circularpatchradius.m

来自「matlab编写的计算圆贴片天线贴片半径的程序。用户指定贴片的工作模式」· M 代码 · 共 45 行

%function for designing radius of circular patch antenna given resonant frequency and substrate dielectric constant and height

function [re, r, fr]=CircularPatchRadius(epr, h, nm, f)
%All lentgh are in mm
%f: resonant frequency in GHz
%nm: n indicates angular mode; m indicates radial mode.
%epr: substrate dielectric constant.
%h: substrate height

%fr: resonant frequencies of the first five modes of the designed antenna
%r: physical radius of the designed antenna
%re: equivalent radius of the designed antenna

switch nm
    case 11
        kai=1.84118;
    case 21
        kai=3.05424;
    case 02
        kai=3.83171;
    case 31
        kai=4.20119;
    case 41
        kai=5.31700;
end

re=kai*3*1e2/(f*2*pi*sqrt(epr));
r0=re;
r=PhysicalRadius(epr, h, re, r0);

fr(1)=1.84118*3*1e2/(re*2*pi*sqrt(epr));
fr(2)=3.05424*3*1e2/(re*2*pi*sqrt(epr));
fr(3)=3.83171*3*1e2/(re*2*pi*sqrt(epr));
fr(4)=4.20119*3*1e2/(re*2*pi*sqrt(epr));
fr(5)=5.31700*3*1e2/(re*2*pi*sqrt(epr));

function r=PhysicalRadius(epr, h, re, r0)

r=re/sqrt(1+2*h*(log(r0/(2*h))+(1.41*epr+1.77)+h*(0.268*epr+1.65)/r0)/(pi*r0*epr));
% r=re/sqrt(1+2*h*(log(pi*r0/(2*h))+1.7726)/(pi*r0*epr));

if abs(r-r0)>0.001
    r=PhysicalRadius(epr, h, re, r);
end