arXiv is now an independent nonprofit! Learn more
License: arXiv.org perpetual non-exclusive license
arXiv:1809.03117v1 [cs.IT] 10 Sep 2018

Mixed-ADC/DAC Multipair Massive MIMO Relaying Systems: Performance Analysis and Power OptimizationThanks: This work was supported in part by the National Natural Science Foundation of China (Nos. 61601020 and 61725101), the Beijing Natural Science Foundation (Nos. 4182049 and L171005), the open research fund of National Mobile Communications Research Laboratory, Southeast University (No. 2018D04), Key Laboratory of Optical Communication and Networks (No. KLOCN2018002), National Key Research and Development Program (No. 2016YFE0200900), Major projects of Beijing Municipal Science and Technology Commission (No. Z181100003218010), National key research and development program under (No. 2016YFE0200900), and the Natural Sciences and Engineering Research Council of Canada (NSERC) through its Discovery program.Thanks: J. Zhang is with the School of Electronic and Information Engineering, Beijing Jiaotong University, Beijing 100044, P. R. China. He is also with National Mobile Communications Research Laboratory, Southeast University, Nanjing 210096, P. R. China (e-mail: jiayizhang@bjtu.edu.cn).Thanks: L. Dai is with Department of Electronic Engineering, Tsinghua University, Beijing 100084, P. R. China (e-mail: daill@tsinghua.edu.cn).Thanks: Z. He was with Department of Electronic Engineering, Tsinghua University, Beijing 100084, P. R. China. He is now with School of Electrical and Computer Engineering, Georgia Institute of Technology, Atlanta, GA 30332 USA (e-mail: heziyan@gatech.edu).Thanks: B. Ai is with the State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University, Beijing 100044, China (e-mail: boai@bjtu.edu.cn). Thanks: O. A. Dobre is with the Faculty of Engineering and Applied Science, Memorial University, St. John’s, NL A1B 3X5, Canada (e-mail: odobre@mun.ca).

Jiayi Zhang    Linglong Dai    Ziyan He Affiliation: Bo Ai,  and Octavia A. Dobre, 
Abstract

High power consumption and expensive hardware are two bottlenecks for practical massive multiple-input multiple-output (mMIMO) systems. One promising solution is to employ low-resolution analog-to-digital converters (ADCs) and digital-to-analog converters (DACs). In this paper, we consider a general multipair mMIMO relaying system with a mixed-ADC/DAC architecture, in which some antennas are connected to low-resolution ADCs/DACs, while the rest of the antennas are connected to high-resolution ADCs/DACs. Leveraging on the additive quantization noise model, both exact and approximate closed-form expressions for the achievable rate are derived. It is shown that the achievable rate can approach the unquantized one by using only 2-3 bits of resolutions. Moreover, a power scaling law is presented to reveal that the transmit power can be scaled down inversely proportional to the number of antennas at the relay. We further propose an efficient power allocation scheme by solving a complementary geometric programming problem. In addition, a trade-off between the achievable rate and power consumption for different numbers of low-resolution ADCs/DACs is investigated by deriving the energy efficiency. Our results reveal that the large antenna array can be exploited to enable the mixed-ADC/DAC architecture, which significantly reduces the power consumption and hardware cost for practical mMIMO systems.

Index Terms: 
Massive MIMO, multipair relay, mixed-ADC/DAC, achievable rate, energy efficiency.

I Introduction

As one of the disruptive technologies for the fifth-generation (5G) wireless communications, massive multiple-input multiple-output (mMIMO) has attracted extensive research interests in recent years [1, 2, 3, 4]. By exploiting quasi-orthogonal random channel vectors between different users, mMIMO can mitigate the inter-user interference to provide high spectral efficiency and energy efficiency via simple linear signal processing, e.g., maximum-ratio combining (MRC) and zero-forcing (ZF) precoding. On the other hand, relaying is an important way of extending coverage and providing uniform service. The inter-user interference of the multiuser relaying system can be suppressed by equipping the relay with a large number of antennas [5, 6, 7].

The practical implementation of an mMIMO relaying system with hundreds or even thousands of antennas is a significant challenge [8, 9, 10]. For example, the perfect synchronization is difficult in mMIMO relaying system. One possible solution is to employ rateless network coding [11]. Typically, each antenna in mMIMO systems is connected to an analog-to-digital converter (ADC) and a digital-to-analog converters (DAC) in the radio frequency (RF) chain, respectively. It is well known that the power consumption and hardware cost of ADCs and DACs linearly increase with the bandwidth and exponentially increase with the number of quantization bits [12]. Thus, high-resolution ADCs and DACs (e.g., 8-12 bits for commercial use) will result in high power consumption and hardware cost in practical mMIMO relaying systems.

To solve this challenging problem, a promising solution is to replace power-hungry high-resolution ADCs and DACs (e.g., 8-12 bits) with low power low-resolution ADCs and DACs (e.g., 1-3 bits) [13, 14]. However, significant signal processing challenges and complex front-end designs (e.g., channel estimation, phase/frequency synchronization, and multiuser detection) inevitably occur due to the strong nonlinear characteristic of coarse quantization [15]. As recently reported in [16], a high signal-to-noise ratio (SNR) channel estimation error floor exists due to the one-bit quantization. Furthermore, the authors in [17] proposed a mixed-ADC architecture, where only a small fraction of ADCs are high-resolution, to facilitate the aforementioned signal processing and the establishment of front-end designs. For example, the CSI for each antenna can be obtained by using the high-resolution ADCs in a round-robin manner, which has been clearly explained in [17]. Moreover, the mixed-ADC architecture is economically beneficial and easier to implement compared with architectures with uniform converter resolution, as it adds some antennas with low-resolution ADCs to the existing high-resolution conventional MIMO system [18].

I-A Related Works

Most of recent works focused on the single-hop mMIMO system with a mixed-ADC architecture. For instance, the mutual information of mixed-ADC mMIMO systems has been investigated in [17, 19], which reveals that the mixed-ADC architecture is able to approach the ideal channel capacity of unquantized systems over both frequency-flat [17] and frequency-selective fading channels [19]. In addition, the achievable rate performance of multi-user mMIMO systems with a mixed-ADC architecture is comparable for Rayleigh [20] and Rician fading channels [21]. This architecture can achieve a better energy-rate trade-off compared with the ideal infinite-resolution and low-resolution ADC architectures. By applying probabilistic Bayesian inference, a family of detectors for mixed-ADC mMIMO systems was developed in [22]. It proves that the available high-resolution ADCs are practically essential since they can effectively eliminate the error floor of a relaxed Bayesian detector. Given the energy constraint at the base station (BS), the sum achievable rate has been maximized in [23], which shows that the optimal rate can be obtained by using only one-bit ADCs in most realistic scenarios. Moreover, [24] considered the downlink mMIMO with both mixed-resolution DACs at the BS and mixed-resolution ADCs at the user side. These important contributions have shown that the power consumption and hardware cost of the single-hop mMIMO system with a mixed-ADC architecture can be considerably reduced while keeping most of the gains in the achievable rate.

In contrast to single-hop systems, very little attention has been paid to the two-hop mMIMO relaying system with both mixed-resolution ADCs and mixed-resolution DACs. Very recently, the authors in [25, 26, 27, 28, 29] investigated the performance of a multipair mMIMO relaying system with low-resolution ADCs and DACs at the relay. The achievable rate of such a system is limited by using very coarse quantization (e.g., one-bit). In this paper, we consider a more general architecture, where ADCs and DACs with arbitrary resolution profile are employed at the relay to achieve a possibly higher rate.

I-B Contributions

In this paper, we focus on a general two-hop mixed-ADC/DAC mMIMO relaying system, where some antennas are connected to low-resolution ADCs/DACs, while the rest of the antennas are connected to high-resolution ADCs/DACs. This study aims to analyse the performance analysis of the multipair mMIMO relaying system with arbitrary quantization noise, which is in contrast to the previous study [25] that employs only one-bit ADCs and DACs. We demonstrate that the achievable rate of the considered system can approach that of the ideal unquantized system. The main contributions of this paper are summarized as follows:

  • Leveraging on AQNM, we present a unified framework to derive the exact closed-form expressions for the achievable rate of mixed-ADC/DAC mMIMO systems. Compared with the Bussgang theorem used in [25], the AQNM can offer analytical tractability for multi-bit ADCs and DACs. Furthermore, approximate achievable rate expressions are derived by using asymptotic arguments. These results can provide insights into the effects of the number of relay antennas, user transmit power, quantization bits, and the fraction of high-resolution quantizers on the achievable rate, respectively.

  • The power-scaling law of the mixed-ADC/DAC architecture is investigated for power saving in the data transmission phase of the relay. Our results reveal that when the number of relay antennas, MM, gets asymptotically large, the transmit power of each antenna can be scaled down by 1/M1/M without any rate loss for the considered system. Moreover, the achievable rate gap between the mixed-ADC/DAC relay system and the unquantized system is a constant in the low power regime.

  • In order to compensate for the rate degradation caused by the coarse quantization, a low-complexity power allocation algorithm is proposed for the considered system. The power allocation problem can be solved by transforming it into a sequence of geometric programming (GP) problems.

  • Finally, using a generic power consumption model, we study the effects of the fraction of high-resolution ADCs and quantization bits on energy efficiency. In order to maximize energy efficiency, the optimal number of quantization bits is derived through numerical computations. Furthermore, our analysis proves that the considered system can significantly reduce power consumption and hardware cost while maintaining considerable performance.

I-C Outline

The remaining parts of the paper are structured as follows. The mixed-ADC/DAC multipair mMIMO relaying system model is briefly introduced in Section II. Both exact and asymptotic achievable rate expressions are derived in Section III. Moreover, the power-scaling laws for the considered system is presented in Section III. In Section V, a simple power allocation scheme is proposed to compensate for the rate loss. Furthermore, numerical results are provided in Section V to illustrate the effect of various system parameters on the achievable rate and energy efficiency. Finally, key findings are concluded in Section VI. Most of the mathematical proofs are given in Appendices A and B.

I-D Notations

In this paper, 𝐱\boldsymbol{\rm x} and 𝐗\boldsymbol{\rm X} in bold typeface are used to represent vectors and matrices, respectively, while scalars are presented in normal typeface, such as xx. We use 𝐗T\boldsymbol{\rm X}^{T} and 𝐗H\boldsymbol{\rm X}^{H} to represent the transpose and conjugate transpose of a matrix 𝐗\boldsymbol{\rm X}, respectively. 𝐈N\boldsymbol{\rm I}_{N} stands for an N×NN\times N identity matrix, and 𝐗F\|\boldsymbol{\rm X}\|_{F} denotes the Frobenius norm of a matrix 𝐗\boldsymbol{\rm X}. Furthermore, 𝔼{}\mathbb{E}\{\cdot\} denotes the expectation operator, and 𝐱𝒞𝒩(𝐦,σ2𝐈)\boldsymbol{\rm x}\sim\mathcal{CN}(\boldsymbol{\rm m},\sigma^{2}\boldsymbol{\rm I}) represents a circularly symmetric complex Gaussian stochastic vector with mean vector 𝐦\boldsymbol{\rm m} and covariance matrix σ2𝐈\sigma^{2}\boldsymbol{\rm I}. Finally, diag(𝐗){{\rm{diag}}\left({\boldsymbol{\rm X}}\right)} denotes a diagonal matrix by keeping only the diagonal elements of matrix 𝐗\boldsymbol{\rm X}.

II System Model

Let us consider a multipair relaying system with KK single-antenna user pairs, denoted as SkS_{k} and DkD_{k}, k=1,,Kk=1,\ldots,K, applying the relay to exchange information with each other. We assume that the direct links between SkS_{k} and DkD_{k} do not exist because of large obstacles or severe shadowing. The large-scale relay is equipped with MM pairs of antennas, namely receive antennas and transmit antennas, with mixed-resolution ADCs and DACs. In the mixed-ADC/DAC architecture, only M0M_{0} pairs of costly high-resolution ADCs and DACs are connected to M0M_{0} relay antennas, while the remaining M1M_{1} (=MM0=M-M_{0}) pairs of less expensive low-resolution ADCs and DACs are connected to M1M_{1} relay antennas. Furthermore, we use κ=ΔM0/M\kappa\mathrel{\mathop{\kern 0.0pt=}\limits^{\Delta}}{{{M_{0}}}\mathord{\left/{\vphantom{{{M_{0}}}M}}\right.\kern-1.2pt}M} (0κ1)\left({0\leq\kappa\leq 1}\right) to denote the fraction of high-resolution ADCs and DACs in the mixed architecture. The low-resolution ADCs bring about severe quantization errors for data reception and the low-resolution DACs lead to obvious signal distortion for data transmission. Therefore, the correlation of the quantization noise is taken into account for the multipair mMIMO relaying system. Furthermore, we assume that the relay operates in half-duplex mode, so it cannot receive and transmit signals simultaneously. Hence, information transmission from SkS_{k} to DkD_{k} is completed in two time slots. In the first time slot, the KK users in the source set SkS_{k} transmit 𝐱SK×1{\boldsymbol{\rm{x}}_{\rm S}}\in\mathbb{C}^{K\times 1} data to the relay independently, and in the next time slot the relay transmits the correlated-quantized signals 𝐱~RM×1\boldsymbol{\tilde{\rm x}}_{\rm R}\in\mathbb{C}^{M\times 1} to KK users in the destination set DkD_{k}. The received signal 𝐲RM×1\boldsymbol{{\rm y}}_{\rm R}\in\mathbb{C}^{M\times 1} at the relay and the received signal 𝐲DK×1\boldsymbol{{\rm y}}_{\rm D}\in\mathbb{C}^{K\times 1} can be respectively given by

𝐲R\displaystyle\boldsymbol{{\rm y}}_{\rm R} =𝐆SR𝐏S1/2𝐱S+𝐧R,\displaystyle={{{\boldsymbol{\rm G}_{\rm{SR}}}\boldsymbol{\rm P}_{\rm S}}^{{1\mathord{\left/{\vphantom{12}}\right.\kern-1.2pt}2}}{\boldsymbol{\rm{x}}_{\rm S}}+{\boldsymbol{\rm{n}}_{\rm{R}}}}, (1)
𝐲D\displaystyle{{\boldsymbol{\rm{y}}}_{D}} =γ𝐆RDT𝐱~R+𝐧D,\displaystyle=\gamma{\boldsymbol{\rm{G}}}_{{\rm{RD}}}^{T}{{\boldsymbol{\rm{\tilde{x}}}}_{{\rm{R}}}}+{{\boldsymbol{\rm{n}}}_{D}}, (2)

where γ\gamma is the normalization factor in order to make the total power at the relay constrained to pRp_{\rm R}, i.e., 𝔼{γ𝐱~R2}=pR{\mathbb{E}}\left\{{{{\left\|{\gamma{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R}}}}\right\|}^{2}}}\right\}={p_{\rm{R}}}. Moreover, 𝐏S\boldsymbol{\rm{P}}_{\rm S} is a diagonal matrix representing the transmit power of the KK source users and its kkth element is given by [𝐏S]kk=pS,k{\left[{{{\boldsymbol{\rm{P}}}_{S}}}\right]_{kk}}={p_{S,k}}. 𝐧R𝒞𝒩(𝟎,𝐈M)\boldsymbol{\rm n}_{\rm R}\sim\mathcal{CN}\left({\boldsymbol{0},{{\boldsymbol{\rm I}}_{M}}}\right), 𝐧D𝒞𝒩(𝟎,𝐈K)\boldsymbol{\rm n}_{\rm D}\sim\mathcal{CN}\left({\boldsymbol{0},{{\boldsymbol{\rm I}}_{K}}}\right) denote the additive white Gaussian noise (AWGN) matrix with independently and identically distributed (i.i.d.) components following the distribution 𝒞𝒩(0,1)\mathcal{CN}\left(0,1\right). We further follow the general assumption that the transmit signal vector 𝐱𝐒\boldsymbol{\rm x_{\rm S}} is Gaussian distributed. The matrices 𝐆SR=[𝐠SR,1,,𝐠SR,K]{{\boldsymbol{\rm{G}}}_{{\rm{SR}}}}=\left[{{\boldsymbol{\rm g}_{{\rm{SR}},1}},\ldots,{\boldsymbol{\rm g}_{{\rm{SR}},K}}}\right] and 𝐆RDT=[𝐠RD,1,,𝐠RD,K]T{{\boldsymbol{\rm{G}}}_{{\rm{RD}}}^{T}}=\left[{{\boldsymbol{\rm g}_{{\rm{RD}},1}},\ldots,{\boldsymbol{\rm g}_{{\rm{RD}},K}}}\right]^{T} refer to the Rayleigh fading channels from the KK sources to the relay with 𝐠SR,k𝒞𝒩(𝟎,βSR,k𝐈M){\boldsymbol{\rm g}_{{\rm{SR}},k}}\sim\mathcal{CN}\left({\boldsymbol{0},{\beta_{{\rm{SR}},k}}{{\boldsymbol{\rm I}}_{M}}}\right) and the channels from the relay to the KK destinations with 𝐠RD,k𝒞𝒩(𝟎,βRD,k𝐈M){\boldsymbol{\rm g}_{{\rm{RD}},k}}\sim\mathcal{CN}\left({\boldsymbol{0},{\beta_{{\rm{RD}},k}{\boldsymbol{\rm I}}_{M}}}\right), respectively. The terms βSR,k\beta_{\rm{SR},k} and βRD,k\beta_{\rm{RD},k} stand for the large-scale fading and are assumed to be known at the relay.

Furthermore, we define 𝐆SR0\boldsymbol{\rm G}_{\rm{SR0}} as the M0×KM_{0}\times K channel matrix from the KK sources to the M0M_{0} relay antennas connected with high-resolution ADCs, and 𝐆SR1\boldsymbol{\rm G}_{\rm{SR1}} as the M1×KM_{1}\times K channel matrix from the KK sources to the remained M1M_{1} relay antennas connected with low-resolution ADCs. Therefore, we have

𝐆SR=[𝐆SR0𝐆SR1].{{\boldsymbol{\rm{G}}}_{{\rm{SR}}}}=\left[{\begin{array}[]{*{20}{c}}{{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}}\\ {{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}}\end{array}}\right]. (3)

Similarly, we can also define 𝐆RD0T\boldsymbol{\rm G}_{\rm{RD0}}^{T} as the K×M0K\times M_{0} channel matrix from the M0M_{0} relay transmit antennas connected with high-resolution DACs to the KK destinations, and 𝐆RD1T\boldsymbol{\rm G}_{\rm{RD1}}^{T} as the M1×KM_{1}\times K channel matrix from the M1M_{1} relay transmit antennas connected with low-resolution DACs to the KK destinations. Then, 𝐆RD\boldsymbol{\rm G}_{\rm{RD}} can be expressed as

𝐆RD=[𝐆RD0𝐆RD1].{{\boldsymbol{\rm{G}}}_{{\rm{RD}}}}=\left[{\begin{array}[]{*{20}{c}}{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}}}\\ {{{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}}}\end{array}}\right]. (4)

With the help of (3) and (4), (1) and (2) can be rewritten as

𝐲R\displaystyle\boldsymbol{{\rm y}}_{\rm R} =[𝐲R0𝐲R1]=[𝐆SR0𝐏S1/2𝐱S+𝐧R0𝐆SR1𝐏S1/2𝐱S+𝐧R1],\displaystyle=\left[{\begin{array}[]{*{20}{c}}\boldsymbol{{\rm y}}_{\rm{R0}}\\ \boldsymbol{{\rm y}}_{\rm{R1}}\end{array}}\right]=\left[{\begin{array}[]{*{20}{c}}{{{\boldsymbol{\rm G}_{\rm{SR0}}}\boldsymbol{\rm P}_{\rm S}}^{{1\mathord{\left/{\vphantom{12}}\right.\kern-1.2pt}2}}{\boldsymbol{\rm{x}}_{\rm S}}+{\boldsymbol{\rm{n}}_{\rm{R0}}}}\\ {{\boldsymbol{\rm G}_{\rm{SR1}}}\boldsymbol{\rm P}_{\rm S}}^{{1\mathord{\left/{\vphantom{12}}\right.\kern-1.2pt}2}}{\boldsymbol{\rm{x}}_{\rm S}}+{\boldsymbol{\rm{n}}_{\rm{R1}}}\end{array}}\right],
𝐲D\displaystyle{{\boldsymbol{\rm{y}}}_{D}} =[𝐲D0𝐲D1]=γ𝐆RD0T𝐱~R0+γ𝐆RD1T𝐱~R1+𝐧D,\displaystyle=\left[{\begin{array}[]{*{20}{c}}{{{\boldsymbol{\rm{y}}}_{{\rm{D0}}}}}\\ {{{\boldsymbol{\rm{y}}}_{{\rm{D1}}}}}\end{array}}\right]=\gamma{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}{{\boldsymbol{\rm{\tilde{x}}}}_{{\rm{R0}}}}+\gamma{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{T}{{\boldsymbol{\rm{\tilde{x}}}}_{{\rm{R1}}}}+{{\boldsymbol{\rm{n}}}_{D}},

where 𝐲R0{{\boldsymbol{\rm{y}}}_{\rm{R0}}} denotes the first M0M_{0} rows of the overall received signals vector 𝐲R{\boldsymbol{\rm{y}}}_{\rm{R}}, and 𝐲R1{{\boldsymbol{\rm{y}}}_{\rm{R1}}} denotes the rest M1M_{1} rows of 𝐲R{{\boldsymbol{\rm{y}}}_{\rm{R}}}. Without loss of generality, the notations 𝐲D0{{\boldsymbol{\rm{y}}}_{\rm{D0}}}, 𝐲D1{{\boldsymbol{\rm{y}}}_{\rm{D1}}}, 𝐧R0{\boldsymbol{\rm{n}}_{\rm{R0}}}, 𝐧R1{\boldsymbol{\rm{n}}_{\rm{R1}}}, 𝐱~R0{{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}} and 𝐱~R1{{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R1}}}}} can also be explained in a similar way.

II-A Quantization with Mixed-Resolution ADCs

For the mixed-ADC architecture, the quantized received signal at the relay can be written as

𝐲~R=[𝐲~R0𝐲~R1]=[𝐲R0(𝐲R1)],{\boldsymbol{\tilde{\rm y}}_{\rm R}}=\left[{\begin{array}[]{*{20}{c}}{\boldsymbol{\tilde{\rm y}}_{\rm R0}}\\ {\boldsymbol{\tilde{\rm y}}_{\rm R1}}\end{array}}\right]=\left[{\begin{array}[]{*{20}{c}}{\boldsymbol{\rm y}_{\rm R0}}\\ {\mathbb{Q}\left({{\boldsymbol{\rm y}_{\rm R1}}}\right)}\end{array}}\right], (11)

where ()\mathbb{Q}\left(\cdot\right) is the scalar quantization function, 𝐲~R0{\boldsymbol{\tilde{\rm y}}_{\rm R0}} denotes the quantized received signals at the output of M0M_{0} high-resolution ADCs, and 𝐲~R1{\boldsymbol{\tilde{\rm y}}_{\rm R1}} is the quantized received signals at the output of M1M_{1} low-resolution ADCs. According to the AQNM [30, Eq. (1)], the quantization operation can be expressed as

𝐲~R1=(𝐲R1)=α𝐲R1+𝐧qa,{{\boldsymbol{\rm{\tilde{y}}}}_{{\rm{R1}}}}=\mathbb{Q}\left({{{\boldsymbol{\rm{y}}}_{{\rm{R1}}}}}\right)=\alpha{{\boldsymbol{\rm{y}}}_{{\rm{R1}}}}+{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}, (12)

where 𝐧qa{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}} refers to the additive Gaussian quantization noise vector which is uncorrelated with 𝐲R1{{\boldsymbol{\rm{y}}}_{{\rm{R1}}}}, and α\alpha denotes a linear gain given by [31, Eq. (13)]

α=1ρ=1𝔼{𝐲~R1𝐲R12}/𝔼{𝐲~R12},\alpha=1-\rho=1-{{\mathbb{E}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{y}}}}}_{\rm R1}}-{{\boldsymbol{\rm{y}}}_{\rm R1}}}\right\|}^{2}}}\right\}}}/{{\mathbb{E}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{y}}}}}_{\rm R1}}}\right\|}^{2}}}\right\}}}, (13)

with ρ\rho as the distortion factor of the low-resolution ADCs. The exact values of ρ\rho are given in Table I with respect to different resolution bits [32]. For large quantization bits (e.g., b>5b>5), the distortion factor ρ\rho can be approximated as ρπ3222b\rho\approx\frac{{\pi\sqrt{3}}}{2}{2^{-2b}} [32]. With the help of (II), (12) and (13), the covariance matrix of 𝐧qa{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}} is expressed as

𝐑𝐧qa=αρdiag(𝐆SR1𝐏S𝐆SR1H+𝐈M1).\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}=\alpha\rho{\rm{diag}}\left({\boldsymbol{\rm G}_{\rm SR1}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}+{\boldsymbol{{\rm I}}_{{M_{1}}}}}\right). (14)

Moreover, (11) can be rewritten as

𝐲~R=[𝐲~R0𝐲~R1]=[𝐆SR0𝐏S1/2𝐱S+𝐧R0α𝐆SR1𝐏S1/2𝐱S+α𝐧R1+𝐧qa].{\boldsymbol{\tilde{\rm y}}_{\rm R}}=\left[{\begin{array}[]{*{20}{c}}{\boldsymbol{\tilde{\rm y}}_{\rm R0}}\\ {\boldsymbol{\tilde{\rm y}}_{\rm R1}}\end{array}}\right]=\left[{\begin{array}[]{*{20}{c}}{{{\boldsymbol{\rm G}_{\rm{SR0}}}\boldsymbol{\rm P}_{\rm S}}^{{1\mathord{\left/{\vphantom{12}}\right.\kern-1.2pt}2}}{\boldsymbol{\rm{x}}_{\rm S}}\!+\!{\boldsymbol{\rm{n}}_{\rm{R0}}}}\\ \alpha{{{\boldsymbol{\rm G}_{\rm{SR1}}}\boldsymbol{\rm P}_{\rm S}}^{{1\mathord{\left/{\vphantom{12}}\right.\kern-1.2pt}2}}{\boldsymbol{\rm{x}}_{\rm S}}\!+\!\alpha{\boldsymbol{\rm{n}}_{\rm{R1}}}\!+\!{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}}\end{array}}\right]. (15)
TABLE I: Distortion Factors For Different Quantization Bits.
bb 1 2 3 4 5
ρ\rho 0.3634 0.1175 0.03454 0.009497 0.002499

II-B Maximum Ratio (MR) Processing at the Relay

We assume that the relay adopts a simple amplify-and-forward (AF)11 1 The AF protocol is considered herein due to its lower implementation complexity compared with the decode-and-forward (DF) protocol [5]. protocol to process the quantized received signals, yielding

𝐱R=𝐖𝐲~R,{\boldsymbol{\rm x}_{\rm R}}=\boldsymbol{\rm W}{\boldsymbol{\tilde{\rm y}}_{\rm R}}, (16)

where 𝐖=𝐆RD𝐆SRH\boldsymbol{\rm W}={\boldsymbol{\rm{G}}}_{{\rm{RD}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR}}}^{H} denotes the MR processing. The MR processing is used at the relay due to its low-complexity, being suitable for the low-cost multipair mMIMO relaying system. Furthermore, according to some research, the MR processing can achieve similar performance as zero-forcing receiver/zero-forcing transmission (ZFR/ZFT) or minimum mean square error (MMSE). By applying (3) and (4), (16) is rewritten as

𝐱R=[𝐆RD0𝐆SR0H𝐲~R0+𝐆RD0𝐆SR1H𝐲~R1𝐆RD1𝐆SR0H𝐲~R0+𝐆RD1𝐆SR1H𝐲~R1].\displaystyle\boldsymbol{{\rm x}}_{\rm R}=\left[{\begin{array}[]{*{20}{c}}{{\boldsymbol{\rm G}_{\rm RD0}^{*}}{\boldsymbol{\rm G}_{\rm SR0}^{H}}{\boldsymbol{\tilde{\rm y}}_{\rm R0}}\!+\!{\boldsymbol{\rm G}_{\rm RD0}^{*}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}{\boldsymbol{\tilde{\rm y}}_{\rm R1}}}\\ {{\boldsymbol{\rm G}_{\rm RD1}^{*}}{\boldsymbol{\rm G}_{\rm SR0}^{H}}{\boldsymbol{\tilde{\rm y}}_{\rm R0}}\!+\!{\boldsymbol{\rm G}_{\rm RD1}^{*}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}{\boldsymbol{\tilde{\rm y}}_{\rm R1}}}\end{array}}\right]. (17)

II-C Quantization with Mixed-Resolution DACs

For simplicity, we assume that the DACs and ADCs have the same resolution. The analysis method can be extended to arbitrary resolution cases. With mixed-DAC architecture at the transmitter, the mixed-ADC/DAC signals from the relay’s transmit antennas can be expressed as

𝐱~R=[𝐱~R0𝐱~R1]=[𝐱R0(𝐱R1)]=[𝐱R0α𝐱R1+𝐧qD],{\boldsymbol{\tilde{\rm x}}_{\rm R}}=\left[{\begin{array}[]{*{20}{c}}{\boldsymbol{\tilde{\rm x}}_{\rm R0}}\\ {\boldsymbol{\tilde{\rm x}}_{\rm R1}}\end{array}}\right]=\left[{\begin{array}[]{*{20}{c}}{\boldsymbol{\rm x}_{\rm R0}}\\ {\mathbb{Q}\left({{\boldsymbol{\rm x}_{\rm R1}}}\right)}\end{array}}\right]=\left[{\begin{array}[]{*{20}{c}}{\boldsymbol{\rm x}_{\rm R0}}\\ \alpha{{\boldsymbol{\rm{x}}}_{{\rm{R1}}}}+{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}\end{array}}\right], (18)

where α\alpha is the distortion factor of the low-resolution DACs, and 𝐧qD{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}} denotes the quantization noise of low-resolution DACs, which is uncorrelated with 𝐱R1{{\boldsymbol{\rm{x}}}_{{\rm{R1}}}}. Note that 𝐱~R0=𝐱R0{\boldsymbol{\tilde{\rm x}}_{\rm R0}}={\boldsymbol{\rm x}_{\rm R0}} because of using high-resolution DACs. Similar as for (14), we can derive the covariance matrix of 𝐧qD{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}} as

𝐑𝐧qD=αρdiag(𝐑𝐱R1),\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}=\alpha\rho{\rm{diag}}\left(\boldsymbol{\rm R}_{{\boldsymbol{\rm{x}}}_{{{\rm{R1}}}}}\right), (19)

