一种干扰信号与语音调制信号混叠的窄带干扰识别方法
【专利摘要】本发明涉及数字信号处理领域的语音信号干扰处理技术，特别涉及语音信号窄带干扰的识别方法。本发明首先对接收语音信号进行分块处理，然后对块内信号做FFT变换，利用语音信号与干扰信号频域响应幅度和宽度上的差异性对语音载频、干扰频率和干扰个数进行预识别，然后再利用前面的块记录中干扰源和干扰频点出现次数最多的组合来对窄带干扰进行最终识别，块间记录的统计信息能进一步的提高正确识别结果。本发明尤其适用于窄带干扰的识别。
【专利说明】一种干扰信号与语音调制信号混叠的窄带干扰识别方法
【技术领域】
[0001] 本发明属于数字信号处理领域的语音信号干扰处理技术，特别涉及语音信号窄带 干扰的识别方法。
【背景技术】
[0002] 无论是军用还是民用，人们都希望利用高科技手段改善语音通信的质量，提高其 清晰度和可懂度。在很多理论和应用中，所用的语音数据大部分都是在近似理想的条件下 采集的。大多数语音识别和语音编码在开始研究时都要在高保真设备上录制语音，尤其要 在无噪环境下录音。然而，当在实际应用时，由于噪声的存在会产生很多问题。比如，环境 噪声的污染使得许多语音处理系统的性能急剧恶化。语音识别正在步入实用阶段，但识别 系统大都是在安静环境中工作的，在噪声环境中尤其是强噪声环境，语音识别系统的识别 率将受到严重影响。而窄带干扰（NBI)普遍存在于各种语音通信系统中，特别是在军事应 用中，许多场合噪声干扰还相当严重，严重影响语音通信的质量。
[0003] 干扰抑制，就是尽可能消除干扰的能量；语音增强，就是从带噪语音中提取尽可 能纯净的语音。语音通信系统中的干扰抑制属于语音增强的范畴。语音增强目前已发展成 为语音信号处理的一个重要分支。语音增强技术在国外发展较早，近年来国内也做了大量 的工作。在近20年，语音增强算法的发展可以概括为以下几个方面：（1)频域研究以谱相 减为代表的幅度谱估计算法，围绕减小音乐噪声和残留噪声，减小失真等性能，发展和产 生了许多改进的谱相减算法；（2)时域研究围绕滤波器设计，研究诸如自适应陷波、最小均 方误差（LMS)等自适应滤波算法；(3)其他变换域或处理空间的算法，例如小波变换和奇 异矩阵分解等。
[0004] 上述语音增强技术的实现是基于干扰信号的识别技术，首先通过干扰信号识别技 术将语音信号与干扰信号区分，然后采用频域谱相减方法、时域滤波方法或其他变换域如 小波变换等方法去除干扰信号，提取尽可能纯净的语音信号。现有的干扰信号的识别技术 都是假设已知语音信号的接收频率，通过信号与干扰不同的频点来区分信号与干扰。而在 很多军用或民用场合，特别是如汽车、火车和飞机等高速移动的语音信号接收平台上，多普 勒频移严重，语音信号的接收频率相对于发射频率有了较大的频率偏移。由于无法准确预 估多普勒频移，因此不能通过语音信号的发射频率来推断语音信号的接收频点，也就无法 通过信号与干扰的频率差异性来区分信号与干扰。
【发明内容】
[0005] 本发明的目的，就是针对现有语音通信中干扰信号的识别技术需要预知语音信号 的接收频率，无法适应高速移动的语音信号接收平台的缺陷，提出一种干扰信号与语音调 制信号混叠的窄带干扰识别方法。
[0006] 本发明的技术方案：一种干扰信号与语音调制信号混叠的窄带干扰识别方法，其 特征在于，包括以下步骤：
[0007] a.信号分块：对接收的信号进行采样，然后将采样的信号进行分块；
[0008] b.窄带干扰预识别，具体包括以下步骤：
[0009] bl.将步骤a中分块后所得的块数据内时域信号进行傅里叶频域变换（FFT)变换 后得到接收信号的频域响应，然后与预设的接收信号频域响应幅度值门限进行比较，得出 大于门限值的频点；这些频点即为载频和干扰信号所在频点，低于门限值的频点为噪声和 语音信号的低频部分；频域响应幅度值门限的设定依赖于接收语音信号频域幅度值的统计 特征，首先统计块内信号频域幅度值的标准差，然后将频域响应幅度值门限设定为该标准 差的s倍，其中s的设置和信号与干扰特性有关。语音信号成分的频率幅度值往往小于载频 和窄带干扰的频率幅度值，通过该门限值即可将载频和干扰频点从语音信号中区分出来。 语音信号中的窄带干扰一般为单频干扰，其频谱宽度很窄。虽然载频信号一般也是单频信 号，但是经过语音信号调制过的载频信号，其频谱宽度得到扩展。因此可比较门限值以上的 频谱宽度，最大宽度的频谱对应的频率即为载频信号的频率，其余频谱对应的频率为干扰 信号频率。在某些情况下，由于FFT的栅栏效应，可能出现存在两个以上相同宽度的最大频 谱，此时通过频谱宽度的判别方法已无法判别，将在后续的窄带干扰最终识别阶段进行判 别。
[0010] b2.将大于门限值的频点进行分组，其中每一组的频点同属于一个干扰源或载频； 计算属于同一个干扰源或者载频的频点个数即可得到该干扰源或者载频的频谱宽度；根据 所得的频谱宽度，判断是否存在多个相同宽度的最大频谱，若是，则进入步骤c，若否，则将 所得的最大宽度的频谱对应的频率识别为载频信号的频率，其余频谱对应的频率为干扰信 号频率，大于门限值的频点分组的个数减去载频个数即为干扰源个数；
[0011] c.窄带干扰识别：重复步骤b，依次对步骤a中信号分块所得的每一块数据进行预 识别，每次完成步骤b后记录该次预识别的块数据的载频频率、干扰频率和干扰源个数；为 根据搜索已经记录的所有块数据的干扰源个数，将出现次数最多的干扰源个数作为本次干 扰源个数识别结果；干扰频率为根据搜索已经记录的所有块数据的干扰频率，出现次数最 多的干扰频率即为本次干扰信号的频率。由于FFT的栅栏效应，某些情况下在前面窄带干 扰预识别阶段利用频谱宽度的识别方法可能无效或识别错误，此时需要利用前面的窄带干 扰识别结果记录来修正最终识别结果。窄带干扰一般是由于接受设备造成的，在短时内具 有平稳性，即在几秒时间内，干扰源的个数和频率是固定的，因此利用前面的块记录中干扰 源和干扰频点出现频率最高的组合来对窄带干扰进行最终识别是可行的
[0012] 具体的，步骤a中对接收的信号进行采样的采样频率至少为信号载频的2倍，将采 样的信号进行分块后的块内样本数为2的指数倍。
[0013] 按照奈奎斯特采样定理，采样频率需为信号载频的2倍以上，为方便后续频域变 换，块内样本数需为2的指数倍，如128、256、512和1024等。块内样本数的设置需要根据 时延和性能要求进行，块内样本数越大，识别效果越好，但其造成的语音时延也越大。如采 样频率为31. 75KHz，每块样本数为512,则语音时延为16ms。
[0014] 具体的，步骤a中对接收的信号进行采样的采样频率至少为信号载频的3倍。
[0015] 然而高速移动的语音信号接收平台存在多普勒频移，本方案为了抑制多普勒频移 对采样的影响，采样频率需设定为信号载频的3倍以上。
[0016] 本发明的有益效果为，本发明首先对接收语音信号进行分块处理，然后对块内信 号做FFT变换，利用语音信号与干扰信号频域响应幅度和宽度上的差异性对语音载频、干 扰频率和干扰个数进行预识别，然后再利用前面的块记录中干扰源和干扰频点出现次数最 多的组合来对窄带干扰进行最终识别，块间记录的统计信息能进一步的提高正确识别结 果。由于本发明利用语音信号与干扰信号频域响应幅度和宽度上的差异性而不是频率值的 差异性来区分语音信号与干扰，该方法无需获知语音信号接收频率和语音接收设备的背景 先验信息，适应于语音信号与窄带干扰识别的所有场合（高速或低速接收平台），计算方法 简单，计算量小，易于实现，是一种低成本、高性能的盲语音通信窄带干扰识别技术。
【专利附图】
【附图说明】
[0017] 图1为本发明的方法流程示意图。
【具体实施方式】
[0018] 下面结合附图和实施例对本发明进行详细的描述
[0019] 如图1所示，本发明的方法首先将接收信号进行采样，按照奈奎斯特采样定理，采 样频率需为信号载频的2倍以上，然后将采样后的信号按时延和性能要求进行分块，为方 便后续频域变换，块内样本数需为2的指数倍，如128、256、512和1024等；对块内时域信号 进行傅里叶频域变换（FFT)得到接收信号的频域响应，根据设定的频域响应幅度值门限， 找到大于该门限值的频点。这些频点即为载频和干扰信号所在频点，低于门限值的频点为 噪声和语音信号的低频部分。测量高于门限值频点的宽度，最大宽度的频谱对应的频率即 为载频信号的频率，其余频谱对应的频率为干扰信号频率。若存在两个以上相同宽度的最 大频谱，则不作处理；记录上一步中估计的该块数据的载频频率、干扰频率、和干扰源个数， 搜索之前记录，出现次数最多的干扰源个数η作为本次干扰源个数识别结果，再次搜索之 前记录，出现次数最多的η个干扰频率即为干扰信号的频率，作为当前块数据干扰信号的 识别结果。
[0020] 本发明方法主要包括一下几个步骤：
[0021] Α、信号分块：首先将接收信号进行采样，然后将采样后的信号进行分块。按照奈奎 斯特采样定理，采样频率需为信号载频的2倍以上，然而高速移动的语音信号接收平台存 在多普勒频移，为了抑制多普勒频移对采样的影响，采样频率需设定为信号载频的3倍以 上。为方便后续频域变换，块内样本数需为2的指数倍，如128、256、512和1024等。块内 样本数的设置需要根据时延和性能要求进行，块内样本数越大，识别效果越好，但其造成的 语音时延也越大。如采样频率为31. 75ΚΗζ，每块样本数为512,则语音时延为16ms。
[0022] B、窄带干扰预识别：首先将块内时域信号进行FFT变换得到接收信号的频域响 应，然后根据接收信号的设定的频域响应幅度值门限，找出大于该门限值的频点。这些频点 即为载频和干扰信号所在频点，低于门限值的频点为噪声和语音信号的低频部分。测量高 于门限值频点的宽度，最大宽度的频谱对应的频率即为载频信号的频率，其余频谱对应的 频率为干扰信号频率。若存在两个以上相同宽度的最大频谱，则不作处理。频域响应幅度 值门限的设定依赖于接收语音信号频域幅度值的统计特征，首先统计块内信号频域幅度值 的标准差，然后将频域响应幅度值门限设定为该标准差的s倍，其中s的设置和信号与干扰 特性有关。语音信号成分的频率幅度值往往小于载频和窄带干扰的频率幅度值，通过该门 限值即可将载频和干扰频点从语音信号中区分出来。语音信号中的窄带干扰一般为单频干 扰，其频谱宽度很窄。虽然载频信号一般也是单频信号，但是经过语音信号调制过的载频信 号，其频谱宽度得到扩展。因此可比较门限值以上的频谱宽度，最大宽度的频谱对应的频率 即为载频信号的频率，其余频谱对应的频率为干扰信号频率。在某些情况下，由于FFT的栅 栏效应，可能出现存在两个以上相同宽度的最大频谱，此时通过频谱宽度的判别方法已无 法判别，将在后续的窄带干扰最终识别阶段进行判别。
[0023] C、窄带干扰最终识别：记录上一步中估计的该块数据的载频频率、干扰频率、和干 扰源个数，搜索之前记录，出现次数最多的干扰源个数η作为本次干扰源个数识别结果，再 次搜索之前记录，出现次数最多的η个干扰频率即为干扰信号的频率，作为当前块数据干 扰信号的识别结果。由于FFT的栅栏效应，某些情况下在前面窄带干扰预识别阶段利用频 谱宽度的识别方法可能无效或识别错误，此时需要利用前面的窄带干扰识别结果记录来修 正最终识别结果。窄带干扰一般是由于接受设备造成的，在短时内具有平稳性，即在几秒时 间内，干扰源的个数和频率是固定的，因此利用前面的块记录中干扰源和干扰频点出现频 率最高的组合来对窄带干扰进行最终识别是可行的。
[0024] 本发明的具体的处理流程为：
[0025] 将接受的语音信号以fs的采样频率进行采样，将采样后的时域语音信号分块，块 内样本数为N，第t个块内信号为 ：
[0026] yt = [y(t-1)N+1，…，ytN]T
[0027] 对块内时域信号进行FFT，将信号从时域变换到频域，频域输出为：
[0028] Yt = [Y1，…，YN] T
[0029] 统计块内信号频域幅度值的标准差：
【权利要求】
1. 一种干扰信号与语音调制信号混叠的窄带干扰识别方法，其特征在于，包括以下步 骤： a. 信号分块：对接收的信号进行采样，然后将采样的信号进行分块； b. 窄带干扰预识别，具体包括以下步骤： bl.将步骤a中分块后所得的块数据内时域信号进行FFT变换，得到接收信号的频域响 应，然后与预设的接收信号频域响应幅度值门限进行比较，得出大于门限值的频点； b2.测量大于门限值频点的宽度，判断是否存在多个相同宽度的最大频谱，若是，则进 入步骤c，若否，则将所得的最大宽度的频谱对应的频率识别为载频信号的频率，其余频谱 对应的频率为干扰信号频率； c. 窄带干扰识别：重复步骤b依次对步骤a所得的每一块数据进行预识别，每次完成 步骤b后记录该次预识别的块数据的载频频率、干扰频率和干扰源个数；其中，干扰源个数 为根据搜索已经记录的所有块数据的干扰源个数，将出现次数最多的干扰源个数作为本次 干扰源个数识别结果；干扰频率为根据搜索已经记录的所有块数据的干扰频率，出现次数 最多的干扰频率即为本次干扰信号的频率。
2. 根据权利要求1所述的一种干扰信号与语音调制信号混叠的窄带干扰识别方法，其 特征在于，步骤a中对接收的信号进行采样的采样频率至少为信号载频的2倍，将采样的信 号进行分块后的块内样本数为2的指数倍。
3. 根据权利要求1所述的一种干扰信号与语音调制信号混叠的窄带干扰识别方法，其 特征在于，步骤a中对接收的信号进行采样的采样频率至少为信号载频的3倍。
【文档编号】G10L15/20GK104361892SQ201410477810
【公开日】2015年2月18日 申请日期:2014年9月18日 优先权日:2014年9月18日
【发明者】黄际彦, 徐保根, 陆川, 万群, 饶中初, 王碧雯 申请人:电子科技大学, 同方电子科技有限公司, 电子科技大学成都研究院