1. Field of the Invention
The present invention relates to Code Division Multiple Access (CDMA) for mobile communications, more particularly to a method and device for estimating carrier frequency offset (CFO) in subscriber terminals in Time Division Synchronous Code Division Multiple Access (TD-SCDMA) system.
2. Description of the Related Art
In the TD-SCDMA system, both base stations and terminals are supposed to send and receive data with nominal carrier frequency. However, there is always some deviation between the actual transmitting and receiving frequency and the nominal frequency because of the technology limitation of the devices. In a base station, the frequency accuracy of an oscillator, in which a carrier frequency offset in a base station required by a TD-SCDMA system is less than 0.05 PPM, may meet to the standard requirement, since limitations subjected by a base station from temperature, size, power consumption, cost and the like are smaller. In a subscriber terminal, however, a carrier frequency offset in a subscriber terminal required by the TD-SCDMA system is less than 0.1 PPM, the frequency accuracy of a selected oscillator may not meet the standard requirement, due to the limitations subjected by the subscriber terminal from temperature, size, power consumption, cost and the like are greater. Consequently, carrier frequency offsets need to be estimated in a subscriber terminal by using a known synchronous downlink pilot signal (SYNC_DL) sent by base stations, and according to the estimated value of carrier frequency offsets to perform automatic frequency control (AFC), so that carrier frequency offsets in the subscriber terminal may meet the standard system requirements.
In a method for estimating carrier frequency offset in prior art, SYNC_DL signals only send by a single base station are utilized in a subscriber terminal to estimate the carrier frequency offset. The base station has one or more antennas to transmit signals with the same carrier frequency.
In the TD-SCDMA system, a dedicated downlink pilot time slot (DwPTS) is set, a base station sends SYNC_DL signals at DwPTS of each frame. Codes by which a SYNC_DL signal is consisted of are used to distinguish adjacent cells using the same carrier frequency, that is, the adjacent cells with the same carrier frequency are distinguished by the SYNC_DL codes. There are 32 different SYNC_DL codes defined in the standard, length of each set of SYNC_DL codes is 64 chips. During the initial search of the cells, the subscriber terminal search the SYNC_DL codes at a specific frequency point used in a cell where it currently expects to resided in to implement the downlink synchronization, and then estimate the carrier frequency offset based on the SYNC_DL codes and adjust the carrier frequency according to the estimated carrier frequency offset value to meet the system requirements.
Provided that the system requires for carrier frequency offset estimation and rough adjusting, adjusted carrier frequency offset in a subscriber terminal is less than 1 kHz, then the subscriber terminal has only received, at a corresponding frequency point, the SYNC_DL signals transmitted by one base station or received SYNC_DL signals of the cell currently expected to resided in are much greater than the SYNC_DL signals of other cells using a same frequency. If the signal power ratio is over 6 dB, it is possible for the adjusted carrier frequency offset of a subscriber terminal to meet to the system requirement. However, if the power difference between the simultaneously received SYNC_DL signals respectively in two cells with the same frequency is small, the SYNC_DL signals transmitted by one cell will interfere those transmitted by another cell. For example, when the power ratios of multiple SYNC_DL signals received by a subscriber terminal at the edge of several cells are within 3 dB, there is about 10-20% probability with the existing methods for carrier frequency offset estimation and adjustment that the carrier frequency offset of the subscriber terminal can't meet the requirement of the system.
In the TD-SCDMA system, because of the short length of a SYNC_DL signal, which is only 64 chips (50 μs), the carrier frequency offset after adjustment permitted by the standard is still greater, for example, less than 1 kHz. Even so, as shown in
In view of this, an object of the present invention is to provide a method and a device for estimating carrier frequency offset in a subscriber terminal to make the carrier frequency offset in a subscriber terminal to meet the system requirements and thus improve the probability of success for the initiate search for a cell.
To achieve the above object, the present invention provides a method for estimating carrier frequency offset in a subscriber terminal, said method comprises steps:
Said determining number of effective base stations from which more than one signals are received by a subscriber terminal in step A comprises steps:
Said signals are synchronous downlink pilot signals, and said step A1 further comprises steps:
Said step A2 further comprises steps
The method further comprises a step before said step A: reading vector data of 128 chips while receiving synchronous downlink pilot signals at the beginning of a downlink pilot time slot.
The method further comprises a step before said step B: multi-path combining signals of each base station.
Said step of multi-path combining signals of each base station comprises steps:
The invention also provides a device for estimating carrier frequency offset, which comprises at least:
Said device further comprises a multi-path combining module for multi-path combining the signals of each base station, and then outputting the multi-path combined signal to the combining module, if the base station number is greater than 1
According to a solution of the invention, when a subscriber terminal simultaneously receives SYNC_DL signals sent by a plurality of base stations (more than one) and their power differences are small, it is possible to make the carrier frequency offset in a subscriber terminal, after automatic frequency control, to meet the system requirement and thus improve the probability of success for the initiate search for a cell.
The method for estimating carrier frequency offset provided by the invention is especially suitable to a condition that the strengths and SIR of SYNC_DL signals received by a mobile user from several base stations are approximate among them, i.e., a mobile user is located at the joint edge of several cells, and with the method for carrier frequency offset estimation provided by the invention, the probability of success for the initiate search for a cell can be improved greatly. Even when the SYNC_DL signal received by a mobile user from one base station is much stronger and those received from other base stations are not so strong, the accuracy of adjustment for carrier frequency offset also can be improved by using the method for estimating carrier frequency offset provided by the invention. With respect to a method utilizing the SYNC_DL signals from one base station, the closer of the SYNC_DL signal strength and SIR from several base stations, the more improvement the accuracy of carrier frequency offset adjustment.
Hereinafter, detailed description will be made to a method and device for estimating carrier frequency offset in a subscriber terminal in TD-SCDMA system according to the embodiments of the invention, with reference to accompanying drawings, so that the invention will be well understood by those skilled in the art.
The first step, in step 201, determine the SYNC_DL signals simultaneously received by a subscriber terminal from several effective base stations with the same frequency and the respective arriving time of their main path.
Then, in step 302, shift complex correlation with the data received at DwPTS to 32 groups of local SYNC_DL codes respectively, shifting from 0 to 63, i.e., calculate the correlation power values of each point.
In step 303, the correlated power values for K frames are averaged (K≦1).
The signal powers for two successive frames arrived at a subscriber terminal may be different because of fading, etc,. To ensure the accuracy of the data obtained in a subscriber terminal, it is preferable to calculate the correlation power by using more frames of data and then average them.
Then, in step 304, averaged correlation power values corresponding to SYNC_DL codes of each group in 32 groups are fined out; the arriving time for the SYNC_DL signals corresponding to the SYNC_DL codes are determined by the SYNC_DL codes, i.e., the arriving time for the main path of the SYNC_DL signals corresponding to the SYNC_DL codes; and, shift complex correlation with data r to the SYNC_DL codes to obtain the peak value of correlation power, which indicates the synchronization of received signals and the corresponding SYNC_DL codes.
In step 305, the peak values of the averaged correlation power are ordered in a falling sequence, and the first Nmax ones are selected, denoted as Pmax(n), n=1, 2, . . . , Nmax. Nmax is a predefined value, and represents the max number of base stations in which a subscriber terminal receives the SYNC_DL signals from different base stations simultaneously.
In step 306, since a subscriber terminal is not able to know how many SYNC_DL signals are received simultaneously from base stations when it receives SYNC_DL signals, thus it is needed to determine the number N of effective base stations from which the SYNC_DL signals have been received simultaneously, based on a threshold, set the threshold VT, and firstly decide whether the Nmaxth peak value of the averaged correlation power meets following condition:
If yes, N=Nmax; otherwise, continue to decide whether the (Nmax−1)th peak value of the averaged correlation power values meets following condition:
If yes, N=Nmax−1; otherwise, the rest may be deduced by analogy, until the number N of effective base stations, from which a subscriber terminal receives the SYNC_DL signals simultaneously, is determined, where 1≦N≦Nmax. The SYNC_DL codes sequence number and the peak values point corresponding to each of the N base stations are recorded as Sync_IDn and Peakn respectively , where n=1, 2, . . . , N.
In the embodiment of the invention, Nmax=3, and VT=4.
The second step, in step 202, combine multipath signals of SYNC_DL received from each effective base station respectively.
Then, the complex SYNC_DL codes numbered as Sync_IDn are operated, and the complex values of SYNC_DL codes are marked as sn, where the kth element is sn,k, k=0, 1, . . . , 63. 64+2M points of data are read beginning from M chips ahead of Peakn in the DwPTS of the received signals, marked as a vector e, where the kth element is ek, k=1, 2, . . . , 63+2M.
e is shift complex correlated with sn, and the result corn is obtained:
then the SYNC_DL signals of each paths are combined and a vector Uu is obtained
where, Sn,k* indicates a conjugate of sn,k, and corn,i* indicates a conjugate of corn,i.
The physical meaning of the formula is to perform Max Ratio Combining of multi-path SYNC_DL signals with different time delay after eliminating the phase differences both between the symbols and between the time delay paths, in which, received signals e multiplied by the conjugate of sn corresponds to eliminating the phase difference between symbols, and multiplied by the conjugate of corn corresponds to eliminating the phase difference between time delay paths.
The process of step 202 of multi-path combining SYNC_DL signals received from a subscriber terminal will be further described with reference to
In step 401, 64+2M points of data are read beginning from M chips ahead of Peakn in the SYNC_DL signals received at DwPTS, i.e., 64+2M received data is read beginning from M chips ahead of the position of the main path, and a complex sequence of 64+2M is obtained, e(k), k=0, 1, . . . , 63+2M , marked as a vector e , where 64 is the length of the SYNC_DL code.
In step 402, the SYNC_DL codes numbered SYNC_ID are complex operated and the complex SYNC_DL codes are obtained , i.e., the local complex SYNC_DL code sequence s(k), k=0, 1, . . . , 63, marked as sn.
Then, in step 406, conjugating operation of each element of the sequence s(k) is performed, and a sequence s* (k)is obtained.
Because e(k+1) includes a phase shift, in order to eliminate the phase shift between subsequent two phases, in step 403, let d(0,k)=e(k), k=0, 1, . . . , 63; in step 404, let d(1,k)=e(k+1), k=0, 1, . . . , 63, . . . , and in step 405, let d(2M,k)=e(k+2M), k=0, 1, . . . , 63.
In step 407, each element of the sequence d(i,k) (i=0,1, . . . , 2M , k=0, 1, . . . , 63) obtained in steps 403˜405 is multiplied by the corresponding element of the sequence s* (k) obtained in step 406 respectively.
In step 408 , all elements of the sequence obtained in step 407 are summed up and then conjugated. In steps 403 ˜step 408, shift multiplex correlation of e and sn* are performed, and a result corn is obtained as follows:
here, ei+k is the d(i,k) and the e(i+k).
In step 409, the result of shift multiplex correlation is multiplied by all elements of the sequence, then, in step 410, all 2M+1 sequences, i.e., multi-path combining the SYNC_DL signal of each path, are summed up and a vector Un is calculated as follows.
The value of M is integrally dependent on both the time delay extension of the received signals and the effect of noise. In the embodiment of the invention, M=4.
The third step, in step 203, multi-path combined SYNC_DL signals of each base station are further combined.
combining the multi-path combined SYNC_DL signals of each effective base station is to combine the calculation results of the SYNC_DL codes, and equal gain combination is used:
and a weighting combination may also be used
where the weighting coefficients Pn can be obtained by the following formula:
Because there is coherent property among the SYNC_DL signals from each effective base station while noncoherent property among noises, it is possible to improve the SNR (or SIR) of available signals by combining the multi-path combined SYNC_DL signals from each effective base station so that the SYNC_DL signals received by a subscriber terminal may be efficiently utilized.
Based on theory formula, the sequence combined by the SYNC_DL signals from each effective base station can be approximately represented as follows:
Uk=Ak·exp(j·2π·Δf·k·TC+φ0)+Nk,k=0,1, . . . , 63,
where, Uk is the kth symbol value in the combined SYNC_DL signal sequence, k=0, 1, . . . , 63;
Ak is magnitude of combined SYNC_DL signals valid in the kth symbol;
exp(x) indicates exponent operation of complex x;
Δf is carrier frequency offset, which is the exact value in theory;
TC is pulse width of each chip;
φ0 is initial phase of the SYNC_DL signals received by a subscriber terminal, range of which is between 0˜2π;
Nk is the complex noise value or interference noise value in the kth symbol.
The fourth step, in step 204, rough estimation values of carrier frequency offset are calculated by using the combined signals.
Based on the theoretical formula in step 203, it can be seen from the expression of Uk that in the combined SYNC_DL signal sequence, a variation rate of phase between symbols is associated with a carrier frequency offset Δf. Therefore, a carrier frequency offset can be roughly estimated by using the combined SYNC_DL signal sequence. In one embodiment of the invention,
Where Δ{circumflex over (f)} is a rough estimation of the carrier frequency offset, TC is an pulse width of each chip, and angle (x) is an operation of the phase angle of complex x. Values P and L will meet formula: P+L≦64, and P≧1, L≧1.
Above formula indicates the effect of noise (or interference) can be minimized by using a phase difference between two symbols spaced by L to estimate the carrier frequency offset, and then averaging the estimated values for P times.
In another embodiment of the invention,
Where Δ{circumflex over (f)} is rough estimated value of the carrier frequency offset, TC is pulse width of each chip, and angle (x) presents the operation of the phase angle of complex x. Values P and L will meet formula: P+L≦64, and P≧1, L>1.
In order to minimize calculating amount of repeatedly calculating phase angle, the above formula indicates to firstly sum up the phase differences between two symbols spaced L, and then calculate phase angles to estimate the carrier frequency offset.
Due to the phases of two consecutive symbols are close together, an angle jitter of the multiplied phase angle will be greater. Therefore, it is preferable to multiply two symbols spaced L to avoid the angle jitter to be more great.
In the embodiment of the invention, P=L=32.
The invention also provides an device for estimating carrier frequency offset in a subscriber terminal, which corresponding to the steps shown in
UE carrier frequency offset is roughly adjusted by using the method provided by the invention. Parameters are selected as follows: K=4, VT=4, M=4. The simulation outcome is shown in
In summery, the description above is only the preferred embodiments of the invention but will not limit the protection scope of the invention.
Number | Date | Country | Kind |
---|---|---|---|
03149619.9 | Aug 2003 | CN | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/CN04/00898 | 8/4/2004 | WO | 10/18/2006 |