技术纖 本发明涉及回声消除技术领域， 更为具体地， 涉及一种抑制残留回 声的方法和装置。

发明背景 语音通讯中，受话器（又称听筒、 receiver、SPK：、 EAR和 EARPHONE) 的信号经过线路反射和声学反射， 会串入送话器 （又称麦克风、 咪头、 微音器、 拾音器） 接收信号中， 并馈送到远端， 使远端听到回声。 回声 会给通话双方造成极大干扰， 影响通话质量， 严重时甚至引起啸叫， 不 仅完全无法通话， 还可能伤害通讯设备。 因此， 为了保证通话质量和设 备安全， 通常需要在语音通讯中对回声进行抑制。

回声信号是受话器信号经由电声转换为实际声 音信号播放出来、 再 经过环境反射产生的。 受话器的电声转换过程和环境反射都可以视为 滤 波过程， 因此回声信号可以看作是受话器经过一个特定 的滤波器产生的 声音信号。

简单的回声抑制是使信道变成半双工， 任何时刻信道内只有单向的 信号在传送， 抑制回声。 但在远近端双讲时， 远端在讲话时听不到近端 的声音， 通话的流畅性受到干扰。

随着对语音通讯要求的不断提升， 一般会采用自适应回声消除技术 来抑制回音， 同时保护近端语音， 保证双工效果。

自适应回声消除技术以自适应滤波器作为回声 路径消除滤波器来 进行回声消除， 滤波器通过比较受话器信号和送话器信号， 自动实时跟 踪回声反射环境的变化， 得到准确的回声路径， 消除回声。 常规的可以 去除大部分回声， 不会影响双工性能。 但是在实际系统中， 由于受话器 本身有失真， 常规的自适应回声消除技术无法消除全部回声 ， 而会有部 分的回声残留。 现有的抑制或消除残留回声技术主要有两种方 式， 一种是监测残留 回声强度， 限制受话器的功率， 使残留回声低于限定水平， 但该种方法 会导致受话器信号强度有较大起伏， 影响近端使用者的听感。 另一种方 法是在残留回声较强时， 将信道转为半双工， 显而易见， 这种方法在抑 制残留回声时会对近端语音造成损伤， 特别是对于一些便携的小型免提 通话设备， 由于受话器的失真较大， 而且由于结构紧凑， 送话器和受话 器距离很近， 残留回声很强时， 为抑制残留回声会使信道变为半双工， 导致通话的流畅性受到影响。 发明内容 鉴于上述问题， 本发明结合阵列处理与回声抑制， 充分利用小型免 提通话设备的声学结构， 提出一种抑制残留回声的方法和装置， 用以在 降低残留回声时减轻对近端语音的损伤， 提升双工性能。

本发明公开了一种抑制残留回声的方法， 该方法适用于包括 M个送 话器和 1个受话器的通话设备， 其中， M为大于 1的自然数， 所述 M个送 话器一字排开形成阵列， 该方法包括：

将 M路送话器信号分别与受话器信号进行自适应� �波处理得到 M路 自适应滤波信号；

将除第一路以外的其他 M-1路自适应滤波信号经过各自对应的阵列 滤波器处理后得到 M-1路阵列滤波器输出信号； 其中， 根据受话器和自 适应滤波信号对应的送话器之间的相对位置确 定该自适应滤波信号所 对应的阵列滤波器的权值；

将第一路自适应滤波信号分别与所述 M-1路阵列滤波器输出信号进 行信号相减处理得到 M-1路差值信号;将该 M-1路差值信号分别进行时频 转换处理， 并从中选择能量最小的一路频域信号；

将第一路和第 M路自适应滤波信号进行时频转换处理后进行� �音概 率滤波处理得到一路频域语音概率信号；

将所述频域语音概率信号和所选择的能量最小 的一路频域信号进 行相乘处理， 并将相乘处理结果进行频时转换处理后作为抑 制残留回声 后的语音输出信号。

在上述方法中， 在将所述频域语音概率信号和所选择的能量最 小的 一路频域信号进行相乘处理之后， 该方法进一步包括：

将相乘处理结果与时频转换后的受话器信号进 行频谱滤波处理， 再 将相乘处理结果与频谱滤波后的结果进行相减 处理， 将相减处理后的结 果进行频时转换处理后作为抑制残留回声后的 语音输出信号。

在上述方法中， 所述根据受话器与自适应滤波信号对应的送话 器以 及第一路送话器之间的相对位置确定该自适应 滤波信号所对应的阵列 滤波器的权值包括：

对于除第一路以外的其他 M-1路自适应滤波信号中的每一个当前自 适应滤波信号， 根据如下公式确定阵列滤波器的权值：

其中， 2是阵列滤波器， 是第一路自适应滤波信号， 是 所述的当前自适应滤波信号， 1¾是受话器和第一路送话器之间 的距离， 1) ₂ 是受话器和当前自适应滤波信号对应的送 话器之间 的距离， C为声速。

