% Szerda - 2024.03.20 
% Ketkapu karakterisztika meghatarozasa 
R = 50;
r = 40;
% x = [U1, U2, Ua, Ia, Ub, Ib] ; s = [I1, I2]
A = [2/R -1/(2*R) -1/R 0 0 0;...
  -1/R 0 1/R 1 0 0;...
  -1/(2*R) 3/(2*R) 0 0 0 1;...
  0 0 1 0 0 r;...
  0 0 0 -r 1 0;...
  0 1 0 0 -1 0]
B = [1 0;0 0;0 1;0 0;0 0;0 0]
% Ismeretlenek ertekeknekm meghatarozasa
x = A \ B
% Ennek felso ket sora (U1,U2) adja az (I1,I2)-tol fuggo reszt
RR = x(1:2,1:2)

%% Atteres masik karakterisztikara
% Terjunk at G-re 
GG = inv(RR)

% Terjunk at K-ra (i1,u2) = K*(u1,i2) 
% A es B kozott kell oszlopokat cserelni 
% A(:,1) <--> B(:,1) 
% Vigyazni, hogy a csere (-1)-es szorzast jelent
AK = [-B(:,1) A(:,2:end)]
BK = [-A(:,1) B(:,2)]
xk = AK\BK
KK = xk(1:2,1:2)

%% ELlenorizzuk a G-t!
AG = [-B A(:,3:end)]
BG = -A(:,1:2)
xg = AG\BG;
GG2 = xg(1:2,:)
%% Lezaras hatasa
% Ketkapukent oldjuk meg mint a multoran


%% Lezaras hatasa 
% Teljes egyenletrendszerkent


