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]])

以上这篇python/sympy求解矩阵方程的方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持脚本之家。

python解一元三次方程_python/sympy求解矩阵方程的方法相关推荐

  1. python计算矩阵方程_python/sympy求解矩阵方程的方法

    sympy版本:1.2 假设求解矩阵方程 AX=A+2X 其中 求解之前对矩阵方程化简为 (A−2E)X=A 令 B=(A−2E) 使用qtconsole输入下面程序进行求解 In [26]: fro ...

  2. python解非线性方程组_python scipy求解非线性方程的方法(fsolve/root)

    使用scipy.optimize模块的root和fsolve函数进行数值求解线性及非线性方程,下面直接贴上代码,代码很简单 from scipy.integrate import odeint imp ...

  3. 【python求一元三次方程】

    python求一元三次方程 # 方程求解 import sympy as sy import sympy import sympy as sp # 导入sympy包 import matha=math ...

  4. 问题三十七:C++怎么解一元四次方程?(2)——怎么解一元三次方程

    37.2 怎么解一元三次方程? 用"盛金公式"求解: ----------------------------------------------main.cpp -------- ...

  5. 盛金公式解一元三次方程_一元三次方程解法(卡尔丹公式法盛金公式法)

    卡尔丹公式法 特殊型一元三次方程 X^3+pX+q=0 (p.q∈R) 判别式Δ=(q/2)2+(p/3)3 卡尔丹公式 X1=(Y1)(1/3)+(Y2)(1/3) X2= (Y1)(1/3)ω+( ...

  6. 解一元三次方程方程问题细节探讨

    解一元三次方程 虽然我是一个新手,按"规矩"应该要经常在博客里面把自己的错误都写出来做个总结,但是我比较保守,总是觉得博客是一个社交平台,一些相对不太成熟的东西还是写个word文档 ...

  7. Python解一元一次方程和一元二次方程,输出数学形式和小数形式

    Python解一元一次方程和一元二次方程,输出数学形式和小数形式 一元一次方程 一元二次方程 完整代码 运行结果 创建函数: def eqt(c,b,a=0): 其中,a.b.c分别是二次项系数.一次 ...

  8. python解复杂方程_Python数据处理篇之Sympy系列(五)---解方程

    前言 sympy不仅在符号运算方面强大,在解方程方面也是很强大. 本章节学习对应官网的:Solvers 官方教程 (一)求解多元一次方程-solve() 1.说明: 解多元一次方程可以使用solve( ...

  9. python解二元一次方程_Python数据处理篇之Sympy系列(五)---解方程

    前言 sympy不仅在符号运算方面强大,在解方程方面也是很强大. 本章节学习对应官网的:Solvers 官方教程 https://docs.sympy.org/latest/tutorial/solv ...

最新文章

  1. Lync与Exchange 2013 UM集成:Exchange 配置
  2. 前端代码规范-CSS
  3. 前后端完全分离出现跨域、无法访问后台解决方案
  4. 容斥原理应用(求1~r中有多少个数与n互素)
  5. java学生签到系统视频教程_手把手教你做一个Java web学生信息、选课、签到考勤、成绩管理系统附带完整源码及视频开发教程...
  6. C/C++遍历进程和进程ID的小工具
  7. 灰色预测模型matlab_Matlab数据分析,2020研究生报名人数灰色预测
  8. 【CodeForces】947 D. Picking Strings
  9. oracle正则表达式
  10. 淘宝中间件借阿里云对外 日均调千亿次
  11. java null转string_java null强转string
  12. 计算机4级学那些课程,计算机四级考试科目是什么
  13. windows cmd命令行查看结束进程
  14. “宅一族”大数据报告:喜欢学习,喜欢健身,生活同样很精彩
  15. 两位末代皇帝的宠物文史小语新浪博客
  16. Kafka(一):概述
  17. 国外客户最讨厌的开发信
  18. 数据分析系列--常见的评估方法
  19. GitHub也能CI/CD了 如何使用GitHub的Action?
  20. 创建deployment

热门文章

  1. webpack5 模块联邦 single-spa基座微前端 systemjs
  2. 【kotin】统计字符串中的字母数(拓展函数)
  3. 使用Pass提高效率
  4. 前端学习-flex布局
  5. latex 插图解释_仅使用一些插图和视频即可解释Big O符号
  6. 组件用.vue还是.js_如何使用Vue.js 2.0构建灵活的图像上传器组件
  7. lighthouse_如何提高您的渐进式网络应用的Google Lighthouse得分(最高可达100分)
  8. Java-时间复杂度和空间复杂度
  9. Python-关于正则表达式的总结
  10. Python的10个神奇的技巧