[发明专利]一种带宽可变低通数字滤波器的设计方法

摘要

本发明提供了一种带宽可变低通数字滤波器的设计方法，通过将Farrow结构的FIR数字滤波器各子滤波器的系数C作为变量，在设计准则下，通过迭代运算，以得到最优的各子滤波器的系数C。这样在保持低通数字滤波器不变的前提下，只要修改输入期望截止频率ωc_des，就可以改变低通数字滤波器截止频率参数ωc即带宽可变，提高了低通数字滤波器使用的灵活性，避免了在线设计或存储大量滤波器系数表的复杂性。

主权项

1.一种带宽可变低通数字滤波器的设计方法，其特征在于，包括以下步骤：（1）、带宽可变低通数字滤波器设计准则准则1：截止频率准则Δc=maxωL<ωc<ωu|ωc-ωc_des|≤δc---(4)]]>准则2：幅度准则Δap=max0<ω<ωc-ωp||H(C,ejω,ωc)|-1|≤δp]]>Δas=maxωc+ωs<ω<π|H(C,ejω,ωc)|≤δs---(5)]]>其中，ω_L为截止频率下限、ω_H为截止频率的上限，ω_c为实际截止频率，ω_{c_des}为输入的期望截止频率，δ_c为截止频率的误差限；H(C,e^jω,ω_c)为带宽可变低通数字滤波器的传递函数，C为带宽可变低通滤波器各子滤波器的系数，ω为数字频率，δ_p为频带[0,ω_c-ω_p]内幅度响应误差限，δ_s为阻带内幅度响应误差限，ω_c-ω_p为通带增益第一次低于幅度响应误差限δ_p时所对应的频率点；（2）、带宽可变低通数字滤波器阶数P₂的确定a)、子滤波器C₀(z)在频带[0,ω_c-ω_p]内幅度响应为1，阻带内幅度响为0，期望截止频率为ω_{c_des}=ω_L，设置子滤波器C₀(z)的初始阶数N₀；b)、利用频率抽样，根据幅度以及相位确定子滤波器C₀(z)的N₀个频响样值，再利用IFFT（傅里叶反变换），计算得到子滤波器C₀(z)的系数c₀(n),n=0,...,N₀-1；c)、根据系数c₀(n),n=0,...,N₀-1，计算子滤波器C₀(z)的实际频响即传递函数H(C,e^jω,ω_c)以及截止频率参数ω_c；d)、如果子滤波器C₀(z)未满足约束条件式（4）、（5）即准则1、2的要求，则更新子滤波器C₀(z)阶数N₀=N₀+1，重复步骤b)、c)，直至满足为止，则此时子滤波器C₀(z)阶数N₀即为带宽可变低通数字滤波器阶数P₂；（3）、带宽可变低通滤波器各子滤波器的系数C的确定a)、子滤波器的个数P₃=0；b)、子滤波器的个数加1即P₃=P₃+1，根据步骤（2）得到的带宽可变低通数字滤波器阶数P₂，按照FIR数字滤波器的泰勒展开式：Hdes(z,ωc_des)≈Σn=0P2-1[Σl=0P3-1cl(n)(ωc_des-ωL)l]z-n]]>=Σl=0P3-1[Σn=0P2-1cl(n)z-n](ωc_des-ωL)l---(6)]]>=Σl=0P3-1Cl(z)(ωc_des-ωL)l]]>写出传输函数H_des(z,ω_{c_des})，然后根据傅里叶变换，得到传递函数H(C,e^jω,ω_c)，其中，带宽可变低通滤波器各子滤波器的系数C为：C=[c0(0),c0(1),...,c0(P2-1),]]>c1(0),c1(1),...,c1(P2-1),...]]>cP3-1(0),cP3-1(1),...,cP3-1(P2-1)]]]>c)、频带分成I等份，即ω_i∈[0,π],i=1,2,...,I；截止频率变化范围为其中的J等份，ω_cj∈[ω_L,ω_U],j=1,2,...,J，将约束条件式（4）、（5）即准则1、2写为：ϵc=max0≤j≤J-1|ωcj-ωcj_des|-δc≤0---(7)]]>ε_a=(Δ_ap-δ_p)+(Δ_as-δ_s)≤0    （8）其中：为取整算子；使用Minimax（最大误差最小化）非线性优化算法，在非线性约束条件式（8）控制下，利用序列二次规划算法(sequential quadratic programming，SQP)，寻找使得式（7）最小的带宽可变低通滤波器各子滤波器的系数C，即性能函数为J=min{max0≤j≤J-1|ωcj-ωcj_des|-δc}---(10)]]>如果最小的带宽可变低通滤波器各子滤波器的系数C存在，且满足约束条件式（7）、（8），则实现了带宽可变低通滤波器，否则返回步骤(b)，增加子滤波器的个数，重新寻找带宽可变低通滤波器各子滤波器的系数C。