基于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)为码减载波差分值，也称为卫星信号伪码-载波偏离度；
步骤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表示以自然常数为底的指数函数，参数α用于控制sigmoid函数的形状，其取值决定函数曲线上升速度，参数β用于控制sigmoid函数的取值范围。
根据大量实验，在保证最优滤波性能基础上，分别设置α＝1，β＝2。
步骤3中，根据大量实验，在保证完好性虚警概率最小基础上，利用高斯膨胀法确定膨胀系数后，设置检测门限值为0.02m/s。
有益效果：
1、本发明提出了在进行卫星信号电离层异常检测前首先对计算得到的卫星信号伪码-载波偏离度进行平滑预处理，从而提高了检测灵敏度；
2、本发明提出了一种改进的sigmoid函数，利用改进的sigmoid函数调节步长因子与预测误差之间的映射关系，lms自适应滤波器能够同时兼顾稳态误差和收敛速度两项性能指标；
3、本发明提出了一种基于lms自适应滤波的卫星信号伪码-载波偏离度高频噪声抑制方法，通过对输入卫星信号的码减载波差分值进行自适应滤波以获得较为准确的电离层延时变化率，提高了电离层异常检测精度与灵敏度。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述或其他方面的优点将会变得更加清楚。
图1是基于lms自适应滤波的gbas电离层异常检测方法流程。
图2是改进sigmoid函数和sigmoid函数对比曲线。
图3是单通道lms自适应滤波器结构图。
图4a为基于二阶线性时不变低通滤波的电离层延时估计曲线。
图4b为基于lms自适应滤波的电离层延时估计曲线。
图5是不同方法检测响应时间。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
如图1所述，本发明提供了一种基于lms自适应滤波的gbas电离层异常检测方法，图1为基于lms自适应滤波的gbas电离层异常检测方法流程，包括以下步骤：
步骤1，地基增强系统(gbas)地面站接收机输出伪距和载波相位观测值计算伪码-载波偏离度，并对计算得伪码-载波偏离度进行平滑预处理，抑制伪距测量、多路径误差等噪声；
步骤2，以步骤1中计算得到的卫星信号伪码-载波偏离度作为检测滤波器的输入，利用最小均方(lms)自适应滤波器抑制伪码-载波偏离度高频噪声，以获取被检测卫星信号电离层延时变化率
步骤3，将步骤2中检测得到的卫星信号电离层延时变化率与检测门限值进行比较，如果电离层延时变化率大于检测门限值，说明电离层梯度异常，需要隔离异常卫星，如果电离层延时变化率不大于检测门限值，说明电离层较为平稳，保留该颗卫星。
步骤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基准站gnss接收机输出的伪距观测值，φ(k)为k时刻gbas基准站gnss接收机输出的载波相位观测值，ρif(k)为k时刻无电离层影响的伪距，r(k)为k时刻星站真实的几何距离，c为光速，δtu(k)和δts(k)分别为k时刻接收机钟差和卫星钟差，i(k)和t(k)分别为k时刻电离层延迟和对流层延迟，λ为波长，n为整周模糊度，ερ(k)和εφ(k)分别为k时刻接收机伪距和相位观测误差，包括多路径误差等。
对码减载波进行差分是一种常用的评估电离层梯度大小方法，根据如下公式进行码减载波计算为：
cmc(k)＝ρ(k)-φ(k)
＝2i(k)+ερ(k)-εφ(k)-λ·n，
其中，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)为码减载波差分值，也称为卫星信号伪码-载波偏离度；
图2为改进sigmoid函数和sigmoid函数对比曲线，由图2可以看出，当误差为2时，sigmoid函数与改进sigmoid函数曲线都已经达到或者接近步长因子的最大值，说明此时滤波算法的收敛速度较快，有利于滤波输出结果快速向误差减小的方向变化。对比sigmoid函数与改进sigmoid函数曲线可知，当系统接近收敛时，改进sigmoid函数步长因子变化更为缓慢，引起的振荡较小，因此系统稳定性能更好。
图3为单通道lms自适应滤波器结构，由图3可以得到步骤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表示以自然常数为底的指数函数，参数α用于控制sigmoid函数的形状，其取值决定函数曲线上升速度，参数β用于控制sigmoid函数的取值范围。
根据大量实验，在保证最优滤波性能基础上，分别设置α＝1，β＝2。
步骤3中，根据大量实验，在保证完好性虚警概率最小基础上，利用高斯膨胀法确定膨胀系数后，设置检测门限值为0.02m/s。
实施例
下面以中国电科某研究所建立的试验型gbas基准站为实验对象，验证本发明提出的gbas电离层异常检测方法的正确性。该试验型gbas地面系统4个基准站、数据处理中心、多载体导航信号模拟器、vdb(veryhighfrequency,vhfdatabroadcasting)电台以及卫星天线。
利用真实卫星信号进行验证实验，采集了多组实验数据，选择其中一组数据进行了详细分析，实验时间为2018年3月21日。由于电离层异常出现频率较低，因此本文在14号卫星伪距与载波相位观测值上仿真一段电离层异常情形，假设在整个观测时段8000s至9000s之间存在梯度值为0.024m/s的电离层异常现象，持续时间为1000s。
图4a使用二阶线性时不变低通滤波器对伪码-载波偏离度进行滤波，图4b使用本发明所提出的lms自适应滤波方法进行滤波。
为进一步验证本发明所提出自适应滤波方法的有效性，在14号卫星观测时段8000s至9000s秒之间模拟不同大小的电离层梯度，电离层梯度分别设置为0.01m/s、0.015m/s、0.02m/s、0.025m/s、0.03m/s、0.035m/s以及0.04m/s，然后利用已有检测方法以及本发明所提出的方法进行异常检测，检测响应时间如图5所示。
本发明提供了基于lms自适应滤波的gbas电离层异常检测方法，具体实现该技术方案的方法和途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
- 还没有人留言评论。精彩留言会获得点赞！