q711.m

MATLAB及在电子信息课程中的应用例题索引与程序

%《MATLAB在电子信息课程中的应用》第七章例7.11程序q711
% 时域采样与重构
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clear; close all;
A=444.128; a=50*sqrt(2)*pi; b=a;
for k=1:2
   if k==1 Fs=400;
   elseif k==2 Fs=1000;end
   T=1/Fs;   dt=T/3;   %每个采样间隔T上g(t)取三个样点
   Tp=0.03;  				% 重构时间区间为[0,0.03s]
   t=0:dt:Tp;				% 生成序列t
   n=0:Tp/T;						% 生成序列n
   TMN=ones(length(n),1)*t- n'*T*ones(1,length(t));	% 生成TNM矩阵
   x=A*exp(-a.*n*T).*sin(b*n*T);		% 生成模拟信号采样序列x(n)
   xa=x*sinc(Fs*TMN);  					%　内插公式
   subplot(2,1,k);plot(t,xa);hold on
   axis([0,max(t),min(xa)-10,max(xa)+10]);
   st1=sprintf('由Fs= %d',Fs);st2='Hz的采样序列x(n)重构的信号';
   ylabel('xa(t)');
   st=[st1,st2];title(st)
   xo=A*exp(-a.*t).*sin(b*t);		%以3Fs对原始模拟信号采样
   stem(t,xo,'.');line([0,max(t)],[0,0])
   emax2=max(abs(xa-xo))
end