CN111865850A - PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection - Google Patents

PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection Download PDF

Info

Publication number
CN111865850A
CN111865850A CN202010660213.2A CN202010660213A CN111865850A CN 111865850 A CN111865850 A CN 111865850A CN 202010660213 A CN202010660213 A CN 202010660213A CN 111865850 A CN111865850 A CN 111865850A
Authority
CN
China
Prior art keywords
branch
late
early
frequency
pcm
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
CN202010660213.2A
Other languages
Chinese (zh)
Other versions
CN111865850B (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.)
Institute of Electronic Engineering of CAEP
Original Assignee
Institute of Electronic Engineering of CAEP
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 Institute of Electronic Engineering of CAEP filed Critical Institute of Electronic Engineering of CAEP
Priority to CN202010660213.2A priority Critical patent/CN111865850B/en
Publication of CN111865850A publication Critical patent/CN111865850A/en
Application granted granted Critical
Publication of CN111865850B publication Critical patent/CN111865850B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/0044—Control loops for carrier regulation
    • H04L2027/0063—Elements of loops
    • H04L2027/0069—Loop filters

Abstract

The invention belongs to the technical field of broadband signal sampling, and particularly relates to a PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection.

Description

PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection
Technical Field
The invention belongs to the technical field of signal processing, and particularly relates to a PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection.
Background
In a wireless communication scene with large dynamic, the received PCM/FM signals have large Doppler frequency offset and first-order and second-order Doppler change rates, and show large dynamic characteristics. The demodulation performance of the telemetry signal is affected by the excessive doppler frequency offset, so the receiving end usually needs to estimate the doppler frequency, complete the capture and tracking of the carrier frequency, and then perform subsequent signal demodulation.
The carrier frequency synchronization of existing PCM/FM telemetry signals typically uses estimation methods based on power spectrum characteristics, such as the power spectrum maxima method, the power spectrum medians method, and the power spectrum centroids method. The first two algorithms complete the estimation of the carrier frequency by solving the frequency corresponding to the maximum amplitude of the power spectrum or estimating the center frequency of the power spectrum in the signal bandwidth after performing FFT on the signal, but the estimation accuracy of the two methods is poor because the maximum amplitude of the power spectrum and the estimation of the center frequency are greatly affected by noise. Especially when the signal-to-noise ratio is low, the method is difficult to realize high-performance carrier synchronization; the third algorithm, the power spectrum gravity center method, utilizes the symmetry of the signal power spectrum to realize the estimation of the carrier frequency. When the modulation code elements '0' and '1' at the transmitting end keep relatively balanced quantity, the center of gravity of the PCM/FM signal power spectrum is always at the carrier frequency, and then the algorithm can realize the carrier frequency synchronization with higher precision. In addition, due to the fact that the utilized power spectrum information is increased, the anti-noise performance of the algorithm is greatly improved compared with the former two algorithms, but unfortunately, in practical application, for most PCM/FM transmitters, due to the non-randomness of the data frame structure, the number of modulation symbols of the transmitter is not balanced, so that the center of gravity of the signal power spectrum is not always at the carrier frequency, which causes the estimation deviation of the synchronization method, and therefore, several carrier frequency estimation methods based on the power spectrum characteristics have certain limitations in practical application.
Besides the estimation method based on frequency domain characteristics, the carrier frequency synchronization of PCM/FM signals also has estimation based on the frequency discrimination demodulation resultThe method, as the Chinese patent publication with publication number CN104486288A, which discloses a method for suppressing carrier frequency offset applicable to PCM/FM telemetering receiver with the publication time of 2015, 4.1.A method for suppressing carrier frequency offset applicable to PCM/FM telemetering receiver is disclosed, which utilizes the statistical characteristics of formed code elements, performs histogram estimation through frequency-discriminated signals to realize modulation frequency estimation corresponding to positive and negative code elements, further averages the modulation frequency estimation values to obtain an estimated value of carrier frequency offset, subtracts the estimated value of carrier frequency offset from the frequency-discriminated signals to achieve the effect of suppressing carrier frequency offset, is less influenced by noise, is applicable to the conditions of large frequency offset and unequal numbers of positive and negative code elements, has small calculated amount, is simple to implement and is easy to implement in real time, and the method obtains the quantity of '0' and '1' through statistical judgment in a period of time, estimating carrier frequency position f corresponding to modulation symbol ' 0 ' and symbol ' 10And f1And by looking for f0And f1The method solves the problem of imbalance of modulation code elements '0' and '1' of the transmitter; however, the estimation accuracy for the decision depends on the accuracy of the decision result, and the decision result of the frequency discrimination demodulation of the PCM/FM signal under the low signal-to-noise ratio is low in reliability, so that the frequency estimation accuracy under the low signal-to-noise ratio is affected. In addition, in engineering practice, after a PCM/FM signal passes through an actual radio frequency channel and a wireless channel, due to the amplitude-frequency response or filter design and the like, it is difficult to ensure that the amplitudes of frequency discrimination results corresponding to positive and negative code elements are equivalent, so that a certain residual frequency offset exists when the method is used for carrier frequency synchronization.
The carrier frequency synchronization based on the early-late gate is the most classical signal synchronization structure, and the signal synchronization of a plurality of modulation systems can realize the carrier frequency synchronization by adopting the thought of the early-late gate. The traditional early-late gate synchronization makes the output results of the early branch and the late branch different, and the carrier frequency estimation is finished iteratively by judging the polarity of the difference and selecting the step of adjusting the frequency towards the positive direction or the negative direction, but the step usually determines the speed of frequency adjustment and the estimation precision of a feedback loop, namely the larger step can realize the faster frequency offset capture, but is not beneficial to the fine tracking of the carrier; a smaller step may enable higher accuracy carrier frequency tracking, but may make the acquisition process in the earlier stage too lengthy.
The adoption of the step-size-changing processing of 'big first and small second' is a concept for solving the problems of acquisition speed and tracking accuracy of an early-late gate, but the complexity of loop control is increased, and particularly when signal interruption needs rapid recapture, the control of carrier synchronization becomes complicated and even the re-compensation fails.
In summary, carrier frequency synchronization of PCM/FM signals needs to overcome the problems of low signal-to-noise ratio, unbalanced modulation symbols "0" and "1", fast acquisition and tracking of large dynamic doppler frequency, and the like, but the current algorithms do not solve these problems well.
Disclosure of Invention
The invention aims to solve the defects and shortcomings in the prior art, realize frequency error estimation by utilizing multi-symbol detection, and realize the PCM/FM signal early-late loop frequency synchronization method for quickly capturing and tracking the carrier frequency of a large dynamic PCM/FM signal by adopting an early-late loop carrier synchronization structure added with a loop filter.
In order to achieve the above object, the present invention discloses a PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection, which is characterized by comprising the following steps:
step 1, according to the estimated value of the frequency deviation of the current carrier wave
Figure BDA0002578207730000031
Frequency interval value of sum branch and late branchfCalculating a frequency deviation value, respectively controlling Numerically Controlled Oscillators (NCO) of an early branch (an advanced branch) and a late branch (a lagging branch) according to the frequency deviation value, and respectively carrying out frequency correction on a received signal (a complex baseband after down-conversion) s (t) to obtain an early branch signal sel(t) and signal s of the late branchlt(t)；
The frequency deviation values for controlling the numerical controlled oscillators of the early branch (leading branch) and the late branch (lagging branch) are respectively
Figure BDA0002578207730000032
And
Figure BDA0002578207730000033
i.e. the early branch signal sel(t) late branch signal
Figure BDA0002578207730000034
slt(t) is
Figure BDA0002578207730000035
The frequency interval value fThe method is optimized and selected through an S curve of frequency error estimation, and the selection principle is that the larger the slope at the origin (frequency offset is 0) of the S curve, the higher the estimation accuracy is, so that the slope of the S curve at the origin should be selected as much as possiblef(ii) a In addition, since the intersection (also referred to as the second zero-crossing point) of the S-curve with the frequency offset axis (amplitude of 0) at a position other than the origin determines the fast capture band of the loop, the second zero-crossing point (farther from the origin) should be selected to allow a large loop capture rangefCombining the above principles, selecting for PCM/FM signalsf0.16Rb (Rb is the symbol rate).
Step 2, the receiver receives the ith PCM/FM modulation symbol sent by the transmitter currently, and then the receiver respectively carries out incoherent multi-symbol detection on the signals of the early branch and the late branch to obtain the module value of the matched filtering output result;
the noncoherent multi-symbol detection is to combine the received signals in the ith symbol interval of the two branches of the early branch and the late branch with L groups of local waveforms (matched filter)
Figure BDA0002578207730000036
L is conjugate correlated symbol by symbol to obtain a modulus of the matched filter output
Figure BDA0002578207730000037
The matched filter is a fixed sampling point in each symbol interval and corresponds to the sampling point of the received signal one by one;
Figure BDA0002578207730000038
Figure BDA0002578207730000039
Wherein "s (t)" is the conjugate of the complex number s (t),
Figure BDA0002578207730000041
is a matched filter.
The matched filter
Figure BDA0002578207730000042
Is NxT, T is the symbol period of the PCM/FM signal, N is the number of symbols required for one detection, N is>1；
Then, the ith local waveform is:
Figure BDA0002578207730000043
Figure BDA0002578207730000044
wherein the content of the first and second substances,
Figure BDA0002578207730000045
for N symbols to be selected
Figure BDA0002578207730000046
(i is 1,2, …, N and
Figure BDA0002578207730000047
) The constituent symbol vectors, i.e.
Figure BDA0002578207730000048
Kf is the modulation index of PCM/FM and takes 0.7, g (t) is the pre-adjusting filter of PCM/FM signal, usually takes the rectangular shaping filter, that is
Figure BDA0002578207730000049
The matched filter
Figure BDA00025782077300000410
The number L of (A) is in a relationship of L-2 with NN。
Preferably, in the present invention, N is 3, so L is 8, and the number of matched filters is reduced by 75% compared to N5 (L is 32) which is commonly used for incoherent detection.
Step 3, respectively solving the maximum value of the N MSD measurements obtained by the early branch and the late branch to obtain the maximum measurement of the two branches
Figure BDA00025782077300000411
Every M symbols, respectively calculating the average value of the maximum measurement of the two branches
Figure BDA00025782077300000412
The symbol number M of the accumulation interval determines the granularity of carrier frequency estimation, the larger M is, the more symbols are used for each estimation, so that the higher the precision of single estimation is, the smaller the estimation jitter is, but the larger M can increase the time delay of loop feedback, thereby slowing down the loop capturing speed. In addition, doppler acceleration may cause certain frequency offset changes to occur within multiple symbols, and thus a larger M may reduce the ability of the loop to adapt to doppler acceleration. On the other hand, the smaller M, the larger the jitter of a single estimation, so that an excessively small value of M may reduce the estimation accuracy of the loop in the tracking phase. Therefore, M needs to be chosen according to the specific requirements of the communication system for acquisition time and estimation accuracy.
Step 4, taking the difference of the maximum measurement cumulant to obtain the estimation result of the carrier frequency error of the time
e(k)＝λEL(k)-λLT(k) Passing e (k) through the second order shown in FIG. 3Loop filter for obtaining the frequency estimation result
Figure BDA0002578207730000051
The second-order loop filter can track the Doppler frequency with the second-order change rate not being 0, wherein the loop filter parameters K1 and K2 can be selected according to the system requirement;
step 5, using the frequency estimation result in step 4
Figure BDA0002578207730000052
And updating the carrier frequency offset estimation value in the step 1, and repeating the steps 1-5 to complete the early-late loop frequency synchronization of the feedback closed loop structure.
Compared with the prior art, the technical scheme of the invention utilizes the two-path multi-symbol detector to respectively calculate the maximum measurement for the received signals subjected to different frequency offsets, and uses the measurement difference of the two branches as the estimated value of the carrier frequency error so as to drive the frequency-locked loop comprising the second-order loop filter, thereby completing the fast capture and high-precision tracking of the carrier frequency offset.
In addition, in the scheme, the multi-symbol detector for carrier frequency estimation adopts simplified 3-symbol matched filtering, compared with a 5-symbol MSD detector commonly used for incoherent demodulation, the scale of a matched filter is reduced by 75%, the complexity of hardware realization is reduced, and the loss of the demodulation performance of a PCM/FM signal can be ignored after the carrier synchronization of the invention is carried out; the frequency interval of the early-late branch is optimized by using the S curve of the frequency error estimation, the estimation precision of the loop can be improved and the loop capture range is considered through the optimized frequency interval, and for PCM/FM signals, the optimized frequency interval is 0.16 times of the symbol rate. The multi-symbol detector for estimating the error adopts the idea of block estimation, and each block calculates the frequency error once by using the average value of continuous M times of multi-symbol detection output, thereby reducing the jitter of frequency error estimation and improving the error estimation precision.
In summary, compared with the prior art, the technical scheme of the invention has the following advantages:
1. the technical scheme of the invention combines a multi-symbol detector, adopts the processing process of matched filtering to realize frequency error estimation, thereby fully utilizing modulation information, and in contrast, the frequency offset estimation algorithm based on the power spectrum characteristic generally takes the received signal as the modulation signal of modulation code elements of '0' and '1', and the like, and does not utilize the modulation information, so that the carrier frequency estimation precision of the invention is higher compared with the estimation algorithm based on the power spectrum characteristic.
2. Another advantage of the present invention is that the carrier frequency estimation is performed by using a multi-symbol detector, and the estimation process detects modulation symbols "0" and "1" separately, so that the detection process is not affected no matter how the symbols of the modulation signal are distributed. Therefore, the invention can solve the problem that the estimation performance is influenced by the asymmetry of the power spectrum caused by the uneven number of the modulation symbols '0' and '1'.
3. The technical scheme of the invention is based on multi-symbol detection and demodulation, the multi-symbol detection is always acknowledged to be the algorithm with the best performance for the incoherent detection of PCM/FM signals, and the demodulation performance is improved by about 2.5dB compared with the performance of the algorithm based on frequency discrimination, so that the performance of the multi-symbol detection is better.
4. The technical scheme of the invention adopts the early-late path error estimation combined with the feedback synchronization structure of the loop filter to replace the traditional method that the early-late gate only utilizes the detection error polarity to adjust the frequency stepping (usually fixed stepping) in the forward direction or the reverse direction. The traditional early-late gate takes the polarity of the difference value of the leading branch and the lagging branch, is a 'hard decision' processing mode, and loses part of detection information. The error is processed by the loop filter, and the 'soft decision' information of the error estimation of the early-late branch is utilized, so that the rapid accumulation of the detection error in the capturing stage can be realized, and the capturing speed is accelerated; in the tracking stage, the detection error oscillates around 0, and the jitter of feedback adjustment can be reduced through the loop filter, so that the precision of carrier frequency estimation is improved, and therefore, compared with the carrier frequency synchronization of the traditional early-late gate structure, the method has better frequency acquisition speed and estimation precision.
Drawings
The foregoing and following detailed description of the invention will be apparent when read in conjunction with the following drawings, in which:
FIG. 1 is a schematic circuit diagram of an early-late loop frequency synchronization structure based on multi-symbol detection according to the present invention;
FIG. 2 is a schematic diagram of the frequency error estimation S-curve of different frequency interval values of the early-late branch according to the present invention;
Fig. 3 is a schematic diagram of a second-order loop filter circuit structure suitable for the present invention.
Detailed Description
The technical solutions for achieving the objects of the present invention are further illustrated by the following specific examples, and it should be noted that the technical solutions claimed in the present invention include, but are not limited to, the following examples.
Example 1
As a specific implementation scheme of the invention, as shown in FIG. 1, a high-precision carrier frequency synchronization method of PCM/FM telemetering signals is disclosed, wherein a loop filter is added in an early-late loop carrier synchronization structure, and the loop filter is connected between two ends of the early-late loop carrier synchronization structure, so that carrier frequency estimation and correction of PCM/FM signals under high dynamic and low signal-to-noise ratio can be realized, and the method belongs to the field of wireless communication. Specifically, the method comprises the following steps:
step 1, according to the estimated value of the frequency deviation of the current carrier wave
Figure BDA0002578207730000061
Frequency interval value of sum branch and late branchfCalculating a frequency deviation value, respectively controlling Numerically Controlled Oscillators (NCO) of an early branch (an advanced branch) and a late branch (a lagging branch) according to the frequency deviation value, and respectively carrying out frequency correction on a received signal (a complex baseband after down-conversion) s (t) to obtain an early branch signal sel(t) and signal s of the late branch lt(t)；
Step 2, the receiver receives the ith PCM/FM modulation symbol sent by the transmitter currently, and then the receiver respectively carries out incoherent multi-symbol detection on the signals of the early branch and the late branch to obtain the module value of the matched filtering output result;
step 3, respectively solving the maximum value of the N MSD measurements obtained by the early branch and the late branch to obtain the maximum measurement of the two branches
Figure BDA0002578207730000071
Every M symbols, respectively calculating the average value of the maximum measurement of the two branches
Figure BDA0002578207730000072
The symbol number M of the accumulation interval determines the granularity of carrier frequency estimation, the larger M is, the more symbols are used for each estimation, so that the higher the precision of single estimation is, the smaller the estimation jitter is, but the larger M can increase the time delay of loop feedback, thereby slowing down the loop capturing speed. In addition, doppler acceleration may cause certain frequency offset changes to occur within multiple symbols, and thus a larger M may reduce the ability of the loop to adapt to doppler acceleration. On the other hand, the smaller M, the larger the jitter of a single estimation, so that an excessively small value of M may reduce the estimation accuracy of the loop in the tracking phase. Therefore, M needs to be chosen according to the specific requirements of the communication system for acquisition time and estimation accuracy.
Step 4, taking the difference of the maximum measurement cumulant to obtain the estimation result of the carrier frequency error of the time
e(k)＝λEL(k)-λLT(k) Passing e (k) through a second order loop filter shown in FIG. 3 to obtain the frequency estimation result of this time
Figure BDA0002578207730000073
The second-order loop filter can track the Doppler frequency with the second-order change rate not being 0, wherein the loop filter parameters K1 and K2 can be selected according to the system requirement;
step 5, using the frequency estimation result in step 4
Figure BDA0002578207730000074
Updating the carrier frequency offset estimation value in the step 1, repeating the steps 1-5, and finishing the feedback closed loopThe early-late loop frequency of the structure is synchronized.
Example 2
As a preferred implementation of the present invention, on the basis of the technical solution of the above example 1, further, as shown in fig. 1, it is assumed that the current carrier frequency offset estimation value is
Figure BDA0002578207730000075
The working steps are as follows:
by a frequency offset value
Figure BDA0002578207730000076
And
Figure BDA0002578207730000077
respectively controlling Numerically Controlled Oscillators (NCO) of a leading branch (early branch) and a lagging branch (late branch), respectively carrying out frequency correction on a received signal (complex baseband after down-conversion) s (t) to obtain an early branch signal sel(t) and signal s of the late branchlt(t)：
Figure BDA0002578207730000081
Figure BDA0002578207730000082
Wherein the content of the first and second substances,ffor the frequency interval values of the early branch and the late branch, the parameter is optimally selected by an S-curve of frequency error estimation, and the S-curves at different frequency interval values are shown in fig. 2.
The selection principle is as follows: the larger the slope of the S-curve at the origin (frequency offset is 0), the higher the estimation accuracy, and therefore, the one that results in the larger slope of the S-curve at the origin should be selected as much as possiblef(ii) a In addition, since the intersection (also referred to as the second zero-crossing point) of the S-curve with the frequency offset axis (amplitude of 0) at a position other than the origin determines the fast capture band of the loop, the second zero-crossing point (farther from the origin) should be selected to allow a large loop capture rangef。
Combining the above principles toPCM/FM signal selectionf0.16Rb (Rb is the symbol rate).
Then, at the ith symbol received, the signals of the two branches are respectively subjected to non-coherent Multi-symbol Detection (MSD), i.e. are respectively matched with L groups of local waveforms
Figure BDA0002578207730000083
Figure BDA0002578207730000083
Figure BDA0002578207730000083
1, 2.. L (i.e., matched filter) is complex correlated, obtaining the modulus of the matched filter output:
Figure BDA0002578207730000084
Figure BDA0002578207730000085
where "s x (t)" is the conjugate of the complex number s (t), matched filters
Figure BDA0002578207730000086
Is NxT, T is the symbol period of PCM/FM signal, N is the number of symbols required for one detection, generally N>1。
Then, the ith local waveform is:
Figure BDA0002578207730000087
Figure BDA0002578207730000088
wherein the content of the first and second substances,
Figure BDA0002578207730000089
for N symbols to be selected
Figure BDA00025782077300000810
(i is 1,2, …, N and
Figure BDA00025782077300000811
) Formed symbol vectorI.e. by
Figure BDA00025782077300000812
Kf is the modulation index of PCM/FM and takes 0.7, and g (t) is a rectangular formed pre-adjusting filter
Figure BDA0002578207730000091
The relation between the number L of matched filters and N is L-2N。
In this embodiment, if N is 3, L is 8, and the number of matched filters is reduced by 75% compared to N is 5(L is 32) which is commonly used for incoherent detection.
Further, the maximum values of the N MSD measurements obtained by the early branch and the late branch are respectively obtained to obtain the maximum measurements of the two branches
Figure BDA0002578207730000092
Every M symbols, respectively calculating the average value of the maximum measurement of the two branches
Figure BDA0002578207730000093
The number of accumulated interval symbols M determines the granularity of carrier frequency estimation, and the larger M is, the more symbols are used for each estimation, so that the higher the precision of single estimation is, the smaller the estimation jitter is.
But too large M increases the delay of the loop feedback, slowing down the loop acquisition speed. In addition, doppler acceleration may cause certain frequency offset changes to occur within multiple symbols, and thus a larger M may reduce the ability of the loop to adapt to doppler acceleration.
On the other hand, the smaller M, the larger the jitter of a single estimation, so that an excessively small value of M may reduce the estimation accuracy of the loop in the tracking phase.
Therefore, M needs to be chosen according to the specific requirements of the communication system for acquisition time and estimation accuracy.
The maximum measurement cumulant is subtracted to obtain the estimation result of the carrier frequency error of the time
e(k)＝λEL(k)-λLT(k) E (k) is subjected to loop filtering as shown in FIG. 3A wave filter for obtaining the current frequency estimation result
Figure BDA0002578207730000094
The second order loop filter can track the Doppler frequency with the second order change rate different from 0, wherein the loop filter parameters K1 and K2 can be selected according to the system requirement.
Finally, utilize
Figure BDA0002578207730000095
And updating the carrier frequency offset estimation value in the step 1, and repeating the steps 1-5 to complete the early-late loop frequency synchronization of the feedback closed loop structure.

Claims (5)

1. A PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection is characterized by comprising the following steps:
step 1, according to the estimated value of the frequency deviation of the current carrier wave
Figure RE-FDA0002639611080000011
Frequency interval value of sum branch and late branchfCalculating a frequency deviation value, respectively controlling Numerically Controlled Oscillators (NCO) of an early branch (an advanced branch) and a late branch (a lagging branch) according to the frequency deviation value, and respectively carrying out frequency correction on a received signal (a complex baseband after down-conversion) s (t) to obtain an early branch signal sel(t) and signal s of the late branchlt(t)；
Step 2, the receiver receives the ith PCM/FM modulation symbol sent by the transmitter currently, and then the receiver respectively carries out incoherent multi-symbol detection on the signals of the early branch and the late branch to obtain the module value of the matched filtering output result;
Step 3, respectively solving the maximum value of the N MSD measurements obtained by the early branch and the late branch to obtain the maximum measurement of the two branches
Figure RE-FDA0002639611080000012
Every M symbols, respectively obtaining the maximum measurement of two branchesAverage value of (2)
Figure RE-FDA0002639611080000013
The number of symbols M of the accumulation interval is selected according to the specific requirements of the communication system on acquisition time and estimation precision;
step 4, making difference between the maximum measurement cumulant to obtain the estimation result e (k) lambda of the carrier frequency error of this timeEL(k)-λLT(k) E (k) is processed by the second order loop filter to obtain the frequency estimation result of this time
Figure RE-FDA0002639611080000014
The selection of the loop filter parameters K1 and K2 can be set according to the system requirements;
step 5, using the frequency estimation result in step 4
Figure RE-FDA0002639611080000015
And updating the carrier frequency offset estimation value in the step 1, and repeating the steps 1-5 to complete the early-late loop frequency synchronization of the feedback closed loop structure.
2. The method of claim 1 for multi-symbol detection based early-late loop frequency synchronization of PCM/FM signals, wherein: the frequency deviation values of the numerical control oscillators used for controlling the early branch (leading branch) and the late branch (lagging branch) in the step 2 are respectively
Figure RE-FDA0002639611080000016
And
Figure RE-FDA0002639611080000017
i.e. the early branch signal sel(t) signal s of the late branchlt(t) is
Figure RE-FDA0002639611080000021
Figure RE-FDA0002639611080000022
3. A PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection as claimed in claim 1 or 2, wherein: the frequency interval value fIs optimally selected by an S curve of frequency error estimationfRb is the symbol rate, 0.16 Rb.
The method of claim 1 for multi-symbol detection based early-late loop frequency synchronization of PCM/FM signals, wherein: step 3, the incoherent multi-symbol detection is to respectively compare the received signals in the ith symbol interval of the two branches of the early branch and the late branch with L groups of local waveforms
Figure RE-FDA0002639611080000023
L is conjugate correlated symbol by symbol to obtain a modulus of the matched filter output
Figure RE-FDA0002639611080000024
Figure RE-FDA0002639611080000025
Figure RE-FDA0002639611080000026
Wherein "s (t)" is the conjugate of the complex number s (t),
Figure RE-FDA0002639611080000027
is a matched filter.
4. The method of claim 4, wherein the PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection is characterized in that: the matched filter
Figure RE-FDA0002639611080000028
Is NxT, T is the symbol period of the PCM/FM signal, N is the number of symbols required for one detection, N is>1；
Then, the ith local waveform is:
Figure RE-FDA0002639611080000029
Figure RE-FDA00026396110800000210
wherein the content of the first and second substances,
Figure RE-FDA00026396110800000211
for N symbols to be selected
Figure RE-FDA00026396110800000212
(i is 1,2, …, N and
Figure RE-FDA00026396110800000213
) The constituent symbol vectors, i.e.
Figure RE-FDA00026396110800000214
Kf is the modulation index of PCM/FM and takes 0.7, and g (t) is the pre-adjusting filter of PCM/FM signal
Figure RE-FDA00026396110800000215
The matched filter
Figure RE-FDA00026396110800000216
The number L of (A) is in a relationship of L-2 with NN。
5. The method of claim 5 for multi-symbol detection based early-late loop frequency synchronization of PCM/FM signals, wherein: the number N of the symbols to be selected is equal to 3, so the matched filter
Figure RE-FDA0002639611080000031
The number L of (c) is 8.
CN202010660213.2A 2020-07-10 2020-07-10 PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection Active CN111865850B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010660213.2A CN111865850B (en) 2020-07-10 2020-07-10 PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010660213.2A CN111865850B (en) 2020-07-10 2020-07-10 PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection

Publications (2)

Publication Number Publication Date
CN111865850A true CN111865850A (en) 2020-10-30
CN111865850B CN111865850B (en) 2022-08-02

Family

ID=73152651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010660213.2A Active CN111865850B (en) 2020-07-10 2020-07-10 PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection

Country Status (1)

Country Link
CN (1) CN111865850B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5711009A (en) * 1996-04-01 1998-01-20 The United States Of America As Represented By The Secretary Of The Navy Method and apparatus for minimizing the effects of frequency selective fading on a PCM/FM data signal
US20060291592A1 (en) * 2004-10-15 2006-12-28 Perrins Erik S Multi-symbol noncoherent CPM detector
CN103457680A (en) * 2013-08-20 2013-12-18 重庆邮电大学 Satellite communication timing synchronization error detection method based on full-digital receiving
CN104486288A (en) * 2014-12-01 2015-04-01 北京理工大学 Carrier frequency deviation suppressing method suitable for PCM/FM (pulse-code modulation/frequency modulation) telemetering receiver
CN107094065A (en) * 2017-03-28 2017-08-25 西安电子科技大学 A kind of remote measurement PCM/FM system transmission methods based on MIMO technology
CN108768604A (en) * 2018-05-08 2018-11-06 北京航空航天大学 A kind of low complex degree bit synchronization method for PCM/FM multiple-symbol detections

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5711009A (en) * 1996-04-01 1998-01-20 The United States Of America As Represented By The Secretary Of The Navy Method and apparatus for minimizing the effects of frequency selective fading on a PCM/FM data signal
US20060291592A1 (en) * 2004-10-15 2006-12-28 Perrins Erik S Multi-symbol noncoherent CPM detector
CN103457680A (en) * 2013-08-20 2013-12-18 重庆邮电大学 Satellite communication timing synchronization error detection method based on full-digital receiving
CN104486288A (en) * 2014-12-01 2015-04-01 北京理工大学 Carrier frequency deviation suppressing method suitable for PCM/FM (pulse-code modulation/frequency modulation) telemetering receiver
CN107094065A (en) * 2017-03-28 2017-08-25 西安电子科技大学 A kind of remote measurement PCM/FM system transmission methods based on MIMO technology
CN108768604A (en) * 2018-05-08 2018-11-06 北京航空航天大学 A kind of low complex degree bit synchronization method for PCM/FM multiple-symbol detections

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
付刚等: "再入遥测环境下的PCM/FM信号软解调研究", 《科学技术与工程》 *
王潇等: "高码率LDPC与MSD技术在PCM/FM遥测体制下性能研究", 《遥测遥控》 *
谢顺钦等: "大动态PCM/FM信号的载波频率同步", 《电讯技术》 *
陈大海等: "遥测接收机的一种载波频偏抑制方法", 《电子科技大学学报》 *

