CN104459743B - Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component - Google Patents

Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component Download PDF

Info

Publication number
CN104459743B
CN104459743B CN201410704995.XA CN201410704995A CN104459743B CN 104459743 B CN104459743 B CN 104459743B CN 201410704995 A CN201410704995 A CN 201410704995A CN 104459743 B CN104459743 B CN 104459743B
Authority
CN
China
Prior art keywords
navigation signal
signal
carrier phase
component
signal component
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
Application number
CN201410704995.XA
Other languages
Chinese (zh)
Other versions
CN104459743A (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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201410704995.XA priority Critical patent/CN104459743B/en
Publication of CN104459743A publication Critical patent/CN104459743A/en
Application granted granted Critical
Publication of CN104459743B publication Critical patent/CN104459743B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry

Abstract

The present invention proposes carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component, the assessment for New Complex modulation system navigation signal is realized by building plural number local code, simultaneously building plural number local code when, add intermodulation component.Compared with component of signal carrier phase bias estimation method compare, method disclosed by the invention is applicable not only to the simple modulation signals such as BPSK, QPSK, it is also applied for the coherent multicarrier modulation navigation signal included containing single sideband modulated signal component, such as AltBoc, Td AltBoc, TMBOC, DualQPSK etc., can quantitatively provide carrier phase deviation between the component of signal of various modulation system navigation signal, and there is higher precision.

Description

Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component
Technical field
The present invention relates to field of satellite navigation, relate to carrier phase deviation between a kind of coherent multicarrier modulated signal component true Determine method.
Background technology
In the development process of aeronautical satellite payload, correctness, investigation signal for checking navigation signal generated The distortion that journey introduces, needs to test property indices, and wherein, component of signal phase deviation is that satellite navigation effectively carries The important indicator of lotus, it be related to the phase relation of the carrier phase of each road component of signal of same frequency whether with the phase of nominal Position relationship consistency.
At present, the assessment of component of signal phase deviation index can be carried out by the following method:
Method 1, the method utilizing commonality vector signal analyzer;
Method 2, utilizes receiver to measure component of signal phase deviation;
Method 3, the method for planisphere.
But, there are the following problems for said method:
(1) method 1 is only capable of assessment simple modulation signal, such as QPSK, component of signal phase deviation, but for relevant The navigation signal of multi-carrier modulation, such as AltBoc, Td-AltBoc signal, it is impossible to be estimated.
(2) method 2 uses receiver to measure the carrier phase of each road component of signal, is calculated signal on this basis The carrier phase deviation of component, the method is suitable only for the center frequency point of component of signal and whole navigation signal center frequency point phase Same situation, it is impossible to measure the situation that the center frequency point of component of signal is different with the center frequency point of whole navigation signal.
(3) problem that method 3 exists: planisphere method reflects component of signal carrier phase deviation more intuitively, but Index cannot be quantified.
Summary of the invention
Present invention solves the technical problem that and be: overcome the deficiencies in the prior art, it is provided that a kind of coherent multicarrier modulated signal Various types of navigation can be estimated by carrier phase offset determination methods between component, are particularly suited for coherent multicarrier and adjust The novel navigation signal of system.
The technical scheme is that carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component, Step is as follows:
1) sample frequency f is sets, navigation signal to be detected is sampled, the navigation signal after being sampled;
2) navigation signal after sampling is carried out quadrature demodulation, obtain I road branch road baseband signal in the same direction and Q road quadrature branch Baseband signal;
3) obtain the prior information of the navigation signal modulation system of outside input, build the plural number of each navigation signal component Local code;
4) according to step 3) the plural local code of each navigation signal component that obtains and step 2) obtain orthogonal The I road that obtains after demodulation, Q roadbed band signal, use maximal possibility estimation, it is thus achieved that each navigation signal component starts in sampling The carrier phase value in moment;
5) by step 4) carrier phase of all navigation signal components that obtains is the poorest, obtains each navigation signal component Between carrier phase relationship, as actual measurement carrier phase relationship;
6) the actual measurement carrier phase relationship between each navigation signal component deducts nominal carrier phase relation, obtains navigation letter Carrier phase deviation between number component.
Step 1) in time navigation signal to be detected is sampled, need to meet following five requirements simultaneously:
Wherein BnavFrequency band range shared by navigation signal, n is positive integer;
2. sample frequency is not less than 500MHz;
3. sample frequency can not be the integral multiple of the spreading rate of pseudo-code in navigation signal;
4. the sampling time is more than the twice of PN-code capture length;
5. the clock that sampling clock and signal generate meets and has identical clock source and lead to.
Step 3) described in build each navigation signal component plural local code, including following three kinds of situations:
Situation one: for the center frequency point of this navigation signal component identical with whole navigation signal center frequency point in the case of, Then utilize the I road of this signal and Q road spread-spectrum pseudorandom codes as the real part of plural number local code and imaginary part, build plural number local code;
Situation two: the center frequency point for this navigation signal component is positioned at the symmetrical of whole navigation signal center frequency point The situation on both sides, then utilize the pseudo-code of this navigation signal be multiplied with subcarrier as plural number local code real part, imaginary part is zero, structure Build plural number local code;
Situation three: the center frequency point for this navigation signal component is different with whole navigation signal center frequency point, and not with Symmetrical situation centered by it, then utilize the pseudo-code of this signal and the long-pending real part of the multiple subcarrier of single sideband modulation and imaginary part Respectively as real part and imaginary part, build plural number local code.
Step 3) described in build each navigation signal component plural local code, for containing intermodulation component navigation believe Number, the plural local code building process of each road navigation signal component adds intermodulation component.
Present invention advantage compared with prior art is:
(1) method disclosed by the invention realizes for New Complex modulation system navigation letter by building plural number local code Number assessment, be applicable not only to the simple modulation signals such as BPSK, QPSK, be also applied for including containing single sideband modulated signal component Coherent multicarrier modulation navigation signal, such as AltBoc, Td-AltBoc, TMBOC, DualQPSK etc., can quantitatively be given various Carrier phase deviation between the component of signal of modulation system navigation signal;
(2) the method is building plural number local code when, adds intermodulation component, leads for coherent multicarrier modulation Between the component of signal of boat signal, the estimation of carrier phase deviation is higher.
Accompanying drawing explanation
Fig. 1 is the inventive method flow chart;
Fig. 2 follows the tracks of the carrier phase tracking knot of Galileo E5 frequency AltBoc (15,10) signal respectively for utilizing receiver Really;
Fig. 3 is the planisphere of Galileo E5 frequency AltBoc (15,10) ideal signal;
Fig. 4 is the planisphere of Galileo E5 frequency AltBoc (15,10) measured signal (signal to noise ratio 54dB);
Fig. 5 is the planisphere of Galileo E5 frequency AltBoc (15,10) measured signal (signal to noise ratio 30dB);
Fig. 6 is the cost function figure during estimation carrier phase.
Detailed description of the invention
Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component, it is characterised in that step is as follows:
1) sample frequency f is sets, navigation signal to be detected is sampled, the navigation signal after being sampled;Treat The navigation signal of detection is when sampling, and need to meet following five requirements simultaneously:
Wherein BnavFrequency band range shared by navigation signal, n is positive integer;
2. sample frequency is not less than 500MHz;
3. sample frequency can not be the integral multiple of the spreading rate of pseudo-code in navigation signal;
4. the sampling time is more than the twice of PN-code capture length;
The most same source sampling, i.e. sampling clock and signal generate clock meet have identical clock source lead to.
2) navigation signal after sampling is carried out quadrature demodulation, obtain I roadbed band signal and Q roadbed band signal;
3) obtain the prior information of the navigation signal modulation system of outside input, build the plural number of each navigation signal component Local code;Including following three kinds of situations:
Situation one: for the center frequency point of component of signal identical with whole navigation signal center frequency point in the case of, such as QPSK Signal, then utilize the I road of QPSK signal and Q road spread-spectrum pseudorandom codes as the real part of plural number local code and imaginary part, build plural number Local code;
Situation two: the center frequency point of component of signal is positioned to the symmetrical both sides of whole navigation signal center frequency point Situation, such as BOC, CBOC, TMBOC signal, then utilizes pseudo-code to be multiplied with subcarrier the real part as plural number local code, and imaginary part is Zero, build plural number local code;
Situation three: the center frequency point for component of signal is different with whole navigation signal center frequency point, and with it be not The symmetrical situation of the heart, such as single sideband modulated signals such as AltBoc, Td-AltBoc, then utilizes pseudo-code and single sideband modulation multiple The long-pending real part of subcarrier and imaginary part, respectively as real part and imaginary part, build plural number local code.
The described plural local code building each component of signal, for the navigation signal containing intermodulation component, each The plural local code building process of road navigation signal component adds intermodulation component.
4) according to step 3) the plural local code of each navigation signal component that obtains and step 2) obtain orthogonal The I road that obtains after demodulation, Q roadbed band signal, use maximal possibility estimation, it is thus achieved that each navigation signal component starts in sampling The carrier phase value in moment;
5) by step 4) carrier phase of all navigation signal components that obtains is the poorest, obtains each navigation signal component Between carrier phase relationship, as actual measurement carrier phase relationship;
6) the actual measurement carrier phase relationship between each navigation signal component deducts nominal carrier phase relation, obtains navigation letter Carrier phase deviation between number component.
For making the purpose of the present invention, technical scheme and advantage clearer, below with reference to accompanying drawing and with AltBoc As a example by (15,10), the present invention is described in more detail.
As a example by Galileo E5 frequency AltBoc (15,10), receiver is utilized to follow the tracks of Galileo E5 frequency AltBoc respectively The carrier phase tracking result of (15,10) signal is as shown in Figure 4, it can be seen that due to E5 frequency AltBoc (15,10) signal Four tunnel component of signals Ad, Ap, Bd, Bp heart frequency the most in the signal, therefore use receiver to receive four tunnel component of signals respectively and obtain To carrier phase only have between Ad, Ap of same frequency or between Bd, Bp, there is comparability, it is impossible to compare between Ad and Bd Carrier phase relationship;
Galileo E5 frequency AltBoc (15,10) ideal signal, the measured signal of signal to noise ratio 54dB and signal to noise ratio 30dB The planisphere of measured signal is respectively shown in Fig. 3, Fig. 4, Fig. 5, it can be seen that planisphere can characterize the phase place of component of signal and close System, but cannot quantitative result, and when signal to noise ratio is relatively low, planisphere seriously obscures;
For Galileo E5 frequency AltBoc (15,10) signal, the operating procedure of method disclosed by the invention is as follows:
(1) sample frequency f is sets=1GHz, carries out same source sampling, sampling time 2ms, letter before sampling to signal to be assessed Number center frequency point be 1191.795MHz, the center frequency point of the signal after sampling is 191.795MHz;
(2) carry out quadrature demodulation to the navigation signal after sampling, obtain I roadbed band signal IR(t) and Q roadbed band signal QR(t);
(3) utilizing the prior information of AltBoc (15,10), the plural number building tetra-component of signals of Ad, Ap, Bd, Bp is local Code, as a example by Ad, the center frequency point of this component of signal is different with whole navigation signal center frequency point, and not left and right centered by it Symmetry, utilize the pseudo noise code of Ad and the long-pending real part of the multiple subcarrier of single sideband modulation and imaginary part respectively as real part and imaginary part, Structure plural number local code:
SLocal(t)=IL(t)+jQL(t)=CAd(t)×{Scos(t)+jSsin(t)}
Wherein, SLocalT () represents plural number local code, ILT () represents the real part of plural number local code, QLT () represents that plural number is originally The imaginary part of ground code, CAdT () represents the pseudo noise code of Ad, Scos(t) and Ssin(t) be respectively AltBoc (15,10) single-side belt adjust The real part of the multiple subcarrier of system and imaginary part;
(4) according to the I road obtained after the plural local code of each component of signal and quadrature demodulation, Q roadbed band signal, estimate Counting each component of signal carrier phase value in sampling start time, use maximal possibility estimation, its cost function is
Wherein,Represent cost function, ILT () represents the real part of plural number local code, QLT () represents plural number local code Imaginary part, IRT () represents the real part after sampled signal quadrature frequency conversion, QRT () represents the void after sampled signal quadrature frequency conversion Portion;Sin and cos represents sinusoidal and cos operation respectively, d andDuring representing pseudo-code time delay and quadrature frequency conversion respectively The local oscillator initial phase used;
WhenDuring acquirement maximumIt is the carrier phase estimated value of component of signal.For improving estimated accuracy, WillLogarithm operation is carried out again, as shown in Figure 6 after inverted.
(5) utilize cost function, use optimization method to estimateValue, uses during i.e. can get quadrature frequency conversion Local oscillator initial phase, using this value as the carrier phase of this road signal, after estimating the carrier phase of each component of signal, institute The carrier phase having component of signal is the poorest, calculates the actual measurement carrier phase relationship between each component of signal;
(6) the actual measurement carrier phase relationship between each component of signal deducts nominal carrier phase relation and obtains between component of signal Carrier phase deviation.
The non-detailed description of the present invention is known to the skilled person technology.

