Filtering Methods for Efﬁcient Dynamic Access Control in 5G Massive Machine-Type Communication Scenarios

: One of the three main use cases of the ﬁfth generation of mobile networks (5G) is massive machine-type communications (mMTC). The latter refers to the highly synchronized accesses to the cellular base stations from a great number of wireless devices, as a product of the automated exchange of small amounts of data. Clearly, an efﬁcient mMTC is required to support the Internet-of-Things (IoT). Nevertheless, the method to change from idle to connected mode, known as the random access procedure (RAP), of 4G has been directly inherited by 5G, at least, until the ﬁrst phase of standardization. Research has demonstrated the RAP is inefﬁcient to support mMTC, hence, access control schemes are needed to obtain an adequate performance. In this paper, we compare the beneﬁts of using different ﬁltering methods to conﬁgure an access control scheme included in the 5G standards: the access class barring (ACB), according to the intensity of access requests. These ﬁltering methods are a key component of our proposed ACB conﬁguration scheme, which can lead to more than a three-fold increase in the probability of successfully completing the random access procedure under the most typical network conﬁguration and mMTC scenario.


Introduction
The first phase of standardization for the future 5th generation of mobile networks (5G) was concluded in the Release 15 of the 3GPP Evolved Universal Terrestrial Radio Access (E-UTRA) specifications [1].Previous mobile generations, 3G and 4G, have merely focused on providing a faster and better mobile broadband service.Instead, 5G relies on three main use cases: massive machine-type communication (mMTC), enhanced mobile broadband (eMBB), and ultra-reliable low-latency communication (URLLC).The support for mMTC and URLLC are the main distinctive characteristics of 5G, and are expected to enable a great number of novel applications.These include, among others, the Internet-of-Things (IoT), where every daily-life object is connected to the Internet to empower a wide range of functionalities.
Nevertheless, the Release 15 of the 3GPP specifications introduced few advances in terms of mMTC support when compared to Release 14.One of the most critical requirements to efficiently support mMTC in 5G was to introduce enhancements to the random access procedure (RAP) of 4G.The RAP is the method used by user equipments (UEs) to switch from idle to connected mode.In other words, it is the method used by UEs for initial access to the cellular network.
Several studies have demonstrated the inefficacy of the RAP, as defined by the 3GPP [1], to support mMTC in 4G [2][3][4][5].Specifically, severe congestion can occur when a bulk of UE access attempts occur in a highly synchronized manner.The latter, which is a typical behavior in mMTC applications, is likely to exceed the signaling capacity of the cellular system.Since the RAP remains the same in 5G, the same limitations and problems persist.
Enhancements are expected to be introduced to the RAP in the following standardization phases of 5G.At the present time, it remains uncertain whether these will be sufficient to support mMTC applications efficiently.On the other hand, access control schemes present an appealing solution to congestion.Because of this, a wide range of access control schemes have been proposed in the literature [6][7][8][9][10][11].
One of the most promising access control schemes is the access class barring (ACB); hence, it has been included in the LTE-A Radio Resource Control (RRC) specifications [1].The ACB scheme redistributes the UE access attempts through time.For this, each UE may delay the beginning of its RAP according to the barring parameters: barring rate and mean barring time.If the barring parameters are correctly configured, sporadic periods of congestion can be relieved in exchange for a longer access delay.This is true even if barring parameters are selected beforehand and remain fixed throughout the congestion period [5,10,12,13].However, congestion is a transitory phenomenon, so an adequate performance under several conditions can only be obtained by continuously adapting the barring parameters to the intensity of access attempts.Doing otherwise would increase the access delay of UEs during intervals in which no congestion occurs, but also fail to relieve congestion episodes that are more severe than expected.Needless to say, correctly configuring the barring parameters in real-time is a complicated task.
Numerous methods to continuously adapt the barring parameters have been proposed in the literature [9,[14][15][16][17][18]; we refer to these methods simply as ACB configuration (ACBC) schemes.In other words, we use the term ACBC scheme to refer to any method or approach to select the barring parameters.In theory, some of the above-mentioned schemes can lead to a near-optimal performance in terms of success probability (i.e., probability of successfully completing the RAP) and access delay.However, most of these ACBC schemes were not designed for 3GPP cellular networks.Specifically, such a high performance is normally achieved by assuming an idealized ACB scheme and by setting extremely precise barring parameters, for which complex processes and bold assumptions are oftentimes needed.As a consequence, most of the ACBC schemes proposed in the literature cannot be implemented in the cellular systems [9,11,14,17].
Tello-Oquendo et al. presented one of the few ACBC schemes that were specifically designed for 3GPP cellular networks [18].In particular, a reinforcement learning approach was utilized.However, results show that the performance achieved under the most typical mMTC scenario with the proposed ACBC scheme is lower than that of correctly selecting the fixed barring parameters beforehand.
In our previous work, we proposed an ACBC scheme that incorporates a filtering block [19].In this block, the well-known least-mean square (LMS) algorithm was implemented.In particular, two different applications of the LMS were investigated: the adaptive line enhancer (ALE) and a novel twist on the latter, referred to as "pulling" ALE (PALE) [19]; preliminary work on the benefits of the PALE can be found in [8].During previous work, we found our ACBC scheme is one of the few that can be implemented in 3GPP networks and that can provide a near-optimal performance, which greatly surpasses that of the ACB with fixed parameters, under different access request intensities [8,19].
The main motivation of the present study is to investigate the potential benefits of our ACBC scheme with different filtering methods, namely, the recursive least-squares (RLS) algorithm or fixed finite-duration impulse response (FIR) filters instead of the LMS.Each of these presents different characteristics in terms of the rate of convergence, numerical stability, and simplicity.Therefore, the main contributions of this paper are described in the following.
1. Two different filtering methods are investigated, namely the RLS adaptive algorithm and a FIR filter.The benefits of these methods are compared to those of the LMS.The latter was the only algorithm considered in our previous work [8,19].
2. An in-depth analysis of the tracking capabilities of the LMS and RLS adaptive algorithms in transient environments is presented.
We observed several filtering methods lead to an adequate operation of our ACBC scheme under the most typical mMTC scenario.However, the main difference between implementing an adaptive algorithm and a simple FIR filter is that, for the latter method, an additional mechanism is needed to reduce the access delay of UEs in mMTC scenarios that can be handled by the RAP itself.That is, when the intensity of access attempts is low when compared to the capacity of the RAP.Clearly, increasing the access delay of UEs in these cases is not beneficial.On the other hand, such a mechanism is not needed if some adaptive algorithm configurations are implemented.
The rest of the paper is organized as follows.The RAP and the ACB scheme are described in Section 2. Next, our ACBC scheme, the studied filtering methods, and the methodology used to assess the benefits of these methods are explained in Section 3. Our main results are presented in Section 4 and Section 5 concludes the paper.

Random Access Procedure (RAP) and Access Class Barring (ACB) Scheme
This section describes the random access (RA) as defined in the 3GPP standards.It comprises the access control and the contention-based RAP [1,20].Throughout this study, the ACB, as defined in the 3GPP standards [1], is studied and we assume no other access control scheme is implemented.Figure 1 provides a brief description of the RA in 3GPP networks with the ACB scheme.The first step for a UE to switch from idle to connected mode is to acquire the network configuration parameters.Once these have been obtained, the UE is subject to the ACB scheme and, finally, it performs the RAP.The network operates in a time-slotted channel in which the minimum unit for scheduling is the subframe.
The configuration parameters are broadcast by the cellular base station (BS) through the System Information Blocks (SIBs).The SIB1 includes, among other parameters, the periodicity of other SIBs in the parameter si-Periodicity ∈ {80, 160, . . ., 5120} ms [1].The SIB2 includes some basic parameters such as the periodicity of the time/frequency resources in which preamble transmissions are allowed t rao (known as random access opportunities, RAOs) and the barring parameters [1,4,5,21].
The jth SIB2 transmission includes a barring rate p acb (j) and a mean barring time t acb (j) [1].Only special categories of high-priority UEs may not be subject to the ACB scheme.The rest of the UEs subject to the ACB scheme must perform a barring check before initiating the RAP (i.e., before the transmission of their first preamble) as described in Algorithm 1 [1,22].
Algorithm 1 Access class barring (ACB) scheme.Acquire the mean barring time t acb (j) and barring rate p acb (j) broadcast in the jth SIB2.
UEs that succeed in a barring check are no longer subject to the ACB scheme [23] and proceed to perform the RAP as follows.
Preamble: Each UE randomly selects one out of the r available preambles and sends it towards the BS in a RAO through the random access channel (RACH).Multiple UEs can access the BS in the same RAO using different preambles due to their orthogonality.As a result, BS decodes the preambles transmitted without (wireless channel) errors by exactly one UE.
In this study, we assume the BS cannot decode preambles transmitted by multiple UEs simultaneously; we refer to these as collided preambles.This assumption goes in line with the 3GPP recommendations for the performance analysis of the RACH [2] and with most of the literature [3,14,15,[24][25][26][27]. The interested reader may refer to [12] for details on the two main assumptions related to collisions during the RAP and their impact on performance.
Next, the UE waits for a predefined time window, known as the RA response (RAR) window, to receive a RAR message from the BS.

RAR:
The BS computes an identifier for each successfully decoded preamble and sends the RAR message through the physical downlink control channel (PDCCH).It includes, among other data, uplink grants for the transmission of the next message, Msg3.Up to one RAR message can be sent per subframe, but it may contain several uplink grants; each of which is associated with a successfully decoded preamble.
The PDCCH resources are limited, so a maximum of n ug uplink grants can be sent within a RAR window.Since up to three uplink grants can be sent per subframe in a RAR message, the number of available uplink grants per RAR window depends on its length.
Connection request and contention resolution: After receiving the corresponding uplink grant, the UEs schedule the transmission of the connection request message (Msg3) to the BS through dedicated resources.Then, the BS responds to each connection request with a contention resolution message (Msg4).The transmission of both of these messages is protected by a hybrid automatic repeat request (HARQ) mechanism.
The maximum number of allowed preamble transmissions for each UE is broadcast by the BS through the SIB2 [1].Access failure can occur, for example, due to a preamble collision or to wireless channel errors during contention resolution.Whenever an access failure occurs, and if the maximum number of preamble transmissions has not been reached, the UE waits for a random backoff time (determined by the backoff indicator); then, it randomly selects and transmits a new preamble at the next RAO.

Access Class Barring Configuration (ACBC) Scheme
This section describes the structure of our ACBC scheme and the configuration of the selected filtering methods.It is important to emphasize our ACBC scheme, first introduced in [8], was designed to operate under the 3GPP specifications, hence, it can be directly implemented in 4G and 5G networks.
Figure 2 illustrates the block diagram that describes the operation of the RAP with our ACBC scheme.Two main blocks are clearly identified in the figure: the RA and the ACBC scheme.The RA, depicted in the upper part of Figure 2, operates as described in Section 2.

ACB configuration (ACBC)
Random access (RA) ACB Let a(i) be the number of UEs that attempt to switch from idle to connected mode for the first time at the ith RAO; these UEs are subject to the ACB scheme.Hereafter we refer to a(i) simply as the UE arrivals.
Next, let f (i) be the number of UEs whose barring check at the ith RAO is successful.These UEs begin the RAP at the ith RAO; hence, f (i) is also the number of UEs whose first preamble is transmitted at the ith RAO.Finally, s (i) is the number of UEs that will receive an uplink grant at the ith RAR window.That is, as a response to a preamble transmitted at the ith RAO.Therefore, the discrete time index i stands for the epoch number.The epoch duration is one RAO since the RA is performed at each and every RAO.
The lower part of Figure 2 depicts the structure of our ACBC scheme, where the BS calculates the barring parameters: mean barring time t acb (j) and barring rate p acb (j) that will be broadcast through the jth SIB2.The SIB2 is broadcast once every t si RAOs, hence, these parameters are adapted according to the perceived access intensity throughout this period.As such, the discrete time index j stands for the epoch number when the epoch duration is t si RAOs.Therefore, the ACBC scheme operates at a time scale that is t si times greater than that of the RA block.Specifically, the jth SIB2 is broadcast at the i = (jt si )th RAO.Therefore, the barring parameters remain constant from the (jt si + 1)th until the ((j + 1) t si )th RAO.Throughout this study we assume t acb (j) is selected beforehand and remains fixed throughout the operation of the network.Hence, the latter is simply denoted as t acb .In the following, we describe the process to calculate p acb (j).
The BS calculates the ratio of idle to available resources immediately before the jth SIB2 transmission.For this, let n(i) be the number of contending UEs (i.e., the total number of preamble transmissions) at the ith RAO.Please observe the number of successful preambles at an arbitrary RAO i is a random variable.From there, the theoretical PRACH (uplink) capacity can be defined as the maximum expected value of the latter random variable for a given number of available preambles r and for any n(i) ∈ R ≥0 .The PRACH capacity is achieved when the number of contending UEs is and is calculated by Lin et al. [15] as Next, let c r, n ug be the theoretical capacity of the channels involved in the RAP: PRACH and PDCCH.The capacity of the PDCCH is exclusively determined by the number of available uplink grants per RAR window n ug [2].Hence, the capacity of the RAP is the minimum between the PRACH capacity and the PDCCH capacity, c r, n ug = min c (r) , n ug . (4) We refer to c r, n ug simply as the RAP capacity.With this information, the ratio of idle to available resources is calculated immediately before the jth SIB2 transmission (i.e., at the (jt si )th RAO) by means of the scaling filter block and the comparator shown in the bottom right corner of Figure 2 as where h = h 0 , h 1 , . . ., h t si −1 is the column vector of coefficients of the scaling filter and s(j) = [s(jt si ), s(jt si − 1), . . ., s(jt si − t si )] is the vector that contains the number of uplink grants assigned at each of the t si RAOs within the SIB2 period.Next, u(j) serves as the input to the filtering process, whose output y(j) is used to calculate the barring rate for the jth SIB2 broadcast interval p acb (j) = min {y(j), 1} .
Hence, p acb (j) decreases with the ratio of idle to available resources.This increases the probability of delaying the beginning of the RAP when most of the resources have been utilized.
The filtering block in our ACBC scheme (at the bottom of Figure 2) is generic, which means different filtering methods can be implemented at this point.In the following, we describe three possibilities: two adaptive filtering algorithms and FIR filters with fixed coefficients.

Recursive Least-Squares (RLS) Algorithm
The RLS is an adaptive algorithm that minimizes the total squared error between its output and the desired response.In particular, it is a recursive implementation of the method of least squares.The RLS can be derived from the least-mean square (LMS) algorithm, described in the next subsection, and is renowned for its fast rate of convergence.For instance, the rate of convergence of the RLS can be up to an order of magnitude faster than that of the LMS.In exchange, its numerical complexity is in the order of O M 2 , which is higher than that of the LMS.Nevertheless, the mathematical formulation and implementation of the RLS is relatively simple [28].
The internal structure of the RLS algorithm is illustrated in Figure 3.A buffer has been incorporated to clearly illustrate that the ratio of idle resources during the last M SIB2 intervals u(j) = [u(j), u(j − 1), . . ., u(j − M + 1)] ( 8) is the input to the algorithm.That is, a single value of u(j) is the input to the buffer as indicated by the thin arrow, and the output of the buffer is a vector, as indicated by the thick arrows.The RLS is summarized in Algorithm 2, and consists of a filtering and an adaptive process.In the filtering process, the output of the FIR filter y(j) is computed from u(j).In the adaptive process, the output y(j) is compared to the desired response d(j) to adapt the weights of the FIR filter The parameter ε(j) is the a priori error, which represents an estimate of the desired response d(j) based on the old least-squares estimate of w (j − 1).It is essential to emphasize that ε(j) is different to the a posteriori estimation error e(j), calculated in the LMS algorithm (see next subsection).
Regarding the configuration of the RLS in our ACBC scheme, the desired response is set to be d(j) = 1 for all j.By doing so, the RLS algorithm is set to minimize the distance between the output y(j) and its maximum value, d(j) = 1, which in turn "pulls" the output toward the aforementioned maximum value.This increases the barring rate p acb (j) and, in turn, reduces the access delay of UEs.This effect will be described in detail on page 16.
The RLS algorithm has two parameters that must be selected empirically, the regularization parameter δ and the forgetting factor λ. The latter determines the rate at which the effect of the previous inputs on the algorithm is minimised.There is no exact method to select δ nor λ.However, the following recommendations exist [28].

•
Regularization parameter δ: Select a small positive value when the input (in our case u(j)) is relatively high with respect to its sudden variations.Select a large positive value otherwise.

•
Forgetting factor λ: Select a positive value that is close to, but less than, 1.This ensures that past the values of the input u(j) are forgotten by the algorithm.On the other hand, λ = 1 corresponds to the method of least squares.
The beginning of Section 4 focuses on finding adequate values for δ and λ.

Least-Mean Square (LMS) Algorithm
The LMS algorithm is widely used because of its simplicity and numerical robustness [28].Concretely, the complexity of the LMS scales linearly with the filter length M. The LMS is an application of the method of stochastic gradient descent, hence, it belongs to a different family of algorithms to that of the RLS.The internal structure of the LMS is depicted in Figure 4 and Algorithm 3 summarizes its operation.

Algorithm 3 LMS adaptive algorithm.
Require: the number of filter coefficients M Require: the adaptation step size µ 1: Initialize the vector of filter coefficients w (0) and the input vector u (0) as 2: for all j = 1, 2, . . .do

6: end for
As with the RLS, the LMS comprises a filtering and an adaptive process.An important difference between these two algorithms is that the output y(j) is compared to the desired response d(j) to obtain the error e(j) in the LMS, whereas the a priori error ε is calculated in the RLS.
The correct selection of the rate of adjustment µ is essential for the LMS to be stable.In particular, the value of µ must be selected to satisfy 0 < 1 2 µ u(j) 2 < 1, for all u(j) (16) where • is the Euclidean norm operator.Two different configurations of the LMS algorithm are considered: an adaptive line enhancer (ALE) and a "pulling" ALE (PALE).In the ALE, the primary input is a delayed version of u(j), namely u(j − 1), while the desired output is set to d(j) = u(j).In the PALE, the primary input is simply u(j) and the desired output is set to d(j) = 1 for all j.The latter configuration has a similar effect on the filter output as the one described above for the RLS configuration.

Simple Low-Pass FIR Filter
Besides investigating the benefits of our ACBC scheme with the two previously described adaptive filter algorithms, we also consider the inclusion of a simple FIR filter with fixed weights w m .We refer to this as the low-pass filtering (LPF) method.
The filter weights must be chosen to ensure the following two conditions are met.
For this, let ρ be the degree of the decreasing polynomial function that determines the fixed weights of the filter.Naturally, the rapidness of the response of the fixed filter to sudden changes in the input u(j) increases with ρ, but so does its sensitivity.In other words, the selection of ρ creates a trade-off between the responsiveness and the stability of the filter output.
In particular, we choose the filter weights as where ∑ M =1 ρ is the normalization factor.

Methodology
In this study, we aim to compare the tracking capabilities of the RLS adaptive algorithm to those of the LMS and simple fixed FIR filters in a transient environment: the traffic model 2 defined by the 3GPP [2].The latter represents a highly synchronized mMTC scenario, in which the UE arrivals follow a Beta (3, 4) distribution over 10 s.Several studies have confirmed the traffic model 2 leads to severe congestion if no access control scheme is implemented [2][3][4][5]12,13].
Results presented in the following section were obtained by simulation.In particular, a C-based simulator is used, which closely replicates the arrival process of the UEs, the ACB scheme, and the RAP as described in the specifications [1,22].In each simulation, the adaptive algorithm is initialized and the filter weights are stabilized (if applicable).Then, n = 30,000 UE arrivals are scheduled within the distribution period t dist which begins at i = 0.The jth SIB2 is broadcast at the (jt si + i r )th RAO, where i r is a discrete random time shift, selected from a uniform distribution between 0 and t si − 1.A simulation run is concluded when every UE has terminated the RA procedure.The number of simulation runs is set to the smallest number that ensures the cumulative results obtained up to the last simulation differ from those obtained up to the previous simulation by less than 0.01 percent.
The first step in our analyses is to use a non-synchronized mMTC scenario, the traffic model 1, to identify an adequate configuration of the RLS algorithm.The traffic model 1 can be seen as a stationary scenario of finite duration with a low access intensity, in which UE arrivals follow a uniform distribution over 60 s.As such, the traffic model 1 (stationary scenario) is used to "train" the RLS algorithm to provide an adequate response in a transient scenario.For this, we set the UEs to ignore the ACB scheme and assess two main characteristics: the rate of convergence of the output of the filter and the relative amplitude of its variations with respect to the mean value.It is important to emphasize that similar trial and error approaches to set the configuration of adaptive algorithms are the most common in the practice, especially in transient scenarios [29].The main characteristics of both traffic models are listed in Table 1.
Table 1.Characteristics of the different traffic models defined by the 3GPP for the performance evaluation of the RAP [2].

Parameter
Traffic In our experiments, we assume the most typical configuration of the RACH, as suggested by the 3GPP [2].The number of available preambles is set to r = 54; this is the most typical value for this parameter.Furthermore, it allows us to evaluate the RAP in a scenario in which the number of successful accesses is greatly affected by r, but also by the number of available uplink grants n ug .In addition, we assume the shortest possible period for the SIB transmissions, 80 ms.These values are listed in Table 2, along with other values selected for the most relevant configuration parameters.To assess the benefits of the RLS algorithm, we test several values for parameter δ and λ, whereas other configuration parameters remain fixed.For example, M = 32 is a widely used filter length that led to satisfactory results in our preliminary tests with the RLS algorithm.Furthermore, it is the filter length that maximizes the benefits of the LMS algorithm [19].Therefore, it is selected throughout the remainder of the paper for every filter.
In addition, we select an intuitive value for the mean barring rate t acb = 1 s among those available in the SIB2; as it will be seen throughout the following section, this value leads to a sufficiently high P s ≥ 0.95 for all the selected filtering methods under the traffic model 2. Certainly, an optimal value for this parameter could be selected, but it is out of the scope of this paper to optimize the performance of each of the studied filtering methods under the selected scenario.Furthermore, selecting a value higher than t acb = 1 only causes a slight increase in the access delay when compared to the results presented in this paper.The interested reader may refer to [19], where our ACBC scheme with the LMS algorithm was optimized, and apply an analogous procedure for the RLS algorithm.The configuration parameters for the filtering methods studied in this paper are listed in Table 3.
In particular, we are set to find configurations of the RLS and fixed FIR filters that successfully reduce the variations of u(j) with the fastest possible convergence toward E [u(j)].Under the traffic model 1, n = 30,000 UE accesses are uniformly distributed within 60 s.Since RAOs occur once every t rao = 5 subframes and the subframe duration is 1 ms, there are i dist = 12,000 RAOs in the distribution period.
Let A be the random variable that defines the number of UE arrivals at an arbitrary RAO within the distribution period.Hence, a(i), shown in Figure 2, is the outcome of a single experiment for random variable A at the ith RAO.Please observe A ∼ B (n, 1/i dist ), hence, we have E [A] = n/i dist = 2.5.Since the intensity of access attempts under the traffic model 1 is considerably low when compared to the RAP capacity, the number of successful accesses per RAO s(i) highly resembles that of UE arrivals (i.e., s(i) ≈ a(i), for all i).As a consequence, E [u(j)] can be approximated by substituting s(i) with a, a column vector of length t si whose entries are E [A], in (5).
which gives E [u(j)] ≈ 5/6 for r = 54.This value has been confirmed by simulation and by an analytical model [12].

Results
In this section, we assess the benefits of our ACBC scheme with different filtering methods.As a starting point, we find adequate values for the regularization parameter δ and the forgetting factor λ, which determine the response of the RLS algorithm.Next, we briefly study the response of the LPF method (fixed FIR filters) and of the LMS algorithm.Next, we assess the benefits provided by the different filtering methods under the traffic model 2, when compared to an ACB scheme with fixed parameters (static ACB) and to our ACBC scheme with no filtering (i.e., where y(j) = u(j)).The section concludes with the comparison of the response of the different filtering methods under the traffic model 2.
Adequate values for parameters δ and λ are selected by observing the response of the algorithm.The latter is illustrated in Figure 5 for δ ∈ {0.001, 0.01, 0.1, 1} and λ ∈ {0.99, 0.999, 1}, given r = 54, along with the calculated u(j).Clearly, the response is not adequate with λ = 0.99 as p acb (j) presents a drastic drop around j = M, hence, the response of the RLS is not stable.On the other hand, the difference in the response between selecting λ = 0.999 and λ = 1 is minimal; both of these values of λ greatly reduce the variations of p acb (j) given a sufficiently low value of δ is selected, for example, δ = 0.001.Figure 5 shows δ = 0.01 is also an acceptable candidate for the regularization parameter, whereas higher values of δ should be avoided, as the amplitude of variations of p acb (j) is high when compared to its mean.It has been observed in the literature that parameter λ plays a similar role in the RLS algorithm as the complement of parameter µ in the LMS algorithm [28].In our previous studies we observed µ = 1/800 is an adequate value for the adaptation step size in our ACBC scheme.Hence, results presented in Figure 5 confirm the relation between λ and µ, as an adequate response is obtained with λ = 0.999 ≈ 1 − 1/800 = 1 − µ.The latter value of λ is selected throughout the remainder of the paper, in combination of δ = 0.001.As it can be seen, p acb (j) with this combination converges toward E [u(j)] around j = M.
Regarding the LPF method, it is immediate to see in Figure 6 that p acb (j) converges toward E [u(j)] at exactly j = M for ρ = 0, and the rate of convergence increases slightly with ρ.It is also interesting to observe that the p acb (j) calculated with the RLS, given λ = 0.999 and δ = 0.1, shown in Figure 5b, is highly similar to that of the LPF method.On the other hand, the p acb (j) calculated with the LMS presents a different behavior to both the RLS and LPF methods.For example, p acb (M) < E [u(j)] and converges toward E [u(j)] approximately until j = 100.Next, we investigate whether the difference in the response of the filters provides different benefits to our ACBC scheme.For this, we consider three key performance indicators (KPIs):

•
Success probability, defined as the probability to successfully complete the RAP.The minimum success probability considered as acceptable is 0.95.Hence, the main goal of our ACBC scheme is to guarantee that more than 95% of the UEs successfully complete the RAP.

•
Access delay, defined as the time elapsed since the arrival of a specific UE and the successful completion of its RAP.This KPI is assessed in terms of the 95th percentile.That is, the delay of 95% of the UEs that successfully complete the RAP is lower than, or equal, to the reported value.

•
Number of preamble transmissions to successfully complete the RAP.This KPI is assessed in terms of its expected value.
It is important to emphasize that the success probability under the traffic model 2 with no implemented ACB scheme is only 0.313, which is poor.In addition, the benefits of our ACBC scheme are minor when no filtering is performed and y(j) = u(j).Specifically, the success probability with t acb = 1 s and no filtering is only 0.419, which is still unacceptable.The reason for such a poor performance is that the UE arrivals and the number of successful accesses per RAO present high variations with respect to their mean value.These variations remain present in u(j) and, naturally, in p acb (j).On the other hand, filtering away these variations in u(j) leads to a stable p acb (j) that is closely related to the intensity of UE accesses in a larger time scale to that of the SIB2 period.
Figure 7 shows the 95th percentile of access delay and the expected number of preamble transmissions for the studied filtering methods.The success probability obtained with our ACBC scheme with each of the filtering methods included in Figure 7 is higher than 0.985, hence, the difference between them, in this KPI, is negligible.In contrast, the success probability under the traffic model 2 without ACB scheme is only 0.313.For benchmark purposes, Figure 7 also includes the KPIs obtained with the optimal configuration of a static ACB.The optimal configuration for the latter is defined as the combination of fixed barring parameters that leads to a success probability of 0.95 or higher with the shortest 95th percentile of access delay.In particular, we found the optimal configuration of the static ACB to be p acb = 0.31 and t acb = 1.75 s, which yields a success probability of 0.95.  Figure 7. 95th percentile of access delay and expected number of preamble transmissions obtained with our ACBC scheme and with each of the studied filtering methods; the KPIs obtained with the optimal configuration of an ACB scheme with fixed parameters (static ACB) are also included.The success probability for these is higher than 0.95.
As it can be seen, regardless of the selected filtering method, the 95th percentile of access delay and the expected number of preamble transmissions are reduced by more than 26% and 12%, respectively, when compared to the optimal static ACB implementation.It can also be seen in Figure 7 that the different values of ρ for the LPF method have a minimal impact on performance.Specifically, the difference in each of the two KPIs shown in the figure for the different valued of ρ is less than 4%.
On the other hand, there is a significant difference between the KPIs obtained with each of the adaptive algorithm configurations.In particular, the ALE leads to the longest access delay but to the lowest number of transmissions, whereas the opposite occurs with the PALE.The RLS leads to an access delay that is in between both LMS configurations, but to the highest number of preamble transmissions.This difference is caused by the "pulling" effect of the PALE and RLS algorithms, which is absent from the ALE.
Figure 8 clearly illustrates the difference in the p acb (j) calculated with the ALE, PALE, RLS, and LPF methods.The number of UE arrivals and successful accesses at each RAO with the RLS algorithm are also illustrated as a reference in Figure 8a.From this figure, it can be seen the number of UE accesses and successful accesses per RAO are highly variable.Furthermore, the number of successful accesses per RAO closely follows that of UE arrivals in the interval 0 ≤ i ≤ 343.Then, for i ≥ 343, the number of successful accesses approaches n ug .This behavior is highly desirable, as it maximizes the utilization of resources and, in turn, reduces the access delay of UEs.
Next, Figure 8b,c exhibit that the LMS algorithm with the ALE configuration operates as a low-pass filter with output p acb (j).That is, the calculated p acb (j) with the ALE configuration highly resembles that with the LPF method with both ρ ∈ {0, 2}.On the other hand, the PALE and the RLS "pull" the latter toward d(j) = 1.As a consequence, the p acb (j) calculated with the PALE and RLS are higher than those with the ALE and LPF methods.However, it is important to observe that the "pulling" effect is greater for the PALE than for the RLS.For instance, p acb (j) increases more rapidly for the PALE than for the RLS as the number of UE arrivals drops (i.e., after the 2000th RAO, which corresponds to the 150th SIB2 transmission).The reason for this effect is the rapid response of the RLS, which converges more rapidly toward E [u(j)], calculated with the M most recent values of u(j).As a result, the benefits provided by the RLS are in between those provided by the ALE and the PALE.On the other hand, the calculated p acb (j) with both ρ ∈ {0, 2} is closely similar in the period j ∈ [50, 175], and the difference is observed when rapid changes in the number of successful accesses occur.Naturally, the rate of change of p acb (j) with respect of s(i) is more rapid for ρ = 2 than for ρ = 0.However, this difference is not evident in the 95th percentile of access delay, illustrated in Figure 7.

Conclusions
In this paper, we investigated the impact of implementing different filtering methods to our ACBC scheme under a typical mMTC scenario: the traffic model 2. The latter leads to severe congestion when no access control scheme is implemented.
Based on the behavior described in Section 4, it is safe to claim that our ACBC scheme can operate efficiently with several filtering methods, and its performance is compromised only if this process is omitted.That is, when compared to the optimal configuration of a static ACB scheme, our ACBC scheme can reduce the 95th percentile of access delay by more than 26%, while increasing the success probability and reducing the number of preamble transmissions with numerous filtering methods.

Figure 2 .
Figure 2.Block diagram of the RAP with our ACBC scheme.The random access is performed at each RAO, whereas the ACBC can only be performed once every t si RAOs.

Figure 4 .
Figure 4. Block diagram of the LMS adaptive algorithm.

Figure 6 .
Figure 6.Ratio of idle to available resources u(j) and barring rate p acb (j) calculated at the jth SIB2 for a single simulation run and r = 54 under the traffic model 1 for the (a) LPF with ρ ∈ {0, 1, 2} and (b) the LMS ALE algorithm.The horizontal line is located at E [u(j)].

Figure 8 .
Figure 8.(a) Number of UE arrivals a(i) and successful accesses s(i) per RAO with the RLS algorithm, and barring rate p acb (j) calculated at the jth SIB2 for a single simulation run and r = 54 under the traffic model 2 with the (b) RLS, ALE LMS, and PALE LMS and (c) LPF with ρ ∈ {0, 2}.Results derived from a single simulation.

Table 2 .
Parameters for the selected network configuration.

Table 3 .
Parameters for the selected filter configurations.