q719.m

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

%《MATLAB在电子信息课程中的应用》第七章例7.19程序q719
% 用DFT作谱分析
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clear;close all
fs=400; T=1/fs;	%采样频率为400Hz
Tp=0.04; N=Tp*fs;	%采样点数N
N1=[N, 4*N, 8*N];			% 设定三种截取长度供调用
st=['|X1(jf)|';'|X4(jf)|';'|X8(jf)|'];	% 设定三种标注语句供调用
%矩形窗截断
for m=1:3
   n=1:N1(m);
   xn=cos(200*pi*n*T)+sin(100*pi*n*T)+cos(50*pi*n*T);%产生采样序列x(n)
   Xk=fft(xn,4096);  %4096点DFT，用FFT实现
   fk=[0:4095]/4096/T;
   subplot(3,2,2*m-1)
   plot(fk,abs(Xk)/max(abs(Xk)));ylabel(st(m,:))
   if m==1 title('矩形窗截取');end
end
%加hamming窗改善谱间干扰
for m=1:3
   n=1:N1(m);
   wn=hamming(N1(m)); %调用工具箱函数hamming产生N长hamming窗序列wn
   xn=(cos(200*pi*n*T)+sin(100*pi*n*T)+cos(50*pi*n*T)).*wn';
   Xk=fft(xn,4096);  %4096点DFT，用FFT实现
   fk=[0:4095]/4096/T;
   subplot(3,2,2*m)
   plot(fk,abs(Xk)/max(abs(Xk)));ylabel(st(m,:))
   if m==1 title('Hamming窗截取');end
end