在上述方法中， 所述将第一路和第 M路自适应滤波信号进 行时频转换处理后进行语音概率滤波处理得到 一路频域语音概 率信号包括：

根据第一路和第 M路自适应滤波信号进行时频转换处理后 的两路结果计算声信号在每个频点上的到达送 话阵器列的到达 角度 «( )；

然后根据如下公式计算语音概率 pFi :

PF(f) = \ - a(f) / 9 ；

c/o为所述得到的频域语音概率信号。

在上述方法中， 所述将相乘处理结果与时频转换后的受话 器信号进行频谱滤波处理包括：

将全频率划分为 N个子带， 子带边界为 〜 ₊₁ ， 在每个子 带内进行如下计算：

计算匹配函数 f _M C :

其中， 为子带编号， ₂ )为所述相乘处理结果， 为 受话器信号经过时频转换后的结果；

计算匹配回声 (/)： Y _M {f) = X{f)H _M {F) ;

将匹配回声！ ^ )乘以一定的因子 得到频谱滤波后的结 则所述再将相乘处理结果与频谱滤波后的结果 进行相减处

本发明还公开了一种抑制残留回声的装置， 该装置适用于包括 M个 送话器和 1个受话器的通话设备， 其中， M为大于 1的自然数， 所述 M个 送话器一字排开形成阵列， 该装置包括： M个自适应滤波组件， M-1个 阵列滤波组件， 1个比较选择器， 1个语音概率估计组件， 1个乘法器和 1 个时频变换器； 其中：

M个自适应滤波组件， 用于将 M路送话器信号分别与受话器信号进 行自适应滤波处理， 输出 M路自适应滤波信号；

M-1个阵列滤波组件分别与除第一路以外的其他 M-1路自适应滤波 信号一一对应， 且每个阵列滤波组件包括： 1个阵列滤波器、 1个减法器 和 1个时频变换器；

其中： 阵列滤波器， 用于对所对应的自适应滤波信号进行阵列滤波 处理后输出给减法器； 减法器， 用于将第一路自适应滤波信号与阵列滤 波器输出的信号进行相减处理后输出给时频变 换器； 时频变换器， 用于 对接收的信号进行时频变换处理后输出频域信 号给比较选择器； 比较选择器， 用于接收 M-l个阵列滤波组件输出的 M-l个频域信号， 从中选择能量最小的一路频域信号输出给乘法 器；

语音概率估计组件包括： 2个时频变换器和 1个语音概率估计器； 其 中：所述 2个时频变换器，用于分别对第一路和第 M路自适应滤波信号进 行时频转换处理后输出给语音概率估计器； 语音概率估计器， 用于根据 所接收的 2路信号进行语音概率滤波处理输出一路频域� �音概率信号给 乘法器；

乘法器， 用于对所接收的 2路频域信号进行相乘处理后输出给频时 变换器；

频时变换器， 用于对所接收信号进行频时转换处理后作为抑 制残留 回声后的语音输出信号。

该装置在所述乘法器和频时变换器之间进一步 包括： 1个频谱滤波 组件；

频谱滤波组件包括： 1个时频变换器， 1个频谱滤波器和 1个减法器； 其中： 时频变换器， 用于对受话器信号进行时频转换处理后输出给 频谱 滤波器； 频谱滤波器， 用于对乘法器的输出信号和时频变换器的输出 信 号进行频谱滤波处理后输出给减法器； 减法器， 用于将乘法器的输出信 号和频谱滤波器的输出信号进行减法处理后输 出给频时变换器。

在该装置中， 所述阵列滤波器的权值由如下公式确定：

其中， 2是阵列滤波器， 是第一路自适应滤波信号， 是 该阵列滤波器对应的自适应滤波信号， A是受话器和第一路送 话器之间的距离， ₂ 是受话器和该阵列滤波器对应的送话器之 间的距离， C为声速。

在该装置中， 所述语音概率估计器， 用于根据所接收的 2 路信号计算声信号在每个频点上的到达送话器 阵列的到达角度

«(/);并根据如下公式计算语音概率 ^ (/): pF(f) = l-a(f)/9 ； 将^ σ M乍为频域语音概率信号输出给减法器。

在该装置中， 所述频谱滤波器， 用于将全频率划分为 Ν个 子带， 子带边界为 Si〜 ₊₁ ， 在每个子带内进行如下计算：

计算匹配函数 f _M C :

其中， 为子带编号， ₂ )为乘法器的输出信号， X 为 时频变换器的输出信号；

计算匹配回声 (/)： Y _M {f) = X{f)H _M {F);

