ch1example2prg2.m

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

% ch1example2prg2.m
g=9.8;      % 重力加速度
v0=0;       % 初始速度
y0=1;       % 初始位置
m=1;        % 小球质量
t0=0;       % 起始时间
K=0.85;     % 弹跳的损耗系数
N=5000;     % 仿真的总步进数
dt=0.005;   % 仿真步长
v=v0;       % 初状态
y=y0;   
for k=1:N 
    if y >0     % 小球在空中的动力方程计算
        v =v -g*dt;
        y =y +v*dt;
    else        % 碰击瞬间的计算
        y =-K.*v*dt;  
        v =-K.*v-g*dt;   
    end 
    plot(0,y,'o'); 
    axis([-2 2 0 1]);               % 坐标范围固定
    set(gcf,'DoubleBuffer','on');   % 双缓冲避免作图闪烁
    drawnow;
end