CN103713299A - Method and device for capturing Galileo signals - Google Patents

Method and device for capturing Galileo signals Download PDF

Info

Publication number
CN103713299A
CN103713299A CN201310746897.8A CN201310746897A CN103713299A CN 103713299 A CN103713299 A CN 103713299A CN 201310746897 A CN201310746897 A CN 201310746897A CN 103713299 A CN103713299 A CN 103713299A
Authority
CN
China
Prior art keywords
road
ranging code
local
pilot tone
galileo
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
CN201310746897.8A
Other languages
Chinese (zh)
Other versions
CN103713299B (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.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN201310746897.8A priority Critical patent/CN103713299B/en
Publication of CN103713299A publication Critical patent/CN103713299A/en
Application granted granted Critical
Publication of CN103713299B publication Critical patent/CN103713299B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30Acquisition or tracking or demodulation of signals transmitted by the system code related
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system

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 method and device for capturing Galileo signals, and belongs to the technical field of satellite positioning. The method includes the steps that product operation is performed on local data circuit ranging codes selected every time in the current selection period and local pilot frequency circuit ranging codes selected every time in the current selection period so that local codes used for operation can be acquired, and square operation is performed on a Galileo E1 signal received every time the local ranging codes are selected in the current selection period so that the square of the Galileo E1 signal can be acquired; the local codes used for operation and the squares of the Galileo E1 signals are input into a multiplying unit to perform multiplication so that the result of each calculation can be obtained, the result of each calculation in the current selection period is input into a coherent integrator, the correlation values between the local data circuit and the Galileo E1 signals as well as between the pilot frequency circuit and the Galileo E1 signals in the current selection period are calculated, and if the correlation values are larger than or equal to a preset threshold value, it is determined that the Galileo E1 signals are successively captured. With the method and device, the operation amount is reduced, and the system expense is saved.

Description

A kind of method that galileo signals is caught and device
Technical field
The present invention relates to satellite positioning tech field, method and device that particularly a kind of galileo signals is caught.
Background technology
Along with popularizing of satellite positioning tech, make people can obtain more convenient, more in time, more accurately desired position information.In position fixing process, what first receiver needed to catch satellite launch passes through CBOC(Composite Binary Offset Carrier, combination binary offset carrier) the Galileo E1 signal that modulation system is modulated, find the code phase of correlation maximum in the ranging code of ranging code in the current Doppler frequency of this Galileo E1 signal and Galileo E1 signal and receiver this locality, just can follow the tracks of and text demodulation realizes positioning function Galileo E1 signal.Wherein, to catching of Galileo E1 signal, need to Doppler and ranging code code phase, search for respectively, find the Doppler frequency and the ranging code code phase that make correlation maximum.In addition, CBOC modulation system can make the autocorrelation function in Galileo E1 signal produce peak, a plurality of limit, can cause receiver acquisition false alarm rate to rise, and causes the error capture to Galileo E1 signal, affects the positioning performance of receiver.
The mode of carrying out Galileo E1 signal capture in prior art can be: by changing local carrier Doppler frequency and local ranging code phase place, Galileo E1 signal is carried out to two-dimensional search, obtain the ranging code that makes in Galileo E1 signal and local ranging code correlation Doppler frequency and local ranging code phase place when maximum.The autocorrelation function solving in prior art in Galileo E1 signal produces peak, a plurality of limit, cause the mode of the problem of receiver error capture to be: by SCPC(Sub Carrier Phase Cancellation, sub-carrier phase method of offset) mode, produce one group and Galileo E1 signal homophase subcarrier He Yi road orthogonal sub-carriers, carry out phase multiplication with Galileo E1 signal respectively, and respectively the data road ranging code of Galileo E1 signal and pilot tone road ranging code are carried out to two-dimensional search.
In realizing process of the present invention, inventor finds that prior art at least exists following problem:
1, prior art need to be carried out two-dimensional search to catching of Galileo E1 signal, travels through the code phase under each Doppler frequency, and search time is longer;
2, in the ranging code of the ranging code of the mode by SCPC in finding Galileo E1 signal and receiver this locality in the process of the code phase of correlation maximum, need to carry out phase multiplication with Galileo E1 signal respectively to two groups of subcarriers, operand is large, the system resource overhead of receiver is large, causes the efficiency of catching Galileo E1 signal of receiver low.
Summary of the invention
In order to solve the problem of prior art, method and device that the embodiment of the present invention provides a kind of galileo signals to catch.Described technical scheme is as follows:
On the one hand, a kind of method that provides galileo signals to catch, described method comprises:
The first step, arranges the access code phase place in the ranging code queue of local data road and the ranging code queue of local pilot tone road;
Second step, from described access code phase place, in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, according to the default time interval, choose one by one local data road ranging code and local pilot tone road ranging code, until the local data road ranging code in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone and local pilot tone road ranging code complete one, choose the cycle;
The 3rd step, current local data road ranging code and the current local pilot tone road ranging code of at every turn choosing in the cycle of choosing at every turn chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and choose in the cycle each result of calculating and input to coherent integration device current, calculate the current correlation of choosing local data road and pilot tone road and Galileo E1 signal in the cycle;
The 4th step, if described correlation is more than or equal to predetermined threshold value, determine the current code phase of the local data road that arranges in the cycle and pilot tone road ranging code of choosing identical with the code phase of data road in described Galileo E1 signal and pilot tone road ranging code, described Galileo E1 signal capture success;
The 5th step, if described correlation is less than predetermined threshold value, changes the access code phase place in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, returns to second step.
Preferably, after the access code phase place of the described local ranging code of described change, described method also comprises:
Numerical value in described coherent integration device is carried out to zero clearing processing.
Preferably, the access code phase place of the described local ranging code of described change, comprising:
A described access code phase delay default value chip to Set For Current in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, and the described local data road ranging code queue of the code phase after postponing in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone after change and the access code phase place of described local pilot tone road ranging code queue.
Preferably, the cycle of choosing described in is data road and the integral multiple in pilot tone road ranging code cycle in described Galileo E1 signal.
On the other hand, the device that provides a kind of galileo signals to catch, described device comprises:
Module is set, for the access code phase place of the ranging code queue of local data road and the ranging code queue of local pilot tone road is set;
Choose module, for from described access code phase place, in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, according to the default time interval, choose one by one local data road ranging code and local pilot tone road ranging code, until the local data road ranging code in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone and local pilot tone road ranging code complete one, choose the cycle;
Computing module, for choose local data road ranging code and the current local pilot tone road ranging code of choosing in the cycle of choosing that the cycle chooses and carry out product calculation acquisition for the local code of computing current at every turn at every turn, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and choose in the cycle each result of calculating and input to coherent integration device current, calculate the current correlation of choosing local data road and pilot tone road and Galileo E1 signal in the cycle;
Trapping module, if be more than or equal to predetermined threshold value for described correlation, determine the current code phase of the local data road that arranges in the cycle and pilot tone road ranging code of choosing identical with the code phase of data road in described Galileo E1 signal and pilot tone road ranging code, described Galileo E1 signal capture success;
Change module, if be less than predetermined threshold value for described correlation, changes the access code phase place in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, returns and chooses module.
Preferably, described device also comprises:
Zero clearing processing module, for carrying out zero clearing processing by the numerical value of described coherent integration device.
Preferably, described change module specifically for:
A described access code phase delay default value chip to Set For Current in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, and the described local data road ranging code queue of the code phase after postponing in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone after change and the access code phase place of described local pilot tone road ranging code queue.
Preferably, the cycle of choosing described in is data road and the integral multiple in pilot tone road ranging code cycle in described Galileo E1 signal.
The beneficial effect that the technical scheme that the embodiment of the present invention provides is brought is:
By current local data road ranging code and the current local pilot tone road ranging code of choosing in the cycle of choosing chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing at every turn at every turn, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and the current result of choosing each calculating in the cycle is inputed to coherent integration device, calculate the current algorithm of choosing the correlation on local data road and pilot tone road and Galileo E1 signal in the cycle, can be by the process of each Doppler frequency of traversal in calculating the process of correlation eliminate during in computing, therefore without carrying out this step, reduce operand, saved system overhead; In addition, in the ranging code queue of local data road and the ranging code queue of local pilot tone road, complete at most one and complete choose circulation and both can capture Galileo E1 signal, further reduced operand, saved system overhead.
Accompanying drawing explanation
In order to be illustrated more clearly in the technical scheme in the embodiment of the present invention, below the accompanying drawing of required use during embodiment is described is briefly described, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skills, do not paying under the prerequisite of creative work, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is the method flow diagram that galileo signals that the embodiment of the present invention one provides is caught;
Fig. 2 is the method flow diagram that galileo signals that the embodiment of the present invention two provides is caught;
Fig. 3 is the apparatus structure schematic diagram that galileo signals that the embodiment of the present invention three provides is caught.
Embodiment
For making the object, technical solutions and advantages of the present invention clearer, below in conjunction with accompanying drawing, embodiment of the present invention is described further in detail.
Embodiment mono-
A kind of method that the embodiment of the present invention provides galileo signals to catch, referring to Fig. 1, method flow comprises:
101: the access code phase place in the ranging code queue of local data road and the ranging code queue of local pilot tone road is set;
102: from described access code phase place, in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, according to the default time interval, choose one by one local data road ranging code and local pilot tone road ranging code, until the local data road ranging code in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone and local pilot tone road ranging code complete one, choose the cycle;
103: current local data road ranging code and the current local pilot tone road ranging code of at every turn choosing in the cycle of choosing at every turn chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and choose in the cycle each result of calculating and input to coherent integration device current, calculate the current correlation of choosing local data road and pilot tone road and Galileo E1 signal in the cycle;
104: if described correlation is more than or equal to predetermined threshold value, determine the current code phase of the local data road that arranges in the cycle and pilot tone road ranging code of choosing identical with the code phase of data road in described Galileo E1 signal and pilot tone road ranging code, described Galileo E1 signal capture success;
105: if described correlation is less than predetermined threshold value, change the access code phase place in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, return to step 102.
The embodiment of the present invention is by current local data road ranging code and the current local pilot tone road ranging code of choosing in the cycle of choosing chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing at every turn at every turn, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and the current result of choosing each calculating in the cycle is inputed to coherent integration device, calculate the current algorithm of choosing the correlation on local data road and pilot tone road and Galileo E1 signal in the cycle, can be by the process of each Doppler frequency of traversal in calculating the process of correlation eliminate during in computing, therefore without carrying out this step, reduce operand, saved system overhead; In addition, in the ranging code queue of local data road and the ranging code queue of local pilot tone road, complete at most one and complete choose circulation and both can capture Galileo E1 signal, further reduced operand, saved system overhead.
Embodiment bis-
A kind of method that the embodiment of the present invention provides galileo signals to catch, referring to Fig. 2, method flow comprises:
201: the access code phase place in the ranging code queue of local data road and the ranging code queue of local pilot tone road is set.
Wherein, the local data road of setting is identical with the access code phase place of local pilot tone road ranging code.
Further, the access code phase place of choosing first can be chosen a position at random in the ranging code queue of local data road and the ranging code queue of local pilot tone road, can start the first position from the ranging code queue of local data road and the ranging code queue of local pilot tone road, the initial position of access code phase place does not limit at this yet.
202: from described access code phase place, in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, according to the default time interval, choose one by one local data road ranging code and local pilot tone road ranging code, until the local data road ranging code in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone and local pilot tone road ranging code complete one, choose the cycle.
Satellite is when generating Galileo E1 signal, the time interval that changes ranging code corresponding to access code phase place in the ranging code queue on data road and ranging code queue presets, therefore in local receiver can by this default time interval one by one from access code phase place, in the ranging code queue of local data road and the ranging code queue of local pilot tone road, choose local data road ranging code and local pilot tone road ranging code.
Wherein, in the ranging code queue of local data road and the ranging code queue of local pilot tone road, preserved respectively a same code phase ranging code, this ranging code loops to be chosen.For example: in the ranging code queue of local data road, preserve 4092 ranging codes, the access code phase place of supposing Set For Current is 1, from the first place of this local data road ranging code queue, choose, every Preset Time interval, choose successively, choose successively 2,3 ... to 4092, after next Preset Time interval, the first place from the ranging code queue of local data road is that code phase 1 starts again to choose.Wherein, by access code phase place 1, started until 4092, in the ranging code queue of local data road, to have completed once circulation and chosen, so this circulation is chosen process once and is called one and chooses the cycle.
Preferably, at some in particular cases, for example: due to weather problems, or receiver is when the medium non-spacious location of building, receiver receives the poor effect of Galileo E1 signal, therefore one need to be chosen in the cycle and starts according to access code phase place in current period, in the ranging code queue of local data road and the ranging code queue of local pilot tone road, carry out two or more complete cycle, the current cycle of choosing correspondence is chosen data road and the pilot tone road ranging code cycle in two or more Galileo E1 signals.The cycle of choosing is data road and the integral multiple in pilot tone road ranging code cycle in Galileo E1 signal.
203: current local data road ranging code and the current local pilot tone road ranging code of at every turn choosing in the cycle of choosing at every turn chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and choose in the cycle each result of calculating and input to coherent integration device current, calculate the current correlation of choosing local data road and pilot tone road and Galileo E1 signal in the cycle.
Current, choose in the cycle, in the Galileo E1 signal receiving while at every turn choosing local ranging code, the satellite ranging code using while mediating is also chosen again, therefore by the local ranging code of choosing and the Galileo E1 signal receiving at that time at every turn, calculate, finally calculate the correlation that each chooses local ranging code corresponding to cycle and Galileo E1 signal.
Wherein, can being expressed as of Galileo E1 signal:
S ( t ) = 1 2 ( D ( t ) × C E 1 - B ( t ) × Sub E 1 - B ( t ) - C sec ( t ) × C E 1 - C ( t ) × Sub E 1 - C ( t ) ) × cos ( 2 π f E 1 t )
Wherein, D (t) is the navigation message signals on data road, C e1-B(t) be the ranging code on data road, Sub e1-B(t) be the subcarrier on data road, C sec(t) be the secondary code on pilot tone road, C e1-Cthe ranging code on pilot tone road, Sub e1-C(t) be the subcarrier on pilot tone road, f e1for Galileo E1 signal imitation carrier frequency, containing Doppler frequency.
Galileo E1 signal is carried out to square operation:
S square ( t ) = 1 2 ( D ( t ) × C E 1 - B ( t ) × Sub E 1 - B ( t ) - C sec ( t ) × C E 1 - C ( t ) × Sub E 1 - C ( t ) ) × ( D ( t ) × C E 1 - B ( t ) × Sub E 1 - B ( t ) - C sec ( t ) × C E 1 - C ( t ) × Sub E 1 - C ( t ) ) × 1 2 ( 1 + cos ( 4 π f E 1 t ) ) = ( 1 2 ( ( C E 1 - B ( t ) ) 2 × ( Sub E 1 - B ( t ) ) 2 ) + 1 2 ( ( C E 1 - C ( t ) ) 2 × ( Sub E 1 - C ( t ) ) 2 ) - ( D ( t ) × C sec ( t ) × C E 1 - B ( t ) × C E 1 - C ( t ) × V cons ) ) × 1 2 ( 1 + cos ( 4 π f E 1 t ) )
Wherein, Sub e1-Band Sub (t) e1-C(t) amass as certain value V cons.While D (t) and C sec(t) value is only+1 and-1.
Wherein, currently choose the local data road ranging code at every turn chosen in the cycle and the current local pilot tone road ranging code of at every turn choosing in the cycle of choosing carries out product calculation acquisition and for the computation process of the local code of computing is:
C E1-B×C E1-C
By current, choose in the cycle, while choosing local code, carry out above-mentioned two formula calculates at every turn, and the value obtaining is inputed to multiplier, carry out multiplying and obtain operation result, and the result of each computing is inputed to coherent integration device carry out coherent accumulation computing, until one obtains the current local ranging code in cycle and the correlation of Galileo E1 signal chosen while choosing end cycle.
Wherein, if when the code phase of the local code phase of local ranging code of choosing and the ranging code of satellite is identical, containing cos (4 π f e1t) after every integration of high frequency content, result is 0, remains in every and contains C e1-B* C e1-Cits integral result of item of the factor is also 0, and the operation result after input coherent integration device is only:
( D ( t ) × C sec ( t ) × C E 1 - B ( t ) × C E 1 - C ( t ) × V cons ) × ( C E 1 - B ( t ) × C E 1 - C ( t ) ) = D ( t ) × C sec ( t ) × C E 1 - B 2 ( t ) × C E 1 - C 2 ( t ) × V cons
Wherein,
Figure BDA0000450263400000082
with
Figure BDA0000450263400000083
two perseverances are 1, and above formula is converted into:
D(t)×C sec(t)×V cons
Final coherent integration result is only that a constant and navigation message and secondary code are long-pending, has completed the work that carrier wave and Doppler peel off simultaneously.And final correlation can be according to being worth greatly after each operation result coherent accumulation calculating.Therefore, the predetermined threshold value of setting can, for a larger value, judge final correlation.
When if the code phase of local ranging code of choosing local is not identical with the code phase of the ranging code of satellite, all contain C in every in coherent integration process e1-B* C e1-Cthe factor, its integral result is 0, therefore final correlation is very little value or 0.
204: if described correlation is more than or equal to predetermined threshold value, determine the current code phase of the local data road that arranges in the cycle and pilot tone road ranging code of choosing identical with the code phase of data road in described Galileo E1 signal and pilot tone road ranging code, described Galileo E1 signal capture success.
205: if described correlation is less than predetermined threshold value, change the access code phase place in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, return to step 202.
If correlation is less than predetermined threshold value the current code phase that arranges in access code phase place that the cycle arranges and satellite asynchronous of choosing is described, need to re-start and choose.Owing to having completed one in the ranging code queue of local data road and the ranging code queue of local pilot tone road, choose circulation, therefore reset access code phase place, and re-start the process of step 201-step 203.Carry out the next one and choose the acquisition procedure in cycle.
Wherein, in step 205, the process of changing the access code phase place of local ranging code is:
A described access code phase delay default value chip to Set For Current in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, and the described local data road ranging code queue of the code phase after postponing in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone after change and the access code phase place of described local pilot tone road ranging code queue.
206: the numerical value in coherent integration device is carried out to zero clearing processing.
This step is the process of carrying out after step 205 is carried out.For the next one, choose the cycle and carry out the preparation that computing is done, at the next one, choose in the cycle, coherent integration device needs independent computing and does not participate in an operation result of choosing in the cycle.
In embodiments of the present invention, by changing the access code phase place of local ranging code, local data road ranging code in local ranging code and local pilot tone road ranging code are arranged ranging code corresponding to each code phase in the ranging code queue of local data road and the ranging code queue of local pilot tone road at most once as access code phase place, and complete corresponding choosing with computing of the cycle of choosing and just can correctly capture Galileo E1 signal.
The embodiment of the present invention is by current local data road ranging code and the current local pilot tone road ranging code of choosing in the cycle of choosing chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing at every turn at every turn, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and the current result of choosing each calculating in the cycle is inputed to coherent integration device, calculate the current algorithm of choosing the correlation on local data road and pilot tone road and Galileo E1 signal in the cycle, can be by the process of each Doppler frequency of traversal in calculating the process of correlation eliminate during in computing, therefore without carrying out this step, reduce operand, saved system overhead; In addition, in the ranging code queue of local data road and the ranging code queue of local pilot tone road, complete at most one and complete choose circulation and both can capture Galileo E1 signal, further reduced operand, saved system overhead.
Embodiment tri-
The device that the embodiment of the present invention provides a kind of galileo signals to catch, referring to Fig. 3, this device comprises:
Module 301 is set, for the access code phase place of the ranging code queue of local data road and the ranging code queue of local pilot tone road is set;
Choose module 302, for from described access code phase place, in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, according to the default time interval, choose one by one local data road ranging code and local pilot tone road ranging code, until the local data road ranging code in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone and local pilot tone road ranging code complete one, choose the cycle;
Computing module 303, for choose local data road ranging code and the current local pilot tone road ranging code of choosing in the cycle of choosing that the cycle chooses and carry out product calculation acquisition for the local code of computing current at every turn at every turn, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and choose in the cycle each result of calculating and input to coherent integration device current, calculate the current correlation of choosing local data road and pilot tone road and Galileo E1 signal in the cycle;
Trapping module 304, if be more than or equal to predetermined threshold value for described correlation, determine the current code phase of the local data road that arranges in the cycle and pilot tone road ranging code of choosing identical with the code phase of data road in described Galileo E1 signal and pilot tone road ranging code, described Galileo E1 signal capture success;
Change module 305, if be less than predetermined threshold value for described correlation, changes the access code phase place in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, returns and chooses module 302.
Wherein, described device also comprises:
Zero clearing processing module 306, for carrying out zero clearing processing by the numerical value of described coherent integration device.
Wherein, described change module 305 specifically for:
A described access code phase delay default value chip to Set For Current in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, and the described local data road ranging code queue of the code phase after postponing in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone after change and the access code phase place of described local pilot tone road ranging code queue.
Wherein, the cycle of choosing described in is data road and the integral multiple in pilot tone road ranging code cycle in described Galileo E1 signal.
The embodiment of the present invention is by current local data road ranging code and the current local pilot tone road ranging code of choosing in the cycle of choosing chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing at every turn at every turn, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and the current result of choosing each calculating in the cycle is inputed to coherent integration device, calculate the current algorithm of choosing the correlation on local data road and pilot tone road and Galileo E1 signal in the cycle, can be by the process of each Doppler frequency of traversal in calculating the process of correlation eliminate during in computing, therefore without carrying out this step, reduce operand, saved system overhead; In addition, in the ranging code queue of local data road and the ranging code queue of local pilot tone road, complete at most one and complete choose circulation and both can capture Galileo E1 signal, further reduced operand, saved system overhead.
The invention described above embodiment sequence number, just to describing, does not represent the quality of embodiment.
One of ordinary skill in the art will appreciate that all or part of step that realizes above-described embodiment can complete by hardware, also can come the hardware that instruction is relevant to complete by program, described program can be stored in a kind of computer-readable recording medium, the above-mentioned storage medium of mentioning can be ROM (read-only memory), disk or CD etc.
The foregoing is only preferred embodiment of the present invention, in order to limit the present invention, within the spirit and principles in the present invention not all, any modification of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (8)

1. the method that galileo signals is caught, is characterized in that, described method comprises:
The first step, arranges the access code phase place in the ranging code queue of local data road and the ranging code queue of local pilot tone road;
Second step, from described access code phase place, in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, according to the default time interval, choose one by one local data road ranging code and local pilot tone road ranging code, until the local data road ranging code in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone and local pilot tone road ranging code complete one, choose the cycle;
The 3rd step, current local data road ranging code and the current local pilot tone road ranging code of at every turn choosing in the cycle of choosing at every turn chosen in the cycle chosen carried out to product calculation acquisition for the local code of computing, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and choose in the cycle each result of calculating and input to coherent integration device current, calculate the current correlation of choosing local data road and pilot tone road and Galileo E1 signal in the cycle;
The 4th step, if described correlation is more than or equal to predetermined threshold value, determine the current code phase of the local data road that arranges in the cycle and pilot tone road ranging code of choosing identical with the code phase of data road in described Galileo E1 signal and pilot tone road ranging code, described Galileo E1 signal capture success;
The 5th step, if described correlation is less than predetermined threshold value, changes the access code phase place in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, returns to second step.
2. method according to claim 1, is characterized in that, after the access code phase place of the described local ranging code of described change, described method also comprises:
Numerical value in described coherent integration device is carried out to zero clearing processing.
3. method according to claim 1, is characterized in that, the access code phase place of the described local ranging code of described change, comprising:
A described access code phase delay default value chip to Set For Current in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, and the described local data road ranging code queue of the code phase after postponing in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone after change and the access code phase place of described local pilot tone road ranging code queue.
4. method according to claim 1, is characterized in that, described in to choose the cycle be data road and the integral multiple in pilot tone road ranging code cycle in described Galileo E1 signal.
5. the device that galileo signals is caught, is characterized in that, described device comprises:
Module is set, for the access code phase place of the ranging code queue of local data road and the ranging code queue of local pilot tone road is set;
Choose module, for from described access code phase place, in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, according to the default time interval, choose one by one local data road ranging code and local pilot tone road ranging code, until the local data road ranging code in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone and local pilot tone road ranging code complete one, choose the cycle;
Computing module, for choose local data road ranging code and the current local pilot tone road ranging code of choosing in the cycle of choosing that the cycle chooses and carry out product calculation acquisition for the local code of computing current at every turn at every turn, and by current choose the Galileo E1 signal receiving while at every turn choosing local ranging code in the cycle carry out square operation obtain described Galileo E1 signal square; Square input multiplier of the described local code for computing and described Galileo E1 signal is carried out to multiplying and obtain each result of calculating, and choose in the cycle each result of calculating and input to coherent integration device current, calculate the current correlation of choosing local data road and pilot tone road and Galileo E1 signal in the cycle;
Trapping module, if be more than or equal to predetermined threshold value for described correlation, determine the current code phase of the local data road that arranges in the cycle and pilot tone road ranging code of choosing identical with the code phase of data road in described Galileo E1 signal and pilot tone road ranging code, described Galileo E1 signal capture success;
Change module, if be less than predetermined threshold value for described correlation, changes the access code phase place in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, returns and chooses module.
6. device according to claim 5, is characterized in that, described device also comprises:
Zero clearing processing module, for carrying out zero clearing processing by the numerical value of described coherent integration device.
7. device according to claim 5, is characterized in that, described change module specifically for:
A described access code phase delay default value chip to Set For Current in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone, and the described local data road ranging code queue of the code phase after postponing in the ranging code queue of described local data road and the road ranging code queue of described local pilot tone after change and the access code phase place of described local pilot tone road ranging code queue.
8. device according to claim 5, is characterized in that, described in to choose the cycle be data road and the integral multiple in pilot tone road ranging code cycle in described Galileo E1 signal.
CN201310746897.8A 2013-12-30 2013-12-30 A kind of method that galileo signals is caught and device Active CN103713299B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310746897.8A CN103713299B (en) 2013-12-30 2013-12-30 A kind of method that galileo signals is caught and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310746897.8A CN103713299B (en) 2013-12-30 2013-12-30 A kind of method that galileo signals is caught and device

Publications (2)

Publication Number Publication Date
CN103713299A true CN103713299A (en) 2014-04-09
CN103713299B CN103713299B (en) 2016-03-02

Family

ID=50406424

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310746897.8A Active CN103713299B (en) 2013-12-30 2013-12-30 A kind of method that galileo signals is caught and device

Country Status (1)

Country Link
CN (1) CN103713299B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105842714A (en) * 2016-06-07 2016-08-10 北京邮电大学 Satellite signal acquisition method and device
CN112698364A (en) * 2020-12-10 2021-04-23 成都国星通信有限公司 Compatible modern GNSS signal ranging code generation method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101262246A (en) * 2007-03-07 2008-09-10 联发科技股份有限公司 Method and apparatus for spreading signal acquisition and tracking
US20080288168A1 (en) * 2005-11-02 2008-11-20 Thales Method for Determining the Position of Satellites in a Navigation System
CN102411148A (en) * 2011-10-10 2012-04-11 北京大学 Method for rapidly capturing navigation satellite signal and apparatus thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080288168A1 (en) * 2005-11-02 2008-11-20 Thales Method for Determining the Position of Satellites in a Navigation System
CN101262246A (en) * 2007-03-07 2008-09-10 联发科技股份有限公司 Method and apparatus for spreading signal acquisition and tracking
CN102411148A (en) * 2011-10-10 2012-04-11 北京大学 Method for rapidly capturing navigation satellite signal and apparatus thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
祝雪芬等: "伽利略E1-BOC信号的改进捕获算法及性能分析", 《东南大学学报(自然科学版)》, vol. 39, no. 5, 30 September 2009 (2009-09-30), pages 941 - 945 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105842714A (en) * 2016-06-07 2016-08-10 北京邮电大学 Satellite signal acquisition method and device
CN105842714B (en) * 2016-06-07 2018-01-05 北京邮电大学 Satellite signal tracking method and device
CN112698364A (en) * 2020-12-10 2021-04-23 成都国星通信有限公司 Compatible modern GNSS signal ranging code generation method
CN112698364B (en) * 2020-12-10 2023-11-14 成都国星通信有限公司 Compatible modern GNSS signal ranging code generation method

Also Published As

Publication number Publication date
CN103713299B (en) 2016-03-02

Similar Documents

Publication Publication Date Title
CN102662183B (en) Method and system for global position system (GPS) signal capture
US10027371B2 (en) Spreading signal generating method, generating device, receiving method and receiving device
CN104570017B (en) Avoiding and stripping method for Neumann-Hoffman codes in navigation messages of Beidou navigation satellite system D1
CN104898136A (en) Capturing method of Beidou second-generation B1 signals and system thereof
CN103926604B (en) Based on the weak signal catching method of overlapping difference circulation coherent integration
CN102508269B (en) Satellite navigation pilot signal acquisition method, pseudo random sequence stripping method and device
CN104181556A (en) BOC modulating signal capturing method based on overlapped difference circulation coherent integration
CN101977091B (en) Method and system for monitoring electromagnetic spectrum
CN103926605B (en) GPS weak signal catching method based on difference circulation coherent integration
CN107888230B (en) Multiphase coherent integration acquisition method applied to weak signal environment
CN111147102B (en) Low signal-to-noise ratio code capturing method based on frequency modulation signal
CN106526632A (en) GNSS capture correlator and signal capture method
CN102355279B (en) Method and system for diversity maximum likelihood spread spectrum communication bit synchronization
CN104035109A (en) Weak signal capturing method based on 1/5 bit differential circulation coherent integration
CN103439718A (en) Unambiguous tracking unit of high-order BOC modulation signals
CN103713299B (en) A kind of method that galileo signals is caught and device
US8665997B1 (en) AIS demodulator for unknown carrier phase and baud alignment
CN113885057A (en) Satellite internet of things signal processing method based on coherent-incoherent spread spectrum modulation
CN103969664A (en) GPS weak signal capturing method based on differential circular coherent integration
CN109031360B (en) Method and device for capturing weak GPS signal
CN109633709A (en) Practical and efficient frequency deviation estimating method in a kind of satellite communication system
CN108957492B (en) L1C/A and L1C combined capturing method of GPS
CN103616699A (en) Binary coded symbol (BCS) optimization and modulation method based on minimum shift keying (MSK) pulse
CN103969668B (en) It is a kind of to be based on frequency compensated GPS bit synchronization methods
CN102340325B (en) Time-domain capturing method for secondary direct sequence spread spectrum signal

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant