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)文章
這篇文章主要介紹了Python的Django框架中從url中捕捉文本的方法,以及URLconf搜索的一些相關(guān)情況,需要的朋友可以參考下2015-07-07
python調(diào)用pytorch的clip模型時(shí)報(bào)錯(cuò)原因及解決方法
這篇文章主要介紹了python調(diào)用pytorch的clip模型時(shí)報(bào)錯(cuò),本文給大家分享問(wèn)題原因及解決方法,本文結(jié)合實(shí)例代碼給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2023-08-08
python 在sql語(yǔ)句中使用%s,%d,%f說(shuō)明
這篇文章主要介紹了python 在sql語(yǔ)句中使用%s,%d,%f說(shuō)明，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-06-06
Python網(wǎng)絡(luò)爬蟲(chóng)信息提取mooc代碼實(shí)例
這篇文章主要介紹了python網(wǎng)絡(luò)爬蟲(chóng)與信息提取mooc,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-03-03
這篇文章主要介紹了Mac下Anaconda的安裝和使用教程，非常不錯(cuò)，具有一定的參考價(jià)借鑒價(jià)值,需要的朋友可以參考下2018-11-11
matplotlib運(yùn)行時(shí)配置(Runtime Configuration，rc)參數(shù)rcParams解析
這篇文章主要介紹了matplotlib運(yùn)行時(shí)配置(Runtime Configuration，rc)參數(shù)rcParams解析，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2021-01-01
這篇文章主要介紹了keras 自定義loss層+接受輸入實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-06-06
這篇文章主要介紹了利用Python實(shí)現(xiàn)批量解壓&壓縮文件夾的示例代碼，文中的實(shí)現(xiàn)步驟講解詳細(xì)，感興趣的小伙伴快跟隨小編一起動(dòng)手試一試2022-04-04
正則表達(dá)式（RE）是一種小型的、高度專(zhuān)業(yè)化的編程語(yǔ)言，它內(nèi)嵌在Python中，并通過(guò)re模塊實(shí)現(xiàn)。下面我們就來(lái)詳細(xì)探討下Python中正則表達(dá)式的使用2016-12-12
大家好，本篇文章主要講的是Python中字典的相關(guān)操作介紹，感興趣的同學(xué)趕快來(lái)看一看吧，對(duì)你有幫助的話記得收藏一下2022-02-02