CN105182376A - GNSS signal mixer-free parallel code phase capturing method - Google Patents
GNSS signal mixer-free parallel code phase capturing method Download PDFInfo
- Publication number
- CN105182376A CN105182376A CN201510637983.4A CN201510637983A CN105182376A CN 105182376 A CN105182376 A CN 105182376A CN 201510637983 A CN201510637983 A CN 201510637983A CN 105182376 A CN105182376 A CN 105182376A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- sequence
- right arrow
- rightarrow
- gnss signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a GNSS signal mixer-free parallel code phase capturing method. The method conducts FFT on signals received at N points, the transformed sequence is recorded as y1<right arrow>, and N is a power of 2. The method then conducts FFT of local code signals at N points, obtains a conjugated sequence of the transformed sequence and records the obtained sequence as y2<right arrow>. The method rotates the y2<right arrow> right for m bit positions, multiplies the y2<right arrow> with corresponding points of the y1<right arrow> and obtains a sequence y3<right arrow>, conducts inverse FFT on the y3<right arrow>, and outputs a sequence y4<right arrow>. The y4<right arrow> is rotated left for m bit positions, and an output sequence Sm<right arrow> is a coherent integration result of N different code phases at the mth frequency point. The maximum value of the envelope of the Sm<right arrow> is compared with a threshold. If the maximum value exceeds the threshold, a signal is considered to exist, and a corresponding code phase and Doppler frequency are output. The method can complete Doppler frequency searching with no need of frequency conversion of the signal first. Besides, the search precision is high, the calculation is simple, and the method is easy for implementation.
Description
Technical field
The present invention relates to a kind of GNSS signal to exempt from mixing and to walk abreast phase acquisition method.
Background technology
Common GNSS signal catching method comprises: serial search methods, parallel frequency search method and parallel code phase search method.Serial search methods searches for different code phases and Doppler frequency successively, and implementation complexity is low, but the search time needed is longer; Parallel frequency search method, has walked abreast the search FFT of different frequency, is equivalent to utilize FFT to carry out spectrum analysis; Parallel code phase search method, has come the coherent integration operation convolution theorem of different code phase.Traditional code phase parallel search method needs first to carry out frequency conversion to signal, and to complete the search of Doppler frequency, computation complexity is high.
The process of parallel code phase search can represent as shown in Figure 1.For convenience of description, first do some hypothesis without loss of generality: suppose that coherent integration is counted as N (in order to carry out FFT computing, limiting N is the power of 2), coherent integration time NT
sfor the one-period of pseudo-random code; Note baseband signal is r
n, in corresponding coherent integration time, the baseband sampling signal phasor of input is
the local pseudo-random code produced of note is c
n, corresponding local pseudo-random code vector is
remember that the coherent integration results under different code phase is s
n, corresponding coherent integration results vector is
represent the coherent integration results in N kind code phase situation respectively; Note is in coherent integration process, and the local carrier vector that the frequency of use is f is:
represent that frequency is the local carrier signal vector of f, be used for being multiplied with baseband signal, to eliminate Doppler frequency remaining in baseband signal, T
srepresent signal sampling interval, N represents that coherent integration is counted.
Summary of the invention
Technical matters to be solved by this invention is, not enough for prior art, provides a kind of GNSS signal to exempt from mixing and to walk abreast phase acquisition method.
For solving the problems of the technologies described above, the technical solution adopted in the present invention is: a kind of GNSS signal is exempted from mixing and to be walked abreast phase acquisition method, comprising:
1) carry out FFT conversion to N point Received signal strength, the sequence after conversion is designated as
wherein N is the power of 2;
2) carry out FFT conversion to N point local code signal, get conjugation to the sequence after conversion, the sequence obtained is designated as
3) will
behind ring shift right m position with
corresponding point are multiplied, and obtain sequence
right
do FFT inverse transformation, the sequence of output is designated as
4) will
ring shift left m position, the sequence of output
be the coherent integration results of the different code phase of N kind under m frequency;
5) will
envelope maximal value and thresholding contrast, if exceed thresholding, think that satellite navigation signals to be captured exists, and export corresponding code phase and Doppler frequency.
Described step 1) in,
Wherein,
P
mfor the column permutation matrix of N × N dimension; T represents transposition; M=0,1 ..., N-1; SHIFT
mrepresent the m position that to be moved to left by sequence loops; F
nfor the matrix for calculating Fourier transform,
for baseband signal.
wherein, conj represents and gets conjugation;
for local pseudo-random code vector.
Described step 5) in, thresholding
wherein, σ
nfor noise mean square root, P
fafor single detects, (word detects and refers to: under a certain group of given pseudo-code phase and Doppler frequency, the satellite navigation signals of reception is reappeared signal with this locality and carries out related operation, the envelope value of operation result and the thresholding of setting contrast, if be greater than thresholding, think that satellite navigation signals to be captured exists, if be less than thresholding, think that signal does not exist) false-alarm probability, ln is for getting natural logarithm.
Compared with prior art, the beneficial effect that the present invention has is: method of the present invention is without the need to carrying out to signal wire the search that frequency conversion can complete Doppler frequency, and search precision is high, calculates simple, easily realizes.
Accompanying drawing explanation
Fig. 1 is the schematic diagram of parallel code phase search;
Fig. 2 is embodiment of the present invention process flow diagram.
Embodiment
When catching, frequency search stepping is set to
frequency f to be searched has N number of, can be expressed as:
m=0,1,...,N-1;
Corresponding local carrier signal vector can be expressed as:
According to the Shifting Property of discrete Fourier transform (DFT), can obtain Received signal strength be multiplied with local carrier signal vector after the Fourier transform of sequence be shown below
Wherein, SHIFT
mrepresent the m position that to be moved to left by sequence loops, subscript T represents transposition.F
nfor the matrix for calculating Fourier transform, its expression formula is as follows:
Local pseudo-random code vector
fFT get conjugation after expression formula be:
Wherein, conj represents and gets conjugation.
Will
with
after corresponding point are multiplied, then carrying out inverse Fourier transform, can to obtain the expression formula of coherent integration results as follows:
According to the character of discrete Fourier transformation, above formula can be written as following form:
Wherein, P
mfor the column permutation matrix (being obtained by every a line ring shift left m position of unit matrix) of N × N dimension, it makes every a line ring shift left m time of the vector of premultiplication with it or matrix.Can obtain exempting from mixing according to above formula to walk abreast phase acquisition structure.
As shown in Figure 2, job step of the present invention is as follows:
1) carry out FFT conversion to N point Received signal strength, N is the power of 2, and the sequence after conversion is designated as
2) carry out FFT conversion to N point local code signal, the sequence obtained after getting conjugation to the sequence after conversion is designated as
3) will
behind ring shift right m position with
corresponding point are multiplied, and obtain sequence
right
do FFT inverse transformation, the sequence of output is designated as
4) will
ring shift left m position, the sequence of output
be the coherent integration results of the different code phase of N kind under m frequency.
Calculate by above-mentioned steps
after, will
envelope maximal value and thresholding contrast, if exceed thresholding, think that signal exists, and export corresponding code phase and Doppler frequency.
Claims (5)
1. GNSS signal is exempted from mixing and to be walked abreast a phase acquisition method, it is characterized in that, comprising:
1) carry out FFT conversion to N point Received signal strength, the sequence after conversion is designated as
wherein N is the power of 2;
2) carry out FFT conversion to N point local code signal, get conjugation to the sequence after conversion, the sequence obtained is designated as
3) will
behind ring shift right m position with
corresponding point are multiplied, and obtain sequence
right
do FFT inverse transformation, the sequence of output is designated as
4) will
ring shift left m position, the sequence of output
be the coherent integration results of the different code phase of N kind under m frequency;
5) will
envelope maximal value and thresholding v
tcontrast, if exceed thresholding, thinks that satellite navigation signals to be captured exists, and exports corresponding code phase and Doppler frequency.
2. GNSS signal according to claim 1 is exempted from mixing and to be walked abreast phase acquisition method, it is characterized in that, described step 1) in,
Wherein,
P
mfor the column permutation matrix of N × N dimension; T represents transposition; M=0,1 ..., N-1; SHIFT
mrepresent the m position that to be moved to left by sequence loops; F
nfor the matrix for calculating Fourier transform,
for baseband signal.
3. GNSS signal according to claim 2 is exempted from mixing and to be walked abreast phase acquisition method, it is characterized in that, described step 2) in,
wherein, conj represents and gets conjugation;
for local pseudo-random code vector.
4. GNSS signal according to claim 3 is exempted from mixing and to be walked abreast phase acquisition method, it is characterized in that, described step 4) in,
5. GNSS signal according to claim 4 is exempted from mixing and to be walked abreast phase acquisition method, it is characterized in that, described step 5) in, thresholding
wherein, σ
nfor noise mean square root, P
fafor the false-alarm probability that single detects, ln is for getting natural logarithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510637983.4A CN105182376A (en) | 2015-09-30 | 2015-09-30 | GNSS signal mixer-free parallel code phase capturing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510637983.4A CN105182376A (en) | 2015-09-30 | 2015-09-30 | GNSS signal mixer-free parallel code phase capturing method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105182376A true CN105182376A (en) | 2015-12-23 |
Family
ID=54904572
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510637983.4A Pending CN105182376A (en) | 2015-09-30 | 2015-09-30 | GNSS signal mixer-free parallel code phase capturing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105182376A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106980125A (en) * | 2016-01-18 | 2017-07-25 | 北京信息科技大学 | The drop amount of calculation deblurring processing method that XFAST is captured in a kind of satellite navigation |
CN109581433A (en) * | 2018-11-28 | 2019-04-05 | 西安开阳微电子有限公司 | A kind of L5 signal acquisition methods, device and computer storage medium |
CN110068841A (en) * | 2019-05-06 | 2019-07-30 | 西安开阳微电子有限公司 | A kind of satellite-signal high sensitivity quick capturing method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070206663A1 (en) * | 2006-03-02 | 2007-09-06 | Qualcomm Incorporated | Method And Apparatus For Code Space Search In A Receiver |
JP2008160642A (en) * | 2006-12-26 | 2008-07-10 | Seiko Epson Corp | Circuit and method for controlling correlation operation |
CN101788677A (en) * | 2010-01-14 | 2010-07-28 | 上海交通大学 | Globe navigation satellite system signal acquisition scanning system and acquisition method thereof |
EP2511734A2 (en) * | 2011-04-14 | 2012-10-17 | Thales | Dual-frequency receiver for satellite positioning and associated reception method |
-
2015
- 2015-09-30 CN CN201510637983.4A patent/CN105182376A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070206663A1 (en) * | 2006-03-02 | 2007-09-06 | Qualcomm Incorporated | Method And Apparatus For Code Space Search In A Receiver |
JP2008160642A (en) * | 2006-12-26 | 2008-07-10 | Seiko Epson Corp | Circuit and method for controlling correlation operation |
CN101788677A (en) * | 2010-01-14 | 2010-07-28 | 上海交通大学 | Globe navigation satellite system signal acquisition scanning system and acquisition method thereof |
EP2511734A2 (en) * | 2011-04-14 | 2012-10-17 | Thales | Dual-frequency receiver for satellite positioning and associated reception method |
Non-Patent Citations (1)
Title |
---|
徐晓: "GNSS软件接收机的研究与实现", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106980125A (en) * | 2016-01-18 | 2017-07-25 | 北京信息科技大学 | The drop amount of calculation deblurring processing method that XFAST is captured in a kind of satellite navigation |
CN109581433A (en) * | 2018-11-28 | 2019-04-05 | 西安开阳微电子有限公司 | A kind of L5 signal acquisition methods, device and computer storage medium |
CN109581433B (en) * | 2018-11-28 | 2023-08-18 | 西安开阳微电子有限公司 | L5 signal capturing method and device and computer storage medium |
CN110068841A (en) * | 2019-05-06 | 2019-07-30 | 西安开阳微电子有限公司 | A kind of satellite-signal high sensitivity quick capturing method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103645483B (en) | Beidou signal capturing method in weak signal environment | |
CN105577229B (en) | CPU auxiliary GPU realizes spread-spectrum signal quick capturing method | |
CN100553168C (en) | A kind of frequency folding time frequency parallel searching method of long-period spread spectrum code | |
CN106772471B (en) | A kind of long code segmentation overlay local correlation catching method based on GPU | |
CN101246210B (en) | GPS locating method | |
CN104360357A (en) | Quick Beidou satellite signal capturing method and system based on circulation mode | |
CN105182376A (en) | GNSS signal mixer-free parallel code phase capturing method | |
CN107153206A (en) | A kind of Big Dipper satellite signal quick capturing method based on FFT | |
CN101561484A (en) | Method for acquiring pseudo code of GNSS signal | |
CN103744061A (en) | Iterative least square method-based MIMO (multiple input multiple output) radar DOA (direction-of-arrival) estimation method | |
CN105553506B (en) | A kind of quick capturing method and device of long code spread-spectrum signal | |
CN103675850B (en) | A kind of double-frequency GPS signal acquisition methods based on segmentation split-radix Fourier transform | |
CN105717522A (en) | Second-generation BeiDou B1 frequency band weak signal capturing method | |
CN102353968A (en) | GPS (Global Positioning System) signal acquisition method based on FPGA (Field Programmable Gate Array) and GPS signal acquisition system | |
CN106680842B (en) | A kind of weak signal catching method for satellite navigation system secondary coding | |
CN102279402B (en) | Frequency domain coherent joint capture method of GNSS data/pilot mixed signal | |
CN104459734B (en) | Beidou satellite navigation signal capturing method based on NH code element jumping detection | |
CN104199060A (en) | Satellite navigation signal capturing method and device based on blind aliasing and blind separation | |
CN100553169C (en) | Short-period spreading code time frequency parallel searching method based on frequency folding | |
CN103344968B (en) | INS auxiliary Beidou signal capturing method | |
CN106291610A (en) | A kind of parallel correlation module of compression for GNSS signal compression capture processing means and its implementation | |
CN103197328A (en) | High dynamic satellite receiver fast capture method based on 2 ms data storage | |
CN103926604A (en) | Weak signal capturing method based on overlapping difference cycle coherent integration | |
CN104035109B (en) | Weak signal catching method based on overlapping 1/5 bit difference circulation coherent integration | |
CN101063717B (en) | Method for capturing GPS signal by FFT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151223 |