发明名称 一种用于船用导航雷达的弱小目标检测方法
摘要 本发明公开了一种用于船用导航雷达的弱小目标检测方法,具体包括训练过程和检测过程,训练过程包括海杂波数据相空间重构和灰色神经网络训练,检测过程具体为海杂波目标检测。本发明的方法基于海杂波的回波内在特性和目标的内在特性不同,利用不含目标的海杂波数据训练灰色神经网络,在纯海杂波情况下总误差均方根值或对消误差趋近零;含有目标的总误差均方根值和对消误差就会很大,因而可以进行弱小目标检测的。本发明的方法相对传统的恒虚警目标检测法,可以在强海杂波背景下检测出弱小目标;相比RBF神经网络进行检测的方法,训练速度相对较快,需要的样本数据信息较少,在海杂波背景下有良好的检测弱小目标的性能。
申请公布号 CN102621533B 申请公布日期 2014.06.18
申请号 CN201210098333.3 申请日期 2012.04.05
申请人 电子科技大学 发明人 卢宁;陈华;唐伟;韩世雄;许宏志
分类号 G01S7/40(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 周永宏
主权项 1.一种用于船用导航雷达的弱小目标检测方法,其特征在于,包括训练过程和检测过程,具体的,所述训练过程包括如下步骤:利用海杂波训练数据进行海杂波相空间的重构,具体包括如下分步骤:S11:计算嵌入维数m;S12:计算延迟时间τ;S13:根据延迟时间τ和嵌入维数m,把海杂波训练数据从时间序列构造成奇异吸引子轨迹向量X<sub>j</sub>=(x<sub>j</sub>,x<sub>j+τ</sub>,…,x<sub>j+(m-1)τ</sub>),j=1,…,n,其中,x<sub>j</sub>表示海杂波训练数据中第j个采样点的海杂波数据,n根据海杂波训练数据总的数目和嵌入维数m得到;利用重构得到的奇异吸引子轨迹向量对灰色神经网络进行训练,所述灰色神经网络具体包括如下参数,输入参数序号t;网络输入参数y<sub>2</sub>(t),…,y<sub>n</sub>(t);网络权值w<sub>11</sub>,w<sub>21</sub>,w<sub>22</sub>,…,w<sub>2n</sub>,w<sub>31</sub>,w<sub>32</sub>,…,w<sub>3n</sub>;网络预测值y<sub>1</sub>;灰色神经网络的四层:第一层LA、第二层LB、第三层LC、第四层LD,具体包括如下分步骤:S21.根据灰色神经网络特征初始化网络结构;S22.根据网络结构调节权值w<sub>11</sub>,w<sub>21</sub>,w<sub>22</sub>,…,w<sub>2n</sub>,w<sub>31</sub>,w<sub>32</sub>,…,w<sub>3n</sub>;S23.对每一个输入序列(t,S<sub>t</sub>),t=1,2,3,…,n,计算每层输出:LA层:a=w<sub>11</sub>t;LB层:<maths num="0001"><![CDATA[<math><mrow><mi>b</mi><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi></mrow></msup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>LC层:c<sub>1</sub>=bw<sub>21</sub>,c<sub>2</sub>=y<sub>2</sub>(t)bw<sub>22</sub>,c<sub>3</sub>=y<sub>3</sub>(t)bw<sub>23</sub>,…,c<sub>n</sub>=y<sub>n</sub>(t)bw<sub>2n</sub>,y<sub>2</sub>(t),…,y<sub>n</sub>(t)与S<sub>t</sub>,t=2,3,…,n一一对应;LD层:d=w<sub>31</sub>c<sub>1</sub>+w<sub>32</sub>c<sub>2</sub>+…+w<sub>3n</sub>c<sub>n</sub>-θ,θ表示第四层LD网络的预测阈值;S24.计算网络预测输出与期望输出的误差,并根据误差调整权值和阈值:LD层误差:δ=d-y<sub>1</sub>;LC层误差:<maths num="0002"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mn>1</mn></msub><mo>=</mo><mi>&delta;</mi><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi></mrow></msup><mo>)</mo></mrow><mo>,</mo><msub><mi>&delta;</mi><mn>2</mn></msub><mo>=</mo><mi>&delta;</mi><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi></mrow></msup><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>&delta;</mi><mi>n</mi></msub><mo>=</mo><mi>&delta;</mi><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi></mrow></msup><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>LB层误差:<maths num="0003"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi></mrow></msup></mrow></mfrac><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi></mrow></msup><mo>)</mo></mrow><mo>*</mo><mrow><mo>(</mo><msub><mi>w</mi><mn>21</mn></msub><msub><mi>&delta;</mi><mn>1</mn></msub><mo>+</mo><msub><mi>w</mi><mn>22</mn></msub><msub><mi>&delta;</mi><mn>2</mn></msub><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>w</mi><mrow><mn>2</mn><mi>n</mi></mrow></msub><msub><mi>&delta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>调整LB到LC的权值,w<sub>21</sub>=-y<sub>1</sub>(0),w<sub>22</sub>=w<sub>22</sub>-μ<sub>1</sub>δ<sub>2</sub>b,…,w<sub>2n</sub>=w<sub>2n</sub>-μ<sub>n-1</sub>δ<sub>n</sub>b,μ<sub>i</sub>表示第二层与第三层的第i+1个节点之间权值调节系数;调整LA层到LB层的权值:w<sub>11</sub>=w<sub>11</sub>+atδ<sub>n+1</sub>;调整阈值:<maths num="0004"><![CDATA[<math><mrow><mi>&theta;</mi><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>w</mi><mn>11</mn></msub><mi>t</mi></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><msub><mi>w</mi><mn>22</mn></msub><mn>2</mn></mfrac><msub><mi>y</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>w</mi><mn>23</mn></msub><mn>2</mn></mfrac><msub><mi>y</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><mfrac><msub><mi>w</mi><mrow><mn>2</mn><mi>n</mi></mrow></msub><mn>2</mn></mfrac><msub><mi>y</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>S25.判断LD层误差δ是否大于预设的期望误差值,若大于,则返回S23,继续训练,否则训练结束;所述检测过程具体包括如下步骤:将待检测海杂波数据输入到已经训练好的灰色神经网络中,得到第四层的输出,即是预测值<img file="FDA0000423335520000022.GIF" wi="81" he="62" />计算对消误差或总误差均方根;对消误差计算公式如下:<img file="FDA0000423335520000023.GIF" wi="341" he="69" />其中,x<sub>n</sub>表示待检测海杂波数据的实际值;总误差均方根计算公式如下:<img file="FDA0000423335520000024.GIF" wi="431" he="141" />若对消误差或总误差均方根趋近零,则表示待检测海杂波数据不含弱小目标;否则,则表示待检测海杂波数据含有弱小目标。
地址 611731 四川省成都市高新区(西区)西源大道2006号