CN114296117B - 一种精密进近下地基增强系统完好性评估方法及系统 - Google Patents

一种精密进近下地基增强系统完好性评估方法及系统 Download PDF

Info

Publication number
CN114296117B
CN114296117B CN202210226260.5A CN202210226260A CN114296117B CN 114296117 B CN114296117 B CN 114296117B CN 202210226260 A CN202210226260 A CN 202210226260A CN 114296117 B CN114296117 B CN 114296117B
Authority
CN
China
Prior art keywords
simulation
monitor
airplane
pseudo
parameters
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210226260.5A
Other languages
English (en)
Other versions
CN114296117A (zh
Inventor
刘海颖
孙颢
李铁香
徐毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Applied Mathematics Center
Original Assignee
Nanjing Applied Mathematics Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Applied Mathematics Center filed Critical Nanjing Applied Mathematics Center
Priority to CN202210226260.5A priority Critical patent/CN114296117B/zh
Publication of CN114296117A publication Critical patent/CN114296117A/zh
Application granted granted Critical
Publication of CN114296117B publication Critical patent/CN114296117B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种精密进近下地基增强系统完好性评估方法,本方法利用蒙特卡洛采样进近过程中的参数,并仿真飞机精密进近过程中的GBAS系统,利用设计的多个监视器对于进近过程中的参数进行检验,最终得到单个监视器的漏检率并进行组合,通过加权最终得到飞机精密进近过程中的完好性风险评估结果。本方法能够模拟极限异常情况下的飞机精密进近过程,消除潜在的完好性风险。

Description

一种精密进近下地基增强系统完好性评估方法及系统
技术领域
本发明涉及卫星导航技术领域,特别涉及一种精密进近下地基的地基增强系统完好性评估方法及系统。
背景技术
为每个可见星提供测距误差校正信息,这些校正极大地提高了飞机的定位精度,并支持精密进近和自动着陆。在飞机精密进近过程中,电离层异常事件对于飞机的定位及精密进近过程有着重大的影响。
GBAS(地基增强系统)的工作原理是通过地面参考站的形式将差分改正信息和误差信息通过甚高频数据广播播发给用户,从而对用户的完好性进行增强。在GBAS(地基增强系统)开发过程中,已经评估过异常电离层梯度对 GBAS 的影响。CAT-I GBAS(地基增强系统一类精密进近)中的现有技术涉及进行详尽的模拟,以在最坏情况的电离层参数和卫星、飞机的最坏几何形状的所有可能组合中搜索最大的差分距离误差。 但是在CAT-I GBAS仿真过程中并未涉及危险的电离层事件。电离层是部分电离的大气区域,完全电离的大气区域称磁层。电离层的范围从离地面约50公里开始一直伸展到约1000公里高度的地球高层大气空域,且其中存在相当多的自由电子和离子,能使无线电波改变传播速度,发生折射、反射和散射,而由于卫星都处于大气层外,卫星信号传播到地球上必须要经过大气层中的电离层,于是必然要受到电离层的影响。由于电离层在不同位置的水平分布较为不均匀,受时间影响的变化复杂,且还存在一定的异常变化特征,给GNSS(全球导航卫星系统)性能提升和建立全球或区域范围内的电离层异常模型带来了挑战。虽然这种情况在概率上微不足道,但其存在会导致非常保守的评估。因此,需要对异常梯度造成的电离层威胁进行更真实的描述。
完好性反映了导航系统在不能用于导航服务时及时向用户提出告警的能力,作为地基增强系统的重要性能指标,完好性评估与风险识别的可用性仍有很大的研究空间,本发明利用蒙特卡洛采样仿真系统模拟极端情况下的飞机精密进近过程,以评估在极端情况中地基增强系统在最严格的运作标准下完好性指标能否满足要求。
发明内容
为了解决上述问题,本发明公开了一种精密进近下地基增强系统完好性评估系统与方法,利用蒙特卡洛采样方法和多完好性监视器组合,以飞机精密进近过程为背景,解决了现有技术中对潜在完好性风险识别能力差的问题。
技术方案如下,一种精密进近下地基增强系统完好性评估方法,包括如下步骤:
步骤1,设置仿真飞机进近参数:包括固定参数与随机参数,固定参数通过系统设定,随机参数通过蒙特卡洛采样得到,由随机参数组成电离层仿真参数;
步骤2,基于进近参数进行地基增强系统下的飞机精密进近仿真:确定飞机的相对位置,测量得到伪距,并通过差分校正该伪距得到伪距误差,判断电离层是否对伪距误差有影响,并最终确认进近飞机的伪距误差;
步骤3,基于步骤1的仿真飞机进近参数及步骤2的仿真结果,进行完好性参数仿真,仿真期间内每秒生成监视器监测量/检验统计量;
步骤4,当信号异常时监视器的监测量/检验统计量会超出预设的阈值,将监视器的监测量/检验统计量与阈值做比较以得到监测超限情况,将该超限情况结合完好性参数仿真结果得到漏检概率P md ,即记录下对应于每组完好性参数下的漏检时间占单次精密进近仿真时间的百分比;
步骤5,将各个监视器的漏检概率P md 进行组合得到组合漏检概率P com ,以消除各个监视器间的相关性,并利用组合漏检概率P com 解出危险误导信息概率P HMI 以判断在飞机精密进近过程中的完好性。
本发明还公开了一种精密进近下地基增强系统完好性评估系统,包括参数仿真模块和测量模块,所述参数仿真模块对飞机进行精密进近参数仿真采样,测量模块整合现场采集的辅助数据,提前测量辅助数据并将其加入到仿真过程中以提升仿真系统的适用性与可靠性;基于测量模块测量的数据进行精密进近参数仿真采样,并将仿真采样结果发送到差分处理模块和监视器模块;所述差分处理模块对飞机精密进近过程差分定位仿真,所述监视器模块在仿真过程中采集每个监视器的测试参数,将每个监视器在进近期间漏检概率的最小值发送至数据记录模块进行记录,最终将采集到的每个监视器的漏检概率进行组合,以计算在极限情况下飞机精密进近过程中的GBAS的完好性风险。
有益效果:
(1)本发明公开的方法通过仿真场景参数随机采样与飞机精密进近仿真相结合,并通过测距误差监测和完好性监视器来判断飞机精密进近过程中异常电离层事件导致的完好性风险,探索了GBAS在精密进近过程中极限情况下的完好性,以满足最严格的航空运作要求;
(2)本发明利用蒙特卡洛仿真,在地基增强系统下飞机精密进近的各个环节都进行了参数采样与仿真,以模拟在最极端的异常事件影响下地基增强系统提供服务的能力,为地基增强系统的实际运作提供理论基础;
(3)本发明通过组合多个监测器在大量仿真后得到的漏检概率,可以得到比传统方法更加有效可靠的危险误导概率信息,这大大提升了计算结果的可用性,提高了仿真系统的可靠性。
附图说明
图1为本发明一个实施例的精密进近下地基增强系统完好性评估方法的流程图;
图2为本发明一个实施例的电离层波峰侧视图;
图3为本发明一个实施例的精密进近下地基增强系统完好性评估系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明公开了一种地基增强系统精密进近下完好性评估方法及系统,通过仿真场景参数随机采样结合飞机精密进近仿真,最终通过测距误差监测和完好性监视器来判断在飞机精密进近过程中异常电离层事件导致的完好性风险。
由于电离层在不同位置的水平分布较为不均匀,受时间影响的变化复杂,且还存在一定的异常变化特征,因此电离层的此类特征给GBAS性能提升和建立全球和区域内的电离层异常梯度模型带来了困难。在本发明中,如图1所示,本专利使用的电离层模型为有线性变化的楔形锋面,该电离层模型为以某个采样得到的速度移动的电离层封面模型的一种。
如图1所示,一种地基增强系统精密进近下完好性评估方法及系统,包括如下步骤:
步骤1,设置仿真飞机进近参数,包括站点参数、完好性监视器参数和飞机运动参数,并将其分为固定参数与随机参数,在仿真系统中事先设定固定参数,通过蒙特卡洛采样得到随机参数,由该随机参数组成电离层仿真参数。
其中,站点参数包括:飞机进近机场编号、GBAS地面站设备与着陆阈值点的距离、电离层仿真参数(电离层锋面宽度、坡度与速度、锋面移动方向)、卫星星座仿真;
完好性监视器参数包括:CCD监测平滑时间常数、相关对码片间距半长、相关峰最小检测误差的拟合参数;飞机运动参数包括:飞机进场速度、飞机进近跑道方向、飞机着陆方向角;
飞机运动参数包括:飞机进场速度、飞机进近跑道方向、飞机着陆方向角。
固定参数与随机参数分配如表1和表2所示,
Figure 954052DEST_PATH_IMAGE001
随机参数由蒙特卡洛采样得到,蒙特卡洛采样步骤为:
步骤1.1,确认仿真时间:蒙特卡洛仿真每一分钟提供一次卫星构型变化数据,一天内共有1440个历元。
步骤1.2,进行卫星星座仿真,利用仿真卫星对飞机进行定位。
步骤1.3,仿真飞机进近方向:从-90°到90°每隔10°模拟飞机的进近方向。
步骤1.4,仿真飞机进近参数中的飞机运动参数、站点参数以及完好性监视器参数均匀分布并随机组合生成异常事件威胁模型,并基于步骤1.1及步骤1.2确认的仿真时间及进近方向对该模型进行随机采样。
Figure 654548DEST_PATH_IMAGE002
采用均匀分布方法生成随机参数并随机组合,其中由于飞机的电离层穿刺点(IPP) 和地面设施在 GBAS 服务区域内在空间上都很近,则可以假设飞机的电离层穿刺点(IPP)与地面电离层穿刺点 (IPP)的相对速度是飞机的进近速度。
步骤2,进行地基增强系统下的飞机精密进近仿真,仿真过程包括飞机运动参数仿真、差分校正过程、误差计算过程、完好性监测过程,具体为:
步骤2.1,通过仿真参数来确定飞机的位置:利用步骤1中的随机参数中的电离层仿真参数解出电离层的位置,并利用采样得到的飞机进近机场编号、进近跑道方向以及事先设定的地面参考站位置,并进行卫星星座仿真,结合该星座仿真结果得到着陆阈值点(LTP),通过反推得到飞机初始仿真位置,此时地面站、飞机、电离层的相对位置便通过仿真得到。
步骤2.2,设定飞机进近过程为抛物线轨迹,基于步骤2.1确定的飞机位置,测量得到载波相位测量值和伪距,并通过载波相位测量的变化值平滑伪距,以减少伪距中的快变误差。
在GBAS中,载波相位测量值平滑伪距的方法用于减轻飞机和地面设施中接收机噪声和多径误差的影响,其本质上是利用载波相位的高精度对于伪距测量值进行辅助,影响伪距测量值。由
Figure 231023DEST_PATH_IMAGE003
表示的载波平滑后的伪距被定义为:
Figure 157391DEST_PATH_IMAGE004
其中,ρ(t)φ(t)分别表示给定历元t的单频伪距和载波相位测量值,,M=τ/T其中τ为载波平滑滤波时间常数,T为原始观测值的采样间隔,通常为0.5秒。在该仿真中,假设伪距和载波相位测量值仅由电离层延迟造成,且飞机和地面设施都使用30秒的平滑时间常数。
对接收到的伪距和载波相位测量值同已知的真实距离做差分,分别计算平滑伪距校正值和载波相位校正值。
步骤2.3,通过蒙特卡洛仿真得到的电离层仿真参数计算电离层对于伪距误差的影响,在仿真中,异常电离层梯度以GBAS的标准方式建模为楔形波峰。如图2所示,为一架飞机和一个朝向GBAS地面设施移动的电离层波峰的侧视图。波峰斜率代表电离层延迟从高到低的线性变化。梯度和电离层锋面宽度的乘积为总的电离层延迟。具体为:确定电离层穿刺点,根据电离层穿刺点的位置决定其对于伪距测量是否有影响,若穿刺点位置在异常电离层范围内,则需要根据电离层梯度及锋面情况计算异常电离层引起的伪距误差并与步骤2.2中的校正误差结合作为进近飞机的伪距误差。
步骤3,进行完好性参数仿真,基于上述参数,在仿真期间内每秒生成四个监视器的测距误差、检验统计量。
步骤4,当信号异常时监视器的监测量/检验统计量会超出预设的阈值,将监视器的监测量/检验统计量与阈值做比较以得到监测超限情况。在每个仿真历元,计算每个监视器的监测量超限情况并记录在LTP点的差分距离误差,最终结合仿真情况,在约1010数量级的仿真样本得到每个监视器的漏检概率P md
四个监视器包括码/载波分歧(CCD)监视器、载波平滑码更新(Inno)监测器、相关峰(correlation peak)监测器和双平滑电离层梯度(DSIGMA)监视器。具体监测过程如下:
码/载波分歧(CCD)监视器:CCD监测通常用于探测电离层风暴是否发生,并且保证所有给定卫星的CCD不发生异常。通常采用几何滑动平均方式来估计CCD,其计算方式如下所示:
Figure 502922DEST_PATH_IMAGE005
式中,CCD (m,n) (k)为码载波分歧值的监测量,其中(m,n)代表对应的第m颗卫星和第n台接收机,τ d 为平均持续时间,是一个值为90秒的常量,T s 是地基增强系统的测量值的更新频率,一般为0.5秒;Z(k)k时刻的码减载波(Code minus carrier, CMC)观测量,通过下式得到:
Z(k)=ρ(k)-φ(k)
Figure 59805DEST_PATH_IMAGE006
其中,ρ(k)表示k时刻的原始伪距观测量,即
Figure 608729DEST_PATH_IMAGE007
时刻未经平滑的原始伪距观测量,φ(k)代表对应于k时刻的载波观测量,M(k)代表多径效应的影响,N(k)代表整周模糊度、λ代表载波波长,
Figure 338788DEST_PATH_IMAGE008
为电离层梯度,由步骤1中仿真得到
Figure 538825DEST_PATH_IMAGE009
通过将计算得到的CCD监测量与阈值相比,可以得到超限的CCD观测量。
相关峰(correlation peak)监测器利用仿真虚拟相关峰来模拟相关峰监测中出现的异常,通过拟合函数模拟相关峰监测的检验统计量来对信号形变异常进行仿真,最终判断潜在的完好性风险。即:利用步骤1中随机生成的拟合参数,结合步骤2中飞机相对位置确定仿真,通过观测数据解算的卫星仰角,通过仿真虚拟相关对来模拟相关峰监测流程中出现的异常情况,通过拟合函数模拟相关峰监测的检验统计量来对信号形变异常进行仿真,最终来判断潜在的完好性风险,具体方式如下:
相关峰的检验统计量为最小可监测误差(Minimum Detective Error, MDE)是指接收机在保证一定的虚警概率(probability of False Alarm, FA)和危险错误信息概率(probability of Hazardly Misleading Information, HMI)条件下,可检测到的最小误差,也是Δ算法的判决阈值。相关峰监测中的最小可监测误差MDE可由如下公式计算:
Figure 328926DEST_PATH_IMAGE010
其中,式中,K ffd 表示由自由虚警概率决定的膨胀系数,K md 表示由漏检概率决定的膨胀系数,飞机在不同等级的精密终端导航和精密进近过程中会有不同的膨胀系数要求。假设飞机接收的卫星信号统计误差满足高斯分布。
Figure 863682DEST_PATH_IMAGE011
表示待测参数的标准差,这里的待测参数是指“从相关对”测量值( 1 2 C-1)与“主相关对”测量值( ref )的差。
相关峰拟合函数是基于信号质量接收机的检验统计量R建模的相关峰最小检测误差MDE关于仰角的函数,在不同码片间距下具有不同的系数,根据函数随机生成当前的MDE即相关峰判决门限,检验统计量R的传统计算方法为:
Figure 131852DEST_PATH_IMAGE012
其中I prompt 为相关时间间隔为0时相关器的输出,对于正常的C/A码,该时间对应的为峰值,I offset 为相关器时间间隔非0时的相关器输出,负号代表滞后,正号代表超前。当信号发生异常时,统计量会超出阈值。
最终对R监测量进行拟合,最终的拟合函数形式为:
Figure 389658DEST_PATH_IMAGE013
Figure 350661DEST_PATH_IMAGE014
表示卫星仰角,R角标代表码片间距半长,式中的系数范围见下表:
Figure 123445DEST_PATH_IMAGE015
通过MDE的拟合函数仿真相关峰监测器的监测值,记录超限部分的阈值用以计算该监视器的P md
载波平滑码更新(Inno)监测器:利用步骤2中计算的伪距和载波相位测量值计算载波平滑码更新监测量,用以检测原始伪距测量中误差,主要检查其中的脉冲和阶跃误差,检验统计量为:
Figure 680459DEST_PATH_IMAGE016
ρ s (k-1)代表(k-1)时刻的平滑后伪距观测量,φ(k)代表k时刻的载波观测量。
双平滑电离层梯度(DSIGMA)监视器:利用设定不同平滑常数的滤波器,在载波平滑伪距过程中产生的校正伪距的差,通过设定阈值对于电离层异常事件进行监测,DSIGMA监测器由两个分别具有平滑时间常数为100秒和30秒的滤波器,其本质为一个二阶的带通滤波器,其检验统计量为:
Figure 855088DEST_PATH_IMAGE017
P c_100 P c_30 分别代表平滑时间常数为100秒和30秒的当前时间的校正伪距。
需要注意的是,监视器中的阈值建模均以高斯分布依据历史数据及RTCA DO-246E标准搭建,在着陆阈值点的伪距误差阈值依据ICAO SARPs(标准及建议措施)设定为2.75m。
步骤5,通过将各个监视器解算出的P md 进行组合以消除各个监视器间的相关性,并利用组合误警率解出危险误导信息P HMI 概率以判断在飞机精密进近过程中的完好性,漏检概率使用如下方式组合:
P com =min[P md,CP ×P md,CCD P md,CP ×P md,DSIGMA P md,CP ×P md,INNO ]
其中P md,-- 代表不同监视器的漏检概率,P md,CP P md,CCD P md,DSIGMA P md,INNO 分别表示相关峰监测器的漏检率、码/载波分歧监视器的漏检率、双平滑电离层梯度监视器的漏检率、载波平滑码更新监测器的漏检率。
P com 加权得到危险误导信息概率P HMI ,具体为:
Figure 986992DEST_PATH_IMAGE018
其中N为仿真参数样本的总数,i为仿真参数样本的部分测距误差超限子集,完好性评估系统会在飞机经过事先设定着陆阈值点测距误差超限的情况下记录下这些样本并用以加权。一旦测距误差超过定义的限值,系统便会记录一次,此时将仿真结束后所有记录到的测距误差超限的仿真次数即为i。
利用该仿真系统得到的危险误导信息概率,可以得到在该环境下的地基增强系统在极端环境下出现危险的概率,通过对比机场性能要求标准,可以评估该地基增强系统在该区域的完好性。
在实际使用过程中,仿真系统利用当地实测的历史数据对于各个监测量建模并设定阈值,并利用蒙特卡洛采样对于各个参数进行设定生成并组合以决定单次精密进近仿真前的参数配置,接下来基于仿真出来的多组参数进行飞机精密进近仿真,记录部分关键测试参数并基于此来解算飞机在单次精密进近过程中的漏检概率,最后将其进行组合以得到飞机危险误导信息概率,进而评估飞机在极端情况下的完好性风险。
如图3所示,本发明还公开了一种地基增强系统精密进近下完好性评估系统,包括参数仿真模块和测量模块,参数仿真模块对飞机进行精密进近参数仿真采样,测量模块整合现场采集的辅助数据,例如机场的具体位置、接收机的观测精度,将这些数据提前测量并将其加入到仿真过程中能够提升仿真系统的适用性与可靠性。基于测量模块测量的参数进行精密进近参数仿真采样,并将仿真数据发送到差分处理模块和监视器模块,差分处理模块对飞机精密进近过程差分定位仿真,监视器模块在仿真过程中采集每个监视器的测试参数,将每个监视器在进近期间漏检概率的最小值发送至数据记录模块进行记录,最终将采集到的每个监视器的漏检概率进行组合,以计算在极限情况下飞机精密进近过程中的GBAS的完好性风险。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种精密进近下地基增强系统完好性评估方法,其特征在于:包括以下步骤:
步骤1,设置仿真飞机进近参数:包括固定参数与随机参数,固定参数通过系统设定,随机参数通过蒙特卡洛采样得到,由随机参数组成电离层仿真参数;
步骤2,基于进近参数进行地基增强系统下的飞机精密进近仿真:确定飞机的相对位置,测量得到伪距,并通过差分校正该伪距得到伪距误差,判断电离层是否对伪距误差有影响,并最终确认进近飞机的伪距误差;
步骤3,基于步骤1的仿真飞机进近参数及步骤2的仿真结果,进行完好性参数仿真,仿真期间内每秒生成监视器监测量/检验统计量;
所述监视器包括码/载波分歧监视器、载波平滑码更新监测器、相关峰监测器和双平滑电离层梯度监视器,其中,所述相关峰监测器利用仿真虚拟相关峰来模拟相关峰监测中出现的异常,以判断潜在的完好性风险,检验统计量为最小可监测误差MDE,对监测到的检验统计量R进行拟合,拟合函数形式为:MDE=a3θ3+a2θ2+a1θ+a0,θ表示卫星仰角,a3、a2、a1、a0为与R相关的系数;
所述双平滑电离层梯度监视器通过设定不同平滑常数的滤波器,在载波相位测量值平滑伪距过程中产生伪距误差,通过设定阈值对于电离层异常事件进行监测,其检验统计量TDSIGMA表示如下:TDSIGMA=Pc_100(t)-Pc_30(t),Pc_100和Pc_30分别代表平滑时间常数为100秒和30秒的当前时间的校正伪距;
步骤4,当信号异常时监视器的监测量/检验统计量会超出预设的阈值,将监视器的监测量/检验统计量与阈值做比较以得到监测超限情况,将该超限情况结合完好性参数仿真结果得到漏检概率Pmd,即在每个仿真历元,计算每个监视器的监测量超限情况并记录在LTP点的差分距离误差,最终结合仿真情况,在约1010数量级的仿真样本得到每个监视器的漏检概率Pmd,记录下对应于每组完好性参数下的漏检时间占单次精密进近仿真时间的百分比;
步骤5,将各个监视器的漏检概率Pmd进行组合得到组合漏检概率Pcom,以消除各个监视器间的相关性,并利用组合漏检概率Pcom解出危险误导信息概率PHMI以判断在飞机精密进近过程中的完好性。
2.根据权利要求1所述的精密进近下地基增强系统完好性评估方法,其特征在于:所述步骤2具体为:
步骤2.1,通过步骤1中的电离层仿真参数解出电离层的位置,并通过采样以得到飞机进近机场编号、进近跑道方向以及事先设定的地面参考站位置,结合卫星星座仿真结果得到着陆阈值点,通过反推得到飞机初始仿真位置,及其与地面站、飞机、电离层的相对位置;
步骤2.2,设定飞机进近过程为抛物线轨迹,基于步骤2.1确定的飞机位置,通过测量得到载波相位测量值和伪距,并通过载波相位测量值平滑伪距,以减少伪距中的快变误差;对接收到的载波平滑后的伪距测量值同已知的真实距离做差分,得到伪距误差;
步骤2.3,根据电离层仿真参数确定电离层穿刺点,根据该穿刺点的位置判断其对于伪距测量是否有影响,若穿刺点位置在异常电离层范围内则判断为有影响,根据电离层梯度及锋面情况计算异常电离层引起的伪距误差并与步骤2.2中的校正误差结合作为进近飞机的伪距误差,否则便不考虑电离层对于定位过程的影响,将步骤2.2中的校正误差作为进近飞机的伪距误差。
3.根据权利要求2所述的精密进近下地基增强系统完好性评估方法,其特征在于:所述码/载波分歧监视器用于监测电离层爆发,以保证伪码测距和载波相位测距的差异在限定范围内,其计算方式如式:
Figure FDA0003646168390000021
式中CCD(k)为码载波分歧值的监测量,τd为平均持续时间,Ts为地基增强系统的测量值的更新频率,Z(k)为k时刻的码减载波观测量。
4.根据权利要求3所述的精密进近下地基增强系统完好性评估方法,其特征在于:所述载波平滑码更新监测器计算载波平滑码更新监测量,用以检测原始伪距测量中误差,检验统计量为:
Figure FDA0003646168390000022
其中,ρ(k)表示k时刻未平滑的伪距,ρs(k-1)代表(k-1)时刻的平滑后的伪距,
Figure FDA0003646168390000023
分别代表k、k-1时刻的载波相位测量值。
5.根据权利要求1或4所述的精密进近下地基增强系统完好性评估方法,其特征在于:步骤4中的组合漏检概率Pcom为:
Pcom=min[Pmd,CP×Pmd,CCD,Pmd,CP×Pmd,DSIGMA,Pmd,CP×Pmd,INNO],
其中Pmd,CP,Pmd,CCD,Pmd,DSIGMA,Pmd,INNO分别表示相关峰监测器的漏检概率、码/载波分歧监视器的漏检概率、双平滑电离层梯度监视器的漏检概率、载波平滑码更新监测器的漏检概率。
6.根据权利要求5所述的精密进近下地基增强系统完好性评估方法,其特征在于:所述危险误导信息概率PHMI由组合漏检概率Pcom加权得到,具体为
Figure FDA0003646168390000031
其中N为仿真参数样本的总数,i为仿真参数样本的部分测距误差超限子集。
7.根据权利要求1所述的一种精密进近下地基增强系统完好性评估方法的一种精密进近下地基增强系统完好性评估系统,其特征在于:包括参数仿真模块和测量模块,所述参数仿真模块对飞机进行精密进近参数仿真采样,测量模块整合现场采集的辅助数据,提前测量辅助数据并将其加入到仿真过程中以提升仿真系统的适用性与可靠性;基于测量模块测量的数据进行精密进近参数仿真采样,并将仿真采样结果发送到差分处理模块和监视器模块;所述差分处理模块对飞机精密进近过程差分定位仿真,所述监视器模块在仿真过程中采集每个监视器的测试参数,将每个监视器在进近期间漏检概率的最小值发送至数据记录模块进行记录,最终将采集到的每个监视器的漏检概率进行组合,以计算在极限情况下飞机精密进近过程中的GBAS的完好性风险。
CN202210226260.5A 2022-03-09 2022-03-09 一种精密进近下地基增强系统完好性评估方法及系统 Active CN114296117B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210226260.5A CN114296117B (zh) 2022-03-09 2022-03-09 一种精密进近下地基增强系统完好性评估方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210226260.5A CN114296117B (zh) 2022-03-09 2022-03-09 一种精密进近下地基增强系统完好性评估方法及系统

Publications (2)

Publication Number Publication Date
CN114296117A CN114296117A (zh) 2022-04-08
CN114296117B true CN114296117B (zh) 2022-06-21

Family

ID=80978451

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210226260.5A Active CN114296117B (zh) 2022-03-09 2022-03-09 一种精密进近下地基增强系统完好性评估方法及系统

Country Status (1)

Country Link
CN (1) CN114296117B (zh)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102103210A (zh) * 2009-12-17 2011-06-22 中国石油大学(北京) 一种卫星导航系统性能评估系统
CN104849728B (zh) * 2015-05-12 2017-07-25 北京航空航天大学 地基增强系统的完好性评估方法
CN110988929B (zh) * 2019-12-21 2020-09-22 中国电子科技集团公司第二十研究所 电离层影响下的gbas系统性能评估方法及装置
CN112130177B (zh) * 2020-09-08 2023-09-29 南京航空航天大学 一种基于稳定分布的地基增强系统完好性监测方法
CN112198533B (zh) * 2020-10-19 2022-11-08 中国电子科技集团公司第二十八研究所 一种多假设下的地基增强系统完好性评估系统及方法

Also Published As

Publication number Publication date
CN114296117A (zh) 2022-04-08

Similar Documents

Publication Publication Date Title
Datta-Barua et al. Ionospheric threat parameterization for local area global-positioning-system-based aircraft landing systems
RU2501039C2 (ru) Устройство и способ контроля целостности в реальном времени спутниковой навигационной системы
US10473791B1 (en) GBAS integrity performance evaluation method based on pseudorange error distribution model
Jung et al. Long-term ionospheric anomaly monitoring for ground based augmentation systems
Datta-Barua et al. Using WAAS ionospheric data to estimate LAAS short baseline gradients
KR101925624B1 (ko) 지역전리층지도 생성 장치 및 방법
Pesonen A framework for Bayesian receiver autonomous integrity monitoring in urban navigation
EP2037291A1 (en) Integrity monitoring method for GNSS navigation based on historical information
CN110988929B (zh) 电离层影响下的gbas系统性能评估方法及装置
EP3570073A1 (en) Multi frequency monitor for detecting ionospheric and tropospheric disturbances
Jin et al. Ionospheric correlation analysis and spatial threat model for SBAS in China region
Caamano et al. Network‐based ionospheric gradient monitoring to support GBAS
Caamano et al. Performance Evaluation of the Ionospheric Threat Mitigation Strategies in Dual-Frequency Multi-Constellation GBAS
CN111007541B (zh) 一种卫星导航地基增强系统仿真性能评估方法
Marini-Pereira et al. Advanced warning of threatening equatorial plasma bubbles to support GBAS in low latitudes
Pervan et al. Integrity monitoring for precision approach using kinematic GPS and a ground‐based pseudolite
CN114296117B (zh) 一种精密进近下地基增强系统完好性评估方法及系统
Osechas et al. Carrier‐Phase Acceleration RAIM for GNSS Satellite Clock Fault Detection
Lee et al. Results from automated ionospheric data analysis for Ground-Based Augmentation Systems (GBAS)
Yoon et al. Multi-dimensional verification methodology of ionospheric gradient observation during plasma bubble events in the Brazilian region
Lee LAAS position domain monitor analysis and test results for CAT II/III operations
CN111175789B (zh) 地基增强系统的电离层异常监测方法、装置以及系统
Hu et al. Multiple reference consistency check algorithm in GBAS based on S-values auxiliary
Konno Design of an aircraft landing system using dual-frequency GNSS
Caamano et al. Using a wide area receiver network to support GBAS ionospheric monitoring

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant