MU-MIMO Downlink Capacity Analysis and Optimum Code Weight Vector Design for 5 G Big Data Massive Antenna Millimeter Wave Communication

1School of Electrical and Electronic Engineering, North China Electric Power University, Beijing 102206, China 2National Key Laboratory of Electromagnetic Environment, China Research Institute of Radiowave Propagation, Qingdao 266107, China 3School of Information and Communication Engineering, Beijing University of Posts and Telecommunications, Beijing 100876, China 4College of Computer Science, Zhejiang University, Hangzhou, Zhejiang 310027, China


Introduction
Millimeter wave (mmWave) communication which explores shorter propagation distance in frequency band of 20-40 GHz is a key enabler for the fifth generation (5G) mobile communication systems [1].The mmWave communication also provides significant benefits to a variety of applications such as vehicular communication, wire-able networks, and autonomous robots [2].In downlink transmission for a MU-MIMO system, a base station (BS) serves multiple users simultaneously in the same frequency and time slot.Thus the throughput can be enhanced by spatial multiplexing.However, the cochannel interference (CCI) becomes a dominant factor in capacity due to nonorthogonal signalling [3].Hence, minimising CCI as much as possible at the end users is of importance.CCI can be suppressed by using linear precoders and decoders at transmitter and receiver sides [4].Additionally, the channel state information (CSI) as well as partial information on the transmitter side can be used to improve the system's performance [5].For all users served by the BS, the CSI should be known in advance at the BS to support the CCI minimisation which is not available at user's end.Explicitly, the overhead on the system is reduced in this approach because the channel information feedback is not required [6,7].To perfectly cancel CCI at each end user, the restriction on the system configuration is necessary.For instance, the number of antenna elements at the BS should be larger than the total number of antennas at the end users [7].This assumption is usually valid for single antenna per user which is generally adopted in the MU-MIMO system.
There are several studies that have been undertaken in MU-MIMO downlink systems on how to minimise the CCI problem related to CSI and analyze the achievable sum-rate per user as well as the average system capacity.In [8], an in-depth capacity analysis for nonorthogonal multiple access (NOMA) mmWave massive MIMO systems was provided.A simplified mmWave channel model was also explored by extending the uniform random single-path (UR-SP) model with the angle-of-arrival (AOA).Furthermore, the capacity analysis was divided into high and low signal to noise ratio (SNR) regimes, where the dominant factors for the signal to interference plus noise ratio (SINR) were determined as interference and noise.A multicast beamforming approach was proposed in [9], where users located within the proximity of the BS can receive two different data streams simultaneously, while those away can receive only one data stream.In [7], a model was designed to transmit beamforming vectors to maximise the signal to leakage ratio (SLR), minimising transmission power which may cause interference to other users.The model is not restricted to a certain number of transmission antenna elements and can be extended to more general scenarios.An asymptotic deterministic SLNR optimisation approach for regularized zero-forcing (RZF) considering perfect CSI and antenna correlation was proposed [10].It was found that when the users are homogeneously distributed and the number of antenna elements is large enough ( → ∞), the SLNR is asymptotically equal to SINR.In [11], an efficient statistical Eigen-mode space division multiple access (SE-SDMA) scheme for downlink ergodic sum-rate analysis based on SLNR was proposed.To maximise the approximate ergodic sum-rate capacity, authors in [12] designed an optimal beamforming vector for each user based on a three-dimensional beamforming algorithm.The authors in [3] focused on imperfect CSI scenario and developed a robust SLNR approach for compensating performance degradation caused by random CSI errors.In [13], a precoding scheme for heterogeneous networks (Het-Nets) was proposed based on SLNR under imperfect CSI scenario.In this case, regularization parameter was used to "weight" the precoding information for other cells (BSs).
In this paper based on the SLNR concept and Ricean fading channel model with perfect CSI in MU-MIMO systems, the statistical Eigen-mode (SE) and zero-forcing (ZF) models are derived.In addition, the average achievable sumrate by each user and the overall average system capacity with different number of users are analyzed.Besides, the expression of dirty paper coding (DPC) method as a function of number of users () and number of antenna elements () is addressed.Moreover, the impact of  on the average achievable sum-rate is demonstrated.The Ricean fading channel with the line-of-sight (LOS) and scattering components as channel model are adopted.We find that the ergodic capacity of DPC is approximately free from channel matrix which mainly depends on the number of BS antenna elements and average transmission power.
Furthermore, by exploring the leakage signal criteria, we propose a new design solution for precoding based on the proposed power iteration technique.The technique finds the optimum weight vector which maximises the SLNR, that is, enhances the system capacity by obtaining the dominant eigenvector for minimising the CCI.Finally, we compare the proposed method with the conventional solutions such as SVD.The numerical results demonstrate that significant throughput performance can be achieved with the proposed technique.
A comparative study of the proposed and SVD methods regarding computation complexity and storage is also carried out.The simplicity of the proposed algorithm is observed in approximating only one eigenvalue of a matrix in a sequence which is considerably more efficient as the number of iterations increase.In comparison, the SVD has to calculate all eigenvectors and pick the one with the most significant eigenvalue (the maximum eigenvalue), which requires larger computation capability and storage capacity.
The contributions of this paper are in twofold: (1) The SE and ZF beam former models in Ricean fading channel are examined based on SLNR instead of the commonly used SINR and orthogonal condition of the beam weight vector in MU-MIMO systems.The linear achievable sum-rate capacity is investigated using these models to validate the proposed method.Moreover, in nonlinear achievable sum-rate capacity, the DPC technique regarding the impact of the number of users and antenna elements on the ergodic capacity is addressed.The results show that the number of antenna elements has a significant contribution to the ergodic capacity, and the DPC sum-rate capacity logarithmically increases with the number of users.
(2) A new optimum weight vector is developed based on the proposed power iteration method, which allows each user to maximise the SLNR and minimise the CCI from other users.The proposed method compared with the conventional SVD method shows that the proposed method can achieve better performance with relatively low cost.
This paper is organized as follows.In Section 2, the MU-MIMO downlink system is introduced.The beamforming and leakage signal approaches are presented in Section 3. In Section 4, the results and discussion of the study are presented.Finally, the conclusion is drawn in Section 5.The significant proofs of the algorithm are addressed in the appendix.
Leakage from U1 to others Mobile users (K users)

System Model
Consider the MU-MIMO system model shown in Figure 1.
The BS is equipped with  uniform linear transmission antennas that serve  users simultaneously.Assume that each user is equipped with a single antenna element [3].
Additionally, the number of the BS antenna elements is assumed to be larger than that of users ( ≥ ); and the power is equally allocated among all users [11,12].

Signal Model.
From the previous assumptions, the received signal   at th user is given by The expression in (1) can be rewritten in a more simplified form as where h  ∈ C 1× is the flat Ricean downlink channel vector between the BS and the th user.The weight vector w  ∈ C ×1 is unit normalised beamforming vector of user  and satisfies ‖w  ‖ = 1.Furthermore, the vector   is the transmitted data symbol of user  with [|  | 2 ] = 1.  → (0,   ) is unit normalised ( 2  = 1) complex additive white Gaussian noise (AWGN).Moreover,   is the transmission power of the th user under the constraint ∑  =1   ≤ .Since the total transmission power () is divided among all the users equally, the average power can be written as   = / for each user [3,11,12].Consequently, the signal magnitude of each user can be written as |h   w  | while the interference of the th user is determined by the sum of interference from all other users.Hence, the interference can be denoted as ∑  =1, ̸ = |h   w  | 2 .Thus the SINR of the th user can be expressed as in (3) [12,14]: (Note: for simplicity we drop the "subscript " from   in the rest of the paper).
The corresponding average achievable data rate for user  can be written as By substituting (3) into (4), the average achievable sum rate per user is given as in (5) [15,16].
Subsequently, the overall system ergodic sum-rate capacity can be written in bps/Hz as

Channel Model.
In this subsection, we describe the channel model as in [11,12].In the most recent studies on mmWave communication, the Rayleigh fading channel model is employed which faces difficulty in capturing the fading variation in the presence of the LOS component.Moreover, mmWave MU massive antenna system can be used mainly in 5G hotpot scenarios in the LOS case.To overcome the capturing problem in Rayleigh fading channel, the Ricean fading channel is applied in this paper.Considering the Ricean fading channel model, the channel vector has two components: LOS and Rayleigh random distribution components.Thus, the channel vector can be expressed as [11] where   is the power ratio between the LOS component and scattering components, also called Ricean -factor of the Wireless Communications and Mobile Computing th user.Moreover, the channel vector h , is the random component, and the entries are independent and identical distribution complex Gaussian random variables.Vector h  ( = 1, 2, . . ., ) is the deterministic component of the channel mean vector of the th user.Based on the assumption that the BS is equipped with  uniform linear array antennas, the deterministic channel vector for each user can be written as [17,18] where  is the space between two antenna elements,  is the wavelength of the signal, and   is the angle of departure (AoD) for the th user,   ∈ [−, ].

Beamforming and Leakage Signals
It is assumed that all users have a perfect instantaneous knowledge of their channel vectors while the BS knows all users channel vector including the channel mean vector.Additionally, assuming  users, it is difficult to maximise the average system capacity due to a couple of  calculations (capacity of each user) required.Therefore, it is hard to use SINR directly to obtain the optimum beamforming vector (w  ).To overcome this problem, the SLNR can be easily implemented to control the leakage signals from a specific user to other users as shown in Figure 1.The SLNR of the th user can be written as in ( 9) by assuming unit normalised Gaussian noise [4,16]: Equations ( 3) and ( 9) have different denominators to define SINR and SLNR.In (3), the user channel vector h  is used with other users' beamforming vectors w  |  =1 ,  ̸ =  to calculate SLNR of the th user.While in (9) the th user beamforming vector w  is used with other users' channel vectors h  |  =1 ,  ̸ =  to calculate SLNR of the th user.Besides, the lower bound (LB) on the average SLNR of the th user can be represented as {SLNR  } ≥ {SLNR  } LB demonstrated in [11,16]: In (10), we can have R  ≜ {h   h  }: where R  is a Hermitian matrix (channel correlation matrix) which can be constructed as in (12), where R  ≜ h   h  , and h  is calculated in (8).Hence, R  of the th user can be decomposed regarding the deterministic component of channel vector h  after normalisation by a factor √  [6].Therefore, we can use the normalised version h  / √  in the decomposition of the channel correlation matrix as in where Ũ  is the orthogonal subspace of h   / √ .Then we can obtain the optimum beamforming vector which maximises the lower bound of the SLNR as follows: On the other hand, to minimise the denominator of (9) we have where vector w opt  is orthogonal to (1/ √ )h   .Accordingly, from (12) we have The maximum value of the lower bound can be achieved if . By substituting (15) in (10), we have . (16) In (16), the lower bound is directly affected by the number of antenna elements , the Ricean factor , and the number of users .Assuming that the number of users is increased as  = 1 → , taking  = 0 and  → ∞, the lower bounds of SLNR are expressed approximately as (/(1 + ( − 1)))| =0 and  * | →∞ , respectively.Therefore, it is seen that the lower bound is independent of  when  = 0.

Analysis of Linear Achievable
Rates by SE and ZF.In this section, we evaluate the Ergodic sum-rate capacity of SE obtained by ( 13) and ( 14) and assume that ZF capacity has perfect CSI when calculating the mean gap loss between  SE and  ZF .Based on the orthogonal beamforming condition in ( 13) and ( 14), the achievable sum-rate capacity obtained by the th user is given by The achievable sum-rate capacity by ZF is expressed as where w ZF  is the unit normalised beamforming vector, which is selected as the th column of the normalised matrix to suppress the interference from all other users for |h  w ZF  | = 0, if  ̸ = .Accordingly, the mean gap loss Δ ZF−SE with perfect instantaneous CSI is given by From the orthogonal condition (1/ √ )h  h   = 0|  ̸ = which is addressed in (14), we neglect the interference from the other users to the th user regarding the signal component.In this case the sum-rate capacity  SE used in (14) When the number of users is large, the DPC achievable sumrate capacity is approximately expressed as [19,20] From (23), we can observe that the sum-rate capacity increases linearly with the number of antenna elements.However, the approximate achievable sum-rate capacity is a nonlinear function of .Therefore, the DPC is regarded as a nonlinear model.

Optimum Weight Vector Formulation.
In this subsection, we find the solution for optimum weight vector w opt , in which the maximum SLNR can be achieved to guarantee that all users have the ability to access the limited resources.For equal transmission of power to all users in a given case, we need to design w  ,  = 1, 2, . . .,  to maximise the SLNR for each user at the same time minimise interference.Hence, the optimum weight vector w opt  = arg max(SLR  ) can be written as follows: where h is congregate interfering channel transfer matrix (CICTM) of the th user, which is an extended channel matrix that excludes h  only for the th user.Then h  is written as where h ∈ C ×(−1) and h  ∈ C 1× .The optimum weight code vector can be derived from (28) as follows [21]: where I  is  ×  square identity matrix.
The optimum weight code vector in (30) is called SVD solution, which is determined from the matrix in (29).Sometimes, when calculating the maximum eigenvalues for the given matrix, the dominant eigenvector does not correspond to the optimum code weight vector.Hence, the eigenvalues are put in ascending or descending order; then an optimum vector corresponding to the most significant eigenvalue is easily selected but at the expense of extra overhead to the system [4,7].To enhance the system performance and minimise the cost, a new approximation method expected to improve the average user achievable sum-rate capacity is proposed.

Proposed Optimum Weight Vector Design.
A new method is developed based on power iteration method to find the optimum weight vector.The method requires only one weight vector during implementation.Equation ( 30) is modified and adopted as the main matrix in the proposed power approximation method as follows: 3.4.1.Power Approximation Method.The power approximation method generates a sequence of vectors A  x 0 , where x 0 is a nonzero initial selected vector.By normalising these sequence vectors under conditions stated in Section 3.4.3, the vector converges to the dominant eigenvector corresponding to the most significant eigenvalue.The normalisation is used to ensure that the most significant component of the given iteration is equal to one [22].The steps in the proposed power iteration method are given in Algorithm 1.
This algorithm is summarised in a sequence of iterations as follows: where   is the component (magnitude) of vector Ax −1 .
Before proceeding to the next iteration, it is necessary to scale down (normalised) the sequence vectors in each approximation so as to keep the largest component at unity.

The Speed of Convergence. Referring to (A.6) in
Appendix A, we observe that the coefficient of V  in the sequence X  goes to zero, which is proportional to ( j / 1 )  ; the convergence speed of sequence {X  } to V 1 is governed by ( 2 / 1 )  .Therefore, the rate of convergence and the convergence of the constants {  } to  1 are linear.For any linear convergent sequence {  }, we can use the Aitken Δ 2 technique which is used to make the linearly convergent sequences fast.The new convergence of the sequence {  } can be written as The convergence property of the algorithm is given by the following.
eorem 1. Assume that there is one and only one eigenvalue  1 of A, and  1 is semisimple; then either the initial vector x 0 which has no component in the invariant subspace associated with  1 or the sequence of vectors generated by the algorithm converges to the eigenvector associated with  1 and  and converges to  1 .The proof is shown in Appendix B. Definition 2. Eigenvalue  1 of matrix A is semisimple with the degree of  if it has a geometric multiplicity one and algebraic multiplicity .The vector A  x 0 is normalised by a specific scalar   to make the most significant component of the vector unity.The initial vector x 0 is decomposed as follows: where   ,  = 1, 2, . . ., , are the spectral projectors associated with the eigenvalues   ,  = 1, 2, . . ., .
From the formula A  =   (  Ι + D  ), the power can be written as A    =   (  Ι + D  )  , where Ι is the identity matrix and D  is the diagonal matrix.Consequently, we have Referring to the Definition 2, D 1 = 0 because  1 is a semisimple eigenvalue.Thus, we obtain The spectral radius of each operator (    + D  ) −1 < 1 and |   −1 1 | < 1 means that the th power will converge to zero.Theorem 1 is true when  1 x 0 = 0.However, if  1 x 0 ̸ = 0, x  converge to normalise  1 x 0 so that the most significant component is one.Meanwhile, the scalar   converges to the eigenvalue  1 , which is an immediate consequence in the form Thus, the sequence of vectors x  is proven to converge to the optimum weight vector which maximises the SLNR in (28).The optimum weight vector is equal to the vector x  : The result in (38), which is the main objective of this paper, maximises the th user SLNR when substituted in (28) to find SLNR  .

Cost Analysis.
In this subsection, the cost function is analyzed regarding storage space required for the proposed power iteration technique and SVD method.Referring to (30), the value (I  + h *  h ) −1 h *  h  results in a new square matrix with dimension ×, where  denotes the number of antenna elements.Therefore, finding w opt  for each user by using the solution of SVD in (30), we need to find at least two matrices,  ×  left eigenvectors matrix,  ×  right eigenvectors matrix, and  diagonal eigenvalues from the original  ×  matrix.On the other hand, the required storage memory for the th user is equal to the size of three matrices in addition to the original square matrix.Hence, the total storage becomes (3 2 + ) *  (full SVD), where  is the number of users.Then, the determination of the  optimum weight vector can be from the right eigenvector or left eigenvector based on the corresponding maximum eigenvalue.Furthermore, we can reduce the size by taking only one side eigenvectors matrix (left or right) to save  2 redundant storage size.Thus the reduced new storage capacity can be written as (2 2 + ) *  (reduced SVD).
Similarly, only the original square matrix with dimension  ×  and a vector with size  × 1 is required in the proposed power approximation method.In this case, there is no need for storing the previous iteration result as it can be overwritten up to the last iteration or reaches the determined tolerance.Thus, the required storage for the power iteration technique is ( 2 + ) *  which is much less than full and reduced SVDs.
Figure 2 shows the storage capacity required in the proposed power iteration technique and SVD method.From Figure 2, it is observed that the number of antenna elements has a significant effect on the cost function.However, the proposed method has much less storage space compared to SVD hence recommended for a massive MIMO with large antennas.

Simulation Results and Discussion
In this section, the numerical results are presented with the following assumptions: the equal power allocation strategy (  = /) and noise effect (  = ) for all users are the same [23].Additionally, the AoD for each user is in a horizontal direction with the uniform distribution in the range [−, ].Furthermore, for each user the Ricean -factor has uniform distribution in the range of [ min ,  max ], where  min = 10 dB and  max = 30 dB.Moreover, the number of users is increased up to  = 32.For SE and ZF analysis  = 16, 32, 64, and 128 are selected as the numbers of transmission antenna elements.The ergodic sum-rate capacities of SE and ZF are obtained in Figures 3 and 4. In Figure 3.It is observed that the average system capacity of SE is, respectively, directly and inversely proportional to the number of antenna elements and number of users.Moreover, when the number of users is small, the impact of  on the SE capacity is more significant and conversely when the numbers are large, the impact is less significant.The SE achieves average capacities approximately in the range of 1.6-6.2bps/Hz.At maximum number of users ( = 32) with  = 16, the corresponding SE sumrate capacity is less than 1.6 bps/Hz.This is because each user faces CCI interference which is considered as sum of leakage signals from all other cochannel users.Referring to (17), it is validated that as the number of users increases, the overall contributed interference to the th user also increases.Therefore, the increase in the denominator in (17) causes a reduction in the sum-rate capacity per user resulting in an average lower system capacity.
In Figure 4, it is shown that the ergodic sum-rate capacities obtained by ZF are directly proportional to the number of antenna elements.By referring to (18) the interference from the other users is perfectly suppressed by the orthogonal beamforming for user channel and weight vectors.These vectors are orthogonal to each other based on the condition |h  w ZF  | = 0 in which  ̸ = .Moreover, the number of antenna elements  has much contribution on the average system capacity with ZF than SE as comparatively shown in Figures 4 and 3.The maximum achievable capacities are 105 bps/Hz and 6.2 bps/Hz for ZF and SE, respectively, because of the interference effects on SE.
Figure 5 shows the upper bound of the mean-rate capacity gap loss ( ZF −  SE ) between ZF and SE.Since the number of users tends to be large, the difference in the mean rate capacity is most likely the ZF capacity, because the capacity obtained  by ZF is close to 15 times that obtained by SE, which precisely affects the mean-gap loss.Nevertheless, when the number of users increases, the difference tends to be very big.It is explicit that when the number of users is large under any number of antenna elements, the achievable sum-rate capacity obtained by SE almost tends to be zero.However, the ergodic sum-rate capacity achieved by ZF positively increases with increase in the number of antenna elements and tends to a fixed level with large number of users.Meanwhile, the achievable sumrate capacity is not profoundly affected by the number of users.
As shown in Figures 3-5, the simulation results depict that the ZF seems ideal due to a perfect cancellation of CCI.However, in practice SE is dominant and realisable because it is challenging to ensure that the interference is perfectly cancelled by the orthogonal condition between the channel vector and optimum weight vector.Thus, carefully focusing the beams alignment between the BS and the users a high system capacity can be achieved.
In Figure 6, the antenna elements  = 16, 32, 64, and 128 are considered with number of users up to 32 to satisfy  ≥  and show the effects of  <  ( = 16).The results show that the ergodic sum-rate capacity achieved by DPC is nonlinearly increased with the number of users.Moreover, if the number of users is fixed to 25, at the number of antenna elements 16, 32, 64, and 128, respectively, the corresponding ergodic sum-rate capacities are 80, 105, 175, and 250 bps/Hz, respectively.As the number of users is increased, the system capacity also increases.Thus, comparing the DPC with ZF and SE achieves much higher capacity with large number of antenna elements and users; hence it can provide multiuser diversity gain.
The results in Figures 7-9 are used to comparatively evaluate the proposed power iteration method in (38) and the SVD method in (30).In Figure 7, it is observed that the SLNR in the range −10.5 dB to 37 dB achieves the system capacities in the range 5 bps/Hz to 11 bps/Hz for the proposed method.On the other hand, the SVD method has the SLNR in the range −15 dB to 31 dB that achieves system capacities in the range 1 bps/Hz to 6 bps/Hz.It is clear that the proposed method has the average throughput nearly two times that of SVD.The main reason is due to the efficiency of the weight vector obtained in the proposed method that is capable of minimising the CCI much more than in the SVD.It is expected that the system capacity can be enhanced in case the BS has more antenna elements but at a higher simulation time.Additionally, the SVD method has a drawback in finding the exact eigenvector corresponding to the dominant eigenvalue in which the eigenvalues are not always in an orderly way (ascending or descending), which limits the SVD performance.Figures 8 and 9    SLNR of the proposed and SVD methods, respectively.It is observed that the proposed method achieves a much better SLNRs than the SVD method.

Conclusion
In this paper, based on MU-MIMO with massive antennas, the ergodic capacity for linear sum-rate analysis including SE and ZF beamforming based on SLNR technique is investigated.In addition, the nonlinear sum-rate analysis using DPC in Ricean fading channels based on SLNR was undertaken.A new method is proposed to find an optimum beam weight vector by exploring the power iteration method using eigenvector approximation.The number of antenna elements is selected up to 128 dramatically increasing the system capacity.By way of simulation, the most significant dominant eigenvector to maximise SLNR as well as minimise the CCI is obtained.By comparing with the SVD method, the proposed method achieves higher performance regarding mean achievable sum-rate capacity per user.The throughput of the proposed method is in the range of 5 bps/Hz-11 bps/Hz, while that of SVD in the range of 1 bps/Hz-6 bps/Hz.Therefore, the proposed method can provide significant system capacity enhancement.In future work, the proposed power iteration technique is recommended for a 5G MU massive antenna system.

Figure 1 :
Figure 1: Schematic diagram for leakage signals in MU-MIMO system model from the first user to all other users.

Figure 2 :
Figure 2: Storage capacity required by the proposed technique and the SVD (full SVD and reduced SVD) method regarding the number of antenna elements.

Figure 3 :
Figure 3: Sum-rate capacities for SE on the number of antenna elements.

Figure 4 :
Figure 4: Achievable sum-rate capacities for ZF beamforming on the number of antenna elements.

Figure 5 :
Figure 5: Mean gap loss capacity ( ZF −  SE ) which is the difference between ZF and SE capacities.

Figure 6 :
Figure 6: Achievable sum-rate capacities for DPC on the number of antenna elements.
are the Cumulative Distribution Function (CDF) and Probability Density Function (PDF) of Wireless Communications and Mobile Computing

Figure 7 :
Figure 7: Sum-rate capacities for the proposed power iteration technique and SVD method.

Figure 8 :
Figure 8: CDFs of the SLNRs for the proposed power iteration and SVD methods.

Figure 9 :
Figure 9: PDFs of the SLNRs for the proposed power iteration and SVD methods.
[19]ges to the loose bound sum-rate capacity of the th user expressed as[11]Nonlinear Achievable Rate Analysis by DPC.The BS can serve all users simultaneously and achieve maximum system capacity as much as possible as the CSI is known at the BS.In MIMO systems, the system capacity obtained by downlink strategy is called DPC.Practically, it is challenging to implement DPC because the encoding and decoding have high computations which is ineffective for large number of users.The achievable sum-rate capacity of the DPC can be written as[19] = Ax −1 , and x  = y  /  ; (normalisation) "the value of x  is the scaled version of y  " if |  −  −1 | ≤ Tol, then  =  + 1; arg max (     h  w      2 ∑  =1, ̸ =      h  w       2 ) (24) SLR  =     h  w      2 ∑  =1, ̸ =      h  w       2 .(25) 3.4.2.Convergence of the Algorithm.The matrix A is  ×  square matrix with  eigenvalues [ 1 ,  2 ,  3 , . . .,   ], which are in descending order as | 1