DOI QR코드

DOI QR Code

Grant-Free Random Access in Multicell Massive MIMO Systems with Mixed-Type Devices: Backoff Mechanism Optimizations under Delay Constraints

  • Yingying, Fang (Nanjing University of Posts and Telecommunications) ;
  • Qi, Zhang (Nanjing University of Posts and Telecommunications)
  • Received : 2022.04.24
  • Accepted : 2023.01.07
  • Published : 2023.01.31

Abstract

Grant-free random access (GFRA) can reduce the access delay and signaling cost, and satisfy the short transmission packet and strict delay constraints requirement in internet of things (IoT). IoT is a major trend in the future, which is characterized by the variety of applications and devices. However, most existing studies on GFRA only consider a single type of device and omit the effect of access delay. In this paper, we study GFRA in multicell massive multipleinput multiple-output (MIMO) systems where different types of devices with various configurations and requirements co-exist. By introducing the backoff mechanism, each device is randomly activated according to the backoff parameter, and active devices randomly select an orthogonal pilot sequence from a predefined pilot pool. An analytical approximation of the average spectral efficiency for each type of device is derived. Based on it, we obtain the optimal backoff parameter for each type of devices under their delay constraints. It is found that the optimal backoff parameters are closely related to the device number and delay constraint. In general, devices that have larger quantity should have more backoff time before they are allowed to access. However, as the delay constraint become stricter, the required backoff time reduces gradually, and the device with larger quantity may have less backoff time than that with smaller quantity when its delay constraint is extremely strict. When the pilot length is short, the effect of delay constraints mentioned above works more obviously.

Keywords

1. Introduction

Machine-type communication (MTC) is one of the major tends for future wireless networks and the driving factor Internet of Things (IoT) [1]. Massive multiple-input multiple-output (MIMO), where each BS equipped with dozens or even hundreds of antennas serve tens of users, is regarded as a key support technology for MTC [2]. Massive MIMO can greatly increase the system capacity and has high robustness and low energy consumption [3,4]. It is also charactered by low cost and interference tolerance [5].

MTC has vast amount of devices that is larger than the available pilots. However, the conventional massive MIMO usually assumes a fully-loaded of devices, which is obviously not appropriate in MTC [6,7]. Therefore, we need to design an effective access scheme for devices in massive MIMO. The existing random access (RA) schemes of massive MIMO systems can be divided into two categories: grant-based and grant-free RA (GFRA) schemes [8-10]. In grant-based RA, the devices need to wait for the permission of the BS to transmit data, which will obviously increase the access delay [11]. In the GFRA scheme, the device directly transmits the pilot sequence and data to the BS without waiting for permission from the BS [12]. The detailed process of GFRA is given in the 3rd generation partnership project (3GPP) [13]. Obviously, compared with the grant-based RA scheme, GFRA scheme can reduce the transmission delay and signaling cost, which can satisfy the low lost and short delay requirement of MTC, and becomes popular in current research [14,15]. In [16], it proposes a GFRA scheme to reduce satellite access delay. In [17], it proposes a RA scheme of preamble-data superposition in massive MIMO system, which significantly improved access success rate and performance. These studies show that GFRA significantly improves the performance of massive MIMO system. Therefore, lots of works began to focus on the combination of massive MIMO and GFRA. The preamble code design scheme is proposed to achieve high success rate of device detection and channel estimation in massive MIMO GFRA system in [18]. In [19], it confirms the existence of an optimal pilot length that maximizes spectral efficiency in GFRA of massive MIMO systems. In [20], it analyzes the parameters that affect system capacity in massive MIMO GFRA system, such as the pilot length. However, all of the above studies considered single-cell scenarios. In [21], GFRA in multicell massive MIMO is studied. The backoff mechanism is usually applied on RA in 3GPP [22]. As the backoff parameter can be adjusted flexibly, some works focus on obtaining the optimal backoff parameter for GFRA [23, 24]. However, these works regard devices that select the same pilot as indistinguishable. In multicell systems, due to the large coverage of network, devices selecting the same pilot may located far away from each other, and the BS can detect the device under its coverage based on the distinction of received powers. Hence, a more precise analysis method is needed for the random access in multicell massive MIMO systems. Moreover, to enhance the device access success probability, we adopt the pilot reuse scheme which makes devices selecting the same pilot can be detected due to the low-level of interferences. Moreover, previous works on GFRA only consider one type of devices. In IoT, various kinds of applications result in different types of devices. Therefore, customized RA strategies should be provided for each type of devices according to their unique configurations and requirements. In addition, IoT devices usually have strict requirement on access delay. Therefore, we take the access delay into account, and set a delay constraint when optimizing the backoff parameter. This has also not been performed in existing literatures.

In this paper, we investigate GFRA in multicell massive MIMO systems with mixed-type devices with delay constraints. With backoff mechanism, each device is randomly activated, and the activated device randomly selects orthogonal pilot in the orthogonal pilot pool to send frequently to the BS. We deduce the average spectral efficiency of each type of devices and the average spectral efficiency of all devices in each cell, and carry out simulation verification from all aspects. Based on it, we obtain the optimal backoff parameter for each type of devices under their delay constraints. It is found that the optimal backoff parameters are closely related to the devices number and delay constraints. In general, devices that have larger quantity should have more backoff time before they are allowed to access. However, as delay constraints become stricter, the required backoff time of devices reduces gradually, and the device with larger quantity may have less backoff time than that with smaller quantity when its delay constraint is extremely strict. Moreover, when the pilot length is short, the effect of delay constraints mentioned above works more obviously.

For ease of understanding, some of the symbols in this article are summarized in Table 1.

Table 1. Significance of symbols

E1KOBZ_2023_v17n1_185_t0001.png 이미지

2. System Model

2.1 Channel Model

Consider a multicell MIMO wireless network with L cells as illustrated in Fig. 1, where each cell consists of a BS equipped with M antennas and N uniformly distributed devices with single antenna. There are K types of devices in the network, and the number of the k-th type of devices is Nk. That is, Σk=1K Nk = N. The n-th device of the k-th type in the l-th cell is denoted as Unkl. Devices transmit their data in the same time-frequency resource. The M × 1 channel vector ginkl between the i-th BS and Unkl is given by

\(\begin{aligned}\mathbf{g}_{i n k l}=\mathbf{h}_{i n k l} \sqrt{\beta_{i n k l}}\end{aligned}\),       (1)

E1KOBZ_2023_v17n1_185_f0001.png 이미지

Fig. 1. Illustration of the multicell MIMO system with type of devices and pilot reuse factor α = 3 .

where βinkl = qinkl/dinklfis the large-scale fading coefficient between the i-th BS and Unkl which models the path loss and shadow fading, while is the distance between them, f is decay exponent, and qinkl represents a log-normal random variable hinkl ~ 𝒞𝒩(0, IM) is the fast fading vector between the i-th BS and Unkl.

2.2 Random Access

In RA, since the BS does not know which device requests access at this time slot, devices need to start the access by randomly selecting a pilot sequence and send it to the BS. In this process, different device may select the same pilot, which may result in the access fail. To improve the access success probability, we use the pilot reuse to extend the distance of devices selecting the same pilot, so that the pilot interference becomes negligible and the colliding devices from different cells can be detected.

We use the backoff mechanism to reduce the device conflict. The BS requires to calculate the device backoff parameters, and then broadcasts the backoff parameters to devices. The device randomly generates a backoff time according to the received backoff parameters broadcast to postpone the activation. When RA starts, the active device randomly selects a pilot sequence from the orthogonal pilot pool and sends it to the BS. Each pilot sequence consists of τ symbols, and thus the total number of orthogonal pilot sequences is also τ. Assuming the pilot reuse factor is α , then the number of pilots in each cell is τ/α . Fig. 1 shows the system model diagram with α = 3, we call the cells that use the same set of pilots as co-pilot cells. For example, cells with the same color in Fig. 1 are co-pilot cells. We assume that the channel coherence time is T symbols long, and T > τ. In each channel coherence interval, the active device directly transmits the orthogonal pilot along with data to the BS.

We call two devices who choose the same pilot as colliding devices. Colliding devices within the same cell cannot be distinguished by the BS. Hence, their access fail at this RA slot. However, if two colliding devices belong to different cells, even they select the same pilot, due to the pilot reuse, they locate far away from each other, and the BS can distinguish them based on the distinction on received power strength. Therefore, only devices does not conflict with other devices in the same cell can be successfully detected by the BS, and we call them as intra-cell non-colliding devices.

3. Performance Analysis

3.1 Channel Estimation

Due to the backoff mechanism, the k-th type of devices are randomly activated according to the backoff probability pb,k. Let γnkl ∈ {0,1}, where γnkl = 1 denotes that Unkl is active and γnkl = 0 that it is inactive. Let \(\begin{aligned}\sqrt{\tau} \varphi_{n k l} \in \mathbb{C}^{\tau \times 1}\end{aligned}\) be the pilot sequence used by Unkl. If Unkl chooses the same pilot as Un'k'l', we have φnklHφn'k'l' = 1. Otherwise, φnklHφn'k'l' = 0. Then, the M × τ received matrix at the i-th BS can be written as

\(\begin{aligned}\mathbf{Y}_{i}=\sum_{l=1}^{L} \sum_{k=1}^{K} \sum_{n=1}^{N_{k}} \sqrt{\tau P_{n k l}} \mathbf{g}_{i n k l} \boldsymbol{\varphi}_{n k l}^{T} \gamma_{n k l}+\mathbf{W}_{i}\end{aligned}\),      (2)

where Pnkl is the transmit power of Unkl , and Wi is additive white Gaussian noise (AWGN) matrix, in which each element of Wi follows an independent distribution 𝒞𝒩(0,1).

We assume that Umk'i is the targeted device and does not collide with other devices in i-th cell. Then, the i-th BS estimates the channel of Umk'i using Yiφmk'i* as follows

