存在无源探测观测站位置误差的辐射源时差定位方法技术

本发明专利技术公开了一种存在无源探测观测站位置误差的辐射源时差定位方法，其实现过程是：(1)提取时差；(2)构造时差定位直接线性方程组；(3)构造无源探测观测站位置误差扰动的线性映射矩阵；(4)构造时差测量误差扰动的线性映射矩阵；(5)构造误差扰动的线性映射矩阵；(6)估计辐射源位置。本发明专利技术采用了时差定位直接线性模型，考虑了无源探测观测站存在的位置误差，可用于无源探测观测站对辐射源的无源定位。

【技术实现步骤摘要】

本专利技术属于通信
，更进一步涉及到无源定位
中的无源探测观测站存在位置误差时的辐射源时差定位。本专利技术可用于飞机、卫星、无人机等运动平台对辐射源的无源定位。
技术介绍
对辐射源的无源定位被广泛应用于飞机，卫星和无人机等等。现代无源定位系统总是将无源探测观测站放置在位置不确知的运动平台上，从而导致无源探测观测站存在位置误差。微小的无源探测观测站位置误差会导致很大的辐射源位置估计误差，因此考虑无源探测观测站位置误差的辐射源位置估计方法得到了广泛的应用。现有对辐射源时差定位的技术主要采用引入辅助变量的时差定位伪线性方程，需要额外的步骤从辅助变量中提取辐射源位置。陈少昌等人在论文“传感器位置误差条件下的约束总体最小二乘时差定位算法”(航空学报，2013，34(5)：1165-1173)提出了一种考虑无源探测观测站位置误差的时差定位方法。采用利用约束总体最小二乘法，全面考虑无源探测观测站位置误差和时差测量误差估计得到辐射源位置坐标。该方法的具体步骤为：首先通过引入辅助变量将非线性时差定位方程转化为伪线性方程，再利用约束总体最小二乘法，全面考虑伪线性方程的无源探测观测站位置误差和时差测量误差，在此基础上推导了定位方程的目标函数，最后根据牛顿迭代方法进行数值迭代，得到辐射源位置坐标。该方法克服了在无源探测观测站位置误差较大时辐射源位置估计精度性能发散的缺陷。但是，该方法存在的不足是，由于采用了引入辅助变量的伪线性方程组，上述方法需要额外的步骤从额外的辅助变量中提取辐射源位置，限制了辐射源位置的估计效率。西安电子科技大学在其申请的专利“一种基于多运动接收站的时差定位方法”(专利申请号201410171377.3，公开号CN103969622A)中公开了一种基于多运动无源探测观测站的时差定位方法。该方法利用目标位置坐标的初始估计值和加权最小二乘估计得到校正后的目标位置坐标。该方法的具体步骤为：首先，得到目标到各个无源探测观测站的时差定位方程组，转换为引入辅助变量的伪线性方程组，然后利用辅助变量的伪线性方程组和辅助变量估计值得到辅助变量估计值的估计偏差，得到目标位置坐标的初始估计偏差的线性模型的矩阵表达式，最后根据目标位置坐标的初始估计值和加权最小二乘估计得到校正后的目标位置坐标。该方法解决了定位结果具有模糊性以及定位偏差会随着测量误差的增大而迅速增大即定位性能不稳定的问题。但是，该方法存在的不足是，由于采用了引入辅助变量的伪线性方程组，需要额外的步骤从额外的辅助变量中提取辐射源位置，限制了辐射源位置的估计效率。另外，由于未考虑无源探测观测站运动时存在的位置误差，使得无源探测观测站存在位置误差时，辐射源位置估计精度不理想。
技术实现思路
本专利技术的目的在于克服上述已有技术的不足，提出一种存在无源探测观测站位置误差的辐射源时差定位方法。为实现上述目的，本专利技术具体实现步骤包括如下：(1)提取时差：(1a)各无源探测观测站截获辐射源传播的信号，测量辐射源传播的信号到达无源探测观测站各自位置的时间，得到到达时间；(1b)任意一个无源探测观测站，接收其他无源探测观测站转发的辐射源传播的信号到达无源探测观测站各自位置的时间，组成到达时间向量；(1c)构造时差向量的提取矩阵；(1d)将时差向量的提取矩阵与到达时间向量相乘，得到时差向量；(2)构造时差定位直接线性方程组：(2a)将时差向量分成两个时差子向量；(2b)根据两个时差子向量，构造矩阵形式的时差定位直接线性方程组；(3)构造无源探测观测站位置误差扰动的线性映射矩阵：(3a)对时差定位直接线性方程组进行最小二乘求解，得到最小二乘解；(3b)将最小二乘解代入无源探测观测站位置误差扰动的线性映射矩阵的公式，得到无源探测观测站位置误差扰动的线性映射矩阵；(4)构造时差测量误差扰动的线性映射矩阵：(4a)将最小二乘解分别代入两个时差子向量的测量误差扰动的线性映射矩阵的公式，得到两个时差子向量的测量误差扰动的线性映射矩阵；(4b)将两个时差子向量的测量误差扰动的线性映射矩阵相加，得到时差测量误差扰动的线性映射矩阵；(5)构造误差扰动的线性映射矩阵：将无源探测观测站位置误差扰动和时差测量误差扰动的线性映射矩阵代入误差扰动的线性映射矩阵的公式，得到误差扰动的线性映射矩阵；(6)按照下式，估计辐射源位置： u ^ = [ A T ( HRH T ) - 1 A ] - 1 A T ( HRH T ) - 1 b ]]>其中，表示辐射源位置估计，A表示定位方程组系数矩阵，H表示误差扰动的线性映射矩阵，R表示时差测量误差和无源探测观测站位置误差的协方差矩阵，b表示定位方程组系数向量。本专利技术与现有技术相比具有以下优点：第一，本专利技术在构造利用时差求解辐射源位置的方程组的过程中，采用了时差定位直接线性模型，得到了时差定位直接线性方程组，克服了现有技术需要较多的步骤和较长的计算时间才能得到辐射源位置的缺点，使得本专利技术能够提高辐射源位置的估计效率。第二，本专利技术在对定位直接线性方程组进行求解的过程中考虑了无源探测观测站位置误差，克服了现有技术在无源探测观测站位置有误差情形下的性能局限，使得本专利技术能够在无源探测观测站存在位置误差情况下提高辐射源位置的估计精度。附图说明图1是本专利技术的流程图；图2是本专利技术与现有技术对辐射源位置估计的均方根误差随无源探测观测站位置误差变化曲线图；图3是本专利技术与现有技术对辐射源位置估计的平均耗时随无源探测观测站数目变化曲线图。具体实施方式下面结合附图对本专利技术做进一步的描述。参照附图1，本专利技术的具体步骤如下。步骤1，提取时差。各无源探测观测站截获辐射源传播的信号，测量辐射源传播的信号到达无源探测观测站各自位置的时间，得到到达时间。任意一个无源探测观测站，接收其他无源探测观测站转发的辐射源传播的信号到达无源探测观测站各自位置的时间，组成到达时间向量。构造时差向量的提取矩阵。时差向量的提取矩阵是按照如下公式构造的：F＝[-EM-1,IM-1]其中，F表示时差向量的提取矩阵，ΕM-1表示(M-1)×1维元素均为1的列向量，IM-1表示M-1维单位矩阵，M表示无源探测观测站数目，M-1表示M个无源探测观测站得到的时差的数目。将时差向量的提取矩阵与到达时间向量相乘，得到时差向量。步骤2，构造时差定位直本文档来自技高网...

