clc;clear all;
%parametros
n=5;
R=50;
f1=1.8e9;
f2=2.6e9;
fc=(f1+f2)/2;
fbw=(f2-f1)/fc;
Er=4.7;
h=1.6;
%valores g
g=[0.6180,1.6180,2,1.6180,0.6180,1];
%Equações de admitancia
ZoJ1=sqrt((pi*fbw)/(2*g(1)));
ZoJn=(pi*fbw)./(2*sqrt(g(2:n).*g(1:(n-1))));
ZoJf=sqrt((pi*fbw)/(2*g(n)*g(n+1)));
ZoJ=[ZoJ1,ZoJn,ZoJf];
%Impedâncias caracteristicas
ZoE=R*(1+ZoJ+(ZoJ).^2);
ZoO=R*(1-ZoJ+(ZoJ).^2);
%Impedâncias equivalentes a linha unica
ZoSE=ZoE/2;
ZoSO=ZoO/2;
%cálculo de w/h SE
A=(ZoSE/60).*sqrt((Er+1)/2)+((Er-1)/(Er+1))*(0.23+0.11/Er);
B=377*pi./(2*ZoSE*sqrt(Er));
Wh_SE=(8*exp(A))./(exp(2*A)-2);
% Wh_SE=(2/pi).*(B-1-log(2*B-1)+(((Er-1)/(2*Er))*(log(B-1)+0.39-0.61/Er)));
%cálculo de w/h SO
A=(ZoSO/60).*sqrt((Er+1)/2)+(((Er-1)/(Er+1))*(0.23+(0.11/Er)));
B=(377*pi)./(2*ZoSO*sqrt(Er));
Wh_SO=(8*exp(A))./(exp(2*A)-2);
% Wh_SO=(2/pi).*(B-1-log(2*B-1)+(((Er-1)/(2*Er))*(log(B-1)+0.39-0.61/Er)));
%cálculo de s/h
num=cosh((pi/2).*Wh_SE)+cosh((pi/2).*Wh_SO)-2;
den=cosh((pi/2).*Wh_SO)-cosh((pi/2).*Wh_SE);
sh = (2/pi).*acosh(num./den);
%cálculo de w/h
syms wh_n;
wh_n_num=0;
wh = zeros(1,length(sh));
for i=1:length(sh)
sh_n = sh(i);
Wh_SE_n = Wh_SE(i);
g = cosh(pi*sh_n/2);
d = cosh((pi*wh_n)+(pi*sh_n/2));
wh_n_num = eval(solve(Wh_SE_n == (2/pi)*acosh((2.*d-g+1)/(g+1)),wh_n));
wh(i) = wh_n_num(1);
end
%Ereff
Erff=((Er+1)./2)+(((Er-1)./2).*(1./sqrt((1+12.*(1./wh)))))
%Dimmensions in mm
Wmm = wh*h
Smm = sh*h
Lmm=3e8./(4.*fc.*sqrt(Erff))*1000