smss1.m

关于电机的.m程序

%                     Synchronous machine steady-state characteristics
clear,clc,close all
disp(date)
disp('           <<<  SYNCHRONOUS MACHINE CHARACTERISTICS  >>> ')
disp('     MACHINE PARAMETERS in per unit:                                  ')
%          Machine parameters (pu) :  
                    Ra=.005; Xd=1.15 ; Xq=0.75 ;
                    Vt=1.0; fs=1.0;
                    Ef=0.5;
disp('    Ra        Xd          Xq ')
disp([Ra Xd Xq])
j=sqrt(-1); den=Xd*Xq+Ra*Ra/(fs*fs); 
 Ef=[0 0.25 .5 .75 1 1.5 2];
figure('Units','normalized','Position',[.256 .05 .489 .32],...
    'Name','Steady-state characteristics','numbertitle','on','menubar','none');   
  % figure(1)
   bar(Ef),title('Ef as a parameter'),grid
for p=1:4
  for k=1:7
     del=-180:2:180;
     delta=del*pi/180;
     vd=-Vt*sin(delta); vq=Vt*cos(delta);
     id=(Xq*(Vt/fs*cos(delta)-Ef(k))-Ra*Vt*sin(delta)/(fs*fs))/den;
     iq=(Xd*Vt*sin(delta)/fs+Ra*(Vt*cos(delta)/fs-Ef(k)))/den;
     Te=Ef(k)*iq+(Xd-Xq)*id.*iq;
     P=id.*vd+iq.*vq ;
     Ia=id+j*iq; V=vd+j*vq; Imag=abs(Ia);
     theta=angle(Ia)-angle(V);
     Ireal=Imag.*cos(theta); Iimag=Imag.*sin(theta);
% ----------------------   PLOTS  ----------------------------------
   if p==1
    figure(2)
 %figure('Units','normalized','Position',[.45 .05 .49 .5],'Name','Steady-state characteristics','numbertitle','on');      
     axis([-180 180 -2 2])
    plot(del,Te,'m'),hold on
    xlabel('TORQUE ANGLE [o]'),ylabel('TORQUE [pu]')
    set(gca,'Xtick',-180:60:180),grid
   end
   if p==2
   figure(3)   
    axis([-180 180 -2 2]);
    plot(del,P,'c'),hold on
    xlabel('TORQUE ANGLE [o]'),ylabel('POWER [pu]')
    set(gca,'Xtick',-180:60:180),grid
   end
   if p==3
    figure(4)  
    axis([-180 180 0 3]);
    plot(del,Imag,'r'),hold on
    xlabel('TORQUE ANGLE [o]'),ylabel('CURRENT')
    set(gca,'Xtick',-180:60:180),grid
    end
   if p==4
      figure(5)
            plot(0,0,'ok',Ireal,Iimag,'r'),grid,hold on
            compass(Vt,'b'), 
            text(Vt-.2,-.1,'Vt'),title('CURRENT PHASOR LOCI')
    axis([-2 2 -3 1]),axis square,axis equal
    end
  end
 end