CN113162652B - Method and module for detecting frequency of spread spectrum signal - Google Patents

Method and module for detecting frequency of spread spectrum signal Download PDF

Info

Publication number
CN113162652B
CN113162652B CN202110471627.5A CN202110471627A CN113162652B CN 113162652 B CN113162652 B CN 113162652B CN 202110471627 A CN202110471627 A CN 202110471627A CN 113162652 B CN113162652 B CN 113162652B
Authority
CN
China
Prior art keywords
signal
frequency
spread spectrum
detected
carrier 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.)
Active
Application number
CN202110471627.5A
Other languages
Chinese (zh)
Other versions
CN113162652A (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.)
SHAANXI LINGYUN ELECTRONICS GROUP CO LTD
Original Assignee
SHAANXI LINGYUN ELECTRONICS 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 SHAANXI LINGYUN ELECTRONICS GROUP CO LTD filed Critical SHAANXI LINGYUN ELECTRONICS GROUP CO LTD
Priority to CN202110471627.5A priority Critical patent/CN113162652B/en
Publication of CN113162652A publication Critical patent/CN113162652A/en
Application granted granted Critical
Publication of CN113162652B publication Critical patent/CN113162652B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/709Correlator structure
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/20Monitoring; Testing of receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/309Measuring or estimating channel quality parameters
    • H04B17/336Signal-to-interference ratio [SIR] or carrier-to-interference ratio [CIR]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • H04L2027/0024Carrier regulation at the receiver end

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Quality & Reliability (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

The invention discloses a method for detecting frequency of spread spectrum signals, which comprises the following steps: receiving a spread spectrum signal to be detected, wherein the spread spectrum signal to be detected comprises modulation information; carrying out square operation on the spread spectrum signal to be detected to eliminate the modulation signal and obtain a frequency estimation signal; sequentially carrying out down-conversion and frequency mixing on the frequency estimation signal, and carrying out narrow-band filtering on the frequency-mixed signal; and performing DFT operation on the filtered signal function to detect the carrier frequency deviation. The frequency domain processing is directly carried out on the signal correlation value by the spread spectrum signal frequency detection method, the modulation phase ambiguity is eliminated, a frequency discriminator is not needed, only a DFT algorithm is needed to carry out frequency spectrum analysis, and the frequency discrimination range is greatly expanded. In addition, the DFT algorithm only needs to use little resource consumption, and the frequency discrimination resolution can be greatly improved.

Description

Method and module for detecting frequency of spread spectrum signal
Technical Field
The invention belongs to the field of signal processing, and particularly relates to a method and a module for detecting frequency of a spread spectrum signal.
Background
In the conventional spread spectrum communication system and satellite navigation receiving application, due to the influence of doppler frequency shift and local carrier error, a received signal may generate slow frequency drift, so that the phase drift of the received signal occurs, which may have a bad influence on a communication system using coherent demodulation, and thus the receiving performance may be degraded. In an all-digital receiver, a local oscillation signal of a digital down-conversion (DDC) is limited, and digital system resources, working frequency and the like cannot slide in a larger range, so that the tracking bandwidth of a loop cannot be large for optimizing the performance and designing a larger carrier loop equivalent bandwidth.
Thus, a corresponding estimation and compensation of the carrier frequency is required before accurate phase tracking. Conventional frequency estimation algorithms fall into two categories, open-loop estimation algorithms and closed-loop estimation algorithms. The traditional closed-loop estimation algorithm is simple to implement, low in calculation requirement and short in acquisition time, but the tracking range is small, the acquisition performance can be close to the best at high signal-to-noise ratio, and at low signal-to-noise ratio, a symbol error caused by noise is transmitted to a feedback loop through a decision device, so that the loop performance is rapidly deteriorated. As a first estimation parameter for demodulation, such performance is clearly not satisfactory.
Disclosure of Invention
In order to solve the above problems in the prior art, the present invention provides a method and a module for detecting a frequency of a spread spectrum signal. The technical problem to be solved by the invention is realized by the following technical scheme:
a method of frequency detection of a spread spectrum signal, comprising:
receiving a spread spectrum signal to be detected, wherein the spread spectrum signal to be detected comprises modulation information;
carrying out square operation on the spread spectrum signal to be detected to eliminate the modulation signal and obtain a frequency estimation signal;
sequentially carrying out down-conversion and frequency mixing on the frequency estimation signal, and carrying out narrow-band filtering on the frequency-mixed signal;
the DFT operation is performed on the filtered signal function to detect the carrier frequency offset.
In a specific embodiment, the expression of the spread spectrum signal to be detected is:
r(t)=A*P(t)*e j(2π(f+Δf)t+θ) +n(t),
wherein A is signal power, P (t) is shaping pulse shape, f is carrier frequency, Δ f is carrier frequency offset, and θ is modulation information phase of current symbol,
Figure BDA0003045557560000021
n (t) is noise.
In one embodiment, the frequency estimation signal expression is:
r 2 '(t)=A 2 *e j(2π(2*f+2*Δf)t) +n 2 (t),
wherein, A is signal power, f is carrier frequency, Δ f is carrier frequency offset, and n (t) is noise.
In one embodiment, the frequency of the mixed signal is: f2=2 x F2-fs/2, where F2 is the frequency of the signal of interest before mixing and fs is the system frequency.
The invention also provides a spread spectrum signal frequency detection module, which can be independently arranged and can be integrated into a system to be used as a functional module of the system, and the module comprises:
the system comprises a signal receiving unit, a signal processing unit and a signal processing unit, wherein the signal receiving unit is used for receiving a spread spectrum signal to be detected, and the spread spectrum signal to be detected comprises modulation information;
the signal modulation unit is used for carrying out square operation on the spread spectrum signal to be detected so as to eliminate the modulation signal and obtain a frequency estimation signal;
the filtering unit is used for sequentially carrying out down-conversion and frequency mixing on the frequency estimation signal and carrying out narrow-band filtering on the frequency-mixed signal;
and the detection unit is used for performing DFT operation on the filtered signal function so as to detect the carrier frequency deviation.
In a specific embodiment, the expression of the spread spectrum signal to be detected is:
r(t)=A*P(t)*e j(2π(f+Δf)t+θ) +n(t),
wherein A is signal power, P (t) is shaping pulse shape, f is carrier frequency, Δ f is carrier frequency offset, and θ is modulation information phase of current symbol,
Figure BDA0003045557560000031
n (t) is noise.
In one embodiment, the frequency estimation signal expression is:
r 2 '(t)=A 2 *e j(2π(2*f+2*Δf)t) +n 2 (t),
wherein, A is signal power, f is carrier frequency, Δ f is carrier frequency offset, and n (t) is noise.
In one embodiment, the frequency of the mixed signal is: f2=2 x F2-fs/2, where F2 is the frequency of the signal of interest before mixing and fs is the system frequency.
The invention has the beneficial effects that:
the frequency domain processing is directly carried out on the signal correlation value by the spread spectrum signal frequency detection method, the modulation phase ambiguity is eliminated, a frequency discriminator is not needed, only a DFT algorithm is needed to carry out frequency spectrum analysis, and the frequency discrimination range is greatly expanded. In addition, the DFT algorithm only needs to use little resource consumption, and the frequency discrimination resolution can be greatly improved.
The present invention will be described in further detail with reference to the accompanying drawings and examples.
Drawings
Fig. 1 is a schematic flowchart of a method for detecting a frequency of a spread spectrum signal according to an embodiment of the present invention;
fig. 2 is an output when the signal-to-noise ratio, Δ f = +7.5KHz in a simulation example of the embodiment of the present invention;
FIG. 3 is a graph of the output at high SNR, Δ f = -7.5KHz in a simulation example of an embodiment of the present invention;
FIG. 4 shows the output when the signal-to-noise ratio = -1dB, and Δ f = +7.5KHz in a simulation example of the embodiment of the present invention;
FIG. 5 shows the output when the signal-to-noise ratio = -1dB, and Δ f = -7.5KHz in a simulation example of the embodiment of the present invention;
fig. 6 is a block diagram of a spread spectrum signal frequency detection module according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to specific examples, but the embodiments of the present invention are not limited thereto.
Example one
Referring to fig. 1, fig. 1 is a schematic flow chart of a method for detecting a frequency of a spread spectrum signal according to an embodiment of the present invention, including:
receiving a spread spectrum signal to be detected, wherein the spread spectrum signal to be detected comprises modulation information;
carrying out square operation on the spread spectrum signal to be detected to eliminate the modulation signal and obtain a frequency estimation signal;
sequentially carrying out down-conversion and frequency mixing on the frequency estimation signal, and carrying out narrow-band filtering on the frequency-mixed signal;
the DFT operation is performed on the filtered signal function to detect the carrier frequency offset.
In a specific embodiment, the expression of the spread spectrum signal to be detected is:
r(t)=A*P(t)*e j(2π(f+Δf)t+θ) +n(t),
wherein A is signal power, P (t) is shaping pulse shape, f is carrier frequency, Δ f is carrier frequency offset, and θ is modulation information phase of current symbol,
Figure BDA0003045557560000041
n (t) is noise.
To more clearly illustrate the solution of the present embodiment, the following exemplary derivation is given:
assuming that the shaping pulse shape is an ideal rectangle, P (t) =1, and the above equation becomes:
r(t)=A*e j(2π(f+Δf)t+θ) +n(t)
squaring r (t) yields:
r 2 (t)=A 2 *e j(2π(2*f+2*Δf)t+2*θ) +n 2 (t)
since θ ∈ (0, π), then 2 ∈ (0, 2 π), from the periodicity, the above equation can be transformed into:
r 2 '(t)=A 2 *e j(2π(2*f+2*Δf)t) +n 2 (t)
the above formula shows that: the integration of the square operation to eliminate the correlation value eliminates the modulation information, and the output signal contains 2 frequency multiplication components [2 x f,2 x delta f ] of the carrier and the frequency deviation, thereby converting the carrier frequency deviation estimation problem of the modulation signal into the frequency estimation of the dot frequency component.
Further, the frequency estimation signal expression is:
r 2 '(t)=A 2 *e j(2π(2*f+2*Δf)t) +n 2 (t),
wherein, A is signal power, f is carrier frequency, Δ f is carrier frequency offset, and n (t) is noise.
In one embodiment, the frequency of the mixed signal is: f2=2 × F2-fs/2, where F2 is the signal frequency of interest before mixing and fs is the system frequency.
The smaller the sampling rate, the smaller the original signal center frequency and bandwidth that can be recovered without aliasing after sampling. Therefore, the frequency estimation signal needs to be down-converted to move it to the vicinity of the baseband (low frequency). The lower the frequency after down-conversion is, the smaller the sampling frequency can be used, the finer the frequency resolution is, or the same resolution is maintained, the smaller the number of points of DFT (N is smaller) can be used, thereby reducing the operation amount of DFT and saving precious hardware resources.
The digital domain mixing and filtering needs a large amount of multiplication and addition operations, particularly, a multiplier is a precious and limited hardware resource, and the fs/2 frequency conversion without multiplication is a mixing method without multiplication to realize the frequency conversion of fixed frequency points.
For a set of discrete sequences:
A(n)=e jπn =cos(πn)=(-1) n
a (n) includes (1, -1), A (n) = (1, -1, \8230;). When multiplying a (n) as a mixing sequence by a time-domain signal sequence, which is equivalent to multiplying this signal sequence by a cosine sampling sequence cos (n), the frequency of a (n) is fs/2 because its cosine value repeats every two points. Although the (-1) n sequence is multiplied by the input signal sequence, in a digital signal processor (such as an FPGA), the multiplication by 1 or-1 is only simple to keep the original value or to invert (the complementary code is inverted and added by one), and multiplication is not actually used, so that the shifting of the input signal in the frequency shift is simply completed.
After mixing, the present application does not perform spectrum analysis by using the existing FFT method, but performs analysis by using DFT (discrete fourier transform).
Specifically, for a finite sequence x (N) of length N, its discrete fourier transform is:
Figure BDA0003045557560000061
in the formula, W N =e -j*2π/N
If the formula is directly used for spectrum analysis, N times of complex multiplication and N (N-1) times of complex addition are required, so that the calculation amount is huge and the realization is difficult. A general solution is to use an optimized FFT kernel for near real-time processing. The FFT implementation method divides the spectrum range of the whole [0:
Figure BDA0003045557560000062
the spectrum analysis realized by the above formula is the DFT method. The calculation amount of DFT analysis is L x N complex multiplication and L (N-1) complex addition, and the smaller the target analysis frequency interval [ L: L-1] is, the smaller the calculation amount is.
The DFT method has variable parameters: spectral analysis interval, spectral analysis precision. When the algorithm initially works, in a larger analysis interval, a signal which possibly exists is rapidly detected by a larger frequency stepping value (poorer analysis precision); when the existence of the signal is detected, the suspected signal is set as the center, and fine search is carried out in a smaller frequency step within a smaller range around the frequency center to obtain accurate signal frequency.
The present embodiment is explained with reference to a specific example in which the system frequency fs =62MHz, and the signal frequency of interest f2 is at 2 × (15.58- Δf ) At MHz, the f2 bandwidth is 40KHz, as known from the carrier tracking bandwidth (+ -10 KHz).
According to the sampling theorem, the minimum sampling frequency of f2 can be effectively recovered to be 80KHz (2 x f 2), namely the resampling frequency fs/D after extraction is more than 80KHz when the Zoom-DFT calculation is carried out, and 775 times of extraction (62 MHz/80 KHz) can be carried out at the maximum, which is an ideal condition for completely converting the signal to the baseband. When using the fs/2 frequency conversion method without multiplication, f2 is multiplied by a cosine sequence with frequency (fs/2 = 31) MHz, then f2 after mixing becomes:
F2=2*(15.58- Δf )-fs/2=31.16-2* Δf -31=0.16-2* Δf (MHz),
the above formula shows that the center frequency of the F4 after mixing is 0.16MHz and the bandwidth is 40KHz.
In the specific implementation of the carrier frequency offset estimation, by setting the 3dB bandwidth of the narrow-band filter to 40KHz, the decimation multiple D =40, and the dft point number N =2048, a spectrogram with a resolution of 756.8Hz (62M/(40 × 2048)) can be obtained. The peak detection output is 0.16 MHz-2- Δf The spectrogram contains 2 × Δf Component, so that the output frequency offset can be actually detected Δf The resolution of (1) is 756.8Hz/2=378.4Hz, i.e. the frequency offset estimation algorithm can accurately distinguish the carrier frequency offset of 378.4 Hz.
To more clearly show the advantages of this embodiment, the simulation is implemented in the FPGA. Please refer to fig. 2-5.
In fig. 2 and 3, no noise is actively added to the input signal, the signal-to-noise ratio is high, the DFT output spectrum is clean, the peak is obvious, and the amplitude is high. The abscissa is the DFT output sequence number, one step represents the frequency offset of the signal carrier 4 af, and the ordinate is the frequency amplitude. The frequency offset of 7.5KHz brings 20 serial number offsets, which is accurate consistent with a resolution of about 378.4Hz for af in the above example analysis.
In fig. 4 and 5, after the signal is mixed with noise with a certain signal-to-noise ratio, the DFT output peak is obviously decreased, and the bottom noise is obviously increased.
When the signal-to-noise ratio of the signal is = -1dB, the module can still accurately estimate the carrier frequency offset, and the stability is good.
The method of the embodiment directly carries out frequency domain processing on the signal correlation value, eliminates the modulation phase ambiguity, does not need to use a frequency discriminator, only needs to use a DFT algorithm to carry out spectrum analysis, and greatly expands the frequency discrimination range. In addition, the DFT algorithm only needs to use little resource consumption, and the frequency discrimination resolution can be greatly improved.
Referring to fig. 6, fig. 6 is a block diagram of a spread spectrum signal frequency detection module provided in the present invention, where the module may be independently disposed, or may be integrated into a system as a functional module of the system, and includes:
the system comprises a signal receiving unit, a signal processing unit and a signal processing unit, wherein the signal receiving unit is used for receiving a spread spectrum signal to be detected, and the spread spectrum signal to be detected comprises modulation information;
the signal modulation unit is used for carrying out square operation on the spread spectrum signal to be detected so as to eliminate the modulation signal and obtain a frequency estimation signal;
the filtering unit is used for sequentially carrying out down-conversion and frequency mixing on the frequency estimation signal and carrying out narrow-band filtering on the frequency-mixed signal;
and the detection unit is used for performing DFT operation on the filtered signal function so as to detect the carrier frequency deviation.
In a specific embodiment, the expression of the spread spectrum signal to be detected is:
r(t)=A*P(t)*e j(2π(f+Δf)t+θ) +n(t),
wherein A is signal power, P (t) is shaping pulse shape, f is carrier frequency, Δ f is carrier frequency offset, and θ is modulation information phase of current symbol,
Figure BDA0003045557560000091
n (t) is noise.
In one embodiment, the frequency estimation signal expression is:
r 2 '(t)=A 2 *e j(2π(2*f+2*Δf)t) +n 2 (t),
wherein, A is signal power, f is carrier frequency, delta f is carrier frequency deviation, and n (t) is noise.
In one embodiment, the frequency of the mixed signal is: f2=2 × F2-fs/2, where F2 is the signal frequency of interest before mixing and fs is the system frequency.
As will be appreciated by one of skill in the art, embodiments of the present application may be provided as a method, module (apparatus, device), or computer program product. Accordingly, this application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects that may all generally be referred to herein as a "module" or "system. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein. A computer program stored/distributed on a suitable medium supplied together with or as part of other hardware, may also take other distributed forms, such as via the Internet or other wired or wireless telecommunication systems.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (devices) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing is a further detailed description of the invention in connection with specific preferred embodiments and it is not intended to limit the invention to the specific embodiments described. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (2)

1. A method for frequency detection of a spread spectrum signal, comprising:
receiving a spread spectrum signal to be detected, wherein the spread spectrum signal to be detected comprises modulation information;
carrying out square operation on the spread spectrum signal to be detected to eliminate the modulation information and obtain a frequency estimation signal;
sequentially carrying out down-conversion and frequency mixing on the frequency estimation signal, and carrying out narrow-band filtering on the frequency-mixed signal;
performing Discrete Fourier Transform (DFT) operation on the filtered signal function x (n) to detect carrier frequency deviation;
the expression of the to-be-detected spread spectrum signal is as follows:
r(t)=A*P(t)*e j(2π(f+Δf)t+θ) +n(t),
wherein A is signal power, P (t) is shaping pulse shape, f is carrier frequency, delta f is carrier frequency deviation, theta is modulation information phase of current symbol,
Figure FDA0003932218990000011
n (t) is noise;
the frequency estimation signal expression is:
r 2 '(t)=A 2 *e j(2π(2*f+2*Δf)t) +n 2 (t),
wherein, A is signal power, f is carrier frequency, delta f is carrier frequency offset, and n (t) is noise;
the frequency of the mixed signal is: f2=2 × F2-fs/2, where F2 is the signal frequency of interest before mixing and fs is the system frequency;
the frequency spectrum range of the whole [0:
Figure FDA0003932218990000012
in the formula, W N =e -j*2π/N
2. A spread spectrum signal frequency detection module, comprising:
the system comprises a signal receiving unit, a signal processing unit and a signal processing unit, wherein the signal receiving unit is used for receiving a spread spectrum signal to be detected, and the spread spectrum signal to be detected comprises modulation information;
the signal modulation unit is used for carrying out square operation on the spread spectrum signal to be detected so as to eliminate the modulation information and obtain a frequency estimation signal;
the filtering unit is used for sequentially carrying out down-conversion and frequency mixing on the frequency estimation signal and carrying out narrow-band filtering on the frequency-mixed signal;
a detection unit for performing Discrete Fourier Transform (DFT) operation on the filtered signal function x (n) to detect a carrier frequency deviation;
the expression of the spread spectrum signal to be detected is as follows:
r(t)=A*P(t)*e j(2π(f+Δf)t+θ) +n(t),
wherein A is signal power, P (t) is shaping pulse shape, f is carrier frequency, Δ f is carrier frequency offset, and θ is modulation information phase of current symbol,
Figure FDA0003932218990000021
n (t) is noise;
the frequency estimation signal expression is:
r 2 '(t)=A 2 *e j(2π(2*f+2*Δf)t) +n 2 (t),
wherein, A is signal power, f is carrier frequency, delta f is carrier frequency deviation, and n (t) is noise;
the frequency of the mixed signal is: f2=2 × F2-fs/2, where F2 is the frequency of the signal of interest before mixing and fs is the system frequency;
the frequency spectrum range of the whole [0:
Figure FDA0003932218990000022
in the formula, W N =e -j*2π/N
CN202110471627.5A 2021-04-29 2021-04-29 Method and module for detecting frequency of spread spectrum signal Active CN113162652B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110471627.5A CN113162652B (en) 2021-04-29 2021-04-29 Method and module for detecting frequency of spread spectrum signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110471627.5A CN113162652B (en) 2021-04-29 2021-04-29 Method and module for detecting frequency of spread spectrum signal

Publications (2)

Publication Number Publication Date
CN113162652A CN113162652A (en) 2021-07-23
CN113162652B true CN113162652B (en) 2023-01-06

Family

ID=76872196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110471627.5A Active CN113162652B (en) 2021-04-29 2021-04-29 Method and module for detecting frequency of spread spectrum signal

Country Status (1)

Country Link
CN (1) CN113162652B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5982809A (en) * 1997-01-23 1999-11-09 Cwill Telecommunications, Inc. Method and apparatus for carrier recovery and compensation in spread spectrum communications
EP1662736A2 (en) * 2004-11-25 2006-05-31 NTT DoCoMo, Inc. Time synchronization and frequency offset estimation in OFDM systems
CN104901718A (en) * 2015-05-22 2015-09-09 哈尔滨工程大学 Doppler estimation method based on measurement of carrier frequency of direct sequence spread spectrum signal

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3296421B2 (en) * 1998-09-28 2002-07-02 日本電気株式会社 Frequency estimation device
JP3633497B2 (en) * 2001-03-22 2005-03-30 三菱電機株式会社 Receiver for performing frequency error estimation and method for estimating frequency error
US7639731B1 (en) * 2006-05-03 2009-12-29 Rockwell Collins, Inc. Carrier frequency estimation via symbol rate estimation
CN202374360U (en) * 2011-12-16 2012-08-08 陕西凌云电器集团有限公司 Multi-mode digital-analog integrated tuner
CN103929391B (en) * 2013-01-15 2017-10-20 电信科学技术研究院 A kind of transmitting frequency calibration method and device
CN110764067A (en) * 2019-10-15 2020-02-07 哈尔滨工程大学 Fourier transform LFM-BPSK composite modulation radar signal parameter estimation method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5982809A (en) * 1997-01-23 1999-11-09 Cwill Telecommunications, Inc. Method and apparatus for carrier recovery and compensation in spread spectrum communications
EP1662736A2 (en) * 2004-11-25 2006-05-31 NTT DoCoMo, Inc. Time synchronization and frequency offset estimation in OFDM systems
CN104901718A (en) * 2015-05-22 2015-09-09 哈尔滨工程大学 Doppler estimation method based on measurement of carrier frequency of direct sequence spread spectrum signal

Also Published As

Publication number Publication date
CN113162652A (en) 2021-07-23

Similar Documents

Publication Publication Date Title
CN101776752B (en) Precise tracking and measuring method of high dynamic signal of air fleet link
JP7228507B2 (en) Method and apparatus for reception of DS/FH signals
JP5310333B2 (en) Reception device, signal processing method, and program
CN106603451B (en) High dynamic Doppler frequency offset and frequency offset change rate estimation method based on time delay autocorrelation
US8462616B2 (en) Apparatus and method for estimating a frequency shift and a time shift
CN110071738A (en) Spread-spectrum signal based on more pseudo-code branch receivers de-spreads tracking
US8611472B2 (en) Carrier frequency offset estimation scheme, for digital standards with MPSK modulated preamble
US7620125B1 (en) Frequency estimation apparatus and related method
CN110300079B (en) MSK signal coherent demodulation method and system
JP2005031073A (en) Search method of gps correlated peak signal and system therefor
CN109067680B (en) Carrier frequency offset estimation method and device of baseband signal
JP2014228536A (en) Apparatus and methods for determining status of tracking loop
CN105610745B (en) It is a kind of for the Fast carrier offset estimation of fsk signal and bearing calibration
CN105516041A (en) Adaptive digital demodulation system at low signal to noise ratio
CN114531329A (en) Method, system and application for estimating carrier frequency of multi-path MSK (minimum shift keying) signals
CN113162652B (en) Method and module for detecting frequency of spread spectrum signal
KR101015890B1 (en) Signal acquisition method and apparatus of GNSS receiver
CN107733824A (en) A kind of carrier synchronization method based on AFC loops
JP4700518B2 (en) Synchronization timing detection device and receiver
Rice et al. Maximum likelihood carrier phase synchronization in FPGA-based software defined radios
CN114301552B (en) Digital modulation signal testing method and system
JP2002208974A (en) Differential phase demodulator incorporating quaternary coherent phase tracking
CN110958035A (en) Frequency secondary capturing device and method based on triangular envelope fitting
JPH05256883A (en) Digital method and device for analyzing frequency
JP2006217601A (en) Method for acquiring positioning signal of geographic localization system, receiver for geographic localization system and computer data carrier comprising program instruction for carrying out the method

Legal Events

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