where 𝐑𝐱R1\boldsymbol{\rm R}_{{\boldsymbol{\rm{x}}}_{{{\rm{R1}}}}} is the covariance matrix of 𝐱R1{{\boldsymbol{\rm{x}}}_{{\rm{R1}}}}. According to (17), 𝐑𝐱R1\boldsymbol{\rm R}_{{\boldsymbol{\rm{x}}}_{{{\rm{R1}}}}} can be written as

𝐑𝐱R1\displaystyle{\boldsymbol{\rm R}_{{\boldsymbol{\rm{x}}}_{{{\rm{R1}}}}}} =𝐆RD1𝐆SR0H𝐑𝐲~R0𝐲~R0𝐆SR0𝐆RD1T\displaystyle={\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{T} (20)
+𝐆RD1𝐆SR0H𝐑𝐲~R0𝐲~R1𝐆SR1𝐆RD1T\displaystyle+{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{T}
+𝐆RD1𝐆SR1H𝐑𝐲~R1𝐲~R0𝐆SR0𝐆RD1T\displaystyle+{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{T}
+𝐆RD1𝐆SR1H𝐑𝐲~R1𝐲~R1𝐆SR1𝐆RD1T,\displaystyle+{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{T},

where

𝐑𝐲~R0𝐲~R0\displaystyle{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}} =𝐆SR0𝐏S𝐆SR0H+𝐈M0,\displaystyle={\boldsymbol{\rm G}_{\rm SR0}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR0}^{H}}+{\boldsymbol{{\rm I}}_{{M_{0}}}}}, (21)
𝐑𝐲~R0𝐲~R1\displaystyle{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}} =α𝐆SR0𝐏S𝐆SR1H,\displaystyle=\alpha{\boldsymbol{\rm G}_{\rm SR0}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}}, (22)
𝐑𝐲~R1𝐲~R0\displaystyle{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}} =α𝐆SR1𝐏S𝐆SR0H,\displaystyle=\alpha{\boldsymbol{\rm G}_{\rm SR1}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR0}^{H}}}, (23)
𝐑𝐲~R1𝐲~R1\displaystyle{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}} =α2(𝐆SR1𝐏S𝐆SR1H+𝐈M1)\displaystyle=\alpha^{2}\left({\boldsymbol{\rm G}_{\rm SR1}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}+{\boldsymbol{{\rm I}}_{{M_{1}}}}}\right)
+αρdiag(𝐆SR1𝐏S𝐆SR1H+𝐈M1).\displaystyle+\alpha\rho{\rm{diag}}\left({\boldsymbol{\rm G}_{\rm SR1}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}+{\boldsymbol{{\rm I}}_{{M_{1}}}}}\right). (24)

Consequently, the normalization factor γ\gamma can be expressed as

γ=pR/𝔼{𝐱~R2}.\gamma=\sqrt{{{{p_{R}}}}/{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R}}}}\right\|}^{2}}}\right\}}}. (25)
Lemma 1.

For mixed-ADC/DAC multipair mMIMO relaying systems, the expectation of the total transmit power at the relay can be expressed as

𝔼{𝐱~R2}=μ(M0+αM1),\displaystyle{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R}}}}\right\|}^{2}}}\right\}}=\mu\left({{M_{0}}+\alpha{M_{1}}}\right), (26)

where μ\mu is given by

μ=α(1α)M1k=1KpS,kβSR,k2βRD,k.\displaystyle\mu=\alpha\left({1-\alpha}\right){M_{1}}\sum\limits_{k=1}^{K}{{p_{{S},k}}\beta_{{\rm{SR}},k}^{2}{\beta_{{\rm{RD}},k}}}. (27)
Proof:

Please refer to Appendix A. ∎

Substituting (26) into (25), the normalization factor γ\gamma can be obtained as

γ=pR/μ(M0+αM1).\gamma=\sqrt{{{{p_{R}}}}/{{\mu\left({{M_{0}}+\alpha{M_{1}}}\right)}}}. (28)

With the normalization factor in hand, we can derive the achievable rate in the following section.

III Achievable Rate Analysis

III-A Exact Achievable Rate Analysis

It is assumed that the destination DkD_{k} applies only statistical CSI to decode the signal, due to the reason that instantaneous CSI leads to excessive high computational complexity for large antenna arrays in a practical mMIMO system. Combining (1), (2), (15), (17) and (18), the received signal at the destination DkD_{k} can be expressed as

yD,k=pS,k𝔼{Tk,k}xS,kdesired signal+n~D,keffective noise,{y_{{D},k}}=\underbrace{\sqrt{{p_{{S},k}}}{\mathbb{E}}\left\{{{T_{k,k}}}\right\}{x_{{S},k}}}_{\text{desired signal}}+\underbrace{\tilde{n}_{{D},k}}_{\text{effective noise}}, (29)

where

Ti,j=γ𝐠RD0,iT𝐆RD0𝐆SR0H𝐠SR0,j+γα𝐠RD0,iT𝐆RD0𝐆SR1H𝐠SR1,j\displaystyle{\displaystyle T_{i,j}}\!=\!\gamma{\boldsymbol{\rm g}_{{\rm{RD0}},i}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR0}},j}}}\!+\!\gamma\alpha\boldsymbol{\rm g}_{{\rm{RD0}},i}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR1}},j}}
+γα𝐠RD1,iT𝐆RD1𝐆SR0H𝐠SR0,j+γα2𝐠RD1,iT𝐆RD1𝐆SR1H𝐠SR1,j,\displaystyle+\gamma\alpha\boldsymbol{\rm g}_{{\rm{RD1}},i}^{T}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR0}},j}}\!+\!\gamma{{\alpha^{2}}\boldsymbol{\rm g}_{{\rm{RD1}},i}^{T}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR1}},j}}},

and

n~D,k=pS,k(Tk,k𝔼{Tk,k})xS,kestimation error+ikTk,ixS,iinter-pair interference\displaystyle{{\tilde{n}}_{D,k}}=\underbrace{\sqrt{{p_{{S},k}}}\left({{T_{k,k}}-{\mathbb{E}}\left\{{{T_{k,k}}}\right\}}\right){x_{{S},k}}}_{\text{estimation error}}+\underbrace{\sum\limits_{i\neq k}{{T_{k,i}}{x_{{S},i}}}}_{\text{inter-pair interference}}
+γ(𝐠RD0,kT𝐆RD0𝐆SR0H𝐧R0+α𝐠RD1,kT𝐆RD1𝐆SR0H𝐧R0)noise of high-resolution quantization at the relay\displaystyle+\underbrace{\gamma\left({{\boldsymbol{\rm{g}}}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{n}}}_{{\rm{R0}}}}+\alpha{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{T}{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{n}}}_{{\rm{R0}}}}}\right)}_{\text{noise of high-resolution quantization at the relay}}
+γ(α𝐠RD0,kT𝐆RD0𝐆SR1H𝐧R1+α2𝐠RD1,kT𝐆RD1𝐆SR1H𝐧R1)noise of low-resolution quantization at the relay\displaystyle\!+\!\underbrace{\gamma\left({\alpha{\boldsymbol{\rm{g}}}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{n}}}_{{\rm{R1}}}}\!+\!{\alpha^{2}}{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{T}{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{n}}}_{{\rm{R1}}}}}\right)}_{\text{noise of low-resolution quantization at the relay}}
+γ(𝐠RD0,kT𝐆RD0𝐆SR1H𝐧qa+α𝐠RD1,kT𝐆RD1𝐆SR1H𝐧qa)quantization noise of ADCs\displaystyle+\underbrace{\gamma\left({{\boldsymbol{\rm{g}}}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}+\alpha{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{T}{\boldsymbol{\rm{G}}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}}\right)}_{\text{quantization noise of ADCs}}
+γ𝐠RD1,kT𝐧qDquantization noise of DACs+nD,knoise,\displaystyle+\underbrace{\gamma{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{T}{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}}_{\text{quantization noise of DACs}}+\underbrace{{{{n}}_{{D},k}}}_{\text{noise}}, (30)

where 𝐧D,k{{\boldsymbol{\rm{n}}}_{{D},k}} is the kkth element of the vector 𝐧D{{\boldsymbol{\rm{n}}}_{{D}}}. We can derive the signal-to-interference-plus-noise ratio (SINR) expression by using [33, Eq. (18)]. Since the “desired signal” and the “effective noise” in (29) are uncorrelated, the exact achievable rate for the kk-th destination is given in Theorem 1.

Theorem 1.

For mixed-ADC/DAC multipair mMIMO relaying systems and using the capacity bound in [34], the exact closed-form achievable rate of the kk-th destination is given as

Rk=τc2τp2τclog2(+AkBk+Ck+Dk+Ek+Fk+Gk+1),{R_{k}}\!=\!\frac{{{\tau_{c}}\!-\!2{\tau_{p}}}}{{2{\tau_{c}}}}{\log_{2}}\left({1\!+\!\frac{{{A_{k}}}}{{{B_{k}}\!+\!{C_{k}}\!+\!{D_{k}}\!+\!{E_{k}}\!+\!{F_{k}}\!+\!{G_{k}}\!+\!1}}}\right), (31)

where τc\tau_{c} denotes the length (in symbols) of each coherence interval, τp\tau_{p} represents the length of the mutually orthogonal pilot sequences, and

Ak=pS,kγ2(M0+αM1)4βSR,k2βRD,k2,\displaystyle{A_{k}}={p_{{S},k}}{\gamma^{2}}{\left({{M_{0}}+\alpha{M_{1}}}\right)^{4}}\beta_{{\rm{SR}},k}^{2}\beta_{{\rm{RD}},k}^{2}, (32)
Bk=pS,kγ2(M0+α2M1)βSR,kβRD,k[2(M0+αM1)2\displaystyle{B_{k}}={p_{{S},k}}{\gamma^{2}}\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right){\beta_{{\rm{SR}},k}}{\beta_{{\rm{RD}},k}}\big[2{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}^{2}}
×βSR,kβRD,k+(M0+α2M1)m=1KβSR,mβRD,m],\displaystyle\times{\beta_{{\rm{SR}},k}}{\beta_{{\rm{RD}},k}}\!+\!\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right)\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}\big], (33)
Ck=γ2(M0+α2M1)ikpS,i\displaystyle{C_{k}}={\gamma^{2}}\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right)\sum\limits_{i\neq k}{{p_{{\rm S},i}}}
×[(M0+αM1)2(βSR,kβRD,k2βSR,i+βRD,kβSR,i2βRD,i)\displaystyle\times\Big[{{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{2}}\left({\beta_{{\rm SR},k}}\beta_{{\rm RD},k}^{2}{\beta_{{\rm SR},i}}+{\beta_{{\rm RD},k}}\beta_{{\rm SR},i}^{2}{\beta_{{\rm RD},i}}\right)}
+(M0+α2M1)βRD,kβSR,im=1KβSR,mβRD,m],\displaystyle+\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right){\beta_{{\rm RD},k}}{\beta_{{\rm SR},i}}\sum\limits_{m=1}^{K}{{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}\Big], (34)
Dk=γ2M0(M0+αM1)2βSR,kβRD,k2\displaystyle{D_{k}}={\gamma^{2}}{M_{0}}{\left({{M_{0}}+\alpha{M_{1}}}\right)^{2}}{\beta_{{\rm SR},k}}\beta_{{\rm RD},k}^{2}
+γ2M0(M0+α2M1)βRD,km=1KβSR,mβRD,m,\displaystyle+{\gamma^{2}}{M_{0}}\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right){\beta_{{\rm RD},k}}\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}, (35)
Ek=γ2α2M1(M0+αM1)2βSR,kβRD,k2\displaystyle{E_{k}}={\gamma^{2}}{\alpha^{2}}{M_{1}}{\left({{M_{0}}+\alpha{M_{1}}}\right)^{2}}{\beta_{{\rm SR},k}}\beta_{{\rm RD},k}^{2}
+γ2α2M1(M0+α2M1)βRD,km=1KβSR,mβRD,m,\displaystyle+{\gamma^{2}}{\alpha^{2}}{M_{1}}\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right){\beta_{{\rm RD},k}}\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}, (36)
Fk=αργ2M1βRD,k{(M0+α2M1)m=1K[βSR,mβRD,m\displaystyle{F_{k}}=\alpha\rho{\gamma^{2}}{M_{1}}{\beta_{{\rm RD},k}}\left\{\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right)\sum\limits_{m=1}^{K}\Big[{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}\right.
×(i=1KpS,iβSR,i+pS,mβSR,m+1)]+(M0+αM1)2\displaystyle\times\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm SR},i}}}+{p_{{\rm S},m}}{\beta_{{\rm SR},m}}+1}\right)\Big]+{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{2}}
×βSR,kβRD,k(i=1KpS,iβSR,i+pS,kβSR,k+1)},\displaystyle\times\left.{{\beta_{{\rm SR},k}}{\beta_{{\rm RD},k}}\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm SR},i}}}+{p_{{\rm S},k}}{\beta_{{\rm SR},k}}+1}\right)}\right\}, (37)
Gk=αργ2M1(M0+αM1)βRD,k{βSR,kβRD,k\displaystyle{G_{k}}=\alpha\rho{\gamma^{2}}{M_{1}}\left({{M_{0}}+\alpha{M_{1}}}\right){\beta_{{\rm RD},k}}\Big\{{{\beta_{{\rm SR},k}}{\beta_{{\rm RD},k}}}
×(i=1KpS,iβSR,i++(M0+αM1)pS,kβSR,k)+m=1KβSR,m\displaystyle\times\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm SR},i}}}\!+\!1\!+\!\left({{M_{0}}\!+\!\alpha{M_{1}}}\right){p_{{\rm S},k}}{\beta_{{\rm SR},k}}}\right)\!+\!\sum\limits_{m=1}^{K}{\beta_{{\rm SR},m}}
×βRD,m(i=1KpS,iβSR,i+1+(M0+αM1)pS,mβSR,m)}\displaystyle\times{\beta_{{\rm RD},m}}\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm SR},i}}}+1+\left({{M_{0}}+\alpha{M_{1}}}\right){p_{{\rm S},m}}{\beta_{{\rm SR},m}}}\right)\Big\}
+α2ρ2γ2M12βRD,k(i=1KpS,iβSR,i2βRD,i+pS,kβSR,k2βRD,k).\displaystyle\!+\!{\alpha^{2}}{\rho^{2}}{\gamma^{2}}M_{1}^{2}{\beta_{{\rm RD},k}}\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}\beta_{{\rm SR},i}^{2}{\beta_{{\rm RD},i}}}\!+\!{p_{{\rm S},k}}\beta_{{\rm SR},k}^{2}{\beta_{{\rm RD},k}}}\right). (38)
Proof:

Please refer to Appendix B. ∎

 


