发明名称 一种频域自适应回声抵消方法
摘要 提供了一种频域自适应回声抵消的方法,在所述方法中:从扬声器输入信号分块,变换到频域,频域滤波,然后反变换回时域,与传声器信号相减,得到了误差信号,利用误差信号进行最优步长的估计以及规整系数的估计,从而更新滤波器的系数。以较小的计算量,实现了回声抵消的功能。
申请公布号 CN101888455A 申请公布日期 2010.11.17
申请号 CN201010142853.0 申请日期 2010.04.09
申请人 熔点网讯(北京)科技有限公司 发明人 吴文焘
分类号 H04M9/08(2006.01)I;G10L21/02(2006.01)I;G10L19/02(2006.01)I 主分类号 H04M9/08(2006.01)I
代理机构 代理人
主权项 1.一种频域自适应回声抵消方法,包括:步骤一:把输入信号x(n)进行分块,块大小为M,对分块的信号进行FFT变换,把信号变换到频域,得到变换后的频域信号X<sub>m</sub>(k);步骤二:对变换后的频域信号X<sub>m</sub>(k)进行自适应滤波处理,得到估计信号<img file="FSA00000062336500011.GIF" wi="98" he="59" />把估计信号<img file="FSA00000062336500012.GIF" wi="74" he="58" />与参考信号y(n)进行比较,得到误差信号e(n);步骤三:结合输入信号x(n)和误差信号e(n),估计最优步长,具体步骤为:首先把误差信号e(n)与估计信号<img file="FSA00000062336500013.GIF" wi="75" he="62" />变换到频域,得到E(k)与<img file="FSA00000062336500014.GIF" wi="102" he="55" />然后计算E(k)与<img file="FSA00000062336500015.GIF" wi="81" he="55" />功率的互相关<img file="FSA00000062336500016.GIF" wi="606" he="45" />和<img file="FSA00000062336500017.GIF" wi="82" he="55" />功率的自相关<img file="FSA00000062336500018.GIF" wi="637" he="45" />其中<img file="FSA00000062336500019.GIF" wi="93" he="43" />为<img file="FSA000000623365000110.GIF" wi="82" he="55" />的功率,P<sub>E</sub>(k)为E(k)的功率,γ为平滑系数,<img file="FSA000000623365000111.GIF" wi="111" he="45" />和<img file="FSA000000623365000112.GIF" wi="114" he="45" />取随机数初始化;接着计算线形回归系数<img file="FSA000000623365000113.GIF" wi="247" he="162" />和回声误差信号的方差<img file="FSA000000623365000114.GIF" wi="221" he="54" />其中<img file="FSA000000623365000115.GIF" wi="649" he="81" />β为平滑系数,<img file="FSA000000623365000116.GIF" wi="124" he="63" />初始化取随机数;计算误差信号方差<img file="FSA000000623365000117.GIF" wi="630" he="69" />其中<img file="FSA000000623365000118.GIF" wi="122" he="61" />初始化取随机数;最后计算最优步长的估计:<img file="FSA000000623365000119.GIF" wi="337" he="111" />步骤四:利用变换后的频域信号X<sub>m</sub>(k),估计优化步长规整系数:P<sub>m</sub>(k)=E{max{ξ<sub>m</sub>(k),|X<sub>m</sub>(k)|<sup>2</sup>}},其中<img file="FSA000000623365000120.GIF" wi="627" he="80" />滤波的输出,<img file="FSA000000623365000121.GIF" wi="869" he="168" />步骤五:利用输入信号、误差信号、最优步长、步长规整系数对滤波器权重W<sub>m</sub>(k)进行更新:<img file="FSA00000062336500021.GIF" wi="651" he="92" />其中X<sub>m</sub><sup>*</sup>(k)表示输入信号频率点k的值的共轭。
地址 100081 北京市海淀区中关村南大街3号海淀科技大厦901室