% Prechodne deje realizovane na transformatoru a vedeni

% Vypocet chodu naprazdno
Tau0=Lh/Rk1
Z0=sqrt(Rk1^2+Xh^2)
I0m=Un1/sqrt(3)*sqrt(2)/Z0
Psi0=atan(Xh/Rk1)
% Rozdil faze oproti ryze induktivnímu proudu
OdchylkaPsi0=(Psi0-pi/2)/(2*pi)*0.02
% Integracni konstanta
C0=I0m*sin(Psi0)

t=0:0.001:0.1;
I0t=C0*exp(-t/Tau0) + I0m*sin(omega*t-Psi0);
plot(t,I0t);
grid on;

% Kontrola pres komplexni vyjadreni I0
I0c=Un1/sqrt(3)/(Rk1+i*Xh)
% Amplituda alternativne
I0m2=sqrt(2)*abs(I0c)
% Faze alternativne 
Psi02=angle(I0c)
