CN110824520B - Method and system for calculating ionospheric phase flicker index in real time - Google Patents

Method and system for calculating ionospheric phase flicker index in real time Download PDF

Info

Publication number
CN110824520B
CN110824520B CN201911155320.3A CN201911155320A CN110824520B CN 110824520 B CN110824520 B CN 110824520B CN 201911155320 A CN201911155320 A CN 201911155320A CN 110824520 B CN110824520 B CN 110824520B
Authority
CN
China
Prior art keywords
carrier
phase
doppler
calculating
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911155320.3A
Other languages
Chinese (zh)
Other versions
CN110824520A (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.)
National Space Science Center of CAS
Original Assignee
National Space Science Center of CAS
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 National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN201911155320.3A priority Critical patent/CN110824520B/en
Publication of CN110824520A publication Critical patent/CN110824520A/en
Application granted granted Critical
Publication of CN110824520B publication Critical patent/CN110824520B/en
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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections

Abstract

The invention discloses a method and a system for calculating an ionospheric phase flicker index in real time, wherein the method comprises the following steps: obtaining a measured value of a carrier phase from a receiver, and calculating a Doppler frequency measured value of a signal; based on the Doppler frequency measurement value of the signal, obtaining a carrier phase for removing phase change generated by relative motion of the satellite and the receiver; filtering the carrier phase by using a six-order Butterworth high-pass filter to remove satellite clock error, receiver clock error, flow error and fitted residual error; obtaining a filtered carrier phase; and calculating the ionospheric phase flicker index based on the filtered carrier phase and by combining the thermal noise correction value. The method can provide high-precision ionospheric phase scintillation index measurement values.

Description

