# Efficient Time-Domain Ray-Tracing Technique for the Analysis of Ultra-Wideband Indoor Environments including Lossy Materials and Multiple Effects

Authors: F. Saez de Adana; O. Gutiérrez; M. A. Navarro; A. S. Mohan
Journal: International Journal of Antennas and Propagation
 (2009)

---

## Abstract

This paper presents an efficient application of the Time-Domain Uniform Theory of Diffraction (TD-UTD) for the analysis of Ultra-Wideband (UWB) mobile communications for indoor environments. The classical TD-UTD formulation is modified to include the contribution of lossy materials and multiple-ray interactions with the environment. The electromagnetic analysis is combined with a ray-tracing acceleration technique to treat realistic and complex environments. The validity of this method is tested with measurements performed inside the Polytechnic building of the University of Alcala and shows good performance of the model for the analysis of UWB propagation.


---

## Body

## 1. Introduction

Ultra-wideband (UWB)  technology has developed rapidly in the past several years. This technology is especially attractive in high data rate and short-range wireless communications. These applications make UWB technology suitable for indoor mobile communications applications, such as wireless personal-area networks (WPANs). This interest has motivated the study of the propagation of the UWB signals in indoor environments as an important task for the implementation of the WPANs.In the last decades, significant effort has been focused on the characterization of the indoor channel for narrowband systems. The statistical [1–4] and deterministic [5–8] models have been most frequently used in these studies. The deterministic models are mostly based on the ray-tracing techniques [5–7] to predict the multipath phenomena, and the Uniform Theory of Diffraction (UTD) technique [9] to calculate the received power or the propagation losses. However, the features of the UWB systems (with bandwidth in the range of the GHz) render the conventional narrowband propagation models, both statistical and deterministic, not applicable. These models are mostly based on frequency-domain analysis while UWB requires a time-domain analysis, due to its wide bandwidth. Therefore, special models must be used to predict the signal propagation in UWB systems. Although statistical models have been proposed [10, 11], they are not adequate, as multipath propagation produces a strong fast-fading on the indoor propagation for UWB systems which is not accurately predicted by the statistical models. Regarding the deterministic models, the frequency domain UTD can be applied, performing an analysis at several frequencies and obtaining the time response using an Inverse Fourier Transform. However, this procedure is computationally inefficient by comparison to a direct analysis in time domain.Instead, Time-Domain Uniform Theory of Diffraction (TD-UTD) was developed to obtain the solution in the time domain for the reflection and the diffraction of a transient electromagnetic wave. The inclusion of the multipath phenomena in this theory, and the analysis in the TD makes this technique suitable for UWB applications. TD-UTD was firstly developed by Vemttipong  and Kouyoumjian [12] who applied the inverse Laplace transform to the frequency-domain UTD formulation. Later, Rousseau and Pathak [13] presented closed-form solutions for the diffraction by an edge, by modifying the formulation presented in [12]. The results obtained in [13] can be directly applied to develop a method for the calculation of the indoor propagation in UWB systems.The formulation presented in [13] has been modified in this paper to introduce the contribution of lossy materials, which form the indoor environment in the reflection, transmission and diffraction. The idea is to obtain the reflection, transmission or diffraction coefficients using an Analytical Time Transform (ATT) from their expressions in the frequency domain. In addition, multiple interactions are also considered in the approach. These interactions included multiple reflections and transmissions, and the interactions between reflected and diffracted rays. Hence, both reflection-diffraction and diffraction-reflection interactions are included. These interactions, which are obtained in the frequency domain from the product of the coefficients involved in the propagation mechanisms, can also be computed in the time domain by convolving those coefficients instead.The ray-tracing associated with the computation of the transient response does not change with respect to the ray-tracing in the frequency domain. The technique developed by the authors in [5] applied to the frequency domain is also used in this approach. The main advantage of that method is the inclusion of a novel ray-tracing acceleration technique called the Angular Z-Buffer (AZB). This technique, as can be shown in [5, 14], reduces the computation time associated to the ray-tracing by several orders of magnitude. Therefore, the analysis of complex and realistic environments can now be computed much more quickly.Other TD-UTD models have also been proposed in the past. In [15], Yao proposed a TD-UTD model that included the contributions of direct, reflected and diffracted rays, in addition to the transmitted ray in [16]. A simulator based on TD-UTD for indoor UWB propagation was presented in [17]. The work in [18] took into account the polarization of both the transmitted and the received antennas by investigating a model consisting of two cube-like objects. The method proposed in this work improves on these previous approaches by generalizing the TD-UTD formulation for the analysis of realistic environments. This aim is achieved by including the multiple-ray contributions. Moreover, the electromagnetic formulation is combined with a previously-developed ray-tracing acceleration technique to reduce the computational cost for analysis of these realistic environments. Finally, with the aim of further reducing this computation time, the convolution of the analytical impulse response with the excitation is performed in an efficient way by expressing the excitation waveform as a sum of simple expansion functions thus allowing the convolution to be performed in closed form.Taking these into account, the analysis of a real site has been performed to prove the validity of the model and its ability to analyze realistic environments. Some experimental measurements have been performed and comparisons with the predictions of the proposed model are presented in this paper. The analysis of such a complex site was made possible due to the improvements included in our approach.


