Wideband and Ultrawideband Channel Models in WorkingMachine Environment

We present statistical models for wideband and ultrawideband (UWB) radio channels in a working machine cabin environment. Based on a set of measurements, it was found that such a small and confined space causes mostly diffuse multipath scattering rather than specular paths. The amplitude of the channel impulse responses in the wideband case is mostly Rayleigh distributed small-scale fading signal, with only a few paths exhibiting Ricean distributions, whereas the ones in the UWB case tend to be log-normally distributed. For the path amplitude, we suggest an exponential decay profile, which has a constant slope in dB scale, with the corresponding parameters for the UWB case. For the wideband case, a twofold exponential decay profile provides excellent fits to the measured data. It was also noted that the root-mean-square (RMS) delay spread is independent of the line-ofsight/obstructed line-of-sight situations of the channel. The multipath components contributing significant energy play a major role in such a small environment if compared to the direct path. In addition, the radio channel gains are attenuated with the presence of a driver inside the cabin.


Introduction
A wireless sensor network (WSN), which is a wireless network consisting of spatially distributed autonomous devices using sensors to cooperatively monitor physical or environmental conditions [1], today covers very widely in many industrial and civilian application areas, including industrial process monitoring and control, health monitoring, environment and habitat monitoring, home automation, and traffic control.Wireless intra-vehicle communications system has been proposed as a new type of communications system that can provide onboard passengers with bidirectional highspeed data exchange services for both communications and entertainment [2,3].For the purpose of safety, comfort and convenience, WSNs are being deployed in the new models of automotives to collect such information as temperature, speed, pressure, and so on.This idea of exploiting WSNs is not limited to applications in a working machine only.The information and communication system used in a working machine makes it easier for driver to control the operation, for example, in a wood processing guiding the trunk cutting according to orders input from the system.This can guarantee quick processing during logging, and high productivity.This information system is expected to be wirelessly controlled by the driver.Using wireless communication links instead of wired link can save a big amount of installing costs.In addition, measuring vibrations acting on the human body, so-called "human vibration," is the prevention of health risks and the evaluation of comfort, for example, in vehicles [4].Wireless vibration sensing application embedded to wireless body area networks (WBANs) can also be included in the future system.
The design of information and communication systems requires a good understanding of the corresponding radio wave propagation channel.The characteristics of an indoor radio propagation is found to be its site-specific behavior.Many existing radio channel measurement campaigns and the corresponding channel models have been restricted to typical residential and office environments, which are highly limited by the attenuation of walls [5,6].Only few radio channel measurements were done inside a vehicle [7,8].To our knowledge, there is no work done in a working machine environment, which can have the different radio channel characteristics due to its small and confined space.This paper discusses on experimental radio channel measurements that are carried out in a working machine environment within two different frequency bands, that is, the industrial, scientific, and medical (ISM) radio frequency band (2.4-2.4835GHz) and the ultrawideband (UWB) frequency band (3.1 and 10 GHz).The statistical radio channel models based on the measurement results are then developed.The novelty of this paper comes from the new environment whose radio channel characteristics are modeled.
Section 2 of the paper is devoted to the description of the measurement setup including the scenarios.Section 3 addresses the characterization of the measured radio channels.The average channel impulse responses for each scenario are analyzed and discussed.The parameters of the statistical channel model are also presented in this section.Section 4 deals with the implementation of the channel model and its evaluation.The conclusions of this work are given in Section 5.

Measurement Campaigns
2.1.Measurement Setup.The measurements were carried out in a frequency domain using Agilent 8720ES S-parameter vector network analyzer (VNA) [9].The network analyzer is operated in a transfer function measurement mode, where port 1 and port 2 are the transmitting and the receiving ports, respectively.Thus, the measurements result as a frequency response of the channel.
The measurements were carried out in the frequency band between 2.4-2.4835GHz, which is in the ISM band, and is used, that is, by the standards IEEE 802.11 and 802.15.4 dedicated for wireless local area network (WLAN) communication.The measurements in an UWB frequency band between 3.1-10 GHz were also conducted.Therefore, the measured bandwidths B were 83.5 MHz and 6.9 GHz, respectively.Consequently, the corresponding delay resolutions are 11.97 ns and 0.14 ns.The maximum number of frequency points per sweep M is 1601, which can then be used to calculate the maximum detectable delay τ max of the channel as Using (1), the maximum detectable delay, τ max of the channels for two different bandwidths are 19.16μs and 231.88 ns, which correspond to 5.75 km and 69.56 m, respectively, in a free-space distance.Maximum detectable delay is just a measurement system-based parameter.It is introduced to show that the measurement system is able to detect all propagation paths having significant contribution to the received power.However, the usable signal energy is in the shorter delays than the maximum detectable delay as can be seen later on.
The antennas for the measurements in WLAN band 2.4-2.4835GHz were HGA7S 7 dBi high gain antennas [10].The SkyCross SMT-3TO10 M-A antennas [11] were used for the measurements in UWB band.Both types of the antennas are linearly polarized.The SkyCross SMT-3TO10 M-A antennas

