名师推荐第7章IIR数字滤波器设计.ppt

1、第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 第第7章章 IIR数字滤波器设计数字滤波器设计 7.1 IIR滤波器原理滤波器原理 7.2 使用使用DSP Builder设计设计IIR滤波器滤波器 7.3 在在Quartus II中使用中使用IIR滤波器滤波器IP核核 第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 7.1 IIR滤波器原理滤波器原理 IIR滤波器一般采用递归方式来实现。也就是说，滤波器当前输出是输入序列和以前各输出值、的函数，这可以用下列差分方程来表示：NiiMiiinyainxbny10)()()(7-1)第第7 7章章 IIRIIR数字滤波器设计数字滤

2、波器设计 其中MN，N是IIR滤波器的阶数。相应地，IIR滤波器的系统函数可以表示为 NiiiMiiizazbzH101)(7-2)第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 在上一章已经提及到，FIR滤波器的系统函数只有零点。而IIR滤波器除了具有极点以外，一般还存在零点。由于极点的存在，IIR滤波器用递归结构来实现较为简单。实现IIR滤波器的基本结构共有三中：直接型、级联型和并联型。下面简单介绍前两种IIR滤波器的结构。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 1.直接型直接型 利用公式(7-2)，可以直接导出直接I型的IIR滤波器结构，可用下式来表示：Nii

