CN101420411B - Fast carrier capture method with low signal-noise ratio - Google Patents

Fast carrier capture method with low signal-noise ratio Download PDF

Info

Publication number
CN101420411B
CN101420411B CN2008102389625A CN200810238962A CN101420411B CN 101420411 B CN101420411 B CN 101420411B CN 2008102389625 A CN2008102389625 A CN 2008102389625A CN 200810238962 A CN200810238962 A CN 200810238962A CN 101420411 B CN101420411 B CN 101420411B
Authority
CN
China
Prior art keywords
frequency
nco
data
sampled
module
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
CN2008102389625A
Other languages
Chinese (zh)
Other versions
CN101420411A (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.)
Space Star Technology Co Ltd
Original Assignee
Space Star Technology 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 Space Star Technology Co Ltd filed Critical Space Star Technology Co Ltd
Priority to CN2008102389625A priority Critical patent/CN101420411B/en
Publication of CN101420411A publication Critical patent/CN101420411A/en
Application granted granted Critical
Publication of CN101420411B publication Critical patent/CN101420411B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Noise Elimination (AREA)

Abstract

A method for quickly catching the low signal-to-noise ratio carrier comprises the following steps: (1) multiplying an input signal with a cos branching data generated by a first NCO module; (2) executing half band filtering; (3) executing down sampling to a first accumulated down sampling module, temporarily storing into an exterior SDRAM through an FIFO-IN module, and the data speed fsl and the maximum capture bandwidth after the first accumulated down sampling satisfy the relationship: fsl>2*f, [f1, f2] is the capture range, (4) reading the data in the exterior SDRAM through an FIFO-OUT module, and multiplying data of two branches of sin and cos generated by a second NCO; (5) entering a second down-sampling through the LPF; (6) executing FFT procession; and (7) adding a second NCO localoscillating frequency, repeating the steps (4)-(6) for n times to find out a maximum value in the results of step (6) and obtain the final frequency offset result. The method according to the invention has the advantages of adaptability for low signal-to-noise ratio, wide capture range and quick capture on the base that the modulation degree characteristic of the universal controlling signal is satisfied.

Description

