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

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

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

% Vedeni
delka=100
Xv1=delka*0.4*i
Xv=Xv1/2

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

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

% Vypocet predporuchoveho stavu, pocatecni zatizeni 75%
% Zvolit bud' velikost E, nebo CosFiS aktualni
% Zde volen CosFiS
Png0=0.9*Png
Pt=Png0
Ig=Png0/CosFi/(3*Unf)
Igcinny =Ig*CosFi
Igjalovy=Ig*sin(acos(CosFi))

% 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)
E=abs(Ekomplexni)
Theta=angle(Ekomplexni)

% Prechodne vnitrni el.-mot. napeti E'
Eckomplexni=Unf + i*Xsumac*(Igcinny-i*Igjalovy)
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=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
