On the Statistical Properties of Deterministic ... - Semantic Scholar

14 downloads 0 Views 580KB Size Report
kinds of distribution functions, such as the Nakagami and Weibull distributions. ... where the quantities cn, fn; and n are called Doppler coe cients, discrete ...
On the Statistical Properties of Deterministic Simulation Models for Mobile Fading Channels MATTHIAS PA TZOLD, MEMBER, IEEE, ULRICH KILLAT, MEMBER, IEEE, FRANK LAUE, YINGCHUN LI Technical University of Hamburg-Harburg Department of Digital Communication Systems Postfach 90 10 52 Denickestr. 17 D-21071 Hamburg Tel. : +4940-7718-3149 Fax : +4940-7718-2941 E-Mail : [email protected]

Abstract

Rice's sum of sinusoids can be used for an ecient approximation of coloured Gaussian noise processes and is therefore of great importance to the software and hardware realization of mobile fading channel models. Although several methods can be found in the literature for the computation of the parameters characterizing a proper set of sinusoids, only less is reported about the statistical properties of the resulting (deterministic) simulation model. This is the topic of the present paper; whereby not only the simulation model's amplitude and phase probability density function will be investigated, but also higher order statistics (e.g. level-crossing rate and average duration of fades). It will be shown that due to the deterministic nature of the simulation model analytical expressions for the probability density function of amplitude and phase, autocorrelation function, level-crossing rate, and average duration of fades can be derived. We also propose a new procedure for the determination of an optimal set of sinusoids, i.e., the method results for a given number of sinusoids in an optimal approximation of Gaussian, Rayleigh, and Rice processes with given Doppler power spectral density properties. It will be shown that the new method can easily be applied to the approximation of various other kinds of distribution functions, such as the Nakagami and Weibull distributions. Moreover, a quasi optimal parameter computation method is presented.

Submitted to

IEEE Transactions on Vehicular Technology September 1995 Final Version 1

I. INTRODUCTION In general, simulation models for mobile radio channels are realized by employing at least two or more coloured Gaussian noise processes. For instance, for the realization of a Rayleigh or Rice process two real coloured Gaussian noise processes are required; whereas the realization of a Suzuki process [1], which is a product process of a Rayleigh and a lognormal process, is based on three real coloured Gaussian noise processes. Such processes (Rayleigh, Rice, Suzuki) are often used as appropriate stochastical models for describing the fading behaviour of the envelope of the signal received over land mobile radio channels, which are often classi ed as frequency non-selective channels. A further example is given by modelling a (n-path) frequency selective mobile radio propagation channel by using the n-tap delay line model [2]. This requires the realization of 2n real coloured Gaussian noise processes. All these examples show us that an ecient design method for the realization of coloured Gaussian noise processes is of particular importance in the area of mobile radio channel modelling. A well known method for the design of coloured Gaussian noise processes is to shape a white Gaussian noise (WGN) process by means of a lter that has a transfer function, which is equal to the square root of the Doppler power spectral density of the fading process (see Fig. 1(a)). Another method, which is recently coming more and more in use, is based on Rice's sum of sinusoids [3, 4]. In this case a coloured Gaussian noise process is approximated by a nite sum of weighted and properly designed sinusoids. The resulting signal ow diagram of such an approximated Gaussian noise process is shown in Fig. 1(b), where the quantities cn , fn ; and n are called Doppler coecients, discrete Doppler frequencies, and Doppler phases, respectively, and N denotes the number of sinusoids. In the meanwhile, diverse methods have been developed for the derivation of the relevant model parameters (Doppler coecients cn and discrete Doppler frequencies fn). For example, the method of equal distances [5, 6] and the mean square error method [6]. The characteristic of both methods is just the same as Rice's original method [3, 4], namely, the di erence between two adjacent discrete Doppler frequencies is equidistant; but the Doppler coecients are adapted to a given Doppler power spectral density in a method speci c manner. The disadvantage of these methods is that (due to the equidistant discrete Doppler frequencies) the approximated coloured Gaussian noise process ~(t) is periodical. Another method is known as the method of equal areas [5,6]. That procedure results in a relativly satisfactory approximation of the desired statistics for the Jakes Doppler power spectral density even for a small number of sinusoids but fails or requires a comparatively large number of sinusoids for other types of Doppler power spectral densities such as those with Gaussian shapes. Widespread in use is the so called Monte Carlo method [7,8]. But unfortunately the performance of that procedure is poor [9]. The reason is that the discrete Doppler frequencies are random variables and therefore, the moments of the Doppler power spectral density of the designed Gaussian noise process ~(t) are also random variables, which can di er widely from the desired (ideal) moments even for a large number of sinusoids N , say N = 100 [6,10]. A further design method was derived by Jakes [11]. For the Jakes method an ecient hardware implementation, using the Intel 8088 microprocessor, is described in [12]. Although all the parameter computation methods can be classi ed in deterministic and statistic procedures, the resulting simulation model is in any case a deterministic one; because all prameters of the simulation model have to be determined during the simulation set up phase and thus, they are known and constant quantities during the total simula2

tion run phase. In this paper, the statistical properties of such deterministic simulation models will be investigated. The paper recapitulates in Section II the statistical properties of Rice processes with given spectral characteristics for the underlying Gaussian processes. Section III introduces the concept of deterministic simulation models and discusses the performance of three di erent parameter computation methods, whereby two of them are new. One of the new procedures is for a given number of sinusoids optimal in the following sense: First, the Doppler coecients cn are calculated such that the probability density of the approximated coloured Gaussian process ~(t) approximates optimally the aspired Gaussian probability density function; and second, the discrete Doppler frequencies fn are optimized so that the autocorrelation function of ~(t) approximates optimally the desired autocorrelation function of a given coloured Gaussian noise process. Moreover, analytical expressions for the deterministic simulation model's amplitude and phase probability density function, autocorrelation function, and cross-correlation function will be derived. Section IV investigates higher order statistical properties of deterministic simulation models such as the level-crossing rate, average duration of fades, and the probability density function of the fading intervals. It is shown analytically and experimentally that the level-crossing rate of properly designed and dimensioned deterministic simulation models for Rayleigh and Rice processes is extremely close to the theory even if the number of sinusoids is small, say seven or eight. Section V concludes the paper.

3

II. A STOCHASTIC ANALYTICAL MODEL FOR RICE PROCESSES In order to compare the statistical properties of the deterministic simulation model as will be introduced in Section III with the desired (exact) statistics, a reference model is required. As reference model, we consider in this section an analytical model for Rice processes. In cases, where the in uence of a direct line of sight component can not be neglected, Rice processes are often used as suitable stochastical models for characterizing the fading behaviour of the received signal level of land mobile terrestrial channels as well as of land mobile satellite channels.

A. Description of the Analytical Model A Rice process, (t), is de ned by taking the absolute value of a non-zero mean complex Gaussian process

(t) = (t) + m(t)

(1)

according to

(t) = j(t)j :

(2)

In (1) all the scattered components in the received signal are represented by a zero mean complex Gaussian noise process

(t) = 1 (t) + j2(t)

(3)

with uncorrelated real components i (t); i = 1; 2; and variance V arf(t)g = 2V arfi(t)g = 202, whereas the in uence of a direct (line of sight) component is taken into account here by a time variant mean value of the form

m(t) = m1(t) + jm2(t) = ej(2ft+ ) :

(4)

Amplitude, Doppler frequency, and phase of the direct component are denoted by ; f, and , respectively. Observe that for f = 0 the mean value m(t) = m = ej is obviously time invariant, which directly corresponds to an orthogonality between the direction of the incoming direct wave component at the receiving antenna of the vehicle and the direction of the vehicle movement. Throughout the paper, we consider the more general case where f 6= 0 and hence, the direct component m(t) depends on time. A typical and often assumed shape for the Doppler power spectral density (psd) of the complex Gaussian noise process (t); S(f ), is for mobile fading channel models given by the Jakes psd [11]

