vdip.m

来自「This file is hysterisis two level contro」· M 代码 · 共 63 行

function [y]=vdip(vector)

% Desbalanceo por disminucion de amplitud en la fase b
t=vector(1); % tiempo

% Datos de las senales
Vllrmsa1=400;
Vllrmsb1=400;
Vllrmsc1=400;
f1=50;
phasea1=0;
phaseb1=-120*pi/180;
phasec1=120*pi/180;
Vllrmsa2=400*0.250;
Vllrmsb2=400*0.25;
Vllrmsc2=400*0.25;
f2=50;
phasea2=0;
phaseb2=-120*pi/180;
phasec2=120*pi/180;



if (t <= 0.3) 
    y(1)=Vllrmsa1*sqrt(2/3)*cos(2*pi*f1*t+phasea1);
    y(2)=Vllrmsb1*sqrt(2/3)*cos(2*pi*f1*t+phaseb1);
    y(3)=Vllrmsc1*sqrt(2/3)*cos(2*pi*f1*t+phasec1); 
elseif (t > 0.3) && (t <= 0.31)  % Cambio en medio ciclo
    y(1)=(((Vllrmsa2-Vllrmsa1)*(t-0.3)/0.01)+Vllrmsa1)*sqrt(2/3)*cos(2*pi*f2*t+phasea2);
    y(2)=(((Vllrmsb2-Vllrmsb1)*(t-0.3)/0.01)+Vllrmsb1)*sqrt(2/3)*cos(2*pi*f2*t+phaseb2);
    y(3)=(((Vllrmsc2-Vllrmsc1)*(t-0.3)/0.01)+Vllrmsc1)*sqrt(2/3)*cos(2*pi*f2*t+phasec2); 
elseif (t > 0.31) && (t <= 0.4) 
    y(1)=Vllrmsa2*sqrt(2/3)*cos(2*pi*f1*t+phasea2);
    y(2)=Vllrmsb2*sqrt(2/3)*cos(2*pi*f1*t+phaseb2);
    y(3)=Vllrmsc2*sqrt(2/3)*cos(2*pi*f1*t+phasec2); 
elseif (t > 0.4) && (t <= 0.41)  % Cambio en medio ciclo
    y(1)=(((Vllrmsa1-Vllrmsa2)*(t-0.4)/0.01)+Vllrmsa2)*sqrt(2/3)*cos(2*pi*f2*t+phasea2);
    y(2)=(((Vllrmsb1-Vllrmsb2)*(t-0.4)/0.01)+Vllrmsb2)*sqrt(2/3)*cos(2*pi*f2*t+phaseb2);
    y(3)=(((Vllrmsc1-Vllrmsc2)*(t-0.4)/0.01)+Vllrmsc2)*sqrt(2/3)*cos(2*pi*f2*t+phasec2); 
elseif (t > 0.41) && (t <= 0.75) 
    y(1)=Vllrmsa1*sqrt(2/3)*cos(2*pi*f1*t+phasea1);
    y(2)=Vllrmsb1*sqrt(2/3)*cos(2*pi*f1*t+phaseb1);
    y(3)=Vllrmsc1*sqrt(2/3)*cos(2*pi*f1*t+phasec1); 
elseif (t > 0.75) && (t <= 0.76)  % Cambio en medio ciclo
    y(1)=(((Vllrmsa2-Vllrmsa1)*(t-0.75)/0.01)+Vllrmsa1)*sqrt(2/3)*cos(2*pi*f2*t+phasea2);
    y(2)=(((Vllrmsb2-Vllrmsb1)*(t-0.75)/0.01)+Vllrmsb1)*sqrt(2/3)*cos(2*pi*f2*t+phaseb2);
    y(3)=(((Vllrmsc2-Vllrmsc1)*(t-0.75)/0.01)+Vllrmsc1)*sqrt(2/3)*cos(2*pi*f2*t+phasec2); 
elseif (t > 0.76) && (t <= 0.85) 
    y(1)=Vllrmsa2*sqrt(2/3)*cos(2*pi*f1*t+phasea2);
    y(2)=Vllrmsb2*sqrt(2/3)*cos(2*pi*f1*t+phaseb2);
    y(3)=Vllrmsc2*sqrt(2/3)*cos(2*pi*f1*t+phasec2); 
elseif (t > 0.85) && (t <= 0.86)  % Cambio en medio ciclo
    y(1)=(((Vllrmsa1-Vllrmsa2)*(t-0.85)/0.01)+Vllrmsa2)*sqrt(2/3)*cos(2*pi*f2*t+phasea2);
    y(2)=(((Vllrmsb1-Vllrmsb2)*(t-0.85)/0.01)+Vllrmsb2)*sqrt(2/3)*cos(2*pi*f2*t+phaseb2);
    y(3)=(((Vllrmsc1-Vllrmsc2)*(t-0.85)/0.01)+Vllrmsc2)*sqrt(2/3)*cos(2*pi*f2*t+phasec2); 
else
    y(1)=Vllrmsa1*sqrt(2/3)*cos(2*pi*f1*t+phasea1);
    y(2)=Vllrmsb1*sqrt(2/3)*cos(2*pi*f1*t+phaseb1);
    y(3)=Vllrmsc1*sqrt(2/3)*cos(2*pi*f1*t+phasec1); 
end