Показать параметры на точках

 %Определение координат, Ам и фазы для плоской ромбической антенной решётки
lym=30;                                 
%длина волны
theta=10.*pi/180;
fi=10.*pi/180;
delx=0.5;
dely=0.5;
% Пъедестал по x и y
Nx=11;                                   
%количество излучателей по оси x
Ny=11;                                   
%количество излучателей по оси y
dx=30;                                   
%расстояние между излучателями по оси x
dy=30;                                  
%расстояние между излучателями по оси y
p1=22.5*pi/180;  
%число задающее дискрет фазирования (delta psi)
k=2*pi/lym;                              
% волновое число
n=Nx;
A=ones(n);
m=(n-1)/2;
A=triu(A,-m);
A=tril(A,m);
A=fliplr(A);
A=triu(A,-m);
A=tril(A,m);
A=reshape(A',1,[]);
 
 
dpsix=(k*dx*sin(theta).*cos(fi))*(-1);          
%фазовый сдвиг между излучателями по оси x
dpsiy=(k*dy*sin(theta).*sin(fi))*(-1);          
%фазовый сдвиг между излучателями по оси y
N=Nx*Ny;                                 
x=0:dx:dx*(Nx-1);
y=0:dy:dy*(Ny-1);
[xe ye]=meshgrid(x,y);
%Координаты x, y и z для каждого излучателя
x=reshape(xe',1,[]); 
y=reshape(ye',1,[]);
z=(0:1:N-1)*0;
% Сдвиг фазы для каждого излучателя
if dpsix==0 ;
   Px1=zeros(1,Nx);
 else
   Px1=0:dpsix:dpsix*(Nx-1);
 end
 if dpsiy==0
   Py1=zeros(1,Ny);
 else 
   Py1=0:dpsiy:dpsiy*(Ny-1);
 end
 [P1 P2]=meshgrid(Px1,Py1);
 Px1=reshape(P1',1,[]);
 Py1=reshape(P2',1,[]);
 K1=(Px1+Py1);
 Px = round((K1/p1)+0.5).*p1;
%определение амплитуды тока на элементах
dx1=0:dx:(dx*(Nx-1));%координаты х элементов решетки
dy1=0:dy:(dy*(Ny-1));%координаты у элементов решетки
dx0=mean(dx1);%координаты центра решетки
dy0=mean(dy1);   
[dx dy]=meshgrid(dx1,dy1);%координаты элементов решетки
flag==1;
ar1=((delx+(1-delx).*cos(pi.*(dx-dx0)./(max(dx1)-min(dx1)+eps))).*(dely+(1-dely).*cos(pi.*(dy-dy0)./(max(dy1)-min(dy1)+eps))));%косинус на пьедестале
ar2=((delx+(1-delx).*cos(pi.*(dx-dx0)./(max(dx1)-min(dx1)+eps)).^2).*(dely+(1-dely).*cos(pi.*(dy-dy0)./(max(dy1)-min(dy1)+eps)).^2));%косинус квадрат на пьедестале
ar1=reshape(ar1',1,[]); 
ar2=reshape(ar2',1,[]);
% Результаты для сохранения документа txt.
x=x.*A;
y=y.*A;
ar1=ar1.*A;
Px=Px.*A;
K1=K1.*A;
z=z.*A;
T=[x;y;z;ar1;Px];
T=T';
plot(T(:,1),y,'dr');
axis equal
Px=Px*(-1);
h = [{'baseline'} {'top'}];
for ik=1:N
str=['(',num2str((floor((ar1(ik))*100))/100),',',num2str(Px(ik)*180/pi),')'];
text(T(ik,1),y(ik),str,'VerticalAlignment',h{mod(ik,2)+1},'FontSize', 7)
end