CN115685271A - Two-stage rapid signal capture method of time division navigation signal under large Doppler - Google Patents

Two-stage rapid signal capture method of time division navigation signal under large Doppler Download PDF

Info

Publication number
CN115685271A
CN115685271A CN202211702895.4A CN202211702895A CN115685271A CN 115685271 A CN115685271 A CN 115685271A CN 202211702895 A CN202211702895 A CN 202211702895A CN 115685271 A CN115685271 A CN 115685271A
Authority
CN
China
Prior art keywords
code phase
signal
pseudo
result
doppler
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211702895.4A
Other languages
Chinese (zh)
Other versions
CN115685271B (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 University of Defense Technology
Original Assignee
National University of Defense Technology
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 University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202211702895.4A priority Critical patent/CN115685271B/en
Publication of CN115685271A publication Critical patent/CN115685271A/en
Application granted granted Critical
Publication of CN115685271B publication Critical patent/CN115685271B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The application relates to a two-stage rapid signal capturing method of time division navigation signals under large Doppler. The method comprises the following steps: acquiring a digital baseband signal received by a navigation receiver, and sampling the digital baseband signal to obtain a first sampling signal; traversing a first code phase searching range, and performing code phase capturing on a time division pulse sequence corresponding to the first sampling signal to obtain a first pseudo code phase capturing result; after waiting for a preset time interval, sampling the digital baseband signal to obtain a second sampling signal, and obtaining a second code phase search range according to the pre-estimated maximum Doppler of the digital baseband signal; traversing a second code phase searching range, and performing code phase capturing on a time division pulse sequence corresponding to a second sampling signal to obtain a second pseudo code phase capturing result; and obtaining carrier Doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result. By adopting the method, the time division spread spectrum signal with large Doppler frequency offset can be rapidly captured.

Description

Two-stage rapid signal capture method of time division navigation signal under large Doppler
Technical Field
The application relates to the technical field of satellite navigation, in particular to a two-stage rapid signal acquisition method of time division navigation signals under large Doppler.
Background
The Beidou navigation system has the characteristics of global coverage, all-day work and the like, and the Beidou navigation system is widely applied to various industries, including the fields of vehicle and aircraft navigation, personal navigation, bridge monitoring, modern agriculture, precise mapping and the like. In a GPS receiver, acquisition of a signal is a prerequisite for signal tracking and demodulation of navigation data bits. During acquisition, carrier Doppler frequency shift search and C/A code initial phase search need to be completed for each satellite in a GPS constellation at the same time. The Beidou system improves the measurement performance of the system by adding the measurement link of the Ka frequency point. The Ka frequency point signal has high carrier frequency and large Doppler dynamic range, and adopts a time division system signal, namely, a time division pulse control sequence is modulated on the basis of the traditional continuous navigation signal, and the continuous signal is changed into a discontinuous signal. On the other hand, the low-orbit satellite is adopted to broadcast the navigation enhancement signal, so that the performance of the Beidou system can be further improved, but the low-orbit satellite has large dynamic state, so that the Doppler variation range of the enhancement signal is large, and meanwhile, in order to solve the problem of receiving and transmitting isolation of the low-orbit satellite, the enhancement signal can also adopt a time division system, so that the low-orbit navigation enhancement signal received on the ground is also a time division signal with large Doppler.
In the conventional method, to complete signal acquisition, a search needs to be performed in an uncertain region of a pseudo code and a doppler frequency, and for a time division navigation signal with large doppler, the conventional continuous signal acquisition method suffers from large performance loss and needs to overcome the problem of a sharp increase of a search space caused by large doppler. The method using the assistance information can compress the search range during signal acquisition, but needs to know the ephemeris of the satellite and the approximate position of the satellite in advance.
Disclosure of Invention
In view of the foregoing, it is necessary to provide a two-stage fast signal acquisition method for time division navigation signals under large doppler.
A two-stage fast signal acquisition method for time-division navigation signals under large doppler, the method comprising:
acquiring a digital baseband signal received by a navigation receiver, and sampling the digital baseband signal to obtain a first sampling signal; the digital baseband signal comprises a time division spread spectrum signal with large Doppler frequency offset;
traversing a preset first code phase searching range, performing code phase capturing on the time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase capturing result in the first code phase searching range, and obtaining a first pseudo code phase capturing result according to the maximum value of the pseudo code phase capturing result;
after waiting for a preset time interval at an initial sampling point of the first sampling signal, sampling the digital baseband signal to obtain a second sampling signal, and calculating according to a pre-estimated maximum Doppler corresponding to the digital baseband signal to obtain a second code phase search range;
traversing the second code phase search range, and performing code phase acquisition on the time division pulse sequence corresponding to the second sampling signal to obtain a corresponding second pseudo code phase acquisition result;
and obtaining carrier Doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result.
In one embodiment, the method further comprises the following steps: the mathematical model of the digital baseband signal is as follows:
Figure 674252DEST_PATH_IMAGE002
wherein,
Figure 122551DEST_PATH_IMAGE003
to navigate the digital baseband signals received by the receiving device,
Figure 558080DEST_PATH_IMAGE004
amplitude of carrier wave,
Figure 484448DEST_PATH_IMAGE005
In order to provide a satellite navigation message,
Figure 705345DEST_PATH_IMAGE006
in order to sample the points of interest,
Figure 58966DEST_PATH_IMAGE007
in order to delay the transmission of the signal,
Figure 477397DEST_PATH_IMAGE008
is the carrier-wave doppler of the signal,w(k) In the case of baseband noise, the noise is,
Figure 207455DEST_PATH_IMAGE009
for the purpose of the spreading codes in the navigation signal,
Figure 282859DEST_PATH_IMAGE010
is a time-division pulse signal, and has a pulse width,
Figure 807381DEST_PATH_IMAGE011
is a symbol of an imaginary unit of a number,
Figure 951923DEST_PATH_IMAGE012
is the radio frequency carrier initial phase.
In one embodiment, the method further comprises the following steps: performing incoherent accumulation on a plurality of coherent integration results obtained after coherent integration is performed on each time division pulse of the first sampling signal and a local spread spectrum code corresponding to the current code phase to obtain an incoherent accumulation result corresponding to the current code phase; the current code phase is in a preset first code phase searching range; and traversing the first code phase searching range to obtain a non-coherent accumulation result corresponding to each code phase in the first code phase searching range, and obtaining a pseudo code phase capturing result in the first code phase searching range according to the magnitude relation between each non-coherent accumulation result and a decision threshold.
In one embodiment, the method further comprises the following steps: when the incoherent accumulation result is larger than a decision threshold, successfully acquiring the code phase, and obtaining an acquisition result according to the code phase corresponding to the current local spread spectrum code; and when the incoherent accumulation result is less than or equal to a decision threshold, failing to acquire the code phase, sliding a preset pseudo code phase search interval backwards in the first code phase search range, calculating an incoherent accumulation result corresponding to the slid code phase, iterating the process until each code phase in the first code phase search range is searched, stopping iteration, and outputting an acquisition result in the first code phase search range.
In one embodiment, the method further comprises the following steps: performing non-coherent accumulation on a plurality of coherent integration results obtained by performing coherent integration on each time division pulse of the first sampling signal and the local spread spectrum code corresponding to the current code phase, and obtaining a non-coherent accumulation result corresponding to the current code phase as follows:
Figure 220093DEST_PATH_IMAGE013
wherein,
Figure 150003DEST_PATH_IMAGE014
is a code phase of
Figure 111006DEST_PATH_IMAGE015
The result of the corresponding non-coherent accumulation,
Figure 775468DEST_PATH_IMAGE016
is the phase of the baseband digital signal and,
Figure 581750DEST_PATH_IMAGE017
is a base-band complex signal and is,
Figure 631746DEST_PATH_IMAGE018
is as follows
Figure 498070DEST_PATH_IMAGE019
The position of the starting sampling point of each pulse signal,
Figure 961413DEST_PATH_IMAGE020
for local replication of phase is
Figure 696019DEST_PATH_IMAGE021
The spreading code of (a) is used,
Figure 131680DEST_PATH_IMAGE022
in order to be a short-time correlation length,
Figure 434485DEST_PATH_IMAGE023
the number of pulses corresponding to the first sampling signal.
In one embodiment, the method further comprises the following steps: and calculating according to the pre-estimated maximum Doppler corresponding to the digital baseband signal to obtain a second code phase search range as follows:
Figure 542381DEST_PATH_IMAGE024
wherein,
Figure 956045DEST_PATH_IMAGE025
for the second code phase search range,
Figure 980632DEST_PATH_IMAGE026
for the maximum doppler corresponding to the digital baseband signal,
Figure 454339DEST_PATH_IMAGE027
Figure 626695DEST_PATH_IMAGE028
is the radio frequency of the signal and,
Figure 968683DEST_PATH_IMAGE029
in order to be at the pseudo-code rate,
Figure 972411DEST_PATH_IMAGE030
which represents a rounding-up operation, is performed,
Figure 492385DEST_PATH_IMAGE031
the time interval is a time interval of,
Figure 152037DEST_PATH_IMAGE032
the interval is searched for the pseudo code phase.
In one embodiment, the method further comprises the following steps: according to the first pseudo code phase capturing result and the second pseudo code phase capturing result, carrier Doppler is obtained as follows:
Figure 907503DEST_PATH_IMAGE033
wherein,
Figure 126257DEST_PATH_IMAGE034
in the form of a carrier wave doppler, the doppler,
Figure 941766DEST_PATH_IMAGE035
is the ratio of the signal radio frequency to the pseudo code rate,
Figure 760818DEST_PATH_IMAGE036
for the second pseudo-code phase acquisition result,
Figure 319975DEST_PATH_IMAGE037
for the first pseudo-code phase acquisition result,
Figure 422929DEST_PATH_IMAGE038
are time intervals.
In one embodiment, the method further comprises the following steps: the second pseudo code phase acquisition result is:
Figure 143761DEST_PATH_IMAGE039
wherein,
Figure 778004DEST_PATH_IMAGE040
for the second pseudo-code phase acquisition result,
Figure 750639DEST_PATH_IMAGE041
is as follows
Figure 849045DEST_PATH_IMAGE042
Get maximum time corresponding
Figure 361017DEST_PATH_IMAGE015
The value of (a) is,
Figure 279295DEST_PATH_IMAGE043
is a code phase of
Figure 383517DEST_PATH_IMAGE044
The result of the corresponding non-coherent accumulation,
Figure 680637DEST_PATH_IMAGE045
for the first pseudo-code phase acquisition result,
Figure 274430DEST_PATH_IMAGE046
for the second code phase search range,
Figure 539058DEST_PATH_IMAGE047
the interval is searched for a pseudo code phase.
In one embodiment, the method further comprises the following steps: the first pseudo code phase acquisition result is:
Figure 446971DEST_PATH_IMAGE048
wherein,
Figure 988811DEST_PATH_IMAGE049
the result is captured for the first pseudo-code phase,
Figure 363291DEST_PATH_IMAGE050
is the phase of the code and is,
Figure 256161DEST_PATH_IMAGE051
is as follows
Figure 125022DEST_PATH_IMAGE052
Taking maximum time to correspond to
Figure 521368DEST_PATH_IMAGE053
The value of (a) is,
Figure 394647DEST_PATH_IMAGE054
is a code phase of
Figure 384599DEST_PATH_IMAGE055
The result of the corresponding non-coherent accumulation,
Figure 165473DEST_PATH_IMAGE056
the interval is searched for the pseudo-code phase,
Figure 72118DEST_PATH_IMAGE057
the number of cells is searched for the pseudo code phase.
In one embodiment, the method further comprises the following steps: and obtaining the carrier Doppler estimation precision according to the time interval and the pseudo code phase search interval as follows:
Figure 913036DEST_PATH_IMAGE058
wherein,
Figure 718180DEST_PATH_IMAGE059
is the accuracy of the estimation of the carrier doppler,
Figure 178112DEST_PATH_IMAGE060
the interval is searched for the pseudo-code phase,
Figure 814629DEST_PATH_IMAGE061
the ratio of the signal radio frequency to the pseudo code rate; adjusting the time interval and the pseudo code phase search interval to adjust an acquisition accuracy of carrier doppler.
A two-stage fast signal acquisition apparatus for time-division navigation signals under large doppler, the apparatus comprising:
the signal acquisition module is used for acquiring a digital baseband signal received by the navigation receiver and sampling the digital baseband signal to obtain a first sampling signal; the digital baseband signal comprises a time division spread spectrum signal with large Doppler frequency offset;
the first-stage searching module is used for traversing a preset first code phase searching range, performing code phase capturing on the time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase capturing result in the first code phase searching range, and obtaining a first pseudo code phase capturing result according to the maximum value of the pseudo code phase capturing result;
the searching range determining module is used for sampling the digital baseband signal after waiting for a preset time interval at an initial sampling point of the first sampling signal to obtain a second sampling signal, and obtaining a second code phase searching range according to the pre-estimated maximum Doppler corresponding to the digital baseband signal;
the second-stage searching module is used for traversing the second code phase searching range and performing code phase capturing on the time division pulse sequence corresponding to the second sampling signal to obtain a corresponding second pseudo code phase capturing result;
and the carrier Doppler estimation module is used for obtaining carrier Doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
acquiring a digital baseband signal received by a navigation receiver, and sampling the digital baseband signal to obtain a first sampling signal; the digital baseband signal comprises a time division spread spectrum signal with large Doppler frequency offset;
traversing a preset first code phase searching range, performing code phase capturing on the time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase capturing result in the first code phase searching range, and obtaining a first pseudo code phase capturing result according to the maximum value of the pseudo code phase capturing result;
after waiting for a preset time interval at an initial sampling point of the first sampling signal, sampling the digital baseband signal to obtain a second sampling signal, and calculating according to a pre-estimated maximum Doppler corresponding to the digital baseband signal to obtain a second code phase search range;
traversing the second code phase search range, and performing code phase acquisition on the time division pulse sequence corresponding to the second sampling signal to obtain a corresponding second pseudo code phase acquisition result;
and obtaining carrier Doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
acquiring a digital baseband signal received by a navigation receiver, and sampling the digital baseband signal to obtain a first sampling signal; the digital baseband signal comprises a time division spread spectrum signal with large Doppler frequency offset;
traversing a preset first code phase searching range, performing code phase capturing on the time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase capturing result in the first code phase searching range, and obtaining a first pseudo code phase capturing result according to the maximum value of the pseudo code phase capturing result;
after waiting for a preset time interval at an initial sampling point of the first sampling signal, sampling the digital baseband signal to obtain a second sampling signal, and calculating according to a pre-estimated maximum Doppler corresponding to the digital baseband signal to obtain a second code phase search range;
traversing the second code phase searching range, and performing code phase acquisition on the time division pulse sequence corresponding to the second sampling signal to obtain a corresponding second pseudo code phase acquisition result;
and obtaining carrier Doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result.
The two-stage rapid signal capturing method of the time division navigation signal under the large Doppler realizes the capturing of the Doppler of the signal by performing two times of serial time dimension searching on the received digital baseband signal, and particularly converts two-dimensional time-frequency searching into one-dimensional time searching and one-dimensional frequency estimation, when the one-dimensional time searching is performed, the pseudo code phase is captured twice, the first stage of capturing can perform rapid searching on the code phase of the signal in a larger time range under the condition that the signal has large Doppler frequency offset, the initial capturing of the code phase of the signal is completed, the second stage of capturing is based on the first stage of capturing, the time searching range of the signal can be greatly compressed, a more accurate phase searching result is obtained, after the two-stage pseudo code phase capturing is realized, the carrier Doppler is calculated according to the two-time pseudo code phase capturing results, the searching space of the Doppler of the signal is greatly compressed, and the capturing efficiency of the time division spread spectrum signal under the large Doppler frequency offset is improved.
Drawings
FIG. 1 is a flow chart illustrating a two-stage fast signal acquisition method for time-division navigation signals under large Doppler in one embodiment;
FIG. 2 is a schematic diagram of an embodiment of a time-division spread-spectrum signal;
FIG. 3 is a diagram illustrating correlation peak results for a first search and a second search, under an embodiment;
FIG. 4 is a block diagram of a two-stage fast signal acquisition device for time-division navigation signals under large Doppler in one embodiment;
FIG. 5 is a diagram of the internal structure of a computer device in one embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
In one embodiment, as shown in fig. 1, a two-stage fast signal acquisition method for time-division navigation signals under large doppler is provided, which includes the following steps:
step 102, acquiring a digital baseband signal received by the navigation receiver, and sampling the digital baseband signal to obtain a first sampling signal.
The GPS signal is a spread spectrum signal that has been subjected to direct sequence spread spectrum modulation and carrier modulation. For spread spectrum systems, acquisition means that the phase difference between the local reference code and the receiving code is less than one symbol width, the clock frequencies of the receiving and transmitting codes are basically consistent, and simultaneously, the carriers are mutually aligned to realize the synchronization of the input signal and the local signal. The coarse synchronization process of the pseudo code phase and the carrier frequency in the GPS system is pseudo code capturing, an effective pseudo code capturing method is the core of research of a high dynamic GPS receiver, and the performance of the system is improved by shortening the capturing time. When capturing the signal of the transmitting end, the receiver can confirm the successful capturing only after the pseudo code phase and the carrier Doppler frequency which are locally reproduced by the receiver are matched with the received signal, thereby realizing the bit synchronization of the received and transmitted signal. In the present invention, the digital baseband signal includes a time division spread spectrum signal with a large doppler frequency offset, the first sampling signal includes a multi-segment time division pulse signal, as shown in the schematic structural diagram of the time division spread spectrum signal shown in fig. 2, compared with the traditional acquisition of the spread spectrum signal, the time division spread spectrum signal adds a time division pulse sequence, which is a discontinuous signal,
Figure 452546DEST_PATH_IMAGE062
for the pulse period, the time-divided pulse sequence being of pulse width
Figure 276146DEST_PATH_IMAGE063
Has signal pulses therein
Figure 867664DEST_PATH_IMAGE064
There is no signal in time, so when the time division pulse sequence is subjected to code phase acquisition, the time division pulse sequence is multiplied by a PN code (pseudo random code), and a signal in the presence of a signal pulse can be processed.
And 104, traversing a preset first code phase search range, performing code phase capturing on the time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase capturing result in the first code phase search range, and obtaining a first pseudo code phase capturing result according to the maximum value of the pseudo code phase capturing result.
The first code phase search range is a preset signal capture time range, the duration of a part of time division spread spectrum signals with signals is short, and the Doppler tolerance range is large, so that only a pseudo code phase can be captured in a capture link, code phase capture is carried out on a time division pulse sequence corresponding to a first sampling signal to serve as first-stage search, a first pseudo code phase capture result is obtained through the first-stage search, pseudo code Doppler is obtained through two pseudo code phase capture results, carrier Doppler is further obtained, and rapid capture of the signals can be still maintained under the condition that a search space is greatly reduced.
And 106, after waiting for a preset time interval at the initial sampling point of the first sampling signal, sampling the digital baseband signal to obtain a second sampling signal, and calculating according to the maximum Doppler corresponding to the digital baseband signal estimated in advance to obtain a second code phase search range.
The principle of Doppler acquisition is to have a time interval of
Figure 968476DEST_PATH_IMAGE065
The received signal is captured twice, the pseudo code Doppler is calculated according to the code phase capture difference of the two times, then the pseudo code Doppler is converted to the carrier Doppler, the carrier Doppler is estimated by utilizing the difference result of the two code phase captures, the second code phase search range is calculated according to the maximum Doppler corresponding to the estimated digital baseband signal, and compared with the first-stage search, the search range of the second-stage search is greatly reduced.
And step 108, traversing the second code phase search range, and performing code phase acquisition on the time division pulse sequence corresponding to the second sampling signal to obtain a corresponding second pseudo code phase acquisition result.
And performing code phase acquisition on the time division pulse sequence corresponding to the second sampling signal as a second-stage search, and obtaining a second pseudo code phase acquisition result through the second-stage search. The second-level searching method is the same as the first-level searching method, but the searching range of the second-level searching is reduced, and the searching precision is improved.
And step 110, obtaining carrier Doppler according to the first pseudo code phase capturing result and the second pseudo code phase capturing result.
In a general acquisition method, because of many units for signal searching, the corresponding signal acquisition time is also long, which is not allowed in some receivers with high real-time requirement. The two-stage capturing algorithm provided by the invention can capture the pseudo code phase twice, converts two-dimensional time-frequency search into one-dimensional time search and one-dimensional frequency estimation, and greatly improves the efficiency of the algorithm.
In the two-stage fast signal capturing method for time division navigation signals under large Doppler, the signal Doppler is captured by performing two-time serial time dimension search on received digital baseband signals, specifically, two-dimensional time-frequency search is converted into one-dimensional time search and one-dimensional frequency estimation, when one-dimensional time search is performed, a pseudo code phase is captured twice, the first-stage capture can perform fast search on the code phase of the signal in a larger time range under the condition that the signal has large Doppler frequency offset, initial capture of the code phase of the signal is completed, the second-stage capture is based on the first-stage capture, the time search range of the signal can be greatly compressed, a more accurate phase search result is obtained, after two-time pseudo code phase capture is performed, carrier Doppler is calculated according to the two-time pseudo code phase capture results, the Doppler search space of the signal is greatly compressed, and the capturing efficiency of the time division spread spectrum signal under large Doppler frequency offset is improved.
In one embodiment, the mathematical model of the digital baseband signal is:
Figure 416774DEST_PATH_IMAGE066
wherein,
Figure 117883DEST_PATH_IMAGE067
for digital baseband signals received at the navigation receiver device,
Figure 44251DEST_PATH_IMAGE068
is the amplitude of the carrier wave and,
Figure 327465DEST_PATH_IMAGE069
in order to provide a satellite navigation message,
Figure 822031DEST_PATH_IMAGE070
in order to sample the points of interest,
Figure 620223DEST_PATH_IMAGE071
in order to delay the transmission of the signal,
Figure 241959DEST_PATH_IMAGE072
is the carrier-wave doppler of the signal,w(k) In the case of baseband noise, the noise is,
Figure 176417DEST_PATH_IMAGE073
for the purpose of the spreading codes in the navigation signal,
Figure 169781DEST_PATH_IMAGE074
in order to time-division pulse signals,
Figure 330635DEST_PATH_IMAGE075
is a symbol of an imaginary unit of a number,
Figure 598805DEST_PATH_IMAGE076
is the radio frequency carrier initial phase. In this embodiment, the signal acquisition process can be described as using a known spreading code sequence and time division pulse sequence of the navigation signal to receive the received signal
Figure 43562DEST_PATH_IMAGE077
In the course of a sequence
Figure 4565DEST_PATH_IMAGE078
And
Figure 387136DEST_PATH_IMAGE079
two parameters are estimated. Compared with the traditional capture of spread spectrum signals, the time division spread spectrum signals are added with time division pulse sequences, so that the capture of the signals by the time division pulse sequences is also needed in the capture process.
In one embodiment, traversing a preset first code phase search range, and performing code phase acquisition on a time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase acquisition result in the first code phase search range includes: performing incoherent accumulation on a plurality of coherent integration results obtained by performing coherent integration on each time division pulse of the first sampling signal and a local spread spectrum code corresponding to the current code phase to obtain an incoherent accumulation result corresponding to the current code phase; the current code phase is in a preset first code phase searching range; traversing the first code phase searching range to obtain a non-coherent accumulation result corresponding to each code phase in the first code phase searching range, and obtaining a pseudo code phase capturing result in the first code phase searching range according to the magnitude relation between each non-coherent accumulation result and a judgment threshold; performing incoherent accumulation on a plurality of coherent integration results obtained by coherently integrating each time division pulse of the first sampling signal and the local spreading code corresponding to the current code phase, wherein obtaining the incoherent accumulation result corresponding to the current code phase comprises: performing incoherent accumulation on a plurality of coherent integration results obtained by performing coherent integration on each time division pulse of the first sampling signal and the local spread spectrum code corresponding to the current code phase, wherein the incoherent accumulation result corresponding to the current code phase is obtained by:
Figure 662259DEST_PATH_IMAGE080
wherein,
Figure 836888DEST_PATH_IMAGE081
is a code phase of
Figure 323453DEST_PATH_IMAGE082
Corresponding non-coherent accumulation result,
Figure 849112DEST_PATH_IMAGE083
Is the phase of the baseband digital signal and,
Figure 334451DEST_PATH_IMAGE084
is a base-band complex signal and is,
Figure 363587DEST_PATH_IMAGE085
is as follows
Figure 604075DEST_PATH_IMAGE086
The position of the initial sampling point of each pulse signal,
Figure 476085DEST_PATH_IMAGE087
for local replication of phase is
Figure 889749DEST_PATH_IMAGE088
The spreading code of (a) is used,
Figure 648758DEST_PATH_IMAGE089
in order to be a short-time correlation length,
Figure 122464DEST_PATH_IMAGE090
the number of pulses corresponding to the first sampling signal.
In this embodiment, in the first level search, a selection is made
Figure 983235DEST_PATH_IMAGE091
The signal is captured in a time-division manner,
Figure 138273DEST_PATH_IMAGE091
the duration of the segment time signal is
Figure 142001DEST_PATH_IMAGE092
Figure 396396DEST_PATH_IMAGE093
And
Figure 118364DEST_PATH_IMAGE094
from time-hopping pulse patterns
Figure 998465DEST_PATH_IMAGE095
Determine that
Figure 59961DEST_PATH_IMAGE091
Segment time division signal and
Figure 875471DEST_PATH_IMAGE096
performing short-time correlation and post-accumulation operation to obtain
Figure 960101DEST_PATH_IMAGE097
Corresponding non-coherent accumulation result
Figure 253680DEST_PATH_IMAGE098
Short time correlation and pulse width
Figure 169683DEST_PATH_IMAGE099
Remain in agreement, i.e.
Figure 782192DEST_PATH_IMAGE100
Figure 213173DEST_PATH_IMAGE101
Is the signal sampling rate. When in use
Figure 451388DEST_PATH_IMAGE102
Figure 18635DEST_PATH_IMAGE103
A preset decision threshold), the code phase is successfully captured, the corresponding local spread spectrum code phase is the captured result, otherwise, the search is continued until the end or the related result exceeds the thresholdTh
In one embodiment, the step of obtaining the pseudo code phase acquisition result in the first code phase search range according to the magnitude relationship between each non-coherent accumulation result and the decision threshold includes: when the incoherent accumulation result is greater than the judgment threshold, successfully acquiring the code phase, and obtaining an acquisition result according to the code phase corresponding to the current local spread spectrum code; when the incoherent accumulation result is less than or equal to the decision threshold, the code phase acquisition fails, a preset pseudo code phase search interval slides backwards in a first code phase search range, the incoherent accumulation result corresponding to the code phase after sliding is calculated, the process is iterated until each code phase in the first code phase search range is searched, the iteration is stopped, and the acquisition result in the first code phase search range is output; the first pseudo code phase acquisition result is:
Figure 379210DEST_PATH_IMAGE104
wherein,
Figure 422121DEST_PATH_IMAGE105
for the first pseudo-code phase acquisition result,
Figure 57502DEST_PATH_IMAGE106
in order to be able to determine the code phase,
Figure 620201DEST_PATH_IMAGE107
is as follows
Figure 213993DEST_PATH_IMAGE108
Get maximum time corresponding
Figure 980086DEST_PATH_IMAGE109
The value of (a) is set to (b),
Figure 684737DEST_PATH_IMAGE110
is a code phase of
Figure 429839DEST_PATH_IMAGE111
The result of the corresponding non-coherent accumulation,
Figure 804320DEST_PATH_IMAGE112
the interval is searched for the pseudo-code phase,
Figure 431610DEST_PATH_IMAGE113
the number of cells is searched for the pseudo code phase. In this embodiment, the pseudo code phase capturing result of the signal can be obtained by traversing all the possibilities of the pseudo code initial phases and taking the phase corresponding to the maximum value of all the search results.
In one embodiment, the second pseudo-code phase acquisition result is:
Figure 533427DEST_PATH_IMAGE114
wherein,
Figure 195353DEST_PATH_IMAGE115
for the second pseudo-code phase acquisition result,
Figure 6314DEST_PATH_IMAGE116
is as follows
Figure 120901DEST_PATH_IMAGE117
Get maximum time corresponding
Figure 105037DEST_PATH_IMAGE118
The value of (a) is set to (b),
Figure 265146DEST_PATH_IMAGE119
is a code phase of
Figure 371642DEST_PATH_IMAGE120
The result of the corresponding non-coherent accumulation,
Figure 848891DEST_PATH_IMAGE121
for the first pseudo-code phase acquisition result,
Figure 433456DEST_PATH_IMAGE122
for the second code phase search range,
Figure 7657DEST_PATH_IMAGE123
searching an interval for a pseudo code phase; obtaining a second code phase search by maximum Doppler calculation corresponding to the pre-estimated digital baseband signalThe cable range includes: and calculating according to the maximum Doppler corresponding to the pre-estimated digital baseband signal to obtain a second code phase search range as follows:
Figure 144109DEST_PATH_IMAGE124
wherein,
Figure 233288DEST_PATH_IMAGE125
for the second code phase search range,
Figure 496910DEST_PATH_IMAGE126
for the maximum doppler corresponding to the digital baseband signal,
Figure 722355DEST_PATH_IMAGE127
Figure 108337DEST_PATH_IMAGE128
is the radio frequency of the signal and,
Figure 310911DEST_PATH_IMAGE129
in order to be at the pseudo-code rate,
Figure 502857DEST_PATH_IMAGE130
which represents a rounding-up operation on the upper part,
Figure 458175DEST_PATH_IMAGE131
in the form of a time interval,
Figure 77375DEST_PATH_IMAGE132
the interval is searched for a pseudo code phase. In this embodiment, the second code phase search range includes
Figure 78829DEST_PATH_IMAGE133
A code phase.
In one embodiment, obtaining carrier doppler based on the first and second pseudo-code phase acquisition results comprises: according to the first pseudo code phase capturing result and the second pseudo code phase capturing result, carrier Doppler is obtained as follows:
Figure 667943DEST_PATH_IMAGE134
wherein,
Figure 133559DEST_PATH_IMAGE135
in the form of a carrier wave doppler, the doppler,
Figure 799027DEST_PATH_IMAGE136
is the ratio of the signal radio frequency to the pseudo code rate,
Figure 818935DEST_PATH_IMAGE137
the result is captured for the second pseudo-code phase,
Figure 555947DEST_PATH_IMAGE138
for the first pseudo-code phase acquisition result,
Figure 236589DEST_PATH_IMAGE139
are time intervals.
In one embodiment, the method further comprises: the estimation precision of the carrier Doppler obtained according to the time interval and the pseudo code phase search interval is as follows:
Figure 197592DEST_PATH_IMAGE140
wherein,
Figure 845742DEST_PATH_IMAGE141
is the accuracy of the estimation of the carrier doppler,
Figure 917603DEST_PATH_IMAGE142
the interval is searched for the pseudo-code phase,
Figure 29916DEST_PATH_IMAGE143
is the ratio of the radio frequency of the signal to the rate of the pseudo code; adjusting time intervals and pseudo code phase search intervals to adjust acquisition accuracy of carrier doppler. In the embodiment, the estimation of carrier Doppler is realized by using the difference result of two times of code phase acquisition, and the carrier Doppler is prolonged
Figure 286454DEST_PATH_IMAGE144
Or reduce
Figure 546534DEST_PATH_IMAGE145
All can improve the capture precision of the signal carrier frequency
Figure 297452DEST_PATH_IMAGE146
Figure 326588DEST_PATH_IMAGE147
A code piece,
Figure 567077DEST_PATH_IMAGE148
The accuracy of the carrier frequency estimation is 620Hz.
In a specific embodiment, the present invention provides a two-stage capturing algorithm for a large doppler time division signal, which comprises the following specific steps, as shown in the following table, first obtaining a digital baseband signal to be captured, sampling the digital baseband signal, then performing a first stage search on a first sampled signal, where the total number of code phases searched in a first code phase search range is
Figure 940551DEST_PATH_IMAGE149
Judging whether the time search is finished or not by the pseudo code phase, if so, acquiring carrier Doppler, and if not, acquiring successfully, otherwise, acquiring unsuccessfully; if the time search of the first-stage search is not finished, sliding the pseudo code phase search interval of the local signal backwards, carrying out short-time correlation accumulation operation on the first sampling signal, the local spread spectrum code and the time-hopping pulse pattern to obtain a correlation value corresponding to each pseudo code phase, comparing the correlation value with a threshold, if the correlation value is larger than the threshold, recording the phase at the moment, otherwise, judging whether the time search is finished again, outputting a first pseudo code phase capture result according to the maximum captured phase, and waiting for the acquisition
Figure 354215DEST_PATH_IMAGE150
And sampling the digital baseband signal after time to obtain a second sampling signal, calculating a second code phase search range in the first code phase search range according to the pre-estimated maximum Doppler, performing second-stage search on the second sampling signal in the second code phase search range to obtain a second pseudo code phase acquisition result, and finally calculating to obtain carrier Doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result to finish signal acquisition.
TABLE 1 two-stage acquisition algorithm for large Doppler time division signals
Figure 378803DEST_PATH_IMAGE151
In a specific embodiment, as shown in fig. 3, a schematic diagram of a correlation peak result of a first-stage search and a second-stage search is provided, as can be seen from fig. 3, a range of the second-stage search is greatly reduced compared with a range of the first-stage search, and by capturing two stages of code phases, doppler frequency offset of a signal can be directly calculated without searching doppler of the signal, which means that the method of the present invention can greatly simplify complexity of capturing a time-division navigation signal with large doppler frequency offset.
It should be understood that, although the steps in the flowchart of fig. 1 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least a portion of the steps in fig. 1 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternately with other steps or at least a portion of the sub-steps or stages of other steps.
In one embodiment, as shown in fig. 4, there is provided a two-stage fast signal acquisition apparatus for time-division navigation signals under large doppler, comprising: a signal acquisition module 402, a first stage search module 404, a search range determination module 406, a second stage search module 408, and a carrier doppler estimation module 410, wherein:
a signal obtaining module 402, configured to obtain a digital baseband signal received by the navigation receiver, and sample the digital baseband signal to obtain a first sampling signal; the digital baseband signal comprises a time division spread spectrum signal with large Doppler frequency offset;
the first-stage searching module 404 is configured to traverse a preset first code phase searching range, perform code phase acquisition on a time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase acquisition result within the first code phase searching range, and obtain a first pseudo code phase acquisition result according to a maximum value of the pseudo code phase acquisition result;
a search range determining module 406, configured to sample the digital baseband signal after waiting for a preset time interval at an initial sampling point of the first sampling signal to obtain a second sampling signal, and obtain a second code phase search range according to a pre-estimated maximum doppler corresponding to the digital baseband signal;
the second-stage searching module 408 is configured to traverse a second code phase searching range, and perform code phase acquisition on the time division pulse sequence corresponding to the second sampling signal to obtain a corresponding second pseudo code phase acquisition result;
and a carrier doppler estimation module 410, configured to obtain carrier doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result.
In one embodiment, the signal acquisition module 402 is further configured to mathematically model the digital baseband signal as:
Figure 586930DEST_PATH_IMAGE152
wherein,
Figure 680657DEST_PATH_IMAGE153
receiving digits for navigationThe baseband signal is then used to generate a baseband signal,
Figure 898012DEST_PATH_IMAGE154
is the amplitude of the carrier wave and,
Figure 511527DEST_PATH_IMAGE155
in order to provide a satellite navigation message,
Figure 359397DEST_PATH_IMAGE156
in order to sample the points of interest,
Figure 815786DEST_PATH_IMAGE157
in order to delay the transmission of the signal,
Figure 722651DEST_PATH_IMAGE158
is the carrier-wave doppler of the signal,w(k) In the case of baseband noise, the noise is,
Figure 456251DEST_PATH_IMAGE159
for the purpose of the spreading codes in the navigation signal,
Figure 537340DEST_PATH_IMAGE160
in order to time-division pulse signals,
Figure 418708DEST_PATH_IMAGE161
is a symbol of an imaginary unit of numbers,
Figure 102499DEST_PATH_IMAGE162
is the radio frequency carrier initial phase.
In one embodiment, the first-stage search module 404 is further configured to perform non-coherent accumulation on a plurality of coherent integration results obtained after performing coherent integration on each time division pulse of the first sampling signal and the local spreading code corresponding to the current code phase to obtain a non-coherent accumulation result corresponding to the current code phase; the current code phase is in a preset first code phase searching range; and traversing the first code phase searching range to obtain a non-coherent accumulation result corresponding to each code phase in the first code phase searching range, and obtaining a pseudo code phase capturing result in the first code phase searching range according to the magnitude relation between each non-coherent accumulation result and the decision threshold.
In one embodiment, the first-stage searching module 404 is further configured to, when the incoherent accumulation result is greater than the decision threshold, successfully acquire the code phase, and obtain an acquisition result according to the code phase corresponding to the current local spreading code; and when the incoherent accumulation result is less than or equal to the judgment threshold, failing to acquire the code phase, sliding a preset pseudo code phase search interval backwards in the first code phase search range, calculating the incoherent accumulation result corresponding to the slid code phase, iterating the process until each code phase in the first code phase search range is searched, stopping iteration, and outputting the acquisition result in the first code phase search range.
In one embodiment, the first-stage search module 404 is further configured to perform non-coherent accumulation on a plurality of coherent integration results obtained after performing coherent integration on each time division pulse of the first sampling signal and the local spreading code corresponding to the current code phase, and obtain a non-coherent accumulation result corresponding to the current code phase as:
Figure 815240DEST_PATH_IMAGE163
wherein,
Figure 942596DEST_PATH_IMAGE164
is a code phase of
Figure 373578DEST_PATH_IMAGE165
The result of the corresponding non-coherent accumulation,
Figure 96945DEST_PATH_IMAGE166
is the phase of the baseband digital signal and,
Figure 867455DEST_PATH_IMAGE167
is a base-band complex signal and is,
Figure 24767DEST_PATH_IMAGE168
is as follows
Figure 818411DEST_PATH_IMAGE169
The position of the initial sampling point of each pulse signal,
Figure 719371DEST_PATH_IMAGE170
for local replication of phase is
Figure 531338DEST_PATH_IMAGE171
The spreading code of (a) is used,
Figure 62813DEST_PATH_IMAGE172
in order to be a short-time correlation length,
Figure 202807DEST_PATH_IMAGE173
the number of pulses corresponding to the first sampling signal.
In one embodiment, the search range determining module 406 is further configured to obtain a second code phase search range according to the pre-estimated maximum doppler corresponding to the digital baseband signal, where the second code phase search range is:
Figure 782825DEST_PATH_IMAGE174
wherein,
Figure 324664DEST_PATH_IMAGE175
for the second code phase search range,
Figure 449877DEST_PATH_IMAGE176
for a pre-estimated maximum doppler corresponding to the digital baseband signal,
Figure 77168DEST_PATH_IMAGE177
Figure 195297DEST_PATH_IMAGE178
is the radio frequency of the signal and,
Figure 857222DEST_PATH_IMAGE179
in order to be at the pseudo-code rate,
Figure 730500DEST_PATH_IMAGE180
which represents a rounding-up operation on the upper part,
Figure 969720DEST_PATH_IMAGE181
the time interval is a time interval of,
Figure 16174DEST_PATH_IMAGE182
the interval is searched for the pseudo code phase.
In one embodiment, the carrier doppler estimation module 410 is further configured to obtain, according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result, carrier doppler as:
Figure 407972DEST_PATH_IMAGE183
wherein,
Figure 452151DEST_PATH_IMAGE184
is a carrier wave doppler (doppler) wave,
Figure 788455DEST_PATH_IMAGE185
the ratio of the signal rf frequency to the pseudo code rate,
Figure 264698DEST_PATH_IMAGE186
for the second pseudo-code phase acquisition result,
Figure 635636DEST_PATH_IMAGE121
the result is captured for the first pseudo-code phase,
Figure 522821DEST_PATH_IMAGE187
are time intervals.
In one embodiment, the second stage search module 408 is further configured to obtain the second pseudo-code phase acquisition result as:
Figure 815262DEST_PATH_IMAGE188
wherein,
Figure 937939DEST_PATH_IMAGE189
for the second pseudo-code phase acquisition result,
Figure 819176DEST_PATH_IMAGE190
is that when
Figure 1896DEST_PATH_IMAGE191
Get maximum time corresponding
Figure 516053DEST_PATH_IMAGE192
The value of (a) is,
Figure 317787DEST_PATH_IMAGE193
is a code phase of
Figure 928897DEST_PATH_IMAGE194
The result of the corresponding non-coherent accumulation,
Figure 902758DEST_PATH_IMAGE195
for the first pseudo-code phase acquisition result,
Figure 700949DEST_PATH_IMAGE196
for the second code phase search range,
Figure 634270DEST_PATH_IMAGE197
the interval is searched for a pseudo code phase.
In one embodiment, the first stage search module 404 is further configured to obtain the following result:
Figure 709674DEST_PATH_IMAGE198
wherein,
Figure 499775DEST_PATH_IMAGE199
the result is captured for the first pseudo-code phase,
Figure 909897DEST_PATH_IMAGE200
is the phase of the code and is,
Figure 443646DEST_PATH_IMAGE201
is that when
Figure 373556DEST_PATH_IMAGE202
Taking maximum time to correspond to
Figure 272242DEST_PATH_IMAGE203
The value of (a) is set to (b),
Figure 45026DEST_PATH_IMAGE204
is a code phase of
Figure 8565DEST_PATH_IMAGE205
The result of the corresponding non-coherent accumulation,
Figure 917615DEST_PATH_IMAGE206
the interval is searched for the pseudo-code phase,
Figure 924885DEST_PATH_IMAGE207
the number of cells is searched for the pseudo code phase.
In one embodiment, the estimation accuracy for carrier doppler obtained from the time interval and the pseudo-code phase search interval is:
Figure 122649DEST_PATH_IMAGE208
wherein,
Figure 998201DEST_PATH_IMAGE209
for the accuracy of the estimation of the carrier doppler,
Figure 886391DEST_PATH_IMAGE210
the interval is searched for the pseudo-code phase,
Figure 454776DEST_PATH_IMAGE211
is the ratio of the radio frequency of the signal to the rate of the pseudo code; adjusting time intervals and pseudo code phase search intervals to adjust carrier dopplerThe capture accuracy of (2).
For specific limitations of the two-stage fast signal capturing apparatus for time division navigation signals under large doppler, reference may be made to the above limitations of the two-stage fast signal capturing method for time division navigation signals under large doppler, which are not described herein again. All or part of each module in the two-stage rapid signal acquisition device for the time division navigation signal under the large Doppler can be 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 is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 5. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operating system and the computer program to run on the non-volatile storage medium. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to realize a two-stage fast signal acquisition method of time division navigation signals under large Doppler. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 5 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In an embodiment, a computer device is provided, comprising a memory storing a computer program and a processor implementing the steps of the method in the above embodiments when the processor executes the computer program.
In an embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which computer program, when being executed by a processor, carries out the steps of the method in the above-mentioned embodiments.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database or other medium used in the embodiments provided herein can include non-volatile and/or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), rambus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
All possible combinations of the technical features in the above embodiments may not be described for the sake of brevity, but should be considered as being within the scope of the present disclosure as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is specific and detailed, but not to be understood as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent application shall be subject to the appended claims.

Claims (10)

1. A two-stage fast signal acquisition method for time division navigation signals under large Doppler is characterized by comprising the following steps:
acquiring a digital baseband signal received by a navigation receiver, and sampling the digital baseband signal to obtain a first sampling signal; the digital baseband signal comprises a time division spread spectrum signal with large Doppler frequency offset;
traversing a preset first code phase searching range, performing code phase capturing on a time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase capturing result in the first code phase searching range, and obtaining a first pseudo code phase capturing result according to the maximum value of the pseudo code phase capturing result;
after waiting for a preset time interval at an initial sampling point of the first sampling signal, sampling the digital baseband signal to obtain a second sampling signal, and calculating according to a pre-estimated maximum Doppler corresponding to the digital baseband signal to obtain a second code phase search range;
traversing the second code phase search range, and performing code phase acquisition on the time division pulse sequence corresponding to the second sampling signal to obtain a corresponding second pseudo code phase acquisition result;
and obtaining carrier Doppler according to the first pseudo code phase acquisition result and the second pseudo code phase acquisition result.
2. The method of claim 1, wherein the mathematical model of the digital baseband signal is:
Figure DEST_PATH_IMAGE001
wherein,
Figure DEST_PATH_IMAGE002
to navigate the digital baseband signals received by the receiving device,
Figure DEST_PATH_IMAGE003
is the amplitude of the carrier wave and,
Figure DEST_PATH_IMAGE004
in order to provide a satellite navigation message,
Figure DEST_PATH_IMAGE005
in order to sample the points of interest,
Figure DEST_PATH_IMAGE006
in order to delay the transmission of the signal,
Figure DEST_PATH_IMAGE007
is the carrier-wave doppler of the signal,w(k) In the case of baseband noise, the noise is,
Figure DEST_PATH_IMAGE008
for the purpose of the spreading codes in the navigation signal,
Figure DEST_PATH_IMAGE009
in order to time-division pulse signals,
Figure DEST_PATH_IMAGE010
is a symbol of an imaginary unit of a number,
Figure DEST_PATH_IMAGE011
is the radio frequency carrier initial phase.
3. The method of claim 1, wherein traversing a preset first code phase search range, and performing code phase acquisition on a time division pulse sequence corresponding to the first sampling signal to obtain a pseudo code phase acquisition result in the first code phase search range comprises:
performing non-coherent accumulation on a plurality of coherent integration results obtained after coherent integration is performed on each time division pulse of the first sampling signal and a local spread spectrum code corresponding to the current code phase to obtain a non-coherent accumulation result corresponding to the current code phase; the current code phase is in a preset first code phase searching range;
and traversing the first code phase searching range to obtain a non-coherent accumulation result corresponding to each code phase in the first code phase searching range, and obtaining a pseudo code phase capturing result in the first code phase searching range according to the magnitude relation between each non-coherent accumulation result and a decision threshold.
4. The method of claim 3, wherein the step of obtaining the pseudo code phase acquisition result in the first code phase search range according to the magnitude relationship between each non-coherent accumulation result and the decision threshold comprises:
when the incoherent accumulation result is larger than a judgment threshold, successfully capturing the code phase, and obtaining a capture result according to the code phase corresponding to the current local spread spectrum code;
and when the incoherent accumulation result is smaller than or equal to the decision threshold, failing to acquire the code phase, sliding a preset pseudo code phase search interval backwards in the first code phase search range, calculating the incoherent accumulation result corresponding to the slid code phase, iterating the process until each code phase in the first code phase search range is searched, stopping iteration, and outputting the acquisition result in the first code phase search range.
5. The method of claim 3, wherein the performing non-coherent accumulation on a plurality of coherent integration results obtained by performing coherent integration on each time division pulse of the first sampling signal and a local spreading code corresponding to a current code phase to obtain a non-coherent accumulation result corresponding to the current code phase comprises:
performing incoherent accumulation on a plurality of coherent integration results obtained by performing coherent integration on each time division pulse of the first sampling signal and the local spread spectrum code corresponding to the current code phase, wherein the incoherent accumulation result corresponding to the current code phase is obtained by:
Figure DEST_PATH_IMAGE012
wherein,
Figure DEST_PATH_IMAGE013
is a code phase of
Figure DEST_PATH_IMAGE014
The result of the corresponding non-coherent accumulation,
Figure DEST_PATH_IMAGE015
is the phase of the baseband digital signal and,
Figure DEST_PATH_IMAGE016
is a base-band complex signal and is,
Figure DEST_PATH_IMAGE017
is as follows
Figure DEST_PATH_IMAGE018
The position of the initial sampling point of each pulse signal,
Figure DEST_PATH_IMAGE019
for local replication of phase is
Figure DEST_PATH_IMAGE020
The spreading code of (a) is used,
Figure DEST_PATH_IMAGE021
is related for a short timeThe length of the first and second support members,
Figure DEST_PATH_IMAGE022
the number of pulses corresponding to the first sampling signal.
6. The method of claim 1, wherein the calculating a second code phase search range according to the pre-estimated maximum doppler corresponding to the digital baseband signal comprises:
calculating according to the pre-estimated maximum Doppler corresponding to the digital baseband signal to obtain a second code phase search range as follows:
Figure DEST_PATH_IMAGE023
wherein,
Figure DEST_PATH_IMAGE024
for the second code phase search range,
Figure DEST_PATH_IMAGE025
for the maximum doppler corresponding to the pre-estimated digital baseband signal,
Figure DEST_PATH_IMAGE026
Figure DEST_PATH_IMAGE027
is the radio frequency of the signal and,
Figure DEST_PATH_IMAGE028
in order to be at the pseudo-code rate,
Figure DEST_PATH_IMAGE029
which represents a rounding-up operation on the upper part,
Figure DEST_PATH_IMAGE030
in the form of a time interval,
Figure DEST_PATH_IMAGE031
the interval is searched for a pseudo code phase.
7. The method of claim 1, wherein obtaining carrier doppler based on the first and second pseudo-code phase acquisition results comprises:
according to the first pseudo code phase capturing result and the second pseudo code phase capturing result, carrier Doppler is obtained as follows:
Figure DEST_PATH_IMAGE032
wherein,
Figure DEST_PATH_IMAGE033
is a carrier wave doppler (doppler) wave,
Figure DEST_PATH_IMAGE034
is the ratio of the signal radio frequency to the pseudo code rate,
Figure DEST_PATH_IMAGE035
for the second pseudo-code phase acquisition result,
Figure DEST_PATH_IMAGE036
the result is captured for the first pseudo-code phase,
Figure DEST_PATH_IMAGE037
are time intervals.
8. The method of any of claims 1 or 7, wherein the second pseudo-code phase acquisition result is:
Figure DEST_PATH_IMAGE038
wherein,
Figure DEST_PATH_IMAGE039
for the second pseudo-code phase acquisition result,
Figure DEST_PATH_IMAGE040
is as follows
Figure DEST_PATH_IMAGE041
Get maximum time corresponding
Figure DEST_PATH_IMAGE042
The value of (a) is,
Figure DEST_PATH_IMAGE043
is a code phase of
Figure DEST_PATH_IMAGE044
The result of the corresponding non-coherent accumulation,
Figure DEST_PATH_IMAGE045
for the first pseudo-code phase acquisition result,
Figure DEST_PATH_IMAGE046
for the second code phase search range,
Figure DEST_PATH_IMAGE047
the interval is searched for the pseudo code phase.
9. The method of any of claims 1 or 7, wherein the first pseudo-code phase acquisition result is:
Figure DEST_PATH_IMAGE048
wherein,
Figure DEST_PATH_IMAGE049
for the first pseudo-code phase acquisition result,
Figure DEST_PATH_IMAGE050
in order to be able to determine the code phase,
Figure DEST_PATH_IMAGE051
is as follows
Figure DEST_PATH_IMAGE052
Get maximum time corresponding
Figure DEST_PATH_IMAGE053
The value of (a) is,
Figure DEST_PATH_IMAGE054
is a code phase of
Figure DEST_PATH_IMAGE055
The result of the corresponding non-coherent accumulation,
Figure DEST_PATH_IMAGE056
the interval is searched for the pseudo-code phase,
Figure DEST_PATH_IMAGE057
the number of cells is searched for the pseudo code phase.
10. The method of claim 7, further comprising:
and obtaining the carrier Doppler estimation precision according to the time interval and the pseudo code phase search interval as follows:
Figure DEST_PATH_IMAGE058
wherein,
Figure DEST_PATH_IMAGE059
for the accuracy of the estimation of the carrier doppler,
Figure DEST_PATH_IMAGE060
the interval is searched for the pseudo-code phase,
Figure DEST_PATH_IMAGE061
is the ratio of the radio frequency of the signal to the rate of the pseudo code;
adjusting the time interval and the pseudo code phase search interval to adjust an acquisition accuracy of carrier doppler.
CN202211702895.4A 2022-12-29 2022-12-29 Two-stage rapid signal capture method of time division navigation signal under large Doppler Active CN115685271B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211702895.4A CN115685271B (en) 2022-12-29 2022-12-29 Two-stage rapid signal capture method of time division navigation signal under large Doppler

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211702895.4A CN115685271B (en) 2022-12-29 2022-12-29 Two-stage rapid signal capture method of time division navigation signal under large Doppler

