sam_inverse_sub.m

很多matlab的源代码

% 此文件用于产生并绘制连续时间信号的抽样信号及由抽样信号重建原信号
% type=1 sinc函数内插
% type=2 ZOH
% type=3 FOH
% type=4 三次样条函数内插

Ts=1/Fs;                   % 抽样周期
Dt=0.00005;                  % 连续信号时间间隔 
t=-0.005:Dt:0.005;         % 连续信号时间坐标
xa=exp(-1000*abs(t));      % 连续时间信号
n=-0.005/Ts:1:0.005/Ts;    % 离散坐标
nTs=n*Ts;       
x=exp(-1000*abs(nTs));     % 抽样后的时间信号
subplot(2,1,1);
plot(t*1000,xa,'r')
hold on
stem(nTs*1000,x);
axis([-5,5,0,1]);
hold off;
legend('Continuous-time Signal','Descrete-time Signal');

subplot(2,1,2);
stem(nTs*1000,x);
axis([-5,5,0,1]);
hold on
switch type
case 1      % sinc函数内插
   x_inv=x*sinc(Fs*(ones(length(n),1)*t-nTs'*ones(1,length(t))));
   plot(t*1000,x_inv,'r')
case 2      % 零阶保持
   stairs(nTs*1000,x,'r');
case 3      % 一阶保持
   plot(nTs*1000,x,'r');
case 4      % 三次样条函数内插
   x_inv=spline(nTs,x,t);
   plot(t*1000,x_inv,'r');
end
legend('Descrete-time Signal','Reconstructed Signal');
hold off