CN109901204B - 一种基于伪距误差分布模型的gbas完好性性能评估方法 - Google Patents
一种基于伪距误差分布模型的gbas完好性性能评估方法 Download PDFInfo
- Publication number
- CN109901204B CN109901204B CN201910236621.2A CN201910236621A CN109901204B CN 109901204 B CN109901204 B CN 109901204B CN 201910236621 A CN201910236621 A CN 201910236621A CN 109901204 B CN109901204 B CN 109901204B
- Authority
- CN
- China
- Prior art keywords
- error
- pseudo
- fault
- gbas
- integrity
- 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
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 32
- 101000972822 Homo sapiens Protein NipSnap homolog 2 Proteins 0.000 title claims abstract 14
- 102100022564 Protein NipSnap homolog 2 Human genes 0.000 title claims abstract 14
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000000926 separation method Methods 0.000 claims description 14
- 230000001133 acceleration Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 5
- 238000005315 distribution function Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 5
- 239000011159 matrix material Substances 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000003416 augmentation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/08—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing integrity information, e.g. health of satellites or quality of ephemeris data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
- G01S19/15—Aircraft landing systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/20—Integrity monitoring, fault detection or fault isolation of space segment
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/40—Correcting position, velocity or attitude
- G01S19/41—Differential correction, e.g. DGPS [differential GPS]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Aviation & Aerospace Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种基于伪距误差分布模型的GBAS完好性性能评估方法包括:机载接收端,通过采集伪距误差样本数据进行GBAS完好性性能评估,包括如下方法步骤:a)对伪距误差样本数据进行分组;b)针对每一组伪距误差样本分别建立高斯核二次高斯多项式尾分布模型,c)将每组伪距误差分布模型进行加权求和,得到总体伪距误差分布模型;d)伪距域误差投影到定位域误差;e)计算定位域误差大于告警限的概率,得到完好性风险概率值;f)评估GBAS系统完好性性能。本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法提高了完好性风险值的计算精确度,改善了完好性评估的有效性。
Description
技术领域
本发明涉及卫星导航技术领域,特别涉及一种基于伪距误差分布模型的GBAS完好性性能评估方法。
背景技术
地基增强系统(Ground Based Augmentation,GBAS)通过建立位置已知的地面参考站,对测距值进行差分处理实现了区域范围内定位精度的增强和导航信号完好性的实时监测。评估所安装的GBAS系统在局部空域内(以GBAS地面站为中心的50KM范围内)提供的导航服务是否满足相应的飞行性能需求需要从精度、连续性、完好性和可用性四个方面进行考察。其中完好性对于保障用户安全至关重要,是民用航空导航领域最为关注的性能需求。
完好性指机载用户利用导航系统进行相应操作的过程中,当导航系统无法用于某些预定工作时及时向用户发出告警的能力。完好性的核心指标包括告警时间、告警限和完好性风险概率值。其中完好性风险概率值指的是系统出现故障无法用于导航且无法在规定告警时间内向用户告警的概率。
然而,现有技术中对于完好性性能的评价存在着样本强相关性的问题,在面对样本数量较少的情况下,会降低完好性风险概率值的精确性,评估的有效性较差。
因此,为了解决现有技术中的一些问题,需要一种基于伪距误差分布模型的GBAS完好性性能评估方法。
发明内容
本发明的目的是提供一种基于伪距误差分布模型的GBAS完好性性能评估方法,数据处理器采集伪距差样本数据进行GBAS完好性性能评估,包括如下方法步骤:
a)对伪距误差样本数据进行分组;
b)针对每一组伪距误差样本分别建立高斯核二次高斯多项式尾分布模型,包括:
分别计算每组误差样本的标准差σi,利用均值μi为零,标准差为σi的高斯分布对误差分布核进行拟合,概率密度分布函数为:
利用高斯核分布的二次多项式对误差分布尾部进行建模,概率密度函数为:
误差样本分布满足:
c)对每组伪距误差分布模型进行加权求和处理,得到总体伪距误差分布模型;
d)伪距域误差投影到定位域误差;
e)计算定位域误差大于告警限的概率,得到完好性风险概率值;
f)评估GBAS完好性性能。
优选地,数据处理器,以1s的时间间隔采集数据。
优选地,对伪距误差样本数据进行分组,包括如下方法步骤:
优选地,利用偏差平方和最小的准则,估计伪距误差分布尾部概率密度函数的未知参数a、b的值。
优选地,给定Ef时计算第i组的最小偏差平方和Qi(Ef),由其最小值min(Qi(Ef))确定核尾分界点Ef。
优选地,评估H0假设下GBAS完好性性能通过Prisk=P(|y|>AL)Pfree-fault计算,其中,
y为定位域误差,AL为告警限,Pfree-fault为参考接收机和测距源均无故障的概率,取值为1,Prisk为完好性风险概率值。
优选地,评估H1假设下GBAS完好性性能时,在步骤a)之前还包括利用GNSS模拟器和故障仿真程序模拟单参考接收机故障,其中,
GNSS模拟器产生卫星信号,通过射频连接到参考接收机,故障仿真程序接收来自参考接收机的测距信息,按照故障参数配置,对来自参考接收机的测距信息添加斜坡误差,模拟单参考接收机故障。
优选地,评估H1假设下GBAS完好性性能通过Prisk=P(|y|>AL)Preceiver-fault计算,其中,
y为定位域误差,AL为告警限,Preceiver-fault为只有一台参考接收机故障的概率,取值为10-5,Prisk为完好性风险概率值。
优选地,评估H2假设下GBAS完好性性能时,在步骤a)之前还包括利用GNSS模拟器和故障仿真程序仿真信号畸变、低信号功率、码载波分离、星历误差过大和超加速度五个测距源故障,
其中,GNSS模拟器产生卫星信号,通过射频连接到参考接收机,故障仿真程序接收来自参考接收机的测距信息,设置故障参数,添加斜坡误差模型模拟信号畸变、低信号功率、码载波分离和星历误差过大四个测距源故障,添加平方误差模型模拟超加速度故障。
优选地,评估H2假设下GBAS完好性性能通过Prisk=P(|y|>AL)Pfault计算,其中,
y为定位域误差,AL为告警限,Pfault为信号畸变、低信号功率、码载波分离、星历误差过大和超加速度五个测距源故障的发生概率,Prisk为完好性风险概率值。
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,利用高斯核分布的二次多项式对误差尾部建模,以偏差平方和最小为准则选取核尾分界点,将核尾分开考虑的误差分布模型既可以确保核部分误差分布模型的准确性,更接近真实分布,又可以降低尾部分误差分布模型的统计不确定性。
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,考虑到GBAS系统伪距误差具有厚尾性和相关性的特征,利用混合高斯核二次高斯多项式尾分布建立伪距误差真实分布模型,结合定位域误差和伪距域误差的投影关系,计算出完好性风险概率值,评估系统完好性,提高了核部分误差分布的精确度,改善了尾部分误差分布包络紧度。误差分布模型更接近误差真实分布,完好性风险概率值更精确,同时增加了系统完好性评估的有效性。
由于GBAS系统伪距误差存在时间相关性,使得独立样本量受限,导致误差分布模型的不准确性,为解决该问题,本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,按一定时间间隔将误差样本分组,再选取每组中时间序号相同的误差样本重组成新的误差样本组,分别对新的误差样本组建立误差分布模型。
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,充分利用误差样本,将各组误差样本的分布模型进行加权求和处理得到混合误差分布模型,即总的误差分布模型
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,在样本受限的情况下减少时间相关性对误差分布建模的影响,降低模型参数的统计不确定性,获得更精确的误差分布模型,完好性风险概率值更准确,从而达到对GBAS系统完好性性能进行有效评估的目的。
由于误差分布核部分的样本较多而尾部分样本较少,为使核部分模型更精确,尾部分布模型相对保守,本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,将误差分布模型的核尾分开考虑,用高斯分布逼近核部分误差分布,用二次高斯多项式分布近似尾部分误差分布。
应当理解,前述大体的描述和后续详尽的描述均为示例性说明和解释,并不应当用作对本发明所要求保护内容的限制。
附图说明
参考随附的附图,本发明更多的目的、功能和优点将通过本发明实施方式的如下描述得以阐明,其中:
图1示意性示出了高斯分布、二次高斯多项式分布和指数分布比较。
图2示出了本发明一个实施例中采集伪距误差样本的框图。
图3示出了本发明一个实施例中H0假设下GBAS完好性性能评估的流程框图。
图4示出了本发明另一个实施例中H1假设下GBAS完好性性能评估的流程框图。
图5示出了本发明再一个实施例中H2假设下GBAS完好性性能评估的流程框图。
具体实施方式
通过参考示范性实施例,本发明的目的和功能以及用于实现这些目的和功能的方法将得以阐明。然而,本发明并不受限于以下所公开的示范性实施例;可以通过不同形式来对其加以实现。说明书的实质仅仅是帮助相关领域技术人员综合理解本发明的具体细节。
在下文中,将参考附图描述本发明的实施例。在附图中,相同的附图标记代表相同或类似的部件,或者相同或类似的步骤。下面通过具体的实施例对本发明的内容进行说明,下面通过具体的实施例对本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法进行说明。由于GBAS系统误差样本具有厚尾特性,真实误差分布模型尾部分介于高斯分布和指数分布之间。若利用高斯分布近似误差分布尾部则会低估系统完好性风险概率值,额外引入完好性风险事件。若利用指数分布逼近误差分布尾部则会过于保守,高估系统完好性风险概率值,降低系统可用性。为保守地估计完好性风险概率值又不会对系统的可用性造成太大的影响,本发明提出二次高斯多项式分布模型近似误差分布尾部,提高了完好性风险概率值的精确度。将高斯分布概率密度函数(probability density function,PDF)作为二次多项式的自变量即可得到二次高斯多项式分布。可将指数分布看作高斯分布的平方根,类似于自变量取值范围为[0,1]时二次多项式取值、一次式取值和平方根式取值的关系,选取适当的分布参数,可使二次高斯多项式分布PDF介于高斯分布PDF和指数分布PDF之间,如图1所示高斯分布、二次高斯多项式分布和指数分布比较,图中曲线a为指数分布,曲线b为二次高斯多项式分布,c为高斯分布。
依据一次式、二次多项式和平方根式在自变量等于0和自变量等于1处的斜率大小关系并假设二次多项式在自变量为0和自变量为1处取值为0,确定二次高斯多项式分布未知参数的约束条件。二次高斯多项式分布PDF为:其中a+b=1,fG(x)是高斯分布的PDF。
值得注意的是,本发明提出的二次高斯多项式分布仅对误差分布尾部进行恰当包络,保守地描述误差分布尾部真实情况,所以不必要求二次高斯多项式分布分布函数最大值为1。
为了使本发明的内容得以更加清晰的说明,首先对数据来源进行说明,如图2所示本发明一个实施例中采集伪距误差样本的框图,地面参考站100向布置于飞行器200上的机载接收端201提供差分数据,卫星300向布置于飞行器200上的机载接收端201提供星历信息,机载接收端根据地面参考站100提供的差分数据和来自卫星的星历信息计算得到伪距误差样本数据。
应当理解,通过差分数据和星历信息计算伪距误差样本数据可以采用本领域技术人员所熟知的方法得到,凡是能够实现伪距误差的计算的方法均可,本发明不再进行具体限定。
数据处理器400由机载接收端采集伪距误差数据。
由于接收机故障概率及其他故障情况发生概率很小,所以实际环境中采集的样本多是无故障情况下的样本,即H0假设,H0假设是指所有的参考接收机和测距源均正常工作没有异常,在该假设下,造成完好性风险事件的原因可能是信号的多径效应和接收机噪声。为了获取H1、H2情况下的误差样本数据,需利用GNSS模拟器结合故障仿真程序设置故障模型,采集故障情况下伪距误差数据。H1假设是指有且仅有一个参考接收机发生故障,在该假设下,发生的故障可能未被地面子系统立即检测出来,因此影响了广播增强信息的有效性,引起机载系统定位误差。H2假设是除了非H0、H1的所有情况,本发明主要考虑与测距源相关的故障情况:信号畸变、低信号功率、码载波分离、超加速度和星历误差过大五个测距源。
H0假设下GBAS完好性性能评估
如图3所示本发明一个实施例中H0假设下GBAS完好性性能评估的流程框图,利用采集的伪距误差样本数据进行GBAS完好性性能评估,包括如下方法步骤:
步骤S101、对伪距误差样本数据进行分组。
由于采集的伪距误差样本数据有限,样本之间具有强相关性且完好性风险是小概率事件,从统计角度计算完好性风险值不易实现,所以本发明通过等时间间隔采样,将样本划分为若干组,然后分别在伪距域建立误差分布模型,再将各分组误差分布模型混合得到混合高斯核二次高斯多项式尾误差分布模型,最后利用伪距域和定位域的投影关系,外推出定位域误差大于告警限的概率。
根据本发明的实施例,假设有m颗可见卫星,以1s的时间间隔在机载端采集数据。对采集的伪距误差样本数据进行分组包括如下方法步骤:
步骤S102、针对每一组伪距误差样本分别建立高斯核二次高斯多项式尾分布模型。
针对每一组伪距误差样本分别建立高斯核二次高斯多项式尾分布模型,包括:
分别计算每组误差样本的标准差σi,利用均值μi为零,标准差为σi的高斯分布对误差分布核进行拟合,概率密度分布函数为:
利用高斯核分布的二次多项式对误差分布尾部进行建模,概率密度函数为:
具体地,通过如下方法利用偏差平方和最小,估计未知参数二次项系数a和一次项系数b的值:
其中,Qi(Ef)为给定Ef时第i组的最小偏差平方和,pj为观测误差xj对应的样本频率,nf为第i组误差样本尾部样本量。
给定Ef时计算第i组的最小偏差平方和Qi(Ef),由其最小值min(Qi(Ef))确定核尾分界点Ef,并将最小偏差平方和的最小值对应的分布模型作为伪距误差分布尾部模型。
故误差样本分布满足:
步骤S103、对每组伪距误差分布模型进行加权求和处理,得到总体伪距误差分布模型。
由于单个分组样本量过少,为保证所建立的误差分布模型更接近真实误差分布,将所有样本汇合,对每个伪距误差分布模型进行加权求和处理,得到总体伪距误差分布模型。
因为相邻样本组之间所包含的信息量相近,故利用相近的样本组建立的误差分布模型相似。为了描述各个误差分布模型对总误差分布模型的贡献,为每个误差分布模型定义权重因子。将所有误差分布模型划分为四组,每组中的第一个误差分布模型权重为即
故对应第k个可见星,总体误差分布模型为:
其中,k=1,..,m。
步骤S104、伪距域误差投影到定位域误差。
第i颗卫星和接收机之间的伪距rk(k=1,..,m)为:rk=dk+ctr+xk,其中dk是第k颗卫星和接收机之间的真实距离,tr是接收机钟差,c是光速,xk是伪距误差。
将伪距域误差投影到定位域误差,则:
Y=AX,其中,
Y=(Δx,Δy,Δz,ctr)T,Δx,Δy,Δz分别是x,y,z方向的定位误差。
X=(x1,...,xm)T表示m颗卫星对应的伪距误差向量。
A=(HTW-1H)-1HTW-1,H是m×4维的系数矩阵,W是m×m维的伪距误差协方差矩阵。
故定位域误差可以表示为伪距误差X=(x1,...,xm)T的函数:
垂直定位误差yL=azX,
其中,az是A矩阵的第三行向量。
步骤S105、计算定位域误差大于告警限的概率。
水平定位误差大于告警限的概率:
垂直定位误差大于告警限的概率:
其中,yL=Δz为垂直定位误差。
定位误差大于告警限的概率为水平定位误差超限概率与垂直定位误差超限概率之和:
P(|y|>AL)=P(|yH|>HAL)+P(|yL|>LAL)。
步骤S106、评估GBAS完好性性能。
根据本发明,本实施例中,在H0假设下,参考接收机和测距源均无故障。评估GBAS完好性性能通过Prisk=P(|y|>AL)Pfree-fault计算,其中,
y为定位域误差,AL为告警限,Pfree-fault为参考接收机和测距源均无故障的概率,取值为1,Prisk为完好性风险概率值。
H1假设下GBAS完好性性能评估
图4示出了本发明另一个实施例中H1假设下GBAS完好性性能评估的流程框图,利用采集的伪距误差样本数据进行GBAS完好性性能评估,包括如下方法步骤:
步骤S201、利用GNSS模拟器和故障仿真程序模拟单参考接收机故障,其中,
GNSS模拟器产生卫星信号,通过射频连接到参考接收机,故障仿真程序接收来自参考接收机的测距信息,按照故障参数配置,对来自参考接收机的测距信息添加斜坡误差,模拟单参考接收机故障。
斜坡误差模型描述为:
ρ(t)=k(t-t0)U(t-t0),
其中ρ(t)是单参考接收机故障导致的误差,k是斜坡误差斜率,t0是故障加入时刻,U(t-t0)是阶跃函数。
将通过上述方法得到的单参考接收机故障导致的伪距误差样本数据,重复实施例H0假设下GBAS完好性性能评估步骤S101至步骤S105的过程。
步骤S202、评估GBAS完好性性能。
根据本发明,本实施例中,在H1假设下,只有一台参考接收机故障的概率。评估GBAS完好性性能通过Prisk=P(|y|>AL)Preceiver-fault计算,其中,
y为定位域误差,AL为告警限,Preceiver-fault为只有一台参考接收机故障的概率,取值为10-5,Prisk为完好性风险概率值。
H2假设下GBAS完好性性能评估
如图5所示本发明再一个实施例中H2假设下GBAS完好性性能评估的流程框图。利用采集的伪距误差样本数据进行GBAS完好性性能评估,包括如下方法步骤:
步骤S301、模拟测距源故障。
利用GNSS模拟器和故障仿真程序仿真信号畸变、低信号功率、码载波分离、星历误差过大和超加速度五个测距源故障,
其中,GNSS模拟器产生卫星信号,通过射频连接到参考接收机,故障仿真程序接收来自参考接收机的测距信息,设置故障参数,添加斜坡误差模型模拟信号畸变、低信号功率、码载波分离和星历误差过大四个测距源故障,添加平方误差模型模拟超加速度故障。添加斜坡误差模型模拟信号畸变、低信号功率、码载波分离和星历误差过大四个测距源故障与H1假设下添加斜坡误差过程相同,这里不再赘述。
平方误差模型描述为:
ρ(t)=A(t-t0)2U(t-t0),
其中,ρ(t)是超加速度导致的误差,A是二次项系数,t0是故障加入时刻,U(t-t0)是阶跃函数。
通过上述方法得到信号畸变、低信号功率、码载波分离、星历误差过大和超加速度五个测距源故障情况下伪距误差样本数据,重复实施例H0假设下GBAS完好性性能评估步骤S101至步骤S105的过程。
步骤S302、
根据本发明的实施例,在H2假设下,考虑信号畸变、低信号功率、码载波分离、星历误差过大和超加速度故障。评估GBAS完好性性能通过Prisk=P(|y|>AL)Pfault计算,其中,
y为定位域误差,AL为告警限,Pfault为信号畸变、低信号功率、码载波分离、星历误差过大和超加速度五个测距源故障的发生概率,Prisk为完好性风险概率值。
通过H0、H1、H2假设下的实测完好性风险概率值分别与规定分配的完好性风险需求值相比较,评估已安装的GBAS系统完好性性能是否达标。若实测完好性风险概率值小于风险需求值,则说明系统完好性性能良好,若实测完好性风险大于风险需求值,说明系统完好性性能不达标。
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,H0假设下,采集实测数据,分组建立高斯核二次高斯多项式尾误差分布模型。将各组样本的误差分布模型进行加权求和获得混合误差分布模型,结合定位域误差和伪距域误差投影关系,得到H0情况下完好性风险概率值。将其与H0情况所分配的完好性风险概率值比较,评估H0假设下GBAS系统完好性性能。
H1假设下,利用GNSS信号模拟器结合故障仿真程序,添加斜坡故障模拟单参考接收机故障情况,采集用户端数据,分组建立高斯核二次高斯多项式尾误差分布模型。将各组样本的误差分布模型进行加权求和获得混合误差分布模型,结合定位域误差和伪距域误差投影关系,得到H1假设下完好性风险概率值。将其与H1假设所分配的完好性风险概率值比较,评估H1假设下GBAS系统完好性性能。
H2假设下,利用GNSS信号模拟器结合故障仿真程序,添加斜坡故障模拟信号畸变、低信号功率、码载波分离、星历误差过大故障,添加平方误差故障模型模拟超加速度故障,采集用户端数据,分组建立高斯核二次高斯多项式尾误差分布模型。将各组样本的误差分布模型进行加权求和获得混合误差分布模型,结合定位域误差和伪距域误差投影关系,得到H2假设下各故障情况对应的完好性风险概率值。分别与各故障情况所分配的完好性风险值比较,评估H2假设下GBAS系统完好性性能。
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,先将误差样本按一定时间间隔进行分组,再取每组时间序号相同的样本重组成新的样本组,对每个新样本组分别建立误差分布模型。根据每个误差分布模型对总体误差模型的贡献和分布模型间信息量的关系,对每个误差分布模型赋予权重。
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,为提高伪距误差分布模型的准确度和对GBAS系统完好性性能评估的有效性,提出了二次高斯多项式分布的概念。为提高完好性风险概率值的计算精确度,利用所提出的二次高斯多项式分布近似伪距误差真实分布尾部。利用高斯核二次高斯多项式尾分布近似总体伪距误差分布,结合定位域误差和伪距域误差的投影关系,计算出完好性风险概率值。
本发明一种基于伪距误差分布模型的GBAS完好性性能评估方法,在给定误差分布核尾分界点的前提下根据偏差平方和最小的准则和未知参数的约束条件,估计伪距误差分布尾部未知参数,并求出最小偏差平方和的值。将最小偏差平方和表示为伪距误差分布模型核尾分界点的函数,求此函数的最小值来确定核尾分界点,并取最小偏差平方和的最小值对应的分布模型参数作为伪距误差分布尾部的参数值。
结合这里披露的本发明的说明和实践,本发明的其他实施例对于本领域技术人员都是易于想到和理解的。说明和实施例仅被认为是示例性的,本发明的真正范围和主旨均由权利要求所限定。
Claims (4)
1.一种基于伪距误差分布模型的GBAS完好性性能评估方法,其特征在于,数据处理器采集伪距差样本数据进行GBAS完好性性能评估,包括如下方法步骤:
a)对伪距误差样本数据进行分组;
b)针对每一组伪距误差样本分别建立高斯核二次高斯多项式尾分布模型,包括:
分别计算每组误差样本的标准差σi,利用均值μi为零,标准差为σi的高斯分布对误差分布核进行拟合,概率密度分布函数为:
利用高斯核分布的二次多项式对误差分布尾部进行建模,概率密度函数为:
误差样本分布满足:
c)对每组伪距误差分布模型进行加权求和处理,得到总体伪距误差分布模型;
d)伪距域误差投影到定位域误差;
e)计算定位域误差大于告警限的概率,得到完好性风险概率值;
f)评估GBAS完好性性能,其中,
评估H0假设下GBAS完好性性能通过Prisk=P(|y|>AL)Pfree-fault计算,其中,
y为定位域误差,AL为告警限,Pfree-fault为参考接收机和测距源均无故障的概率,取值为1,Prisk为完好性风险概率值;
评估H1假设下GBAS完好性性能时,在步骤a)之前还包括利用GNSS模拟器和故障仿真程序模拟单参考接收机故障,其中,
GNSS模拟器产生卫星信号,通过射频连接到参考接收机,故障仿真程序接收来自参考接收机的测距信息,按照故障参数配置,对来自参考接收机的测距信息添加斜坡误差,模拟单参考接收机故障;
评估H1假设下GBAS完好性性能通过Prisk=P(|y|>AL)Preceiver-fault计算,其中,
y为定位域误差,AL为告警限,Preceiver-fault为只有一台参考接收机故障的概率,取值为10-5,Prisk为完好性风险概率值;
评估H2假设下GBAS完好性性能时,在步骤a)之前还包括利用GNSS模拟器和故障仿真程序仿真信号畸变、低信号功率、码载波分离和星历误差过大和超加速度五个测距源故障,
其中,GNSS模拟器产生卫星信号,通过射频连接到参考接收机,故障仿真程序接收来自参考接收机的测距信息,设置故障参数,添加斜坡误差模型模拟信号畸变、低信号功率、码载波分离和星历误差过大四个测距源故障,添加平方误差模型模拟超加速度故障;
评估H2假设下GBAS完好性性能通过Prisk=P(|y|>AL)Pfault计算,其中,
y为定位域误差,AL为告警限,Pfault为信号畸变、低信号功率、码载波分离、星历误差过大和超加速度五个测距源故障的发生概率,Prisk为完好性风险概率值。
3.根据权利要求1所述的方法,其特征在于,利用偏差平方和最小的准则,估计伪距误差分布尾部概率密度函数的未知参数a、b的值。
4.根据权利要求1所述的方法,其特征在于,给定Ef时计算第i组的最小偏差平方和Qi(Ef),由其最小值min(Qi(Ef))确定核尾分界点Ef。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910236621.2A CN109901204B (zh) | 2019-03-27 | 2019-03-27 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
US16/430,847 US10473791B1 (en) | 2019-03-27 | 2019-06-04 | GBAS integrity performance evaluation method based on pseudorange error distribution model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910236621.2A CN109901204B (zh) | 2019-03-27 | 2019-03-27 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109901204A CN109901204A (zh) | 2019-06-18 |
CN109901204B true CN109901204B (zh) | 2020-12-04 |
Family
ID=66953550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910236621.2A Active CN109901204B (zh) | 2019-03-27 | 2019-03-27 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10473791B1 (zh) |
CN (1) | CN109901204B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110596737B (zh) * | 2019-08-24 | 2023-04-14 | 西南交通大学 | 一种gnss虚拟参考站自适应建站方法 |
CN110646820B (zh) * | 2019-09-20 | 2021-11-30 | 广州市中海达测绘仪器有限公司 | Rtk定位数据的质检方法、装置、设备和存储介质 |
US11988752B2 (en) | 2020-01-30 | 2024-05-21 | Rohde & Schwarz Gmbh & Co. Kg | Method of analyzing a ground-based augmentation system signal and test system for testing a ground-based augmentation system |
CN111366956B (zh) * | 2020-05-26 | 2020-08-14 | 天津七一二通信广播股份有限公司 | 一种应用于gbas系统的数据处理设备 |
CN111913192B (zh) * | 2020-07-08 | 2023-02-03 | 北京航空航天大学 | 一种基于关键星的gbas完好性风险分配方法 |
US11645928B2 (en) | 2020-09-28 | 2023-05-09 | Honeywell International Inc. | Method of using multiple approach guidance systems to provide higher integrity with improved performance and availability |
CN112596080B (zh) * | 2020-11-25 | 2024-04-16 | 中国人民解放军93216部队 | 无人机差分北斗起降引导系统完好性指标的测试方法 |
CN113406670A (zh) * | 2021-06-17 | 2021-09-17 | 哈尔滨工程大学 | 一种北斗非机动型广播星历故障实时监测方法 |
CN113985457B (zh) * | 2021-11-01 | 2022-07-22 | 北京航空航天大学 | 一种gbas保护级计算方法 |
CN113917495B (zh) * | 2021-12-14 | 2022-03-11 | 天津七一二通信广播股份有限公司 | 基于北斗gbas多频点多星座高可靠性自主监测方法及设备 |
CN113917510B (zh) * | 2021-12-15 | 2022-03-08 | 腾讯科技(深圳)有限公司 | 数据处理方法、装置、设备、存储介质及计算机程序产品 |
CN115941129B (zh) * | 2022-10-17 | 2024-04-16 | 西安电子科技大学 | 基于线性拟合误差矩阵的高速链路眼图实现方法 |
CN115913993B (zh) * | 2022-12-12 | 2024-06-04 | 千寻位置网络有限公司 | Gnss接收机观测数据故障建模与仿真方法及其装置 |
CN116699659B (zh) * | 2023-04-25 | 2024-03-15 | 哈尔滨工程大学 | 一种基于贝叶斯模型的连续性风险加速测试方法 |
CN116413752B (zh) * | 2023-06-08 | 2023-08-18 | 北京航空航天大学 | 基于稳定分布参数概率密度估计的gbas伪距误差包络方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545967A (zh) * | 2009-04-16 | 2009-09-30 | 北京航空航天大学 | 卫星导航完好性参数的解算方法和监测系统 |
CN101776762A (zh) * | 2009-12-30 | 2010-07-14 | 北京航空航天大学 | 基于多地基增强系统的完好性监测方法、装置与系统 |
CN102096075A (zh) * | 2010-12-23 | 2011-06-15 | 中国航天科工信息技术研究院 | 基于Galileo系统完好性概念的多模用户完好性评估方法 |
CN104200217A (zh) * | 2014-08-07 | 2014-12-10 | 哈尔滨工程大学 | 基于一种复合核函数的高光谱分类方法 |
CN104880713A (zh) * | 2015-04-10 | 2015-09-02 | 北京航空航天大学 | Gbas完好性评估处理方法和系统 |
EP2930533A1 (en) * | 2014-04-07 | 2015-10-14 | Honeywell International Inc. | Systems and methods for a code carrier divergence high-pass filter monitor |
CN105866802A (zh) * | 2016-03-24 | 2016-08-17 | 北京航空航天大学 | 一种对流层延迟完好性监测方法与装置 |
CN105911567A (zh) * | 2016-05-14 | 2016-08-31 | 四川中卫北斗科技有限公司 | 一种gnss系统完好性评估方法和装置 |
CN106324622A (zh) * | 2016-08-05 | 2017-01-11 | 西安希德电子信息技术股份有限公司 | 一种局域增强系统完好性监测及实时定位增强方法 |
CN107390237A (zh) * | 2017-06-02 | 2017-11-24 | 北京航空航天大学 | 基于无人机的北斗gbas测试方法及系统 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2928741B1 (fr) * | 2008-03-11 | 2020-06-26 | Thales | Dispositif et procede de surveillance d'integrite en temps reel d'un systeme de navigation par satellite |
ES2348778T3 (es) * | 2008-04-30 | 2010-12-14 | Gmv Aerospace And Defence S.A. | Procedimiento de determinacion autonoma de niveles de proteccion para posicionamiento gnss basado en residuos de navegacion y en una relacion de confianza isotropica. |
JP5424338B2 (ja) * | 2010-03-18 | 2014-02-26 | 日本電気株式会社 | 衛星測位システムの異常値検出装置、異常値検出方法及び異常値検出プログラム |
ES2498720T3 (es) * | 2010-12-01 | 2014-09-25 | European Space Agency | Método y aparato para determinar un parámetro de indicación de integridad que indica la integridad de la información de posicionamiento determinada en un sistema de posicionamiento global |
FR2989163B1 (fr) * | 2012-04-06 | 2014-04-11 | Thales Sa | Systeme embarque d'aide au pilotage d'un aeronef, base sur un systeme gnss, a architecture rebondante et dissimilaire pour niveau d'integrite eleve |
EP3009860B1 (en) * | 2014-10-16 | 2019-12-18 | GMV Aerospace and Defence S.A. | Method for computing an error bound of a Kalman filter based GNSS position solution |
CN107064961B (zh) * | 2017-03-24 | 2018-05-22 | 北京航空航天大学 | 对卫星导航系统完好性监测性能进行测试的方法及装置 |
CN107807371B (zh) * | 2017-09-14 | 2019-10-01 | 北京航空航天大学 | 基于北斗卫星观测数据的数据相关性处理方法 |
US11119222B2 (en) * | 2017-12-18 | 2021-09-14 | Korea Advanced Institute Of Science And Technology (Kaist) | Method and system for local-area differential GNSS for UAV navigation, and for generating optimal protection level and geometry screening therefor |
CN108828627A (zh) * | 2018-07-25 | 2018-11-16 | 中国电子科技集团公司第二十八研究所 | 一种基于高斯膨胀法的gbas完好性示警阈值估计方法 |
-
2019
- 2019-03-27 CN CN201910236621.2A patent/CN109901204B/zh active Active
- 2019-06-04 US US16/430,847 patent/US10473791B1/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545967A (zh) * | 2009-04-16 | 2009-09-30 | 北京航空航天大学 | 卫星导航完好性参数的解算方法和监测系统 |
CN101776762A (zh) * | 2009-12-30 | 2010-07-14 | 北京航空航天大学 | 基于多地基增强系统的完好性监测方法、装置与系统 |
CN102096075A (zh) * | 2010-12-23 | 2011-06-15 | 中国航天科工信息技术研究院 | 基于Galileo系统完好性概念的多模用户完好性评估方法 |
EP2930533A1 (en) * | 2014-04-07 | 2015-10-14 | Honeywell International Inc. | Systems and methods for a code carrier divergence high-pass filter monitor |
CN104200217A (zh) * | 2014-08-07 | 2014-12-10 | 哈尔滨工程大学 | 基于一种复合核函数的高光谱分类方法 |
CN104880713A (zh) * | 2015-04-10 | 2015-09-02 | 北京航空航天大学 | Gbas完好性评估处理方法和系统 |
CN105866802A (zh) * | 2016-03-24 | 2016-08-17 | 北京航空航天大学 | 一种对流层延迟完好性监测方法与装置 |
CN105911567A (zh) * | 2016-05-14 | 2016-08-31 | 四川中卫北斗科技有限公司 | 一种gnss系统完好性评估方法和装置 |
CN106324622A (zh) * | 2016-08-05 | 2017-01-11 | 西安希德电子信息技术股份有限公司 | 一种局域增强系统完好性监测及实时定位增强方法 |
CN107390237A (zh) * | 2017-06-02 | 2017-11-24 | 北京航空航天大学 | 基于无人机的北斗gbas测试方法及系统 |
Non-Patent Citations (2)
Title |
---|
GBAS基准站布设方案设计与评估方法;张悦 等;《北京航空航天大学学报》;20181231;第44卷(第12期);全文 * |
The Impact of Tropospheric Anomalies on Sea-Based JPALS Integrity;Yue Zhang et al.;《sensors》;20180807;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109901204A (zh) | 2019-06-18 |
US10473791B1 (en) | 2019-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109901204B (zh) | 一种基于伪距误差分布模型的gbas完好性性能评估方法 | |
CN104483678B (zh) | 一种空地协同的多星座卫星导航完好性多级监测方法 | |
RU2501039C2 (ru) | Устройство и способ контроля целостности в реальном времени спутниковой навигационной системы | |
JP5122066B2 (ja) | 積率限界を利用してデータ処理システムの動作のインテグリティを評価する方法及び装置 | |
CN110007317B (zh) | 一种选星优化的高级接收机自主完好性监测方法 | |
EP1965219A1 (en) | Method and system for predicting the performance of satellite navigation systems | |
CN109782303B (zh) | 一种地基增强系统的性能测试方法、装置及系统 | |
Choi et al. | Demonstrations of multi-constellation advanced RAIM for vertical guidance using GPS and GLONASS signals | |
CN110988930B (zh) | 一种先进的接收机自主完好性监测仿真评估方法及装置 | |
Lee et al. | Integrity assurance of Kalman-filter based GNSS/IMU integrated systems against IMU faults for UAV applications | |
CN110988929B (zh) | 电离层影响下的gbas系统性能评估方法及装置 | |
CN106125101A (zh) | 基于北斗的araim可用性的评估处理方法及装置 | |
Ochieng et al. | An assessment of the RAIM performance of a combined Galileo/GPS navigation system using the marginally detectable errors (MDE) algorithm | |
Rife et al. | Formulation of a time-varying maximum allowable error for ground-based augmentation systems | |
CN112198533B (zh) | 一种多假设下的地基增强系统完好性评估系统及方法 | |
CN111007541B (zh) | 一种卫星导航地基增强系统仿真性能评估方法 | |
Bang et al. | Effect of user antenna group delay variation error on advanced raim | |
Felux et al. | Towards airborne multipath models for dual constellation and dual-frequency GNSS | |
Kannemans | The Generalized Extreme Value statistical method to determine the GNSS integrity performance | |
Lee | LAAS position domain monitor analysis and test results for CAT II/III operations | |
Yoon et al. | Multi-dimensional verification methodology of ionospheric gradient observation during plasma bubble events in the Brazilian region | |
Neri et al. | GBAS NSE Model for Cat Ii/Iii Autoland Simulations | |
Simonetti et al. | GNSS Autonomous Integrity Monitoring with Barometric Pressure Measurements and Weather Data | |
Gopalakrishnan et al. | Simulation of Integrity and Availability of multi-constellation and multi-frequency GNSS augmentation systems | |
Sakai et al. | Avoiding improper modeling in SBAS ionospheric correction with shrunk observations |
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 |