ЛР3 / scriptAD

%===================parameters=======================
P2n=55*10^3; %Номинальная мощность
U1n=220; %Номинальное напряжение
w0=2*pi*50; %rated frequency of rotation of the magnetic field
Zp=2; %number of pole pairs
J=2; %момент инерции ротора с барабаном
KPD=0.75; %Efficiency
cosf=0.9; %Power factor
Xm0=4.2; %Inductive Magnetization Resistance
R1g=0.027; %stator winding resistance
X1g=0.086; %stator winding inductance
R2g=0.015; %rotor winding resistance
X2g=0.14; %rotor winding inductance
%===========calculation of the parameters of the T-shaped model==========
%stator rated phase current
In1=(P2n)/(3*U1n*KPD*cosf);
%inductance T-shaped equivalent circuit
X1t=(2*Xm0*X1g)/(Xm0+sqrt((Xm0^2)+4*Xm0*X1g));
%conversion factor
c=X1g/X1t;
%stator winding resistance (T-shaped)
R1t=R1g/c;
%rotor winding resistance (T-shaped)
R2t=R2g/(c^2);
%rotor winding inductance (T-shaped)
X2t=X2g/(c^2);
%Base resistance
Zbaz=U1n/In1;
%===============conversion to named units================================
%stator winding resistance
R1=R1t*Zbaz;
%rotor winding resistance
R2=R2t*Zbaz;
%Inductive Magnetization Resistance
Xm=Xm0*Zbaz;
%Stator Scatter Resistance
X1s=X1t*Zbaz;
%Rotor Scatter Resistance
X2s=X2t*Zbaz;
%Stator impedance
X1=X1s+Xm;
%Rotor impedance
X2=X2s+Xm;
%Magnetization loop inductance
Lm=Xm/w0;
%Stator Inductance
L1s=X1s/w0;
%Rotor Inductance
L2s=X2s/w0;
%Total stator inductance
L1=Lm+L1s;
%Total rotor inductance
L2=Lm+L2s;
%scattering coefficient
sig=1-((Lm^2)/(L1*L2));
%Electromagnetic Rotor Time Constant
T2=L2/R2;
%The coefficient of the mathematical model
L2st=L1-((Lm^2)/L2);
k2=Lm/L2;
r=R1+(k2^2)*R2;
T1=L2st/r;
Tg=1; %Постоянная времени жидкости в барабане
%Коэффициенты полинома
A=0;
B=0.5;
C=1;