Statistical Analysis of Radio Propagation Channel in Ruins Environment

The cellphone based localization system for search and rescue in complex high density ruins has attracted a great interest in recent years, where the radio channel characteristics are critical for design and development of such a system. This paper presents a spatial smoothing estimation via rotational invariance technique (SS-ESPRIT) for radio channel characterization of high density ruins. The radio propagations at three typical mobile communication bands (0.9, 1.8, and 2GHz) are investigated in two different scenarios. Channel parameters, such as arrival time, delays, and complex amplitudes, are statistically analyzed. Furthermore, a channel simulator is built based on these statistics. By comparison analysis of average excess delay and delay spread, the validation results show a good agreement between the measurements and channel modeling results.


Introduction
As mobile phone is becoming more and more widely used nowadays, a phone-based localization system for search and rescue in emergency situations, such as earthquakes, avalanches, and landslides, has received considerable attention recently [1][2][3][4].To develop and implement such a system, a comprehensive study on the radio propagation behaviors in disaster scenarios is required, especially for system configuration and localization.However, there is a lack of comprehensive study on this issue so far.For pedestrian localization of victims in ruins, the radio propagation characteristics between the cellphone and handheld detection device were studied in [5].In [6], the base station (BS) to mobile station (MS) radio propagation channel was investigated in two simple artificial ruins, where the rubbles covered only a small area around transmitter in an open field.In actual disaster scenarios, large area of ruins is very common, and the disaster channel suffers extremely complicated fading.To evaluate the influence of large amounts of scatterers and to ensure effective implementation of the phone-based localization system, further researches are required to characterize the radio propagation behaviors in ruins, and investigation of channel modeling method for such complex scenarios is very critical.
For characterization of the radio propagation channel, wideband measurement can be conducted to derive the impulse responses by parameter estimation.Traditionally, the inverse Fourier transform method is utilized for this application.However, the order of Fourier modeling is too large and less efficient for simulation [6,7].Over the past years, conventional spectral estimation methods [8] have been applied for channel modeling.In [7], a low-order autoregressive model was employed to model the indoor radio propagation channel.In [9], three parameter estimation methods, including minimum norm [10], autoregressive model, and frequency response matching method, were compared.In the dense ruins environment, a lot of interferences show up due to poor channel conditions.Under this circumstance, the aforementioned parameter estimation methods are no longer effective.In this study, the ESPRIT method is applied for modeling of the radio propagation channel in dense ruins.As a high efficiency spectral estimation approach, the ESPRIT method has been utilized for resolution enhancement in several specific situations, such as in [11,12], where it was applied to estimate the time delays for ground penetrating radar.
In this paper, the SS-ESPRIT approach is utilized to derive the parameters of channel impulse responses in ruins.Through wideband measurements, the parameters are statistically analyzed in two typical ruins scenarios at three communication bands: 0.9, 1.8, and 2 GHz.Based on the statistics, a channel simulator is built and the statistical properties of two key channel parameters, including average excess delay and delay spread, are analyzed.
where ,   ,   , and   are the number, amplitudes, phases, and delays of multipath components, respectively, and () denotes additive white Gaussian noise.

ESPRIT-Based
where , and Let and ; then we can derive where The ESPRIT estimation is a nonparametric method, which is based on eigenvalue decomposition of a specific correlation matrix, and the largest  eigenvalues correspond to signal subspace, whereas the remaining eigenvalues correspond to noise subspace.The   in Φ can be derived directly by three ESPRIT-based approaches: total least squares [13], time-shift [14], and spatial smoothing.In this study, the SS-ESPRIT method shows a superior performance compared to the other two, which will be illustrated in Section 4. For channel modeling, the number of multipath components, , should be determined first and the minimum description length method in [15] is recommended for estimation, which is given as follows: where  is the independent samples,  is the total number of eigenvalues of the correlation matrix,  1 () and  2 () are the arithmetic and geometric averages of the smallest  −  eigenvalues, respectively, and  3 (, ) = 0.5(2 − ) log().
After estimating  and   in ( 1), the values of   and   are determined by solving the least squares solution of the equation H = Ac, which is referred to as extended Prony's method [9].

Model Validation by Simulation.
To validate the SS-ESPRIT method, a simulation is performed, and it is compared with the other two parameter estimation methods: Burg-autoregressive method and minimum norm.The normalized root mean square error (NRMSE) is utilized to evaluate the performance of these methods, which is calculated by the following expression: where Ĥ(  ) is the reconstructed responses based on estimated parameters and (  ) is the measured responses.The simulated channel frequency responses are at the 900 MHz frequency band with a bandwidth of 400 MHz.Moreover, the number of multipath components is set to be 5 and the delays and amplitudes are uniformly distributed in the range between [100 ns, 150 ns] and [0, 0.6], respectively.Based on these assumptions, the "measured" channel frequency responses are constructed according to (1), which is referred to as "original data" in the following description.By 200 times of simulation, the average NRMSEs and model orders at different signal-to-noise-ratios (SNRs) are summarized in Table 1 these methods can perfectly estimate the channel parameters at high SNRs, but at low SNRs, the SS-ESPRIT approaches show a lower model order and a better antinoise performance.Figure 1 depicts a comparison example between the original channel frequency response and the SS-ESPRIT reconstructed responses, in which a good agreement can be found.Furthermore, another advantage of SS-ESPRIT is that almost no disturbance multipath components show up in the estimated channel impulse response at low SNRs.

Measurement Campaign
A chaotic yard of collapsed buildings and pieces of furniture, with an area of about 100 m × 20 m, was selected for measurements.The channel frequency responses were measured in two different scenarios along a designated path from 10 m to 28 m away from the transmitting (Tx) antenna at regular spatial intervals of 0.2 m, as shown in Figure 2, in which the circle represents measurement point.At each distance, three locations at 0.5 m spacing were measured to gain enough spatial independent samples.Thus, a total of about 270 spatial samples at each scene were obtained.To simulate the mobile base station, the receiving (Rx) antenna was mounted on a 4 m high pole, which was moved from one measurement point to another during the measurement.In the two scenarios, the Tx antenna was placed behind ruins and then covered by stratified media.To reduce the interference, all the mobile phones around were turned off, and an average of 200 snapshots for each sample was taken by vector network analyzer (VNA).
During the measurements, the channel responses from 700 MHz to 2200 MHz with an interval of 0. frequency interval is the time length of signal, a 2 s time window can be obtained for 0.5 MHz interval, which is long enough to analyze the impulse responses in this study.Two zero gain omnidirectional biconical antennas, which were vertically polarized, were utilized to transmit and receive signal.

Results and Discussion
Radio channel in the test field is a static problem, which is in line with the scenarios after disaster.The time-varying factors, such as moving rescuers, are not taken into consideration in this study.As the disaster channel is complicated and time-invariant, it is viewed as a wide sense stationary uncorrelated scattering channel.Based on the experimental data, three mobile communication frequency bands are studied: 0.8-1 GHz, 1.7-1.9GHz, and 1.9-2.1 GHz.A bandwidth of 200 MHz is utilized for parameter estimation according to analysis of the measured data.deviation of the peaks mainly results from spectral leakage and interferences.Besides, it is found that a lot of small spurious interferences show up in the estimated channel impulse responses due to complex multipath fading, which leads to the bad performance of conventional modeling methods.
In order to reduce the unreasonable results, multipath components after 400 ns, whose amplitudes are significantly small in all cases, are removed from the measured channel frequency responses by 15th order Butterworth filter.Consequently, the SS-ESPRIT estimated impulse responses are truncated and the model orders are reduced from about 200 to 35.To characterize the radio propagation behavior in ruins, statistical properties of the other three channel parameters,   ,   , and   , are analyzed in the next part.

Statistical Results of the Delays and Amplitudes.
In our investigation, the delays, amplitudes, and phases of multipath components, show no regularity, which are viewed as random variables, and statistical analysis is performed for characterization.Distribution functions (normal, exponential, Rayleigh, Nakagami, Rice, and Weibull) are applied to fit the aforementioned parameters, and the Akaike Information Criterion is utilized to compare the fitting performance.Besides, the best fitting distribution function for each parameter is validated by the Kolmogorov-Smirnov goodness-of-fit test.As the spatial samples are measured at different locations, the delays and amplitudes of multipath components should be normalized first.The normalized delays   , namely, excess delay, are obtained by subtracting the line-of-sight propagation delay from the measured   , whereas the amplitudes   are normalized by dividing the theoretical values of free space propagation.
The excess delay of the first arrival  fa is a very important channel parameter, especially for localization methods based on arrival time.As shown in Figure 5, statistical analysis indicates that exponential distribution fits the measured cumulative distribution function (CDF) of  fa best; that is,  fa ∼ ().Besides, the parameters of the fitted distributions are summarized in the third row in Table 2.The mean value of exponential distribution is exactly the parameter  and it is found that the mean of  fa is quite small.In general, if a lot of large obstacles show up, the arrival time would become longer.Furthermore, the normalized amplitudes of the first arrival  fa is normally distributed; that is,  fa ∼ (, ).And the parameters of the fitted normal distribution are summarized as well in the last two rows in Table 2.In this investigation, no explicit relationship between the normalized amplitudes and  fa is found and the expected phenomenon in which the amplitudes decrease with increasing  fa does not show up.Moreover, the signal energy is larger around the arrival time and it tends to reduce after arrival.Therefore, to describe the behavior of the normalized   elaborately, the 400 ns time window is divided into five groups according to the difference between the excess delay and  fa .Among these groups, given in Table 3, Group 2 is the arrival time interval.Besides, the time length of Group 4 is set very large as subdivided groups in this range obey almost the same distribution.As shown in Figure 6, the best description of the normalized   in each group is Rayleigh distribution and Table 3 summarizes the corresponding distribution parameters under different conditions.Because of more serious blocking effect, the amplitudes of Group 2 in Scenario 1 are much smaller than in Scenario 2.
In addition, the phases of multipath components-  in each group are uniformly distributed in the range [0, 2] in this study.Besides, the excess delay of each multipath component obeys uniform distribution as well in the observed time interval.

