基于极值理论的卫星导航地基增强系统的误差包络方法
基于极值理论的卫星导航地基增强系统的误差包络方法
【技术领域】
[0001] 本发明涉及卫星导航技术，尤其涉及一种基于极值理论的卫星导航地基增强系统 的误差包络方法。
【背景技术】
[0002] 地基增强系统（Ground-Based Augmentation Systems，简称为：GBAS)采用差分技 术提高飞机定位精度。目前GBAS已可以满足最高到CAT I的民航运行需求，可作为主用甚 至唯一导航系统使用，从而使卫星导航系统替代传统的陆基无线电导航系统成为可能。
[0003] 为保证飞行安全，GBAS的机载设备在实时进行差分定位的同时，还要计算指定风 险概率下定位误差的上限，称为保护级。在GBAS中，地面站实时计算视界内每颗卫星的伪 距校正值，同时，地面站假设伪距校正误差服从零均值高斯分布，并估计其标准差〇 每颗卫星的伪距校正值及校正误差标准差被广播给飞机。飞机假设地面站发送的伪距校正 值的误差是零均值高斯分布的，标准差为〇 ^gnd，以此来计算保护级。
[0004] 但在实际中，地面反射多径等引起的误差可能是非高斯的、非零均值的，或者没有 足够的数据来验证实际误差是高斯分布的，导致实际误差的标准差超过了 0 Pund值，造成 潜在的完好性风险。因此，为了补偿假设的误差概率分布与真实误差概率分布之间的差，必 须找到一定的方法来处理误差的这些分布特性，以保证保护级的可靠性，并且这个方法不 需要误差必须是高斯分布、方差已知的。
[0005] 目前的GBAS中普遍使用了一种称为包络的技术来解决此问题，首先根据实际观 测值计算误差标准差的估计值〇Pnd_est，然后计算放大因子（Inflation Factor)kinf，并根 据公式：kinfX o ，使得计算出的〇卩^值能够包络（Overbound)实际误差， 从而机载接收机根据〇 Pqnd计算的保护级可以满足完好性需求。
[0006] 然而，实际中可用于计算误差包络的独立样本数量与误差包络所要保护的风险概 率相比极其有限，即使在真实误差满足零均值高斯分布的条件下，也仍然需要考虑由于使 用的样本数量有限所带来的不确定性。此外，由于真实误差来源于具有不同标准差的总体 分布、处理过程中导致的误差混合和不同参考接收机数据间的相关性等问题，导致真实误 差呈现厚尾分布，并且其真实分布未知，而现有方法对分布尾进行了保守的假设，导致所计 算出的放大因子过大，从而导致系统的连续性降低。
【发明内容】
[0007] 本发明实施例提供一种基于极值理论的卫星导航地基增强系统的误差包络方法， 以克服现有方法对分布尾进行了保守的假设，导致所计算出的放大因子过大，从而导致系 统的连续性降低的问题。
[0008] 本发明第一方面提供一种基于极值理论的卫星导航地基增强系统的误差包络方 法，包括：
[0009] 获取至少两组伪距误差校正样本值；
[0010] 对每组伪距误差校正样本值执行下述操作：
[0011] 确定本组伪距误差校正样本值对应的具有指定风险的第一分位数和第一标准 差；
[0012] 根据地基增强系统GBAS总的运行风险在标准高斯分布下所述GBAS总的运行风险 对应的第二分位数；
[0013] 根据所述第一分位数、所述第二分位数和所述第一标准差确定本组伪距误差校正 样本值对应的放大因子。
[0014] 结合第一方面，在第一方面的第一种可能的实现方式中，所述确定本组伪距误差 校正样本值对应的具有指定风险的第一分位数，包括：
[0015] 确定所述本组伪距误差校正样本值对应的阈值；
[0016] 根据重采样法获得所述本组伪距误差校正样本值对应的至少两组重采样样本 值；
[0017] 根据所述本组伪距误差校正样本值对应的阈值和所述本组伪距误差校正样本值 对应的所述至少两组重采样样本值中的每组重采样样本值确定每组重采样样本值中的超 出量样本值，所述超出量样本值为所述每组重采样样本值中大于所述阈值的值与所述阈值 的差值；
[0018] 根据所述超出量样本值确定所述第一分位数。
[0019] 结合第一方面的第一种可能的实现方式，在第一方面的第二种可能的实现方式 中，所述根据所述超出量样本值确定所述第一分位数，包括：
[0020] 根据本组重采样样本值对应的超出量样本值确定所述本组重采样样本值中超过 所述阈值的第一概率；
[0021] 确定所述本组重采样样本值对应的超出量样本值的第二标准差；
[0022] 根据所述第一概率和所述第二标准差确定所述第一分位数。
[0023] 结合第一方面的第二种可能的实现方式，在第一方面的第三种可能的实现方式 中，所述根据所述第一概率和所述第二标准差确定所述第一分位数，包括：
[0024] 根据所述第一概率确定指定置信度的第一参数置信限值；
[0025] 根据所述第二标准差确定指定置信度的第二参数置信限值；
[0026] 根据所述第一参数置信限值和所述第二参数置信限值确定所述第一分位数。
[0027] 结合第一方面的第一种可能的实现方式，在第一方面的第四种可能的实现方式 中，所述确定所述本组伪距误差校正样本值对应的阈值，包括：
[0028] 根据平均超出量函数MEF确定所述本组伪距误差校正样本值对应的阈值。
[0029] 结合第一方面，在第一方面的第五种可能的实现方式中，所述获取至少两组伪距 误差校正样本值，包括：
[0030] 根据预设时间间隔获取至少两个伪距误差校正样本值；
[0031] 对所述至少两个伪距误差校正样本值进行分组得到至少两组伪距误差校正样本 值。
[0032] 本发明中，获取至少两组伪距误差校正样本值，并确定各组伪距误差校正样本值 对应的具有指定风险的第一分位数和第一标准差，然后根据地基增强系统GBAS总的运行 风险确定在标准高斯分布下GBAS总的运行风险对应的第二分位数，进而根据第一分位数、 第二分位数和第一标准差确定放大因子，最后，根据放大因子确定本组伪距误差校正样本 值对应的伪距校正值，通过上述方法可以准确计算出放大因子，从而计算的伪距校正值能 够包络实际误差，有效提高系统的连续性。
【附图说明】
[0033] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发 明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以 根据这些附图获得其他的附图。
[0034]图1为本发明实施例提供的基于极值理论的卫星导航地基增强系统的误差包络 方法的流程图。
【具体实施方式】
[0035] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0036]图1为本发明实施例提供的基于极值理论的卫星导航地基增强系统的误差包络 方法的流程图，如图1所示，本实施例的方法可以包括：
[0037] 步骤101 :获取至少两组伪距误差校正样本值。
[0038] 具体的，根据预设时间间隔获取至少两个伪距误差校正样本值，然后，对至少两个 伪距误差校正样本值进行分组得到至少两组伪距误差校正样本值。
[0039] 在实际应用中，通常，GBAS地面参考站每间隔0.5秒计算得到1个当前时刻的伪 距差分校正值，对应可以得到一个伪距校正误差的样本值。然后获得全部时刻的至少两个 伪距校正误差的样本值，对全部时刻的伪距校正误差样本值进行筛选，以得到独立的伪距 校正误差样本值。
[0040] 典型的筛选方法为：首先按100秒的时间间隔对全部时刻的伪距校正误差样本值 进行采样，然后对采样后的伪距校正误差样本值按照其对应的仰角进行分组，即，仰角为0 度到5度的伪距校正误差样本值为一组，5度到10度的伪距校正误差样本值为一组，以此类 推，共得到18个分组的伪距校正误差样本值。
[0041] 步骤102 :对每组伪距误差校正样本值执行下述操作：也即，对每一组伪距校正误 差样本值分别执行步骤103至步骤105中的方法。
[0042] 例如：根据预设时间间隔获取至少50000个伪距误差校正样本值，然后对至少 50000个伪距误差校正样本值进行分组得到18组伪距误差校正样本值。
[0043] 步骤103 :确定各组伪距误差校正样本值对应的具有指定风险的第一分位数和第 一标准差。
[0044] 具体的，首先需要确定本组伪距误差校正样本值对应的阈值，设{XJ，i = 1，…， N，为经筛选的一个分组的伪距误差样本值，N为本组伪距误差样本值的个数，利用平均超出 量函数（Mean Excess Function，简称为：MEF)确定本组伪距误差校正样本值对应的阈值。
[0045] 具体的确定阈值的方法为：
[0046] 对{XJ中的伪距误差样本值按照从大到小的顺序进行排列，得到次序统计量 X#X2彡…彡X N，其中，MEF定义为：
[0048] 其中Nu为超过阈值u的样本数量，X/为伪距误差样本值中大于u的样本值，E(u) 代表大于阈值的样本的超出量的数学期望。
[0049] 在二维坐标轴中按照上述MEF的公式绘制点（u，E(u ))构成的曲线。通过选取充 分大的叫，使得当x彡 U(l时E(x)为近似正斜率的线性函数，取U(l为阈值。
[0050] 值得注意的是，通常选取的阈值叫在1. 96至2. 58倍样本标准差0 samp之间。
[0051] 然后，根据重采样法获得每组伪距误差校正样本值对应的至少两组重采样样本 值。
[0052] 根据GBAS伪距误差样本，采用重采样方法获得多组伪距误差的重采样样本；
[0053] 重采样方法可采用自助法（Bootstrap)等方法。
[0054] 以自助法为例，对{XJ进行B次重采样，得到B组重采样样本：
[0055] {X*bi}，i = 1，…，N ;b = 1，…，B。
[0056] 进而，根据该组伪距误差校正样本值对应的阈值和每组重采样样本值确定每组重 采样样本值中的超出量样本值，超出量样本值为每组重采样样本值中大于阈值的值与阈值 的差值，
[0057] 可根据如下公式计算超出量样本值，Y = X\_U(I，X、）>U(I，其中，Y为超出 量样本值。
[0058] 最后，根据超出量样本值确定第一分位数。
[0059] 可选的，根据超出量样本值确定第一分位数，包括：
[0060] 根据本组重采样样本值对应的超出量样本值确定本组重采样样本值中超过阈值 的第一概率；
[0061] 确定本组重采样样本值对应的超出量样本值的第二标准差；
[0062] 根据对GBAS伪距误差的现有研宄，GBAS伪距误差的概率分布的核服从高斯分布， 概率分布的尾未知，但处于高斯分布和拉普拉斯分布之间。根据极值理论，不论GBAS伪距 误差X的分布为何种形式，其极大值服从I型极值分布（Gumbel分布）。进一步的，其超出 量Y服从I型帕累托（Pareto)分布：
[0063] Fe(y) = Gjy) = 1-exp {_y/〇 J
[0064] 其中，Fe(y)为超出量Y的概率分布，〇丨为超出量Y的第二标准差。
[0065] 超出量Y的概率分布Fe (y)与伪距误差X的概率分布Fr (x)间的关系为：
[0067] 因此：根据上面公式可知：
[0068] Fr(x) = (l-Ptail)+PtailFe(x_u)，x>u
[0069] 其中 ptail= l-Fr(u)〇
[0070] 上述模型中需要估计的参数为第一概率Ptail和〇 ^
[0071] 而Ptail的参数估计值为每一个重采样样本的经验分布：
[0072] 也即，ptail=m/n
[0073] 采用极大似然法估计〇 i，其估计值为：
[0074] 〇 !=E(Y)
[0075] 即可根据估计出的第一概率Ptail和第二标准差〇 :确定第一分位数。
[0076] 可选的，根据第一概率和第二标准差确定第一分位数，包括：
[0077] 对估计出的PtaijP 〇的上限进行置信度检验。
[0078] 令根据B组重采样样本估计的参数值rtail，b}和{ <}的经验a分位数为Ptail 和〇的置信度为a的上限。
[0079] 根据参数置信限值flu(㈨和参数置信限值 岣确定放大因子kinf。
[0080] 进一步的，将估计出的第一概率Ptail对应的第一参数置信限值P#tail， B(a )和第二 标准差〇 i对应的第一参数置信限值〇 a )带入：
[0081] Fr(x) = (l-Ptail)+PtailFe(x_u)，x>u
[0082] 得到GBAS伪距误差的概率分布：
[0084] 最后根据第一参数置信限值P#tail，B(a)和第二参数置信限值〇\(a)确定第一分 位数。
[0085] 其中，根据下述公式计算第一分位数：
[0087] 其中，Qr，risk为第一分位数，u。为阈值，P' ail，B(a )为第一参数置信限值，〇 '(a ) 为第二参数置信限值，P&k为GBAS对应的指定风险。
[0088] 步骤104 :根据地基增强系统GBAS总的运行风险确定标准高斯分布对应GBAS总 的运行风险对应的第二分位数。
[0089] 步骤105 :根据第一分位数、第二分位数和第一标准差确定本组伪距误差校正样 本值对应的放大因子。
[0090] 具体的，根据第一分位数、第二分位数和第一标准差确定放大因子kinf，包括：
[0091] 根据以下公式确定放大因子kinf:
[0092] kinf - Q r, risk/kffmd/0samp
[0093] 其中，Quisk为第一分位数、kffDld为第二分位数，a SMP为第一标准差。
[0094] 本发明实施例提供一种基于极值理论的卫星导航地基增强系统的误差包络方法， 包括：获取至少两组伪距误差校正样本值，并确定各组伪距误差校正样本值对应的具有指 定风险的第一分位数和第一标准差，然后根据地基增强系统GBAS总的运行风险确定在标 准高斯分布下GBAS总的运行风险对应的第二分位数，进而根据第一分位数、第二分位数和 第一标准差确定本组伪距误差校正样本值对应的放大因子，最后，根据放大因子确定本组 伪距误差校正样本值对应的伪距校正值，通过上述方法可以准确计算出放大因子，从而计 算的伪距校正值能够包络实际误差，有效提高系统的连续性。
[0095] 本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通 过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程 序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：只读存储记忆 体（Read-Only Memory，ROM)、随机存储记忆体（Random Access Memory，RAM)、磁碟或者光 盘等各种可以存储程序代码的介质。
[0096] 最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制； 尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其 依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征 进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技 术方案的范围。
【主权项】
1. 一种基于极值理论的卫星导航地基增强系统的误差包络方法，其特征在于，包括： 获取至少两组伪距误差校正样本值； 对每组伪距误差校正样本值执行下述操作： 确定本组伪距误差校正样本值对应的具有指定风险的第一分位数和第一标准差； 根据地基增强系统GBAS总的运行风险确定在标准高斯分布下所述GBAS总的运行风险 对应的第二分位数； 根据所述第一分位数、所述第二分位数和所述第一标准差确定本组伪距误差校正样本 值对应的放大因子。2. 根据权利要求1所述的方法，其特征在于，所述确定本组伪距误差校正样本值对应 的具有指定风险的第一分位数，包括： 确定所述本组伪距误差校正样本值对应的阈值； 根据重采样法获得所述本组伪距误差校正样本值对应的至少两组重采样样本值； 根据所述本组伪距误差校正样本值对应的阈值和所述本组伪距误差校正样本值对应 的所述至少两组重采样样本值中的每组重采样样本值确定每组重采样样本值中的超出量 样本值，所述超出量样本值为所述每组重采样样本值中大于所述阈值的值与所述阈值的差 值； 根据所述超出量样本值确定所述第一分位数。3. 根据权利要求2所述的方法，其特征在于，所述根据所述超出量样本值确定所述第 一分位数，包括： 根据本组重采样样本值对应的超出量样本值确定所述本组重采样样本值中超过所述 阈值的第一概率； 确定所述本组重采样样本值对应的超出量样本值的第二标准差； 根据所述第一概率和所述第二标准差确定所述第一分位数。4. 根据权利要求3所述的方法，其特征在于，所述根据所述第一概率和所述第二标准 差确定所述第一分位数，包括： 根据所述第一概率确定指定置信度的第一参数置信限值； 根据所述第二标准差确定指定置信度的第二参数置信限值； 根据所述第一参数置信限值和所述第二参数置信限值确定所述第一分位数。5. 根据权利要求2所述的方法，其特征在于，所述确定所述本组伪距误差校正样本值 对应的阈值，包括： 根据平均超出量函数MEF确定所述本组伪距误差校正样本值对应的阈值。6. 根据权利要求1所述的方法，其特征在于，所述获取至少两组伪距误差校正样本值， 包括： 根据预设时间间隔获取至少两个伪距误差校正样本值； 对所述至少两个伪距误差校正样本值进行分组得到至少两组伪距误差校正样本值。
【专利摘要】本发明提供一种基于极值理论的卫星导航地基增强系统的误差包络方法，包括：获取至少两组伪距误差校正样本值，并确定各组伪距误差校正样本值对应的具有指定风险的第一分位数和第一标准差，然后根据地基增强系统GBAS总的运行风险确定在标准高斯分布下GBAS总的运行风险对应的第二分位数，进而根据第一分位数、第二分位数和第一标准差确定本组伪距误差校正样本值对应的放大因子，通过上述方法可以准确计算出放大因子，从而计算的伪距校正值能够包络实际误差，有效提高系统的连续性。
【IPC分类】G01S19/41
【公开号】CN104898140
【申请号】CN201510284335
【发明人】薛瑞, 张军, 朱衍波, 王志鹏
【申请人】北京航空航天大学
【公开日】2015年9月9日
【申请日】2015年5月28日
【技术领域】
[0001] 本发明涉及卫星导航技术，尤其涉及一种基于极值理论的卫星导航地基增强系统 的误差包络方法。
【背景技术】
[0002] 地基增强系统（Ground-Based Augmentation Systems，简称为：GBAS)采用差分技 术提高飞机定位精度。目前GBAS已可以满足最高到CAT I的民航运行需求，可作为主用甚 至唯一导航系统使用，从而使卫星导航系统替代传统的陆基无线电导航系统成为可能。
[0003] 为保证飞行安全，GBAS的机载设备在实时进行差分定位的同时，还要计算指定风 险概率下定位误差的上限，称为保护级。在GBAS中，地面站实时计算视界内每颗卫星的伪 距校正值，同时，地面站假设伪距校正误差服从零均值高斯分布，并估计其标准差〇 每颗卫星的伪距校正值及校正误差标准差被广播给飞机。飞机假设地面站发送的伪距校正 值的误差是零均值高斯分布的，标准差为〇 ^gnd，以此来计算保护级。
[0004] 但在实际中，地面反射多径等引起的误差可能是非高斯的、非零均值的，或者没有 足够的数据来验证实际误差是高斯分布的，导致实际误差的标准差超过了 0 Pund值，造成 潜在的完好性风险。因此，为了补偿假设的误差概率分布与真实误差概率分布之间的差，必 须找到一定的方法来处理误差的这些分布特性，以保证保护级的可靠性，并且这个方法不 需要误差必须是高斯分布、方差已知的。
[0005] 目前的GBAS中普遍使用了一种称为包络的技术来解决此问题，首先根据实际观 测值计算误差标准差的估计值〇Pnd_est，然后计算放大因子（Inflation Factor)kinf，并根 据公式：kinfX o ，使得计算出的〇卩^值能够包络（Overbound)实际误差， 从而机载接收机根据〇 Pqnd计算的保护级可以满足完好性需求。
[0006] 然而，实际中可用于计算误差包络的独立样本数量与误差包络所要保护的风险概 率相比极其有限，即使在真实误差满足零均值高斯分布的条件下，也仍然需要考虑由于使 用的样本数量有限所带来的不确定性。此外，由于真实误差来源于具有不同标准差的总体 分布、处理过程中导致的误差混合和不同参考接收机数据间的相关性等问题，导致真实误 差呈现厚尾分布，并且其真实分布未知，而现有方法对分布尾进行了保守的假设，导致所计 算出的放大因子过大，从而导致系统的连续性降低。
【发明内容】
[0007] 本发明实施例提供一种基于极值理论的卫星导航地基增强系统的误差包络方法， 以克服现有方法对分布尾进行了保守的假设，导致所计算出的放大因子过大，从而导致系 统的连续性降低的问题。
[0008] 本发明第一方面提供一种基于极值理论的卫星导航地基增强系统的误差包络方 法，包括：
[0009] 获取至少两组伪距误差校正样本值；
[0010] 对每组伪距误差校正样本值执行下述操作：
[0011] 确定本组伪距误差校正样本值对应的具有指定风险的第一分位数和第一标准 差；
[0012] 根据地基增强系统GBAS总的运行风险在标准高斯分布下所述GBAS总的运行风险 对应的第二分位数；
[0013] 根据所述第一分位数、所述第二分位数和所述第一标准差确定本组伪距误差校正 样本值对应的放大因子。
[0014] 结合第一方面，在第一方面的第一种可能的实现方式中，所述确定本组伪距误差 校正样本值对应的具有指定风险的第一分位数，包括：
[0015] 确定所述本组伪距误差校正样本值对应的阈值；
[0016] 根据重采样法获得所述本组伪距误差校正样本值对应的至少两组重采样样本 值；
[0017] 根据所述本组伪距误差校正样本值对应的阈值和所述本组伪距误差校正样本值 对应的所述至少两组重采样样本值中的每组重采样样本值确定每组重采样样本值中的超 出量样本值，所述超出量样本值为所述每组重采样样本值中大于所述阈值的值与所述阈值 的差值；
[0018] 根据所述超出量样本值确定所述第一分位数。
[0019] 结合第一方面的第一种可能的实现方式，在第一方面的第二种可能的实现方式 中，所述根据所述超出量样本值确定所述第一分位数，包括：
[0020] 根据本组重采样样本值对应的超出量样本值确定所述本组重采样样本值中超过 所述阈值的第一概率；
[0021] 确定所述本组重采样样本值对应的超出量样本值的第二标准差；
[0022] 根据所述第一概率和所述第二标准差确定所述第一分位数。
[0023] 结合第一方面的第二种可能的实现方式，在第一方面的第三种可能的实现方式 中，所述根据所述第一概率和所述第二标准差确定所述第一分位数，包括：
[0024] 根据所述第一概率确定指定置信度的第一参数置信限值；
[0025] 根据所述第二标准差确定指定置信度的第二参数置信限值；
[0026] 根据所述第一参数置信限值和所述第二参数置信限值确定所述第一分位数。
[0027] 结合第一方面的第一种可能的实现方式，在第一方面的第四种可能的实现方式 中，所述确定所述本组伪距误差校正样本值对应的阈值，包括：
[0028] 根据平均超出量函数MEF确定所述本组伪距误差校正样本值对应的阈值。
[0029] 结合第一方面，在第一方面的第五种可能的实现方式中，所述获取至少两组伪距 误差校正样本值，包括：
[0030] 根据预设时间间隔获取至少两个伪距误差校正样本值；
[0031] 对所述至少两个伪距误差校正样本值进行分组得到至少两组伪距误差校正样本 值。
[0032] 本发明中，获取至少两组伪距误差校正样本值，并确定各组伪距误差校正样本值 对应的具有指定风险的第一分位数和第一标准差，然后根据地基增强系统GBAS总的运行 风险确定在标准高斯分布下GBAS总的运行风险对应的第二分位数，进而根据第一分位数、 第二分位数和第一标准差确定放大因子，最后，根据放大因子确定本组伪距误差校正样本 值对应的伪距校正值，通过上述方法可以准确计算出放大因子，从而计算的伪距校正值能 够包络实际误差，有效提高系统的连续性。
【附图说明】
[0033] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发 明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以 根据这些附图获得其他的附图。
[0034]图1为本发明实施例提供的基于极值理论的卫星导航地基增强系统的误差包络 方法的流程图。
【具体实施方式】
[0035] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0036]图1为本发明实施例提供的基于极值理论的卫星导航地基增强系统的误差包络 方法的流程图，如图1所示，本实施例的方法可以包括：
[0037] 步骤101 :获取至少两组伪距误差校正样本值。
[0038] 具体的，根据预设时间间隔获取至少两个伪距误差校正样本值，然后，对至少两个 伪距误差校正样本值进行分组得到至少两组伪距误差校正样本值。
[0039] 在实际应用中，通常，GBAS地面参考站每间隔0.5秒计算得到1个当前时刻的伪 距差分校正值，对应可以得到一个伪距校正误差的样本值。然后获得全部时刻的至少两个 伪距校正误差的样本值，对全部时刻的伪距校正误差样本值进行筛选，以得到独立的伪距 校正误差样本值。
[0040] 典型的筛选方法为：首先按100秒的时间间隔对全部时刻的伪距校正误差样本值 进行采样，然后对采样后的伪距校正误差样本值按照其对应的仰角进行分组，即，仰角为0 度到5度的伪距校正误差样本值为一组，5度到10度的伪距校正误差样本值为一组，以此类 推，共得到18个分组的伪距校正误差样本值。
[0041] 步骤102 :对每组伪距误差校正样本值执行下述操作：也即，对每一组伪距校正误 差样本值分别执行步骤103至步骤105中的方法。
[0042] 例如：根据预设时间间隔获取至少50000个伪距误差校正样本值，然后对至少 50000个伪距误差校正样本值进行分组得到18组伪距误差校正样本值。
[0043] 步骤103 :确定各组伪距误差校正样本值对应的具有指定风险的第一分位数和第 一标准差。
[0044] 具体的，首先需要确定本组伪距误差校正样本值对应的阈值，设{XJ，i = 1，…， N，为经筛选的一个分组的伪距误差样本值，N为本组伪距误差样本值的个数，利用平均超出 量函数（Mean Excess Function，简称为：MEF)确定本组伪距误差校正样本值对应的阈值。
[0045] 具体的确定阈值的方法为：
[0046] 对{XJ中的伪距误差样本值按照从大到小的顺序进行排列，得到次序统计量 X#X2彡…彡X N，其中，MEF定义为：
[0048] 其中Nu为超过阈值u的样本数量，X/为伪距误差样本值中大于u的样本值，E(u) 代表大于阈值的样本的超出量的数学期望。
[0049] 在二维坐标轴中按照上述MEF的公式绘制点（u，E(u ))构成的曲线。通过选取充 分大的叫，使得当x彡 U(l时E(x)为近似正斜率的线性函数，取U(l为阈值。
[0050] 值得注意的是，通常选取的阈值叫在1. 96至2. 58倍样本标准差0 samp之间。
[0051] 然后，根据重采样法获得每组伪距误差校正样本值对应的至少两组重采样样本 值。
[0052] 根据GBAS伪距误差样本，采用重采样方法获得多组伪距误差的重采样样本；
[0053] 重采样方法可采用自助法（Bootstrap)等方法。
[0054] 以自助法为例，对{XJ进行B次重采样，得到B组重采样样本：
[0055] {X*bi}，i = 1，…，N ;b = 1，…，B。
[0056] 进而，根据该组伪距误差校正样本值对应的阈值和每组重采样样本值确定每组重 采样样本值中的超出量样本值，超出量样本值为每组重采样样本值中大于阈值的值与阈值 的差值，
[0057] 可根据如下公式计算超出量样本值，Y = X\_U(I，X、）>U(I，其中，Y为超出 量样本值。
[0058] 最后，根据超出量样本值确定第一分位数。
[0059] 可选的，根据超出量样本值确定第一分位数，包括：
[0060] 根据本组重采样样本值对应的超出量样本值确定本组重采样样本值中超过阈值 的第一概率；
[0061] 确定本组重采样样本值对应的超出量样本值的第二标准差；
[0062] 根据对GBAS伪距误差的现有研宄，GBAS伪距误差的概率分布的核服从高斯分布， 概率分布的尾未知，但处于高斯分布和拉普拉斯分布之间。根据极值理论，不论GBAS伪距 误差X的分布为何种形式，其极大值服从I型极值分布（Gumbel分布）。进一步的，其超出 量Y服从I型帕累托（Pareto)分布：
[0063] Fe(y) = Gjy) = 1-exp {_y/〇 J
[0064] 其中，Fe(y)为超出量Y的概率分布，〇丨为超出量Y的第二标准差。
[0065] 超出量Y的概率分布Fe (y)与伪距误差X的概率分布Fr (x)间的关系为：
[0067] 因此：根据上面公式可知：
[0068] Fr(x) = (l-Ptail)+PtailFe(x_u)，x>u
[0069] 其中 ptail= l-Fr(u)〇
[0070] 上述模型中需要估计的参数为第一概率Ptail和〇 ^
[0071] 而Ptail的参数估计值为每一个重采样样本的经验分布：
[0072] 也即，ptail=m/n
[0073] 采用极大似然法估计〇 i，其估计值为：
[0074] 〇 !=E(Y)
[0075] 即可根据估计出的第一概率Ptail和第二标准差〇 :确定第一分位数。
[0076] 可选的，根据第一概率和第二标准差确定第一分位数，包括：
[0077] 对估计出的PtaijP 〇的上限进行置信度检验。
[0078] 令根据B组重采样样本估计的参数值rtail，b}和{ <}的经验a分位数为Ptail 和〇的置信度为a的上限。
[0079] 根据参数置信限值flu(㈨和参数置信限值 岣确定放大因子kinf。
[0080] 进一步的，将估计出的第一概率Ptail对应的第一参数置信限值P#tail， B(a )和第二 标准差〇 i对应的第一参数置信限值〇 a )带入：
[0081] Fr(x) = (l-Ptail)+PtailFe(x_u)，x>u
[0082] 得到GBAS伪距误差的概率分布：
[0084] 最后根据第一参数置信限值P#tail，B(a)和第二参数置信限值〇\(a)确定第一分 位数。
[0085] 其中，根据下述公式计算第一分位数：
[0087] 其中，Qr，risk为第一分位数，u。为阈值，P' ail，B(a )为第一参数置信限值，〇 '(a ) 为第二参数置信限值，P&k为GBAS对应的指定风险。
[0088] 步骤104 :根据地基增强系统GBAS总的运行风险确定标准高斯分布对应GBAS总 的运行风险对应的第二分位数。
[0089] 步骤105 :根据第一分位数、第二分位数和第一标准差确定本组伪距误差校正样 本值对应的放大因子。
[0090] 具体的，根据第一分位数、第二分位数和第一标准差确定放大因子kinf，包括：
[0091] 根据以下公式确定放大因子kinf:
[0092] kinf - Q r, risk/kffmd/0samp
[0093] 其中，Quisk为第一分位数、kffDld为第二分位数，a SMP为第一标准差。
[0094] 本发明实施例提供一种基于极值理论的卫星导航地基增强系统的误差包络方法， 包括：获取至少两组伪距误差校正样本值，并确定各组伪距误差校正样本值对应的具有指 定风险的第一分位数和第一标准差，然后根据地基增强系统GBAS总的运行风险确定在标 准高斯分布下GBAS总的运行风险对应的第二分位数，进而根据第一分位数、第二分位数和 第一标准差确定本组伪距误差校正样本值对应的放大因子，最后，根据放大因子确定本组 伪距误差校正样本值对应的伪距校正值，通过上述方法可以准确计算出放大因子，从而计 算的伪距校正值能够包络实际误差，有效提高系统的连续性。
[0095] 本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通 过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程 序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：只读存储记忆 体（Read-Only Memory，ROM)、随机存储记忆体（Random Access Memory，RAM)、磁碟或者光 盘等各种可以存储程序代码的介质。
[0096] 最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制； 尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其 依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征 进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技 术方案的范围。
【主权项】
1. 一种基于极值理论的卫星导航地基增强系统的误差包络方法，其特征在于，包括： 获取至少两组伪距误差校正样本值； 对每组伪距误差校正样本值执行下述操作： 确定本组伪距误差校正样本值对应的具有指定风险的第一分位数和第一标准差； 根据地基增强系统GBAS总的运行风险确定在标准高斯分布下所述GBAS总的运行风险 对应的第二分位数； 根据所述第一分位数、所述第二分位数和所述第一标准差确定本组伪距误差校正样本 值对应的放大因子。2. 根据权利要求1所述的方法，其特征在于，所述确定本组伪距误差校正样本值对应 的具有指定风险的第一分位数，包括： 确定所述本组伪距误差校正样本值对应的阈值； 根据重采样法获得所述本组伪距误差校正样本值对应的至少两组重采样样本值； 根据所述本组伪距误差校正样本值对应的阈值和所述本组伪距误差校正样本值对应 的所述至少两组重采样样本值中的每组重采样样本值确定每组重采样样本值中的超出量 样本值，所述超出量样本值为所述每组重采样样本值中大于所述阈值的值与所述阈值的差 值； 根据所述超出量样本值确定所述第一分位数。3. 根据权利要求2所述的方法，其特征在于，所述根据所述超出量样本值确定所述第 一分位数，包括： 根据本组重采样样本值对应的超出量样本值确定所述本组重采样样本值中超过所述 阈值的第一概率； 确定所述本组重采样样本值对应的超出量样本值的第二标准差； 根据所述第一概率和所述第二标准差确定所述第一分位数。4. 根据权利要求3所述的方法，其特征在于，所述根据所述第一概率和所述第二标准 差确定所述第一分位数，包括： 根据所述第一概率确定指定置信度的第一参数置信限值； 根据所述第二标准差确定指定置信度的第二参数置信限值； 根据所述第一参数置信限值和所述第二参数置信限值确定所述第一分位数。5. 根据权利要求2所述的方法，其特征在于，所述确定所述本组伪距误差校正样本值 对应的阈值，包括： 根据平均超出量函数MEF确定所述本组伪距误差校正样本值对应的阈值。6. 根据权利要求1所述的方法，其特征在于，所述获取至少两组伪距误差校正样本值， 包括： 根据预设时间间隔获取至少两个伪距误差校正样本值； 对所述至少两个伪距误差校正样本值进行分组得到至少两组伪距误差校正样本值。
【专利摘要】本发明提供一种基于极值理论的卫星导航地基增强系统的误差包络方法，包括：获取至少两组伪距误差校正样本值，并确定各组伪距误差校正样本值对应的具有指定风险的第一分位数和第一标准差，然后根据地基增强系统GBAS总的运行风险确定在标准高斯分布下GBAS总的运行风险对应的第二分位数，进而根据第一分位数、第二分位数和第一标准差确定本组伪距误差校正样本值对应的放大因子，通过上述方法可以准确计算出放大因子，从而计算的伪距校正值能够包络实际误差，有效提高系统的连续性。
【IPC分类】G01S19/41
【公开号】CN104898140
【申请号】CN201510284335
【发明人】薛瑞, 张军, 朱衍波, 王志鹏
【申请人】北京航空航天大学
【公开日】2015年9月9日
【申请日】2015年5月28日
最新回复(0)