ms15r21.m

《数字信号处理实验(MATLAB版)》这

dt=0.1;
t=-4:dt:4;
f=sinc(t);                %建立原连续信号
subplot(4,1,1),plot(t,f,'k');
title('原连续信号');
Tm=1;fm=1/Tm;             %由原连续信号周期,计算最高频率
fs=3*fm;Ts=1/fs;          %取fs>2fm
n1=-4:Ts:4;
f1=sinc(n1);
subplot(4,1,2),stem(n1,f1,'filled','k');
title('fs>2fm抽样信号');
fs=2*fm;Ts=1/fs;          %取fs=2fm
n2=-4:Ts:4;
f2=sinc(n2);
subplot(4,1,3),stem(n2,f2,'filled','k');
title('fs=2fm抽样信号');
fs=1.2*fm;Ts=1/fs;        %取fs<2fm
n3=-4:Ts:4;
f3=sinc(n3);
subplot(4,1,4),stem(n3,f3,'filled','k');
title('fs<2fm抽样信号');