The present invention relates to the field of mobile communications, especially to a coarse timing synchronization acquisition method in a mobile communication system.
In order to provide high-speed multimedia transmission services, it is necessary in a broadband wireless access system to adopt advanced techniques which counteract the frequency-selective fading effects of wireless channels. Many researchers have appreciated the use of Orthogonal Frequency Division Multiplexing (OFDM), which is well known as an excellent approach for alleviating the detrimental effects of frequency selective fading, and which is a preferred physical transmission technique employed in various wireless standards, and is to be standardized and industrialized.
Mobile multimedia broadcast is an important application of the OFDM technique. An OFDM transmission system requires both frequency synchronization and timing synchronization. Timing synchronization further comprises coarse timing synchronization and fine timing synchronization. The coarse timing synchronization searches in received signals to acquire the timing of the transmission of an OFDM signal, and so as to determine the start time instant of an OFDM symbol. OFDM systems are vulnerable to timing errors as they can introduce inter-subcarrier interference by destroying the orthogonality between subcarriers, and finally lead to the performance degradation at the receiver. Existing coarse timing synchronization methods not only require the receiver to make use of a high-precision crystal oscillator but also impose a heavy computational burden on the timing synchronization module.
The present invention aims to solve the above-mentioned technical shortcomings by providing a coarse timing synchronization acquisition method in a mobile communication system, applied in a mobile communication system that transmits the same signal block repeatedly, especially in a mobile multimedia broadcast system, with simple, effective and robust implementation in practical operation. The method provides a synchronization prediction function that can be used for predicting a sample-level timing offset that can be resulted in by an inaccurate sampling clock. The method permits the receiver to employ a low-precision crystal oscillator. In addition, the technical scheme provided by the present invention is flexible; it can be broadly applied in various systems that implement timing synchronization by transmitting the same signal block repeatedly. The predicted result provides frequency offset information as well, part of which can be put to use in correcting estimates of frequency offset.
To solve the above-mentioned problems, the present invention provides a coarse timing synchronization acquisition method in a mobile communication system, which is applied for timing synchronization acquisition of a mobile communication system that transmits the same signal block repeatedly. The method comprises the following processing steps:
Further, in said step C:
reliable timing synchronization positions consecutively approaching the predicted value of synchronization position of the frame may mean:
the timing synchronization positions of Lacq consecutive frames are reliable, and the difference between the timing synchronization positions and the predicted timing synchronization position of the current frame is less than or equal to Tacq; where Lacq is the minimum number of consecutive reliable frames of timing synchronization acquisition; and where Tacq is the threshold for timing synchronization acquisition decision;
determining the timing synchronization position of last frame to be the acquired timing synchronization position of current frame may mean:
the acquired timing synchronization position is the timing synchronization position of the last frame of the Lacq consecutive frames.
Further, this method may be applied in mobile multimedia broadcast systems.
Further, said step A may further comprise: computing the synchronization metric of DCA values;
in said step B, the position that corresponds to the maximum DCA value may be the position at which the synchronization metric is the local peak in the time slot.
Further, said synchronization metric may be a power P of the DCA value.
Further, said step B may further comprise: outputting a reliability indicator of timing synchronization position of current frame; the reliability indicator being used to signify whether the timing synchronization position of the current frame is reliable;
in said step C, all timing synchronization position of Lacq consecutive frames may being reliable means: reliability indicators of all the Lacq consecutive frames show that timing synchronization positions of these frames are reliable.
Further, in said step B, a timing synchronization position and a reliability indicator may be obtained according to the following formula:
the timing synchronization position of current frame is:
Vk is an estimate of timing synchronization position in a time slot, namely, U is an estimate of timing synchronization position corresponding to the maximum peak among the M local peaks.
A reliability indicator of timing synchronization position of the current frame is:
where Tdec is the threshold for timing synchronization reliability decision.
Further, in said step B, a timing synchronization position and a reliability indicator may be obtained according to the following formula:
a timing synchronization position of current frame is:
U=V
└M/2┘;
a reliability indicator of a timing synchronization position of current frame is:
where, local peak positions are normalized; Tdec is the threshold for timing synchronization reliability decision, V is a timing synchronization position estimate, and max(•) is a “maximum” operand.
Further, in said step B, a correlation envelope detection algorithm or maximum correlation value detection algorithm may be adopted to search for a local peak of each frame and the timing synchronization position estimate corresponding to the local peak.
Further, in step A, calculating the DCA value may mean:
calculating the initial value of DCA first:
where, r(k) is a digital baseband received signal, N is a length of DCA, ND is a number of delayed samples of delay correlation, k0 is a beginning time instant of timing synchronization acquisition, and r*(•) denotes a complex conjugate operation;
based on the initial value of DCA, starting from the beginning time instant of timing synchronization acquisition, a DCA value at each sample is computed iteratively:
θi=θi−1+r(k0+i+N−1)r*(k0+i+ND+N−1)−r(k0+i−1)r*(k0+i+ND−1),
where, i=1, 2, 3, . . . K.
Further, in said step A, DCA value may be calculated according to the following formula:
where, r(k) is a digital baseband received signal, N is a length of DCA, ND is a number of delayed samples of delay correlation, k0 is a beginning time instant of timing synchronization acquisition, and r*(•) denotes a complex conjugate operation; i=1, 2, 3, . . . K.
Further, in said step C, a predicted value of timing synchronization position of the next frame being obtained by using the L timing synchronization positions may mean:
adopting a robust least-square linear fitting algorithm to predict the timing synchronization position of the (L+1)th frame, by utilizing said L timing synchronization positions.
Compared with existing techniques, the method of the present invention is characterized by the flexible design of timing acquisitions, by a less demanding computation for a timing synchronization search, high accuracy of synchronization acquisition results, good robustness of a receiver system. The method completely avoids complex modulo or division operations that would be required in computation in an existing technique, and can be conveniently carried out in either hardware or software.
The technical scheme of the present invention will be illustrated in further details with reference to the drawing and embodiments hereinafter.
A core idea of the present invention is to search for correlation peaks within a certain time range by utilizing identical signal blocks at each time slot of physical frames and adopting a delay correlation method, so as to acquire the estimate of the receiver's coarse timing synchronization. In order to make the coarse synchronization result more robust and reliable, a timing synchronization result at a current time slot is predicted, which can greatly decrease the probability of the appearance of abnormal values in coarse timing synchronization results, and which can increase the accuracy of the receiver's timing acquisition.
Main technical characteristics of the present invention are as follows: (1) a delay correlation technique is adopted; (2) a related accumulation value carries frequency offset information; (3) a simplified method of power computation is adopted; (4) search results are stable and reliable; (5) a synchronization prediction technique is adopted; (6) a receiver can work even if it has a crystal oscillator of low precision. Combining the above-mentioned characteristics, the present invention, when it uses a low-precision crystal oscillator, can acquire the coarse timing synchronization of the mobile communication system transmitting the same signal block repeatedly with low complexity in real time, so that the requirement of timing synchronization raised by this type of mobile communication systems is well satisfied.
In order to make the technical scheme and technical advantages of the present invention clearer, a terminal receiver of a China Mobile Multimedia Broadcasting system will be discussed, thus illustrating in detail a technical scheme of the coarse timing synchronization acquisition provided by the present invention, with reference to
Exemplary values of main design parameters related to the present embodiment are as shown in Table 1.
In practical application, the parameters are not limited to the above-mentioned values, which can be chosen based on actual situations.
A, calculating a timing synchronization metric
Step A comprises the following steps specifically:
A1, calculating a DCA value
starting from the beginning time instant of timing synchronization acquisition, iteratively computing the DCA values at each sample.
In order to reduce the converging time of timing synchronization search, an initial value of delay correlation accumulation can be computed first:
where, r(k) is a digital baseband received signal, N is a length of DCA, ND is a number of delayed samples of delay correlation, k0 is a beginning time instant of timing synchronization acquisition, and r*(•) denotes a complex conjugate operation.
Based on the initial DCA value, starting from the beginning time instant of timing synchronization acquisition, a DCA value at each sample is computed iteratively:
θi=θi−1+r(k0+i+N−1)r*(k0+i+ND+N−1)−r(k0+i−1)r*(k0+i+ND−1),
where, i=1, 2, 3, . . . K, and K is a number of DCA values at each time slot.
Of course, in order to reduce storage space needed for the computation, the directly related formula below can also be employed:
To be noted, a disadvantage of a direct computing method is the relatively large amount of computation.
A2, calculating a timing synchronization metric
A DCA value is a complex number and, for the convenience of searching for timing synchronization, it is needed to construct a timing synchronization metric, which is a function of a DCA value.
In the present embodiment, powers of the DCA values are taken as the timing synchronization metric, the power is:
P
i=θiθ*i.
The following simplified formula can also be adopted in order to reduce computational burden:
P
i=αmax(|real(θi)|,|imag(θi)|)+βmin(|real(θi)|,|imag(θi)|),
where, α+1,β+½ in this embodiment, real(•) is the real part of the complex number while imag(•) is the imaginary part thereof, and max(•) and min(•) are the maximum and minimum of their respective operands, respectively.
In a practical application, it is not required that the power be employed as timing synchronization metric.
B, updating synchronization data
Step B comprises the following steps specifically:
B1, searching for a local peak
M time slots are searched in each frame of timing synchronization, and M timing synchronization position estimates can be obtained. When the normal timing synchronization position is searched, said DCA value appears to be the maximum one, and said timing synchronization metric value is the peak point, namely the maximum or the minimum value. In this embodiment, a local peak is the power in step A2, and the position corresponding to the local peak of this power serves as the estimate of timing synchronization position. Supposing that there is one timing synchronization position among every K samples, a local peak searching unit needs to find one timing synchronization position estimate within the range of every K sample points according to the power.
The present invention supports but is not limited to any one of the two algorithms of searching for local peak, where one is a correlation envelope detection algorithm, and the other one is a maximum correlation value detection algorithm.
If the second detection algorithm is adopted, the result of a timing synchronization position estimate corresponding to the local peak of local peak searching is:
It is assumed that the local peaks searched within 10 consecutive frames in step B1 are as shown in Table 2.
B2, determining current timing synchronization
In each frame, a local peak searching unit obtains a plurality of timing synchronization position estimates, selects the most reliable one from the multiple synchronization position estimates as the timing synchronization position U of current frame, decides whether the position U is reliable at the same time, and can be but is not limited to calculating the reliability indicator f of the current timing synchronization position, wherein this reliability indicator is used to indicate whether the timing synchronization position of the current frame is reliable. A reliability indicator of “1” means the timing synchronization position is reliable, while a reliability indicator of “0” means the timing synchronization position is not reliable. The most reliable estimate value can be selected according to the actual situation, such as selecting the one at a middle position. If only one time slot is searched in a frame, selection is unneeded. The present invention supports but is not limited to two decision algorithms, one being a peak-based detection algorithm, and the other one being a peak position-based determining algorithm.
a, peak-based detection algorithm
A timing synchronization position of current frame is
U is a timing synchronization position Vk of maximum peak among M local peaks, and PVk is the power corresponding to position Vk.
A reliability indicator of a timing synchronization position of the current frame is
where, PU is the power corresponding to position U.
b, peak position-based determining algorithm
Supposing that timing synchronization position estimates Vk corresponding to local peaks of several time slots in current frame are sorted in ascending order, the timing synchronization position of the current frame is
U=V
└M/2┘,
where, └M/2┘ is the largest integer not greater than M/2.
A reliability indicator of a timing synchronization position of the current frame is
(Note: local peaks need to be normalized hereby.)
A peak position-based determining algorithm as in b is adopted in the present embodiment. According to the partial peak list shown in Table 2, timing synchronization position and reliability indicator are calculated, and results of 10 consecutive frames are shown in Table 3.
C, acquiring synchronization results
Step C comprises the following steps specifically as follows:
C1, predicting synchronization results
After L frames are searched, L timing synchronization positions are obtained. A timing synchronization position of current frame is the timing synchronization position U corresponding to frame synchronization positions obtained in step B2. A timing synchronization predicting unit will, by making use of these L timing synchronization positions, predict the timing synchronization position pL of the (L+1)th frame. Supposing that the L timing synchronization positions are respectively s0, s1, . . . , sL−1, a fitting curve can be obtained by using robustness least-square linear fitting algorithm, that is:
p
i
=ax
i
+b,
where, xi=i, i=0, 1, . . . . Then the predicted value of timing synchronization position of the (L+1)th frame is:
p
L
=aL+b.
In accordance with the above-mentioned prediction algorithm, and according to the synchronization positions from frame 0 to frame 7 shown in Table 3, a fitting curve can be obtained as follows:
p
i=196.1i+12260,
then the predicted value of a synchronization position of frame 8 is:
p
L
p
8=196.1×8+12260=13829.
According to the synchronization positions from frame 1 to frame 8 shown in Table 3, another fitting curve can be obtained as follows:
p
i=197i+12250.
Then the predicted value of a synchronization position of frame 9 is:
p
L+1
=p
9=197×9+12250=14023.
In practical applications, one is not limited to a least-square linear fitting algorithm to obtain a fitting curve; indeed any linear fitting algorithm that can eliminate abnormal values can be adopted for this purpose.
C2, calculating the synchronization results
when there are consecutive reliable timing synchronization positions close to the predicted values of timing synchronization positions of current frame, it is determined that the acquired timing synchronization position is the timing synchronization position of the last frame.
Specifically, when several frames are searched, a predicted value sequence of timing synchronization position sL, sL+1, . . . , a reliable indicator sequence fL, fL+1, . . . , and a timing synchronization prediction sequence pL, pL+1 . . . are obtained, when Lacq consecutive frames satisfy these conditions:
|sL+q−pL+q|≦Tacq,
and
f
L+q=1, where, q=0, 1 . . . , Lacq−1,
it is considered that a reliable timing synchronization position is acquired, and the outputting acquired timing synchronization position is that of the last frame: sL+q. Otherwise, the synchronization acquisition process will be continued.
Of course, it is not required to output the reliability indicator in practical applications, and when Lacq consecutive frames satisfy the following condition:
|sL+q−pL+q|≦Tacq,
and timing synchronization position of the Lacq consecutive frames are reliable, the outputting acquired timing synchronization position is that of the last frame of the Lacq consecutive frames: sL+q.
According to the above method, synchronization acquisition conditions in the present embodiment and data in the above tables, are calculated as:
|sL−pL|=|s8−p8|=5≦Tacq=32, and f8=1;
|sL+1−pL+1|=|s9−p9|=28≦Tacq=32, and f9=1;
then it holds that the receiver has acquired the timing synchronization position, which is sL+1=14051.
The above discussion, with reference to specific preferable embodiments, is a further illustration of the present invention. The invention should not be construed, however, as limited to the embodiments illustrated above. It should be understood by those skilled in the field that, without departing from the spirit of the present invention, any simple deduction or substitute is deemed to be included within the scope of the present invention, as defined by the claims.
The present invention can, under the condition that low-precision crystal oscillator is adopted, take full use of structural characteristics of physical frames of the mobile communication system transmitting the same signal block repeatedly, and can acquire coarse timing synchronization of such mobile communication systems fast, in real time and with low complexity of calculation, so as to satisfy well the requirements of timing synchronization raised by such mobile communication systems.
Number | Date | Country | Kind |
---|---|---|---|
200710124223.9 | Oct 2007 | CN | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/CN07/03944 | 12/29/2007 | WO | 00 | 3/19/2010 |