实验三窗函数法设计FIR数字滤波器..pptVIP

实验三窗函数法设计FIR数字滤波器 实验目的 掌握用窗函数法设计FIR数字滤波器的原理和方法 熟悉线性相位FIR数字滤波器特性 了解各种窗函数对滤波特性的影响 实验原理 如果所希望的滤波器的理想频率响应函数为 ，则其对应的单位脉冲响应为 窗函数设计法的基本原理是用有限长单位脉冲响应 逼近 。由于 往往是无限长序列，且是非因果的，所以用窗函数 将 截断，并进行加权处理，得到： 实验原理 就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数 为 式中，Ｎ为所选窗函数 的长度(N-1阶数)。 实验原理 如果要求线性相位特性，则 还必须满足： 根据上式中的正、负号和长度 的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类。 窗函数基本参数 窗函数 旁瓣峰值幅度/dB 过渡带宽 阻带最小衰减/dB 矩形窗 -13 4π/N? -12 三角形窗 -25 8π/N -25 汉宁窗 -31 8π/N -44 哈明窗 -41 8π/N -53 不莱克曼窗 -57 12π/N -74 凯塞窗( α=7.865) -57 10π/N -80 实验原理 以高通滤波器为例，其常用的设计指标有： 1、通带截止频率 2、阻带截止频率 3、通带波动δp，单位为 dB 4、阻带波动αs，单位为 dB 数字频率 = 模拟频率/采样频率 实验用MATLAB函数介绍 fir1(); fft(); freqz(); abs(); angle(); text(); 实验用MATLAB函数介绍 b=fir1(n,Wn,’ftype’,window)， 窗函数法线性相位滤波器设计，b包含n阶滤波的n+1个系数，n 为阶数、Wn 是截止频率，如果输入是形如[W1 W2]的矢量时，本函数将设计带通/带阻滤波器。 默认选用hamming窗（window缺省）； type=high, 设计高通FIR滤波器 type=stop,设计带阻FIR滤波器 实验用MATLAB函数介绍 [h,w] = freqz(hd,n)：离散系统频响特性 abs(X) ：绝对值（复数的幅值） P = angle(Z) ：相位角 text（x,y,‘s’）：在图面（x,y）位置处书写字符注释。 实验用MATLAB函数介绍 MATLAB中的窗函数 （1）矩形窗（Rectangle Window） ??? 调用格式：w=boxcar(n)，根据长度 n 产生一个矩形窗 w。 （2）三角窗（Triangular Window）??? 调用格式：w=triang(n)，根据长度 n 产生一个三角窗 w。 （3）汉宁窗（Hanning Window）?????? 调用格式：w=hanning(n)，根据长度 n 产生一个汉宁窗 w。 （4）海明窗（Hamming Window）?????? 调用格式：w=hamming(n)，根据长度 n 产生一个海明窗 w。 （5）布拉克曼窗（Blackman Window）? 调用格式：w=blackman(n)，根据长度 n 产生一个布拉克曼窗 w。 （6）恺撒窗（Kaiser Window）??????? 调用格式：w=kaiser(n,beta)，根据长度 n 和影响窗函数旁瓣的β参数产生一个恺撒窗w。 实验用MATLAB函数介绍 例1、设计一个长度为 8 截止频率为0.4π的线性相位 FIR 滤波器。 ? 用矩形窗: Window=boxcar(8); b=fir1(7,0.4,Window); freqz(b,1) ? 实验用MATLAB函数介绍 ? 实验用MATLAB函数介绍 例1、设计一个长度为 8 截止频率为0.4π的线性相位 FIR 滤波器。 ? 用blackman窗: Window=blackman(8); b=fir1(7,0.4,Window); freqz(b,1) 实验用MATLAB函数介绍 实验用MATLAB函数介绍 例2、 设计线性相位带通滤波器，其长度N=15，上下边带截止频率分别为W1= 0.3π，W2=0.5π Window=blackman(15); b=fir1(14,[0.3 0.5],Window); freqz(b,1) 实验用MATLAB函数介绍 实验用MATLAB函数介绍 例3?、用窗函数法设计一FIR带通滤波器： 阻带下截止频率： 通带下截止频率： 通带上截止频率： 阻带上截止频率： 通带最大衰减： 阻带最小衰减： 实验用MATLAB函