CN104215981A - Adaptive tracking method for receiver in high-dynamic environment - Google Patents

Adaptive tracking method for receiver in high-dynamic environment Download PDF

Info

Publication number
CN104215981A
CN104215981A CN201410433179.XA CN201410433179A CN104215981A CN 104215981 A CN104215981 A CN 104215981A CN 201410433179 A CN201410433179 A CN 201410433179A CN 104215981 A CN104215981 A CN 104215981A
Authority
CN
China
Prior art keywords
carrier
loop
fll
phase
frequency
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
CN201410433179.XA
Other languages
Chinese (zh)
Other versions
CN104215981B (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.)
Sichuan Jiuzhou Electric Group Co Ltd
Original Assignee
Sichuan Jiuzhou Electric Group Co Ltd
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 Sichuan Jiuzhou Electric Group Co Ltd filed Critical Sichuan Jiuzhou Electric Group Co Ltd
Priority to CN201410433179.XA priority Critical patent/CN104215981B/en
Publication of CN104215981A publication Critical patent/CN104215981A/en
Application granted granted Critical
Publication of CN104215981B publication Critical patent/CN104215981B/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/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/29Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, 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
    • 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/23Testing, monitoring, correcting or calibrating of receiver elements

Abstract

The invention discloses an adaptive tracking method for a receiver in a high-dynamic environment. According to the adaptive tracking method, bandwidths of loop filters for a frequency-locked loop, a phase-locked loop and an auxiliary phase-locked loop of the frequency-locked loop are adjusted by detecting the dynamic property and the signal-to-noise property of tracking signals of the receiver in real time, so that optimum values of the bandwidths are kept; meanwhile, a dynamic application range of a tracking loop is widened by adjusting an auxiliary coefficient of the auxiliary phase-locked loop of the frequency-locked loop, so that stability of the tracking loop in the high-dynamic environment is improved. The adaptive tracking method for the receiver in the high-dynamic environment has the advantages that the adaptive tracking method is capable of enabling the receiver to adjust own parameters timely in a dynamic change environment, errors of tracking satellite navigation signals are kept within a small range, and the shortcomings of large tracking errors and proneness to lock loss of the receiver in the prior art are overcome.

Description

Adaptive tracking method under a kind of receiver high dynamic environment
Technical field
The present invention relates to satellite communication receiver, particularly relate to adaptive tracking method under a kind of receiver high dynamic environment.
Background technology
Beidou satellite navigation system (hereinafter referred to as BD) is that China independently builds, independent operating, and the GPS (Global Position System) shared with other satellite navigation system compatibilities of the world.By 5 geostationary orbit (GEO) satellites, circle Earth's orbit (MEO) satellite and 3 inclination geostationary orbit (IGSO) satellites form in 27.Current Beidou satellite navigation system has succeeded in sending up 16 satellites, defines regional service ability, provides the services such as passive location, navigation and time service towards China and periphery most area.
BD satellite-signal has B1, B2 and B3 tri-frequencies, and wherein B1 and B2 is two civil signal.B1 and B2 signal is made up of on carrier wave " ranging code+navigation message " orthogonal modulation of I, Q two branch roads.
The base band signal process of BD receiver mainly comprise catch, follow the tracks of, the synchronous step such as (frame synchronization and bit synchronization) and navigation calculation.Carrier signal, after the carrier signal being received corresponding frequency by antenna, is become digital intermediate frequency signal through down coversion and analog to digital conversion, is supplied to trapping module by BD receiver.Trapping module obtains the carrier frequency of all visible satellites and the coarse value of ranging code phase place by reproduction intermediate frequency carrier and ranging code and the parameter that adjusts them.The phase place of information to the carrier frequency of visible satellite and ranging code that track loop obtains based on trapping module carries out further accurately measuring.When track loop success locking satellite signal, just can carry out bit synchronization, and frame synchronization operation, the textual information that last basis obtains carries out position, speed, time (PVT) resolves.
In BD receiver, track loop is by carrying out accurate tracking to the carrier frequency of BD signal and ranging code phase place, and obtain pseudo-range information according to the ranging code phase place that code tracking loop is measured, from in-phase branch, demodulate navigation message information, thus provide required pseudorange and telegraph text data for navigation calculation.Therefore, the pseudorange accuracy measured of track loop and navigation message demodulation quality directly determine the positioning precision of BD receiver.
In track loop, the performance of loop filter decides the performance of track loop, and determines exponent number and the noise bandwidth of the mainly loop filter of loop filter performance.High exponent number can follow the tracks of the Dynamic Signal of high order.Narrow noise bandwidth contributes to noise in filtering loop, improves the tracking accuracy of loop.Wide noise bandwidth following range is comparatively large, is applicable to following the tracks of dynamically higher signal, but also can increases the noise in loop simultaneously.At present, existing BD receiver tracking loop circuit at least comprises FLL, phaselocked loop and FLL and assists one in phaselocked loop, its exponent number and loop bandwidth set already, limit the dynamic scope of application of loop, under high dynamic environment, tracking error is easily caused to strengthen, the even phenomenon of loop losing lock.Particular problem is as follows:
1) although phaselocked loop is less to carrier tracking error, dynamic property is poor, follows the tracks of carrier phase difficulty comparatively greatly, easy losing lock under high dynamic environment;
2) although FLL has larger following range to Dynamic Signal, tracking error is larger;
3) although FLL assists phaselocked loop can improve the dynamic property of track loop to a certain extent, its fixing initial loop bandwidth limits the ability of loop dynamics and anti-noise jamming.And although can tracking signal preferably based on the track loop of Kalman Filter Technology, its complexity be high, and software simulating cost is high.
Summary of the invention
For the problems referred to above, the present invention proposes adaptive tracking method under a kind of new receiver high dynamic environment, it comprises the following steps:
The intermediate frequency carrier frequency that S100, basis are caught and ranging code phase place, initialization carrier generator and ranging code generator produce the correlation parameter of signal;
S200, the signal pair frequency carrier wave digital signal utilizing carrier generator and ranging code generator to produce carry out carrier wave stripping and ranging code stripping, extract baseband digital signal;
S300, baseband digital signal carried out to phase demodulation, frequency discrimination and carrier-to-noise ratio and calculate;
S400, according to phase demodulation, frequency discrimination and carrier-to-noise ratio result of calculation, ask for according to the mathe-matical map relation preset the auxiliary coefficient that new frequency-locked loop bandwidth sum phase-locked loop bandwidth and FLL assist phaselocked loop;
S500, the auxiliary coefficient of phaselocked loop is assisted to reset the bandwidth of loop filter according to new frequency-locked loop bandwidth sum phase-locked loop bandwidth and FLL;
S600, according to the Output rusults of loop filter adjustment carrier generator and ranging code generator, return step S200.
According to embodiments of the invention, in above-mentioned steps S100, described correlation parameter comprises the phase increment control word of carrier generator and the initial phase of ranging code generator;
Can the phase increment control word △ P of initialization carrier generator according to the following formula f:
Δ P f = 2 M · f carr f clk
In formula, M is the word length of the phase accumulator of carrier generator, f clkfor the clock of carrier generator, f carrfor the intermediate frequency carrier frequency values of catching;
Can the initial phase C of initialization ranging code generator according to the following formula pha0:
C pha0=C pha
In formula, C phafor the ranging code phase value of catching.
According to embodiments of the invention, in above-mentioned steps S300, described phase demodulation comprises yard ring phase demodulation and carrier wave ring phase demodulation, and described frequency discrimination comprises carrier wave ring frequency discrimination.
According to embodiments of the invention, preferably, in above-mentioned steps S300, first can carry out integration to baseband digital signal and add up, and then the calculating of phase demodulation, frequency discrimination and carrier-to-noise ratio is carried out to integration accumulation result.
According to embodiments of the invention, in above-mentioned steps S400, described mathe-matical map relation comprises time interval T kthe average of interior carrier-to-noise ratio result of calculation is as the carrier-to-noise ratio of the tracked satellite-signal of t
CNR ‾ t = 1 K Σ k = 0 K CNR ( k )
In formula, CNR (k) is a kth carrier-to-noise ratio result of calculation, and K is moment t-T kto the number of carrier-to-noise ratio result of calculation between moment t.
According to embodiments of the invention, in above-mentioned steps S400, described mathe-matical map relation comprises the loop bandwidth Bf calculating new i rank FLL according to following formula i':
Bf i ′ = π 2 T 2 · CNR ‾ t · ( V ( i + 1 ) ) 2 · ( λ i ) 2 i 9 F · ( 1 + 1 T · CNR ‾ t ) 2 i + 1
In formula, i represents FLL exponent number, V (i+1)represent the i+1 order derivative of receiver bearer rate, λ ifor the scale-up factor of i rank FLL loop bandwidth and characteristic frequency, for the carrier-to-noise ratio of the tracked satellite-signal of t, T is integration accumulation interval, and F is one and carrier-to-noise ratio relevant parameter.
According to embodiments of the invention, in above-mentioned steps S400, described mathe-matical map relation comprises the loop bandwidth Bp calculating new j rank phaselocked loop according to following formula j':
B p j ′ = 4 CNR ‾ t · ( μ j ) 2 j · ( V ( j ) ) 2 9 ( 1 + 1 2 T · CNR ‾ t ) 2 j + 1
In formula, j represents phaselocked loop exponent number, V (j)represent the j order derivative of receiver bearer rate, μ jfor the scale-up factor of j rank PLL loop bandwidth and characteristic frequency, for the carrier-to-noise ratio of the tracked satellite-signal of t, T is integration accumulation interval.
According to embodiments of the invention, in above-mentioned steps S500, described mathe-matical map relation comprises:
Calculate the sighting distance acceleration a of receiver carrier t twith sighting distance acceleration a t';
a t = 1 K Σ k = 0 K f e ( k ) · ( B f i / λ i ) , a t ′ = 1 K Σ k = 0 K f e ( k ) · ( B f i / λ i ) 2
In formula, f ek () is a kth carrier wave ring frequency discrimination result, K is moment t-T kto the number of the carrier wave ring frequency discrimination result exported between moment t, λ ifor the scale-up factor of i rank FLL loop bandwidth and characteristic frequency, Bf ifor the loop bandwidth that i rank FLL is current;
If sighting distance acceleration and sighting distance acceleration are all less than corresponding threshold value, the auxiliary coefficient α >1 that FLL assists phaselocked loop is set, 0< β <1;
If sighting distance acceleration is greater than corresponding threshold value and sighting distance acceleration is less than corresponding threshold value, the auxiliary coefficient 0< α <1 that FLL assists phaselocked loop is set, β >1;
If sighting distance acceleration and sighting distance acceleration are all greater than corresponding threshold value, the auxiliary coefficient 0< α <0.5 that FLL assists phaselocked loop is set, β >1.
According to embodiments of the invention, further, in above-mentioned steps S500, described loop filter comprises FLL, phaselocked loop and FLL and assists cycle of phase-locked loop ring wave filter and Loop filter, according to the loop bandwidth of new FLL, the loop bandwidth of phaselocked loop, and the product of identified result that auxiliary coefficient α and carrier wave ring phase detector export, the product of the frequency discrimination result that auxiliary coefficient β and carrier wave ring frequency discriminator export resets the bandwidth that FLL, phaselocked loop and FLL assist cycle of phase-locked loop ring wave filter.
According to embodiments of the invention, further, further, in above-mentioned steps S600, the Output rusults of cycle of phase-locked loop ring wave filter is assisted according to FLL, carrier generator adjustment produces the correlation parameter of signal, and according to filtered code ring identified result, the adjustment of ranging code generator produces the correlation parameter of signal.
Compared with prior art, the technical scheme that the present invention proposes has the following advantages:
1, dynamic perfromance and noise characteristic by detecting receiver in real time carry out the bandwidth of adjustment loop, make track loop can under the environment of dynamic change, especially the bandwidth of timely adjustment loop wave filter under the environment of high dynamic change, make it keep optimal value, tracking error is remained in less scope.
2, assist the auxiliary coefficient of phaselocked loop to increase the dynamically adapting scope of track loop further by adjustment FLL, thus enhance track loop under the environment of dynamic change, the stability especially under the environment of high dynamic change.
3, the present invention proposes technical scheme to overcome in prior art receiver tracking error under high dynamic environment large, the shortcoming of loss of lock easily occurs, and complexity is low, can realize easily by software.
Other features and advantages of the present invention will be set forth in the following description, and partly become apparent from instructions, or understand by implementing the present invention.Object of the present invention and other advantages realize by structure specifically noted in instructions, claims and accompanying drawing and obtain.
Accompanying drawing explanation
Accompanying drawing is used to provide a further understanding of the present invention, and forms a part for instructions, with embodiments of the invention jointly for explaining the present invention, is not construed as limiting the invention.In the accompanying drawings:
Fig. 1 is the structured flowchart of the track loop of the BD receiver adopted in the embodiment of the present invention;
Fig. 2 is the implementing procedure figure of the adaptive tracking method of BD receiver in the embodiment of the present invention.
Embodiment
Below in conjunction with drawings and Examples, technical scheme of the present invention is described in detail.
Fig. 1 shows the structured flowchart of the track loop of the BD receiver adopted in the embodiment of the present invention.Wherein, carrier wave NCO and code ring NCO represents carrier generator and ranging code generator respectively, is digital controlled oscillator.FLL assists PLL loop ring wave filter to be the abbreviation that FLL assists cycle of phase-locked loop ring wave filter.
This BD receiver, can by the control errors of tracking satellite navigation signal in less scope according to the process flow diagram work shown in Fig. 2.
The intermediate frequency carrier frequency that S100, basis are caught and ranging code phase place, initialization carrier generator and ranging code generator produce the correlation parameter of signal.Particular content is as follows.
After BD receiver to receive the carrier signal of corresponding frequency by antenna, through down coversion and analog to digital conversion, carrier signal is become digital intermediate frequency signal, be supplied to trapping module.Trapping module reproduction intermediate frequency carrier and ranging code, obtain intermediate frequency carrier frequency f carrwith the coarse value C of ranging code phase place pha.
According to the intermediate frequency carrier frequency values f obtained carrwith ranging code phase value C pha, initialization carrier generator and ranging code generator produce the correlation parameter of signal.In the present embodiment, these correlation parameters comprise the phase increment control word △ P of carrier generator f, and the initial phase C of ranging code generator pha0.Wherein:
&Delta; P f = 2 M &CenterDot; f carr f clk C pha0=C pha
Wherein, M is the word length of the phase accumulator of carrier generator, f clkfor the clock of carrier generator.
S200, the signal pair frequency digital signal utilizing carrier generator and ranging code generator to produce carry out carrier wave stripping and ranging code stripping, extract baseband digital signal.Particular content is as follows.
Carrier generator and ranging code generator produce corresponding signal respectively and carry out carrier wave stripping and ranging code stripping to the digital intermediate frequency signal obtained through analog to digital conversion, thus obtain following six road baseband digital signals:
Advanced code i e(n) q e(n) orthogonal two-way;
Instantaneous code i p(n) q p(n) orthogonal two-way;
Delayed code i l(n) q l(n) orthogonal two-way;
Wherein n=0,1,2 ... for sampled point sequence number.
S300, carries out phase demodulation, frequency discrimination and carrier-to-noise ratio to baseband digital signal and calculates.Particular content is as follows.
In embodiments of the present invention, preferably, also integration is carried out to above-mentioned six roadbed band signals and adds up, obtain following six road low speed base band datas:
I E ( k ) = &Sigma; n = kN ( k + 1 ) N - 1 i E ( n ) , Q E ( k ) = &Sigma; n = kN ( k + 1 ) N - 1 q E ( n ) ,
I P ( k ) = &Sigma; n = kN ( k + 1 ) N - 1 i P ( n ) , Q P ( k ) = &Sigma; n = kN ( k + 1 ) N - 1 q P ( n ) ,
I L ( k ) = &Sigma; n = kN ( k + 1 ) N - 1 i L ( n ) , Q L ( k ) = &Sigma; n = kN ( k + 1 ) N - 1 q L ( n ) ,
Wherein, k=0,1,2 ... for cumulative rear data point sequence number, T is integration accumulation interval, and N is original intermediate frequency data sampled point number in time T.
Then six road low speed base band datas are inputed to a yard ring phase detector, carrier wave ring phase detector, carrier wave ring frequency discriminator and carrier-to-noise ratio computing unit, carry out corresponding phase demodulation, frequency discrimination and carrier-to-noise ratio and calculate (as shown in Figure 1).Wherein, the identified result that code ring phase detector exports is designated as φ ek (), the identified result that carrier wave ring phase detector exports is designated as θ ek (), the frequency discrimination result that carrier wave ring frequency discriminator exports is designated as f ek (), the result of calculation that carrier-to-noise ratio computing unit exports is designated as CNR (k).
Should be understood that, in the present embodiment, carry out integration add up to baseband digital signal, then cumulative to integration result is carried out phase demodulation, frequency discrimination and carrier-to-noise ratio and is calculated, and is to reduce noise further.But this does not get rid of, directly baseband digital signal is carried out to the feasibility of phase demodulation, frequency discrimination and carrier-to-noise ratio calculating.
S400, according to phase demodulation, frequency discrimination and carrier-to-noise ratio result of calculation, ask for according to the mathe-matical map relation preset the auxiliary coefficient that new frequency-locked loop bandwidth sum phase-locked loop bandwidth and FLL assist phaselocked loop.Particular content is as follows.
In the present embodiment, preferred acquisition time interval T kthe result that interior carrier-to-noise ratio computing unit exports, gets the carrier-to-noise ratio of its average as the tracked satellite-signal of t
CNR &OverBar; t = 1 K &Sigma; k = 0 K CNR ( k )
Wherein, CNR (k) is a kth carrier-to-noise ratio result of calculation, and K is moment t-T kto the number of the result that carrier-to-noise ratio computing unit between moment t exports.
The loop bandwidth Bf of new (i.e. subsequent time) i rank FLL is calculated according to following formula i':
Bf i &prime; = &pi; 2 T 2 &CenterDot; CNR &OverBar; t &CenterDot; ( V ( i + 1 ) ) 2 &CenterDot; ( &lambda; i ) 2 i 9 F &CenterDot; ( 1 + 1 T &CenterDot; CNR &OverBar; t ) 2 i + 1
Wherein, i represents FLL exponent number, V (i+1)represent the i+1 order derivative of receiver bearer rate, λ ifor the scale-up factor of i rank FLL loop bandwidth and characteristic frequency, T is integration accumulation interval.F be one with the carrier-to-noise ratio of tracked satellite-signal relevant parameter, works as carrier-to-noise ratio when being greater than given threshold value, F gets 1, works as carrier-to-noise ratio when being less than given threshold value, F gets 2.
The loop bandwidth Bp of new (i.e. subsequent time) j rank phaselocked loop is calculated according to following formula j':
B p j &prime; = 4 CNR &OverBar; t &CenterDot; ( &mu; j ) 2 j &CenterDot; ( V ( j ) ) 2 9 ( 1 + 1 2 T &CenterDot; CNR &OverBar; t ) 2 j + 1
Wherein, j represents phaselocked loop exponent number, V (j)represent the j order derivative of receiver bearer rate, μ jfor the scale-up factor of j rank PLL loop bandwidth and characteristic frequency, T is integration accumulation interval.
Based on the frequency discrimination result f that frequency discriminator exports ek (), calculates the sighting distance dynamic value of receiver carrier t, also namely calculate sighting distance acceleration a twith sighting distance acceleration a t':
a t = 1 K &Sigma; k = 0 K f e ( k ) &CenterDot; ( B f i / &lambda; i ) , a t &prime; = 1 K &Sigma; k = 0 K f e ( k ) &CenterDot; ( B f i / &lambda; i ) 2
Wherein, f ek () is a kth frequency discrimination result, K is moment t-T kto the number exporting frequency discrimination result between moment t; Bf ifor the loop bandwidth of current time i rank FLL.
In the present embodiment, FLL assists auxiliary coefficient α and β of phaselocked loop according to sighting distance acceleration a twith sighting distance acceleration a t' size, reset according to following mechanism:
The threshold value supposed sighting distance given and accelerate with the threshold value of sighting distance acceleration
When time, auxiliary coefficient value α >1,0< β <1;
When time, auxiliary coefficient value 0< α <1, β >1;
When time, auxiliary coefficient value 0< α <0.5, β >1.
S500, the auxiliary coefficient of phaselocked loop is assisted to reset the bandwidth of loop filter according to new frequency-locked loop bandwidth sum phase-locked loop bandwidth and FLL.Particular content is as follows:
According to the loop bandwidth Bf of new FLL i', the loop bandwidth Bp of new phaselocked loop j', and the identified result θ that new auxiliary coefficient α and carrier wave ring phase detector export ethe product of (k), the frequency discrimination result f that new auxiliary coefficient β and carrier wave ring frequency discriminator export ek the product of () resets the bandwidth that FLL, phaselocked loop and FLL assist cycle of phase-locked loop ring wave filter.
S600, according to the Output rusults of loop filter adjustment carrier generator and ranging code generator, return step S200.
Particularly, as shown in Figure 1, assist the Output rusults of cycle of phase-locked loop ring wave filter according to FLL, carrier generator adjustment produces the correlation parameter of signal, and produces new signal accordingly, and the carrier wave carrying out a new round is peeled off; Simultaneously according to the Output rusults of Loop filter (also i.e. filtered identified result φ e(k)), the adjustment of ranging code generator produces the correlation parameter of signal, and produces new signal accordingly, and the ranging code carrying out a new round is peeled off.
Adjusted by this iterative cycles, the error of tracking satellite navigation signal remains in less scope by BD receiver, achieves the accurate measurement of carrier frequency to satellite and ranging code phase place.
The tracking that the present invention proposes, residing for detection receiver, the dynamic perfromance of environment and the noisiness of signal adjust the bandwidth of track loop loop filter and the auxiliary coefficient of FLL and phaselocked loop in real time, track loop is made to adjust inherent parameters in time under the environment of dynamic change, tracking error is remained in less scope, also increase the dynamically adapting scope of track loop simultaneously, improve the stability of track loop, solve the problem that existing BD receiver track loop error under high dynamic environment is large, losing lock easily occurs.In addition, the mathe-matical map relation set in above-described embodiment can by writing corresponding software simulating in the original central control unit of receiver, and without the need to increasing extra cost, implementation complexity is low.
Last it is noted that above each embodiment is only in order to illustrate technical scheme of the present invention, be not intended to limit; Although with reference to foregoing embodiments to invention has been detailed description, those of ordinary skill in the art is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein some or all of technical characteristic; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of various embodiments of the present invention technical scheme.

Claims (10)

1. an adaptive tracking method under receiver high dynamic environment, comprises the following steps:
The intermediate frequency carrier frequency that S100, basis are caught and ranging code phase place, initialization carrier generator and ranging code generator produce the correlation parameter of signal;
S200, the signal pair frequency carrier wave digital signal utilizing carrier generator and ranging code generator to produce carry out carrier wave stripping and ranging code stripping, extract baseband digital signal;
S300, baseband digital signal carried out to phase demodulation, frequency discrimination and carrier-to-noise ratio and calculate;
S400, according to phase demodulation, frequency discrimination and carrier-to-noise ratio result of calculation, ask for according to the mathe-matical map relation preset the auxiliary coefficient that new frequency-locked loop bandwidth sum phase-locked loop bandwidth and FLL assist phaselocked loop;
S500, the auxiliary coefficient of phaselocked loop is assisted to reset the bandwidth of loop filter according to new frequency-locked loop bandwidth sum phase-locked loop bandwidth and FLL;
S600, according to the Output rusults of loop filter adjustment carrier generator and ranging code generator, return step S200.
2. adaptive tracking method as claimed in claim 1, is characterized in that:
In step S100, described correlation parameter comprises the phase increment control word of carrier generator and the initial phase of ranging code generator;
The phase increment control word △ P of initialization carrier generator according to the following formula f:
&Delta; P f = 2 M &CenterDot; f carr f clk
Wherein, M is the word length of the phase accumulator of carrier generator, f clkfor the clock of carrier generator, f carrfor the intermediate frequency carrier frequency values of catching;
The initial phase C of initialization ranging code generator according to the following formula pha0:
C pha0=C pha
Wherein, C phafor the ranging code phase value of catching.
3. adaptive tracking method as claimed in claim 1, is characterized in that:
In step S300, described phase demodulation comprises yard ring phase demodulation and carrier wave ring phase demodulation, and described frequency discrimination comprises carrier wave ring frequency discrimination.
4. the adaptive tracking method as described in claims 1 to 3 any one, is characterized in that:
Further, in step S300, first integration is carried out to baseband digital signal and add up, and then the calculating of phase demodulation, frequency discrimination and carrier-to-noise ratio is carried out to integration accumulation result.
5. adaptive tracking method as claimed in claim 4, it is characterized in that, in step S400, described mathe-matical map relation comprises time interval T kthe average of interior carrier-to-noise ratio result of calculation is as the carrier-to-noise ratio of the tracked satellite-signal of t
CNR &OverBar; t = 1 K &Sigma; k = 0 K CNR ( k )
Wherein, CNR (k) is a kth carrier-to-noise ratio result of calculation, and K is moment t-T kto the number of carrier-to-noise ratio result of calculation between moment t.
6. adaptive tracking method as claimed in claim 4, it is characterized in that, in step S400, described mathe-matical map relation comprises the loop bandwidth Bf calculating new i rank FLL according to following formula i':
Bf i &prime; = &pi; 2 T 2 &CenterDot; CNR &OverBar; t &CenterDot; ( V ( i + 1 ) ) 2 &CenterDot; ( &lambda; i ) 2 i 9 F &CenterDot; ( 1 + 1 T &CenterDot; CNR &OverBar; t ) 2 i + 1
Wherein, i represents FLL exponent number, V (i+1)represent the i+1 order derivative of receiver bearer rate, λ ifor the scale-up factor of i rank FLL loop bandwidth and characteristic frequency, for the carrier-to-noise ratio of the tracked satellite-signal of t, T is integration accumulation interval, and F is one and carrier-to-noise ratio relevant parameter.
7. adaptive tracking method as claimed in claim 4, it is characterized in that, in step S400, described mathe-matical map relation comprises the loop bandwidth Bp calculating new j rank phaselocked loop according to following formula j':
B p j &prime; = 4 CNR &OverBar; t &CenterDot; ( &mu; j ) 2 j &CenterDot; ( V ( j ) ) 2 9 ( 1 + 1 2 T &CenterDot; CNR &OverBar; t ) 2 j + 1
Wherein, j represents phaselocked loop exponent number, V (j)represent the j order derivative of receiver bearer rate, μ jfor the scale-up factor of j rank PLL loop bandwidth and characteristic frequency, for the carrier-to-noise ratio of the tracked satellite-signal of t, T is integration accumulation interval.
8. adaptive tracking method as claimed in claim 4, it is characterized in that, in step S500, described mathe-matical map relation comprises:
Calculate the sighting distance acceleration a of receiver carrier t twith sighting distance acceleration a t';
a t = 1 K &Sigma; k = 0 K f e ( k ) &CenterDot; ( B f i / &lambda; i ) , a t &prime; = 1 K &Sigma; k = 0 K f e ( k ) &CenterDot; ( B f i / &lambda; i ) 2
Wherein, f ek () is a kth carrier wave ring frequency discrimination result, K is moment t-T kto the number of the carrier wave ring frequency discrimination result exported between moment t, λ ifor the scale-up factor of i rank FLL loop bandwidth and characteristic frequency, Bf ifor the loop bandwidth of current time i rank FLL;
If sighting distance acceleration and sighting distance acceleration are all less than corresponding threshold value, the auxiliary coefficient α >1 that FLL assists phaselocked loop is set, 0< β <1;
If sighting distance acceleration is greater than corresponding threshold value and sighting distance acceleration is less than corresponding threshold value, the auxiliary coefficient 0< α <1 that FLL assists phaselocked loop is set, β >1;
If sighting distance acceleration and sighting distance acceleration are all greater than corresponding threshold value, the auxiliary coefficient 0< α <0.5 that FLL assists phaselocked loop is set, β >1.
9. the adaptive tracking method as described in claims 1 to 3 any one, is characterized in that:
Further, in step S500, described loop filter comprises FLL, phaselocked loop and FLL and assists cycle of phase-locked loop ring wave filter and Loop filter, according to the loop bandwidth of new FLL, the loop bandwidth of phaselocked loop, and the product of identified result that auxiliary coefficient α and carrier wave ring phase detector export, the product of the frequency discrimination result that auxiliary coefficient β and carrier wave ring frequency discriminator export resets the bandwidth that FLL, phaselocked loop and FLL assist cycle of phase-locked loop ring wave filter.
10. the adaptive tracking method as described in claims 1 to 3 any one, is characterized in that:
Further, in step S600, assist the Output rusults of cycle of phase-locked loop ring wave filter according to FLL, carrier generator adjustment produces the correlation parameter of signal, and according to filtered code ring identified result, the adjustment of ranging code generator produces the correlation parameter of signal.
CN201410433179.XA 2014-08-28 2014-08-28 Adaptive tracking method under a kind of receiver high dynamic environment Active CN104215981B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410433179.XA CN104215981B (en) 2014-08-28 2014-08-28 Adaptive tracking method under a kind of receiver high dynamic environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410433179.XA CN104215981B (en) 2014-08-28 2014-08-28 Adaptive tracking method under a kind of receiver high dynamic environment

Publications (2)

Publication Number Publication Date
CN104215981A true CN104215981A (en) 2014-12-17
CN104215981B CN104215981B (en) 2017-03-15

Family

ID=52097677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410433179.XA Active CN104215981B (en) 2014-08-28 2014-08-28 Adaptive tracking method under a kind of receiver high dynamic environment

Country Status (1)

Country Link
CN (1) CN104215981B (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105699992A (en) * 2016-01-13 2016-06-22 上海交通大学 High-dynamic GNSS carrier wave signal tracking method and system thereof
CN105891854A (en) * 2016-04-01 2016-08-24 北京邮电大学 Dynamic threshold carrier ring
CN105891855A (en) * 2016-04-13 2016-08-24 南京理工大学 High dynamic GPS receiver carrier tracking method based on fuzzy control
CN105929419A (en) * 2016-04-13 2016-09-07 南京理工大学 BP artificial neural network-based GPS carrier tracking method
CN106019324A (en) * 2016-05-16 2016-10-12 河北科技大学 High-dynamic satellite navigation receiver loop tracker and tracking method thereof
CN106353780A (en) * 2016-08-31 2017-01-25 成都天奥信息科技有限公司 Method for high-dynamic adaption of satellite navigation receiver under condition of low signal-to-noise ratio
CN107070832A (en) * 2016-12-22 2017-08-18 中国电子科技集团公司第三十研究所 A kind of high dynamic UAV TT & C signal carrier tracking method
CN107305253A (en) * 2016-04-21 2017-10-31 大唐半导体设计有限公司 A kind of receiver tracking device and the method for realizing receiver tracking
CN108072883A (en) * 2017-12-24 2018-05-25 航天天绘科技有限公司 A kind of signal multimode of more spacecraft compatible interoperations receives locking processing method
CN108490416A (en) * 2018-05-21 2018-09-04 北京理工大学 A kind of high dynamic carrier phase tracking method of pseudo-code modulation continuous wave radar echo
CN108931792A (en) * 2017-05-22 2018-12-04 上海华测导航技术股份有限公司 A kind of method of loop parameter dynamic weighting
US10771030B2 (en) 2018-08-24 2020-09-08 Analog Devices International Unlimited Company Phase-locked loop with adjustable bandwidth
CN111812683A (en) * 2020-09-11 2020-10-23 四川九洲电器集团有限责任公司 Low-orbit satellite real-time tracking method, device, equipment and medium
CN113098497A (en) * 2021-03-23 2021-07-09 上海交通大学 All-digital phase-locked loop low-noise digital phase discriminator based on self-adaptive Kalman filter
CN113872569A (en) * 2021-12-06 2021-12-31 北京北斗华大科技有限公司 Smoothing filtering method and filter based on NWPR carrier-to-noise ratio algorithm

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107748376B (en) * 2017-09-26 2021-02-02 南京航空航天大学 Self-adaptive bandwidth tracking loop design method for satellite positioning receiver of unmanned vehicle

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6313789B1 (en) * 1998-06-10 2001-11-06 Topcon Positioning Systems, Inc. Joint tracking of the carrier phases of the signals received from different satellites
CN102621561A (en) * 2012-04-16 2012-08-01 南京航空航天大学 Loop self-regulation method of satellite navigation receiver
CN103163534A (en) * 2013-03-11 2013-06-19 西北工业大学 Adaptive noise bandwidth carrier loop tracking method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6313789B1 (en) * 1998-06-10 2001-11-06 Topcon Positioning Systems, Inc. Joint tracking of the carrier phases of the signals received from different satellites
CN102621561A (en) * 2012-04-16 2012-08-01 南京航空航天大学 Loop self-regulation method of satellite navigation receiver
CN103163534A (en) * 2013-03-11 2013-06-19 西北工业大学 Adaptive noise bandwidth carrier loop tracking method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陆煦: "跟踪环路在高动态GNSS接收机中的研究和应用", 《中国优秀硕士学位论文全文数据库 信息科技辑 》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105699992A (en) * 2016-01-13 2016-06-22 上海交通大学 High-dynamic GNSS carrier wave signal tracking method and system thereof
CN105699992B (en) * 2016-01-13 2017-12-15 上海交通大学 High-dynamic GNSS carriers signal trace method and system
CN105891854A (en) * 2016-04-01 2016-08-24 北京邮电大学 Dynamic threshold carrier ring
CN105891854B (en) * 2016-04-01 2018-03-30 北京邮电大学 A kind of dynamic threshold carrier wave ring
CN105929419B (en) * 2016-04-13 2019-01-18 南京理工大学 A kind of GPS carrier tracking based on BP artificial neural network
CN105891855A (en) * 2016-04-13 2016-08-24 南京理工大学 High dynamic GPS receiver carrier tracking method based on fuzzy control
CN105929419A (en) * 2016-04-13 2016-09-07 南京理工大学 BP artificial neural network-based GPS carrier tracking method
CN107305253A (en) * 2016-04-21 2017-10-31 大唐半导体设计有限公司 A kind of receiver tracking device and the method for realizing receiver tracking
CN106019324A (en) * 2016-05-16 2016-10-12 河北科技大学 High-dynamic satellite navigation receiver loop tracker and tracking method thereof
CN106353780A (en) * 2016-08-31 2017-01-25 成都天奥信息科技有限公司 Method for high-dynamic adaption of satellite navigation receiver under condition of low signal-to-noise ratio
CN107070832B (en) * 2016-12-22 2019-09-24 中国电子科技集团公司第三十研究所 A kind of high dynamic UAV TT & C signal carrier tracking method
CN107070832A (en) * 2016-12-22 2017-08-18 中国电子科技集团公司第三十研究所 A kind of high dynamic UAV TT & C signal carrier tracking method
CN108931792A (en) * 2017-05-22 2018-12-04 上海华测导航技术股份有限公司 A kind of method of loop parameter dynamic weighting
CN108072883A (en) * 2017-12-24 2018-05-25 航天天绘科技有限公司 A kind of signal multimode of more spacecraft compatible interoperations receives locking processing method
CN108490416A (en) * 2018-05-21 2018-09-04 北京理工大学 A kind of high dynamic carrier phase tracking method of pseudo-code modulation continuous wave radar echo
CN108490416B (en) * 2018-05-21 2019-04-02 北京理工大学 A kind of high dynamic carrier phase tracking method of pseudo-code modulation continuous wave radar echo
US10771030B2 (en) 2018-08-24 2020-09-08 Analog Devices International Unlimited Company Phase-locked loop with adjustable bandwidth
CN111812683A (en) * 2020-09-11 2020-10-23 四川九洲电器集团有限责任公司 Low-orbit satellite real-time tracking method, device, equipment and medium
CN111812683B (en) * 2020-09-11 2020-12-15 四川九洲电器集团有限责任公司 Low-orbit satellite real-time tracking method, device, equipment and medium
CN113098497A (en) * 2021-03-23 2021-07-09 上海交通大学 All-digital phase-locked loop low-noise digital phase discriminator based on self-adaptive Kalman filter
CN113872569A (en) * 2021-12-06 2021-12-31 北京北斗华大科技有限公司 Smoothing filtering method and filter based on NWPR carrier-to-noise ratio algorithm

Also Published As

Publication number Publication date
CN104215981B (en) 2017-03-15

Similar Documents

Publication Publication Date Title
CN104215981A (en) Adaptive tracking method for receiver in high-dynamic environment
CN107765270B (en) Satellite navigation receiver tracking loop based on Kalman filtering
CN104316941B (en) Vector tracking method based on carrier frequency assisted phase
CN102435999B (en) Baseband module of GPS (global positioning system) receiver and GPS signal acquiring and tracing method
CN111308520B (en) Tracking loop, method and receiver of satellite signal under high dynamic low signal-to-noise ratio
JPH02105746A (en) Differental phase transition keying modulator
CN102486539B (en) Improved navigation satellite signal tracking method
CN102621563A (en) GPS (global positioning system) software receiver signal tracking method based on FPGA (field programmable gate array) and system of GPS software receiver signal tracking method
CN103916201B (en) A kind of aerial signal initial phase difference, time delay and difference on the frequency estimation unit and method
CN104345323A (en) GPS satellite signal carrier loop tracking method and device
CN108199761B (en) USB transponder baseband digitization method
CN105388498A (en) Combined incoherent integral vector tracking method based on spatial domain
Weill A high performance code and carrier tracking architecture for ground-based mobile GNSS receivers
CN104765052B (en) GEO navigation satellite high-sensitivity carrier tracking method
US8094697B2 (en) Method and device for the demodulation of satellite radio navigation signals
CN109283557B (en) Double-pass pseudo code auxiliary carrier high-precision inter-satellite ranging system and method
Kozhaya et al. Blind receiver for LEO beacon estimation with application to UAV carrier phase differential navigation
CN103339526B (en) Equipment and method
Li et al. High dynamic carrier tracking using Kalman filter aided phase-lock loop
Shachi et al. Coherent BPSK demodulator using Costas loop and early-late gate synchronizer
CN104614742A (en) Beidou space-based high-precision differential information receiving implementation method
Patel et al. Faster methods for GPS signal acquisition in frequency domain
CN108169773A (en) A kind of satellite navigation signals tracking based on maximum likelihood coherent integration
Phyo et al. Implementation and analysis of signal tracking loops for software defined GPS receiver
CN115380224A (en) Positioning method, positioning chip and terminal equipment

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