Question

用matlab解二元二阶微分方程（用ode45）

有一个一元的 但是二元的我不会了
下面是那个一元的





Answer 1

function dy=dong22(t,y)

m=1;

k=2;

g=    %注意输入g l0

l0=

dy=[y(2);

    y(1)*y(4)^2+g*cos(y(3))-(k/m)*(y(1)-l0);

    y(4);

    -(2/y(1))*y(2)*y(4)-(g/y(1))*sin(y(3))];

    然后和一次微分方程的一样。

追问

那我在命令框里面输入
y0=[1;0;5*3.14/180;0];
[t,y]=ode45('dong22',[0,200],y0];
这里就开始说有错误了
还有一句还没输入
plot（t,y(:,1))
请帮我看一下哪里错误了 谢谢

追答

function dy=dong22(t,y)
里面的g和l0要给出的·。

追问

我已经给出了，可不可以加我qq373605738 谢谢 感激不尽

追答

已经加了。