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

yipeiwu_com5年前Python基础

sympy版本:1.2

假设求解矩阵方程

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]])

将结果验证一下:

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 ,如果矩阵方程变为

XA=A+2X

那么求解结果为:

In [35]: X=A*B.inv()

In [36]: X
Out[36]: 
Matrix([
[ 3, -8, -6],
[ 2, -9, -6],
[-2, 12, 9]])

将结果验证一下:

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求解矩阵方程的方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持【听图阁-专注于Python设计】。

相关文章

python 执行终端/控制台命令的例子

如下所示: import os os.system() os.popen().read().strip() #上面2种方法 是python 执行终端/控制台 命令的常见方法 #os...

Python3标准库总结

Python3标准库 操作系统接口 os模块提供了不少与操作系统相关联的函数。 >>> import os >>> os.getcwd()...

浅谈Python3识别判断图片主要颜色并和颜色库进行对比的方法

浅谈Python3识别判断图片主要颜色并和颜色库进行对比的方法

【更新】主要提供两种方案: 方案一:(参考网上代码,感觉实用性不是很强)使用PIL截取图像,然后将RGB转为HSV进行判断,统计判断颜色,最后输出RGB值 方案二:使用opencv库函数...

python操作xlsx文件的包openpyxl实例

Python扩展库openpyxl,可以操作07版以上的xlsx文件。可以创建工作簿、选择活动工作表、写入单元格数据,设置单元格字体颜色、边框样式,合并单元格,设置单元格背景等等。 需要...

用Python配平化学方程式的方法

用Python配平化学方程式的方法

不知不觉已经毕业多年了,不知道大家是否还记得怎么配平化学方程式呢?反正小编我是已经记不太清了,所以今天的文章除了分享如何用python配平化学方程式,顺带着还会复习 一些化学方程式的知识...