fg_03_12.m

英文书《Digital Signal Processing with Examples in MATLAB》附带的MATLAB实例

% fg_03_12  samples and waveforms.
N=8;
t=[0:2:14];
tt=linspace(0,10,500);
c=[4 5 5 6];
x=[1 c 6 7 1000; 8 c 1 0 -1000; 3 c 5 5 5];
y=zeros(3,500);
% Construct 3 waveforms.
for k=1:3
   p=polyfit(t,x(k,:),N);
   y(k,:)=polyval(p,tt);
end
% Plot
sp_fig(1);
axis([0 7 0 10]);
set(gca,'fontname','times','fontsize',16)
set(gca,'position',[.13 .15 .8 .6]);
line(t(2:5)-2,c,'linestyle','none','marker','o','color','k');
for k=1:3
   line(tt-2,y(k,:),'color','k'); grid;
end
xlabel('t'); ylabel('x(t)');
title('T=2');