本文根据参考资料整理出一套设计内嵌馈线式(Inset-Fed Patch Antenna)的公式,并用Octave描述出公式,只要输入必要的参数如板材厚度、介质常数、工作频率等即可生成了一个微带天线的初始值。然而这个初始值只作参考,在仿真软件内按照初始值画出天线后,还需要 调参。
以下是正文内容。
第一部分:公式
0. 选择基质材料,查资料获得介电常数以及选择合适的基质厚度
基质材料选择常用的RO4350,该材料介电常数为3.66,基质厚度决定为0.256mm。
下图是内嵌馈线式的微带天线形状。
宽度定为W,长度为L,Wf为四分之一阻抗变换器的宽度,y0是四分之一阻抗变换器的嵌入距离。
1. 天线的宽度W
书本公式:
2 计算等效介电常数
3 由于磁场边缘效应造成的长度差
4 计算实际长度L
5 计算微带天线阻抗
计算G1书本提供了两个公式。第一个:
第二个是近似值
上式中的J()函数是第一类贝塞尔函数。J0表示零阶第一类贝塞尔函数。
其表达式是:
而后可以算出输入阻抗了:
这里的正负可以取正(我刚入门书本例子[参考资料2]用正号我就用正号……)
6 馈线的嵌入位置y0计算
而
因此
7 四分之一阻抗转换器的宽度Wf
用[3]计算出微带天线的输入阻抗Zi,输出阻抗为了阻抗匹配,需要置为50欧姆。通过下式计算出馈线的阻抗:
然后用工具Si9000计算出相应阻抗下的线宽。
8 四分之一阻抗转换器的长度
四分之一就是指波长的四分之一。
这样这几个步骤,就可以得到天线的初始值了。但是这样的参数还是需要优化的,这个优化过程,可以在仿真软件里面实现。借用知乎上面的一句话,实验最靠谱,仿真比较靠谱,
理论公式完全不靠谱。
第二部分,Octave代码
%计算得到微带天线的长L和宽W:
%先输入Substrate: 介质常数Er、高度h、工作频率fr、光速v0
%reference book:Antenna Theory-Analysis and Design (3rd Edition).pdf
% 0. configure parameters
Er = 2.2;
fr = 10;%GHz
h=0.1588;%cm
v0=30;%30e8
% 1. Width
W = (v0/2/fr)*sqrt(2/(Er+1));%cm
% 2. the effective dielectric constant
Ereff = (Er+1)/2 + (Er-1)/2/sqrt(1+12*h/W);
% 3. The extended incremental length of the patch delta_L
delta_L = h*0.412* (Ereff+0.3)*(W/h+0.264) / (Ereff-0.258)/(W/h+0.8);%cm
% 4. The actual length L of the patch
L = 30/2/fr/sqrt(Ereff)-2*delta_L;%cm
% 5. The input impedance
lambda0 = v0/fr;
%G1=0.00328;
%k0=sqrt((1-120*lambda0*G1/W)*24/h^2);%ref:Antenna Theory-Analysis and Design (3rd Edition).pdf page826
% 先说明一点,这里由于书本没有写清楚k0的取值,给出了公式唯独缺少这个参数的说明,于是我根据例子,把k0这个未知数算了出来。
k0=2.0427;
G1=W/(120*lambda0)*(1-1/24*(k0*h)^2);%formatation(14-8a)
if W<lambda0
G1= 1/90*(W/lambda0)^2;%f(14-13)
else
G1= 1/120*(W/lambda0);%f(14-13)
end
stepsize=0.01;
t=0:stepsize:pi;
epoch = 300;
G12=0;
for m=0:1:epoch
%J0=(-1)^m/4^m/ factorial(m)^2;
G12_p = (sin(k0*W/2*cos(t))./cos(t)).^2 *(-1)^m/4^m/ factorial(m)^2 *(k0*L)^(2*m) .* sin(t).^(2*m) .*sin(t).^3;
G12 = G12 + sum(G12_p*stepsize)/ 120/pi^2;
end
Rin= 1/2/(G1+G12);
% 6. The position of the inset feed point where the input impedance is 50 ohms
y0=acos(sqrt(50/Rin))/pi*L;%cm
% 7. Z1 = sqrt(Rin*50);
% 8. WL = 3e8/fr/sqrt(Ereff)/4;
以上代码我在Octave上面已经测试过,用书本的例子测试同时24GHz天线的长宽计算结果和Infineon的开发板资料比较接近。
题外话:支持开源社区的发展,支持Octave。今天2018年12月30日发现csdn的代码编辑器已经不支持Matlab了,我也已不敢随意用Matlab。
参考资料:
1. RECTENNA: INSET-FED AND EDGE-FED “PATCH” ANTENNAS WITH RECTIFYING CIRCUIT
2. Antenna Theory-Analysis and Design (3rd Edition).pdf
3. http://www.emtalk.com/mpacalc.php