Ⅰ 三自由度質量彈簧系統,用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見附件