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字符串加密解密的三种方法分享(base64 win32com)
Jan 19 Python
Python使用matplotlib实现在坐标系中画一个矩形的方法
May 20 Python
Python中函数参数设置及使用的学习笔记
May 03 Python
python安装oracle扩展及数据库连接方法
Feb 21 Python
python实现决策树、随机森林的简单原理
Mar 26 Python
使用anaconda的pip安装第三方python包的操作步骤
Jun 11 Python
Python Pandas 获取列匹配特定值的行的索引问题
Jul 01 Python
详解python调用cmd命令三种方法
Jul 08 Python
Python 格式化打印json数据方法(展开状态)
Feb 27 Python
python闭包、深浅拷贝、垃圾回收、with语句知识点汇总
Mar 11 Python
使用keras框架cnn+ctc_loss识别不定长字符图片操作
Jun 29 Python
2021年值得向Python开发者推荐的VS Code扩展插件
Jan 25 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 函数学习简单小结
2010/07/08 PHP
安装ImageMagick出现error while loading shared libraries的解决方法
2014/09/23 PHP
Yii入门教程之目录结构、入口文件及路由设置
2014/11/25 PHP
PHP生成图片验证码功能示例
2017/01/12 PHP
windows下的WAMP环境搭建图文教程(推荐)
2017/07/27 PHP
nodejs初步体验篇
2015/11/23 NodeJs
图解Sublime Text3使用技巧
2015/12/21 Javascript
JavaScript中各种引用类型的常用操作方法小结
2016/05/05 Javascript
使用requirejs模块化开发多页面一个入口js的使用方式
2017/06/14 Javascript
基于javascript 显式转换与隐式转换(详解)
2017/12/15 Javascript
JavaScript命名空间模式实例详解
2019/06/20 Javascript
jquery实现两个div中的元素相互拖动的方法分析
2020/04/05 jQuery
解决vue+elementui项目打包后样式变化问题
2020/08/03 Javascript
jQuery实现动态操作table行
2020/11/23 jQuery
wxPython 入门教程
2008/10/07 Python
Django1.9 加载通过ImageField上传的图片方法
2018/05/25 Python
django连接mysql配置方法总结(推荐)
2018/08/18 Python
Django之Mode的外键自关联和引用未定义的Model方法
2018/12/15 Python
Python实现批量执行同目录下的py文件方法
2019/01/11 Python
Python使用numpy模块实现矩阵和列表的连接操作方法
2019/06/26 Python
解决python中导入win32com.client出错的问题
2019/07/26 Python
Python time库基本使用方法分析
2019/12/13 Python
python 实现Flask中返回图片流给前端展示
2020/01/09 Python
浅谈tensorflow中Dataset图片的批量读取及维度的操作详解
2020/01/20 Python
Python列表操作方法详解
2020/02/09 Python
如何使用python切换hosts文件
2020/04/29 Python
Python利用matplotlib绘制散点图的新手教程
2020/11/05 Python
HTML5之消息通知的使用(Web Notification)
2018/10/30 HTML / CSS
一套C++笔试题面试题
2012/06/06 面试题
员工团队活动方案
2014/08/28 职场文书
防灾减灾标语
2014/10/07 职场文书
新郎婚礼答谢词
2015/01/04 职场文书
婚宴父母致辞
2015/07/27 职场文书
2016年党员公开承诺书范文
2016/03/24 职场文书
【超详细】八大排序算法的各项比较以及各自特点
2021/03/31 Python
使用vuex-persistedstate本地存储vuex
2022/04/29 Vue.js