% ch2example1prg25.m
% 采用較低速率仿真，并用樣條插值使得結(jié)果精細(xì)化
simurate=3000;                          % 較低仿真速率
mysimopts = simset('Solver','ode5');    % 設(shè)置仿真求解器為ode5的
mysimopts = simset(mysimopts,'FixedStep',1./simurate); % 設(shè)置仿真步進(jìn)
timespan=[0, 0.01];                     % 仿真時間段參數(shù)
sim('ch2example25',timespan,mysimopts); % 設(shè)置仿真參數(shù)，仿真時間段并執(zhí)行仿真
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');                    % 畫出插值結(jié)果
% 采用較高速率仿真
simurate=30000;                         % 較高仿真速率
mysimopts = simset(mysimopts,'FixedStep',1./simurate); % 設(shè)置仿真步進(jìn)
sim('ch2example25',timespan,mysimopts); % 設(shè)置仿真參數(shù)，仿真時間段并執(zhí)行仿真
plot(simout.time,simout.signals.values,'b');% 畫出仿真結(jié)果
legend('較低速率仿真輸出波形','樣條插值精細(xì)化的結(jié)果','較高速率仿真輸出波形');