% Prechodne deje realizovane na transformatoru a vedeni

% Vypocet chodu nakratko
TauK=Ls/Rk
Ikm=Un1*Ukp/100/sqrt(3)*sqrt(2)/Zk
PsiK=atan(Xs/Rk)
% Rozdil faze oproti ryze induktivnímu proudu
OdchylkaPsiK=(PsiK-pi/2)/(2*pi)*0.02
% Integracni konstanta
Ck=Ikm*sin(PsiK)

t=0:0.0001:0.1;
Ikt=Ck*exp(-t/TauK) + Ikm*sin(omega*t-PsiK);
plot(t,Ikt);
grid on;

% Kontrola pres komplexni vyjadreni Ik
Ikc=Un1*Ukp/100/sqrt(3)/(Rk+i*Xs)
% Amplituda alternativne
Ikm2=sqrt(2)*abs(Ikc)
% Faze alternativne 
PsiK2=angle(Ikc)

% Kontrola pres jmenovity proud
Ikm3=In*sqrt(2)
I0m3=I0p/100*In*sqrt(2)
