%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 por separado
sinThetaDelta = sin(theta);
kddiv2 = k*d/2;
numerador = sin(N*kddiv2*sinThetaDelta);
denominador = sin(kddiv2*sinThetaDelta);

rho = numerador./denominador/N;

%también podemos hacerlo todo en la misma expresión
%rho = abs(sin((N*k*d/2)*sin(theta)) ./ sin((k*d/2)*sin(theta)) / N);

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