将匹配回声！ ^ )乘以一定的因子 得到频谱滤波后的结 果后输出给减法器。 根据本发明的抑制残留回声的方法和装置， 借助新的回声消除技术 将残留回声滤波设计与阵列空间滤波相结合， 充分利用小型免提通话设 备的声学特性以及送话器和受话器的位置信息 ， 利用回声和近端语音到 送话器的时间和相位差异， 提取语音概率信息， 将语音区域和回声区域 划分开来， 有效滤除残留回声， 并保护近端语音。 为了实现上述以及相关目的， 本发明的一个或多个方面包括后面将 详细说明并在权利要求中特别指出的特征。 下面的说明以及附图详细说 明了本发明的某些示例性方面。 然而， 这些方面指示的仅仅是可使用本 发明的原理的各种方式中的一些方式。 此外， 本发明旨在包括所有这些 方面以及它们的等同物。 附图简要说明 通过参考以下结合附图的说明及权利要求书的 内容， 并且随着对本 发明的更全面理解， 本发明的其它目的及结果将更加明白及易于理 解。 在附图中：

图 1 为本发明所应用的受话器和送话器阵列相对位 置实施例示意 图；

图 2 为本发明所应用的小型免提通话设备使用场景 的实施例示意 图；

图 3为本发明实施例的送话器阵列回声消除系统� �方框示意图； 图 4为送话器信号与近端语音成分的能量对比示� �图；

图 5为根据本发明实施例的阵列滤波器的形状示� �图；

图 6为本发明实施例的经阵列滤波器滤波后的效� �示意图； 图 7为根据本实施例的语音概率抑制后的信号能� �对比示意图； 图 8为受话器信号和残留回声频谱的对比示意图� �

图 9为匹配回声和残留回声频谱的对比示意图；

图 10为本发明进行回声滤波后的最终效果比较示� ��图；

图 11 为根据本发明实施例的送话器阵列回声消除方 法的流程示意 图；

图 12为本发明实施例的残留回声滤波的流程示意� ��；

图 13 是本发明实施例中的抑制残留回声的装置的通 用结构的示意 图；

图 14是本发明实施例中的一种抑制残留回声的方� ��通用流程图。 在所有附图中相同的标号指示相似或相应的特 征或功能。

实施本发明的方式 以下将结合附图对本发明的具体实施例进行详 细描述。

本发明采用送话器阵列消除残余回声， 因此适用于有多个送话器和 单个受话器的小型免提通话设备。 图 1和图 2分别为本发明所应用的受 话器和送话器阵列相对位置以及小型免提通话 设备使用场景的实施例 示意图， 如图 1和图 2所示， 在这种通话设备中， 一般受话器放置在送 话器阵列 90度方向， 使用者在送话器阵列 0度方向。 这与某些应用场 景一致， 比如车载免提 caridt, 使用者正对送话器阵列， 位置较为固定。 图 1和图 2演示的是双送话器的例子， 两个以上多送话器阵列的摆放与 之类似。

对于使用者来说， 使用者到送话器的距离接近相等， 也就是说， 送 话器阵列中的每一个送话器所接收到的使用者 发出的语音信号是基本 相同的。 但是对于受话器而言， 受话器到各个送话器的距离不等， 以双 送话器阵列为例， 如果受话器和送话器 1的距离是 D 与送话器 2的距 离是 D ₂ ， 并且 0^>=20 ₂ ， 例如 D^ m, D ₂ =4cm。 由于受话器发出的回 声， 到两个送话器之间有相位差， 而使用者发出的语音到两个送话器是 同相的， 因此根据相位关系的差异， 可以区分语音和回声， 达到分离语 音和回声的目的。 本发明正是利用这一相位关系的差异来进行语 音和回 声的分离。

图 3为本发明实施例的送话器阵列回声消除系统� �方框示意图。 如图 3所示， 本发明所提供的送话器阵列回声消除系统 300主要由 自适应回声滤波单元 320和残留回声滤波单元 340两部分组成， 其中自 适应回声滤波单元 320和残留回声滤波单元 340在结构上为级联关系， 自适应回声滤波单元 320的输入是受话器信号和送话器阵列的两路送 话 器信号， 残留回声滤波单元 340的输入信号是受话器信号和自适应回声 滤波单元 320的两路输出， 残留回声滤波单元 340的输出信号是送话器 阵列回声消除系统的输出， 即分离了回声的语音信号。

从送话器阵列回声消除系统 300的整体应用来说， 本发明所提出的 送话器阵列回声消除系统连接在送话器和受话 器之间， 受话器信号 X和 送话器信号 d (包括送话器阵列的所有信号 d d ₂ . . .d _M ) 是系统输入， 送话器信号 d由回声信号 y以及近端语音信号 V组成。 两路送话器信号 经过自适应回声滤波单元 320， 回声成分被大部分滤除， 仅有部分的残 留回声进入残留回声滤波单元 340， 经残留回声滤波单元 340后， 残留 回声也被滤除， 最后得到的是只有近端语音信号 V， 作为分离了回声的 语音信号被送话器阵列回声消除系统 300输出送到语音通讯的远端。

图 4为送话器信号与近端语音成分的能量对比示� �图。如图 4所示， 在典型应用中， 实线表示送话器信号 d 虚线表示送话器信号中的近端 语音成分 _Vl ， 点线表示回声成分 _yi 。 可以看到， 送话器信号中， 近端语 音成分的能量很低， 完全被淹没在回声中。

