e319.m

matlab例子

Ts=0.0002;n=-25:1:25;Fs=1/Ts;nTs=n*Ts;
x=exp(-1000*abs(n*Ts));
Dt=0.00005;t=-0.005:Dt:0.005;
xa=x*sinc(Fs*(ones(length(n),1)*t-nTs'*ones(1,length(t))));
subplot(2,1,1);plot(t*1000,xa);hold on
stem(n*Ts*1000,x);hold off;
error=max(abs(xa-exp(-1000*abs(t))))