# Distributed Optimization for Coordinated Beamforming in Multi-Cell Multigroup Multicast Systems: Power Minimization and SINR Balancing

## Abstract

This paper considers coordinated multicast beamforming in a multi-cell multigroup multiple-input single-output system. Each base station (BS) serves multiple groups of users by forming a single beam with common information per group. We propose centralized and distributed beamforming algorithms for two different optimization targets. The first objective is to minimize the total transmission power of all the BSs while guaranteeing the user-specific minimum quality-of-service targets. The semidefinite relaxation (SDR) method is used to approximate the non-convex multicast problem as a semidefinite program (SDP), which is solvable via centralized processing. Subsequently, two alternative distributed methods are proposed. The first approach turns the SDP into a two-level optimization via primal decomposition. At the higher level, inter-cell interference powers are optimized for fixed beamformers while the lower level locally optimizes the beamformers by minimizing BS-specific transmit powers for the given inter-cell interference constraints. The second distributed solution is enabled via an alternating direction method of multipliers, where the inter-cell interference optimization is divided into a local and a global optimization by forcing the equality via consistency constraints. We further propose a centralized and a simple distributed beamforming design for the signal-to-interference-plus-noise ratio (SINR) balancing problem in which the minimum SINR among the users is maximized with given per-BS power constraints. This problem is solved via the bisection method as a series of SDP feasibility problems. The simulation results show the superiority of the proposed coordinated beamforming algorithms over traditional non-coordinated transmission schemes, and illustrate the fast convergence of the distributed methods.

## 1Introduction

The performance of a communication system can be improved with a signal processing technique called transmit beamforming (or equivalently precoding). Given a certain quality of channel state information, properly designed advanced multi-antenna beamforming techniques have a potential of significantly improving spectral efficiency. However, inter-cell interference (ICI) may be a limiting performance factor without proper interference coordination between neighboring cells. To address this issue, coordinated beamforming has been adopted in the literature as a powerful interference management method. The idea of the approach is to involve the inter-cell interference in the beamforming optimization design, which improves the system performance especially at cell-edge areas [1]. In coordinated beamforming, each data stream is linearly precoded in the spatial domain and transmitted from a single base station (BS). More specifically, the aim is to achieve a network performance target with some user-specific quality of service constraints, by jointly designing precoded data transmissions among BSs. Coordinated beamforming schemes rely on the availability of channel state information (CSI) at the BSs. Depending on the implementation, the total amount of required CSI is different. In general, a centralized algorithm refers to a case where global CSI is required, which essentially means the knowledge of the channels between all the BSs and all the users in the system. On the other hand, in the distributed methods, each BS has to know only its channels towards all the users in the coordinating network to perform the optimization. In the literature, this is usually called local CSI. Throughout this paper, the acquired global or local CSI is assumed to be perfect. Generally, decentralized schemes may be easier to implement in practice since they reduce the signaling overhead and require lower computational requirements per processing unit. In the recent wireless communications research, coordinated beamforming has received an extensive attention for various system design objectives, such as sum power minimization [2], minimum SINR maximization (also called SINR balancing problem in the literature) [3], sum rate maximization [4] and energy efficiency maximization [5]. The sum power minimization problem [2] is to minimize the sum power of all the BSs with certain user-specific SINR targets. This system design objective is of practical interest for wireless applications which have stringent data rate and delay constraints. For the sum power minimization problem, the works of [2] and [2] have proposed centralized and distributed beamforming designs, respectively. The dual and primal decomposition was applied in [11] and [12], respectively, while the work of [14] proposed a robust alternating direction method of multipliers (ADMM)-based approach. On the other hand, the problem of SINR balancing targets at maximizing the minimum SINR among the users. In this objective, the available power is used to provide the best possible Quality of Service (QoS) for the users in a fair manner. Coordinated beamforming for this optimization target has been studied, e.g., in [15].

#### Related Work

The previously presented literature assumes an independent data stream for each user, which is known as unicast beamforming. Another interesting technique which has also been recognized in LTE standards is called multicast beamforming, where one stream is intended to a group of users. It has a great potential to address the nature of future traffic demands, e.g., to support demanding video broadcasting applications. Furthermore, physical layer multicasting has applications in satellite communications [18]. For the power minimization target, a physical layer multicasting problem was originally proposed in [20], proven NP-hard and accurately approximated by the semidefinite relaxation (SDR) and Gaussian randomization techniques. Bornhorst *et al.* [21] proposed a distributed successive convex approximation method to solve the multicasting problem in relay networks. In [22], a unified framework was derived for physical layer multigroup multicasting, where independent sets of common data are transmitted to different interfering groups of users. They proved the NP-hardness of the sum power minimization and the minimum SINR maximization problems in a multicast multigroup system with a sum power constraint (also known as the QoS and the max-min fair problems), and proposed accurate approximations to solve the problems. The works of [23] derived a consolidated solution for the weighted max-min fair multigroup multicast beamforming under per-antenna power constraints. This work was extended for the sum rate maximization problem in [25], where the important aspect of user scheduling was also considered. The coordinated multicast beamforming problem with a single group per cell was considered in [26], where the authors derived a distributed and centralized method for the QoS and max-min fair problems, respectively. In [27], the authors proposed a centralized algorithm for the energy efficiency maximization problem in a multi-cell multigroup system. In the literature, however, there is a lack of generic centralized and distributed algorithms for the sum power minimization and SINR balancing problems in a multi-cell multicast system with multiple groups per cell.

#### Contributions

In this paper, we propose centralized and distributed beamforming designs for multi-cell multigroup multicast systems. Two different optimization targets are considered, i.e., power minimization and SINR balancing. We first generalize the methods of [22] from single-cell to a multi-cell system model (i.e., the centralized methods) and [26] to a multigroup system model (i.e., centralized and primal decomposition -based distributed method for sum power minimization and a centralized method to SINR balancing problem). To be more specific, in the first design, we aim at minimizing the total transmission power of the system while guaranteeing predefined minimum SINRs for all the active users. The resulting problem is non-convex and difficult to tackle as such. We propose to use SDR to approximate the problem as a semidefinite program (SDP) [22]. If the SDR results in solutions where all the covariance matrices are rank-1, then the algorithms are optimal for the original problem. Otherwise, the Gaussian randomization method is utilized to provide a sub-optimal, but feasible, beamforming solution. However, the simple power scaling as in [26] cannot be applied to a multi-group setting due to multiple power variables in each cell. Thus, differently from [26], we propose a linear programming method to perform Gaussian randomization. We further propose two alternative methods to obtain a distributed beamforming implementation. Firstly, we generalize the primal decomposition method of [26] to a multigroup system model by reformulating the one-level SDP into two optimization levels. In the higher level, upper bounding levels for inter-cell interference powers are optimized while the lower level is in charge of optimizing the beamformers for a given set of interference power constraints. Secondly, differently from [26], we propose another distributed solution via an alternating direction method of multipliers, where the inter-cell interference optimization is divided into local and global optimization by forcing the equality via consistency constraints. The distributed methods require only local CSI at each BS and low-rate backhaul scalar information exchange. In the second problem formulation, we consider a multicast SINR balancing problem to obtain fairness between the users. More specifically, the aim is to maximize the minimum SINR of all the users in the network. In the proposed centralized design, we again generalize the method of [22] from a single-cell to a multi-cell and [26] to a multigroup system model by reformulating the original nonconvex max-min fractional program based on SDR. Then, the problem is reformulated to find a bisection type algorithm where the SDP feasibility problem is solved in each iteration. Subsequently, differently from the earlier works, by limiting the maximum inter-cell interference towards each user at each BS, we propose a simple distributed beamforming design which can be implemented without any information exchange among the cells. The superiority of the proposed algorithms over the conventional transmission schemes is demonstrated via numerical examples. We also compare the convergence rates of the distributed methods and illustrate when the proposed algorithms are optimal, i.e., produce rank-1 solutions.

Parts of this paper have been published in the previous conference publication [28], i.e., the centralized approach and the primal decomposition based distributed implementation for the sum power minimization problem with preliminary simulation results. The following additional contributions can be found in this paper which have not been in the previous work. Here, we propose an alternative distributed beamforming design based on the ADMM method for the power minimization problem. We further provide a centralized and a distributed algorithm for the SINR balancing problem, and show how the Gaussian randomization is performed in this formulation. Finally, we have extended the simulation model to take into account varying interference conditions, i.e., we have added a cell separation parameter into simulations. Using this, we have provided a more extensive set of simulations to illustrate the effectiveness of the proposed methods.

#### Organization and Notation

The paper is organized as follows. In Sections Section 2 and Section 3, the multi-cell multigroup multicast system is introduced, and the optimization problems are formulated, respectively. In sections Section 4 and Section 4.2, the centralized and distributed beamforming algorithms for the power minimization problem are derived, respectively. Section 5 presents the centralized and distributed algorithm for the SINR balancing problem. The performance of the proposed algorithms is examined in Section 6 via numerical examples. Finally, conclusions are drawn in Section 7.

The following notations are used. Bold face lower and upper case characters denote column vectors and matrices, respectively. The operators and correspond to the conjugate transpose and the trace operator. Calligraphic letters are for sets and notation means that is positive semidefinite.

## 2System Model

The multi-cell multigroup multicasting system consists of BSs, groups and users. Each BS has transmit antennas, whereas each user is equipped with a single receive antenna. We denote the sets of BSs, groups and users by , and , respectively. Each group of users is served by an independent data stream transmitted from a single BS. Therefore, there exists inter-group interference between the groups of the same cell and between the groups belonging to different cells. The former interference type is known as intra-cell interference and the latter as inter-cell interference. The number of groups served by BS is given by , and the corresponding set of groups is denoted by . The number of users in group is denoted by , and the corresponding set of users is given by . Note that each user belongs only to a single group. Thus, the sets of users belonging to different groups are disjoint. This can be mathematically expressed as , . The signal received by user is given by

where is the channel vector from th BS to th user, is the unnormalized beamformer of group , is the corresponding normalized information symbol and is the complex realization of white Gaussian noise with zero mean and variance .

Gaussian codebook is assumed to be used for each data stream. Thus, the rate for user can be expressed as

where denotes the SINR of user . The mathematical expression of is written as

where user belongs to a group .

## 3Problem Formulation

We consider two different optimization targets. In the first one, the (network-level) optimization goal is to minimize the sum transmission power of the BSs while satisfying the user-specific rate targets for active users. The mathematical expression of the problem is given by

where is a predefined rate target for user . In the multiple-input single-output (MISO) problem setting in , there exists a direct mapping between the user-specific rate and SINR. Hence, the rate targets can be equivalently changed into SINR targets . The resulting problem is given by

where is the fixed SINR target. Note that may be infeasible for some system settings and/or channel conditions. This can happen in interference limited scenarios, e.g., if the number of antennas at some BS is smaller than the total number of users in neighboring cells plus the sum of the user groups served by BS . Nevertheless, admission control is responsible for handling infeasibility issues by alleviating the system-specific requirements such as reducing the number of simultaneously scheduled users or decreasing the minimum SINR values. [29]. In [9] and [26], the feasibility issue was considered for unicast and multicast beamforming concepts, respectively. A joint admission control and multicast beamforming problem was solved and discussed in [30]. The feasibility of is assumed in the rest of this paper.

Problem is difficult to solve as such because it is non-convex and NP-hard. In fact, a conceptually simple single-cell multicast problem was shown to be NP-hard in [20], and is a more generic multi-cell version of the problem in [20]. It is worth mentioning that in case of heterogeneous power requirements between the BSs, BS-specific priority weights could be easily included in problem . For the ease of notation, however, they are left out in this paper. On the other hand, we could easily include different power constraints in , e.g., BS-specific, antenna-specific or antenna group -specific power constraints. All of these constraints admit an SDP formulation and they are readily separable between the base stations. By adding the given power constraints will reduce the size of the feasible solution set of .

The second problem of interest is a coordinated multicast SINR balancing problem where we aim at providing fairness among the users by maximizing the minimum SINR of the users, subject to BS-specific power constraints. The problem can be cast as

where is a maximum transmit power allowed for BS . Similarly to the power minimization problem, is hard to tackle as such due to its non-convexity and NP-hardness. Note that we could add the priority weights for each user in problem if some users are more important than others. However, the weights are left out for the ease of notation. The algorithms proposed in this paper can be straightforwardly generalized to different power constraints, e.g., for antenna-specific or antenna group -specific power constraints.

## 4Power Minimization Problem

Here we aim at minimizing the total transmit power of all the BSs with a given user-specific SINR constraints. We first derive a centralized algorithm, which is followed by the primal decomposition -based and ADMM-based distributed methods, respectively. Subsequently, the distributed Gaussian randomization method is presented for the cases where the solution is not rank-1. Finally, practical aspects of the distributed methods are discussed, e.g., the effects of limiting the number of iterations and signaling overhead caused by the algorithms.

### 4.1Centralized Beamforming Design

Due to the non-convexity of problem , a more tractable formulation which can be solved efficiently is desired. To this end, the SDR method can be used to obtain a convex approximation of problem . Specifically, the term is replaced with a semidefinite matrix , . The resulting convex SDP can be stated as

where . In order to solve , global CSI is required at a central controlling unit or at each BS. If all the optimal covariance matrices produced by are rank-1, the solution is also optimal for the original problem . In this case, the eigenvalue decomposition can be used to extract the optimal beamformers from . The resulting beamformers are given by , , where and are the principal eigenvalue and eigenvector of . However, due to the relaxation, the rank of can be larger than one. Consequently, an optimal solution of is not necessarily optimal for .

For some specific problems in the literature the SDR provides optimum solutions, e.g., the optimal unicast beamforming in [8]. Generally, however, the unit rank solution of the relaxed problem cannot be guaranteed due to the NP-hardness of the multicast problem. To this end, a general approach is to apply a rank-one approximation over the higher rank solution. In case of multicast beamforming, the Gaussian randomization method has been shown to result in the highest accuracy [31]. Let us denote by the optimal symmetric positive semidefinite matrices of the relaxed problem. In the Gaussian randomization, we generate a rank-one beamforming approximation of the original problem as a complex Gaussian vector with zero mean and covariance equal to , i.e. , . Vector is called a candidate beamforming solution. For the candidate beamformers, however, the feasibility of the solution to the original problem is not guaranteed. Thus, the feasibility of each candidate beamformer has to be checked by solving a power minimization problem. The aim is to minimize the transmit powers with the user-specific SINR targets , for each fixed candidate beamformer (normalized to have unit power). This yields a linear program (LP)

After solving , the beamforming solution for each set of candidate beamformers is obtained as , , where is the optimal power associated with fixed candidate beamformer . Finally, after generating a predetermined number of candidate solutions, the one that yields the lowest objective value of the original problem is chosen. Note that the beamformers are sub-optimal, but feasible, for the original problem. The accuracy of this approximate solution is measured by the distance of the approximate objective value and the optimal value of the relaxed problem. This accuracy increases with the increasing number of Gaussian randomizations. The proposed centralized multicast approach is summarized in *Algorithm ?*. The algorithm requires the knowledge of the global CSI either in the central controlling unit (if the controller performs the optimization) or in each BS (if the CSI is shared between the BSs). In the latter case, the BSs run *Algorithm ?* in parallel, i.e., at the same time. It is worth mentioning that an alternative iterative solution based on iterative convex approximation without introducing the SDR was proposed in [32]. However, the advantages of the SDR approach lie in the fact that only a single problem needs to be solved if the solution is rank-1, and even the Gaussian randomization requires solving a linear program which can be implemented efficiently. The method of [32] can be more useful in the cases where the ranks of the covariance matrices would be far from rank-1.

### 4.2Distributed Beamforming Designs

The algorithm derived in the previous section requires centralized processing due to the coupling ICI terms in . In this section, we derive distributed methods which decouple the problems so that each BS can perform beamformer optimization based on local channel state information and low-rate scalar information exchange. Two alternative formulations are proposed, one based on primal decomposition and the other based on alternating direction method of multipliers. We derive the algorithms, explain the distributed implementations, summarize the details in step-by-step algorithms and discuss some practical aspects. The Gaussian randomization method is presented in case the solution is not rank-1.

To apply distributed optimization, we first reformulate . In this respect, we separate interference power to intra-cell and inter-cell terms, and add auxiliary variables to denote the ICI terms. Now, the coupling is transferred from beamformers to ICI variables. The reformulated problem is expressed as

where is the user set served by BS , is the ICI from BS to user , and the vector consists of all ICI variables. Since the inequality constraints are met with equality at the optimal solution, and are equivalent. Problem still requires central processing but it will be the basis for deriving the distributed methods in the next sections.

#### Primal Decomposition -based Approach

Primal decomposition is applicable to decouple an optimization problem which has coupling constraints [33]. Specifically, the specific problem decouples by fixing the coupling constraints. Primal decomposition divides the one-level optimization problem into two levels, i.e., the lower level subproblems and the higher level master problem. Herein, we generalize the method of [26] to a multigroup system model.

In the primal decomposition method, is divided into a two-level optimization. Firstly, by fixing the inter-cell interference levels, the beamformers are solved from the BS-specific subproblems by locally minimizing the BS-specific transmit powers. Secondly, the interference levels are optimized by solving a network wide master problem. The local subproblem for BS is formulated as an SDP

which can be optimally solved using any SDP solver. The network wide master problem can be cast as

where I is a set of positive -dimensional real vectors, and the optimal objective value of for given is denoted by . The vector with length includes the inter-cell interference terms which are associated with BS . Problem to find the inter-cell interference variables can be solved by using the projected subgradient method

where is the projection onto a positive orthant, is the iteration index, is the step-size and is the subgradient of at point . Since problem is convex, we can find the subgradient via the dual problem of as in [13]. The resulting subgradient at point is given by , where is the dual variable associated with in the SINR constraint of user at its serving BS (i.e., in subproblem ) and is the dual variable associated with in the inter-cell interference constraint of user at the interfering BS (i.e., in subproblem ). Because is an SDP, solving it using standard SDP solvers returns the optimal dual variables as side information (i.e., a certificate for optimality). Another method to find the dual variables is to formulate and solve the dual problem of .

The step-size of the projected subgradient method has to be properly chosen [33] in order to solve the master problem optimally. One example of the valid step size rule is, e.g., a nonsummable diminishing step size, satisfying and . If all the optimal covariance matrices returned by the proposed method have unit ranks, then this solution is also optimal for the original problem . In this case, the optimal beamformers are obtained from by applying the eigenvalue decomposition, i.e., , , . In case at least one of has a rank higher than one, then the SDR is not optimal. Then we can use the Gaussian randomization method presented in Section ? to find feasible rank-one beamformers.

#### Distributed Implementation

When local CSI is available at each BS and the BSs are allowed to exchange scalar information via low-rate backhaul links, we can develop the following distributed beamforming implementation. The subproblem in and the corresponding part of the master problem in , i.e., the update of , are solved independently at BS , for all in parallel. At subgradient iteration , BS signals the dual variables associated with the SINR constraints, i.e., , to all the interfering BSs using the backhaul links. On the other hand, each BS signals the dual variables associated with the ICI constraints, i.e., , only to the serving BS of each user , since each ICI variable couples only two BSs, the serving BS of user and the interfering BS . The total amount of the required backhaul signaling at each subgradient iteration is which is the sum of the real-valued terms exchanged between the coupled BS pairs, assuming a fully connected network and an equal number of users and groups at each cell (i.e., , ). Although the primal decomposition method results in optimal solution for , the BSs only know their own covariance matrices, while the other BSs may have higher rank covariance matrices. To this end, each BS can send a one-bit feedback to indicate if the rank is higher than 1. If Gaussian randomization is required, the BS-specific powers for each Gaussian randomization instance need to be shared among other BSs to select the best one in a distributed manner. The overall distributed approach is summarized in *Algorithm ?*, which is performed at BS , for all in parallel. The proposed algorithm applies a subgradient method to update the ICI variables for which the convergence to the optimal solution (for convex problems) has been analyzed, e.g., in [34], [35].

#### ADMM-based Approach

Next we present the proposed ADMM-based distributed solution to . As a first step, we equivalently reformulate as

where we have introduced two new variables and for each [11]. In the first constraint, each BS has replaced the ICI variable (i.e., the ICI that BS causes to the user of BS ) with , i.e., this interference is locally optimized. On the other hand, in the second constraint, each BS has replaced with . Now, let us consider a user which is served by BS . Then, the ICI which BS causes to user is but this value is locally optimized by both BS (this variable is ) and also BS (this variable is ). The last constraint of (13) then makes sure that . The introduction of the new variables is required to arrive at a form tractable to the application of ADMM as presented shortly. Intuitively, we can interpret as the local version of optimized by BS . A vector is composed of all the local variables which are associated with BS , i.e., the interference from BS to the users of other cells and the interference from other BSs towards the users of cell . This means that each interference term couples only two base stations. Similarly, involves the global variables associated with BS in a same order as in . With the new formulation in , let us define

to be a local feasible set for cell . As a result, problem can be compactly expressed as

At this point, it is observed that admits a form of standard global consensus problem which lends itself to the application of ADMM [36]. Specifically, the objective function and the constraints are separable, and then the consensus constraints force the equality between the local variables and the global variables . Let us define the local variables of BS as , and define

To apply the ADMM method, we first write

which is a partial augmented Lagrangian of . In , the dual variables related to the interference equality constraints are denoted by . The final term in with parameter penalizes for the violation of the equality constraints . The variable updates in ADMM can be written as:

where is the iteration index. From , the local variables are found from the optimal solution of the following convex problem

and then the update is written as . To find the global variables from , the updates are written as , where is the optimal solution of the quadratic unconstrained convex program

Problem is separable in and can be solved in parallel when the local variables have been exchanged between the coupling BSs. By solving the point of the zero gradient of gives

where and are the serving and interfering BS of user , respectively. Equation can be simplified as , because inserting to results in

From the above expressions, we can see that is actually a complement of , i.e., . The ADMM-based distributed method is compactly written in Algorithm ?. The convergence proof of ADMM to optimal solution, applied to the convex problem of type , can be found, e.g., in [36].

#### Distributed Implementation

As in the primal decomposition based approach, the distributed implementation is enabled if each BS acquires local CSI and scalar information exchange between the BSs is allowed via low-rate backhaul links. The first local subproblems are completely decoupled and can be solved independently in parallel at each BS, with fixed global variables . The update of global variables requires the knowledge of the local copies as shown in . At ADMM iteration , BS signals the local copies , to all the interfering BSs. Since each of inter-cell interference term couples only two BSs, the number of exchanged scalars per ADMM iteration is given by which is the same as in the primal decomposition based approach. Finally, each BS can again locally update the dual variables in with the updated from the previous step. Similarly to the primal decomposition -based method, a one-bit feedback can be transmitted to indicate the rank of the covariance matrices and possible Gaussian randomization requires the BS-specific powers to be shared between the BSs. The overall distributed approach is summarized in *Algorithm ?*, which is performed at BS , for all in parallel.

