Snt=100
frekv=50
omega=2*pi*frekv
dP0p=0.1
dPkp=1
Ukp=10
I0p=1
Un1=110
Un2=22

>> Znt=Un1^2/Snt
Znt =  121
>> Rk=(dPkp/100)*(Un1^2/Snt)
Rk =  1.2100
>> Rk1=Rk/2
Rk1 =  0.60500
>> Zk=(Ukp/100)*(Un1^2/Snt)
Zk =  12.100
>> Xs=sqrt(Zk^2-Rk^2)
Xs =  12.039
>> Ls=Xs/omega
Ls =  0.038322
>> Ls1=Ls/2
Ls1 =  0.019161
>> Gfe=(dP0p/100)*(Snt/Un1^2)
Gfe =  0.0000082645
>> Rfe=1/Gfe
Rfe =  121000
>> Y0=(I0p/100)*(Snt/Un1^2)
Y0 =  0.000082645
>> Xh=1/sqrt(Y0^2-Gfe^2)
Xh =  12160.95756
>> Lh=Xh/omega
Lh =  38.709

>> Tau0=Lh/Rk1
Tau0 =  63.983
>> Z0=sqrt(Rk1^2+Xh^2)
Z0 =  12160.95758
>> I0m=Un1/sqrt(3)*sqrt(2)/Z0
I0m =  0.0073855
>> Psi0=atan(Xh/Rk1)
Psi0 =  1.5707
>> C0=I0m*sin(Psi0)
C0 =  0.0073855

>> t=[0:0.001:0.1]';
>> I0t=C0*exp(-t/Tau0) + I0m*sin(omega*t-Psi0);
>> plot(t,I0t);

>> In=Snt/Un1/sqrt(3)
In =  0.52486

>> TauK=Ls/Rk
TauK =  0.031671
>> Ikm=Un1*Ukp/100/sqrt(3)*sqrt(2)/Zk
Ikm =  0.74227
>> PsiK=atan(Xs/Rk)
PsiK =  1.4706
>> Ck=Ikm*sin(PsiK)
Ck =  0.73855

>> Ikt=Ck*exp(-t/TauK) + Ikm*sin(omega*t-PsiK);
>> plot(t,Ikt);
