CN105467406B - 多测站广播星历确定用户测距精度上界值的方法 - Google Patents
多测站广播星历确定用户测距精度上界值的方法 Download PDFInfo
- Publication number
- CN105467406B CN105467406B CN201510966601.2A CN201510966601A CN105467406B CN 105467406 B CN105467406 B CN 105467406B CN 201510966601 A CN201510966601 A CN 201510966601A CN 105467406 B CN105467406 B CN 105467406B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- range accuracy
- user range
- ura
- 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
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)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种多测站广播星历确定用户测距精度上界值(URA_UB)的方法,涉及卫星导航系统性能的监测与评估。该方法采用频次加权模式识别确定单个测站用户测距精度(URA)表示形式,再根据多个测站广播星历中的用户测距精度(URA)值确定对应时刻、对应卫星的用户测距精度上界值(URA_UB)。本发明合成的用户测距精度上界值可用于全球卫星导航系统空间信号的完好性监测与评估,为提升卫星导航系统性能提供更科学的参考依据。
Description
技术领域
本发明涉及卫星导航系统性能监测与评估,更具体涉及合成GPS广播星历时,确定用户测距精度(User Range Accuracy-URA)参数统一表达形式以及上界值的方法,该方法适用于GPS卫星导航系统完好性监测与评估的广播星历合成。
背景技术
精度、完好性、连续性与可用性是卫星导航系统四个重要性能指标,其中完好性对民航等涉及生命安全的应用领域至关重要,其性能监测与评估备受关注。而GPS空间信号(Signal in Space–SIS)完好性的监测与评估是GPS系统完好性保障中的重要组成部分,它和广播星历中的用户测距精度(URA)指标以及卫星健康标记密切相关:GPS现代化之前,当卫星健康标记为0时,用户测距误差(User Range Error-URE)超过4.42×URA_UB(URA上界值),但未及时报警的概率小于10-5/h(称之为完好性风险);GPS现代化以后,当完好性状态标记(Integrity Status Flag-ISF)为0时,完好性风险与现代化前的相同,当完好性标记(ISF)为1(即为“增强”模式)时,用户测距误差(URE)超过5.73倍的用户测距精度上界值(URA Upper Bound value–URA_UB),但未及时报警的概率小于10-8/h。由此可见,用户测距精度上界值(URA_UB)在卫星导航系统的完好性保障、完好性监测与评估中具有关键作用。
然而,现在合成的广播星历中还没有一种统一的用户测距精度上界值可用于GPS空间信号完好性监测与评估。广播电文中的用户测距精度是四位二进制数表示的指数形式(URAIndex–URA_IX),除此外,还有用户测距精度上界值(URA_UB),用户测距精度下界(URALowerBound value-URA_LB)、用户测距精度名义值(URA Normal Value-URA_NV)两种形式。四种表示形式对应的值如下表所示(GPS Navstar Joint Program Office,NavstarGPS Space Segment/Navigation User Interfaces,IS-GPS-200,Revision H,GPSNavstarJoint Program Office,El Segundo,CA,Sep.7,2013):
表2用户测距精度的四种不同表示形式
表2中,当用户测距精度指数等于15时,表示无可用的预测准确度,用户使用该卫星定位时,风险自担。
各接收机厂家采用的处理方式不同,导致其输出的URA值含义不尽相同,例如国际GNSS服务组织(IGS)有些测站获取的广播星历中的用户测距精度采用名义值形式,如测站dgar的接收机(类型:ASHTECH UZ-12)、测站dgav的接收机(类型:JAVAD TRE_G3TH),而有些采用指数形式,如测站nist的接收机(类型:NOV OEM4-G2)、测站mat1的接收机(类型:TRIMBLE 4000SSI),此外还有些接收机采用下界值、上界值作为用户测距精度的表示形式。而国际GNSS服务组织(IGS)提供的合成广播星历(brdc)中的用户测距精度(URA)参数表达形式也不统一,比如2015年9月17日(年积日260天)至2015年10月16日(年积日289天)间,国际GNSS服务组织(IGS)提供的GPS合成广播星历(brdc)中的用户测距精度值主要是{0,2,2.8,2.9,4,5.7,8,4096}中的元素,比如2015年9月20日(年积日263天)的用户测距精度值有:0,2,2.8,2.9,4,5.7,由这些值很难判断是哪种用户测距精度的表达形式,2015年9月16日(年积日262天)中的用户测距精度值甚至有负值,即使卫星健康标记为不健康,这也是不合理的。由此可见,国际GNSS服务组织(IGS)提供的GPS合成广播星历是不能用于该系统的完好性监测与评估的。
此外,GPS卫星可能受到周围环境、星上载荷工作状态影响,导致播发的数据与所注入数据不同,即使GPS运行控制系统的地面监测站能够监测到异常,一般需要45分钟,甚至更长的时间作出反应。因此地面运控端的注入星历也不能作为GPS系统完好性评估的基础数据。
在解决这类难题的方法中,有学者(Heng L.,Gao G.X.,WalterT.,et al.(2010).GPS Ephemeris Error Screening and Results for 2006–2009[C].ION2010International Technical Meeting January 25-27,2010,San Diego,CA.;Heng L.,Gao G.X.,Walter T.,et al.(2011).Digging into GPS Integrity-Charting theEvolution of Signal-in-Space Performance by Data Mining 400,000,000Navigation Messages[J].GPS World,22(11):44-49.;Heng L.,Gao G.X.,Walter T.,etal.(2012).GPS Signal-in-Space Integrity Performance Evolution in the LastDecade[J].IEEE Transactions on Aerospace and Electronic Systems,48(4):2932-2946)提出采用基于模式识别的方法分五步判定导航星历文件中提取的用户测距精度(URA)数据属于哪种表示形式,并对2004年至2010年的广播星历数据进行了合成。虽对同一台接收机来说,其解码得到的所有卫星URA值表示形式是统一的,这也易于采用模式识别方法判定,但是仍有少数国际GNSS服务组织(International GNSS Service–IGS)跟踪站的导航星历文件中的用户测距精度(URA)值并不是上述四种表示形式中的任一个值。而上述文献所提出的模式识别法将这类数值认为是除用户测距精度指数(URA_IX)表达形式外的其它形式,并将其归入差值绝对值最近的用户测距精度(URA)表示形式。这种处理方法忽略了导航文件中用户测距精度(URA)不同数值出现的频次,实际上异类用户测距精度(URA)值出现的次数很少。因此,前述文献提出的用户测距精度(URA)表示形式判定的模式识别法在实际应用中存在不足。
发明内容
本发明的目的是:提供一种多测站广播星历确定用户测距精度上界值(URA_UB)的方法,该方法采用频次加权模式识别确定单个测站用户测距精度(URA)表示形式,再根据多个测站广播星历中的用户测距精度(URA)值确定对应时刻、对应卫星的用户测距精度上界值(URA_UB)。该方法合成的用户测距精度上界值可用于全球卫星导航系统空间信号的完好性监测与评估,为提升卫星导航系统性能提供更科学的参考依据。
为达到上述目的,本发明采用了如下技术方案:
1、多测站广播星历确定用户测距精度上界值的方法,该方法是在有N个可以提供广播星历数据的GNSS测站的条件下,确定用户测距精度的上界值,具体步骤如下:
a、给每个测站一个序号,第m个测站是指序号为m的测站,m=1,2,…,N;
a1、令m=1;
b、如果m>N,执行步骤m,否则读取第m个测站的广播星历数据,存入集合中,其中Mi表示该测站第i组广播星历数据,i=1,2,…,nm,nm表示第m个测站的广播星历组数;
c、将第m个测站广播星历中的用户测距精度存入向量Sm中,其中si为第i组广播星历Mi中的用户测距精度,i=1,2,…,nm;
d、构造向量S',向量S'的第一个元素为s'1,s'1=s1;
d1、令i=2,l=1;
d2、读取Sm中的第i个元素si,令j=1;
d3、读取S'中的第j个元素s'j;
d4、当si=s'j,令i=i+1,如果i>nm执行步骤e,否则执行步骤d2,当si≠s'j,令j=j+1,如果j>l,执行步骤d5,否则执行步骤d3;
d5、令l=l+1,S'的第l个元素s'l=si;
e、统计S'的元素s'i在Sm中出现的次数,记为fi,其中i=1,2,…,l,组成向量f=[f1,f2,…,fl];
f、对f进行归一化处理,得到f'
g、令L1=URA_IX=[0,1,…,15],称URA_IX为用户测距精度指数向量,
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为用户测距精度下界向量,
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为用户测距精度名义向量,
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],称URA_UB为用户测距精度上界向量;
h、辨识第m个测站采用的用户测距精度表示形式,令k=1;
h1、构造b=[b1,b2,…,bl],其中
h2、令wk=f'bT,其中bT表示b的转置向量;
h3、令k=k+1
h4、当k≤4时,返回步骤h1,否则执行步骤i;
i、令wmax=max{w1,w2,w3,w4}
i1、如果wmax<0.9,将令m=m+1,执行步骤b;
i2、否则,令k=4;
i3、如果wmax=wk,则用户测距精度采用的是Lk对应的用户测距精度表达形式,执行步骤j;
i4、如果wmax≠wk,令k=k-1,执行步骤i3;
j、将第m个测站的用户测距精度替换为上界值,令k0=1;
j1、读取Sm中的第k0个元素sk0;
j2、如果sk0>6144.00,令执行步骤j5,否则,令k1=1,执行步骤j3;
j3、如果sk0=Lk(k1),令执行步骤j5,其中Lk(k1)表示向量Lk的第k1个元素,k的取值在执行步骤j前已确定,1≤k≤4,否则,令k1=k1+1,执行步骤j4;
j4、如果k1≤15,执行步骤j3,否则,如果k≤2且k1≥16,或者k>2且k1>15,令执行步骤j5;
j5、将广播星历Ek0中的用户测距精度替换为
j6、将广播星历Ek0中6个调和改正项幅值Cuc、Cus、Crc、Crs、Cic、Cis、平动点偏差Δn、参考历元的平近点角M0、轨道偏心率e、按参考历元计算的升交点赤经Ω0、升交点赤经变化率近地点角ω、轨道半长轴平方根以及参考历元的轨道倾角i0,共14个参数作为稳健星历参数向量ER,k0,而将卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB作为非稳健星历参数向量即
卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB在中的顺序依次为1,2,3,4,5;
j7、将集合Mm第k0个广播星历Ek0替换为
j8、令k0=k0+1,如果k0≤nm,执行步骤j1,否则执行步骤k;
k、令m=m+1,当m≤N,执行步骤b;
m、将N个测站所有广播星历Mm,m=1,2,…,N,存入集合M中,
M={E1,E2,……,Ep},其中,nm表示第m个测站有nm组广播星历数据,p为M的元素总数;
n、将第一组广播星历数据E1的稳健参数ER,1存入ER中,非稳健参数存入中,ER,1=ER,1,
n1、令k=2,q=1;
n2、令i=1;
n3、读取M中的广播星历数据Ek;
n4、如果Ek的稳健星历数据部分ER,k与ER中的第i组稳健星历数据ER,i相同,则
令k=k+1,当k≤p时,执行步骤n2,当k>p时,执行步骤o;
n5、否则,令i=i+1,如果i≤q,执行步骤n4,如果i>q执行步骤n6;
n6、令q=q+1,将ER,k、分别作为稳健、非稳健星历参数集合ER、的第q个元素,即ER,q=ER,k,
步骤n4中ER,k与ER,i相同是指ER,k和ER,i中所有对应的广播星历数据差的绝对值均不大于相应的阈值,相应的阈值如表1所示
表1
o、构造集合以及向量Y;
o1、令为空集;
o2、令k=1;
o3、令p=1;
o4、读取集合中的第k个元素
o5、如果的行数小于3,令k=k+1,当k≤r时,其中r为中元素的个数,执行步骤o4,当k>r时,执行步骤p;
o6、否则,当的行数为偶数时,删除的最后一行;
o7、取中第1至4列的中位数以及第5列的最大值按顺序存入向量T,令其中表示矩阵的行数;
o8、如果为空集,令的第一个元素向量Y的第一个元素Y1=y,令k=k+1,当k≤r时,执行步骤o4,当k>r时,执行步骤p;
o9、当为非空集合时,令i=1;
o10、读取中的第i个元素
o11、如果的第一、第二个元素分别与T的第一个、第二个元素均相等,执行步骤o12,否则令i=i+1,当i≤p时,执行步骤o10,当i>p时,执行步骤o14;
o12、如果Y的第i个元素Yi<y,令Yi=y,执行步骤o13;
o13、令k=k+1,当k≤r时执行步骤o4,当k>r时,执行步骤p;
o14、令p=p+1,将T、y分别作为Y的第p个元素,即Yp=y,执行步骤o15;
o15、令k=k+1,当k≤r时执行步骤o4,当k>r时,执行步骤p;
p、通过上述步骤得到用于卫星导航系统完好性监测评估的广播星历数据集合它的每个元素包括卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB。
本发明的优点是:提出了一种基于频次加权识别单站用户测距精度(URA)表达形式的方法,并利用多测站用户测距精度上界值(URA_UB)合成最终的用户测距精度上界值(URA_UB),与IGS提供的brdc合成星历以及相关文献提供的合成结果相比,本方法合成的用户测距精度上界值(URA_UB)更为合理、准确,为提升卫星导航系统性能提供更科学的参考依据。
附图说明
图1是多测站广播星历确定用户测距精度上界值(URA_UB)的方法流程示意图。
具体实施方式
以下结合附图及附表,对本发明作进一步的详细说明。
实施例一
多测站广播星历确定用户测距精度上界值(URA_UB)的方法,该方法包含下列步骤:
多测站广播星历确定用户测距精度上界值的方法,该方法是在有N个可以提供广播星历数据的GNSS测站的条件下,确定用户测距精度的上界值,具体步骤如下:
a、给每个测站一个序号,第m个测站是指序号为m的测站,m=1,2,…,N;
a1、令m=1;
b、如果m>N,执行步骤m,否则读取第m个测站的广播星历数据,存入集合中,其中Mi表示该测站第i组广播星历数据,i=1,2,…,nm,nm表示第m个测站的广播星历组数;
c、将第m个测站广播星历中的用户测距精度存入向量Sm中,其中si为第i组广播星历Mi中的用户测距精度,i=1,2,…,nm;
d、构造向量S',向量S'的第一个元素为s'1,s'1=s1;
d1、令i=2,l=1;
d2、读取Sm中的第i个元素si,令j=1;
d3、读取S'中的第j个元素s'j;
d4、当si=s'j,令i=i+1,如果i>nm执行步骤e,否则执行步骤d2,当si≠s'j,令j=j+1,如果j>l,执行步骤d5,否则执行步骤d3;
d5、令l=l+1,S'的第l个元素s'l=si;
e、统计S'的元素s'i在Sm中出现的次数,记为fi,其中i=1,2,…,l,组成向量f=[f1,f2,…,fl];
f、对f进行归一化处理,得到f'
g、令L1=URA_IX=[0,1,…,15],称URA_IX为用户测距精度指数向量,
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为用户测距精度下界向量,
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为用户测距精度名义向量,
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],称URA_UB为用户测距精度上界向量;
h、辨识第m个测站采用的用户测距精度表示形式,令k=1;
h1、构造b=[b1,b2,…,bl],其中
h2、令wk=f'bT,其中bT表示b的转置向量;
h3、令k=k+1
h4、当k≤4时,返回步骤h1,否则执行步骤i;
i、令wmax=max{w1,w2,w3,w4}
i1、如果wmax<0.9,将令m=m+1,执行步骤b;
i2、否则,令k=4;
i3、如果wmax=wk,则用户测距精度采用的是Lk对应的用户测距精度表达形式,执行步骤j;
i4、如果wmax≠wk,令k=k-1,执行步骤i3;
j、将第m个测站的用户测距精度替换为上界值,令k0=1;
j1、读取Sm中的第k0个元素sk0;
j2、如果sk0>6144.00,令执行步骤j5,否则,令k1=1,执行步骤j3;
j3、如果sk0=Lk(k1),令执行步骤j5,其中Lk(k1)表示向量Lk的第k1个元素,k的取值在执行步骤j前已确定,1≤k≤4,否则,令k1=k1+1,执行步骤j4;
j4、如果k1≤15,执行步骤j3,否则,如果k≤2且k1≥16,或者k>2且k1>15,令执行步骤j5;
j5、将广播星历Ek0中的用户测距精度替换为
j6、将广播星历Ek0中6个调和改正项幅值Cuc、Cus、Crc、Crs、Cic、Cis、平动点偏差Δn、参考历元的平近点角M0、轨道偏心率e、按参考历元计算的升交点赤经Ω0、升交点赤经变化率近地点角ω、轨道半长轴平方根以及参考历元的轨道倾角i0,共14个参数作为稳健星历参数向量ER,k0,而将卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB作为非稳健星历参数向量即
卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB在中的顺序依次为1,2,3,4,5;
j7、将集合Mm第k0个广播星历Ek0替换为
j8、令k0=k0+1,如果k0≤nm,执行步骤j1,否则执行步骤k;
k、令m=m+1,当m≤N,执行步骤b;
m、将N个测站所有广播星历Mm,m=1,2,…,N,存入集合M中,
M={E1,E2,……,Ep},其中,nm表示第m个测站有nm组广播星历数据,p为M的元素总数;
n、将第一组广播星历数据E1的稳健参数ER,1存入ER中,非稳健参数存入中,ER,1=ER,1,
n1、令k=2,q=1;
n2、令i=1;
n3、读取M中的广播星历数据Ek;
n4、如果Ek的稳健星历数据部分ER,k与ER中的第i组稳健星历数据ER,i相同,则
令k=k+1,当k≤p时,执行步骤n2,当k>p时,执行步骤o;
n5、否则,令i=i+1,如果i≤q,执行步骤n4,如果i>q执行步骤n6;
n6、令q=q+1,将ER,k、分别作为稳健、非稳健星历参数集合ER、的第q个元素,即ER,q=ER,k,
步骤n4中ER,k与ER,i相同是指ER,k和ER,i中所有对应的广播星历数据差的绝对值均不大于相应的阈值,相应的阈值如表1所示
表1
o、构造集合以及向量Y;
o1、令为空集;
o2、令k=1;
o3、令p=1;
o4、读取集合中的第k个元素
o5、如果的行数小于3,令k=k+1,当k≤r时,其中r为中元素的个数,执行步骤o4,当k>r时,执行步骤p;
o6、否则,当的行数为偶数时,删除的最后一行;
o7、取中第1至4列的中位数以及第5列的最大值按顺序存入向量T,令其中表示矩阵的行数;
o8、如果为空集,令的第一个元素向量Y的第一个元素Y1=y,令k=k+1,当k≤r时,执行步骤o4,当k>r时,执行步骤p;
o9、当为非空集合时,令i=1;
o10、读取中的第i个元素
o11、如果的第一、第二个元素分别与T的第一个、第二个元素均相等,执行步骤o12,否则令i=i+1,当i≤p时,执行步骤o10,当i>p时,执行步骤o14;
o12、如果Y的第i个元素Yi<y,令Yi=y,执行步骤o13;
o13、令k=k+1,当k≤r时执行步骤o4,当k>r时,执行步骤p;
o14、令p=p+1,将T、y分别作为Y的第p个元素,即Yp=y,执行步骤o15;
o15、令k=k+1,当k≤r时执行步骤o4,当k>r时,执行步骤p;
p、通过上述步骤得到用于卫星导航系统完好性监测评估的广播星历数据集合它的每个元素包括卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB。
实施例二
多测站广播星历确定用户测距精度上界值的方法的应用。
选取国际卫星导航系统服务组织(IGS)网站上提供的2004年至2010年多天的全球跟踪站的广播星历,采用实施例一的步骤得到当天合成星历的用户测距精度上界值(URA_UB)。所选择的实验数据年份、年积日,以及当天有广播星历数据的测站数见表3。
表3
以2004年、年积日202天,即本年第202天,测站“albh”为例,执行步骤b至步骤f,得到S'={2,2.8,4,64,128},f'={0.5374,0.3878,0.0612,0.0068,0.0068},执行至步骤i,得到{w1,w2,w3,w4}={0.5986,0,1,0},wmax=1,执行至步骤j时,得到k=3,即该测站用户测距精度采用的是名义向量形式,完成步骤p后,将合成广播星历的用户测距精度上界值有{2.4,3.4,4.85,48,96,192},而国际GNSS服务组织提供的同一天的广播星历数据中的用户测距精度值有{0,1,2,2.8,4,7,8,64,128},显然采用本发明步骤合成的用户测距精度均是用户测距精度上界值向量中的元素,而国际GNSS服务组织提供的广播星历中的用户测距精度值并不全是任一Lk中的元素,k=1,2,3,4,即国际GNSS服务组织提供的广播星历的用户测距精度不统一,因而不能用于全球卫星导航系统的完好性监测与评估。为进一步说明存在的问题,提取了国际GNSS服务组织提供的、表3对应时间的广播星历中的用户测距精度数据,并提取了相同时间、采用本发明步骤合成的用户测距精度上界值,见表4。
表4
经实际检验本方法合成的广播星历用户测距精度上界值(URA_UB)与(Heng L.,Gao G.X.,Walter T.,et al.(2012).GPS Signal-in-Space Integrity PerformanceEvolution in the Last Decade[J].IEEE Transactions on Aerospace and ElectronicSystems,48(4):2932-2946)中的合成结果大部分相同,篇幅所限在此省略,仅针对部分不同的计算结果进行分析,并验证本文确定的用户测距精度上界值(URA_UB)的正确性,详情见表5,关于该表有以下几点说明:
⑴异常时刻为卫星用户测距误差(URE)开始超限的时刻;
⑵持续时间表示卫星用户测距误差(URE)超限至卫星告警之间的时间长度;
⑶URA_UB_HL一列为(Heng L.,Gao G.X.,Walter T.,et al.(2012).GPS Signal-in-Space Integrity Performance Evolution in the Last Decade[J].IEEETransactions on Aerospace and Electronic Systems,48(4):2932-2946)中合成星历中的用户测距精度上界值(URA_UB);
⑷多测站用户测距精度上界值(URA_UB)是指本专利申请方法在用户测距精度上界值(URA_UB)确定过程中获得的上界值,例如:3.40(29)+2.40(1)表示总共有29个测站的用户测距精度上界值(URA_UB)为3.40,仅有1个测站的用户测距精度上界值(URA_UB)为2.40,根据本专利申请方法确定的原则应选3.40作为用户测距精度上界值(URA_UB)的最终值;
⑸TOE表示与异常时刻对应的有效星历参考历元时刻。
从表5中可以看出:本专利申请的方法合成的用户测距精度上界值
(URA_UB)与(Heng L.,Gao G.X.,WalterT.,et al.(2012).GPS Signal-in-SpaceIntegrityPerformance Evolution in the Last Decade[J].
IEEE Transactions onAerospace and Electronic Systems,48(4):2932-2946)中的值有差异,但不难发现我们合成的用户测距精度上界值(URA_UB)是大多数测站共同观测到的,因而是较为可信的。
表5
Claims (1)
1.多测站广播星历确定用户测距精度上界值的方法,其特征在于,该方法是在有N个可以提供广播星历数据的GNSS测站的条件下,确定用户测距精度的上界值,具体步骤如下:
a、给每个测站一个序号,第m个测站是指序号为m的测站,m=1,2,…,N;
a1、令m=1;
b、如果m>N,执行步骤m,否则读取第m个测站的广播星历数据,存入集合中,其中Mi表示该测站第i组广播星历数据,i=1,2,…,nm,nm表示第m个测站的广播星历组数;
c、将第m个测站广播星历中的用户测距精度存入向量Sm中,其中si为第i组广播星历Mi中的用户测距精度,i=1,2,…,nm;
d、构造向量S',向量S'的第一个元素为s1',s1'=s1;
d1、令i=2,l=1;
d2、读取Sm中的第i个元素si,令j=1;
d3、读取S'中的第j个元素s'j;
d4、当si=s'j,令i=i+1,如果i>nm执行步骤e,否则执行步骤d2,当si≠s'j,令j=j+1,如果j>l,执行步骤d5,否则执行步骤d3;
d5、令l=l+1,S'的第l个元素s′l=si;
e、统计S'的元素s′i在Sm中出现的次数,记为fi,其中i=1,2,…,l,组成向量f=[f1,f2,…,fl];
f、对f进行归一化处理,得到f'
<mrow>
<msup>
<mi>f</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>f</mi>
<mn>1</mn>
<mo>&prime;</mo>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>f</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>f</mi>
<mi>l</mi>
<mo>&prime;</mo>
</msubsup>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>l</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
g、令L1=URA_IX=[0,1,…,15],称URA_IX为用户测距精度指数向量,
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为用户测距精度下界向量,
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为用户测距精度名义向量,
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],称URA_UB为用户测距精度上界向量;
h、辨识第m个测站采用的用户测距精度表示形式,令k=1;
h1、构造b=[b1,b2,…,bl],其中
h2、令wk=f'bT,其中bT表示b的转置向量;
h3、令k=k+1
h4、当k≤4时,返回步骤h1,否则执行步骤i;
i、令wmax=max{w1,w2,w3,w4}
i1、如果wmax<0.9,将置为空集,令m=m+1,执行步骤b;
i2、否则,令k=4;
i3、如果wmax=wk,则用户测距精度采用的是Lk对应的用户测距精度表达形式,执行步骤j;
i4、如果wmax≠wk,令k=k-1,执行步骤i3;
j、将第m个测站的用户测距精度替换为上界值,令k0=1;
j1、读取Sm中的第k0个元素sk0;
j2、如果sk0>6144.00,令执行步骤j5,否则,令k1=1,执行步骤j3;
j3、如果sk0=Lk(k1),令执行步骤j5,其中Lk(k1)表示向量Lk的第k1个元素,k的取值在执行步骤j前已确定,1≤k≤4,否则,令k1=k1+1,执行步骤j4;
j4、如果k1≤15,执行步骤j3,否则,如果k≤2且k1≥16,或者k>2且k1>15,令执行步骤j5;
j5、将广播星历Ek0中的用户测距精度替换为
j6、将广播星历Ek0中6个调和改正项幅值Cuc、Cus、Crc、Crs、Cic、Cis、平动点偏差Δn、参考历元的平近点角M0、轨道偏心率e、按参考历元计算的升交点赤经Ω0、升交点赤经变化率近地点角ω、轨道半长轴平方根以及参考历元的轨道倾角i0,共14个参数作为稳健星历参数向量ER,k0,而将卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB作为非稳健星历参数向量即
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>R</mi>
<mo>,</mo>
<mi>k</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>u</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>u</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>r</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>&Delta;</mi>
<mi>n</mi>
<mo>,</mo>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mi>e</mi>
<mo>,</mo>
<msub>
<mi>&Omega;</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mover>
<mi>&Omega;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>,</mo>
<mi>&omega;</mi>
<msqrt>
<mrow>
<mi>A</mi>
<mo>,</mo>
</mrow>
</msqrt>
<msub>
<mi>i</mi>
<mn>0</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>k</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>P</mi>
<mi>R</mi>
<mi>N</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>o</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
<mi>t</mi>
<mi>h</mi>
<mo>_</mo>
<mi>i</mi>
<mi>f</mi>
<mo>,</mo>
<mi>U</mi>
<mi>R</mi>
<mi>A</mi>
<mo>_</mo>
<mi>U</mi>
<mi>B</mi>
<mo>&rsqb;</mo>
</mrow>
卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB在ER,k0中的顺序依次为1,2,3,4,5;
j7、将集合Mm第k0个广播星历Ek0替换为
j8、令k0=k0+1,如果k0≤nm,执行步骤j1,否则执行步骤k;
k、令m=m+1,当m≤N,执行步骤b;
m、将N个测站所有广播星历Mm,m=1,2,…,N,存入集合M中,
M={E1,E2,……,Ep},其中,nm表示第m个测站有nm组广播星历数据,p为M的元素总数;
n、将第一组广播星历数据E1的稳健参数ER,1存入ER中,非稳健参数存入中,ER,1=ER,1,
n1、令k=2,q=1;
n2、令i=1;
n3、读取M中的广播星历数据Ek;
n4、如果Ek的稳健星历数据部分ER,k与ER中的第i组稳健星历数据ER,i相同,则
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>E</mi>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>E</mi>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
2
令k=k+1,当k≤p时,执行步骤n2,当k>p时,执行步骤o;
n5、否则,令i=i+1,如果i≤q,执行步骤n4,如果i>q执行步骤n6;
n6、令q=q+1,将ER,k、分别作为稳健、非稳健星历参数集合ER、的第q个元素,即ER,q=ER,k,
步骤n4中ER,k与ER,i相同是指ER,k和ER,i中所有对应的广播星历数据差的绝对值均不大于相应的阈值,相应的阈值如表1所示
表1
o、构造集合以及向量Y;
o1、令为空集;
o2、令k=1;
o3、令p=1;
o4、读取集合中的第k个元素
o5、如果的行数小于3,令k=k+1,当k≤r时,其中r为中元素的个数,执行步骤o4,当k>r时,执行步骤p;
o6、否则,当的行数为偶数时,删除的最后一行;
o7、取中第1至4列的中位数以及第5列的最大值按顺序存入向量T,令其中表示矩阵的行数;
o8、如果为空集,令的第一个元素向量Y的第一个元素Y1=y,令k=k+1,当k≤r时,执行步骤o4,当k>r时,执行步骤p;
o9、当为非空集合时,令i=1;
o10、读取中的第i个元素
o11、如果的第一、第二个元素分别与T的第一个、第二个元素均相等,执行步骤o12,否则令i=i+1,当i≤p时,执行步骤o10,当i>p时,执行步骤o14;
o12、如果Y的第i个元素Yi<y,令Yi=y,执行步骤o13;
o13、令k=k+1,当k≤r时执行步骤o4,当k>r时,执行步骤p;
o14、令p=p+1,将T、y分别作为Y的第p个元素,即Yp=y,执行步骤o15;
o15、令k=k+1,当k≤r时执行步骤o4,当k>r时,执行步骤p;
p、通过上述步骤得到用于卫星导航系统完好性监测评估的广播星历数据集合它的每个元素包括卫星编号PRN、星历参考历元toe、卫星钟差参考历元toc、卫星健康状况标记health_if、用户测距精度上界值URA_UB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510966601.2A CN105467406B (zh) | 2015-12-18 | 2015-12-18 | 多测站广播星历确定用户测距精度上界值的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510966601.2A CN105467406B (zh) | 2015-12-18 | 2015-12-18 | 多测站广播星历确定用户测距精度上界值的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105467406A CN105467406A (zh) | 2016-04-06 |
CN105467406B true CN105467406B (zh) | 2017-10-24 |
Family
ID=55605304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510966601.2A Active CN105467406B (zh) | 2015-12-18 | 2015-12-18 | 多测站广播星历确定用户测距精度上界值的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105467406B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106873009B (zh) * | 2017-01-12 | 2019-08-02 | 中国科学院测量与地球物理研究所 | 利用非组合ppp辅助的长基线模糊度解算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010053395A1 (en) * | 2008-11-05 | 2010-05-14 | Joint Stock Company "Russian Space Systems" (Jsc "Rss") | Method and device of satellite navigation system's integrity control |
EP2244097B1 (en) * | 2009-04-20 | 2014-02-26 | Astrium GmbH | A method of reducing the impact of space vehicle frequency jumps on a global navigation device |
CN104215974B (zh) * | 2014-09-15 | 2017-01-11 | 中国航天标准化研究所 | 一种卫星导航系统的完好性监测可用性确定方法 |
CN104614734B (zh) * | 2014-12-22 | 2017-03-29 | 中国科学院国家授时中心 | 北斗天基高精度差分电文生成方法 |
CN104749586A (zh) * | 2015-03-27 | 2015-07-01 | 中国电子科技集团公司第二十研究所 | 一种卫星导航伪距异常实时检测方法 |
-
2015
- 2015-12-18 CN CN201510966601.2A patent/CN105467406B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105467406A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102353936B (zh) | 带有时钟差建模的多gnss信号处理 | |
CN109313268A (zh) | 具有宽巷偏差校正值和窄巷偏差校正值的导航卫星的轨道和低延时时钟的确定 | |
England et al. | Late Quaternary to decadal velocity fields in Asia | |
CN104502935A (zh) | 一种基于非差非组合模型的网络rtk模糊度解算方法 | |
CN109154664A (zh) | 具有低延时时钟校正值的导航卫星轨道和时钟的确定 | |
CN103728643B (zh) | 附有宽巷约束的北斗三频网络rtk模糊度单历元固定方法 | |
CN109219762A (zh) | 导航卫星宽巷偏差的确定和超范围调整系统和方法 | |
CN103675835A (zh) | 一种北斗三频信号载波相位整周模糊度单历元确定方法 | |
Lee | New Advanced RAIM with Improved Availability for Detecting Constellation‐wide Faults, Using Two Independent Constellations | |
Li et al. | Undifferenced cycle slip estimation of triple-frequency BeiDou signals with ionosphere prediction | |
CN105158778B (zh) | 多系统联合实施载波相位差分故障卫星剔除方法及其系统 | |
CN108196284A (zh) | 一种进行星间单差模糊度固定的gnss网数据处理方法 | |
CN108549095A (zh) | 一种区域cors网非差并行增强方法及系统 | |
Dimitrova et al. | High-resolution view of active tectonic deformation along the Hikurangi subduction margin and the Taupo Volcanic Zone, New Zealand | |
CN106017475A (zh) | 一种航迹更新方法及装置 | |
CN105842720A (zh) | 一种基于载波相位的广域精密实时定位方法 | |
CN110398758A (zh) | 实时钟差估计中的粗差探测方法、装置、设备及存储介质 | |
CN107219537A (zh) | 一种融合选星与完好性检测的多系统兼容定位方法 | |
CN105511481B (zh) | 一种星载定轨优化方法 | |
Mohanty et al. | Learning GNSS positioning corrections for smartphones using graph convolution neural networks | |
CN105467406B (zh) | 多测站广播星历确定用户测距精度上界值的方法 | |
CN108873023A (zh) | 一种提高定位精度的观测量的处理方法 | |
CN111142125A (zh) | 一种卫星完好性监测方法和系统 | |
Choi et al. | Jason-1 precision orbit determination by combining SLR and DORIS with GPS tracking data | |
CN107194161A (zh) | 基于用户需求分级的araim可用性预测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |