【发明授权】改进的两通道IIR的QMFB设计方法_杭州电子科技大学_201911167427.X 

申请/专利权人:杭州电子科技大学

申请日:2019-11-25

发明/设计人:王浩;李祥振;赵知劲;赵晨子;李伟琪;靳一

公开(公告)日:2020-09-15

代理机构:杭州君度专利代理事务所(特殊普通合伙)

公开(公告)号:CN111010144B

代理人:杨舟涛

主分类号:H03H17/02(20060101)

地址:310018 浙江省杭州市下沙高教园区2号大街

分类号:H03H17/02(20060101)

优先权:

专利状态码:有效-授权

法律状态:2020.09.15#授权;2020.05.08#实质审查的生效;2020.04.14#公开

摘要:本发明公开了一种改进的两通道IIR的QMF组设计方法,本发明针对全通滤波器的相位和整体失真传输函数的相位进行最小最大化问题,通过确定相应的相位误差和合适的加权值,对得到的非线性优化目标进行一阶泰勒展开,转换为线性问题,以求得最优的全通滤波器的系数。本发明最大的改进在于不再仅仅只考虑全通滤波器的性能或者只考虑分析‑合成滤波器的性能,而是把两者的性能综合考虑,保证了信号经过的每一个子滤波器都具有更近似线性的相位,减少相位失真的可能性,达到设计预期目标。而且本发明并未对幅度讨论,只对全通滤波器和整体传输函数的相位进行了讨论,而幅度完全由相位控制。

主权项:1.改进的两通道IIR的QMFB设计方法,其特征在于,该方法具体包括以下步骤:步骤一、根据设计要求,确定全频带上的频率点数L、两个全通滤波器的阶数N1和N2、分析滤波器中的低通滤波器的通带截止频率ωp,阻带截止频率ωs,令迭代初始系数k=0,第k次的全通滤波器系数aik=0,初始加权值Wi=1,i=1,2;其中N1=N2+1;步骤二、确定全通滤波器的实际相位误差及整体失真传输函数Tejω的实际相位误差;2.1.确定全通滤波器的理想相位全通滤波器的理想相位在ω∈[0,π]的频率带上满足:当IIR滤波器的相位满足ω=0时,相位是0;当ω=π时,相位满足θπ=-Nπ;此时全通滤波器是稳定的;全通滤波器的理想相位分别是θd1ω=-N1ω+0.25ω和θd2ω=-N2ω-0.25ω;2.2.求全通滤波器的实际相位误差通过1、2、3式分别得到全通滤波器的表达式、全通滤波器的实际相位表达式和全通滤波器的实际相位误差,其中ω∈[0,2ωp],i=1,2; θeiω=θiω-θdiω3其中ain表示滤波器系数ai的第n个元素,n=1,2,…,Ni;2.3.求整体失真传输函数Tejω的实际相位误差θeTω式4、式5、式6分别为整体失真传输函数的理想相位、实际相位和实际相位误差,其中ω∈[0,ωp];θdTω=θd12ω+θd22ω-ω4θTω=θ12ω+θ22ω-ω5θeTω=θTω-θdTω6其中θ12ω、θ22ω、θd12ω、θd22ω分别表示为式7,式8,式9和式10 θd12ω=-2N1ω+0.5ω9θd22ω=-2N2ω+0.5ω10步骤三、求解得出第k次迭代下全通滤波器的滤波器系数aik;3.1.得到目标优化函数表示为minW1||θe1ω||+W2||θe2ω||+||θeTω||11式11中的θe1ω、θe2ω和θeTω依次是指全通滤波器A1ejω、A2ejω和整体失真传输函数Tejω的实际相位误差;该目标优化函数是最大最小化问题,同时也是一个非线性的问题,步骤3.2可将该非线性问题转换为线性问题;3.2.计算第k-1次迭代所得的θeiω关于系数aik-1的一阶偏导数如式12所示; 3.3.计算第k次迭代时全通滤波器Aiejω的相位误差θeiω如式13所示; 式13中,△i代表第k次迭代时全通滤波器Aiejω的系数增量,Δi=aik-aik-1;计算第k次迭代时整体失真传输函数Tejω的实际相位误差θeTω如式14所示; 3.4.第k次迭代时优化目标函数表示为凸优化问题进行求解,确定第k次迭代的系数增量△1和△2; 步骤四、若式16不成立,则转到步骤五;若式16成立,则将aik作为最终设计出的全通滤波器Aiejω的系数,迭代结束; 式16中,Ek是第k次迭代中目标函数的最大值,Ek-1是第k-1次迭代中目标函数的最大值;μ为设定阈值;Ek=maxW1|θe1a1k,ω|+maxW2|θe2a2k,ω|+max|θeTak,ω|17步骤五、根据群延迟误差的包络计算加权值Wi首先计算出第k次迭代的群延迟误差其中gdi表示第k次迭代的全通滤波器Aiejw的实际群延迟,τdi表示全通滤波器Aiejω的理想群延迟;然后计算出的包络如果不满足使然后令k=k+1,返回步骤三;否则的话直接令k=k+1,返回步骤三;其中ε为设定阈值。

全文数据:

权利要求:

百度查询: 杭州电子科技大学 改进的两通道IIR的QMFB设计方法