Claims (3)

1. carrier phase offset determination methods between a coherent multicarrier modulated signal component, it is characterised in that step is as follows:
1) sample frequency f is sets, navigation signal to be detected is sampled, the navigation signal after being sampled;To be detected Navigation signal when sampling, need to meet following five requirements simultaneously:
Wherein BnavFrequency band range shared by navigation signal, n is positive integer;
2. sample frequency is not less than 500MHz;
3. sample frequency can not be the integral multiple of the spreading rate of pseudo-code in navigation signal;
4. the sampling time is more than the twice of PN-code capture length;
5. the clock that sampling clock and signal generate meets and has identical clock source and lead to;
2) navigation signal after sampling is carried out quadrature demodulation, obtain I road branch road baseband signal in the same direction and Q road quadrature branch base band Signal;
3) obtaining the prior information of the navigation signal modulation system of outside input, the plural number building each navigation signal component is local Code;
4) according to step 3) the plural local code of each navigation signal component that obtains and step 2) quadrature demodulation that obtains After obtain I road, Q roadbed band signal, use maximal possibility estimation, it is thus achieved that each navigation signal component sampling start time Carrier phase value;
5) by step 4) carrier phase of all navigation signal components that obtains is the poorest, obtains between each navigation signal component Carrier phase relationship, as actual measurement carrier phase relationship;
6) the actual measurement carrier phase relationship between each navigation signal component deducts nominal carrier phase relation, obtains navigation signal and divides Carrier phase deviation between amount.
Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component the most according to claim 1, its Be characterised by: step 3) described in build each navigation signal component plural local code, including following three kinds of situations:
Situation one: for the center frequency point of this navigation signal component identical with whole navigation signal center frequency point in the case of, then profit By the I road of this signal and Q road spread-spectrum pseudorandom codes as the real part of plural number local code and imaginary part, build plural number local code;
Situation two: the center frequency point of this navigation signal component is positioned to the symmetrical both sides of whole navigation signal center frequency point Situation, then utilize the pseudo-code of this navigation signal to be multiplied with subcarrier the real part as plural number local code, imaginary part is zero, builds multiple Number local code;
Situation three: the center frequency point for this navigation signal component is different with whole navigation signal center frequency point, and with it is not The symmetrical situation in center, then utilize the pseudo-code of this signal and the long-pending real part of the multiple subcarrier of single sideband modulation and imaginary part respectively As real part and imaginary part, build plural number local code.
Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component the most according to claim 1, its Be characterised by: step 3) described in build each navigation signal component plural local code, for the navigation containing intermodulation component Signal, adds intermodulation component in the plural local code building process of each road navigation signal component.
CN201410704995.XA 2014-11-27 2014-11-27 Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component Active CN104459743B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410704995.XA CN104459743B (en) 2014-11-27 2014-11-27 Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410704995.XA CN104459743B (en) 2014-11-27 2014-11-27 Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component

Publications (2)

Publication Number Publication Date
CN104459743A CN104459743A (en) 2015-03-25
CN104459743B true CN104459743B (en) 2017-01-04

Family

ID=52906092

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410704995.XA Active CN104459743B (en) 2014-11-27 2014-11-27 Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component

Country Status (1)

Country Link
CN (1) CN104459743B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105116425B (en) * 2015-08-21 2017-07-28 西安空间无线电技术研究所 A kind of parallel AltBOC navigation signals intermediate frequency generation method
CN105549043B (en) * 2015-12-16 2018-02-09 西安空间无线电技术研究所 Carrier phase relationship detection method between a kind of each component of satellite navigation signals
CN106226793B (en) * 2016-07-29 2019-01-15 北京空间飞行器总体设计部 A kind of in-orbit navigation signal IQ phase equalization scaling method
CN109507703B (en) * 2018-12-17 2023-06-16 成都国星通信有限公司 GNSS carrier phase calculation method
CN112578422A (en) * 2019-09-27 2021-03-30 清华大学 Composite navigation signal receiving method and receiver
CN111521887B (en) * 2020-05-07 2022-06-07 西安热工研究院有限公司 380V power supply-based high-voltage station service power supply nuclear phase system and method
CN115801524B (en) * 2023-01-30 2023-05-16 中国电子科技集团公司第五十四研究所 Low signal-to-noise ratio carrier recovery method based on multi-symbol correlation and maximum likelihood selection

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5343209A (en) * 1992-05-07 1994-08-30 Sennott James W Navigation receiver with coupled signal-tracking channels
CN101039126B (en) * 2006-03-16 2012-01-25 湖北众友科技实业股份有限公司 Method and system of phase deviation measurement for testing TD-SCDMA terminal
CN101567709B (en) * 2009-05-27 2012-10-03 西华大学 Method and device for weakening the influence of multipath on positioning accuracy of receiver antenna
CN103278825B (en) * 2013-05-02 2015-05-27 西安空间无线电技术研究所 Method for determining satellite navigation signal quality evaluation parameters

