运动观测站误差下外辐射源雷达时差与频差协同定位方法与流程

本发明属于雷达数据处理领域，具体涉及一种运动观测站误差下外辐射源雷达tdoa/fdoa协同定位方法。

背景技术：

外辐射源雷达不主动发射信号，依靠第三方辐射源(如电视广播信号、电视信号、手机信号、卫星导航信号等)探测目标，具有隐蔽性好、探测低空与超低空目标、低成本、抗干扰能力强等特点。外辐射源雷达接收站(观测站)接收发射源直达波信号和经目标散射的回波信号，通过信号处理得到目标到达角度(doa)、到达时差(tdoa)和到达频差(fdoa)的量测信息。外辐射源雷达系统作为一种双/多基地结构的传感器组网系统，通过对量测值进行数据融合处理，实现目标无源定位。对于运动目标，通常联合tdoa和fdoa获得目标的位置和速度估计。

目前，现有联合tdoa/fdoa的定位算法主要针对目标辐射源定位系统，基于外辐射源的tdoa/fdoa定位算法研究较少。赵勇胜等针对单站外辐射源提出一种基于极大似然估计的tdoa/fdoa联合定位算法，上述外辐射源tdoa/fdoa定位问题未考虑观测站位置和速度存在误差的情况，而实际问题中观测站常常被安装在卫星、飞机、舰艇或地面车辆等运动平台上，尽管可由定位平台上的导航系统获得接收机位置，但是仍不可避免的含有随机误差。忽视误差的影响会导致目标定位估计性能严重下降，甚至产生虚假目标。因此，外辐射源tdoa/fdoa联合定位和误差校正是外辐射源雷达系统数据处理的一项关键技术。

技术实现要素：

本发明考虑观测站位置和速度误差的影响，针对多发单收外辐射源雷达网tdoa/fdoa定位问题，提出了一种基于两步迭代加权最小二乘估计算法，通过估计运动目标状态(位置和速度)，实现观测站位置误差的校正和运动目标精确定位。

本发明方法的具体步骤是：

步骤1.外辐射源雷达观测站(观测站)接收来自目标散射第三方辐射源发射的信号，得到目标tdoa和fdoa的量测信息；

步骤2.忽略量测噪声和观测站位置以及速度误差的影响，对tdoa量测信息，构造辅助变量rp，将tdoa非线性量测方程转化为伪线性估计方程；

步骤3.忽略量测噪声、观测站位置误差和观测站速度误差的影响，对tdoa伪线性方程对时间求导，构造辅助变量rp和获得fdoa的伪线性估计方程；

步骤4.联立tdoa和fdoa伪线性估计方程，选择运动目标状态为估计向量x，构造线性估计方程z＝hx；

步骤5.考虑量测误差和观测站位置误差对系数矩阵h和z的影响，构造线性方程ε1＝z1-h1x1＝a1n+b1δsr；

步骤6.根据观测站位置和速度误差、tdoa和fdoa量测误差设计权重w1，采用加权最小二乘估计算法得到目标位置xwls＝(h1^tw1h1)^-1h1^tw1z1；

步骤7.在上述估计结果的基础上，考虑待求变量之间的关联性，采用关联最小二乘估计算法对步骤6的估计值进行改进。

本发明的有益效果：

1.考虑观测站位置和速度误差对目标定位性能的影响，根据观测站位置和速度误差以及tdoa和fdoa量测噪声设计优化指标权重，从而降低误差对目标定位性能的影响，提高目标定位精度。

2.通过引入中间变量，合理将多基外辐射源雷达强非线性量测模型转化为伪线性估计方程，在保证估计性能的前提下降低外辐射源定位的复杂度。

3.考虑辅助变量与待求变量之间的关联性，设计关联最小二乘算法，进一步减小估计误差。

具体实施方式：

运动观测站误差下外辐射源雷达时差与频差协同定位方法，该方法包括以下步骤：

步骤1：在多发单收外辐射源雷达网中，包括m个外辐射源和一个观测站。观测站的真实位置位于原点，真实速度为观测站实际位置为sr＝[x0,y0]^t，实际速度为且δsr为观测站的位置误差向量，为观测站的速度噪声向量，并假设为独立的高斯零均值白噪声，其协方差分别为e[δsrδsr^t]＝qs和第m个发射源的坐标向量为p个目标，第p个目标的坐标向量为速度向量为则tdoa和fdoa量测为

