Maximal Ratio Combining Using Channel Estimation in Chaos Based Pilot-Added DS-CDMA System with Antenna Diversity

This paper presents an adaptive multiuser transceiver scheme for DS-CDMA systems in which pilot symbols are added to users’ data to estimate complex channel fading coefficients.The performance of receiver antenna diversity with maximal ratio combining (MRC) technique is analyzed for imperfect channel estimation in flat fading environments. The complex fading coefficients are estimated using least mean square (LMS) algorithm and these coefficients are utilized by themaximal ratio combiner for generating the decision variable. Probability of error in closed form is derived. Further, the effect of pilot signal power on bit error rate (BER) and BER performance of multiplexed pilot and data signal transmission scenario are investigated. We have compared the performance of added and multiplexed pilot-data systems and concluded the advantages of both systems. The proposed CDMA technique uses the chaotic sequence as spreading sequence. Assuming proper synchronization, the computer simulation results demonstrate the better bit error rate performance in the presence of channel estimator in the chaotic based CDMA system and the receiver antenna diversity technique further improves the performance of the proposed system. Also, no channel estimator is required if there is no phase distortion to the transmitted signal.


Introduction
The fundamental phenomenon which makes reliable wireless transmission difficult is time varying multipath fading.It is this phenomenon which makes wireless transmission a challenge when compared to fiber, coaxial cable, line-ofsight microwave, or even satellite transmission.Severe rapid fluctuation in received signal due to multipath environment makes it extremely difficult for the receiver to determine the transmitted signal [1].In multipath fading environment, BER improvement of order 10 −2 to 10 −3 using typical modulation and coding schemes may require SNR improvement up to 10 dB.Since there is power transmission limitation in wireless communication systems, hence BER improvement cannot be achieved by increasing the transmission power [2].Spread spectrum techniques, channel estimation algorithms, and diversity techniques [3] are widely used methods for mitigating the effect of fading without increasing the transmission power.
CDMA system has received much attention in wireless communication system because the available spectrum is shared by all the users.In CDMA system each user is distinguished from others by uniquely defined spread sequence.Chaotic and Walsh sequences are widely used for identifying the correct user.Since interuser interference degrades the performance of CDMA system, the multiuser channel estimation gains a significant interest in recent years.Several methods have been proposed for estimating the channel based on short and long spreading codes [4][5][6][7][8][9][10].
In pilot assisted channel estimation [11][12][13][14] technique a predefined data sequence, known as pilot, is transmitted along with the actual message signal.These pilot sequences are used for estimating the fading coefficients and synchronization.Although there are various ways for transmitting the pilot signal, for example, sending pilot signal at the beginning or middle of each packet, in our work we have assumed that the pilot signal and actual users messages are transmitted simultaneously using separate spreading signals, which is a standard method in CDMA-2000 system.In this research work, it is assumed that the pilot and users data have the same power and they are not multiplexed but added to each other, which removes the assumption that pilot and users data have the same fading and noise effects.Additionally, a separate mathematical analysis has been done for different pilot and data powers and multiplexed pilot-data case.
Pilot signal and user's message have to be separated for estimating the channel coefficients and user's data; therefore spreading sequences should have impulse-like autocorrelation function and zero cross correlation values.Chaotic sequences possess good properties for pilot aided CDMA system.They are aperiodic sequences with mixed deterministicstochastic nature and are very sensitive to initial conditions.Therefore, infinite number of chaotic signals can be generated just from one chaotic map with impulse-like autocorrelation function with low cross correlation values [15].Theoretical analysis of the chaos based CDMA systems has been done in many research works [16][17][18][19][20].In these works channel coefficients are assumed to be real and positive.Complex channel gains are considered in [21][22][23] and perfect estimation of fading channel coefficients is assumed.Performance of these systems can degrade if fading coefficients are complex and the correct channel coefficients are not known at the receiver.Various adaptive equalization algorithms and corresponding BER performance for CDMA system with binary spreading codes have been investigated in [24][25][26][27][28] but none of these were used in chaos based CDMA system.Further, complex fading coefficients are estimated using Bayesian estimator in [14,29] and with least mean square (LMS) in [30,31] for chaos based CDMA systems.
Further, in [30,31], the received pilot signal is modified before using it for estimating the channel coefficients using LMS estimator.This paper shows analytically that the inclusion of spreading sequences increases the value of minimum mean square error for the LMS estimator, which can be minimized by normalizing the received pilot signal.Also, only added pilot-data systems are studied in [30,31], whereas this paper compares, mathematically as well as using simulation results, the added and multiplexed pilot-data system.Moreover, perfect estimator and no estimator cases are studied in this paper, which demonstrate the need of channel estimator in the systems which have complex fading coefficients.Lastly, the first-order Markov channel model includes the Doppler effect in the channel; and we have shown the effects of Doppler frequency on BER performance of the system using simulation results.
Diversity is principally used to combat fading.Multipleantenna diversity is the most widely used diversity technique to mitigate the fading effect [32].The method of fading mitigation using antenna diversity in chaos based CDMA system is discussed in [33] with known channel coefficients.The basic idea is that if different copies of the same signal are available then there is a high probability that at least one of them is of a good quality.Choosing the best copy and rejecting all others are not the optimal solution.This brings us to the problem of choosing the best way to combine all of them [11,34].MRC [35] technique is the most effective technique for combining.Channel estimator plays vital role in MRC because all the received signals must have the same phase before combining.Therefore channel estimators play dual role here; firstly they remove the fading effect from the received signals which are used for generating decision variable at each branch.Secondly they ideally maintain the zero phase signal at every branch; hence all the input signals to maximal ratio combiner have the same phase, to generate final decision variable.
The objective of this paper is to show the effect of imperfect channel estimation on the performance of chaos based DS-CDMA system in flat fading environment which can be extended for binary spreading sequence also.We have derived generalized BER expressions in closed form for imperfect channel estimation with receiver antenna diversity using MRC to optimize the decision.BER expression shows that the probability of error is inversely proportional cosine of the estimated phase error.LMS algorithm [36] is used for estimating the channel coefficients with the help of chaotic signals in CDMA system with receiver antenna diversity.All the received signals from different antennas are combined together using MRC method to generate BER curves and the results are compared with perfect channel estimation case.We assume a block fading channel; that is, length of data signal is less than the coherence time of the channel.Under these assumptions the estimation algorithm works in two modes, that is, training mode and tracking mode.In training mode channel coefficients are estimated with the help of known pilot signals.Then in tracking mode actual message is detected with the help of estimated channel coefficients.Finally we have shown that the BER expression derived in [16] is the special case of BER expression derived in this paper.During this complete process proper synchronization is assumed.Robust synchronization techniques for chaos based DS-CDMA system are shown in [37,38] and effect of Rayleigh fading coefficients on synchronization process is analyzed in [17].
This paper is organized as follows.Section 2 presents the baseband representation of proposed chaos based DS-CDMA system with channel estimation in the presence of complex noise and fading.Section 3 deals with the generation of decision variables.The analytical performance of MRC technique is shown in Section 4. Simulation results are presented in Section 5. Section 6 reports some concluding remarks.

System Model
Figure 1 shows the baseband representation of the chaos based CDMA system.The baseband signal  ()  for th user at th chip instant is given by where  ()  is th user symbol at th time instant, is th chip of th user chaotic spreading sequence or chip within an information bit (i.e.,  = 1, . . ., 2), and 2 is the spreading factor.
Similarly baseband pilot signal is represented by User "1" data Impulse response

Accumulating over spreading factor chips
Received data symbol

Received pilot symbol
Accumulating over spreading factor chips where is the pilot symbol at th time instant and is the chaotic spreading sequence for pilot.Pilot symbol is usually the series of ones; that is, Here we assume that spreading sequence of pilot and data have the same chip power   .Bit error rate results for different pilot and data power for added pilot-data and multiplexed pilot-data case have been included in Appendix C.
The flat fading channel between transmitter and receiver antenna  at time  is denoted by complex variable ℎ , which includes the amplitude variation ( , ) and phase variations ( , ) caused by the channel, expressed as We assume that noise is affecting the signal on chip level, whereas fading is happening on bit level.Then th chip of received signal

𝑚(𝑔) 𝑘
at antenna  of th user at time  from  transmission users is given by where  , denotes the complex additive white Gaussian noise for th receiving antenna with power spectral density equal to  0 .Clearly in the above equation every 2 chip is representing one bit duration.

