ch1example2prg1.m

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

% 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]);