PSMPC – Prediktivní regulátor typu pulse-step

Symbol blokuPotřebná licence: ADVANCED
Qt SVG Document Exported by REXYGEN Studio sp pv tv hv MAN mv dmv de SAT pve iE PSMPC

Popis funkce
Funkční blok PSMPC (Pulse Step Model Predictive Control) je určen pro realizaci vysoce kvalitních regulátorů pro obtížně regulovatelné lineární časově invariantní soustavy s omezením akční veličiny (např. soustavy s dopravním zpožděním nebo s neminimální fází). Zvlášť výhodný je pro případy, kdy je požadován velmi rychlý přechod z jedné hodnoty regulované veličiny na druhou bez překmitu. Regulátor PSMPC však může být obecně použit všude tam, kde je běžně nasazován standardní PID regulátor a kde žádáme vysokou kvalitu regulace.

image/svg+xml0 Ts 2Ts 3Ts 4Ts N Ts cas y(t) g(1)g(3)g(2)g(4)g(N)h(1)h(2)h(3)h(4)

PSMPC je prediktivní regulátor s explicitně zadaným intervalovým omezením akční veličiny. Pro účely predikce je použit model ve tvaru diskrétní přechodové charakteristiky g(j),j = 1,,N. Na obrázku výše je naznačen způsob, jakým lze tuto posloupnost získat ze spojité přechodové charakteristiky. Poznamenejme, že N musí být zvoleno dostatečně velké, aby přechodová charakteristika byla popsána až do ustáleného stavu (N TS > t95, kde TS je perioda vzorkování regulátoru a t95 je doba ustálení na 95 % konečné hodnoty).

Pro stabilní, lineární a t-invariantní systémy s monotónní přechodovou charakteristikou je alternativně možné použít momentový množinový model [6] a popsat systém pouze třemi charakteristickými čísly κ, μ a σ2 , které je možno určit z jednoduchého pulzního experimentu. Řízený systém pak aproximujeme buď přenosem prvého řádu s dopravním zpožděním

FFOPDT (s) = K τs + 1 eDs,κ = K,μ = τ + D,σ2 = τ2 (7.1)

nebo přenosem druhého řádu s dopravním zpožděním

FSOPDT (s) = K τs + 12 eDs,κ = K,μ = 2τ + D,σ2 = 2τ2 (7.2)

se stejnými charakteristickými čísly. Typ aproximace se zadává parametrem imtype.

Pro zjednodušení on-line optimalizace v otevřené smyčce je množina přípustných posloupností řízení omezena pouze na posloupnosti ve tvaru "pulz-skok"zobrazené na obrázku níže.

image/svg+xmlN c n 2 n 1 u u + u n 1 n 2 N c u u + u p 0 =1p 0 =0

Poznamenejme, že každá taková posloupnost je jednoznačně určena jen třemi čísly n1,n2 0,,NC a u u,u+, kde NC 0,1, je horizont řízení a u,u+ označují po řadě zadanou dolní a horní mez akční veličiny regulátoru. On-line optimalizace (vzhledem k n1,n2 a u) spočívá v minimalizaci kritéria

I = i=N1N2 ê(k + i|k)2 + λ i=0NC Δû(k + i k)2 min, (7.3)

kde ê(k + i|k) je v kroku k predikovaná regulační odchylka na intervalu predikce i N1,N2, Δû(k + i|k) jsou diference řídicího signálu na intervalu i 0,NC a λ je koeficient penalizace změn akční veličiny. Pro nalezení optima úlohy (7.3) je použita kombinace metody nejmenších čtverců a hrubé síly. Hodnota u je určena pro všechny přípustné kombinace p0, n1 a n2 a následně je z nich vybrána optimální řídicí sekvence pro řízení v otevřené smyčce. Ve skutečnosti je však vždy aplikován pouze prvý krok této řídicí sekvence a v další vzorkovací periodě je celý optimalizační proces zopakován. Tím se řídicí strategie mění na zpětnovazební řízení.

Parametry prediktivního regulátoru, kromě modelu řízené soustavy a omezení jejího vstupu, jsou horizont řízení NC, horizont predikce N1, N2 a koeficient λ. Pouze poslední uvedený parametr λ je určen pro ruční doladění kvality regulace při rutinním uvádění do provozu. V případě použití modelu soustavy ve tvaru přenosu (7.1) nebo (7.2) jsou parametry N1, N2 zvoleny automaticky na základě charakteristických čísel μ,σ2. Regulátor potom může být efektivně laděn „ručně“ pouze seřizováním charakteristických čísel procesu κ, μ, σ2 .

Varování
Při použití bloku PSMPC pro simulaci v systému Matlab/Simulink je třeba zajistit, aby pole sr bylo dostatečně velké, tak aby jím definovaný buffer pojmul interně vygenerovanou přechodovou charakteristiku určenou z FOPDT nebo SOPDT modelu. V opačném případě blok hlásí chybu.

Tento blok propaguje kvalitu signálu. Více informací je uvedeno v sekci 1.4.

Vstup

sp

Požadovaná hodnota (setpoint)

Double (F64)

pv

Řízená veličina

Double (F64)

tv

Veličina pro vysledování

Double (F64)

hv

Hodnota výstupu v manuálním režimu

Double (F64)

MAN

Manuální nebo automatický režim

Bool

off ..

automatický režim

on ...

manuální režim

Parametr

nc

Délka horizontu řízení  5

Long (I32)

np1

Začátek koincidenčního intervalu  1

Long (I32)

np2

Konec koincidenčního intervalu  10

Long (I32)

lambda

Koeficient penalizace změn řízení  0.05

Double (F64)

umax

Horní mez akčního zásahu regulátoru  1.0

Double (F64)

umin

Dolní mez akčního zásahu regulátoru  -1.0

Double (F64)

imtype

Typ modelu řízené soustavy  3

Long (I32)

1 ....

model prvního řádu

2 ....

model druhého řádu

3 ....

přechodová charakteristika

kappa

Statické zesílení  1.0

Double (F64)

mu

Míra zpoždění soustavy  20.0

Double (F64)

sigma

Míra délky odezvy soustavy  10.0

Double (F64)

nmax

Rezervovaná paměť pro pole   10  10000 32

Long (I32)

sr

Diskrétní přechodová charakteristika  [0 0.2642 0.5940 0.8009 0.9084 0.9596 0.9826 0.9927 0.9970 0.9988 0.9995]

Double (F64)

Výstup

mv

Akční zásah regulátoru (manipulated variable)

Double (F64)

dmv

Rychlostní výstup regulátoru (diference)

Double (F64)

de

Regulační odchylka

Double (F64)

SAT

Saturace

Bool

off ..

lineární zákon řízení

on ...

výstup regulátoru je saturován

pve

Odhad řízené veličiny

Double (F64)

iE

Kód chyby

Long (I32)

0 ....

bez chyby

1 ....

nesprávný FOPDT model

2 ....

nesprávný SOPDT model

3 ....

chyba v zadání přechodové charakteristiky

2024 © REX Controls s.r.o., www.rexygen.com