解密回声消除技术之一(理论篇)

简介:
一、前言
因为工作的关系,笔者从 2004 年开始接触回声消除 (Echo Cancellation) 技术,而后一直在某大型通讯企业从事与回声消除技术相关的工作,对回声消除这个看似神秘、高端和难以理解的技术领域可谓知之甚详。
要了解回声消除技术的来龙去脉,不得不提及作为现代通讯技术的理论基础—— 数字 信号处理理论。 首先,数字信号处理理论里面有一门重要的分支,叫做自适应信号处理。而在经典的教材里面,回声消除问题从来都是作为一个经典的自适应信号处理案例来讨论的。既然回声消除在教科书上都作为一种经典的具体的应用,也就是说在理论角度是没有什么神秘和新鲜的,那么回声消除的难度在哪里?为什么提供回声消除技术(不管是芯片还是算法)的公司都是来自国外?回声消除技术的神秘性在哪里?
 
二、回声消除原理
从通讯回音产生的原因看,可以分为声学回音( Acoustic Echo )和线路回音( Line Echo ),相应的回声消除技术就叫声学回声消除( Acoustic Echo Cancellation AEC )和线路回声消除( Line Echo Cancellation, LEC )。声学回音是由于在免提或者会议应用中,扬声器的声音多次反馈到麦克风引起的(比较好理解);线路回音是由于物理电子线路的二四线匹配耦合引起的(比较难理解)。
回音的产生主要有两种原因:
1.   由于空间声学反射产生的声学回音(见下图):


图中的男子说话,语音信号( speech1 )传到女士所在的房间,由于空间的反射,形成回音 speech1(Echo) 重新从麦克风输入,同时叠加了女士的语音信号( speech2 )。此时男子将会听到女士的声音叠加了自己的声音,影响了正常的通话质量。此 时在 女士所在房间应用回音抵消模块,可以抵消掉男子的回音,让男子只听到女士的声音。
2.   由于 2-4 线转换引入的线路回音(见下图):

 
ADSL Modem 和交换机上都存在 2-4 线转换的电路,由于电路存在不匹配的问题,会有一部分的信号被反馈回来,形成了回音。如果在交换机侧不加回音抵消功能,打电话的人就会自己听到自己的声音。
不管产生的原因如何,对语音通讯终端或者语音中继交换机需要做的事情都一样:在发送时,把不需要的回音从语音流中间去掉。
试想一下,对一个至少混合了两个声音的语音流,要把它们分开,然后去掉其中一个,难度何其之大。就像一瓶蓝墨水和一瓶红墨水倒在一起,然后需要把红墨水提取出来,这恐怕不可能了。所以回声消除被认为是神秘和难以理解的技术也就不奇怪了。诚然,如果仅仅单独拿来一段混合了回音的语音信号,要去掉回音也是不可能的(就算是最先进的盲信号分离技术也做不到)。但是,实际上,除了这个混合信号,我们是可以得到产生回音的原始信号的,虽然不同于回音信号。
我们看下面的 AEC 声学回声消除框图(本图片转载)。
 
Figure   Acoustic Echo Cancellation in a voice communication terminal
 
其中,我们可以得到两个信号:一个是蓝色和红色混合的信号 1 ,也就是实际需要发送的 speech 和实际不需要的 echo 混合而成的语音流;另一个就是虚线的信号 2 ,也就是原始的引起回音的语音。那大家会说,哦,原来回声消除这么简单,直接从混合信号 1 里面把把这个虚线的 2 减掉不就行了?请注意,拿到的这个虚线信号 2 和回音 echo 是有差异的,直接相减会使语音面目全非。我们把混合信号 1 叫做近端信号 ne ,虚线信号 2 叫做远端参考信号 fe ,如果没有 fe 这个信号,回声消除就是不可能完成的任务,就像“巧妇难为无米之炊”。
虽然参考信号 fe echo 不完全一样,存在差异,但是二者是高度相关的,这也是 echo 称之为回音的原因。至少,回音的语义和参考信号是一样的,也还听得懂,但是如果你说一句,马上又听到自己的话回来一句,那是比较难受的。既然 fe echo 高度相关, echo 又是 fe 引起的,我们可以把 echo 表示为 fe 的数学函数: echo=F fe )。函数 F 被称之为回音路径。在声学回声消除里面,函数 F 表示声音在墙壁,天花板等表面多次反射的物理过程;在线路回声消除里面,函数 F 表示电子线路的二四线匹配耦合过程。很显然,我们下面要做的工作就是求解函数 F 。得到函数 F 就可以从 fe 计算得到 echo ,然后从混合信号 1 里面减掉 echo 就实现了回声消除。
 
