首页 / 专利库 / 数学与统计 / 自相关矩阵 / 一种适用于矢量水听器的线谱检测方法

一种适用于矢量听器的线谱检测方法

阅读:770发布:2020-05-08

专利汇可以提供一种适用于矢量听器的线谱检测方法专利检索,专利查询,专利分析的服务。并且本 发明 提出一种适用于矢量 水 听器 的线谱检测方法,其包括步骤:步骤一,对矢量阵列各 传感器 通道的时域 信号 做波束形成处理;步骤二,对时域波束形成的结果做 带通滤波 和降 采样 处理;步骤三,通过对多种自适应线谱增强 算法 的对比分析;步骤四,使用多重自相关法,利用信号与噪声在统计特性上的可区分性,把淹没于噪声中的线谱信号提取出来;步骤五,使用分段周期图平均法对通过自适应线谱增强和多重自相关处理的时间序列做 功率谱 估计;步骤六,对 门 限检测的结果做时间统计。该线谱检测方法充分利用矢量水听器的空间增益和时间增益来提高线谱信号的 信噪比 ,自动识别与提取出淹没在宽带干扰背景中的目标线谱。,下面是一种适用于矢量听器的线谱检测方法专利的具体信息内容。

1.一种适用于矢量听器的线谱检测方法,其特征在于,包括如下步骤:
步骤一,对矢量阵列各传感器通道的时域信号做波束形成处理;
步骤二,对时域波束形成的结果做带通滤波和降采样处理,为自适应滤波增强处理效果和提高处理速度做准备;
步骤三,通过对多种自适应线谱增强算法的对比分析,简化的谐波簇自适应线谱增强算法在低信噪比情况下能有效地抑制背景噪声;
步骤四,使用多重自相关法,利用信号与噪声在统计特性上的可区分性,把淹没于噪声中的线谱信号提取出来;
步骤五,使用分段周期图平均法对通过自适应线谱增强和多重自相关处理的时间序列做功率谱估计;
步骤六,对限检测的结果做时间统计。
2.如权利要求1所述的适用于矢量水听器的线谱检测方法,其特征在于,所述步骤一中的波束形成处理采用的计算公式为B=W(θ)*A,其中 为校准因子,θ为扫描度在0 360°范~
围内取值,A为阵列流形,
在每个扫描角度下对数据做步骤一到步骤五的处理,最终得到过门限谱的频率方位矩阵 ,对信号处理频带内的每一个频率 对应的N个处理结果 取最大值 作为后续时间统计的输入。
3.如权利要求2所述的适用于矢量水听器的线谱检测方法,其特征在于,在所述步骤二中,为了避免降采样后出现频谱混叠,根据采样定理的要求,带通滤波的高截止频率fh和降采样后的采样率Fs的关系应该满足Fs>2*fh的条件。
4.如权利要求3所述的适用于矢量水听器的线谱检测方法,其特征在于,在所述步骤三中,简化的谐波簇的自适应线谱增强算法以一对频率为 的正交信号为参考信号,的最大取值为声纳的最大工作频率 ,参考信号与权系数w相乘求和得到Y(k)即为算法的输出,使用输入p(k)与输出Y(k)之差来更新权系数使得自适应过程持续下去,权系数更新迭代公式如下:
………………1
………………2
………………3
当正交参考信号的频率 与输入信号中的线谱频率对应时,输入信号的线谱分量即可被增强,自适应步长 影响着算法的噪声抑制效果和响应速度,应根据实际情况选择合适的值。
5.如权利要求4所述的适用于矢量水听器的线谱检测方法,其特征在于,在所述步骤四中,设信号为 ,其中 为高斯白噪声,自相关
函数为: ,
其中 为正弦信号的自相关函数,通过公式推导可得 ,
、 、 分别为噪声与信号、噪声与噪声的相关函数,
若噪声为标准白噪声,则噪声与信号不自相关,但实际测量中,观测时间是有限的,且噪声白化程度不一定十分理想,从而导致自相关函数 、 、 不为零,可将
当作新的噪声Rn(t),其幅值与噪声n(t)相比很小,多重自相关法是将 当作 重复上述步骤多次。
6.如权利要求5所述的适用于矢量水听器的线谱检测方法,其特征在于,在所述步骤五中,将Ns的数据分为2*(N-1)+1段,每段包含1s的数据,各段子序列重叠率50%,计算出每段的周期图,然后加权平均得到功率谱估计,获得功率谱估计之后采用双向α滤波获取功率谱的连续谱趋势,去趋势化后设定门限与直接设定固定门限的方法相比,可以更好的剔除连续谱背景对线谱检测的影响,避免因固定门限设置不合理带来的大片连续谱信息残留。
7.如权利要求6所述的适用于矢量水听器的线谱检测方法,其特征在于,在所述步骤六中,设count(fi)为动态累计数组,若当前时刻频率fi对应的功率谱值p(fi)通过门限检测后不为零则对应的count(fi)= count(fi)+1,否则count(fi)=count(fi)-1,然后,设置门限值为C,对动态累计结果count(fi)满足大于等于C的条件,则说明存在线谱,其频率为fi。

