imssmod3.m

来自「关于电机的.m程序」· M 代码 · 共 23 行

% Induction machine steady-state model used in imss3.mdl
function [f,x0]=imssmod3(t,x,u,flag,Rs,Rr,Lsl,Lrl,Lm,wm)
switch flag,
    case 3,
         %              Outputs
         j=sqrt(-1);
    Zr=Rr./((u(2)-wm)+eps)+j*Lrl;
    Zm=j*Lm.*Zr./(j*Lm+Zr);
    Zin=Rs/u(2)+j*Lsl+Zm;
    Req=real(Zm);
    I=u(1)./(u(2)*Zin);
    Ir=real(I);Ii=imag(I);
    Te=I.*conj(I).*Req;
       f=[wm Te Ir Ii] ;
    case 0,
       m=4*length(wm);
       f=[0,0,m,2,0,1];
       x0=[];
    case {1,2,4,9},
         f=[];
       otherwise
         error(['unhandled flag = ',num2str(flag)]);
  end