% ultrasonido
clear
N=20;%arreglo cuadrado de NxN antenas
lam=34000/40000;
d=.1*lam, %distancia entre antenas
a=0:.03:pi;
ng=length(a);
p=2*pi*d*cos(a);
P=p'*ones(1,ng);
g=abs(sin(N*P/2)./sin(P/2)/N).*abs(sin(N*P'/2)./sin(P'/2)/N);
A=a'*ones(1,ng);B=A';
K=0;%ajuste de duagrama de radiacion de cada antena
go=(abs(sin(A)).*abs(sin(B))).^K;
%go=1;%surf(180/pi*a,180/pi*a',go)
surf(180/pi*a,180/pi*a',g.*go),
for i=1:180 view(i-1,i-1),pause(.01),end
surfl(180/pi*a,180/pi*a',g.*go),
shading interp,colormap(gray);
for i=1:180 view(i-1,i-1),pause(.01),end