PERFORMANCE OF DIFFERENTIALLY DETECTED GMSK IN A MOBILE RADIO CHANNEL by W I L L I A M CHI H U N G M O K B . E N G . A.C.G.I. (Electrical and Electronics Engineering), Imperial College, University of London, United Kingdom, 1995 A THESIS SUBMITTED IN PARTIAL F U L F I L L M E N T OF T H E REQUIREMENTS FOR T H E D E G R E E OF MASTER OF APPLIED SCIENCE in T H E FACULTY OF GRADUATE STUDffiS DEPARTMENT OF E L E C T R I C A L AND COMPUTER ENGINEERING We accept this thesis as conforming to the required standard T H E UNIVERSITY O F BRITISH C O L U M B I A April 1998 © William C . H. Mok, 1998 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of The University of British Columbia Vancouver, Canada DE-6 (2/88) Abstract As the demand for current cellular voice-oriented networks to include other services such as data and video applications is increasing rapidly, the reliability of information transmission over the mobile radio environment is important. Performance evaluation in terms of the bit error rate (BER) and block (packet) error rate (BKER) of modulation schemes in such an environment is of considerable research interest. In this thesis, we study the B E R as well as the B K E R performance of Gaussian-filtered Minimum Shift Keying (GMSK) modulation, currently used in the Global System for Mobile Communication (GSM) standard, in additive white Gaussian noise (AWGN) as well as flat and frequency selective Rayleigh fading. 1-bit and 2-bit differential detectors are used for demodulation. A detailed analysis is carried out to examine the IF filtering effect on the B E R performance over the fading channels. The B K E R performance is mainly studied through simulations. Since the premodulation Gaussian filter at the G M S K transmitter introduces significant Intersymbol Interference (ISI), we also examine the performance of a decision feedback (DF) scheme which partially removes the ISI, over such channels. The B E R results in both flat and frequency selective (FS) Rayleigh fading show that 2-bit conventional (with no decision feedback) G M S K performs slightly better than 1-bit conventional G M S K when E /N b Q < 25 dB. Above this value, the latter performs better and has a lower irreducible error. Results also show that 2-bit conventional G M S K is more sensitive to the IF filtering distortion than 1-bit conventional G M S K in Rayleigh fading. The percentage improvement in B K E R of 1-bit D F G M S K over 1-bit conventional G M S K is almost equal to that of the corresponding BERs. In the case of 2-bit ii Ill D F G M S K , the B K E R performance shows an increased improvement over 2-bit conventional G M S K when compared to the corresponding B E R performance, owing to increased error propagation. Since fading gives rise to bursty bit errors, the actual B K E R is lower than the predicted assuming independent bit errors. In all fading cases, 1-bit conventional G M S K has a better B K E R performance over 2-bit conventional G M S K when E /N b decision feedback, similar performance trends are observed. Q >25 dB. In the case of Table of Contents Abstract ii List of Tables vi List of Figures vii Acknowledgement x Chapter 1 Introduction 1 1.1 Motivation and Objectives 1 1.2 Outline of the Thesis 2 Background 4 2.1 Continuous Phase Modulation 7 2.2 Rayleigh Fading 14 2.3 Coherent Detection 17 2.4 Non-coherent Detection 18 2.5 Intersymbol Interference 19 B E R Performance of G M S K in Rayleigh Fading 20 BER Performance of 1-bit Detection in Flat Rayleigh fading channel 21 3.1.1 Unfiltered case 23 3.1.2 Unfiltered Case with Decision Feedback 35 3.1.3 Filtered Case 37 B E R Performance of 1-bit Detection in FS Rayleigh fading channel 46 3.2.1 Unfiltered case 46 3.2.2 Filtered case 52 Chapter 2 Chapter 3 3.1 3.2 iv V 3.3 BER Performance of 2-bit Detection in flat Rayleigh fading channel 58 3.3.1 Unfiltered case 59 3.3.2 Unfiltered Case with Decision Feedback 66 3.3.3 Filtered case 68 BER Performance of 2-bit Detection in FS Rayleigh fading channel 71 3.4.1 Unfiltered case 71 3.4.2 Filtered case 75 Results and Discussions 81 4.1 Simulation Model 81 4.2 Error Performance in a Non Fading Environment 85 4.3 Comparison of BER Performances between 'filtered' and 'unfiltered' cases in Rayleigh Fading 89 4.4 B E R Performance in Slow Rayleigh Fading 94 4.5 B K E R Performance in Rayleigh Fading 103 Conclusions 118 3.4 Chapter 4 Chapter 5 Glossary 121 Bibliography 126 Appendix A 130 Appendix B 132 Appendix C 135 Appendix D 141 List of Tables Table D . l Phase shifts 6 _ - (in degrees) for B T= 0.25 in 1-bit detection 141 Table D.2 Delayed phase shifts 0/_y(in degrees) for B T= 0.25 with x = 0.2T in 1-bit detection 141 Differential phase angles Att and Act^ (in degrees) of the one-bit detector corresponding to various input data combinations for B T= 0.25. The contributions of bj_ and bj are ignored 141 Table D.4 Phase shifts y _ • (in degrees) for B T = 0.25 in 2-bit detection 142 Table D.5 Delayed phase shift <jlr_ - (in degrees) for B T= 0.25 with x = 0.2T in 2-bit detection 142 Differential phase angles Afy and Aty (in degrees) of the one-bit detector corresponding to various input data combinations for B T= 0.25. The contributions of Cj_ and Cj are ignored 142 Table D.3 ; ; t t 5 ; t t 2 Table D.6 +2 [ t i t sl dt t 3 +2 vi List of Figures Figure 2.1 Block diagram of a digital communication system 4 Figure 2.2 General block diagram of continuous phase modulation 8 Figure 2.3 Information carrying phase pattern of (a) non-differentially encoded G M S K (b) differentially encoded G M S K 13 Figure 2.4 Response of Gaussian lowpass filter to a rectangular pulse of duration T. 13 Figure 2.5 A diagram to illustrate the differences between flat and frequency selective fading Two Ray Model: the block diagram of a frequency selective Rayleigh fading 16 channel 17 Figure 3.1 Block diagram of a 1-bit G M S K differential receiver 22 Figure 3.2 Block diagram of a 1-bit DF G M S K receiver 35 Figure 3.3 Phase-state diagram of 1-bit differential detection: (a) before feedback (b) after Figure 2.6 feedback , 36 Figure 3.4 Block diagram of 2-bit G M S K transmitter (for differential detection) 58 Figure 3.5 Figure 3.6 Block diagram of a 2-bit G M S K differential receiver Block diagram of a 2-bit DF G M S K receiver 59 66 Figure 3.7 Phase-state diagram of 2-bit differential detection: (a) before feedback (b) after feedback 67 Figure 4.1 BER of 1 and 2-bit conventional G M S K in a non fading environment with AWGN. For 2-bit 'detection', 'thres = 0' denotes the case where a decision threshold of zero is used 86 Figure 4.2 B K E R of 1-bit conventional G M S K in a non fading environment with AWGN...87 Figure 4.3 B K E R of 2-bit conventional G M S K in a non fading environment with AWGN. Note that the threshold is optimized for every E /N 88 b Figure 4.4 Q C D F of the number of bit errors in a block at = 16 dB in a non fading environment for (a) 1-bit conventional G M S K (b) 2-bit conventional G M S K with optimum threshold 89 vii Figure 4.5 Average BER as a function of f T for 'filtered' and 'unfiltered' cases of 1 and 2bit conventional G M S K in (a) flat fading (b) FS fading with CDR = 20 dB, x = 0.27 91 Figure 4.6 Average BER as a function of f T for 'filtered' and 'unfiltered' cases of 1 and 2bit conventional G M S K in FS fading with (a) CDR = 20 dB, x = 0.47 (b) CDR = 10 dB, x = 0.27" (c) CDR = 10 dB, x = OAT 92 Figure 4.7 Average BER as a function of CDR for 'filtered' and 'unfiltered' cases of 1 and 2bit conventional G M S K in FS fading for f T = 0.002 with (a) x = 0.2T (b) x = 0.47" 93 D D D Figure 4.8 Average BER as a function of E /N inflatfading Figure 4.9 Average B E R as a function of E /N x = 0.27" in FS fading with CDR = 20 dB, Figure 4.10 Figure 4.11 Figure 4.12 b 0 b a 98 Average BER as a function of E /N x = 0.47" in FS fading with CDR = 20 dB, Average BER as a function of E /N x = 0.27 in FS fading with CDR = 10 dB, b 0 b 95 99 Q 100 Average BER as a function of E /N in FS fading with CDR = 10 dB, x = 0.27. Note that d is chosen to be 0.5, instead of 1 for the simulated 1-bit DF G M S K curve and d is chosen to be 1.5, instead of 2 for the simulated 2-bit DF G M S K curve 101 b a x 2 Figure 4.13 Average BER as a function of E /N x = 0.47 b in FS fading with CDR = 10 dB, 0 102 Figure 4.14 Average B K E R as a function of E /N for (a) 1-bit conventional G M S K and 1bit D F G M S K in flat fading (b) 2-bit conventional G M S K and 2-bit D F G M S K in flat fading (c) 1-bit conventional G M S K and 1-bit D F G M S K in FS fading with CDR = 20 dB, x = 0.27 (d) 2-bit conventional G M S K and 2-bit D F G M S K in FS fading with CDR = 20 dB x = 0.27 104 Figure 4.15 Average B K E R as a function of E /N for (a) 1-bit conventional G M S K and 1bit D F G M S K in FS fading with CDR = 20 dB, x = 0.47 (b) 2-bit conventional G M S K and 2-bit D F G M S K in FS fading with CDR = 20 dB, x = 0.47 (c) 1bit conventional G M S K and 1-bit D F G M S K in FS fading with CDR = 10 dB, x = 0.27 (d) 2-bit conventional G M S K and 2-bit D F G M S K in FS fading with C D R = lOdB, x = 0.27 106 b b viii 0 Q Figure 4.16 Average B K E R as a function of E /N for (a) 1-bit conventional G M S K and 1bit DF G M S K in FS fading with CDR = 10 dB, x = 0.47/ (b) 2-bit conventional G M S K and 2-bit DF G M S K in FS fading with CDR = 10 dB,x = OAT 108 Figure 4.17 Average B K E R as a function of E /N for 1-bit conventional G M S K and 2-bit conventional G M S K in (a) flat fading (b) FS fading with CDR = 20 dB, x = 0.27/ (c) FS fading with CDR = 20 dB, x = OAT (d) FS fading with CDR = 10 dB, x = 0.27/ 110 Figure 4.18 Average B K E R as a function of E /N for 1-bit conventional G M S K and 2-bit conventional G M S K in FS fading with CDR = 10 dB, x = OAT Ill Figure 4.19 C D F of the number of bit errors in a block at E /N = 40 dB for (a) 1-bit conventional G M S K and 1-bit D F G M S K in flat fading (b) 2-bit conventional G M S K and 2-bit DF G M S K in flat fading 112 Figure 4.20 C D F of the number of bit errors in a block at E /N = 20 dB for 2-bit conventional G M S K and 2-bit DF G M S K in flat fading b b b a 0 Q b Q b 0 113 Figure 4.21 C D F of the number of bits errors in a block at E /N conventional G M S K in flat fading Figure 4.22 C D F of the number of bit errors in a block at E /N = 40 dB for (a) 1-bit conventional G M S K and 1-bit DF G M S K in FS fading with CDR = 20 dB, x = 0.27 (b) 2-bit conventional G M S K and 2-bit D F G M S K in FS fading with CDR = 20 dB, x = 0.27/ (c) 1-bit conventional G M S K and 1-bit D F G M S K in FS fading with CDR = 20 dB, x = 0.4T (d) 2-bit conventional G M S K and 2-bit DF G M S K in FS fading with CDR = 20 dB, x = 0.47/ 114 Figure 4.23 C D F of the number of bit errors in a block at E /N = 40 dB for (a) 1-bit conventional G M S K and 1-bit DF G M S K in FS fading with CDR = 10 dB, x = 0.27/ (b) 2-bit conventional G M S K and 2-bit D F G M S K in FS fading with CDR = 10 dB, x = 0.27/ (c) 1-bit conventional G M S K and 1-bit D F G M S K in FS fading with CDR = 10 dB, x = 0.47/ (d) 2-bit conventional G M S K and 2-bit D F G M S K in FS fading with CDR = 10 dB, x = 0.47/ 115 Figure 4.24 C D F of the number of bit errors in a block at E /N = 40 dB for (a) 1-bit and 2bit conventional G M S K in FS fading with CDR = 20 dB, x = 0.2T (b) 1-bit and 2-bit conventional G M S K in FS fading with CDR = 20 dB, x = 0.47/ (c) 1-bit and 2-bit conventional G M S K in FS fading with CDR = 10 dB, x = 0.27/ (d) 1-bit and 2-bit conventional G M S K in FS fading with CDR = 10 dB, x = 0.47/ 116 b b b b ix a = 40 dB for 1-bit and 2-bit 113 0 0 0 Acknowledgment I would like to express my gratitude to my research supervisor, Dr. Cyril Leung, for his continuous guidance and valuable advice throughout the course of the research work. Also, the sincere assistance from my colleagues is much appreciated. This work was partially supported by NSERC grant OGP0001731. Finally, I would like to thank my parents and my sister for their constant encouragement and support throughout my graduate study at U B C . Chapter 1 Introduction Today, the area of wireless communications is experiencing unprecedented growth, with the development of a cellular integrated services network, known as the Universal Mobile Telecommunication System (UMTS) [1] being under way. Currently, the Global System for Mobile Communication (GSM) forms the standard for most mobile voice networks worldwide [2]. The introduction of General Packet Radio Service (GPRS) [3], intended to support frequent 1 non continuous transmission of bursty packet data up to data rate of 200 kbits/s in the G S M networks, will provide a proving ground for the integrated services to be supported in future U M T S networks. Therefore, G S M is now evolving from a voice-oriented network that generally supports circuit switching, to include a packet switching platform which improves bandwidth utilization for packet data transmission. Gaussian-filtered Minimum Shift Keying (GMSK) is the modulation scheme currently used in the G S M standard. It belongs to a modulation class, known as Continuous Phase Modulation (CPM) which has been closely examined due to its improved spectral characteristics and power efficiency relative to that of Phase shift keying (PSK) and Frequency modulation (FM). 1.1 Motivation and Objectives In connection with the current development trend in the G S M standard, it is worthwhile to study the bit error rate (BER) as well as the block (packet) error rate (BKER) performances of G M S K in a mobile radio channel. Differential detection is to be used, because it eliminates the need for precise regeneration of the reference carrier, difficult to achieve in fading conditions, 1 Before the introduction of GPRS, G S M is already supporting packet switching using the Short Message Service (SMS) intended to provide infrequent, non continuous transmission of very short data bursts of length up to 160 characters (seven bit data). 1 Chapter 1 Introduction 2 otherwise. The BER performance of G M S K in various channels has been investigated in [5], [9]-[20], [22],[23]. However, the B K E R performance has rarely been studied. In this thesis, a simulation model of G M S K employing 1-bit and 2-bit differential detection is developed to study the BER and B K E R performance in non fading additive white Gaussian noise (AWGN) channel, flat and frequency selective (FS) Rayleigh fading channels. A number of papers which analyse differential detection of G M S K in flat and FS Rayleigh fading channels, assume that the intermediatefrequency (IF) filter at the receiver has a negligible effect on the signal spectrum but little justification has been given [9],[10],[13],[19],[21]. Therefore, a detailed analysis is also presented to compare the B E R performance of G M S K when the IF filtering effect is taken into account and when it is not. The numerical results from the analysis are verified by comparison with the simulation results. Since the premodulation Gaussian filter at the transmitter for G M S K introduces significant ISI, a decision feedback scheme was proposed in [22] to partially remove the effect of destructive ISI. The performance of this scheme was simulated in [22] for a non fading environment and it was analyzed in [13], [23] and simulated in [16] for flat fading channels; but not in FS fading channels. In this thesis, the BER and B K E R performances of this scheme in flat and FS Rayleigh fading channels are investigated. In particular, the B E R performances over such channels will be analysed and compared with the simulation results. The B K E R performance is used to look at the bit error distribution under this scheme. 1.2 Outline of the Thesis This thesis is organized as follows. In Chapter 2, some background material including the mathematical modelling of mobile radio Chapter 1 Introduction 3 channels is reviewed. In Chapter 3, we analyze the B E R performance of differentially detected G M S K with and without decision feedback in flat and frequency selective Rayleigh fading channels. In Chapter 4, the BER and B K E R performances in AWGN, flat and FS Rayleigh fading channels are presented. In particular, the numerical B E R results from chapter 3 are compared with the simulation results. In Chapter 5, the conclusion and some recommendations for future work are given. Chapter 2 Background Figure 2.1 shows the functional blocks of a digital communication system [6]. The input source may be either an analog or a digital signal. The source encoder maps the analog signal into a sequence of binary digits with as little redundancy as possible. The channel encoder encodes the binary information sequence so as to introduce some redundancy into the information. The added redundancy can assist the receiver to reconstruct the transmitted data from the received signal which is corrupted by noise and distortions. The output of the channel encoder is modulated onto a carrier. The focus in this thesis is on the modulator/demodulator section of the system. digital signal source Source Encoder transmitted signal Modulator Channel Encoder Transmitter earner (^Demodulator Transmission Channel Channel Decoder received signal Noise, interference and distortion Figure 2.1 Receiver Clock (synch.) carrier reference (coherent) Block diagram of a digital communication system 4 Source Decoder output digital signal Chapter 2 Background 5 There are a number of considerations that enter into the choice of a modulation technique for use in a wireless application. In a communication medium or channel with additive white Gaussian noise, Shannon derived a formula for the capacity, C (bits/s), as a function of the bandwidth of the channel, W (Hz), and the average received signal to noise power ratio, SNR [6] C = Wlog (l + S N R ) 2 (2.1) where P 2 is the transmitted signal power, E s b (postdetection) noise power, R b is the transmitted energy per bit, G is the received n is the transmission rate in bits/s and N 0 is the one-sided power spectral density (PSD) of AWGN. The significance of the channel capacity formula given above is that it puts an upper limit on R for reliable communication over a noisy channel. Specifically, one can transmit informab tion, with as small an error probability as desired, as long as R is smaller than C. Using (2.1), we b can plot R /W b versus E /N , b 0 giving the bandwidth efficiency plane [7]. From the plane, we can see that there is a trade-off to be made between the transmitted power P and the channel s bandwidth W in maximizing the channel capacity and achieving a specified level of probability of bit error (BER). In designing a digital communication system, the objectives are to minimize the 2 P is denned as the desired (average) signal power at the point before AWGN is added. In this case, it is equal to the transmitted signal power. Chapter 2 Background 6 B E R and the parameters P and W; but also to maximize R . However, these objectives are in s b conflict with each other. For example, if we reduce P , with R s b and W remaining constant in the system, then the B E R is degraded. The ratio Rj/W is a measure of the bandwidth efficiency. For example, m-ary phase shift keying (MPSK) requires a bandwidth at IF (double-sided bandwidth) equal to the symbol rate, R [7]: s W = R (2.2) s where R v 1 = R b log m . Therefore, the bandwidth efficiency of MPSK modulated signals can be ex- 2 pressed as R — = log m. (2.3) 2 On the other hand, non-coherent m-ary orthogonal frequency shift keying (MFSK) modulation requires a bandwidth equal to mR corresponding to the m different orthogonal waveforms. s Hence, the bandwidth efficiency of M F S K modulation is equal to R h \og m 0 w Vv = —m <-> 2 4 At a bit rate of 9600 bits/s and a symbol rate of 2400 symbols/s, M P S K has a bandwidth efficiency of 2 (hence m=4) and requires an E /N b whereas 4-FSK requires an E /N b 0.5 [7]. a of 17.5 dB to achieve a B E R of 10~ , 5 a of only 8.1 dB but can only achieve a bandwidth efficiency of Chapter 2 Background 7 2.1 Continuous Phase Modulation Most mobile radio products are designed with Class-C power amplifiers, which offer the highest power efficiency among the common types of power amplifiers. However, this class of amplifiers are highly nonlinear, which requires the signal to be amplified to have a constant envelope. Any amplitude fluctuations in the input signal will result in a spectral widening of the output signal, which gives rise to adjacent channel interference (ACI) as explained in the next section [2]. Continuous phase modulation (CPM) [8] is a class of modulation schemes that can achieve a good compromise between power and bandwidth efficiency while maintaining a specified level of B E R at the expense of only a reasonable increase in system design complexity. In C P M , the phase is constrained to be continuous, making memory being present in the input to the phase or frequency modulator. The RF signal after modulation has a constant envelope which contributes to a couple of advantages. Firstly, constant envelope signals are less sensitive to amplifier nonlinearities due to the lack offluctuationsin the signal amplitude, resulting in a more compact power spectrum. Sidelobe spreading, as compared with that of non-constant envelope signals is reduced hence improving the bandwidth efficiency . 3 Secondly, they are more resistant to adjacent channel interference (ACI), interference due to signals of adjacent channels. This is due to their compact power spectrum so that more of these energy is contained in the channel bandwidth [6],[8]. In C P M , the transmitted signal s(t) is given as If we define bandwidth as the frequency band that contains 99.9% of the signal power, then the smaller the normalized bandwidth (bandwidth-bit duration product) of a modulation scheme, the higher the bandwidth efficiency it can achieve. Chapter 2 Background s(t) = 2£„ l^ os(2nf t C c + Q>(t,b)) (2.1.1) bjgiy-jT) \dv (2.1.2) where ' f $(t;b) = 2nh\ ]T (2.1.3) h = 2f T d where E is the transmitted energy per symbol, T is the symbol duration, h is the modulation s index,/^ is the peak frequency deviation, and g(t) will be described later. Premodulation Filter d(t) FM Modulator HJf) s(t) NRZ Sequence 2nh Figure 2.2 General block diagram of continuous phase modulation Assuming the transmitted power P = E /T = 0.5 , (2.1.1) can be written as s s(t) = cos(2nf t + §(t;b)) c (2.1.4) The transmitted information is contained in the phase which is used to modulate a carrier Chapter 2 Background 9 in such a way that its phase is continuous. In Figure 2.2, {bj} denotes the non-return-to-zero (NRZ) sequence of amplitudes obtained by mapping fc-bit blocks of binary digits from the information sequences {a.j} into the amplitude levels ± 1 , ± 3 , ... ± (ra - 1) where ra = 2 . After k the premodulation filter with frequency response H (f), t the output signal d(t) can be expressed as oo d(t) = X bjg(t-jT) (2.1.5) j = -°° with g(t) = h {t) ® p(t), where ® denotes convolution, p(t) is a rectangular pulse of duration T t and unity amplitude, h (t) is the impulse response of the premodulation filter and g(t) is the t response of the transmitted rectangular pulse to the pre-modulationfilter.If g(t) = 0 for t > T, the C P M signal is called full response C P M . If g(t) ^ 0 for t >T, then the C P M signal is called partial response C P M [8]. From the definition of F M modulation, one can observe that the phase response described in (2.1.2) is a function of the integral of the pulse g(t). The shape of g(t) determines the smoothness of the transmitted information carrying phase. The rate of change of the phase is proportional to the parameter h, which is the modulation index. The pulse g(t) is normalized such that oo J g(t)dt = 0.5 [8]. —oo A variety of C P M signals can be obtained by varying the three parameters, h, g(t) and w. A few popular schemes include Tamed Frequency Modulation (TFM), Raised Cosine (L-RC), Chapter 2 Background 10 Rectangular (L-REC) and Gaussian Minimum Shift Keying (GMSK). M S K is a member of the L - R E C scheme which has a rectangular pulse as its pulse response: 0<t<LT 1 g(t) = 2LT 0 . (2.1.6) otherwise In M S K , L is 1, m = 2 and h is 0.5. It should also be pointed out that the premodulation filter is not required in the implementation of M S K . Following from (2.1.2), the phase of the carrier in the interval nT < t < (n + 1)7/ [6] is W;b) n= \% £ „ 1 bj+ ; =— 1 , Kb q(t-nT) n (2.1.7) ft-nT where ne» = f 1 X ; *y ^(r) = \g(%)dx. Substituting (2.1.7) into (2.1.4) gives (2-1-8) = -oo (2.1.9) Chapter 2 Background 11 1 , ft-nT + Q + ^nb ^ s(t) = cos(2nf t c n n " T = cos {**{fc + ^fK} ~ \™b (2.1.10) + G„) n nT < t < ( » + 1)7" The equation above indicates that M S K signal is a sinusoid having one of the two possible frequencies in the interval nT <t<(n + 1)T .If we define these frequencies as fi=fc-jf fl (2.1.H) = fc + 4T then equation (2.1.10) can be written as s (t) t = cos(27u/.r + e + i m t ( - l ) ' j , _ 1 n The frequency separation A / = f -f\ 2 nT<t<(n+\)T, = i = 1,2. (2.1.12) is the minimum frequency separation that is necessary to ensure orthogonality of the signals s^t) and s (t) over a signalling interval of 2 length T and explains the name Minimum Shift Keying. One of the advantages of M S K is that it does not introduce any intersymbol interference (ISI). The transmitted pulse corresponding to each symbol is confined within the bit duration, so that there is no adjacent interference between symbols. M S K allows a simple implementation as in the F M form depicted in Figure 2.2. It can also be realized in the form of 4-phase PSK as it is the only linear modulation in the C P M family [2]. Gaussian-filtered Minimum Shift Keying (GMSK) was invented in an attempt to improve the M S K power spectrum [27],[28]. In GMSK, the power spectrum is made more compact, with Chapter 2 Background 12 spectral sidelobes further reduced by applying a lowpass filter (i.e. the premodulation filter as shown in Figure 2.2) to the data stream prior to modulation in its F M form. The constant envelope property of C P M is preserved as filtering is performed before modulation. The premodulation filter chosen is a Gaussian lowpass filter because it satisfies certain properties suitable for reshaping the output power spectrum, hence the name GMSK. The impulse response of the premodulation Gaussian filter, h (t) is given by t (2.1.13) where K is an arbitrary constant (the choice of K is elaborated in Chapter 3) and B is the 3-dB t bandwidth of the filter. Figure 2.3 shows the information carrying phase §(t;b) of M S K and G M S K for a particu- lar data sequence. Note that in differentially encoded G M S K , the data sequence {bj} shown in Figure 2.2, is a differentially encoded version of the actual transmitted information sequence. Note that the phase changes for G M S K are smooth due to pulse shaping by premodulation filtering, whereas M S K has linear, sharp phase changes during data transitions. To illustrate the bandwidth efficiency of C P M signals in general, we can define the bandwidth as the frequency band that contains 99.9% of the signal power. The 3-dB bandwidth of the premodulation Gaussian lowpass filter is denoted by B . For M S K where B T = °°, the normalized bandwidth is 2.76 and t t for G M S K with B T= 0.25, it is approximately 1 [27]. Hence, G M S K (with B T < o o ) achieves a t t higher bandwidth efficiency than M S K at the expense of the introduction of ISI. In general, as B T t decreases, the bandwidth efficiency of the signal increases but the spreading of the pulse also Chapter 2 Background 13 becomes more severe, causing more ISI. (a) Figure 2.3 (b) Information carrying phase pattern of (a) non-differentially encoded GMSK (b) differentially encoded GMSK Figure 2.4 illustrates the response of the Gaussian filter to a rectangular pulse of period T and the effect of variation of B T on pulse spreading. Both M S K and G M S K can be coherently or t non-coherently detected, as further described in Sections 2.3 and 2.4. 1 1 1 1 1 1i 1 1 B T = ~(MSK) ( B T = 0.5 t - A h B T = 0.25 ( jX 0 Figure 2.4 1 2 3 4 5 y 6 B T = 0 1 t 7 8 Response of Gaussian lowpass filter to a rectangular pulse of duration T. 9 Chapter 2 Background 14 2.2 Rayleigh Fading In most mobile radio channels, the transmitted wave travels to the receiver via multiple paths, reflected from objects such as buildings or scattered from objects close to the receiver. The channel in which the transmitted signal can propagate via multiple paths is called a multipath channel. Associated with each path, there exists an attenuation factor and a propagation delay; the former accounts for the amplitude variations in each multipath component and the latter accounts for the phase variations. Due to the constructive and destructive interferences of these components in the vincinity of the receiver, there will be amplitude and phase variations in the received signal. Fading occurs when the receiver moves through an area of destructive interference [2], [6]. A set of channel correlation functions and power spectra that define the characteristics of a fading multipath channel can be developed. The characterization involves the following parameters: the time variations of the channel can be measured by the Doppler spread and the coherence time. The Doppler spread f D is the width of the Doppler spectrum, which represents the strength of the Doppler shift at different frequencies, caused by the movement of the mobile. The coherence time t is approximately equal to the reciprocal off c D and is a measure of the time interval over which the channel conditions remain approximately constant. The frequency variations of the channel can be measured by the coherence bandwidth and the multipath delay spread. The coherence bandwidth Af c is a measure of the bandwidth over which the signals spaced in frequency, will be similarly affected by the channel. The reciprocal of Af c is approxi- mately equal to the multipath delay spread , x , which is a measure of the time dispersion of a 4 d 4 Another measurement of delay spread is the root mean squared (rms) delay spread [42]. Chapter! Background 15 transmitted signal caused by the channel [6], [29], [30]. Fading channels are characterized asflatif all frequencies of the signalling bandwidth are similarly affected by the channel; otherwise they are frequency selective, where different frequencies are subjected to different fading. Therefore from the above definition, if the transmitted signal bandwidth is small in comparison with A / , c the channel is said to be flat or frequency non-selective; otherwise it is said to be frequency selective (FS). The sum of these multipath components or phasors, is a complex Gaussian process. The magnitude of the process has a Rayleigh distribution. It will change into a Ricean distribution in the presence of a line-of-sight component [2]. In a Rayleigh fading channel, the transmitted signal is randomly modulated with a complex Gaussian process denoted as Z (t), where s Z (t) s = x (t) + jy (t) s (2.2.1) s and x (t), y (t) are independent, zero mean Gaussian lowpass processes. s s A model for mobile radio channels has been proposed by Clarke [31]. The Doppler power spectrum derived from this model is given by \f\<f. D 1/1 >/D 0 where f D (2.2.2) is the maximum Doppler frequency (the Doppler spread) given by • - ZD ~ X- ll c V (2.2.3) Chapter 2 Background 16 In Equation (2.2.3), v denotes the speed of the mobile unit in metres/second (m/s),/ is the carrier c frequency in Hz and c is the speed of light which is equal to 3 x 1 0 m/s. The implementation of 8 the Rayleigh fading channel in the software simulation will be further described in Chapter 4. Since the multipath delay spread is approximately equal to the reciprocal of A / , we can c also categorize a channel as flat or frequency selective according to the delay spread. Figure 2.5 illustrates the differences between flat and FS fading. In general, multipath signal components of which their path delays are closely distributed about a mean delay would be regarded as one Rayleigh faded signal (via the Central Limit Theorem). The term 'closely distributed' means the difference in path delays are small relative to a symbol duration. The mobile channel in which only one such path exists is known as a flat Rayleigh fading channel (path 1). Signal components with path delays closely distributed about a mean delay different from that of path 1, would be regarded as a second delayed Rayleigh faded path, represented by path 2 in Figure 2.5. If the \ Receiver Figure 2.5 A diagram to illustrate the differences between fiat and frequency selective fading Chapter 2 Background 17 difference in the mean path delays, T , is significant relative to a symbol duration, then there exists frequency selectivity. A frequency selective channel model contains multiple paths of which their mean path delay differences are significant with respect to a symbol duration. In this thesis, we will use the two ray model (represented by path 1 and path 2 in Figure 2.5) which is simple and adequate enough to simulate frequency selectivity in mobile radio channels. The model can be implemented as shown in Figure 2.6. It has been reported in [38] that channels with rms delay spread less than about 0.1 of a symbol duration are considered flat; otherwise they are considered as frequency selective. PATH 1 AWGN signal at output of transmitter Rayleigh faded signal Z (t) and Zj(t) s are independent complex Gaussian processes PATH 2 Figure 2.6 Two Ray Model: the block diagram of a frequency selective Rayleigh fading channel 2.3 Coherent Detection In coherent detection systems, the receiver requires an estimate of the carrier frequency and phase in order to recover the transmitted data in the received signal. A majority of methods in combating fading in coherent systems involve the transmission of pilot tones [33],[34] or pilot sequences [35],[36]. The pilot scheme involves the transmission of a pilot tone at a convenient frequency in the data spectrum. The pilot tone is extracted at the receiver and the channel impairments can then be deduced from the tone. The information from the channel can assist in the Chapter 2 Background 18 reconstruction of a coherent reference signal. The disadvantage of this scheme is that it resembles a double sideband modulation scheme which is bandwidth inefficient because the sidebands need to have zero response around dc, to allow bandwidth for pilot transmission. A n improvement to this scheme is proposed in [34]. 2.4 Non-coherent Detection In contrast to coherent systems, non-coherent detection does not need to generate a carrier reference at the receiver, hence allowing simple implementation. Non-coherent detection can be subdivided into limiter/discriminator [4],[11],[15],[17] and differential detection [9],[10],[13],[16],[18]. Limiter/discriminator consists of a hard limiter to remove the amplitude variations of the received wave, followed by a bandpass filter [26]. The resulting wave is a constant-envelope sinusoidal and is demodulated by a discriminator. However, this detection scheme suffers from the F M threshold effect when the signal carrier to noise ratio falls below a certain level. In differential detection, information, instead of being transmitted using absolute carrier phases, is encoded into carrier phase differences [24]. Information is therefore recovered at the receiver by obtaining the difference in phase of received signal at current time and at some time instant in the past, usually at multiples of the symbol period. Since information is transmitted and detected through phase differences of the carrier, the absolute carrier phase is not important and so the need for carrier recovery is eliminated. Differential detection does not suffer from the threshold effect and is affected little by severe channel delay distortions [21]. Its performance in A W G N is worse than that of coherent detection, but it has an improved performance in fading channels [9],[10],[27]. This is because it is difficult to precisely regenerate the reference carrier in coherent detection under fading conditions. Moreover, differential detection improves performance by cancelling the phase distortion between adjacent symbols during phase subtraction. In the Chapter 2 Background 19 case of G M S K , the inherent ISI due to pulse shaping can be partially eliminated using decision feedback [22]. The differential detection of G M S K is described in Chapter 3. 2.5 Intersymbol Interference Intersymbol interference (ISI) introduced by the premodulation filter in the transmitted G M S K signals was described in Section 2.1 and in this section, we will describe the ISI induced by the channel. In any communication systems, it is desirable to transmit information at a rate as high as possible while maintaining the reliability of transmission. Due to non ideal channel frequency response characteristics, a succession of pulses transmitted at rates on the order of the channel bandwidth W will be distorted and they overlap each other so that each of the pulses can no longer be distinguished independently. The overlapping affects decision making at sampling instants since there is interference contributed from other adjacent symbols, giving rise to ISI. A signal pulse can be designed that allows us to transmit at symbol rates equal to or smaller than 2W the Nyquist rate, with zero or a controlled amount of ISI. Non-constant envelope digital modulation schemes make use of filters to satisfy Nyquist's criteria for zero ISI or duobinary coding which introduces ISI in a controlled manner in order to achieve a transmission rate equal to 2W [6],[25], assuming that the channel's response characteristics are known a priori. If the channel conditions are not known a priori, we have to employ equalization methods to eliminate or reduce ISI [6]. In the case of GMSK, it is not possible to use filters for pulse shaping to avoid ISI. This is because the pulse shaping in G M S K is performed using the premodulation Gaussian lowpass filter which can improve on the bandwidth efficiency as well as maintaining the constant envelope property of G M S K signals. In this thesis, we will only look at the effect of ISI in G M S K signals caused by the premodulation filter and consider the use of decision feedback [22] to combat ISI. Chapter 3 BER Performance of GMSK in Rayleigh Fading A careful analysis of BER performance of a digital communication system generally takes into account the following factors which cause signal distortions during transmission: additive white Gaussian noise, the statistical characteristics of the flat fading channel, the frequency selectivity in the channel, Doppler spread in the carrier frequency due to the motion of the mobile unit, intersymbol interference due to the premodulation filter, intersymbol interference due to bandlimiting by the IF filter, the co-channel interference from a user in another cell and the adjacent channel interference which arises between channels whose frequency bands are in proximity. Some of these factors have been described in the previous chapter. In this chapter, we will analyse the IFfilteringeffect on the BER performances of 1 and 2bit differentially detected G M S K (hereafter called 1 and 2-bit G M S K ) in flat and FS Rayleigh fading channels. We will also look at the BER performances when a decision feedback scheme is employed on such channels. As mentioned in Chapter 1, the B E R performances of 1 and 2-bit G M S K have been analysed in [9], [10], [13] under the assumption that the IF filtering effect is negligible (the 'unfiltered' case) in flat Rayleigh fading channels but the IFfilteringeffect on the BER performances has not been analysed. There are also papers on the B E R performances of 1 and 2-bit G M S K in FS Rayleigh fading channels [14], [18], [20] where the 'unfiltered' case is assumed in [14], and the filtering effect has been taken into account (the 'filtered' case) in [18] and [20]. However, there has not been a direct comparison between the B E R performances for the 'unfiltered' and the 'filtered' case. 20 Chapter 3 BER Performance of GMSK in Rayleigh Fading 21 In [9] and [10], the modulation index, h, of the G M S K signal is chosen so as to minimize the probability of error. In this thesis, we will use h = 0.5 which allows G M S K to be demodulated using coherent detection as that employed in MSK. Hence, choosing such a value of h will allow us to compare the performance of differential detection with results from other papers in which coherent detection is looked at. We will use the approach developed in [10] which deals with the performance of 2-bit G M S K in flat Rayleigh fading channel to evaluate the B E R performances of 1 and 2-bit G M S K in both flat and FS Rayleigh fading channels. This approach eliminates the need for computing a double integral which appears in [20],[21]. In the case of a non-fading channel, A W G N with two-sided power spectral density (PSD) N(/2, is simply added to the transmitted G M S K signal before the input of the receiver. The BER performance of 1 and 2-bit differentially detected G M S K in A W G N channel is analysed in [5]. The 2-bit detection performs better than the 1-bit detection and the performance advantage is greater for lower values of the bandwidth-time product, B T of the Gaussian premodulation filter. t The improvement in the BER performance for the 2-bit case is due to the asymmetric eye opening of the 2-bit detector eye patterns, enabling one to bias the detection threshold towards the larger opening of the eye, as compared to that of the 1-bit case. This makes 2-bit detection less sensitive to errors caused by AWGN. 3.1 B E R Performance of 1-bit detection in Flat Rayleigh fading channel Figure 2.2 in Chapter 2 illustrates a G M S K transmitter, consisting of a premodulation Gaussian lowpass filter and a F M modulator with modulation index, h = 0.5 giving a peak frequency deviation/^ equal to 1/(47/). The output of the premodulation filter, d(t) and the Chapter 3 BER Performance of GMSK in Rayleigh Fading 22 transmitted signal at the output of the modulator, s(t) are given by (2.1.5) and (2.1.4) respectively. Figure 3.1 shows the receiver for the 1-bit detection case. The output of the 1-bit detector V'(t) is o obtained from multiplying the output of the IF filter, r(t) by a 1-bit period delayed, 90 phase shifted version of itself. The product is then passed through a lowpass filter to produce which is sampled at t = IT to give V(IT), V(t), where / is the Ith sampling interval: le {...,-1,0, 1,2, ...} . A hard decision is then made on V(IT) to decide whether a "1" or a "0" was sent. In an ideal channel (no fading and no noise), the receiver output is given by (3.1.23). The operation of the 1-bit differential detector will be described mathematically in Section 3.1.1. Detector Ht) 90° Delay T n(t) phase shift V(IT) Sample Sj( t): signal after fading with average power o 2 Hard Decision "0" or "V n(t): AWGN which produces narrowband noise with 2 average noise power <s at the IFfilteroutput Figure 3.1 Block diagram of a 1-bit GMSK differential receiver Chapter 3 BER Performance of GMSK in Rayleigh Fading 23 3.1.1 Unfiltered case After propagation through the flat fading mobile radio channel, the signal after fading, 2 sj(t), with average power <J , is given by S s (t) f = Re[(x (t)+jy (t))-&xpM (t) s s + (O t)] s c = / J « [ { c o s ( c o r + (|),(0) + j s i n ( ( D r + <|>,(0)} • {x (t) + jy (t)}] = x,(f)cos(co f + ^(r)) - y (t) sin(co r + 0,(0) c c c s s (3.1.1) s c where x (t)and y (t) are independent stationary zero mean random Gaussian lowpass processes s s that define the Rayleigh fading channel and § (t) is the information carrying phase, as described s in Chapter 2. We first consider the 'unfiltered' case in which the filter bandlimits the noise but does not distort the signal. Assuming the noise is Gaussian and broadband prior to IF filtering, thefilterednoise, n(t) can be written as n(t) = x (t)cos® t-y (t)sm(£> t n c n c (3.1.2) where x (t) and y (t) are independent stationary zero mean Gaussian lowpass processes. We n n denote r(t) as the sum of Sj(t) and n(t): r(t) = s (t) + n(t). f (3.1.3) The received signal, denoted by r(i) at the output of the IF filter is the sum of the filtered noise and the filtered signal after fading , which should have been distorted by the filter. In the 5 The signal after fading will be referred to as signal, hereafter. Chapter 3 BER Performance of GMSK in Rayleigh Fading 24 'unfiltered' case, we assume that the IF filter does not distort the signal Sj(t) and hence the received signal is equal to r(t). r(r) = r(0 = x (t)cos((a t s + x (r)cosco ? n n s s s + s c + ^ (t)) c s -y (t)smty (t)]cos(i) t s + x (t)cos(d t s c = [x (t)costy (t) - [x (t)sm§ (t) s y (t)smo) t c n + ^ (t))-y (t)sm((i) t c s c y (t)cos$ (t)]sin(O t s s c - y {t) sinco r n c = X(r)cosco r- Y(t)s'm(O t c = c Re{Z(t)exp(ja> t)} c where Z(0 = X(t) + jY(t) The processes x (t), (3.1.5) X(t) = x (t) cos^CO-y,(0sin<|),(0 + x (t) , (3.1.6) Y(t) = x (t)sm<$> (t) + y (t)cos<b (t) + y (t) . (3.1.7) y (t), s , s s n s x (t) n s and y (t) n s s n are independent and of zero mean, hence the follow- ing relations can be obtained: x (t) = = y (t) s s Xi(t) (t) yi = y (r) = 0 B = yjijxjf) = Xi(t)Xj(t) = yjfjyjf) = 0 *,-(0y/0 = y,.(f)*;(0 = 0 x (t)y (t i j + ^ ) = Xiit + fyyjit) = {t) {t xJiTfyy~(7) 3 yi Xj + t) = y,(r+ £)*,•( f) = 0 = y.(t + Z )x (t) = 0 > *,•(* +$)*,•(') = {t yi where i, j = s,n i + Z) (t) yi = a?p.(cj j which will be used later in this section. In (3.1.8), the 'bar' denotes statistical average. Chapter 3 BER Performance of GMSK in Rayleigh Fading 25 Using (3.1.5), (3.1.6) and (3.1.7), Z(f) can be rewritten as Z(0 = x (t)cos(<b (t))-y (t)sm($ (t)) s s s (3.1.9) s + j[x (t) sin(<|>,(0) + y (f)cos(<MO)] nW + 7>„(0 = Z^OexpEy^COl+Z^O + x s 5 where Z (t) and Z (t) are independent complex zero mean Gaussian processes given by s n = Z(t) z n(0 = x.(t)+jy.(t), (3.1.10) x (t)+jy (t). (3.1.11) n n The output of the detector V'(t) (see Figure 3.1) is given by V'(t) = r(t)-r(t-T)\ n (3.1.12) positive phase shift = Re{Z(t)exp(j«> t)} • c Re{Z(t-T)exp(j(d (t-T))ex (jn/2)} c V where Z(r) is defined in (3.1.5) and r(t-T)\ = n Re{Z(t-T)exp(ja (t-T))} c - positive phase shift positive phase shift = /?e{Z(r-7/)}cosco (f-7/)| c I - positive phase shift Im{Z(t-T)}smm(t-T) (3.1.13) positive phase shift = Re{Z(t-T)cos(oJt-T = + n/2)}-Im{Z(t-T)sin(oJt-T + n/2)} Re{Z(t-T)expU® (t-T))txp(jn/2)} c Representing Z(t) by its real and imaginary parts as in (3.1.5) and multiplying out (3.1.12) Chapter 3 BER Performance of GMSK in Rayleigh Fading 26 yields, V\t) = Re{Z(t)expU® t)} • c Re{Z(t-T)exv)(j(£> (t-T))exip(jn/2)} c = (X(r)cosco ?-r(r)sinco O . (-X(t - T) sin® (t c c = - -X(t)Y(t-T)[cos(4nf t-2nf T) + l c - T) - Y(t - T)cos(& (t - T)) c c c cos(2nf T)] c (3.1.14) 1 + ^Y{t)X(t- T)[-cos{AKf t-2nf T) c - -X(t)X(t- T)[sm(4nf tl c + + c 2nf T) c c - sin(2nf T)] + sm(2Kf T)] -Y(t)Y(t-T)[sin(4llf t-2Kf T) l c cos(2nf T)] c c c In (3.1.14), the product contains signal components at around 2f , which are rejected by c the lowpass filter to yield V(r) (see Figure 3.1). Assuming that f c 2nf T c is chosen such that = (O T = 2%m where m is a positive integer, the output of the lowpass filter is c V(t) = = ±Re[-j(X(t) = \x(t)Y(t-T)+ -Y{t)X(t-T) l + jY(t))(X(t-T)- jY(t- T))} -Re[-jZ{t)Z*(t-T)] l (3.1.15) Chapter 3 BER Performance of GMSK in Rayleigh Fading 27 Expanding (3.1.15) and substituting (3.1.9) for Z(t) gives V(t) = ±Re{-j[Z (t)exp(Jt(t))+Z (t)]s n [Z*(t-T)exp(-j<$> (t-T)) + s = Z *(t-T)]} n -Re{-jZ (t)Z*{t-T)zx j($ {t)-$ {t-T)) + Z {t)} X s V \ = -Re{ [- jx {t)x {t -T)- s jy (t)y (t -T)- l s s s s N (3.1.16) T)] • x (t)y (t -T) + y (t)x (t- s s s s s [cos A(|> (r) + jsin A^(0] + Z (t)} 5 N = ±[x (t)x (t s -T) + y (t)y (t- s s + \[y (t)x (ts s T)-x (t)y (t- s s s T)] sinA$ (t) s T)]cosA^ (t) s +n\t) where Z (t) represents all terms associated with the noise processes Z (t) and Z *(t - T), and N n n n'(t) is the real part of Z (t). N The sampled output of the receiver V(IT) (see Figure 3.1), is obtained by sampling V(t) at time t=lT. The phase change over the 1-bit period is A(t> = J;/ 4>,(zr)-4>(zr-r) a oo = 2nf d oo IT X bj J g(v-jT)dv j = -oc (l-l)T = 2izf d IT-jT (3.1.17) £ bj J g(v)dv j = -oo a-\)T-jT where g(t) is the response of the premodulation Gaussian filter with impulse response h (t) t given by (2.1.13) to a transmitted rectangular pulse and bj is the jth transmitted symbol, i.e. 8(t) = (3.1.18) Chapter 3 BER Performance of GMSK in Rayleigh Fading 28 Equation (3.1.17) can be written as oo Ws,i= 1 j i-j j = -~ b (3- - ) Q 1 19 where 0/_ is equal to ; i-j = ^rf Q J ZlJ(t-l)T 8(u-JT)du. (3.1.20) When 0.2 < BfT < 0.3, 9,_ -is negligible for | / - j\ > 3 [22] (Table D . l in Appendix D shows the values of 0, *Ki • for B T= 0.25), therefore A$ t s l = b^Q.i +b i ^ + b f i o can be simplified to +b ^ + b , ^ Table D.3 in Appendix D shows the differential phase angles Aty s l (3.1.21) for all possible input sequence combinations. The transmitted signal s(t) after propagation through an ideal channel where there is no fading and no noise, is just the signal at the output of the modulator, defined in equation (2.1.4) in Chapter 2 and repeated here for convenience: s(t) = cos(co r + (j)(0)c (3.1.22) 5 Comparing (3.1.1) with (3.1.22), it is clear that the lowpass Gaussian processes x (t) and y (t) in s (3.1.1) is equal to 1 and 0, respectively for all values of t, when the channel is not fading. s Chapter 3 BER Performance of GMSK in Rayleigh Fading Hence, in an ideal channel, the receiver output V will be 29 6 V = ^sinA^ . (3.1.23) In Table D.3 (Appendix D), it is noted that Aty (or equivalently sinAc^ ) is positive when the st t transmitted bit b is T and negative when b is '0'. Consequently, the decision rule at the t l receiver is: decide the transmitted bit is a T if V > 0 and vice versa. In an ideal channel, correct decisions are always made and therefore the decision b\ made at time t = IT is always equal to the corresponding transmitted bit, denoted as b . t After propagation through the mobile channel, the transmitted bandpass signal spectrum is being shaped by the Doppler spectrum defined in (2.2.2) in Chapter 2 and the bandpass noise spectrum is being shaped by the Gaussian IF filter with 3-dB bandwidth B , and frequency t response H^{f) given by H (f) if where H(f) = H(f-f )+H(-f-f ) c (3.1.24) c is the equivalent baseband frequency response of the filter: tf(/) = Kexp{-(^J.f}. (3.1.25) In (3.1.25), K is an arbitrary constant. Assuming that the receiver has an omnidirectional antenna [9], the equivalent baseband spectra of the unmodulated carrier fading signal and noise, 7 Hereafter, we omit the time dependence of V on T, assuming it to be understood. Chapter 3 BER Performance of GMSK in Rayleigh Fading 30 denoted as W (f) and W (f), respectively, can be obtained by performing a frequency translas tion of / n (the carrier frequency) towards the baseband, on their bandpass spectra: 2 ^ W,(/) = i— D \J \f\<f \—Ju 0 \f\>f D 2 W (f) = —^=exp\ n -[ ) ln2 \ (3.1.27) Wln2 where//) is the maximum Doppler frequency. In (3.1.27), W (f) is obtained by multiplying the n one-sided noise spectral density N by H(f) defined in (3.1.25) and choosing K, where a a 2 n = NB 0 (3.1.28) n and B is the equivalent noise bandwidth given by n B„ = (3.1.29) The probability density function (p.d.f.) of the sampled output of the receiver V can be derived by taking the inverse Fourier Transform of the characteristic function F(v) of Vas follows [10]. Define a column matrix Z 7 The unmodulated carrier fading signal is equal to Re{Z (t)exp(ja> t)} , where Z^(f) is defined in (3.1.10). s c Chapter 3 BER Performance of GMSK in Rayleigh Fading z= 31 Z{IT) (3.1.30) Z((l-\)T)_ z 2 Its covariance matrix C is denned to be z C (3.1.31) = i(Z-Z)*(Z-Z) 2 where the bar denotes taking the statistical average. Using the relations defined in (3.1.8), it is found that Z j = Z = 0 and hence C can be written as 2 z C ^ = 2 7 * ^1 7 Zj z * z,*z Z *Zj 2 =a 2 2 z 2 (3.1.32) 1 p* IP 1. where p is the correlation coefficient between Z and Z j [9],[21], and is given by 2 a p = a p (r)exp(;A(l)^) 2 2 5 In (3.1.32), a 5 + ^p„(r) (3.1.33) represents the average power of the received signal, which is given by: 2 2 2 2 2 (3.1.34) In (3.1.34) and (3.1.33), a , G„ are the average powers of the signal, Sj-(t), and noise, respecs Chapter 3 BER Performance of GMSK in Rayleigh Fading 32 tively; p (t) and p„(r) are the normalized autocorrelation functions of the unmodulated carrier s fading signal and noise, respectively. We start by obtaining p (T) and p (T) , which we find by s n taking the inverse Fourier Transform of the spectra defined in (3.1.26) and (3.1.27), respectively. This yields [9] p (T) = J (2nf T) s 0 , D p„(T) = exp (3.1.35) (3.1.36) ln2 2 2 Defining the average signal-to-noise ratio (SNR) to be S = <5 /<5 , p is evaluated as s p = ^^y ( ^D^ 2 o e x PO' A ( l>, /) ) +( 5 ^ ) e x P{-(^„ : n r ) 2 / l n 2 } • C - - ) 3 1 3 7 The characteristic function F(v) of the sampled output V is defined as a Fourier Transform on the p.d.f. of V and can be shown to be equal to [39] F(v) = [det(I-2jvC *Q)] (3.1.38) _ 1 where I is a 2 x 2 identity matrix and Q is the 2 x 2 Hermitian matrix of the quadratic form of V. Using (3.1.38), F(v) is evaluated in Appendix A as F(v) = [l-2va (p-p*) + 4vV(l-|p| )] 2 2 1 . (3.1.39) Solving the quadratic equation in the denominator of F(v), the two poles of the character- Chapter 3 BER Performance of GMSK in Rayleigh Fading 33 istic function can be found to be Kpj±Ji^p ) r ^1.2 2 (3.1.40) , ,2 2a(l-|p| ) 2 where p and py are the real and imaginary parts of p , respectively. The p.d.f. of the receiver r output V is equal to P(V) = 1 (3.1.41) ^JF(v)e- dv JvV Using (3.1.41), it is shown in Appendix B that the conditional p.d.f.'s of V are given by 1 exp 2a (l-|p| ) 2 y<o 2 (3.1.42) P{V\b = ±\) = l ( 1 2 !"Pr exp 2cr (l-|p| ) 2 y>o 2 In (3.1.42), p and py are the real and imaginary parts of p (see (3.1.37)), respectively, with p r being a function of Aty 8 s1 which depends on the input bit, b , as described earlier. t According to the decision rule described earlier, error occurs if V < 0 when b\ = 1 or vice versa. The probability of errors for each of the cases are given as follows: In the FS fading case, p is also a function of A<() . , which depends on b,, as defined in Section 3.2. Chapter 3 BER Performance of GMSK in Rayleigh Fading P = P ( V < 0 | 6 = 1) = \P(V\b,= el Z 34 \)dv = \ 1 - (3.1.43) ( P e2 = P(V>0|A, = - 1 ) = J/>(V|i, = - l ) A = I (3.1.44) 1+—^ The unconditional probability of error is (3.1.45) Substituting (3.1.37) into (3.1.43) and (3.1.44), (3.1.45) is equal to P e X 7 ( 2 7 t / r ) s i n ( 0 + A(p ) o 5 t / l - [ X / ( 2 7 t / r ) c o s ( 0 + A(p, ) + Fexp{-(7i5 T) /ln2}] 2 A + D o = 2 f 1 1 + 2 V D o o >/ i : (3.1.46) X 7 ( 2 7 c / r ) s i n ( - 0 + A(p ,) D o 1 - [XJ (2nf T)cos(0 D O 5 6 + Acp, {) + yexp{-(7Cfi T)Vln2}] 0 I S 1 where X = -—-, Y = -—-, 0 represents the signal and A(p represents the ISI terms: o -f* 1 O o "H 1 A( P,,/ = / 5t ' 2-2 + ^ 10-l+Vl l+V2 2e & + 0 + e (3.1.47) Chapter 3 BER Performance of GMSK in Rayleigh Fading 35 3.1.2 Unfiltered Case with Decision Feedback HO Lowpass Filter * F Figure 3.2 Delay T 1 90° F Sample F F V(t) V(t) Hard Decision V(IT) or V F phase shift Delay T Block diagram of a 1-bit DF GMSK receiver Decision feedback receiver improves the system performance by partially removing the effect of ISI. At the sampling instant t = IT, when b is to be decided, the receiver has already t made a decision, bi_\on the transmitted bit A phase shift X equal to t>i_\Q x introduced to eliminate the ISI contribution due to b _ § . l in states 4 and 5 since b ^b _ l [ x x can be The effect of ISI is always destructive and therefore these two states are located closest to the thresh- l old, as shown in Figure 3.3(a). With decision feedback, we move these two states away from the threshold. The new differential angles for these two states are as follows: _ J rV/ + ~ \A$ -9J A( D F ^' A< s 1 S z 0 l for state 3 for state 4 (3.1.48) Chapter 3 BER Performance of GMSK in Rayleigh Fading 36 In doing so, some of the states are moved closer to the threshold (see Figure 3.3(b) and Table D.3 in Appendix D). However, this only causes minimal B E R performance degradation because states 3 and 4 are the critical states in determining the performance. Therefore, the unconditional probability of error is identical to (3.1.46) with A(j) 5t DF replaced by A$ s l as follows: = bfi. + A ^ (3.1.49) where %j A = b l + Q_ 2 2 + b Q_ +b ^Q l+l 1 l l + b _ Q -b - Q l 2 2 [ l l (3.1.50) Chapter 3 BER Performance of GMSK in Rayleigh Fading 37 3.1.3 Filtered Case In the previous section, we considered the 'unfiltered' case where the IF filter's effect is only to bandlimit the noise and not to distort the signal. In reality, the IF filter at the receiver distorts the signal's phase and amplitude. Its effect can be taken into account by convolving the received signal with the filter's impulse response h^(t) given by h (t) = 2Re{h(t)exv}(j(i) t)} if (3.1.51) c where h(t) is the equivalent lowpass impulse response of the IF filter, given by the inverse Fourier Transform of the filter's baseband frequency response H(f) defined in (3.1.25) 2n (B/2) t ~ 2 2 2 (3.1.52) hT2 where B has been defined in Section 3.1.1. i The 'filtered' received signal from the output of the IF filter is denoted as r(t) which we can express in a similar form as r(t) defined by equation (3.1.4) for the "unfiltered" case (Section 3.1.1): r{t) = X(r)cosco r-Y(t)sm(£> t c = c , (3.1.53) Re{Z (t)exp(j(D t)} if c where Z^(r) is the complex envelope of r(t). The 'filtered' signal r(t) is related to the 'unfiltered' signal r(t) by the convolution integral: Chapter 3 BER Performance of GMSK in Rayleigh Fading ~r(t) = 38 r(t)®h (t) if oo = J r(u)h {t - u)du • if (3.1.54) —oo oo = j Re{Z(u) exp(y'co w)} • Re{h(t - u)exp(j(d (t - u))}du c c —oo Assuming that the center frequency of the IF filter is the carrier frequency, it can be shown that the filtering effect is the same regardless of whether the signal is at baseband or at IF. That is, the bandpass system can be equivalently represented by the complex envelope Z(f) given in (3.1.9) and the lowpass function h(t) [26]. The complex envelope Z^(t) of r(t) is related to the complex envelope Z(t) of r(t) by the following convolution: Z (t) if = Z{t)®h{t) (3.1.55) oo = J Z(u)h(t - u)du —oo where <8> denotes convolution. It is only required to perform the convolution integral defined in (3.1.55) over the interval, denoted as T , for which the value of h(t) is not negligible. Assuming that the fading rate is t chosen such that it is much less than the bit rate, i.e. | / | « \/T , D Z (t) s h then the complex amplitude defined in (3.1.10) is essentially constant over the interval T . In this case, Z^(t) in the absence of noise, denoted by Z^Jt) is given by i Chapter 3 BER Performance of GMSK in Rayleigh Fading Z' (t) = if 39 \Z {u)ex [j{$ {u))]h{t-u)du s V s = [x (t) + jy (t)] j exp[j(^ (u))]h(t s s = x (t) j cos (u))h(t s - u)du s - u)du s Ti (3.1.56) - y (t) J s'm($ (u))h(t - u)du s s Ti + jx (t) j sin(§ (u))h(t s - u)du s Ti + jy (t) j cos(§ (u))h(t s - u)du s Consequently using (3.1.56), the received signal r(t) in (3.1.54) can be expressed as ~r{t) = Re{Z {t)&xxy(j(a t)} if c = Re{(Z' {t) + Z {t))ex (ja t)} if n V c n = Re{[{x (t)A {t)-y (t)A {t) s c j(x (t)A (t) s s + s = Re{(X(t) 1 + x (t))+ s n y (t)A (t)+y (t))]expU(O t)} s + c n c jY(t))ex Ua t)} V c where Z [t) is the zero mean complex Gaussian process for noise. n From above, it can be seen that Z (t) = if X(t) Y(t) +j = Z (t)A (t)+Z (t) s C0 n where X(t) = x (t)A (t)-y (t)A (t)+x (t), s c s s Y(t) = x (t)A (t)+y (t)A (t) s A(t) s s (3.1.59) n c = A(t) + jA(t), + y (t), n (3.1.60) (3.1.61) Chapter 3 BER Performance of GMSK in Rayleigh Fading 40 A (t) = cos<|> (0®MO c 5 (3.1.62) = j" cos (u))h(t - u)du s T, A (t) s = sin<^(r) ® h(t) (3.1.63) = j sin((j) (M))/i(r - u)du 5 Evaluating (3.1.57) yields: Re{(X(t)+jY(t)) x (j<» t)} r(0 = = e V c X(r)cosco r-T(f)sina> r c c = x M) A (t) cos a t-xM) AM) sin a t - [y (t)A (t)cos(£> t + y (t)A (t)sin(£> t] c s s c c . c s c c (3.1.64) + x (t) cos a> t - y (t) sin (£> t n c n c The term xM)AM)cos(OJ-xM)A(t)sm(o t (3.1.65) can be written as x (t)A(t)cos§ (t)cos(dj-x (t)A(t)sin$ (t)sin(xi t s s s s ^ \ 66) c = xM)A(t)cos((i) t + is(*)) where A M) and A (t) can be expressed in terms of A(t) and § (t); s AM) = A(0cos^(r), (3.1.67) AM) = A(r)sin<t>,(r); (3.1.68) Chapter 3 BER Performance of GMSK in Rayleigh Fading 41 and ty (t) = tan" s •sin<|> (0®VO5 = tan A (r)J A (t) = A (t)+A (t) 2 2 [cos^(t)®h (t)] 2 2 if (3.1.69) Lcos^r)® c + [sin^(t)®h (t)] 2 if . (3.1.70) Similarly, the term y„(f)A (f)cosco r+ y_(r)A (r)sina> f (3.1.71) c can be written as y (f)A(r)sin(^ (r)cos co f + y (f)A(r)cos0 (Osinco r J i c J iS (3.1.72) c = y (0A(0sin(co f + <j),(0) s c Hence using (3.1.66) and (3.1.72), (3.1.64) can be written as ~r{t) = x (t)A(t)cos(<d t + s c ty (t))-y (t)A(t)sm(G) t s s + x (t) cos a t - y (t) sin on t n c n c + <k (t)) s (3.1.73) c The output of the 1-bit differential detector is given by (3.1.15), where Z(f) is now replaced with Z (t) if defined in (3.1.58): Chapter 3 BER Performance of GMSK in Rayleigh Fading V(t) = 42 ±Re{-jZ (t)Z *(t-T)} if = if -Re{-j(Z (t)A (t) + Z (t))(Z*(t- l s co T)A *(t-T)+ n = -Re{- j(Z (t)Z*(t - T)A (t)A *(t- l s C0 Z*(t- co T))} T)) + Z (t)} C0 N 1 (3 1 74) = -{-AM)Ap-T)x {t)x (t-T) s - A (t)A (tc c + s T)x (t)y (ts c - A (t)A (tc where A {t), co s s s s s c s - T)x (t)y (t s s s + A (t)A (t - T)y (t)x (t -T) + A (t)A (t- T)y (t)y (t s T)y (t)y (t s T) - A (t)A (t s + A (t)A {t-T)y (t)x {t-T) c A {t)A (t-T)x (t)x (t-T) s s s s c s s s s - T) - T) -T) + ri{t)} A (t) and A (t) are defined in (3.1.61), (3.1.67) and (3.1.68), respectively. It is c s noted that A (t)A (t-T)-A (t)A {t-T) s c = c s ( 3 1 7 5 ) A(t)A(t-r)sin(A<i),(0) and A (t)A (t-T) c + c = A (t)A (t-T) s s (3.1.76) A(t)A(t-T)cos(A$ (t)) s where A M O = i (t)-Ut-T). (3.1.77) s Hence using (3.1.75) and (3.1.76), V(r)from (3.1.74) can be expressed as V(t) = l{A(t)A(t-T)[x (t)x (t-T)+y (t)y (t-T)]sm(A^ ) s s s s s . 2 + A(t)A(t-T)[y(t)x(t-T)-x(t)y(t-T)]cos(A^ ) s + n'(t)} (3.1.78) Chapter 3 BER Performance of GMSK in Rayleigh Fading 43 The phase change over the 1-bit period from t = ( / - 1)7" to t = IT at the Ith sampling interval, denoted as A§ s A^, z t is now given by = MIT) - MIT ~T) (3.1.79) J J sin(<)) (u))h(lT - T sin(§ (u))h(lT-u)du s s J) = tan J -u)du Ti tan J cos((j) (u))h(lT - T - u)du cos (u))h(lT-u)du s s where the integration is taken over the interval T for which h(.) is not negligible. i As in Section (3.1.1), the probability of error is given by equations (3.1.43), (3.1.44) and (3.1.45); where p and p are the real and imaginary parts of p , which is now given by r ; Gjfjjp = A A a p ( r ) e x p ( ; A ( ^ ) + c p „ ( 7 ' ) 2 1 2 J 2 / .2 a, = Co = A] = A 2 The parameters p (T) s and p (T) 2 2 .2 2 2°s °n> A 2 2 + A (IT-T) (3.1.80) (3.1.81) (3.1.82) 2 = (3.1.83) A (IT) 2 are given in (3.1.35) and (3.1.36). The average powers of the n 2 signal and noise are denoted as o s 2 and c , respectively. Equations (3.1.80), (3.1.81)and (3.1.82) n are derived as shown in Appendix C. It is noted that in the unfdtered case, the IF fdter can be 44 Chapter 3 BER Performance of GMSK in Rayleigh Fading considered to be an ideal bandpass filter with equivalent noise bandwidth B . In this situation, n A(t) in (3.1.70) will be equal to 1 for all values of t and the signal's phase is undistorted i.e. A<j) / = ;. This brings (3.1.80) back to (3.1.33) and (3.1.81) will be equal to (3.1.82), which 5; is identical to (3.1.34). Therefore, the unconditional probability of error can be obtained from (3.1.45) as r P e = 2 f \_ 2 U b,= \ 1 ) -(W , + exp{-(7c5 r) /ln2}) ; 2 V 2 bi= 2 M f 1 1 + 2 + V A (3.1.84) /(a a /a ) -(W 2 1 2 2 _ + exp{-(7t5 7') /ln2}) ; 2 fci = 1 2 n where (3.1.85) (3.1.86) b, = ±i = 5A A y (2K/ r)cosA(|), , w 1 2 0 D ;Z ,2 and S = G2/o . The subscripts in U and W refer to the corresponding input bit b for the s l n parameter A<|) /, which is a function of b . 5> l With decision feedback, the new differential angles are given as ;DF (3.1.87) where bi-\ is the decision made by the receiver on the transmitted bit b _ . The unconditional t x Chapter 3 BER Performance of GMSK in Rayleigh Fading 45 probability of error in the case of decision feedback is identical to (3.1.84) with A ^ / replaced by s> In the sections which deal with the performance of G M S K with D F in a FS channel only, the B E R analysis for 'filtered' G M S K with D F will be considered. This is because the 'filtered' analysis is more realistic as it takes into account the filtering effect. Moreover, both analytical and simulation results show that the differences in B E R between the 'filtered' and 'unfiltered' cases with D F is very small. Chapter 3 BER Performance of GMSK in Rayleigh Fading 46 3.2 B E R Performance of 1-bit detection in FS Rayleigh fading channel 3.2.1 Unfiltered case In a commonly used frequency selective fading channel [42], the received signal is the sum of two incoming signals, each with independent Rayleigh fading characteristics: the desired signal Sjr(t) given in (3.1.1) and a delayed signal J^(f) which also undergoes Rayleigh fading that is independent of that for s^(t). This is known as the two-ray model as described in the previous chapter. The delayed signal d^(t) is given by d (t) = x (t-x)cos(a> (t-x) f d + ^ (t-T))-y (t-x)cos(oi (t-x) c s d c + ^ (t-%)) s . (3.2.1) In (3.2.1), xjt) and y (t) are independent stationary zero mean Gaussian lowpass processes. The d time delay difference between the two signals is denoted as t. In the 'unfiltered' case, the received signal denoted by r(t) at the output of the IFfilteris equal to r(t), which is given by the sum of Sf(t), df(t) and the noise n(t) defined in (3.1.2). r(t) = s (t) + d (t) + n(t) f f = x (0cos(oy + <|>(0) ~ y (0sin(cB r + <|>(0) 5 5 c 4 5 + x (t - x)cos(fo (r - T) + § (t - x)) - y (t - %) cos(co (r - T) + § (t - x)) d c s d c s • (3-2.2) + x (t) cos co t -y (t) sin ® t n c n c In (3.2.2), if the fixed phase shift of co x is included in the random process x (t-%) and c d y (t -1), then r(t) can be expressed as d r(f) = X(0cosco r-y(f)sinco ? c c = /?«{Z(0exp(y(0 0} c (3.2.3) Chapter 3 BER Performance of GMSK in Rayleigh Fading 47 where Z(r) = X(t)+jY(t) X(t) = , (3.2.4) x (t)cos§ (t)-y (t)sm<b (t) s s s s + x (t - x)cos())(r - x) - y (t - x) sin d - x) + x (t) d J n Y{t) = ^(Osin^CO+y^Ocos^CO + x (t - x) sin <!>,(* - x) + y (t - x) cos0 (? - x) + y (t) d d n 5 We can re-write (3.2.4) using complex envelope notation as Z(r) = ^(Ocos^CO-y^Osin^CO + x (t - x)cos(j)(r - x) - y (t - x) s i n ^ f - x) + j[x (t) sin 0,(0 + ^(Ocos^CO + x (t - x) sin^Cr - x) + y (t - x) cos - x)] + x (t) + jy (t) = Z (t) exp [y<|>(0] + Z (t - T)exp[./<t> (f - *)] + Z (t) d d J s d d n n s d 5 5 n where Z (t), Z {t) and Z (t) are independent zero mean complex Gaussian processes. The s d n processes Z (t) and Z (r)are given by (3.1.10) and (3.1.11), respectively and Z^(0is given by s n d(t) = x (t) + jy (t). Z d d (3.2.8) The output of the lowpass filter V(t) in the 1-bit differential receiver is given in (3.1.15). Substituting Z(r) from (3.2.7) yields Chapter 3 BER Performance of GMSK in Rayleigh Fading V(t) = ^Re{-j[Z (t)cx (j^ (t)) s V + Z (t-x) xp(j^{t-x)) s d [Z*(t-T) xp(-j^ (t-T)) e l s d + Z (t)]n Z *(t-x-T)cxp(-j^(t-x-T))+Z *(t-T)]} d n - T)ex j(ty (t) -6) (t- T)) V s x - T) e x p y X K O + Z (t)Z *(ts e + s = -Re{- j[Z (t)Z*(t 48 s -ty (t-x-T)) s + Z (t - x)Z*(t - 7>xpy(<|>,(' -1) - W ~ T)) d + Z (t-x)Z *(td x - T)expM (t-x) d - 4>,(f - x - T))] + Z {t)} s N (3.2.9) where Z (t) represents all terms associated with the noise processes Z (t) and Z *(tN n n T). Substituting (3.1.10) and (3.2.8) for Z (t) and Z (t) into (3.2.9) gives s V(t) = -{(x (t)x (t-T) l s s d + y (t)y (t-T)) mA$ (t) s s S s + (x (t- x)x (t- x - T) + y (t-x)y {t-x d d + d - T))s\nA$ (t) d d (-x (t)y (t-T)+y (t)x (t-T))co A^ (t) s s s + (- x (t - x)y (t -x-T)+ d s s d d x d y (t)y (ts s x)x (t - x - 7j)cosA(^(r) d x-T)+ y (t)x (t- d + {x (t)x {t -x-T)+ s S y (t- d + (- x (t)y (t- s (3.2.10) - r))cosA(|)'(0 x - T)) sin Af(r) d + (-x (t - x)y (t -T) + y {t - x)x {t- T))cosAc>"(0 d s d s + (x (t- x)x (t -T) + y (t - x)y {t - T)) sinA<|>"(f) + n\t)} d s d s where A ^ O = <$> (t)-<b (t-T), s A<M0 = s $ (t-T)-$ (t-i-T), s s (3.2.11) (3.2.12) Chapter 3 BER Performance of GMSK in Rayleigh Fading A<t>'(0 = <b (t)-<$> (t-T-T), (3.2.13) A<|)"(0 = ( ^ - t ) - ^ - ! ' ) , (3.2.14) s and n'(t) is the real part of 49 s Z (t). N The sampled output of the receiver V is obtained by sampling V(t) at time t = IT. The 2 and the correlation coefficient p can be inferred from (3.1.34) and (3.1.33) signal power c s evaluated for the flat fading 'unfiltered' case (Section 3.1.1), but with the effect of second delayed signal dj-(t) also taken into account, as described earlier in this section. This yields: a 2 = a + o^ + a , 2 (3.2.15) 2 a p = a p,(r)exp(;A(t), ) + a ^ ( r ) e x p ( 7 A ^ ) + a p „ , 2 2 2 2 (3.2.16) 2 )/ / 2 where o , o and c are the average powers of the desired signal, delayed signal and noise, s d n respectively and p (r), p (r) and p (r) are the normalized autocorrelation functions of the unmodus d n lated carrier of the desired fading signal, the unmodulated carrier of the delayed fading signal 9 and noise, respectively. Assuming that the desired signal and the delayed signal experience the same Doppler spreading, then p (T) is equal to pd(T) = jQ(2nf T). s D We can observe that the probability of error defined in (3.1.45) in the FS case is related with the desired signal phase change A ^ given in (3.1.17) and the delayed signal phase change l The unmodulated carrier of the delayed fading signal is equal to Re{Z (t)exp(j(O t)}, where Z (t) is denned in d (3.2.8). c d Chapter 3 BER Performance of GMSK in Rayleigh Fading 50 A<j)^ , which is the phase difference over the 1-bit period from (l-l)T-T to IT-i. A§ ; d t is given by A^, = ^(IT-T)-MIT-T-T) . oo = 2nf d IT-X X bj j = -oo oo J (/_i)7-_ g(v-jT)dv = 2nf £ d bj J j = -°° t (3.2.17) IT-x-jT g(v)dv (l-i)T-x-jT Equation (3.2.17) can be written as oo Wd,i = 2 bfii-j (-- ) 3 2 18 where 0* is equal to 6/_ - = g(u-jT)du 7 ^ and is the/ fc (3.2.19) 2IJ(l-\)T-x transmitted symbol. As with Aty , A$ s t d l can be simplified to $d, l = l + 6-2 + A b 2 + 10-1 + M o + b _ j 0! + b _ 0 . l l 2 (3.2.20) 2 2 2 We define the average signal-to-noise ratio (SNR) to be S = o /o s n and the average 2 2 desired (carrier) signal-to-delayed signal ratio (CDR) to be S and D into (3.2.16), p is given as D = c /cj . By direct substitution of s d Chapter 3 BER Performance of GMSK in Rayleigh Fading P /) + [ [DS + S + ) o(^fDT)^V(m , = J D s 51 )j (2Kf T)exp(jAh, S DS S+ + D 0 D /) [DS fs + ^Vi-^Bjf/iXnl)} + (3-2.21) The probability of error is given in (3.1.45). By substituting (3.2.21) for p (hence p and r p •), (3.1.21) for Aty and (3.2.17) for A$ sl we can derive the unconditional probability of d error in FS Rayleigh fading channel for the unfiltered case as P„ = (Pel + Pel) (3.2.22) r f 1 X sin(0 + A(p ) + X s i n ( 0 + Aq ,) } 1 Jl- o s z 2 o d [ X c o s ( 0 + A(pi;,) + X c o s ( 0 + Acp^,) + Fexp{-(7tfl,7/) V l n 2 } ] 1 o 2 o f 1 Xj sin(-0 + A(p ,) + X s i n ( - 0 + Acp^,) 1 + 2 . / l - t Z j C o s ^ O o + Acp^ ,) + X c o s ( - 0 + A(p ) + 7exp{-(7t5 .r) /ln2}] ; V o + s> 2 O 2 2 y = = D S + D s D' ° J + i 2 n f D T ) ' A ( ? s > 1 i s § i v e n m + 2 2 z+ rf> z ( (3-1-47) and P</,/ = b[ 0 _ + & ,0_i + A( O fc §i M + ^/_ 02 • 2 (3.2.23) 2 Chapter 3 BER Performance of GMSK in Rayleigh Fading 52 3.2.2 Filtered case We can apply a similar analysis described in the 'filtered' case for flat fading channel (section 3.1.3) by including the IF filtering effect on the signal and the delayed signal components in a FS channel. The 'filtered' received signal from the output of the IF filter is denoted as r(r) which is given by r(t) = Re{Z (t)expUG> t)} if where Z^(t) (3.2.24) c is the complex envelope of r(t). As described in the 'filtered' case for flat fading, the complex envelope Z y(r) is related to the complex envelope Z(f) of the 'unfiltered' received ( signal r(t) by equation (3.1.55), where Z(t) is given by (3.2.4). The received signal r{t) will have a similar form as that derived for the flat 'filtered' case in Section 3.1.3, with the terms associated with the delayed signal also taken into account. Consequently, the received signal given in (3.2.24) can be expressed as ~r(t) = = Re{Z (t)exp(j® t)} Re{[x (t)A (t)-y (t)A (t) if c s c s + x (t - x)A (t -x)-y (td + c d s x)A (t — x)+ x (t)] s (3-2-25) n j[x (t)A (t)+y (t)A (t) s s s c + x (t - x)A (t — x) + y (t - x)A (t — x)+ y (t)] exp(;0) r)} d s = Re{(X(t) + where Z (t) n d c n jY(t))exp(j(O t)} c is the zero mean complex Gaussian lowpass process for noise. From above, it can be seen that c Chapter 3 BER Performance of GMSK in Rayleigh Fading 53 Z (t) = X(t) + jY(t) = Z,(t)A (t) + Z (t - T)A (t- X ) + Z„(0 (3.2.26) if co d co where X(t) = x (t)A (t)-y (t)A (t) + x (t- x)A {t- x) -y (ts c s d s c d s s d d , (3.2.28) s c s (3.2.27) x)A (t - x) + xM) •7(0 = x (t)A (t) + y (t)A (t) + x (t - x)A (t - x) + y (t - x)A {t- x) + y {t) s , c n and A (t) is defined in (3.1.61), A (t) in (3.1.62) and A (t) in (3.1.63). co c s Alternatively, r(t) can be expressed as ~r{t) = x (t)A(t) cos ((a t + ^(0)-y (r)A(r)sin(co f + ^(r)) s c 5 c + x (t - x)A(t- x)cos(co r + § (t - x)) - y (t - x)A(t - x)cos(oo f + <j>(f - x)) d c s d c + x (t)cos(x) t - y (t) sinco r n c n 5 (3.2.29) c where the filtered signal's phase (j>(f) and the receivefilteredsignal amplitude A(t) are given by s (3.1.69) and (3.1.70), respectively. The terms A (t) and A (t) can be expressed in terms of A{t) c s and ^ ( f ) , as given in (3.1.67) and (3.1.68). The output of the 1-bit differential detector is given by (3.1.15), where Z(t) is now replaced with Z y(r) defined in (3.2.26): ( Chapter 3 BER Performance of GMSK in Rayleigh Fading V(t) = ±Re{-jZ 't)Z *(t-T)} if = if -Re{-j(Z (t)A (t) s co d T)A *(t-x-T)+Z*{t-T))} co (3.2.30) ^Re{-j(Z (t)Z*(t-T)A (t)A *(t-T) s co co + Z (t - x)Z*(t - T)A (t - x)A *(t- + Z (t)Z *(t C0 d co s d co - x - T)A (t)A *(tC0 + Z (t-x)Z/(t-x- T) x - T) T)A (t-x)A *(t-x-T))+ d where Z (t) co -T) + Z/(t-x- C0 1 + Z„(f)) • + Z (t-x)A (t-x) l (Z*(t - T)A *(t = 54 co Z (t)} co N represents all terms associated with noise processes Z (t) N n Substituting (3.1.10) for Z (t), (3.2.8) for Z (t), s d and Z *(tn T). (3.1.11) for Z „ ( f ) , (3.1.61) for A (t) co and multiplying out (3.2.30) gives V(t) = |[^(0^(f-r) + y (0y (r-7')]A(0A(r-7')sin(A^(0) J 5 + [x (t - x)x (t -x-T) d + y (t- d + x)y {t -x-T)]A(t-x)A(t-x- d d T) sin(A^(0) [y (t)x (t-T)-x (t)y (t-T)]A(t)A(t-T)cos(AUt)) s s s + [x (t- x)x (t- T) + y (td s d d s T) + y (t)y (t-x- d s + [y (t)x (t-x-T)-x (t)y (t-xs d s x {t - x)y (t- T)]A(t- s + [x (t)x (t-xs x)y (t - T)]A(t- d + [y (t- T)x (t -T)- s s d d x)A{t- T)sin(A^"(0) x)A(t - 7)cos(A^"(0) T)]A(t)A(t-x- r)sin(A^'(0) T)]A(t)A(t — x — T)cos(A<j>;(0) + [y (t-x)x (t-x-T)-x (t-x)y (t-x-T)]A(t-x)A(t-x-T)cos(A$ (t)) d d d d where A(j) (r) is given by (3.1.77) and s (3.2.31) d + n'(t) Chapter 3 BER Performance of GMSK in Rayleigh Fading 55 (3.2.32) A<t>'(0 = 6)s(t)-$ (t-i--T), (3.2.33) A<|>"(f) = faft-T)-$ (t-T). (3.2.34) s s At the Ith sampling interval, the desired differential phase change over the 1-bit period, Aty I is given by (3.1.79) and the delayed differential phase angle, A ^ / is given as s ; (3.2.35) J sin(()) (u))h(lT - x - u)du s = tan J cos(§ {u))h(lT - x - u)du s | sm(§ (u))h(lT - x - T - u)du s - tan -l Jj J cos(§ (u))h(lT - x - T - u)du s T, where T- denotes the time interval over which h(,) is not negligible. The probability of error is given in (3.1.45) where p is given by a , o p = A A a p,(r)expOAt>, )+A A a^(7;)exp(jA^,/) + a 2 2 and 1 2 J Tl t2 2 p„(r) (3.2.36) Chapter 3 BER Performance of GMSK in Rayleigh Fading a 2 .2 2 .2 2 .2 2 .2 2 x2 = A\ = A (IT) A = A (lT-T (3.2.37) 2 (3.2.38) +a. , d A\ n A {IT-T) 2 2 = A (lT + x) 2 2 (3.2.39) + x) 2 xX A\ s s 2 2 = A a +A a 2 2 p (T) 2 56 and p (T) are given in (3.1.35) and (3.1.36). The average powers of the signal, Sjr(t), and n 2 noise are denoted as o 2 and G , respectively. In Appendix C, the derivations of (3.2.37), (3.2.38) s n and (3.2.39) are shown. Therefore, the total probability of error can be obtained from (3.1.45) as P. = x r f 1 12 L V Ui, i + bl= — ) -(W 1 ) f c ( = 1 +W 2 ; U 2ibi=1 , ( = 1 exp{-(7u5.T)/ln2}) + J ( 3 2 4 Q ) f + 1 1+ , 2 ) " (W 2 lt where f/j fc _ ± 1 and _ fc ± 1 b i = _ + exp{-(7t5.r) /ln2}) ; _ +W x 2 2> b i = 2 x are given in (3.1.85) and (3.1.86), respectively and 2,b, = ±l = (3.2.41) N \ 2 o( fD ) ^d,h U D W2,b, = ±i = A A J 2K T sinA (3.2.42) D A A J (2Kf T)cosAt^ i. In (3.2.41) and (3.2.42), the parameter A§ N d> l 2 0 D di / depends on the corresponding input bit b and D l N Chapter 3 BER Performance of GMSK in Rayleigh Fading 2 is the average delayed signal to noise ratio, defined as 57 2 c /a . d n With decision feedback, the new differential angles are given as ~ DF ~ ~ A^/ = A^z-d,^.^!, = A i ^ - d ^ e p where d x (3.2.43) (3.2.44) is a variable that can be adjusted to provide the optimum feedback angle so that the B E R is minimized. This is because in FS fading, the phase shift £>/_ i0] is not always optimum owed to the presence of the delayed signal. The total probability of error in this case is identical to ~ (3.2.40) with A(j) ~DF iS> ~ ~DF / replaced by A ^ ; and A(j)^ / by Aty^ /. Chapter 3 BER Performance of GMSK in Rayleigh Fading 58 3.3 B E R Performance of 2-bit detection in flat Rayleigh fading Figure 3.4 and Figure 3.5 depict the transmitter and receiver for the 2-bit detection case. The transmitter is similar to that for the 1-bit detection case except that it includes a differential encoder before the premodulation filter. This is because if the differential encoder is absent, the decisions made by the detector will be a differentially encoded version of the input data sequence and not the original input. This will become clearer in the subsequent section, where the mathematical analysis is shown. From Figure 3.4, we have bj = Cj • C y _ j since bj is either T or r 1 Cj-l Figure 3.4 ' — • Gaussian lowpass filter FM modulator s(t) T Block diagram of 2-bit GMSK transmitter (for differential detection) The output of the 2-bit detector V'(t) is obtained from multiplying the output of the IF filter, r(t) by a 2-bit period delayed version of itself. The product is then passed through a lowpass filter to produce V(t), which is assumed to be sampled at time t = IT to give V(IT), where / denotes the Ith sampling interval: le {...,-1,0, 1,2, ...}. A hard decision is then made on V(IT) to decide whether a '1' or a '0' was sent. In an ideal channel (no fading and no noise), the receiver output is given by (3.3.9). The operation of the 2-bit differential detector will be described mathematically in section 3.3.1. Chapter 3 BER Performance of GMSK in Rayleigh Fading 59 Detector Delay 2T V(0 V(IT) Hard Decision Sj(t): signal after fading with average power ^ s "0" or "!' n(t): AWGN which produces narrowband noise with 2 average power °„ at the IF filter output Figure 3.5 Block diagram of a 2-bit GMSK differential receiver 3.3.1 Unfiltered case In the 'unfiltered' case, the received signal is equal to r{t) which is given by (3.1.4) and the complex envelope Z(t) of r(t) is given by (3.1.5) or equivalently (3.1.9). Since we assume the IF Gaussian filter has no effect on the received signal, except for bandlimiting the noise, therefore the signal's phase is undistorted at the output of the filter. The output of the detector V'(t) (see Figure 3.5) is given by Chapter 3 BER Performance of GMSK in Rayleigh Fading V(t) = [r(t)-r(t-2T)] = Re{Z(t)exp(j(ti t)} 60 • Re{Z{t-2T)ex {j(Si {t-2T))} c V • c ( 3 3 > 1 ) Multiplying out (3.3.1) yields V"(0 = Re{Z(t)exp(j(£> t)} • c Re{Z(t-2T)exp(ja> (t-2T))} c = (X(t) cos a y - 7(0 sin a y ) • (X(t - 2T) cosco (f -2T)-Y(t-2T) sinco (f - 2T)) c = ±X(t)X(t-2T)[cos(4nf t - 2nf T) c c + c cos(2nf T)] c (3 3 2^ -±Y(t)X(t-2T)[sm(4%f t-2nf T) c + c -±X(t)Y(t-2T)[sm(4nf t-2Kf T)-sm(2nf T)] + -Y{t)Y(t-2T)[- cos(4nf t-2Kf T) c c c c + X c sm(2nf T)] c cos(2nf T)] c The lowpass filter rejects the signal components shown in (3.3.2) at around Assuming that f c is chosen such that 2nf T c 2f . c = (£> T = 2%m where ra is a positive integer, the C output of the lowpass filter is V(0 = ^[X(t)X(t-2T) + Y(t)Y(t-2T)] = -Re[(X(t) + jY(t))(X(t-2T)-jY(t-2T))] l = ±Re{Z(t)Z*(t-2T)} Substituting (3.1.9) for Z(t) into (3.3.3) gives • (3.3.3) Chapter 3 BER Performance of GMSK in Rayleigh Fading V(t) = 61 -Re{[Z (t)cxp(JW)+Z (t)]- l s [Z*(t - 2T)exv(-j<$> (t- IT)) + Z*(t n - 2T)]} s (3.3.4) - 27>xpy(<M0 " <M'" = ±Re{Z (t)Z *(t s = s -[x (t)x (t-2T) s s + \ [x (t)y (ts where Z (t) N + l » + Z ^)> y (t)y (t-2T)]cosA<b (t) s s 2T) - y (t)x(t- s 2 r s 2T)] sin A^(f) + n'(t) s represents all terms associated with the noise processes Z (t) n and n'{t) is the real part of and Z *(tn 2T), Z (t). N The sampled output of the receiver V(IT) (denoted as V hereafter) is obtained by sampling V(t) at t = IT. The phase change over the 2-bit period from t = (l-2)T to t = IT, denoted as A(j) where 5l A(|> l ' = <b (lT)-ty (lT-2T) s s = fd 2K (3.3.5) IT-jT oo X j c J S(v)dv 0-2)T-jT ; = -•» and g(r) is the response of the premodulation Gaussian filter to a transmitted rectangular pulse, given in (3.1.18). Equation (3.3.5) can be written as oo A t,i = X JVI-J C (3.3.6) Chapter 3 BER Performance of GMSK in Rayleigh Fading 62 where V(/,_ is equal to ; IT Vl.j = 2^1 8(u-jT)du. (3.3.7) (/-2)7 When 0.2 < 5-7 < 0.3 , \|/,__. is negligible for |/ - ; | > 3 [22] (Table D.4 in Appendix D shows the values of 4>,,/ A for 5 T = 0.25), therefore Aty ( = c / + 2¥-2 + can be simplified to sl / i¥-i+c ¥ c + / 0 + c _ \ l / + c _ \)/ + z 1 1 / 2 2 . c _3\|/3 / (3.3.8) As in the 1-bit case, the transmitted signal s(t) after propagation through an ideal channel is defined by (3.1.22). Comparing (3.1.22) with (3.3.4), the lowpass Gaussian processes x (t) and y (t) s s in (3.3.4) is equal to 1 and 0, respectively for all values of t, when the channel is not fading. Hence, in an ideal channel, the receiver output V at time IT will be V = \cosA$ . (3.3.9) Stl It is noted from Table D.4 that when 0 < |A<t>^. | < 7 i / 2 , the product c _ c z l c are of opposite signs), and when 7t/2 < |A(t) [\<n, c _ c z 5 signs). Since we have b = c _ c l l l l l x l l l is -1 (c _ l l and is 1 (c _ j and c are of equal t t , the decision rule at the receiver is: decide the transmitted bit b, is a ' 1' when V < 0 and vice versa. Therefore, Chapter 3 BER Performance of GMSK in Rayleigh Fading $s,l = ± ( V o + V i ) A <bs,l A = ± ( V 0 + A( rV/ = ±2\|/ - V l ) + A ( P 5 , / = 0 + A ( 0 63 given b = 1 given b - -1 p5,/ t %,l + A (3.3.10) l where \\f is equal to \\f (from Table D.4) and 0 x (3.3.11) Following a similar procedure as in the 1-bit detection case, we define a column matrix Z : z= 2," Z(IT) z Z((/-2)T) (3.3.12) 2 Its covariance matrix C is given by (3.1.32), namely 1 p* IP (3.3.13) 1. where p is the correlation coefficient between Z and Z j , and is given by: 2 ap =G 2 2 ( 2 7 > x p a A ^ ,) + a p„(27;). 2 Pj The average power of the received signal denoted as a 2 (3.3.14) given by (3.1.34). The average 2 powers of the signal, Sjr(t) and noise are denoted as G , G , respectively; and the normalised s n autocorrelation functions of the unmodulated carrier fading signal and noise are denoted as p (t) s and p „ ( r ) , respectively. The parameters p (2T) and p „ ( 2 T ) are evaluated as [10]: s p (2T) = J (4nf T) s Q D , (3.3.15) Chapter 3 BER Performance of GMSK in Rayleigh Fading p„(27/) = exp {-(InBiT) 64 Vln2} (3.3.16) 2 , 2 Defining the average signal-to-noise ratio (SNR) to be S = G /G , s P = [^]/ (47t/ r)exp(jA(t), ) + (^ 0 D ;/ n p is evaluated to be )exp{-(27Ufi„r) /(ln2)} 2 T . (3.3.17) The characteristic function F(v) of the sampled output V is given in (3.1.38). For the 2bit detection case, it is evaluated to be: F(v) = [ l - 2 > G ( p + p*) + 4 v o - ( l - | p | ) f . 2 2 4 2 1 (3.3.18) Solving the quadratic equation in the denominator of F(v), the two poles of the characteristic function can be found to be: _ 7(P ±7T-p]) r P l < 2 ~ 2 (3.3.19) , ,2 2a (l-|p| ) 2 2 where p and p - are the real and imaginary parts of p , respectively. The p.d.f. of the receiver r y output V, P(V), is given by (3.1.41). Following the analysis shown in Appendix B for 1-bit detection, P(V) for 2-bit detection can be evaluated to be exp P{V\b = ±\) = 4a V<0 2a (l-|p| ) 2 . 2 (3.3.20) l exp 4 a 2 ^ 2 •V V>0 , 2a (l-|p| ) 2 2 Chapter 3 BER Performance of GMSK in Rayleigh Fading 65 In (3.3.20), p and pj are the real and imaginary parts of p (see (3.3.17)), respectively with p r being a function of A$ 1 0 s 1 which depends on the input bit, b , as shown in (3.3.10). t According to the decision rule described earlier, error occurs if V > 0 when b = 1 and l vice versa, the probability of errors for each of the cases are given as follows: Pr f P el = \ 1- = P(V<0|fc, = - l ) = \P{V\b^-l)dv A (3.3.21) P V oo Pel = P(V>0\b =l) = jP(V\b^l)dv [ =\ (3.3.22) V The unconditional probability of error is given by P„ = (Pel + Pel) = \ e\+\P (3.3.23) \{( el\^s, l = V P 2 since for the case when b = T , A§ l s t + 0 A( P*, /) \( el\ ^, + P A / = ~ Vo + <P,, /) 2 A can have either of the two values, as shown in (3.3.10). Substituting (3.3.14) into (3.3.21) and (3.3.22), (3.3.23) yields In the FS fading case, p is also a function of A<(> . , which depends on b, , defined in Section 3.4. Chapter 3 BER Performance of GMSK in Rayleigh Fading 66 X7 (47t/ r)cos(A(p ) + yexp{-(2nB .r) V l n 2 } D 0 5 I Jl-[XJ (4nf T)sm(A<? )] 2 0 + D s XJ (4nf T)cos(2\\f 1 + 1 2 0 D Q 5 l Jl - [XJ (4nf T)sin(2i|/ 0 D + A(p )] 0 (3.3.24) 2 5 X7 (47t/ r)cos(-2v|/ + Aq> ) + 7exp{-(2jcB,-r) V l n 2 } 1 1 + 2 Jl - [XJ (4nf T)sm(-2y + Acp )] V 0 + + Aq> ) + rexp{-(27ifl T) V l n 2 } D 5 0 2 0 where X = 5+1 ,Y = D 0 5 5+1 3.3.2 Unfiltered Case with Decision Feedback < 1 F Delay 2T F Sample Lowpass Filter V(0 Hard F F V(IT) V(0 Delay T A, "t-l N Delay T "1-2 Figure 3.6 Block diagram of a 2-bit DF GMSK receiver Decision feedback can also be applied to 2-bit detection. Observe from Figure 3.7(a) that states 17 and 24 are subjected to the greatest amount of ISI, therefore the objective is to move these two states and also other states away as far as possible from the threshold. The phase shift X, can be implemented as follows [22]: Chapter 3 BER Performance of GMSK in Rayleigh Fading 0 67 ifc _ = C _ / / 2 1 It is noted that the biased threshold is only employed in a non-fading channel with A W G N , as described in the beginning of this chapter. In the case of fading, we will only consider the zerolevel threshold because there is little performance improvement in a fading channel with a biased threshold. Decision Threshold Decision Threshold 24 22,23 22,24 9,10 ! 21 20 ol r 21,23 ! \ T 1 / 18,19 Figure 3.7 1 1 1 1 I ^ 17 (a) 1 / k W 18,20 j 11,12, 13,14 10,11, 14,15 15,16 12,16 1 biased threshold (b) Phase-state diagram of 2-bit differential detection: (a) before feedback (b) after feedback Therefore, the new probability of error is identical to (3.3.24) with A$ s t replaced by DF A$ s l as follows: DF A where C z = /Vo + / - i ¥ i + C / c c A (3.3.26) Chapter 3 BER Performance of GMSK in Rayleigh Fading A ( P 5 , / l 2^-2 c = + Cl- V2 + + 68 l-3^3-^ • C 2 (3.3.27) 3.3.3 Filtered case The input to the receiver, r(t) is given by (3.1.64) in section 3.1.3. It is also given by (3.1.57) when expressed in complex envelope notation where Z^(t) is given by (3.1.58). The output of the 2-bit differential detector is given by (3.3.3), where Z(r) is replaced with Z^(r): V(t) = ±Re{Z (t)Z *(t-2T)} if if = ±Re{(Z (t)A (t) + Z (t))(Z*(t = ±Re{Z (t)Z*(t- 2T)A (t)A *(t- s C0 s = co -Re{A {t)A (tc s - A (t)A (t- 2T)x (t)y (t + A (t)A (t + A (t)A (t c c c where A (t), co s s c A (t) c -2T) + Z „ * ( f ) ) } co 2T) + Z (t)} co 2T)x (t)x (t l c - 2T)A *(t n N - 2T) + A (t)A (t s s - 2T)x (t)x (t s s s - 2T) + A (t)A (t - 2T)x (t)y {t - 2T) - 2T)y (t)x (t- 2T) - A (t)A (t - 2T)y (t)x (t -2T) - 2T)y (t)y (t- 2 7 ) + A (t)A (t- s s s s s s s and A (t) c s s c s s s s 2T)y (t)y (t- s s s - 2T) { 3 3 2 S ) 2T) + n'(t)} are defined in (3.1.61), (3.1.62) and (3.1.63), respectively. It is s noted that A (t)A (t c c = - 2T) + A (t)A (ts s 2T) (3.3.29) A(t)A(t-2T)cos(A$ (t)) s and A {t)A {t-2T)-A {t)A {t-2T) s = c c s A(t)A(t-27/)sin(A(j> (r)) i (3.3.30) 69 Chapter 3 BER Performance of GMSK in Rayleigh Fading where A^(0 = Ut)-Ut-2T). (3.3.31) Hence using (3.3.29) and (3.3.30), V(t) from (3.3.28) can be expressed as V(t) = i { A ( 0 A ( r - 2 r ) [ x ( 0 ^ ( f - 2 r ) + y , ( 0 y s ( r - 2 r ) ] c o s ( A i ( 0 ) J + A(t)A(t-2T)[x s J (t)y(t-2T)-y (t)x(t-2T)]sin(AUt)) s ( 3 3 3 2 ) + n'(t)} The phase change over the 2-bit period from t = (I - 2)T to t = IT at the Ith sampling interval, denoted as A(j) , is now given by 5 A<|>.z = l UIT)-UIT-2T) (3.3.33) J sm(§ (u))h(lT = tan J sin(§ (u))h(lT - 2T - u)du -u)du s T s h - tan J J cos (0 {u))h(lT -2T cos(§ (u))h(lT-u)du s s - u)du L.T, where the integration is taken over the interval T for which h(.) is not negligible. i The analogue of (3.1.80) for the 2-bit detection case is given as a,a p = A A a p ( 2 r ) e x p O A ^ ) + a p„(27') , 2 2 1 2 i 2 > / (3.3.34) Chapter 3 BER Performance of GMSK in Rayleigh Fading 70 where 2 .2 2 a, = A a x s (3.3.35) n 2 °~2 = .2 2 2 20s+°n> (3.3.36) A\ = A (IT-2T) (3.3.37) A = A (IT) 2 The parameters p (2T) 2 + <5 , s A 2 2 2 and p (2T) are given in (3.3.15) and (3.3.16). Equations (3.3.34), n (3.3.35) and (3.3.36) can be derived in a similar manner as that for 1-bit detection, shown in Appendix C. As in Section 3.3.1, the probability of error is given by equations (3.3.21), (3.3.22) and (3.3.23), where p and py are the real and imaginary parts of p , which is given by (3.3.34). r Therefore, the unconditional probability of error can be obtained as f W 1+— b + exp{-(27i5 r) /ln2} — 2 = 1 (<*l°2 n) - b,= /CT w 1 1 2 V fc/ = + A n l U 2 2 J (3.3.38) \ _ + exp{-(27Cfi r)Vln2} 1 li (G^/a ) '-U 2 2 _ 2 b i = where b, = ±i = SA^jQiAnf^sinAtysj, u W b, = ±i = 5 A A 7 ( 4 j c / r ) c o s A ^ / , 1 2 0 o i (3.3.39) (3.3.40) Chapter 3 BER Performance of GMSK in Rayleigh Fading 2 2 and S = a /a . s 71 n The subscripts in U and W refer to the corresponding input bit b for the l parameter A§ /, which is a function of c and b = c • c _ . With decision feedback, the new differential angle is given as Si { t t l x A ^ / = A^j-l F (3.3.41) where A. is given in (3.3.25). The total probability of error in the case of decision feedback is identical to (3.3.38) with Acj)^ / replaced by A^j. 3.4 BER Performance of 2-bit Detection in FS Rayleigh fading channel 3.4.1 Unfiltered case In a commonly used frequency selective fading channel model, the received signal is the sum of two incoming signals, each with independent Rayleigh fading characteristics. It is given by (3.2.2) or equivalently by (3.2.3) if expressed in complex envelope form. The complex envelope Z(r) of r(r) is given by (3.2.7). The output of the 2-bit differential detector is given by (3.3.3). Substituting (3.2.7) for Z(t) into (3.3.3) yields Chapter 3 BER Performance of GMSK in Rayleigh Fading 72 V(t) = -Re{Z{t)Z*{t-2T)} l = ±Re{(Z (t)expj($ (t))+Z (t-x)expM (t-x)) s s d •(Z M-2T)txp(-j$ (t-2T)) s -Re{Z {t)Z*{t- n + Z *(t-x-2T)exp{-j§M-x-2T)) s = + Z (t)) s 2r)exp;((t),(0 - l s + d 2 T )) (3.4.1) 2 x - 2T)expM (t) + Z (t)Z *(ts d - $ (t - T - 27)) s s + Z (t - x)Z*(t- 2T)cxpM (t d -x)- s x)Z *(t - x - 2T)expj($ (t + Z (td Z*{t-2T))} d ^(f - 2D) -x)- s 27)) + Z (t)} <$> (t -xs N Substituting (3.1.10) and (3.2.8) for Z (t) and Z (t) into (3.4.1) gives s V(t) = ±{(x (t)x (t-2T) s + y (t)y (t-2T))cosA<$> (t) s s + (x (t)y (t - 2T) - y (t)x (t- + (x (t)x (t -x-2T)+ s s s s s + (x (t)y (t-x s d s s 2D) sinA^CO s y (t)y (ts x- s - 2D -y (t)x (t-x- s s 2T))cosA^' (t) d 2D)sinA$ (t) s d (3.4.2) + (x (t- x)x (t - 2 D + y (t - x)y (t - 2T)) cos A^" (t) d s d s d + (x (t - x)y (t - 2 D - y {t - x)x (td s d 2D) sinA$' (t) s d + (x (t - x)x (t - X-2T) + y (t - x)y {t -xd d d d 2D)cos A<^(r) + (x (t- x)y {t- x — 2T) — y {t - x)x {t -xd d d d 2D) sin A ^ O + n\t)} where A<|>(0 = $ (t)-$ {t-2T), 5 s A0 (O = <|>,(f d A$\t) (3.4.3) s x) — — <|>(f 5 X - 2D , = <|> (0-<Mf-T-2D, 5 (3-4.4) (3.4.5) Chapter 3 BER Performance of GMSK in Rayleigh Fading 73 A<|>"(0 = t(t-x)-Mt-2T), (3.4.6) and ri(t) is the real part of Z (t). N As in the 1-bit FS case (section 3.2.1), the probability of error depends on both the desired signal phase change A<j) given in (3.3.8) and the delayed signal phase change A(|) 5 A0 rf l = d i / t , where § {lT-x)-$ {lT-x-2T) s IT-x oo X = *fd 2 s j c J oo S(v-jT)dv = 2nf d U-2)T-t y = -°° (3.4.7) IT-x-jT ]T ^ J j = -°° g(v)dv (l-2)T-x-jT Equation (3.4.7) can be written as $d,i = X A where I S th (--) 3 4 8 g{u-jT)du T (3.4.9) transmitted differentially encoded symbol. As with A<|> i , A(|) s rf l can be simplified to <^d,l = l lV-2 + l \V-\ A $l-j equal to ,-=^r~ and cy is the j c c c + + + $0 + l-\V\ c C + l-lVl c + In (3.4.10), unlike the case of A§ , , \ j / is not equal to s 0 l-lVl- c (3.4.10) . The values of • for Chapter 3 BER Performance of GMSK in Rayleigh Fading 74 \l - j\ < 3 with B T = 0.25 and x = 0.2 T is given in Table D.5 in Appendix D. t 2 We define the average signal-to-noise ratio (SNR) to be 5 = cs/on 2 2 2 desired (carrier) signal-to-delayed signal ratio (CDR) to be D = a /(5 . s d and the average The analogue of (3.2.16) for the 2-bit case is given by a p = a, p (2r)expO-A^ ) + a p (2T)expaA^ ) + a p„. 2 2 2 5 2 rf / / (3.4.11) By direct substitution of 5 and D into (3.4.11), p is given as P = [DS + S S + D) D j Q ( 4 n f D T ) e X 1 ? U A \1 DS + S + D ^> l ) + [DS + S + j J o ( 7 /D?>xp(7A(^,) 4 7 I 2 ( e x p { - ( 2 7 t 5 „ T ' ) V(ln2)} The probability of error is given in (3.3.23) where P el and P e2 3 A 1 2 > are given by (3.3.21) and (3.3.22), respectively. By substituting (3.4.12) for p (hence p and p •), (3.3.8) for A<j> and r (3.4.10) for A(j) . In this instance, the unconditional probability of error is derived as rf ; 5>/ Chapter 3 BER Performance of GMSK in Rayleigh Fading P X,cos(2\)/ + Acp^j) +X cos(\j/ + vj/! + Acp^ ,) + Y 0 1+ =r L V 2 1 f V 0 j]/ 2 0 0 Jl-[X 5/ 2 0 0 5> V 1 2 rfi# 0 rf> X c o s A ( p + Z cos(\i/ -\}/ +Aq ) 1 5/ 2 0 1 cU +Y (3.4.13) ^ -[XjSinCAcp^;) + X s i n ( \ j / - \ j / + A<p ,)]' 2 1 12 d>/ sin(-2\|/ + Aq> ,) + X s i n ( - \ j / - ijjr, + A(p ,)]' l f + 1 X,cos(-2\|/ + A(p ) +X cos(-\j/ -\j/ + A(p ) + F 1_ ' 2 1 2 0 r 1 1+ 2 + 2 Jl -[A sin(2\(/ + A(p ) + A: sin(\j/ + vi/ +Aq> )] . 2 + 75 0 1 d> Z,cosA(p j + XjCost-ij/o + iij/, + A(p^ /) + T 5 Jl-[X sin(A(p ,) + X s i n ( - y + vj/, + A<p ,)]' l S) 2 0 d> where (3.4.14) (3.4.15) D DS+S+D v A( and 5 = c]/a ,D 2 n P</,/ =Z c + + 2V-2 c (2*W J Z l¥-l+C/_ V2 + (3.4.16) In 2 2 + c /-3 r 3' li / (3.4.17) = a)/a]. 3.4.2 Filtered case The input of the receiver, r(t) is given by (3.2.29) in section 3.2.2. It is also given by (3.2.25) when expressed in complex envelope notation where Z^(t) is given by (3.2.26). The Chapter 3 BER Performance of GMSK in Rayleigh Fading 76 output of the 2-bit differential detector is given by (3.3.3), where Z(r) is replaced with V(t) = ±Re{Z (t)Z *(t-2T)} if if = -Re{(Z (t)A (t) +Z (t-x)A {t- l s (Z*(t- C0 d 2T)A *(t-2T) s n C0 2T))} - x - 2T)A (t)A *(t co ( x)Z*(t - 2T)A (t - x)A *(t + Z (t- x)Z *(t - x - 2T)A (t - x)A *(t - x - 2T) + Z (t)} d AM) d co co and A (t) s 1 8 ) -2T) C0 N are defined in (3.1.61), (3.1.62) and (3.1.63), respectively, and represents all terms associated with the noise processes ZM) Substituting (3.1.10) for ZM), using (3.3.29) and (3.3.30), gives A -x-2T) co co 3 c + Z (td N + Z *(t- co -Re{ZM)Z *{t-2T)A {t)A *(t-2T) d Z (t) 2T)A *(t-x-2T) l + ZM)Z *(t C x) + ZM)) • co + Z/(t-x- co = where A M), ZM): (3.2.8) for Z (t), d and Z *(t - T). n (3.1.11) for Z (t) n into (3.4.18) and Chapter 3 BER Performance of GMSK in Rayleigh Fading V(t) = -A(t)A(t-2T)[x (t)x (t-2T) s + + s 77 y (t)y (t-2T)]cos(AUt)) s s A(t)A(t-2T)[xM)y (t-2T)-y (t)x (t-2T)]sin{AUt)) s + A(t)A(t- x - 2T)[x (t)x (ts s s x — 2T) + y (t)y (t - x - 27/)] cos(A(j>' (r)) d s d + A(t)A(t - x - 2T)[x (t)y (t — x — 2T) — y (t)x (t - x - 2T)] sin(A0'(O) s d s d + A(t - x)A(t- 2T)[x {t - x)x {t - 2T) + y (t - x)y (t - 27/)]cos(A<K'(0) d + s d s ( 3 A 1 9 ) A(t-x)A(t-2T)[x (t-x)y (t-2T)-y (t-x)x (t-2T))sm(A^(t)) d s d s + A(t - x)A(t -x- 2T)[x (t- x)x (t- x-2T) + y (t- x)y (t- x - 27/)]cos(A<jv(r)) d d d + A(t - x)A{t - x - 2T)[x (t- x)y (t -X-2T)d d d y (t- x)x (t — x — 2T)) sin(A^(0) d d where A<() (£) is given in (3.3.31) and 5 AWO = Mt-x)-Mt-x-2T) (3.4.20) (3.4.21) Af(0 = Ut-D-Ut-2T) (3.4.22) At the Ith sampling interval, the desired differential phase over the 2-bit period, A ^ / is given by (3.3.33) and the delayed differential phase angle, Atj)^ / is given as 78 Chapter 3 BER Performance of GMSK in Rayleigh Fading A(k/ = <$>S(IT-T)-UIT-T-2T) J sin (§ (u))h(lT - x - u)du s (3.4.23) Is = tan j cos(§ (u))h(lT— x — u)du s J sin(§ (u))h(lT s - tan-l - x- 2T - u)du Ii | cos (§ (u))h(lT -x-2T s - u)du where T is defined earlier. t The probability of error is given in (3.3.23) where p is given by a,a p = A A a p ( 2 r ) e x p ( y A ^ ) + A A a ^ ( 2 r ) e x p ( ; A ^ ) + a p „ ( 2 r ) 2 2 1 2 2 5 ( / T l x 2 ) / (3.4.24) and 2 .2 2 °1 = \ ^ s A + 2 .2 2 ° 2 = 2°s i2 d A A 2 A\ A xl A 2 2 +A = A .2 2 2 ^d n > .2 a (S 2 2 + °n » (3.4.25) (3.4.26) A (IT-2T) 2 = A (IT) 2 (3.4.27) = A (lT-2T 2 2 + = A (lT 2 + x) + x) p (2T) and p (27) are given in (3.3.15) and (3.3.16). The average powers of the signal, Sj(t) s w Chapter 3 BER Performance of GMSK in Rayleigh Fading 2 and noise are denoted as a 5 79 2 and c , respectively. Equations (3.4.24), (3.4.25) and (3.4.26) are n derived in a similar manner, as in 1-bit detection (Section 3.2.2). Similar to the analysis done for the 1-bit detection and using (3.3.21), (3.3.22), (3.3.23), the unconditional probability of error in the 2-bit case is obtained as Wi - i + W 1+— bl P. = r + 1 2 ! + exp{-(27ufi„r) /ln2} 2 2b Wi, b, = -i + 2, b, = - l + exp{-(27tZ?„r) /ln2} W 2 2~2 1 (c a /a ) x 2 2 — n (3.4.28) (C/i, fc/ = _ i -i- where ^ i , * , = ±i = SA A J X 2 {4nf T)smA$ , Q 2b, = ±\ = D (3.4.29) sJ (3.4.30) N l 2 o( fD ) ' $d,h U D A A J 4n T S mA W i,b, = ±i = 5 A A 7 ( 4 7 t / r ) c o s A ^ , 1 W 2,b, = ±l = 2 2 2 2 D 0 / N l 2 o( fD ) D A A J 4n T <bd,h COsA (3.4.31) (3.4.32) 2 and S = os/Gn, DN = od/an (the average delayed signal to noise power ratio). In (3.4.29) to (3.4.32), the parameters A$Sj / and A ( j ) ^ ; depend on the corresponding input bit . Chapter 3 BER Performance of GMSK in Rayleigh Fading 80 With decision feedback, the new differential angles are given as A<j)^ = Aty i - d X, (3.4.33) A$7,i = A4d,i-d X, (3.4.34) Si 2 2 where d is a variable that can be adjusted to provide the optimum feedback angle so that the 2 BER is minimized, as in 1-bit detection. The probability of error in the case of decision feedback ~ DF ~ ~ ~DF is identical to (3.4.28) with A(j) / replaced by A ^ , / and A<t>^ / by A ^ /. S) Chapter 4 Results and Discussions Chapter 4 Results and Discussions In this chapter, the B E R formulae derived in Chapter 3 are illustrated with numerical examples. In order to validate the formulae, comparisons with simulation results are also presented. 4.1 Simulation Model The transmitter consists of a random binary source and a G M S K modulator. The source output data bits are generated using a pseudo-random number generator. It should be pointed out that in 2-bit detection, the data sequence requires differential encoding. The simulation written in C, is implemented as an equivalent baseband model, which has been described in Chapter 3. The baseband transmitted signal is given in (3.1.9), with Z (t) s and Z (t) n set to 1 and 0, respectively. The received signal, distorted by fading and noise, is then passed through the equivalent baseband Gaussian IF filter and then to the receiver implemented as a 1-bit differential detector for 1-bit detection or as a 2-bit differential detector for 2-bit detection. The zero mean random Gaussian noise process is generated using the method described in [44] and the generation of the random Gaussian processes (the in-phase and quadrature components) that characterize the Rayleigh fading waveform is described below. The Rayleigh fading waveform that characterizes multipath propagation is generated using the method described by Jakes [32]. The Jakes model is a discrete, time averaged approximation to the desired Doppler spectrum defined in (2.2.2) in the form of N os low-frequency oscillators. As described in the previous chapters, the in-phase and quadrature components, denoted as x (t) s and y (t) s respectively, are given by: 81 82 x {t) = 2 ^ cosP cos(co O + j2cos((a t) s n n= y (t) s (4.1.1) D n 1 = 2 2 sinP„cos(ay) (4.1.2) n= 1 where B„ = (N +l) n is to provide the appropriate gains for the sinusoids, and os (0„ = G) cos(2nn/N ). d In (4.1.3), N s = 2(2N 0S (4.1.3) s + 1) and a> is the maximum Doppler frequency in radians/sec. D The above technique can be extended to provide multiple independent fading processes [32]. Two independent Rayleigh fading waveforms can be generated and used to implement the FS two-ray model as described in Chapter 2. The following is a list of parameters being used in obtaining the numerical or simulation results: (1) sampling frequency (f ) s The 'sampling frequency' is the number of samples per second. It is set to 8 • 10 samples per 5 second in the simulations. (2) bit rate (R ) b The 'bit rate' is the transmission rate in number of bits per second. In the simulations, the bit rate is chosen to be 16000 bits per second. Hence, this gives 50 samples per bit using the value of f s 83 stated above. (3) signal power (P ) s The 'signal power' is the desired signal power at the point before A W G N is added. When the channel is non fading with only AWGN, the signal power is given by the transmitted power of the signal, s(t), defined in (2.1.4) and is set to 0.5; whereas in the case of Rayleigh fading, it is given by the average power of the signal, Sy(r), defined in (3.1.1) and is equal to 0.5(x + y ). s s 2 (4) postdetection noise power (o" ) n In this thesis, error rates are usually expressed as a function of E /N , b P/R b Q where E is equal to b in a non fading A W G N channel and it is an average value in the case of fading because the signal power (after fading) varies; and N is the one-sided PSD of A W G N . The received signal is 0 passed through the IF filter and the noise is therefore being bandlimited. The postdetection noise 2 . . , power <5 is given by n (4.1.4) n where B is the equivalent noise bandwidth of the IF filter and is equal to 2 A / l n ^ ^ ' ^ ' k * § ^ 1 n e n e 3-dB bandwidth of the IF filter. Throughout the simulations, B is chosen to be 0.9R because it i b is optimum for a varieties of channel parameters (in both fading and non fading channels) and B E R performances [5],[18],[20], hence B n is equal to 15328.8 Hz. 84 (5) normalised 3-dB bandwidth of the pre-modulation filter (B T) t The parameter B T is chosen to be 0.25. The reasons for choosing this value is three-fold: (1) this t value is assumed in a number of papers and thus allows a direct comparison between our results and previous works. (2) G M S K is made more spectrally compact by having a smaller B T (in the t G S M standard, B T = 0.3), at the cost of significant ISI. (3) the effectiveness of the decision t feedback scheme can be more clearly demonstrated, given the amount of inherent ISI. (6) number of low-frequency oscillators in the Jakes model in simulating Rayleigh fading (N ): os The number N os is set to 8. By choosing p n = — ' powers of x (t) s and y (t), denoted as x (7) doppler frequency s s — in (4.1.1) and (4.1.2), the average os ) and y , are given by N s os and N os + 1, respectively. (f ) D The Doppler frequency has been defined in (2.2.3). The Doppler frequency-time product fT D (also known as the normalized Doppler spread) is used as a parameter in the simulations, where T = l/R b is the bit duration. The values that are chosen for f T D are given in later sections. (8) time delay difference in a FS channel (x) The parameter x is the time delay difference between the main (undelayed) and the secondary (delayed) signal in a FS channel, as described in Chapter 3. Two values, 0.2Tand 0.47; are used in the simulations. 85 (9) decision feedback coefficients (d and d ) x 2 The decision feedback coefficient, d , for 1-bit DF is chosen to be 1 [22] as it minimizes the BER x for most of the cases of fading considered in this thesis. The corresponding parameter d , for 22 bit D F is chosen to be 2 in all the curves presented, except that in Figure 4.12. The simulation results presented in this thesis have 99% confidence intervals of approximately ± 1 % in a non fading AWGN channel, 95% confidence intervals of approximately ± 5 % in flat Rayleigh fading and 95% confidence intervals of approximately ± 1 0 % in FS Rayleigh fading. We will denote conventional (no decision feedback) G M S K as 'conventional' and G M S K with decision feedback as DF. In the figure legends, 't' stands for theoretical results, 's' for simulation results, 'c' for conventional detection and 'DF' for detection with decision feedback. 4.2 Error Performance in a Non Fading Environment In this section, the simulated BER and B K E R results of 1 and 2-bit 'conventional' in a non fading channel with A W G N are shown. A numerical procedure for calculating the B E R is given in [5]. For our purposes, it is more convenient to obtain the BER, which is used later to calculate the B K E R , from simulations. The simulation values were verified to agree with the values in [5]. Figure 4.1 shows the BER of 1 and 2-bit 'conventional' G M S K in a non fading environment with A W G N . In the 2-bit detection case, two curves are shown: the first curve results if the detection threshold is optimized for every E /N b a and the second curve is for a fixed detection threshold of zero [5]. As mentioned in Chapter 3, the B E R performance of 2-bit detection can be much improved when the threshold is optimized. The advantage is apparent when E /N b D is greater 86 than about 10 dB, as shown in Figure 4.1. A fixed zero threshold 2-bit 'conventional' shows a better performance than 1-bit 'conventional' when E /N b E /N b 0 is less than 17 dB. However, for greater than 17 dB, 1-bit 'conventional' has a better BER performance. 0 Eb/N (dB) 0 Figure 4.1 BER of 1 and 2-bit conventional GMSK in a non fading environment with AWGN. For 2-bit detection, 'thres = 0' denotes the case where a decision threshold of zero is used. The block error rate (BKER) is the probability that there is at least one bit error in a packet of N bits. Assuming independent bit errors, the B K E R is given by l-(\-P y e where P e (4.2.1) is the BER. The probability of more than M bit errors in a block of length N bits, denoted as P(M, N), is important in the performance evaluation when using block coding in an error control scheme. The CDF is equal to the probability that a block has less than or equal to M errors, given by 87 M Pr{number of errors < M} = ^ P inumber n=Q M -x n = of errors = n} r 0 (4.2.2) {y.w and Pr{number of errors <M} = 1 - P(M, N). (4.2.3) Figure 4.2 and Figure 4.3 show the B K E R for 1-bit and 2-bit 'conventional' in a non fading environment with AWGN, respectively. Both of the figures show that the B K E R for /V = 255 bits obtained from simulation is just slightly lower than that obtained from (4.2.1), using the 16 18 20 Eb/N Figure 4.2 0 22 (dB) BKER of 1-bit conventional GMSK in a non fading environment with AWGN corresponding simulated B E R . This is because for differential detection, the bit errors are correlated (rather than independent) as each decision depends on the current phase and the phase at 1-bit period or 2-bit periods earlier (depending on whether it is 1-bit or 2-bit differential detection). It is also observed that as N increases, the difference between the B K E R from (4.2.1) 88 and the simulated B K E R decreases. This observation agrees with that in [45] which investigates the block (packet) error rates of differentially encoded coherent BPSK, in which the problem of correlation between bit errors also arises. Eb/N (dB) 0 Figure 4.3 BKER of 2-bit conventional GMSK in a non fading environment with AWGN. Note that the threshold is optimized for every E^/N' . Figure 4.4 shows the C D F of the number of bit errors in a block at E /N b 0 = 16 dB in a non fading environment for 1-bit and 2-bit 'conventional', with N = 255, 511 and 1023 bits. In both cases, the CDFs obtained from simulations agree very well with the corresponding theoretical CDFs obtained from (4.2.2). This indicates that (4.2.2) (and (4.2.1)) is valid for approximating the block error statistics in a non fading environment, despite the fact that differential detection introduces some correlation between the bit errors. In summary, the performance of 2-bit detection can be much improved when the threshold is optimized. However, if a fixed zero threshold is used in 2-bit detection, it only shows a better performance than 1-bit detection for E /N b 0 < 17 dB. In both cases of detection, small correlation 89 (a) Figure 4.4 (b) CDF of the number of bit errors in a block at E^/N - 16 dB in a non fading environment for (a) g 1-bit conventional GMSK (b) 2-bit conventional GMSK with optimum threshold. exists between the bit errors, introduced by differential detection. 4.3 Comparison of BER Performances between 'filtered' and 'unfiltered' cases in Rayleigh Fading In a fading environment, the received signal power varies and therefore the parameter E /N b 0 refers to the ratio of average energy (after fading) per bit to the noise PSD. In this section, the average B E R performances are obtained from the numerical formulae derived in Chapter 3. For 1-bit detection, the 'unfiltered' and 'filtered' BERs in flat fading are obtained by (3.1.46) and (3.1.84), respectively. In FS fading, the corresponding BERs are given by (3.2.22) and (3.2.40). For 2-bit detection, the 'unfiltered' and 'filtered' BERs in flat fading are obtained by (3.3.24) and (3.3.38), respectively. In FS fading, they are given by (3.4.13) and (3.4.28). 90 Figure 4.5(a) compares the BER of 1-bit and 2-bit 'conventional' for both 'filtered' and 'unfiltered' cases in flat fading, with the Doppler frequency-time product f T as a variable D parameter for E /N b Q at 20 dB and at 60 dB. For both cases, we can see that at E /N b the B E R is almost insensitive to f T. However, at E /N D f T. D b This is because at lower values of E /N , b values of E /N b 0 0 Q = 20 dB, a = 60 dB, the B E R increases with the BER is noise-limited. Conversely, at high (>40 dB) the noise has a negligible effect on the BERs, which are now dominated by the effect of the Doppler spread. In subsequent sections, results that show irreducible error floors being present at E /N >40 b o dB with f T D > 0.001, are displayed. Error floors occur when the B E R cannot be further lowered by increasing E /N , b due to the effect of Doppler 0 spread. When the 'filtered' and the 'unfiltered' cases are compared: at E /N b Q = 60 dB, the 'filtered' B E R is on average about 25% higher than the 'unfiltered' BER of 1-bit detection for all values of f T D shown; and this difference is about 35% for 2-bit detection. This suggests that the BER for 2-bit detection is more sensitive to the IF filtering effect than that for the 1-bit case and can be explained by the fact that the B E R expressions depends on the parameter where n = 1, 2 , for 1-bit and 2-bit detection, respectively. Since J (.) 0 of its argument and a smaller / (.) 0 J (2nf nT) 0 D is a decreasing function means an increased B E R (see equations (3.1.84) and (3.3.38)), the B E R degradation through IFfilteringis more apparent in the 2-bit case. Figure 4.5(b) shows the B E R of 1 and 2-bit 'conventional' for both 'unfiltered' and 91 'filtered' cases as a function of f T, with average carrier to delayed signal power ratio D (CDR) = 20 dB and time delay difference x = 0.27 in a FS channel for E /N b 60 dB. In both cases, we see that at E /N b the flat fading case. At E /N b small values of f T D (f T D at 20 dB and at a = 20 dB, the BER is almost insensitive to f T, 0 D as in = 60 dB, the BER increases much more gradually, especially for 0 < 0.001) when compared with the flat fading case. This is because in a FS channel, the BER is controlled by the combined effect of f T, D CDR and x. 10'' - D 1 - b i t filtered - A 2-bit, filtered • 1 -bit, unfiltered | A - A Z - b i t , unfiltered 10' cr E./N b o E./N b o = 20dB LU LU CO CO = 60dB 10 ' - r 0.000 0.001 1-bit, filtered & 2 - b i t , filtered • 1-bit, unfiltered A 2 - b i t , unfiltered 0.002 0.003 (a) Figure 4.5 10"' 0.000 0.001 0.002 0.003 (b) Average BER as a function of f T for 'filtered' and 'unfiltered' cases of 1 and 2-bit D conventional GMSK in (a) flat fading (b) FS fading with CDR = 20 dB, x = 0.2T . When E /N b a = 60 dB and f T D < 0.001, the effects of Doppler and noise are negligible and the parameters CDR and x, which are fixed for the results shown in Figure 4.5(b), overwhelms the 92 Doppler effect on the BER performance. This explains why the B E R is almost insensitive at low values of f T. A comparison of the 'filtered' and the 'unfiltered' B E R shows that the B E R D 10"' 10" 1 - H 1 -bit, filtered - A 2-bit, filtered JQ - EH-bit, unfiltered IA - A2-bit, unfiltered -€H-biL filtered - A 2-bit, filtered In - EI1 4ML unfiltered |A - A 2 - b i L unfiltered 10' - E./N = 20dB b o cr UJ CO LU CB IO" 2 H E./N b o = 20dB E./N = 60dB b o 10" J E./N b 10"' 0.000 = 60dB o _L 0 . 0 0 1 _L 10" 0 . 0 0 0 ; 0 . 0 0 2 0 . 0 0 3 0 . 0 0 1 0 . 0 0 2 fnT 10" 1 E./N b o = 20dB -Ol-bitfWtered -A2-bit,Wtered • 1 -bit, unfiltered IA - A 2-bit unfiltered ciur 10 o 2 E./N = 60dB b o IO"- 0 . 0 0 1 0 . 0 0 0 Figure 4.6 0 . 0 0 2 0 . 0 0 3 Average BER as a function of / T for 'filtered' and 'unfiltered' cases of 1 and 2-bit D conventional GMSK in FS fading with (a) CDR = 20 dB, i = OAT (b) CDR = 10 dB, T = 0.2T (c) CDR = 10 dB, T = 0.4T . 0 . 0 0 3 93 degradation due to IF filtering is negligible when f T is small and it increases with f T D E /N b 0 for D = 60 dB, whereas it remains almost constant independent of f T D when E /N b 0 = 20 dB. 2-bit detection is also more sensitive to IFfiltering.Similar reasons as stated earlier are accounted for the observations. Similar observations can also be deduced from Figure 4.6. In the case of flat fading (Figure 4.5(a)), it is observed that at both E /N b = 20 dB and 60 dB, the percentage differ- 0 ence between the 'filtered' and the 'unfiltered' BER is almost constant for the values of CDR D CDR (a) Figure 4.7 f T. (b) Average BER as a function of CDR for 'filtered' and 'unfiltered' cases of 1 and 2-bit conventional GMSK in FS fading for f T = 0.002 with (a) D T = 0.2T (b) x = 0.4T . From Figure 4.7, it can be seen that at both E /N b a = 20 dB and 60 dB, the percentage difference between the B E R in the 'unfiltered' and 'filtered' cases increases with CDR. This is because for a small CDR, the relatively high delayed signal power tends to mask the effect of the IF filtering; but as C D R increases, the IF filtering effect becomes more apparent. At 94 E /N b = 6 0 d B , with CDR = 60dB for both x = 0.2T and x = OAT, the B E R s for the 0 'unfiltered' and 'filtered' cases approach the corresponding BERs of the flat fading case. This is expected because at such a high CDR, the delayed signal power is much weaker than that of the carrier signal power and can be considered to be negligible regardless of x, resembling the case of flat fading. In the case when E /N b = 20 dB, the BERs stay almost constant as the C D R o increases. This section investigated the effect of IF filtering on the BER performances of 1-bit and 2bit conventional G M S K as a function of f T D the BER is almost insensitive to f T and CDR. For both 'filtered' and 'unfiltered' cases, at low E /N D b but at high E /N , Q b by the effect of the Doppler spread, thus it increases with f T.A D 0 the BER is dominated comparison of the 'filtered' and the 'unfiltered' BER in FS fading shows that the BER degradation due to IFfilteringis negligible when f T D is small and it increases with f T D almost the same for all values of f T D observed that for a given E /N , b 0 for high E /N , b when E /N b 0 0 whereas the BER degradation is is small. In the case of flat fading, it is the percentage difference between the 'filtered' and the 'unfiltered' B E R is almost constant, independent of the values of f T D considered. In all the fading cases that were considered, 2-bit detection is more sensitive to IF filtering than 1-bit detection. 4.4 BER Performance in Slow Rayleigh Fading With slow fading, the signal is assumed to be constant for the duration of over at least one bit period, i.e. the product f T D should be much less than 1. In this case, f D is chosen to be 32 95 Hz, giving f T= D 0.002 where T is given by the reciprocal of the bit rate R . b Figure 4.8 shows the average B E R in a flat Rayleigh fading channel. Irreducible error floors start to appear at E /N b 0 ~ 45 dB . It can be seen from the figure that for 1 and 2-bit 'conventional' and 1-bit DF, the simulation results are in good agreement with the theoretical 30 Eb/N (dB) 0 Figure 4.8 Average BER as a function of E /N b g in flat fading. results as given by (3.1.84), (3.3.38) and (3.1.84)*, respectively. The '*' indicates that for the decision feedback schemes, the BER formulae are based on those of the corresponding 'conventional' detection, but with the pertinent values of differential angles substituted, as described in Chapter 3. In the case of 2-bit DF, the simulations B E R values are substantially higher than the theoretical results from (3.3.38)*. This is because (3.3.38)* assumes that all decisions being fed back are correct in order to simplify the analysis. However in the simulation, true decisions, some of which incorrect, are being fed back. In the 2-bit D F case, the phase shifting rule looks at two consecutive decisions which are then differentially encoded in order to determine the feedback 96 angle. Therefore, if one decision is incorrect, after differential encoding error propagation occurs to give double-bit errors. This explains why the simulated BER is higher than the theoretical BER from (3.3.38)*. To verify this explanation, a simulation was run in which all of the decisions being fed back are correct. In this case, it was found that the average BER for 2-bit DF is in close agreement with the analytical result from (3.3.38)*. The improvement in BER for 1-bit DF over 1-bit 'conventional' are about 55% and 60% at E /N b 0 = 30 dB and 60 dB, respectively. For 2-bit DF, the corresponding percentages are about 50% and 40%, respectively. It is observed that 2-bit 'conventional' performs slightly better than 1bit 'conventional' when E /N b 0 <25 dB. For example, 2-bit 'conventional' provides an approxi- mate 1 dB gain over 1-bit 'conventional' at a B E R = 1 0 . When E /N _1 b a ~ 25 d B , there is a switch over between the two performances: 1-bit 'conventional' now performs better than 2-bit 'conventional' and has a lower irreducible error floor. In the cases of FS fading, there also exist similar trends, as shown in subsequent figures. This is expected because the BER performances of 1-bit and 2-bit 'conventional' (with zero-level decision threshold) in a non fading A W G N channel shows a similar trend (see Figure 4.1). In the case of decision feedback, there is almost no difference between 1-bit D F and 2-bit D F when E /N b a <30 dB. However, beyond E /N b a = 30 dB, 1-bit D F has a better performance and has a lower error floor. We can compare the BER performances of 1 and 2-bit conventional G M S K with that of coherent detection obtained in [27] 11 in flat Rayleigh fading (known as the dynamic B E R 97 performance). When f T D at E /N b = 0.0025 1 2 (f D = 40 Hz), the BER obtained in [27] is about l x l O - 3 = 60 dB (the error floor), whereas they are 1.68 x 10" (from (3.1.84)) and 4 0 5.72 x 10" (from (3.3.38)) for 1 and 2-bit conventional G M S K , respectively. When f T 4 D = 2.5 x 1 0 " ( / = 4 H z ) , the BER from [27] is 1 x 10~ and the corresponding BERs for 1 and 24 5 D bit 'conventional' are 4.25 x 10 and 8.15 x 10 , respectively. Hence, in both cases, differen- 6 6 tial detection achieves a lower BER over coherent detection in fiat Rayleigh fading. The improvement comes from the fact that differential detection cancels the phase distortion between adjacent symbols during phase comparison. In FS fading, the theoretical BERs for 1 and 2-bit 'conventional'; 1 and 2-bit D F G M S K are given by (3.2.40), (3.4.28), (3.2.40)* and (3.4.28)*, respectively. Figure 4.9 shows the average B E R in FS Rayleigh fading channel with x = 0.2T and CDR = 20 dB. Irreducible error floors start to appear at around E /N b a = 40 dB. As in the flat fading case, the simulation results for 1 and 2-bit conventional G M S K and 1-bit D F G M S K are in good agreement with the theoretical results. In the case of 2-bit D F G M S K , the average B E R from simulation is about 50% higher than the theoretical results from (3.4.28)* for E /N b 0 > 30 dB. The degree of disagreement is less than that in the flat fading case. The improvements in BER of 1-bit D F G M S K over 1-bit conventional G M S K are about 50% and 15% at E /N b Q = 30 dB and 60 dB, respectively. In 2-bit D F G M S K , the corresponding percentages are about 45% and 20%, respectively. 1 1 The experimental results in [27] were reproduced in [12] by software simulation. 1 2 The bit duration Tis chosen to be 6.25 10 s , which is also assumed in this thesis. 5 98 Eb/N Figure 4.9 Average BER as a function of E /N' 0 (dB) in FS fading with CDR = 20 dB, x = 0.2T. b Figure 4.10 shows the average B E R in FS Rayleigh fading with x = OAT and CDR = 20 dB . Irreducible error floors occur at around E /N b a = 40 dB. There is good agreement between the simulation results and the theoretical results for 1 and 2-bit 'conventional' and 1-bit DF. In the case of 2-bit DF, the average B E R from simulation is about 30% higher than the corresponding theoretical result for E /N b > 30 dB. The improvement in B E R with 1-bit D F over 0 1-bit 'conventional' is about 40% at E /N b 0 = 30 dB. However at E /N b 0 = 60 dB, decision feedback does not improve the BER at all. The improvement in B E R obtained with 2-bit D F are about 40% and 30% at E /N b 0 = 30 dB and 60 dB, respectively. Figure 4.11 shows the average B E R in FS Rayleigh fading with x = 0.2T and CDR = 10 dB . It shows that irreducible error floors start to occur at about E /N =A0 b o dB. Similar to the previous fading cases described earlier, the simulation results agree with the 99 theoretical results for 1 and 2-bit 'conventional' and 1-bit DF. With 2-bit DF, the average B E R from simulations are about 30% and 25% higher than the corresponding theoretical results for E /N b at 30 dB and 60 dB, respectively. 0 10° I 0 i i i i 10 i 20 i i 30 i 40 i i 50 i I 60 Eb/N (dB) 0 Figure 4.10 Average BER as a function of E /N b in FS fading with CDR = 20 dB, x = 0.47* The improvement in B E R with 1-bit D F over 1-bit 'conventional' is about 30% at E /N b a = 30 dB . At E /N b 0 = 60 dB, decision feedback even degrades the B E R by about 20%. This is because the feedback coefficient, d,, that we assume throughout the simulations (d x = 1) is not optimal in this case. For 2-bit DF, the improvement in B E R over 2-bit 'conventional' is about 20% at E /N b a = 30 dB. At E /N b a = 60 dB, the B E R obtained with decision feedback is about 5% higher than that of 'conventional' though theoretical result shows a 25% improvement in BER. This is because the analysis does not take into account of error propagation, as explained earlier. 100 Figure 4.12 shows the average BER in FS fading at E /N b 0 from 40 dB to 60 dB for the same CDR and T as in Figure 4.11, but with new optimal feedback coefficients. In 1-bit DF, both theoretical and simulation results show that the optimal feedback coefficient, d x is equal to 0.5. In this case, 1-bit D F exhibits a slight improvement in BER when compared to 1-bit 'conventional'. With 2-bit DF, the optimal feedback coefficient , d , from theory is 2 whereas it is 1.5 13 2 from simulation. The difference is again due to the fact that 'error propagation' is not considered in the theoretical analysis. The simulated BER from 2-bit DF is approximately the same as that of 2-bit 'conventional' when the new coefficient is used. Figure 4.13 shows the average B E R in FS Rayleigh fading with x = OAT and 1 3 For both 1-bit DF and 2-bit DF, the optimum coefficient is chosen in steps of 0.25 from the original value; the step size beingfineenough to give an accurate optimum value. 101 • 1-bit conventional, s 1-bit DF, s • 2-bit conventional, s * 2-bit DF, s 1-bit conventional, t 1-bit DF, t 2-bit conventional, t 2-bit DF, t A DC 111 00 ro > < chosen to be 0.5, instead of 1 for the simulated 1-bit DF GMSK curve and d 2 is chosen to be 1.5, instead of 2 for the simulated 2-bit DF GMSK curve. CDR = 10 dB . In this case, irreducible error floors occur at around E /N b a = 30 dB, which is lower than the corresponding values for the other cases of fading previously described. There is good agreement between the simulation results and the theoretical results for 1 and 2-bit 'conventional' and 1-bit DF. In the case of 2-bit DF, the average B E R from simulation is about 30% higher than the theoretical results for both E /N b >30 dB. The improvement in B E R with 1-bit 0 DF over 1-bit 'conventional' is about 20% at E /N b at all when E /N b 0 0 = 30 dB. However, the BER is not improved = 60 dB. The improvement in BER obtained with 2-bit D F over 2-bit 'conven- tional' is about 25% at both E /N b 0 = 30 dB and 60 dB. This section presented BER results of differentially detected G M S K in Rayleigh fading. For 1 and 2-bit conventional G M S K and 1-bit D F G M S K , the simulation results are in good agreement with the theoretical results. In the case of 2-bit D F GMSK, the simulation BER values are higher than the theoretical results due to the fact that error propagation is not taken into account in the theoretical results. The improvement in BER using decision feedback is greatest in 102 20 30 E b / N (dB) 40 0 Figure 4.13 Average BER as a function of E /N b in FS fading with CDR = 10 dB, T = 0.47. flat fading. There is almost no BER improvement with decision feedback for the FS channel when CDR = 10 dB and x = 0.2T. In the different cases of fading that have been described, 2-bit conventional G M S K performs slightly better than 1-bit conventional G M S K when E /N b 0 < 25 d B . However, at a higher E /N , b 0 the latter performs better and has a lower irreduc- ible error floor. In the case of decision feedback, there is almost no difference in performance between 1-bit D F G M S K and 2-bit D F G M S K when E /N b E /N b Q 0 <30 dB. However, when > 30 d B , 1-bit D F G M S K has a better performance and a lower error floor. It is also important to note that the corresponding BER performances for 1-bit and 2-bit conventional G M S K and their D F schemes, degrade as we go from the case of flat fading shown in Figure 4.8 to FS fading with CDR = 10 dB and x = 0.4T shown in Figure 4.13. 103 4.5 BKER Performance in Rayleigh Fading The block lengths used in the simulation are N=255, 511 and 1023 bits. As in Section 4.2, f is chosen to be 32 Hz, the corresponding f T D D T B is the duration of a block. Since f T D B B are 0.51, 1.022 and 2.046, respectively; and is not much less than 1 in all of the cases, the signal cannot be considered to be constant over a block and therefore the environment is not slow fading with respect to B K E R performance in this circumstance. Figure 4.14(a), Figure 4.14(c), Figure 4.15(a), Figure 4.15(c) and Figure 4.16(a) show the B K E R performances of 1-bit 'conventional' and 1-bit DF for block lengths N=255, 511 and 1023 bits, in flat and FS Rayleigh fading with CDR = 20 dB, x = 0.27/; C D R = 20 dB, x = OAT, C D R = 10 dB, x = 0.2T and C D R = 10 dB, x = OAT, respectively. Figure 4.14(b), Figure 4.14(d), Figure 4.15(b), Figure 4.15(d) and Figure 4.16(b) show the corresponding B K E R performances of 2-bit 'conventional' and 2-bit D F for the same block lengths and fading parameters, respectively. In general, we see that irreducible error floors appear at E /N b Q in the range of 40 dB to 50 dB. With 1-bit detection: in Figure 4.14(a) for flat fading, Figure 4.14(c) and Figure 4.15(a) for FS fading where the C D R = 20 dB, the B K E R approximately doubles when N doubles at E /N b Q > 50 dB. In Figure 4.15(c) and Figure 4.16(a) for FS fading where the CDR = 10 dB, the B K E R increases with N as expected but the increase is less than double when TV doubles. In the case of 2-bit detection, the B K E R performances show similar trends, as in 1-bit detection. 104 Figure 4.14 Average BKER as a function of E /N b o for (a) 1-bit conventional GMSK and 1-bit DF GMSK in flat fading (b) 2-bit conventional GMSK and 2-bit DF GMSK inflatfading (c) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 20 dB, % = 0.2T (d) 2-bit conventional GMSK and 2-bit DF GMSK in FS fading with CDR = 20 dB , T = 0.2T. 105 Figure 4.14(a) shows the B K E R performance of 1-bit 'conventional' and 1-bit D F in flat fading. With DF, at E /N = b 60 dB, the improvement in B K E R over 1-bit 'conventional' is Q approximately 55% for the values of N considered. The percentage improvement is almost equal to that of the corresponding BER. Figure 4.14(c) shows the B K E R performance of 1-bit 'conventional' and 1-bit D F in FS Rayleigh fading with x = 0.2T and C D R = 20dB . The B K E R obtained from 'conventional' is improved by about 15% when D F is applied, for the values of N considered at E /N b a = 60 dB. The percentage improvement is roughly equal to that of the corresponding BER performance. Figure 4.15(a) shows the B K E R performance of 1-bit 'conventional' and 1-bit D F in FS Rayleigh fading with x = 0.4T and CDR = 20 dB. With decision feedback, the improvement in B K E R is about 40% at E /N b 0 = 30 dB but there is almost no improvement at E /N b Q = 60 dB, for the values of N considered. The corresponding BER performance shows a similar percentage in improvement. Figure 4.15(c) shows the B K E R performance of 1-bit 'conventional' and 1-bit D F in FS Rayleigh fading with x = 0.2T and CDR = 10 dB. With decision feedback, the improvement in B K E R is about 30% at E /N b = 30 dB for N = 255 bits and it is about 25% when N= 1023 bits. 0 However, the B K E R is slightly degraded at E /N b a = 60 dB, for all three values of N. As with the B E R performance, the degradation is due to the fact that the feedback coefficient has not been chosen appropriately at E /N b 0 = 60 dB. S> 10to CU > < — A 1-bit c, block length N=511bits E 10"' * 40 50 60 40 50 60 — • 2-bit c, block length N=255bits CD > < — • 1-bit c, block length N=255bits • B — 4 2-bit c, block length N=511 bits — • 2-bit c, block length N=1023bhs —•1-bit c, block length N=1023bits • - - B 2-bit DF, block length N=255bits - - • 1-bit DF, block length N=255brts | A - -A2-bitDF, block length N=511bits * - - *2-bit DF, block length N=1023bits - -A1-bitDF,block length N=511bits - -*1-bitDF, block length N=1023Ms J 10 j l 20 30 I L 50 i 40 10'' 60 10 20 Eb/N (dB) 30 Eb/N (dB) 0 0 (a) (b) 10 rr u cr LU LU CO 0) CD CO CD > < CO i CD 1 -bit c, block length N=255bils > 1 -bit c, block length N=511 bits «1-bit c, block length N=1023bits < B 1-bit DF, block length N=255bits — • 2-bit c. block length N=255bHs 10" lA - - A1-bit DF, block length N=511bits | * - - • 1-bit DF, block length N=1023bits — i 2-bit c, block length N=511btts — • 2-bit c. block length N=1023Ms • - - • 2-bit DF, block length N=255bits |A - - A 2 - M DF, block length N=511bils \ » - - * 2-bit DF, block length N=1023bits 10"' 10 20 30 40 50 10 20 30 Eb/N (dB) Eb/N (dB) 0 0 (C) Figure 4.15 60 (d) Average BKER as a function of E /N B Q for (a) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 20 dB, x = 0.42" (b) 2-bit conventional GMSK and 2-bit DF GMSK in FS fading with CDR = 20 dB , T = 0.4F (c) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 10 dB, x = Q.2T (d) 2-bit conventional GMSK and 2-bit DF GMSK in FS fading with CDR = 10 dB, x = 0.27/. 107 Figure 4.14(b) shows the BKER performance of 2-bit 'conventional' and 2-bit DF in flat fading. When DF is applied, the improvement in BKER over 2-bit 'conventional' is approximately 70% for the values of N considered. In the case of BER performance, the percentage improvement is only about 40%. Better improvement with DF shown by the BKER performance is due to the fact that with 2-bit DF, the bit errors tend to occur in pairs owing to error propagation whereas the number of blocks with an odd number of errors is reduced. Figure 4.14(d) shows the BKER performance of 2-bit 'conventional' and 2-bit DF in FS Rayleigh fading with x = 0.27/ and CDR = 20dB . At E /N = b Q 60 dB, the BKER obtained with DF is approximately a 50% improvement over that of 'conventional', for the values of N considered. In the case of BER performance, the improvement is only about 20%. Figure 4.15(b) shows the BKER performance of 2-bit 'conventional' and 2-bit DF in FS Rayleigh fading with x = OAT and CDR = 20 dB. With decision feedback, the percentage improvement in BKER is about 45% at E /N b Q = 60 dB, for the values of N considered. The corresponding percentage improvement in the BER performance is about 30%. Figure 4.15(d) shows the BKER performance of 2-bit 'conventional' and 2-bit DF in FS Rayleigh fading with x = 0.2T and CDR = 10 dB. When decision feedback is applied, the improvement in BKER is about 25% at E /N b 0 = 60 dB for N = 255 bits but it decreases to about 18% when N = 1023 bits. However, we have previously mentioned that the BER is slightly degraded (see Figure 4.11) because the decision feedback coefficient that is assumed in other fading cases, is not optimal in this case. The reason that there exists an improvement in BKER is due to error propagation, as described earlier. 108 Figure 4.16(a) shows the BKER performances of 1-bit 'conventional' and 1-bit DF in FS Rayleigh fading with x = OAT and CDR = 10 dB. At E /N = 60 dB, decision feedback does B 0 not improve the BKER, for all three values of N. This is similar to the trend shown by the BER performance. Figure 4.16(b) shows the corresponding BKER performance of 2-bit 'conventional' and 2-bit DF. When decision feedback is applied, the improvement in BKER is about 35% at E /N = 60 dB for N = 255 bits and it is about 15% when N = 1023 bits. B A 10 rr u rr LU LU DO CD O) CO m CD D> 2 cu 3 — • 2-bit c, block length N=255bits 3 — a 1-bit c, block length N=255bils — A 2-bit c, block length N=511bits — • 2-bit c, block length N=1023bils — * 1-bit c. block length N=511bils • - - •2-bit DF, block length N=255bits — • 1-bit c, block length N=1023bits | A - - A2-bH DF, block length N=511bits |Q - - • 1-bit DF, block length N=255bils - - * 2 - M DF, block length N=1023bits |A - - Al-bit DF, block length N=511 bits I* - - *1-brt DF, block length N=1023bits 10" 10 20 30 40 50 Eb/N (dB) 0 60 10" _1_ 10 20 _1_ 30 40 Eb/N (dB) 50 60 0 (b) (a) Figure 4.16 Average BKER as a function of E /N b for (a) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 10 dB, T = OAT (b) 2-bit conventional GMSK and 2-bit DF GMSK in FS fading with CDR = 10 dB , t = 0.47. Figure 4.17 and Figure 4.18 show the BKER performances of 1-bit and 2-bit 'conventional' for blocks lengths N= 255, 511 and 1023 bits in fiat and FS fading with the same parameters (CDR and x) as for Figure 4.14 to Figure 4.16. In all cases, 1-bit 'conventional' has a better 109 performance over 2-bit 'conventional' when E /N b lower irreducible error floor for high E /N , b 0 0 >25 dB because it has a lower B E R and a For flat and FS fading with CDR = 20 dB at E /N b a = 60 dB, the improvement is about the same for all values of N. For FS fading with C D R = 10 dB, the improvement drops as N increases. When the channel goes from flat fading to FS fading, the improvement also decreases for all values of N. For example, there is around 70% improvement with flat fading but it is down to about 50% for FS fading with x = OAT and CDR = 10 dB, for N = 255 bits. The cumulative density function (CDF) of the number of bit errors in a block for 1-bit and 2-bit detection in flat fading at E /N b 0 = 40 dB for N = 255, 511 and 1023 bits are shown in Figure 4.19 and Figure 4.21, respectively. Figure 4.19(a) compares the CDFs between 1-bit 'conventional' and 1-bit DF. It can be seen that the Pr{errors <M} for 1-bit 'conventional' is always lower than that for 1-bit DF. This is because the BER for 1-bit 'conventional' is always higher than that for 1-bit DF as shown in Figure 4.8. Figure 4.19(b) compares the CDFs between 2-bit 'conventional' and 2-bit DF. The C D F for 2-bit D F shows a stair-case-like shape: for example, the Prjerrors < 1} is just slightly higher (about 0.2% higher) than 7V{errors <0} , but the Pr{errors <2} is about 3.5% higher than Pr{ errors < 1} . This implies that the number of blocks that contain only one bit error is much less than those that have double-bit errors. The predomination of double-bit errors was discussed in Section 4.4. rr LU m 0 CD « i — • 1-bit c, block length N=255bits CD > < — * 1-bit c, block length N=S11 bits —•1-bit c, block length N=1023bits !• - - B 2-bit c, block length N=255bits |A - - A2-bit c, block length N=511bits | » - - » 2-bit c block length N=1023bits J 20 10 30 40 Eb/N (dB) 0 Average BKER as a function of E^/N' I I l_ 30 40 Eb/N (dB) (d) 50 0 (C) Figure 4.17 20 for 1-bit conventional GMSK and 2-bit conventional GMSK in (a) flat fading (b) FS fading with CDR = 20 dB, x = 0.27 (c) FS fading with CDR = 20 dB, x = 0.47 (d) FS fading with CDR = 10 dB, x = 0.2T. 60 Ill tr LU CO CD cn co * CD 3 20 30 40 Eb/N (dB) 0 Figure 4.18 Average BKER as a function of E^/N' for 1-bit conventional GMSK and 2-bit conventional GMSK in FS fading with CDR = 10 dB, t = 0.41. Generally, the probability that a block contains an odd number of errors is less than that of an even number of errors. This is clearly illustrated in Figure 4.20 which shows the CDFs at E /N b 0 = 20 dB . We see that in the C D F of 2-bit DF, 'flat tops' appear when M goes from even to odd and they are joined by 'slopes' when M changes from odd to even. Figure 4.21 shows the CDFs of 1-bit and 2-bit 'conventional' in flat fading for N = 255, 511 and 1023 bits at E /N b a = 40 dB. It can be seen that for M < 2, 1-bit 'conventional' shows a better performance than 2-bit 'conventional', for all three values of N. This is indicated by its higher 7 r{errors < M} over 2-bit 'conventional'. However, when M > 2, the performance of 1J bit 'conventional' is equivalent to or slightly worse than that of 2-bit 'conventional', although 1bit 'conventional' has a lower B E R . If we compare the simulated B K E R which is equal to 112 l-Pr{errors = 0} in Figure 4.21 and that obtained with (4.2.1) using the corresponding average BER in flat fading, we found that the simulated B K E R is lower than that obtained in (4.2.1) for 14 both 1-bit and 2-bit 'conventional' with all three values of N. This is because the assumption of independent bit errors made in (4.2.1) is not realistic in the case of fading, where error bursts occur during fades and hence correlation exists between bit errors. (a) Figure 4.19 CDF of the number of bit errors in a block at E^/N (b) o - 40 dB for (a) 1-bit conventional GMSK and 1-bit DF GMSK in flat fading (b) 2-bit conventional GMSK and 2-bit DF GMSK inflatfading. Both Figure 4.22 and Figure 4.23 compare the CDFs of 1-bit and 2-bit 'conventional' with their corresponding decision feedback schemes in FS fading, for the appropriate values of CDR and x, when N = 255, 511 and 1023 bits. Similar trends as in Figure 4.21 can be observed. As an example, for N= 255 bits, the BKER obtained from (4.2.1) is 0.08653 and 0.138468, for 1-bit and 2-bit 'conventional', respectively. Figure 4.21 CDF of the number of bits errors in a block at E^/N' = 40 dB for 1-bit and 2-bit conventional GMSK in flat fading 114 x 0.95 block length N=255bits x o block length N=511 bits * block length N=1023bits block length M=255bi1s o block length N=511bits * block length N=1023bits 2-bit c 2-bit DF f 0.9 1-bit c 1-bitDF £ 0.85 N 0.75 (a) (b) x block length N=255bits x o block length N=511bits o block length N=511bits * block length N=1023bits x block length N=1023bits block length N=255bits 2-bit c 2-bit DF 1-bit c 1-bitDF 0.65 15 M (c) Figure 4.22 20 25 30 0.3' 0 5 10 15 M 20 25 (d) CDF of the number of bit errors in a block at E^/N = 40 dB for (a) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 20 dB, x = 0.2T (b) 2-bit conventional GMSK and 2-bit DF GMSK in FS fading with CDR = 20 dB, x = 0.2T (c) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 20 dB, x = OAT (d) 2-bit conventional GMSK and 2-bit DF GMSK in FS fading with CDR = 20 dB, x = OAT 30 115 (d) (c) Figure 4.23 CDF of the number of bit errors in a block at E^/N o = 40 dB for (a) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 10 dB, X = 0.2T (b) 2-bit conventional GMSK and 2bit DF GMSK in FS fading with CDR = 10 dB, x = 0.2T (c) 1-bit conventional GMSK and 1-bit DF GMSK in FS fading with CDR = 10 dB , X = 0.47 (d) 2-bit conventional GMSK and 2-bit DF GMSK in FS fading with CDR = 10 dB, x = OAT. 116 Figure 4.24 compares the CDFs between 1-bit and 2-bit 'conventional' in FS fading (for CDR = 20 d B , x = 0.27/; C D R = 20 d B , x = 0.47/; C D R = 10 d B , x = 0.27 and 0.3 0 1 1 2 1 1 4 1 6 1 8 1 10 M 12 ' 1 14 1 16 1 18 0 1 20 1 0 (c) 1 10 1 20 1 30 1 M 40 1 50 1 60 (d) Figure 4.24 CDF of the number of bit errors in a block at E /N b Q = 40 dB for (a) 1-bit and 2-bit conventional GMSK in FS fading with CDR = 20 dB, x = 0.27 (b) 1-bit and 2-bit conventional GMSK in FS fading with CDR = 20 dB, x = 0.47 (c) 1-bit and 2-bit conventional GMSK in FS fading with CDR = 10 dB, x = 0.27 (d) 1-bit and 2-bit conventional GMSK in FS fading with CDR = 10 dB, x = 0.47. 70 117 CDR = 10 dB , x = 0.47") with N = 255, 511 and 1023 bits at E /N b a = 40 dB. In all of the cases, 1-bit 'conventional' performs better than 2-bit 'conventional' when M is small and the advantage disappears when M increases. This section examined the BKER performances of 1-bit and 2-bit differential detection of GMSK in flat and FS Rayleigh fading. In general, we see that irreducible errorfloorsappear at E /N b in the range of 40 dB to 50 dB. The percentage improvement in BKER of 1-bit DF 0 GMSK over 1-bit conventional GMSK is almost equal to that of the corresponding BERs presented in Section 4.4. In the case of 2-bit DF GMSK, the BKER performance shows an increased improvement over 2-bit conventional GMSK when compared to the corresponding BER performance. This is due to the presence of increased error propagation in 2-bit DF GMSK. In all cases, 1-bit conventional GMSK has a better BKER performance over 2-bit conventional GMSK when E /N b a >25 dB. This is because 1-bit conventional GMSK has a lower BER and a lower irreducible errorfloorin this E /N b 0 range. In the case of 2-bit DF GMSK, for a given block, its probability of containing an odd number of errors is less than that of an even number of errors, due to error propagation. When examining the CDFs of 1-bit and 2-bit conventional GMSK especially at E /N b 0 >20 dB, the simulated BKER is always lower than the corresponding BKER obtained with (4.2.1) which assumes independent bit errors. This is because fading causes the bit errors to occur in bursts, thus reducing the BKER. Chapter 5 Conclusions In this thesis, we examined the BER and B K E R performances of G M S K employing 1-bit and 2-bit differential detection in AWGN, flat and FS Rayleigh fading channels through analysis and computer simulations. A detailed analysis was carried out to look at the IF filtering effect on the BER performance over fading channels. The B K E R performance, which gives information on the bit error distributions, was mainly studied through simulations. Since the premodulation Gaussian filter at the G M S K transmitter introduces significant ISI, we also examined the performance of a decision feedback scheme which partially removes the effect of destructive ISI, over such channels. The effect of IF filtering on the B E R performance over Rayleigh fading channels was examined. For both 'filtered' and 'unfiltered' cases, the BER is almost insensitive to f T at low D E /N b but at high E /N , a b increases with f T. Q the B E R is dominated by the effect of the Doppler spread, and it A comparison of the 'filtered' and the 'unfiltered' BER in FS fading, shows D that the BER degradation due to IF filtering is negligible when f T D fT D for high E /N , b when E /N b 0 Q is small but it increases with whereas the BER degradation is almost the same for all values of is small. In the case of flat fading, it is observed that for the same E /N , b Q fT D the percentage difference between the 'filtered' and the 'unfiltered' B E R is almost constant, independent of f T. D The B E R results of differentially detected G M S K as a function of E /N b 0 in Rayleigh fading were then presented. For 1 and 2-bit conventional G M S K and 1-bit D F G M S K , the simula- 118 Chapter 5 Conclusions 119 tion results are in good agreement with the theoretical results. In the case of 2-bit D F G M S K , the simulation BER values are higher than the theoretical results. This is because error propagation is not taken into account in the analysis. The improvement in B E R using decision feedback is greatest in flat fading. In the different cases of fading that were considered, 2-bit conventional G M S K performs slightly better than 1-bit conventional G M S K for E /N b E /N , b Q a < 25 dB. For higher the latter has a lower irreducible error floor than the former. This is expected because the BER performances of 1-bit and 2-bit conventional G M S K (with a fixed zero decision threshold) in a non fading A W G N channel also shows a similar trend. In the case of decision feedback, similar performance trends are observed. As expected, the corresponding BER performances for 1-bit and 2-bit conventional G M S K and their D F counterparts, degrade as the CDR decreases. In FS fading for the same CDR, the BER performance also degrades when the time delay difference, x, is increased. For the purpose of computing the B K E R in a non fading A W G N channel, it was found that the bit errors can be assumed to be independent, though differential detection introduces some correlation between them. The B K E R results of 1-bit and 2-bit differential detection of G M S K in flat and FS Rayleigh fading show that irreducible error floors generally appear at E /N b 0 in the range of 40 dB to 50 dB. The percentage improvement in B K E R of 1-bit D F G M S K over 1-bit conventional G M S K is almost equal to that of the corresponding BERs at high E /N . b a In the case of 2-bit D F G M S K , the B K E R performance shows an increased improvement over 2-bit conventional G M S K when compared to the corresponding BER performance, due to the increased presence of error propagation in 2-bit D F GMSK. Moreover, due to error propagation, there is a higher probability that a given block contains an even number of errors than that of an odd Chapter 5 Conclusions 120 number of errors. Since fading gives rise to bursty bit errors, the actual B K E R is better than the predicted assuming independent bit errors. Some suggestions for related future research work are given as follows. • To derive formulae for B K E R performance in very slow fading, taking into account bit error correlation introduced by differential detection, and attempt to extend the analysis to the slow fading case. • Since G M S K and C P M in general have memory such that their information carrying phase follows a certain trellis path, Viterbi decoding can be applied for error correction after differential detection [6]. The BER and B K E R performances of differentially detected C P M can be analyzed with the application of Viterbi algorithm over the mobile radio channel. • To investigate the BER and B K E R performances in Nakagami fading used to model fading in microcellular systems [38]. Glossary ACI - Adjacent channel interference {dj} - Input information sequence of logical ' 1' and '0' AWGN - Additive White Gaussian Noise B - 3-dB bandwidth of the premodulation Gaussian filter B - 3-dB bandwidth of the IF Gaussian filter B - Equivalent noise bandwidth of the IF Gaussian filter {bj} - Non-return-to-zero (NRZ) signalling representation of {a^} P - Gain of the nth oscillator in the Jake's fading model t i n n BER - Bit Error Rate BKER - Block (Packet) Error Rate g c - Speed of light equal to 3 x 10 m/s C - Channel capacity {c .} - Differentially encoded sequence of {bj} CDF - Cumulative density function CDR - Average carrier (desired signal) to delayed signal power ratio CPM - Continuous Phase Modulation d and d - Decision feedback coefficients for 1-bit and 2-bit DF GMSK, respectively dj-(t) D(f) - Delayed (secondary) signal after fading - Doppler spectrum x 2 121 122 DF - Decision feedback A/ c - Coherence bandwidth A0 Jf l and A<|> ,; A<j> / and d> 4i A<jv, / - 'Unfiltered' differential phase angle (phase change) for the desired signal and the delayed signal, respectively. The 'tilde' denotes the 'filtered' case. U W Ac^ , .DF A ~,DF K ~DF and A ^ , ; A ^ / and A ^ / z - 'Unfiltered' differential phase angle (phase change) for the desired signal and the delayed signal after decision feedback, respectively. The 'tilde' denotes the 'filtered' case. - (Average) energy per bit or (average) energy per symbol, as specified fc a n d °>c - Peak frequency deviation fd f D - Carrier frequency in Hz and radians/s, respectively a n d D ffl - Doppler spread or maximum Doppler frequency in Hz and radians/s, respectively - Sampling frequency F(v) - Characteristic function of the sampled output V FM - Frequency Modulation FS - Frequency selective FSK - Frequency Shift Keying g(t) - Response of the premodulation filter to a transmitted rectangular pulse GMSK - Gaussian-filtered Minimum Shift Keying GPRS - General Packet Radio Service GSM - Global System for Mobile Communication 123 h - Modulation index n,y(r) - Impulse response of IF filter h(t) - Equivalent lowpass impulse response of IF filter h (t) - Impulse response of the premodulation filter Hi/if) - Frequency response of IF filter H(f) - Equivalent baseband frequency response of IF filter H (f) - Frequency response of the premodulation filter IF - Intermediate frequency ISI - Intersymbol Interference J (.) - Zeroth order Bessel function of the first kind A, - Decision feedback angle m - Number of bits per symbol in modulation MFSK - m-ary FSK MPSK - m-ary PSK MSK - Minimum Shift Keying n(t) - B andlimited noise N - Block length in bits t t a iV 0 N os - One-sided power spectral density of Gaussian noise - Number of low-frequency oscillators in the Jakes Model co - Frequency of the nth oscillator in the Jakes Model p.d.f. - Probability density function n 124 P - BER (for a non fading AWGN channel) or average BER (for fading) P(M, N) - Probability of more than M bit errors in a block of length N bits e $(t;b) or Q> (t) s - Transmitted information carrying phase § (t) - Filtered information carrying phase P - Desired (average) signal power at the point before A W G N is added PSD - Power spectral density or power spectrum PSK - Phase Shift Keying r(0 - 'Unfiltered' received signal r(t) - 'Filtered' received signal R - Bit rate s s b R - Symbol rate RF - Radio frequency s p,(.), p ( . ) a n d p „ ( . ) d - Normalized autocorrelation functions of the unmodulated carrier of the desired fading signal, the delayed fading signal and noise, respectively s(t) - Transmitted signal Sjr(t) - Desired (main or undelayed) signal after fading 2 G d 2 G n - Average power of dj(t) - Postdetection noise power, i.e. power of n(t) 125 C - Average power of sy(r) SNR - (Average) signal-to-noise ratio, it is an average value in fading t - Coherence time x - Time delay difference between the two signals in a FS channel s c % - Multipath delay spread T - Bit or symbol duration d T - Block duration UMTS - Universal Mobile Telecommunication Systems v - Velocity of the mobile unit V - Sampled output of lowpass filter in the receiver W - Channel bandwidth B W (t) - Equivalent baseband spectrum of noise W (t) - Equivalent baseband spectrum of the unmodulated carrier fading signal n s x (t), y (t), x (t), y {t), x (t) and y (t) s s d d n n - Independent zero mean, lowpass Gaussian processes Z(t) - Complex envelope of r(t) Z f(t) - Complex envelope of r(t) Z (t) - Equal to x (t) + jy (t) ZM) -Equal to Z (t) - Equal to x (t) + jy (t) t d s d d x (t)+jy (t) n s n s Bibliography [I] S. T. S. Chia and D. Grille-, "UMTS - mobile communications beyond the year 2000: Requirements, architecture and system options", IEE Electronics and Communication Engineering Journal, Vol. 4, pp. 331-340, Oct. 1992. [2] K. Pahlavan and A. H . Levesque, Wireless Information Networks, New York: John Wiley & Sons, Inc., 1995. [3] M . Mouly and M - B . Pautet, "Current Evolution of the G S M Systems", IEEE Personal Communications, Vol. 2, pp. 9-19, Oct. 1995. [4] M . K. Simon and C. C. Wang, "Differential versus Limiter-Discriminator Detection of Narrowband F M " , IEEE Trans, on Communications, Vol. 31, No. 11, pp. 1227-1234, Nov. 1983. [5] M . K. Simon and C. C. Wang, "Differential Detection of Gaussian M S K in a Mobile Radio Environment", IEEE Trans, on Vehicular Technology, Vol. 33, No. 4, pp. 307-320, Nov. 1984. [6] J. G. Proakis, Digital Communications, Third edition, New York: McGraw-Hill, Inc., 1995. [7] B. Sklar, "Defining, Designing, and Evaluating Digital Communication Systems", IEEE Communications Magazine, Vol. 31, pp. 91-101, Nov. 1993. [8] C - E Sundberg, "Continuous Phase Modulation, IEEE Communications Magazine, Vol. 24, No. 4, pp. 25-38, April 1986. [9] S. M . Elnoubi, "Analysis of G M S K with Differential Detection in Land Mobile Radio Channels", IEEE Trans, on Vehicular Technology, Vol. 35, No. 4, pp. 162-167, Nov. 1986. [10] S. M . Elnoubi, "Analysis of G M S K with Two-bit Differential Detection in Land Mobile Radio Channels", IEEE Trans, on Communications, Vol. 35, No. 2, pp. 237-240, Feb. 1987. [II] S. M . Elnoubi, "Analysis of G M S K with Discriminator Detection in Mobile Radio Channels", IEEE Trans, on Vehicular Technology, Vol. 35, No. 2, pp. 71-76, May 1986. [12] P. P. S. Carter and A . M . D. Turkmani, "Performance Evaluation of Rayleigh and Lognormal G M S K signals in the presence of Co-channel Interference", IEE Proceedings-Part I, Vol. 139, No. 2, pp. 156-164, April 1992. [13] I. Korn, "GMSK with Differential Phase Detection in the Satellite Mobile Channel", IEEE Trans, on Communications, Vol. 38, No. 11, pp. 1980-1986, Nov. 1990. 126 127 [14] I. Korn, "GMSK with Frequency Selective Rayleigh Fading and Co-channel Interference", IEEE Journal on Selected Areas in Communications, Vol. 10, No. 3, pp. 506-515, April 1992. [15] I. Korn, "GMSK with Limiter Discriminator Detection in Satellite Mobile Channel", IEEE Trans, on Communications, Vol. 39, No. 1, pp. 94-101, Jan. 1991. [16] S. S. Shin and P. T. Mathiopoulos, "Differentially Detected G M S K Signals in CCI Channels for Mobile Cellular Telecommunication Systems", IEEE Trans, on Vehicular Technology, Vol. 42, No. 3, pp. 289-293, August 1993. [17] P. Varshney and S. Kumar, "Performance of G M S K in a Land Mobile Radio Channel", IEEE Trans, on Vehicular Technology, Vol. 40, No. 3, pp. 607-614, August 1991. [18] P. Varshney, J. E . Salt and S. Kumar, "BER Analysis of G M S K with Differential Detection in a Land Mobile Channel", IEEE Trans, on Vehicular Technology, Vol. 42, No. 4, pp. 683689, Nov. 1993. [19] C. Tellambura and V. K. Bhargava, "Performance of G M S K in Frequency selective Rayleigh Fading and Multiple Cochannel Interferers", IEEE 45th Vehicular Technology Conference, Count Down to the Wireless Twenty-First Century, Chicago, IL, USA. pp. 211-215, July 1995. [20] J. Horikoshi and S. Shimura, "Multipath Distortion of Differentially Encoded G M S K with 2-b Differential Detection in Bandlimited Frequency-Selective Mobile Radio Channel", IEEE Trans, on Vehicular Technology, Vol. 39, No. 4, pp. 308-315, Nov. 1990. [21] K. Hirade, M . Ishizuka, F. Adachi and K. Ohtani, "Error-Rate Performance of Digital F M with Differential Detection in Land Mobile Radio Channels", IEEE Trans, on Vehicular Technology, Vol. 28, No. 3, pp. 204-212, August 1979. [22] A. Yongacoglu, D. Makrakis and K. Feher, "Differential Detection of G M S K using Decision Feedback", IEEE Trans, on Communications, Vol. 36, No. 6, pp. 641-649, June 1988. [23] S. Elnoubi, "Comments on " G M S K with Differential Phase Detection in the Satellite Mobile Channel"" , IEEE Transactions on Communications, Vol. 40, No. 4, pp. 666-669, April 1992. [24] R. R. Anderson, W. R. Bennett, J. R. Davey and J. Salz, "Differential Detection of Binary F M " , The Bell System Technical Journal, Vol. 44, pp. 111-159, Jan. 1965. [25] B. P. Lathi, Modern Digital and Analog Communication Systems, New York: Holt, Rinehart and Winston, 1983. [26] S. Haykin, Communication Systems, Second Edition, New York: John Wiley & Sons, 1983. 128 [27] K. Murota and H. Hirade, "GMSK Modulation for Digital Mobile Radio Telephony", IEEE Trans, on Communications, Vol. 29, No. 7, pp. 1044-1050, July 1981. [28] K. Murota, "Spectrum Efficiency of G M S K Land Mobile Radio", IEEE Trans, on Vehicular Technology, Vol. 34, No. 2, pp. 69-75, May 1985. [29] W. C. Y. Lee, Mobile Cellular Telecommunications Systems, New York: McGraw Hill Inc., 1990. [30] C. A . Siller, Jr., "Multipath Propagation", IEEE Communications Magazine, Vol. 22, No. 2, pp. 6-15, Feb. 1984. [31] R. H . Clarke, "A Statistical Theory of Mobile-Radio Reception", Bell System Technical Journal, Vol. 47, pp. 957-1000, July-August 1968. [32] W. C. Jakes, Jr., Microwave Mobile Communications, New York: John Wiley & Sons Inc., 1974. [33] F. Davarian, "Mobile Digital Communication via Tone Calibration", IEEE Trans, on Vehicular Technology, Vol. 36, No. 2, pp. 55-61, May 1987. [34] M . K. Simon, "Dual-pilot Tone Calibration Technique", IEEE Technology, Vol. 35, No. 2, pp. 63-70, May 1986. Trans, on Vehicular [35] R. Srinivasan, "Feeback Communications over Fading Channels", IEEE Communications, Vol. 29, No. 1, pp. 50-57, Jan. 1981. Trans, on [36] L . Hanzo, R. Steele and P-M. Fortune, "A Subband coding, B C H coding and 16-QAM System for Mobile Radio Speech Communications", IEEE Trans, on Vehicular Technology, Vol. 39, No. 4, pp. 327-339, Nov. 1990. [37] J. K . Cavers, "An Analysis of Pilot Symbol Assisted Modulation for Rayleigh Fading Channels", IEEE Trans, on Vehicular Technology, Vol. 40, No. 4, pp. 686-693, Nov. 1991. [38] J. Andersen, T. Rappaport and S. Yoshida, "Propagation Measurements and Models for Wireless Communications Channels", IEEE Communications Magazine, Vol. 33, pp. 42-49, Jan. 1995. [39] M . Schwartz, Communication Systems and Techniques, New York: McGraw-Hill, Inc., 1966. [40] N . Mohanty, Random Signals Estimation and Identification (Analysis and Applications), New York: Van Nostrand Reinhold Company, 1986. [41] E . Kreyszig, Advanced Engineering Mathematics 6th Edition, New York: John Wiley & Sons, Inc., 1988. 129 [42] B. Sklar, "Rayleigh Fading Channels in Mobile Digital Communication Systems: Part I and Part IT, IEEE Communications Magazine, Vol. 35, pp. 90-111, July 1997. [43] M . C. Jeruchim, "Techniques for Estimating the Bit Error rate in the Simulation of Digital Communication Systems", IEEE Journal on Selected Areas in Communications, Vol. S A C - 2 , N o . l,pp. 153-170, Jan. 1984. [44] W. H. Press, S. A . Teukolsky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in C: The Art of Scientific Computing 2nd Edition, Cambridge: Cambridge University Press, 1997. [45] J. A . Roberts, "Packet Error Rates for DPSK and Differentially Encoded Coherent BPSK", IEEE Transactions on Communications, Vol. 42, No. 2/3/4, pp. 1441-1444, Feb./March/ April 1994. Appendix A Evaluation of the Characteristic function of the Detector output The problem to be considered is to find the p.d.f. of the detector output V(IT) (denoted as V hereafter) through its characteristic function F(v) [39] F(v) = [ ^ f ( I - 2 y v C * Q ) ] _1 (A.l) where Q is the 2x2 Hermitian matrix of the quadratic form and it satisfies V = \z *QZ. T (A.2) 1-Bit Detection For 1-bit detection, Vis given by V = \Re{-jz z *} x (A.3) 2 where z\ = Z(IT) and z = Z((l-1)T). V can be readily shown to be Hermitian since Q is Hermitian. 2 Equating (A.2) with (A.3) and ignoring the scalar factor, Q can be obtained as Q = 0 j -j q 130 . (A.4) Appendix A 131 The covariance matrix C is given in (3.1.31). Therefore F(v) is obtained as z 2 F(v) = det 1 + 2VG p -l 2 -2va 2 2 (A.5) p 1 - 2VG 2VG = [l-2va (p-p*)+4vV(l-|p| )] 2 2 1 The derivation of the covariance matrix C for filtered 1-bit detection in flat and FS Rayleigh z fading channels is given in Appendix C. 2-Bit Detection Using the above procedure, Q can be obtained as Q = 0 1 (A.6) 1 0 for the 2-bit case. The corresponding characteristic function F(v) is found to be 2 F(v) = det 1-2jvG p -2JVG 2 2 -2jvo 2 1-2JVG p (A.7) = [ 1 - 2;va (p + p*) + 4 v V ( 1 - |p| )f 2 2 1 Appendix B 132 Appendix B Evaluation of the p.d.f. of the 1-bit detector output from its Characteristic function The residue theorem can be applied to evaluate the following improper integral, as a contour integral. The contour integral is being taken in the anti-clockwise sense around a closed path [41]. Therefore, the characteristic function F(v) can be taken to be |F(v ) e ->^v = j ' -'' ''" " l-2nj(res[F(z)e 2 I where the residue of a simple pole at z =PQ r« The ( [ F ( z ) e J f t ] ) ]) jVz V ° V>0 < <B.l) for a function f(z), is defined to be / ( z ) = lim (z-p )f(z) 0 (B.2) characteristic function F(v) for one-bit detection is derived in Appendix A and is given by (A.5). The two poles of F(v) are found to be P *> ~ 2 2 , ,2 2a(l-|p| ) 2 (B.3) Appendix B Consider the case when V < 0, res[F(z)e ] jVz lim (z-p )F(z)e j V z 1 2 ^ Pi lim z^ (z-p ){z-p ) Pl l 1 4 2 4 G 1 2 ( 1 e~ jVz = *lim -^>4o- (l-|p| ) 4 _|p| ) (z-p ) 2 2 j ; f J ~P r Pj J exp<{ — - — - i - V ^ 4a ^T^ l2a (l-|p| ) = l 2 + J ! 2 2 2 J Therefore the conditional p.d.f. of V when V < 0, is given by P(V\b = ± 1 ) = • l J F(v)e ; v V Jv ^•27U7-r 5[F(z)e" ' ] 7 Vz E [ 7 -Pr 1 4a /rrp~ 2 / where is the /r/z input bit. 1 2eXP + Pj l2a (l-|p| ) 2 2 Appendix B 134 For the case V > 0, the analysis is similar. res[F(z)e = ] jVz lim (z- p )F(z)e -jVz 2 (B.6) -jVz = lim ^^4a (l-|p| ) 4 2 (z-P\) V 2a (l-|p| ) 2 2 Therefore, the conditional p.d.f. of Vwhen V > 0, is given by P(V\b l = ±l) = dv ^z\F(v) = ~ (B.7) • 2nj • res[F(z)e- ] jVz 4a JTp" 2 A/1 - P r - P j ' exp-! 2 ^ 2o (i-ipi ); 2 2 V V>0 Appendix C 135 Appendix C Derivation of the Covariance matrix C for Filtered 1-Bit Detection in Flat and FS Rayleigh z Fading Channels Flat Rayleigh Fading Channel The filtered received waveform at the receiver output r(t) is given as r(t) = X(r)cosco r-Y(t) sin® t c (C.l) c where X(t) = x (t)A (t) - y (t)A (t) + x {t), n (C.2) Y(t) = x (t)A (t) + y (t)A (t) + y (t), (C.3) s s c s s s s c n and A (t) = cos(|) (0® h (t), c if s , (CA) A (t) = sm$ (t)®h (t). s s (C.5) if T We define a vector a = [x (t), y (t), x (t), y (t)\, s s n n where x (t), y (t), x (t) and y (t) are s s n n mutually independent zero mean Gaussian processes. Thus, the following relations are obtained: Appendix C 136 x (t) = y (t) s s = x (t) = y {t) n n = 0 (Oy,-(0 = y,-(0*,-(0 = *,-(0*,-(0 = y,-(0y,-(0 = o (Oy/0 = y,-(0*/0 = o (C.6) x (f + $)y,(0 = y,-(f + ^,.(0 = 0 where i , j = s, n j where the bar denotes taking the ensemble average. Hence, the covariance matrix of a is given by _ C a = — T [a-a][a-a] (C.7) a! 0 0 0 0 a; 0 0 0 0 a„ 0 0 0 0 a; 2 2 where a and a are the average received powers of the signal after fading and noise. Define a vector e = [X Y , X , Y ], where s M v x 2 2 Xj(0 = X ( 0 , r , ( 0 = Y(t),X (t) 2 = X(t-T),Y (t) 2 = Y(t-T) (C.8) and X- = X (lT) t , F- = y,.(/r) , ( i = 1,2;/is an integer). (C.9) Appendix C 137 Its covariance matrix H is defined as e H = e (CIO) [ e - e ] [ e - e f where the bar denotes taking the statistical average. H is computed to be e H = e aT 0 0 a 0"lPl 2 -c p a 2p2 a lPl a pj 2 2 o"iPi " 2 p 2 a °2P2 2 2 2 a 0 (C.ll) 0 a 2 where 2 °" a 2 2 .22 2 .22 2 (C.12) (C.13) = 2°,+<V A i P i = i 2 ?P,(7 )cos(A^ ) + A G A , ; / °1>P2 = -^1^2 ?P,( ) a r s i n a p„(r), 2 (C.14) (C.15) ( ^,/)' A A (t) = A (t) + A (t) , 2 A<j> = 5(/ and 2 (C.l 6) 2 tan" (A (/r)/A (/r))-tan" (A,(/r-r)/A (/r-7')) 1 1 5 c c ) (C.17) Appendix C 138 A] = A (IT-T) (C.18) 2 A\ = A (IT) 2 Define the vector Z to be Z = [Z(IT) Z(IT-T)] (C.19) where Z(t) = Xi(t) + yYi(t) and its covariance matrix is defined as C = z -[Z-Zf[Z-Z] X (C.20) T The elements of C can be computed from H as e z G p +jO 0 2 2 1 l 2 O^P 2 a,a p a 2 lPl -7" 2p2 CT a 2 (C.21) 2 where 0"jO p is given as 2 a,a p = o p - ; a p 2 2 2 1 2 = A j A ^ p ^ e x p C / A ^ / ) + CT p (T) 2 n (C.22) F S Rayleigh Fading Channel In the case of FS fading, the filtered received waveform at the receiver output is given in (C.l), where X(t) = x (t)A (t)-y (t)A (t) + x (t- x)A (t -x)-y (ts d c s c (C.23) s d x)A (t -x) + x (t) s n Appendix C 139 Y(t) = x (t)A (t) + y (t)A (t) + x (t- x)A (t - x) + y (t- x)A (t -x) + y (t) s s s d c s d c . (C.24) n T We define a vector a = [x , y , x , y , x , y ]. Following the same procedure as in the s s d d n n filtered FS fading case and using the relations given in (C.6) in which i,j now denote the subscripts s,d or n and i ^ j, we can obtain the corresponding covariance matrix H as given in e (C.l 1) where Oj = A o +A G 1 xl s +c d , n (C.25) <h = Ws + \i^d + «. A a 2 k = A A a p,(7')cos(Ac> ) + A , A , a J 2 P l 1 ap 2 2 J;/ 1 2 P r f = -[A A CT p (r)sin(A^ ) + A A 1 5 2 i/ T l c 26 (r)cos(A^ ) + a ; 2 2 (- ) a T 2 2 p„(r), a^(r)sin(A^ )], ) / (C.27) (C.28) with Acj)^ / being given in (C.17) and A<j>4/ = tan (A (IT-x)/A {IT-x)) -tan {A (IT- T-x)/A {lT- X S 1 C S c T-x)). (C.29) Equation (C.l 8) defines the variables A j and A and 2 A ! = A (lT-T 2 A 2 = A (/7 + 2 2 2 + x) . ) (C30 T) The corresponding covariance matrix C can be similarly obtained as that in thefilteredflatfading z 2 2 case. C is given in (C.20) with a , G given by (C.25) and (C.26), respectively and z x 2 Appendix C 140 2 2 o,a p = o , p - 7 a P 2 2 = A 1 l 2 5 p*( ) A a 2 r _ 2 e X P O ' < k /) + x\ z2°dPd( ) V(J $d, A A A T ex A l) + C p (T) 2 n n (C.3i) Appendix D 141 Appendix D Differential Phase Angles of 1-bit and 2-bit Detection BT t 0.25 9-3 0-2 0-1 0.0009 0.6 18.2 0 01 O 52.4 02 18.2 0.6 0. 0.0009 Phase shifts 6^ _ . (in degrees) for B,T = 0.25 in 1-bit detection. BT t 0.25 0-3 0-2 0-1 00 01 02 03 0.00017 0.2 11.3 50.3 26.7 1.5 0.0042 Table D.2 Delayed phase shifts (in degrees) for B T= 0.25 with x - 0.27 in 1-bit detection. t Bit Combinations State ''/-I h l /+/ /; 1 1 1 1 88.8 88.3 1 1 -1 2 52.4 34.9 -1 1 1 3 52.4 65.7 -1 1 -1 4 16.0 12.3 1 -1 1 5 -16.0 -12.3 1 -1 -1 6 -52.4 -65.7 -1 -1 1 7 -52.4 -34.9 -1 -1 -1 8 -88.8 -88.3 Table D.3 Differential phase angles A§ ^ and Aty^ ^ (in degrees) of the one-bit detector corresponding to various input data combinations for B T- 0.25. The contributions of and b[ are ignored. s t +2 Appendix D 142 0.25 ¥-3 ¥-2 ¥-1 ¥ 0.00095 0.6 18.8 70.6 0 ¥l ¥ 70.6 18.8 0.6 2 • ¥ 3 " Table D.4 Phase shifts \|/ _ • (in degrees) for B,T - 0.25 in two-bit detection. ; B,T 0.25 ¥-3-' ¥-2 _ ¥-i ¥o '. ¥i ¥2 ¥3 0.0002 0.2 11.5 61.5 77.0 28.2 1.5 Table D.5 Delayed phase shift _ • (: degrees) for B,T - 0.25 with x = 0.2 T , in two-bit detection. Suic Bit Combinations ' l-2 l-l c l c 1 1 1 - 1 l+l c - 1 1 i c 1 1 1 . i - . 9 37.6 24.2 10 37.6 55.2 11 0.0 -32.2 . - 1 1 - 1 12 0.0 -1.2 - 1 1 - 1 1 13 0.0 1.2 - 1 - 1 1 1 14 0.0 32.2 .1 1 -1 -1 15 -37.6 -55.2 -1 -1 1 -1 16 -37.6 -24.2 17 -103.6 -98.8 1 - . - 1 1 . -. -1 -1 18 -141.2 -155.2 .. -1 -1 1 19 -141.2 -121.8 -1 -1 -1 -1 20 -178.8 -178.2 21 178.8 178.2 1 1 1 1 . 1 1 - 1 22 141.2 121.8 - 1 1 1 1 23 141.2 155.2 24 103.6 98.8 - 1 1 1 . 1 Table D.6 Differential phase angles A^^ and A())^ (in degrees) of the two-bit detector corresponding to various input data combinations for B T- 0.25. The contributions of C/_3 and c; are ignored. ; ; t +2
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Performance of differentially detected GMSK in a mobile...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Performance of differentially detected GMSK in a mobile radio channel Mok, William Chiu Hung 1998
pdf
Page Metadata
Item Metadata
Title | Performance of differentially detected GMSK in a mobile radio channel |
Creator |
Mok, William Chiu Hung |
Date Issued | 1998 |
Description | As the demand for current cellular voice-oriented networks to include other services such as data and video applications is increasing rapidly, the reliability of information transmission over the mobile radio environment is important. Performance evaluation in terms of the bit error rate (BER) and block (packet) error rate (BKER) of modulation schemes in such an environment is of considerable research interest. In this thesis, we study the BER as well as the BKER performance of Gaussian-filtered Minimum Shift Keying (GMSK) modulation, currently used in the Global System for Mobile Communication (GSM) standard, in additive white Gaussian noise (AWGN) as well as flat and frequency selective Rayleigh fading. 1-bit and 2-bit differential detectors are used for demodulation. A detailed analysis is carried out to examine the IF filtering effect on the BER performance over the fading channels. The BKER performance is mainly studied through simulations. Since the premodulation Gaussian filter at the GMSK transmitter introduces significant Intersymbol Interference (ISI), we also examine the performance of a decision feedback (DF) scheme which partially removes the ISI, over such channels. The BER results in both flat and frequency selective (FS) Rayleigh fading show that 2-bit conventional (with no decision feedback) GMSK performs slightly better than 1-bit conventional GMSK when Eb/NQ < 25 dB. Above this value, the latter performs better and has a lower irreducible error. Results also show that 2-bit conventional GMSK is more sensitive to the IF filtering distortion than 1-bit conventional GMSK in Rayleigh fading. The percentage improvement in BKER of 1-bit DF GMSK over 1-bit conventional GMSK is almost equal to that of the corresponding BERs. In the case of 2-bit DF GMSK, the BKER performance shows an increased improvement over 2-bit conventional GMSK when compared to the corresponding BER performance, owing to increased error propagation. Since fading gives rise to bursty bit errors, the actual BKER is lower than the predicted assuming independent bit errors. In all fading cases, 1-bit conventional GMSK has a better BKER performance over 2-bit conventional GMSK when Eb/NQ >25 dB. In the case of decision feedback, similar performance trends are observed. |
Extent | 5090233 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-05-19 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0065240 |
URI | http://hdl.handle.net/2429/7962 |
Degree |
Master of Applied Science - MASc |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of Electrical and Computer Engineering, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 1998-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_1998-0278.pdf [ 4.85MB ]
- Metadata
- JSON: 831-1.0065240.json
- JSON-LD: 831-1.0065240-ld.json
- RDF/XML (Pretty): 831-1.0065240-rdf.xml
- RDF/JSON: 831-1.0065240-rdf.json
- Turtle: 831-1.0065240-turtle.txt
- N-Triples: 831-1.0065240-rdf-ntriples.txt
- Original Record: 831-1.0065240-source.json
- Full Text
- 831-1.0065240-fulltext.txt
- Citation
- 831-1.0065240.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0065240/manifest