f=100
w=2*pi*f
Uf=100
Zn=30*i
Za=30+10*i
Zb=20
Zc=20+10*i
a=-0.5+0.866*i
t=0:0.0001:0.03

ua=Uf*sqrt(2)*sin(w*t)
ub=Uf*sqrt(2)*sin(w*t-2*pi/3)
uc=Uf*sqrt(2)*sin(w*t+2*pi/3)
f=100
w=2*pi*f
Uf=100
Zn=30*i
Za=30+10*i
Zb=20
Zc=20+10*i
a=-0.5+0.866*i
t=0:0.0001:0.03
ia=11.5*sqrt(2)*sin(w*t-0.16*pi)
ib=11.13*sqrt(2)*sin(w*t-0.655*pi)
ic=15.7*sqrt(2)*sin(w*t+0.58*pi)
ua=Uf*sqrt(2)*sin(w*t)
ub=Uf*sqrt(2)*sin(w*t-2*pi/3)
uc=Uf*sqrt(2)*sin(w*t+2*pi/3)

plot(t,ua)
hold on
plot(t,ub)
hold on
plot(t,uc)
hold on
plot(t,ia)
hold on
plot(t,ib)
hold on
plot(t,ic)
hold on
Ea=Uf
Eb=a*a*Uf
Ec=a*Uf

Ya=1/Za
Yb=1/Zb
Yc=1/Zc
Yo=1/Zn

Uoo=(Ea*Ya+Eb*Yb+Ec*Yc)/(Ya+Yb+Yc)


%Układ gwiazda brak przewodu zerowego
Ia=(Ea-Uoo)*Ya
Ib=(Eb-Uoo)*Yb
Ic=(Ec-Uoo)*Yc
faza_a=atan(imag(Ia)/real(Ia))
faza_b=atan(imag(Ib)/real(Ib))
faza_c=atan(imag(Ic)/real(Ic))

%Jest przewód zerowy
Uoo=(Ea*Ya+Eb*Yb+Ec*Yc)/(Ya+Yb+Yc+Yo)

Ia=(Ea-Uoo)*Ya
Ib=(Eb-Uoo)*Yb
Ic=(Ec-Uoo)*Yc
Ioo=Ia+Ib+Ic


faza_a=atan(imag(Ia)/real(Ia))
faza_b=atan(imag(Ib)/real(Ib))
faza_c=atan(imag(Ic)/real(Ic))


% Trójkąt
Ia=(Ea-Eb)*Ya-(Ec-Ea)*Yc
Ib=(Eb-Ec)*Yb-(Ea-Eb)*Ya
Ic=(Ec-Ea)*Yc-(Eb-Ec)*Yb


faza_a=atan(imag(Ia)/real(Ia))
faza_b=atan(imag(Ib)/real(Ib))
faza_c=atan(imag(Ic)/real(Ic))