ditongchouyang.m

低通抽样定理程序有课程设计的有需求的自己

clear all 
close all 
dt=0.01;   
t=0:dt:10; 
xt=0.1*cos(0.15*pi*t)+1.5*sin(2.5*pi*t)+0.5*cos(4*pi*t);
[f,xf]=T2F(t,xt); 
% 抽样信号，抽样速率为4Hz 
fs=4; 
sdt=1/fs; 
t1=0:sdt:10; 
st=0.1*cos(0.15*pi*t1)+1.5*sin(2.5*pi*t1)+0.5*cos(4*pi*t1);
 [f1,sf]=T2F(t1,st); 
% 恢复原始信号 
t2=-50:dt:50; 
gt=sinc(fs*t2); 
stt=sigexpand(st,sdt/dt); 
xt_t=conv(stt,gt);
figure(1) 
subplot(311); 
plot(t,xt); 
title('原始信号') 
subplot(312);
stem(t1,st); 
title('抽样信号') 
subplot(313); 
t3=-50:dt:60+sdt-dt; 
plot(t3,xt_t); 
title('抽样信号恢复') 
axis([0 10 -4 4]) 
*****************************************************  
   程序中用到二个函数，分别是：T2F.m和sigexpand.m。 
************************************************************
 function [f,sf]=T2F(t,st) 
dt=t(2)-t(1); 
T=t(end);  
df=1/T; 
N=length(st);  
f=-N/2*df:df:N/2*df-df; 
sf=fft(st); 
sf=T/N*fftshift(sf); 
************************************************************
*********************************************************** 
function [out]=sigexpand(d,M)  
N=length(d); 
out=zeros(M,N); 
out(1,:)=d; 
out=reshape(out,1,M*N);
***********************************************************