一种洋流干扰中微弱目标信号提取方法
【专利摘要】本发明提出一种洋流干扰中微弱目标信号提取方法，是按照下述方式实现的：获取海洋磁场变化信息f0i；对f0i进行多尺度分解得到信号f1i；对f1i采取随机共振方法；对非线性随机共振系统利用四阶龙格?库塔方法进行高精度数值化求解；找到最优的随机共振参数；求得最佳的共振输出数据序列f2i；对f2i进行多尺度分解；在目标信号所在的主要频带区间选取不同的阈值，重构出经过阈值处理的信号f3i；对f3i进行特征提取；把提取的特征参数，与水下目标运动的数据库进行匹配比对。本发明能从复杂海洋背景电磁干扰中甄别出异常的目标扰动磁场特征，从而间接探测、发现和识别过境的水下动目标，有效提高水下探测、预警能力。
【专利说明】
-种洋流干扰中微弱目标信号提取方法
技术领域
[0001] 本发明属于水下动目标探测技术领域，具体设及一种利用随机共振与小波分析技 术实现对水下动目标信号的提取与识别。
【背景技术】
[0002] 水下动目标的探测和定位不论在民用还是在军事领域，一直都是备受关注的问 题。在不同环境下针对不同的水下动目标，各种探测方法各有其优势和不足。例如，目前最 基本的方法是声纳探测，但是它的虚警率高，且在探测海洋中的目标时受到很大限制。光学 探测技术容易受到恶劣环境气象条件的影响。磁场异常探测技术受到环境影响较小，可W 弥补其他探测技术留下的探测盲区，但磁异信号探测距离有限。
[0003] 近年来，随着减振降噪技术和消磁、消电技术被大量运用，使得水下动目标本身的 声、磁、电信号更难捕捉，寻找一种新的信号源和信号提取识别方法变得十分重要和迫切。 海水是一种弱导电介质，水下动目标在航行过程中会对海水造成端流扰动，运种扰动与地 球磁力线切割会产生很微弱电磁场变化。通常，目标运动引起的微弱电磁场变化与海洋本 身的电磁场变化叠加在一起，难W分离，被当成干扰和噪声信号。从复杂海洋电磁干扰中提 取出微弱目标信号十分困难，受到数据提取处理分析能力的限制，很长时间只是停留在理 论阶段，需不断探索研究新机制、新方法。
[0004] 通常的信号检测方法都认为噪声是有害的，一般将噪声进行抑制，采用各种方式 抑制或滤除噪声，然后把有用信号提取出来;新兴的随机共振方法，通过系统非线性，利用 噪声增强微弱信号，把噪声变害为宝，实现部分噪声能量向特征信号能量转移，使得在削弱 噪声的同时，又强化了待测信号，从而实现微弱信号的有效检测。小波分析方法是一种新型 的信号分析工具，它将时域和频域有效结合起来，具有多分辨率、方向选择性和自动聚焦等 特点，在信号检测、特征提取等领域有广泛应用。
[0005] 本发明巧妙利用随机共振与小波分析方法各自的优点，建立了基于随机共振机制 及小波分析方法相结合的目标洋流干扰磁信号提取与识别方法，在强干扰噪声中微弱信号 检测方面现出独特的优势和巨大魅力，开拓了一种水下目标检测与识别的新思路和新方 法。
【发明内容】
[0006] 本发明的目的是针对现有技术的不足，提出一种洋流干扰中微弱目标信号提取方 法，能够准确的需找海洋中的信号源，提取相应的信号信息。
[0007] 本发明的洋流干扰中微弱目标信号提取方法，是按照下述方式实现的：
[000引步骤1，通过置于水底的磁场收集天线获取海洋磁场变化信息，采集海浪波动、洋 流流动、地球磁场变化和动目标经过产生的磁场扰动信息，采集到的数据表示为f01，其中i =1 ???N，N为数据点个数，采样频率为fS Hz ;
[0009]步骤2,采用小波分析技术对步骤1采集到的数据进行多尺度分解，去除采集仪器 本身高频噪声及环境工频干扰噪声之后得到信号f 11;
[0010] 步骤3,对步骤2得到的信号fii采取随机共振方法，利用随机共振系统的非线性特 性，将海洋背景电磁噪声的部分能量转化为水下动目标扰动信号的能量，削弱噪声的同时， 强化待测信号，提高动目标信号的信噪比和检出率；
[0011] 步骤4,对步骤3中的非线性随机共振系统，利用四阶龙格-库塔方法进行高精度数 值化求解，求解步骤如下：
[0017] 式（1.4)-(1.8)中扣和9。分另懐示共振输出义(*)和输入信号的第11个采样值，口。= fln，h=l/fs为时间步长，即数据采样频率fs的倒数;a，b为随机共振系统的系统结构参数； a, b的取值范围与检测信号中的噪声类型W及信号强度等相关，通常取值范围可选择[Ο? ι.0]。
[0018] 步骤5,利用步骤4的方法，通过调节随机共振系统的系统结构参数a, b，找到最优 的随机共振参数a*和b*;基于最优共振参数a*和b*，利用步骤4的方法，求得最佳的共振输出 数据序列，记为f 21;
[0019] 步骤6,采用小波分析多尺度分解方法对步骤5随机共振处理之后的数据f2i进行多 尺度分解，依次把f2i中的各种频率成分从高向低逐步分离为不同的频带系数fj-乂 W-〉 fj-2->---->fl ；
[0020] 步骤7,在目标信号所在的主要频带区间[fa-fb]选取不同的阔值，对小于阔值的 系数置零，对大于阔值的系数进行非线性处理，而后利用小波重构方法重构出经过阔值处 理的信号，重构后的数据记为f3i;
[0021] 步骤8,对步骤7得到的数据f3i进行特征提取，得到信号波动周期、信号衰减长度、 信号峰值信息；
[0022] 步骤9,把步骤8提取的特征参数，与水下目标运动的数据库进行匹配比对，从而识 别出最相似的目标类型，实现对目标的分类和特征识别；
[0023] 步骤10，在数据处理终端对探测到的目标进行预警提示，告知目标的位置、深度、 大小、类型参数。
[0024] 步骤3采取的随机共振系统是由郎之万方程描述： 闺
（1.0
[0029] 式（1.3)中X为系统输出，s(t)为目标扰动引起的周期驱动信号，丫（t)为噪声信 号，a, b为随机共振系统的系统结构参数。
[0030] 步骤2中去除采集仪器本身高频噪声及环境工频干扰噪声的方式是将高于50化的 信号分量设为零，剩余信号分量之和记为f 11。
[0031] 本发明通过置于水底的磁场信号收集天线对海洋的磁场变化进行长时间不间断 监测，采用随机共振与小波分析相结合的信号提取技术，能从复杂海洋背景电磁干扰中甄 别出异常的目标扰动磁场特征，从而间接探测、发现和识别通过的水下动目标。本发明可用 于近海海岸防御、港口防御及海岛防御等，有效提高水下探测、预警能力。
【附图说明】
[0032] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W 根据运些附图获得其他的附图。
[0033] 图1为本发明的流程图。
[0034] 图2(a)真实海洋环境中磁场收集天线的检测数据foi。
[0035] 图2(b)为去除仪器噪声及工频干扰后的信号数据fii。
[0036] 图3(a)为采用第一组随机共振系数的信号数据f2i。
[0037] 图3(b)为采用第二组随机共振系数的信号数据f2i。
[0038] 图3(c)为采用第Ξ组随机共振系数的信号数据f2i。
[0039] 图3(d)为采用第四组随机共振系数的信号数据f2i。
[0040] 图4为经过多尺度的小波重构法之后得到的信号数据f3i。
【具体实施方式】
[0041] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0042] 本发明针对水下动目标探测防御的需求，可W巧妙利用海洋本身的电磁干扰噪 声，实现噪声与目标信号的随机共振，把噪声能量用来增强信号幅度，从洋流干扰中挖掘出 有用的目标信号，用于鉴别目标的存在，确定目标特征参数等，从而实现磁隐身、弱磁甚至 无磁水下目标的间接探测。本发明可W有效拓展我国水下侦测预警手段能力，提升我国水 下监测水平。
[0043] 如图1所示，本发明的洋流干扰中微弱目标信号提取方法，是按照下述方式实现 的：
[0044] 步骤1，通过置于水底的磁场收集天线长时间、不间断地获取海洋磁场变化信息， 采集海浪波动、洋流流动、地球磁场变化和可能的动目标经过产生的磁场扰动信息，采集到 的数据表示为foi，其中i = 1 ···Ν，Ν为数据点个数，采样频率为fS Hz;
[0045] 步骤2,采用小波分析技术对步骤1采集到的数据进行多尺度分解，将高于50化的 信号分量设为零，剩余信号分量之和记为f 11 α=1 ???N); f 11 α=1 ???N)即为去除采集仪器本 身高频噪声及环境工频干扰噪声之后的信号;如图2(b)所示；
[0046] 步骤3,对步骤2得到的信号fii采取随机共振方法，利用随机共振系统的非线性特 性，将海洋背景电磁噪声的部分能量转化为水下动目标扰动信号的能量，削弱噪声的同时， 强化待测信号，提高动目标信号的信噪比和检出率；
[0047] 随机共振系统可W采用多种方程描述，如按照下述方式采用郎之万化angevin)方 程描述：
[0052] 式（1.3)中X为系统输出，s(t)为目标扰动引起的周期驱动信号，丫（t)为噪声信 号，a, b为随机共振系统的系统结构参数。
[0053] 步骤4,对步骤3中的非线性随机共振系统，利用四阶龙格-库塔(Runge-Kutta)方 法进行高精度数值化求解，求解步骤如下：
[0059] 式（1.4)-(1.8)中Xn和Pn分别表示共振输出x(t)和输入信号的第η个采样值，Pn = fln，h=l/fs为时间步长，即数据采样频率fs的倒数;a，b为随机共振系统的系统结构参数；
[0060] 步骤5,利用步骤4的方法，通过调节随机共振系统的系统结构参数a, b，找到最优 的随机共振参数a*和b*;基于最优共振参数a*和b*，利用步骤4的方法，求得最佳的共振输出 数据序列，记为f2i，其中a、b-般在0.01~1.0之间，在选择最优随机共振参数a*和b*的时 候，可W选择将0.01~1.0之间的数字带入上述方程（1.3)内进行比较，得至Ija滿b*。步骤6， 采用小波分析多尺度分解方法对步骤5随机共振处理之后的数据f 21进行多尺度分解，依次 把f2i中的各种频率成分从高向低逐步分离为不同的频带系数fj-乂 乂 W-〉···一〉fl; j代表频带的个数
[0061] 步骤7,在目标信号所在的主要频带区间[fa-fb]选取不同的阔值，fa代表该频带 的最小频率，fb代表该频带的最大频率，对小于阔值的系数置零，对大于阔值的系数进行非 线性处理，而后利用小波重构方法重构出经过阔值处理的信号，重构后的数据记为f3i;
[0062] 步骤8,对步骤7得到的数据f3i进行特征提取，得到信号波动周期、信号衰减长度、 信号峰值等信息；
[0063] 步骤9,把步骤8提取的特征参数，与多种类型水下目标运动的数据库进行匹配比 对，从而识别出最相似的目标类型，实现对目标的分类和特征识别;水下目标运动的数据库 是现有的数据库，是经过多年计算和监测得到的参考数据库。
[0064] 步骤10，在数据处理终端对探测到的目标进行预警提示，告知目标的位置、深度、 大小、类型等参数。
[0065] 通常，噪声被认为是"讨厌"的东西，认为它破坏系统的有序行为，是微弱信号检测 的最大障碍。本发明采用了随机共振技术，解决了传统线性分析方法不能解决的微弱信号 检测问题，能够检测出与噪声信号同频的信号;随机共振采用四阶的龙格-库塔方法数据求 解，计算速度快，能够实时检测分析。基于随机共振数据的小波分析技术，可W实现对目标 扰动信号的窄带检测。基于小波分析技术提取出的目标特征，与水下运动目标特征数据库 对比，就可W得到水下动目标的特征参数识别。
[0066] 本发明基于随机共振和小波分析技术，有效提取出淹没在海洋电磁干扰中的动目 标信号，并与特征数据库进行匹配比对，能够实现水下运动目标的间接探测和预警;利用水 下目标对海水的扰动产生的间接信号，可W间接探测到电磁场隐身的水下目标;利用目标 引起的洋流扰动与海洋本身洋流的非线性作用特性，采用随机共振技术，巧妙利用系统非 线性特性，将海洋背景电磁噪声能量向信号能量的转化，实现海洋背景电磁噪声的"变废为 宝"。
[0067] 实施例：如图2(a)所示，我们测试真实海洋环境中磁场收集天线检测到的数据，包 含海洋扰动、地磁场变化、工频信号、仪器高频噪声和一个小型动目标经过产生的扰动，对 应步骤（1)中的foi;图2(b)为去除仪器高频噪声及工频干扰后的信号，对应步骤（2)中的 hi。
[0068] 如图3(a)、3(b)、3(c)和3(d)所示，显示了步骤(5)中a*和b*的计算方法，如图3(a) 取曰=0.01，6 = 0.02，图3化)取曰=0.02,6 = 0.05、图3(。)取曰=0.01，6 = 0.02，图3((1)取曰= 0.02，b = 0.01，通过上述四组典型数据和图的比较，可W知道在运个动目标中，a = 0.02，b = 0.01为最优的结果，图3(d)中显示的结果对应步骤(5)中的f2i。
[0069] 之后再利用多尺度小波重构及阔值处理方法得到步骤(7)中的信号f3i，如图4所 示。将信号f3i与数据库中的参考数据库进行匹配，得到通过目标的速度约为0.6米/秒，深度 约7米，识别出的目标类型为小型慢速目标，与实验测得的结果相符合。
[0070] W上所述仅为本发明的较佳实施例而已，并不用W限制本发明，凡在本发明的精 神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种洋流干扰中微弱目标信号提取方法，其特征在于是按照下述方式实现的： 步骤1，通过置于水底的磁场收集天线获取海洋磁场变化信息，采集海浪波动、洋流流 动、地球磁场变化和动目标经过产生的磁场扰动信息，采集到的数据表示为fQl，其中i=1… N，N为数据点个数，采样频率为fs Hz; 步骤2,采用小波分析技术对步骤1采集到的数据进行多尺度分解，去除采集仪器本身 高频噪声及环境工频干扰噪声之后得到信号f η ; 步骤3,对步骤2得到的信号fu采取随机共振方法，利用随机共振系统的非线性特性，将 海洋背景电磁噪声的部分能量转化为水下动目标扰动信号的能量，削弱噪声的同时，强化 待测信号，提高动目标信号的信噪比和检出率； 步骤4,对步骤3中的非线性随机共振系统，利用四阶龙格-库塔方法进行高精度数值化 求解，求解步骤如下：式（1.4)-(1.8)中&和?"分别表示共振输出奴〇和输入信号的第11个采样值，？11 =心11，11 = l/fs为时间步长，即数据采样频率匕的倒数;a，b为随机共振系统的系统结构参数;a，b的 取值范围与检测信号中的噪声类型以及信号强度相关，取值范围大于0小于等于1; 步骤5,利用步骤4的方法，通过调节随机共振系统的系统结构参数a，b，找到最优的随 机共振参数a*和b*;基于最优共振参数a*和b*，利用步骤4的方法，求得最佳的共振输出数据 序列，记为f2l; 步骤6,采用小波分析多尺度分解方法对步骤5随机共振处理之后的数据f2l进行多尺度 分解，依次把f2i中的各种频率成分从高向低逐步分离为不同的频带系数fj一〉fj-Ι一〉fj-2- >---->f 1 ； 步骤7,在目标信号所在的主要频带区间[fa+fb]选取不同的阈值，对小于阈值的系数 置零，对大于阈值的系数进行非线性处理，而后利用小波重构方法重构出经过阈值处理的 信号，重构后的数据记为f3l; 步骤8,对步骤7得到的数据f3l进行特征提取，得到信号波动周期、信号衰减长度、信号 峰值信息； 步骤9,把步骤8提取的特征参数，与水下目标运动的数据库进行匹配比对，从而识别出 最相似的目标类型，实现对目标的分类和特征识别； 步骤10，在数据处理终端对探测到的目标进行预警提示，告知目标的位置、深度、大小、 类型参数。2. 根据权利要求1所述的洋流干扰中微弱目标信号提取方法，其特征在于:步骤3采取 的随机共振系统是由郎之万方程描述：把式(1.2)带入式(1.1)，得到式（1.3)中X为系统输出，s(t)为目标扰动引起的周期驱动信号，γ (t)为噪声信号，a，b 为随机共振系统的系统结构参数。3. 根据权利要求1所述的洋流干扰中微弱目标信号提取方法，其特征在于:步骤2中去 除采集仪器本身高频噪声及环境工频干扰噪声的方式是将高于50Hz的信号分量设为零，剩 余信号分量之和记为f η。
【文档编号】G01D3/032GK105823492SQ201610157984
【公开日】2016年8月3日
【申请日】2016年3月18日
【发明人】易忠, 张绍华, 孟立飞, 王斌, 唐小金, 代佳龙, 张超
【申请人】北京卫星环境工程研究所