8 < fmaxp12?02(f=fmax)2 ; jf j  fmax ; S(f ) = : 0 ; jf j > fmax ; 4

(5)

where fmax is the maximum Doppler frequency. From this equation it follows for the corresponding autocorrelation function (acf) r (t), i.e. the inverse Fourier transform of S (f ), the following relation

r(t) = 202J0(2fmaxt) ; (6) where Jo() denotes the zeroth order Bessel function of the rst kind. The resulting structure of the analytical model for Rice processes (t) with underlying Jakes psd characteristic is shown in Fig. 2(a). Another important shape for the Doppler psd S (f ) is given by the frequency shifted Gaussian psd

r





2

2

2 f ?pfs (7) S(f ) = 2f0 ln2 e? fc= ln 2 ; c where fs and fc are the frequency shift and the 3-dB cuto frequency, respectively. Superpositions of two Gaussian psd functions with di erent values for 02; fs, and fc have been speci ed for example in [13] for describing the typical Doppler psd of large echo delays, which occur in the frequency selective mobile radio channel of the pan-European cellular GSM system. Furthermore, frequency non-shifted Gaussian psd functions, i.e. fs = 0 in (7), have been introduced for describing the Doppler psd of aeronautical channels [14, 15]. The acf of the Gaussian psd S (f ) is given by pfc r(t) = 22 e?  ln 2 t

ej2fst : (8) Observe that the acf r(t) is complex if fs 6= 0. Thus, it follows from the properties of the Fourier transform that in such cases the real Gaussian processes 1(t) and 2(t) are correlated. Fig. 2(b) shows the resulting structure of the analytical model for Rice processes (t) with underlying frequency shifted Gaussian psd characteristic. 0

B. Statistical Properties of the Analytical Model In this subsection we study besides the amplitude and phase probability density function, also the level-crossing rate, and the average duration of fades of Rice processes with given spectral shapes for the underlying Gaussian processes. We refer to these (ideal) statistical properties in Section III, where we investigate the statistics of deterministic simulation models. The probability density function (pdf) of the Rice process (t), p (x), is given by [16] 8 x ? x2+2  x  < e 2 o Io ; x0 ; p (x) = : o (9) o 0 ; x p1  ff ; > < 0 0max 1 => 1  @ f. ? fs A ; > p > : 0 fc pln 2 (

=

2 0(fmaxp)2 ; 2 2 0(fc= ln 2) ;

Jakes psd;

(20a)

Gaussian psd;

(20b)

Jakes psd; Gaussian psd; 7

(21a) (21b)

where

0 = 02.

From the above it is clear now that the higher order statistics considered here (lcr N(r) and adf T? (r)) are completely determined by the amplitude and Doppler frequency of the direct component, i.e.  and f , as well as by the shape of the acf r (t) and its time derivative of rst and second order at t = 0.

8

III. A DETERMINISTIC SIMULATION MODEL FOR RICE PROCESSES In this section, we investigate the statistical properties of deterministic simulation models for Rice processes with desired spectral characteristics for the underlying Gaussian processes. Therefore, we rst of all introduce the reader to the concept of deterministic simulation models and subsequently, we present a new and optimal method for the computation of the simulation model parameters. For reasons of comparison, we also investigate the statistical properties of the simulation system for two other parameter computation methods. Let us proceed by considering the two real functions ~1(t) and ~2(t) which will be expressed as follows

~i(t) =

Ni X n=1

ci;n cos(2fi;nt + i;n) ; i = 1; 2 ;

(22)

where Ni denotes the number of sinusoids of the function ~i(t). The quantities ci;n, fi;n, and i;n are simulation model parameters, which are adapted to the desired Doppler psd function and are therefore called Doppler coecients, discrete Doppler frequencies, and Doppler phases, respectively. It is worth mentioning that the simulation model parameters ci;n, fi;n, and i;n have to be computed during the simulation set up phase, e.g. by one of the methods described in the next subsections. Afterwards these parameters are known quantities and are kept constant during the whole simulation run phase. From the fact that all parameters are known quantities, it follows that ~i(t) itself is determined for all time and thus, ~i(t) can be considered as a deterministic function. We will see that the statistical properties of ~i (t) approximate closely the statistical properties of coloured Gaussian stochastic processes. In order to emphasize that property, we will denote henceforth the deterministic function ~i(t) as deterministic (Gaussian) process. Consequently, the interpretation of ~i(t) as a deterministic process (function) allows us to compute an analytical expression for the III.3 corresponding acf r~i i (t) via the de nition 1 r~ii (t) := Tlim !1 o 2To

ZTo

?To

~i( )~i (t +  )d :

(23)

Substituting (22) in (23) and taking thereupon the Fourier transform, we nd the following analytical expressions for the acf r~i i (t) and the corresponding psd S~ii (f )

r~ii (t) =

  S~ii (f ) = for i = 1; 2.

Ni c2 X i;n n=1

2 cos(2fi;nt) ;

(24a)

Ni c2 X i;n 4 [(f ? fi;n) + (f + fi;n )] ; n=1

9

(24b)

Notice, that the deterministic processes ~1(t) and ~2(t) are uncorrelated if and only if f1;n 6= f2;m for all n = 1; 2; : : : ; N1 and m = 1; 2; : : : ; N2. But if f1;n = f2;m is valid for some or all n; m, then ~1(t) and ~2(t) are correlated and the following cross-correlation function, r~12 (t), can be derived

8 < PN c1;n c2;m cos(2f1;nt ? 1;n  2;m) ; r~12 (t) = : n=1 2 0 ;

if f1;n = f2;m ; if f1;n 6= f2;m ;

(25)

where N indicates the minimum value of N1 and N2, i.e. N = minfN1; N2g. Fig. 3 shows the resulting general structure of the deterministic simulation model for Rice processes in its continuous-time representation. Therefrom, the discrete-time fading simulation model can immediately be obtained by replacing the time variable t by t = kT , where k is a natural number and T denotes the sampling interval. For the parameters ci;n and fi;n appearing in (22) we will present in the next two subsections a new computation procedure, which is for a given number of sinusoids Ni optimal in the following two senses: i) The pdf of the deterministic process ~i (t), named by p~i (x), is according to the L2-norm an optimal approximation of the ideal (desired) pdf pi (x) of the stochastic process i (t). ii) The acf of the deterministic process ~i(t), named by r~i i (t), is according to the L2-norm an optimal approximation of the ideal (desired) acf ri i (t) of the stochastic process i (t) within a distinct time interval.

A. Determination of Optimal Doppler Coecients In this subsection we compute the Doppler coecients ci;n in such a way that the pdf p~i (x) of the deterministic process ~i (t) results for a given number of sinusoids Ni in an optimal approximation of the ideal Gaussian pdf pi (x) of the stochastic process i (t). In the rst step, we derive the pdf of ~i(t) (see (22)). For that purpose let us consider a single sinusoid

~i;n (t) = ci;n cos(2fi;nt + i;n) :

(26)

