% ELE-8401 Machines et entraînements électriques. Hiver 2018. Jesus Morales.
close all; clear; clc

dt = 1e-4;   tmax = 5/60;   t = 0:dt:tmax;   Ns = length(t);

f = 60;             w = 2*pi*f;

vabc = 10*[ cos(w*t); cos(w*t-2*pi/3); cos(w*t+2*pi/3) ]; 

figure; plot(t,vabc); legend('v_a','v_b','v_c')

%____________________________________________________________________
% Transformation au cadre qd0

vqd0 = zeros(3,Ns);   
alpha = 0;
for n = 1:Ns
    
    th = w*t(n) + alpha;
    K = 2/3*[ cos(th)   cos(th-2*pi/3)   cos(th+2*pi/3)
              sin(th)   sin(th-2*pi/3)   sin(th+2*pi/3)
                1/2          1/2              1/2       ];
 
    vqd0(:,n) = K * vabc(:,n);
    
end

figure; plot(t,vqd0); legend('v_q','v_d','v_0')