% Definimos el paso y el vector que recorra el plano de coordenadas.

delta = (2*pi)/300;
phi = 0:delta:(2*pi);
n = [1, 2, 3, 4, 5, 10, 25, 50:343];
n(50<n) = 0;

% Definimos la función a representar.

i = 1;

if i == 1
    funcion = (sin(2.5*pi*sin(phi))./(2.5*pi*sin(phi))).^2;
elseif i == 2
    funcion = (sin(7*pi*sin(phi))./(7*pi*sin(phi))).^2;
else
    funcion = (cos(phi)).^n;
end

% Dividimos la función entre su valor máximo para normalizar entre 0 y 1.

funcion=funcion./(max(funcion));

% Expresamos la función en escala logarítmica entre 0 y 30 dB.

funcion = 10*log10(funcion);
funcion(funcion<-30) = -30;
funcion = funcion+30;

% Y representamos la función con los ejes normalizados entre -30 y 0 dB.

pantalla = get(0, 'ScreenSize');
ventana = [pantalla(3)/2-300, pantalla(4)/2-300, 600, 600];
figure('Name', 'Diagrama de radiación en 2D', 'Position', ventana);

polar(phi, funcion);