免提通话系统中的稳健回声抵消方法
【专利摘要】一种免提通话系统中的稳健回声抵消方法,其特征在于：该方法包括以下步骤：步骤一、获得一个采样序列，步骤二、获得误差信号：步骤三、权系数的更新：步骤四、对上述步骤一至步骤三重复执行50~100次，然后计算所有滤波器输出信号的平均值，将该平均值作为最终消除掉回声的纯净信号。其能够准确估计出回声路径，得到纯净的语音信号，且计算复杂度低。
【专利说明】免提通话系统中的稳健回声抵消方法
【技术领域】 [0001] ：本发明涉及一种针对免提通话系统的回声消除技术，基于稳健估计理 论对具有脉冲特征的噪声进行非线性抑制，采用自适应滤波方法实现对通话系统路径的辨 识，属于信号处理。
【背景技术】 [0002] ：回声抵消算法是会议电视、车载电话等免提通话系统中实现全双工通 话的关键技术。通常由扬声器、封闭环境和麦克风所构成的回声路径冲击响应是未知、复杂 且随时间变化的，所以无法设计和使用固定的回声抵消器得到对所有环境都适合的回声抵 消器。根据回声路径冲击响应的特点，有效地阻止和控制回声的方法是自适应回声抵消，其 基本思想是：用一个估计的回声路径冲击响应去代替实际的回声路径冲击响应，以产生与 回声具有相同特性的信号，然后从实际信号中将它减去，实现抵消。目前，己经提出了多种 自适应回声抵消方法，并有基于高速数字信号处理芯片的实现方案。
[0003] 近年来，随着无线接入网、卫星网等的广泛使用，通话的距离越来越长，信号传输 延迟也大大增加，从而使得现代免提通话系统具有一种特殊的特性--稀疏性，即系统的 大部分能量只集中在脉冲响应的一小部分，多数权系数为〇或很小的数值，只有少数权系 数较大。系统的这种稀疏特性增加了基于一般自适应滤波算法的回声抵消器的有效处理时 间。
[0004] 但目前绝大多数方法不断努力追求的目标主要是如何提高收敛速度，增强其对环 境稀疏性的适应性，减少计算复杂度和失调量。而对影响语言质量的另一重要因素一随 机噪声只作了简单的高斯假设，即只考虑了噪声概率密度函数轻薄拖尾的情况，从而其代 价函数或是基于NLMS的均方误差准则，或是基于最小二乘准则，或是基于距离定义。当然， 根据中心极限定理，这个假设是合理的。实际上，正是这个假设给信号处理分析带来了极大 的方便。比如：信号处理模型存在简易的封闭表达式，而且这些表达式通常都是线性的；只 需要均值和方差两个统计量就足以表达等等。
[0005] 实际上，在很多免提通话系统中，与回声同时存在且影响语音质量的背景噪声并 不呈现高斯特性。例如，在车内使用免提电话时，存在着汽车发动机点火引起的电磁噪声、 车外气流噪声、轮胎摩擦噪声等。而在电视会议系统中，室内回声可能混杂来自空调运行噪 声等。这些噪声都表现出比高斯噪声显著得多的脉冲性，因此对环境噪声仅仅做出高斯分 布的假设并不完全符合真实情况。
[0006] 由稳健统计理论可知，前面所述算法中对代价函数的运算都可以归结为二范数的 范畴，它们对不符合高斯分布的脉冲噪声和信号具有放大作用，因而使得自适应回声抵消 算法中的梯度估计偏差很大，继而影响权向量的搜索方向和相关矩阵的估计。因此，存在于 免提通话系统中的非高斯脉冲性干扰和噪声必然对现有基于高斯分布假设的算法造成不 良影响，从而限制其实用性。
【发明内容】
：
[0007] 发明目的：本发明提供一种免提通话系统中的稳健回声抵消方法，其目的是解决 以往的方式所存在的不能准确估计出回声路径和语音信号，且计算复杂度高的问题。
[0008] 技术方案：本发明是通过以下技术方案来实现的：
[0009] -种免提通话系统中的稳健回声抵消方法，其特征在于：该方法包括以下步骤：
[0010] 步骤一、获得一个采样序列X (η)，其中η = 1，2,…，N，N为采样序列X (η)的采样 点数量；所述的采样序列X (η)是一维信号，且其中包含N个采样点；
[0011] 步骤二、获得误差信号：采用所得到的采样序列x(n)，获得理想的输出信号，并计 算实际接收信号和理想接收信号的误差信号e (η);
[0012] 步骤三、权系数的更新：应用零吸引-最小平均Ip范数算法，在每一离散时间点上 进行迭代，对未知的回声路径h进行估计和更新调整，计算出自适应滤波器的输出士/7);当 自适应滤波器收敛后，该输出信号即为回声信号的一个复制，将其从期望信号d(n)中减去 便可消除回声；
[0013] 步骤四、对上述步骤一至步骤三重复执行50?100次，然后计算所有滤波器输出 信号的平均值，将该平均值作为最终消除掉回声的纯净信号。
[0014] 步骤三中应用零吸引-最小平均Ip范数算法是在每一离散时间点上进行迭代，其 更新过程为：
[0015] e(n) = din)-%1 (ti)h(n -1) (I)
[0016] h (η) = h (n_l) - κ sgn {h (η)} + μ I e (n) I p-1sgn [e (η) ] x (η) (2)
[0017] 其中，?〇)是对未知的回声路径h进行估计和更新调整，其中K是一个比较小的 正数，0 < μ < 1是每个独立自适应滤波器的步长，e(n)是误差信号，d(n)是期望信号，L 是滤波器长度，x(n) = [x(n)，x(n-l),…，x(n-L+l)]T。
[0018] 优点及效果：
[0019] 本发明提供一种免提通话系统中的稳健回声抵消方法，其能够准确估计出回声路 径，得到纯净的语音信号，且计算复杂度低。
[0020] 本发明与现有技术相比，具有如下优点：
[0021] (1)、本发明采用Ip范数作为代价函数，可以大大降低脉冲噪声的冲击幅度，减少 其对回声消除效果的影响。
[0022] (2)、本发明充分应用了系统的稀疏特性，利用零吸引算子对系统权系数进行约 束，使得本发明在收敛速度和估计精度方面达到了更好的折中，并大大减少了计算复杂度。
[0023] (3)、本发明能够准确有效地估计出免提通话系统参数，消除各种回声，并通过合 理的参数选择有效抑制各类脉冲噪声，能够有效提高免提通话系统的稳健性和适应性
【专利附图】
【附图说明】：
[0024] 图1是本发明的免提通话系统的自适应回声消除器的原理框图；
[0025] 图2是本发明的应用流程图；
[0026] 图3是本发明实施例中用到的回声路径，（a)稀疏度为0. 69的实测室内信道（b) 稀疏度为〇. 49的仿真信道；
[0027] 图4是在非高斯噪声条件下，本发明的性能图；
[0028] 图5是在高斯噪声条件下，本发明的性能图。
【具体实施方式】 [0029] ：下面结合附图对本发明做进一步的描述：
[0030] 如图1所示，本发明提供一种免提通话系统中的稳健回声抵消方法，其利用稳健 估计思想，将原代价函数中的I 2范数用Ip，通过对参数P的选择，实现了对各种脉冲程度不 同的噪声的抑制。同时，在迭代过程中，采用零吸引算子实现对稀疏系统中非零值的有效利 用，达到了减少计算量的效果。具体地，本发明包括如下步骤：
[0031] 步骤一、获得一个采样序列X (η)，其中η = 1，2,…，N，N为采样序列X (η)的采样 点数量；所述的采样序列X (η)是一维信号，且其中包含N个采样点；
[0032] 步骤二、获得误差信号：采用所得到的采样序列χ(η)，获得理想的输出信号，并计 算实际接收信号和理想接收信号的误差信号e (η);
[0033] 步骤三、权系数的更新：应用零吸引-最小平均Ip范数算法，在每一离散时间点上 进行迭代，对未知的回声路径h进行估计和更新调整，计算出自适应滤波器的输出士/7);当 自适应滤波器收敛后，该输出信号即为回声信号的一个复制，将其从期望信号d(n)中减去 便可消除回声；
[0034] 应用零吸引-最小平均Ip范数算法是在每一离散时间点上进行迭代，其更新过程 为：
[0035] e(n) = d(n)- X1 (n)h(n -1) (I)
[0036] h (η) = h (n_l) - κ sgn {h (η)} + μ I e (η)卜1sgn [e (η) ] x (η) (2)
[0037] 按照公式：
【权利要求】
1. 一种免提通话系统中的稳健回声抵消方法，其特征在于：该方法包括以下步骤： 步骤一、获得一个采样序列X(n)，其中n= 1，2,…，N，N为采样序列X(n)的采样点数 量；所述的采样序列x(n)是一维信号，且其中包含N个采样点； 步骤二、获得误差信号：采用所得到的采样序列x(n)，获得理想的输出信号，并计算实 际接收信号和理想接收信号的误差信号e(n); 步骤三、权系数的更新：应用零吸引-最小平均lp范数算法，在每一离散时间点上进行 迭代，对未知的回声路径h进行估计和更新调整，计算出自适应滤波器的输出士/7);当自适 应滤波器收敛后，该输出信号即为回声信号的一个复制，将其从期望信号d(n)中减去便可 消除回声； 步骤四、对上述步骤一至步骤三重复执行50?100次，然后计算所有滤波器输出信号 的平均值，将该平均值作为最终消除掉回声的纯净信号。
2. 根据权利要求1所述的免提通话系统中的稳健回声抵消方法，其特征在于： 步骤三中应用零吸引-最小平均lp范数算法是在每一离散时间点上进行迭代，其更新 过程为：
h(n) =h(n-1) -ksgn{h(n)} +u|e(n) |p_1sgn[e(n) ]x(n) (2) 其中，是对未知的回声路径h进行估计和更新调整，其中k是一个比较小的正数， 0 <ii< 1是每个独立自适应滤波器的步长，e(n)是误差信号，d(n)是期望信号，L是滤 波器长度，x(n) = [x(n)，x(n-l)，…，x(n-L+l)]T。
【文档编号】H04M9/08GK104410762SQ201410659854
【公开日】2015年3月11日 申请日期:2014年11月18日 优先权日:2014年11月18日
【发明者】郭莹, 白艳梅 申请人:沈阳工业大学