.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(..
29 'SearchMethod', 'auto');
42
....
Níže uveden příklad kódu, jeho podrobné vysvětlení výsledky aproximace pro
topení obývacím pokoji poruchový vstup venkovní teploty. řádu,
přičemž původní model bez zkrácení dynamiky měl řád 228.
26 'InitialCondition', 'zero', .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. 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.4 Aproximace
Pro návrh řízení teploty bylo nutné zjednodušit původní model systému analyzovat
chování jednotlivých zón.
28 'Display', 'on', .4.. Pro každou nich byla
vykreslena frekvenční charakteristika, která následně sloužila jako podklad pro apro-
ximaci modelem nižšího řádu.Z{1}) počet nul
14 length(Fs_cont_tf_mr_bezI....
27 'EnforceStability', true,..Tímto způsobem
bylo získáno celkem přenosových funkcí reprezentujících tepelný systém 196. 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', .
Výpis 4
