%% Delutani szekcio - multkori folytatasa
% 1. Elso pelda - halozati egyenletek megoldasa  
clear all
clc
R = 2.4; % kOhm
U0 = 10; % V
alfa = 0.78; % erosites
A = [3 -1 0 -2*R*alfa;-1 5 -2 2*R*alfa;0 -1 4 0;0 -1 1 -R]
B = [0;0;U0;0]
x = A \ B
% helyettesito generator parametereinek szamitasa
% szakadasi feszultseg - UT
A2 = [3 -1 0 -2*R*alfa;-1 5 -2 2*R*alfa;0 -1 2 0;0 -1 1 -R]
x2 = A2\B

% rovidzarasi aram - IN
A3 = [3 -1 0 -2*R*alfa;-1 5 0 2*R*alfa;0 0 R R;0 1 0 R]
det(A)
B = [0;0;U0;0]
x3 = A3 \ B
% belso ellenallas szamitasa
RB = 5.7627/4.1667
% maximalis teljesitmeny, ha Rt=RB
5.7627^2/(4*1.383)

% tetszoleges lezaras eseten a halozat megoldasa
% csak az az egyenlet(ek) valtozik, amely(ek) az Rt-hez kapcsolodo
% csomopont(ok)ra vonatkoznak
Rt=2; 
A = [3 -1 0 -2*R*alfa;-1 5 -2 2*R*alfa;0 -1/R 2/R+1/Rt 0;0 -1 1 -R];
x = A\B

%% 2. Ketkapu (lezart) megoldasa 
%  RR = [2 7;5 8]; 
% primer oldal UT=10 V, Rg=2 kOhm; 
% szekunder oldal Rt=10 kOhm
% valtozok : u1,i1,u2,i2 
format compact
A = [1 -2 0 0 ;0 0 1 -10;1 -2 0 -7;0 -5 1 -8]
B=[10;0;0;0]
x=A\B

%% 2/b. reciprocitas vizsgalata
% primer oldal U0=10 mA, szekunder U2 = 0 (rz)
A = [1 0 0 0 ;0 0 1 0;1 -2 0 -7;0 -5 1 -8]
B =[10;0;0;0];
x1 = A \ B
% megforditott eset 
% primer : rz (U=0), szekunder : gerj. (U=10)
A = [1 0 0 0 ;0 0 1 0;1 -2 0 -7;0 -5 1 -8];
B=[0;10;0;0];
x2 = A\B
% reciprok, ha x1(4) == x2(2)
%% 2/c. Reciprok ketkapu (R12=R21)
A = [1 0 0 0 ;0 0 1 0;1 -2 0 -5;0 -5 1 -8]
x2 = A\B
B=[10;0;0;0]
x1 = A \ B

%% 2/d. Szimmetrikus ketkapu (R12=R21 ES R11=R22)
A = [1 0 0 0 ;0 0 1 0;1 -8 0 -5;0 -5 1 -8]
x1 = A \ B
B=[0;10;0;0]
x2 = A\B