## 2. TD-UTD Model for Indoor Propagation

As shown in Figure1, the classical UTD in the frequency domain obtains the field at an observation point inside an indoor environment as the sum of the contribution of different rays. These rays that started from the source (S) reached that observation point (O) either directly or after one or several reflections, diffractions, transmissions or the multiple combinations of these effects. Accordingly, the TD-UTD analytical impulse response in that environment can be obtained from an ATT, which consisted of  a one-sided Inverse Fourier Transform (IFT) of the frequency response, as can be seen in [12, 13], and can be written as:Figure 1
Ray-tracing mechanisms between a given source (S) and an observer (O). The figure shows the complexity of the multipath phenomenon in an indoor environment.(1)f+I(t)=fIi+(t)+fIr+(t)+fIt+(t)+fId+(t)+fImr+(t)+fImt+(t)+fIrd+(t)+fIdr+(t),
where fIi+(t), fIr+(t), fIt+(t), fId+(t), fImr+(t), fImt+(t), fIrd+(t), and fIdr+(t) are the analytical signal representations for the direct field fIi+(t), reflected field fIr+(t), transmitted field fIt+(t), diffracted field fId+(t), multiply reflected field fImr+(t), multiply transmitted fImt+(t), reflected-diffracted-field fIrd+(t) and diffracted-reflected field fIdr+(t), respectively. Therefore, the impulse response as shown in (1) included all the multipath phenomena as mentioned previously. Each term in (1) will be described in the following sections.


### 2.1. Direct Field

The contribution of the direct field to the impulse response is obtained as the ATT of the usual Geometrical Optics (GO) incident field and can be expressed by the following equation [13]:(2)fIi+(t)=E0i|Ai(si)|δ+(t-sic)
where E0i is the initial field value, which is constant with time and frequency, Ai(si) is the spreading (or spatial divergence) factor for the direct ray, and si is the distance between the source and the observation point. The spreading factor for the direct ray is given by [13]:(3)Ai(si)=ρ1iρ2i(ρ1i+si)(ρ2i+si)
where ρ1i and ρ2i are the principal radii of curvature of the incident wavefront at the observation point (see Figure 2).Figure 2
Incident ray wavefront.


### 2.2. Reflected Field

Similar to the case of the direct field, the contribution of the reflected field to the impulse response is obtained from the ATT of the classical GO expression in the frequency domain by the following equation:(4)fIr+(t)=E0iR̅+(t-sic-src)|Ai(si)Ar(sr)|,
where si is, in this case, the distance between the source and the reflection point, sr is the distance between the reflection point and the observation point, Ar(sr) is the spreading (or spatial divergence) factor for the reflected ray which is expressed as [13](5)Ar(sr)=ρ1rρ2r(ρ1r+sr)(ρ2r+sr)
where  ρ1r and ρ2r are the principal radii of curvature of the reflected wavefront at the observation point.In (4), R̅+(t) is the TD dyadic reflection coefficient which must be expressed in terms of its vertical and parallel components. These coefficients are obtained as the ATT of the classical Fresnel reflection coefficients [15] for a reflected surface composed by a lossy material. By performing this ATT, the parallel component can be expressed as(6)R+∥(t)=Kδ+(t)+4κ1-κ2e-att∑n=1∞(-1)n+1nKnIn(at),
where In is the modified Bessel function of order n and(7)K=1-κ1+κ,κ=εr-cos2ϕεrsinϕ,a=120πσc2εr
and the perpendicular component is as follows(8)R+⊥(t)=-[Kδ+(t)+4κ1-κ2e-att∑n=1∞(-1)n+1nKnIn(at)],
where in this case,
(9)κ=sinϕεr-cos2ϕ,a=120πσc2εr(1-cos2ϕεr)
and the rest of parameters are the same.


