CN111679294B - 一种星基增强系统完好性数据的统计检验方法 - Google Patents

一种星基增强系统完好性数据的统计检验方法 Download PDF

Info

Publication number
CN111679294B
CN111679294B CN202010723220.2A CN202010723220A CN111679294B CN 111679294 B CN111679294 B CN 111679294B CN 202010723220 A CN202010723220 A CN 202010723220A CN 111679294 B CN111679294 B CN 111679294B
Authority
CN
China
Prior art keywords
udre
data
satellite
satellite broadcast
test
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
Application number
CN202010723220.2A
Other languages
English (en)
Other versions
CN111679294A (zh
Inventor
尹继东
李烁星
何宁
黄智�
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Academy Of Aerospace Science Technology And Communications Technology Co ltd
Original Assignee
Academy Of Aerospace Science Technology And Communications Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Academy Of Aerospace Science Technology And Communications Technology Co ltd filed Critical Academy Of Aerospace Science Technology And Communications Technology Co ltd
Priority to CN202010723220.2A priority Critical patent/CN111679294B/zh
Publication of CN111679294A publication Critical patent/CN111679294A/zh
Application granted granted Critical
Publication of CN111679294B publication Critical patent/CN111679294B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • 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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/08—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing integrity information, e.g. health of satellites or quality of ephemeris data
    • 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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/10—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals
    • G01S19/11—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals wherein the cooperating elements are pseudolites or satellite radio beacon positioning system signal repeaters
    • 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 Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Security & Cryptography (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种星基增强系统完好性数据的统计检验方法，包括以下步骤：S1、通过监测站观测数据获取卫星导航电文信息和SBAS差分修正信息；S2、通过卫星导航电文信息和SBAS差分修正信息，经伪距误差修正得到实测UDRE数据；S3、从监测站观测数据中解析出卫星广播UDRE′数据；S4、对实测UDRE与卫星广播UDRE′数据进行方差齐性检验；S5、判定星基增强系统伪距域的完好性数据是否可信，系统工作状态是否出现异常。与现有技术相比，本发明提供的统计检验方法，基于假设检验统计思想，从卫星端的角度进行统计检验，能直接检测出故障卫星，对多故障情况的检测能力较快。

Description

一种星基增强系统完好性数据的统计检验方法
技术领域
本发明涉及星基增强技术领域，具体是一种星基增强系统完好性数据的统计检验方法。
背景技术
由于卫星导航系统的定位精度并不能满足其在某些应用领域的需求，因此各航天大国争相发展各自的星基增强系统(SBAS)。星基增强系统是高精度卫星导航系统的重要组成部分，它通过观测站计算有关GNSS的完好性和修正数据，并借助地球静止轨道(GEO)卫星上的透明转发器向用户播发完好性信息和修正信息(包括星历误差、卫星钟差、电离层延迟等)给SBAS用户，同时提供测距信号增强卫星信号，从而来提高原有导航系统的定位精度及各项性能。完好性参数包括用户差分伪距误差(UDRE)和电离层格网点垂直改正误差(GIVE)，其中GIVE与电离层延迟有关，更新周期为120s；UDRE与卫星星历与钟差改正有关，更新周期为6s。从数据更新周期来看，UDRE对定位误差的影响更敏感，所以SBAS的完好性主要通过UDRE的检验结果来反应。完好性反应了通过观测值解算出定位结果的可靠性，在用户安全性方面，有着比精度更重要的地位。目前，SBAS用户的完好性主要通过计算用户端位置域定位误差，通过完好性阈值与告警限值作比较来判定，但是位置域检验方法存在计算量大、多卫星故障检测能力差等缺点。
发明内容
本发明的目的在于克服现有技术通过计算用户端位置域定位误差，通过完好性阈值与告警限值作比较来判定，但是位置域检验方法存在计算量大、多卫星故障检测能力差等缺点的不足，提供了一种星基增强系统完好性数据的统计检验方法，基于假设检验统计思想，从卫星端的角度出发对SBAS完好性参数进行统计检验，通过单颗卫星的伪距域UDRE参数检验结果，判定SBAS完好性数据是否可信，系统工作状态是否正常，能直接检测出故障卫星，相比基于位置域的完好性数据检验方法，其对多故障情况的检测能力较快。
本发明的目的主要通过以下技术方案实现：
一种星基增强系统完好性数据的统计检验方法，包括以下步骤：S1、通过监测站观测数据获取卫星导航电文信息和SBAS差分修正信息；S2、通过卫星导航电文信息和SBAS差分修正信息，经伪距误差修正得到实测UDRE数据；S3、从监测站观测数据中获取卫星广播UDRE′数据；S4、对实测UDRE数据与卫星广播UDRE′数据进行方差齐性检验；S5、若实测UDRE数据与卫星广播UDRE′数据方差检验结果齐性，则对实测UDRE数据与卫星广播UDRE′数据的均值差进行T检验，并判断实测UDRE数据和卫星广播UDRE′数据的误差是否在卫星广播误差估计值限值内，若是则判定星基增强系统伪距域的完好性数据可信，反之，则认为完好性数据不可信，系统工作出现异常；若实测UDRE数据与卫星广播UDRE′数据方差检验结果非齐性，则使用校正T检验对实测UDRE数据与卫星广播UDRE′数据的均值差进行检验，并判断实测UDRE数据和卫星广播UDRE′数据的误差是否在卫星广播误差估计值限值内，若是则判定星基增强系统伪距域的完好性数据可信，反之，则认为完好性数据不可信，系统工作出现异常。
现有技术中星基增强系统用户的完好性主要通过计算用户端位置域定位误差，通过完好性阈值与告警限值作比较来判定，但是位置域检验方法存在计算量大、多卫星故障检测能力差等缺点。针对现有技术从通过计算用户端位置域定位误差的不足之处，本技术方案首先通过监测站观测数据获取GPS电文与SBAS差分校正信息，再经伪距误差修正计算得到实测UDRE，同时直接解析SBAS广播UDRE′参数，其次基于统计检验方法，先对监测站实测UDRE与卫星广播UDRE′数据进行方差齐性检验，在方差检验结果齐性时再进行T检验，方差检验结果非齐性时进行校正T检验，最后在既定显著性水平及置信区间内判定星基增强系统伪距域的完好性，通过T检验联合方差齐性检验，使数据更加满足T检验需求，提高检验结果的可靠性；本技术方案通过监测站观测数据获取GPS电文与SBAS差分校正信息，通过差分校正的方法可以降低甚至消除某些误差，从而提高导航性能，通过导航电文信息对伪距观测值进行伪距误差修正得到实测UDRE数据，能够使得到的实测UDRE与卫星广播UDRE′为两组相互独立的观察值，通过将实测UDRE与卫星广播UDRE′进行方差齐性检验确定实测UDRE与卫星广播UDRE′数据总体平均数的差异，以判断实测UDRE与卫星广播UDRE′的差异，再通过T检验对实测UDRE与卫星广播UDRE′进行分析，判断实测UDRE与卫星广播UDRE′的符合性，得到伪距误差的可信度，以判定星基增强系统伪距域的完好性数据是否可信，从而判断卫星完好性，相比于位置域的完好性检验，对多故障情况有快速的检测能力，本技术方案中的星基增强系统完好性数据指星基增强系统伪距域的完好性数据。本技术方案基于假设检验统计思想，从卫星端的角度出发对SBAS完好性参数进行统计检验，通过单颗卫星的伪距域UDRE参数检验结果，判定SBAS完好性数据是否可信，系统工作状态是否正常，能直接检测出故障卫星，相比基于位置域的完好性数据检验方法，能对所有可见卫星逐颗进行检验，其对多故障情况的检测能力较快。本技术方案中的卫星导航电文信息包括卫星星历、钟差修正、电离层延迟修正等信息，SBAS差分修正信息包括快钟修正、完好性信息及长期钟差修正等。
进一步的，S2包括如下步骤：S2.1、通过卫星导航电文信息对伪距观测值进行伪距误差修正，得到修正后的伪距观测量；S2.2、对修正后的伪距观测量进行SBAS差分修正，得到修正后的伪距值；S2.3、用修正后的伪距值与真距值做差计算得到实测UDRE的值。
本技术方案提供了通过导航电文信息对伪距观测值进行伪距误差修正，得到实测UDRE数据的具体方法，通过对修正后的伪距观测量进行SBAS差分修正得到修正后的伪距值，通过差分校正的方法可以降低甚至消除某些误差，从而提高导航性能，通过导航电文信息对伪距观测值进行伪距误差修正得到实测UDRE数据，能够使得到的实测UDRE与卫星广播UDRE′为两组相互独立的观察值，提高检验结果的可靠性。需要说明的是，本技术方案采用的SBAS差分修正采用现有技术。
进一步的，S2.1中采用的伪距误差修正公式为：P＝ρ-Δt(s)+Δtu-dr+dion+dtrop+εp，其中，P为修正后的伪距观测值，ρ为卫星与用户接收机的几何距离，Δt(s)为广播星历引起的伪距修正值，Δtu为监测站接收机钟差引起的伪距修正值，dr为相对论效应伪距修正值，dion为电离层误差修正值，dtrop为对流层误差修正值，εp为GPS误差修正后的残差。
本技术方案提供了伪距误差修正的具体方法，其中伪距误差修正公式中的参数Δt(s)、dr、dion、dtrop可通过导航电文解析获取，参数ρ通过确定用户接收机的三维坐标信息直接获得，Δtu的值可通过用户接收机获取，伪距误差修正的计算过程采用现有的方法。
进一步的，S4中方差齐性检验包括如下步骤：S4.1、设X1,X2,…,Xn是来自实测UDRE数据X的一个样本，Y1,Y2,…,Yn是来自广播UDER′数据Y的一个样本，X与Y相互独立，且
Figure BDA0002600767160000031
两个样本的均值分别为
Figure BDA0002600767160000032
与
Figure BDA0002600767160000033
样本方差分别为
Figure BDA0002600767160000034
与
Figure BDA0002600767160000035
参数μ1,μ2,
Figure BDA0002600767160000036
均未知；S4.2、方差齐性检验假设：
Figure BDA0002600767160000037
即认为实测UDRE与卫星广播UDRE′方差无显著性差异；
Figure BDA0002600767160000038
即认为实测UDRE与卫星广播UDRE′方差存在显著性差异；S4.3、当假设H0成立时，检验统计量为：
Figure BDA0002600767160000039
n1和n2为给定的样本容量；S4.4、在给定显著性水平α、样本容量n1与n2时，根据样本值计算s1 2,s2 2的值，与检验统计量的观测值F作以下判断：若F≥Fα/2(n1-1,n2-1)或F≤F1-α/2(n1-1,n2-1)，则拒绝H0；若F1-α/2(n1-1,n2-1)＜F＜Fα/2(n1-1,n2-1)，则接受H0。
本技术方案提供了对实测UDRE与卫星广播UDRE′数据进行方差齐性检验的具体方法，通过方差齐性检验假设，并对假设进行检验，来确定实测UDRE和广播UDER′的值是否满足方差齐性，通过T检验联合方差齐性检验，使数据更加满足T检验需求，提高检验结果的可靠性。
进一步的，当实测UDRE与卫星广播UDRE′数据方差检验结果齐性时，对实测UDRE与卫星广播UDRE′数据的均值差进行T检验的方法如下：S5.1、设X1,X2,…,Xn是来自实测UDRE数据X的一个样本，Y1,Y2,…,Yn是来自广播UDER′数据Y的一个样本，X与Y相互独立，且
Figure BDA0002600767160000041
两个样本的均值分别为
Figure BDA0002600767160000042
与
Figure BDA0002600767160000043
样本方差分别为
Figure BDA0002600767160000044
与
Figure BDA0002600767160000045
参数μ1,μ2,
Figure BDA0002600767160000046
均未知；S5.2、检验假设：H3：μ1-μ2＝δ≤0，即认为实测UDRE的值比卫星广播UDRE′的值要低，即真实误差在卫星广播误差估计值限值内，卫星播发的校正参数可信；H4：μ1-μ2＝δ＞0，即认为实测UDRE的值比卫星广播UDRE′的值高，卫星播发的修正参数不可信；S5.3、在H3为真时，T统计量为：
Figure BDA0002600767160000047
其中，
Figure BDA0002600767160000048
n1和n2为给定的样本容量；S5.4、在给定显著性水平α、样本容量n1与n2时，根据样本值计算
Figure BDA0002600767160000049
sW的值，对检验统计量的观测值作以下判断：若t≥ta(n1+n2-2)，则拒绝H3，即星基增强系统伪距域的完好性数据不可信；若t＜tα(n1+n2-2)，则接受H3，即星基增强系统伪距域的完好性数据可信。
本技术方案提供了当实测UDRE与卫星广播UDRE′数据方差检验结果齐性时，对实测UDRE与卫星广播UDRE′数据的均值差进行T检验方法，通过检验统计量的观测值进行判断，确定星基增强系统伪距域的完好性数据是否可信，相比于位置域的完好性检验，对多故障情况有快速的检测能力。需要说明的是，本技术方案中的参数μ1,
Figure BDA00026007671600000410
是实测UDRE数据X服从正态分布的均值和方差，μ2,
Figure BDA00026007671600000411
是卫星广播UDRE′数据Y服从正态分布的均值和方差，为总体指标，统计思想用样本数据指标检验总体指标。
进一步的，当实测UDRE数据与卫星广播UDRE′数据的方差检验结果非齐性时，使用校正T检验对实测UDRE数据与卫星广播UDRE′数据的均值差进行检验的方法为：对实测UDRE数据与卫星广播UDRE′数据进行校正，对校正后的实测UDRE与卫星广播UDRE′数据进行方差齐性检验，若校正后的实测UDRE与卫星广播UDRE′数据方差检验结果齐性，则进行T检验；若校正后实测UDRE与卫星广播UDRE′数据的方差检验结果非齐性，则调整给定显著性水平α，再进行T检验。
本技术方案中当实测UDRE与卫星广播UDRE′数据方差检验结果非齐性时，使用校正T检验对实测UDRE与卫星广播UDRE′数据的均值差进行检验，通过对实测UDRE与卫星广播UDRE′的数据进行校正，若校正后的数据方差检验结果仍为非齐性，则调整给定显著性水平α，再进行T检验。本技术方案中通过调整给定显著性水平α，通常设定显著性水平α取值为0.05，置信度为95％，可通过调整显著性水平α，使置信度提高，以使获得的星基增强系统伪距域的完好性数据判定结果可靠。优选的，在调整给定显著性水平α后，对校正后的实测UDRE与卫星广播UDRE′数据进行方差齐性检验，若校正后的实测UDRE与卫星广播UDRE′数据方差检验结果齐性，则进行T检验，若检验结果为非齐性，则继续调整给定显著性水平α直至校正后的实测UDRE与卫星广播UDRE′方差检验结果齐性或近似于齐性。
进一步的，采用处理数据缺失值和剔除数据的异常值对实测UDRE与卫星广播UDRE′数据进行校正。
本技术方案通过处理数据缺失值和剔除数据的异常值校正数据，处理数据缺失值包括删除存在缺失值的个案和缺失值插补，可以采用现有的简单删除法和权重法删除含有缺失值的个案，采用现有的均值插补、利用同类均值插补、极大似然估计或多重插补，通过处理数据缺失值能够有效减少由于缺失数据造成的估计量偏倚，以使校正后的实测UDRE与卫星广播UDRE′数据方差检验结果齐性；通过剔除数据的异常值，将数据中存在错误或重复的数据进行剔除，可以采用现有的图形法，如箱线图、正态分布图，和建模法，如线性回归、聚类算法、K近邻算法，以使校正后的实测UDRE与卫星广播UDRE′数据方差检验结果齐性或近似于齐性。
综上所述，本发明与现有技术相比具有以下有益效果：本发明通过监测站观测数据获取GPS电文与SBAS差分校正信息，通过差分校正的方法可以降低甚至消除某些误差，从而提高导航性能，通过导航电文信息对伪距观测值进行伪距误差修正得到实测UDRE数据，能够使得到的实测UDRE与卫星广播UDRE′为两组相互独立的观察值，通过将实测UDRE与卫星广播UDRE′进行方差检验结果齐性检验确定实测UDRE与卫星广播UDRE′的总体平均数的差异，以判断实测UDRE与卫星广播UDRE′的差异，再通过T检验对实测UDRE与卫星广播UDRE′进行分析，判断实测UDRE与卫星广播UDRE′的符合性，得到伪距误差的可信度，以判定星基增强系统伪距域的完好性数据是否可信，从而判断卫星完好性，相比于位置域的完好性检验，对多故障情况有快速的检测能力。本发明基于假设检验统计思想，从卫星端的角度出发对SBAS完好性参数进行统计检验，通过单颗卫星的伪距域UDRE参数检验结果，判定SBAS完好性数据是否可信，系统工作状态是否正常，能直接检测出故障卫星，相比基于位置域的完好性数据检验方法，能对所有可见卫星逐颗进行检验，其对多故障情况的检测能力较快。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为一种星基增强系统完好性数据的统计检验方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1：
如图1所示，本实施例包括以下步骤：
S1、通过监测站观测数据获取卫星导航电文信息和SBAS差分修正信息；S2、通过卫星导航电文信息和SBAS差分修正信息，通过伪距误差修正得到实测UDRE数据；
S3、从监测站观测数据中获取卫星广播UDRE′数据；
S4、对实测UDRE数据与卫星广播UDRE′数据进行方差齐性检验；
S5、若实测UDRE数据与卫星广播UDRE′数据方差检验结果齐性，则对实测UDRE数据与卫星广播UDRE′数据的均值差进行T检验，并判断实测UDRE数据和卫星广播UDRE′数据的误差是否在卫星广播误差估计值限值内，若是则判定星基增强系统伪距域的完好性数据可信，反之，则认为完好性数据不可信，系统工作出现异常；若实测UDRE数据与卫星广播UDRE′数据方差检验结果非齐性，则使用校正T检验对实测UDRE数据与卫星广播UDRE′数据的均值差进行检验，并判断实测UDRE数据和卫星广播UDRE′数据的误差是否在卫星广播误差估计值限值内，若是则判定星基增强系统伪距域的完好性数据可信，反之，则认为完好性数据不可信，系统工作出现异常。
优选的，S2包括如下步骤：
S2.1、通过卫星导航电文信息对伪距观测值进行伪距误差修正，得到修正后的伪距观测量；
S2.2、对修正后的伪距观测量进行SBAS差分修正，得到修正后的伪距值；
S2.3、用修正后的伪距值与真距值做差计算得到实测UDRE的值。
优选的，S2.1中采用的伪距误差修正公式为：P＝ρ-Δt(s)+Δtu-dr+dion+dtrop+εp，其中，P为修正后的伪距观测值，ρ为卫星与用户接收机的几何距离，Δt(s)为广播星历引起的伪距修正值，Δtu为监测站接收机钟差引起的伪距修正值，dr为相对论效应伪距修正值，dion为电离层误差修正值，dtrop为对流层误差修正值，εp为GPS误差修正后的残差。
优选的，S4中方差齐性检验包括如下步骤：
S4.1、设X1,X2,…,Xn是来自实测UDRE数据X的一个样本，Y1,Y2,…,Yn是来自广播UDER′数据Y的一个样本，X与Y相互独立，且
Figure BDA0002600767160000071
两个样本的均值分别为
Figure BDA0002600767160000072
与
Figure BDA0002600767160000073
样本方差分别为
Figure BDA0002600767160000074
与
Figure BDA0002600767160000075
参数μ1,μ2,
Figure BDA0002600767160000076
均未知；
S4.2、方差齐性检验假设：
Figure BDA0002600767160000077
即认为实测UDRE与卫星广播UDRE′方差无显著性差异；
Figure BDA0002600767160000078
即认为实测UDRE与卫星广播UDRE′方差存在显著性差异；
S4.3、当假设H0成立时，检验统计量为：
Figure BDA0002600767160000079
n1和n2为给定的样本容量；
S4.4、在给定显著性水平α、样本容量n1与n2时，根据样本值计算
Figure BDA00026007671600000710
的值，与检验统计量的观测值F作以下判断：若F≥Fα/2(n1-1,n2-1)或F≤F1-α/2(n1-1,n2-1)，则拒绝H0；若F1-α/2(n1-1，n2-1)＜F＜Fα/2(n1-1，n2-1)，则接受H0。
优选的，当实测UDRE与卫星广播UDRE′数据方差检验结果齐性时，对实测UDRE与卫星广播UDRE′数据的均值差进行T检验的方法如下：
S5.1、设X1,X2,…,Xn是来自实测UDRE数据X的一个样本，Y1,Y2,…,Yn是来自广播UDER′数据Y的一个样本，X与Y相互独立，且
Figure BDA00026007671600000711
两个样本的均值分别为
Figure BDA00026007671600000712
与
Figure BDA00026007671600000713
样本方差分别为
Figure BDA00026007671600000714
与
Figure BDA00026007671600000715
参数
Figure BDA00026007671600000716
均未知；
S5.2、检验假设：H3：μ1-μ2＝δ≤0，即认为实测UDRE的值比卫星广播UDRE′的值要低，即真实误差在卫星广播误差估计值限值内，卫星播发的校正参数可信；H4：μ1-μ2＝δ＞0，即认为实测UDRE的值比卫星广播UDRE′的值高，卫星播发的校正参数不可信；
S5.3、在H3为真时，T统计量为：
Figure BDA0002600767160000081
其中，
Figure BDA0002600767160000082
n1和n2为给定的样本容量；
S5.4、在给定显著性水平α、样本容量n1与n2时，根据样本值计算
Figure BDA0002600767160000083
的值，对检验统计量的观测值作以下判断：若t≥tα(n1+n2-2)，则拒绝H3，即星基增强系统伪距域的完好性数据不可信；若t＜tα(n1+n2-2)，则接受H3，即星基增强系统伪距域的完好性数据可信。
实施例2：
本实施例在实施例1的基础上，还包括：当实测UDRE与卫星广播UDRE′方差检验结果非齐性时，使用校正T检验对实测UDRE与卫星广播UDRE′数据的均值差进行检验的方法为：对实测UDRE与卫星广播UDRE′的数据进行校正，对校正后的实测UDRE与卫星广播UDRE′数据进行方差齐性检验，若校正后的实测UDRE与卫星广播UDRE′数据方差检验结果齐性，则进行T检验；若校正后实测UDRE的与卫星广播UDRE′数据方差检验结果非齐性，则调整给定显著性水平α，再进行T检验。优选的，在调整给定显著性水平α后，对校正后的实测UDRE与卫星广播UDRE′进行方差齐性检验，若校正后的实测UDRE与卫星广播UDRE′数据方差检验结果齐性，则进行T检验，若方差检验结果为非齐性，则继续调整给定显著性水平α的值直至校正后的实测UDRE与卫星广播UDRE′方差检验结果齐性或近似于齐性。所述T检验与实施例1中当实测UDRE与卫星广播UDRE′方差齐性时采用的T检验方法相同。
优选的，采用处理数据缺失值和剔除数据的异常值对实测UDRE与卫星广播UDRE′的数据进行校正。
优选的，对校正后的实测UDRE与卫星广播UDRE′进行方差齐性检验包括如下步骤：
(1)设X1,X2,…,Xn是来自实测UDRE数据X的一个样本，Y1,Y2,…,Yn是来自广播UDER′数据Y的一个样本，X与Y相互独立，且
Figure BDA0002600767160000084
两个样本的均值分别为
Figure BDA0002600767160000085
与
Figure BDA0002600767160000086
样本方差分别为
Figure BDA0002600767160000087
与
Figure BDA0002600767160000088
参数μ1,μ2,
Figure BDA0002600767160000089
均未知；
(2)方差齐性检验假设：
Figure BDA00026007671600000810
即认为实测UDRE与卫星广播UDRE′方差无显著性差异；
Figure BDA0002600767160000091
即认为实测UDRE与卫星广播UDRE′方差存在显著性差异；
(3)当假设H0成立时，检验统计量为：
Figure BDA0002600767160000092
n1和n2为给定的样本容量；
(4)在给定显著性水平α、样本容量n1与n2时，根据样本值计算
Figure BDA0002600767160000093
的值，与检验统计量的观测值F作以下判断：若F≥Fα/2(n1-1,n2-1)或F≤F1-α/2(n1-1,n2-1)，则拒绝H0；若F1-α/2(n1-1，n2-1)＜F＜Fα/2(n1-1，n2-1)，则接受H0。
本实施例基于假设检验统计思想，从卫星端的角度出发对SBAS完好性参数进行统计检验，通过单颗卫星的伪距域UDRE参数检验结果，判定SBAS完好性数据是否可信，系统工作状态是否正常，能直接检测出故障卫星，相比基于位置域的完好性数据检验方法，能对所有可见卫星逐颗进行检验，其对多故障情况的检测能力较快。本实施例通过监测站观测数据获取GPS电文与SBAS差分校正信息，通过差分校正的方法可以降低甚至消除某些误差，从而提高导航性能，通过导航电文信息对伪距观测值进行伪距误差修正得到实测UDRE数据，能够使得到的实测UDRE与卫星广播UDRE′为两组相互独立的观察值，通过将实测UDRE与卫星广播UDRE′进行方差齐性检验确定实测UDRE与卫星广播UDRE′的总体平均数的差异，以判断实测UDRE与卫星广播UDRE′的差异，再通过T检验对实测UDRE与卫星广播UDRE′进行分析，判断实测UDRE与卫星广播UDRE′的符合性，得到伪距误差的可信度，以判定星基增强系统伪距域的完好性数据是否可信，从而判断卫星完好性，相比于位置域的完好性检验，对多故障情况有快速的检测能力。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (4)

1.一种星基增强系统完好性数据的统计检验方法，其特征在于，包括以下步骤：
S1、通过监测站观测数据获取卫星导航电文信息和SBAS差分修正信息；
S2、通过卫星导航电文信息和SBAS差分修正信息，经伪距误差修正得到实测UDRE数据；
S3、从监测站观测数据中获取卫星广播UDRE′数据；
S4、对实测UDRE数据与卫星广播UDRE′数据进行方差齐性检验；
S5、若实测UDRE数据与卫星广播UDRE′数据方差检验结果齐性，则对实测UDRE数据与卫星广播UDRE′数据的均值差进行T检验，并判断实测UDRE数据和卫星广播UDRE′数据的误差是否在卫星广播误差估计值限值内，若是则判定星基增强系统伪距域的完好性数据可信，反之，则认为完好性数据不可信，系统工作出现异常；若实测UDRE数据与卫星广播UDRE′数据方差检验结果非齐性，则使用校正T检验对实测UDRE数据与卫星广播UDRE′数据的均值差进行检验，并判断实测UDRE数据和卫星广播UDRE′数据的误差是否在卫星广播误差估计值限值内，若是则判定星基增强系统伪距域的完好性数据可信，反之，则认为完好性数据不可信，系统工作出现异常；
S4中方差齐性检验包括如下步骤：
S4.1、设X1,X2,…,Xn是来自实测UDRE数据X的一个样本，Y1,Y2,…,Yn是来自广播UDER′数据Y的一个样本，X与Y相互独立，且
Figure FDA0002929682080000011
两个样本的均值分别为
Figure FDA0002929682080000012
与
Figure FDA0002929682080000013
样本方差分别为
Figure FDA0002929682080000014
与
Figure FDA0002929682080000015
参数μ1,μ2,
Figure FDA0002929682080000016
均未知；
S4.2、方差齐性检验假设：
Figure FDA0002929682080000017
即认为实测UDRE数据与卫星广播UDRE′数据的方差无显著性差异；
Figure FDA0002929682080000018
即认为实测UDRE数据与卫星广播UDRE′数据的方差存在显著性差异；
S4.3、当假设H0成立时，检验统计量为：
Figure FDA0002929682080000019
n1和n2为给定的样本容量；
S4.4、在给定显著性水平α、样本容量n1与n2时，根据样本值计算
Figure FDA00029296820800000110
的值，与检验统计量的观测值F作以下判断：若F≥Fα/2(n1-1,n2-1)或F≤F1-α/2(n1-1,n2-1)，则拒绝H0；若F1-α/2(n1-1，n2-1)＜F＜Fα/2(n1-1，n2-1)，则接受H0；
当实测UDRE数据与卫星广播UDRE′数据的方差检验结果齐性时，对实测UDRE数据与卫星广播UDRE′数据的均值差进行T检验的方法如下：
S5.1、设X1,X2,…,Xn是来自实测UDRE数据X的一个样本，Y1,Y2,…,Yn是来自广播UDER′数据Y的一个样本，X与Y相互独立，且
Figure FDA0002929682080000021
两个样本的均值分别为
Figure FDA0002929682080000022
与
Figure FDA0002929682080000023
样本方差分别为
Figure FDA0002929682080000024
与
Figure FDA0002929682080000025
参数μ1,μ2,
Figure FDA0002929682080000026
均未知；
S5.2、检验假设：H3：μ1-μ2＝δ≤0，即认为实测UDRE的值比卫星广播UDRE′的值要低，即真实误差在卫星广播误差估计值限值内，卫星播发的校正参数可信；H4：μ1-μ2＝δ＞0，即认为实测UDRE的值比卫星广播UDRE′的值高，卫星播发的校正参数不可信；
S5.3、在H3为真时，T统计量为：
Figure FDA0002929682080000027
其中，
Figure FDA0002929682080000028
n1和n2为给定的样本容量；
S5.4、在给定显著性水平α、样本容量n1与n2时，根据样本值计算
Figure FDA0002929682080000029
sW的值，对检验统计量的观测值作以下判断：若t≥tα(n1+n2-2)，则拒绝H3，即星基增强系统伪距域的完好性数据不可信；若t＜tα(n1+n2-2)，则接受H3，即星基增强系统伪距域的完好性数据可信；
当实测UDRE数据与卫星广播UDRE′数据的方差检验结果非齐性时，使用校正T检验对实测UDRE数据与卫星广播UDRE′数据的均值差进行检验的方法为：对实测UDRE数据与卫星广播UDRE′数据进行校正，对校正后的实测UDRE与卫星广播UDRE′数据进行方差齐性检验，若校正后的实测UDRE与卫星广播UDRE′数据方差检验结果齐性，则进行T检验；若校正后实测UDRE与卫星广播UDRE′数据的方差检验结果非齐性，则调整给定显著性水平，再进行T检验。
2.如权利要求1所述的一种星基增强系统完好性数据的统计检验方法，其特征在于，S2包括如下步骤：
S2.1、通过卫星导航电文信息对伪距观测值进行伪距误差修正，得到修正后的伪距观测量；
S2.2、对修正后的伪距观测量进行SBAS差分修正，得到修正后的伪距值；
S2.3、用修正后的伪距值与真距值做差计算得到实测UDRE的值。
3.如权利要求2所述的一种星基增强系统完好性数据的统计检验方法，其特征在于，S2.1中采用的伪距误差修正公式为：P＝ρ-Δt(s)+Δtu-dr+dion+dtrop+εp，其中，P为修正后的伪距观测值，ρ为卫星与用户接收机的几何距离，Δt(s)为广播星历引起的伪距修正值，Δtu为监测站接收机钟差引起的伪距修正值，dr为相对论效应伪距修正值，dion为电离层误差修正值，dtrop为对流层误差修正值，εp为GPS误差修正后的残差。
4.如权利要求1所述的一种星基增强系统完好性数据的统计检验方法，其特征在于，采用处理数据缺失值或剔除数据的异常值对实测UDRE与卫星广播UDRE′数据进行校正。
CN202010723220.2A 2020-07-24 2020-07-24 一种星基增强系统完好性数据的统计检验方法 Active CN111679294B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010723220.2A CN111679294B (zh) 2020-07-24 2020-07-24 一种星基增强系统完好性数据的统计检验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010723220.2A CN111679294B (zh) 2020-07-24 2020-07-24 一种星基增强系统完好性数据的统计检验方法

Publications (2)

Publication Number Publication Date
CN111679294A CN111679294A (zh) 2020-09-18
CN111679294B true CN111679294B (zh) 2021-03-19

Family

ID=72438280

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010723220.2A Active CN111679294B (zh) 2020-07-24 2020-07-24 一种星基增强系统完好性数据的统计检验方法

Country Status (1)

Country Link
CN (1) CN111679294B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113406673B (zh) * 2021-06-17 2023-03-24 哈尔滨工程大学 一种广域差分增强系统实时修正服务完好性监测方法
CN113238257A (zh) * 2021-07-12 2021-08-10 航天科工通信技术研究院有限责任公司 基于单接收机载波相位差分的gnss欺骗干扰检测方法
CN113534205B (zh) * 2021-09-16 2021-12-17 长沙海格北斗信息技术有限公司 卫星导航信号的异常判定方法、卫星导航方法及接收机
CN115372919B (zh) * 2022-08-30 2023-06-09 中国船舶集团有限公司第七二三研究所 一种基于t检验的雷达目标回波模拟性能评估方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101299063A (zh) * 2008-06-30 2008-11-05 北京航空航天大学 多星座sbas系统时差改正方法
CN101866009A (zh) * 2010-06-10 2010-10-20 北京航空航天大学 一种适用于卫导系统缺星情况下的完好性监测方法
CN103870679A (zh) * 2014-02-24 2014-06-18 苏州热工研究院有限公司 一种夏比冲击试样重组技术中最小插入段尺寸的判断方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107766424B (zh) * 2017-09-13 2020-09-15 深圳市宇数科技有限公司 一种数据探索管理方法、系统、电子设备及存储介质

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101299063A (zh) * 2008-06-30 2008-11-05 北京航空航天大学 多星座sbas系统时差改正方法
CN101866009A (zh) * 2010-06-10 2010-10-20 北京航空航天大学 一种适用于卫导系统缺星情况下的完好性监测方法
CN103870679A (zh) * 2014-02-24 2014-06-18 苏州热工研究院有限公司 一种夏比冲击试样重组技术中最小插入段尺寸的判断方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Estimated t test and F test with Interval-Censored Normal Data";Taichi KARIYA;《KAWASAKI MED.J.》;19841231;第10卷(第3期);第197-206页 *
"基于t检验的星基增强系统完好性数据核验";王娇;《电讯技术》;20171031;见正文第2-3节及图1-2 *
"方差非齐及小样本下总体均值差检验的探讨";林正大 等;《上海师范大学学报（自然科学版）》;19951231;正文第2-3节 *

Also Published As

Publication number Publication date
CN111679294A (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
CN111679294B (zh) 一种星基增强系统完好性数据的统计检验方法
CN109521443B (zh) 一种探测星历异常的方法
EP2544024B1 (en) Satellite navigation system fault detection based on threshold derived from test statatistics
JP6119838B2 (ja) 状態検出方法、補正値処理装置、測位システム、および状態検出プログラム
CN109359270B (zh) 北斗地基增强系统完好性风险监测的阈值模型建立方法
US11409002B2 (en) Method for operating a plurality of GNSS receivers for detecting satellite signal deformation
CN101419275B (zh) 基于多接收机的局域机场监视方法和系统
JP2011517770A (ja) 衛星ナビゲーションシステムのインテグリティのリアルタイム監視用の装置および方法
CN111965668B (zh) 一种面向卫星多故障的raim方法
CN107807368B (zh) 基于卫星观测量的定位方法和定位系统
CN109031356B (zh) 特征斜率加权的最小二乘残差接收机自主完好性监测方法
CN105301617A (zh) 一种卫星导航系统中的整周模糊度有效性检验方法
CN111856513A (zh) 卫星观测值获取方法、装置、计算机设备和存储介质
CN111435167B (zh) 一种基于北斗三号的接收机数据质量分析方法和装置
CN114859387A (zh) 静态站周跳探测的方法、装置、设备及计算机存储介质
CN115201872A (zh) 一种基于选星的gnss伪距粗差探测方法及系统
CN114417552A (zh) 一种模糊度确认方法、存储介质以及电子设备
CN113835105B (zh) 一种基于gnss模拟器的gbas完好性监测方法
CN116719073B (zh) 面向gnss解算域的粗差探测与剔除方法
CN105204047B (zh) 一种卫星导航系统中观测量单个粗差的探测与修复方法
CN105652298A (zh) 一种bds三频伪距相位组合的周跳探测修复方法及装置
CN112083463B (zh) 模糊度是否固定正确的检测方法及装置、定位终端
Zhang et al. Improved method for single and multiple GNSS faults exclusion based on consensus voting
Wang et al. Real-time assessment of GNSS observation noise with single receivers
WO2023123147A1 (zh) 一种数据预处理方法、数据预处理装置和芯片

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant