Spherical Harmonic Expansion of Fisher-Bingham Distribution and 3D Spatial Fading Correlation for Multiple-Antenna Systems Yibeltal F
Total Page:16
File Type:pdf, Size:1020Kb
IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY SUBMISSION 1 Spherical Harmonic Expansion of Fisher-Bingham Distribution and 3D Spatial Fading Correlation for Multiple-Antenna Systems Yibeltal F. Alem, Student Member, IEEE, Zubair Khalid, Member, IEEE and Rodney A. Kennedy, Fellow, IEEE Abstract—This paper considers the 3D spatial fading correla- closed-form expressions for evaluating the SFC function have tion (SFC) resulting from an angle-of-arrival (AoA) distribution been developed in the existing literature (e.g., [10], [16]–[25]). that can be modelled by a mixture of Fisher-Bingham distribu- The elliptic (directional) nature of the Fisher-Bingham tions on the sphere. By deriving a closed-form expression for the spherical harmonic transform for the component Fisher-Bingham distribution offers flexibility in modelling the distribution of distributions, with arbitrary parameter values, we obtain a normalized power or AoA of the multipath components for closed-form expression of the 3D-SFC for the mixture case. The most practical scenarios. Exploiting this fact, a 3D spatial 3D-SFC expression is general and can be used in arbitrary multi- correlation model has been developed in [9], where the dis- antenna array geometries and is demonstrated for the cases of tribution of AoA of the multipath components is modelled by a 2D uniform circular array in the horizontal plane and a 3D regular dodecahedral array. In computational aspects, we use a positive linear sum of Fisher-Bingham distributions, each recursions to compute the spherical harmonic coefficients and with different parameters. Such modelling has been shown give pragmatic guidelines on the truncation size in the series to be useful in a sense that it fits well with the multi-input representations to yield machine precision accuracy results. The multi-output (MIMO) field data and allows the evaluation results are further corroborated through numerical experiments of the SFC as a function of the angular spread, ovalness to demonstrate that the closed-form expressions yield the same results as significantly more computationally expensive numerical parameter, azimuth, elevation, and prior contribution of each integration methods. cluster. Although the SFC function presented in [9] is general in a sense that it is valid for any arbitrary antenna array Index Terms—Fisher-Bingham distribution, spherical har- monic expansion, spatial correlation, MIMO, angle of ar- geometry, it has not been expressed in closed-form and has rival (AoA). been only evaluated using numerical integration techniques, which can be computationally intensive to produce sufficiently accurate results. If the spherical harmonic expansion of the I. INTRODUCTION Fisher-Bingham distribution is given in a closed-form, the The Fisher-Bingham distribution, also known as the Kent SFC function can be computed analytically using the spherical distribution, belongs to the family of spherical distributions harmonic expansion of the distribution of AoA of the multipath in directional statistics [1]. It has been used in applications components [21]. To the best of our knowledge, the spherical in a wide range of disciplines for modelling and analysing harmonic expansion of Fisher-Bingham distribution has not directional data. These applications include sound source lo- been derived in the existing literature. calization [2], [3], joint set identification [4], modelling protein In the current work, we present the spherical harmonic structures [5], 3D beamforming [6], classification of remote expansion of Fisher-Bingham distribution and a closed-form sensing data [7], modelling the distribution of AoA in wireless expression that enables the analytic computation of the spher- communication [8], [9], to name a few. ical harmonic coefficients. We also address the computational arXiv:1501.04395v1 [cs.IT] 19 Jan 2015 In this work, we focus our attention to the use of Fisher- considerations required to be taken into account in the evalua- Bingham distribution for modelling the distribution of angle tion of the proposed closed-form. Using the proposed spherical of arrival (AoA) (also referred as the distribution of scatterers) harmonic expansion of the Fisher-Bingham distribution, we in wireless communication and computation of spatial fading also formulate the SFC function experienced between two correlation (SFC) experienced between elements of multiple- arbitrary points in 3D-space for the case when the distribution antenna array systems. Modelling the distribution of scatterers of AoA of the multipath components is modelled by a mixture and characterising the spatial correlation of fading channels (positive linear sum) of Fisher-Bingham distributions. The is a key factor in evaluating the performance of wireless SFC presented here is general in a sense that it is expressed as communication systems with multiple antenna elements [10]– a function of arbitrary points in 3D-space and therefore can be [15]. It has been an active area of research for the past two used to compute spatial correlation for any 2D and 3D antenna decades or so and a number of spatial correlation models and array geometries. Through numerical analysis, we also validate The authors are with the Research School of Engineering, College the correctness of the proposed spherical harmonic expansion of Engineering and Computer Science, The Australian National Univer- of Fisher-Bingham distribution and the SFC function. In this sity, Canberra, ACT 2601, Australia. Emails: fyilbeltal.alem, zubair.khalid, paper, our main objective is to employ the proposed spherical [email protected]. This work was supported by the Aus- tralian Research Council’s Discovery Projects funding scheme (Project no. harmonic expansion of the Fisher-Bingham distribution for DP150101011). computing the spatial correlation. However, we expect that IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY SUBMISSION 2 ` the proposed spherical harmonic expansion can be useful in where dm;m0 (·) denotes the Wigner-d function of degree ` and various applications where the Fisher-Bingham distribution is orders m and m0 [26]. used to model and analyse directional data (e.g., [2]–[7]). By completeness of spherical harmonics, any finite energy The remainder of the paper is structured as follows. We function f(xb) on the 2-sphere can be expanded as review the mathematical background related to signals defined 1 ` on the 2-sphere and spherical harmonics in Section II. In X X m m f(xb) = (f)` Y` (xb); (4) Section III, we define the Fisher-Bingham distribution and `=0 m=−` present its application in modelling 3D spatial correlation. where (f)m is the spherical harmonic coefficient given by We present the spherical harmonic expansion of the Fisher- ` Z Bingham distribution and analytical formula for computing m m (f)` = f(xb)Y` (xb) ds(xb): (5) the spherical harmonic coefficients in Section IV, where we 2 S also address computational issues. We derive a closed-form The signal f is said to be band-limited in the spectral domain expression for the 3D SFC between two arbitrary points in 3D- m at degree L if (f)` = 0 for ` > L. For a real-valued function space when the AoA of an incident signal follows the Fisher- f(x), we also note the relation Bingham probability density function (pdf) in Section V. b −m m m In Section VI, we carry out a numerical validation of the (f)` = (−1) (f)` ; (6) proposed results and provide examples of the SFC for uniform which stems from the conjugate symmetry property of spher- circular array (2D) and regular dodecahedron array (3D) ical harmonics [26]. antenna elements. Finally, the concluding remarks are made in Section VII. C. Rotation on the Sphere II. MATHEMATICAL PRELIMINARIES The rotation group SO(3) is characterized by Euler angles A. Signals on the 2-Sphere ('; #; !) 2 SO(3), where ' 2 [0; 2π), # 2 [0; π] and ! 2 [0; 2π). We define a rotation operator on the sphere We consider complex valued square integrable functions 2 D('; #; !) that rotates a function on the sphere, according defined on the 2-sphere, S . The set of such functions forms to 0zyz0 convention, in the sequence of ! rotation around the L2( 2) a Hilbert space, denoted by S , that is is equipped with z-axis, # rotation around the y-axis and ' rotation around f g the inner product defined for two functions and defined z-axis. A rotation of a function f(θ; φ) on sphere is given by on S2 as [26] −1 Z (D('; #; !)f)(xb) , f(R xb); (7) hf; gi , f(x)g(x) ds(x); (1) 2 b b b S where R is a 3 × 3 real orthogonal unitary matrix, referred which induces a norm kfk hf; fi1=2. Here, (·) as rotation matrix, that corresponds to the rotation operator , D('; #; !) and is given by denotes the complex conjugate operation and xb , [sin θ cos θ; sin θ sin φ, cos θ]T 2 2 ⊂ 3 represents a point S R R = R (')R (#)R (!); on the 2-sphere, where [·]T represents the vector transpose, z y z (8) θ 2 [0; π] and φ 2 [0; 2π) denote the co-latitude and longitude where 2cos ' − sin ' 03 respectively and ds(xb) = sin θ dθ dφ is the surface measure on the 2-sphere. The functions with finite energy (induced norm) Rz(') = 4sin ' cos ' 05 ; are referred as signals on the sphere. 0 0 1 2 cos # 0 sin #3 B. Spherical Harmonics Ry(#) = 4 0 1 0 5 : Spherical harmonics serve as orthonormal basis functions − sin # 0 cos # for the representation of functions on the sphere and are Here R (') and R (#) characterize individual rotations by defined for integer degree ` ≥ 0 and integer order jmj ≤ ` as z y ' along z-axis and # along y-axis, respectively. m m m m imφ Y` (xb) ≡ Y` (θ; φ) , N` P` (cos θ)e ; III. PROBLEM FORMULATION with s 2` + 1 (` − m)! A. Fisher-Bingham Distribution on Sphere N m ; (2) ` , 4π (` + m)! Definition 1 (Fisher-Bingham Five-Parameter (FB5) Distribu- tion).