e_plane1.m

patchantenna matlab code

function [E_plane, phi]= e_plane1(fr, width, Leff, h, directdB);

% This function was created by Gazi K ASadullah
% fr= Resonant Frequency in GHz.
% w= Width of the microstrip
% h= Thickness of the microstrip
% Er = Di electric Constant
% The function will calculate the Length of the microstrip Antenna.
% References: 
% Microwave Engineering - Third edition- David Pozar
% Microwave filters, Impedence-matching, Networks, and Coupling structures - Mathaei Young Jones
% Passive and Active filters- Theory and implementations - Wai-kai-Chen.
% Digital Signal Processing using MATLAB- Vinay K Ingle.
% Stripline Circuit design- Harlan Howe
% Antenna THeory: Analysis and Design -Balanis (chapter 14)



r=1;
phi = (-pi/2):(pi/1000):(pi/2);
fr=fr*(1e+9);
c=(2.9986e+8)*1000; %mm/sec
lamda_o=c/fr; %mm
ko=(2*pi)/lamda_o;
const=((j)*width*ko*(exp((-j)*ko*r)))/(pi*r);
B=sin(((ko*h)/2)*cos(phi));
C= (((ko*h)/2)*cos(phi));
%B= (sin(((ko*h)/2)*cos(phi)))/(((ko*h)/2)*cos(phi));
D=B/C;
%E=cos(((ko*Leff)/2)*cos(phi));
E= cos(((ko*Leff)/2)*sin(phi));
%E_plane= 20*log10(abs(const*D*E));
plane1=abs(const*D*E);
plane_max=max(plane1);
E_plane = (20*log10(plane1/plane_max)) + directdB; %normalizing


%e_plane=(E_plane);
%[THETA,RHO] = cart2pol(phi,E_plane);
%[E_plane, phi]= e_plane(fr, width, Leff, h);
%max(E_plane+e_plane_max);
%min(E_plane+e_plane_max);
%polar(phi, E_plane, 'r'), xlabel ('Angle (\phi)'), ylabel('GAIN (dB)'), title ('E-Plane E(\phi)' ); title ('E-Plane E(\phi)' );
polarhg(phi',E_plane','rlim',[-30 10],'rtick',[-30 -20 -10 0 10], 'tstep',60,'color','b'), title ('E-Plane E(\phi) (normalized)' );
plotedit on;