Characterisation of the Satellite-to-Indoor Propagation Channel

Please download to get full document.

View again

of 6
All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.
Document Description
Reliable propagation models are required for the design of satellite broadcast to mobile or satellite navigation systems at S-band in the indoor environment. In particular, these models will support the development of hybrid satellite systems
Document Share
Document Tags
Document Transcript
  Characterisation of the Satellite-to-Indoor Propagation Channel F. Perez-Fontan * , V. Hovinen † , M. Schönhuber  ¤ , R. Prieto-Cerdeira § , J. Rivera-Castro § , P. Valtr  * , J. Kyröläinen † , F. Teschl ¤   * University of Vigo, Spain, e-mail:, † Centre for Wireless Communications, Oulu, Finland, e-mail:, ¤  Joanneum Research, Graz, Austria, e-mail: § European Space Agency, ESTEC, Noordwijk, The Netherlands, e-mail:,  Abstract   — Reliable propagation models are required for the design of satellite broadcast to mobile or satellite navigation systems at S-band in the indoor environment. In particular, these models will support the development of hybrid satellite systems including complementary ground components. Apart from the effects of shadowing and multipath encountered in LMS (land mobile satellite) urban and suburban environments, the satellite-to-indoor channel is affected by high building penetration losses, but also lower Doppler spread due to static or quasi-static situation of the user terminal. Results from a measurement campaign simulating a satellite-to-indoor link are presented. Cluster-based wideband measurements were carried out with a channel sounder placed inside several buildings (airports, shopping mall, office buildings…), receiving the signal from the transmitter located in a helicopter hovering over the building at different (azimuth, elevation) combinations in igloo configuration. The measured data has been analysed with direct methods and iterative super-resolution algorithm in order to derive building penetration losses, time spreading and angle-of-arrival information. Super-resolution techniques make it possible to identify rays and clusters of rays, and thus give added value to the development of models. The super-resolution information is effectively used for evaluation of inter-cluster and intra-cluster azimuth angle of arrival distributions. The availability of spatial information in MIMO and SIMO measurements helps developing and validating site-specific models based on ray-tracing techniques given that the measured echoes are fully characterized in terms of amplitude, phase, delay and angle of arrival. Furthermore, the parameters corresponding to a directional statistical channel model have been extracted for each scenario thanks to the availability of this type of information. The results obtained contribute to a better characterisation of the satellite-to-indoor propagation channel.  Satellite-to-indoor, wideband statistical channel model, entry loss, angle-of-arrival, time-spreading. I.   I  NTRODUCTION The development of new mobile satellite systems, including  both satellite interactive communication systems and satellite  broadcast systems to mobile users, intend to provide coverage from the satellite to users located on a wide range of environments, including indoor. ESA supports the development of precise prediction methods for satellite to indoor propagation environment, which lacks experimental data. As a first step, some statistical  propagation models have been developed for narrowband and wideband signals transmitted from a geostationary satellite. Other prediction propagation models have to be developed for non-geo satellites or to assess the coverage performance obtained with satellites plus complementary ground components. The indoor environment presents strong attenuation values introduced by building structures. However, penetration through windows or the roof still can provide sufficient signal strength inside the building. The information of multipath spreading and angle of arrival is also of relevant importance. A well recognized channel model, due to Saleh and Valenzuela [1], has been used as the baseline model for describing time dispersion effects. A modification of the above model by Spencer et al. [2] has been used for describing the spatial scattering effects, i.e. angles of arrival. II.   T HE SATELLITE - TO - INDOOR CHANNEL Two main issues are of importance when modelling the satellite-to-indoor channel, namely, (a) the entry loss and (b) the delay and angle spreading of the received energy. Entry loss is more critical of the two, given its very large values, while the delays are in the order of tens of ns. Echoes produced  by interactions (reflections, diffractions) with other buildings will, in general, be negligible in comparison to the energy entering the building by direct illumination of its façade. Entry loss is the result of very complex processes including transmission through building materials and paths through different external and internal walls, the roof, the floor, through windows, glass façades, etc. where also window glass can be metal tinted, thus giving rise to strong attenuations. Moreover, satellite paths are very much dependent on the incidence angle. 978-1-4244-2204-3/08/$25.00 ©2008 IEEE    Figure 1. Assumed cluster structure in the Saleh-Valenzuela model Channel models can be classified as empirical, statistical or deterministic / analytical. Here an empirical approach has been selected for the entry loss while a statistical model has been used for the spreading effects. The developed model and associated simulator are based on a combination of two existing models. One adapted to satellite-wireless LAN interference studies by Glazunov-Berg [3] and based on earlier COST 231 work [4]. The other model is due to Saleh-Valenzuela [1], including its directional version [2] (Modified Saleh-Valenzuela model due to Spencer et al.). The model's outcome, in its most general form, is a received echo structure (impulse response / scattering matrix) consisting of a maximum of 5 echo ensembles (four external walls and the roof). From this general result, other products: statistics, narrow-band complex time series, tapped delay-line complex time-series for link level simulations, can be generated through the use in the post-processing phase. The above assumption (up to five building entry paths) is a very general one, only rarely valid for the terrestrial case where the signal levels used are much larger than those in the satellite case. Thus, out of those five multipath ensembles, a maximum of two are expected to be observed above the measurement noise floor. Note that these ensembles are further composed of clusters of echoes and individual echoes or rays. The S&V model makes the assumption that echoes arrive in clusters, the overall impulse response being ) ττ ( δ ) jexp( β ) τ ( kl l kl l k kl   T h −−= ∑∑  φ   (1) where index l   indicates the cluster number and index k   indicates the echo number within a cluster. The amplitude (representing the received voltage) of each echo is β kl  (Fig. 1). Each external wall and the roof gives rise to an echo ensemble consisting of one or several clusters of individual echoes spread out in time (delay axis). The overall received echo structure is given by the sum of the individual complex impulse responses, h ( τ ), for all external walls and the roof, i.e., ) τ () τ () τ () τ () τ () τ ( Roof 4-Wall3-Wall2-Wall1-WallTot  hhhhhh ++++=  (2) where one, or at most two terms are significant. The S&V wideband model was extended in [2] to include Angles of Arrival, AoA. From measurements carried out in indoor environments it was found how ray clusters are not only spread out in time/delay but in angle of arrival. Thus, the srcinal Saleh-Valenzuela model was extended to include angles of arrival,  θ , as indicated by the equation ) θ ( δ ) τ ( δ ) jexp( β ) θτ ,( kl l kl l kl  l k kl   T h  θ τ φ  −Θ−−−= ∑∑  (3) where, again, index l   indicates cluster number and index kl   indicates ray number k   within a cluster l  . The model for the angles of arrival, θ , is a zero-mean Laplacian distribution with standard deviation σ , and where Θ l   is uniformly distributed  between 0 and 2 π . III.   D ATA A  NALYSIS To clarify the data analysis process, the different data sets and pre-processing steps taken before the actual analysis are summarized. The data set consists of a number of "products". The Propsound channel sounder [5] from Elektrobit was used,  providing an "instantaneous" data set consisting of one "complex" (I and Q) channel response v. delay per patch antenna and polarization. These instantaneous measurements are called cycles. The channel sounder transmits a pseudo random sequence where different code lengths are possible. Depending on the selected code length and code rate, the update rate of such set (full scan over all antennas and both  polarizations) varies giving rates ranging from several tens of cycles per second to a few hundreds. The transmit antenna was circularly polarized (RHCP) while the receive antenna consisted of a set of patch antennas (SIMO) with two orthogonal linear polarizations covering a surface that approximates a semi-sphere. From each of the two linearly cross-polarized measurements, both the received copolar (RHCP) and cross-polar (LHCP) components can be calculated. Thus, one instantaneous measurement or cycle consists of an ensemble of instantaneous individual antenna complex channel responses obtained through cross-correlation between a transmitted pseudo-random sequence and an identical, internally generated sequence in the receiver. The channel sounding process tries to measure the channel impulse response, h ( τ ). In channel modelling, it is usual to assume a sufficiently large bandwidth so that the impulse response is made up of complex deltas (ideal channel response). The band limiting effects are considered later when the impulse response is converted into a tapped delay line, TDL, fitting the bandwidth requirements of whatever system has to be simulated. In general, the results obtainable from the analyses reported should be valid for systems with RF  bandwidths smaller than that of the channel sounder, i.e., 200 MHz (chip rate 100 Mcps).  -50 0 50 100 150-50-49-48-47-46-45-44-43-42-41Excess delay (ns)    N  o  r  m  a   l   i  z  e   d  p  o  w  e  r   (   d   B   )   Figure 2. Measured and simulated / fitted power delay profile: convolution of the ideal impulse response and the response of a correlation receiver. Graz Airport, Gate 10, Elevation angle 60 degree. model TABLE I. S UMMARY OF MEASUREMENTS AND RELATED ( PRE -)  PROCESSED PRODUCTS  Instantaneous, single patch antenna, single polarization, complex channel response v. delay Cycle : ensemble of all (16 antennas × 2 polarizations) instantaneous, single patch antenna, single polarization, complex channel responses v. delay Accumulation of 5 cycles: input to super-resolution processing Averaged power delay profile, APDP Moreover, the impulse response is not constant in time but time-variant, i.e., h ( t  , τ ). This imposes some constraints on the sounding process (the duration of one full instantaneous measurement or cycle) whereby sounding has to be carried out at a rate consistent with the rate of chance of the channel (its coherence time). Assuming that the sounding process is fast enough, it is possible to consider the channel impulse response to be constant over the duration of one cycle. To sound the channel, pseudo random, noise-like sequences were used. Ideally, noise has a delta autocorrelation function. Pseudo random sequences show very narrow autocorrelation peaks: a base of the order of ± T  chip , the chip duration, and amplitude equal to the sequence length, m . The measurement of the channel impulse response is carried out by computing the cross-correlation between the received signal and an identical sequence synchronized with the one at the transmitter. The measurement process does not exactly provide the ideal channel impulse response but is approximately the result of the convolution between the channel impulse response and the code autocorrelation pulse. Given the narrowness of the cross-correlation pulse, this  processing permits a good approximation of h ( τ ). The in-phase and quadrature parts are obtained in this process. The so called  power delay profile, PDP, is in fact the most frequently used characterization parameter for the wideband channel and is given by 2 |) τ (*) τ (|log10) τ (  wh P  ≈  (4) where w ( τ ) is the channel sounder correlation pulse. Further  pre-processing of the row data includes the computation of so-called averaged power delay profiles, APDPs, encompassing a sufficient number of cycles (several tens of cycles). This is a standard procedure (ITU-R Rec. P.1407) [8] for smoothing out fading effects throughout the delay range in the channel response due to possible in-phase or out-of-phase combinations of echoes with similar delays. Table I summarizes the various measurement products. In super-resolution approach instantaneous measurements or cycles are passed on to a SAGE-based (space alternating generalized expectation-maximization) [6] ISIS (Initialization and Search Improved SAGE) super-resolution algorithm [7] for extraction of individual echoes or rays given in terms of amplitude, phase, delay, angle of arrival and departure, ) ψδ ( ψ ) θδ ( θ ) ττ ( δ ) j(exp iiiii a −−− φ   (5) The ISIS algorithm needs at least a 15 dB dynamic range to operate correctly. For instantaneous channel data suffering from low dynamic range, a coherent combination of a few cycles before applying ISIS may be used to improve the signal-to-noise ratio with respect to that of a single cycle. While ISIS still is producing good estimates for delay spread after moderate data combining, it however easily suffers from the phase distortion that most dramatically is reflected in the angle spread results. In this study, the relatively fast cycle rate allowed accumulation of data over five consecutive cycles while still preserving the phase information of the major echo contributions. Measurements were performed at fixed positions within selected buildings, with a helicopter carrying the transmitter attempting at being stationary during each measurement. At each measurement location, a sufficiently large number of instantaneous measurements (cycles) was recorded: about one minute. This allows the accumulation of sufficient data for the two types of cycle combinations pointed out above: every five cycles for super-resolution ISIS analyses and several tens of cycles for computing APDPs. Two kinds of parameters were extracted from the measurements; the first group has to do with the entry loss and the other with the time and angle dispersion effects. As an example, Fig. 2 illustrates a measured averaged power delay  profile and a simulated/fitted power delay profile resulting from the convolution of an ideal (deltas) impulse response and the response of a correlation receiver. Also in the figure, individual cluster responses are shown. For extracting the entry loss, the measured averaged power delay profiles, APDPs, were compared with a reference measurement carried out outside each building. In the beginning, to measure the entry loss, the sum of all samples above the noise floor was calculated both for the reference outdoor measurement,  P  Outdoor  , and the indoor    Figure 3. Measuring the building entry loss model -50 0 50 100 150 200 250 300-50-48-46-44-42-40-38-36-34-32-30Excess Delay (ns)    N  o  r  m  a   l   i  z  e   d   P  o  w  e  r   (   d   B   )   Figure 4. Measured and simulated / fitted power delay profile: convolution of the ideal impulse response and the response of a correlation receiver. Graz Airport, Meeting Room, Elevation angle 15 degree. measurements  P  Indoor  , i.e., ∑∑ == τ 10/) τ (Outdoor  τ 10/) τ (Indoor  Outdoor Indoor  10log10and10log10  P  P   P  P  (6, 7) This criterion was later changed to take into account only the direct, LOS contribution as illustrated in Fig. 3. These parameters are not exactly the total received power, but rather some sort of power density (dB/ns), but this does not matter since the interest is in their difference (relative measurement). The entry loss is given by ( ) Indoor Indoor Outdoor Outdoor Entry ),(  P d d  P  L −∆+=  (8) where )/log(20),( Indoor Outdoor Indoor Outdoor   d d d d  =∆  (9) is a distance correction factor needed when the reference measurement has been carried out at a different distance than that used in the indoor measurement.  P  Outdoor  ( τ ) and  P  Indoor  ( τ ) are the measured APDPs. Similarly, the averaged PDPs were normalized with respect the outdoor reference parameter  P  Outdoor  , i.e., ( ) ),() τ () τ ( Indoor Outdoor Outdoor Indoor Indoor   d d  P  P  P  ∆+−= ′  (10) This implies that the direct, LOS contribution power outdoors is assumed to be 0 dB, thus the extracted model  parameters (see next section) will be referred to the conditions  just outside the building. Their transformation into absolute values is straightforward. IV.   M ODEL PARAMETERS AND RESULTS  In this section, entry loss results for all buildings as a function of the enervation angle are provided, Table II, and Figs. 5 and 6. More details are given for two different rooms at Graz Airport, namely, Gate 10 and a conference room, are  presented. Gate 10 is a room on the ground floor of the terminal building with external glass walls, with one floor above it. On the other hand, the conference room is on the first floor with a ceiling giving directly to the terminal's roof. The extracted parameter ranges (min-max) are provided in Tables III to VII. Fig. 4 illustrates a fitted APDP corresponding to the conference room. Here the time spreading is larger than at Gate 10, Fig. 2. V.   S UMMARY  This paper presents results of a measurement campaign simulating the satellite to indoor channel at S-Band where a helicopter has been used to carry a wideband, directive, SIMO channel sounder. A model initially based on the Saleh-Valenzuela model has been developed. Preliminary wideband model parameters have also been reported. A CKNOWLEDGMENT This work has been performed under ESA / ESTEC Contract No. 19769/06/NL/GLC. R  EFERENCES   [1]   A. A. M. Saleh and R. A. Valenzuela, "A statistical model for indoor  multipath propagation," IEEE Journal on selected areas in communications, vol. 5, no. 2, pp. 128-137, Feb. 1987. [2]   Q. H. Spencer, B. D. Jeffs, M. A. Jensen and A. L. Swindlehurst, "Modeling the Statistical Time and Angle of Arrival Characteristics of  an Indoor Multipath Channel," IEEE Journal on Selected Areas in Communications, vol. 18, no. 3, pp. 347-360, March 2000. [3]   A. A. Glazunov and J. E. Berg, “Building-shielding loss modelling”, Proc. Vehicular technology Conference 2000, vol. 3, pp. 1835-1839, 2000. [4]   European Cooperation in the Field of Scientific and Technical Research (COST) Telecommunications. Digital mobile radio toward future generation systems. COST-231 Final Report, European Commission, 1999. [5]   Elektrobit. Propsopund. [6]   B. H. Fleury, M. Tschudin, R. Heddergott, D. Dahlhaus, and K. L. Pedersen, "Channel parameter estimation in mobile radio environments using the SAGE algorithm," IEEE Journal on Selected Areas in Communications, vol. 17, pp. 434–450, Mar. 1999. [7]   A. Stucki and P. Jourdan, "MIMO Radio Channel Parameter Estimation Using the Initialization and Search Improved SAGE (ISIS) Algorithm," Proceedings of Virginia Tech's Thirteenth Symposium on Wireless Personal Communications, Blacksburg, VA, June 4-6, 2003. [8]   ITU-R Rec. P.1407-3 (2007). Multipath propagation and  parameterization.  Figure 5. Entry losses for all buildings Figure 6. Entry loss values for all buildings as a function of the elevation angle TABLE II. E  NTRY LOSSES FOR ALL BUILDINGS  Max Average Min Std Graz Airport. Conference room 21.5 14.5 7.3 3.9 Graz Airport. Gate10 23.3 19.5 15.6 2.1 Graz. Shopping Center 38.3 31.2 24.2 3.7 Graz. Residential Home 33.8 29.7 23.1 2.4 Vienna. International Airport 27.3 21.7 16.7 2.4 Vienna. Millenium Tower floor 22 33.1 28.0 22.1 2.6 Vienna. Millenium Tower floor 44 27.6 22.2 14.2 3.1 Vienna. Office building 28.3 22.2 15.4 3.2 TABLE III. S UMMARY OF OTHER EXTRACTED PARAMETERS AT G RAZ AIRPORT  Elev. (deg) T  max  (ns) # of clusters Delay spread (ns) Sp (dB/ µ s) µ  (ns) 15 100-250 6-12 30.8-82.2 14.7-95.9 2-12 30 60-170 4-8 24.6-44.7 80.4-212.5 10-25 45 90-210 4-11 27.0-61.7 28.0-144.6 2-10 Graz Airport Gate 10 (ground floor) 60 105-140 4-12 27.5-68.8 60.1-115.7 2-15 15 180-360 9-17 55.6-100.2 19.2-62.9 3-12 30 305-400 14-20 75.7-95.6 35.1-50.0 2-20 45 310-350 12-15 66.2-87.6 26.1-51.2 8-15 Graz Airport Meeting room (top floor) 60 350-390 13-16 74.5-89.1 32.3-44.2 8-15
Similar documents
View more...
Search Related
We Need Your Support
Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks