CN110988928A - 一种基于长度约束的卫星导航着陆完好性监测方法 - Google Patents

一种基于长度约束的卫星导航着陆完好性监测方法 Download PDF

Info

Publication number
CN110988928A
CN110988928A CN201911290143.XA CN201911290143A CN110988928A CN 110988928 A CN110988928 A CN 110988928A CN 201911290143 A CN201911290143 A CN 201911290143A CN 110988928 A CN110988928 A CN 110988928A
Authority
CN
China
Prior art keywords
node
difference
satellite
pseudo
double
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.)
Granted
Application number
CN201911290143.XA
Other languages
English (en)
Other versions
CN110988928B (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.)
CETC 20 Research Institute
Original Assignee
CETC 20 Research Institute
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 CETC 20 Research Institute filed Critical CETC 20 Research Institute
Priority to CN201911290143.XA priority Critical patent/CN110988928B/zh
Publication of CN110988928A publication Critical patent/CN110988928A/zh
Application granted granted Critical
Publication of CN110988928B publication Critical patent/CN110988928B/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/13—Receivers
    • G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明提供了一种基于长度约束的卫星导航着陆完好性监测方法，先任意选择主节点与基准卫星，然后通过计算多节点间双差值和伪距双差残差，从而进行多节点综合判断。本发明采用双差的方式进行完好性监测，由于双差观测量消除了大部分公共误差且转化得到的距离理论双差值精度高，因此能够更好的对故障进行监测识别。在现有算法运行之前执行，并不影响现有算法体系，可与现有算法同时运行，提升***的完好性监测能力。仅涉及双差、坐标转换、比对等运算，算法简单，占用运算资源较小。

Description

一种基于长度约束的卫星导航着陆完好性监测方法
技术领域
本发明涉及卫星导航定位领域，是卫星导航着陆***中一种利用地面基准站内多个基准接收机天线之间固定的基线向量作为约束信息的完好性监测方法。
背景技术
完好性监测是卫星导航着陆***的一个关键环节，决定着***所能提供的服务等级，关系着飞机的着陆安全。在卫星导航着陆***中，地面基准站设备通常配置的4部基准接收机用于完好性监测处理，对可能的故障源(包括卫星钟、卫星星历、电离层、对流层、多径和接收机)进行分类的监视。完好性监测算法利用3～4台基准接收机信号测量提供的信息通过信号质量件监视(SQM)、量测质量监视(MQM)、数据质量监视(DQM)等一系列方法完成故障监测和告警，实现故障监测完好性。SQM对卫星导航信号功率异常和相关峰畸变情况进行监测；DQM连续检查卫星星历和星钟参数，验证卫星导航数据足够可信；MQM主要对于生成的各种观测质量的好坏和状态加以判断和鉴别。上述一系列算法仅通过对各台基准接收机数据进行独立运算分析以及多台基准接收机的数据冗余性比对实现故障的监测。
上述算法并未利用天线间准确基线关系进行完好性监测，而在卫星导航着陆***中，地面站基准接收机天线布设在固定位置，且每台基准接收机天线位置事先进行了标定，同时多天线间的准确基线关系可通过测量的方式获得。由于通过测定多天线间准确基线关系可达到毫米精度，同时天线间准确基线关系能够转化为“卫星-基准接收机”间的距离理论双差值，通过将基准接收机测量获得的伪距值进行双差运算获得的伪距双差值与转化获得的距离双差值进行比对，即可对伪距的异常进行监测，实现完好性监测的目的。
发明内容
为了克服现有技术的不足，本发明提供一种基于长度约束的卫星导航着陆完好性监测方法。本发明介绍一种基于长度约束的卫星导航着陆完好性监测方法，利用多个节点(节点数不少于3个)卫星导航基准接收机天线之间固定的基线向量作为约束信息进行完好性监测。基于长度约束的卫星导航着陆完好性监测方法总体框图如图1所示。其中多节点包括节点1、节点2....节点n(n＝1,2,3…N，N≥3)，共N个节点，其中节点M为在节点n中选定的主节点(M＝1,2,3…N，N≥3)。由于N个节点间任意两个节点间的基线矢量关系可以事先测定，因此在选定主节点M后，可以节点1-节点M、节点2-节点M....节点n-节点M(M,n＝1,2,3…N，N≥3，n≠M)间的基线关系作为约束信息进行完好性监测。
本发明解决其技术问题所采用的技术方案的具体步骤如下：
步骤1：主节点与基准卫星任意选择；
选择节点1、节点2、节点3、....节点n中的一个节点作为主节点，记作M，其中，n＝1,2,3…N，N≥3；选择节点1、节点2、节点3、....节点n的共视卫星集合Q(共L颗卫星)中的一颗卫星作为基准卫星，记作S；
步骤2：多节点间双差值计算；
在选定主节点M以及基准卫星S后，分别计算除去M节点后的其他节点与节点M间的伪距双差观测值与伪距双差理论值；
伪距双差观测值
Figure BDA0002318745100000021
通过对不同节点间、不同卫星间基准接收机测量获得的伪距观测量的进行两次作差获得：
Figure BDA0002318745100000022
其中n代表节点编号，n＝1,2,3...N,n≠M；k代表N个节点中共视卫星编号，k≠S；
Figure BDA0002318745100000023
为构建的伪距双差观测量；
Figure BDA0002318745100000024
为节点n与节点M间卫星k与卫星S真实距离的双差观测值；
Figure BDA0002318745100000025
为伪距双差观测量测量噪声；
伪距双差理论值
Figure BDA0002318745100000026
通过将事先测量获得的基线向量BnM进行转换，转换为不同节点间、不同卫星间基准接收机所对应的双差值获得：
Figure BDA0002318745100000027
其中，BnM为事先测量的节点n与节点M间的基线向量，ek、eS分别为节点n与节点M到卫星k与卫星S的方向矢量，n代表节点编号，n＝1,2,3...N,n≠M；k代表N个节点中共视卫星编号，k≠S；
Figure BDA0002318745100000028
为节点n与节点M间卫星k与卫星S的伪距双差理论值；
步骤3：伪距双差残差计算；
根据步骤2计算得到的伪距双差理论值、伪距双差观测值，并将其作差得到双差残差值
Figure BDA0002318745100000031
Figure BDA0002318745100000032
N≥3，n≠M，
Figure BDA0002318745100000033
代表由伪距双差理论值与伪距双差观测值构建的伪距双差残差值；
步骤4：多节点综合判断；
根据步骤3计算得到的伪距双差残差值
Figure BDA0002318745100000034
构建伪距双差残差平均值
Figure BDA0002318745100000035
Figure BDA0002318745100000036
分别计算伪距双差残差值与伪距双差平均值的距离
Figure BDA0002318745100000037
Figure BDA0002318745100000038
Figure BDA0002318745100000039
记
Figure BDA00023187451000000310
中最大值为
Figure BDA00023187451000000311
将
Figure BDA00023187451000000312
与预设门限
Figure BDA00023187451000000313
进行比较，若距离
Figure BDA00023187451000000314
小于门限，则判定通过；若
Figure BDA00023187451000000315
大于等于门限，则判定节点n的卫星k发生故障；按照99.7％置信概率对预设门限进行
Figure BDA00023187451000000316
设置；
Figure BDA00023187451000000317
Figure BDA00023187451000000318
步骤1中为保证选择的基准卫星具有较好的数据质量，进行以下限定：
1)该颗卫星载噪比在所有共视卫星中最大；
2)该颗卫星仰角大于45度。
步骤4中为了更加真实的评估伪距双差观测量测量噪声
Figure BDA00023187451000000319
采用仰角加权模型对
Figure BDA00023187451000000320
进行取值。由于双差观测量消除了大部分公共误差，仅剩接收机热噪声以及多路径误差影响，因此参照卫星导航着陆***性能标准，选取模型为
Figure BDA00023187451000000321
其中参数a、b、c为可变参数，θ为卫星k的仰角，按照标准规定的中精度等级参数取值，设置a、b、c分别为0.16、1.07、15.5，因此
Figure BDA00023187451000000322
本发明的有益效果在于：
1)本发明采用双差的方式进行完好性监测，由于双差观测量消除了大部分公共误差且转化得到的距离理论双差值精度高，因此能够更好的对故障进行监测识别。
2)在现有算法运行之前执行，并不影响现有算法体系，可与现有算法同时运行，提升***的完好性监测能力。
3)本发明仅涉及双差、坐标转换、比对等运算，算法简单，占用运算资源较小。
附图说明
图1为本发明基于长度约束的卫星导航着陆完好性监测方法框图。
图2为本发明的基于长度约束的卫星导航着陆完好性监测方法流程图。
图3为本发明基线向量与双差理论值转换关系示意图。
图4为本发明天线位置与基线关系示意图。
图5为本发明完好性监测结果示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明是一种利用地面基准站内多个基准接收机天线之间固定的基线向量作为约束信息的卫星导航着陆***完好性监测的方法，具体步骤如图2所示：
步骤一：主节点与基准卫星选择
主节点与基准卫星选择任意。选择节点1、节点2、节点3、....节点n(n＝1,2,3…N，N≥3)中的一个节点作为主节点，记作M；选择节点1、节点2、节点3、....节点n(n＝1,2,3…N，N≥3)的共视卫星集合Q(共L颗卫星)中的一颗卫星作为基准卫星，记作S。要求选择的S卫星的载噪比在所有共视卫星中最大且该颗卫星仰角大于45度。
步骤二：多节点间双差值计算
在选定主节点M以及基准卫星S后，分别计算除去M节点的其他节点与节点M间的伪距双差观测值与伪距双差理论值，方法如下：
(1)伪距双差观测值计算
对不同节点间、不同卫星间基准接收机测量获得的伪距观测量的两次作差，得到伪距双差观测值如下：
Figure BDA0002318745100000041
式中n代表节点编号，n＝1,2,3...N,n≠M；k代表N个节点中共视卫星编号，k≠S；
Figure BDA0002318745100000042
为构建的伪距双差观测量；
Figure BDA0002318745100000043
为节点n与节点M间卫星k与卫星S真实距离的双差观测值；
Figure BDA0002318745100000044
为伪距双差观测量测量噪声。
(2)伪距双差理论值计算
将事先测量获得的基线关系进行转换，转换为不同节点间、不同卫星间基准接收机所对应的双差值，得到伪距双差理论值如下：
如图3所示，BnM为事先测量的节点n与节点M间的基线向量，ek、eS分别为节点n与节点M指向卫星k与卫星S的方向矢量。根据转换关系，可得：
Figure BDA0002318745100000051
式中n代表节点编号，n＝1,2,3...N,n≠M；k代表N个节点中共视卫星编号，k≠S；
Figure BDA0002318745100000052
为节点n与节点M间卫星k与卫星S理论距离的双差观测值。
步骤三：伪距双差残差计算
根据步骤二计算得到的伪距双差理论值、伪距双差观测值，并将其作差得到双差残差值
Figure BDA0002318745100000053
Figure BDA0002318745100000054
Figure BDA0002318745100000055
代表由伪距双差理论值与伪距双差观测值构建的伪距双差残差值，根据
Figure BDA0002318745100000056
发生跳变的情况即可判断卫星k是否发生异常。但由于
Figure BDA0002318745100000057
包含节点n、节点M的异常信息，因此
Figure BDA0002318745100000058
发生跳变时无法辨别是节点n与节点M的异常情况，因此需利用N个节点生成的N-1个
Figure BDA0002318745100000059
信息进行综合判断确定。
步骤四：多节点综合判断
根据步骤三计算得到的伪距双差残差值
Figure BDA00023187451000000510
构建伪距双差残差平均值
Figure BDA00023187451000000511
Figure BDA00023187451000000512
分别计算伪距双差残差值与伪距双差平均值的距离
Figure BDA00023187451000000513
Figure BDA00023187451000000514
记
Figure BDA00023187451000000515
中最大值为
Figure BDA00023187451000000516
将
Figure BDA00023187451000000517
与预设门限
Figure BDA00023187451000000518
进行比较，若距离
Figure BDA00023187451000000519
小于门限，则通过；若
Figure BDA00023187451000000520
大于等于门限，则判定节点n的卫星k发生故障。式中按照99.7％置信概率对预设门限进行
Figure BDA00023187451000000521
设置。
Figure BDA00023187451000000522
采用仰角加权模型对
Figure BDA00023187451000000523
进行取值，
Figure BDA00023187451000000524
θ为卫星k的仰角。
为验证该方法的实际性能，开展了如下试验。架设4副卫星导航接收天线，其中天线1与天线2基线距离100.08m，天线1与天线3基线距离145.798m，天线1与天线4基线距离99.634m。天线位置与基线位置关系如图4所示。
试验中，选择天线1作为主节点，选择10号卫星作为基准卫星，在46秒、91秒时刻在节点1的1号星伪距上分别设置1m、2m的故障，完好性监测结果如图5所示。右图看出在第46秒时，节点1的伪距双差残差值与伪距双差平均值的距离最大且超出门限值；在第91秒时，节点1、3、4的伪距双差残差值与伪距双差平均值的距离均超出门限，但节点1的伪距双差残差值与伪距双差平均值的距离最大。因此，在46秒、91秒分别判定节点1的1号卫星发生故障。因此基于长度约束的卫星导航着陆完好性监测方法能够将有效的对卫星导航着陆***的故障进行监测识别，降低***的完好性风险。

Claims (3)

1.一种基于长度约束的卫星导航着陆完好性监测方法，其特征在于包括下述步骤：
步骤1：主节点与基准卫星任意选择；
选择节点1、节点2、节点3、….节点n中的一个节点作为主节点，记作M，其中，n＝1,2,3…N，N≥3；选择节点1、节点2、节点3、….节点n的共视卫星集合Q(共L颗卫星)中的一颗卫星作为基准卫星，记作S；
步骤2：多节点间双差值计算；
在选定主节点M以及基准卫星S后，分别计算除去M节点后的其他节点与节点M间的伪距双差观测值与伪距双差理论值；
伪距双差观测值
Figure FDA0002318745090000011
通过对不同节点间、不同卫星间基准接收机测量获得的伪距观测量的进行两次作差获得：
Figure FDA0002318745090000012
其中n代表节点编号，n＝1,2,3...N,n≠M；k代表N个节点中共视卫星编号，k≠S；
Figure FDA0002318745090000013
为构建的伪距双差观测量；
Figure FDA0002318745090000014
为节点n与节点M间卫星k与卫星S真实距离的双差观测值；
Figure FDA0002318745090000015
为伪距双差观测量测量噪声；
伪距双差理论值
Figure FDA0002318745090000016
通过将事先测量获得的基线向量BnM进行转换，转换为不同节点间、不同卫星间基准接收机所对应的双差值获得：
Figure FDA0002318745090000017
其中，BnM为事先测量的节点n与节点M间的基线向量，ek、eS分别为节点n与节点M到卫星k与卫星S的方向矢量，n代表节点编号，n＝1,2,3...N,n≠M；k代表N个节点中共视卫星编号，k≠S；
Figure FDA0002318745090000018
为节点n与节点M间卫星k与卫星S的伪距双差理论值；
步骤3：伪距双差残差计算；
根据步骤2计算得到的伪距双差理论值、伪距双差观测值，并将其作差得到双差残差值
Figure FDA0002318745090000019
Figure FDA00023187450900000110
M,n＝1,2,3…N，N≥3，n≠M，
Figure FDA00023187450900000111
代表由伪距双差理论值与伪距双差观测值构建的伪距双差残差值；
步骤4：多节点综合判断；
根据步骤3计算得到的伪距双差残差值
Figure FDA00023187450900000112
构建伪距双差残差平均值
Figure FDA00023187450900000113
Figure FDA0002318745090000021
分别计算伪距双差残差值与伪距双差平均值的距离
Figure FDA0002318745090000022
Figure FDA0002318745090000023
(M,n＝1,2,3…N，N≥3，n≠M)；记
Figure FDA0002318745090000024
中最大值为
Figure FDA0002318745090000025
将
Figure FDA0002318745090000026
与预设门限
Figure FDA0002318745090000027
进行比较，若距离
Figure FDA0002318745090000028
小于门限，则判定通过；若
Figure FDA0002318745090000029
大于等于门限，则判定节点n的卫星k发生故障；按照99.7％置信概率对预设门限进行
Figure FDA00023187450900000210
设置；
Figure FDA00023187450900000211
Figure FDA00023187450900000212
M,n＝1,2,3…N，N≥3，n≠M。
2.根据权利要求1所述的一种基于长度约束的卫星导航着陆完好性监测方法，其特征在于：
步骤1中基准卫星进行以下限定：
1)基准卫星载噪比在所有共视卫星中最大；
2)基准卫星仰角大于45度。
3.根据权利要求1所述的一种基于长度约束的卫星导航着陆完好性监测方法，其特征在于：
步骤4中为了更加真实的评估伪距双差观测量测量噪声
Figure FDA00023187450900000213
采用仰角加权模型对
Figure FDA00023187450900000214
进行取值，选取模型为
Figure FDA00023187450900000215
其中参数a、b、c为可变参数，θ为卫星k的仰角，设置a、b、c分别为0.16、1.07、15.5，因此
Figure FDA00023187450900000216
CN201911290143.XA 2019-12-16 2019-12-16 一种基于长度约束的卫星导航着陆完好性监测方法 Active CN110988928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911290143.XA CN110988928B (zh) 2019-12-16 2019-12-16 一种基于长度约束的卫星导航着陆完好性监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911290143.XA CN110988928B (zh) 2019-12-16 2019-12-16 一种基于长度约束的卫星导航着陆完好性监测方法

Publications (2)

Publication Number Publication Date
CN110988928A true CN110988928A (zh) 2020-04-10
CN110988928B CN110988928B (zh) 2023-02-14

Family

ID=70094091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911290143.XA Active CN110988928B (zh) 2019-12-16 2019-12-16 一种基于长度约束的卫星导航着陆完好性监测方法

Country Status (1)

Country Link
CN (1) CN110988928B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112526563A (zh) * 2020-11-30 2021-03-19 北京航空航天大学 一种gnss信号质量监测方法及***
CN113253321A (zh) * 2021-04-14 2021-08-13 长沙学院 一种适用于接收机时间误差为秒级的agps定位方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090037105A1 (en) * 2007-08-01 2009-02-05 Astrium Gmbh Positioning Device and Positioning Method for a Satellite Navigation System
CN104570011A (zh) * 2014-12-26 2015-04-29 中国电子科技集团公司第二十研究所 一种卫星导航相对定位装置及其载波相位周跳修复方法
CN106597514A (zh) * 2016-12-15 2017-04-26 中国电子科技集团公司第二十研究所 一种基于多节点数据融合的卫星导航相对定位装置及方法
CN108732601A (zh) * 2018-04-16 2018-11-02 哈尔滨工业大学 基于差分卫星导航的垂直返回运载器着陆段导航方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090037105A1 (en) * 2007-08-01 2009-02-05 Astrium Gmbh Positioning Device and Positioning Method for a Satellite Navigation System
CN104570011A (zh) * 2014-12-26 2015-04-29 中国电子科技集团公司第二十研究所 一种卫星导航相对定位装置及其载波相位周跳修复方法
CN106597514A (zh) * 2016-12-15 2017-04-26 中国电子科技集团公司第二十研究所 一种基于多节点数据融合的卫星导航相对定位装置及方法
CN108732601A (zh) * 2018-04-16 2018-11-02 哈尔滨工业大学 基于差分卫星导航的垂直返回运载器着陆段导航方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
戴伟 等: "北斗卫星导航***伪距差分定位技术研究", 《全球定位***》, vol. 40, no. 2, 30 April 2015 (2015-04-30), pages 36 - 39 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112526563A (zh) * 2020-11-30 2021-03-19 北京航空航天大学 一种gnss信号质量监测方法及***
CN112526563B (zh) * 2020-11-30 2022-05-03 北京航空航天大学 一种gnss信号质量监测方法及***
CN113253321A (zh) * 2021-04-14 2021-08-13 长沙学院 一种适用于接收机时间误差为秒级的agps定位方法

Also Published As

Publication number Publication date
CN110988928B (zh) 2023-02-14

Similar Documents

Publication Publication Date Title
CN109359270B (zh) 北斗地基增强***完好性风险监测的阈值模型建立方法
US7994974B2 (en) Device and method for the real-time monitoring of the integrity of a satellite navigation system
CN103970997B (zh) 一种无人直升机传感器故障快速诊断方法
CN110068840B (zh) 一种基于伪距测量特征值提取的araim故障检测方法
AU2019204184B2 (en) A method for operating a plurality of GNSS receivers for detecting satellite signal deformation
US20090273511A1 (en) Method of operating a satellite navigation receiver
CN110988928B (zh) 一种基于长度约束的卫星导航着陆完好性监测方法
US10732289B1 (en) LDACS-based air-ground cooperative positioning and integrity monitoring method
CN102426350B (zh) 一种星载阵列天线测向通道幅相误差的确定方法
CN102135621B (zh) 一种多星座组合导航***的故障识别方法
CN103592656A (zh) 一种适用于星载导航接收机的自主完好性监测方法
WO2021052511A1 (zh) 北斗机载设备raim性能符合性测试方法及***
CN108682140B (zh) 一种基于压缩感知和自回归模型的增强型异常检测方法
CN101419275A (zh) 基于多接收机的局域机场监视方法和***
KR101477041B1 (ko) Ｄｇｎｓｓ 기준국 위성신호 이상 모니터링 시스템 및 이를 이용한 모니터링 방법
CN104331602A (zh) 基于马尔科夫链的gnss完好性检测率估计方法
CN109308518A (zh) 一种基于概率神经网络的监测***及其平滑参数优化方法
CN106443577A (zh) 一种针对星间射频相对测量的多径误差检测与消除方法
KR101503001B1 (ko) 지상설비 안테나와 기저선 길이 예측값을 이용한 위성항법시스템의 고장 판단 시스템 및 그 방법
CN111027721B (zh) 一种***故障定位方法
CN115453579A (zh) 基于北斗ppp-rtk的合成星历a类故障完好性监测方法及装置
CN112230247A (zh) 一种用于城市复杂环境下gnss完好性监测方法
CN113835105B (zh) 一种基于gnss模拟器的gbas完好性监测方法
US20110169693A1 (en) Integrity communication in a satellite navigation system
KR101040053B1 (ko) 의사거리 측정잡음 감시기반 위성전파항법시스템 위성시계 고장검출 및 고장위성 식별 방법

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