当前位置: 首页 > 专利查询>吉林大学专利>正文

一种基于黎曼流形的单基地MIMO雷达目标检测方法技术

技术编号:22882750 阅读:21 留言:0更新日期:2019-12-21 06:52
本发明专利技术涉及一种基于黎曼流形的单基地MIMO雷达目标检测方法,属于多输入多输出雷达目标检测领域。本发明专利技术通过研究基于黎曼流形的MIMO雷达目标检测方法,在单快拍下利用正则burg递推法得到接收信号和噪声的统计协方差矩阵,将接收信号统计协方差矩阵与噪声统计协方差矩阵黎曼均值之间的黎曼距离作为检测统计量,利用判决规则判断信号是否存在,该方法属于黎曼几何方法,提高了在低信噪比、单快拍时的MIMO雷达目标检测性能。

A target detection method of single base MIMO Radar Based on Riemannian manifold

【技术实现步骤摘要】
一种基于黎曼流形的单基地MIMO雷达目标检测方法
本专利技术属于多输入多输出(MIMO)雷达目标检测领域,涉及到基于黎曼流形的单基地MIMO雷达目标检测方法,适用于低信噪比、单快拍情况下单基地MIMO雷达的目标检测。
技术介绍
多输入多输出(MIMO)雷达是一种采用多个天线发射和接收的新型体制雷达,即多个发射天线同时发射相互正交的信号,再用多个接收天线接收信号,并利用波形分集等技术对接收到的信号进行处理分析,从而提高雷达目标检测的性能和估计精度。黎曼几何是19世纪中期由德国数学家黎曼提出的一种几何学理论,其中黎曼流形是一种微分流形,近年来,其在水声通信、物理学、神经网络、通信编码、图像处理等领域的应用日益广泛,已成为学术界的一个新的研究热点。在MIMO雷达目标检测的过程中,当快拍数较小时,样本协方差矩阵不能代替统计协方差矩阵;而且,在传统欧氏空间中计算点与点之间的距离,是直接利用欧氏距离进行矢量的模运算,会导致在计算流形上两点间的距离时出现偏差,从而影响MIMO雷达目标检测性能。
技术实现思路
本专利技术提供一种基于黎曼流形的单基地MIMO雷达目标检测方法,适用于低信噪比、单快拍情况下的单基地MIMO雷达目标检测;在实际噪声的背景下,提高单基地MIMO雷达的目标检测性能。本专利技术采取的技术方案是,包括下列步骤:步骤一,单基地MIMO雷达接收端接收的信号为:其中,ar(θ)=[1,e-j2πsin(θ)d/λ,…,e-j2πsin(θ)(N-1)d/λ]T为接收导向矢量,at(θ)=[1,e-j2πsin(θ)d/λ,…,e-j2πsin(θ)(M-1)d/λ]T为发射导向矢量,θ为波达方向角,λ为信号波长,d为阵元间距离,(·)T表示矩阵的转置;β为单快拍下的目标散射系数;S为M个发射波形组成的归一化正交矩阵,且SSH=IM,(·)H表示矩阵的共轭转置;Z是高斯白噪声矩阵;对接收到的信号X进行匹配滤波和矢量化,处理后得到单快拍下的接收信号矢量为:y=aβ+z其中,为MN×1维发射-接收联合导向矩阵,为Kronecker积;z为处理后的MN×1维的噪声向量;步骤二:在单快拍下,基于黎曼流形的单基地MIMO雷达目标检测方法的二元假设模型为:H0:y=zH1:y=aβ+z其中,H0表示仅存在噪声的情况,H1表示信号和噪声同时存在的情况;步骤三:构造单快拍下接收信号的统计协方差矩阵,采用正则burg递推法求解接收信号的统计协方差矩阵,过程如下:令单快拍下的接收信号矢量为:y=[y1,…,yMN]T接收信号矢量的统计协方差矩阵R可表示为一个MN×MN维的正定的托普利兹、厄密特矩阵,即:其中,(·)H表示矩阵的共轭转置,矩阵中的每个元素采用如下正则burg递推方法进行估计:r1=P1,r2=-a2P1,其中,Rk-1表示从R截取的(k-1)×(k-1)维子阵,根据AR模型的burg递推法,可得到预测误差功率P1、P2、Pk-1以及burg模型系数a1、a2、ak-1;步骤四:推导噪声统计协方差矩阵时,首先获得每个参考单元的噪声矢量,表示为:z=[z1,…,zMN]T与计算接收信号统计协方差矩阵过程相同,利用正则burg递推法得到噪声矢量的统计协方差矩阵Rz为:Rz=E[zzH]步骤五:假设黎曼流形上存在任意两个满足正定和厄米特对称的矩阵Rm和Rn,那么这两个矩阵Rm和Rn之间的黎曼距离可表示如下:其中,tr(·)表示矩阵的迹;步骤六:用表示从第l个参考单元噪声数据通过正则burg递推计算得到的噪声统计协方差矩阵Rz,收集总共L个参考单元的噪声数据,计算噪声统计协方差矩阵的黎曼均值MR,过程如下:首先求出L个噪声统计协方差矩阵的算术平均矩阵:根据黎曼流形的等距映射,定义函数:M1=M1/2对进行奇异值分解,并将左、右奇异值矩阵U和V与噪声协方差矩阵相乘,得到:此时,黎曼均值为:M'=M1'M1'H设置局部阈值为ε,精度h为:h=d(M',M)其中,d(M',M)表示M'与M之间的黎曼距离,若h>ε,用M'代替M继续迭代;若h<ε,迭代中止,此时得出黎曼均值为MR=M';步骤七:利用基于黎曼流形的目标检测方法,将接收信号的统计协方差矩阵R与噪声黎曼均值矩阵MR之间进行黎曼距离运算,并将其作为检测统计量为T,即:T=d(R,MR)步骤八:根据虚警概率以及噪声z的统计分布计算阈值,定义虚警概率Pfa与阈值rα的关系式为:根据虚警概率Pfa值求出阈值rα,其中α为置信度;步骤九:基于黎曼流形的单基地MIMO雷达目标检测的判决规则为:本专利技术所述步骤三中接收信号矢量的统计协方差矩阵R的具体推到步骤如下:单快拍下的接收信号矢量为:y=[y1,…,yMN]T对单快拍下接收信号y的统计协方差矩阵的每个元素进行正则burg计算,根据AR模型的burg递推法,首先设置前、后向预测误差以及预测误差功率的初始值:f0=b0=y定义前向误差fm和后向误差bm的递推关系表示如下:fm=fm-1+kmbm-1其中,m=1,2,…,MN为AR模型阶次,(·)*表示复数共轭,km为反射系数,表达式为:其中,分别为fm、bm从第p行开始截取的子向量,根据反射系数及Levinson-Durbin公式,得出正则burg递推法模型系数的表达式为:am=km计算此时的预测误差功率为:Pm=(1-|km|2)Pm-1利用如下正则burg递推方法估计y的统计协方差矩阵中每个元素:r1=P1,r2=-a2P1,其中,Rk-1表示从R截取的(k-1)×(k-1)维子阵,然后得到单快拍下对应的接收信号矢量的统计协方差矩阵:其中,(·)H表示矩阵的共轭转置。本专利技术首先收集测试单元的观测信号和L个参考单元的观测噪声,利用正则burg递推法分别推导出它们的统计协方差矩阵,计算L个噪声协方差矩阵的黎曼均值,建立基于黎曼流形的MIMO雷达目标检测二元假设检验模型,根据虚警概率以及噪声的统计分布计算阈值,并通过判决规则检测是否存在目标,在条件相同的情况下,比较虚警概率对检测性能的影响。本专利技术采用基于黎曼流形的检测方法,即:将经过正则burg递推法得到的接收信号统计协方差矩阵与噪声统计协方差矩阵的黎曼均值之间的黎曼距离作为检测统计量,利用判决规则判断信号是否存在,该方法属于黎曼几何方法,可以在低信噪比、单快拍的情况下进行目标检测。本专利技术的优点是:通过研究基于黎曼流形的MIMO雷达目标检测方法,在单快拍下利用正则burg递本文档来自技高网...

【技术保护点】
1.一种基于黎曼流形的单基地MIMO雷达目标检测方法,其特征在于,包括下列步骤:/n步骤一,单基地MIMO雷达接收端接收的信号为:/n

【技术特征摘要】
1.一种基于黎曼流形的单基地MIMO雷达目标检测方法,其特征在于,包括下列步骤:
步骤一,单基地MIMO雷达接收端接收的信号为:



其中,ar(θ)=[1,e-j2πsin(θ)d/λ,…,e-j2πsin(θ)(N-1)d/λ]T为接收导向矢量,at(θ)=[1,e-j2πsin(θ)d/λ,…,e-j2πsin(θ)(M-1)d/λ]T为发射导向矢量,θ为波达方向角,λ为信号波长,d为阵元间距离,(·)T表示矩阵的转置;β为单快拍下的目标散射系数;S为M个发射波形组成的归一化正交矩阵,且SSH=IM,(·)H表示矩阵的共轭转置;Z是高斯白噪声矩阵;
对接收到的信号X进行匹配滤波和矢量化,处理后得到单快拍下的接收信号矢量为:
y=aβ+z
其中,为MN×1维发射-接收联合导向矩阵,为Kronecker积;z为处理后的MN×1维的噪声向量;
步骤二:在单快拍下,基于黎曼流形的单基地MIMO雷达目标检测方法的二元假设模型为:
H0:y=z
H1:y=aβ+z
其中,H0表示仅存在噪声的情况,H1表示信号和噪声同时存在的情况;
步骤三:构造单快拍下接收信号的统计协方差矩阵,采用正则burg递推法求解接收信号的统计协方差矩阵,方法过程如下:
令单快拍下的接收信号矢量为:
y=[y1,…,yMN]T
接收信号矢量的统计协方差矩阵R可表示为一个MN×MN维的正定的托普利兹、厄密特矩阵,即:



其中,(·)H表示矩阵的共轭转置,矩阵中的每个元素采用如下正则burg递推方法进行估计:
r1=P1,
r2=-a2P1,



其中,Rk-1表示从R截取的(k-1)×(k-1)维子阵,根据AR模型的burg递推法,可得到预测误差功率P1、P2、Pk-1以及burg模型系数a1、a2、ak-1;
步骤四:推导噪声统计协方差矩阵时,首先获得每个参考单元的噪声矢量,表示为:
z=[z1,…,zMN]T
与计算接收信号统计协方差矩阵过程相同,利用正则burg递推法得到噪声矢量的统计协方差矩阵Rz为:
Rz=E[zzH]
步骤五:假设黎曼流形上存在任意两个满足正定和厄米特对称的矩阵Rm和Rn,那么这两个矩阵Rm和Rn之间的黎曼距离可表示如下:



其中,tr(·)表示矩阵的迹;
步骤六:用表示从第l个参考单元噪声数据通过正则burg递推计算得到的噪声统计协方差矩阵Rz,收集总共L个参考单元的噪声数据,计算噪声统计协方差矩阵...

【专利技术属性】
技术研发人员:姜宏周美含孙帅郭帛洋赵安琪张铭航曲良东莫秀玲
申请(专利权)人:吉林大学
类型:发明
国别省市:吉林;22

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1