Channel Measurement Description.
The measurements were conducted inside a cabin of a working machine with the dimensions as illustrated in Figure 1.The transmit (Tx) antenna and the receive (Rx) antenna are, in general, set to face each other, and the line-of-sight (LOS) path always exists except when the links are blocked by the seat.The arrangement of the antennas' positions is shown in Figure 1.The position of the Tx antenna is fixed at the seat on 1.08 m above the cabin floor.The positions of the Rx antenna were moved to 17 different spots, namely, from the position Rx1 till the position Rx17.At each position, 100 consecutive measurements of the channel were taken to improve statistical reliability.This number of typically used in static cases.In order to study radio channels inside a cabin in a working machine for real situations, the following three scenarios are investigated: (i) scenario 1: without a driver, when the engine is off.
(ii) scenario 2: with a driver sitting, when the engine is off.
(iii) scenario 3: with a driver sitting, when the engine is on.

Channel Measurement Results and Data Analysis
Measured S 21 -parameters, that is, channel transfer functions, are converted to the time domain, that is, to channel impulse responses h(τ), using an inverse fast Fourier transform (IFFT).A Hamming window is used to reduce sidelobes.The channel impulse response is generally given by a tapped delay line model [12] as where a l is the amplitude of the lth path, τ l is the arrival time of the lth path, φ l is the phase for the lth path, and L is the number of paths.δ(t) stands for the Dirac function, which is the impulse symbol and is sometimes called as "Dirac's delta function."The phase φ l is modeled by a uniform distribution U{0, 2π}.The statistical parameters defining the channel model are presented in the following paragraphs.

Average Channel Impulse
Response.One hundred individual realizations of the channel impulse response were averaged for each position.The reason for using one hundred realizations is since this is typically used to guarantee statistical reliability.This number of sweeps is recorded for all the individual links measured, so the total number of sweeps carried out during the measurement campaign is much higher.When using frequency domain measurement technique, the environment needs to be static, which means that one hundred consecutive measurements per link is enough.
To proceed further analysis, all measured channel impulse responses are first truncated above the noise threshold, set to the noise level adding 10 dB (about −98 dBm).The only paths above the noise threshold are taken into account.Some examples are presented here in order to explain the analysis.The large-scale features (i.e., path loss) and also small-scale features (due to the change in pulse shape) of the channels are taken into account.The large-scale models are necessary for network planning and link budget design while small-scale models are necessary for efficient receiver design and performance analysis.This kind of models is not publicly available at open literature for working machine's cabin environment so far.
As preliminarily shown in [13], the measured channel impulse responses for the ISM band have two constant slopes in dB scale, which correspond to a typical exponential decay profile in a linear scale.This shows that there are many reflections rather than a significant specular path.This is excellently suited to characterize diffuse scattering [14].
The average channel impulse responses for the UWB band (6.9 GHz bandwidth) in scenario 1, without a driver, when the engine is off shown in Figure 2 were also represented in [13].The UWB measurements allow us to check the presence of many paths arriving at the receiver within the delay resolution and the maximum delay.Firstly, the channel impulse responses have a constant slope in dB, which corresponds to a typical exponential decay profile in a linear scale.The other advantage of the UWB measurements is the high temporal resolution, which makes it possible to see the first path more precisely.The first path arrival time can approximately translate into the distance between Tx and Rx antennas.Therefore, the distances between Tx and Rx1, Rx3, Rx8 as well as Rx17 are 1.17 m, 1.35 m, 1.05 m, and 0.81 m, respectively (3.9 ns, 4.5 ns, 3.5 ns, and 2.7 ns, resp.).The Rx1 and Rx3 are in the obstructed line-of-sight (OLOS) case; thus the first path is not the strongest path.Depending on the position, the constructive multipaths can create a stronger path to the receiver.The Rx8 and Rx17 are LOS cases; therefore, the first path is also the strongest path.Figure 3 illustrates the average measured channel impulse responses for the ISM band, Rx5 (OLOS), in three different scenarios.When the driver is inside the cabin, the radio waves in all multipaths from the Tx close to the body of the driver (≈3 cm) are attenuated.In addition, the first path arrival time is longer than the one in the case of no driver inside the cabin.Note that in the following, the analysis has been carried out against several distributions, and the best fits have been selected.Moreover, all the parameters used in the models were extracted from the measured data during the postprocessing.From the extracted parameters, we develop measurements-based deterministic channel models for both ISM band and UWB band cases.

Exponential Decaying Factor (a l ).
As explained above, the path amplitude a l in ( 2) is approximately modeled by an exponential decay profile with a Ricean factor γ 0 and an exponential decaying factor Γ for UWB case as ( For ISM case, the channel impulse responses have two constant slopes in dB, and thus the path amplitude a l in ( 2) is approximately modeled by twofold exponential decay profile with two Ricean factors γ 01 and γ 02 as well as two exponential decaying factors Γ 1 and Γ 2 as where l 1 is the last path of the first part.An example from each frequency band is presented in Figure 4 showing the average channel impulse response for the ISM band and for Rx17 (LOS) in the first scenario  (without a driver, the engine is off) and the twofold least square (LS) linear fitting curve in dB scale.The corresponding factors γ 01 and γ 02 are −44 dB and −81 dB, respectively, as well as the exponential decaying factor Γ 1 and Γ 2 are 2.33 and 24.39, respectively.Figure 5 shows the average channel impulse response for the UWB band and for Rx17 (LOS) and the LS linear fitting in dB.The corresponding factor γ 0 and the exponential decaying factor Γ are −66 dB and 21.74.

Amplitude Variation (χ).
This part represents the amplitude variation χ over the mean value of the path amplitude.The cumulative distribution functions (CDFs) in Figures 6 and 7 depict that the amplitude variations of the  average channel impulse responses are well fitted by the lognormal distribution with zero mean and standard deviation σ χ of 6.69 dB and 4.50 dB, respectively.

Number of Paths (L).
The CDFs of the number of paths L (above the noise threshold) are excellently fitted by Poisson distribution with a certain mean μ L .An example for Rx14 (LOS) in the UWB band shown in Figure 8, respectively.The number of paths decreases with the presence of the driver.In the case of the ISM band, the CDFs of the number of paths L are fitted by Gaussian distribution with certain mean μ L and standard deviation σ L as an example in Figure 9 for Rx14 (LOS).The number of paths presented in the paper is for channel modeling purpose.Nonetheless, they are much more than the typical number of paths in system design, for example, the number of paths within 10 dB of the strongest path and the number of paths containing 85 percent of the energy, which give approximately 35-50 paths in the UWB band case.As can be seen, the impact of the vibration (of the engine) on the number of the paths is clear, especially, in the case of OLOS unlike the other channel characteristics.

Delay Dispersion.
In order to compare different multipath channels and to develop general design guidelines, parameters which grossly quantify the multipath channel are used.The mean excess delay τ m and root-mean-square (RMS) delay spread τ RMS are commonly used to imply the time delay dispersive properties of a wideband multipath channel parameters that can be determined from a power delay profile.The received signal will suffer spreading in time if compared to the transmitted signal.This dispersion can also imply frequency selective fading and intersymbol interference (ISI).No serious ISI is likely to occur if the symbol duration is longer than, say, ten times the RMS delay spread [14].The RMS delay spread is the standard deviation (or root-mean-square) value of the delay of reflections, weighted proportional to the energy in the reflected waves defined as [14] where τ m is expressed as Figures 10 and 11 depict the CDFs of the measured RMS delay spreads for LOS and OLOS situations, respectively, for the UWB band measurements.All values are in the range of 1.7-5.8ns, which are very small if compared to typical indoor RMS delay spread varying between 14 and 18 ns [14].We can see that the RMS delay spread is independent of the LOS/OLOS situations of the channel.Neither of the distance, this is because the cabin space is very small.The multipath components contributing significant energy play a major role in such a small environment if compared to the direct path.This is very important on the selection of the antenna position and also to the receiver structure.An example for the case of the cabin with a driver, when the Rx antenna is at Rx12 (LOS), is shown in Figure 12.The RMS delay spread decreases by approximately 1.2 ns, since the presence of people in the cabin attenuates the path gains.Corresponding CDFs for OLOS case (Rx4) are shown in Figure 13.The standard deviations σ of the RMS delay spreads, when the engine is on (cabin vibration), are larger than the ones, when the engine is off.This shows that the radio channels differ considerably for the case with and without driver.

Amplitude Distribution.
In the ISM band, the power received for each delay tap is the vectorial sum of many multipath components arriving in the corresponding delay bin.The central limit theorem is therefore valid, and the amplitudes of the delay bins exhibit a Rayleigh (in nonLOS or OLOS case) or a Ricean (in LOS case) distribution [15].On the contrary, only few multipath components (MPCs) fall into each resolvable delay in the UWB band.Therefore, the central limit theorem is not applicable and the amplitudes of the delay bins consequently do not exhibit a Rayleigh or Rice distributions.
In the ISM band case, in all scenarios, the first delay bins tend to be Ricean distribution due to their high K factors [16], whereas the other delay bins can be assumed to be Rayleigh distributed due to their very small K factors.
For LOS case in the UWB band, as explained above, only few MPCs fall into each resolvable delay and the central limit theorem is therefore not applicable.The measured amplitudes are rather log-normally distributed [17].

Large-Scale Channel Characterization.
Path loss (PL) is an important large-scale characteristic for link budget calculation and system design.In the conventional narrowband channel, Friis transmission formula is used to model the received signal power, P r in free space given by [18] P r (d where P t is the transmitted power, and G t and G r are the Tx and the Rx antenna gains, respectively.λ is the wavelength, and d is the Tx-Rx separation distance.The Friis transmission formula shows that the received signal power falls off as the square of the Tx-Rx separation distance.From (7), the term PL is given by (λ/4πd) 2 which predicts that signal power will decrease with the square of increasing frequency and distance.These show the presence of both distance and frequency dependency in path loss.Nevertheless, the frequency is constant within a bandwidth of interest and is thus negligible.However, for UWB systems, the frequency dependency in PL can be significant due to the large bandwidth (>500 MHz).The channel would distort the signal spectrum of the signal which is directly related to the signal distortions, that is, pulse shape distortion.PL can be obtained directly from the measured channel transfer function H( f ).The distance dependency of the path loss is calculated by time snapshot and frequency averaging over the total bandwidth at each Rx position as where H( f i , t j ; d) denotes the jth channel transfer function time snapshot at frequency f i in a distance d.N and M denote the number of sample points within the band and the number of snapshots, respectively.The mean PLs for both LOS and OLOS cases are calculated from (8) and averaged over the results from all Rx positions in each case.The results for each case and each scenario are summarized in Table 2.As expected, the PL increases in the OLOS case and with the presence of a driver.Moreover, it is noted that the PL in the case when the engine is on is not very different than in the case when the engine is off.This follows the initial assumption.the engine is on and off, only the number of the paths are different, especially, in the case of OLOS; that is, the vibration (of the engine) has only the impact on the number of paths, whereas the impact on other channel properties is trivial.
In both frequency ranges, the exponential decay profile was modeled for the path with the corresponding parameters.For the amplitude distribution, many MPCs arrive in the corresponding delay bin in the ISM band, the central limit theorem is therefore valid, and the amplitudes of the delay bins exhibit a Rayleigh (in non-LOS or OLOS case) or Ricean (in LOS case) distributions.On the contrary, only few MPCs fall into each resolvable delay in the UWB band.Therefore, the central limit theorem is not applicable and the amplitudes of the delay bins consequently do not exhibit a Rayleigh or Ricean distribution, but rather log-normal distribution.For the delay dispersion, the RMS delay spread was independent of the LOS/OLOS situations of the channel.Neither of the distance, this is because the cabin space is very small.The MPCs contributing significant energy play a major role in such a small environment if compared to the direct path.This is very important on the selection of the antenna position and the receiver design.For the path loss, the frequency dependency of the path loss in the UWB band was firstly observed.Secondly, the mean path loss increases in the OLOS case and with the presence of the driver.Moreover, it was noted that the mean path loss in the case, when the engine is on, is not very different than in the case, when the engine is off.The measurements obtained were used to estimate the corresponding channel parameters in order to build statistical simulated channel models.
The developed radio channel models give a great additional information to wireless communications designs in a working machine cabin environment.

Figure 2 :
Figure 2: Average channel impulse responses for the UWB band in scenario 1: without a driver, when the engine is off.The positions of Rx antenna are Rx1 and Rx8.

Figure 3 :
Figure 3: Average measured channel impulse responses for the ISM band, Rx5 (OLOS), in three different scenarios.

Figure 4 :
Figure 4: Average channel impulse responses for the ISM band for Rx17 (LOS) and the least square linear fitting curve in dB scale.

Figure 5 :
Figure 5: Average channel impulse responses for the UWB band for Rx17 (LOS) and the least square linear fitting curve in dB scale.

Figure 6 :
Figure 6: CDFs of the amplitude variation over the mean value of the path amplitude and the log-normal fitting curve for the ISM band and for Rx17 (LOS).

Figure 7 :
Figure 7: CDFs of the amplitude variation over the mean value of the path amplitude and the log-normal fitting curve for the UWB band and for Rx17 (LOS).

Figure 8 :Figure 9 :
Figure 8: CDFs of the number of paths L for Rx14 (LOS) in the UWB band measurements.

Figure 10 :
Figure 10: CDFs of the measured RMS delay spreads for LOS in the UWB band measurements.

Figure 11 :Figure 12 :
Figure 11: CDF of the measured RMS delay spreads for OLOS in the UWB band measurements.

2 .Figure 13 :
Figure 13: CDFs of the measured RMS delay spreads for Rx4 (OLOS) in the UWB band measurements in all three scenarios.

Figure 16 :Figure 17 :
Figure 16: Average measured and simulated channel impulse responses for Rx17 (LOS) in the ISM band.

Table 2 :
Average path loss (PL) within the measured band.