python解一元三次方程_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]])
以上这篇python/sympy求解矩阵方程的方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持脚本之家。
python解一元三次方程_python/sympy求解矩阵方程的方法相关推荐
- python计算矩阵方程_python/sympy求解矩阵方程的方法
sympy版本:1.2 假设求解矩阵方程 AX=A+2X 其中 求解之前对矩阵方程化简为 (A−2E)X=A 令 B=(A−2E) 使用qtconsole输入下面程序进行求解 In [26]: fro ...
- python解非线性方程组_python scipy求解非线性方程的方法(fsolve/root)
使用scipy.optimize模块的root和fsolve函数进行数值求解线性及非线性方程,下面直接贴上代码,代码很简单 from scipy.integrate import odeint imp ...
- 【python求一元三次方程】
python求一元三次方程 # 方程求解 import sympy as sy import sympy import sympy as sp # 导入sympy包 import matha=math ...
- 问题三十七:C++怎么解一元四次方程?(2)——怎么解一元三次方程
37.2 怎么解一元三次方程? 用"盛金公式"求解: ----------------------------------------------main.cpp -------- ...
- 盛金公式解一元三次方程_一元三次方程解法(卡尔丹公式法盛金公式法)
卡尔丹公式法 特殊型一元三次方程 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)ω+( ...
- 解一元三次方程方程问题细节探讨
解一元三次方程 虽然我是一个新手,按"规矩"应该要经常在博客里面把自己的错误都写出来做个总结,但是我比较保守,总是觉得博客是一个社交平台,一些相对不太成熟的东西还是写个word文档 ...
- Python解一元一次方程和一元二次方程,输出数学形式和小数形式
Python解一元一次方程和一元二次方程,输出数学形式和小数形式 一元一次方程 一元二次方程 完整代码 运行结果 创建函数: def eqt(c,b,a=0): 其中,a.b.c分别是二次项系数.一次 ...
- python解复杂方程_Python数据处理篇之Sympy系列(五)---解方程
前言 sympy不仅在符号运算方面强大,在解方程方面也是很强大. 本章节学习对应官网的:Solvers 官方教程 (一)求解多元一次方程-solve() 1.说明: 解多元一次方程可以使用solve( ...
- python解二元一次方程_Python数据处理篇之Sympy系列(五)---解方程
前言 sympy不仅在符号运算方面强大,在解方程方面也是很强大. 本章节学习对应官网的:Solvers 官方教程 https://docs.sympy.org/latest/tutorial/solv ...
最新文章
- Lync与Exchange 2013 UM集成:Exchange 配置
- 前端代码规范-CSS
- 前后端完全分离出现跨域、无法访问后台解决方案
- 容斥原理应用(求1~r中有多少个数与n互素)
- java学生签到系统视频教程_手把手教你做一个Java web学生信息、选课、签到考勤、成绩管理系统附带完整源码及视频开发教程...
- C/C++遍历进程和进程ID的小工具
- 灰色预测模型matlab_Matlab数据分析,2020研究生报名人数灰色预测
- 【CodeForces】947 D. Picking Strings
- oracle正则表达式
- 淘宝中间件借阿里云对外 日均调千亿次
- java null转string_java null强转string
- 计算机4级学那些课程,计算机四级考试科目是什么
- windows cmd命令行查看结束进程
- “宅一族”大数据报告:喜欢学习,喜欢健身,生活同样很精彩
- 两位末代皇帝的宠物文史小语新浪博客
- Kafka(一):概述
- 国外客户最讨厌的开发信
- 数据分析系列--常见的评估方法
- GitHub也能CI/CD了 如何使用GitHub的Action?
- 创建deployment
热门文章
- webpack5 模块联邦 single-spa基座微前端 systemjs
- 【kotin】统计字符串中的字母数(拓展函数)
- 使用Pass提高效率
- 前端学习-flex布局
- latex 插图解释_仅使用一些插图和视频即可解释Big O符号
- 组件用.vue还是.js_如何使用Vue.js 2.0构建灵活的图像上传器组件
- lighthouse_如何提高您的渐进式网络应用的Google Lighthouse得分(最高可达100分)
- Java-时间复杂度和空间复杂度
- Python-关于正则表达式的总结
- Python的10个神奇的技巧