www.pudn.com > DSP3_MATLAB.rar > exa030702.m

%---------------------------------------------------------------------------- 
% exa030703, for example 3.7.3 
% to discuss the effect of padding with zeros for DTFT;  
%---------------------------------------------------------------------------- 
clear; 
 
% 计算长度为N的原始信号的DTFT 
 f1=2.67;f2=3.75;f3=6.75;fs=20;w=2*pi/fs; 
 N=16; 
 x=sin(w*f1*(0:N-1))+sin(w*f2*(0:N-1)+pi/2)+sin(w*f3*(0:N-1)); 
 f=0:fs/N:fs/2-1/N; 
 X=fft(x);  
 X=abs(X); 
 f=fs/N*(0:N/2-1); 
 subplot(221) 
 stem(f,X(1:N/2),'.');grid on; 
 xlabel('Hz') 
  
% 在数据末补N个零  
 x(N:2*N-1)=0; 
 X=fft(x); X=abs(X); 
 f=fs*(0:N-1)/(2*N); 
 subplot(222) 
 stem(f,X(1:N),'.');grid on; 
 xlabel('Hz') 
  
% 在数据末补7*N个零 
 x(N:8*N-1)=0; 
 X=fft(x); X=abs(X); 
 f=fs*(0:4*N-1)/(8*N); 
 subplot(223) 
 stem(f,X(1:4*N),'.');grid on; 
 xlabel('Hz') 
  
% 在数据末补29*N个零 
 x(N:30*N-1)=0; 
 X=fft(x); X=abs(X); 
 f=fs*(0:15*N-1)/(30*N); 
 subplot(224) 
 plot(f,X(1:15*N));grid on; 
 xlabel('Hz')