Time-Varying Ultra-Wideband Channel Modeling and Prediction

This paper considers the channel modeling and prediction for ultra-wideband (UWB) channels. The sparse property of UWB channels is exploited, and an efficient prediction framework is developed by introducing two simplified UWB channel impulse response (CIR) models, namely, the windowing-based on window delay (WB-WD) and the windowing-based on bin delay (WB-BD). By adopting our proposed UWB windowing-based CIR models, the recursive least square (RLS) algorithm is used to predict the channel coefficients. By using real CIR coefficients generated from measurement campaign data conducted in outdoor environments, the modeling and prediction performance results and the statistical properties of the root mean square (RMS) delay spread values are presented. Our proposed framework improves the prediction performances with lower computational complexity compared with the performance of the recommended ITU-R UWB-CIR model. It is shown that our proposed framework can achieved 15% lower prediction error with a complexity reduction by a factor of 12.


Introduction
The Federal Communications Commission (FCC) defines ultra-wideband (UWB) communications [1,2] as systems that operate within the fractional bandwidth of B/ f c ≥ 20%, where B is the transmission bandwidth and f c is the center frequency of the band, or with a total bandwidth of more than 500 MHz [3,4]. The UWB systems are currently deployed on an unlicensed frequency band of 3.1 GHz to 10.6 GHz [5]. However, under the FCC definition, practical UWB communication systems may include some existing systems, such as long-term evolution (LTE) and LTE advanced (LTE-A) utilizing the single carrier (SC) frequency division multiplexing access (FDMA) scheme [6][7][8][9], 60 GHz WLAN [10,11] and future 5G wireless communications [12][13][14] operating at millimetre wave (mm-wave). According to the Mobile and Wireless Communication Enablers for the Twenty-Twenty Information Society (METIS 2020) project, the bandwidth for 5G systems ranges from 0.5 GHz to 2 GHz and operates within the candidate frequency band range of 0.45 GHz to 85 GHz [15].
The key to harvesting this high bandwidth advantage of UWB wireless communication systems is understanding their channels' characteristics and behaviors. Many works have been dedicated to studying the propagation channels of large bandwidth systems in different environments [16][17][18][19][20][21][22][23][24]. In [16], the path loss and multipath parameters were investigated based on propagation measurement using UWB signals in a forest environment. Based on the measurements conducted in an indoor environment, the work of [17] modified the CLEAN algorithm to estimate the time-of-arrival, angle-of-arrival, and waveform shape. In [18], the fundamental differences between UWB and other wireless channels, e.g., narrowband and wideband, were discussed. A three-dimensional (3D) prediction can be achieved. Second, based on our proposed models, the UWB-CIR prediction using the RLS algorithm is analyzed. The complexity reduction analysis when adopting this model is also presented. Third, the modeling and prediction performance results and the statistical properties of the root mean square (RMS) delay spread are presented. These results are generated by utilizing the proposed windowing-based models on the measured UWB-CIRs from an outdoor measurement campaign. The remainder of this paper is organized as follows. In Section 2, the UWB channel impulse response is presented. The two proposed window-based UWB-CIR models are developed in Section 3. The CIR tap prediction using the RLS algorithm for the proposed window-based CIR is presented in Section 4. Section 5 presents the performance evaluation, and Section 6 provides the computation complexity analysis. The modeling and prediction performance results are presented in Section 7. Finally, the paper is concluded in Section 8.

Channel Model
The static multipath UWB channels, assuming no pulse distortion for simplicity, can be represented using an L-tapped delay line model: where δ(·) corresponds to a Dirac delta function, and the coefficients α i and m i are the i-th tap gain and delay, respectively. For time-varying channels, (1) can be modeled as where t is either the time or the spatial location. Assuming that the transmitter moves at a constant velocity away from the receiver, the conversion between the time and spatial location parameters is straightforward [34]. Hence, the two parameters are used interchangeably in this work. For real communication systems, the space and delay must be sampled. Using T s to denote the sampling period over delay and ∆d to denote the sampling distance between the Rx and Tx antennas, the discrete form of (2) when applying the transformations of τ = mT s and t = p∆d is where m ∈ {1, . . . , N} is the delay index and p ∈ [1, . . . , P] is the position or location index for Tx or Rx antenna movement. Using the shorthand τ = mT s and p = p∆d, the discrete UWB channel model in (3) becomes

UWB Channel Measurement
In this work, real UWB-CIRs were collected from the time-domain UWB channel measurements to evaluate the performance of the proposed prediction method. The time-domain measurements were conducted in an outdoor environment considering line-of sight (LOS) scenarios using a pair of PulsON R 410 transceivers. Each transceiver had a vertically polarized omnidirectional wide-band dipole antenna mounted on it with operating frequencies of 3.1 to 10.6 GHz. The heights of the transmitter (Tx) and receiver (Rx) antennas were set to 2 m and 1.7 m, respectively. In our measurement setup, a pulse bandwidth of 2.2 GHz over the frequency range of 3.1-5.3 GHz was used. The Rx was locked at 61 ps to capture the signal transmitted by the Tx. The measurements occurred in front of the Wireless Communication Center, Universiti Teknologi Malaysia, as depicted in Figure 1. This outdoor environment consisted of a concrete floor, wooden column roofing, and a parking space with several cars. The choice of this environment was done to illustrate an infostation scenario [37,38], where users can download data as they move along the corridor. Throughout this measurement, the Rx was located at a fixed position in front of the metallic sheet. This metallic sheet was placed between the person operating the setup and the setup itself to ensure a stationary environment. The Tx was then positioned along a straight line away from the Rx, which is marked by the white tape in Figure 1. The channel realization was captured starting at a Tx-Rx separation distance of d = 3 m and at each d = 1 cm transmitter movement away from the Rx. The final channel realization was captured at a separation distance of d = 4.26 m. As a result, the measurement campaign contained P = 127 sets of measurement data corresponding to all the Tx-Rx separation distances.

UWB Channel Impulse Response Extraction
In this section, we elaborate the process of extracting the UWB channel impulse response (CIR) from the measurement waveforms. q(p) denotes the template waveform sent by the Tx at the p location index, and the delayed received waveform r(p, τ) after propagating through the channel h(p, τ) of (4) becomes where ⊗ denotes the convolution operator. The channel analysis tool (CAT) integrated into the PulsON R 410 platform was used to log and capture the received UWB radio frequency (RF) waveforms in CAT logfile format (.csv file), which is compatible for processing using MATLAB.
To obtain a high-resolution CIR, the received waveform can be further processed using the CLEAN algorithm [39]. The CLEAN algorithm is capable of resolving dense multipath components (MPCs) that are usually irresolvable by the conventional inverse filter, by deploying an iterative, high-resolution, and subtractive deconvolution procedure. To obtain h(p, τ) from r(p, τ), the CLEAN algorithm correlates the received waveform r(p, τ) with a template waveform q(p), finds the highest correlation, assigns a weight, and then subtracts the template waveform from the received waveform [16,40]. This process is repeated until the maximum iterations expire or when an MPC with gain of a certain threshold value below the strongest MPC is detected. A 20 dB threshold value is recommended by the International Telecommunication Union-Radiocommunication (ITU-R) [41].

Window-Based UWB Channel Impulse Response Proposed Model
The significant transmission bandwidth of UWB systems enables a fine-grained delay resolution at the receiver. In many environments, only some of the resolvable delay bins carry significant multipath energy, yielding a sparse channel structure. To exploit this sparse structure, our model starts by grouping the delay bins of the channel impulse response using W non-overlapping windows. Next, a channel tap selection procedure, named the maximum tap gain selection, is used to choose the channel taps within the delay bins of each window. By implementing the proposed selection procedure, each window is represented using only one tap, hence, reducing and simplifying the original L taps sparse CIR of h(p, τ) in (4) with an approximated CIR of lower W taps.

Window Selection
We further elaborated the process of grouping the delay bins (i.e., channel taps) of UWB CIR using W non-overlapping windows, which can be written as In general, the first delay bin of h(p, τ) in (4) contains the LOS component (i.e., LOS channel path), while the remaining L − 1 delay bins carry the MPCs. Hence, h(p, τ) can be divided into To maintain the channel information of the LOS channel component, the size of the first window is always set to one delay bin to represent h LOS (p, τ). Regarding the remaining L − 1 delay bins, by grouping S delay bins into one window that generate the W − 1 MPCs windows, where W − 1 = L−1 S , h MP (p, τ) can be divided using W − 1 windows given by where h MPC (p, τ) for the last window contains the remaining L − 1 − S(W − 2) channel taps. Using (7) and (8) in (6), the CIR of UWB channels of (4) can be grouped into W non-equal and non-overlapping windows as follows: where Here, the first multipath LOS component is placed in the first window as given in (10). In (11), the next multipath components are grouped using windows 2 to W-1, where each window has the same bin size of S. As an example, if (11) is used for window 3 (i.e., w = 3) and bin size of 10, the multipath components is Finally, the remaining L − 1 − S(W − 2) multipath components are grouped using the last window, given by (12).

Channel Tap Selection
As previously explained, once the delay bins have been grouped using windows, the next step is to approximate the channel taps within each window by using only one tap. This process requires the generation of two new parameters to represent the one channel tap of each window, i.e., the tap gainᾱ w (p) and delayτ w (p) for w ∈ {1, 2, . . . , W}. Using these parameters, the CIR of (4) can be approximated using In this subsection, the tap selection approach is proposed to computeᾱ w (p) andτ w (p). This approach is called maximum tap gain selection, where a significant path with maximum tap gain is selected from each of the windows to approximate the sparse CIR of UWB channels. Under this approach the maximum tap gainᾱ w (p) for each window is given by The delay of proposed window models is presented in two approaches. The first approach is to use the delay bin of the selected channel path with the maximum tap gain in each window, termedτ w (p). The proposed tap selection is called the windowing-based on bin delay (WB-BD) approach. For further simplification, in the second approach, the delay parameter is set to the delay of the last bin in each window and termedτ w,1 (p). We call this approach the windowing-based on window delay (WB-WD) approach, where the delay parameter is the window size delay defined as Representing the delay parameters with the window size delay standardizes the approximated CIRs at all location indexes p with similar bin sizes, hence simplifying the prediction of the delay parameters for future CIRs.

Channel Impulse Response Tap Prediction Algorithms
This section explains the time-domain prediction strategies that are used to predict the evolution of the UWB CIR tap. Specifically, the RLS algorithm is used to predict the simplified CIR tap of UWB channels. The aim of the prediction strategies is to predict the CIR taps associated with the future transmitted symbols, given the previous CIR taps is considered. Hence, the problem of predicting the CIR at the index p + 1 given the measured CIRs for spatial locations at indices 1, 2, . . . , p. Based on the proposed and simplified CIR using windows, two parameters must be predicted to obtain the CIR at index p + 1, i.e., the tap gainᾱ w (p + 1) and the delayτ w (p + 1) for the WB-BD approach. However, the delay parameters do not need to be predicted for the WB-WD approach. In this work, the fading process of the simplified CIR is modeled as an autoregressive process (AP) of the order p 1 ≤ p, and the estimate of both tap gain and delay parameters can be obtained usinḡ where Cᾱ ,w (p) = (Cᾱ ,w [1], . . . , Cᾱ ,w [p 1 ]) T is the prediction filter coefficient vector for the channel tap gain, and Cτ ,w (p) = (Cτ ,w [1], . . . , Cτ ,w [p 1 ]) T is the filter coefficient for the channel delay. We definē The superscripts T and H denote the transpose and Hermitian transpose operators, respectively. Here, the RLS algorithm is considered to adaptively compute the coefficients of the prediction filters in (14) and (15).
Under the RLS algorithm, the predictor filter C A,w (p) for A ∈ {ᾱ,τ} can be computed by minimizing the following cost function: Here, λ A is the forgetting factor. In this work, we set λ A to 1. The resulting RLS update equation for C A,w (p) is where the prediction error e * A,w (p) can be computed by Then, the RLS gain vector in (17) is where matrix R A,w (p) denotes the inverse of the w-th CIR taps of p 1 × p 1 dimensional sample covariance matrix, which can be computed recursively using where I p is the p × p identity matrix.

Evaluation Criterion
The proposed windowing-based UWB-CIR models and prediction techniques were evaluated by comparing the mean square error (MSE) of the waveform obtained from measurement campaigns with the waveform reconstructed using the modeled and predicted CIRs. The reconstructed waveform was generated by convolving the template pulse q(p) with the modeled and predicted CIRs. By denoting the measured received waveform by r(p, τ) in (4), the reconstructed waveform from the modeled CIR by r m (p, τ) and from the predicted CIR by r p (p, τ), the MSE between the measured and reconstructed waveform is given by for i ∈ {m, p}, and K is the total number of samples in the received waveform. To evaluate the proposed method at L location indices, the average MSE is given by The proposed UWB-CIRs models were also evaluated from the perspective of time dispersion characteristics given by the RMS delay spread parameter [42] where τ rms is the RMS delay spread that is defined by the square root of the second central moment of a power delay profile (PDP), P l is the power for the l-th path, τ l is the arrival time of the MPCs, τ b is the first path arrival time, and τ m is the mean excess delay that can be represented by the first moment of the PDP as

Complexity Analysis
This section discusses the computational complexity of our proposed windowing-based CIR models under RLS prediction algorithms, expressed by the number of basic operations (i.e., addition, division, and multiplication). Table 1 presents the computational cost of the RLS algorithm per iteration [43]. For the standard UWB-CIR with L delay bins and M prediction order, the prediction requires L(M 2 + 5M + 1) multiplication, L(M 2 + 3M) addition, and L division operations per iteration. Here, prediction is required for every delay bin, even though it does not contain any MPCs. By utilizing the sparsity of the UWB-CIR, the proposed UWB-CIR models use a lower number of W bins, i.e., W < L. For the WB-BD model, the prediction of the channel tap gain and its delay parameter is required. As a result, the WB-BD approach requires 2W(M 2 + 5M + 1) multiplication, 2W(M 2 + 3M) addition, and 2W division operations per iteration. For the WB-WD model, prediction for the delay parameter is not required as this value is further approximated to the window bin delay. Thus, the complexity of the WB-WD model is halved compared to that of WB-BD model. Here, both WB-BD and WB-WD models use fewer W bins to achieve a reduction in complexity by factors of L/(2W) and L/W, respectively, over the standard CIR model with L bins. Note that, even though the complexity analysis is performed using the RLS algorithm, the reduction in prediction complexity is achieved by adopting simplified UWB-CIR models. Hence, the models will also simplify prediction methods utilizing different algorithms. Table 1. Computational cost of the recursive least square (RLS) algorithm per iteration [43].

Results and Discussion
This section presents the modeling and prediction MSE performances along with the statistical properties of the RMS delay spread for the proposed windowing-based CIR models.

