%Calcula Distorsion Armonica Total (THD)
wb = waitbar(0,'Calculando THD...');
perc=passo;
npup=ms*p-1;

TF=2*abs(fft(vo(np-npup:np)/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/8),0.5)