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标准日志模块logging的使用方法
Nov 01 Python
python中常用的九种预处理方法分享
Sep 11 Python
Django 如何获取前端发送的头文件详解(推荐)
Aug 15 Python
浅谈Matplotlib简介和pyplot的简单使用——文本标注和箭头
Jan 09 Python
matplotlib.pyplot绘图显示控制方法
Jan 15 Python
python实现猜数字游戏
Mar 25 Python
Python使用进程Process模块管理资源
Mar 05 Python
浅析Django 接收所有文件,前端展示文件(包括视频,文件,图片)ajax请求
Mar 09 Python
python脚本监控logstash进程并邮件告警实例
Apr 28 Python
Python实现汇率转换操作
May 03 Python
Pycharm创建文件时自动生成文件头注释(自定义设置作者日期)
Nov 24 Python
python 实现一个简单的线性回归案例
Dec 17 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
PHP配置心得包含MYSQL5乱码解决
2006/11/20 PHP
用PHP函数解决SQL injection
2006/12/09 PHP
PHP strncasecmp字符串比较的小技巧
2011/01/04 PHP
配置Nginx+PHP的正确思路与过程
2016/05/10 PHP
thinkphp中U方法按路由规则生成url的方法
2018/03/12 PHP
php微信开发之关注事件
2018/06/14 PHP
laravel框架中视图的基本使用方法分析
2019/11/23 PHP
js函数中onmousedown和onclick的区别和联系探讨
2013/05/19 Javascript
jQuery实现友好的轮播图片特效
2015/01/12 Javascript
jquery实现网页的页面平滑滚动效果代码
2015/11/02 Javascript
javascript先序遍历DOM树的方法
2016/02/27 Javascript
JavaScript学习笔记整理_简单实现枚举类型,扑克牌应用
2016/09/19 Javascript
原生js封装的一些jquery方法(详解)
2016/09/20 Javascript
Vue数据驱动模拟实现5
2017/01/13 Javascript
JS实现的Unicode编码转换操作示例
2017/04/28 Javascript
详解Angular 4.x NgIf 的用法
2017/05/22 Javascript
深入学习JavaScript 高阶函数
2019/06/11 Javascript
javascript设计模式 ? 命令模式原理与用法实例分析
2020/04/20 Javascript
JavaScript实现拖动对话框效果的实现代码
2020/10/12 Javascript
python数据库操作常用功能使用详解(创建表/插入数据/获取数据)
2013/12/06 Python
Python科学计算环境推荐——Anaconda
2014/06/30 Python
python爬虫入门教程之点点美女图片爬虫代码分享
2014/09/02 Python
Python获取Windows或Linux主机名称通用函数分享
2014/11/22 Python
Python基础篇之初识Python必看攻略
2016/06/23 Python
Python实现二分查找与bisect模块详解
2017/01/13 Python
django 使用 request 获取浏览器发送的参数示例代码
2018/06/11 Python
opencv python Canny边缘提取实现过程解析
2020/02/03 Python
python3:excel操作之读取数据并返回字典 + 写入的案例
2020/09/01 Python
CSS3中Animation动画属性用法详解
2016/07/04 HTML / CSS
使用CSS3制作版头动画效果
2020/12/24 HTML / CSS
运动会四百米广播稿
2014/01/19 职场文书
领导调研接待方案
2014/02/27 职场文书
Python djanjo之csrf防跨站攻击实验过程
2021/05/14 Python
gateway网关接口请求的校验方式
2021/07/15 Java/Android