1.基于误差校正的联合DOA与TOA单站无源定位方法，其特征在于：该方法包括以下步骤：步骤1：外辐射源雷达系统接收站接收来自目标散射第三方辐射源发射的信号，得到目标到达方位角DOA和到达时间TOA的量测信息，具体如下其中，接收站位于原点(0,0)，发射源个数为M，第m个发射源的坐标向量为qm＝(xm,ym)；

P个目标，第p个目标的坐标向量为 θp为目标p到接收站的到达方位角量测值，Δθp为到达方位角系统误差， 为到达方位角量测误差，服从均值为0，方差为 的高斯分布；τm,p为发射源m和接收源构成的外辐射源雷达系统对目标p测量的到达时间量测值，Δτm为该到达时间的系统误差， 为到达时间量测误差，服从均值为0，方差为 的高斯分布；c为光速，c＝3×108m/s；

为简化将式(2)两边乘以c，可以得到其中，um,p＝cτm,p为发射源m和接收源构成的外辐射源雷达系统对目标p的双基距量测值；δm＝cΔτm为对应的双基距量测系统误差， 为双基距量测误差，服从均值为0，方差为 的高斯分布；

步骤2：引入辅助变量 忽略量测噪声影响，将非线性方程式(1)和式(3)转化为伪线性方程；具体形式如下其中,

建立目标位置状态与接收站位置状态之间的伪线性模型，将式(4)和(5)写成矩阵形式，具体如下其中，

步骤3：采用迭代最小二乘法获得目标位置和系统误差的估计值；

步骤3.1：初始化，令迭代次数k＝0，系统误差初值 和目标距接收站的距离初值

步骤3.2：令k＝k+1，采用最小二乘法求得 得到目标位置估计值以及角度误差估计值 p＝1,2,...,P和双基距误差估计值

步骤3.3：判断 和 若不满

足，则转步骤3.2；否则，迭代停止，得到目标位置和系统误差的初始估计值 ε1和ε2为阈值；

步骤4：考虑量测误差对 和 的影响，将伪线性方程即式(6)中 和 噪声分量提取出来，构造基于约束总体最小二乘估计的形式；

步骤4.1：同时考虑量测误差对 和 的影响，对式(6)进行一阶泰勒展开，并忽略二阶及以上误差项，得到式中，

Ea为E第a项分量，L为X向量的长度；

步骤4.2：构建误差配准模型为一个约束总体最小二乘估计模型，具体如下：步骤4.3：将有约束的优化问题式(9)转化为如下无约束优化问题形式其中，f(X)为目标函数， C＝[H Z]，xc为X向量重第c个分量；

步骤5：利用迭代牛顿算法对式(10)求解，得到其中，n为迭代次数，X(n)为第n次迭代的待估参数；Hn为f(X)的海塞矩阵，Tn为f(X)的梯度向量，具体求解如下：

步骤5.1：将步骤3得到的目标位置和系统误差估计值作为初始值 n＝0，n表示迭代次数；

步骤5.2：令n＝n+1，采用迭代牛顿算法求解公式(11)，得到新的目标位置和系统误差估计值X(n)，即目标位置估计值 距离估计值 以及系统误差估计值 和

步骤5.3：当 和 迭代停止，

得到目标位置和系统误差的估计值 否则，转步骤5.2，ε3和ε4为阈值；

步骤6：考虑辅助变量 和目标位置 的相关性，采用相关最小二乘算法对步骤5的估计值 进行改进，具体如下：步骤6.1：选择目标位置平方项 角度误差和双基距误差估计值 作为变量，构建系统误差配准模型，将基于CTLS的估计值 带入系数矩阵B'和Z'，得到Z'＝H'X'+B'ΔX                            (12)其中，

步骤6.2：对式(12)采用加权最小二乘，得到估计值如下X'LS＝(H'TW'-1H')-1H'TW'-1Z'                   (13)其中，

步骤6.3：根据步骤6.2获得的估计值X'LS，得到目标位置的平方项对其开方，在二维空间得到目标位置 的4个候选值，从中选择最佳目标位置估计值 具体如下

其中，最佳目标位置估计值 Zp为接收站接收M个外辐射源对目标p的量测T T值，Zp＝[θp u1,p...uM,p]；系统误差Δ＝[Δθp δ1...δM]；

dp表示接收站获得目标p的方位角和双基距的估计值；

最终，得到目标位置的估计 从而实现系统误差下联合DOA与TOA信息的单站无源定位。