Method and system for calculating ionospheric phase flicker index in real time
Technical Field
The invention relates to the field of satellite communication and navigation positioning, in particular to a method and a system for calculating an ionospheric phase scintillation index in real time.
Background
Ionospheric scintillation refers to the electromagnetic wave amplitude fading and phase jitter caused by the uneven structure of the ionosphere, which brings serious harm to satellite communication, navigation positioning and the like. Therefore, the real-time monitoring of ionospheric scintillation is of great significance. Ionospheric phase scintillation is one of the important observables to measure ionospheric scintillation. Besides the flicker phase, the carrier phase measured by the receiver also contains the change caused by the relative motion of the satellite receiver, the clock error of the receiver and the satellite, the phase jitter caused by the troposphere and the thermal noise.
The phase trend caused by the satellite motion needs to be removed when calculating the ionospheric scintillation phase. A common method is to fit a 4 th order polynomial to the phase trend and subtract, but this method is only applied in post-hoc data processing.
Disclosure of Invention
The invention aims to overcome the technical defects and provides a method for calculating the ionospheric phase scintillation index in real time. The method can effectively remove the influence of other factors and accurately measure the ionospheric scintillation phase.
In order to achieve the above object, the present invention provides a method for calculating an ionospheric phase flicker index in real time, the method comprising:
obtaining a measured value of a carrier phase from a receiver, and calculating a Doppler frequency measured value of a signal;
based on the Doppler frequency measurement value of the signal, obtaining a carrier phase for removing phase change generated by relative motion of the satellite and the receiver;
filtering the carrier phase by using a six-order Butterworth high-pass filter to remove satellite clock error, receiver clock error, troposphere error and fitted residual error; obtaining a filtered carrier phase;
and calculating the ionospheric phase flicker index based on the filtered carrier phase and by combining the thermal noise correction value.
As an improvement of the above method, the carrier phase is obtained by removing phase changes caused by relative motion between the satellite and the receiver based on doppler frequency measurements of the signals; the method specifically comprises the following steps:
calculating observation quantity of polynomial fitting according to the time recorded by the Doppler frequency and the Doppler frequency;
solving coefficients of the polynomial by using a least square method;
calculating the fitted carrier Doppler size according to the current time;
subtracting the fitted carrier Doppler and the measured carrier Doppler to obtain a carrier frequency value with the movement trend removed;
and integrating the carrier frequency values without the motion tendency to obtain carrier phases.
As an improvement of the above method, the calculating of the observation of polynomial fitting from the time of doppler frequency recording and the doppler frequency; the method specifically comprises the following steps:
establishing a relation between the time of Doppler frequency recording and the Doppler frequency:
T(N)P(N)＝C(N)
wherein T (N) is a time matrix with a variable at the power of time; p (N) is the coefficient of the polynomial; c (N) is a carrier Doppler vector; n represents the number of data used to perform the fitting;
the time matrix T (N) is:
Figure GDA0003204828940000021
wherein t (i) represents the time corresponding to the ith sampling point, and i is more than or equal to 1 and less than or equal to N;
the carrier doppler vector c (n) is:
Figure GDA0003204828940000022
wherein, carph (i) represents the carrier phase doppler corresponding to the ith sampling point;
the polynomial coefficient vector P (N) is:
Figure GDA0003204828940000023
Pithe ith coefficient of the polynomial.
As an improvement of the above method, the coefficients of the polynomial are found by a least square method; the method specifically comprises the following steps:
the curve equation is normalized by multiplying the two sides of T (N) -P (N) -C (N) by the transpose of T (N) ((N))TObtaining:
U(N)P(N)＝D(N)
wherein U (N) is:
U(N)＝T(N)TT(N)
wherein D (N) is:
D(N)＝T(N)TC(N)
the increment Δ U (n +1) from U (n) to U (n +1) is:
Figure GDA0003204828940000031
wherein N is more than or equal to 1 and less than or equal to N-1;
the recurrence formula for U (n) is therefore:
U(n+1)＝U(n)+ΔU(n+1)
calculating by using the recursion formula and the initial value U (1) to obtain U (N);
the increments of D (n) to D (n +1) are:
Figure GDA0003204828940000032
the recurrence formula for D (N) is therefore:
D(N+1)＝D(N)+ΔD(N+1)
calculating by using the recursion formula and the initial value D (1) to obtain D (N);
p (n) is obtained by the least square method based on u (n) p (n) ═ d (n).
As an improvement of the above method, the fitted carrier doppler magnitude is calculated according to the current time; the method specifically comprises the following steps:
the fitted carrier doppler magnitudes are:
Figure GDA0003204828940000033
wherein t is the current time.
As an improvement of the above method, calculating an ionospheric phase flicker index based on the filtered carrier phase in combination with a thermal noise correction value specifically includes:
phase sigma due to thermal noisepllComprises the following steps:
Figure GDA0003204828940000034
wherein, BnIs the loop bandwidth, C/N, of the receiver0Representing the carrier-to-noise ratio, T is the loop integration time;
ionospheric phase flicker index corrected for thermal noise error
Figure GDA0003204828940000035
Comprises the following steps:
Figure GDA0003204828940000041
wherein the content of the first and second substances,
Figure GDA0003204828940000042
is the filtered carrier phase.
The invention also provides a system for calculating the ionospheric phase scintillation index in real time, which comprises:
the Doppler frequency measurement value calculation module is used for acquiring the measurement value of the carrier phase from the receiver and calculating the Doppler frequency measurement value of the signal;
the carrier phase correction module is used for obtaining a carrier phase for removing phase change generated by relative motion of the satellite and the receiver based on the Doppler frequency measurement value of the signal;
the carrier phase filtering module is used for filtering the carrier phase by using a six-order Butterworth high-pass filter to remove satellite clock error, receiver clock error, troposphere error and fitted residual error; obtaining a filtered carrier phase;
and the ionospheric phase flicker index calculation module is used for calculating the ionospheric phase flicker index based on the filtered carrier phase and combined with the thermal noise correction value.
The invention has the advantages that:
1. the method of the invention removes the phase trend by using a polynomial fitting method and can run on the GNSS ionosphere scintillator in real time, and in addition, the method also considers the measurement error brought by the loop and corrects the measurement error;
2. the method has stable operation and accurate result, and can provide high-precision ionospheric phase scintillation index measurement values.
Drawings
Fig. 1 is a flowchart of a method for calculating an ionospheric phase scintillation index in real time according to the present invention.
Detailed Description
The technical solution of the present invention will be described in detail below with reference to the accompanying drawings.
In order to calculate the phase jitter caused by ionospheric flicker, the phase jitter caused by other factors needs to be filtered. The phase change due to the relative motion of the satellite and the receiver is removed first. The phase change caused by the relative motion between the satellite and the receiver can be calculated by fitting a polynomial to the carrier phase, and subtracted from the measured carrier phase. In order to reduce the calculation amount, the method fits the carrier Doppler of the signal through a second-order polynomial, and then subtracts the fitted carrier Doppler from the measured carrier Doppler so as to obtain the carrier Doppler with the movement trend removed. Therefore, the carrier phase to remove the relative motion between the satellite and the receiver can be calculated by the following formula:
Figure GDA0003204828940000043
coefficient p of polynomial in the formulaiCalculated by the least square method, see formula (2), t is time,
Figure GDA0003204828940000051
is the carrier doppler measured by the receiver;
UP＝D (2)
where U is the input time series normalization matrix, P is the column vector of polynomial coefficients, and D is the product of carrier doppler and time series. The polynomial coefficient P can be found by multiplying both sides of the equation by the inverse of the U matrix.
The satellite clock error, the receiver clock error and the phase change caused by neutral atmosphere can be filtered by a 6-order Butterworth high-pass filter, and the cut-off frequency of the filter is 0.1 Hz. The parameters of the filter are shown in table 1,
table 1: butterworth filter coefficient
Figure GDA0003204828940000052
The phase due to thermal noise can be estimated using the carrier-to-noise ratio, the loop bandwidth and the integration time, and is given by the following formula,
Figure GDA0003204828940000053
in the formula BnIs the loop bandwidth, C/N, of the receiver0Representing the carrier-to-noise ratio, T is the loop integration time; sigmapllPhase due to thermal noise; the error caused by the thermal noise can be corrected by using the estimated noise jitter and combining the following formula:
Figure GDA0003204828940000054
Figure GDA0003204828940000055
is the ionospheric phase flicker index corrected for errors.
As shown in fig. 1, embodiment 1 of the present invention provides a method for calculating an ionospheric phase scintillation index in real time, including:
step 1) obtaining a measured value of a carrier phase from a receiver;
step 2) carrying out differential operation on the carrier phase to obtain a Doppler frequency measurement value of the signal;
step 3) calculating the observed quantity of polynomial fitting according to the time recorded by the Doppler frequency and the Doppler frequency;
the specific calculation method is as follows;
T(N)P(N)＝C(N) (5)
wherein T (N) is a matrix with a variable at the power of time;
p (N) is the coefficient of the polynomial;
c (N) is the measured carrier Doppler;
n represents the number of data used to perform the fit.
The time matrix T (N) is:
Figure GDA0003204828940000061
where t (i) represents the time corresponding to the ith sampling point,
the carrier measurement vector c (n) is:
Figure GDA0003204828940000062
carph (i) represents the sampled value of the ith carrier phase doppler;
the polynomial coefficient vector P (N) is:
Figure GDA0003204828940000063
Pirepresents the ith coefficient of the polynomial;
the polynomial coefficients are solved by using the least square method, firstly, the curve equation is normalized, namely coordinates on two sides of the equation are multiplied by the transposition of T (N),
U(N)P(N)＝D(N) (9)
wherein U (N) is a group of,
U(N)＝T(N)TT(N) (10)
wherein D (N) is a group of,
D(N)＝T(N)TC(N) (11)
since the receiver resources are limited and all past values cannot be stored, the recursion formulas of u (n) and d (n) need to be found, and u (n) and d (n) need to be found in an iterative manner.
The increments of U (N) to U (N +1) are:
Figure GDA0003204828940000071
the recurrence formula of u (n) is,
U(N+1)＝U(N)+ΔU(N+1) (13)
the increments of D (N) to D (N +1) are:
Figure GDA0003204828940000072
the recurrence formula for D (N) is therefore:
D(N+1)＝D(N)+ΔD(N+1) (15)
step 4) solving coefficients of the polynomial by using a least square method and calculating the fitted carrier Doppler size according to the current time;
step 5) subtracting the fitted carrier Doppler and the measured carrier Doppler to obtain a carrier frequency value without the movement trend;
step 6) integrating the carrier frequency values without the motion tendency to obtain carrier phases;
step 7) filtering the carrier phase by using a six-order Butterworth high-pass filter to remove satellite clock error, receiver clock error, troposphere error and fitted residual error;
and 8) calculating the ionospheric phase flicker index by combining the thermal noise correction value.
Example 2
Embodiment 2 of the present invention provides a system for calculating an ionospheric phase flicker index in real time, the system including:
the Doppler frequency measurement value calculation module is used for acquiring the measurement value of the carrier phase from the receiver and calculating the Doppler frequency measurement value of the signal;
the carrier phase correction module is used for obtaining a carrier phase for removing phase change generated by relative motion of the satellite and the receiver based on the Doppler frequency measurement value of the signal;
the carrier phase filtering module is used for filtering the carrier phase by using a six-order Butterworth high-pass filter to remove satellite clock error, receiver clock error, troposphere error and fitted residual error; obtaining a filtered carrier phase;
and the ionospheric phase flicker index calculation module is used for calculating the ionospheric phase flicker index based on the filtered carrier phase and combined with the thermal noise correction value.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention and are not limited. Although the present invention has been described in detail with reference to the embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (5)

1. A method of calculating ionospheric phase scintillation indices in real-time, the method comprising:
obtaining a measured value of a carrier phase from a receiver, and calculating a Doppler frequency measured value of a signal;
based on the Doppler frequency measurement value of the signal, obtaining a carrier phase for removing phase change generated by relative motion of the satellite and the receiver;
filtering the carrier phase by using a six-order Butterworth high-pass filter to remove satellite clock error, receiver clock error, troposphere error and fitted residual error; obtaining a filtered carrier phase;
calculating an ionospheric phase flicker index based on the filtered carrier phase in combination with the thermal noise correction value;
the carrier phase which removes the phase change generated by the relative motion of the satellite and the receiver is obtained based on the Doppler frequency measurement value of the signal; the method specifically comprises the following steps:
calculating observation quantity of polynomial fitting according to the time recorded by the Doppler frequency and the Doppler frequency;
solving coefficients of the polynomial by using a least square method;
calculating the fitted carrier Doppler size according to the current time;
subtracting the fitted carrier Doppler and the measured carrier Doppler to obtain a carrier frequency value with the movement trend removed;
integrating the carrier frequency values without the motion tendency to obtain carrier phases;
the observed quantity of polynomial fitting is calculated according to the time recorded by the Doppler frequency and the Doppler frequency; the method specifically comprises the following steps:
establishing a relation between the time of Doppler frequency recording and the Doppler frequency:
T(N)P(N)＝C(N)
wherein T (N) is a time matrix with a variable at the power of time; p (N) is the coefficient of the polynomial; c (N) is a carrier Doppler vector; n represents the number of data used to perform the fitting;
the time matrix T (N) is:
Figure FDA0003204828930000011
wherein t (i) represents the time corresponding to the ith sampling point, and i is more than or equal to 1 and less than or equal to N;
the carrier doppler vector c (n) is:
Figure FDA0003204828930000021
wherein, carph (i) represents the carrier phase doppler corresponding to the ith sampling point;
the polynomial coefficient vector P (N) is:
Figure FDA0003204828930000022
2. the method according to claim 1, wherein the coefficients of the polynomial are found by a least squares method; the method specifically comprises the following steps:
normalizing the curve equation by multiplying the two sides of T (N) P (N) C (N) by the transpose of T (N) ((N))TObtaining:
U(N)P(N)＝D(N)
wherein U (N) is:
U(N)＝T(N)TT(N)
wherein D (N) is:
D(N)＝T(N)TC(N)
the increment Δ U (n +1) from U (n) to U (n +1) is:
Figure FDA0003204828930000023
wherein N is more than or equal to 1 and less than or equal to N-1;
the recurrence formula for U (n) is therefore:
U(n+1)＝U(n)+ΔU(n+1)
calculating by using the recursion formula and the initial value U (1) to obtain U (N);
the increments of D (n) to D (n +1) are:
Figure FDA0003204828930000024
the recurrence formula for D (N) is therefore:
D(N+1)＝D(N)+ΔD(N+1)
calculating by using the recursion formula and the initial value D (1) to obtain D (N);
p (n) is obtained by the least square method based on u (n) p (n) ═ d (n).
3. The method of claim 2, wherein the fitted carrier doppler magnitude is calculated from a current time; the method specifically comprises the following steps:
the fitted carrier doppler magnitudes are:
Figure FDA0003204828930000031
wherein t is the current time.
4. The method of claim 3, wherein calculating the ionospheric phase flicker index based on the filtered carrier phase in combination with the thermal noise correction value comprises:
phase sigma due to thermal noisepllComprises the following steps:
Figure FDA0003204828930000032
wherein, BnIs the loop bandwidth, C/N, of the receiver0Representing the carrier-to-noise ratio, T is the loop integration time;
ionospheric phase flicker index corrected for thermal noise error
Figure FDA0003204828930000033
Comprises the following steps:
Figure FDA0003204828930000034
wherein the content of the first and second substances,
Figure FDA0003204828930000035
is the filtered carrier phase.
5. A system for calculating an ionospheric phase scintillation index in real time, the system comprising:
the Doppler frequency measurement value calculation module is used for acquiring the measurement value of the carrier phase from the receiver and calculating the Doppler frequency measurement value of the signal;
the carrier phase correction module is used for obtaining a carrier phase for removing phase change generated by relative motion of the satellite and the receiver based on the Doppler frequency measurement value of the signal;
the carrier phase filtering module is used for filtering the carrier phase by using a six-order Butterworth high-pass filter to remove satellite clock error, receiver clock error, troposphere error and fitted residual error; obtaining a filtered carrier phase;
the ionospheric phase flicker index calculation module is used for calculating the ionospheric phase flicker index based on the filtered carrier phase and combined with the thermal noise correction value;
the specific processing of the carrier phase correction module includes:
calculating observation quantity of polynomial fitting according to the time recorded by the Doppler frequency and the Doppler frequency;
solving coefficients of the polynomial by using a least square method;
calculating the fitted carrier Doppler size according to the current time;
subtracting the fitted carrier Doppler and the measured carrier Doppler to obtain a carrier frequency value with the movement trend removed;
integrating the carrier frequency values without the motion tendency to obtain carrier phases;
the observed quantity of polynomial fitting is calculated according to the time recorded by the Doppler frequency and the Doppler frequency; the method specifically comprises the following steps:
establishing a relation between the time of Doppler frequency recording and the Doppler frequency:
T(N)P(N)＝C(N)
wherein T (N) is a time matrix with a variable at the power of time; p (N) is the coefficient of the polynomial; c (N) is a carrier Doppler vector; n represents the number of data used to perform the fitting;
the time matrix T (N) is:
Figure FDA0003204828930000041
wherein t (i) represents the time corresponding to the ith sampling point, and i is more than or equal to 1 and less than or equal to N;
the carrier doppler vector c (n) is:
Figure FDA0003204828930000042
wherein, carph (i) represents the carrier phase doppler corresponding to the ith sampling point;
the polynomial coefficient vector P (N) is:
Figure FDA0003204828930000043
CN201911155320.3A 2019-11-22 2019-11-22 Method and system for calculating ionospheric phase flicker index in real time Active CN110824520B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911155320.3A CN110824520B (en) 2019-11-22 2019-11-22 Method and system for calculating ionospheric phase flicker index in real time

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911155320.3A CN110824520B (en) 2019-11-22 2019-11-22 Method and system for calculating ionospheric phase flicker index in real time

Publications (2)

Publication Number Publication Date
CN110824520A CN110824520A (en) 2020-02-21
CN110824520B true CN110824520B (en) 2021-11-09

Family

ID=69558458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911155320.3A Active CN110824520B (en) 2019-11-22 2019-11-22 Method and system for calculating ionospheric phase flicker index in real time

Country Status (1)

Country Link
CN (1) CN110824520B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111323796B (en) * 2020-03-18 2021-11-09 中国科学院国家空间科学中心 GNSS receiver high-sampling clock error resolving method
CN111983654B (en) * 2020-08-24 2021-03-19 中国矿业大学 Method for constructing ionosphere phase scintillation factor in arctic region based on GNSS
CN113031036B (en) * 2021-03-01 2021-09-24 中国矿业大学 Ionosphere phase flicker factor construction method based on GNSS 30s sampling frequency data
WO2023224508A1 (en) * 2022-05-16 2023-11-23 Limited Liability Company "Topcon Positioning Systems" Method and apparatuses of measuring phase ionosphere scintillations

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012128979A2 (en) * 2011-03-22 2012-09-27 Trimble Navigation Limited Gnss signal processing with ionospheric bridging for reconvergence
CN106134462B (en) * 2011-12-06 2014-07-09 北京遥测技术研究所 A kind of spaceborne ionospheric scintillation monitoring device that utilizes gps signal phase place to change
CN107991690A (en) * 2017-06-28 2018-05-04 江苏省测绘工程院 A kind of computational methods of the scintillation index based on Beidou navigation
CN108983258A (en) * 2018-05-30 2018-12-11 南京信息工程大学 A kind of GNSS ionospheric scintillation and TEC monitoring device
CN109188474A (en) * 2018-08-14 2019-01-11 中国科学院光电研究院 A kind of ionospheric scintillation monitoring system for significantly ionospheric scintillation
CN109444935A (en) * 2018-10-17 2019-03-08 桂林电子科技大学 A kind of Doppler's detection and reparation for cycle slips method of low sampling rate

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012128979A2 (en) * 2011-03-22 2012-09-27 Trimble Navigation Limited Gnss signal processing with ionospheric bridging for reconvergence
CN106134462B (en) * 2011-12-06 2014-07-09 北京遥测技术研究所 A kind of spaceborne ionospheric scintillation monitoring device that utilizes gps signal phase place to change
CN107991690A (en) * 2017-06-28 2018-05-04 江苏省测绘工程院 A kind of computational methods of the scintillation index based on Beidou navigation
CN108983258A (en) * 2018-05-30 2018-12-11 南京信息工程大学 A kind of GNSS ionospheric scintillation and TEC monitoring device
CN109188474A (en) * 2018-08-14 2019-01-11 中国科学院光电研究院 A kind of ionospheric scintillation monitoring system for significantly ionospheric scintillation
CN109444935A (en) * 2018-10-17 2019-03-08 桂林电子科技大学 A kind of Doppler's detection and reparation for cycle slips method of low sampling rate

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
卫星导航系统电离层闪烁监测技术研究;徐昌元;《中国优秀硕士学位论文全文数据库信息科技辑》;20190115;第20-25页、第33-34页 *

