% Napeti
i=sqrt(-1)
Un=400
Unf=Un/sqrt(3)

% Generator
Sng=100
CosFi=0.8
Png=Sng*CosFi
Xd= (120/100)*(Un^2/Sng)*i
Xdc=( 30/100)*(Un^2/Sng)*i

% Transformator T1
Snt1=100
Ukt1=10
Xt1=(Ukt1/100)*(Un^2/Snt1)*i

% Vedeni
delka=50
Xv1=delka*0.3*i/110^2*Un^2
Xv=Xv1/2

% Transformator T2
Snt2=100
Ukt2=10
Xt2=(Ukt2/100)*(Un^2/Snt2)*i

% Sit
Sk=10000
KoefC=1
Xs=KoefC*(Un^2/Sk)*i

% Vypocet predporuchoveho stavu, pocatecni zatizeni 75%
% Zvolit bud' velikost E, nebo CosFiS aktualni, 
% Treti moznost je navolit zatezny uhel
% Zde volen zatezny uhel Theta
Png0=0.75*Png
Pt=Png0
Theta=60/180*pi
Thetac=Theta

% Celkova vazebni reaktance
Xsuma =abs(Xd +Xt1+Xv+Xt2+Xs)
Xsumac=abs(Xdc+Xt1+Xv+Xt2+Xs)

% Vypocet potrebneho nabuzeni E a E'
E= Png0/3/Unf*Xsuma/sin(Theta)
Ec=Png0/3/Unf*Xsumac/sin(Theta)

Pgkontrola=3*E*Unf/Xsuma*sin(Theta)
Pgckontrola=3*Ec*Unf/Xsumac*sin(Thetac)

% Zkrat v polovine vedeni Xv1

% Nahrada trojuhelnika hvezdou
Xa=Xv1*Xv1/2/(Xv1+Xv1/2+Xv1/2)
Xb=Xv1*Xv1/2/(Xv1+Xv1/2+Xv1/2)
Xc=Xv1/2*Xv1/2/(Xv1+Xv1/2+Xv1/2)

Xl=Xdc+Xt1+Xa
Xr=Xb+Xt2+Xs
Xsumazkrat=abs(Xl + Xr + Xl*Xr/Xc)

% Odpojeni poruchy
Xsumaodpojeno=abs(i*Xsumac -Xv + Xv1)
Xsumaodpojeno2=abs(Xdc + Xt1 + Xv1 + Xt2 + Xs)

% Konecne alternativy maximalnich dodavanych cinnych vykonu
Pmax1=3*Ec*Unf/Xsumac
Pmax2=3*Ec*Unf/Xsumazkrat
Pmax3=3*Ec*Unf/Xsumaodpojeno

% Kontrolni admitance
Y12=        1/Xsumac
Y12zkrat=   1/Xsumazkrat
Y12odpojeno=1/Xsumaodpojeno
