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

% Generator
Png=225
CosFi=0.9
Sng=Png/CosFi
Xd= (200/100)*(Un^2/Sng)*i
Xdc=( 30/100)*(Un^2/Sng)*i

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

% Vedeni
delka=30
Xv1=delka*0.4*i/110^2*Un^2
Xv=Xv1/2

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

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

% Vypocet predporuchoveho stavu, pocatecni zatizeni 100%
% Zvolit bud' velikost E, nebo CosFiS aktualni
% Zde volen CosFiS
Png0=1.00*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
Ey=    Igcinny *Xsuma
Ex=Unf+Igjalovy*Xsuma
E=sqrt(Ex^2+Ey^2)
Theta=atan(Ey/Ex)

% Prechodne vnitrni el.-mot. napeti E'
Eyc=    Igcinny *Xsumac
Exc=Unf+Igjalovy*Xsumac
Ec=sqrt(Exc^2+Eyc^2)
Thetac=atan(Eyc/Exc)

% Alternativni vypocet komplexnimy cisly
Ekomplexni=Unf + i*Xsuma*(Igcinny-i*Igjalovy)
Ekontrola=abs(Ekomplexni)
Thetakontrola=angle(Ekomplexni)

Eckomplexni=Unf + i*Xsumac*(Igcinny-i*Igjalovy)
Eckontrola=abs(Eckomplexni)
Thetackontrola=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
