 Research
 Open Access
 Published:
Dualbased bounds for resource allocation in zeroforcing beamforming OFDMASDMA systems
EURASIP Journal on Wireless Communications and Networking volume 2013, Article number: 51 (2013)
Abstract
We consider multiantenna base stations using orthogonal frequencydivision multiple access and space division multiple access techniques to serve singleantenna users. Some users, called realtime users, have minimum rate requirements and must be served in the current time slot while others, called non realtime users, do not have strict timing constraints and are served on a besteffort basis. The resource allocation (RA) problem is to find the assignment of users to subcarriers and the transmit beamforming vectors that maximize the total user rates subject to power and minimum rate constraints. In general, this is a nonlinear and nonconvex program and the zeroforcing technique used here makes it integer as well, exact optimal solutions cannot be computed in reasonable time for realistic cases. For this reason, we present a technique to compute both upper and lower bounds and show that these are quite close for some realistic cases. First, we formulate the dual problem whose optimum provides an upper bound to all feasible solutions. We then use a simple method to get a primalfeasible point starting from the dual optimal solution, which is a lower bound on the primal optimal solution. Numerical results for several cases show that the two bounds are close so that the dual method can be used to benchmark any heuristic used to solve this problem. As an example, we provide numerical results showing the performance gap of the wellknown weight adjustment method and show that there is considerable room for improvement.
1 Introduction
With the ubiquitous use of smart phones, tablets, laptops and other devices, traffic demand on wireless access networks is increasing geometrically. Multiantenna base stations using orthogonal frequencydivision multiple access (OFDMA) and space division multiple access (SDMA) can transmit at the same time to different sets of users on multiple subcarriers. This diversity increases the system throughput by assigning transmitting resources to users with good channel conditions. High data rates are thus made possible by exploiting the degrees of freedom of the system in time, frequency and space dimensions. OFDMASDMA is also supported by WiMAX and LTEAdvanced systems, which are the technologies currently being set up to implement the fourth generation (4G) cellular networks [1, 2].
Due to these increased degrees of freedom, it is critical to use a dynamic and efficient scheduling and resource allocation (RA) mechanism that takes full advantage of all OFDMASDMA transmitting resources [3]. In this mechanism, the scheduler selects the users that are served at each frame and the RA algorithm allocates to these users the transmission resources required to meet the quality of service (QoS) requested from the upper layers.
In this article, we focus on the RA part of the problem for an OFDMASDMA system supporting real time traffic with minimum rate requirements. We propose a dualbased method to get an upper bound to all feasible solutions. Because the dual solution is not necessarily primal feasible, we also present an algorithm to get a feasible point from the dual solution, which is a lower bound to the optimal primal solution. The importance of these bounds is that they give us limits on the duality gap and we can use them to estimate how far the solution given by any heuristic method is from the optimum.
1.1 State of the art
The general RA is a nonlinear, nonconvex program so that it is almost impossible to solve it directly for any realistic number of subcarriers, users and antennas. For this reason, most research study focuses on developing heuristic algorithms. In this context, an important question is always how accurate are the results. For the RA problem with rate constraints, it turns out that there are very few results of that kind, as we shall see.
Traffic in the system can be divided into two main groups: Delaysensitive real time (RT) services and delayinsensitive non real time (nRT) services. Early studies on OFDMA systems focused on solving the RA problem for nRT services only, where the objective is to maximize the total throughput with only power constraints and possibly minimum BER constraints. In [4], the complete optimization problem is divided into two subproblems: Selection of users for each carrier and then power allocation to these users, which are both solved by a heuristic. A similar approach using zeroforcing (ZF) beamforming is reported in [5]. The study of [4, 5] does not solve the complete optimization problem; instead, it separates it into uncoupled subproblems that provide suboptimal solutions.
There is a definite need to benchmark the performance of these heuristic algorithms. For the RA problem with nRT traffic only, several methods have been proposed to compute nearoptimal solutions. For example, genetic algorithms are proposed in [6], while [7–10] provide methods to compute a nearoptimal solution based on dual decomposition. In addition to providing a benchmark, nearoptimal algorithms can also lead to the design of efficient RA methods as shown in [10], where heuristic algorithms derived from the dual decomposition method are proposed.
Several methods have been used to solve the RA problem for OFDMASDMA systems with both RT and nRT traffic. In [11], the objective is to maximize the sum of the user rates subject to peruser minimum rate constraints that model the priority assigned to each user at each frame. The optimization problem is solved approximately for each frame by minimizing a cost function representing the increase in power needed when increasing the number of users or the modulation order. The advantages of this approach are that it handles user scheduling and RA together and supports RT and nRT traffic. Its weaknesses are that no comparison is made against a nearoptimal solution and the method used to determine user priorities at every frame is very complex.
In [12], both RT and nRT traffic are supported. Priorities are set according to the remaining deadline time for RT users and to the difference between the achieved rate and the desired rate required for nRT users. The user with the highest priority is paired with the subchannel with the highest vector norm and semiorthogonal users are multiplexed on the same subchannel. Comparisons against the algorithm of [11] show that the packet drop rate for RT users is significantly reduced. The algorithm’s complexity is also reduced because of the semiorthogonality criteria used to add users. However, as in [11], a performance comparison with a nearoptimal solution is not provided.
In [13], the objective is to maximize a utility function without any hard minimum rate constraints for the RT users. The channel quality information is added to the utility function to favor users with good channel conditions and priorities are set by increasing user weights in the utility function. The advantage is that the perframe optimization problem has only a power constraint and no rate constraints, which makes it simpler to solve. The disadvantage with this reactive approach is that RT users with poor channel conditions are backlogged until their delay is close to the deadline, increasing the average delay and jitter.
In [14], a heuristic algorithm is proposed for the sum rate maximization problem with proportional rates among the user data rates, i.e., the ratio among allocated user rates is predetermined. The criteria to form user groups includes semiorthogonality as in [12], but also fairness through proportional rate constraints. This method is extended to include hard minimum rates in [15] which attempts to solve exactly the same problem we formulate in Section 2. Again, there is not reported method to evaluate the accuracy of these heuristics, except by comparing them to each other.
Another approach is to maximize the sum rate subject to constraints on the average rate delivered to a user [16]. However, unlike the study presented in [17], where an optimal solution is provided for the single antenna OFDMA RA problem with average rate constraints, the algorithm presented in [16] is a heuristic approximation. Note also that with average rate constraints, RT users tend to be served when they have good channel conditions which can create unwanted delay violations and jitter.
1.2 Article contribution and organization
None of the previous study provide a nearoptimal solution to the OFDMASDMA RA problem with minimum rate requirements. This is important not only as a benchmark for any heuristic algorithms, but also to get a better insight into the problem and to help devise efficient heuristics.
The main contribution of this article is a method that provides an upper bound to the following OFDMASDMA RA problem for mixed RT and nRT traffic: for a given time slot, find the user selection and beamforming vectors that maximize a linear function of the users rates, given a total transmit power constraint and minimum rate constraints for RT users. The user weights in the linear utility function are arbitrary and can be the result of a prioritization or fairness policy by the scheduler. We solve the RA problem by Lagrange decomposition and we show, for small cases where we can find the optimal solution, that the duality gap is small.
A second contribution is a simple offline heuristic algorithm to compute a feasible point based on the dual solution. This point is a lower bound for the optimal solution. This lower bound, in conjunction with the upper bound from the dual, can be used to limit the optimality gap in larger cases where an optimal solution is not available.
We then study several cases where we compare the performance of the upper and lower bounds. The results show that the two bounds are tight when the number of RT users is small and that their difference increases for larger values but that it stays quite small. Thus, the dual method provides a good approximation to the optimal solution. We also compare the performance of the weight adjustment algorithm versus the solution provided by the two bounds. The results indicate that adjusting the user weights to prioritize RT users can lead to significantly suboptimal solutions.
We describe the system and formulate the optimization problem we want to solve in Section 2, where we also briefly discuss a direct enumeration method to find the optimal solution for small problems. We present in Section 3 the dual method and two algorithms: One that finds the dual solution (the upper bound) and the other that finds a dualbased primal feasible solution (the lower bound). In Section 4, we present numerical results showing the accuracy of the upper and lower bounds and of the weight adjustment algorithm for different scenarios. Finally, we present our conclusions in Section 5.
2 System description and problem formulation
We consider the RA problem for the downlink transmission in a MISO system with a single base station. There are K users, some of which have RT traffic with minimum rate requirements while the others have nRT traffic that can be served on a besteffort basis. The BS is equipped with M transmit antennas and each user has one receive antenna. In this configuration, the BS can transmit in the downlink data to different users on each subcarrier by performing linear transmit beamforming precoding. At each OFDM symbol, the BS can change the beamforming vector for each user on each subcarrier to maximize some performance function. In this article, we assume that we use a channel coding that reaches the channel capacity. The data rate are in units of bits per OFDM symbol, or equivalently bits per second per Hertz (bps/Hz).
2.1 Signal model
First, we describe the model used to compute the bit rate received by each user. Define
${\stackrel{~}{s}}_{k,n}$
the symbol transmitted by the BS to user k on subcarrier n. We assume that the ${\stackrel{~}{s}}_{k,n}$ are independently identically distributed (i.i.d) random variables with ${\stackrel{~}{s}}_{k,n}\sim \mathcal{C}\mathcal{N}(0,1)$.
w_{k,n} an Mcomponent column vector representing the beamforming vectors for user k on subcarrier n. Unless otherwise noted, we denote w the vector made up by the column stacking of the vectors w_{k,n}.
x_{ n } an Mcomponent column vector representing the signal sent by the array of M antennas at the BS for each subcarrier n.
h_{k,n} an Mcomponent row vector representing the channel between the M antennas at the BS and the receive antenna at user k for each subcarrier n.
z_{k,n} the additive white gaussian noise at the receiver for user k on subcarrier n. The z_{k,n} are i.i.d. Gaussian random variables and, without loss of generality, we assume that they have unit variance, that is ${z}_{k,n}\sim \mathcal{C}\mathcal{N}(0,1)$.
y_{k,n} the signal received by user k on subcarrier n.
${r}_{k,n}^{0}$the rate of user k on subcarrier n in bps/Hz.
The signal vector x_{ n } is built by a linear precoding scheme which is a linear transformation of the information symbols ${\stackrel{~}{s}}_{k,n}$:
The signal received by user k on subcarrier n is then given by
The second and third terms in the righthand side of (2) correspond to the interference and noise terms. Since the signals and noise are Gaussian, their sum is also Gaussian and the data rate of user k for subcarrier n is given by the Shannon channel capacity for an additive white Gaussian noise channel:
2.2 Rate maximization problem
The general rate maximization problem corresponding to the OFDMASDMA RA problem with mixed RT and nRT traffic is to find a set of beamforming vectors w_{k,n} that will maximize the weighted sum of user rates. This is limited by the total power available for the transmission at the base station and some users with real time QoS requirements must receive a minimum rate. More precisely, we assume that we know
K Number of users in the cell.
$\mathcal{K}$ Set of users in the cell: {1, …, K}.
$\mathcal{D}$ A subset of $\mathcal{K}$ containing the users that have minimum rate requirements. We define $D=\left\mathcal{D}\right$.
${\stackrel{\u02c7}{d}}_{k}$ Minimum rate requirement for user K.
M Number of antennas at the BS.
N Number of subcarriers available.
$\stackrel{\u030c}{P}$Total power available at the base station for transmitting over all channels.
c_{ k } Weight of the user rates in the objective function. These could be computed by the scheduler to implement prioritization or fairness.
We then want to solve the following optimization problem to obtain the RA
The total transmit power is represented by the sum of the L_{2} norms of the beamforming vectors in constraint (5). The achievable rate over all subcarriers should be higher or equal than the required minimum rate per user as in (6).
Problem (4–6) is a nonconvex, nonlinear optimization problem. Using an exact algorithm to find a global optimal solution is very hard considering the size of a typical problem where there can be up to a hundred users and hundreds of subchannels. Another option is to use a standard nonlinear program (NLP) solver to compute a local optimal solution and use different starting points in the hope of finding a good global solution. The problem with this approach is that (1) we don’t know how close we are to the true optimum and (2) the technique is quite timeconsuming.
Nevertheless, we tried this approach for some small cases and observed that most users end up with a zero beamforming vector and only a small subset of users actually get some rate, often no more than M. Furthermore, in accordance to what was reported for the SDMA problem in [18], we observed that at high SNR, a socalled zeroforcing (ZF) solution is very close to the local optimum. This ZF solution can be easily computed by channel diagonalization and waterfilling power allocation and is nearoptimal compared to the general beamforming solution in the high SNR regime. Moreover, due to multiuser diversity, there is a good chance of finding users with high SNR channels when the total number of users increase. Therefore, in a multiuser Rayleigh fading environment, the ZF beamforming technique can provide results close to the general beamforming solution, even in the moderate SNR regime. For these reasons, we now turn to the socalled ZeroForcing beamforming strategy.
2.3 Zeroforcing beamforming
In general, user k is subject to the interference from other users which reduces its bit rate, as indicated by the denominator in (3). Zeroforcing beamforming is a strategy that completely eliminates interference from other users. For each subcarrier n, we choose a set s of g ≤ M users which are allowed to transmit. This is called an SDMA set. We then impose the condition that for each user k in this set, the beamforming vector of user k must be orthogonal to the channel vectors of all the other users of the set. This amounts to adding the orthogonality constraints
and the user k data rate for subcarrier n then simplifies to:
With ZF beamforming, the problem is now made up of two parts. We need to select a SDMA set s(n) for each subcarrier n and, for each selected SDMA set, we must compute the beamforming vectors in such a way that the total rate received by all users is maximized. Because of this, we need to add another set of decision variables
α_{k,n} a binary variable that is 1, if we allow user k to transmit on subcarrier n and zero otherwise. We denote the collection of α_{k,n} by the vector α.
This results in the following ZF problem
where A and B are some large positive constants. Constraint (12) guarantees that we do not choose more than M users for each subcarrier, constraint (13) guarantees that if we have chosen two users k and j, they meet the ZF constraints and is redundant for other users, and constraint (14) guarantees that the beamforming vector is zero for users that are not chosen. Problem (9–14) is a nonlinear mixed integer program (MIP) and these are known to be very hard to solve exactly.
In this article, whenever we need to get an exact solution, we use a complete enumeration of the binary variables α satisfying the constraints (12). For each α, we then find the beamforming vectors maximizing the utility function given the power, minimum rate and ZF constraints. The optimal solution is the vector α that maximizes the utility function. The problem of finding the optimal w variables for a particular α is a relatively small nonconvex problem. Using the pseudoinverse approach to satisfy the ZF constraint (see Section 3.3), it can also be transformed into an approximate convex problem and solved using a dual algorithm similar to the one described in Section 3.4.
It would seem that the zeroforcing model is not improving things much: We have gone from a nonconvex nonlinear program to a nonconvex mixed nonlinear program that can be solved by picking the global solution from a large collection of small convex problems. However, as we explain in Section 3, this allows us to design an efficient and accurate algorithm.
3 Dualbased solution method
Obviously, we cannot solve Problem (9–15) fast enough to use it in a realtime system. Not only is it NPcomplete [4] but the actual computation time becomes quickly prohibitive for realistic sizes, even for offline computations. Still, we need to compute solutions so that we can use them as benchmarks to evaluate the quality of realtime heuristic approximations. We now present two offline solution techniques that are tractable for problems of realistic size based on the Lagrange relaxation of the primal. Because Problem (9–15) is not convex, there will often be a strictly positive duality gap at the solution of the dual. However, if it is small enough, we can use the solution provided by the dual method as a useful benchmark to check the accuracy of heuristic methods. Results in Section 4 show that in many cases, the duality gap is in fact less than a few percent.
Solving the ZF problem will require some form of search over the α_{k,n} variables. Note that this ranges over all subsets with a number of users smaller than or equal to M, so that the search space is going to be fairly large. Our first transformation is thus to separate the problem into singlesubcarrier subproblems. For this, we dualize the constraints (10) and (11) since they are the ones that couple the subcarriers. Define the dual variables
λ Lagrange multiplier for power constraint (10).
μ_{ k } Lagrange multipliers for minimum rate constraint (11) of user K. The collection of μ_{ k } is denoted μ.
In order to simplify the derivation, we also define the dual variables μ_{ k } for all users $k\in \mathcal{K}$. For users with no minimum rate requirements ($k\notin \mathcal{D}$), we have μ_{ k } = 0 and ${\stackrel{\u02c7}{d}}_{k}=0$. In what follows, we use the standard form of Lagrangian duality which is expressed in terms of minimization with inequality constraints of the form ≤. Under these conditions, the multipliers λ, μ ≥ 0. We get the partial Lagrangian
with constraints (12–15). The value of the dual function Θ at some point (λ, μ) is obtained by minimizing the Lagrange function over the primal variables
and the dual problem is
which we can solve by the well known subgradient algorithm [19]. From now on, we concentrate on the calculation of the subproblem (17).
3.1 Subchannel subproblem
Because of the relaxation of the carriers coupling constraints (10–11), the subproblems in (17) decouple by subcarrier since the objective (16) is separable in n and so are constraints (12–14). Computing the dual function then requires the solution of N independent subproblems. For each subcarrier n, this has the form
where w_{ n } is the vector made up by the column stacking of the vectors w_{k,n} for subcarriers n and α_{ n } denote the collection of α_{k,n} for subcarrier n. Problem (20–23) is still a mixed NLP, albeit of a smaller size.
3.2 SDMA subproblem
A simple solution procedure is to enumerate all possible choices for α_{k,n} that meet constraint (12). This is called the extensive formulation of the problem. Each choice defines a SDMA set s and κ = s and here, we present a solution technique for the subproblems with given s. For each s, the problem separates into κ independent problems to compute the optimal beamforming vector w_{k,n,s} for each user k ∈ s. For each user k ∈ s we know the set of channel vectors for the other members of s and we collect these vectors in the (κ  1) × M matrix H_{k,n,s}. Problem (20–23) can then be rewritten as
where ${f}_{k,n,s}^{\ast}$ is given by the solution of the optimization
where ${c}_{k}^{\prime}={c}_{k}+{\mu}_{k}$, w_{k,n,s} is the beamforming vector for user k on subcarrier n for SDMA set s, and w_{n,s} is the vector made up by the column stacking of the vectors w_{k,n,s} for the κ users in s. Note that constraint (21) is automatically satisfied by the construction of s, constraint (23) simply drops out since w_{k,n,s} = 0 for k ∉ s and constraint (22) remains only for k ∈ s, but we write it as (27) because we are considering only users that belong to SDMA sets.
This is certainly not a feasible realtime algorithm, but for realistic values of K and M the number of SDMA sets is still manageable and the optimization subproblem (24–27) is a small nonlinear program with M variables and κ  1 linear constraints. It can thus be solved quickly by a number of techniques. Still, the overall computation load can be quite large. There will be κ such problems to solve for each SDMA set, and there are $S=\sum _{i=1}^{M}\left(\genfrac{}{}{0.0pt}{}{K}{i}\right)$ such sets for each of the N subcarriers so that we have to solve the problem κ × S × N times, and this for each iteration of the subgradient algorithm. Clearly, any simplification of the beamforming subproblem can reduce the overall computation time significantly.
3.3 Approximate solution to the beamforming problem
This can be done by the following construction. Instead of searching in the whole orthogonal subspace of H_{k,n,s} as defined by (27), we pick a direction vector in that subspace and search only on its support. This will give a good approximation to the extent that the direction vector is close to the optimal vector. The choice of direction is motivated by the fact that the objective function depends only on the product h_{ k }w_{k,n,s}. We then introduce a new independent variable
and because this variable is not independent of w_{k,n,s}, we add (28) as a constraint. We then get the equivalent problem
Constraints (30) and (31) can then be rewritten in the standard form G_{k,n,s}w_{k,n} = b_{k,n,s} where the G_{k,n,s} matrix is the concatenation of h_{ k } and H_{k,n,s} and b_{k,n,s} = [q_{k,n,s}, 0, 0 … 0]^{T}.
Since we are proposing to transform the constrained optimization over the κ variables into an unconstrained optimization over q_{k,n,s} only, we must be able to express w_{k,n,s} as a function of q_{k,n,s}. The linear system being underdetermined, this is obviously not unique. We propose to use ${\mathbf{G}}_{k,n,s}^{+},$ the pseudoinverse of G_{k,n,s}, for the backtransformation ${\mathbf{w}}_{k,n,s}={\mathbf{G}}_{k,n,s}^{+}{\mathbf{b}}_{k,n,s}.$ The pseudoinverse picks the vector of minimum norm compatible with the linear system. In other words, choosing this transformation will minimize ∥w_{k,n,s}∥ so that it is minimizing the power term in the objective function in (29). Because λ ≥ 0, this has the effect of contributing to the maximization of ${f}_{k,n,s}^{\ast}$. Note that this technique provides only an approximate solution of the beamforming problem; we cannot invoke Theorem 1 from [20] which shows that in certain cases, the pseudoinverse transformation is optimal. A strong assumption for the theorem is that the objective function depends only on the q_{k,n,s} variable, which is not the case here since (29) also depends on ∥w_{k,n,s}∥^{2}. However, we observed from numerical results that the difference between the pseudoinverse solution and the optimal solution is small. With this approximation we fix the direction of the beamforming vectors to
where ${\left[{\mathbf{G}}_{k,n,s}^{+}\right]}_{1}$ denotes the first column of ${\mathbf{G}}_{k,n,s}^{+}$. Now, we can obtain a problem formulation in terms of the user powers only by replacing the following expression in (29):
where ${\gamma}_{k,n,s}=\parallel {\left[{\mathbf{G}}_{k,n,s}^{+}\right]}_{1}\parallel $ and p_{k,n,s} = q_{k,n,s}^{2}. Adding the constraint p_{k,n,s} ≥ 0, we get the equivalent problem
which has the solution
so that the computation time is basically that for computing ${\mathbf{G}}_{k,n,s}^{+}$. Also note that using ${\mathbf{G}}_{k,n,s}^{+}$ we can also find the optimal beamforming vectors for all users in s, the only difference being that γ_{k,n,s} is computed using the column of ${\mathbf{G}}_{k,n,s}^{+}$ corresponding to the channel vector of this user.
3.4 Solving the dual problem
To summarize, the dual function Θ(λ, μ) is obtained for the current values of the multipliers by finding for each subcarrier $n=1,\dots ,N$ the optimal SDMA set s^{∗}(n) to the minimization problem in (24), where
and p_{k,n,s} is given by (35). Substituting back in (17), the dual function is
with f_{n,s} given by (36). For any value of the dual variables (λ, μ) we can determine the optimal primal variables (α, w); α is obtained by the optimal subcarrier assignment vector s(n) after performing the minimization over s in (37), and the optimal beamforming vectors ${\mathbf{w}}_{n,k}^{\ast}$ for the users k ∈ s^{∗}(n) are given by
The largest part of the computation to evaluate the dual function is the calculation of ${\mathbf{G}}_{k,n,s}^{+}$ which has to be done for each subchannel and each possible SDMA set. The number of evaluations can become quite large but the size of each matrix is relatively small so that the calculation remains feasible for mediumsize networks. Furthermore, while solving the dual problem requires multiple subgradient iterations, the calculation of the pseudoinverses is independent of the value of the multipliers. This means that the calculation of ${\mathbf{G}}_{k,n,s}^{+}$ can be done only once in the initialization step of the subgradient procedure.
Algorithm 1. Calculation of the dual problem
Note that if the algorithm finds an optimal solution, the corresponding primal computed from the optimization of the Lagrange function will be feasible since the subgradients ${g}_{\mu}^{\left(k\right)}$ and g_{ λ } are small. We use the (somewhat inappropriate) expression “dual feasible” to denote such a solution. If, on the other hand, the algorithm stops before reaching the optimum, the primal will generally not be feasible.
Algorithm 1 finds the optimal dual variables (λ^{∗}, μ^{∗}) that solve the dual problem (18) using the subgradient method [19] with a fixed step size δ and therefore yields an upper bound to problem (9–15). In Algorithm 1, the pseudoinverse matrices ${G}_{k,n,s}^{+}$ can be computed by any number of well known algorithms. Here, we have used the Matlab command pinv[21]. Note that this algorithm can be used to solve the beamforming problem or equivalently the power allocation problem for a fixed SDMA set assignment. The only difference is that (24) becomes trivial since there is only one possible SDMA set per subcarrier.
It is not possible to evaluate the overall complexity of Algorithm 1 because we don’t have a bound on the number of dual iterations. We can nevertheless give a bound on the complexity of one iteration since it is dominated by the pseudoinverse matrices computation. We have to compute N S ∼ N K^{M} matrix inversions each one with complexity O(M^{3}), giving a total computational complexity O(N K^{M}M^{3}). Computing the subgradient vectors and updating the dual variables have much lower complexity than the pseudoinverse matrices computation.
3.5 Analysis of the dual function
Lets denote Θ^{∗} the maximum of the dual function Θ(λ, μ) over (λ, μ) ≥ 0. If U^{1} is the weighted sumrate objective function achieved by any feasible point in the primal problem and U^{∗} its optimum value, then the following inequalities hold [22]
The value Θ^{∗}, or any feasible approximation Θ(λ, μ) to it, is thus an upper bound to the optimum value of the primal problem U^{∗}. Thus, Θ^{∗} can be used to benchmark any other solution method.
Figure 1 shows a contour plot of the dual function Θ(λ, μ) for the case of one RT user. The diamond marker shows the maximum.
We can get some insight on the shape of the dual function from Figure 2, where we plotted the function with respect to μ for a fixed value of λ. The solid line curve corresponds to the same dual function as in Figure 1, where the rate constraint is active, ${\stackrel{\u02c7}{d}}_{1}=50$ bps/Hz. We see that the dual function goes through a maximum at μ = 0.24. We have also shown the case where we increase the minimum rate constraint so much that the problem becomes infeasible, e.g., we make ${\stackrel{\u02c7}{d}}_{1}=100$ bps/Hz. As expected from duality theory, the dual function has no maximum since lim_{μ→∞}Θ(λ, μ) = ∞ as shown by the dashdotted curve. Finally, the dashed curve at the bottom corresponds to ${\stackrel{\u02c7}{d}}_{1}=0$ bps/Hz such that the constraint is inactive and the solution where the maximum occurs is located at μ_{1} = 0.
For completeness we show in Figure 3 the dual function as a function of λ when the rate constraints are feasible. The dual function increases rapidly and reaches a maximum at λ = 1.83.
3.6 Dualbased primal feasible method
The SDMA set selection and beamforming vectors found by Algorithm 1 do not always provide a primal feasible solution. The rate or power constraints might be violated whenever the algorithm stops because the number of iterations has been reached before the convergence rule is met. In this subsection, we propose a simple procedure to obtain a feasible point to problem (9–15) from the dual solution found with Algorithm 1. This point is not optimal but because we start from the dual optimal solution, we expect that it will be close to the optimal solution. Obviously this will give us a lower bound to the optimal primal solution.
Algorithm 2. Calculating a feasible point from the dual solution
Algorithm 2 summarizes this method. The algorithm begins by solving the dual problem (18–19) using Algorithm 1. If the solution is not feasible either directly or by recomputing the power allocation for the SDMA set assignment found in the dual problem, the algorithm performs a search by increasing the dual variables associated to the users whose QoS constraints are not met until a new SDMA set assignment is found. It then solves the power allocation problem for this new SDMA set assignment and checks if the minimum rate constraints are met. The search for new SDMA sets continues using this method until a feasible SDMA set assignment is found or a maximum number of iterations is reached.
Algorithm 2 invokes Algorithm 1 which has complexity O(N K^{M}M^{3}). Assuming that the maximum number of iterations is fixed independently of the problem parameters, the computational complexity of the feasible point search is O(N K^{M}) for the subcarrier reassignment and O(N(D+K)) for the power allocation. These are lower than the computational complexity of Algorithm 1. Therefore, the computational complexity of Algorithm 2 is also O(N K^{M}M^{3}).
In contrast to the enumeration method described in Section 2.3 which performs an enumeration of all possible SDMA set assignments, the dualbased Algorithm 2 is a method that finds new candidate SDMA set assignments close to the dual optimal and then uses them to solve a simple power allocation problem until the rate and power constraints are met. This makes the search for a nearoptimal feasible point much faster than finding the exact solution.
4 Performance analysis
In this section, we present some numerical results on the performance of the dualbased algorithm and the accuracy of the upper and lower bounds. To show how they can be used to evaluate heuristic algorithms, we also compare them with the solution provided by a weight adjustment method described in Section 4.2.
4.1 Convergence of the dual algorithm
First we present in Figure 4 the value of the dual function and Lagrange multipliers as a function of the number of iterations for a given channel realization. The corresponding transmit power and the rate received by the RT user are shown in Figure 5. The parameters used for the calculation are listed in the figure titles.
We see that the algorithm converges very quickly to a solution that is both close to the minimum value and feasible. This is typical of several other configurations, except that the number of iteration increases with the number of RT users.
4.2 Weight adjustment heuristic
As discussed in Section 1.1, several RA algorithms provide support for users with RT traffic by increasing the user weights in the utility function until they receive enough transmission resources. In this section, we describe a generic weight adjustment method which will be used to show that this technique leaves much room for improvement.
Algorithm 3. Weight adjustment algorithm
We want to find a set of weights in the utility function (9) such that the rate requirements of the RT users are met when we solve problem (9–15) without the rate constraints (11). Also, the set of weights must not be very different from one user to the other in order to maximize the multiuser diversity gain. Algorithm 3 implements a generic method that can do this. It increases the user weights for RT users until enough resources are allocated to meet the minimum rate requirements. The parameter ϵ controls how much the weights are increased with respect to the rate bounds. The rates achieved by Algorithms 2 and 3 are different since they solve different problems. Algorithm 3 can be seen as solving problem (9–15) by using a linear penalty method for constraints (11) of the form
The modified objective function is then
At each iteration of the penalty method, whenever rate constraints are active, the solution of (41) cannot be smaller than that of (9–15) since it is a relaxation. Notice that problem (41) is quite simple since it has a single constraint for the transmit power but it has to be solved many times to adjust the weights of the real time users. In weight adjustment algorithms such as [13], the user weights are increased at each time slot using an increasing function of the packets delay, so the computation task is distributed over time. However, this distributed approach does not guarantee that the rate requirements are met in a given time slot which can lead to delay violations and jitter.
4.3 Parameter setup and methodology
We now present the method and parameter values used to compare the performance of the different methods to solve Problem (9). We used a Rayleigh fading model to generate the user channels such that each component of the channel vectors h_{k,n} are i.i.d. random variables distributed as $\mathcal{C}\mathcal{N}(0,1)$. We also assumed independent fading between users, antennas and subcarriers. Unless otherwise noted, we used a configuration with M = 3 antennas, K = 16 users and N = 16 subcarriers. We have only one RT user when we examine the effect of various parameters and we also look at the impact of increasing the number of RT users. The minimum rate constraint was set at 40 bps/Hz unless otherwise stated. We also fixed the power constraint to $\stackrel{\u030c}{P}=20$ dBm and used a largescale attenuation of 0 dB for all users. The user weights in (9) were set to c_{ k } = 1 for all users. The results are the average over the feasible cases from 100 independent channel realizations.
We compared the performance of the different methods for various scenarios where we increased the resource requirements for the RT users until the minimum rate requirements can no longer be met for all RT users. For each scenario and channel realization, the upper bound was computed from the dual solution using Algorithm 1 described in Section 3.4. For small systems, we also found the exact solution using the primal enumeration method given in Section 2.3. We also computed the lower bound given by dualbased primal feasible Algorithm 2 and the heuristic solution provided by the weight adjustment Algorithm 3 described in Sections 3.6 and 4.2. We use the upper bound given by the dual optimal solution as the reference point when computing the gap when the exact solution is not available.
4.4 Single user, increasing minimum rate
In this first scenario, we have a single RT user and we increase its minimum rate ${\stackrel{\u02c7}{d}}_{1}$. First we consider a small system with K = 4 users and N = 2 subcarriers where it is possible to compute an exact solution. We present in Table 1 the average gap in percent between the three methods used to find feasible solutions against the dual upper bound for a small system. As the required minimum rate increases from 13.33 to 20 bps/Hz, the upper bound decreases as more resources need to be assigned to the RT user until the problem is no longer feasible. For this small configuration, we see that all methods give excellent results and the duality gap is very small.
In the remaining results, we use a larger system with K = 16 users and N = 16 subcarriers. With these values, it is no longer possible to use the primal enumeration method and we compute the gap relative to the dual upper bound. We present in Table 2 the difference in percentage between the bound and the solutions of the dual feasible and the weight adjustment algorithms. The dual feasible algorithm gives a lower bound within 0.25 % of the dual upper bound while the weight adjustment solutions difference can be almost 10 %. As discussed in Section 4.2, this is due to the fact that the weight adjustment algorithm stops as soon as it finds a feasible solution and does not have the option of finding a better assignment. As a result, the objective does not change much when the minimum rate is increased. This can be seen from Figure 6 which shows the sum rate achieved by the dual feasible algorithm and the weight modification method against the minimum rate requirement.
4.5 Single user, increasing attenuation
As the user moves away from the BS and the channel attenuation increases, the RA algorithm dedicates more resources to the RT user until the problem is unfeasible. Figure 7 shows the average total rate when the largescale channel attenuation of the RT user varies from 0 to 15 dB. The results show that for all SNR, the dual lower bound provides a tight solution with the upper bound while the weight adjustment method shows a large performance gap. Table 3 shows the error in percentage between the objective and the upper bound. For an attenuation of 15 dB, neither method is able to find a feasible solution; the problem is feasible because the dual upper bound is around 140, but the algorithms cannot find a solution. That is, Algorithm 1 for the dual upper bound finds a solution as long as the primal problem is feasible. On the other hand, Algorithm 2 cannot find a solution when the solution point is close to the infeasibility region. In that case, the algorithm exits declaring that a feasible solution could not be found. Also, note that the proposed dualbased algorithms provides a nearoptimal solution to the ZF Problem (9–14) for all SNRs where a feasible solution can be found by the algorithm.
Figure 7 shows average sum rates. The average is performed over the channel realizations that provide feasible points. When the attenuation of the RT user is low, most of the channel realizations produce feasible points. When it is high, some of the channel realizations do not produce feasible points and are discarded. For example, in Figure 7 when the attenuation is 15 dB, none or very few of the channel realizations produce feasible points. Therefore, Figure 7 does not show a sum rate for that point.
4.6 Increasing number of RT users
Finally, Figure 8 shows the upper dual bound, the lower bound and the solution given by weight adjustment methods as a function of the number of RT users. Table 4 lists the performance gap against the dual bound in percentage. The dual feasible lower bound is again very close to the upper bound. Meanwhile, we can see that the performance of the weight adjustment method quickly degrades when the number of RT users increases. It cannot find feasible points with 6 or 7 RT users while the dual algorithm yields solutions for these values within 3.52 % of the upper bound.
For a single RT user, we have seen in Tables 2 and 3 that the difference between the upper and lower solution is small. In Figure 8, we see that this difference increases for three or more RT users. Still, this growth is not large and we can consider that a 3.52 % is an acceptable error tolerance. Based on this, we can claim that it is possible to find a quasioptimal solution to problem (9–15) with the proposed method, albeit with an offline algorithm.
Furthermore, the results show that the weight modification method has a large performance gap which becomes more significant as the number of RT users increase. Also, the dual method can find feasible solutions for cases where the weight adjustment method cannot. This shows that the weight modification method should be used carefully for RA in OFDMASDMA systems with RT users and that more efficient heuristics should be developed to approach the performance of the dualbased feasible solution.
5 Conclusion
In this study, we proposed a method to compute the beamforming vectors and the user selection in an OFDMASDMA MISO system with minimum rate requirements for some RT users. We gave a rigorous mathematical formulation of the zeroforcing model and then used a Lagrangian relaxation of the power and rate constraints to solve the dual problem using a subgradient algorithm. The Lagrange decomposition yields subproblems separated per subcarrier, SDMA sets and users which substantially reduces the computational complexity. We obtained a simple expression of the dual function for the beamforming problem for a given SDMA set based on a pseudoinverse condition on the beamforming vectors. The dual optimum can then be used as a benchmark to compare against any other solution methods and heuristics. The dual function also gives us a better understanding of the problem. Its shape is related to the rate constraint activation and problem feasibility, and it also justifies the splitting of the subcarrier assignment and power allocation processes used in several heuristic methods.
We then proposed an algorithm which finds a feasible point by starting from the dualbased optimum solution and searching among the dual variables of the rate constraints. Numerical results indicate that the two bounds are tight so that the feasible point is nearoptimal solution. As a point of comparison, we also evaluated the performance of a weight adjustment method which uses weight adjustments in the objective function to achieve the required rates. Our results show that the performance gap of this approach is large and grows when the SNR of a single RT user increases or when the number of RT users increases.
In addition, the weight adjustment method requires many time slots to adjust the weights and schedule real time users. The dualbased method explicitly includes the minimum rate constraints which allows RT users to be scheduled in the current slot, which decreases the overall packet delay and jitter.
The significant gap between the weight adjustment algorithm and the optimal RA solution, suggests that there is a need to find better heuristics. The dual approach looks promising to guide the design of efficient novel heuristics. To implement the RA algorithm in real time, we also need to design fast methods to reduce the number of SDMA sets to be searched. The design of these heuristic algorithms is outside the scope of this article and is part of our current efforts. Finally, the upper and lower bounds also provide a very useful benchmark to compare the performance of any heuristic method.
Abbreviations
 BER:

Bit error rate
 BS:

Base station
 CPU:

Central processing unit
 LTE:

Long term evolution
 RA:

Resource allocation
 MIP:

Mixed integer program
 MISO:

Multiple input single output
 NLP:

Non linear program
 NP:

Nondeterministic polynomial time
 nRT:

non real time
 OFDMA:

Orthogonal frequency division multiplexing access
 QoS:

Quality of service
 4G:

Fourth generation
 RT:

Real time
 SDMA:

Spatial division multiple access
 SISO:

Single input single output
 SNR:

Signal to noise ratio.
References
 1.
3GPP: Further advancements for EUTRA: physical layer aspects Rel. 9. 2010.
 2.
IEEE: Draft amendment to IEEE standard for local and metropolitan area networks part 16: air interface for fixed and mobile broadband wireless access systems: Multihop relay specification. 2009.
 3.
Letaief K, Zhang Y: Dynamic multiuser resource allocation and adaptation for wireless systems. IEEE Wirel. Commun. Mag 2006, 13(4):3847. 10.1109/MWC.2006.1678164
 4.
Bartolomé D, PérezNeira A: Practical implementation of bit loading schemes for multiantenna multiuser wireless OFDM systems. IEEE Trans. Commun 2007, 55(8):15771587.
 5.
Maciel TF, Klein A: A resource allocation strategy for SDMA/OFDMA systems. In Proc. of 16th IST Mobile and Wireless Communications Summit. Budapest; 2007:429433.
 6.
Ozbek B, D Le Ruyet: Adaptive resource allocation for SDMAOFDMA systems with genetic algorithm. In 6th International Symposium on Wireless Communication Systems, ISWCS’09. Tuscany; 2009:438442.
 7.
Tsang Y, Cheng R: Optimal resource allocation in SDMA/multiinputsingleoutput/OFDM systems under QoS and power constraints. In Proc. IEEE Wireless Communications and Networking Conference, WCNC. Atlanta; 2004:15951600.
 8.
Chan P, Cheng R: Capacity maximization for zeroforcing MIMOOFDMA downlink systems with multiuser diversity. IEEE Trans. Wirel. Commun 2007, 6(5):18801889.
 9.
Xingmin L, Hui T, Qiaoyun S, Lihua L: Utility based scheduling for downlink OFDMA/SDMA systems with multimedia traffic. In Proc. IEEE Wireless Communications and Networking Conference, WCNC. Sydney; 2010:130134.
 10.
PereaVega D, Frigon J, Girard A: Nearoptimal and efficient heuristic algorithms for resource allocation in MISOOFDM systems. In IEEE International Conference on Communications ICC. Cape Town; 2010:16.
 11.
Tsai C, Chang C, Ren F, Yen C: Adaptive radio resource allocation for downlink OFDMA/SDMA systems with multimedia traffic. IEEE Trans. Wirel. Commun 2008, 7(5):17341743.
 12.
Chung W, Wang L, Chang C: A lowcomplexity beamformingbased scheduling for downlink OFDMA/SDMA systems with multimedia traffic. In Proc. of Global Telecommunications conference IEEE GLOBECOM. Honolulu; 2009:15.
 13.
Huang W, Sun K, Bo T: A new weighted proportional fair scheduling algorithm for SDMA/OFDMA systems. In Proc. 3rd Int. Conf. on Communications and Networking in China, ChinaCom. Hangzhou; 2008:538541.
 14.
Papoutsis SKV, Fraimis I: User selection and resource allocation algorithm with fairness in MISOOFDMA. IEEE Commun. Lett 2010, 14(5):411413.
 15.
Papoutsis V, Kotsopoulos S: Resource allocation algorithm for MISOOFDMA systems with QoS provisioning. In Seventh International Conference on Wireless and Mobile Communications, ICWMC. Luxembourg; 2011:711.
 16.
Tralli V, Henarejos P, PerezNeira A: A low complexity scheduler for multiuser MIMOOFDMA systems with heterogeneous traffic. In Proc. International Conference on Information Networking, ICOIN. (Kuala Lumpur; 2011:251256.
 17.
Wang X, Giannakis G: Ergodic capacity and average rateguaranteed scheduling for wireless multiuser OFDM systems. In International Symposium on Information Theory ISIT. Toronto; 2008:16911695.
 18.
Yoo T, Goldsmith A: On the optimality of multiantenna broadcast scheduling using zeroforcing beamforming. IEEE J. Sel. Areas Commun 2006, 24(3):528541.
 19.
Bertsekas D: Convex Analysis and Optimization. Belmont: Athena Scientific; 2003.
 20.
Wiesel A, Eldar Y, Shamai S: Optimal generalized inverses for zero forcing precoding. In Proc. of 41st Annual Conf. on Information Sciences and Systems. Baltimore; 2007:130134.
 21.
Matlab 2012.http://www.mathworks.com/help/matlab/ref/pinv.html MoorePenrose pseudoinverse of matrix. Matlab documentation center. [Online]. Available:
 22.
Boyd S, Vandenberghe L: Convex Optimization. Cambridge: Cambridge University Press; 2004.
Acknowledgements
This research project was partially supported by NSERC grant CRDPJ 33593406.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
PereaVega, D., Girard, A. & Frigon, JF. Dualbased bounds for resource allocation in zeroforcing beamforming OFDMASDMA systems. J Wireless Com Network 2013, 51 (2013). https://doi.org/10.1186/16871499201351
Received:
Accepted:
Published:
Keywords
 Resource allocation
 OFDMA/SDMA
 QoS
 Minimum rate constraints
 Zeroforce beamforming