CN107942299A - 一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法 - Google Patents

一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法 Download PDF

Info

Publication number
CN107942299A
CN107942299A CN201711115256.7A CN201711115256A CN107942299A CN 107942299 A CN107942299 A CN 107942299A CN 201711115256 A CN201711115256 A CN 201711115256A CN 107942299 A CN107942299 A CN 107942299A
Authority
CN
China
Prior art keywords
frog
radar
value
sea clutter
group
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.)
Pending
Application number
CN201711115256.7A
Other languages
English (en)
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201711115256.7A priority Critical patent/CN107942299A/zh
Publication of CN107942299A publication Critical patent/CN107942299A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Computing Systems (AREA)
  • Mathematical Analysis (AREA)
  • Software Systems (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Evolutionary Computation (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法，系统由雷达、数据库和上位机依次相连组成，雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库，所述的上位机包括数据预处理模块、鲁棒预报模型建模模块、智能寻优模块、目标检测模块、模型更新模块以及结果显示模块：本发明针对海上目标检测的复杂特性，对雷达杂波数据进行重构，并对目标数据进行检测，引入改进混合蛙跳算法，从而提供一种实现在线检测、智能性高的雷达海上目标检测系统及方法。

Description

一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及 方法
技术领域
本发明涉及雷达数据处理领域，特别地，涉及一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法。
背景技术
海杂波,即来自于海面的雷达后向散射回波。近几十年来，随着对海杂波认识的深入，德国、挪威等国家相继尝试利用雷达观测海杂波获取雷达海浪图像来反演海浪信息，以获得关于海洋状态的实时信息，如海浪的波高、方向和周期等，从而进一步对海上微小目标进行检测，这对海上活动具有十分重要的意义。
海上目标检测技术具有重要的地位，提供准确的目标判决是对海雷达工作的重要任务之一。雷达自动检测系统依据判决准则在给定的检测阈值下做出判决，而强海杂波往往成为微弱目标信号的主要干扰。如何处理海杂波将直接影响到雷达在海洋环境下的检测能力：1)识别导航浮标、小片的冰，漂浮在海面的油污，这些可能会对导航带来潜在的危机；3)监测非法捕鱼是环境监测的一项重要的任务。
在传统的目标检测时,海杂波被认为是干扰导航的一种噪声被去掉。然而，在雷达对海观测目标时，微弱的运动目标回波常常湮没在海杂波中，信杂比较低，雷达不易检测到目标，同时海杂波的大量尖峰还会造成严重虚警，对雷达的检测性能产生较大影响。对于各种对海警戒和预警雷达而言，研究的主要目标是提高海杂波背景下目标的检测能力。因此，不仅具有重要的理论意义和实际意义，而且也是国内外海上目标检测的难点和热点。
发明内容
为了克服已有雷达海上目标检测方法无法实现在线检测、智能性较差的不足，本发明提供一种实现在线检测、智能性强的基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法。
本发明解决其技术问题所采用的技术方案是：
一种基于改进混合蛙跳算法的智能雷达海上目标检测系统，包括雷达、数据库以及上位机，雷达、数据库和上位机依次相连，所述雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库，所述的上位机包括：
数据预处理模块，用以进行雷达海杂波数据预处理，采用如下过程完成：
(1)雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库；
(2)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本，i＝1,...,N；
(3)对训练样本进行归一化处理，得到归一化幅值
其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；
(4)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：
其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；
所述鲁棒预报模型建模模块，用以建立预报模型，采用如下过程完成：
将数据预处理模块得到的X、Y代入如下线性方程：
其中
权重因子vi由下式计算：
其中是误差变量ξi标准差的估计，c1,c2为常量；
求解得待估计函数f(x)：
其中，M是支持向量的数目，1v＝[1,...,1]T,上标T表示矩阵的转置，是拉格朗日乘子，b*是偏置量，K＝exp(-||xi-xj||/θ2)，其中i＝1,…,M，j＝1,…,M，和exp(-||x-xi||/θ2)均为支持向量机的核函数，xj为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；
智能寻优模块，采用改进混合蛙跳算法对鲁棒预报模型的核参数θ和惩罚系数γ进行优化，采用如下过程完成：
步骤1：初始化青蛙群体参数，设种群青蛙数为P,最大迭代次数Maxgen，局部搜索的迭代次数Mmax，最大更新长度Dmax，分组数m以及每组青蛙数n，由于模型有两个参数需要优化，所以位置pi的维度为2维，随机生成每个青蛙的位置pi＝(pi1,pi2)，置初始迭代次数k＝0；
步骤2：计算所有青蛙的适应度值，并进行排序，分组，选取种群最优青蛙pg；
步骤3：按照下式对子群最差青蛙对子群中最差的青蛙进行更新，接着在子群内重新排序，然后对子群最差青蛙进行更新；重复该局部搜索过程Mmax次；
D(k+1)＝W×(r1×D(k)+r2×D(k-1)+r3×His(pw))+r4×(pb-pw)
p'w(k+1)＝pw(k)+D,-Dmax≤D≤Dmax
W＝We+(Ws-We)×a
a＝exp(-30×(k/Maxgen)s)
其中，pw为子群最差青蛙，pb为子群最优青蛙，Dmax为最大变异尺度，p'w为更新后的青蛙，r1、r2、r3、r4是0～1之间的随机数，k是当前迭代次数，W是权重因子，Ws、We是权重因子的起始值和最终值，s为大于1的整数。首先利用子群最优青蛙去更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则种群最优青蛙代替子群最优青蛙进行更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则随机生成一个青蛙取代原来子群最差青蛙。
步骤4：当所有子群的局部搜索都已经完成，对所有青蛙进行混合、排序和分组，选取种群最优青蛙pg；
步骤5：k＝k+1，若k<Maxgen，则转至步骤3；否则输出种群最优青蛙xg，为鲁棒预报模型的最优参数，算法终止；
其中，初始种群大小为200，分组数10，每组子群数20，种群最大迭代数100，子群最大迭代数10，最大更新长度5，权重因子的起始值和最终值分别为0.9和0.4，s为3。
目标检测模块，用以进行目标检测，采用如下过程完成：
1)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[xt-D+1,…,xt]，xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，xt表示第t采样时刻的海杂波回波信号幅值；
2)进行归一化处理；
3)代入鲁棒预报模型建模模块得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。
4)计算海杂波预报值与雷达回波实测值的差值e，计算控制限Qα：
其中，α是置信度，θ1,θ2,θ3,h0是中间变量，λj i表示协方差矩阵的第j个特征值的i次方，k是样本维数，Cα是正态分布置信度为α的统计；
5)进行检测判断：当e2差值大于控制限Qα时，该点存在目标，否则没有目标。
模型更新模块，用以按设定的采样时间间隔，采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。
结果显示模块，用以将目标检测模块的检测结果在上位机显示。
一种基于改进混合蛙跳算法的智能雷达海上目标检测系统所使用的雷达海上目标检测方法，所述的方法包括以下步骤：
(1)雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库；
(2)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本，i＝1,...,N；
(3)对训练样本进行归一化处理，得到归一化幅值xi：
其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；
(4)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：
其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；
(5)将得到的X、Y代入如下线性方程：
其中
权重因子vi由下式计算：
其中是误差变量ξi标准差的估计，c1,c2为常量；
求解得待估计函数f(x)：
其中，M是支持向量的数目，1v＝[1,...,1]T,上标T表示矩阵的转置，是拉格朗日乘子，b*是偏置量，K＝exp(-||xi-xj||/θ2)，其中i＝1,…,M，j＝1,…,M，和exp(-||x-xi||/θ2)均为支持向量机的核函数，xj为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；
(6)用混合蛙跳算法对步骤4)的核参数θ和惩罚系数γ进行优化，采用如下过程完成：
(6.1)初始化青蛙群体参数，设种群青蛙数为P,最大迭代次数Maxgen，局部搜索的迭代次数Mmax，最大更新长度Dmax，分组数m以及每组青蛙数n，由于模型有两个参数需要优化，所以位置pi的维度为2维，随机生成每个青蛙的位置pi＝(pi1,pi2)，置初始迭代次数k＝0；
(6.2)计算所有青蛙的适应度值，并进行排序，分组，选取种群最优青蛙pg；
(6.3)按照下式对子群最差青蛙对子群中最差的青蛙进行更新，接着在子群内重新排序，然后对子群最差青蛙进行更新；重复该局部搜索过程Mmax次；
D(k+1)＝W×(r1×D(k)+r2×D(k-1)+r3×His(pw))+r4×(pb-pw)
p'w(k+1)＝pw(k)+D,-Dmax≤D≤Dmax
W＝We+(Ws-We)×a
a＝exp(-30×(k/Maxgen)s)
其中，pw为子群最差青蛙，pb为子群最优青蛙，Dmax为最大变异尺度，p'w为更新后的青蛙，r1、r2、r3、r4是0～1之间的随机数，k是当前迭代次数，W是权重因子，Ws、We是权重因子的起始值和最终值，s为大于1的整数。首先利用子群最优青蛙去更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则种群最优青蛙代替子群最优青蛙进行更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则随机生成一个青蛙取代原来子群最差青蛙。
(6.4)当所有子群的局部搜索都已经完成，对所有青蛙进行混合、排序和分组，选取种群最优青蛙pg；
(6.5)k＝k+1，若k<Maxgen，则转至步骤3；否则输出种群最优青蛙xg，为鲁棒预报模型的最优参数，算法终止；
其中，初始种群大小为200，分组数10，每组子群数20，种群最大迭代数100，子群最大迭代数10，最大更新长度5，权重因子的起始值和最终值分别为0.9和0.4，s为3。
(7)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[xt-D+1,…,xt]，TX表示海杂波从第t-D+1采样时刻到第t采样时刻的信号幅值矩阵，xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，xt表示第t采样时刻的海杂波回波信号幅值；
(8)进行归一化处理；
(9)代入步骤(5)得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。
(10)计算海杂波预报值与雷达回波实测值的差值e，计算控制限Qα：
其中，α是置信度，θ1,θ2,θ3,h0是中间变量，λj i表示协方差矩阵的第j个特征值的i次方，k是样本维数，Cα是正态分布置信度为α的统计；
(11)进行检测判断：当e2差值大于控制限Qα时，该点存在目标，否则没有目标。
(12)按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。
本发明的技术构思为：本发明针对雷达海杂波的混沌特性，对雷达海杂波数据进行重构，并对重构后的数据进行非线性拟合，建立雷达海杂波的预报模型，计算雷达海杂波的预报值和实测值的差，有目标存在时的误差会显著大于没有目标时，引入改进混合蛙跳算法，从而实现海杂波背景下的强智能目标检测。
本发明的有益效果主要表现在：1、可在线检测海上目标；2、所用的检测方法只需较少样本；3、智能性强、受人为因素影响小。
附图说明
图1是本发明所提出的系统的硬件结构图；
图2是本发明所提出的上位机的功能模块图。
具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。
实施例1
参照图1、图2，一种基于改进混合蛙跳算法的智能雷达海上目标检测系统，包括雷达1、数据库2、及上位机3，雷达1、数据库2和上位机3依次相连，所述雷达1对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库2，所述的上位机3包括：
数据预处理模块4，用以进行雷达海杂波数据预处理，采用如下过程完成：
(1)雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库；
(2)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本，i＝1,...,N；
(3)对训练样本进行归一化处理，得到归一化幅值
其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；
(4)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：
其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；
所述鲁棒预报模型建模模块，用以建立预报模型，采用如下过程完成：
将数据预处理模块得到的X、Y代入如下线性方程：
其中
权重因子vi由下式计算：
其中是误差变量ξi标准差的估计，c1,c2为常量；
求解得待估计函数f(x)：
其中，M是支持向量的数目，1v＝[1,...,1]T,上标T表示矩阵的转置，是拉格朗日乘子，b*是偏置量，K＝exp(-||xi-xj||/θ2)，其中i＝1,…,M，j＝1,…,M，和exp(-||x-xi||/θ2)均为支持向量机的核函数，xj为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；
智能寻优模块6，用以采用改进混合蛙跳算法对鲁棒预报模型的核参数θ和惩罚系数γ进行优化，采用如下过程完成：
步骤1：初始化青蛙群体参数，设种群青蛙数为P,最大迭代次数Maxgen，局部搜索的迭代次数Mmax，最大更新长度Dmax，分组数m以及每组青蛙数n，由于模型有两个参数需要优化，所以位置pi的维度为2维，随机生成每个青蛙的位置pi＝(pi1,pi2)，置初始迭代次数k＝0；
步骤2：计算所有青蛙的适应度值，并进行排序，分组，选取种群最优青蛙pg；
步骤3：按照下式对子群最差青蛙对子群中最差的青蛙进行更新，接着在子群内重新排序，然后对子群最差青蛙进行更新；重复该局部搜索过程Mmax次；
D(k+1)＝W×(r1×D(k)+r2×D(k-1)+r3×His(pw))+r4×(pb-pw)
p'w(k+1)＝pw(k)+D,-Dmax≤D≤Dmax
W＝We+(Ws-We)×a
a＝exp(-30×(k/Maxgen)s)
其中，pw为子群最差青蛙，pb为子群最优青蛙，Dmax为最大变异尺度，p'w为更新后的青蛙，r1、r2、r3、r4是0～1之间的随机数，k是当前迭代次数，W是权重因子，Ws、We是权重因子的起始值和最终值，s为大于1的整数。首先利用子群最优青蛙去更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则种群最优青蛙代替子群最优青蛙进行更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则随机生成一个青蛙取代原来子群最差青蛙。
步骤4：当所有子群的局部搜索都已经完成，对所有青蛙进行混合、排序和分组，选取种群最优青蛙pg；
步骤5：k＝k+1，若k<Maxgen，则转至步骤3；否则输出种群最优青蛙xg，为鲁棒预报模型的最优参数，算法终止；
其中，初始种群大小为200，分组数10，每组子群数20，种群最大迭代数100，子群最大迭代数10，最大更新长度5，权重因子的起始值和最终值分别为0.9和0.4，s为3。
目标检测模块7，用以进行目标检测，采用如下过程完成：
1)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[xt-D+1,…,xt]，xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，xt表示第t采样时刻的海杂波回波信号幅值；
2)进行归一化处理；
3)代入鲁棒预报模型建模模块得到的函数f(x)得到采样时刻(t+1)的海杂波预报值；
4)计算海杂波预报值与雷达回波实测值的差值e，计算控制限Qα：
其中，α是置信度，θ1,θ2,θ3,h0是中间变量，λj i表示协方差矩阵的第j个特征值的i次方，k是样本维数，Cα是正态分布置信度为α的统计；
5)进行检测判断：当e2差值大于控制限Qα时，该点存在目标，否则没有目标。
模型更新模块8，用以按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。
结果显示模块9，用以将目标检测模块的检测结果在上位机显示。
所述上位机3的硬件部分包括：I/O元件，用于数据的采集和信息的传递；数据存储器，存储运行所需的数据样本和运行参数等；程序存储器，存储实现功能模块的软件程序；运算器，执行程序，实现指定的功能；显示模块，显示设置的参数和检测结果。
实施例2
参照图1、图2，一种基于改进混合蛙跳算法的智能雷达海上目标检测方法，所述的方法包括以下步骤：
(1)雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库；
(2)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本，i＝1,...,N；
(3)对训练样本进行归一化处理，得到归一化幅值
其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；
(4)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：
其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；
(5)将得到的X、Y代入如下线性方程：
其中
权重因子vi由下式计算：
其中是误差变量ξi标准差的估计，c1,c2为常量；
求解得待估计函数f(x)：
其中，M是支持向量的数目，1v＝[1,...,1]T,上标T表示矩阵的转置，是拉格朗日乘子，b*是偏置量，K＝exp(-||xi-xj||/θ2)，其中i＝1,…,M，j＝1,…,M，和exp(-||x-xi||/θ2)均为支持向量机的核函数，xj为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；
(6)用混合蛙跳算法对步骤4)的核参数θ和惩罚系数γ进行优化，采用如下过程完成：
(6.1)初始化青蛙群体参数，设种群青蛙数为P,最大迭代次数Maxgen，局部搜索的迭代次数Mmax，最大更新长度Dmax，分组数m以及每组青蛙数n，由于模型有两个参数需要优化，所以位置pi的维度为2维，随机生成每个青蛙的位置pi＝(pi1,pi2)，置初始迭代次数k＝0；
(6.2)计算所有青蛙的适应度值，并进行排序，分组，选取种群最优青蛙pg；
(6.3)按照下式对子群最差青蛙对子群中最差的青蛙进行更新，接着在子群内重新排序，然后对子群最差青蛙进行更新；重复该局部搜索过程Mmax次；
D(k+1)＝W×(r1×D(k)+r2×D(k-1)+r3×His(pw))+r4×(pb-pw)
p'w(k+1)＝pw(k)+D,-Dmax≤D≤Dmax
W＝We+(Ws-We)×a
a＝exp(-30×(k/Maxgen)s)
其中，pw为子群最差青蛙，pb为子群最优青蛙，Dmax为最大变异尺度，p'w为更新后的青蛙，r1、r2、r3、r4是0～1之间的随机数，k是当前迭代次数，W是权重因子，Ws、We是权重因子的起始值和最终值，s为大于1的整数。首先利用子群最优青蛙去更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则种群最优青蛙代替子群最优青蛙进行更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则随机生成一个青蛙取代原来子群最差青蛙。
(6.4)当所有子群的局部搜索都已经完成，对所有青蛙进行混合、排序和分组，选取种群最优青蛙pg；
(6.5)k＝k+1，若k<Maxgen，则转至步骤3；否则输出种群最优青蛙xg，为鲁棒预报模型的最优参数，算法终止；
其中，初始种群大小为200，分组数10，每组子群数20，种群最大迭代数100，子群最大迭代数10，最大更新长度5，权重因子的起始值和最终值分别为0.9和0.4，s为3。
(7)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[xt-D+1,…,xt]，TX表示海杂波从第t-D+1采样时刻到第t采样时刻的信号幅值矩阵，xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，xt表示第t采样时刻的海杂波回波信号幅值；
(8)进行归一化处理；
(9)代入步骤(5)得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。
(10)计算海杂波预报值与雷达回波实测值的差值e，计算控制限Qα：
其中，α是置信度，θ1,θ2,θ3,h0是中间变量，λj i表示协方差矩阵的第j个特征值的i次方，k是样本维数，Cα是正态分布置信度为α的统计；
(11)进行检测判断：当e2差值大于控制限Qα时，该点存在目标，否则没有目标。
(12)按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。
由以上实施例可见，本发明建立了智能雷达海上目标检测系统和方法，可以在线检测雷达目标；且所用的检测方法只需较少样本即可；另外，减少了人为因素的影响，智能性高，鲁棒性强。

