Embodiment
The present invention carries out the QR resolution process by the embodiment of the invention to channel matrix, and according to the value of the elements in a main diagonal that obtains the upper triangular matrix R that obtains after QR decomposes, determines the order of this channel matrix; Thereby it is large to have solved the method operand that has the order of determining channel matrix now, and the low problem of reliability of the order of determining.
Channel matrix can obtain by channel situation being carried out channel estimating in the embodiment of the invention, with LTE (Long Term Evolution, Long Term Evolution) system is example, UE (User Equipment, subscriber equipment) can carry out channel estimating by cell common pilots (CRS) or channel information measurement pilot tone (CSI RS), in certain physical resource scope, obtain channel matrix.
Below in conjunction with Figure of description the embodiment of the invention is described in further detail.
As shown in Figure 1, the embodiment of the invention is determined the method for the order of channel matrix, may further comprise the steps:
Step 101, channel matrix is carried out the QR resolution process, obtain the value of the elements in a main diagonal of the upper triangular matrix R that obtains after QR decomposes;
Wherein, channel matrix H is N
r* N
tDimension matrix, wherein N
rThe number of expression reception antenna, N
tThe number of expression transmitting antenna carries out obtaining unitary matrice Q and upper triangular matrix R after QR decomposes to channel matrix H;
Step 102, according to the value of the elements in a main diagonal in the R matrix, determine the order of channel matrix.
In specific implementation process, determine in the step 102 that the order of channel matrix H comprises following dual mode:
Mode one: according to the value of the elements in a main diagonal in the R matrix and the threshold value of setting, determine the order of channel matrix; Wherein threshold value is empirical value, can set as required;
In specific implementation process, determine that the order of channel matrix comprises:
Determine that R matrix the elements in a main diagonal intermediate value is greater than the number of the elements in a main diagonal of the threshold value of setting;
With the number of the elements in a main diagonal the determined order as channel matrix.
Suppose that threshold value is r_th, after channel matrix H was decomposed through QR, main diagonal element was r in the R matrix that obtains
I, i, wherein i ∈ 1 ..., N
t, then from i=1, judge one by one r
I, iValue whether greater than threshold value r_th, until i=N
tIf the result who judges is: have N main diagonal element must be worth greater than threshold value r_th, determine that then the order of channel matrix H is N;
Preferably, before determining the order of channel matrix, also comprise in the step 102: with the processing of sorting of the value of the elements in a main diagonal of R matrix.
In specific implementation process, if the value to the elements in a main diagonal is carried out descending, then from first the elements in a main diagonal, whether the value of the elements in a main diagonal after definite descending is greater than the threshold value of setting successively, after the value of determining the elements in a main diagonal that first is not more than threshold value, determine to be positioned at the value of number of all the elements in a main diagonals before this elements in a main diagonal as the order of channel matrix;
Suppose N
t=4, the elements in a main diagonal of the R matrix that obtains after the process step 101 is respectively: r
1,1, r
2,2, r
3,3And r
4,4, with r
1,1, r
2,2, r
3,3And r
4,4Value carry out descending, make r
1,1, r
2,2, r
3,3And r
4,4Value satisfy: r
1,1>r
2,2>r
3,3>r
4,4, after judging, if r
1,1>r_th, r
2,2>r_th, r
3,3<r_th then need to not judge r
4,4Value and the relation of r_th, the order of determining channel matrix H is 2;
If the value of the elements in a main diagonal is carried out ascending order arranges, then from first the elements in a main diagonal, determine that successively whether the value of the elements in a main diagonal after ascending order is arranged is greater than threshold value, after determining first value greater than the elements in a main diagonal of threshold value, the value of number of determining to comprise this elements in a main diagonal and being positioned at all the elements in a main diagonals after this elements in a main diagonal is as the order of channel matrix;
Suppose N
t=4, the elements in a main diagonal of the R matrix that obtains after the process step 101 is respectively: r
1,1, r
2,2, r
3,3And r
4,4, with r
1,1, r
2,2, r
3,3And r
4,4Value carry out ascending order and arrange, make r
1,1, r
2,2, r
3,3And r
4,4Value satisfy: r
1,1<r
2,2<r
3,3<r
4,4, after judging, if r
1,1<r_th, r
2,2>r_th then need to not judge r
3,3And r
4,4Value and the relation of r_th, the order of determining channel matrix H is 3.
Mode two, with the processing of sorting of the value of the elements in a main diagonal of R matrix, the ratio of the value of adjacent two the elements in a main diagonals and these adjacent two threshold values that the elements in a main diagonal is corresponding of setting after processing according to ordering are determined the order of channel matrix;
In specific implementation process, if the value of the elements in a main diagonal of R matrix is carried out descending, then from first the elements in a main diagonal r
1,1Beginning is determined the ratio of the value of the value of current the elements in a main diagonal and adjacent next the elements in a main diagonal, i.e. Ratio successively
i=r
I, i/ r
I+1, i+1, the ratio R atio that is determining
iGreater than these adjacent two threshold value Ratio_th that the elements in a main diagonal is corresponding
iThe time, the order of determining channel matrix H is current the elements in a main diagonal r
I, iLocation label i, two of arbitrary neighborhood the elements in a main diagonals corresponding threshold value Ratio_th all wherein
i, this threshold value is empirical value, can set as required;
Need to prove, if all Ratio
iAll greater than the threshold value Ratio_th of correspondence
i, determine that then the order of channel matrix H is N
t
Suppose N
t=4, the elements in a main diagonal of the R matrix that obtains after the process step 101 is respectively: r
1,1, r
2,2, r
3,3And r
4,4, with r
1,1, r
2,2, r
3,3And r
4,4Value carry out descending, make r
1,1, r
2,2, r
3,3And r
4,4Value satisfy: r
1,1>r
2,2>r
3,3>r
4,4, adjacent r
1,1With r
2,2Corresponding threshold value is Rati
o 1, adjacent r
2,2With r
3,3Corresponding threshold value is Ratio_th
2, adjacent r
3,3With r
4,4Corresponding threshold value is Ratio_th
3, then determine the order of channel matrix H by following process:
Make i=1, calculate Ratio
1=r
1,1/ r
2,2
If Ratio
1≤ Ratio_th
1, then calculate Ratio
2=r
2,2/ r
3,3If Ratio
1>Ratio_th
1, determine that then the order of channel matrix H is 1;
If Ratio
1≤ Ratio_th
1And Ratio
2≤ Ratio_th
2, then calculate
If Ratio
1≤ Ratio_th
1And Ratio
2>Ratio_th
2, determine that then the order of channel matrix H is 2;
If Ratio
1≤ Ratio_th
1, and Ratio
2≤ Ratio_th
2, and Ratio
3≤ Ratio_th
3, determine that then the order of channel matrix H is 4; If Ratio
1≤ Ratio_th
1, and Ratio
2≤ Ratio_th
2, and Ratio
3>Ratio_th
3, determine that then the order of channel matrix H is 3.
If being carried out ascending order, the value of the elements in a main diagonal of R matrix arranges, then from last the elements in a main diagonal
The ratio R atio of the value of the value of current the elements in a main diagonal and an adjacent upper the elements in a main diagonal is determined in beginning successively
i=r
I, i/ r
I-1, i-1, the ratio R atio that is determining
iThreshold value Ratio_th greater than correspondence
iThe time, the order of determining channel matrix is the number N of the elements in a main diagonal in the R matrix
tWith a upper the elements in a main diagonal r
I-1, i-1The difference of location label i-1, namely the order of channel matrix H is: N
t-(i-1).
In specific implementation process, in the step 101 channel matrix is carried out the QR resolution process and comprise the following steps:
For i column vector in Q matrix and the R matrix, be listed as N according to current i in the Q matrix
tThe mould value of column vector, to the processing of reordering of current i column vector in Q matrix and the R matrix, i=1 wherein ..., N
t, N
tNumber for transmitting antenna;
Be listed as N according to current i in the Q matrix that reorders after processing
tThe value of element in the capable vector of current i in the value of element and the R matrix determines that i is listed as N in the Q matrix in the column vector
tThe value of each element in the column vector, and current i is listed as N in the replacement Q matrix
tThe value of each element in the column vector, and determine in the R matrix value of each element in the capable vector of i, and replace in the R matrix value of each element in the capable vector of current i.
In specific implementation process, to the processing of reordering of current i column vector in Q matrix and the R matrix, comprise the following steps: in the step 101
Determine that current i is listed as N in the Q matrix
tThe mould value that column vector is corresponding;
The column vector that current i column vector in the Q matrix is corresponding with maximum mould value exchanges processing, and will determine in the R matrix column vector identical with column vector position corresponding to maximum mould value in the Q matrix, and current i column vector in the R matrix and the column vector of determining are exchanged processing.
Be elaborated below in conjunction with the implementation process of in 2 pairs of embodiment of the invention steps 101 of accompanying drawing channel matrix H being carried out the QR resolution process.
Step 201, Q matrix and R matrix are carried out initialization, the initial Q matrix is channel matrix H, and the Initial R matrix is N
t* N
tThe null matrix of dimension is set a 1 * N
tDimensional vector k=[1,1 ... 1];
Step 202, for i column vector in i column vector and the R matrix in the Q matrix, judge that whether i is less than N
t, i=1 wherein ..., N
t, N
tNumber for transmitting antenna;
If then execution in step 203;
If not, the elements in a main diagonal of then process ends, and output R matrix;
Step 203, according to formula 1, calculate i row the~the N in the Q matrix
tThe performance number P of column vector
j, wherein j ∈ i ..., N
t, with each performance number P
jStore among the vectorial E, namely
And make the 1st value to (i-1) individual element among the vectorial E be set to-inf (negative infinite);
P
j=|| Q (:, j) ||
2Formula 1;
Wherein, Q (:, j) j column vector in the expression Q matrix, || Q (:, j) || the mould value of j column vector in the expression Q matrix, j ∈ i ..., N
t;
If i=1 then calculates the 1st row the~the N in the Q matrix
tThe performance number P of column vector
j, namely be not set among the vectorial E-element of inf;
Step 204, the value of each element among the vectorial E is sorted, the label j that therefrom element of selected value maximum is corresponding, and make i element k (i) among the vectorial k=j;
Step 205, i column vector and k (i) column vector of Q matrix exchanged, i column vector and k (i) column vector of R matrix exchanged;
Step 206, according to the value of the element of the i column vector of Q matrix, calculate the value r of i the elements in a main diagonal in the R matrixes according to formula 2
I, i,
r
I, i=|| Q (:, i) || formula 2;
Wherein, Q (:, i) i column vector in the expression Q matrix, || Q (:, i) || the mould value of i column vector in the expression Q matrix;
The value r of i the elements in a main diagonal in the R matrix that step 207, basis are determined
I, iAnd the value of the element of the current i column vector of Q matrix, calculate the value of each element of i column vector of Q matrixes and the value of replacing current each element of i column vector according to formula 3;
Q (m, i)=Q (m, i)/r
I, iFormula 3;
Wherein, the value of the capable i column element of m in Q (m, i) the expression Q matrix, m ∈ 1 ..., N
r, N
rThe number of expression reception antenna;
Step 208, according to the value of each element of Q matrix, capable according to i among the formula 4 compute matrix R in i+1 be listed as N
tThe value of element of row, and replace among the matrix R current i capable in i+1 be listed as N
tThe value of the element of row;
r
I, n=Q (:, i)
HQ (:, n)/|| Q (:, i) ||
2Formula 4;
Wherein, Q (:, i)
HThe i column vector is carried out the vector behind the conjugate transpose in the expression Q matrix, Q (:, n) n column vector in the expression Q matrix, n ∈ i+1 ..., N
t;
Step 209, capable according to i in the R matrix in i+1 be listed as N
tCurrent i+1 is listed as N in the value of the element of row and the Q matrix
tThe element value of column vector is listed as N according to i+1 in the formula 5 calculating Q matrixes
tThe element value of column vector, and current i+1 is listed as N in the replacement Q matrix
tThe element value of column vector;
Q (:, n)=Q (:, n)-r
I, nQ (:, i) formula 5;
Wherein, Q (:, n) n column vector in the expression Q matrix, r
I, nThe element value of the capable n row of i in the expression R matrix, Q (:, i) i column vector in the expression Q matrix;
Step 210, make i=i+1, and execution in step 202 is to step 209.
Need to prove, the embodiment of the invention is carried out QR to channel matrix and is decomposed, and the value of the elements in a main diagonal of the R matrix that obtains after decomposing according to QR is determined the order of channel matrix, the QR decomposition of wherein channel matrix being carried out can adopt the method for existing QR decomposition to obtain the elements in a main diagonal of R matrix, can adopt also in the embodiment of the invention that the row to channel matrix reorder in the QR decomposable process, satisfy descending so that obtain the value of the elements in a main diagonal of R matrix as far as possible, make the reliability of order of the channel matrix of determining according to the value of the elements in a main diagonal of R matrix higher.
In practical communication system, it is 2,4,8 configuration that base station and terminal adopt the number of transmitting antenna usually, situation than 2 transmitting antennas, in the situation that 4 transmitting antennas and 8 transmitting antennas, adopt the method for the embodiment of the invention to determine that the reduction of complexity of order of channel matrix is more obvious.
Based on same inventive concept, a kind of device of order of definite channel matrix also is provided in the embodiment of the invention, because the principle of this device solves problem is similar to the method for the order of above-mentioned definite channel matrix, so the enforcement of this device can referring to the enforcement of method, repeat part and repeat no more.
As shown in Figure 3, the embodiment of the invention is determined the device of the order of channel matrix, comprising:
Decomposing module 30 is used for channel matrix is carried out the QR resolution process, obtains the value of the elements in a main diagonal of the upper triangular matrix R that obtains after QR decomposes;
Determination module 31 is used for the value according to the elements in a main diagonal of R matrix, determines the order of channel matrix.
In specific implementation process, determination module 31 is determined the order of channel matrix according to the value of the elements in a main diagonal of R matrix and the threshold value of setting.
Concrete, determination module 31 is determined the elements in a main diagonal intermediate value of R matrixes greater than the number of the elements in a main diagonal of the threshold value of setting, and with the number of the elements in a main diagonal the determined order as channel matrix.
Preferably, before determining the order of channel matrix, the value that determination module 31 will obtain the elements in a main diagonal of the R matrix processing of sorting.
In specific implementation process, if the value of the elements in a main diagonal of 31 pairs of R matrixes of determination module is carried out descending, then determination module 31 is from first the elements in a main diagonal, whether the value of the elements in a main diagonal after definite descending is greater than threshold value successively, after the value of determining the elements in a main diagonal that first is not more than threshold value, determine to be positioned at the value of number of all the elements in a main diagonals before this elements in a main diagonal as the order of channel matrix;
If carrying out ascending order, the value of the elements in a main diagonal of 31 pairs of R matrixes of determination module arranges, then determination module 31 is from first the elements in a main diagonal, determine that successively whether the value of the elements in a main diagonal after ascending order is arranged is greater than threshold value, after determining first value greater than the elements in a main diagonal of threshold value, the value of number of determining to comprise this elements in a main diagonal and being positioned at all the elements in a main diagonals after this elements in a main diagonal is as the order of channel matrix.
Preferably, determination module 31 the value with the elements in a main diagonal of R matrix sort process after, can also according to the ratio of the value of adjacent two the elements in a main diagonals and these adjacent two threshold values that the elements in a main diagonal is corresponding of setting after the ordering processing, determine the order of channel matrix.
In specific implementation process, if the value of the elements in a main diagonal of 31 pairs of R matrixes of determination module is carried out descending, then determination module 31 is from first the elements in a main diagonal, determine successively the ratio of the value of the value of current the elements in a main diagonal and adjacent next the elements in a main diagonal, during greater than the threshold value of correspondence, the order of determining channel matrix is the location label of current the elements in a main diagonal at the ratio of determining;
If carrying out ascending order, the value of the elements in a main diagonal of 31 pairs of R matrixes of determination module arranges, then determination module 31 is from last the elements in a main diagonal, determine successively the ratio of the value of the value of current the elements in a main diagonal and an adjacent upper the elements in a main diagonal, at the ratio of determining during greater than the threshold value of correspondence, the order of determining channel matrix is the difference of the location label of the number of the elements in a main diagonal in the R matrix and a upper the elements in a main diagonal.
Preferably, decomposing module 30 specifically is used for: for Q matrix and R matrix i column vector, be listed as N according to current i in the Q matrix
tThe mould value of column vector, to the processing of reordering of current i column vector in Q matrix and the R matrix, i=1 wherein ..., N
t, N
tNumber for transmitting antenna; Be listed as N according to current i in the Q matrix that reorders after processing
tThe value of element in the capable vector of current i in the value of element and the R matrix determines that i is listed as N in the Q matrix in the column vector
tThe value of each element in the column vector, and current i is listed as N in the replacement Q matrix
tThe value of each element in the column vector, and determine in the R matrix value of each element in the capable vector of i, and replace in the R matrix value of each element in the capable vector of current i.
In the specific implementation process, for i column vector in Q matrix and the R matrix, decomposing module 30 determines that current i is listed as N in the Q matrix
tThe mould value that column vector is corresponding; Current i column vector and j column vector in the Q matrix are exchanged processing, and current i column vector and j column vector in the R matrix are exchanged processing, wherein the j column vector is maximum column vector corresponding to mould value in the Q matrix, j ∈ i ..., N
t.
Those skilled in the art should understand, embodiments of the invention can be provided as method, system or computer program.Therefore, the present invention can adopt complete hardware implementation example, complete implement software example or in conjunction with the form of the embodiment of software and hardware aspect.And the present invention can adopt the form of the computer program of implementing in one or more computer-usable storage medium (including but not limited to magnetic disc store, CD-ROM, optical memory etc.) that wherein include computer usable program code.
The present invention is that reference is described according to flow chart and/or the block diagram of method, equipment (system) and the computer program of the embodiment of the invention.Should understand can be by the flow process in each flow process in computer program instructions realization flow figure and/or the block diagram and/or square frame and flow chart and/or the block diagram and/or the combination of square frame.Can provide these computer program instructions to the processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing device producing a machine, so that the instruction of carrying out by the processor of computer or other programmable data processing device produces the device that is used for realizing in the function of flow process of flow chart or a plurality of flow process and/or square frame of block diagram or a plurality of square frame appointments.
These computer program instructions also can be stored in energy vectoring computer or the computer-readable memory of other programmable data processing device with ad hoc fashion work, so that the instruction that is stored in this computer-readable memory produces the manufacture that comprises command device, this command device is realized the function of appointment in flow process of flow chart or a plurality of flow process and/or square frame of block diagram or a plurality of square frame.
These computer program instructions also can be loaded on computer or other programmable data processing device, so that carry out the sequence of operations step producing computer implemented processing at computer or other programmable devices, thereby be provided for realizing the step of the function of appointment in flow process of flow chart or a plurality of flow process and/or square frame of block diagram or a plurality of square frame in the instruction that computer or other programmable devices are carried out.
Although described the preferred embodiments of the present invention, in a single day those skilled in the art get the basic creative concept of cicada, then can make other change and modification to these embodiment.So claims are intended to all changes and the modification that are interpreted as comprising preferred embodiment and fall into the scope of the invention.
The embodiment of the invention is carried out the QR resolution process to channel matrix, and according to the value of the elements in a main diagonal that obtains the upper triangular matrix R that obtains after QR decomposes, determines the order of this channel matrix; Because the QR of channel matrix decomposes than the Eigenvalues Decomposition of the correlation matrix of channel matrix, complexity is lower, thereby has reduced operand; The embodiment of the invention is decomposed the elements in a main diagonal of the R matrix of output by QR, therefore the uncontrollable problem of precision that the characteristic value of existing Eigenvalues Decomposition output therefore can not occur adopts the reliability of order of the channel matrix that the method for the embodiment of the invention determines higher.
Obviously, those skilled in the art can carry out various changes and modification to the present invention and not break away from the spirit and scope of the present invention.Like this, if of the present invention these are revised and modification belongs within the scope of claim of the present invention and equivalent technologies thereof, then the present invention also is intended to comprise these changes and modification interior.