dfim.m

来自「一个用S函数写成的在dq坐标系中的双馈异步电机模型。」· M 代码 · 共 43 行

function [sys,x0,str,ts] =DFIM2(t,x,u,flag,p,Rs,Ls,Rr,Lr,Lm)
switch flag
  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
  case 1,
    sys=mdlDerivatives(t,x,u,Rs,Ls,Rr,Lr,Lm);%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  = 4;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 5;
sizes.NumInputs      = 6;
sizes.DirFeedthrough = 0;%输出方程中用到了输入量
sizes.NumSampleTimes = 1;  
sys = simsizes(sizes);
x0  = [0 0 0 0];
str = [];
ts  = [-1 0];
function sys=mdlDerivatives(t,x,u,Rs,Ls,Rr,Lr,Lm)
l1=[-Ls 0 Lm 0;0 -Ls 0 Lm;-Lm 0 Lr 0;0 -Lm 0 Lr];
D=inv(l1);
R=[Rs -u(5)*Ls 0 u(5)*Lm;u(5)*Ls Rs -u(5)*Lm 0;0 -u(6)*Lm -Rr u(6)*Lr;u(6)*Lm 0 -u(6)*Lr -Rr];
C=D*R;
sys=C*x(1:4)+D*u(1:4);

function sys=mdlOutputs(t,x,u,Lm,p)
x(5)=1.5*p*Lm*(x(1)*x(4)-x(2)*x(3));%求Te
sys=x;

%在这里有四个状态变量[isd,isq,ird,irq],六个输入变量[Usd,Usq,Urd,Urq,W1,Ws,],
%五个输出变量[isd,isq,ird,irq,Te]
%