# Downlink Goodput Analysis for D2D Underlaying Massive MIMO Networks

The performance of downlink massive multiple-input-multiple-output networks with co-channel device-to-device communications is investigated in this paper. Specifically, we consider a cellular network with sufficient number of antennas at the base station and typical hexagonal cell coverage, where the cell users and device-to-device transmitters are randomly and uniformly distributed. To obtain the analytical expressions of system-level performance, the asymptotic signal-to-interference ratios for both downlink and device-to-device links are first obtained, which depend on the pathloss and small-scale fading of the interference channels. Since these information may not be available at the service base station or device-to-device transmitters, there exists a chance of packet outage. Therefore, we continue to derive the closed-form approximation of the average goodput, which measures the average number of information bits successfully delivered to the receiver. Hence, the system design trade-off between downlink and co-channel device-to-device communications can be investigated analytically. Moreover, the performance region in which the co-channel device-to-device communications could lead to better overall spectral efficiency can be obtained. Finally, it is shown by simulations that the analytical results matches the actual performance very well.

## Authors

• 8 publications
• 2 publications
• 176 publications
• 74 publications
• ### Downlink Performance of Dense Antenna Deployment: To Distribute or Concentrate?

Massive multiple-input multiple-output (massive MIMO) and small cell den...
10/20/2019 ∙ by Mounia Hamidouche, et al. ∙ 0

• ### Optimized Power Control for Massive MIMO with Underlaid D2D Communications

In this paper, we consider device-to-device (D2D) communication that is ...
03/06/2019 ∙ by Amin Ghazanfari, et al. ∙ 0

• ### Spectral Efficiency Analysis of Multi-Cell Massive MIMO Systems with Ricean Fading

This paper investigates the spectral efficiency of multi-cell massive mu...
08/26/2018 ∙ by Pei Liu, et al. ∙ 0

• ### Cell-Free Massive MIMO with Underlaid D2D Communications and Low Resolution ADCs

05/20/2020 ∙ by Hamed Masoumi, et al. ∙ 0

• ### Joint Transmit and Circuit Power Minimization in Massive MIMO with Downlink SINR Constraints: When to Turn on Massive MIMO?

In this work, we consider the downlink of a multi-cell multiple-input mu...
11/28/2017 ∙ by Kamil Senel, et al. ∙ 0

• ### Underlay Radar-Massive MIMO Spectrum Sharing: Modeling Fundamentals and Performance Analysis

In this work, we investigate underlay radar-massive MIMO cellular coexis...
08/04/2020 ∙ by Raghunandan M. Rao, et al. ∙ 0

• ### Average SEP-Optimal Precoding for Correlated Massive MIMO with ZF Detection: An Asymptotic Analysis

This paper investigates the symbol error probability (SEP) of point-to-p...
01/02/2019 ∙ by Zheng Dong, et al. ∙ 0

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## I Introduction

Massive multiple-input multiple-output (MIMO) is an efficient technique to boost the spectral efficiency. However, as elaborated in the existing literature, the issue of pilot contamination [1], which refers to the undiminished inter-cell interference caused by pilot reuse, may severely degrade its performance. Although there have been significant research efforts devoted to address the pilot contamination issue [2, 3] and some techniques proposed to mitigate the pilot contamination [4, 5], most of the works handle the inter-cell interference from the conventional cellular network point of view, instead of novel network topologies.

Recently, the deployment of wireless cache nodes in cellular networks has drawn significant research attentions [6, 7], where the device-to-device (D2D) links has to be introduced into cellular network. The co-channel deployment of D2D communications and cellular networks was studied in some literature [8, 9]

, where it is shown that the D2D transmission reusing the cellular spectrum may cause severe interference. To alleviate the interference, D2D underlay massive MIMO cellular networks has been proposed by exploiting spatial degrees of freedom at the base station (BS). For example, in

[10], the authors proposed a pilot reuse strategy for D2D receivers and a novel interference-aided minimum mean square error (MMSE) detector to suppress the D2D-to-cellular interference. In [11], a novel revised graph coloring-based pilot allocation (RGCPA) algorithm was proposed for pilot allocation, and an iterative scheme was adopted to minimize the transmission power of D2D links. In order to evaluate the overall effect of D2D links on cellular network, a system-level performance analysis is necessary. The interplay between massive MIMO uplink transmission and co-channel D2D transmission has been studied in [12]. However, the conclusion of its analysis cannot be directly applied in downlink. Moreover, the authors assume that the antenna number at the BS is infinity, which is not practical111In massive MIMO systems, the typical number of antennas at the BS is up to a few hundred, which is large but not infinity.. As a result, performance analysis of massive MIMO downlink transmission with co-channel D2D links is still open.

## Ii System Model

### Ii-a Network Model

We consider a typical cellular network with hexagonal cells where the radius of each cell is , as illustrated in Figure. 1. Each cell consists of a base station equipped with antennas, active single-antenna downlink users and active D2D transmitters. BSs and D2D transmitters transmit with constant powers and , respectively. The D2D links may refer to the data delivery between the downlink users and wireless cache nodes. For example, the desired data of downlink users is found in one cache node nearby, and then a direct D2D communication link is established. Since the focus of this paper is on the physical-layer SIR and throughput analysis, the establishment of D2D links is outside the scope of this work. Since massive MIMO technology is considered, is sufficiently large, e.g., a few hundred. The downlink users and D2D transmitters are uniformly and independently distributed. Without loss of generality, we investigate the performance of the first cell while other cells are all interfering cells. The -th downlink user of the -th cell is referred to as the -th downlink user.

The massive MIMO network is working in time-division duplex (TDD) mode. Thus it is assumed that the downlink channel of downlink users is estimated from their uplink pilot transmission within the same coherent fading block. Moreover, in order to improve the overall spectrum efficiency, we consider the co-channel deployment of D2D and downlink transmission, i.e., the D2D transmitters use the same spectrum as the cellular network. Note that the coexistence issue of massive MIMO uplink and D2D communications has been investigated in

[12], the focus of this paper is put on sharing the downlink transmission opportunities with D2D links. All D2D transmitters and receivers are equipped with single antenna. The -th D2D link (transmitter or receiver) of the -th cell is referred to as the -th D2D link (transmitter or receiver). The notations of downlink and D2D transmissions are summarized below.

• represent the downlink channel vectors from the

-th BS to the -th downlink user and -th D2D receiver in the -th cell, respectively. Each component of and

is complex Gaussian with mean zero and variance

and respectively. and are the pathloss from the -th BS to the -th downlink user and -th D2D receiver, where and are the distances from the -th BS to the -th downlink user and -th D2D receiver. is the pathloss exponent between BSs and users. and are the aggregation of and within one cell.

• represents the downlink channel vector from the -th D2D transmitters to the -th D2D receiver and -th downlink user in the -th cell, respectively. and are complex Gaussian with mean zero and variance and , respectively. and are the pathloss from the -th D2D transmitter to the -th D2D receiver and -th downlink user, where and are the distances from the -th D2D transmitter to the -th D2D receiver and -th downlink user. is the pathloss exponent between users. and are aggregation of and .

• is the pilot sequence of the -th downlink user. We have and where is the transmission power of each mobile user and represents the pilot length in the uplink. is the aggregation of pilot sequences from active downlink users in the -th cell. .

###### Remark 1.

The coexistence SIR analysis of massive MIMO uplink transmission and D2D transmission in [12] cannot be applied in downlink scenario, as the source of interference is completely different. In this paper, we also propose a new analytical framework to evaluate the asymptotic goodput performance with sufficiently large (but finite) number of antennas at BS . Note that the approach introduced in [12]

is for infinite M, which may not be accurate when M is only a few hundred. Moreover, we use Gaussian approximation to obtain a simple closed-form expression of the cumulative distribution function (CDF) of SIR, based on which we also derive the average goodput as the performance metric so that the potential packet outage can be counted. These results cannot be obtained with the approach in

[12].

### Ii-B Channel Model

Since D2D links share the downlink transmission opportunities, they are silent in the uplink subframe. Thus in the channel estimation phase (as illustrated in Fig. 2), the received signal of the -th BS is given by

 Ypi=(Hii)HXpi+∑∀l≠i(Hil)HXpl.

With match filter, the estimated uplink channel can be written as

 (^Hii)H = Ypi(Xpi)H/LpPu (1) = (Hii)H+∑∀l≠i(Hil)H/√Lp(ΔHii)H,

where is the channel estimate error of the -th BS. The first term is the desired CSI and the second term is the inter-cell interference due to pilot reuse (pilot contamination).

In the downlink transmission, the aggregated received signal of the users in the -th cell (as illustrated in Fig. 1) is

 Yi = HiiXdi+∑l≠iHliXdl+∑∀lUliSl, (2)

where is due to the interference from D2D transmitters. and are the downlink signals for the -th downlink user and the -th D2D receiver in the -th cell. and are the aggregations of downlink data blocks for downlink and D2D receivers in the -th cell with average power and .

Moreover, the aggregated received signal of the D2D receivers in the -th cell during the downlink subframe can be written as

 YD2Di=∑∀lGliSl+VliXdl. (3)

We neglect the effect of noise because the noise is much smaller than the interference.

Since the channel is estimated at the BS before downlink transmission, the zero-forcing precoder based on imperfect channel state information (CSI) is designed as

 Pi=(^Hii)H(^Hii(^Hii)H)−1

Denoting the aggregated downlink data block before precoding as , the downlink signal for downlink users in the -th cell is

 Yi = HiiPiWdi+∑l≠iHliPlWdl+∑∀lUliSl (4) = ∑∀lHli(^Hll)H(^Hll(^Hll)H)−1Wdl+∑∀lUliSl

and the SIR of the -th downlink user can be given by

 γi,k=1∥∥Δhii,kPi∥∥2+∑l≠i∥∥hli,kPl∥∥2+∑∀l,m∥∥ul,mi,k∥∥2PdPb, (5)

where is the channel estimate error of . Similarly, the downlink signal for D2D receivers in the -th cell can be rewritten as

 YD2Di=∑∀lGliSdl+VliPlWdl. (6)

It is assumed that the D2D transmitter has no knowledge of CSI, but the receiver has perfect CSI. Hence the SIR of the -th D2D receiver can be given by

 γD2Di,k=∥∥gi,ki,k∥∥2∑l∥∥vli,kPl∥∥2PbPd+∑(l,m)≠(i,k)∥∥gl,mi,k∥∥2. (7)

Note that the D2D transmitter does not know , it has to determine the data rate according to the statistics of . Thus goodput is introduced in the following section as the performance metric.

### Ii-C Average Goodput

It can be observed from (5) that the downlink SIR is determined by channel estimation error and CSI of interference channel , which may be unknown to the service BS. Hence the link capacity becomes random to the service BS, and it is possible that the scheduled downlink data rate may be larger than the channel capacity, which leads to the packet outage. In order to take the potential packet loss into consideration, we use the average goodput as the performance metric [13]. Given a scheduled downlink data rate for the -th downlink user, the goodput is defined as

 \emphgi,k=ri,kI(ri,k≤log2(1+γi,k)), (8)

where is is an indicator function with value if the event is true and otherwise. The the average goodput spanning all possible channel realization is given by

 ¯\emphgi,k=EH[\emphgi,k]=ri,kPr(ri,k≤log2(1+γi,k)). (9)

In this paper, we consider the downlink transmission with a target outage probability

, thus

 Pr(ri,k≤log2(1+γi,k))=1−ε.

Define the SIR Threshold determined by as , where

 Pr(γi,k≤Ti,k)=ε. (10)

The average goodput of the -th downlink user becomes

 ¯\emphgi,k=(1−ε)×L×log2(1+Ti,k), (11)

where is the number of total downlink symbols within one subband of a frame.

Similarly, for the -th D2D receiver, given the outage probability , the average goodput becomes

 ¯\emphgD2Di,k=(1−ε)×L×log2(1+TD2Di,k), (12)

where satisfies

 Pr(γD2Di,k≤TD2Di,k)=ε. (13)

## Iii Goodput Performance analysis

In this section, we first derive the asymptotic SIR expression for both downlink users and D2D receivers, and then provide the approximated expressions of goodput. First of all, we have the following lemma on asymptotic SIR.

###### Lemma 1 (Aysmptotic Downlink SIR).

When the number of BS antennas is sufficiently large, the asymptotic SIR of the -th downlink user is given by

 γ1,k=1∑∀i≠1,j[(hi1,k(hii,j)HMρii,j)2+1Lp(ρi1,kρii,j)2]+∑∀i,jui,j1,k(ui,j1,k)HPdPb (14)

and the asymptotic SIR of the -th D2D receiver is

 γD2D1,k=g1,k1,k(g1,k1,k)H∑∀i,j(vi1,k(hii,j)HMρii,j)2PbPd+∑∀(i,j)≠(1,k)gi,j1,k(gi,j1,k)H. (15)
###### Proof.

Please refer to the Appendix A. ∎

From equations (14) and (15), we can observe that the downlink performance is related to downlink users’ and D2D transmitters’ locations as well as the channel variation. Some of them are unknown to the service BS or D2D transmitters. In this case, the service BS and D2D transmitter cannot ensure that the data rate is below the channel capacity. However, if we know the distribution of the SIR, the outage probability can be controlled by setting the data rate appropriately. Therefore, we derive the CDF of downlink SIR for both downlink users and D2D receivers in the following.

###### Lemma 2 (CDF of SIRs).

Given that the number of BS antennas M and the number of interfering downlink users and D2D transmitters are sufficiently large, the CDF of the -th downlink user’s downlink SIR can be approximated as

 Pr[γ1,k≤τ]≈Q⎡⎢ ⎢ ⎢ ⎢⎣1τ−∑i≠1μi,BK−∑∀iμi,DD√∑i≠1σ2i,BK+∑∀iσ2i,DD⎤⎥ ⎥ ⎥ ⎥⎦, (16)

where the

-function is the tail probability of stantard normal distribution.

and are denoted as the expectation of and , while and are their variances.

The CDF of the -th D2D receiver’s downlink SIR can be approximated as

 Pr[γ1,k≤τ] ≈ F⎡⎣τ⎛⎝∑i≠1(μi,B′K+μi,DD)+μ1,D(D−1)⎞⎠⎤⎦λ=ρ1,k1,k, (17)

where

-function represents the tail probability of exponential distribution with expectation

and is the expectation of .

###### Proof.

Please refer to the Appendix B. ∎

###### Remark 2.

Based on the uniform distribution of cellular users and D2D transmitters, and can be calculated as follows

 μi,B = limζ→∞∫R√32R12xcos−1(√3R2x)(πζ2−3√32)R2⎡⎢⎣ρi1,kMρii,j+1Lp(ρi1,kρii,j)2⎤⎥⎦dx +limζ→∞∫ζRR2πx(πζ2−3√32)R2⎡⎢⎣ρi1,kMρii,j+1Lp(ρi1,kρii,j)2⎤⎥⎦dx,
 σ2i,B = limζ→∞∫R√32R12xcos−1(√3R2x)(πζ2−3√32)R2(ρi1,k2Mρii,j2+ρi1,k4Lpρii,j4+2ρi1,k3MLpρii,j3)dx +limζ→∞∫ζRR2πx(πζ2−3√32)R2(ρi1,k2Mρii,j2+ρi1,k4Lpρii,j4+2ρi1,k3MLpρii,j3)dx.

and can be obtained Similarly. Note that the distance between D2D transmitters and receivers is much smaller than the cell radius, the distribution of D2D users can also be considered uniform, thus , and , can be calculated.

As a result, we have the following theorem on average goodput of downlink users and D2D receivers.

###### Theorem 1 (Average Goodput).

With outage probability , the average goodput of the -th downlink user is given as

 ¯\emphg1,k=(1−ε)×L×log2(1+1Q−1(ε)σC1+μC1), (18)

where

 σC1=√∑i≠1σ2i,BK+∑∀iσ2i,DD,μC1=∑i≠1μi,BK+∑∀iμi,DD.

Similarly, for the -th D2D receiver the goodput data rate can be presented as

 (19)

where . Then the overall average goodput of the -st cell is given as

 ¯\emphg1=KEU[¯\emphg1,k]+DEU[¯\emphgD2D1,k], (20)

where and are expectations based on the distribution of locations for single downlink user’s average goodput and single D2D receiver’s average goodput.

###### Proof.

Theorem 1 is directly derived from lemma 2. ∎

## Iv Simulation Results

In the simulation, we consider a network with 19 hexagonal cells each with radius . The number of downlink users in each cell is 10, and the antennas number of BS is 250. The length of pilot sequence is and the length of downlink symbols within one frame is . The pathloss exponents are and . The transmitting power of each BS and each D2D transmitter is 46 dBm and 23 dBm.

In Fig. 3, there are 10 D2D transmitters in each cell and the transmission distance of each D2D link is . For both downlink users and D2D receivers, we give the numerical results as well as the asymptotic results. It can be observed that the average goodput of downlink users decreases with respect to its distance to the service BS. This is because of stronger inter-cell interference. On the other hand, the average goodput of D2D links increases with respect to its distance to BS. This demonstrates the impact of interference from BS. Even there are a large number of antennas at the BS, its interference to D2D receiver is still strong and dominant.

Fig. 3 shows that the analytical results fit the actual performance quite well. Thus we can study the effect of D2D number with analytical expressions instead of complicated simulations. In Fig. 4, the effect of D2D number is demonstrated through analytical results, which show that the cell overall average goodput increases almost linearly with the D2D number, while the average goodput of all downlink users declines. Thus we reveal the trade-off between cell overall average goodput and downlink users’ average goodput. Therefore, to achieve decent cell overall average goodput and fair average goodput for downlink users simultaneously, the number of D2D links per cell can be determined through our analytical results.

## V Conclusion

In this paper, we focus on the average goodput performance analysis of a D2D underlay downlink massive MIMO system. It is assumed that the cell coverage is hexagonal and the distribution of downlink users and D2D transmitters is independent and uniform. The asymptotic SIR expressions for both downlink transmission and D2D links are firstly derived. To take the potential packet outage, we continue to derive the approximated expressions of average goodput for both downlink and D2D links, which measures the number of information bits successfully delivered to the receiver. Based on it, the trade-off between the two types of links can be studied. Through simulation we show that the analytical results fit the numerical results quite well.

## Appendix A Proof of Lemma 1

In massive MIMO systems, when is sufficiently large, we usually utilize asymptotic orthogonality of channel as

 hij,k(hij,k)HM→ρij,k,hij,k(hlm,n)HM→0,(i,j,k)≠(l,m,n)

It is also the same for . Thus the interference term caused by neighboring BSs in Equation (4) can be simplified as

 ∑∀l≠iHliPlWdl≈∑∀l≠iHli(^Hll)HRl−1Wdl ≈ ∑l≠iHli(Hli)H(Xpi(Xpl)HPuLp)Rl−1Wdl+∑l≠iHli(Hll)HRl−1Wdi.

The first term remains because . The second term remains because . Similarly, for D2D receivers in the -th cell, the interference from neighboring BSs can be simplified as

 ∑∀l≠iVliPlW≈∑l≠iVli(Hll)HRl−1Wdi.

Therefore, the received signal of the -th downlink user can be presented as

 y1,k = √Pbwd1,k+∑∀i≠1,j√Pbρi1,kxp1,k(xpi,j)Hρii,j−1wdi,j/(PuLp) +∑∀i≠1,j√Pb(hi1,k)Hhii,jwdi,jMρii,j+∑∀i,j√Pdui,j1,ksdi,j,

where the intra-cell interference for downlink users is neglected for it is much smaller than inter-cell interference. Similarly, the received signal of the -th D2D receiver is

 yD2D1,k≈∑∀i,j√Pb(vi1,k)Hhii,jwdi,jMρii,j+∑∀i,j√Pdgi,j1,ksdi,j.

In this paper, we assume the transmit signals are Gaussian with zero mean and unit variance. Then asymptotic SIR expressions as (14) and (15) can be straightforward obtained.

## Appendix B Proof of lemma2

We first consider the downlink users. The asymptotic SIR expression is presented in (14). In the dominator, are independent variables with expectation and variance . Meanwhile, are also independent variables with expectation and variance . When the number of interfering cellular users and D2D transmitters is sufficiently large, we have and

, which satisfies requirements of the Lindeberg-Feller Central Limit Theorem

[14]. According to the theorem, converges to a Gaussian variable when is sufficiently large. Thus, we apply Gaussian approximation on . Then the probability that the downlink SIR of the -th cellular user is less than (coverage outage probability) can be written as

 Pr[γ1,k≤τ]≈Q⎡⎢ ⎢⎣1τ−μi,BK−μi,DD√σ2i,BK+σ2i,DD⎤⎥ ⎥⎦.

The probability that the SIR of the -th D2D link is less than can be presented as

 Pr[γ1,k≤τD2D]=Pr[g1,k1,kg1,k1,kH−τD2D(IBS1,k+ID2D1,k)≤0],

where , . Due to D2D receivers are quite close to their associated transmitters, the variation of the interference is relatively much smaller than the signal variation. Thus we can regard the interference as a constant. Since the downlink channel is a complex Gaussian variable, the signal power is an exponential variable, where the coverage outage probability can be simplified as

 Pr[γ1,k≤τD2D]=Pr[g1,k1,kg1,k1,kH−τD2DE[IBS1,k+ID2D1,k]≤0]

Thus (2) can be easily obtained. This finishes the proof.

## References

• [1] T. L. Marzetta, “Noncooperative cellular wireless with unlimited numbers of base station antennas,” IEEE Transactions on Wireless Communications, vol. 9, pp. 3590–3600, November 2010.
• [2] F. Fernandes, A. Ashikhmin, and T. L. Marzetta, “Inter-cell interference in noncooperative TDD large scale antenna systems,” IEEE Journal on Selected Areas in Communications, vol. 31, pp. 192–201, February 2013.
• [3] M. A. Teeti, R. Wang, Y. Liu, and Q. Ni, “Pilot optimization in multicell massive mimo,” in 2016 IEEE International Conference on Communication Systems (ICCS), pp. 1–6, Dec 2016.
• [4] R. Wang, Y. Chen, and H. Tan, “Data-assisted channel estimation for uplink massive mimo systems,” in 2014 IEEE Global Communications Conference, pp. 3766–3771, Dec 2014.
• [5] R. Wang, Y. Chen, H. Tan, and Q. Zhang, “Data-assisted massive mimo uplink transmission with large backhaul cooperation delay: Scheme design and system-level analysis,” in 2015 IEEE Global Communications Conference (GLOBECOM), pp. 1–7, Dec 2015.
• [6] S. W. Jeon, S. N. Hong, M. Ji, G. Caire, and A. F. Molisch, “Wireless multihop device-to-device caching networks,” IEEE Transactions on Information Theory, vol. 63, pp. 1662–1676, March 2017.
• [7] N. Golrezaei, A. F. Molisch, A. G. Dimakis, and G. Caire, “Femtocaching and device-to-device collaboration: A new architecture for wireless video distribution,” IEEE Communications Magazine, vol. 51, pp. 142–149, April 2013.
• [8] X. Wang, X. J. Li, H. Y. Shwe, M. Yang, and P. H. J. Chong, “Interference-aware resource allocation for device-to-device communications in cellular networks,” in 2015 10th International Conference on Information, Communications and Signal Processing (ICICS), pp. 1–5, Dec 2015.
• [9] X. Luan, J. Wu, C. Piao, Y. Cheng, and H. Xiang, “Cooperative transmission based on multi-relay device-to-device communications in cellular networks,” in 11th International Conference on Wireless Communications, Networking and Mobile Computing (WiCOM 2015), pp. 1–5, Sept 2015.
• [10] X. Liu, Y. Li, X. Li, L. Xiao, and J. Wang, “Pilot reuse and interference-aided MMSE detection for D2D underlay massive MIMO,” IEEE Transactions on Vehicular Technology, vol. 66, pp. 3116–3130, April 2017.
• [11] H. Xu, N. Huang, Z. Yang, J. Shi, B. Wu, and M. Chen, “Pilot allocation and power control in D2D underlay massive MIMO systems,” IEEE Communications Letters, vol. 21, pp. 112–115, Jan 2017.
• [12] X. Lin, R. W. Heath, and J. G. Andrews, “The interplay between massive MIMO and underlaid D2D networking,” IEEE Transactions on Wireless Communications, vol. 14, pp. 3337–3351, June 2015.
• [13] R. Wang and V. K. N. Lau, “Cross layer design of downlink multi-antenna OFDMA systems with imperfect csit for slow fading channels,” IEEE Transactions on Wireless Communications, vol. 6, pp. 2417–2421, July 2007.
• [14] Y. S. Chow and H. Teicher, Central Limit Theorems, pp. 290–323. New York, NY: Springer US, 1978.