clear all
close all
clc

Uref= 231; %napätie referenčné

Z= [2 3 2 1 0.5 5.5 1.5 2.5]; % impedancie vo vetvách vetiev

Y= inv(diag(Z)); % diagonálna matica vetvových admitancií

Kr= [ 1 0 -1 0 0 0;  % redukovaná prvá incidenčná matica
 1 0 0 -1 0 0; 
0 0 -1 0 0 0; 
0 0 -1 0 1 0;
0 0 0 -1 1 0;
0 0 0 0 -1 1;
0 0 0 1 0 -1;
0 1 0 -1 0 0;
;];

Izat= [ -5; % stĺpcová matica uzlových zaťažových prúdov
 -55;
 -39;
 -40;];

Yadm= transpose(Kr)*Y*Kr % uzlová admitančná matica
Y11= Yadm(1:2,1:2);     %vzdy podla Yadm (pocet nap vs odber)
Y12= Yadm(1:2,3:6);
Y21= Yadm(3:6,1:2);
Y22= Yadm(3:6,3:6);

Unap= [231-Uref;
    225-Uref;
       
       ]; % stĺpcová matica uzlových napätí v napájacích uzloch

Uzat= inv(Y22) * (Izat - Y21*Unap) % napätia v zaťažných uzloch

U= [Unap; % matica uzlových napätí
Uzat; ]

%Inap= Y11 * Unap + Y12 *inv(Y22) * (Izat-Y21*Unap) % prúdy v zaťažných uzloch

%Icelkove= [Inap; % matica uzlových prúdov
%Izat; ]
%Un= inv(Yadm)*Icelkove

deltaU= Kr*U % matica úbytkov napätí vo vetvách
I= Y*deltaU % matia prúdov vo vetvách