### 2.3. Transmitted Field

The impulse response for the transmitted field is analogous to the response for the reflected field and can be written as:(10)fIr+(t)=E0iT̅+(t-sic-stc)|Ai(si)At(st)|,
where si is, in this case, the distance between the source and the transmission point, st is the distance between the transmission point and the observation point, At(st) is the spreading (or spatial divergence) factor for the transmitted ray that can be obtained using expression (5) by replacing the radii of curvature of the reflected wavefront by those of the transmitted wavefront, and T̅+(t) is the TD dyadic transmission coefficient, which must be also expressed in terms of its vertical and parallel components. The transmission coefficient can be easily obtained by considering that the relationship between the reflection and the transmission coefficients is the same in the time domain than as in in the frequency domain, and is given by(11)T̅+(t)=I+R̅+(t),
where I is the identity matrix.


### 2.4. Diffracted Field

In the case of diffraction, its contribution to the impulse response is given by the ATT of the UTD expression for the frequency domain as follows:(12)fId+(t)=E0iD̅+(t-sic-sdc)|Ai(si)Ad(sd)|,
where si is, in this case, the distance between the source and the diffraction point, sd is the distance between the diffraction point and the observation point, Ad(sd) is the spreading (or spatial divergence) factor for the reflected ray given by [13](13)Ad(sd)=ρs(ρs+sd)·sd,
where  ρs is the principal radius of curvature of the diffracted wavefront at the observation point.On the other hand,D̅+(t) is the TD dyadic diffraction coefficient, which must be expressed in terms of its components with respect to the edge-fixed system. The diffraction coefficients for a PEC wedge are obtained as the sum of four terms as  in [13]. If the lossy materials are included in the formulation, the last two terms corresponding to the contribution of the reflected shadow boundary must include the effect of the lossy reflection coefficient. In the frequency domain, this is performed by a product. Therefore, the convolution between the reflection coefficient and these two terms of the diffraction coefficient must instead be performed in the time domain. Performing this convolution results in the following expression for the diffraction coefficients:(14)D̅+β,ϕ(t)=D+1(t)+D+2(t)+R̅+⊥,∥*[D+3(t)+D+4(t)].The expressions forD+1(t), D+2(t), D+3(t) and D+4(t) are given as follows [13]:
(15)D+1(t)=A0cot[π+(ϕ-ϕ′)2n]F+[Lia+(ϕ-ϕ′)],D+2(t)=A0cot[π-(ϕ-ϕ′)2n]F+[Lia-(ϕ-ϕ′)],D+3(t)=A0cot[π+(ϕ+ϕ′)2n]F+[Lrna+(ϕ+ϕ′)],D+4(t)=A0cot[π-(ϕ+ϕ′)2n]F+[Lroa-(ϕ+ϕ′)],
with A0=-1/(2n2πsinβ0i) where β0i is the angle between the direction of incidence and the vector of the edge and F+(x,t)=x/π(jt+x/c))/t(t+x/c).TheLi are distance parameters associated with the incident shadow boundaries and are the same in the frequency domain. These parameters are given by(16)Li=sd(ρei+sd)ρ1iρ2iρei(ρ1i+sd)(ρ2i+sd)sin2β0i,(17)Lro,n=sd(ρero,n+sd)ρ1ro,nρ2ro,nρero,n(ρ1ro,n+sd)(ρ2ro,n+sd)sin2β0i,
where ρ1i is the radii of curvature 1 of the incident wavefront at the edge, ρ2i is the radii of curvature 2 of the incident wavefront at the edge, ρei is the radii of curvature of the incident wavefront at edge-fixed plane of incidence, ρ1ro,rn is the radii of curvature 1 of the reflected wavefront from the o and n faces, respectively, ρ2ro,rn is the radii of curvature 2 of the reflected wavefront from the o and n faces, respectively, and ρero,rn is the radii of curvature of the reflected wavefront from in the plane containing the reflected ray and the edge.The functiona± of the expressions (15) is given by
(18)a±(β±)=2cos2(2nπN±-β±2),β±=ϕd±ϕi,2πnN+-(ϕd±ϕi)=π,2πnN--(ϕd±ϕi)=-π.The geometrical parameters involved in the calculation of the diffraction coefficients are the same as in the frequency domain as shown in Figure3 and are explained in [9].Figure 3
Geometrical parameters and edge-fixed system for the diffracted ray. In the figureê is the vector of the edge, Qd is the diffraction point, ŝi and ŝd are the incident and diffraction directions, respectively, β0i, β0d, ϕi and ϕd (β0i=β0d=β0) are the diffraction angles and β̂0i,ϕ̂i and β̂0d,ϕ̂d are the vectors of the edge-fixed system necessary to apply the diffraction coefficients. All these parameters are explained in [9], because they are the same to the frequency domain case.


