python怎么解微分方程组_python能解微分方程吗

python怎么解微分方程组_python能解微分方程吗SympySympy是一个数学符号库(sym代表了symbol,符号),包括了积分,微分方程等各种数学运算方法,为python提供了强大的数学运算支持。对于图像来说,虽然都是做离散的计算,操作最多的还是numpy里的数组,但实际上,这个库包含了积分微分,三角等最基本的数学运算,可以说是工科最基本的,用起来媲美matlab。sympy安装sudopipinstallsympy或者sudo…

Sympy

Sympy是一个数学符号库(sym代表了symbol,符号),包括了积分,微分方程等各种数学运算方法,为python提供了强大的数学运算支持。对于图像来说,虽然都是做离散的计算,操作最多的还是numpy里的数组,但实际上,这个库包含了积分微分,三角等最基本的数学运算,可以说是工科最基本的,用起来媲美matlab。

python怎么解微分方程组_python能解微分方程吗

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能解微分方程吗

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。
如需转载请保留出处:https://bianchenghao.cn/38339.html

(0)
编程小号编程小号

相关推荐

发表回复

您的电子邮箱地址不会被公开。 必填项已用*标注