Building and Validation of the Channel
Simulator.Simulation of radio channel is always a research hotspot, which can help configure the localization system and evaluate the performance of base station to mobile station communication as well as localization.Based on the statistical analysis of the parameters of channel impulse responses, a channel simulator is built and here are the steps.
(1) First, the normalized   and excess delay of the first arrival are determined as illustrated in Table 2.
(2) Then, a set of normalized   in each interval group is generated, according to the Rayleigh parameters in Table 3.
(3) Afterwards, a total of  excess delays, which are uniformly distributed from 0 ns to about 400 ns, are (5) Finally, the "normalized" channel impulse response is constructed based on the multipath channel model in (1).
To characterize the radio propagation in ruins and validate the channel model, two important channel parameters  "Measured" means that the data were obtained by experiments whereas "regenerated" means that they were obtained by channel simulator.
is within about 10%.Generally, the error is acceptable for simulation analysis of channel characteristics.In addition, the delay spread values in the two scenarios are both about 100 ns.

Conclusion
In this paper, a SS-ESPRIT based characterization of radio propagation channel in high density ruins is presented.By measurements in a large yard full of collapsed buildings, the radio propagation channel has been investigated at three mobile communication bands (0.9, 1.8, and 2 GHz) and the results show the SS-ESPRIT method is applicable for such analysis.Furthermore, the delays, amplitudes, and phases of multipath components are statistically analyzed and different distribution functions are applied for description of these parameters.Based on these results, a channel simulator is established and it is validated through the analysis of two key channel parameters (average excess delay and delay spread) where the experimental data show an acceptable agreement with the channel modeling results.As a conclusion, this paper provides an effective method for modeling of communication channel in high density ruins, and the analysis results reported in the paper can be utilized for design of the phonebased localization system.

2 International
Journal of Antennas and Propagation

FrequencyFigure 3 :Figure 4 :
Figure 3: Comparison of the measured data and the SS-ESPRIT reconstructed signal in Scenario 1.

Figure 5 :
Figure 5: Empirical cumulative distribution and fitted exponential distribution of the excess delay of the first arrival at 2 GHz in Scenario 1.

Figure 8 :
Figure 8: Comparison of the empirical distribution and fitted normal distribution of delay spread at 900 MHz in Scenario 1."Measured" means that the data were obtained by experiments whereas "regenerated" means that they were obtained by channel simulator.

Table 1 :
Comparison of the NRMSEs of different methods at different SNRs.
Channel Parameter Estimation.Through measurements, we can obtain the discrete channel frequency responses (  ) with  = 1, 2, . . .,  and   =  0 +  ⋅ Δ, where  0 is the start frequency and Δ is the sampling interval.Afterwards, parameter estimation method is applied to determine the four parameters in (1): ,   ,   , and   .
. It is observed from Table 1 that almost all Figure 1: Comparison of the simulated channel response and the SS-ESPRIT reconstructed responses at 10 dB SNR.

Table 2 :
Distribution parameters of the excess delay and the normalized amplitude of the first arrival.

Table 3 :
The normalized   , the parameters of the fitted Rayleigh distributions in two scenarios.

Table 4 :
Average excess delay and delay spread, the parameters of the fitted normal distributions.

Table 4 .
For average excess delay, the measured parameters show a good agreement with the regenerated parameters.From Table4, it is found that the average excess delay in Scenario 1 is about 130 ns, which is slightly higher than Scenario 2 due to more serious blocking effect.As to delay spread, a little deviation of the regenerated data (about several nanoseconds) is found in Table4, which