imic3.m

关于电机的.m程序

%                <<<   INDUCTION MOTOR INITIALIZATION >>>
%  Motor is initialized at a specified operating condition established by speed
%clc,close all,clear
disp('                                                            imic3.m')
disp('        <<<    Induction motor transients  >>>         ')
disp('               [ Sinusoidal excitation  ]              ')
disp('               ( per-unit formulation  )               ')
disp('          Parameters of induction machine in per unit    ')
%     DATA
R_s=findobj('Tag','R_s');
Rs=eval(get(R_s,'String'));
R_r=findobj('Tag','R_r');
Rr=eval(get(R_r,'String'));
L_s=findobj('Tag','L_s');
Lsl=eval(get(L_s,'String'));
L_r=findobj('Tag','L_r');
Lrl=eval(get(L_r,'String'));
L_m=findobj('Tag','L_m');
Lm=eval(get(L_m,'String'));
H_i=findobj('Tag','H_i');
Hi=eval(get(H_i,'String'));
v_s=findobj('Tag','v_s');
vs=eval(get(v_s,'String'));
w_s=findobj('Tag','w_s');
ws=eval(get(w_s,'String'));
fo_s=findobj('Tag','fo_s');
freq0=eval(get(fo_s,'String'));
n0_i=findobj('Tag','n0_i');
n0=eval(get(n0_i,'String'));
%%%%%%%%%%%%%%%%%
disp('                                                         ')
disp('    Rs        Lsl       Rr        Lrl       Lm       vs         ws')
param=[Rs Lsl Rr Lrl Lm vs ws];
disp(param)
wo=2*pi*freq0;
Ls=Lsl+Lm ; Lr=Lrl+Lm;
j=sqrt(-1);
    Zr=Rr/(ws-n0)+j*Lrl; Zm=j*Lm*Zr/(j*Lm+Zr);  Zin=Rs/ws+j*Lsl+Zm;
    Is=vs/(ws*Zin);  Ir=-Is*j*Lm/(j*Lm+Zr);
    Req=real(Zm);    Tl0=Is*conj(Is)*Req;
    phis=Ls*Is+Lm*Ir;  phir=Lm*Is+Lr*Ir;    phids=real(phis);  phiqs=imag(phis);
    phidr=real(phir);  phiqr=imag(phir);
    op=[n0  Tl0];    disp('       n      Tl  ')  
    disp(op) 
        L=[Lsl+Lm   0      Lm     0   
          0   Lsl+Lm    0     Lm
          Lm    0     Lrl+Lm  0
          0     Lm      0     Lrl+Lm]
ic=[phids phiqs phidr phiqr n0 0]'              %initial conditions 
% ic=zeros(6,1)