python的scipy里的odeint这个求微分方程的函数怎么用啊 scipy.integrate.odeint(func,y0,t,args=(),Dfun=None,col_deriv=0,full_output=0,ml=None,mu=None,rtol=None,atol=None,tcrit=None,h0=0.0,hmax=0.0,hmin=0.0,ixpr=0,mxstep=0,mxhnil=0,mxordn=12,mxords=5,printmessg=0)实际使用中,还是主要使用前三个参数,即微分方程的描写函数、初值和需要求解函数值对应的的时间点。接收数组形式。这个函数,要求微分方程必须化为标准形式,即dy/dt=f(y,t,)。fromscipyimportodeint y=odeint(dy/dt=r*y*(1-y/k),y(0)=0.1,t)对于微分方程全还给老师了,python的scipy里的odeint这个求微分方程的函数怎么用啊? odeint实际是用来解微分方程组的。令z=y’,可以把方程化为方程组:y'=zz'=-b*y-a*z将y,z分别替换为y[0],y[1]就得到了程序里的returnarray([y[1],a*y[0]+b*y[1]]),这个程序把a,b前面的符号放到参数赋值里了。搞明白上面的就能用来解方程了,下面画图部分无关紧要hold('on')是用来保持之前画的曲线legend()显示曲线的标签python中如何求微分方程的特解?谢谢 #y\"+a*y'+b*y=0from scipy.integrate import odeintfrom pylab import*def deriv(y,t):#返回值是y和y的导数组成的数组a=-2.0b=-0.1return 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的sympy模块解高阶微分方程,ytho是一个很有趣的语言,可以在命令行窗口运行。ytho中有很多功能强大的模块,这篇经验告诉你,如何使用ytho的ymy模块求解高阶。求助关于一个用Python求微分方程并算出极值 所说所有的变量都是对象。对象在python里,其实是一个指针,指向一个数据结构,数据结构里有属性,有方法。对象通常就是指变量。从面向对象OO的概念来讲,对象是类的一个。如何使用python计算常微分方程? 常用形式odeint(func,y0,t,args,Dfun)一般这种形式就够用了。下面是官方的例子,求解的是D(D(y1))-t*y1=0为了方便,采取D=d/dt。如果我们令初值y1(0)=1.0/3*(2.0/3.0)/gamma(2.0/3.0)D(y1)(0)=-1.0/3*(1.0/3.0)/gamma(1.0/3.0)这个微分方程的解y1=airy(t)。令D(y1)=y0,就有这个常微分方程组。D(y0)=t*y1D(y1)=y0Python求解该微分方程。from scipy.integrate import odeintfrom scipy.special import gamma,airyy1_0=1.0/3*(2.0/3.0)/gamma(2.0/3.0)y0_0=-1.0/3*(1.0/3.0)/gamma(1.0/3.0)y0=[y0_0,y1_0]def func(y,t):return[t*y[1],y[0]]def gradient(y,t):return[[0,t],[1,0]]x=arange(0,4.0,0.01)t=xychk=airy(x)[0]y=odeint(func,y0,t)y2=odeint(func,y0,t,Dfun=gradient)print ychk[:36:6][0.355028 0.339511 0.324068 0.308763 0.293658 0.278806]print y[:36:6,1][0.355028 0.339511 0.324067 0.308763 0.293658 0.278806]print y2[:36:6,1][0.355028 0.339511 0.324067 0.308763 0.293658 0.278806]得到的解与精确值相比,误差相当小。args是额外的参数。用法请参看下面的例子。这是一个洛仑兹曲线的求解,并且用matplotlib绘。有没有专门画微分方程相图或方向场的软件,或者用Python如何实现? 有一个三变量(满足$x+y+z=1$)的微分方程组,有没有专门画相图或方向场的软件呢?另,我用python画了图…python fipy 求解偏微分方程 jingyan.baidu.com/article/46650658d71872f548e5f852.html 热量的传递由微分方程给出,可以简单地理解为热量按照温度降低最快的方向进行传递。公式右边是温度的梯度,左边。
随机阅读
- 诸位大神跪求答案!!靠谱卧式蒸汽锅炉厂家价格,卧式蒸汽锅炉上档次么?? 节能卧式常压蒸气锅炉
- 镇巴永乐镇领导班 石泉自驾游攻略?
- 泰迪狗几个月剃胎毛好 刚剃光毛吃什么毛发好 泰迪狗几个月毛能长的好看
- 野战军歼敌数 国共内战死亡人数谁知道啊
- 火山版水中花 水中花这首歌词是什么意思
- 口袋妖怪白金光水之石
- 南京市汉中路155号离夫子庙远吗 南京市中山东路305号到建邺区汉中路155号远吗?
- 赵丽颖参加过什么综艺节目? 高梓淇赵丽颖参加的综艺节目
- 河北衡水市橡胶城 衡水市区橡胶公司
- 大连到瓦房店高速几个收费站 从阜新到瓦房店走高速收费得用多少钱
- 哪个牌子的花生好吃? 锦州贵香源
- 增值税一般纳税人分别有:17%、11%、6%、3%税率分别适用哪些行业? 增值税的纳税人有哪些
- 山楂罐头木耳 银耳和山楂能在一起做吗?
- 从南昌火车站该坐几路公交车到南昌大桥呢? 南昌大桥到火车站
- 《爱与自由》读后感 阿马蒂亚森 《理性与自由》读后感
- 革命圣地延安在哪个省 延安市圣地路家和酒店
- 二硫化钼是否属于危险化学品 三氧化钼pH值多少
- 张亮麻辣烫(新天地店)怎么样 张亮麻辣烫怎么样
- 停止进程服务
- 小龙大肉面和路子铺 南京现在有什么网红面馆?