ch1example1prg1.m

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

% ch1example1prg1.m
g=9.8;  % 重力加速度
v=0;    % 设定初始速度条件
s=0;    % 设定初始位移条件
t=0;    % 设定起始时间
dt=0.1; % 设置计算步长
N=20;   % 设置仿真递推次数. 仿真时间等于N与dt的乘积
for k=1:N
  v=v+g*dt;                 % 计算新时刻的速度
  s(k+1)=s(k)+v*dt;         % 新位移
  t(k+1)=t(k)+dt;           % 时间更新
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,s,'o', t_theory,s_theory, '-');
xlabel('时间 t'); ylabel('位移 s');
legend('仿真结果','理论结果');