【技术保护点】
一种存在无源探测观测站位置误差的辐射源时差定位方法，包括如下步骤：(1)提取时差：(1a)各无源探测观测站截获辐射源传播的信号，测量辐射源传播的信号到达无源探测观测站各自位置的时间，得到到达时间；(1b)任意一个无源探测观测站，接收其他无源探测观测站转发的辐射源传播的信号到达无源探测观测站各自位置的时间，组成到达时间向量；(1c)构造时差向量的提取矩阵；(1d)将时差向量的提取矩阵与到达时间向量相乘，得到时差向量；(2)构造时差定位直接线性方程组：(2a)将时差向量分成两个时差子向量；(2b)根据两个时差子向量，构造矩阵形式的时差定位直接线性方程组；(3)构造无源探测观测站位置误差扰动的线性映射矩阵：(3a)对时差定位直接线性方程组进行最小二乘求解，得到最小二乘解；(3b)将最小二乘解代入无源探测观测站位置误差扰动的线性映射矩阵的公式，得到无源探测观测站位置误差扰动的线性映射矩阵；(4)构造时差测量误差扰动的线性映射矩阵：(4a)将最小二乘解分别代入两个时差子向量的测量误差扰动的线性映射矩阵的公式，得到两个时差子向量的测量误差扰动的线性映射矩阵；(4b)将两个时差子向量的测量误差扰动的线性映射矩阵相加，得到时差测量误差扰动的线性映射矩阵；(5)构造误差扰动的线性映射矩阵：将无源探测观测站位置误差扰动和时差测量误差扰动的线性映射矩阵代入误差扰动的线性映射矩阵的公式，得到误差扰动的线性映射矩阵；(6)按照下式，估计辐射源位置：u^=[AT(HRHT)-1A]-1AT(HRHT)-1b]]>其中，表示辐射源位置估计，A表示定位方程组系数矩阵，H表示误差扰动的线性映射矩阵，R表示时差测量误差和无源探测观测站位置误差的协方差矩阵，b表示定位方程组系数向量。...