下面分别对自适应回声滤波单元 320和残留回声滤波单元 340的组 成结构和信号处理过程进行示例性说明。

自适应回声滤波单元 320接收受话器信号和送话器阵列信号， 基于 所接收的受话器信号分别对送话器阵列中的每 一个送话器信号进行回 声滤波， 以得到回声滤波信号阵列。 在图 3所示的实施例中， 自适应回 声滤波单元 320的输入有三路， 受话器信号 X和两路送话器信号 d d ₂ ， 输出为自适应滤波后的回声滤波信号阵列， 在图 3所示的实施例中为 _ei 与 e ₂ 。 在本发明中， 自适应回声滤波单元 320部分的工作原理与通用的 自适应回声滤波相似， 可以采用时域或者频域或者时频域混合的滤波 形 式。

在本实施例中， 自适应回声滤波单元 320包括与送话器阵列中的每 一送话器分别对应设置的滤波器 321、 滤波控制器 322和加法器 323构 成。 其中， 滤波器和滤波控制器分别为自适应滤波器和自 适应滤波控制 器， 通过比较受话器信号和送话器信号的相似性， 自适应匹配出回声信 号并从借助加法器 323送话器信号中滤除回声信号， 得到相应的回声滤 波信号。

残留回声滤波单元 340与自适应回声滤波单元 320级联， 基于所接 收的受话器信号， 对自适应回声滤波单元 320输出的回声滤波信号阵列 进行残留回声滤波。 也就是说， 残留回声滤波单元 340的输入有三路， 受话器信号 X和自适应回声滤波单元 320输出的回声滤波信号阵列 _ei 和 e ₂ ， 输出为去除所有回声成分的语音信号 e。 _M ,。本实施例中的残留回声滤 波单元 340主要包括阵列滤波器 341、第一加法器 342、时频变换器 343、 语音概率估计器 344和乘法器 345， 另外还可以包括频谱滤波器 346、 第二加法器 347以及频时变换器 348。

在残留回声滤波单元 340 收到受话器信号和自适应回声滤波单元 320输出的回声滤波信号阵列后， 时频变换器 343会将回声滤波信号阵 列中的每一路回声滤波信号 _ei 和 分别从时域变换到频率域，以便进行 频率域的语音概率估计。 时频变换可以采用傅立叶变换实现， 也可以用 改进离散数字余弦变换等方式实现。

其中， 阵列滤波器 341与回声滤波信号阵列中的 相卷积， 以利用 卷积结果初步消除 _ei 中的残留回声成分。

由于受话器和送话器阵列之间的位置差异， 送话器 1中的回声， 近 似于送话器 2中的回声与一个阵列滤波器相卷积， 因此， 将这样一个阵 列滤波器与送话器 2的回声滤波信号 卷积， 再从送话器 1中减去， 就 会减小送话器 1中的残留回声成分， 获得第一残留回声滤波信号 em 该第一残留回声滤波信号 e _mi 在每一路回声滤波信号 _ei 和 e ₂ 的基础上初 步滤除了部分残留回声。下面对阵列滤波器 341的权值做如下分析推导: 在如图 1所示的通话设备中， 送话器和受话器的位置是确定的， 根 据声传播的物理特性， 回声到两个送话器之间的相对关系也是确定的 ， 若受话器发声时， 受话器中心的声音为 则传播到两个送话器上的信 号 近似

其中， 0和0 ₂ 分别为送话器到受话器 1和受话器 2之间的距离， t 表示当前时间， c为声速。 为能量衰减因子， 由受话器电声特性决定， 为定值。

两路残留回声的关系也与之类似， 符合近似关系：

因此可以由受话器与送话器之间的相对位置确 定阵列滤波器的权 值， 即:幅度为 D ₂ /D 时延为 (D _r D ₂ )/c。

图 5为根据本发明实施例的阵列滤波器的形状示� �图。如图 5所示， 以受话器距离两个送话器的距离 Df n^ D ₂ =4cm为例， 滤波器的峰值 位置， 其时延是 (D _r D ₂ )/c， 其幅度近似为 D ₂ /D _{l t3} 此阵列滤波器的权值可 以根据具体的应用需求事先离线算出并固定。

图 6为本发明实施例的经阵列滤波器滤波后的效� �示意图。 如图 6 所示， 点线部分表示自适应回声滤波后的信号 E 实线部分表示阵列滤 波后的信号 E _mi ， 虚线部分表示近端语音成分 ¼， 从图 6的三部分信号 的能量对比可以看出， 经阵列滤波器 341阵列滤波前后残留回声能量变 化， 平均降低约 6dB。 通过效果对比示意图能够充分说明利用阵列滤 波 器 341能够在一定程度上滤除部分回声信号。

