ch1example2prg2.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

% 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