代码：
% Discrete-time Signal x1(n)
% Ts = 0.0002; n = -25:1:25; nTs = n*Ts; Fs = 1/Ts; x = exp(-1000*abs(nTs));
Ts = 0.001; n = -5:1:5; nTs = n*Ts; Fs = 1/Ts; x = exp(-1000*abs(nTs)); % Analog Signal
Dt = 0.00005; t = -0.005:Dt:0.005;
xa = x * sinc(Fs*(ones(length(n),1)*t - nTs'*ones(1,length(t)))) ; %% --------------------------------------------------------------------
%%
%% --------------------------------------------------------------------
figure('NumberTitle', 'off', 'Name', '<<DSP MATLAB>> Example3.21');
set(gcf,'Color','white');
subplot(1,1,1); plot(t*1000,xa); grid on; %axis([0,1,0,1.5]);
% title('Reconstructed Signal from x1(n) using sinc function');
title('Reconstructed Signal from x2(n) using sinc function');
xlabel('t in msec units.'); ylabel('xa(n)'); hold on;
% stem(n*Ts*1000,x); gtext('Ts=0.2 msec'); hold off;
stem(n*Ts*1000,x); gtext('Ts=1 msec'); hold off; % subplot(2,1,2); plot(w/pi, X); grid on; % axis([-1,1,-1,1]);
% title('Discrete-time Fourier Transform');
% xlabel('frequency in \pi units'); ylabel('X1(w)'); %% -------------------------------------------------------------------
%%
%% ------------------------------------------------------------------- % check
error = max(abs(xa - exp(-1000*abs(t))))
运行结果：
