CN106101039B - One kind assisting the adjustable frequency deviation estimating method of precision based on data - Google Patents

One kind assisting the adjustable frequency deviation estimating method of precision based on data Download PDF

Info

Publication number
CN106101039B
CN106101039B CN201610398391.6A CN201610398391A CN106101039B CN 106101039 B CN106101039 B CN 106101039B CN 201610398391 A CN201610398391 A CN 201610398391A CN 106101039 B CN106101039 B CN 106101039B
Authority
CN
China
Prior art keywords
function
frequency deviation
algorithm
sebolic addressing
symbol sebolic
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
CN201610398391.6A
Other languages
Chinese (zh)
Other versions
CN106101039A (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.)
Peking University
Original Assignee
Peking University
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 Peking University filed Critical Peking University
Priority to CN201610398391.6A priority Critical patent/CN106101039B/en
Publication of CN106101039A publication Critical patent/CN106101039A/en
Application granted granted Critical
Publication of CN106101039B publication Critical patent/CN106101039B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0024—Carrier regulation at the receiver end
    • H04L2027/0026—Correction of carrier offset

Abstract

The invention discloses one kind to assist the adjustable frequency deviation estimating method of precision based on data.The method include the steps that by the angle difference of the adjacent x R function of calculating come estimation frequency deviation；Wherein, R function is the function of the receiver carrying angular deflection signal that symbol sebolic addressing is calculated based on the received.Calculate the formula of the exemplary frequency deviation values are as follows:2≤x≤N-1 and be even number；3≤x≤N-1 is odd number；1≤m≤N, N are the positive integer no more than L-1.This method has better accuracy under the conditions of flat fading channel, than M&M algorithm, has bigger estimation range than L&R algorithm, Fitz algorithm.

Description

