close all
% Importovat z přílohy č. 2 z druhého listu sloupce A a B
% Prejmenovat na struct s hodnotami ze sloupce A na T_ven_cas, ze sloupce B
% na T_ven_hodnota
% Spustit před simulací v simulinku

% % -------------------------
T_poc = 20;                                 %[°C]
T_poc_steny = 21;                           %[°C]
T_poz = [21 21 22 22 20 20 23 23];          %[°C]
Vykon_rad = 2061;
Step_time = 100;

% Regulátor
Kr = 11.612;
Ti = 1062.269;
I=0.1;

% Vzduch uvnitř místnosti
V_vzduch = 75;                                  %[m3] 
ro_vzduch = 1.2047;                             %[kg/m3]
c_vzduch = 1010;                                %[J/(kg*K)]
m_vzduch = ro_vzduch*V_vzduch;                  %[kg]

% Okno
d_okno = 0.044;                                 %[m]
S_okno = 1.5*1.5;                               %[m2]
V_okno = d_okno*S_okno;                         %[m3]
m_okno = 29;                                    %[kg]
c_okno = 850;                                   %[J/(kg*K)]
lambda_okno = 0.078;                            %[W/(m*K)]
U_okno = lambda_okno/d_okno;                    %[W/(m2*K)]

% Dveře vnitřní
d_dvere = 0.05;                                 %[m]
m_dvere = 30;                                   %[kg]
c_dvere = 1300;                                 %[J/(kg*K)]
lambda_dvere = 0.32;                            %[W/(m*K)]
S_dvere = 2.1105;                               %[m2]
U_dvere = lambda_dvere/d_dvere;                 %[W/(m2*K)]

% Strop beton
d_beton = 0.25;                                 %[m]
ro_beton = 2300;                                %[kg/m3]
c_beton = 1200;                                 %[J/(kg*K)]
lambda_beton = 1.36;                            %[W/(m*K)]
V_beton = 5*5*d_beton;                          %[m3]
m_beton = V_beton*ro_beton;                     %[kg]
U_beton = lambda_beton/d_beton;                 %[W/(m2*K)]

% Strop polystyren
d_poly = 0.15;                                  %[m]
lambda_poly = 0.044;                            %[W/(m*K)]
V_poly = 25*d_poly;                             %[m3]
U_poly = lambda_poly/d_poly;                    %[W/(m2*K)]

% Strop
m_strop = m_beton;                              %[kg]
c_strop = c_beton;                              %[J/(kg*K)]
S_strop = 25;                                   %[m2]
U_strop = 1/(d_poly/lambda_poly+d_beton/lambda_beton);%[W/(m2*K)]

% Stena vnejsi
d_zed_ven = 0.44;                               %[m]
V_zed_ven = 10*3*d_zed_ven-S_okno*d_zed_ven;    %[m3]
lambda_zed_ven = 0.165;                         %[W/(m*K)]
c_zed_ven = 1000;                               %[J/(kg*K)]
ro_zed_ven = 750;                               %[kg/m3]
m_zed_ven = ro_zed_ven*V_zed_ven;               %[kg]
S_zed_ven = 3*10-S_okno;                        %[m2]
U_zed_ven = lambda_zed_ven/d_zed_ven;           %[W/(m2*K)]

% Stena vnitrni
d_zed_vni = 0.14;                               %[m]
V_zed_vni = 5*3*d_zed_vni-S_dvere*d_zed_vni;    %[m3]
lambda_zed_vni = 0.26;                          %[W/(m*K)]
c_zed_vni = 1000;                               %[J/(kg*K)]
ro_zed_vni = 750;                               %[kg/m3]
m_zed_vni = ro_zed_ven*V_zed_vni;               %[kg]
S_zed_vni = 3*5-S_dvere;                        %[m2]
U_zed_vni = lambda_zed_vni/d_zed_vni;           %[W/(m2*K)]
