CN107390166A - 一种自适应干扰源定位飞行校验方法 - Google Patents

Abstract

本发明公开了一种自适应干扰源定位飞行校验方法，可以实现对地面各类无线电干扰源的监测和定位，保证民航无线电安全。飞行校验无线电干扰定位的基本原理是借助于机载无线电侦测设备获取干扰源来波方位，结合校验飞机自身的位置、速度、飞行姿态、高度等飞行参数，应用单站交叉定位原理，采用一种自适应滤波方法，能够准确的求得干扰源的位置信息。本发明的主要步骤包括：建立校验飞机干扰侦测系统的状态方程；建立校验飞机干扰侦测系统的非线性量测方程；离散化系统状态方程和量测方程；利用自适应干扰源定位算法求解离散后的方程，输出干扰源的位置。本发明所采用的自适应干扰源定位算法，提高了定位算法的定位精度和定位过程的稳定性。

Description

一种自适应干扰源定位飞行校验方法

技术领域

本发明涉及一种自适应干扰源定位飞行校验方法，应用于航空飞行校验技术领域和无线电侦测技术领域。

背景技术

随着我国无线电台站总量的不断增加,以及民用航空机场、航线、航班数量的日益增多，民用航空无线电专用频率遭受电磁干扰的事件逐年呈上升趋势。无线电干扰对民航通信系统和空中交通管制系统构成了严重威胁，轻则使飞机改变原来的航向或使航班延误而带来巨大经济损失，重则威胁飞机及广大乘客的安全。对于民用航空无线电专用频率遭受电磁干扰的问题，需要加强对影响民航无线电专用频率安全使用相关环节的管理，提高无线电监测与干扰侦测的能力，保证必要的装备和手段。如何根据当前的形式，大力加强民航无线电干扰侦测系统设备建设，及时有效地查处干扰源，是民航飞行校验面临的重要课题。

飞行校验干扰定位问题是一种典型的单站交叉定位问题，其本质上是一类参数估计问题，属于非线性滤波问题。传统的滤波方法有最小二乘法、扩展卡尔曼滤波(EKF)和无迹卡尔曼滤波(UKF)。EKF在非线性系统中应用最广泛，通过泰勒级数展开将非线性滤波问题转化为近似线性化问题，但是由于线性化忽略了泰勒级数最高项，引入了截断误差；同时过程噪声协方差和量测噪声协方差以及初始状态均采用预先估计的方式确定，如果估计得不太准确，会导致滤波发散。

发明内容

本发明的目的是为了解决上述问题，提出一种自适应干扰源定位飞行校验方法，实现对地面干扰源的精确测定。

本发明的一种飞行程序校验干扰源定位的方法，应用于飞行校验，可以实现对地面各类无线电干扰源的监测和定位，保证民航无线电安全。

飞行校验无线电干扰定位的基本原理是借助于机载无线电侦测设备获取干扰源来波方位，结合校验飞机自身的位置、速度、飞行姿态、高度等飞行参数，应用单站交叉定位原理，采用一种自适应滤波方法，能够准确求得干扰源的位置信息。本发明的主要步骤包括：建立校验飞机干扰侦测系统的状态方程；建立校验飞机干扰侦测系统的非线性量测方程；离散化系统状态方程和量测方程；利用自适应干扰源定位算法求解离散后的方程，输出干扰源的位置。本发明所采用的自适应干扰源定位算法，提高了定位算法的定位精度和定位过程的稳定性。

本发明的优点在于：

(1)能够实时进行干扰源的测向和定位；

(2)通过对初值自适应选择和量测误差的自适应修正，提高了卡尔曼滤波的稳定性，从而提高定位的稳定性。

(3)通过卡尔曼滤波过程中的迭代处理，减小扩展卡尔曼滤波的截断误差，从而提高定位的精确度。

附图说明

图1为现有民航干扰侦测的流程图；

图2为本发明的总体流程图；

图3为自适应定位算法流程图。

具体实施方式

下面将结合附图和实施例对本发明作进一步的详细说明。

本发明的一种自适应干扰源定位飞行校验方法，包括两部分：一是通过总体递推最小二乘法获得干扰源粗略位置；二是将获得的干扰源粗略位置作为新息迭代扩展卡尔曼滤波的初始值进行求解，得到最终的干扰源位置，如图3所示，具体说明如下：

步骤一：通过总体递推最小二乘法获得干扰源位置

在飞行校验过程中，利用飞行校验系统中的数据分析设备可以得到飞机飞行位置的一系列位置点即航迹，同时又通过机载无线电侦测设备可以获得干扰源的来波方位角，计算各个航迹点，通过交叉定位可以求得干扰源的位置。

总体递推最小二乘法算法描述如下：

a)设定干扰源位置初始估计值最右奇异向量的初始值v_m+1(0)和矩阵F的初始值F(0)；

b)通过量测数据，更新矩阵F，其中：F(k)表示矩阵F的第k个递推，d_k＝[sinφ_k -cosφ_k x_k sinφ_k-y_k cosφ_k]^T，φ_k、(x_k,y_k)分别为飞机在第k个测向点测得的方位角和飞机的位置；

c)更新最右奇异向量v_m+1(k)，w(k)＝F(k)v_m+1(k-1)，其中，v_m+1(k)表示最右奇异向量的第k个递推，表示v_m+1(k)的估计值；

d)获取干扰源粗略位置：其中，为第k次递推的位置估计值，表示的第i个分量(i∈[1,m+1])；

步骤二：扩展卡尔曼滤波

在假设系统为高斯白噪声的前提下，基于总体最小二乘法的新息迭代卡尔曼滤波算法如下：

a)将步骤一中求得的目标位置估计值与飞机X_o的相对位置作为EKF的初始值设定估计误差协方差阵P_0|0，过程噪声方差Q₀和量测噪声方差R₀，迭代次数为N；

b)k＝1时，

(1)计算状态一步预测其中为状态向量X_k-1的估计值，u_k-1为k-1次递推的控制输入，A∈R^n×n、B∈R^n×n为状态向量和控制输入向量的常值系数矩阵，为状态向量X_k-1的一步预测估计值；

(2)计算协方差一步预测P_k|k-1＝AP_k-1|k-1A^T+Q₀，其中P_k-1|k-1为估计值的误差协方差阵，P_k|k-1为的误差协方差阵；

(3)计算量测预测值其中y(·)为观测方程的函数形式，为观测值的一步预测；

(4)计算增益其中为观测函数的一阶偏导数；

(5)计算状态更新方程其中，Y_k为第k次观测量的值；

(6)计算协方差更新方程I为单位矩阵；

c)k＝[2,N]时,

在每一步进行迭代，i＝[1,M]；

(1)计算状态一步预测

(2)计算协方差一步预测P_k|k-1i＝AP_k-1|k-1,iA^T+Q₀；

(3)计算量测预测值

(4)计算增益其中，

(5)计算状态更新方程

(6)计算协方差更新方程

最终，得到相对位置的精确估计值得到干扰源的精确位置X_oN为飞机的位置。

本发明应用于飞行校验，利用机载无线电侦测设备和由机载设备获取的飞机飞行参数，采用单站交叉定位，通过一种自适应卡尔曼滤波算法对地面各类无线电干扰源进行实时监测和定位。

本发明是一种自适应干扰源定位飞行校验方法，现有民航干扰侦测的一般流程如图1所示，本发明主要研究的是利用校验飞机进行干扰侦测时数据处理和干扰源解算的问题。如图2所示，由本发明的干扰源定位方法的总流程图，主要包括以下步骤：

步骤100、配置监听接收机，监测是否存在干扰，如果存在干扰，通过监听接收机获得信号的调制方式、频谱和中心频率信息。

步骤200、根据读取的监听接收机参数，根据测向接收机手册配置测向接收机，并对干扰信号进行跟踪，获取干扰信号的来波方位角、信号质量观测量。

步骤300、对观测量进行选择，并导入来自GPS接收机、惯性导航系统、VOR、DME、大气数据系统的飞机飞行参数：飞机位置、速度、航向、姿态信息。

自适应卡尔曼滤波算法对干扰源进行实时定位解算流程如图3所示，有如下步骤：