### 2.5. Multiple Reflected and Multiple Transmitted Fields

The expressions for them-order reflections and transmissions are easily derived recursively from the first order effects. For instance, the second-order reflection will be a single reflection where the source is set as the first reflection point and the incident field is the simple reflected field. Using this recursion, an m-order reflection that reaches the observation point would contribute the following term to the impulse response:(19)fImr+(t)=E0iR̅+M|At(st)|,
where(20)st=si+∑j=1msrj,At(st)=Ai(si)∏j=1mAr(srj),R̅+M=R̅+1(t-sic-sr1c)*⋯*R̅+m(t-sic-∑i=1msrjc),withsrj being the distance between the jth the (j+1)th reflection points (or between the mth reflection point and the observation point in the case of the last reflection).Analogously, them-order transmitted field will be(21)fImt+(t)=E0iT̅+M|At(st)|,
where T̅M+(t)=I+R̅M+(t).


### 2.6. Reflected-Diffracted and Diffracted-Reflected Fields

Following the same procedure for multiple reflections, the contribution of the interaction between an edge and a reflecting surface to the impulse response can be written as(22)fIrd+(t)=E0i[R̅+(t-sic-src)*D̅+(t-sic-src-sdc)]·|Ai(si)Ar(sr)Ad(sd)|
for the case of reflection-diffraction interaction and(23)fIdr+(t)=E0i[D̅+(t-sic-sdc)*R̅+(t-sic-src-sdc)]·|Ai(si)Ad(sd)Ar(sr)|
for the case of diffraction-reflection interaction.The meaning of these parameters is analogous to the previous effects.


## 2.1. Direct Field

The contribution of the direct field to the impulse response is obtained as the ATT of the usual Geometrical Optics (GO) incident field and can be expressed by the following equation [13]:(2)fIi+(t)=E0i|Ai(si)|δ+(t-sic)
where E0i is the initial field value, which is constant with time and frequency, Ai(si) is the spreading (or spatial divergence) factor for the direct ray, and si is the distance between the source and the observation point. The spreading factor for the direct ray is given by [13]:(3)Ai(si)=ρ1iρ2i(ρ1i+si)(ρ2i+si)
where ρ1i and ρ2i are the principal radii of curvature of the incident wavefront at the observation point (see Figure 2).Figure 2
Incident ray wavefront.


## 2.2. Reflected Field

Similar to the case of the direct field, the contribution of the reflected field to the impulse response is obtained from the ATT of the classical GO expression in the frequency domain by the following equation:(4)fIr+(t)=E0iR̅+(t-sic-src)|Ai(si)Ar(sr)|,
where si is, in this case, the distance between the source and the reflection point, sr is the distance between the reflection point and the observation point, Ar(sr) is the spreading (or spatial divergence) factor for the reflected ray which is expressed as [13](5)Ar(sr)=ρ1rρ2r(ρ1r+sr)(ρ2r+sr)
where  ρ1r and ρ2r are the principal radii of curvature of the reflected wavefront at the observation point.In (4), R̅+(t) is the TD dyadic reflection coefficient which must be expressed in terms of its vertical and parallel components. These coefficients are obtained as the ATT of the classical Fresnel reflection coefficients [15] for a reflected surface composed by a lossy material. By performing this ATT, the parallel component can be expressed as(6)R+∥(t)=Kδ+(t)+4κ1-κ2e-att∑n=1∞(-1)n+1nKnIn(at),
where In is the modified Bessel function of order n and(7)K=1-κ1+κ,κ=εr-cos2ϕεrsinϕ,a=120πσc2εr
and the perpendicular component is as follows(8)R+⊥(t)=-[Kδ+(t)+4κ1-κ2e-att∑n=1∞(-1)n+1nKnIn(at)],
where in this case,
(9)κ=sinϕεr-cos2ϕ,a=120πσc2εr(1-cos2ϕεr)
and the rest of parameters are the same.


