您好,登錄后才能下訂單哦!
這篇文章主要講解了基于python檢查矩陣計(jì)算結(jié)果的方法,內(nèi)容清晰明了,對(duì)此有興趣的小伙伴可以學(xué)習(xí)一下,相信大家閱讀完之后會(huì)有幫助。
鑒于最近復(fù)習(xí)線性代數(shù)計(jì)算量較大,且1800答案常常忽略一些逆陣、行列式的計(jì)算答案,故用Python寫出矩陣的簡單計(jì)算程序,便于檢查出錯(cuò)的步驟。
1、行列式
可自行更改階數(shù)
from numpy import * # 求行列式 ,建議:取小數(shù)點(diǎn)前整數(shù) A = array([[3, 1, 1, 1], [1, 3, 1, 1], [1, 1, 3, 1], [1, 1, 1, 3]]) B = linalg.det(A) print(B) # 48.000000000000014 正確答案:48
2、矩陣相乘
注意要內(nèi)標(biāo)相同
from numpy import * # 求矩陣相乘 A = array([[1, -1, 1], [1, 1, 0], [-1, 0, 1]]) B = array([[3, 0, 0], [0, 0, 0], [0, 0, 0]]) # N=AB N = dot(A, B) # N=BA,則 N = dot(B, A) print(N) # 正確答案: # [ 3 0 0] # [ 3 0 0] # [-3 0 0]
3、逆矩陣
自行判斷|A|≠0,這里 A∗ = A−1 · |A|
from numpy import * # 求逆矩陣 ,建議:取小數(shù)點(diǎn)后一位化為分?jǐn)?shù) A = mat([[1, -1, 1], [1, 1, 0], [-1, 0, 1]]) B = A.I print(B) # [ 0.33333333 0.33333333 -0.33333333] # [-0.33333333 0.66666667 0.33333333] # [ 0.33333333 0.33333333 0.66666667] # 0.333≈ 1/3 ,0.667≈ 2/3
看完上述內(nèi)容,是不是對(duì)基于python檢查矩陣計(jì)算結(jié)果的方法有進(jìn)一步的了解,如果還想學(xué)習(xí)更多內(nèi)容,歡迎關(guān)注億速云行業(yè)資訊頻道。
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。