CN108873024B - Method for assisting loop by using clock drift - Google Patents

Method for assisting loop by using clock drift Download PDF

Info

Publication number
CN108873024B
CN108873024B CN201710666685.7A CN201710666685A CN108873024B CN 108873024 B CN108873024 B CN 108873024B CN 201710666685 A CN201710666685 A CN 201710666685A CN 108873024 B CN108873024 B CN 108873024B
Authority
CN
China
Prior art keywords
clock drift
model
loop
auxiliary
value
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
CN201710666685.7A
Other languages
Chinese (zh)
Other versions
CN108873024A (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.)
Shanghai Huace Navigation Technology Ltd
Original Assignee
Shanghai Huace Navigation Technology Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Huace Navigation Technology Ltd filed Critical Shanghai Huace Navigation Technology Ltd
Publication of CN108873024A publication Critical patent/CN108873024A/en
Application granted granted Critical
Publication of CN108873024B publication Critical patent/CN108873024B/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/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/29Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
    • 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/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30Acquisition or tracking or demodulation of signals transmitted by the system code related

Abstract

The invention relates to a method for utilizing a clock drift auxiliary loop, which comprises the following steps: s1, repairing and removing the generated pseudo range rate, and distinguishing the carrier waves with cycle slip; s2, establishing a velocity measurement equation according to the pseudo-range change rate, fitting parameters of the self-stabilized motion model by the calculated clock drift value, establishing a clock drift model by the measured clock drift, obtaining each parameter of the model through least square, comparing the latest clock drift value with the clock drift value predicted by the model, and stopping output if the difference value is overlarge; and step S3, obtaining an original clock drift, converting different systems into auxiliary parameters of different frequency points, and associating the auxiliary parameters with a carrier ring and a code ring to eliminate Doppler influence caused by crystal oscillator drift. The method of the invention has simple and reliable algorithm and can not cause huge burden on resources. The method can adapt to various levels of crystal oscillators, and the sensitivity and the reliability of data are improved.

Description

Method for assisting loop by using clock drift
Technical Field
The invention relates to the field of navigation, in particular to a method for assisting a loop by using a clock drift.
Background
High-precision and high-sensitivity positioning is widely applied to the fields of various surveying and mapping and carrier navigation. The stability of the data directly affects the performance index of the product, especially the index of the data utilization rate. At present, products of various manufacturers have basic functions, but the method becomes a hot spot of next research on how to improve data quality and product indexes. The manufacturers at home and abroad have gone to the front in this respect, and because of technical blockade and commercial value and the like, the manufacturers at home and abroad do not disclose the technology for improving the data reliability, stability and the like, and the thesis and the patent almost have no description in this respect. The method is based on a self-developed B380 multi-mode multi-frequency high-precision board card, and explains how to utilize clock drift to assist a carrier loop.
Disclosure of Invention
The present invention provides a method for improving loop dynamics and sensitivity. The method aims to eliminate Doppler influence caused by clock drift in the loop tracking process, thereby improving the reliability of the loop. The working idea and the processing flow of the invention are as follows:
s1, repairing and removing the generated pseudo range rate, and distinguishing the carrier waves with cycle slip;
s2, establishing a velocity measurement equation according to the pseudo-range change rate, fitting parameters of the self-stabilized motion model by the calculated clock drift value, establishing a clock drift model by the measured clock drift, obtaining each parameter of the model through least square, comparing the latest clock drift value with the clock drift value predicted by the model, and stopping output if the difference value is overlarge;
and step S3, obtaining an original clock drift, converting different systems into auxiliary parameters of different frequency points, and associating the auxiliary parameters with a carrier ring and a code ring to eliminate Doppler influence caused by crystal oscillator drift.
In the method using the clock drift auxiliary loop, the clock drift value obtained according to the clock drift model is firstly converted into a Doppler frequency word required by channel tracking and then is superposed into an NCO frequency word of the carrier loop.
In the above method using the clock drift auxiliary loop, the code loop is assisted by the doppler value of the carrier loop, and the doppler generated by the clock drift is not directly assisted.
The method for using a clock drift auxiliary loop, wherein the step of generating the auxiliary value corresponding to the loop by the clock drift comprises: the clock drift value is divided by the wavelength of the corresponding loop times the integration time to obtain the doppler value due to the clock drift during one setting of the NCO.
The method using the clock drift auxiliary loop, wherein the step S2 further includes:
and converting the processed clock difference value into auxiliary Doppler required by the loop, dividing the clock difference by the wavelength to obtain Doppler, determining an auxiliary word under the time length according to the setting time of the NCO, and superposing the auxiliary word into the carrier loop, wherein the auxiliary of the code loop is provided by the carrier loop.
The method of the invention has simple and reliable algorithm and can not cause huge burden on resources. The method can adapt to various levels of crystal oscillators, and the sensitivity and the reliability of data are improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the embodiments or the prior art descriptions will be briefly described 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 inventive exercise.
FIG. 1 is a flow chart of a process for utilizing a clock drift assist loop according to an embodiment of the present invention;
FIG. 2 is a flow chart of the steps for setting the NCO in accordance with the present invention.
Detailed Description
In the following description, numerous specific details are set forth in order to provide a more thorough understanding of the present invention. It will be apparent, however, to one skilled in the art, that the present invention may be practiced without one or more of these specific details. In other instances, well-known features have not been described in order to avoid obscuring the invention.
In the following description, for the purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the present invention. The following detailed description of the preferred embodiments of the invention, however, the invention is capable of other embodiments in addition to those detailed.
At present, no relevant patent is provided for describing the influence and processing of the clock drift on the loop in detail, and some loop parameters are simply improved from parameters such as integration time, loop bandwidth and the like. So that the data quality can be improved upwards.
The invention provides a method for utilizing a clock drift auxiliary loop, which is shown in a figure 1 and a figure 2 and comprises the following steps:
s1, repairing and removing the generated pseudo range rate, and distinguishing the carrier waves with cycle slip;
s2, establishing a velocity measurement equation according to the pseudo-range change rate, fitting parameters of an auto-stable motion model (ARMA) by the calculated clock drift value, establishing a clock drift model by the measured clock drift, obtaining each parameter of the model through least square, comparing the latest clock drift value with the clock drift value predicted by the model, and stopping output if the difference value is too large;
and S4, obtaining an original clock drift, converting different systems into auxiliary parameters of different frequency points, and associating the auxiliary parameters with a carrier ring and a code ring to eliminate Doppler influence caused by crystal oscillator drift.
In an alternative embodiment of the present invention, the clock drift value obtained according to the clock drift model is first converted into a doppler frequency word required for channel tracking, and then added to the NCO frequency word of the carrier loop.
The observed quantity required before the velocity measurement equation is established is the pseudo-range change rate, the pseudo-range change rate obtained by the carrier phase is used, the velocity measurement precision and the clock drift precision of the change rate are higher, but the cycle slip is easy to occur in the carrier phase, and the condition needs to be identified and repaired.
Further, the clock difference processed in step S2 needs to be converted into an auxiliary doppler needed by the loop, and the calculation idea is to divide the clock difference by the wavelength to obtain the doppler, then determine the auxiliary word in the time length according to the setting time of the NCO, and then superimpose the auxiliary word on the carrier loop, where the code loop is provided by the carrier loop.
In an alternative embodiment of the invention, the code loop is assisted by the doppler of the carrier loop, and the doppler generated by the clock drift is not directly assisted.
In an alternative embodiment of the present invention, the step of generating the auxiliary value corresponding to the loop by the clock drift comprises: the clock drift value is divided by the wavelength of the corresponding loop times the integration time to obtain the doppler value due to the clock drift during one setting of the NCO.
The method of the invention has simple and reliable algorithm and can not cause huge burden on resources. The method can adapt to various levels of crystal oscillators, and the sensitivity and the reliability of data are improved.
The above description is of the preferred embodiment of the invention. It is to be understood that the invention is not limited to the particular embodiments described above, in that devices and structures not described in detail are understood to be implemented in a manner common in the art; those skilled in the art can make many possible variations and modifications to the disclosed embodiments, or modify equivalent embodiments to equivalent variations, without departing from the spirit of the invention, using the methods and techniques disclosed above. Therefore, any simple modification, equivalent change and modification made to the above embodiments according to the technical essence of the present invention are still within the scope of the protection of the technical solution of the present invention, unless the contents of the technical solution of the present invention are departed.

Claims (1)

1. A method of using a clock drift assist loop, comprising the steps of:
s1, repairing and removing the generated pseudo range rate, and distinguishing the carrier waves with cycle slip;
step S2, a velocity measurement equation is established according to the pseudo-range change rate, the calculated clock drift value is used for fitting parameters from a stable motion model, the measured clock drift is used for establishing a clock drift model, each parameter of the model is obtained through least square, then the latest clock drift value is compared with the clock drift value predicted by the model, if the difference value is too large, the output is cut off, then,
dividing a clock by a wavelength to obtain Doppler, then determining an auxiliary word under the time length according to the setting time of NCO, and then superimposing the auxiliary word into a carrier loop, wherein the auxiliary of a code loop is provided by the carrier loop;
s3, obtaining an original clock drift, converting different systems into auxiliary parameters of different frequency points, and associating the auxiliary parameters with a carrier ring and a code ring to eliminate Doppler influence caused by crystal oscillator drift;
firstly, converting a clock drift value obtained according to a self-stabilizing motion clock drift model into a Doppler frequency word required by channel tracking, and then superposing the Doppler frequency word to an NCO frequency word of a carrier loop;
the assistance of the code ring depends on the Doppler value of the carrier ring, and the Doppler generated by the clock drift is not directly assisted.
CN201710666685.7A 2017-05-15 2017-08-07 Method for assisting loop by using clock drift Active CN108873024B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2017103391265 2017-05-15
CN201710339126 2017-05-15

Publications (2)

Publication Number Publication Date
CN108873024A CN108873024A (en) 2018-11-23
CN108873024B true CN108873024B (en) 2022-01-28

Family

ID=64325440

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710666685.7A Active CN108873024B (en) 2017-05-15 2017-08-07 Method for assisting loop by using clock drift

Country Status (1)

Country Link
CN (1) CN108873024B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995269A (en) * 2014-05-26 2014-08-20 南京理工大学 Method for inertial information assisting GNSS tracking loop
CN106199652A (en) * 2016-06-24 2016-12-07 南京理工大学 A kind of self adaptation vector tracking method of GPS
CN106199655A (en) * 2016-06-24 2016-12-07 南京理工大学 A kind of vector tracking method based on federated filter
CN106501832A (en) * 2016-12-16 2017-03-15 南京理工大学 A kind of fault-tolerant vector tracking GNSS/SINS deep integrated navigation methods

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070211791A1 (en) * 2005-09-07 2007-09-13 Suman Ganguly Coherent processing using composite codes

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995269A (en) * 2014-05-26 2014-08-20 南京理工大学 Method for inertial information assisting GNSS tracking loop
CN106199652A (en) * 2016-06-24 2016-12-07 南京理工大学 A kind of self adaptation vector tracking method of GPS
CN106199655A (en) * 2016-06-24 2016-12-07 南京理工大学 A kind of vector tracking method based on federated filter
CN106501832A (en) * 2016-12-16 2017-03-15 南京理工大学 A kind of fault-tolerant vector tracking GNSS/SINS deep integrated navigation methods

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GNSS_INS标量深组合跟踪技术研究与原型系统验证;张提升;《中国优秀博士学位论文全文数据库 信息科技辑》;20160115(第1期);I136-83 *
Improvement in tracking loop threshold of high dynamic GNSS receiver by installation of crystal oscillator on gyroscopic mounting;Maryam Abedi等;《2015 Joint Conference of the IEEE International Frequency Control Symposium & the European Frequency and Time Forum》;20150702;1-12 *
双系统高动态导航接收机基带技术研究与实现;陈小鱼;《中国优秀硕士学位论文全文数据库 基础科学辑》;20121015(第10期);A008-76 *

Also Published As

Publication number Publication date
CN108873024A (en) 2018-11-23

Similar Documents

Publication Publication Date Title
CN104410413B (en) Atomic frequency standard frequency Correction Method, device and atomic frequency standard
WO2016050055A1 (en) Signal capture method and device, and computer storage medium
CN104485948B (en) The control method and time standard device of a kind of time standard equipment
CN104716904B (en) A kind of method of crystal oscillator frequency compensation
CN103576167B (en) Based on the cycle-slip detection and repair method of HHT and support vector machine
CN104485954B (en) The control method and time device of a kind of time device
CN103344971A (en) Optimization method suitable for GNSS real-time data processing
CN103018555A (en) High-precision electric power parameter software synchronous sampling method
CN109085616A (en) A kind of satellite timing method, device and storage medium
CN104702214A (en) Crystal oscillator frequency compensation method
CN104300969A (en) High-precision synchronization clock realization method based on full digital phase-locked loop
CN203812005U (en) Gps time service system
CN104199061A (en) Method for establishing carrier phase frequency standard of GPS (global position system) and BDS (BeiDou Navigation Satellite system)
CN107566105B (en) Time synchronization equipment compensation method, device, storage medium and computer equipment thereof
CN104062645A (en) Method for measuring phase difference between parametric array difference-frequency wave and same-frequency small-amplitude wave signal
CN206684760U (en) Carrier signal generative circuit for radio-frequency identification transponder equipment
CN108873024B (en) Method for assisting loop by using clock drift
CN107976700A (en) A kind of stabilization output method of satellite navigation receiver pulse per second (PPS)
CN104820225B (en) The anti-interference monitoring system of receiver and method with carrier phase is monitored based on markers
CN106371092A (en) Deformation monitoring method based on GPS and strong-motion seismograph observation adaptive combination
CN112415542A (en) Deformation monitoring and resolving method based on combination of Beidou and INS
CN103884910A (en) Electric system phasor calculating method suitable for frequency deviation
CN203366002U (en) Atomic clock
JP2007278708A (en) Satellite navigation system
CN104990616A (en) Multi-path desynchronous pulse counting compensation method based on cascading adaptive trap wave device

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