Open access peer-reviewed chapter

Propagation Models for the Characterization of the Indoor UWB Channel

By Francisco Saez de Adana

Submitted: November 7th 2010Reviewed: March 23rd 2011Published: August 9th 2011

DOI: 10.5772/20171

Downloaded: 1759

1. Introduction

Ultra-wideband (UWB) technology has developed rapidly over 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 communication applications, such as wireless personal-area networks (WPAN). This interest has motivated the study of the propagation of the UWB signals in indoor environments as an important task for the implementation of WPANs.

In the last decades, significant effort has been focused on the characterization of the indoor channel for narrowband systems. Statistical (Motley & Keenan, 1990; Saleh & Valenzuela, 1987; Seidel & Rappaport, 1992; Tornevik et al., 1993) and deterministic (Lauer et al., 1984; Saez de Adana et al., 2000; Tarng et al., 1997; Whitman et al., 1995) models have been used most frequently in these studies. The statistical models are based on the obtention of closed formulas to characterize the propagation channel. These formulas are derived from the data obtained from measurement campaigns in different environments. Alternatively, the deterministic models are based mostly on the use of ray-tracing techniques (Saez de Adana et al., 2000; Tarng et al., 1997; Whitman et al., 1995) to predict the multipath phenomena and the Uniform Theory of Diffraction (UTD) technique (Kouyoumjian & Pathak, 1974) to calculate the received power or the propagation losses. However, the features of the UWB systems (with bandwidth in the range of GHz) render the conventional narrowband propagation models, both statistical and deterministic, inapplicable. These models are based primarily 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 the number of statistical models developed for UWB systems is not as extensive as that for narrowband systems, some recent examples can be found in the literature (Cassioli et al., 2001, 2002, Dabin et al., 2006; Molisch et al., 2006). These models have been obtained, as in the case of narrowband systems, by obtaining closed expressions that fit the behavior of the received signal measured in several locations in a measurement campaign.

Regarding the deterministic models, 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 in comparison to direct analysis in the time domain. Instead, the Time-Domain Uniform Theory of Diffraction (TD-UTD) was developed to obtain a 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 first developed by Veruttipong and Kouyoumjian (Veruttipong & Kouyoumjian, 1979), who applied the inverse Laplace transform to the frequency-domain UTD formulation. Later, Rousseau and Pathak (Rousseau & Pathak, 1975) presented closed-form solutions for the diffraction by an edge by modifying the formulation presented in (Veruttipong & Kouyoumjian, 1979). The results obtained in (Rousseau & Pathak, 1975) can be directly applied to develop a method for the calculation of the indoor propagation in UWB systems.

In this chapter, the formulation of both a statistical approach and a deterministic approach are presented. The statistical approach has been selected from the available literature because it seems very suitable for the case of UWB systems. The deterministic formulation was developed by the author of this chapter and consists of modifying the formulation presented in (Rousseau & Pathak, 1975) to introduce the contribution of lossy materials present in the indoor environment to reflection, transmission and diffraction. The goal 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 this approach. These interactions include multiple reflections and transmissions and the interactions between reflected and diffracted rays. Thus, 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 deterministic approach is completed, including the analysis of a real site, which proves the validity of the model and its ability to analyze realistic environments. Some experimental measurements and comparisons with the predictions of the proposed model are presented in this chapter.

2. Statistical approach

The statistical models are obtained based on the statistical analysis of the experimental data. Some measurements are performed for a given scenario and a propagation model is obtained based on these results after the statistical analysis. In this chapter, one of these statistical models is presented. This model was obtained from measurements performed in a typical modern office building (Cassioli et al., 2002).

This propagation model provides all the parameters and distributions necessaries to obtain the power delay profile. First, the attenuation of the received power satisfies the following expression, dependent on the distance:

L(dB)={20.4log(dd0)d11m56+74log(dd0)d11mE1

The small-scale average-power delay profile (SSA-PDP) accomplishes the following expression:

g(τ)=k=1NbinsGkδ(ττk)E2

where g(τ)is the average energy received with delay , normalized to the total energy at one meter, and Nbins is the number of bins selected in the observation window. The idea of quantizing the delay axis into bins comes from (Rappaport et al., 1991). In that way, the local PDP is obtained by integrating the received power within each bin. This local PDP is expressed in terms of the pairs{Gk,τk}, where Gk is the ratio between the energy received at a distance d and the total energy received at one meter. This value is called the energy gain and is obtained over a bin with size Δτbeginning with a delayτk=(k1)Δτ. In the model proposed in (Cassioli et al., 2002), Δτ=2 ns.

An exponential decay from the second bin can be assumed. Therefore,

g(τ)=G1δ(ττ1)+k=2NbinsG2δ(ττk)e(τkτ2)εE3

where is the decay constant of the SSA-PDP.

The total average energy received in the observation interval T is

Gtot=0Tg(τ)dτ=G1+k=2NbinsG2e(τkτ2)εE4

Because the second term of expression (4) is a geometric series,

Gtot=G1[1+rF(ε)]E5

where r=G2G1is the power ratio and

F(ε)=1e(Nbins1)Δτε1eΔτεE6

A lognormal distribution around the mean value of the path loss can be considered:

Gtot=LN(L;4.3)E7

The average energy gains are obtained by inverting (5)

Gk={Gtot1+rF(ε)fork=1Gtot1+rF(ε)re(τkτ2)εfork=2,....,NbinsE8

and therefore

g(τ)=Gtot1+rF(ε){δ(ττ1)+k=2Nbins[re(τkτ2)ε]δ(ττk)}E9

This model characterizes the local PDP by the pairs{Gk,τk}, as mentioned above, with τk=(k1)Δτand Gk generated by a superposition of large- and small-scale statistics.

The first step is to obtain Gtot from (12). Next, the power ratio r and the decay factor are generated. These two values are treated as stochastic variables, modeled as lognormal variables according to the experience for narrowband systems shown in (Hashemi, 1979). The shape of the distribution is obtained from the measured data. Therefore,

ε=LN(16.1;1.27)E10
ε=LN(16.1;1.27)E11

and the width of the observation window is T=5.

With these parameters, SSA-PDP is completely specified by (9) and the local PDP can be obtained by obtaining the normalized energy gains Gk(i)of every bin k and every location i. Gk(i)are considered to be gamma-distributed independent variables with an average given by (8) and values mk generated as independent truncated Gaussian variables

mk=TN(μm(τk);σm2(τk))E12

with

μm(τk)=3.5τk73E13
σm2(τk)=1.84τk160E14

3. Deterministic approach

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 reached that observation point either directly or after one or several reflections, diffractions, transmissions or combinations of these effects. Accordingly, the TD-UTD analytical impulse response in that environment can be obtained from an ATT, which consists of a one-sided Inverse Fourier Transform (IFT) of the frequency response, as can be seen in (Rousseau & Pathak, 1975; Veruttipong & Kouyoumjian, 1979), and can be written as

f+I(t)=fIi+(t)+fIr+(t)+fIt+(t)+fId+(t)+fImr+(t)+fImt+(t)+fIrd+(t)+fIdr+(t)E15

wherefIi+(t), fIr+(t), fIt+(t), fId+(t), fImr+(t), fImt+(t), fIrd+(t)and fIdr+(t)are the analytical signal representations for the direct fieldfIi+(t), reflected fieldfIr+(t), transmitted fieldfIt+(t), diffracted fieldfId+(t), multiply reflected fieldfImr+(t), multiply transmittedfImt+(t), reflected-diffracted-field fIrd+(t)and diffracted-reflected fieldfIdr+(t), respectively. Therefore, the

impulse response shown in equation (15) includes all multipath phenomena, as mentioned previously. Each term in equation (15) will be described in the following sections.

3.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 (Rousseau & Pathak, 1975):

fIi+(t)=E0i|Ai(si)|δ+(tsi/c)E16

where E0iis 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 siis the distance between the source and the observation point. The spreading factor for the direct ray is given by (Rousseau & Pathak, 1975):

Ai(si)=ρ1iρ2i(ρ1i+si)(ρ2i+si)E17

where ρ1iand ρ2iare the principal radii of curvature of the incident wavefront at the observation point (see Figure 1).

Figure 1.

Incident-ray wavefront

3.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:

fIr+(t)=E0iR¯+(tsi/csr/c)|Ai(si)Ar(sr)|E18

where siis, in this case, the distance between the source and the reflection point and sris 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 (Rousseau & Pathak, 1975)

Ar(sr)=ρ1rρ2r(ρ1r+sr)(ρ2r+sr)E19

where ρ1rand ρ2rare the principal radii of curvature of the reflected wavefront at the observation point.

In equation (18), 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 (Yao et al., 2003) for a reflected surface composed of a lossy material. By performing this ATT, the parallel component can be expressed as

R+||(t)=Kδ+(t)+4κ1κ2eattn=1(1)n+1nKnIn(at)E20

where In is the modified Bessel function of order n and

K=1κ1+κκ=εrcos2ϕεrsinϕa=120πσc2εrE21

The perpendicular component is as follows:

R+(t)=[Kδ+(t)+4κ1κ2eattn=1(1)n+1nKnIn(at)]E22

where, in this case,

κ=sinϕεrcos2ϕa=120πσc2εr(1cos2ϕεr)E23

and the rest of parameters are the same.

3.3. Transmitted field

The impulse response for the transmitted field is analogous to that response for the reflected field and can be written as

fIr+(t)=E0iT¯+(tsi/cst/c)|Ai(si)At(st)|E24

where siis, in this case, the distance between the source and the transmission point, stis the distance between the transmission point and the observation point, At(st)is the spreading (or spatial divergence) factor for the transmitted ray, which can be obtained using expression (19) 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 also be expressed in terms of its vertical and parallel components. The transmission coefficient can be obtained easily by considering that the relationship between the reflection and the transmission coefficients are the same in the time domain as in the frequency domain and is given by

T¯+(t)=I+R¯+(t)E25

where I is the identity matrix.

3.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:

fId+(t)=E0iD¯+(tsi/csd/c)|Ai(si)Ad(sd)|E26

where siis, in this case, the distance between the source and the diffraction point, sdis 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 (Rousseau & Pathak, 1975)

Ad(sd)=ρs(ρs+sd)sdE27

where ρsis 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 (Rousseau & Pathak, 1975). 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 inclusion 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:

D¯+β,φ(t)=D+1(t)+D+2(t)+R¯+,||*[D+3(t)+D+4(t)]E28

The expressions forD+1(t), D+2(t), D+3(t)and D+4(t)are given as follows (Rousseau & Pathak, 1975):

D+1(t)=A0cot[π+(ϕϕ')2n]F+[Lia+(ϕϕ')]E29
D+2(t)=A0cot[π(ϕϕ')2n]F+[Lia(ϕϕ')]E30
D+3(t)=A0cot[π+(ϕ+ϕ')2n]F+[Lrna+(ϕ+ϕ')]E31
D+4(t)=A0cot[π(ϕ+ϕ')2n]F+[Lroa(ϕ+ϕ')]E32

withA0=12n2πsinβ0i, where β0iis the angle between the direction of incidence and the vector of the edge andF+(x,t)=x/π(jt+x/c))t(t+x/c).

The Li are distance parameters associated with the incident shadow boundaries and are the same as in the frequency domain. These parameters are given by

Li=sd(ρei+sd)ρ1iρ2iρei(ρ1i+sd)(ρ2i+sd)sin2β0iE33
Lro,n=sd(ρero,n+sd)ρ1ro,nρ2ro,nρero,n(ρ1ro,n+sd)(ρ2ro,n+sd)sin2β0iE34

where ρ1iis the radius of curvature 1 of the incident wavefront at the edge, ρ2iis the radius of curvature 2 of the incident wavefront at the edge, ρeiis the radius of curvature of the incident wavefront at edge-fixed plane of incidence, ρ1ro,rnis the radius of curvature 1 of the reflected wavefront from the o and n faces, respectively, ρ2ro,rnis the radius of curvature 2 of the reflected wavefront from the o and n faces, respectively, and ρero,rnis the radius of curvature of the reflected wavefront from the plane containing the reflected ray and the edge.

The function a in the expressions (27)-(30) is given by

a±(β±)=2cos2(2nπN±-β±2)β±=ϕd±ϕi2πnN+-(ϕd±ϕi)=π2πnN--(ϕd±ϕi)=-πE35

The geometrical parameters involved in the calculation of the diffraction coefficients are the same as in the frequency domain as shown in Figure 2 and are explained in (Kouyoumjian & Pathak, 1974).

3.5. Multiple reflected and multiple transmitted fields

The expressions for the m-order reflections and transmissions are easily derived recursively from the first-order effects. For instance, the second-order reflection is 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:

fImr+(t)=E0iR¯+M|At(st)|E36

where

st=si+j=1msrjE37
At(st)=Ai(si)j=1mAr(srj)E38
R¯+M=R¯+1(tsi/csr1/c)*.....*R¯+m(tsi/ci=1msrj/c)E39

srjis 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).

Figure 2.

Geometrical parameters and an edge-fixed system for the diffracted ray. In the figure, e ^ is the vector of the edge, Qd is the diffraction point, s ^ i and s ^ d are the incident and diffraction directions, respectively, β 0 i , β 0 d , ϕ i and ϕ d ( β 0 i = β 0 d = β 0 ) are the diffraction angles, and β ^ 0 i , ϕ ^ i and β ^ 0 d , ϕ ^ d are the vectors of the edge-fixed system necessary to apply the diffraction coefficients. All parameters are explained in (Kouyoumjian & Pathak, 1974) because they are identical to those in the frequency domain case

Analogously, the m-order transmitted field is

fImt+(t)=E0iT¯+M|At(st)|E40

where.

T¯M+(t)=I+R¯M+(t)E41

3.6. Reflected-diffracted and diffracted-reflected fields

Following the same procedure as for multiple reflections, the contribution of the interaction between an edge and a reflecting surface to the impulse response can be written as

fIrd+(t)=E0i[R¯+(tsi/csr/c)*D¯+(tsi/csr/csd/c)]·|Ai(si)Ar(sr)Ad(sd)|E42

for the case of reflection-diffraction interactions and

fIdr+(t)=E0i[D¯+(tsi/csd/c)*R¯+(tsi/csr/csd/c)]·|Ai(si)Ad(sd)Ar(sr)|E43