Also Published As

Publication number Publication date
CN111865850B (en) 2022-08-02

Similar Documents

Publication Publication Date Title
US11277254B2 (en) Receiver with enhanced clock and data recovery
US7961823B2 (en) System and method for adjusting multiple control loops using common criteria
CN102136850A (en) Method and device for realizing automatic frequency control
WO2018129843A1 (en) Artm cpm demodulation and synchronization method with low implementation complexity
CN108270715A (en) It is suitble to the carrier recovery system and method for high-order 4096-QAM
CN112399551B (en) High-precision synchronization method for short-time burst signals
US11139949B2 (en) Equalizer adaptation based on eye monitor measurements
CN108494467B (en) Physical layer self-adaptive ACM synchronization equipment based on satellite communication
CN105356993B (en) The channel time delay control method of PCM/FM polarizations synthesis
CN107302409B (en) Automatic gain control method based on signal-to-noise ratio estimation of over-sampled signal
RU2321955C2 (en) Device and method for tracking radio frequency and capturing it
CN114697941A (en) Low-power consumption Bluetooth baseband receiving method
US6154510A (en) Symbol timing recovery based on adjusted, phase-selected magnitude values
CN111865850B (en) PCM/FM signal early-late loop frequency synchronization method based on multi-symbol detection
CN110880964B (en) Bit synchronization tracking system based on data conversion tracking loop
CN110493152B (en) Self-adaptive equalization circuit based on frequency spectrum balancing method
CN114513393A (en) Correction method of analog equalizer, control chip, receiver and storage medium
CN112468281A (en) High-precision symbol synchronization system
CN110943752A (en) OQPSK self-adaptive variable rate digital transceiver based on adjacent space link protocol
CN114465691A (en) Low-complexity constant envelope phase modulation signal sampling deviation estimation and compensation method and system
CN112929310A (en) Carrier recovery method based on high-order QAM
CN114050822A (en) Self-adaptive phase tracking method and system
US20030081699A1 (en) Phase detector
US20050129148A1 (en) Control algorithm in QAM modems
CN113242200B (en) Method for dynamically calculating optimal decision threshold based on 4FSK signal soft demodulation

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