%% Parametry obiektu
TwN = 20;                   %temp. wewn. budynku
TzewN = -20;                %temp. na zewn. budynku
TgzN = 90;                  %temp. wpływająca do grzejnika
TgpN = 70;                  %temp. wypływająca z grzejnika

Cpp = 1000;                 %ciepło właściwe powietrza
Cpw = 4175;                 %ciepło włlaściwe wody
ep = 1.2;                   %gęstość powietrza
ew = 1000;                  %gęstość wody

Vp = 3 * 4 * 2.3;           %objętość pomieszczenia
Vw = 2 * 0.7 *1;            %objętość grzejnika
Vk = 1 * 1 * 1;             %objętość kotła

QN = 5000;                  %ciepło oddawane przez grzejnik
QdN = 0;
PkN = 3 * QN;               %moc kotła na 3 pomieszczenia

T0=1000;                    %opóźnienie 
Tdelay = 1000;

%% warunki początkowe (nominalne)
Tzew0=TzewN;
Tgz0=TgzN;
Tkz0=TkzN;
Tkp0=TkpN;

F0 = FN;
Qd0=QdN;
Pk0 = PkN;

%% zakłócenia
dTzew0=0;
dTgz0=0;

dF1=0;
dF2=0;
dF3=0;

dQd1=0;
dQd2=500;
dQd3=0;

dPk0=0;

%% Obliczenia dla wartości nominalnych
cvg=Cpw*ew*Vw;           %pojemność pieplna grzejnika
cvw=Cpp*ep*Vp;           %pojemność pieplna pokoju
cvk=Cpw*ew*Vk;           %pojemność pieplna kotła

TgsrN = (TgzN + TgpN) / 2;
TgpN = 2 * TgsrN - TgzN;
TkpN=TgpN;
TkzN=TgzN;
FN = QN / (2 * Cpw * ew * (TgzN - TgsrN));
FkN = 3 * FN;
kg = QN / (TgsrN - TwN);
kstr = (QN + QdN) / (TwN - TzewN);

%% Obliczenia cd.
Tgsr0 = ( 2 * Cpw * ew * F0 * Tgz0 * (kg + kstr) + kg * (kstr*Tzew0 + Qd0) ) / ( kg * kstr + 2 * Cpw * ew * F0 * (kg + kstr) );
Tw0 = (Qd0 + kg * Tgsr0 + kstr * Tzew0) / (kg + kstr);



