1、%用窗函数法设计线性相位低通滤波器clear all;wp=0.5*pi;ws=0.6*pi;wd=ws-wp;%hanning 窗N_hann=ceil(8*pi/wd);wc=(wp+ws)/2;if mod(N_hann,2)=0N_hann=N_hann+1;endhann_win=hanning(N_hann);b_hann=fir1(N_hann-1,wc/pi,hann_win);freq_axis=0:pi/512:pi-pi/512;freqn=0:511/512;H_hann=freqz(b_hann,1,512);figure(1);subplot(211);plot(f
2、reqn,20*log10(abs(H_hann);xlabel(归一化频率 w/pi);ylabel(幅度谱（dB）);title(采用 hanning 窗设计 fir 滤波器幅度谱);subplot(212);plot(freqn,angle(H_hann);xlabel(归一化频率 w/pi);ylabel(相位谱);title(采用 hanning 窗设计 fir 滤波器相位谱);%hamming 窗N_hamm=ceil(8*pi/wd);wc=(wp+ws)/2;if mod(N_hamm,2)=0N_hamm=N_hamm+1;endhamm_win=hamming(N_hamm
3、);b_hamm=fir1(N_hamm-1,wc/pi,hamm_win);H_hamm=freqz(b_hamm,1,512);figure(2)subplot(211);plot(freqn,20*log10(abs(H_hamm);xlabel(归一化频率 w/pi);ylabel(幅度谱（dB）);title(采用 hamming 窗设计 fir 滤波器幅度谱);subplot(212);plot(freqn,angle(H_hamm);xlabel(归一化频率 w/pi);ylabel(相位谱);title(采用 hamming 窗设计 fir 滤波器相位谱);N2_hamm=2*
4、N_hamm;wc=(wp+ws)/2;if mod(N2_hamm,2)=0N2_hamm=N2_hamm+1;endhamm_win=hamming(N2_hamm);b_hamm=fir1(N2_hamm-1,wc/pi,hamm_win);H_hamm=freqz(b_hamm,1,512);figure(3)subplot(211);plot(freqn, 20*log10(abs(H_hamm) );xlabel(归一化频率 w/pi);ylabel(幅度谱（dB）);title(窗长增大一倍，采用 hamming 窗设计 fir 滤波器幅度谱);subplot(212);plot( freqn,angle(H_hamm); xlabel(归一化频率 w/pi);ylabel(相位谱);title(窗长增大一倍，采用 hanning 窗设计 fir 滤波器相位谱);
