[go: up one dir, main page]

CN104993844A - Method and device for searching frequency domain - Google Patents

Method and device for searching frequency domain Download PDF

Info

Publication number
CN104993844A
CN104993844A CN201510370038.2A CN201510370038A CN104993844A CN 104993844 A CN104993844 A CN 104993844A CN 201510370038 A CN201510370038 A CN 201510370038A CN 104993844 A CN104993844 A CN 104993844A
Authority
CN
China
Prior art keywords
frequency
frequency domain
sub
stage
coherent integration
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
CN201510370038.2A
Other languages
Chinese (zh)
Other versions
CN104993844B (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.)
UNICORE COMMUNICATIONS (BEIJING) Inc
Original Assignee
UNICORE COMMUNICATIONS (BEIJING) Inc
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 UNICORE COMMUNICATIONS (BEIJING) Inc filed Critical UNICORE COMMUNICATIONS (BEIJING) Inc
Priority to CN201510370038.2A priority Critical patent/CN104993844B/en
Publication of CN104993844A publication Critical patent/CN104993844A/en
Application granted granted Critical
Publication of CN104993844B publication Critical patent/CN104993844B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Complex Calculations (AREA)

Abstract

The invention provides a method and a device for searching a frequency domain. The method comprises the steps as follows: performing first-order down-conversion to a sampling point sequence according to a frequency domain search range; segmenting a coherent integration time according to a first preset value so as to obtain sub coherent integration times; calculating a first-order correlation value, corresponding to each sub coherent integration time, of a first result,; performing second-order down-conversion to each first order correlation value; performing second-order coherent integration calculation to each obtained second result; performing discrete Fourier transform DFT to each obtained second-order correlation value; and determining the frequency of a capturing moment according to each frequency value obtained after performing the DFT transformation. According to the technical solution, the first-order serial search is finished by performing the first-order down-conversion to the sampling point sequence according to the frequency domain search range, and the second-order serial search is finished by performing the second-order down-conversion, and the third-order parallel search is finished by using the DFT, so that the method and the device of the invention could be used for accurately and fast searching the frequency domain in wide range.

Description

Frequency domain searching method and device
Technical Field
The present invention relates to satellite navigation technologies, and in particular, to a frequency domain searching method and apparatus based on a spread spectrum receiver.
Background
In a direct sequence spread spectrum communication system, the acquisition of a spreading code by a receiver is a precondition for realizing signal demodulation. Due to the relative motion of the transmitting and receiving parties, a doppler frequency offset exists between the received signal and the local signal. Theoretical research shows that, during time domain search, the existence of doppler frequency offset affects the correlation detection amplitude, as shown in formula (1):
<math> <mrow> <msub> <mi>Z</mi> <mi>d</mi> </msub> <mo>=</mo> <mi>Z</mi> <mo>&times;</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mo>|</mo> <mfrac> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>&pi;Lf</mi> <mi>d</mi> </msub> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>&pi;f</mi> <mi>d</mi> </msub> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>|</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein Z is the correlation detection amplitude without frequency deviation, fdIs Doppler frequency offset (hereinafter referred to as frequency offset), ZdAt a frequency offset of fdAmplitude of time-dependent detection, L being the spreading ratio, TcIs the duration of each chip. As can be seen from equation (1), in the presence of frequency offset, the correlation detection amplitude decays.
When the frequency offset is small, the attenuation of the relevant detection amplitude value can be ignored, and the receiver can complete the acquisition of the spread spectrum code by pure time domain search. However, a large frequency offset causes a significant loss of the time-domain correlated signal-to-noise ratio, and therefore, two-dimensional search in the frequency domain and the time domain is required to acquire the spreading code.
The time domain searching is to align the phase of the local spreading code generator with the phase of the spreading code in the received signal to realize the capture of the spreading code; the frequency domain searching is to estimate the doppler frequency offset and compensate the received signal, so as to reduce the signal-to-noise ratio loss caused by the frequency offset of the received signal and ensure the success of the time domain searching.
In a Global Navigation Satellite System (GNSS), a Global Positioning System (GPS) and a beidou Satellite System use a code division multiple access technology to distinguish different satellites by different pseudo random codes, and a GLONASS (GLONASS) System uses a frequency division multiple access technology to distinguish different satellites by different carrier frequencies. In practice, these systems spread the modulated data using direct sequence spreading. In order to receive the navigation data of a certain satellite, pseudo-random codes for modulating the data must be reproduced, and the reproduced pseudo-codes and the input pseudo-codes are subjected to correlation operation on different phase errors to synchronize the reproduced pseudo-codes and the input pseudo-codes, so that the despreading of the navigation data is completed, and the process is called pseudo-code acquisition. On the other hand, due to the radial movement between the GNSS satellite moving at high speed and the receiver, the carrier wave generates the doppler shift effect. Therefore, after analog down-conversion (radio frequency module) and digital down-conversion of the satellite signal, the frequency value is not zero, but a doppler shift is added on the basis of zero intermediate frequency. In order to demodulate navigation data of a certain satellite, the value of the doppler shift generated by the corresponding satellite must be searched, and the process is called carrier acquisition or frequency domain search.
The existing frequency domain searching method has two different methods of serial and parallel. The advantage of serial search is that the algorithm is simple to implement, has low requirements for system resources, but has a slow capture speed, and is not suitable for application in a long-period code and high-dynamic environment. The parallel search method generally uses a Fast Fourier Transform (FFT) parallel search method, which has a fast acquisition speed, but the search frequency range is limited, which results in a higher signal-to-noise ratio loss due to an excessively large frequency interval.
Disclosure of Invention
In order to solve the technical problem, the invention provides a frequency domain searching method and a frequency domain searching device, which can realize large-range, accurate and quick searching of a frequency domain.
In order to achieve the object of the present invention, the present invention provides a frequency domain searching method, comprising:
performing first-stage down-conversion on the sampling point sequence according to the frequency domain searching range;
segmenting the coherent integration time according to a first preset value to obtain sub-coherent integration time; respectively calculating first-stage correlation values of the first results corresponding to the sub-coherent integration time;
performing second-stage down-conversion on each calculated first-stage correlation value; performing second-stage coherent integration operation on each obtained second result;
performing Discrete Fourier Transform (DFT) on each obtained second-stage correlation value; and determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion.
Further, the performing a first stage down-conversion on the sampling point sequence according to the frequency domain search range includes:
determining a first center frequency of the frequency domain searching range according to the frequency domain searching range;
performing first-stage down-conversion on the sampling point sequence according to the determined first central frequency to obtain a first result;
wherein the first center frequency is a center frequency value corresponding to the frequency domain search range.
Further, the performing of the second-stage down-conversion on each calculated first-stage correlation value includes:
segmenting the determined frequency domain searching range according to a first preset value;
respectively determining second center frequencies of the sub-frequency domain search ranges obtained after segmentation;
performing second-stage down-conversion on each first-stage correlation value according to each determined second center frequency to obtain a second result;
wherein each second center frequency is a center frequency value corresponding to each first sub-frequency domain search range.
Further, the performing a second-stage coherent integration operation on the obtained second result includes:
and respectively carrying out second-stage coherent integration operation on each obtained second result according to a second sub-coherent integration time to obtain each second-stage correlation value, wherein the second sub-coherent integration time is sub-coherent integration time L, and L is a positive integer.
Further, the performing DFT on each obtained second-stage correlation value includes:
segmenting each sub-frequency domain searching range according to a second preset value; to obtain each second sub-frequency domain search range;
determining a third center frequency of each obtained second sub-frequency domain searching range;
performing DFT on each obtained second-stage correlation value according to each determined third center frequency;
wherein each third center frequency is a center frequency value corresponding to the each second sub-frequency domain search range.
Further, the determining the frequency at the capturing time according to the obtained frequency values after DFT conversion includes:
judging each frequency value obtained after DFT conversion and a judgment threshold;
and when the frequency value after DFT conversion is larger than the judgment threshold, determining the frequency value after DFT conversion as the frequency of the capture moment.
The invention also provides a frequency domain searching method, which is characterized by comprising the following steps: the device comprises a first module, a calculation module, a second module and a determination module; wherein,
the first module is used for carrying out first-stage down-conversion on the sampling point sequence according to the frequency domain searching range;
the calculation module is used for segmenting the coherent integration time according to a first preset value so as to obtain sub-coherent integration time; respectively calculating first-stage correlation values of the first results corresponding to the sub-coherent integration time;
the second module is used for carrying out second-stage down-conversion on each calculated first-stage correlation value; performing second-stage coherent integration operation on each obtained second result;
the determining module is used for performing Discrete Fourier Transform (DFT) on each obtained second-stage correlation value; and determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion.
Further, the first module is specifically configured to:
determining a first center frequency of the frequency domain searching range according to the frequency domain searching range;
performing first-stage down-conversion on the sampling point sequence according to the determined first central frequency to obtain a first result;
wherein the first center frequency is a center frequency value corresponding to the frequency domain search range.
Further, the second module is specifically configured to:
segmenting the determined frequency domain searching range according to a first preset value;
respectively determining second center frequencies of the sub-frequency domain search ranges obtained after segmentation;
performing second-stage down-conversion on each first-stage correlation value according to each determined second center frequency to obtain a second result;
wherein each second center frequency is a center frequency value corresponding to each first sub-frequency domain search range.
Further, the second module is further specifically configured to:
and respectively carrying out second-stage coherent integration operation on each obtained second result according to a second sub-coherent integration time to obtain each second-stage correlation value, wherein the second sub-coherent integration time is sub-coherent integration time L, and L is a positive integer.
Further, the determining module is specifically configured to:
segmenting each sub-frequency domain searching range according to a second preset value; to obtain each second sub-frequency domain search range;
determining a third center frequency of each obtained second sub-frequency domain searching range;
performing DFT on each obtained second-stage correlation value according to each determined third center frequency;
wherein each third center frequency is a center frequency value corresponding to the each second sub-frequency domain search range.
Further, the determining module is further specifically configured to:
judging each frequency value obtained after DFT conversion and a judgment threshold;
and when the frequency value after DFT conversion is larger than the judgment threshold, determining the frequency value after DFT conversion as the frequency of the capture moment.
The technical scheme of the invention comprises the following steps: performing first-stage down-conversion on the sampling point sequence according to the frequency domain searching range; segmenting the coherent integration time according to a first preset value to obtain sub-coherent integration time; respectively calculating first-stage correlation values of the first results corresponding to the sub-coherent integration time; performing second-stage down-conversion on each calculated first-stage correlation value; performing second-stage coherent integration operation on each obtained second result; performing Discrete Fourier Transform (DFT) on each obtained second-stage correlation value; and determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion. In the technical scheme of the invention, the first-stage down-conversion is carried out on the sampling point sequence through the frequency domain searching range to complete the first-stage serial search, the second-stage serial search is completed through the second-stage down-conversion, and the third-stage parallel search is completed through DFT, so that the large-range, accurate and rapid search of the frequency domain is realized.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the invention without limiting the invention. In the drawings:
FIG. 1 is a flow chart of a frequency domain searching method of the present invention;
fig. 2 is a schematic structural diagram of the frequency domain searching apparatus according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be noted that the embodiments and features of the embodiments in the present application may be arbitrarily combined with each other without conflict.
Fig. 1 is a flowchart of a frequency domain searching method of the present invention, as shown in fig. 1, including the following steps:
step 101: and carrying out first-stage down-conversion on the sampling point sequence according to the frequency domain searching range. The method specifically comprises the following steps:
determining a first center frequency of the frequency domain searching range according to the frequency domain searching range;
performing first-stage down-conversion on the sampling point sequence according to the determined first central frequency to obtain a first result;
wherein the first center frequency is a center frequency value corresponding to the frequency domain search range.
For example, when the search range is ± 5KHz, the first center frequency is a center frequency value of 0KHz of ± 5KHz, and when the search range is 0KHz to 20KHz, the first center frequency is 10 KHz. The following calculation methods of the second center frequency and the third center frequency are the same as the calculation method of this step, and will not be described again.
The step also includes determining the frequency domain search range before the step, and the specific implementation is well known to those skilled in the art, and is not used to limit the protection scope of the present invention, and will not be described herein again.
It should be noted that how to perform the down-conversion is a conventional technical means known to those skilled in the art, and is not intended to limit the present invention, and will not be described herein.
Step 102: segmenting the coherent integration time according to a first preset value to obtain sub-coherent integration time; first-stage correlation values of the first results corresponding to the respective sub-coherent integration times are respectively calculated.
The first result is the result of performing the first stage down-conversion on the sampling point sequence according to the first center frequency.
For example, the coherent integration time is known as TcohAssuming that each sub-coherent integration time to be acquired is T1Then can be divided intoSub-coherent integration time, the first preset value in step 102 may be M, where M is a positive integer greater than or equal to 1, and M first-stage correlation values may be obtained in this step and may be recorded as x0,x1,..,xM-1
It should be noted that how to perform the coherent integration operation on the first result to obtain the first-order correlation value belongs to the conventional technical means known to those skilled in the art, and is not described herein again.
Step 103: performing second-stage down-conversion on each calculated first-stage correlation value; and carrying out second-stage coherent integration operation on each obtained second result.
And the second result is the result of carrying out second-stage down-conversion on each calculated first-stage correlation value.
Performing a second-stage down-conversion on each calculated first-stage correlation value comprises:
segmenting the determined frequency domain searching range according to a first preset value;
respectively determining second center frequencies of the sub-frequency domain search ranges obtained after segmentation;
performing second-stage down-conversion on each first-stage correlation value according to each determined second center frequency to obtain a second result;
wherein each second center frequency is a center frequency value corresponding to each first sub-frequency domain search range.
And performing a second stage of coherent integration operation on the obtained second result comprises:
and respectively carrying out second-stage coherent integration operation on each obtained second result according to a second sub-coherent integration time to obtain each second-stage correlation value, wherein the second sub-coherent integration time is sub-coherent integration time L, and L is a positive integer.
For example, the second result can be obtained by using equation (2):
<math> <mrow> <msub> <mover> <mi>x</mi> <mo>~</mo> </mover> <mi>m</mi> </msub> <mo>=</mo> <msub> <mi>x</mi> <mi>m</mi> </msub> <mo>&times;</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mn>2</mn> <mi>&pi;</mi> <mi>j</mi> <mo>&times;</mo> <msub> <mi>f</mi> <mn>2</mn> </msub> <mo>&times;</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>&times;</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein M is 0, 1, M-1, f2Is the second center frequency, T1For the sub-coherent integration time, j is an imaginary unit, e.g., in complex number a + bi, a is the real part and bi is the imaginary part. And performing second-stage coherent integration operation on the second result according to the second sub-coherent integration time (i.e. accumulating each L points into one point) to obtain N-M/L second-stage correlation values which are recorded as y0,y1,...,yN-1
Step 104: performing Discrete Fourier Transform (DFT) on each obtained second-stage correlation value; and determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion.
Performing DFT on each obtained second-level correlation value includes:
segmenting each sub-frequency domain searching range according to a second preset value; to obtain each second sub-frequency domain search range;
determining a third center frequency of each obtained second sub-frequency domain searching range;
performing DFT on each obtained second-stage correlation value according to each determined third center frequency;
wherein each third center frequency is a center frequency value corresponding to each second sub-frequency domain search range.
And the number of the first and second groups,
determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion, comprising:
judging each frequency value obtained after DFT conversion and a judgment threshold;
and when the frequency value after DFT conversion is larger than the judgment threshold, determining the frequency value after DFT conversion as the frequency of the capture moment.
By way of example, following the above example, y is plotted according to equation (3)0,y1,...,yN-1The DFT is performed and the result is obtained,
<math> <mrow> <msub> <mi>Y</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>y</mi> <mi>n</mi> </msub> <mo>&times;</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mn>2</mn> <mi>&pi;</mi> <mi>j</mi> <mo>&times;</mo> <msub> <mi>f</mi> <mn>3</mn> </msub> <mo>&times;</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>&times;</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein N is 0, 1, N-1, k is 0, 1, N-1f3Is the third center frequency, T2Is the sub-coherent integration time, the second sub-coherent integration time, T2=T1XL, j is the unit of imaginary number.
In the method, the first-stage down-conversion is carried out on the sampling point sequence through the frequency domain searching range to complete the first-stage serial search, the second-stage serial search is completed through the second-stage down-conversion, and the third-stage parallel search is completed through DFT, so that the wide-range, accurate and quick search of the frequency domain is realized.
Fig. 2 is a schematic structural diagram of the frequency domain searching apparatus of the present invention, as shown in fig. 2, including: the device comprises a first module, a calculation module, a second module and a determination module. Wherein,
and the first module is used for carrying out first-stage down-conversion on the sampling point sequence according to the frequency domain searching range.
Wherein, the first module is specifically configured to:
determining a first center frequency of the frequency domain searching range according to the frequency domain searching range;
performing first-stage down-conversion on the sampling point sequence according to the determined first central frequency to obtain a first result;
wherein the first center frequency is a center frequency value corresponding to the frequency domain search range.
The calculation module is used for segmenting the coherent integration time according to a first preset value so as to obtain sub-coherent integration time; first-stage correlation values of the first results corresponding to the respective sub-coherent integration times are respectively calculated.
The first result is the result of performing the first stage down-conversion on the sampling point sequence according to the first center frequency.
The second module is used for carrying out second-stage down-conversion on each calculated first-stage correlation value; and carrying out second-stage coherent integration operation on the obtained second result.
Wherein, the second module is specifically configured to:
segmenting the determined frequency domain searching range according to a first preset value;
respectively determining second center frequencies of the sub-frequency domain search ranges obtained after segmentation;
performing second-stage down-conversion on each first-stage correlation value according to each determined second center frequency to obtain a second result;
wherein each second center frequency is a center frequency value corresponding to each first sub-frequency domain search range.
The second module is further specifically configured to:
and respectively carrying out second-stage coherent integration operation on each obtained second result according to a second sub-coherent integration time to obtain each second-stage correlation value, wherein the second sub-coherent integration time is sub-coherent integration time L, and L is a positive integer.
A determining module, configured to perform a Discrete Fourier Transform (DFT) on each obtained second-level correlation value; and determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion.
The determining module is specifically configured to:
segmenting each sub-frequency domain searching range according to a second preset value; to obtain each second sub-frequency domain search range;
determining a third center frequency of each obtained second sub-frequency domain searching range;
performing DFT on each obtained second-stage correlation value according to each determined third center frequency;
wherein each third center frequency is a center frequency value corresponding to each second sub-frequency domain search range.
And the determining module is further specifically configured to:
judging each frequency value obtained after DFT conversion and a judgment threshold;
and when the frequency value after DFT conversion is larger than the judgment threshold, determining the frequency value after DFT conversion as the frequency of the capture moment.
It will be understood by those skilled in the art that all or part of the steps of the above methods may be implemented by instructing the relevant hardware through a program, and the program may be stored in a computer readable storage medium, such as a read-only memory, a magnetic or optical disk, and the like. Alternatively, all or part of the steps of the above embodiments may be implemented using one or more integrated circuits. Accordingly, each module/unit in the above embodiments may be implemented in the form of hardware, and may also be implemented in the form of a software functional module. The present application is not limited to any specific form of hardware or software combination.
The above description is only a preferred example of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (12)

1. A frequency domain searching method, comprising:
performing first-stage down-conversion on the sampling point sequence according to the frequency domain searching range;
segmenting the coherent integration time according to a first preset value to obtain sub-coherent integration time; respectively calculating first-stage correlation values of the first results corresponding to the sub-coherent integration time;
performing second-stage down-conversion on each calculated first-stage correlation value; performing second-stage coherent integration operation on each obtained second result;
performing Discrete Fourier Transform (DFT) on each obtained second-stage correlation value; and determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion.
2. The method of claim 1, wherein the first stage down-converting the sequence of sample points according to the frequency domain search range comprises:
determining a first center frequency of the frequency domain searching range according to the frequency domain searching range;
performing first-stage down-conversion on the sampling point sequence according to the determined first central frequency to obtain a first result;
wherein the first center frequency is a center frequency value corresponding to the frequency domain search range.
3. The method of claim 1, wherein said second down-converting each calculated first-stage correlation value comprises:
segmenting the determined frequency domain searching range according to a first preset value;
respectively determining second center frequencies of the sub-frequency domain search ranges obtained after segmentation;
performing second-stage down-conversion on each first-stage correlation value according to each determined second center frequency to obtain a second result;
wherein each second center frequency is a center frequency value corresponding to each first sub-frequency domain search range.
4. The method of claim 1, wherein performing a second stage coherent integration operation on the obtained second result comprises:
and respectively carrying out second-stage coherent integration operation on each obtained second result according to a second sub-coherent integration time to obtain each second-stage correlation value, wherein the second sub-coherent integration time is sub-coherent integration time L, and L is a positive integer.
5. The method of claim 1, wherein performing a DFT on each obtained second-stage correlation value comprises:
segmenting each sub-frequency domain searching range according to a second preset value; to obtain each second sub-frequency domain search range;
determining a third center frequency of each obtained second sub-frequency domain searching range;
performing DFT on each obtained second-stage correlation value according to each determined third center frequency;
wherein each third center frequency is a center frequency value corresponding to the each second sub-frequency domain search range.
6. The method according to any one of claims 1 to 5, wherein determining the frequency at the time of capture from the obtained DFT-transformed individual frequency values comprises:
judging each frequency value obtained after DFT conversion and a judgment threshold;
and when the frequency value after DFT conversion is larger than the judgment threshold, determining the frequency value after DFT conversion as the frequency of the capture moment.
7. A frequency domain searching apparatus method, comprising: the device comprises a first module, a calculation module, a second module and a determination module; wherein,
the first module is used for carrying out first-stage down-conversion on the sampling point sequence according to the frequency domain searching range;
the calculation module is used for segmenting the coherent integration time according to a first preset value so as to obtain sub-coherent integration time; respectively calculating first-stage correlation values of the first results corresponding to the sub-coherent integration time;
the second module is used for carrying out second-stage down-conversion on each calculated first-stage correlation value; performing second-stage coherent integration operation on each obtained second result;
the determining module is used for performing Discrete Fourier Transform (DFT) on each obtained second-stage correlation value; and determining the frequency of the capturing moment according to the obtained frequency values after DFT conversion.
8. The apparatus of claim 7, wherein the first module is specifically configured to:
determining a first center frequency of the frequency domain searching range according to the frequency domain searching range;
performing first-stage down-conversion on the sampling point sequence according to the determined first central frequency to obtain a first result;
wherein the first center frequency is a center frequency value corresponding to the frequency domain search range.
9. The apparatus of claim 7, wherein the second module is specifically configured to:
segmenting the determined frequency domain searching range according to a first preset value;
respectively determining second center frequencies of the sub-frequency domain search ranges obtained after segmentation;
performing second-stage down-conversion on each first-stage correlation value according to each determined second center frequency to obtain a second result;
wherein each second center frequency is a center frequency value corresponding to each first sub-frequency domain search range.
10. The apparatus of claim 9, wherein the second module is further specifically configured to:
and respectively carrying out second-stage coherent integration operation on each obtained second result according to a second sub-coherent integration time to obtain each second-stage correlation value, wherein the second sub-coherent integration time is sub-coherent integration time L, and L is a positive integer.
11. The apparatus of claim 7, wherein the determining module is specifically configured to:
segmenting each sub-frequency domain searching range according to a second preset value; to obtain each second sub-frequency domain search range;
determining a third center frequency of each obtained second sub-frequency domain searching range;
performing DFT on each obtained second-stage correlation value according to each determined third center frequency;
wherein each third center frequency is a center frequency value corresponding to the each second sub-frequency domain search range.
12. The apparatus according to any one of claims 7 to 11, wherein the determining module is further specifically configured to:
judging each frequency value obtained after DFT conversion and a judgment threshold;
and when the frequency value after DFT conversion is larger than the judgment threshold, determining the frequency value after DFT conversion as the frequency of the capture moment.
CN201510370038.2A 2015-06-29 2015-06-29 A kind of frequency domain search method and device Active CN104993844B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510370038.2A CN104993844B (en) 2015-06-29 2015-06-29 A kind of frequency domain search method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510370038.2A CN104993844B (en) 2015-06-29 2015-06-29 A kind of frequency domain search method and device

Publications (2)

Publication Number Publication Date
CN104993844A true CN104993844A (en) 2015-10-21
CN104993844B CN104993844B (en) 2017-10-13

Family

ID=54305607

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510370038.2A Active CN104993844B (en) 2015-06-29 2015-06-29 A kind of frequency domain search method and device

Country Status (1)

Country Link
CN (1) CN104993844B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106998586A (en) * 2017-06-01 2017-08-01 电子科技大学 The synchronization acquiring method of wireless communication system in a kind of high dynamic environment
CN113820731A (en) * 2021-11-24 2021-12-21 湖南跨线桥航天科技有限公司 Method and system for capturing non-periodic long code signal of satellite navigation receiver

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101056139A (en) * 2007-05-24 2007-10-17 重庆重邮信科(集团)股份有限公司 Method for realizing the downlink rough synchronization in the initial search process of TD-SCDMA system cell
CN101174852A (en) * 2006-11-02 2008-05-07 中兴通讯股份有限公司 A method and system for WCDMA cell search time slot synchronization
CN101286753A (en) * 2007-04-11 2008-10-15 中兴通讯股份有限公司 A preamble sending and receiving method of a broadband wireless communication system
CN101447823A (en) * 2007-11-27 2009-06-03 杰脉通信技术(上海)有限公司 TD-SCDMA system down link synchronisation method
CN101741800A (en) * 2009-10-14 2010-06-16 苏州恩巨网络有限公司 Synchronous searching method
CN101789807A (en) * 2008-07-18 2010-07-28 俊茂微电子(上海)有限公司 Frequency deviation evaluation method, device and communication device
CN102162852A (en) * 2010-12-07 2011-08-24 清华大学 Method and system for capturing weak GNSS (Global Navigation Satellite System) signal under condition of large-scale frequency deviation
CN102315859A (en) * 2011-09-23 2012-01-11 华为技术有限公司 Multi-path search method and device
CN102422664A (en) * 2009-05-05 2012-04-18 意法爱立信有限公司 PPS (pulse Per second) detection method of LTE (Long term evolution) communication network
CN102655419A (en) * 2012-05-08 2012-09-05 中国人民解放军国防科学技术大学 Calculation method of decision function for captured spread spectrum signals
CN102655420A (en) * 2011-03-02 2012-09-05 中兴通讯股份有限公司 Subsection frequency domain method and device for realizing multi-diameter search
CN103616703A (en) * 2013-11-28 2014-03-05 中国航天科工信息技术研究院 Base band data processing method and system of navigation receiver
CN104601512A (en) * 2015-01-05 2015-05-06 江苏指南针导航通信技术有限公司 Method and system for detecting carrier frequency offset of phase-modulated signals

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101174852A (en) * 2006-11-02 2008-05-07 中兴通讯股份有限公司 A method and system for WCDMA cell search time slot synchronization
CN101286753A (en) * 2007-04-11 2008-10-15 中兴通讯股份有限公司 A preamble sending and receiving method of a broadband wireless communication system
CN101056139A (en) * 2007-05-24 2007-10-17 重庆重邮信科(集团)股份有限公司 Method for realizing the downlink rough synchronization in the initial search process of TD-SCDMA system cell
CN101447823A (en) * 2007-11-27 2009-06-03 杰脉通信技术(上海)有限公司 TD-SCDMA system down link synchronisation method
CN101789807A (en) * 2008-07-18 2010-07-28 俊茂微电子(上海)有限公司 Frequency deviation evaluation method, device and communication device
CN102422664A (en) * 2009-05-05 2012-04-18 意法爱立信有限公司 PPS (pulse Per second) detection method of LTE (Long term evolution) communication network
CN101741800A (en) * 2009-10-14 2010-06-16 苏州恩巨网络有限公司 Synchronous searching method
CN102162852A (en) * 2010-12-07 2011-08-24 清华大学 Method and system for capturing weak GNSS (Global Navigation Satellite System) signal under condition of large-scale frequency deviation
CN102655420A (en) * 2011-03-02 2012-09-05 中兴通讯股份有限公司 Subsection frequency domain method and device for realizing multi-diameter search
CN102315859A (en) * 2011-09-23 2012-01-11 华为技术有限公司 Multi-path search method and device
CN102655419A (en) * 2012-05-08 2012-09-05 中国人民解放军国防科学技术大学 Calculation method of decision function for captured spread spectrum signals
CN103616703A (en) * 2013-11-28 2014-03-05 中国航天科工信息技术研究院 Base band data processing method and system of navigation receiver
CN104601512A (en) * 2015-01-05 2015-05-06 江苏指南针导航通信技术有限公司 Method and system for detecting carrier frequency offset of phase-modulated signals

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106998586A (en) * 2017-06-01 2017-08-01 电子科技大学 The synchronization acquiring method of wireless communication system in a kind of high dynamic environment
CN113820731A (en) * 2021-11-24 2021-12-21 湖南跨线桥航天科技有限公司 Method and system for capturing non-periodic long code signal of satellite navigation receiver
CN113820731B (en) * 2021-11-24 2022-06-10 湖南跨线桥航天科技有限公司 Method and system for capturing non-periodic long code signal of satellite navigation receiver

Also Published As

Publication number Publication date
CN104993844B (en) 2017-10-13

Similar Documents

Publication Publication Date Title
US8462616B2 (en) Apparatus and method for estimating a frequency shift and a time shift
CN104536016B (en) A new GNSS system signal acquisition device and method
US10042056B2 (en) Signal processing method for ultra-fast acquisition and tracking of severely attenuated spread spectrum signals with doppler frequency and apparatus thereof
CN104765052B (en) GEO navigation satellite high-sensitivity carrier tracking method
CN100553169C (en) Short-period spreading code time frequency parallel searching method based on frequency folding
CN108196274A (en) Be applicable in BOC (n, n) signal without fuzziness catching method and device
CN104765050A (en) A New Beidou Signal Second Acquisition Algorithm
Tamazin et al. Robust fine acquisition algorithm for GPS receiver with limited resources
Kim et al. Design of FFT-based TDCC for GNSS acquisition
CN103760578B (en) A kind of GNSS satellite navigation signal without fuzzy tracking method
CN104993844B (en) A kind of frequency domain search method and device
CN109581434B (en) B2a signal capturing method and device
US20110050496A1 (en) Energy Domain Based Peak Reconstruction Methods And Apparatuses
CN104765053B (en) The catching method and device of pseudo-code in a kind of GNSS receiver
CN104422939B (en) Method and apparatus for being associated with the signal received from navigational satellite system
Rui et al. Blind Estimation of GPS M‐Code Signals under Noncooperative Conditions
CN111708058B (en) Signal capturing method and device
CN112764063A (en) Method for realizing capture processing and receiver
CN103339526A (en) Apparatus and method
CN107132552B (en) Parallel code phase searching device and method for realizing parallel code phase searching
Rakipi et al. GPS signal acquisition and sensitivity analysis using different algorithms on a software receiver
CN107132555B (en) A kind of parallel code phase search device and the method for realizing parallel code phase search
Beldjilali et al. A new proposed GPS satellite signals acquisition algorithm based on the fast Fourier transform
Zhao et al. A novel PN-code acquisition method based on local frequency folding for BeiDou system
Su et al. An unfuzzy acquisition algorithm based on matched filtering for BOC (n, n)

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant