Next Article in Journal
Hidden Variable Models in Text Classification and Sentiment Analysis
Previous Article in Journal
Fuzzy PID Control Design of Mining Electric Locomotive Based on Permanent Magnet Synchronous Motor
Previous Article in Special Issue
On the Quasistationarity of the Ambient Electromagnetic Field Generated by Wi-Fi Sources
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Outdoor-to-Indoor mmWave Relaying with Massive MIMO: Impact of Imperfect Channel Estimation †

by
Nawal Bounouader
*,
Houda Chafnaji
and
Mustapha Benjillali
Communication Systems Department, INPT (Institut National des Postes et Télécommunications), Rabat 10100, Morocco
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in 6th International Conference on Advanced Communication Technologies and Networking (CommNet), Rabat, Morocco, 11–13 December 2023.
Electronics 2024, 13(10), 1857; https://doi.org/10.3390/electronics13101857
Submission received: 4 April 2024 / Revised: 6 May 2024 / Accepted: 8 May 2024 / Published: 10 May 2024
(This article belongs to the Special Issue New Trends and Methods in Communication Systems, 2nd Edition)

Abstract

:
Assuming incomplete knowledge of the channel state information (CSI), we investigate two scenarios involving millimeter wave (mmWave) relaying to support outdoor-to-indoor communications. We proceed to derive the average signal-to-noise ratio (SNR) expressions for two relaying scenarios and quantify the asymptotic SNR. The performance of the two relaying scenarios is evaluated using the outage probability—for which we have derived closed-form equations—the end-to-end channel capacity, and the energy efficiency. The obtained results are compared with those derived assuming complete knowledge of the CSI. The effect of the imperfect CSI is therefore assessed in relation to the reference of perfect CSI. In these scenarios, an outside base station (BS) in an urban cellular network serves several indoor users. In the context of a two-hop full-duplex (FD) relaying scheme, we initially suggest a method in which the base station (BS) utilizes zero-forcing (ZF) precoding, and we take into account the overall channel response. Furthermore, we make the assumption that the base station (BS) engages in precoding only depending on the response of the channel in the first hop; in this second design, the relay precodes (using the response of the second-hop channel), amplifies, and sends the signals. Both techniques utilize massive multiple-input–multiple-output (mMIMO) arrays to permit transmission. We also present Monte Carlo simulation results to assess the accuracy of our analytical results. Finally, the two systems are compared in terms of channel estimation and precoding complexity, the number of antennas, as well as the number of users. Practical deployment recommendations are formulated at the end of this work.

1. Introduction

A preliminary version of this work was published in [1], and the results obtained previously will be considered as a reference to compare with the new results. Hence, the present work presents an extended version wherein the important and more practical impact of channel estimation is taken into consideration.
The millimeter wave (mmWave) spectrum is an essential element of fifth-generation (5G) and future (6G) wireless communications, as stated in [2]. The mmWave technology offers the benefit of increased capacity [3] and more focused beams [4,5], but it comes with the drawback of greater loss of signal when penetrating buildings. The aforementioned disadvantage is often reduced by employing relays to connect indoor users to outdoor transceivers, along with the implementation of scalable massive multiple-input–multiple-output (mMIMO) configurations [6,7] to enhance signal coverage [8,9], utilizing nearly optimal linear processing techniques (such as zero-forcing (ZF) and maximum-ratio (MR)), and reducing channel estimation complexity (using achieved channel hardening, among other methods). Many mMIMO detection algorithms have been proposed in the literature; ref. [10] provides a survey of those algorithms.
However, full-duplex (FD) relaying, as described in [11], involves nodes that may broadcast and receive at the same time and on the same frequency resource. This is a significant advancement towards fulfilling the spectrum efficiency and system capacity needs of future wireless networks. Although there are still underlying concerns that need to be addressed [12], the directional characteristics of mMIMO (massive multiple-input–multiple-output) technology operating in the millimeter wave frequency range help to mitigate the challenges associated with implementing full-duplex (FD) relaying.
Power considerations are of utmost importance in this context, particularly because precoding is necessary to prepare the signals for the mMIMO broadcasts. To account for the power limitations at the base station (BS) and the relay node, it is necessary to normalize the precoding matrix by adjusting its power [11,13]. In [14], the authors present a mathematical formula for calculating the trace of the zero-forcing (ZF) precoding matrix at the base station. On the other hand, [15] provides the expression of the normalization coefficient at the relay node while adopting different relaying schemes.
There are two main approaches to normalization: vector normalization (VN) and matrix normalization (MN). The authors in [16] present a comparative approach that examines the differences between VN and MN for various linear precoding techniques. Channel estimation is a crucial concern that has been addressed in numerous studies exploring the problem of mMIMO-based mmWave relaying (e.g., [17,18]). According to [14], the complexity of channel estimation can be decreased in a mMIMO system due to the occurrence of channel hardening. Within this specific situation, it is only necessary to use uplink pilot symbols to calculate the channel using reciprocity; in this case, the residual channel estimation error can be disregarded when determining the power normalization coefficient [19].
This paper examines two outdoor-to-indoor millimeter wave (mmWave) full-duplex (FD) relaying strategies considering both perfect [1] and imperfect channels. The performance of these strategies is evaluated based on outage probability, capacity, and energy efficiency while also considering the complexity of channel estimation and precoding for each scheme. In the first design, the base station (BS) employs zero-forcing (ZF) precoding: utilizing the two-hop channel response from end to end. In this scenario, the relay performs the functions of normalizing, amplifying, and forwarding the signals to the indoor users. The base station (BS) predicts the comparable uplink channel, and the downlink channel state information (CSI) is derived by leveraging the principle of channel reciprocity [20].
In contrast, the second technique involves the BS performing precoding only based on the channel knowledge of the first hop. In this second scheme, the relay performs the functions of normalizing, precoding (using the second-hop channel response), amplifying, and forwarding the signals. The relay estimates the second-hop channel, while the first-hop channel is estimated at the BS. The power normalization coefficients are gradually approached as a limit, enabling us to calculate the average and asymptotic signal-to-noise ratios (SNRs) in the presence of both perfect and imperfect signal channel information (SCI). The outage probability expressions are derived for each precoding scheme.
Furthermore, we offer estimations of the level of complexity involved in channel estimation and precoding. We also deduce a relationship between the number of antennas in each scheme to guarantee equal complexity. Therefore, the subsequent performance comparison framework will consider two cases: (i) analyzing performance without taking complexity into account and (ii) choosing different numbers of antennas in both schemes to ensure equal complexity and enable a fair performance comparison.
The rest of the paper is structured as follows. The global system model is introduced in Section 2. Section 3 discusses and analyzes the first scheme, whereas Section 4 provides the same analysis for the second scheme. Section 5 displays the numerical and simulation results that were obtained; it also highlights the performance that was obtained when taking into account different system parameters and constraints. The paper is concluded in Section 6.

Mathematical Notation

In the rest of the paper, we use the following notational conventions: ( · ) T and ( · ) H , respectively, for the transpose and the Hermitian transpose, E [ · ] for the mathematical expectation, · for vectors norm, tr [ · ] for the trace of a matrix, and ( · ) + for the Moore–Penrose pseudo-inverse of a matrix.

2. System and Signal Models

As illustrated in Figure 1, we consider N u single-antenna users in a dual-hop downlink system in communication with a BS B equipped with N B antennas via an FD relay node R equipped with N R 1 receive antennas and N R 2 transmit antennas.
At time instant i, R receives a data vector y RR ( i ) and transmits the previously received data, which has been processed, normalized, and amplified y RT ( i 1 ) . To alleviate the notation, we drop the time index in the remainder of the analysis. Let us note that H 1 C ( N R 1 , N B ) and H 2 C ( N u , N R 2 ) are the channel matrix between the BS and the relay receiver and the channel matrix between the relay transmitter and the users, respectively. Both nodes B and R transmit with large numbers of antennas. Assuming that all the paths between the BS and the receiver side of R have the same large-scale fading statistic β 1 , we denote H 1 = β 1 Z 1 , and since a Rayleigh fading channel model is considered in this work, the elements of Z 1 are independent circularly symmetric Gaussian with zero mean and unit-normalized variance. The second hop is modeled similarly, where H 2 denotes the channel matrix of the second hop, and all paths between the relay transmitter and the individual users have the same large-scale fading statistic β 2 , i.e., H 2 = β 2 Z 2 , where the elements of Z 2 are independent circularly symmetric Gaussian with zero mean and unit-normalized variance. We denote by q the N u × 1 vector of user symbols, with E [ q q H ] = I N u . The transmitted vector x is then given by
x = P F 1 q ,
where P is the source power budget for each user, and F 1 is the precoding matrix performed such that E [ x 2 ] = N u P . The received signal at R is given by:
y RR = H 1 x + H RR y RT + n 1 ,
where H RR C ( N R 1 , N R 2 ) is the matrix representing the self-interference channel (taking into consideration the advantages of mMIMO antenna array directivity, we neglect the self-interference term in the remainder of the analysis [21]), y RT is the data vector transmitted by R during the same time interval (containing data received in the previous time interval, which has been processed, normalized, and amplified), and n 1 represents zero-mean additive white Gaussian noise (AWGN) at the relay’s receiver side with variance σ 1 2 .
The processing at the relay consists of a linear signal processing unit, which is denoted by the relaying matrix F 2 and results in a processed version expressed as
y RP = F 2 y RR = P F 2 H 1 F 1 q + F 2 n 1 .
The relay then normalizes and amplifies the signal before forwarding it to the indoor users. Let α denote the normalization coefficient (in the remainder of the text, to avoid confusion, we denote the normalization coefficients by α 1 and α 2 when analyzing the first and second, relaying scheme, respectively) given by
α = tr y RP y RP H = tr F 2 ( P H 1 F 1 F 1 H H 1 H + n 1 n 1 H ) F 2 H .
With a power budget of P R at the relay, the transmitted signal is then expressed as
y RT = P R α y RP .
Finally, user k receives
y ( u , k ) = P R P α h ( 2 , k ) F 2 H 1 F 1 q + P R α h ( 2 , k ) F 2 n 1 + h ( D , k ) F 1 q + n ( 2 , k ) ,
where n ( 2 , k ) represents the zero-mean AWGN with variance σ 2 2 at user k, h ( 2 , k ) and h ( D , k ) denote, respectively, the k-th lines of matrices H 2 and H D , and H D C ( N u , N B ) represents the direct link channel matrix. Due to the blockage of mmWave signals by the outdoor-to-indoor separation, the direct link is subsequently neglected [22].

3. All-Pass Relay with Zero-Forcing Precoding at BS

3.1. Case of Perfect CSI

In this case, the relay is simply designed as an all-pass amplify-and-forward (AF) unit, i.e., N R 2 = N R 1 = N R and F 2 = I N R , while the BS performs precoding based on the end-to-end channel. Let us denote H eq = H 2 × H 1 = N R β 1 β 2 Z eq . According to the central limit theorem, the elements of Z eq are complex Gaussian with zero mean and unit variance. As in [14], the ZF beamforming matrix is given by F 1 = N B N u Z eq + , where Z eq + = Z eq H ( Z eq Z eq H ) 1 .
To characterize the complexity of this first scheme, we consider the complexity analysis framework presented in [23] wherein only multiplications are taken into account. Therefore, the precoding complexity of Scheme 1 can be directly approximated as O ( 2 N B N u 2 + N u 3 ) .

3.1.1. Asymptotic SNR Analysis

In the following, we first give the asymptotic value of α in the large number of antennas regime; then, we derive the instantaneous, average, and asymptotic SNR expressions.
With Scheme 1, the normalization coefficient is expressed as α 1 = tr [ ( H 1 x + n 1 ) ( H 1 x + n 1 ) H ] , which, after a few manipulations, yields
α 1 = P ( N B N u ) · tr H 1 Z eq H ( Z eq Z eq H ) 1 q q H ( ( Z eq Z eq H ) 1 ) H Z eq H 1 H + tr n 1 n 1 H .
Based on [24] (Lemma 1), when N R = N B = N , we obtain
H 1 H 1 H N N β 1 I N , H 2 H 2 H N N β 2 I N u , Z eq Z eq H N N I N u , and n 1 n 1 H N N σ 1 2 I N .
Therefore, the normalization coefficient α 1 asymptotically goes to
α 1 N N P β 1 N u + σ 1 2 ,
and since the signal received by user k is expressed as
y ( u , k ) = P ( N B N u ) P R N R β 1 β 2 α 1 q k + P R α 1 h ( 2 , k ) n 1 + n ( 2 , k ) ,
the end-to-end instantaneous SNR can be expressed as
Γ 1 inst = P R P N R ( N B N u ) β 1 β 2 | q k | 2 P R | h ( 2 , k ) n 1 | 2 + α 1 | n ( 2 , k ) | 2 ;
the average SNR can be expressed as
Γ 1 avg = P R P N R ( N B N u ) β 1 β 2 P R N R β 2 σ 1 2 + α 1 σ 2 2 ,
and, in the large number of antennas regime, the end-to-end SNR at user k can be expressed as
Γ 1 asym = P R P β 1 β 2 P R β 2 σ 1 2 N + P β 1 N u σ 2 2 N + σ 1 2 σ 2 2 N .
Equation (12) shows that the asymptotic SNR for Scheme 1 goes to infinity as N . This means that the undesirable effects from relay and user noise disappear when N becomes large, and only the useful signal is dominant.
It is well known that one advantage of using mMIMO is that power scaling can improve the energy efficiency while maintaining a desired capacity due to the large diversity gain of the large antenna array. Therefore, we propose and analyze a few typical power scaling options for Scheme 1. The asymptotic SNRs with power scaling ( P = E / N a , P R = E R / N b , 0 a , b 1 ) are summarized in Table 1. Note that E and E R are fixed regardless of N. Finally, we can straightforwardly express the system’s capacity and asymptotic capacity, respectively, as
C 1 = log 2 ( 1 + N u Γ 1 avg ) ,
C 1 asym = log 2 ( 1 + N u Γ 1 asym ) .

