python/sympy求解矩阵方程的方法
作者：落叶_小唱 发布时间：2023-07-04 07:06:34
标签：python,sympy,矩阵,方程
sympy版本:1.2
假设求解矩阵方程
AX=A+2X
其中
求解之前对矩阵方程化简为
(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]])
来源：https://blog.csdn.net/ouening/article/details/82711761
0
投稿猜你喜欢
- 代码如下:function checkip(checkstring)'用正则判断IP是否合法 dim re1 set re1=new
0. 我们如何通过邮件系统完成远程控制电脑（关机、重启等）？实现思路：需要有两个邮箱：接收指令邮箱（A）发送指令邮箱（B）被控制的电脑（查看- 一.思路我们通过网页版的微信公众平台的图文消息中的超链接获取到我们需要的接口从接口中我们可以得到对应的微信公众号和对应的所有微信公众号文章。
- 上期回顾:亚马逊购物用户体验分析 (二)方便的导航元素任何网上商店的成功，至关重要的一点，就是用户可以简单轻松地使用导航条。基本店铺分类，用
- 但是如果是让你接手一个二等残废的网站，并让你在上面改版，而且不能推翻式改版，只能逐步替换旧的程序，那么你会非常痛苦，例如我遇到的问题： 问题
- 这篇文章主要介绍了python框架django项目部署相关知识详解,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价
- 本文实例讲述了python处理csv数据的方法。分享给大家供大家参考。具体如下：Python代码：#coding=utf-8__author
- 一、说明：urllib库是python内置的一个http请求库，requests库就是基于该库开发出来的，虽然requests库使用更方便，
- 导语：使用 python-plotly 模块来进行压测数据的绘制，并且生成静态 html 页面结果展示。不少小伙伴在开发过程中都有对模块进行
- 废话不多说了，关键代码如下所示：<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 T
- 学习网络爬虫难免遇到使用代理的情况，下面介绍一下如何使用requests设置代理：如果需要使用代理，你可以通过为任意请求方法提供 proxi
- javascript可以根据输入值自动搜索显示相关的select列表，对于列表很长时可以很方便的查找到要的值。js代码：<script
- 关于数据库的逻辑设计，是一个很广泛的问题。本文主要针对开发应用中遇到在MS SQL Server上进行表设计时，对表的主键设计应注意的问题以
- 图片提取为了方便技术展示，我们选取素材为演员杨紫的一段演讲视频，用例仅为技术交流演示使用，不针对任何指定人。为达到我们AI换脸的目的，我们首
- 1. 预估执行计划 - Explain PlanExplain plan以SQL语句作为输入，得到这条SQL语句的执行计划，并将执行计划输出
- ctrl + r => 输入drivers回车 => etc/hosts , 用记事本打开它,在 127.0.0.1 local
- Jaslabs的Justin Silverton列出了十条有关优化MySQL查询的语句，我不得不对此发表言论，因为这个清单非常非常糟糕。另外
- 新一代W3C，xhtml代码规范，大家在设计网站的时候务必遵循这一规范 ，这将对网站的优化，网站的推广，搜索引擎的友好
- 前言许多任务程序如果为其构造为一个命令行界面，就可以通过接受不同的参数来改变它的工作方式。例如，在爬虫程序中，不同 URL&nbs
- 首先让我们看下 YUI 是如何处理的：var toObject = function(a) { var o = {
