python怎么解微分方程组_python能解微分方程吗
Sympy
Sympy是一个数学符号库(sym代表了symbol,符号),包括了积分,微分方程等各种数学运算方法,为python提供了强大的数学运算支持。对于图像来说,虽然都是做离散的计算,操作最多的还是numpy里的数组,但实际上,这个库包含了积分微分,三角等最基本的数学运算,可以说是工科最基本的,用起来媲美matlab。
sympy安装
sudo pip install sympy 或者 sudo pip3 install sympy
如果你用anaconda则 conda install sympy
相关推荐:《python视频》
odeint()函数
odeint()函数是scipy库中一个数值求解微分方程的函数
odeint()函数需要至少三个变量,第一个是微分方程函数,第二个是微分方程初值,第三个是微分的自变量。
例子:#y"+a*y'+b*y=0
from scipy.integrate import odeint #使用odeint之前,首先从scipy.integrate中调用它from pylab import *
def deriv(y,t): # 返回值是y和y的导数组成的数组
a = -2.0
b = -0.1
return array([ y[1], a*y[0]+b*y[1] ])
time = linspace(0.0,50.0,1000)
yinit = array([0.0005,0.2]) # 初值
y = odeint(deriv,yinit,time)
figure()
plot(time,y[:,0],label='y') #y[:,0]即返回值的第一列,是y的值。label是为了显示legend用的。
plot(time,y[:,1],label="y'") #y[:,1]即返回值的第二列,是y’的值
xlabel('t')
ylabel('y')
legend()
show()
输出结果如下:
python怎么解微分方程组_python能解微分方程吗相关推荐
- 用python怎样解偏微分方程组_Python能解偏微分方程吗
fipy 菲克定律是指在不依靠宏观的混合作用发生的传质现象时,描述分子扩散过程中传质通量与浓度梯度之间关系的定律.菲克定律是阿道夫·菲克(Adolf Fick)于1855年提出. 由菲克第二定律可以得 ...
- 龙格库塔公式法解微分方程组初值问题实例
用四阶龙格库塔公式(取h=0.1)解下列微分方程组初值问题: { y 1 ′ = 120 − 2 y 1 + 2 y 2 ( 0 ≤ x ≤ 1.0 ) , y 2 ′ = 2 y 1 − 5 y 2 ...
- matlab解二阶微分方程组ode,MATLAB解含参数方程、矩阵方程、二阶微分方程组
1.如下公式: 其中x.y.m.n为参数,a.b为未知数,利用MATLAB求解方程: syms x y m n a b; [a,b]=solve('x=m*cos(a)+n*cos(a+b)','y= ...
- 第二十五讲 用线性代数解微分方程组
一,上一讲的例题,如图: 设x=T1x=T_1x=T1,y=T2y=T_2y=T2 方程组为·:{x′=−2x+2yy′=2x−5y\left\{\begin{matrix}{x}'=-2x+2y ...
- python求解微分方程组_python – SymPy / SciPy:求解具有不同变量的常微分方程组...
我是SymPy和Python的新手,我目前正在使用Python 2.7和SymPy 0.7.5,其目标是: a)从文本文件中读取微分方程组 b)解决系统问题 我已经阅读了this question和t ...
- python怎么算一元二次方程_python如何解一元二次方程
本文收集整理关于python如何解一元二次方程的相关议题,使用内容导航快速到达. 内容导航: Q1:python 使用try except解一元二次方程问题 123456789101112131415 ...
- python求解微分方程组_python – SymPy / SciPy:求解具有不同变量的常微分方程组
我是SymPy和 Python的新手,我目前正在使用Python 2.7和SymPy 0.7.5,其目标是: a)从文本文件中读取微分方程组 b)解决系统问题 我已经阅读了this question和 ...
- 龙格库塔法解微分方程组的matlab程序,MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc...
MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc MATLAB实例源码教程龙格库塔法求解微分方程组源代码实例题目用经典 Runge-Kutta方法求下列一阶微分方程组的近似解y1 ...
- matlab解微分方程组_MATLAB编程入门 求解常微分方程 通解 特解 数值解
一.求微分方程的通解 clear all; f1=dsolve('Dy-y=sin(x)') %默认情况下以t为自变量求解 求出的C2为积分常数 f2=dsolve('Dy-y=sin(x ...
- python求解方程组_python如何解方程组
解方程的基本思想是:aa = solve([f],[x]) f是方程,x是变量,这个代码解的是关于x的方程f=0 既然可以解方程(组),就一定可以解微分方程(组).from sympy import ...
最新文章
- dockerfile构建镜像的命令_编写Dockerfile的最佳实践
- java保存数字_Java:JFormattedTextField保存数字
- i2c--ioctl--主机控制器驱动(i2c_adapter)--外设驱动(i2c_driver)
- eclipse mybatis插件_Archetype插件的介绍和使用
- 第四十三节,文件、文件夹、压缩包、处理模块shutil
- ubuntu10.04下audacious2.4源码编译过程(解决2.3cue的bug)
- java hanlp分词_Hanlp分词实例:Java实现TFIDF算法
- kendo Grid json解析的问题
- openeim再被刺也忍着痛继续吃
- 骁龙8性能巅峰旗舰!黑鲨5系列获3C认证:配备120W快充
- VS2010中添加WebService注意的几个地方
- 西门子精彩SMART触摸屏下载程序时常见问题分析与处理对策
- 用phpstudy安装DVWA
- Python添加或去掉百分号
- Linuxptp安装部署
- 2012人类将从“微博体”过渡到“微媒体”。至于你信不信,反正我信了
- Excel 里筛选手机号码所属运营商
- 2020全国普通高校大学生竞赛排行榜出炉!
- 鼠标灵敏度调多少合适
- 创业者必学的六种植物[转]