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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

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

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

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