#### Distributed Gaussian Randomization

As mentioned earlier, the solution of the SDR is not optimal to the original problem if at least one of produced by distributed methods has a rank higher than one. In this case, we propose to use the Gaussian randomization method to find feasible beamformers. As in the centralized Gaussian randomization presented in Section 4, the optimal covariance matrices returned by the distributed methods are used to generate candidate (normalized) beamforming solutions , , , which are Gaussian random variables with zero mean and covariance . To guarantee the feasibility of the beamforming solution, each BS solves the following LP

After solving , BS can define its beamformers by , , where is the optimal power associated with the candidate beamformer . The resulting beamformers are sub-optimal for the original problem. After generating a predefined number of candidate solutions, the one that gives the lowest objective value of the original problem is selected. To find the lowest network-wide objective, the BS-specific power values have to be exchanged between the BSs. Note that the Gaussian randomization can be implemented in a fully distributed manner, since the ICI variables are fixed in . The inter-cell interference values are obtained from the optimal solution of and for primal decomposition based method and ADMM-based method, respectively.

#### Practical Considerations

To acquire optimal performance, *Algorithms ? and ?* need to be run until convergence (provided that the obtained covariance matrices are all rank-one). However, this is somewhat impractical since the more iterations are run, the higher the signaling/computational load and the longer the caused delay. In this respect, *Algorithms ? and ?* naturally lend themselves to a practical design where they can be stopped after a limited number of iterations to reduce delay and signaling load. Limiting the number of iterations comes at the cost of increased sum power. In the primal decomposition, since the inter-cell interference levels are fixed at each iteration, feasible beamformers can be computed via the eigenvalue decomposition or the Gaussian randomization procedure, depending on the rank properties of the covariance matrices. In case of ADMM, the local and global interference variables (i.e., and ) are different during the iterative processing. As a result, the feasibility of the SINR constraints in the original problem is not necessarily guaranteed before the convergence. However, since the global inter-cell interference terms in are inherently feasible, we can always get feasible beamformers even if the local inter-cell interference terms would be infeasible. Specifically, in order to guarantee a feasible beamforming solution at any iteration, each BS can fix the local interference terms to be equal to the global ones , and solve by keeping fixed.

The overhead of backhaul signaling for the centralized and distributed algorithms are compared in Table 1. In the centralized algorithm, we assume that global CSI is made available for each BS by exchanging local CSI of each BS via backhaul links. Thus, the total backhaul signaling load in terms of scalar-valued channel coefficients is given by if we assume equal number of users and groups per cell. Here, one complex channel coefficient is considered as two real-valued coefficients. For the distributed algorithm, the total backhaul signaling load is presented per subgradient iteration or ADMM iteration, because the total backhaul signaling per iteration is the same for the two methods. In Table 1, the values inside the brackets denote the percentage of the signaling load that the distributed methods require compared to the centralized algorithm. It can be observed that significant backhaul signaling savings are achieved with the distributed methods and the savings are significant even when multiple iterations are performed. To give an example about the signaling overhead between the centralized and distributed approaches, let us consider the case where (i.e., the last row of Table 1). In this case, we could run 64 iterations in the distributed methods and the signaling overhead would be still less than in the centralized method. The gap between the methods increases with the network size. The backhaul signaling overhead can be significantly reduced by limiting the number of iterations. It is worth observing that the signaling overhead of the centralized method scales both with the number of BS antennas and the number of users , while only with for the distributed methods. Thus, if the number of antennas is large, the distributed methods require less signaling even when the channels are static for longer time. Also, if the number of antennas grows high, it is obvious that the signaling overhead becomes a problem and it is relevant to reduce the signaling overhead. However, the benefits of the distributed methods become more important in the time-varying channels [11] because the BSs can easily acquire local channel information, and they only need to exchange the user-specific scalars. The centralized method would require sharing all the channel information every time when the channels change, causing significant signaling overhead and delays. In practice, building very high rate and low-latency optical fibers for backhauling could enable centralized methods to be practical, but this would result in very high costs for the operators [37]. However, with smaller backhaul overhead, wider range of capacity-limited backhaul technologies can be applicable, especially when the network dimensions increase [39].

