ms15r23.m

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

%用内插公式重建信号
fm=1;
dt=0.1;
t=-4:dt:4;
x=sinc(t); 
subplot(4,1,1),plot(t,x,'k');grid
axis([min(t) max(t) 1.1*min(x) 1.1*max(x)]);
title('用时域卷积重建抽样信号');
for i=1:3;
    fs=i*fm;Ts=1/fs;         %确定采样频率和周期      
     n=-4:Ts:4; 
     t1=-4:Ts:4;   
     x1=sinc(n*fm/fs);              %生成抽样信号
     T_N=ones(length(n),1)*t1-n'*Ts*ones(1,length(t1));%生成t-nT矩阵
    xa=x1*sinc(fs*pi*T_N);                      %内插公式
    subplot(4,1,i+1),plot(n,xa,'k');grid
    axis([min(n) max(n) 1.1*min(xa) 1.1*max(xa)]);
end