### Abstract

We consider the multiuser beamforming problem for a multi-input single-output downlink channel that takes into account the errors in the channel state information at the transmitter side (CSIT). By modeling the CSIT errors as elliptically bounded uncertainty regions, this problem can be formulated as minimizing the transmission power subject to the worst-case signal-to-interference-plus-noise ratio constraints. Several methods have been proposed to solve this nonconvex optimization problem, but none can guarantee a global optimal solution. In this article, we consider a semidefinite relaxation (SDR) for this multiuser beamforming problem, and prove that the SDR method actually solves the robust beamforming problem to global optimality as long as the channel uncertainty bound is sufficiently small or when the transmitter is equipped with at most two antennas. Numerical examples show that the proposed SDR approach significantly outperforms the existing methods in terms of the average required power consumption at the transmitter.

### Introduction

Consider the transmit beamforming problem in a multi-user multi-input single-output (MISO) downlink channel where the transmitter equipped with multiple-antennas sends independent data steams to a number of receivers, each having a single receive antenna. For such a system, a popular measure of quality of service (QoS) is the signal-to-interference-plus-noise ratio (SINR) [1] and optimal downlink transmit beamforming can be formulated as minimizing the transmission power subject to the SINR constraints.

It is well known that if the perfect channel state information is available at the transmitter (CSIT), the optimal transmit beamforming problem can be solved efficiently via convex optimization [1-4]. However, perfect CSIT is impossible to obtain in practice due to the finite length of training signal and/or limited feedback bandwidth [5,6]. This is unfortunate, because we cannot simply ignore the CSIT imperfection in the design of downlink beamformers, for otherwise the resulting solution might violate the SINR constraints for some users.

As a remedy, various researchers have proposed robust beamforming formulations by taking the CSIT imperfection into consideration. There are two major types of robust downlinks beamforming design: the stochastic robust design and the worst-case robust design. For the stochastic robust beamforming design (see [7-13]), the CSIT errors are modeled as random variables with known statistical properties. In this case, the design goal of transmit beamforming is to minimize transmission power while ensuring that all the receive SINR requirements are satisfied with high probability. In contrast, the worst-case robust beamforming design (see [14-20]) assumes that the CSIT errors belong to some known bounded uncertainty sets. The robust beamforming vectors are designed to satisfy the SINR requirements for all possible channel realizations in the uncertainty regions. For instance in [18,19], the authors formulated the robust beamforming problem as minimizing the transmission power subject to the worst-case SINR constraints, while assuming the CSIT errors lie in some elliptically bounded regions. Unfortunately the resulting robust downlink beamforming problem is nonconvex and has infinitely many constraints. It is difficult to solve in practice.

To overcome the computational difficulties associated with the robust downlink beamforming, researchers have proposed various optimization methods to approximate the problem by restricting the non-convex feasible region to a smaller but convex set [18,19]. Due to convex restriction, these methods typically can not find the globally minimum transmission power although they do guarantee the satisfaction of the QoS constraints. In fact, the computed beamforming solutions can be highly suboptimal.

In this article, we consider a semidefinite programming relaxation (SDR) method to solve the worst-case robust downlink beamforming problem for the MISO channel system. The SDR method has been extensively used in addressing the robust beamforming problem [21,22]. A critical issue of SDR is whether it can exactly solve the original problem, which is equivalent to whether the method can produce a rank-one solution. The key result of this article is that the robust beamforming problem can be globally solved by SDR method (i.e., SDR is tight) under elliptically bounded CSI errors, provided that the size of channel errors is sufficiently small or if the transmitter is equipped with no more than two antennas. We give a computable bound which ensures the tightness of the SDR if the radius of the error region is less than this bound. In the simulations, this computable bound is always reasonably large, which makes the tightness result valuable. Moreover, the SDR method significantly outperforms the existing robust beamforming methods based on the worse-case design [18,19].

The organization of the rest of the article is as follows. In Section “System model and problem formulation”, we present the system model and the problem formulation. Section “relaxation and its dual” proposes the SDR relaxation of robust beamforming problem and its dual. We provide our main results in Section “Main results”. The cases that the base station (BS) is equipped with two antennas and three antennas are considered in Sections “Special case: two transmit antennas” and “Special case: three transmit antennas”, respectively. Simulation results are presented in Section “Simulation results” and conclusions are given in Section “Conclusions”.

*Notations*: Throughout this article, we use uppercase and lowercase boldface letters to represent
matrices and vectors, respectively. (·)^{H} denotes the Hermitian conjugate of a matrix or a vector, and ∥·∥ denotes the spectral
norm of a matrix or the Euclidean norm of a vector. The notation
**A **≽ **B** means that **A **− **B** is a positive semidefinite matrix. We denote identity matrix by **I **while we use **0** to represent all-zero vectors or matrices. We denote
*x* is circularly symmetric complex Gaussian distributed with mean *μ* and covariance matrix **Q**.

### System model and problem formulation

Consider a *K*-user downlink MISO system where the BS is equipped with *M* antennas. The BS intends to transmit *K* independent data streams *s*_{1}*s*_{2},…,*s*_{K} to *K* users, respectively. Each *s*_{k}is circularly symmetric complex Gaussian distributed with zero mean and unit variance.
As we are considering the MISO scenario, each user is equipped with only one antenna.
The BS uses linear beamforming vector
*s*_{k}, *k*=1,2,…,*K*. The received signal at receiver *k* is

where
*k*-th user, and
*k*. Therefore, the SINR of the *k*-th user is

The BS is interested in minimizing its total transmission power (
**h**_{k}’s are known exactly at the BS (perfect CSIT), the problem can be formulated as minimizing
total transmission power subject to SINR constraints (non-robust formulation)

where *γ*_{k }> 0 is the SINR requirement for receiver *k*. Although problem (P_{NR}) is nonconvex, it can be transformed to a convex problem and solved globally and
efficiently
[1-4].

Because of the finite length of training signal and/or limited feedback bandwidth,
perfect CSIT is not available in practice. Imperfect CSIT can dramatically impact
the system performance and the users’ QoS. To maintain the desired QoS for all users,
we consider robust transmit beamforming problem under imperfect channel knowledge.
In particular, let us assume that the channel vector **h**_{k} lies within a ball with radius *ε*_{k} around the estimated channel vector

for all *k*=1,2,…,*K*, where
*ε*_{k}. In addition, we suppose that the desired QoS constraints are given in the form of
SINR constraints. Hence, the robust QoS-constrained beamforming problem can be formulated
as follows:

In the following, a SDR method is proposed for solving Problem (*P*_{R}).

### SDP relaxation and its dual

In this section, we will first present the SDR formulation of problem (*P*_{R}) and its dual problem, and then present three properties of the optimal solutions
of problem (*P*_{R}) and its dual problem.

#### SDP relaxation and its dual problem

According to the definition of
_{R}) is equivalent to

where

and

are both quadratic functions of variable *δ*_{k}. Using the *S*-Procedure
[23,24] this problem with infinitely many constraints can be reformulated as an SDP with
rank constraints.

#### Lemma 1

(S-procedure) Let

holds true if and only if there exists a nonnegative *λ*such that

Using the above *S*-procedure, we can rewrite the problem (*P*_{R}) as

The nonlinear constraint
**G**_{k} ≽ 0 and rank(**G**_{k}) = 1. Because of the rank constraint the above problem is nonconvex and difficult
to solve. If we drop this rank constraint, we obtain the following convex SDP problem
which can be efficiently solved.

We introduce the following dual variables to dualize the corresponding constraints (Table 1).

**Table 1.** Dual variables and their corresponding constraints

Then the dual of this SDP can be formulated as

Notice that the channel uncertainty bounds
*(Pε)* and its dual *(Dε)*.

#### Properties of primal problem *(Pε)* and dual problem *(Dε)*

The following three properties of the optimal solutions of the primal problem *(Pε)* and dual problem *(Dε)* will be useful throughout the rest of the article.

