Simulácia riadenia TITO systému adaptívnym STC regulátorom

Riadenie dvojrozmerného (TITO - Two Inputs Two Outputs) systému pomocou STC regulátora vychádza z nasledujúcej principiálnej schémy:


Táto schéma znázorňuje zapojenie uzavretého regulačného obvodu, v ktorom blok "Regulovaný systém" má v tomto prípade dva vstupy, dva výstupy a nemá poruchu v. K tomuto systému je paralelne zapojený blok priebežnej rekurzívnej identifikácie, ktorý v každom kroku riadenia odhaduje parametre riadeného TITO systému a tie zapisuje do bloku "Výpočet parametrov regulátora". Vypočítané parametre regulátora q sa v bloku "Regulátor" použijú na výpočet akčného zásahu, teda výstupu regulátora. To prebieha cyklicky pre každý krok riadenia.
Simulačná schéma riadenia TITO systému STC regulátorom je potom zostavená podľa tejto principiálnej schémy v jazyku Simulink z blokov knižnice STCLIB podľa nasledujúceho obrázku:


V tejto simulačnej schéme blok STC regulátora "PP_TITO_a" obsahuje blok "Výpočet parametrov regulátora" a "Regulator" znázornené v predošlej schéme.
Spustením simulácie tejto schémy v jazyku Simulink pri zadaných nastaveniach blokov získajú sa priebehy jednotlivých veličín obvodu, ktoré sa znázorňujú v grafoch.
Príklad riadenia TITO systému STC regulátorom
Podmienky simulácie:
Riadený systém je popísaný maticami:


• Perióda vzorkovania Tvz = 1 [s].
• Doba trvania simulácie Stop_Time = 200 [s].
• Diagonálna matica M(z-1), ktorou sa rozmiestnia póly STC regulátora, je:

kde:

• Počiatočný odhad parametrov ARX modelu TITO systému pre RMNŠ je:

• Použitá je identifikačná metóda so smerovým zabúdaním:
C(0) = 103, j(0) = 1, l(0) = 0.001, u(0) = 10-6, r(0) = 0.99
• Obmedzenie akčného zásahu regulátora je +10 a -10.
• Riadiaca hodnota w1(k), w2(k) je generovaná so skokovými zmenami a so sínusovým priebehom (ako je zrejmé z obrázkov simulácií).
• Na vstup systému pôsobí skoková a vysokofrekvenčná sínusová porucha podľa obrázka.



• A na výstupe systému pôsobí aditívny, dynamikou systému filtrovaný, biely šum zobrazený na obrázku.



• V k = 180 boli zmenené parametre systému: b2 z 0.4 na hodnotu 0,5 a b3 z 0.9 na 0,8 a b6 z hodnoty 0,2 na 0.25.

Získané priebehy:







Nasledujúce formuláre umoňňujú zmeniť niektoré parametre simulácie a simuláciu znova spustiť. Jej výsledky sa zobrazia na vyššie zobrazených grafoch.

Parametre STC regulátora TITO systému - " PP_TITO_a "

Koeficienty matice M :
Obmedzenie výstupu u1(k) regulátora - " Saturation1 "
Horné : Dolné :
Obmedzenie výstupu u2(k) regulátora - " Saturation2 "
Horné : Dolné :
V bloku STC regulátora možno nastaviť parametre matice predpísaných pólov M a to tak, že v hranatých zátvorkách do vstupného poľa vložíte 8 koeficientov ( [m1 m2 ... m8] ). Matica M je v tvare:

Parametre diskrétneho modelu systému - " System TITO "

Koeficienty matice A pred zmenou parametra :
Koeficienty matice B pred zmenou parametra :
Koeficienty matice A po zmene parametra :
Koeficienty matice B po zmene parametra :
Čas, kedy má nastať zmena parametra : [s]
V bloku diskrétneho systému sa matice A a B zadávajú pomocou 8-ich koeficientov zadaných v hranatých zátvorkách ( [ a1 a2 ... a8 ] , [ b1, b2, ... b8 ] ). Matice A a B sú:


V nastavenom čase simulácie sa zmenia parametre v maticiach A, B na hodnoty uvedené v maticiach určujúcich koeficienty A, B po zmene parametra.

Parametre simulácie

Perióda vzorkovania : [s]
Čas začiatku simulácie (Start time) : [s]
Čas ukončenia simulácie (Stop time) : [s]
Perióda vzorkovania platí pre celú simulačnú schému.

Parametre rekurzívnej identifikačnej metódy najmenších štvorcov
- " Identifikacia TITO "

Rád identifikovaného modelu :
Typ identifikačnej metódy :
Počiatočný odhad parametrov modelu - On1 :
Počiatočný odhad parametrov modelu - On2 :
Počiatočná kovariančná matica - Pn0 :
Počiatočná hodnota zabúdania - ex0 :
Počiatočná lambda - lambda0:
Počiatočná teta - teta0:
Pomocný skalár - ro:
Parametrom "Rád identifikovaného modelu" možno určiť, či sa budú identifikovať len parametre na diagonále matice A, alebo všetky.
"Typ identifikačnej metódy" môže byť: LSM ez - RMNŠ s exponenciálnym zabúdaním, alebo LSM až s adaptívnym smerovým zabúdaním.
On1 sa zadáva v hranatých zátvorkách a je v tvare: [ a1 a2 a3 a4 b1 b2 b3 b4 ], tieto koeficienty predstavujú prvý riadok matíc A, B uvedených vyššie.
On2 sa zadáva podobne: [ a5 a6 a7 a8 b5 a6 a7 a8 ], ako druhý riadok matíc A, B.

Parametre generátora signálu vstupnej poruchy - " porucha "

Generátor sínusového signálu 1 - " Sine Wave1 "
Amplitúda : Frekvencia : [Hz]
Pripojenie poruchy k 1. vstupu systému - " Manual Switch1 "
--------------------------------------------------------------------------------------------------
Generátor sínusového signálu 2 - " Sine Wave2 "
Amplitúda : Frekvencia : [Hz]
Pripojenie poruchy k 2. vstupu systému - " Manual Switch2 "
Posledný priebeh vstupnej poruchy:


Vnútorná schéma bloku vstupnej poruchy - blok "porucha":

Parametre generátora signálu výstupnej poruchy - " sum+porucha "

Úroveň bieleho šumu 1 - " Gain1 "
Zosilnenie :
Pripojenie poruchy k 1. výstupu systému - " Manual Switch1 "
--------------------------------------------------------------------------------------------------
Úroveň bieleho šumu 2 - " Gain2 "
Zosilnenie :
Pripojenie poruchy k 2. výstupu systému - " Manual Switch2 "
Posledný priebeh výstupnej poruchy:


Vnútorná schéma bloku výstupnej poruchy - blok "sum+porucha":
Blok - "System TITO" je automaticky nastavený podľa nastavenia matíc A pred a po zmene parametra, pričom matica B je jednotková. Tento systém tu slúži ako filter. Ak je táto porucha pripojená k výstupu riadeného systému, identifikuje sa ARX model systému.

Teraz môžete spustiť nový simulačný výpočet modelu riadenia TITO systému STC regulátorom stlačením tlačidla "Spust výpocet".