Weighted sum-rate maximization for multi-user SIMO multiple access channels in cognitive radio networks

In this article, an efficient distributed and parallel algorithm is proposed to maximize the sum-rate and optimize the input distribution policy for the multi-user single input multiple output multiple access channel (MU-SIMO MAC) system with concurrent access within a cognitive radio (CR) network. The single input means that every user has a single antenna and multiple output means that base station(s) has multiple antennas. The main features are: (i) the power distribution for the users is updated by using variable scale factors which effectively and efficiently maximize the objective function at each iteration; (ii) distributed and parallel computation is employed to expedite convergence of the proposed distributed algorithm; and (iii) a novel water-filling with mixed constraints is investigated, and used as a fundamental block of the proposed algorithm. Due to sufficiently exploiting the structure of the proposed model, the proposed algorithm owns fast convergence. Numerical results verify that the proposed algorithm is effective and fast convergent. Using the proposed approach, for the simulated range, the required number of iterations for convergence is two and this number is not sensitive to the increase of the number of users. This feature is quite desirable for large scale systems with dense active users. In addition, it is also worth noting that the proposed algorithm is a monotonic feasible operator to the iteration. Thus, the stop criterion for computation could be easily set up.


Introduction
The radio spectrum is a precious resource that demands for efficient utilization as the currently licensed spectrum is severely underutilized [1]. Cognitive radio (CR) [2][3][4], which adapts the radios operating characteristics to the real-time conditions, is the key technology that allows flexible, efficient and reliable spectrum utilization in wireless communications. This technology exploits the underutilized licensed spectrum of the primary user(s) (PU) and introduces secondary user(s) (SU) to operate on the spectrum that is either opportunistically being available or concurrently being shared by the PU and the SU. Under this situation and according to the definition of a cognitive (radio) network [5], opportunistically *Correspondence: phe@ee.ryerson.ca 1 Department of Electrical and Computer Engineering, Ryerson University, Ontario,M5B 2K3, Canada Full list of author information is available at the end of the article utilizing the spectrum means that the SUs may fill the spectrum gaps or holes left by the PUs; whereas concurrently utilizing the spectrum means that the SUs transmit over the same spectrum as the PUs, in a way that the interference from the transmitting SUs does not violate the quality requirement from the PUs. This article focuses on the latter case. Furthermore, the multipleinput multiple-output (MIMO) technology uses multiple antennas at either the transmitter or the receiver to significantly increase data throughput and link range without additional bandwidth or transmitted power. Thus it plays an important role in wireless communications today. In infrastructure-supported networks, such as the widely used cellular network, base stations are typically shared by a large number of users. Within the scope of this article, it is therefore assumed that the base station under consideration is shared by multiple PUs and multiple SUs. In this article, a MIMO-enhanced CR network is considered http://asp.eurasipjournals.com/content/2013/1/80 to fully ensure the quality of service (QoS) of the PUs as well as to maximize the weighted sum-rate of the SUs. We consider multiple SUs accessing the base station, referred to as a multiple-access channel (MAC).
The weighted sum-rate maximization problem is to compute the "best" achievable rate vector in the capacity region [6][7][8] by specifying the working point at the boundary of the capacity region. This optimality problem is of the Pareto meaning under multi-objective optimization.
For the non-CR cases, the sum-rate maximization problem has been intensively explored for both Gaussian broadcast channel (BC) [9,10] and Gaussian MAC [11]. Typical approaches include iterative water-filling algorithms [9,11] and dual-decomposition [10]. The conventional water-filling algorithm [12] which is an efficient resource allocation algorithm needs to be used inside each of the iterations as an inner loop operation. In addition, the setup of the well known duality between the Gaussian BC and the sum-power constrained Gaussian dual MAC [13][14][15] facilities the transform of BC sum-rate problems into its dual MAC problem. As for the weighted sum-rate problem, it is easily seen that as the weighted coefficients all being unity, the weighted sum-rate problem is reduced into a sum-rate optimization problem. Thus, solving the weighted sum-rate problem is more general. However, due to the more complicated problem structure, the conventional water-filling algorithm [12] is not able to compute its solution. For computing the maximum weighted sum-rate for a class of the Gaussian singleinput multiple-output (SIMO) BC systems or equivalent dual MAC systems [16], has presented some algorithms using a cyclic coordinate ascent algorithm to provide the max-stability policy.
For the CR cases, besides the individual power constraints to the SUs, the total interference power from the SUs needs to be included into the constraints of the target problem. Since single-antenna mobile users are quite common and compose a major served group due to the size and cost limitations of mobile terminals, this article is confined to a single input multiple output multiple access channel (SIMO-MAC) in the CR network. Earlier study [17,18] investigated the sum-rate problem and the weighted sum-rate problem in CR-SIMO-MAC cases, respectively. In addition, for the ergodic sum capacity of single input single out (SISO) system [19], studied the maximum (non-weighted) sum-rate problem with a simple form of the objective function.
In this article, by exploiting the structure of the weighted sum-rate optimization problem, we propose an efficient iterative algorithm to compute the optimal input policy and to maximize the weighted sum-rate, via solving a generalized water-filling problem in each of the iterations. The water-filling machinery is experiencing continuous development [12,[20][21][22][23]. In this article, we propose a generalized weighted water-filling algorithm (GWWFA) to form a fundamental step (inner loop algorithm) for the target problem. In the inner loop, the weighted sum-rate problem is decomposed into a series of generalized water-filling problems. With this decomposition, a decoupled system with each equation of the decoupled system containing only a scalar variable is formed and solved. Any one of the equations is solved by the GWWFA with a finite number of loops. To speed up the computation of the solution to each of the equations, we also specify the intervals the solution belongs to.
For the outer loop of the algorithm, a variable scale factor is applied to update the covariance vector of the users. The optimal scale factor is obtained by maximizing the target objective value (i.e., the weighted sum-rate) in the scalar variable β to expedite convergence of the proposed algorithm. In order to achieve this purpose, we determine an optimal scale factor by searching in a range which consists of a few discrete values. As a result, parallel operation can be used to expedite the search and to avoid the requirement of another nested loop. This parallel operation can be distributed to and carried out by multiple processors (for example, four processors).
Compared with earlier study [18], the main difference of our study is that: (i) in [18], the dual-decomposition approach [10] is used. In our study, we apply the iterative water-filling algorithm [9] and extend the algorithm to solve the target problem. The advantage of the iterative water-filling algorithm is that it is a monotonic feasible operator to the iteration. That is to say, the proposed algorithm generates a sequence composed of feasible points in its iterations. The objective function values, corresponding to this point sequence, are monotonically increasing. Hence, the stop criterion for computation might be easily set up. However, the regular primal-dual method used in [18] is not a feasible point method; (ii) for the constraints of the target problem, we make the individual power constraints more strict and more reasonable, due to the values of signal powers being assumed to be greater than or equal to zero; (iii) the convergence rate is improved significantly. In the numerical example illustrated by Figure 1 of [18], the convergence of the weighted sum-rate is obtained after 110 iterations for a system with 3 SUs and 2 PUs. However, with our proposed algorithm, we achieve the weighted sum-rate convergence with two iterations with the simulated range (number of SUs up to 110). In addition, even if the PUs and SUs are served by different base stations, it is easy to see that the proposed machinery can be used with some minor modifications.
In the remainder of this article, the system model for a CR-SIMO-MAC system and its weighted sum-rate are described in Section 2. Section 3 discusses the proposed http://asp.eurasipjournals.com/content/2013/1/80 algorithm to solve the maximal weighted sum-rate problem through an inner loop algorithm. The optimality proof of the inner loop algorithm (GWWFA) is presented in Section 3.1. Then the outer loop algorithm (AWCR) and its implementation are presented in Section 3.2. Section 4 provides the convergence proof of the AWCR. Section 5 presents numerical results and some complexity analysis to show the effectiveness of the proposed algorithm.
Key notations that are used in this article are as follows: |A| and Tr (A) give the determinant and the trace of a square matrix A, respectively; E[ X] is the expectation of the random variable X; the capital symbol I for a matrix denotes the identity matrix with the corresponding size. A square matrix B 0 means that B is a positive semi-definite matrix. Further, for two arbitrary positive semi-definite matrices B and C, the expression B C means the difference of B − C is a positive semi-definite matrix. In addition, for any complex matrix, its superscripts † and T denote the conjugate transpose and the transpose of the matrix, respectively.

SIMO-MAC in CR network and its weighted sum-rate
For a SIMO-MAC in a CR network, as shown in Figure 1, assume that there are one base-station (BS) with N r antennas, and K SUs and N PUs, each of which is equipped with one single antenna. The received signal y ∈ C N r ×1 at the BS is described as where the jth entry x j of x ∈ C K×1 is a scalar complex input signal from the jth SU and x is assumed to be a Gaussian random vector having zero mean with independent entries. The jth entryx j ofx is a scalar complex input signal from the jth PU andx is assumed to be a Gaussian random vector having zero mean with independent entries. The noise term, z ∈ C N r ×1 is an additive Gaussian noise random vector, i.e., z ∼ N(0, σ 2 I). The channel input,x, x, and z are also assumed to be independent. Furthermore, the jth SU's transmitted power can be expressed as Note that S j , ∀j, is non-negative. The mathematical model of the weighted sum-rate optimization problem for the SIMO-MAC in the CR network can be stated as follows (refer to (2.16) in [6] and the references therein): Given a group of weights {w k } K k=1 which is assumed to be in decreasing order (users can be arbitrarily renumbered to satisfy this condition) with the achievable rate of the secondary user k, the weighted sum-rate is organized by where, for the MAC cases, the peak power constraint on the kth SU exists and is denoted by a group of positive numbers: P k , k = 1, . . . , K; the power threshold to ensure the QoS of the PUs is denoted by the positive number P t . Further, when no confusion is possible, f wmac is simply written as f. For convenience, we define η k by w k − w k+1 for k = 1, . . . , K − 1; and η K by w K , as a group of nonnegative real numbers, and assume at least one of them to be non-zero. Further, the term g k = h k h † k , ∀k, is the channel power gain of the kth SU to the BS. Also, we denote the covariance matrix of the random vector N j=1ĥ † jx j + z by C 0 . It is easy to see that the matrix C 0 is positive definite.
The constraint K k=1 g k S k ≤ P t is called the sum-power constraint with gains. The constraint is obtained in the following analysis. Let Thus, the received signal at BS is y =Ĥx + (Hx + z), where Hx + z can be regarded as the additive interference and noise to the transmitted signal Hx from the PUs. To guarantee the QoS for the PUs, the power of the interference and noise should be less than a threshold value, P TH . This condition can be expressed as It can be written as where the power constraint value P t is the interference and noise threshold subtracted by the Gaussian noise power.
As an alternative, to guarantee the QoS for each of the PUs, individually, the power of the interference and noise should be less than a threshold value, P TH (i), ∀i. Similarly, it is obtained that That is to say, the condition above is equivalent to Name min i {P t (i)} as P t ; then the target model can still cover the case that the QoS for each of the PUs is considered individually. Note that at the base station with multiple antennas, the received signals can be regarded as a stochastic vector or point in a Hilbert space and the received signal powers are abstracted into the norm squared of the vector. The transmitted powers of the PUs have been taken into account by forming C 0 and P t mentioned above, which appear in (3).
It is seen that the sequence {η k } K k=1 stems from the vector of weights used in the multi-user information theory [6]. The parameter or item η k , ∀k, in the sequence is called the weighted coefficient without confusion.
A more strict weighted sum-rate model can also be obtained that reflects the essence of the issue for the CR-SIMO-MAC. Along a similar way mentioned above, we may choose the power thresholds P t,i to limit the impact from the SUs on each of the antennas of the BS. Thus the sum-power constraint with the gains is evolved into It is seen that such a weighted sum-rate problem with more power constraints can be solved by solving a similar problem in (3). Therefore, the proposed article aims at computing the solution to the problem (3). Note that if (3), we remove the user i 0 and then the number of the users is reduced to K − 1. Thus, we can assume that For the important special case of the sum-rate problem, which is included in (3), assume that M = rank(H). Applying the QR decomposition, H = QR, let Q = [q 1 , . . . , q M ] ∈ C N r ×M have orthogonal and normalized column vectors. R ∈ C M×K is an upper triangle matrix with r i,j denoting the (i, j)th entry of the matrix R. Q † is regarded as an equalizer to the received signal by the BS. Thus, the ith SU should have the rate: It is easy to see that the rate just mentioned comes from the expression: i.e., C 0 = I in this case.

Algorithm AWCR
The proposed algorithm for solving the weighted sumrate problem in the cognitive radio network, denoted by AWCR, is an iterative algorithm. It consists of two layers of loops. Inside the inner loop, a generalized weighted water-filling algorithm is proposed and used. Due to special problem structure and the complexity of the weighted sum-rate problem, the proposed water-filling is more general than regular weighted water-filling. It is discussed in Section 3.1. For the outer loop of AWCR, a variable scale factor with parallel computation is applied to expedite the convergence. This discussion is presented in Section 3.2. http://asp.eurasipjournals.com/content/2013/1/80

Generalized weighted water-filling algorithm (GWWFA)
Being a fundamental block of the optimum resource allocation problem for the CR-SIMO-MAC systems, the generalized water-filling problem is abstracted as follows.
For a multiple receiving antenna system, it is given that P t > 0, as the total power or volume of the water; K is the total number of the users; the allocated power and the propagation path (non-negative) gains for the ith user are given as S i for i = 1 . . . K, and {a ij } K j=i , respectively. The generalized weighted water-filling problem under consideration then reads where the set Due to the specific CR SIMO MAC setup considered as well as the inclusion of arbitrary weights {η j }, the problem structure (9) is novel. It is easy to see that if a ij = 0, as i = j, and P i >> 0, ∀i, then the problem (9) is reduced into the conventional weighted water-filling problem. Further, if equal weights are chosen, it is reduced into the conventional water-filling problem, which can be solved by the conventional water-filling algorithm [12].
To find the solution to the more complicated generalized problem above, the generalized weighted waterfilling algorithm (GWWFA) is presented as follows. Let Utilize a permutation operation π on {λ i } such that It is easy to see that the function J i (s i ) is strictly monotonically decreasing and continuous over the interval The steps of the GWWFA can be described as below.
Then increase the iteration from n to n + 1. Repeat the procedure in (14) until the point s (n) Otherwise, go to step (2). Remarks 3.1. Note, in (1) of the GWWFA, that the initial λ min may be chosen as λ π(K+1) , and λ max may be chosen as λ π(1) .
In (3), for the initialization of s (0) π(k) , first, we may choose an interval, such as [ 0, P π(k) ] , and use the secant method or the bisection method [24] over the interval to compute, in parallel, an approximate solution to the system J π(k) (s π(k) ) − λ = 0, ∀k. Hence, only through a few loops (≤ log 2 P π(k) + 1 loops), |e 0 |, as an absolute error between the accurate solution and the approximate solution obtained by the method above, is less than 0.5. The initialization of s (0) π(k) , for k = 1, . . . , i, is assigned by the above approximate solution.
Since the problem in Proposition 3.1 is a differentiable convex optimization problem with linear constraints, not only is the KKT condition mentioned above sufficient, but it is also necessary for optimality. Note that it is easily seen that the constraint qualification (the CQ) of the optimization problem (9) holds. Proposition 3.1 hence is proved.

Remarks 3.2.
To decouple the variables in the objective function of the problem (3), a sum expression is acquired by adding the objective function, just mentioned, K times. Then the sum expression is operated, by one variable being selected as an optimized variable with respect to the others being fixed. Thus, from the expression (3), the problem (20), is implied as follows: Since where S k , ∀k, is fixed and the optimization problem is equivalent to the problem below: If the CR SIMO weighted case is generalized to the CR MIMO weighted case, it is still an open question whether there exists a fast water-filling solution like the algorithm mentioned above.
(1) Generate effective channels where the superscript with a pair of bracket, (n), represents the number of iterations.
(2) Treating these effective channels as parallel, noninterfering channels, the new covariances are generated by the GWWFA under the sum power constraint P t . That is to say, is the optimal solution to (24): Note that (24) is similar to (20), only S (n−1) i and G (n) ij in the former take place of S i and G il in the latter, respectively, for any i, j, l. Let as the innovation, where the function f has been defined in (3). Then, the covariance update step is The updated covariance is a convex combination of the newly obtained covariance and the immediate past covariance.
(4) Increase the iteration from n to n + 1. Go to (1) until convergence.
Note that the new algorithm employs variable weighting factors, which are obtained to maximize the objective function and to update the covariance.
In this section, the optimality of S (n) is the solution to (20), by Proposition 3.1.

Remarks 3.3.
Due to the objective function f βγ (n) + (1 − β)p (n−1) in step (3) of Algorithm AWCR being (upper) convex, i.e., being concave, in the scalar variable β, for computing the maximum solution to the corresponding optimization problem, we can choose finite searching steps with even fewer evaluations of the objective function. Without loss of generality, the objective function in step (3) is evaluated at the four points β = 1 K , 1 K + 1 3 1 − 1 K , 1 K + 2 3 1 − 1 K and 1 by parallel computation to determine β * . That is to say, this parallel operation can be distributed to and carried out by multiple processors (for example, four processors) at the base station, in order to expedite convergence of the proposed algorithm. Finally, the obtained satisfying solution is then distributed or returned to the corresponding secondary users.

Convergence of Algorithm AWCR
There are two methods by either of which convergence of the proposed algorithm can be proved. The first method is to utilize convergence of Algorithm AWCR with β * = 1 4 (refer to [25]) and the innovation, as a spacer step, by Zangwill's convergence theorem B ( [26], p. 128). However, we will then still need to prove Algorithm AWCR with β * = 1 4 , as a basic mapping, to satisfy the closedness condition of Zangwill's convergence theorem B. This point requires much explanation and an abstract proof. As an alternative, the second method which is more intuitive than the first method is used. The fixed point approach proposed in this article could also be generalized to solve other problems.
In this section, utilizing results from Section 3.1, convergence of Algorithm AWCR will be strictly proved under a weaker assumption. Note that, due to the power constraint being coupled between the optimization stages of (20) with the weighted coefficients in the objective function while being decoupled between the optimization stages of the MIMO-MAC case without the weighted coefficients, usage of the waterfilling principle in the former is different from that of the latter.

Convergence proof of the proposed algorithm
In this article, as a more general model, we eliminate the assumption in [9] that the optimal solution is unique to prove convergence of the proposed algorithm. To the best knowledge of the authors, this is one of the proposed novelty for convergence of this class of algorithm http://asp.eurasipjournals.com/content/2013/1/80 with the spacer step ( [26], p. 125). Since our convergence proof is based on more general functions including an objective function and a few constraint functions, it will also enrich the optimization theory and methods. It is assumed that a mapping projects a point to a set. First, two concepts are introduced. The first concept is of an image of a mapping (or algorithm) that projects a point to a set; the second one is of a fixed point under the mapping (algorithm). Then, two lemmas are proposed, followed by the convergence proof of the proposed algorithm.

Definition 4.1. (Image under mapping or Algorithm
A) (see e.g., [26], p. 84). Assume that X and Y are two sets. Let A be a mapping or an algorithm from X to Y, which projects from a point in X to a set of points in Y. If the point in X is denoted by x and the set of the points in Y is denoted by A(x), then A(x) is called the image of x under A.

Definition 4.2. (Fixed point under mapping or Algorithm A). Let A be a mapping or an algorithm from
Note that (20) can be changed into a general form: due to the condition of the optimal solution uniqueness being removed. Further, corresponding to this change, step (2) of Algorithm AWCR will be carried out in this way: given a feasible point , its image under step (2) of Algorithm AWCR is a set of points. A point in this set is chosen arbitrarily as the next point S generated by Algorithm AWCR. Thus, Algorithm AWCR can generate a point sequence under this change. In the following, we will still call this algorithm Algorithm AWCR despite the changes. The feasible set is denoted by V d .
For any convergent subsequence, whose limit is denoted by (S 1 , . . . , S K ), generated by Algorithm AWCR, we may use the following lemma to prove that the limit is a fixed point under Algorithm AWCR, when Algorithm AWCR is regarded as a mapping. Proof. See Appendix 2.
Based on the lemmas above, we obtain the conclusion that Algorithm AWCR is convergent. At the same time, step (3) of Algorithm AWCR is then regarded as a computation for a point. With these lines of proofs, Algorithm AWCR generates a point sequence and every point of the point sequence consists of the K non-negative numbers, e.g., S (n) 1 , . . . , S (n) K . The details are described below. Proof. Due to compactness of the set of feasible solutions for the problem in (3), the point sequence generated by Algorithm AWCR already includes a convergent subsequence. For every convergent subsequence, according to Lemma 1, the convergent subsequence must converge to a fixed point under Algorithm AWCR. Then, according to Lemma 2, it converges to one of the optimal solutions to the problem in (3).
In addition, conversely, as stated by the sufficient and necessary conditions of Lemmas 1 and 2, for any optimal solution to the problem in (3), there is a point sequence generated by Algorithm AWCR such that the point sequence converges to that optimal solution.
With Algorithm AWCR generating the point sequence, the definition of Algorithm AWCR and (30) in Appendix 1 imply that the sequence of the objective values, obtained by evaluating the objective function at the point sequence, monotonically increases to the optimal objective value. This is due to (30) and any convergent subsequence of the point sequence converging to one of the optimal solutions to the optimization problem in (3).
To reduce the cost of computation, (20) and (25) in Section 3 may utilize the Fibonacci search. To improve the performance of the algorithm and reduce the cost of the computation, the objective function in step (3) of the AWCR can be evaluated at the four points mentioned in Remark 3.3, by parallel computation to find the estimate of β * of (25). http://asp.eurasipjournals.com/content/2013/1/80

Numerical results and complexity analysis
In this section, numerical examples are provided to illustrate the effectiveness of the proposed algorithm. For comparison purpose, a regular feasible direction method utilizing the gradient [27] in the optimization is chosen. It is denoted as Algorithm AFD. Note that, as a benchmark and a feasible direction method, the Algorithm AFD can also generate a sequence of feasible points (as a feasible point algorithm). It is easy to set up a stop criterion of computation for a feasible point algorithm, especially for a monotonic feasible point algorithm like the proposed one. Due to the feasible set being a convex polygon, the recently developed AFD algorithm is used as a reference. We didn't select [18] for comparison since the primal-dual algorithm used in [18] is not a feasible point method; in addition, the assumption of the constrains is different and system model is different, too. Figures 2 and 3 show the evolution of the weighted sum-rate values versus the number of iterations for AWCR and AFD for some choices of the number of users (K). In the calculation, the number of antennas at the base station (m) is set to be 4. Channel gain vectors are generated randomly using random m × 1 vectors with each entry drawn independently from the standard Gaussian distribution. {P k } is the set of randomly chosen positive numbers. The sum power constraint is P t = 10 dB. A group of different weights are also generated randomly. In these figures, the cross markers and the diamond markers represent the results of our proposed Algorithms AWCR and AFD, respectively. These results show that the proposed Algorithm AWCR exhibits much faster convergence rate, especially with an increasing number of users.
Let f * be the maximum sum-rate, f (n) the sumrate at the nth iteration and |f (n) − f * | the error in the sum-rate. Figures 4 and 5 show the corresponding error in the sum-rate versus the number of iterations. Note it is easy to see that using the fixed-point theory of the proposed Lemma 2 one can determine the maximum sum-rate f * mentioned. As shown in these figures, the algorithms converge linearly. The proposed algorithm exhibits a much larger slope in the sum-rate error function, which translates to a faster convergence rate.
We can further observe that the convergence rate of the proposed algorithm is not sensitive to the increase of the number of users. For clearly understanding, we define where the point {(j, f (j) )} is generated by the AWCR and = 10 −3 without loss of generality; N AFD is similarly defined but generated by the AFD. Each of these numbers can be regarded as the required number of iterations for the corresponding algorithm. We simulate different selection of K, and list the corresponding N AWCR and N AFD in Table 1. We can observe that in the simulated range, using the proposed algorithm, the required number of iterations for convergence is about 2, whereas for the AFD, the required number of iterations is much larger.
Since the AFD and the proposed algorithm use the same matrix inverse operations, which consist of the most significant part of the computation, to compute the gradient of the objective values, both algorithms have similar computational complexity O(m 3 ) in each of the iterations (refer to [28]). This is because for a m × m square matrix, its inverse needs m(m 2 − 1) + m(m − 1) 2 , i.e., O(m 3 ), arithmetic operations; its determinant needs 2 3 m 3 + m, i.e., O(m 3 ), operations (the Cholesky decomposition approach is used for efficiency and our objects). Thus, since these operations are used with finite times, it is easily seen that, for each iteration, computational complexities for both AFD and AWCR are O(m 3 ).
Also for conveniently checking the algorithms, deterministic instances are chosen as η k = k K k=1 k , ∀k, P t = 10 dB and P i = 9 dB, ∀i, and the channel gains are randomly generated as for K = 4 and K = 5, respectively. Let the normalized covariances of C 0 be the identity matrix. The calculated weighted sum-rate is plotted as a function of the iterations in Figure 6. It is shown that N AWCR = 1, keeping the same least value for both cases; and N AFD = 10 and 12 as K = 4 and K = 5, respectively.

Conclusion
The proposed algorithm AWCR, as a class of iterative water-filling algorithms, is used to solve the problem of the weighted sum-rate for the MIMO-MAC in a CR network. By exploiting the concept of variable weighting factor for covariance update, together with the machinery of distributed and parallel computation, the proposed AWCR algorithm can greatly speed up the convergence rate of the weighted sum-rate maximization computation. The required number of iterations for convergence exhibits non-sensitivity to the increase of the number of the users. Furthermore, a novel GWWFA, as a fundamental block of the proposed algorithm, is proposed.
Convergence of the proposed algorithm is strictly proved by the designed fixed point theory. We present an equivalent optimality condition by Lemma 2, i.e., a point is one of the optimal solutions to the problem of maximum weighted sum-rate for the MISO-MAC in the CR network if and only if the point is a fixed point of the AWCR. In the derivation, for more general problems, the assumption used in [9] that the optimal solution is unique to prove the convergence could be eliminated. Numerical examples are presented to demonstrate the effectiveness of the proposed algorithm. In the simulated range, the required number of iterations for convergence is shown to be fixed at two, which is a significant reduction compared with the conventional algorithms.

Proof of Lemma 1
Note that in the following proof, we use the notation n to stand for the number of iterations for convenience.
The necessity is proved first. For the limit (S 1 , . . . , S K ) of any convergent subsequence, there is a convergent sub- for any n and (S 1 , . . . , S K ) ∈ V d . Replacing n with n k , we obtain: We have the following relationships:  Note that the set arg max (S 1 ,...,S K )∈V d K i=1 f (S 1 , . . . , S i−1 , S i , S i+1 , . . . , S K ) does not need to be a singlepoint set. However, we may choose (S 1 , . . . , S K ) as one of the optimal solutions to the problem max (S 1 ,...,S K )∈V d K i=1 f (S 1 , . . . , S i−1 , S i , S i+1 , . . . , S K ). This corresponds to step (2) of Algorithm AWCR. Further, (S 1 , . . . , S K ) = β * (S 1 , . . . , S K ) + (1 − β * )(S 1 , . . . , S K ), based on the choice of the optimal solution mentioned above. This corresponds to step (3) of Algorithm AWCR. Therefore, resulting from the two correspondences mentioned above and the definition of Algorithm AWCR, it is true that (S 1 , . . . , S K ) is a fixed point under Algorithm AWCR, which is viewed as a mapping.
The sufficiency will be proved as follows: If (S 1 , . . . , S K ) is a fixed point under Algorithm AWCR, it is seen that if (S Note that in the proving process above, we do not have the following assumption:

Proof of Lemma 2
The necessity is proved first.