[发明专利]大规模MIMO阵列天线远场辐射场的精确快速计算方法

说明书

技术领域

本发明属于电磁数值计算领域，具体涉及一种采用迭代散射算法求解的大规模MIMO阵列天线远场辐射的精确快速分析方法。

背景技术

大规模天线阵列系统(即Massive MIMO)被认为是未来5G最具潜力的传输技术，它是现有4G网络中MIMO技术的扩展和延伸。Massive MIMO系统在增加系统容量、提高通信质量和针对复杂环境下设备的通用性方面都有很好的作用。它顺应了当前射频部件更贴近天线的发展趋势，降低维护成本和能源成本，同时进一步提高网络性能和部署灵活性，在第五代移动通信系统中具有重要的作用。Massive MIMO有源阵列天线技术或分布式天线技术将成为5G无线传输关键性技术之一，从理论分析上，这种天线技术有着非常明显的优势，但是，在实际应用中将面临着很多问题。一方面由于无线通信环境非常复杂，存在很多不确定的干扰因素，从而对天线性能产生影响；另一方面，Massive MIMO天线阵列和智能控制部分本身的性能也有待在实际使用中进行验证。为了对Massive MIMO天线的性能进行全面考核，需要研究该类型天线性能的精确快速分析方法。

当阵列规模较大时，可以根据经典的方向图乘积定理，忽略单元间耦合的影响，通过阵因子和单元方向图的乘积快速地计算出阵列天线的方向图。然而，由于阵元间存在互耦效应，结果往往与真实情况存在较大差异。

对于阵列天线辐射方向图的精确计算，通常借助于矩量法、有限元法、时域有限差分等电磁数值算法或者基于这些电磁数值算法的商业电磁仿真软件来实现。但是，由于受单台计算机的计算能力限制，当阵列天线规模过大时，任何全波仿真方法，即使是目前高效率的混合矩量法(Hybrid MoM Solution)，也因内存和计算时间的巨大消耗而无法快速精确地实施天线辐射场的计算。

发明内容

本发明针对现有技术通常借助于矩量法、有限元法、时域有限差分等电磁数值算法或者基于这些电磁数值算法的商业电磁仿真软件来实现。但是，由于受单台计算机的计算能力限制，当阵列天线规模过大时，任何全波仿真方法，即使是目前高效率的混合矩量法(Hybrid MoM Solution)，也因内存和计算时间的巨大消耗而无法快速精确地实施天线辐射场的计算等问题，提出一种大规模MIMO阵列天线远场辐射的精确快速分析方法。

本发明的技术方案是：大规模MIMO阵列天线远场辐射场的精确快速计算方法，该计算方法包括以下步骤：

步骤一：确定M×N元平面阵列天线的阵元间距参数d_x,d_y，其中d_x是平面阵列天线x方向的阵元间距，d_y是平面阵列天线y方向的阵元间距；

步骤二：求解单元天线入射场与散射场关系，得到关系式为：

步骤三：根据单元天线间的互耦特性，在阵列环境条件下选择提取单元的子阵形式和规模，取子阵规模为p×q；

步骤四：对于p×q个单元的天线子阵，求解子阵中单元远场辐射方向图；

步骤五：计算M×N元平面阵列天线远场辐射方向图。

所述的大规模MIMO阵列天线远场辐射场的精确快速计算方法，所述步骤二的具体方法为：

考虑平面波入射，将天线单元的入射场和散射场按照球面矢量波形式展开，则有

其中M_mn、N_mn表示球矢量波函数，散射场的球面矢量波系数可以通过单元天线的边界条件确定；

ψ_mn表示球坐标中标量波函数，表达式如下：

其中，是球面第二类汉克尔函数，是勒让德多项式；

空间总电场可以表示为基于矩量法模型，将天线辐射贴片划分为若干网格单元，在每个网格区域外表面应用理想导体边界条件，即可求得单元天线入射场与散射场之间的系数关系，即有

其中，下标j表示天线辐射贴片剖分后网格的编号。

所述的大规模MIMO阵列天线远场辐射场的精确快速计算方法，所述步骤四的具体方法为：

假设单元天线的个数为p×q，取(除去第i个天线外的)p×q-1个天线的散射场作为第i个天线的入射场，即有以下关系式：

该入射场产生新的散射场，其表达式如下：

将新产生的散射场作为下一次迭代的入射场，则第v次迭代的入射场和散射场的表达式如下：