Ширина ДН и УБЛ

clear all
clc
 
lym=1;
% lym- длина волны
m=7;
n=7;
% где m-число излучателей по оси x, n-число излучателей по оси y
N = m*n;
% N-кол-во излучателей
dx2=0.5*lym;
dy2=0.5*lym;
% где dx,dy шаг решётки по оси x и y
k=2*pi/lym;
% k-волновое число
deltapsix2=0.5;
deltapsiy2=0.5;
ksix2=deltapsix2/(k*dx2);
ksiy2=deltapsiy2/(k*dy2);
% где ksix,ksiy-коэффициенты замедления
theta = (-90:0.58:90)/180*pi;
fi1 = 0; fi2= 90*pi/180;
p3=k*dx2/2; p4=k*dy2/2;
alfa3=p3*((sin(theta).*cos(fi1))-ksix2);
alfa4=p3*((sin(theta).*cos(fi2))-ksix2);
beta3=p4*((sin(theta).*sin(fi1))-ksiy2);
beta4=p4*((sin(theta).*sin(fi2))-ksiy2);
% где alfa и beta -обобщённые аргументы
Fp3 = (sin(m*alfa3)./sin(alfa3));
Fp4 = (sin(m*alfa4)./sin(alfa4)).*(sin(n*beta4)./sin(beta4));
Fp3=abs(Fp3);
Fp4=abs(Fp4);
% где Fp1,2,3,4,5,6 - множитель направленности антенной решётки
F1=cos(theta);
% где F1-комплексная векторная хар-ка направленности изотропного излучателя
Dn3=F1.*Fp3; 
Dn4=F1.*Fp4; 
DnN3 = Dn3/max(max(Dn3));
DnN4 = Dn4/max(max(Dn4));
DnN3=20*log10(DnN3);
DnN4=20*log10(DnN4);
theta = theta/pi*180;
plot(theta,DnN3,theta,DnN4); grid
xlim([-90 90])
ylim([-25 0])
xlabel('theta')
ylabel('F(theta)')
title('Диаграмма направленности')