% PRACTICA 3: ANALISIS FRECUENCIAL COMPLEJO DEL VECTOR DE PARK
% SCRIPT:     filtro_fir_1.m
% DESCRIPCION: 
%     Para distintos ensayos sobre un motor asincrono se analizan las corrientes 
%     ir(t), is(t), it(t)=-ir(t)-is(t) en el motor utilizando el enfoque del vector 
%     de Park
%
%                       i = ir(t) + a*is(t) + a^2*it(t)
%
%     Se describe aqui el analisis frecuencial de esta secuencia compleja.
%
% AUTOR
%       Ignacio Diaz Blanco. 
%       Area de Ingenieria de Sistemas y Automatica.
%       Universidad de Oviedo. 
%       Enero 2002

close all;
load dataicann;
z   = z  ([1 2 3 7 4 5 6 8 9]);
nom = nom([1 2 3 7 4 5 6 8 9]);


fm = 5000;


for Registro=1:9                    % Para cada ensayo
    N  = length(z{Registro});       % Determinar la longitud del registro
    
    % Determinar las corrientes ir(t), is(t), it(t)
    ir = z{Registro}(:,4); 
    is = z{Registro}(:,5); 
    it = -ir-is;
    
    % Eliminar sus valores medios (debidos a "bias" en acelerometros, etc.)
    ir = ir - mean(ir);
    is = is - mean(is);
    it = it - mean(it);

    % Generar el vector de Park
    a  = exp(j*2*pi/3);
    ipark = ir + a*is + a^2*it;    
    figure(1);
    subplot(3,3,Registro);
    plot(ipark);
    axis equal;
    axis([-2 2 -2 2]);
    title(nom(Registro));
    
    % Analisis frecuencial    
    f = linspace(-fm/2,fm/2,N+1); 
    f = f(1:end-1);
    figure(2);
    subplot(3,3,Registro);
    plot(f,fftshift(abs(fft(ipark))));
    title(nom(Registro));
    
end;
