北京邮电大学《数字信号处理》门爱东dsp04iir 数字滤波器设计和实现

IIR数字滤波器设计和实现北京邮电大学《数字信号处理》第四章｜主讲：门爱东Contents本章内容概览数字信号处理—IIR数字滤波器设计和实现01引言：数字滤波器设计基础02模拟滤波特性的逼近03脉冲响应不变法与数字仿真04双线性变换法05数字滤波器变换与实现结构06优化设计与MATLAB实践CHAPTER01引言：数字滤波器设计基础理解数字滤波器的设计流程、系统函数表达与方法选择IIRFILTERDESIGN数字滤波器设计的三个基本步骤数字滤波器设计是一个从需求定义到数学建模再到工程实现的系统过程，需依次完成性能指标确定、系统函数逼近和有限精度运算三个关键阶段。STEP01确定性能要求根据实际应用需求明确滤波器类型（低通/高通/带通/带阻）及各项频率响应指标关键参数包括通带截止频率ωp、阻带起始频率ωs、通带最大衰减Rp和阻带最小衰减Rs这些指标直接决定滤波器的选择范围和设计复杂度ωp·ωs·Rp·RsSTEP02系统函数逼近用因果稳定的系统函数H(z)逼近目标性能，可选IIR（有理函数）或FIR（多项式）两类IIR系统函数为z的有理分式，可用较低阶数实现陡峭过渡带，但需保证极点位于单位圆内常用逼近方法包括巴特沃斯、切比雪夫和椭圆函数等H(z)STEP03有限精度实现选择运算结构（直接型、级联型、并联型等），不同结构对量化误差敏感度差异显著确定运算字长和有效数字处理方法，控制有限字长效应对滤波器频率特性的影响级联型结构通常具有更好的数值稳定性直接·级联·并联Chapter04·IIRDigitalFiltersIIR系统函数的数学表达与设计方法IIR数字滤波器的系统函数是z的有理函数，分子分母多项式的系数决定了滤波器的零极点分布和频率响应特性。设计方法包括零极点累试法、模拟滤波器转换法和优化技术三种路径，其中模拟转换法因理论成熟、设计流程规范而最为常用。01H(z)有理分式IIR系统函数H(z)表示为z的有理分式，分子阶数M和分母阶数N决定滤波器复杂度和频率选择性M/N阶数02零极点累试法通过在z平面上配置零极点位置逐步逼近目标频率响应，适用于简单滤波器设计Z平面配置03模拟滤波器转换法先设计Ha(s)再映射为H(z)，利用成熟的模拟滤波器理论实现高效设计Ha(s)→H(z)04优化技术通过计算机辅助迭代搜索系统函数的最优系数，适用于复杂指标下的精确设计迭代搜索DIGITALFILTERCOMPARISONIIR与FIR数字滤波器的核心差异IIR和FIR是数字滤波器的两大体系，在系统结构、相位特性、稳定性和计算效率方面存在根本性差异。IIR滤波器特征STRUCTURE系统函数为有理分式，存在反馈路径，冲激响应理论上无限延续，故称无限冲激响应EFFICIENCY可用较低阶数（如4-8阶）实现陡峭的幅频特性过渡带，计算量和存储需求相对较小4–8阶STABILITY极点必须位于z平面单位圆内才能保证稳定，设计时需额外验证稳定性条件FIR滤波器特征STRUCTURE系统函数仅有零点（无分母多项式），冲激响应在有限时间内衰减为零，天然绝对稳定PHASE可通过对称系数设计实现严格的线性相位，避免信号通过滤波器后产生相位失真LINEARPHASECOST实现同等选择性通常需要比IIR高5-10倍的阶数，在实时处理场景中计算开销较大5–10×阶数CHAPTER02模拟滤波特性的逼近掌握巴特沃思与切比雪夫经典滤波器的幅频特性、设计方法与MATLAB实现FILTERCLASSIFICATION滤波器分类与频带基本划分滤波器按处理信号类型、有源/无源属性和通频带特性可分为多个类别。按通频带划分为低通、高通、带通和带阻四种基本类型，每种滤波器的幅频特性均包含通带、阻带和过渡带三个关键区域，过渡带宽度是衡量滤波器频率选择性的核心指标。按信号类型与应用分类01模拟滤波器输入输出均为模拟量，数字滤波器处理离散数字信号，两者设计理论相互关联02有源滤波器包含运算放大器等有源器件可提供增益，无源滤波器仅由电阻电容电感组成模拟/数字按通频带特性分类01低通滤波器允许低频信号通过并抑制高频，高通滤波器反之，两者互为频率变换的对偶关系02带通滤波器仅允许特定频带通过，带阻滤波器抑制特定频带，可通过低通原型经频率变换获得四种基本类型幅频特性的三个区域01通带内信号衰减小于容许值Rp，阻带内衰减大于要求值Rs，过渡带连接通带和阻带02理想滤波器过渡带宽度为零（矩形响应），实际滤波器的设计目标是在给定阶数下最小化过渡带过渡带宽度IIRFilterDesign巴特沃思(Butterworth)滤波器的幅频特性巴特沃思滤波器以通带内最大平坦的振幅特性著称，阶数N决定过渡带陡峭程度，N越大越接近理想矩形响应。振幅平方函数|Ha(jΩ)|²=1/(1+(Ω/Ωc)2N)Ωc为3dB截止频率，N为滤波器阶数通带最大平坦所有导数在Ω=0处为零，这是Butterworth区别于其他逼近方式的核心特征单调下降特性正频率范围内幅频特性单调下降，N越大过渡带越陡峭，阻带衰减速率为每倍频程6NdB极点分布2N个极点均匀分布在s平面半径为Ωc的圆周上，左半平面N个极点构成稳定传递函数不同阶数Butterworth滤波器幅频响应N越大过渡带越陡峭，越接近理想矩形响应IIRFilterDesign巴特沃思滤波器设计方法与步骤巴特沃思滤波器的设计流程包括指标转换、阶数计算、查表获取归一化传递函数和频率去归一化四个步骤。阶数N由通带和阻带的衰减要求共同决定，计算结果向上取整确保实际性能不低于设计指标，是工程设计中保证可靠性的关键原则。STEP01根据通带截止频率Ωp、阻带起始频率Ωs及衰减指标Rp、Rs，计算频率比λ=Ωs/ΩpSTEP02利用公式计算最小阶数N并向上取整，确保滤波器阶数满足通带与阻带的衰减要求N≥lg[(10Rs/10−1)/(10Rp/10−1)]/(2lgλ)STEP03根据阶数N查归一化巴特沃思多项式表，获得归一化传递函数Ha(s')的极点与系数STEP04通过频率去归一化s=s'·Ωc将截止频率调整至实际值，得到最终系统函数Ha(s)MATLABFunctions巴特沃思滤波器的MATLAB实现MATLAB通过buttord和butter两个核心函数实现Butterworth滤波器的自动化设计，buttord计算最小阶数和截止频率，butter生成系统函数参数。阶数计算·buttord[n,Wn]=buttord(Wp,Ws,Rp,Rs,'s')计算给定指标下最小Butterworth模拟滤波器阶数n和3dB截止频率Wn。n,Wn零极点增益形式[z,p,k]=butter(n,Wn,'s')返回归一化滤波器的零点z、极点p和增益k，可直接构建系统函数。z,p,k传递函数系数[b,a]=butter(n,Wn,'s')返回传递函数的分子系数b和分母系数a，便于频响分析和仿真验证。b,a多类型滤波器支持带通/带阻设计中Wp和Ws为两元素向量分别指定低频和高频边缘，buttord同样适用于高通和带阻场景。LP·HP·BP·BSIIRFilterDesign切比雪夫(Chebyshev)滤波器的频率特性切比雪夫I型滤波器以通带内等波纹为代价换取比同阶Butterworth更快的过渡带衰减速率，通带波纹幅度由参数ε控制。01幅度平方函数|Ha(jΩ)|²=1/(1+ε²·CN²(Ω/Ωc))，CN为N阶切比雪夫多项式，ε控制通带波纹幅度。02通带内等波纹起伏，波纹峰值个数与阶数N成正比；阻带内幅频特性单调下降，衰减速率优于同阶Butterworth。03参数ε由通带最大衰减Rp确定：ε=√(10Rp/10−1)，ε越小通带越平坦但过渡带越宽。04当ε=1时，Ωc对应通带的3dB截止频率；通常ε<1，Ωc表示通带边缘频率而非3dB频率。Chebyshev与Butterworth幅频响应对比(N=4)同阶数下Chebyshev过渡带衰减显著快于Butterworth，但通带存在等波纹POLEDISTRIBUTION切比雪夫滤波器极点分布与设计方法N阶Chebyshev滤波器的2N个极点分布在s平面的椭圆上，成复共轭对出现且关于虚轴对称。极点坐标由两个辅助圆通过作图法确定——纵坐标取自大圆、横坐标取自小圆，仅选取左半平面的N个极点构建稳定传递函数，这一分布规律是理解Chebyshev设计流程的关键。01极点分布特征幅度平方函数在s平面上有2N个极点，成复共轭对出现且关于虚轴对称，无极点落在虚轴上2N极点02坐标确定方法极点纵坐标由大圆（半径bΩc）上均匀分布的2N个点确定，横坐标由小圆（半径aΩc）对应点确定大圆·小圆03形状控制参数参数a和b由波纹系数ε决定：a=(1/N)·sh⁻¹(1/ε)，b=(1/N)·ch⁻¹(1/ε)，控制椭圆形状a,b←ε04设计流程设计步骤与Butterworth类似：频率归一化求λp和λs，查归一化曲线确定阶数N，再查表获取传递函数归一化→查表IIRFilterDesign切比雪夫滤波器MATLAB实现与选型对比MATLAB通过cheb1ord和cheby1函数实现ChebyshevI型滤波器的自动化设计。在工程选型中，Butterworth适合通带平坦度要求高的场景，Chebyshev适合阶数受限且需要陡峭过渡带的场景，两者互为补充构成模拟滤波器逼近的基础工具集。ImplementationMATLAB实现[n,Wn]=cheb1ord(Wp,Ws,Rp,Rs,'s')计算最小ChebyshevI型模拟滤波器阶数，Rp为通带波纹（dB）[z,p,k]=cheby1(n,Rp,Wn,'s')生成滤波器零点、极点和增益，支持低通、高通、带通、带阻四种类型cheb1ord→cheby1SelectionCriteria选型对比原则通带平坦度优先→选Butterworth；过渡带陡峭度优先→选Chebyshev，相同指标下Chebyshev阶数更低Chebyshev的通带波纹Rp是设计输入参数，可根据应用容忍度灵活调整以优化阶数和过渡带的平衡Butterworth:FlatChebyshev:SteepChapter03脉冲响应不变法与数字仿真掌握时域采样映射原理，理解脉冲响应不变法的推导过程、频率混叠问题及适用条件PRINCIPLESOFDIGITALSIMULATION模拟系统数字仿真的基本原理模拟滤波器到数字滤波器的转换本质是s平面到z平面的映射，必须满足两个核心约束：虚轴映射到单位圆、左半平面映射到单位圆内。01虚轴→单位圆s平面虚轴（jΩ）映射到z平面单位圆（ejω），实现模拟频率响应与数字频率响应的一一对应关系，保证频域特性的一致性。02左半平面→圆内s平面左半平面（Re[s]<0）映射到z平面单位圆内（|z|<1），确保模拟系统的稳定性在数字域中得以保持，极点位置决定系统稳定。03脉冲响应不变法通过时域采样建立映射关系：h(n)=ha(nT)，保持冲激响应的采样对应关系，适用于要求时域特性匹配的设计场景。04双线性变换法通过代数变换建立映射：s=(2/T)·(1−z⁻¹)/(1+z⁻¹)，在频域建立非线性频率对应，避免混叠现象，是工程中最常用的方法。MATHEMATICALDERIVATION脉冲响应不变法的数学推导脉冲响应不变法通过对模拟冲激响应进行等间隔采样建立模拟与数字滤波器的对应关系，模拟极点pk被指数映射为数字极点e^(pkT)，映射关系简洁且具有明确的物理意义。01基本原理：h(n)=ha(nT)=ha(t)|t=nT，数字滤波器的脉冲响应等于模拟冲激响应的等间隔采样序列，采样周期T决定数字域的频率分辨率。02部分分式展开：Ha(s)=Σ(Ak/(s−pk))，将高阶模拟系统分解为若干一阶子系统的并联叠加，每个子系统仅含一个极点。03极点对应映射：每个一阶项Ak/(s−pk)→Ak/(1−epkT·z⁻¹)，模拟极点pk通过指数函数映射为数字极点epkT。04稳定性保持：左半平面极点Re[pk]＜0映射后|epkT|＜1，数字极点严格位于z平面单位圆内，因果稳定性自动得到保证。IIRFILTERDESIGN脉冲响应不变法的频率混叠与适用条件脉冲响应不变法在时域采样的同时导致频域频谱的周期延拓叠加，产生不可消除的频率混叠。混叠程度取决于模拟滤波器高频衰减速度，因此该方法仅适用于高频衰减充分的低通和带通滤波器设计，不适合高频能量较大的高通和带阻滤波器。脉冲响应不变法的频谱混叠示意频谱周期延拓导致过渡带和阻带出现混叠，高频衰减不足时混叠尤为严重01数字频响H(ejω)=(1/T)·ΣHa(j(ω/T−2πk/T))，是模拟频响的无限移位叠加，存在固有混叠02当模拟滤波器高频衰减充分时（如低通），相邻频谱副本重叠小，混叠误差可控制在工程容许范围内03高通和带阻滤波器在高频段幅值不衰减甚至为常数，频谱移位叠加后混叠严重，频率响应严重失真04减小采样周期T可提高采样频率、拉开频谱副本间距从而减轻混叠，但无法从根本上消除CHAPTER04双线性变换法掌握s-z代数映射原理、频率预畸变技术及无混叠的完整IIR滤波器设计流程IIR数字滤波器设计双线性变换法的映射原理双线性变换通过代数替换s=(2/T)(1-z⁻¹)/(1+z⁻¹)将s平面虚轴整体压缩映射到z平面单位圆，从根本上消除频率混叠。整个左半平面映射到单位圆内保证稳定性。代价是模拟频率Ω与数字频率ω之间呈正切非线性关系，高频段频率刻度被压缩，需通过预畸变校正。01核心映射公式s=(2/T)·(1-z⁻¹)/(1+z⁻¹)，将连续域有理函数直接转换为离散域有理函数02虚轴→单位圆s平面虚轴(jΩ,-∞到+∞)被压缩映射到z平面单位圆一周(e^jω,-π到π)，完整覆盖无遗漏03稳定性保持s左半平面(Re[s]<0)映射为|z|<1（单位圆内），s右半平面映射为|z|>1，稳定性严格保持04频率非线性Ω=(2/T)·tan(ω/2)，模拟频率与数字频率呈正切非线性关系，高频段压缩明显IIRDIGITALFILTERDESIGN双线性变换的频率预畸变技术双线性变换的非线性频率映射会导致数字域截止频率偏移，频率预畸变技术通过在设计前将数字频率指标经Ω=(2/T)tan(ω/2)反变换为模拟频率，补偿非线性压缩效应，确保最终数字滤波器的关键频率点精确落在设计指标要求的位置。非线性频率关系Ω=(2/T)tan(ω/2)导致低频段近似线性、高频段严重压缩，直接使用数字频率设计会产生偏差预畸变操作：将数字通带截止频率ωp和阻带起始频率ωs代入Ω=(2/T)tan(ω/2)得到预畸变后的模拟频率用预畸变后的模拟频率指标设计模拟滤波器Ha(s)，再经双线性变换得到H(z)，关键频率点自动对齐采样周期T的选择不影响最终数字滤波器特性（T在预畸变和变换中互相抵消），通常取T=2简化计算Ω=(2/T)tan(ω/2)频率映射关系低频段近似线性，高频段模拟频率急剧增大，非线性压缩效应显著IIRFILTERDESIGN双线性变换法完整设计流程双线性变换法通过"数字指标→预畸变→模拟设计→双线性变换"的四步流程实现无混叠的IIR数字滤波器设计。该方法适用于低通、高通、带通和带阻所有类型的滤波器，因其流程规范、无混叠风险而成为工程实践中最常用的IIR滤波器设计方法。01确定数字滤波器性能指标ωp·ωs·Rp·Rs确定数字通带截止频率ωp、阻带起始频率ωs、通带衰减Rp和阻带衰减Rs四项核心参数。02频率预畸变Ω=(2/T)tan(ω/2)将数字频率ωp和ωs通过Ω=(2/T)tan(ω/2)转换为模拟频率Ωp和Ωs，补偿非线性映射。03设计模拟原型滤波器N+Ha(s)根据Ωp、Ωs、Rp、Rs确定滤波器阶数N，获取模拟系统函数Ha(s)。04双线性变换s→H(z)将s=(2/T)(1-z⁻¹)/(1+z⁻¹)代入Ha(s)，整理得到数字系统函数H(z)。CHAPTER05数字滤波器变换与实现结构掌握频率变换方法、直接型/级联型/并联型实现结构及有限字长效应分析IIRFILTERDESIGN数字滤波器的频率变换方法数字域频率变换通过z平面上的变量替换将低通原型转换为高通、带通和带阻滤波器。变换函数必须是全通有理函数以保证稳定性，每种变换有特定的参数计算公式，工程实践中先设计低通原型再经变换获得目标滤波器类型，实现设计流程的标准化复用。01低通→低通z⁻¹→(z⁻¹−α)/(1−αz⁻¹)，通过参数α调整截止频率，α=sin[(ωc'−ωc)/2]/sin[(ωc'+ωc)/2]02低通→高通z⁻¹→−(z⁻¹+α)/(1+αz⁻¹)，频率轴翻转实现低高通转换，α=cos[(ωc'+ωc)/2]/cos[(ωc'−ωc)/2]03低通→带通用二阶全通函数替换z⁻¹，将单截止频率扩展为上下两个截止频率定义的通带区间，实现频谱搬移与带宽控制，适用于需要保留特定频段的应用场景04低通→带阻与带通变换对偶，将通带转换为阻带，变换参数的计算基于目标阻带的上下边缘频率，用于滤除特定干扰频段同时保留两侧有用信号FILTERSTRUCTURESIIR滤波器的直接型实现结构直接型结构是IIR滤波器最基础的实现方式，直接I型需M+N个延时单元，直接II型通过合并共享延时链将延时减少到max(M,N)个，但对系数量化误差敏感。TYPEI直接I型结构按H(z)分子分母直接实现：先前馈（零点）部分再反馈（极点）部分，需M+N个延时单元。结构直观易于理解，适合教学演示。延时单元数量较多，硬件实现成本较高，在资源受限的嵌入式系统中需权衡。零极点独立实现，便于分别调试和验证系统响应特性。M+N延时单元TYPEII·CANONIC直接II型（正准型）交换零极点子系统顺序并合并共享延时链，延时单元减少至max(M,N)个，是最经济的直接实现方案。对系数量化误差敏感度高于级联型，高阶时可能导致极点偏移甚至不稳定，需配合高精度运算。节省存储资源，在DSP芯片和FPGA实现中广泛应用，是工程实践的首选结构。max(M,N)延时单元IIRFILTERSTRUCTURESIIR滤波器的级联型与并联型结构级联型将高阶系统分解为二阶节串联，量化误差被局限在单个二阶节内不扩散，是DSP芯片中最常用的实现结构。并联型将系统展开为二阶节并联求和，各支路可并行运算适合高速场景。两者均显著优于直接型的量化误差性能，是工程实现的推荐方案。级联型结构01将H(z)分解为若干二阶节级联：H(z)=H₁(z)·H₂(z)···Hₖ(z)，每节独立实现一对零极点02量化误差局限在单个二阶节内不扩散，零极点可独立调整，是实际工程中最常用的IIR实现结构03二阶节顺序和增益分配影响整体动态范围，需优化排列以获得最佳信噪比并联型结构01将H(z)展开为部分分式之和：H(z)=H₁(z)+H₂(z)+···+Hₖ(z)，各支路独立运算后求和02各支路可并行执行，总延迟等于单个二阶节的延迟，适合对处理速度要求高的实时应用场景03极点量化误差被局限在各支路内，但零点由各支路共同决定，零点精度不如级联型FINITEWORDLENGTHEFFECTS有限字长效应对IIR滤波器实现的影响有限字长效应包括系数量化导致的零极点偏移、运算舍入产生的量化噪声和定点溢出三种主要形式。IIR滤波器因存在反馈结构对量化误差尤为敏感，系数量化可能使极点移出单位圆导致系统不稳定。采用级联/并联结构、增加运算字长和引入溢出保护是工程中的主要对策。系数量化效应有限位表示导致滤波器系数偏离设计值，极点位置偏移可能使系统不稳定（极点移出单位圆）极点偏移运算舍入噪声每次乘法运算后的舍入操作等效于注入白噪声，经反馈结构放大后影响输出信噪比量化噪声定点溢出问题累加器数值超出寄存器范围导致严重非线性失真，需采用饱和运算或增加保护位防止溢出非线性失真缓解策略选用级联/并联型结构降低敏感度、增加运算字长、合理安排二阶节顺序和增益分配级联/并联CHAPTER06优化设计与MATLAB实践掌握IIR滤波器优化策略、MATLAB综合设计流程及性能分析方法IIRFilterOptimizationIIR滤波器的优化设计策略IIR滤波器优化设计通过计算机辅助迭代算法在满足频率响应指标的前提下寻求最低阶数或最优系数组合。常用方法包括零极点位置优化和系数直接优化两类路径，优化目标可以是最大误差最小化(Minimax)或加权均方误差最小化，适用于标准逼近方法无法满足特殊约束的场景。零极点优化法在z平面上直接调整零极点位置，以频率响应误差为目标函数进行迭代搜索最优配置z平面迭代系数优化法直接优化系统函数的分子分母系数，常用Fletcher-Powell共轭梯度法和准牛顿法共轭梯度Minimax准则最小化通带和阻带中的最大偏差（最坏情况），保证频率响应在所有频点满足指标最坏情况加权优化对通带和阻带误差施加不同权重，灵活平衡各频段的逼近精度以满足差异化需求差异化权重IIRFilterDesignMATLABIIR滤波器综合设计实例通过MATLAB实例串联IIR滤波器设计全流程——从指标定义、阶数计算到