实验报告(五)
——FIR 滤波器的设计
一、实验任务:
一、利用窗函数法设计一个线性相位低通 FIR 数字滤波器,wp=0.4pi,wst=0.6pi,阻
带最小衰减-50dB。
二、利用频率抽样法,设计一个低通 FIR 数字滤波器,wc=0.5pi,抽样点数为 N=33,
要求滤波器具有线性相位。
(1)增加一点过渡带进行优化,过渡带抽样点值为 0.5;
(2)增加抽样点数至 N=65,加两点过渡带优化,过渡带抽样点值为 0.5886,0.1065;
三、实验源代码:
一、
wp=0.4*pi;wst=0.6*pi;
wc=(wp+wst)/2;
N=ceil(6.6*pi/(wst-wp)); --------------------------
b=fir1(N,wc/pi,hamming(N+1));----------------
[h,w]=freqz(b,1,500);----------------------------------
plot(w/pi,20*log10(abs(h)),'-');xlabel('H');
ylabel('20log|H|');grid on;axis([0 1 -100 10]);
%二、
N=33;
Hw=[ones(1,9),zeros(1,16),ones(1,8)];-----------
k=[0:N-1];
angH=-2*pi/N*k*(N-1)/2;----------------------------
H=Hw.*exp(j*angH);---------------------------------
b=ifft(H,N);--------------------------------------------
[H1,w1]=freqz(b,1,500);
(1)增加一点过渡带进行优化,
过渡带抽样点值为 0.5;
N=33;
Hw=[ones(1,9),zeros(1,16),ones(1,8)];-------------
Hw(10)=0.5;Hw(25)=0.5;
k=[0:N-1];
angH=-2*pi/N*k*(N-1)/2;---------------------------
H=Hw.*exp(j*angH);-----------------------------------
b=ifft(H,N);-------------------------------------------
[H2,w2]=freqz(b,1,500); ----------------------