尽管回声消除是非常复杂的技术,但我们可以简单的描述这种处理方法:
   1 、房间 A 的音频会议系统接收到房间 B 中的声音
   2 、声音被采样,这一采样被称为回声消除参考
   3 、随后声音被送到房间 A 的音箱和声学回声消除器中
   4 、房间 B 的声音和房间 A 的声音一起被房间 A 的话筒拾取
   5 、声音被送到声学回声消除器中,与原始的采样进行比较,移除房间 B 的声音
 
求解回音路径函数 F 的过程恐怕就是比较难以表达的数学公式了。鉴于通俗表达数学公式的难度比发现数学公式还难,笔者就不费力解释了。 下面这段表达了利用自适应滤波器原理求解 函数 F 的过程。(以下可以跳过
 
自适应滤波器
自适应滤波器是以输入和输出信号的统计特性的估计为依据,采取特定算法自动地调整滤波器系数,使其达到最佳滤波特性的一种算法或装置。自适应滤波器可以是连续域的或是离散域的。离散域自适应滤波器由一组抽头延迟线、可变加权系数和自动调整系数的机构组成。附图表示一个离散域自适应滤波器用于模拟未知离散系统的信号流图。自适应滤波器对输入信号序列 x (n) 的每一个样值,按特定的算法,更新、调整加权系数,使输出信号序列 y (n) 与期望输出信号序列 d (n) 相比较的均方误差为最小,即输出信号序列 y (n) 逼近期望信号序列 d (n)

 
以最小均方误差为准则设计的自适应滤波器的系数可以由维纳 - 霍甫夫方程解得。
B. 维德罗提出的一种方法,能实时求解自适应滤波器系数,其结果接近维纳-霍甫夫方程近似解。这种算法称为最小均方算法或简称  LMS 法。这一算法利用最陡下降法,由均方误差的梯度估计从现时刻滤波器系数向量迭代计算下一个时刻的系数向量
   
式中 ks 为一负数,它的取值决定算法的收敛性,  V ε 2 (n) 】为均方误差梯度估计,
自适应滤波器应用于通信领域的自动均衡、回声消除、天线阵波束形成,以及其他有关领域信号处理的参数识别、噪声消除、谱估计等方面。对于不同的应用,只是所加输入信号和期望信号不同,基本原理则是相同的。(以上部分可以跳过
 
上面这段话表明,需要求解的回音路径函数 F 就是一个自适应滤波器 W (n) 收敛的过程。所加输入信号 x (n) fe ,期望信号是 echo ,自适应滤波器收敛后的 W (n) 就是回音路径函数 F   收敛之后,当实际回音发生,我们把 fe 通过函数 W (n) ,就可以得到一个很准确的 echo ,把混合信号直接减去 echo ,得到实际需要发送的语音 speech ,完成回声消除任务。
值得注意的两点:
1、            自适应滤波器收敛阶段,期望信号是完全的 echo ,不能混杂有 speech 。因为 speech fe 是没有关系的,会扰乱 W (n) 的收敛过程。也就是说要求回声消除算法开始运转后收敛要非常快,最好对方还来不及说话,你一说就收敛好了;收敛好之后,如果对方开始说话,也就是有 speech 混合过来,这个 W (n) 系数就不要变化了,需要稳定下来。
2、             回音路径可能是变化的,一旦出现变化,回声消除算法要能判断出来,因为自适应滤波器学习要重新开始,也就是 W (n) 需要一个新的收敛过程,以逼近新的回音路径函数 F
基本上来说,上面这两点是两难的,一个需要自适应滤波器收敛后保持系数稳定,以保证不受 speech 说话干扰,另一个需要自适应滤波器随时保持更新状态,以保证能够追踪变化的 回音路径。这样一来,仅从数学算法层面,回声消除已经是难上加难!简单地说,回声消除自适应滤波器的设计具有两个互相矛盾的特性,也就是快速收敛和高度的稳定性,如何同时实现这两项特性,正是设计上的主要挑战。
经过上面的分析,相信大家对回声消除的原理和技术有了深刻的理解,这是一门即容易理解又难以实现的技术。
目录
相关文章
信道建模流程 | 带你读《大规模天线波束赋形技术原理与设计 》之二十八
本节将详细介绍衰落信道的整体建模流程,内容上与 3D 信道模 型 3GPP TR36.873 7.3 节和 3GPP TR38.901 的 7.5 节对应。两者在内容上大体相同,前者的目标为6GHz以下的信道建模(记为模型1),后者为0.5~100GHz 的信道建模(记为模型 2)。对于 6GHz 以下的信道建模,两者均可以使用, 在下文的描述中,两者不同的地方均会列出。
信道建模流程  | 带你读《大规模天线波束赋形技术原理与设计 》之二十八
|
机器学习/深度学习
小尺度信道建模 | 带你读《大规模天线波束赋形技术原理与设计 》之二十六
小尺度衰落是指无线电信号在短时间或短距离(若干波长)传播后其幅度、 相位或多径时延的快速变化。这种衰落是由于同一传输信号沿不同的路径传播, 由不同时刻(或相位)到达接收机的信号互相叠加所引起的,这些不同路径到 达的信号称为多径信号,接收机的信号强度取决于多径信号的强度、相对到达 时延以及传输信号的带宽。
 小尺度信道建模  | 带你读《大规模天线波束赋形技术原理与设计 》之二十六
大尺度信道建模 | 带你读《大规模天线波束赋形技术原理与设计 》之二十四
本节将介绍 3D 信道模型中 不同传输场景的大尺度衰落模型,包括路损计算、穿透损耗、直射径概率、阴 影衰落等。
大尺度信道建模  | 带你读《大规模天线波束赋形技术原理与设计 》之二十四
|
机器学习/深度学习 传感器 算法
毫米波系统中混合预编码的交替最小化附matlab代码
毫米波系统中混合预编码的交替最小化附matlab代码
一种改善调制宽带变换器的有意混叠方法(Matlab代码实现)
一种改善调制宽带变换器的有意混叠方法(Matlab代码实现)
|
算法
MATLAB|主动噪声和振动控制算法——对较大的次级路径变化具有鲁棒性
MATLAB|主动噪声和振动控制算法——对较大的次级路径变化具有鲁棒性
107 0
|
机器学习/深度学习 传感器 算法
【通信-SCMA】基于串行干扰消除的多用户检测算法附matlab代码
【通信-SCMA】基于串行干扰消除的多用户检测算法附matlab代码
|
算法 量子技术
从实际代码出发了解量子相位估计算法的原理
从实际代码出发了解量子相位估计算法的原理
309 0
|
机器学习/深度学习 算法 数据挖掘
自适应共振理论网络-1|学习笔记
快速学习自适应共振理论网络-1
自适应共振理论网络-1|学习笔记
|
机器学习/深度学习 存储 算法
自适应共振理论网络-2 | 学习笔记
快速学习自适应共振理论网络-2 。
自适应共振理论网络-2 | 学习笔记