example1_21.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 28 行

%例1-21程序，example1_21
%由x1恢复xa
Ts=0.0002;n=-25:1:25;nTs=n*Ts;
x1=exp(-1000*abs(nTs));
subplot(2,2,1);stem(nTs,x1)
title('x1');
%模拟信号恢复
Dt=0.00005;t=-0.005:Dt:0.005;
Fs=1/Ts;
xa=x1*sinc(Fs*(ones(length(n),1)*t-nTs'*ones(1,length(t))));
subplot(2,2,2);plot(t,xa)
title('由x1恢复xa');
%校验
error1=max(abs(xa-exp(-1000*abs(t))))

%由x2恢复xa
Ts=0.001;n=-5:1:5;nTs=n*Ts;
x=exp(-1000*abs(nTs));
subplot(2,2,3);stem(nTs,x)
title('x2');
%模拟信号恢复
Dt=0.00005;t=-0.005:Dt:0.005;
Fs=1/Ts;
xa=x*sinc(Fs*(ones(length(n),1)*t-nTs'*ones(1,length(t))));
subplot(2,2,4);plot(t,xa)
title('由x2恢复xa');
%校验
error2=max(abs(xa-exp(-1000*abs(t))))