说明书全文

一种适用于矢量听器的线谱检测方法

[0001]

技术领域

[0002] 本发明涉及一种适用于矢量水听器的线谱检测方法,属于水声信号处理技术领域。
[0003]

背景技术

[0004] 矢量水听器是接收换能器的一种。在水声领域,通常,将传感器称为换能器,接收换能器主要包括标量传感器和矢量传感器,也叫标量水听器和矢量水听器。在声场测量中,传统的方法是采用标量水听器,只能测量声场中的标量参数。矢量水听器能测量声场中的矢量参数,它的应用有助于获得声场的矢量信息,对声呐设备的功能扩展具有极为关键的意义。
[0005] 舰船辐射噪声一般包含线谱和连续谱,其中线谱主要是由于舰艇螺旋桨的周期性击水以及叶片共振产生的。因产生线谱的声源的功率和惯性都相当大,工作条件比较稳定,所以线谱具有较高的能量稳定性,使用线谱检测技术可以显著提高被动声呐的作用距离。近年来,现代舰艇的隐身性能得到了高度重视,随着减振降噪技术在舰艇上的大量应用,目标辐射噪声大幅降低,出现了很多安静型目标。从强海洋背景噪声中检测出微弱的线谱信号,成为近年来水声信号处理的一个重要的研究课题。一些理论成熟且行之有效的方法已经被应用在微弱线谱检测的研究中。
[0006] 波束形成器可以获得空间处理增益,为检测弱信号创造条件。波束形成器是一种空域滤波器,它将各通道数据进行加权求和得到输出,通过设计合适的加权值,可以有选择地增强来自某一指定方向的信号,抑制其他方向到达的信号,提高输出信噪比
[0007] 自适应线谱增强是自适应信号处理的一个重要应用,它具有在低信噪比条件下较好地提取线谱的能,且无需噪声参考信号就能自适应地将线谱信号从宽带噪声中分离出来。但是该算法的固有缺陷严重制约了其工程实际应用。
[0008] 因此,如何利用矢量水听器实现自动识别与提取淹没在强海洋背景噪声中的远距离或安静型目标辐射的微弱线谱信号,还有待进一步地研究。
[0009]

发明内容

