ФАР и ДН антенны

clear
clc
i=sqrt(-1);
% мнимая единица
m = 1:1:200;
%вектор используемый для формирования вектора % значений ДН в следующей строке – интервал от 90 
% градусов до -90
gr=-90+1.8*m;
 
for m2=1:200 %Цикл по количеству значений ДН
      tetta = -pi/2+ pi/100*m2;    
      % формирование угла прихода эталонного
      % сигнала
      % вектор элементов АР
            n=1:1:10;
            f=cos(pi*sin(tetta)*(n-1));
            g=sin(pi*sin(tetta)*(n-1));
        % h - вектор комплесных отсчетов эталонного
         % сигнала в элементах АР
               h=f+i*g;
               % вычисление модуля ДН
                 symma=0;
                  for lm=1:10
                    symma = symma+h(lm);
                  end     
               symma1=abs(symma);
               diagr(m2)=symma1;
end
% нахождение максимального значения ДН
max1=max(diagr);
diagr=diagr/max1; 
diagr1=20*log10(diagr);% вывод значений в децибеллах
%subplot(2,2,4);
plot (gr,diagr)
%subplot(2,2,2); plot (gr,diagr1)
%subplot(2,2,[1 3]); polar(gr,diagr)