## 2.3. Transmitted Field

The impulse response for the transmitted field is analogous to the response for the reflected field and can be written as:(10)fIr+(t)=E0iT̅+(t-sic-stc)|Ai(si)At(st)|,
where si is, in this case, the distance between the source and the transmission point, st is the distance between the transmission point and the observation point, At(st) is the spreading (or spatial divergence) factor for the transmitted ray that can be obtained using expression (5) by replacing the radii of curvature of the reflected wavefront by those of the transmitted wavefront, and T̅+(t) is the TD dyadic transmission coefficient, which must be also expressed in terms of its vertical and parallel components. The transmission coefficient can be easily obtained by considering that the relationship between the reflection and the transmission coefficients is the same in the time domain than as in in the frequency domain, and is given by(11)T̅+(t)=I+R̅+(t),
where I is the identity matrix.


## 2.4. Diffracted Field

In the case of diffraction, its contribution to the impulse response is given by the ATT of the UTD expression for the frequency domain as follows:(12)fId+(t)=E0iD̅+(t-sic-sdc)|Ai(si)Ad(sd)|,
where si is, in this case, the distance between the source and the diffraction point, sd is the distance between the diffraction point and the observation point, Ad(sd) is the spreading (or spatial divergence) factor for the reflected ray given by [13](13)Ad(sd)=ρs(ρs+sd)·sd,
where  ρs is the principal radius of curvature of the diffracted wavefront at the observation point.On the other hand,D̅+(t) is the TD dyadic diffraction coefficient, which must be expressed in terms of its components with respect to the edge-fixed system. The diffraction coefficients for a PEC wedge are obtained as the sum of four terms as  in [13]. If the lossy materials are included in the formulation, the last two terms corresponding to the contribution of the reflected shadow boundary must include the effect of the lossy reflection coefficient. In the frequency domain, this is performed by a product. Therefore, the convolution between the reflection coefficient and these two terms of the diffraction coefficient must instead be performed in the time domain. Performing this convolution results in the following expression for the diffraction coefficients:(14)D̅+β,ϕ(t)=D+1(t)+D+2(t)+R̅+⊥,∥*[D+3(t)+D+4(t)].The expressions forD+1(t), D+2(t), D+3(t) and D+4(t) are given as follows [13]:
(15)D+1(t)=A0cot[π+(ϕ-ϕ′)2n]F+[Lia+(ϕ-ϕ′)],D+2(t)=A0cot[π-(ϕ-ϕ′)2n]F+[Lia-(ϕ-ϕ′)],D+3(t)=A0cot[π+(ϕ+ϕ′)2n]F+[Lrna+(ϕ+ϕ′)],D+4(t)=A0cot[π-(ϕ+ϕ′)2n]F+[Lroa-(ϕ+ϕ′)],
with A0=-1/(2n2πsinβ0i) where β0i is the angle between the direction of incidence and the vector of the edge and F+(x,t)=x/π(jt+x/c))/t(t+x/c).TheLi are distance parameters associated with the incident shadow boundaries and are the same in the frequency domain. These parameters are given by(16)Li=sd(ρei+sd)ρ1iρ2iρei(ρ1i+sd)(ρ2i+sd)sin2β0i,(17)Lro,n=sd(ρero,n+sd)ρ1ro,nρ2ro,nρero,n(ρ1ro,n+sd)(ρ2ro,n+sd)sin2β0i,
where ρ1i is the radii of curvature 1 of the incident wavefront at the edge, ρ2i is the radii of curvature 2 of the incident wavefront at the edge, ρei is the radii of curvature of the incident wavefront at edge-fixed plane of incidence, ρ1ro,rn is the radii of curvature 1 of the reflected wavefront from the o and n faces, respectively, ρ2ro,rn is the radii of curvature 2 of the reflected wavefront from the o and n faces, respectively, and ρero,rn is the radii of curvature of the reflected wavefront from in the plane containing the reflected ray and the edge.The functiona± of the expressions (15) is given by
(18)a±(β±)=2cos2(2nπN±-β±2),β±=ϕd±ϕi,2πnN+-(ϕd±ϕi)=π,2πnN--(ϕd±ϕi)=-π.The geometrical parameters involved in the calculation of the diffraction coefficients are the same as in the frequency domain as shown in Figure3 and are explained in [9].Figure 3
Geometrical parameters and edge-fixed system for the diffracted ray. In the figureê is the vector of the edge, Qd is the diffraction point, ŝi and ŝd are the incident and diffraction directions, respectively, β0i, β0d, ϕi and ϕd (β0i=β0d=β0) are the diffraction angles and β̂0i,ϕ̂i and β̂0d,ϕ̂d are the vectors of the edge-fixed system necessary to apply the diffraction coefficients. All these parameters are explained in [9], because they are the same to the frequency domain case.


