发明名称 一种雷达海上目标检测系统及方法
摘要 一种雷达海上目标检测系统,包括雷达、数据库以及上位机;雷达、数据库和上位机依次相连,所述雷达对所检测海域进行照射,并将雷达海杂波数据储存到所述的数据库,所述的上位机包括数据预处理模块、最小二乘预测模型建模模块、目标检测模块、模型更新模块以及结果显示模块。以及提出了一种雷达海上目标检测方法。本发明提供一种响应速度快、在线检测效率高的雷达海上目标检测系统及方法。
申请公布号 CN102183746A 申请公布日期 2011.09.14
申请号 CN201110051096.0 申请日期 2011.03.03
申请人 浙江大学 发明人 刘兴高;闫正兵
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;王利强
主权项 1.一种雷达海上目标检测系统,包括雷达、数据库以及上位机,雷达、数据库和上位机依次相连,其特征在于:所述雷达对所检测海域进行照射,并将雷达海杂波数据储存到所述的数据库,所述的上位机包括:所述的数据预处理模块,用以进行雷达海杂波数据预处理,采用如下过程完成:1)从数据库中采集N个雷达海杂波回波信号幅值x<sub>i</sub>作为训练样本,i=1,...,N;2)对训练样本进行归一化处理,得到归一化幅值<img file="FDA0000048706210000011.GIF" wi="56" he="47" /><maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>&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>x</mi></mrow><mrow><mi>max</mi><mi>x</mi><mo>-</mo><mi>min</mi><mi>x</mi></mrow></mfrac></mrow></math>]]></maths>其中,minx表示训练样本中的最小值,maxx表示训练样本中的最大值;3)将归一化后的训练样本重构,分别得到输入矩阵X和对应的输出矩阵Y:<maths num="0002"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub></mtd><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub></mtd><mtd><mi>L</mi></mtd><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>D</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub></mtd><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>3</mn></msub></mtd><mtd><mi>L</mi></mtd><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>D</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mi>O</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>N</mi><mo>-</mo><mi>D</mi></mrow></msub></mtd><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>N</mi><mo>-</mo><mi>D</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>L</mi></mtd><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mi>Y</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>D</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>D</mi><mo>+</mo><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><mi>M</mi></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>N</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,D表示重构维数,D为自然数,且D<N,D的取值范围为50-70;最小二乘预报模型建模模块,用以建立预报模型,采用如下过程完成:将得到的X、Y代入如下线性方程:<maths num="0004"><![CDATA[<math><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><mi>K</mi><mrow><mo>+</mo><msup><mi>&gamma;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow><mi>I</mi></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>b</mi></mtd></mtr><mtr><mtd><mi>&alpha;</mi></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></math>]]></maths>求解得到待估计函数f(x):<maths num="0005"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><msub><mi>&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>&theta;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><msup><mi>b</mi><mo>*</mo></msup></mrow></math>]]></maths>其中,M是支持向量的数目,1<sub>v</sub>=[1,...,1]<sup>T</sup>,<img file="FDA0000048706210000017.GIF" wi="380" he="57" />K=exp(-‖x<sub>i</sub>-x<sub>j</sub>‖/θ<sup>2</sup>),I为单位阵,上标-1表示矩阵的逆,上标T表示矩阵的转置,<img file="FDA0000048706210000018.GIF" wi="48" he="57" />是拉格朗日乘子,b<sup>*</sup>是偏置量,i=1,…,M,j=1,…,M,<img file="FDA0000048706210000019.GIF" wi="407" he="62" />和exp(-‖x-x<sub>i</sub>‖/θ<sup>2</sup>)均为支持向量机的核函数,x<sub>j</sub>为第j个雷达海杂波回波信号幅值,θ是核参数,x表示输入变量,γ是惩罚系数;目标检测模块,用以进行目标检测,采用如下过程完成:1)在采样时刻t采集D个海杂波回波信号幅值得到TX=[x<sub>t-D+1</sub>,K,x<sub>t</sub>],x<sub>t-D+1</sub>表示第t-D+1采样时刻的海杂波回波信号幅值,x<sub>t</sub>表示第t采样时刻的海杂波回波信号幅值;2)进行归一化处理;<maths num="0006"><![CDATA[<math><mrow><mover><mi>TX</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><mi>TX</mi><mo>-</mo><mi>min</mi><mi>x</mi></mrow><mrow><mi>max</mi><mi>x</mi><mo>-</mo><mi>min</mi><mi>x</mi></mrow></mfrac></mrow></math>]]></maths>3)代入最小二乘预报模型建模模块得到的待估计函数f(x)得到采样时刻(t+1)的海杂波预报值;4)计算海杂波预报值与雷达回波实测值的差值e,计算控制限Q<sub>α</sub>:<maths num="0007"><![CDATA[<math><mrow><msub><mi>Q</mi><mi>&alpha;</mi></msub><mo>=</mo><msub><mi>&theta;</mi><mn>1</mn></msub><msup><mrow><mo>[</mo><mfrac><mrow><msub><mi>C</mi><mi>&alpha;</mi></msub><msub><mi>h</mi><mn>0</mn></msub><msqrt><mn>2</mn><msub><mi>&theta;</mi><mn>2</mn></msub></msqrt></mrow><msub><mi>&theta;</mi><mn>1</mn></msub></mfrac><mo>+</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>&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>&theta;</mi><mn>2</mn></msub></mfrac><mo>]</mo></mrow><mfrac><mn>1</mn><msub><mi>h</mi><mn>0</mn></msub></mfrac></msup></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>&lambda;</mi><mi>j</mi><mi>i</mi></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3</mn></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msub><mi>h</mi><mn>0</mn></msub><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mn>2</mn><msub><mi>&theta;</mi><mn>1</mn></msub><msub><mi>&theta;</mi><mn>3</mn></msub></mrow><mrow><mn>3</mn><msubsup><mi>&theta;</mi><mn>2</mn><mn>2</mn></msubsup></mrow></mfrac></mrow></math>]]></maths>其中,α是置信度,θ<sub>1</sub>,θ<sub>2</sub>,θ<sub>3</sub>,h<sub>0</sub>是中间变量,<img file="FDA0000048706210000025.GIF" wi="50" he="61" />表示协方差矩阵的第j个特征值的i次方,k是样本维数,C<sub>α</sub>是正态分布置信度为α的统计;5)进行检测判断:当e<sup>2</sup>差值大于控制限Q<sub>α</sub>时,该点存在目标,否则没有目标。
地址 310027 浙江省杭州市西湖区浙大路38号