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


Posted in Python onNovember 08, 2018

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制作一个桌面便签软件
Aug 09 Python
Python3.6正式版新特性预览
Dec 15 Python
Python实现通过文件路径获取文件hash值的方法
Apr 29 Python
python+pandas生成指定日期和重采样的方法
Apr 11 Python
python实现下载pop3邮件保存到本地
Jun 19 Python
Python时间序列处理之ARIMA模型的使用讲解
Apr 02 Python
使用Python打造一款间谍程序的流程分析
Feb 21 Python
解决pycharm中opencv-python导入cv2后无法自动补全的问题(不用作任何文件上的修改)
Mar 05 Python
PHP基于phpqrcode类库生成二维码过程解析
May 28 Python
DRF框架API版本管理实现方法解析
Aug 21 Python
pycharm配置QtDesigner的超详细方法
Jan 25 Python
Django如何创作一个简单的最小程序
May 12 Python
python生成lmdb格式的文件实例
Nov 08 #Python
python实现嵌套列表平铺的两种方法
Nov 08 #Python
python用列表生成式写嵌套循环的方法
Nov 08 #Python
在Python中实现shuffle给列表洗牌
Nov 08 #Python
python实现RabbitMQ的消息队列的示例代码
Nov 08 #Python
对Python 3.5拼接列表的新语法详解
Nov 08 #Python
Python使用random.shuffle()打乱列表顺序的方法
Nov 08 #Python
You might like
总集篇&特番节目先行播出!《SAO Alicization War of Underworld》第2季度TV动画4月25日放送!
2020/03/06 日漫
解决phpmyadmin中缺少mysqli扩展问题的方法
2013/05/06 PHP
下载站控制介绍字数显示的脚本 显示全部 隐藏介绍等功能
2009/09/19 Javascript
通过JS 获取Mouse Position(鼠标坐标)的代码
2009/09/21 Javascript
jquery1.4 教程二 ajax方法的改进
2010/02/25 Javascript
19个很有用的 JavaScript库推荐
2011/06/27 Javascript
Js数组的操作push,pop,shift,unshift等方法详细介绍
2012/12/28 Javascript
JS实现多物体缓冲运动实例代码
2013/11/29 Javascript
jQuery实现的在线答题功能
2015/04/12 Javascript
javascript转换静态图片,增加粒子动画效果
2015/05/28 Javascript
js实现三张图(文)片一起切换的banner焦点图
2015/08/25 Javascript
Augularjs-起步详解
2016/07/08 Javascript
Vue.js 表单校验插件
2016/08/14 Javascript
浅谈angular4 ng-content 中隐藏的内容
2017/08/18 Javascript
es6系列教程_ Map详解以及常用api介绍
2017/09/25 Javascript
AngularJS实现注册表单验证功能
2017/10/16 Javascript
node.js+express+mySQL+ejs+bootstrop实现网站登录注册功能
2018/01/12 Javascript
解决js ajax同步请求造成浏览器假死的问题
2018/01/18 Javascript
详解vue-cli 本地开发mock数据使用方法
2018/05/29 Javascript
layer弹出的iframe层在执行完毕后关闭当前弹出层的方法
2018/08/17 Javascript
elementUI 设置input的只读或禁用的方法
2018/10/30 Javascript
jQuery实现小火箭返回顶部特效
2020/02/03 jQuery
[01:10:24]DOTA2-DPC中国联赛 正赛 VG vs Aster BO3 第一场 2月28日
2021/03/11 DOTA
django通过ajax发起请求返回JSON格式数据的方法
2015/06/04 Python
Python实现批量下载图片的方法
2015/07/08 Python
详解Python3中的Sequence type的使用
2015/08/01 Python
python添加模块搜索路径方法
2017/09/11 Python
Flask框架信号用法实例分析
2018/07/24 Python
django解决跨域请求的问题
2018/11/11 Python
对DataFrame数据中的重复行,利用groupby累加合并的方法详解
2019/01/30 Python
python字典key不能是可以是啥类型
2020/08/04 Python
Surfdome西班牙:世界上最受欢迎的生活方式品牌
2019/02/13 全球购物
yy结婚证婚词
2014/01/10 职场文书
教师自我剖析材料范文
2014/09/30 职场文书
给老婆的保证书怎么写
2015/05/08 职场文书
SQL SERVER中常用日期函数的具体使用
2021/04/08 SQL Server