[发明专利]一种基于ADS-B的三维椭球交叉GNSS干扰源定位方法

权利要求书

1.一种基于ADS-B的三维椭球交叉GNSS干扰源定位方法，其特征在于，包括如下步骤：

步骤1：获取ADS-B接收机的所有飞机的ADS-B数据信息；

步骤2：搜索ADS-B数据信息中有位置信息丢失的航线；

步骤3：分析每条位置信息丢失航线的位置丢失点和位置重获取点；

步骤4，利用线性化降维及最小二乘法实现三维椭球GNSS干扰源定位：

对于第k条有位置丢失点和重获取点的航线，建立以下方程：

其中，x_J、y_J、z_J为干扰源的三维坐标，x_Lk、y_Lk、z_Lk为第k条有位置丢失点和重获取点航线的位置丢失点的三维坐标，x_Rk、y_Rk、z_Rk为第k条有位置丢失点和重获取点航线的位置重获取点的三维坐标；

步骤41：对方程(1)进行去根号，可得：

步骤42：对方程(2)进行线性化降维，记：

通过寻找GNSS干扰源的初始位置(x_J,0,y_J,0,z_J,0)开始线性化，更新的干扰源的坐标为：

其中，Δx_J1、Δy_J1、Δz_J1为增量，f_k(x_J,1,y_J,1,z_J,1)的泰勒展开式为：

经过这种变化后，就实现了方程(2)的线性化方程：

步骤43：对于有位置丢失点和重获取点的航线，都建立方程(6)，当多于三条及三条以上这样的航线时，可建立求解干扰源位置的方程组：

对于该方程组，采用最小二乘进行求解，最小二乘是给出一个无解的方程Ax＝b，A有m行和n列，mn；其中观测值b₁,...,b_m的个数大于求解参数的个数，称之为最优解应使误差向量的长度最小，则等价于使得||e||²＝(b-Ax)^T(b-Ax)最小，最小化该二次方程得到：

采用最小二乘可求得方程组(7)的解Δx_J1、Δy_J1、Δz_J1；

步骤44：求得Δx_J1、Δy_J1、Δz_J1后，由方程(4)求得干扰源的新位置，将该新位置作为新的初始位置(x_J,0,y_J,0,z_J,0)，从方程(5)到方程(7)循环计算，直到增量Δx_J1、Δy_J1、Δz_J1收敛到设定阈值为止。