Abstract

Compared with the line-of-sight (LOS) condition, the multipath effect is more serious in the non-line-of-sight (NLOS) condition. Therefore, the LOS and NLOS identification is necessary for the multipath analysis of signal propagation. The commonly used method is the support vector machine (SVM) method with high computational complexity. To tackle this problem, this paper adopts the SVM classifier based on fewer selected features of the normalized power delay profile (PDP). Therein, the PDP can be obtained using the sliding correlation method. The results show that the SVM-based classifier can achieve high accuracy on LOS and NLOS identification. We then analyze the impact of the signal-to-noise ratio (SNR) and transmitting-receiving (Tx-to-Rx) distance on distinguishable multipaths under LOS and NLOS conditions. According to statistical measurement results, a function of distinguishable multipath numbers is established. Finally, we investigate the multipath power and delay parameters of average delay spread and root mean square (RMS) delay spread based on multipath results. The outcomes of this paper provide a useful support for analyzing signal propagation characteristics.

1. Introduction

Communication systems can be developed in various environments such as tunnels instead of being limited to free space [1, 2]. Therefore, it is necessary to explore signal propagation characteristics in narrow spaces [3, 4].

Due to the confined space of a tunnel environment, internal wireless signals are susceptible to be interfered caused by non-line-of-sight (NLOS) effects, resulting in amplitude attenuation and phase shift of signal propagation [5, 6]. Therefore, it is significant to identify line-of-sight (LOS) and NLOS conditions. The ranging-based estimation methods are investigated in [79], and the power-based methods are adopted in [10, 11]. Additionally, deep learning methods are used for identification. For example, an overview of opportunities and challenges of the deep learning-based classification method is presented in [12]. The multimodal deep learning method is proposed to design the mobile encrypted traffic classifier based on automatically extracted features [13].

Moreover, the LOS and NLOS conditions are identified in several systems. For example, for the ultrawideband (UWB) system, the support vector machines (SVMs) are developed in [1416], with the results showing that the NLOS condition can be effectively identified based on UWB signal features. In [17], the multipath condition is added to be considered, i.e., the machine learning techniques are compared to identify LOS, NLOS, and multipath conditions, with the results showing the effectiveness of the machine learning method in the field of the UWB indoor localization. For the vehicle-to-vehicle (V2V) system, in [18], the SVM-based method is used via combining the angular information and conventional channel features to identify LOS and NLOS conditions. In [19], the machine learning-based LOS identification methods are compared based on angular properties of the dynamic channel impulse response (CIR), which effectively improves the identification accuracy. Inspired by these, we aim to investigate the LOS and NLOS identification method for the tunnel environment.

In addition to LOS and NLOS identification, a considerable amount of research has been conducted on the characteristics of wireless channels in tunnels. For example, the 2.4 GHz band signal is used for a subway in Madrid, and the arrival time distribution of multipath delay clusters is analyzed under LOS, NLOS, and far-line-of-sight (FLOS) conditions [20]. The measurement results of [21] provide insights into the large-scale fading characteristics in actual curved subway tunnels with distances of 300 m and 500 m at 920 MHz, 2.4 GHz, and 5.7 GHz. Zhou et al. [22] investigated measurement data on radio frequency (RF) propagation in tunnel and mine environments for vertically, horizontally, and cross-polarized signals; their results showed the importance of tunnel size, frequency, polarization, and electrical characteristics of the tunnel wall on signal propagation characteristics. Radio channel characteristics such as path loss, root mean square (RMS) delay spread, channel stationarity, Doppler shift, and channel capacity are investigated in a subway tunnel [23].

In this paper, we adopt the SVM-based method to identify LOS and NLOS conditions for the tunnel environment. Compared with the used SVM in [14, 16], we only select variance, mean excess delay, root mean square (RMS) delay spread, and kurtosis of the normalized power delay profile (PDP) as identification features. The main contributions of our works are outlined as follows:(i)We only select four features to identify LOS and NLOS conditions with the SVM method for the tunnel environment, which significantly reduces the computational complexity.(ii)We analyze the change in distinguishable multipaths by comparing the normalized PDP of measurement data under LOS and NLOS conditions. Moreover, a function of distinguishable multipath numbers is established according to the statistical results. Additionally, the multipath power and delay parameters are analyzed under the LOS and NLOS conditions.(iii)All measurement data are available to be downloaded at https://github.com/sunjc-cumt/measurement-data. To the best of our knowledge, this is valuable to provide downloadable actual measurement data for the tunnel environment.

The rest of this paper is organized as follows. In Section 2, the system model is presented. In Section 3, the SVM classifier and identification method are presented. In Section 4, the design of the measurement system and the measurement scenarios are described. In Section 5, the measurement results are analyzed. Finally, conclusions are presented in Section 6.