3、iMiiizazbzH1011)(7-3)第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 可以看出，和反馈环节 中的延时单元不能共用，需要个延 时单元。为了减少延时单元，对于线性系统，公式(7-3)也可以写成：Miiizb0Niiiza111MiiiNiiizbzazH0111)(第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 由此可以导出直接型的滤波器结构。图7-1表示了一个直接型IIR滤波器的结构，该滤波器的阶数是2阶，共存在两个反馈环节。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 z1z1x(n)a1a2b0b1b2y(n)图7-1 直接型IIR滤波器结

4、构第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 2.级联型级联型 对式(7-2)的分子、分母进行因式分解，由于中的系数都为实数，的极、零点只可能是实数或者复共轭对，对于复共轭对因子，可以复合成二阶因子：212112312111231211)1()1()1()1()(NiiiNiiMiiiMiizazaazbzbbAzH第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 式中系数都为实数。如果把实数因子(一阶因子)看成是二次项系数为0的二阶因子，则上式可以写成：101022112211)(11)(ssNiiNiiiiizHAzazazbzbAzH2211221111)(zaza

5、zbzbzHiiiii其中 第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 由此可见，是一个2阶的IIR滤波器，N阶IIR滤波器可以看成是由多个2阶IIR滤波器级联而成的，即IIR滤波器可以采用下列结构来表示，也就是级联型结构(见图7-2)：第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 z11a11b10b12x(n)G1a21b20a12z1b11y1(n)a22z1z1y(n)b21b22图7-2 级联型IIR滤波器结构第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 7.2 使用DSP Builder设计IIR滤波器 7.2.1 4阶直接阶直接型型IIR滤波

6、器设计滤波器设计 在此需要设计一个4阶IIR滤波器。给定滤波器的系统函数如下：4321432112.072.074.195.1105.020.031.020.005.0)(zzzzzzzzzH第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 即滤波器系数为 a0=1b0=0.05 a1=-1.95b1=0.20 a2=1.74 b2=0.31 a3=-0.72b3=0.20 a4=0.12 b4=0.05 可以看出，这是一个4阶IIR滤波器，我们采用直接型结构来实现该滤波器，其具体结构可由信号流图(见图7-3)来表示。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 1a1b0

7、 x(n)a2a3a3z1z1z1z1b1b2b3b4y(n)图图7-3 直接直接型型4阶阶IIR滤波器信号流图滤波器信号流图第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 在DSP Builder中，建立一个新模型，按照图7-3调用DSP Builder模块进行设计。设计完的4阶IIR滤波器模型如图7-4所示。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-4 4阶阶IIR滤波器模型滤波器模型第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 该模型中各个模块的参数设置如下：X模块：(Altbus)库：Altera DSP Builder中Bus Manip

8、ulation库参数“Bus Type”设为“signed Fractional”(有符号小数)参数“Node Type”设为“Input port”参数“number of bits.”设为“2”参数“.number of bits”设为“8”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 Y模块：(Altbus)库：Altera DSP Builder中Bus Manipulation库参数“Bus Type”设为“signed Fractional”参数“Node Type”设为“Output port”参数“number of bits.”设为“4”参数“.number of

9、 bits”设为“23”BusConversion模块：(BusConversion)库：Altera DSP Builder中Bus Manipulation库第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 参数“Input Bus Type”设为“signed Fractional”参数“Input number of bits.”设为“4”参数“Input.number of bits”设为“23”参数“Output Bus Type”设为“signed Fractional”参数“Output number of bits.”设为“2”参数“Output.number of

10、bits”设为“8”FeedBack Adder模块：(Parallel Adder Subtractor)库：Altera DSP Builder中Arithmetic库参数“Number of Inputs”设为“2”参数“Add(+)Sub(-)”设为“+-”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 FeedBack Adder1模块：(Parallel Adder Subtractor)库：Altera DSP Builder中Arithmetic库参数“Number of Inputs”设为“4”参数“Add(+)Sub(-)”设为“+”FeedForward Add

11、er模块：(Parallel Adder Subtractor)库：Altera DSP Builder中Arithmetic库参数“Number of Inputs”设为“5”参数“Add(+)Sub(-)”设为“+”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 Delay、Delay1、Delay2、Delay3模块：(Delay)库：Altera DSP Builder中Storage库参数“Depth”设为“1”参数“Clock Phase Selection”设为“1”A1、A2、A3、A4、B0、B1、B2、B3、B4模块：(Gain)库：Altera DSP Buil

12、der中Arithemtic库参数“Gain Value”按照设计要求中指定的系数设置(直接输入)参数“Map Gain Value to Bus Type”设为“Signed Fraction”参数“Gain value number of bits.”设为“2”参数“.Gain value number of bits”设为“12”参数“Number of Pipeline Levels”设为“0”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 另外，用于Simulink仿真的模块的参数设置如下：Pulse Generator模块：(Pulse Generator)库：Simul

13、ink中Sources库参数“Pulse type”设为“Sample based”参数“Amplitude”设为“1”参数“Period(number of samples)”设为“1000”参数“Pulse width(number of samples”设为“1”参数“Phase delay(number of samples)”设为“10”参数“Sample time”设为“1”“Interpret vector parameters as 1-D”设为“”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 Scope模块：(Scope)库：Simulink中sinks库 参数“

14、Number of Axes”为“2”在这个实现中，使用了增益(Gain)模块来完成常数系数的乘法。这里使用了有符号小数(Signed Fractional)来简化IIR滤波器系数的代入过程。值得注意的是，参数“number of bits.”设置的是整数部分的位宽，其中最高位为符号位；参数“.number of bits”设置了小数部分的位宽。假如把参数“number of bits.”设置为“2”，那就意味着该小数的整数部分只能是0、1、-1、-2(分别对应00、01、10或10、10)，而且整数部分要为-2，必须是小数部分为全零。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计

15、建立IIR滤波器模型后就可以进行仿真了。Simulink仿真模块“Pulse Generator”(脉冲发生器)模拟了一个单位冲激函数。通过“Scope”模块来观察该直接型IIR滤波器的冲激响应。修改Simulink的仿真参数设置，设置“Stop time”仿真停止时间为“50”，“Solver options”中“Type”为“Fixed-Step”。启动仿真。仿真结果见图7-5。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-5 IIR滤波器仿真结果滤波器仿真结果第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 7.2.2 4阶级联型阶级联型IIR滤波器设计滤波器

16、设计 1.建立模型建立模型 参照图7-2，建立一个4阶的级联型IIR滤波器模型，该模型共由两节2阶直接型IIR滤波器构成，见图7-6。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-6 4阶级联型阶级联型IIR滤波器滤波器第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 模型中各个模块的参数设置如下：X模块：(Altbus)库：Altera DSP Builder中Bus Manipulation库参数“Bus Type”设为“signed Fractional”(有符号小数)参数“Node Type”设为“Input port”参数“number of bits.”

17、设为“2”参数“.number of bits”设为“8”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 Y模块：(Altbus)库：Altera DSP Builder中Bus Manipulation库参数“Bus Type”设为“signed Fractional”参数“Node Type”设为“Output port”参数“number of bits.”设为“4”参数“.number of bits”设为“23”BusConv模块：(BusConversion)库：Altera DSP Builder中Bus Manipulation库参数“Input Bus Type”设

18、为“signed Fractional”参数“Input number of bits.”设为“4”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 参数“Input.number of bits”设为“18”参数“Output Bus Type”设为“signed Fractional”参数“Output number of bits.”设为“2”参数“Output.number of bits”设为“15”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 BusConv1、BusConv4模块：(BusConversion)库：Altera DSP Builder中Bus

19、Manipulation库参数“Input Bus Type”设为“signed Fractional”参数“Input number of bits.”设为“3”参数“Input.number of bits”设为“15”参数“Output Bus Type”设为“signed Fractional”参数“Output number of bits.”设为“2”参数“Output.number of bits”设为“15”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 BusConv2、BusConv3、BusConv5模块：(BusConversion)库：Altera DSP

20、Builder中Bus Manipulation库参数“Input Bus Type”设为“signed Fractional”参数“Input number of bits.”设为“4”参数“Input.number of bits”设为“25”参数“Output Bus Type”设为“signed Fractional”参数“Output number of bits.”设为“2”参数“Output.number of bits”设为“15”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 FeedBackAdder、FeedBackAdder2模块：(Parallel Adde

21、r Subtractor)库：Altera DSP Builder中Arithmetic库参数“Number of Inputs”设为“2”参数“Add(+)Sub(-)”设为“+-”FeedBackAdder1、FeedBackAdder3模块：(Parallel Adder Subtractor)库：Altera DSP Builder中Arithmetic库参数“Number of Inputs”设为“2”参数“Add(+)Sub(-)”设为“+”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 FeedForwardAdder、FeedForwardAdder1模块：(Para

22、llel Adder Subtractor)库：Altera DSP Builder中Arithmetic库参数“Number of Inputs”设为“3”参数“Add(+)Sub(-)”设为“+”Delay、Delay1、Delay2、Delay3模块：(Delay)库：Altera DSP Builder中Storage库参数“Depth”设为“1”参数“Clock Phase Selection”设为“1”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 G、A11、A12、A21、A21、B10、B11、B12、B20、B21、B22模块：(Gain)库：Altera DSP

23、 Builder中Arithemtic库参数“Map Gain Value to Bus Type”设为“Signed Fraction”参数“Gain value number of bits.”设为“2”参数“.Gain value number of bits”设为“10”参数“Number of Pipeline Levels”设为“0”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 2.IIR滤波器系数的计算滤波器系数的计算假设4阶IIR滤波器的设计指标如下：采样频率Fs为50MHz；滤波器类型为高通(HighPass)；滤波器Fc为5MHz；滤波器设计类型为Butterw

24、orth。打开MATLAB中的FDATool进行滤波器系数计算，见图7-7。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-7 计算计算IIR滤波器系数滤波器系数第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 设置阶数为4，输入相关参数指标，点击“Design Filter”按钮完成IIR滤波器设计，接着打开相关分析。图7-8显示的是IIR滤波器的相频特性，图7-9显示了IIR滤波器的阶跃响应。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-8 相频特性相频特性第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-9 阶跃响应阶跃响应第

25、第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 3.设置系数设置系数 选择FDATool的“File”菜单中的“Export”，导出IIR滤波器系数。对于级联型IIR滤波器，导出的系数分成两个部分：SOS矩阵和G增益，其中对于多个级联的IIR滤波器，SOS矩阵为一个二维矩阵。在此先把SOS矩阵和G导出到Matlab的工作区(Workspace)，设置上面级联型IIR滤波器中的各个增益模块的“Gain Value(增益值)”。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 G模块的“Gain Value”是G；B10模块为SOS(1,1)，B11模块为SOS(1,2)，B12模

26、块为SOS(1,3)；A11模块为SOS(1,5)，A12模块为SOS(1,6)；B20模块为SOS(2,1)，B21模块为SOS(2,2)，B22模块为SOS(2,3)；A21模块为SOS(2,5)，A22模块为SOS(2,6)。设置完增益值后更新一下模型，IIR滤波器模型的系数就全部设置完成了，见图7-6。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 4.Simulink仿真仿真 设置Simulink仿真参数“Stop time(停止时间)”为1e-6，类型为“Fixed-step”。模型中的Simulink仿真模块的参数设置如下：第第7 7章章 IIRIIR数字滤波器设计数字

27、滤波器设计 Step模块：(Step)库：Simulink中Sources库参数“Step time”设为“1e-7”参数“Initial value”设为“0”参数“Final value”设为“1”参数“Sample time”设为“2e-8”“Interpret vector parameters as 1-D”设为“”第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 Scope模块：(Scope)库：Simulink中sinks库参数“Number of Axes”为“2”启动仿真。由于Step模块生成了一个阶跃函数，故仿真结果为IIR滤波器的阶跃响应，见图7-10。5.由由S

28、imulink模型转成模型转成VHDL 打开SingalCompiler，选定对应器件，把模型转成VHDL文件，并在SignalCompiler中选择QuartusII进行综合。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图7-10 Simulink仿真结果一阶跃响应 第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 6.适配下载适配下载 在QuartusII中打开SignalCompiler建立的Quartus项目文件，选择具体器件，锁定管脚，完成适配后下载至FPGA DSP开发板中。在FPGA DSP开发板上加入高频信号源，验证IIR滤波器的高通滤波效果。第第7 7章

29、章 IIRIIR数字滤波器设计数字滤波器设计 7.3 在QuartusII中使用IIR滤波器IP核 从上面的流程来看，虽然借助DSP Builder中丰富的DSP功能模块和MATLAB强大的计算环境，不再需要关心DSP模块的具体实现过程，也就不会陷入复杂的滤波器设计算法，从而大大简化了IIR滤波器的设计过程。但是，以上述设计方法得到的IIR滤波器，在资源利用上和性能上往往不是最优的。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 同FIR Filter Core一样，Altera以及其它FPGA厂商也提供了IIR Filter Core。IIR Filter Core同FIR Fil

30、ter Core一样，与DSP Builer相结合，可以在Simulink环境中使用(图7-11)。IIR滤波器核在DSP Builder中的使用，与FIR滤波器核大致是相同的，因此本节就不再一一详述了。IIR滤波器核除了与DSP Builder可以衔接开发之外，也可以在QuartusII中直接使用。由于在前文中没有详细提及在QuartusII中使用IP Core，因而在这里加以介绍。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-11 IIR Core模块模块第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 7.3.1 配置配置QuartusII以便使用以便使用II

31、R滤波器核滤波器核 首先，必须确认IIR滤波器核(IIR Compiler)和QuartusII是否安装。打开QuartusII集成环境，新建一个项目。选择菜单“Assignments”“Settings”，打开QuartusII的设置对话框，在对话框的左侧选择“User Libraries(用户库)”，见图7-12。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-12 设置用户库设置用户库第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 7.3.2 使用使用IIR滤波器核滤波器核 首先完成IIR滤波器核的安装配置，然后按照第2章给出的相类似的方法，定制使用IIR滤波

32、器核。具体步骤如下：1.打开打开MegaWizard Plug-In Manager 选择菜单“Tools”“MegaWizard Plug-In Manager”，接着打开MegaCore向导。在窗口的左侧选中“IIR Compiler”，在右侧选择实现时的HDL类型，并输入生成IIR滤波器核的名字(见图7-13)。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-13 选择选择IIR Core第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 2.选择选择IIR滤波器结构滤波器结构 在前文中已经提到，常用的IIR滤波器结构共有三种：直接型(Direct Form)，级

33、联型(Cascaded)和并联型(Parallel)。Altera的IIR滤波器核对这三种IIR结构都是支持的。图7-14是IIR滤波器结构设置对话框，在这里我们选择了级联型进行实现，滤波器阶数设为8。图7-15、图7-16分别表示了另外两种滤波器结构的计算公式和结构。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图7-14 IIR滤波器结构第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-15 直接直接型型IIR结构结构第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图7-16 并联型IIR结构 第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计

34、3.导入滤波器系数导入滤波器系数 由于IIR滤波器系数的计算要比FIR滤波器复杂得多，对于Altera的IIR滤波器核，其系数计算需要通过其它工具进行，比如可以使用MATLAB的FDATool来完成，再通过图7-17的对话框导入(通过文件)。系数导入后，IIR Compiler可以显示当前系数的IIR滤波器的幅频响应。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-17 确定确定IIR系数系数第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 对于IIR滤波器而言，还可以通过导入零极点值来确定IIR滤波器的系数。IIR Compiler也提供了导入零极点来确定IIR系数

35、的方式。图7-18显示了导入后IIR滤波器的零极点图。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-18 IIR滤波器的零极点滤波器的零极点第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 4.IIR滤波器系数的量化滤波器系数的量化 同FIR滤波器一样，在FPGA上实现，IIR系数必须进行量化处理，见图7-19。我们在这里设置了输入信号的位宽为8位，系数位宽为16位。采用自动方式让IIR Compiler具体决定系数的归一化问题。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-19 量化量化IIR滤波器系数滤波器系数第第7 7章章 IIRIIR数字

36、滤波器设计数字滤波器设计 5.IIR滤波器的具体实现滤波器的具体实现 为了提高IIR滤波器的性能，同时减少IIR滤波器使用的资源，需要对IIR滤波器实现时的某些细节进行调整，如图7-20所示。然后点击“Finish”完成IIR滤波器核的设置。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图7-20 IIR滤波器实现的调整 第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 6.调用设计好的调用设计好的IIR滤波器核滤波器核 设计好的IIR滤波器核可以通过多种方式调用。图7-21显示了在QuartusII的原理图方式设计中，对设计好的IIR滤波器核作为一个元件(Symbol)进行调用。当然，采用VHDL调用设计好的IIR滤波器核也是相当方便的。IIR Compiler在设计时提供了调用模板文件供设计者参考。另外，IIR Compiler也能生成一些仿真测试文件来简化对核的测试过程。第第7 7章章 IIRIIR数字滤波器设计数字滤波器设计 图图7-21 调用设计好的调用设计好的IIR滤波器核滤波器核