## 2.5. Multiple Reflected and Multiple Transmitted Fields

The expressions for them-order reflections and transmissions are easily derived recursively from the first order effects. For instance, the second-order reflection will be a single reflection where the source is set as the first reflection point and the incident field is the simple reflected field. Using this recursion, an m-order reflection that reaches the observation point would contribute the following term to the impulse response:(19)fImr+(t)=E0iR̅+M|At(st)|,
where(20)st=si+∑j=1msrj,At(st)=Ai(si)∏j=1mAr(srj),R̅+M=R̅+1(t-sic-sr1c)*⋯*R̅+m(t-sic-∑i=1msrjc),withsrj being the distance between the jth the (j+1)th reflection points (or between the mth reflection point and the observation point in the case of the last reflection).Analogously, them-order transmitted field will be(21)fImt+(t)=E0iT̅+M|At(st)|,
where T̅M+(t)=I+R̅M+(t).


## 2.6. Reflected-Diffracted and Diffracted-Reflected Fields

Following the same procedure for multiple reflections, the contribution of the interaction between an edge and a reflecting surface to the impulse response can be written as(22)fIrd+(t)=E0i[R̅+(t-sic-src)*D̅+(t-sic-src-sdc)]·|Ai(si)Ar(sr)Ad(sd)|
for the case of reflection-diffraction interaction and(23)fIdr+(t)=E0i[D̅+(t-sic-sdc)*R̅+(t-sic-src-sdc)]·|Ai(si)Ad(sd)Ar(sr)|
for the case of diffraction-reflection interaction.The meaning of these parameters is analogous to the previous effects.


## 3. The Angular Z-Buffer (AZB) Algorithm

The AZB technique was presented in [5, 14] and is summarized here due to its applicability to the time-domain formulation to reduce the computational cost associated with the analysis of realistic environments. The method consisted of dividing the space from the source point into spherical sectors called anxels. These anxels are defined by the spherical coordinates, theta (θ) and phi (ϕ) for a fixed coordinate system with the origin set at the source.The anxel belonging to each facet is computed from the spherical coordinates of its vertices. This information is stored in the so-called AZB matrix. This information depends, exclusively, on the source point and on the environmental model. The facets of each anxel are sorted according to their distances from the source (S). For instance, Figure4 shows a simple 2D outdoor scene where a space partitioning in 12 anxels has been accomplished. Table 1 shows the storage of the facets in the AZB matrix. The procedure is easily generalized for the 3D case by including the second spherical coordinate.Table 1
Facets storage in the AZB of the scene of Figure4.
AnxelFacets114, 15, 16, 18, 19, 20, 21214, 16, 17, 27, 28, 18, 21, 31, 32, 30, 33, 34314, 13, 17, 27, 29414, 13, 23, 24, 25, 26512, 11, 23, 22, 26610, 12, 9, 11710, 9891, 2101, 2111, 3, 2, 4, 5, 6124, 7, 8, 6Figure 4
Example of a 2D outdoor scene divided in 12 “anxels”.To perform the shadowing test for a given observation point (O), O is located in the corresponding anxel by means of its spherical coordinates. Then facets placed in the anxel with distances to S less than the distance S-O are included. This inclusion allows reduction of the computation time for the shadowing test calculation by several orders of magnitude. As this test, as explained in [14],  encompassed most of the computational cost of the ray-tracing techniques, the complexity of the analysis in general is drastically reduced by the application of the AZB algorithm.


## 4. Discussion of Results and Practical Considerations