式中，和分别为tdoa和fdoa的真实值，c为信号的传播速度c＝3×10⁸m/s，fm为外辐射源m的频率，||·||为欧几里得距离；为目标到观测站位置距离，目标到到观测站位置距离变化率为目标到发射源位置距离，为目标到发射源位置距离变化率；为外辐射源到观测站位置距离，为外辐射源到观测站位置距离变化率；和分别为tdoa和fdoa的量测误差，服从高斯分布。

由于外辐射源的位置和频率先验已知,因此tdoa和fdoa转化为距离和差um,p和距离和差变化率ρm,p

式中，分别为距离和差真实值、距离和差变化率真实值，

为距离和差量测噪声，服从均值为零，方差为qu的高斯分布；为距离和差变化率量测噪声，服从均值为零，方差为qρ的高斯分布。

步骤2：在双基距量测模型中引入中间变量rp，忽略量测噪声和δsr的影响，将上述非线性方程(3)转化为伪线性方程，形式如下

其中，

步骤3：将式(5)等式两边同时对时间求导，得

其中，

步骤4：将目标位置目标速度辅助变量rp和作为待求变量，联立式(5)和(6)，构造线性估计方程

z＝hx(7)

式中，

采用最小二乘估计值获得目标的估计值

步骤5：考虑距离和差量测误差和距离和差变化率以及观测站位置误差δsr和速度误差对h和z的影响，将距离和差量测伪线性方程式(5)和距离和差变化率量测伪线性方程(6)中h和z噪声分量提取出来，构造目标位置伪线性估计方程。将和带入式(5)，展开可得

其中，

将和带入式(6)，展开可得

其中，

联立式(9)与式(10)写成矩阵形式：

ε1＝z1-h1x1＝a1n+b1δs(11)

式中：

b11＝diag(b11(1,1),…,b11(m,p))，

b12＝diag(b12(1,1),…,b12(m,p))，

步骤6：根据距离和差量测误差和距离和差变化率以及观测站位置误差δsr和速度误差设计权重，采用加权最小二乘估计算法得到目标位置的估计值。

步骤6.1：初始化。令迭代次数k＝0，将式(8)获得的最小二乘估计值作为目标初始估计值

步骤6.2：由估计值计算系数矩阵h1，z1，a1和b1。根据观测站位置和速度误差、距离和差以及距离和差变化率量测噪声设计优化指标权重w1，则为量测噪声协方差矩阵，为观测站位置和速度误差的协方差矩阵。

步骤6.3：令k＝k+1，采用加权最小二乘估计获得目标的位置估计值与目标速度估计值与以及中间变量与

步骤6.4：判断其中η1，η2，η3，η4为阈值；若满足条件算法迭代停止，得到目标的位置加权最小二乘估计值否则，转步骤6.2。

步骤7：考虑辅助变量和与目标位置和速度关联性，设计关联最小二乘算法对步骤6的估计值xwls进行改进，具体如下：

步骤7.1：构建关联最小二乘估计模型

ε2＝z2-h2x2＝a2δx1+b2δs(12)

其中，

δx1＝[δx1(1)^t…δx1(p)^t]^t，a2＝blkdiag(a2(1),…,a2(p))

步骤7.2：根据观测站位置误差和速度误差，以及目标状态x1的估计误差协方差设计权重w2＝e[ε2ε2^t]＝(a2cov(x1)a2^t+b2qβb2^t)^-1，cov(x1)＝(h1^tw1h1)^t为目标状态x1的估计误差协方差。

步骤7.3：采用加权最小二乘法估计得到

步骤7.4：x2中变量包含目标位置与观测站位置之差的平方项以及目标速度与观测站速度之差的平方项，要求获得目标的位置需要对x2开根号，目标的位置具体公式如下：

其中，π＝diag{sgn(x1(3p-2)-x0)sgn(x1(3p-1)-y0)}，sgn(·)为符号函数；

目标的速度公式为

获得目标的位置估计值和目标速度的估计值