s
1
=
s
1,1
,s
2,1
,...,s
C,1
T
Cx1 vector
r
m
=
s
1
,s
2
,...,s
K
h
1
0 .. 0
0 h
2
.. 0
: : .. :
0 0 .. h
K
b
1,m
b
2,m
:
b
K,m
+
n
1
n
2
:
n
K
r
m
= [s
1
h
1
,s
2
h
2
,...,s
K
h
K
]b
m
+ n
m
(7)
Equation (7) can be represented in a more compact
form:
r
m
= Gb
m
+ n
m
(8)
where theCxK matrix G is assumed full rank. We can
see the similarity between the IDMA model of equa-
tion (8) and the N-ICA model of equation (3). The
goal of the Noisy-ICA based IDMA detection is to
recover the symbol vector b
m
for each user k without
knowing the parametric form of G which depends on
the channel coefficients. The objective is to estimate
the filter weight w such that the variable at the output
of the filter is one of the ICs(source signal):
y
m
= w
T
r
m
(9)
If BPSK modulation is used, the symbol of desired
user k can be obtained by using this decision formula:
ˆ
b
k,m
= sgn(w
T
r
m
) (10)
4 N-ICA ESTIMATION
ALGORITHM
The proposed system is a hybrid structure composed
of two parts where a classical IDMA receiver is com-
bined with a N-ICA block as shown in figure 1. Block
IDMA, described in the previous section, works for a
number of iterations (it) after which the block N-ICA
takes over. The proposed N-ICA will act as a post pro-
cessor attached to an IDMA receiver in the presence
of noise. The aim of our N-ICA block is to avoid con-
tinuous tracking of channel state information. In this
section, we will derive estimation algorithms for the
proposed N-ICA post processor in IDMA context.
4.1 Principal Component Analysis
based Processing
The Principal Component Analysis (PCA) based part
of the model consists of whitening the input signals.
This step of processing is achieved by using the PCA
in (Davies, 2004) to extract the Principal Components
(PCs). It is based on the diagonalization concept of
the input signals covariance matrix. This can be done
for the noiseless case as follows
Y = Λ
−1/2
U
T
GB (11)
where the matrix U corresponds to the Eigen vec-
tor of the data covariance matrix C and the diago-
nal matrix Λ that contains the related Eigen values:
Λ
−1/2
= diag[λ
−1/2
1
,λ
−1/2
2
,...,λ
−1/2
n
]. This PCA pro-
cessing can be extended to noisy data using bias re-
movaltechnique (Ekici,2004). In the regular ICA pro-
cess, the covariance matrix of the noise free data r
(nf)
m
can be given by:
C = E{r
(nf)
m
(r
(nf)
m
)
T
} = GG
T
(12)
On the other hand, the covariance matrix of the ob-
served noisy data can be written as:
Γ = E{r
m
(r
T
m
)} = GG
T
+ σ
2
I = C+ C
n
(13)
where σ
2
is the noise power and C
n
is the diagonal
noise covariance matrix. In the noise bias removal
technique, the Eigen values and vectors of matrix Γ−
C
n
is used for whitening instead of matrix Γ which
is called quasi-whitening (Hyvarinen,1999). In fact,
quasi whitening can be performed on the noisy data
as follows:
z = (Λ− σ
2
I)
−1/2
U
T
r
m
(14)
The covariance matrix of quasi white data is then
given by :
E{zz
T
} = I + σ
2
(Λ− σ
2
I)
−1
(15)
From (15), we notice that the covariance matrix is
different from the identity matrix. Therefore, we have
to take into account the non-whiteness of the data.
This is achieved by using the fast ICA algorithm that
is presented in the next subsection.
4.2 Fast ICA Algorithm
The purpose of this work is to establish a new
scheme in which the system can take into account
such random deformations in the detection step.
To improve the performance, the presence of the
noise should be reduced to the minimum using the
extracted PCs without additional prior knowledge of
their statistical properties. This is the purpose of the
ICA based part of the model. Therefore, the ICA
model should include a noise term as well in its linear
transform matrix. We have used two algorithms for
detecting and separating the received signals: IDMA
algorithm (Schoeneich,2005), and the fastICA in
IDMA in (Hamza, 2009). The second ICA approach
that we present here is our contribution to take into
account the noise in the ICA model. This means that
the bias due to noise should be removed, or at least
reduced. The N-ICA algorithm performs as follows:
Let k be the desired user, r
m
,m = 1,..,M the re-
ceived block data and
ˆ
b denotes the estimate of b.
BLIND DETECTION IN IDMA SYSTEMS
139