[0010] 为解决现有技术存在的不足,本发明旨在提出一种适用于矢量水听器的线谱检测方法,其充分利用矢量水听器的空间增益和时间增益来提高线谱信号的信噪比,自动识别与提取出淹没在宽带干扰背景中的目标线谱。
[0011] 为实现上述目的,本发明的适用于矢量水听器的线谱检测方法包括如下步骤:步骤一,对矢量阵列各传感器通道的时域信号做波束形成处理;
步骤二,对时域波束形成的结果做带通滤波和降采样处理,为自适应滤波增强处理效果和提高处理速度做准备;
步骤三,通过对多种自适应线谱增强算法的对比分析,简化的谐波簇自适应线谱增强算法在低信噪比情况下能有效地抑制背景噪声;
步骤四,使用多重自相关法,利用信号与噪声在统计特性上的可区分性,把淹没于噪声中的线谱信号提取出来;
步骤五,使用分段周期图平均法对通过自适应线谱增强和多重自相关处理的时间序列做功率谱估计;
步骤六,对限检测的结果做时间统计。
[0012] 进一步地,步骤一中的波束形成处理采用的计算公式为B=W(θ)*A,其中 为校准因子,θ为扫描度在0 360°范围内取值,A为阵列流形,~
在每个扫描角度下对数据做步骤一到步骤五的处理,最终得到过门限谱的频率方位矩阵 ,对信号处理频带内的每一个频率 对应的N个处理结果 取最大值 作为后续时间统计的输入。
[0013] 进一步地,在步骤二中,为了避免降采样后出现频谱混叠,根据采样定理的要求,带通滤波的高截止频率fh和降采样后的采样率Fs的关系应该满足Fs>2*fh的条件。
[0014] 进一步地,在步骤三中,简化的谐波簇的自适应线谱增强算法以一对频率为 的正交信号为参考信号,的最大取值为声纳的最大工作频率 ,参考信号与权系数w相乘求和得到Y(k)即为算法的输出,使用输入p(k)与输出Y(k)之差来更新权系数使得自适应过程持续下去,权系数更新迭代公式如下:………………1
………………2
………………3
当正交参考信号的频率 与输入信号中的线谱频率对应时,输入信号的线谱分量即可被增强,自适应步长 影响着算法的噪声抑制效果和响应速度,应根据实际情况选择合适的值。
[0015] 进一步地,在步骤四中,设信号为 ,其中 为高斯白噪声,自相关函数为: ,
其中 为正弦信号的自相关函数,通过公式推导可得 ,
、 、 分别为噪声与信号、噪声与噪声的相关函数。若噪声为标准白噪声,则噪声与信号不自相关,但实际测量中,观测时间是有限的,且噪声白化程度不一定十分理想,从而导致自相关函数 、 、 不为零,可将 当作新的噪声Rn
(t),其幅值与噪声n(t)相比很小,多重自相关法是将 当作 重复上述步骤多次。
[0016] 进一步地,在步骤五中,将Ns的数据分为2*(N-1)+1段,每段包含1s的数据,各段子序列重叠率50%,计算出每段的周期图,然后加权平均得到功率谱估计,获得功率谱估计之后采用双向α滤波获取功率谱的连续谱趋势,去趋势化后设定门限与直接设定固定门限的方法相比,可以更好的剔除连续谱背景对线谱检测的影响,避免因固定门限设置不合理带来的大片连续谱信息残留。
[0017] 进一步地,在步骤六中,设count(fi)为动态累计数组,若当前时刻频率fi对应的功率谱值p(fi)通过门限检测后不为零则对应的count(fi)= count(fi)+1,否则count(fi)=count(fi)-1,然后,设置门限值为C,对动态累计结果count(fi)满足大于等于C的条件,则说明存在线谱,其频率为fi。
[0018] 本发明的适用于矢量水听器的线谱检测方法在波束形成的基础上,结合自适应线谱增强、多重自相关、自动门限检测等算法,充分利用矢量水听器的空间增益和时间增益来提高线谱信号的信噪比,自动识别与提取出淹没在宽带干扰背景中的目标线谱。
[0019]附图说明
[0020] 下面结合附图对本发明作进一步描写和阐述。
[0021] 图1是本发明的适用于矢量水听器的线谱检测方法的流程图
[0022] 图2是图1中线谱检测方法步骤三中简化的基于谐波簇的自适应线谱增强算法流程图。
[0023] 图3是图1中线谱检测方法步骤五中双向α滤波流程图。
[0024]

具体实施方式