One kind assisting the adjustable frequency deviation estimating method of precision based on data
Technical field
The invention belongs to digital communication technology fields, and in particular to digital receiver docking is collected mail in carrier wave communication system Number frequency offset estimation.
Background technique
Nonlinear Transformation in Frequency Offset Estimation is about one of the main problem in digital receiver research, and common method is using dedicated lock Phase ring core piece construct hardware module method, but hardware approach there are reaction time length, Measurement bandwidth relative narrower, frequency with The problems such as track accuracy and unstable precision.More than the method that numeric field realizes offset estimation can be solved the problems, such as effectively, Wherein by whether training sequence being needed to assist, data auxiliary (data-aided, DA) and unbound nucleus (non-can be divided into Data-aided, NDA), under data subsidiary conditions classical algorithm include maximum likelihood function method, L&R algorithm, Fitz algorithm, M&M algorithm etc., they can be applied under Gaussian white noise channel and flat fading channel conditions.Wherein maximum likelihood function method Calculate more complex, then 3 kinds of algorithms have preferable engineering application value.Wherein, 3 kinds of algorithm essences of this under the conditions of awgn channel Degree is very close, and the estimation range of M&M algorithm is much larger than other two kinds.But under the conditions of Rayleigh channel, tied according to emulation Fruit, M&M algorithm estimated accuracy substantially decay compared with other two algorithm.
The present invention is directed to above situation, has carried out corresponding modification to M&M algorithm, modified algorithm reduces estimation appropriate Range but still in advantageous situation, adjusting parameter can obtain more preferable estimated accuracy according to actual needs.
Summary of the invention:
For the technical problems in the prior art, the purpose of the present invention is to provide one kind assists precision based on data Adjustable frequency deviation estimating method.
The principle of the data-aided frequency excursion algorithms such as L&R algorithm, Fitz algorithm, M&M algorithm is converted in A/D Afterwards, it is handled by carrying out the autocorrelative mathematical operation of class to reception sequence known to symbol, to estimate that angle caused by frequency deviation is inclined Turn value.
Assuming that sending the symbol sebolic addressing C that length is Lk={ c1,c2,...,ck,...,cn, and
Then receiving symbol is x (k)=ckej(2πfkT+θ)+ n (k), wherein ckFor k-th of sequence of symhols of transmission, n (k) is to make an uproar Sound item.
The function that class auto-correlation computation obtains carrying angular deflection signal is done to symbol is received, does class auto-correlation fortune twice It calculates, obtains
Wherein 1≤m≤N
Wherein N is the positive integer no more than L-1, then the angle value of R function is the angle information for carrying frequency deflection, right R function asks angle to have
arg{R(m)}≈[2πmfT+λ(m)]2π
Wherein f is required frequency values, and λ (m) is noise item, and L&R algorithm, Fitz algorithm are by calculating N number of R (m) The angle value of function simultaneously does appropriate weighted average and does to frequency departure and estimate.
L&R algorithm:
Fitz algorithm:
M&M algorithm is by calculating the angle value of two adjacent R functions and doing smoothing processing to frequency by window function Deviation is estimated.
M&M algorithm:
T is to receive symbol period in algorithm expression formula above, and arg { f (x) } is to seek f (x) angle in [0,2 π] range Operation, w (m) are smooth window function, and the angle difference that M&M algorithm calculates two neighboring R function includes single f information；And it is another Two kinds of algorithms are up to N number of f information, and the angle value range that can be estimated is determined by N, therefore, the theoretical estimation model of Fitz algorithm Enclose forL&R algorithm isAnd M&M algorithm isM&M algorithm estimation range is much larger than other two kinds.
But since single f information is less, thus be also easier to it is affected by noise, under the conditions of flat fading channel, M&M calculate Method precision is poor.Therefore, the present invention proposes the modification of algorithm, and frequency is estimated by seeking the angle difference of adjacent x R function Rate deviation, expression formula are as follows:
2≤x≤N-1 and be even number
3≤x≤N-1 is odd number
In above formulaAs estimate that frequency deviation value, x are positive integer and 2≤x≤N-1, the theoretical estimation range of above formula is It is still greater than L&R algorithm and Fitz algorithm.
BecauseIt can be seen that identical In the case where f information, modified algorithm is affected by noise smaller.
Compared with prior art, the positive effect of the present invention are as follows:
Under the conditions of flat fading channel, there is better accuracy than M&M algorithm, have more than L&R algorithm, Fitz algorithm Big estimation range.
The value that simulation result as shown in Figure 2,3 shows by adjusting x, can be the case where suitably reducing estimation range Lower acquisition better estimated accuracy under flat fading Rayleigh channel.X is bigger, and estimated accuracy is better, but estimation range subtracts therewith It is small；When x is maximized (N-1), algorithm estimation range and precision are close to L&R algorithm.
Detailed description of the invention
Fig. 1 is program circuit schematic diagram；
Fig. 2 is when x takes 3, and simulation estimate of the three kinds of algorithms of this algorithm and other under flat fading Rayleigh channel is square Poor effect is with signal-to-noise ratio variation diagram.
Fig. 3 is when x takes N-1, and simulation estimate of the three kinds of algorithms of this algorithm and other under flat fading Rayleigh channel is equal Variance effect is with signal-to-noise ratio variation diagram.
Specific embodiment
1, assume that sending the symbol sebolic addressing that length is L is Ck={ c1,c2,...cn, and
2, receiver A/D is represented by x (k)=c to the sampled symbols sequence obtained after the symbol sebolic addressing conversion receivedkej (2πfkT+θ)+ n (k), wherein k is symbol ordinal number, ckFor the sequence of symhols in the symbol sebolic addressing of transmission, f is frequency deviation value to be estimated, and T is Symbol period, θ are symbol-modulated angle, and n (k) is the noise item of k-th of reception symbol.
3, the R function that class auto-correlation computation twice obtains carrying angular deflection signal is done to sampled symbols sequence x (k), such as Under
Wherein 1≤m≤N
Wherein m, N are positive integer less than L, according to theory analysis and simulation result,When estimation effect it is best, Therefore it takes
5, estimated by calculating R function angle value and doing appropriate weighted average and can be done to frequency departure, the present invention proposes Algorithm arrangement, expression formula is as follows:
2≤x≤N-1 and be even number
3≤x≤N-1 is odd number
In above formulaAs estimate that frequency deviation value, x are positive integer and 2≤x≤N-1, theoretical estimation range isIt will connect Symbol is received to removeCaused by angular deflection, realize Frequency Synchronization.

Claims (2)

1. one kind assists the adjustable frequency deviation estimating method of precision based on data, which is characterized in that by calculating adjacent x R function Angle difference carry out estimation frequency deviation；Wherein, R function is the receiver carrying that symbol sebolic addressing is calculated based on the received The function of angular deflection signal；
Calculate the formula of the exemplary frequency deviation values are as follows: It and is even number；For odd number；Wherein,For the exemplary frequency deviation values, x is positive integer；L is symbol sebolic addressing Length, T is the period of symbol sebolic addressing, and 1≤m≤N, N are positive integer no more than L-1；
The method for obtaining the R function are as follows: 1) receiver obtains sampled symbols sequence after converting to the symbol sebolic addressing received；Its In, symbol sebolic addressing Ck={ c1,c2,...ck... cn, andSampled symbols sequence is x (k)=ckej(2πfkT+θ)+n (k), ckFor k-th of sequence of symhols of transmission, f is frequency deviation value to be estimated, and θ is symbol-modulated angle, and n (k) is k-th of reception symbol Noise item；2) class auto-correlation computation is done to the sampled symbols sequence according to the symbol sebolic addressing, obtains carrying angular deflection signal Function；3) class auto-correlation computation construction R function is done to the function:
2. the method as described in claim 1, which is characterized in that
CN201610398391.6A 2016-06-07 2016-06-07 One kind assisting the adjustable frequency deviation estimating method of precision based on data Active CN106101039B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610398391.6A CN106101039B (en) 2016-06-07 2016-06-07 One kind assisting the adjustable frequency deviation estimating method of precision based on data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610398391.6A CN106101039B (en) 2016-06-07 2016-06-07 One kind assisting the adjustable frequency deviation estimating method of precision based on data

Publications (2)

