如何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
追问能不能把程序发一下谢谢追答可以的，已发私信请查收。
