clc

clear all
format long

disp('design rectangular microstrip antenna by transmission line model using

matlab');
disp('============================================================================'
);
er=input('enter the dielectric constant:');
h=input('enter the substrate thickness (in mm):');
f=input('enter the frequency (GHz):');
z=input('enter the input impedence(ohm):');
disp('calculating. Please wait.');
f=f*1e9;

%calculate the Width

Wid=(3e8/(sqrt((er+1)/2)*2*f))*1000; %in mm

%calculate the effective dielectric constant

e_eff=((er+1)/2)+((er-1)/2)*(1+((12*h)/Wid))^-0.5;

%calculate the extension of length L

del_l=(((e_eff+0.3)*((Wid/h)+0.264))/((e_eff-0.258)*((Wid/h)+0.8)))*(0.412*h);
%in mm

%calculate the effective length

l_eff=(3e8/(2*f*sqrt(e_eff)))*1000;

%calculate the actual length

L=l_eff-(2*del_l);

la=(3e8/f)*1000;
k=(2*pi)/la;
x=k*(Wid);
i1=-2+cos(x)+(x*sinint(x))+(sin(x)/x);
g1=i1/(120*pi*pi); %conductance

a=@(th)(((sin((x./2).*cos(th))./cos(th)).^2).*(besselj(0,
(k.*L.*sin(th)))).*(sin(th)).^3);
a1=quad(a,0,pi);
g12=a1/(120*pi*pi); %in siemens
r_in=1/(2*(g1+g12)); %in ohms

inset=(L/pi)*(acos(sqrt(z/r_in))); %in mm

Lg_min=6*h+L; %based on Funit s. Nakar (2004) book

Wg_min=6*h+Wid ;
B=60*pi*pi/(z*sqrt(er));
m1=2*2*B-1;
m=log(m1);
n1=B-1;
n=log(n1);

W=(2*h/pi)*(B-1-m+(((er-1)/(2*er))*(n+(0.39*0.61)/er))); %based on Liang J.

(2004) book

g=(3e8*4.65e-9)/(sqrt(2*e_eff)*f*10^-9); %based on journal by M A Matin

(2010)

disp('RECTANGULAR PATCH')
disp(['The width of the patch (Wp) is:',num2str(Wid),'mm'])
disp(['The length of the patch (Lp) is:',num2str(L),'mm'])
disp(['The inset feed patch (Fi) is:',num2str(inset),'mm'])
disp(['The width of the feed line (wf) is:',num2str(W),'mm'])
disp(['The gap of the feed line (Gpf) is:',num2str(g),'mm'])
disp(['The minimum length of ground plane is:',num2str(Lg_min),'mm'])
disp(['The minimum width of ground plane is:',num2str(Wg_min),'mm'])