% Definice transformace
i=sqrt(-1)
A=-0.5+sqrt(3)/2*i
F=[ 1   1   1
    A*A A   1
    A   A*A 1 ]
pause

% Generator 10.5 kV
Png=200
CosFi=0.9
Sng=Png/CosFi
Xdc=(30/100)*(10.5^2/Sng)*(110^2/10.5^2)*i
pause 

% Transformator T1 10.5/110 kV
Snt1=250
Ukt1=10
Xt1=(Ukt1/100)*(110^2/Snt1)*i
pause 

% Vedeni
delka=10
Xv=delka*0.3*i

% Transformator T2 110/400 kV
Snt2=300
Ukt2=10
Xt2=(Ukt2/100)*(110^2/Snt2)*i
pause 

% Sit 400 kV
Sk=1000
Xk=(400^2/Sk)*(110^2/400^2)*i
pause

% Nahrada trojuhelnika hvezdou
Z12=(Xv/2)*Xv/(Xv/2 +Xv/2+Xv)
Z22=(Xv/2)*(Xv/2)/(Xv/2 +Xv/2+Xv)
Z32=(Xv/2)*Xv/(Xv/2 +Xv/2+Xv)
Z2=Z22+(Xdc+Xt1+Z12)*(Xk+Xt2+Z32)/((Xdc+Xt1+Z12)+(Xk+Xt2+Z32))
pause

% Protoze veskere impedance vedeni pocitame jako 3-nasobek
% hodnot pro prevod trojuhelnik hvezda ze zpetne soustavy
Z10=3*Z12
Z20=3*Z22
Z30=3*Z32
Z0=Z20+(Xt1+Z10)*(Xk+Xt2+Z30)/((Xt1+Z10)+(Xk+Xt2+Z30))
pause

% Pro 1-f zemni zkrat je zapojeni nahradnich reaktanci seriove
Zekv=Z2+Z0
pause

% Nahradni reaktance sousledne slozky je stejna jako u zpetne
Z11=Z12
Z21=Z22
Z31=Z32
Z1=Z2
pause

Unf=110000/sqrt(3)
I1=1.1*Unf/(Z1+Zekv)
pause

% Proudy ve zpetne a netocive slozkove soustave jsou stejne
I2=I1
I0=I1
pause

% Podil proudu zleva a zprava v sousledne soustave
Ia1=I1*(Xk+Xt2+Z31)/((Xdc+Xt1+Z11)+(Xk+Xt2+Z31))
Ib1=I1*(Xdc+Xt1+Z11)/((Xdc+Xt1+Z11)+(Xk+Xt2+Z31))
pause

% Podil proudu zleva a zprava v zpetne soustave
Ia2=I2*(Xk+Xt2+Z32)/((Xdc+Xt1+Z12)+(Xk+Xt2+Z32))
Ib2=I2*(Xdc+Xt1+Z12)/((Xdc+Xt1+Z12)+(Xk+Xt2+Z32))
pause

% Podil proudu zleva a zprava v netocive soustave
Ia0=I0*(Xk+Xt2+Z30)/((Xt1+Z10)+(Xk+Xt2+Z30))
Ib0=I0*(Xt1+Z10)/((Xt1+Z10)+(Xk+Xt2+Z30))
pause

% Napeti v miste poruchy
U1=I1*Zekv
U2=-I2*Z2
U0=-I0*Z0
pause

% Napeti na svorkach generatoru
Ug1=1.1*Unf-Xdc*Ia1
Ug2=-Xdc*Ia2
Ug0=0
pause

% Napeti na svorkach transformatoru
Ut11=1.1*Unf-(Xt1+Xdc)*Ia1
Ut12=-(Xt1+Xdc)*Ia2
Ut10=-Xt1*Ia0
pause

% Transformace zpetne napeti v miste poruchy
Upor=[U1 U2 U0]
UporN=F*Upor.'
pause

% Transformace zpetne napeti v miste svorek generatoru
Ug=[Ug1 Ug2 Ug0]
UgN=F*Ug.'
pause

% Transformace zpetne napeti v miste svorek transformatoru 1
Ut1=[Ut11 Ut12 Ut10]
Ut1N=F*Ut1.'
pause

