Abstract
The present paper investigates the tradeoffs between the energy efficiency (EE) and spectral efficiency (SE) in the fullduplex (FD) multiuser multiinput multioutput (MUMIMO) cloud radio access networks (CRANs) with simultaneous wireless information and power transfer (SWIPT). In the considered network, the central unit (CU) intends to concurrently not only transfer both energy and information toward downlink (DL) users using power splitting structures but also receive signals from uplink (UL) users. This communication is executed via FD radio units (RUs) which are distributed nearby users and connected to the CU through limited capacity fronthaul (FH) links. In order to unveil interesting tradeoffs between the EE and SE metrics, we first introduce three conventional singleobjective optimization problems (SOOPs) including (i) system sum rate maximization, (ii) total power minimization, and (iii) fractional energy efficiency maximization. Then, by making use of the multiobjective optimization (MOO) framework, the MOO problem (MOOP) with the objective vector of the achievable rate and power consumption is addressed. All considered problems are nonconvex with respect to designing variables comprising precoding matrices, compression matrices, and DL power splitting factors; thus, it is extremely intractable to solve these problems directly. To overcome these issues, we develop iterative algorithms by utilizing the sequential convex approximation (SCA) approach for the first two SOO problems and the SCAbased Dinkelbach method for the fractional EE problem. Regarding the MOOP, we first rewrite it as an SOOP by applying the modified weighted Tchebycheff method and, then, propose the iterative algorithmbased SCA to find its optimal Pareto set. Various numerical simulations are conducted to study the system performance and appealing EESE tradeoffs in the considered system.
1. Introduction
In the context of explosive growth in mobile data traffic and in wireless connectivities, the cloud radio access networks (CRANs) have been suggested as a modern architecture for the next generation of wireless communication systems. It is promising to overcome the strict requirements for concurrently enhancing the 1000fold data rate as well as reducing by half the power consumption [1–3]. Specifically, a typical CRAN composes of a central unit (CU), a set of distributed radio units (RUs), and terminal devices [1]. The CU is responsible for cooperatively processing baseband signals of all terminals and, then, communicates them to RUs through FH connections. These RUs which are inherently distributed nearby user equipment (UE) are responsible for transmitting or/and receiving data messages toward or/and from terminal devices. Based on this structure, the baseband signals are centrally processed at the CU, and thus, CRANs provide not only a cooperative data transmission but also a better interference management among RUs as compared with the conventional ones, for example, multicell models [4, 5] and coordinated multipoint (CoMP) schemes [6, 7]. In addition, another remarkable advantage of CRANs is the replacement of the highcost highconsuming power base stations by the lowcost lowconsuming power RUs that can significantly reduce the system power consumption.
On the other hand, with potentially increasing the degrees of freedom in the design, the multiple antenna transmission techniques can thus support to convey more data streams over a same radio resource and promisingly boost the system throughput. Hence, the multiantenna techniques have been integrated and investigated in the context of CRANs. The CRANs implemented with multipleantenna RUs and singleantenna terminals are named as MISO CRANs and investigated in [8–13]. Meanwhile, the ones deployed with multipleantenna RUs and multipleantenna UEs were referred to as MIMO CRANs and studied in [14–17]. Regarding the MISO CRAN deployment, [8] proposed the joint design of RU activation and beamformers to minimize power consumption in downlink (DL) MISO CRAN systems. In [9], the authors first proposed a spatialcompressionandforward strategy for uplink (UL) MISO CRAN systems and, then, addressed the beamforming and compressing design to maximize the minimum user's quality of service (QoS). Alternatively, the work in [10] developed a nonorthogonal multiple access (NOMA) transmission in DL MISO CRANs and, then, proposed a solution to handle the maximin individual user rate subjected to transmit power budget and limited FH capacity constraints. In addition, a wirelesspowered MISO CRAN system was developed in [11] where the terminal devices conduct both information decoding (ID) and energy harvesting (EH) in the first IDEH DL phase and each user then transmits its own data in the UL data transmission phase using harvested energy. The UL rate maximization problem is tackled under minimum DL rate requirements, power budget, and FH capacity constraints. Several novel techniques including hybrid digitalanalog beamforming and rate splitting multiple access were addressed in DL MISO CRANs in [12, 13], respectively. As the result, the system spectral efficiency (SE) is thus considerably improved. On the other hand, by considering the MIMO CRAN scheme, the joint design of precoding and compressing matrices was investigated for the UL MIMO CRAN in [14, 15] and for the DL MIMO CRAN in [16], respectively. Additionally, the authors in [17] took multicluster interference into consideration and developed an effective algorithm to handle them in the multicluster DL MIMO CRAN.
It is worth noting that the aforementioned works separately studied the UL CRAN channels or DL CRAN ones and can be referred to as halfduplex (HD) transmission that leads to inefficient spectrum usage. Recently, thanks to the advanced analog and digital selfinterference (SI) cancellation strategies, the fullduplex (FD) radio is able to concurrently transmit and receive signals over the same radio resources. Then, the FD techniques potentially provides the SE approximately twice as much as that of the conventional HD counterparts [18, 19]. The integration between the FD radio approach and CRAN scheme has recently drawn significant attention [20–26]. Reference [20] investigated a classical Wyner model based on FDaided CRANs and illustrated that the achievable rate of the FD CRAN system is significantly improved in comparison to both HD CRANs and singlecell processing ones. Reference [21] developed a twostage algorithm to jointly design RU beamformers, user power allocation, and user selection for the FD MISO CRAN. Furthermore, the system data rate for the multiuser scheme was analysed and optimal power allocation for the singleuser case was investigated in the FD CRANs [22]. Additionally, by considering a hybrid RU set which comprises both multiantenna and singleantenna RUs, [23] aimed for joint resource allocation and interference cancellation in order to maximize the DL rate subjected to power and UL rate QoS constraints in the FD CRANs. By exploiting the difference of the convex (DC) form in the system data rate and FH capacities, the authors in [24] proposed an iterative algorithm based on the majorizationminimization method for the FD MIMO CRANs which jointly optimizes precoding and compressing matrices under transmit power and FH rate constraints. Alternatively, the work in [25] developed a lowcomplexity design for FD MIMO CRANs by making the use of the connection between weighted minimum mean squared error (WMMSE) problem and weighted sum rate (WSR) problem. The precoders and compressing matrices are then alternatively optimized by solving the convex problem in each iteration. Most recently, in [26], the simultaneous wireless information and power transfer (SWIPT) protocol was investigated on the FD MISO CRANs. The joint design of DL beamforming, EHpower splitting ratio, UL power allocation, and RU postcoding was addressed. Reference [27] investigated the energythroughput tradeoff by designing the beamforming for SWIPT CRANs. Alternatively, the authors in [28] developed the hybrid timepower splitting protocol for downlink user information and uplink IDEH sensor in FDaided massive MIMO wireless systems.
However, the aforementioned works were concentrated only on designing either the sum rate maximization (referred to as SE maximization) or the power consumption minimization in the considered CRAN. This might lead to ineffective and negative impacts on the system EE performance [29]. The EE issue was, therefore, taken into account in the transceiver design for CRAN networks [30–36]. Specifically, reference [30] investigated the resource allocation in order to minimize the entire system power consumption for both data sharing and data compressingbased DL CRAN models. The authors in [31] first generalized these models to the multihop FH links MIMO CRANs and, then, proposed a successive convex quadratic programming (SCQP) algorithm to jointly design precoders, userRU association, and RU activation to optimize the system fractional EE defined as a ratio of the system sum rate to power consumption. In order to further reduce power consumed by ineffective RUs, the issues of turning RU on/off and userRU association were studied in various recent works. Particularly the joint design of precoding matrices, turning RUs on/off was investigated for both SE and EE optimization problems in DL MIMO CRANbased multihop FH link networks [32]. In [33], the EE maximization problem constrained with FH capacity limitation, transmit power budget, and user’s QoS requirements was tackled for a DL MISO CRAN. Alternatively, the EE beamforming designbased RU selecting and RUUE associating designs for DL MISO CRAN were proposed either with virtual computing resource models in [34] or with FH ratedependent power consumption models in [35]. To take advantage of FD techniques, reference [36] proposed an efficient algorithm for optimizing the system EE in FD MISO CRANbased data sharing systems. However, the FH capacity constraints and FH power consumption model were not taken into consideration.
We would like to emphasize that, in the aforementioned works, only one specific system performance metric is considered as a singleobjective function in design problems while the others are taken into account in the constraints. This approach can be referred to as singleobjective optimization (SOO) design and can lead to an unfairness between the conflicting system measures, such as systems SE, EE and power consumption. More recently, the multiobjective optimization (MOO) framework has been introduced as an effective mathematical tool to overcome this issue as well as to provide a balance among the system performance metrics on designing the wireless communication systems [37, 38]. As a result, the MOO framework was recently exploited in many studies to address various tradeoffs between two or more system metrics [39–48]. Particularly, the harvested energy and information transfer tradeoffs in SWIPTenabled wireless system were studied either in cognitive FD MUMIMO systems [39] or in MIMO broadcast channels [40, 41]. In [42], the authors studied the optimal Pareto set for DL and UL transmit power fairness in the FD MISO networks. Furthermore, the balance between transmit power and achievable rate was investigated in [43, 44] for secure DL MUMISO networks and DL MISO CRANs, respectively. In addition, the EESE tradeoffs for massive MIMO systems were addressed in [45, 46]. On the other hand, the works in [47, 48] have developed algorithms using MOO frameworks to study the balance among three essential and conflicting metrics. Particularly, the balance among three inconsistent design objectives of the MSE minimization, harvested energy maximization, and physical security optimization was revealed for SWIPTaided MIMO interference channels in [47]. Furthermore, the work in [48] addressed the fairness among total transmit power minimization, EH efficiency maximization, and interference leakagetotransmit power ratio minimization for SWIPTbased MISO cognitive radio networks.
Motivated by the above viewpoints, this paper studies EESE tradeoffs in an FD MUMIMO CRANenabling SWIPT network (shown in Figure 1). The considered CRAN model includes that a set of FD RU stations, which is connected to the CU through the limited capacity FH links, concurrently communicates to the sets of EHenabled DL users (EHDLUs) and UL users (ULUs). Specifically, regarding DL channels, the baseband signals associating with DLUs are first centrally processed at the CU and, then, conveyed to the RUs via FH links. After radio frequency (RF) steps, the distributed RUs transmit these signals toward DLUs through wireless access links. At the DLU, the received signal is divided into two proportions using power splitting (PS) receiver structures, which are then performed for ID and EH purposes. The total amount of harvested energy can be stored at the DLUs to prolong their lifetime or to utilize for signal processing operations as discussed in previous works [26, 39, 47, 48]. In the UL channels, the individualreceived baseband signals at RUs are sent to the CU via FH links for central processing. It is noted that, in order to guarantee the rate over finite capacity FH links, the baseband signals are quantized and compressed by using FH compression strategies [14–16]. Hence, we aim to jointly design precoding matrices, FH compressing covariance matrices, and power splitting parameters. Moreover, to investigate the tradeoffs between two conflicting EE and SE metrics, we introduce three SOO problems (SOOPs) and one MOO problem (MOOP); they are all under constraints of transmit power budget, minimum data rate, maximum FH link capacity, and minimum amount of harvested energy. Specifically, the SOOPs include (i) sum rate maximization, (ii) power consumption minimization, and (iii) fractional energy efficiency maximization whereas the MOOP is the simultaneous sum rate and power consumption optimization one. All interesting problems are the challenging nonconvex optimization ones and hard to be directly handled. This is due not only to the nonconcavity in data rate functions and nonconvexity in FH link rates caused by the coupled of UL and DL precoder variables but also to the inverse convexity in EH constraints. To overcome these issues, we exploit the efficient inner convex approximation and SCQP approaches to develop iterative algorithms with low computational complexity.
The main contributions of this work are summarized as twofold as follows: (1)This paper attempts to comprehensively study the SEEE tradeoff in the FD MIMO CRANs with SWIPT. Compared to previous related works in the literatures [24–26, 36], we highlight the novelty of the present work as follows: (i)Although the works in [24, 25] have addressed the transceiver design for FD CRANs in SE perspectives, the present paper first extends the system model for enabling SWIPT at DLUs and, then, studies the additional system EE metric using two approaches, namely, the maximization of the ratio between the sum rate over the power consumption based on the SOO programming and simultaneous sum rate maximization and power consumption minimization based on MOO programming.(ii)On the other hand, it is worth noting that the work [26] has also attempted to minimize power beamforming under the QoS of the data rate and harvested energy in SWIPTenabled FD MISO CRANs while the authors in [36] have also proposed an algorithm to optimize the overall system EE in FD MISO CRANs. However, both works assumed singleantenna users and unlimitedcapacity FH links, therefore, those approaches cannot be directly applied for FD MIMO CRANs with SWIPT considered in this paper. Thus, the present work is distinguished itself and can be referred to as the more comprehensive investigation(2)We also develop a computationally efficient iterative algorithm based on the SCQP approach, in which the logdeterminant expressions in both objective function and FH rate constraints are recast by the quadratic convex ones within each iteration. As the result, the computational cost, which is essentially caused by logdeterminant expressions, is significantly reduced as compared with the methods in [24, 25].
The rest of this paper is organized as follows. Section 2 introduces the system model and transmission signaling and presents the problem formulations for SOOPs, which include sum rate maximization, power consumption minimization, and energy efficiency maximization, as well as introduces an EESE tradeoff MOOP. Section 3 develops algorithms for considered SOOPs and MOOP. Finally, the simulation results and conclusions are provided in Section 4 and Section 5, respectively.
Notations. Bold lower and uppercased letters represent vectors and matrices, respectively. , , , and are the transpose, Hermitian transpose, trace, and determinant of , respectively. The notation (, respectively) means that is a positive semidefinite (positive definite, respectively) matrix. represents the identical matrix. We denote and . Operator is the real part of the complex number and denotes the expectation operator.
2. System Model
This paper considers an FD MIMO CRANenabled SWIPT as shown in Figure 1, in which two sets of UL and DL users communicate with a CU through a pool of RUs. It is assumed that the limitedcapacity FH links are used in CURU connections. Hence, to meet the bottleneck rate constraint, the data is compressed before being transferred over FH links. For notation simplicity, the sets of RUs, DLUs, and ULUs are denoted as , , and , respectively. The number of transmit and receive antennas at RU are and , respectively. Thus, the numbers of all transmit and all receive antennas at RUs are and , respectively. In addition, the numbers of transmit antennas of ULU and receive antennas of DLU are denoted as and , respectively. Moreover, the channel state information (CSI) of all links is assumed to be fully known at the CU by exploiting the joint UL and DL channel estimation algorithm [49]. It should be noticed that such an assumption of the perfect CSI is widely adopted in the recent related works [9–17, 23–26]. In practice, the CSI can be achieved through the channel estimation, feedback, or reverse links, and thus, the imperfect CSI is available at the CU. In such cases, the results in this work can be referred to as the upperbound performances.
At the same time and over the same frequency resource, the CU intends to transmit the message for DLU in the DL channels while the ULU also wishes to convey its own message to the CU, with and representing the numbers of data streams of DLU and ULU , respectively. All data symbols are assumed with zero mean and unit powerindependent distribution, i.e., with and with an appropriate dimension for all . In addition, to take advantages of multiantenna transmission in spatial diversity enhancement, the DL and UL signals are first beamformed by multiplying with its corresponding precoding matrices, namely, for DLU and for ULU . Thus, in the UL channels, the transmit signal at ULU is given as . Then, the power constraints at ULUs, i.e., , are determined by
Concerning the DL channels, it is emphasized that the precoding matrix for DLU is constructed at the CU in the form [16, 17], in which is the precoder matrix for DLU at RU . Hence, the stacked DL data vector at RUs is given by and the DL data vector corresponding to RU is given as with
Prior to being conveyed throughout the limited capacity FH link, the DL signal for RU is quantized and compressed. Adopting a Gaussian test channel, the compressed DL signal for RU is given as where the quantization noise of DL FH link , , is independent of precoded data stream and obeys as in [16, 17, 24, 25]. This process is to guarantee the rate constraint over the FH link. Hence, the power constraint at RUs, i.e., with , is equivalently expressed as
2.1. DL Channels with Information and Energy Transferring
From (3), the noiseless FH rate of the link from CU to RU is determined as [16, 17]
In addition, the received signal vector at DLU , , is given as where and denote the channel matrix from RUs to DLU with being the complex channel matrix between RU and DLU while represents the CCI channel coefficient to DLU from ULU and represents the additive Gaussian noise vector at DLU with zero mean and covariance matrix .In order to realize SWIPT via received RF signals, this paper adopts power splitting receiver structures at DLUs. Specially, we exploit perantenna power splitting (PAPS) structures in EHDLUs. Denoting as the signal proportion factor for ID at antenna of DLU with , the PAPS coefficients of DLU are represented by a diagonal matrix which satisfies and . Then, the received signal for ID is where represents the additive circuitry Gaussian noise vector at DLU with zero mean and covariance matrix . Accordingly, the residual amount of signals is used for scavenging energy. Therefore, the signal for the EH purpose is given as
From (7), the DL rate at DLU can be expressed as the function of precoders, compression matrices, and PAPS matrices and given by where denotes the covariance matrix of the desired downlink signal and which is the covariance matrix of multiuser, cochannel, and compression interferences plus additive noise at DLU , where we have defined .Moreover, upon (8), the amount of harvested energy at DLU is expressed as follows where the DL signal covariance mapping matrix . In this paper, we consider a practical nonlinear model of energy harvester which is mathematically formulated as [50] where with . Meanwhile, denotes the maximum power harvested at DLU at the saturation point of EH circuit and and are the nonlinear EH parameters by hardware limitations as illustrated in [50]. This quantity of energy, , can be stored either in a battery or in a super capacitor and later exploited to prolong its own operation such as processing/decoding signals or transmitting control signals to the RUs [26, 39, 47, 48, 50, 51].
2.2. UL Channels with Information Transferring
As the transmission strategy discussed in the early subsection, the received baseband signal at RU in UL channels is given as where is the channel gain between ULU and RU whereas with is the interference link from RU transmission to RU reception. Specifically, if , interference is referred to as SI at RU ; otherwise, it is an interRU interference link between receive antennas at RU and transmit antennas at RU . The additive Gaussian noise at RU is denoted as with zero mean and covariance matrix .Then, by adopting a Gaussian test channel compression model, the compressed version of the UL signal for the FH link is where is the quantization noise which is randomized and independent from the received data stream [14, 15, 24, 25]. The rate corresponding to UL FH link is thus expressed as where
Define as the channel link from ULU to RUs. represents the aggregative link from transmit to receive antennas between RUs. In addition, , , and are the stacked versions of UL, DL quantization noise, and additive noise at RUs. Then, the stacked UL baseband signal vector received at all RUs, i.e., , can be expressed as
Similar to [24, 52], the optimal minimum mean square error and successive interference cancellation (MMSESIC) technique is used at the CU to decode the UL messages. Based on the MMSESIC decoding strategy, the user’s message with the strongest channel condition is retrieved first whereas that with the worst channel link is decoded last. After detecting any user data, the CU then removes it from the received signal for interference mitigation. As the result, the later decoding ULU suppresses from the multipleaccess interferences of all earlier decoding ULUs. Therefore, the data rate of ULU is given as [14, 24] where . Then, the total achievable rate of UL channels, i.e., , can be formulated as where is the received covariance matrix of the desired UL signal at RUs. Meanwhile, the covariance matrix of SI and compression noise plus additive noise at RUs is with being the diagonally stacking matrix of covariance compression noise of UL FH links.
2.3. Power Consumption Model
The whole transmit power consumed by RUs and ULUs is given by where and are the transmitted power at ULU and RU and given in (1) and (4), respectively. In addition, are the power amplifier efficiency factors of RU and ULU , respectively. Otherwise, the total circuitry power is given as where and are the circuitry power consumed at RU and ULU , respectively. Additionally, the power consumption of the FH links is defined as the ratedependent model; hence, the total power consumption over all DL and UL FH links is given as [30, 31] where and are constant scaling factors. and are FH rates given by (5) and (14), respectively. Meanwhile, and are the maximum rate to be able to convey through the DL and UL FH link , respectively. Nevertheless, and are the whole power consumption values when the FH link is active at the maximum rate [30, 31].
2.4. Formula of Considered Optimization Problems
Firstly, we introduce some system requirements including minimum harvested energy, minimum user’s data rate, and limited FH capacity requirements. In order to ensure the operating necessity of the energy harvester, the total amount of harvested energy at DLUs must satisfy the following condition where is the required harvested energy at DLU . The value of is chosen properly providing that the harvested energy is useful and can activate the EH circuitry [50]. Hence, this value should be set up to be larger than or at least equal to the minimum sensitive level of its energy harvester [26, 39, 48].
On the other hand, to meet the maximum limited FH capacity, the DL and UL FH rates are thus constrained as [24, 25]
Finally, the QoS management is predefined as the following achievable data rate constraints: where and denote the minimum QoS requirement of DL and UL data rates, respectively.
For notation simplicity, we define the compact set as
In this paper, we consider three SOOPs and an MOOP which are defined below.
2.4.1. Problem SOOP1
The sum rate maximization problem is written as a minimization one where is the overall UL and DL rate.
2.4.2. Problem SOOP2
The power consumption minimization problem can be expressed as where is the whole power consumed for data transmission, FH links, and circuitry.
2.4.3. Problem SOOP3
The conventional EE maximization problem is where is the whole system EE.
Note that (1) and (4) present the transmit power constraints at ULUs and RUs, respectively. The nonconvex constraints in (23) guarantee the minimum harvested energy at each EHDLU while the nonconvex constraints of DL and UL FH rate are given in (24) and (25), respectively. As a result, the three above problems are highly nonconvex due not only to nonconvexities in objective functions but also to nonconvex constraints. Thus, these problems are all extremely intractable to directly find their own optimum points. In the following sections, we develop the iterative algorithms using sequential convex approximation (SCA) methods to solve these nonconvex problems.
2.4.4. Problem MOOP
The energyspectral efficiency optimization using the MOO framework is given as [37, 38] in which the objective function is a tworow vector whose elements are the SOOP1 and SOOP2 objective functions. Once again, it is very challenging to jointly find the optimal solution to these nonconvex multicriteria problems (32a) and (32b). Therefore, in the next section, we introduce an approach based on the modified weighted Tchebycheff and SCA methods to find the optimal Pareto set.
3. Proposed Iterative Algorithm
This section represents an approach to develop algorithms for considered nonconvex problems using SCA methods. Firstly, we provide three inner convex approximations of the data rate, FH rate, and harvested energy.
3.1. Convex Inner Approximation Analysis
In this subsection, we provide inner approximations in order to convexify the nonconvex problems.
3.1.1. Concave Quadratic Lower Approximation for the DL and UL Rate
Following Theorem 1 of [52], [53], and Appendix A that the DLU rate in (9) is lower bounded by the following concave quadratic expression, where and
Moreover, the concave quadratic lower bound of UL rate (18) is given as where , , , and
Accordingly, the nonconvex QoS rate constraints (26) and (27) are innerly approximated, respectively, as the following convex constraints
3.1.2. Convex Inner Approximation for DL and UL FH Rate Constraint
The convex upper bound of the DL FH rate is given by where with
In addition, the UL FH rate is convex upper bounded by where , with being found by substituting , , and into (15) instead. The complete proof is given in Appendix B.
Finally, the inner convex approximation of nonconvex FH rate constraints according to the DL and UL in (24) and (25), respectively, is given as
3.1.3. Convex Inner Approximation for DLU EH Constraints
Firstly, the compact constraint in (23) is equivalently rewritten after some simple mathematical manipulations as where . Next, by introducing additional diagonal matrix variables with , which satisfies the constraints the harvested energy at DLU is lower, approximated as (see Appendix C for detail) with
Consequently, the nonconvex constraint (43) is recast as the following convex inner approximation below
3.2. Iterative Algorithm for Considered Optimization Problems
Upon the above analyses, the nonconvex SOOP1 in (29a) is innerly approximated at iteration as the following convex one where and with .Additionally, the power consumption minimization problems (30a) and (30b) can be expressed as convex one at iteration where and given by where and are given in (39) and (41), respectively.
These two recast problems are now convex ones and can be efficiently tackled with existing convex solvers, such as CVX packages [54]. The whole procedure is summarized in Algorithm 1.

On the other hand, in order to find the Pareto set, the MOO problems (32a) and (32b) are first transformed to the SOO one by utilizing the modified weighted Tchebycheff method as [37] where and are the optimums of SOOP1 and SOOP2, respectively. The weighted factor with represents the tradeoff between two considered metrics. In addition, the small constant should be chosen to obtain the Pareto optimality [19, 37].
It is important to notice that problems (53a) and (53b) are still challenging to be directly solved due to the nonsmooth nonconvex property of the objective function and nonconvex constraints. Therefore, we exploit the analyses in the previous subsections and SCA method to develop an algorithm to iteratively handle it. Following inner convex approximations in the previous subsection, problems (53a) and (53b) are approximated at iteration as For given and , problems (54a) and (54b) are a convex one and can be handled with existing convex solvers, e.g., CVX packages [54]. The entire set of optimal Pareto points is achieved by varying value under the condition of . The total procedure to find Pareto points is summarized in Algorithm 2.

3.3. DinkelbachBased Iterative Algorithm for the Conventional Fractional EE Maximization Problem
Next, we turn our attention to the traditional fractional EE optimization problem. With the bound expressions derived in the previous subsection, the EE maximization problem can be recast at iteration as where for all with and being already defined below as (50a), (50b), (51a), and (51b), respectively. Problems (55a) and (55b) are now a concaveconvex fractional program with convex constraints. In specific, the objective is a pseudoconcave function that inspires us to exploit the nonlinear fractional programming approach [55]. Based on the Dinkelbach algorithm, we first reformulate (55a) and (55b) into a parameterized polynomial subtractive form with respect to as follows
Following [56], it is emphasized on parametrized problem (56) that is a monotonically decreasing and continuous function with respect to . This implies that is fulfilled with a unique solution which is given as On the other hand, the optimal solution of (56) at is also an optimal solution of (55a) and (55b) and is the optimal objective value. The detailed explanation can be refereed in [56]. It implies that one can first fix to find the optimal solution of convex problems (55a) and (55b). Afterward, we tackle identity to get the unique . Finally, the whole procedure to achieve the optimal solution is summarized in Algorithm 3.

