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

clear, clc
% параметры синхронного генератора
Tr=2;
md=0.87;
xq=0.83;
xd=2;
xds=0.26;
%параметры нагрузки
Xn1=0.6;
%Параметры регулятора напряжения
kk=25;
Tk=0.25; %u0
Un=1;
ku=1.25;
ki=1.25;

%Количество шагов
N=6000;
dt=0.001;
t=0;
%начальные условия
iy=0.3;
Fr=1.005;
for k=0:N
t=t+dt;
t0=1;t01=3.5;
if (t<t0 || t>t01)
Xn=10000;
else
Xn=Xn1;
end
id=Fr/(Xn+xds);
Uq=Fr-xds*id;
ur=ku*Uq+ki*xd*id-iy;
if ur>4
ur=4;
end
if ur<0
ur=0;
end
ir=Uq+xd*id;
pFr=(ur-ir)/Tr;Fr=Fr+pFr*dt;
piy=(-iy+kk*(Uq-Un))/Tk;iy=iy+piy*dt;
if iy<0
iy=0;
end
if iy>1
iy=1;
end

T(k+1)=t;
T1(k+1)=t;
T2(k+1)=t;
T3(k+1)=t;
T4(k+1)=t;
UQ(k+1)=Uq;
Ir(k+1)=ir;
ID(k+1)=id;
fr(k+1)=Fr;
UR(k+1)=ur;
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');
%вывод на печать значений Uq fr id Ur при t=3.51c
disp(UQ(3501));disp(fr(3501));disp(ID(3501));disp(UR(3501));