Method for improving the reception of a CDMA receiver, and corresponding CDMA receiver

Abstract: The invention relates in general to systems for digital transmission, for example digital mobile radio systems such as the UMTS (Universal Mobile Telecommunications System) and in particular to the receiving-end estimation of parameters which are required for retrieving the data. An object of the invention resides, in particular, in indicating an efficient way which, particularly in the case of CDMA systems with power control, ensures estimates of higher quality than with methods of the prior art, and simultaneously largely avoids a substantial increase in the complexity of system implementation. For this purpose, the invention proposes a method for improving the reception of a CDMA receiver based on a parameter estimation method, as well as a correspondingly adapted CDMA receiver, in which the estimates are produced by processing received values (y i [k]) provided with correction factors (s), the correction factors (s) serving to equalize different power levels of the received signal.

Description: The invention relates in general to systems for digital transmission, for example digital mobile radio systems such as the UMTS (Universal Mobile Telecommunications System), and in particular to a method for improving the receiving-end estimation of parameters which are required for retrieving the data, and a receiver of a CDMA (Code Division Multiple Access) system adapted for carrying out the method. As the person skilled in the art in this field is well aware, examples of parameter estimation methods are methods for estimating a frequency offset between modulator and demodulator, which is subsequently compensated, and/or channel estimation methods which are required, for example to set the coefficients for a rake receiver so as to eliminate channel distortions. For the purpose of better comprehension of a rake receiver, reference may be made, in particular, to the publication by A.J. Viterbi, "CDMA", Addison-Wesley, Reading, Mass., 1995. Particularly in the case of digital transmission over dispersive channels, for example via a mobile radio channel, the transmitted signal is distorted and disturbed by noise. The channel distortions necessitate special concepts for retrieving the transmitted data; for example a receiving-end sequence estimation with the aid of the Viterbi algorithm is carried out in a TDMA (Time Division Multiple Access) transmission, as is used in the case of GSM systems. When transmission is done using the CDMA (Code Division Multiple Access) method, as happens in third-generation mobile radio systems, for example in UMTS systems, retrieval is possible by means of a rake receiver, for example. Methods which are more complicated, but also more powerful, are based, in particular, on a multiuser detection, supplementary reference being made for this purpose to the publication by S. Verdu, "Multiuser Detection", Cambridge University Press, Cambridge, 1998. The reception methods conventionally require estimates of the channel pulse response. In most CDMA systems known pilot symbols whose corresponding received values can be used to carry out the estimation are transmitted to the receiver for this purpose. A multiplicity of channel estimation methods in the literature are based on this concept. The power efficiency of the rake receiver and/or of the multiuser detection is greatly influenced by the quality of the channel estimates. Consequently, CDMA systems are strongly dependent on the use of an adaptive power control in order to achieve a desired capacity of the overall system and transmission quality. The term "power control" in this case covers both the fast power control for reaction to fading, and the longer-term adaptation of the transmit power which is required, for example in the case of a change in data rate. With all currently existing CDMA systems it is usual to undertake such power control in blockwise fashion, such that the transmit power is held constant in each case over individual blocks of symbols and is changed discontinuously between individual blocks. Since each channel estimate is thereby generally influenced by received values from different blocks, different power levels being assigned to the various blocks in systems with power control, this results in the need to combine symbols with different transmit power levels. if these different transmit power levels are not taken into account explicitly, this can therefore lead to a lower quality of the channel estimation results. However, this problem is ignored in systems of the prior art, and/or a palpable degradation of the channel estimate caused by this discontinuity is accepted. International Patent Publication W0 00/33472 discloses method for improving the reception of a CDMA receiver based on a parameter estimation method, in which the estimates are produced by processing received values (y i [k]) provided with correction factors (s), the correction factors (s) serving to equalize different power levels of the received signal. The present invention is characterised over the disclosure of W0 00/33472 in that a discontinuity of transmit power levels and the correction factors (s) resulting therefrom are estimated by statistical hypothesis testing and are assigned to the received values (y i [k]). Preferred embodiments of the invention advantageously avoid the indicated problems of the prior art, in particular the degradation described. Preferred embodiments of the invention also advantageously provides an efficient way which, particularly in the case of CDMA systems with power control, ensures estimates of higher quality than when use is made of methods according to the prior art, and simultaneously largely avoids a substantial increase in the complexity of system implementation. The solution according to the invention is supplied by a method, a CDMA receiver and an implementation program having the features of Claims 1, 10 and 14, respectively. Advantageous and/or preferred embodiments and/or developments are the subject-matter of the respective subclaims. For the purpose of parameter estimation methods which are as efficient as possible, in particular for the purpose of estimating a channel pulse response in the case of a CDMA transmission, the invention provides for estimating a correction factor with which different power levels of various receiving blocks can be equalized for the parameter estimate. Discontinuous variations in the received signal of the parameter estimate can therefore be avoided by estimating and compensating a power control factor. Since power control information is transmitted, for example in the case of UMTS systems, only with a relatively low reliability, it is provided in a preferred development to use methods of the statistical hypothesis test to estimate power control factors. An advantage of this is that analytical expressions can be specified for the method on the basis of the description of the attenuation of transmission channels in mobile radio, and also of superimposed noise as complex Gaussian processes. The estimate of the power control factor is, moreover, preferably based on a maximum a-posteriori approach, a decision going to the power control factor for which the probability of occurrence is a maximum in terms of measurement of prescribable observed variables. According to the invention, different observed variables are specified for individual preferred estimation methods, it being possible to estimate advantageously both with the aid of received values which are assigned the pilot symbols, and with the aid of received values which are assigned to unknown data. Furthermore, the method according to the invention is also suitable for profitable application in estimating further parameters such as, for example, a frequency offset between a transmitter and a receiver. The method according to the invention, and the receiver adapted according to the invention for carrying out the method therefore ensure a rise in the quality of parameter estimation at the receiving end, such as, in particular, channel estimation in CDMA systems with power control. Since, particularly with the use of a rake receiver to reconstruct the transmitted data, the detection is substantially based on channel estimates and/or estimates of other parameters such as, for example, the frequency offset, the improvement in parameter estimation therefore leads to a reduction in the bit error rate of the transmission system, which therefore entails a possible saving in the required transmit power, and thus a higher overall capacity. By contrast, in practical implementation this leads only to a slight extra outlay. The invention is described in detail below by way of example with the aid of preferred exemplary embodiments and with reference to the attached drawings, in which: Figure 1 shows a simplified block diagram of a rake receiver with maximum ratio combining, taking explicit account of the power control, Figure 2 shows a block diagram corresponding to Figure 1, but without the account of the power control, Figure 3a shows an input signal for a parameter estimate without compensation of the power control, Figure 3b shows an input signal for a parameter estimate with compensation of the power control, and Figure 4 shows an input signal of a rake finger for estimating a power control factor. For a better understanding of the approach according to the invention, and of the system model on which the following description is based, reference is firstly made to Figure 2, which illustrates a simplified block diagram of a rake receiver for a transmission with DS-CDMA (Direct Sequence Code Division Multiple Access) for a linearly distorting channel. As is known per se to the person skilled in the art, after traversing a receiving filter 1 an input signal coming from an antenna is processed for this purpose in paths, the so-called fingers F 0 , F i , F L-1 of the rake receiver, comprising a plurality of devices, in particular for delay 2, despreading 3 and channel estimation 4. The time-discrete output signal y[k] of the rake receiver in accordance with Figure 2 is given in this case by L denoting the number of the fingers F 0 , F i , F L-1 of the rake receiver. If propagation times which occur are ignored in order to simplify the representation, the output signal y i [k] of the ith finger F i in accordance with Figure 2 is given by y i [k] = h i [k]·a[k]·s+n i [k], i∈{0, 1,..., L-1) The sequence a[k] in this case includes the transmitted symbols which, depending on the modulation method used, can either be truly real or complex. The parameters h i [k] constitute time-dependent weighting factors relating to the ith rake finger F i , which are given by the channel pulse response, the parameter h and i [k] in equation (1) denoting the estimated channel weighting factor relating to the ith rake finger F i , which is determined by the channel estimating device 4 in Figure 2. The parameters n i [k] denote interfering signals of the ith rake finger F i , which are conventionally composed of different components, for example white noise and MAI (Multiple Access Interference). The factor s constitutes a weighting factor which influences the power of the received signal, that is to say exerts power control (fast power control or longer term power matching), and is constant for the receiving blocks of duration K. As is further known to the person skilled in the art in this field, both h and i [k], and n i [k] can be assumed to a good approximation in many mobile radio applications to be complex average less Gaussian processes, this assumption forming the basis of the following description. In the case of the UMTS systems, for example, the power control factor s is suddenly communicated to the receiver, but the corresponding reception symbols are not highly reliable, since the requirement for a short delay renders it impossible to use a powerful channel coding for these symbols. In general, the receiver therefore has no precise knowledge of the factor s, and so it is therefore necessary, rather, essentially to consider all possibilities for s. Their (a-priori) probabilities can be determined with the aid of the reception symbols containing the power control information, and of the corresponding symbol error probabilities. For the sake of simplicity, the symbols of each block are enumerated with a constant power control factor s at discrete times k∈{0, 1, ..., K-1}. It is further assumed that each block K p includes pilot symbols, known to the receiver, at the positions k 0 , k 1 , ..., k Kp-1 . Channel estimation for determining the factors h i [k] can be undertaken with the aid of the received values y i [k 0 ], y i [k 1 ], ..., y i [k Kp-1 ], in which case various methods based, for example, on filtering and/or interpolation are known for this purpose from the literature, for example from publications by J.K. Cavers, "An Analysis of Pilot Symbol Assisted Modulation for Rayleigh Fading Channels", Transactions on Vehicular Technology, VT-40:686-693, Nov. 1991; G. Auer et al., "Adaptive Mobile Channel Prediction for Decision Directed Rake Receivers", IEE Colloguium on Adaptive Signal Processing for Mobile Communications Systems, pages 13/1-13/5, Oct. 1997; H. Andoh et al., "Channel Estimation Filter Using Time-Multiplexed Pilot Channel for Coherent Rake Combining in DS-CDMA Mobile Radio", IEICE Trans. Commun., E81-B(7): 1517-1526, July 1998. It is also possible with these methods to incorporate into the channel estimate received values which correspond to data symbols, doing so by applying the principles of decision feedback. This results in an increase in the quality of estimation. It is known to the person skilled in the art that received values from different blocks generally influence the current channel estimate, the consequence of which is that symbols are combined with various power control factors s, as may be seen from Figure 3b) for example. As already mentioned above, this can result in a degradation in the adaptation of the channel estimating filter or in a lower quality of the result of channel estimation. Furthermore, a similar problem also occurs, for example in frequency offset estimation methods which operate with estimates for the channel autocorrelation sequence, reference being made for this purpose by way of example to the publication by W.-Y. Kuo et al. "Frequency Offset Compensation of Pilot Symbol Assisted Modulation in Frequency Flat Fading", Transactions on Communications, COM-45: 1412-1416, Nov. 1997. In a departure from the parameter estimation methods according to the prior art, in which these problems are ignored or a resulting degradation is accepted, the approach of the method according to the invention is to estimate the power control factor s and thereafter compensate it in order, in particular, to avoid abrupt discontinuities in the input signal of the parameter estimate, as shown in Figure 3b). The estimation to be carried out for this purpose is described below with reference to Figures 1 and 4. It is generally impossible to make direct use of the factor s communicated to the receiver by the transmitter, since the symbols which include this information are mostly not sufficiently reliable. It is therefore provided in order to estimate the factor s to apply methods of the statistical hypothesis test as disclosed, for example, in the publication by A.D. Whalen, "Detection of Signals in Noise", Academic Press, New York, 1971. In accordance with Figure 1, following estimation, carried out in a device 5, of the correction factor s, that is to say substantially of the power control factor, values 1/s·y i [k], in which abrupt changes in power no longer occur, are made available to the device 4' for parameter estimation (for example channel estimation), and this improves the result of parameter estimation. In order to carry out estimation of the factor s, it is preferred to consider the first M output values of each rake finger F 0 , F i , F L-1 per block n, n+1. These values are combined in reception vectors y i = [ y i [0] y i [1] ... y i [ M - 1]] T , i ∈ {0,1,..., L - 1}, (.) T denoting transposition. On the basis of this, interest now attaches, in particular, to the overall density of y i the case of M = 2 being described below by way of example. Since h i [k] and n i [k] can be taken as complex Gaussian processes, the result for given values of s, a[0] and a[1] (compare equation (2)) as conditional density for y i is a two-dimensional complex Gaussian density with the probability density function of p ( y i | s , a [0], a [1]) = 1 π 2 |det( M i )| exp(- y H i M -1 i y i ), (.) H denoting Hermitian transposition of the matrix, and the autocorrelation matrix M i being given by The following definitions hold here: where and σ 2 / i, σ 2 / h,i and σ 2 / n,i respectively representing the variance of the output signal of the ith rake finger F i , the variance of the channel weighting factor relating to the ith rake finger F i and the noise variance at the ith rake finger F i and the ϕ h i h i [κ] and ϕ n i n i [κ] representing the autocorrelation sequences of h i [k] and n i [k] respectively, where and The determination of the autocorrelation functions and variances σ 2 / h,i and σ 2 / n,i can be performed in the receiver, for example with the aid of the channel estimates. It is to be noted that σ 2 / i is a function of the power control factor s to be tested, and the correlation coefficient ρ i is, furthermore, a function of a[0] and a[1]. Furthermore, consideration will be given for what follows to the case of white noise n i [k], that is to say, ϕ n i n i [κ] = 0, κ ≠ 0 which case mostly obtains in practice to a good approximation. This leads to ρ i = ϕ h i h i [1]a[0]a[1]s 2 σ 2 i . A vector (or else scalar) u i = ƒ (y i ) is formed before estimation of the power control factor s. It is generally known to be necessary to assume that no pilot symbols are assigned to the time steps 0, 1, ..., M-1 under consideration, and that the corresponding transmission symbols are therefore unknown for the estimation. In other words, no use can be made of decided data, since obtaining them requires knowledge of channel estimates, whereas a is estimated upstream of the channel estimation in the current block. The vectorial or scalar function f(.) serves the purpose of keeping the influence of the unknown data symbols as slight as possible, it also being possible, as described at a later junction, to cover the case of known pilot symbols by suitable definition of f(.). Although it is also possible in principle to use criteria known to the person skilled in the art to estimate the factor s, a maximum a-posteriori approach is used below for the exemplary description. The estimated factor is therefore determined in accordance with p (s|u 0 , u 1 ,..., u L-1 ) denoting the conditional probability density function of s with u 0 , u 1 , ..., u L-1 being known. Bayes' theorem is used to obtain the expression which is accessible for numerical evaluations, that is to say: p ( s |u 0 ,u 1 ,...,u L-1 ) = p (u 0 ,u 1 ,...,u L-1 |s) Pr (s) p (u 0 ,u 1 ,...,u L-1 ) . Equation (13) is to be evaluated in this case for all possibilities, that is to say hypotheses for the factor s. Here, p (u 0 ,u 1 ,...,u L-1 ) can be neglected as an irrelevant factor for determining the most likely value, and the calculations of the a-priori probabilities Pr (s) can be carried out in a simple way with the aid of the error probabilities relating to the symbols containing the power control information. Furthermore, the overall density p (u 0 ,u 1 ,...,u L-1 |s) to be calculated can be simplified. Since the channel coefficients h i [k] for various rake fingers F 0 , F i , F L-1 can be modelled in a fashion statistically independent of one another, as can the interference processes n i [k], it is possible with a suitable selection of the function f(.) to assume the validity of the factorization p (u 0 ,u 1 ,...,u L-1 |s) = p (u 0 | s )· p (u 1 |s)·...· p (u L-1 |s). The function f(.) must be known in order to determine the density functions p (u i |s). Three preferred possibilities for forming f(.) are discussed below. A first possibility, without the need to know training frequencies on the part of the receiver, consists in forming the squares of absolute values component by component, that is to say Furthermore, the number of the symbols is selected as M=2, and binary antipodal modulation is assumed, such as is used, for example, in UMTS systems, that is to say a [k] ∈ {±1}. The representation is selected for calculating p(u i |s). The calculation of an overall probability density function of the squares of absolute values of two correlated complex Gaussian random variables, the calculation of p (u i =f(y i )|s, a[0], a[1]) in the case under consideration, proves to be a standard problem for the case of a two-dimensional complex Gaussian density for y i for given s, a[0] and a[1]; its solution is known from the literature, for example from publications by M. Schwartz et al. "Communication Systems and Techniques", McGraw-Hill, New York 1966 and by C.W. Helstrom, "The Resolution of Signals in White, Gaussian Noise", Proceedings of the IRE, 43:1111-1118, Sept. 1955. With σ 2 / i and ρ i in accordance with equation (6) and equation (10) respectively, and the modified Bessel function of first kind of zero order the result is compare, for example, the abovementioned publication by M. Schwartz et al. In accordance with equation (10), the absolute value of ρ i is independent of a[0] and a[1] for binary antipodal modulation and white noise. This results in accordance with equations (16) and (17) in where σ 2 i = σ 2 h,i s 2 + σ 2 n,i , and |ρ i | = ϕ h i h i [1] s 2 σ 2 i . A simplification can be undertaken given low signal-to-noise ratios, which occur as a rule at the outputs of the individual rake fingers F 0 , F i , F L-1 . Use may be made in this case of the approximation ln (I 0 (χ)) ≈ χ 2 /4 which is valid for small arguments and which leads to supplementary reference being made for this purpose by way of example to the publication by D. Raphaeli, "Noncoherent Coded Modulation", Transactions on Communications, COM-44:172-183, Feb. 1996. Finally, using equation (13) and equation (14) the estimated power control factor s is yielded as that value which maximizes the expression As a second possibility, the squares of the absolute values of the individual received values can be summed to form the function f(.), that is to say M = 2 is chosen again by way of example, and the case of binary antipodal modulation is considered. The conditional density function of u i =f(y i ) for a given s can be derived starting from the publication by M. Schwartz already referred to. Once again, the expression independent of a[0] and a[1] is yielded for the conditional density p (f(u i ) |s,a[0],a[1]), leading finally to equations (20) and (21) holding, respectively, in turn for σ 2 / i and |ρ i |. Starting therefrom, it is therefore possible to use equations (13) and (14) to specify directly the target function to be maximized as a function of s, more specifically A further preferred possibility consists in directly using y i to estimate s, that is to say, u i = y i in the case of the occurrence of training symbols, known at the receiving end, at the points k∈{0,1,...,M-1}. The case of M = 2 and a[0] = a[1] = 1 will be considered as an example. The resulting conditional density function is where σ 2 i = σ 2 h,i s 2 + σ 2 n,i , ρ i = ϕ h i h i [1]s 2 + ϕ n i n i [1] σ 2 i . Equations (13) and (14) can therefore be used to set up the target function, which is to be maximized as a function of s, as Whereas the above exclusively considered the case of considering only signal values of the new receiving block n+1 to estimate s, that is to say after a change in the power control factor s, an estimate can also be made by additionally considering received values of the old block n. As an example, the use of received values from an old and a new receiving block n and n+1, respectively, are considered (Figure 4) below for this purpose, that is to say in the case of M=2 the received vector relating to the ith rake finger F i is formed as It is assumed, furthermore, that the compensation of the power control factor s has already been undertaken for the old received value. A Gaussian density in accordance with equation (4) is yielded once again as conditional density p (y i |s, a[-1], a[0]), the result for the autocorrelation matrix now being Here, s 0 denotes the already determined power control factor of the old block n. As a person skilled in the art in this field will recognize, it is therefore possible to use this result once again to derive target functions for estimating s for the three possibilities discussed above which are based on function f(.) and binary antipodal modulation. In a practical implementation, in particular inside a receiver of a CDMA-based mobile radio system, the method entails only a slight extra outlay and is started in relation to any period of the power control. After the storage of the samples y i [k] for all L rake fingers, the latter are therefore additionally preprocessed with the formation of the function f(.) based, in particular, on one of the three preferred possibilities described above. Consequently, it is subsequently possible in conjunction with consideration of the a-priori probability Pr(s) for the occurrence of a power control factor s which is known per se from a transmitted power control command and the uncoded error rate then to set up a target function for all possible power control factors s, it being possible to use samples from a receiving block n+1 or, alternatively, received values from different blocks n, n+1 (Figure 4): for example, when selecting the number of symbols to be considered in relation to M=2, either the samples y i [0] and y i [1] from the new receiving block n+1 or, alternatively, the samples y i [-1] and y i [0] from different receiving blocks. The estimate for the power control factor s actually occurring is thereupon selected as that for which the target function can be maximized. The method according to the invention for estimating the correction factor can therefore be performed both with the aid of received values which are assigned to pilot symbols, and with the aid of received values which are assigned to unknown data. The method according to the invention ensures enhancement of the quality of the parameter estimation, thus ensuring a reduction in the bit error rate of the transmission system and/or in a possible saving of necessary transmit power.

