一种基于线性关系的两步估计高精确时差定位算法
【专利摘要】基于我们申请并已公布的专利201110197395.5（一种基于线性关系的高精确时差定位算法），本发明进一步公开了一种基于线性关系的两步估计高精确时差定位算法。首先利用已公布专利201110197395.5中的方法完成辐射源位置的初步估计。在关键的第二步估计中，基于第一步估计得到的辐射源位置得到另一线性方程组，该方程组以辐射源位置修正量为变量。基于最小均方原则，第二步估计得到辐射源位置修正量。综合两步估计结果，得到辐射源位置最终估计值。由于两步估计中皆用线性关系，所得解唯一。通过第二步估计技术对已公布专利201110197395.5所得结果进行修正，本发明方法大幅度提高了估计精度和鲁棒性。
【专利说明】一种基于线性关系的两步估计高精确时差定位算法
【技术领域】
[0001] 本发明涉及时差定位算法，尤其是基于线性关系的高精确高鲁棒性时差定位算 法。
【背景技术】
[0002] 时差定位属于双曲定位范畴，又称"反罗兰定位法"。用三个以上已知相对距离 的感知器，测量同一辐射源发出的信号，根据信号到达各接收点的时间差，确定发射源的位 置。时差定位技术在机场进近着陆、移动车辆及其它辐射目标定位、嵌入声纳或雷达阵列构 成定位系统等方面都具有重要应用。此外，在军事上，可以对对目标进行高精度、实时、隐蔽 定位，发展该技术可有效提高卫星对地目标侦查能力，掌握信息化战场的主动性。
[0003] 时差定位的基础是精确获得接收辐射源信号的不同传感器间时差信号，可以通过 最大似然估计、联合估计及相关估计等方法求得，也可以通过一些信号模型实现。获得时差 信号后，时差定位算法对估计辐射源位置具有关键作用。目前所用时差定位算法主要有泰 勒级数法、分而治之法、球面交叉法、球面插值法，Fang方法，Chan方法以及约束条件下的 全局最小方差定位方法。最近有结合时差和频差进行定位的报道，但时差定位本身尚需进 一步研究。
[0004] 泰勒级数法首先对双曲定位方程同时线性化，用泰勒级数估计这些线性代数方程 的最小均方误差解。泰勒级数法最大的问题是其收敛性没有保障，所求解是否是辐射源位 置与初始估计密切相关。相关估计方法的实质在于寻找相关函数的"尖峰"。而当SNR达不 到要求时，特别是处理窄带信号时，相关峰"宽胖"不明显，在其周围还有若干干扰峰存在， 导致所谓的"模糊"现象，影响估计的准确性及其方差性能。当时差值分组合理时，分而治 之法是无偏的。否则该方法产生的均方误差超过Cramer-Rao下界。因此，如何恰当地给时 差值分组是十分关键的。球面交叉法利用交叉球面确定辐射源位置，其利用的代数关系为 辐射源距离的二次式。基于建立在辐射源距离和辐射源深度基础上的误差方程组，球面插 值法把距离差转化为位置估计。在最大似然意义下，Fang方法是最优的。但是这种方法只 适合二维定位。通过把三维定位问题分解为多个二维定位问题，Fang方法也可用于三维定 位。Y. T. Chan和K. C. Ho在1994年提出的Chan方法，对于固定辐射源给出了近似最大 似然估计的闭式解，已成为公认的解决双曲线定位的最好方法之一，并得到了广泛的应用。 尽管给出的方法在理论上可以趋于精确值，但是在运算过程中需要对大地坐标系与空间坐 标系的进行转换，并且迭代的数值需要近似，这些都会对结果造成一定的误差和模糊。
[0005] 除上述问题以外，所有以上方法都还具有一共同缺陷：即使时差值足够，辐射源位 置也不能唯一求出。球面交叉法和Fang方法求解一个二次式，一般都产生两个解。简单认 为其中一解在合理范围，而摒弃另外一解。这在一些极端情况下并不适用。球面插值法运 行一次得到辐射源距离和辐射源深度信息，但其必须运行两次，所以也会产生两个解。Chan 方法利用的几何关系是辐射源位置向量的二次式，即使在二维空间中使用四个甚至更多的 感知器，仍然会产生两个解，甚至会产生复数解。所以，在实际应用中，这些已有方法都必须 结合一定先验知识，判定各所求解的合理性。
[0006] 2012年发明人刘怡光等在已公布专利201110197395. 5 (-种基于线性关系的高 精确时差定位算法）中提出了一种新型求解辐射源位置的线性方法
【权利要求】
1. 基于线性关系的两步估计高精确时差定位算法，该方法首先利用已公布专利 201110197395. 5中的方法完成辐射源位置的初步估计；在关键的第二步估计中，基于第一 步估计得到的辐射源位置得到另一线性方程组，该方程组以辐射源位置修正量为变量；基 于最小均方原则，第二步估计得到辐射源位置修正量；综合两步估计结果，得到辐射源位置 最终估计值；由于两步估计中皆用线性关系，所得解唯一；通过第二步估计技术对已公布 专利201110197395. 5所得结果进行修正，本发明方法大幅度提高了估计精度和鲁棒性。
2. 根据权利要求1所述的基于线性关系的两步估计高精确时差定位算法，其特征在于 它的步骤如下： 步骤一：利用已公布专利201110197395. 5中的方法完成辐射源位置的初步估计； 步骤二：基于第一步估计得到的辐射源位置得到另一线性方程组，该方程组以辐射源 位置修正量为变量； 步骤三：基于最小均方原则，第二步估计得到辐射源位置修正量； 步骤四：综合两步估计结果，得到辐射源位置最终估计值。
3. 根据权利要求2所述的基于线性关系的两步估计高精确时差定位算法，其特征在于 步骤一：基于我们申请并已公布的专利201110197395. 5,构建线性方程组，以辐射源位置 向量为变量，各接收传感器的位置向量和时差为系数或偏移向量；在定位传感器足够时，用 最小均方误差法从线性定位方程组初步估计辐射源位置 r = (ArAy1ArB (1) 矩阵J中第i行定义为,矩阵f中第i行定义为 之=+ IWw-! + w
4. 根据权利要求2所述的基于线性关系的两步估计高精确时差定位算法，其特征在于 步骤二：基于（1)式的估计所得r，可得到信号从r传播到弋的时间||r -乃I/r与实际量 测所得时间?,之差为 Sti = - ||r - Ti I / c (2) 假设r修正为^ 时满足 ||r + fir - JriI = Cti (3) 则由⑵和（3)得到 ||r + fir - r,.|| - |r - ^.|] = CSti (4) 进一步本发明得到以辐射源位置修正量为变量的另一组线性方程组为 Κδτ = L (5) 其中矩阵JT的第i行定义为2(^.+1 - ，矩阵i；的第i行定义为 2ψ?】.|η. - r| - 5ii+1|ri+1 - _r||]+ (CvJij.)2 - (c5ii+J。
5. 根据权利要求2所述的基于线性关系的两步估计高精确时差定位算法，其特征在 于步骤三：基于最小均方原则，第二步估计得到辐射源位置修正量 &' = If1Kj1KjL (6) 〇
6.根据权利要求2所述的基于线性关系的两步估计高精确时差定位算法，其特征在于 步骤四：综合两步估计结果，得到辐射源位置最终估计值 Trjt = r + Sr (7) 以上基于最小均方原则求解线性关系的两步估计高精确时差定位算法，仅为本发明的 基于线性关系的两步估计高精确时差定位算法的较佳实例之一，并非用于限定本发明的实 质技术内容范围，实质内容已广泛地定义于权利要求书中，任何他人所完成的技术实体和 方法，若是与权利要求书中所定义者完全相同，或同一等效变更(如采用其它等效方法推导 出的线性关系)，均将被视为涵盖于此专利保护范围之内。
【文档编号】G01S5/00GK104375118SQ201410636004
【公开日】2015年2月25日 申请日期:2014年11月13日 优先权日:2014年11月13日
【发明者】刘怡光, 曹丽萍 申请人:四川大学