Как рассчитать направленность антенны в Matlab

Question

Я пытаюсь вычислить направленность 20 элементов равномерной линейной антенной решетки в Matlab, но у меня не получилось.

Я уже рассчитал коэффициент массива, нормализовал его, преобразовал его в дБ и построил с помощьюкод ниже.

N = 10;                          % Number of array elements
Amplitude = ones(1,N);           % Excitation amplitude of array elements
d = 0.25 : 0.5 : ((N*0.5)-0.25); % Distance between array elements (0.5 lambda)
Phase = zeros(1,N);              % Excitation phase of array elements
Theta = -90:0.25:90 ;            % Angle

ArrayFactor = 0;

for n = 1 : N
    ArrayFactor = ArrayFactor + (Amplitude(n) * cos(2*pi*d(n)*sin(deg2rad(Theta)) + (deg2rad(Phase(n)))));
end

ArrayFactor = ArrayFactor./max(ArrayFactor); % Normalization
ArrayFactor_db = 20 * log10(abs(ArrayFactor)); % Converting to dB

plot(Theta, ArrayFactor_db);
xlim([-90 90]);
ylim([-80 0]);
set(gca,'XTick',(-80:10:80),'YTick',(-80:10:0));

Мой вопрос: как я могу рассчитать направленность этой антенны, используя ArrayFactor_db и Theta