CN103326755A - Method and device determining rank of channel matrix - Google Patents

Method and device determining rank of channel matrix Download PDF

Info

Publication number
CN103326755A
CN103326755A CN201210078419XA CN201210078419A CN103326755A CN 103326755 A CN103326755 A CN 103326755A CN 201210078419X A CN201210078419X A CN 201210078419XA CN 201210078419 A CN201210078419 A CN 201210078419A CN 103326755 A CN103326755 A CN 103326755A
Authority
CN
China
Prior art keywords
value
elements
matrix
main diagonal
order
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201210078419XA
Other languages
Chinese (zh)
Other versions
CN103326755B (en
Inventor
吴凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Telecommunications Technology CATT
Datang Mobile Communications Equipment Co Ltd
Original Assignee
China Academy of Telecommunications Technology CATT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Academy of Telecommunications Technology CATT filed Critical China Academy of Telecommunications Technology CATT
Priority to CN201210078419.XA priority Critical patent/CN103326755B/en
Publication of CN103326755A publication Critical patent/CN103326755A/en
Application granted granted Critical
Publication of CN103326755B publication Critical patent/CN103326755B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

The invention relates to the technical field of wireless communication, in particular to a method and device determining the rank of a channel matrix. The method and device determining the rank of the channel matrix is used for solving the problems that an existing method determining the rank of the channel matrix is large in operand, and low in reliability of a determined rank. The method determining the rank of the channel matrix comprises the following steps that QR decomposition processing is conducted on the channel matrix and the value of a main diagonal element of an upper rectangular matrix R obtained after the QR decompression is obtained; the rank of the channel matrix is determined according to the value of the main diagonal element of the R matrix. The method determining the rank of the channel matrix carries out QR decompression on the channel matrix and is low in complexity, and therefore operand is reduced; in addition, due to the fact that the rank of channel matrix is determined according to the value of the main diagonal element of the R matrix, the reliability of the determined rank is higher.

Description

A kind of method of order of definite channel matrix and device
Technical field
The present invention relates to wireless communication technology field, particularly a kind of method of order of definite channel matrix and device.
Background technology
In multiple antenna communication, the wireless channel simultaneously number of parallel transmission data flow depends on the number of the vector that the channel matrix neutral line is irrelevant, the i.e. size of the order of channel matrix.When the number of the independent data stream of wireless channel transmission during greater than the order of channel matrix, receiving terminal is difficult to acquisition and detects preferably performance; Otherwise, when the number of the independent data stream of wireless channel transmission during less than the order of channel matrix, might cause the waste of space channel resource.Therefore, to the accurate judgement of the order of channel matrix, be to guarantee that multi-aerial radio communication system can obtain preferably transmission performance, and take full advantage of the important supplementary means of space channel resource.
At present, generally to determine the order of channel matrix, the method is by the Eigenvalues Decomposition to channel correlation matrix, to obtain a plurality of characteristic values of channel correlation matrix based on the method for Eigenvalues Decomposition in employing, pass through again the numerical analysis to a plurality of characteristic values, thereby determine the order of channel matrix.Specifically comprise the steps:
Step 1, (wherein channel matrix H is N according to channel matrix H r* N tThe dimension matrix, N rThe number of expression reception antenna, N tThe number of expression transmitting antenna), calculate channel correlation matrix R Xx, R wherein Xx=H HH;
Step 2, to channel correlation matrix R XxCarry out Eigenvalues Decomposition, obtain M eigenvalue λ 1, λ 2..., λ M, λ wherein 1>λ 2>...>λ M, M=min{N r, N t;
Step 3, to eigenvalue λ 1, λ 2..., λ MNumerical values recited analyze: begin to carry out iteration from i=1, iterative formula is: ratio ii/ λ I+1, wherein, i=1 ..., M;
If ratio i>ratio_th determines that then the order of channel matrix is current i, wherein the threshold value of ratio_th for setting;
If ratio i≤ ratio_th, then i=i+1, and continuation iteration.
In sum, adopt the order of determining channel matrix based on the method for Eigenvalues Decomposition, need to carry out preliminary treatment to channel correlation matrix, increased operand; Because the complexity of the Eigenvalues Decomposition of channel correlation matrix is higher, and owing to whenever determining a characteristic value and all need through iterative operation repeatedly, so that the precision of the characteristic value of obtaining is uncontrollable, thereby make the reliability of order of definite channel matrix low.
Summary of the invention
The embodiment of the invention provides a kind of method and device of order of definite channel matrix, and be used for solving the method that has the order of determining channel matrix now and exist operand large, and the low problem of reliability of the order of determining.
The embodiment of the invention provides a kind of method of order of definite channel matrix, comprising:
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;
According to the value of described the elements in a main diagonal, determine the order of described channel matrix.
The embodiment of the invention provides a kind of device of order of definite channel matrix, comprising:
Decomposing module 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 is used for the value according to described the elements in a main diagonal, determines the order of described channel matrix.
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.
Description of drawings
Fig. 1 is the method flow diagram that the embodiment of the invention is determined the order of channel matrix;
Fig. 2 is the embodiment of the invention is carried out the implementation process of QR resolution process to channel matrix flow chart;
Fig. 3 is the structural representation of the device of the embodiment of the invention order of determining channel matrix.
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
Figure BDA0000145949940000051
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
Figure BDA0000145949940000061
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
Figure BDA0000145949940000071
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.

Claims (18)

1. the method for the order of a definite channel matrix is characterized in that, the method comprises:
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;
According to the value of described the elements in a main diagonal, determine the order of described channel matrix.
2. the method for claim 1 is characterized in that, determines that the order of described channel matrix comprises:
According to the value of described the elements in a main diagonal and the threshold value of setting, determine the order of described channel matrix.
3. method as claimed in claim 2 is characterized in that, determines that the order of described channel matrix comprises:
Determine that described 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 described channel matrix.
4. method as claimed in claim 3 is characterized in that, determines also to comprise before the order of described channel matrix: with the processing of sorting of the value of described the elements in a main diagonal.
5. method as claimed in claim 4 is characterized in that, determines that the order of described channel matrix comprises:
If the value to described the elements in a main diagonal is carried out descending, from first the elements in a main diagonal, whether the value of the elements in a main diagonal after definite descending is greater than described threshold value successively, after the value of determining the elements in a main diagonal that first is not more than described 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 described channel matrix;
If the value of described the elements in a main diagonal is carried out ascending order arranges, 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 described threshold value, after determining first value greater than the elements in a main diagonal of described 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 described channel matrix.
6. the method for claim 1 is characterized in that, determines also to comprise before the order of described channel matrix: with the processing of sorting of the value of described the elements in a main diagonal;
The order of determining described channel matrix also comprises:
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 described channel matrix.
7. method as claimed in claim 6 is characterized in that, determines that the order of described channel matrix comprises:
If the value to described the elements in a main diagonal is carried out descending, 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 described channel matrix is the location label of current the elements in a main diagonal at described ratio;
If the value of described the elements in a main diagonal is carried out ascending order arranges, 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 described ratio during greater than the threshold value of correspondence, the order of determining described 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.
8. such as the arbitrary described method of claim 1~7, it is characterized in that, described channel matrix carried out the QR resolution process, comprising:
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.
9. method as claimed in claim 8 is characterized in that, for i column vector in Q matrix and the R matrix, to the processing of reordering of current i column vector in Q matrix and the R matrix, comprising:
Determine 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.
10. the device of the order of a definite channel matrix is characterized in that, this device comprises:
Decomposing module 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 is used for the value according to described the elements in a main diagonal, determines the order of described channel matrix.
11. device as claimed in claim 10 is characterized in that, described determination module specifically is used for:
According to the value of described the elements in a main diagonal and the threshold value of setting, determine the order of described channel matrix.
12. device as claimed in claim 11 is characterized in that, described determination module specifically is used for:
Determine that described 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; Reach the number of the elements in a main diagonal that will determine as the order of described channel matrix.
13. device as claimed in claim 12 is characterized in that, described determination module also is used for:
Before determining the order of described channel matrix, with the processing of sorting of the value of described the elements in a main diagonal.
14. device as claimed in claim 13 is characterized in that, described determination module specifically is used for:
If the value to described the elements in a main diagonal is carried out descending, from first the elements in a main diagonal, whether the value of the elements in a main diagonal after definite descending is greater than described threshold value successively, after the value of determining the elements in a main diagonal that first is not more than described 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 described channel matrix;
If the value of described the elements in a main diagonal is carried out ascending order arranges, 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 described threshold value, after determining first value greater than the elements in a main diagonal of described 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 described channel matrix.
15. device as claimed in claim 10 is characterized in that, described determination module also is used for:
With the processing of sorting of the value of described the elements in a main diagonal; 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 described channel matrix.
16. device as claimed in claim 15 is characterized in that, described determination module specifically is used for:
If the value to described the elements in a main diagonal is carried out descending, 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 described channel matrix is the location label of current the elements in a main diagonal at described ratio;
If the value of described the elements in a main diagonal is carried out ascending order arranges, 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 described ratio during greater than the threshold value of correspondence, the order of determining described 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.
17. such as the arbitrary described device of claim 10~16, it is characterized in that, described decomposing module specifically is used for:
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.
18. device as claimed in claim 17 is characterized in that, described decomposition specifically is used for:
For i column vector in Q matrix and the R matrix, determine 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.
CN201210078419.XA 2012-03-22 2012-03-22 A kind of method and device of the order determining channel matrix Active CN103326755B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210078419.XA CN103326755B (en) 2012-03-22 2012-03-22 A kind of method and device of the order determining channel matrix

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210078419.XA CN103326755B (en) 2012-03-22 2012-03-22 A kind of method and device of the order determining channel matrix

Publications (2)

Publication Number Publication Date
CN103326755A true CN103326755A (en) 2013-09-25
CN103326755B CN103326755B (en) 2016-12-14

Family

ID=49195314

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210078419.XA Active CN103326755B (en) 2012-03-22 2012-03-22 A kind of method and device of the order determining channel matrix

Country Status (1)

Country Link
CN (1) CN103326755B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101166047A (en) * 2006-10-17 2008-04-23 中国科学院上海微系统与信息技术研究所 Multi-antenna communication system transmitting device for channel geometric mean decomposition, receiving device, transmitting method and receiving method
CN101627672A (en) * 2007-01-05 2010-01-13 高通股份有限公司 Be used for method and apparatus in wireless station deal with data
CN101783697A (en) * 2009-01-15 2010-07-21 北京信威通信技术股份有限公司 Transmit diversity method of OFDMA multiple-antenna system
US20110199888A1 (en) * 2010-02-12 2011-08-18 Fujitsu Limited Signal separation device and signal separating method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101166047A (en) * 2006-10-17 2008-04-23 中国科学院上海微系统与信息技术研究所 Multi-antenna communication system transmitting device for channel geometric mean decomposition, receiving device, transmitting method and receiving method
CN101627672A (en) * 2007-01-05 2010-01-13 高通股份有限公司 Be used for method and apparatus in wireless station deal with data
CN101783697A (en) * 2009-01-15 2010-07-21 北京信威通信技术股份有限公司 Transmit diversity method of OFDMA multiple-antenna system
US20110199888A1 (en) * 2010-02-12 2011-08-18 Fujitsu Limited Signal separation device and signal separating method

Also Published As

Publication number Publication date
CN103326755B (en) 2016-12-14

Similar Documents

Publication Publication Date Title
CN111867049B (en) Positioning method, positioning device and storage medium
CN102404055B (en) Method, system and device for measuring interference
CN104541571B (en) Data transmission method for uplink and device
CN108243451A (en) A kind of method and system positioned in TD-LTE grids based on MR data
CN110351010A (en) A kind of channel blind detection method, method for transmitting signals and relevant device
CN112836806B (en) Data format adjustment method, device, computer equipment and storage medium
CN102256365A (en) Method and device for reducing inter-cell interference
EP4184981A1 (en) Processing method and processing apparatus for saving energy of base station
CN104640222A (en) Pilot frequency scheduling method for multi-input multi-output system and synergetic equipment
CN103249160B (en) A kind of resource allocation methods in lte-a system under CoMP
CN102263579B (en) Method and device for judging whether wireless communication user devices are paired, and method and device for pairing wireless communication user devices
CN108632849A (en) A kind of method of adjustment of antenna-feed parameter, device and equipment
CN105379352A (en) System, device, and method for optimizing network
EP2395782B1 (en) Radio resource manager and method for radio resource allocation of physical uplink control channel
CN107404344A (en) Communication means, the network equipment and terminal device
CN105165049B (en) Communication control unit, communication control method and information processing equipment
CN101127547A (en) Diversity switching
CN114009112A (en) Method, device, equipment and medium for determining control channel detection capability
CN108632088B (en) Method for processing business, device and server
CN105182389A (en) Method for collecting positioning data by using user equipment and user equipment
CN107306175A (en) Information processing method, apparatus and system
CN110971320B (en) Beam quality reporting method, device, equipment and storage medium
WO2017148246A1 (en) Data configuration method and device
CN117014057A (en) Network slice resource allocation method, device and storage medium
CN103326755A (en) Method and device determining rank of channel matrix

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

CP01 Change in the name or title of a patent holder
TR01 Transfer of patent right

Effective date of registration: 20210601

Address after: 100085 1st floor, building 1, yard 5, Shangdi East Road, Haidian District, Beijing

Patentee after: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

TR01 Transfer of patent right