dfig.m

这是一个用S函数写成的在altha/beta坐标系下的双馈风力发电机模型。

function [sys,x0,str,ts] =DFIG(t,x,u,flag,p,Rs,Ls,Rr,Lr,Lm,J,np)
Lm=0.06931;Lr=0.002;Ls=0.002;Rs=0.435;Rr=0.816;p=2;J=0.001;np=314;
switch flag
  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
  case 1,
    sys=mdlDerivatives(t,x,u,Rs,Ls,Rr,Lr,Lm,np,J);%p代表极对数
  case 2,
    sys=[];
  case 3,
    sys=mdlOutputs(t,x,u,flag,Lm);
  case {4,9}
    sys=[];
otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end

function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 5;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 5;
sizes.NumInputs      = 5;
sizes.DirFeedthrough = 1;%输出方程中用到了输入量
sizes.NumSampleTimes = 1;  
sys = simsizes(sizes);
x0  = [0 0 0 0 0];
str = [];
ts  = [-1 0];
function sys=mdlDerivatives(t,x,u,Rs,Ls,Rr,Lr,Lm,np,J)
l1=[Ls 0 Lm 0;0 Ls 0 Lm;Lm 0 Lr 0;0 Lm 0 Lr];
D=inv(l1);
R=[-Rs 0 0 0;0 -Rs 0 0;0 -x(5)*Lm -Rr -x(5)*Lr;x(5)*Lm 0 x(5)*Lr -Rr];
C=D*R;
sys(1:4)=C*x(1:4)+D*u(1:4);
sys(5)=np/J*(1.5*np*Lm*(x(2)*x(3)-x(1)*x(4))-u(5))
function sys=mdlOutputs(t,x,u,Lm,p)
sys=x;

%在这里有五个状态变量[isd,isq,ird,irq，wr],五个输入变量[Usd,Usq,Urd,Urq,T],
%五个输出变量[isd,isq,ird,irq,wr]
%