#### Property 1

Problem *(Dε)* is always strictly feasible.

#### Proof

See Appendix 1. □

#### Property 2

Suppose that problem *(Pε)* is feasible, then, strong duality holds true for problem *(Pε)* and its dual problem *(Dε)*, and they can both attain their optimal values.

#### Proof

Since problem *(Pε)* is feasible and by Property 1, Slater’s conditions
[23] always hold for problem *(Dε)*. Hence, the problem *(Pε)* can attain its minimum and strong duality holds
[23]. In addition, according to weak duality theorem we know that problem *(Dε)* has a finite maximum value, i.e., *α*_{k}is bounded above. This means **Y**_{k }and **y**_{k }are bounded. Hence, the feasible set of problem *(Dε)* is bounded and closed, which means that problem *(Pε)* can attain its maximum value. □

#### Property 3

Suppose {**Y**_{k},**y**_{k},*α*_{k}} is any optimal solution of Problem *(Dε)*, then

(a) *α*_{k }> 0 for all *k*=1,2,…,*K*;

(b)
*k*.

#### Proof

If there exists one *α*_{k} = 0, we could increase its value without compromising dual feasibility, which would
then contradict the optimality of *α*_{k}’s. This implies that (a) holds true. Notice that
**Y**_{k}, **y**_{k} and *α*_{k} are dual optimal solutions of dual problem *(Dε)*. Hence,

which proves (b). □

### Main results

In this section, we will show that for *ε*_{k}’s lying in a region upper bounded by some specific value, the solution of the problem
*(Pε)* must be of rank-one, which means that the SDP-relaxation *(Pε)* is tight in this case.

To this end, we first study the relationship between the robust beamforming problem
(*P*_{R}) and its non-robust version (*P*_{NR}). Obviously, if problem (*P*_{NR}) with
*P*_{NR}) is feasible.

#### Lemma 2

Let
*P*_{NR}) with

for all *k* and

Then, for any ** ε** ∈

*Ψ*, the problem (

*P*

_{R}) must have a feasible solution.

#### Proof

Based on the results of
[1], all the constraints of the problem (P_{NR}) should be active at the optimal solution, i.e.,

which is equivalent to

Using the definition of *η*_{k} (c.f. (6)), we can easily see that for any *ε*_{k} ∈[0*η*_{k}),

For each *k*=1,2,…,*K*, we further define

and

Combining the fact that *t* ≥ *t*_{k }with equality (9) and using the inequality in (8), we can show

Now for any

where the second equality is a direct consequence of (7), and the first inequality
is due to Cauchy–Schwartz inequality and the matrix norm inequality. The last inequality
follows from (10). Hence,
*P*_{R}). □

#### Remark 1

Lemma 2 shows that Problem (*P*_{R}) is feasible when uncertainty bounds are not too large. The region of uncertainty
bounds that guarantee feasibility of Problem (*P*_{R}) can be calculated from the optimal solution of Problem (*P*_{NR}).

In order to prove the results in our article we need the following key lemma. For the sake of legibility, the proof of this lemma is relegated to Appendix 2.

#### Lemma 3

Suppose that the strong duality holds for problem *(Pε)* and its dual problem *(Dε)*. Suppose *γ*_{k} ≥ 1, *k*=1,2,…,*K*. Let {*λ*_{k}}, {**G**_{k}}, {**Y**_{k}}, {**y**_{k}}, {*α*_{k}} be any primal and dual optimal solutions. Then, the optimal solution **G**_{k }of problem *(Pε)* and **Z**_{k }of problem *(Dε)* satisfy

Now, we present the main result of the article, which states that the solution of
Problem (*P*_{R}) is of rank 1 when the uncertainty bounds *ε*_{k}’s are small enough.

#### Theorem 1

Suppose that for some choice of uncertainty bounds

Then, for any vector of uncertainty bounds
*(Pε)* is feasible, and each optimal solution {**G**_{k}} must be rank one, i.e., rank(**G**_{k})=1 for all *k*.

#### Proof

For any

(13) implies that any feasible point of problem
*(Pε)*. Hence, by Property 2, the strong duality holds and the optimal values of *(Pε)* and its dual *(Dε)* are attained.

For any
*V*(** ε**) denote the optimal value of

*(Pε)*(or its dual

*(Dε)*). Then, it is obvious that

Let {*λ*_{k}}, {**G**_{k}}, {**Y**_{k}}, {**y**_{k}}, and {*α*_{k}} be any primal and dual optimal solutions satisfying complementary slackness for
the SDP problem *(Pε)* and its dual *(Dε)*. As we proved in Lemma 3 in Appendix 2, we have rank(**G**_{k}) ≥ 1. To establish rank(**G**_{k}) ≤ 1, we use the complementarity slackness condition Tr(**G**_{k}**Z**_{k})=0 and prove rank(**Z**_{k})=*M*−1. By Property 3, we have *α*_{k} > 0 and:

In addition,

Hence,

Now we can compute rank(**Z**_{k}):

where the last inequality holds true due to the fact that

#### Remark 2

Combining the result of Theorem 1 and Lemma 2, we can easily find a set of uncertainty
bounds for which solving problem *(Pε)* will find a solution to (*P*_{R}). In other words, the SDR is tight.

#### Remark 3

It is easy to extend the result of the above theorem to the case where the uncertainty bounds form an ellipsoid [20] as follows:

for all *k*=1,2,…,*K*, where **B**_{k} is a positive definite matrix,
*δ*_{k}, and *ε*_{k} are defined as before. By a similar derivation as in Theorem 1, we can prove that,
the solution of the problem *(Pε)* must be rank one when *ε*_{k}s are small enough such that

for all *k*.

### Special case: two transmit antennas

As we proved in the previous section, the SDR relaxation is tight if the error is
small enough. In this section, and the following one we will see that in the special
cases of two and three transmit antenna case, we can further improve the results.
Suppose the BS is equipped with two antennas and all the predefined SINR thresholds
*γ*_{k}’s are larger than or equal to 1, the following theorem shows that one can globally
solve the problem (*P*_{R}) by showing that all optimal solutions of *(Pε)* are rank-one solutions.

#### Theorem 2

For given positive constants *σ*_{k}, *γ*_{k}, *ε*_{k}, and channel vectors
*k*=1,2,…,*K*, we assume that problem *(Pε)* is feasible, and *γ*_{k} ≥ 1, *k*=1,2,…,*K*, then, any optimal solution **G**_{k}of problem *(Pε)* is rank one.

#### Proof

As we assume that problem *(Pε)* is feasible, by Property 2, we know strong duality holds for problem *(Pε)* and *(Dε)*, and their optimal values are both attained. Suppose {*λ*_{k}}, {**G**_{k}}, {**Y**_{k}}, {**y**_{k}} and {*α*_{k}} are any primal and dual optimal solutions with zero duality gap. Then, applying
Lemma 3, we complete the proof. □

### Special case: three transmit antennas

For the case that the BS is equipped with three antennas and all the predefined thresholds
*γ*_{k}’s of SINR are larger than or equal to 1, if {**G**_{k}} is solution of problem *(Pε)*, the following Theorem 3 shows that there are at least *K*−1**G**_{i}’s that must be of rank one.

#### Theorem 3

For given positive constants *σ*_{k}, *γ*_{k}, *ε*_{k}, and vector
*k*=1,2,…,*K*, we assume that problem *(Pε)* is feasible and *γ*_{k} ≥ 1, for all *k*. Then, for any set of optimal solutions
*(Pε)*, there are at least *K*−1**G**_{k}’s with rank one.

#### Proof

As we assume that problem *(Pε)* is feasible, by Property 2, we know strong duality holds for problem *(Pε)* and *(Dε)*, and their optimal value are both attained. Suppose {*λ*_{k}}, {**G**_{k}}, {**Y**_{k}}, {**y**_{k}} and {*α*_{k}} are any primal Using the fact that *γ*_{k} ≥ 1 and
*s*,*j* ≤ *K*, *s*≠*j*we have

which implies that

From Lemma 3, for any **Z**_{k}, we have

From (18) and (19), any two of **Z**_{k}’s can not be simultaneously of rank one. Hence, we can conclude that there is at
most one
**Z**_{k}’s are of rank two. Since we assume the primal problem is feasible and thus the strong
duality holds, then, we can conclude from the complementarity condition that there
are at least *K*−1 of **G**_{k}’s that are rank one. □

### Simulation results

In this section we first provide some numerical examples to show that the SDR method indeed gives rank one solution when the uncertainty region is small enough. In addition, the simulations indicate that when the uncertainty regions are in the form of balls, the SDR methods always gives rank one solution (as far as the problem is feasible), even if the uncertainty bound is large. This is an interesting property which can motivate future works on this problem.

In the rest of this section, we will provide some examples to show that the SDR method
does not always generate a rank one solution. These examples are generated using ellipsoid
uncertainty regions which we discussed in Remark 3. We now present simulation results
to corroborate the result of Theorem 1 and to demonstrate the effectiveness of the
SDR method. In the simulations, the BS has three antennas (i.e., *M*=3) and serves three users (i.e., *K*=3); the channel coefficients are independent and identically distributed (i.i.d)
complex Gaussian random variables with zero mean and unit variance; the channel uncertainty
bounds *ε*_{k}’s are all 0.1; the noise variances of all users are set equally to be 0.1.

To demonstrate the tightness of SDR method we compute the rank one test (ROT) ratio defined as

where *λ*_{i}(*G*_{k}) denotes the *i*-th largest eigenvalue of *G*_{k}. The simulation result is shown in Figure
1 where each data point represents the maximum ROT ratio over 2000 Monte Carlo simulations.
From this figure, we can see that the ROT ratio is very close to zero, suggesting
that the solution of the SDR method is indeed always of rank one.

**Figure 1.** **ROT ratio versus the target SINR.**

We also compare the SDR method with the convex restriction method proposed in [18]. Figure 2 presents the average transmission power versus the target SINR level, where each data point is averaged over 2000 Monte Carlo simulations. Note that, to guarantee the feasibility of the problem, we compare two methods only when the convex restriction method is feasible. From the figure, it can be observed that the SDR method is more power-efficient than the convex restriction method (which is abbreviated as ‘cvxRes’ in the figure).

**Figure 2.** **Average transmission power versus the target SINR.**

Now we present three counter examples in which the uncertain bound does not satisfy inequalities (18) and the solution is not rank one accordingly. In the examples, we consider the channel uncertainty regions as (18). We set

in (18).

#### Example 1

Three channel vectors are

with SINR requirement *γ*_{1} = 0.8174, *γ*_{2} = 0.6475 and *γ*_{3} = 0.7893. and ellipsoid norm bound

The solution of problem
**G**_{3} has eigenvalues

Our proposed rank one bound is

and Figure
3 reveals the rank distribution with respect to uncertain bound *ε*_{3} when

#### Example 2

Three channel vectors are

**Figure 3.** **ROT ratio of ****G**_{3 }**Vs. the uncertain bound ****ε**_{3}**.**

The SINR requirements are the same as previous example and ellipsoid norm bound is

and Figure
4 reveals the rank distribution with respect to uncertain bound *ε*_{3} when

**Figure 4.** **ROT ratio of ****G**_{3 }**versus the uncertain bound ***ε*_{3}**.**

#### Example 3

Three channel vectors are

with SINR requirement *γ*_{1} = 1.2174, *γ*_{2} = 0.6475 and *γ*_{3} = 0.7893 and ellipsoid norm bound

and Figure
5 illustrates the change in the rank with respect to the increase in the uncertainty
bound *ε*_{3} when

**Figure 5.** **ROT ratio of ****G**_{3 }**Vs. the uncertain bound ***ε*_{3}**.**