%《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