








Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Community
Ask the community for help and clear up your study doubts
Discover the best universities in your country according to Docsity users
Free resources
Download our free guides on studying techniques, anxiety management strategies, and thesis advice from Docsity tutors
Clarke model for multipath and working
Typology: Lecture notes
1 / 14
This page cannot be seen from the preview
Don't miss anything!
Abstract Clarke’s classical model of mobile radio reception assumes isotropic rich scattering around the mobile receiver antenna. The assumption of isotropic scattering is valid only in limited circumstances. In this contribution we develop a generalized Clarke model, which is applicable to mobile radio reception in general scattering environments. We give expressions for the autocorrelation and power spectral density (PSD) of the channel fading process and demonstrate the generality of the model by applying it to different non-isotropic scattering scenarios. Using the generalized model, we analyze the effect of mobile direction of travel and the non-isotropicity on the statistics of the channel fading process. We also show that if the mobile direction of travel is equiprobable in all directions, a non-isotropic scattering environment on average is as good as an isotropic scattering environment.
Index Terms Non-Isotropic scattering, Rayleigh fading, Autocorrelation, Power Spectral Density, Power Azimuth Spec- trum, mobile fading channel
I. INTRODUCTION In real world communication scenarios, the transmitter and/or the receiver may be in motion. In a mobile- radio situation in which the transmitter is fixed in position while the receiver is moving, the direct line between the transmitter and receiver may be obstructed by buildings. At ultra high frequencies and above, therefore, the mode of propagation of the electromagnetic energy from transmitter to receiver is largely by way of scattering [1]. The amplitude fluctuations of the received signal have been shown to follow Rayleigh distribution in this communication scenario. By assuming a 2-dimensional (2D) uniform power azimuth spectrum (PAS) which corresponds to uniform distribution of angle-of-arrival (AOA) around the receive antenna, the autocorrelation of the received signal was shown to be strictly real valued and the resulting power spectral density (PSD) was U-shaped symmetric [1]. It has been argued and experimentally demonstrated (see [2] and references therein) that the scattering encountered in many suburban and rural environments is non-isotropic i.e., the distribution of AOA of waves is not uniform as assumed in [1]. The use of a directional antenna with non-uniform gain pattern at the receiver also results in non-isotropic PAS as seen by the antenna.
The authors are with the Applied Signal Processing group, Department of Information Engineering, Research School of Information Sciences and Engineering, the Australian National University, Canberra ACT 0200, Australia, email:{rauf.iqbal, thushara.abhayapala, tharaka.lamahewa}@anu.edu.au.
The assumption of uniform PAS resulted in closed form expressions for the statistics of the channel in [1]. For non-uniform PAS, the resulting expressions for the statistics of the channel are not in closed form [3]. A quadratic form for the probability distribution function (pdf) of the AOA proposed in [4] resulted in a closed form expression for the correlation of the complex envelope of the received signal in a non-isotropic scattering environment. In [5], von Mises distribution was used to model the pdf of AOA and effects of non-isotropic scattering on the correlation properties and velocity estimation in a Rician fading channel were discussed. In practice, PAS under non-isotropic conditions has also been shown to be well modelled by truncated Laplacian, truncated Cosine, and truncated Gaussian distribution [6]. Some of the other work related to mobile fading channels are [7]–[9]. The assumption of uniform PAS introduces small errors on the first order statistics of the received signal but a significant error on the second order statistics [10], like correlation function or, equivalently, PSD, and level crossing rates or, equivalently, the fading rate. There are certain communication system parameters like the estimation of vehicle velocity [5] for handoff decisions and the achievable information rates [11] without channel state information (CSI) that depend on the correlational properties of the received signal. It is, therefore, of some interest to develop a model that accurately models the statistics of the mobile Rayleigh fading channel in isotropic and non-isotropic scattering environments. In this sequel we use the term generalized Rayleigh fading to denote Rayleigh fading in general scattering environments. The main contribution of this paper is the generalization of the Clarke model of mobile radio reception to generalized Rayleigh fading. This model can be used for the accurate performance prediction and evaluation of a communication system in generalized Rayleigh fading. Using the proposed generalized model, the impact on channel statistics of different parameters like the mobile direction of travel, mean angle of arrival (AOA) and the angular spread of the scattering environment can be accurately determined. The rest of the paper is organized as follows. In Section II, we develop a generalized discrete time Rayleigh fading model and derive expressions for the autocorrelation and the PSD. In Section III, we show that the Clarke’s isotropic scattering model is a special case of this generalized model. In Section III, we also analyze the effects of non-isotropicity and mobile velocity on the channel statistics. Some of the possible applications of the model are described in Section IV. Finally, the conclusions are drawn in Section V. Throughout the paper, the following notation will be used: Bold lower (upper) letters denote vectors (matrices). The superscript ∗^ denotes the conjugate transpose. <(·) represents the real part operator. The notation E {·} denotes the mathematical expectation and the symbol ∮^ represents integration over a circle.
II. CHANNEL MODEL We consider a downlink transmission system where the transmitter is stationary while the receiver is moving with some speed |v| at an angle of φv with respect to the x-axis, where v ≡ (|v|, φv ) is the velocity of the mobile as shown in Fig. 1. We assume that the scatterers are distributed in the far-field from the transmitter and receiver antennas. We also assume that the channel between the transmitter and the mobile receiver is a strictly bandlimited, flat-fading (frequency non-selective), wide-sense stationary, zero mean circularly symmetric Gaussian fading process.
continuous-time and discrete-time models ((1) and (2) respectively) are related through the following:
x[j] = √^1 T s
∫ (^) (j+1)Ts jTs x(t) dt, (3)
y[j] = √^1 Ts
∫ (^) (j+1)Ts jTs y(t) dt, (4)
n[j] = √^1 T s
∫ (^) (j+1)Ts jTs n(t) dt, (5)
hd[j] = √∫ (^) Ts^1 0
∫ (^) Ts 0 Φ(s^ −^ t)^ dt ds
∫ (^) (j+1)Ts jTs hc(t) dt, (6)
where Φ(·) in (6) is the covariance of the continuous channel process, hc(t). Equations (3)−(6) are based on the facts that matched filter is an integrator, the output of the filter is sampled at the end of the each symbol interval Ts and discrete-time input, output, noise and channel processes must be normalized to have the same (co)variances as their continuous-time counterparts. It is to be noted that while most of the material on channel modeling deals with continuous-time model, we have chosen to work in the discrete-time domain, firstly, because the two domains are equivalent as long as the matched filter output is sampled at least at the Nyquist rate corresponding to rate of channel variation (for an idea of equivalence between continuous-time and discrete-time models please see [17]). Secondly, whenever it comes to practical implementation of the model, discrete-time models are computationally more efficient. There are some instances in literature of discrete-time modeling e.g., see [18] and references therein. Equation (2) can be obtained as a special case of the discrete-time triply (i.e., time-frequency-space) selective model of [18]. With no loss of generality, let the mobile be at some arbitrary point ‘O’ at the signaling instant j′. Thus, at the signaling interval j = j′^ + k, the mobile will be at the point (|v|ηkTs, φv ) with respect to ‘O’, where Ts is the symbol duration, η = 2π/λ is the free space phase constant, λ is the wavelength, and k is an integer, which represents the time lag. Assume that the scattered signals are impinging on the mobile receiver from all directions on the 2D (horizontal) plane. Let f (β) be the scattering gain at the origin due to the signals arriving from the direction β with respect to the x-axis (see Fig. 1). Then we write the channel gain as
h[j] =
f (β) exp
iηjTs v · βˆ
dβ, (7)
where βˆ ≡ (1, β) represents a unit vector along the direction β, ‘·’ is the scalar product between two vectors and i = √− 1. Note that the factor exp(iηjTsv · βˆ) in ((7)) reflects the phase delay of the incoming signal from the direction β at the mobile receiver with respect to the origin. It is crucial to highlight an important conceptual difference between our approach and that used in [1] for modelling mobile radio reception. In [1], firstly, a probability of arrival of waves is associated with each direction in the azimuth. Secondly, the complex scattering gain, f (β), from a certain direction of arrival is also random. In this contribution, we essentially assume that there is no probability distribution associated with AOA i.e., waves are assumed to be impinging on the mobile receive antenna from all the directions in the azimuth. Only the complex random scattering gain, f (β), is assumed to be random and has, associated with it, a probability distribution. An azimuth direction with zero associated probability of arrival of wave is equivalent to a direction with waves assumed to be impinging from that particular direction but zero complex scattering gain. In other
(j′)
v
(j)
O x
β φv
plane wave
Fig. 1. Illustration of the key parameters: direction of mobile travel φv , mobile velocity v, direction of signal (wave) arrival β, time instances j′^ and j, and the origin O of the co-ordinate system.
words, our approach is a simplified yet equivalent form of that employed in [1] but proves more convenient, as we would see, for arriving at a generalized Rayleigh fading model for mobile radio reception.
A. Autocorrelation of the Channel Fading Process
Using (7) we write the correlation between the channel gain at the signaling intervals j′^ and j as Φ(j, j′) = E {hj h∗ j′^ } =
E {f (β)f ∗(β′)} exp
iηTs(j v · βˆ − j′^ v · βˆ′
dβdβ′, (8)
where E {·} stands for mathematical expectation. Assuming that the scattering gain from two distinct directions are uncorrelated and are zero mean [3], i.e.,
E {f (β)f ∗(β′)} = E
|f (β)|^2
δ(β − β′). (9)
where δ is the Dirac delta function. Making use of (9) in (8) , we can write
Φ(j, j′) = Φ(j − j′) = Φ(k) =
Ψ(β) exp
iηTsk v · βˆ
dβ (10)
where k = j − j′^ and
Ψ(β) = E {|f (β)|^2 }^ , (11)
normally termed as the angular power distribution of the received signal. Thus Ψ(β) is the average power received from the direction β.
where ω ∈ [−π, π] is the continuous radian frequency variable. We simplify (16) in Appendix A to obtain
Φ(ω) = (^) ω^1 D
m=−∞
γm eimφv^ Fm
( (^) ω ωD
where ωD = ωdTs is the maximum Doppler spread normalized by the symbol rate, and fD = ωD/ 2 π is the normalized maximum Doppler frequency which is also called the normalized fading rate. It is easy to see that (17) is equivalent to the following
Φ(ω) = (^) ω^1 D
1 − (ω/ωD )^2 + 2
m=
< (γm eimφv^ )^ Fm
( (^) ω ωD
where <(·) is the real part of the argument. Equation (18) clearly shows that the PSD is real-valued as expected.
C. Truncation of series expansions
Equations (15) and (18) involve summation over infinite number of terms. However it is almost always possible in cases of practical interest to safely truncate the respective series up to a finite number of terms in the light of the following two facts: i) For a fixed order m, Jm(x) starts small and reaches its maximum at argument x ≈ O(m) before it starts decaying slowly. It was shown in [22] that Jm(x) ≈ 0 for |m| > 2 dx/ 2 e+1 with e = 2. 7183.. .. Since the argument of Jm(kωD) depends on the lag variable k, the approximation would also depend on k for a fixed normalized fading rate. ii) The Fourier coefficients γm must decay with m for Fourier series to be convergent (see e.g., [23]). The rate of decay depends on the smoothness of the function which is, in turn, related to the number of continuous derivatives of the function. In fact, the Fourier coefficients of an analytic (infinitely differentiable) function decay exponentially with m. Fourier coefficients of a Gaussian distribution, for example, decay exponentially with m, and decay polynomially for Laplace distribution [24]. The Fourier coefficients for a uniform distribution decay as 1 /m. In all these cases, the Fourier coefficients tend to zero^5 with m. The rapid the decay of the Fourier series coefficients, the less the number of Fourier modes (γm) with significant contribution and vice versa. For small k, the fact (i) above is useful in approximating (15) due to the presence of Jm(kωD). When k = 0, for example, Jm(kωD) = 1 for m = 0 and zero otherwise. In other words, none of the Fourier modes γm except m = 0 contribute to the autocorrelation of the channel process. For k = 1, Jm(kωD) ≈ 0 after 6 |m| > (2dkωD/ 2 e + 1). As k increases, the number of Fourier coefficients that are “allowed” to contribute also increases. In the limit k → ∞, the number of Fourier modes that contribute to (15) also approaches infinity. This is where the fact (ii) comes into effect. Due to the decay (e.g., exponential for a Gaussian distribution) of γm with m, the values of γm become increasingly small so that there must exist some finite m 0 such that γm ≈ 0 for all |m| > m 0. We can , therefore, truncate the infinite summation in (15) to a summation over |m| = m 0 terms.
(^5) The Delta distribution is an exception to this behavior. The Fourier coefficients of a Delta distribution do not tend to zero at all, indicative of the fact that it is not an ordinary function and its Fourier series does not converge in the standard sense. (^6) For a normalized fading rate fD = 0. 05 and maximum time lag of k = 10 symbols, this approximation requires only |m| = 3 terms.
∆r ∆r
φv x
β 0 v
O
Fig. 2. Illustraion of an uniform-limited scattering scenario where the scattered power is uniformly distributed with magnitude 1 /2∆r over a part of the azimuth with a mean angle β 0 and a maximum deviation of ∆r on each side of the mean. The direction of the mobile travel φv is also shown.
The approximation of the PSD of the fading process in (18) seems difficult due to the presence of the factor^7 Fm(ω/ωD) ((27) with x = ω/ωD) which, for all m, sharply increases as ω approaches ωD becoming infinite at ωD. The approximation, therefore, would depend on the the decay of the Fourier coefficients of a particular scattering distribution. As long as the decay of γm is sufficiently fast, we can use the fact (ii) to approximate (18) to some finite |m| = m′ 0 with the approximation error that approaches zero as m 0 approaches infinity. The accuracy of this approximation, moreover, depends on the angular spread^8 of the scattering distribution and the mobile direction of travel.
III. EFFECT OF NON-ISOTROPICITY AND MOBILE VELOCITY ON CHANNEL STATISTICS Equations (15) and (18) give the second order channel statistics, namely autocorrelation and power spectral density (PSD) for any 2D scattering environment around the receiver. For illustration purposes, throughout this section, we consider the so-called uniform limited scattering scenario, i.e., scattered waves are arriving uniformly from an angular sector as shon in Fig. 2. For this case, we have
F (β) =
1 /(2 (^4) r ), if |β − β 0 | ≤ 4r 0 , otherwise
where β 0 is the mean angle of arrival, β is any other angle of arrival, and (^24) r is the angular spread of the uniform limited scattering arrival signals as illustrated in Fig. 2. The values of γm for this distribution were derived in closed form in [?] and are given as
γm = exp (−imβ 0 ) sinc(m∆r ). (20)
(^7) Notice that the Fourier transformation of Jm(kωD) in (16), that gave rise to the factor Fm(ω/ωD ), involves summation over k = ±∞ which implies that infinite number of modes are allowed to contribute to the PSD of the fading process. (^8) The angular spread, Λ, is defined as the standard deviation of scattering distribution (PAS).
Using (20) in (15) and (18), we can write the autocorrelation and the PSD respectively, as Φ(k) =
m=−∞
im^ sinc(m∆r ) Jm(ωD k) exp (−im (β 0 − φv )) , (21a)
Φ(ω) = (^) ω^1 D
m=−∞
sinc(m∆r ) Fm
( (^) ω ωD
exp (−im (β 0 − φv )). (21b) In the rest of this section, we first show that the classical Clarke model is a special case of the generalized Clarke model introduced in this contribution. We then explore the effect of the non-isotropicity and mobile velocity on channel autocorrelation and PSD.
A. Clarke’s model as a special case
When ∆r = π, i.e., when the scattered power is uniformly distributed over the full azimuth plane around the mobile receiver, it can easily be verified that (21a) and (21b) collapse to the following equations for autocorrelation and PSD, respectively,
Φ(k) = J 0 (ωDk), (22) Φ(ω) = (^) ωD (1 − 2 (ω/ωD) (^2) ) , (23)
which are the well-known Clarke’s model [1] for 2D isotropic scattering around the receive antenna. Thus Clarke’s model is a special case of the generalized model developed in this contribution. Notice that, in general, the autocorrelation in (21a) is complex valued unlike that given in (22) for isotropic scattering environment which is strictly real valued.
B. Effect of mobile velocity
For a uniform-limited scattering scenario, the effect of changing the mobile direction of travel on the auto- correlation and PSD of the received signal has been shown in Figures 3 and 4 respectively. The autocorrelation and PSD for isotropic case have also been plotted for comparison. A marked deviation from the isotropic case can be observed. The skewness of the PSD is easily observed: If the mobile is moving into the non-isotropic scattering environment, the Doppler spectrum becomes (emphasized and) concentrated towards positive Doppler frequency axis. On the other hand, the Doppler spectrum is skewed towards negative Doppler frequency axis if the mobile moves away from the scatterers. The spectrum is symmetric about the mean only when the mobile moves at right angles to the mean scattering angle. The above discussion of the autocorrelation and PSD implicitly assumes that the direction of mobile travel is perfectly known which is usually not the case in practice. It may be of some interest to find out the autocorrelation and PSD if the mobile direction of travel is unknown at the receiver. Suppose the mobile direction of travel is equiprobable in all directions i.e., p (φv ) = 1/ 2 π, then, from (21a), the average autocorrelation, Φavg(k) is given by
Φavg(k) =
m=−∞
im^ sinc(m∆r ) Jm(ωDk) exp (−imβ 0 )
∫ (^2) π 0 exp (imφv ) dφv , (24)
and using (21b), the average PSD, Φavg(ω) , is given by
Φavg(ω) = (^) ω^1 D
m=−∞
sinc(m∆r )Fm
( (^) ω ωD
exp (−imβ 0 )
∫ (^2) π 0 exp (imφv ) dφv , (25)
It is not hard to see that, irrespective of ∆r and the β 0 , the integrals in (24) and (25) are zero for all m 6 = 0, and these two equations converge, respectively, to (22) and (23) i.e., the Clarke’s isotropic case. In other words, if the mobile direction of travel is equiprobable in all directions, a non-isotropic scattering environment on average is as good as an isotropic scattering environment.
C. Effect of non-isotropicity
For a fixed direction of mobile travel, the effect of changing the degree of non-isotropicity on the autocorre- lation and PSD of the fading process has been plotted in Figures 5 and 6 respectively. It can be seen from Fig. 5 that the channel fading process can have significantly higher correlation over time in non-isotropic scattering environments as compared to the isotropic scattering environment. With increasing ∆r (or, in other words, decreasing the non-isotropicity) the correlation curves tend towards those of isotropic case. Figure 6 shows the effect of changing the degree of non-isotropicity on the PSD of the channel process. The normalized Doppler spread seems to be directly proportional to ∆r. For a fixed carrier frequency, the normalized fading rate fD depends directly on the mobile speed. In our case none of the parameters, except the scattering environment (∆r ), is being changed. It can, therefore, be concluded that changing the degree of non-isotropicity is actually equivalent to changing the normalized fading rate to some effective value, fDeff. Since the second order statistics are directly related to the normalized fading rate, this verifies the point of view of [10] that the assumption of isotropic distribution of scattered power in a non-isotropic distribution of power introduces significant errors in the second order statistics.
IV. APPLICATIONS OF THE GENERALIZED CLARKE MODEL The generalized Clarke model presented in this contribution can be used for accurate performance evaluation and prediction of a communication system in Rayleigh fading under any scattering environment. Among many possible applications, following are some important applications of the generalized model:
APPENDIX A PROOF OF (17) The following identity exists for continuous time Fourier transform (CTFT) of the Bessel function of the first kind and integer order μ [3] ∫ (^) ∞ −∞ Jμ(ωd∆t) e−iω∆td∆t =
Fμ
( (^) ω ωd
iμωd^ ,^ (26) where
Fμ(x) , 2 cos^
(μ cos− (^1) (x)) √ 1 − x^2.^ (27)
We know from the basic Fourier theory that the DTFT of a sampled process is essentially a magnitude and frequency scaled version of CTFT of the continuous-time process with 2 π periodicity [20]. Making use of (26) and (27), the DTFT of the sampled Bessel function of the first kind and integer order μ is given by
∑^ ∞ p=−∞
Jμ(ωdpTs) e−iωpTs^ = (^) T^1 s
`=−∞
Fμ
ωdTs (ω^ −^2 π`)
iμωd^ ,
=
`=−∞
Fμ
ωD (ω^ −^2 π`)
iμωD^ ,^ (28) where ωD = ωdTs is the maximum Doppler spread normalized by the symbol rate. Since we are working at the baseband level, only the low pass term (` = 0) in (28) is of interest to us, i.e., ∑^ ∞ p=−∞
Jμ(ωdpTs) e−iωpTs^ =
Fμ
( (^) ω ωD
iμωD^.^ (29) Using (27) and (29), we can rewrite (16) to get the desired expression (17) for PSD
Φ(ω) = (^) ω^1 D
m=−∞
γm eimφv^ Fm
( (^) ω ωD
[1] R. H. Clarke, “A statistical theory of mobile-radio reception,” in Bell System Technical Journal, vol. 47, 1968, pp. 957–1000. [2] A. Abdi, J. A. Barger, and M. Kaveh, “A parametric model for the distribution of the angle of arrival and the associated correlation function and power spectrum at the mobile station,” vol. 51, no. 3, pp. 425–434, May 2002. [3] J. S. Sadowsky and V. Kafedziski, “On the correlation and scattering functions of the WSSUS channel for mobile communications,” in IEEE Transactions on Vehicular Technology, vol. 47, Feb 1998, pp. 270–282. [4] K. Anim-Appiah, “Complex envelope correlations for non-isotropic scattering,” IEE Electronics Letters, vol. 34, no. 9, pp. 1039–1052, Apr 1998. [5] C. Tepedelenlioglu and G. Giannakis, “On velocity estimation and correlation properties of narrow-band communication channels,” IEEE Trans. Vehic. Technol., vol. 50, no. 4, pp. 1039–1052, July 2001. [6] P. D. Teal, T. D. Abhayapala, and R. A. Kennedy, “Spatial correlation for general distributions of scatterers,” IEEE Signal Processing Letters, vol. 9, no. 10, pp. 305–308, Oct. 2002. [7] T. Aulin, “A modified model for the fading signal at a mobile radio channel,” IEEE Trans. on Vehicular Technology, vol. VT-28, no. 3. [8] J. D. Parsons, The mobile radio propagation channel. London, UK: Pentech Press, 1992. [9] J. D. Parsons and A. M. D. Turkmani, “Characterisation of mobile radio signals: model description,” in IEE Proceedings-I, vol. 138, no. 6, 1991, pp. 549–556. [10] W. C. Y. Lee, “Finding the aproximate angular probability density function of wave arrival by using a directional antenna,” in IEEE transactions on Antennas and Propagation, vol. AP-21, 1973, pp. 328–334. [11] X. Deng and A. M. Haimovich, “Information rates of time varying Rayleigh fading channels,” in Proc. IEEE International Conference on Communication , ICC’2004, Paris, France, 2004. [12] D. Torrieri, Principles of Spread-Spectrum Communication Systems. Springer, 2005. [13] D. Tse and P. Viswanath, Fundamentals of wireless communication. New York, NY, USA: Cambridge University Press, 2005. [14] F. D. Neeser and J. L. Massey, “Proper complex random processes with applications to information theory,” vol. 39, pp. 1293–1302, Jul. 1993. [15] T. A. Lamahewa, T. D. Abhayapala, R. A. Kennedy, and J. T. Y. Ho, “Space-time cross correlation and space-frequency cross spectrum in non-isotropic scattering environment,” International Conference on Acoustics, Speech, and Signal Processing, vol. 4, May 2006. [16] J. R. Barry, D. G. Messerschmitt, and E. A. Lee, Digital Communication: Third Edition. Norwell, MA, USA: Kluwer Academic Publishers, 2003. [17] V. K. N. Lau and Y. R. Kwok, Channel-Adaptive Technologies and Cross-Layer Designs for Wireless Systems with Multiple Antennas: Theory and Applications (Wiley Series in Telecommunications and Signal Processing). Wiley-Interscience, 2006. [18] C. Xiao, J. Wu, S. Leong, Y. R. Zheng, and K. Letaief, “A discrete-time model for triply selective MIMO Rayleigh fading channels,” IEEE Trans. on Wireless Comm., vol. 3, no. 5, pp. 1678–1688, Sep. 2004. [19] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory, 2nd ed. New York: Springer, 1998, vol. 93. [20] A. V. Oppenheim and R. Shafer, Eds., Discrete-Time Signal Processing, 2nd ed. New Jersey: Prentice Hall, 1989. [21] D. C. Champeney, A Handbook of Fourier Transforms. U.K.: Cambridge University Press, 1987. [22] H. M. Jones, R. A. Kennedy, and T. D. Abhayapala, “On dimensionality of multipath fields: Spatial extent and richness,” in Proc. IEEE Int. Conf. Acoust., Speech, Signal Processing, ICASSP’02, vol. 3, Orlando, Florida, May 2002, pp. 2837–2840. [23] G. B. Folland, Real Analysis, 2nd ed. New York: Wiley-Interscience, 1972. [24] F. Comte, Y. Rozenholc, and M. L. Taupin, “Penalized contrast estimator for density deconvolution,” in The Candian Journal of Statistics, vol. 34, no. 3, 2006. [25] A. A. Farid, Z. Q. Luo, and Z. Ding, “Blind channel equalization based on second order statistics,” International Conference on Acoustics, Speech, and Signal Processing, vol. 3, pp. 557–560, 2005. [26] L. Tong, G. Xu, B. Hassibi, and T. Kailath, “Blind channel identification based on second-order statistics: A frequency-domain approach,” IEEE Trans. on Inform. Theory, pp. 329–334, 1995.