一种基于特征模理论的天线阵列快速分析方法与流程

[0001]
本发明属于天线阵列快速分析技术领域，特别是一种基于特征模理论的天线阵列快速分析方法。


背景技术：

[0002]
无线通信技术不断在发展，天线是通信系统中必要的一部分，通信技术要求的天线性能也在不断提高，大规模的阵列天线散射与辐射特性是目前的研究热点之一。
[0003]
目前针对多尺度、大规模、宽频带的阵列天线仿真与分析主要还是采用体面积分矩量法来分析，这种方法对天线进行整体剖分，导体面和介质体分别采用面剖分和体剖分，用矩量法求解。该方法适用范围广，理论上可以计算任意曲面、电小尺寸共形微带天线阵列，且计算过程简单，但用该方法计算大尺寸阵列时会产生庞大的未知量，甚至容易产生内存溢出的现象。


技术实现要素：

[0004]
本发明的目地在于提供一种适用于大规模天线阵列快速分析、占用内存低、资源消耗低、计算效率高的天线阵列快速分析方法。
[0005]
实现本发明目的的技术解决方案为：一种基于特征模理论的天线阵列快速分析方法，包括以下步骤：
[0006]
步骤1、使用体面积分矩量法分析天线阵元，根据阻抗矩阵方程推导出广义特征方程；
[0007]
步骤2、求解广义特征方程，得到天线阵元的特征值和特征电流，进而得到天线阵元的模式显著性参数以及模式权重系数，选择天线阵元模式显著性参数的模式中的n个主要模式，使用特征电流组成的特征电流矩阵和权重系数表示实际电流；
[0008]
步骤3、用特征模作为天线阵元的全域基函数和整个天线阵列的子域基函数，对矩量法形成的整个天线阵列进行降阶，完成天线阵列快速分析。
[0009]
本发明与现有技术相比，其显著优点在于：(1)将特征模式作为天线阵元的全域基函数和整个天线的子域基函数，改善了矩阵性态；(2)减少了需要求解的未知量，减少系统资源消耗，提高了计算效率。
附图说明
[0010]
图1是本发明基于特征模理论的天线阵列快速分析方法的流程示意图。
具体实施方式
[0011]
本发明为了解决内存时间资源消耗过多的情况，基于特征模理论来降低传统基于体面积分矩量法分析天线阵列的阻抗矩阵维数。特征模方法属于典型的全域基函数方法，
此类方法的本质是使用低阶局域基函数的线性组合来逼近全域或子域全基函数。从区域上看，用低阶基函数聚合起来逼近全域基函数的定义范围，相当于把原来全域基函数定义区域进行网格划分，在每一个网格上定义低阶基函数。从函数关系上来看，用低阶基函数的线性组合来逼近全域基函数，需要通过采用适当的方法决定基函数的聚合形式。特征模是一种非常灵活的全域或子域全基函数形式，具有全域或子域全基函数的特征，但又不受电磁系统的形状与结构的限制。
[0012]
任意形状物体表面的实际电流可以由若干特征模式线性叠加得到，天线的特征模式只与其结构和材料有关，与单元的位置、外加激励无关，每个模式对应着天线本身固有的属性。对于相同的天线单元，对应的特征模式一样，并且谐振状态的模式占主导地位，只需要少数特征模式就可以很好地描述实际电流。在本发明方法中，将特征模式作为天线阵元的全域基函数和整个天线的子域基函数，不仅能够减少求解未知量，且改善了矩阵性态。
[0013]
本发明一种基于特征模理论的天线阵列快速分析方法，包括以下步骤：
[0014]
步骤1、使用体面积分矩量法分析天线阵元，根据阻抗矩阵方程推导出广义特征方程；
[0015]
步骤2、求解广义特征方程，得到天线阵元的特征值和特征电流，进而得到天线阵元的模式显著性参数以及模式权重系数，选择天线阵元模式显著性参数的模式中的n个主要模式，使用特征电流组成的特征电流矩阵和权重系数表示实际电流；
[0016]
步骤3、用特征模作为天线阵元的全域基函数和整个天线阵列的子域基函数，对矩量法形成的整个天线阵列进行降阶，完成天线阵列快速分析。
[0017]
作为一种具体示例，步骤1所述的使用体面积分矩量法分析天线阵元，根据阻抗矩阵方程推导出广义特征方程，具体如下：
[0018]
步骤1.1、由于天线阵元是金属介质混合结构，所以选取体面积分方程矩量法作为基本方法，当电磁波入射到金属介质混合目标时，根据金属表面与介质的边界条件得：
[0019][0020][0021]
其中e
inc
表示入射场，表示体散射场，表示面散射场，下标tan表示切向分量；e
total
表示电场总场；
[0022]
同时由麦克斯韦方程推导出体面积分方程的基本形式：
[0023][0024][0025]
其中下标s表示金属面，v表示介质体，d(r)和ε(r)分别为电位移矢量和介电常数关于距离r的函数，a
s
和a
v
分别是金属面的矢量位函数和介质体的矢量位函数，φ
s
和φ
v
分别是金属面的矢量位函数和介质体的标量位函数；j为复数单位，ω为角频率，表示梯度；
[0026]
步骤1.2、设定天线阵元的金属表面为s，介质体为v，分别用体电流和面电流表示体散射场和面散射场，并且将面电流j
s
用rwg基函数展开，将体电流j
v
和电位移矢量d
v
用swg
基函数展开，进一步将体面积分方程写为矩阵方程形式：
[0027]
z
·
i＝v
ꢀꢀ
(5)
[0028]
其中，为阻抗矩阵，z
ss
表示源和场都是rwg基函数时对应的阻抗矩阵，z
sv
表示源是swg基函数、场是rwg基函数时对应的阻抗矩阵，z
vs
表示源是rwg基函数、场是swg基函数时对应的阻抗矩阵，z
vv
表示源和场都是swg基函数时对应的阻抗矩阵，为待求电流系数，为电压向量；
[0029]
步骤1.3、由阻抗矩阵z的实部r和虚部x建立广义特征方程：
[0030]
xj
n
＝λ
n
rj
n
ꢀꢀ
(6)
[0031]
式中j
n
是特征电流向量，λ
n
是特征值，下标n表示所对应的第n个模式。
[0032]
作为一种具体示例，步骤2所述的求解广义特征方程，得到天线阵元的特征值和特征电流，进而得到天线阵元的模式显著性参数以及模式权重系数，选择天线阵元模式显著性参数的模式中的n个模式，使用特征电流组成的特征电流矩阵和权重系数表示实际电流，具体如下：
[0033]
步骤2.1、求解式(6)，得到天线阵元的特征值和特征电流，进而得到天线阵元的模式显著性参数以及模式权重系数；
[0034]
步骤2.2、根据特征模理论中模式显著性的定义ms＝1/|1+jλ
n
|，将天线阵元的模式显著性参数的模式归为谐振模式，即主要模式，反之显著性参数较小的模式则归为高次模式；
[0035]
步骤2.3、设单个天线阵元未知量的个数为n，选取m个主要特征模式，m远小于n，将原始求解的特征电流矩阵从n
×
n维，降为m
×
m维，用特征电流组成的特征电流矩阵和权重系数表示实际电流：
[0036]
i＝j0i
α
ꢀꢀ
(7)
[0037]
其中j0＝[j1,j2,
…
,j
m
]是由阵元的特征模组成的列向量矩阵，i
α
＝[α1,α2,
…
,α
m
]是阵元的特征模式对应的权重系数；
[0038]
步骤2.4、将式(7)代入式(5)并左乘使得用特征模作为全域基函数的阻抗矩阵保持对称性：
[0039][0040]
作为一种具体示例，步骤3所述的用特征模作为天线阵元的全域基函数和整个天线阵列的子域基函数，对矩量法形成的整个天线阵列进行降阶，完成天线阵列快速分析，具体如下：
[0041]
步骤3.1、若天线阵列中有p个天线阵元，则运用矩量法形成的整个天线阵列的阻抗矩阵z
p
维数为np
×
np，如下式：
[0042][0043]
步骤3.2、用特征模作为天线阵元的全域基函数和整个天线阵列的子域基函数，对矩量法形成的整个天线阵列的阻抗矩阵进行降阶：
[0044][0045]
其中用于降阶的j由j0构成块对角矩阵，如式(11)所示：
[0046][0047]
降阶后形成的矩阵维数降为mp
×
mp。
[0048]
作为一种具体示例，步骤3中所述的降阶过程中，相同天线阵元的自作用矩阵及相对位置相同的两天线阵元互耦作用矩阵只需计算一次，能够复用；降阶过程中的特征模式对于结构相同的天线阵元也能够复用。
[0049]
作为一种具体示例，步骤3中所述的降阶过程中，对式(9)整个天线阵列形成的阻抗矩阵降阶过程，转换为对z
p
中的子矩阵的循环降阶，每次开辟的内存大小仅为子矩阵大小，处理完立即释放。
[0050]
下面结合附图及具体实施例对本发明做进一步详细介绍。
[0051]
实施例
[0052]
结合图1，本发明一种基于特征模理论的天线阵列快速分析方法，包括以下步骤：
[0053]
步骤1、使用体面积分矩量法分析天线阵元，根据阻抗矩阵方程推导出广义特征方程，具体如下：
[0054]
步骤1.1、首先求解一个天线阵元的特征模式，由于一般天线是金属介质混合结构，所以选取体面积分方程矩量法作为基本方法，当电磁波入射到金属介质混合目标时，根据金属表面与介质的边界条件可得：
[0055][0056][0057]
其中e
inc
表示入射场，表示体散射场，表示面散射场，下标tan表示切向分量；
[0058]
同时由麦克斯韦方程可推导出体面积分方程的基本形式：
[0059][0060]
[0061]
其中下标s表示金属面，v表示介质体，d为电位移矢量，a和φ分别为矢量位和标量位函数；
[0062]
步骤1.2、考虑天线阵元为金属介质复合结构，金属表面为s，介质体为v，分别用体电流和面电流表示体散射场和面散射场，并且将面电流j
s
用rwg基函数展开，将体电流j
v
和电位移矢量d
v
用swg基函数展开，进一步可将体面积分方程写为矩阵方程形式：
[0063]
z
·
i＝v
ꢀꢀ
(5)
[0064]
其中，为阻抗矩阵，z
ss
表示源和场都是rwg基函数时对应的阻抗矩阵，z
sv
表示源是swg基函数、场是rwg基函数时对应的阻抗矩阵，z
vs
表示源是rwg基函数、场是swg基函数时对应的阻抗矩阵，z
vv
表示源和场都是swg基函数时对应的阻抗矩阵，为待求电流系数，为电压向量；
[0065]
步骤1.3、由阻抗矩阵z的实部r和虚部x可建立广义特征方程：
[0066]
xj
n
＝λ
n
rj
n
ꢀꢀ
(6)
[0067]
式中j
n
是特征电流向量，λ
n
是特征值，下标n表示所对应的第n个模式。
[0068]
步骤2、求解广义特征方程，得到天线阵元的特征值和特征电流，进而得到天线阵元的模式显著性参数以及模式权重系数，选择天线阵元模式显著性参数的模式中的n个主要模式，使用特征电流组成的特征电流矩阵和权重系数表示实际电流，具体如下：
[0069]
步骤2.1、求解式(6)，得到天线阵元的特征值和特征电流，进而得到天线阵元的模式显著性参数以及模式权重系数；
[0070]
步骤2.2、根据特征模理论中模式显著性的定义ms＝1/|1+jλ
n
|，将天线阵元模式显著性参数的模式归为谐振模式，即主要模式，反之显著性参数较小的模式则归为高次模式；
[0071]
步骤2.3、设单个天线阵元未知量的个数为n，选取m(m<<n)个主要特征模式，将原始求解的特征电流矩阵从n
×
n维，降为m
×
m维，用特征电流组成的特征电流矩阵和权重系数表示实际电流：
[0072]
i＝j0i
α
ꢀꢀ
(7)
[0073]
其中j0＝[j1,j2,
…
,j
m
]是由阵元的特征模组成的列向量矩阵，i
α
＝[α1,α2,...,α
m
]是阵元的特征模对应的权重系数；
[0074]
步骤2.4、将式(7)代入式(5)并左乘使得用特征模作为全域基函数的阻抗矩阵保持对称性：
[0075][0076]
步骤3、用特征模作为天线阵元的全域基函数和整个天线阵列的子域基函数，对矩量法形成的整个天线阵列进行降阶，具体如下：
[0077]
步骤3.1、若天线阵列中有p个天线阵元，则运用矩量法形成的矩阵维数为np
×
np，如下式：
[0078][0079]
步骤3.2、用特征模作为天线阵元的全域基函数和整个天线阵列的子域基函数，对矩量法形成的整个天线阵列的阻抗矩阵进行降阶：
[0080][0081]
其中用于降阶的j可由j0构成块对角矩阵，如式(11)所示：
[0082][0083]
降阶后形成的矩阵维数降为mp
×
mp。
[0084]
上述降阶过程，其相同天线阵元的自作用矩阵及相对位置相同的两天线阵元互耦作用矩阵只需计算一次，可复用；且降阶过程中的特征模对于结构相同的天线阵元也可进行复用，大大减少计算量。此外，式(9)整个天线阵列形成的阻抗矩阵降阶过程可转换为对z
p
中的子矩阵的循环降阶，每次开辟的内存大小仅为子矩阵大小，处理完立即释放，进一步减少内存消耗。