CN115118564A - Carrier frequency deviation estimation method and device - Google Patents

Carrier frequency deviation estimation method and device Download PDF

Info

Publication number
CN115118564A
CN115118564A CN202210699989.4A CN202210699989A CN115118564A CN 115118564 A CN115118564 A CN 115118564A CN 202210699989 A CN202210699989 A CN 202210699989A CN 115118564 A CN115118564 A CN 115118564A
Authority
CN
China
Prior art keywords
sequence
carrier frequency
phase
phase difference
frequency deviation
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
CN202210699989.4A
Other languages
Chinese (zh)
Other versions
CN115118564B (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.)
Hunan Econavi Technology Co Ltd
Original Assignee
Hunan Econavi 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 Hunan Econavi Technology Co Ltd filed Critical Hunan Econavi Technology Co Ltd
Priority to CN202210699989.4A priority Critical patent/CN115118564B/en
Publication of CN115118564A publication Critical patent/CN115118564A/en
Application granted granted Critical
Publication of CN115118564B publication Critical patent/CN115118564B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2689Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2695Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
    • 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
    • H04L2027/0026Correction of carrier offset
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

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

Abstract

The invention discloses a carrier frequency deviation estimation method and a device, and the method comprises the following steps: s01, inputting a sampling sequence received by a digital intermediate frequency receiver, obtaining a complex baseband sequence after down-conversion, and calculating a phase angle of each sampling point in the complex baseband sequence; s02, calculating phase differences between adjacent sampling points according to the calculated phase angles; s03, correcting according to the phase difference between the adjacent sampling points obtained through calculation to remove phase folding, and obtaining a corrected phase difference sequence; and S04, calculating carrier frequency deviation output according to the corrected phase difference sequence. The invention has the advantages of simple realization method, low complexity, small calculation amount, compatibility with various linear modulation patterns, wide application range and the like.

Description

