clear;clc; format compact; close all;
disp(' ');
disp('             A   H Á R O M F Á Z I S Ú   F E L A D A T         ');
disp('         A generátor csillag, a fogyasztó delta kapcsolású     ');
disp(' ');

% kiindulasi adatok
ug=231;ug1=ug*exp(j*2*pi/3);ug2=ug;ug3=ug*exp(j*4*pi/3);f=50;w=2*pi*f;
zgv=.06+.1*j;L=.052;C=100e-6;R=27;zab=R+j*w*L;zbc=2*R;zca=R+1/(j*w*C);

disp(' ');
disp(' A hurokáramok módszerével az ágáramok [A] ');
disp(' ');

z=[2*zgv+zab    -zgv     -zab
    -zgv      2*zgv+zbc  -zbc
    -zab        -zbc   zab+zbc+zca];
u=[ug1-ug2;ug2-ug3;0]; % A pontosvessző ; oszlopvektort eredményez!
h=inv(z)*u;
iab=h(1)-h(3);ibc=h(2)-h(3);ica=-h(3);%a fázisáramok
ia=h(1);ib=h(2)-h(1);ic=-h(2);%a vonaláramok
uab=iab*zab;ubc=ibc*zbc;uca=ica*zca;%a fogyasztó fázis- és vonalfeszültségei
ifh=[iab;ibc;ica]
ifhabs=abs(ifh)
ivh=[ia;ib;ic]
ivhabs=abs(ivh)

disp(' ');disp(' ');
disp(' A csomóponti potenciálok módszerével az ágáramok [A] ');
disp(' ');

y=[1/zgv+1/zca+1/zab    -1/zab        -1/zca
     -1/zab      1/zgv+1/zab+1/zbc    -1/zbc
     -1/zca             -1/zbc     1/zgv+1/zbc+1/zca];
i=[ug1/zgv;ug2/zgv;ug3/zgv];
up=inv(y)*i;
ifp=[(up(1)-up(2))/zab;(up(2)-up(3))/zbc;(up(3)-up(1))/zca]
ifpabs=abs(ifp)
ivp=[(ug1-up(1))/zgv;(ug2-up(2))/zgv;(ug3-up(3))/zgv]
ivpabs=abs(ivp)

% fazorábrák

Uab=ug1-ug2; Ubc=ug2-ug3; Uca=ug3-ug1;
Ugfv=[ug1;ug2;ug3;Uab;Ubc;Uca;uab;ubc;uca];
figure;
compass(Ugfv)
text(real(ug1),imag(ug1),'ug1');text(real(ug2),imag(ug2),'ug2');
text(real(ug3),imag(ug3),'ug3');text(real(Uab),imag(Uab),'Uab');
text(real(Ubc),imag(Ubc),'Ubc');text(real(Uca),imag(Uca),'Uca');
text(real(uab),imag(uab),'uab');text(real(ubc),imag(ubc),'ubc');
text(real(uca),imag(uca),'uca');
title('A generátor fázis és vonali-, továbbá a fogyasztó fázis- feszültsége [V]');
zoom on;
grid on;

fazorv=[uab/20;ubc/20;uca/20;ia;ib;ic];
figure;
compass(fazorv);
text(real(uab/20),imag(uab/20),'uab/20');text(real(ubc/20),imag(ubc/20),'ubc/20');
text(real(uca/20),imag(uca/20),'uca/20');
text(real(ia),imag(ia),'ia');text(real(ib),imag(ib),'ib');
text(real(ic),imag(ic),'ic');
title('A vonali feszültségek [V] és áramok [A] a fogyasztóra');
zoom on;
grid on;

fazorf=[uab/20;ubc/20;uca/20;iab;ibc;ica];
figure;
compass(fazorf);
text(real(uab/20),imag(uab/20),'uab/20');text(real(ubc/20),imag(ubc/20),'ubc/20');
text(real(uca/20),imag(uca/20),'uca/20');
text(real(iab),imag(iab),'iab');text(real(ibc),imag(ibc),'ibc');
text(real(ica),imag(ica),'ica');
title('A fázis feszültségek [V] és áramok [A] a fogyasztóra');
zoom on;
grid on;

% további eredmények

disp(' ');disp(' ');
disp(' A fogyasztó fázis feszültsége [V] ');
disp(' ');
uf=[uab;ubc;uca]
ufabs=abs(uf)

disp(' ');disp(' ');
disp(' A generátor fázis feszültsége [V] ');
disp(' ');
ugf=[ug1;ug2;ug3]
ugfabs=abs(ugf)

disp(' ');disp(' ');
disp(' A generátor vonali feszültsége [V] ');
disp(' ');
Ugv=[Uab;Ubc;Uca]
Ugvabs=abs(Ugv)

disp(' ');disp(' ');
disp(' A generátor vonali feszültségének időeltolása a fázis feszültségéhez képest [s] (késés = +) ');
disp(' ');
Tfv=[(angle(ug1)-angle(Uab))/2/pi/f;(angle(ug2)-angle(Ubc))/2/pi/f;(angle(ug3)-angle(Uca))/2/pi/f]

disp(' ');disp(' ');
disp(' A fugyasztó feszültségének időeltolása a generátor vonali feszültségéhez képest [s] (késés = +) ');
disp(' ');
Tvv=[(angle(Uab)-angle(uab))/2/pi/f;(angle(Ubc)-angle(ubc))/2/pi/f;(angle(Uca)-angle(uca))/2/pi/f]

disp(' ');disp(' ');
disp(' A fugyasztó vonali áramának időeltolása a feszültségéhez képest [s] (késés = +) ');
disp(' ');
Tiv=[(angle(uab)-angle(ia))/2/pi/f;(angle(ubc)-angle(ib))/2/pi/f;(angle(uca)-angle(ic))/2/pi/f]

disp(' ');disp(' ');
disp(' A fugyasztó fázis áramának időeltolása a feszültségéhez képest [s] (késés = +) ');
disp(' ');
Tif=[(angle(uab)-angle(iab))/2/pi/f;(angle(ubc)-angle(ibc))/2/pi/f;(angle(uca)-angle(ica))/2/pi/f]

disp(' ');disp(' ');
disp(' A teljesítmények [W] ');
disp(' ');
Pfogy=((abs(iab)^2)*R+(abs(ibc)^2)*2*R+(abs(ica)^2)*R)/2    % "/2" a csúcsértékek miatt!
Pveszt=(real(zgv)*((abs(ia)^2)+(abs(ib)^2)+(abs(ic)^2)))/2  % "/2" a csúcsértékek miatt!
PR=Pfogy+Pveszt
Pg=-(real(ug1*(conj(ia))+ug2*(conj(ib))+ug3*(conj(ic))))/2  % "/2" a csúcsértékek miatt!