for the case of diffraction-reflection interaction.

The meaning of these parameters is analogous to the previous effects.

4. Deterministic analysis of a realistic environment

The results of the deterministic approach presented in this chapter compared with measurements are shown in this section. The measurements were performed in a complex realistic site to investigate the validity of the approach. The measurements were done in the corridor of the second floor of the Polytechnic Building of the University of Alcala. Figure 3 shows the plan schematic of the measurement site. The dimensions of the scenario are 7.9x20.7 m. 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 I lists the electrical properties of these three materials considered in our model.

Materialεrσ (S/m)
Concrete4.50.01
Wood210-5
Glass6.510-12

Table 1.

Electrical characteristics of the materials for the measurement scenario

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. Figure 3 illustrates the position of the transmitter and the receiver in both cases. The coordinates of the transmitter were (1.60, 5.55, 1.10). The coordinates of 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 meters.

Measurements were conducted in the frequency domain using the network analyzer (VNA) Agilent E8362B. A linearly polarized double-ridged waveguide antenna was used as the

Figure 3.

Plan of the Polytechnic Building corridor of Alcala University and TX and RX placements for the measurement antennas

transmitter (TX) and as the receiver (RX). The frequency range of these antennas was 1 to 18 GHz. In this range, the average VSWR ratio was lower than 1.5.

Figure 4 shows an overview of the measurement setup. The VNA was set to transmit 3201 tones uniformly distributed over the 1-18 GHz frequency range. This transmission 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 4.

Overview of the measurement 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 with analytical signal representations that 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, the antenna transfer function is included in this representation. Equations (38) and (39) show the representation of the input signal and the closed form for the convolution, respectively.

i+(t)=jπn=1nAnt+jτnE44
o+(t)=12i+(t)f+(t)=n=1NAnf+(t+jτn)E45

The comparisons between the measurements and our approach for the normalized PDP in the LOS and the NLOS cases are shown in Figures 5 and 6, respectively. As can be seen, a good agreement between calculation and measurement is obtained in both cases. The mean errors were 3.5 dB for the LOS case and 4.6 dB for the NLOS case, which are very good for UWB applications.

Figure 5.

Comparison between measurements and simulation for the LOS case

Figure 6.

Comparison between measurements and simulation for the NLOS case.

Acknowledgments

This work has been supported, in part, by the University of Alcalá, project UAH GC2010-004.

© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Francisco Saez de Adana (August 9th 2011). Propagation Models for the Characterization of the Indoor UWB Channel, Ultra Wideband Communications: Novel Trends - Antennas and Propagation, Mohammad Matin, IntechOpen, DOI: 10.5772/20171. Available from:

chapter statistics

1759total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Frequency UWB Channel

By Gonzalo Llano, Juan C. Cuellar and Andres Navarro

Related Book

First chapter

Overview of Wireless Sensor Network

By M.A. Matin and M.M. Islam

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us