CN105353367A - 一种双基地mimo雷达空间机动目标跟踪方法 - Google Patents
一种双基地mimo雷达空间机动目标跟踪方法 Download PDFInfo
- Publication number
- CN105353367A CN105353367A CN201510835721.9A CN201510835721A CN105353367A CN 105353367 A CN105353367 A CN 105353367A CN 201510835721 A CN201510835721 A CN 201510835721A CN 105353367 A CN105353367 A CN 105353367A
- Authority
- CN
- China
- Prior art keywords
- target
- matrix
- represent
- state
- sampling instant
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种双基地MIMO雷达的空间机动目标跟踪方法,能够实现对空间机动目标的三维位置、速度等运动参数的准确估计与预测;该方法将双基地MIMO雷达均匀圆阵体制应用于空间机动目标跟踪,在双基地MIMO雷达均匀圆阵角度多普勒频率估计的基础上,基于目标状态空间,在发射均匀圆阵与接收均匀圆阵的测量直角坐标系下,根据测元与目标运动状态变量之间的函数关系,进而建立测量方程,并采用适用于机动目标跟踪的“当前”统计模型建立描述目标运动参数随时间演化的状态方程,利用适用于非线性状态滤波的无迹卡尔曼滤波算法实时估计目标的运动参数,并可对目标运动参数进行预测,实现对空间机动目标的跟踪。
Description
技术领域
本发明属于雷达数据处理领域,具体涉及一种双基地MIMO雷达空间机动目标跟踪方法。
背景技术
双基地雷达是采用接收机和发射机***分置结构实现的。这种结构的主要特点是:发射机置于后方,接收机置于前方,可以接近目标进行隐蔽侦查,同时避免了雷达电磁波双程传播带来的威力损失,提高了目标的信噪比。传统的双基地米波雷达利用目标相对于接收端的角度及距离和对空间目标进行定位,由于双基地雷达的接收机与发射机难以满足精确的时间同步,且其角度分辨率与测距精度均比较低,导致对目标定位的精度低。采用MIMO(多输入多输出)技术的双基地雷达,即双基地MIMO雷达,可以在接收端获得发射端的角度信息,在不需要时间同步和目标距离的情况下,能够对目标进行精确定位,获得目标的位置坐标。
西安电子科技大学申请的专利技术“多输入多输出雷达***目标定位方法”(申请号200810150754.X,公开号CN101349748A)中公开了一种多输入多输出MIMO雷达目标定位的方法。该方法只能确定雷达目标的二维位置坐标,无法估计空间目标的三维位置坐标。
哈尔滨工程大学申请的专利技术“双基地多输入多输出雷达多目标定位方法”(申请号201110001351.0,公开号CN102135617A)中公开了一种双基地多输入多输出雷达多目标定位方法。该方法也仅适用于二维目标定位,无法估计空间目标的三维位置坐标。
西安电子科技大学申请的专利技术“双基地米波雷达目标三维精确定位方法”(专利申请号:201218001807.9)中公开了一种双基地多输入多输出雷达的多目标三维定位方法,将双基地雷达配置为发射阵为均匀圆阵、接收阵均匀线阵,估计出目标相对于均匀圆阵的方位角与俯仰角、目标相对于均匀线阵的接收角,利用这三个角度以及基线长度通过几何解算方法得到目标的三维位置坐标。该方法虽然可用于空间目标的三维定位,但仅利用了三个角度信息,且采用的是几何解算方法,当目标距离较远时,测角误差放大效应会使目标定位精度明显下降,无法对空间机动目标进行高精度跟踪。
西安电子科技大学申请的专利技术“双基地MIMO雷达均匀圆阵角度多普勒频率估计方法”(专利申请号:201510487104.4)中提出了一种双基地MIMO雷达均匀圆阵角度多普勒频率估计方法,将双基地雷达配置为发射阵、接收阵均为均匀圆阵,可同时估计出多个目标相对于接收阵与发射阵的方位角、俯仰角以及归一化多普勒频率。但该方法未涉及目标运动状态估计及预测跟踪问题。
发明内容
有鉴于此,本发明提供了一种双基地MIMO雷达空间机动目标的跟踪方法,其目的基于现有技术的不足,综合利用目标相对于双基地MIMO雷达的角度、多普勒频率估计信息,实现对空间机动目标的三维位置、速度等运动参数的准确估计与预测。
实施本发明的具体方案如下:
一种双基地MIMO雷达空间机动目标跟踪方法,包括以下步骤:
步骤1,分别将双基地MIMO雷达的发射机配置为M个阵元的均匀圆阵、接收机配置为N个阵元的均匀圆阵,并使发射机中M个阵元发射相互正交的波形信号;其中,M表示发射机阵元个数,N表示接收机阵元个数,且M、N均为自然数;
步骤2,利用所述发射机中M个阵元发射相互正交的波形信号,接收机中的N个阵元分别接收所述发射机中M个阵元发射的相互正交的波形信号,并进行匹配滤波,获得匹配滤波后的NM×1维雷达回波信号s和L次快拍积累得到的NM×L维雷达回波信号矩阵S,进而获得矩阵S中第n个接收阵元的M×L维切片矩阵形式Sn,其中,L表示快拍次数,n为自然数,n=1,2,...,N;
步骤3,根据N个接收阵元的M×L维切片矩阵S1~SN,利用平行因子算法分别获得发射方向向量的估计接收方向向量的估计和归一化多普勒频率方向向量的估计
步骤4,根据所述发射方向向量的估计获得目标相对于发射机的方位角估计值和俯仰角估计值根据所述接收方向向量的估计获得目标相对于接收机的方位角估计值和俯仰角估计值
步骤5,根据所述归一化多普勒频率方向向量的估计利用最小二乘算法获得目标的归一化多普勒频率估计值
步骤6,利用地心直角坐标系建立目标的运动状态空间,基于目标的运动状态空间,采用CS模型建立目标状态方程其中,为目标运动状态变量,包括目标三维位置[x,y,z]、速度加速度k表示目标状态的采样时刻序号,k=1,2,3,...,X(k-1)和X(k)分别表示第k-1个和第k个采样时刻的目标状态,Φ表示状态转移矩阵,U表示控制矩阵,表示第k个采样时刻的目标机动加速度均值,W(k)表示第k个采样时刻的目标状态噪声向量;
步骤7,基于目标状态空间,在发射均匀圆阵与接收均匀圆阵的测量直角坐标系下,根据测元与目标运动状态变量X之间的函数关系,建立测量方程Y(k)=h(X(k))+v(k),其中,v(k)表示各测元的测量噪声向量;Y(k)表示第k个采样时刻的目标测元,h(X(k))表示第k个采样时刻的目标测元与第k个采样时刻的目标运动状态之间的映射关系;
步骤8,基于测量方程,根据初始两个采样时刻获取的目标测元Y(1),Y(2),计算X(1)=[x(1),y(1),z(1)]T和X(2)=[x(2),y(2),z(2)]T,进而计算目标运动状态变量的初值和对应的协方差矩阵P(0),得到目标运动状态的滤波初值;
步骤9,获得滤波初值后,采用UKF算法,利用第k个采样时刻的测元Y(k)、目标在第k-1个采样时刻的状态估计与协方差矩阵P(k-1)获得目标在第k个采样时刻的状态估计与协方差矩阵P(k),依次获得目标在各采样时刻的运动状态,完成对空间机动目标的跟踪。
进一步地,在L次快拍积累中第n个接收阵元的M×L维切片矩阵Sn,其表达式为:
Sn=AR(n)ATBT+Wn,n=1,2,...,N
其中,AR(n)表示接收方向向量AR的第n个元素,Wn表示第n个接收阵元L次快拍积累中的噪声矩阵;AT表示发射方向向量,B表示归一化多普勒频率方向向量,其中,ρ为目标雷达波反射系数,fd为目标的归一化多普勒频率。
进一步地,在步骤3中,所述发射方向向量的估计所述接收方向向量的估计和所述归一化多普勒频率方向向量的估计的具体子步骤为:
a)计算第m个发射阵元的第n个接收阵元的第l次快拍的平行因子三线性模型形式sm,n,l,其具体表达式为
sm,n,l=AT(m)AR(n)B(l)+wm,n,l,m∈{1,…,M},n∈{1,…,N},l∈{1,…,L}
其中,AT(m)表示发射方向向量AT的第m个元素,AR(n)表示接收方向向量AR的第n个元素,B(l)表示归一化多普勒频率方向向量B的第l个元素,wm,n,l表示sm,n,l的三维噪声数据集;
b)根据第m个发射阵元的第n接收阵元的第l次快拍的平行因子三线性模型形式sm,n,l的对称性,对sm,n,l的第二维和第三维分别进行切片,分别获得第m个发射阵元的L×N维切片矩阵形式Fm、第l次快拍的N×M维切片矩阵形式Zl,其表达式分别为:
Fm=BAT(m)AR T+Wm,m∈{1,2,...,M}
Zl=ARB(l)AT T+Wl,l∈{1,2,...,L}
其中,Wm表示第m个发射阵元噪声,Wl表示第l次快拍噪声;
c)将N个接收阵元的M×L维切片矩阵形式S1~SN按列形成NM×L维矩阵S,将M个发射阵元的L×N维切片矩阵形式F1~FM按列形成ML×N维矩阵F,将L次快拍的N×M维切片矩阵形式Z1~ZL按列形成NL×M维矩阵Z,具体地,
S=[ARοAT]BT+WS=[S1,...,SN]T
F=[ATοB]AR T+WF=[F1,…,FM]T
Z=[BοAR]AT T+WZ=[Z1,…,ZL]T
其中,ο表示Khatri-Rao积,WS表示NM×L维矩阵S的噪声矩阵,WF表示ML×N维矩阵F的噪声矩阵,WZ表示NL×M维矩阵Z的噪声矩阵;
d)根据NM×L维矩阵S、ML×N维矩阵F和NL×M维矩阵Z,利用平行因子算法得到发射方向向量的估计接收方向向量的估计和归一化多普勒频率方向向量的估计具体地,
其中,[·]+表示取伪逆。
进一步地,在步骤8中,完成滤波初值的计算,具体为:
a)利用Y(1)中的4个角度计算目标在第1个采样时刻的位置参数X(1)=[x(1),y(1),z(1)]T;利用Y(2)中的4个角度计算目标在第2个采样时刻的位置参数X(2)=[x(2),y(2),z(2)]T;
根据4个角度θt,φt,θr,φr解算目标位置x,y,z的具体过程如下:
θt表示目标相对于发射机的方位角,φt表示目标相对于发射机的俯仰角,θr表示目标相对于接收机的方位角,φr表示目标相对于接收机的俯仰角;
i)将4个角度从发射阵与接收阵测量直角坐标系转至测量地心直角系为:
其中,Ωt T表示地心直角系到发射阵测量直角系的转换矩阵的转置,Ωr T表示地心直角系到接收阵测量直角系的转换矩阵的转置,(lt,mt,nt)为发射阵中心至目标连线的方向向量,θ′t、φ′t分别表示测量地心直角系下目标相对于发射机的方位角和俯仰角;(lr,mr,nr)为接收阵中心至目标连线的方向向量,θ′r、φ′r分别表示测量地心直角系下目标相对于接收机的方位角和俯仰角;
ii)采用水平投影法解算目标在地心直角系下的位置
其中
其中,(Xt0,Yt0,Zt0)为发射阵中心在地心直角坐标系下的位置;ΔX1为发射阵中心至目标在X轴方向的距离,Xr0为接收阵中心在地心直角坐标系X轴方向的位置,Zr0为接收阵中心在地心直角坐标系Z轴方向的位置;当k=1时,获得采样时刻的位置参数X(1)=[x(1),y(1),z(1)]T,当k=2时,获得采样时刻的位置参数X(2)=[x(2),y(2),z(2)]T;k取自然数;
b)对前两个采样时刻的位置参数进行差分计算目标在第2个采样时刻的近似速度
其中,T表示采样周期;
c)根据上述计算结果构造滤波初值,
P(0)=106I9
其中,I9表示9维的单位矩阵。
进一步地,在步骤9中,所述完成对空间机动目标跟踪的具体步骤为:
a)状态向量采样
对状态估计向量进行sigma采样,k取自然数,χp(k-1)表示第p个状态采样点,p=0,…,2nx,nx=9;采样点为:
其中,κ为比例参数,表示矩阵的第i列,i=1,2,…,nx
b)时间更新
利用CS模型对目标状态采样点进行一步预测得到χp(k|k-1),一步预测状态转移矩阵为Φa,根据目标测元与目标运动状态之间的映射函数h(·)计算观测一步预测的采样点Υp(k|k-1),再对其加权求和,计算观测量的一步预测
χp(k|k-1)表示第p个状态采样点从k-1时刻到k时刻的预测值;Υp(k|k-1)表示第p个观测一步预测的采样点;表示第p个采样点的均值权重;为第k个采样时刻的观测量的一步预测;
χp(k/k-1)=Φaχp(k-1)
Υp(k|k-1)=h(χp(k|k-1))
c)状态更新
P(k)=P(k|k-1)-G(k)Pyy(k|k-1)GT(k)
G(k)表示第k个采样时刻目标状态更新增益矩阵;P(k|k-1)表示目标状态协方差矩阵的一步预测;Pyy(k|k-1)表示观测协方差矩阵的一步预测;和分别表示第k-1到第k个采样时刻的目标状态一步预测和测元一步预测的均值;GT(k)表示第k个采样时刻的目标状态更新增益矩阵的转置。
有益效果:
(1)本发明采用发射机、接收机均为均匀圆阵的配置,不但能够估计目标的四个角度和一个多普勒频率信息,还能够对空间目标进行三维定位。
(2)在进行目标运动状态参数解算时,传统的几何解算方法只能利用角度测元,无法利用多普勒频率测元,造成了测元信息浪费,而本发明采用了基于CS模型的UKF算法估计目标的运动状态参数,一方面提高了对测量噪声的抑制能力,另一方面能够融合利用所有测元信息,这两方面共同提高了目标运动状态参数的解算精度。
(3)本发明采用了适用于机动目标跟踪的CS模型,基于此可以对空间机动目标的运动状态进行一定时间的预测,实现对空间机动目标的有效稳定跟踪。
附图说明
图1为本发明的一种双基地MIMO雷达空间机动目标跟踪方法实现流程图
图2为本发明的双基地MIMO雷达配置示意图。
图3为用本发明方法仿真信噪比20dB情况下某临近空间运动目标相对于发射机的俯仰角和方位角的估计误差图。
图4为用本发明方法仿真信噪比20dB情况下某临近空间运动目标相对于发射机的俯仰角和方位角的估计误差图。
图5为用本发明方法仿真信噪比20dB情况下某临近空间运动目标多普勒频率的估计误差图。
图6为用本发明方法仿真信噪比为20dB时的目标位置滤波误差图。
图7为用本发明方法仿真信噪比为20dB时的目标速度滤波误差图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
图1为本发明的一种双基地MIMO雷达空间机动目标跟踪方法实现流程图,实现本发明的具体步骤如下:
步骤1,分别将双基地MIMO雷达的发射机配置为M个阵元的均匀圆阵、接收机配置为N个阵元的均匀圆阵,并使发射机中M个阵元发射相互正交的波形信号;其中,M表示发射机阵元个数,N表示接收机阵元个数,且M、N均为自然数;
图2为本发明的双基地MIMO雷达配置示意图。坐标系采用地心直角坐标系OG-XGYGZG,原点OG为地球中心,OGXG轴的正向为沿起始天文子午面与地球赤道平面的交线指向外方向,OGZG轴沿地球自转轴指向北极,OGYG轴与OGXG、OGZG轴构成右手直角坐标系。其中,图2中Ot为发射圆阵的中心点,Or为接收圆阵的中心点。
发射圆阵与接收圆阵均布设在当地水平面内。设发射圆阵的圆点Ot的大地坐标(经度Lt,纬度Bt,高程ht)为(Lt,Bt,ht),接收圆阵的圆点Or的大地坐标(经度Lr,纬度Br,高程hr)为(Lr,Br,hr)。发射圆阵半径为rt,其阵元个数M=2floor(2πrt/λ)+1,floor(·)表示向下取整运算,λ表示发射波的波长;接收圆阵半径为rr,其阵元个数N=2floor(2πrr/λ)+1。
假设空间有一目标P,在地心系下的位置坐标为(x,y,z)、速度坐标为发射圆阵的圆心与目标的连线OtP与发射圆阵平面的夹角为俯仰角,记为φt,OtP在发射圆阵平面的投影与当地正北方向的夹角为方位角,记为θt;接收圆阵的圆心至目标连线OrP与接收圆阵平面的夹角为俯仰角,记为φr,OrP在接收圆阵平面的投影与当地正北方向的夹角为方位角,记为θr;由目标运动引起的接收信号多普勒频率为fd。
步骤2,利用所述发射机中M个阵元发射相互正交的波形信号,接收机中的N个阵元分别接收所述发射机中M个阵元发射的相互正交的波形信号,并进行匹配滤波,获得匹配滤波后的NM×1维雷达回波信号s和L次快拍积累得到的NM×L维雷达回波信号矩阵S,进而获得矩阵S中第n个接收阵元的M×L维切片矩阵形式Sn,其中,L表示快拍次数,n为自然数,n=1,2,...,N;
具体地,利用发射机中M个阵元发射相互正交的波形信号,接收机中的N个阵元分别接收该发射机中M个阵元的发射信号,并进行匹配滤波,依次获得匹配滤波后的NM×1维雷达回波信号s,考察第l个脉冲匹配滤波后的回波信号s(l),其表达式为
其中,ρ为目标雷达波反射系数,假定目标对不同正交发射信号的反射系数在相同脉冲内是完全相关的,在不同脉冲间是相互独立的,ar(θr,φr)为目标接收方向向量,at(θt,φt)为目标发射方向向量,表示Kronecker积,为第l个脉冲匹配滤波后的接收噪声,l=1,2,…,L;
其中,βrn为第n个接收阵元的方位角,βrn=2π(n-1)/N,n=1,2,…,N,符号表示“简记为”;arn(θr,φr)表示第n个接收阵元的接收方向因子。
其中,βtm为第m个发射阵元的方位角,βtm=2π(m-1)/M,m=1,2,…,M;αtm(θt,φt)表示第m个发射阵元的的发射方向因子。
经过L次快拍积累得到NM×L维矩阵S,在L次快拍内认为目标运动状态参数不变,S的表达式为
其中, 上标T表示转置。表示由形成的矩阵,B表示归一化多普勒频率方向向量。
进而得到L次快拍积累中第n个接收阵元的M×L维切片矩阵Sn,其表达式为:
Sn=AR(n)ATBT+Wn,n=1,2,...,N
其中,AR(n)表示向量AR的第n个元素,Wn表示第n个接收阵元L次快拍积累中的噪声矩阵;AT表示发射方向向量,AR表示接收方向向量。
步骤3,根据N个接收阵元的M×L维切片矩阵形式S1~SN,利用平行因子算法分别获得发射方向向量的估计接收方向向量的估计和多普勒频率方向向量的估计
步骤3的具体子步骤为:
3a)第n个接收阵元的M×L维切片矩阵Sn,根据N个接收阵元的M×L维切片矩阵S1~SN,获得M×N×L的三维数据集,进而得到第m个发射阵元的第n个接收阵元的第l次快拍的平行因子三线性模型形式sm,n,l,其具体表达式为
sm,n,l=AT(m)AR(n)B(l)+wm,n,l,m∈{1,…,M},n∈{1,…,N},l∈{1,…,L},
其中,AT(m)表示发射方向向量AT的第m个元素,AR(n)表示接收方向向量AR的第n个元素,B(l)表示多普勒频率方向向量B的第l个元素,wm,n,l表示三维噪声数据集;
3b)根据sm,n,l,分别得到第m个发射阵元的L×N维切片矩阵形式Fm、第l次快拍的N×M维切片矩阵形式Zl,进而分别获得M个发射阵元的L×N维切片矩阵形式F1~FM,L次快拍的N×M维切片矩阵形式Z1~ZL。
具体地,根据第m个发射阵元的第n接收阵元的第l次快拍的平行因子三线性模型形式sm,n,l的对称性,对其第二维和第三维分别进行切片,分别获得第m个发射阵元的L×N维切片矩阵形式Fm、第l次快拍的N×M维切片矩阵形式Zl,其表达式分别为:
Fm=BAT(m)AR T+Wm,m∈{1,2,...,M}
Zl=ARB(l)AT T+Wl,l∈{1,2,...,L}
3c)将N个接收阵元的M×L维切片矩阵形式S1~SN按列形成NM×L维矩阵S,将M个发射阵元的L×N维切片矩阵形式F1~FM按列形成ML×N维矩阵F,将L次快拍的N×M维切片矩阵形式Z1~ZL按列形成NL×M维矩阵Z,具体地,
S=[ARοAT]BT+WS=[S1,...,SN]T
F=[ATοB]AR T+WF=[F1,…,FM]T
Z=[BοAR]AT T+WP=[Z1,...,ZL]T
其中,ο表示Khatri-Rao积,WS表示NM×L维矩阵S的噪声矩阵,WF表示ML×N维矩阵F的噪声矩阵,WP表示NL×M维矩阵P的噪声矩阵。
3d)根据NM×L维矩阵S、ML×N维矩阵F和NL×M维矩阵P,利用平行因子算法得到发射方向向量的估计接收方向向量的估计和多普勒频率方向向量的估计具体地,
其中,[·]+表示取伪逆。
步骤4,根据所述发射方向向量的估计获得目标相对于发射机的方位角估计值和俯仰角估计值根据所述接收方向向量的估计获得目标相对于接收机的方位角估计值和俯仰角估计值
具体地,发射方向向量的估计为接收方向向量的估计的为其表达式分别为:
一般,βr1=0,βt1=0,中的每一项都除以第一项然后去掉其第一项,获得发射新向量a1,再取a1对数的虚部得到a1=Im(lna1);中的每一项都除以第一项然后去掉其第一项,得到接收新向量a2,再取a2对数的虚部得到a′2=Im(lna2),a′1和a′2的表达式分别为:
其中,ξt=2πrt/λ,ξr=2πrr/λ;
a′1中的第m项除以(cosβtm′-1),m′∈{2,3,...,M},得到β1,
a′2中的第n项除以(cosβrn′-1),n′∈{2,3,...,N},得到β2,则β1和β2的表达式分别为:
其中,ct0=ξtsinφtcosθt,ct1=ξtsinφtsinθt,cr0=ξrsinφrcosθr,cr1=ξrsinφrsinθr;
根据 可知 的求解是一个标准的参数估计问题,可以用最小二乘估计得到 的估计值
其中,
根据 可知 的求解是一个标准的参数估计问题,可以用最小二乘估计,得到 的估计值
其中,
进而获得雷达目标的发射方位角估计值发射俯仰角估计值接收方位角估计值和接收俯仰角估计值分别为 发射方位角估计值发射俯仰角估计值的估计误差如图3所示,接收方位角估计值和接收俯仰角估计值的估计误差如图4所示。
步骤5,根据所述归一化多普勒频率方向向量的估计利用最小二乘算法获得目标的归一化多普勒频率估计值
具体地,根据归一化多普勒频率方向向量令B的每一项都除以第一项ρ,得到再取对数的虚部,得到其相位
其中,ρ表示目标的雷达波发射系数,angle(·)表示取相位。
根据 可知 均求解是一个标准的参数估计问题,可以用最小二乘估计,得到其估计
其中,
归一化多普勒频率估计值的估计误差如图5所示。
步骤6,利用地心直角坐标系建立目标的运动状态空间,基于目标的运动状态空间,采用CS模型建立目标状态方程其中,为目标运动状态变量,包括目标在地心直角坐标系下的位置[x,y,z]、速度加速度k表示目标状态的采样时刻序号,k=1,2,3,...,X(k-1)和X(k)分别表示第k-1和第k个采样时刻的目标状态,Φ表示状态转移矩阵,U表示控制矩阵,表示第k个采样时刻的目标机动加速度均值,W(k)表示第k个采样时刻的目标状态噪声向量;
具体地,
其中,表示p1方向上的正向极限加速度,表示p1方向上的负向极限加速度,p1=x,y,z,表示p1方向上的加速度均值。
q12=0.5α-4(e-2αT+1-2e-αT+2αTe-αT-2αT+α2T2)
q13=0.5α-3(1-e-2αT-2αTe-αT)
q22=0.5α-3(4e-αT-3-e-2αT+2αT)
q23=0.5α-2(e-2αT+1-2e-αT)
q33=0.5α-1(1-e-2αT)
步骤7,基于目标状态空间,在发射圆阵与接收圆阵的测量坐标系下,根据测元与目标运动状态变量X之间的函数关系建立测量方程Y(k)=h(X(k))+v(k),其中,v(k)表示各测元测量噪声构成的向量;
其中,v(k)=[ε1,ε2,ε3,ε4,ε5]T为测量噪声向量,ε1,ε2,ε3,ε4,ε5分别为对应的噪声,满足测量噪声协方差矩阵 diag(·)表示取向量构成的对角矩阵,分别表示各种测元的测量噪声方差;Δxt,Δyt,Δzt,vt分别表示目标在发射圆阵测量直角坐标系下的三维位置坐标和径向速率,Δxr,Δyr,Δzr,vr分别表示目标在接收圆阵测量直角坐标系下的三维位置坐标和径向速率,c表示空气介质中的光速,f0表示雷达载波频率,PRF表示雷达脉冲重复频率。
其中,Ωt表示地心直角系到发射阵测量直角系的转换矩阵,Ωr表示地心直角系到接收阵测量直角系的转换矩阵,Rx(α),Ry(α),Rz(α)分别表示绕x,y,z轴逆时针旋转α角的旋转矩阵,具体为
其中,表示大地坐标到地心直角坐标的转换函数,的具体形式为
其中a为地球参考椭球半长轴,e2为参考椭球第一偏心率的平方。
步骤8,基于测量方程,根据初始两个采样时刻获取的目标测元Y(1),Y(2),计算X(1)=[x(1),y(1),z(1)]T和X(2)=[x(2),y(2),z(2)]T,进而计算目标运动状态变量的初值和对应的协方差矩阵P(0),得到滤波初值;具体为:
8a)利用Y(1)中的4个角度计算目标在第1个采样时刻的位置参数x(1),y(1),,z(1);利用Y(2)中的4个角度计算目标在第2个采样时刻的位置参数x(2),y(2),z(2);
根据4个角度θt,φt,θr,φr解算目标位置x,y,z的具体过程如下:
i)将4个角度从发射阵与接收阵的测量直角坐标系转至测量地心直角系下(其坐标原点与测量直角坐标系的原点相同,各坐标轴与地心直角坐标系平行)
其中,Ωt T表示地心直角系到发射阵测量直角系的转换矩阵的转置,Ωr T表示地心直角系到接收阵测量直角系的转换矩阵的转置,(lt,mt,nt)为发射阵中心至目标连线的方向向量,θ′t、φ′t分别表示测量地心直角系下目标相对于发射机的方位角和俯仰角;(lr,mr,nr)为接收阵中心至目标连线的方向向量,θ′r、φ′r分别表示测量地心直角系下目标相对于接收机的方位角和俯仰角;
ii)采用水平投影法解算目标在地心直角系下的位置
其中
其中,(Xt0,Yt0,Zt0)为发射阵中心在地心直角坐标系下的位置;ΔX1为发射阵中心至目标在X轴方向的距离,Xr0为接收阵中心在地心直角坐标系X轴方向的位置,Zr0为接收阵中心在地心直角坐标系Z轴方向的位置;当k=1时,获得采样时刻的位置参数X(1)=[x(1),y(1),z(1)]T,当k=2时,获得采样时刻的位置参数X(2)=[x(2),y(2),z(2)]T;k取自然数;
8b)对前两个采样时刻的位置参数进行差分计算目标在第2个采样时刻的近似速度
其中,T表示采样周期;
8c)根据计算结果构造滤波初值,
P(0)=106I9
这里I9表示9维的单位阵。
步骤9,在获得滤波初值后,采用UKF算法,利用第k个采样时刻的测元Y(k)、目标在第k-1个采样时刻的状态估计与协方差矩阵P(k-1)获得目标在第k个采样时刻的状态估计与协方差矩阵P(k),依次获得目标在各采样时刻的运动状态,完成对空间机动目标的跟踪。
9a)状态向量采样
对状态向量进行sigma点采样,k取自然数,χp(k-1)表示第p个状态采样点,p=0,…,2nx,nx=9,公别表示第p个采样点的均值权重与方差权重,κ为设计参数:
其中,表示矩阵的第i列,i=1,2,…,nx
采样后的向量:
9b)时间更新
利用CS模型对目标状态采样点进行一步预测,采用加速度均值自适应算法,有
χp(k|k-1)=Φaχp(k-1)
式中χp(k|k-1)表示第p个采样点χp(k-1)从k-1时刻到k时刻的预测值。
表示p1方向上机动加速度的一步预测,p1=x,y,z。
根据状态Sigma采样点的一步预测,利用测元与目标状态变量之间的映射函数h(·)计算观测一步预测的Sigma点,Υp(k|k-1)表示观测一步预测的第p个Sigma点,再对其加权求和,计算观测量的一步预测
Υp(k|k-1)=h(χp(k|k-1))
9c)状态更新
P(k)=P(k|k-1)-G(k)Pyy(k|k-1)GT(k)
G(k)表示第k个采样时刻目标状态更新增益矩阵;Pyy(k|k-1)表示观测协方差矩阵的一步预测,Pyy(k|k-1)表示状态观测互协方差矩阵的一步预测。和分别表示第k-1到第k个采样时刻的目标状态一步预测和测元一步预测的均值;GT(k)表示第k个采样时刻的目标状态更新增益矩阵的转置。
上述式中Q为状态噪声阵,R为测量噪声阵,κ为比例参数,用于调节sigma点和的距离,仅影响二阶之后的高阶矩带来的偏差,一般取κ=0;α为正值的比例缩放因子,控制sigma点分布的范围,通常取(0,1)的一个较小的值;β为引入f(·)高阶项信息的参数,通常取β=2。目标位置跟踪误差如图6所示和目标速度跟踪误差如图7所示。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种双基地MIMO雷达空间机动目标跟踪方法,其特征在于,包括以下步骤:
步骤1,分别将双基地MIMO雷达的发射机配置为M个阵元的均匀圆阵、接收机配置为N个阵元的均匀圆阵,并使发射机中M个阵元发射相互正交的波形信号;其中,M表示发射机阵元个数,N表示接收机阵元个数,且M、N均为自然数;
步骤2,利用所述发射机中M个阵元发射相互正交的波形信号,接收机中的N个阵元分别接收所述发射机中M个阵元发射的相互正交的波形信号,并进行匹配滤波,获得匹配滤波后的NM×1维雷达回波信号s和L次快拍积累得到的NM×L维雷达回波信号矩阵S,进而获得矩阵S中第n个接收阵元的M×L维切片矩阵形式Sn,其中,L表示快拍次数,n为自然数,n=1,2,...,,N;
步骤3,根据N个接收阵元的M×L维切片矩阵S1~SN,利用平行因子算法分别获得发射方向向量的估计接收方向向量的估计和归一化多普勒频率方向向量的估计
步骤4,根据所述发射方向向量的估计获得目标相对于发射机的方位角估计值和俯仰角估计值根据所述接收方向向量的估计获得目标相对于接收机的方位角估计值和俯仰角估计值
步骤5,根据所述归一化多普勒频率方向向量的估计利用最小二乘算法获得目标的归一化多普勒频率估计值
步骤6,利用地心直角坐标系建立目标的运动状态空间,基于目标的运动状态空间,采用CS模型建立目标状态方程其中,为目标运动状态变量,包括目标三维位置[x,y,z]、速度加速度k表示目标状态的采样时刻序号,k=1,2,3,…,X(k-1)和X(k)分别表示第k-1个和第k个采样时刻的目标状态,Φ表示状态转移矩阵,U表示控制矩阵,表示第k个采样时刻的目标机动加速度均值,W(k)表示第k个采样时刻的目标状态噪声向量;
步骤7,基于目标状态空间,在发射均匀圆阵与接收均匀圆阵的测量直角坐标系下,根据测元与目标运动状态变量X之间的函数关系,建立测量方程Y(k)=h(X(k))+v(k),其中,v(k)表示各测元的测量噪声向量;Y(k)表示第k个采样时刻的目标测元,h(X(k))表示第k个采样时刻的目标测元与第k个采样时刻的目标运动状态之间的映射关系;
步骤8,基于测量方程,根据初始两个采样时刻获取的目标测元Y(1),Y(2),计算X(1)=[x(1),y(1),z(1)]T和X(2)=[x(2),y(2),z(2)]T,进而计算目标运动状态变量的初值和对应的协方差矩阵P(0),得到目标运动状态的滤波初值;
步骤9,获得滤波初值后,采用UKF算法,利用第k个采样时刻的测元Y(k)、目标在第k-1个采样时刻的状态估计与协方差矩阵P(k一1)获得目标在第k个采样时刻的状态估计与协方差矩阵P(k),依次获得目标在各采样时刻的运动状态,完成对空间机动目标的跟踪。
2.如权利要求1所述一种双基地MIMO雷达空间机动目标跟踪方法,其特征在于,在L次快拍积累中第n个接收阵元的M×L维切片矩阵Sn,其表达式为:
Sn=AR(n)ATBT+Wn,n=1,2,...,N
其中,AR(n)表示接收方向向量AR的第n个元素,Wn表示第n个接收阵元L次快拍积累中的噪声矩阵;AT表示发射方向向量,B表示归一化多普勒频率方向向量,其中,ρ为目标雷达波反射系数,fd为目标的归一化多普勒频率。
3.如权利要求1所述一种双基地MIMO雷达空间机动目标跟踪方法,其特征在于,在步骤3中,所述发射方向向量的估计所述接收方向向量的估计和所述归一化多普勒频率方向向量的估计的具体子步骤为:
a)计算第m个发射阵元的第n个接收阵元的第l次快拍的平行因子三线性模型形式sm,n,l,其具体表达式为
sm,n,l=AT(m)AR(n)B(l)+wm,n,l,m∈{1,…,M},n∈{1,…,N},l∈{1,…,L}
其中,AT(m)表示发射方向向量AT的第m个元素,AR(n)表示接收方向向量AR的第n个元素,B(l)表示归一化多普勒频率方向向量B的第l个元素,wm,n,l表示sm,n,l的三维噪声数据集;
b)根据第m个发射阵元的第n接收阵元的第l次快拍的平行因子三线性模型形式sm,n,l的对称性,对sm,n,l的第二维和第三维分别进行切片,分别获得第m个发射阵元的L×N维切片矩阵形式Fm、第l次快拍的N×M维切片矩阵形式Zl,其表达式分别为:
Fm=BAT(m)AR T+Wm,m∈{1,2,…,M}
Zl=ARB(l)AT T+Wl,l∈{1,2,...,L}
其中,Wm表示第m个发射阵元噪声,Wl表示第l次快拍噪声;
c)将N个接收阵元的M×L维切片矩阵形式S1~SN按列形成NM×L维矩阵S,将M个发射阵元的L×N维切片矩阵形式F1~FM按列形成ML×N维矩阵F,将L次快拍的N×M维切片矩阵形式Z1~ZL按列形成NL×M维矩阵Z,具体地,
S=[ARοAT]BT+WS=[S1,…,SN]T
F=[ATοB]AR T+WF=[F1,…,FM]T
Z=[BοAR]AT T+WZ=[Z1,…,ZL]T
其中,ο表示Khatri-Rao积,WS表示NM×L维矩阵S的噪声矩阵,WF表示ML×N维矩阵F的噪声矩阵,WZ表示NL×M维矩阵Z的噪声矩阵;
d)根据NM×L维矩阵S、ML×N维矩阵F和NL×M维矩阵Z,利用平行因子算法得到发射方向向量的估计接收方向向量的估计和归一化多普勒频率方向向量的估计具体地,
其中,[·]+表示取伪逆。
4.如权利要求1所述一种双基地MIMO雷达空间机动目标跟踪方法,其特征在于,在步骤8中,完成滤波初值的计算,具体为:
a)利用Y(1)中的4个角度计算目标在第1个采样时刻的位置参数X(1)=[x(1),y(1),z(1)]T:利用Y(2)中的4个角度计算目标在第2个采样时刻的位置参数X(2)=[x(2),y(2),z(2)]T;
根据4个角度θt,φt,θr,φr解算目标位置x,y,z的具体过程如下:
θt表示目标相对于发射机的方位角,φt表示目标相对于发射机的俯仰角,θr表示目标相对于接收机的方位角,φr表示目标相对于接收机的俯仰角;
i)将4个角度从发射阵与接收阵测量直角坐标系转至测量地心直角系为:
其中,Ωt T表示地心直角系到发射阵测量直角系的转换矩阵的转置,Ωr T表示地心直角系到接收阵测量直角系的转换矩阵的转置,(lt,mt,nt)为发射阵中心至目标连线的方向向量,θ′t、φ′t分别表示测量地心直角系下目标相对于发射机的方位角和俯仰角;(lr,mr,nr)为接收阵中心至目标连线的方向向量,θ′r、φ′r分别表示测量地心直角系下目标相对于接收机的方位角和俯仰角;
ii)采用水平投影法解算目标在地心直角系下的位置
其中
其中,(Xt0,Yt0,Zt0)为发射阵中心在地心直角坐标系下的位置;ΔX1为发射阵中心至目标在X轴方向的距离,Xr0为接收阵中心在地心直角坐标系X轴方向的位置,Zr0为接收阵中心在地心直角坐标系Z轴方向的位置;当k=1时,获得采样时刻的位置参数X(1)=[x(1),y(1),z(1)]T,当k=2时,获得采样时刻的位置参数X(2)=[x(2),y(2),z(2)]T;k取自然数;
b)对前两个采样时刻的位置参数进行差分计算目标在第2个采样时刻的近似速度
其中,T表示采样周期;
c)根据上述计算结果构造滤波初值,
P(0)=106I9
其中,I9表示9维的单位矩阵。
5.如权利要求1所述一种双基地MIMO雷达空间机动目标跟踪方法,其特征在于,在步骤9中,所述完成对空间机动目标跟踪的具体步骤为:
a)状态向量采样
对状态估计向量进行sigma采样,k取自然数,χp(k-1)表示第p个状态采样点,p=0,…,2nx,nx=9;采样点为:
其中,κ为比例参数,表示矩阵的第i列,i=1,2,…,nx
b)时间更新
利用CS模型对目标状态采样点进行一步预测得到χp(k|k-1),一步预测状态转移矩阵为Φa,根据目标测元与目标运动状态之间的映射函数h(·)计算观测一步预测的采样点γp(k|k-1),再对其加权求和,计算观测量的一步预测
χp(k|k-1)表示第p个状态采样点从k-1时刻到k时刻的预测值;γp(k|k-1)表示第p个观测一步预测的采样点;表示第p个采样点的均值权重;为第k个采样时刻的观测量的一步预测;
χp(k/k-1)=Φaχp(k-1)
γp(k|k-1)=h(χp(k|k-1))
c)状态更新
P(k)=P(k|k-1)-G(k)Pyy(k|k-1)GT(k)
G(k)表示第k个采样时刻目标状态更新增益矩阵;P(k|k-1)表示目标状态协方差矩阵的一步预测;Pyy(k|k-1)表示观测协方差矩阵的一步预测;和分别表示第k-1到第k个采样时刻的目标状态一步预测和测元一步预测的均值;GT(k)表示第k个采样时刻的目标状态更新增益矩阵的转置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510835721.9A CN105353367A (zh) | 2015-11-26 | 2015-11-26 | 一种双基地mimo雷达空间机动目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510835721.9A CN105353367A (zh) | 2015-11-26 | 2015-11-26 | 一种双基地mimo雷达空间机动目标跟踪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105353367A true CN105353367A (zh) | 2016-02-24 |
Family
ID=55329368
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510835721.9A Pending CN105353367A (zh) | 2015-11-26 | 2015-11-26 | 一种双基地mimo雷达空间机动目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105353367A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107728138A (zh) * | 2017-09-15 | 2018-02-23 | 电子科技大学 | 一种基于当前统计模型的机动目标跟踪方法 |
CN108490425A (zh) * | 2018-03-07 | 2018-09-04 | 西安电子科技大学 | 一种双基地mimo雷达的测角方法 |
CN108646238A (zh) * | 2018-03-06 | 2018-10-12 | 中国船舶重工集团公司第七二四研究所 | 一种基于副瓣对消系数映射的干扰源跟踪方法 |
CN109143223A (zh) * | 2018-08-14 | 2019-01-04 | 中国电子科技集团公司第三十八研究所 | 一种双基地雷达的空间目标跟踪滤波装置及方法 |
CN109257085A (zh) * | 2018-03-30 | 2019-01-22 | 北京润科通用技术有限公司 | 人造卫星与飞行设备间多普勒频移的获得方法及装置 |
CN110031802A (zh) * | 2019-04-04 | 2019-07-19 | 中国科学院数学与***科学研究院 | 具有未知量测零偏的双红外传感器的融合定位方法 |
CN110115823A (zh) * | 2018-02-06 | 2019-08-13 | 英飞凌科技股份有限公司 | 跑步机和跑步机上的非接触式感测方法 |
WO2020056586A1 (zh) * | 2018-09-18 | 2020-03-26 | 深圳市大疆创新科技有限公司 | 高度确定方法、装置、电子设备和计算机可读存储介质 |
US11914070B2 (en) | 2020-05-29 | 2024-02-27 | Rohde & Schwarz Gmbh & Co. Kg | Radar target simulator front end and method for simulating |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102156279A (zh) * | 2011-05-12 | 2011-08-17 | 西安电子科技大学 | 基于mimo的双基地雷达地面动目标检测方法 |
EP2496958A1 (en) * | 2009-11-06 | 2012-09-12 | Saab AB | Radar system and method for detecting and tracking a target |
CN104020480A (zh) * | 2014-06-17 | 2014-09-03 | 北京理工大学 | 一种带自适应因子的交互式多模型ukf的卫星导航方法 |
CN104793201A (zh) * | 2015-05-04 | 2015-07-22 | 哈尔滨工业大学 | 一种跟踪临近空间高超声速目标的修正变结构网格交互多模型滤波方法 |
CN105068068A (zh) * | 2015-08-10 | 2015-11-18 | 西安电子科技大学 | 双基地mimo雷达均匀圆阵角度多普勒频率估计方法 |
-
2015
- 2015-11-26 CN CN201510835721.9A patent/CN105353367A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2496958A1 (en) * | 2009-11-06 | 2012-09-12 | Saab AB | Radar system and method for detecting and tracking a target |
CN102156279A (zh) * | 2011-05-12 | 2011-08-17 | 西安电子科技大学 | 基于mimo的双基地雷达地面动目标检测方法 |
CN104020480A (zh) * | 2014-06-17 | 2014-09-03 | 北京理工大学 | 一种带自适应因子的交互式多模型ukf的卫星导航方法 |
CN104793201A (zh) * | 2015-05-04 | 2015-07-22 | 哈尔滨工业大学 | 一种跟踪临近空间高超声速目标的修正变结构网格交互多模型滤波方法 |
CN105068068A (zh) * | 2015-08-10 | 2015-11-18 | 西安电子科技大学 | 双基地mimo雷达均匀圆阵角度多普勒频率估计方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107728138B (zh) * | 2017-09-15 | 2020-11-17 | 电子科技大学 | 一种基于当前统计模型的机动目标跟踪方法 |
CN107728138A (zh) * | 2017-09-15 | 2018-02-23 | 电子科技大学 | 一种基于当前统计模型的机动目标跟踪方法 |
CN110115823A (zh) * | 2018-02-06 | 2019-08-13 | 英飞凌科技股份有限公司 | 跑步机和跑步机上的非接触式感测方法 |
CN110115823B (zh) * | 2018-02-06 | 2022-03-08 | 英飞凌科技股份有限公司 | 跑步机和跑步机上的非接触式感测方法 |
CN108646238A (zh) * | 2018-03-06 | 2018-10-12 | 中国船舶重工集团公司第七二四研究所 | 一种基于副瓣对消系数映射的干扰源跟踪方法 |
CN108490425A (zh) * | 2018-03-07 | 2018-09-04 | 西安电子科技大学 | 一种双基地mimo雷达的测角方法 |
CN109257085B (zh) * | 2018-03-30 | 2020-11-03 | 北京润科通用技术有限公司 | 人造卫星与飞行设备间多普勒频移的获得方法及装置 |
CN109257085A (zh) * | 2018-03-30 | 2019-01-22 | 北京润科通用技术有限公司 | 人造卫星与飞行设备间多普勒频移的获得方法及装置 |
CN109143223A (zh) * | 2018-08-14 | 2019-01-04 | 中国电子科技集团公司第三十八研究所 | 一种双基地雷达的空间目标跟踪滤波装置及方法 |
CN109143223B (zh) * | 2018-08-14 | 2023-03-07 | 中国电子科技集团公司第三十八研究所 | 一种双基地雷达的空间目标跟踪滤波装置及方法 |
WO2020056586A1 (zh) * | 2018-09-18 | 2020-03-26 | 深圳市大疆创新科技有限公司 | 高度确定方法、装置、电子设备和计算机可读存储介质 |
CN110031802B (zh) * | 2019-04-04 | 2020-10-09 | 中国科学院数学与***科学研究院 | 具有未知量测零偏的双红外传感器的融合定位方法 |
CN110031802A (zh) * | 2019-04-04 | 2019-07-19 | 中国科学院数学与***科学研究院 | 具有未知量测零偏的双红外传感器的融合定位方法 |
US11914070B2 (en) | 2020-05-29 | 2024-02-27 | Rohde & Schwarz Gmbh & Co. Kg | Radar target simulator front end and method for simulating |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105353367A (zh) | 一种双基地mimo雷达空间机动目标跟踪方法 | |
CN102981144B (zh) | 空中运动平台对目标的三维无源定位方法 | |
CN105445730B (zh) | 一种基于角度分集的海洋流场反演星载sar***及其方法 | |
CN106353744A (zh) | 基于双基地fda‑mimo雷达的多参数联合估计方法 | |
CN101915911B (zh) | 基于相消积累空时谱的空间任意构型分布式sar动目标参数估计方法 | |
CN104698453B (zh) | 基于合成孔径天线阵列的雷达信号被动定位方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN107526089B (zh) | 一种基于时延二次差分的非共视雷达信号无源定位方法 | |
CN105044667A (zh) | 一种运动目标的双星跟踪方法、装置和*** | |
US20220082707A1 (en) | Techniques for Determining Geolocations | |
US20210011109A1 (en) | Method and apparatus for determining the direction of arrival of radio or acoustic signals, and for transmitting directional radio or acoustic signals | |
CN103576137A (zh) | 一种基于成像策略的多传感器多目标定位方法 | |
CN105974386A (zh) | 一种多基地雷达多目标成像定位方法 | |
CN108872971A (zh) | 一种基于运动单阵列的目标定位方法与装置 | |
CN105353355A (zh) | 一种基于稀疏重构和投影成像的多基地雷达多目标定位方法 | |
CN104678418A (zh) | 一种基于多星gnss-r海面目标定位模糊消除方法 | |
CN104267389A (zh) | 一种mimo天波超视距雷达信号处理方法 | |
CN104678386A (zh) | 一种利用gnss海面反射信号相关功率探测目标的方法 | |
CN104267420A (zh) | 一种星载对运动目标的三维定位方法、装置和*** | |
CN104820221B (zh) | 多基合成孔径雷达的目标三维定位方法 | |
Luo et al. | Target location and height estimation via multipath signal and 2D array for sky-wave over-the-horizon radar | |
CN105372652A (zh) | 基于接收线阵的mimo雷达空间机动目标跟踪方法 | |
CN105548959A (zh) | 一种基于稀疏重建的多传感器多目标的定位方法 | |
CN106569189A (zh) | 一种高分辨率sar卫星距离模糊度性能分析方法 | |
Fois et al. | DOPSCAT: A mission concept for a Doppler wind-scatterometer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160224 |
|
WD01 | Invention patent application deemed withdrawn after publication |