%Calcula Distorsion Armonica Total (THD)
wb = waitbar(0,'Calculando THD...');
perc=passo;
ntpuntos=length(Vo);    % numero total de puntos  
npup=(ntpuntos/nc);     % numero de puntos de 1 periodo

TF=2*abs(fft(Vo(np-(np/nc):np-1)/Vbase))/(npup);
THC(2)=TF(2)^2;
THC(3)=TF(3)^2;
fun=TF(2);

for q=4:(npup/2)
   THC(q)=TF(q)^2+THC(q-1);
   if q==(npup*perc/2)     % si alcanzo el porcentaje de actualizacion...
      waitbar(perc)         % actualiza la barra de progresion
      perc=perc+passo;      % próximo porcentaje de actualizacion
   end
end
close(wb);

RHC=sqrt(THC(q));
THD=100*RHC/fun
figure
bar(TF(2:ms/2),0.5)