Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com
|
customercare@cbspd.com
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.
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)d≤11m−56+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=(k−1)Δτ. 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=G2G1 is the power ratio and
F(ε)=1−e−(N−bins1)Δτε1−e−Δτε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)
This model characterizes the local PDP by the pairs{Gk,τk}, as mentioned above, with τk=(k−1)Δτ 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
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
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)|δ+(t−si/c)E16
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 (Rousseau & Pathak, 1975):
Ai(si)=ρ1iρ2i(ρ1i+si)(ρ2i+si)E17
where ρ1i and ρ2i are 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¯+(t−si/c−sr/c)|Ai(si)Ar(sr)|E18
where si is, in this case, the distance between the source and the reflection point and 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 (Rousseau & Pathak, 1975)
Ar(sr)=ρ1rρ2r(ρ1r+sr)(ρ2r+sr)E19
where ρ1r and ρ2r are 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
The impulse response for the transmitted field is analogous to that response for the reflected field and can be written as
fIr+(t)=E0iT¯+(t−si/c−st/c)|Ai(si)At(st)|E24
where si is, 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¯+(t−si/c−sd/c)|Ai(si)Ad(sd)|E26
where si is, 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 ρ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 (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:
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 β0i is 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
where ρ1i is 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,rn is 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
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:
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^iand s^d are the incident and diffraction directions, respectively, β0i, β0d, ϕiand ϕ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 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
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)
Concrete
4.5
0.01
Wood
2
10-5
Glass
6.5
10-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.
This work has been supported, in part, by the University of Alcalá, project UAH GC2010-004.
References
1.CassioliD.WinM.MolischA.2001A Statistical Model for the UWB Indoor Channel, Proceedings of 20015 53rd Vehicular Technology Conference,11590-78036-728-6Greece,
2.CassioliD.WinM.MolischA.2002The Ultra-Wide Band Indoor Channel: from Statistical Model to Simulations, IEEE Journal on Selected Areas in Telecommunications,206124712570733-8716
3.DabinJ. A.HaimovichA. M.GrebelH.2006A Statistical Ultra-Wide Band Indoor Channel Model and the Effects of Antenna Directivity on Path Loss and Multipath Propagation, IEEE Journal on Selected Areas in Telecommunications2447527580733-8716
4.HashemiH.1979Simulation of the Urban Radio Propagation Channel. IEEE Transactions on Vehicular Technology,2832132250018-9545
5.KouyoumjianR. G.PathakP. H.1974A Uniform Geometrical Theory of Diffraction for an Edge in a Perfectly Conducting Surface, Proceedings of the IEEE,6211144814610018-9219
6.Lauer, A.; Bahr, A. & Wolff, I. (1984). FDTD Simulations of Indoor Propagation, Proceedings of 1994 44th Vehicular Technology Conference, Stockholm, Sweden0-78031-927-3883 EOF
7.MolischA.CassioliD.ChongC.EmamiS.FortA.KannanB.KaredalJ.KunischJ.SchantzH. G.SiwiakK.WinM. Z.2006A Comprenhensive Standardized Model for Ultrawideband Propagation Channels, IEEE Transactions on Antennas and Propagation,5411315131660001-8926X.
8.MotleyA. J.KeenanJ. M. P.1990Radio Coverage in Buildings, British Telecom Technological Journal, Special Issue Mobile Communications,8119241358-3948
9.RappaportT. S.SiedelS. Y.TakamizawaK.1991Statistical Channel Impulse Response Models for Factory and Open Plan Building Radio Communication System Design. IEEE Transactions on Communications,3957948070090-6778
10.RousseauP. R.PathakP. H.1975Time-Domain Uniform Theory of Diffraction for a Curved Wedge, IEEE Transactions on Antennas and Propagation,4312137513820001-8926X.
11.Saez de Adana, F.; Gutiérrez, O.; González, I.; Pérez, J. & Cátedra, M.F. (2000). Propagation Model Based on Ray-Tracing for the Design of Personal Communications, IEEE Transactions on Vehicular Technology, 0018-954564921052112
12.SalehA. A.ValenzuelaR. A.1987A Statistical Model for Indoor Multipath Propagation, IEEE Journal on Selected Areas on Communication,521281370733-8716
13.SeidelS. Y.RappaportT. S.1992MHz Path Loss Prediction Models for Indoor Wireless Communications in Multifloored Buildings, IEEE Transactions on Antennas and Propagation,4022072170001-8926X.
14.TarngJ. H.ChangW. R.HsuB. J.1997Three-Dimensional Modelling of 900 MHz and 2.44 GHz Radio Propagation in Corridors, IEEE Transactions on Vehicular Technology,4625195260018-9545
15.Tornevik, C.; Berg, J.E., Lotse, F. & Madfors, M. (1993). Propagation Models, Cell Planning and Channel Allocation for Indoor Applications of Cellular Systems, Proceedings of 1993 43rd Vehicular Technology Conference, pp. 867-870, ISBN 0-7803-1267-8, Secaucus, New Jersey, USA0-78031-267-8867870
16.Veruttipong, T. & Kouyoumjian, R.G. (1979). Early-Time Responses of Currents and Charges on Wedges and Strips, Proceedings of the 1979 Antennas and Propagation Society International Symposium, pp. 590-593, Seattle, Washington, USA590593
17.WhitmanG. M.KimK. S.NiverE.1995A Theoretical Model for Radio Signal Attenuation inside Buildings, IEEE Transactions on Vehicular Technology,4436216290018-9545
18.Yao, R.; Chen, Z. & Zhu, W. (2003). An Efficient Time-Domain Ray Model for UWB Indoor Multipath Propagation Channel, Proceedings of 2003 58th Vehicular Technology Conference, pp. 1293, ISBN 0-7803-7954-3, Orlando, Florida, USA0-78037-954-31293
Written By
Francisco Saez de Adana
Submitted: 07 November 2010Published: 09 August 2011