## Noname manuscript No.

(will be inserted by the editor)

## Adaptive Blind CDMA Receivers Based on ICA Filtered Structures

Zaid Albataineh · Fathi M. Sale

the date of receipt and acceptance should be inserted later

Abstract Code Division Multiple Access (CDMA) is a channel access method adopted by various radio technologies world-wide. In particular, CDMA is used as an access method in many mobile standards such as CDMA2000 and WCDMA. We address the problem of blind multiuser equalization in the wideband CDMA systems in the noisy multipath propagation environment. Herein, we propose three new blind receiver schemes based on variations of Independent Component Analysis (ICA) within several filtering structures. These adaptive blind CDMA (ABC) receivers do not require knowledge of the propagation parameters or spreading code sequences of the users-they primarily exploit the natural assumption of statistical independence among the symbol signals. We also develop three semi-blind adaptive detectors by incorporating new adaptive methods into the standard Rake receiver structure. Extensive comparative case-studies, based on Bit error rate (BER) performance are carried out for as a function of (i) the number of users, (ii) the number of symbols per user, and (iii) the signal to noise ratio (SNR).The conventional detectors include the baseline Linear Minimum mean squared error (LMMSE) detector. The results show that the proposed methods outperform other detectors in estimating the symbol signals from the received mixed CDMA signals. Moreover, the new blind detectors mitigate the multi access interference (MAI) in CDMA.

Keywords Direct Sequence Code Division Multiple Access (DS-CDMA) systems · Wide-band CDMA (W-CDMA) · Independent Component Analysis (ICA) · Robust ICA · Linear Minimum mean squared error (LMMSE) · Multi Access Interference (MAI) · Rake detector · Principle Component Anaylisis (PCA) · FAST ICA · Bit error rate (BER) · signal to noise ratio (SNR).

Address(es) of author(s) should be given

## 1 Introduction

Code Division Multiple Access (CDMA) is a channel access method ubiquitously used in various modalities and platforms worldwide. It is based on spread-spectrum technology as is found, e.g., in third-generation (3G) cellular telephony, terrestrial and satellite communications systems, and indoor wireless networks [1-2], [9]. Although, LTE (4G) is utilized by several cellular companies inside and outside the U.S., their networks are still not fully built, and LTE coverage is still not universal. Thus, most of the older 2G and 3G systems are ubiquitous and exist in parallel with the newer 4G systems worldwide. In the U.S., companies like AT&amp;T and T-Mobile use GSM/WCDMA/HSPA while Verizon, Sprint, and MetroPCS use CDMA2000/EV-DO [3-5]. Moreover, the newer LTE wireless interface is incompatible with the 2G and 3G networks, so that it must be operated on a separate wireless spectrum. While 4G technology is intended to eventually replace the 3G technologies, it is now evident that it will take some time before LTE coverage is fully developed and widely adopted even in the developed countries [26-27].

As with any radio communication system, CDMA-based systems also suffer from various types of interferences. Specifically, they suffer from (i) an internal multiple access interference (MAI) due to the non-ideal cross-correlations among the users spreading sequences, (ii) narrow-band inter-symbol interference (ISI), and (iii) background noise at the receiver. These drawbacks affect the performance of a CDMA system. The conventional detectors most frequently utilized to counteract CDMA interference is based on second order statistics. In highly loaded systems, conventional detectors are not considered a suitable choice. Most of the conventional detectors suffer from external interference sources and treat all interferences as a lumped background noise. In CDMA-based systems, however, the primary source of interference is MAI. This has motivated the development of numerous interference rejection techniques to overcome the MAI and the near-far problem in conventional receivers [1, 7]. Several state-of-the-art approaches have been proposed in the literature to overcome this challenge, e.g., using pilot signals and training [40].

In CDMA-based systems, multiuser detection is desirable in order to enhance channel capacity and mitigate MAI [10, 21]. Multiuser detection has been introduced to obtain an optimum multiuser detector for multi-Gaussian channels in [1]. Several suboptimal detectors have also been proposed in [6-8], to overcome the computational complexity in realizing optimal detectors. In [1] and [32-36], training pilot sequence techniques have been used to present suboptimal detectors, namely an adaptive linear detector and a zero-forcing detector.

In [20, 21], Wang and Poor proposed the blind minimum mean square error (MMSE) and the blind de-correlating detectors. The suboptimal detector based on the linear minimum mean square error (LMMSE) method has been described in [32].In [31-36], adaptive blind detectors were proposed based on incorporating the minimum output energy with constrained optimization methods. Several subspace approaches were proposed in the literature, e.g., in [20],

[23], [36]. In [10], several types of group-blind linear detectors were proposed in order to enhance the performance for the uplink and downlink channels. The key idea of these detectors is to take advantage of the cross-correlation matrix which was constructed by exploiting the correlation between successive samples of received signals. These detectors, however, are too complex to be practically implemented, especially at the mobile unit. Also, they require information regarding signal timing and the spreading codes of the desired user.

The aforementioned techniques periodically require the base station to send a training sequence that must be known by the mobile receiver in order to enable the latter in estimating the parameters of the channel propagation model. These parameters also attempt to capture the multiple reflections of the radio waves on encountered obstacles, i.e. buildings, cars, trees, etc. Furthermore, according to [42], it has been reported that 20% of the bandwidth in GSM, and up to 40% in UMTS CDMA, is devoted to the training sequence. In spite of the good performance of the training sequence techniques, the cost tends to be significantly large in terms of bandwidth. Adaptive signal processing techniques, on the other hand, provide more efficient methods for CDMA systems in the presence of high dynamic conditions as a result of the receiver mobility, the short channel codes and the fortuitous channel access. In particular, the desire to ensure a high communication rate has made blind adaptive techniques a hot topic, driven by their potential to eliminate/reduce training sessions. Moreover, blind techniques help recover symbol signals in other situations e.g., i) eavesdropping, where using the training sequence is not available, and ii) tracking, when the receiver fails to keep the desired user locked in track. It is also noted that the underlying user symbol sequences are reasonably assumed to be statistically independent. Therefore, statistical independence, or near independence, is a key assumption that makes a CDMA system suitable for the blind techniques, e.g., using information maximization [1] or minimum mutual information [6]. In [6-8], typical CDMA based systems are represented by wide stationary slowly fading multipath environment and are expressed by a linear multi-channel convolution model. Thus, the received signals in a CDMA mobile can be considered as signals generated by the linear convolutive model of statistically independent components of independent users as shown in [6], [10], [31-36]. The adaptive LMMSE detector has been originally proposed to overcome the necessary complex matrix inversion operation [38]; however, it still requires the spreading codes of all users. While the LMMSE detector maybe suitable for the uplink to the base station, as computational resources are usually abundant, it is less practical in the downlink to the receiver as computational resources are scarcer.

This paper aims at recovering the source symbol sequences from the linear convolutive received mixture without any knowledge of the user short channelizing codes and in the absence of explicit channel identification. In essence, the paper proposes improved blind adaptive detections, based on the state space approach[37, 6],using the natural gradient method for multipath channels of CDMA-based systems. Three update-laws are derived for various

filtering structures [6] and then three adaptive blind CDMA detectors are introduced for more effective MAI, ISI suppression and symbol estimation. The second contribution of the paper is three semi-blind adaptive stochastic gradient algorithms fused into the conventional Rake receiver. Specifically, we fuse algorithms based on, respectively, FastICA, RobustICA, and principle component analysis (PCA). Furthermore, higher order statistics (HOS) are exploited in order to make the proposed methods robust and secure against incomplete cross-correlation and the near-far problem in conventional detectors [42]. Extensive Monte Carlo simulations have been carried out to verify and evaluate the effectiveness of the proposed methods in estimating the users symbols. In summary, we provide metric comparisons in the bit error-rate (BER) as a function of (i) the number of users, and (ii) the number of symbols per user, and (iii) the signal-to-noise (SNR). The comparisons include the proposed methods with existing and conventional ones in terms of BER performance and computational complexity.

We now set the notation used throughout the paper. Lower case letters denote scalars, bold lower case letters denote vectors, and bold upper case letters denote matrices. Moreover, the following symbols are used:

- * ( ) . T refers to the transpose operator;
- * ( ) . H refers to the Hermitian transpose operator;
- * trace ( ) refers to the trace operator; .
- * j = √ -1 refers to the imaginary symbol;
- * diag ( ) refers to the standard diagonal of a matrix; .
- * Diag ( ) . refers to the diagonal of a block matrix, where elements may be block matrices themselves;
- * sgn ( ) refers to the sign operator; .
- * E . [ ] refers to the statistical expectation operator.

The remainder of the paper is organized as follows. In Section II, brief descriptions and derivations of synchronous CDMA signal models in multipath fading are presented. The conventional Rake receiver model is described in Section III. Section IV is dedicated to the derivation of adaptive update laws and to the proposed new detection schemes. The comparative simulations with summary results and conclusions are given in Section V and Section VI, respectively.

## 2 CDMA SIGNAL MODEL

We now briefly present two signal models for a CDMA based system using one layer of channel spreading codes. Specifically, we describe the DS-CDMA signal and WCDMA signal models in a typical synchronous CDMA system usually employed, e.g., for cellphones, indoor ATM, and certain ad hoc wireless networks [1], [3].

## 2.1 A DS-CDMA Receiver Signal Model

In a DS-CDMA system, several users share the medium simultaneously by using unique individualized code signatures. We refer to Fig. 1 below for a typical system schematic block diagram. In this paper, we assume the data transmission to be quaternary phase shift keying (QPSK). At the mobile unit receiver, assume a total of K active users in an L multipath environment and M transmitted symbols during the observation frame time. The simplest downlink received signal model r t ( ), at the time sample t over a single symbol interval is given by [21]

$$r \left ( t \right ) = \sum _ { m = 1 } ^ { M } \sum _ { k = 1 } ^ { K } \sum _ { l = 0 } ^ { L - 1 } \alpha _ { l m } b _ { k, m } s _ { k } \left ( t - m T _ { b } - d _ { l } T _ { c } \right ) + n \left ( t \right ) \text{ \quad \ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text$$

where

- -l,k,m are path, user and symbol indices, respectively.
- -α lm is the path gain - in the downlink model the path gain is assumed to be the same among users because all users' signals are transmitted together. Thus, the path gain α lm and propagation delay factor d l do not depend on the user k.
- -b k,m is the kth user m symbol.
- -s k ( ) is the kth user spreading code (chip sequence). .
- -d l is the propagation delay factor.
- -t, T b , T c are time, symbol, and chip duration, respectively.
- -n t ( ) is the channel additive white Gaussian noise (AWGN) with zero mean and covariance equals q.

The system is assumed to be time-invariant, over a small duration, which means that the channel parameters are much slower than the frequency of transmitted symbol data. Let us assume that G is the number of chips per symbol, K is the number of users, and L is the number of paths. Thus, the scalar form of Equation (1) can be transformed to a vector form [6, 21] as:

$$r = H S b + n$$

where r is a received (G1) dimensional vector signal; H is a (G1) x G matrix with G1 ≥ G+L -1, which represents the multipath propagation coefficients; S is a G x K block-diagonal matrix; b is a K dimensional vector, which represents the users data symbols; and n is the ( G 1) dimensional channel noise vector with covariance matrix, say, Q . This standardized model of received signals has been used in deriving the conventional detectors, e.g., Match filter, Rake filter, blind LMMSE and other blind detectors [21]. We shall use it in our development as well. In addition, an alternative two-tap symbol signal model is given by [10]:

$$r _ { n } = H _ { 0 } b _ { n } + H _ { 1 } b _ { n - 1 } + n _ { n } = \bar { H } \bar { b } _ { n } + n _ { n }$$

where

Fig. 1 Signal generation model for a typical QPSK DS-CDMA system

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-3ecfe392-475e-4c63-9393-1fb11b111042.jpg)

- * r n is the total received user's signal vector;
- * H 0 = [ h 1 , . . . , h k ] is the signature matrix of the current symbol vectors of all users including MAI, specifically,

$$\max \pi, \, \text{perm} ( \Omega ), \\ h _ { k } = \begin{bmatrix} 0 \\ h _ { k } ( 0 ) \\ \cdot \\ \cdot \\ h _ { k } ( G - D _ { l } - 1 ) \end{bmatrix} \quad \begin{bmatrix} 4 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 1 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ \cdot \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 3 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 2 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 6 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 5 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 9 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ } 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 } \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 4 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } & } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } $ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }; } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }} } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \\ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }$$

- * H 1 = [ ¯ h 1 , . . . , ¯ h k ] is the signature matrix of the previous symbol vectors of all users including ISI, where,

$$\mathbb { m } _ { 5 } \omega, \, \text{max}, \\ \mathring { h } _ { k } = \begin{bmatrix} \ h _ { k } \left ( G - D _ { l } \right ) \\. \\. \\ \hbar { \, } \left ( G + L - 1 \right ) \\ 0 \end{bmatrix} \quad & \quad ( 5 ) \\ \bullet \, \cdot \, \cdot \,. \,. \,. \,. \,. \,. \,.$$

D l ∈ { 0 1 , , . . . , G - } 1 is the delay in chip periods.

- * ¯ H = [ H H 0 1 ] is the signature matrix of all users;
- * b n = [ b 1 ( n , . . . , ) b K ( n )] T are the current symbols of all users;
- * b n -1 = [ b 1 ( n -1) , . . . , b K ( n -1)] T are the previous symbols of all users;
- * ¯ b n = [ b T n , b T n -1 ] T are the augmented two-tap symbols of all users;
- * n n = [ n ( nG ,..., ) n ( nG + G -1)] T is the independent white composite Gaussian Noise vector. We defer further details to [6, 10].

In the asynchronous uplink CDMA systems, one can assume that the columns of H 0 and H 1 are mutually independent. Therefore, ¯ H is a full rank matrix. Whereas for the synchronous downlink CDMA communication, ¯ H is full-rank with some restrictions. The main focus in this paper is on the synchronous downlink CDMA communication system, although our proposed algorithms work well in the uplink asynchronous CDMA systems [10], [30].

## 2.2 WCDMA Receiver Signal Model

One difference between a WCDMA system and a DS-CDMA system is the presence of scrambling codes. The main cause of the MAI in WCDMA systems is the intra-cell multiple user signals sharing the same multipath channels. Fig. 2 depicts a block diagram that shows the additional code scrambling before transmission through the air interface. In Fig. 2, DPDCH stands for Dedicated Physical Data CHannel which is a term adopted in UMTS (Universal Mobile Telecommunications Systems) and a S/P block stands for serial to parallel Converter. Consequently, the basic received signal model r(t) is given by [6]:

