[信息与通信]第5章有限长单位脉冲响应FIR滤波器的设计方法.ppt

第5章 有限长单位脉冲响应（FIR）滤波器的设计方法 5.1 线性相位FIR数字滤波器的特性 5.2 窗口设计法（时间窗口法） 5.3 频率取样法 5.1 线性相位FIR数字滤波器的特性  3. h(n)奇对称，N为奇数，h(n)=-h(N-1-n)    四种线性相位FIR滤波器   5.2 窗口设计法（时域）     5.3 频率采样法  三、?约束条件  5.凯塞窗 以上四种窗函数，都是以增加主瓣宽度为代价来降低旁瓣。凯塞窗则可自由选择主瓣宽度和旁瓣衰减。 I0(x)是零阶修正贝塞尔函数，参数β可自由选择，决定主瓣宽度与 旁瓣衰减。β越大，w(n)窗越窄，其频谱的主瓣变宽，旁瓣变小。一般取 4β9。 β=5.44 接近汉明 β=8.5 接近布莱克曼 β=0 为矩形 图2 凯塞窗函数 图1 零阶修正贝塞尔函数 I0(x) x 0 1 而 当MN时， hM(n)≈hd(n) 零阶贝塞尔函数 窗口设计法的主要工作是计算hd(n)和w(n)，当 较为复杂时，hd(n)不容易由反付里叶变换求得。这时一般可用离散傅里叶变换代替连续傅里叶变换，求得近似值： 令 ??过渡带宽At阻带最小衰减为 例2 用凯塞窗设计一FIR低通滤波器，低通边界频率 ，阻带边界频率 ，阻带衰减At不小于50dB。 解： wn=kaiser(30,4.55); nn=[0:29]; alfa=(30-1)/2; hd=sin(0.4*pi*(nn-alfa))./(pi*(nn-alfa)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1))); axis([0,1,-80,10]); grid; xlabel(归一化频率/pi); ylabel(幅度/dB); clear; w1=0.3*pi; w2=0.5*pi; wc=(w1+w2)/2; dw=w2-w1; At=50; if At=50 bata=0.1102*(At-8.7) elseif At21 bata=0.5842*(At-21)^0.4+0.07886*(At-21) else bata=0 end N=ceil((At-8)/(2.286*dw)) N=ceil((At-8)/(2.286*dw)) alfa=(N-1)/2; M=1024; for mm=0:M-1 w=2*pi/M*mm; if w=wc Hd(mm+1)=exp(-j*alfa*w); else Hd(mm+1)=0; end end hd=real(ifft(Hd)); wn=kaiser(N,bata); nn=[0:N-1]; hd=real(hd(1:N)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1))); axis([0,1,-80,10]); grid; xlabel(归一化频率/pi); ylabel(幅度/dB); 工程上，常给定频域上的技术指标，所以采用频域设计更直接。 一、基本思想 使所设计的FIR数字滤波器的频率特性在某些离散频率点上的值准确地等于所需滤波器在这些频率点处的值，在其它频率处的特性则有较好的逼近。 内插公式 第二种FIR滤波器（h(n)偶对称，N为偶数）， ，即 是 的零点，为单根； 第四种FIR滤波器（h(n)奇对称，N为偶数）， ，即 是 的零点，为单根； 第三种FIR滤波器（h(n)奇对称，N为奇数）， ，即 是 的零点，为单根。 如果希望得到的滤波器的理想频率响应为 ，那么 FIR滤波器的设计就在于寻找一个传递函数 去逼近 ，逼近方法有三种： 窗口设计法（时域逼近） 频率采样法（频域逼近） 最优化设计（等