Power allocation and precoding matrix computation method in a wireless communication system

Abstract

This invention presents methods for power allocation for each data stream on each transmitting antennas in MU-MIMO wireless communication systems comprising the BS computing the temporary beamforming matrix, constructing special matrices and vectors with the maximum transmitting power on each antenna, the allocated power to each data stream, the channel quality of each data stream, and the amplitude of each element of the temporary beamforming matrix, calculating the power allocated to each data stream on each antenna with the calculated special matrices and vectors, and adjusting each element of the temporary beamforming matrix to obtain the final beamforming matrix.

Description

FIELD OF THE INVENTION

The disclosed inventions relate generally to wireless communication, and in particular, to the mechanism for a Base Station (BS) to allocate power and precode the signal before it is transmitted to the User Equipment (UE) in massive Multiple-Input Multiple-Output (MIMO) communication systems.

BACKGROUND

Large-scale MIMO systems or massive MIMO systems were firstly introduced in [1] where each BS is equipped with dozens to several hundreds of transmit antennas. One main advantage of such systems is the potential capability to offer linear capacity growth without increasing power or bandwidth [1][4], which is realized by employing Multi-User MIMO (MU-MIMO) to achieve the significant higher spatial multiplexing gains than conventional systems. In this system, the BS groups UEs at each scheduling slot and transmits data to them on the same time and frequency resource.

It has been proved that Zero-Forcing (ZF) precoding with a total transmitting power constraint is almost the best choice to maximize the sum rate for large-scale MU-MIMO systems [2]. However, in practice, the power of each antenna is restricted instead of the total power. It means that maximizing the power utilization requires the sum power of all users at each antenna to be the same. Unfortunately, it is generally not the case in practice, because of the randomness of the ZF precoding matrix. As a result, it causes a dilemma to the BS: on the one hand, ZF precoding could not fully use the transmit power, which leads to throughput loss; on the other hand, full power utilization means that there exists residual interference among the grouped users, since the ZF precoding matrix is violated, which also results in throughput loss. Conjugate Beamforming (CB) is another practical precoding method for MU-MIMO precoding in large-scale MIMO communication systems because of its simplicity for implementation. Similarly to ZF, CB also faces the optimal power allocation problem when the power of each antenna is restricted. Therefore, more sophisticated power allocation methods are needed to maximum the sum rate of MU-MIMO systems. Due to the aforementioned reasons, this invention provides four different methods to allocate the power to each data stream based on two different optimization objectives when ZF precoding is employed by the BS. In addition, a simple power allocation method is also offered when CB is employed by the BS. The advantages of this invention include: 1. when ZF precoding is employed, two of the four power allocation methods have better performance than the rest two in the low Signal-to-Noise Ratio (SNR) region and vice versa, so different power allocation methods could be employed in different SNR regions to achieve the maximum sum rate of MU-MIMO systems; 2. when CB is employed, a very simple power allocation method could be employed with little sum rate loss; 3. the sum rate losses of all methods provided in this invention are negligible compared to the case where the total transmitting power instead of the per-antenna power is constrained; 4. most importantly, these methods are not affected by a scaling factor of each channel vector so channel estimation with an arbitrary scaling factor would be sufficient, which alleviates the accuracy requirement of channel measurement in massive MIMO systems.

SUMMARY OF THE INVENTION

This application provides several methods to complete power allocation and precoding matrix computation in MU-MIMO systems. For ZF precoding, two methods are provided to maximize the power utilization, where one is based on orthogonal projection while the other one is based on iterative searching the optimal solution in the constraint domain. In addition, two more methods are provided to minimize the inter-user interference among UEs in the MU-MIMO user group, where one is based on linear scaling while the other one is based on iterative searching. Among the four methods, the first two methods are the better choices in the low SNR region while the latter two methods are the better choices in the high SNR region. For CB, a simple but rate-lossless method is provided where the power of each antenna could be totally consumed.

BRIEF DESCRIPTION OF THE DRAWINGS

The aforementioned implementation of the invention as well as additional implementations would be more clearly understood as a result of the following detailed description of the various aspects of the invention when taken in conjunction with the drawings. Like reference numerals refer to corresponding parts throughout the several views of the drawings.

FIG. 1 shows the block diagram illustrating the iteration process of searching the optimal solution to (12).

FIG. 2 shows the block diagram illustrating the iteration process of searching the optimal solution to (19).

FIG. 3 shows a block diagram illustrating the process of power allocation and precoding matrix computation.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

For a large-scale MIMO system, suppose that the number of antennas at the BS side is M, and the number of multiplexed data streams is K in the downlink on a Resource Element (RE) such as a subcarrier, or a Resource Block (RB), etc. Note that the K data streams belong to N users, where N≤K. Suppose that the channel matrix corresponding to the K data streams is H=[h₁ h₂ . . . h_K]^T, which may be acquired by BS through uplink channel measurement or uplink feedback channel.

The K modulated signals on the current RE are precoded by a matrix W before being further processed, where W has the dimension of M×K and the (m,k)th element is w_mk, where m=1, 2, . . . , M, and k=1, 2, . . . , K.

If the ZF precoding method is employed, the BS firstly computes a temporary matrixH^inv=H^H(HH^H)⁻¹,  (1)then the (m,k)th element of the matrix H^inv can be written ash_mk^inv=|h_mk^inv|e^jθmk.  (2)

Let P₁^ant, P₂^ant, . . . , P_M^ant denote the power allocated to the current RE belonging to the M antennas respectively, then the power allocated to the kth data stream by the BS is P_k, which satisfies Σ_k=1^K P_k=Σ_m=1^M P_m^ant. One possible example is

$\begin{array}{cc}{P}_k=\frac{\sum _{m=1}^M{P}_m^{\mathrm{ant}}}{K}.& \left(3\right)\end{array}$

In order to complete the power allocation, four methods belonging to two categories are provided in this invention, where the first category is based on maximizing the power utilization, while the second one is based on minimizing the inter-user interference.

Category-1: Maximizing the Power Utilization.

In this category, the BS constructs an (M+K−1)×MK matrix A by deleting the k_dth row vector of the following temporary matrix

$\begin{array}{cc}{A}_T=\left[\begin{array}{ccccccccccccc}1& 1& \dots & 1& 0& 0& \dots & 0& \dots & 0& 0& \dots & 0\\ 0& 0& \dots & 0& 1& 1& \dots & 1& \dots & 0& 0& \dots & 0\\ ⋮& ⋮& \ddots & ⋮& ⋮& ⋮& \ddots & ⋮& \ddots & ⋮& ⋮& \ddots & ⋮\\ 0& 0& \dots & 0& 0& 0& \dots & 0& \dots & 1& 1& \dots & 1\\ 1& 0& \dots & 0& 1& 0& \dots & 0& \dots & 1& 0& \dots & 0\\ 0& 1& \dots & 0& 0& 1& \dots & 0& \dots & 0& 1& \dots & 0\\ ⋮& ⋮& \ddots & ⋮& ⋮& ⋮& \ddots & ⋮& \ddots & ⋮& ⋮& \ddots & ⋮\\ 0& 0& \dots & 1& 0& 0& \dots & 1& \dots & 0& 0& \dots & 1\end{array}\right],& \left(4\right)\end{array}$ where the (i,j)th element of A_T satisfies the conditions

$\begin{array}{cc}{a}_{\mathrm{ij}}=\{\begin{array}{cc}1,& \mathrm{if}1\le i\le K,j=\left(i-1\right)M+l,l=1,\dots ,M,\\ 0,& \mathrm{if}1\le i\le K,j\ne \left(i-1\right)M+l,l=1,\dots ,M,\\ 1,& \mathrm{if}K+1\le i\le K+M,j=\mathrm{lM}+\left(i-K\right),l=0,\dots ,K-1,\\ 0,& \mathrm{if}K+1\le i\le K+M,j\ne \mathrm{lM}+\left(i-K\right),l=0,\dots ,K-1\end{array},& \left(5\right)\end{array}$ and k_d∈{1, . . . , MK} may be any one of the MK possible values.

The BS constructs an (M+K−1)×1 vector b by deleting the k_dth element of the temporary vectorb_T=[P₁. . . P_KP₁^ant. . . P_M^ant]^T.  (6)

An MK×1 vector r is constructed asr=[r₁₁. . . r_M1r₁₂. . . r_M2. . . r_1K. . . r_MK]^T  (7)where the elements of r satisfy

$\begin{array}{cc}\frac{r_{\mathrm{mk}}}{r_{\mathrm{lk}}}=\frac{{h_{\mathrm{mk}}^{\mathrm{inv}}}^2}{{h_{\mathrm{lk}}^{\mathrm{inv}}}^2},l,m=1,\dots ,M,k=1,\dots ,K,e.g.,r_{\mathrm{mk}}={h_{\mathrm{mk}}^{\mathrm{inv}}}^2.& \left(8\right)\end{array}$

With the matrix A and vectors b and r, two possible methods could be used to compute the power allocated to each data stream on each antenna.

Method-1: Orthogonal Projection.

In this method, the BS projects the vector r into the solution space of the equation Ax=b firstly by{tilde over (p)}=[I−Ã^T(ÃÃ^T)⁻¹Ã]{tilde over (r)},  (9)where Ã=[A b] and {tilde over (r)}=[r^T−a]^T with a being a positive real number. Then, the elements of power allocation matrix are computed as

$\begin{array}{cc}{p}_{\mathrm{mk}}=\frac{\tilde{p}\left(\mathrm{mk}\right)}{\tilde{p}\left(\mathrm{MK}+1\right)},m=1,\dots ,M,k=1,\dots ,K.& \left(10\right)\end{array}$

Method-2: Iterative Searching.

In this method, the power allocation vector is computed by solving the following problemmin∥p−r∥₂²,s·t·Ap−b=0,−p≤0.  (11)

The problem (11) can be solved by iterative searching in the constraint domain. One possible solution is to firstly transform (11) into an equivalent problemmin f(t,p)=−t∥p−r∥₂²−Σ_i=1^MKp_i,s·t·Ap=b,  (12)then the iterative searching process in FIG. 1 is used to find the optimal p. Specifically, after the searching process starts 1, the input parameters are initialized as t>0, ε_o>0, ε_i>0, p₀∈{Ap=b, p_i≥0}, μ>0, and γ>0, where ε_o and ε_i are endurable errors for the outer and inner searching cycles respectively, P₀ is an initial power allocation matrix, while μ and γ are two adjusting parameters 2. Then, the outer searching cycle runs while

$-\frac{1}{t}\sum _{i=1}^{\mathrm{MK}}{p}_i>{\varepsilon }_o$ 3, and the inner searching cycle runs while λ(p)/2>ε_i 4. Inside the inner cycle, the first step is to solve the equation

$\left[\begin{array}{cc}{\nabla }^2{f}_p& A^T\\ A& 0\end{array}\right]\left[\begin{array}{c}\Delta p_t\\ \chi \end{array}\right]=\left[\begin{array}{c}-\nabla f_p\\ 0\end{array}\right],$ where χ is an adjusting parameter 5. Then, the Newton decrement is calculated as λ(p)=Δp_t^T∇²f_pΔp_t 6. After that, the variable p is updated as p←p+γΔp_t 7. After the inner cycle ends 8, the parameter t is updated as t←μt 9. After the outer cycle ends 10, the whole process ends 11. Finally, the vector p is reshaped to a matrix with elements p_mk, m=1, . . . , M, and k=1, . . . , K.

With p_mk, the elements w_mk of the precoding matrix W can be computed asw_mk=√{square root over (p_mk)}e^jθmk,m=1, . . . ,M,k=1, . . . ,K.  (13)

Category-2: Minimizing the Inter-User Interference.

In this category, two possible methods are provided to minimize the inter-user interference.

Method-1: Linear Scaling.

In this method, the BS computes a temporary power allocation matrix with elements

${\tilde{p}}_{\mathrm{mk}}=\frac{P_k{h_{\mathrm{mk}}^{\mathrm{inv}}}^2}{\sum _{m=1}^M{h_{\mathrm{mk}}^{\mathrm{inv}}}^2}$ firstly, where h_mk^inv is the same as in (2), then computes the power consumed on each antenna as Q_m=Σ_k=1^K {tilde over (p)}_mk, m=1, . . . , M. After that, the BS chooses the maximum value of Q_m, which is denoted as Q_max Finally, the power allocation matrix is computed as

$\begin{array}{cc}{p}_{\mathrm{mk}}=\frac{{\tilde{p}}_{\mathrm{mk}}{P}_m^{\mathrm{ant}}}{Q_{\max }},m=1,\dots ,M,k=1,\dots ,K.& \left(14\right)\end{array}$

Method-2: Iterative Water-Filling Method.

In this method, the BS constructs an (M+K)×K matrix A with a form of

$\begin{array}{cc}A=\left[\begin{array}{ccc}{a}_{11}& \dots & a_{1K}\\ ⋮& \ddots & ⋮\\ a_{M1}& \dots & a_{\mathrm{MK}}\\ -1& \dots & 0\\ ⋮& \ddots & ⋮\\ 0& \dots & -1\end{array}\right],& \left(15\right)\end{array}$ where the elements a_mk, m=1, . . . , M, and k=1, . . . , K, satisfy

$\begin{array}{cc}{a}_{\mathrm{mk}}=\{\begin{array}{c}\frac{{h_{\mathrm{mk}}^{\mathrm{inv}}}^2}{\sum _{m=1}^M{h_{\mathrm{mk}}^{\mathrm{inv}}}^2},m=1,\dots ,M,k=1,\dots ,K,\\ 1,m=M+k,\\ 0,\mathrm{otherwise}\end{array}.& \left(16\right)\end{array}$

The BS constructs an (M+K)×1 vector b asb=[P₁^ant. . . P_M^ant0 . . . 0]^T.  (17)

Let the power allocation vector p^s be p^s=[p₁^s . . . P_K^s]^T, where P_k^s, k=1, . . . , K, are the power allocated to the kth data stream. Then, p^s can be obtained by solving the following optimization problemmin−Σ_k=1^K log(1+g_kP_k^s),s·t·Ap^s−b≤0,  (18)where

$g_k=\frac{{h_k}^2}{{\sigma }_{\mathrm{NI}}^2}$ denotes the signal-to-Interference-plus-Noise Ratio (SINR) of the kth stream with σ_N1² being the power of the noise and interference.

Problem (18) can be solved by iterative searching in the constraint domain. One possible solution is to firstly transform (18) into an equivalent problemmin G_t(p^s)=min[−tΣ_k=1^K log(1+g_kP_k^s)−Σ_k=1^K+M log f_i(p^s)],  (19)where f_i(p^s)=a_i^Tp^s−b_i and a_i^T, is the ith column vector of A. Then, the iterative searching process in FIG. 2 is used to find the optimal p^s. Specifically, after the searching process starts 12, the input parameters are initialized as t>0, ε_o>0, ε_i>0, p^s∈{Ap^s−b≤0, p_i≥0}, μ>0, and γ>0, where ε_o and ε_i are endurable errors for the outer and inner searching cycles respectively, while μ and γ are two adjusting parameters 13. Then, the outer searching cycle runs while

$-\frac{1}{t}\sum _{i=1}^{M+K}\log f_i\left(p^s\right)>{\varepsilon }_o$ 14, and the inner searching cycle runs while λ(p^s)/2>ε_i 15. Inside the inner cycle, the first step is to calculate the decrement Δp^s=−∇²G_t(p^s)⁻¹∇G_t(p^s) 16. After that, the Newton decrement is calculated as λ=∇G_t (p^s)^T∇²G_t(p^s)⁻¹∇G_t(p^s) 17. Then, the variable p^s is updated as p^s←p^s+γΔp^s 18. After the inner cycle ends 19, the parameter t is updated as t←μt 20. After the outer cycle ends 21, the whole process ends 22.

With the solution of problem (19), the power allocation matrix can be computed as

$\begin{array}{cc}{p}_{\mathrm{mk}}=\frac{p_k^s{h_{\mathrm{mk}}^{\mathrm{inv}}}^2}{\sum _{m=1}^M{h_{\mathrm{mk}}^{\mathrm{inv}}}^2},m=1,\dots ,M,k=1,\dots ,K.& \left(20\right)\end{array}$

With p_mk in (20), the elements w_mk of the precoding matrix W can be computed by w_mk=√{square root over (p_mk)}e^jθmk, m=1, . . . , M, and k=1, . . . , K.

If CB is employed by the BS, it firstly computes the phases of the elements of precoding matrix W asØ_mk=−θ_mk,m=1, . . . ,M,k=1, . . . ,K,  (21)then it computes the elements of the precoding matrix W as

$\begin{array}{cc}{w}_{\mathrm{mk}}=\sqrt{\frac{P_k}{M}}{e}^{j{\varnothing }_{\mathrm{mk}}},m=1,\dots ,M,k=1,\dots ,K.& \left(22\right)\end{array}$

The process of power allocation and precoding matrix computation is illustrated in FIG. 3. Specifically, after the process starts 23, the BS determines the precoding method first 24. Then, the BS computes the phase of each element of the precoding matrix θ_mk 25. Next, the BS computes the power of each data stream on each antenna p_mk 26. After that, the BS computes each elements of the precoding matrix as w_mk=p_mke^jØmk 27, before the process ends 28.

With the precoding matrix W, the signals belonging to these K data streams are precoded, further processed, and sent by the M antennas.

Although the foregoing descriptions of the preferred embodiments of the present inventions have shown, described, or illustrated the fundamental novel features or principles of the inventions, it is understood that various omissions, substitutions, and changes in the form of the detail of the methods, elements or apparatuses as illustrated, as well as the uses thereof, may be made by those skilled in the art without departing from the spirit of the present inventions. Hence, the scope of the present inventions should not be limited to the foregoing descriptions. Rather, the principles of the inventions may be applied to a wide range of methods, systems, and apparatuses, to achieve the advantages described herein and to achieve other advantages or to satisfy other objectives as well.

Claims

1. A method for power allocation in a Multi-User Multiple Input Multiple Output (MU-MIMO) Base Station (BS) with an array of M antennas comprising computing a temporary precoding matrix based on a channel matrix to generate K multiplexed data streams in a downlink in a Resource Element (RE) where M>K;computing a vector with elements that are ratios of a magnitude of elements of the temporary precoding matrix;using the vector to calculate power allocation to each data stream on each of said antennas to maximize power utilization in the MU-MIMO BS;generating a precoding matrix using the temporary precoding matrix and the calculated power allocation to generate K multiplexed data streams in the downlink in the RE, wherein the RE is a subcarrier or a resource block;allocating the calculated power to each data stream on each of said antennas using the precoding matrix; andtransmitting the generated K multiplexed data streams via the M antennas.

2. The method in claim 1, wherein calculating power allocation to maximize power utilization is by orthogonal projection.

3. The method in claim 1, wherein calculating power allocation to maximize power utilization is by iterative searching.

4. The method in claim 1, wherein calculating power allocation to maximize the power utilization comprising constructing special matrices and vectors using numbers and values related to the antennas, an allocated power to each data stream, and a channel quality, and calculating the power allocated to each data stream on each antenna using the constructed special matrices and vectors, and using results from said calculating to obtain a final beamforming matrix to maximize the power utilization.

5. The method in claim 4, further comprising calculating the power allocated to each data stream on each antenna by orthogonal projection using the constructed special matrices and vectors.

6. The method in claim 4, further comprising calculating the power allocated to each data stream on each antenna by iterative searching using the constructed special matrices and vectors.