A Wiener filtering Approach to Subcarrier combining for Orthogonal Multi-Carrier CDMA

by Jean-Paul Linnartz

PIMRC 94

 

The application of Wiener filtering to the detection of Multi-Carrier Code Division Multiple Access (MC-CDMA) signals is analyzed. An analytical model is developed to evaluate an MC-CDMA system with a very large number of subcarriers. This model is compared with simulations for 8 and 64 subcarriers. Simulation results are also presented on the effect of correlated fading (in the frequency domain) at the subcarriers.

 

Background

With orthogonal MC-CDMA, each data symbol is simultaneously transmitted at N binary phase shift-keying (BPSK) narrowband subcarriers, each separated by F/Tb Hz where F is an integer.Due to its special signal structure, MC-CDMA signals will not experience significant linear distortion in fading channels where the symbol duration, , is much larger than the delay spread, . As shown in the transmitter model of Fig. 1, each of the N subcarrier waveforms is modulated (multiplied) by a single (time constant) chip belonging to a spreading code of length N. Different users transmit at the same set of subcarriers but with a different spreading code in the frequency domain. In contrast to Direct-Sequence CDMA. (or other wideband) signals, MC-CDMA signals are not significantly affected by delay spreads due to their narrowband composition. This modulation technique should not be confused with transmitting multiple DS-CDMA signals at different frequencies. Note that the signal structure of MC-CDMA is similar to that of Orthogonal Frequency Division Multiplexing (OFDM). but the manner in which the signals are used is very different.

While the general shape of each individual subcarrier will not experience significant linear distortion after transmission, the flat fading that each subcarrier with experience results in an amplitude scaling. Unequal scaling a different subcarriers distorts the orthogonality between users. While the conventional detection methods of equal gain combining (EGC) and maximal ratio combining (MRC) may be sufficient in an additive white Gaussian noise (AWGN) channel, these detection methods do not directly address the issues of orthogonality and of interference cancellation. Consequently, these detection methods do not perform as well in interference limited channels.

In this paper, we will apply Wiener filtering to the detection of MC-CDMA signals. Wiener filtering is optimal in a mean-squared error sense with respect to both the noise and the interference. In most fading channels, the determination of the Wiener coefficients is performed adaptively. In this paper, we are primarily concerned with the theoretical limits and not the implementation aspects of how to track the channel fluctuations. Thus, it is assumed that accurate estimates of the complete channel state information (i.e., the fading at the subcarriers) are available and that the Wiener coefficients are chosen accordingly.

Channel Model

Denote the vector of the data symbols by

where corresponds to the data symbol of the mth user. If user m is inactive, then . Define the code matrix C to be

where the mth column vector of C corresponds to the spreading code, , of the mth user. Using the vector and matrix notation described above, an equivalent discrete vector representation of the transmitted signal described in Fig. 1 is

where si[k] represents the signal components of all users at the ith subcarrier for the kth symbol interval. The actual continuous-time transmitted signal corresponding the kth data symbol of the mth user, am[k], is

where is the carrier frequency and is a unit pulse that is non-zero for and zero otherwise. Note that in a distortionless channel, that data symbol of the mth user can be recovered using the orthogonality of the codes as indicated by the following operation

where . denotes the inner product between two vectors and

is a vector containing the spreading code of the mth user.

In this paper, we will focus on downlink transmissions, i.e., from the base station to the mobiles. Assuming that the delay spread is much smaller than the symbol duration, the effect of the channel at the ith subcarrier may be approximated by a constant amplitude scaling, pi[k], and a constant phase offset, , over the symbol duration. Using this assumption, the index k will be suppressed throughout the rest of this paper. Applying the received signal to the receiver model shown in Fig. 2, the equivalent discrete representation of the received signal is

where represents the component of the received signal at the ith subcarrier, the channel matrix H is defined to be

,

and is a vector containing the corresponding AWGN terms with representing the noise term at the ith subcarrier with power . Obviously, the discrete representation of the demodulated signal is an approximation. However, under the condition , each subcarrier faces a relatively flat channel.

 

Linear Receivers

 

The application of Wiener filtering to this received signal involves linearly combining the different subcarrier diversity components to yield the decision variable

where the vector represents the optimal weighting coefficients. These coefficients can also be viewed as an amplitude equalization to compensate for the fading at the subcarriers. Arbitrarily choosing m = 0 as the desired signal, the optimal choice of the equalization vector D in the mean-squared sense can be determined to be

where

represents the autocorrelation matrix of the received vector Y and represents the cross-correlation vector between the desired symbol, a0, and the received signal vector, Y. In the formulation of this method, it is assumed that estimates of the channel amplitudes, , are available and are consequently treated as deterministic constants.

A closed form solution is difficult to obtain except for the case of a full load, i.e., when all N users are active. The corresponding equalization coefficient at the ith subcarrier for a full load is

.

As expected, this indicates that the linear combination of the subcarrier components should include the inner product of the desired user's spreading code vector with the received signal vector. For small , the equalization coefficient should be small to avoid excessive amplification of the noise. For large , the correction factor should be proportional to taking the inverse of the channel, , in order to restore the orthogonality between users.

For sufficiently large values of N, the average bit error rate (BER) for the case of a full load may be approximated using the Central Limit Theorem.

where

and

.

 

In obtaining this expression, independent and identically distributed (IID) fading at the subcarriers was assumed.

To study the effect of correlated fading through simulations, the following algorithm for generating correlated Rayleigh r.v.s was used. The correlation between two frequencies separated by radians in a complex gaussian fading channel is given in Jakes Jr., Microwave Mobile Communications, New York: Wiley, 1974. to be

.

where and is the local-mean power of the mth user. Defining as a vector of dimension N whose elements, hi for i = 0, 1, ..., N-1, are the complex correlated gaussian fading variables at the different subcarriers, the autocorrelation matrix of , where denotes the conjugate transpose of , can be determined. The vector can be generated by taking linear combinations of independent, idenctically distributed (i.i.d.) unit normal gaussian random variables, using

where G is a vector of dimension N containing i.i.d. unit normal r.v.s. A possible choice of the transformation matrix A is

where E is an matrix with the eigenvectors of as its column vectors and is a diagonal matrix with the eigenvalues of as its entries. The relationship between the fading variables of See , and is

.

Simulation and Numerical Results

Due to the lengthy time consumption of the simulations, results for only spreading factors of N = 8 and N = 64 were obtained. Shown in Fig. 3 is a plot of the average BER versus the signal-to-noise ratio (SNR) in a Rayleigh fading channel for a full load with N = 8. IID fading at the subcarriers is assumed. It is also assumed that all users are received with equal power. Note that the analytical results for Wiener filtering differ significantly from the simulation results for large SNR. This discrepancy is due to the application of the Central Limit Theorem to small spreading factors. Shown in Fig. 4 are the results for N = 64. Examining Fig. 3 and 4, it can be seen that See is more accurate for large N. Also included in Fig. 3 are the curves for a single narrowband subcarrier with AWGN, and with and without fading. Note that the performance of MC-CDMA is better than a single narrowband transmission due to frequency diversity even with a fully loaded system.

The average BER versus the SNR for correlated fading was also examined for the case of N = 8 (see Fig. 5). The degree of correlation was measured in terms of

.

Small a's correspond to highly correlated fading at the subcarriers while large a's correspond to IID fading. Ideally, the larger the value of a, the higher the degree of frequency diversity. In a DS-CDMA system, the number of resolvable paths is equal to

 

.

Note that for a = 0.125 (L = 2) most of the benefits in performance from frequency diversity have already been achieved.

Shown in Fig. 6 is a plot of the simulation results for the average BER versus the number of interferes. The SNR was chosen to be 10dB. Because of the time consuming nature of the simulator, a relatively small spreading factor of N = 8 was used. Besides the curves for Wiener filtering, the simulation results of conventional EGC and MRC detection are also included See N. Yee, J. P. Linnartz and G. Fettweis, "Multi-Carrier CDMA in Indoor Wireless Radio Networks," Proceedings PIMRC '93, Yokohama, Japan, 1993, pp. 109-113. . As expected, Wiener filtering outperforms EGC and MRC substantially in combating interference.

Conclusion

Wiener (or MMSE) filtering was applied to the detection process of MC-CDMA. In contrast to conventional diversity combining techniques, this detection technique directly addresses the effects of the interference on the BER. It was shown that even with a full load, MC-CDMA with Wiener filtering outperforms single narrowband transmissions in a Rayleigh fading channel. It was also shown that Wiener filtering outperforms conventional detection methods substantially. Analytical approximations for appeared reasonably accurate for SNRs up to 12dB and for N = 64 subcarriers. As the number of subcarriers is increased, the accuracy of the analytical results improves.. It was also shown that most of the benefits of frequency diversity can be achieved even in a channel with a significant degree of correlated fading at the subcarriers.

 

 

References

  1. N. Yee, J. P. Linnartz and G. Fettweis, "Multi-Carrier CDMA in Indoor Wireless Radio Networks," Proceedings PIMRC '93, Yokohama, Japan, 1993, pp. 109-113.
  2. G. Fettweis, "Multi-Carrier Code Division Multiple Access (MC-CDMA): Basic Idea," Technekron Communication Systems - Internal Report.
  3. R. L. Pickholtz, D. L. Schilling and L.B. Milstein, "Theory of Spread-Spectrum Communications - A Tutorial," IEEE Transactions on Communications, COM-30, May 1982, pp. 855-884.
  4. S. Kondo and L. B. Milstein, "MultiCarrier CDMA System with Co-channel Interference Cancellation for an Asynchronous Forward Link", 44th IEEE Veh Tech Conf, June 7-10 1994, Stockholm.
  5. Leonard J. Cimini, Jr., "Analysis and Simulation of a Digital Mobile Channel Using Orthogonal Frequency Division Multiplexing," IEEE Transactions on Communications, vol. COM-33, pp. 665-6.
  6. G. Fettweis, K. Anvari, and A. Shaikh Bahai, "On MultiCarrier Code Division Multiple Access (MC-CDMA) Modem Design", 44th IEEE Veh Tech Conf, June 7-11 1994, Stockholm.
  7. W. C. Jakes Jr., Microwave Mobile Communications, New York: Wiley, 1974.