h、辨识第m个测站采用的用户测距精度表示形式，令k = l;
[0032] hi、构造 b = [bi，b2,…，bi]，其中
[0033]
[0034] 112、令￥1{ = ；^'131'，其中131'表示13的转置向量；
[0035] h3、令 k = k+l
[0036] h4、当k < 4时，返回步骤hi，否则执行步骤i ;
[0037] ? N'4'Wmax = max{wi ,W2 ,W3 ,W4}
[0038] il、如果Wmax〈0.9,将置为空集，令m=m+l，执行步骤b;
[0039] i2、否则，令 k = 4;
[0040] i3、如果Wmax=Wk，则用户测距精度采用的是Lk对应的用户测距精度表达形式，执行 步骤j;
[0041 ] i4、如果Wmax关Wk，令k = k_l，执行步骤i3;
[0042] j、将第m个测站的用户测距精度替换为上界值，令k0 = l;
[0043] jl、读取Sm中的第k0个元素 sk0;
[0044] j2、如果skQ>6144· 00，令T = 101%执行步骤j5，否则，令kl = 1，执行步骤j3;
[0045] j3、如果sk〇 = Lk(kl)，令J = .，执行步骤j5,其中Lk(kl)表示向量Lk的第kl个 元素，k的取值在执行步骤j前已确定，1 < k < 4,否则，令kl =kl+l，执行步骤j4;
[0046] j4、如果kl < 15,执行步骤j3,否则，如果2且kl 2 16,或者k>2且kl>15,令 歹=1015'，执行步骤j5;
[0047] j5、将广播星历EkQ中的用户测距精度替换为:Γ;
[0048] j6、将广播星历Ek〇中6个调和改正项幅值0>。、(^、(^、(^、(^、(^、平动点偏差八11、 参考历元的平近点角Mo、轨道偏心率e、按参考历元计算的升交点赤经Ω 〇、升交点赤经变化 率_￡^、近地点角ω、轨道半长轴平方根以及参考历元的轨道倾角i〇,共Η个参数作为稳 健星历参数向量ER,kQ，而将卫星编号PRN、星历参考历元t〇 e、卫星钟差参考历元t。。、卫星健康
[0049] 状况标记health_if、用户测距精度上界值URA_UB作为非稳健星历参数向量，即
[0050]
[0051] 卫星编号PRN、星历参考历元toe、卫星钟差参考历元t。。、卫星健康状况标记 health_if、用户测距精度上界值1^_1^在％:;1〇中的顺序依次为1，2，3，4，5;
[0052] j7、将集合I第k0个广播星历Ek〇替换为=[五^。五_]:;
[0053] j8、令k0 = k0+l，如果ko < n，执行步骤jl，否则执行步骤k;
[0054] k、令 m=m+l，当m<N，执行步骤 b;
[0055] m、将N个测站所有广播星历Mm，m = 1，2，…，N，存入集合Μ中，
[0056]
b，其中，nm表示第m个测站有^组广播星历数据，ρ 为Μ的元素总数；
[0057] η、将第一组广播星历数据E1的稳健参数Eiu存入爲中，非稳健参数士4存入^中，
[0058] nl、令 k = 2，q = l;
[0059] n2、令i = l;
[0060] n3、读取Μ中的广播星历数据Ek;
[0061] n4、如果Ek的稳健星历数据部分中的第i组稳健星历数据馬~相同，则
[0062]
[0063] 令k = k+l，当k < p时，执行步骤n2,当k>p时，执行步骤〇;
[0064] n5、否贝I】，令i = i+l，如果i < q，执行步骤n4,如果i>q执行步骤n6;
[0065] n6、令q = q+l，将ER.kJu.分别作为稳健、非稳健星历参数集合晃、馬的第q个 元素，即
[0066] 步骤n4中ER,k与^^相同是指ER,k和ffl,中所有对应的广播星历数据差的绝对值均 不大于相应的阈值，相应的阈值如表1所示
[0067] 表 1
[0068]
[0069]
[0070] ο、构造集合巧，·以及向量Y;
[0071] 〇1、令巧为空集；
[0072] 〇2、令 k = l;
[0073] 〇3、令 p = l;
[0074] 〇4、读取集合&中的第k个元素巧.s ;
[0075] 〇5、如果巧^的行数小于3,令k = k+l，当时，其中中元素的个数，执行 步骤〇4,当k>r时，执行步骤p;
[0076] 〇6、否则，当馬』的行数为偶数时，删除·的最后一行；
[0077] 〇7、取^^中第1至4列的中位数以及第5列的最大值按顺序存入向量T，令 y - row ?馬ul，其中 row (鳥」表示矩阵馬的行数；
[0078] 〇8、如果g为空集，令的第一个元素向量γ的第一个元素 Yl = y，令k = 1^+1，当1^<1'时，执行步骤〇4，当1^>1'时，执行步骤口；
[0079] 〇9、当<^为非空集合时，令i = l;
[0080] 〇10、读取1?中的第i个元素 g,,;
[0081] 〇11、如果巧^的第一、第二个元素分别与T的第一个、第二个元素均相等，执行步骤 〇12,否则令i = i+l，当i < p时，执行步骤〇10,当i>p时，执行步骤〇14;
[0082] 〇12、如果Y的第i个元素￥々，令^%:=了小=7，执行步骤〇13;
[0083] 〇13、令k = k+l，当k < r时执行步骤〇4,当k>r时，执行步骤p;
[0084] 〇14、令p = p+l，将T、y分别作为g、Y的第p个元素，即，_Yp = y，执行步骤 〇15；
[0085] 〇15、令k = k+l，当k < r时执行步骤〇4,当k>r时，执行步骤p;
[0086] p、通过上述步骤得到用于卫星导航系统完好性监测评估的广播星历数据集合 A，它的每个元素包括卫星编号PRN、星历参考历元toe、卫星钟差参考历元t。。、卫星健康状 况标记heal th_if、用户测距精度上界值URAJJB。
[0087] 本发明的优点是:提出了一种基于频次加权识别单站用户测距精度(URA)表达形 式的方法，并利用多测站用户测距精度上界值(URAJJB)合成最终的用户测距精度上界值 (URAJJB)，与IGS提供的brdc合成星历以及相关文献提供的合成结果相比，本方法合成的用 户测距精度上界值(URA_UB)更为合理、准确，为提升卫星导航系统性能提供更科学的参考 依据。
【附图说明】
[0088] 图1是多测站广播星历确定用户测距精度上界值(URA_UB)的方法流程示意图。
【具体实施方式】
[0089] 以下结合附图及附表，对本发明作进一步的详细说明。
[0090] 实施例一
[0091] 多测站广播星历确定用户测距精度上界值(URA_UB)的方法，该方法包含下列步 骤：
[0092] 多测站广播星历确定用户测距精度上界值的方法，该方法是在有N个可以提供广 播星历数据的GNSS测站的条件下，确定用户测距精度的上界值，具体步骤如下：
[0093] a、给每个测站一个序号，第m个测站是指序号为m的测站，m=l，2,…，N;
[0094] al、$m=l;
[0095] b、如果m>N，执行步骤m，否则读取第m个测站的广播星历数据，存入集合 氣,=_1,财2,-,^^中，其中1表示该测站第1组广播星历数据，1 = 1，2，~，11，11表示第111个 测站的广播星历组数；
[0096] c、将第m个测站广播星历中的用户测距精度存入向量5"中，5"=[幻，82，"_，~]，其 中81为第i组广播星历吣中的用户测距精度，i = l，2,…，η;
[0097] d、构造向量S'，向量S'的第一个元素为=
[0098] dl、令 i = 2，l = l;
[0099] d2、读取S中的第i个元素 Sl，令j = l;
[0100] d3、读取S'中的第j个元素 s'J;
[0101] d4、当Si = s'j，令i = i+l，如果i>n执行步骤e，否则执行步骤d2,当Si关s'j，令j = j +1，如果j>l，执行步骤d5，否则执行步骤d3;
[0102] d5、令 1 = 1+1，5'的第1个元素8/1 = &;
[0103] e、统计S'的元素 Y i在Sm中出现的次数，记为fi，其中i = l，2,…，1，组成向量f = [fi,f2,··· ,fi];
[0104] f、对f进行归一化处理，得到f·'
[0105]
[0106] g、令1^ = URA_IX = [0，1，…，15]，称URA_IX为用户测距精度指数向量，
[0107] L2 = URA_LB= [0.00,2.40,3.40,4.85,6.85,9.65,13.65,24.00,48.00,96.00, 192 · 00,384· 00,768 · 00,1536 · 00,3072 · 00,6144· 00]，称URA_LB为用户测距精度下界向量，
[0108] L3 = URA_NV= [2.00,2.83,4.00,5.66,8.00,11.31,16.00,32.00,64.00,128.00, 256 · 00,512 · 00,1024· 00,2048 · 00,4096 · 00]，称URA_NV为用户测距精度名义向量，
[0109] L4 = URA_UB= [2.40,3.40,4.85,6.85,9.65,13.65,24.00,48.00,96.00,192.00, 384· 00,768 · 00,1536 · 00,3072 · 00,6144· 00]