快捷導(dǎo)航

Python創(chuàng)建對稱矩陣的方法示例【基于numpy模塊】

更新時間：2017年10月12日 10:30:51 作者：Inside_Zhang

這篇文章主要介紹了Python創(chuàng)建對稱矩陣的方法,結(jié)合實例形式分析了Python基于numpy模塊實現(xiàn)矩陣運算的相關(guān)操作技巧,需要的朋友可以參考下

本文實例講述了Python創(chuàng)建對稱矩陣的方法。分享給大家供大家參考，具體如下：

對稱（實對稱）矩陣也即：

step 1：創(chuàng)建一個方陣

>>> import numpy as np
>>> X = np.random.rand(5**2).reshape(5, 5)
>>> X
array([[ 0.26984148, 0.25408384, 0.12428487, 0.0194565 , 0.91287708],
  [ 0.31837673, 0.35493156, 0.74336268, 0.31810561, 0.04409245],
  [ 0.06644445, 0.8967897 , 0.10990936, 0.05036292, 0.72581982],
  [ 0.94758512, 0.21375975, 0.36781736, 0.1633904 , 0.36070709],
  [ 0.53263787, 0.18380491, 0.0225521 , 0.91239367, 0.75521585]])

step 2：保留其上三角部分

>>> X = np.triu(X)
# 保留其上三角部分
>>> X
array([[ 0.26984148, 0.25408384, 0.12428487, 0.0194565 , 0.91287708],
  [ 0.  , 0.35493156, 0.74336268, 0.31810561, 0.04409245],
  [ 0.  , 0.  , 0.10990936, 0.05036292, 0.72581982],
  [ 0.  , 0.  , 0.  , 0.1633904 , 0.36070709],
  [ 0.  , 0.  , 0.  , 0.  , 0.75521585]])

step 3：將上三角”拷貝”到下三角部分

>>> X += X.T - np.diag(X.diagonal())
>>> X
array([[ 0.26984148, 0.25408384, 0.12428487, 0.0194565 , 0.91287708],
  [ 0.25408384, 0.35493156, 0.74336268, 0.31810561, 0.04409245],
  [ 0.12428487, 0.74336268, 0.10990936, 0.05036292, 0.72581982],
  [ 0.0194565 , 0.31810561, 0.05036292, 0.1633904 , 0.36070709],
  [ 0.91287708, 0.04409245, 0.72581982, 0.36070709, 0.75521585]])

注意，要減去一次對角線上的元素。因為上三角cov，和下三角cov.T在進行相加時會把主對角線上的元素相加兩次。

step 4：測試

>>> X.T == X
array([[ True, True, True, True, True],
  [ True, True, True, True, True],
  [ True, True, True, True, True],
  [ True, True, True, True, True],
  [ True, True, True, True, True]], dtype=bool)

更多關(guān)于Python相關(guān)內(nèi)容可查看本站專題：《Python數(shù)學(xué)運算技巧總結(jié)》、《Python正則表達(dá)式用法總結(jié)》、《Python數(shù)據(jù)結(jié)構(gòu)與算法教程》、《Python函數(shù)使用技巧總結(jié)》、《Python字符串操作技巧匯總》及《Python入門與進階經(jīng)典教程》

希望本文所述對大家Python程序設(shè)計有所幫助。

您可能感興趣的文章:

相關(guān)文章

Python運行異常管理解決方案
這篇文章主要介紹了Python運行異常管理解決方案,文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友可以參考下
2020-03-03
Python中進度條tqdm包使用方法及特性介紹
這篇文章主要介紹了Python中進度條tqdm包使用方法及特性,具有很好的參考價值,希望對大家有所幫助,如有錯誤或未考慮完全的地方,望不吝賜教
2023-08-08
python對DICOM圖像的讀取方法詳解
DICOM（Digital Imaging and Communications in Medicine）即醫(yī)學(xué)數(shù)字成像和通信，是醫(yī)學(xué)圖像和相關(guān)信息的國際標(biāo)準(zhǔn)（ISO 12052）。下面這篇文章主要給大家介紹了關(guān)于python對DICOM圖像讀取的相關(guān)資料，需要的朋友可以參考借鑒，下面來一起看看吧。
2017-07-07
python采用getopt解析命令行輸入?yún)?shù)實例
這篇文章主要介紹了python采用getopt解析命令行輸入?yún)?shù)實例,非常具有實用價值,需要的朋友可以參考下
2014-09-09
Python循環(huán)語句之break與continue的用法
這篇文章主要介紹了Python循環(huán)語句之break與continue的用法,是Python入門學(xué)習(xí)中的基礎(chǔ)知識,需要的朋友可以參考下
2015-10-10
LRUCache的實現(xiàn)原理及利用python實現(xiàn)的方法
LruCache 是 Android 的一個內(nèi)部類，提供了基于內(nèi)存實現(xiàn)的緩存，而下面這篇文章主要給大家介紹了關(guān)于LRUCache的實現(xiàn)原理以及利用python實現(xiàn)的方法，文中通過示例代碼介紹的非常詳細(xì)，需要的朋友可以參考借鑒，下面來一起看看吧。
2017-11-11
json 轉(zhuǎn) mot17數(shù)據(jù)格式的實現(xiàn)代碼 (親測有效)
這篇文章主要介紹了json 轉(zhuǎn) mot17數(shù)據(jù)格式的實現(xiàn)代碼 (親測有效),本文給大家介紹的非常詳細(xì)，對大家的學(xué)習(xí)或工作具有一定的參考借鑒價值，需要的朋友可以參考下
2021-03-03
python實現(xiàn)Floyd算法
這篇文章主要為大家詳細(xì)介紹了python實現(xiàn)Floyd算法，具有一定的參考價值，感興趣的小伙伴們可以參考一下
2018-01-01
numpy如何按條件給元素賦值np.where、np.clip
這篇文章主要介紹了numpy如何按條件給元素賦值np.where、np.clip問題，具有很好的參考價值，希望對大家有所幫助。如有錯誤或未考慮完全的地方，望不吝賜教
2023-06-06
Python輸出列表(List)不帶中括號和引號的問題及解決方法
這篇文章主要介紹了Python輸出列表(List)不帶中括號和引號的問題及解決方法，本文給大家介紹的非常詳細(xì)，對大家的學(xué)習(xí)或工作具有一定的參考借鑒價值,需要的朋友可以參考下
2023-02-02