The measurements in our experiments have been performed in a complex realistic site to investigate the validity of the approach. The measurements have been done in the corridor of the second floor of the Polytechnic Building of the University of Alcala. Figure5 shows the plan schematic of the measurement site. The dimensions of the scenario are 7.9×20.7 metres. A 3D plane-facets model has been designed to represent the realistic environment composed of 77 facets. The material composition of the elements of the site was concrete for the walls, wood for the doors and glass for the windows. Table 2 lists the electrical properties of these three materials considered in our model.Table 2
Electrical characteristics of the materials of the measurements scenario.
Materialεrσ (S/m)Concrete4.50.01Wood210-5Glass6.510-12Figure 5
Plant of the Polytechnic Building corridor of the Alcala University and TX and RX placements for the measurement antennas.Several measurements were performed on the site. Examples of one Line of Sight (LOS) and one Non Line of Sight (NLOS) case will be shown in this section. Figure5 illustrates the position of the transmitter and the receiver in both cases. The coordinates of the transmitter were (1.60, 5.55, 1.10) and the receivers were (5.0, 1.78, 1.10) for the LOS case and (9.55, 1.3, 1.10) for the NLOS case. All the coordinates are given in metres.Measurements have been conducted in the frequency domain using the network analyzer (VNA) Agilent E8362B. Two Double-Ridged Waveguide antennas linearly polarized were used as transmitter (TX) and receiver (RX), respectively, similar to the approach in [19]. The frequency range of these antennas was 1 to 18 GHz. In this range, the average VSWR ratio was lower than 1.5.Figure6 shows an overview of the measurements setup. The VNA was set to transmit 3201 tones uniformly distributed over the 1–18 GHz frequency range. This gave an excess delay of ~188 ns and a maximum distance of 56.4 m. The temporal resolution for the 17 GHz frequency was 59 ps.Figure 6
Overview of the measurements setup.The input signal is expressed as a sum of a small number of simple expansion functions for a more efficient convolution with the TD-UTD impulse response. In this approach the input signal is represented as the sum of waveforms whose analytical signal representations are simple poles in the complex-time plane. This representation allows the convolution to be expressed in a closed form, thus speeding up computation. Moreover, in this representation is included the antenna transfer function. Equations (24) and (25) show the representation of the input signal and the closed form for the convolution, respectively:(24)i+(t)=jπ∑n=1nAnt+jτn,(25)o+(t)=12i+(t)*f+(t)=∑n=1NAnf+(t+jτn).The comparison between the measurements and our approach for the normalized PDP in the LOS and the NLOS cases is shown in Figures7 and 8, respectively. As can be seen, a good agreement between calculation and measurement is obtained in both cases. The mean errors were only 3.5 dB for the LOS case and 4.6 dB for the NLOS case, which are very good for UWB applications. If the lossy components of the materials are not included the mean errors increase to 7.2 dB and 10.1 dB, respectively. Obviously, neglecting the lossy components in the model decreases the accuracy of the prediction. The computation time was 1 minute 23 seconds. for the LOS case and 1 minute 42 seconds. for the NLOS case considering up to sixth order reflections. Both results were obtained using an AMD Athlon microprocessor at 1.99 GHz. Hence, the method is very suitable for analysing realistic environments at high speed. This high computation speed was due to the combination of the TD-UTD formulation with the AZB ray-tracing acceleration technique. It is important to mention that the computation time includes the creation of the AZB matrices. Therefore, considering that the AZB matrices only have to be created once, the computation time for the NLOS case is only 35 seconds after the computation of the LOS case. Therefore, the calculation of additional point would only take an amount of time of that order of magnitude.Figure 7
Comparison between measurements and simulation for the LOS case.Figure 8
Comparison between measurements and simulation for the LOS case.


## 5. Conclusions

An accurate, high-speed method based on TD-UTD to analyze indoor propagation for UWB systems has been presented. The time domain analysis used is especially suitable for the problem under study, but it is also important to include in the TD-UTD formulation the contribution of lossy materials and multiple effects for this kind of analysis. Experimental measurements have been performed at the University of Alcala and the agreement of the presented model with the measured results is excellent. The main advantage of the approach is the possibility of analyzing UWB propagation in realistic environments in a reasonable amount of time. This advantage is afforded by the combination of the electromagnetic analysis with acceleration ray-tracing techniques such as the AZB algorithm.



---

*Source: 390782-2010-02-22.xml*