
%%analyticky
%Model B
clear
syms p Ui

R6 = 2;
C1 = 0.04;
R1 = 50;
L2 = 2;
Res2 = 110;
C2 = 0.0005;
R2 = 1300;
L3 = 0.95;
Res3 = 60;
C3 = 0.0012;
R3 = 550;
Z1 = R6+(R1*(1/(p*C1)))/(R1+(1/(p*C1)));
Z2 = ((Res2*(p*L2))/(Res2+(p*L2)))+((R2*(1/(p*C2)))/(R2+(1/(p*C2))));
Z3 = ((Res3*(p*L3))/(Res3+(p*L3)))+((R3*(1/(p*C3)))/(R3+(1/(p*C3))));
Z4 = (Z2*Z3)/(Z2+Z3);
Z = (Z1*Z4)/(Z1+Z4);
I = Ui/Z;
I1 = I*(Z1/(Z1+Z4));
I2 = I1*(Z2/(Z3+Z2));
I3 = I2*((1/(p*C3))/((1/(p*C3))+R3));
U3 = R3*I3;
%Pøenos ze srdce na nohu
noha_srdce = U3/Ui;
noha_srdce = simplify(noha_srdce);
pretty(noha_srdce)

%Pøenos ze srdce na ruku
I20 = I1-I2;
I21 = I20*((1/(p*C2))/((1/(p*C2))+R2));
U2 = I21*R2;
ruka_srdce = U2/Ui;
ruka_srdce = simplify(ruka_srdce);
pretty(ruka_srdce)

%Pøenos z ruky (U2) na nohu (U3)
I20 = I1-I2;
I21 = I20*((1/(p*C2))/((1/(p*C2))+R2));
U2 = I21*R2;
ruka_noha = U3/U2;
ruka_noha = simplify(ruka_noha);
pretty(ruka_noha)


%Model E
%Parametry jednotlivých prvkù
clear
syms p Ui % p = j*omega, Ui = vstupni napeti
Ra = 5;
C1 = 7e-3;
R1 = 120;
R2 = 15;
L2 = 3.6;
C2 = 1e-3;
R21 = 1035;
Rb = 5;
R3 = 45;
C3 = 1e-2;
R4 = 60;
L4 = 0.8;
C4 = 4e-3;
R41 = 465;

%Impedance frekvenènì závislých prvkù prvkù v Laplaceove oblasti
RC1 = 1/(p*C1);
RL2 = p*L2;
RC2 = 1/(p*C2);
RC3 = 1/(p*C3);
RL4 = p*L4;
RC4 = 1/(p*C4);

%Impedance jednotlivých smyèek -    metoda postupného zjednodušování
Z1 = (R1*RC1)/(R1+RC1);
Z21 = (R2*RL2)/(R2+RL2);
Z22 = (R21*RC2)/(R21+RC2);
Z2  = Z21+Z22; 
Z3 = (R3*RC3)/(R3+RC3);
Z41 = (R4*RL4)/(R4+RL4);
Z42 = (R41*RC4)/(R41+RC4);
Z4 = Z41+Z42;
Z5 = (Z1*Z2)/(Z1+Z2);
Z6 = (Z3*Z4)/(Z3+Z4);
Z7 = Rb+Z6;
Z8 = (Z5*Z7)/(Z5+Z7);
Z = Ra+Z8;
I = Ui/Z;

%Dìlièe proudu v jednotlivých uzlech - krouzkovano cervene
I1 = I*(Z5/(Z5+Z7));
I2 = I1*(Z3/(Z3+Z4));
I3 = I2*(RC4/(RC4+R41)); %proud tekoucí odporem R41

%Napìtí na prvku R41 - místo kde hledám pøenos vùèi Uin
UR41 = I3*R41;

%Pøenosová funkce
noha_srdce = UR41/Ui;
noha_srdce = simplify(noha_srdce);
pretty(noha_srdce) 

%%pomoci funkce tfest
tfdata = iddata(vstup,vystup);
sys = tfest(tfdata,rad_modelu)

%%analyza pomoci Linear analysis je dominantou Simulinku
