CN102183750A - 一种鲁棒雷达海杂波预报系统及方法 - Google Patents
一种鲁棒雷达海杂波预报系统及方法 Download PDFInfo
- Publication number
- CN102183750A CN102183750A CN 201110051117 CN201110051117A CN102183750A CN 102183750 A CN102183750 A CN 102183750A CN 201110051117 CN201110051117 CN 201110051117 CN 201110051117 A CN201110051117 A CN 201110051117A CN 102183750 A CN102183750 A CN 102183750A
- Authority
- CN
- China
- Prior art keywords
- overbar
- radar
- clutter
- training sample
- extra large
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000010606 normalization Methods 0.000 claims description 12
- 238000006467 substitution reaction Methods 0.000 claims description 12
- 238000012706 support-vector machine Methods 0.000 claims description 6
- 230000017105 transposition Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 3
- 230000003760 hair shine Effects 0.000 claims description 2
- 230000001678 irradiating effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
一种鲁棒雷达海杂波预报系统,包括雷达、数据库以及上位机;雷达、数据库和上位机依次相连:所述雷达对所检测海域进行照射,并将雷达海杂波数据储存到所述的数据库,所述的上位机包括数据预处理模块、鲁棒预报模型建模模块、海杂波预报模块、判别模型更新模块以及结果显示模块。以及提出了一种鲁棒雷达海杂波预报方法。本发明提供一种实现在线预测雷达海杂波、具有鲁棒性的雷达海杂波预报系统及方法。
Description
技术领域
本发明涉及雷达数据处理领域,特别地,涉及一种鲁棒雷达海杂波预报系统及方法。
背景技术
海杂波,即来自于被雷达发射信号照射的一片海面的后向散射回波。由于海杂波对来自海面或接近海面的“点”目标,如航海浮标及漂浮在海上的冰块等目标的雷达回波的可检测性形成严重制约,因此海杂波的研究对海洋背景中轮船等目标的探测性能具有十分重要的影响从而具有重要理论意义和实用价值。
习惯上海杂波被视为单一随机过程,如对数正态分布、K分布等。然而这些模型在实际应用中都有其特定的局限性,其中一个重要原因是海杂波看似随机的波形,实际上并不具有随机分布特性。
发明内容
为了克服传统雷达海杂波预报方法鲁棒性能较差的不足,本发明提供一种鲁棒雷达海杂波预报系统及方法。
本发明解决其技术问题所采用的技术方案是:
一种鲁棒雷达海杂波预报系统,包括雷达、数据库以及上位机,雷达、数据库和上位机依次相连,所述雷达对所检测海域进行照射,并将雷达海杂波数据储存到所述的数据库,所述的上位机包括:
数据预处理模块,用以进行雷达海杂波数据预处理,采用如下过程完成:
1)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本,i=1,...,N;
其中,minx表示训练样本中的最小值,maxx表示训练样本中的最大值;
3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
鲁棒预报模型建模模块,用以建立预报模型,采用如下过程完成:
将得到的X、Y代入如下线性方程:
其中
权重因子vi由下式计算:
求解得待估计函数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)在采样时刻t采集D个海杂波回波信号幅值得到TX=[xt-D+1,K,xt],xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值,xt表示第t采样时刻的海杂波回波信号幅值;
2)进行归一化处理;
3)代入鲁棒预报模型建模模块得到的函数f(x)得到采样时刻(t+1)的海杂波预报值。
作为优选的一种方案:所述上位机还包括:判别模型更新模块,用以按设定的采样时间间隔采集数据,将得到的实测数据与模型预报值比较,如果相对误差大于10%,则将新数据加入训练样本数据,更新预报模型。
作为优选的另一种方案:所述上位机还包括:结果显示模块,用以将海杂波预报模块计算得到的预报值在上位机显示。
一种鲁棒雷达海杂波预报系统所使用的海杂波预报方法,所述的方法包括以下步骤:
(1)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本,i=1,...,N;
其中,minx表示训练样本中的最小值,maxx表示训练样本中的最大值;
(3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
(4)将得到的X、Y代入如下线性方程:
其中
权重因子vi由下式计算:
求解得待估计函数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表示输入变量,γ是惩罚系数;
(5)在采样时刻t采集D个海杂波回波信号幅值得到TX=[xt-D+1,K,xt],xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值,xt表示第t采样时刻的海杂波回波信号幅值;
(6)进行归一化处理;
(7)代入步骤(4)得到的待估计函数f(x)得到采样时刻(t+1)的海杂波预报值。
作为优选的一种方案:所述的方法还包括:
(8)、按设定的采样时间间隔采集数据,将得到的实测数据与模型预报值比较,如果相对误差大于10%,则将新数据加入训练样本数据,更新预报模型。
作为优选的另一种方案:在所述的步骤(7)中,将计算得到的海杂波预报值在上位机显示。
本发明的技术构思为:本发明针对雷达海杂波的混沌特性,对雷达海杂波数据进行重构,并对重构后的数据进行非线性拟合,引入鲁棒方法,从而建立雷达海杂波的鲁棒预报模型。
本发明的有益效果主要表现在:1、建立了雷达海杂波预报模型,可以在线预测雷达海杂波;2、所用的建模方法只需较少样本即可;3、鲁棒性强。
附图说明
图1是本发明所提出的系统的硬件结构图;
图2是本发明所提出的上位机的功能模块图;
具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例1
参照图1、图2,一种鲁棒雷达海杂波预报系统,包括雷达1连接的数据库2、及上位机3,雷达1、数据库2和上位机3依次相连,所述雷达1对所检测海域进行照射,并将雷达海杂波数据储存到所述的数据库2,所述的上位机3包括:
数据预处理模块4,用以进行雷达海杂波数据预处理,采用如下过程完成:
1)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本,i=1,...,N;
其中,minx表示训练样本中的最小值,maxx表示训练样本中的最大值;
3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
鲁棒预报模型建模模块5,用以建立预报模型,采用如下过程完成:
将得到的X、Y代入如下线性方程:
其中
权重因子vi由下式计算:
求解得待估计函数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)在采样时刻t采集D个海杂波回波信号幅值得到TX=[xt-D+1,K,xt],xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值,xt表示第t采样时刻的海杂波回波信号幅值;
2)进行归一化处理;
3)代入鲁棒预报模型建模模块得到的待估计函数f(x)得到采样时刻(t+1)的海杂波预报值。
所述的上位机3还包括:判别模型更新模块8,按设定的采样时间间隔,采集数据,将得到的实测数据与模型预报值比较,如果相对误差大于10%,则将新数据加入训练样本数据,更新预报模型。
所述上位机3还包括:结果显示模块7,用于将海杂波预报模块计算得到的预报值在上位机显示。
所述上位机3的硬件部分包括:I/O元件,用于数据的采集和信息的传递;数据存储器,存储运行所需的数据样本和运行参数等;程序存储器,存储实现功能模块的软件程序;运算器,执行程序,实现指定的功能;显示模块,显示设置的参数和运行结果。
实施例2
参照图1、图2,一种鲁棒雷达海杂波预报方法,所述的方法包括以下步骤:
(1)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本,i=1,...,N;
其中,minx表示训练样本中的最小值,maxx表示训练样本中的最大值;
(3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
(4)将得到的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表示输入变量,γ是惩罚系数;
(5)在采样时刻t采集D个海杂波回波信号幅值得到TX=[xt-D+1,K,xt],xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值,xt表示第t采样时刻的海杂波回波信号幅值;
(6)进行归一化处理;
(7)代入步骤(4)得到的待估计函数f(x)得到采样时刻(t+1)的海杂波预报值。
所述的方法还包括:(8)、按设定的采样时间间隔,采集数据,将得到的实测数据与模型预报值比较,如果相对误差大于10%,则将新数据加入训练样本数据,更新预报模型。
所述的方法还包括:在所述的步骤(7)中将计算得到的海杂波预报值在上位机显示。
Claims (6)
1.一种鲁棒雷达海杂波预报系统,包括雷达、数据库以及上位机,雷达、数据库和上位机依次相连,其特征在于:所述雷达对所检测海域进行照射,并将雷达海杂波数据储存到所述的数据库,所述的上位机包括:
数据预处理模块,用以进行雷达海杂波数据预处理,采用如下过程完成:
1)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本,i=1,...,N;
2)对训练样本进行归一化处理,得到归一化幅值
其中,minx表示训练样本中的最小值,maxx表示训练样本中的最大值;
3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
鲁棒预报模型建模模块,用以建立预报模型,采用如下过程完成:
将得到的X、Y代入如下线性方程:
其中
权重因子vi由下式计算:
求解得待估计函数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)在采样时刻t采集D个海杂波回波信号幅值得到TX=[xt-D+1,K,xt],xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值,xt表示第t采样时刻的海杂波回波信号幅值;
2)进行归一化处理;
3)代入鲁棒预报模型建模模块得到的函数f(x)得到采样时刻(t+1)的海杂波预报值。
2.如权利要求1所述的鲁棒雷达海杂波预报系统,其特征在于:所述上位机还包括:判别模型更新模块,用以按设定的采样时间间隔采集数据,将得到的实测数据与模型预报值比较,如果相对误差大于10%,则将新数据加入训练样本数据,更新预报模型。
3.如权利要求1或2所述的鲁棒雷达海杂波预报系统,其特征在于:所述上位机还包括:结果显示模块,用以将海杂波预报模块计算得到的预报值在上位机显示。
4.一种如权利要求1所述的鲁棒雷达海杂波预报系统所使用的海杂波预报方法,其特征在于:所述的方法包括以下步骤:
(1)从数据库中采集N个雷达海杂波回波信号幅值xi作为训练样本,i=1,...,N;
(2)对训练样本进行归一化处理,得到归一化幅值
其中,minx表示训练样本中的最小值,maxx表示训练样本中的最大值;
(3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:
其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;
(4)将得到的X、Y代入如下线性方程:
其中
权重因子vi由下式计算:
求解得待估计函数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表示输入变量,γ是惩罚系数;
(5)在采样时刻t采集D个海杂波回波信号幅值得到TX=[xt-D+1,K,xt],xt-D+1表示第t-D+1采样时刻的海杂波回波信号幅值,xt表示第t采样时刻的海杂波回波信号幅值;
(6)进行归一化处理;
(7)代入步骤(4)得到的待估计函数f(x)得到采样时刻(t+1)的海杂波预报值。
5.如权利要求4所述的鲁棒雷达海杂波预报方法,其特征在于:所述的方法还包括:
(8)、按设定的采样时间间隔采集数据,将得到的实测数据与模型预报值比较,如果相对误差大于10%,则将新数据加入训练样本数据,更新预报模型。
6.如权利要求4或5所述的鲁棒雷达海杂波预报方法,其特征在于:在所述的步骤(7)中,将计算得到的海杂波预报值在上位机显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100511179A CN102183750B (zh) | 2011-03-03 | 2011-03-03 | 一种鲁棒雷达海杂波预报系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100511179A CN102183750B (zh) | 2011-03-03 | 2011-03-03 | 一种鲁棒雷达海杂波预报系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102183750A true CN102183750A (zh) | 2011-09-14 |
CN102183750B CN102183750B (zh) | 2012-07-25 |
Family
ID=44569950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100511179A Expired - Fee Related CN102183750B (zh) | 2011-03-03 | 2011-03-03 | 一种鲁棒雷达海杂波预报系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102183750B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101140324A (zh) * | 2007-10-11 | 2008-03-12 | 上海交通大学 | 海域合成孔径雷达图像点目标提取的方法 |
WO2008112361A2 (en) * | 2007-02-08 | 2008-09-18 | Raytheon Company | Methods and apparatus for log-ftc radar receivers having enhanced sea clutter model |
CN101806887A (zh) * | 2010-03-19 | 2010-08-18 | 清华大学 | 基于空间跟踪滤波器的海杂波抑制和目标检测方法 |
CN101881826A (zh) * | 2009-05-06 | 2010-11-10 | 中国人民解放军海军航空工程学院 | 扫描模式海杂波局部多重分形目标检测器 |
CN101887119A (zh) * | 2010-06-18 | 2010-11-17 | 西安电子科技大学 | 基于子带anmf海杂波中动目标检测方法 |
-
2011
- 2011-03-03 CN CN2011100511179A patent/CN102183750B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008112361A2 (en) * | 2007-02-08 | 2008-09-18 | Raytheon Company | Methods and apparatus for log-ftc radar receivers having enhanced sea clutter model |
CN101140324A (zh) * | 2007-10-11 | 2008-03-12 | 上海交通大学 | 海域合成孔径雷达图像点目标提取的方法 |
CN101881826A (zh) * | 2009-05-06 | 2010-11-10 | 中国人民解放军海军航空工程学院 | 扫描模式海杂波局部多重分形目标检测器 |
CN101806887A (zh) * | 2010-03-19 | 2010-08-18 | 清华大学 | 基于空间跟踪滤波器的海杂波抑制和目标检测方法 |
CN101887119A (zh) * | 2010-06-18 | 2010-11-17 | 西安电子科技大学 | 基于子带anmf海杂波中动目标检测方法 |
Non-Patent Citations (1)
Title |
---|
《舰船电子对抗》 20100430 郭锦成 对海雷达目标检测性能测试方法 第70-71,75页 1-6 第33卷, 第2期 2 * |
Also Published As
Publication number | Publication date |
---|---|
CN102183750B (zh) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102147463B (zh) | 一种群智雷达海杂波预报系统及方法 | |
CN102147464B (zh) | 一种鲁棒智能雷达海杂波预报系统及方法 | |
CN102183749B (zh) | 一种自适应雷达海上目标检测系统及方法 | |
CN102147465B (zh) | 一种混沌寻优雷达海上目标检测系统及方法 | |
Read et al. | Process‐guided deep learning predictions of lake water temperature | |
Watson‐Parris et al. | ClimateBench v1. 0: A benchmark for data‐driven climate projections | |
CN102183745B (zh) | 一种智能雷达海杂波预报系统及方法 | |
CN107942312A (zh) | 一种基于差分进化入侵杂草优化算法的智能雷达海上目标检测系统及方法 | |
CN102147466B (zh) | 一种敏捷雷达数据处理系统及方法 | |
CN102183746B (zh) | 一种雷达海上目标检测系统及方法 | |
Li et al. | A hybrid forecasting model of carbon emissions with optimized VMD and error correction | |
CN102183751B (zh) | 一种智能雷达海上目标检测系统及方法 | |
CN102183752B (zh) | 一种自适应雷达海杂波预报系统及方法 | |
Gaidai et al. | Multivariate risks assessment for complex bio-systems by Gaidai reliability method | |
CN102183754B (zh) | 一种鲁棒智能雷达海上目标检测系统及方法 | |
CN107656251A (zh) | 一种基于改进入侵杂草优化算法的智能雷达海杂波预报系统及方法 | |
CN102183753B (zh) | 一种混沌寻优雷达海杂波预报系统及方法 | |
CN102183750B (zh) | 一种鲁棒雷达海杂波预报系统及方法 | |
CN102183748B (zh) | 一种雷达海杂波预报系统及方法 | |
CN102183744B (zh) | 一种群智雷达海上目标检测系统及方法 | |
CN102183747B (zh) | 一种敏捷雷达目标检测系统及方法 | |
CN107942303A (zh) | 一种基于改进人工蜂群算法的智能雷达海杂波预报系统及方法 | |
CN202119905U (zh) | 敏捷雷达数据处理装置 | |
CN102156278B (zh) | 一种鲁棒雷达海上目标检测系统及方法 | |
Zuur et al. | Violation of independence–Part II |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120725 Termination date: 20130303 |