Níže uveden příklad kódu, jeho podrobné vysvětlení výsledky aproximace pro
topení obývacím pokoji poruchový vstup venkovní teploty..4 Aproximace
Pro návrh řízení teploty bylo nutné zjednodušit původní model systému analyzovat
chování jednotlivých zón.Z{1}) počet nul
14 length(Fs_cont_tf_mr_bezI.
27 'EnforceStability', true,.1: Ukázkový kód aproximaci přenosů jazyce MATLAB
1 Fs_cont d2c(Fs);
2 Převod transfer funkci
3 Fs_cont_tf zpk(Fs_cont);
4
5 Minimální realizace
6 tol 1e-4; tolerance
7 Fs_cont_tf_mr minreal(Fs_cont_tf,tol)
8
9 Odebrání integrátoru
10 Fs_cont_tf_mr_bezI minreal(Fs_cont_tf_mr*tf('s'),tol);
11
12 Výpočet relativního řádu
13 length(Fs_cont_tf_mr_bezI...Tímto způsobem
bylo získáno celkem přenosových funkcí reprezentujících tepelný systém 196.
28 'Display', 'on', ..
Výpis 4....
29 'SearchMethod', 'auto');
42
.p{1}) počet polů
15 np-nz
16
17 Výpočet frekvenční odezvy původního systému bez integrátoru
18 wstart 10e-8;
19 wstop 10e-2;
20 logspace(log10(wstart), log10(wstop), 5000);
21 sysgFs frd(Fs_cont_tf_mr_bezI, w);
22
23 Nastavení tfest
24 opttfest tfestOptions(. Aby bylo možné nahradit původní modely vysokého
řádu vhodnějšími aproximacemi, byly kombinovány funkce minreal (pro minimali-
zaci řádu), frd (frekvenční reprezentace) tfest (odhad parametrického modelu)...
25 'InitializeMethod', 'all', .4. Pomocí nástroje Model Linearizer prostředí MATLAB
Simulink byly stanoveny operátorové přenosy čtyř tepelných zón jednotlivé ří-
dicí (topení) poruchové (venkovní teplota, sluneční záření) vstupy.
26 'InitialCondition', 'zero', . Pro každou nich byla
vykreslena frekvenční charakteristika, která následně sloužila jako podklad pro apro-
ximaci modelem nižšího řádu.. řádu,
přičemž původní model bez zkrácení dynamiky měl řád 228
