CN101599789B - 一种信道矩阵秩的测量方法、装置、终端 - Google Patents

Abstract

本发明提供了一种信道矩阵秩的测量方法、装置、终端，其中的方法包括：根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；利用埃尔米特矩阵A构造具有非负实数的特征值的方程组；计算方程组的根，从而得到信道矩阵H的特征值；根据特征值判断信道矩阵H的秩，以使基站根据秩来决定向终端发送的数据流的个数。本发明可用于具有四发四收天线的无线通信***的终端，能够测得四阶矩阵的秩，这样基站可以根据当前信道矩阵秩的情况决定并行发送数据流的个数，从而提高了在同一频率同时传输多个并行数据流的个数。本发明的方法同时也适用于具有两发两收、四发两收和两发四收天线的无线通信***的终端。

Description

一种信道矩阵秩的测量方法、装置、终端

技术领域

本发明涉及无线通信技术领域，尤其涉及一种应用于3GPP LTE(LongTerm Evolution)***终端的信道矩阵秩的测量方法、装置、终端。

背景技术

为提高无线通信***中的数据传输速率，在未来的无线通信***-LTE***中引入了MIMO(Multiple-Input Multiple-Output，多输入多输出)技术，MIMO技术是根据多个天线之间传输信道的不相关特性来提高传输速率的。这就意味着MIMO技术可以在同一频率同时传输多个并行的数据流，并行数据流之间通过信道矩阵的不相关性来区别，因此MIMO***中能同时并行传输的数据流个数取决于信道矩阵中线性无关的向量的个数，也即取决于信道矩阵的秩的大小。

MIMO技术中，发送信号，接收信号和信道矩阵的关系可用如下公式简单表示：

Y_R×1＝H_R×TX_T×1

式中，R表示接收天线个数，T表示发送天线个数。H为信道矩阵，X为发送的信号向量，Y为接收的信号向量。

在无线信道中，信道情况是实时变化的，因此在使用MIMO技术传输的情况下，终端必须时时将信道矩阵秩的信息反馈给基站，这样基站就可以根据当前信道矩阵秩的情况决定并行发送数据流的个数。***中，这一信息被称作RI(Rank Index)。

目前的无线通信***终端一般采用单根发送天线和单根接收天线来进行数据的传输，这种收发方式最多只能传输一个数据流。但是，对于采用多根天线的终端来说，如何准确传输更多的数据流成为亟待解决的问题。

发明内容

本发明的目的是提供一种信道矩阵秩的测量方法、装置、终端，以提高在同一频率同时传输的并行数据流的个数。本发明的信道矩阵秩的测量方法包括步骤：

根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；

利用所述埃尔米特矩阵A构造具有非负实数的特征值的方程组；

计算所述方程组的根，从而得到所述信道矩阵H的特征值；

根据所述特征值判断所述信道矩阵H的秩，以使基站根据所述秩来决定向所述终端发送的数据流的个数。

所述无线通信***的天线配置具体包括：四发四收或两发两收或两发四收或四发两收。

所述无线通信***的天线为四发四收时，所述生成埃尔米特矩阵A的步骤具体为：

A_4×4＝H_4×4 ^H×H_4×4，

其中，H_4×4 ^H为H_4×4的共轭转置矩阵；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述生成埃尔米特矩阵A的步骤具体为：根据所述信道矩阵H计算维数较小的埃尔米特矩阵A，即，如果R＜T，则A_2×2＝H_R×T×H_R×T ^H；

否则，A_2×2＝H_R×T ^H×H_R×T；

其中，

R为接收天线个数；

T为发送天线个数。

所述无线通信***的天线配置为四发四收时，所述构造具有非负实数的特征值的方程组的步骤具体为：利用|λI-A|＝0，构造一个一元四次方程组，其中I为4×4的单位阵，假设A_4×4为：

$A_{4\×4}=\left[\begin{array}{cccc}{r}_{11}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}\end{array}\right]$

则|λI-A|＝0具体为：

$\left|\begin{array}{cccc}{r}_{11}-\mathrm{\λ}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}-\mathrm{\λ}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}-\mathrm{\λ}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}-\mathrm{\λ}\end{array}\right|=0,$

即所述一元四次方程组为：λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0，其中方程各次项a₀～a₃按照常规的行列式计算得到；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述构造具有非负实数的特征值的方程组的步骤具体为：利用|λI-A|＝0，构造一个一元二次方程，其中I为2×2的单位阵，假设A_2×2为：

$A_{2\×2}=\left[\begin{array}{cc}{r}_{11}& r_{12}\\ r_{21}& r_{22}\end{array}\right]$

则|λI-A|＝0具体为：

$\left|\begin{array}{cc}{r}_{11}-\mathrm{\λ}& r_{12}\\ r_{21}& r_{22}-\mathrm{\λ}\end{array}\right|=0,$

即所述一元二次方程为λ²+bλ+c＝0，

其中，b＝-(r₁₁+r₂₂)，c＝(r₁₁r₂₂-r₁₂r₂₁)。

所述无线通信***的天线配置为四发四收时，所述得到所述信道矩阵H的特征值的步骤具体为：

采用费拉里解法对一元四次方程λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0进行求解，得到矩阵A的四个特征值分别为λ₀，λ₁，λ₂，λ₃；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述得到所述信道矩阵H的特征值的步骤具体为：

根据韦达定理解出矩阵A的两个特征值λ′₀，λ′₁，

其中，

${{\mathrm{\λ}}_0}^{\′}=\frac{-b+\sqrt{b^2-4c}}{2}$

${{\mathrm{\λ}}^{\′}}_1=\frac{-b-\sqrt{b^2-4c}}{2}.$

所述无线通信***的天线配置为四发四收时，所述判断所述信道矩阵H的秩的步骤具体为：

将四个特征值λ₀，λ₁，λ₂，λ₃由大到小排序得到λ_i0，λ_i1，λ_i2，λ_i3，按照如下方法进行秩的判断：

如果 $\left(\right|\frac{{\mathrm{\λ}}_{i0}}{{\mathrm{\λ}}_{i1}}|>\mathrm{\Δ}),$ 则所述信道矩阵H的秩为1；

否则，如果 $\left(\right|\frac{{\mathrm{\λ}}_{i0}}{{\mathrm{\λ}}_{i2}}|>\mathrm{\Δ}),$ 则所述信道矩阵H的秩为2；

否则，如果 $\left(\right|\frac{{\mathrm{\λ}}_{i0}}{{\mathrm{\λ}}_{i3}}|>\mathrm{\Δ}),$ 则所述信道矩阵H的秩为3；

否则，所述信道矩阵H的秩为4；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述判断所述信道矩阵H的秩的步骤具体为：

如果 $\left(\right|\frac{{{\mathrm{\λ}}^{\′}}_0}{{{\mathrm{\λ}}^{\′}}_1}|>\mathrm{\Δ}),$ 则所述信道矩阵H的秩为1；

否则，所述信道矩阵H的秩为2；

其中，Δ为所述埃尔米特矩阵A的条件数门限。

利用信道矩阵计算其对应的埃尔米特矩阵，并进一步利用埃尔米特矩阵构造一元四/二次方程，并利用费拉里算法/韦达定理解出该方程的根，即信道矩阵的特征值，最后利用特征值来判断矩阵的秩。

本发明还提供了一种信道矩阵秩的测量装置，包括：

埃尔米特矩阵生成模块，用于根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；

方程生成模块，用于利用所述埃尔米特矩阵A生成具有非负实数的特征值的方程组；

特征值计算模块，用于根据所述方程组计算所述信道矩阵H的特征值；

矩阵秩判断模块，用于根据所述特征值判断信道矩阵H的秩，以使基站根据所述秩来决定向终端发送的数据流的个数。

所述无线通信***的天线配置具体包括：四发四收或两发两收或两发四收或四发两收。

所述无线通信***的天线配置为四发四收时，计算所述信道矩阵H的特征值的方法为费拉里解法；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，计算所述信道矩阵H的特征值的方法为韦达定理。

本发明还提供了一种终端，包括：信道矩阵秩的测量装置，所述测量装置包括：

埃尔米特矩阵生成模块，用于根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；

方程生成模块，用于利用所述埃尔米特矩阵A生成具有非负实数的特征值的方程组；

特征值计算模块，用于根据所述方程组计算所述信道矩阵H的特征值；

矩阵秩判断模块，用于根据所述特征值判断信道矩阵H的秩，以使基站根据所述秩来决定向所述终端发送的数据流的个数。

与现有技术相比，本发明具有以下有益效果：

本发明利用信道矩阵构造出一个特征值均为非负实数的矩阵，并通过对该矩阵求解，得到该矩阵的秩。本发明可用于具有四发四收天线的无线通信***的终端，能够测得四阶矩阵的秩，这样基站可以根据当前信道矩阵秩的情况决定并行发送数据流的个数，从而提高了在同一频率同时传输多个并行数据流的个数。本发明的方法同时也适用于具有两发两收、四发两收和两发四收天线的无线通信***的终端。

附图说明

图1为本发明的信道矩阵秩的测量方法流程图；

图2为本发明的信道矩阵秩的测量装置示意图。

具体实施方式

矩阵的秩等于该矩阵不为0的特征值的个数，因此只要计算出矩阵的特征值就可以计算出矩阵的秩。而对于特征值的计算可以转换为对复数方程一元N次方程|λI-H|＝0求解的问题。值得注意的是对于一元四次和四次以下的一般实数(不是复数)代数方程的求根公式都是存在的，例如对于一元四次方程可采用费拉里解法，对于一元三次方程可采用卡当公式，对于一元二次方程则可以用大家所熟知的韦达定理。对于四次以上的代数方程则没有固定的求根公式。上面所提及的三种方程解法都是针对实数而言的，而信道矩阵却是一个复数矩阵，并不能直接应用上述求根公式。

本发明根据矩阵H与A＝H^HH具有相同的秩，且A的特征值都是非负实数这一原理，可将信道矩阵H秩的计算转化为对实数的一元N次方程|λI-A|＝0根的求解。

下面结合附图对本发明的具体实施方式作进一步详细说明。

参考图1，图1为本发明的信道矩阵秩的测量方法流程图，包括步骤：

步骤1，根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；

步骤2，利用埃尔米特矩阵构造具有非负实数的特征值的方程组；

步骤3，计算方程组的根，从而得到信道矩阵H的特征值；

步骤4，根据特征值判断信道矩阵H的秩，以使基站根据所述秩来决定向终端发送的数据流的个数。

得到信道矩阵H的秩之后，基站就可以根据信道矩阵H的秩决定并行传输的数据流的个数。

下面分别针对LTE***中，MIMO收发天线配置的两发两收、四发两收、两发四收和四发四收四种情况，对上述步骤的实施过程进行详细说明。这里的m发n收是指：m是基站侧的发送天线个数，n是指终端侧的接收天线个数。

步骤1：

对于四发四收这种情况，按如下方式操作：

假设信道矩阵 $H_{4\×4}=\left[\begin{array}{ccc}{h}_{11}& ...& h_{14}\\ ...& ...& ...\\ h_{41}& ...& h_{44}\end{array}\right]\text{,}$ 则A_4×4＝H_4×4 ^H×H_4×4’当然，矩阵A_4×4也可以利用H_4×4×H_4×4 ^H得到，即A_4×4＝H_4×4×H_4×4 ^H。

这里(·)^H表示矩阵(·)的共轭转置。

设矩阵A_4×4由如下元素构成： $A_{4\×4}=\left[\begin{array}{cccc}{r}_{11}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}\end{array}\right]$

对于两发两收、四发两收、两发四收这三种情况，按照如下方式操作：

假设信道矩阵 $H_{R\×T}=\left[\begin{array}{ccc}{h}_{11}& ...& h_{1T}\\ ...& ...& ...\\ h_{R1}& ...& h_{\mathrm{RT}}\end{array}\right],$

如果R＜T，则A_2×2＝H_R×T×H_R×T ^H

否则，A_2×2＝H_R×T ^H×H_R×T

这里(·)^H表示矩阵(·)的共轭转置，min(R，T)＝2。

设矩阵A_2×2由如下元素构成：

$A_{2\×2}=\left[\begin{array}{cc}{r}_{11}& r_{12}\\ r_{21}& r_{22}\end{array}\right]$

步骤2：

对于四发四收这种情况，按如下方式操作：

根据|λI-A|＝0来构造一个一元四次方程组。其中I为4×4的单位阵，具体为：

$\left|\begin{array}{cccc}{r}_{11}-\mathrm{\λ}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}-\mathrm{\λ}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}-\mathrm{\λ}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}-\mathrm{\λ}\end{array}\right|=0,$ 即λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0，其中方程各次项a₀～a₃可按照常规的行列式计算得到。

对于两发两收、四发两收、两发四收这三种情况，按照如下方式操作：

根据|λI-A|＝0来构造一个一元二次方程组。其中I为2×2的单位阵，具体为：

$\left|\begin{array}{cc}{r}_{11}-\mathrm{\λ}& r_{12}\\ r_{21}& r_{22}-\mathrm{\λ}\end{array}\right|=0,$ 即λ²+bλ+c＝0，其中b＝-(r₁₁+r₂₂)，c＝(r₁₁r₂₂-r₁₂r₂₁)

步骤3：

对于四发四收这种情况，按如下方式操作：

采用费拉里解法计算方程λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0的四个根。具体包括以下两个步骤：

第一步，通过解一元三次方程计算配平参数α。

计算如下一元三次方程的根。

$y^3-a_2{y}^2+(a_3{a}_1-{4a}_0)y+a_0({4a}_2-a_3^2)-a_1^2=0$

为表达简便起见，将上式改写成：

y³+b₂y²+b₁y+b₀＝0

其中，b₂＝-a₂，b₁＝(a₃a₁-4a₀)， $b_0=a_0({4a}_2-a_3^2)-a_1^2.$

计算变量p和q：

$p=b_1-\frac{b_2^2}{3}$

$q=b_0+\frac{{2b}_2^3}{27}-\frac{b_2{b}_1}{3}$

计算变量u，v：

$u=\sqrt[3]{\frac{-q+\sqrt{q^3+\frac{{4p}^3}{27}}}{2}}$

$v=\sqrt[3]{\frac{-q-\sqrt{q^3+\frac{{4p}^3}{27}}}{2}}$

u和v均取实数的三次方根。

$y_0=u+v-\frac{b_2}{3},$ $y_1=\mathrm{ue}+{\mathrm{ve}}^2-\frac{b_2}{3},$ $y_2=\mathrm{ve}+{\mathrm{ue}}^2-\frac{b_2}{3}$

其中， $e=e^{j\frac{2\mathrm{\π}}{3}}$

第二步，构造两个一次方程。

将上述一元三次方程的任一实根带入如下两个方程：

${\mathrm{\λ}}^2+\frac{a_3+\sqrt{4y+a_3^2-{4a}_2}}{2}\mathrm{\λ}+(\frac{y}{2}+\frac{a_{3}y-\frac{a_1}{2}}{2\sqrt{4y+a_3^2-{4a}_2}})=0$

${\mathrm{\λ}}^2+\frac{a_3-\sqrt{4y+a_3^2-{4a}_2}}{2}\mathrm{\λ}+(\frac{y}{2}-\frac{a_{3}y-\frac{a_1}{2}}{2\sqrt{4y+a_3^2-{4a}_2}})=0$

利用韦达定理可解出方程的两个一元二次方程的四个根分别为λ₀，λ₁，λ₂，λ₃。

对于两发两收、四发两收、两发四收这三种情况，可直接根据韦达定理解出特征值λ′₀、λ′₁：

${{\mathrm{\λ}}^{\′}}_0=\frac{-b+\sqrt{b^2-4c}}{2}$

${{\mathrm{\λ}}^{\′}}_1=\frac{-b-\sqrt{b^2-4c}}{2}$

步骤4：在实际应用中进行秩的判断时，由于信号不可避免地受到噪声影响，所以不会出现某(几)个特征值恰好为0的情况，这时使用条件数来进行秩的判断往往更有效。

对于四发四收这种情况，按如下方式操作：

将四个特征值由大到小排序得到λ_i0，λ_i1，λ_i2，λ_i3，按照如下方法进行秩的判断：

如果 $\left(\right|\frac{{\mathrm{\λ}}_{i0}}{{\mathrm{\λ}}_{i1}}|>\mathrm{\Δ}),$ 则矩阵A的秩为1，即信道矩阵H的秩为1；

否则，如果 $\left(\right|\frac{{\mathrm{\λ}}_{i0}}{{\mathrm{\λ}}_{i2}}|>\mathrm{\Δ}),$ 则矩阵A的秩为2，即信道矩阵H的秩为2；

否则，如果 $\left(\right|\frac{{\mathrm{\λ}}_{i0}}{{\mathrm{\λ}}_{i3}}|>\mathrm{\Δ}),$ 则矩阵A的秩为3，即信道矩阵H的秩为3；

否则，矩阵A的秩为4，即信道矩阵H的秩为4。

对于两发两收、四发两收、两发四收这三种情况，按照如下方式操作：

如果 $\left(\right|\frac{{{\mathrm{\λ}}^{\′}}_0}{{{\mathrm{\λ}}^{\′}}_1}|>\mathrm{\Δ}),$ 则矩阵A的秩为1，即信道矩阵H的秩为1；

否则，矩阵A的秩为2，即信道矩阵H的秩为2。

这里Δ为矩阵A的条件数门限。条件数越高表示矩阵越不稳定，也就是越接近于奇异矩阵，该门限应该是一个比较大的值。

本发明还提供了一种信道矩阵秩的测量装置，参考图2所示，图2为本发明的信道矩阵秩的测量装置示意图，包括：

埃尔米特矩阵生成模块，用于根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；

方程生成模块，用于利用埃尔米特矩阵A构造具有非负实数的特征值的方程组；

特征值计算模块，用于根据方程组计算信道矩阵H的特征值；

矩阵秩判断模块，根据特征值判断信道矩阵H的秩，以使基站根据所述秩来决定向终端发送的数据流的个数。

每一个模块的具体实现方法分别对应于上文中的步骤1～4。

总之，本发明利用信道矩阵生成具有非负实数的特征根的新矩阵，并根据生成的新矩阵构造了一元四次或一元二次方程，然后根据方程的系数确定信道矩阵的秩。从而使得基站可以根据一元四次或一元二次方程确定的秩来决定终端并行发送数据流的个数，提高了同时传输并行数据流的个数。本发明不仅适用于具有四发四收天线的无线通信***的终端，也同样适用于具有四发两收，两发四收，两发两收的情况，适用范围较广泛。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种信道矩阵秩的测量方法，其特征在于，包括以下步骤：

根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；

利用所述埃尔米特矩阵A构造具有非负实数的特征值的方程组；

计算所述方程组的根，从而得到所述信道矩阵H的特征值；

根据所述特征值判断所述信道矩阵H的秩，以使基站根据所述秩来决定向所述终端发送的数据流的个数；

其中，所述无线通信***的天线配置具体包括：四发四收或两发两收或两发四收或四发两收；

所述无线通信***的天线为四发四收时，所述生成埃尔米特矩阵A的步骤具体为：

A_4×4＝H_4×4 ^H×H_4×4，

其中，H_4×4 ^H为H_4×4的共轭转置矩阵；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述生成埃尔米特矩阵A的步骤具体为：根据所述信道矩阵H计算维数较小的埃尔米特矩阵A，即，如果R＜T，则A_2×2＝H_R×T×H_R×T ^H；

否则，A_2×2＝H_R×T ^H×H_R×T；

其中，

R为接收天线个数；

T为发送天线个数；

所述无线通信***的天线配置为四发四收时，所述构造具有非负实数的特征值的方程组的步骤具体为：利用|λI-A|＝0，构造一个一元四次方程组，其中I为4×4的单位阵，假设A_4×4为：

$A_{4\×4}=\left[\begin{array}{cccc}{r}_{11}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}\end{array}\right]$

则|λI-A|＝0具体为：

$\left|\begin{array}{cccc}{r}_{11}-\mathrm{\λ}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}-\mathrm{\λ}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}-\mathrm{\λ}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}-\mathrm{\λ}\end{array}\right|=0,$

即所述一元四次方程组为：λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0，其中方程各次项a₀～a₃按照常规的行列式计算得到；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述构造具有非负实数的特征值的方程组的步骤具体为：利用|λI-A|＝0，构造一个一元二次方程，其中I为2×2的单位阵，假设A_2×2为：

$A_{2\×2}=\left[\begin{array}{cc}{r}_{11}& r_{12}\\ r_{21}& r_{22}\end{array}\right]$

则|λI-A|＝0具体为：

$\left|\begin{array}{cc}{r}_{11}-\mathrm{\λ}& r_{12}\\ r_{21}& r_{22}-\mathrm{\λ}\end{array}\right|=0,$

即所述一元二次方程为λ²+bλ+c＝0，

其中，b＝-(r₁₁+r₂₂)，c＝(r₁₁r₂₂-r₁₂r₂₁)；

所述无线通信***的天线配置为四发四收时，所述得到所述信道矩阵H的特征值的步骤具体为：

采用费拉里解法对一元四次方程λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0进行求解，得到矩阵A的四个特征值分别为λ₀，λ₁，λ₂，λ₃；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述得到所述信道矩阵H的特征值的步骤具体为：

根据韦达定理解出矩阵A的两个特征值λ′₀，λ′₁，

其中，

${{\mathrm{\λ}}_0}^{\′}=\frac{-b+\sqrt{b^2-4c}}{2}$

${{\mathrm{\λ}}^{\′}}_1=\frac{-b-\sqrt{b^2-4c}}{2};$

所述无线通信***的天线配置为四发四收时，所述判断所述信道矩阵H的秩的步骤具体为：

将四个特征值λ₀，λ₁，λ₂，λ₃由大到小排序得到λ_i0，λ_i1，λ_i2，λ_i3，按照如下方法进行秩的判断：

如果

则所述信道矩阵H的秩为1；

否则，如果

则所述信道矩阵H的秩为2；

否则，如果则所述信道矩阵H的秩为3；

否则，所述信道矩阵H的秩为4；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述判断所述信道矩阵H的秩的步骤具体为：

如果

则所述信道矩阵H的秩为1；

否则，所述信道矩阵H的秩为2；

其中，Δ为所述埃尔米特矩阵A的条件数门限；

利用信道矩阵计算其对应的埃尔米特矩阵，并进一步利用埃尔米特矩阵构造一元四/二次方程，并利用费拉里算法/韦达定理解出该方程的根，即信道矩阵的特征值，最后利用特征值来判断矩阵的秩。

2.一种信道矩阵秩的测量装置，其特征在于，包括：

埃尔米特矩阵生成模块，用于根据无线通信***的终端天线收到的来自基站的相关信息计算得到的信道矩阵H生成埃尔米特矩阵A；

方程生成模块，用于利用所述埃尔米特矩阵A生成具有非负实数的特征值的方程组；

特征值计算模块，用于根据所述方程组计算所述信道矩阵H的特征值；

矩阵秩判断模块，用于根据所述特征值判断信道矩阵H的秩，以使基站根据所述秩来决定向终端发送的数据流的个数；

其中，所述无线通信***的天线配置具体包括：四发四收或两发两收或两发四收或四发两收；

所述无线通信***的天线为四发四收时，所述埃尔米特矩阵生成模块采用下述方法生成埃尔米特矩阵A：

A_4×4＝H_4×4 ^H×H_4×4，

其中，H_4×4 ^H为H_4×4的共轭转置矩阵；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述埃尔米特矩阵生成模块采用下述方法生成埃尔米特矩阵A：

根据所述信道矩阵H计算维数较小的埃尔米特矩阵A，即，如果R＜T，则A_2×2＝H_R×T×H_R×T ^H；

否则，A_2×2＝H_R×T ^H×H_R×T；

其中，

R为接收天线个数；

T为发送天线个数；

所述无线通信***的天线配置为四发四收时，所述方程生成模块采用下述方法生成具有非负实数的特征值的方程组：

利用|λI-A|＝0，构造一个一元四次方程组，其中I为4×4的单位阵，假设A_4×4为：

$A_{4\×4}=\left[\begin{array}{cccc}{r}_{11}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}\end{array}\right]$

则|λI-A|＝0具体为：

$\left|\begin{array}{cccc}{r}_{11}-\mathrm{\λ}& r_{12}& r_{13}& r_{14}\\ r_{21}& r_{22}-\mathrm{\λ}& r_{23}& r_{24}\\ r_{31}& r_{32}& r_{33}-\mathrm{\λ}& r_{34}\\ r_{41}& r_{42}& r_{43}& r_{44}-\mathrm{\λ}\end{array}\right|=0,$

即所述一元四次方程组为：λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0，其中方程各次项a₀～a₃按照常规的行列式计算得到；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述方程生成模块采用下述方法生成具有非负实数的特征值的方程组：

利用|λI-A|＝0，构造一个一元二次方程，其中I为2×2的单位阵，假设A_2×2为：

$A_{2\×2}=\left[\begin{array}{cc}{r}_{11}& r_{12}\\ r_{21}& r_{22}\end{array}\right]$

则|λI-A|＝0具体为：

$\left|\begin{array}{cc}{r}_{11}-\mathrm{\λ}& r_{12}\\ r_{21}& r_{22}-\mathrm{\λ}\end{array}\right|=0,$

即所述一元二次方程为λ²+bλ+c＝0，

其中，b＝-(r₁₁+r₂₂)，c＝(r₁₁r₂₂-r₁₂r₂₁)；

所述无线通信***的天线配置为四发四收时，所述特征值计算模块采用下述方法计算所述信道矩阵H的特征值：

采用费拉里解法对一元四次方程λ⁴+a₃λ³+a₂λ²+a₁λ+a₀＝0进行求解，得到矩阵A的四个特征值分别为λ₀，λ₁，λ₂，λ₃；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述特征值计算模块采用下述方法计算所述信道矩阵H的特征值：

根据韦达定理解出矩阵A的两个特征值λ′₀，λ′₁，

其中，

${{\mathrm{\λ}}_0}^{\′}=\frac{-b+\sqrt{b^2-4c}}{2}$

${{\mathrm{\λ}}^{\′}}_1=\frac{-b-\sqrt{b^2-4c}}{2};$

所述无线通信***的天线配置为四发四收时，所述矩阵秩判断模块采用下述方法判断信道矩阵H的秩：

将四个特征值λ₀，λ₁，λ₂，λ₃由大到小排序得到λ_i0，λ_i1，λ_i2，λ_i3，按照如下方法进行秩的判断：

如果

则所述信道矩阵H的秩为1；

否则，如果

则所述信道矩阵H的秩为2；

否则，如果

则所述信道矩阵H的秩为3；

否则，所述信道矩阵H的秩为4；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，所述矩阵秩判断模块采用下述方法判断信道矩阵H的秩：

如果

则所述信道矩阵H的秩为1；

否则，所述信道矩阵H的秩为2；

其中，Δ为所述埃尔米特矩阵A的条件数门限；

利用信道矩阵计算其对应的埃尔米特矩阵，并进一步利用埃尔米特矩阵构造一元四/二次方程，并利用费拉里算法/韦达定理解出该方程的根，即信道矩阵的特征值，最后利用特征值来判断矩阵的秩。

3.如权利要求2所述的测量装置，其特征在于，所述无线通信***的天线配置为四发四收时，计算所述信道矩阵H的特征值的方法为费拉里解法；

所述无线通信***的天线配置为两发两收或两发四收或四发两收时，计算所述信道矩阵H的特征值的方法为韦达定理。

Images