Publications (2)

Publication Number Publication Date
CN115685271A true CN115685271A (en) 2023-02-03
CN115685271B CN115685271B (en) 2023-03-17

Family

ID=85055686

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211702895.4A Active CN115685271B (en) 2022-12-29 2022-12-29 Two-stage rapid signal capture method of time division navigation signal under large Doppler

Country Status (1)

Country Link
CN (1) CN115685271B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6735243B1 (en) * 2000-11-08 2004-05-11 Nokia Corporation Method, apparatus and system for fast acquisition of a spread spectrum signal
JP2006254500A (en) * 2006-05-29 2006-09-21 Sony Corp Correlation detecting apparatus and method and receiver
CN102520423A (en) * 2011-12-06 2012-06-27 成都金本华科技有限公司 Circuit for rapidly capturing long period pseudo random spread spectrum code of satellite navigation receiver and capturing method thereof
WO2013115481A1 (en) * 2012-01-30 2013-08-08 한국과학기술원 Compressed sensing method and device for quickly obtaining global navigation satellite system (gnss) and spread spectrum signals
CN103592664A (en) * 2013-10-17 2014-02-19 中国科学院光电研究院 Coarse acquisition and fine acquisition spread spectrum signal synchronization method
CN105049081A (en) * 2015-06-05 2015-11-11 重庆大学 Long-code spread spectrum signal rapid capturing method adaptive to high dynamic environment
CN106291615A (en) * 2016-07-28 2017-01-04 西安空间无线电技术研究所 A kind of two benches catching method of high dynamic Doppler frequency deviation
CN106707309A (en) * 2016-11-23 2017-05-24 北京自动化控制设备研究所 Precise trapping method for carrier and pseudo code two-dimensional searching
CN111399004A (en) * 2020-04-07 2020-07-10 北京理工大学 High-dynamic high-sensitivity GNSS signal capturing method
CN113972929A (en) * 2021-10-26 2022-01-25 上海无线电设备研究所 Method for capturing spread spectrum signal under high dynamic Doppler

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6735243B1 (en) * 2000-11-08 2004-05-11 Nokia Corporation Method, apparatus and system for fast acquisition of a spread spectrum signal
JP2006254500A (en) * 2006-05-29 2006-09-21 Sony Corp Correlation detecting apparatus and method and receiver
CN102520423A (en) * 2011-12-06 2012-06-27 成都金本华科技有限公司 Circuit for rapidly capturing long period pseudo random spread spectrum code of satellite navigation receiver and capturing method thereof
WO2013115481A1 (en) * 2012-01-30 2013-08-08 한국과학기술원 Compressed sensing method and device for quickly obtaining global navigation satellite system (gnss) and spread spectrum signals
CN103592664A (en) * 2013-10-17 2014-02-19 中国科学院光电研究院 Coarse acquisition and fine acquisition spread spectrum signal synchronization method
CN105049081A (en) * 2015-06-05 2015-11-11 重庆大学 Long-code spread spectrum signal rapid capturing method adaptive to high dynamic environment
CN106291615A (en) * 2016-07-28 2017-01-04 西安空间无线电技术研究所 A kind of two benches catching method of high dynamic Doppler frequency deviation
CN106707309A (en) * 2016-11-23 2017-05-24 北京自动化控制设备研究所 Precise trapping method for carrier and pseudo code two-dimensional searching
CN111399004A (en) * 2020-04-07 2020-07-10 北京理工大学 High-dynamic high-sensitivity GNSS signal capturing method
CN113972929A (en) * 2021-10-26 2022-01-25 上海无线电设备研究所 Method for capturing spread spectrum signal under high dynamic Doppler

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GUO YICHAO等: ""A code Doppler compensation algorithm in acquisition for high dynamic spread spectrum signals"" *
WANG LE等: ""High Dynamic Spread Spectrum Signal Acquisition Algorithm based on Delay Autocorrelation"" *
刘晓明: ""高动态环境下长码扩频信号快捕算法"" *
毕卫红: ""高动态下长周期扩频信号快速捕获技术的研究"" *

