一种基于虚源理论的6极化MIMO信道建模方法

摘要

本发明公开了一种基于虚源理论的6极化MIMO信道建模方法，包括以下步骤：计算自由空间信道中6极化信道矩阵H

说明书

技术领域

本发明涉及MIMO无线通信领域，尤其一种基于虚源理论的6 极化MIMO信道建模方法。

背景技术

随着移动多媒体和移动互联网等多种高数据率无线业务的飞速 发展，在无线通信网络中对带宽的需求不断增加。然而无线频谱是 不可再生的资源，因此有限的带宽将严重制约未来高速无线通信系 统的发展。无线通信的研究者不得不对电磁波的信息传输能力进行 重新审视，以获得频谱效率更高的无线传输技术。MIMO技术由于 在频谱利用率和网络覆盖等方面极具潜力而成为多种新一代无线信 息传输系统中不可或缺的核心技术，如蜂窝系统中的3GPP LTE， WiMAX，无线局域网系统中的IEEE 802.11n等。然而，在单极化MIMO 系统中，对于散射丰富的环境，获得独立衰落支路所需的发射、接 收天线阵元间距至少要达到半个波长，而对于多径角度扩展较小的 环境，则需要5到10倍波长甚至更大的阵元间距，因此MIMO系统 在小尺寸的便携收发设备中很难实现。因此近年来人们逐渐把目光 转向可以缩小收发阵孔径的多极化MIMO技术。

一个MIMO系统所能支持的独立子信道数可以表示为系统的自 由度。与单极化MIMO系统中利用收发阵元的空间距离来获得MIMO 自由度不同，多极化MIMO技术通过对电磁波极化特性的利用来进 一步提高无线通信系统的信息传输能力，使无线通信系统的容量在 时域、频域和空域之外又获得了另一个维度。更重要的是多极化 MIMO系统可利用空间共点的多极化发射和接收天线来实现，从而使 MIMO技术可以在一些设备尺寸受限的场合得到应用。

虽然在传统的抗衰落技术中，也考虑通过不同的极化状态来获 得分集增益，但这些系统多限于使用水平和垂直两种极化，而6极化 MIMO技术可以对电磁场中E_x，E_y，E_z和H_x，H_y，Y_z共6个极化状 态进行利用。最近一些研究表明通过发射和接收端对这6种极化状态 的利用可以获得一个6×6的MIMO系统，这个系统的自由度为2到6 之间，主要受环境的散射特性影响。由于多极化MIMO系统的信道 特性对环境中散射体的分布和材料等非常敏感，因此建立符合实际 环境的信道模型对于评估多极化MIMO系统性能有着至关重要的意 义。虽然关于单极化MIMO信道建模方法的研究已趋于成熟，但人 们对多极化MIMO信道建模的原理和方法研究在世界范围内还处在 起步阶段。与单极化MIMO相比，多极化MIMO信道建模的难度更 大，因为要考虑信道中散射体对不同极化状态信号间的耦合作用。 与常用的2极化相比，6极化MIMO系统不但要考虑不同极化的电场 信号间的耦合，还要考虑各种极化的磁场信号间以及电场和磁场信 号间的相互耦合，而这些耦合作用具有复杂的物理机制，因此目前 还没有公认的6极化MIMO信道模型。此外，已有考虑极化的MIMO 信道建模多针对2极化或3极化，而关于更多重极化的MIMO信道模 型还很缺乏。在这些考虑极化的MIMO信道建模中，多采用基于统 计分析的建模方法。这些模型主要可以归纳为两类，一类基于散射 矩阵，另一类基于极化信号的相关性。虽然这种基于统计分析的信 道建模方法可以使我们从宏观上对多极化MIMO信道的统计特性有 所认识，但从微观角度这种信道模型无法从物理上体现多极化MIMO 信道特性的细节。在基于散射矩阵的多极化MIMO信道建模中，一 般认为散射后信号的相位在[0，2π)内均匀分布，而功率的变化则通过 XPD(cross-polarization discrimination)来体现，这些XPD一般被描述 为正态分布的随机变量。因为这些信道模型中一般考虑散射体均匀 分布的丰富散射环境，因此认为经过散射后信号的相位变化均匀分 布，但实际信道中可能具有稀疏多径的情况，散射后信号的幅度和 相位变化与散射体的分布以及收发阵与散射体的相对位置因素等都 有密切的关系，而这些关系很难通过这种基于统计的信道模型来体 现。在基于极化信号相关性的多极化MIMO信道建模中，通过多径的 AOD(angle of departure)和AOA(angle of arrival)的不同分布来描述 不同的信道散射环境，多针对这些角度的分布进行假设，然后得到 它们与不同极化状态信道响应相关系数的关系，用随机相关矩阵模 型来进行信道建模。然而，这种建模方法中没有体现信道散射过程 对极化耦合的作用，也无法直接体现MIMO信道容量与其传播环境 间的关系以及收发阵的结构和位置等参数对MIMO信道特性的影响。 此外，这种建模方法中多针对AOD和AOA的连续、均匀分布情况， 而对强直达信号和稀疏多径等信道情况下的应用较为困难。此外， 在已有多极化MIMO信道分析中多进行远场平面波的假设，而这种 假设会使多极化MIMO信道矩阵的计算受到不可忽略的影响。

发明内容

本发明的目的在于解决6极化MIMO系统信道建模时，自由空间 信道中各种极化状态的电、磁场信号间的耦合以及环境中的散射对 这些耦合作用的影响问题。

为了解决以上问题，提出了一种基于虚源理论的6极化MIMO信 道建模方法，包括以下步骤：

计算自由空间信道中6极化信道矩阵H_los；

确定每个虚源的位置，根据边界条件确定每个虚源对于接收电磁场 各分量的贡献，计算相应于每个虚源的信道矩阵H_i(i＝1，...，N)；

将自由空间中直射信号的贡献和每个虚源的贡献相加，得到总的信 道矩阵。

进一步，作为优选方案，所述自由空间中H_los的计算方法如下：

$H_{\mathrm{los}}=\frac{e^{-j{k}_{0}r}}{4\pi r^2}\left(\begin{array}{cccccc}{a^x}_1& {a^x}_2& {a^x}_3& 0& {a^x}_4& {a^x}_5\\ {a^y}_3& {a^y}_1& {a^y}_2& {a^y}_5& 0& {a^y}_4\\ {a^z}_2& {a^z}_3& {a^z}_1& {a^z}_4& {a^z}_5& 0\\ 0& -{a^x}_4& -{a^x}_5& {a^x}_1& {a^x}_2& {a^x}_3\\ -{a^y}_5& 0& -{a^y}_4& {a^y}_3& {a^y}_1& {a^y}_2\\ -{a^z}_4& -{a^z}_5& 0& {a^z}_2& {a^z}_3& {a^z}_1\end{array}\right)$

其中k₀＝2π/λ，r是发射和接收天线间的距离，

a^x₁＝2cos²θ_x(1+ξ^-1)-sin²θ_x(1+ξ+ξ^-1)

${a^x}_2=\sin 2{\theta }_x\cos {\phi }_x(1+{\xi }^{-1})+\frac{1}{2}\sin 2{\theta }_x\cos {\phi }_x(1+\xi +{\xi }^{-1})$

${a^x}_3=\sin 2{\theta }_x\sin {\phi }_x(1+{\xi }^{-1})+\frac{1}{2}\sin 2{\theta }_x\sin {\phi }_x(1+\xi {+\xi }^{-1})$

a^x₄＝-sinφ_xsinθ_x(1+ξ)

a^x₅＝cosφ_xsinθ_x(1+ξ)

ξ＝jk₀r

其中aⁱ_l，l∈{1，2，...，5}的上标，以及θ_i与φ_i的下标，i∈{x，y，z}表示相 应于x，y和z轴的极化，a^y_l和a^z_l具有与a^x_l相似的表达式，除了其 中的θ_x，φ_x被替换成相应的θ_y，φ_y，θ_z，φ_z；

所述相应于每个虚源的信道矩阵H_i的计算方法为：

可以采用与上式相类似的方法得到，同时根据各个虚源的位置 来确定上式中的各种参数，另外对于不同极化方向的电偶极子和磁 偶极子的虚源方向要满足边界条件。

本发明中将虚源理论应用到6极化MIMO信道建模中。虚源理论 指出，当一个天线发射的电磁波遇到一个均匀、无限大平面时，这 个平面会产生反射，反射信号可以看成由一个关于这个平面对称的 虚拟天线发出，因此，在这个平面以上空间的场可以看成来自真实 的辐射源和来自平面下的虚拟辐射源所产生的电磁场之和。

附图说明

当结合附图考虑时，通过参照下面的详细描述，能够更完整更 好地理解本发明以及容易得知其中许多伴随的优点，但此处所说明 的附图用来提供对本发明的进一步理解，构成本发明的一部分，本 发明的示意性实施例及其说明用于解释本发明，并不构成对本发明 的不当限定，其中：

图1由水平和垂直反射面构成的信道示意图；

图2x轴极化的电偶极子在由水平和垂直PEC反射面构成的信道 中虚源的位置和方向示意图；

图3收发天线高度相等时极化自由度与通信距离的关系 (h_t＝h_r＝2m)示意图；

图4收发天线高度不相等时极化自由度与通信距离的关系 (h_t＝20m，h_r＝2m)示意图。

具体实施方式

以下参照图1-4对本发明的实施例进行说明。

为使上述目的、特征和优点能够更加明显易懂，下面结合附图 和具体实施方式对本发明作进一步详细的说明。

当一个无线信道中存在多个反射面从而产生多个虚源时，我们将 6极化信道矩阵H表示为：

$H=H_{\mathrm{los}}+\underset{i=1}{\overset{N}{\Sigma }}{H}_i---\left(1\right)$

其中H_los来自于自由空间LOS(line-of-sight)信号的贡献，H_i来自 于第i(i＝1，...，N，N为虚源的个数)个虚源信号的贡献，它们均为6×6 的矩阵。

一种基于虚源理论的6极化MIMO信道建模方法，包括以下步骤： 计算自由空间信道中6极化信道矩阵H_los；

确定每个虚源的位置，根据边界条件确定每个虚源对于接收电磁场 各分量的贡献，计算相应于每个虚源的信道矩阵H_i(i＝1，...，N)； 将自由空间中LOS信号的贡献和每个虚源的贡献相加，得到总的信 道矩阵。

现以如图1所示的由一个水平反射面和一个垂直反射面构成的 信道中的6极化MIMO信道矩阵的计算作为一个示例性实施例进行详 细的说明。

为了去除空间分集效应对多极化MIMO系统性能的影响，我们 考虑如图1所示的共点、垂直3个电偶极子和3个磁偶极子构成的收 发阵。对于m或n＝{1，2，3}，信道矩阵的元素h_mn分别代表沿着x，y和 z轴发射或接收的电偶极子；类似的，对于m或n＝{4，5，6}，信道矩 阵的元素h_mn分别代表沿着x，y和z轴发射或接收的磁偶极子。

H_los的计算：

对于一个x轴极化的自由空间中的无限小电偶极子，我们可以得 到球坐标中Maxwell方程的精确解：

$E_r=\left(I_{E}l\right){\eta }_0\frac{\cos \theta }{2\pi r^2}(1+\frac{1}{j{k}_{0}r})e^{-j{k}_{0}r}$

$E_{\theta }=j\left(I_{E}l\right){\eta }_0\frac{k_0\sin \theta }{4\mathrm{\pi r}}(1+\frac{1}{j{k}_{0}r}-\frac{1}{{\left(k_{0}r\right)}^2})e^{-j{k}_{0}r}$

$H_{\phi }=j\left(I_{E}l\right)\frac{k_0\sin \theta }{4\mathrm{\pi r}}(1+\frac{1}{j{k}_{0}r})e^{-j{k}_{0}r}---\left(2\right)$

E_φ＝H_r＝H_θ＝0

通过电磁场的对偶性，也可得到自由空间中x轴极化的无限小磁 偶极子的电磁场：

$H_r=\frac{\left(I_{m}l\right)}{{\eta }_0}\frac{\cos \theta }{2\pi r^2}(1+\frac{1}{j{k}_{0}r})e^{-j{k}_{0}r}$

$H_{\theta }=j\frac{\left(I_{M}l\right)}{{\eta }_0}\frac{k_0\sin \theta }{4\mathrm{\pi r}}(1+\frac{1}{j{k}_{0}r}-\frac{1}{{\left(k_{0}r\right)}^2})e^{-j{k}_{0}r}$

$E_{\phi }=-j\left(I_{M}l\right)\frac{k_0\sin \theta }{4\mathrm{\pi r}}(1+\frac{1}{j{k}_{0}r})e^{-j{k}_{0}r}---\left(3\right)$

E_r＝E_θ＝H_φ＝0

(式(2)和(3)来自于H.A.Haus，and J.R.Melcher，Electromagnetic  Fields and Energy，NJ：Prentice-Hall，1989)

与(2)和(3)类似，我们也可以获得自由空间中沿y和z轴极化的电 偶极子和磁偶极子的电磁场表达式。然而有一点需要注意的是，对 于沿y和z轴极化的情况，如果采用与(2)和(3)相同的表达式，则其中 θ和φ的定义也要进行相应的调整。举个例子，对于沿x轴极化的电 偶极子而言，θ是与x轴正向的夹角，φ是在y-z平面的投影与y 轴正向的夹角，其中是连接收发之间的位置矢量。

这样，如果考虑发射和接收端均为空间共点的3个无限小电偶 极子和磁偶极子天线，对它们在自由空间中球坐标系的电磁场(2)和 (3)进行直角坐标系下的坐标变换，并采用如下的归一化方法：I_Ml＝ η₀(I_El)，其中η₀是自由空间的波阻抗。2)接收磁场分量被乘以η₀。我 们可以得到如下的信道矩阵：

$H_{\mathrm{los}}=\frac{e^{-j{k}_{0}r}}{4\pi r^2}\left(\begin{array}{cccccc}{a^x}_1& {a^x}_2& {a^x}_3& 0& {a^x}_4& {a^x}_5\\ {a^y}_3& {a^y}_1& {a^y}_2& {a^y}_5& 0& {a^y}_4\\ {a^z}_2& {a^z}_3& {a^z}_1& {a^z}_4& {a^z}_5& 0\\ 0& -{a^x}_4& -{a^x}_5& {a^x}_1& {a^x}_2& {a^x}_3\\ -{a^y}_5& 0& -{a^y}_4& {a^y}_3& {a^y}_1& {a^y}_2\\ -{a^z}_4& -{a^z}_5& 0& {a^z}_2& {a^z}_3& {a^z}_1\end{array}\right)---\left(4\right)$

其中k₀＝2π/λ，r是发射和接收天线间的距离。

a^x₁＝2cos²θ_x(1+ξ^-1)-sin²θ_x(1+ξ+ξ^-1)

${a^x}_2=\sin 2{\theta }_x\cos {\phi }_x(1+{\xi }^{-1})+\frac{1}{2}\sin 2{\theta }_x\cos {\phi }_x(1+\xi +{\xi }^{-1})$

${a^x}_3=\sin 2{\theta }_x\sin {\phi }_x(1+{\xi }^{-1})+\frac{1}{2}\sin 2{\theta }_x\sin {\phi }_x(1+\xi {+\xi }^{-1})---\left(5\right)$

a^x₄＝-sinφ_xsinφ_x(1+ξ)

a^x₅＝cosφ_xsinθ_x(1+ξ)

ξ＝jk₀r

式(5)中aⁱ_l，l∈{1，2，...，5}的上标，以及θ_i与φ_i的下标，i∈{x，y，z} 表示相应于x，y和z轴的极化。a^y_l和a^z_l具有与a^x_l相似的表达式， 除了其中的θ_x，φ_x被替换成相应的θ_y，φ_y，θ_z，φ_z。这样我们就获得了 在自由空间LOS信道中的6×6信道矩阵。

虚源位置和方向的确定

在图2中以x轴极化的电偶极子为例，描述了信道中虚源的位置 和方向。在如图1所示的由水平和垂直PEC(perfectly electric  conducting)反射面构成的信道中，为了满足边界条件(在PEC面上电 场的切向分量和磁场的法向分量为0)，对于发射源O_t存在3个虚源， 分别是来自于水平反射面的虚源O^ih_t，来自垂直反射面的虚源O^iv_t，它 们的方向与O_t相反，另外一个虚源是O^iv_t来自于水平反射面的虚源O^il_t(或看成O^ih_t来自于垂直反射面的虚源)，它的方向与O_t相同。这样在 收发之间就存在4个多径信号，这些来自于虚源的信号对信道矩阵的 贡献可以按照前面自由空间中LOS信号的信道矩阵相类似的方法计 算，首先根据这些虚源的位置确定这些多径信号的水平方位角和俯 仰角，然后代入式(4)和(5)中就可以实现。需要注意的是，图2中仅以 x轴极化的电偶极子为例给出了虚源的位置和方向，而对于沿其它方 向极化的电偶极子和磁偶极子所对应的虚源的位置与图2所示相同， 但为了满足边界条件，这些虚源的方向要进行相应的调整。

多极化MIMO自由度计算结果

自由度是衡量MIMO系统性能的一个重要参数，它表示此系统所 能支持的有效子信道数。一般采用如下的定义：

$\mathrm{DOF}={\left(\underset{i=1}{\overset{6}{\Sigma }}{\sigma }_i\right)}^2/\underset{i=1}{\overset{6}{\Sigma }}{{\sigma }_i}^2---\left(6\right)$

其中σ_i，i∈{1，2，...，6}是矩阵H的奇异值。下面利用上述信道建 模方法，我们在图3和图4中对自由空间信道和图1所示的存在水平 和垂直反射面的信道中6极化MIMO系统的自由度随收发阵距离的变 化进行了比较。其中工作频率为800MHz，dl_t＝50m，λ为波长，曲线 1表示自由空间信道，曲线2表示图1所示的存在水平和垂直反射面 信道。

首先，在图3中，我们考虑了发射和接收天线高度相等的情况。 图3(a)中的结果显示，在近场情况下，对于这两种信道，6极化MIMO 系统均可以获得很高的极化自由度，约为5.33。然而，随着收发间距 的增加，6极化MIMO系统的自由度值对于自由空间信道，很快下 降为2，而对于存在水平和垂直反射面的信道，并不会单调下降，在 横坐标位于0到1之间会出现一个接近3的峰，然后随距离的增加极 化自由度在3和1之间振荡。图3(b)中的结果同样显示对于收发阵高 度相等的情况，在近场中，6极化MIMO自由度在这两种信道中基本 相等，而在远场会存在较大差异。图4中显示，角度ζ对多极化自 由度的影响也很大，当h_t＝20m，h_r＝2m时，ζ＝π/2时自由空间信道中6 极化MIMO系统的自由度随距离的增加没什么变化，对于存在水平 和垂直反射面的信道中，在横坐标位于2和3之间会出现一个接近3 的峰。ζ＝π/4时，对于自由空间信道，近场6极化MIMO自由度为2 而远场可以达到3.5，当信道中存在水平和垂直反射面时会使近场的 极化自由度增加至2.5。

虽然以上描述了本发明的具体实施方式，但是本领域的技术人 员应当理解，这些具体实施方式仅是举例说明，本领域的技术人员 在不脱离本发明的原理和实质的情况下，可以对上述方法和系统的 细节进行各种省略、替换和改变。例如，合并上述方法步骤，从而 按照实质相同的方法执行实质相同的功能以实现实质相同的结果则 属于本发明的范围。因此，本发明的范围仅由所附权利要求书限定。