Claims: Method for reception in a CDMA receiver based on a parameter estimation method, in which the estimates are produced by processing received values (y i [k]) provided with correction factors (s), the correction factors (s) serving to equalize different power levels of the received signal, characterized in that a discontinuity of transmit power levels and the correction factors (s) are estimated by statistical hypothesis testing and are assigned to the received values (y i [k]). Method according to Claim 1, characterized in that a-priori probability (Pr (s)) for the presence of a specific discontinuity of a transmit power level being used in the statistical hypothesis testing to improve the estimation of the current discontinuity of the transmit power level of the received signal (y i [k]). Method according to one of the preceding claims, characterized in that the parameter estimation method undertakes to estimate the discontinuity of the transmit power between two blocks (n, n+1) respectively assigned to a transmit power. Method according to one of Claims 1 to 3, characterized in that in order to estimate the correction factor, the parameter estimation method uses a vector from component-by-component squares of absolute values of the received signal (y i [k]) to carry out the statistical hypothesis test. Method according to one of Claims 1 to 4, characterized in that sums of squares of absolute values are used to carry out the statistical hypothesis test. Method according to one of the preceding claims, characterized in that the statistical hypothesis test is carried out during the reception of received values (y i [k]) corresponding to training symbols. Method according to one of Claims 1 to 4, characterized in that only received values (y i [k]) from the block (n+1) belonging to the power control factor (s) currently to be estimated are used to estimate the power level. Method according to one of Claims 1 to 7, characterized in that received values (y i [k]) from the block (n+1) belonging to the power control factor (s) currently to be estimated, and from the preceding block (n) are used to estimate the power level. Method according to one of the preceding claims, characterized in that in the event of a discontinuity, the estimation of discontinuities of the power control factors (s) signals further estimation methods and is used as a basis for these further estimation methods. CDMA receiver comprising a parameter estimation device for carrying out a parameter estimation method in accordance with one of Claims 1 to 9. CDMA receiver according to Claim 10, characterized in that the receiver is designed as a rake receiver. CDMA receiver according to Claim 11, characterized in that more than one finger (F 0 , F i , F L-1 ) of the rake receiver are assigned a parameter estimation device (5). CDMA receiver according to Claim 11 or 12, characterized in that signals from more than one finger (F 0 , F i , F L-1 ) of the rake receiver are used to carry out the statistical hypothesis test. Implementation program for carrying out a method according to one of Claims 1 to 9, in particular inside a CDMA receiver according to one of Claims 10 to 13.

Inventors: BACHL RAINER, GERSTACKER WOLFGANG HELMUT, RAU RICHARD, BACHL, RAINER, GERSTACKER, WOLFGANG HELMUT, RAU, RICHARD
Applicants: LUCENT TECHNOLOGIES INC, LUCENT TECHNOLOGIES INC.
Application year: 2001
Main classification: H04B17
IPC codes: H04B1, H04B7, H04L25, H04J13, H04B17, H04W52