CN103326755B - A kind of method and device of the order determining channel matrix - Google Patents
A kind of method and device of the order determining channel matrix Download PDFInfo
- Publication number
- CN103326755B CN103326755B CN201210078419.XA CN201210078419A CN103326755B CN 103326755 B CN103326755 B CN 103326755B CN 201210078419 A CN201210078419 A CN 201210078419A CN 103326755 B CN103326755 B CN 103326755B
- Authority
- CN
- China
- Prior art keywords
- value
- matrix
- elements
- 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.)
- Active
Links
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
The present embodiments relate to wireless communication technology field, particularly to the method and device of a kind of order determining channel matrix, for solving the method for the existing order determining channel matrix, to there is operand big, and the problem that the reliability of the order determined is low.The method of the order of the determination channel matrix of the embodiment of the present invention includes: channel matrix carries out QR resolution process, obtains the value of the elements in a main diagonal of the upper triangular matrix R obtained after QR decomposes；And the value of the elements in a main diagonal according to R matrix, determine the order of channel matrix.The embodiment of the present invention carries out QR decomposition to channel matrix, and complexity is relatively low, thus decreases operand, and determines the order of channel matrix according to the value of the elements in a main diagonal of R matrix, makes the reliability of order determined higher.
Description
Technical field
The present invention relates to wireless communication technology field, particularly to the method and device of a kind of order determining channel matrix.
Background technology
In multiple antenna communication, wireless channel the number of simultaneously transmitting data in parallel stream can depend on channel matrix
The size of the order of the number of middle linearly independent vector, i.e. channel matrix.Number when the independent data stream of wireless channel transmission
During more than the order of channel matrix, receiving terminal is difficult to acquisition and preferably detects performance；Otherwise, when the independent data of wireless channel transmission
When the number of stream is less than the order of channel matrix, it is possible to the waste of spatial channel resources can be caused.Therefore, the order to channel matrix
Accurate judgement, be to ensure that multi-aerial radio communication system is obtained in that preferable transmission performance, and make full use of space channel
The important supplementary means of resource.
At present, the general method using feature based value to decompose is to determine the order of channel matrix, and the method is by channel
The Eigenvalues Decomposition of correlation matrix, to obtain multiple eigenvalues of channel correlation matrix, then by the numerical value to multiple eigenvalues
Specificity analysis, so that it is determined that the order of channel matrix.Specifically include following steps:
Step 1, according to channel matrix H, (wherein channel matrix H is Nr×NtDimension matrix, NrRepresent the number of reception antenna,
NtRepresent the number of transmission antenna), calculate channel correlation matrix Rxx, wherein Rxx=HHH；
Step 2, to channel correlation matrix RxxCarry out Eigenvalues Decomposition, obtain M eigenvalue λ1,λ2,…,λM, wherein λ1＞
λ2＞ ... ＞ λM, M=min{Nr,Nt}；
Step 3, to eigenvalue λ1,λ2,…,λMNumerical values recited be analyzed: proceeding by iteration from i=1, iteration is public
Formula is: ratioi=λi/λi+1, wherein, i=1 ..., M；
If ratioi＞ ratio_th, it is determined that the order of channel matrix is current i, wherein ratio_th is the door set
Limit value；
If ratioi≤ ratio_th, then i=i+1, and continue iteration.
In sum, the method using feature based value to decompose determines the order of channel matrix, needs channel correlation matrix
Carry out pretreatment, add operand；Owing to the complexity of the Eigenvalues Decomposition of channel correlation matrix is higher, and due to every true
A fixed eigenvalue is required for operating through successive ignition so that the precision of the eigenvalue of acquisition is uncontrollable, so that determine
The reliability of the order of channel matrix is low.
Summary of the invention
Embodiments provide the method and device of a kind of order determining channel matrix, be used for solving existing determining letter
It is big to there is operand in the method for road rank of matrix, and the problem that the reliability of the order determined is low.
A kind of method embodiments providing order determining channel matrix, including:
Channel matrix is carried out QR resolution process, obtains the elements in a main diagonal of the upper triangular matrix R obtained after QR decomposes
Value；
According to the value of described the elements in a main diagonal, determine the order of described channel matrix.
Embodiments provide the device of a kind of order determining channel matrix, including:
Decomposing module, for channel matrix being carried out QR resolution process, the upper triangular matrix R's obtained after obtaining QR decomposition
The value of the elements in a main diagonal；
Determine module, for the value according to described the elements in a main diagonal, determine the order of described channel matrix.
The embodiment of the present invention carries out QR resolution process to channel matrix, and according to obtaining upper three angular moments obtained after QR decomposes
The value of the elements in a main diagonal of battle array R, determines the order of this channel matrix；Owing to the QR of channel matrix decomposes compared to channel matrix
The Eigenvalues Decomposition of correlation matrix, complexity is relatively low, thus decreases operand；The embodiment of the present invention decomposes output by QR
The elements in a main diagonal of R matrix, asks thus without the precision of the eigenvalue output of existing Eigenvalues Decomposition occur is uncontrollable
Topic, the reliability therefore using the order of the channel matrix that the method for the embodiment of the present invention determines is higher.
Accompanying drawing explanation
Fig. 1 is the method flow diagram that the embodiment of the present invention determines the order of channel matrix；
Fig. 2 is the flow chart that the embodiment of the present invention carries out the implementation process of QR resolution process to channel matrix；
Fig. 3 is the structural representation that the embodiment of the present invention determines the device of the order of channel matrix.
Detailed description of the invention
The present invention carries out QR resolution process by the embodiment of the present invention to channel matrix, and obtains according to obtaining after QR decomposes
The value of the elements in a main diagonal of upper triangular matrix R, determine the order of this channel matrix；Thus solve and existing determine channel matrix
The method operand of order big, and the problem that the reliability of the order determined is low.
In the embodiment of the present invention, channel matrix can be estimated to obtain, with LTE (Long by channel situation carries out channel
Term Evolution, Long Term Evolution) as a example by system, UE (User Equipment, subscriber equipment) can be public by community
Pilot tone (CRS) or channel information measurement pilot tone (CSI RS) carry out channel estimation, obtain in the range of certain physical resource
Channel matrix.
Below in conjunction with Figure of description, the embodiment of the present invention is described in further detail.
As it is shown in figure 1, the method that the embodiment of the present invention determines the order of channel matrix, comprise the following steps:
Step 101, the QR resolution process that carries out channel matrix, the master obtaining the upper triangular matrix R obtained after QR decomposes is right
The value of diagonal element；
Wherein, channel matrix H is Nr×NtDimension matrix, wherein NrRepresent the number of reception antenna, NtRepresent transmission antenna
Number, after channel matrix H is carried out QR decomposition, obtains unitary matrice Q and upper triangular matrix R；
Step 102, according to the value of the elements in a main diagonal in R matrix, determine the order of channel matrix.
In specific implementation process, step 102 determining, the order of channel matrix H includes following two mode:
Mode one: according to the threshold value of the value of the elements in a main diagonal in R matrix Yu setting, determine the order of channel matrix；Its
Middle threshold value is empirical value, can set as required；
In specific implementation process, determine that the order of channel matrix includes:
Determine the R matrix the elements in a main diagonal intermediate value number more than the elements in a main diagonal of the threshold value set；
Using the number of the elements in a main diagonal that determines as the order of channel matrix.
Assuming that threshold value is r_th, channel matrix H is after QR decomposes, and in the R matrix obtained, main diagonal element is ri,i,
Wherein i ∈ 1 ..., Nt, then from the beginning of i=1, judge r one by onei,iValue whether more than threshold value r_th, until i=Nt；If sentencing
Disconnected result is: have N number of main diagonal element must be worth more than threshold value r_th, it is determined that the order of channel matrix H is N；
Preferably, in step 102 before determining the order of channel matrix, also include: by the elements in a main diagonal of R matrix
Value is ranked up processing.
In specific implementation process, if the value of the elements in a main diagonal is carried out descending, then from first leading diagonal
Element starts, and determines that the value of the elements in a main diagonal after descending, whether more than the threshold value set, is determining first successively
After the value of the individual the elements in a main diagonal being not more than threshold value, determine all leading diagonals unit being positioned at before this elements in a main diagonal
The value of the number of element is as the order of channel matrix；
Assume Nt=4, the elements in a main diagonal of the R matrix obtained after step 101 is respectively as follows: r1,1、r2,2、r3,3With
r4,4, by r1,1、r2,2、r3,3And r4,4Value carry out descending, make r1,1、r2,2、r3,3And r4,4Value meet: r1,1＞ r2,2＞
r3,3＞ r4,4, after judging, if r1,1> r_th, r2,2> r_th, r3,3< r_th then need not judging r4,4Value and r_th
Relation, determines that the order of channel matrix H is 2；
If the value of the elements in a main diagonal is carried out ascending order arrangement, then, from the beginning of first the elements in a main diagonal, determine successively
The value of the elements in a main diagonal after ascending order arrangement, whether more than threshold value, is determining first leading diagonal unit more than threshold value
After the value of element, determine include this elements in a main diagonal and be positioned at all the elements in a main diagonals after this elements in a main diagonal
The value of number is as the order of channel matrix；
Assume Nt=4, the elements in a main diagonal of the R matrix obtained after step 101 is respectively as follows: r1,1、r2,2、r3,3With
r4,4, by r1,1、r2,2、r3,3And r4,4Value carry out ascending order arrangement, make r1,1、r2,2、r3,3And r4,4Value meet: r1,1＜ r2,2＜
r3,3＜ r4,4, after judging, if r1,1< r_th, r2,2> r_th, then need not judging r3,3And r4,4The relation of value and r_th,
The order determining channel matrix H is 3.
Mode two, the value of the elements in a main diagonal of R matrix is ranked up process, according to sequence process after adjacent two masters
The ratio of the value of diagonal entry and threshold value corresponding to this adjacent two the elements in a main diagonal of setting, determine channel matrix
Order；
In specific implementation process, if the value of the elements in a main diagonal of R matrix is carried out descending, then from first master
Diagonal entry r1,1Start, determine the value of current the elements in a main diagonal and the value of next adjacent the elements in a main diagonal successively
Ratio, i.e. Ratioi=ri,i/ri+1,i+1, in ratio R atio determinediThe threshold corresponding more than this adjacent two the elements in a main diagonal
Value Ratio_thiTime, determine that the order of channel matrix H is current the elements in a main diagonal ri,iLocation label i, wherein arbitrary neighborhood
All corresponding threshold value Ratio_th of two the elements in a main diagonalsi, this threshold value is empirical value, can set as required；
If it should be noted that all of RatioiThe most corresponding threshold value Ratio_thi, it is determined that channel matrix H
Order is Nt。
Assume Nt=4, the elements in a main diagonal of the R matrix obtained after step 101 is respectively as follows: r1,1、r2,2、r3,3With
r4,4, by r1,1、r2,2、r3,3And r4,4Value carry out descending, make r1,1、r2,2、r3,3And r4,4Value meet: r1,1＞ r2,2＞
r3,3＞ r4,4, adjacent r1,1With r2,2Corresponding threshold value is Ratio_th1, adjacent r2,2With r3,3Corresponding threshold value is Ratio_
th2, adjacent r3,3With r4,4Corresponding threshold value is Ratio_th3, then the order of channel matrix H is determined by procedure below:
Make i=1, calculate Ratio1=r1,1/r2,2；
If Ratio1≤Ratio_th1, then Ratio is calculated2=r2,2/r3,3；If Ratio1＞ Ratio_th1, it is determined that letter
The order of road matrix H is 1；
If Ratio1≤Ratio_th1And Ratio2≤Ratio_th2, then Ratio is calculated3=r3,3/r4,4；If Ratio1≤
Ratio_th1And Ratio2＞ Ratio_th2, it is determined that the order of channel matrix H is 2；
If Ratio1≤Ratio_th1, and Ratio2≤Ratio_th2, and Ratio3≤Ratio_th3, it is determined that channel
The order of matrix H is 4；If Ratio1≤Ratio_th1, and Ratio2≤Ratio_th2, and Ratio3＞ Ratio_th3, it is determined that
The order of channel matrix H is 3.
If the value of the elements in a main diagonal of R matrix is carried out ascending order arrangement, then from last the elements in a main diagonal
Start, determine the value of current the elements in a main diagonal and ratio R atio of the value of an adjacent upper the elements in a main diagonal successivelyi=
ri,i/ri-1,i-1, in ratio R atio determinediMore than corresponding threshold value Ratio_thiTime, determine that the order of channel matrix is R matrix
The number N of middle the elements in a main diagonaltWith a upper the elements in a main diagonal ri-1,i-1The difference of location label i-1, i.e. channel matrix H
Order be: Nt-(i-1)。
In specific implementation process, step 101 carries out QR resolution process to channel matrix and comprises the following steps:
For Q matrix and the i-th column vector in R matrix, according to the i-th row current in Q matrix to NtThe modulus value of column vector is right
Q matrix carries out process of reordering, wherein i=1 with current i-th column vector in R matrix ..., Nt, NtNumber for transmission antenna；
According to the i-th row current in the Q matrix after processing that reorders to NtIn column vector, the value of element is current with R matrix
The value of element in i-th row vector, determines in Q matrix that the i-th row are to NtThe value of each element in column vector, and replace in Q matrix when
Front i-th arranges NtThe value of each element in column vector, and determine in R matrix the value of each element in the i-th row vector, and replace R square
The value of each element in current i-th row vector in Zhen.
In specific implementation process, reorder in step 101 place to the i-th column vector current in Q matrix and R matrix
Reason, comprises the following steps:
Determine that in Q matrix, current i-th row are to NtThe modulus value that column vector is corresponding；
The column vector that the i-th column vector current in Q matrix is corresponding with maximum modulus value is swapped process, and R will be determined
The column vector that column vector position corresponding to modulus value maximum with Q matrix in matrix is identical, and by the i-th row current in R matrix to
Amount swaps process with the column vector determined.
Channel matrix H is carried out the enforcement of QR resolution process by 2 pairs of embodiment of the present invention steps 101 below in conjunction with the accompanying drawings
Journey is described in detail.
Step 201, initializing Q matrix and R matrix, initial Q matrix is channel matrix H, and Initial R matrix is Nt×
NtThe null matrix of dimension, sets a 1 × NtDimensional vector k=[1,1 ... 1]；
Step 202, for the i-th column vector in the i-th column vector in Q matrix and R matrix, it is judged that whether i less than Nt, wherein i
=1 ..., Nt, NtNumber for transmission antenna；
The most then perform step 203；
If it is not, then terminate flow process, and export the elements in a main diagonal of R matrix；
Step 203, according to formula 1, calculate the i-th row～N in Q matrixtPerformance number P of column vectorj, wherein j ∈ i ...,
Nt, by each performance number PjStore in vector E, i.e. E=[P1,P2,…PNt], and make vector E in the 1st individual to (i-1)
The value of element is set to-inf (negative infinite)；
Pj=| | Q (:, j) | |2Formula 1；
Wherein, Q (:, j) represent jth column vector in Q matrix, | | Q (:, j) | | represent the modulus value of jth column vector in Q matrix,
j∈{i,…,Nt}；
If i=1, then calculate the 1st row～N in Q matrixtPerformance number P of column vectorj, i.e. vector E is not set to-
The element of inf；
Step 204, the value of each element in vector E is ranked up, the label that element that therefrom selected value is maximum is corresponding
J, and make i-th element k (the i)=j in vector k；
Step 205, the i-th column vector of Q matrix is swapped with kth (i) column vector, by the i-th column vector of R matrix with
Kth (i) column vector swaps；
Step 206, value according to the element of the i-th column vector of Q matrix, calculate i-th master couple in R matrix according to formula 2
Value r of diagonal elementi,i,
ri,i=| | Q (:, i) | | formula 2；
Wherein, Q (:, i) represent the i-th column vector in Q matrix, | | Q (:, i) | | represent the modulus value of the i-th column vector in Q matrix；
Step 207, according to value r of i-th the elements in a main diagonal in the R matrix determinedi,iAnd Q matrix current i-th row to
The value of the element of amount, calculates the value of the i-th each element of column vector of Q matrix and replaces the current i-th each element of column vector according to formula 3
Value；
Q (m, i)=Q (m, i)/ri,iFormula 3；
Wherein, Q (m, i) value of m row the i-th column element in expression Q matrix, m ∈ 1 ..., Nr, NrRepresent reception antenna
Number；
Step 208, value according to each element of Q matrix, calculate in matrix R i+1 in the i-th row according to formula 4 and arrange Nt
The value of element of row, and replace in matrix R i+1 in current i-th row and arrange NtThe value of the element of row；
ri,n=Q (:, i)HQ(:,n)/||Q(:,i)||2Formula 4；
Wherein, Q (:, i)HVector after the i-th column vector carries out conjugate transpose in expression Q matrix, Q (:, n) represent Q matrix
In the n-th column vector, n ∈ i+1 ..., Nt}；
Step 209, arrange N according to i+1 in the i-th row in R matrixtCurrent i+1 in the value of the element of row and Q matrix
Arrange NtThe element value of column vector, calculates i+1 in Q matrix according to formula 5 and arranges NtThe element value of column vector, and replace Q
In matrix, current i+1 arranges NtThe element value of column vector；
Q (:, n)=Q (:, n)-ri,nQ (:, i) formula 5；
Wherein, Q (:, n) represent the n-th column vector in Q matrix, ri,nThe element value of i-th row the n-th row in expression R matrix, Q (:,
I) the i-th column vector in Q matrix is represented；
Step 210, make i=i+1, and perform step 202 to step 209.
It should be noted that the embodiment of the present invention carries out QR decomposition to channel matrix, and according to the R square obtained after QR decomposition
The value of the elements in a main diagonal of battle array determines the order of channel matrix, and the QR wherein carried out channel matrix decomposes can use existing QR
The method decomposed is to obtain the elements in a main diagonal of R matrix, it would however also be possible to employ right in QR catabolic process in the embodiment of the present invention
The row of channel matrix reorder, so that the value obtaining the elements in a main diagonal of R matrix meets descending as far as possible, make basis
The reliability of the order of the channel matrix that the value of the elements in a main diagonal of R matrix determines is higher.
In practical communication system, base station and terminal generally use the configuration that number is 2,4,8 of transmission antenna, compared to
The situation of 2 transmission antennas, in the case of 4 transmission antennas and 8 transmission antennas, the method using the embodiment of the present invention is true
The reduction of the complexity determining the order of channel matrix becomes apparent from.
Based on same inventive concept, the embodiment of the present invention additionally provides the device of a kind of order determining channel matrix, by
The principle solving problem in this device is similar to the method for the above-mentioned order determining channel matrix, and therefore the enforcement of this device can be joined
The enforcement of square method, repeats no more in place of repetition.
As it is shown on figure 3, the embodiment of the present invention determines the device of the order of channel matrix, including:
Decomposing module 30, for channel matrix carries out QR resolution process, obtains the upper triangular matrix R obtained after QR decomposes
The value of the elements in a main diagonal；
Determine module 31, for the value of the elements in a main diagonal according to R matrix, determine the order of channel matrix.
In specific implementation process, determine the threshold value of the module 31 value according to the elements in a main diagonal of R matrix and setting,
Determine the order of channel matrix.
Concrete, determine that module 31 determines the elements in a main diagonal intermediate value main diagonal angle more than the threshold value set of R matrix
The number of line element, and using the number of the elements in a main diagonal that determines as the order of channel matrix.
Preferably, before determining the order of channel matrix, determine that module 31 will obtain the value of the elements in a main diagonal of R matrix
It is ranked up processing.
In specific implementation process, however, it is determined that module 31 carries out descending to the value of the elements in a main diagonal of R matrix, then
Determine that module 31, from the beginning of first the elements in a main diagonal, determines that the value of the elements in a main diagonal after descending is the biggest successively
In threshold value, after the value determining first the elements in a main diagonal being not more than threshold value, determine and be positioned at this elements in a main diagonal
The value of the number of all the elements in a main diagonals before is as the order of channel matrix；
If it is determined that module 31 carries out ascending order arrangement to the value of the elements in a main diagonal of R matrix, it is determined that module 31 is from first
Individual the elements in a main diagonal starts, and determines that the value of the elements in a main diagonal after ascending order arrangement, whether more than threshold value, is determining successively
After first more than the value of the elements in a main diagonal of threshold value, determine and include this elements in a main diagonal and be positioned at this leading diagonal unit
The value of the number of all the elements in a main diagonals after element is as the order of channel matrix.
Preferably, determine that module 31 is after being ranked up the value of the elements in a main diagonal of R matrix processing, it is also possible to according to
After sequence process, the ratio of the value of adjacent two the elements in a main diagonals and these adjacent two the elements in a main diagonals of setting are corresponding
Threshold value, determines the order of channel matrix.
In specific implementation process, however, it is determined that module 31 carries out descending to the value of the elements in a main diagonal of R matrix, then
Determine that module 31, from the beginning of first the elements in a main diagonal, determines the value of current the elements in a main diagonal and next adjacent master successively
The ratio of the value of diagonal entry, when the ratio determined is more than the threshold value of correspondence, determines that the order of channel matrix is for current main right
The location label of diagonal element；
If it is determined that module 31 carries out ascending order arrangement to the value of the elements in a main diagonal of R matrix, it is determined that module 31 is from finally
One the elements in a main diagonal starts, and determines the value of current the elements in a main diagonal and the value of an adjacent upper the elements in a main diagonal successively
Ratio, when the ratio determined is more than corresponding threshold value, determine that the order of channel matrix is the number of the elements in a main diagonal in R matrix
The difference of the location label of mesh and a upper the elements in a main diagonal.
Preferably, decomposing module 30 specifically for: for the i-th column vector in Q matrix and R matrix, according to current in Q matrix
I-th arranges NtThe modulus value of column vector, reorder process, wherein i=to Q matrix and current i-th column vector in R matrix
1,…,Nt, NtNumber for transmission antenna；According to the i-th row current in the Q matrix after processing that reorders to NtUnit in column vector
The value of element in current i-th row vector in value and the R matrix of element, determines in Q matrix that the i-th row are to NtEach element in column vector
Value, and replace in Q matrix current i-th row to NtThe value of each element in column vector, and determine in R matrix each in the i-th row vector
The value of element, and replace in R matrix the value of each element in current i-th row vector.
In specific implementation process, for the i-th column vector in Q matrix and R matrix, decomposing module 30 determines in Q matrix current
I-th arranges NtThe modulus value that column vector is corresponding；The i-th column vector current in Q matrix and jth column vector are swapped process, and
The i-th column vector current in R matrix and jth column vector are swapped process, and wherein in Q matrix, jth column vector is maximum mould
The column vector that value is corresponding, j ∈ i ..., Nt}。
Those skilled in the art are it should be appreciated that embodiments of the invention can be provided as method, system or computer program
Product.Therefore, the reality in terms of the present invention can use complete hardware embodiment, complete software implementation or combine software and hardware
Execute the form of example.And, the present invention can use at one or more computers wherein including computer usable program code
The upper computer program product implemented of usable storage medium (including but not limited to disk memory, CD-ROM, optical memory etc.)
The form of product.
The present invention is with reference to method, equipment (system) and the flow process of computer program according to embodiments of the present invention
Figure and/or block diagram describe.It should be understood that can the most first-class by computer program instructions flowchart and/or block diagram
Flow process in journey and/or square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided
Instruction arrives the processor of general purpose computer, special-purpose computer, Embedded Processor or other programmable data processing device to produce
A raw machine so that the instruction performed by the processor of computer or other programmable data processing device is produced for real
The device of the function specified in one flow process of flow chart or multiple flow process and/or one square frame of block diagram or multiple square frame now.
These computer program instructions may be alternatively stored in and computer or other programmable data processing device can be guided with spy
Determine in the computer-readable memory that mode works so that the instruction being stored in this computer-readable memory produces and includes referring to
Make the manufacture of device, this command device realize at one flow process of flow chart or multiple flow process and/or one square frame of block diagram or
The function specified in multiple square frames.
These computer program instructions also can be loaded in computer or other programmable data processing device so that at meter
Perform sequence of operations step on calculation machine or other programmable devices to produce computer implemented process, thus at computer or
The instruction performed on other programmable devices provides for realizing at one flow process of flow chart or multiple flow process and/or block diagram one
The step of the function specified in individual square frame or multiple square frame.
Although preferred embodiments of the present invention have been described, but those skilled in the art once know basic creation
Property concept, then can make other change and amendment to these embodiments.So, claims are intended to be construed to include excellent
Select embodiment and fall into all changes and the amendment of the scope of the invention.
The embodiment of the present invention carries out QR resolution process to channel matrix, and according to obtaining upper three angular moments obtained after QR decomposes
The value of the elements in a main diagonal of battle array R, determines the order of this channel matrix；Owing to the QR of channel matrix decomposes compared to channel matrix
The Eigenvalues Decomposition of correlation matrix, complexity is relatively low, thus decreases operand；The embodiment of the present invention decomposes output by QR
The elements in a main diagonal of R matrix, asks thus without the precision of the eigenvalue output of existing Eigenvalues Decomposition occur is uncontrollable
Topic, the reliability therefore using the order of the channel matrix that the method for the embodiment of the present invention determines is higher.
Obviously, those skilled in the art can carry out various change and the modification essence without deviating from the present invention to the present invention
God and scope.So, if these amendments of the present invention and modification belong to the scope of the claims in the present invention and equivalent technologies thereof
Within, then the present invention is also intended to comprise these change and modification.
Claims (14)
1. the method for the order determining channel matrix, it is characterised in that the method includes:
Channel matrix is carried out QR resolution process, obtains the value of the elements in a main diagonal of the upper triangular matrix R obtained after QR decomposes；
According to the value of described the elements in a main diagonal, determine the order of described channel matrix；
Wherein it is determined that the order of described channel matrix, including: according to value and the threshold value of setting of described the elements in a main diagonal, really
The order of fixed described channel matrix；Or
It is ranked up the value of described the elements in a main diagonal processing, the value of adjacent two the elements in a main diagonals after processing according to sequence
Ratio and threshold value corresponding to this adjacent two the elements in a main diagonal of setting, determine the order of described channel matrix.
2. the method for claim 1, it is characterised in that determine that the order of described channel matrix includes:
Determine the described the elements in a main diagonal intermediate value number more than the elements in a main diagonal of the threshold value set；
Using the number of the elements in a main diagonal that determines as the order of described channel matrix.
3. method as claimed in claim 2, it is characterised in that before determining the order of described channel matrix, also include: by described
The value of the elements in a main diagonal is ranked up processing.
4. method as claimed in claim 3, it is characterised in that determine that the order of described channel matrix includes:
If the value of described the elements in a main diagonal is carried out descending, from the beginning of first the elements in a main diagonal, determine fall successively
The value of the elements in a main diagonal after sequence arrangement, whether more than described threshold value, is determining first master being not more than described threshold value
After the value of diagonal entry, determine that the value of the number being positioned at all the elements in a main diagonals before this elements in a main diagonal is as institute
State the order of channel matrix；
If the value of described the elements in a main diagonal is carried out ascending order arrangement, from the beginning of first the elements in a main diagonal, determine liter successively
The value of the elements in a main diagonal after sequence arrangement, whether more than described threshold value, is determining that first master being more than described threshold value is right
After the value of diagonal element, determine and include this elements in a main diagonal and be positioned at all leading diagonals unit after this elements in a main diagonal
The value of the number of element is as the order of described channel matrix.
5. the method for claim 1, it is characterised in that determine that the order of described channel matrix includes:
If the value of described the elements in a main diagonal is carried out descending, from the beginning of first the elements in a main diagonal, determine successively and work as
The value of front the elements in a main diagonal and the ratio of the value of next adjacent the elements in a main diagonal, be more than corresponding threshold value at described ratio
Time, determine the location label that order is current the elements in a main diagonal of described channel matrix；
If the value of described the elements in a main diagonal is carried out ascending order arrangement, from the beginning of last the elements in a main diagonal, determine successively
The value of current the elements in a main diagonal and the ratio of the value of an adjacent upper the elements in a main diagonal, at described ratio more than corresponding threshold
During value, determine that the order of described channel matrix is the number of the elements in a main diagonal and the position of a upper the elements in a main diagonal in R matrix
The difference of label.
6. the method as described in Claims 1 to 5 is arbitrary, it is characterised in that described channel matrix is carried out QR resolution process, bag
Include:
For Q matrix and the i-th column vector in R matrix, according to the i-th row current in Q matrix to NtThe modulus value of column vector, to Q matrix
Current i-th column vector carries out process of reordering, wherein i=1 with R matrix ..., Nt, NtNumber for transmission antenna；
According to the i-th row current in the Q matrix after processing that reorders to NtThe value of element and current i-th row in R matrix in column vector
The value of element in vector, determines in Q matrix that the i-th row are to NtThe value of each element in column vector, and replace in Q matrix current i-th
Arrange NtThe value of each element in column vector, and determine in R matrix the value of each element in the i-th row vector, and replace in R matrix when
The value of each element in front i-th row vector.
7. method as claimed in claim 6, it is characterised in that for the i-th column vector in Q matrix and R matrix, to Q matrix and R
In matrix, current i-th column vector carries out process of reordering, including:
Determine that in Q matrix, current i-th row are to NtThe modulus value that column vector is corresponding；
The i-th column vector current in Q matrix and jth column vector are swapped process, and by the i-th column vector current in R matrix with
Jth column vector swaps process, the column vector that modulus value that wherein in Q matrix, jth column vector is maximum is corresponding, j ∈ i ...,
Nt}。
8. the device of the order determining channel matrix, it is characterised in that this device includes:
Decomposing module, for channel matrix carries out QR resolution process, the master obtaining the upper triangular matrix R obtained after QR decomposes is right
The value of diagonal element；
Determine module, for the value according to described the elements in a main diagonal, determine the order of described channel matrix；
Wherein, described determine that module determines the order of described channel matrix, including: according to value and the setting of described the elements in a main diagonal
Threshold value, determine the order of described channel matrix；Or
It is ranked up the value of described the elements in a main diagonal processing, the value of adjacent two the elements in a main diagonals after processing according to sequence
Ratio and threshold value corresponding to this adjacent two the elements in a main diagonal of setting, determine the order of described channel matrix.
9. device as claimed in claim 8, it is characterised in that described determine module specifically for:
Determine the described the elements in a main diagonal intermediate value number more than the elements in a main diagonal of the threshold value set；And the master that will determine
The number of diagonal entry is as the order of described channel matrix.
10. device as claimed in claim 9, it is characterised in that described determine that module is additionally operable to:
Before determining the order of described channel matrix, it is ranked up the value of described the elements in a main diagonal processing.
11. devices as claimed in claim 10, it is characterised in that described determine module specifically for:
If the value of described the elements in a main diagonal is carried out descending, from the beginning of first the elements in a main diagonal, determine fall successively
The value of the elements in a main diagonal after sequence arrangement, whether more than described threshold value, is determining first master being not more than described threshold value
After the value of diagonal entry, determine that the value of the number being positioned at all the elements in a main diagonals before this elements in a main diagonal is as institute
State the order of channel matrix；
If the value of described the elements in a main diagonal is carried out ascending order arrangement, from the beginning of first the elements in a main diagonal, determine liter successively
The value of the elements in a main diagonal after sequence arrangement, whether more than described threshold value, is determining that first master being more than described threshold value is right
After the value of diagonal element, determine and include this elements in a main diagonal and be positioned at all leading diagonals unit after this elements in a main diagonal
The value of the number of element is as the order of described channel matrix.
12. devices as claimed in claim 9, it is characterised in that described determine module specifically for:
If the value of described the elements in a main diagonal is carried out descending, from the beginning of first the elements in a main diagonal, determine successively and work as
The value of front the elements in a main diagonal and the ratio of the value of next adjacent the elements in a main diagonal, be more than corresponding threshold value at described ratio
Time, determine the location label that order is current the elements in a main diagonal of described channel matrix；
If the value of described the elements in a main diagonal is carried out ascending order arrangement, from the beginning of last the elements in a main diagonal, determine successively
The value of current the elements in a main diagonal and the ratio of the value of an adjacent upper the elements in a main diagonal, at described ratio more than corresponding threshold
During value, determine that the order of described channel matrix is the number of the elements in a main diagonal and the position of a upper the elements in a main diagonal in R matrix
The difference of label.
13. as arbitrary in claim 8～12 as described in device, it is characterised in that described decomposing module specifically for:
For Q matrix and the i-th column vector in R matrix, according to the i-th row current in Q matrix to NtThe modulus value of column vector, to Q matrix
Current i-th column vector carries out process of reordering, wherein i=1 with R matrix ..., Nt, NtNumber for transmission antenna；According to weight
In Q matrix after sequence process, current i-th arranges NtThe value of element and unit in current i-th row vector in R matrix in column vector
The value of element, determines in Q matrix that the i-th row are to NtThe value of each element in column vector, and replace in Q matrix current i-th row to NtRow
The value of each element in vector, and determine in R matrix the value of each element in the i-th row vector, and replace in R matrix current i-th row to
The value of each element in amount.
14. devices as claimed in claim 13, it is characterised in that described decomposition specifically for:
For the i-th column vector in Q matrix and R matrix, determine that in Q matrix, current i-th row are to NtThe modulus value that column vector is corresponding；By Q
In matrix, current i-th column vector and jth column vector swap process, and the i-th column vector current in R matrix and jth are arranged to
Amount swaps process, the column vector that modulus value that wherein in Q matrix, jth column vector is maximum is corresponding, j ∈ i ..., Nt}。
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 CN103326755A (en) | 2013-09-25 |
| CN103326755B true 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 (3)
| 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 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5446959B2 (en) * | 2010-02-12 | 2014-03-19 | 富士通株式会社 | Signal separation device and signal separation method |
-
2012
- 2012-03-22 CN CN201210078419.XA patent/CN103326755B/en active Active
Patent Citations (3)
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103326755A (en) | 2013-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104380780B (en) | Transmission method and communication equipment | |
| CN106301506A (en) | A kind of method and device of codebook subset constraint | |
| CN106559365A (en) | Self-interference removing method and equipment | |
| CN109152005A (en) | A kind of uplink beam indicating means, UE, base station and storage medium | |
| CN110351010A (en) | A kind of channel blind detection method, method for transmitting signals and relevant device | |
| CN103220780A (en) | Indoor positioning method, device and system | |
| CN104915855B (en) | The predictor method and device of order competition for orders rate | |
| CN109152042A (en) | Configuration, recognition methods and the base station of lead code, terminal, readable storage medium storing program for executing | |
| CN101464205B (en) | Test modal analysis method based on reduced-basis method | |
| CN107783028A (en) | A kind of chip enters the control method and system of test pattern | |
| CN103326755B (en) | A kind of method and device of the order determining channel matrix | |
| CN106254050A (en) | Extensive mimo system dynamic pilot allocation method based on large scale information | |
| CN110233715A (en) | A kind of method of acquisition of information, terminal, base station and system | |
| CN102855278B (en) | A kind of emulation mode and system | |
| KR20140002299A (en) | Apparatas and method for selecting device to device communication mode according to signal interference | |
| CN102137484B (en) | Method and device for measuring RSCP (received signal code power) of adjacent cells | |
| CN109144840A (en) | A kind of test method, device and its equipment | |
| CN101854438B (en) | Method for searching for contact persons in mobile terminal, mobile terminal and system | |
| CN104104963B (en) | A kind of vector candidate's list generation method and device | |
| CN106325606A (en) | Signal scanning method and device of infrared touch equipment and touch equipment | |
| WO2016184305A1 (en) | System and method for synthetic trace model | |
| CN113128696A (en) | Distributed machine learning communication optimization method and device, server and terminal equipment | |
| CN1588834A (en) | Frequency deviation evaluation device and method | |
| CN1863408A (en) | Method and system for implementing common frequency measurement in TD-SCDMA system | |
| CN106060918A (en) | Power control method and base station |
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 |