CN104880713A - Gbas完好性评估处理方法和*** - Google Patents
Gbas完好性评估处理方法和*** Download PDFInfo
- Publication number
- CN104880713A CN104880713A CN201510169712.0A CN201510169712A CN104880713A CN 104880713 A CN104880713 A CN 104880713A CN 201510169712 A CN201510169712 A CN 201510169712A CN 104880713 A CN104880713 A CN 104880713A
- Authority
- CN
- China
- Prior art keywords
- subsystem
- unknown
- unknown subsystem
- time
- receiver
- 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.)
- Granted
Links
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/13—Receivers
- G01S19/20—Integrity monitoring, fault detection or fault isolation of space segment
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Security & Cryptography (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种GBAS完好性评估处理方法和***,包括:获取GBAS的时间相关误差,并根据时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H;分别对第一未知子***G1和第二未知子***G2进行时序模型拟合,获取第一未知子***G1的模型参数和拟合残差,以及第二未知子***G2的模型参数和拟合残差;根据第一未知子***G1的模型参数和拟合残差,第二未知子***G2的模型参数和拟合残差,以及已知子***H的单位脉冲响应,获取时间相关误差的包络;根据时间相关误差的包络,以获取GBAS的保护级;比较保护级和预设的告警限,并在保护级超过预设的告警线时,生成告警信号。
Description
技术领域
本发明涉及卫星导航技术,尤其涉及一种GBAS完好性评估处理方法和***。
背景技术
目前,卫星导航地基增强***(Ground Based Augmentation System,简称GBAS)广泛使用差分定位技术来提高定位精度,同时,利用完好性监测算法,对***提供误差校正信息的正确性进行度量和监测。GBAS通过计算保护级(Protection Level,简称PL)并和规定的告警限(Alarm Limit,简称AL)比较,在超过告警限时告警,从而实现对完好性的监测。
在计算保护级的过程中,始终假设误差是服从高斯分布的,然而实际误差并非服从高斯分布,且没有足够的样本确定误差的真实分布。因此,一种称为“包络”的方法被引入到计算保护级的过程中。
要对GBAS伪距校正误差进行包络,首先面临的问题是误差样本的时间相关性。传统包络方法基于独立样本,但由于已知和未知的原因,误差样本存在强烈的时间相关性。传统方法将误差样本采样时间间隔从Hatch滤波的0.5秒增大到200秒,认为这样采集的误差样本独立。但同时却带来了一个新的问题:过大的采样间隔导致误差样本匮乏。因此,在有限的样本采集时间内无法兼顾包络计算的高置信度和***的高可用性。
发明内容
本发明提供一种GBAS完好性评估处理方法和***,通过采样效率的提升,从根本上解决了GBAS的高置信度和***高可用性兼顾的问题。
本发明提供一种GBAS完好性评估处理方法,包括:
获取GBAS的时间相关误差,并根据所述时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H;
分别对所述第一未知子***G1和所述第二未知子***G2进行时序模型拟合,获取所述第一未知子***G1的模型参数和拟合残差,以及所述第二未知子***G2的模型参数和拟合残差;
根据所述第一未知子***G1的模型参数和拟合残差,所述第二未知子***G2的模型参数和拟合残差,以及所述已知子***H的单位脉冲响应,获取所述时间相关误差的包络;
根据所述时间相关误差的包络,获取所述时间相关误差的置信区间,并通过所述置信区间投影到预设定位域的目标置信区间,以获取所述GBAS的保护级;
比较所述保护级和预设的告警限,并在所述保护级超过所述预设的告警线时,生成告警信号。
本发明还提供一种GBAS完好性评估处理***,包括:
分解模块,用于获取GBAS的时间相关误差,并根据所述时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H;
拟合模块,用于分别对所述第一未知子***G1和所述第二未知子***G2进行时序模型拟合,获取所述第一未知子***G1的模型参数和拟合残差,以及所述第二未知子***G2的模型参数和拟合残差;
获取模块,用于根据所述第一未知子***G1的模型参数和拟合残差,所述第二未知子***G2的模型参数和拟合残差,以及所述已知子***H的单位脉冲响应,获取所述时间相关误差的包络;
所述获取模块,还用于根据所述时间相关误差的包络,获取所述时间相关误差的置信区间,并通过所述置信区间投影到预设定位域的目标置信区间,以获取所述GBAS的保护级;
比较模块,用于比较所述保护级和预设的告警限,并在所述保护级超过所述预设的告警线时,生成告警信号。
本发明提供的GBAS完好性评估处理方法和***,通过对引入的线性***进行分解,获取第一未知子***G1和第二未知子***G2,以及已知子***H,然后分别对第一未知子***G1和第二未知子***G2进行时序模型拟合,获取第一未知子***G1的模型参数和拟合残差,以及第二未知子***G2的模型参数和拟合残差,从而进一步计算所述引入的线性***的时间相关误差的包络,不再需要增大样本采样间隔,从而使得样本采样频率大幅提升,兼顾了包络计算的高置信度和GBAS的高可用性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明第一实施例的GBAS完好性评估处理方法的流程图;
图2为本发明第二实施例的GBAS完好性评估处理方法的流程图;
图3为本发明第二实施例的线性***分解示意图;
图4为本发明第三实施例的GBAS完好性评估处理方法的流程图;
图5为本发明第三实施例的Panit-Wu方法的流程图;
图6为本发明第四实施例的GBAS完好性评估处理方法的流程图;
图7为本发明第五实施例的GBAS完好性评估处理***的流程图;
图8为本发明第六实施例的GBAS完好性评估处理***的流程图;
图9为本发明第七实施例的GBAS完好性评估处理***的流程图;
图10为本发明第八实施例的GBAS完好性评估处理***的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明第一实施例的GBAS完好性评估处理方法的流程图;如图1所示,本实施例的GBAS完好性评估处理方法,具体可以包括如下步骤:
步骤S100、获取GBAS的时间相关误差,并根据时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H。
在本实施例中,由于考虑到GBAS中的Hatch滤波器引入的相关性以及Hatch滤波器的输入也是相关的且原因未知。因此,本发明将整个引入相关性的线性***分成两个串联的子***进行分析,其中第一未知子***G1和第二未知子***G2对应未知原因的部分,其***函数未知,需要进行拟合求得;其中,已知子***H是Hatch滤波器的***函数。
具体地,τ为Hatch滤波器的时间常数;T为采样时间间隔;且Hatch滤波器的***函数H(z)满足:
步骤S200、分别对第一未知子***G1和第二未知子***G2进行时序模型拟合,获取第一未知子***G1的模型参数和拟合残差,以及第二未知子***G2的模型参数和拟合残差。
步骤S300、根据第一未知子***G1的模型参数和拟合残差,第二未知子***G2的模型参数和拟合残差,以及已知子***H的单位脉冲响应,获取时间相关误差的包络。
步骤S400、根据时间相关误差的包络,获取时间相关误差的置信区间,并通过置信区间投影到预设定位域的目标置信区间,以获取GBAS的保护级。
具体地,在本实施例中,目标置信区间的上限作为GBAS的保护级。
步骤S500、比较保护级和预设的告警限,并在保护级超过预设的告警线时,生成告警信号。
本发明提供的GBAS完好性评估处理方法,通过获取GBAS的时间相关误差,并根据时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H;分别对第一未知子***G1和第二未知子***G2进行时序模型拟合,获取第一未知子***G1的模型参数和拟合残差,以及第二未知子***G2的模型参数和拟合残差;根据第一未知子***G1的模型参数和拟合残差,第二未知子***G2的模型参数和拟合残差,以及已知子***H的单位脉冲响应,获取时间相关误差的包络;根据时间相关误差的包络,获取时间相关误差的置信区间,并通过置信区间投影到预设定位域的目标置信区间,以获取GBAS的保护级,从而实现GBAS完好性监测。其中,本实施例根据时间相关误差,对获取的线性***进行分解成为第一未知子***G1和第二未知子***G2,以及已知子***H,通过对时间相关误差计算公式的等价变换,将第一未知子***和第二未知***的输入及***函数由未知变为已知,然后与已知的子***根据球对称包络理论得到时间相关误差包络。本发明可以直接利用时间相关样本进行分析,而不再需要增大样本采样间隔,获得时间无关样本,从而使的样本采样的频率大幅提升,兼顾了包络计算的高置信度和GBAS的高可用性。
图2为本发明第二实施例的GBAS完好性评估处理方法的流程图,在本发明第一实施例的技术方案的基础上,如图2所示,步骤S100,具体包括以下步骤:
步骤S101、根据第一接收机t时刻输出的伪距校正值PRCt 1,以及第二接收机t时刻输出的伪距校正值PRCt 2,采用公式
在本实施例中,第一接收机和第二接收机位于GBAS的地面检测站中,一般来说,GBAS的地面监测站由3~4台接收机组成。GBAS中每个接收机无法获得伪距校正误差值,因此只能通过第一接收机与第二接收机的伪距校正值之差的一半作为第一接收机或者第二接收机的伪距校正误差。
可以理解的是,由公式(1)得到的时间相关误差yt也可以作为第二接收机的时间相关误差,本实施例中,为了分析方便只针对第一接收机进行分析,第二接收机与第一接收机的分析方式相同,在此不再赘述。
第一接收机和第二接收机为GBAS中的任意两个接收机,且第一接收机和第二接收机的距离小于预设的接收机距离;GBAS中的接收机的误差均服从相同分布。一般的,预设的接收机距离都比较小,也就是说,第一接收机和第二接收机位置相近。
具体地,PRCt 1和PRCt 2是根据下列公式得到的:
其中,和分别为第一接收机和第二接收机t时刻的载波平滑码伪距,和分别为第一接收机和第二接收机t时刻到卫星的真实距离,和分别为第一接收机和第二接收机t时刻的钟差的估计值,和是从第一接收机直接获取,和是从第二接收机直接获取。
步骤S102、根据方程组
通过解方程,计算得到t时刻第一未知子***G1的输出ε1,t和t时刻第二未知子***G2的输出ε2,t。
其中,τ1为第一接收机的第一滤波时间常数;τ2为第一接收机的第二滤波时间常数;T为采样时间间隔;y1,t为t时刻第一接收机在第一滤波时间常数τ1时的时间相关误差;y2,t为t时刻第一接收机在第二滤波时间常数τ2时的时间相关误差,y1,t和y2,t是根据公式(1)计算得到的。
在本实施例中,方程组(2)是经过下列步骤推导得到的:由于载波平滑码伪距是经过Hatch滤波器平滑得到的,因此,满足如下公式:
其中,ρt为第一接收机t时刻的载波码伪距,为第一接收机t-1时刻的载波平滑码伪距,为第一接收机t时刻的载波相位观测量,为第一接收机t-1时刻的载波相位观测量,τ为Hatch滤波器的时间常数,T为采样时间间隔。
由步骤S101可知,其中,为接收机t时刻的载波平滑码伪距,Rt为接收机t时刻到卫星的真实距离,为接收机t时刻的钟差的估计值,其中的接收机可以为第一接收机也可以为第二接收机。
令εc,t为接收机t时刻的钟差的估计误差,则有
其中,Ct为接收机t时刻的钟差的真实值。
将代入 得到:
令算子d表示t时刻与t-1时刻值的差,则上式可以写成
令Δ为第一接收机和第二接收机的差分,则
由于同时定义εc,t=Δ(ρt-Rt-Ct)/2和则上式可以写成:
另重新定义 则有
因此,对于时间相关误差yt,已知子***H的输入为该输入也是第一未知子***G1和第二未知子***G2的输出,其中,ε1,t=ερ,t+εc,t。
令ε1,t为t时刻第一未知子***G1的输出,ε2,t为t时刻第二未知子***G2的输出,同时使用不同的Hatch滤波时间常数τ1和τ2进行滤波,即可得到方程组(2),通过解方程组(2)将ε1,t和ε2,t分离。
根据上述分析,可知整个线性***可以分解成如图3所示的形式,图3为本发明第二实施例的线性***分解示意图,如图3所示,待计算的白噪声和分别通过第一未知子***G1和第二未知子***G2得到ε1,t和ε2,t,ε2,t乘以(τ-T)/τ后与ε1,t融合得到已知子***H的输入,在经过已知子***H的平滑,最终得到时间相关误差yt。
在本实施例中,根据造成误差的不同原因,将线性***分解成为两个未知子***和已知子***,通过采用不同的Hatch滤波时间常数进行获得两个未知子***的输出,使得GBAS具有更广泛的适应性,能够应对GBAS中Hatch滤波时间参数的变化的情况,大大缩短GBAS投入使用的时间,同时,本发明采用采样时间间隔小的样本,不再需要增大样本采样间隔,因此,样本的采集的效率得到大幅提升,为实现高置信度的包络提供了足够的样本量,同时兼顾了GBAS的高可用性。
图4为本发明第三实施例的GBAS完好性评估处理方法的流程图,在本发明第一实施例或者第二实施例的技术方案的基础上,如图4所示,步骤S200,具体可以包括如下步骤:
步骤S201、根据第一未知子***G1的输出ε1,t和第二未知子***G2的输出ε2,t,分析第一未知子***G1的输出ε1,t的自相关函数和偏自相关函数和第二未知子***G2的输出ε2,t的自相关函数和偏自相关函数,获取第一未知子***G1对应的第一时序模型类型,以及第二未知子***G2对应的第二时序模型类型。
在本实施例中,可选地,第一未知子***G1采用自回归滑动平均(Auto-Regressive and Moving Average,简称ARMA)模型进行拟合较为合适,而第二未知子***G2则可采用更为简单的滑动平均(Moving Average,简称MA)模型进行拟合。
具体的,ARMA模型定义如下:如果随机序列{zt,t=0,±1,±2,....}满足其中序列nt为白噪声序列,zt为随机变量,t为时间参数则称zt满足ARMA(p,q)模型。
其中,(p,q)为模型的阶数,β=(φ1,φ2,...,φp,θ1,θ2,...,θq)和为模型的参数,为序列nt的标准差。特别的,当p=0时,ARMR模型退化为MR模型。
步骤S202、采用基于似然比检验的Panit-Wu方法,分别估计第一时序模型类型对应的第一参数和第二时序模型类型对应的第二参数。
在本实施例中,由于本发明处理的样本的数据量大,实时性高,需要利用计算机自动进行估计,因此,采用无需认为参与且计算量较小的基于似然比检验的Panit-Wu方法对时序模型的阶数和参数进行估计。图5为本发明第三实施例的Panit-Wu方法的流程图,如图5所示,Panit-Wu方法具体包括以下步骤:
S1:开始。
S2:n=1。
S3:拟合ARMA(2n,2n-1)。
S4:拟合ARMA(2n+2,2n+1)。
S5:将ARMA(2n+2,2n+1)拟合的残差和ARMA(2n,2n-1)拟合的残差作F检验,判断F值是否显著,若F值显著,则执行S6;若F值不显著,则执行S7。
其中,“拟合”指的是在给定模型阶数条件下利用最大似然估计方法估计模型参数,F检验,是一种似然比检验方法,检验统计量F服从F(s,N-r)分布
其中,Q0和Q1分别为嵌套的模型ARMA(p,q)和ARMA(p*,q*)的拟合残差平方,N为样本容量,r=p+q,s=r-(p*+q*),其中p*≤p,q*≤q。
S6:n=n+1;返回继续执行S3。
S7:检验φ2n和θ2n-1是否显著为零;若为否,则执行S8;若为是,则继续执行S9。
S8:确定模型ARMA(2n,2n-1)。
S9:拟合ARMA(2n-1,2n-2)。
S10:将ARMA(2n-1,2n-2)拟合的残差和S3中ARMA(2n,2n-1)拟合的残差作F检验,判断F值是否显著;若F值不显著,则执行S11;若F值显著,则执行S12。
S11:舍去ARMA(2n-1,2n-2)模型中较小的MA参数,确定模型ARMA(2n-1,m)。
其中,m<2n-2。
S12:舍去ARMA(2n,2n-1)模型中较小的MA参数,确定模型ARMA(2n,m)。
其中,m<2n-2。
若是MA模型,则从n从1开始逐步增加,利用F检验比较MA(n)和MA(n+1)。当随着n的增加,F检验不再显著,则停止增加n选择当前的MA(n)模型。
步骤S203、根据第一参数和第一未知子***G1的输出ε1,t、第二参数和第二未知子***G2的输出ε2,t,分别获取第一未知子***G1对应的第一拟合残差α1,t和第二未知子***G2对应的第二拟合残差α2,t。
在本实施例中,第一拟合残差和第二拟合残差实际上是模型的白噪声输入。具体地,根据第一参数和第一未知子***,计算对应的第一拟合残差α1,t,和根据第二参数和第二未知子***,计算对应的第二拟合残差α2,t都是发展成熟的技术,这里不再赘述。
图6为本发明第四实施例的的GBAS完好性评估处理方法的流程图;在本发明第一实施例、第二实施例和第三实施例的技术方案的基础上,如图6所示,步骤S300,具体可以包括如下步骤:
步骤S301、根据第一拟合残差α1,t和第二拟合残差α2,t,分别获取第一未知子***G1的拟合残差α1,t对应的第一包络分布密度函数和第二未知子***G2的拟合残差α2,t对应的第二包络分布密度函数
在本实施例中,由于第一拟合残差α1,t和第二拟合残差α2,t是***的输入,且样本独立,因此,可以利用传统的误差包络方法对第一拟合残差α1,t和第二拟合残差α2,t进行包络,分别获取第一包络分布密度函数和第二包络分布密度函数具体流程如下:
(1)、建立时间相关误差分布模型(Pseudorange Error DistributionModel,简称PEDM)。
具体的,PEDM是对误差分布模型的假设,通常认为时间相关误差服从一个零均值高斯核拉普拉斯尾的分布。
(2)、建立时间相关误差包络模型(Pseudorange Error OverboundModel,简称PEOM)。
具体的,PEOM是一个高斯分布,其尾部累积概率大于PEDM。
根据上述流程可知,获取的第一包络分布密度函数和第二包络分布密度函数是独立同分布高斯的形式。
本实施例中,建立时间相关误差分布和建立时间相关误差包络模型都是已发展成熟的技术,这里不再赘述。
步骤S302、根据第一包络分布密度函数和第一未知子***G1的单位脉冲响应,采用公式
计算第一未知子***G1的输出ε1,t的包络分布密度函数fε1(t)。
其中,g1为第一未知子***G1的单位脉冲响应;为第二未知子***G1的输出ε1,t的包络;“>”表示包络关系。
步骤S303、根据第二包络分布密度函数和第二未知子***G2的单位脉冲响应,采用公式
计算第二未知子***G2的输出ε2,t的包络分布密度函数
其中,g2为第二未知子***G2的单位脉冲响应;为第二未知子***G2的输出ε2,t的包络;“>”表示包络关系。
在本实施例中,对于步骤S302和步骤S303,和是根据球对称包络理论得到的。球对称理论是发展成熟的技术,这里不再赘述。
步骤S304、根据第一未知子***G1的输出ε1,t的包络分布密度函数第二未知子***G2的输出ε2,t的包络分布密度函数和已知子***H的单位脉冲响应,采用公式
计算时间相关误差yt的包络分布密度函数
其中,τ为第一接收机的滤波时间常数,T为采样时间间隔,h为已知子***H的单位脉冲响应,可由滤波时间常数τ和采样时间间隔T得到,为时间相关误差的包络。
具体的,公式(4)是根据如下推导得到的,由步骤S102可知,已知子***H的输入εt满足:
而包络关系在线性组合的过程中保持稳定,即有
包络已知子***H的输入εt。由上可知,类似公式(3)和公式(4),时间相关误差yt的包括如公式(5)所示。
在本实施例中,根据造成误差的不同原因,将线性***分解成为两个未知子***和已知子***,通过采用不同的Hatch滤波时间常数进行获得两个未知子***的输出,使得GBAS具有更广泛的适应性,能够应对GBAS中Hatch滤波时间参数的变化的情况,大大缩短GBAS投入使用的时间,同时,本发明采用采样时间间隔小的样本,不再需要增大样本采样间隔,因此,样本的采集的效率得到大幅提升,为实现高置信度的包络提供了足够的样本量,同时兼顾了GBAS的高可用性。
图7为本发明第五实施例的的GBAS完好性评估处理***的流程图;如图7所示,本实施例的GBAS完好性评估处理***,具体可以包括:分解模块400、拟合模块500、获取模块600和比较模块700;其中,
分解模块400,用于获取GBAS的时间相关误差,并根据时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H。
拟合模块500,用于分别对第一未知子***G1和第二未知子***G2进行时序模型拟合,获取第一未知子***G1的模型参数和拟合残差,以及第二未知子***G2的模型参数和拟合残差。
获取模块600,用于根据第一未知子***G1的模型参数和拟合残差,第二未知子***G2的模型参数和拟合残差,以及已知子***H的单位脉冲响应,获取时间相关误差的包络。
获取模块600,还用于根据时间相关误差的包络,获取时间相关误差的置信区间,并通过置信区间投影到预设定位域的目标置信区间,以获取GBAS的保护级。
比较模块700,用于比较保护级和预设的告警限,并在保护级超过预设的告警线时,生成告警信号。
上述***实施例中各模块,对应地,可用于执行图1方式实施例的技术方案,其实现原理和技术效果类似,在此不再赘述。
图8为本发明第六实施例的GBAS完好性评估处理***的流程图,在本发明第五实施例的技术方案的基础上,如图8所示,分解模块400,包括:第一获取单元401和第一计算单元402,其中,
第一获取单元401,用于根据第一接收机t时刻输出的伪距校正值PRCt 1,以及第二接收机t时刻输出的伪距校正值PRCt 2,采用公式
获取第一接收机时间相关误差yt。
第一计算单元402,用于根据方程组
通过解方程,计算得到t时刻第一未知子***G1的输出ε1,t和t时刻第二未知子***G2的输出ε2,t。
其中,τ1为第一接收机的第一滤波时间常数;τ2为第一接收机的第二滤波时间常数;T为采样时间间隔;y1,t为t时刻第一接收机在第一滤波时间常数τ1时的时间相关误差;y2,t为t时刻第一接收机在第二滤波时间常数τ2时的时间相关误差。
第一接收机和第二接收机为GBAS中的任意两个接收机,且第一接收机和第二接收机的距离小于预设的接收机距离;GBAS中的接收机的误差均服从相同分布。
上述***实施例中各模块,对应地,可用于执行图2方法实施例的技术方案,其实现原理和技术效果类似,在此不再赘述。
图9为本发明第七实施例的GBAS完好性评估处理***的流程图,在本发明第五实施例、第六实施例的技术方案的基础上,如图9所示,拟合模块500,包括分析单元501、估计单元502和第二获取单元503;其中,
分析单元501,用于根据第一未知子***G1的输出ε1,t和第二未知子***G2的输出ε2,t,分析ε1,t的自相关函数和偏自相关函数和ε2,t的自相关函数和偏自相关函数,获取第一未知子***G1对应的第一时序模型类型,以及第二未知子***G2对应的第二时序模型类型。
估计单元502,用于采用基于似然比检验的Panit-Wu方法,分别估计第一时序模型类型对应的第一参数和第二时序模型类型对应的第二参数。
第二获取单元503,用于根据第一参数和第一未知子***G1的输出ε1,t,第二参数和第二未知子***G2的输出ε2,t,分别获取第一未知子***G1对应的第一拟合残差α1,t和第二未知子***G2对应的第二拟合残差α2,t。
其中,第一时序模型类型为ARMA模型,第二时序模型类型为MA模型。
上述***实施例中各模块,对应地,可用于执行图4方法实施例的技术方案,其实现原理和技术效果类似,在此不再赘述。
图10为本发明第八实施例的GBAS完好性评估处理***的流程图,在本发明第五实施例、第六实施例和第七实施例的技术方案的基础上,如图10所示,获取模块600,包括:第三获取单元601和第二计算单元602,其中,
第三获取单元601,用于根据第一拟合残差α1,t和第二拟合残差α2,t,分别获取第一未知子***G1的拟合残差α1,t对应的第一包络分布密度函数和第二未知子***G2的拟合残差α2,t对应的第二包络分布密度函数
第二计算单元602,用于根据第一包络分布密度函数和第一未知子***G1的单位脉冲响应,采用公式
获取第一未知子***G1的输出ε1,t的包络分布密度函数
其中,g1为第一未知子***G1的单位脉冲响应;为第一未知子***G1的输出ε1,t的包络。
第二计算单元602,还用于根据第二包络分布密度函数和第二未知子***G2的单位脉冲响应,采用公式
获取第二未知子***G2的输出ε2,t的包络分布密度函数
其中,g2为第二未知子***G2的单位脉冲响应;为第二未知子***G2的输出ε2,t的包络。
第二计算单元602,用于根据第一未知子***G1的输出ε1,t的包络分布密度函数fε1(t),第二未知子***G2的输出ε2,t的包络分布密度函数和已知子***H的单位脉冲响应,采用公式
计算时间相关误差yt的包络分布密度函数
其中,h为已知子***H的单位脉冲响应;τ为第一接收机的滤波时间常数;T为采样时间间隔;为时间相关误差的包络,且满足:
上述***实施例中各模块,对应地,可用于执行图6方法实施例的技术方案,其实现原理和技术效果类似,在此不再赘述。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种GBAS完好性评估处理方法,其特征在于,包括:
获取GBAS的时间相关误差,并根据所述时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H;
分别对所述第一未知子***G1和所述第二未知子***G2进行时序模型拟合,获取所述第一未知子***G1的模型参数和拟合残差,以及所述第二未知子***G2的模型参数和拟合残差;
根据所述第一未知子***G1的模型参数和拟合残差,所述第二未知子***G2的模型参数和拟合残差,以及所述已知子***H的单位脉冲响应,获取所述时间相关误差的包络;
根据所述时间相关误差的包络,获取所述时间相关误差的置信区间,并通过所述置信区间投影到预设定位域的目标置信区间,以获取所述GBAS的保护级;
比较所述保护级和预设的告警限,并在所述保护级超过所述预设的告警线时,生成告警信号。
2.根据权利要求1所述的方法,其特征在于,所述获取GBAS的时间相关误差,并根据所述时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H,具体包括:
根据第一接收机t时刻输出的伪距校正值PRCt 1,以及第二接收机t时刻输出的伪距校正值PRCt 2,采用公式获取所述第一接收机的时间相关误差yt;
根据方程组 通过解方程,计算得到t时刻所述第一未知子***G1的输出ε1,t和t时刻所述第二未知子***G2的输出ε2,t;
其中,τ1为所述第一接收机的第一滤波时间常数;τ2为所述第一接收机的第二滤波时间常数;T为采样时间间隔;y1,t为t时刻所述第一接收机在所述第一滤波时间常数τ1时的时间相关误差;y2,t为t时刻所述第一接收机在第二滤波时间常数τ2时的时间相关误差;
所述第一接收机和所述第二接收机为所述GBAS中的任意两个接收机,且所述第一接收机和所述第二接收机的距离小于预设的接收机距离;所述GBAS中的接收机的误差均服从相同分布。
3.根据权利要求1或2所述的方法,其特征在于,所述分别对所述第一未知子***G1和所述第二未知子***G2进行时序模型拟合,获取所述第一未知子***G1的模型参数和拟合残差,以及所述第二未知子***G2的模型参数和拟合残差,具体包括:
根据所述第一未知子***G1的输出ε1,t和所述第二未知子***G2的输出ε2,t,分析ε1,t的自相关函数和偏自相关函数和ε2,t的自相关函数和偏自相关函数,获取所述第一未知子***G1对应的第一时序模型类型,以及所述第二未知子***G2对应的第二时序模型类型;
采用基于似然比检验的Panit-Wu方法,分别估计所述第一时序模型类型对应的第一参数和所述第二时序模型类型对应的第二参数;
根据所述第一参数和所述第一未知子***G1的输出ε1,t,所述第二参数和所述第二未知子***G2的输出ε2,t,分别获取所述第一未知子***G1对应的第一拟合残差α1,t和所述第二未知子***G2对应的第二拟合残差α2,t。
4.根据权利要求3所述的方法,其特征在于,所述根据所述第一未知子***G1的模型参数和拟合残差,所述第二未知子***G2的模型参数和拟合残差,以及所述已知子***H的单位脉冲响应,获取所述时间相关误差的包络,具体包括:
根据所述第一拟合残差α1,t和所述第二拟合残差α2,t,分别获取所述第一未知子***G1的拟合残差α1,t对应的第一包络分布密度函数和所述第二未知子***G2的拟合残差α2,t对应的第二包络分布密度函数
根据所述第一包络分布密度函数和所述第一未知子***G1的单位脉冲响应,采用公式计算所述第一未知子***G1的输出ε1,t的包络分布密度函数其中,g1为所述第一未知子***G1的单位脉冲响应;为所述第二未知子***G1的输出ε1,t的包络;
根据所述第二包络分布密度函数和所述第二未知子***G2的单位脉冲响应,采用公式:计算所述第二未知子***G2的输出ε2,t的包络分布密度函数其中,g2为所述第二未知子***G2的单位脉冲响应;为所述第二未知子***G2的输出ε2,t的包络;
根据所述第一未知子***G1的输出ε1,t的包络分布密度函数所述第二未知子***G2的输出ε2,t的包络分布密度函数和已知子***H的单位脉冲响应,采用公式计算所述时间相关误差yt的包络分布密度函数其中,h为已知子***H的单位脉冲响应;τ为所述第一接收机的滤波时间常数;T为采样时间间隔;y为时间相关误差的包络。
5.根据权利要求3所述的方法,其特征在于,所述第一时序模型类型为ARMA模型,所述第二时序模型类型为MA模型。
6.一种GBAS完好性评估处理***,其特征在于,包括:
分解模块,用于获取GBAS的时间相关误差,并根据所述时间相关误差,对获取的线性***进行分解,以获取第一未知子***G1和第二未知子***G2,以及已知子***H;
拟合模块,用于分别对所述第一未知子***G1和所述第二未知子***G2进行时序模型拟合,获取所述第一未知子***G1的模型参数和拟合残差,以及所述第二未知子***G2的模型参数和拟合残差;
获取模块,用于根据所述第一未知子***G1的模型参数和拟合残差,所述第二未知子***G2的模型参数和拟合残差,以及所述已知子***H的单位脉冲响应,获取所述时间相关误差的包络;
所述获取模块,还用于根据所述时间相关误差的包络,获取所述时间相关误差的置信区间,并通过所述置信区间投影到预设定位域的目标置信区间,以获取所述GBAS的保护级;
比较模块,用于比较所述保护级和预设的告警限,并在所述保护级超过所述预设的告警线时,生成告警信号。
7.根据权利要求6所述的***,其特征在于,所述分解模块,包括:
第一获取单元,用于根据第一接收机t时刻输出的伪距校正值PRCt 1,以及第二接收机t时刻输出的伪距校正值PRCt 2,采用公式获取所述第一接收机的时间相关误差yt;
第一计算单元,用于根据方程组 通过解方程,计算得到t时刻所述第一未知子***G1的输出ε1,t和t时刻所述第二未知子***G2的输出ε2,t;
其中,τ1为所述第一接收机的第一滤波时间常数;τ2为所述第一接收机的第二滤波时间常数;T为采样时间间隔;y1,t为t时刻所述第一接收机在所述第一滤波时间常数τ1时的时间相关误差;y2,t为t时刻所述第一接收机在所述第二滤波时间常数τ2时的时间相关误差;
所述第一接收机和所述第二接收机为所述GBAS中的任意两个接收机,且所述第一接收机和所述第二接收机的距离小于预设的接收机距离;所述GBAS中的接收机的误差均服从相同分布。
8.根据权利要求5或6所述的***,其特征在于,所述拟合模块,包括:
分析单元,用于根据所述第一未知子***G1的输出ε1,t和所述第二未知子***G2的输出ε2,t,分析ε1,t的自相关函数和偏自相关函数和ε2,t的自相关函数和偏自相关函数,获取所述第一未知子***G1对应的第一时序模型类型,以及所述第二未知子***G2对应的第二时序模型类型;
估计单元,用于采用基于似然比检验的Panit-Wu方法,分别估计所述第一时序模型类型对应的第一参数和所述第二时序模型类型对应的第二参数;
第二获取单元,用于根据所述第一参数和所述第一未知子***G1的输出ε1,t,所述第二参数和所述第二未知子***G2的输出ε2,t,分别获取所述第一未知子***G1对应的第一拟合残差α1,t和所述第二未知子***G2对应的第二拟合残差α2,t。
9.根据权利要求7所述的***,其特征在于,所述获取模块,还包括:
第三获取单元,用于根据所述第一拟合残差α1,t和所述第二拟合残差α2,t,分别计算获取所述第一未知子***G1的拟合残差α1,t对应的第一包络分布密度函数和所述第二未知子***G2的拟合残差α2,t对应的第二包络分布密度函数
第二计算单元,用于根据所述第一包络分布密度函数和所述第一未知子***G1的单位脉冲响应,采用公式计算所述第一未知子***G1的输出ε1,t的包络分布密度函数其中,g1为所述第一未知子***G1的单位脉冲响应;为所述第一未知子***G1的输出ε1,t的包络;
所述第二计算单元,还用于根据所述第二包络分布密度函数和所述第二未知子***G2的单位脉冲响应,采用公式:计算所述第二未知子***G2的输出ε2,t的包络分布密度函数其中,g2为所述第二未知子***G2的单位脉冲响应;为所述第二未知子***G2的输出ε2,t的包络;
所述第二计算单元,还用于根据所述第一未知子***G1的输出ε1,t的包络分布密度函数所述第二未知子***G2的输出ε2,t的包络分布密度函数和已知子***H的单位脉冲响应,采用公式计算所述时间相关误差yt的包络分布密度函数其中,h为已知子***H的单位脉冲响应;τ为第一接收机的滤波时间常数;T为采样时间间隔;为时间相关误差的包络。
10.根据权利要求8所述的***,其特征在于,所述第一时序模型类型为ARMA模型,所述第二时序模型类型为MA模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169712.0A CN104880713B (zh) | 2015-04-10 | 2015-04-10 | Gbas完好性评估处理方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510169712.0A CN104880713B (zh) | 2015-04-10 | 2015-04-10 | Gbas完好性评估处理方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104880713A true CN104880713A (zh) | 2015-09-02 |
CN104880713B CN104880713B (zh) | 2017-07-25 |
Family
ID=53948279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510169712.0A Active CN104880713B (zh) | 2015-04-10 | 2015-04-10 | Gbas完好性评估处理方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104880713B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107064961A (zh) * | 2017-03-24 | 2017-08-18 | 北京航空航天大学 | 对卫星导航***完好性监测性能进行测试的方法及装置 |
CN107390237A (zh) * | 2017-06-02 | 2017-11-24 | 北京航空航天大学 | 基于无人机的北斗gbas测试方法及*** |
CN109061684A (zh) * | 2018-08-01 | 2018-12-21 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 民航卫星导航完好性增强***误差的实时包络方法 |
CN109901204A (zh) * | 2019-03-27 | 2019-06-18 | 北京航空航天大学 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331618A (zh) * | 2014-10-29 | 2015-02-04 | 北京航空航天大学 | 地基增强***完好性评估方法和*** |
-
2015
- 2015-04-10 CN CN201510169712.0A patent/CN104880713B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331618A (zh) * | 2014-10-29 | 2015-02-04 | 北京航空航天大学 | 地基增强***完好性评估方法和*** |
Non-Patent Citations (5)
Title |
---|
MATHIEU JOERGER: "Carrier Phase Relative RAIM Algorithms and Protection Level Derivation", 《JOURNAL OF NAVIGATION》 * |
张倩倩: "AR模型异常值探测的Bayes方法及其在GNSS周跳探测中的应用", 《中国优秀硕士论文全文数据库基础科学辑》 * |
张倩倩等: "不同类型识别变量的自回归模型异常值探测的Bayers方法", 《测绘学报》 * |
朱衍波等: "民航GPS地基区域完好性监视***设计与实现", 《北京航空航天大学学报》 * |
淡志强等: "基于改进的包络模型的GBAS完好性评估方法", 《现代导航》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107064961A (zh) * | 2017-03-24 | 2017-08-18 | 北京航空航天大学 | 对卫星导航***完好性监测性能进行测试的方法及装置 |
CN107064961B (zh) * | 2017-03-24 | 2018-05-22 | 北京航空航天大学 | 对卫星导航***完好性监测性能进行测试的方法及装置 |
CN107390237A (zh) * | 2017-06-02 | 2017-11-24 | 北京航空航天大学 | 基于无人机的北斗gbas测试方法及*** |
CN109061684A (zh) * | 2018-08-01 | 2018-12-21 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 民航卫星导航完好性增强***误差的实时包络方法 |
CN109901204A (zh) * | 2019-03-27 | 2019-06-18 | 北京航空航天大学 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
CN109901204B (zh) * | 2019-03-27 | 2020-12-04 | 北京航空航天大学 | 一种基于伪距误差分布模型的gbas完好性性能评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104880713B (zh) | 2017-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104880713A (zh) | Gbas完好性评估处理方法和*** | |
US10168158B2 (en) | Navigation system and method for error correction | |
EP2930533B1 (en) | Systems for a code carrier divergence high-pass filter monitor | |
CN110687555B (zh) | 一种导航卫星原子钟弱频率跳变在轨自主快速检测方法 | |
EP2618181A1 (en) | Satellite navigation augmentation system and satellite navigation augmentation method | |
EP2897001A3 (en) | Multi-level/multi-threshold/multi-persistency GPS/GNSS atomic clock monitoring | |
CN112346120B (zh) | 一种槽波频散曲线的变异度计算方法及异常体的识别方法 | |
CN104766291A (zh) | 多台摄像机标定方法和*** | |
DE102014211175A1 (de) | Verfahren und System zur Initialisierung eines Sensorfusionssystems | |
KR101433908B1 (ko) | Gnss 데이터 품질 평가 방법 및 시스템 | |
CN104280741A (zh) | 电离层异常检测方法 | |
CN113654899B (zh) | 输送管道的压力分析方法、装置、设备及存储介质 | |
CN105115573A (zh) | 一种洪水流量预报的校正方法和装置 | |
CN105116419A (zh) | 一种gnss接收机、双路载波相位伪距平滑方法 | |
CN105387801A (zh) | 一种地铁隧道管片错台量检测方法 | |
CN105578587B (zh) | 一种时间频率比对方法和设备 | |
CN105093251A (zh) | Gnss接收机静态模式下的高精度相对定位方法 | |
CN104331618A (zh) | 地基增强***完好性评估方法和*** | |
CN106371092B (zh) | 一种基于gps与强震仪观测自适应组合的形变监测方法 | |
CN103454653B (zh) | 一种基于gps***的野值替换方法及装置 | |
EP2239598A1 (de) | Verfahren und Vorrichtung zur Abschätzung des Integritätsrisikos in einem Satellitennavigationssystem | |
CN107884791B (zh) | 一种gnss定位数据的滤波方法 | |
KR101596323B1 (ko) | 위치 추정 시스템, 방법, 및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 매체 | |
US20240255649A1 (en) | Network rtk anti-ionospheric disturbance positioning method, apparatus, system, device and storage medium | |
CN116481416B (zh) | 基于北斗导航的桥梁挠度监测方法、电子设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |