%Representación visual del aliasing.
%Efecto de ir aumentando la velocidad de giro hasta q se sobrepasa la frecuencia de Nyquist. En este momento se
%produce un cambio de velocidad, pasando a girar a velocidad contraria y aumentando paulatinamente su velocidad.

%creacción del vector giratorio
% v = exp(-j*w*t);
%
%vector de tiempos del sistema muestreado
t = linspace(0,1,101); % <-- fm = 1/100;
%vector de tiempos de la señal real ( para simular la señal real, instroducimos una resolución muy alta)
tr = linspace(0,1,1001); % <-- fm = 1/1000;
%frecuencia real
fr = 1/(tr(2)-tr(1));
%frecuencia de muestreo (100 Hz)
fm = 1/(t(2)-t(1));
%frecuencia de Nyquist (Hz)
FN = fm/2;
%velocidad angular (Hz) <---- PARÁMETRO PARA COMPROBAR EL ALIASING
w = 50;
%vector giratorio muestreado
v = exp(-j*2*pi*w*t);
%vector giratorio real
vr = exp(-j*2*pi*w*tr);

%dibujamos el vector de manera animada
for cc = 1:length(t)
   figure(1);
   pause((t(2)-t(1)));
   clf;
   polar([0 angle(v(cc))],[0 abs(v(cc))],'ro-');hold on;
   polar([0 angle(vr(cc))],[0 abs(vr(cc))],'o-');hold off;
   title(strcat('Velocidad de giro: ',sprintf('%d Hz.\n',w),'Azul: Señal en cámara de alta velocidad. Roja: Señal en cámara normal'));
end

%POSIBILIDADES DE AMPLIACIÓN:
%
%	- Hacer un bucle for, en el cual se vaya cambiando la frecuencia de giro de manera automática, observando el proceso
%	  desde un correcto muestreo hasta que se produce aliasing.
%
%	- Dibujar simultaneamente la señal girando en el sentido correcto
%
%	- Comprobar que ocurre el mismo efecto tanto si cambiamos la frecuencia de muestreo como si se cambia la velocidad
%	  de giro.
%
