Question

matlab数值法求解三元一阶微分方程组

如图所示的方程，参数可以随便带入，求出三个未知量的数值解并画图（自变量是t），求matlab程序代码

Answer 1

这是一个常规问题, 给出代码参考吧:

m=1; b=1; r=1; g=9.8; w=1; F0=1; % 设为自己需要的数.

fun=@(t,y) [y(2);b/m*y(2)-g/r*sin(y(1))+F0*cos(y(3))/m/r;w];

[t,y]=ode45(fun,[0,10],[1,1,1]);

plot3(y(:,1),y(:,2),y(:,3))

grid on

xlabel \theta

ylabel y

zlabel z

追问

其中的y(:,1),y(:,2),y(:,3)分别对应三个变量？

追答

分别对应θ,y,z, 编程的代码中已经体现出来了.