файлы матлаб / 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));