imrunsif.m

关于电机的.m程序

% Induction machine model
function df=imrunsif(t,f)
global Rs Rr L J Vr fs fk wo km Tl p 
% ------------------ Step input excitations -------------
      Vs=Vr  ;
    if t<0.2   Vs=Vr ;
       elseif t>=0.2 & t<=0.6 Vs=Vr;Tl=11.90;
       elseif t>0.6 & t<=1.0  Tl=0.12;
      % else    Tl=5.95   ; end 
       else    fs=5/6 ; Vs=5/6*Vr ;Tl=0.12;
    end
% ------------------------------------------------------    
     vas=Vs*cos(f(7));    vbs=Vs*sin(f(7));
     vds=vas*cos(f(6))+vbs*sin(f(6))   ; vqs=-vas*sin(f(6))+vbs*cos(f(6));
     flux=[f(1) f(2) f(3) f(4)]';
% --------------------- Dynamic equations of induction machine ----
      if km==2, fk=f(5)/wo; end
           cur=L\flux;
       df =[vds-Rs*cur(1)+fk*wo*f(2)
            vqs-fk*wo*f(1)-Rs*cur(2)
           -Rr*cur(3)+f(4)*(wo*fk-f(5))
           -Rr*cur(4)-f(3)*(wo*fk-f(5))
            0.5*p*((cur(2)*f(1)-cur(1)*f(2))*0.75*p-Tl)/J
            wo*fk
            wo*fs];
% ------------------------------------------------------------------