recMicrostripAtenna.m

clear all;
clc;
% Given values
% To calculate W and L use this link:
% https://goo.gl/kYUHc3
% W=11.86*10^-3;
% L=9.06*10^-3;
eps_r = 4.5;
fr = 2.4*10^9;
h = 1.6*10^(-3);
c=3*10^8;
eps_0=8.854*10^(-12);
mu_0= 4*pi*10^-7;
%Calculate W, delta L, L and eps_reff
 W = vpa((c/(2*fr))*sqrt(2/(eps_r+1)))
deltaL = 0.412*h*((eps_r+0.3)*(W/h+0.264))/((eps_r-0.258)*(W/h+0.8));
eps_reff = (eps_r+1)/2+(eps_r-1)*(1+12*h/W)^(-1/2)/2;
 L = vpa( 1/(2*fr*sqrt(eps_reff*mu_0*eps_0))-2*deltaL )
%calculate G12
syms theta t
lambda = c/fr;
% W = 0.03768891807;
% L = 0.02909632537;
k_0=2*pi/lambda;
X=k_0*W;
G_12 = (1/(120*pi.^2))*(sin((k_0*W*cos(theta))/2)/cos(theta))^2*(sin(theta))^3*besselj(0,k_0*L*sin(theta));
G_12 = vpaintegral(G_12,theta,[0 pi])

%calculate G1
S_i = vpaintegral(sin(t)/(t),t,[0 X]);
I_1=double(-2+cos(X)+X*S_i+ sin(X)/X);
G_1=(I_1/(120*pi.^2))

%Calculate Rin
R_in = 1/(2*(G_1+ G_12))

%Calculate the x and y-position of the feed point
x_fed = vpa((L/pi)*asin(sqrt(50/R_in)))
y_fed = vpa((L/pi)*acos(sqrt(50/R_in)))

%