CN113225138A - 一种消除水声强宽带干扰的方法 - Google Patents

一种消除水声强宽带干扰的方法 Download PDF

Info

Publication number
CN113225138A
CN113225138A CN202110446496.5A CN202110446496A CN113225138A CN 113225138 A CN113225138 A CN 113225138A CN 202110446496 A CN202110446496 A CN 202110446496A CN 113225138 A CN113225138 A CN 113225138A
Authority
CN
China
Prior art keywords
tap
signal
filter
broadband interference
underwater acoustic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110446496.5A
Other languages
English (en)
Other versions
CN113225138B (zh
Inventor
陈劼
韩冰
刘洁
马绪峰
安杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202110446496.5A priority Critical patent/CN113225138B/zh
Publication of CN113225138A publication Critical patent/CN113225138A/zh
Application granted granted Critical
Publication of CN113225138B publication Critical patent/CN113225138B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B11/00Transmission systems employing sonic, ultrasonic or infrasonic waves
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208Noise filtering
    • G10L21/0216Noise filtering characterised by the method used for estimating noise
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B13/00Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Human Computer Interaction (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Noise Elimination (AREA)

Abstract

本发明属于水声通信技术领域,具体涉及一种消除水声强宽带干扰的方法。本发明为了解决水声通信中现有的强宽带干扰消除方案的不足,提出了一种适用于水下声通信系统的自适应噪声消除方法。它可以指导噪声消除自适应滤波器的重要信道抽头,使信道均衡器不必跟踪整个信道,从而快速收敛并避免了目标信号的失真。本发明的有益效果为:可以在水声通信环境中,使滤波器抽头快速收敛,将水听器接收信号中的声波宽带干扰的消除,并有效避免目的信号失真,提高水下声通信可靠性。

Description

一种消除水声强宽带干扰的方法
技术领域
本发明属于水声通信技术领域,具体涉及一种消除水声强宽带干扰的方法。
背景技术
水声通信是目前常用的水下无线数据传输技术。但是,在一些水声通信场景中,信号是在恶劣的环境中记录的,例如封闭的港口或噪声较大的船只附近。如果接收机收到与目标信号有相同频带的强宽带干扰,系统信噪比会很低,通信的可靠性会大大下降。目前大多数的水声强宽带干扰的干扰消除解决方案主要考虑了对强干扰信号的缓解,并没有设计来避免目标信号的失真。此外,由于需要对水声信道进行跟踪,且水声传播具有延迟较长的特点,现有宽带干扰消除方案需要较长的收敛时间。
发明内容
针对上述问题,本发明为了解决水声通信中现有的强宽带干扰消除方案的不足,提出了一种适用于水下声通信系统的自适应噪声消除方法。它可以指导噪声消除自适应滤波器的重要信道抽头,使信道均衡器不必跟踪整个信道,从而快速收敛并避免了目标信号的失真。
本发明的技术方案是:
一种消除水声强宽带干扰的方法,其特征在于,包括以下步骤:
S1、对水听器的接收信号和参考向量进行预处理:
设Nf为自适应滤波器的阶数,NS为参考缓冲区中的样品数,并且N=NS/Nf是一个整数,将信号和参数向量划分为N个时间窗,基于递归最小二乘自适应滤波器,计算第n个时间窗的基础方程:
Figure BDA0003037107980000011
其中,rn为参考缓冲区的第n块,Covn为误差的协方差矩阵,λ为选定的标量;令sn为信号缓冲区的第n块,第n个时间窗的自适应滤波器权重wn间的关系式为:
Figure BDA0003037107980000021
S2、将水听器的接收信号通过归一化匹配滤波器,与抽头阈值比较,挑选出与水声宽带干扰信号相关性较高的抽头的索引:
初始化设置抽头索引,组成集合Tw={t1,t2,…tL},其中t1,t2,…tL为所标识的有效抽头的位置;基于正则非完全β函数计算出虚警概率Pt,得出合适的抽头门限tThr
Figure BDA0003037107980000022
令水听器的接收信号通过归一化匹配滤波器,实现与干扰信号匹配,结果序列为M(i),1≤i≤L,选择绝对值高于门限的抽头的索引,得到有效抽头索引序列为Tv={…,ti,…,tj,…},1≤i,j≤L;
S3、根据得到的有效抽头的索引,初始化自适应滤波器抽头权重为:
Figure BDA0003037107980000023
S4、基于步骤S1中提出的改进的递归最小二乘法的基础方程,对滤波器抽头进行迭代计算,直至抽头权重序列收敛:
Figure BDA0003037107980000024
水听器接收到的信号经过滤波器抽头的递归过滤,得到完整的目的信号。
本发明的有益效果为:可以在水声通信环境中,使滤波器抽头快速收敛,将水听器接收信号中的声波宽带干扰的消除,并有效避免目的信号失真,提高水下声通信可靠性。
附图说明
图1为本发明的流程示意图。
具体实施方式
下面结合附图,对本发明技术方案进行详细描述:
如图1所示,本发明的方法具体包括:
S1、对水听器的接收信号和参考向量进行预处理,划分时间窗,计算出递归最小二乘法的基本方程。具体方法为:
设Nf为自适应滤波器的阶数,NS为参考缓冲区中的样品数。假设N=NS/Nf是一个整数值,将信号和参考向量划分为N个时间窗。基于递归最小二乘自适应滤波器,计算第n个时间窗的基础方程:
Figure BDA0003037107980000031
其中rn为参考缓冲区的第n块,Covn为误差的协方差矩阵,λ为选定的标量。sn为信号缓冲区的第n块,第n个时间窗的自适应滤波器权重wn间的关系式为:
Figure BDA0003037107980000032
S2、将水听器的接收信号通过归一化匹配滤波器,与抽头阈值比较,挑选出与水声宽带干扰信号相关性较高的抽头的索引。具体方法为:
初始化设置抽头索引,组成集合Tw={t1,t2,…tL},其中t1,t2,…tL为所标识的有效抽头的位置。基于正则非完全β函数计算出虚警概率Pt,得出合适的抽头门限tThr
Figure BDA0003037107980000033
之后让水听器的接收信号通过归一化匹配滤波器,实现与干扰信号匹配,结果序列为M(i),1≤i≤L。只选择绝对值高于门限的抽头的索引,使抽头可以更针对性地消除干扰信号,同时降低对目标信号的失真影响。得到有效抽头索引序列为Tv={…,ti,…,tj,…},1≤i,j≤L。
S3、根据得到的有效抽头的索引,并结合信号经归一化匹配滤波器后得到的与干扰信号强相关的重要路径值,初始化滤波器抽头。并将其它相关度弱的索引处的抽头值设为零,降低对目标信号的失真影响。具体方法为:
基于归一化匹配滤波器得出的与干扰信号强相关的抽头索引,定制化滤波器的抽头权重。抓住滤波器的抽头系数集合Tw与误差协方差矩阵Cov1有比例关系的特点,利用上述步骤得出的与干扰信号强相关的抽头索引,初始化自适应滤波器抽头权重为:
Figure BDA0003037107980000041
S4、基于步骤S1中提出的改进的递归最小二乘法的基础方程,迭代计算修正滤波器抽头,至抽头逐渐收敛进入稳定状态。具体方法为:
依照改进的递归最小二乘法,对滤波器抽头进行迭代计算,至抽头权重序列收敛。
Figure BDA0003037107980000042
最终,水听器接收到的信号经过滤波器抽头的递归过滤,得到完整的目的信号。在此过程中抽头序列通过不断递归迅速收敛,提高了消除水声宽带干扰信号的速度,并避免了目的信号的失真。

Claims (1)

1.一种消除水声强宽带干扰的方法,其特征在于,包括以下步骤:
S1、对水听器的接收信号和参考向量进行预处理:
设Nf为自适应滤波器的阶数,NS为参考缓冲区中的样品数,并且N=NS/Nf是一个整数,将信号和参数向量划分为N个时间窗,基于递归最小二乘自适应滤波器,计算第n个时间窗的基础方程:
Figure FDA0003037107970000011
其中,rn为参考缓冲区的第n块,Covn为误差的协方差矩阵,λ为选定的标量;令sn为信号缓冲区的第n块,第n个时间窗的自适应滤波器权重wn间的关系式为:
Figure FDA0003037107970000012
S2、将水听器的接收信号通过归一化匹配滤波器,与抽头阈值比较,挑选出与水声宽带干扰信号相关性较高的抽头的索引:
初始化设置抽头索引,组成集合Tw={t1,t2,…tL},其中t1,t2,…tL为所标识的有效抽头的位置;基于正则非完全β函数计算出虚警概率Pt,得出合适的抽头门限tThr
Figure FDA0003037107970000013
令水听器的接收信号通过归一化匹配滤波器,实现与干扰信号匹配,结果序列为M(i),1≤i≤L,选择绝对值高于门限的抽头的索引,得到有效抽头索引序列为Tv={…,ti,…,tj,…},1≤i,j≤L;
S3、根据得到的有效抽头的索引,初始化自适应滤波器抽头权重为:
Figure FDA0003037107970000014
S4、基于步骤S1中的基础方程,对滤波器抽头进行迭代计算,直至抽头权重序列收敛:
Figure FDA0003037107970000021
水听器接收到的信号经过滤波器抽头的递归过滤,得到完整的目的信号。
CN202110446496.5A 2021-04-25 2021-04-25 一种消除水声强宽带干扰的方法 Active CN113225138B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110446496.5A CN113225138B (zh) 2021-04-25 2021-04-25 一种消除水声强宽带干扰的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110446496.5A CN113225138B (zh) 2021-04-25 2021-04-25 一种消除水声强宽带干扰的方法

Publications (2)

Publication Number Publication Date
CN113225138A true CN113225138A (zh) 2021-08-06
CN113225138B CN113225138B (zh) 2022-01-25

Family

ID=77088778

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110446496.5A Active CN113225138B (zh) 2021-04-25 2021-04-25 一种消除水声强宽带干扰的方法

Country Status (1)

Country Link
CN (1) CN113225138B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101567863A (zh) * 2008-04-24 2009-10-28 魏昕 浅海水声通信系统的间接自适应均衡方法
CN106712781A (zh) * 2016-12-08 2017-05-24 浙江大学 从数字领域和模拟领域抑制干扰的带内全双工水声通信机
WO2017117732A1 (zh) * 2016-01-05 2017-07-13 北京司响无限文化传媒有限公司 一种基于信道特征的空气声波通信信道均衡方法和装置
CN108900443A (zh) * 2018-06-05 2018-11-27 中国科学院声学研究所 一种水声通信中的水声信道干扰消除方法
CN111092833A (zh) * 2019-12-30 2020-05-01 哈尔滨工程大学 一种全双工水声通信期待信道与自干扰信道联合估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101567863A (zh) * 2008-04-24 2009-10-28 魏昕 浅海水声通信系统的间接自适应均衡方法
WO2017117732A1 (zh) * 2016-01-05 2017-07-13 北京司响无限文化传媒有限公司 一种基于信道特征的空气声波通信信道均衡方法和装置
CN106712781A (zh) * 2016-12-08 2017-05-24 浙江大学 从数字领域和模拟领域抑制干扰的带内全双工水声通信机
CN108900443A (zh) * 2018-06-05 2018-11-27 中国科学院声学研究所 一种水声通信中的水声信道干扰消除方法
CN111092833A (zh) * 2019-12-30 2020-05-01 哈尔滨工程大学 一种全双工水声通信期待信道与自干扰信道联合估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王安义等: "水声软频率复用网络的干扰缓解与资源分配", 《科学技术与工程》 *

Also Published As

Publication number Publication date
CN113225138B (zh) 2022-01-25

Similar Documents

Publication Publication Date Title
CN107993670B (zh) 基于统计模型的麦克风阵列语音增强方法
CN112735456B (zh) 一种基于dnn-clstm网络的语音增强方法
US8014230B2 (en) Adaptive array control device, method and program, and adaptive array processing device, method and program using the same
CN106301289A (zh) 利用自适应滤波算法消除泥浆脉冲信号中的泵冲噪声的方法
CN109977724B (zh) 一种水下目标分类方法
CN101876715A (zh) 一种拖曳声学阵列的拖船噪声抑制方法
CN106059531B (zh) 一种非负自适应滤波器
CN105702262A (zh) 一种头戴式双麦克风语音增强方法
CN113225138B (zh) 一种消除水声强宽带干扰的方法
CN111090089B (zh) 一种基于两类辅助数据的空时自适应检测方法
Peyvandi et al. SONAR systems and underwater signal processing: classic and modern approaches
CN107635181B (zh) 一种基于信道学习的多址感知源反馈优化方法
CN111239814B (zh) 基于同相轴分频追踪平滑的浅剖数据机械干扰压制方法
CN107976663A (zh) 一种基于子空间投影的外辐射源雷达对目标压缩检测方法
EP2045620A1 (en) Acoustic propagation delay measurement
CN112329523A (zh) 一种水声信号种类识别方法、系统及设备
CN116520333A (zh) 一种基于复杂干扰下的多基浮标融合检测方法和系统
CN116047488A (zh) 一种针对压制式干扰源的主动声纳干扰抑制方法
Lampl et al. Noise cancellation with LMS, NLMS and RLS filtering algorithms to improve the fault detection of an industrial measurement system
CN108667538B (zh) 复杂噪声环境下时延和幅度衰减的联合估计方法
RU2616357C1 (ru) Способ регистрации малошумного морского объекта с использованием медианной фильтрации
CN114035193A (zh) 一种两阶可调谐空时自适应检测方法
Yu et al. The application of improved variable step-size LMS algorithm in sonar signal processing
CN109407078B (zh) 一种双自适应拷贝相关器方法
CN113794489B (zh) 一种通信抗强相关干扰的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant