CN111766616A - Beidou second-order time transfer satellite-side multipath error correction method - Google Patents

Beidou second-order time transfer satellite-side multipath error correction method Download PDF

Info

Publication number
CN111766616A
CN111766616A CN202010539934.8A CN202010539934A CN111766616A CN 111766616 A CN111766616 A CN 111766616A CN 202010539934 A CN202010539934 A CN 202010539934A CN 111766616 A CN111766616 A CN 111766616A
Authority
CN
China
Prior art keywords
satellite
beidou
pseudo
range
observation
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.)
Pending
Application number
CN202010539934.8A
Other languages
Chinese (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.)
People's Liberation Army 61081 Unit
Original Assignee
People's Liberation Army 61081 Unit
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 People's Liberation Army 61081 Unit filed Critical People's Liberation Army 61081 Unit
Priority to CN202010539934.8A priority Critical patent/CN111766616A/en
Publication of CN111766616A publication Critical patent/CN111766616A/en
Pending legal-status Critical Current

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/396Determining accuracy or reliability of position or pseudorange measurements
    • 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/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07Cooperating 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
    • G01S19/073Cooperating 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 involving a network of fixed stations
    • 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/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/10Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals
    • G01S19/11Cooperating 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
    • G01S19/115Airborne or satellite based pseudolites or repeaters
    • 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
    • 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/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

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

The invention provides a method for correcting multipath errors of a Beidou second time transmission satellite terminal, which is used for solving the problem of multipath errors of the Beidou second time transmission satellite terminal. The Beidou second time transmission satellite-side multipath error correction method selects a MGEX observation station as a data source, after grouping, the MP sequence is used for segmenting and dividing each grouped data to generate a pseudo-range correction information table, and therefore during BDS CV and/or BDS PPP time transmission, pseudo-range correction node information is searched according to satellite types and signal frequency points and is fitted to correct a pseudo-range observation value. According to the invention, satellite-end multipath error correction is introduced in Beidou time transmission, the pseudo range observation value of the Beidou second satellite is corrected based on the observation data of the global MGEX observation station, and the pseudo range is directly corrected systematically at the user terminal, so that the influence of satellite-end multipath on time transmission is weakened, the problem of satellite-end multipath error transmission of the Beidou second time is solved, and the precision of the Beidou second time transmission is improved.

Description

Beidou second-order time transfer satellite-side multipath error correction method
Technical Field
The invention belongs to the field of navigation, and particularly relates to a Beidou second time transfer satellite-side multipath error correction method.
Background
The chinese BeiDou Navigation Satellite System (BDS) is a global Satellite Navigation System developed by the self in china, and formally provides Navigation, positioning and time service to asia-pacific region from 12 months in 2012, and is one of the four global Satellite Navigation systems. The BDS consists of a space section, a ground section and a user section, the Beidou No. two satellite which runs in orbit comprises 5 GEO (geostationary earth orbit) satellites (C01-C05), 7 IGSO (unsaturated geosynchronous satellite orbit) satellites (C06-C10, C13 and C16) and 3 MEO (medium earth orbit) satellites (C11, C12 and C14), and the Beidou No. three system is expected to be built in 6 months in 2020.
The time-keeping laboratories participating in the international time alignment mainly perform the remote time alignment through the united states global navigation satellite System (GPS) and the russian global navigation satellite System (GLONASS). The frequency stability of the GPS precision Point Positioning technology (PPP) in a day with a distance of hundreds or even thousands of kilometers can reach an order of E-15 to E-16, and the time transfer method thereof includes GPS Common-View time transfer (GPS Common-View, GPS cv), GPS All-View time transfer (GPS All-View, GPS AV), and GPS Carrier Phase time transfer (GPS Carrier-Phase, GPS CP). There is a significant reliability risk in relying on a single system for time transfer. In the long term, the international time comparison method gradually develops from a single system to a multi-system fusion time comparison direction. In 2015, a CGGTTS V2E co-view standard was formally released, and a co-view method of GPS, GLONASS, european union global navigation satellite system Galileo, BDS, and QZSS was clearly defined in the standard. Preliminary tests prove that the BDS CV time comparison standard deviation (STD) is about 2-3 ns, and the number of common vision satellites is about 2-3. Due to the limited service range of the Beidou II, the time transfer of the BDS PPP is limited to the Asia-Pacific region, and the STD of the BDS PPP time comparison is better than 0.5 ns.
The big dipper second satellite has satellite-end multipath error, and can generate system error of several decimeters to 1 meter for pseudo range. The pseudo-range deviation of the Beidou second satellite is closely related to the altitude angle and the signal frequency of the satellite, and the deviation is unrelated to the antenna type and the observation time period of the receiver. Therefore, the pseudorange bias is preliminarily concluded to be derived from Beidou satellite multipath. In the prior art, the research and application of the satellite-end multipath influence of the Beidou second satellite mainly focuses on satellite positioning, and the research of the influence of systematic errors on time transfer is lacked.
Disclosure of Invention
The embodiment of the invention provides a method for correcting the time transfer satellite-end multipath errors of a Beidou second satellite, which aims to solve the problem of the time transfer satellite-end multipath errors of the Beidou second satellite.
In order to achieve the above purpose, the embodiment of the invention adopts the following technical scheme:
a Beidou second time transfer satellite-end multipath error correction method comprises the following steps:
step S1, selecting data of a multi-system observation network MGEX observation station of the International Global navigation satellite service organization as a data source, and grouping the data source according to a satellite number and a signal frequency point to obtain grouped data;
step S2, segmenting each grouped data by utilizing a pseudo-range multi-path combined MP sequence, constructing a segmented linear model, and generating pseudo-range correction information tables of different types of satellites and different frequency points;
and step S3, searching pseudo-range correction node information in the pseudo-range correction information table according to the satellite type and the signal frequency point during time transmission, fitting, calculating a pseudo-range correction value, and correcting the pseudo-range observation value.
In the foregoing solution, in step S2, each packet of data is segmented, and segmentation is further performed according to the change of altitude, where an objective function of the segmented linear model is represented as:
Figure BDA0002538546100000021
s.t.fj(ej)-fj+1(ej)=0 (1)
in formula (1), j is 1,2, …, m-1, m represents the total number of segmentation nodes; i is 1,2 …, n represents the length of the MP sequence of the corresponding segment; f represents a piecewise function linearly related to the altitude; MP represents pseudorange multipath combination observations.
In the above scheme, the MP pseudorange and multipath combination observed value is obtained by combining a single-frequency pseudorange and a dual-frequency carrier phase, and is represented as:
Figure BDA0002538546100000031
in formula (2), i and j (i, j ≠ 1,2,3, i ≠ j) represent different frequencies; MP represents pseudo-range multi-path combination observed value; p and L represent pseudorange and carrier phase observations, respectively; f represents a carrier frequency; m represents a carrier multipath error; b comprises phase ambiguity and hardware delay deviation; representing the observed noise.
In the foregoing solution, the step S2 specifically includes the following steps:
step S21, cycle slip detection is carried out on the carrier phase observed value of each grouped data, and arc segment division is carried out on each grouped data according to whether cycle slip occurs or not;
step S22, calculating the average value of MP combination observation values of each arc segment of each grouped data, and subtracting the average value from all MP observation values;
and step S23, dividing MP combined observed values in the range of 5-85 degrees at preset intervals, and performing piecewise linear fitting to generate pseudo-range correction information tables of different types of satellites and different frequency points.
In the foregoing solution, in step S23, the predetermined interval is a 5 ° elevation angle, 17 parameters to be estimated are obtained by dividing MP combined observed values in a range of 5 ° to 85 °, and for any observed value [ e, MP ], the observation equation is expressed as:
Figure BDA0002538546100000032
in the formula (3), x represents the value of the parameter MP to be estimated corresponding to the node of the segment, and j represents the segment corresponding to the elevation angle e.
In the above scheme, the time transfer is the BDS CV common view of the Beidou satellite and/or the BDS PPP time transfer of the precise single-point positioning of the Beidou satellite.
In the foregoing solution, in step S3, the correcting the pseudorange observation value during time transfer specifically includes the following steps:
step S31, the Beidou broadcast ephemeris takes the frequency point B3 as reference, and the broadcast ephemeris group delay GD is corrected:
Figure BDA0002538546100000033
in the formula (4), f1And f2B1 and B2, respectivelyGD1And TGD2Respectively representing the group delay of B1 and B2 frequency points;
at step S32, the satellite co-view equation is expressed as:
Figure BDA0002538546100000034
in the formula (5), the reaction mixture is,
Figure BDA0002538546100000041
representing a dual-frequency ionosphere-free combination code,
Figure BDA0002538546100000042
which represents the coordinates of the satellite or satellites,
Figure BDA0002538546100000043
representing the antenna phase center coordinate calculated by using the dual-frequency ionosphere-free combined observation value, S representing the earth rotation effect, △ trelRepresenting relativistic effects △ ttropRepresenting tropospheric delay; GD is the broadcast ephemeris group delay.
In the above scheme, the pseudorange observation value is corrected and directly performed at the user side.
The invention has the following beneficial effects:
according to the method for correcting the satellite-end multipath error in the Beidou time transmission, the satellite-end multipath error correction is introduced in the Beidou time transmission, based on the observation data of the global MGEX observation station, a piecewise linear correction model is constructed by using the MP combined observation value and the altitude angle, a pseudo-range correction information table is obtained, the pseudo-range observation value of the Beidou second satellite is corrected, the pseudo-range is corrected systematically at a user terminal directly according to the altitude angle, the influence of the satellite-end multipath on the time transmission is weakened, the problem of the satellite-end multipath error in the Beidou second time transmission is solved, and the precision of the Beidou second time transmission is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for correcting errors of a Beidou second time transmission satellite terminal in the embodiment of the invention.
Detailed Description
The technical problems, aspects and advantages of the invention will be explained in detail below with reference to exemplary embodiments. The following exemplary embodiments are merely illustrative of the present invention and are not to be construed as limiting the invention. It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
The embodiment of the invention provides a method for correcting multipath errors of a Beidou second time transfer satellite terminal. The Beidou second satellite comprises 5 GEO (geostationary earth orbit) satellites (C01-C05), 7 IGSO (inclinedogenic geosynchronous satellite) satellites (C06-C10, C13 and C16) and 3 MEO (medium earth orbit) satellites (C11, C12 and C14). The Beidou second time transmission satellite-end multipath error correction method is realized by constructing a Beidou second time transmission satellite-end multipath error correction model.
Fig. 1 is a schematic flow chart of the Beidou second time transfer satellite-side multipath error correction method. As shown in fig. 1, the method for correcting the satellite-side multipath error in the second beidou time transfer is described by taking the satellite-side multipath error correction of the IGSO satellite and the MEO satellite in the second beidou as an example, and is also applicable to other satellites in the second beidou. The satellite-end multipath error correction method comprises the following steps:
step S1, selecting Multi-System observation network (meg) observation station data of a Global Navigation Satellite System (GNSS) Service organization (International GNSS Service, IGS) as a data source, and grouping the data source according to a Satellite number and a signal frequency point to obtain grouped data.
In this embodiment, because IGSO satellites mainly cover asia-pacific regions, observation angles of other regions are low, pseudo-range multi-path combination (MP) errors are large, which easily "pollutes" the overall observation value and reduces the model estimation accuracy. Therefore, when modeling the IGSO pseudo-range bias, the asia-pacific observation station should be selected to ensure a complete observation arc. For MEO satellites, monitoring stations distributed over a wide area should be selected as much as possible to ensure adequate observation over all ranges of elevation angles.
And step S2, segmenting each grouped data by using the MP sequence, constructing a segmented linear model, and generating pseudo-range correction information tables of different types of satellites and different frequency points.
Preferably, in this step, segmentation is performed according to the change of the altitude angle, and the objective function of the linear model after segmentation is represented as:
Figure BDA0002538546100000051
s.t.fj(ej)-fj+1(ej)=0 (1)
in formula (1), j is 1,2, …, m-1, m represents the total number of segmentation nodes; i is 1,2 …, n represents the length of the MP sequence of the corresponding segment; f represents a piecewise function linearly related to the altitude; MP represents the multipath combination observation.
When the Beidou second satellite is used for positioning, the measured distance contains the clock error of a user receiver and an atomic clock of the satellite and atmospheric refraction delay, but is not the true distance, and the true distance is called pseudo distance. In this embodiment, a pseudo-range multi-path combination (MP) is obtained by combining a single-frequency pseudo-range and a dual-frequency carrier phase, and is represented as:
Figure BDA0002538546100000061
in formula (2), i and j (i, j ≠ 1,2,3, i ≠ j) represent different frequencies; MP represents pseudo-range multi-path combination observed value; p and L represent pseudorange and carrier phase observations, respectively; f represents a carrier frequency; m represents a carrier multipath error; b comprises phase ambiguity and hardware delay deviation; representing the observed noise. The MP combined observation eliminates geometric distance, tropospheric delay, ionospheric delay, clock error and the like, and remains carrier phase ambiguity, hardware delay, multipath and observation noise. The hardware delay is kept unchanged in a short time, and the carrier multipath error is much smaller than the pseudo-range multipath error. So as long as no cycle slip occurs in the carrier phase, the MP combination mainly reflects the effects of pseudorange multipath effects. During actual calculation, cycle slip detection is firstly carried out on original observation data, then the MP mean value of each complete arc segment is calculated, and finally the MP mean value of the corresponding arc segment is subtracted from all MP values, so that the multipath error is obtained.
Further, step S2 specifically includes the following steps:
step S21, cycle slip detection is carried out on the carrier phase observed value of each grouped data, and arc segment division is carried out on each grouped data according to whether cycle slip occurs or not;
step S22, calculating the average value of MP combination observation values of each arc segment of each grouped data, and subtracting the average value from all MP observation values;
and step S23, dividing MP combined observed values in the range of 5-85 degrees at preset intervals, and performing piecewise linear fitting to generate pseudo-range correction information tables of different types of satellites and different frequency points.
Since the MP error is large when the altitude angle is low, which may affect the model estimation accuracy, in this embodiment, MP observation values of more than 5 ° are selected, and preferably, 5 ° are used as intervals. Furthermore, considering that only a few stations have MP observations close to 90 °, to prevent the low observations from making the linear parameter estimation of the segment inaccurate, the elevation angle of the MP sequence is less than 85 °. Preferably, the linear fitting is performed by using a whole least squares method.
In this step, when the elevation angle is 5 ° as an interval, the number of parameters to be estimated is 17, and for any observed value [ e, MP ], the observation equation is expressed as:
Figure BDA0002538546100000062
in the formula (3), x represents the value of the parameter MP to be estimated corresponding to the node of the segment, and j represents the segment corresponding to the elevation angle e.
Because the observation data quantity is large, in order to improve the calculation efficiency, a small number of observation values are selected to calculate initial MP node parameters and covariance matrixes thereof; then, using the node parameters and the covariance matrix as virtual observed quantities to construct a new observation equation; and finally, sequentially calculating final MP node parameters by adopting a sequential adjustment method.
Table 1 is an example of the pseudo-range correction information table of B1I/B2I based on IGSO/MEO generated in this embodiment. It should be noted that, in practical application, the node parameters are recalculated according to historical observation data.
TABLE 1 IGSO/MEO B1I/B2I pseudorange correction values
Figure BDA0002538546100000071
And step S3, searching pseudo-range correction node information in the pseudo-range correction information table according to the satellite type and the signal frequency point and fitting when the BDS CV and/or BDS PPP time is transferred, and correcting the pseudo-range observed value according to the fitting result.
In this step, taking the common view BDS CV of the big dipper satellite as an example, the satellite-end multipath error correction is explained, and the following calculation process is also applicable to the precise single-point positioning PPP time transfer of the big dipper, and specifically includes:
step S31, the Beidou broadcast ephemeris takes the frequency point B3 as reference, and the broadcast ephemeris group delay GD is corrected:
Figure BDA0002538546100000072
in the formula (4), f1And f2B1 and B2, respectivelyGD1And TGD2Respectively representing the group delay of B1 and B2 frequency points;
at step S32, the satellite co-view equation is expressed as:
Figure BDA0002538546100000081
in the formula (5), the reaction mixture is,
Figure BDA0002538546100000082
representing a dual-frequency ionosphere-free combination code,
Figure BDA0002538546100000083
which represents the coordinates of the satellite or satellites,
Figure BDA0002538546100000084
representing the antenna phase center coordinate calculated by using the dual-frequency ionosphere-free combined observation value, S representing the earth rotation effect, △ trelRepresenting relativistic effects △ ttropRepresenting tropospheric delay; GD is the broadcast ephemeris group delay.
In the step, the multipath error of the second Beidou satellite terminal is directly related to the altitude angle through the satellite coordinates
Figure BDA0002538546100000085
And receiver antenna phase center coordinates
Figure BDA0002538546100000086
And calculating the altitude angle of each satellite, performing linear interpolation according to the node information of the pseudo-range correction information table corresponding to the satellite type, the signal frequency point and the altitude angle, calculating a pseudo-range correction value, and correcting the pseudo-range observation value.
According to the technical scheme, satellite multipath error correction is introduced into Beidou time transmission, based on global MGEX observation station observation data, a piecewise linear correction model is constructed by using MP combined observation values and altitude angles, a pseudo range correction information table is obtained, the pseudo range observation values of the Beidou second satellite are corrected, the pseudo ranges are corrected systematically directly at a user terminal according to the altitude angles, the influence of satellite multipath on the time transmission is weakened, the problem of satellite multipath error transmission of the Beidou second time transmission is solved, and the precision of the Beidou second time transmission is improved.
While the foregoing is directed to the preferred embodiment of the present invention, it is understood that the invention is not limited to the exemplary embodiments disclosed, but is made merely for the purpose of providing those skilled in the relevant art with a comprehensive understanding of the specific details of the invention. It will be apparent to those skilled in the art that various modifications and adaptations of the present invention can be made without departing from the principles of the invention and the scope of the invention is to be determined by the claims.

Claims (8)

1. A Beidou second time transfer satellite-side multipath error correction method is characterized by comprising the following steps:
step S1, selecting data of a multi-system observation network MGEX observation station of the International Global navigation satellite service organization as a data source, and grouping the data source according to a satellite number and a signal frequency point to obtain grouped data;
step S2, segmenting each grouped data by utilizing a pseudo-range multi-path combined MP sequence, constructing a segmented linear model, and generating pseudo-range correction information tables of different types of satellites and different frequency points;
and step S3, searching pseudo-range correction node information in the pseudo-range correction information table according to the satellite type and the signal frequency point during time transmission, fitting, calculating a pseudo-range correction value, and correcting the pseudo-range observation value.
2. The Beidou second time-transfer satellite-side multipath error correction method according to claim 1, wherein in the step S2, each packet data is segmented and further segmented according to altitude angle change, and an objective function of a segmented linear model is represented as:
Figure FDA0002538546090000011
s.t.fj(ej)-fj+1(ej)=0 (1)
in formula (1), j is 1,2, …, m-1, m represents the total number of segmentation nodes; i is 1,2 …, n represents the length of the MP sequence of the corresponding segment; f represents a piecewise function linearly related to the altitude; MP represents pseudorange multipath combination observations.
3. The Beidou second time-transfer satellite-terminal multipath error correction method according to claim 2, characterized in that the MP pseudorange multipath combination observed value is obtained by combining a single frequency pseudorange and a dual frequency carrier phase, and is expressed as:
Figure FDA0002538546090000012
in formula (2), i and j (i, j ≠ 1,2,3, i ≠ j) represent different frequencies; MP represents pseudo-range multi-path combination observed value; p and L represent pseudorange and carrier phase observations, respectively; f represents a carrier frequency; m represents a carrier multipath error; b comprises phase ambiguity and hardware delay deviation; representing the observed noise.
4. The Beidou second time-transfer satellite-side multipath error correction method according to any one of claims 1 to 3, wherein the step S2 specifically comprises the steps of:
step S21, cycle slip detection is carried out on the carrier phase observed value of each grouped data, and arc segment division is carried out on each grouped data according to whether cycle slip occurs or not;
step S22, calculating the average value of MP combination observation values of each arc segment of each grouped data, and subtracting the average value from all MP observation values;
and step S23, dividing MP combined observed values in the range of 5-85 degrees at preset intervals, and performing piecewise linear fitting to generate pseudo-range correction information tables of different types of satellites and different frequency points.
5. The Beidou second time-transfer satellite-side multipath error correction method according to claim 4, wherein the predetermined interval in the step S23 is a 5 ° elevation angle, 17 parameters to be estimated are obtained by dividing MP combined observation values in a range of 5 ° to 85 °, and for any observation value [ e, MP ], the observation equation is expressed as:
Figure FDA0002538546090000021
in the formula (3), x represents the value of the parameter MP to be estimated corresponding to the node of the segment, and j represents the segment corresponding to the elevation angle e.
6. The Beidou second time transfer satellite-side multipath error correction method according to claim 5, wherein the time transfer is Beidou satellite common view BDS CV and/or Beidou satellite precise single point positioning BDS PPP time transfer.
7. The method for correcting the multipath error at the Beidou second time-transfer satellite end of claim 5, wherein in the step S3, the pseudo-range observed value is corrected during the time transfer, and the method specifically comprises the following steps:
step S31, the Beidou broadcast ephemeris takes the frequency point B3 as reference, and the broadcast ephemeris group delay GD is corrected:
Figure FDA0002538546090000022
in the formula (4), f1And f2B1 and B2, respectivelyGD1And TGD2Respectively representing the group delay of B1 and B2 frequency points;
at step S32, the satellite co-view equation is expressed as:
Figure FDA0002538546090000023
in the formula (5), the reaction mixture is,
Figure FDA0002538546090000024
representing a dual-frequency ionosphere-free combination code,
Figure FDA0002538546090000025
which represents the coordinates of the satellite or satellites,
Figure FDA0002538546090000026
representing the antenna phase center coordinate calculated by using the dual-frequency ionosphere-free combined observation value, S representing the earth rotation effect, △ trelRepresenting relativistic effects △ ttropRepresenting tropospheric delay; GD is the broadcast ephemeris group delay.
8. The method of claim 7, wherein the correction of the pseudorange observations is performed directly at the user end.
CN202010539934.8A 2020-06-15 2020-06-15 Beidou second-order time transfer satellite-side multipath error correction method Pending CN111766616A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010539934.8A CN111766616A (en) 2020-06-15 2020-06-15 Beidou second-order time transfer satellite-side multipath error correction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010539934.8A CN111766616A (en) 2020-06-15 2020-06-15 Beidou second-order time transfer satellite-side multipath error correction method

Publications (1)

Publication Number Publication Date
CN111766616A true CN111766616A (en) 2020-10-13

Family

ID=72721252

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010539934.8A Pending CN111766616A (en) 2020-06-15 2020-06-15 Beidou second-order time transfer satellite-side multipath error correction method

Country Status (1)

Country Link
CN (1) CN111766616A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230254A (en) * 2020-10-29 2021-01-15 新乡医学院三全学院 Method and device for correcting multipath errors of GPS carrier phase
CN112924992A (en) * 2021-01-25 2021-06-08 中国科学院国家授时中心 GEO orbit precision evaluation method and device, electronic equipment and storage medium
CN114355421A (en) * 2021-12-22 2022-04-15 杭州电子科技大学 Flood detection method based on Beidou satellite L4 and CMC combined observation value
CN114814903A (en) * 2022-04-26 2022-07-29 中国计量科学研究院 Common-view data processing method based on Beidou third-generation navigation satellite
CN115308781A (en) * 2022-09-02 2022-11-08 中国人民解放军国防科技大学 BDGIM assistance-based phase smoothing pseudorange high-precision time transfer method
CN116299618A (en) * 2023-03-24 2023-06-23 中国科学院精密测量科学与技术创新研究院 Carrier phase satellite common view time transfer method based on PPP (point-to-point protocol) calculation parameters

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009139281A (en) * 2007-12-07 2009-06-25 Xanavi Informatics Corp Positioning method
CN105607089A (en) * 2016-01-18 2016-05-25 西安测绘研究所 Pseudo-range-assisted-carrier-phase-based beidou common-view time transmission method
CN107544082A (en) * 2016-06-29 2018-01-05 武汉大学 The step modeling of Big Dipper IGSO/MEO satellite pseudorange codes deviation one
CN108919634A (en) * 2018-08-13 2018-11-30 中国科学院国家授时中心 A kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou
CN110146906A (en) * 2019-05-23 2019-08-20 中国科学院国家授时中心 Remote time transmission method based on single poor carrier phase observation data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009139281A (en) * 2007-12-07 2009-06-25 Xanavi Informatics Corp Positioning method
CN105607089A (en) * 2016-01-18 2016-05-25 西安测绘研究所 Pseudo-range-assisted-carrier-phase-based beidou common-view time transmission method
CN107544082A (en) * 2016-06-29 2018-01-05 武汉大学 The step modeling of Big Dipper IGSO/MEO satellite pseudorange codes deviation one
CN108919634A (en) * 2018-08-13 2018-11-30 中国科学院国家授时中心 A kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou
CN110146906A (en) * 2019-05-23 2019-08-20 中国科学院国家授时中心 Remote time transmission method based on single poor carrier phase observation data

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
GUOJUN LI ET AL.: "Analysis on BDS Satellite-Induced Multipath and Its Impact on Time Transfer", CHINA SATELLITE NAVIGATION CONFERENCE (CSNC) 2019 PROCEEDINGS, pages 518 - 530 *
孔祥元 等: "控制测量学 下册", 30 November 2015, 武汉大学出版社, pages: 310 - 314 *
常军;党海龙;李昕;: "北斗卫星多路径系统偏差改正的研究", 测绘地理信息, no. 01 *
李昕;张小红;曾琪;潘林;朱锋;: "北斗卫星伪距偏差模型估计及其对精密定位的影响", 武汉大学学报(信息科学版), no. 10 *
汪捷;何锡扬;: "北斗卫星伪距多路径系统偏差的修正方法及其对单频PPP的影响分析", 测绘学报, no. 07 *
王阿昊;陈俊平;张益泽;王解先;王彬;: "基于分区综合改正数的北斗卫星导航系统和GPS组合的动态精密单点定位", 同济大学学报(自然科学版), no. 03 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230254A (en) * 2020-10-29 2021-01-15 新乡医学院三全学院 Method and device for correcting multipath errors of GPS carrier phase
CN112230254B (en) * 2020-10-29 2024-05-14 新乡医学院三全学院 Correction method and device for GPS carrier phase multipath error
CN112924992A (en) * 2021-01-25 2021-06-08 中国科学院国家授时中心 GEO orbit precision evaluation method and device, electronic equipment and storage medium
CN112924992B (en) * 2021-01-25 2022-11-04 中国科学院国家授时中心 GEO (geosynchronous orbit) track precision evaluation method and device, electronic equipment and storage medium
CN114355421A (en) * 2021-12-22 2022-04-15 杭州电子科技大学 Flood detection method based on Beidou satellite L4 and CMC combined observation value
CN114355421B (en) * 2021-12-22 2023-08-18 杭州电子科技大学 Flood detection method based on Beidou satellite L4 and CMC combined observation value
CN114814903A (en) * 2022-04-26 2022-07-29 中国计量科学研究院 Common-view data processing method based on Beidou third-generation navigation satellite
CN114814903B (en) * 2022-04-26 2022-12-30 中国计量科学研究院 Common-view data processing method based on Beidou third-generation navigation satellite
CN115308781A (en) * 2022-09-02 2022-11-08 中国人民解放军国防科技大学 BDGIM assistance-based phase smoothing pseudorange high-precision time transfer method
CN116299618A (en) * 2023-03-24 2023-06-23 中国科学院精密测量科学与技术创新研究院 Carrier phase satellite common view time transfer method based on PPP (point-to-point protocol) calculation parameters
CN116299618B (en) * 2023-03-24 2024-03-19 中国科学院精密测量科学与技术创新研究院 Carrier phase satellite common view time transfer method based on PPP (point-to-point protocol) calculation parameters

Similar Documents

Publication Publication Date Title
CN111766616A (en) Beidou second-order time transfer satellite-side multipath error correction method
CN111190203B (en) PPP-RTK positioning method based on non-difference observation model and rank deficiency elimination
CN109581452B (en) GNSS reference station carrier phase integer ambiguity resolution method
US10281587B2 (en) Navigation satellite system positioning involving the generation of correction information
CN111045034B (en) GNSS multi-system real-time precise time transfer method and system based on broadcast ephemeris
US10078140B2 (en) Navigation satellite system positioning involving the generation of advanced correction information
CN109520512A (en) A kind of precision orbit determination method and device
CN102998681B (en) A kind of high-frequency clock error estimation method of satellite navigation system
CN114966760B (en) Ionosphere weighted non-differential non-combination PPP-RTK technology implementation method
CN111694030A (en) BDS local difference method and system based on grid virtual observation value
CN114994727B (en) Equipment for realizing high-precision time calibration and satellite positioning
CN115933356B (en) High-precision time synchronization system and method for virtual atomic clock
CN116148909B (en) Multi-frequency multi-mode non-combination precise single-point positioning instantaneous narrow-lane ambiguity fixing method
CN115993617B (en) GNSS system time deviation monitoring method
CN115225245B (en) Non-differential non-combined PPP-RTK long baseline real-time transfer method
Dai et al. Improving real-time clock estimation with undifferenced ambiguity fixing
CN113325446A (en) Multi-mode common-frequency GNSS carrier phase time transfer method and system
CN116047555A (en) Positioning method and system based on Beidou PPP-B2B
CN117111112A (en) GNSS and LNSS data fusion processing method and system
Li et al. Assessment and analysis of the four-satellite QZSS precise point positioning and the integrated data processing with GPS
CN116125513A (en) GNSS ambiguity fixed solution checking method and system
CN111766614B (en) Sub-nanosecond real-time comparison method for joint networking of time keeping laboratory
CN115308781A (en) BDGIM assistance-based phase smoothing pseudorange high-precision time transfer method
CN112485813A (en) Method and system for correcting frequency offset of non-combined ranging codes between GLONASS measuring stations
CN110542910B (en) Method, device and equipment for correcting hardware delay deviation in satellite positioning

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