Correlation Filter and Decision Variable
In Figure 1, the received signal is processed at every antenna to separate the pilot and data symbols.Then the extracted pilot symbols are used to estimate the channel coefficients.These estimated channel coefficients are utilized by the maximal ratio combiner for generating the decision variable.[39] is used to extract desired signal from the noisy signal.In Figure 1, received signal is multiplied by the locally generated pilot chaotic sequence and th user chaotic sequence to extract the pilot and th user data, respectively.After multiplication of chaotic sequences, the resultant sequences are accumulated over the 2 interval.Hence the two correlator outputs are given by

Pilot and Data Extraction. Correlation filter
where  () ,data and  () ,pilot are the extracted data and pilot in th antenna branch at time  of th user.

Channel Estimation and Decision
Variable.This extracted pilot is used as the input to the LMS filter and the resultant algorithm can be represented as follows: where â,   φ, are the estimated channel coefficients which are initialized with zero,  is the step size, and   is known as error function in the LMS filter. () ,pilot and ∑ =2 =2(−1)+1 ( ()  ) 2 are known as desired signal and known signal in estimation theory [25].The minimum mean square error (MMSE)  min for this algorithm is (see Appendix A for derivation) Channel coefficients can not be estimated properly with such an enormous amount of minimum mean square error.This value of  min can be reduced by normalizing the desired signal and known signal by ∑ In this case the new LMS algorithm is represented by The decision variable is obtained by multiplying the extracted data symbol  () ,data with the conjugate of the estimated channel coefficients.Hence the decision variable  ()  after processing the signal at antenna  at time  is given by where * denotes the conjugation operation.Putting value of  () ,data from (5) to the above equation and after separating mathematical equations of various terms, that is, desired signal, multiuser interference, and channel noise, we get the above equation in the following form:

Performance Analysis of MRC Technique
In MRC technique all intermediate decision variables  ()  ,  = 1, . . .,  are added together to generate the final decision variable, which is given by where () is the final decision variable used for symbol detection and Re is the real part of complex function.
Putting value of  ()  from (11) to the above equation, we get final decision variable as follows: Im ( , ) where Re( , ) and Im( , ) are the real and imaginary parts of complex AWGN noise.Sign of this decision variable is the estimated symbol; that is, The probability of error ( () (  )) for th bit of th user is given by [15]  () where Pr and erfc are the probability operations and complementary error function, respectively.Solving the above equation (see Appendix B for derivation) and omitting the subscript  that is related to bit under investigation, we have the final equation of  () () as follows: where  = var[( and   = 2  .Hence the probability of error depends on the number of users (), spreading factor (2), estimated channel coefficients (â  , φ ), the signal to noise ratio (  / 0 ) and total number of receiving antennas ().Clearly probability of error increases as the number of users in the system increases and decreases as the spreading factor, signal to noise ratio, and/or number of receiving antennas increases.Further, the error probability depends on the cosine value of the difference between estimated phase and actual phase of the complex channel coefficients.Therefore the performance of the system is greatly dependent on the performance of the channel estimation algorithms.From this expression it is clear that the probability of error is minimum for perfect channel estimation where the difference between actual and estimated phase is zero and then increases gradually as the difference increases.
For perfect channel estimation case (17) reduces to If we assume that fading coefficients are real and positive (i.e.,   = 0 with no channel estimator in the system (i.e., ĥ, = 1 0 ), then we can neglect the pilot processing part in Figure 1 and total numbers of interuser interferers are reduced to ( − 1) in the absence of pilot signal.In this case (17) will become Thus the above BER expression which is derived in [16] and used in research works [16,17] is a special case of (17).From (17) and (19) we can see that the performance of the system deteriorates significantly in the presence of complex channel coefficients.Removing pilot interference term from (18), that is, changing  to ( − 1) and assuming only one receiving antenna, then from ( 18) and ( 19) we can observe that, for real fading coefficients, that is, when channel is not causing any phase distortion, the BER expression of the no estimator case is the same as the perfect estimator case.Therefore no channel estimator is required for real fading coefficients.
In general, the BER of the system ( ff+AWGN ) can be expressed as [40] where the subscript ff denotes the flat fading and () is the probability density function of Rayleigh distribution.

Simulation Results
Several cases are simulated to demonstrate the performance of the proposed pilot-added CDMA system with BPSK signaling.In the simulation we have used the first-order Markov process to model the fading process of channel [5,41,42], which is described as where V , is the complex Gaussian process for antenna  at time .In the simulation it is assumed that fading processes are independent to each other with unit variances.  is the correlation coefficient that depends on maximum Doppler frequency   and is defined as where  0 (⋅) is the Bessel function of first kind and zerothorder and   is the signaling rate.
In the simulation, the value of the spreading factor 2 is 50.In LMS algorithm, the order of the filter is one and the step size  is set to 0.9.We have chosen the following Chebyshev polynomial function to generate the chaotic sequence because this chaotic generator outperforms many other chaotic sequence generators [43]. where denotes the th chip value of th user and pilot.In the simulation results SISO (single input, single output) and SIMO (single input, multiple output) represent the one receiving and two receiving antennas, respectively, with one transmitting antenna.We have used the MRC technique to generate the final BER curve in SIMO case throughout the simulation.Similarly "LMS BER," "True channel BER," and "Simple BER" represents the bit error rate with LMS channel estimation, perfect channel estimation, and no channel estimation case respectively.
Figure 2 compares the performance of no estimator case with perfect estimator case for real fading coefficients with one user and pilot signal and having fading rate of     = 0.01.For simulation, the real fading coefficients are generated by taking absolute value of (21).For this case the performance of the system used in [16,17] is matched with the perfect estimator case.Also there is a significant gain in the BER performance with two receiving antennas' diversity system over single antenna system.Figure 3 demonstrates the BER comparison between various cases with one and two receiving antennas for     = 0.01 in the presence of one user and pilot signal with complex fading coefficients.In this case the BER for no estimator case is 0.5 for complete SNR range, which was quite good for real fading coefficients.This is happening because, for complex fading coefficients, no estimator case is not considering the phase change due to channel and simply giving a random guess about the transmitted bit, which has the probability of 0.5.From the plot it is clear that the channel estimator improves the performance of the proposed CDMA system.In this figure "LMS SISO-channel 1" and "LMS SISO-channel 2" denote the BER due to individual antenna processing.Almost 5 dB performance improvement is achieved using MRC technique at SNR = 15 dB with two receiving antennas as compared to one antenna system.This performance is further improved for higher SNR conditions.
Figure 4 shows the performance of two receiving antennas' diversity system for different fading rates for 1 user and pilot system.There is almost 2 dB performance decrement in BER performance with increase in fading rate from 0.01 to 0.1.
Figure 5 depicts the performance of LMS estimator in two receiving antennas' diversity system for different number of users.The performance of the system degrades as the number of users increases.
Figure 6 illustrates that the performance of two receiving antennas' diversity system improves with higher values of spreading factor.On comparing this result with Figure 3, we can observe that the performance of system with 10 users with higher spreading factor, that is, 300, is equivalent to performance of the system with monouser with lower spreading factor value, that is, 50.
Effect of higher power of pilot is shown in Figure 7.Here we assumed that pilot power is 6 dB higher than the data and     = 0.1.It can be seen in the BER curves of estimation cases that there is an intersection point between equal pilotdata and higher pilot-data power case.Before intersection point, that is, at low SNR conditions, the performance of the higher pilot power case has better performance than equal  power case and the performance degrades after the intersection points.From (C.1) it can be observed that pilot power increases the interuser interference term, but at the same time performance of the channel estimator improves because of higher pilot power.This performance improvement compensates the interuser interference caused by higher pilot power.But at higher SNR, the estimation algorithm can achieve minimum mean square value at low pilot power; therefore higher power is only producing the interuser interference.Since the higher pilot power transmission is not required at high SNR condition, hence this performance degradation can be improved by controlling the power transmission.Finally in Figure 8, multiplexed pilot-data scenario for two receiving antennas is simulated with four users and pilot signal for     = 0.1.If we assume that pilot and data have the same fading and noise effect, then the performance of this system improves with increase in the pilot power transmission because the pilot is not creating any interuser interference.Similarly for perfect channel estimator case, this system has the same BER curves for higher power transmission and equal power transmission cases.Also, for higher pilot power transmission the BER performance of the system with LMS estimator is nearly equal to the perfect estimator case.But this performance improvement can degrade if the pilot and data do not have the same fading and noise effect.

Conclusion
In this paper, we propose the pilot-added transceiver for CDMA system in flat fading channel in which chaotic sequences are used for spreading the data symbols of various users.In simulation, performance of LMS channel estimation algorithm is compared with the perfect estimator case and with simple CDMA system in which no estimation algorithm is used, for the same pilot and data powers.There is huge BER performance improvement in the proposed CDMA system using channel estimation algorithm.However, for real fading coefficients, that is, when channel is not causing any phase distortion, then the BER expression of no estimator case is the same as the perfect estimator case as shown in simulation results.Therefore no channel estimator is required for real fading coefficients.The performance of the proposed system is then compared with the two receiver antennas' diversity system.Simulation results show that there is further improvement in the system using diversity technique.We have also verified that BER performance is slightly better for low Doppler spreading.It is shown that there is an increase in bit error rate with increase in pilot power, but this performance degradation is compensated by the performance improvement of the estimator at low SNR conditions.Since the high transmission pilot power is not required for high SNR conditions, therefore if the transmission power is controlled for high SNR conditions, then the proposed system's performance can be improved for high SNR case also.Separate mathematical equation is derived for the multiplexed pilot-data case.In this case we assume that pilot and data have the same fading effects.For this system the BER equation does not depend directly on the pilot power, but the performance of such system increases with increase in pilot power because the performance of estimator increases with power, which results in improvement in BER.Further the simulation results show that the BER performance of multiplexed data-pilot system using LMS channel estimator is closely matched to the perfect channel estimation case when pilot power is 6 dB higher than the data signal power.In final conclusion we can say that, for multiplexed pilot-data case, if pilot and data symbols do not have the same fading and noise effects then the added pilot-data system performs better than the multiplexed system and vice versa.16) and after rearranging the equation we have (17).

Figure 2 :
Figure 2: Simulated BER comparison of no channel estimator case and perfect estimator case for real fading coefficients in 1 user and pilot system with added pilot-data case,     = 0.01 and 2 = 50.

Figure 3 :
Figure3: Simulated BER comparison of various methods in chaos based CDMA system for complex fading coefficients in 1 user and pilot system with added pilot-data case,     = 0.01 and 2 = 50.

1 Figure 4 :Figure 5 :
Figure4: Simulated BER comparison of different fading rates with two receiving antennas' diversity in chaos based CDMA system for 1 user and pilot system with added pilot-data case, 2 = 50.

Figure 6 :
Figure 6: Simulated BER comparison of different values of spreading factor (SF) for 10 users and pilot for two receiving antennas' diversity in chaos based CDMA system for added pilot-data case,     = 0.1.

Figure 7 :
Figure 7: Simulated BER comparison of different pilot power with two receiving antennas' diversity in chaos based CDMA system for 4 users and pilot system with added pilot-data case,     = 0.1 and 2 = 50.

Figure 8 :
Figure 8: Simulated BER comparison of different pilot power with two receiving antennas' diversity in chaos based CDMA system for 4 users and pilot system with multiplexed pilot-data case,     = 0.1 and 2 = 50.