Chapter 6 GEOMETRIC MODELS FOR ANGLE AND TIME OF ARRIVAL 1. INTRODUCTION Array antennas are sometimes employed at the receiver for a variety of purposes such as to combat fading or to reduce co-channel interference.To evaluate the performance of a wireless communication system using antenna arrays it becomes necessary to have spatial channel models that describe the angle of arrival and time of arrival of the multipath components.Among the most widely used radio propagation models is the single bounce scattering geometric model,where propagation between the transmitting and receiving antennas is assumed to take place via single scattering from an intervening ob- stacle.The scattering cross section of the scatterers is assumed to be isotropic, but their density is umed to vary fro om locatio n to locat Currently,sev eral geometric models are available such as the ring model 53 discrete uni form model [118],circular scattering model [30],elliptical scattering model [981.Gaussian angle of arrival model [931.etc.with each model being ap- plicable to a s ecific environment type For example,the circular scatter ble to macro cell type of environments,whereas the elliptic model is suitable for a micro-or a pico-cell types of environment.Most of these models are 2D in nature in that they assume radio propagation to take place in a plane containing the transmitter and the receiver.In this Chapter we will consider the circular scattering model,the elliptical scattering model and the Gaussian density model and show statistics for the angle and time of arrival
Chapter 6 GEOMETRIC MODELS FOR ANGLE AND TIME OF ARRIVAL 1. INTRODUCTION Array antennas are sometimes employed at the receiver for a variety of purposes such as to combat fading or to reduce co-channel interference. To evaluate the performance of a wireless communication system using antenna arrays it becomes necessary to have spatial channel models that describe the angle of arrival and time of arrival of the multipath components. Among the most widely used radio propagation models is the single bounce scattering geometric model, where propagation between the transmitting and receiving antennas is assumed to take place via single scattering from an intervening obstacle. The scattering cross section of the scatterers is assumed to be isotropic, but their density is assumed to vary from location to location. Currently, several geometric models are available such as the ring model [53], discrete uniform model [118], circular scattering model [30], elliptical scattering model [98], Gaussian angle of arrival model [93], etc., with each model being applicable to a specific environment type. For example, the circular scattering model is applicable to macro cell type of environments, whereas the elliptical model is suitable for a micro- or a pico-cell types of environment. Most of these models are 2D in nature in that they assume radio propagation to take place in a plane containing the transmitter and the receiver. In this Chapter we will consider the circular scattering model, the elliptical scattering model, and the Gaussian density model and show statistics for the angle and time of arrival
164 RADIO WAVE PROPAGATION SMART ANTENNAS 2. GENERAL FORMULATION Figure 6.1 shows the geometry of scatterers distributed around the mobile station which is transmitting to a base station.The scatter region is considered convex in shape and can possibly extend to infinity.It is assumed that propa- gation takes place in the horizontal plane containing the tip of the transmitting and rec ving antennas. The global in this plane r edenoted b (x,y).The line joining the transmitter and receiver makes an angle with the x-axis and the coordinates of a point with respect to this line are (, The polar coordinates with respect to the base station are(r,B).The mobile station is assumed to be located at a distance D from the b estation along this line. The coordin ates of a poin with respect to the mobile station are denoted by (m,ym).The polar coordinates with respect to the mobile station are (rm,m).The relationships between these various coordinates are 6=TocosB=xcos+ysin;Im=rm cos om=-D (6.la) yo ro sin B=-xsin+ycos; hm=Tmsin中m=b (6.1b) r2=r话+D2-2rDcos 中=6+. (6.1c) The density of scatterers about the mobile station is given as a probabil- m Scatter Region y%↑ X BS D Figure 6.1.Geometrically based single bounce scattering model. ity density function pi(m,ym).Waves from the transmitter undergo single scattering from the scatter region before arriving at the receiver.The delay variable is related to the total path length traveled+m and the speed of light c via CT To+rm (6.2a) =r%+Vr号+D2-2 D cos B. (6.2b)
164 RADIO WAVE PROPAGATION & SMART ANTENNAS 2. GENERAL FORMULATION Figure 6.1 shows the geometry of scatterers distributed around the mobile station which is transmitting to a base station. The scatter region is considered convex in shape and can possibly extend to infinity. It is assumed that propagation takes place in the horizontal plane containing the tip of the transmitting and receiving antennas. The global coordinates in this plane are denoted by (x, y). The line joining the transmitter and receiver makes an angle with the x-axis and the coordinates of a point with respect to this line are The polar coordinates with respect to the base station are The mobile station is assumed to be located at a distance D from the base station along this line. The coordinates of a point with respect to the mobile station are denoted by The polar coordinates with respect to the mobile station are The relationships between these various coordinates are The density of scatterers about the mobile station is given as a probabilFigure 6.1. Geometrically based single bounce scattering model. ity density function Waves from the transmitter undergo single scattering from the scatter region before arriving at the receiver. The delay variable is related to the total path length traveled and the speed of light c via
Geometric Models for Angle and Time of Arrival 165 The locus of all points having a constant value of is an ellipse with the base station and mobile station at its foci as shown in Figure 6.2.Waves arriving Ip +rm=ct Scatter Region A Figure 6.2.Waves arriving with delay less thanr will arrive from the intersection of the ellipse and the scatter region. from the scattering region that is enclosed by this ellipse will have delays less than Maximum delay TM occurs when the ellipse just encloses all points of the scatter region.If the scattering region extends to infinity,the maximun delay becomes infinity. The objective is to determine the marginal p.d.f.p(B)of the angle of arrival (AOA)and the marginal p.d.f.p()of the time of arrival (TOA)as seen from the given is assumed that the radiatio patterns of the transmitting and receiving antennas are omnidirectional in the xy plane.This is not a serious limitation as the transmitting antenna pattern could always be included n the definition of p(m,ym)and the receiving p().The ttem could always of random variables and p e problem from t lan sys tem (xm,ym)to the polar coordinate system (ro,B)or to (T,B).The joint density in (B)can be obtained from the standard Jacobian of the transfor- mation between Cartesian and polar coordinates relating the elemental area damdym==rodrodB.The result is p(Tb,B)=rop1(xm,ym) cing-D (6.3) A.PDFfor AOA. The marginal p.d.f.p(B)is obtained by including contribution from all points along the line B=constant.If this line intersects the scatter region at rbm()and roM()with0≤rm≤r%≤rbM<oo,then TAM(B) p()= p(r6,B)dra rom(B)
Geometric Models for Angle and Time of Arrival 165 The locus of all points having a constant value of is an ellipse with the base station and mobile station at its foci as shown in Figure 6.2. Waves arriving Figure 6.2. Waves arriving with delay less than will arrive from the intersection of the ellipse and the scatter region. from the scattering region that is enclosed by this ellipse will have delays less than Maximum delay occurs when the ellipse just encloses all points of the scatter region. If the scattering region extends to infinity, the maximum delay becomes infinity. The objective is to determine the marginal p.d.f. of the angle of arrival (AOA) and the marginal p.d.f. of the time of arrival (TOA) as seen from the base station given It is assumed that the radiation patterns of the transmitting and receiving antennas are omnidirectional in the xy plane. This is not a serious limitation as the transmitting antenna pattern could always be included in the definition of and the receiving antenna pattern could always be included in and The problem boils down to transformation of random variables from the Cartesian system to the polar coordinate system or to The joint density in can be obtained from the standard Jacobian of the transformation between Cartesian and polar coordinates relating the elemental area The result is A. PDF for AOA: The marginal p.d.f. is obtained by including contribution from all points along the line = constant. If this line intersects the scatter region at and with then
166 RADIOWAVE PROPAGATION SMART ANTENNAS ToM(B) ropi(ro CosB-D,rosinB)dro (6.4) rom(B) It is of interest to detemine the mean angand the variancefor a given distribution.It is reasonable to assume that the scatter density is symmetrically distributed about the line joining the base station and the mobile.In this case the p.d.f.p()will be an even function in B and the mean angle =0.The 2=B2p()d (6.5) Equation (6.5)can be cast in a different form by making use of Parseval's theorem for Fourier series.Expressing Band p()as a Fourier series in the interval (,)equation (6.5)can be expressed as 3+m∑(-1)马 (6.6) where the even p.d.f.p()is expressed as Fourier cosine series 间=云+B,on0 (6.7) n=1 with coefficients Bn 2 p(B) (6.8) Equation (6.6)is very convenient if the Fourier cosine coefficients of p(B) can be found in a closed form.The r.m.s.angular spread in a scattering environment is defined as 2(-B). In the study of array antennas,it is also of interest to determine the complex spatial correlation p(kod)=(e-jkodco (6.9) between the field received at two points located along the x-axis and sepa- rated by a distance d.The quantity ko is the usual wavenumber in free-space and angle brackets denote expectation with respect to the angular variable =+B.In Appendix B.it is shown that the expression for (kod)for
166 RADIOWAVE PROPAGATION & SMART ANTENNAS It is of interest to determine the mean angle and the variance for a given distribution. It is reasonable to assume that the scatter density is symmetrically distributed about the line joining the base station and the mobile. In this case the p.d.f. will be an even function in and the mean angle The variance for a symmetric distribution is Equation (6.5) can be cast in a different form by making use of Parseval’s theorem for Fourier series. Expressing and as a Fourier series in the interval equation (6.5) can be expressed as where the even p.d.f. is expressed as Fourier cosine series with coefficients Equation (6.6) is very convenient if the Fourier cosine coefficients of can be found in a closed form. The r.m.s. angular spread in a scattering environment is defined as In the study of array antennas, it is also of interest to determine the complex spatial correlation between the field received at two points located along the x–axis and separated by a distance d. The quantity is the usual wavenumber in free-space and angle brackets denote expectation with respect to the angular variable In Appendix B, it is shown that the expression for for
Geometric Models for Angle and Time of Arrival 167 symmetric distributions,with p(-)=p(),can also be put in the form: p(kod)= dco(dsin sinB)p(B)d (6.10) Jo(kod)+2>(-j)"Jn(kod)Bn cosno, (6.11) n-1 where Jn()is the Bessel function of the first kind of order n. B.PDF for TOA: The marginal p.d.f.for the tune of arrival could,in principle,be obtained by first finding the joint density p(,B)and integrating over all possible values of B(r).However,a simpler approach is to first find the c.d.f.P(r)and then to take its derivative w.r.t.r to get p(r)[301.To determine the c.d.f.for a to this.Denoting the intersection of the scatter region with the ellipse as Ar. the c.d.f and the p.d.f.for the time of arrival are P(r)= (6.12) a p(r)= aT JJA. p(cm,hm)dzmdym (6.13) The mean delay and the second moment of the delay 72 can be determined from 干=rp()d (6.14) 不=2pr)dn (6.15) 3.ELLIPTICAL SCATTERING MODEL In this model,first described in [78],the scatter density is assumed to be estation and the mobile. he ellipse SCTM/2 the semi-minor axis6is equal to The area of the ellipse is Ae=mab.The equation of the ellipse in polar coordinates is 4a2-D2 Tbe=4a-2D Cos B (6.16)
Geometric Models for Angle and Time of Arrival 167 symmetric distributions, with can also be put in the form: where is the Bessel function of the first kind of order n. B. PDF for TOA: The marginal p.d.f. for the tune of arrival could, in principle, be obtained by first finding the joint density and integrating over all possible values of However, a simpler approach is to first find the c.d.f. and then to take its derivative w.r.t. to get [30]. To determine the c.d.f. for a particular it is essential to consider all delays that are less than All points of the scatter region that lie within the ellipse will contribute to this. Denoting the intersection of the scatter region with the ellipse as the c.d.f and the p.d.f. for the time of arrival are The mean delay and the second moment of the delay can be determined from 3. ELLIPTICAL SCATTERING MODEL In this model, first described in [78], the scatter density is assumed to be uniform inside an elliptical region of maximum delay with foci at the base station and the mobile. The semi-major axis a of the ellipse is and the semi-minor axis 6 is equal to The area of the ellipse is The equation of the ellipse in polar coordinates is