基于python檢查矩陣計(jì)算結(jié)果
鑒于最近復(fù)習(xí)線性代數(shù)計(jì)算量較大,且1800答案常常忽略一些逆陣、行列式的計(jì)算答案,故用Python寫出矩陣的簡(jiǎn)單計(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=ABN = 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.Iprint(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
以上就是本文的全部?jī)?nèi)容,希望對(duì)大家的學(xué)習(xí)有所幫助,也希望大家多多支持好吧啦網(wǎng)。
相關(guān)文章:
1. 以PHP代碼為實(shí)例詳解RabbitMQ消息隊(duì)列中間件的6種模式2. html小技巧之td,div標(biāo)簽里內(nèi)容不換行3. PHP字符串前后字符或空格刪除方法介紹4. 將properties文件的配置設(shè)置為整個(gè)Web應(yīng)用的全局變量實(shí)現(xiàn)方法5. nestjs實(shí)現(xiàn)圖形校驗(yàn)和單點(diǎn)登錄的示例代碼6. AspNetCore&MassTransit Courier實(shí)現(xiàn)分布式事務(wù)的詳細(xì)過程7. XML入門的常見問題(一)8. jsp cookie+session實(shí)現(xiàn)簡(jiǎn)易自動(dòng)登錄9. css進(jìn)階學(xué)習(xí) 選擇符10. Echarts通過dataset數(shù)據(jù)集實(shí)現(xiàn)創(chuàng)建單軸散點(diǎn)圖
