CN106253946B - A kind of satellite navigation signals data/pilot branch spread spectrum code optimization matching method - Google Patents
A kind of satellite navigation signals data/pilot branch spread spectrum code optimization matching method Download PDFInfo
- Publication number
- CN106253946B CN106253946B CN201610654757.1A CN201610654757A CN106253946B CN 106253946 B CN106253946 B CN 106253946B CN 201610654757 A CN201610654757 A CN 201610654757A CN 106253946 B CN106253946 B CN 106253946B
- Authority
- CN
- China
- Prior art keywords
- code
- pilot
- data
- tracking
- branch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
- H04B1/7075—Synchronisation aspects with code phase acquisition
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
- H04B1/7085—Synchronisation aspects using a code tracking loop, e.g. a delay-locked loop
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7097—Interference-related aspects
- H04B1/711—Interference-related aspects the interference being multi-path interference
- H04B1/7115—Constructive combining of multi-path signals, i.e. RAKE receivers
- H04B1/7117—Selection, re-selection, allocation or re-allocation of paths to fingers, e.g. timing offset control of allocated fingers
Abstract
The invention discloses a kind of satellite navigation signals data/pilot branches to spread code optimization matching method, by being redistributed to data branch spreading code and pilot tone branch spreading code, ensure under certain related spacing (subtracting space lag in advance), code tracking mean error minimum caused by cross-correlation interference in channel when independent TPS Tracking Pilot Signal.The specific steps are：First, data and the pilot spreading code race of signal to be optimized are defined；Then, code tracking error caused by cross-correlation interference in the channel under any one data/pilot spreading code pairing relationship is calculated；Then, code tracking mean error under different data/pilot spreading code pairing relationship is obtained according to code tracking error matrix；Finally, it by minimizing code tracking mean error, obtains data/pilot and spreads code optimization matching side case.Present invention reduces the cross correlations between data branch and pilot tone branch, improve the tracking performance of satellite navigation signals.
Description
Technical field
The invention belongs to GNSS signal system design technical fields, and in particular to a kind of satellite navigation signals data/pilot
Branch spreads code optimization matching method.
Background technology
Pseudo noise code (or spreading code) is not only the foundation that receiver distinguishes different satellite-signals, is also to discriminate between same signal
The foundation of different components (such as data and pilot tone).Especially for the pseudo noise code in data and pilot channel, since they are same
When (code phase is consistent), with phase (carrier phase is identical) the characteristics of, the size of cross correlation will directly affect GNSS and connect
The code tracking performance of receipts machine.With the propulsion of GPS modernization and Galileo system Constructions, the research work of pseudo-code design aspect
Also it has made great progress.To improve signal performance, Galileo E1 OS introduce completely new Random codes.With traditional puppet
Random code is different, and Random codes cannot be generated by shift register.Meanwhile GPS L1C also select new pseudo-code sequence, i.e. base
In the Weil codes of Legendre sequence.The appearance of Weil codes and Random codes is related with the new proposition of pseudo-code criterion, such as balance, zero
Autocorrelation sidelobe characteristic, even odd correlation etc..The pseudo noise code of GPS L1C and Galileo E1 OS is mainly examined in the design phase
The characteristic (such as even odd correlation) for considering code book body the spies such as coexists for specific modulation system, multiplex mode, data/pilot
Point considers less.Theoretically, Random codes and Weil codes are only only when data and pilot signal use orthogonal modulation optimal
's.But in fact, either Galileo E1 OS or GPS L1C, their data are all with mutually transmitting with pilot signal
, this largely exacerbates the cross correlation between data and pilot tone branch.Cross correlation between data and pilot tone
The tracking performance of satellite navigation signals can be influenced, and then influences the range accuracy and positioning accuracy of user, it is therefore necessary to be studied
Influence and pairing optimization method of the data/pilot signal cross correlation to navigation signal performance.
Invention content
The present invention proposes satellite navigation signals data/pilot branch spread spectrum code optimization matching method, by data and
Pilot spreading code is redistributed, and ensures that under the conditions of certain related spacing (subtracting space lag in advance), individually tracking is led
Code tracking mean error caused by cross-correlation interference is minimum in channel when frequency signal.
The present invention is achieved through the following technical solutions：A kind of satellite navigation signals data/pilot branch spread spectrum code optimization
Matching method includes the following steps：
(1) under default related spacing, code tracking error caused by all pairing relationships is calculated, code tracking error moments are obtained
Battle array；The pairing relationship be data branch spreading code and the spreading code of pilot tone branch between one-to-one relationship, code with
Track error is determined by individually tracking pilot tone tributary signal time-code phase demodulation curve offset；
(2) the code tracking error matrix for utilizing step (1), calculates the code tracking mean error under different pairing relationships；
(3) match the minimum corresponding pairing relationship of code tracking mean error as data/pilot branch spread spectrum code optimization
To scheme.
Wherein, the related spacing is to subtract space lag in advance.
Wherein, the result of calculation of the code phase demodulation curve offset described in step (1) is related with code phase phase discriminator type, code phase
Position phase discriminator is relevant advanced to subtract retarded type, incoherent advanced subtract lag encapsulating type, incoherent subtract after-power type, quasi- phase in advance
Do long-pending power-type or relevant dot product type.
Wherein, using different related spacing or different types of code phase phase discriminator, different optimization pairings can be obtained
Scheme.
The present invention is compared to the advantages of background technology：
This method is ensured by being redistributed to data branch spreading code and pilot tone branch spreading code certain related
Under spacing condition, code tracking mean error caused by cross-correlation interference is minimum in channel when independent TPS Tracking Pilot Signal.It reduces
Cross correlation between data branch and pilot tone branch, improves the tracking performance of satellite navigation signals, and then improves user
Range accuracy and positioning accuracy.
Description of the drawings
Fig. 1 is the method for the present invention flow chart of steps；
Fig. 2 is the S- curve offset errors of Galileo E1 OS signals before optimization；
Fig. 3 is Galileo E1 OS signal datas/pilot spreading code pairing relationship after optimization.
Fig. 4 is the front and back Galileo E1 OS signal datas/pilot spreading code pairing relationship comparison diagram of optimization.
Fig. 5 is the front and back Galileo E1 OS signal code tracking mean error comparison diagrams of optimization.
Specific implementation mode
It elaborates to the present invention below in conjunction with drawings and examples.The scope of the present invention is not by these embodiments
Limitation, the scope of the present invention propose in detail in the claims.
The spread spectrum race of the spread spectrum race and pilot tone branch of given data branch, is denoted as a respectivelydAnd ap, it is defined as：
Wherein, N represents the spread spectrum code sequence number that data branch and pilot tone branch include,Represent kth
The spread spectrum code sequence of a data branch,The spread spectrum code sequence of k-th of pilot tone branch is represented, T represents transposition.
The present invention keeps data branch spreading code, the grouping of pilot tone branch spreading code constant, only changes the spread spectrum of pilot channel
The code relations of distribution.That is, the front and back data branch spreading code race of optimization and pilot tone branch spreading code race remain unchanged, only it is pair
It should be related to and be changed.After pairing relationship optimization, when ensureing independent TPS Tracking Pilot Signal, in channel caused by cross-correlation interference
Code tracking mean error is minimum.
For the data/pilot spreading code pairing relationship optimized, when needing to calculate independent TPS Tracking Pilot Signal, Suo Youke
Code phase demodulation curve (i.e. S- curves) offset under the spreading code pairing relationship of energy.Code tracking defines error matrix BN×N：
Wherein,Represent spreading code pairing (i.e.As data branch spreading code,As leading
Frequency branch spreading code),S- curves caused by cross-correlation interference are inclined in channel when representing independent TPS Tracking Pilot Signal
It moves.
The code phase delay of receiver S- curve zero crossings, as S- curve offsets.The concrete outcome of S- curve offsets with
It is related to receive the code phase phase discriminator type used, specific code phase phase discriminator includes the relevant advanced retarded type, incoherent of subtracting
It is advanced to subtract lag encapsulating type, incoherent advanced subtract after-power type, accurate relevant dot product power-type, relevant dot product type etc..Below with
It is concerned with and illustrates the circular of S- curve offsets for subtracting retarded type in advance.Phase discriminator is lagged for relevant advanced subtract
It says, S- curves D (ε) is represented by：
Wherein,It is to receive signal spread-spectrum code c1With receiver local spreading code c2Cross-correlation function.It is corresponding
Zero crossing (i.e. keyed end) is defined as：
D(εbias)=0 (4)
Wherein, Δ is to subtract lag spacing (i.e. related spacing) in advance, and ε is code phase delay, εbiasIt is exactly S- curve offsets.
Code tracking error matrix BN×NIt is independent tracking pilot tone branch as a result, thereforeIt is that receiver local spreading code isWhen εbiasSpecial case.
The best pairing mode of data/pilot is selected to ensure that code tracking mean error is minimum, can be converted into and minimize generation
Valence function：
Wherein
Select the restrictive condition of pairing relationship as follows：
(1)That is code tracking error matrix BN×NEach row in one and only one element energy
Enough it is selected；
(2)That is code tracking error matrix BN×NEvery a line in one and only one element energy
Enough it is selected.
In simple terms, which is actually in code tracking error matrix BN×NOften row and each column in choose and
An element (shared N number of) is only selected, ensures the sum of their absolute value minimum.
It solves after minimizing cost function z, according to vij, i, j=1,2 ... the value of N can obtain data/pilot branch
Result is matched in the optimization of road spreading code.Even vij=1, then it representsIt is selected one group of data/pilot branch spreading code
It is right, otherwiseJust not in optimum results.
Below the data/pilot spreading code of the E1 OS signals of Galileo navigation system match using the present invention excellent
Change.According to downlink signal Interface Control File, Galileo E1 OS include 50 data and pilot spreading code sequence.It is specific real
Apply that steps are as follows：
The first step defines Galileo E1 OS data spreading codes race adWith pilot spreading code race ap：
Second step when independent TPS Tracking Pilot Signal, calculates the S- curve offsets of all possible spreading code pair, obtain code with
Track error matrix：
Assuming that receiver select it is relevant subtract lag phase discriminator in advance, and subtract the chip of lag separation delta=1 in advance, then specifically
Code tracking error matrix it is as shown in Figure 2.
Third walks, and minimizes cost function z.The most intuitive solution of the optimization problem is exactly enumerative technique.But due to N values
Larger, the operand of enumerative technique is excessive, should not use.In fact, pairing optimization here is considered as the warp in " operational research "
Allusion quotation problem --- assignment problem (Assignment Problem).Assignment problem can pass through " Hungary " algorithm Efficient Solution.
The solution procedure of " Hungary " algorithm is as follows：
Process one：Find out absolute code tracking error matrix | BN×N| the least member in per a line；Where least member
Row subtracts the least member from each element, builds a new matrix；For this new matrix, find out in each row
Least member.In the row where least member, the least member is subtracted from each element, the new matrix of structure one (claims
To reduce error matrix)；
Process two：Draw the straight line (horizontal, vertical or both) of minimum number so that they can cover reduction error moments
All 0 in battle array；If necessary to N straight line, then optimal solution can be found out in capped 0 of matrix；If needed
The straight line wanted is less than N items, then is transferred to process three；
Process three：In reducing error matrix, find out not by the non-zero element of minimum of drawn straight line covering in process two
(its value is denoted as ζ)；It is not covered from each of reduction error matrix now and subtracts ζ in element, to the member covered by two straight lines
Element adds ζ；It is then back to process two.
4th step, according to vij, i, j=1,2 ... the value of N can obtain the optimization pairing of data/pilot branch spreading code
As a result.Even vij=1, then it representsIt is selected one group of data/pilot branch spreading code pair, otherwiseJust not
In optimum results.
Galileo E1 OS data/pilot spreading code pairing relationships after optimization are as shown in Figure 3.Wherein, " satellite number/number
According to " row provides is satellite number in Galileo E1 OS Interface Control File and data spreading code number.Due to not changing
Data spreading code race, therefore keep satellite number, data spreading code number consistent in Interface Control File here." pilot tone " row
What is provided is the pilot spreading code number in Galileo E1 OS Interface Control File.The spreading code and its lastrow same column
Data spreading code forms new spreading code pair.For example, the data/pilot spreading code pairing relationship of satellite 1 is (Isosorbide-5-Nitrae 7) after optimization,
Wherein, " 1 " represents data spreading code sequence number in Interface Control File, and " 47 " represent pilot tone in Interface Control File and spread
Code sequence number.Obviously, optimize the spreading code of front and back PRN 6,33 and 40 to not changing.The front and back Galileo E1 of optimization
Fig. 4 is shown in the comparison of OS data/pilot spreading code pairing relationships.The front and back Galileo E1 OS code tracking mean errors comparison of optimization is shown in
Fig. 5.Obviously, after the optimization of this patent method, the code tracking performance of Galileo E1 OS will be obviously improved.
Claims (4)
1. a kind of satellite navigation signals data/pilot branch spreads code optimization matching method, which is characterized in that including following step
Suddenly：
（1）Under default related spacing, code tracking error caused by all pairing relationships is calculated, code tracking error matrix is obtained；
The pairing relationship is the one-to-one relationship between the spreading code and the spreading code of pilot tone branch of data branch, and code tracking misses
Difference is determined by individually tracking pilot tone tributary signal time-code phase demodulation curve offset；
（2）Utilize step（1）Code tracking error matrix, calculate the code tracking mean error under different pairing relationships；
（3）Minimum code tracking mean error is obtained using cost function is minimized, minimum code tracking mean error is corresponded to
Pairing relationship as data/pilot branch spread code optimization matching side case.
2. a kind of satellite navigation signals data/pilot branch according to claim 1 spreads code optimization matching method, special
Sign is：The related spacing is to subtract space lag in advance.
3. a kind of satellite navigation signals data/pilot branch according to claim 1 spreads code optimization matching method, special
Sign is：Step（1）The result of calculation of the code phase demodulation curve offset is related with code phase phase discriminator type, code phase phase demodulation
Device is relevant advanced to subtract retarded type, incoherent advanced subtract lag encapsulating type, incoherent advanced subtract after-power type, accurate relevant dot product
Power-type or relevant dot product type.
4. a kind of satellite navigation signals data/pilot branch according to claim 3 spreads code optimization matching method, special
Sign is：Using different related spacing or different types of code phase phase discriminator, different optimization pairing schemes can be obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610654757.1A CN106253946B (en) | 2016-08-11 | 2016-08-11 | A kind of satellite navigation signals data/pilot branch spread spectrum code optimization matching method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610654757.1A CN106253946B (en) | 2016-08-11 | 2016-08-11 | A kind of satellite navigation signals data/pilot branch spread spectrum code optimization matching method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106253946A CN106253946A (en) | 2016-12-21 |
| CN106253946B true CN106253946B (en) | 2018-07-24 |
Family
ID=58078251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610654757.1A Active CN106253946B (en) | 2016-08-11 | 2016-08-11 | A kind of satellite navigation signals data/pilot branch spread spectrum code optimization matching method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106253946B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109495130A (en) * | 2017-09-10 | 2019-03-19 | 北京遥感设备研究所 | A kind of missile-borne band spread receiver time-division single channel code phase discrimination method |
| CN110045398B (en) * | 2019-04-30 | 2020-11-10 | 中国人民解放军国防科技大学 | Optimal correlation interval-based code phase zero crossing point deviation suppression method |
| CN113691344B (en) * | 2021-10-25 | 2021-12-28 | 长沙海格北斗信息技术有限公司 | Multi-path Weil code generation method and device |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102736520A (en) * | 2011-04-13 | 2012-10-17 | 北京邮电大学 | Satellite navigation system principle simulation method and satellite signal simulator |
| EP2745141A1 (en) * | 2011-08-16 | 2014-06-25 | European Space Agency | A navigation system using spreading codes based on pseudo-random noise sequences |
-
2016
- 2016-08-11 CN CN201610654757.1A patent/CN106253946B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102736520A (en) * | 2011-04-13 | 2012-10-17 | 北京邮电大学 | Satellite navigation system principle simulation method and satellite signal simulator |
| EP2745141A1 (en) * | 2011-08-16 | 2014-06-25 | European Space Agency | A navigation system using spreading codes based on pseudo-random noise sequences |
Non-Patent Citations (3)
| Title |
|---|
| 《卫星导航定位系统中伪随机码的研究》;卢晓春等;《时间频率学报》;20040630;全文 * |
| 《基于卫星配对的北斗系统扩频码优化分配》;朱建锋等;《导航定位学报》;20130930;全文 * |
| 《面向卫星导航接收机应用的自动增益控制环路设计》;王会;《硕士论文》;20071231;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106253946A (en) | 2016-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106253946B (en) | A kind of satellite navigation signals data/pilot branch spread spectrum code optimization matching method | |
| CA2578018C (en) | Point positioning apparatus and point positioning method | |
| CN106932787A (en) | Beidou satellite navigation signal large-scale parallel real-time Transmission method and system | |
| CN103424754B (en) | Based on synthesizing the MBOC modulation signal of related function without fuzzy multipaths restraint method | |
| Yao et al. | Quadrature multiplexed BOC modulation for interoperable GNSS signals | |
| CN103558612B (en) | Phaselocked loop combines phase demodulation track loop with subcarrier ring | |
| CN105607096A (en) | Dual-satellite time difference and frequency difference positioning method and positioning device | |
| CN103558615B (en) | FLL combines frequency discrimination track loop with subcarrier ring | |
| CN106842263A (en) | Based on the satellite navigation method that navigation signal space time processing is combined with vector tracking | |
| US20150063426A1 (en) | Systems and methods for pseudo-random coding | |
| CN107167820A (en) | A kind of digital demultiplexing satellite navigation signal simulator, method and detecting system | |
| CN103869337A (en) | Signal tracking method of global navigation satellite system, signal tracking apparatus and satellite signal reciever having the same | |
| Hodgart et al. | Improvements and additions to the double estimation technique | |
| CN108008422A (en) | Signal capture apparatus and method when pseudo satellite, pseudolite is jumped | |
| Fantino et al. | Performance analysis of MBOC, AltBOC and BOC modulations in terms of multipath effects on the carrier tracking loop within GNSS receivers | |
| CN105372678B (en) | A kind of sine BOC modulated signal without fuzzy tracking method | |
| CN102012499A (en) | Chinese terrestrial digital television single frequency network-based positioning method and system | |
| US9784815B2 (en) | Separating ranging and data signals in a wireless positioning system | |
| CN106526630B (en) | Half cycle fuzziness removing method and device | |
| US20080191932A1 (en) | Method and Device for Correlation Detection in Spread Spectrum Transmission Systems by Fast Fourier Transformation | |
| JP5765902B2 (en) | Code tracking device for satellite signals | |
| JP5607606B2 (en) | Multipath detection apparatus and GNSS reception apparatus | |
| CN102375149B (en) | Multi-associated-branch mixed spacing delay lock loop and application method thereof | |
| CN103439715A (en) | BOC signal anti-multipath method based on two-dimensional tracking structure | |
| CN107479071B (en) | Pseudolite receiver and Pseudolite signal method of reseptance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |