Joint Precoder and Combiner Design for MMSE Distributed Beamforming with PerAntenna Power Constraints ^{†}^{†}thanks: This work was supported by the Defense Advanced Research Projects Agency (DARPA) under contract D17PC00006. The views, opinions and/or findings expressed are those of the authors and should not be interpreted as representing the official views or policies of the Department of Defense or the U.S. Government.
Abstract
We consider minimum mean square error (MMSE) joint precoder and combiner design for single and multi carrier distributed beamforming systems with nonuniform perantenna transmit power constraints. We show that, similar to the maximumgain problem, an iterative GaussSeidel algorithm can be used for minimizing MSE which alternately optimizes the transmitter and receiver coefficients. In a single carrier system the optimum transmit coefficients are obtained by a simple projection of the effective MISO channel. In the multicarrier case with a sumMSE objective, the GaussSeidel approach is once again applicable, but the transmit coefficients must be found by solving a quadratically constrained quadratic problem for which we apply a dual gradient algorithm. A numerical example is presented which shows improvement of 0.7 dB in carrier signaltonoise ratio (SNR) relative to a projected eigenvector method for a multicarrier DBF system with Rayleighfaded multipath channels.
I Introduction
Distributed beamforming (DBF), is a promising technology for increasing communication range, providing improved signaltonoise ratio (SNR), or reducing probability of intercept and detection for geographically separated groups of communicating nodes [1]. In a DBF system, multiple transmitters with potentially heterogeneous hardware form an antenna array in an adhoc manner. Each transmit radio is power constrained and since the group is heterogeneous, this imposes a set of nonuniform perantenna power constaints (PAPC) on the transmit array. On the receiving end, another group of radios forms a receive array with nonuniform noise profile, again due to heterogeneous hardware.
Several studies of multiple input multiple output (MIMO) communications with uniform and nonuniform PAPC have appeared recently for various scenarios and objectives. For example, in [2] the authors consider precoder design for a multiple input single output (MISO) system with an outage probability objective. In [3] a zeroforcing precoder is designed for a broadcast channel with sumrate objective. In [4], transmit beamformers are designed for a multiuser scenario with a signaltoleakageplusnoise criterion. In [5] multiuser sum rate is maximized, again for MISO channels. In [6] directional beamforming under PAPC is considered for MISO channels. In [7] a downlink cellular maxmin problem is formulated for SINR over a set of users. Joint transmit and receive beamforming optimization with PAPC has also been studied. Optimum precoders (in terms of beamforming gain) for various combining strategies with uniform PAPC, also called equal gain transmission (EGT) were presented in [8]. Zheng, et. al. [9] then proposed a cyclic (GaussSeidel) algorithm for joint optimization of the precoder and combiner weights, again with a gain objective. Since then, the GaussSeidel approach has been applied for multiuser MIMO cellular downlinks with sumMSE objectives in [10] and for multicarrier systems with arithmetic error probability objectives in [11]. In this paper we consider joint precoder and combiner optimization for single and multicarrier beamforming systems with PAPC and MSE objectives. We note that besides DBF systems, conventional (colocated) MIMO systems often have PAPC too and our formulation is equally applicable in these cases.
It is well known that in a narrowband MIMO beamforming system (i.e., where the precoder and combiner matrices consist of only one column), maximum gain is achieved when the precoder and combiner weight vectors lie in the directions of the channel’s dominant right and left eigenvectors, respectively. This principle is used to guide beamforming weight selection with a totaltransmitpower (across all antennas) constraint. However, when perantenna constraints are imposed, the gainmaximizing transmit and receive weights cannot be found in closed form. Instead one can appeal to the GaussSeidel approach which involves alternately optimizing the transmitter and receiver weights until convergence [9, 10, 12, 11]. With EGT constraints, and maxgain objective, the transmit weights are found by projecting the effective MISO channel onto the set of vectors with unitmagnitude components [8], while the receive weights satisfy the maximum ratio combining (MRC) principle and lie in the direction of the effective single input multiple output (SIMO) channel.
In this paper, we show that a GaussSeidel approach can be applied to find the MMSE transmit and receive weights for both narrowband and wideband (multicarrier) systems. We make the following contributions:

We show that, similar to the maximum gain case, the singlecarrier MMSE transmit weights under PAPC are obtained by a simple projection of the MISO channel but only when this channel satisfies a certain norm condition. We then show that a GaussSeidel algorithm can be used to jointly optimize the transmit and receive weights for MSE under PAPC.

We find the optimal transmit weights with PAPC for a multicarrier system with a sumMSE objective using a dual gradient algorithm which we then embed in a GaussSeidel algorithm to find the jointly optimal transmit and receive weights with PAPC.
We use the following notation. Boldfaced lowercase symbols represent complex vectors and boldfaced uppercase symbols represent complex matrices. The superscripts and denote transpose and conjugate (Hermitian) transpose of a matrix or vector, respectively and indicates the complex conjugate. The quantity is a diagonal matrix with on the diagonal and is a blockdiagonal matrix with on the diagonal. Statistical expectation is denoted by . The magnitude (modulus) of a complex number is denoted and the angle by . Vector norms are represented by with the subscript denoting the type of norm, and projection operators are denoted with the superscript indicating the projection set. Finally, and are the positive and nonnegative orthants of respectively.
The remainder of the paper is organized as follows. Section II presents the MIMO beamforming system model and introduces the projection and norm operators which facilitate the analysis in the sequel. Section III summarizes the maxgain problem which provides useful comparisons to the MMSE systems which we develop later. Next in Section IV we derive the jointly optimal MMSE transmit and receive weights for single carrier systems and Section V extends this to multicarrier systems with sumMSE objective. A numerical example is presented in Section VI which shows the utility of the proposed algorithms. Finally, the paper is concluded in Section VII.
Ii System Model
Iia PAPC MIMO Beamforming System
The transmitter, receiver MIMO beamforming system is described by the following equation
(1) 
where and are the complex transmitted and equalized symbols, respectively, is the receive noise vector, is the complex MIMO channel matrix, and and are the transmit and receive beamforming weight vectors, respectively. The multicarrier case is modeled by
where the subscript indicates the th carrier. We assume a distributed system in which some node, which we call the fusion center in keeping with standard DBF terminology, has full channel state information (CSI). The fusion center computes the jointly optimal transmit and receive weight vectors and feeds this information to the other nodes. We consider perantenna transmitter power constraints of the form
(2) 
where is the th component of the transmit weight vector and is the maximum allowable power of the th element. In the multicarrier case, we constrain the total power across carriers: where is the transmit weight of the th carrier on antenna .
IiB Projection and Norm
Define the constraint vector and matrix and the total power . Define as the set of feasible transmit vectors which meet the constraints (2) with equality (i.e., the boundary of the feasible set):
It can be easily shown that the closest to an arbitrary nonzero is
(3) 
and thus is a projection operator onto the set [13], which we call the projection. Equation (3) is valid for any norm as well as the following “weighted” norm [14], which we call the norm, and which will facilitate further analysis.
We note a few important facts regarding the projection and the norm. First, when , its norm is . When the power constraints are all unity, i.e., , the norm is equivalent to the norm. From (3), the projection is independent of the magnitudes of the components of and thus for any nonzero real . As with any norm, we have for any real . Finally, for any nonzero , we have
IiC Constrained MMSE Problems
In this paper, our goal is to minimize MSE subject to PAPC. From the narrowband model (1). The normalized MSE is given by
(4)  
where and is the noise covariance. When , the transmit power of the th antenna is and we will make this assumption henceforth. The singlecarrier constrained MMSE problem is
subject to  (5) 
With multiple carriers, we minimize the sumMSE across carriers and constrain the total power at each antenna. SumMSE has been used as an objective for multiuser systems [10], spatial multiplexing systems [15], and multicarrier systems [16]. The joint optimization problem is
subject to  (6) 
Iii Max Gain Beamforming with PAPC
The problem of maximizing the gain under EGT constraints has been studied extensively, e.g., [8], [9] and the focus of the present paper is MSE minimization. However, a brief summary of the results on gain maximization is in order here as there are many parallels to the MMSE case. Thus, in this section we provide this summary and in so doing, we also extend the EGT results to the case of nonuniform PAPC.
Iiia MaxGain Precoder and Combiner
The optimum unitnorm receive weight with transmit weight is the MRC vector
(7) 
where is an arbitrary phase shift. With no individual transmit power constraints, the optimum (maximal ratio transmission) unitnorm transmit weight vector for receive weight is
(8) 
where is again arbitrary. The jointly optimal weights are where and are the dominant left and right eigenvectors of .
IiiB Max Gain Precoder and Combiner with PAPC
With PAPC on the transmitter, the optimum transmit vector with receive vector is
(9) 
with arbitrary . Note that the optimum EGT weight vector given in [8] is a special case of the above equation, with . Assuming MRC receive weights, the optimum transmit phase vector solves (also a generalization of a result from [8]).
IiiC Joint Gain Maximization with GaussSeidel Algorithm
While a closedform solution for the jointly optimal weights has not been found, [9] proposed a GaussSeidel (or “cyclic”) algorithm wherein the transmit and receive weights are alternately updated according to (7) and (9). GaussSeidel algorithms are also referred to as block coordinate descent algorithms and their convergence is studied in [17] and [18]. They are appealing since at every iteration, the objective is guaranteed to improve or stay the same. However, convergence can only be to a local minimum (for objectives which are nonconvex in the joint set of decision variables), and even this convergence requires somewhat restrictive conditions [17, 10]. Nevertheless, the GaussSeidel approach has been applied successfully in several MIMO precoder/combiner optimization problems [9, 10, 11]. In the maxgain case, the algorithm starts with initializing , for example to . Then equations (9) and (7) are alternately applied until convergence.
Iv SingleCarrier MMSE Beamforming with PAPC
We now turn to the MMSE problem, which is the focus of this paper. In this section we study the singlecarrier case in depth and in the next section we study the multicarrier case. Our goal is to minimize the MSE subject to PAPC, i.e., to solve problem (5). Note that the MSE (4) is convex in for a fixed and vice versa, but it is not convex in [16], which provides some motivation for a GaussSeidel approach. We first summarize the MMSE problem without PAPC.
Iva MMSE Precoder and Combiner
The MMSE receive weight for a fixed transmit weight vector is
(10) 
The optimal transmit weight vector for fixed is
which is the same as (8) up to a constant. Joint optimization of and is accomplished by expressing the resultant MSE with the optimal receive vector (10) as
(11) 
which is minimized for a maximum total transmit power by where is the (unit norm) dominant eigenvector of [16].
IvB MMSE Precoder and Combiner with PAPC
IvB1 Constrained Problem
To find the jointly optimal transmit and receive weights for the MMSE problem with PAPC (5) we could try to maximize the denominator of the RHS of (11) under perantenna power constraints, but the problem would not be convex. Instead we attempt to find a formula analogous to (9) which optimizes the transmitter weights for fixed receive weights with MSE objective. This, along with (10) can then be used in a cyclic algorithm to jointly optimize and for MSE. From (4), we see that minimizing with respect to is equivalent to minimizing , where
(12) 
is the effective MISO channel. The transmitter optimization problem is
subject to  (13) 
where , and is the th standard basis vector of . In Appendix A we show that this quadratically constrained quadratic problem (QCQP) exhibits strong duality. The solutions are given in the following three propositions, which are proved in Appendix A. These solutions depend on the MISO channel given in (12).
Proposition 1
If , then is an optimal powerconstrained transmit beamforming weight vector, i.e., a solution to (13), with nonnegative Lagrange multipliers , .
Proposition 2
If , then is a solution to (13) with Lagrange multipliers , where .
Proposition 3
If for all , then is a solution to (13) with Lagrange multipliers , where .
Of the solutions above, only the first (Proposition 1) has active constraints. This solution is also nearly identical to the maxgain under PAPC solution (9). However this solution is only valid when the norm of the effective MISO channel is less than unity. We show next, however, that in the context of a GaussSeidel algorithm for joint optimization of and , this condition is easily managed and the solution of Proposition 1 can be used throughout the algorithm’s iterations.
Proposition 4
If , then is the optimal powerconstrained transmit beamforming vector for receive vector
and the MSE with and is lower than any achievable MSE with receive vector . That is .
IvB2 Joint Transmit / Receive Optimization
From Proposition 4, if , we know that the jointly optimal solution cannot include . From Proposition 1, if , is optimal for . Thus the optimal pair must satisfy , and . Furthermore, when is equal to the projection, we have . Thus the joint optimization problem is
(14)  
subject to 
This problem is nonconvex but has a simple solution in the MISO case (), with receiver noise variance , channel , and scalar receive weight . The MISO problem is independent of and convex in with solution
with arbitrary.
IvB3 Shadow Prices
The resultant MSE with is given above in (14) and is differentiable as a function of with . The “shadow price” interpretation of is the reduction in MSE that can be realized per unit of power relaxation of the th power constraint (for small relaxations) [13]. The largest shadow price occurs for the element that maximizes . Thus, a lowpower element which experiences large gain to the beamformed receiver would do well to increase its power beyond its constraint, if possible.
IvC GaussSeidel MMSE Algorithm
Recall that the cyclic algorithm attempts to find a jointly optimal pair by sequentially optimizing and . Optimization of can be carried out using (10) and can be optimized using Proposition 1. However the situation may arise where . But from Proposition 4 we know that in such a case the MSE can be reduced by updating to , and then updating to . But since , there is no need to update to . Thus the cyclic algorithm only needs to sequentially update and using and (10). The transmit vector can be initialized to (see Section IVA).
V Multicarrier PAPC Beamforming with SumMSE Objective
In the multicarrier case, the goal is to solve problem (6). Once again, joint optimization of transmit and receive coefficients is a nonconvex problem and we appeal to the GaussSeidel method. Optimization of the receive weights given a set of transmit weights is a separable problem. Each carrier’s receive weight can be computed using (10). Thus, this portion of the GaussSeidel algorithm is easy. It remains to find the optimum PAPC transmit weights given a set of receive weights. The multicarrier precoder optimization problem is
subject to  (15) 
This problem is a convex QCQP and admits strong duality, as in the single carrier case, but the constraints introduce coupling in the variables. A primal gradient projection algorithm could be used, but the projection step would be difficult. On the other hand, the dual problem has feasible region , for which projection is simple. Thus we solve the dual problem with a gradient projection algorithm. The dual function, derived in Appendix B, is
where . The dual problem is
The dual function is concave and differentiable on and projection onto the orthant is simple. Primal variables are calculated with the Lagrangian minimizers given in (21) and (22).
We now have the necessary tools to jointly optimize . In Algorithm 1 optimization steps for and are alternated. First the transmit weights are initialized to the projected dominant eigenvectors of each carrier and then assigned equal power (Step 1). Then the cyclic iterations begin. Optimization of the receive weights is a separable problem, where each carrier’s weight is found using as in (10). Next, the dual variables are initialized and updated using gradient ascent on the dual function followed by projection on the orthant . The step size is found using a line search. Finally, the updated transmit weights are found from the Lagrangian minimization equation (21) as shown in Step 9, (or (22) if lies on the boundary of the orthant). Note .
Vi Numerical Example
In this section we demonstrate the utility of the cyclic algorithm through a numerical example. We consider two distributed groups of nodes with transmitters and receivers. The transmit power constraints (in Watts) are chosen uniformly in and the receive noise variances (in dBW) are uniform with spread 10 dB. The noise variance is frequency independent so and the noise between receivers is uncorrelated so is diagonal. The system uses carriers over MHz centered at GHz. We assume the receivers are spaced sufficiently far apart such that for each transmitter/receiver pair, we have an independent Rayleigh channel with exponential intensity profile and delay spread of 4.
We compare the multicarrier cyclic algorithm to three suboptimal approaches and to the optimum totalpowerconstrained weights [16, eqn 12]. The latter approach provides a lower bound on MSE since its feasible region contains the feasible region of (15). The first suboptimal approach uses . That is, for each carrier we find the optimum unconstrained (without PAPC) weight, which is the dominant eigenvector of , and project this weight on the feasible set. No attempt is made to optimally allocate power across carriers, and instead all carriers are allocated equal power. Note that this choice of transmit vectors is also used to initialize the cyclic algorithm (Step 1 of Algorithm 1). In the second suboptimal method, the singlecarrier cyclic algorithm from Section IVC is run independently for each carrier and once again equal power is assigned to each carrier. Thus this approach should yield better performance than the first, since each carrier’s transmit and receive weights are jointly optimized. But once again the transmit power is not allocated optimally across carriers. Thirdly we consider a naive approach where we find and for every antenna in violation of its perantenna constraint, scale the weight magnitude (equally across carriers) so that the constraint is met with equality. In implementing Algorithm 1, we terminate the cyclic loop after 20 iterations and the dual gradient loop after 200 iterations.
In Figure 1 we plot the empirical cumulative distribution functions (CDFs) of the carrier SNR, with 400 MonteCarlo trials, using the approaches described above. The figure shows that the singlecarrier cyclic approach is only marginally better than the projected eigenvector approach. This is consistent with results reported for the maxgain case which suggest that the projected dominant eigenvector is nearly optimal for the singlecarrier problem [12]. Returning to the figure, we see that the cyclic multicarrier algorithm improves median SNR by approximately 0.7 dB relative to the first two suboptimal approaches. The naive method of scaling the out of tolerance weights of performs very poorly.
The totalpowerconstrained approach performs 0.2 dB better than the multicarrier cyclic algorithm, but this marginal improvement comes at a steep price in terms of perantenna constraint violations. Figure 2 sheds light on these violations by showing the CDFs for the number of antennas in violation and for the maximum percent violation. There are always at least 7 antennas in violation and the median value of the maximum percent violation is well over 200%. Thus this method is not feasible with perantenna constraints.
Vii Conclusion
We have presented a method for joint transmit and receive beamforming optimization with MSE objective and nonuniform PAPC for both single and multicarrier systems, which have applicability for distributed beamforming systems. We showed that the optimum transmit weights under PAPC are nearly identical to the gain maximizing weights under a certain norm condition on the effective MISO channel. We also developed a GaussSeidel algorithm for joint MSE optimization of the transmitter and receiver weights, which is not encumbered by this norm condition. Optimum transmit weights for multicarrier systems with sumMSE objective were found using a dual gradient algorithm which solves a QCQP and another GaussSeidel algorithm was developed for the multicarrier case. Finally, we showed through numerical example the benefits of the method by comparison to several suboptimal approaches. The numerical example also demonstrated the disadvantage of a totalpowerconstrained approach.
Appendix A Proofs of Propositions 1–3
Proof of Proposition 1
First we note that the problem (13) can be converted into an equivalent problem with all real variables as in [19]. Define the operator which maps complex vectors in to real vectors in , and complex matrices in to real matrices in according to:
It is easy to show that is positive semidefinite if and only if is positive semidefinite Now (13) can be converted to the following real problem
subject to  (16) 
Since and are positive semidefinite, and are also positive semidefinite. Problem (16) has a convex objective and convex inequality constraints. Therefore it is a convex problem [13]. The Lagrangian function is given by
(17)  
where and are the Lagrange multipliers associated with the magnitude constraints. Finally and .
Since (16) is a convex problem and a strictly feasible exists (for example, the vector with ), Slater’s condition is satisfied, the problem admits strong duality, and the following KarushKuhnTucker (KKT) conditions are thus sufficient for the primaldual pair to be optimal [13]:


satisfies all inequality constraints in (16)


From (17) the zerogradient condition is . Thus, in terms of complex variables, KKT condition 1 can be written , or
which is equivalent to
(18) 
It is easily verified that this equation is satisfied for
(19) 
Furthermore, this choice of has magnitude . Thus the first three KKT conditions hold. By the assumption , the Lagrange multipliers are nonnegative. Thus the solution (19) is optimal.
Proof of Proposition 2
Proof of Proposition 3
Apply KKT conditions again.
Appendix B Derivation of Multicarrier Dual Function
The primal problem given in (15) can be converted to a real problem using the method in Appendix A. This problem has Lagrangian
which is quadratic in . The dual function is the infimum of the Lagrangian over all transmit vectors. The matrix coefficient is positive definite when all Lagrange multipliers are positive and positive semidefinite when at least one Lagrange multiplier is zero. (We assume dual feasibility so no Lagrange multipliers are negative). The Lagrangian is minimized when for all , or, in terms of the complex variables, when
(20) 
In the positive definite case () the minimizers are
(21)  
Inserting these values into the Lagrangian gives the dual function for
For the positive semidefinite case, assume there is a zerovalued Lagrange multiplier: . Then either of the following sets of vectors satisfy (20) and thus minimize the Lagrangian
(22) 
(where is the th component of and is the th standard basis vector of ). The dual function is thus
when is on the boundary of the nonnegative orthant . The dual function is differentiable on since the minimizers in (21) are unique [20]. However, we cannot claim differentiability when is on the boundary of the orthant since the minimizers in (22) are not unique.
References
 [1] R. Mudumbai, D. R. Brown III, U. Madhow, and H. V. Poor, “Distributed transmit beamforming: challenges and recent progress,” IEEE Commun. Mag., vol. 47, no. 2, pp. 102–110, Feb. 2009.
 [2] M. Medra and T. N. Davidson, “Lowcomplexity robust MISO downlink precoder design with perantenna power constraints,” IEEE Trans. Signal Process., vol. 66, no. 2, pp. 515–527, Jan. 2018.
 [3] T. M. Pham, R. J. Farrell, J. Dooley, E. Dutkiewicz, D. N. Nguyen, and L.N. Tran, “Efficient zeroforcing precoder design for weighted sumrate maximization with perantenna power constraint,” IEEE Transactions on Vehicular Technology, 2017.
 [4] H. Shen, W. Xu, A. L. Swindlehurst, and C. Zhao, “Transmitter optimization for perantenna power constrained multiantenna downlinks: An SLNR maximization methodology,” IEEE Trans. Signal Process., vol. 64, no. 10, pp. 2712–2725, May 2016.
 [5] P. L. Cao, T. J. Oechtering, R. F. Schaefer, and M. Skoglund, “Optimal transmit strategy for MISO channels with joint sum and perantenna power constraints,” IEEE Trans. Signal Process., vol. 64, no. 16, pp. 4296–4306, Aug. 2016.
 [6] M. Ding, H. Luo, and W. Chen, “Polyblock algorithmbased robust beamforming for downlink multiuser systems with perantenna power constraints,” IEEE Trans. Wireless Commun., vol. 13, no. 8, pp. 4560–4573, Aug. 2014.
 [7] G. Dartmann, X. Gong, W. Afzal, and G. Ascheid, “On the duality of the max–min beamforming problem with perantenna and perantennaarray power constraints,” IEEE Transactions on Vehicular Technology, vol. 62, no. 2, pp. 606–619, 2013.
 [8] D. J. Love and R. W. Heath, “Equal gain transmission in multipleinput multipleoutput wireless systems,” IEEE Transactions on Communications, vol. 51, no. 7, pp. 1102–1110, 2003.
 [9] X. Zheng, Y. Xie, J. Li, and P. Stoica, “MIMO transmit beamforming under uniform elemental power constraint,” IEEE Transactions on Signal Processing, vol. 55, no. 11, pp. 5395–5406, 2007.
 [10] S. Shi, M. Schubert, N. Vucic, and H. Boche, “MMSE optimization with perbasestation power constraints for network MIMO systems,” in Communications, 2008. ICC’08. IEEE International Conference on. IEEE, 2008, pp. 4106–4110.
 [11] Q. Huang, M. Ghogho, Y. Li, D. Ma, and J. Wei, “Transmit beamforming for MISO frequencyselective channels with perantenna power constraint and limitedrate feedback,” IEEE Transactions on Vehicular Technology, vol. 60, no. 8, pp. 3726–3735, 2011.
 [12] S.H. Tsai, “Transmit equal gain precoding in rayleigh fading channels,” IEEE Transactions on Signal Processing, vol. 57, no. 9, pp. 3717–3721, 2009.
 [13] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge University Press, 2004.
 [14] E. J. Candes, M. B. Wakin, and S. P. Boyd, “Enhancing sparsity by reweighted minimization,” Journal of Fourier analysis and applications, vol. 14, no. 56, pp. 877–905, 2008.
 [15] A. J. Paulraj and C. B. Papadias, “Spacetime processing for wireless communications,” IEEE Signal Processing Magazine, vol. 14, no. 6, pp. 49–83, 1997.
 [16] D. P. Palomar, J. M. Cioffi, M. A. Lagunas, and A. P. Iserte, “Convex optimization theory applied to joint beamforming design in multicarrier MIMO channels,” in Communications, 2003. ICC’03. IEEE International Conference on, vol. 4. IEEE, 2003, pp. 2974–2978.
 [17] D. P. Bertsekas, Nonlinear Programming. Athena Scientific Belmont, 1999.
 [18] P. Stoica and Y. Selen, “Cyclic minimizers, majorization techniques, and the expectationmaximization algorithm: a refresher,” IEEE Signal Processing Magazine, vol. 21, no. 1, pp. 112–114, 2004.
 [19] Z.Q. Luo, W.K. Ma, A. M.C. So, Y. Ye, and S. Zhang, “Semidefinite relaxation of quadratic optimization problems,” IEEE Signal Processing Magazine, vol. 27, no. 3, pp. 20–34, 2010.
 [20] M. S. Bazaraa, H. D. Sherali, and C. M. Shetty, Nonlinear Programming: Theory and Algorithms. John Wiley & Sons, 2013.