Carrier frequency deviation estimation method and device
Technical Field
The present invention relates to the field of non-cooperative communication technologies, and in particular, to a carrier frequency offset estimation method and apparatus.
Background
In cooperative communication, carrier frequency deviation is mainly caused by doppler shift and unequal transceiver oscillator frequencies. In non-cooperative communication, because it is difficult for a receiver to know an accurate value of carrier frequency of a transmitter, carrier frequency deviation after down-conversion is often very large and even can be compared with signal bandwidth. The carrier frequency deviation can rotate the signal constellation point, so that carrier synchronization is required for correct demodulation to enable the carrier at the receiving end to have the same frequency as that at the sending end, namely, the corresponding frequency deviation needs to be estimated and compensated, otherwise, the problems of performance reduction of a demodulator, large demodulation loss and the like can be caused, and even correct information data can not be demodulated.
If there is no modulation information, the down-converted received signal is a single frequency signal with a frequency equal to the carrier frequency offset. If the modulation information can be removed from the received signal, a frequency estimation method of a single frequency signal in noise can be used to estimate the carrier frequency offset. In fact, in the data-aided method of the training sequence in the prior art, the main role of the training sequence is to remove the influence of the modulation information on the received signal. For non-cooperative communication, the data-aided approach is not applicable since no a priori information is available, and the pilot information or synchronization header is often unknown. Therefore, in the demodulation process of the non-cooperative signal, only a non-data aided method, generally a method of non-linear transformation, such as a square law method, a cyclic accumulation method, etc., can be used to remove the modulation information in order to obtain the carrier frequency deviation of the signal.
In the prior art, the square law method is to perform M power on a signal correspondingly according to the modulation order M of the signal, discrete spectral lines appear in a frequency spectrum at M times of a carrier frequency, and carrier frequency estimation can be completed by searching frequency points corresponding to the spectral lines, but the calculation amount is large, and the problem of noise power amplification after square or high power exists, so that the square law method is not suitable for high-order modulation signals. The cyclic accumulation method is to extract carrier frequency information by using cyclic moment and cyclic accumulation, has high estimation precision, but has the problems of large calculation amount and high algorithm complexity for high-order modulation signals. For carrier frequency offset estimation, the prior art usually adopts differential operation between symbols, but this kind of method only uses the optimal sampling point of each data symbol for estimation, and must require the receiver to complete symbol synchronization before estimation, and the amount of calculation is large and the complexity is high.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: aiming at the technical problems in the prior art, the invention provides the carrier frequency deviation estimation method and the carrier frequency deviation estimation device which have the advantages of simple implementation method, low complexity, small calculation amount, compatibility with various linear modulation patterns and wide application range.
In order to solve the technical problems, the technical scheme provided by the invention is as follows:
a carrier frequency offset estimation method, comprising the steps of:
s01, inputting a sampling sequence received by a digital intermediate frequency receiver, obtaining a complex baseband sequence after down-conversion, and calculating a phase angle of each sampling point in the complex baseband sequence;
s02, calculating phase differences between adjacent sampling points according to the calculated phase angles;
s03, correcting according to the phase difference between the adjacent sampling points obtained through calculation to remove phase folding, and obtaining a corrected phase difference sequence;
and S04, calculating carrier frequency deviation output according to the corrected phase difference sequence.
Further, in step S01, a cordic algorithm is used to calculate the instantaneous phase angle of the complex baseband sequence, and the calculation expression is:
Figure BDA0003704021100000021
wherein,
Figure BDA0003704021100000022
for instantaneous phase angle sequences, n denotes the time of day, phi a (n) phase at time n of data symbol, Im [ x (n)]Is the imaginary part of the complex base band sequence x (n), Re [ x (n)]Is the real part, T, of the complex baseband sequence x (n) s For the sampling interval time, Δ f is the carrier frequency deviation, Φ N And (n) is phase disturbance caused by channel noise.
Further, in step S02, the instantaneous phase difference between adjacent sampling points is specifically calculated according to the following formula:
Figure BDA0003704021100000023
wherein,
Figure BDA0003704021100000024
is the instantaneous phase angle at time n,
Figure BDA0003704021100000025
is the phase angle at the time n-1,
Figure BDA0003704021100000026
for a sequence of instantaneous phase differences, T s For the sampling interval time, Δ f is the carrier frequency deviation, Φ a (n) is the phase at time n of the data symbol, phi N And (n) is phase disturbance caused by channel noise at the time n.
Further, in step S03, the phase difference between the adjacent sampling points is corrected specifically according to the following formula:
Figure BDA0003704021100000027
wherein,
Figure BDA0003704021100000028
is a sequence of phase angles and is,
Figure BDA0003704021100000029
in order to be a sequence of instantaneous phase differences,
Figure BDA00037040211000000210
for the corrected instantaneous phase difference sequence, n represents the time of day.
Further, after the step S03 and before the step S04, low-pass filtering is performed on the corrected phase difference sequence to obtain a filtered signal output.
Further, in step S04, the carrier frequency offset is specifically calculated according to the following formula:
Figure BDA00037040211000000211
wherein,
Figure BDA00037040211000000212
low-pass filtered signal, T, for said corrected phase difference sequence s Is the sampling interval time.
A carrier frequency offset estimation apparatus comprising:
the phase angle calculation module is used for inputting a sampling sequence received by the digital intermediate frequency receiver, obtaining a complex baseband sequence after down-conversion, and calculating the phase angle of each sampling point in the complex baseband sequence;
the phase difference calculation module is used for calculating the phase difference between adjacent sampling points according to each phase angle obtained by calculation;
the phase folding removal module is used for correcting according to the phase difference between the adjacent sampling points obtained through calculation so as to remove phase folding and obtain a corrected phase difference sequence;
and the frequency deviation calculation module is used for calculating carrier frequency deviation output according to the corrected phase difference sequence.
Further, a low-pass filtering module is connected between the phase folding removal module and the frequency deviation calculation module, and is used for performing low-pass filtering on the corrected phase difference sequence to obtain a filtered signal output.
A computer device comprising a processor and a memory for storing a computer program, the processor being adapted to execute the computer program to perform the method as described above.
A computer-readable storage medium having stored thereon a computer program which, when executed, implements the method as described above.
Compared with the prior art, the invention has the advantages that:
1. the carrier frequency deviation estimation method and the carrier frequency deviation estimation device do not need to remove modulation information in signals, carry out carrier frequency deviation estimation on baseband signals directly, have low complexity and small calculated amount, have an estimation process irrelevant to modulation orders, can be compatible with various linear modulation modes such as FSK/MSK/MPSK/DMPSK/MAPSK/MQAM and the like, simultaneously estimate by utilizing all sampling data sequences without paying attention to sampling position deviation, realize carrier frequency deviation estimation before symbol synchronization and have a very large adaptive carrier frequency deviation range.
2. The carrier frequency deviation estimation method and the carrier frequency deviation estimation device can be suitable for non-cooperative communication and high-order modulation signals, the non-cooperative communication does not need to be identified in a modulation mode in advance, the high-order modulation signals do not need to be subjected to high-order nonlinear calculation, the problem of noise power amplification does not exist, the complexity is not increased, and the realization complexity and the calculation amount can be greatly reduced.
3. The carrier frequency deviation estimation method and the device realize carrier frequency deviation estimation by combining differential operation of adjacent sampling points, can remove modulation information by the differential operation, do not need to remove the modulation information by utilizing prior information such as training sequences and the like or by nonlinear transformation like the traditional estimation method, further combine a narrow-band low-pass filter, can filter out phase difference distortion values at jump positions of adjacent symbols, do not need prior information such as the training sequences and the like or remove the modulation information by nonlinear transformation, can quickly and efficiently realize accurate carrier frequency deviation estimation, have low complexity and are easy to realize engineering.
Drawings
Fig. 1 is a schematic flow chart of a carrier frequency offset estimation method according to this embodiment.
Fig. 2 is a schematic structural diagram of the carrier frequency offset estimation device of the present embodiment.
Fig. 3 is a schematic flow chart of the carrier frequency offset estimation implemented in the embodiment of the present invention.
Detailed Description
The invention is further described below with reference to the drawings and specific preferred embodiments of the description, without thereby limiting the scope of protection of the invention.
As shown in fig. 1 and 2, the carrier frequency offset estimation method of the present embodiment includes the steps of:
s01, inputting a sampling sequence received by a digital intermediate frequency receiver, obtaining a complex baseband sequence after down-conversion, and calculating a phase angle of each sampling point in the complex baseband sequence;
s02, calculating phase differences between adjacent sampling points according to the calculated phase angles;
s03, correcting according to the phase difference between the adjacent sampling points obtained through calculation to remove phase folding, and obtaining a corrected phase difference sequence;
and S04, calculating carrier frequency deviation output according to the corrected phase difference sequence.
Considering that the sampling frequency in the digital intermediate frequency receiver is far greater than the symbol rate, the embodiment implements carrier frequency deviation estimation by combining the differential operation of adjacent sampling points on the basis that the used sampling frequency is far greater than the symbol rate, and since most of the adjacent sampling points are in the same data symbol, the modulation information of the adjacent sampling points can be considered to be the same, the modulation information can be removed by the differential operation without using prior information such as training sequences or removing the modulation information by nonlinear transformation as in the conventional method. And because all the sampled data sequences are used for estimation, rather than only the optimal sampling point of each data symbol as in the traditional estimation mode, the estimation does not need to pay attention to the sampling position deviation, and the estimation of the carrier frequency deviation is not influenced even if the symbol synchronization is not finished by signals, the estimation of the carrier frequency deviation can be realized before the symbol synchronization, the carrier synchronization is finished, so that the selection of a receiver symbol synchronization algorithm is wider, and the realization complexity is greatly reduced.
For the digital intermediate frequency receiver, the sampling sequence r (n) of the received signal is subjected to digital down-conversion to obtain a complex baseband signal sequence x (n), which can be expressed as:
Figure BDA0003704021100000041
where A is the signal amplitude, a (mTs) is the transmitted data symbol, g [ (n-m) T s ]For shaping the filtered impulse response, M is the shaping filter order, T s For the sampling interval, Δ f is the carrier frequency deviation, θ is the initial phase of the carrier, N (nT) s ) The signal is sampled for noise.
In step S01 of this embodiment, a cordic algorithm is specifically adopted to calculate the phase of each sample point in the complex baseband signal sequence x (n) according to the following formula, so as to obtain an instantaneous phase sequence
Figure BDA0003704021100000043
Figure BDA0003704021100000042
Wherein, Im [ x (n)]Is the imaginary part of the complex base band sequence x (n), Re [ x (n)]Is the real part, Φ, of the complex baseband sequence x (n) a (n) is the phase at time n of the data symbol, related only to the transmitted data symbol and the shaping filter, phi N And (n) is phase disturbance caused by channel noise.
The amplitude and phase angle of the complex baseband signal sequence x (n) satisfy the relation:
Figure BDA0003704021100000051
wherein a (n) is the amplitude value of the complex baseband signal sequence x (n),
Figure BDA0003704021100000052
is the phase angle of the complex baseband signal sequence x (n).
Calculating instantaneous phase sequence as described above
Figure BDA0003704021100000053
Thereafter, step S02 of the present embodiment further calculates the instantaneous phase sequence
Figure BDA0003704021100000054
In every two adjacent miningPhase difference between the samples. The frequency is the first difference of the instantaneous phase sequence, i.e.:
Figure BDA0003704021100000055
wherein,
Figure BDA0003704021100000056
is the instantaneous phase angle at time n,
Figure BDA0003704021100000057
is the phase angle at the time instant n-1,
Figure BDA0003704021100000058
is the nth instantaneous phase difference.
Calculating the phase difference between two adjacent sampling points of the complex baseband signal sequence x (n)
Figure BDA0003704021100000059
The specific expression of (A) is as follows:
Figure BDA00037040211000000510
since the phase of x (n) is calculated modulo 2 pi, there is a phase fold that needs to be corrected for the instantaneous phase difference of adjacent samples. In step S03, the present embodiment corrects the phase difference between adjacent sampling points according to the following formula:
Figure BDA00037040211000000511
wherein,
Figure BDA00037040211000000512
is a corrected instantaneous phase difference sequence.
In this embodiment, after step S03 and before step S04, the method further includes performing low-pass filtering on the corrected phase difference sequence to obtain a filtered signal output.
Sampling frequency f in digital intermediate frequency receiver s Will be much larger than the data symbol rate R s I.e. the sampling interval time Ts is much smaller than the width of the transmitted data symbol, the phases of two adjacent sampling points of the transmitted signal can be considered equal except near the transition point of the transmitted data symbol, i.e.:
Φ a (n)≈Φ a (n-1) (7)
and [ phi ] because the transmitted data symbols are random a (n)-Φ a (n-1)]Has a mean value of zero, [ phi ] N (n) is the phase perturbation due to channel noise, thus [ phi ] N (n)-Φ N (n-1)]Should also be zero. Calculated by equation (5)
Figure BDA00037040211000000513
Passing through a narrow band (the passband bandwidth is much less than the symbol rate R s ) After the low-pass filtering, the phase distortion value near the data symbol jumping point and the phase disturbance caused by channel noise can be sent, namely the last two terms ([ phi ] of the right side of the equation (5) are filtered a (n)-Φ a (n-1 and. PHINn-1). Then to the corrected instantaneous phase difference sequence
Figure BDA00037040211000000514
After narrow-band low-pass filtering, the following components are obtained:
Figure BDA00037040211000000515
for the corrected instantaneous phase difference sequence
Figure BDA00037040211000000516
After narrow-band low-pass filtering, phase disturbance caused by data symbols and noise can be filtered out to obtain filtered signals
Figure BDA00037040211000000517
The passband of the low-pass filter is much smaller than the symbol rate R of the signal s Start of stop bandFrequency is also less than R s
In step S04 of this embodiment, the carrier frequency deviation is calculated according to the following formula specifically based on the relationship between the frequency and the phase:
Figure BDA0003704021100000061
wherein,
Figure BDA0003704021100000062
low-pass filtered signal, T, for the corrected phase difference sequence s Is the sampling period.
In the embodiment, by adopting a differential operation method of adjacent sampling points, modulation information in signals does not need to be removed, carrier frequency offset estimation is directly performed on baseband signals, and a narrow-band low-pass filter is further combined to filter out a phase difference distortion value at a jump position of adjacent symbols, so that accurate carrier frequency offset estimation can be quickly and efficiently realized.
The adaptive carrier frequency deviation range is an important index for estimating the carrier frequency deviation, and the sampling frequency f of the invention is s Much larger than the symbol rate R s And the difference operation between adjacent sampling points is adopted, and the carrier frequency offset range is ((-f) s /4,f s 4) greater than the symbol rate R of the signal s In a non-cooperative communication scene, the method can be used for directly carrying out carrier frequency offset estimation on the signal without pretreatment such as carrier frequency rough estimation and the like.
In a specific application embodiment, as shown in fig. 3, a detailed procedure for implementing carrier frequency offset estimation includes:
step 1: the amplitude and phase angle of the complex baseband sequence x (n) are calculated using cordic algorithm.
Figure BDA0003704021100000063
Figure BDA0003704021100000064
Figure BDA0003704021100000065
Wherein a (n) is the amplitude value of x (n),
Figure BDA00037040211000000612
is the phase angle of x (n).
Step 2: calculating the instantaneous phase difference of adjacent sampling points:
Figure BDA0003704021100000066
and step 3: removing phase folding, and correcting the instantaneous phase difference of adjacent sampling points:
Figure BDA0003704021100000067
and 4, step 4: for corrected instantaneous phase difference sequence
Figure BDA0003704021100000068
Narrow-band low-pass filtering is carried out to filter out phase disturbance caused by data symbols and noise to obtain filtered signals
Figure BDA0003704021100000069
The passband of the low-pass filter is much smaller than the symbol rate R of the signal s The starting frequency of the stop band is less than R s
And 5: according to frequency and phase relationship, according to filtered
Figure BDA00037040211000000610
Calculating a frequency value to obtain a carrier frequency deviation:
Figure BDA00037040211000000611
as shown in fig. 2, the carrier frequency offset estimation device of the present embodiment includes:
the phase angle calculation module is used for inputting a sampling sequence received by the digital intermediate frequency receiver, obtaining a complex baseband sequence after down-conversion, and calculating the phase angle of each sampling point in the complex baseband sequence;
the phase difference calculation module is used for calculating the phase difference between adjacent sampling points according to each phase angle obtained by calculation;
the phase folding removal module is used for correcting according to the phase difference between the adjacent sampling points obtained through calculation so as to remove phase folding and obtain a corrected phase difference sequence;
and the frequency deviation calculation module is used for calculating the carrier frequency deviation output according to the corrected phase difference sequence.
In this embodiment, a low-pass filtering module is further connected between the phase folding removal module and the frequency deviation calculation module, so as to perform low-pass filtering on the corrected phase difference sequence, and obtain a filtered signal output.
As shown in fig. 2, the phase angle calculation module includes a down-conversion circuit composed of two multipliers and a phase angle calculation circuit, the down-conversion circuit inputs a sampling sequence r (n), outputs two paths of intersecting signals i (n) and q (n) through down-conversion, and calculates the phase angle of a complex baseband sequence x (n) by the phase angle calculation circuit according to a cordic algorithm; the phase difference calculation module samples a differential circuit and calculates an instantaneous phase difference sequence according to a formula (4)
Figure BDA0003704021100000071
Outputting the filtered signal after correction and low-pass filtering
Figure BDA0003704021100000072
And calculating the final frequency deviation delta f and outputting the final frequency deviation delta f through a frequency deviation calculation module according to the formula (8).
The carrier frequency offset estimation apparatus of the present embodiment corresponds to the carrier frequency offset estimation method one by one, and is not described herein again.
The present embodiment also provides a computer device comprising a processor and a memory, the memory being configured to store a computer program, the processor being configured to execute the computer program to perform the method as described above.
The present embodiment is a computer-readable storage medium storing a computer program, which when executed implements the method described above.
The invention can quickly and efficiently realize the carrier frequency offset estimation of the baseband signal, has very wide modulation modes, can be compatible with various linear modulation modes such as FSK/MSK/MPSK/DMPSK/MAPSK/MQAM and the like, does not need to change the realization structure and parameter configuration, is independent of the modulation order, does not need to carry out modulation mode identification in advance if being applied to non-cooperative communication, does not need to carry out high-order nonlinear calculation when being applied to high-order modulation signals, does not have the problem of noise power amplification, does not increase the complexity, and can greatly reduce the realization complexity and the calculated amount.
The foregoing is considered as illustrative of the preferred embodiments of the invention and is not to be construed as limiting the invention in any way. Although the present invention has been described with reference to the preferred embodiments, it is not intended to be limited thereto. Therefore, any simple modification, equivalent change and modification made to the above embodiments according to the technical essence of the present invention shall fall within the protection scope of the technical solution of the present invention, unless the technical essence of the present invention departs from the content of the technical solution of the present invention.

Claims (10)

1.一种载波频率偏差估计方法,其特征在于,步骤包括:1. a carrier frequency deviation estimation method, is characterized in that, step comprises: S01.输入数字中频接收机接收到的采样序列,经过下变频后得到复基带序列,计算所述复基带序列中每个采样点的相位角;S01. Input the sampling sequence received by the digital intermediate frequency receiver, obtain a complex baseband sequence after down-conversion, and calculate the phase angle of each sampling point in the complex baseband sequence; S02.根据计算得到个各所述相位角计算相邻采样点之间的相位差;S02. Calculate the phase difference between the adjacent sampling points according to the calculated phase angles; S03.根据计算得到的所述相邻采样点之间的相位差进行校正,以去除相位折叠,得到校正后的相位差序列;S03. Correct according to the calculated phase difference between the adjacent sampling points, to remove the phase folding, and obtain the corrected phase difference sequence; S04.根据所述校正后的相位差序列计算载波频率偏差输出。S04. Calculate the carrier frequency deviation output according to the corrected phase difference sequence. 2.根据权利要求1所述的载波频率偏差估计方法,其特征在于,所述步骤S01中采用cordic算法计算复基带序列的瞬时相位角,计算表达式为:2. carrier frequency deviation estimation method according to claim 1, is characterized in that, adopts cordic algorithm to calculate the instantaneous phase angle of complex baseband sequence in described step S01, and calculation expression is:
Figure FDA0003704021090000011
Figure FDA0003704021090000011
其中,
Figure FDA0003704021090000012
为瞬时相位角序列,n表示时刻,Φa(n)为数据符号的n时刻的相位,Im[x(n)]为复基带序列x(n)的虚部,Re[x(n)]为复基带序列x(n)的实部,Ts为采样间隔时间,Δf为载波频率偏差,ΦN(n)为信道噪声引起的相位扰动。
in,
Figure FDA0003704021090000012
is the instantaneous phase angle sequence, n represents time, Φ a (n) is the phase of the data symbol at time n, Im[x(n)] is the imaginary part of the complex baseband sequence x(n), Re[x(n)] is the real part of the complex baseband sequence x(n), T s is the sampling interval time, Δf is the carrier frequency deviation, and Φ N (n) is the phase disturbance caused by channel noise.
3.根据权利要求1所述的载波频率偏差估计方法,其特征在于,所述步骤S02中,具体按照下式计算相邻采样点的瞬时相位差:3. carrier frequency deviation estimation method according to claim 1, is characterized in that, in described step S02, specifically calculates the instantaneous phase difference of adjacent sampling point according to following formula:
Figure FDA0003704021090000013
Figure FDA0003704021090000013
其中,
Figure FDA0003704021090000014
为n时刻的瞬时相位角,
Figure FDA0003704021090000015
为n-1时刻的相位角,
Figure FDA0003704021090000016
为瞬时相位差序列,Ts为采样间隔时间,Δf为载波频率偏差,Φa(n)为数据符号的n时刻的相位,ΦN(n)为n时刻信道噪声引起的相位扰动。
in,
Figure FDA0003704021090000014
is the instantaneous phase angle at time n,
Figure FDA0003704021090000015
is the phase angle at time n-1,
Figure FDA0003704021090000016
is the instantaneous phase difference sequence, T s is the sampling interval, Δf is the carrier frequency deviation, Φ a (n) is the phase of the data symbol at time n, and Φ N (n) is the phase disturbance caused by channel noise at time n.
4.根据权利要求1所述的载波频率偏差估计方法,其特征在于,所述步骤S03中,具体按照下式对所述相邻采样点之间的相位差进行校正:4. The carrier frequency deviation estimation method according to claim 1, wherein in the step S03, the phase difference between the adjacent sampling points is corrected according to the following formula:
Figure FDA0003704021090000017
Figure FDA0003704021090000017
其中,
Figure FDA0003704021090000018
为瞬时相位角序列,
Figure FDA0003704021090000019
为瞬时相位差序列,
Figure FDA00037040210900000110
为校正后的瞬时相位差序列,n表示时刻。
in,
Figure FDA0003704021090000018
is the instantaneous phase angle sequence,
Figure FDA0003704021090000019
is the instantaneous phase difference sequence,
Figure FDA00037040210900000110
is the corrected instantaneous phase difference sequence, and n represents the time.
5.根据权利要求1~4中任意一项所述的载波频率偏差估计方法,其特征在于,所述步骤S03后、步骤S04前还包括对所述校正后的相位差序列进行窄带低通滤波,得到滤波后信号输出。5. The carrier frequency deviation estimation method according to any one of claims 1 to 4, characterized in that, after the step S03 and before the step S04, it further comprises performing narrow-band low-pass filtering on the corrected phase difference sequence , to get the filtered signal output. 6.根据权利要求5所述的载波频率偏差估计方法,其特征在于,所述步骤S04中,具体按照下式计算载波频率偏差:6. The carrier frequency deviation estimation method according to claim 5, wherein in the step S04, the carrier frequency deviation is calculated according to the following formula:
Figure FDA0003704021090000021
Figure FDA0003704021090000021
其中,
Figure FDA0003704021090000022
为所述校正后的相位差序列进行低通滤波后信号,Ts为采样间隔时间。
in,
Figure FDA0003704021090000022
The signal after low-pass filtering is performed for the corrected phase difference sequence, and T s is the sampling interval time.
7.一种载波频率偏差估计装置,其特征在于,包括:7. A carrier frequency deviation estimation device, characterized in that, comprising: 相位角计算模块,用于输入数字中频接收机接收到的采样序列,经过下变频后得到复基带序列,计算所述复基带序列中每个采样点的相位角;The phase angle calculation module is used to input the sampling sequence received by the digital intermediate frequency receiver, obtain a complex baseband sequence after down-conversion, and calculate the phase angle of each sampling point in the complex baseband sequence; 相位差计算模块,用于根据计算得到个各所述相位角计算相邻采样点之间的相位差;a phase difference calculation module, configured to calculate the phase difference between adjacent sampling points according to the calculated phase angles; 相位折叠去除模块,用于根据计算得到的所述相邻采样点之间的相位差进行校正,以去除相位折叠,得到校正后的相位差序列;a phase folding removing module, configured to perform correction according to the calculated phase difference between the adjacent sampling points to remove the phase folding and obtain a corrected phase difference sequence; 频率偏差计算模块,用于根据所述校正后的相位差序列计算载波频率偏差输出。The frequency deviation calculation module is used for calculating the carrier frequency deviation output according to the corrected phase difference sequence. 8.根据权利要求7所述的载波频率偏差估计装置,其特征在于,所述相位折叠去除模块与所述频率偏差计算模块之间还连接有低通滤波模块,以用于对所述校正后的相位差序列进行低通滤波,得到滤波后信号输出。8. The carrier frequency deviation estimation device according to claim 7, characterized in that, a low-pass filtering module is also connected between the phase folding removal module and the frequency deviation calculation module, so as to be used for the correction after the correction. The phase difference sequence is low-pass filtered to obtain the filtered signal output. 9.一种计算机设备,包括处理器以及存储器,所述存储器用于存储计算机程序,其特征在于,所述处理器用于执行所述计算机程序以执行如权利要求1~6中任意一项所述方法。9. A computer device, comprising a processor and a memory, wherein the memory is used to store a computer program, wherein the processor is used to execute the computer program to execute the computer program according to any one of claims 1 to 6 method. 10.一种存储有计算机程序的计算机可读存储介质,其特征在于,所述计算机程序执行时实现如权利要求1~6中任意一项所述的方法。10 . A computer-readable storage medium storing a computer program, wherein, when the computer program is executed, the method according to any one of claims 1 to 6 is implemented. 11 .
CN202210699989.4A 2022-06-20 2022-06-20 A carrier frequency deviation estimation method and device Active CN115118564B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210699989.4A CN115118564B (en) 2022-06-20 2022-06-20 A carrier frequency deviation estimation method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210699989.4A CN115118564B (en) 2022-06-20 2022-06-20 A carrier frequency deviation estimation method and device

Publications (2)

Publication Number Publication Date
CN115118564A true CN115118564A (en) 2022-09-27
CN115118564B CN115118564B (en) 2023-12-29

Family

ID=83328726

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210699989.4A Active CN115118564B (en) 2022-06-20 2022-06-20 A carrier frequency deviation estimation method and device

Country Status (1)

Country Link
CN (1) CN115118564B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115801518A (en) * 2022-11-30 2023-03-14 湖南智领通信科技有限公司 Frequency offset estimation method and device based on probability distribution statistics and computer equipment
CN116389210A (en) * 2023-03-21 2023-07-04 湖南艾科诺维科技有限公司 Blind Estimation Method, System and Medium of FSK Carrier Frequency Offset Based on L&R Algorithm

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020057733A1 (en) * 2000-08-09 2002-05-16 Sullivan Mark C. Frequency translator using a cordic phase rotator
US6567480B1 (en) * 1999-08-10 2003-05-20 Lucent Technologies Inc. Method and apparatus for sampling timing adjustment and frequency offset compensation
US6771715B1 (en) * 2000-03-30 2004-08-03 Adtran, Inc. Demodulator using cordic rotator-based digital phase locked loop for carrier frequency correction
CN1549622A (en) * 2003-05-23 2004-11-24 乐金电子(中国)研究开发中心有限公 Method and apparatus for carrier deviation estimation in mobile communication system
US20060176802A1 (en) * 2005-02-04 2006-08-10 Samsung Electronics Co., Ltd. Apparatus and method for compensating for frequency offset in wireless communication system
US20090135978A1 (en) * 2007-11-28 2009-05-28 Chi-Tung Chang Apparatus and method for estimating and compensating sampling frequency offset
US7702040B1 (en) * 2006-04-12 2010-04-20 Sirf Technology, Inc. Method and apparatus for frequency discriminator and data demodulation in frequency lock loop of digital code division multiple access (CDMA) receivers
CN102223344A (en) * 2011-07-20 2011-10-19 四川虹微技术有限公司 Method and system for estimating residual fractional frequency offset
CN102647381A (en) * 2012-03-22 2012-08-22 中兴通讯股份有限公司 Method and device for elevating frequency deviation in MPSK (Mary Phase Shift Keying) coherent light communication system
CN104601512A (en) * 2015-01-05 2015-05-06 江苏指南针导航通信技术有限公司 Method and system for detecting carrier frequency offset of phase-modulated signals
CN106411802A (en) * 2016-09-06 2017-02-15 北京理工大学 Frequency offset correction method based on unmanned aerial vehicle data link transmission
CN106559105A (en) * 2016-11-05 2017-04-05 北京晓程科技股份有限公司 Broadband power carrier OFDM system sampling frequency deviation is estimated and bearing calibration
CN112291013A (en) * 2020-09-25 2021-01-29 湖南工业大学 A Blind Estimation Method of Sampling Clock Frequency Offset for Optical OFDM Systems

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6567480B1 (en) * 1999-08-10 2003-05-20 Lucent Technologies Inc. Method and apparatus for sampling timing adjustment and frequency offset compensation
US6771715B1 (en) * 2000-03-30 2004-08-03 Adtran, Inc. Demodulator using cordic rotator-based digital phase locked loop for carrier frequency correction
US20020057733A1 (en) * 2000-08-09 2002-05-16 Sullivan Mark C. Frequency translator using a cordic phase rotator
CN1549622A (en) * 2003-05-23 2004-11-24 乐金电子(中国)研究开发中心有限公 Method and apparatus for carrier deviation estimation in mobile communication system
US20060176802A1 (en) * 2005-02-04 2006-08-10 Samsung Electronics Co., Ltd. Apparatus and method for compensating for frequency offset in wireless communication system
US7702040B1 (en) * 2006-04-12 2010-04-20 Sirf Technology, Inc. Method and apparatus for frequency discriminator and data demodulation in frequency lock loop of digital code division multiple access (CDMA) receivers
US20090135978A1 (en) * 2007-11-28 2009-05-28 Chi-Tung Chang Apparatus and method for estimating and compensating sampling frequency offset
CN102223344A (en) * 2011-07-20 2011-10-19 四川虹微技术有限公司 Method and system for estimating residual fractional frequency offset
CN102647381A (en) * 2012-03-22 2012-08-22 中兴通讯股份有限公司 Method and device for elevating frequency deviation in MPSK (Mary Phase Shift Keying) coherent light communication system
CN104601512A (en) * 2015-01-05 2015-05-06 江苏指南针导航通信技术有限公司 Method and system for detecting carrier frequency offset of phase-modulated signals
CN106411802A (en) * 2016-09-06 2017-02-15 北京理工大学 Frequency offset correction method based on unmanned aerial vehicle data link transmission
CN106559105A (en) * 2016-11-05 2017-04-05 北京晓程科技股份有限公司 Broadband power carrier OFDM system sampling frequency deviation is estimated and bearing calibration
CN112291013A (en) * 2020-09-25 2021-01-29 湖南工业大学 A Blind Estimation Method of Sampling Clock Frequency Offset for Optical OFDM Systems

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
3GPP: ""GP-052574 Feasibility study 45912 V030"", 3GPP SPECS\\45_SERIES *
JUN WOO KIM, YONG SU LEE: "Carrier Frequency Offset Estimation for OFDM System with Large Oscillator Phase Noise", 2021 INTERNATIONAL CONFERENCE ON INFORMATION AND COMMUNICATION TECHNOLOGY CONVERGENCE *
唐比: "基于光频梳的相干光通信系统中联合载波恢复研究", 中国优秀硕士学位论文数据库 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115801518A (en) * 2022-11-30 2023-03-14 湖南智领通信科技有限公司 Frequency offset estimation method and device based on probability distribution statistics and computer equipment
CN115801518B (en) * 2022-11-30 2024-04-09 湖南智领通信科技有限公司 Frequency offset estimation method and device based on probability distribution statistics and computer equipment
CN116389210A (en) * 2023-03-21 2023-07-04 湖南艾科诺维科技有限公司 Blind Estimation Method, System and Medium of FSK Carrier Frequency Offset Based on L&R Algorithm

Also Published As

Publication number Publication date
CN115118564B (en) 2023-12-29

Similar Documents

Publication Publication Date Title
US9270390B2 (en) Frequency and phase offset compensation of modulated signals with symbol timing recovery
US20050259768A1 (en) Digital receiver and method for processing received signals
CN110912847B (en) GMSK signal demodulation method
CN115118564B (en) A carrier frequency deviation estimation method and device
US8571139B2 (en) Method for tracking the phase of a signal modulated by continuous phase modulation and synchronization device implementing said method
WO2001082546A1 (en) Blind carrier offset detection for quadrature modulated digital communication systems
CN106603217B (en) Sampling frequency offset suppression method for Bluetooth signal of wireless comprehensive tester
CN108353066B (en) Device and method for carrier frequency offset correction and storage medium thereof
CN111935046A (en) Low-complexity frequency shift keying signal symbol rate estimation method
CN107579942B (en) Demodulation method of APSK (amplitude phase Shift keying) modulation signal
CN105049174A (en) Carrier and clock combined synchronization method for OQPSK modulation
Gudovskiy et al. A novel nondata-aided synchronization algorithm for MSK-type-modulated signals
CN108011649B (en) Construction Method of Symbol Synchronization Circuit in Bluetooth EDR Receiver
KR100760793B1 (en) Correction of Orthogonal and Gain Errors in Synchronous Receivers
CN114499646B (en) VDE-SAT channel physical layer signal time domain pilot frequency equalization method and equipment
US7778359B2 (en) Method and device for clock synchronisation with a vestigial-sideband-modulated transmitted signal
TW201631939A (en) Maximum likelihood sequence detection in the phase domain
CN115603888B (en) A method for estimating intra-symbol delay of wireless communication system and recovering signal
US7869534B2 (en) Method and apparatus for clock-pulse synchronization with an offset QPSK-modulated transmission signal
JP3873078B2 (en) Timing extracting apparatus and method, and demodulating apparatus including the timing extracting apparatus
CN116073979B (en) A symbol synchronization method and system for a multi-ary differential phase shift keying system
CN115632923A (en) OQPSK (offset quadrature phase shift keying) -based unmanned aerial vehicle and satellite ultra-wideband communication method and related equipment
JP4335125B2 (en) Timing synchronization circuit
CN106789797B (en) A Symbol Synchronization Method for PSK Signals in Wireless Burst Communication
JPH02230845A (en) Psk signal demodulation system

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