%?ch1problem2.m
g=9.8;???%?重力加速度
k=-1;????%?空氣阻力系數(shù)
dt=0.1;??%?設(shè)置計(jì)算步長
N=30;????%?設(shè)置仿真遞推次數(shù).?仿真時(shí)間等于N與dt的乘積
for?m=[1?2?10]??%?三種落體質(zhì)量
????v=0;????%?設(shè)定初始速度條件
????s=0;????%?設(shè)定初始位移條件
????t=0;????%?設(shè)定起始時(shí)間
????for?i=1:N
??????a=g+k/m*v;????????????????%?計(jì)算加速度
??????v=v+a*dt;?????????????????%?計(jì)算新時(shí)刻的速度
??????s（i+1）=s（i）+v*dt;?????????%?新位移
??????t（i+1）=t（i）+dt;???????????%?時(shí)間更新
????end
plot（ts‘o‘）;
hold?on;
end
%?理論計(jì)算?以便與仿真結(jié)果對照
t_theory=0:0.01:N*dt;???????%?設(shè)置解析計(jì)算的時(shí)間點(diǎn)
v_theory=g*t_theory;????????%?解析計(jì)算的瞬時(shí)速度
s_theory=1/2*g*t_theory.^2;?%?解析計(jì)算的瞬時(shí)位移
%?作圖:?仿真結(jié)果與解析結(jié)果對比
t=0:dt:N*dt;
plot（t_theorys_theory?‘-‘）;
xlabel（‘時(shí)間?t‘）;?ylabel（‘位移?s‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????805??2009-08-05?12:01??ch1problem2.m

-----------?---------??----------?-----??----

??????????????????805????????????????????1