Massive multiple-input multiple-output or massive MIMO system has great potential for 5th generation (5G) wireless communication systems as it is capable of providing game-changing enhancements in area throughput and energy efficiency (EE). This work proposes a realistic and practically implementable EE model for massive MIMO systems while a general and canonical system model is used for single-cell scenario. Linear processing schemes are used for detection and precoding, i.e., minimum mean squared error (MMSE), zero-forcing (ZF), and maximum ratio transmission (MRT/MRC). Moreover, a power dissipation model is proposed that considers overall power consumption in uplink and downlink communications. The proposed model includes the total power consumed by power amplifier and circuit components at the base station (BS) and single antenna user equipment (UE). An optimal number of BS antennas to serve total UEs and the overall transmitted power are also computed. The simulation results confirm considerable improvements in the gain of area throughput and EE, and it also shows that the optimum area throughput and EE can be realized wherein a larger number of antenna arrays at BS are installed for serving a greater number of UEs.

In existing 3G and 4G standards, BS allows only up to 8 antenna ports. However, the BS antennas for 5G standard are increased to support data transmission rate up to several GBs/seconds. The enhancements are forecasted to be achieved through the techniques, i.e., network densification (extra BS nodes), amplified bandwidth (mm-wave spectrum) or by using massive MIMO systems. Massive MIMO systems are multicarrier systems with “L” cells that use time division duplex operation protocol. In this system, to realize channel hardening, BS is prepared with number of antennas where each BS communicate with number of UEs instantaneously. Realizing unbounded EE is impossible because the system model does not consider the power expended by analog circuits (for radio frequency (RF) and baseband processing) and signal processing that raises linearly with and . Hence, it can be taken as constant only in massive MIMO setups where values of and are comparatively small, while its changeability can be observed in massive MIMO systems in which, (). The definition of massive MIMO in [1] also supposed the ration of , while in [2, 3], the ratio has been taken as a small constant value. All BS process its signals independently by using linear transmit precoding and received combining. We have considered a canonical massive MIMO system model that is implementable in real-time testbeds. Massive MIMO system shows many advantages over MIMO systems low power consumption because of enlarged antenna aperture and improved network capacity [2]. So far, efficient modeling of power consumption is a protuberant apprehension of these systems. Especially, in mobile UEs, the battery technology is not improving as per the massively rising demand for multimedia communication [3]. Analyzing a massive MIMO network for realistic modeling of total power consumption during uplink and downlink communications as it plays a vital role to achieve optimal area throughput, and maximum EE is a prominent aim of this work. Area throughput and EE can be defined as Definitions 1 and 2, respectively.

Definition 1. Total bits successfully transmitted per kilometer in one second is called area throughput measured as bits-per-Joule [4].

Definition 2. Total bits successfully transmitted by consuming a Joule of energy is called EE measured as bits-per-Joule [4]. In massive MIMO systems, EE is dependent on many factors, i.e., spectral efficiency, network architecture, power consumption by the entire system, and transmission protocol [49]. Existing literature provides mathematical modeling of mentioned parameters to enhance area throughput and EE. In [10], authors have provided an exhaustive power consumption modeling for EE optimization. Moreover, the circuit power utilization impact on transmitting antenna is deliberated by authors in [1117]. Specifically, [13] emphasizes on uplink communication and power consumption of the massive MIMO network in uplink communication. Authors determined the optimized EE while the antennas of UEs are turned off. In [1416], authors examine the downlink communication of a system and concluded EE as a concave function of total antennas at BS. In [1721], authors provide the optimum value of to attain the maximum EE of a multicell scenario for a specified number of UEs. Authors in [2226] modeled the EE by efficient power allocation for uplink and downlink by using MMSE processing scheme. Realization of infinite throughput and EE is perceptibly unrealistic because the power required for signal processing and analog circuits increases with increasing values of [913]. However, existing literature claims that it is possible to attain infinite EE because most of the modeling has been done by neglecting the impact of an increasing number of UEs.

The primary concentration of this paper is to mathematically model the impact of , , and the total power utilization on EE of a massive MIMO architecture by using linear processing schemes, i.e., MMSE, MRT/MRC, and ZF. We have also illustrated how energy hoarding help to improve the area throughput and EE at the network level. The second objective is to deliver an accurate power dissipation model that encompasses the power consumed by the entire system during uplink and downlink communication. Additionally, the linear precoding and received combining schemas are used including MMSE, ZF (perfect/imperfect CSI), and MRT/MRC at BS and UEs. The resulting expression (38) gives understandings of the impact of atmospheric effects during propagation, structural parameters, and the hardware apparatuses cast off to originate the power utilization modeling. Prescribed objectives are well achieved and outcomes are summarized as follows. (i)MMSE and ZF schemes carry the max EE gain; however, ZF is better as compared with min square root error and MRT because of less complexity and better interference mitigation practice(ii)Total circuit power and transmit power increase by increasing BS antennas(iii)BS has hundreds of antennas that is an utmost efficient way to serve UE with optimal EE and area throughput(iv)The ratio of reducing cell size is also applied as provided in [26] expression (40). According to the expression cell, radius and user density values are set accordingly. Reducing cell radius decreases the network capacity of the system; however, it increases the EE.

Listed results are appropriate to authenticate the massive MIMO architecture that it is an accurate method to attain the maximum EE in 5th generation cellular networks. The remaining sections of this paper are structured as follows. Uplink and Downlink Massive MIMO System Model illustrates the canonical massive MIMO system model for both uplink and downlink communication. Additionally, it also describes the user distribution and linear precoding schemes. Later on, Methodology and Calculations is the main section as it includes all significant mathematical modeling and computations of data rate and power utilization model. Optimization of , , power, area throughput, and EE is articulated in the same Methodology and Calculations. Moreover, Results and Discussion stretches results to validate the theoretic analysis and draws the contrasts of the proposed model by using prescribed linear processing schemes with existing works. Finally, key insinuation conclusions are drawn in Conclusions. Table 1 shows the symbolic presentations used in this paper.

The following segment represents the specifications of massive MIMO network covering the distribution of users, linear processing scheme, and channel model. A single-cell case is cogitated in which numbers of transmit antennas are employed to attend single antenna UEs. Furthermore, the flat-fading channel is considered for transmission which is delimited by coherence time-frequency blocks as characterized by . Time division duplex is used for an operation mode because the BS and UEs are copiously synchronized to each other. Figure 1 provides a demonstration of TDD coherence block frame for bidirectional transmissions.

The uplink pilots are used at the BS to estimate the channels of UEs, whereas the downlink pilots are used at UE for estimating the interference and at the channel by using linear precoding schemes [11]. Prearranged channel massive MIMO model is considered because it enables the allocation of antenna correlation and path loss component. Additionally, it permits considerable antenna correlation of antenna either due to less antenna spacing or the absence of scattering f4or large antenna arrays [21, 22]. Various types of antenna configuration are proposed for massive MIMO systems in the literature including cylinder configuration, linear configuration, and rectangular configuration [2325]. We have considered the cylinder configuration. Moreover, MMSE, MRT, and ZF processing and detection techniques are operational in these networks that help to reduce power consumption significantly; however, the overall circuit power dissipation upsurges progressively by growing [19]. Received vector at BS is denoted with , given as

The matrix between BS and UE is specified by in equation (2). The channel coefficients among UE and BS antenna are identified as (). The , where ( represents the average dispersed power) represents the () symbols vectors that are transmitted simultaneously by UEs. Moreover, the is the Additive White Gaussian Noise (AWGN) vector. The written as , l ranges from {}, and is the flat-fading coefficient from UE and BS antenna. It symbolizes the signal fading because of obstacles such as large building and moving vehicles, and also due to the extension distance signals propagation from BS to UE, where represents channel vector given as . The () are the channel vector entries that define the propagation channel among the UE and BS antenna. is the shadow fading and attenuation coefficient and signifies the identity matrix. The three precoding schemes including MMSE are used to carry out the downlink tasks of uplink data detection and data precoding. The LS-MIMO transmitter reduces interuser interference by using the precoding like spatial signal processing methods. Supposedly, the BS is capable of obtaining perfect channel information by uplink pilot sequences. Therefore, the results of ZF are gotten for single-cell scenario and also (with perfect CSI) for MMSE and MRT/C. The represents the uplink linear received matrix in which the the column is assigned to the UE. As the ZF, MMSE, and MRT/MRC schemes are employed, therefore, the precoding matrix classification in uplink can be taken as in (3) [20]:

In (3), represent the Hermiston matrix and inverse transpose, individually. Further, includes the channel values for all the UEs, and the power vector for uplink is . The ZF precoding is capable of attaining the efficiency and the capacity of the network, if receiver has exemplary channel information (CSI) for a large number of UEs. Contrarily, the defective CSI at the BS the completion of ZF declines sustainably to the precision of the channel information [24]. To get the full multiplexing gain, ZF demands immense overhead addressing the signal-to-noise ratio (SNR). Defective CSI outcomes as of the turn-out loss due to interuser interference. The interference continues until the users are unable to form an independent beam for independent communication. For downlink communication, linear processing schemes are characterized as . The matrix is modeled as in (4):

According to the equation [19], to reduce the calculus complexity. The MRC is basically a diversity combining precoding scheme that combines the signals from each channel and then adds them, and alternatively, it is known as a predetection combining scheme [27]. For independent AWGN, it is used as a combiner and also reinstates a signal to its original form. The user terminal (UE/UTs) is actually located on () measured in meters.

All users are randomly distributed in the cell. The is calculated w.r.t BS with the Round Robin (RR) selection; UEs are wisely selected for the communication, and the UE locations are reserved with arbitrary variables from the user distribution as shown in Figure 2. The distances from BS to UEs are comparatively greater than the range among BS antennas of the array configuration. Suppose that () denotes the max distance of a UE and () represents the min distance a UE has from BS. All UEs are spread around the BS evenly. Further, UE’s location is designated by modeled as (5) [19].

Additionally, represents the large-scale fading given as . path loss exponent and constant symbolizes the attenuation factor in the channel at . Lastly, in modeling, the invers--channel attenuation can be calculated as follows [19]:

3. Methodology and Calculations

By applying MMSE, ZF, and MRT/C precoders, the improvement in EE of a massive MIMO system and the modeling of the optimized power consumption are the prominent objectives of this work. We have explored the methodology of energy consumption in which energy is saved by installing low-power BS according to the traffic requirements. This is a distinguished approach to reduce the consumption of power, by cutting the infrastructure cost and managing the spectrum efficiently [20]. The power consumed at the BS depends upon the total number of dynamics transmit antennas against the number of UEs obliged simultaneously.

Calculations performed in this section are shown in Figure 3. The uplink and downlink communication data rates are calculated in the first and second steps. Moreover, the right side of the flow diagram displays computations of power consumption by power amplifiers and circuitry at the BS and UEs. Further, the data rate totaled as an average. The rate in the step of total rate and the power consumption calculations are added to find the total power expended.

Lastly, the seventh step delivers the ration of average data rate and the overall power dissipation used to compute EE. Subsequent portion provides the detailed computations given in Figure 3.

3.1. Average Data Rate

The upper left portion of Figure 3 contains two steps for the computation of the average rate. Conferring to (7), the average rate is a summation of achievable data rate downlink and in uplink communications. We aimed to calculate the total achievable data rate for UE is given as where is achievable data rate in uplink and represents the downlink data rate. By the Gaussian codebooks, the uplink achievable data rate is given as (8) [3]: where is the uplink transmission fraction and represents the uplink pilot length. Moreover, is the pilot overhead value, and the gross rate (excluding overhead factors) is given as , where is calculated as (9): where shows the noise power and is the uplink power allocation. In [22], the power distribution vector in uplink is , and () are elements of given as (10):

According to [21], if zero-forcing precoder is applied where , and , where is proportional to the received SINR that is a designed parameter. According to the Gaussian codebooks, the downlink rate for perfect CSI is given as in (11)L where is the downlink transmission fraction and is the downlink-pilot length; the pilot overhead can be modeled as , and the gross rate is given as . The value of is (12): where is the downlink power vector and () is the component of . Fraction values of are cast off to compute the realizable rate for forward link.

The results of (8) and (9) are plugged in (7) to calculate the achievable rate for as (14), where is taken as the gross rate, and it is computed by totaling up the achievable data rates for bidirectional communications.

3.2. Computations of Power Consumptions

The total power consumption is the summation of power dissipated by a power amplifier (PA) and also consumed by circuit components (of UEs and BS) [24], as our goal is the proposition of a power distribution model to elevate the EE based upon data rates. It is stimulating because (1) of the difficulty level of associated iterative algorithms and (2) uncontrolled positioning and a massive number of BS antennas that make it difficult to define an unapproachable cluster (set of BS antennas) without focusing to the external interference [25]. Accordingly, we have executed a more reliable, more detailed, and realistic model to calculate the power consumption of the massive MIMO system that includes all aspects of energy consumption including transmit power and component power consumption at UEs and BS, modeled as (15):

In (15), is the power consumption by PA and is represented as circuit power consumption, correspondingly. In previous works, the calculations of total circuit power are considered as (), where the is taken as a constant value of power consumption calculated for load-independent backhaul, control signaling, and the power required for cooling the system and baseband processor [4, 17]. This way of calculating circuit power is mistaken. As a system in which the value of grows consequently, the circuit power also shows proportional growth with (that is taken constant in existing literature). The existing circuit power consumption claims that achieving an infinite EE is conceivable by increasing ; however, it is impossible to attain infinite EE. An average power consumed by power amplifier is delineated by [19] as . Henceforth, in accordance with the power consumption by PA during uplink, transmission is intended as in (16):

Similarly, the downlink PA’s power can be modeled as (17)

The overall power consumption by PA is calculated as that is the sum of (16) and (17). The () represents the efficiency of PA given as .

The is the sum of power used in analog components and by several signal processing [20]. Mathematical model of total circuit power consumption is provided (19):

In (19), represents the power consumed in transceiver chains. It can be computed as a summation of power needed for the BS antenna array, UEs, and synthesizers. Additionally, the power consumption while coding/decoding () can be modeled as

The load-dependent backhaul is represented by (used for data transmission between BS and core network), and it is directly proportional to the average sum data rate. Hence, it can be further modeled as the sum of average rate during uplink and downlink communications as (21), whereas is required by backhaul traffic power (measured in Watt/bits/sec).

We have considered coherence block/seconds, and the estimation of the channel is also executed once for each coherence block. BS performs channel estimation by getting () matrix from UE and multiplying it by pilot-sequence length [3]. Moreover, represents the power used in the estimation of channel state given as and further modeled in (22)

Moreover, is used to represent the power disbursed while applying linear processing schemes at BS. Two categories of power are measured in linear processing: the first type is required while calculating and , and the other is needed for matrix-vector multiplication by each data symbol [23]. The precoding/decoding matrices are calculated once in every coherent block while its complexity is contingent on the type of scheme. Hence, for MRT/MRC, power consumption is given as (23)

Moreover, the power consumed for ZF is approximately modeled as (24) where and are computational efficiency symbols of UE and BS [19]. Lastly, the power consumed while using MMSE is modeled as (25): where is supposed as a fixed power required for baseband, cooling, and control signaling.

3.3. Optimization of , , and

The optimized value of can be symbolized as () and further intended as (26): where the values and is specified as PA’s power consumption calculated as in (18), where and are power consumed by components at UE and BS given as (27) and (28):

The theorem in (26) shows the selection by the prescribed system whereas (26) displays that the optimal value is independent of rate-dependent power including the power required for the backhaul traffic and the power consumed in channel coding/decoding. Conversely, it increases by increasing the rate-independent powers such as power required by synthesizers, components of UE, and fixed power for control signaling. The lower bound of optimal can be calculated as (29):

The optimized transmit power is represented by () and calculated by using the Lambert function in (30) where and are given in (31) and (32), respectively.

The is also independent of the rate-dependent power defined previously, such as power consumption during coding and decoding and power required by the backhaul traffic; however, it grows with fixed power and also with the power required by UEs and BS components for site cooling and control signaling, etc. If the large circuit powers are used, in that case, the greater PA power can be afforded by the network, meanwhile has a slight influence on the overall power dissipation. The authors in [6, 10, 11] illustrate that the TDD-based system allows a power drop that can be directly proportional to ( or imperfect CSI), however, keeping the nonzero data rates for the value that approaches to infinity. The lower bound of can be intended as (33).