$$r \left ( t \right ) = \sum _ { m = 1 } ^ { M } \sum _ { k = 1 } ^ { K } \sum _ { l = 0 } ^ { L } \alpha _ { l m } b _ { k, m } c _ { k } \left ( t - \ d _ { l } T _ { c } \right ) s _ { k } \left ( t - m T _ { b } - d _ { l } T _ { c } \right ) + n \left ( t \right ) \quad ( 6 )$$

where, in addition to the previous parameters, one adds c k (t) ∈ {± 1 ± } j , the complex cell-specific scrambling sequences. The remaining variables are defined in model (1). The received signal at the mobile unit is passed through a chip-matched filter and sampled at the chip rate. The received discrete vector r in this case can be expressed as [6, 10, 21].

$$r = \text{HCSb} + n$$

where C is the GxG complex diagonal scrambling matrix with CC H = I GxG and the remaining variables are defined similarly as in (2). The form of C is given by:

$$C \, = \, \text{diag} ( \mathbf c _ { 1 } \, \ \mathbf c _ { 2 } \dots \ \ \mathbf c _ { G } \, ) \quad \quad \quad ( 8 ) \\ \text{where} \, \mathbf c _ { i } \in \{ \pm 1 \, \pm j \} \quad \forall \ 1 \leq i \, \leq G$$

## 3 CONVENTIONAL BLIND LINEAR MULTIUSER DETECTORS

We briefly describe the baseline conventional linear multiuser detectors such as the Match Filter (MF), the Rake receiver and the LMMSE detector in multipath environments. For further details, see [1, 9]. '

Fig. 2 Signal generation based on the proposed 3GPP UMTS FDD standard

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-2a1c2dbc-c4f3-4ed0-8c26-66fb339c7131.jpg)

## 3.1 Single user detector (SUD)

The SUD is a standard MF detector which exploits the users code signature to provide an estimate of the users symbol sequence from the received data. This detector completely ignores the presence of MAI due to other users. One can express the MF detector for the ith user in the DS-CDMA system as follows:

$$b _ { i, M F } ^ { D } = S _ { i } ^ { H } r$$

where S i = Diag ( ¯ s i , ¯ s i , . . . , ¯ s i ) , ¯ s i = [0 0 . . . s i . . . 0]. s i is the ith users signature code, r is the received discrete signal vector, and b D i , MF is the estimated DS-CDMA ith symbol vector.

## 3.2 Rake Detector

Perhaps, the most popular linear user detection is the Rake detector, which consists of multiple parallel chip-delayed SUD fingers. In this paper, we implement the Rake detector with the estimated known channel gain coefficients, but not the channel delays. One can express the Rake detector for the DSCDMA system mathematically as follows:

$$b _ { i, \text{Rake} } ^ { D } = S _ { i } ^ { H } H ^ { H } r$$

where H represents the estimated channel matrix, and b D Rake i , is the estimated ith users symbol vector.

## 3.3 LMMSE Detector

Conventional linear detectors based on the Least Square (LS), Zero-Force (ZF) and BLUE algorithms [9] perform poorly especially in the presence of colored noise. The LMMSE detector, however, is considered to be one of the best linear detectors for DS-CDMA systems. Mathematically, one can express the LMMSE as follows:

$$b _ { i, L M M S E } ^ { D } = S _ { i } ^ { H } \mathbf H ^ { H } \left ( \sigma ^ { 2 } \mathbf H \mathbf H ^ { H } + Q \right ) ^ { - 1 } r$$

where ( σ 2 HH H + Q ) = R = E [ rr H ] is the auto-correlation of the received data at the mobile unit, and σ 2 is the average power of the received signal. There are several drawbacks in the implementation of the LMMSE receiver. The main drawback is that the computation of the auto-correlation R is very expensive. If possible, one may use eigen-structure decomposition instead of inverting the auto-correlation matrix R directly to obtain

$$b _ { i, L M M S E } ^ { W } = S _ { i } ^ { H } H ^ { H } \left ( V _ { s } D _ { s } ^ { - 1 } V _ { s } ^ { H } \right ) r$$

where V s is the estimated eigen-vector matrix of the auto-correlation matrix R , and D s is the corresponding diagonal eigenvalue matrix. Additionally, one can use adaptive algorithms to estimate the LMMSE users symbols as in [32].

## 4 THE PROPOSED ADAPTIVE BLIND DETECTION SCHEMES

In this section, we introduce new blind detection strategies for the filtering structures. We propose three blind multiuser detectors based on (i) a feedforward structure, (ii) a feedback structure I, and (iii) a feedback structure II, as in [6]. These filtering structures are depicted in figures 3, 4 and 5, respectively.

To that end, one recalls the discrete received signal model (3), namely,

$$r _ { n } = H _ { 0 } b _ { n } + H _ { 1 } b _ { n - 1 } + n _ { n }$$

The aim here is to detect the symbol vector b n from the received data vector r n , over the discrete index n, under the following assumptions:

- * AS1) the G1xK matrices H 0 and H 1 are of full column rank.
- * AS2) the symbol signal vector series, b n , have statistically independent components and are identically distributed (i.i.d).
- * AS3) the Additive Noise vector n n is white, Gaussian, and independent of the symbol source signals.
- * AS4) the power of the transmitted symbol signals are normalized to be unity.
- * AS5) the maximum lag in the entire multipath channels is smaller than the spreading gain G of the CDMA codes.

- * AS6) the CDMA system is not over-saturated, which means the number of users (K) is less than the number of the spreading gain (G).
- * AS7) the channel is assumed to be a slowly fading wide sense stationary.

For methodical convenience, each detector algorithm involves two steps: first, a preprocessing stage; second, the (matrix) rotation stage based on the filtering structures. In the next subsection, we will present the common preprocessing stage (i.e., whitening processes), and then we will derive each of the three algorithms based on each filtering structure in individual subsections.

## 4.1 Step1: Preprocessing (i.e. Data Whitening)

The outcome of this step is that the symbol signals are detected up to a unitary rotational matrix. This step uses second order statistics (SOS) in order to normalize the variance (or power) of the received discrete signal vector. It may also be used to eliminate redundancy in the data based on PCA. Under Assumptions AS1-AS4, the G1xG1 covariance matrix, say ( Cov ), of the noiseless received discrete signal vector can be expressed as

$$C o v = E \left [ r _ { n } r _ { n } ^ { H } \right ] - q I _ { G 1 }$$

We will now consider the two-tap signal model. Then we may generalize it using induction techniques. Under Assumptions AS1-AS7, substituting r n from eqn (3) into (13) results in the following covariance matrix:

$$\begin{smallmatrix} \text{Cov} = \text{H} _ { 0 } E \left [ b _ { n } b _ { n } ^ { H } \right ] \text{H} _ { 0 } ^ { H } + \text{H} _ { 1 } E \left [ b _ { n - 1 } b _ { n - 1 } ^ { H } \right ] \text{H} _ { 1 } ^ { H } \\ \text{Cov} = \text{H} _ { 0 } \text{H} _ { 0 } ^ { H } + \text{H} _ { 1 } \text{H} _ { 1 } ^ { H } = \left [ \text{H} _ { 0 } \text{ } \text{H} _ { 1 } \right ] \text{[H} _ { 0 } \text{ } \text{H} _ { 1 } \right ] ^ { H } \end{smallmatrix}$$

Observe that under AS2, E [ b b n H n ] = I K and [ b n -1 b H n -1 ] = I K . Without loss of generality, we shall briefly proceed with the basic algebraic procedure by adopting the eigen-structure decomposition for the symmetric square matrix Cov and use it to obtain a singular value decomposition for the combined matrix [ H 0 H 1 ] Thus, let

