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 19 Python
深入理解NumPy简明教程---数组1
Dec 17 Python
Pycharm技巧之代码跳转该如何回退
Jul 16 Python
python使用KNN算法手写体识别
Feb 01 Python
基于anaconda下强大的conda命令介绍
Jun 11 Python
Django Admin中增加导出Excel功能过程解析
Sep 04 Python
Python单元测试与测试用例简析
Nov 09 Python
Python高阶函数、常用内置函数用法实例分析
Dec 26 Python
解决python脚本中error: unrecognized arguments: True错误
Apr 20 Python
如何利用Python识别图片中的文字
May 31 Python
Python headers请求头如何实现快速添加
Nov 03 Python
基于Pytorch版yolov5的滑块验证码破解思路详解
Feb 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编程之高级技巧——利用Mysql函数
2006/10/09 PHP
Linux下进行MYSQL编程时插入中文乱码的解决方案
2007/03/15 PHP
遭遇php的in_array低性能问题
2013/09/17 PHP
php实现求相对时间函数
2015/06/15 PHP
深入浅析php中sprintf与printf函数的用法及区别
2016/01/08 PHP
php版微信小店调用api示例代码
2016/11/12 PHP
jQuery.Autocomplete实现自动完成功能(详解)
2010/07/13 Javascript
javascript 获取所有id中包含某关键字的控件的实现代码
2010/11/25 Javascript
jQuery UI的Dialog无法提交问题的解决方法
2011/01/11 Javascript
javascript制作loading动画效果 loading效果
2014/01/14 Javascript
jquery设置按钮停顿3秒不可用
2014/03/07 Javascript
javascript变量声明实例分析
2015/04/25 Javascript
javascript文本模板用法实例
2015/07/31 Javascript
谈谈我对JavaScript原型和闭包系列理解(随手笔记8)
2015/12/24 Javascript
微信小程序 获取设备信息 API实例详解
2016/10/02 Javascript
canvas 弹幕效果(实例分享)
2017/01/11 Javascript
Vue.2.0.5过渡效果使用技巧
2017/03/16 Javascript
Vue中计算属性computed的示例解读
2017/07/26 Javascript
js实现滑动滑块验证登录
2020/07/24 Javascript
JQuery基于FormData异步提交数据文件
2020/09/01 jQuery
[50:20]DOTA2上海特级锦标赛主赛事日 - 5 总决赛Liquid VS Secret第四局
2016/03/06 DOTA
[03:19]2016国际邀请赛中国区预选赛第四日TOP10镜头集锦
2016/07/01 DOTA
Win10下Python环境搭建与配置教程
2016/11/18 Python
Python登录并获取CSDN博客所有文章列表代码实例
2017/12/28 Python
Python中GIL的使用详解
2018/10/03 Python
python实现飞机大战游戏(pygame版)
2020/10/26 Python
python实现贪吃蛇游戏源码
2020/03/21 Python
Python tkinter实现简单加法计算器代码实例
2020/05/13 Python
详解Django关于StreamingHttpResponse与FileResponse文件下载的最优方法
2021/01/07 Python
分享CSS3制作卡片式图片的方法
2016/07/08 HTML / CSS
浅谈CSS3特性查询(Feature Query: @supports)功能简介
2017/07/31 HTML / CSS
linux面试题参考答案(11)
2012/05/01 面试题
国贸专业自荐信范文
2014/03/02 职场文书
导游词之上饶龟峰
2019/10/25 职场文书
修改并编译golang源码的操作步骤
2021/07/25 Golang
Mysql将字符串按照指定字符分割的正确方法
2022/05/30 MySQL