Claims (2)

1.一种基于改进混合蛙跳算法的智能雷达海上目标检测系统，包括雷达、数据库以及上位机，雷达、数据库和上位机依次相连，其特征在于：所述雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库，所述的上位机包括数据预处理模块、鲁棒预报模型建模模块、智能寻优模块、目标检测模块、模型更新模块以及结果显示模块：
所述数据预处理模块，用以进行雷达海杂波数据预处理，采用如下过程完成：
(1)雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库；
(2)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本，i＝1,...,N；
(3)对训练样本进行归一化处理，得到归一化幅值
<mrow> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> <mrow> <mi>max</mi> <mi> </mi> <mi>x</mi> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> </mfrac> </mrow>
其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；
(4)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：
<mrow> <mi>Y</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>D</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>D</mi> <mo>+</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mi>N</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；
所述鲁棒预报模型建模模块，用以建立预报模型，采用如下过程完成：
将数据预处理模块得到的X、Y代入如下线性方程：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msubsup> <mn>1</mn> <mi>v</mi> <mi>T</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msub> <mn>1</mn> <mi>v</mi> </msub> </mtd> <mtd> <mrow> <mi>K</mi> <mo>+</mo> <msub> <mi>V</mi> <mi>&amp;gamma;</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msup> <mi>b</mi> <mo>*</mo> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mi>&amp;alpha;</mi> <mo>*</mo> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mi>Y</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中
权重因子vi由下式计算：
其中是误差变量ξi标准差的估计，c1,c2为常量；
求解得待估计函数f(x)：
<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msup> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <mo>*</mo> </msup> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mo>|</mo> <mo>|</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>/</mo> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>b</mi> <mo>*</mo> </msup> </mrow>
其中，M是支持向量的数目，1v＝[1,...,1]T,上标T表示矩阵的转置，是拉格朗日乘子，b*是偏置量，K＝exp(-||xi-xj||/θ2)，其中i＝1,…,M，j＝1,…,M，和exp(-||x-xi||/θ2)均为支持向量机的核函数，xj为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；
所述智能寻优模块，用以采用改进混合蛙跳算法对鲁棒预报模型的核参数θ和惩罚系数γ进行优化，采用如下过程完成：
(A)初始化青蛙群体参数，设种群青蛙数为P,最大迭代次数Maxgen，局部搜索的迭代次数Mmax，最大更新长度Dmax，分组数m以及每组青蛙数n，由于模型有两个参数需要优化，所以位置pi的维度为2维，随机生成每个青蛙的位置pi＝(pi1,pi2)，置初始迭代次数k＝0；
(B)计算所有青蛙的适应度值，并进行排序，分组，选取种群最优青蛙pg；
(C)按照下式对子群最差青蛙对子群中最差的青蛙进行更新，接着在子群内重新排序，然后对子群最差青蛙进行更新；重复该局部搜索过程Mmax次；
D(k+1)＝W×(r1×D(k)+r2×D(k-1)+r3×His(pw))+r4×(pb-pw)
p'w(k+1)＝pw(k)+D,-Dmax≤D≤Dmax
W＝We+(Ws-We)×a
a＝exp(-30×(k/Maxgen)s)
其中，pw为子群最差青蛙，pb为子群最优青蛙，Dmax为最大变异尺度，p'w为更新后的青蛙，r1、r2、r3、r4是0～1之间的随机数，k是当前迭代次数，W是权重因子，Ws、We是权重因子的起始值和最终值，s为大于1的整数。首先利用子群最优青蛙去更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则种群最优青蛙代替子群最优青蛙进行更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则随机生成一个青蛙取代原来子群最差青蛙。
(D)当所有子群的局部搜索都已经完成，对所有青蛙进行混合、排序和分组，选取种群最优青蛙pg；
(E)k＝k+1，若k<Maxgen，则转至步骤3；否则输出种群最优青蛙xg，为鲁棒预报模型的最优参数，算法终止；
其中，初始种群大小为200，分组数10，每组子群数20，种群最大迭代数100，子群最大迭代数10，最大更新长度5，权重因子的起始值和最终值分别为0.9和0.4，s为3。
所述目标检测模块，用以进行目标检测，采用如下过程完成：
(a)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[xt-D+1,…,xt]，xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，xt表示第t采样时刻的海杂波回波信号幅值；
(b)进行归一化处理；
<mrow> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mi>T</mi> <mi>X</mi> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> <mrow> <mi>max</mi> <mi> </mi> <mi>x</mi> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> </mfrac> </mrow>
(c)代入鲁棒预报模型建模模块得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。
(d)计算海杂波预报值与雷达回波实测值的差值e，计算控制限Qα：
<mrow> <msub> <mi>Q</mi> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <msup> <mrow> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mi>&amp;alpha;</mi> </msub> <msub> <mi>h</mi> <mn>0</mn> </msub> <msqrt> <mrow> <mn>2</mn> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> </mrow> </msqrt> </mrow> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <msub> <mi>h</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> </mfrac> <mo>&amp;rsqb;</mo> </mrow> <mfrac> <mn>1</mn> <msub> <mi>h</mi> <mn>0</mn> </msub> </mfrac> </msup> </mrow>
<mrow> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>&amp;lambda;</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> </mrow>
<mrow> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> </mrow> <mrow> <mn>3</mn> <msubsup> <mi>&amp;theta;</mi> <mn>2</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow>
其中，α是置信度，θ1,θ2,θ3,h0是中间变量，λj i表示协方差矩阵的第j个特征值的i次方，k是样本维数，Cα是正态分布置信度为α的统计；
(e)进行检测判断：当e2差值大于控制限Qα时，该点存在目标，否则没有目标。
所述模型更新模块，用以按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。
所述结果显示模块，用以将目标检测模块的检测结果在上位机显示。
2.一种权利要求1所述的基于混合蛙跳算法的智能雷达海上目标检测系统所使用的雷达海上目标检测方法，其特征在于，所述的方法包括以下步骤：
(1)雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库；
(2)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本，i＝1,...,N；
(3)对训练样本进行归一化处理，得到归一化幅值
<mrow> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> <mrow> <mi>max</mi> <mi> </mi> <mi>x</mi> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> </mfrac> </mrow>
其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；
(4)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：
<mrow> <mi>Y</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>D</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>D</mi> <mo>+</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mi>N</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；
(5)将得到的X、Y代入如下线性方程：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msubsup> <mn>1</mn> <mi>v</mi> <mi>T</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msub> <mn>1</mn> <mi>v</mi> </msub> </mtd> <mtd> <mrow> <mi>K</mi> <mo>+</mo> <msub> <mi>V</mi> <mi>&amp;gamma;</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msup> <mi>b</mi> <mo>*</mo> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mi>&amp;alpha;</mi> <mo>*</mo> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mi>Y</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中
权重因子vi由下式计算：
其中是误差变量ξi标准差的估计，c1,c2为常量；
求解得待估计函数f(x)：
<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msup> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <mo>*</mo> </msup> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mo>|</mo> <mo>|</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>/</mo> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>b</mi> <mo>*</mo> </msup> </mrow>
其中，M是支持向量的数目，1v＝[1,...,1]T,上标T表示矩阵的转置，是拉格朗日乘子，b*是偏置量，K＝exp(-||xi-xj||/θ2)，其中i＝1,…,M，j＝1,…,M，和exp(-||x-xi||/θ2)均为支持向量机的核函数，xj为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；
(6)用混合蛙跳算法对步骤4)的核参数θ和惩罚系数γ进行优化，采用如下过程完成：
(6.1)初始化青蛙群体参数，设种群青蛙数为P,最大迭代次数Maxgen，局部搜索的迭代次数Mmax，最大更新长度Dmax，分组数m以及每组青蛙数n，由于模型有两个参数需要优化，所以位置pi的维度为2维，随机生成每个青蛙的位置pi＝(pi1,pi2)，置初始迭代次数k＝0；
(6.2)计算所有青蛙的适应度值，并进行排序，分组，选取种群最优青蛙pg；
(6.3)按照下式对子群最差青蛙对子群中最差的青蛙进行更新，接着在子群内重新排序，然后对子群最差青蛙进行更新；重复该局部搜索过程Mmax次；
D(k+1)＝W×(r1×D(k)+r2×D(k-1)+r3×His(pw))+r4×(pb-pw)
p'w(k+1)＝pw(k)+D,-Dmax≤D≤Dmax
W＝We+(Ws-We)×a
a＝exp(-30×(k/Maxgen)s)
其中，pw为子群最差青蛙，pb为子群最优青蛙，Dmax为最大变异尺度，p'w为更新后的青蛙，r1、r2、r3、r4是0～1之间的随机数，k是当前迭代次数，W是权重因子，Ws、We是权重因子的起始值和最终值，s为大于1的整数。首先利用子群最优青蛙去更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则种群最优青蛙代替子群最优青蛙进行更新，若新得到的青蛙优于原来子群最差青蛙，则取代之；否则随机生成一个青蛙取代原来子群最差青蛙。
(6.4)当所有子群的局部搜索都已经完成，对所有青蛙进行混合、排序和分组，选取种群最优青蛙pg；
(6.5)k＝k+1，若k<Maxgen，则转至步骤3；否则输出种群最优青蛙xg，为鲁棒预报模型的最优参数，算法终止；
其中，初始种群大小为200，分组数10，每组子群数20，种群最大迭代数100，子群最大迭代数10，最大更新长度5，权重因子的起始值和最终值分别为0.9和0.4，s为3。
(7)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[xt-D+1,…,xt]，TX表示海杂波从第t-D+1采样时刻到第t采样时刻的信号幅值矩阵，xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，xt表示第t采样时刻的海杂波回波信号幅值；
(8)进行归一化处理；
<mrow> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mi>T</mi> <mi>X</mi> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> <mrow> <mi>max</mi> <mi> </mi> <mi>x</mi> <mo>-</mo> <mi>min</mi> <mi> </mi> <mi>x</mi> </mrow> </mfrac> </mrow>
(9)代入步骤(5)得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。
(10)计算海杂波预报值与雷达回波实测值的差值e，计算控制限Qα：
<mrow> <msub> <mi>Q</mi> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <msup> <mrow> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mi>&amp;alpha;</mi> </msub> <msub> <mi>h</mi> <mn>0</mn> </msub> <msqrt> <mrow> <mn>2</mn> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> </mrow> </msqrt> </mrow> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <msub> <mi>h</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> </mfrac> <mo>&amp;rsqb;</mo> </mrow> <mfrac> <mn>1</mn> <msub> <mi>h</mi> <mn>0</mn> </msub> </mfrac> </msup> </mrow>
<mrow> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>&amp;lambda;</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> </mrow>
<mrow> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> </mrow> <mrow> <mn>3</mn> <msubsup> <mi>&amp;theta;</mi> <mn>2</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow>
其中，α是置信度，θ1,θ2,θ3,h0是中间变量，λj i表示协方差矩阵的第j个特征值的i次方，k是样本维数，Cα是正态分布置信度为α的统计；
(11)进行检测判断：当e2差值大于控制限Qα时，该点存在目标，否则没有目标。
(12)按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。
CN201711115256.7A 2017-11-13 2017-11-13 一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法 Pending CN107942299A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711115256.7A CN107942299A (zh) 2017-11-13 2017-11-13 一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711115256.7A CN107942299A (zh) 2017-11-13 2017-11-13 一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法

