MATLAB对简谐运动的分析

一、画图函数

t=0:pi/500:6;
x0=10sin(2t);
x1=10sin(2t+pi/3);
x2=10sin(2t+2pi/3);
plot(t,x0,’-r’);
hold on;
plot(t,x1,’-g’);
hold on;
plot(t,x2,’-b’);
hold on;
xlabel(‘t’);
ylabel(‘x’);
gtext('x0=10sin(2t)’);
gtext(‘x1=10sin(2t+pi/3)’);
gtext('x2=10sin(2t+2*pi/3)’)

**

二、运用for循环对步长进行检索

**
（最合适T=m.max步长，此处步长t=t的递增值0.1）*

T=10;
F1=0;F2=0;F3=0;F4=0;
t=0:0.1:20;
for m=1:2:7;
move=(4/pi)(1/m)sin((2pim/T)t);
F1=F1+move;
end;
for m=1:2:9;
move=(4/pi)(1/m)sin((2pim/T)t);
F2=F2+move;
end;
for m=1:2:19;
move=(4/pi)(1/m)sin((2pim/T)t);
F3=F3+move;
end;
for m=1:2:99;
move=(4/pi)(1/m)sin((2pi*m/T)*t);
F4=F4+move;
end;
plot(t,F1,‘y’,t,F2,‘r’,t,F3,‘b’,t,F4,‘g’)