samp7_2.m

matlab采用512个频率点绘制各种窗函数的幅频特性

%Samp7_2
clf;Nf=512;                  %窗函数复数频率特性的数据点数
Nwin=20;                      %窗函数数据长度
figure(1)
for ii=1:4
   switch ii
   case 1
      w=boxcar(Nwin);              %矩形窗
      stext='矩形窗';
   case 2
      w=hanning(Nwin);              %汉宁窗
      stext='汉宁窗';
   case 3
      w=hamming (Nwin);              %哈明窗
      stext='哈明窗';
   case 4
      w=bartlett(Nwin);                % Bartlett窗
      stext='Bartelett窗';
   end
   [y,f]=freqz(w,1,Nf); %求解窗函数的幅频特性，窗函数相当于一个数字滤波器
   mag=abs(y);%求得窗函数幅频特性
   posplot=['2,2,' int2str(ii)];
   subplot(posplot);
   plot(f/pi,20* log10(mag))%/max(mag)));            %绘制窗函数的幅频特性
   xlabel('归一化频率');ylabel('振幅/dB');
   title(stext);grid on;
end
figure(2)
for ii=1:4
   switch ii
   case 1
      w=blackman(Nwin);                 %Blackman 窗
      stext='Blackman窗';
   case 2
      w=triang(Nwin);                  %三角窗
      stext='三角窗';
   case 3
      w=kaiser(Nwin,4);               %Kaiser窗
      stext='Kaiser窗(Beta=4)';
   case 4
      w=chebwin(Nwin,40);            %Chebyshev 窗
      stext='Chebyshev窗(r=40)';
   end
   [y,f]=freqz(w,1,Nf);               %以 Nf点数求解窗函数的幅频响应
   mag=abs(y);                        %求得窗函数幅频响应
   posplot=['2,2,' int2str(ii)]; 
   subplot(posplot);plot(f/pi,20* log10(mag/max(mag))); %绘制幅频响应
   xlabel('归一化频率');ylabel('振幅/dB');
   title(stext);grid on;
end