对于送话器多于两个的送话器阵列， 阵列滤波器 341的计算方式与 该具有两个送话器的送话器阵列的计算方式类 似。设有 M个送话器， 自 适应回声滤波的输出为 E E ₂ 〜E _M ， 通过送话器 1和其它送话器 k之间的 相对位置，可以计算送话器 1和其他送话器 k之间的阵列滤波器 /¾(l<k ^M), 这样通过 (M-l)个阵列滤波器， 就可以得到 (M-1)个阵列滤波器输 出 E _mi — _k ， 取其中能量最小的输出即可作为最终的 E _mi 。

阵列滤波器的计算需满足最小均方差准则， 设此阵列滤波器为 h， 则 h需要满足

dE

= 0

dh

其中 E[.]为取平均值运算。 如果设阵列滤波器输出为 em

em、 = e、一 h * e 和

L

em _l (n) = e _l (n) - ^e ₂ (n

k=l

其中 L为阵列滤波器长】 通过上述阵列滤波器和加法器的回声滤波处理 ， 所得到的第一残留 回声滤波信号中仍然残留有一定的回声信号， 因此， 在获得的第一残留 回声滤波信号 的基础上， 本发明对其中的残留回声进行进一步的滤 波处理， 在这一层的回声滤波处理中， 主要采用语音概率估计器和乘法 器相结合进行。

首先，在经阵列滤波器 341通过卷积滤除 _ei 中的部分回声信号输出 611 ₁ 后， 需要利用时频变换器 343将该信号 em A时域变换到频率域， 即从 e _mi 变换到 E _mi 。

语音概率估计器 344通过比较两路回声滤波信号 _ei 和 经时频域变 换后的两路信号 和 E ₂ 的时间和相位关系， 得出指示两路回声滤波信 号中的语音和回声分布在哪些频率区域的频域 语音概率信息 pF。

语音概率估计的工作原理是， 两路送话器信号 E ₂ 中都含有残留 回声和语音信号， 残留回声是从 90度方向来的， 存在相位差， 语音信 号则是从 0度方向来的， 同相位。 因此， 当近端语音越强， 语音的成分 越多， £ ₁ 和 E ₂ 的相位就越接近。 当近端语音很弱， 语音成分越少， 和£ ₂ 的相位差就越明显。 通过比较各个频点上 Ε^Π Ε^相位， 可以获 得送话器阵列获得的语音信号在频率上的分布 。

具体实现方法是， 根据 E ₁ 和 ¾计算声信号在空间上的到达角度， 也即声信号是从哪个方向来， 如果来自 90度方向， 则信号是回声信号， 语音概率为 0， 如果来自 0度， 则是语音方向， 语音概率为 1。 如果来 自于 0和 90度之间的方向 α，则语音概率在 0和 1之间，具体值为 1-α/90。

经过语音概率判断后， 可以计算得到每个频点上的语音概率： PF(f) = l - a(f) / 9 其中 /为频率。 多于两个的送话器阵列，可以选取第 1个和第 M个送话器的自适应 滤波器输出 和 E _M 计算语音概率， 方法与双送话器阵列相同。

将阵列滤波后的信号 e _mi 经时频转换后变为 Em 与语音概率 pF 相乘， 则可以得到输出 Em ₂ 。 Em ₂ (f) = Em _l (f) * pF(f)

图 7为根据本实施例的语音概率抑制后的信号能� �对比示意图。 如 图 7所示， 点线部分表示阵列滤波后信号 E _mi ， 实线部分表示语音概率 抑制后信号 Em ₂ ， 虚线部分表示近端语音成分 ¼， 从图 7的三部分信号 的能量对比可以看出， 经过阵列滤波和语音概率相乘后， 残留回声信号 进一步降低， 而语音成分则被保留。 1¾的语音-回声比可以比 提高 10dB以上， 从而进一步抑制残留回声， 得到更纯净的近端语音。

另外， Em ₂ 中的残余回声还可以通过频谱滤波器 346进一步去除。 残余回声从能量上看， 有很多能量峰值， 将峰值抑制下去， 就可以进一 步去除回声。

图 8为受话器信号和残留回声频谱的对比示意图� � 如图 8所示， 实 现部分表示残留回声， 点线部分表示受话器信号， 由于回声是由受话器 信号产生的， 回声和受话器信号的谐波峰值位置相同或相近 ， 只是总体 起伏形状和信号能量不同。 因此， 可以将受话器信号和残余回声信号做 频谱包络形状匹配， 再乘以一定的因子 Ag做能量匹配， 并与残余回声 信号做相位匹配， 得到匹配回声， 然后再将匹配回声从残留回声信号中 减去， 可以去除残留回声。 Ag是不小于 1 的因子， 一般取 1〜8之间， 根据残留回声能量而定， 为了能将残留回声滤除干净， 残留回声越强， Ag 则应越大， 频谱滤波的力度越强， 但对近端语音伤害也就越大。 而 经过语音概率估计的处理后， 残留回声已经较弱， 因此匹配回声可以取 得较小， 不太强的频谱滤波可以将残留回声滤除， 相应的， 近端语音则 可以较好保留。

具体地， 频谱滤波器 346的频谱滤波过程如下：

先将受话器信号和残余回声信号做频谱包络形 状匹配， 得到回声估 计。 可以用以下方式进行：

将全频率划分为 M个子带， 子带边界为 ^〜¾^， M可以为 32或 16。 Em2和 X在每个子带内求能量， 将能量相除， 得到能量匹配函数 H _M 。 受话器信号 X乘以 H _M 得到匹配回声 Y _M 。

图 9为匹配回声和残留回声频谱的对比示意图。 如图 9所示， 实现 部分表示残留回声， 点线部分表示匹配回声， 匹配效果比较明显。

计算匹配函数的方法如下面的公式，对于属于 第 M个子带的频 其匹配函数为 f-B _M

∑ Em ₂ ² (f)

H _M 、f、二^^ ， _l < ^ f j < B

∑ x f)

