CN113341686B - Single-satellite multi-calendar timing method and device, computer equipment and storage medium - Google Patents
Single-satellite multi-calendar timing method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN113341686B CN113341686B CN202110752475.6A CN202110752475A CN113341686B CN 113341686 B CN113341686 B CN 113341686B CN 202110752475 A CN202110752475 A CN 202110752475A CN 113341686 B CN113341686 B CN 113341686B
- Authority
- CN
- China
- Prior art keywords
- epoch
- local clock
- sliding window
- satellite
- difference
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000006243 chemical reaction Methods 0.000 claims abstract description 13
- 238000012937 correction Methods 0.000 claims description 39
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000013075 data extraction Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a single-star multi-calendar timing method, a single-star multi-calendar timing device, computer equipment and a storage medium, wherein the method comprises the following steps: acquiring a single satellite number, and extracting data of all epoch single satellites from a pre-stored one-way time service data file according to the satellite number; calculating the local clock error of each epoch single star according to the data of all epoch single stars; the method comprises the steps of carrying out unit conversion on local clock differences of all epochs to obtain the starting time and the ending time of a time period with the same frequency difference, obtaining the frequency difference of each epoch in a preset sliding window width through a sliding window according to data of all epoch singles, the converted local clock differences of each epoch and each epoch in the preset sliding window width, and obtaining the local clock differences of each epoch in the time period with the same frequency difference in consideration according to the frequency difference, the starting time and the ending time of each epoch in the preset sliding window width, each epoch in the preset sliding window and the converted local clock differences of each epoch. Improve the local clock error accuracy.
Description
Technical Field
The invention belongs to the technical field of satellite navigation receiving, and particularly relates to a single-satellite multi-epoch time correction method and device, computer equipment and a storage medium.
Background
With the rapid development of modern society and the continuous improvement of the requirement on the time synchronization accuracy, the satellite time service research and application become an important trend for the development of the time frequency field in China. The satellite time service has an important function in the field of basic research and also has wide application in national defense and national economic construction. The single-satellite time service is a simpler time service method in the satellite time service technology, and is a premise of viewing some time service methods such as time service and the like together, so that the problem of precision of one-way time service is researched, the time service precision of the Beidou satellite can be improved, and the follow-up development work of the Beidou satellite time service is facilitated.
The single-satellite time service means that high-precision system time can be obtained only by receiving on the premise that the user equipment does not transmit signals. The purpose of single-satellite time service is to obtain the time difference between a local clock and system time, namely the local clock difference, and then directly or indirectly synchronize the local time to the system time. The Beidou satellite unidirectional time service obtains the time deviation between the local time and the Beidou satellite. Therefore, the accurate level of the local clock difference directly influences the accurate level of the single-satellite time service.
Disclosure of Invention
In view of the above technical problems, the present invention provides a single-satellite multi-epoch timing method, apparatus, computer device, and storage medium, which can improve the accuracy level of local clock error.
The technical scheme adopted by the invention for solving the technical problems is as follows:
in one embodiment, a single-star multi-epoch timing method includes the steps of:
step S100: acquiring a single satellite number, and extracting data of all epoch single satellites from a pre-stored one-way time service data file according to the satellite number;
step S200: calculating to obtain the local clock error of each epoch single satellite without considering the frequency error according to the data of all epoch single satellites;
step S300: the method comprises the steps of carrying out unit conversion on the local clock difference of each epoch unit satellite without considering the frequency difference to obtain the starting time and the ending time of a time period with the same frequency difference, obtaining the frequency difference of each epoch in the preset sliding window width according to the data of all epoch unit satellites, the converted local clock difference of each epoch unit satellite without considering the frequency difference and each epoch in the preset sliding window width through a sliding window type, and completing time correction according to the frequency difference, the starting time and the ending time of each epoch in the preset sliding window width, the frequency difference of each epoch in the preset sliding window and the converted local clock difference of each epoch unit satellite without considering the frequency difference.
Preferably, the data for all epoch single stars includes raw pseudoranges, satellite positions, satellite clock offsets, ionospheric delay corrections, tropospheric delay corrections, autorotation of the earth, observation point positions.
Preferably, step S200 is specifically:
wherein, Δ tRjThe local clock offset without taking the frequency difference into account for each epoch, ρ' is the raw pseudorange, (X, y, z) is the position of the observation point, and X, y, z are the positions of the observation point in the X, y and z axes, respectively, (Xj,Yj,Zj) As satellite position, Xj、Yj、ZjPosition of the satellite in the x, y and z axes, respectively, Δ tjIs the satellite clock error, deltaionFor ionospheric delay correction, deltatroFor tropospheric delay correction, Δ ρ is the earth rotation correction.
Preferably, in step S300, the obtaining, through a sliding window, the frequency difference of each epoch in the preset sliding window width according to the data of all epoch singles, the converted local clock difference of each epoch without considering the frequency difference, and each epoch in the preset sliding window width specifically includes:
X＝(HTH)-1HTY
wherein,
Δ t in H matrixjRepresenting the time interval from the current epoch to the first epoch in the current sliding window, Δ f in the X matrix representing the local frequency difference calculated by the least squares method, Δ t in the Y matrixR' local clock difference, Δ t in Y matrix, when the frequency difference is not considered for the first epoch in the current sliding windowRj' indicating that each epoch takes into account the corresponding satelliteThe local clock error after the satellite clock error, the ionospheric delay correction, the tropospheric delay correction and the earth rotation correction.
Preferably, in step S300, the local clock offset after frequency offset is considered for each epoch in a time period with the same frequency offset according to the frequency offset of each epoch in the preset sliding window width, the start time and the end time, each epoch in the preset sliding window, and the local clock offset without frequency offset being considered for each epoch after conversion, which is specifically:
ΔtRj*＝ΔtR+(Δt×Δf)
wherein, Δ tRjRepresenting the local clock offset, Δ t, of each epoch after conversion without taking into account the frequency offsetRThe local clock difference after the frequency difference is considered for each epoch, and Δ t represents the time difference between the current epoch and the first epoch in the current sliding window.
Preferably, step S300 is followed by:
step S400: and extracting a known clock error from a pre-stored unidirectional time service data file, and obtaining a clock error according to the local clock error and the known clock error after the frequency difference is considered.
Preferably, step S400 is specifically:
D_value＝ΔtR-ΔtR*
wherein D _ value represents a clock error, Δ tRLocal clock difference after frequency difference, Δ t, is taken into account for each epochRDenotes the known clock error for each epoch.
In one embodiment, a single-star multi-calendar timing device comprises:
the data extraction module is used for acquiring a single satellite number and extracting data of all single epochs from a pre-stored one-way time service data file according to the satellite number;
the local clock error calculation module is used for calculating and obtaining the local clock error of each epoch single satellite without considering the frequency error according to the data of all epoch single satellites;
the local clock difference calculation module after considering the frequency difference is used for performing unit conversion on the local clock difference of each epoch monaural without considering the frequency difference to obtain the start time and the end time of the time period with the same frequency difference, the preset sliding window width comprises a preset number of epochs, the frequency difference of each epoch in the preset sliding window width is obtained according to the data of all epoch monaural, the converted local clock difference of each epoch monaural without considering the frequency difference and each epoch in the preset sliding window width through a sliding window, and the local clock difference of each epoch in the time period with the same frequency difference after considering the frequency difference is obtained according to the frequency difference, the start time and the end time of each epoch in the preset sliding window width, each epoch in the preset sliding window and the converted local clock difference of each epoch monaural without considering the frequency difference, so that the time correction is completed.
In an embodiment, a computer device comprises a memory and a processor, the memory storing a computer program, the processor implementing the steps of the above method when executing the computer program.
In an embodiment, a computer-readable storage medium has stored thereon a computer program which, when being executed by a processor, carries out the steps of the above-mentioned method.
According to the single-satellite multi-epoch timing method, the single-satellite multi-epoch timing device, the computer equipment and the storage medium, because the frequency difference is unchanged within a period of time, on the basis of obtaining the local clock difference of each epoch single satellite without considering the frequency difference through calculation, the frequency difference of each epoch is obtained through calculation by using the preset sliding window, and then the influence of the frequency difference among the epochs is eliminated through the sliding window, so that more accurate local clock difference can be obtained, the timing of the local clock difference of the single-satellite multi-epoch is realized, and the accuracy level of the local clock difference is effectively improved.
Drawings
FIG. 1 is a flowchart of a single-star multi-epoch timing method according to an embodiment of the present invention;
FIG. 2 is a flowchart of a single-star multi-epoch timing method according to another embodiment of the present invention;
FIG. 3 is a single-satellite multi-epoch local clock error versus time plot;
FIG. 4 is a single-satellite single-epoch local clock error versus time plot;
FIG. 5 is a comparison graph of local clock error for single-satellite single epoch and multi-epoch;
FIG. 6 is a flow chart of verifying that single-satellite multi-epoch is more accurate than single-epoch.
Detailed Description
In order to make the technical solutions of the present invention better understood, the present invention is further described in detail below with reference to the accompanying drawings.
In one embodiment, as shown in fig. 1, a single-star multi-calendar timing method includes the following steps:
step S100: and acquiring a single satellite number, and extracting the data of all the epoch single satellites from a pre-stored one-way time service data file according to the satellite number.
Specifically, the single satellite time service means that high-precision system time can be obtained only by receiving on the premise that the user equipment does not transmit signals. The purpose of single-satellite time service is to obtain the time difference between a local clock and system time, namely the local clock difference, and then directly or indirectly synchronize the local time to the system time. The Beidou satellite unidirectional time service obtains the time deviation between the local time and the Beidou satellite. Therefore, the accurate level of the local clock difference directly influences the accurate level of the single-satellite time service.
Further, the data for all epoch singles includes raw pseudoranges, satellite positions, satellite clock offsets, ionospheric delay corrections, tropospheric delay corrections, autorotation of the earth corrections, and observation point positions.
Step S200: and calculating the local clock error of each epoch single satellite without considering the frequency error according to the data of all epoch single satellites.
Specifically, calculating the local clock offset using a single satellite and multiple epochs means that the local clock offsets of multiple epochs are calculated sequentially according to the data of a single satellite. Since the coordinates x, y, z of the observation point location are known, it is assumed that at time epoch t, the satellite S is observed by the receiver at the observation point location with known coordinatesjThe corresponding pseudo-range observation equation is:
wherein,ΔtRjlocal clock error for each epoch single star, unit: meter (converted by multiplying by the speed of light), ρ' is the raw pseudorange, (x, y, z) is the position of the observation point, and x, y, z are the positions of the observation point on the x, y, and z axes, in units: rice, (X)j,Yj,Zj) As satellite position, Xj、Yj、ZjThe position of the satellite in the x, y and z axes, respectively, in units: rice, Delta tjIs the satellite clock error, unit: meter (converted by the speed of light), deltaionFor ionospheric delay correction, deltatroFor tropospheric delay correction, Δ ρ is the earth rotation correction.
Because the earth's fixed coordinate system corresponding to the time of satellite signal transmission and the earth's fixed coordinate system corresponding to the time of signal reception by the receiver are different due to the influence of earth rotation, the original pseudorange should be corrected by adding earth rotation.
Therefore, the local clock offset of the receiver at epoch t under the premise of considering the influence of other errors is:
wherein, Δ tRjFor each epoch, the local clock offset without taking frequency offset into account, ρ' is the original pseudo range, X, y and z are the observation point positions, and Xj、Yj、ZjAs the satellite position, Δ tjIs the satellite clock error, deltaionFor ionospheric delay correction, deltatroFor tropospheric delay correction, Δ ρ is the earth rotation correction.
Step S300: the method comprises the steps of carrying out unit conversion on the local clock difference of each epoch unit satellite without considering the frequency difference to obtain the starting time and the ending time of a time period with the same frequency difference, obtaining the frequency difference of each epoch in the preset sliding window width according to the data of all epoch unit satellites, the converted local clock difference of each epoch unit satellite without considering the frequency difference and each epoch in the preset sliding window width through a sliding window type, and completing time correction according to the frequency difference, the starting time and the ending time of each epoch in the preset sliding window width, the frequency difference of each epoch in the preset sliding window and the converted local clock difference of each epoch unit satellite without considering the frequency difference.
Specifically, the unit of the local clock difference needs to be converted from meters to seconds. The frequency difference means that there exists a frequency difference between the frequency received at the receiver and the frequency of the transmitted signal, the magnitude of the frequency difference is related to the moving direction of the object relative to the receiver, and the higher the moving speed, the larger the frequency difference. Therefore, the influence of the frequency difference is eliminated every time.
In one embodiment, in step S300, unit conversion is performed on the local clock offset of each epoch without considering the frequency offset, specifically:
ΔtRj*＝ΔtRj÷c
wherein, Δ tRjLocal clock offset, Δ t, for each epoch without taking into account frequency offsetRjAnd c is the speed of light, and takes 299792458.458 m/s.
In one embodiment, in step S300, the frequency difference of each epoch in the preset sliding window width obtained through the sliding window according to the data of all epoch singles, the converted local clock difference of each epoch without considering the frequency difference, and each epoch in the preset sliding window width is specifically:
X＝(HTH)-1HTY
wherein,
Δ t in H matrixjRepresenting the time interval from the current epoch to the first epoch in the current sliding window, Δ f in the X matrix representing the local frequency difference calculated by the least squares method, Δ t in the Y matrixR' denotes the local clock difference when the frequency difference is not taken into account for the first epoch in the current sliding window, Δ t in the Y matrixRj' means that each epoch takes into account the corresponding satellite clock error, ionospheric delay correction, tropospheric delay correction and local clock error after earth rotation correction.
Specifically, on the basis of single-satellite single-epoch local clock offset calculation, frequency offset influence between epochs is eliminated, and more accurate local clock offset, namely the local clock offset after frequency offset is considered, can be obtained. Since the frequency difference is constant in a period of time, the frequency difference of each epoch in the sliding window is calculated by using the sliding window with the width N, and further, although the frequency difference is generally constant in a period of time, but slightly changes, the frequency difference of each epoch in the current sliding window is recalculated every time the sliding window slides backwards, so that the obtained frequency difference is more accurate.
In one embodiment, in step S300, the local clock offset after frequency offset is considered for each epoch in a time period with the same frequency offset according to the frequency offset of each epoch in the preset sliding window width, the start time and the end time, each epoch in the preset sliding window, and the local clock offset without frequency offset being considered for each epoch after conversion, which is specifically:
ΔtRj*＝ΔtR+(Δt×Δf)
wherein, Δ tRjRepresenting the local clock offset, Δ t, of each epoch after conversion without taking into account the frequency offsetRThe local clock difference after the frequency difference is considered for each epoch, and Δ t represents the time difference between the current epoch and the first epoch in the current sliding window.
Specifically, after the frequency difference is obtained, the local clock difference after the frequency difference is considered by the second epoch in the current sliding window is calculated according to the frequency difference, the sliding window moves backwards, the local clock difference after the frequency difference is considered by the second epoch in the current sliding window is solved again, and similarly, the local clock difference after the frequency difference is considered by each epoch in the time period with the same frequency difference is solved in sequence, and further, in the embodiment, 30 epochs are the width of the sliding window, the frequency difference of the thirty epochs is calculated by using the 30 epochs according to a frequency difference calculation method, the frequency offset is used to eliminate the influence of the frequency offset of the second epoch in the 30 epochs, then the sliding window is moved backwards, the frequency offset of the next 30 epochs is continuously calculated, and the method is used for eliminating the influence caused by the frequency difference of the second epoch in the next 30 epochs until the epochs in the time period with the same frequency difference are calculated, thereby completing the time correction.
In one embodiment, as shown in fig. 2, step S300 is followed by:
step S400: and extracting a known clock error from a pre-stored unidirectional time service data file, and obtaining a clock error according to the local clock error and the known clock error after the frequency difference is considered.
In one embodiment, step S400 specifically includes:
D_value＝ΔtR-ΔtR*
where D _ value represents the clock error, Δ tRLocal clock difference after frequency difference, Δ t, is taken into account for each epochRDenotes the known clock error for each epoch. The corresponding local clock error-time diagram of single-satellite multi-epoch is shown in fig. 3, and the local clock error of single-satellite single-epoch can be obtained according to the known clock error corresponding to each epoch and the local clock error after the frequency difference is not considered for each epoch, as shown in fig. 4.
Specifically, it can be seen that the pseudorange measurement accuracy of the receiver, the position accuracy known by the observation stations, the position accuracy of the satellites, the accuracy of the satellite clock error, the accuracy of the tropospheric delay correction and the ionospheric delay correction, and the like are main factors that influence the accuracy level of the receiver clock error.
In one embodiment, as shown in fig. 5 and 6, the local clock error of the single-satellite single epoch and the local clock error of the single-satellite multi-epoch are calculated by using a single-satellite single-epoch local clock error calculation principle and a single-satellite multi-epoch local clock error calculation principle, and the single-satellite single-epoch local clock error and the single-satellite multi-epoch local clock error are respectively obtained according to the local clock error of the single-satellite single epoch, the single-satellite multi-epoch local clock error and the known clock error in the same time period, and compared with the error, the accuracy of the single-satellite multi-epoch is higher than that of the single-satellite single epoch.
According to the single-satellite multi-epoch timing method, the frequency difference is unchanged within a period of time, so that on the basis of obtaining the local clock difference of each epoch single satellite without considering the frequency difference through calculation, the frequency difference of each epoch in the current sliding window is obtained through calculation by using the preset sliding window, and then the frequency difference influence among the epochs is eliminated through the sliding window, so that more accurate local clock difference can be obtained, the timing of the local clock difference of the single-satellite multi-epoch is realized, and the accuracy level of the local clock difference can be effectively improved.
In one embodiment, the single-star multi-calendar timing device comprises a data extraction module, a local clock difference calculation module and a local clock difference calculation module after frequency difference is considered,
the data extraction module is used for acquiring a single satellite number and extracting data of all single epochs from a pre-stored one-way time service data file according to the satellite number;
the local clock error calculation module of each epoch single satellite is used for calculating and obtaining the local clock error of each epoch single satellite without considering the frequency error according to the data of all epoch single satellites;
the clock difference calculation module after considering the frequency difference is used for carrying out unit conversion on the local clock difference of each epoch unit satellite without considering the frequency difference to obtain the starting time and the ending time of a time period with the same frequency difference, a preset number of epochs are contained in the preset sliding window width, the frequency difference of each epoch in the preset sliding window width is obtained through a sliding window according to the data of all epoch unit satellites, the converted local clock difference of each epoch unit satellite without considering the frequency difference and each epoch in the preset sliding window width, and the local clock difference of each epoch in the time period with the same frequency difference after considering the frequency difference is obtained according to the frequency difference, the starting time and the ending time of each epoch in the preset sliding window width, each epoch in the preset sliding window and the converted local clock difference of each epoch unit satellite without considering the frequency difference, so that the time correction is completed.
For specific limitations of the single-satellite multi-epoch timing device, reference may be made to the above limitations of the single-satellite multi-epoch timing method, which are not described herein again. The modules in the single-satellite multi-calendar timing device can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device includes a memory storing a computer program and a processor implementing the steps of a single-star multi-epoch timing method when the processor executes the computer program.
In one embodiment, a computer readable storage medium has stored thereon a computer program which, when executed by a processor, performs the steps of a single-star multi-epoch timing method.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-readable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein. The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The single-star multi-epoch timing method, the single-star multi-epoch timing device, the computer device and the storage medium provided by the invention are described in detail above. The principles and embodiments of the present invention are explained herein using specific examples, which are presented only to assist in understanding the core concepts of the present invention. It should be noted that, for those skilled in the art, it is possible to make various improvements and modifications to the present invention without departing from the principle of the present invention, and those improvements and modifications also fall within the scope of the claims of the present invention.
Claims (10)
1. A single-satellite multi-calendar timing method is characterized by comprising the following steps:
step S100: acquiring a single satellite number, and extracting data of the single satellite of all epochs from a pre-stored unidirectional time service data file according to the satellite number;
step S200: calculating to obtain the local clock error of the single satellite of each epoch without considering the frequency error according to the data of the single satellite of all epochs;
step S300: the unit conversion is carried out on the local clock offset of the single satellite of each epoch without considering the frequency offset, the starting time and the ending time of the time period with the same frequency offset are obtained, the preset sliding window width contains the preset number of epochs, obtaining the frequency difference of each epoch in the preset sliding window width through a sliding window according to the data of the single satellite of all epochs, the converted local clock difference of the single satellite of each epoch without considering the frequency difference and each epoch in the preset sliding window width, and obtaining the local clock offset of each epoch after frequency offset is considered in the time period with the same frequency offset according to the frequency offset of each epoch in the preset sliding window width, the starting time and the ending time, each epoch in the preset sliding window and the converted local clock offset of each epoch of the single satellite without frequency offset, and finishing time correction.
2. The method of claim 1, wherein all epochs of the single-satellite data comprise raw pseudoranges, satellite positions, satellite clock error, ionospheric delay corrections, tropospheric delay corrections, earth rotation corrections, observation point positions.
3. The method according to claim 2, wherein step S200 is specifically:
wherein, Δ tRjThe local clock offset without taking the frequency difference into account for each epoch, ρ' is the raw pseudorange, (X, y, z) is the position of the observation point, and X, y, z are the positions of the observation point in the X, y and z axes, respectively, (Xj,Yj,Zj) As satellite position, Xj、Yj、ZjPosition of the satellite in the x, y and z axes, Δ tjIs the satellite clock error, deltaionFor ionospheric delay correction, deltatroFor tropospheric delay correction, Δ ρ is the earth rotation correction.
4. The method according to claim 3, wherein the step S300 of obtaining the frequency difference of each epoch in the preset sliding window width through a sliding window according to the data of the single star of all epochs, the converted local clock difference of each epoch without considering the frequency difference, and each epoch in the preset sliding window width specifically comprises:
X＝(HTH)-1HTY
wherein,
Δ t in H matrixjRepresenting the time interval from the current epoch j to the first epoch in the current sliding window, Δ f in the X matrix representing the local frequency difference calculated by the least square method, Δ t in the Y matrixR' local clock difference, Δ t in Y matrix, when the frequency difference is not considered for the first epoch in the current sliding windowRj' means that each epoch takes into account the corresponding satellite clock error, ionospheric delay correction, tropospheric delay correction and local clock error after earth rotation correction.
5. The method according to claim 4, wherein in step S300, the local clock offset after frequency offset is considered for each epoch in a time period with the same frequency offset is obtained according to the frequency offset of each epoch in the preset sliding window width, the start and end times, each epoch in the preset sliding window, and the local clock offset without frequency offset being considered for each scaled epoch, and specifically:
ΔtRj*＝ΔtR+(Δt×Δf)
wherein, Δ tRjRepresenting the local clock offset, Δ t, of each epoch after conversion without taking into account the frequency offsetRThe local clock difference after the frequency difference is considered for each epoch, and Δ t represents the time difference between the current epoch and the first epoch in the current sliding window.
6. The method of claim 1, wherein step S300 is further followed by:
step S400: and extracting a known clock difference from the pre-stored unidirectional time service data file, and obtaining a clock difference error according to the local clock difference after the frequency difference is considered and the known clock difference.
7. The method according to claim 6, wherein step S400 is specifically:
D_value＝ΔtR-ΔtR*
where D _ value represents the clock error, Δ tRLocal clock difference after frequency difference, Δ t, is taken into account for each epochRDenotes the known clock error for each epoch.
8. A single-satellite multi-epoch timing device, comprising:
the data extraction module is used for acquiring a single satellite number and extracting the data of the single satellite of all epochs from a pre-stored one-way time service data file according to the satellite number;
the local clock error calculation module is used for calculating and obtaining the local clock error of each epoch, which does not consider the frequency error, of the single satellite according to the data of the single satellite of all epochs;
a local clock difference calculation module with frequency difference considered, configured to perform unit conversion on the local clock difference of the single satellite of each epoch without frequency difference considered, obtain start and end times of a time period with the same frequency difference, where a preset sliding window width includes a preset number of epochs, obtain a frequency difference of each epoch within the preset sliding window width through a sliding window according to data of the single satellite of all epochs, the converted local clock difference of the single satellite without frequency difference considered of each epoch and each epoch within the preset sliding window width, obtain the local clock difference of each epoch within the time period with the same frequency difference considered of each epoch within the preset sliding window width according to the frequency difference of each epoch within the preset sliding window width, the start and end times, each epoch within the preset sliding window, and the converted local clock difference of the single satellite without frequency difference considered of each epoch, and finishing time correction.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the method of any of claims 1 to 7.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110752475.6A CN113341686B (en) | 2021-07-02 | 2021-07-02 | Single-satellite multi-calendar timing method and device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110752475.6A CN113341686B (en) | 2021-07-02 | 2021-07-02 | Single-satellite multi-calendar timing method and device, computer equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113341686A CN113341686A (en) | 2021-09-03 |
| CN113341686B true CN113341686B (en) | 2022-05-13 |
Family
ID=77482401
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110752475.6A Active CN113341686B (en) | 2021-07-02 | 2021-07-02 | Single-satellite multi-calendar timing method and device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113341686B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116125512B (en) * | 2023-04-13 | 2023-07-18 | 中国科学院精密测量科学与技术创新研究院 | PPP self-adaptive clock difference model estimation method considering clock frequency time-varying characteristics |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08146166A (en) * | 1994-11-24 | 1996-06-07 | Advantest Corp | Reference frequency generating device |
| CN101035038A (en) * | 2007-04-24 | 2007-09-12 | 北京邮电大学 | Self-adapted method for detecting the clock resetting based on the mode recognition |
| JP2008051529A (en) * | 2006-08-22 | 2008-03-06 | Seiko Epson Corp | Time correcting device, electronic watch with time correcting device, and method for correcting time |
| CN104181550A (en) * | 2014-08-19 | 2014-12-03 | 北京无线电计量测试研究所 | Common view time and frequency transmitting method based on BeiDou navigation satellite system |
| CN104678408A (en) * | 2015-02-05 | 2015-06-03 | 清华大学 | Satellite borne navigation receiver time service method, time service type satellite borne navigation receiver and satellite borne navigation application system |
| CN106019333A (en) * | 2016-05-16 | 2016-10-12 | 中国人民解放军国防科学技术大学 | Beidou navigation signal vector tracking method based on incoherent discriminator |
| KR20170062208A (en) * | 2015-11-27 | 2017-06-07 | 한전케이디엔주식회사 | Time synchronized system |
| CN108931915A (en) * | 2018-05-08 | 2018-12-04 | 和芯星通科技(北京)有限公司 | Utilize time service method and device, the computer readable storage medium of navigation satellite |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19914355A1 (en) * | 1999-03-30 | 2000-10-05 | Wolfgang Schaefer | Method for synchronizing remote clocks with central clock via satellite |
| JP4600480B2 (en) * | 2008-01-16 | 2010-12-15 | セイコーエプソン株式会社 | Electronic clock |
| CN101866009B (en) * | 2010-06-10 | 2012-07-25 | 北京航空航天大学 | Integrity monitoring method suitable for satellite lacking situation of global navigation satellite system (GNSS) |
| CN103149571B (en) * | 2013-02-18 | 2014-12-10 | 桂林电子科技大学 | GNSS (Global Navigation Satellite System)-based signal aided time frequency difference comprehensive correction method |
| CN103226324B (en) * | 2013-03-06 | 2016-12-28 | 中国计量科学研究院 | A kind of high-accuracy temporal frequency source taming time and frequency standard in real time |
| CN110850450A (en) * | 2019-12-03 | 2020-02-28 | 航天恒星科技有限公司 | Adaptive estimation method for satellite clock error parameters |
-
2021
- 2021-07-02 CN CN202110752475.6A patent/CN113341686B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08146166A (en) * | 1994-11-24 | 1996-06-07 | Advantest Corp | Reference frequency generating device |
| JP2008051529A (en) * | 2006-08-22 | 2008-03-06 | Seiko Epson Corp | Time correcting device, electronic watch with time correcting device, and method for correcting time |
| CN101035038A (en) * | 2007-04-24 | 2007-09-12 | 北京邮电大学 | Self-adapted method for detecting the clock resetting based on the mode recognition |
| CN104181550A (en) * | 2014-08-19 | 2014-12-03 | 北京无线电计量测试研究所 | Common view time and frequency transmitting method based on BeiDou navigation satellite system |
| CN104678408A (en) * | 2015-02-05 | 2015-06-03 | 清华大学 | Satellite borne navigation receiver time service method, time service type satellite borne navigation receiver and satellite borne navigation application system |
| KR20170062208A (en) * | 2015-11-27 | 2017-06-07 | 한전케이디엔주식회사 | Time synchronized system |
| CN106019333A (en) * | 2016-05-16 | 2016-10-12 | 中国人民解放军国防科学技术大学 | Beidou navigation signal vector tracking method based on incoherent discriminator |
| CN108931915A (en) * | 2018-05-08 | 2018-12-04 | 和芯星通科技(北京)有限公司 | Utilize time service method and device, the computer readable storage medium of navigation satellite |
Non-Patent Citations (2)
| Title |
|---|
| GPS接收机模糊搜索粗时定位方法;靖守让;《国防科技大学学报》;20150630;第37卷(第3期);34-38 * |
| 单星频差无源定位技术原理及仿真;蒋东旭;《舰船电子对抗》;20170430;第40卷(第2期);16-22 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113341686A (en) | 2021-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109581452B (en) | GNSS reference station carrier phase integer ambiguity resolution method | |
| WO2022156481A1 (en) | Ephemeris forecasting method and apparatus | |
| RU2749667C1 (en) | Method and system for fast and accurate positioning | |
| CN111045034B (en) | GNSS multi-system real-time precise time transfer method and system based on broadcast ephemeris | |
| EP4258024A1 (en) | Clock error predicting method and device | |
| CN106168672B (en) | A kind of GNSS multimode single-frequency RTK Cycle Slips Detection and device | |
| AU2007300586B2 (en) | Method for using three GPS frequencies to resolve whole-cycle carrier-phase ambiguities | |
| RU2432584C2 (en) | Method of determining coordinates of satellite radio navigation system (srns) mobile receiver | |
| CN107656300B (en) | Satellite/inertia ultra-tight combination method based on Beidou/GPS dual-mode software receiver | |
| CN1930488A (en) | Method for backup dual-frequency navigation during brief periods when measurement data is unavailable on one of two frequencies | |
| CN110146906A (en) | Remote time transmission method based on single poor carrier phase observation data | |
| CN108845339B (en) | A kind of GNSS localization method and GNSS positioning device | |
| CN105607089A (en) | Pseudo-range-assisted-carrier-phase-based beidou common-view time transmission method | |
| CN110389364A (en) | Localization method, device, computer equipment and storage medium | |
| CN107607971B (en) | Time frequency transmission method based on GNSS common-view time comparison algorithm and receiver | |
| CN104898145B (en) | A kind of fuzziness fixing means and system based on half cycle fuzziness | |
| CN102197317A (en) | Satellite navigation device | |
| CN114994727B (en) | Equipment for realizing high-precision time calibration and satellite positioning | |
| CN112799107B (en) | Dynamic time service system and method | |
| CN113341686B (en) | Single-satellite multi-calendar timing method and device, computer equipment and storage medium | |
| CN114779301B (en) | Satellite navigation real-time precise single-point positioning method based on broadcast ephemeris | |
| CN113325446A (en) | Multi-mode common-frequency GNSS carrier phase time transfer method and system | |
| CN114935770A (en) | Method and device for accelerating precision single-point positioning convergence speed by multiple calendars | |
| CN117388881A (en) | Method and system for tracing satellite-borne atomic clock of low-orbit satellite to UTC (k) | |
| JP2009222438A (en) | Positioning device for movable body |
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 | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210903 Assignee: Hunan Xinglian Aerospace Technology Co.,Ltd. Assignor: CHANGSHA University Contract record no.: X2024980011685 Denomination of invention: A single star multi epoch timing method, device, computer equipment, and storage medium Granted publication date: 20220513 License type: Common License Record date: 20240814 Application publication date: 20210903 Assignee: Hunan Lingyan Huixing Technology Co.,Ltd. Assignor: CHANGSHA University Contract record no.: X2024980011680 Denomination of invention: A single star multi epoch timing method, device, computer equipment, and storage medium Granted publication date: 20220513 License type: Common License Record date: 20240814 |