[0025] 下面将结合附图、通过对本发明的优选实施方式的描述,更加清楚、完整地阐述本发明的技术方案。
[0026] 本发明提出一种适用于矢量水听器的线谱检测方法,其基本思路是:首先,在时域上利用波束形成得到对线谱目标的空间增益;然后,使用自适应线谱增强、多重自相关等算法利用信号相关而噪声不相关的特点,提取出淹没在噪声中的相关成分;最后,使用频谱估计、双向α滤波等手段实现自动门限检测提取出相关成分中的线谱,根据线谱存在的时长来判断环境中是否有线谱声目标,并提取出目标的频率。
[0027] 如图1所示,本发明的适用于矢量水听器的线谱检测方法包括如下布置:步骤一,对矢量阵列各传感器通道的时域信号做波束形成处理。采用常规波束形成方法实现波束的计算,计算公式为B=W(θ)*A,其中 为校准因子,θ为扫描角度在0 360°范围~
内取值,为阵列流形。
[0028] 在每个扫描角度下对数据做步骤一到步骤五的处理,最终得到过门限谱的频率方位矩阵 。对信号处理频带内的每一个频率 对应的N个处理结果 取最大值作为后续时间统计的输入。
[0029] 步骤二,对时域波束形成的结果B做带通滤波和降采样处理,为自适应滤波增强处理效果和提高处理速度做准备。为了避免降采样后出现频谱混叠,根据采样定理的要求,带通滤波的高截止频率fh和降采样后的采样率Fs的关系应该满足Fs>2*fh的条件。
[0030] 步骤三,通过对多种自适应线谱增强算法的对比分析,发现简化的谐波簇自适应线谱增强算法在低信噪比情况下仍可以有效地抑制背景噪声。简化的谐波簇的自适应线谱增强算法的原理结构如图2所示,算法以一对频率为 的正交信号为参考信号,的最大取值为声纳的最大工作频率 ,参考信号与权系数w相乘求和得到Y(k)即为算法的输出,使用输入p(k)与输出Y(k)之差来更新权系数使得自适应过程持续下去,权系数更新迭代公式如下。
[0031] ………………1………………2
………………3
当正交参考信号的频率 与输入信号中的线谱频率对应时,输入信号的线谱分量即可被增强。自适应步长 影响着算法的噪声抑制效果和响应速度,应根据实际情况选择合适的值。
[0032] 步骤四,使用多重自相关法,利用信号与噪声在统计特性上的可区分性,把淹没于噪声中的线谱信号提取出来。
[0033] 设信号为 ,其中 为高斯白噪声,自相关函数为: 。
[0034] 其中 为正弦信号的自相关函数,通过公式推导可得 ,尽管与原正弦信号的幅值和相位不同,频率却没有发声变化。
[0035] 、 、 分别为噪声与信号、噪声与噪声的相关函数。若噪声为标准白噪声,则噪声与信号不自相关,但实际测量中,观测时间是有限的,且噪声白化程度不一定十分理想,从而导致自相关函数 、 、 不为零,可将 当作新的噪声Rn(t),其幅值与噪声n(t)相比很小。
[0036] 可见,通过自相关处理提高了线谱信号的信噪比,但改善程度有限,多重自相关法是将 当作 重复上述步骤多次,理论上自相关次数越多,信噪比提高的越多,因此可检测出淹没于噪声中的微弱线谱信号。
[0037] 步骤五,使用分段周期图平均法对通过自适应线谱增强和多重自相关处理的时间序列做功率谱估计。将Ns的数据分为2*(N-1)+1段,每段包含1s的数据,各段子序列重叠率50%,计算出每段的周期图,然后加权平均得到功率谱估计。获得功率谱估计之后采用双向α滤波获取功率谱的连续谱趋势,去趋势化后设定门限与直接设定固定门限的方法相比,可以更好的剔除连续谱背景对线谱检测的影响,避免因固定门限设置不合理带来的大片连续谱信息残留。双向α滤波流程如图3所示。
[0038] 步骤六,对门限检测的结果做时间统计,设count(fi)为动态累计数组,若当前时刻频率fi对应的功率谱值p(fi)通过门限检测后不为零则对应的count(fi)= count(fi)+1,否则count(fi)=count(fi)-1。
[0039] 然后,设置门限值为C,对动态累计结果count(fi)满足大于等于C的条件,则说明存在线谱,其频率为fi。
[0040] 上述具体实施方式仅仅对本发明的优选实施方式进行描述,而并非对本发明的保护范围进行限定。在不脱离本发明设计构思和精神范畴的前提下,本领域的普通技术人员根据本发明所提供的文字描述、附图对本发明的技术方案所作出的各种变形、替代和改进,均应属于本发明的保护范畴。本发明的保护范围由权利要求确定。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