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
投稿猜你喜欢
之前的文章介绍了python抓取网页数据并将数据保存到本地excel文件，后续可以将数据保存到数据库（SqlServer、mysql等）中，- 此货很干,跟上脚步!!!Cookiecookie是什么东西?小饼干?能吃吗?简单来说就是你第一次用账号密码访问服务器服务器在你本机硬盘上设置
- 自动发送邮件功能是我们经常要用到的，比如每天定时统计报表信息，然后自动发送给运营人员，协助运营人员进行业务数据分析。本文是用Python写的
- 对其中的参数进行解释plt.subplots(figsize=(9, 9))设置画面大小，会使得整个画面等比例放大的sns.heapmap(
- 前言表是数据库存储数据的基本单位。一个表包含若干个字段或记录。表的操作包括创建新表、修改表和删除表。这些操作都是数据库管理中最基本，也是最重
- 在没介绍正文之前先给大家补充点go语言基本知识及实例。Go 语言教程Go 是一个开源的编程语言，它能让构造简单、可靠且高效的软件变得容易。G
- 如果当前绝对定位的元素需要透明（没啥内容、且不设置背景），背景元素有内容透出来的时候，IE6/IE7响应的不是期望的当前元素，而是背景元素。
- 本文实例为大家分享了PHP实现简单注册登录系统的具体代码，供大家参考，具体内容如下目录结构如下，其中function文件夹下包含两个函数文件
- 一、Array数组的更改array数组中修改的API示例如下：//创建大量相同元素的数组//创建有10个String类型元素的数组，并且每个
- String含义：String是一个封装char[]数组的对象,字符串不可变String str = “abc&
- 简而言之就是，nn.Sequential类似于Keras中的贯序模型，它是Module的子类，在构建数个网络层之后会自动调用forward(
- Tuple 叫做 tuple，用小括号、或者无括号来表述，是一连串有顺序的数字。a_tuple = (12, 3, 5, 15 , 6)an
- Python实现对网易云音乐的数据进行一个数据清洗和可视化分析对音乐数据进行数据清洗与可视化分析关于数据的清洗，实际上在上一一篇文章关于抓取
- cgi.h #ifndef CGI_H#define CGI_H#include <stdio.h>#include <s
- 1、截取GB2312字符用的函数 PHP代码 <?php //截取中文字符串 function mysubstr($str, $sta
- torchvision包 包含了目前流行的数据集，模型结构和常用的图片转换工具。torchvision.datasets中包含了以下数据集M
- 本文实例讲述了PHP缓存集成库phpFastCache用法。分享给大家供大家参考。具体分析如下：phpFastCache是一个开源的PHP缓
- 前言之前一直写不出来，这周周日花了一下午终于弄懂了， 顺便放博客里，方便以后忘记了再看看。 要实现的是输入一张 图，起点，终点，输出起点和终
- 问题你想反方向迭代一个序列解决方案使用内置的 reversed() 函数，比如：>>> a = [1, 2, 3, 4]&
- 目录prometheus通过exporter监控mysql，并用grafana图表展示1、测试机器 2、配置mysql host0
