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

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

%параметры нагрузки
Xn1=1;
%Параметры регулятора напряжения
kk=50;
Tk=0.1;
Un=1;
ku=1.3;
ki=1.3;
%Количество шагов
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');
%вывод на печать значений w mv i при t=3.51c
disp(UQ(3001));disp(fr(3001));disp(ID(3001));disp(UR(3001));