fg_05_08.m

英文书《Digital Signal Processing with Examples in MATLAB》附带的MATLAB实例

% fg_05_08  Kaiser power spectra.
N=127; dB0=-120;
beta=[0 sqrt(3)*pi 0 2*sqrt(2)*pi];
name=str2mat('Hamming','Kaiser; \beta=5.44','Blackman',...
             'Kaiser; \beta=8.89');
nu=linspace(0,.5,500);
sp_fig(1,9,6);
type=[5 7 6 7];
clr=['r' 'r' 'c' 'c'];
for i=1:4
   w=window(N,type(i),beta(i));
   H=gain_f(w,1,nu);
   dB=20*log10(max(abs(H/H(1)),-log10(-dB0/2)));
   subplot(2,2,i);
   axis([0 .5 dB0 0]);
   line(nu,dB,'color',clr(i));
   grid;
   set(gca,'fontsize',16);
   if(i==3|i==4)
      xlabel('\nu (Hz-s)');
   end
   if(i==1|i==3)
      ylabel('dB');
   end
   ht(1)=text(.21,-10,name(i,:));
   set(ht,'fontsize',16);
end