python/sympy求解矩陣方程的方法
更新時(shí)間：2018年11月08日 11:36:06 作者：落葉_小唱
今天小編就為大家分享一篇python/sympy求解矩陣方程的方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧
sympy版本:1.2
假設(shè)求解矩陣方程
AX=A+2X
其中
求解之前對(duì)矩陣方程化簡(jiǎn)為
(A−2E)X=A
令
B=(A−2E)
使用qtconsole輸入下面程序進(jìn)行求解
In [26]: from sympy import * In [27]: from sympy.abc import * In [28]: A=Matrix([[4,2,3],[1,1,0],[-1,2,3]]) In [29]: A Out[29]: Matrix([ [ 4, 2, 3], [ 1, 1, 0], [-1, 2, 3]]) In [30]: B=A-2*diag(1,1,1) In [31]: B Out[31]: Matrix([ [ 2, 2, 3], [ 1, -1, 0], [-1, 2, 1]]) In [32]: B.inv()*A Out[32]: Matrix([ [ 3, -8, -6], [ 2, -9, -6], [-2, 12, 9]])
將結(jié)果驗(yàn)證一下：
In [38]: X=B.inv()*A In [39]: X Out[39]: Matrix([ [ 3, -8, -6], [ 2, -9, -6], [-2, 12, 9]]) In [40]: A*X-A-2*X Out[40]: Matrix([ [0, 0, 0], [0, 0, 0], [0, 0, 0]])
求解矩陣方程過(guò)程中注意的問(wèn)題是左乘還是右乘問(wèn)題，在此例中是B.inv()*A ，如果矩陣方程變?yōu)?/p>
XA=A+2X
那么求解結(jié)果為：
In [35]: X=A*B.inv() In [36]: X Out[36]: Matrix([ [ 3, -8, -6], [ 2, -9, -6], [-2, 12, 9]])
將結(jié)果驗(yàn)證一下：
X=A*B.inv() X Out[36]: Matrix([ [ 3, -8, -6], [ 2, -9, -6], [-2, 12, 9]]) X*A-A-2*X Out[37]: Matrix([ [0, 0, 0], [0, 0, 0], [0, 0, 0]])
以上這篇python/sympy求解矩陣方程的方法就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。
相關(guān)文章
這篇文章主要為大家詳細(xì)介紹了python定向爬取淘寶商品價(jià)格，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2018-02-02
使用pymysql查詢(xún)數(shù)據(jù)庫(kù),把結(jié)果保存為列表并獲取指定元素下標(biāo)實(shí)例
這篇文章主要介紹了使用pymysql查詢(xún)數(shù)據(jù)庫(kù),把結(jié)果保存為列表并獲取指定元素下標(biāo)實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-05-05
在學(xué)習(xí)Python爬蟲(chóng)部分，需要已經(jīng)學(xué)過(guò)Python基礎(chǔ)和前端的相關(guān)知識(shí)，本文對(duì)python爬蟲(chóng)概念及原理給大家詳細(xì)介紹，需要的朋友跟隨小編一起看看吧2021-04-04
今天小編就為大家分享一篇python讀取文件名并改名字的實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-01-01
Python無(wú)法用requests獲取網(wǎng)頁(yè)源碼的解決方法
爬蟲(chóng)獲取信息,很多時(shí)候是需要從網(wǎng)頁(yè)源碼中獲取鏈接信息的,下面這篇文章主要給大家介紹了關(guān)于Python無(wú)法用requests獲取網(wǎng)頁(yè)源碼的解決方法,文中通過(guò)示例代碼介紹的非常詳細(xì),需要的朋友可以參考下2022-07-07
這篇文章主要介紹了Python中的zip函數(shù)使用示例,Python中的zip函數(shù)不是用來(lái)壓縮文件的,而是用來(lái)傳入?yún)?shù)返回一個(gè)tuple列表的,需要的朋友可以參考下2015-01-01
這篇文章主要介紹了python 發(fā)送get請(qǐng)求接口詳解,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2020-11-11
今天小編就為大家分享一篇PyCharm配置mongo插件的方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2018-11-11
聊聊boost?python3依賴(lài)安裝問(wèn)題
這篇文章主要介紹了boost?python3依賴(lài)安裝,本文給大家介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值，需要的朋友可以參考下2021-12-12
python利用google翻譯方法實(shí)例（翻譯字幕文件）
這篇文章主要給大家介紹了關(guān)于python利用google翻譯（翻譯字幕文件）的相關(guān)資料，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2020-09-09