CN100438399C - Timing tracking apparatus, receiver, timing tracking and regulating method - Google Patents
Timing tracking apparatus, receiver, timing tracking and regulating method Download PDFInfo
- Publication number
- CN100438399C CN100438399C CNB2005100735390A CN200510073539A CN100438399C CN 100438399 C CN100438399 C CN 100438399C CN B2005100735390 A CNB2005100735390 A CN B2005100735390A CN 200510073539 A CN200510073539 A CN 200510073539A CN 100438399 C CN100438399 C CN 100438399C
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- data
- mfrac
- msup
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000001105 regulatory effect Effects 0.000 title abstract 2
- 238000005070 sampling Methods 0.000 claims abstract description 61
- 238000012549 training Methods 0.000 claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000001360 synchronised effect Effects 0.000 claims description 16
- 238000012935 Averaging Methods 0.000 claims description 11
- 238000005259 measurement Methods 0.000 claims description 2
- 238000000926 separation method Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
Abstract
The present invention provides a timing tracking device, a receiving device and a timing tracking and regulating method. The timing tracking device uses the data of an early sampling point and a late sampling point to calculate a timing error. The timing tracking device comprises a data separating device, a G (tau/Tc) calculator, an equalizer and a timing error calculator, wherein the data separator is used for separating the data of the early sampling point and the late sampling point from receiving data, the early sampling point and the late sampling point are used for a training sequence, the calculator is used for calculating G (tau/Tc) according to the training sequence data separated from the data separating device, the equalizer is used for calculating the equalizing value of G (tau/Tc) of n sub-frames, and the timing error calculator is used for calculating a timing error according to the equalizing value output by the equalizer.
Description
Technical Field
The present invention relates to a timing tracking device, a receiving device, a timing tracking and adjusting method in a wireless communication system, and more particularly, to a timing tracking device, a receiving device, a timing tracking and adjusting method in a time division synchronous code division multiple access (TD-SCDMA) mobile terminal.
Background
In order to correctly receive information transmitted from a Base Station (BS) in a digital mobile communication system, a mobile terminal (UE) must periodically sample and perform data frame analysis on a received signal with correct clock information, that is, a locally generated PN code needs to be synchronized with a PN code in the received signal, that is, downlink synchronization with the base station needs to be achieved. The downlink synchronization of the mobile terminal is divided into two stages of synchronous establishment and synchronous tracking.
(1) Synchronization establishment phase
The mobile terminal does not know whether the base station transmits a signal at the beginning, and therefore, a search process is required, namely, a synchronous training sequence transmitted by the base station is searched and captured within a certain frequency and time range, and the position of a data frame and the phase of a sampling clock are determined according to the synchronous training sequence. This phase is also referred to as synchronization establishment. That is, the phase difference between the signal transmitted from the other party and the local signal is within the synchronization maintaining range, i.e., within one chip.
(2) Synchronous tracking phase
Once this stage is completed, the synchronous tracking process is entered, i.e. the synchronization is continuously maintained, and the synchronization is not lost due to external influence. That is, no matter what kind of factors the phases at both ends shift, as long as the difference between the shifted phases is within the synchronization holding range, the synchronization shift can be estimated by the synchronization tracking, and the phase of the terminal sampling clock can be adjusted in real time, so that the PN code at the receiving end tracks the change of the PN code at the transmitting end, i.e. the difference between the phases is kept less than a fraction of a chip.
The commonly used synchronous tracking methods are both delay locked loops and tau-jitter loops based on early-late gate timing error detectors. They all belong to the advance-retard type of phase-locked loops. The phase-locked loop functions by correlating the received signal with two locally generated phase-shifted (early and late) signals. The delay-locked loop uses two independent correlators, while the tau-dither loop uses a single correlator for time division. The disadvantage of this conventional synchronous tracking method is that the speed and phase jitter are a contradictory quantity. If the adjustment speed is required to be high, large phase jitter exists; if less phase jitter is required, the corresponding speed of adjustment will also be slower. Therefore, it is necessary to provide a synchronous tracking device having a high speed, a small jitter, and the like.
Disclosure of Invention
In view of the above-described drawbacks, it is an object of the present invention to provide a timing tracking apparatus, a receiving apparatus, a timing tracking and adjusting method, which are fast and have small phase jitter.
An aspect of the present invention is to provide a timing tracking apparatus for calculating a timing error using data of an early sampling point and a late sampling point, comprising a data separator for separating data of the early sampling point and the late sampling point of a training sequence from received data, respectivelyA calculator for calculating the time of the measurement,
<math>
<mrow>
<mover>
<mi>G</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</math>
Pefor early sampling of the power of the data, PlFor the power of the late sample data, for estimating the power of the early sample point data and the power of the late sample point data, respectively, based on the training sequence data of the early sample point data and the late sample point data separated by the data separator, and calculatingA mean value for calculating n sub-framesAnd a timing error calculator for calculating a timing error based on the mean value output by the averager.
Another aspect of the present invention is to provide a receiving apparatus having the timing tracking apparatus as described above, which adjusts timing using a timing error calculated by the timing tracking apparatus.
Another aspect of the present invention is to provide a timing tracking method, which includes a training sequence data separating step of separating early sample point data and late sample point data of a training sequence from received data, respectively, using a data separator,a calculation step of usingA calculator for estimating power of the early sample point data and power of the late sample point data respectively according to the training sequence data of the early sample point data and the late sample point data separated by the data separator, and calculating
<math>
<mrow>
<mover>
<mi>G</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</math>
A value of (b), wherein PeFor early sampling of the power of the data, PlIs the power of the late sampled data; averaging step, using an averager to calculate n sub-framesAnd a timing error calculation step of calculating a timing error based on a result of the averaging step using a timing error calculator.
Another aspect of the present invention is to provide a timing adjustment method including a step of adjusting a timing using a timing error calculated by the timing tracking method as described above.
According to the invention, because the early sampling point and the late sampling point have different power/amplitude when the timing error is not 0, once the timing tracking loop converges within the allowable range of timing precision, that is, the timing error is smaller than the required precision, the power/amplitude difference between the early sampling point and the late sampling point is smaller than a fixed threshold value. Therefore, the method has the advantages of high speed, small phase jitter and the like, and can effectively keep the downlink synchronization within the requirement of a precision range, thereby effectively avoiding the system communication quality reduction and even communication interruption caused by the synchronization adjustment and tracking failure of the receiver.
Drawings
Fig. 1 is a power-time delay diagram for a multipath channel.
Fig. 2 is a block diagram of a receiving apparatus according to the present invention.
Fig. 3 is a block diagram of a timing tracking device.
Fig. 4 is a structure diagram of TD-SCDMA subframe.
Fig. 5 is a structure diagram of TD-SCDMA time slot.
FIG. 6 is a structural diagram of DwPTS.
FIG. 8 is a diagram of a correlation-based implementation methodStructural block diagram of the calculator.
Detailed Description
As shown in fig. 2, the receiving apparatus according to the present invention includes an rf front end 21, an automatic gain controller 22, an analog-to-digital converter 23, a cosine transform filter 24, a receiver 26 and a timing tracking apparatus 25. The structure and function of the rf front-end 21, the agc 22, the adc 23, the radiconic cosine filter 24 and the receiver 26 are the same as in the prior art and their construction will not be described in detail here. The timing tracking means 25 receives the two paths of data of the late and early sampling points from the root cosine filter 24, and estimates a timing error based on the received data. Here, the late and early sampling points are for the ideal sampling point. The timing error obtained by the timing tracking means 25 is inputted to the analog-to-digital converter 23 to adjust the timing. The estimation of timing error needs two paths of over-sampled data, and the receiver receives normal data demodulation and needs one path of sampled data, so the receiving device in the invention at least needs three times of over-sampled data.
The principle of timing tracking in the present invention is that when the timing error is not 0, the early and late sampling points have different powers/amplitudes, so that the different powers/amplitudes of the early and late sampling points can be used to estimate the timing error. Once the timing tracking loop converges within the timing accuracy tolerance, i.e. the timing error is less than the required accuracy, the power/amplitude difference between the early and late sampling points will be less than a fixed threshold. The estimated timing error is for the sample that is intermediate the late sample and the early sample.
In an actual communication system, signals reach a receiving end from a transmitting end through different reflection paths to form multipath signals. Different propagation paths of a multipath signal cause different propagation delays, i.e., multipath effects. Fig. 1 shows a power-delay diagram of a multipath channel, wherein the multipath signal includes three significantly different paths, hereinafter referred to as effective path components. The timing tracking method can track any effective path, but preferably tracks the first effective path, so that the receiver can apply the effective path information as much as possible during data demodulation, and the data demodulation performance of the receiver is improved. In the following, only the tracking of the first effective path will be described.
Fig. 3 shows a block diagram of the timing tracking device 25 according to the present invention. As shown in FIG. 3, the timing estimation device 25 includes a data separator 251, a A calculator 252, an averager 253, and a timing error calculator 254.
The data separator 251 separates the data of the early sampling point and the data of the late sampling point of the training sequence from the received data of the whole sub-frame, and obtains training sequence data of the early sampling point and training sequence data of the late sampling point. Since in most digital mobile communication systems, channel estimation can be performed by a training sequence, in TD-SCDMA systems (subframe structure is shown in fig. 4), the training sequence for downlink data can include midamble or SYNC-DL data of DwPTS portion. Therefore, the midamble or the early and late sampling data of SYNC-DL can be directly applied to perform timing error estimation in the TD-SCDMA system. Fig. 5 shows the time slot structure of TD-SCDMA, midamble is 144 chips in the TD-SCDMA time slot. SYNC-DL is located in DwPTS shown in fig. 6.
Therefore, the training sequence data obtained by the data separator 251 may be midamble data or SYNC-DL data. The preferred embodiment of the present invention is that the data separator 251 separates the last 128 chips of midamble of 144 chips in the time slot. In the following, an implementation of separating midamble by the data separator 251 will be described as an example, which is also applicable to the implementation of SYNC-DL.
The calculator 252 estimates powers of the early and late sample data based on the early and late sample data of the midamble part input from the data separator 251, and calculates the power based on the estimated powersThe early and late sampling data of the midamble part are respectively marked as emid_e(n), n-1, 2, …, 128 and emid_l(n),n=1,2,…,128。Can be formulated as:
here, TcAnd τ are the chip period and timing error, respectively. Delta is the same time difference between the ideal sample point and the late sample point, respectively, the early sample point.Is the result of the estimation of R (tau). For the root-raised cosine filter, when a is 0.22,
there are two ways to implement calculator 252. Implementation method 1 is based on FFT implementation method to calculateThe method is called as an FFT-Based timing tracking (FFT-Based timing tracking) method, and a structural block diagram Based on the method is shown in FIG. 7. The implementation method 2 is based on the calculation of the related implementation methodThis is called Correlation-Based Timing Tracking (Correlation-Based Tracking) method in the present invention, and FIG. 8 shows a structural block Based on this methodFigure (a). These two implementations will be described separately in detail below.
As shown in fig. 7, the FFT-based implementation is calculatedIs/are as followsThe calculator 252 includes a channel estimator 71 for calculating the power of the early-sampled data, a first effective path decimator 72 and a power calculator 73, a channel estimator 71 for calculating the power of the late-sampled data, a first effective path decimator 72 and a power calculator 73, an adder 75, a subtractor 74, and a divider 76.
When the early-sampled data is input to the channel estimator 71, the channel estimator 71 calculates a channel impulse response from the received midamble and the local midamble. In the channel estimator 71, the channel impulse response heCan be calculated by the following formula,
he=IFFT[FFT(emid_e)./FFT(mid)] (3)
here, he=[he(1),he(2),...,he(128)]Is the estimated channel impulse response, emid_e=[emid_e(1),emid_e(2),...,emid_e(128)]Is the early sampled data portion, L, of the received midamble datamid_e=[lmid(1),lmid(2),...,lmid(128)]Is the local midamble data. FFT (), IFFT () respectively represent fast fourier transform and inverse transform of the data sequence in brackets.
The channel impulse response calculated by the channel estimator 71 is input to a first effective path decimator 72. A first effective path extractor 72 extracts a first effective path in the channel profile window (the first value in the channel profile window is defaulted to the first effective path) based on the known user code information (i.e., any code in the user slot). In the present invention, the default value of the code channel number used is the smallest code channel of the user in the user time slot. The extracted position can be calculated by:
Nextact=112-(K-kuser)*W (4)
here, K is the number of the largest midamble shift values in one cell. W is the shift value of two adjacent miambles. k is a radical ofuserIndicating the sequence number of the smallest code channel in the user time slot. When the number of the largest midamble shift values is K,(represents the largest integer less than ·). If we consider using the beacon code track of TSO to implement FFT-based timing tracking, Nextact0 is a fixed value.
The power calculator 73 calculates the power value of the extracted first effective diameter:
Pe=|he(Nextact)|2 (5)
similarly, the power of the late sampled data is also calculated by a channel estimator 71, a first effective path decimator 72 and a power calculator 73, in the same configuration as that used for the calculation of the power of the early sampled data. The difference is that the received data used in the calculation is late sampled midamble data. The late sample data power is noted as:
Pl=|hl(Nextact)|2 (6)
the adder 75, the subtractor 75 and the divider 76 calculate based on the power of the early sampled data and the power of the late sampled dataIs formulated as:
as shown in fig. 8, the calculation is based on the implementation of correlationOfThe calculator 252' includes a correlator 81 and a power calculator 82 for calculating the power of early-sampled data, a correlator 81 and a power calculator 82 for calculating the power of late-sampled data, a shift enable generator 86, an adder 84, a subtractor 83, and a divider 85.
The shifted midamble generator 86 generates shifted midamble codes from the basic midamble codes and the code channel information of the user. The generation process comprises the generation of a complex midamble and the generation of a shifted midamble.
The generation of complex midambles is a process of converting a basic midamble into a complex midamble code, which is described as:
m(n)=jn*m(n);n=1,2,...,128 (8)
here, m (n) is a binary representation of the basic midamble.
here, n is 1, 2, …, 128. K is the number of the largest midamble shift values in a cell. W is the shift value of two adjacent midambles. k is a radical ofuserIndicating the sequence number of the smallest code channel in the user time slot. When the number of the largest midamble shift values is K,(represents the largest integer less than ·). If we consider using the beacon code track of TSO to realize the timing tracking based on correlation, (K-K)user) W-112 is a fixed value.
The shifted midamble code generated by the shifted midamble generator 86 is input to the correlator 81. The correlator 81 calculates the correlation result of the shifted midamble and the received early sampled midamble. Correlation result R of early sampled dataeCan be expressed as:
the correlation result ReInto the power calculator 82. The power calculator 82 calculates the power value of the extracted first effective diameter from the input correlation result:
Pe=|Re|2 (11)
the late sampled data power is also calculated by a correlator 81 and a power calculator 82, which are the same in structure as those employed for the calculation of the early sampled data power. The difference is that the received data used in the calculation is late sampled midamble data.
Correlation result R of late sampled datalCan be expressed as:
the power of the late sampled data is noted as:
Pl=|Rl|2 (13)
the adder 84, the subtractor 83 and the divider 85 calculate based on the power of the early sampled data and the power of the late sampled dataIs formulated as:
calculated by calculator 252And input to the averager 253. The averager 253 calculates n sub-framesThe result can be expressed as:
where n is a numberNumber of subframes of the mean. The minimum value of n may be 1. When the timing variation is within the allowable accuracy requirement (e.g., 1/8 chips), the larger the value of n,the closer the mean is to the statistical mean.
The timing error calculator 254 calculates a timing error from the result of the averager 253. There are two ways to implement this calculation:
a first way to implement this is for the tracking loop to make a fixed step adjustment of the timing. Here, the adjusted timing error value is calculated based on the threshold value T and the adjustment step value α. Timing error adjustment valueCan be calculated using the following equation:
A second way to implement this is for the tracking loop to make a non-fixed step adjustment of the timing.
Here, the timing error adjustment value is based on the threshold value T, the adjustment step value alpha and the oversampling multiple Sn(i.e., the sampling clock frequency) are calculated jointly.
Here, ,maximum integer representing absolute value smaller than · sign [ ·]The symbol of. T is a threshold value, which can be defined by two of the above-mentionedA method for determining; k isThe slope of the function curve, as described above.
When the relation between the oversampling multiple of the system and the required precision is 1/SnThe second implementation is fully equivalent to the first implementation when α.
When the relation between the oversampling multiple of the system and the required precision is 1/SnWhen the value is less than alpha, the convergence effect of the timing error of the second implementation method is better than that of the first implementation method.
For the two methods, the threshold value and the adjustment step value may be determined according to the accuracy requirement of the timing, and taking 1/8 chips as an example, the adjustment step may be set to α -1/8.
There are two calculation methods for the threshold value T.
The first method can be directly based on the precision requirementThe slope of the function curve at 0 point is set, for example, in TD-SCDMA system, the filter is root raised cosine filter, soThe slope of the function curve is K1.402, and if 1/8 chips of accuracy are required, the threshold T is equal to:
the second method can be directly obtained by using a formula:
T=G(Trequired/Tc) (19)
here, TrequiredIs the required progress accuracy of the synchronization tracking,
<math>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>/</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>,</mo>
</mrow>
</math>
taking a root raised cosine filter with a being 0.22 as an example, R (τ) is shown in formula (2), and the calculation method of the present invention is also suitable for other filters.
In the following, a timing tracking method according to the present invention is described. The timing tracking method comprises the following steps:
a training sequence data separating step of separating early sample point data and late sample point data of the training sequence from the received data by using a data separator 251,calculation step of usingA calculator 252 for calculating training sequence data according to the training sequence data separated by the data separatorAn averaging step of calculating n sub-frames using the averager 253And a timing error calculation step of calculating a timing error based on the result of the averaging step using the timing error calculator 254.
In the training sequence data separation step, the separated training sequence data is midamble data or SYNC-DL data.
There are two implementation methods for the calculation step of (a) and (b), namely an FFT-based implementation method and a correlation-based implementation method.
The FFT-based implementation method comprises calculating the power of the early sampled data according to the training sequence data of the early sampled point, the local training sequence and the user code channel informationStep of calculating power of late sampling data according to training sequence data of late sampling point, local training sequence and user code channel information, and calculating power of early sampling data and power of late sampling data according to calculated power of early sampling dataIn the step (2) of (a),
<math>
<mrow>
<mover>
<mi>G</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</math>
Pefor early sampling of the power of the data, PlIs the power of the late sampled data.
The related realization method comprises a step of calculating the power of early sampling data according to the training sequence data of the early sampling point and the shifting training sequence, a step of calculating the power of late sampling data according to the training sequence data of the late sampling point and the shifting training sequence, and a step of calculating the power of the early sampling data and the power of the late sampling data according to the calculated power of the early sampling data and the calculated power of the late sampling dataIn the step (2) of (a),
<math>
<mrow>
<mover>
<mi>G</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</math>
Pefor early sampling of the power of the data, PlIs the power of the late sampled data.
The timing error adjustment value calculation step has two implementation methods. The first method for realizing the timing error adjustment value calculation step is to calculate the timing error adjustment value according to the threshold value T and the adjustment step value alpha
The second method for realizing the timing error adjustment value calculation step is to calculate the adjustment step value alpha and the oversampling multiple S according to the threshold value TnCalculating a timing error adjustment value
Here, ,as a result of the averaging step, the average value,maximum integer representing absolute value smaller than · sign [ ·]The symbol of.
In both implementations of the timing error calculation step described above, the threshold value T is equal toThe slope of the function curve at point 0 is multiplied by the required accuracy value, or by the formula T ═ G (T)required/Tc) And (4) calculating. Here, TrequiredIs the required accuracy of the synchronous tracking progress,
<math>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>/</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>,</mo>
</mrow>
</math>
Tcand τ is the chip period and timing error, respectively, and Δ is the same time difference between the ideal sample and the late and early samples, respectively.
The timing adjustment method according to the present invention adjusts the timing using the timing error calculated by the above-described timing tracking method.
Additional advantages and modifications will readily appear to those skilled in the art. Therefore, the present invention is not limited to the embodiments described in the specification. Accordingly, various modifications may be made without departing from the spirit or scope of the invention as defined by the appended claims and their equivalents.
Claims (22)
1. A timing tracking apparatus for calculating a timing error using data at an early sample point and a late sample point, comprising
A data separator for separating the early sampling point and late sampling point data of the training sequence from the received data,
aA calculator for calculating the time of the measurement,
<math>
<mrow>
<mover>
<mi>G</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</math>
Pefor early sampling of the power of the data, PlFor late sampling of the power of the data, saidThe calculator is used for respectively estimating the power of the early sampling point data and the power of the late sampling point data according to the training sequence data of the early sampling point data and the late sampling point data separated by the data separator, and calculatingThe value of (a) is,
And the timing error calculator is used for calculating the timing error according to the average value output by the averager.
2. The timing tracking device of claim 1,
the training sequence data separated by the data separator is midamble data or SYNC-DL data.
3. The timing tracking device of claim 1,
the calculator being based on FFTA calculator comprising
A channel estimator for calculating power of the early sampled data based on the training sequence data of the early sampled point, the local training sequence and the user code channel information, an effective path extractor and a power calculator,
a channel estimator for calculating power of late sampled data in accordance with training sequence data of late sampled points, local training sequence and user code channel information, an effective path extractor and a power calculator, and a power calculator for calculating power P of early sampled data in accordance with calculated power P of late sampled dataeAnd power P of the late sampled datalComputingThe computing unit of (1).
4. The timing tracking device of claim 1,
A correlator and a power calculator for calculating the power of the early sample data in accordance with the training sequence data of the early sample point, the shifted training sequence,
a correlator and a power calculator for calculating the power of the late sampled data in accordance with the training sequence data of the late samples, the shifted training sequence, and
5. The timing tracking device of claim 1,
the timing error calculator calculates a timing error adjustment value according to the threshold value T and the adjustment step value alpha
6. The timing tracking device of claim 1,
the timing error calculator adjusts the step value alpha and the oversampling multiple S according to the threshold value TnCalculating a timing error adjustment value
9. The timing tracking device of claim 5,
T=G(Trequired/Tc)
here, TrequiredIs the required accuracy of the synchronous tracking progress,
<math>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>/</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>,</mo>
</mrow>
</math>
Teand τ is the chip period and timing error, respectively, and Δ is the same time difference between the ideal sample and the late and early samples, respectively.
10. The timing tracking device of claim 6,
T=G(Trequired/Tc)
here, TrequiredIs the required accuracy of the synchronous tracking progress,
<math>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>/</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>,</mo>
</mrow>
</math>
Tcand τ is the chip period and timing error, respectively, and Δ is the same time difference between the ideal sample and the late and early samples, respectively.
11. A receiving device having a timing tracking device according to any one of claims 1 to 10 which adjusts the timing using the timing error calculated by the timing tracking device.
12. A timing tracking method, comprising
A training sequence data separation step, wherein a data separator is used to separate the early sampling point data and the late sampling point data of the training sequence from the received data respectively,
calculating stepProcedure of usingA calculator for estimating power of the early sample point data and power of the late sample point data respectively according to the training sequence data of the early sample point data and the late sample point data separated by the data separator, and calculating
<math>
<mrow>
<mover>
<mi>G</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</math>
A value of (b), wherein PcFor early sampling of the power of the data, PlIn order to late-sample the power of the data,
averaging step, using an averager to calculate n sub-framesIs a mean value of
And a timing error calculation step of calculating a timing error according to the result of the averaging step by using a timing error calculator.
13. The timing tracking method of claim 12, wherein
The separated training sequence data is midamble data or SYNC-DL data.
14. The timing tracking method of claim 12,
calculating the power of the early sampling data according to the training sequence data of the early sampling point, the local training sequence and the user code channel information,
calculating the power of the late sampled data based on the training sequence data of the late samples, the local training sequence and the user code channel information, an
According to the calculated power P of the early sampling dataeAnd power P of the late sampled datalComputingThe step (2).
15. The timing tracking method of claim 12,
Calculating the power of the early sample data by shifting the training sequence according to the training sequence data of the early sample point,
calculating the power of the late sampled data by shifting the training sequence in accordance with the training sequence data of the late sampled points, an
16. The timing tracking method of claim 12,
the timing error calculation step is to calculate the timing according to the threshold value T and the adjustment step value alphaTime error adjustment value
Wherein,the result obtained in the averaging step.
17. The timing tracking method of claim 12,
the timing error calculation step is based on the threshold value T, the adjustment step value alpha and the oversampling multiple SnCalculating a timing error adjustment value
20. The timing tracking method of claim 16,
T=G(Trequired/Tc)
here, TrequiredIs the required accuracy of the synchronous tracking progress,
<math>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>/</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>,</mo>
</mrow>
</math>
Tcand τ is the chip period and timing error, respectively, and Δ is the same time difference between the ideal sample and the late and early samples, respectively.
21. The timing tracking method of claim 17,
T=G(Trequired/Tc)
here, TrequiredIs the required accuracy of the synchronous tracking progress,
<math>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>τ</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>/</mo>
<mo>[</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>-</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>τ</mi>
<mo>+</mo>
<mi>Δ</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>]</mo>
<mo>,</mo>
</mrow>
</math>
Tcand τ is the chip period and timing error, respectively, and Δ is the same time difference between the ideal sample and the late and early samples, respectively.
22. A timing adjustment method comprises the steps of,
a step of adjusting the timing by using the timing error calculated by the timing tracking method according to any one of claims 12 to 21.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100735390A CN100438399C (en) | 2005-06-02 | 2005-06-02 | Timing tracking apparatus, receiver, timing tracking and regulating method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100735390A CN100438399C (en) | 2005-06-02 | 2005-06-02 | Timing tracking apparatus, receiver, timing tracking and regulating method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1688120A CN1688120A (en) | 2005-10-26 |
CN100438399C true CN100438399C (en) | 2008-11-26 |
Family
ID=35306160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100735390A Expired - Fee Related CN100438399C (en) | 2005-06-02 | 2005-06-02 | Timing tracking apparatus, receiver, timing tracking and regulating method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100438399C (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101534426B (en) * | 2008-03-11 | 2010-12-01 | 卓胜微电子(上海)有限公司 | National standard timing tracking method of digital television terrestrial broadcasting |
CN102196553A (en) * | 2011-03-31 | 2011-09-21 | 重庆重邮信科通信技术有限公司 | Timing synchronization correcting method of mobile terminal |
CN103368873B (en) * | 2012-03-27 | 2016-09-28 | 北京电信技术发展产业协会 | A kind of method obtaining sampling dithering |
CN108809869B (en) * | 2017-04-28 | 2020-12-29 | 深圳市中兴微电子技术有限公司 | A method and device for controlling sampling time |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1267978A (en) * | 1999-11-10 | 2000-09-27 | 深圳市中兴通讯股份有限公司 | Precise PN code synchronizing method and device for wide-band CDMA system |
CN1322073A (en) * | 2000-07-26 | 2001-11-14 | 深圳市中兴通讯股份有限公司 | Code tracking and despreading and descramble method of wideband code division multiple access system and its device |
CN1327312A (en) * | 2000-06-07 | 2001-12-19 | 华为技术有限公司 | Forming method of training sequence in channel evaluation |
WO2003028238A1 (en) * | 2001-09-03 | 2003-04-03 | Infineon Technologies Ag | Device and method for fine synchronization on the sampling of spread-coded received signals |
US20040085918A1 (en) * | 2002-10-21 | 2004-05-06 | Sanyogita Shamsunder | Method and apparatus for block-based chip timing estimation in a code division multiple access communication system |
WO2004047326A1 (en) * | 2002-11-15 | 2004-06-03 | Telecom Italia S.P.A. | Method and device for fine synchronization of a digital telecommunication receiver |
-
2005
- 2005-06-02 CN CNB2005100735390A patent/CN100438399C/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1267978A (en) * | 1999-11-10 | 2000-09-27 | 深圳市中兴通讯股份有限公司 | Precise PN code synchronizing method and device for wide-band CDMA system |
CN1327312A (en) * | 2000-06-07 | 2001-12-19 | 华为技术有限公司 | Forming method of training sequence in channel evaluation |
CN1322073A (en) * | 2000-07-26 | 2001-11-14 | 深圳市中兴通讯股份有限公司 | Code tracking and despreading and descramble method of wideband code division multiple access system and its device |
WO2003028238A1 (en) * | 2001-09-03 | 2003-04-03 | Infineon Technologies Ag | Device and method for fine synchronization on the sampling of spread-coded received signals |
US20040085918A1 (en) * | 2002-10-21 | 2004-05-06 | Sanyogita Shamsunder | Method and apparatus for block-based chip timing estimation in a code division multiple access communication system |
WO2004047326A1 (en) * | 2002-11-15 | 2004-06-03 | Telecom Italia S.P.A. | Method and device for fine synchronization of a digital telecommunication receiver |
Also Published As
Publication number | Publication date |
---|---|
CN1688120A (en) | 2005-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7068593B2 (en) | Apparatus and method for synchronizing frequency in orthogonal frequency division multiplexing communication system | |
JP3125644B2 (en) | Demodulator | |
TWI247549B (en) | Code tracking loop with automatic power normalization | |
CN101189847B (en) | Fine timing acquisition | |
EP1618696B1 (en) | Frequency synchronization apparatus and frequency synchronization method | |
USRE41607E1 (en) | Data aided symbol timing system for precoded continuous phase modulated signals | |
US7298806B1 (en) | Method and system for data-aided timing offset estimation for frequency selective fading channels | |
CN108449298B (en) | Timing synchronization and frequency offset estimation method suitable for high-speed broadband communication | |
JP2004532587A (en) | CDMA searcher with time lag compensation | |
CN109756968B (en) | Precise synchronization timing method and device for single carrier spread spectrum system | |
US8000415B2 (en) | Method and device for detecting a synchronization signal in a communication system | |
CN112953593B (en) | LoRa advanced receiver | |
US8005157B2 (en) | Correlation calculation unit and correlation calculation apparatus | |
TW201310952A (en) | OFDM receivers, compensating devices and methods for detecting and compensating for a sampling clock offset in a receiver | |
CN100361424C (en) | Data Transmission Method in Short Range Wireless Network | |
CN103441969B (en) | The carrier frequency bias estimation of pilot frequency space variable | |
CN100438399C (en) | Timing tracking apparatus, receiver, timing tracking and regulating method | |
CN114465691A (en) | Low-complexity constant envelope phase modulation signal sampling deviation estimation and compensation method and system | |
KR20170096725A (en) | method for processing WCDMA signal timing offset for signal analyzing equipment | |
CN109561042B (en) | A Timing Frequency Synchronization Method of OFDM System Receiver | |
US7643602B2 (en) | Method and system for estimating frequency offsets | |
KR102574425B1 (en) | Apparatus for receiving Orthogonal Frequency Division Multiplexing signal based on baseband and Method for estimating error of sample clock | |
CN101095291B (en) | Fine time tracking with improved resolution | |
CN113162717A (en) | Feedforward timing synchronization method and system for burst communication | |
CN1925470B (en) | Multi-path judging method for time-division duplex mobile communication system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20081126 Termination date: 20180602 |
|
CF01 | Termination of patent right due to non-payment of annual fee |