3.1.2. Outage Probability Analysis

In this subsection, we derive the outage probability of Scheme 1 using (10) and (8). The instantaneous SNR can be rewritten as
Γ 1 inst = A λ 1 X + λ 2 Y ,
where A = P R P N R ( N B N u ) β 1 β 2 , λ 1 = P R N R β 2 σ 1 2 / 2 , λ 2 = α 1 σ 2 2 / 2 , and X and Y are two independent chi-squared distributed variates with two degrees of freedom [25]. Note that λ 1 X and λ 2 Y are gamma variables with shape parameter 1 and respective scale parameters a 1 = 1 / 2 λ 1 and a 2 = 1 / 2 λ 2 . Then, by following the same approach as in [25], we obtain the cumulative distribution function (CDF) of the approximate SNR in Scheme 1 as
P 1 out ( Γ th ) = a 1 exp a 2 A Γ th a 2 exp a 1 A Γ th a 1 a 2 .
This yields the outage probability expression that we will evaluate numerically later in Section 5.

3.2. Case of Imperfect CSI

In this subsection, we consider the case of imperfect CSI, for which the channels are estimated using a minimum mean square error (MMSE) approach, and we re-derive the SNR and the outage probability expressions obtained in the previous subsection. For Scheme 1, we assume that the relay amplifies and forwards the received pilots, and only the BS estimates the equivalent end-to-end channel.

3.2.1. Pilot Symbol Transmission

We denote by τ c the number of samples that can be sent during a channel coherence interval, where τ c = B · T , with T being the duration of the coherence interval in which the channel is considered to remain time-invariant, and B is the bandwidth of the waveform, for which the frequency response of the channel is considered to be flat. Each coherence interval hosts N u orthogonal pilot waveforms of length τ P , where N u τ P τ c [14]. Note that the BS estimates the equivalent end-to-end channel for each user based on the received uplink pilots.
Let us denote by Φ a τ P × N u unitary matrix that contains the pilot symbols, such that
Φ H Φ = I N u ,
and all the users simultaneously transmit N u signals of duration τ P and of the form [14]
X P = τ P Φ H .
The relay receives the pilot symbols and amplifies and retransmits the signal to the BS, where the received signal can be expressed as
Y PR = H 2 T X R + N 3 ,
where N 3 represents the zero-mean AWGN matrix at the relay, the elements of which all have a variance σ 3 2 .
Finally, let us denote by α P = tr ( Y PR Y PR H ) the power normalization coefficient, and the relay hence transmits
Y PT = 1 / α P Y PR .
We can note that α P has the same structure as α 1 and can be approximated as:
α P N N τ P β 2 N u + σ 3 2 ,

3.2.2. Despreading at the BS

During this phase, the BS uses the pilot symbol matrix Φ for despreading the signal. Since the received pilot signal at the BS is given by
Y PB = H 1 T Y PT + N 4 ,
where N 4 represents the zero-mean AWGN matrix at the BS, for which the elements all have a variance σ 4 2 , this yields a received signal with the form
Y P = Y PB Φ = τ P / α P H eq T + ( 1 / α P H 1 T N 3 + N 4 ) Φ .

3.2.3. MMSE Channel Estimation

Let us denote by N a ( 1 / α P H 1 T N 3 + N 4 ) Φ the aggregate noise. The elements of N a are independent and identically distributed (i.i.d.) circularly symmetric Gaussian variables, i.e., N a ( i , j ) CN ( 0 , N R β 1 σ 3 2 / α P + σ 4 2 ) , where i and j denote the i-th line and j-th column, respectively, of N a ( i , j ) . Referring to [26], the estimated channel between the BS and the users is thus given by
H ^ eq T = τ P α P β e Δ Y P ,
with Δ = ( τ P β e + N R β 1 σ 3 2 ) / α P + σ 4 2 , and β e = N R β 1 β 2 . It is worth noting that the variance of the elements of H ^ eq T can be expressed as
ζ e 2 = τ P α P Δ β e 2 ,
and denoting by E H ^ eq H eq the error estimation matrix, the variance of the elements of E is β e ζ e 2 , and 0 μ e ζ e 2 β e 1 denotes the channel estimation reliability.
The computational complexity of channel estimation was studied and approximated for several channel estimation methods in [27]. Referring to it, we approximate the MMSE-based channel estimation complexity as O N B 3 τ p 3 . This complexity will be added to the precoding complexity already computed in Section 3.1.

3.2.4. Asymptotic SNR Analysis

