发明名称 基于噪声先验的联合高斯模型雷达目标稳健识别方法
摘要 本发明提出了一种基于噪声先验的联合高斯模型雷达目标稳健识别方法,主要解决现有雷达飞机目标识别技术中的统计模型对噪声不稳健的问题。其实现过程是:对雷达高分辨距离像数据做预处理;将预处理之后的数据取模获得其时域特征;确定雷达训练目标高分辨距离像数据每一帧的联合高斯模型的均值、加载矩阵和噪声协方差矩阵;统计雷达测试高分辨距离像数据非信号支撑区中噪声的方差,并以此计算联合高斯模型的均值和噪声协方差矩阵的修正值;计算雷达测试目标高分辨距离像数据时域特征对应每一个雷达训练目标的后验概率值;确定雷达测试目标高分辨距离像数据的类别属性。本发明具有对噪声稳健的优点,可用于对雷达飞机目标的稳健识别。
申请公布号 CN102628938B 申请公布日期 2013.08.14
申请号 CN201210131026.0 申请日期 2012.04.29
申请人 西安电子科技大学 发明人 刘宏伟;潘勉;杜兰;张学峰;冯博;王鹏辉
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.基于噪声先验的联合高斯模型雷达目标稳健识别方法,包括如下步骤: (1)对雷达训练目标高分辨距离像数据进行角域分帧、帧内对齐和信号能量归一化的预处理,对雷达测试目标高分辨距离像数据进行信号能量归一化的预处理; (2)对预处理后的雷达训练目标高分辨距离像数据和雷达测试目标高分辨距离像数据取模得到它们的时域特征; (3)用雷达训练目标高分辨距离像数据的每一帧时域特征,来训练联合高斯模型,先求得雷达训练目标高分辨距离像数据每一帧联合高斯模型的均值μ,并通过最大期望EM算法确定雷达训练目标高分辨距离像数据每一帧联合高斯模型的加载矩阵Σ<sub>x</sub>和噪声协方差矩阵ψ<sub>noise</sub>; (4)统计预处理后的雷达测试目标高分辨距离像中非信号支撑区的噪声方差<img file="FDA00003112672500012.GIF" wi="66" he="67" />;(5)根据雷达测试目标高分辨距离像中非信号支撑区的噪声方差<img file="FDA00003112672500013.GIF" wi="64" he="65" />,对上述联合高斯模型的均值μ和噪声协方差ψ<sub>noise</sub>进行修正,得到雷达训练目标高分辨距离像数据每一帧联合高斯模型的均值修正值μ<sub>0</sub>和噪声协方差修正值ψ<sub>0</sub>:(5a)计算对应于第c个目标第m帧的雷达训练目标的辅助函数f(x<sup>(c,m)</sup>): <img file="FDA00003112672500011.GIF" wi="1647" he="199" />其中,上标(c,m)表示参数对应于第c个目标的第m帧,w和v为用于积分的变量,x<sup>(c,m)</sup>(l)表示雷达训练目标高分辨距离像x<sup>(c,m)</sup>第l个距离单元,μ<sup>(c,m)</sup>(l)表示第c个目标第m帧联合高斯模型的均值中的第l个元素,<img file="FDA00003112672500014.GIF" wi="105" he="68" />(l)表示第c个目标第m帧联合高斯模型的加载矩阵中第l行,(·)<sup>T</sup>表示矩阵的转置操作,L表示雷达训练目标高分辨距离像时域特征中距离单元个数;(5b)根据辅助函数f(x<sup>(c,m)</sup>),计算第c个目标第m帧联合高斯模型的均值修正值<img file="FDA00003112672500015.GIF" wi="108" he="70" />:5b1)根据辅助函数f(x<sup>(c,m)</sup>),计算辅助向量α<sup>(c,m)</sup>,α<sup>(c,m)</sup>是一个L×1维向量,它的第l个元素α<sup>(c,m)</sup>(l)表示为: <img file="FDA00003112672500021.GIF" wi="1166" he="106" />5b2)计算第c个目标第m帧联合高斯模型的均值修正值<img file="FDA00003112672500027.GIF" wi="108" he="71" />:<img file="FDA00003112672500022.GIF" wi="447" he="78" />其中,μ<sup>(c,m)</sup>表示第c个目标第m帧联合高斯模型的均值; (5c)根据辅助函数f(x<sup>(c,m)</sup>)计算第c个目标第m帧联合高斯模型的噪声协方差的修正值<img file="FDA00003112672500028.GIF" wi="120" he="69" />:5c1)根据辅助函数f(x<sup>(c,m)</sup>),计算第c个目标第m帧的雷达训练目标的辅助矩阵φ<sup>(c,m)</sup>,φ<sup>(c,m)</sup>是L×L维的矩阵,它的第l<sub>1</sub>行l<sub>2</sub>列的元素φ<sup>(c,m)</sup>(l<sub>1</sub>,l<sub>2</sub>)表示为: <img file="FDA00003112672500023.GIF" wi="1584" he="100" />其中α<sup>(c,m)</sup>(l)为辅助向量α<sup>(c,m)</sup>中的第l个的元素; 5c2)根据辅助矩阵φ<sup>(c,m)</sup>,计算第c个目标第m帧联合高斯模型的噪声协方差的修正值<img file="FDA00003112672500029.GIF" wi="111" he="71" />:<img file="FDA00003112672500024.GIF" wi="852" he="97" />其中diag(·)表示令矩阵非对角线元素置零操作,<img file="FDA00003112672500025.GIF" wi="116" he="81" />表示雷达训练目标高分辨距离像数据时域特征x<sup>(c,m)</sup>对应于第c个目标第m帧联合高斯噪声协方差矩阵;(6)根据雷达训练目标高分辨距离像数据每一帧联合高斯模型的均值μ、均值的修正值μ<sub>0</sub>、加载矩阵Σ<sub>x</sub>、噪声协方差矩阵ψ<sub>noise</sub>及噪声协方差修正值ψ<sub>0</sub>,计算雷达测试目标高分辨距离像数据时域特征x<sub>test</sub>与每一个雷达训练目标相对应的后验概率值p(x<sub>test</sub>|c): <img file="FDA00003112672500026.GIF" wi="1634" he="130" />其中,<img file="FDA00003112672500031.GIF" wi="675" he="87" />,它表示雷达测试目标高分辨距离像数据时域特征x<sub>test</sub>对应于第c个目标第m帧联合高斯模型协方差矩阵,m=1,...,M<sub>c</sub>,c=1,...,N<sub>c</sub>,<img file="FDA00003112672500032.GIF" wi="127" he="80" />表示雷达训练目标高分辨距离像数据时域特征x<sup>(c,m)</sup>对应于第c个目标第m帧联合高斯噪声协方差矩阵,上标(c,m)表示参数对应于第c个目标的第m帧,det(·)表示求矩阵的行列式操作,(B<sup>(c,m)</sup>)<sup>-1</sup>表示求协方差矩阵B<sup>(c,m)</sup>的逆,M<sub>c</sub>表示第c个目标的包含的帧数,N<sub>c</sub>表示雷达训练目标总数;(7)取雷达训练目标相对应的后验概率值p(x<sub>test</sub>|c)中最大值所对应的雷达训练目标的类别属性,作为雷达测试目标高分辨距离像数据的类别属性。 
地址 710071 陕西省西安市太白南路2号