发明名称 基于稀疏特征的雷达目标识别方法
摘要 本发明提供了一种基于稀疏特征的雷达目标识别方法,该方法以雷达目标图像的稀疏特征作为雷达目标图像训练样本和待测雷达目标的识别特征,帮助提升雷达目标识别的数据针对性并减少数据处理运算量,再通过训练样本的稀疏特征构建待测雷达目标的稀疏线性方程,然后通过贝叶斯压缩感知算法求解稀疏线性方程,借助压缩感知理论实现对待测雷达目标的识别,不需要借助目标方位角估计,降低了识别复杂程度,避免了识别准确性对目标方位角估计的依赖,同时基于压缩感知理论进行雷达目标识别在噪声环境下也具有良好的识别性能,从而解决了现有技术中雷达目标识别系统较为复杂、识别准确性有限的问题,达到提高雷达目标识别的处理效率和识别准确性的目的。
申请公布号 CN103226196B 申请公布日期 2015.05.13
申请号 CN201310184862.X 申请日期 2013.05.17
申请人 重庆大学 发明人 张新征;吴奇政;秦建红;谭熠峰
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 重庆博凯知识产权代理有限公司 50212 代理人 张先芸
主权项 基于稀疏特征的雷达目标识别方法,其特征在于,包括如下步骤:A)针对多类不同的已知雷达目标,分别采集多个已知雷达目标图像作为训练样本,并分别提取各个类别中各个训练样本的稀疏特征矩阵;该步骤具体为:a1)针对多类不同的已知雷达目标,分别采集多个已知雷达目标图像作为训练样本;所采集的第i类已知雷达目标的第n<sub>i</sub>个训练样本<img file="FDA0000612842130000011.GIF" wi="72" he="68" />为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>r</mi><msub><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow><mi>i</mi></msub></msub><mo>=</mo><msub><mi>C</mi><msub><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow><mi>i</mi></msub></msub><mo>&times;</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></msub><mo>+</mo><msub><mi>w</mi><msub><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow><mi>i</mi></msub></msub><mo>;</mo></mrow>]]></math><img file="FDA0000612842130000012.GIF" wi="463" he="78" /></maths>其中,<img file="FDA0000612842130000013.GIF" wi="86" he="71" />表示第i类已知雷达目标第n<sub>i</sub>个训练样本的目标反射率图像矩阵;<img file="FDA0000612842130000014.GIF" wi="90" he="71" />表示第i类已知雷达目标第n<sub>i</sub>个训练样本的叠加噪声;<img file="FDA0000612842130000015.GIF" wi="91" he="78" />表示第i类已知雷达目标第n<sub>i</sub>个训练样本的观测转换矩阵;所述观测转换矩阵<img file="FDA0000612842130000016.GIF" wi="90" he="78" />为已知参数,所述叠加噪声<img file="FDA0000612842130000017.GIF" wi="80" he="70" />和目标反射率图像矩阵<img file="FDA0000612842130000018.GIF" wi="78" he="69" />为未知参数;1≤i≤K,K表示所采集的已知雷达目标的类别总数;1≤n<sub>i</sub>≤N<sub>i</sub>,N<sub>i</sub>表示针对第i类已知雷达目标所采集的训练样本的个数;a2)采用投影梯度下降算法对下式求解,得到第i类已知雷达目标第n<sub>i</sub>个训练样本的目标反射率图像矩阵<img file="FDA0000612842130000019.GIF" wi="74" he="66" />的最小取值<img file="FDA00006128421300000110.GIF" wi="106" he="84" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>u</mi><mo>^</mo></mover><mrow><mi>i</mi><mo>,</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></msub><mo>=</mo><mi>arg</mi><munder><mi>min</mi><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></msub></munder><mo>[</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></msub><mo>-</mo><msub><mi>C</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></msub><mo>&times;</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>+</mo><mi>&lambda;</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>]</mo><mo>;</mo></mrow>]]></math><img file="FDA00006128421300000111.GIF" wi="910" he="120" /></maths>并以该最小取值<img file="FDA00006128421300000112.GIF" wi="78" he="83" />作为第i类已知雷达目标第n<sub>i</sub>个训练样本的稀疏特征矩阵<img file="FDA00006128421300000113.GIF" wi="106" he="84" />即取<img file="FDA00006128421300000114.GIF" wi="224" he="80" />由此方式分别提取各个类别中各个训练样本的稀疏特征矩阵;其中,<img file="FDA00006128421300000115.GIF" wi="228" he="112" />表示使得目标反射率图像矩阵<img file="FDA00006128421300000116.GIF" wi="76" he="67" />的目标函数为最小值时<img file="FDA00006128421300000117.GIF" wi="76" he="70" />的取值;||·||<sub>1</sub>为L1范数运算符;||·||<sub>2</sub>为L2范数运算符,λ为L2范数规范参数,且0.1≤λ≤0.7;B)将步骤A提取的每个训练样本的稀疏特征矩阵作为一个稀疏特征训练样本,从而由步骤A提取的各个类别各个训练样本的稀疏特征矩阵的集合构成稀疏特征训练样本集;C)针对待测雷达目标,采集待测雷达目标图像,提取待测雷达目标图像的稀疏特征矩阵;该步骤具体为:c1)针对待测雷达目标,采集待测雷达目标图像r<sub>z</sub>:r<sub>z</sub>=C<sub>z</sub>×u<sub>z</sub>+w<sub>z</sub>;其中,u<sub>z</sub>表示待测雷达目标图像的目标反射率图像矩阵;w<sub>z</sub>表示待测雷达目标图像的叠加噪声;C<sub>z</sub>表示待测雷达目标图像的观测转换矩阵;所述观测转换矩阵C<sub>z</sub>为已知参数,所述叠加噪声w<sub>z</sub>和目标反射率图像矩阵u<sub>z</sub>为未知参数;c2)采用投影梯度下降算法对下式求解,得到待测雷达目标图像的目标反射率图像矩阵u<sub>z</sub>的最小取值<img file="FDA0000612842130000021.GIF" wi="79" he="73" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>u</mi><mo>^</mo></mover><mi>z</mi></msub><mo>=</mo><mi>arg</mi><munder><mi>min</mi><msub><mi>u</mi><mi>z</mi></msub></munder><mo>[</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>z</mi></msub><mo>-</mo><msub><mi>C</mi><mi>z</mi></msub><mo>&times;</mo><msub><mi>u</mi><mi>z</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>+</mo><mi>&lambda;</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>u</mi><mi>z</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>]</mo><mo>;</mo></mrow>]]></math><img file="FDA0000612842130000022.GIF" wi="775" he="108" /></maths>并以该最小取值<img file="FDA0000612842130000023.GIF" wi="56" he="72" />作为待测雷达目标图像的稀疏特征矩阵b<sub>z</sub>,即取<img file="FDA0000612842130000024.GIF" wi="172" he="72" />其中,<img file="FDA0000612842130000025.GIF" wi="226" he="101" />表示使得目标反射率图像矩阵u<sub>z</sub>的目标函数为最小值时u<sub>z</sub>的取值;||·||<sub>1</sub>为L1范数运算符;||·||<sub>2</sub>为L2范数运算符,λ为L2范数规范参数,且0.1≤λ≤0.7;D)利用稀疏特征训练样本集中个各个稀疏特征训练样本建立稀疏线性方程,对待测雷达目标图像的稀疏特征矩阵进行稀疏线性表示,并采用贝叶斯压缩感知算法求解得到该稀疏线性方程的系数向量;E)分别提取步骤D所得系数向量对应于稀疏特征训练样本集中每一类已知雷达目标的稀疏系数向量,然后分别计算对应于每一类已知雷达目标的稀疏系数向量的L2范数值,将L2范数值最大的稀疏系数向量所对应的一个已知雷达目标类别判定为待测雷达目标所属的雷达目标类别,实现对待测雷达目标的识别。
地址 400044 重庆市沙坪坝区沙正街174号