此示例演示了使用无限阵列分析来模拟单个元素的行为 - 嵌入阵列中的单位单元[1] - [3]。假设阵列具有两个维度的无限程度,位于XY平面中。
这单位细胞指向无限数组中的单个元素。单元单元需要一个接地面。没有地平面的天线需要有反射器作为支撑。每种情况的一个典型例子是由反射器和微带贴片天线支持的偶极子。本例将使用反射器支持的偶极子,并分析在10ghz时的单元阻抗行为。单位细胞将有一个 X 横截面。
FREQ = 10E9;vp = physconst('LightSpeed');λ= vp /频率;ucdx = 0.5 *λ;ucdy = 0.5 *λ;
定义单独的元素创造一个长度略微小的薄脂 并将其作为尺寸反射器的激励器分配 在每一边。
d =偶极;d.length = 0.495 * lambda;d.width = lambda / 160;D.Tilt = 90;d.tiltaxis = [0 1 0];r =反射器;r.exciter = D;r.spacing = lambda / 4;R.groundPlaneLength = UCDX;R.groundPlaneWidth = UCDY; figure show(r)
创建无限阵列,并指定反射镜支持的偶极子作为元素并查看它。
infArray = infiniteArray;infArray。元素= r;infArrayFigure =图;显示(infArray)
通过计算扫描阻抗分析无限阵列的阻抗特性。扫描阻抗是单元单元在单频率下的阻抗随扫描角度的变化。使用无限阵列上的扫描角度属性ScanAzimuth和ScanElevation来定义扫描行为。在这里,我们计算了单平面的扫描阻抗,方位角= 0度,仰角在1度步长内从0到90度变化。
%扫描平面定义阿兹= 0;%E-Planeel = 0:1:90;%海拔%计算和情节scanZ =南(1,元素个数(el));infArray。ScanAzimuth =阿兹;为了i = 1:numel(el)infarray.scanelevation = el(i);scanz(i)=阻抗(infarray,freq);结尾图绘制(el,实际(scanZ), el,图像放大(scanZ),'行宽',2);网格在传奇(“抵抗”那'actance')包含('扫描海拔(DEG)') ylabel (“阻抗(ω\)”)标题(扫描阻抗在az = 'num2str(az)的平面度])
无限阵列分析取决于周期性的绿色功能,包括无限的双重求和。有关此信息的更多信息,请参阅文档页面(InfiniteArray.
)。该双重求和中的术语数对结果的融合产生了影响。增加总结术语的数量,以提高收敛。执行所示的命令将使术语总数增加到101(每个用于负数和正面指数的50个术语,第0项术语为第0项)。
numsummationterms(infarray,50);
术语数量较多的结果分别在3平面,方位角= 0,45和90°中示出了扫描阻抗。2.4 GHz机器上每次扫描平面约为100秒,具有32 GB内存。
Az = [0 45 90];% E D h面加载scanzdata.
E-plane
图绘图(EL,Real(ScanZ50Terms(1,:)),El,Imag(ScanZ50Terms(1,:)),'行宽',2);网格在传奇(“抵抗”那'actance')包含('扫描海拔(DEG)') ylabel (“阻抗(ω\)”)标题(扫描阻抗在az = 'num2str(阿兹(1))的平面度])
d平面
图绘图(EL,Real(ScanZ50Terms(2,:)),El,Imag(ScanZ50Terms(2,:)),'行宽',2);网格在传奇(“抵抗”那'actance')包含('扫描海拔(DEG)') ylabel (“阻抗(ω\)”)标题(扫描阻抗在az = 'num2str(阿兹(2))的平面度])
h面
图绘制(el真实(scanZ50terms (3,:)), el,图像放大(scanZ50terms (3:)),'行宽',2);网格在传奇(“抵抗”那'actance')包含('扫描海拔(DEG)') ylabel (“阻抗(ω\)”)标题(扫描阻抗在az = 'num2str(阿兹(3))的平面度])
将扫描角度固定到一个特定的值,并扫描频率以观察这个单元单元的阻抗行为。
az_scan = 0;el_scan = 45;percent_bw =含量;bw = percent_bw *频率;Fmin = freq - bw/2;Fmax = freq + bw/2;infArray。ScanAzimuth = az_scan;infArray。ScanElevation = el_scan; figure impedance(infArray,linspace(fmin,fmax,51));
使用从无限阵列分析得到的扫描阻抗数据来导出扫描单元图(在有限阵列情况下也称为嵌入/阵列单元图)。如[1]-[4]所示,使用隔离单元图和阻抗来计算它。通过分析由一个无限反射器支持的偶极子,并计算其在10ghz的功率模式和阻抗来实现这一点。
R.groundPlaneLength = Inf;R.groundPlaneWidth = INF;Giso = NaN(Numel(AZ),Numel(EL));Gisodb = NaN(Numel(AZ),Numel(EL));为了i = 1:numel(az)giso(i,:) =模式(r,freq,az(i),el,“类型”那'力量');gisodB(我:)= 10 * log10 (giso(我,:));if (i,:) = if (i,:) - max(i,:);结尾Ziso =阻抗(R,FREQ);
扫描元素模式的计算要求我们定义发电机阻抗。在这里,我们选择它是广泛的扫描阻力。
Rg = 185;Xg = 0;Zg = Rg + 1i*Xg;gs =南(元素个数(az),元素个数(el));gsdB =南(元素个数(az),元素个数(el));为了i = 1:numel(az)gs(i,:)= 4 * rg *真实(ziso)。* giso(i,:) ./(abs(scanz50terms(i,:) + zg))。^ 2;GSDB(i,:) = 10 * log10(gs(i,:));GSDB(i,:) = GSDB(I,:) - Max(GSDB(I,:));结尾数字;绘图(EL,GSDB(1,:),EL,GSDB(2,:),EL,GSDB(3,:),'行宽', 2.0)网格在Axis ([0 90 -20 0]) xlabel(“扫描海拔(度)。”) ylabel ('电源模式(DB)'')标题(strcat ('E-Plane(AZ = 0°)电源模式')) 传奇('az = 0 deg'那'az = 45 deg'那'az = 90度'那“位置”那'最好的')
[1] J.Illen,“扫描偶极子阵列的增益和阻抗变化”,天线和传播的IRE交易,Vol.10,No.5,PP.56-572,1962年9月。
[2] R. C. Hansen,相控阵天线,第7章和第8章,John Wiley & Sons公司,第二版,1998。
[3] R. J. Mailloux,“相控阵天线手册”,Artech House,第二版,2005
[4] W. Stutzman, G. Thiele,“天线理论与设计”,John Wiley & Sons Inc.,第3版,2013。