線性相位

matlab用窗函數(shù)法設(shè)計線性相位fir低通濾波器

%用窗函數(shù)法設(shè)計線性相位低通濾波器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窗設(shè)計fir濾波器幅度譜');ubplot(212);

plot(freqn,angle(H_hann));

某label('歸一化頻率w/pi');ylabel('相位譜');title('采用

hanning窗設(shè)計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窗設(shè)計fir濾波器幅度譜');ubplot(212);

plot(freqn,angle(H_hamm));

某label('歸一化頻率w/pi');ylabel('相位譜');title('采用

hamming窗設(shè)計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窗設(shè)計fir濾波器幅度譜

');ubplot(212);

plot(freqn,angle(H_hamm));

某label('歸一化頻率w/pi');ylabel('相位譜');

title('窗長增大一倍，采用hanning窗設(shè)計fir濾波器相位譜');