$$C o v = V D V ^ { H }$$

where V is a G1xG1 matrix of orthogonal eigenvectors satisfying

$$V V ^ { H } = V ^ { H } V = I _ { G 1 }$$

and D is the corresponding G1xG1 diagonal eigen-matrix containing its eigenvalue entries along the diagonal. Thus, from (14), the GxG H 0 and H 1 matrices can be represented respectively as

$$\begin{smallmatrix} { \mathbf H } _ { 0 } = { \mathbf V } _ { 0 } { \Lambda } _ { 0 } { \mathbf U } _ { 0 } ^ { \mathbf H } \\ { \mathbf H } _ { 1 } = { \mathbf V } _ { 1 } { \Lambda } _ { 1 } { \mathbf U } _ { 1 } ^ { \mathbf H } \end{smallmatrix}$$

where V 0 and V 1 are composed of orderly non-overlapping columns of the G1xG1 unitary matrix V U . 0 and U 1 are constant but unknown KxK right

singular-value unitary matrices with U U 0 0 H = U U 1 1 H = I K , and Λ 0 and Λ 1 are the appropriate GxK singular value matrices. We note that the whitening or algebraic PCA procedure can (i) estimate the noise power in eqn(13) and (ii) reduce the whitened signal dimension to the signal subspace, in this case K. Now, we process the received data to obtain the (whitened) data, specifically, we define:

$$r _ { n } ^ { w } = \Lambda ^ { + } \mathbf V ^ { H } r _ { n }$$

where the KxG1 matrix Λ + denotes the pseudo-inverse of the singular value matrices. One simplifies (18) to eventually obtain:

$$r _ { n } ^ { w } = U _ { 0 } ^ { \, H } b _ { n } + U _ { 1 } ^ { \, H } b _ { n - 1 } + \left ( \Lambda ^ { + } V ^ { H } \right ) n _ { n }$$

Thus, the whitening step renders the whitened data expressed in (18) or (19) as having a reduced dimension to the symbol space and a covariance matrix equal to the identity. That is E [ r w wH n r n ] = I K .

Note that, after the preprocessing step, the detection of the symbol signal ˆ b n reduces to determining or compensating for the unknown K x K (rotation) unitary matrices U 0 and U 1 . Next, we proceed with the development and derivations of the three proposed adaptive filtering structures, based on (i) feed-forward structure (FF), (ii) feedback structure I (FB-I) and (iii) feedback structure II (FB-II) [6].

Remark: For the purposes of the adaptive filtering to be discussed next, we shall re-label these unknown (but fixed) unitary matrices as the starred values for the environment. Specifically, in eqn(19), we set

∗

$$U _ { 0 } = U _ { 0 } ^ { * } \\ U _ { 1 } = U _ { 1 } ^ { * }$$

The developed adaptive filtering will have parameter matrices that, when adaptation is successful, will converge to (approximately) these fixed starred environment parameters.

4.2 Step 2a: Determining the rotation unitary matrix U for the feedforward structure

The output from the FF structure, as depicted in Fig. 3, is expressed as

$$y _ { n } = U _ { 0 } r _ { n } ^ { w } + \sum _ { k = 1 } ^ { K } U _ { k } r _ { n - k } ^ { w }$$

For simplicity of presentation, we begin with a two-tap model; thus the two-tap FF structure becomes

$$y _ { n } = U _ { 0 } r _ { n } ^ { w } + U _ { 1 } r _ { n - 1 } ^ { w }$$

$$\Phi _ { 0 } r _ { n } ^ { w } + U _ { 1 } r _ { n - 1 } ^ { w }$$

Fig. 3 Feed-Forward (FF) Demixing Structure

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-aee2d3db-cb81-4e48-91d4-3ff23dfd2742.jpg)

The goal for a successful adaptive algorithm is to bring about the convergence of the parameter matrices to the (starred) environment parameters. Specially, the adaptive algorithm succeeds when its parameter matrices converge to U ∗ 0 , and U ∗ 1 respectively.

We now proceed with the development. One can re-write this convolutive filter (21) as the following (static) map

$$\left [ \begin{smallmatrix} y _ { n } \\ r _ { n - 1 } ^ { w } \end{smallmatrix} \right ] = \left [ \begin{smallmatrix} U _ { 0 } & & U _ { 1 } \\ 0 & & I \end{smallmatrix} \right ] \left [ \begin{smallmatrix} r _ { n } ^ { w } \\ r _ { n - 1 } ^ { w } \end{smallmatrix} \right ]$$

Then, one defines, respectively, the new augmented output, the static map, and the augmented input as

$$\tilde { Y } & = \begin{bmatrix} y _ { n } \\ r _ { n - 1 } ^ { w } \end{bmatrix} \\ \tilde { U } & = \begin{bmatrix} U _ { 0 } & 0 \\ U _ { 1 } & 1 \end{bmatrix} \\ \tilde { R } & = \begin{bmatrix} r _ { n } ^ { w } \\ r _ { n - 1 } ^ { w } \end{bmatrix} \\ \? 2 \end{bmatrix} \text{becomes the static ma}$$

Thus, the expression in (22) becomes the static map

$$\tilde { \mathbf Y } = \tilde { U } ^ { H } \tilde { \mathbf R }$$

Based on the natural gradient approach [43, 44], the update law for the columns of the augmented de-mixing matrix ˜ U can be expressed as

$$u ^ { + } = u - \mu E \left [ \tilde { R } \left ( g \left ( u ^ { H } \tilde { R } \right ) \right ) \right ]$$

where u , respectively u + , is the current, respectively next, value of one column vector of ˜ , U µ is the step size and g is the chosen score function. Noting the structure of the de-mixing matrix in (23), one decomposes the column vector as

$$u = \begin{bmatrix} u _ { 0 } \\ u _ { 1 } \end{bmatrix}$$

Hence, the update law is correspondingly decomposed as (Note, we have suppressed the E . [ ] operator):

$$\begin{bmatrix} u _ { 0 } ^ { + } \\ u _ { 1 } ^ { + } \end{bmatrix} = \begin{bmatrix} u _ { 0 } \\ u _ { 1 } \end{bmatrix} - \mu \, \begin{bmatrix} r _ { n } ^ { w } \\ r _ { n - 1 } ^ { w } \end{bmatrix} g \left ( y _ { n } \right )$$

where u 0 , u 1 are the column vectors of U 0 and U 1 in (22), respectively. Therefore, the update laws for the individual (sub-) columns are

$$u _ { 0 } ^ { + } = u _ { 0 } - \mu r _ { n } ^ { w } g \left ( y _ { n } \right )$$

$$u _ { 1 } ^ { + } = u _ { 1 } - \mu r _ { n - 1 } ^ { w } g \left ( y _ { n } \right )$$

Now, by induction, the update law for the kth lag element u k is

$$u _ { k } ^ { + } = u _ { k } - \mu r _ { n - k } ^ { w } g \left ( y _ { n } \right )$$

4.3 Step 2b: Determining the rotation unitary matrix U based on feedback structure I (FB-I)

The output of FB-I, as depicted in Fig. 4, results in the filtering expression

$$y _ { n } = U _ { 0 } ^ { - 1 } \left ( r _ { n } ^ { w } - \sum _ { k = 1 } ^ { K } U _ { k } y _ { n - k } \right )$$

Consider now just two taps of FB-I , i.e.,

$$y _ { n } = U _ { 0 } ^ { - 1 } \left ( r _ { n } ^ { w } - U _ { 1 } y _ { n - 1 } \right )$$

One can re-write this convolutive filter into the following augmented static form

$$\begin{bmatrix} r _ { n } ^ { w } \\ y _ { n - 1 } \end{bmatrix} = \begin{bmatrix} U _ { 0 } & & U _ { 1 } \\ 0 & & I \end{bmatrix} \begin{bmatrix} y _ { n } \\ y _ { n - 1 } \end{bmatrix}$$

Or

$$\begin{bmatrix} y _ { n } \\ y _ { n - 1 } \end{bmatrix} = \begin{bmatrix} U _ { 0 } & & U _ { 1 } \\ 0 & & I \end{bmatrix} ^ { - 1 } \begin{bmatrix} r _ { n } ^ { w } \\ y _ { n - 1 } \end{bmatrix}$$

Fig. 4 Feedback Demixing Structure I FB-I

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-1ce54ab7-ca8b-4ef4-8c00-8d9a7fdb2e79.jpg)

$$\begin{bmatrix} y _ { n } \\ y _ { n - 1 } \end{bmatrix} = \begin{bmatrix} U _ { 0 } ^ { - 1 } & - U _ { 0 } ^ { - 1 } U _ { 1 } \\ 0 & I \end{bmatrix} \begin{bmatrix} r _ { n } ^ { w } \\ y _ { n - 1 } \end{bmatrix}$$

Thus, in this case, one defines the augmented output, de-mixing matrix, and input as follows:

$$\tilde { Y } & = \begin{bmatrix} y _ { n } \\ y _ { n - 1 } \end{bmatrix} \\ \tilde { U } & = \begin{bmatrix} U _ { 0 } & 0 \\ U _ { 1 } & I \end{bmatrix} \\ \tilde { R } & = \begin{bmatrix} r _ { n } ^ { w } \\ y _ { n - 1 } \end{bmatrix} \\ \} ) \, \text{into the compact equatic}$$

One then re-expresses (32) into the compact equation

$$\tilde { \mathbf R } = \tilde { U } ^ { H } \tilde { \mathbf Y }$$

Again, using the natural gradient approach, the update law for a column of the de-mixing matrix˜ U is

$$u ^ { + } = u - \mu E \left [ \tilde { Y } \left ( g \left ( u ^ { H } \tilde { Y } \right ) \right ) \right ]$$

As before, u , respectively u + , is the current, respectively next, value of one column vectors of ˜ , U µ is the step size and g is the chosen score function.

One can exploit the block matrix structure of the de-mixing matrix and simplify the update law. To that end, consider the block matrix

$$u ^ { 0 } = \begin{bmatrix} u _ { 0 } ^ { 0 } \\ u _ { 1 } ^ { 0 } \end{bmatrix}$$

Thus, the update laws can be calculated to produce

$$\begin{bmatrix} u _ { 0 } ^ { + } \\ u _ { 1 } ^ { + } \end{bmatrix} = \begin{bmatrix} u _ { 0 } \\ u _ { 1 } \end{bmatrix} - \ \mu \begin{bmatrix} y _ { n } \\ y _ { n - 1 } \end{bmatrix} \, g \left ( u _ { 0 } y _ { n } + u _ { 1 } y _ { n - 1 } \right ) \quad \quad ( 3 7 )$$

Similarly, the next block matrices can be defined as

$$u ^ { 1 } = \left [ \begin{smallmatrix} 0 \\ i ^ { 1 } \end{smallmatrix} \right ]$$

This leads to the specialized form

$$\left [ \begin{matrix} 0 ^ { + } \\ i ^ { + } \end{matrix} \right ] = \left [ \begin{matrix} 0 \\ i \end{matrix} \right ] - \mu \, \left [ \begin{matrix} y _ { n } \\ y _ { n - 1 } \end{matrix} \right ] g \left ( y _ { n - 1 } \right )$$

Thus, the update laws for the individual columns are

$$u _ { 0 } ^ { + } = u _ { 0 } - \mu y _ { n } g \left ( r _ { n } ^ { w } \right )$$

and

$$u _ { 1 } ^ { + } = u _ { 1 } - \mu y _ { n - 1 } g \left ( r _ { n } ^ { w } \right )$$

Analogously, by induction, the update law for the kth lag element, say u k is

$$u _ { k } ^ { + } = u _ { k } - \mu y _ { n - k } g \left ( r _ { n } ^ { w } \right )$$

4.4 Step 2c: Determining the rotation unitary matrix U based on Feedback Structure II

The output of FB-II, as depicted in Fig. 5, is expressed as:

$$y _ { n } = U _ { 0 } r _ { n } ^ { w } + \sum _ { k = 1 } ^ { K } U _ { k } y _ { n - k }$$

Again, consider two taps of FB-II, i.e.,

$$y _ { n } = U _ { 0 } r _ { n } ^ { w } - U _ { 1 } y _ { n - 1 } & & ( 4 4 )$$

Hence, one re-writes this convolutive filter in the following augmented static form

$$\begin{bmatrix} y _ { n } \\ y _ { n - 1 } \end{bmatrix} = \begin{bmatrix} U _ { 0 } & & - U _ { 1 } \\ 0 & & I \end{bmatrix} \begin{bmatrix} r _ { n } ^ { w } \\ y _ { n - 1 } \end{bmatrix}$$

Similarly, define the augmented entities as

$$\tilde { Y } = \begin{bmatrix} y _ { n } \\ y _ { n - 1 } \end{bmatrix}$$

Fig. 5 Feedback Demixing Structure II (FB-II)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-5192a62d-0c58-4240-a59e-24fa94e39aaa.jpg)

$$\tilde { U } & = \begin{bmatrix} U _ { 0 } & 0 \\ - U _ { 1 } & I \end{bmatrix} \\ \tilde { R } & = \begin{bmatrix} r _ { n } ^ { w } \\ y _ { n - 1 } \end{bmatrix} \\. &. &.$$

Thus one re-writes (45) into the compact mapping

$$\tilde { \mathbf Y } = \tilde { U } ^ { H } \tilde { \mathbf R }$$

Using the natural gradient approach, the update laws for a weight column of the de-mixing matrix ˜ U is expressed as

$$\mathfrak { u } ^ { + } = \mathfrak { u } - \mu \mathbb { E } \left [ \tilde { R } \left ( g \left ( \mathfrak { u } ^ { H } \tilde { R } \right ) \right ) \right ]$$

where, as before, u , respectively u + , is the current, respectively next, value of one column vectors of ˜ , U µ is the step size and g ( ) is the chosen score function. . One can appropriately decompose a column vector in order to simplify the update expressions as:

$$u = \begin{bmatrix} u _ { 0 } \\ u _ { 1 } \end{bmatrix}$$

Then the update law becomes decomposed as follows:

$$\left [ \begin{matrix} u _ { 0 } ^ { + } \\ u _ { 1 } ^ { + } \end{matrix} \right ] = \left [ \begin{matrix} u _ { 0 } \\ u _ { 1 } \end{matrix} \right ] - \mu \, \left [ \begin{matrix} r _ { n } ^ { w } \\ y _ { n - 1 } \end{matrix} \right ] g \left ( y _ { n } \right )$$

Thus, the update laws for the individual sub-columns are

$$u _ { 0 } ^ { + } = u _ { 0 } - \mu r _ { n } ^ { w } g ( y _ { n } )$$

and

$$u _ { 1 } ^ { + } = u _ { 1 } - \mu y _ { n - 1 } g \left ( y _ { n } \right )$$

Finally, by induction, the update law for the kth lag element u k is

$$u _ { k } ^ { + } = u _ { k } - \mu y _ { n - k } g \left ( y _ { n } \right )$$

4.5 The proposed adaptive Rake-based detectors

While the previous filtering structures constitute new adaptive filters, one can augment the existing conventional Rake detectors to improve its performance adaptively. We now develop three adaptive modifications of the conventional Rake detector based on, respectively, Independent Component Analysis (ICA), Robust ICA and Principle Component Analysis (PCA). Recalling the Rake detector's structure as given in (10), one can mathematically express the adaptive modified Rake detector for DS-CDMA systems as follows:

$$b _ { i, \text{Rake} } ^ { D } = S _ { i } ^ { H } W H ^ { H } r$$

where, as before, H is the crudely estimated (inverse) channel matrix usually based only on time-delays of the finger Rake filter, S i is a vector associated with the ith user's signature code, and b D Rake i , is the estimated ith user's symbol. A GxG matrix W is inserted which will adaptively augment and improve the estimate of the channel inverse. In the following, we summarize the process in Algorithms 1, 2, and 3 to adaptively estimate the matrix W using the FastICA, Robust ICA and PCA algorithms, respectively.

## Algorithm 1 Adaptive Rake based FastICA method

1: procedure Initialization

2:

r ← M × N matrix of realization

3:

W = I G ← Initial demixing matrix

4: Itr ← number of iterations

5: γ ← Step Size

6: H ← the estimated channel matrix

7: g ( y ) = y 3 ← the nonlinear fcn

8:

Pre-Whitening

:

9:

r

←

V

∗

r

=

Λ

((

-

1

)

/

2

)

E

T

r

10:

For Loop

:

11:

i

←

1

. . . N

12:

r

←

WH r

H

(:

,

i

)

13:

W

+

←

E g

[[

(

Wr

)]

T

]

-

E g

[

′

(

Wr W

)]

14: Normalization :

15: W + ← W norm W / ( )

16: b D ICA i , (: , i ) ← S H i W¯ H H r

17: goto For Loop .

18:

close

;

19: Output :

20: b D ICA i , ← the estimated Symbols

## Algorithm 2 Adaptive Rake based RICA method

- 1: procedure Initialization 2: r ← M × N matrix of realization 3: W = I G ← Initial demixing matrix 4: Itr ← number of iterations 5: µ ← Step Size 6: H ← the estimated channel matrix 7: g ( y ) ← the gradient of the Kurtosis 8: Pre-Whitening : 9: r ← V ∗ r = Λ (( -1 ) / 2 ) E T r 10: For Loop : 11: i ← 1 . . . N 12: r ← WH r H (: , i ) 13: W + ← W + µ ( I G -g r ( ) ∗ g r ( ) H ) W 14: Normalization : 15: W + ← W norm W / ( ) 16: b D RICA i , (: , i ) ← S H i W¯ H H r 17: goto For Loop . 18: close ; 19: Output : 20: b D RICA i , ← the estimated Symbols

## Algorithm 3 Adaptive Rake based PCA method

```
Algorithm 3  Adaptive Rake based PCA n
==============================================================================
 1: procedure  INITIALIZATION
 2:     r <- M    x   N matrix of realization
 3:      W = I _ { G } <- Initial demixing matrix
 4:     Itr <- number of iterations
 5:     \gamma <- Step Size
 6:     H <- the estimated channel matrix
 7: Pre-Whitening:
 8:     r <- V * r = A((-1)/2) E T  r
 9: For Loop:
 10:     i <- 1 \dots N
 11:     r <- WHH ^ r (:, i)
 12:     W + <- W + \gamma (I _ { G } - r * r ^ H ) W
 13:  Normalization:
 14:     W + <- W /norm (W)
 15:     b _ { I _ { PCA } (:, i) <- S _ { WHH ^ r } }
 16:     goto For Loop.
 17:     close;
 18:  Output:
 19:     b _ { I _ { PCA } <- the estimated Symbols
 }
```

## 5 SIMULATION RESULTS

A series of extensive simulations are carried out in order to verify and evaluate the performance of the proposed adaptive filters and algorithms in the multipath downlink DS-CDMA system in the presence of AWGN. We summarize the case study results as follows. We assume a constant spreading gain, which is G = 63 for Gold Codes and G = 64 for Orthogonal Variable Spreading Factor (OVSF) codes. The received CDMA signal experiences five multipath

channels L=5 with delays of 0,1,2,3 ,4 chips, respectively. Also, we set the complex attenuation coefficients to represent the multipath channels, specifically, h 0 = 0 3684 + 0 5364 , . . i h 1 = 0 1982 + 0 0187 , . . i h 2 = 0 0237 + 0 5683, . . h 3 = 0 1112 + 0 0835 , and . . i h 4 = 0 2203 + 0 2756 , respectively. We use the . . i following model function for sub-Gaussian sources for which the source signals have a negative kurtosis sign:

$$g ^ { \text{SUB} } \left ( \hat { \mathbf b } \right ) = \hat { \mathbf b } - \left ( \tanh \left ( \text{Re} \left \{ \hat { \mathbf b } \right \} \right ) + j \tanh \left ( \text{Im} \left \{ \hat { \mathbf b } \right \} \right ) \right ) \quad \ \ ( 5 4 )$$

Monte Carlo Simulations have been run to verify the validity of the algorithms. We also use the signal-to-noise ratio (SNR) as a figure of merit which represents the ratio of the energy per symbol and the power spectral density (PSD) of the noise. Moreover, all the user symbols are assumed to be transmitted with the same power. Fig. 6 (a) and (b) show the simulation results of BER vs. SNR for the proposed detectors in contrast to the existing and conventional ones for the number of users K=30 and K=50, respectively. The other parameters were set as (i) number of symbols M=1000 and (ii) number of paths L=5, with the values of SNR in the range of -10 dB to 30dB.

Fig. 6 shows that the proposed algorithms improve the performance of the CDMA system. One observes that the blind multiuser detection based on FBII has resulted in the lowest BER, and thus it outperforms all other detectors. One also observes that the proposed algorithms work even in cases which cause difficulties for the LMMSE receiver, as in the high SNR ratio, and when the sample set is fairly small. Moreover, the performance of the blind multiuser detection degrades as the number of users increases as comparatively seen in Fig. 6 (b).

Furthermore, we have also evaluated the effect of the OVSF codes as depicted in Fig. 7. As in Fig. 7, it is generally the case that using the OVSF codes enhances the performance of the proposed methods.

In the WCDMA System case, we assume that the channel coefficients are h 0 = 0 3684 + 0 5364 , . . i h 1 = 0 1982 + 0 0187 , . . i h 2 = 0 0237 + 0 5683, . . h 3 = 0 1112+0 0835 , and . . i h 4 = 0 2203+0 2756 , respectively. Also, all user-specific . . i codes use two types of spreading codes, namely, Gold codes with spreading gain G=63 and OVSF (or Walsh-Hadamard) codes with spreading gain G=64.

In Fig. 8 and 9, we document and demonstrate the performance of the various methods in terms of BER for the WCDMA downlink scenario. We observe that the LMMSE is slightly better than some presented detectors under good SNR conditions. However, the proposed algorithm based on FB-II outperforms all detectors over all SNR depicted ranges and has again produced the lowest BER when compared to all other methods.

It is also worthwhile to compare the presented algorithms with a relatively large data sample set. Thus, Fig. 10 and Fig. 11 present the performance of the various detectors with fairly long sample set, namely, M=30 000 in each of the DS-CDMA and WCDMA systems. It is noted that the benchmark LMMSE detector performs much better for high SNR. It is plausible to assume that the LMMSE detector becomes better than other detectors under good SNR

Fig. 6 Average BER as a function of SNR for DS-CDMA downlink. Using Gold codes G=63. (a) Using 30 users (b) Using 50 users

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-748f8223-0ab1-43f4-8f96-8f9aae8a7751.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-5ed3253a-9cbd-4a0d-ba4e-acffc7455a60.jpg)

conditions. However, the proposed algorithm based on FB-II has exceeded the LMMSE detector at all SNRs less than 22 dB.

Finally, we evaluate the effect of the number of users and the size of the sample set on the performance of the proposed FB-II method in Figs. 12 and 13, respectively. In Fig. 12, the simulation results show the BER vs. SNR

Fig. 7 Average BER as a function of SNR for DS-CDMA downlink. Using OVSF codes G=64. (a) Using 30 users (b) Using 50 users

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-04d8f750-dd90-4f95-b2a6-c058c87bb8fb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-2d47e64c-caf7-4334-80cd-d8664073b587.jpg)

with various K users at 500 symbols for each user for blind multiuser detection based on the FB II detector. As expected, Fig. 12 shows that the FB-II detector decreases in performance as K, the number of users, is increased. Moreover, Fig. 13 shows the simulation results of BER vs. SNR with 30 users (K=30) for various data samples (M). The proposed FB-II algorithm appears robust

Fig. 8 Average BER as a function of SNR for WCDMA downlink. Using Gold codes G=63. (a) Using 30 users (b) Using 50 users

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-32847891-b2cb-47d4-83f6-bf63579fe6de.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-1744de61-9163-4d6c-8a21-1b8faac6d0af.jpg)

and performs resonably well, and it is obvious that its performance improves more consistently as M increases by mitigating the MIA.

Overall, the proposed variant detectors and algorithms perform well in solving the symbol estimation problem in the DS/WCDMA downlink system, especially when the size of the sample set is reasonably small.

Fig. 9 Average BER as a function of SNR for WCDMA downlink. Using OVSF codes G=64. (a) Using 30 users (b) Using 50 users

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-41374f60-8155-40c1-91c5-3846e0103014.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-02b0b867-4df6-4a59-829f-4c2c00b7fcb4.jpg)

## 6 CONCLUSION

We have presented formulations, derivations, and subsequent extensive simulations of various filtering algorithms within various structures for multi-user detection in CDMA based systems. We have developed three blind multiuser

Fig. 10 Average BER as a function of SNR for DS-CDMA downlink. For 30 users (a) Using Gold codes G=63. (b) Using OVSF codes G=64.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-d50df44a-5d2b-4200-befd-5509197fafe2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-3c036959-dc37-4cf6-83bc-225397e6e184.jpg)

detectors based on different filtering structures and three algorithms, namely, ICA, RICA and PCA. The results appear to show that the proposed structures perform well in the symbol estimation problem in DS/CDMA systems and more generally outperform all other detectors, including the LMMSE detector. Our results also show that MAI can be mitigated by the proposed

Fig. 11 Average BER as a function of SNR for WCDMA downlink. For 30 users (a) Using Gold codes G=63. (b) Using OVSF codes G=64.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-bd73a410-9654-4a1e-9f04-f8dd43594221.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-7bc2b9bc-3e70-4498-893d-25e780576ab0.jpg)

algorithms, particularly the proposed FB-II detector. Although the FB-II detector further improves as the size of the sample set increases, the results show that it performs well even when the sample sets are small- unlike the LMMSE detector. Finally, the proposed algorithms, unlike the adaptive LMMSE detector, do not require the spreading codes of the interfering users. While these

Fig. 12 Average BER as a function of SNR for various number of users K

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-4bf6c2a8-4141-4629-b51d-f2dcf2149521.jpg)

Fig. 13 Average BER as a function of SNR for various sample sets M

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0196v2-figure-19677584-2427-4fb2-a1aa-4aee2244d2b0.jpg)

detectors are more suitable for the downlink case, they can also be used in the uplink case as well.

## Acknowledgement

This research was supported in part by the NSF Grant ECCS-1549517.

## References

- 1. X. D. Wang and H. V. Poor, Wireless Communication Systems: Advanced Techniques for Signal Reception. Prentice Hall, 2004.
- 2. P. Comon, C. Jutten (eds.), Handbook of Blind Source Separation Independent Component Analysis and Applications (Academic Press, Oxford, 2010).
- 3. 'An Introduction to LTE'. 3GPP LTE Encyclopedia. Retrieved December 3, 2010.
- 4. LTE An End-to-End Description of Network Architecture and Elements. 3GPP LTE Encyclopedia. 2009.
- 5. 'Verizon launches its first LTE handset'. Telegeography.com. 2011-03-16. Retrieved 201203-15.
- 6. K. Waheed, F. Salem, Blind information-theoretic multiuser detection algorithms for DSCDMA and WCDMA downlink systems, IEEE Trans. Neural Netw., vol. 16, no. 4, pp. 937948, Jul. 2005.
- 7. R.C. de Lamare, Diniz, P.S.R., Blind Adaptive Interference Suppression Based on SetMembership Constrained Constant-Modulus Algorithms With Dynamic Bounds ', IEEE Trans. Signal Process., vol. 56, no. 5, pp. 1288 1301 2013
- 8. L. Hu, X. Zhou, L. Zhang, 'Blind Multiuser Detection Based on Tikhonov Regularization', IEEE Communications Letters, vol. 15, no. 5, May 2011.
- 9. S. Verdu, Multiuser Detection. Cambridge U.K.: Cambridge Univ. Press, 1998.
- 10. G. B. G. Zhang and L. Zhang 'Group-blind intersymbol multiuser detection for downlink cdma with multipath', IEEE transaction on Wireless Communications, vol. 4, no. 2, pp.434 -443 2005
- 11. R. C. de Lamare , M. Haardt and R. Sampaio-Neto 'Blind adaptive constrained reducedrank parameter estimation based on constant modulus design for CDMA interference suppression', IEEE Trans. Signal Process., vol. 56, no. 6, pp.2470 -2482 2008
- 12. Z. Albataineh, F. Salem, New blind multiuser detection DS-CDMA algorithm based on Extension of Efficient FAST Independent Component Analysis (EF-ICA). 4th International Conference on Intelligent Systems Modelling and Simulation (ISMS), pp. 543-548, 2013.
- 13. Z. Albataineh, F. Salem, New blind multiuser detection DS-CDMA algorithm using simplified fourth order cumulant matrices IEEE International Symposium on Circuit and System (ISCAS), pp. 1946-1949, 2013.
- 14. Jean-Franois Cardoso, High-order contrasts for independent component analysis, Neural Computation, vol. 11, no 1, pp. 157-192, Jan. 1999.
- 15. X. Liu, R. B. Randall, A New Efficient Independent Component Algorithm: Joint Approximate Diagonalization of Simplified Cumulant Matrices. The 16th National Congress of Australian Institute of Physics, Canberra, Australia, 2005.
- 16. V. Zarzoso and P. Comon, Robust Independent Component Analysis by Iterative Maximization of the Kurtosis Contrast With Algebraic Optimal Step Size, IEEE Transactions on Neural Networks, vol. 21, no. 2, pp. 248-261, 2010.
- 17. Hyvarinen. A, 'Fast and robust fixed-point algorithm for independent component analysis'. IEEE Transactions on Neural Network, Vol. lo, no. 3, pp. 626-634, May 1999.
- 18. E. Oja and Z. Yuan 'The fastICA algorithm revisited: Convergence-analysis', IEEE Trans. Neural Netw., vol. 17, no. 6, pp.1370 -1381 2006

