一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法技术

本发明专利技术公开了一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法。步骤1：通过GNSS宽带干扰定位系统中获取的测量值特征，分别构建到达角AOA观测模型和差分到达时间TD0A观测模型；步骤2：基于步骤1中的到达角AOA观测模型，建立到达角AOA定位模型；步骤3：基于步骤1中的差分到达时间TDOA观测模型，建立差分到达时间TDOA定位模型；步骤4：基于步骤2和步骤3的到达角AOA定位模型和差分到达时间TDOA定位模型，建立面向GNSS宽带干扰的到达角AOA与差分到达时间TDOA松组合定位模型，实现最终定位。针对现有AOA/TDOA联合定位方法的不足。

【技术实现步骤摘要】
一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法
本专利技术属于卫星导航
；具体涉及一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法。
技术介绍
由于个人隐私设备的急剧增加，对于安全至关重要的GNSS用户如机场和海港等，其当务之急是要随时了解本地GNSS干扰，其中宽带干扰发射器是较为有效且常见的干扰手段。宽带干扰与频率无关，并且会导致附近的GNSS用户持续丢失GNSS信号。因此，当宽带GNSS干扰器对GNSS用户构成更大危险时，必须能够识别宽带GNSS干扰器并对其进行地理位置定位。针对GNSS宽带干扰定位需求，主要采用两种技术对其进行地理定位：AOA方法和TDOA方法，而这两种定位方法都各自有其优势和不足。因此，相比于单一定位方法的不足，将所有站点的AOA和TDOA信息结合起来进行宽带干扰定位具有较高的实际应用价值。目前，针对到达角(AngleofArrival，AOA)方法和差分到达时间(TimeDifferenceofArrival，TDOA)方法的研究已经取得了较多成果，而AOA/TDOA联合定位方法研究还不成熟。一些学者采用了非加权算法进行联合定位，但常规算法没有考虑异方差下AOA和TDOA测量值的融合，并且还没有公平地表征它们相对于现有方法的优点。
技术实现思路
本专利技术提出一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，针对现有AOA/TDOA联合定位方法的不足。本专利技术通过以下技术方案实现：一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，所述增强型GNSS宽带干扰定位方法包括以下步骤：步骤1：通过GNSS宽带干扰定位系统中获取的测量值特征，分别构建到达角AOA观测模型和差分到达时间TD0A观测模型；步骤2：基于步骤1中的到达角AOA观测模型，建立到达角AOA定位模型；步骤3：基于步骤1中的差分到达时间TDOA观测模型，建立差分到达时间TDOA定位模型；步骤4：基于步骤2和步骤3的到达角AOA定位模型和差分到达时间TDOA定位模型，建立面向GNSS宽带干扰的到达角AOA与差分到达时间TDOA松组合定位模型，实现最终定位。进一步的，所述步骤2中在各基站到达角AOA估计具有不等方差情况下，求解仅到达角AOA定位方法的干扰机坐标。进一步的，所述步骤3中在各基站差分到达时间TDOA估计具有不等方差情况下，求解仅差分到达时间TDOA定位方法的干扰机坐标。进一步的，所述步骤1具体包括以下步骤：步骤1.1：利用多个基站天线系统获取到达角AOA的测量值，步骤1.2：选取的到达角AOA观测模型为正态分布，步骤1.3：根据步骤1.1与步骤1.2定义到达角AOA观测模型包含高斯噪声，以此建立到达角AOA的观测模型，其中，θl0表示真实的AOA，表示AOA的测量值，表示正态分布的方差，l表示基站序号，θ表示中间变量，N表示正态分布；步骤1.4：获得各基站观测到的差分到达时间TDOA，步骤1.5：选取的差分到达时间TDOA观测模型为多元高斯分布，步骤1.6：根据步骤1.4与步骤1.5定义多元高斯分布的均值τ0和协方差分量Στ，以此建立差分到达时间TDOA的观测模型，式中：其中，表示是从l站到第1站的真实TDOA，l∈{1,2...P}，表示第P站的TDOA测量方差，P表示用于TDOA测量基站个数。进一步的，所述步骤2所包含以下步骤：步骤2.1：根据基站观测到的到达角AOA测量值，建立基站坐标与干扰坐标间的数学模型具体为，其中，θl是在站点l处观测到的宽带干扰源AOA，l∈{1,2...L}，L代表用于AOA测量的基站个数，(xu,yu)是干扰源的二维东北笛卡尔坐标，而(xl,yl)是第l个站点的笛卡尔坐标；当θl＝±90°时，tan(θl)容易出现数值不稳定，从而导致逼近无穷大，所以将两边都乘以cos(θl)：因此，干扰机坐标(xu,yu)的普通最小二乘估计为：然而，干扰信号传播到每个不同的接收站时受到不同路径损耗和局部多径以及天线阵列本身的作用，导致AOA估计具有不相等的方差，上述AOA定位方法将不再是最佳线性无偏估计器，需进行修正；步骤2.2：根据步骤2.1的数学模型构建各基站到达角AOA估计具有不等方差情况下的到达角AOA测量误差协方差矩阵Σθ∈RN×N，其中Σθ非对角元素为零，对角元素[Σθ]l,l对应于第l个基站的AOA误差方差；步骤2.3：根据步骤2.2的测量误差协方差矩阵求解干扰坐标对到达角AOA测量值的偏导数，其中，rl表示第l个基站到干扰源的距离；确定到达角AOA测量向量的雅可比矩阵JA，步骤2.4：根据步骤2.3的雅可比矩阵推导高斯-牛顿迭代法，求解到达角AOA定位模型下干扰机坐标其中，[Δxu,Δyu]T表示由角度测量误差Δθl引起的干扰源位置估计误差，Σθ表示到达角AOA测量误差协方差矩阵，表示第l个基站到达角AOA的估计值。进一步的，所述步骤3所包含以下步骤：步骤3.1：对于宽带GNSS干扰器，使用互相关信号处理方法测量在几个站之间到达的差分到达时间TDOA，建立差分到达时间TDOA测量与宽带干扰坐标(xu,yu)的数学模型，其中，||||是向量的欧几里得距离，c为电磁波传播速度，τij表示第i个基站与第j个基站间的差分到达时间TDOA，ri表示第i个基站到宽带干扰机的距离，rj表示第j个基站到宽带干扰机的距离，P表示TDOA测量基站的个数；步骤3.2：根据步骤3.1的数学模型构建各差分到达时间TDOA测量具有不等方差情况下的差分到达时间TDOA误差协方差矩阵Στ，其中Στ非对角元素为零，对角元素[Στ]l,l对应于第l个基站到第1个基站的TDOA误差方差；步骤3.3：根据步骤3.2的测量误差协方差矩阵求解干扰坐标对差分到达时间TDOA测量值的偏导数，确定差分到达时间TDOA测量向量的雅可比矩阵JT，其中，τN,1表示第P个站与第1个基站间的差分到达时间TDOA；步骤3.4：根据步骤3.3的雅可比矩阵推导高斯-牛顿迭代法，求解差分到达时间TDOA定位模型下干扰机坐标，其中，表示第i站到第1个基站的差分到达时间TDOA的估计值，Δτi1表示第i个站到第1个基站的差分到达时间TDOA的测量误差。进一步的，所述步骤4所包含以下步骤：步骤4.1：分别获取到达角AOA定位估计坐标[xA，yA]T与差分到达时间TDOA定位估计坐标[xT，yT]T及其位置误差协方差ΣA和ΣA；步骤4.2：将到达角AOA和差分到达时间TDOA获得的位置信息进行松组合，步骤4.3：将步骤4.2得到的松组合通过加权最小二乘求解到达角AOA和差分到达本文档来自技高网...

【技术保护点】
1.一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，其特征在于，所述增强型GNSS宽带干扰定位方法包括以下步骤：/n步骤1：通过GNSS宽带干扰定位系统中获取的测量值特征，分别构建到达角AOA观测模型和差分到达时间TD0A观测模型；/n步骤2：基于步骤1中的到达角AOA观测模型，建立到达角AOA定位模型；/n步骤3：基于步骤1中的差分到达时间TDOA观测模型，建立差分到达时间TDOA定位模型；/n步骤4：基于步骤2和步骤3的到达角AOA定位模型和差分到达时间TDOA定位模型，建立面向GNSS宽带干扰的到达角AOA与差分到达时间TDOA松组合定位模型，实现最终定位。/n

【技术特征摘要】
1.一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，其特征在于，所述增强型GNSS宽带干扰定位方法包括以下步骤：
步骤1：通过GNSS宽带干扰定位系统中获取的测量值特征，分别构建到达角AOA观测模型和差分到达时间TD0A观测模型；
步骤2：基于步骤1中的到达角AOA观测模型，建立到达角AOA定位模型；
步骤3：基于步骤1中的差分到达时间TDOA观测模型，建立差分到达时间TDOA定位模型；
步骤4：基于步骤2和步骤3的到达角AOA定位模型和差分到达时间TDOA定位模型，建立面向GNSS宽带干扰的到达角AOA与差分到达时间TDOA松组合定位模型，实现最终定位。


