数字信号处理：几种常见数字滤波器实现原理详解引用CSDN1.https://blog.csdn.net/qq_42947060/article/details/116915181
数字信号处理中的数字滤波器是信号处理领域的重要工具，用于从含有噪声的数据记录中估计出信号的某些特征或信号本身。本文将详细介绍几种常见的数字滤波器实现原理，重点讲解巴特沃斯滤波器的设计过程，包括模拟滤波器设计和数字化转换的关键步骤。
现代滤波器理论研究的主要内容是从含有噪声的数据记录中估计出信号的某些特征或者信号本身，估计出的信号的信噪比将比原信号的高，将信号和噪声都视为随机信号，利用其统计特征导出一套最佳的估计算法。
设计的步骤：
给出所需的滤波器的技术指标；
设计一个H(z)使其逼近所需的技术指标；
实现所设计的H(z)，其中步骤2是本章和下一章所讨论的主要内容。
按功能分类按相位响应分为线性相位和非线性相位；按特殊要求分为最小相位滞后滤波器、梳状滤波器等；按冲激响应分为无限冲激响应数字滤波器(IIR：有反馈，无限长)，有限冲激响应数字滤波器(FIR：无反馈，有限长)
示例：巴特沃斯滤波器
滤波器设计
目的：找到一个满足技术指标要求的可实现的因果稳定的数字滤波器来逼近理想的滤波器幅度特性。
思路：间接法-先设计模拟滤波器，然后通过双线性变换法来制作数字低通滤波器、数字高通滤波器、数字带通滤波器、数字带阻滤波器
模拟滤波器设计：通带截止频率wp，阻带截止频率ws， 通带最大衰减1-δ1（Rp(dB)）， 阻带最小衰减δ2（As(dB)） 。
1、确定技术指标：Ωp or s=2∗π∗fp or s，1−δ1,δ2，Ωc（3dB频率）Ωc≥lg(100.1As−1100.1Ap−1)=ΩcpΩc≤lg(Ωs100.1As−1)=Ωcs
2、计算滤波器所需阶数N，CN=lg(100.1As−1100.1Ap−1)/(2∗lg(λs))C2=10αp/10−1
3、查表确定归一化低通滤波器系统函数G(p)G(p)=1(p−p1)(p−p2)(p−p3)...(p−pN)
因为几点成对共轭出现，所以可以按N的奇偶性如下表示：N为偶G(p)=∏K=1N/2Gk(p)N为偶G(p)=K=1∏N/2Gk(p)N为奇G(p)=1p+1∏K=1(N−1)/2Gk(p)N为奇G(p)=p+11K=1∏(N−1)/2Gk(p)
Gk(p)=1p2−2pcos(2k+N−12Nπ)+1极点：pk=exp(j2k+N−12Nπ) k=1,..N
4、将Has(s)转换为所需类型的低通滤波器系统函数G(s)，因为p=sΩc，所以G(s)=G(sΩc)
关于G§还有一种查表计算法：G(p)=1(p−p1)(p−p2)(p−p3)...(p−pN)=11+a1p+a2p2+...+aNpN
模拟滤波器数字化方法有冲激响应不变法、阶跃响应不变法和双线性变换法等。冲激响应不变法(适合低通和带通);双线性变换法：s->z
简单说来，将水平面的转移函数映射到z平面，就是数字化，变换公式在思路流程图中。
实际应用滤波的过程就是解常系数线性差分方程的过程。y(n)为滤波后信号，x(n)为滤波前信号y(n)=∑m=0Mbm x(n−m)−∑k=1Nak y(n−k)
其中ak,bm分别为系统函数H(z)分母与分子的系统数组
滤波器的素质转模拟技术指标的规则
低通截止频率fp,阻带起始频率fs,通带最大衰减αp,阻带最小衰减αs，由此计算N值λp=1λs=(fs/fp)p=sΩp
高通截止频率fp,阻带起始频率fs,通带最大衰减αp,阻带最小衰减αs，由此计算N值λp=1λs=(fp/fs)p=sΩp
带通Ω=ωf通带下限截止频率ω1通带上限截止频率ω3下阻带截止频率ωzl上阻带截止频率ωzh通频带宽ωBW=ω3−ω1定义中心频率ω22=ω1ω3
将频率归一化处理：ηzl=ωzl/ωBW;ηzh=ωzh/ωBW;η=ω1/ωBW;η=ω3/ωBW
因为η3−η1=1=>λp=1λs=ηzh2−η22ηzh−λs=ηzl2−η22ηzl最终λs取上绝对值的最小值p=s2+Ω1Ω3s(Ω3−Ω1)
带阻λs=ηzhηzh2−η22−λs=ηzlηzl2−η22p=s(Ω3−Ω1)s2+Ω1Ω3最终λs取上绝对值的最小值
G(p)公式
巴特沃斯(Butterworth)∣G(jΩ)∣2=11+C2(Ω2)NG(p)=1(p−p1)(p−p2)...(p−pN)C=10αp/10−1N=lg(100.1As−1100.1Ap−1)/(2∗lg(λs))
切尔雪夫I型(Chebyshev-I)∣G(jΩ)∣2=11+ε2Cn2(Ω)G(p)=1ε×2n−1∏k=1n(p−pk)ε2=10αp/10−1a2=100.1As−1100.1Ap−1N=arcoshaarcoshλs注：arcosh(x)=ln(x+sqrt(x2−1))
切尔雪夫II型(Chebyshev-II)∣G(jΩ)∣2=11+ε2[Cn2(Ωs)Cn2(Ωs/Ω)]2
椭圆滤波器∣G(jΩ)∣2=11+ε2⋃n2(Ω)