Link. r exp (jo)=z exp W O ) + w exp (j4), z, w>o (1) From this, p(r) is given by. p(r) = r/bo 56 exp [- (r2 +z2)/2bo]zo(rz/b~)p(z) dz.

Transcription

1 22 IEEE TRANSACTIOKS ON VEHICULAR TECHNOLOGY, VOL. VT-3, NO. 3, AUGUST 985 A Statistical Model for a Land Mobile Satellite Link Abstract-A statistical model for a land mobile satellite link is LI. STATISTICAL MODEL described. The model assumes that the amplitude of the line-of-sight The mathematical derivations required to describe the component under foliage attenuation (shadowing) is lognormally distributed and the received multipatb interference has a Rayleigh distribution. model are given in this section. Expressions for the level crossing rate and average fade duration are given. Comparison between results calculated from the model and from A. Probability Density Function measurements show reasonably good agreement. The model should be useful for designing communications systems and for simulating propaga- Consider the sum of a lognormally distributed random tion effects in the laboratory. phasor and a Rayleigh phasor [lo], I. INTRODUCTION STATISTICAL model for received signal envelope A fading encountered on a land mobile satellite channel is useful for predicting communications systems performance with various modulation schemes. The model will also facilitate the design of propagation simulator. There have been many papers published on the modeling of mobile radio channels in urban environments, a few of which are given in []-[9]. The model which is of interest here is for a land mobile satellite channel in a rural environment, where, for most of the time, a line-of-sight signal component is available at the receiver. This paper describes a statistical model in terms of its probability density function, level crossing rate (LCR) and average fade duration (AFD). These characteristics were obtained under the hypothesis that foliage not only attenuates but also scatters radio waves at UHF and at L-band frequencies. The model assumes thathe line-of-sight component under foliage attenuation (shadowing) is lognormally distributed and that the multipath effect is Rayleigh distributed. Additionally, these two random processes are assumed to be correlated. Although no attempt is made at physical justification of the model, results derived from the model are compared with those derived from measurements and are found to be good engineering approximations. The probability density function of the model facilitates the calculation of fade margins in the designs of communication systems. The results of LCR and AFD provide information relevant to modem designs and error correction techniques. Section I of the paper describes the statistical model and numerical results are given in Section. Conclusions are presented in Section IV. Manuscript received February 2, 98, revised January 7, 985. The author is with the Department of Communications, Communications Research Centre, P.O. Box 9. Station H, Ottawa, Ontario, Canada K2H 8S2. Telephone (63) r exp (jo)= exp O ) + w exp (j),, w>o () where the phases o and are uniformly distributed between and 2a, is lognormally distributed, and w has a Rayleigh distribution. If is temporarily kept constant, then the conditional probability density function of r is simply that of a Rician vector [ : where bo represents the average scattered power due to multipath, and Zo ( a ) is the modified Bessel function of eroth order. Applying the theorem of total probability, one obtains From this, p(r) is given by p(r) = r/bo 56 exp [- (r2 +2)/2bo]Zo(r/b~)p() d. () It has been assumed that p() is lognormal, given by where &and p are the standar deviation and mean, respectively. Thus, substituting (5) into (), one obtains p(r)=r/(bo=) jm l/ exp [-(ln It has been shown [lo] that r is lognormal for large values OO8-955/85/8-22$.OO 985 IEEE

2 LOO: LAND MOBILE SATELLITE LINK 23 and Rayleigh distributed for small values as given by (7). where B = bob2 - bl, is the amplitude of the sine wave and the b are the moments [ 3. exp [- r- pc>2/2d for % 6 (7) For mathematical simplicity, assume that the noise (multipath) spectrum isymmetrical abouthe line of sight compo- r/bo exp [ - r2/2bo] for r a &. nent frequency. Then bl =O, B= bob2, For intermediate values of r, the complicated integral given by (6) must be evaluated numerically. The following relation and from 39 the joint Probability becomes can be used to allow integration over [, R] instead of [R, 3. p(r> R) = lm p(r) dr R R = - p(r) dr. lr (8) r(2a) p(r, f)= bg 2 exp [ - (r2 + 2)/2bo] exp [ - i2/2b2] exp [r cos B/bo] db. () --* The integral in () can be recognied as the modified Bessel B. Level Crossing Rate and Average Fade Duration function of eroth order and is given by 2a (r/bo). Thus the In this subsection, expressions for the LCR and AFD of the expression given received signal due to the combined effect of shadowing and multidath on a land mobile satellite link are derived. RiEe [I2 hashown that for any wide sense stationary and random process, NR (which is the LCR), defined as the expected rate at which the envelope crosses a specified signal level R with positive slope is given by the following integral: i) di, where the overdot indicates time derivative and p(r, i) is the joint probability density function of i and r at r = R. If the lognormal variation of the direct sign,al component is ignored, the remaining signal envelope would have a Rician distribution, and it is easy to show under simplifying assumptions, based on previous work [ 3, that the resultant received signal envelope and its rate of change are statistically independent. Similarly, if the multipath component is ignored, the same statistical independence can be shown [9]. This gives justification for the reduction of (9) into a more tractable form. In addition, based on the referenced analyses, it can be shown that the rate of change of the resultant envelope is a Gaussian random process under both fading conditions. This can be used for the determination of LCR s and the AFD of the combined received signal. First, to show the desired properties of the received signal under Rician fading conditions, let the lognormal component () remain temporarily constant as in Section -A. Then the received signal can be modelled as a sine wave plus narrowband Gaussian noise, and from Rice [3], in ( ) can be separated into two parts: p(i) = - exp [ - i2/2b2]. ( 3) 5& From (2) and (3), it is clear that the rate of change of the (9) signal envelope i is uncorrelated with the envelope, r; that is, p(r, i) = p(r) p(i). Further, (3) shows that the probability density function of i is Gaussian. If one now considers the lognormal component in isolation, it can be shown [9] that the rate of change of r, i, is distributed normally. The derivation is rewritten here for convenience. Following [9], the joint probability density p(r, i) is given by p(r, i)= 2-r exp [-(In r - ~ )~/2d~ - i2/2d2] () where p and do are the mean. and variance of the envelope r, and d2 is the variance of i. One can write () as where and p(r) = - exp [- (In r - ~ )~/2d~] (5) s r p(i) = - exp [ - i2/2d2]. This shows that r and i are also uncorrelated for lognormal fading. Equations (5) and (6) are lognormal and normal it exp - [r2/2bo- r cos B/bo+2/2b --* respectively. distribution, It has been shown for both Rician fading and lognormal variation of the line-of-sight signal, that envelope r is () uncorrelated with the rate of change of the envelope i. The + boi2/2b + 2bl sin B/B + b:2 sin2 B/2Bbo] db,

3 2 IEEE TRANSACTIONS VEHICULAR TECHNOLOGY, VOL. VT-3, NO. 3, AUGUST 985 expression in (9) for LCR can thus be rewritten as NR=p(r) ip(i) di, ( 7) Then substituting (22) into (7), the expression of LCR is given by where p(r) is the probability density function of the received envelope due to the combined effects of multipath and shadowing, given by (6). In (7), one needs to determine the probability density function of the rate of change of the envelope, p(i). Then to obtain values of the LCR, one must evaluate the integral. It has been shown in previous paragraphs that the probability density function of the rate of change of envelope, i, for both multipath and shadowing is Gaussian. Therefore, it seems reasonable to assume that the probability density function of the rate of change of received signal envelope, p(i), comes from two correlated Gaussian random processes. This allows one to obtain results with varying degrees of correlation, for comparison with results derived from measurements. Let i combined rate of change of the envelope; x rate of change of the envelope due to multipath with standard deviation &and mean = ; j rate of change of envelope due to shadowing with standard deviation &%-and mean = ; ith the assumption that the fading spectrum due to shadowing and multipath is symmetrical and Gaussian [ 3, that is b2 = bo(2tfm) d2 = do(2tfm) (2) where fm is the maximum Doppler-shift frequency, (23) can be normalied with respect to fm to make results independent of vehicle velocity, giving (25) The average fade duration, TV is related to the LCR [8] by the equation then i = x + y. The correlated Gaussian random process of two variables is given by [] as exp [ -- (X2/b2-2pkj/JbZd +j2/d2 (8) 2( - P2) where p is the correlation coefficient. After substituting 2 = i - 3 into (8) and some algebraic manipulations, the following integral can be evaluated. Crn P(i) = [2~(l - p2)(b2 + 2 p G where a2 in (2) is given by Using (2), the integral given in + d2)] (7) can be evaluated, as m ip(i) di= 6 Jb2+2pJbdfd CY This completes the mathematical derivation required to obtain expressions for the probability density function, LCR, and AFD due to the combined effect of multipath and shadowing. III. NUMERICAL RESULTS This section presents numerical values for the probability density function, LCR, and AFD which define the model given in this paper. In addition, a comparison with similar functions computed from measured data is discussed. Equations (6),(7), and (8) were used to calculate the probability density function, p(r). Many calculations with different values for bo, do, and p were carried out with the objective of fitting results from application of the model to those derived from measurements made. by colleagues on a simulated satellite path. During these measurements a helicopter was flown at a fixed elevation angle of 5 with respect to the receiver [ll]. The measurement site was a rural area with about 35 percent tree cover. Numerical evaluation of the probability density function in (6) was time consuming, but allowed computation of the LCR and the AFD using (25) and (26). Fig. shows a comparison of the cumulative distribution function for the received envelope calculated using (8) and that computed from measurements for infrequent light shadowing (sparse tree cover), frequent heavy shadowing (dense tree cover) and overall results (the combined results of light and heavy shadowing). For the case of infrequent light shadowing, the model shows the best fit around the median region and some deviation near the tails of the distribution. Generally, results of the model showed a slightly higher shadowing effect than those from measured data. For the case of frequent heavy shadowing, the results of the model show reasonably good

4 LOO: LAND MOBILE SATELLITE LINK ' I l l I I I I I I I I I I I I I PERCENT OF TIME RECEIVED SIGNAL IS GREATER THAN ORDINATE Fig.. A comparison of measured and calculated values of the probability distribution of signal, Conditions Infrequent light shadowing Frequent heavy shadowing Overal I results TABLE I MODEL PARMETERS Standard Deviation logio(/bo) Mean loglo(p) Multipath Power loglo (bo) -a. -2. agreement with measured values throughout the fading range. For the combined results, the fit was poor about the median but reasonably good in the weak signal range which is most important for fade margin calculations. The model parameters were obtained by trial and error to fit measured values. These parameters are given in Table I. Fig. 2 shows values for the LCR calculated using (6) and (25) for the case of infrequent light shadowing. The results show that the maximum LCR occurs when there is no correlation between the rate of change of the envelope due to multipath and that due to shadowing. Measured values [l for this case are also shown. These results indicate that the model gives a good indication of the LCR, when high correlation, p =.5 to.9, between multipath and shadowing is assumed. The figure also indicates that there is a bound on the LCR when the signal is large. Fig. 3 shows a plot of the AFD calculated from the model and from measurements [ll]. Application of the model results in a good approximation to the AFD anticipated on a mobile satellite link. For low signal -6 w t..o w B Y d. 3 5 P... CALCULATED - MEASURED I I I I I I Fig. 2. Level crossing rate-infrequent light shadowing. i

5 26. - IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY, VOL. VT-3, NO. 3, AUGUST CALCULATED - MEASURED 2 s.- 9 CI 2 e. - I I' I I I I Fig. 3. Average fade duration-infrequent light shadowing. 'o""l. lo Fig. 5. Average fade duration-frequent heavy shadowing..o g v) 8 g. J E 5. LEVEL RELATIVE TO LINE-OF-SIGHT LEVEL (db) Fig.. Level crossing rate-frequent heavy shadowing..' I I I I I I Fig. 6. Level crossing rate-overdl results.

6 LOO: LAND MOBILE S.TELLITE LIKK 27 ~ Y Y a u w a - I I I I I I I Fig. 7. Average fade duration-overall results. levels, the two curves compare reasonably well. These results also show that the best fit is obtained when p is large. Fig. shows a plot of normalied LCR versus relative signal level for the case of frequent heavy shadowing. Measured values are also shown. The values calculated from the model for p between.5 and.8 bound the values calculated from measured data throughout the relative signal level range. Fig. 5 shows AFD for the case of frequent heavy shadowing; the calculated values, with p =.7 give good agreement when compared with measurements. This indicates that the correlation between the rate of change of envelope due to multipath and foliage attenuation is high for heavy shadowing as well as for light shadowing. Figs. 6 and 7 show the LCR and AFD for the case of the combined light and heavy shadowing with correlation coefficient, p as a parameter. There were no measured values available for this case. Using measured values for the previous two cases as a guide, however, one can estimate that the measured curve would be closest to the case of p =.8, but this remains to be verified. A comparison between the probability distribution of signal level calculated from (8) and that calculated from measured data for the two extreme cases, infrequent light shadowing and frequent heavy shadowing, has been described. It has been found that the model can provide a good indication of LCR and AFD throughout the range of relative signal levels indicated in the figures. The model should thus be useful for predicting the performance of communication systems. IV. CONCLUSION A statistical model for a land mobile satellite link has been described in the paper. Expressions for the LCR and AFD were obtained. Comparisons between values calculated from the model and those calculated from measured data show reasonably good agreement. It is also shown that the assumption of high correlation between the effects of multipath and those caused by shadowing yields results that compare well with values derived from measurements. This is true for both light and heavy shadowing. The model should be useful in estimating the performance of communications systems with various modulation schemes and it should also facilitate the development of a propagation simulator. REFEREKES S. H. Lin, Statistical behavior of a fading signal, Bell Syst. Tech. J., vol. 5, no., pp , Dec. 97. G. L. Turin et al., A statistical model of urban multipath propagation, IEEE Trans. Veh. Technol., vol. VT-2, no., pp. -9, Feb H. Suuki, A statistical model for urban radio propagation, IEEE Trans. Commun., vol. COM-25, no. 7, pp , July 977. F. Hansen and F. I. Meno, Mobile fading-rayleigh and lognormal superimposed, IEEE Trans. Veh. Technol., vol. VT-26, no., pp , NOV R. C. French, The effect of fading and shadowing on channel reuse in mobile radio, IEEE Trans. Veh. Technol., vol. VT-28, no. 3, pp. 7-8, Aug T. Aulin, A modified model for the fading signal at a mobile radio channel, IEEE Trans. Veh. Technol., vol. VT-28, no. 3, pp , Aug J. Zander, A stochastical model of the urban UHF radio channel, IEEE Trans. Veh. Technol., vol. VT-3, no., pp. 5-55, Nov C. Jakes, Jr.. Ed., MicrowaveMobile Communications. New York: iley, 97, pp C. Y. Lee, Mobile Communications Engineering. New York: McGraw-Hill, 982, pp , P. Beckmann, Probability in Communication Engineering. New York: Harcourt, Brace and orld, Inc., 967. pp J. S. Butterworth and E. E. Matt, The characteriation of propagation effects for land mobile satellite services, Inst. Elec. Eng. Cony. Proc. Satellite Syst. Mobile Commun. and Navigations, June 983. S.. Rice, Mathematical analysis of random noise, Bell Syst. Tech. J., vol. 23, pp , July 9 and vol. 2, pp. 6-56, Jan , Statistical properties of a sine wave plus random noise, Bell Syst. Tech. J., vol. 27, pp. 9-7: Jan. 98. A. Papoulis, Probability, Random Variables and Stochastic Processes. New York: McGraw-Hill, 965, pp Chun Loo ( 62) received the B. Eng. degree in electrical engineering from McGill University, Montreal, PQ, Canada, in 962, and the M.Sc. degree in electrical engineering from the University of Ottawa, Ontario, Canada, in 966. From 962 to 965, he was employed by Computing Devices of Canada, Ottawa, where he participated in the development of avionic equipment. He was with Northern Electric Research and Development Labs (now Bell Northern Research) from 966 to 97 where he was engaged in the theoretical studies of FM distortions, FM threshold characteristics and in computer simulation of communications systems. Since 97, hehasbeen with the Communications Research Centre, Department of Communications, Ottawa. He did systems analysis under the CTS (Hermes) project a forerunner of direct broadcast satellite, and the ANIK B Pilot project. He is the author of several papers relating to FM and digital transmission through a nonlirlear satellite channel. He is currently doing system studies on various modulation schemes and on channel characteriation for mobile satellite communications applicable to the MSAT program.