Background
In recent years, as mobile operators, equipment manufacturers, mobile phone manufacturers and government departments of all countries in the world invest considerable funds and energy for the construction and operation of 3 rd generation mobile communication networks (3G), the computer industry has introduced a new broadband wireless access technology and named WiMAX. WiMAX is an abbreviation for Worldwide Interoperability for microwave Access, which is commonly translated into "Worldwide Interoperability for microwave Access".
WiMAX technology is based on the 802.16 series of standards of IEEE. In the 802.16 series standard, the technical requirements of the air interface between the base station bs (base station) and the subscriber station ss (subscriber station) are specified in detail, and particularly, the frame structure requirements of the physical layer, the system design parameters, and the like are specified in detail. The invention is mainly designed for timing synchronization of a receiving end of a WiMAX system (IEEE 802.16-2004) base station based on an OFDM physical layer.
OFDM, an abbreviation of Orthogonal Frequency Division Multiplex, in chinese, means Orthogonal Frequency Division multiplexing. The OFDM technology is based on orthogonal multiple carriers, which is a multi-carrier spread spectrum technology. The great advantage of OFDM is against frequency selective fading or narrow-band interference, which can cause the entire communication link to fail in a single carrier system, but only a small fraction of the carriers are disturbed in a multi-carrier system, and error correction codes can be used for these sub-channels for error correction. In the OFDM system, the carriers of each sub-channel are orthogonal to each other, and the frequency spectrums are overlapped with each other, so that the mutual interference among the sub-carriers is reduced, and the frequency spectrum utilization rate is improved.
In the design of WiMAX system based on OFDM, the requirement for synchronization is quite high in order to ensure orthogonality between subcarriers. Once out of synchronization, orthogonality between subcarriers is affected, thereby severely affecting system performance. Therefore, the superiority and inferiority of the synchronization algorithm can lead to the performance of the whole system. A high-performance and practical synchronization algorithm is designed, and becomes a key link of the whole system design.
The synchronization algorithm includes time domain synchronization and frequency domain synchronization. Frame synchronization is a precondition and a basis for the whole WiMAX system to correctly receive data. The physical layer of the 802.16-2004 protocol based on the OFDM mode supports frame-based transmission. One frame includes one downlink subframe and one uplink subframe. The downlink subframe consists of one downlink physical layer PDU and the uplink subframe consists of contention slots for initial search and bandwidth request purposes and one or more uplink physical layer bursts transmitted by different SSs. The Burst of each uplink SS contains a short preamble.
The downlink physical layer PDU is started by a long preamble for physical layer synchronization. The uplink physical layer PDU (protocol data unit) is started by a short preamble for synchronization of uplink bursts.
An OFDM symbol is a useful symbol time T from the time domainbAnd length T of cyclic prefix CPgSum Ts. The CP is a copy of the length of the end of the useful symbol time to collect multipath and maintain the orthogonality of the subcarriers. The length of the CP should be greater than the maximum multipath delay length. The OFDM symbol is shown in fig. 1.
According to the stipulation of the protocol, the value of the uplink short lead code in the frequency domain sequence of the WiMAX system based on OFDM is determined by the following formula:
<math> <mrow> <msub> <mi>P</mi> <mn>128</mn> </msub> <mo>=</mo> <mtable> <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msqrt> <mn>2</mn> </msqrt> <msub> <mi>P</mi> <mi>ALL</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>·</mo> <mo>·</mo> <mo>·</mo> <msub> <mi>k</mi> <mrow> <mi>mod</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mn>0</mn> <mo>;</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> <mo>·</mo> <mo>·</mo> <mo>·</mo> <msub> <mi>k</mi> <mrow> <mi>mod</mi> <mn>2</mn> </mrow> </msub> <mo>≠</mo> <mn>0</mn> <mo>;</mo> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtable> </mrow> </math>
whereinThe factor is related to the 3dB gain. PALLThe frequency domain sequences are uniformly specified by the protocol and have fixed values. Thus, the preamble received by the base station to the SS (subscriber station) appears in the time domain as two repetitions of 128 samples, preceded by the cyclic prefix of the preamble. As shown in fig. 2.
According to the design of practical system, the base station receives bursts of physical layer of the up-going SS (subscriber station), it must first detect the short preamble in front of the Burst of each SS for frame timing synchronization, and only after detecting the preamble, the base station can receive and process the data Burst.
In an actual WiMAX system, a data processing flow chart at the receiving end of a base station is shown in fig. 3. The upper part of the figure is the data flow of the transmitter and the lower part is the data flow of the receiver. The transmitter performs serial/parallel conversion on the data after modulation processing (constellation mapping), performs IFFT conversion, performs parallel/serial conversion on the converted data, inserts cyclic prefix, windowing and digital-to-analog conversion (DAC), and transmits the data to a wireless communication environment through radio frequency. The receiver is the inverse process of the transmitter and mainly comprises the steps of radio frequency receiving, analog-to-digital conversion (ADC), cyclic prefix removing, serial/parallel conversion, FFT demodulation and the like.
According to the basic principle of OFDM, when multicarrier transmission is performed by using OFDM, the requirement of the base station receiving end on synchronization is quite high in order to correctly receive the transmitted data symbols and improve the performance of the system. The synchronization process at the receiving end of the WiMAX base station is shown in fig. 4.
Detailed Description
The invention provides a joint synchronization method for frame, symbol timing and sampling timing of a receiving end of a base station of a WiMAX system, which has the following basic principle:
let the expression of the SS transmit signal be: x (t) ej(2πfct+*)Wherein x (t) is the waveform of SS transmitting terminal after D/A conversion, fcTo transmit the carrier frequency, * is the carrier phase.
After passing through the transmission channel, the signals at each point at the receiving end of the base station can be represented as:
wherein,frequency and phase estimates are estimated for the carrier. Δ fcΔ * is the carrier frequency offset and phase offset, and η (t) is the channel noise.
Obtaining a data stream of r after samplingn},
Wherein,for the purpose of the sampling interval estimation,
fourier FFT is carried out, and r isnSubstituting the expression of (1) into (2) for simplification to obtain:
<math> <mrow> <msub> <mi>R</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>r</mi> <mi>n</mi> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>πnk</mi> <mo>/</mo> <mi>N</mi> </mrow> </msup> </mrow> </math>
(1-3)
wherein,
are all due to interference caused by noise, S
mAnd m, N and k are 0, 1, 2 … … and N-1 for the frequency domain sequence before IFFT transformation at the SS sending end.
If the timing deviation deltat of the upper sample is considered againy(corresponding to the reception signal y (t) in <math> <mrow> <mi>t</mi> <mo>=</mo> <mi>n</mi> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>-</mo> <mi>Δ</mi> <msub> <mi>t</mi> <mi>y</mi> </msub> </mrow> </math> Sampling at a moment to obtain rn(ii) a Reconsidering the presence of symbol timing deviation Δ tfEquivalent to following from the cyclic prefix <math> <mrow> <mi>Δ</mi> <msub> <mi>n</mi> <mi>f</mi> </msub> <mo>=</mo> <mi>Δ</mi> <msub> <mi>t</mi> <mi>f</mi> </msub> <mo>/</mo> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> </mrow> </math> The FFT operation is performed for the first N samples of the sample. Here we assume Δ nf,ΔtfTaking a value of positive means leading the optimum symbol sample timing; a negative value indicates a lagging optimum symbol sample timing.
An OFDM symbol that takes into account sample timing offset and symbol timing offset is shown in fig. 5.
The output signal after FFT demodulation at the receiving end while comprehensively considering the carrier synchronization deviation, the sample synchronization deviation, and the symbol synchronization deviation is:
<math> <mrow> <msub> <mi>R</mi> <mi>k</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>S</mi> <mi>m</mi> </msub> <mfrac> <mrow> <mi>sin</mi> <mo>[</mo> <mi>π</mi> <mrow> <mo>(</mo> <mi>m</mi> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <mi>k</mi> <mo>-</mo> <mi>NΔ</mi> <msub> <mi>f</mi> <mi>c</mi> </msub> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mrow> <mi>sin</mi> <mo>[</mo> <mi>π</mi> <mrow> <mo>(</mo> <mi>m</mi> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <mi>k</mi> <mo>-</mo> <mi>NΔ</mi> <msub> <mi>f</mi> <mi>c</mi> </msub> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>N</mi> <mo>]</mo> </mrow> </mfrac> <mo>×</mo> </mrow> </math>
(1-4)
as can be seen from the above formula, only Δ fc,ΔTsNot only causes the phase of the signal to deflect, but also changes the amplitude of the signal, bringing about the reduction of the signal-to-noise ratio (SNR). Other deviations Δ *, Δ ty,ΔnfAnd noise etc. only cause the signal phase to deflect without affecting the signal amplitude.
For the convenience of analysis, the product term between various deviations is considered to be small, and after the product term is ignored, the following formula can be simplified:
<math> <mrow> <msub> <mi>R</mi> <mi>k</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <mi></mi> <msub> <mi>S</mi> <mi>k</mi> </msub> <mfrac> <mrow> <mi>sin</mi> <mo>[</mo> <mi>π</mi> <mrow> <mo>(</mo> <mi>kΔ</mi> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <mi>NΔ</mi> <msub> <mi>f</mi> <mi>c</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mrow> <mi>sin</mi> <mo>[</mo> <mi>π</mi> <mrow> <mo>(</mo> <mi>kΔ</mi> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <mi>NΔ</mi> <msub> <mi>f</mi> <mi>c</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>N</mi> <mo>]</mo> </mrow> </mfrac> <mo>×</mo> </mrow> </math>
<math> <mrow> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>Σ</mi> <munder> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>m</mi> <mo>≠</mo> <mi>k</mi> </mrow> </munder> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>S</mi> <mi>m</mi> </msub> <mfrac> <mrow> <mi>sin</mi> <mo>[</mo> <mi>π</mi> <mrow> <mo>(</mo> <mi>m</mi> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <mi>k</mi> <mo>-</mo> <mi>NΔ</mi> <msub> <mi>f</mi> <mi>c</mi> </msub> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mrow> <mi>sin</mi> <mo>[</mo> <mi>π</mi> <mrow> <mo>(</mo> <mi>m</mi> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <mi>k</mi> <mo>-</mo> <mi>NΔ</mi> <msub> <mi>f</mi> <mi>c</mi> </msub> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>N</mi> <mo>]</mo> </mrow> </mfrac> <mo>×</mo> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mi>S</mi> <mi>k</mi> </msub> <msub> <mi>I</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>+</mo> <munderover> <mi>Σ</mi> <munder> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>m</mi> <mo>≠</mo> <mi>k</mi> </mrow> </munder> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>S</mi> <mi>m</mi> </msub> <msub> <mi>I</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> <mo>+</mo> <msup> <mi>η</mi> <mrow> <mo>′</mo> <mo>′</mo> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein
Indicating interference of other subcarriers with the desired carrier.
First term on right side of equation (1-5)I.e. the desired data symbol term, the second term is the effect on the desired symbol caused by the symbols transmitted on the other subcarriers. It is due to carrier frequency deviation, carrier phase deviation, sample frequency deviation, sample timing deviation, symbol timing deviation, noise, etc. that the desired data symbols are passed through the coefficients Ik,kWeighting and interference on other sub-carriers.
In the practical WiMAX system, after the initial sample and symbol coarse synchronization, the frequency offset fine synchronization and the channel estimation and correction, Δ f is not considered for the moment for the convenience of analysisc、ΔTsThe simplified equation (1-5) can now be:
it follows that, assuming that the phase deviation of the received signal with respect to the transmitted signal is α, then
α=-2πkΔty/(NTs)-2πkΔnf/N-Δ*-α′ (1-7)
Where α' is the phase deflection due to noise, with some randomness. The simplification (1-7) can obtain:
α=-2πk[Δty/(NTs)+Δnf/N]-Δ*-α′ (1-8)
through analysis, the phase deviation caused by the first term on the right side of the expression (1-8) is in direct proportion to the subcarrier number k, the phase deviation of the second term on the right side is the same for all subcarriers, namely, the second term is a quantity which is independent of the carrier number k, and the third term on the right side is the phase deviation caused by other factors such as noise and the like, and the phase deviation has certain randomness.
Aiming at the analysis, the invention provides a combined algorithm of fine timing synchronization and sample value deviation timing synchronization at a receiving end of a base station of a WiMAX system, which comprises the following steps:
and (one) performing time domain frame coarse synchronization before FFT transformation according to a short preamble of a known sequence transmitted by a SS transmitting end, and determining a starting point nO of the OFDM symbol (the coarse synchronization can be determined by performing delay correlation through 128 repeated samples of the short preamble and by an output peak value of a correlator).
(II) the frequency domain sequence value R of the short preamble receiving signal sequence after the FFT transformation after the frequency offset correction and the channel estimationk(k is 0, 1, … … 255) and the known short preamble frequency domain sequence value are subjected to conjugate multiplication to obtain an intermediate variable Rk′=RkSk *In which S isk *Represents a pair SkTaking conjugation, taking Rk' the argument yields the sequence value αk=arg(Rk'), (k ═ 0, 1, … … 255), where arg () denotes the argument.
(III) taking betak=αk+1-αk(k-0, 1, … … 254), converting the sequence β into a sequence βkTaking an average to obtain <math> <mrow> <mover> <mi>β</mi> <mo>‾</mo> </mover> <mi></mi> <mo>=</mo> <mi></mi> <mfrac> <mn>1</mn> <mn>255</mn> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>254</mn> </munderover> <msub> <mi>β</mi> <mi>k</mi> </msub> </mrow> </math> Taking the intermediate variable theta-256 beta/2 pi, the integer part of theta is delta nfThe fractional part is Δ ty/TsAccording to Δ ny、ΔtyAs a result, the OFDM time domain frame timing and sample timing before FFT conversion are controlled, thereby achieving precise timing synchronization and sample timing synchronization for OFDM symbols.
The flow chart of the present invention is shown in fig. 6 (the flow of frequency offset synchronization, channel estimation, etc. is omitted in the figure).
The algorithm is further described below with reference to the algorithm flow chart of fig. 6:
firstly, the WiMAX base station performs analog/digital conversion on the received signal sent by the uplink SS, performs frame timing coarse synchronization (the synchronization method may adopt delay correlation to find the maximum value) through the time domain structure feature of the uplink preamble (the uplink preamble is an OFDM symbol, the front is a cyclic prefix, and then two 128 sample sequences are repeated), and determines the coarse timing synchronization point n0 as the start position of the preamble.
According to the result of frame coarse synchronization timing, we can remove the cyclic prefix sample according to the CP length specified by the system, and perform FFT fourier transform on the OFDM valid data to obtain the frequency domain sequence value of the preamble symbol (note that before and after FFT, it is necessary to perform frequency offset correction and channel estimation according to the result of coarse timing synchronization to complete necessary signal processing and correctly receive signals).
Since the protocol specifies the frequency domain sequence value of the SS short preamble symbol of the WiMAX system, the received short preamble frequency domain sequence value is multiplied by the conjugate of the known frequency domain sequence complex value, Rk′=RkSk *From the previous derivation, we can derive the intermediate variable Rk′,RkThe magnitude of' will be the square of the modulus of the known short preamble frequency domain complex sequence, and its phase will be the sum of the phase deviations due to the various deviations. Including a timing symbol deviation deltan proportional to the subcarrier number kfAnd the sampling timing deviation Δ tyThe phase deviation caused by the carrier frequency is not related to the subcarrier frequency offset number k, and the carrier frequency is not related to the subcarrier frequency offset number k.
By applying to an intermediate variable RkTaking the amplitude angle, we obtain the total deviation sequence alpha of the received short preamble frequency domain sequence caused by various deviationsk(k is 0, 1, … 255), and then alpha is addedkSubtracting the sequence value from the sequence value to obtain another intermediate variable sequence betak=αk+1-αk(k-0, 1, … 254). By the foregoing derivation, weIt can be seen that through the process of this subtraction, various phase deviations caused by various synchronization deviations, which are independent of the subcarrier number k, including the phase deviation caused by the carrier phase deviation Δ *, the residual carrier frequency deviation, the sample frequency deviation, and so on, can be effectively eliminated.
Then beta is mixedkAveraging over all sub-carriers of the whole preamble we get <math> <mrow> <mover> <mi>β</mi> <mo>‾</mo> </mover> <mi></mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>255</mn> </mfrac> <mi></mi> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>254</mn> </munderover> <msub> <mi>β</mi> <mi>k</mi> </msub> <mo>,</mo> </mrow> </math> By this averaging process, random fluctuations of the phase deviation α' caused by various random deviations, mainly random phase deviations caused by various noises, can be effectively eliminated. It is known that one advantage of multi-carrier systems is that they are efficient against frequency selective fading, since only a part of the sub-carriers are affected by noise. This averaging of the phase deviations is also very resistant to random phase shifts caused by random noise.
After processing by the foregoing processes, we obtain a phase offset value that is substantially related only to timing symbol offset and sample timing offset <math> <mrow> <mover> <mi>β</mi> <mo>‾</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mn>255</mn> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>254</mn> </munderover> <msub> <mi>β</mi> <mi>k</mi> </msub> <mo>=</mo> <mo>-</mo> <mn>2</mn> <mi>π</mi> <mo>[</mo> <mi>Δ</mi> <msub> <mi>t</mi> <mi>y</mi> </msub> <mo>/</mo> <mrow> <mo>(</mo> <mi>N</mi> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>Δ</mi> <msub> <mi>n</mi> <mi>f</mi> </msub> <mo>/</mo> <mi>N</mi> <mo>]</mo> <mo>,</mo> </mrow> </math> Where Δ tyIs less than TsError term of, Δ nfIs an error term with integer value, in the WiMAX system, N is 256, TsA constant, representing the sampling time interval, is taken according to the system design, as the inverse of the sampling frequency. Then we multiply β by-256/2 π to get Δ ty/Ts+ΔnfWherein the integer part is Δ nfThe fractional part is Δ ty/TsBy multiplying the fractional part by TsCan obtain Δ tyThe value of (a).
From the previous analysis, we have derived from Δ nfReturning to the step of removing the control cyclic prefix before FFT conversion to obtain accurate timing synchronization according to delta tyThe value of (2) can control the synchronization of the sample timing. Of course, according to the actual system design requirement, in order to reduce the influence of the delay of the system on the performance of the system, the delay can also be based on Δ nfAnd Δ tyTo control cyclic prefix removal and sample timing synchronization of the next OFDM symbol. In particular, if Δ nfTaking a positive value to indicate frame timing advance optimum timing Δ nfCoarse synchronization timing n0 should be delayed by an amount of Δ nfTaking the sample value as a frame synchronization timing point; if Δ nfA value of negative indicates an optimum timing deltan of the frame timing lagfSample number, coarse synchronization timing n0 should be extended forward by anfThe samples serve as frame sync timing points. For Δ tyWith the same operation if Δ tyTaking the value as positive, delaying the sampling moment backward by delta ty(ii) a If Δ tyTaking the value as negative, extending the sampling moment by delta ty。
The traditional timing fine synchronization algorithm is based on pilot frequency carrier waves to carry out accurate timing synchronization, and compared with the traditional fine timing synchronization algorithm, the timing fine synchronization algorithm has the advantages that: based on the short preamble symbol sent by the uplink SS, accurate frame and symbol fine timing synchronization and sample value timing synchronization are carried out on the basis of frame and symbol coarse synchronization, the synchronization timing process of the system is accelerated, and the synchronization timing precision is improved. Meanwhile, according to the design principle of the algorithm, the algorithm has strong resistance to carrier frequency phase fixed deviation and phase random deviation caused by noise, and the anti-interference capability of timing synchronization is improved.