Here, zero-forcing precoding is based on the estimated channel in (24), and the precoding matrix is given by F ^ 1 = N B N u Z ^ eq H ( Z ^ eq Z ^ eq H ) 1 , where Z ^ eq ( 1 / ζ e ) H ^ eq is the normalized estimated equivalent channel. The transmitted signal vector is thus given by
x ^ = P F ^ 1 q ,
and the relay receives the following signal:
y ^ RR = H 1 x ^ + n 1 .
Let us denote α ^ 1 = tr ( y ^ RR y ^ RR H ) , which, after a few manipulations, can be written as
α ^ 1 = P ( N B N u ) · tr H 1 Z ^ eq H ( Z ^ eq Z ^ eq H ) 1 q q H ( ( Z ^ eq Z ^ eq H ) 1 ) H Z ^ eq H 1 H + tr [ n 1 n 1 H ] .
To alleviate the derivations, we write Z ^ eq as Z ^ eq = 1 ζ e ( H eq + E ) . When N R = N B = N , we obtain the following convergences: H 1 Z ^ eq H N N β 1 ζ e H 2 H + 1 ζ e H 1 E H N , H 2 H 2 H N N β 2 I N u , Z ^ eq Z ^ eq H N N I N u , and n 1 n 1 H N N σ 1 2 I N , and since H 1 E H N N 0 N u , the normalization coefficient α ^ 1 asymptotically converges as
α ^ 1 N N P β 1 N u μ e + σ 1 2 .
Referring to (25), we can express μ e 2 as μ e 2 = τ P β 1 β 2 τ P β 1 β 2 + β 1 σ 3 2 + ( τ P β 2 N u + σ 3 2 ) σ 4 2 . This yields an expression of the signal received by the user under the form
y ^ u = P R P N R β 1 β 2 μ e 2 α ^ 1 ( N B N u ) q + P R α ^ 1 ( H 2 n 1 E x ) + n 2 .
Under these notations, we can derive the instantaneous SNR for user k as
Γ ^ 1 inst = P R P ( N B N u ) N R β 1 β 2 μ e 2 P R | E ( k , : ) x | 2 + P R | h ( 2 , k ) n 1 | 2 + α ^ 1 | n ( 2 , k ) | 2 ,
where E ( k , : ) is the k-th line of the channel estimation error matrix. Then, the average SNR can be expressed as
Γ ^ 1 avg = P R P ( N B N u ) N R β 1 β 2 μ e 2 P R P N u β e ( 1 μ e 2 ) + P R N R β 2 σ 1 2 + α ^ 1 σ 2 2 ,
and the end-to-end SNR in the large number of antennas regime at user k can be expressed as
Γ ^ 1 asym = P R P β 1 β 2 μ e 2 P R β 2 σ 1 2 N + P β 1 N u σ 2 2 μ e 2 N + σ 1 2 σ 2 2 N .
From (33), we can see clearly that the main influence of the channel estimation reliability μ e 2 is the decrease in the capacity of Scheme 1 due to the channel estimation error. Some typical power scaling options for Scheme 1 under imperfect CSI are summarized in Table 2. We observe that μ e 2 does not change the power scaling effects on Scheme 1 as N increases, which will also be verified in the numerical results in Section 5.
Finally, we can straightforwardly express the system capacity and asymptotic capacity, respectively, similarly to (13) and (14).

3.2.5. Outage Probability Analysis

In this subsection, we derive the outage probability of Scheme 1 using (31) and (29). The instantaneous SNR can be rewritten as
Γ ^ 1 inst = B λ 1 X + λ 2 Y + λ 3 Z ,
where B = P R P ( N B N u ) N R β 1 β 2 μ e 2 , λ 1 = P R P N u β e ( 1 μ e 2 ) / 2 , λ 2 = P R N R β 2 σ 1 2 / 2 , λ 3 = α ^ 1 σ 2 2 / 2 , and X, Y, and Z are three independent chi-squared distributed variates with two degrees of freedom. Note that λ 1 X , λ 2 Y , and λ 3 Z are gamma variates with shape parameter 1 and respective scale parameters b 1 = 1 / ( 2 λ 1 ) , b 2 = 1 / ( 2 λ 2 ) , and b 3 = 1 / ( 2 λ 3 ) . Then, by following the same approach as in [25], we obtain the CDF of the instantaneous SNR in Scheme 1 under imperfect CSI as
P ^ 1 out ( Γ th ) = b 3 b 2 ( b 3 b 2 ) exp b 1 B Γ th b 3 b 1 ( b 3 b 1 ) exp b 2 B Γ th + b 2 b 1 ( b 2 b 1 ) exp b 3 B Γ th ( b 3 b 2 ) ( b 3 b 1 ) ( b 2 b 1 ) .
This yields the outage probability expression that we will evaluate numerically later in Section 5.

4. Hop-by-Hop Zero-Forcing Precoding at BS and Relay

4.1. Case of Perfect CSI

In this second scheme, hereafter called Scheme 2, the base station precodes the signal using the channel response of the first hop; then, the relay receives, without interference, the users’ signals with N u antennas, precodes them based on the channel response of the second hop, and normalizes, amplifies, and retransmits the signals to the users. In this case, N R 1 = N u and N R 2 = N R . The precoding matrix at the base station is given by F 1 = N B N u Z 1 + , and the relay processing matrix is given by F 2 = N R N u Z 2 + . Therefore, the precoding complexity of Scheme 2 can be approximated as O ( 2 N B N u 2 + N u 3 ) + O ( 2 N R 2 N u 2 + N u 3 ) .
Assuming the same large number of antennas at both the base station and the relay, i.e., N R 2 = N B = N , the complexity can be rewritten as 2 O ( N N u 2 + N u 3 ) . Under these considerations, we derive the normalization coefficient as
α 2 = P β 1 ( N N u ) 2 tr Z 2 + q ( Z 2 + q ) H + ( N N u ) Z 2 + n 1 ( Z 2 + n 1 ) H .

4.1.1. Asymptotic SNR Analysis

Based on [24] (Lemma 1), when N R = N B = N , we obtain
Z 2 Z 2 H N N I N u .
Then, for the normalization coefficient
α 2 N N P β 1 N u .
The received signal at user k is given by
y ( u , k ) = P P R β 1 β 2 ( N B N u ) ( N R N u ) α 2 q k + P R β 2 ( N R N u ) α 2 n ( 1 , k ) + n ( 2 , k ) ,
where n ( 1 , k ) is the k-th element of vector n 1 .
Now, we can derive the instantaneous and average SNRs, respectively, as
Γ 2 inst = P R P ( N R N u ) ( N B N u ) β 1 β 2 P R | n ( 1 , k ) | 2 + α 2 | n ( 2 , k ) | 2 ,
Γ 2 avg = P P R ( N R N u ) ( N B N u ) β 1 β 2 P R ( N R N u ) β 2 σ 1 2 + α 2 σ 2 2 .
When N B = N R = N and taking into consideration (38), the asymptotic SNR for this scenario can be expressed as
Γ 2 asym N P R P N β 1 β 2 P R β 2 σ 1 2 + P β 1 N u σ 2 2 .
As we have observed in (12) for Scheme 1, Equation (42) also shows that the asymptotic SNR in Scheme 2 grows to infinity as N . The asymptotic SNRs for a few typical power scaling options are summarized in Table 3.
In Table 3, we analyze a few typical power scaling options for Scheme 2 under perfect CSI to take advantage of using mMIMO and to improve the energy efficiency while maintaining the desired capacity. Therefore, we propose some practical cases of asymptotic SNR according to specific values of a and b.
From Table 1 and Table 3, we see clearly that Scheme 2 outperforms Scheme 1 for all power scaling options. Moreover, for a large number of antennas, increasing a and b leads to a dramatic reduction in the performance of Scheme 1, i.e., a + b 1 .
Finally, we derive the system capacity and asymptotic capacity, respectively, as
C 2 = log 2 ( 1 + N u Γ 2 avg ) ,
C 2 asym = log 2 ( 1 + N u Γ 2 asym ) .

4.1.2. Outage Probability Analysis

We adopt the same approach as for Scheme 1 and derive the CDF of the approximate SNR in Scheme 2 as
P 2 out ( Γ th ) = c 1 exp c 2 C Γ th c 2 exp c 1 C Γ th c 1 c 2 .
where C = P R P ( N R N u ) ( N B N u ) β 1 β 2 , c 1 = 1 / ( P R ( N R N u ) β 2 σ 1 2 ) , and c 2 = 1 / α 2 σ 2 2 .

4.2. Case of Imperfect CSI

Here, for Scheme 2, we adopt a cascaded channel estimation, wherein the relay estimates the second-hop channel response and transmits pilot symbols to the BS, which estimates the first-hop channel response.

4.2.1. MMSE Channel Estimation

First, the second-hop channel estimation is given by
H ^ 2 T = τ P β 2 τ P β 2 + σ 3 2 Y P , 1 ,
where Y P , 1 is the pilot symbol matrix received by the relay after despreading. Similarly to the perfect CSI case, let us denote by ζ 2 2 the variance of the elements of H ^ 2 T , which is given by
ζ 2 2 = τ P β 2 2 τ P β 2 + σ 3 2 ,
and let E 2 H ^ 2 H 2 denote the channel estimation error matrix over the second hop, for which the elements’ variance is given by β 2 ζ 2 2 , and 0 μ 2 = ζ 2 2 β 2 1 denotes the channel estimation reliability in this case.
Second, and in order to estimate the first-hop channel, the relay sends its own pilot symbols to the BS, which obtains an estimation given by
H ^ 1 T = τ P β 1 τ P β 1 + σ 4 2 Y P , 2 ,
where Y P , 2 is the pilot matrix received by the BS after despreading. Again, ζ 1 2 denotes the variance of the elements of H ^ 1 T and is given by
ζ 1 2 = τ P β 1 2 τ P β 1 + σ 4 2 ,
E 1 H ^ 1 H 1 is the channel estimation error matrix over the first hop, for which the elements’ variance is given by β 1 ζ 1 2 , and 0 μ 1 = ζ 1 2 β 1 1 denotes the channel estimation reliability.
In this case, the channel estimation complexity is approximated as O ( N B 3 τ p 3 + N R 3 τ p 3 ) .

4.2.2. Asymptotic SNR Analysis

In this case, zero-forcing precoding is based on the estimated channels in (46) and (48), with precoding matrices F ^ 1 N B N u Z ^ 1 + and F ^ 2 N R N u Z ^ 2 + , where Z ^ 1 = ( 1 / ζ 1 2 ) H ^ 1 and Z ^ 2 = ( 1 / ζ 2 2 ) H ^ 2 are the normalized estimated channel matrices of the first and the second hops, respectively.
Under these notations, the signal received at the relay side can be expressed as
y ^ RR = ( N B N u ) P ( H ^ 1 E 1 ) Z ^ 1 + q + n 1 .
In this case, α ^ 2 can be expressed as
α ^ 2 = P ( N B N u ) · tr H 1 Z ^ 1 H ( Z ^ 1 Z ^ 1 H ) 1 q × q H ( ( Z ^ 1 Z ^ 1 H ) 1 ) H Z ^ 1 H 1 H + tr [ n 1 n 1 H ] .
To alleviate the derivations, we write Z ^ 1 as Z ^ 1 = 1 ζ 1 2 ( H 1 + E 1 ) . When N R = N B = N , we obtain H 1 Z ^ 1 H N N β 1 ζ 1 2 I N u + H 1 E 1 H N , Z ^ 1 Z ^ 1 H N N I N u , and n 1 n 1 H N N σ 1 2 I N u , and since H 1 E 1 H N N 0 N u , the normalization coefficient α ^ 2 asymptotically converges to
α ^ 2 N N P β 1 N u μ 1 2 .
Thus, the signal received by all the users is given by
y u = P R ( N R N u ) α ^ 2 ( Z ^ 2 E 2 ) Z ^ 2 + P ζ 1 2 ( N B N u ) q P ( N B N u ) E 1 Z ^ 1 + q + n 1 + n 2 .
Similarly to the analysis in Section 3.2.4 and after a few simple manipulations, we can derive the instantaneous and the average SNRs, which are given, respectively, by
Γ ^ 2 inst = P R P ( N B N u ) ( N R N u ) β 1 β 2 μ 1 2 μ 2 2 P ( N B N u ) β 1 μ 1 2 | E 2 ( k , : ) x 2 | 2 + P R ( N R N u ) β 2 μ 2 2 | E 1 ( k , : ) x 1 | 2 + P R ( N R N u ) β 2 μ 2 2 | ( n ( k , 1 ) | 2 + α ^ 2 | n ( k , 2 ) | 2 ,
and
Γ ^ 2 avg = P R P ( N B N u ) ( N R N u ) β 1 β 2 μ 1 2 μ 2 2 P R P N u ( N B N u ) β 1 β 2 μ 1 2 ( 1 μ 2 2 ) + P R P N u ( N R N u ) β 1 β 2 μ 2 2 ( 1 μ 1 2 ) + P R ( N R N u ) β 2 μ 2 2 σ 1 2 + α ^ 2 σ 2 2 ,
where E 1 ( k , : ) and E 2 ( k , : ) are the k-th lines of channel estimation error matrix E 1 and E 2 , respectively, and x 1 = P ( N B N u ) Z ^ 1 + q and x 2 = P R ( N R N u ) Z ^ 2 + q . We also assume that H 2 Z ^ 2 + = ζ 2 2 I N u . The end-to-end SNR in the large number of antennas regime at user k can be expressed as
Γ ^ 2 asym = P R P N β 1 β 2 μ 1 4 μ 2 2 P R β 2 σ 1 2 μ 1 2 μ 2 2 + P β 1 N u σ 2 2 .
Intuitively, Scheme 1 is more sensitive in terms of channel estimation reliability compared to Scheme 2. This is mainly due to the pilot noise amplification at the relay. A few typical power scaling options for Scheme 2 under imperfect CSI are summarized in Table 4.
Finally, we can express the system capacity and asymptotic capacity as in (43) and (44), respectively.

4.2.3. Outage Probability Analysis

Let us now derive the outage probability of Scheme 2 under imperfect channel estimation using (54) and (52). The instantaneous SNR can be rewritten as
Γ ^ 1 inst = D λ 1 X + λ 2 Y + λ 3 W + λ 4 Z ,
where D = P R P ( N B N u ) ( N R N u ) β 1 β 2 μ 1 2 μ 2 2 , λ 1 = P P R ( N B N u ) β 1 β 2 μ 1 2 ( 1 μ 2 2 ) / 2 , λ 2 = P P R N u ( N R N u ) β 1 β 2 μ 2 2 ( 1 μ 1 2 ) / 2 , λ 3 = P R ( N R N u ) β 2 μ 2 2 σ 1 2 / 2 , λ 4 = α ^ 2 σ 2 2 / 2 , and X, Y, W, and Z are independent chi-squared distributed variates with two degrees of freedom. Then, we obtain the CDF of the approximate SNR in Scheme 2 under imperfect CSI as
P ^ 2 out ( Γ th ) = j = 1 4 i = 1 , i j 4 d i i = 1 , i j 4 ( d i d j ) exp d j D Γ th ,
where d 1 = 1 / ( 2 λ 1 ) , d 2 = 1 / ( 2 λ 2 ) , d 3 = 1 / ( 2 λ 3 ) , and d 4 = 1 / ( 2 λ 4 ) . This yields the outage probability expression that we will evaluate numerically later in Section 5.

5. Numerical Results

Using Monte Carlo simulations, we now validate our analytical analysis and give a comparative discussion of both schemes taking into consideration the precoding complexity. Without loss of generality, we assume that σ 1 = σ 2 = σ 3 = σ 4 = 1 , β 1 = β 2 = 1 , and τ p = 10 .
Please note that in all figures that will be presented in this section, the markers represent the results of simulations, while solid lines represent the analytical and asymptotic results under the perfect CSI assumption, and the dashed lines represent the results with imperfect MMSE channel estimation.

5.1. Capacity and Energy Efficiency

In this subsection, we first examine the capacity and the energy efficiency (EE) of both schemes, where EE is defined as EE = C / ( P + P R ) . (For simplicity, we only consider transmission energy here, as it is the major component, and we neglect other energy consumption in the system. An exhaustive energy efficiency framework will complement this work in the future.) Figure 2 and Figure 3 show the simulated capacity together with the average and asymptotic capacity for different values of a and b, while Figure 4 and Figure 5 show the energy efficiency for different values of a and b. All results are given for both perfect and imperfect CSI.
Clearly, the asymptotic values presented in Table 1, Table 2, Table 3 and Table 4 can perfectly predict the performance of Scheme 1 and Scheme 2 under perfect and imperfect CSI. In fact, as N grows infinitely large, the capacity when a = b = 0 and a = b = 0.3 has no upper bound. In the other cases, as N increases, the capacity increases towards a constant asymptotic value—except for a = b = 1 with Scheme 1, where it decreases towards zero. Moreover, as discussed in Section 4.1.1, Scheme 2 outperforms Scheme 1 for all power scaling options when N is large enough. Further, observing Figure 4 and Figure 5, we can conclude that Scheme 2 is the best in terms of energy efficiency when combined with relevant power scaling.
However, in terms of complexity, Scheme 2 is worse than Scheme 1. Hence, for a fair comparison, and in order to get the same complexity, we select different values of N for each scheme in what follows. Let us denote by N 1 and N 2 the numbers of antennas for Scheme 1 and Scheme 2, respectively. Based on the complexity approximations for channel estimation and precoding for each scheme, and in order to get the same complexity for both schemes, the following constraint is respected:
N 1 = 2 N 2 + N u 2 , for perfect CSI ,
and
N 1 2 3 N 2 , for imperfect CSI .
Figure 6, Figure 7, Figure 8 and Figure 9 show the capacity and energy efficiency for both schemes with the same complexity level. We clearly see that with the additional complexity constraint, Scheme 1 now always outperforms Scheme 2 for the perfect CSI condition—except when the power scaling at the relay is very high, i.e., b = 1 . However, under imperfect CSI, Scheme 2 remains a better choice at the same complexity level for all proposed a and b values.

5.2. Outage Probability

In this subsection, we compare the studied schemes in terms of outage probability, and for space considerations, we limit our analysis to the case of power scaling at the relay, i.e., a = 0 and b = 1 . Moreover, in all figures, we assume that both schemes have the same complexity, i.e., the constraints in (59) and (60) are respected. In Figure 10, we first note that the simulation results confirm the accuracy of the analytical expressions obtained in Section 3.1.2, Section 3.2.5, Section 4.1.2, and Section 4.2.3. We conclude that Scheme 1 outperforms Scheme 2 in term of outage probability when assuming the same complexities and a perfect CSI for some numbers of transmitting antennas. Under imperfect CSI, we see that Scheme 2 is always the best.

6. Conclusions

In this paper, we have analyzed the end-to-end capacity, its asymptotic approximation, and the outage probability of an outdoor-to-indoor dual-hop full-duplex mmWave multi-user system for which the transmitting nodes are equipped with massive numbers of antennas while the final users are equipped with single antennas. Two precoding schemes were proposed, and the performance metrics were derived for each scheme under both perfect and imperfect CSI conditions. An approximation of the precoding complexity was also given for each precoding scheme. To respect the power budget allocated to the relay, the normalization coefficient was approximated, capitalizing on the law of large numbers. Monte Carlo simulation results were presented to confirm the accuracy of our analysis. Practical recommendations were formulated at the end of the analysis for both pure performance and for fair performance–complexity trade-off comparisons.

Author Contributions

Conceptualization, N.B., H.C. and M.B.; Data curation, N.B., H.C. and M.B.; Formal analysis, N.B., H.C. and M.B.; Investigation, N.B., H.C. and M.B.; Methodology, N.B., H.C. and M.B.; Project administration, N.B., H.C. and M.B.; Resources, N.B., H.C. and M.B.; Software, N.B., H.C. and M.B.; Supervision, N.B., H.C. and M.B.; Validation, N.B., H.C. and M.B.; Visualization, N.B., H.C. and M.B.; Writing—original draft, N.B., H.C. and M.B.; Writing—review and editing, N.B., H.C. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All obtained data in this work is the result of the mathematical calculations detailed in the paper. Supporting simulation scripts can be shared upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bounouader, N.; Chafnaji, H.; Benjillali, M. Power Normalization and Precoding for Massive MIMO in Outdoor-to-Indoor mmWave Relaying. In Proceedings of the 6th International Conference on Advanced Communication Technologies and Networking (CommNet), Rabat, Morocco, 11–13 December 2023; pp. 1–7. [Google Scholar]
  2. Hong, W.; Jiang, Z.H.; Yu, C.; Hou, D.; Wang, H.; Guo, C.; Hu, Y.; Kuai, L.; Yu, Y.; Jiang, Z.; et al. The role of millimeter-wave technologies in 5G/6G wireless communications. IEEE J. Microw. 2021, 1, 101–122. [Google Scholar] [CrossRef]
  3. Sakaguchi, K.; Yoneda, T.; Iwabuchi, M.; Murakami, T. Mmwave massive Analog relay MIMO. ITU J. Future Evol. Technol. 2021, 2, 43–55. [Google Scholar] [CrossRef]
  4. Nguyen, N.T.; Nguyen, H.N.; Nguyen, N.L.; Le, A.T.; Nguyen, T.N.; Voznak, M. Performance analysis of NOMA-based hybrid satellite-terrestrial relay system using mmWave technology. IEEE Access 2023, 11, 10696–10707. [Google Scholar] [CrossRef]
  5. Li, D.; Xu, G.; Gao, M.; Song, Z.; Zhang, Q.; Zhang, W. Performance Analyses of RIS-Assisted Stochastic UAV mmWave Relay Communication System With Moment Matching Estimation. IEEE Wirel. Commun. Lett. 2024, 13, 1198–1202. [Google Scholar] [CrossRef]
  6. Marzetta, T.L. How much training is required for multiuser MIMO? In Proceedings of the 40-th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 29 October–1 November 2006; pp. 359–363. [Google Scholar]
  7. Marzetta, T.L. Noncooperative Cellular Wireless with Unlimited Numbers of Base Station Antennas. IEEE Trans. Wirel. Commun. 2010, 9, 3590–3600. [Google Scholar] [CrossRef]
  8. Grossman, D. mmWave for Fixed Wireless: What We’ve Learned; Technical Report, Strategy Analytics; Strategy Analytics: Ottawa, ON, Canada, 2022. [Google Scholar]
  9. Ntontin, K.; Verikoukis, C. Relay-aided outdoor-to-indoor communication in millimeter-wave cellular networks. IEEE Syst. J. 2019, 14, 2473–2484. [Google Scholar] [CrossRef]
  10. Albreem, M.A.; Juntti, M.; Shahabuddin, S. Massive MIMO detection techniques: A survey. IEEE Commun. Surv. Tutor. 2019, 21, 3109–3132. [Google Scholar] [CrossRef]
  11. Zhang, Z.; Long, K.; Vasilakos, A.V.; Hanzo, L. Full-Duplex Wireless Communications: Challenges, Solutions, and Future Research Directions. Proc. IEEE 2016, 104, 1369–1409. [Google Scholar] [CrossRef]
  12. Koc, A.; Le-Ngoc, T. Full-duplex mmWave massive MIMO systems: A joint hybrid precoding/combining and self-interference cancellation design. IEEE Open J. Commun. Soc. 2021, 2, 754–774. [Google Scholar] [CrossRef]
  13. Osama, I.; Rihan, M.; Elhefnawy, M.; Eldolil, S.; Malhat, H.A.E.A. A review on Precoding Techniques For mm-Wave Massive MIMO Wireless Systems. Int. J. Commun. Netw. Inf. Secur. 2022, 14, 26–36. [Google Scholar] [CrossRef]
  14. Marzetta, T.L.; Ngo, H.Q. Fundamentals of Massive MIMO; Cambridge University Press: Cambridge, UK, 2016. [Google Scholar]
  15. Chae, C.B.; Tang, T.; Heath, R.W.; Cho, S. MIMO relaying with linear processing for multiuser transmission in fixed relay networks. IEEE Trans. Signal Process. 2008, 56, 727–738. [Google Scholar] [CrossRef]
  16. Lee, C.; Chae, C.B.; Kim, T.; Choi, S.; Lee, J. Network massive MIMO for cell-boundary users: From a precoding normalization perspective. In Proceedings of the IEEE Globecom Workshops, Anaheim, CA, USA, 3–7 December 2012; pp. 233–237. [Google Scholar]
  17. Lin, H.; Zhu, H.; Yang, Y.; Huang, G.; Huang, K.; Chen, Y. Robust Channel estimation for 5G mmWave massive MIMO relay communication systems. In Proceedings of the IEEE 21st International Conference on Communication Technology (ICCT), Tianjin, China, 13–16 October 2021; pp. 289–294. [Google Scholar]
  18. Eltayeb, M.E. Relay-aided channel estimation for mmWave systems with imperfect antenna arrays. In Proceedings of the IEEE International Conference on Communications (ICC), Shanghai, China, 20–24 May 2019; pp. 1–5. [Google Scholar]
  19. Tabataba, F.S.; Sadeghi, P.; Pakravan, M.R. Outage Probability and Power Allocation of Amplify and Forward Relaying with Channel Estimation Errors. IEEE Trans. Wirel. Commun. 2011, 10, 124–134. [Google Scholar] [CrossRef]
  20. Cao, Y.; Liang, X.; Wang, X.; Li, Y.; Fu, Q.; Wang, D.; You, X. Performance of Cell-Free Massive MIMO Under Imperfect Channel State Information and Reciprocity Calibration. IEEE Syst. J. 2023, 17, 4383–4394. [Google Scholar] [CrossRef]
  21. Niu, Y.; Li, Y.; Jin, D.; Su, L.; Vasilakos, A.V. A survey of millimeter wave communications (mmWave) for 5G: Opportunities and challenges. Wirel. Netw. 2015, 21, 2657–2676. [Google Scholar] [CrossRef]
  22. Akdeniz, M.R.; Liu, Y.; Samimi, M.K.; Sun, S.; Rangan, S.; Rappaport, T.S.; Erkip, E. Millimeter Wave Channel Modeling and Cellular Capacity Evaluation. IEEE J. Sel. Areas Commun. 2014, 32, 1164–1179. [Google Scholar] [CrossRef]
  23. Pandey, B.C.; Mohammed, S.K.; Raviteja, P.; Hong, Y.; Viterbo, E. Low complexity precoding and detection in multi-user massive MIMO OTFS downlink. IEEE Trans. Veh. Technol. 2021, 70, 4389–4405. [Google Scholar] [CrossRef]
  24. Cui, H.; Song, L.; Jiao, B. Multi-pair two-way amplify-and-forward relaying with very large number of relay antennas. IEEE Trans. Wirel. Commun. 2014, 13, 2636–2645. [Google Scholar]
  25. Coelho, C.A. The generalized integer Gamma distribution: A basis for distributions in multivariate analysis. J. Multivar. Anal. 1998, 64, 86–102. [Google Scholar] [CrossRef]
  26. Ngo, H.Q.; Larsson, E.G.; Marzetta, T.L. Energy and spectral efficiency of very large multiuser MIMO systems. IEEE Trans. Commun. 2013, 61, 1436–1449. [Google Scholar]
  27. Hawej, M.; Shayan, Y.R. Low-complexity channel estimation for time division duplex massive multi-user multi-input multi-output systems. IET Commun. 2021, 15, 43–50. [Google Scholar] [CrossRef]
Figure 1. Adopted system model with massive-MIMO-based outdoor-to-indoor relaying with multiple indoor users.
Figure 1. Adopted system model with massive-MIMO-based outdoor-to-indoor relaying with multiple indoor users.
Electronics 13 01857 g001
Figure 2. Capacity and asymptotic capacity for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 2. Capacity and asymptotic capacity for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g002
Figure 3. Capacity and asymptotic capacity for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 3. Capacity and asymptotic capacity for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g003
Figure 4. Energy efficiency for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 4. Energy efficiency for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g004
Figure 5. Energy efficiency for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 5. Energy efficiency for different values of a and b, with N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g005
Figure 6. Capacity of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 6. Capacity of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g006
Figure 7. Capacity of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 6 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 7. Capacity of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 6 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g007
Figure 8. Energy efficiency of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 8. Energy efficiency of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g008
Figure 9. Energy efficiency of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Figure 9. Energy efficiency of Scheme 1 and Scheme 2 at equal complexities for different values of a and b when N u = 5 , E = 2 , and E R = 10 under perfect and imperfect CSI.
Electronics 13 01857 g009
Figure 10. Analytical and simulated outage probability assuming the same complexities for a = 0 and b = 1 , E = 10 , E R = 10 , and Γ th = 0.3 .
Figure 10. Analytical and simulated outage probability assuming the same complexities for a = 0 and b = 1 , E = 10 , E R = 10 , and Γ th = 0.3 .
Electronics 13 01857 g010
Table 1. Asymptotic SNRs with different power scaling—Scheme 1 and perfect CSI.
Table 1. Asymptotic SNRs with different power scaling—Scheme 1 and perfect CSI.
Power Scaling ParametersAsymptotic SNR ( N )
a + b < 1
a + b = 1 , a 0 and b 0 E E R β 1 β 2 / σ 1 2 σ 2 2
a = 1 and b = 0 E E R β 1 β 2 / [ ( β 2 E R + σ 2 2 ) σ 1 2 ]
a = 0 and b = 1 E E R β 1 β 2 / [ ( β 1 N u E + σ 1 2 ) σ 2 2 ]
a + b > 1 0
Table 2. Asymptotic SNRs with different power scaling—Scheme 1 and imperfect CSI.
Table 2. Asymptotic SNRs with different power scaling—Scheme 1 and imperfect CSI.
Power Scaling ParametersAsymptotic SNR ( N )
a + b < 1
a + b = 1 , a 0 and b 0 E E R β 1 β 2 μ e 2 / σ 1 2 σ 2 2
a = 1 and b = 0 E E R β 1 β 2 μ e 2 / [ ( β 2 E R + σ 2 2 ) σ 1 2 ]
a = 0 and b = 1 E E R β 1 β 2 μ e 2 / [ ( β 1 N u E μ e 2 + σ 1 2 ) σ 2 2 ]
a + b > 1 0
Table 3. Asymptotic SNRs with different power scaling—Scheme 2 and perfect CSI.
Table 3. Asymptotic SNRs with different power scaling—Scheme 2 and perfect CSI.
Power Scaling ParametersAsymptotic SNR ( N )
0 a , b < 1
a = 1 and 0 b < 1 E β 1 / σ 1 2
0 a < 1 and b = 1 E R β 2 / N u σ 2 2
a = b = 1 E R E β 1 β 2 / ( E R β 2 σ 1 2 + E β 1 N u σ 2 2 )
Table 4. Asymptotic SNRs with different power scaling—Scheme 2 and imperfect CSI.
Table 4. Asymptotic SNRs with different power scaling—Scheme 2 and imperfect CSI.
Power Scaling ParametersAsymptotic SNR ( N )
0 a , b < 1
a = 1 and 0 b < 1 E β 1 μ 1 2 / σ 1 2
0 a < 1 and b = 1 E R β 2 μ 1 4 μ 2 2 / N u σ 2 2
a = b = 1 E R E β 1 β 2 μ 1 4 μ 2 2 ( E R β 2 σ 1 2 μ 1 2 μ 2 2 + E β 1 N u σ 2 2 )
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bounouader, N.; Chafnaji, H.; Benjillali, M. Outdoor-to-Indoor mmWave Relaying with Massive MIMO: Impact of Imperfect Channel Estimation. Electronics 2024, 13, 1857. https://doi.org/10.3390/electronics13101857

AMA Style

Bounouader N, Chafnaji H, Benjillali M. Outdoor-to-Indoor mmWave Relaying with Massive MIMO: Impact of Imperfect Channel Estimation. Electronics. 2024; 13(10):1857. https://doi.org/10.3390/electronics13101857

Chicago/Turabian Style

Bounouader, Nawal, Houda Chafnaji, and Mustapha Benjillali. 2024. "Outdoor-to-Indoor mmWave Relaying with Massive MIMO: Impact of Imperfect Channel Estimation" Electronics 13, no. 10: 1857. https://doi.org/10.3390/electronics13101857

APA Style

Bounouader, N., Chafnaji, H., & Benjillali, M. (2024). Outdoor-to-Indoor mmWave Relaying with Massive MIMO: Impact of Imperfect Channel Estimation. Electronics, 13(10), 1857. https://doi.org/10.3390/electronics13101857

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop