CN105071849A - Method for realizing multi-stream beam forming in TD-LTE-Advanced - Google Patents

Abstract

The present invention discloses a method for realizing multi-stream beam forming in TD-LTE-Advanced. The method comprises acquiring a 4*8 channel matrix A1 from upstream sounding reference signals of a TD-LTE-Advanced base station, decomposing the channel matrix A1 to obtain two 4*4 sub-channel matrixes, carrying out Household transformation on each sub-channel matrix to generate an upper Hessenberg matrix J(1), carrying out Givens rotation on the upper Hessenberg matrixes J(1) to transform the upper Hessenberg matrixes J(1) into diagonal matrixes, the above operation is repeated for at least five times, a matrix obtained by previous calculation serves as a Hessenberg matrix used in the next calculation, and an obtained right multiplication Household transformation matrix and an obtained right multiplication Givens matrix are multiplied to obtain a 4*4 matrix V. A maximum ratio transmit algorithm is utilized to carry out weighting treatment on each column of the generated matrix V, thus to generate an ultimate beam forming weight vector. With adoption of the method, defects of a conventional Eigen-based beamforming (EBB) algorithm can be overcome, the multi-stream beam forming weight vector can be accurately calculated, and an error rate is effectively lowered.

Description

A kind of method realizing multi-flow beam forming in TD-LTE-Advanced

Technical field

The invention belongs to mobile communication technology field, more specifically, relate to a kind of method realizing multi-flow beam forming in TD-LTE-Advanced.

Background technology

Global communication cause is in recent years fast-developing, and the demand of radio communication is increasing, and radio communication cause obtains flourish.But along with people's continuing to increase wireless communication needs, the contradiction between huge communication requirement amount and very limited frequency spectrum resource is more and more outstanding.How to utilize limited frequency spectrum resource efficiently, and the important topic that large-scale raising power system capacity has become radio communication circle urgently to be resolved hurrily under the prerequisite of ensuring the quality of products, the wave beam forming in smart antenna has become the important directions addressed this problem.

Wave beam forming is applied to closely spaced aerial array transmission technology, the strong correlation in space and the principle of interference of ripple is utilized to produce the antenna pattern of directivity, make the main lobe of antenna pattern point to arrival bearing user adaptively, thus improve signal to noise ratio, power system capacity and coverage.

The algorithm of traditional compute beam forming weight vector is the algorithm of the wave beam forming (Eigen-basedBeamforming) that feature based decomposes, although this algorithm realizes simple, but the deficiency of depositing when calculating multi-flow beam forming weight vector both ways: () is very high to the selective dependency of initial vector, if good primary iteration vector can not be selected, may cause being difficult to convergence, wave beam forming weight vector can not be tried to achieve; (2) when channel matrix exists same characteristic features value, the wave beam forming weight vector of single current can only be obtained, cause the higher error rate the most at last.

Summary of the invention

For above defect or the Improvement requirement of prior art, the invention provides a kind of method realizing multi-flow beam forming weight vector in TD-LTE-Advanced, its object is to, overcome the deficiency of existing EBB algorithm, calculate multi-flow beam forming weight vector exactly, and effectively reduce the error rate.

For achieving the above object, according to one aspect of the present invention, provide a kind of method realizing multi-flow beam forming in TD-LTE-Advanced, comprise the following steps:

(1) the channel matrix A of 4 × 8 is obtained from the uplink detection reference signal of TD-LTE-Advanced base station ¹;

(2) by channel matrix A ¹decompose, to obtain two 4 × 4 sub-channel matrix A ^{(p, 1)}, it corresponds respectively to two sub-aerial arrays, and wherein p represents the sequence number of sub antenna array;

(3) to each sub-channel matrix A ^{(p, 1)}carry out Household conversion, to generate upper Hessenberg matrix J ⁽¹⁾

(4) to upper Hessenberg matrix J ⁽¹⁾carry out Givens rotation, with by this matrix J ⁽¹⁾be transformed into diagonal matrix;

(5) computing repeating above-mentioned steps (4) reaches at least 5 times, and wherein the last matrix calculated again can as the Hessenberg matrix used in calculating next time;

(6) taken advantage of on the right side obtained in step (3) right side obtained in Household transformation matrix and step (4) to take advantage of Givens matrix to fold to take advantage of, to obtain the matrix V of a 4*4, namely each row of this matrix V are designated as v ₁, v ₂, v ₃, v ₄.

(7) utilize high specific transmission algorithm to each row v of the matrix V that step (6) generates ₁, v ₂, v ₃, v ₄be weighted process, to generate final wave beam forming weight vector v' ₁, v' ₂, v' ₃, v' ₄:

v′ ₁＝p ₁v ₁

v′ ₂＝p ₂v ₂

v′ ₃＝p ₃v ₃

v′ ₄＝p ₄v ₄

Wherein p is the transmitting power factor.

Preferably, step (3) comprises following sub-step:

(3-1) counter k=1 is set;

(3-2) premultiplication Household matrix L is constructed ^(k), use this matrix to A ^{(p, k)}carry out premultiplication, to obtain matrix A ^{(p, k+1/2)}=L ^(k)a ^{(p, k)};

(3-3) judged whether that k+1=4 sets up, if it is matrix A ^{(p, 1)}become upper Hessenberg matrix J ⁽¹⁾, process terminates, otherwise enters step (3-4);

(3-4) construct the right side and take advantage of Household matrix R ^(k), use the right side to take advantage of Household matrix R ^(k)to A ^{(p, k+1/2)}carry out the right side to take advantage of, to obtain matrix A ^{(p, k+1)}=A ^{(p, k+1/2)}r ^(k);

(3-5) k=k+1 is set, and returns step (3-2);

Preferably, premultiplication Household matrix is constructed in the following ways:

(3-2-1) first zoom factor K is calculated _x, order represent and be transformed matrix A ^{(p, k)}in i-th row jth row element:

$K_x=1+\sqrt{\frac{\underset{i=k}{\overset{N}{\mathrm{\Σ}}}{a}_{i,k}^{\left(k\right)}{a}_{i,k}^{\left(k\right)*}}{a_{k,k}^{\left(k\right)}{a}_{k,k}^{\left(k\right)*}}}$

(3-2-2) column vector { xtemp} is constructed ^k:

(3-2-3) by { xtemp} ^kbe converted into unit vector { x} ^(k):

${\left\{x\right\}}^{\left(k\right)}=\frac{1}{\sqrt{\underset{i=k}{\overset{N}{\mathrm{\Σ}}}{\mathrm{xtemp}}_i^{\left(k\right)}{\mathrm{xtemp}}_i^{\left(k\right)*}}}{\left\{\mathrm{xtemp}\right\}}^k$

(3-2-4) according to { x} ^(k)structure premultiplication Household matrix L ^(k):

L ^(k)=I-2{x} ^(k){ x} ^{(k) *}, wherein I is the unit matrix of 4*4.

Preferably, construct the right side and take advantage of Household matrix in the following ways:

(3-4-1) first zoom factor K is calculated _y, order represent and be transformed matrix A ^{(p, k+1/2)}in row k i-th arrange element

$K_y=1+\sqrt{\frac{\underset{i=k+1}{\overset{N}{\mathrm{\Σ}}}{a}_{k,i}^{(k+1/2)}{a}_{k,i}^{(k+1/2)*}}{a_{k,k+1}^{(k+1/2)}{a}_{k,k+1}^{(k+1/2)*}}}$

(3-4-2) column vector { ytemp} is constructed ^k:

(3-4-3) by { ytemp} ^kturn to unit vector { y} ^(k):

${\left\{y\right\}}^{\left(k\right)}=\frac{1}{\sqrt{\underset{i=k+1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{ytemp}}_i^{\left(k\right)}{\mathrm{ytemp}}_i^{\left(k\right)*}}}{\left\{\mathrm{ytemp}\right\}}^k$

(3-4-4) according to { y} ^(k)structure R ^(k)

R ^(k)=I-2{y} ^(k){ y} ^{(k) *}, wherein I is the unit matrix of 4*4.

Preferably, step (4) comprises following sub-step:

(4-1) counter m=1 is set;

(4-2) dextrorotation Givens matrix Q is constructed ^(m), to J ^(m)carry out the right side to take advantage of, to obtain J ^(m+1/2).That is: J ^(m+1/2)=J ^(m)q ^(m);

(4-3) left-handed Givens matrix P is constructed ^(m), to J ^(m+1/2)carry out premultiplication, obtain J ^(m+1).That is: J ^(m+1)=P ^(m)j ^(m+1/2);

(4-4) counter m=m+1 is set;

(4-5) judge whether m=4, if process terminates, otherwise return step (4-2).

Preferably, construct the right side and take advantage of Givens matrix in the following ways:

(4-2-1) initial vector is set wherein for J ^(m)in the i-th row jth column element;

(4-2-2) by [f ^(m), g ^(m)] calculate Givens dextrorotation matrix Q ^(m), for Q ^(m)element on middle diagonal and l ≠ m and l ≠ m+1;

(4-2-3) Q is made ^(m)middle matrix-block

$\left[\begin{array}{cc}{q}_{m,m}^{\left(m\right)}& q_{m,m+1}^{\left(m\right)}\\ q_{m+1,m}^{\left(m\right)}& q_{m+1,m+1}^{\left(m\right)}\end{array}\right]=\left[\begin{array}{cc}{f}^{\left(m\right)*}/r& -g^{\left(m\right)}{f}^{\left(m\right)*}/(r.R_f)\\ g^{\left(m\right)*}/r& R_f/r\end{array}\right],$

Wherein r is the mould of initial vector, R _ffor f ^(m)mould;

(4-2-4) for all the other elements, 0 is got.

Preferably, premultiplication Givens matrix is constructed in the following ways:

(4-3-1) initial vector is set

${[f^{(m+1/2)},g^{(m+1/2)}]}^T={[j_{m,m}^{(m+1/2)},j_{m+1,m}^{(m+1/2)}]}^T;$ Wherein for J ^(m+1/2)in i-th row jth row element;

(4-3-2) by [f ^(m+1/2), g ^(m+1/2)] ^tcalculate the left-handed matrix P of Givens ^(m), for the element on diagonal and l ≠ m and l ≠ m+1;

(4-3-3) order matrix block

$\left[\begin{array}{cc}{p}_{m,m}^{\left(m\right)}& p_{m,m+1}^{\left(m\right)}\\ p_{m+1,m}^{\left(m\right)}& p_{m+1,m+1}^{\left(m\right)}\end{array}\right]=\left[\begin{array}{cc}{f}^{(m+1/2)*}/r& g^{(m+1/2)*}/r\\ -g^{(m+1/2)}{f}^{(m+1/2)*}/(r\·R_f^{\left(m\right)})& R_f^{(m+1/2)}/r\end{array}\right],$ Wherein r is the mould of initial vector, R _ffor f ^(m+1/2)mould;

(4-3-4) for all the other elements, 0 is got.

Preferably, weighted factor is tried to achieve by following formula:

$\begin{array}{c}{p}_1=\frac{{\mathrm{\σ}}_1^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\\ p_2=\frac{{\mathrm{\σ}}_2^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\\ p_3=\frac{{\mathrm{\σ}}_3^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\\ p_4=\frac{{\mathrm{\σ}}_4^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\end{array}$

Wherein σ ₁, σ ₂, σ ₃, σ ₄be respectively sub-channel matrix A ^{(p, 1)}singular value.

According to another aspect of the present invention, provide a kind of multi-flow beam forming in TD-LTE-Advanced of realizing, comprising:

First module, for obtaining the channel matrix A of 4 × 8 from the uplink detection reference signal of TD-LTE-Advanced base station ¹;

Second module, for by channel matrix A ¹decompose, to obtain two 4 × 4 sub-channel matrix A ^{(p, 1)}, it corresponds respectively to two sub-aerial arrays, and wherein p represents the sequence number of sub antenna array;

3rd module, for each sub-channel matrix A ^{(p, 1)}carry out Household conversion, to generate upper Hessenberg matrix J ⁽¹⁾

Four module, for upper Hessenberg matrix J ⁽¹⁾carry out Givens rotation, with by this matrix J ⁽¹⁾be transformed into diagonal matrix;

5th module, reaches at least 5 times for the computing repeating above-mentioned four module, and wherein the last matrix calculated again can as the Hessenberg matrix used in calculating next time;

6th module, the right side that the right side for the 3rd module being obtained takes advantage of Household transformation matrix and four module to obtain is taken advantage of Givens matrix to fold and is taken advantage of, to obtain the matrix V of a 4*4, namely each row of this matrix V are designated as v ₁, v ₂, v ₃, v ₄.

7th module, for utilizing high specific transmission algorithm to each row v of the matrix V of the 6th CMOS macro cell ₁, v ₂, v ₃, v ₄be weighted process, to generate final wave beam forming weight vector v' ₁, v' ₂, v' ₃, v' ₄:

v′ ₁＝p ₁v ₁

v′ ₂＝p ₂v ₂

v′ ₃＝p ₃v ₃

v′ ₄＝p ₄v ₄

Wherein p is the transmitting power factor.

In general, the above technical scheme conceived by the present invention compared with prior art, can obtain following beneficial effect:

1, the present invention can overcome the very high problem of the selective dependency to initial vector that exists in existing method: because the present invention does not adopt initial vector, therefore the present invention does not rely on the selection of initial vector;

2, the singular value decomposition of channel matrix is achieved owing to present invention employs Householder transformation and Givens rotation.Therefore, the present invention, when channel matrix exists same characteristic features value, also can obtain the wave beam forming weight vector of multithread, this reduces the error rate.

Accompanying drawing explanation

Fig. 1 is the aerial array schematic diagram that the present invention uses.

Fig. 2 is the flow chart that the present invention realizes the method for multi-flow beam forming in TD-LTE-Advanced.

Embodiment

In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with drawings and Examples, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.In addition, if below in described each execution mode of the present invention involved technical characteristic do not form conflict each other and just can mutually combine.

Basic ideas of the present invention are, by carrying out Householder transformation to channel matrix, channel matrix is converted into upper Hessenberg matrix, afterwards Givens rotation is carried out to upper Hessenberg matrix, be translated into diagonal matrix, so just achieve the singular value decomposition of channel matrix, each weight vector arranged as wave beam forming obtaining the right singular value matrix of channel matrix is used.

The present invention adopts 4 × 4 cross polarised antenna arrays in Fig. 1, and this aerial array is made up of two subarrays, and 1,2,3,4 composition subarrays 1, adopt-45 ° of polarization, 5,6, and 7,8 composition subarrays 2, adopt+45 ° of polarization.

This aerial array pattern supports 2.6GHz carrier frequency, and base station side of the present invention configures eight antennas, and user side configures four antennas.

As shown in Figure 2, the present invention realizes the method for multi-flow beam forming in TD-LTE-Advanced, comprises the following steps:

(1) the channel matrix A of 4 × 8 is obtained from the uplink detection reference signal (soundingreferencesignal is called for short SRS) of TD-LTE-Advanced base station ¹;

(2) by channel matrix A ¹decompose, to obtain two 4 × 4 sub-channel matrix A ^{(p, 1)}, it corresponds respectively to two sub-aerial arrays, and wherein p represents the sequence number of sub antenna array;

(3) to each sub-channel matrix A ^{(p, 1)}carry out Household conversion, to generate upper Hessenberg matrix J ⁽¹⁾, this step comprises following sub-step:

(3-1) counter k=1 is set;

(3-2) premultiplication Household matrix L is constructed ^(k), use this matrix to A ^{(p, k)}carry out premultiplication, to obtain matrix A ^{(p, k+1/2)}, that is: A ^{(p, k+1/2)}=L ^(k)a ^{(p, k)}; Wherein construct premultiplication Household matrix in the following ways:

(3-2-1) first zoom factor K is calculated _x, order represent and be transformed matrix A ^{(p, k)}in i-th row jth row element:

$K_x=1+\sqrt{\frac{\underset{i=k}{\overset{N}{\mathrm{\Σ}}}{a}_{i,k}^{\left(k\right)}{a}_{i,k}^{\left(k\right)*}}{a_{k,k}^{\left(k\right)}{a}_{k,k}^{\left(k\right)*}}}$

(3-2-2) column vector { xtemp} is constructed ^k:

(3-2-3) by { xtemp} ^kbe converted into unit vector { x} ^(k):

${\left\{x\right\}}^{\left(k\right)}=\frac{1}{\sqrt{\underset{i=k}{\overset{N}{\mathrm{\Σ}}}{\mathrm{xtemp}}_i^{\left(k\right)}{\mathrm{xtemp}}_i^{\left(k\right)*}}}{\left\{\mathrm{xtemp}\right\}}^k$

(3-2-4) according to { x} ^(k)structure premultiplication Household matrix L ^(k):

L ^(k)=I-2{x} ^(k){ x} ^{(k) *}; Wherein I is the unit matrix of 4*4;

(3-3) judged whether that k+1=4 sets up, if it is matrix A ^{(p, 1)}become upper Hessenberg matrix J ⁽¹⁾, process terminates, otherwise enters step (3-4);

(3-4) construct the right side and take advantage of Household matrix R ^(k), use the right side to take advantage of Household matrix R ^(k)to A ^{(p, k+1/2)}carry out the right side to take advantage of, to obtain matrix A ^{(p, k+1)}, that is: A ^{(p, k+1)}=A ^{(p, k+1/2)}r ^(k); Wherein construct the right side and take advantage of Household matrix in the following ways:

(3-4-1) first zoom factor K is calculated _y, order represent and be transformed matrix A ^{(p, k+1/2)}in row k i-th arrange element

$K_y=1+\sqrt{\frac{\underset{i=k+1}{\overset{N}{\mathrm{\Σ}}}{a}_{k,i}^{(k+1/2)}{a}_{k,i}^{(k+1/2)*}}{a_{k,k+1}^{(k+1/2)}{a}_{k,k+1}^{(k+1/2)*}}}$

(3-4-2) column vector { ytemp} is constructed ^k:

(3-4-3) by { ytemp} ^kturn to unit vector { y} ^(k):

${\left\{y\right\}}^{\left(k\right)}=\frac{1}{\sqrt{\underset{i=k+1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{ytemp}}_i^{\left(k\right)}{\mathrm{ytemp}}_i^{\left(k\right)*}}}{\left\{\mathrm{ytemp}\right\}}^k$

(3-4-4) according to { y} ^(k)structure R ^(k), I is the unit matrix of 4*4:

R ^(k)＝I-2{y} ^(k){y} ^(k)*

(3-5) k=k+1 is set, and returns step (3-2);

(4) to upper Hessenberg matrix J ⁽¹⁾carry out Robin Givens (Givens) to rotate, with by this matrix J ⁽¹⁾be transformed into diagonal matrix, this step comprises following sub-step:

(4-1) counter m=1 is set;

(4-2) dextrorotation Givens matrix Q is constructed ^(m), to J ^(m)carry out the right side to take advantage of, to obtain J ^(m+1/2).That is: J ^(m+1/2)=J ^(m)q ^(m); Wherein construct premultiplication Givens matrix in the following ways:

(4-2-1) initial vector is set wherein for J ^(m)in the i-th row jth column element;

(4-2-2) by [f ^(m), g ^(m)] calculate Givens dextrorotation matrix Q ^(m), for Q ^(m)element on middle diagonal and l ≠ m and l ≠ m+1

(4-2-3) Q is made ^(m)middle matrix-block

$\left[\begin{array}{cc}{q}_{m,m}^{\left(m\right)}& q_{m,m+1}^{\left(m\right)}\\ q_{m+1,m}^{\left(m\right)}& q_{m+1,m+1}^{\left(m\right)}\end{array}\right]=\left[\begin{array}{cc}{f}^{\left(m\right)*}/r& -g^{\left(m\right)}{f}^{\left(m\right)*}/(r.R_f)\\ g^{\left(m\right)*}/r& R_f/r\end{array}\right],$

Wherein r is the mould of initial vector, R _ffor f ^(m)mould.

(4-2-4) for all the other elements, 0 is got.

(4-3) left-handed Givens matrix P is constructed ^(m), to J ^(m+1/2)carry out premultiplication, obtain J ^(m+1).That is: J ^(m+1)=P ^(m)j ^(m+1/2); Wherein construct premultiplication Givens matrix in the following ways:

(4-3-1) initial vector is set

${[f^{(m+1/2)},g^{(m+1/2)}]}^T={[j_{m,m}^{(m+1/2)},j_{m+1,m}^{(m+1/2)}]}^T;$ for J ^(m+1/2)in i-th row jth row element.

(4-3-2) by [f ^(m+1/2), g ^(m+1/2)] ^tcalculate the left-handed matrix P of Givens ^(m), for the element on diagonal and l ≠ m and l ≠ m+1

(4-3-3) order matrix block

$\left[\begin{array}{cc}{p}_{m,m}^{\left(m\right)}& p_{m,m+1}^{\left(m\right)}\\ p_{m+1,m}^{\left(m\right)}& p_{m+1,m+1}^{\left(m\right)}\end{array}\right]=\left[\begin{array}{cc}{f}^{(m+1/2)*}/r& g^{(m+1/2)*}/r\\ -g^{(m+1/2)}{f}^{(m+1/2)*}/(r\·R_f^{\left(m\right)})& R_f^{(m+1/2)}/r\end{array}\right],$ Wherein r is the mould of initial vector, R _ffor f ^(m+1/2)mould.

(4-3-4) for all the other elements, 0 is got.

(4-4) counter m=m+1 is set;

(4-5) judge whether m=4, if process terminates, otherwise return step (4-2).

(5) computing repeating above-mentioned steps (4) reaches at least 5 times, and wherein the last matrix calculated again can as the Hessenberg matrix used in calculating next time;

(6) taken advantage of on the right side obtained in step (3) right side obtained in Household transformation matrix and step (4) to take advantage of Givens matrix to fold to take advantage of, to obtain the matrix V of a 4*4, namely each row of this matrix V are designated as v ₁, v ₂, v ₃, v ₄.

(7) utilize high specific transmission algorithm (MaximumRatioTransmission is called for short MRT) to each row v of the matrix V that step (6) generates ₁, v ₂, v ₃, v ₄be weighted process, to generate final wave beam forming weight vector v' ₁, v' ₂, v' ₃, v' ₄:

v′ ₁＝p ₁v ₁

v′ ₂＝p ₂v ₂

v′ ₃＝p ₃v ₃

v′ ₄＝p ₄v ₄

Wherein p is the transmitting power factor, and it is tried to achieve by following formula:

$\begin{array}{c}{p}_1=\frac{{\mathrm{\σ}}_1^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\\ p_2=\frac{{\mathrm{\σ}}_2^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\\ p_3=\frac{{\mathrm{\σ}}_3^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\\ p_4=\frac{{\mathrm{\σ}}_4^2}{{\mathrm{\σ}}_1^2+{\mathrm{\σ}}_2^2+{\mathrm{\σ}}_3^2+{\mathrm{\σ}}_4^2}\end{array}$

Wherein σ ₁, σ ₂, σ ₃, σ ₄be respectively sub-channel matrix A ^{(p, 1)}singular value.

Generally speaking, the present invention has following beneficial effect:

1, the present invention can overcome the very high problem of the selective dependency to initial vector that exists in existing method: because the present invention does not adopt initial vector, therefore the present invention does not rely on the selection of initial vector;

2, the singular value decomposition of channel matrix is achieved owing to present invention employs Householder transformation and Givens rotation, therefore, the present invention, when channel matrix exists same characteristic features value, also can obtain the wave beam forming weight vector of multithread, this reduces the error rate.

Those skilled in the art will readily understand; the foregoing is only preferred embodiment of the present invention; not in order to limit the present invention, all any amendments done within the spirit and principles in the present invention, equivalent replacement and improvement etc., all should be included within protection scope of the present invention.

Claims (9)

1. realize a method for multi-flow beam forming in TD-LTE-Advanced, it is characterized in that, comprise the following steps:

(1) the channel matrix A of 4 × 8 is obtained from the uplink detection reference signal of TD-LTE-Advanced base station ¹;

(2) by channel matrix A ¹decompose, to obtain two 4 × 4 sub-channel matrix A ^{(p, 1)}, it corresponds respectively to two sub-aerial arrays, and wherein p represents the sequence number of sub antenna array;

(3) to each sub-channel matrix A ^{(p, 1)}carry out Household conversion, to generate upper Hessenberg matrix J ⁽¹⁾

(4) to upper Hessenberg matrix J ⁽¹⁾carry out Givens rotation, with by this matrix J ⁽¹⁾be transformed into diagonal matrix;

(5) computing repeating above-mentioned steps (4) reaches at least 5 times, and wherein the last matrix calculated again can as the Hessenberg matrix used in calculating next time;

(6) taken advantage of on the right side obtained in step (3) right side obtained in Household transformation matrix and step (4) to take advantage of Givens matrix to fold to take advantage of, to obtain the matrix V of a 4*4, namely each row of this matrix V are designated as v ₁, v ₂, v ₃, v ₄.

(7) utilize high specific transmission algorithm to each row v of the matrix V that step (6) generates ₁, v ₂, v ₃, v ₄be weighted process, to generate final wave beam forming weight vector v ' ₁, v ' ₂, v ' ₃, v ' ₄:

v′ ₁＝p ₁v ₁

v′ ₂＝p ₂v ₂

v′ ₃＝p ₃v ₃

v′ ₄＝p ₄v ₄

Wherein p is the transmitting power factor.

2. method according to claim 1, is characterized in that, step (3) comprises following sub-step:

(3-1) counter k=1 is set;

(3-2) premultiplication Household matrix L is constructed ^(k), use this matrix to A ^{(p, k)}carry out premultiplication, to obtain matrix A ^{(p, k+1/2)}=L ^(k)a ^{(p, k)};

(3-3) judged whether that k+1=4 sets up, if it is matrix A ^{(p, 1)}become upper Hessenberg matrix J ⁽¹⁾, process terminates, otherwise enters step (3-4);

(3-4) construct the right side and take advantage of Household matrix R ^(k), use the right side to take advantage of Household matrix R ^(k)to A ^{(p, k+1/2)}carry out the right side to take advantage of, to obtain matrix A ^{(p, k+1)}=A ^{(p, k+1/2)}r ^(k);

(3-5) k=k+1 is set, and returns step (3-2).

3. method according to claim 2, is characterized in that, constructs premultiplication Household matrix in the following ways:

(3-2-1) first zoom factor K is calculated _x, order represent and be transformed matrix A ^{(p, k)}in i-th row jth row element:

(3-2-2) column vector { xtemp} is constructed ^k:

(3-2-3) by { xtemp} ^kbe converted into unit vector { x} ^(k):

(3-2-4) according to { x} ^(k)structure premultiplication Household matrix L ^(k):

L ^(k)=I-2{x} ^(k){ x} ^{(k) *}, wherein I is the unit matrix of 4*4.

4. method according to claim 3, is characterized in that, Household matrix is taken advantage of in the following ways in the structure right side:

(3-4-1) first zoom factor K is calculated _y, order represent and be transformed matrix A ^{(p, k+1/2)}in row k i-th arrange element

(3-4-2) column vector { ytemp} is constructed ^k:

(3-4-3) by { ytemp} ^kturn to unit vector { y} ^(k):

(3-4-4) according to { y} ^(k)structure R ^(k)

R ^(k)=I-2{y} ^(k){ y} ^{(k) *}, wherein I is the unit matrix of 4*4.

5. method according to claim 4, is characterized in that, step (4) comprises following sub-step:

(4-1) counter m=1 is set;

(4-2) dextrorotation Givens matrix Q is constructed ^(m), to J ^(m)carry out the right side to take advantage of, to obtain J ^(m+1/2).That is: J ^(m+1/2)=J ^(m)q ^(m);

(4-3) left-handed Givens matrix P is constructed ^(m), to J ^(m+1/2)carry out premultiplication, obtain J ^(m+1).That is: J ^(m+1)=P ^(m)j ^(m+1/2);

(4-4) counter m=m+1 is set;

(4-5) judge whether m=4, if process terminates, otherwise return step (4-2).

6. method according to claim 5, is characterized in that, Givens matrix is taken advantage of in the following ways in the structure right side:

(4-2-1) initial vector is set wherein for J ^(m)in the i-th row jth column element;

(4-2-2) by [f ^(m), g ^(m)] calculate Givens dextrorotation matrix Q ^(m), for Q ^(m)element on middle diagonal and l ≠ m and l ≠ m+1;

(4-2-3) Q is made ^(m)middle matrix-block

Wherein r is the mould of initial vector, R _ffor f ^(m)mould;

(4-2-4) for all the other elements, 0 is got.

7. method according to claim 6, is characterized in that, constructs premultiplication Givens matrix in the following ways:

(4-3-1) initial vector is set

wherein for J ^(m+1/2)in i-th row jth row element;

(4-3-2) by [f ^(m+1/2), g ^(m+1/2)] ^tcalculate the left-handed matrix P of Givens ^(m), for the element on diagonal and l ≠ m and l ≠ m+1;

(4-3-3) order matrix block

wherein r is the mould of initial vector, R _ffor f ^(m+1/2)mould;

(4-3-4) for all the other elements, 0 is got.

8. method according to claim 7, is characterized in that, weighted factor is tried to achieve by following formula:

Wherein σ ₁, σ ₂, σ ₃, σ ₄be respectively sub-channel matrix A ^{(p, 1)}singular value.

9. one kind realizes multi-flow beam forming in TD-LTE-Advanced, it is characterized in that, comprising:

First module, for obtaining the channel matrix A of 4 × 8 from the uplink detection reference signal of TD-LTE-Advanced base station ¹;

Second module, for by channel matrix A ¹decompose, to obtain two 4 × 4 sub-channel matrix A ^{(p, 1)}, it corresponds respectively to two sub-aerial arrays, and wherein p represents the sequence number of sub antenna array;

3rd module, for each sub-channel matrix A ^{(p, 1)}carry out Household conversion, to generate upper Hessenberg matrix J ⁽¹⁾

Four module, for upper Hessenberg matrix J ⁽¹⁾carry out Givens rotation, with by this matrix J ⁽¹⁾be transformed into diagonal matrix;

5th module, reaches at least 5 times for the computing repeating above-mentioned four module, and wherein the last matrix calculated again can as the Hessenberg matrix used in calculating next time;

6th module, the right side that the right side for the 3rd module being obtained takes advantage of Household transformation matrix and four module to obtain is taken advantage of Givens matrix to fold and is taken advantage of, to obtain the matrix V of a 4*4, namely each row of this matrix V are designated as v ₁, v ₂, v ₃, v ₄.

7th module, for utilizing high specific transmission algorithm to each row v of the matrix V of the 6th CMOS macro cell ₁, v ₂, v ₃, v ₄be weighted process, to generate final wave beam forming weight vector v ' ₁, v ' ₂, v ' ₃, v ' ₄:

v′ ₁＝p ₁v ₁

v′ ₂＝p ₂v ₂

v′ ₃＝p ₃v ₃

v′ ₄＝p ₄v ₄

Wherein p is the transmitting power factor.