最近在写毕业论文,用到了龙格库塔法,可是函数怎么使用遇到了困惑,想求助下大家帮忙,谢谢大家。
下面是 代码
:
tspan = [1,20000]; %%仿真时间
y0 = [0;0;0.1;2.5]; %%微分方程初始值
[t,y] = ode45(@equation,tspan,y0);
figure(1);plot(t,y(:,1));
xlabel('t');ylabel('θ');title('系绳的偏角');grid on
figure(2);plot(t,y(:,2));
xlabel('t');ylabel('ω');title('系绳偏角的导数');grid on
figure(3);plot(t,y(:,3));
xlabel('t');ylabel('L');title('系绳绳长变化');grid on
figure(4);plot(t,y(:,4));
xlabel('t');