3.4. Initiating a Feasible Point
At the beginning of the proposed algorithms, it is hard to find one feasible set. Thus, we consider the following maxmin problem which helps to find an initial feasible set. The initial feasible point is generated by solving the following maxmin problem The solution can be iteratively found by solving a inner convex approximation problem at iteration as where , , and are given in (33), (35), and (45), respectively. The algorithm terminates as the objective value is larger than or equal to .
4. Numerical Results
In this section, various numerical simulation results are conducted to illustrate EESE tradeoffs for the interested FD MUMIMO CRANs with SWIPT.
4.1. Simulation Parameters
We adopt the CRAN model in [24, 25], in which a square cell with a length of 100 meters is divided into four equally small square cells with a length of 50 meters. One RU is located at the center of each small cell and a pair of a DLU and a ULU is in each cell. Specifically, the DLU is randomly distributed within a circle with a radius of meters from RU while the ULU is randomized within its square cell area. In addition, the minimum distance between RU and any user is set as 5 meters. This cell setup implies that 4 DLUs, 4 ULUs, and 4 RUs are simulated. Furthermore, the transmit and receive antennas at each user is set as antennas for all and while each RU possesses antennas for all . The numbers of data streams are and .For the sake of convenience, an identical power budget is set for all ULUs and all RUs, i.e., , and , , respectively. In addition, the power budgets are set as dBm and dBm [24, 25]. The power amplifier factors are and (W) and with dBW for all [31]. The noise floor is set as dBW [51]. Regarding energy harvesters, the nonlinear EH parameters are set (mW) and for all [57]. The required harvested energy is for all and we also set (dBm). In addition, we also set the required data rate of DL and UL channels as for all and . The maximum capable FH capacities and FH power scale factor are set as and , respectively, for all . The required minimum rate is set to (bps/Hz). Without other statements, we also set (W) and (Mbps) with over the channel bandwidth of (MHz), which implies (W/(bps/Hz)) [30–32].
The channel matrix from node to node with distance meters are simulated as with being the path loss exponent. Meanwhile, the SI channel matrix at RU is conducted as where (in decibel) denotes the SI suppression level. We also set and dB [52]. The nominal smallscale fading matrix follows the Rician fading model as , where dB is the Rician factor, the lineofsight (LoS) term is a deterministic matrix, and each element of is an independent and identically distributed complex Gaussian variable with zero mean and unit variance, i.e., . In this paper, we utilize the farfield uniform linear antenna array model [51] , where and refer to the receiving and transmitting antennas, and with antenna spacing