CN105204036A - Tracking loop design method based on multi-path inhibiting signals - Google Patents

Tracking loop design method based on multi-path inhibiting signals Download PDF

Info

Publication number
CN105204036A
CN105204036A CN201510586803.4A CN201510586803A CN105204036A CN 105204036 A CN105204036 A CN 105204036A CN 201510586803 A CN201510586803 A CN 201510586803A CN 105204036 A CN105204036 A CN 105204036A
Authority
CN
China
Prior art keywords
signal
represent
navigation
carrier wave
tracking loop
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201510586803.4A
Other languages
Chinese (zh)
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.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of 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 Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN201510586803.4A priority Critical patent/CN105204036A/en
Publication of CN105204036A publication Critical patent/CN105204036A/en
Pending legal-status Critical Current

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/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/22Multipath-related issues
    • 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/428Determining position using multipath or indirect path propagation signals in position determination

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)
  • Radio Relay Systems (AREA)

Abstract

The invention relates to a method for reducing errors in the message demodulation and navigation process, in particular to a tracking loop design method based on multi-path inhibiting signals. The method comprises the following steps that a satellite navigation signal received through an antenna of a navigation receiver is given, digital down conversion is conducted on the satellite navigation signal received through the navigation receiver, an intermediate frequency signal is obtained, quantification is conducted on the intermediate frequency signal, and a quantification signal is obtained; the quantification signal passes through a filter, autocorrecting computing is conducted on the quantification signal, and a sequence input to an equalizer can be obtained; after the intermediate frequency signal is captured through the navigation receiver, the intermediate frequency signal conducts computing with an I by-pass carrier wave and a Q by-pass carrier wave which are output through a carrier-tracking loop, and carrier wave stripping is completed; the intermediate frequency signal which is subjected to carrier wave stripping conducts channel compensation by means of the equalizer which is input with the sequence, interrelated tracking is conducted on a code tracking loop and the carrier-tracking loop, and message demodulation and navigation are achieved. By means of the tracking loop design method, the tracking effect of a navigation machine is more accurate, and the positioning accuracy is improved.

Description

A kind of track loop method for designing based on suppressing multipath signal
Technical field
The present invention relates to the method reducing error in demodulation navigation message process, being specially a kind of track loop method for designing based on suppressing multipath signal.
Background technology
As one of information industry with fastest developing speed, satellite navigation industry is more and more subject to the attention of country.Satellite navigation signals, as a kind of information resources shared for hundreds of thousands of subscribers, for vast user, needs to determine that the satellite receiver of locating information is caught, track navigation satellite-signal.Therefore, navigation neceiver plays vital effect positioning in the process resolved.
Can draw from current research, multipath effect has become the main source of error that navigation neceiver is measured, its there is the difficult problem of time variation and environmental characteristics.The harm of multipath effect, first can make to be modulated to the ranging code on navigation signal and navigation data distortion, but also the phase place of carrier wave can be caused to distort; The quality that also can affect each measurement data reduces; Situation the most serious is, multipath signal can cause the loop losing lock of receiver.
Summary of the invention
The problem of error can be there is in the present invention under multipath signal interference to solve navigation neceiver track loop, provide one and take balancing technique in track loop, its equalizer coefficients upgrades a kind of track loop method for designing based on suppressing multipath signal taked based on CMA algorithm.
The present invention adopts following technical scheme to realize: a kind of track loop method for designing based on suppressing multipath signal, comprises the following steps:
The satellite navigation signals s that given navigation neceiver antenna receives (i)(t), in formula, i represents the numbering of satellite, and subscript I, Q are respectively two branch roads, and A represents signal amplitude, and C represents ranging code, and D represents the numeric data code that pseudo-code is modulated, f 0represent carrier frequency, represent the initial phase of carrier wave, N'(t) represent noise signal;
Satellite navigation signals s (i)t () obtains intermediate-freuqncy signal s after Digital Down Convert (i)' (t), s ( i ) ′ ( t ) = Σ i = 0 N ∂ i A C ( t - τ i ) D ( t - τ i ) c o s [ 2 π ( f I F + f d ) t + φ i ] + N ′ ( t ) , In formula, i=0 represents satellite direct signal, and i ≠ 0 represents the multipath signal on other roads, and N is the quantity of multipath signal, and A represents signal amplitude, represent signal amplitude fading coefficients, D (t) represents the numeric data code of navigation message, and C (t) represents C/A code, τ irepresent chip delay, f iFrepresent the frequency of intermediate-freuqncy signal, f drepresent Doppler frequency deviation, represent i-th signal phase, N'(t) represent noise signal;
To intermediate-freuqncy signal s (i)' (t) carry out quantification and obtain quantized signal s [n], t in formula sfor sampling interval;
Quantized signal s [n] device and auto-correlation computation after filtering, can obtain the sequence W (n+1) of input equalizer, carry out signal equalization on this basis, W (n+1) iterative formula is as follows:
W ( n + 1 ) = W ( n ) - α f ( n ) z ( n ) [ z ( n ) 2 - | a ^ ( n ) | 2 ] X * ( n ) ,
f ( n ) = 1 , sgn [ z ( n ) 2 - | α ^ ( n ) | 2 ] = sgn [ z ( n ) 2 - R ] 0 , sgn [ z ( n ) 2 - | α ( n ) ^ | 2 ] ≠ sgn [ z ( n ) 2 - R ] ; W (n+1) represents the equalized sequence in n+1 moment, and W (n) represents the equalized sequence in n moment, and α represents equalizing coefficient, and z (n) represents the signal quantization sequence in n moment, represent the modulus value of the signal quantization sequence in n moment, X *n () represents satellite-signal, R represents threshold value;
Navigation neceiver completes after intermediate-freuqncy signal catches, and intermediate-freuqncy signal enters track loop and starts working, and the I branch road carrier wave that intermediate-freuqncy signal exports with carrier tracking loop respectively and Q branch road carrier wave carry out computing respectively, complete stripping carrier wave;
Complete and peel off the intermediate-freuqncy signal of carrier wave and utilize input to have the balanced device of sequence W (n+1) to carry out channel compensation, finally enter code tracking loop respectively and carrier tracking loop carries out correlation tracking, realize demodulation navigation message.
This loop effectively can weaken multipath effect, when there is multi-path jamming in channel, introduce the tracking accuracy of code tracking precision higher than typical early slow DLL loop of this equalization algorithm in track loop, thus improve the measuring accuracy of pseudorange, and then the positioning precision of navigation neceiver can be improved.
Accompanying drawing explanation
The classical code tracking loop block diagram of Fig. 1.
The classical carrier tracking loop block diagram of Fig. 2.
Fig. 3 is inventive receiver track loop block diagram.
Embodiment
Based on the track loop method for designing suppressing multipath signal, comprise the following steps:
The satellite navigation signals s that given navigation neceiver antenna receives (i)(t), in formula, i represents the numbering of satellite, and subscript I, Q are respectively two branch roads, and A represents signal amplitude, and C represents ranging code, and D represents the numeric data code that pseudo-code is modulated, f 0represent carrier frequency, represent the initial phase of carrier wave, N'(t) represent noise signal;
The satellite navigation signals s of the civil navigation receiver reception of single-frequency (i)t () obtains intermediate-freuqncy signal s after Digital Down Convert (i)' (t), s ( i ) ′ ( t ) = Σ i = 0 N ∂ i A C ( t - τ i ) D ( t - τ i ) c o s [ 2 π ( f I F + f d ) t + φ i ] + N ′ ( t ) , In formula, i=0 represents satellite direct signal, and i ≠ 0 represents the multipath signal on other roads, and N is the quantity of multipath signal, and A represents signal amplitude, represent signal amplitude fading coefficients, D (t) represents the numeric data code of navigation message, and C (t) represents C/A code, τ irepresent chip delay, f iFrepresent the frequency of intermediate-freuqncy signal, f drepresent Doppler frequency deviation, represent i-th signal phase, N'(t) represent noise signal;
To intermediate-freuqncy signal s (i)' (t) carry out quantification and obtain quantized signal s [n], t in formula sfor sampling interval;
Quantized signal s [n] device and auto-correlation computation after filtering, can obtain the sequence W (n+1) of input equalizer, carry out signal equalization on this basis, W (n+1) iterative formula is as follows:
W ( n + 1 ) = W ( n ) - α f ( n ) z ( n ) [ z ( n ) 2 - | a ^ ( n ) | 2 ] X * ( n ) ,
f ( n ) = 1 , sgn [ z ( n ) 2 - | α ^ ( n ) | 2 ] = sgn [ z ( n ) 2 - R ] 0 , sgn [ z ( n ) 2 - | α ( n ) ^ | 2 ] ≠ sgn [ z ( n ) 2 - R ] ; W (n+1) represents the equalized sequence in n+1 moment, and W (n) represents the equalized sequence in n moment, αrepresent equalizing coefficient, z (n) represents the signal quantization sequence in n moment, represent the modulus value of the signal quantization sequence in n moment, X *n () represents satellite-signal, R represents threshold value;
Navigation neceiver completes after intermediate-freuqncy signal catches, and intermediate-freuqncy signal enters track loop and starts working, and the I branch road carrier wave that intermediate-freuqncy signal exports with carrier tracking loop respectively and Q branch road carrier wave carry out computing respectively, complete stripping carrier wave;
Complete and peel off the intermediate-freuqncy signal of carrier wave and utilize input to have the balanced device of sequence W (n+1) to carry out channel compensation, finally enter code tracking loop respectively and carrier tracking loop carries out correlation tracking, realize demodulation navigation message.

Claims (1)

1., based on the track loop method for designing suppressing multipath signal, it is characterized in that comprising the following steps:
The satellite navigation signals s that given navigation neceiver antenna receives (i)(t),
in formula, i represents the numbering of satellite, and subscript I, Q are respectively two branch roads, and A represents signal amplitude, and C represents ranging code, and D represents the numeric data code that pseudo-code is modulated, f 0represent carrier frequency, represent the initial phase of carrier wave, N'(t) represent noise signal;
Satellite navigation signals s (i)t () obtains intermediate-freuqncy signal s after Digital Down Convert (i)' (t),
S ( i ) ′ ( t ) = Σ i = 0 N ∂ i A C ( t - τ i ) D ( t - τ i ) c o s [ 2 π ( f I F + f d ) t + φ i ] + N ′ ( t ) , In formula, i=0 represents satellite direct signal, and i ≠ 0 represents the multipath signal on other roads, and N is the quantity of multipath signal, and A represents signal amplitude, represent signal amplitude fading coefficients, D (t) represents the numeric data code of navigation message, and C (t) represents C/A code, τ irepresent chip delay, f iFrepresent the frequency of intermediate-freuqncy signal, f drepresent Doppler frequency deviation, represent i-th signal phase, N'(t) represent noise signal;
To intermediate-freuqncy signal s (i)' (t) carry out quantification and obtain quantized signal s [n],
t in formula sfor sampling interval;
Quantized signal s [n] device and auto-correlation computation after filtering, can obtain the sequence W (n+1) of input equalizer, carry out signal equalization on this basis, W (n+1) iterative formula is as follows:
W ( n + 1 ) = W ( n ) - α f ( n ) z ( n ) [ z ( n ) 2 - | a ^ ( n ) | 2 ] X * ( n ) ,
f ( n ) = 1 , sgn [ z ( n ) 2 - | α ^ ( n ) | 2 ] = sgn [ z ( n ) 2 - R ] 0 , sgn [ z ( n ) 2 - | α ( n ) ^ | 2 ] ≠ sgn [ z ( n ) 2 - R ] , W (n+1) represents the equalized sequence in n+1 moment, and W (n) represents the equalized sequence in n moment, and α represents equalizing coefficient, and z (n) represents the signal quantization sequence in n moment, represent the modulus value of the signal quantization sequence in n moment, X *n () represents satellite-signal, R represents threshold value;
Navigation neceiver completes after intermediate-freuqncy signal catches, and intermediate-freuqncy signal enters track loop and starts working, and the I branch road carrier wave that intermediate-freuqncy signal exports with carrier tracking loop respectively and Q branch road carrier wave carry out computing respectively, complete stripping carrier wave;
Complete and peel off the intermediate-freuqncy signal of carrier wave and utilize input to have the balanced device of sequence W (n+1) to carry out channel compensation, finally enter code tracking loop respectively and carrier tracking loop carries out correlation tracking, realize demodulation navigation message.
CN201510586803.4A 2015-09-15 2015-09-15 Tracking loop design method based on multi-path inhibiting signals Pending CN105204036A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510586803.4A CN105204036A (en) 2015-09-15 2015-09-15 Tracking loop design method based on multi-path inhibiting signals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510586803.4A CN105204036A (en) 2015-09-15 2015-09-15 Tracking loop design method based on multi-path inhibiting signals

Publications (1)

Publication Number Publication Date
CN105204036A true CN105204036A (en) 2015-12-30

Family

ID=54951814

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510586803.4A Pending CN105204036A (en) 2015-09-15 2015-09-15 Tracking loop design method based on multi-path inhibiting signals

Country Status (1)

Country Link
CN (1) CN105204036A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106597481A (en) * 2016-12-12 2017-04-26 太原理工大学 Vector tracking multi-path interference suppression algorithm based on blind equalizer
CN106906728A (en) * 2017-02-14 2017-06-30 北京市道路工程质量监督站 Detection method and device are peeled off between a kind of asphalt pavement structural layer
CN110118979A (en) * 2018-11-26 2019-08-13 太原理工大学 The method of improved differential evolution algorithm estimation multipath parameter based on broad sense cross-entropy
CN111880200A (en) * 2020-08-18 2020-11-03 南京航空航天大学 Multi-path inhibition method based on steepest descent method
CN113297780A (en) * 2021-07-27 2021-08-24 湖南卫导信息科技有限公司 High-speed simulation method, device, medium and equipment for carrier shielding and multipath signals

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6151368A (en) * 1999-03-22 2000-11-21 Sicom, Inc. Phase-noise compensated digital communication receiver and method therefor
CN101232358A (en) * 2007-01-23 2008-07-30 复旦大学 Method and system of channel blind equalization in MIMO system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6151368A (en) * 1999-03-22 2000-11-21 Sicom, Inc. Phase-noise compensated digital communication receiver and method therefor
CN101232358A (en) * 2007-01-23 2008-07-30 复旦大学 Method and system of channel blind equalization in MIMO system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐金标等: ""基于CMA算法的双模式盲均衡算法"", 《通信学报》 *
李春宇等: ""基于盲判决反馈均衡的GPS信号多径干扰抑制方法"", 《电子学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106597481A (en) * 2016-12-12 2017-04-26 太原理工大学 Vector tracking multi-path interference suppression algorithm based on blind equalizer
CN106906728A (en) * 2017-02-14 2017-06-30 北京市道路工程质量监督站 Detection method and device are peeled off between a kind of asphalt pavement structural layer
CN106906728B (en) * 2017-02-14 2019-08-16 北京市道路工程质量监督站 Detection method and device are removed between a kind of asphalt pavement structural layer
CN110118979A (en) * 2018-11-26 2019-08-13 太原理工大学 The method of improved differential evolution algorithm estimation multipath parameter based on broad sense cross-entropy
CN110118979B (en) * 2018-11-26 2023-02-28 太原理工大学 Method for estimating multipath parameters by using improved differential evolution algorithm based on generalized mutual entropy
CN111880200A (en) * 2020-08-18 2020-11-03 南京航空航天大学 Multi-path inhibition method based on steepest descent method
WO2022037309A1 (en) * 2020-08-18 2022-02-24 南京航空航天大学 Multipath suppression method based on steepest descent method
US11716106B1 (en) 2020-08-18 2023-08-01 Nanjing University Of Aeronautics And Astronautics Multipath suppression method based on steepest descent method
CN113297780A (en) * 2021-07-27 2021-08-24 湖南卫导信息科技有限公司 High-speed simulation method, device, medium and equipment for carrier shielding and multipath signals
CN113297780B (en) * 2021-07-27 2021-10-12 湖南卫导信息科技有限公司 High-speed simulation method, device, medium and equipment for carrier shielding and multipath signals

Similar Documents

Publication Publication Date Title
CN105204036A (en) Tracking loop design method based on multi-path inhibiting signals
TWI483592B (en) Continuous sequential scatterer estimation
WO2018107441A1 (en) Signal capturing method and receiver for satellite navigation system
CN101807940B (en) Anti-multipath interference device of GNSS receiving system and method thereof
Wang et al. Doppler estimation and timing synchronization of underwater acoustic communication based on hyperbolic frequency modulation signal
CN101997807B (en) Channel estimation method and device
CN103338024B (en) The complementary Kalman filtering apparatus and method of time delay in antenna array
CN104836759B (en) Sane underwater acoustic channel method of estimation based on joint norm optimization
KR101828790B1 (en) Frequency shift keying signal receiving method and device
CN102307054B (en) Direct sequence spread spectrum signal acquisition method
CN102571652A (en) Method for estimating global navigation satellite system (GNSS) baseband signal
CN103428153A (en) Gaussian minimum shift keying (GMSK) signal receiving method in satellite mobile communication
CN102164110B (en) Method and system for balancing frequency domain
CN104765052A (en) GEO navigation satellite high-sensitivity carrier tracking method
CN111190197A (en) Navigation satellite signal quality on-orbit optimization and maintenance method
CN1256811A (en) Method and apparatus for compensating for click noise in FM receiver
CN101630000B (en) System for estimating influence of interference signal on GPS performance
CN103916199B (en) The time delay of a kind of aerial signal and phase adjusting apparatus and method
CN106330342A (en) Water sound communication doppler factor estimation method with low computation complexity
CN103023830A (en) Turbo iteration detecting device and method directing at big carrier frequency excursion
KR101421156B1 (en) Apparatus and Method for acquiring GPS Signal using Multi-stage Partial Cross-Correlator
CN102565822B (en) Capture method of GPS L5 signal and apparatus thereof
CN109738916A (en) A kind of multipath parameter estimation method based on compressed sensing algorithm
JPWO2007020693A1 (en) GPS positioning method and GPS positioning device
CN110780320A (en) Software and hardware integrated satellite navigation signal processing method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20151230

RJ01 Rejection of invention patent application after publication