【技术特征摘要】
1.一种存在无源探测观测站位置误差的辐射源时差定位方法，包括如下步骤：(1)提取时差：(1a)各无源探测观测站截获辐射源传播的信号，测量辐射源传播的信号到达无源探测观测站各自位置的时间，得到到达时间；(1b)任意一个无源探测观测站，接收其他无源探测观测站转发的辐射源传播的信号到达无源探测观测站各自位置的时间，组成到达时间向量；(1c)构造时差向量的提取矩阵；(1d)将时差向量的提取矩阵与到达时间向量相乘，得到时差向量；(2)构造时差定位直接线性方程组：(2a)将时差向量分成两个时差子向量；(2b)根据两个时差子向量，构造矩阵形式的时差定位直接线性方程组；(3)构造无源探测观测站位置误差扰动的线性映射矩阵：(3a)对时差定位直接线性方程组进行最小二乘求解，得到最小二乘解；(3b)将最小二乘解代入无源探测观测站位置误差扰动的线性映射矩阵的公式，得到无源探测观测站位置误差扰动的线性映射矩阵；(4)构造时差测量误差扰动的线性映射矩阵：(4a)将最小二乘解分别代入两个时差子向量的测量误差扰动的线性映射矩阵的公式，得到两个时差子向量的测量误差扰动的线性映射矩阵；(4b)将两个时差子向量的测量误差扰动的线性映射矩阵相加，得到时差测量误差扰动的线性映射矩阵；(5)构造误差扰动的线性映射矩阵：将无源探测观测站位置误差扰动和时差测量误差扰动的线性映射矩阵代入误差扰动的线性映射矩阵的公式，得到误差扰动的线性映射矩阵；(6)按照下式，估计辐射源位置： u ^ = [ A T ( HRH T ) - 1 A ] - 1 A T ( HRH T ) - 1 b ]]>其中，表示辐射源位置估计，A表示定位方程组系数矩阵，H表示误差扰动的线性映射矩阵，R表示时差测量误差和无源探测观测站位置误差的协方差矩阵，b表示定位方程组系数向量。2.根据权利要求1所述的存在无源探测观测站位置误差的辐射源时差定位方法，其特征在于：步骤(1c)中所述的时差向量的提取矩阵是按照如下公式构造的：F＝[-EM-1,IM-1]其中，F表示时差向量的提取矩阵，ΕM-1表示(M-1)×1维元素均为1的列向量，IM-1表示M-1维单位矩阵，M表示无源探测观测站数目，M-1表示M个无源探测观测站得到的时差的数目。3.根据权利要求1所述的存在无源探测观测站位置误差的辐射源时差定位方法，其特征在于：步骤(2a)中所述的将时差向量分成两个时差子向量的具体步骤如下：第一步，按照下式，构造高投影矩阵和低投影矩阵：Th＝[0M-2,IM-2]T1＝[IM-2,0M-2]其中，Th表示高投影矩阵，0M-2表示(M-2)×1维元素都为0的列向量，IM-2表示M-2维单位矩阵，M表示无源探测观测站数目，M-2表示将M个无源探测观测站得到的时差向量分成两个时差子向量中时差的数目，T1表示低投影矩阵；第二步，按照下式，将时差向量分成两个时差子向量：τH＝ThττL＝T1τ其中，τH和τL分别表示两个时差子向量，Th和T1分别表示高投影矩阵和低投影矩阵，τ表示时差向量。4.根据权利要求1所述的存在无源探测观测站位置误差的辐射源时差定位方法，其...

【专利技术属性】
技术研发人员：董春曦，朱颖童，董阳阳，赵国庆，饶鲜，许锦，唐峥钊，刘明明，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61