f=Bi

而经过频谱包络形状匹配的回声为

Y _M (f) = ^(fW _M (F)

再经过能量匹配和相位匹配， 得到的回声估计信号为

最后利用加法器从 Em ₂ 中减去回声估计信号， 即得到最终的语音 $1刖

在最终消除回声对语音的影响之后， 尚需利用频时变换器将语音信 号从频率域变回到时域， 频时变换可以用逆傅立叶变换实现， 也可以用 逆离散数字余弦变换等其他方式实现。

经过频时转换后， 频域信号 变换为时域信号 ^， 作为送话器 阵列回声消除系统 300的总输出。

图 10为本发明进行回声滤波后的最终效果比较示� ��图。 如图 10所 示， 实线表示送话器信号 dl， 虚线表示送话器信号中的近端语音信号成 分 vl， 点线表示系统输出 eout。 从图 10中的信号能量对比可以看出， 在只有回声的区域，输出信号 eout能量很小，说明回声被压制的很彻底， 在有近端语音的部分， 系统输出和语音成分能量接近， 说明近端语音被 很好保护。

图 11 为根据本发明实施例的送话器阵列回声消除方 法的流程图。 如图 11所示， 该方法始于步骤 S1110, 在送话器阵列接收到包括近端语 音信号 v和受话器信号 X后， 首先基于受话器信号对送话器阵列信号进 行自适应回声滤波， 得到回声滤波信号阵列。 具体地说， 对于具有两个 送话器的送话器阵列而言， 送话器 1接收到的声音信号 dl包括近端语 音信号 vl和受话器产生的回声信号 yl， 即 dl=yl+vl ; 送话器 2接收到 的声音信号 d2包括近端语音信号 v2和受话器产生的回声信号 y2， 即 d2=y2+v2。 分别对这两个送话器所接收的声音信号进行自 适应滤波， 通 常采用自适应滤波器、 自适应滤波控制器和加法器构成的自适应滤波 单 元实现。 自适应滤波后得到回声滤波信号阵列， 在本实施例中即为 el 和 e2。

然后， 在步骤 S1120中， 基于所接收的受话器信号， 对所得到的回 声滤波信号阵列进行残留回声滤波。 残留回声的滤波采用与自适应滤波 单元级联的残留回声滤波单元实现， 其中具体的滤波过程如图 12所示。

图 12示出了本发明实施例的残留回声滤波的流程� ��，如图 12所示， 经过初步的自适应滤波之后， 送话器所接收的信号中还会有部分回声残 留， 此时， 首先从回声滤波信号阵列中的一个回声滤波信 号中减去另一 个回声滤波信号与阵列滤波器相卷积所产生的 信号， 以输出第一残留回 声滤波信号 （步骤 S121 ) , 其中阵列滤波器的权值由受话器与送话器阵 列之间的相对位置确定； 然后将回声滤波信号阵列中的每一回声滤波信 号以及第一残留回声滤波信号分别从时域变换 到频率域 （步骤 S122) , 其中， 将回声滤波信号阵列中的每一回声滤波信号分 别从时域变换到频 率域的处理与步骤 S121 中滤波处理的顺序没有必然的先后， 可以同时 进行， 也可以顺序进行。

在将回声滤波信号阵列中的每一回声滤波信号 分别从时域变换到频 率域之后， 在步骤 S123 中， 通过比较经过时频域变换器变换的回声滤 波信号阵列中的每一回声滤波信号的时间和相 位关系， 确定语音和回声 分布区域的频域语音概率信息。

然后， 在步骤 S124 中， 根据所确定的频域语音概率信息进一步降 低所述从时域变换到频率域的第一残留回声滤 波信号中的残留回声信 号， 得到第二残留回声滤波信号。 这一滤波过程通过使第一残留回声滤 波信号与所确定的频域语音概率相乘实现。

接着， 在步骤 S125 中， 基于所接收的受话器信号和第二残留回声 滤波信号利用频谱滤波器确定回声估计信号。 具体地说， 将受话器信号 和残余回声信号做频谱包络形状匹配， 然后乘以一因子 Ag进行能量匹 配。 并与第二残留回声滤波信号做相位匹配， 得到匹配回声， 然后再根 据所得到的匹配回声确定回声估计信号， Ag的取值为不小于 1的实数， 一般取为 1〜8， 根据残留回声的强度而定。

