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(.4..
26 'InitialCondition', 'zero', .Z{1}) počet nul
14 length(Fs_cont_tf_mr_bezI..
28 'Display', 'on', . 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)..
29 'SearchMethod', 'auto');
42
.
Výpis 4.
Níže uveden příklad kódu, jeho podrobné vysvětlení výsledky aproximace pro
topení obývacím pokoji poruchový vstup venkovní teploty..
25 'InitializeMethod', 'all', .4 Aproximace
Pro návrh řízení teploty bylo nutné zjednodušit původní model systému analyzovat
chování jednotlivých zón... 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.
27 'EnforceStability', true,. 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.Tímto způsobem
bylo získáno celkem přenosových funkcí reprezentujících tepelný systém 196...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
