CN114152962A - Method for determining service range of satellite-based augmentation system - Google Patents

Method for determining service range of satellite-based augmentation system Download PDF

Info

Publication number
CN114152962A
CN114152962A CN202111344000.XA CN202111344000A CN114152962A CN 114152962 A CN114152962 A CN 114152962A CN 202111344000 A CN202111344000 A CN 202111344000A CN 114152962 A CN114152962 A CN 114152962A
Authority
CN
China
Prior art keywords
error
satellite
clock
orbit
clock error
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
CN202111344000.XA
Other languages
Chinese (zh)
Other versions
CN114152962B (en
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 CN202111344000.XA priority Critical patent/CN114152962B/en
Priority claimed from CN202111344000.XA external-priority patent/CN114152962B/en
Publication of CN114152962A publication Critical patent/CN114152962A/en
Application granted granted Critical
Publication of CN114152962B publication Critical patent/CN114152962B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40Correcting position, velocity or attitude

Abstract

The invention provides a method for determining a service range of a satellite-based augmentation system, which comprises the steps of selecting a target research area, resolving current epoch global navigation satellite system satellite augmentation orbit/clock error and user differential ranging error information, calculating orbit/clock error comprehensive errors projected by each satellite on grid points in the current epoch, counting target area grid point UDRE enveloping probability, calculating grid point ionized layer delay error values and grid ionized layer vertical error information, counting target area grid point GIVE enveloping probability, carrying out merging realization on public areas under integrity parameter enveloping requirements, determining a final SBAS service area, and determining a merged overlapping area as the final service area. The invention is based on the requirement to avoid the occurrence of integrity risk events: the integrity parameter corrects the residual error to form an envelope for the corresponding correction number with the envelope probability of 99.9 percent, the processing method is reasonable, and the processing efficiency is high.

Description

Method for determining service range of satellite-based augmentation system
Technical Field
The invention relates to the field of Satellite navigation enhancement, in particular to a method for determining the service range of a Satellite Based Augmentation System (SBAS).
Background
The satellite-based augmentation system is an augmentation system created for improving the positioning accuracy and integrity performance of a basic navigation system, and ephemeris correction numbers and integrity information are mainly broadcast to a user through Geostationary orbit satellites (GEO), so that integrity service is provided for a user side. Currently, SBAS in formal operation in the world include Wide Area Augmentation System (WAAS) in the united states, Geostationary Navigation Augmentation Service System (EGNOS) in europe, multi-function Satellite-based Augmentation System (MSAS) in japan, Global Positioning System (GPS) assisted static orbit Augmentation Navigation System (gan) in japan, and Global Positioning System (GPS) in india, and SBAS systems developed in other countries, such as difference Monitoring System (sdc) for Differential tracking and correction System (sdc) in russia, and enhanced Satellite System (KASS, koion Satellite Augmentation System) in korea, etc.
And the WAAS calculates corresponding protection levels according to different user service level requirements through the integrity parameters, and further gives a specific service range according to the protection levels. EGNOS utilizes the existing reference station, and determines a service range based on the constraint that a Horizontal Navigation System Error (HNSE) is better than 3m and a Vertical Navigation System Error (VNSE) is better than 4m under the effective Navigation positioning Precision Approach (PA) result. Other enhancement systems only give a rough service scope and have no explicit data or reports to support.
Therefore, specific methods for determining the service range are not disclosed abroad, and related documents are not introduced at home.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a method for determining the service range of a satellite-based augmentation system. The invention provides a method for backstepping an SBAS service range based on an SBAS integrity parameter 99.9% envelope probability requirement based on BDSBAS development current situation and innovativeness
The technical scheme adopted for solving the technical problems comprises the following specific steps:
step 1: selecting a target research area as a primary service area of the satellite-based augmentation system;
step 2: resolving information of a current epoch Global Navigation Satellite System (GNSS) Satellite enhanced orbit/clock Error and a User Differential Range Error (UDRE) with a calculation period of 2 minutes, and updating SBAS orbit clock Error correction information once in 2 minutes;
and step 3: dividing a target research area by latitude and longitude differences of 5 degrees multiplied by 5 degrees, calculating the orbit/clock error comprehensive error projected by each satellite on grid points of the current epoch, and judging whether the orbit/clock error comprehensive error of the UDRE is enveloped; the solved UDRE value is larger than the track/clock difference comprehensive error absolute value, namely envelope;
and 4, step 4: counting the UDRE enveloping probability of the target area grid points, and determining the area S with the UDRE enveloping probability more than 99.9 percent(UDRE)
And 5: calculating Grid point Ionospheric delay Error values and Grid Ionospheric Vertical Error (GIVE) information, wherein the calculation period is 5 minutes, the SBAS Ionospheric information is broadcasted once in 5 minutes, and whether the GIVE envelops the Ionospheric delay Error is judged; the solved GIVE value is larger than the absolute value of the ionosphere delay error, namely envelope;
step 6: counting the GIVE enveloping probability of grid points of the target region, and determining the region S with the GIVE enveloping probability more than 99.9 percent(GIVE)
And 7: merging public areas under the requirement of 99.9% envelope of integrity parameters (UDRE, GIVE), and determining the final SBAS service area S(SBAS)And the merged overlapping area is the final service area.
In the step 2, the current epoch GNSS satellite enhanced orbit/clock error and the UDRE information are solved;
the method comprises the following steps of calculating the satellite position and clock error through a received navigation message, and correcting the satellite position and clock error by using the orbit/clock error correction number in the satellite-based enhanced message to obtain an enhanced orbit/clock error result, wherein the method comprises the following specific steps:
2.1, calculating the position and clock error [ X ] of the satellite in the geocentric geostationary coordinate system by using the navigation messageE YE ZE BE];
2.2, calculating the orbit/clock error correction number by utilizing the satellite-based enhanced message:
the slow varying correction information comprises a satellite ephemeris slow varying correction and a satellite clock slow varying correction, and in the single-frequency enhanced message, the slow varying correction information is broadcast by a message 25;
the slowly varying correction of satellite ephemeris is calculated by:
Figure BDA0003353307370000021
[δxk δyk δzk]Tfor the track correction at the current time, [ δ x δ y δ z [ ]]TIn order to be a track deviation, the track deviation,
Figure BDA0003353307370000022
is the rate of change of track deviation;
the satellite clock slow-varying correction is calculated by:
δΔtSV(t)=δaf0+δaf1(t-t0)+δafG0where t is the current time, δ Δ tSV(t) is the current time clock correction, δ af0To clock skew, δ af1Is the rate of change of clock skew, t0To correct the reference time, deltaafG0Parameters were corrected for GLONASS satellites, broadcast in the text 12, which is 0 for non-GLONASS satellites;
2.3, correcting the satellite position and the clock error by using the correction number:
[XC YC ZC BC]T=[XE YE ZE BE]T+[δxk δyk δzk δΔtSV*c]T
wherein, [ X ]C YC ZC BC]TThe corrected enhanced orbit/clock error at the time t;
2.4 UDRE information is obtained via text 6.
In the step 3, the target area is segmented by latitude and longitude differences of 5 degrees multiplied by 5 degrees, the orbit/clock error comprehensive error projected by each satellite of the current epoch on the grid point is calculated, and whether the orbit/clock error comprehensive error is enveloped by the UDRE is judged;
respectively subtracting the corrected enhanced orbit/clock error from a post-operation precise orbit/clock error product to obtain a corrected error, projecting the corrected error to grid points, and judging whether the comprehensive error of the UDRE to the orbit/clock error is enveloped or not; the method comprises the following specific steps:
3.1, taking a difference value between the corrected satellite position and clock error and the precise satellite orbit and clock error;
[ΔX ΔY ΔZ ΔB]T=[XC YC ZC BC]T-[XP YP ZP BP]T
wherein, [ X ]C YC ZC BC]TFor the corrected satellite position and clock offset at time t, [ X [ ]P YP ZP BP]TFor the precise satellite orbit and clock error at time t, [ Δ X Δ Y Δ Z Δ B [ ]]TThe difference value between the satellite position and clock error corrected at the moment t and the precise satellite orbit and clock error;
3.2, calculating the projection of the corrected residual error on the grid point of the target area;
Figure BDA0003353307370000031
wherein, Delta R is the projection of the orbit/clock error synthetic error on grid points, luserIs a unit direction vector from the satellite to the grid point;
and 3.3, judging whether the comprehensive error of the UDRE to the track/clock error is enveloped.
The step 5: calculating the ionized layer delay error value of the grid point and GIVE information (the calculation period is 5 minutes, and the SBAS ionized layer information is broadcasted once in 5 minutes), and judging whether the ionized layer delay error is enveloped by the GIVE;
and taking the ionospheric delay amount calculated by the global ionospheric grid model as a reference, obtaining an ionospheric error by subtracting the ionospheric delay amount from the ionospheric delay at the corresponding grid point provided in the SBAS message, and counting the envelope probability of the GIVE to the difference value. The method comprises the following specific steps:
5.1, calculating enhanced message ionosphere delay IC and GIVE values of grid points of a target region through messages 18 and 26, calculating corresponding grid point global ionosphere grid model ionosphere delay I of the grid points, and if necessary, carrying out plane geometric interpolation.
And 5.2, calculating the difference value delta I between the ionospheric correction number and the ionospheric delay.
ΔI=IC-I
And 5.3, judging whether the ionized layer correction error enveloped by the GIVE is present.
The invention has the beneficial effects that:
based on the requirement to avoid the occurrence of integrity risk events: the integrity parameter corrects the residual error to form an envelope for the corresponding correction number by 99.9 percent of envelope probability, and innovatively provides a method for reversely deducing the service range of the satellite-based augmentation system. The treatment method is reasonable and has high treatment efficiency.
Based on the principle and the steps of the invention, research thinking and theoretical support can be provided for determining the service range of the BDSBAS, and further reference is provided for industrial users to use the BDSBAS service, and the application and popularization of the BDSBAS are promoted.
Drawings
FIG. 1 is a flow chart of a service range determination method of a satellite-based augmentation system of the present invention.
Detailed Description
The invention is further illustrated with reference to the following figures and examples.
The invention relates to a method for determining the service range of a satellite-based augmentation system, which comprises the following specific steps:
step 1: selecting a target research area as a primary service area of the satellite-based augmentation system;
and according to the distribution range of the SBAS monitoring stations, the longitude and latitude are respectively expanded by 10-20 degrees to the periphery to serve as target research areas.
Step 2: resolving current epoch GNSS satellite enhanced orbit/clock error and UDRE information;
the method comprises the following steps of calculating the satellite position and clock error through a received navigation message, and correcting the satellite position and clock error by using the orbit/clock error correction number in the satellite-based enhanced message to obtain an enhanced orbit/clock error result, wherein the method comprises the following specific steps:
2.1, calculating the position and clock error [ X ] of the satellite in the geocentric geostationary coordinate system by using the navigation messageE YE ZE BE];
2.2, calculating the orbit/clock error correction number by utilizing the satellite-based enhanced message:
the slow varying correction information comprises a satellite ephemeris slow varying correction and a satellite clock slow varying correction, and in the single-frequency enhanced message, the slow varying correction information is broadcast by a message 25;
the slowly varying correction of satellite ephemeris is calculated by:
Figure BDA0003353307370000041
[δxk δyk δzk]Tfor the track correction at the current time, [ δ x δ y δ z [ ]]TIn order to be a track deviation, the track deviation,
Figure BDA0003353307370000042
is the rate of change of track deviation;
the satellite clock slow-varying correction is calculated by:
δΔtSV(t)=δaf0+δaf1(t-t0)+δafG0
where t is the current time, δ ΔtSV(t) is the current time clock correction, δ af0To clock skew, δ af1Is the rate of change of clock skew, t0To correct the reference time, deltaafG0Parameters were corrected for GLONASS satellites, broadcast in the text 12, which is 0 for non-GLONASS satellites;
2.3, correcting the satellite position and the clock error by using the correction number:
[XC YC ZC BC]T=[XE YE ZE BE]T+[δxk δyk δzk δΔtSV*c]T
wherein, [ X ]C YC ZC BC]TThe corrected enhanced orbit/clock error at the time t;
2.4, the UDRE information is obtained through a telegraph text 6;
and step 3: dividing the target area by latitude and longitude differences of 5 degrees multiplied by 5 degrees, calculating the orbit/clock error comprehensive error projected by each satellite of the current epoch on the grid point, and judging whether the orbit/clock error comprehensive error of the UDRE is enveloped;
respectively subtracting the corrected enhanced orbit/clock error from a post-operation precise orbit/clock error product to obtain a corrected error, projecting the corrected error to grid points, and judging whether the comprehensive error of the UDRE to the orbit/clock error is enveloped or not; the method comprises the following specific steps:
3.1, taking a difference value between the corrected satellite position and clock error and the precise satellite orbit and clock error;
[ΔX ΔY ΔZ ΔB]T=[XC YC ZC BC]T-[XP YP ZP BP]T
wherein, [ X ]C YC ZC BC]TFor the corrected satellite position and clock offset at time t, [ X [ ]P YP ZP BP]TFor the precise satellite orbit and clock error at time t, [ Δ X Δ Y Δ Z Δ B [ ]]TCorrected satellite position and clock offset and precise satellite orbit and clock offset for time tA difference of (d);
3.2, calculating the projection of the corrected residual error on the grid point of the target area;
Figure BDA0003353307370000051
wherein, Delta R is the projection of the orbit/clock error synthetic error on grid points, IuserIs a unit direction vector from the satellite to the grid point;
3.3, judging whether the comprehensive error of the UDRE to the track/clock error is enveloped;
and 4, step 4: counting the UDRE enveloping probability of the target area grid points, and determining the area S with the enveloping probability more than 99.9 percent(UDRE)
S(UDRE)=S(P(σudre≥ΔR)≥99.9%)
And 5: calculating the ionized layer delay error value of the grid point and GIVE information (the calculation period is 5 minutes, and the SBAS ionized layer information is broadcasted once in 5 minutes), and judging whether the ionized layer delay error is enveloped by the GIVE;
and taking the ionospheric delay amount calculated by the global ionospheric grid model as a reference, obtaining an ionospheric error by subtracting the ionospheric delay amount from the ionospheric delay at the corresponding grid point provided in the SBAS message, and counting the envelope probability of the GIVE to the difference value. The method comprises the following specific steps:
5.1, calculating enhanced message ionosphere delay IC and GIVE values of grid points of a target region through messages 18 and 26, calculating corresponding grid point global ionosphere grid model ionosphere delay I of the grid points, and if necessary, carrying out plane geometric interpolation.
And 5.2, calculating the difference value delta I between the ionospheric correction number and the ionospheric delay.
ΔI=IC-I
5.3, judging whether the ionized layer correction error of the GIVE is enveloped;
step 6: counting the GIVE enveloping probability of the grid points of the target region, and determining the region S with the enveloping probability more than 99.9 percent(GIVE)
S(GIVE)=(P(σGIVE≥ΔI)≥99.9%)
And 7: the public areas under the envelope requirement of 99.9% of the integrity parameters (UDRE, GIVE) are superimposed:
S(SBAS service scope)=S(GIVE)∩S(UDRE)
Determining a final SBAS service area S(SBAS service scope)

Claims (4)

1. A method for determining service range of a satellite-based augmentation system is characterized by comprising the following steps:
step 1: selecting a target research area as a primary service area of the satellite-based augmentation system;
step 2: resolving satellite enhanced orbit/clock error and user differential ranging error information of a current epoch global navigation satellite system, wherein the calculation period is 2 minutes, and the SBAS orbit clock error correction information is updated once in 2 minutes;
and step 3: dividing a target research area by latitude and longitude differences of 5 degrees multiplied by 5 degrees, calculating the orbit/clock error comprehensive error projected by each satellite on grid points of the current epoch, and judging whether the orbit/clock error comprehensive error of the UDRE is enveloped; the solved UDRE value is larger than the track/clock difference comprehensive error absolute value, namely envelope;
and 4, step 4: counting the UDRE enveloping probability of the target area grid points, and determining the area S with the UDRE enveloping probability more than 99.9 percent(UDRE)
And 5: calculating the ionized layer delay error value of the grid point and the vertical error information of the grid ionized layer, wherein the calculation period is 5 minutes, the SBAS ionized layer information is broadcasted once in 5 minutes, and whether the ionized layer delay error is enveloped by the GIVE is judged; the solved GIVE value is larger than the absolute value of the ionosphere delay error, namely envelope;
step 6: counting the GIVE enveloping probability of grid points of the target region, and determining the region S with the GIVE enveloping probability more than 99.9 percent(GIVE)
And 7: merging public areas under the requirement of 99.9% envelope of integrity parameters (UDRE, GIVE), and determining the final SBAS service area S(SBAS)And the merged overlapping area is the final service area.
2. The method of claim 1, wherein the service range of the satellite-based augmentation system is determined by:
in the step 2, the current epoch GNSS satellite enhanced orbit/clock error and the UDRE information are solved;
the method comprises the following steps of calculating the satellite position and clock error through a received navigation message, and correcting the satellite position and clock error by using the orbit/clock error correction number in the satellite-based enhanced message to obtain an enhanced orbit/clock error result, wherein the method comprises the following specific steps:
2.1, calculating the position and clock error [ X ] of the satellite in the geocentric geostationary coordinate system by using the navigation messageE YE ZE BE];
2.2, calculating the orbit/clock error correction number by utilizing the satellite-based enhanced message:
the slow varying correction information comprises a satellite ephemeris slow varying correction and a satellite clock slow varying correction, and in the single-frequency enhanced message, the slow varying correction information is broadcast by a message 25;
the slowly varying correction of satellite ephemeris is calculated by:
Figure FDA0003353307360000011
[δxk δyk δzk]Tfor the track correction at the current time, [ δ x δ y δ z [ ]]TIn order to be a track deviation, the track deviation,
Figure FDA0003353307360000012
is the rate of change of track deviation;
the satellite clock slow-varying correction is calculated by:
δΔtSV(t)=δαf0+δαf1(t-t0)+δαfG0
where t is the current time, δ Δ tSV(t) is the current time clock correction, δ αf0To clock skew, δ αf1Is the rate of change of clock skew, t0When referred to as correction numbersEtching of delta alphafG0Parameters were corrected for GLONASS satellites, broadcast in the text 12, which is 0 for non-GLONASS satellites;
2.3, correcting the satellite position and the clock error by using the correction number:
[XC YC ZC BC]T=[XE YE ZE BE]T+[δxk δyk δzk δΔtSV*c]T
wherein, [ X ]C YC ZC BC]TThe corrected enhanced orbit/clock error at the time t;
2.4 UDRE information is obtained via text 6.
3. The method of claim 1, wherein the service range of the satellite-based augmentation system is determined by:
in the step 3, the target area is segmented by latitude and longitude differences of 5 degrees multiplied by 5 degrees, the orbit/clock error comprehensive error projected by each satellite of the current epoch on the grid point is calculated, and whether the orbit/clock error comprehensive error is enveloped by the UDRE is judged;
respectively subtracting the corrected enhanced orbit/clock error from a post-operation precise orbit/clock error product to obtain a corrected error, projecting the corrected error to grid points, and judging whether the comprehensive error of the UDRE to the orbit/clock error is enveloped or not; the method comprises the following specific steps:
3.1, taking a difference value between the corrected satellite position and clock error and the precise satellite orbit and clock error;
[ΔX ΔY ΔZ ΔB]T=[XC YC ZC BC]T-[XP YP ZP BP]T
wherein, [ X ]C YC ZC BC]TFor the corrected satellite position and clock offset at time t, [ X [ ]P YP ZP BP]TFor the precise satellite orbit and clock error at time t, [ Δ X Δ Y Δ Z Δ B [ ]]TCorrected satellite position and clock error and precise satellite orbit for time tDifference between track and clock;
3.2, calculating the projection of the corrected residual error on the grid point of the target area;
Figure FDA0003353307360000021
wherein, Delta R is the projection of the orbit/clock error synthetic error on grid points, luserIs a unit direction vector from the satellite to the grid point;
and 3.3, judging whether the comprehensive error of the UDRE to the track/clock error is enveloped.
4. The method of claim 1, wherein the service range of the satellite-based augmentation system is determined by:
the step 5: calculating the ionized layer delay error value of the grid point and GIVE information (the calculation period is 5 minutes, and the SBAS ionized layer information is broadcasted once in 5 minutes), and judging whether the ionized layer delay error is enveloped by the GIVE;
taking the ionospheric delay amount calculated by the global ionospheric grid model as a reference, obtaining an ionospheric error by making a difference with the ionospheric delay at the corresponding grid point provided in the SBAS message, and counting the envelope probability of the GIVE to the difference; the method comprises the following specific steps:
5.1, calculating enhanced message ionosphere delay IC and GIVE values of grid points of a target region through messages 18 and 26, calculating corresponding grid point global ionosphere grid model ionosphere delay I of the grid points, and if necessary, carrying out plane geometric interpolation;
5.2, calculating a difference value delta I between the ionospheric correction number and the ionospheric delay;
ΔI=IC-I
and 5.3, judging whether the ionized layer correction error enveloped by the GIVE is present.
CN202111344000.XA 2021-11-14 Method for determining service range of star-based enhanced system Active CN114152962B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111344000.XA CN114152962B (en) 2021-11-14 Method for determining service range of star-based enhanced system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111344000.XA CN114152962B (en) 2021-11-14 Method for determining service range of star-based enhanced system

Publications (2)

Publication Number Publication Date
CN114152962A true CN114152962A (en) 2022-03-08
CN114152962B CN114152962B (en) 2024-05-14

Family

ID=

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115575978A (en) * 2022-11-23 2023-01-06 中国矿业大学(北京) Grid ionosphere delay correction method and device for user side and receiver

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104732085A (en) * 2015-03-23 2015-06-24 北京航空航天大学 Satellite navigation satellite-based augmentation system availability prediction method
CN106154291A (en) * 2016-07-08 2016-11-23 北京航空航天大学 SBAS availability Forecasting Methodology based on satellite geometry structure and device
US20190271782A1 (en) * 2017-07-18 2019-09-05 Wuhan University Ionospheric delay correction method for leo satellite augmented navigation systems
CN112099063A (en) * 2020-08-25 2020-12-18 中国电子科技集团公司第二十研究所 Method for quickly searching Beidou satellite-based enhanced user error maximum projection direction
CN113253303A (en) * 2021-05-13 2021-08-13 中国电子科技集团公司第二十研究所 Method for monitoring performance of single-frequency satellite-based augmentation system in real time

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104732085A (en) * 2015-03-23 2015-06-24 北京航空航天大学 Satellite navigation satellite-based augmentation system availability prediction method
CN106154291A (en) * 2016-07-08 2016-11-23 北京航空航天大学 SBAS availability Forecasting Methodology based on satellite geometry structure and device
US20190271782A1 (en) * 2017-07-18 2019-09-05 Wuhan University Ionospheric delay correction method for leo satellite augmented navigation systems
CN112099063A (en) * 2020-08-25 2020-12-18 中国电子科技集团公司第二十研究所 Method for quickly searching Beidou satellite-based enhanced user error maximum projection direction
CN113253303A (en) * 2021-05-13 2021-08-13 中国电子科技集团公司第二十研究所 Method for monitoring performance of single-frequency satellite-based augmentation system in real time

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《中国公路学报》编辑部;: "中国汽车工程学术研究综述・2017", 中国公路学报, no. 06, 15 June 2017 (2017-06-15) *
崔瑞云;倪育德;王琳琳;刘逸;: "北斗卫星导航系统完好性参数研究", 现代导航, no. 01, 15 February 2015 (2015-02-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115575978A (en) * 2022-11-23 2023-01-06 中国矿业大学(北京) Grid ionosphere delay correction method and device for user side and receiver

Similar Documents

Publication Publication Date Title
EP3035080B1 (en) Navigation satellite system positioning involving the generation of correction information
CN105842720B (en) A kind of wide area precision real-time location method based on carrier phase
EP2260319B1 (en) Improved sbas receiver
US11226416B1 (en) System and method to reduce PPP filter convergence time using LEO frequency band signals
CN111983641B (en) Method for generating Beidou satellite-based augmentation system integrity parameters in real time
CN113253303B (en) Method for monitoring performance of single-frequency star-based enhanced system in real time
CN112099056B (en) Method for predicting Beidou satellite-based enhanced service performance in real time
CN113406673B (en) Method for monitoring real-time correction service integrity of wide area difference enhancement system
CN112099063B (en) Method for quickly searching Beidou satellite-based enhanced user error maximum projection direction
CN103076615A (en) Ionized layer delaying estimation method of Beidou navigation satellite system and corresponding receiving device
CN103698790A (en) Beidou and GPS (global positioning system) double-system wide-lane carrier phase mixed-frequency inter-satellite differential combination method
CN115407371A (en) PPP-B2B-based real-time high-precision time transfer method and device
Tsai et al. Performance assessment on expanding SBAS service areas of GAGAN and MSAS to Singapore region
Euler et al. Improvement of positioning performance using standardized network RTK messages
CN103675858A (en) B1/L1 carrier phase frequency mixing difference method for Beidou system and GPS (Global Positioning System)
CN108254762B (en) Pseudo-range differential positioning method and system
CN114152962B (en) Method for determining service range of star-based enhanced system
CN114152962A (en) Method for determining service range of satellite-based augmentation system
CN115902968A (en) PPP terminal positioning method based on Beidou third GEO broadcast enhancement information
Luo et al. Benefit of sparse reference network in BDS single point positioning with single-frequency measurements
CN115685283A (en) PPP-B2B-based real-time vehicle-mounted positioning method and system
Jin et al. Performance analysis of SBAS ephemeris corrections and integrity algorithms in China region
Paziewski et al. Application of SBAS pseudorange and carries phase signals to precise instantaneous single-frequency positioning
CN109738912B (en) Method for realizing fixed point time service based on GNSS satellite signals
Collins et al. Real-Time WADGPS Corrections from Undifferenced Carrier Phase

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