An APP Algorithm for Fading Channels using Forward-Only Prediction

An APP Algorithm for Fading Channels using Forward-Only Prediction Anders ansson,, or Aulin, and eith M. Chugg Department of Computer Engineering Communication Sciences Institute Chalmers University of echnology University of Southern California 42 96 Göteborg Los Angeles, CA 90089-2565, USA {ahansson, tor}@ce.chalmers.se {hansson, chugg@usc.edu} ABSRAC A multiple-antenna receiver is proposed for achieving a diversity effect that partly overcomes the severity of continuous-time frequency-flat Rayleigh fading channels, and fast fading channels in particular. he error performance is further improved by bit interleaving and channel coding, where the encoders/channel is viewed as a serially concatenated system: a convolutional code (CC constitutes the outer code, whereas a differential encoder and the fading channel (having truncated memory form a joint inner code. In order to obtain a feasible detector structure, it is desirable to perform iterative decoding, which implies the need for an APP (a posteriori probability algorithm that can operate over fading channels. For this purpose, the well-nown BCJR algorithm (due to Bahl, Coce, Jeline, Raviv [] is generalized to handle channels having memory. Numerical results indicate that iterative decoding becomes more powerful when the exploited channel memory depth is extended, but the chief performance gain is due to space diversity provided by the multiple antennas. INRODUCION Shortly after Berrou et al. announced that interleaved parallel concatenated convolutional codes (PCCC, socalled turbo codes [2], yield performance close to the theoretical Shannon capacity limits, serially concatenated codes (SCC were shown to offer comparable (and in some cases even superior performance [3]. he nearoptimal performance over additive white Gaussian noise (AWGN channels is partly due to the interleaver embedded in the code structure, which forms an overall concatenated code with very large memory. he analysis in [3] revealed that the use of a recursive systematic convolutional encoder (RSC as inner encoder guarantees an interleaver gain [3]. his fact further indicates that continuous phase modulation (CPM [4] should be equally suitable as an inner encoder (CPM is a recursive trellis code, achieving bandwidth efficient modulation at the same time. Such a system, SCCPM, has recently proved to give performance comparable to turbo codes both over AWGN channels [5] as well as over piecewise constant memoryless Rayleigh fading channels [6]. o simplify the treatment of the fading channel, the latter study assumes the detector has somehow attained perfect channel state information (CSI. Even without the ideal assumption of perfect CSI, excellent performance has been achieved for piecewise constant memoryless Rayleigh channels [7]. Such channel models have unfortunately a somewhat limited validity [8]. A few results have also been presented for PCCC over slow and fast Rayleigh fading channels [9]. he encoder proposed in [9] has the drawbac of being highly complex (it consists of two RSCs, two interleavers, one deinterleaver, and one differential encoder, leading to excessive decoding complexity. his paper analyzes an SCCPM system over continuous-time fading channels, i.e., over channels where the rate of variation of the fading process cannot be neglected [0]. o somewhat reduce the receiver complexity, only the most elementary non-trivial CPM member, namely MS (minimum shift eying [], [2] is treated. SYSEM AND CANNEL DESCRIPION An all-zero tail is first appended to the information bits, before the bits are fed to an outer CC. his is done in order to provide perfect nowledge of the terminal state of the outer trellis. Next, the coded pacage of bits is bit interleaved, MS modulated, and transmitted over the equivalent lowpass channel, which comprises D diversity lins. In each lin the transmitted sequence s (; t i of waveforms is perturbed by two time-varying random processes, one is the multiplicative fading vector f ( t (dimension D by, while the other is the additive thermal noise vector n ( t (dimension D by. he variable i in s (; t i refers to the hypothesis label, i.e., {,, 2 + i } [2]. he distorted signal is received by D identical, omnidirectional antenna elements, and can be expressed as follows by using a complex baseband notation [3]: [ r ( t L r ( t] = r( t = s ( t; i f( t + n( t, ( D where superscript denotes the transpose. Note that both f ( t and n ( t are complex-valued. he waveform

sequence s (; t i, the fading vector f ( t, and the noise vector n ( t are assumed to be statistically independent. Let us for simplicity assume that the noise vector is composed of independent, stationary, zero-mean, white Gaussian processes. If we in addition assume that noise processes in different diversity lins are equally strong, the noise covariance-function matrix can be written as: { t t+ τ } = N0 E n( n ( I δ( τ, (2 where denotes the ermitian (complex conjugate transpose, N 0 is the double-sided power spectral density, I is the identity matrix, and δ ( denotes Dirac s delta function. Let the propagation geometry be restricted to the horizontal plane, and consider the Fraunhofer region (or the far-field; i.e., assume plane wavefronts. hese two assumptions are reasonable, provided the distance between transmitter and receiver is sufficiently large [4], [5]. If we let x l be the (Euclidean geometrical vector for receiver antenna l, and denotes a unit vector normal to the wavefront impinging from an angle θ, the array response vector (or steering vector can be written as [6]: { j j } = g L g, (3 a ( θ exp η [ η [ D where η is the circular wavenumber, and g denotes scalar product. Now, following [7], the frequency-flat fading vector is zero-mean and has the following covariance-function matrix: { f t f t+ τ } = A θ0 J0( π fmτ s E ( ( (, 2 / θ0 + = ( θ ( θ dθ J0 ( 2 π fmτ / s, 2 a a θ0 where the angle of arrival (AOA has been modeled according to a uniform distribution over [ θ 0, θ 0 + ]. Further, J 0 ( is the zero-order Bessel function of the first ind and f m denotes the maximum Doppler frequency shift normalized to the symbol duration s. Below, we will consider the two fading rates f m = 0.0 and f m = 0., which corresponds to slow and fast fading, respectively. From (4 it is seen that the fading induces infinite memory, i.e., the correlation extends over all time separations. he receiver first extracts N observables in each symbol interval up to, i.e., the received continuous-time process is represented by a sequence of N random variables. he obtained observation vector, say r, could be built up either by sampling or by an orthogonal expansion [0], [8]. In this paper, the observables will be derived by Nyquist sampling the (4 lowpass filtered antenna signals, where the cutoff frequency W of the (ideal lowpass pre-filter is chosen such that the loss of information becomes negligible. he observables are then used for computing soft information that is fed to the inner SISO (soft-input softoutput module [3]. In addition, interleaved extrinsic APPs generated from the outer of the two constituent codes are fed to the inner SISO and used as a priori probabilities. Only in the very first iteration, before APPs have been generated for the outer code, no a priori information is provided. In a similar manner, the inner SISO computes extrinsic APPs that are deinterleaved and used as a priori information in the outer SISO. Information obtained about the inner tail bits is superfluous and simply dropped before deinterleaving. Finally, after a certain number of iterations, the detector arrives at a decision for each information bit. AN APP ALGORIM USING FORWARD-ONLY PREDICION he generation of APPs for the outer trellis is a standard procedure. See [3] for any details. owever, the original BCJR algorithm [] was derived for a discrete memoryless channel (DMC, while the branch metric in our inner tree depends on previously derived observables. his implies the need for a different APP algorithm. he authors of [9] face a similar problem, but their proposal is less general. he APP for edge number l at time index can be written as: b ( i ( r = ( r ( ( l APP( e P e C P e P( e, l @, l, l, l, = C P( b ( i P r b ( i P e b ( i ( r b b (: i e, l b (: i e, l = C P ; ( i = C P ( i P ( i. where C is a normalization constant, b ( i is the bit sequence on hypothesis i, and b ( i : e, l denotes all bit sequences consistent with edge e, l. his quantity can be computed by the following recursive expressions, where the channel memory has been truncated to include only the L + most recent symbol intervals: (@ ( r ; b( r b r r b P i P i ( ( = P( ; ( i P( b ( i P ; ( i L L P( r ; b ( i P( b ( i P r r ; b ( i = P ( i γ (, i MS is not strictly bandlimited. (5 (6

and similarly: (@ ( r + ; b + ( r ; b( ( 2 2 + + r+ ; b+ ( r ; b ( + Pb ( + ( i P( r+ r ; b ( i ( 2 ; 2 ( + ; + r+ b+ r b ( + Pb ( + ( i P( r+ r+ L; b+ L( i P i P i i = P i i P i i + γ + = P ( i (. i Any realizable detector must truncate the (infinite channel memory and is inevitably suboptimal. owever, the correlation is fortunately a decreasing function of the time lag, which causes the approximation entailed in truncating the memory depth to become increasingly accurate as the truncation length is extended. As soon as the memory is truncated, the tree folds into a trellis. Note that the size of the trellis is determined by the truncation length, which is a design parameter. If the channel is truncated to include only the L + most recent symbol intervals, the modulation and the fading can thus described by a joint (time-invariant trellis having 2 L states. A trellis state is now defined as the L most recent data bits b ( fed to the MS modulator, while the L+ i (7 branches are labeled b Li (. Further, it is possible to terminate the inner trellis by appending L modulated zero tail bits to the transmitted waveform sequence, socalled post-interleaver flushing. From (5, the (approximated APP is finally obtained by summing over all possible sequences consistent with the considered edge. Since sequences that merge in the forward (or bacward direction will have identical branch metrics for all future (or past trellis sections, the summation in (5 can most efficiently be computed by means of the well-nown BCJR recursions []. his was also realized by Gertsman and Lodge [9]. We will now give an expression for the branch metric that proves to be equal to the classical MLSD (maximum-lielihood sequence detection metric [20]. Alternatively, using a set of linear prediction filters yield an equivalent result [20], [2], [22] (an approach taen in [9]. Let Z Li ( be a diagonal matrix comprising sam- ples of the unperturbed message sequence s L(; t i, and let F L denote the covariance matrix for samples of the fading process f ( t taen from time L up to, i.e., let F L be a discrete-time representation of (4. Since MS has constant envelope, i.e., ( i ( ( Es Z L i Z L = I, (8 s where E s is the symbol energy, it is possible to factorize the conditional covariance matrix for observables derived during symbol intervals L up to (the variance of each complex-valued noise sample is 2WN0 [23]: { ( i } s 0 ( L L L L C ( i @ E r r b ( = Z L( i F L + 2 WN I Z L( i. Es Note that the determinant of C Li ( does not depend on the hypothesis i (due to the constant envelope, which leads to the following form of the (conditional loglielihood function: { p( r r L b L i } ( ( i ( ( i log ; ( L L L L L L (9 r C ( r r C ( r. (0 It is interesting to note that the forward and bacward recursions use the same branch metric. Recently, socalled bi-directional algorithms have proved to achieve a performance gain for channels where the parametric uncertainty is modeled as a Gauss-Marov process [24]. NUMERICAL RESULS Error performance in terms of bit error rate ( versus signal-to-noise ratio (SNR has been obtained from Monte-Carlo simulations with a minimum bit error limit of 000 errors. Samples of the fading vector f ( t were generated by means of a method similar to the one described in [25] (the correlation coefficients were derived by Cholesy decomposition. Moreover, we investigated systems with a (7,5 convolutional encoder and a 024 o bit interleaver, assuming the mean AOA θ 0 = 45 and the azimuth spread 2 = 90 o. Figure shows the five first iterations for a twoelement array with antenna separation 0.5 carrier wavelengths, operating over a fading channel with bandwidth f m = 0.. Clearly, the receiver gains from increasing the exploited channel memory. Specifically, the iteration gain becomes more significant for large truncation lengths L, but recall from (0 that this performance improvement is achieved at the cost of considerably higher detection complexity. Figure 2 shows the three first iterations and assumes f m = 0.. hree different antenna constellation were investigated: a single antenna, a two-element array with antenna separation 0.5 carrier wavelengths, and a fourelement array with the elements equidistantly placed at a

circle with diameter 0.5 carrier wavelengths. Even though there is little to gain from running more than two iterations, there is obviously much to gain from spacediversity; the error performance of our proposed array receiver is superior to previously published results obtained for fast fading channels [20], [9]. An explanation to the moderate iteration gain is that only hard decisions are circulating between the inner and outer SISO modules after the second iteration. his might be improved by replacing MS be a CPM scheme with larger memory, thus obtaining a stronger inner code. In figure 3, the two fading bandwidths f m = 0.0 and f m = 0. have been investigated for a single-antenna receiver, as well as for a two-element array with antenna separation 0.5 wavelengths. When the fading rate is slow compared with the transmission rate, i.e., when f m = 0.0, iterating more than two times is seen to become somewhat more powerful. Figure : Error performance for various truncation lengths. 0 0 2 0 3 4 4.5 5 5.5 6 6.5 7 7.5 0 0 0 0 2 0 3 0 4 runcation Length L=2 L=4 L=6 D= (Single Antenna D=2 (2 Element Array D=4 (4 Element Array 0 5 0 2 4 6 8 0 2 Figure 2: Error performance for various antenna constellations. 0 0 0 0 2 0 3 0 4 0 5 0 2 4 6 8 0 2 Figure 3: Error performance for various fading bandwidths and antennas. CONCLUSIONS his paper has generalized the BCJR algorithm to be valid for fading channels. Moreover, numerical results indicated that iterative decoding offers a limited performance gain for fast fading channels, because most of the improvement is achieved in the second iteration. With a large truncation length L, the performance is significantly improved but the cost is unfortunately high decoding complexity. owever, a comparable or even superior gain is easily achieved by using additional receiver antenna elements. Stimulated by the good results in [24], the authors are currently woring on generalizing the bi-directional algorithms to handle the more general channel model that was used in this paper. REFERENCES D=, f m =0. D=, f m =0.0 D=2, f m =0. D=2, f m =0.0 [] L. R. Bahl, J. Coce, F. Jeline, and J. Raviv, Optimal decoding of linear codes for minimizing symbol error rate, IEEE rans. Inform. heory, vol. 20, pp. 284 287, Mar. 974. [2] C. Berrou, A. Glavieux, and P. hitimajshima, Near Shannon limit error-correcting coding and decoding: urbo codes, Proc. IEEE International Conf. Commun. (ICC 93, pp. 064 070, Geneva, May 993. [3] S. Benedetto, D. Divsalar, G. Montorsi, and F. Pollara, Serial concatenation of interleaved codes, IEEE rans. Inform. heory, vol. 44, no. 3, pp. 909 926, May 998. [4]. Aulin, CPM A Power and Bandwidth Efficient Digital Constant Envelope Modulation Scheme, Ph. D. dissertation, University of Lund, Sweden, Nov. 979.

[5] P. Moqvist and. Aulin, Power and bandwidth efficient serially concatenated CPM with iterative decoding, Proc. IEEE Global elecommun. Conf. (GLOBECOM 00, San Francisco, Nov./Dec. 2000. [6] C. Brutel and J. Boutros, Serial concatenation of interleaved convolutional codes and M-ary continuous phase modulations, Ann. élécommunications, vol. 54, no. 3/4, pp. 235 240, Mar./Apr. 999. [7] M. Peleg and S. Shamai (Shitz, Efficient Communication over memoryless Rayleigh fading channels with turbo coding/decoding, Proc. IEEE International Conf. Commun. (ICC 99, pp. 83 88, Vancouver, June 999. [8] E. Biglieri, J. Proais, and S. Shamai (Shitz, Fading Channels: Information-theoretic and communication aspects, IEEE rans. Inform. heory, vol. 44, no. 6, pp. 269 2692, Oct. 998. [9] I. D. Marsland and P.. Mathiopoulos, Multiple differential detection of parallel concatenated convolutional (turbo codes in correlated fast Rayleigh fading, IEEE J. Select. Areas Comm., vol. 6, no. 2, pp. 265 275, Feb. 998. [0] U. ansson and. Aulin, Aspects on single symbol signaling on the frequency flat Rayleigh fading channel, IEEE rans. Commun., vol. 47, no. 6, pp. 874 883, June 999. [] S. Pasupathy, Minimum shift eying. A spectrally efficient modulation, IEEE Commun. Mag., pp. 4 22, July 979. [2] B. Rimoldi, A decomposition approach to CPM, IEEE rans. Inform. heory, vol. 34, no. 2, pp. 260 270, Mar. 988. [3] J. G. Proais, Digital Communications, 3 rd ed., McGraw-ill, 995. [4] P. C. F. Eggers, Angular propagation descriptions relevant for base station adaptive antenna operation, luwer Personal Commun., vol., no., pp. 3 29, Oct. 999. [5] C. A. Balanis, Antenna heory: Analysis and Design, arper & Row, 982. [6] S. R. Saunders, Antenna and Propagation for Wireless Communication Systems, John Wiley & Sons, 999. [7] A. F. Naguib, Adaptive Antennas for CDMA Wireless Networs, Ph. D. dissertation, Stanford University, Aug. 996. [8] A. ansson and. Aulin, On antenna array receiver principles for space-time-selective Rayleigh fading channels, IEEE rans. Commun., vol. 48, no. 4, pp. 648 657, Apr. 2000. [9] M. Gertsman and J. Lodge, Symbol-by-symbol MAP demodulation of CPM and PS signals on Rayleigh flat-fading channels, IEEE rans. Commun., vol. 45, no. 7, pp. 788 799, July 997. [20] U. ansson, Efficient Digital Communication over the ime Continuous Rayleigh Fading Channel, Ph. D. dissertation, Chalmers Univ. ech., Dec. 997. [2] J.. Lodge and M. L. Moher, Maximum lielihood sequence estimation of CPM signals transmitted over Rayleigh flat-fading channels, IEEE rans. Commun., vol. 38, no. 6, pp. 787 794, June 990. [22] X. Yu and S. Pasupathy, Innovations-based MLSE for Rayleigh fading channels, IEEE rans. Commun., vol. 43, no. 2/3/4, pp. 534 544, Feb./Mar./Apr. 995. [23] C. Schlegel, rellis Coding, IEEE Press, 997. [24] A. Anastasopoulos,. Chugg, Adaptive soft-input soft-output algorithms for iterative detection with parametric uncertainty, IEEE rans. Commun., vol. 48, no. 0, pp. 638 649, Oct. 2000. [25] A. ansson and. Aulin, Generation of N correlated Rayleigh fading processes for the simulation of space-time-selective radio channels, Proc. European Wireless 99, Munich, pp. 269 272, Oct. 999.