Fast carrier capture method with low signal-noise ratio
Technical field
The present invention relates to a kind of low signal-to-noise ratio high-frequency offset carrier quick capturing method, be mainly used in the fast carrier capture under the standard USB observing and controlling system, solved the problem of fast carrier capture under the lower big offset frequency situation of signal to noise ratio.
Background technology
Unified S frequency range (USB) observing and controlling system is a kind of observing and controlling system of maturation, and it can be finished angle tracking simultaneously, tests the speed, range finding, remote measurement, remote control and function of voice communication.Along with the continuous development of digital communication system, digitlization capturing carrier technology becomes a common technology wherein.For the estimation and the rectification of demodulating process medium frequency side-play amount, existing ripe algorithm, as phase-locked loop (PLL) method, methods such as automatic frequency control (AFC) method and Frequency Estimation.Owing to contain Doppler frequency shift in the intermediate-freuqncy signal, PLL, AFC or some are applicable to that high speed Frequency Estimation method can not finish the task that frequency deviation value is caught fast under low noise, and the condition that also can't adapt to the capture range broad, this opens up new method with regard to needs and solves this problem.On the frequency spectrum of unified USB measurement and control signal, the amplitude of residual carrier is maximum in modulation degree when not being very big, therefore can utilize this specific character to extract the frequency values of carrier wave; Simultaneously the frequency spectrum of measurement and control signal is symmetrically distributed on frequency domain, also can estimate centre frequency according to this symmetrical structure, both frequency values of carrier wave, but the common speed of this method is slower, and the memory space of data is bigger, need take big hardware resource.Article " the method research of USB answering machine fast carrier capture ", Zhang Chunlong etc., 2003 national emulation technology academic meeting paper collection, characteristic according to the frequency symmetry, method with FFT-DFT combination is judged the frequency of carrier wave, the resolution of this method is low, and capture range is narrow, the characteristic that can not satisfy wide capture range simultaneously and catch fast.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiencies in the prior art, a kind of fast carrier capture method with low signal-noise ratio is provided, this method has the advantages that to take into account low signal-to-noise ratio, wide capture range, catch fast on the basis of satisfying general measurement and control signal modulation degree characteristic.
Technical solution of the present invention is: fast carrier capture method with low signal-noise ratio comprises the following steps:
(1) coS that produces of an input signal and a NCO (digital controlled oscillator) props up circuit-switched data and multiplies each other, and wherein the local frequency of a NCO is f NCO1, the frequency values of input signal is f 0+ Δ f, f 0Be the centre frequency of input signal, Δ f is the frequency deviation value of input signal;
(2) step (1) multiplied result is carried out semi-band filtering and handle the filtering high fdrequency component;
(3) result of step (2) is sent into the first down-sampled module that adds up and carry out down-sampled processing, the data after the processing are temporarily stored among the outside SDRAM by the FIFO_IN module, and the storage number of data points is greater than N among the SDRAM FFT* L, wherein L = f s 1 f s 2 Be first down-sampled the counting of adding up, f S1Be first frequency after the down-sampled data that adds up, f S2For being the first down-sampled back frequency that adds up, the first data rate f that adds up after the down-sampled processing S1And should satisfy between the maximum acquisition bandwidth f: f S12 * f, [f 1, f 2] be capture range, f 1Be negative direction minimum frequency value, f 2Be the maximum acquisition bandwidth of positive direction, f=f 2-f 1
(4) by FIFO_OUT the data among the outside SDRAM are read in, and multiply each other the local frequency f of the 2nd NCO with the sin of the 2nd NCO generation and the data of two branch roads of cos NCO2Variable, the local frequency f of the 2nd NCO NCO2Initial value be 0;
(5) step (4) result is passed through LPF (low pass filter) respectively, the filtering high fdrequency component enters the second down-sampled module that adds up then and carries out the down-sampled processing second time, and down-sampled purpose is to improve signal to noise ratio for the second time, the data rate f after for the second time down-sampled S2Local carrier frequency N with the 2nd NCO CO2The step value Δ between relation should satisfy: f S2=2 * Δ;
(6) result of step (5) is carried out N FFTThe FFT of point handles, wherein N FFTBe counting that FFT handles; N wherein FFTValue determined the precision of the frequency deviation value that calculates, N FFTBig more, precision is high more, but while N FFTCan not be too big, because N FFTBig more, the resource that program takies is big more, calculates I simultaneously 2+ Q 2Maximum, and record maximum pairing frequency deviation value f m, wherein the I road is a real part, the Q road is an imaginary part, m ∈ [1,2 ... .., n];
(7) the 2nd NCO local frequency f in the increase step (4) NCO2, repeating step (4)~(6) n time, wherein n = int [ f 2 - f 1 Δ + 1 ] , Find out the maximum f among n repeating step (6) result i, i ∈ [1,2 ... .., n], final frequency deviation value result is:
Figure G2008102389625D00032
The principle of the inventive method: Mathematical Modeling is as follows, k wherein 1, k 2, k 3Be constant, ω cBe the centre frequency of input signal, Δ ω is a frequency offseting value:
If the expression formula of input signal is:
a(t)=k 1sin(ω ct+Δω)    (1)
The signal that the one NCO produces is:
b(t)=k 2cos(ω ct) (2)
The result was respectively after input signal and local carrier multiplied each other:
a ( t ) * ( t ) = k 1 sin ( ω c t + Δω ) × k 2 cos ( ω c t )
(3)
= 1 2 k 1 k 2 [ sin ( 2 ω c t + Δω ) + sin ( Δω ) ]
Through behind the half-band filter, the high fdrequency component in the filtered signal is output as:
s ( t ) = 1 2 k 1 k 2 sin ( Δω ) - - - ( 4 )
The two-way orthogonal signalling that the 2nd NCO module produces are respectively, wherein
Figure G2008102389625D00036
Being the local carrier frequency of the 2nd NCO, is variable, and the each data that produce with the 2nd NCO of signal multiply each other, the local carrier of the 2nd NCO step value Δ that adds up:
c I(t)=k 3cos(φ) (5)
c Q(t)=k 3sin(φ) (6)
Signal through first add up down-sampled after, deposit among the outside SDRAM by FIFO_IN, the two-way orthogonal data that data of reading in from FIFO_OUT and the 2nd NCO produce multiplies each other, and obtains:
Wherein, A = 1 2 k 1 k 2 k 3 .
Through LPF with the high fdrequency component filtering, carry out then second add up down-sampled after, be equivalent to pass through a comb filter, improved signal to noise ratio in the image frequency component the filtering band outside, at this moment data are output as:
Figure G2008102389625D00042
Figure G2008102389625D00043
Give FFT module as I, Q two paths of data respectively with (9), (10), carry out FFT and handle.
The present invention's beneficial effect compared with prior art is:
(1) the present invention has solved the big data quantity storage problem that in the past occurred because of resource-constrained in the FPGA sheet by using the method for outside SDRAM, has improved acquisition speed simultaneously;
(2) the present invention adopts the two-stage down-sampled processing that adds up, and adopts FFT Frequency Estimation method, when reducing storage and counting, has improved signal to noise ratio, can catch accurately under than the low signal-to-noise ratio situation;
(3) the present invention adopts two-stage NCO processing mode, the 2nd NCO local carrier frequency is variable, has increased the frequency acquisition scope, and has solved the problem that needs repeatedly to store data in the frequency sweep facture in the past, only needs storage one secondary data to get final product among the present invention.
(4) the present invention is mainly used in the capturing carrier under the standard USB observing and controlling system, has solved the problem of fast carrier capture under than low signal-to-noise ratio, big offset frequency situation.This method is through actual test, prove its have capture time short, catch advantages such as precision height, capture range be wide.Its capture range can reach-300KHz~+ 300KHz, capture time is less than 1S, C/N 0The minimum 35dB that reaches.
Description of drawings
Fig. 1 is the flow chart of fast carrier capture method with low signal-noise ratio of the present invention.
Fig. 2 is the operation principle block diagram of two fifo modules of the present invention.
Embodiment
With integrated test high speed remote sensing data receiving and processing equipment and many systems detection and control terminal is that example illustrates specific implementation process of the present invention, as shown in Figure 1,
(1) in this method, the hardware system processing clock is 110MHz, and capture range [300,300] KHz ,-300KHz are negative direction minimum frequency value, and 300KHz is the positive direction maximum frequency values, and the local frequency step value Δ of the 2nd NCO is 2.5KHz;
(2) frequency values is f 0The Cos that the input signal of+Δ f and a NCO module produce props up circuit-switched data and multiplies each other, and wherein the local frequency of a NCO is f NCO1Be 39.6975MHz, f 0=40MHz is the centre frequency of input signal, and Δ f=300KHz is the frequency deviation value of input signal;
(3) step (2) multiplied result is carried out semi-band filtering and handle the filtering high fdrequency component;
(4) result of step (3) is sent into the first down-sampled module that adds up and carry out down-sampled processing, the first data rate f that adds up after the down-sampled processing S1And should satisfy between the maximum frequency deviation value f: f S12 * f, so f=f 2-f 1=600KHz considers that down-sampled the counting that add up should be chosen for positive integer, therefore gets f S1=1.718750MHz, down-sampled the counting of promptly adding up for the first time is the 110MHz/1.71875MHz=64 point;
(5) data after step (4) processing are temporarily stored among the outside SDRAM by the FIFO_IN module, the storage number of data points is greater than N FFT* L, wherein L = int [ f s 1 f s 2 + 0.5 ] = 344 Be down-sampled the counting of adding up for the second time, in actual applications, the storage number of data points is 360000 points, and the storage needed time of data is about 0.21s;
(6) by the FIFO_OUT module data among the outside SDRAM are read in, and multiply each other the local frequency f of the 2nd NCO with the sin of the 2nd NCO generation and the data of two branch roads of cos NCO2Variable, the local frequency f of the 2nd NCO NCO2Initial value be 0;
(7) step (6) result is passed through LPF respectively, the filtering high fdrequency component enters the second down-sampled module that adds up then and carries out the down-sampled processing second time, and down-sampled purpose is to improve signal to noise ratio for the second time, the data rate f after for the second time down-sampled S2Local carrier frequency f with the 2nd NCO NCO2The step value Δ between the pass be: f S2=2 * Δ=5KHz, second down-sampled count into L = int [ f s 1 f s 2 + 0.5 ] = 344 , This moment, data clock was reduced to 320KHz by 110M Hz, and data rate is reduced to 5KHz by 1.71875MHz, and clock sampling 1024 point data of 320KHz need the time to be about 3.2ms;
(8) result of step (7) is carried out N FFT=1024 FFT handles, wherein N FFTBe counting that FFT handles; N wherein FFTValue determined the precision of the frequency deviation value that calculates, N FFTBig more, precision is high more, but while N FFTCan not be too big, because N FFTBig more, the resource that program takies is big more, and in sum, choosing counting of FFT in the practical application is 1024 points, and the resolution behind the FFT can reach 4.88Hz;
(9) I, Q two paths of data after the FFT processing are calculated I 2+ Q 2Maximum, and record maximum pairing frequency deviation value f m, wherein the I road is a real part, the Q road is an imaginary part, m ∈ [1,2 ... .., n], it is 29.976us that each FFT handles the needed time;
(10) the 2nd NCO local frequency f in the increase step (6) NCO2, repeating step (6)~(9) n time, wherein n = int [ f 2 - f 1 Δ + 1 ] = 241 , Find out the maximum f among 241 repeating steps (9) result i, i ∈ [1,2 ... .., 241], final frequency deviation value result is:
Δ f=(f 0-f NCO1)-(i * Δ+f i)=299.99997KHz, with the actual value deviation be 0.03Hz;
(11) needed capture time is: 0.21s+241* (29.976us+3.2ms)=0.988s, promptly capture time is less than 1s.
As shown in Figure 2, to be the method that adopts time-sharing multiplex become a big Capacity FIFO Memory with the SDRAM emulation of FPGA outside for FIFO_IN, FIFO_OUT module functions.The 1Kx16bit FIFO that the input and output port of FIFO has FPGA BIock RAM to constitute respectively.FIFO_IN is responsible for the data of the outside input of buffer memory, and gives inner sdram controller with data passes.FIFO_OUT is responsible for the data that the buffer memory sdram controller provides, and these data can be read away by external logic from output port.Sdram controller is responsible for monitoring state and control and the read-write SDRAM of two FIFO.When among the FIFO_IN data being arranged, sdram controller is about to these data and reads and write among the SDRAM, has clearance spaces to be in preserving new data and FIFO_OUT among the SDRAM, and FIFO_OUT is read and write to sdram controller with the data among the SDRAM.Like this, under the time-sharing work of sdram controller, it is a total FIFO that the input and output port of FIFO is modeled into.As long as input and output data transfer rate summation is no more than the interface data rate of SDRAM, can keep this FIFO steady operation.
The frequency deviation value that table 1 obtains for this method Practical Calculation and the table of comparisons of ideal value.
Frequency deviation value (MHz) C/N0(dB) Measured value (KHz) With theoretical value deviation (Hz)
39.88407 35 186.572269 2.269
39.99031 34 292.812512 2.512
39.99031 36 292.812512 2.512
40.30 35 602.49997 0.03
40.0218 35 324.301744 1.744
39.70 35 2.499977 0.023
39.70 36 2.499977 0.023
40.00 35 302.50485 4.85
In sum, it is as follows to reach a conclusion: method for capturing carrier of the present invention can correctly be caught when C/N0 〉=35dB.
The part that the present invention does not elaborate belongs to techniques well known.