在步骤 S126 中， 从第二残留回声滤波信号中减去所确定的回声 估 计信号， 即得到分离的语音信号。

最后， 在步骤 S127中， 将分离的语音信号从频率域变换到时域。 在上述实施例中以送话器阵列包含两个送话器 为例进行了说明。 在 实际应用中送话器阵列中的送话器个数也可以 是更大的数值， 如 3、 4、 5等。 下面对本发明中的通用的抑制残留回声的装置 和方法进行说明。

图 13是本发明实施例中的抑制残留回声的装置的� ��用结构的示意 图。 该装置适用于包括 M个送话器和 1个受话器的通话设备， 其中， M为 大于 1的自然数， 所述 M个送话器一字排开形成阵列， 如图 13所示， 该装 置包括： M个自适应滤波组件， M-1个阵列滤波组件， 1个比较选择器， 1个语音概率估计组件， 1个乘法器和 1个时频变换器； 其中：

M个自适应滤波组件， 用于将 M路送话器信号分别与受话 器信号进行自适应滤波处理， 输出 M路自适应滤波信号；

M- 1个阵列滤波组件分别与除第一路以外的其他 M- 1路自 适应滤波信号一一对应， 且每个阵列滤波组件包括： 1个阵列滤 波器、 1个减法器和 1个时频变换器；

其中： 阵列滤波器， 用于对所对应的自适应滤波信号进行 阵列滤波处理后输出给减法器； 减法器， 用于将第一路自适应 滤波信号与阵列滤波器输出的信号进行相减处 理后输出给时频 变换器； 时频变换器， 用于对接收的信号进行时频变换处理后 输出频域信号给比较选择器； 比较选择器， 用于接收 M- l个阵列滤波组件输出的 M- l个频 域信号， 从中选择能量最小的一路频域信号输出给乘法 器； 语音概率估计组件包括： 2个时频变换器和 1个语音概率估 计器； 其中： 所述 2个时频变换器， 用于分别对第一路和第 M路 自适应滤波信号进行时频转换处理后输出给语 音概率估计器； 语音概率估计器，用于根据所接收的 2路信号进行语音概率滤波 处理输出一路频域语音概率信号给乘法器；

乘法器，用于对所接收的 2路频域信号进行相乘处理后输出 给频时变换器；

频时变换器， 用于对所接收信号进行频时转换处理后作为 抑制残留回声后的语音输出信号。 在图 13所示的装置中， 当 M等于 2时： 所述比较选择器， 用 于接收 1个阵列滤波组件输出的 1个频域信号， 将该频域信号输 出给乘法器； 或者， 图 3所示的装置不包含比较选择器， 阵列滤 波组件中的时频变换器直接将频域信号输出给 乘法器， 此时， 图 13所示的装置变形为 3所示的装置。 如图 13所示， 该装置在所述乘法器和频时变换器之间进一 步包括： 1个频谱滤波组件；

频谱滤波组件包括： 1个时频变换器， 1个频谱滤波器和 1 个减法器； 其中： 时频变换器， 用于对受话器信号进行时频转 换处理后输出给频谱滤波器； 频谱滤波器， 用于对乘法器的输 出信号和时频变换器的输出信号进行频谱滤波 处理后输出给减 法器； 减法器， 用于将乘法器的输出信号和频谱滤波器的输出 信号进行减法处理后输出给频时变换器。 在图 13所示的装置中， 所述阵列滤波器的权值由如下公式 确定：

其中， 2是阵列滤波器， 是第一路自适应滤波信号， 是 该阵列滤波器对应的自适应滤波信号， A是受话器和第一路送 话器之间的距离， ₂ 是受话器和该阵列滤波器对应的送话器之 间的距离， C为声速。 在图 13所示的装置中， 所述语音概率估计器， 用于根据所 接收的 2路信号计算声信号在每个频点上的到达送话� �阵列的 到达角度 a、f、； 并根据如下公式计算语音概率 pF、f、： pF(f) = 1 - a(f) / 9 _; 将 pF f)作为频域语音概率信号输出给减法

在图 13所示的装置中， 所述频谱滤波器， 用于将全频率划 分为 N个子带， 子带边界为 β^β^, 在每个子带内进行如下计 计算匹配函数 f _M C :

f=B _i+ l

∑ Em ₂ ² (f)

H ( f) = ^ f=Bi

其中， 为子带编号， /)为乘法器的输出信号， χη为 时频变换器的输出信号；

计算匹配回声 (/)： Y _M {f) = X{f)H _M {F);

将匹配回声！ ^ /l乘以一定的因子 得到频谱滤波后的结 果后输出给减法器。 图 14是本发明实施例中的一种抑制残留回声的方� ��通用流 程图。 该方法适用于包括 M个送话器和 1个受话器的通话设备， 其中， M为大于 1的自然数，所述 M个送话器一字排开形成阵列， 如图 14所示该方法包括：

