clear;clc; format compact; close all; % Minden törlése

% Bemenő adatok

R=27;L=52e-3;C=100e-6;
Ug=100;Ig=2*exp(-j*5*pi/18);
f=500; w=2*pi*f;
Z1=R+j*w*2*L; Z2=j*w*L;
Z3=1/(j*w*C); Z4=R;
Z5=10*R+j*2*w*L;

% Kirchhoff egyenletekkel

Z=[ 1     0     1    0     0
   -1     1     0    0     1
    0    -1     0   -1     0
   -Z1    0     Z3   0    -Z5
    0    -Z2    0    Z4    Z5];
G=[Ig
    0
  -Ig
   Ug
    0];
I=inv(Z)*G
Iabs=abs(I)'
% VIGYÁZZ !!! A transzponálás ' komplex vektornál conjugálást is eredményez!

% Hurokáramokkal

ZH=[Z1+Z3+Z5       -Z5
      -Z5        Z2+Z4+Z5];
GH=[Z1*Ig+Ug
      Z2*Ig];
HA=inv(ZH)*GH;
IH=[Ig-HA(1)
    Ig-HA(2)
       HA(1)
       HA(2)
 HA(2)-HA(1)]
IHabs=abs(IH)'

% Csomoponti potenciálokkal

YP=[1/Z1+1/Z3     -1/Z1             0
    -1/Z1      1/Z1+1/Z2+1/Z5     -1/Z2
      0           -1/Z2         1/Z2+1/Z4];
GP=[-Ug/Z3+Ig
            0
          -Ig];
P=inv(YP)*GP;
IP=[(P(1)-P(2))/Z1
    (P(2)-P(3))/Z2
      (P(1)+Ug)/Z3
          -P(3)/Z4
           P(2)/Z5]
IPabs=abs(IP)'

% Az áramok megjelenítése

figure;
IP(6)=Ig;
compass(IP,'m');
text(real(IP(1)),imag(IP(1)),'i1');text(real(IP(2)),imag(IP(2)),'i2');
text(real(IP(3)),imag(IP(3)),'i3');text(real(IP(4)),imag(IP(4)),'i4');
text(real(IP(5)),imag(IP(5)),'i5');text(real(Ig),imag(Ig),'Ig');
title('Az egyes ágak áramai');
zoom on;
grid on;