% Prechodne deje relizovane na transformatoru
% a vedeni

% Chod nakratko
Ikm=(Uk/100)*Um/Zk
TauK=Ls/Rk
PsiK=atan(Xs/Rk)

t=[0:0.001:0.1]';
Ikt=Ikm*(sin(PsiK)*exp(-t/TauK)+sin(Omega*t-PsiK));
plot(t,Ikt);
grid on;
