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类参数self使用示例
Feb 17 Python
python实现提取百度搜索结果的方法
May 19 Python
Python中pip安装非PyPI官网第三方库的方法
Jun 02 Python
scrapy爬虫实例分享
Dec 28 Python
mac系统安装Python3初体验
Jan 02 Python
Python 12306抢火车票脚本 Python京东抢手机脚本
Feb 06 Python
Django读取Mysql数据并显示在前端的实例
May 27 Python
numpy.linspace 生成等差数组的方法
Jul 02 Python
使用coverage统计python web项目代码覆盖率的方法详解
Aug 05 Python
python构造函数init实例方法解析
Jan 19 Python
python图形开发GUI库pyqt5的详细使用方法及各控件的属性与方法
Feb 14 Python
django 解决model中类写不到数据库中,数据库无此字段的问题
May 20 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
用PHP调用数据库的存贮过程!
2006/10/09 PHP
phpadmin如何导入导出大数据文件及php.ini参数修改
2013/02/18 PHP
基于Web标准的UI组件 — 树状菜单(2)
2006/09/18 Javascript
javascript编程起步(第七课)
2007/02/27 Javascript
很全的显示阴历(农历)日期的js代码
2009/01/01 Javascript
js运动框架_包括图片的淡入淡出效果
2013/05/11 Javascript
jquery判断密码强度的验证代码
2020/04/22 Javascript
JavaScript实现瀑布流以及加载效果
2017/02/11 Javascript
Node.js 8 中的重要新特性
2017/06/28 Javascript
React Native中Navigator的使用方法示例
2017/10/13 Javascript
浅谈node模块与npm包管理工具
2018/01/03 Javascript
解决angularjs WdatePicker ng-model的问题
2018/09/13 Javascript
Vue.js实现可排序的表格组件功能示例
2019/02/19 Javascript
如何实现iframe父子传参通信
2020/02/05 Javascript
[01:14]英雄,所敬略同——2018完美盛典宣传视频4K
2018/12/05 DOTA
Python多线程和队列操作实例
2015/06/21 Python
教大家使用Python SqlAlchemy
2016/02/12 Python
Python实现调度算法代码详解
2017/12/01 Python
Python之reload流程实例代码解析
2018/01/29 Python
Python3.5.3下配置opencv3.2.0的操作方法
2018/04/02 Python
PyQt5实现无边框窗口的标题拖动和窗口缩放
2018/04/19 Python
python 怎样将dataframe中的字符串日期转化为日期的方法
2019/09/26 Python
基于python3的socket聊天编程
2020/02/17 Python
Python requests上传文件实现步骤
2020/09/15 Python
微信html5页面调用第三方位置导航的示例
2018/03/14 HTML / CSS
世界第一曲奇连锁店:Mrs. Fields Cookies
2017/02/04 全球购物
英国Radley包德国官网:Radley London德国
2019/11/18 全球购物
MATCHESFASHION澳大利亚/亚太地区:英国时尚奢侈品电商
2020/01/14 全球购物
毕业生个人求职自荐信
2014/02/26 职场文书
应届毕业生自荐信
2014/05/28 职场文书
员工趣味活动方案
2014/08/27 职场文书
个人总结怎么写
2015/02/26 职场文书
推销搭讪开场白
2015/05/28 职场文书
开票证明
2015/06/23 职场文书
企业安全隐患排查治理制度
2015/08/05 职场文书
Win11电脑显示本地时间与服务器时间不一致怎么解决?
2022/04/05 数码科技