线性相位

matlab用窗函数法设计线性相位fir低通滤波器

%用窗函数法设计线性相位低通滤波器clearall;wp=0.5某pi;w=0.6

某pi;wd=w-wp;

%hanning窗

N_hann=ceil(8某pi/wd);wc=(wp+w)/2;

ifmod(N_hann,2)==0N_hann=N_hann+1;end

hann_win=hanning(N_hann);

b_hann=fir1(N_hann-1,wc/pi,hann_win);freq_a某

i=[0:pi/512:pi-pi/512];freqn=[0:511]/512;

H_hann=freqz(b_hann,1,512);figure(1);ubplot(211);

plot(freqn,20某log10(ab(H_hann)));

某label('归一化频率w/pi');ylabel('幅度谱（dB）');title('采

用hanning窗设计fir滤波器幅度谱');ubplot(212);

plot(freqn,angle(H_hann));

某label('归一化频率w/pi');ylabel('相位谱');title('采用

hanning窗设计fir滤波器相位谱');%hamming窗

N_hamm=ceil(8某pi/wd);wc=(wp+w)/2;

ifmod(N_hamm,2)==0

N_hamm=N_hamm+1;end

hamm_win=hamming(N_hamm);

b_hamm=fir1(N_hamm-

1,wc/pi,hamm_win);H_hamm=freqz(b_hamm,1,512);figure(2)ubplot(211

);

plot(freqn,20某log10(ab(H_hamm)));

某label('归一化频率w/pi');ylabel('幅度谱（dB）');title('采

用hamming窗设计fir滤波器幅度谱');ubplot(212);

plot(freqn,angle(H_hamm));

某label('归一化频率w/pi');ylabel('相位谱');title('采用

hamming窗设计fir滤波器相位谱');N2_hamm=2某N_hamm;

wc=(wp+w)/2;

ifmod(N2_hamm,2)==0

N2_hamm=N2_hamm+1;end

hamm_win=hamming(N2_hamm);

b_hamm=fir1(N2_hamm-

1,wc/pi,hamm_win);H_hamm=freqz(b_hamm,1,512);figure(3)ubplot(211

);

plot(freqn,20某log10(ab(H_hamm)));

某label('归一化频率w/pi');ylabel('幅度谱（dB）');

title('窗长增大一倍，采用hamming窗设计fir滤波器幅度谱

');ubplot(212);

plot(freqn,angle(H_hamm));

某label('归一化频率w/pi');ylabel('相位谱');

title('窗长增大一倍，采用hanning窗设计fir滤波器相位谱');