clc; close all; clear all;

f1=10;    %frecuencia de rotacion del disco
var = 0; %porcentaje de variacion de velocidad
N =1;    %numero de rayos del disco
Fs=30;   %frecuencia de sampleo
V= 20;   %numero de rotaciones del disco

Ts=1/Fs;
%t=[0:Ts:V/f1];  % numero de vueltas
t=[0:Ts:.5];     % tiempo finito
f_var=f1*var/100*(2*t/t(end)-1)+f1;

h1=figure(1);
k = (0:0.01:1)'*2*pi;
x = sin(k);
y = cos(k);

phi0=[0:1:N-1]'/N*2*pi;
x0=zeros(N*2,1);
y0=zeros(N*2,1);

M=moviein(length(t));
for i=1:length(t)
   clf
   fill(x,y,'w','LineStyle','none')
   x0(2:2:N*2)= sin(2*pi*f_var(i)*t(i) + phi0);
   y0(2:2:N*2)= cos(2*pi*f_var(i)*t(i) + phi0);
   hold on
   plot(x0,y0,'LineWidth',8)
   hold off
   axis square

   texto=sprintf('f_r_o_t = %d rps  F_s = %d Hz',round(f_var(i)),Fs);
   title(texto,'FontSize',14)
   axis off
   M(:,i)=getframe(h1);
end
%movie2avi(M,'rd_n1_f10.avi','FPS',Fs/N/3,'quality',60, 'compression','indeo3')