python/sympy求解矩陣方程的方法
sympy版本:1.2
假設(shè)求解矩陣方程
AX=A+2X
其中
求解之前對矩陣方程化簡為
(A−2E)X=A
令
B=(A−2E)
使用qtconsole輸入下面程序進行求解
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é)果驗證一下:
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]])
求解矩陣方程過程中注意的問題是左乘還是右乘問題,在此例中是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é)果驗證一下:
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求解矩陣方程的方法就是小編分享給大家的全部內(nèi)容了,希望能給大家一個參考,也希望大家多多支持腳本之家。
相關(guān)文章
python調(diào)用pytorch的clip模型時報錯原因及解決方法
這篇文章主要介紹了python調(diào)用pytorch的clip模型時報錯,本文給大家分享問題原因及解決方法,本文結(jié)合實例代碼給大家介紹的非常詳細(xì),對大家的學(xué)習(xí)或工作具有一定的參考借鑒價值,需要的朋友可以參考下2023-08-08Python網(wǎng)絡(luò)爬蟲信息提取mooc代碼實例
這篇文章主要介紹了python網(wǎng)絡(luò)爬蟲與信息提取mooc,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友可以參考下2020-03-03matplotlib運行時配置(Runtime Configuration,rc)參數(shù)rcParams解析
這篇文章主要介紹了matplotlib運行時配置(Runtime Configuration,rc)參數(shù)rcParams解析,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2021-01-01