ch1example2prg1.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 26 行

% ch1example2prg1.m
g=9.8;      % 重力加速度
v0=0;       % 初始速度
y0=1;       % 初始位置
m=1;        % 小球质量
t0=0;       % 起始时间
K=0.85;     % 弹跳的损耗系数
N=5000;     % 仿真的总步进数
dt=0.001;   % 仿真步长
v=v0;       % 初状态
y=y0;       
for k=1:N 
    if (y>0)|(v>0)    % 小球在空中的(含刚刚弹起瞬间)动力方程计算
        v =v -g*dt;
        y =y +v*dt;
    else              % 碰击瞬间的计算
        y =y-K.*v*dt;  
        v =-K.*v-g*dt;   
    end 
    s(k)=y;     % 将当前位移记录到s数组中以便作图
end
t=t0:dt:dt*(N-1);% 仿真时间长度
plot(t,s);
xlabel('时间 t');
ylabel('位移 y(t)');
axis([0 5 0 1.1]);