Publications (1)

Publication Number Publication Date
CN107942299A true CN107942299A (zh) 2018-04-20

Family

ID=61933967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711115256.7A Pending CN107942299A (zh) 2017-11-13 2017-11-13 一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法

Country Status (1)

Country Link
CN (1) CN107942299A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108647491A (zh) * 2018-05-14 2018-10-12 浙江大学 一种智能白血病癌细胞检测仪

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7372393B2 (en) * 2006-07-07 2008-05-13 Mitsubishi Electric Research Laboratories, Inc. Method and system for determining unwrapped phases from noisy two-dimensional wrapped-phase images
CN102147464A (zh) * 2011-03-03 2011-08-10 浙江大学 一种鲁棒智能雷达海杂波预报系统及方法
CN102147463A (zh) * 2011-03-03 2011-08-10 浙江大学 一种群智雷达海杂波预报系统及方法
CN102147465A (zh) * 2011-03-03 2011-08-10 浙江大学 一种混沌寻优雷达海上目标检测系统及方法
CN102183749A (zh) * 2011-03-03 2011-09-14 浙江大学 一种自适应雷达海上目标检测系统及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7372393B2 (en) * 2006-07-07 2008-05-13 Mitsubishi Electric Research Laboratories, Inc. Method and system for determining unwrapped phases from noisy two-dimensional wrapped-phase images
CN102147464A (zh) * 2011-03-03 2011-08-10 浙江大学 一种鲁棒智能雷达海杂波预报系统及方法
CN102147463A (zh) * 2011-03-03 2011-08-10 浙江大学 一种群智雷达海杂波预报系统及方法
CN102147465A (zh) * 2011-03-03 2011-08-10 浙江大学 一种混沌寻优雷达海上目标检测系统及方法
CN102183749A (zh) * 2011-03-03 2011-09-14 浙江大学 一种自适应雷达海上目标检测系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘悦婷: "基于混合蛙跳的数据挖掘模糊聚类算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108647491A (zh) * 2018-05-14 2018-10-12 浙江大学 一种智能白血病癌细胞检测仪