Furthermore, let us assume in this subsection that the time t is a random variable, which ?1 ) (suppose f 6= 0). Then ~ (t) is also a is uniformly distributed in the interval (0; fi;n i;n i;n random variable and the pdf of ~i;n(t) is given by (cf. [16] on p.98)

p~i;n (x) =

(

p1 ci;n 1?(x=ci;n )2 0

; jxj < ci;n ; ; jxj  ci;n :

(27)

If the random variables ~i;n (t) are independent with respective densities p~i;n (x), then the density p~i (x) of their sum

~i(t) = ~i;1(t) + ~i;2(t) + : : : + ~i;Ni (t) 10

(28)

equals the convolution of their corresponding densities

p~i (x) = p~i;1 (x)  p~i;2 (x)  : : :  p~i ;Ni (x) :

(29)

Although (29) provides us with a rule for the computation of the density p~i (x), it is more appropriate to apply the concept of the characteristic function1. The characteristic function ~ i;n ( ) of the random variable ~i;n(t) is obtained by taking the Fourier transform of (27) ~ i;n ( ) = Jo(2ci;n ) : (30) The properties of characteristic functions are essentially the same as the properties of Fourier transforms. Thus, the convolution of the densities (29) reduces to the product of their characteristic functions, i.e. ~ i ( ) = ~ i;1 ( )  ~ i;2 ( )  : : :  ~ i;Ni ( )

(31a)

Jo(2ci;n  ) :

(31b)

=

Ni Y

n=1

Finally, from the above equation the pdf p~i (x) can be obtained by taking the inverse Fourier transform, i.e.

p~i (x) =

Z1

?1

~ i ( )e?j2x d = 2

p

Z1 Y Ni

0 n=1

Jo(2ci;n  ) cos(2x)d :

(32)

Let fi;n 6= 0 and ci;n = 0 2=Ni for all n = 1; 2; : : : ; Ni and i = 1; 2, then the sum ~i(t) (de ned by (28)) is a random variable with zero mean and variance ~02 = (c2i;1 + : : : + c2i;Ni )=2 = 02. The central limit theorem states that in the limit Ni ! 1 the random variable ~i(t) tends to a Gaussian distributed random variable i(t) with zero mean and variance 02, i.e. ? 2x22 1 e 0: lim p~ (x) = pi (x) = p (33) Ni!1 i 20 By taking the Fourier transform of the equation given above, we obtain (34) lim ~ i ( ) = i ( ) = e?2(0)2 ; Ni!1

and it follows by making use of (31b) the noteworthy result

r 2 Ni ?! 2 Ni !1 e?2(0  ) : Jo 20 N  i

 

(35)

The characteristic R 1 j2xfunction of a random variable X is de ned as the statistical average ( ) = e pX (x)dx, where  is a real variable. In spite of the positive sign in the exponential = ?1 ( ) is often called the Fourier transform of pX (x).

1 Efej2X g

11

Of course, for a limited number of sinusoids Ni , we must write p~i (x)  pi (x) and ~ p i ( )  i ( ). By considering Fig. 4(a), where p~i (x) has been evaluated for ci;n = o 2=Ni , we see that the approximation error is small and can be neglected for most practical applications if Ni  7. Now, the following question arises: Does there exist an optimal set of Doppler coecients fci;ng in such a sense that the pdf of the simulation system p~i (x) results for a given number of sinusoids Ni in an optimal approximation of the ideal Gaussian pdf pi (x)? This question can be answered by considering the following L2-norm, which is also known as the square root of the mean square error, i.e.

91=2 8 Z1 = < 2 (2) Epi := : pi (x) ? p~i (x) dx; : ?1

(36)

After substituting (32) and (33) in (36), the above error function can be minimized numerically by using a suitable optimization procedure such as the method described(optin) [18]. The minimization of (36) gives us an optimized set for the Doppler coecients ci;n . Fig. 4(b) shows the resulting pdf of the simulation system p~i (x) by using the optimized (opt) quantities ci;n = ci;n . A proper choice of the initial values of the Doppler coecients is p given by choosing ci;n = o 2=Ni . For these values the evaluation of the error function Ep(2)i is shown in Fig. 5(a). The same gure shows also the resulting error function Ep(2)i for (opt), which are for a given number of sinusoids the optimized Doppler coecients ci;n = ci;n Ni all identical after the minimization of (36). A further result worth mentioning is that the optimization procedure always decreases the initial values for the Doppler coecients, p (opt) i.e. ci;n  0 2=Ni , and therefore, the variance of the deterministic process ~i(t) is for nite Ni always smaller than the variance of the stochastic process i(t) as depicted in Fig. 5(b). But this e ect diminishes if Ni increases and can completely be neglected for reasonably values of Ni, say Ni  7. The consequence from the above is that the improvements made by the optimization procedure are negligible for most practical applications (opt) if Ni  7, and thus, we can say for simplicity that ci;n is approximately given by p ci;n = 0 2=Ni if Ni  7 (37) for all n = 1; 2; : : : ; Ni. Next, let us impose on this optimization approach the power constraint, i.e. ~02 = 02. The power constraint is nothing else but a boundary condition, which can easily be included in the optimization procedure by optimizing the rst Ni ? 1 Doppler coecients ci;1;    ; ci;Ni?1, whereas the remaining Doppler coecient ci;Ni will be computed so that ~02 = 02 yields. It follows that the result (after applying the optimization procedure) is for all Ni identical with (37) even if the initial values are choosen arbitrarily. Thus, (37) ensures for a given number of sinusoids Ni an optimal approximation of ideal Gaussian pdfs pi (x) given the power constraint. Henceforth, we will make use of (37) if we consider the approximation of stochastic processes with Gaussian distribution. In the following, we are interested in the amplitude and phase pdf of the deterministic simulation system for Rice processes as represented in Fig. 3. We only consider the timeinvariant direct component m = m1 + jm2 (see (4) with f = 0), where the pdf of the 12

real and imaginary part are given by pmi (xi) = (xi ? mi); i = 1; 2. Thus, the pdf of the deterministic process ~i (t) can be expressed by

p~i (xi) = p~i (xi)  pmi (xi) = p~i (xi ? mi) :

(38)

Let us assume that the deterministic processes ~1 (t) and ~2 (t) are statistically independent, i.e. f1;n 6= f2;m for all n = 1; 2; : : : ; N1 and m = 1; 2; : : : ; N2, then we can express the joint pdf of ~1 (t) and ~2 (t), p~1 2 (x1; x2), as follows

p~1 2 (x1; x2) = p~1 (x1)  p~2 (x2) :

(39)

Finally, the transformation of the cartesian coordinates (x1; x2) to polar coordinates (z; ) by using

x1 = z cos  ; x2 = z sin  ;

(40a,b)

allows us to derive the joint pdf of the amplitude ~(t) and phase #~(t), p~# (z; ), in the following way

p~#(z; ) = z p~1 2 (z cos ; z sin ) ;

(41a)

p~#(z; ) = z p~1 (z cos )  p~2 (z sin ) ;

(41b)

p~#(z; ) = z p~1 (z cos  ?  cos )  p~2 (z sin  ?  sin ) :

(41c)

From (41c) we can directly compute the desired amplitude pdf p~ (z) and phase pdf p~# () of the deterministic simulation system according to

Z

p~ (z) = z p~1 (z cos  ?  cos )  p~2 (z sin  ?  sin ) d ; ?

p~#() =

Z1 0

zp~1 (z cos  ?  cos )  p~2 (z sin  ?  sin ) dz :

(42a) (42b)

The result of the evaluation of (42a) and (42b) by using (32) is shown in the Figs. 6(a) and 6(b), respectively. Thereby, the Doppler coecients ci;n are computed from (37), and the number of sinusoids Ni has been restricted to N1 = N2 = 7. For the reasons of comparison, we present also in the same gures the behaviour of the ideal probability density functions p (z) and p#() according to the eqs. (9) and (10), respectively. The proposed method not only applies to the approximation of Gaussian processes or stochastic processes derived from Gaussian processes (e.g. Rayleigh processes, Rice processes, and lognormal processes) but can also be used for the approximation of other kinds of distribution functions such as the Nakagami distribution. 13

The Nakagami distribution [19], also denoted as m-distribution, o ers a greater exibility and often tends to t the pdf of the received signal amplitude computed from experimental data of certain urban channels better than the common used Rayleigh or Rice distribution [1]. The Nakagami distribution with parameters and m is de ned by m z 2m?1 e?(m= )z2 2 m p (z) = ; m  1=2 ; z  0 ; ?(m) m

(43)

where ?() is the Gamma function, =< (t)2 > denotes the time average power of the received signal strength, and m =< (t)2 >2 = < ((t)2? < (t)2 >)2 > is the inverse of the normalized variance of (t)2. We remark that for m = 1=2 and m = 1 the Nakagami distribution becomes a one-sided Gaussian and Rayleigh distribution, respectively; and nally, the Rice and lognormal distribution also can be approximated by the Nakagami distribution [19, 20]. For further details on the derivation and simulation of Nakagami fading channel models we refer to [21, 22]. In order to nd a suitable set for the Doppler coecients fci;n g, so that the amplitude pdf of the proposed simulation system p~ (z) ts the Nakagami distribution (43), we proceed in a similar way as we have successfully done for the minimization of the L2-norm (36) in connection with Gaussian distributions. We only have to replace in (36) the Gaussian distribution pi (x) by the Nakagami distribution p (z) as de ned by (43), and p~i (x) has to be replaced by p~ (z) as introduced by (42a). The optimization results for the Nakagami distribution are shown for various values of m in Fig. 7, whereby in each case N1 = N2 = 10 sinusoids have been used. Weibull distributions can be derived from uniform distributions, and on the other hand uniform distributions can be derived from Gaussian distributions. Consequently, the derivation of deterministic Weibull processes is straightforward and will not be investigated in detail.

B. Determination of Optimal Discrete Doppler Frequencies In the previous subsection, we have discussed a method which enables us to compute the Doppler coecients ci;n in such a way that the pdf p~i (x) of the deterministic process ~i(t) results in an optimal approximation of the Gaussian pdf pi (x) of the stochastic process i (t). In the present subsection, we will see that the discrete Doppler frequencies fi;n can be determined in such a way that the acf r~ii (t) of the deterministic process ~i(t) results in an optimal approximation of any desired acf rii (t) of the stochastic process i(t) within an appropriate time interval. The method discussed here is a modi cation of the Lp-norm method [10]. The proposed (modi ed) Lp-norm method assumes that the Doppler coecients ci;n are xed and given by (37), whereas optimal values for the discrete Doppler frequencies fi;n can be found numerically by minimizing the Lp-norm for p = 2, i.e. the square root of the mean square 14

error norm,

2 ZTo 31=2 Er(2)ii := 4 T1 jrii (t) ? r~ii (t)j2 dt5 ; o

(44)

0

where r~i i (t) is given by (24a), ri i (t) can be any desired acf (e.g. resulting from "snap shot" measurements of real world mobile channels), and To is de ned below and denotes an appropriate time interval [0; To] over which the approximation of rii (t) is of interest. Our investigations of the performance of the method presented above will be restricted to both Doppler power spectral densities S (f ) as introduced in Section II (i.e. Jakes and Gaussian). For the Jakes psd and the frequency non-shifted Gaussian psd, we obtain | by using (6) and (8) - for the acf rii (t) the expression

8 > < 02 Jo(2fmaxt) rii (t) = r(t)=2 = >  : 2 e?  pflnc 2 t2 0

(Jakes) ;

(45a)

(Gauss) : (45b) Numerical investigations have shown that for these acfs appropriate values for To, i.e. the upper bound of the integral in (44), are given if To is adapted to the number of sinusoids Ni according to 8 Ni < 2fmax (Jakes) ; To = : N q (46) i ln2 (Gauss) : 4fc 2 The evaluation of the error function Er(2)ii (see (44)) with the resulting optimized discrete (opt) is shown for the Jakes psd and the Gaussian psd in the Doppler frequencies fi;n = fi;n Figs. 8(a) and 8(b), respectively. For the sake of comparison, we also have shown in these gures the results obtained by applying two other ecient parameter computation techniques, which will be consisly described below.

Method of Equal Areas [5, 6]:

The principle of this method is as follows: Select a set of discrete Doppler frequencies ffi;ng in such a manner that in the range of fi;n?1  f < fi;n the area under the (symmetrical) Doppler psd Sii (f ) is equal to 02=(2Ni ) for all n = 1; 2; : : : ; Ni with fi;0 = 0. The application of the method of equal areas to the Jakes psd Sii (f ) = S (f )=2 (see (5)) gives us for the discrete Doppler frequencies  n  (47) fi;n = fmax sin 2N ; i for all n = 1; 2; : : : ; Ni and i = 1; 2. Whereas for the frequency non-shifted Gaussian psd (see (7) with fs = 0) no closedform expression for the discrete Doppler frequencies fi;n exists, and we must calculate the discrete Doppler frequencies fi;n by nding the zeros of   n ? erf fi;n pln 2 = 0 (48) Ni fc 15

for all n = 1; 2; : : : ; Ni and i = 1; 2. The Doppler coecients are exactly the same as introduced by (37). For a comprehensive derivation and discussion of the method of equal areas, we refer to [6], where also some applications and simulation results can be found.

Method of Exact Doppler Spread: This procedure is due to its simplicity and high performance predestinated to the approximation of the acf rii (t) as given by (45a), i.e. when the Jakes psd is of interest. The method can easily be derived by using the integral representation [23]

Z =2 2 Jo(z) =  cos(z sin )d ; 0 and replacing the above integral by a series expansion as follows Ni X 2 Jo(z) = Nlim cos(z sin n ) ; i !1  n=1

(49)

(50)

where n := (2n ? 1)=(4Ni ) and  := =(2Ni). Hence, we can write by using (45a) and (50)

rii (t) = o2 Jo (2fmaxt) = Nlim i !1

o2 Ni

Ni X n=1

cos(2fi;nt) ;

(51)

where

  1 fi;n = fmax sin 2N (n ? 2 ) : i

(52)

Finally, as for nite values of Ni we require rii (t)  r~ii (t) , a comparison of (51) with (24a) leads to

o2 r~ii (t) = N

Ni X

i n=1

cos(2fi;nt) :

(53)

A comparison of (53) with (24a) shows that the discrete Doppler frequencies fi;n are identical with (52), and the Doppler coecients ci;n can immediately be identi ed with those given by (37). The Doppler psd Sii (f ) is often characterized by the Doppler spread Bi , that is the square root of the second central moment (variance) of Si i (f ). In [6] it has been shown that the Doppler spread of deterministic simulation models, B~i , depends directly on the discrete Doppler frequencies fi;n and Doppler coecients ci;n. Consequently, the approximation quality Bi  B~i is mainly a ected by the employed parameter computation method. In Section IV, we will see that for the Jakes psd the proposed method always results in a Doppler spread of the simulation model B~i that is identical with the Doppler spread of the analytical model Bi , i.e. Bi = B~i . Therefore, we have designated the 16

procedure as "method of exact Doppler spread". Clearly, the method of exact Doppler spread is adapted to the Jakes psd, but a comparison of (47) and (52) motivates us - as a heuristic approach - to replace in (48) n by n ? 1=2. This enables the computation of the discrete Doppler frequencies fi;n for the frequency non-shifted Gaussian psd by nding the zeros of





2n ? 1 ? erf fi;n pln 2 = 0 2Ni fc

(54)

for all n = 1; 2; : : : ; Ni and i = 1; 2. The Doppler coecients ci;n are remain unchanged and are further on given by (37). After the calculation of the discrete Doppler frequencies fi;n and the Doppler coecients ci;n by applying one of the three methods described above, the acf of the simulation model r~ii (t) can be computed by means of (24a). The result for r~ii (t) that has been obtained by applying the L2-norm method, the method of equal areas, and the method of exact Doppler spread are presented for the Jakes psd in Fig. 9(a) and for the frequency nonshifted Gaussian psd in Fig. 9(b). In all cases, the number of sinusoids Ni was determined by Ni = 7. Observe that we have also depicted in Fig. 9(a) and Fig. 9(b) the corresponding acfs of the analytical model rii (t) = o2Jo(2fmaxt) and ri i (t) = o2 exp[?(fct)2= ln 2], respectively. This allows us directly to compare and criticize the performance of the various parameter computation methods. The L2-norm method is more complex and requires higher numerical e ort than the method of exact Doppler spread, but surprisingly, we realize from Fig. 8(a) and Fig. 9(a) that these two methods are yielding nearly equivalent and excellent approximation results in the interval 0  t  To. Thus, due to its simplicity and high approximation quality, we prefer for the Jakes psd the method of exact Doppler spread. We gather from the Figs. 8(b) and 9(b) that for the Gaussian psd the method of exact Doppler spread results in essential smaller approximation errors than the method of equal areas. A comparison of the performance of these two methods with the L2-norm method reveals undoubtedly the advantage of that optimization procedure. Therefore, we prefer for the Gaussian psd the L2-norm method. Finally, we mention that the period of the deterministic process ~i(t) (see (22)) is inversely proportional to the greatest common devisor of the discrete Doppler frequencies fi;n. For all discussed methods, the greatest common devisor can be made small enough for realistic values of Ni , say Ni  7, simply by increasing the word length of the discrete Doppler frequencies fi;n; and thus, the period of ~i(t) can be made extremely large.

C. Determination of the Doppler Phases The Doppler phases i;n are uniformly distributed random variables, which can be obtained simply by means of a random generator with uniform distribution over the interval (0; 2]. But for deterministic simulation models, where the Doppler coecients ci;n and the discrete Doppler frequencies fi;n are all determined by using deterministic methods, it is an obvious idea to introduce also for the computation of the Doppler phases i;n a deterministic approach. 17

~ i with Ni deterministic phase Therefore, we consider the following standard phase vector  components   ~i = 2 1 ; 2 2 ; : : :; 2 Ni (55) Ni + 1 Ni + 1 Ni + 1 ; and we group the Doppler phases i;n of the Ni sinusoids to a so called Doppler phase vector ~i in accordance with ~i = (i;1 ; i;2 ; : : :; i;Ni ) (56) for i = 1; 2. Now, the components of the Doppler phase vector ~i can be identi ed with ~ i. the components of the standard phase vector ~ i after permutating the components of  In this way, for a given number of sinusoids Ni one can construct Ni! di erent Doppler phase vectors ~i with equally distributed components. Let us assume that the real deterministic processes ~1(t) and ~2(t) are uncorrelated, then it can be shown that the corresponding Doppler phases 1;n and 2;n have no in uence on the statistical properties of the complex deterministic process ~(t) = ~1(t) + j ~2(t). Consequently, for one and the same sets of Doppler coecients fci;ng and discrete Doppler frequencies ffi;ng it is possible to construct Ni! di erent sets of Doppler phases fi;ng, and thus N1!  N2! di erent complex deterministic processes ~(t) = ~1(t) + j ~2(t) with di erent time behaviour but identical statistics can be realized. This is of great advantage for simulating channels for mobile digital transmission systems with frequency hopping. In such cases, a frequency hop necessitates a new permutation of the Doppler phases. On the other side, if ~1(t) and ~2(t) are correlated, then the Doppler phases i;n have an in uence on the cross-correlation function r~1 2 (t), as can easily be seen by considering (25). In this case, the Doppler phases i;n can be designed such that they control the higher order statistics of ~(t) = j~(t) + m(t)j, e.g. the simulation model's level-crossing rate and the average duration of fades [24].

18

IV. HIGHER ORDER STATISTICAL PROPERTIES OF DETERMINISTIC SIMULATION MODELS In this section, we derive for the deterministic simulation models of Figs. 3(a) and 3(b) analytical expressions for the level-crossing rate N~ (r) and the average duration of fades T~? (r). The analytical expressions do not only allow us to avoid the determination of N~ (r) and T~? (r) from time consuming simulations of the channel output signal ~(t), but also enable us to study the degradation e ects due to the selected parameter computation method and the limited number of sinusoids Ni. Moreover, we investigate the pdf of fading time intervals by using the Jakes and Gaussian psds. In Chapter A of Section III it has been shown that the pdf of the deterministic process ~i(t), p~i (x), is nearly identical with the pdf of the stochastic process i(t), pi (x), if the number of sinusoids Ni is suciently large, say Ni  7. On the assumption that p~i (x) = pi (x) yields, the level-crossing rate of the deterministic simulation model N~(r) is furthermore de ned by (12); one merely has to replace the quantities of the analytical model and (consider therefore (13) and (14)) by the corresponding quantities of the simulation model ~ and ~. Thus ! ~_ 0 q ~  ~ = 2f ? ~ 2 ; (57) 0 _2 ~   ~ = ? ~0 ? ~0 ; 0

(58)

where ~0, ~_ 0 are related by the acf r~11 (t) (see (24a)) and ccf r~12 (t) (see (25)) of the simulation model as follows ~ = d2 r~  (t) = r~  (0) ; (59) 1 1 0 dt2 1 1 t=0 ~_ 0 = d r~1 2 (t) = r~_ 1 2 (0) ; (60) dt t=0 P i c2 =2. Let the Doppler coecients c are as given by (37), and ~0 = r~11 (0) = Nn=1 i;n 1;n then ~0 = 0 yields. The quality of the approximations  ~ and  ~ depends on the underlying Doppler psd and will be investigated separatly for the Jakes and frequency shifted Gaussian psds.

A. Higher Order Statistical Properties by Using the Jakes PSD The Jakes psd (5) is a symmetrical function and thus, it follows from the properties of the Fourier transform that the ccf r12 (t) is zero, i.e. 1 (t) and 2(t) are uncorrelated. 19

Now, we impose the same requirement on the simulation system. From (25) it follows that r~1 2 (t) = 0 can easily be achieved for all design methods discussed here by selecting the number of sinusoids such that N1 = N2 ? 1 yields. Then N1 is unequal to N2 and per de nition the discrete Doppler frequencies are designed in such a way that f1;n 6= f2;m holds for all n = 1; 2; : : : ; N1 and m = 1; 2; : : : ; N2. Thus, it follows from (57) - (60) in connection with (24a) that

q

~ = 2f= 2 ~ and

~ = ? ~0 = 22

(61) N1 X n=1

(c1;nf1;n)2 :

(62)

The evaluation of the above equation by using the method of exact Doppler spread can be performed by substituting (37) and (52) in (62), i.e.

~ = (20fmax)2  N1

N1 X

1 n=1

By taking notice of



0 = 02

sin2

  1 2N1 (n ? 2 ) :

(63)

and making use of the relation



N1 1 X 2  (n ? 1 ) = 1 ; sin N1 n=1 2N1 2 2

(64)

we nally obtain for ~ and ~ the result ~ = p1  ff ; 0 max ~ = 2 0(fmax)2 :

(65) (66)

A comparison of (20a) and (21a) with (65) and (66) reveals that ~ = and ~ = , i.e., the approximation q error is zero, and thepDoppler spreads B~i and Bi are identical too, ~ o=(2) and Bi = = o=(2). because B~i = = In a similar way, we can compute ~ by using the method of equal areas. In this case, we obtain ~ = (1 + 1=Ni ) and ~ tends to when Ni ! 1. Let us for a moment assume that the Doppler frequency of the direct component f is equal to zero, i.e. f = 0, then = ~ = 0 and the lcr N (r) reduces to (17). Now the relative error of N~ (r) can be expressed by

q

p ~ (r) ? N(r) ~ ? N "N~ = : = p N (r)

(67)

20

q

The substitution of ~ = (1 + 1=Ni) in (67) and the use of the approximation 1 + N1i  1 + 2N1 i allows us to express the relative error "N~ directly as function of the number of sinusoids 1 : "N~  2N (68) i

The preceding equation shows that by using the method of equal areas about 25 sinusoids are required for the design of ~1(t) and ~2(t) in order to reduce the relative error of N~(r) to 2 %.

It shoud be noted that for the L2-norm method no closed-form expression for ~ can be derived. The use of that method leads after the numerical evaluation of (62) to the results 2 as depicted in Fig. 10(a). This gure also shows the ~ max for the normalized quantity =f 2 as function of the number of sinusoids Ni by using the method of ~ max behaviour of =f equal areas and the method of exact Doppler spread. Next, we consider the deterministic simulation model of Fig. 3(a), and we calculate directly the lcr N~ (r) by counting the number of crossings per second at which the simulated channel output signal ~(t) crosses a speci ed signal level r with positive slope. Therefore, we have chosen a maximum Doppler frequency fmax (see (5)) of fmax = 91:73 Hz, that corresponds to a mobile speed of 110 km/h and a carrier frequency of 900 MHz. The variance 02 was normalized to 02 = 1=2. The parameters of the deterministic simulation system have been designed according to the method of exact Doppler spread with N1 = 7 and N2 = 8 sinusoids, so that the cross-correlation between the real part ~1(t) and the imaginary part ~2(t) is zero. For the realization of the discrete-time deterministic simulation system we have used a sampling interval T of T = 3  10?4 s and altogether K = 4  105 samples of the deterministic process ~(kT ), k = 0; 1; : : : ; K ? 1, have been simulated for the evaluation of the lcr N~ (r). The result for the normalized lcr of the deterministic simulation model N~ (r)=fmax is presented in Fig. 11(a). This gure demonstrates for a deterministic Rice process with  = 1 the in uence of the Doppler frequency of the direct component f on the behaviour of N~ (r). For the reasons of comparison, we also have shown in Fig. 11(a) the lcr of the (ideal) analytical model N (r) as de ned by (12). It can be seen that the lcr of the simulation model N~ (r) coincides extremely closely with the lcr of the theoretical model N (r) by applying the method of exact Doppler spread.

B. Higher Order Statistical Properties by Using the Gaussian PSD The Gaussian psd (7) is for fs 6= 0 an unsymmetrical function. Consequently, the ccf r1 2 (t) is unequal to zero and thus 1 (t) and 2(t) are correlated. Let the stochastic processes 1(t) and 2(t) (see Fig. 2(b)) are uncorrelated, then we yield after some simple computations the following expressions for the acf r11 (t) and ccf r12 (t)

r11 (t) = 21 [r11 (t) + r22 (t)]cos(2fst) ; r12 (t) = 21 [r11 (t) + r22 (t)]sin(2fs t) ; 21

(69a) (69b)

respectively. Starting with the above equations and using ri i (0) = 02 = 0, we can show that (15) and (16) are resulting for the frequency shifted Gaussian psd in 0 = ?

_0 =

o

h

i

p

2(fc= ln 2)2 + (2fs)2 ;

0  2fs

(70a)

;

(70b)

respectively. In a similar way, the corresponding quantities of the deterministic simulation model can be derived. The result is ~0 = ?

~_0 =

"X N1 2 n=1

0  2fs

(c1;nf1;n)2 +

N2 X n=1

#

(c2;nf2;n)2 ? 0(2fs)2 ;

;

(71a) (71b)

where we have used r~ii (0) = ~0 = 0, i.e., the Doppler coecients are given by ci;n = p 0 2=Ni . A comparison of (70b) with (71b) shows that ~_0 = _0. Now, we are able to compute the quantities ~ and ~, which determine the approximation quality of the deterministic simulation model's level-crossing rate. Therefore, we substitute (71a), (71b) in (57), (58) and obtain nally

s

~ =  2~ (f ? fs ) ; ~ = 2

"X N1 n=1

(c1;nf1;n)2 +

(72) N2 X n=1

#

(c2;nf2;n)2 :

(73)

Eq. (72) allows further observations: when ~ ! , ~ tends to and consequently N~(r) tends to N (r). Obviously, the quantity ~ is responsible for the approximation quality of the higher order statistical properties of the deterministic simulation model. For all parameter computation methods discussed in Section III, ~ can not be expressed in a closed-form for the Gaussian psd, so that (73) has to be evaluated numerically. Fig. 10(b) shows the numerical evaluation results for the normalized quantity ~=fc2, which have been obtained by employing all parameter computation methods investigated in the present paper. Now, let us consider Fig. 3(b) and let us determine for the frequency shifted Gaussian psd the lcr N~ (r) from the simulated sequence ~(kT ), k = 0; 1; : : : ; K ? 1, in the same way as we have already done for the Jakes psd. Therefore, p we have selected the following model parameters: T = 3  10?4 s, K = 4  105 , fc = ln 2  fmax, fmax = 91:73 Hz,  = 1, f = 0, andin order to study the in uence of the frequency shift fs , we have used the values fs 2 0; 32 fmax; fmax . p Observe, by considering of (21), that the 3-dB cuto frequency fc = ln 2  fmax was selected such that the quantity is identical for both power spectral densities (Jakes 22

and Gaussian). The discrete Doppler frequencies fi;n and Doppler coecients ci;n have been determined by the method of exact Doppler spread with N1 = 25 and N2 = 26. The results for the normalized level crossing rates of the deterministic simulation model N~ (r)=fmax as well as of the analytical model N (r)=fmax are shown in Fig. 11(b). A comparison of Fig. 11(a) with Fig. 11(b) shows the same result for both the analytical and the simulation model: identical level-crossing rates are observed | as it was expected by the theory | although the underlying Doppler power spectral densities are quite di erent. In a similar way, an expression for the average duration of fades of the simulation model, T~? (r), can be derived. One merely has to replace in (18) N (r) by N~ (r) and P? (r) by P~? (r). Thus, the analytical and numerical investigations of T~? (r) are straightforward and will not be presented here.

C. Distribution of the Fading Time Intervals The performance of mobile communication systems can mainly be improved by an ecient design of the interleaver/deinterleaver and the channel coding/decoding unit. For the speci cation of these units, not only the knowledge of the level-crossing rate and the average duration of fades is of special interest, but also the distribution of the fading time intervals. The rest of this section is focused on the probability density function of the fading time intervals that will be denoted here by p0? (? ; r). This pdf p0?(? ; r) is the conditional probability density that the stochastic process (t) crosses a certain level r for the rst time at t2 = t1 + ? with positive slope given a crossing with negative slope at time t1. An exact analytical evaluation of p0? (? ; r) is still an unsolved problem even for Rayleigh and Rice processes. So we will determine p~0? (?; r) experimentally by employing deterministic simulation systems. For that purpose, we simulate ~(kT ) with a sampling interval of T = 0:5  10?4 s, and we x  = 0, i.e. ~(kT ) is a deterministic Rayleigh process. Typical measurement results of p~0?(? ; r) for r = 0:1 and r = 1 by evaluating 2  104 fading time intervals ? are presented in Figs. 12(a) and 12(b) for the Jakes psd and in Figs. 12(c) and 12(d) for the Gaussian psd, respectively, whereby we have applied the method of exact Doppler spread by using various values for the pnumber of sinusoids Ni. The remaining model parameters were fmax = 91 Hz, fc = fmax ln 2; o2 = 1=2, and fs = 0. Approximative solutions for p0?(? ; r) can be found for example in [25, 26, 27, 28, 29]. The rst approximation to p0? (? ; r) is p1? (?; r), where p1? (? ; r) is the conditional probability density that the stochastic process (t) passes upwards through a certain level r at time t2 = t1 + ? given that (t) has passed downward through r at time t1. Note, that nothing is said about the behaviour of (t) between t = t1 and t = t2. In [25] an analytical solution for p1?(? ; r) can be found that tends for deep fades to





d 2 I ( 2 )e? u22 ; if r ! 0 ; p1? (?; r) ! ? du uT? 1 u2

(74)

where u = ? =T? (r) and I1() denotes the modi ed Bessel function of rst order. The graph of this function also is shown for the Rayleigh process (t) with underlying Jakes and Gaussian psd characteristics in Figs. 12(a) and 12(c), respectively. These gures impressivly show the excellent accordance between the theoretical approximation p1? (? ; r) 23

and the measurement results p~0? (? ; r) for deep fades, i.e. for small values of r, where the probability that only short fading time intervals ? occur is high. Moreover, Fig. 12 led us to the hypothesis that stochastic processes (t) with di erent Doppler psd characteristics but identical values for the quantities and have identical probability density functions of the fading time intervals p0? (? ; r) for small values of r. Assuming that the hypothesis is true, then some consequences arise for the modelling of real-world mobile channels on the basis of measurements of the statistics of the received signal. We conclude from the above that for the simulation model's Doppler psd it is not necessary to match exactly the measured Doppler psd, but the approximation quality of the respective autocorrelation functions around the time origin is of utmost importance for the modelling of the measured channel statistics.

24

V. CONCLUSION In this paper we have investigated the statistical properties of deterministic simulation models for mobile fading channels. Especially for frequency non-selective deterministic simulation models, we have not only derived analytical expressions for the probability density function of the amplitude and phase, but also higher order statistics like levelcrossing rate and average duration of fades have been investigated analytically. This makes the time consuming evaluation and estimation of the statistics from simulated channel output sequences super uous. Moreover, the derived analytical expressions allow detailed investigations of the degradation e ects due to a limited number of sinusoids or due to the applied computation method for the simulation model's discrete Doppler frequencies and Doppler coecients, and they provide us with a powerful tool when discussing and comparing the eciency of di erent parameter computation methods. A performance study of three ecient parameter computation methods was also a topic of the present paper, whereby two of the considered methods are new. In particular for the often used Jakes and Gaussian Doppler power spectral densities, it turned out by a comparison of the so called method of exact Doppler spread with the method of equal areas that the former allows a drastic reduction in the number of sinusoids, without producing noteworthy deviations from the desired (ideal) statistics. The third introduced method, which we have named L2-norm method, requires a greater numerical expenditure than the other two methods. This method develops its full performance by the approximation of measured autocorrelation functions obtained from real-world mobile radio channels, where other more conventional parameter computation methods tend to fail.

25

References [1] H. Suzuki, \A statistical model for urban radio propagation," IEEE Trans. Commun., vol. COM-25, no. 7, pp. 673{680, July 1977. [2] J.D. Parsons, The Mobile Radio Propagation Channel. London: Pentech Press, 1992, p. 166. [3] S.O. Rice, \Mathematical analysis of random noise," Bell Syst. Tech. Journal, vol. 23, pp. 282{332, July 1944. [4] S.O. Rice, \Mathematical analysis of random noise," Bell Syst. Tech. Journal, vol. 24, pp. 46{156, Jan. 1945. [5] M. Patzold, U. Killat, F. Laue, \A deterministic model for a shadowed Rayleigh land mobile radio channel," in Proc. of the Fifth IEEE International Symposium on Personal, Indoor, and Mobile Radio Communications (PIMRC'94), pp. 1202{1210, The Hague, Netherlands, Sept. 1994. [6] M. Patzold, U. Killat, F. Laue, \A deterministic digital simulation model for Suzuki processes with application to a shadowed Rayleigh land mobile radio channel," IEEE Trans. Veh. Technol., vol. VT-45, no.2, pp. 318-331, May 1996. [7] H. Schulze, \Stochastic models and digital simulation of mobile channels (in German)," in U.R.S.I/ITG Conf. in Kleinheubach 1988, Germany (FR), vol. 32, pp. 473{483, Proc. Kleinheubacher Berichte by the German PTT, Darmstadt, Germany, 1989. [8] P. Hoher, \A statistical discrete-time model for the WSSUS multipath channel," IEEE Trans. Veh. Technol., vol. VT-41, no. 4, pp. 461{468, Nov. 1992. [9] M. Patzold, U. Killat, F. Laue, Y. Li, \On the problems of Monte Carlo method based simulation models for mobile radio channels," in Proc. IEEE 4th Int. Symp. on Spread Spectrum Techniques & Applications (ISSSTA'96), Mainz, Germany, pp. 1214-1220, Sept. 1996. [10] M. Patzold, U. Killat, Y. Shi, F. Laue, \A deterministic method for the derivation of a discrete WSSUS multipath fading channel model," European Transactions on Telecommunications and Related Technologies (ETT), vol. ETT-7, no. 2, pp. 165{175, March{April 1996. [11] W.C. Jakes, Ed., Microwave Mobile Communications. New York: IEEE Press, 1993. [12] E.F. Casas, C. Leung, \A simple digital fading simulator for mobile Radio," in Proceedings of the IEEE Vehicular Technology Conference, pp. 212{217, Sept. 1988. [13] CEPT/COST 207 WG1, \Proposal on channel transfer functions to be used in GSM tests late 1986," COST 207 TD (86)51 Rev. 3, Sept. 1986. [14] A. Neul, J. Hagenauer, W. Papke, F. Dolainsky, F. Edbauer, \Aeronautical channel characterization based on measurement ights," IEEE Conf. GLOBECOM'87, pp. 1654{1659, Tokyo, Japan, Nov. 1987. 26

[15] P.A. Bello, \Aeronautical channel characterization," IEEE Trans. Commun., vol. COM-21, pp. 548{563, May 1973. [16] A. Papoulis, Probabilities, Random Variables, and Stochastic Processes. New York: McGraw-Hill, 3 edition, 1991. [17] M. Patzold, U. Killat, F. Laue, \An extended Suzuki model for land mobile satellite channels and its statistical properties," IEEE Trans. Veh. Technol., submitted and accepted for publication, 1995. [18] R. Fletcher, M.J.D. Powell, \A rapidly convergent descent method for minimization," Computer Journal, vol. CJ-6, no. 2, pp. 163{168, 1963. [19] M. Nakagami, \The m-distribution { A general formula of intensity distribution of rapid fading," in Statistical Methods in Radio Wave Propagation, W.G. Ho man (Ed.), Oxford, England: Pergamon Press, 1960. [20] U. Charash, \Reception through Nakagami fading multipath channels with random delays," IEEE Trans. Commun., vol. COM-27, no. 4, pp. 657{670, Apr. 1979. [21] W.R. Braun, U. Dersch, \A physical mobile radio channel model," IEEE Trans. Veh. Technol., vol. VT-40, no. 2, pp. 472{482, May 1991. [22] U. Dersch, R.J. Ruegg, \Simulations of the time and frequency selective outdoor mobile radio channel," IEEE Trans. Veh. Technol., vol. VT-42, no. 3, pp. 338{344, Aug. 1993. [23] M. Abramowitz, I.A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Washington: National Bureau of Standards, 1972. [24] M. Patzold, U. Killat, F. Laue, Y. Li, \An ecient deterministic simulation model for land mobile satellite channels," in Proc. IEEE Veh. Tech. Conf., Atlanta, Georgia, USA, pp. 1028{1032, Apr./May 1996. [25] S.O. Rice, \Distribution of the duration of fades in radio transmission: Gaussian noise model," Bell Syst. Tech. Journal, vol. 37, pp. 581{635, May 1958. [26] J.A. McFadden, \The axis-crossing intervals of random functions-II," IRE Trans. Inform. Theory, vol. IT-4, pp. 14{24, Mar. 1958. [27] M.S. Longuet-Higgins, \The distribution of intervals between zeros of a stationary random function," Phil. Trans. Royal. Soc., vol. A 254, pp. 557{599, 1962. [28] A.J. Rainal, \Axis-crossing intervals of Rayleigh processes," Bell Syst. Tech. Journal, vol. 44, pp. 1219{1224, 1965. [29] D. Wolf, T. Munakata, J. Wehhofer, \Statistical properties of Rice fading processes," in Signal Processing II: Theories and Applications, Proc. Eusipco-83 Second European Signal Processing Conference, H.W. Schuler (Ed.), Erlangen, Germany, pp. 17{20, Elsevier Science Publishers B.V. (North-Holland), 1983.

27

LIST OF FIGURE CAPTIONS Fig. 1: Models for the realization of coloured Gaussian noise processes

(a) by ltering a white Gaussian noice process and (b) by a nite sum of weighted sinusoids. Fig. 2: Stochastic analytical models for Rice processes with given Doppler psd shape according to p (a) the Jakes psd, where Hi(f ) = 0=[fmax 1 ? (f=fmax )2]1=2, and  p ? 1=4 exp ? ln 2f=f . (b) the frequency shifted Gaussian psd, where Hi(f ) = p0fc ln2 c 

Fig. 3: Deterministic simulation models for Rice processes with given Doppler psd

shape according to (a) the Jakes psd and (b) the frequency shifted Gaussian psd. Fig. 4: The approximated Gaussian pdf p~i (x) for p (a) ci;n = 0 2=Ni and for Ni = 5; 7; 1 and 02 = 1. (b) ci;n = c(opt) i;n

Fig. 5: (a) Error function Ep(2)i and

p

(b) variance ~02 of the deterministic process ~i(t) for ci;n = 0 2=Ni (ooo) and ci;n = c(opt) (***) with 02 = 1. i;n Fig. 6: (a) Amplitude pdfs p (z) (analytical model), p~ (z) (simulaion model), and (b) phase pdfs p#() (analytical model), p~# () (simulation model) for N1 = N2 = 7 (f = 0;  = ; 02 = 1=2). Fig. 7: Approximation of the Nakagami distribution with N1 = N2 = 10 ( = 1).

Fig. 8: Error function Er(2)ii in connection with (a) the Jakes psd and (b) the Gaussian

psd (02 = 1, fmax = fc = 91 Hz; fs = 0). Fig. 9: The ideal autocorrelation functions (a) ri i (t) = 02J0(2fmaxt) (Jakes psd) and (b) rii (t) = 02 exp [?(fct)2= ln 2] (Gaussian psd) in comparison with the respective approximated autocorrelation functions r~i i (t) for Ni = 7 (02 = 1, fmax = fc = 91 Hz, fs = 0). Fig. 10: Characteristic quantities and ~ of Sii (f ) and S~ii (f ), respectively, for: (a) the Jakes psd and (b) the Gaussian psd. Fig. 11: The normalized level-crossing rate of the analytical model N (r)=fmax for 02 = 1=2,  = 1 and of the simulation model for (a) the Jakes psd p (N1 = 7, N2 = 8) and (b) the Gaussian psd (N1 = 25, N2 = 26) with fc = fmax ln 2 and f = 0. Fig. 12: The pdf of the fading time intervals p~0? ( ; r) of the deterministic Rayleigh process ~(t) by using the Jakes psd (02 = 1=2; fmax = 91 Hz;  = 0) if the signal level ~(t) = r is (a) r = 0p :1 and (b) r = 1, as well as by using the Gaussian psd (02 = 1=2; fc = fmax ln 2;  = 0) if (c) r = 0:1 and (d) r = 1. 28

(a)

(b)

Figure 1

29

(a)

(b)

Figure 2

30

(a)

(b)

Figure 3

31

(a) p˜ µ (x) i

0.4 0.35 0.3

c i,n = σ o (2/N i ) 1/2

Ni = ∞ Ni = 7 Ni = 5

0.25 0.2 0.15 0.1 0.05 0 -3

-2

-1

0

1

2

3



x

(b) p˜ µ (x) i

0.4 0.35 0.3

Ni = ∞

c i,n = c (opt) i,n

Ni = 7 Ni = 5

0.25 0.2 0.15 0.1 0.05 0 -3

-2

-1

0

1 →

x

Figure 4

32

2

3

(a) E (2) p

µi

0.02 o

E (2) if c i,n = σ o (2/N i ) 1/2 p

*

E (2) if c i,n = p

µi

0.015 c (opt) i,n

µi

0.01

0.005

0 0

5

10

15

20

25

30

25

30



Ni

(b) ˜ 2o σ 1.1 1 0.9 0.8 0.7

o

˜ 2o if c i,n = σ o (2/N i ) 1/2 σ

*

˜ 2o if c i,n = σ

c (opt) i,n

0.6 0.5 0

5

10

15

20 →

Ni

Figure 5

33

(a) p ξ (z), p˜ ξ (z)

0.9 0.8

ρ= 0

0.7

ρ= 1/2

0.6

ρ= 1

0.5 0.4 0.3 0.2

analy. model

0.1

simul. model

0 0

0.5

1

1.5 →

z

(b)

2

2.5

p ϑ (θ), p˜ ϑ (θ)

0.5

analy. model simul. model

ρ= 1

0.4

ρ=1/2

0.3

ρ= 0

0.2 0.1 0 0

π θ

2π →

Figure 6 34

Nakagami distribution 1.6 1.4

p ξ (z) (analy. model)

1.2

p˜ ξ (z) (simul. model)

1

m=3

0.8

m=2

0.6

m=1

0.4 0.2 0 0

0.5

1

1.5

2 →

z

Figure 7

35

2.5

3

(a) E (2) r

µiµi

0.14

Method of Exact Doppler Spread

0.12

Method of Equal Areas 0.1

L p -Norm Method

0.08 0.06 0.04 0.02 0 0

5

10

15

25

30



Ni

(b)

20

E (2) r

µiµi

0.5

Method of Exact Doppler Spread

0.4

Method of Equal Areas 0.3

L p -Norm Method

0.2

0.1

0 0

5

10

15

20 →

Ni

Figure 8 36

25

30

(a) rµ

iµi

(t), r˜µ

iµi

(t)

Method of Exact Doppler Spread 1 Method of Equal Areas L 2 -Norm Method Analytical Model

0.5

0

-0.5 0

0.025 To t/sec

(b)



iµi

(t), r˜µ

0.05 →

iµi

0.075

(t)

Method of Exact Doppler Spread 1 Method of Equal Areas 0.8 L 2 -Norm Method 0.6 Analytical Model 0.4 0.2 0 -0.2 0

0.002

0.004

0.006 t/sec

Figure 9 37

0.008 →

0.010

To

(a)

β/f 2max , β˜/f 2max

Method of Exact Doppler Spread Method of Equal Areas L 2 -Norm Method

30 -

Analytical Model

20

0

5

10

15

25

30



Ni

(b)

20

β/f 2c, β˜/f 2c

450 400

Method of Exact Doppler Spread

350

Method of Equal Areas

300

L 2 -Norm Method -

250

Analytical Model

200 150 100 50 0 0

5

10

15

20 →

Ni

Figure 10 38

25

30

(a) N ξ (r)/f max 1.2

Analytical Model Simulation (N 1 =7, N 2 =8)

1 0.8

f ρ /f max = 1

0.6

f ρ /f max = 2/3

0.4

f ρ /f max = 0

0.2 0 0

1

2

3

4

r

(b) N ξ (r)/f max 1.2

Analytical Model Simulation (N 1 =25, N 2 =26)

1 0.8

f s /f max = 1

0.6

f s /f max = 2/3

0.4

f s /f max = 0

0.2 0 0

1

2

3 r

Figure 11 39

4

(a)

(b) p˜ 0 (τ - ,r) ⋅ T ξ -

-

N 1 = 4, N 2 = 5 N 1 = 6, N 2 = 7 N 1 =24, N 2 =25 analytical approx. ( p 1 (τ - ,r) )

1 0.8

p˜ 0 (τ - ,r) ⋅ T ξ -

1.8

-

N 1 = 4, N 2 = 5 N 1 = 6, N 2 = 7 N 1 =24, N 2 =25

1.6 1.4 1.2

ρ = 0 fρ = 0 r = 1

1

-

0.6 0.8

ρ = 0 fρ = 0 r = 0.1

0.4

0.6 0.4

0.2 0.2 0 0

1

2 τ - /T ξ

3

0 0

4

1

2 τ - /T ξ

-

(c)

3

4

-

(d) p˜ 0 (τ - ,r) ⋅ T ξ -

-

N 1 =5, N 2 =6 N 1 =25, N 2 =26 N 1 =50, N 2 =51 analytical approx. ( p 1 (τ - ,r) )

1 0.8

p˜ 0 (τ - ,r) ⋅ T ξ -

1

N 1 =5, N 2 =6 N 1 =25, N 2 =26

0.8

N 1 =50, N 2 =51 0.6

-

ρ = 0

0.6 ρ = 0 fs = 0 r = 0.1

0.4

fs = 0

0.4

r = 1 0.2

0.2 0 0

-

1

2 τ - /T ξ

3

0 0

4

-

Figure 12

40

1

2 τ - /T ξ

3 -

4