一种用于局域增强系统完好性监测的门限值确定方法
【专利摘要】本发明提供了一种用于局域增强系统完好性监测的门限值确定方法，根据局域增强系统的虚警率要求计算对应的分位数；然后按照实际使用需要将卫星仰角区间进行划分，并将采集到的观测量数据按区间进行分组；计算每组数据的均值和标准差，并进行归一化处理；然后利用膨胀的高斯分布去包络归一化数据的概率密度分布从而计算出膨胀因子；利用分位数、均值、标准差和膨胀因子计算出每个区间的门限值；最后利用各区间的门限值进行拟合得到门限值曲线。本发明满足了局域增强系统对虚警率的要求；同时，因为是使用系统本身采集的数据来确定门限值，也更加符合系统的实际特性。
【专利说明】-种用于局域增强系统完好性监测的门限值确定方法
【技术领域】
[0001] 本发明属于卫星导航领域，涉及一种用于局域增强系统完好性监测的门限值确定 方法。
【背景技术】
[0002] 卫星导航具有全天候、广覆盖和低成本等优点，已经成为导航体系中最主要的导 航手段。随着导航应用的不断深入，单纯利用卫星导航系统传统的单点定位服务已经不能 满足各种用户所需的精度和完好性等要求。这就需要建立增强系统对导航性能进行增强， 局域增强系统就是其中之一。局域增强系统被广泛应用在机场对飞机的进近引导上，它是 一种地基增强系统，其地面基准站利用多个已知精确位置的基准接收机（通常为3至4个） 接收可视卫星的测距信号，综合计算出每颗卫星的码伪距校正量、载波相位校正量等增强 信息，并将这些信息发送给机载子系统。机载子系统利用差分增强原理，在定位解算中用这 些信息来减小相关误差的影响，可以获得性能得到增强的导航服务。
[0003] 精密进近对定位精度和完好性要求很高，II类、III类精密进近对完好性风险的 要求达到2XKT9/进近，即每10亿次进近中，发生引导信息不可用且机载子系统没有得到 报警的次数不得超过2次。苛刻的完好性风险要求局域增强系统在将各种增强信息发送给 机载子系统之前，必须要进行完好性监测，避免将不可用的引导信息发出。导致引导信息不 可用的因素主要包括地面站设备和空间卫星信号的故障及异常，对应到完好性监测中，就 要进行相应的质量监测，包括信号质量监测，数据质量监测和测量值质量监测。在局域增强 系统中，进行以上质量监测的做法是用地面站基准接收机采集的各种观测量与其相应的门 限值进行比较，若观测量的大小超出了门限值，则认为该观测量发生了异常，从而可以进行 异常值剔除或进一步确定对应的观测通道出现了故障。
[0004] 在进行质量监测时，门限值的确定成为关键。若门限值取得过大，则异常观测量无 法被完全检测出，会造成错误的引导信息；若门限值取得过小，则正常观测量也会被当做异 常值处理，系统的连续性会受到影响。在局域增强系统完好性监测中，观测量质量监测对门 限值的要求是：让正常观测量超出门限值范围的概率小到一定程度（通常要求达到1〇_ 9的 量级）以满足系统整体连续性的要求，此概率称之为虚警率。以往确定门限值的做法主要 有两种：一种是依靠工程实践经验为各种观测量设定门限值，这种做法缺少理论依据而无 法保证虚警率，同时，以往的经验通常都是基于GPS卫星及其接收机的，对于我国北斗导航 系统及其接收机，可能并不适用。另一种是根据各种观测量门限值模型来确定门限值，但由 于实际使用时的各种条件（如地面基准站的位置、接收机的性能、环境的差异、所用的卫星 导航系统）可能与门限值模型建立时有较大的出入，因此按照模型来确定门限值的做法并 不一定正确。通常，这两种做法被结合起来使用，但还是无法保证完全适用。
【发明内容】
[0005] 为了克服现有技术的不足，本发明提供一种用于局域增强系统完好性监测的门限 值确定方法，根据观测量数据的概率密度分布，结合虚警率要求计算出对应的门限值。
[0006] 本发明解决其技术问题所采用的技术方案包括以下步骤：
[0007](1)计算局域增强系统要求的虚警率p对应的分位数a= (p/2)，其中函数Q(s) 为标准正态分布S的分位数函数，s彡0,Q(S)代表S的取值大于s的概率，即Q(S) =P(S>s);
[0008]⑵将卫星仰角按照L度作为一区间进行划分，则得到的各个仰角区间为（0,L]， (L，2L]，（2L，3L]…（90-L，90]，共有M= 90/L个仰角区间，L取5度或10度；
[0009] (3)对采集到的观测量数据，将属于同一个仰角区间的分为一组，得到M组观测量 数据D1, D2, Dm;
[0010] (4)计算每组观测量数据的均值iii，ii2，ii3*"iiM和标准差〇i，〇 2，〇 3*" 〇 M;
[0011] (5)对每一组观测量数据Di的每一个数据点减去该组数据的均值yi后，再除以 该组观测量数据的标准差Oi，得到M组经过归一化处理的数据Dn1,Dn2,Drv"DnM,i= 1, 2,3…M;
[0012] (6)计算M组经过归一化处理的数据Dn1, Dn2, Drv"DnM的膨胀因子f\，f2, ff fM ; 具体做法是将每一组经过归一化处理的数据的概率密度分布Yni和高斯分布Gni?N (0, f2) 画于同一幅图中，并让f的初值为1并以设定步长增长，直到所画出的高斯分布能够刚好包 络住概率密度分布Yni的两侧尾部，此时的f值即为对应的膨胀因子&;
[0013]⑵计算每个仰角区间的门限值Ti = ii 〇i(i=1，2, 3*"M);
[0014] (8)按照各个仰角区间的门限值进行曲线拟合，得到门限值曲线；在实际使用时， 根据卫星仰角在曲线上取对应仰角下的值作为门限值即可。
[0015] 本发明的有益效果是：
[0016](1)本发明提出的门限值确定方法，相比于以往依靠工程经验或门限值模型的做 法，按照虚警率计算对应的分位数，满足了局域增强系统对虚警率的要求；
[0017](2)本发明提出的门限值确定方法，相比于以往依靠工程经验或门限值模型的做 法，使用系统本身采集的统计数据来确定门限值，更符合系统的实际特性。
【专利附图】
【附图说明】
[0018] 图1是本发明进行门限值确定的处理流程图；
[0019] 图2是具体实施例中的原始载波加速度数据示意图；
[0020] 图3是具体实施例中的进行了归一化处理后的载波加速度数据示意图；
[0021] 图4是具体实施例中进行膨胀因子计算的示意图；
[0022] 图5是具体实施例中计算得到的门限值曲线图。
【具体实施方式】
[0023] 下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施 例。
[0024] 本发明利用高斯分布去包络住观测量数据的概率密度分布，并结合虚警率要求计 算出对应的门限值，具体步骤如下：
[0025] (1)计算虚警率对应的分位数
[0026] 局域增强系统要求的虚警率为p，对高斯分布X?N(y，〇2)，取门限值为T= U±a* 〇，其中a为虚警率对应的分位数（a>0)，则X取值超过门限值的概率为P(|X-ii >a*〇) =P(|S| >a) =p，其中S为标准正态分布，S?N(0，1)，则可计算得虚警率对应 的分位数a= (T(p/2)，其中函数Q(S)为标准正态分布S的分位数函数（s>0)，Q(S)代 表S的取值大于s的概率，即Q(S) =P(S>s)，GT1为Q的逆函数。
[0027] (2)卫星仰角区间划分
[0028] 局域增强系统完好性监测中各种观测量的大部分统计特性都与卫星的仰角有关。 仰角越低，多径现象越严重，电离层和对流层的倾斜因子越大，导致观测量误差也越大。因 此在确定门限值时，有必要对卫星仰角区间进行划分。卫星仰角在〇到90度之间变化，将L 度的仰角范围划分为一个区间，则得到的各个仰角区间为（〇，L]，（L，2L]，（2L，3L]…（90-L， 90]，共有M= 90/L个仰角区间，L一般取5度或10度（通常卫星仰角小于5度的数据并 不使用）。
[0029] (3)观测量数据分组
[0030] 对采集到的观测量数据，将属于同一个仰角区间的分为一组，得到M组观测量数 据D1,D2,D3 …Dm。
[0031] (4)计算每组数据的均值和标准差
[0032] 对按照卫星仰角区间分组得到的M组观测量数据D1,D2,D^Dm,计算每组数据的 均值UpU2，Ii 3...IiM和标准差 〇 1，〇 2，〇 3... 〇M。
[0033] (5)归一化处理
[0034] 对按照卫星仰角区间分组得到的M组观测量数据进行归一化处理。即对每一个数 据组DiQ= 1，2,3…M)，将其每一个数据点减去该组数据的均值1^后，再除以该组数据的 标准差〇i，得到M组经过归一化处理的数据Dn1,Dn2,Drv"DnM。
[0035] (6)计算膨胀因子
[0036] 对M组经过归一化处理的数据Dn1,Dn2,Drv"DnM，计算其对应的膨胀因子f\，f2， fV"fM。具体做法是：对每一个数据组DniQ= 1，2, 3…M)，将其概率密度分布Yni和高斯 分布Gni?N(0,f2)画于同一幅图中，并让f?的初值为1并以一定步长增长，直到所画出的 高斯分布能够刚好包络住分布Yni的两侧尾部（越到外侧概率越小，超出门限值范围的数 据都落在此区域），此时的f值即为对应的膨胀因子。此时，在两侧尾部，分布Yni的概率 值小于高斯分布Gni?N(0,A2)的概率值。
[0037] (7)计算每个仰角区间的门限值
[0038] 对按照卫星仰角划分得到的M个仰角区间，每个区间的门限值取为Ti= ii 〇Ji= 1,2,3*"M),其中iii和〇i为步骤⑷中计算出的第i个仰角区间数 据的均值和标准差，a为步骤（1)中计算出的虚警率对应的分位数，&为步骤￠)中计算出 的第i个区间的膨胀因子。经过步骤￠)中膨胀因子的计算，Di的概率密度分布Yi的两侧 尾部能够被高斯分布Gi?N(yi，（A* 〇J2)包络住。此时有Di超过其门限值Ti的概率为 PidDi-iiJ=PJlSl>a) =2Q(a) =p，其中S 和Q的含义与步骤（1)中的相同。即有Pi(IDi-ii」>a*fjOi) <p，满足Di超过其门限 值Ti的概率小于虚警率p的要求。
[0039] (8)门限值曲线拟合
[0040] 经过以上步骤已经得到各个仰角区间的门限值Ti(i= 1，2,3…M)，进行曲线拟合 后得到门限值曲线。在实际使用时，根据卫星仰角在曲线上取对应仰角下的值作为门限值 即可。
[0041] 下面以局域增强系统完好性监测中测量值质量监测的载波加速度监测为例，对本 发明的门限值确定方法的【具体实施方式】进行说明。载波加速度监测属于局域增强系统完 好性监测中对载波相位测量值的一种监测，它利用原始载波相位测量值计算其在各个历元 间变化的加速度，然后与门限值比较，若加速度超出门限，则认为载波相位测量值出现了异 常，从而完成监测。利用本发明提出的方法，在某次实际使用中确定载波加速度监测所使用 的门限值的具体步骤如下：
[0042] (1)计算虚警率对应的分位数
[0043] 局域增强系统的虚警率取为？ = 1父10'则其对应的分位数为3 = (^(?/2)= Q-1(IX10_9/2) 6. 1095〇
[0044] (2)卫星仰角区间划分
[0045]将卫星仰角按照10度一个区间进行划分，则得到9个仰角区间，分别为（0, 10]、 (10,20]、 （20,30]、 （30,40]、 （40,50]、 （50,60]、 （60,70]、 （70,80]和（80,90]。
[0046] (3)观测量数据分组
[0047] 将采集得到的载波相位加速度数据按照步骤（2)中的仰角区间进行分组，得到9 组加速度数据AcciQ= 1，2,3…9)，如图2所示。在图2中，由于实际采集数据的局域增强 系统地面站的地理位置的限制，卫星仰角小于10度的情况基本上观测不到，因此在本实施 例中对于仰角从〇到10度的区间不作处理。
[0048] (4)计算每组数据的均值和标准差
[0049]对每组载波加速度数据AcciQ= 2,3*"9)，计算出其对应的均值iiJi= 1，2,3… 9)和标准差〇i(i= 2, 9)，结果如表1所示。
[0050] 表1各仰角区间加速度数据的均值和标准差
[0051]
【权利要求】
1. 一种用于局域增强系统完好性监测的门限值确定方法，其特征在于包括下述步骤： (1)计算局域增强系统要求的虚警率P对应的分位数a = (Γ (p/2)，其中函数Q (S)为 标准正态分布s的分位数函数，s彡0, Q (s)代表S的取值大于s的概率，即Q (s) = P (S > s); ⑵将卫星仰角按照L度作为一区间进行划分，则得到的各个仰角区间为（0, L]，（L， 2L]，（2L，3L]…（90-L，90]，共有M = 90/L个仰角区间，L取5度或10度； (3) 对采集到的观测量数据，将属于同一个仰角区间的分为一组，得到M组观测量数据 D1, D2, D3 …Dm ; (4) 计算每组观测量数据的均值μ P μ 2, μ 3…μ M和标准差〇 i，σ 2，σ 3··· σ M ; (5) 对每一组观测量数据Di的每一个数据点减去该组数据的均值Ui后，再除以该组 观测量数据的标准差〇 i，得到M组经过归一化处理的数据Dn1, Dn2, Drv"DnM, i = 1，2, 3··· M ； (6) 计算M组经过归一化处理的数据DnpDrvDrvDnM的膨胀因子f2, fV"fM ;具体 做法是将每一组经过归一化处理的数据的概率密度分布Yni和高斯分布Gni?N(0, f2)画 于同一幅图中，并让f的初值为1并以设定步长增长，直到所画出的高斯分布能够刚好包络 住概率密度分布Yn i的两侧尾部，此时的f值即为对应的膨胀因子& ; ⑵计算每个仰角区间的门限值Ti = μ = 1，2,3···Μ); (8)按照各个仰角区间的门限值进行曲线拟合，得到门限值曲线；在实际使用时，根据 卫星仰角在曲线上取对应仰角下的值作为门限值即可。
【文档编号】G01S19/42GK104318089SQ201410537847
【公开日】2015年1月28日 申请日期:2014年10月13日 优先权日:2014年10月13日
【发明者】熊帅, 吕超, 李斌, 王晓旺, 路阳 申请人:中国电子科技集团公司第二十研究所