Modeling Results
Before discussing the prediction performance, we first present results to show the impact of reducing the number of MPCs to the CIR when adopting our proposed windowing-based models in terms of both the MSE and the RMS delay spread parameter. Figure 2 presents the received waveform r(p, τ) at a Tx-Rx separation distance of 4.07 m. Using this r(p, τ) measured waveform, several CIRs were extracted by using different approaches and are depicted in Figure 3. As a benchmark, Figure 3a represents the CIR extracted by the CLEAN algorithm with the ITU-R recommended threshold of 20 dB (CA-ITUR20dB). The CA-ITUR20dB CIR is very sparse; there are 380 resolution delay bins with only 23 non-zero MPCs. Note that, although there are only 23 MPCs, the prediction for CIR coefficients in the next time slot must be performed within all 380 delay bins because the locations of these MPCs are random within the delay bins. Lowering the threshold of the CLEAN algorithm reduced the channel gain gap between the strongest and weakest MPCs of the CIR, hence reducing the number of MPCs. As a result, a CLEAN algorithm with different thresholds (CA-WDT) can be used to generate CIRs with different number of MPCs. Figure 3b,c show two CIRs with reduced 9 and 7 MPCs extracted by the CLEAN algorithm with lower 15 dB and 10 dB thresholds, respectively. Despite the reduction in the number of MPCs, the complexity of channel prediction for these CIRs was similar to the CIR of the CA-ITUR20dB approach, i.e., all zero and non-zero MPCs in 380 delay bins must be predicted because their time delay parameters are random. Finally, Figure 3d shows the CIR with similar 7 MPCs when utilizing our proposed WB-BD model. We can see that the proposed WB-BD CIR in Figure 3d is more structured when compared to the CIR of CA-WDT of 10 dB in Figure 3c; this structure provides benefit to certain applications, such as channel prediction.
Let us now compare the modeling average the MSE results for the windowing-based models with the CA-ITUR20dB and CA-WDT approaches. These results are shown in Figure 4. We see that the CIR from CA-ITUR20dB approach showed the lowest MSE as a result of having the most MPCs. We can also see that the average MSE started to increase when the number of MPCs reduced for the CA-WDT approach and the windowing-model CIRs. For a fair comparison, a similar number of MPCs was set for CIR in both the windowing-based models and the CA-WDT approach. The number of MPCs was determined by setting the threshold of the CLEAN algorithm and the window size (which is defined by the number of bin per window (bpw)) for the CA-WDT approach and windowing-based models, respectively. The CIR of the CA-WDT approach introduced the lowest increment in average MSE when compared to CA-ITUR20dB. However, the average MSE increment for the proposed WB-BD model can still be considered small and acceptable for 25 bpw and lower. In contrast, the CIR of the WB-WD model showed significant increases in the average MSE compared with the other three CIRs. The MSEs of the WB-WD model were approximately 7 and 8 times larger than the benchmark CA-ITUR20dB for 5 and 35 bpw, respectively.    Figure 5 presents the RMS delay spread values for the windowing-based models with different window sizes, generated using the measured CIR at a Tx-Rx separation distance of 3.57 m. The RMS delay spread of 1.98 ns for the CIR of CA-ITUR20dB approach was used as a benchmark to determine how the RMS delay spreads of our proposed models deviate from the original CIR of CA-ITUR20dB. The RMS delay spread of the WB-BD model was slightly higher than the one achieved by CA-ITUR20dB-approximately 2.00 ns for window sizes of 5 to 25 bpw. The RMS delay spread of the WB-BD model then decreased to 1.98 ns at 30 bpw, similar to the result of the benchmarked CA-ITUR20dB. The results show that our proposed WB-BD CIR model with lower number of MPCs has a very small impact on the statistics of RMS delay spread. Hence, the system transmission rate is not impacted when adopting this less complex WB-BD CIRs. Regarding the results of WB-WD model in Figure 5, the RMS delay spread increased with the window size in the range between 2.01 and 2.65 ns for 5 and 35 bpw, respectively. The increased of RMS delay spread value reduced the system transmission rate if it was larger than the inverse of the system bandwidth. In Figure 5, the WB-WD CIR with 35 bpw has a RMS delay spread of 2.65 ns, which is 0.67 ns greater than the RMS delay spread of the benchmark CA-ITUR20dB and higher than the tolerable system inverse bandwidth of 0.45 ns. When adopting the simplified WB-WD CIR with 35 bpw, this increased of RMS delay spread by 0.22 ns above the tolerable inverse system bandwidth reduces the system transmission rate by approximately 30% when compared to the communication system utilizing the CIR of CA-ITUR20dB.
Next, Figure 6 presents the values of RMS delay spread for the windowing-based models with different window sizes generated using another measured CIR at a Tx-Rx separation distance of 3 m. A similar RMS delay spread pattern of values can be observed as the results shown in Figure 5 for a Tx-Rx separation distance of 3.57 m. Here, the RMS delay spread of CA-ITUR20dB was shown to be 1.53 ns, which is comparable to the values achieved by our proposed WB-BD model. For the CIR of the WB-WD model, the RMS delay spread values varied between 1.57 and 2.25 ns for 5 and 35 bpw, respectively. The RMS delay spread for the WB-WD model with 35 bpw was approximately 0.72 ns greater than the CA-ITUR20dB approach and 0.27 ns higher than the tolerable inverse system bandwidth of 0.45 ns. This higher RMS delay spread reduced the system transmission rate by approximately 36% when compared to the communication system utilizing the CIR of CA-ITUR20dB.

Prediction Results
We present the prediction performances, i.e., the average MSE of (22), resulting from the adoption of the windowing-based models. For prediction order M = 5, the CIR at the Tx-Rx location index of 3.05 was predicted by the RLS algorithm utilizing measured CIRs from the 5 previous Tx-Rx location indices of 3.0 m to 3.04, and then, the MSE of (21) was computed. Here, the average MSE of (22) for M = 5 in our plot was obtained by averaging the MSE values of (21) computed from predicting the CIRs at 122 location indices of 3.05 m up to 4.26 m. The average predicted MSEs under the RLS algorithm utilizing CIRs for the benchmarked CA-ITUR20dB are shown in Figure 7. For predictions using the CA-ITUR20dB CIRs, the channel tap gains in all 380 delay bins must be predicted. We see that the average MSE of CA-ITUR20dB model improved when higher prediction orders of the RLS algorithm were applied. Figure 7 also shows the average prediction MSE for our proposed WB-BD models using 5, 15, 25 and 35 bpw. We see that the average prediction MSE value adopting the WB-BD models increased with the window size because increasing the window size reduces the number of MPCs of the WB-BD CIRs. For M = 20, the average MSEs for CIR prediction utilizing WB-BD models of 5, 15, 25 and 35 bpw were approximately 1.25 × 10 −3 , 1.37 × 10 −3 , 1.5 × 10 −3 , and 1.87 × 10 −3 , respectively. In general, the average MSE for CIR prediction adopting the WB-BD model improved with the prediction order; this improvement decreased when a higher window size was used. For example, the average MSE for the window size of 35 bpw remained comparable at approximately 1.87 × 10 −3 , even at higher prediction orders, and for the window size of 5 bpw, the average MSE improved from approximately 1.35 × 10 −3 for M = 5 to 1.25 × 10 −3 for M = 30. The main advantage of the proposed WB-BD model for prediction is the reduction in computational complexity via the lower number of MPCs. As shown in Figure 7, our proposed WB-BD models not only requires a lower prediction complexity, but it also outperforms the benchmarked CA-ITUR20dB CIR at the lower prediction order region. For example, when considering the RLS prediction order M = 5, the WB-BD model with 25 bpw achieved a lower average MSE of 1.5 × 10 −3 , whereas the CA-ITUR20dB model achieved a higher average MSE of approximately 1.77 × 10 −3 . Here, the proposed WB-BD model achieved a 15% lower prediction error rate with a prediction complexity reduction factor of 12, because at 25 bpw, the simplified WB-BD CIR is represented by only 15 MPCs.
Next, Figure 8 presents the average prediction MSEs under the RLS algorithm utilizing CIRs from the WB-WD models for different window sizes of 5, 15, 25, and 35 bpw. By approximating the time delay parameters to the window size using the proposed WB-WD models, the prediction complexity reduced further by a factor of 2 when compared to the WB-BD models; however, this reduced complexity slightly degraded the prediction performances. For the comparison at the prediction order of 20 and a window size of 15 bpw, the average MSE of the WB-WD model in Figure 8 was approximately 7.8 × 10 −3 , and the average MSE for the WB-BD model in Figure 7 was approximately 1.38 × 10 −3 .

CDF of RMS Delay Spread
Next, we considered the RMS delay spread distribution for the proposed WB-BD and WB-WD models. First we considered the distribution of the modeled windowing-based approaches. Figure 9 shows the cumulative density function (CDF) comparing the RMS delay spreads for the WB-BD and CA-ITUR20dB models. We found that the distribution of RMS delay spreads for the WB-BD model using 5, 15, and 25 bpw were comparable with the CIR of CA-ITUR20dB. This finding implies that the simplified WB-BD model maintains the second-order characteristics of the UWB channels. The CDF of the RMS delay spreads for the second WB-WD model is depicted in Figure 10. Note that, by further approximation of the time delay parameter with the window size, the CDF of RMS delay spreads for the WB-WD CIRs deviated from the CA-ITUR20dB model. This deviation increased with an increased window size. Here, 90% of the RMS delay spread values for CIRs of CA-ITUR20dB were lower than 1.9 ns, and 90% of the RMS delay spreads of the WB-WD CIRs were below 2 ns, 2.1 ns, and 2.2 ns for the window sizes of 5, 15, and 25 bpw, respectively.  We then considered the distributions of the RMS delay spreads for the predicted CIRs using the WB-BD and CA-ITUR20B models. Figure 11 shows the CDFs of the RMS delay spreads for the predicted CIRs using the RLS algorithm with M = 10 for the WB-BD and CA-ITUR20dB models. Here, only the CDFs of delay spreads for the WB-BD models with 5 and 15 bpw were considered because they have lower average prediction MSE values than the CA-ITUR20dB model. The delay dispersion has great importance for high-data rate systems, where the error probability due to delay dispersion is proportional to the RMS delay spread. Figure 11 shows that 90% of RMS delay spreads utilizing the proposed WB-BD models with 5 and 15 bpw were less than 1.6 ns. For the CA-ITUR20dB models, 90% of the RMS delay was less than 2.0 ns. The results show that the predicted CIRs using the WB-BD models with 5 and 15 bpw not only improved the prediction performance but also reduced their delay dispersions provided by the RMS delay spreads compared to the predicted CIRs using the CA-ITUR20dB model. Figure 11. CDF of the RMS delay spread for the predicted CIRs using the recursive least square (RLS) algorithm with M = 10.

Conclusions
This paper applied real CIRs collected from measurement campaigns in an outdoor environment to investigate a time-domain based channel prediction for UWB channels. A less complex and simplified CIR prediction method involving the exploitation of the sparse property of UWB channels was proposed by introducing two new models: the WB-BD and WB-WD UWB-CIR models. The new models were generated by grouping channel taps of ITU-R UWB-CIRs with a 20 dB threshold using non-overlapping windows, and consequently, two procedures were introduced to select channel taps within each window. The MSE performances and RMS delay spread values were used to evaluate the proposed methods. Based on the modeling and predicted MSE performances, the WB-BD model outperformed the WB-WD model. Based on RLS algorithms, the channel prediction utilizing the WB-BD model was shown to outperform not only the WB-WD model, but it also achieved a 15% lower prediction error rate with a reduction factor of 12 in the prediction complexity over the ITU-R UWB-CIR. In addition to the lower prediction performance and complexity when using the WB-BD model, the predicted CIRs were found to possess smaller RMS delay spread values. These results highlight the potential of the UWB WB-BD channel model for critical applications such as the future 5G systems requiring both high data rate and low-latency.