2.根据权利要求1所述的一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，其特征在于，所述步骤2中在各基站到达角AOA估计具有不等方差情况下，求解仅到达角AOA定位方法的干扰机坐标。


3.根据权利要求1所述的一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，其特征在于，所述步骤3中在各基站差分到达时间TDOA估计具有不等方差情况下，求解仅差分到达时间TDOA定位方法的干扰机坐标。


4.根据权利要求1所述的一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，其特征在于，所述步骤1具体包括以下步骤：
步骤1.1：利用多个基站天线系统获取到达角AOA的测量值，
步骤1.2：选取的到达角AOA观测模型为正态分布，
步骤1.3：根据步骤1.1与步骤1.2定义到达角AOA观测模型包含高斯噪声，以此建立到达角AOA的观测模型，



其中，θl0表示第l个基站真实的AOA，表示第l个基站AOA的测量值，表示第l个基站的AOA测量方差，θ表示中间变量，N表示正态分布；
步骤1.4：获得各基站观测到的差分到达时间TDOA，
步骤1.5：选取的差分到达时间TDOA观测模型为多元高斯分布，
步骤1.6：根据步骤1.4与步骤1.5定义多元高斯分布的均值τ0和协方差分量Στ，以此建立差分到达时间TDOA的观测模型，



式中：



其中，表示是从l站到第1站的真实TDOA，l∈{1,2...P}，表示第P站的TDOA测量方差，P表示用于TDOA测量基站个数。


5.根据权利要求2所述的一种基于AOA/TDOA联合的增强型GNSS宽带干扰定位方法，其特征在于，所述步骤2所包含以下步骤：
步骤2.1：根据基站观测到的到达角AOA测量值，建立基站坐标与干扰坐标间的数学模型具体为，



其中，θl是在站点l处观测到的宽带干扰源AOA，l∈{1,2...L}，L代表用于AOA测量的基站个数，(xu,yu)是干扰源的二维东北笛卡尔坐标，而(xl,yl)是第l个站点的笛卡尔坐标；当θl＝±90°时，tan(θl)容易出现数值不稳定，从而导致逼近无穷大，所以将两边都乘以cos(θl)：



因此，干扰机坐标(xu,yu)的普通最小二乘估计为：



然而，干扰信号传播到每个不同的接收站时受到不同路径损耗和局部多径以及天线阵列本身的作用，导致AOA估计具有不相等的方差，上述AOA定位方法将不再是最佳线性无偏估计器，需进行修正；
步骤2.2：根据步...

【专利技术属性】
技术研发人员：沈锋，赵雨晴，徐定杰，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23