python/sympy求解矩陣方程的方法

發(fā)布時間：2020-09-30 01:57:40 來源：腳本之家閱讀：280 作者：落葉_小唱欄目：開發(fā)技術(shù)

sympy版本:1.2

假設(shè)求解矩陣方程

AX=A+2X

其中

python/sympy求解矩陣方程的方法

求解之前對矩陣方程化簡為

(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)容了，希望能給大家一個參考，也希望大家多多支持億速云。

向AI問一下細節(jié)

python/sympy求解矩陣方程的方法

猜你喜歡

最新資訊

相關(guān)推薦

相關(guān)標簽