1401， 将 M路送话器信号分别与受话器信号进行自适应� � 波处理得到 M路自适应滤波信号；

1402， 将除第一路以外的其他 M- 1路自适应滤波信号经过 各自对应的阵列滤波器处理后得到 M- 1路阵列滤波器输出信 号； 其中， 根据受话器和自适应滤波信号对应的送话器之 间的 相对位置确定该自适应滤波信号所对应的阵列 滤波器的权值；

1403， 将第一路自适应滤波信号分别与所述 M- 1路阵列滤 波器输出信号进行信号相减处理得到 M- 1路差值信号； 将该 M- 1 路差值信号分别进行时频转换处理， 并从中选择能量最小的一 路频域信号；

1404， 将第一路和第 M路自适应滤波信号进行时频转换处 理后进行语音概率滤波处理得到一路频域语音 概率信号；

1405， 将所述频域语音概率信号和所选择的能量最小 的一 路频域信号进行相乘处理， 并将相乘处理结果进行频时转换处 理后作为抑制残留回声后的语音输出信号。 在图 14所示的方法中， 当所述 M等于 2时：

所述将除第一路以外的其他 M- 1路自适应滤波信号经过各 自对应的阵列滤波器处理后得到 M - 1路阵列滤波器输出信号 为： 将除第一路以外的另一路自适应滤波信号经过 阵列滤波器 处理后得到 1路阵列滤波器输出信号；

所述将第一路自适应滤波信号分别与所述 M- 1路阵列滤波 器输出信号进行信号相减处理得到 M- 1路差值信号； 将该 M- 1 路差值信号分别进行时频转换处理， 并从中选择能量最小的一 路频域信号为： 将第一路自适应滤波信号与另一路阵列滤波器 输出信号进行信号相减处理得到 1路差值信号，将该路差值信号 进行时频转换处理；

所述将所述频域语音概率信号和所选择的能量 最小的一路 频域信号进行相乘处理为： 将所述频域语音概率信号和所述差 值信号进行时频转换处理后的信号进行相乘处 理。 在图 14所示的方法中， 在将所述频域语音概率信号和所选 择的能量最小的一路频域信号进行相乘处理之 后， 该方法进一 步包括：

将相乘处理结果与时频转换后的受话器信号进 行频谱滤波 处理， 再将相乘处理结果与频谱滤波后的结果进行相 减处理， 将相减处理后的结果进行频时转换处理后作为 抑制残留回声后 的语音输出信号。 在图 14所示的方法中， 所述根据受话器与自适应滤波信号 对应的送话器以及第一路送话器之间的相对位 置确定该自适应 滤波信号所对应的阵列滤波器的权值包括：

对于除第一路以外的其他 M- 1路自适应滤波信号中的每一 个当 ，根据如下公式确定阵列滤波器的权值：

其中， 2是阵列滤波器， 是第一路自适应滤波信号， 是 所述的当前自适应滤波信号， 1¾是受话器和第一路送话器之间 的距离， 1) ₂ 是受话器和当前自适应滤波信号对应的送 话器之间 的距离， C为声速。 在图 14所示的方法中， 所述将第一路和第 M路自适应滤波 信号进行时频转换处理后进行语音概率滤波处 理得到一路频域 语音概率信号包括：

根据第一路和第 M路自适应滤波信号进行时频转换处理后 的两路结果计算声信号在每个频点上的到达送 话阵器列的到达 角度 «( )；

然后根据如下公式计算语音概率 pF

PF(f) = \-a(f)/9 ；

c/o为所述得到的频域语音概率信号。 在图 14所示的方法中， 所述将相乘处理结果与时频转换后 的受话器信号进行频谱滤波处理包括：

将全频率划分为 N个子带， 子带边界为 〜 ₊₁ ， 在每个子 带内进行如下计算：

计算匹配函数 f _M C :

f=B _i+ l

2(f)

其中， 为子带编号， 为所述相乘处理结果， χη为 受话器信号经过时频转换后的结果；

计算匹配回声 (/)： Y _M {f) = X{f)H _M {F);

将匹配回声！ ^ )乘以一定的因子 得到频谱滤波后的结 果；

则所述再将相乘处理结果与频谱滤波后的结果 进行相减处 理为： Em ₂ ( | )|

Em,

如上参照附图以示例的方式描述了本发明用 于回声消除的送话器 阵列残留回声消除方法及系统。 显然， 本领域的普通技术人员在不偏离 本发明概念的条件下， 不需要付出创造性的劳动， 就可以对在此公开的 特定装置和技术进行许多应用和修改， 并且可以与在此所公开的特定装 置和技术有所不同。 因此， 本发明应该被理解成包括此处所公开的装置 和技术所提出或拥有的每一个新颖特征和新颖 的特征组合， 但凡本领域 普通技术人员根据本发明所揭示的内容所作的 等效修饰和变化， 皆应纳 入权利要求书中记载的保护范围内。