Similar Documents

Publication Publication Date Title
CN102183749B (zh) 一种自适应雷达海上目标检测系统及方法
Ji et al. Seabed sediment classification using multibeam backscatter data based on the selecting optimal random forest model
CN104297748B (zh) 一种基于轨迹增强的雷达目标检测前跟踪方法
Pennino et al. Modeling sensitive elasmobranch habitats
CN102147465B (zh) 一种混沌寻优雷达海上目标检测系统及方法
US8138963B1 (en) Method for detecting targets using space-time adaptive processing and shared knowledge of the environment
CN108983179A (zh) 一种群智能敏捷的雷达海上目标检测系统
CN105372636A (zh) 一种基于自适应遗忘因子的杂波图更新方法及装置
CN107942312A (zh) 一种基于差分进化入侵杂草优化算法的智能雷达海上目标检测系统及方法
CN108931776A (zh) 一种高精度的匹配场定位方法
CN107656250A (zh) 一种基于人工蜂群算法的智能雷达海上目标检测系统及方法
Yan et al. Floating small target detection in sea clutter using mean spectral radius
CN108983181A (zh) 一种群智寻优的雷达海上目标检测系统
CN108983178A (zh) 一种敏捷自适应的智能雷达海上目标检测系统
CN107942299A (zh) 一种基于改进混合蛙跳算法的智能雷达海上目标检测系统及方法
CN102183746B (zh) 一种雷达海上目标检测系统及方法
CN102183751B (zh) 一种智能雷达海上目标检测系统及方法
CN102183744B (zh) 一种群智雷达海上目标检测系统及方法
CN107942300A (zh) 一种基于改进人工蜂群算法的智能雷达海上目标检测系统及方法
CN110287924A (zh) 一种基于gru-rnn模型的土壤参数分类方法
CN107976652A (zh) 一种基于混合蛙跳算法的智能雷达海上目标检测系统及方法
CN107918117A (zh) 一种基于自适应混合蛙跳算法的智能雷达海上目标检测系统及方法
CN107656248A (zh) 一种基于混合蛙跳算法的智能雷达海杂波预报系统及方法
Dai et al. Learning From Noise: An Unsupervised GPR Data Denoising Scheme based on Generative Adversarial Networks
CN102147466A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180420

RJ01 Rejection of invention patent application after publication