Ⅰ 三自由度质量弹簧系统,用matlab解这个函数,ode45或4阶龙格——库塔法都可以。谢谢
[T,Y]=ode45(@fun,[01],[351246]);
plot(T,Y(:,1))
holdon
plot(T,Y(:,1),'-ro')
plot(T,Y(:,2),'-g*')
plot(T,Y(:,3),'-ys')
plot(T,Y(:,4),'-gs')
plot(T,Y(:,5),'--b*')
plot(T,Y(:,6),'--ko')
legend('Y(:,1)','Y(:,2)','Y(:,3)','Y(:,4)','Y(:,5)','Y(:,6)')
函数fun.m见附件