W ramach przygotowania do mini-projektu należy przeczytać rozdziały ze skryptu R: 1, 2.1, 2.2.1, 2.2.2, 2.2.4.

http://anna.czemplik.staff.iiar.pwr.wroc.pl/images/Dmodele/zz_cz1.pdf

Zapisać: % zmienne wejściowe % zmienne wyjściowe Schemat I-7

zz_cz1.pdf

ob1_skrypt1.m

% ob1_skrypt1.m
%======= I część (identyfikacja)=======
% wartości nominalne
TzewN=-20;   % °C
TwewN=20;    % °C
TgN=45;      % °C
QgN=1500;    % W
Vg=10/1000;           % m3, objetosc grzenika
Vw-2.5*10*10;         % m3, objetosc pokoju
cpp=1000; rop=1.2;    % J/(kgK), kg/m3, powietrze
cpo=2400; roo=1200;   % J/(kgK), kg/m3, olej

%--------------------------------------
% identyfikacja parametrów statycznych
FpN=Vw/ (6*60*60);
Kg=QgN/ (TgN-TwewN);
K1=(QgN-cpp*rop*FpN* (TwewN-TzewN))/(TwewN-TzewN);

%--------------------------------------
% parametry "dynamiczne"
Cvg=cpo*roo*Vg;
Cvw=cpp *rop*Vw;

%======= II część (punkt pracy)========
% warunki początkowe
Tzewo = TzewN + 0;  % 0, 10, 20
Qg0 = QgN 1.0;      % 1.0, 0.7, 0.3
Fp0 = FPN 1.0;      % 1.0, 0.7, 0.3

%--------------------------------------
% stan równowagi
Twew0= Qg0/ (K1+cpp*rop*Fp0) +Tzew0;
Tg0 Qg0/Kg + Twew0;

%======= III część (symulacje)=========
czas = 50000; % czas symulacji
% zakłócenie
czas_skok = 5000;
dTzew=1;  % 0 lub np. 1
dQg=0;    % 0 lub np. 0.05*QgN (5% zakresu)
dFp=0;    % 0 lub np. 0.2*FpN (20% zakresu)

%--------------------------------------
% symulacja
[t]=sim ('obl', czas); % t-wektor czasu

%--------------------------------------
% wykresy
figure, plot (t, aTwew, 'r'), grid on, title ('Reakcja Twew na skok Tzew');
figure, plot (t, aTg, 'r'), grid on, title('Reakcja Tg na skok Tzew');