Notations: the notations , , and denote the scalar, vector, and matrix, respectively. The operator represents the convolution. denotes the transpose of matrix. and denote the absolute value and -norm, respectively. denotes the smallest integer greater than . denotes the expectation operator, and denotes the sign function. .

2. System Model

Due to the confined space of a tunnel environment, obvious refraction and reflection phenomena occur during signal propagation [24], as shown in Figure 1. The transmitted signal travels multiple paths to the receiver. The general formulation of a CIR in a tunnel environment [25] can be given bywhere represents the total distinguishable multipath numbers and , , and represent the amplitude, delay, and phase of the -th path, respectively. The phase is assumed to be a priori, statistically independent uniform random variable distributed over [26].

PDP is important in the analysis of channel characteristics, which can be obtained by the sliding correlation method. The transmitted signal is assumed to be an m-sequence modulated by QPSK. By not considering the phase factor in (1), the received signal can be given aswhere denotes the additive white Gaussian noise (AWGN). The received signal is correlated with the known transmitted signal by using the sliding correlation method. Based on the stable autocorrelation of an m-sequence, we can obtainwhere represents a cross-correlation function between the received and transmitted signals and represents an autocorrelation function of . A correlation function profile consisting of multiple peaks is obtained, which is known as the PDP. Subsequently, the channel characteristics can be represented by PDPs.

3. LOS and NLOS Identification

The SVM classifier is adopted to identify LOS and NLOS conditions based on selected features. For reducing computational complexity, we only select variance, mean excess delay, RMS delay spread, and kurtosis of the normalized PDP as identification features. The form of the SVM classifier is given bywhere

and are unknown parameters which are determined based on , is the number of reference points, and denotes the selected features, i.e., variance, mean excess delay, RMS delay spread, and kurtosis, and is a predetermined function, and with

The SVM classifier aims to maximize the margin with slack between two classes [27], which is presented aswhere are slack variables and is the penalty term.

Generally, the Lagrangian dual of the optimization (7) is a quadratic program (QP), which can be solved with the sequential minimal optimization (SMO) method [28, 29]. However, for complexity consideration, this method is usually replaced by the least squares (LS) method [30]. The LS-SVM is used to turn the inequality constraint into an equality constraint, which facilitates the solution of the Lagrange multiplier. Consequently, the problem is a linear equation in the LS-SVM:where denotes the error.

Based on the principle of the LS-SVM, the method of the identification of LOS and NLOS conditions is given in Algorithm 1. In the training phase, the Lagrangian dual of optimization problem (8) is given bywhere . denotes Lagrange multiplier. Then, by setting partial derivatives of the Lagrangian dual to parameters , , , and to zero, we can obtain values of , , and

Input: training set ,
Testing data .
Output: identification results.
\\Training phase:
(1)Obtain the Lagrangian dual based on equation (9);
(2)Set partial derivatives of the Lagrangian dual to parameters , , , and to zero;
(3)We can obtain
\\Testing phase:
(4)For the given testing data , we calculate the output based on equation (11);
(5)Identify LOS and NLOS conditions based on equation (6).

In the testing phase, based on trained parameters, we can obtain the identification result for the given testing data , which is given bywhere denotes the kernel function. Typically, the linear, poly, Gaussian, and radial basis function (RBF) kernel functions are adopted to solve (11).

4. Measurement System and Environment

4.1. Measurement System

In this section, we describe the measurement system and environment. The designed measurement system is shown in Figure 2.

In this study, we adopt a RF vector signal generator (N5172B), which can generate 2.4 GHz Wi-Fi to transmit signals and adopt a vector signal analyzer (N9020B) to receive signals. Additionally, we use a BNC line as trigger to connect the transmitting and receiving instruments for time synchronization.

The transmitted signal (m-sequence) and trigger signal are imported into N5172B through processor-side Toolkit software and sent through the transmitting antenna after setting transmission power and sample rate. Then, N9020B receives and stores signals through the receiving antenna after time synchronization. It is worth pointing out that the directions of the receiving antenna and the transmitting antenna are always aligned throughout the experiment. Finally, we use the sliding correlation method to obtain PDPs at the processor side. The parameters of the transmitting and receiving antennas are summarized in Table 1.

4.2. Measurement Environment

A typical underground tunnel environment is applied as the measurement environment. The tunnel is arched with closed glass doors at both ends and two ramps of the same size in the middle section.

Figure 3(a) shows a long straight tunnel, wherein the width of both ramps is 5.8 m. The other side of the ramp is an auxiliary road with a width of 1 m; it runs parallel to the tunnel. Figure 3(b) shows an arched cross section with a side height of 2.3 m, a vault height of 3.6 m, and a width of 3 m. The measurement setup is shown in Figure 3(c).

For convenience of position description, each test point in Figure 3(a) is given in the form of two-dimensional coordinates. Tx1 is taken as the reference point, and the coordinate is set to in meters. The coordinates of reference points are listed in Table 2.

5. Analysis of Measurement Results

5.1. LOS and NLOS Identification

We measure received signals at the given Rx points under LOS and NLOS conditions. In order to ensure the accuracy of the analysis, we collect data under LOS conditions and data under NLOS conditions. Then, the SVM classifier is adopted to identify LOS and NLOS conditions. Figure 4 demonstrates the receiver operating characteristic (ROC) curve of four kernel functions on identification. The black dashed line is the identification ability baseline, i.e., there is no identification ability in this line. The result shows that the SVM classifier with four kernel functions has a high accuracy on identifying LOS and NLOS conditions.

Numerically, we present the area under the ROC curve (AUC) values as shown in Table 3. The results intuitively present that the AUC values can be up to 0.96 and 0.97. Therefore, each kernel function can be adopted to aid the SVM classifier for identification and has a high accuracy.

The computational complexity of the SVM classifier is divided into the training phase and testing phase. In the training phase, generally, we have , and the computational complexity is , where is the dimension of the selected feature set, is the number of the support vector, and is the number of reference points. In the testing phase, the computational complexity is , where denotes the number of operations required to evaluate the kernel. Therefore, the total computational complexity is [31]. In comparison to previous works [14, 16], we only select four features in our works, i.e., , and the dimension of the feature set is lower. Therefore, the computational complexity is obviously reduced.

5.2. Analysis of Distinguishable Multipaths

In this section, the effects of various factors on the distinguishable multipaths are analyzed under LOS and NLOS conditions. It is worth pointing out that the secondary peak can be considered as a distinguishable multipath component when the peak value is within 30 dB (refer to recommendation ITU-R P.1407-4 [32]; the threshold used for the identification of the number of multipath components depends on the dynamic range of the measurement equipment; a typical value is 20 dB in the outdoor environment. According to our measurement environment, equipment, and data, the threshold is set to 30 dB in this paper) lower than that of the main peak in each PDP.

5.2.1. LOS Conditions

The receiver can receive direct signals from the transmitter under LOS conditions. Some examples are presented: (i) the transmitting and receiving positions are set at Tx1 and Rx7, respectively, with a distance of 30 m, as shown in Figure 3(a). The transmission power is set to −10 dBm and 0 dBm to obtain two PDPs. The PDPs can be compared as shown in Figure 5(a). (ii) The transmission power is fixed to 0 dBm, and the receiver is placed at Rx7 and Rx8, while the transmitter is placed at Tx1. The change in the distinguishable multipath numbers is analyzed with respect to Tx-to-Rx distances, as shown in Figure 5(b).

As illustrated in Figure 5(a), only the direct signal can be distinguished, and the nondirect component cannot be distinguished when the transmission power is set to −10 dBm with receiving signal-to-noise ratio (SNR) being 8 dB. However, the two components can be distinguished with a delay of 30 ns when the transmission power increases to 0 dBm with receiving SNR being 17 dB. It can be concluded that the nondirect components can be easily distinguished at the receiver with the SNR increasing.

Figure 5(b) shows that the two paths can be distinguished with a delay of 30 ns when the distance is 30 m. The three paths can be distinguished with a delay of 20 ns and 40 ns when the distances increase to 40 m. It can be concluded that when the transmission power is sufficiently high, the distinguishable multipath numbers increase with distances increasing within a certain range.

5.2.2. NLOS Conditions

The receiver can only receive refracted and reflected signals from the transmitter under NLOS conditions. Similarly, several examples are conducted under NLOS conditions: (i) the transmitter and receiver are placed at Tx2 and Rx3, respectively. The PDPs in the LOS and NLOS scenarios are compared at the same distances, as shown in Figure 6(a). (ii) The transmitter is placed at Tx2, and the receiver is placed at Rx5 and Rx6, respectively. A comparison of the PDPs is shown in Figure 6(b) (for convenience, we mark the distance of the x-axis of Tx-to-Rx in Figure 6(b), and the real Tx-to-Rx distance can be calculated by coordinates in Table 2). (iii) The transmitter is placed at Tx1, while the receiver is placed at Rx10 and Rx11, respectively. The measurement paths look like “L” and “U” types, respectively, as shown in Figure 3(a). The multipath effects of two types of measurement conditions are compared and analyzed, as shown in Figure 6(c).

Figure 6 shows that, in the NLOS scenario, the distinguishable multipath numbers reduced compared with the LOS scenario with the same transmission power and distances. Furthermore, the effects of NLOS cause path loss to increase as the distances increase, and the receiving signal is severely distorted. Finally, in “U” type measurement conditions, signal power is drastically attenuated after traveling through two walls, and the multipath signal cannot be clearly distinguished. It can be concluded that the receiver cannot receive the information effectively with an increase in the number of obstructions.

Then, we consider that the transmitter and receiver are placed in a straight tunnel that is fully occluded by the obstacles, i.e., an environment with closed glass doors is shown in Figure 3(a). The measurement setup is shown as follows. (i) The transmitter is placed at Tx1, while the receiver is placed at Rx1 and Rx6, respectively. The obtained PDPs are compared with those obtained under LOS conditions with the transmission power of 10 dBm and distance of 20 m, as shown in Figure 7(a). (ii) The PDPs are analyzed under different transmission powers with obstacles, while the transmitter and receiver are placed at Tx1 and Rx2, respectively, as shown in Figure 7(b). Therein, receiving SNR is 0 dB with the transmission power of −20 dBm, and receiving SNR is 9 dB with the transmission power of 10 dBm.

Figure 7 shows that the signal is distorted as it passes through the obstacles with a low SNR. The distances of signal propagation are shorter than those in the LOS because of the presence of the obstacles. The solution is to increase the SNR to ensure efficient transmission of the signal.

Based on the above measurement results, the distinguishable multipaths in the receiver are related to SNR, Tx-to-Rx distance, and measurement scenarios. The LOS and NLOS conditions are considered when analyzing the change in distinguishable multipaths.

According to statistical results, the distinguishable multipath number can be expressed as a function within the error tolerance by using the fitting method, which is given bywhere in meters denotes the Tx-to-Rx distance, , , and denote coefficients, and denotes the bias.

Based on statistical measurement results, we determine the coefficients and bias as shown in Table 4. Therefore, the distinguishable multipath number can be determined using the given function.

5.3. Analysis of Multipath Power

Based on multipath results, we investigate the CDF of the multipath power in the normalized PDP, which is shown in Figure 8. The overall observation from this figure, in comparison to the NLOS condition, shows that the multipath power of the LOS condition is higher. In detail, at a CDF of , the multipath power is less than −23 dBm under the LOS condition, while that is less than −30 dBm under the NLOS condition. The reason is that there is no direct path, which causes that the signal attenuation is more severe under the NLOS condition.

5.4. Analysis of Multipath Arrival Delay

In channel multipath propagation, signals traveling through different paths arrive at the receiver at different times, resulting in a delay spread of the receiving signal. In general, the propagation is characterized by delay parameters including average delay spread and RMS delay spread .

The average delay spread is defined as the first order of the power delay spectrum:where is the power of the -th distinguishable multipath with the relative delay (the relative delay denotes the delay of the distinguishable multipath component relative to the main peak) .

The RMS delay spread is defined as the square root of the second moment of the power delay spectrum:where

Based on results of the above distinguishable multipath, channel delay parameters are analyzed at different locations under LOS and NLOS conditions, which are shown in Figure 9.

Figure 9 indicates that and of the NLOS are larger than those of the LOS. Moreover, the average delay spread increases as the distance increases, whereas the change in the RMS delay spread is not obvious. This reveals a low correlation between the RMS delay spread and the distances for both LOS and NLOS conditions [33].

6. Conclusion

In this paper, the SVM classifier was adopted to identify LOS and NLOS conditions. The results showed that the SVM with four conventional kernel functions can achieve a high identification accuracy. Moreover, a function of distinguishable multipath numbers was established. The coefficients of the function were determined based on statistical measurements. Furthermore, the analysis of the multipath power reveals that the multipath component signal is severely attenuated under the NLOS condition. The multipath delay parameters indicate that the average delay spread increases as the distance increases, while the RMS delay spread has a low correlation with the distance. The results of this study provided a theoretical basis for developing communication systems for tunnel environments.

For the future avenue of our works, we aim to consider more reliable and low-complexity identification methods, e.g., deep learning methods. Therefore, the multipath of the tunnel environment can be analyzed more effectively.

Data Availability

The data used to support the findings of this study can be downloaded at https://github.com/sunjc-cumt/measurement-data.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work of Shuai Ma and Shiyin Li was supported by the National Natural Science Foundation of China (Grant nos. 61701501 and 61771474) and in part by the Shandong Provincial Major Scientific and Technological Innovation Project (MSTIP) (Grant no. 2019JZZY020505), the Open Fund of IPOC (BUPT), the Open Research Fund of National Mobile Communications Research Laboratory, Southeast University (Grant no. 2021D02), Young Talents of Xuzhou Science and Technology Plan Project (Grant no. KC19051), and the Program for “Industrial IoT and Emergency Collaboration” Innovative Research Team in CUMT (Grant no. 2020ZY002).