Claims (2)

1. fast carrier capture method with low signal-noise ratio is characterized in that step is as follows:
(1) cos that produces of an input signal and a NCO props up circuit-switched data and multiplies each other, and wherein the local frequency of a NCO is f NCO1, the frequency values of input signal is f 0+ Δ f, f 0Be the centre frequency of input signal, Δ f is the frequency deviation value of input signal;
(2) step (1) multiplied result is carried out semi-band filtering and handle the filtering high fdrequency component;
(3) result of step (2) is sent into the first down-sampled module that adds up and carry out down-sampled processing, the data after the processing are temporarily stored among the outside SDRAM by the FIFO_IN module, the first data rate f that adds up after the down-sampled processing S1And should satisfy between the maximum acquisition bandwidth f: f S1>2 * f, [f 1, f 2] be capture range, f 1Be negative direction minimum frequency value, f 2Be the positive direction maximum frequency values, maximum acquisition bandwidth f=f 2-f 1
(4) by the FIFO_OUT module data among the outside SDRAM are read in, and multiply each other the local frequency f of the 2nd NCO with the sin of the 2nd NCO generation and the data of two branch roads of cos NCO2Variable, the local frequency f of the 2nd NCO NCO2Initial value be 0;
(5) step (4) result is passed through LPF respectively, the filtering high fdrequency component enters the second down-sampled module that adds up then and carries out the down-sampled processing second time, the data rate f after for the second time down-sampled S2Local carrier frequency f with the 2nd NCO NCO2The step value Δ between relation should satisfy: f S2=2 * Δ;
(6) result of step (5) is carried out N FFTThe FFT of point handles, wherein N FFTBe counting that FFT handles;
(7) the 2nd NCO local frequency f in the increase step (4) NCO2, repeating step (4)~(6) n time, wherein
Figure FSB00000207856400011
Find out the maximum fi among n repeating step (6) result, i ∈ [1,2 ... .., n], final frequency deviation value result is:
Figure FSB00000207856400012
2. fast carrier capture method with low signal-noise ratio according to claim 1 is characterized in that: the FIFO_IN module in the described step (3) is temporarily stored into number of data points among the outside SDRAM greater than N FFT* L, wherein
Figure FSB00000207856400021
Be down-sampled the counting of adding up for the second time, f S1Be the first down-sampled back data rate that adds up, f S2It is the second down-sampled back data rate that adds up.
CN2008102389625A 2008-12-05 2008-12-05 Fast carrier capture method with low signal-noise ratio Active CN101420411B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008102389625A CN101420411B (en) 2008-12-05 2008-12-05 Fast carrier capture method with low signal-noise ratio

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008102389625A CN101420411B (en) 2008-12-05 2008-12-05 Fast carrier capture method with low signal-noise ratio

