ms15_3.m

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

%2个正弦叠加信号
f0=1;T0=1/f0;dt=0.01;
fm=5*f0;Tm=1/fm;
t=0:dt:3*T0;
x=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t);   %建立原连续信号
subplot(4,1,1),plot(t,x,'k');
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=0:(3*T0)/Ts 
  t1=0:Ts:3*T0;
  x1=sin(2*pi*n*f0/fs)+1/3*sin(6*pi*n*f0/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(t1,xa,'k');
  axis([min(t1) max(t1) 1.1*min(xa) 1.1*max(xa)]);
end