aki={βSR,iβSR,k(M0+αM1)(+αρM1(M0+αM1)2+βSR,iβRD,iβSR,kβRD,k[+αρM1(M0+αM1)2]+1(M0+αM1)m=1KβSR,mβRD,mβSR,kβRD,k),ik1(M0+αM1)(+αρM1(M0+αM1)[+2(M0+αM1)+αρM1(M0+αM1)2]+1(M0+αM1)m=1KβSR,mβRD,mβSR,kβRD,k),i=k\displaystyle{a_{ki}}=\left\{{\begin{array}[]{ll}\frac{\frac{{{\beta_{{\rm SR},i}}}}{{{\beta_{{\rm SR},k}}}}}{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}}\left({1\!+\!\frac{{\alpha\rho{M_{1}}}}{{{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}^{2}}}}\!+\!\frac{{{\beta_{{\rm SR},i}}{\beta_{{\rm RD},i}}}}{{{\beta_{{\rm SR},k}}{\beta_{{\rm RD},k}}}}\left[{1\!+\!\frac{{\alpha\rho{M_{1}}}}{{{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}^{2}}}}}\right]\!+\!\frac{1}{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}}\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{{\beta_{{\rm SR},k}}{\beta_{{\rm RD},k}}}}}\right),&{i\neq k}\\ \frac{1}{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}}\left({2\!+\!\frac{{\alpha\rho{M_{1}}}}{{\left({{M_{0}}+\alpha{M_{1}}}\right)}}\left[{2\!+\!\frac{2}{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}}\!+\!\frac{{\alpha\rho{M_{1}}}}{{{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}^{2}}}}}\right]+\frac{1}{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}}\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{{\beta_{{\rm SR},k}}{\beta_{{\rm RD},k}}}}}\right),&{i=k}\end{array}}\right. (41)

 

R^k,Rkτc2τp2τclog2[1+ESER(α+ρκ)2ES(α+ρκ)i=1KβSR,i2βRD,iβSR,k2βRD,k2+m=1KβSR,mβRD,mβSR,k2βRD,k2+ER(α+ρκ)βSR,k].{\hat{R}_{k}},{R_{k}}\to\frac{{{\tau_{c}}-2{\tau_{p}}}}{{2{\tau_{c}}}}{\log_{2}}\Bigg[{1+\frac{{{E_{\rm S}}{E_{\rm R}}{{\left({\alpha+\rho\kappa}\right)}^{2}}}}{{{E_{\rm S}}\left({\alpha+\rho\kappa}\right)\sum\limits_{i=1}^{K}{\frac{{\beta_{{\rm SR},i}^{2}{\beta_{{\rm RD},i}}}}{{\beta_{{\rm SR},k}^{2}\beta_{{\rm RD},k}^{2}}}}+\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{\beta_{{\rm SR},k}^{2}\beta_{{\rm RD},k}^{2}}}+\frac{{{E_{\rm R}}\left({\alpha+\rho\kappa}\right)}}{{{\beta_{{\rm SR},k}}}}}}}\Bigg]. (50)

With the help of (28) and after some simplifications, we can derive the compact expression for the sum achievable rate as

R=τc2τp2τck=1Klog2(1+νk),{R}=\frac{{{\tau_{c}}-2{\tau_{p}}}}{{2{\tau_{c}}}}\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+{\nu_{k}}}\right)}, (39)

where

νk\displaystyle{\nu_{k}} =pS,k/ξk,\displaystyle={{{p_{{\rm S},k}}}}/{{{\xi_{k}}}}, (40)
ξk\displaystyle{\xi_{k}} =i=1KpS,iaki+pR1(i=1KpS,ibki+ck)+dk,\displaystyle=\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{a_{ki}}}+p_{\rm R}^{-1}\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{b_{ki}}}+{c_{k}}}\right)+{d_{k}}, (41)
bki\displaystyle{b_{ki}} =1(M0+αM1)βSR,i2βRD,iβSR,k2βRD,k2[1+αρM1(M0+αM1)2]\displaystyle=\frac{1}{{\left({{M_{0}}+\alpha{M_{1}}}\right)}}\frac{{\beta_{{\rm SR},i}^{2}{\beta_{{\rm RD},i}}}}{{\beta_{{\rm SR},k}^{2}\beta_{{\rm RD},k}^{2}}}\left[{1+\frac{{\alpha\rho{M_{1}}}}{{{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{2}}}}}\right]
+1(M0+αM1)2βSR,im=1KβSR,mβRD,mβSR,k2βRD,k2,\displaystyle+\frac{1}{{{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{2}}}}\frac{{{\beta_{{\rm SR},i}}\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{\beta_{{\rm SR},k}^{2}\beta_{{\rm RD},k}^{2}}}, (42)
ck\displaystyle{c_{k}} =1(M0+αM1)2m=1KβSR,mβRD,mβSR,k2βRD,k2,\displaystyle=\frac{1}{{{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{2}}}}\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{\beta_{{\rm SR},k}^{2}\beta_{{\rm RD},k}^{2}}}, (43)
dk\displaystyle{d_{k}} =1(M0+αM1)1βSR,k+αρM1(M0+αM1)31βSR,k\displaystyle=\frac{1}{{\left({{M_{0}}+\alpha{M_{1}}}\right)}}\frac{1}{{{\beta_{{\rm SR},k}}}}+\frac{{\alpha\rho{M_{1}}}}{{{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{3}}}}\frac{1}{{{\beta_{{\rm SR},k}}}}
+1(M0+αM1)2m=1KβSR,mβRD,mβSR,k2βRD,k,\displaystyle+\frac{1}{{{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{2}}}}\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{\beta_{{\rm SR},k}^{2}{\beta_{{\rm RD},k}}}}, (44)

with akia_{ki} given by (41) at the bottom of next page.

III-B Asymptotic Analysis

Note that Fk{F_{k}} and Gk{G_{k}} are respectively derived by using exact values of the covariance matrices 𝐑𝐧qa\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}} and 𝐑𝐧qD\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}, which make results in (37) and (38) cumbersome. In order to provide more insights into the effect of various parameters on the achievable rate, we consider a large number of antennas and use the law of large numbers. The covariance matrix 𝐑𝐧qa\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}} in (14) is given by

𝐑𝐧qaαρ(k=1KpS,kβSR,k+1)𝐈M1.\displaystyle\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}\approx\alpha\rho\left({\sum\limits_{k=1}^{K}{{p_{{\rm S},k}}{\beta_{{\rm SR},k}}}+1}\right){\boldsymbol{\rm I}_{{M_{1}}}}. (42)

Similarly, the covariance matrix 𝐑𝐧qD\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}} is approximated as

𝐑𝐧qD\displaystyle\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}} 𝔼{𝐑𝐧qD}=αρμ𝐈M1,\displaystyle\approx{\mathbb{E}}\left\{\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}\right\}=\alpha\rho\mu{\boldsymbol{\rm I}_{{M_{1}}}}, (43)

where μ\mu is given by (27). Note that the expression of the approximate 𝐑𝐧qD\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}} is obtained in (94), in Appendix A. Based on the aforementioned discussion, we can derive a more concise close-form approximation for the achievable rate by simplifying Fk{F_{k}} and Gk{G_{k}} as in (42) and (43), respectively. Similar to Theorem 1, we can derive the approximate achievable rate R^k\hat{R}_{k} by calculating the power expectations of the signal and interference as shown in the following theorem.

Theorem 2.

For mixed-ADC/DAC multipair mMIMO relaying systems, the approximate achievable rate of the kk-th destination is

R^k=τc2τp2τclog2(+AkBk+Ck+Dk+Ek+F^k+G^k+1),{\hat{R}_{k}}\!=\!\frac{{{\tau_{c}}\!-\!2{\tau_{p}}}}{{2{\tau_{c}}}}{\log_{2}}\left({1\!+\!\frac{{{A_{k}}}}{{{B_{k}}\!+\!{C_{k}}\!+\!{D_{k}}\!+\!{E_{k}}\!+\!{\hat{F}_{k}}\!+\!{\hat{G}_{k}}\!+\!1}}}\right), (44)

where G^k=γ2αρμM1βRD,k{\hat{G}_{k}}={\gamma^{2}}\alpha\rho\mu{M_{1}}{\beta_{{\rm RD},k}},

F^k=γ2αρ(k=1KpS,kβSR,k+1)M1βRD,k((M0+αM1)2CLOSE\displaystyle{\hat{F}_{k}}={\gamma^{2}}\alpha\rho\left({\sum\limits_{k=1}^{K}{{p_{{\rm S},k}}{\beta_{{\rm SR},k}}}+1}\right){M_{1}}{\beta_{{\rm RD},k}}\Big({\left({{M_{0}}+\alpha{M_{1}}}\right)^{2}}
×βSR,kβRD,k+(M0+α2M1)m=1KβSR,mβRD,m),\displaystyle\times{\beta_{{\rm SR},k}}\beta_{{\rm RD},k}+\left({{M_{0}}+{\alpha^{2}}{M_{1}}}\right)\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}\Big), (45)

and AkA_{k}, BkB_{k}, CkC_{k}, DkD_{k} and EkE_{k} are given by (32), (33), (34), (35), and (36), respectively.

Proof:

From (111) in Appendix B, we have

F^k=ρα(k=1KpS,kβSR,k+1)Ek.\displaystyle{\hat{F}_{k}}=\frac{{{\rho}}}{{{\alpha}}}\left({\sum\limits_{k=1}^{K}{{p_{{\rm S},k}}{\beta_{{\rm SR},k}}}+1}\right){E_{k}}. (46)

Substituting (42) into (46), we obtain (2). Similarly, we derive G^k=γ2αρμ𝔼{|𝐠RD1,kT𝐠RD1,k|}.{\hat{G}_{k}}={\gamma^{2}}\alpha\rho\mu{\mathbb{E}}\left\{{\left|{{\boldsymbol{\rm g}}_{\rm{RD1},k}^{T}{\boldsymbol{\rm g}}_{\rm{RD1},k}^{*}}\right|}\right\}. Using the fact that 𝔼{|𝐠RD1,kT𝐠RD1,k|}=M1βSR,k{\mathbb{E}}\left\{{\left|{{\boldsymbol{\rm g}}_{\rm{RD1},k}^{T}{\boldsymbol{\rm g}}_{\rm{RD1},k}^{*}}\right|}\right\}={M_{1}}{\beta_{{\rm SR},k}}, we can then derive G^k{\hat{G}_{k}}. ∎

Following a similar reasoning as in the exact achievable rate analysis, we substitute (27) and (28) into Theorem 2 to deduce the compact expression for the approximate achievable rate as

R^k=((τc2τp)/2τc)k=1Klog2(1+ν^k),{\hat{R}_{k}}=\left({\left({{\tau_{c}}-2{\tau_{p}}}\right)}/{{2{\tau_{c}}}}\right)\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+{\hat{\nu}_{k}}}\right)}, (47)

where

ν^k=\displaystyle{\hat{\nu}_{k}}= pS,k/ξ^k,\displaystyle{{{p_{{\rm S},k}}}}/{{{\hat{\xi}_{k}}}}, (48)
ξ^k=\displaystyle{\hat{\xi}_{k}}= i=1KpS,ia^ki+(i=1KpS,ibki+ck)/pR+d^k,\displaystyle\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\hat{a}_{ki}}}+\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{b_{ki}}}+{c_{k}}}\right)/p_{\rm R}+{\hat{d}_{k}}, (49)
d^k=\displaystyle{\hat{d}_{k}}= 1(M0+αM1)βSR,k+1(M0+αM1)2m=1KβSR,mβRD,mβSR,k2βRD,k,\displaystyle\frac{1}{{\left({{M_{0}}+\alpha{M_{1}}}\right)}{{\beta_{{\rm SR},k}}}}+\frac{1}{{{{\left({{M_{0}}+\alpha{M_{1}}}\right)}^{2}}}}\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{\beta_{{\rm SR},k}^{2}{\beta_{{\rm RD},k}}}},

with a^ki\hat{a}_{ki} given as

a^ki=1(M0+αM1)βSR,iβSR,k{+βSR,iβRD,iβSR,kβRD,k\displaystyle{\hat{a}_{ki}}=\frac{1}{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}}\frac{{{\beta_{{\rm SR},i}}}}{{{\beta_{{\rm SR},k}}}}\Big\{1\!+\!\frac{{{\beta_{{\rm SR},i}}{\beta_{{\rm RD},i}}}}{{{\beta_{{\rm SR},k}}{\beta_{{\rm RD},k}}}}
×[+α2ρ2M12(M0+αM1)3]+1(M0+αM1)m=1KβSR,mβRD,mβSR,kβRD,k}.\displaystyle\times\left[{1\!+\!\frac{{{\alpha^{2}}{\rho^{2}}M_{1}^{2}}}{{{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}^{3}}}}}\right]+\frac{1}{{\left({{M_{0}}\!+\!\alpha{M_{1}}}\right)}}\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{{\beta_{{\rm SR},k}}{\beta_{{\rm RD},k}}}}\Big\}.

Note that bkib_{ki} and ckc_{k} have been defined in (42) and (43), respectively.

It is clear to see from (47) that the approximate achievable rate R^k{\hat{R}_{k}} increases with the total power of the relay pRp_{R}, and decreases with the transmit power of other sources. The mixed-ADC/DAC multipair mMIMO relaying system is interference-limited, which is consistent with [25]. Moreover, we find that including more low-resolution ADCs and DACs decreases R^k{\hat{R}_{k}}. This is reasonable since the quantization noise increases.

III-C Power Scaling Law

In this subsection, we investigate the potential for power saving in the data transmission phase due to the deployment of a very large antenna array at the relay. Here, let pS=ES/Mp_{\rm S}={{{E_{\rm S}}}\mathord{\left/{\vphantom{{{E_{\rm S}}}M}}\right.\kern-1.2pt}M} (i.e., the power of all sources is the same, pS=pS,kp_{\rm S}=p_{{\rm S},k}, k=1,,Kk=1,\ldots,K) and pR=ER/Mp_{\rm R}={{{E_{\rm R}}}\mathord{\left/{\vphantom{{{E_{\rm R}}}M}}\right.\kern-1.2pt}M}, where the transmit power of the source ESE_{S} and of the relay ERE_{R} are fixed. As MM\to\infty, the exact and approximate achievable rate for the considered system is provided in the following Corollary.

Corollary 1.

With pS=ES/Mp_{\rm S}={{{E_{\rm S}}}\mathord{\left/{\vphantom{{{E_{\rm S}}}M}}\right.\kern-1.2pt}M}, pR=ER/Mp_{\rm R}={{{E_{\rm R}}}\mathord{\left/{\vphantom{{{E_{\rm R}}}M}}\right.\kern-1.2pt}M} and ESE_{\rm S}, ERE_{\rm R} fixed, the achievable rate limit of mixed-ADC/DAC multipair mMIMO relaying systems is given by (50) at the bottom of this page.

Proof:

We start with the approximate achievable rate R^k{\hat{R}_{k}}. Let pS=ES/Mp_{\rm S}={{{E_{\rm S}}}\mathord{\left/{\vphantom{{{E_{\rm S}}}M}}\right.\kern-1.2pt}M}, pR=ER/Mp_{\rm R}={{{E_{\rm R}}}\mathord{\left/{\vphantom{{{E_{\rm R}}}M}}\right.\kern-1.2pt}M}, and with the help of (48) and (49), we can obtain

ν^k=ESESi=1Ka^ki+ER1(i=1KESMbki+M2ck)+Md^k.{\hat{\nu}_{k}}=\frac{{{E_{{\rm S}}}}}{{{E_{{\rm S}}}\sum\limits_{i=1}^{K}{{\hat{a}_{ki}}}+E_{\rm R}^{-1}\left({\sum\limits_{i=1}^{K}{{E_{{\rm S}}}{M}{b_{ki}}}+{M^{2}}{c_{k}}}\right)+M{\hat{d}_{k}}}}. (51)

As MM\to\infty, the terms related to MM in (51) are derived as

a^ki0,\displaystyle\hat{a}_{ki}\to 0, (52)
MbkiβSR,i2βRD,i/(α+ρκ)βSR,k2βRD,k2,\displaystyle M{b_{ki}}\to{{\beta_{{\rm SR},i}^{2}{\beta_{{\rm RD},i}}}}/{{\left({\alpha+\rho\kappa}\right)}}{{\beta_{{\rm SR},k}^{2}\beta_{{\rm RD},k}^{2}}}, (53)
M2ck1(α+ρκ)2m=1KβSR,mβRD,mβSR,k2βRD,k2,\displaystyle M^{2}{c_{k}}\to\frac{1}{{{{\left({{\alpha+\rho\kappa}}\right)}^{2}}}}\frac{{\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}{{\beta_{{\rm SR},k}^{2}\beta_{{\rm RD},k}^{2}}}, (54)
Md^k1(α+ρκ)1βSR,k.\displaystyle M{\hat{d}_{k}}\to\frac{1}{{\left({\alpha+\rho\kappa}\right)}}\frac{1}{{{\beta_{{\rm SR},k}}}}. (55)

Substituting (52), (53), (54) and (55) into (51), we can derive the limit of R^k\hat{R}_{k} after some simple mathematical manipulations. Following a similar way, the limit of exact achievable rate RkR_{k} can be derived. With pS=ES/Mp_{\rm S}={{{E_{\rm S}}}\mathord{\left/{\vphantom{{{E_{\rm S}}}M}}\right.\kern-1.2pt}M} and pR=ER/Mp_{\rm R}={{{E_{\rm R}}}\mathord{\left/{\vphantom{{{E_{\rm R}}}M}}\right.\kern-1.2pt}M}, (40) can be rewritten as

νk=ESESi=1Kaki+ER1(i=1KESMbki+M2ck)+Mdk.{\nu_{k}}=\frac{{{E_{{\rm S}}}}}{{{E_{{\rm S}}}\sum\limits_{i=1}^{K}{{a_{ki}}}\!+\!E_{\rm R}^{\!-\!1}\left({\sum\limits_{i=1}^{K}{{E_{{\rm S}}}{M}{b_{ki}}}\!+\!{M^{2}}{c_{k}}}\right)\!+\!M{d_{k}}}}. (56)

As MM\to\infty, the terms related to MM in (56) are given by

aki0,\displaystyle a_{ki}\to 0, (57)
Mdk1(α+ρκ)1βSR,k,\displaystyle M{d_{k}}\to\frac{1}{{\left({\alpha+\rho\kappa}\right)}}\frac{1}{{{\beta_{{\rm SR},k}}}}, (58)

and the limits of MbkiM{b_{ki}} and M2ckM^{2}{c_{k}} are given by (53) and (54), respectively. Since the limit of akia_{ki} and dkd_{k} are separately the same as a^ki\hat{a}_{ki} and d^k\hat{d}_{k} with MM\to\infty, RkR_{k} approaches the same constant limit as R^k\hat{R}_{k}. After some simplifications, the proof is concluded by deriving (50). ∎

It is clear from (50) that both exact and approximate results tend to a same constant value with MM\to\infty. We can find that the proportion of the high-resolution ADCs/DACs κ\kappa and the distortion factor of the low-resolution ADCs/DACs ρ\rho have effects on the limit rate when scaling down the transmit power proportion to 1/M{{{1}}\mathord{\left/{\vphantom{{{1}}M}}\right.\kern-1.2pt}M}. More specifically, the limit can be improved by increasing κ\kappa. Adopting the fact that ρκ+α=(1κ)α+κ\rho\kappa+\alpha=\left(1-\kappa\right)\alpha+\kappa is a monotonic increasing function of α\alpha, the limit of achievable rate monotonically increases with α\alpha, which means that we can boost the achievable rate by using higher quantization bits in the M1M_{1} low-resolution ADCs/DACs.

Proposition 1.

With pS=ES/Mp_{\rm S}={{{E_{\rm S}}}\mathord{\left/{\vphantom{{{E_{\rm S}}}M}}\right.\kern-1.2pt}M}, pR=ER/Mp_{\rm R}={{{E_{\rm R}}}\mathord{\left/{\vphantom{{{E_{\rm R}}}M}}\right.\kern-1.2pt}M} and ES0E_{\rm S}\to 0, ERE_{\rm R} fixed, we can derive the factor of the sum rate gap between the mixed-ADC/DAC relay system and the unquantized one as

RkRkp(α+ρκ)2(ERβSR,kβRD,k2+m=1KβSR,mβRD,m)(α+ρκ)ERβSR,kβRD,k2+m=1KβSR,mβRD,m.\frac{R_{k}}{R^{p}_{k}}\to\frac{{{{\left({\alpha+\rho\kappa}\right)}^{2}}\left({{E_{\rm{R}}}{\beta_{{\rm SR},k}}\beta_{{\rm RD},k}^{2}+\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}\right)}}{{\left({\alpha+\rho\kappa}\right){E_{\rm{R}}}{\beta_{{\rm SR},k}}\beta_{{\rm RD},k}^{2}+\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}.
Proposition 2.

With pS=ES/Mp_{\rm S}={{{E_{\rm S}}}\mathord{\left/{\vphantom{{{E_{\rm S}}}M}}\right.\kern-1.2pt}M}, pR=ER/Mp_{\rm R}={{{E_{\rm R}}}\mathord{\left/{\vphantom{{{E_{\rm R}}}M}}\right.\kern-1.2pt}M} and ESE_{\rm S} fixed, ER0E_{\rm R}\to 0, we can derive the factor of the sum rate gap between the mixed-ADC/DAC relay system and the unquantized one as

RkRkp(α+ρκ)2(ESm=1KβSR,m2βRD,m+m=1KβSR,mβRD,m)(α+ρκ)ESm=1KβSR,m2βRD,m+m=1KβSR,mβRD,m.\frac{R_{k}}{R^{p}_{k}}\to\frac{{{{\left({\alpha+\rho\kappa}\right)}^{2}}\left({{E_{S}}\sum\limits_{m=1}^{K}{\beta_{{\rm SR},m}^{2}{\beta_{{\rm RD},m}}}+\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}\right)}}{{\left({\alpha+\rho\kappa}\right){E_{S}}\sum\limits_{m=1}^{K}{\beta_{{\rm SR},m}^{2}{\beta_{{\rm RD},m}}}+\sum\limits_{m=1}^{K}{{\beta_{{\rm SR},m}}{\beta_{{\rm RD},m}}}}}. (59)

From Propositions 1 and 2, it is clear that the sum rate gap between the mixed-ADC/DAC relay system and the unquantized system is a constant factor in the low power regime. The factors in the case ES0E_{\rm S}\to 0 or ER0E_{\rm R}\to 0 are both related to α+ρκ\alpha+\rho\kappa. Using the fact that ρκ+α=(1κ)α+κ\rho\kappa+\alpha=\left(1-\kappa\right)\alpha+\kappa is a monotonic increasing function of α\alpha, we can find that the factors also increase with α\alpha. For the special case of α=1\alpha=1, i.e., α+ρκ=1\alpha+\rho\kappa=1, the achievable rate of the mixed-ADC/DAC system is the same as that of the ideal unquantized system.

IV Power Allocation

In this section, we try to maximize the sum achievable rate of the mixed-ADC/DAC multipair mMIMO relaying system constrained to a given total sum power PTP_{\rm T}, i.e., k=1KpS,k+pRPT\sum\limits_{k=1}^{K}{{p_{{\rm S},k}}}+p_{\rm R}\leq P_{\rm T}, and formulate it as a power allocation problem.

Let us define 𝐩S=[pS,1,,pS,K]T{\boldsymbol{\rm p}_{\rm S}}={\left[{{{p_{{S},1}},\ldots,{p_{{S},K}}}}\right]}^{T}, the power allocation problem can be expressed as

𝒫1:maximize𝐩S,pR\displaystyle\mathcal{P}_{1}:\mathop{{\rm{maximize}}}\limits_{{\boldsymbol{\rm p}_{\rm S}},p_{\rm R}}\quad τc2τp2τck=1Klog2(1+νk)\displaystyle\frac{{{\tau_{c}}-2{\tau_{p}}}}{{2{\tau_{c}}}}\sum\limits_{k=1}^{K}{{{\log}_{2}}\left({1+{\nu_{k}}}\right)} (60)
subjectto\displaystyle\rm{subject}\ \rm{to}\quad k=1KpS,k+pRPT\displaystyle\sum\limits_{k=1}^{K}{{p_{{\rm S},k}}}+p_{\rm R}\leq P_{\rm T} (61)
𝐩S𝟎,pR0.\displaystyle{\boldsymbol{\rm p}_{\rm S}}\geq\boldsymbol{\rm 0},p_{\rm R}\geq 0. (62)

Since log()\log\left({\cdot}\right) is a monotonic increasing function, the problem 𝒫1\mathcal{P}_{1} can be reformulated as

𝒫2:minimize𝐩S,pR\displaystyle\mathcal{P}_{2}:\mathop{{\rm{minimize}}}\limits_{{\boldsymbol{\rm p}_{\rm S}},p_{\rm R}}\quad k=1K(1+νk)1\displaystyle\prod\limits_{k=1}^{K}{\left({1+{\nu_{k}}}\right)}^{-1} (63)
subjectto\displaystyle\rm{subject}\ \rm{to}\quad k=1KpS,k+pRPT\displaystyle\sum\limits_{k=1}^{K}{{p_{{\rm S},k}}}+p_{\rm R}\leq P_{\rm T} (64)
𝐩S𝟎,pR0.\displaystyle{\boldsymbol{\rm p}_{\rm S}}\geq\boldsymbol{\rm 0},p_{\rm R}\geq 0. (65)

Problem 𝒫2\mathcal{P}_{2} is a general nonconvex complementary geometric program (CGP), which can be approximated by solving a sequence of GP problems. After that, we can use standard convex optimization tools (e.g., CVX) to solve the GP problems [35]. The detailed steps of the power allocation algorithm are provided in Algorithm 1. Following the successive approximation algorithm in [25], it is efficient to solve the power allocation problem.

1) Initialization. Define a tolerance ϵ\epsilon and parameter θ\theta. Set j=1j=1 and set the initial value ν~k\tilde{\nu}_{k} according to the SINR expression in Theorem 1 with pS,k=PT2Kp_{S,k}=\frac{P_{T}}{2K} and pR=PT2p_{R}=\frac{P_{T}}{2}. 2) iteration jj. Compute δk=ν~k1+ν~k\delta_{k}=\frac{\tilde{\nu}_{k}}{1+\tilde{\nu}_{k}}. Then solve the GP problem 𝒫3\mathcal{P}_{3}:
𝒫3:minimize𝐩S,pR\displaystyle\mathcal{P}_{3}:\mathop{{\rm{minimize}}}\limits_{{\boldsymbol{\rm p}_{\rm S}},p_{\rm R}}\quad k=1Kνkδk\displaystyle\prod\limits_{k=1}^{K}{\nu_{k}}^{-\delta_{k}} (66)
subjectto\displaystyle\rm{subject}\ \rm{to}\quad θ1ν~kνkθν~k,k=1,,K\displaystyle{\theta^{-1}\tilde{\nu}_{k}\leq\nu_{k}\leq\theta\tilde{\nu}_{k},\ k=1,\cdots,K} (67)
νkpS,k1ξk1,k=1,,K\displaystyle{\nu_{k}{p_{S,k}}^{-1}\xi_{k}\leq 1,\ k=1,\cdots,K} (68)
k=1KpS,k+pRPT\displaystyle{\sum\limits_{k=1}^{K}{{p_{{\rm S},k}}}+p_{\rm R}\leq P_{\rm T}} (69)
𝐩S𝟎,pR0.\displaystyle{\boldsymbol{\rm p}_{\rm S}}\geq\boldsymbol{\rm 0},p_{\rm R}\geq 0. (70)
Denote the optimal solutions by νk(j){\nu_{k}}^{\left(j\right)}, for k=1,,Kk=1,\cdots,K. 3) Stopping criterion. If maxk|νk(j)ν~k|<ϵ{\rm{max}}_{k}\left|{{\nu_{k}}^{\left(j\right)}-{\tilde{\nu}_{k}}}\right|<\epsilon, stop; otherwise, go to step 4). 4) Update initial values. Set νk=νk(j)\nu_{k}={\nu_{k}}^{\left(j\right)}, and j=j+1j=j+1. Go to step 2)
Algorithm 1 Successive approximation algorithm for 𝒫2\mathcal{P}_{2}
Refer to caption
Fig. 1: Sum achievable rate against the number of relay antennas MM for pS=10p_{\rm S}=10 dB, pR=10p_{\rm R}=10 dB, K=10K=10 and κ=1/2\kappa=1/2.

V Numerical Results

In this section, we conduct numerous simulations to verify the accuracy of the analytical results. Apart from that, insights are also provided. Then, the benefit of the proposed power allocation algorithm is demonstrated. Moreover, we investigate the energy efficiency to show the advantage of the mixed-ADC/DAC architecture.

In the Monte Carlo simulation, we assume that the users are distributed in a hexagonal cell with a radius of 1000 meters, while the minimum distance between the users and relay is rmin=100r_{\min}=100 meters. The length of the coherence interval and pilot sequence are set as τc=20K\tau_{c}=20K and τp=K\tau_{p}=K, respectively. Furthermore, the large-scale fading coefficients are arbitrarily generated by βSR,k=zk(rSR,k/rmin)α{\beta_{{\rm SR},k}}={z_{k}}{\left({{{{r_{{\rm SR},k}}}\mathord{\left/{\vphantom{{{r_{{\rm SR},k}}}{{r_{\min}}}}}\right.\kern-1.2pt}{{r_{\min}}}}}\right)^{-\alpha}} and βRD,k=zk(rRD,k/rmin)α{\beta_{{\rm RD},k}}={z_{k}}{\left({{{{r_{{\rm RD},k}}}\mathord{\left/{\vphantom{{{r_{{\rm RD},k}}}{{r_{\min}}}}}\right.\kern-1.2pt}{{r_{\min}}}}}\right)^{-\alpha}}, where zkz_{k} is a log-normal random variable with standard derivation 8 dB, rSR,k{r_{{\rm SR},k}} and rRD,k{r_{{\rm RD},k}} represent the distances from the sources to the relay and destinations to the relay, respectively, and α=3.8\alpha=3.8 denotes the pathloss exponent.

V-A Achievable Rate

In Fig. 1, the simulated achievable rate, the analytical exact result (39) and the approximate result (47) are plotted against the number of relay antennas. It can be seen that the analytical exact and approximate results, as well as simulation results are close to each other, which validates the correctness of our derived expressions. For a small number of antennas at the relay, the sum achievable rates for the cases of b=1,2,b=1,2,\infty matches well with each other. While as the quantization bits bb increase, the gap between the approximate and simulated curves becomes small. Finally, better rate performance is achieved with a larger number of quantization bits (b>1b>1).

In Fig. 2, we investigate the power scaling law of the mixed-ADC/DAC multipair mMIMO relaying system. The fraction of the number of high-resolution ADCs/DACs in the mixed-ADC/DAC architecture is κ=0,1/2\kappa=0,1/2 and 1, respectively. It can be seen that the exact and approximate expressions tend to a constant for MM\to\infty, and a higher value of κ\kappa increases the achievable rate, which agrees with Corollary 1.

Refer to caption
Fig. 2: Sum achievable rate against the number of relay antennas for pS=10p_{\rm S}=10 dB, pR=10p_{\rm R}=10 dB, K=10K=10 and b=1b=1.
Refer to caption
Fig. 3: Sum achievable rate against the number of relay antennas MM for K=10K=10, κ=1/2\kappa=1/2 and PT=10P_{\rm T}=10 dB.

V-B Power Allocation

We show the impact of the efficient power allocation scheme on the sum achievable rate in Fig. 3. The uniform power allocation scheme, i.e., pS=PT2Kp_{\rm S}=\frac{P_{\rm T}}{2K} and pR=PTKp_{\rm R}=\frac{P_{\rm T}}{K}, is also investigated as a benchmark for comparison. It is clear that the proposed optimal power allocation scheme significantly boosts the sum rate compared with the one of the unquantized system with uniform power allocation. This important finding demonstrates the significance of adopting an efficient power allocation scheme in the mixed-ADC/DAC multipair mMIMO relaying system.

V-C Energy Efficiency

Up to now, we have investigated the achievable rate of mixed-ADC/DAC multipair mMIMO relaying systems. As expected, under the same power allocation, the sum rate of unquantized system outperforms the one with mixed-resolution ADCs/DACs, at the cost of expensive hardware and power consumption. There should be a fundamental trade-off between the achievable rate and energy efficiency, and therefore, we also study the energy efficiency of mixed-ADC/DAC multipair mMIMO relaying systems.

According to [21], the energy efficiency can be defined as

ηEE=R×BPtotalbit/Joule,{\eta_{EE}}=\frac{{R\times B}}{{{P_{\rm total}}}}{{\rm bit}\mathord{\left/{\vphantom{{\rm bit}{\rm Joule}}}\right.\kern-1.2pt}{\rm Joule}}, (71)

where RR denotes the sum achievable rate, BB refers to the transmission bandwidth assumed to be 20 MHz, and PtotalP_{\rm total} is the total power consumption. Combining [21, Eq. (43)] and [36, Eq. (9)], PtotalP_{\rm total} can be expressed as

Ptotal=M(Pmix+Pfilt)+2Psyn+M1(cPAGC+PDACL)\displaystyle{P_{{\rm{total}}}}=M\left({{P_{\rm mix}}+{P_{\rm filt}}}\right)+2{P_{\rm syn}}+{M_{1}}\left({c{P_{\rm AGC}}+P_{\rm DAC}^{L}}\right)
+M(PLNA+Pmix+PIFA+Pfilr)+M0(PAGC+PADCH)\displaystyle+M\left({{P_{\rm LNA}}+{P_{\rm mix}}+{P_{\rm IFA}}+{P_{\rm filr}}}\right)+{M_{0}}\left({{P_{\rm AGC}}+P_{\rm ADC}^{H}}\right)
+M1(cPAGC+PADCL)+M0(PAGC+PDACH),\displaystyle+{M_{1}}\left({c{P_{\rm AGC}}+P_{\rm ADC}^{L}}\right)+{M_{0}}\left({{P_{\rm AGC}}+P_{\rm DAC}^{H}}\right), (72)

where PmixP_{\rm mix}, PfiltP_{\rm filt}, PsynP_{\rm syn}, PLNAP_{\rm LNA}, PIFAP_{\rm IFA}, PfilrP_{\rm filr}, PAGCP_{\rm AGC}, PADCHP_{\rm ADC}^{H}, PADCLP_{\rm ADC}^{L}, PDACHP_{\rm DAC}^{H} and PDACLP_{\rm DAC}^{L} are the power consumption values for the mixer, the active filters at the transmitter side, the frequency synthesizer, low-noise amplifiers (LNA), the intermediate frequency amplifier (IFA), the active filters at the receiver side, the automatic gain control (AGC), high-resolution ADCs, low-resolution ADCs, high-resolution DACs and low-resolution DACs, respectively. In addition, cc denotes the flag related to quantization bits of low-resolution ADCs, which is given by

c={0,b=1,1,b>1.c=\left\{{\begin{array}[]{ll}0,\quad b=1,\\ 1,\quad b>1.\end{array}}\right. (73)

According to [37], the power consumed in DACs and ADCs can be respectively expressed in terms of the number of quantization bits as

PDAC\displaystyle{P_{\rm DAC}} =12VddI0(2b1)+bCp(2B+fcor)Vdd2,\displaystyle=\frac{1}{2}{V_{dd}}{I_{0}}\left({{2^{b}}-1}\right)+b{C_{p}}\left({2B+{f_{cor}}}\right)V_{dd}^{2}, (74)
PADC\displaystyle{P_{\rm ADC}} =3Vdd2Lmin(2B+fcor)100.1525b+4.838,\displaystyle=\frac{{3V_{dd}^{2}{L_{\min}}\left({2B+{f_{cor}}}\right)}}{{{{10}^{-0.1525b+4.838}}}}, (75)

where bb denotes the quantization bits, BB is the bandwidth of the original signal assumed to be 20 MHz, VddV_{dd} is the power supply of converter, I0I_{0} is the unit current source corresponding to the least significant bit (LSB), CpC_{p} is the parasitic capacitance of each switch in the converter, LminL_{min} is the minimum channel length for the given CMOS technology, fcorf_{cor} is the corner frequency of the 1/f1/f noise and all those parameters are specifically defined in [37]. (74) holds for binary-weighted current-steering DACs [37] and (75) is established for the complete class of CMOS Nyquist-rate high speed ADCs [38]. In numerical examples, we consider the following classical values: Pmix=30.3P_{\rm mix}=30.3 mW, Pfilt=Pfilr=2.5P_{\rm filt}=P_{\rm filr}=2.5 mW, Psyn=50.0P_{\rm syn}=50.0 mW, PLNA=20P_{\rm LNA}=20 mW, PIFA=3P_{\rm IFA}=3 mW and PAGC=2P_{\rm AGC}=2 mW as in [21] and [36], and the power consumption values of other various circuit blocks have been discussed in [37].

Refer to caption
Fig. 4: Energy efficiency of multipair mMIMO relaying systems with mixed-resolution ADCs/DACs against the number of low-resolution quantization bits bb for pS=10p_{\rm S}=10 dB, pR=10p_{\rm R}=10 dB and K=10K=10.

The energy efficiency of mixed-ADC/DAC multipair mMIMO relaying systems against the quantization bits is illustrated in Fig. 4. It is clear from the figure that the relay adopting pure low-resolution ADCs/DACs attains the best energy efficiency. That means the energy efficiency increases with the proportion of the number of low-resolution ADCs/DACs (1κ1-\kappa) in the mixed-ADC/DAC architecture. Although the pure low-resolution ADC/DAC architecture can achieve better energy efficiency than the mixed-ADC/DAC architecture, as shown in 1, the spectral efficiency of the low-resolution ADC/DAC architecture is much lower than that of the mixed-ADC/DAC architecture. Moreover, the channel estimation in the mixed-ADC/DAC architecture is more tractable than that in the low-resolution ADC/DAC architecture due to the use of partial high-resolution ADCs. Moreover, the pure low-resolution ADCs/DACs has low spectrum efficiency. Fig. 4 indicates that we can achieve a better spectral efficiency by reducing the burden of power consumption considerably by using the mixed-ADC/DAC architecture.

VI Conclusions

In this paper, we investigated the achievable rate of a multipair mMIMO relaying system with mixed-resolution ADCs and DACs at the relay. Both exact and approximate closed-form expressions for the achievable rate were derived. Then, we proved that the transmit power of each users can be scaled down as 1/M{1/M} for the considered system. Despite the rate loss due to the use of low-resolution ADCs and DACs, employing massive antenna arrays still enables high achievable rate and large power saving. Furthermore, we proposed an efficient power allocation scheme, which can compensate for the rate degradation caused by low-resolution ADCs and DACs. Finally, the energy efficiency was investigated, and showed that the mixed-ADC/DAC architecture can attain a considerable rate and energy efficiency simultaneously, which is promising for practical mMIMO relaying systems.

Appendix A Proof of Lemma 1

For 𝔼{𝐱~R2}=𝔼{𝐱~R02}+𝔼{𝐱~R12}{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R}}}}\right\|}^{2}}}\right\}}={{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R0}}}}\right\|}^{2}}}\right\}}+{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}}, we first

A-1 Calculate 𝔼{𝐱~R02}{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R0}}}}\right\|}^{2}}}\right\}}

𝐑𝐱~R0𝐱~R0\displaystyle{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}}} =𝐆RD0𝐆SR0H𝐑𝐲~R0𝐲~R0𝐆SR0𝐆RD0T\displaystyle={\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T} (76)
+𝐆RD0𝐆SR0H𝐑𝐲~R0𝐲~R1𝐆SR1𝐆RD0T\displaystyle+{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}
+𝐆RD0𝐆SR1H𝐑𝐲~R1𝐲~R0𝐆SR0𝐆RD0T\displaystyle+{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}
+𝐆RD0𝐆SR1H𝐑𝐲~R1𝐲~R1𝐆SR1𝐆RD0T\displaystyle+{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}
=𝐐1+𝐐2+𝐐3+𝐐4,\displaystyle={\boldsymbol{\rm Q}}_{1}+{\boldsymbol{\rm Q}}_{2}+{\boldsymbol{\rm Q}}_{3}+{\boldsymbol{\rm Q}}_{4},

where 𝐑𝐲~R0𝐲~R0{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}}, 𝐑𝐲~R0𝐲~R1{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}}, 𝐑𝐲~R1𝐲~R0{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R0}}}}}} and 𝐑𝐲~R1𝐲~R1{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{\tilde{y}}}}}_{{\rm{R1}}}}}} are separately given by (21) to (24), and 𝐐1{\boldsymbol{\rm Q}}_{1}, 𝐐2{\boldsymbol{\rm Q}}_{2} ,𝐐3{\boldsymbol{\rm Q}}_{3} and 𝐐4{\boldsymbol{\rm Q}}_{4} are defined as the four parts of 𝐑𝐱~R0𝐱~R0{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}}}, respectively.

𝔼{𝐐1}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm Q}}_{1}}\right\} =𝔼{𝐆RD0𝐆SR0H𝐆SR0𝐏S𝐆SR0H𝐆SR0𝐆RD0T}\displaystyle=\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm G}_{\rm SR0}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR0}^{H}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}
+𝔼{𝐆RD0𝐆SR0H𝐆SR0𝐆RD0T}.\displaystyle+\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}. (77)

Using 𝔼{𝐠SR0,m4}=M0(M0+1)βSR,m2\mathbb{E}\left\{{{{{\left\|{{\boldsymbol{\rm g}_{{\rm{SR0}},m}}}\right\|}^{4}}}}\right\}=M_{0}\left({M_{0}+1}\right)\beta_{{\rm SR},m}^{2}, we have

𝔼{𝐆RD0𝐆SR0H𝐆SR0𝐆RD0T}=M0m=1KβSR,mβRD,m𝐈M0,\displaystyle{\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}}=M_{0}\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}{{\boldsymbol{\rm I}}_{M0}}, (78)
𝔼{𝐆RD0𝐆SR0H𝐆SR0𝐏S𝐆SR0H𝐆SR0𝐆RD0T}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm G}_{\rm SR0}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR0}^{H}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR0}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}
=M0m=1KβSR,mβRD,m(i=1KpS,iβSR,i+M0pS,mβSR,m)𝐈M0.\displaystyle=M_{0}\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm{SR}},i}}}\!+\!M_{0}{p_{{\rm S},m}}{\beta_{{\rm{SR}},m}}}\right)}{{\boldsymbol{\rm I}}_{M0}}. (79)

Then substituting (A-1) and (78) into (A-1), we directly obtain

𝔼{𝐐1}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm Q}}_{1}}\right\} =M0m=1KβSR,mβRD,m\displaystyle=M_{0}\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}} (80)
×(i=1KpS,iβSR,i+M0pS,mβSR,m+1)𝐈M0.\displaystyle\times\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm{SR}},i}}}+M_{0}{p_{{\rm S},m}}{\beta_{{\rm{SR}},m}}+1}\right){{\boldsymbol{\rm I}}_{M0}}.

Similar to the computation of 𝐐1{\boldsymbol{\rm Q}}_{1}, the expectation of 𝐐2{\boldsymbol{\rm Q}}_{2}, 𝐐3{\boldsymbol{\rm Q}}_{3} and 𝐐4{\boldsymbol{\rm Q}}_{4} can be derived respectively as

𝔼{𝐐2}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm Q}}_{2}}\right\} =αM0M1m=1KpS,mβSR,m2βRD,m𝐈M0,\displaystyle=\alpha M_{0}M_{1}\sum\limits_{m=1}^{K}{{p_{{\rm S},m}}{\beta^{2}_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}{{\boldsymbol{\rm I}}_{M0}}, (81)
𝔼{𝐐3}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm Q}}_{3}}\right\} =αM0M1m=1KpS,mβSR,m2βRD,m𝐈M0.\displaystyle=\alpha M_{0}M_{1}\sum\limits_{m=1}^{K}{{p_{{\rm S},m}}{\beta^{2}_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}{{\boldsymbol{\rm I}}_{M0}}. (82)
𝔼{𝐐4}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm Q}}_{4}}\right\} =α2M1m=1KβSR,mβRD,m\displaystyle=\alpha^{2}M_{1}\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}
×(i=1KpS,iβSR,i+M1pS,mβSR,m+1)𝐈M0\displaystyle\times\left(\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm{SR}},i}}}+M_{1}{p_{{\rm S},m}}{\beta_{{\rm{SR}},m}}+1\right){{\boldsymbol{\rm I}}_{M0}}
+αρ𝔼{𝐆RD0𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1𝐆RD0T}.\displaystyle+\alpha\rho\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}. (83)

where

𝔼{𝐆RD0𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1𝐆RD0T}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}
=𝔼{𝐆RD0𝔼[𝐆SR1Hdiag(𝐆SR1𝐏S𝐆SR1H)𝐆SR1]𝐆RD0T}\displaystyle=\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}\mathbb{E}\left[{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}{\rm{diag}}\left({\boldsymbol{\rm G}_{\rm SR1}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}}\right){{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}\right]{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}
+M1m=1KβSR,mβRD,m𝐈M0.\displaystyle+M_{1}\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}{{\boldsymbol{\rm I}}_{M0}}. (84)

The expectation of the diagonal term can be decomposed as

𝔼{𝐠SR1,iHdiag(𝐆SR1𝐏S𝐆SR1H)𝐠SR1,j}\displaystyle\mathbb{E}\left\{{\boldsymbol{\rm{g}}}_{{\rm{SR1}},i}^{H}{\rm{diag}}\left({\boldsymbol{\rm G}_{\rm SR1}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}}\right){{\boldsymbol{\rm{g}}}_{{\rm{SR1}},j}}\right\} (85)
=\displaystyle= {0,ij,m=1M1pS,n𝔼{|gSR1,mn|4}+m=1M1inKpS,i𝔼{|gSR1,mn|2}𝔼{|gSR1,mi|2},i=j=n.\displaystyle\left\{{\begin{array}[]{ll}0,&i\neq j,\\ {\begin{aligned} &\sum\limits_{m=1}^{{M_{1}}}{{p_{{S},n}}{\mathbb{E}}\left\{{{{\left|{{g_{{\rm{SR1}},mn}}}\right|}^{4}}}\right\}}\\ &\!+\!\sum\limits_{m\!=\!1}^{{M_{1}}}{\sum\limits_{i\neq n}^{K}{{p_{{S},i}}{\mathbb{E}}\left\{{{{\left|{{g_{{\rm{SR1}},{{mn}}}}}\right|}^{2}}}\right\}{\mathbb{E}}\left\{{{{\left|{{g_{{\rm{SR1}},{{mi}}}}}\right|}^{2}}}\right\}}}\end{aligned}},&i\!=\!j\!=\!n.\end{array}}\right.

Applying the fact that 𝔼{|gSR1,mn|4}=2βSR,n2{\mathbb{E}}\left\{{{{\left|{{g_{{\rm{SR1}},{{mn}}}}}\right|}^{4}}}\right\}=2\beta^{2}_{{\rm{SR}},n}, (85) can be expressed as

𝔼{𝐠SR1,iHdiag(𝐆SR1𝐏S𝐆SR1H)𝐠SR1,j}\displaystyle\mathbb{E}\left\{{\boldsymbol{\rm{g}}}_{{\rm{SR1}},i}^{H}{\rm{diag}}\left({\boldsymbol{\rm G}_{\rm SR1}{{\boldsymbol{\rm{P}}}_{S}}{\boldsymbol{\rm G}_{\rm SR1}^{H}}}\right){{\boldsymbol{\rm{g}}}_{{\rm{SR1}},j}}\right\} (86)
=\displaystyle= {0,ij,M1βSR,n(i=1KpS,iβSR,i+pS,nβSR,n),i=j=n.\displaystyle\left\{{\begin{array}[]{ll}0,&i\neq j,\\ {M_{1}}{\beta_{{\rm{SR}},n}}\left({\sum\limits_{i=1}^{K}{{p_{{S},i}}{\beta_{{\rm{SR}},i}}}+{p_{{S},n}}{\beta_{{\rm{SR}},n}}}\right),&i=j=n.\end{array}}\right.

Substituting (86) into (A-1), we can obtain

𝔼{𝐆RD0𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1𝐆RD0T}\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{\boldsymbol{\rm{G}}}_{{\rm{SR1}}}}{\boldsymbol{\rm{G}}}_{{\rm{RD0}}}^{T}}\right\}
=M1m=1KβSR,mβRD,m(i=1KpS,iβSR,i+pS,mβSR,m+1)𝐈M0.\displaystyle=M_{1}\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm{SR}},i}}}\!+\!{p_{{\rm S},m}}{\beta_{{\rm{SR}},m}}\!+\!1}\right){{\boldsymbol{\rm I}}_{M0}}. (87)

With the help of (83) and (A-1), we derive 𝔼{𝐐4}\mathbb{E}\left\{{{\boldsymbol{\rm Q}}_{4}}\right\} as

𝔼{𝐐4}=αM1m=1KβSR,mβRD,m\displaystyle\mathbb{E}\left\{{{\boldsymbol{\rm Q}}_{4}}\right\}=\alpha M_{1}\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}} (88)
×(i=1KpS,iβSR,i+αM1pS,mβSR,m+ρpS,mβSR,m+1)𝐈M0.\displaystyle\times\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm{SR}},i}}}\!+\!\alpha M_{1}{p_{{\rm S},m}}{\beta_{{\rm{SR}},m}}\!+\!\rho{p_{{\rm S},m}}{\beta_{{\rm{SR}},m}}\!+\!1}\right){{\boldsymbol{\rm I}}_{M0}}.

Combining (80), (81), (82) and (88), we can derive

𝔼{𝐑𝐱~R0𝐱~R0}=𝔼{𝐑𝐱R0𝐱R0}=μ𝐈M0.\mathbb{E}\left\{{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}}}\right\}=\mathbb{E}\left\{{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{x}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{x}}}}_{{\rm{R0}}}}}}\right\}=\mu{\boldsymbol{\rm I}}_{M0}. (89)
𝔼{𝐱~R02}=𝔼{𝐱R02}=μM0.{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R0}}}}\right\|}^{2}}}\right\}}={{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{x}}}}_{\rm{R0}}}}\right\|}^{2}}}\right\}}=\mu M_{0}. (90)

A-2 Calculate 𝔼{𝐱~R12}{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}}

𝔼{𝐱~R12}=α2𝔼{𝐱R12}+𝔼{𝐧qDH𝐧qD}.\displaystyle{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}}=\alpha^{2}{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{x}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}}+\mathbb{E}\left\{{{{\boldsymbol{\rm{n}}}^{H}_{{{\rm{q}}_{D}}}}{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}}\right\}. (91)

Similar to the calculation of (89) and (90), 𝔼{𝐑𝐱R0𝐱R0}\mathbb{E}\left\{{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{x}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{x}}}}_{{\rm{R0}}}}}}\right\} and 𝔼{𝐱R12}{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{x}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}} can be respectively expressed as

𝔼{𝐑𝐱R1𝐱R1}=μ𝐈M1,\mathbb{E}\left\{{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{x}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{x}}}}_{{\rm{R1}}}}}}\right\}=\mu{\boldsymbol{\rm I}}_{M1}, (92)
𝔼{𝐱R12}=μM1.{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{x}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}}=\mu M_{1}. (93)

As for 𝔼{𝐧qDH𝐧qD}\mathbb{E}\left\{{{{\boldsymbol{\rm{n}}}^{H}_{{{\rm{q}}_{D}}}}{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}}\right\}, considering (19) and (92), we can derive 𝐑𝐧qD\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}} as

𝔼{𝐑𝐧qD}=αρ𝔼{𝐑𝐱R1𝐱R1}=αρμ𝐈M1.\mathbb{E}\left\{\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}\right\}=\alpha\rho\mathbb{E}\left\{{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{x}}}}_{{\rm{R1}}}}{{{\boldsymbol{\rm{x}}}}_{{\rm{R1}}}}}}\right\}=\alpha\rho\mu{\boldsymbol{\rm I}}_{M1}. (94)

Hence,

𝔼{𝐧qDH𝐧qD}=αρμM1.\mathbb{E}\left\{{{{\boldsymbol{\rm{n}}}^{H}_{{{\rm{q}}_{D}}}}{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{D}}}}}\right\}=\alpha\rho\mu M_{1}. (95)

Substituting (93) and (95) into (91), we can directly obtain

𝔼{𝐱~R12}=μαM1.{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}}=\mu\alpha M_{1}. (96)

Therefore, 𝔼{𝐱~R2}=𝔼{𝐱~R02}+𝔼{𝐱~R12}=μ(M0+αM1){{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R}}}}\right\|}^{2}}}\right\}}={{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R0}}}}\right\|}^{2}}}\right\}}+{{\mathbb{E}}\left\{{{{\left\|{{{{\boldsymbol{\rm{\tilde{x}}}}}_{\rm{R1}}}}\right\|}^{2}}}\right\}}=\mu\left({M_{0}+\alpha M_{1}}\right). The proof is concluded.

Appendix B Proof of Theorem 1

The argument of the log function in the right-hand side of (31) consists of six terms: 1) desired signal power AkA_{k}; 2) estimation error BkB_{k}; 3) inter-pair interference CkC_{k}; 4) noise at the relay DkD_{k} and EkE_{k}; 5) quantization noise of ADCs FkF_{k}; 6) quantization noise of DACs GkG_{k}.

B-1 Compute AkA_{k}

Since

𝔼{Tk,k}=γ(M0+αM1)2βSR,kβRD,k,\displaystyle{{\mathbb{E}}\left\{{T_{k,k}}\right\}}=\gamma\left({M_{0}+\alpha M_{1}}\right)^{2}\beta_{{\rm SR},k}\beta_{{\rm RD},k}, (97)

we have

Ak=pS,kγ2(M0+αM1)4βSR,k2βRD,k2.{A_{k}}={p_{{S},k}}{\gamma^{2}}{\left({{M_{0}}+\alpha{M_{1}}}\right)^{4}}\beta_{{\rm{SR}},k}^{2}\beta_{{\rm{RD}},k}^{2}. (98)

B-2 Compute BkB_{k}

Bk=pS,kVar(Tk,k)=pS,k𝔼{|Tk,k|2}Ak,{B_{k}}={p_{{S},k}}{\rm Var}\left({T_{k,k}}\right)={p_{{S},k}}{\mathbb{E}}\left\{{\left|T_{k,k}\right|^{2}}\right\}-A_{k}, (99)

We define t1t_{1} to t10t_{10} as the decomposed terms of 𝔼{Tk,iTk,iH}{\mathbb{E}}\left\{{T_{k,i}T^{H}_{k,i}}\right\} in order. Note that the undefined terms in tkt_{k} mean that they are included in the expressions if and only if i=ki=k.

t1\displaystyle t_{1} =𝔼{𝐠RD0,kT𝐆RD0𝐆SR0H𝐠SR0,i𝐠SR0,iH𝐆SR0𝐆RD0T𝐠RD0,k}\displaystyle={\mathbb{E}}\left\{{{\boldsymbol{\rm g}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR0}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR0}},i}}{\boldsymbol{\rm G}}_{{\rm{SR0}}}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD0}},k}}}\right\}
=M02βSR,iβRD,k(M0βSR,kβRD,k+m=1KβSR,mβRD,mCLOSE\displaystyle=M^{2}_{0}\beta_{{\rm SR},i}\beta_{{\rm RD},k}({M_{0}\beta_{{\rm SR},k}\beta_{{\rm RD},k}\!+\!\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}}
OPEN+M0βSR,iβRD,i)+M04βSR,k2βRD,k2,\displaystyle\!+\!M_{0}\beta_{{\rm SR},i}\beta_{{\rm RD},i})+{M^{4}_{0}\beta^{2}_{{\rm SR},k}\beta^{2}_{{\rm RD},k}}, (100)
t2\displaystyle t_{2} =𝔼{α2𝐠RD0,kT𝐆RD0𝐆SR1H𝐠SR1,i𝐠SR1,iH𝐆SR1𝐆RD0T𝐠RD0,k}\displaystyle={\mathbb{E}}\left\{{\alpha^{2}{\boldsymbol{\rm g}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR1}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR1}},i}}{\boldsymbol{\rm G}}_{{\rm{SR1}}}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD0}},k}}}\right\}
=α2M0M1βSR,iβRD,k(M0βSR,kβRD,k+m=1KβSR,mβRD,mCLOSE\displaystyle=\alpha^{2}M_{0}M_{1}\beta_{{\rm SR},i}\beta_{{\rm RD},k}({M_{0}\beta_{{\rm SR},k}\beta_{{\rm RD},k}\!+\!\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}}
OPEN+M1βSR,iβRD,i)+α2M02M12βSR,k2βRD,k2,\displaystyle\!+\!M_{1}\beta_{{\rm SR},i}\beta_{{\rm RD},i})+{\alpha^{2}M^{2}_{0}M^{2}_{1}\beta^{2}_{{\rm SR},k}\beta^{2}_{{\rm RD},k}}, (101)
t3\displaystyle t_{3} =𝔼{α2𝐠RD1,kT𝐆RD1𝐆SR0H𝐠SR0,i𝐠SR0,iH𝐆SR0𝐆RD1T𝐠RD1,k}\displaystyle={\mathbb{E}}\left\{{\alpha^{2}{\boldsymbol{\rm g}_{{\rm{RD1}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR0}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR0}},i}}{\boldsymbol{\rm G}}_{{\rm{SR0}}}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD1}},k}}}\right\}
=α2M0M1βSR,iβRD,k(M1βSR,kβRD,k+m=1KβSR,mβRD,mCLOSE\displaystyle=\alpha^{2}M_{0}M_{1}\beta_{{\rm SR},i}\beta_{{\rm RD},k}({M_{1}\beta_{{\rm SR},k}\beta_{{\rm RD},k}+\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}}
OPEN+M0βSR,iβRD,i)+α2M02M12βSR,k2βRD,k2,\displaystyle+M_{0}\beta_{{\rm SR},i}\beta_{{\rm RD},i})+{\alpha^{2}M^{2}_{0}M^{2}_{1}\beta^{2}_{{\rm SR},k}\beta^{2}_{{\rm RD},k}}, (102)
t4\displaystyle t_{4} =𝔼{α4𝐠RD1,kT𝐆RD1𝐆SR1H𝐠SR1,i𝐠SR1,iH𝐆SR1𝐆RD1T𝐠RD1,k}\displaystyle={\mathbb{E}}\left\{{\alpha^{4}{\boldsymbol{\rm g}_{{\rm{RD1}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR1}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR1}},i}}{\boldsymbol{\rm G}}_{{\rm{SR1}}}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD1}},k}}}\right\}
=α4M12βSR,iβRD,k(M1βSR,kβRD,k+m=1KβSR,mβRD,mCLOSE\displaystyle=\alpha^{4}M^{2}_{1}\beta_{{\rm SR},i}\beta_{{\rm RD},k}\Big({M_{1}\beta_{{\rm SR},k}\beta_{{\rm RD},k}+\sum\limits_{m=1}^{K}{{\beta_{{\rm{SR}},m}}{\beta_{{\rm{RD}},m}}}}
OPEN+M1βSR,iβRD,i)+α4M14βSR,k2βRD,k2,\displaystyle+M_{1}\beta_{{\rm SR},i}\beta_{{\rm RD},i}\Big)+{\alpha^{4}M^{4}_{1}\beta^{2}_{{\rm SR},k}\beta^{2}_{{\rm RD},k}}, (103)
t5\displaystyle t_{5} =𝔼{α𝐠RD0,kT𝐆RD0𝐆SR0H𝐠SR0,i𝐠SR1,iH𝐆SR1𝐆RD0T𝐠RD0,k}\displaystyle={\mathbb{E}}\left\{{\alpha{\boldsymbol{\rm g}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR0}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR1}},i}}{\boldsymbol{\rm G}}_{{\rm{SR1}}}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD0}},k}}}\right\}
=αM02M1βSR,i2βRD,k(βRD,i+M0βRD,k),\displaystyle=\alpha M^{2}_{0}M_{1}\beta^{2}_{{\rm SR},i}\beta_{{\rm RD},k}\left({\beta_{{\rm RD},i}+{M_{0}\beta_{{\rm RD},k}}}\right), (104)
t6\displaystyle t_{6} =𝔼{α𝐠RD0,kT𝐆RD0𝐆SR0H𝐠SR0,i𝐠SR0,iH𝐆SR0𝐆RD1T𝐠RD1,k}\displaystyle={\mathbb{E}}\left\{{\alpha{\boldsymbol{\rm g}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR0}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR0}},i}}{\boldsymbol{\rm G}}_{{\rm{SR0}}}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD1}},k}}}\right\}
=αM02M1βSR,kβRD,k2(βSR,i+M0βSR,k),\displaystyle=\alpha M^{2}_{0}M_{1}\beta_{{\rm SR},k}\beta^{2}_{{\rm RD},k}({\beta_{{\rm SR},i}+{M_{0}\beta_{{\rm SR},k}}}), (105)
t7\displaystyle t_{7} =𝔼{α2𝐠RD0,kT𝐆RD0𝐆SR0H𝐠SR0,i𝐠SR1,iH𝐆SR1𝐆RD1T𝐠RD1,k}\displaystyle={\mathbb{E}}\left\{{\alpha^{2}{\boldsymbol{\rm g}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR0}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR0}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR1}},i}}{\boldsymbol{\rm G}}_{{\rm{SR1}}}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD1}},k}}}\right\}
=α2M02M12βSR,k2βRD,k2,\displaystyle={\alpha^{2}M^{2}_{0}M^{2}_{1}\beta^{2}_{{\rm SR},k}\beta^{2}_{{\rm RD},k}}, (106)
t8\displaystyle t_{8} =𝔼{α2𝐠RD0,kT𝐆RD0𝐆SR1H𝐠SR1,i𝐠SR0,iH𝐆SR0𝐆RD1T𝐠RD1,k}\displaystyle={\mathbb{E}}\left\{{\alpha^{2}{\boldsymbol{\rm g}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR1}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR0}},i}}{\boldsymbol{\rm G}}_{{\rm{SR0}}}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD1}},k}}}\right\}
=α2M02M12βSR,k2βRD,k2,\displaystyle={\alpha^{2}M^{2}_{0}M^{2}_{1}\beta^{2}_{{\rm SR},k}\beta^{2}_{{\rm RD},k}}, (107)
t9\displaystyle t_{9} =𝔼{α3𝐠RD0,kT𝐆RD0𝐆SR1H𝐠SR1,i𝐠SR1,iH𝐆SR1𝐆RD1T𝐠RD1,k}\displaystyle={\mathbb{E}}\left\{{\alpha^{3}{\boldsymbol{\rm g}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}{\boldsymbol{\rm g}_{{\rm{SR1}},i}}}{\boldsymbol{\rm g}^{H}_{{\rm{SR1}},i}}{\boldsymbol{\rm G}}_{{\rm{SR1}}}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm g}^{*}_{{\rm{RD1}},k}}}\right\}
=α3M0M12βSR,kβRD,k2(βSR,i+M1βSR,k),\displaystyle=\alpha^{3}M_{0}M^{2}_{1}\beta_{{\rm SR},k}\beta^{2}_{{\rm RD},k}({\beta_{{\rm SR},i}+{M_{1}\beta_{{\rm SR},k}}}), (108)
t10\displaystyle t_{10} =α3M0M12βSR,i2βRD,k(βRD,i+M1βRD,k).\displaystyle=\alpha^{3}M_{0}M^{2}_{1}\beta^{2}_{{\rm SR},i}\beta_{{\rm RD},k}({\beta_{{\rm RD},i}+{M_{1}\beta_{{\rm RD},k}}}). (109)

Substituting (100) into (109), we have

𝔼{Tk,iTk,iH}\displaystyle{\mathbb{E}}\left\{{T_{k,i}T^{H}_{k,i}}\right\} =γ2(t1+t2+t3+t4+2t5+2t6CLOSE\displaystyle=\gamma^{2}\Big(t_{1}+t_{2}+t_{3}+t_{4}+2t_{5}+2t_{6}
OPEN+2t7+2t8+2t9+2t10).\displaystyle+2t_{7}+2t_{8}+2t_{9}+2t_{10}\Big). (110)

After some simplifications with i=ki=k, we can obtain (33).

B-3 Compute CkC_{k}

Similar to the calculation of BkB_{k}, We can derive the expression (35) of CkC_{k}.

B-4 Compute DkD_{k} and EkE_{k}

Similar to the calculation of BkB_{k}, we can obtain (36) and (37).

B-5 Compute FkF_{k}

Fk=γ2𝔼{|𝐠RD0,kT𝐆RD0𝔼{𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1}𝐆RD0T𝐠RD0,k\displaystyle{\displaystyle F_{k}}={\gamma^{2}}{\mathbb{E}}\left\{{\left|{{\boldsymbol{\rm{g}}}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\mathbb{E}}\left\{{{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{{\boldsymbol{\rm G}}}_{{\rm{SR1}}}}}\right\}{{\boldsymbol{\rm G}}}_{{\rm{RD0}}}^{T}{\boldsymbol{\rm{g}}}_{{\rm{RD0}},k}^{*}}\right.}\right. (111)
+α𝐠RD0,kT𝐆RD0𝔼{𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1}𝐆RD1T𝐠RD1,k\displaystyle+\alpha{\boldsymbol{\rm{g}}}_{{\rm{RD0}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD0}}}^{*}{\mathbb{E}}\left\{{{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{{\boldsymbol{\rm G}}}_{{\rm{SR1}}}}}\right\}{{\boldsymbol{\rm G}}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{*}
+α𝐠RD1,kT𝐆RD1𝔼{𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1}𝐆RD0T𝐠RD0,k\displaystyle+\alpha{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{*}{\mathbb{E}}\left\{{{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{{\boldsymbol{\rm G}}}_{{\rm{SR1}}}}}\right\}{{\boldsymbol{\rm G}}}_{{\rm{RD0}}}^{T}{\boldsymbol{\rm{g}}}_{{\rm{RD0}},k}^{*}
+α2𝐠RD1,kT𝐆RD1𝔼{𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1}𝐆RD1T𝐠RD1,k|}.\displaystyle+\left.{\left.{{\alpha^{2}}{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{T}{\boldsymbol{\rm G}}_{{\rm{RD1}}}^{*}{\mathbb{E}}\left\{{{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{{\boldsymbol{\rm G}}}_{{\rm{SR1}}}}}\right\}{{\boldsymbol{\rm G}}}_{{\rm{RD1}}}^{T}{\boldsymbol{\rm{g}}}_{{\rm{RD1}},k}^{*}}\right|}\right\}.

Using results in (86), we have

𝔼{𝐆SR1H𝐑𝐧qa𝐧qa𝐆SR1}=diag(a1,,aK),{\mathbb{E}}\left\{{{\boldsymbol{\rm G}}_{{\rm{SR1}}}^{H}\boldsymbol{\rm R}_{{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}{\boldsymbol{\rm{n}}}_{{{\rm{q}}_{\rm{a}}}}}{{{\boldsymbol{\rm G}}}_{{\rm{SR1}}}}}\right\}={\rm{diag}}\left({a_{1},...,a_{K}}\right), (112)

where an=M1βSR,n(i=1KpS,iβSR,i+pS,nβSR,n+1).a_{n}=M_{1}{\beta_{{\rm{SR}},n}}\left({\sum\limits_{i=1}^{K}{{p_{{\rm S},i}}{\beta_{{\rm{SR}},i}}}+{p_{{\rm S},n}}{\beta_{{\rm{SR}},n}}+1}\right). Substituting (112) into (111), we can obtain (37).

B-6 Compute GkG_{k}

With the help of (20), we can obtain (38) by applying similar approaches in the derivations of 𝔼{𝐑𝐱~R0𝐱~R0}\mathbb{E}\left\{{{\boldsymbol{\rm{R}}}_{{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}{{{\boldsymbol{\rm{\tilde{x}}}}}_{{\rm{R0}}}}}}\right\}.

Combining all derived terms completes the proof.

References

  • [1] V. W. Wong, R. Schober, D. W. K. Ng, and L.-C. Wang, Key Technologies for 5G Wireless Systems. Cambridge University Press, 2017.
  • [2] T. L. Marzetta, E. G. Larsson, H. Yang, and H. Q. Ngo, Fundamentals of Massive MIMO. Cambridge University Press, 2016.
  • [3] A. Yadav and O. A. Dobre, “All technologies work together for good: A glance to future mobile networks,” arXiv:1804.05963, Apr. 2018.
  • [4] A. Morgado, K. M. S. Huq, S. Mumtaz, and J. Rodriguez, “A survey of 5G technologies: Regulatory, standardization and industrial perspectives,” Digital Commun. Netw., vol. 4, no. 2, pp. 87–97, Apr. 2018.
  • [5] J. Zhang, X. Xue, E. Björnson, B. Ai, and S. Jin, “Spectral efficiency of multipair massive MIMO two-way relaying with hardware impairments,” IEEE Wireless Commun. Lett., vol. 7, no. 1, pp. 14–17, Feb. 2018.
  • [6] Q. Wang and Y. Jing, “Performance analysis and scaling law of MRC/MRT relaying with CSI error in multi-pair massive MIMO systems,” IEEE Trans. Wireless Commun., vol. 16, no. 9, pp. 5882–5896, Sept. 2017.
  • [7] S. Jin, X. Liang, K. Wong, X. Gao, and Q. Zhu, “Ergodic rate analysis for multipair massive MIMO two-way relay networks,” IEEE Trans. Wireless Commun., vol. 14, no. 3, pp. 1480–1491, Mar. 2015.
  • [8] E. Björnson, J. Hoydis, and L. Sanguinetti, “Massive MIMO networks: Spectral, energy, and hardware efficiency,” Foundations and Trends® in Signal Processing, vol. 11, no. 3-4, pp. 154–655, 2017.
  • [9] J. Zhang, L. Dai, X. Li, Y. Liu, and L. Hanzo, “On low-resolution ADCs in practical 5G millimeter-wave massive MIMO systems,” IEEE Commun. Mag., vol. 56, no. 7, pp. 205–211, Jul. 2018.
  • [10] A. Esswie, M. El-Absi, O. A. Dobre, S. Ikki, and T. Kaiser, “A novel FDD massive MIMO system based on downlink spatial channel estimation without CSIT,” in Proc. IEEE ICC, May 2017, pp. 1–6.
  • [11] X. Li, T. Jiang, S. Cui, J. An, and Q. Zhang, “Cooperative communications based on rateless network coding in distributed MIMO systems,” IEEE Wireless Commun., vol. 17, no. 3, pp. 60–67, Jun. 2010.
  • [12] H. Lee and C. G. Sodini, “Analog-to-digital converters: Digitizing the analog world,” Proc. IEEE, vol. 96, no. 2, pp. 323–334, Feb. 2008.
  • [13] J. Zhang, L. Dai, S. Sun, and Z. Wang, “On the spectral efficiency of massive MIMO systems with low-resolution ADCs,” IEEE Commun. Lett., vol. 20, no. 5, pp. 842–845, May 2016.
  • [14] L. Fan, S. Jin, C.-K. Wen, and H. Zhang, “Uplink achievable rate for massive MIMO systems with low-resolution ADC,” IEEE Commun. Lett., vol. 19, no. 12, pp. 2186–2189, Dec. 2015.
  • [15] J. Zhang, L. Dai, X. Zhang, E. Björnson, and Z. Wang, “Achievable rate of Rician large-scale MIMO channels with transceiver hardware impairments,” IEEE Trans. Veh. Techonol., vol. 65, no. 10, pp. 8800–8806, Oct. 2016.
  • [16] Y. Li, C. Tao, G. Secogranados, A. Mezghani, A. L. Swindlehurst, and L. Liu, “Channel estimation and performance analysis of one-bit massive MIMO systems,” IEEE Trans. Signal Process., vol. 65, no. 15, pp. 4075–4089, Aug. 2017.
  • [17] N. Liang and W. Zhang, “Mixed-ADC massive MIMO,” IEEE J. Sel. Areas Commun., vol. 34, no. 4, pp. 983–997, Apr. 2016.
  • [18] J. Yuan, S. Jin, C.-K. Wen, and K. Wong, “The distributed MIMO scenario: Can ideal ADCs be replaced by low-resolution ADCs?” IEEE Wireless Commun. Lett., vol. 6, no. 4, pp. 470–473, Aug 2017.
  • [19] N. Liang and W. Zhang, “Mixed-ADC massive MIMO uplink in frequency-selective channels,” IEEE Trans. Commun., vol. 64, no. 11, pp. 4652–4666, Nov. 2016.
  • [20] W. Tan, S. Jin, C. Wen, and Y. Jing, “Spectral efficiency of mixed-ADC receivers for massive MIMO systems,” IEEE Access, vol. 4, pp. 7841–7846, 2016.
  • [21] J. Zhang, L. Dai, Z. He, S. Jin, and X. Li, “Performance analysis of mixed-ADC massive MIMO systems over Rician fading channels,” IEEE J. Sel. Areas Commun., vol. 35, no. 6, pp. 1327–1338, Jun. 2017.
  • [22] T. Zhang, C. Wen, S. Jin, and T. Jiang, “Mixed-ADC massive MIMO detectors: Performance analysis and design optimization,” IEEE Trans. Wireless Commun., vol. 15, no. 11, pp. 7738–7752, Nov. 2016.
  • [23] H. Pirzadeh and A. L. Swindlehurst, “Spectral effciency of mixed-ADC massive MIMO,” IEEE Trans. Signal Process., vol. 66, no. 13, pp. 3599–3613, Jul. 2018.
  • [24] J. Xu, W. Xu, and F. Gong, “On performance of quantized transceiver in multiuser massive MIMO downlinks,” IEEE Wireless Commun. Lett., vol. 6, no. 5, pp. 562–565, Oct. 2017.
  • [25] C. Kong, A. Mezghani, C. Zhong, A. L. Swindlehurst, and Z. Zhang, “Multipair massive MIMO relaying systems with one-bit ADCs and DACs,” IEEE Trans. Signal Process., vol. 66, no. 11, pp. 2984–2997, Jun. 2018.
  • [26] P. Dong, H. Zhang, W. Xu, and X. You, “Efficient low-resolution ADC relaying for multiuser massive MIMO system,” IEEE Trans. Veh. Techonol., vol. 66, no. 12, pp. 11 039–11 056, Dec. 2017.
  • [27] J. Liu, J. Xu, W. Xu, S. Jin, and X. Dong, “Multiuser massive MIMO relaying with mixed-ADC receiver,” IEEE Signal Process. Lett., vol. 24, no. 1, pp. 76–80, Jan. 2017.
  • [28] C. Kong, C. Zhong, S. Jin, S. Yang, H. Lin, and Z. Zhang, “Full-duplex massive mimo relaying systems with low-resolution adcs,” IEEE Trans. Wireless Commun., vol. 16, no. 8, pp. 5033–5047, Aug. 2017.
  • [29] X. Jia, M. Zhou, M. Xie, L. Yang, and H. Zhu, “Optimal design of secrecy massive mimo amplify-and-forward relaying systems with double-resolution adcs antenna array,” IEEE Access, vol. 4, pp. 8757–8774, 2016.
  • [30] O. Orhan, E. Erkip, and S. Rangan, “Low power analog-to-digital conversion in millimeter wave systems: Impact of resolution and bandwidth on performance,” in Proc. ITA Workshop, Feb. 2015, pp. 191–198.
  • [31] A. K. Fletcher, S. Rangan, V. K. Goyal, and K. Ramchandran, “Robust predictive quantization: Analysis and design via convex optimization,” IEEE J. Sel. Topics Signal Process., vol. 1, no. 4, pp. 618–632, Dec. 2007.
  • [32] J. Max, “Quantizing for minimum distortion,” IEEE Trans. Inf. Theory, vol. 6, no. 1, pp. 7–12, Mar. 1960.
  • [33] P. Liu, S. Jin, T. Jiang, Q. Zhang, and M. Matthaiou, “Pilot power allocation through user grouping in multi-cell massive MIMO systems,” IEEE Trans. Commun., vol. 65, no. 4, pp. 1561–1574, Apr. 2017.
  • [34] H. Q. Ngo, L.-N. Tran, T. Q. Duong, M. Matthaiou, and E. G. Larsson, “On the total energy efficiency of cell-free massive MIMO,” IEEE Trans. Green Commun. Netw., vol. 2, no. 1, pp. 25–39, Mar. 2018.
  • [35] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge University Press, 2004.
  • [36] C. He, B. Sheng, P. Zhu, D. Wang, and X. You, “Energy efficiency comparison between distributed and co-located MIMO systems,” Int. J. Commun. Syst., vol. 27, no. 1, pp. 81–94, Mar. 2014.
  • [37] S. Cui, A. J. Goldsmith, and A. Bahai, “Energy-constrained modulation optimization,” IEEE Trans. Wireless Commun., vol. 4, no. 5, pp. 2349–2360, Sept. 2005.
  • [38] E. Lauwers and G. Gielen, “Power estimation methods for analog circuits for architectural exploration of integrated systems,” IEEE Trans. VLSI System, vol. 10, no. 2, pp. 155–162, Apr. 2002.