ch2example25prg1.m

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

% ch2example25prg1.m
% 采用较低速率仿真，并用样条插值使得结果精细化
simurate=3000;                          % 较低仿真速率
mysimopts = simset('Solver','ode5');    % 设置仿真求解器为ode5的
mysimopts = simset(mysimopts,'FixedStep',1./simurate); % 设置仿真步进
timespan=[0, 0.01];                     % 仿真时间段参数
sim('ch2example25',timespan,mysimopts); % 设置仿真参数，仿真时间段并执行仿真
plot(simout.time,simout.signals.values,'o-k');
hold on;axis([0 0.005 -2 2]);
t=timespan(1):1/30000:timespan(2);      % 插值时间序列
intp_y=interp1(simout.time,simout.signals.values,t,'spline');% 样条插值
plot(t,intp_y,'xr');                    % 画出插值结果
% 采用较高速率仿真
simurate=30000;                         % 较高仿真速率
mysimopts = simset(mysimopts,'FixedStep',1./simurate); % 设置仿真步进
sim('ch2example25',timespan,mysimopts); % 设置仿真参数，仿真时间段并执行仿真
plot(simout.time,simout.signals.values,'b');% 画出仿真结果
legend('较低速率仿真输出波形','样条插值精细化的结果','较高速率仿真输出波形');