// An highlighted block
import numpy as np
##定义未知量和待定参数amplitude和detuning
amplitude=20
detuning=20
#郎之万方程经典形式
def dfun(t,x,y,q,p):
dvec=np.zeros((4))
dvec[0]=-0.5*x-(detuning+0.02*q)*y+amplitude
dvec[1]=-0.5*y+(detuning+0.02*q)*x
dvec[2]=20*p
dvec[3]=0.02*(x**2+y**2)-20*q-0.01*p
return dvec
#四阶龙格库塔:dfun=微分方程;t为自变量,x,y,q,p为因变量;求解区间(a,b);x0,y0,q0,p0=初值,h=步长
def main(dfun,a,b,t,x,y,q,p,x0,y0,q0,p0,h):
K1=dfun(t,x,y,q,p)[0]
L1=dfun(t,x,y,q,p)[1
基于python的四阶龙格库塔求解一阶常微分方程组的问题
最新推荐文章于 2023-03-20 17:03:54 发布
![](https://img-home.csdnimg.cn/images/20240711042549.png)