\(\begin{aligned}\mathbf{Y}_{i} \boldsymbol{\varphi}_{m k^{\prime} i}^{*}=\sqrt{\tau P_{m k^{\prime} i}} \mathbf{g}_{m k^{\prime} i}+\sum_{l \in \mathcal{C}_{i}} \sum_{k=1}^{K} \sum_{n=1}^{N_{k}} \sqrt{\tau P_{n k l}} \mathbf{g}_{i n k l} \delta_{n k l}^{\left(m k^{\prime} i\right)}+\boldsymbol{\omega}_{i}\end{aligned}\),       (3)

where 𝒞i is the co-pilot cells of the i-th cell, and δnkl(mk'i) ∈ {0,1} indicates whether Umk'i select the same pilot as Unklnkl(mk'i) = 1) or not (δnkl(mk'i) = 0). Therefore, the estimates of gimk'i can be written as

\(\begin{aligned}\hat{\mathbf{g}}_{i m k^{\prime} i}=\frac{1}{\sqrt{\tau P_{i m k^{\prime} i}}} \mathbf{Y}_{i} \varphi_{m k^{\prime} i}^{*}\end{aligned}\)       (4)

Further, (4) can be represented as (5)

\(\begin{aligned}\hat{\mathbf{g}}_{i m k^{\prime} i}=\mathbf{g}_{i m k^{\prime} i}+\sum_{l \in \mathcal{C}_{i}} \sum_{k=1}^{K} \sum_{n=1}^{N_{k}} \sqrt{\frac{P_{n k l}}{P_{m k^{\prime} i}}} \mathbf{g}_{i n k l} \delta_{n k l}^{\left(m k^{\prime} i\right)}+\frac{1}{\sqrt{\tau P_{i m k^{\prime} i}}} \boldsymbol{\omega}_{i}\end{aligned}\)       (5)

From (5), we can see that the channel estimation of intra-cell non-colliding devices will be interfered not only by Gaussian white noise, but also by devices in the co-pilot cells. Using the estimated channels obtained above, we can then obtain the average spectral efficiency.

3.2 Average Spectral Efficiency

We can get that the received M × 1 signal at the i-th BS is given by (6), where xnkl is the data symbol with 𝔼{|xnkl|2} = 1, and οi is AWGN vector, and oi ~ 𝒞𝒩(0, iM). Note that ri is divided into five independent terms. The first term is the effective signal. The second and third terms indicate interference from other devices in the same cell. The fourth term represents interference from devices in other cells. The fifth term is noise.

\(\begin{aligned} \mathbf{r}_{i}= & \sqrt{P_{m k^{\prime} i}} \mathbf{g}_{i m k^{\prime} i} x_{m k^{\prime} i}+\sum_{n=1, n \neq m}^{N_{k}} \sqrt{P_{n k^{\prime} i}} \mathbf{g}_{i n k^{\prime} i} x_{n k^{\prime} i} \gamma_{n k^{\prime} i} \\ & +\sum_{k=1, k \neq k^{\prime}}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k i}} \mathbf{g}_{i n k i} x_{n k i} \gamma_{n k i}+\sum_{l=1, l \neq i}^{L} \sum_{k=1}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k l}} \mathbf{g}_{i n k l} x_{n k l} \gamma_{n k l}+\mathbf{o}_{i} .\end{aligned}\).       (6)

Since Umk'i does not collide with other devices in i-th cell, we can receive it. We adopt the maximum-ratio-combing (MRC) receiver here and process the data signal, which is given by

\(\begin{aligned} e_{m k^{\prime} i} & =\hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{r}_{i} \\ & =\sqrt{P_{m k^{\prime} i}} \hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i m k^{\prime} i} x_{m k^{\prime} i}^{\prime}+\sum_{n=1, n \neq m}^{N_{k}} \sqrt{P_{n k^{\prime} i}} \hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i n k^{\prime} i} x_{n k^{\prime} \gamma^{\prime}} \gamma_{n k^{\prime} i}+\hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{o}_{i} \\ & +\sum_{k=1, k \neq k^{\prime}}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k i}} \hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i n k i} x_{n k i} \gamma_{n k i}+\sum_{l=1, l \neq i}^{L} \sum_{k=1}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k l}} \hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i n k l} x_{n k l} \gamma_{n k l} .\end{aligned}\).       (7)

Further, (7) can be represented as (8)

\(\begin{aligned} e_{m k^{\prime} i} & =\sqrt{P_{m k^{\prime} i}} \mathbb{E}_{\mathrm{h}}\left\{\hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i m k^{\prime} i}\right\} x_{m k^{\prime} i}+\sum_{n=1, n \neq m}^{N_{k}} \sqrt{P_{n k^{\prime} i}} \hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i n k^{\prime} i} x_{n k^{\prime} i} \gamma_{n k^{\prime} i} \\ & +\sqrt{P_{m k^{\prime} i}}\left[\hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i m k^{\prime} i}-\mathbb{E}_{\mathrm{h}}\left\{\hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i m k^{\prime} i}\right\}\right] x_{m k^{\prime} i}+\hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{o}_{i} \\ & +\sum_{k=1, k \neq k^{\prime}}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k i}} \hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i n k i} x_{n k i} \gamma_{n k i}+\sum_{l=1, l \neq i}^{L} \sum_{k=1}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k l}} \hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i n k l} x_{n k l} \gamma_{n k l},\end{aligned}\),       (8)

where 𝔼h{⦁} is the average of the small-scale fading. We further assume that the channel is ergodic so that an ergodic achievable uplink rate Rmk'i of Umk'i is given b

\(\begin{aligned}R_{m k^{\prime} i}=\log _{2}\left(1+\frac{\sqrt{P_{m k^{\prime} i}}\left|\mathbb{E}_{\mathrm{h}}\left\{\hat{\mathbf{g}}_{i m k^{\prime} i}^{H} \mathbf{g}_{i m k^{\prime} i}\right\}\right|^{2}}{\mathcal{F}_{m k^{\prime} i}}\right)\end{aligned}\),       (9)

where 𝓕mk'i is represented by

\(\begin{aligned} \mathcal{F}_{m k^{\prime} i} & =\sqrt{P_{m k^{\prime} i}} \operatorname{Var}\left\{\hat{\mathbf{g}}_{i m k k^{H} i}^{H} \mathbf{g}_{i m k^{\prime} i}\right\}+\sum_{n=1, n \neq m}^{N_{k}} \sqrt{P_{n k^{\prime} i}} \mathbb{E}_{\mathrm{h}}\left\{\left|\hat{\mathbf{g}}_{i m k k^{\prime} i}^{H} \mathbf{g}_{i n k^{\prime} i}\right|^{2}\right\} \gamma_{n k^{\prime} i} \\ & +\sum_{k=1, k \neq k^{\prime}}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k i}} \mathbb{E}_{\mathrm{h}}\left\{\left|\hat{\mathbf{g}}_{i m k k^{H} i}^{H} \mathbf{g}_{i n k k}\right|^{2}\right\} \gamma_{n k i}+\mathbb{E}_{\mathrm{h}}\left\{\|\left.\hat{\mathbf{g}}_{i m k^{\prime} i}^{H}\right|^{2}\right\} \\ & +\sum_{l=1, l \neq i}^{L} \sum_{k=1}^{K} \sum_{n=1}^{N_{k}} \sqrt{P_{n k l}} \mathbb{E}_{\mathrm{h}}\left\{\left|\hat{\mathbf{g}}_{i m k i}^{H} \mathbf{g}_{i n k l}\right|^{2}\right\} \gamma_{n k l}\end{aligned}\)       (10)

The channel estimation is obtained when Unkl does not collide with other devices in the cell. Let ρnkl ∈ {0,1} , where ρnkl = 1 denotes that Unkl is an intra-cell non-colliding device and ρnkl = 0 that it is not an intra-cell non-colliding device. Then, average spectral efficiency of the k-th type of device in the i-th cell is given by

\(\begin{aligned}S_{k i}=\left(1-\frac{\tau}{T}\right) \mathbb{E}\left\{\sum_{n=1}^{N_{k}} R_{n k i} \rho_{n k i}\right\}\end{aligned}\).       (11)

Since Rnki and ρnki are independent, (11) can be represented as

\(\begin{aligned}S_{k i}=\left(1-\frac{\tau}{T}\right) \sum_{n=1}^{N_{k}} \mathbb{E}\left\{R_{n k i}\right\} \mathbb{E}\left\{\rho_{n k i}\right\}=\left(1-\frac{\tau}{T}\right) \mathcal{R}_{k i} N_{\text {non }, k}\end{aligned}\),       (12)

where 𝓡ki = 𝔼(Rnki} is got by averaging the parameters such as large-scale fading, pilot selection and device activation of Rnki, which eliminates the difference among devices in the same cell and the corresponding subscript is omitted, and Nnon,k = ∑n=1Nk𝔼{ρnki} indicates the average number of intra-cell non-colliding devices of the k-th type device per cell. Hence, the average spectral efficiency per cell is further written as

\(\begin{aligned}S_{i}=\sum_{k=1}^{K} S_{k i}\end{aligned}\).       (13)

Next, we give the expressions for the average number of intra-cell non-colliding devices of the k-th type device per cell.

Lemma 1. The distributions of the indicator γnkl and δnkl(mk'i) are

ℙ[γnkl = 1] = pb,k, ℙ[δnkl(mk'i) = 1|l∈𝒞i] = αpb,k/τ.       (14)

The average number of intra-cell non-colliding devices of the k-th type device per cell can be written as

\(\begin{aligned}N_{\mathrm{non}, k}=\sum_{n=1}^{N_{k}} n\left(\begin{array}{c}N_{K} \\ n\end{array}\right) p_{\mathrm{b}, k}^{n}\left(1-p_{\mathrm{b}, k}\right)^{N_{k}-n}\left(1-\frac{\alpha}{\tau}\right)^{\left(\sum_{k=1}^{K} p_{\mathrm{b}, k} N_{k}-1\right)}\end{aligned}\).       (15)

Proof. See Appendix A:

In the above results, we can see that Nnon,k is related to τ and the parameters of devices, and it is not only affected by its own activation probability and number, but also by the activation probability and number of other types in the cell. When the system model reduce to a single cell with one type of device, (15) becomes \(\begin{aligned}N_{\text {no-colli }}=N_{\text {tatal }} p_{b}\left(1-\frac{1}{\tau}\right)^{N_{\text {tatal }} p_{b}-1}\end{aligned}\), where Ntotal is the total number of devices, and pb is the backoff probability. We can find that Nno-colli is the same as (1) in [18].

We use the power control technique to compensate for the large-scale fading of devices in each cell. The power control factor is given by λ = Pnklβinkl (for any n, k, and l). Then, with it we can substitute Nnon,k from (15) into (12) to get Si.

Theorem 1. The approximation of the average spectral efficiency of the i-th cell is given by

\(\begin{aligned}S_{i}=\sum_{k=1}^{K} S_{k i}\end{aligned}\),       (16)

where Ski is as follows

\(\begin{aligned}S_{k i} \approx\left(1-\frac{\tau}{T}\right) N_{\mathrm{non}, k} \log _{2}\left(1+\frac{\tau M \lambda^{2}}{\mathcal{Q}}\right)\end{aligned}\),       (17)

in which 𝒬 is given by (18). In (18), 𝒫 denotes the average transmitted power of all devices,

\(\begin{aligned} \mathcal{Q} & =(\tau \lambda+1) \lambda\left[p_{\mathrm{b}, k^{\prime}}\left(N_{k^{\prime}}-1\right)+\sum_{k=1, k \neq k^{\prime}}^{K} p_{\mathrm{b}, k} N_{k}\right]+(\tau \lambda+1)(\lambda+1) \\ & +\mathcal{P}^{2} \alpha\left(\sum_{k=1}^{K} p_{\mathrm{b}, k} N_{k}\right)\left[1+p_{\mathrm{b}, k^{\prime}}\left(N_{k^{\prime}}-1\right)+\sum_{k=1, k \neq k^{\prime}}^{K} p_{\mathrm{b}, k} N_{k}\right] \sum_{l \in \mathcal{C}_{i} v_{i k l}} \varepsilon_{i k l}^{2} \\ & +\mathcal{P}^{2} \alpha\left(\sum_{k=1}^{K} p_{\mathrm{b}, k} N_{k}\right)^{2} \sum_{l_{1}=1, l_{1} \neq i l_{2} \in \mathcal{C}_{l} \backslash, l_{2} \neq l_{1}}^{L} \varepsilon_{i k l_{1}} \varepsilon_{i k l_{2}}+\mathcal{P}^{2} M \alpha\left(\sum_{k=1}^{K} p_{\mathrm{b}, k} N_{k}\right) \sum_{l \in \mathcal{C}_{i}, i} \mu_{i k l} \\ & +\mathcal{P}\left(\sum_{k=1}^{K} p_{\mathrm{b}, k} N_{k}\right)\left\{\alpha \tau+\alpha+1+\tau \lambda+\alpha \lambda\left[p_{\mathrm{b}, k^{\prime}}\left(N_{k^{\prime}}-1\right)+\sum_{k=1, k \neq k^{\prime}}^{K} p_{\mathrm{b}, k} N_{k}\right]\right\} \sum_{l \in \mathcal{C}_{i} \mid i} \varepsilon_{i k l} .\end{aligned}\)       (18)

while εikl = 𝔼{βinkl} and μikl = 𝔼{βinkl2} represent the average and average square of largescale fading from a device of the k-th type in the l-th cell to the i-th BS.

Proof. See Appendix B:

The above results are averaged for all randomness, such as device activation, and pilot selection. That is, strong versatility, suitable for most networks. Obviously, Si is not a monotone function of pb,k, so we need to find the optimal value for pb,k. Next, we look for the optimal pb,k to maximize (16).

3.3 Optimal backoff parameters

In this section, we aim to find the optimal pb,k for each type of device that maximize Si. That is, to solve the following optimization problem.

\(\begin{aligned}(P 1) \quad \max _{p_{\mathrm{b}, 1}^{*}, \ldots, p_{\mathrm{b}, K}^{*}} \quad S_{i}\\ \begin{array}{ll}\text { s.t. } & 0<p_{\mathrm{b}, k} \leq 1, \forall k \\ & 0<\mathcal{D}_{k} \leq \mathcal{D}_{k}^{\mathrm{u}}, \forall k\end{array}\\\end{aligned}\).       (19)

Because of the backoff mechanism, it takes time for each user to access, and the time required for user access is called access delay. We use 𝒟k to denote the delay for the k-th type of device. Considering that devices in IoT usually short time requirements, we set delay constraints for the k-th type of devices as 𝒟ku, i.e., the maximum permissible delay for the k-th type of devices. Note that the unit of delay constraints is time slot and will not be specified for simplication. From an average point of view, 𝒟k= 1/pb,k. Therefore, (P1) can be simplified into (P2).

\(\begin{aligned}\begin{array}{l}(P 2) \quad \max _{p_{b, 1,1} \ldots, p_{b, K}} \quad S_{i} \\ \text { s.t. } \quad \frac{1}{\mathcal{D}_{k}^{\mathrm{u}}}<p_{\mathrm{b}, k} \leq 1, \forall k \text {. } \\\end{array}\end{aligned}\).       (20)

For (P2), we use gradient descent. pb,k is iterated continuously until Si converges, and the specific algorithm details can be summarized in algorithm 1, as shown in Table 2.

Table 2. Optimization algorithm

E1KOBZ_2023_v17n1_185_t0002.png 이미지

We also analyze the complexity of algorithm 1. We can see that the complexity of algorithm 1 depends on the gradient of the backoff probability. The main overall complexity is 𝑂(tKNk + tK + tKϑ + 2tK2ϑ + tK2Lϑ), where ϑ is the number of co-pilot cells in the i-th cell [25]. Therefore, we can see that the algorithm complexity mainly depends on the number of device types, the number of devices and the number of co-pilot cells.

4. Numerical Results

In this section, we consider two types of devices exists in the network, denoted as device 1 and 2, respectively. Devices are randomly and uniformly distributed in a cell with a radius rc = 1000 of meters, and there is an area with a radius of rh = 100 meters in the middle where no devices are distributed. The large-scale fading model is expressed as βinkl = zinkl/(rinkl/rh)υ, where zinkl is a log-normal random variable with standard deviation σ, υ is the path loss exponent, and rinkl is the distance between Unkl and the i-th BS. We assume that σ = 8 dB and υ = 3.8. εikl and μikl are got by averaging over 1000 realizations of devices spatial distributions in each cell. We assume that the data is transmitted using an OFDM system and T = 196 symbols. We set λ = 0.1 and α = 3.

4.1 Performance Validation

In Fig. 2, we compare the simulated spectral efficiency in (11) and the analytic approximation in (16). Clearly, in all cases, no matter with a small or large τ, the analytical value and the simulation value are very tight, which indicates the accuracy and rigor of our analysis results. Given the tightness, we use the analysis results of (16) for the following investigations.

E1KOBZ_2023_v17n1_185_f0002.png 이미지

Fig. 2. Average spectral efficiency per cell vs. BS antenna number, where pb,1 = pb,2 = 0.5.

4.2 Performance Analysis

In Fig. 3, we study the relationship between the backoff probability of devices and the average spectral efficiency. We use pb,1 and pb,2 to represent the backoff probability of device 1 and device 2 , respectively. Firstly, we find that with the increase of pb,1, the average spectral efficiency Si keeps growing until it reaches a critical point, and gradually decrease. This is because when backoff probability is small, there are few active devices connected to the system. Therefore, pilots are sufficient for devices and Si will increase with growth of active devices, i.e., the growth of backoff probability. However, when backoff probability increases to a critical point, the number of active devices exceeds the number of pilots. Therefore, the pilot collision is aggravated when active devices grows and Si begin to reduce. Given that, we find that there exists an optimal backoff probability that maximizes Si. Secondly, it can be seen from Fig. 3 that when τ is fixed, the optimal backoff probability decreases with the increment of the device number. Conversely, when the number of devices is fixed, the optimal backoff probability decreases with the decrement of τ. Next, we discuss the variation of the optimal backoff probability with respect to devices number and delay constraints.

E1KOBZ_2023_v17n1_185_f0003.png 이미지

Fig. 3. Average spectral efficiency per cell vs. backoff probability, where M = 200, and pb,1 = pb,2.

Fig. 4 shows the optimal backoff probability with and without delay constraints, where the number of device 2 is fixed (N2 = 140 ). Fig. 4. (a) is the case without delay constraints, and it can be seen that when N1 < N2, pb,1 > pb,2 ; and when N1 > N2, pb,1 < pb,2. This indicates that devices with larger quantity should have more backoff time. When the total number of devices in a cell is large, the type of devices that have larger quantity backoff altogether, that is, the access opportunity should be completely transferred to devices with smaller quantity. Fig. 4. (b) is the case with delay constraints. It can be seen that due to delay constraints, devices with larger quantity cannot completely backoff, and their optimal backoff probability is the lower bound determined by delay constraints. Moreover, when there is a large number of devices in the cell, the lower bound is taken for both devices with larger quantity and devices with smaller quantity.

E1KOBZ_2023_v17n1_185_f0004.png 이미지

Fig. 4. Optimal backoff probability vs. number of device 1, where N2 = 140, and M = 200.

Now, we discuss the influence of delay constraints on the optimal backoff probability. Fig. 5 shows the change of optimal backoff probability of two devices with different numbers when 𝒟2u is fixed (𝒟2u = 10 ). It can be seen that with the increment of 𝒟1u, the optimal pb,1 decreases and is always at the lower limit of its feasible region due to delay constraints, while the optimal pb,2 first lies a the lower bound and then grow up to exceed the pb,1 when 𝒟1u increases. This indicates that a larger number of devices should be allocated more backoff time when delay constraints are relatively loose. However, as delay constraints become stricter, the required backoff time of this type of devices reduces, and may even be less than devices with smaller quantity in the case with extremely strict delay constraints.

E1KOBZ_2023_v17n1_185_f0005.png 이미지

Fig. 5. Optimal backoff probability vs. delay constraints of device 1, where 𝒟2u = 10, and M = 200.

Fig. 6 shows how the optimal backoff probability of the two types of devices varies with respect to delay constraints where 𝒟1u = 𝒟2u. We can see that when the number of device remains unchanged and the number of device increases from 60 to 130, the optimal pb,1 and optimal pb,2 are the same as the lower limit of the feasible domain under strict delay constraints, which confirms the conclusion of Fig. 5. Moreover, we compare the optimal backoff probability under and τ = 90 and τ = 120, and find that in the case with large τ, the optimal backoff probability has a stronger tolerance to delay constraints, and the system can accommodate more devices to access. When τ is small, the optimal backoff probability take the lower limit of the feasible region more often than that with large τ, which indicates that the effect of delay constraints works more obviously when the pilot length is short.

E1KOBZ_2023_v17n1_185_f0006.png 이미지

Fig. 6. Optimal backoff probability vs. delay constraints of device 1, where 𝒟1u = 𝒟2u, and M = 200.

In Fig. 7, we compare the average spectral efficiency Si with full access, i.e., pb,k = 1, with the optimal backoff mechanism we proposed. It can be seen that when N1 increases gradually, Si increases first and then tends to be stable with optimal pb,k, but Si with pb,k = 1 increases first and then decreases. It can be seen that the system with optimal backoff mechanism is obviously better than the system without backoff mechanism. We also find that with the increment of the number of devices, the advantage of the optimal backoff mechanism becomes more obvious, i.e., the gain of Si also increases. Therefore, in large user scenarios, which is a typical scenario in IoT, it is necessary to adopt the optimal backoff mechanism we proposed.

E1KOBZ_2023_v17n1_185_f0007.png 이미지

Fig. 7. Average spectral efficiency per cell vs. number of device 1, where N1 = N2, M = 200.

5. Conclusion

This paper has investigated the GFRA in multicell massive MIMO systems with mixed-type devices with delay constraints. The analytic approximation of the average spectral efficiency is derived by using indicators which captured devices' behaviors uncertainty and verified by simulation. we obtain the optimal backoff parameter for each type of devices under their delay constraints. It is found that the optimal backoff parameters are closely related to the devices number and delay constraints. In general, devices that have larger quantity should have more backoff time before they are allowed to access. However, as delay constraints become stricter, the required backoff time of devices reduces gradually, and the device with larger quantity may have less backoff time than that with smaller quantity when its delay constraint is extremely strict. Moreover, when the pilot length is short, the effect of delay constraints mentioned above works more obviously. In conclusion, the system with backoff mechanism is more practical and has better performance. Especially in the scenario of IoT, it is necessary to adopt the optimal backoff mechanism proposed by us.

Appendix

Acknowledgement

This work was supported in part by the Graduate Research and Innovation Projects of Jiangsu Province under SJCX21_0263. The work of Qi Zhang was supported in part by the National Natural Science Foundation of China under Grant 62171231 and 61801244, in part by the National Key Research & Development Program of China under Grant No. 2020YFB1807202.

References

  1. J. Li, M. Wen and T. Zhang, "Group-Based Authentication and Key Agreement With Dynamic Policy Updating for MTC in LTE-A Networks," IEEE Internet Things J., vol. 3, no. 3, pp. 408-417, Jun. 2016. https://doi.org/10.1109/JIOT.2015.2495321
  2. X. Liu, T. Huang, N. Shlezinger, Y. Liu, J. Zhou and Y. C. Eldar, "Joint transmit beamforming for multiuser MIMO communications and MIMO radar," IEEE Trans. Signal Process., vol. 68, pp. 3929-3944, 2020. https://doi.org/10.1109/tsp.2020.3004739
  3. H. Q. Ngo, E. G. Larsson, and T. L. Marzetta, "Energy and spectral efficiency of very large multiuser MIMO systems," IEEE Trans. Commun., vol. 61, no. 4, pp. 1436-1449, Apr. 2013. https://doi.org/10.1109/TCOMM.2013.020413.110848
  4. M. Y. Ali, T. Hossain, and M. M. Mowla, "A trade-off between energy and spectral efficiency in massive MIMO 5G system," in Proc. of ICECTE, Rajshahi, Bangladesh, pp. 209-212, 2019.
  5. O. Mahmoud and A. El-Mahdy, "Performance evaluation of noncoherent DPSK signal detection algorithms in massive MIMO systems," in Proc. of Conf. TELE, Leiria, Portugal, pp. 1-5, 2021.
  6. T. Wang, Y. Wang, C. Wang, Z. Yang and J. Cheng, "Group-based random access and data transmission scheme for massive MTC networks," IEEE Trans. Commun., vol. 69, no. 12, pp. 8287-8303, Dec. 2021. https://doi.org/10.1109/TCOMM.2021.3111609
  7. J. Choi, "Fast retrial for low-latency connectivity in MTC with two different types of devices," IEEE Wireless Commun. Lett., vol. 9, no. 10, pp. 1786-1789, Oct. 2020. https://doi.org/10.1109/lwc.2020.3004822
  8. Z. Chen, F. Sohrabi, and W. Yu, "Multi-cell sparse activity detection for massive random access: Massive MIMO versus cooperative MIMO," IEEE Trans. Wireless Commun., vol. 18, no. 8, pp. 4060-4074, Aug. 2019. https://doi.org/10.1109/twc.2019.2920823
  9. X. Fang, M. Zhao, S. Zhang, W. Zhou, and J. Zhu, "Grouping-based grant-free random access based on statistical distribution of user rates," in Proc. of ICC, Shanghai, China, pp. 1-6, 2019.
  10. H. Han, L. Fang, W. Lu, K. Chi, W. Zhai and J. Zhao, "A novel grant-based pilot access scheme for crowded massive MIMO systems," IEEE Trans. Veh. Technol., vol. 70, no. 10, pp. 11111-11115, Oct. 2021. https://doi.org/10.1109/TVT.2021.3111627
  11. O. S. Nishimura, J. C. Marinello and T. Abrao, "A grant-based random access protocol in extralarge massive MIMO system," IEEE Wireless Commun. Lett, vol. 24, no. 11, pp. 2478-2482, Nov. 2020. https://doi.org/10.1109/LCOMM.2020.3012586
  12. A. T. Abebe and C. G. Kang, "Multi-sequence spreading random access (MSRA) for compressive sensing-based grant-free communication," IEEE Trans. Commun., vol. 69, no. 11, pp. 7531-7543, Nov. 2021. https://doi.org/10.1109/TCOMM.2021.3103542
  13. Technical specification group services and system aspects, 2-step RACH for NR, document 3GPP TR 21.916 (V16.2.0), 3GPP, June. 2022. [Online]. Available: https://www.3gpp.org/ftp/Specs/archive/21_series/21.916
  14. J. Ding, M. Feng, M. Nemati, and J. Choi, "Performance analysis of massive MIMO assisted semi-grant-free random access," in Proc. of CCNC, Las Vegas, NV, USA, pp. 1-7, 2021.
  15. Y. Liu, M. Zhao, L. Xiao and S. Zhou, "Pilot domain NOMA for grant-free massive random access in massive MIMO marine communication system," China Commun., vol. 17, no. 6, pp. 131-144, June 2020. https://doi.org/10.23919/jcc.2020.06.011
  16. Z. Zhang, Y. Li, C. Huang, Q. Guo, L. Liu, C. Yuen, and Y. L. Guan, "User activity detection and channel estimation for grant-free random access in LEO satellite-enabled internet of things," IEEE Internet Things J., vol. 7, no. 9, pp. 8811-8825, Sept. 2020. https://doi.org/10.1109/jiot.2020.2997336
  17. J. Ding and J. Choi, "Preamble-data superposition random access in massive MIMO systems," IEEE Wireless Commun. Lett., vol. 9, no. 6, pp. 906-910, June 2020. https://doi.org/10.1109/lwc.2020.2974724
  18. H. Jiang, D. Qu, J. Ding, and T. Jiang, "Multiple preambles for high success rate of grant-free random access with massive MIMO," IEEE Trans. Wireless Commun., vol. 18, no. 10, pp. 4779-4789, Oct. 2019. https://doi.org/10.1109/twc.2019.2929126
  19. J. Ding, D. Qu, and H. Jiang, "Optimal preamble length for spectral efficiency in grant-free RA with massive MIMO," in Proc. of ICEIC, Auckland, New Zealand, pp. 1-5, 2019.
  20. H. Wang, X. Zhang, and S. Zhou, "Analysis of system capacity of grant-free access with massive MIMO," in Proc. of ICCC Workshops, Xiamen, China, pp. 410-415, 2021.
  21. Z. Chen, F. Sohrabi, and W. Yu, "Multi-cell sparse activity detection for massive random access: Massive MIMO versus cooperative MIMO," IEEE Trans. Wireless Commun., vol. 18, no. pp. 4060-4074, Aug. 2019. https://doi.org/10.1109/twc.2019.2920823
  22. Medium Access Control (MAC) protocol specification, backoff parameter values, document 3GPP TS 36.321 (V10.3.0), 3GPP, Sept. 2011. [Online]. Available: https://www.3gpp.org/ftp/Specs/archive/36_series/36.321
  23. P. Jiang, H. Wang and M. De Mari, "Optimal dynamic backoff for grant-free NOMA IoT networks: a mean field game approach," in Proc. of ICCC, Sanshui, Foshan, China, pp. 997-1002, 2022.
  24. J. Ding, D. Qu, H. Jiang and T. Jiang, "Success probability of grant-free random access with massive MIMO," IEEE Internet of Things Journal, vol. 6, no. 1, pp. 506-516, Feb. 2019. https://doi.org/10.1109/jiot.2018.2869003
  25. G. Lei, Y. Liu and X. Xiao, "Analysis of joint transmit and receive antenna selection in CPM MIMO systems" KSII Transactions on Internet and Information Systems, vol. 11, no. 3, pp. 1425-1440, Mar. 2017. https://doi.org/10.3837/tiis.2017.03.010