- 19. X. G. Doukopoulos and G. V. Moustakides 'Adaptive power techniques for blind channel estimation in CDMA systems', IEEE Trans. Signal Process., vol. 53, no. 3, pp.1110 -1120 2005
- 20. X. Wand and H. V. Poor, Blind multiuser detection: A subspace approach', IEEE Trans. Inform. Theory, vol. 44, pp.677 -690 1998
- 21. X. Wang and A. Host-Madsen, 'Group-blind multiuser detection for uplink CDMA', IEEE J. Sel. Areas Commun., vol. 17, no. 11, pp.1971 -1984 1999
- 22. S. Cui, M. Kisialiou , Z.-Q. Luo and Z. Ding 'Robust blind multiuser detection against signature waveform mismatch based on second-order cone programming', IEEE Trans. Commun., vol. 4, no. 4, pp.1285 -1291 2005
- 23. A. Kachenoura , L. Albera , L. Senhadji and P. Comon 'ICA: A potential tool for BCI systems', IEEE Signal Process. Mag., vol. 25, no. 1, pp.57 -68 2008
- 24. P. Loubation and E. Moulines, On blind multiuser forward link channel estimation by the subspace method: Identifiability results, IEEE Trans. Signal Process., vol. 48, no. 8, pp. 23662376, Aug. 2000.
- 25. E. Oja and Z. Yuan 'The fastICA algorithm revisited: Convergence-analysis', IEEE Trans. Neural Netw., vol. 17, no. 6, pp.1370 -1381 2006
- 26. http://www.bloomberg.com/slideshow/2013-09-19/countries-with-the-most-4g-mobileusers.html,slide1
- 27. http://gigaom.com/2013/09/20/mapping-out-the-worlds-lte-coverage-its-in-fewerplaces-than-you-think/
- 28. J. W. H. Kao, S. M. Berber, and V. Kecman, Blind multiuser detection of a chaosbased CDMA system using support vector machine, in Proc. IEEE 10th Int. Symp. Spread Spectrum Tech. Applicat., Aug. 2008, pp. 5862.
- 29. T. Yang and B. Hu 'Blind multiuser detection based on kernel approximation', Advances in Neural Networks, vol. 3973, pp.94 -101 2006
- 30. J. Kao, S. Berber, and V. Kecman, 'Blind Multi-User Detector for Chaos-Based CDMA Using Support Vector Machine,' IEEE Neural Networks, vol. 21, pp. 1221 1231, 2010.
- 31. Jiang Zhang, Hang Zhang ; Zhifu Cui, Dual-antenna-based blind joint hostile jamming cancellation and multi-user detection for uplink of asynchronous direct-sequence codedivision multiple access systems, IEEE Communications IET, vol. 7, 2013
- 32. X. Liu, K. C. Teh and E. Gunawan, 'A blind adaptive MMSE multiuser detector over multipath CDMA channels and its analysis,' IEEE Trans. Commun., vol. 7, no. 1, pp. 90-97, Jan 2008.
- 33. M. Hubert, P. J. Rousseeuw, and K. Branden, ROBPCA: A new approach to robust principal component analysis, Technometr., vol. 47, no. 1, pp. 6479, 2005.
- 34. X. X. Zhang and T. S. Qiu, 'Blind multiuser detection based on improved Infomax and Fast ICA,' 2010, 2nd International Conference on Advanced Computer Control, Shenyang, 2010, pp. 476-479.
- 35. Y. Washizawa , Y. Yamashita , T. Tanaka and A. Cichocki 'Blind extraction of global signal from multi-channel noisy observations', IEEE Trans. Neural Netw., vol. 21, no. 9, pp.1472 -1481 2010
- 36. A. Bayati, S. Prakriya and S. Prasad, Semi-blind space-time receiver for multiuser detection of DS/CDMA signals in multipath channels, The Institution of Engineering and Technology, vol. 153, no. 3, 2006.
- 37. K. Waheed and F. M. Salem, Blind Source Recovery: A Framework in the State Space, Journal of Machine Learnin Research, 4 (2003) pp. 1114-1146.
- 38. Koivisto T, Koivunen V. Blind Dispreading of Short-Code DS-CDMA Signals in Asynchronous Multi-user Systems. Signal Processing, 2007, 11(87): 2560-2568.
- 39. L. De Lathauwer and A. de Baynast 'Blind deconvolution of DS-CDMA signals by means of decomposition in rank-(1,L,L) terms', IEEE Trans. Signal Process., vol. 56, no. 4, pp.1562 -1571 2008
- 40. W. N. Yuan , Y. F. Tu and P. Z. Fan 'Optimal training sequences for cyclic-prefix-based single-carrier multi-antenna systems with space-time block-coding', IEEE Trans. Wireless Commun., vol. 7, no. 11, pp.4047 -4050 2008
- 41. L. Hu, X. Zhou, L. Zhang, 'Blind Multiuser Detection Based on Tikhonov Regularization', IEEE Communications Letters, vol. 15, no. 5, May 2011.
- 42. D. Nion and L. De Lathauwer, A Block-Component Model Based Blind DS-CDMA Receiver, IEEE Trans. on Signal Processing, Vol. 56, No. 11, pp. 5567-5579, Nov. 2008.

- 43. S. C. Douglas and M. Gupta 'Scaled natural gradient algorithms for instantaneous and convolutive blind source separation', Proc. ICASSP, vol. II, pp.637 -640 2007
- 44. Shun-Ichi Amari, Natural gradient works efficiently in learning, Neural Computation, v.10 n.2, p.251-276, Feb. 15, 1998
Zaid Albataineh Fathi M. Salem
A Blind Adaptive CDMA Receiver Based on State Space Structures
Code Division Multiple Access (CDMA) is a channel access method, based on
spread-spectrum technology, used by various radio technologies world-wide. In
general, CDMA is used as an access method in many mobile standards such as
CDMA2000 and WCDMA. We address the problem of blind multiuser equalization in
the wideband CDMA system, in the noisy multipath propagation environment.
Herein, we propose three new blind receiver schemes, which are based on state
space structures and Independent Component Analysis (ICA). These blind
state-space receivers (BSSR) do not require knowledge of the propagation
parameters or spreading code sequences of the users they primarily exploit the
natural assumption of statistical independence among the source signals. We
also develop three semi blind adaptive detectors by incorporating the new
adaptive methods into the standard RAKE receiver structure. Extensive
comparative case study, based on Bit error rate (BER) performance of these
methods, is carried out for different number of users, symbols per user, and
signal to noise ratio (SNR) in comparison with conventional detectors,
including the Blind Multiuser Detectors (BMUD) and Linear Minimum mean squared
error (LMMSE). The results show that the proposed methods outperform the other
detectors in estimating the symbol signals from the received mixed CDMA
signals. Moreover, the new blind detectors mitigate the multi access
interference (MAI) in CDMA.