Also Published As

Publication number Publication date
CN104459743A (en) 2015-03-25

Similar Documents

Publication Publication Date Title
CN104459743B (en) Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component
CN103278825B (en) Method for determining satellite navigation signal quality evaluation parameters
JP5607607B2 (en) GNSS receiver
US7443339B2 (en) Cross-correlation suppression technique for position location receivers
USRE49560E1 (en) Method and a measuring device for measuring broadband measurement signals
RU2658625C1 (en) Spread spectrum signal generating method, generating apparatus, receiving method and receiving apparatus
US7816904B2 (en) Modulation signature trigger
US8762096B2 (en) Trigger generation for digital modulation signal analysis
JP2016526150A5 (en)
CN106918822A (en) Calculate the GNSS receiver of the non-fuzzy discriminator for parsing subcarrier tracking fuzziness
CN105891600A (en) Quaternary phase shift keying signal spectrum estimation method
US20080089401A1 (en) Signal Testing System
CN108562918B (en) BOC (n, n) ambiguity-free capturing method and device based on correlation shift
CN102495417A (en) High-robustness carrier tracking system of global navigation satellite system receiver and method
CN103905360B (en) The non-cooperation bpsk signal coding/decoding method for adding " sentencing pole " to operate
CN107272026B (en) A kind of navigation signal component phase test bias method
RU2524673C1 (en) Method of measuring phase shift variations of centre frequency of minimum frequency shift-keyed signals and apparatus for implementing said method
JP6706579B2 (en) Improvement of satellite positioning method
CN108957492A (en) A kind of L1C/A the and L1C joint acquisition method and device of GPS
JP5344825B2 (en) Frequency offset estimation and C / I ratio measurement method using pilot signal by repetition for TDMA signal having offset larger than Nyquist frequency of reference symbol rate
CN111753248A (en) Frequency offset time vernier measurement method and system and computer readable storage medium
Sun et al. Sub cross-correlation shift cancellation technique for unambiguous tracking of BOC (kn, n) signals
CN110441798A (en) The Beidou RDSS weak signal capturing method for integrating and selecting star to assist based on multiplication accumulation
CN102299880A (en) Method for calculating modulation signal phase characteristic
Sui et al. Analysis and Test of Unambiguous Tracking Performance for High-Order BOC Signals

Legal Events

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