Publication Number Publication Date
CN106101039A CN106101039A (en) 2016-11-09
CN106101039B true CN106101039B (en) 2019-01-25

Family

ID=57227340

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610398391.6A Active CN106101039B (en) 2016-06-07 2016-06-07 One kind assisting the adjustable frequency deviation estimating method of precision based on data

Country Status (1)

Country Link
CN (1) CN106101039B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107342960B (en) * 2016-11-29 2019-12-31 重庆邮电大学 Non-data-aided frequency offset estimation method suitable for amplitude phase shift keying
CN107204782B (en) * 2017-04-10 2020-11-20 北京大学 BCH decoder and implementation method of compiler for generating BCH decoder
CN109600330B (en) * 2018-11-23 2021-12-14 周口师范学院 Simplified diagonal cross-correlation carrier frequency offset estimation method
CN113556304B (en) * 2021-06-02 2022-11-04 北京大学 Time-varying frequency offset estimation method, system and medium based on particle filter

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101854229A (en) * 2010-05-14 2010-10-06 中国人民解放军理工大学 Iteration demodulation decoding method of encoded modulation signals based on climax frequency deviation compensation
CN102904843A (en) * 2012-08-02 2013-01-30 京信通信系统（广州）有限公司 Frequency offset estimation method and device
US8755462B2 (en) * 2010-09-17 2014-06-17 Vecima Networks Inc. Frequency offset estimator for upstream cable signals
CN103916348A (en) * 2012-12-30 2014-07-09 重庆重邮信科通信技术有限公司 Calculation methods and systems for phase deviant, timing deviation and frequency deviation

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100568069B1 (en) * 2004-09-02 2006-04-05 한국전자통신연구원 Apparatus and Method of Feedforward Frequency offset Estimation in TDMA System

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101854229A (en) * 2010-05-14 2010-10-06 中国人民解放军理工大学 Iteration demodulation decoding method of encoded modulation signals based on climax frequency deviation compensation
US8755462B2 (en) * 2010-09-17 2014-06-17 Vecima Networks Inc. Frequency offset estimator for upstream cable signals
CN102904843A (en) * 2012-08-02 2013-01-30 京信通信系统（广州）有限公司 Frequency offset estimation method and device
CN103916348A (en) * 2012-12-30 2014-07-09 重庆重邮信科通信技术有限公司 Calculation methods and systems for phase deviant, timing deviation and frequency deviation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
An Improved Frequency Offset Estimator for OFDM Applications;Michele Morelli等;《IEEE Communications Letters》;19990331;全文
Feedforward Frequency Estimation for PSK: a Tbtorial Review;MICHELE MOREL等;《Special Issue》;19980430;正文第5.1、5.4、6.5、6.7节

Also Published As

Publication number Publication date
CN106101039A (en) 2016-11-09

Similar Documents

Publication Publication Date Title
CN106101039B (en) One kind assisting the adjustable frequency deviation estimating method of precision based on data
CN104320361B (en) A kind of serial separation method of spaceborne ais signal based under parameter Estimation
JP4979413B2 (en) Pulse radio receiver
CN108471347A (en) A kind of parallel time synchronization method based on accurate loop filtering
JP4315658B2 (en) Method for decoding incoming pulse signal of ultra-wideband type
US20140118066A1 (en) Low-cost digital predistortion apparatus and method using envelope detection feedback
CN103428153B (en) GMSK signal acceptance method in satellite mobile communication
CN111884685A (en) Synchronous demodulation method and device for digital communication signal
CN109274623B (en) Implementation method for correcting carrier frequency offset based on large-point FFT (fast Fourier transform)
CN111935046B (en) Low-complexity frequency shift keying signal symbol rate estimation method
CN107342960B (en) Non-data-aided frequency offset estimation method suitable for amplitude phase shift keying
CN103338024B (en) The complementary Kalman filtering apparatus and method of time delay in antenna array
CN107483130A (en) One kind joint broader frequency spectrum perceives and angle-of- arrival estimation method
CN105162570A (en) Timing synchronization method and device for signal parallel processing
CN107623647B (en) Carrier synchronization method based on scattered pilot frequency assistance
CN103916199A (en) Device and method for time delay and phase adjustment of antenna signal
US9450745B2 (en) Method and apparatus for radio frequency (RF) pulse synchronization in super regenerative receiver (SRR)
CN102035769B (en) Phase shift keying signal blind detection method based on plural discrete full-feedback neural network
CN105933265B (en) A kind of pair of QAM signal carries out the phase noise blind estimating method of unbound nucleus
CN106026974B (en) A kind of passband global response error constraints airspace matrix filter design method
KR20210103823A (en) Clock and data recovery device using machine learning technology and training method thereof
CN107045121A (en) A kind of near field ultra-broadband signal phase difference ranging method and system
CN110445740A (en) Frequency deviation estimating method and system based on repetitive sequence
CN106341199A (en) Signal to noise ratio determination method based on Manchester code
CN108401582B (en) A kind of GEO satellite mobile communication system initial frequency synchronization method adapting to TDMA/FDMA systems

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