基于LMS自适应滤波的GBAS电离层异常检测方法技术

本发明专利技术公开了基于LMS自适应滤波的GBAS电离层异常检测方法，本发明专利技术通过构造单通道变步长LMS自适应滤波器以抑制伪码‑载波偏离度高频噪声，单通道LMS自适应滤波器是在标准双通道LMS自适应滤波器基础上，利用被检测信号短时相关性及其量化噪声的非相关性，构造一个采用被检测信号延时量作为参考输入的自适应滤波器，同时对Sigmoid函数进行改进，使得自适应滤波器在前期收敛速度快，且待滤波器收敛后保持较高稳定性。本发明专利技术所提出的滤波方法具有更好的高频噪声抑制效果，对电离层慢增长故障检测性能最优。

GBAS ionospheric anomaly detection method based on LMS adaptive filtering

【技术实现步骤摘要】
基于LMS自适应滤波的GBAS电离层异常检测方法
本专利技术属于星基导航增强系统领域，具体涉及基于LMS自适应滤波的GBAS电离层异常检测方法。
技术介绍
目前民用航空导航中普遍使用仪表着陆系统(InstrumentLandingSystem，ILS)为飞机提供精密进近与着陆引导服务，其频道仅为40，严重制约了进近航班量的提升，逐渐成为限制航班量增长的瓶颈。近年来，全球导航卫星系统(GlobalNavigationSatelliteSystem，GNSS)逐步完善，GNSS凭借全球、全天候、高精度等优点得到越来越多的应用，国际民航组织(InternationalCivilAviationOrganization，ICAO)已计划利用GNSS替代ILS，降低航空导航成本、增加航路设计和进场路线的灵活性。然而由于GNSS信号在传输过程中受到电离层延时、对流层延时、多路径误差等影响，基于信号本身的测距精度不能满足ICAO对民用航空进近所需导航性能要求，因此需要对GNSS信号进行增强，其中地基增强系统(ground-basedaugmentationsystems，GBAS)通过建立位置点已知的基准站，采用差分和完好性监视技术，提高了系统完好性和可用性。2009年，Honeywell公司研制的基于全球定位系统(GlobalPositioningSystem，GPS)L1频点的GBAS获得ICAOCATI类许可认证。实验表明，单频点GBAS不能满足CATIII类导航需求，为此相关研究机构起草了多星座、多频点CATIII类GBAS研究规划。目前国内中电20所、中电54所、中电28所等研究机构正在开展CATI类原型系统研制工作，围绕伪距差分校正、卫星信号完好性监测以及机载差分定位等展开研究，北京航空航天大学长期以来一直致力于机载完好性评估与仿真关键技术研究，积极推动了该领域技术进步。电离层是地球大气层中的一个电离区域，卫星信号通过该区域后受到其中自由电子和离子的非线性散射影响，信号传播速度与路径发生了变化，从而会降低导航定位精度，该误差也称为电离层延时，是GNSS主要误差源之一。尽管电离层延时对GNSS精度影响可以通过数学建模或使用双频GNSS卫星信号组合加以削弱，但是在特定使用领域中，由于电离层延时引起的位置误差严重时会影响用户的生命安全，因此必须实时监测电离层状态，例如，用于支持飞机精密进近与着陆的GBAS受到电离层梯度影响后难以满足其导航精度与完好性需求。对于GBAS而言，电离层异常检测精度和灵敏度较为苛刻且具有挑战性，目前如何快速、准确地检测异常电离层是GBAS电离层监测的重要研究方向之一。
技术实现思路
专利技术目的：针对现有技术的不足，本专利技术提出一种基于LMS自适应滤波的GBAS电离层异常检测方法，通过构造单通道变步长最小均方自适应滤波器以抑制被检测卫星信号伪码-载波偏离度高频噪声，从而能够获得较为准确的电离层延时变化率，并与设定的电离层异常检测门限值进行比较，如果检测得到的延时变化率大于异常检测门限值，则说明卫星信号异常。技术方案：本专利技术所述的基于LMS自适应滤波的GBAS电离层异常检测方法，包括以下步骤：步骤1，根据地基增强系统(GBAS)地面站接收机输出的伪距和载波相位观测值计算卫星信号伪码-载波偏离度，并对卫星信号伪码-载波偏离度进行平滑预处理，抑制伪距测量、多路径误差等噪声；步骤2，以步骤1中计算得到的卫星信号伪码-载波偏离度作为检测滤波器的输入，利用单通道变步长最小均方(LMS)自适应滤波器抑制卫星信号伪码-载波偏离度高频噪声，以获取被检测卫星信号电离层延时变化率；步骤3，将步骤2中检测得到的卫星信号电离层延时变化率与检测门限值进行比较，如果卫星信号电离层延时变化率大于检测门限值，说明电离层梯度异常，需要隔离异常卫星，如果卫星信号电离层延时变化率不大于检测门限值，说明电离层较为平稳，保留该颗卫星。步骤1-1，计算被检测卫星信号码减载波值；步骤1-2，由码减载波值计算得卫星信号伪码-载波偏离度，并进行平滑预处理。步骤1-1包括：进行如下计算：ρ(k)＝ρIF(k)+I(k)+ερ(k)，φ(k)＝ρIF(k)-I(k)+λ·N+εφ(k)ρIF(k)＝r(k)+c[δtu(k)-δts(k)]+T(k)，其中，ρ(k)为k时刻GBAS基准站卫星信号接收机输出的伪距观测值，φ(k)为k时刻GBAS基准站卫星信号接收机输出的载波相位观测值，ρIF(k)为k时刻无电离层影响的伪距，r(k)为k时刻星站真实的几何距离，c为光速，δtu(k)和δts(k)分别为k时刻接收机钟差和卫星钟差，I(k)和T(k)分别为k时刻电离层延迟和对流层延迟，λ为波长，N为整周模糊度，ερ(k)和εφ(k)分别为k时刻接收机伪距和相位观测误差；然后计算得到码减载波值。步骤1-1中，根据如下公式计算码减载波值：其中，CMC(k)为k时刻码减载波值。步骤1-2包括：设定一段时间T1内(此时间段T1不固定，需要根据伪距和载波相位输出值进行判别，一般为数小时不等，一旦接收机失锁，则N值会发生跳变，此时一般需要重新初始化dz值)，接收机处于锁定状态，无失锁和失周现象，因而N值大小不变，同时令δε(k)＝ερ(k)-εφ(k)，即δε(k)表示k时刻ερ(k)和εφ(k)的差值，则有：其中，δε(k-(k0+1))表示k-(k0+1)时刻ερ(k-(k0+1))和εφ(k-(k0+1))的差值，dz(k)为码减载波差分值，也称为卫星信号伪码-载波偏离度；为电离层延时速率；Ts＝1s，表示相邻两个历元间隔时间；k0＝20，表示平滑窗口时间长度。步骤2包括如下步骤：步骤2-1，初始化W(M)＝0；步骤2-2，计算输出值：k＝M+1,M+2,…，y(k)＝WT(k)X(k)；步骤2-3，计算预测误差：e(k)＝dz(k)-y(k)；步骤2-4，调整步长因子μ(k)；步骤2-5，更新权重：其中，k表示观测历元时间；M为LMS自适应滤波窗口长度；y(k)为滤波器输出值，即被检测卫星信号电离层延时变化率；X(k)表示LMS自适应滤波器k时刻输入信号矢量，X(k)＝[dz(k-1),dz(k-2),…,dz(k-M)]T，其中，dz(k-1),dz(k-2),...,dz(k-M)分别表示k时刻码减载波差分值延时1至M个历元输入量；W(k)表示LMS自适应滤波器k时刻权重更新系数，W(k)＝[w1,k,w2,k,…,wM,k]T，其中，wM,k表示k时刻LMS自适应滤波器窗口长度为M时的第M个权重系数值，WT(k)为W(k)的转置，μ(k)表示k时刻单通道变步长最小均方滤波器步长因子，e(k)表示k时刻预测误差。步骤2-4中，采用改进的Sigmoid函数调节步长因子μ与预测误差e之间的映射关系：其中，exp本文档来自技高网...

【技术保护点】
1.基于LMS自适应滤波的GBAS电离层异常检测方法，其特征在于，包括以下步骤：/n步骤1，计算卫星信号伪码-载波偏离度，并对卫星信号伪码-载波偏离度进行平滑预处理；/n步骤2，以卫星信号伪码-载波偏离度作为检测滤波器的输入，利用单通道变步长最小均方自适应滤波器抑制卫星信号伪码-载波偏离度高频噪声，以获取被检测卫星信号电离层延时变化率；/n步骤3，将步骤2得到的卫星信号电离层延时变化率与检测门限值进行比较，如果卫星信号电离层延时变化率大于检测门限值，说明电离层梯度异常，需要隔离异常卫星，如果卫星信号电离层延时变化率不大于检测门限值，说明电离层平稳，保留该颗卫星。/n

【技术特征摘要】
1.基于LMS自适应滤波的GBAS电离层异常检测方法，其特征在于，包括以下步骤：
步骤1，计算卫星信号伪码-载波偏离度，并对卫星信号伪码-载波偏离度进行平滑预处理；
步骤2，以卫星信号伪码-载波偏离度作为检测滤波器的输入，利用单通道变步长最小均方自适应滤波器抑制卫星信号伪码-载波偏离度高频噪声，以获取被检测卫星信号电离层延时变化率；
步骤3，将步骤2得到的卫星信号电离层延时变化率与检测门限值进行比较，如果卫星信号电离层延时变化率大于检测门限值，说明电离层梯度异常，需要隔离异常卫星，如果卫星信号电离层延时变化率不大于检测门限值，说明电离层平稳，保留该颗卫星。


2.根据权利要求1所述的方法，其特征在于，步骤1包括：
步骤1-1，根据如下公式计算被检测卫星信号码减载波值：
ρ(k)＝ρIF(k)+I(k)+ερ(k)，
φ(k)＝ρIF(k)-I(k)+λ·N+εφ(k)
ρIF(k)＝r(k)+c[δtu(k)-δts(k)]+T(k)，
其中，ρ(k)为k时刻GBAS基准站卫星信号接收机输出的伪距观测值，φ(k)为k时刻GBAS基准站卫星信号接收机输出的载波相位观测值，ρIF(k)为k时刻无电离层影响的伪距，r(k)为k时刻星站真实的几何距离，c为光速，δtu(k)和δts(k)分别为k时刻接收机钟差和卫星钟差，I(k)和T(k)分别为k时刻电离层延迟和对流层延迟，λ为波长，N为整周模糊度，ερ(k)和εφ(k)分别为k时刻接收机伪距和相位观测误差；
根据如下公式进行码减载波计算：



其中，CMC(k)为k时刻码减载波值。
步骤1-2，由码减载波值计算得卫星信号伪码-载波偏离度，并进行平滑预处理：
设定一段时间T1内，接收机处于锁定状态，无失锁和失周现象，则N值大小不变，同时令δε(k)＝ερ(k)-εφ(k)，即δε(k)表示k时刻ερ(k)和εφ(k)的差值，则有：



其中，δε(k-(k0+1))表示k-(...

【专利技术属性】
技术研发人员：胡杰，严勇杰，谢晓妤，石潇竹，李贺，单尧，
申请(专利权)人：中国电子科技集团公司第二十八研究所，
类型：发明
国别省市：江苏;32