It has been intended by finding the roots of the following polynomial of (35) [16]. The calculation of optimal “” for the optimal value of EE is our target. Hence, to provide critical tractability, let us suppose the summation of (“SINR” and in that way the power of PA) and the total numbers of BS antennas for each UE, retained as constants (N) and with () and (). Therefore, the gross rate value can be taken as , and the optimal value of the number of users can be taken as

In (34), describes the positive real roots of (35). where

As described earlier, the optimal “” is the real roots of the polynomial specified in (35). Quartic polynomials have four roots, and there are basic closed-form root expressions for a polynomial [27]. Nevertheless, these closed-form expressions are overlong, that is why these are seldom used, and also, as an alternative to these expressions, many simple algorithms are used to calculate the roots with greater statistical precision [27]. The number of UEs () are affected by different parameters, to check the impact; suppose and are both nearly equal to 0. This case is mainly relevant as and , basically, diminish with the computational efficiencies that are predictable to rise speedily in the future. For the prescribed case of supposing two types of power consumption equal to zero, the will become (38):

The equation is a decreasing function of the power needed by BS and UE circuit components that are always growing by increasing the values of and . The previous subsections deliver modest closed-form expressions of , , and that enable us to get the maximization in the area throughput gains and EE. The eventual objective is to calculate joint-global optimum. As and are numerals, so the global optimum value can be found by a comprehensive exploration of overall equitable combinations of the pair () and evaluate the optimal power appropriation. In extensive MIMO systems, there is a trade-off between total power consumption and achievable rate [22].

3.4. Modeling of EE

According to the definition of EE, it is in a bit-per-Joule [2] and computed as a ratio of avg data rate and average power dissipation. In massive MIMO networks, metric calculation of total EE holds the form of (40). It depicts a decrease in overall power used by the system that carries enhanced results, and it can be accomplished by dropping the static power consumption of the system. The EE for the prearranged system model can be evaluated by applying the proposed power consumption model. EE can be modeled as the following optimization problem.

An optimized value of EE single-cell scenarios is calculated by the solution of (39). We have plugged in all previously computed values, i.e., the sum rate in (14), the power consumed by the circuit components (16), and the power of power amplifiers as in (18) to outline the EE in (40).

The relationship shows a trade-off between total power dissipation and data rate [21]. EE can be improved by increasing the antennas where the total power dissipates, and it also grow proportionally.

4. Results and Discussion

This section provides the performed simulation results and discussion on the predefined model. We have considered an appropriate number for the EE calculations with ZF, MMSE, and MRC/MRT precoding schemes. The UEs are dispersed inconstantly; thus, the large-scale fading is demonstrated using where values are given in Table 2.

The efficiency of PA for BS and UEs is given as 0.4 and 0.3, respectively. The remaining parameters used for simulation purpose are listed in Table 1, consistently stirred by the previous task in [19]. We have examined an appropriate approach to model power consumption illustrated in Methodology and Calculations. After the initialization of simulation parameters provided in Tables 1 and 2, Algorithm 1 provides the simulation algorithm that gives computations of flow chart given in Methodology and Calculations and optimized values of the parameters to draw the figures.

Proposed algorithm:
Step 1: According to Table 1, initialize the simulation parameters
Step 2: Calculate the orthogonal pilot sequences and inverse-channel attenuation
Step 3: If
     MRT = true || MMSE = true
Step 4: For
     M =1: max;
     Calculate RF power; Circuit Power; SINR and opt EE
Step 5: if ZF algorithm = true
     Compute circuit power; RF power; opt EE; and SNIR
      opt, opt, and optimum power.
Step 6: Plot of figures
   EE for single and multicell; RF; antenna power, and area throughput.
4.1. EE Performance Comparison

In Figure 4, the results of EE vs. transmit antennas in the single-cell (MRT/C with perfect CSI, ZF with perfect and imperfect CSI, and MMSE with perfect CSI) are shown. By comparing these results with the results of previous work in [16, 19], the purposed model demonstrates improved performance.

Table 3 demonstrates the contrast of all the outcomes in which the peak value of EE contrary to the number of BS antennas in a single-cell scenario (with perfect/imperfect channel state information) gives 15% improved results with a smaller . The results publicized in Figure 4 also portray that we can be able to make a more energy-efficient system by minimizing the cell size and .

Moreover, the PA’s power dissipation that exploits EE for optimal with is given in Figure 5. The MMSE has power dissipation of 50 Watt, and ZF with perfect and imperfect CSI is 55 and 70 Watts, respectively, whereas for MRT/MRC, the value is 5 Watt that is quite small as compared to the MMSE and ZF. For the ZF scheme, increasing the transmit power by the growing value of is an extra energy-efficient approach.

It is similar to Corollary 8 in [19]; however, it is dissimilar as compared with the results of [9, 10] that claim that transmission power is decreased with decreasing . Further, the power consumed by an antenna at BS is inversely proportional to as shown in Figure 6. According to this figure, ZF is 100 (m-Watt/antenna) and the transmit power with MSSE where MRT/C is responsible is about 200 (m-Watt/antenna) that is lesser than calculated by [21].

4.2. Results for Area Throughput

Figure 7 reveals a real-time advancement in the area throughput as compared with existing literature. The maximum of the improvement is achieved in the case of MMSE and ZF with imperfect CSI.

It approves that the massive MIMO with applicable interference mitigation linear precoding schemes can achieve the outstanding area throughput and EE. In sharp contrast, the deployments of a massive number of transmit antennas while using MRT/C processing are severely restraining for both tasks either the EE or the area throughput.

Comparing the numerical values of plots in Figure 7 and the results of the previous work, the given plots demonstrate that our proposed model gives noteworthy progress in terms of area throughput. Figure 8 provides a 3D plot to show the results of an appropriate number of antennas and users equipment to realize optimal values of EE. The optimal EE is achieved at that serves around UE efficiently. Table 4 shows the comparison of area throughput of our work existing literature. The improvement of area throughput can be significantly observed from the generated numerical results.

5. Conclusions

In this paper, we have examined massive MIMO systems to get optimal EE and the area throughput gains. As a first step, we have computed the average data rates of the system for uplink and downlink communications and then suggested an essentially applicable, less complex, and an energy-efficient power consumption model. We have recognized an exceptionally optimal power allocation arrangement that can achieve considerable improvement in EE and area throughput in a single-cell scenario. The simulation results reveal noteworthy implications.

The study was fundamentally founded upon uncorrelated fading, where every user can have an exclusive channel covariance matrix. The results of MMSE and MRT/MRC are produced by the Monte Carlo simulations. Although the MMSE is the optimum for throughput gain, it is observed that ZF processing provides higher efficiency gain. It is because of the complex computations of MMSE; however, the variance in results is relatively small. The MMSE scheme has the benefit of handling the system where the value of . The ZF precoding (with imperfect CSI) has the same performance as MMSE and ZF (with perfect CSI). The results with MRT/C are not that much good but it operates under high IUI that is why the rate/UE is less. In this scheme, complexity of signal processing is also lesser than the other schemes for the same values of and ; however, the power saving is not large enough to recompense the low data rates. In order to attain the same rates as ZF, the MRT/MRC needs a system where the values of . Conversely, it will significantly increase the circuit power dissipation.

It is predicted that the circuit consumption will reduce with time, indicating that the higher EE can be achieved by using fewer UEs, smaller number of BS antennas, reduced power, and progressive processing schemes. The findings of this work specify that the massive MIMO system can be built by employing a low-power transceiver (consumer grade) apparatus at coverage-tier BSs as an alternative for conservative industrial-grade high-power consuming apparatus. The decrement in power dissipation brings improved results, and it is possible by abbreviating the circuit power consumption, as the results are demonstrated in Figure 4. The ratio of reducing cell size is also applied as provided in [26] expression (40). Reducing cell radius decreases the network capacity of the system; however, it increases the EE. Conclusively, we can say the reduction in the max distance of a cell helps to expand the EE, and the reason is when we reduce the cell radius essentially, we are reducing the transmit power that results in terms of improved EE of a massive MIMO system.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare no conflict of interest.