clear;
clc;
close all;
%% Basic Electromagnetic Parameters
Frequency = 10e9;
Lightspeed = physconst('LightSpeed');
Wavelength = Lightspeed/Frequency;
Wavenumber = 2*pi/Wavelength;
%% Array Parameters
N =6;
X = (1:N)*Wavelength/2;
I = ones(1,N);
alpha = zeros(1,N);
%% ArrayFactor Samping
Ns =1000;% Sampling number
theta = linspace(-60,60,Ns);
E =zeros(1,Ns);
for num = 1:Ns
E(num)=sum(I.*exp(1j*(Wavenumber*X*sind(theta(num))+alpha)));
end
%% plot figure
E_dB = db(E)-max(db(E));
figure()
plot(theta,E_dB,'LineWidth',2);%normalized
xlabel('\theta(\circ)');ylabel('dB');
ylim([-40,0]);
grid on
set(gca,'Fontsize',19)
%%
[hpbw,ang_max] = cal_hpbw_2d(E_dB,theta);