senales1.m

来自「This file is used to develop Phase locke」· M 代码 · 共 48 行

function [y]=senales1(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;
Vllrmsb2=350;
Vllrmsc2=400;
f2=50;
phasea2=0;
phaseb2=-120*pi/180;
phasec2=120*pi/180;



if (t <= 0.2) 
    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.2) && (t <= 0.21)  % Cambio en medio ciclo
    y(1)=(((Vllrmsa2-Vllrmsa1)*(t-0.2)/0.01)+Vllrmsa1)*sqrt(2/3)*cos(2*pi*f2*t+phasea2);
    y(2)=(((Vllrmsb2-Vllrmsb1)*(t-0.2)/0.01)+Vllrmsb1)*sqrt(2/3)*cos(2*pi*f2*t+phaseb2);
    y(3)=(((Vllrmsc2-Vllrmsc1)*(t-0.2)/0.01)+Vllrmsc1)*sqrt(2/3)*cos(2*pi*f2*t+phasec2); 
elseif (t > 0.21) && (t <= 0.3) 
    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.3) && (t <= 0.31)  % Cambio en medio ciclo
    y(1)=(((Vllrmsa1-Vllrmsa2)*(t-0.3)/0.01)+Vllrmsa2)*sqrt(2/3)*cos(2*pi*f2*t+phasea2);
    y(2)=(((Vllrmsb1-Vllrmsb2)*(t-0.3)/0.01)+Vllrmsb2)*sqrt(2/3)*cos(2*pi*f2*t+phaseb2);
    y(3)=(((Vllrmsc1-Vllrmsc2)*(t-0.3)/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