% EJEMPLO DE APLICACION DE LA FFT
%
%
% Descripcion:
% En esta practica se aplica la FFT para el analisis de la corriente de inducido
% de un motor de CC de 6000 Kw utilizado en laminacion. La corriente de
% excitacion presenta, ademas de su nivel de continua, un armonico a 600 Hz
% (rizado de 12 rizos por ciclo, correspondiente a señal trifasica rectificada 
% procedente de dos transformadores a 30 grados entre si) y otro armonico
% de 300 Hz que aparece cuando los dos puentes de tiristores asociados a
% ambos transformadores trabajan de forma desigual...
%
% AUTOR
%       Ignacio Diaz Blanco. 
%       Prof. Titular de Universidad
%       Area de Ingenieria de Sistemas y Automatica.
%       Universidad de Oviedo. 
%       Febrero 2004

load datosmotorlaminacion;

fm = 5000;                      % Frecuencia de muestreo
tm = inv(fm);                   % Periodo de muestreo
ttot = 0:tm:tm*(length(p)-1);   % Vector de tiempos de la señal a analizar

N = 4096;                       % Tamaño de la ventana de datos a analizar
f = linspace(0,fm,N+1);         % Vector de frecuencias asociado
f = f(1:N);                     % 


offset   = 10000;               % Posicion de la ventana de datos dentro del total
idx      = offset + (1:N);
variable = 1;                   % Variable a analizar (1=Ia, 2=Va, 3=If, 4=w)



figure(1);
for k = 1:4,
    subplot(4,1,k);
    plot(ttot,p(k,:));
    title(nombres(k));
    if k == variable
        hold on;
        plot(ttot(idx(1))*ones(1,20),linspace(min(p(k,:)),max(p(k,:)),20),'r');
        plot(ttot(idx(end))*ones(1,20),linspace(min(p(k,:)),max(p(k,:)),20),'r');
        hold off;
    end
end






y = p(variable,idx);

figure(2);
subplot(211);
plot(ttot(idx),y);
axis tight;
title(nombres(variable));
subplot(212);
plot(f,abs(fft(y))/N);

soundsc(y,fm);