close all
clear all
clc

Zimp= [2 2 1 0.5 5.5 1.5 3];   %Impedancia
Zv= diag(Zimp)      %diag. matica vetových impedancií


C= [-1 -1 -1 -1 -1 0 1;
    1 0 0 0 0 0 0;
    0 0 -1 -1 -1 0 1;
    0 0 0 0 -1 -1 1;
    0 0 0 -1 0 1 0;
    0 0 0 0 0 -1 0;
    0 0 0 0 0 0 -1;
%rozsírena druhá incidencna mat.
   
   ]
       
Zs= transpose(C)*Zv*C %matica sluckovych impedancii

Z11= Zs(1:1,1:1)
Z12= Zs(1:1,2:5)
Z13= Zs(1:1,6:7)
Z21= Zs(2:5,1:1)
Z22= Zs(2:5,2:5)
Z23= Zs(2:5,6:7)
Z31= Zs(6:7,1:1)
Z32= Zs(6:7,2:5)
Z33= Zs(6:7,6:7)

Izat=[-5; %matica zatazovych prudov
      -39;
      -40;
      -55;]

Unap= [231-231;
        ] %matica uzlových napajacích napätí

%matica napajacích prúdov
Inap= inv(1-inv(Z11)*Z13*inv(Z33)*Z31)*inv(Z11)*(Unap-(Z12-Z13*inv(Z33)*Z32)*Izat)

Ins= inv(Z33)*(-Z31*Inap-Z32*Izat) %matica prúdov v lin. nezávislých sluckach

Ics= [Inap; %celkova matica sluckovych prúdov
      Izat;
      Ins;]
  
Iv= C*Ics %prúdy vo vetvách

deltaU= Zv*Iv  %úbytky napätí vo vetvách