A distributed approach allows some special case designs where the number of optimization variables is reduced, leading to a lower computational load and even a further decreased signaling overhead. These special case designs come at the cost of somewhat decreased performance. Some of the possible special cases are presented below:

Common interference constraint: .

Fixed interference constraints: , , where is a predefined constant. Does not require any backhaul signaling.

Inter-cell interference nulling, i.e., , . Does not require any backhaul signaling.

Centralized | |||
---|---|---|---|

256 | 16 (6.3) | 160 (63) | |

1728 | 48 (2.8) | 480 (28) | |

6144 | 96 (1.6) | 960 (16) | |

## 5SINR Balancing Problem

### 5.1Centralized Beamforming Design

Problem is a nonconvex max-min fractional program and difficult to solve as such. To find a more tractable formulation, we first write the equivalent transformation of as

Problems and are equivalent since the first constraint of holds with equality at optimum. By applying the SDR, we get

The above problem is an epigraph form of a quasiconvex optimization problem which can be optimally solved by a bisection search over [41]. Specifically, in the beginning, we initialize lower bound and upper bound and set . For fixed , we solve problem

which is an SDP feasibility problem. If the problem is feasible, we set . Otherwise, we set . Then we calculate new as and repeat until , where is a small threshold. It is worth mentioning that the same optimality/rank analysis as in power minimization problem applies also to the SINR balancing problem so that if the solution is not rank 1, we can use Gaussian randomization to find a suboptimal solution. In this case, instead of solving as in the power minimization problem, we solve the quasi-linear program

for a given set of normalized candidate beamformers and choose the one that gives the maximum . Problem can again be solved via the bisection method. The centralized algorithm is summarized in Algorithm ?.

### 5.2Distributed Beamforming Design

Finding a distributed solution for the SINR balancing problem is more challenging because in addition to the inter-cell interference terms as in the power minimization problem, global variable couples all the cells. Here we propose a simple distributed design where each cell maximizes its own minimum SINR with the constraint that the interference towards all the other cells are kept below a certain threshold. This method can be realized without any backhaul information exchange. Specifically, each cell locally solves the problem

where is the maximum predefined inter-cell interference from BS to user . The second constraint makes sure that the interference is kept below the limit . In this method, each cell gets slightly different values for the minimum SINR but it is simple and can achieve good performance if the interference levels are chosen reasonably. In practice, these values could be averaged over longer time and we could then easily use some average values. Testing e.g., long-term values for different system parameters in offline and using look-up tables to pick up good average values could be a practical option. In the distributed approach, the Gaussian randomization can be performed locally because the maximum inter-cell interference levels are fixed. Specifically, each BS uses the bisection method to solve the problem

for a given set of normalized candidate beamformers and the one that gives the maximum is chosen at each BS . The distributed algorithm is summarized in Algorithm ?.

## 6Simulation Results

This section numerically evaluates the performances of the proposed algorithms. In the simulation model, we assume a network of BSs with transmit antennas per each BS. The total number of groups in the network is denoted by , and the groups are then equally divided between the BSs. The number of users per each group is given by , where is the total number of users in the network. The channel conditions are assumed to be frequency-flat Rayleigh fading, and the channel coefficients between the antennas are uncorrelated. We define cell separation parameter to represent the average path loss between BS and users of neighboring cells of BS . This means that the interference power towards the users of other cells is attenuated by the value of . More specifically, the observed channel from BS to user () is , where each element of