CN102621533B - 一种用于船用导航雷达的弱小目标检测方法 - Google Patents

一种用于船用导航雷达的弱小目标检测方法 Download PDF

Info

Publication number
CN102621533B
CN102621533B CN201210098333.3A CN201210098333A CN102621533B CN 102621533 B CN102621533 B CN 102621533B CN 201210098333 A CN201210098333 A CN 201210098333A CN 102621533 B CN102621533 B CN 102621533B
Authority
CN
China
Prior art keywords
layer
error
extra large
delta
large clutter
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.)
Expired - Fee Related
Application number
CN201210098333.3A
Other languages
English (en)
Other versions
CN102621533A (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 CN201210098333.3A priority Critical patent/CN102621533B/zh
Publication of CN102621533A publication Critical patent/CN102621533A/zh
Application granted granted Critical
Publication of CN102621533B publication Critical patent/CN102621533B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种用于船用导航雷达的弱小目标检测方法,具体包括训练过程和检测过程,训练过程包括海杂波数据相空间重构和灰色神经网络训练,检测过程具体为海杂波目标检测。本发明的方法基于海杂波的回波内在特性和目标的内在特性不同,利用不含目标的海杂波数据训练灰色神经网络,在纯海杂波情况下总误差均方根值或对消误差趋近零;含有目标的总误差均方根值和对消误差就会很大,因而可以进行弱小目标检测的。本发明的方法相对传统的恒虚警目标检测法,可以在强海杂波背景下检测出弱小目标;相比RBF神经网络进行检测的方法,训练速度相对较快,需要的样本数据信息较少,在海杂波背景下有良好的检测弱小目标的性能。

Description

一种用于船用导航雷达的弱小目标检测方法
技术领域
本发明属于船用导航雷达技术领域,具体涉及其中的弱小目标检测方法。
背景技术
船用导航雷达是船舶不可缺少的导航设备之一,但是船用雷达的目标检测性能常常会受到海杂波的影响。海杂波即为来自海洋表面的雷达反射回波,是船舰导航雷达的主要杂波源。通常情况下它不像地杂波是一种分布散射现象,而是表现出更强的动态特性,因此海杂波的存在严重的影响了雷达对海上目标的检测性能。而随着近些年海运事业、港口贸易经济的发展,船舶对导航雷达的目标检测的速度、精度和抗干扰能力的要求越来越高,以用来能够准确、及时地发现海面目标,减小事故发生几率。所以研究海杂波的特性,把海杂波抑制掉对船舶的安全航行及民用海上救生等方面的应用是很有意义的。
目前处理海杂波主要有两种方法:一个是以传统力学为基础,建立的海杂波分布模型,如瑞利(Rayleigh)分布、对数正态(Log-Normal)分布、韦布尔(Weibul)分布和K分布;另一个方法是基于海杂波混沌动力学特性,建立海杂波模型。
然而经过对海杂波的深入研究,发现海杂波并非完全是一种随机信号,包含了很大程度上的确定性因素。基于传统统计学的建模方法,仅仅只是模拟海杂波的表面特性,无法模拟海杂波的内在特性,而且当海杂波幅度较大,目标较小的时候,如(1)掠海飞行的小型飞机、导弹;(2)小型舰船、导航浮标、小片的冰,漂浮在海面的油污,这些可能会对导航带来潜在的危机等等。基于统计特性建立的模型,在强海杂波和低虚警概率条件下发现弱小目标的能力有限,因此如果能够在一定程度上消除或减少海杂波的干扰,将有可能在很大程度上提高海防预警雷达的目标检测性能。
传统的海杂波背景下目标检测方法,主要是基于统计特性利用其表面看似随机过程而建立起来的分布模型上进行目标检测的,认为海杂波是线性的,目标的回波远大于海杂波回波,设定检测门限,因此这种方法只能检测出较大的目标。然而假如目标较小,海杂波较强的情况下;恒虚警目标检测方法是束手无策的,很难检测出弱小目标,甚至导致漏检的发生。
目前基于混沌特性的目标检测方法主要是利用径向基函数神经网络(Radical BasisFunction)进行检测,虽然此方法确实可行,但也存在着“维数灾难”这个缺点,而恰恰网络收敛的速率对输入空间维数的依赖正是这个“维数灾难”。在RBF神经网络中,逼近函数的固有复杂度与输入空间维数成指数增长关系,根本原因是定义在高维空间的函数很可能远远比低维空间上的函数复杂得多,而且不易区分这些复杂的东西,另外RBF网络的训练耗时也是相对很大。
发明内容
本发明的目的是为了解决现有的船用导航雷达的弱小目标检测中存在的上述问题,提出了一种用于船用导航雷达的弱小目标检测方法。
本发明的技术方案是:一种用于船用导航雷达的弱小目标检测方法,包括训练过程和检测过程,具体的,所述训练过程包括如下步骤:
利用海杂波训练数据进行海杂波相空间的重构,具体包括如下分步骤:
S11:计算嵌入维数m;
S12:计算延迟时间τ;
S13:根据延迟时间τ和嵌入维数m,把海杂波训练数据从时间序列构造成奇异吸引子轨迹向量Xj=(xj,xj+τ,…,xj+(m-1)τ),j=1,…,n,其中,xj表示海杂波训练数据中第j个采样点的海杂波数据,n根据海杂波训练数据总的数目和嵌入维数m得到;
利用重构得到的奇异吸引子轨迹向量对灰色神经网络进行训练,所述灰色神经网络具体包括如下参数,输入参数序号t;网络输入参数y2(t),…,yn(t);网络权值w11,w21,w22,…,w2n,w31,w32,…,w3n;网络预测值y1;灰色神经网络的四层:第一层LA、第二层LB、第三层LC、第四层LD,具体包括如下分步骤:
S21.根据灰色神经网络特征初始化网络结构;
S22.根据网络结构调节权值w11,w21,w22,…,w2n,w31,w32,…,w3n
S23.对每一个输入序列(t,St),t=1,2,3,…,n,计算每层输出:
LA层:a=w11t;
LB层: b = f ( w 11 t ) = 1 1 + e - w 11 t ;
LC层:c1=bw21,c2=y2(t)bw22,c3=y3(t)bw23,…,cn=yn(t)bw2n,y2(t),…,yn(t)与St,t=2,3,…,n一一对应;
LD层:d=w31c1+w32c2+…+w3ncn-θ,θ表示第四层LD网络的预测阈值;
S24.计算网络预测输出与期望输出的误差,并根据误差调整权值和阈值:
LD层误差:δ=d-y1
LC层误差: δ 1 = δ * ( 1 + e - w 11 t ) , δ 2 = δ * ( 1 + e - w 11 t ) , · · · , δ n = δ * ( 1 + e - w 11 t ) ;
LB层误差: δ n + 1 = 1 1 + e - w 11 t * ( 1 - 1 + e - w 11 t ) * ( w 21 δ 1 + w 22 δ 2 + · · · + w 2 n δ n ) ;
调整LB到LC的权值,w21=-y1(0),w22=w221δ2b,…,w2n=w2nn-1δnb,μi表示第二层与第三层的第i+1个节点之间权值调节系数;
调整LA层到LB层的权值:w11=w11+atδn+1
调整阈值: θ = ( 1 + e - w 11 t ) ( w 22 2 y 2 ( t ) + w 23 2 y 3 ( t ) + · · · + w 2 n 2 y n ( t ) - y 1 ( 0 ) ) ;
S25.判断LD层误差δ是否大于预设的期望误差值,若大于,则返回S23,继续训练,否则训练结束;
所述检测过程具体包括如下步骤:
将待检测海杂波数据输入到已经训练好的灰色神经网络中,得到第四层的输出,即是预测值
Figure GDA0000423335530000033
计算对消误差或总误差均方根。
对消误差计算公式如下:
Figure GDA0000423335530000034
其中,xn表示待检测海杂波数据的实际值;
总误差均方根计算公式如下: R = 1 n Σ i = 1 n ( x ‾ n - x n ) 2 ,
若对消误差或总误差均方根趋近零,则表示待检测海杂波数据不含弱小目标;否则,则表示待检测海杂波数据含有弱小目标。
本发明的有益效果:本发明的方法基于海杂波的回波内在特性和目标的内在特性不同,利用不含目标的海杂波数据训练灰色神经网络,在纯海杂波情况下总误差均方根值或对消误差趋近零;含有目标的总误差均方根值和对消误差就会很大,因而可以进行弱小目标检测的。本发明的方法相对传统的恒虚警目标检测法,可以在强海杂波背景下检测出弱小目标;相比RBF神经网络进行检测的方法,避免了RBF具有“维数灾难”的问题,训练速度相对较快,需要的样本数据信息较少,在海杂波背景下有良好的检测弱小目标的性能。
附图说明
图1为本发明弱小目标检测方法的流程示意图。
图2为灰色神经网络的结构示意图。
图3为实施例所用的第#26组海杂波数据示意图。
图4为纯海杂波数据中不含小目标的示意图,
图5对海杂波数据中在采样点数为60000时含一个小球目标的示意图。
图6对图4所示的数据经过本发明方法的处理得到的对消误差示意图。
图7为对图5所示的数据经过本发明方法的处理得到的对消误差示意图。
图8为对图4所示的数据经过本发明方法的处理得到的总均方误差根示意图
图9为对图5所示的数据经过本发明方法的处理得到的总均方误差根示意图。
具体实施方式
下面结合附图和具体实施方式做进一步的说明。
如图1所示,本发明的用于船用导航雷达的目标检测方法,包括训练过程和检测过程,具体的,所述训练过程包括如下步骤:
利用海杂波训练数据进行海杂波相空间的重构,具体包括如下分步骤:
S11:计算嵌入维数m;
S12:计算延迟时间τ;
S13:根据延迟时间τ和嵌入维数m,把海杂波训练数据从时间序列构造成奇异吸引子轨迹向量Xj=(xj,xj+τ,…,xj+(m-1)τ),j=1,…,n,其中,xj表示海杂波训练数据中第j个采样点的海杂波数据,n根据海杂波训练数据总的数目和嵌入维数m得到;
这里,计算嵌入维数m和计算延迟时间τ属于本领域的现有知识,不再详细说明。
利用重构得到的奇异吸引子轨迹向量对灰色神经网络进行训练,灰色神经网络的结构示意图如图2所示,所述灰色神经网络具体包括如下参数,输入参数序号t;网络输入参数y2(t),…,yn(t);网络权值w21,w22,…,w2n,w31,w32,…,w3n;网络预测值y1;灰色神经网络的四层:第一层LA、第二层LB、第三层LC、第四层LD,具体包括如下分步骤:
S21.根据灰色神经网络特征初始化网络结构;
S22.根据网络结构调节权值w11,w21,w22,…,w2n,w31,w32,…,w3n
S23.对每一个输入序列(t,St),t=1,2,3,…,n,计算每层输出:
LA层:a=w11t;
LB层: b = f ( w 11 t ) = 1 1 + e - w 11 t ;
LC层:c1=bw21,c2=y2(t)bw22,c3=y3(t)bw23,…,cn=yn(t)bw2n,y2(t),…,yn(t)与St,t=2,3,…,n一一对应;
LD层:d=w31c1+w32c2+…+w3ncn-θ,θ表示第四层LD网络的预测阈值;
S24.计算网络预测输出与期望输出的误差,并根据误差调整权值和阈值:
LD层误差:δ=d-y1
LC层误差: δ 1 = δ * ( 1 + e - w 11 t ) , δ 2 = δ * ( 1 + e - w 11 t ) , · · · , δ n = δ * ( 1 + e - w 11 t ) ;
LB层误差: δ n + 1 = 1 1 + e - w 11 t * ( 1 - 1 + e - w 11 t ) * ( w 21 δ 1 + w 22 δ 2 + · · · + w 2 n δ n ) ;
调整LB到LC的权值,w21=-y1(0),w22=w221δ2b,…,w2n=w2nn-1δnb,μi表示第二层与第三层的第i+1个节点之间权值调节系数;
调整LA层到LB层的权值:w11=w11+atδn+1
调整阈值: θ = ( 1 + e - w 11 t ) ( w 22 2 y 2 ( t ) + w 23 2 y 3 ( t ) + · · · + w 2 n 2 y n ( t ) - y 1 ( 0 ) ) ;
S25.判断LD层误差δ是否大于预设的期望误差值,若大于,则返回S23,继续训练,否则训练结束;
所述检测过程具体包括如下步骤:
将待检测海杂波数据输入到已经训练好的灰色神经网络中,得到第四层的输出,即是预测值
Figure GDA0000423335530000054
计算对消误差。
对消误差计算公式如下:
Figure GDA0000423335530000055
其中,xn表示待检测海杂波数据的实际值;
若对消误差趋近零,则表示待检测海杂波数据不含弱小目标;否则,则表示待检测海杂波数据含有弱小目标。
这里的也可以采用总误差均方根进行评价,总误差均方根计算公式如下:
Figure GDA0000423335530000056
若总误差均方根趋近零,则表示待检测海杂波数据不含弱小目标;否则,则表示待检测海杂波数据含有弱小目标。
这里的,趋近零可以理解为小于0.01,当海杂波数据含有弱小目标时,检测得到的对消误差或总误差均方根会有较大的很明显的波动,这里对趋近零的理解不影响本发明的实现。
在这里采用加拿大McMaster大学S.Haykin教授采用IPIX雷达在Dartmouth地区海域实测带有目标的海杂波数据对本发明所述的方法进行验证,具体用第#26组海杂波数据,如图3所示,此数据总共有14个距离单元,每个距离15m,每个距离单元有131072个采样点数,此数据在第七距离单元第60000个采样点数附近有一个用铝箔包的小球目标。
图4是回波数据为纯海杂中不含小目标的示意图,图5是回波数据为在采样点数为60000时含一个小球目标。图6为对图4所示的数据经过本发明方法的处理得到的对消误差示意图,图7为对图5所示的数据经过本发明方法的处理得到的对消误差示意图,图8为对图4所示的数据经过本发明方法的处理得到的总均方误差根示意图,图9为对图5所示的数据经过本发明方法的处理得到的总均方误差根示意图。
可以看出,本发明的方法基于海杂波的回波内在特性和目标的内在特性不同,利用不含目标的海杂波数据训练灰色神经网络,在纯海杂波情况下总误差均方根值或对消误差趋近零;含有目标的总误差均方根值和对消误差就会很大,因而可以进行弱小目标检测的。本发明的方法相对传统的恒虚警目标检测法,可以在强海杂波背景下检测出弱小目标;相比RBF神经网络进行检测的方法,避免了RBF具有“维数灾难”的问题,训练速度相对较快,需要的样本数据信息较少,在海杂波背景下有良好的检测弱小目标的性能。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (1)

1.一种用于船用导航雷达的弱小目标检测方法,其特征在于,包括训练过程和检测过程,具体的,所述训练过程包括如下步骤:
利用海杂波训练数据进行海杂波相空间的重构,具体包括如下分步骤:
S11:计算嵌入维数m;
S12:计算延迟时间τ;
S13:根据延迟时间τ和嵌入维数m,把海杂波训练数据从时间序列构造成奇异吸引子轨迹向量Xj=(xj,xj+τ,…,xj+(m-1)τ),j=1,…,n,其中,xj表示海杂波训练数据中第j个采样点的海杂波数据,n根据海杂波训练数据总的数目和嵌入维数m得到;
利用重构得到的奇异吸引子轨迹向量对灰色神经网络进行训练,所述灰色神经网络具体包括如下参数,输入参数序号t;网络输入参数y2(t),…,yn(t);网络权值w11,w21,w22,…,w2n,w31,w32,…,w3n;网络预测值y1;灰色神经网络的四层:第一层LA、第二层LB、第三层LC、第四层LD,具体包括如下分步骤:
S21.根据灰色神经网络特征初始化网络结构;
S22.根据网络结构调节权值w11,w21,w22,…,w2n,w31,w32,…,w3n
S23.对每一个输入序列(t,St),t=1,2,3,…,n,计算每层输出:
LA层:a=w11t;
LB层: b = f ( w 11 t ) = 1 1 + e - w 11 t ;
LC层:c1=bw21,c2=y2(t)bw22,c3=y3(t)bw23,…,cn=yn(t)bw2n,y2(t),…,yn(t)与St,t=2,3,…,n一一对应;
LD层:d=w31c1+w32c2+…+w3ncn-θ,θ表示第四层LD网络的预测阈值;
S24.计算网络预测输出与期望输出的误差,并根据误差调整权值和阈值:
LD层误差:δ=d-y1
LC层误差: δ 1 = δ * ( 1 + e - w 11 t ) , δ 2 = δ * ( 1 + e - w 11 t ) , · · · , δ n = δ * ( 1 + e - w 11 t ) ;
LB层误差: δ n + 1 = 1 1 + e - w 11 t * ( 1 - 1 + e - w 11 t ) * ( w 21 δ 1 + w 22 δ 2 + · · · + w 2 n δ n ) ;
调整LB到LC的权值,w21=-y1(0),w22=w221δ2b,…,w2n=w2nn-1δnb,μi表示第二层与第三层的第i+1个节点之间权值调节系数;
调整LA层到LB层的权值:w11=w11+atδn+1
调整阈值: θ = ( 1 + e - w 11 t ) ( w 22 2 y 2 ( t ) + w 23 2 y 3 ( t ) + · · · + w 2 n 2 y n ( t ) - y 1 ( 0 ) ) ;
S25.判断LD层误差δ是否大于预设的期望误差值,若大于,则返回S23,继续训练,否则训练结束;
所述检测过程具体包括如下步骤:
将待检测海杂波数据输入到已经训练好的灰色神经网络中,得到第四层的输出,即是预测值
Figure FDA0000423335520000022
计算对消误差或总误差均方根;
对消误差计算公式如下:
Figure FDA0000423335520000023
其中,xn表示待检测海杂波数据的实际值;
总误差均方根计算公式如下:
Figure FDA0000423335520000024
若对消误差或总误差均方根趋近零,则表示待检测海杂波数据不含弱小目标;否则,则表示待检测海杂波数据含有弱小目标。
CN201210098333.3A 2012-04-05 2012-04-05 一种用于船用导航雷达的弱小目标检测方法 Expired - Fee Related CN102621533B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210098333.3A CN102621533B (zh) 2012-04-05 2012-04-05 一种用于船用导航雷达的弱小目标检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210098333.3A CN102621533B (zh) 2012-04-05 2012-04-05 一种用于船用导航雷达的弱小目标检测方法

Publications (2)

Publication Number Publication Date
CN102621533A CN102621533A (zh) 2012-08-01
CN102621533B true CN102621533B (zh) 2014-06-18

Family

ID=46561561

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210098333.3A Expired - Fee Related CN102621533B (zh) 2012-04-05 2012-04-05 一种用于船用导航雷达的弱小目标检测方法

Country Status (1)

Country Link
CN (1) CN102621533B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103577877B (zh) * 2013-11-19 2016-05-25 北京航空航天大学 一种基于时频分析和bp神经网络的船舶运动预报方法
CN107238823A (zh) * 2017-06-13 2017-10-10 电子科技大学 真假目标的状态特征提取方法
CN108009469B (zh) * 2017-10-24 2020-11-24 中国科学院电子学研究所苏州研究院 一种基于结构递归神经网络的海上油气平台检测方法
JP7182869B2 (ja) * 2017-12-28 2022-12-05 古野電気株式会社 物標検出装置
CN108388732B (zh) * 2018-02-27 2021-11-05 中国人民解放军空军工程大学 掠海目标多径散射特性仿真计算方法及系统
CN108711312B (zh) * 2018-05-24 2020-09-01 大连海事大学 基于bp神经网络的船舶与静态物标碰撞风险预判方法
CN110263646A (zh) * 2019-05-21 2019-09-20 华中科技大学 一种基于卷积神经网络的海面弱目标探测方法和系统
CN112101250B (zh) * 2020-09-18 2022-07-15 电子科技大学 一种基于上下文语义感知的近岸舰船目标检测方法
CN112068085B (zh) * 2020-10-16 2022-05-06 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于深度学习的雷达海杂波原始数据快速预处理方法
CN112711026B (zh) * 2020-12-19 2022-09-13 西北工业大学 基于变微分步长相对变化熵的目标检测与信号频率估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101887116A (zh) * 2010-04-14 2010-11-17 电子科技大学 一种用于微弱目标的改进型随机Hough变换检测方法

Also Published As

Publication number Publication date
CN102621533A (zh) 2012-08-01

Similar Documents

Publication Publication Date Title
CN102621533B (zh) 一种用于船用导航雷达的弱小目标检测方法
Lipa et al. HF radar detection of tsunamis
CN107607092B (zh) 一种基于无人机遥感的波浪与浮体运动预报系统
Al-Ani et al. Deterministic sea waves prediction using mixed space–time wave radar data
Gopalakrishnan et al. Assimilation of HF radar-derived surface currents on tidal-timescales
CN108051784A (zh) 基于先验知识的静止杂波环境动目标检测的波形优化方法
CN103308910B (zh) 利用高频地波雷达检测海上非航行状态船只目标的方法
CN103809170A (zh) 一种多普勒频偏修正方法及系统
CN103487793B (zh) 一种基于简正波理论的宽带混响波形仿真方法
CN105589061A (zh) 一种岸基雷达的信号处理算法
Mas-Soler et al. Estimating on-site wave spectra from the motions of a semi-submersible platform: an assessment based on model scale results
Brushett et al. Evaluation of met-ocean forecast data effectiveness for tracking drifters deployed during operational oil spill response in Australian waters
Naaijen et al. Reducing operational risks by on-board phase resolved prediction of wave induced ship motions
Smith et al. Methods of associating CODAR SeaSonde vessel detection data into unique tracks
Piscopo et al. Measurement of sea waves
Li et al. Small boat detection for radar image datasets with yolo V3 network
Boltryk et al. An ultrasonic transducer array for velocity measurement in underwater vehicles
Vesecky et al. Identifying ship echoes in CODAR HF radar data: A Kalman filtering approach
Rätsep et al. Surface vessel localization from wake measurements using an array of pressure sensors in the littoral zone
Laws et al. Monitoring coastal vessels for environmental applications: Application of Kalman filtering
Caiazzo et al. The limitations of statistical low runs prediction in rough seas: A study based on real wave data
Radhakrishnan et al. A Methodology for Tuning of Computational Vessel Models Utilizing Wave Measurements From X-Band Marine Radar and Wave Buoy
Zhang et al. Prediction method of ocean wave spectrum based on an echo state network
Wang et al. Pre-processing of images from shipborne X-band wave measuring radar
Панасенкова et al. On Iceberg Drift in the Barents Sea

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140618

Termination date: 20170405

CF01 Termination of patent right due to non-payment of annual fee