Publications (2)

Publication Number Publication Date
CN101420411A CN101420411A (en) 2009-04-29
CN101420411B true CN101420411B (en) 2011-02-09

Family

ID=40631021

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008102389625A Active CN101420411B (en) 2008-12-05 2008-12-05 Fast carrier capture method with low signal-noise ratio

Country Status (1)

Country Link
CN (1) CN101420411B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101741424B (en) * 2009-12-24 2013-08-28 航天恒星科技有限公司 Method for rapidly capturing multi-mode high dynamic spread spectrum signal
CN109379314B (en) * 2018-12-10 2022-02-08 北京卫星信息工程研究所 High speed burst digital demodulation method and apparatus

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1018827A1 (en) * 1999-01-08 2000-07-12 Sony International (Europe) GmbH Synchronisation structure for OFDM system
WO2002076056A2 (en) * 2001-03-16 2002-09-26 Rohde & Schwarz Gmbh & Co. Kg Method for the frequency and time synchronization of an ofdm receiver
CN1859345A (en) * 2005-04-30 2006-11-08 中国科学院上海微系统与信息技术研究所 Quick frequency catching method based on OFDM system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1018827A1 (en) * 1999-01-08 2000-07-12 Sony International (Europe) GmbH Synchronisation structure for OFDM system
WO2002076056A2 (en) * 2001-03-16 2002-09-26 Rohde & Schwarz Gmbh & Co. Kg Method for the frequency and time synchronization of an ofdm receiver
CN1859345A (en) * 2005-04-30 2006-11-08 中国科学院上海微系统与信息技术研究所 Quick frequency catching method based on OFDM system

Also Published As

Publication number Publication date
CN101420411A (en) 2009-04-29

Similar Documents

Publication Publication Date Title
CN109495410B (en) High dynamic PCM/FM signal carrier frequency accurate estimation method
CN101726746B (en) Intermediate frequency direct sequence spread spectrum receiver for satellite ranging
CN102520423B (en) Rapid capturing method for Circuit for rapidly capturing long period pseudo random spread spectrum code of satellite navigation receiver and capturing method thereof
CN108055058B (en) High-precision measurement method for carrier Doppler and change rate thereof
CN201607527U (en) Intermediate-frequency direct-sequence spread spectrum receiver
CN102508200B (en) Plural spread spectrum signal TDOA (time difference of arrival) measuring method
CN102540217A (en) System for rapidly acquiring GLONASS satellite signals in field programmable gate array (FPGA)
CN103837740A (en) High-precision digital instantaneous frequency measurement method and device
CN110113278B (en) Modulation mode identification method based on all-digital receiver
CN104777495A (en) QPSK (quadrature phase shift keying) modulated I/Q branch orthogonality test method based on distribution histogram
CN105187348A (en) Any-rate CPFSK (Continuous Phase Frequency Shift Key) signal timing synchronization method
CN102508269B (en) Satellite navigation pilot signal acquisition method, pseudo random sequence stripping method and device
CN105429918B (en) A kind of quick capturing method for deep space TT&C signal
CN101420411B (en) Fast carrier capture method with low signal-noise ratio
CN110336583A (en) A kind of quick catching device and method of Direct Sequence Spread Spectrum Signal
CN102495280B (en) Anti-noise wide-range frequency measurement method and phase locking frequency meter
CN111796306A (en) Navigation satellite signal receiving method and receiver
CN107222226A (en) The method and apparatus of harmonic carcellation interference
CN101174849B (en) Spread-spectrum code chip synchronous catching and tracing method and device of wireless sensing net node
CN106253945A (en) A kind of spread-spectrum signal broad sense carrier synchronization system and method for super large dynamic environment
CN109067676B (en) High-precision time domain performance evaluation method for satellite navigation signals
CN106453188A (en) Rapid and accurate frequency synchronization method applicable for MPSK demodulation
CN105717519A (en) Satellite navigation signal channel SCB characteristic detecting method
CN110018500A (en) A kind of Big Dipper satellite signal capture method based on circular shifting
CN102624418B (en) Hydroacoustic biphase modulation direct sequence spread spectrum signal carrier frequency estimation 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
C14 Grant of patent or utility model
GR01 Patent grant