如何Matlab解三元二阶非线性微分方程组？

发布网友发布时间：2022-03-25 19:59

共1个回答

热心网友时间：2022-03-25 21:28

三元二阶非线性微分方程组用Matlab解，可以按下列思路求解。

首先，自定义微分方程组函数，其内容：

dy(1)=x(2);%降阶

dy(2)=0.6615-(0.3*x(2)+2848*x(1)+0.0025*x(3)+1.227e3*(x(1)^3+x(1))/(x(1)^2-1e-4))/0.0675

dy(3)=(0.025*x(2)-x(3)/2e5)/2.135

其二，确定初值，x0=[0,0,0];

其三，确定时间变量范围，t【0，10】

其四，用ode45函数求解。其命令格式

[t,x]=ode45(@func,tspan,x0)

其五，用plot函数绘制t—x(t)、t—dx/dt、t—V(t)

执行结果，可以发现，该方程似乎为振动型方程，其x的振幅为±3.8，Vx速度为±800，Vx速度为±4e-3

追问能不能把程序发一下谢谢

追答可以的，已发私信请查收。