使用Python怎么求逆矩陣?很多新手對(duì)此不是很清楚,為了幫助大家解決這個(gè)難題,下面小編將為大家詳細(xì)講解,有這方面需求的人可以來(lái)學(xué)習(xí)下,希望你能有所收獲。
站在用戶(hù)的角度思考問(wèn)題,與客戶(hù)深入溝通,找到江陵網(wǎng)站設(shè)計(jì)與江陵網(wǎng)站推廣的解決方案,憑借多年的經(jīng)驗(yàn),讓設(shè)計(jì)與互聯(lián)網(wǎng)技術(shù)結(jié)合,創(chuàng)造個(gè)性化、用戶(hù)體驗(yàn)好的作品,建站類(lèi)型包括:成都網(wǎng)站建設(shè)、網(wǎng)站建設(shè)、企業(yè)官網(wǎng)、英文網(wǎng)站、手機(jī)端網(wǎng)站、網(wǎng)站推廣、國(guó)際域名空間、網(wǎng)頁(yè)空間、企業(yè)郵箱。業(yè)務(wù)覆蓋江陵地區(qū)。import numpy as np kernel = np.array([1, 1, 1, 2]).reshape((2, 2)) print(kernel) print(np.linalg.inv(kernel))
注意,Singular matrix奇異矩陣不可求逆
補(bǔ)充:python+numpy中矩陣的逆和偽逆的區(qū)別
對(duì)于矩陣A,如果存在一個(gè)矩陣B,使得AB=BA=E,其中E為與A,B同維數(shù)的單位陣,就稱(chēng)A為可逆矩陣(或者稱(chēng)A可逆),并稱(chēng)B是A的逆矩陣,簡(jiǎn)稱(chēng)逆陣。(此時(shí)的逆稱(chēng)為凱利逆)
矩陣A可逆的充分必要條件是|A|≠0。
偽逆矩陣是逆矩陣的廣義形式。由于奇異矩陣或非方陣的矩陣不存在逆矩陣,但可以用函數(shù)pinv(A)求其偽逆矩陣。
基本語(yǔ)法為X=pinv(A),X=pinv(A,tol),其中tol為誤差,pinv為pseudo-inverse的縮寫(xiě):max(size(A))*norm(A)*eps。
函數(shù)返回一個(gè)與A的轉(zhuǎn)置矩陣A' 同型的矩陣X,并且滿(mǎn)足:AXA=A,XAX=X.此時(shí),稱(chēng)矩陣X為矩陣A的偽逆,也稱(chēng)為廣義逆矩陣。
pinv(A)具有inv(A)的部分特性,但不與inv(A)完全等同。
如果A為非奇異方陣,pinv(A)=inv(A),但卻會(huì)耗費(fèi)大量的計(jì)算時(shí)間,相比較而言,inv(A)花費(fèi)更少的時(shí)間。
代碼如下:
import numpy as np a = np.array([[1, 2], [3, 4]]) # 初始化一個(gè)非奇異矩陣(數(shù)組) print(np.linalg.inv(a)) # 對(duì)應(yīng)于MATLAB中 inv() 函數(shù) # 矩陣對(duì)象可以通過(guò) .I 求逆,但必須先使用matirx轉(zhuǎn)化 A = np.matrix(a) print(A.I)
import numpy as np # 定義一個(gè)奇異陣 A A = np.zeros((4, 4)) A[0, -1] = 1 A[-1, 0] = -1 A = np.matrix(A) print(A) # print(A.I) 將報(bào)錯(cuò),矩陣 A 為奇異矩陣,不可逆 print(np.linalg.pinv(A)) # 求矩陣 A 的偽逆(廣義逆矩陣),對(duì)應(yīng)于MATLAB中 pinv() 函數(shù)
這就是矩陣的逆和偽逆的區(qū)別
截至2020/10/4,matrix函數(shù)還可以使用,但已經(jīng)過(guò)時(shí),應(yīng)該是mat函數(shù)這種。
看完上述內(nèi)容是否對(duì)您有幫助呢?如果還想對(duì)相關(guān)知識(shí)有進(jìn)一步的了解或閱讀更多相關(guān)文章,請(qǐng)關(guān)注創(chuàng)新互聯(lián)行業(yè)資訊頻道,感謝您對(duì)創(chuàng)新互聯(lián)網(wǎng)站建設(shè)公司,的支持。
網(wǎng)頁(yè)題目:使用Python怎么求逆矩陣-創(chuàng)新互聯(lián)
鏈接URL:http://aaarwkj.com/article0/dpjjio.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供網(wǎng)站設(shè)計(jì)公司、軟件開(kāi)發(fā)、營(yíng)銷(xiāo)型網(wǎng)站建設(shè)、網(wǎng)站內(nèi)鏈、面包屑導(dǎo)航、網(wǎng)站導(dǎo)航
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶(hù)投稿、用戶(hù)轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如需處理請(qǐng)聯(lián)系客服。電話(huà):028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來(lái)源: 創(chuàng)新互聯(lián)
猜你還喜歡下面的內(nèi)容