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

% Generator
Png=700
CosFi=0.95
Sng=Png/CosFi
Xd= (200/100)*(Un^2/Sng)
Xdc=( 80/100)*(Un^2/Sng)

% Transformator T1
Snt1=700
Ukt1=8
Xt1=(Ukt1/100)*(Un^2/Snt1)

% Vedeni
delka=50
Xv1=delka*0.4
Xv=Xv1/2

% Transformator T2
Snt2=700
Ukt2=5
Xt2=(Ukt2/100)*(Un^2/Snt2)

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

% Vypocet predporuchoveho stavu, pocatecni zatizeni 100%
Pt=Png
Ig=Png/CosFi/(3*Unf)
Igcinny =Ig*CosFi
Igjalovy=Ig*sin(acos(CosFi))
Is=(Png-i*Png*tan(acos(CosFi)))/sqrt(3)/Un

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

% Vypocet potrebneho nabuzeni
% Ustalane vnitrni el.-mot. napeti E dle fazoroveho diagramu
Ekomplexni=Unf + i*Xsuma*(Igcinny-i*Igjalovy)
Ekomplexni=Unf + i*Xsuma*Is
E=abs(Ekomplexni)
Theta=angle(Ekomplexni)

% Prechodne vnitrni el.-mot. napeti E'
Eckomplexni=Unf + i*Xsumac*(Igcinny-i*Igjalovy)
Eckomplexni=Unf + i*Xsumac*Is
Ec=abs(Eckomplexni)
Thetac=angle(Eckomplexni)

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=Xl + Xr + Xl*Xr/Xc

% Odpojeni poruchy
Xsumaodpojeno=Xsumac -Xv + Xv1
Xsumaodpojeno2=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