步骤301、建立以地心坐标系，并将参数转换到地心坐标系之下。

步骤302、设定目标位置初始估计值和最右奇异向量的初始值v_m+1(0)、矩阵F的初始值P(0)及迭代次数。

步骤303、由量测数据，更新矩阵F，其中d_k＝[sinφ_k -cosφ_k x_k sinφ_k-y_k cosφ_k]^T，φ_k、(x_k,y_k)分别飞机在第k个测向点测得的方位角和飞机的位置；更新最右奇异向量v_m+1(k)，w(k)＝F(k)v_m+1(k-1)，估计目标位置得到总体最小二乘法下干扰源的位置

步骤304、将步骤303求得的目标位置估计值与飞机位置求向量差后作为EKF的初始值设定估计均方误差P_0|0，过程噪声方差Q₀和量测噪声方差R₀。

步骤305、k＝1时，计算状态一步预测计算协方差一步预测P_k|k-1＝AP_k-1|k-1A^T+Q₀；计算量测预测值计算增益计算状态更新方程计算协方差更新方程

步骤306、k＝[2,N]时,i＝[1,M]，计算状态一步预测计算协方差一步预测P_k|k-1i＝AP_k-1|k-1,iA^T+Q₀；计算量测预测值计算增益计算协方差更新方程计算状态更新方程计算状态更新方程

步骤307、经过N次测向定位，得到相对位置的精确估计值最终得到干扰源的精确位置

步骤400、输出定位数据及定位误差，结束飞行，生成侦测报告。

Claims (1)

1.一种自适应干扰源定位飞行校验方法，包括以下几个步骤：

步骤一：通过总体递推最小二乘法获得干扰源粗略位置；

具体的：

(1)设定干扰源位置初始估计值最右奇异向量的初始值v_m+1(0)和矩阵F的初始值F(0)；

(2)通过量测数据，更新矩阵F，其中：F(k)表示矩阵F的第k个递推，d_k＝[sinφ_k -cosφ_k x_ksinφ_k -y_kcosφ_k]^T，φ_k、(x_k,y_k)分别为飞机在第k个测向点测得的方位角和飞机的位置；

(3)更新最右奇异向量v_m+1(k)，w(k)＝F(k)v_m+1(k-1)，其中，v_m+1(k)表示最右奇异向量的第k个递推，表示v_m+1(k)的估计值；

(4)获取干扰源粗略位置：其中，为第k次递推的位置估计值，表示的第i个分量(i∈[1,m+1])；

步骤二：将获得的干扰源粗略位置作为新息迭代扩展卡尔曼滤波的初始值进行求解，得到最终的干扰源位置；

设系统为高斯白噪声，基于总体最小二乘法的新息迭代卡尔曼滤波算法如下：

(1)将步骤一中求得的干扰源粗略位置与飞机位置X_o的相对位置作为EKF的初始值设定估计误差协方差阵初值P_0|0，过程噪声方差Q₀和量测噪声方差R₀；设定迭代次数为N；

(2)当k＝1时，

①计算状态一步预测其中为状态向量X_k-1的估计值；u_k-1为k-1次递推的控制输入，A∈R^n×n、B∈R^n×n为状态向量和控制输入向量的常值系数矩阵，为状态向量X_k-1的一步预测估计值；

②计算协方差一步预测P_k|k-1＝AP_k-1|k-1A^T+Q₀，其中P_k-1|k-1为估计值的误差协方差阵；P_k|k-1为的误差协方差阵；

③计算量测预测值其中y(·)为观测方程的函数形式，为观测值的一步预测；

④计算增益其中为观测函数的一阶偏导数；；

⑤计算状态更新方程其中，Y_k为第k次观测量的值；

⑥计算协方差更新方程I为单位矩阵；

(3)当k＝[2,N]时,

在每一步进行迭代，i＝[1,M]；

①计算状态一步预测

②计算协方差一步预测P_k|k-1i＝AP_k-1|k-1,iA^T+Q₀；

③计算量测预测值

④计算增益其中，⑤计算状态更新方程

⑥计算协方差更新方程

最终，得到相对位置的精确估计值得到干扰源的精确位置X_oN为飞机的位置。