файлы матлаб / lab5

clear,clc

% Параметры синхронного генератора
Tr = 2;
md = 0.87;
xq = 0.83;
xd = 2;
xds = 0.26;

% Параметры нагрузки
Xn1 = 0.6;
Rn = 0.8;

% Параметры регулятора напряжения
kk = 25;
Tk = 0.25;
Un = 1; % u0
ku = 1.25;
ki = 1.25;
Tj = 1.5;
Tp = 0.25;
kp = 33;
omega0 = 1;
Mxx = 0.2;

% Количество шагов
N = 6000;
dt = 0.001;
t = 0;

% Начальные условия
omega = 0;
mp = 1,3;
iy = 0.3;
Fr = 1.005;
dt = 0.001;
t = 0;

for k = 0:N,
t = t+dt; % шаг для значений
t0 = 1; % включение выключателя
t01 = 3.5; % выключение
if (t<t0 | t>t01)
Xn = 10000;
else
Xn = Xn1;
end

% Расчетная модель для исследования
id = (Xn+xq)*Fr/((Xn+xds)*(Xn+xq)+Rn^2);
iq = (Rn*id)/(Xn+xq);
i = sqrt(id^2+iq^2);
uq = Fr-xds*id; % Fr - это Psi_f
ud = xq*iq;
ur = ku*uq+ki*xd*id-iy;
if ur>4
ur = 4;
end
if ur<0
ur = 0;
end
ir = Fr+md*xd*id; % md - мю д
Me = uq*iq-ud*id; % Me - это Мэ
pFr = (ur-ir)/Tr; % уравнение p*Psi_f
Fr = Fr+pFr*dt;
piy = (-iy+kk*(uq-Un))/Tk;
iy = iy+piy*dt; % уравнение 17
if iy<0
iy = 0;
end
if iy>1
iy = 1;
end
pomega = (mp-Me-Mxx)/Tj;
omega = omega+pomega*dt;
pmp = (-mp-kp*(omega-omega0))/Tp;
mp = mp+pmp*dt;
if mp<0
mp = 0;
end
if mp>1.3
mp = 1.3;
end

T(k+1)=t;
T1(k+1)=t;
T2(k+1)=t;
T3(k+1)=t;
T4(k+1)=t;
T5(k+1)=t;

UQ(k+1)=uq;
Ir(k+1)=ir;
ID(k+1)=id;
fr(k+1)=Fr;
UR(k+1)=ur;
Omega(k+1)=omega;
end

plot(T,UQ,T1,1.05*fr,T2,0.5*Ir,T3,2*ID,T4,0.5*UR)
grid on;
xlabel('Время');
legend('UQ','fr','Ir','ID','UR');


%Вывод на печать значений при t=3.51c
disp(UQ(3001));
disp(fr(3001));
disp(ID(3001));
disp(UR(3001));
disp(Ir(3001));