Původně byla tato inicializace realizována
přivedením záporné sumy proporcionální diferenční složky.J.3). však třeba
poznamenat, simulační model použitý při optimalizaci zahrnuje určité vnitřní
zpoždění: konkrétně zpoždění 0. tohoto důvodu byla inicializace upravena přivedení záporné sumy pouze
proporcionální složky hodnoty, která byla naladěna ručně metodou pokus–omyl
53
. Při výpočtu ampli-
tudy přenos |𝐺(𝑗𝜔řez)| popisuje systém bez integrační složky.
Výpis 5.slx");
4 out.3)
Touto rovnicí nastavujeme zisk regulátoru tak, aby otevřená smyčka měla přesně
jednotkové zesílení místě, kde jsme nastavili frekvenci řezu. každém kroku spustí simulační model zahrnující systém bez modelova-
ného dopravního zpoždění, jelikož při návrhu regulátoru byl kódu použit model se
zpožděním pro výpočet frekvenční odezvy systému (viz Obrázek 5. kódu této apro-
ximace dosaženo násobením 𝜔řez, což uvedeném vztahu již neprojevuje.Zesílení regulátoru vypočítáno podle vztahu:
𝐾 =
1
|𝐺(𝑗𝜔řez)|
·
√︁
𝜀2𝜔2
řez +1
√︁
𝑇2
1 𝜔2
řez +1·
√︁
𝑇2
2 𝜔2
řez +1
(5. Ukázalo však, že
tento přístup vede nevyvážené dynamice regulátoru diferenční složka díky své
rychlé dynamice velmi rychle odezní, zatímco sumační složka zůstane ovlivněna po-
čátečním impulzem. Dále
probíhá optimalizace parametru který ovlivňuje chování regulátoru při změně žá-
dané hodnoty.2: Ukázkový kód funkci pro výpočet ITAE kritéria pro optimalizaci v
jazyce MATLAB
1 function opt_beta(x)
2 assignin("base","beta",x);
3 out sim("regulator_beta_optimalizace. Po
určení všech parametrů regulátor sestaven tvaru ZPK (nuly, póly, zesílení) a
následně převeden standardního tvaru PID parametry 𝐾𝑝, 𝑇𝑖, odvozenou
back-calculation sledovací konstantou =
√
𝑇𝑖 ·𝑇𝑑 pro ošetření wind-up jevu. Pro zrychlení
reakce regulátoru snížení kmitů průběhu akčního zásahu byla vstup sumace
přidána inicializace počáteční hodnoty. Hodnocení kvality
řízení založeno ITAE kritériu, které hodnotí regulační plochu vážením podle
času. Výsledkem pak výrazný překmit teploty nesmyslných hod-
not.Data;
5 end
Limitační blok pro akční zásah nastaven rozsahu kW. optimalizaci používá funkce "fminsearch", která hledá takovou
hodnotu která zajistí minimální překmit přechodného děje. Model vypočte hodnotu kritéria ta
se numericky minimalizuje prostřednictvím hodnoty 𝛽.5𝑇𝑠 způsobené blokem Zero-Order Hold dále zpož-
dění 1𝑇𝑠 vzniklé použitím bloku Unit Delay
