主要求证了一个计算天线阵列导纳的方法

2.求证方法

首先算出Y矩阵。在计算单个半波偶极子的导纳时，对其馈电为1伏（原文中是一个单位，进行过归一化处理），因此yjkij，ij是第j个端口的电流（k1,,N81），由此得到Y矩阵。在此过程中将Ek,储存起来。

然后计算计算任意负载下的输入阻抗和合成辐射场。认为1端口是馈电点，2端口是由N个端口组成的可重构单元，由i=Yv得到

i1y11 i2y21y12v1 （1） Y22v2对Y矩阵进行分割，如上式所示。k1终端口的导纳为yL,k,k1,,N1，所以有

i2=-YLv2，YL是个对角矩阵，yL,kkyL,k，与（1）式联立得到

v2(Y22+YL)y21v1-1 （2）

（3）

-1i1yy(Y+Y)y21v1111222Lyinyin是从1端口看入的可重构单元的输入导纳。v1可以由以下式子得到

NE(,)vi1jEj(,) （4）

可重构单元上的电压可由（2）式得到。

3. 3元阵例子

取波长为200M，半波振子半径为0.3mm，剖分成21节，振元间距离为一个波长。

10.1278-4.42877i 10.1278-4.42877i 10.1278-4.42877iY=9.95196-4.76398i 9.95196-4.76398i 9.95196-4.76398i9.95196-4.76398i 9.95196-4.76398i 9.95196-4.76398i13.9 0YL

0 13.9

单位为ms

 0.3095 - 0.0537iv计算出2

 0.3095 - 0.0537i下面是几组全波仿真结果和网络分析结果的数据比较   E(,)/v 全波仿真结果 95 60 140 160 55 60 50 60 0.4403 0.3613 0.2446 0.1011 网络分析结果 0.4115 0.3489 0.2320 0.0485 结论：辐射大的方向基本比较吻合，辐射小的地方相差较大，网路分析结果略小于全波仿真真结果。

仿真时，阵子之间的距离影响准确性。 实际上5元，9元的结果相差比较大，我加上匹配负载试了一下，结果也相差比较大，所以没有粘贴结果，但是普遍网络分析结果小于全波仿真结果。