Notes

close all
clear all
clc

%% EMII assignment 2

steps = 1000;
phi = 0;
theta = linspace(-pi,pi,steps);
f = 7.6*10^9;
omega = 2*pi*f;
c = 3*10^8;
k = omega/c;
a = 0.02;
b = 0.01;

r = 1; % Not relevant since the results are normalized
E0 = 1; % Not relevant since the results are normalized

% Question 6
% Plot the radiation pattern in the xz-plane for the center frequency
kx = k*sin(theta)*cos(phi);
Ephi = (1i*k*E0*exp(-1i*k*r)*b*cos(a*kx/2))./(r*a*((pi/a)^2-kx.^2));
EphiN = Ephi./Ephi(steps/2); % Normalization

figure()
polarplot(theta,EphiN.^2)
title('Normalized H-plane radiation pattern')

% Question 7
% Plot the radiation pattern in the yz-plane for the center frequency
phi = pi/2;
theta = linspace(-pi,pi,steps);
kx = k*sin(theta)*cos(phi);
ky = k*sin(theta)*sin(phi);
Etheta = (1i*k*E0*exp(-1i*k*r)*b*sinc(b*ky/2).*cos(a*kx/2).*cos(theta)*sin(phi))./(r*a*((pi/a)^2-kx.^2));
EthetaN = Etheta./Etheta(steps/2); % Normalization

figure()
polarplot(theta,EthetaN.^2)
title('Normalized E-plane radiation pattern')