%control 02 pregunta 03
clearvars;
close all;
clc;


%arreglo de theta (angulo)
theta = 0:2*pi/3600:2*pi;

%variables importantes del grafico polar
N = 9;
f = 5000;
c = 344;
k = 2*pi*f/c;
d = 0.05;

%calculamos la parte de la función que no depende de delta
sinThetaDelta = sin(theta);
kddiv2 = k*d/2;
numerador = sin(N*kddiv2*sinThetaDelta);
denominador = sin(kddiv2*sinThetaDelta);

%graficamos
for delta=0:2*pi/360:2*pi
  rho = abs( sin(N*k*d/2*sin(theta-delta) )./(sin(k*d/2*sin(theta-delta) ))  )/N;

  figure(1)
  polar(theta,rho)
  title('patron polar arreglo lineal');
  grid on;
  box on;

  pause(0.01);
end