Also Published As

Publication number Publication date
CN115685271B (en) 2023-03-17

Similar Documents

Publication Publication Date Title
US7567636B2 (en) Method and apparatus for performing signal correlation using historical correlation data
US8472503B2 (en) Method and apparatus for performing frequency synchronization
US8351486B2 (en) Parallel correlator implementation using hybrid correlation in spread-spectrum communication
CN100580477C (en) GPS position measuring method and GPS position measuring system
US8634504B2 (en) Correlation calculating method, satellite signal capturing method, and correlation calculating circuit
JP2016524710A (en) Receiver and method for direct sequence spread spectrum signals
CN104360357A (en) Quick Beidou satellite signal capturing method and system based on circulation mode
EP1749349A1 (en) Performing an acquisition in a receiver
CN105553506B (en) A kind of quick capturing method and device of long code spread-spectrum signal
CN104065397A (en) Method and device for synchronously capturing pseudo codes in real time
CN106291610B (en) The compression correlation module and implementation method of GNSS signal compression acquisition equipment
JP2007520100A (en) GPS receiver using differential correlation
CN102798871B (en) Pseudo code capturing method and device based on pseudo code reconstruction
JP2004340855A (en) Method and system of satellite positioning
CN106291618B (en) The recovery module and implementation method of GNSS signal compression acquisition equipment
US7995682B2 (en) Method and apparatus for performing signal processing using historical correlation data
JP2001223672A (en) Synchronization method for receiver, positioning system, receiver and electronic device
CN115685271B (en) Two-stage rapid signal capture method of time division navigation signal under large Doppler
CN109633707B (en) Variable coefficient matched filtering rapid capturing method based on pre-averaging processing
CN104993844B (en) A kind of frequency domain search method and device
van der Merwe et al. Efficient multi-hypothesis tests for extended coherent acquisition using modified linear feedback shift register secondary codes
CN100438360C (en) Method and equipment for realizing parallel correlator utilizing block integral in spread spectrum communication
CN112764063A (en) Method for realizing capture processing and receiver
US20040141574A1 (en) Determination of the code phase between a code modulated signal and a replica code sequence
Wei et al. Simulation and analysis of GPS software receiver

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