Also Published As

Publication number Publication date
CN110824520A (en) 2020-02-21

Similar Documents

Publication Publication Date Title
CN110824520B (en) Method and system for calculating ionospheric phase flicker index in real time
JP3455266B2 (en) GPS code measurement value smoothing method, GPS receiver differential position determination method and apparatus
CN107678050B (en) GLONASS phase inter-frequency deviation real-time tracking and precise estimation method based on particle filtering
CN109581856B (en) Time synchronization and time keeping method based on high-performance crystal oscillator frequency calibration
CN109507704B (en) Double-satellite positioning frequency difference estimation method based on mutual ambiguity function
US9621288B2 (en) Method and apparatus for estimating the current signal-to-noise ratio
CN110307780A (en) PGC phase demodulating real-time error compensation method based on iterative calculation
TW201003106A (en) Digital front end for a satellite navigation receiver
US20140354475A1 (en) Mitigation of spurious signals in gnss receivers
CN109613568A (en) A kind of satellite navigation receiver frequency discriminator measurement noise removing method
JP6203608B2 (en) GLONASS receiver
CN114050867B (en) Signal compensation method, device, equipment and storage medium
TWI384771B (en) Signal receiver and method for estimating residual doppler frequencies
EP3936898A1 (en) Method and device for group delay variation compensation
CN111753248B (en) Frequency offset time vernier measurement method, system and computer readable storage medium
CN113253310A (en) Method and system for measuring group delay of high-precision GNSS receiver of carrier phase
CN110716215B (en) Method for improving Doppler precision
CN102307050B (en) Digital mirror image inhibition system and method thereof
CN116260690B (en) Carrier phase tracking method for QAM weak signals
Gunst et al. Application of digital wide band mismatch calibration to an I/Q receiver
CN115118304B (en) Frequency hopping signal parameter measurement method for vector signal analyzer
CN115598675B (en) Method and device for extracting ionosphere delay by Beidou three-frequency phase conversion
CN112162153B (en) Sine wave signal frequency estimation method based on phase straight line fitting
CN113640834A (en) Method for improving satellite double-difference pseudo range positioning accuracy
CN115561486A (en) Laser Doppler frequency spectrum correction method based on power spectrum analysis

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