发明名称 一种用于消除声学回声的自适应滤波器及滤波方法
摘要 本发明提出了一种用于消除声学回声的自适应滤波器,所述自适应滤波器采用基于以下目标函数的自适应滤波算法,并采用限制系数对滤波器更新进行限制,且本发明的自适应滤波器每迭代R次对δ(k,m)进行一次更新。本发明克服已有的变步长控制方法的不足,提出一种收敛性能较好、稳态误差小且对双讲情况鲁棒的步长控制方法。该方法基于频域滤波器实现,通过对每次迭代过程中滤波器系数更新进行约束,获得较小的稳态误差和较好的双讲鲁棒性能,同时在频域对约束条件进行更新,提高了滤波器的收敛速度和跟踪速度。
申请公布号 CN104883462A 申请公布日期 2015.09.02
申请号 CN201410073711.1 申请日期 2014.02.28
申请人 中国科学院声学研究所;北京中科信利技术有限公司 发明人 吴超;付强;颜永红
分类号 H04M9/08(2006.01)I;G10L21/0232(2013.01)I 主分类号 H04M9/08(2006.01)I
代理机构 北京法思腾知识产权代理有限公司 11318 代理人 杨小蓉
主权项 一种用于消除声学回声的自适应滤波器,所述自适应滤波器采用基于以下目标函数的自适应滤波算法:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>W</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mrow><mi>W</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>Z</mi></mrow></munder><msup><mrow><mo>|</mo><msub><mi>E</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mo>|</mo><mi>W</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>-</mo><mi>W</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mo>&le;</mo><mfenced open=' ' close=''><mtable><mtr><mtd></mtd></mtr><mtr><mtd></mtd></mtr></mtable></mfenced><mfenced open='{' close=''><mtable><mtr><mtd><mi></mi><mn>2</mn><mi></mi><msup><mi>&delta;</mi><mi>l</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mtd><mtd><mi>m</mi><mo>&le;</mo><msub><mi>N</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mi></mi><mn>2</mn><mi>&delta;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo></mtd><mtd><mi>m</mi><mo>></mo><msub><mi>N</mi><mn>0</mn></msub></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000471276440000011.GIF" wi="1294" he="299" /></maths>其中,m表示频率点,δ(k,m)和δ<sup>l</sup>(k,m)是滤波器系数限制参数,k表示数据块数,arg表示定义域,E<sub>P</sub>()表示频域后验误差,|E<sub>P</sub>(k,m)|<sup>2</sup>表示第k块第m频点的后验误差,Z表示复空间,s.t.表示“需满足条件”;W(k,m)表示更新前滤波器权值,W(k+1,m)表示更新后滤波器权值;X(k,m)表示频域参考信号;N<sub>0</sub>为整数,且满足<img file="FDA0000471276440000012.GIF" wi="336" he="138" />fs表示声卡采样率且单位为Hz,[]表示取整数运算。
地址 100190 北京市海淀区北四环西路21号