29.01.2015 Views

Diplomová práce - Ústav automatizace a informatiky - Vysoké učení ...

Diplomová práce - Ústav automatizace a informatiky - Vysoké učení ...

Diplomová práce - Ústav automatizace a informatiky - Vysoké učení ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Diplomová práce<br />

Použití frekvenčních charakteristik při<br />

analýze a syntéze regulačních obvodů<br />

Vypracoval:<br />

Vedoucí práce:<br />

Obor:<br />

Specializace:<br />

2006<br />

Miroslav Kij<br />

Ing. Olga Davidová, Ph. D<br />

Inženýrská informatika a <strong>automatizace</strong><br />

Automatizace


Strana 3<br />

Vysoké učení technické v Brně, Fakulta strojního inženýrství<br />

Ústav: Automatizace a informatika<br />

Akademický rok 2005/2006<br />

ZADÁNÍ DIPLOMOVÉ PRÁCE<br />

pro:<br />

Miroslava KIJE<br />

který/která studuje v magisterském studijním programu<br />

obor: M3917 Inženýrská informatika a <strong>automatizace</strong><br />

Ředitel ústavu Vám v souladu se zákonem č.111/1998 o vysokých školách a se Studijním a<br />

zkušebním řádem VUT v Brně určuje následující téma diplomové práce:<br />

POUŽITÍ FREKVENČNÍCH CHARAKTERISTIK PŘI ANALÝZE A SYNTÉZE<br />

REGULAČNÍCH OBVODŮ<br />

a anglickém jazyce:<br />

Analysis and synthesis of control systems using frequency characteristic<br />

Stručná charakteristika problematiky úkolu diplomové práce:<br />

K řešení analýzy a syntézy regulačních obvodů lze použít také metody, které využívají<br />

frekvenční charakteristiky. Cílem diplomové práce je provést jednak ucelený přehled<br />

těchto metod a jednak je aplikovat na příkladech.<br />

Doporučená osnova práce:<br />

1. Popsat frekvenční charakteristiky spojitých a diskrétních systémů<br />

2. Zpracovat metody stability regulačních obvodů využívající frekvenční charakteristiky<br />

3. Zpracovat metody syntézy regulačních obvodů využívající frekvenční charakteristiky<br />

4. Vybrané metody uvést na příkladech


Strana 4<br />

Rozsah grafických prací: ---<br />

Rozsah průvodní zprávy: cca 60 stran<br />

Seznam odborné literatury:<br />

Balátě, J. Automatické řízení. Praha : Nakladatelství BEN-technická literatura, 2003.<br />

664s. ISBN 80-7300-020-2.<br />

Kachaňák, A. Teória automatického riadenia II. Bratislava : Edičné středisko SVŠT v<br />

Bratislavě, 1987. 334s.<br />

Kubík, S., Kotek, Z., Strejc, V. & štacha, J. Teorie automatického řízení I - Lineární a<br />

nelineární systémy. Praha : SNTL Praha, 1982. 528s.<br />

Švarc, I. Teorie automatického řízení I. Brno : Vysoké učení technické v Brně, 1992.<br />

210s. ISBN 80-214-0516-3.<br />

Švarc, I. Teorie automatického řízení II. Brno : Vysoké učení technické v Brně, 1993.<br />

231s. ISBN 80-214-0550-3.<br />

Švarc, I. Automatizace - Automatické řízení. Brno : Vysoké učení technické v Brně, 2005.<br />

262s. ISBN 80-214-2943-7.<br />

Švec, J. Šiška, I. & Vavřín, P. Teorie řízení I - Lineární systémy. Brno : Vysoké učení<br />

technické v Brně, 1982. 208s.<br />

Zítek, P., Hofreiter, M. & Hlava, J. Automatické řízení. Praha : ČVUT Praha, 2001.<br />

148s. ISBN 80-01-02044-4.<br />

Vedoucí diplomové práce: Ing. Olga Davidová, PhD.<br />

Termín odevzdání diplomové práce je stanoven časovým plánem akademického roku<br />

2005/2006.<br />

V Brně dne 21. listopadu 2005<br />

L.S.<br />

Doc. RNDr. Ing. Miloš Šeda, PhD.<br />

ředitel ústavu<br />

Prof. Ing. Josef Vačkář, CSc.<br />

děkan


Strana 5<br />

ANOTACE<br />

Diplomová práce je zaměřena na analýzu a syntézu regulačních obvodů s využitím<br />

frekvenčních charakteristik. Zabývá se zjišťováním stability spojitých a diskrétních<br />

regulačních obvodů a rovněž návrhem optimálních parametrů regulátoru z hlediska kvality<br />

regulace, a to vše za pomocí frekvenčních charakteristik.<br />

ANNOTATION<br />

The diploma work is focused on the analysis and synthesis of control systems using<br />

frequency characteristic. It is dealt with pinpoint stability of the uninterrupted and<br />

unobtrusive control systems and also the proposition of the optimum parameters of the<br />

regulator from the aspect of the regulation quality, and that all by means of the frequency<br />

characteristics.


Strana 7<br />

PODĚKOVÁNÍ<br />

Na úvod mé diplomové práce bych rád poděkoval mé vedoucí práce Ing. Olze<br />

Davidové, Ph.D, za poskytnutí studijních podkladů, absolvování přínosných konzultací s<br />

cennými radami a připomínkami, čímž mě vytvořila optimální podmínky pro její<br />

vypracování.


Strana 9<br />

PROHLÁŠENÍ<br />

Prohlašuji, že jsem tuto diplomovou práci vypracoval samostatně dle rad a pokynů<br />

vedoucího diplomové práce, s použitím odborné literatury a článků uvedených na<br />

internetových stránkách.<br />

23.5.2006 ………………………………………...


Strana 11<br />

OBSAH<br />

Seznam použitých symbolů .........................................................................................13<br />

1 Úvod .....................................................................................................................15<br />

2 Analýza a syntéza................................................................................................17<br />

3 Použití frekvenčních charakteristik u spojitých systémů ...............................21<br />

3.1 Frekvenční popis spojitých systémů ............................................................21<br />

3.1.1 Frekvenční přenos............................................................................21<br />

3.1.2 Frekvenční charakteristika v komplexní rovině ..............................22<br />

3.1.3 Logaritmické frekvenční charakteristiky.........................................24<br />

3.2 Analýza spojitých systémů frekvenčními metodami...................................26<br />

3.2.1 Stabilita regulačních obvodů ...........................................................26<br />

3.2.2 Michajlov-Leonhardovo kritérium ..................................................27<br />

3.2.3 Nyquistovo kritérium.......................................................................29<br />

3.3 Syntéza spojitých systémů frekvenčními metodami....................................36<br />

3.3.1 Ukazatele kvality regulace...............................................................36<br />

3.3.2 Frekvenční metody syntézy.............................................................41<br />

3.3.2.1 Volba zesílení rozpojeného regulačního obvodu....................41<br />

3.3.2.2 Sériové korekční členy ...........................................................42<br />

3.3.2.3 Metoda typizované logaritmické frekvenční charakteristiky .42<br />

4 Použití frekvenčních charakteristik u diskrétních systémů ...........................45<br />

4.1 Frekvenční popis diskrétních systémů.........................................................45<br />

4.1.1 Frekvenční přenos............................................................................45<br />

4.1.2 Frekvenční charakteristika v komplexní rovině ..............................46<br />

4.1.3 Logaritmická frekvenční charakteristika .........................................46<br />

4.2 Analýza diskrétních systému frekvenčními metodami................................47<br />

4.2.1 Stabilita regulačních obvodů ...........................................................47<br />

4.2.2 Michajlov-Leonhardovo kritérium ..................................................49<br />

4.2.3 Nyquistovo kritérium.......................................................................51<br />

4.3 Syntéza diskrétních systémů frekvenčními metodami.................................51<br />

4.3.1 Typizovaná logaritmická amplitudová frekvenční charakteristika .52<br />

4.3.2 Metoda transformovaných frekvenčních charakteristik ..................55<br />

5 Použití vybraných metod na příkladech...........................................................59<br />

5.1 Příklady pro spojité regulační obvody.........................................................59<br />

5.1.1 Michajlov-Leonhardovo kritérium ..................................................59<br />

5.1.2 Nyquistovo kritérium.......................................................................61<br />

5.2 Příklady pro diskrétní regulační obvody......................................................62<br />

5.2.1 Michajlov-Leonhardovo kritérium ..................................................62<br />

5.2.2 Nyquistovo kritérium.......................................................................64<br />

5.2.3 Metoda transformovaných frekvenčních charakteristik ..................66<br />

6 Závěr ....................................................................................................................69<br />

Seznam použité literatury............................................................................................71


Strana 13<br />

SEZNAM POUŽITÝCH SYMBOLŮ<br />

a 0 činitel autoregulace<br />

a i koeficient jmenovatele přenosu<br />

A w modul (amplituda) kmitočtového přenosu řízení uzavřeného regulačního obvodu<br />

A r rezonanční převýšení<br />

A amplituda<br />

e regulační odchylka<br />

e T (t) trvalá regulační odchylka<br />

G 0 přenos rozpojeného regulačního obvodu<br />

G R přenos regulátoru<br />

G S přenos regulované soustavy<br />

k 0 zesílení rozpojeného regulačního obvodu<br />

k 0opt optimální zesílení rozpojeného regulačního obvodu<br />

n stupeň charakteristického polynomu zavřeného regulačního obvodu<br />

r 0 proporcionální složka regulátoru<br />

r 1 derivační složka regulátoru<br />

r -1 integrační složka regulátoru<br />

t r doba regulace<br />

t max čas<br />

T perioda kmitů přechodové charakteristiky<br />

T 95 doba, za kterou přechodová charakteristika regulované soustavy dosáhne 95%<br />

ustálené hodnoty<br />

T d dopravní zpoždění<br />

T D derivační časová konstanta<br />

T I integrační časová konstanta<br />

v poruchová veličina<br />

w řídicí veličina<br />

y regulovaná veličina<br />

y hom homogenní rovnice<br />

y max maximální překmit<br />

y part partikulární integrálγ fázová bezpečnost<br />

δ amplitudová bezpečnost<br />

ϕ fáze frekvenčního přenosu<br />

ω úhlová frekvence používaná pro spojité systémy<br />

ω d šířka pásma používaná pro spojité systémy<br />

ω Td fázové zpoždění<br />

ω ř úhlová frekvence řezu používaná pro spojité systémy<br />

ω r rezonanční frekvence používaná pro spojité systémy<br />

ω h šířka pásma propustnosti používaná pro spojité systémy<br />

Ω úhlová frekvence používaná pro diskrétní systémy<br />

frekvence řezu používaná pro diskrétní systémy<br />

Ω ř


Strana 15<br />

1 ÚVOD<br />

Bavíme-li se o regulačních obvodech, pak nedílnou součástí se stává posuzování<br />

(analýza) jejich stability, tj. schopnosti obvodu ustálit se a vykazovat tak správnou funkci.<br />

Pro zjišťování této stability využíváme několika metod. V této práci se zaměřím na metody<br />

využívající frekvenční charakteristiky. Aby byla jasná podstata, je nejprve potřeba rozebrat<br />

frekvenční charakteristiky pro spojité i diskrétní systémy a teprve poté se zaměřit na<br />

jednotlivé metody. Rovněž zde uvedu, jakým způsobem se nastavují parametry regulátoru<br />

(syntéza).<br />

Pro analýzu regulačního obvodu je vypracována řada metod, díky kterým lze<br />

snadno určit stabilitu regulačního obvodu, aniž bychom museli provádět složité výpočty,<br />

na které v dnešní době využíváme téměř výhradě výpočetní techniku.<br />

Pro návrh regulačního obvodu nejsou prozatím vypracovány exaktní metody, které<br />

by jednoznačně vedly k cíli. Stále zde tedy hraje důležitou roli intuice plynoucí ze<br />

zkušenosti a citu navrhovatele. V drtivé většině spočívá metoda návrhu v přímé extrapolaci<br />

nebo interpolaci z již realizovaných a vyzkoušených řešení.<br />

V oblasti syntézy je však zpracována celá řada dobře využitelných metod.<br />

Rozhodující je, aby bylo dosaženo převedení požadavků na regulační obvod, které<br />

formuluje provozovatel, konstruktér a projektant regulovaného projektu na matematickou<br />

formulaci požadavků, kritérií a cílů, vhodných pro další zpracování. Dosud v praxi<br />

převládají jednoduché metody syntézy bez velkých matematických nároků.<br />

V poslední době se rozvíjejí metody syntézy regulačních obvodů, u kterých<br />

určování vlastností řízeného systému a jeho řízení probíhá podle jediného algoritmu, nebo<br />

se vyvíjejí algoritmy řízení, které nevyžadují výchozí přesnou identifikaci řízeného<br />

systému a tyto způsoby řízení označujeme jako adaptivní řízení.<br />

Náplní této diplomové práce je využití frekvenčních charakteristik při analýze a<br />

syntéze regulačních obvodů. Vybrané metody analýzy a syntézy pak budou aplikovány na<br />

příkladech, kde jednak bude zkontrolována stabilita regulačního obvodu a také budou<br />

navrženy parametry jednotlivých regulátorů (P, I, PI, PD, PID). Výsledky budou<br />

prezentovány formou tabulek a grafů a následně vysvětlujícím textem se závěrem o<br />

stabilitě či nestabilitě regulačního obvodu a s výsledkem návrhu parametrů regulátoru.


Strana 17<br />

2 ANALÝZA A SYNTÉZA<br />

I přesto, že je v úvodu (kap. 1) naznačen popis pojmů analýza a syntéza, zdá se být<br />

na místě uvést podrobnější popis, který nás hlouběji vtáhne do problematiky a pomůže nám<br />

pochopit níže popsané skutečnosti.<br />

Analýza a syntéza patří mezi základní a nejčastěji užívané vědecké metody.<br />

Původní význam řeckého slova analýza znamenalo rozložení nějakého komplexu na části a<br />

syntéza měla význam spojení rozmanitostí k jednotě v celku. Z metodologického hlediska<br />

jsou tato slova používána ve smyslu metod k získávání nových poznatků, nebo ve smyslu<br />

metody výkladu poznatků. Syntéza je proti analýze proces opačný, nebo doplňující. Jde o<br />

sjednocování, složení nějakého předmětu, jevu či procesu z jeho základních prvků ať již<br />

myšlenkově, či fakticky v nějaký celek. Toto sjednocování nemusí být jen u jednotlivých<br />

částí, které byly předtím vyděleny analýzou. Syntéza má však jako metodologický princip<br />

analýzu vždy doplňovat. Tím syntéza umožňuje poznání předmětu v jeho úplnosti. Pomocí<br />

syntézy nalézáme vztahy nějakého jevu k jiným jevům, zařazujeme jev, nebo proces do<br />

většího celku a objasňujeme vztahy a mechanismus funkcí u tohoto jevu. Syntézou lze<br />

rozumět také takový proces, při němž hledáme spojováním části v celek takovou strukturu,<br />

která by měla námi předem požadované chování. V tomto případě syntéza není pouhou<br />

skladbou jednotlivých jevů čí procesů, ale je to zároveň kreace nových celků, případně<br />

jejich proměna. Syntéza tedy může být hledáním nejvhodnější varianty dosahované<br />

kombinací jednotlivých prvků a jejich vlastností. Jednoduše řečeno pojmem "syntéza"<br />

rozumíme stanovení takové struktury a parametrů regulačního obvodu, aby byly splněny<br />

požadavky, které klademe na regulační pochod.<br />

Syntéza je první etapou návrhu regulačního obvodu. Dalšími etapami návrhu<br />

regulačního obvodu je volba jednotlivých členů regulačního obvodu. Při návrhu<br />

regulačního obvodu vycházíme z provozních podmínek, kladených na regulační obvod. Do<br />

provozních podmínek zahrnujeme např. požadavky na rozměry a hmotnost zařízení<br />

(zvláště u létajících objektů - letadel, družic apod.), dále požadavky na pracovní prostředí<br />

(vlhkost, agresivita, nevýbušnost), režim provozu (dlouhodobý, krátkodobý), požadavky na<br />

typizaci s ohledem na jednoduchost údržby celého zařízení a také požadavky, zda<br />

regulační obvod má být složen z přístrojů elektrických, pneumatických nebo<br />

hydraulických.<br />

Pro tyto fáze návrhu nejsou zatím vypracovány žádné exaktní metody, které by<br />

jednoznačně vedly k cíli. Zde stále ještě hraje důležitou roli intuice, zkušenost a<br />

konstruktérský cit navrhovatele. V praxi metoda návrhu spočívá především v tom, že<br />

přímo extrapolujeme z již realizovaných a vyzkoušených řešení.<br />

V otázkách syntézy regulačního obvodu se však můžeme opřít o řadu dobře<br />

vypracovaných metod, kterými se zde budeme zabývat.<br />

Při syntéze regulačního obvodu je velmi důležité převést požadavky na regulační<br />

obvod, které formuluje provozovatel regulovaného objektu (technolog u technologického<br />

procesu, pilot letounu apod.), na matematickou formulaci požadavků, kritérií a cílů,<br />

vhodnou pro další zpracování. Dosud v praxi převládají jednoduché metody syntézy bez<br />

velkých matematických nároků.<br />

Výchozí předpoklady pro syntézu mohou být zadány různě:<br />

1. Můžeme volit libovolně strukturu i parametry regulačního obvodu a jsme<br />

omezeni jen splněním podmínek fyzikální realizovatelnosti.


Strana 18<br />

Analýza a syntéza<br />

2. Je zadána část struktury i část parametrů regulačního obvodu<br />

3. Je plně zadána struktura a jsou zadány některé parametry regulačního obvodu<br />

S případem 1 se v praxi setkáváme zřídka. Vyskytuje se např. při syntéze filtrů.<br />

Prakticky všechny úlohy syntézy regulačního obvodu lze zahrnout pod body 2 a 3.<br />

Pod bod 3 patří velké množství průmyslových regulací, u kterých lze regulační<br />

obvod rozdělit na regulátor a regulovaný systém. Regulátory pro průmyslové regulace se<br />

vyrábějí jako univerzální regulátory s pevnou strukturou - obvykle spojitý regulátor PID.<br />

Úloha syntézy se zde redukuje pouze na určení nastavitelných parametrů regulátoru.<br />

Vhodnost správné volby typu regulátoru ověříme po provedené syntéze regulačního<br />

obvodu jeho simulací na matematickém modelu regulačního obvodu a později provozními<br />

zkouškami na regulovaném objektu přímo v provozu.<br />

Pod bod 2 patří regulační obvody, u kterých neprovádíme rozdělení na regulátor a<br />

řízený systém. Jsou to např. servomechanismy, což jsou regulační obvody sloužící k vlečné<br />

regulaci polohy nebo jejích derivací. U těchto regulačních obvodů navrhujeme jak jejich<br />

strukturu, tak i jejich parametry.<br />

ϕ 1<br />

+<br />

ϕ<br />

ovladač<br />

servomotor<br />

ϕ 2<br />

−<br />

ϕ 2<br />

Obr. 2.1 - Polohový servomechanismus<br />

Servomechanismus je tedy regulační obvod, který je tvořen tzv. ovládačem a<br />

servomotorem s převodovkou (obr. 2.1). Ovladač působící na servomotor je tvořen<br />

výkonovým zesilovačem, předzesilovačem s korekcemi a měřícím členem.<br />

Při návrhu servomechanismu navrhujeme všechny jeho členy. Servomotor a jeho<br />

převodovku navrhujeme podle výkonu požadovaného na výstupní hřídeli. Typem<br />

servomotoru je určen výkonový zesilovač, kterým u elektrických servomotorů bývá<br />

obvykle tyristorový měnič. Volba měřícího členu je ovlivněna požadovanou přesností<br />

servomechanismu. Předzesilovač s jeho korekčními členy navrhujeme podle požadavků na<br />

dynamické vlastnosti servomechanismu.<br />

Při syntéze servomechanismu vycházíme obvykle z daných vlastností výkonového<br />

zesilovače, servomotoru s převodovkou a měřícího členu a určujeme vlastnosti<br />

předzesilovače s korekcemi.<br />

Při syntéze regulačního obvodu potřebujeme znát<br />

1. vlastnosti regulovaného objektu,<br />

2. předpokládaný průběh řídicí veličiny,<br />

3. předpokládané průběhy poruchových veličin a místa jejich vstupu do řízeného<br />

systému,<br />

4. omezení akčních veličin,<br />

5. požadavky na kvalitu regulace.<br />

Vlastnosti regulovaného objektu určujeme buď analýzou objektu, nebo rozborem<br />

experimentálně získaných průběhů veličin v objektu. Obě tyto metody identifikace vedou<br />

na sestavení matematického modulu - řízeného systému. Přitom se ukazuje, že některá


Analýza a syntéza Strana 19<br />

zjednodušení, která provádíme při analýze vlastností objektu, se v uzavřeném regulačním<br />

obvodu nepříznivě projeví tím, že model a reálný objekt se v uzavřeném regulačním<br />

obvodu chovají odlišně. Proto jsou výhodné ty metody identifikace, při kterých je<br />

identifikovaný objekt zapojen přímo v regulačním obvodu (např. identifikace pomocí<br />

adaptivního modelu).<br />

Uvedené potíže vedly k tomu, že se v poslední době rozvíjejí syntézy regulačního<br />

obvodu, u kterých určování vlastností objektu a jeho řízení probíhá současně podle<br />

jediného algoritmu, nebo se vyvíjejí algoritmy řízení, které nevyžadují přesnou identifikaci<br />

objektu. Tyto způsoby řízení označujeme jako adaptivní řízení.<br />

Rozborem fyzikálních jevů, které nastanou během regulačního pochodu, určíme<br />

podstatu regulované veličiny i případných poruchových veličin a odhadneme jejich průběh<br />

v čase. Fyzikální podstata regulované veličiny bude ovlivňovat volbu čidla, které převádí<br />

regulovanou veličinu na signál vhodný k dalšímu zpracování. Časová funkce, podle které<br />

se bude měnit regulovaná veličina, bude především ovlivňovat požadavky na kvalitu a<br />

přesnost regulace, tj. požadavky na dynamiku regulačního pochodu.<br />

Někdy požadujeme, aby regulovaná veličina přesně sledovala řídicí veličinu. Např.<br />

u polohových servomechanismů (obr. 2.1) požadujeme, aby výstupní poloha ϕ 2 hřídele<br />

servomechanismu co nejvěrněji sledovala průběh žádané hodnoty ϕ 1 a vliv poruchových<br />

veličin (zatěžovacích momentů) často ani neuvažujeme. Naopak při regulaci teploty,<br />

napětí, síťového kmitočtu apod. je často žádaná hodnota trvale konstantní a úkolem<br />

regulace je kompenzovat poruchy vstupující do regulovaného objektu.<br />

Ve skutečnosti můžou mít žádané hodnoty regulovaných veličin i poruchy<br />

vstupující do regulovaného objektu zcela obecný průběh. Pro zjednodušení výpočtu<br />

uvažujeme jako vstupní veličiny tzv. typizované funkce, jejichž matematické vyjádření je<br />

snadné a z odezvy regulačního obvodu na tyto funkce můžeme soudit na přesnost a kvalitu<br />

regulace. Nejčastěji používané typizované funkce jsou jednotkový skok, Diracův impuls,<br />

skok rychlosti a zrychlení vstupního průběhu a harmonický průběh.<br />

Akční veličina je výstupní veličinou regulátoru a zároveň vstupní veličinou<br />

řízeného systému. Velikost akční veličiny je vždy omezena. Omezení akční veličiny je<br />

způsobeno jednak tím, že výstupní signál z regulátoru nemůže nabývat libovolně velké<br />

hodnoty, a jednak tím, že vstupní signál do regulovaného objektu se může měnit pouze v<br />

dovoleném rozsahu, který je dán fyzikální podstatou objektu i ekonomickými omezeními.<br />

Omezení jsou v podstatě dvojího druhu. Pro názornost je označíme jako omezení<br />

typu "skála" a "propast". Omezení typu "skála" je omezení na dorazech; toto omezení<br />

nemůžeme nikdy přesáhnout. Např. regulační ventil má krajní polohy, kdy je plně otevřen<br />

nebo uzavřen. Omezení typu "propast" je omezení, jehož překročení může způsobit<br />

poruchu zařízení. Proto musíme zajistit, abychom hodnotu tohoto omezení nikdy<br />

nepřekročili. Např. napětí nebo proud kotvy motoru nesmí nikdy přesáhnout velikost<br />

plynoucí z izolační pevnosti vinutí, popř. jeho přípustného oteplení.,<br />

Respektování omezení veličin značně komplikuje úlohu syntézy. Proto nejčastěji<br />

provádíme syntézu regulačního obvodu bez respektování omezení a potom kontrolujeme,<br />

zda při dané velikosti řídicích a poruchových veličin nastane omezení signálů a jak se<br />

omezení projeví na dynamických vlastnostech celého regulačního obvodu. [Kubík, 1982]<br />

Jak již jsme uvedli výše, syntéza regulačního obvodu ve frekvenční oblasti vychází<br />

ze struktury regulačního obvodu, kde není možné respektovat místo a tvar vstupující<br />

poruchové veličiny, ale je možno dosáhnout pouze požadovaných frekvenčních vlastností<br />

obvodu.


Strana 21<br />

3 POUŽITÍ FREKVENČNÍCH CHARAKTERISTIK U<br />

SPOJITÝCH SYSTÉMŮ<br />

K řízení reálných objektů se rozvíjí metody automatické regulace. Aby tyto metody<br />

byly použitelné pro širokou třídu reálných objektů, jsou vytvořeny abstraktní modely<br />

reálných objektů, které se nazývají systémy. Abstrakcí velmi široké třídy reálných modelů<br />

vznikly spojité systémy, u nichž jsou všechny veličiny funkcemi času t. [Kubík, 1982]<br />

3.1 Frekvenční popis spojitých systémů<br />

V této části třetí kapitoly se budu zabývat frekvenčním popisem spojitých<br />

regulačních obvodů. Vysvětlím zde, které pojmy a výpočtové vztahy využiji při analýze a<br />

syntéze regulačních obvodů. Důležité rovněž bude ukázat si postup při konstrukci<br />

frekvenčních charakteristik.<br />

3.1.1 Frekvenční přenos<br />

Frekvenční přenos se získá tak, že je na vstup systému přiveden harmonický signál.<br />

Typickým harmonickým signálem je sinusový průběh<br />

u t)<br />

= u sin ωt<br />

(<br />

0<br />

(3.1)<br />

amplituda vstupního<br />

signálu<br />

úhlová<br />

frekvence<br />

Na výstupu systému se dostane podle obr. 3.1 (po odeznění přechodového jevu)<br />

opět sinusový signál ovšem s jinou amplitudou, stejnou úhlovou frekvencí a fázově proti<br />

vstupnímu signálu posunutý<br />

u(t)<br />

y(t)<br />

y<br />

( t) = y sin( ω t + ϕ )<br />

0<br />

(3.2)<br />

T<br />

u 0<br />

t<br />

ϕ<br />

T<br />

y 0<br />

t<br />

Lépe se ale jeví vyjádřit vstupní i<br />

výstupní funkci v komplexním tvaru<br />

jωt<br />

( t) = u e ;<br />

0<br />

u(t)<br />

S<br />

y(t)<br />

j ( ω t +ϕ )<br />

( )<br />

t = y0e<br />

(3.3)<br />

Obr. 3.1


Strana 22<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

To jsou v komplexní rovině vektory, které se otáčí úhlovou rychlostí ω. Poměr<br />

těchto vektorů nám definuje frekvenční přenos<br />

j(<br />

ωt+<br />

ϕ )<br />

( t)<br />

y0e<br />

G( jω)<br />

= = =<br />

jωt<br />

( t)<br />

u e<br />

0<br />

y<br />

u<br />

0<br />

0<br />

e<br />

jϕ<br />

(3.4)<br />

kde<br />

y<br />

u<br />

0<br />

0<br />

je poměr amplitud a ϕ je fázové posunutí.<br />

<br />

Základem všeho je důležitý vzorec pro výpočet přenosu z diferenciální rovnice<br />

bms<br />

G(<br />

s)<br />

=<br />

a s<br />

n<br />

m<br />

n<br />

+ ... + b s + b<br />

1<br />

+ ... + a s + a<br />

1<br />

0<br />

0<br />

(3.5)<br />

Pro výpočet frekvenčního přenosu z koeficientů diferenciální rovnice lze odvodit<br />

následující vztah<br />

b<br />

G(<br />

jω)<br />

=<br />

a<br />

m<br />

n<br />

( jω)<br />

( jω)<br />

m<br />

n<br />

+ ... + b1<br />

jω<br />

+ b<br />

+ ... + a jω<br />

+ a<br />

1<br />

0<br />

0<br />

(3.6)<br />

Vztah je formálně stejný jako vztah (3.5) pro přenos G(s), pouze místo komplexní<br />

proměnné s v něm figuruje výraz jω. Tím je zároveň dána relace mezi přenosem a<br />

frekvenčním přenosem, která spočívá ve formální záměně s za jω eventuálně naopak<br />

G ( jω)<br />

= G(<br />

s)<br />

G ( s)<br />

= G(<br />

jω)<br />

s=<br />

jω;<br />

jω<br />

= s;<br />

(3.7)<br />

Zavedení frekvenčního přenosu má velký praktický význam pro řešení regulačních<br />

problémů. Frekvenční přenos je základem pro používání frekvenčních metod. Znázornění<br />

frekvenčního přenosu ve tvaru frekvenčních charakteristik umožní řešit otázky stability<br />

regulačních obvodů, kvalitu regulace i syntézu regulačních obvodů. Také je možno<br />

používat experimentálně zjištěné a naměřené frekvenční charakteristiky.<br />

3.1.2 Frekvenční charakteristika v komplexní rovině<br />

Použití frekvenčních charakteristik v komplexní rovině se dá považovat za<br />

nejčastější způsob při určování stability regulačních obvodů.<br />

Frekvenční charakteristika je grafické vyjádření frekvenčního přenosu G(jω)<br />

v komplexní rovině, když se za úhlovou frekvenci ω dosazují hodnoty 0 až ∞. [Švarc,<br />

2002]


Použití frekvenčních charakteristik u spojitých systémů Strana 23<br />

Při praktickém sestrojování frekvenční charakteristiky si frekvenční přenos G(jω)<br />

ještě v obecném tvaru (před dosazením hodnot ω) upravím na složkový tvar komplexního<br />

čísla (rozšířením zlomku číslem komplexně sdruženým ke jmenovateli).<br />

[ G(<br />

jω)<br />

] j Im[ G(<br />

j )]<br />

G ( jω)<br />

= Re + ω<br />

Při pohledu na obr. 3.2a je možné vidět sestrojení frekvenční charakteristiky ( jω)<br />

(3.8)<br />

G podle<br />

(3.8), kde za ω byla dosazena zvolená čísla.<br />

Druhým způsobem, jak lze sestrojit frekvenční charakteristiku v komplexní rovině<br />

je z exponenciálního tvaru komplexního čísla (obr. 3.2b). Číslo a + jb se vyjádří ve<br />

složkovém, goniometrickém nebo exponenciálním tvaru. Po úpravách díky Eulerova<br />

vztahu [Švarc, 2002] pro převod goniometrického tvaru na exponenciální se dostaneme<br />

vztah<br />

G(<br />

jω)<br />

=<br />

A(<br />

ω)<br />

⋅ e<br />

jϕ (ω )<br />

.<br />

(3.9)<br />

Odvozování Eulerových vztahů zde uvádět nebudu. Uvedu zde pouze dva vztahy.<br />

První (3.10) v podstatě vyplývá z Pythagorovy věty a slouží k určení velikosti<br />

amplitudy A . Druhým vzorcem (3.11) zjistím fázi ϕ . Oba údaje se uvádí do tabulky při<br />

sestrojování frekvenční charakteristiky v komplexní rovině.<br />

2<br />

A = a +<br />

b<br />

b<br />

ϕ = arctg<br />

a<br />

2<br />

(3.10)<br />

(3.11)<br />

a) Im<br />

b)<br />

Im<br />

ω = ∞<br />

ω = 20<br />

ω = 5<br />

ω = 2<br />

a = Acosα<br />

b = Asinα<br />

ω =1<br />

Re<br />

ω = 0<br />

ω = 20<br />

ω = 0,5<br />

ω = 5<br />

ω = ∞<br />

G ( jω)<br />

G( jω)<br />

ω = 2<br />

A<br />

ϕ<br />

ω =1<br />

Re<br />

ω = 0<br />

ω = 0,5<br />

a + jb<br />

Obr. 3.2 - Frekvenční charakteristika v komplexní rovině<br />

Pro názornost takovou tabulku uvedu a k ní pak sestrojím frekvenční<br />

charakteristiku v komplexní rovině. Ke zvoleným hodnotám ω na kalkulačce dopočítám<br />

hodnoty Re a Im a výsledky doplním do tabulky 3.1.


Strana 24<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

Tabulka 3.1<br />

ω Re(ω) Im(ω) A(ω) ϕ(ω)<br />

0 1,5 0 1,5 0°<br />

0,1 1.39 -0,43 1,46 -17°01'<br />

0,2 1,14 -0,74 1,37 -33°07'<br />

0,3 0,83 -0,91 1,23 -47°40'<br />

0,4 0,54 -0,95 1,09 -60°28'<br />

0,5 0,3 -0,90 0,95 -71°34'<br />

0,7 0,01 -0,71 0,71 -89°27'<br />

0,8 -0,07 -0,62 0,62 -96°49'<br />

1,0 -0,15 -0,45 0,47 -108°26'<br />

1,5 -0,16 -0,21 0,26 -127°52'<br />

2 -0,12 -0,11 0,16 -139°24'<br />

10 -0,01 -0,01 0,1 -171°33'<br />

ω 0 0 0 -180°<br />

Podle této tabulky pak frekvenční charakteristiku zkonstruuji. Hodnoty ω jsem<br />

dosazoval do rovnice, kterou zde neuvádím, ale pro názornost je dobré si ukázat, jak<br />

taková tabulka vypadá a hlavně jak vypadá frekvenční charakteristika (obr. 3.3) sestrojená<br />

z těchto hodnot.<br />

Im<br />

2<br />

1,5<br />

1<br />

0,8<br />

ω = ∞<br />

G( jω)<br />

Re<br />

ω = 0<br />

0,1<br />

0,2<br />

0,3<br />

0,5<br />

0,4<br />

Obr. 3.3 - Frekvenční charakteristika v komplexní rovině sestrojená podle tabulky 3.1<br />

3.1.3 Logaritmické frekvenční charakteristiky<br />

Frekvenční charakteristiku v komplexní rovině se také někdy nazývá amplitudofázovou<br />

frekvenční charakteristikou. Z jednoho bodu této charakteristiky se pro danou<br />

frekvenci odečte amplitudu A i fázi ϕ frekvenčního přenosu (3.9)


Použití frekvenčních charakteristik u spojitých systémů Strana 25<br />

G(<br />

jω)<br />

=<br />

A(<br />

ω)<br />

⋅ e<br />

jϕ<br />

( ω )<br />

Tím pádem je možno tuto charakteristiku rozdělit na dvě charakteristiky, amplitudovou<br />

A=A (ω) a fázovou ϕ=ϕ (ω), jak je ukázáno na obr. 3.4.<br />

Vhledem k úzkému frekvenčnímu pásmu se lineární souřadnice, uvedené na obr.<br />

3.4, příliš nepoužívají. Pokud bych chtěl toto pásmo rozšířit, pak nejdůležitější část<br />

charakteristiky s podstatnou změnou amplitudy by byla nahuštěna v úzkém rozsahu<br />

frekvencí [Švarc, 2002]. Také pracnost konstrukce těchto charakteristik je poměrně velká<br />

(stejně jako předchozích charakteristik v komplexní rovině).<br />

amplitudo-fázová<br />

frekvenční charakteristika<br />

v komplexní rovině<br />

Im<br />

Re<br />

A<br />

ϕ<br />

ω = 0,15<br />

A<br />

[–]<br />

ϕ<br />

[°]<br />

amplitudová v lineárních<br />

souřadnicích A<br />

20<br />

[dB]<br />

10<br />

0 0,1 0,2 0,3<br />

-20°<br />

-40°<br />

-60°<br />

ω[1/s]<br />

fázová v lineárních<br />

souřadnicích<br />

ϕ<br />

[°]<br />

0<br />

-45°<br />

-90°<br />

-135°<br />

-180°<br />

amplitudová v<br />

logaritmických<br />

souřadnicích<br />

0,1 1 10 100<br />

ω[1/s]<br />

fázová<br />

v logaritmických<br />

souřadnicích<br />

Obr. 3.4 - Frekvenční charakteristiky v amplitudových a fázových souřadnicích<br />

Proto se začalo upřednostňovat použití těchto charakteristik v logaritmických<br />

souřadnicích. Na vodorovné ose obou charakteristik, amplitudové i fázové, je vynesena<br />

frekvence v logaritmickém měřítku. Tím se dosáhne velkého rozmezí frekvencí ω.<br />

U amplitudové frekvenční charakteristiky v logaritmických souřadnicích se na<br />

svislou osu vynáší amplituda frekvenčního přenosu G(jω) a to v jednotkách decibel [dB.]<br />

A(<br />

ω)<br />

=<br />

G(<br />

jω)<br />

=<br />

y<br />

u<br />

0<br />

0<br />

e<br />

jϕ<br />

=<br />

y<br />

u<br />

0<br />

0<br />

(3.12)<br />

Decibel je dvacetinásobek dekadického logaritmu zesílení. Amplituda A je podíl amplitud<br />

výstupního a vstupního sinusového signálu y 0 /u 0 tedy zesílení označené A[dB] vyjádřené<br />

vzorcem<br />

y<br />

(3.13)<br />

0<br />

A [ dB]<br />

= 20log A[<br />

−]<br />

= 20log<br />

u<br />

U fázových frekvenčních logaritmických charakteristik je fáze vynášena na svislou osu v<br />

lineárním měřítku (ve stupních nebo radiánech).<br />

0


Strana 26<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

3.2 Analýza spojitých systémů frekvenčními metodami<br />

V této kapitole budou analyzovány vlastnosti spojitých řízených systémů. Budou<br />

analyzovány různé popisy řízených systémů a uvedu zde některé metody určení vlastností<br />

a stability systému. Nejprve zde vysvětlím pojem stabilita regulačního obvodu a poté se<br />

zaměřím na jednotlivá kritéria stability.<br />

3.2.1 Stabilita regulačních obvodů<br />

Stabilita je jedním ze základních požadavků, které se kladou na regulační obvod.<br />

Regulační obvod je stabilní, jestliže po vychýlení regulačního obvodu z rovnovážného<br />

stavu a odeznění vnějších sil, které tuto odchylku způsobily, se regulační obvod během<br />

času znovu vrátí do původního rovnovážného stavu. Jinak řečeno je stabilita vlastnost<br />

regulačního obvodu udržet se v okolí rovnovážného stavu nebo se do něj vrátit po<br />

odeznění vnějších působících sil.<br />

Z hlediska stability se regulační obvod rozlišuje na stabilní, na mezi stability a<br />

nestabilní (obr. 3.5). Regulační obvod na mezi stability se obecně považuje za stabilní a<br />

vždy se vyžaduje, aby regulační obvod byl za všech okolností stabilní. Zatímco parametry<br />

a dynamické vlastnosti regulované soustavy jsou dány konstrukcí soustavy,<br />

technologickým procesem apod. a nemůžeme je tudíž měnit, můžeme měnit dynamické<br />

vlastnosti regulátoru nastavováním volitelných parametrů regulátoru. Tím lze dosáhnout<br />

stability (a dalších vlastností) regulačního obvodu.<br />

y hom<br />

( t)<br />

y hom<br />

( t)<br />

( t)<br />

a) b) c)<br />

y hom<br />

t<br />

t<br />

t<br />

stabilní obvod<br />

obvod na hranici<br />

stability<br />

nestabilní obvod<br />

Obr. 3.5 - Určení stability regulačního obvodu<br />

Jak je uvedeno v [Balátě, 2004], nutnou a postačující podmínkou pro stabilitu<br />

uzavřeného lineárního regulačního obvodu je, aby všechny kořeny charakteristické<br />

rovnice obvodu měly zápornou reálnou část. Z této definice jasně vyplývá, že aby byl<br />

regulační uzavřený lineární obvod stabilní, musí všechny kořeny charakteristické rovnice<br />

ležet v levé polorovině komplexní roviny "s" (obr. 3.6).<br />

Vzhledem ke skutečnosti, že samotné určení stability je pracné i s použitím<br />

výpočetní techniky, protože je zapotřebí vyčíslení kořenů charakteristické rovnice vyššího


Použití frekvenčních charakteristik u spojitých systémů Strana 27<br />

než druhého stupně, byla sestavena<br />

matematická kritéria, tzv. kritéria stability,<br />

která umožňují z charakteristické rovnice<br />

určit, zdali jsou její kořeny se zápornou<br />

reálnou částí nebo ne, a tím stabilitu<br />

obvodu, aniž by se musela daná rovnice<br />

řešit.<br />

Kritéria stability lze rozdělit na<br />

algebraická a frekvenční. Vzhledem k<br />

tématu mé práce se budu věnovat jen<br />

frekvenčním, mezi něž patří dvě<br />

nejpoužívanější, a to Michajlov-<br />

Leonhardovo kritérium a Nyquistovo<br />

kritérium.<br />

nestabilní<br />

oblast<br />

Im<br />

Obr. 3.6<br />

stabilní<br />

oblast<br />

Re<br />

hranice<br />

stability<br />

3.2.2 Michajlov-Leonhardovo kritérium<br />

Jedná se o frekvenční kritérium, které vychází z charakteristické rovnice<br />

uzavřeného regulačního obvodu<br />

a s<br />

n ... 0<br />

(3.14)<br />

n<br />

+ + a1 s + a0<br />

=<br />

Pro označení Michajlov-Leonhardovy křivky se používají symboly N nebo H. V<br />

této práci budu používat symbol H.<br />

Kritérium hodnotí stabilitu podle křivky, kterou opíše koncový bod<br />

charakteristického vektoru H(jω) v komplexní rovině při změně frekvence ω od 0 do ω.<br />

Vektor H(jω) vznikne z charakteristické funkce dosazením s = jω<br />

H<br />

n<br />

( jω)<br />

= an ( jω) + ... + a1( jω) + a0<br />

(3.15)<br />

Tato křivka se nazývá křivkou H(jω) nebo také Michajlovov-Leonhardovou<br />

křivkou. (A.V. Michajlov, ruský matematik, jeho práce uveřejněna v roce 1938; A.<br />

Leonhard, německý technik, práce uveřejněna v roce 1943).<br />

Křivku H(jω) není nutné vždy kreslit celou, postačí jen vypočítat polohu jejich<br />

průsečíků se souřadnými osami. V tom případě se reálná a imaginární část výrazu H(jω)<br />

položí rovna nule a z toho se vypočítají frekvence zmíněných průsečíků. Z frekvencí se<br />

pak určí jejich poloha a z polohy snadno určíme průběh celé charakteristiky. [Švarc, 2002]<br />

Definice Michajlov-Leonhardova kritéria stability:<br />

Uzavřený regulační obvod je stabilní tehdy a jen tehdy, když Michajlova<br />

charakteristika H(jω) začíná na kladné reálné poloose [H(0) = a 0 > 0] a při změně<br />

úhlového kmitočtu ω od 0 do ∞ postupně v kladném smyslu (tj. proti směru pohybu<br />

hodinových ručiček) projde n kvadranty (ke n je stupeň charakteristického polynomu<br />

zavřeného regulačního obvodu).


Strana 28<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

Začíná-li Michajlovova charakteristika H(jω) v počátku souřadnic, pak<br />

charakteristický mnohočlen uzavřeného regulačního obvodu H(s) má nejméně jeden<br />

nulový kořen a regulační obvod je na nekmitavé mezi stability. Na obr. 3.7 je uveden<br />

případ, kdy n = 3 . Zde je možné pozorovat průběh křivky H ( jω)<br />

pro stabilní či nestabilní<br />

obvod nebo pro obvod na hranici stability.<br />

Im<br />

a) b)<br />

Im<br />

ω = ∞<br />

H ( jω)<br />

ω = 0<br />

Re<br />

ω = ∞<br />

H ( jω)<br />

ω = 0<br />

Re<br />

c)<br />

Im<br />

d)<br />

ω = ∞<br />

H ( jω)<br />

Im<br />

H ( jω)<br />

ω = 0<br />

Re<br />

ω = 0<br />

Re<br />

ω = ∞<br />

e)<br />

ω = ∞<br />

H ( jω)<br />

Im<br />

ω = 0<br />

Re<br />

Obr. 3.7 - a) stabilní obvod; b) obvod na hranici stability; c), d), e) nestabilní obvod


Použití frekvenčních charakteristik u spojitých systémů Strana 29<br />

3.2.3 Nyquistovo kritérium<br />

Budu-li rozebírat toto kritériu, je na úvod potřeba říci, že mezi jeho výhody patří to,<br />

že není třeba znát přenos nebo diferenciální rovnici rozpojeného regulačního obvodu, ale<br />

stačí vycházet z experimentálně zjištěné frekvenční charakteristiky rozpojeného<br />

regulačního obvodu. Na rozdíl od algebraických kritérií, kde se zkoumá regulační obvod<br />

pouze z pohledu stability, lze za pomoci Nyquistova kritéria určit také kvalitu regulačního<br />

obvodu, jednoduše řečeno jak moc je obvod stabilní.<br />

Na základě frekvenční charakteristiky rozpojeného regulačního obvodu umožňuje<br />

Nyquistovo kritérium stability ověřovat stabilitu uzavřeného regulačního obvodu.<br />

Frekvenční charakteristika může být k dispozici i v podobě grafu či tabulky získané<br />

experimentálně.<br />

Nyquistovo kritérium vychází z přenosu rozpojeného regulačního obvodu, který si<br />

vyjádřím ve tvaru podílu polynomů<br />

G<br />

0<br />

( s) = G ( s) G ( s)<br />

S R<br />

=<br />

(3.16)<br />

Na obr. 3.8a je ukázka rozpojeného regulačního obvodu. Pokud se na některém<br />

místě obvod fiktivně rozpojí, získají se tím samostatné dva regulační členy a také člen s<br />

otáčením znaménka v sériovém zapojení. Na vstup se přivede sinusový signál. Na výstupu<br />

lze pozorovat stejný signál, který má ovšem jinou amplitudu A a je fázově posunut o<br />

hodnotu ϕ. Pokud se fiktivně rozpojený obvod opět spojí, pak kmity, které byly poslány na<br />

vstup, se udrží a to i bez opětovného přivedení sinusového signálu. A právě teď se<br />

rozhoduje o stabilitě regulačního obvodu. Pokud se totiž tyto kmity po určité době zmírní,<br />

jedná se o stabilní regulační obvod. Pokud se naopak zesílí, jedná se o nestabilní regulační<br />

obvod. Rozhodujícím případem je posouzení regulačního obvodu na hranici jeho stability.<br />

Tento případ nastane, pokud je na vstup fiktivně rozpojeného regulačního obvodu přiveden<br />

sinusový signál a na výstupu se objeví přesně ten stejný. Jak je vidět na obr. 3.8b, regulační<br />

obvod je v sériovém zapojení s G 0 (s) a s členem otáčejícím znaménko. Na výstupu z G 0 (s)<br />

je sinusový signál stejný jako na vstupu, se stejnou amplitudou, ale fázově posunutý o<br />

180°. Poté je signál přiveden do členu otáčejícího znaménko, kde je znovu posunut o 180°<br />

a tím pádem je na výstupu totožný, a to co se týče amplitudy i fázového posunu.<br />

Budu-li uvažovat regulační obvod bez členu otáčejícího znaménko a vstupní i<br />

výstupní signál bude mít stejnou amplitudu o fázové posunutí o 180°, pak frekvenční<br />

charakteristika G 0 (s) musí procházet tzv. kritickým bodem [-1, 0].<br />

Jestliže frekvenční charakteristika rozpojeného regulačního obvodu prochází<br />

kritickým bodem [-1, 0], pak se tedy jedná o obvod na hranici stability. Jestliže při<br />

frekvenci, kde je na výstupu signál posunut o 180° oproti vstupu a amplituda výstupních<br />

kmitů je větší než vstupních, pak nedochází ke zmírnění signálu a obvod je nestabilní.<br />

M<br />

N<br />

( s)<br />

( s)


Strana 30<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

a)<br />

v=0 y<br />

G R (s)<br />

G S (s)<br />

w=0<br />

m 1 m 2<br />

b)<br />

m 1<br />

G S (s)<br />

G R (s)<br />

y -y<br />

m 2<br />

G 0 (s)<br />

Obr. 3.8 - a) rozpojený regulační obvod; b) regulační obvod v sériovém zapojení s<br />

G 0 (s) a s členem otáčejícím znaménko<br />

Im<br />

ω = ∞<br />

ω = 0<br />

Re<br />

nestabilní<br />

−1<br />

0<br />

1<br />

na hranici stability<br />

stabilní<br />

( ) jω G 0<br />

Obr. 3.9 - Průběhy amplitudo-fázové frekvenční charakteristiky rozpojeného<br />

regulačního obvodu G 0 (jω)


Použití frekvenčních charakteristik u spojitých systémů Strana 31<br />

Pro stabilní rozpojený regulační obvod lze již nyní zformulovat Nyquistovo<br />

kritérium stability:<br />

Aby byl rozpojený regulační obvod stabilní, pak jeho frekvenční<br />

charakteristika G 0 (jω) musí ležet vlevo od kritického bodu [-1, 0], a to pro frekvence<br />

ω od 0 do ∞ (obr. 3.9).<br />

Toto kritérium pojaté pro rozpojený regulační obvod se může použít pouze pro<br />

stabilní obvod.<br />

U uzavřeného regulačního obvodu nastává určité omezení použitelnosti Nyquistova<br />

kritéria. Jmenovatel přenosu N(s) nesmí obsahovat kladné kořeny (kořeny v pravé<br />

komplexní polorovině [Švarc, 2002]). Pak je možné napsat, že je-li rozpojený regulační<br />

obvod nestabilní, uzavřený regulační obvod stabilní být může. Tuto stabilitu lze vyšetřit<br />

jen zobecněným Nyquistovým kritériem.<br />

Řešení stability obvodu s dopravní zpožděním použitím zjednodušeného Nyquistova<br />

kritéria<br />

Beru-li v úvahu obvody s dopravním zpožděním, je dobré si stručně vysvětli<br />

charakter těchto obvodů. Existují situace, kdy na vstup regulačního obvodu přivedu nějaký<br />

signál, ale na výstupu se nic neobjeví. Teprve až po uplynutí určité doby se začne výstupní<br />

veličina měnit. Tomuto jevu se říká dopravní zpoždění T d .<br />

Charakteristická rovnice uzavřeného regulačního obvodu s dopravním zpožděním<br />

má tvar<br />

− sTd<br />

1 + G0 ( s)<br />

e = 0.<br />

(3.17)<br />

Dopravní zpoždění způsobuje fázové zpoždění o úhel ω Td . Modul G 0 (s) se nemění,<br />

s = . Vyjádřím-li si frekvenční přenos rozpojeného regulačního<br />

GT<br />

d<br />

protože modul ( ) 1<br />

obvodu (3.20) ve smyslu rovnic (3.18 a 3.19)<br />

( jω) = Re[ G( jω)<br />

] j G( jω)<br />

G +<br />

0<br />

jϕ<br />

(<br />

( ) ( ) ) 0 ω<br />

jω<br />

A ω e ,<br />

G =<br />

0<br />

Im[ ],<br />

j ( )<br />

( ) ( ) ϕ ω<br />

j arg G( jω<br />

jω<br />

A ω e G( jω) e<br />

) ,<br />

G =<br />

(3.18)<br />

(3.19)<br />

(3.20)<br />

potom frekvenční proces rozpojeného regulačního obvodu s dopravním zpožděním bude<br />

G<br />

T<br />

0<br />

d<br />

[ ] .<br />

j ϕ ( ω ) ϕ ( ω )<br />

( ω) ( ω) ( ω) ( ω) 1<br />

O + T d<br />

j = G j G j = A ⋅ ⋅ e<br />

0<br />

Td<br />

(3.21)<br />

Pro zjednodušené Nyquistovo kritérium je podmínkou stability uzavřeného<br />

regulačního obvodu bez dopravního zpoždění splnění podmínky<br />

0


Strana 32<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

(3.22a)<br />

(3.22b)<br />

Z podmínky (3.22b) se určí úhlová frekvence ω (obr. 3.10), při které amplitudová a<br />

fázová frekvenční charakteristika G 0 (jω) protínají reálnou osu roviny "G 0 (jω)".<br />

Při výskytu dopravního zpoždění v regulačním obvodu se vztahy (3.22) změní<br />

analogicky na<br />

(3.23)<br />

Při určitém dopravním zpoždění projde frekvenční charakteristika G T ( jω)<br />

d<br />

0<br />

kritickým bodem [-1, 0] a uzavřený obvod bude na hranici stability, jak je naznačeno na<br />

obr. 3.10.<br />

Při dalším zvětšování T d je obvod již nestabilní, frekvenční charakteristika<br />

neponechává bod [-1, 0] vlevo od svého průběhu. Hodnoty ω a T d , při kterých frekvenční<br />

charakteristika G T ( jω)<br />

d<br />

Re<br />

Im<br />

Re<br />

Im<br />

[ G0( jω)<br />

] > −1,<br />

[ G ( jω)<br />

] = 0 ⇒ ω.<br />

0<br />

− jωTd<br />

[ GR<br />

( jω) GS<br />

( jω)<br />

e ]<br />

− jωTd<br />

G ( jω) G ( jω)<br />

e<br />

> −1,<br />

[ ] = 0 ⇒ ω.<br />

R<br />

S<br />

0<br />

prochází kritickým bodem [-1, 0], se nazývá kritickými hodnotami<br />

s označením ω k a T dk .<br />

Obr. 3.10 - K řešení stability uzavřeného regulačního obvodu s dopravním<br />

zpožděním zjednodušeným Nyquistovým kritériem<br />

Kritický modul má tedy vztah<br />

A ≡ 1 =<br />

k<br />

Re<br />

2<br />

2<br />

[ G ( jω )] + Im G ( jω<br />

)<br />

o<br />

k<br />

[ ],<br />

o<br />

k<br />

(3.24)<br />

a pro kritickou fázi platí: − π = ϕ ( ω ) + ϕ ( ),<br />

k k T<br />

ω<br />

d k


Použití frekvenčních charakteristik u spojitých systémů Strana 33<br />

Im<br />

− π = arctg<br />

Re<br />

[ G0<br />

( jωk<br />

)]<br />

[ G ( jω<br />

)]<br />

0<br />

k<br />

− ω T<br />

k<br />

dk<br />

.<br />

(3.25)<br />

Ze vztahu (3.24) vypočítám kritickou úhlovou frekvenci ω k a z kritické fáze<br />

ϕ (3.25) kritické dopravní zpoždění T dk .<br />

( )<br />

k<br />

ω k<br />

Jak už jsem uvedl v kapitole 3.1.3, při zjišťování stability regulačních obvodů<br />

využívám i frekvenčních charakteristik v logaritmických souřadnicích. Tím se zabývá i<br />

Nyquistovo kritérium, které zde uvedu nyní, ale ve zjednodušené formě.<br />

Zjednodušené Nyquistovo kritérium v logaritmických souřadnicích<br />

Verzi zjednodušeného Nyquistova kritéria zde uvedu v logaritmických<br />

souřadnicích. Předpokladem je, že přenos rozpojeného regulačního obvodu nemá žádný<br />

pól v pravé polorovině roviny kořenů "s".<br />

Je-li přenos rozpojeného regulačního obvodu ve tvaru<br />

1<br />

1<br />

tj. s<br />

1<br />

= − < 0,<br />

s<br />

2<br />

= − < 0,<br />

s = 0 3<br />

,<br />

T<br />

T<br />

1<br />

G<br />

o<br />

( s)<br />

2<br />

=<br />

s<br />

k<br />

O<br />

( T s + )( T s 1) ,<br />

1<br />

1 2<br />

+<br />

(3.26)<br />

(nulový pól je zahrnut do levé poloroviny), přitom lze uvažovat relaci koeficientů přenosu<br />

(zesílení) k o1 < k o3 < k o2 .<br />

Amplitudové a fázové frekvenční charakteristiky přenosu (3.26) jsou zobrazeny v<br />

části obr. 3.10a. Pro zesílení k 01 přenosu rozpojeného regulačního obvodu je uzavřený<br />

regulační obvod stabilní, pro k 03 je na mezi stability a pro k O2 je uzavřený regulační obvod<br />

nestabilní. Rozhodl jsme tak podle polohy bodu [-1, 0] vzhledem k průběhu jednotlivých<br />

amplitudových a fázových frekvenčních charakteristik rozpojeného regulačního obvodu.<br />

Z obr. 3.11 vyplývá souvislost průběhů frekvenčních charakteristik rozpojeného<br />

regulačního obvodu v komplexní rovině a v logaritmických souřadnicích. V části obr.<br />

3.11b je amplitudová logaritmická frekvenční charakteristika, v části obr. 3.11c je<br />

zakreslena fázová logaritmická frekvenční charakteristika.<br />

Frekvenčním bodem pro rozhodování o stabilitě zjednodušeným Nyquistovým<br />

kritériem v komplexní rovině byl bod [-1, 0]. V tomto bodě modul frekvenčního přenosu<br />

rozpojeného regulačního obvodu má hodnotu<br />

G O<br />

( jω) = 1,<br />

(3.27)<br />

a fáze má hodnotu ϕ ( ω) = −180° .<br />

V logaritmických souřadnicích si zobrazím odpovídající kritický bod výpočtem<br />

[ ] = 20 log1 = 0.<br />

A dB<br />

(3.28)


Strana 34<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

Jednotková kružnice v komplexní rovině se tedy zobrazí do přímky s hodnotou 0<br />

dB (viz obr. 3.11b).<br />

Fáze má hodnotu stejnou, zůstává nezměněna. Kritickým bodem pro rozhodování o<br />

stabilitě zjednodušeným Nyquistovým kritériem v logaritmických souřadnicích je bod v<br />

amplitudové logaritmické frekvenční charakteristice s amplitudou A[dB] = 0 a nazývá se<br />

úhlovým kmitočtem řezu ω ř .<br />

Při praktické aplikaci teorie automatického řízení je možné narazit na systémy, u<br />

kterých je splněna silná podmínka fyzikální realizovatelnosti, tj. n > m. Z toho plyne, že<br />

stupeň jmenovatele přenosu rozpojeného regulačního obvodu je vyšší než stupeň čitatele a<br />

tedy fáze ϕ má zápornou hodnotu. To je případ i přenosu (3.26), který představuje<br />

integrační člen se setrvačností 2. řádu. Jeho fázová logaritmická frekvenční charakteristika<br />

je zobrazena na obr. 3.11c. Porovnáním hodnot fáze kritického bodu v částech obr. 3.11c a<br />

obr. 3.11a je vidět, že je možné fázi kritického bodu stanovit na hodnotu<br />

ϕ<br />

( ω) = −180° .<br />

(3.29)<br />

Potom pro ( ω) > −180°<br />

pro ( ω) < −180°<br />

ϕ (tj. např. pro − 160° , …) bude uzavřený obvod stabilní a naopak<br />

ϕ (tj. např. pro − 200 ° , …) bude uzavřený regulační obvod nestabilní.<br />

Definice stability podle zjednodušeného Nyquistova kritéria v logaritmických<br />

souřadnicích<br />

Uzavřený regulační obvod je stabilní, jestliže pro úhlovou frekvenci řezu ω ř , při<br />

které; amplitudová logaritmická frekvenční charakteristika rozpojeného regulačního<br />

obvodu protíná osu 0 dB, bude ϕ ( ω) > −180°<br />

. Důsledkem je, že amplitudová logaritmická<br />

frekvenční charakteristika protíná osu 0 dB vlevo od úhlové frekvence, při němž má fáze<br />

hodnotu ϕ ( ω) = −180°<br />

(viz obr. 3.11).<br />

Z uvedeného je zřejmá výhoda používání frekvenčních charakteristik v<br />

logaritmických souřadnicích. Jestliže vím, že místa lomu amplitudové charakteristiky jsou<br />

definována převrácenou hodnotou časových konstant přenosu, snadno poznám, kterou z<br />

časových konstant přenosu rozpojeného regulačního obvodu musím změnit, aby uzavřený<br />

regulační obvod byl stabilní.


Použití frekvenčních charakteristik u spojitých systémů Strana 35<br />

a)<br />

b)<br />

c)<br />

Obr. 3.11 - Souvislost zjednodušeného Nyquistova kritéria v komplexní rovině a<br />

v logaritmických souřadnicích: a) v komplexní rovině; b) amplitudová<br />

logaritmická frekvenční charakteristika; c) fázová logaritmická<br />

frekvenční charakteristika


Strana 36<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

3.3 Syntéza spojitých systémů frekvenčními metodami<br />

Byla vypracována celá řada metod syntézy regulačního obvodu. Jedny z metod<br />

tvoří klasické metody syntézy pomocí frekvenčních charakteristik. Tyto metody se<br />

používají pro syntézu struktury i parametrů regulátorů v regulačních obvodech.<br />

Při syntéze regulačního obvodu pomocí frekvenčních charakteristik se vychází z<br />

dané frekvenční charakteristiky rozpojeného regulačního obvodu nebo jeho části. Syntéza<br />

se provádí v komplexní rovině nebo v logaritmických souřadnicích.<br />

Frekvenční metody syntézy v podstatě využívají toho, že nevhodné póly i nuly<br />

přenosu rozpojené smyčky se kompenzují nulami a póly přenosu korekčních členů<br />

(korekčními členy bývají nejčastěji standardní lineární regulátory typu P, PI, PD, PID<br />

apod.).<br />

3.3.1 Ukazatele kvality regulace<br />

Kvalitu regulace je možné posuzovat podle ukazatelů v časové, frekvenční nebo<br />

obrazové oblasti [Davidová, 2004]. Z průběhu přechodové charakteristiky se kvalita<br />

regulace v časové oblasti určuje z odezvy regulačního obvodu na skokový signál. Právě z<br />

této přechodové charakteristiky se určí kvantitativní ukazatele, díky kterým lze hodnotit<br />

chování regulačního obvodu. Mezi kvantitativní ukazatele patří maximální překmit y max<br />

regulované veličiny, doba regulace t<br />

r<br />

, trvalá regulační odchylka e<br />

T<br />

a perioda kmitů<br />

T přechodové charakteristiky. Pro určení kvality regulace v obrazové oblasti je zásadní<br />

poloha pólů uzavřeného regulačního obvodu v komplexní rovině.<br />

Kvantitativní ukazatele, které jsem zmínili výše, jsou zobrazeny na obr. 3.12 a ve<br />

stručnosti si zde uvedeme jejich popis.<br />

• maximální překmit y max regulované veličiny = maximální hodnota, kterou<br />

dosáhne regulovaná veličina<br />

• doba regulace t r = doba, za kterou trvale klesne regulační odchylka pod 5%<br />

ustálené hodnoty<br />

• trvalá regulační odchylka e T (t) = rozdíl mezi požadovanou a skutečnou hodnotou<br />

regulované veličiny<br />

• perioda kmitů T přechodové charakteristiky.<br />

Někdy se také používají další ukazatele kvality regulace, jako např.:<br />

• čas t max = doba, kdy dochází k největšímu překmitu<br />

• počet kmitů.<br />

Tyto ukazatelé kvality dávají přímou informaci o časovém průběhu regulačního<br />

pochodu, což je jejich výhodou, ale zároveň platí podmínka, že přechodová charakteristika<br />

musí mít kmitavý průběh, nikoliv aperiodický, jinak maximální překmit ani periodu kmitů<br />

nelze určit.


Použití frekvenčních charakteristik u spojitých systémů Strana 37<br />

y(t)<br />

+5% y(∞)<br />

y(∞)<br />

y max<br />

T<br />

–5% y(∞)<br />

t max<br />

t<br />

t r<br />

Obr. 3.12 - Ukazatele kvality regulace na přechodové charakteristice<br />

Pro kvalitu regulace v obrazové oblasti je zásadní poloha pólů uzavřeného<br />

regulačního obvodu v komplexní rovině. Mezi všemi uvedenými ukazateli kvality existuje<br />

jednoznačná souvislost.<br />

Především jsou důležité dvě hodnoty, a to maximální překmit y max regulované<br />

soustavy a doba regulace t r . Je dobré si představit, že jsou regulovány např. otáčky<br />

nějakého točivého stroje, pak y max je hodnotou, na kterou musí být stroj ještě pevnostně<br />

dimenzován a otáčky y max nesmí např. padnout do kritických otáček stroje. Dalším<br />

požadavkem je, aby doba regulace t r nebyla příliš dlouhá a perioda kmitů T také někdy<br />

nesmí padnout do určitého pásma.<br />

Při regulačním pochodu v podstatě dochází k výměně energie. V oblasti<br />

podregulování má regulovaná soustava nedostatek energie a je zapotřebí jí energii dodávat,<br />

aby se hodnota regulované veličiny zvýšila na požadovanou hodnotu. V oblasti<br />

přeregulování má naopak přebytek energie a tu předává okolí, aby se hodnota regulované<br />

veličiny snížila na požadovanou hodnotu. To pochopitelně není možné vzhledem např. k<br />

setrvačným hmotám apod., tedy je alespoň požadováno, aby výměna energie byla co<br />

nejmenší – aby regulační plocha byla minimální.<br />

Regulační plocha, zmíněná v předešlém odstavci, je podle obr. 3.13 vyšrafovaná<br />

plocha nad a nebo pod průběhem regulované veličiny.<br />

Minimum regulační plochy znamená také kompromis mezi protichůdnými<br />

požadavky na minimální y max a minimální t r . Tam se totiž ukazují protiklady požadavků na<br />

tyto dvě hodnoty. Čím větší je např. zesílení regulátoru r 0 , tím kratší je sice doba regulace,<br />

ale tím větší je maximální překmitnutí. Optimum mezi těmito dvěma protichůdnými<br />

požadavky je minimum regulační plochy. V praxi se kvalita regulace hodnotí podle<br />

charakteristických parametrů regulačního pochodu (viz obr. 3.13).


Strana 38<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

y(t)<br />

y(∞)<br />

+ +<br />

– –<br />

t<br />

Obr. 3.13 - Ukázka regulační plochy<br />

Vzhledem k tématu mé diplomové práce, ve které se zabývám frekvenčními<br />

metodami, bude dále má pozornost upřena na ukazatele kvality regulace ve frekvenční<br />

oblasti.<br />

Ukazatelé kvality regulace ve frekvenční oblasti<br />

Pro hodnocení kvality řízení budu na frekvenční charakteristice uzavřené smyčky<br />

definovat následující míry:<br />

• rezonanční převýšení A r = maximální hodnota zesílení (maximální hodnota<br />

amplitudy frekvenčního přenosu). Velké rezonanční převýšení znamená velký<br />

překmit na přechodové charakteristice.<br />

• rezonanční frekvence ω r = frekvence, při níž frekvenční charakteristika dosahuje<br />

hodnoty rezonančního převýšení A r a rozhoduje o relativní stabilitě uzavřeného<br />

obvodu<br />

• šířka pásma propustnosti ω h = frekvence, na níž poklesne zesílení o 3 dB oproti<br />

zesílení na nízkých frekvencích. Širší propustné pásmo znamená rychlejší odezvu<br />

systému, tj. kratší dobu náběhu přechodové charakteristiky (dobu, za kterou přejde<br />

výstup z 10% na 90% ustálené hodnoty). Na druhou stranu vyšší mezní frekvence<br />

však znamená, že systém může reagovat i na vysokofrekvenční rušení zpravidla<br />

přítomné na vstupech systému.<br />

Tyto ukazatele kvality jsou zobrazeny na obr. 3.14, kde je vynesena typická<br />

frekvenční charakteristika uzavřeného regulačního obvodu.


Použití frekvenčních charakteristik u spojitých systémů Strana 39<br />

Existuje souvislost mezi rezonančním převýšením A r amplitudové charakteristiky a<br />

maximálním překmitem y max na přechodové charakteristice. I přes to, že platí čím větší je<br />

rezonanční převýšení, tím větší je maximální překmit, tak většina regulačních obvodů se v<br />

praxi navrhuje tak, aby A<br />

r<br />

∈ ( 1.1,1.5 ), resp. A rdB<br />

∈ ( 1dB,<br />

3dB)<br />

, a to z důvodu urychlení<br />

odezvy. Stejně tak existuje souvislost mezi dobou regulace t r a šířkou pásma propustnosti<br />

ω h , kde platí nepřímá úměrnost. Čím větší je šířka pásma propustnosti, tím rychlejší je<br />

přechodový děj a tím pádem je kratší doba regulace. Doba regulace je omezena vztahem,<br />

neboli nerovností<br />

π 4π<br />

< t<br />

r<br />

< .<br />

ω ω<br />

(3.30)<br />

ř<br />

Pro připomenutí ω ř je úhlová frekvence řezu, při níž je absolutní hodnota<br />

frekvenčního přenosu rovna jedné.<br />

ř<br />

A<br />

[dB]<br />

A r<br />

ω r<br />

ω h<br />

3dB<br />

ω<br />

Obr. 3.14 - Ukazatele kvality regulace z frekvenční charakteristiky<br />

uzavřeného regulačního obvodu<br />

Pro frekvenční návrh uzavřeného regulačního obvodu jsou na frekvenční<br />

charakteristice rozpojeného regulačního obvodu definovány ukazatele kvality regulace<br />

uvedené v [Davidová, 2004]. Jsou to tyto:<br />

• amplitudová bezpečnost δ říká, kolikrát se ještě může zvětšit zesílení v rozpojené<br />

smyčce, než se zpětnovazební systém dostane na mez stability<br />

• fázová bezpečnost γ je dána úhlem frekvenční charakteristiky γ otevřeného<br />

regulačního obvodu při frekvenci ω ř , při kterém se amplituda frekvenční<br />

charakteristiky rovná jedné. Fázová bezpečnost je záporně vzatá změna fáze<br />

otevřeného obvodu, která přivede uzavřený regulační obvod na hranici stability.<br />

Tyto ukazatele kvality regulace jsou znázorněny na obr. 3.15, kde je vyobrazena<br />

frekvenční charakteristika rozpojeného regulačního obvodu.


Strana 40<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

γ<br />

a) b)<br />

Im<br />

A<br />

[dB]<br />

1/δ<br />

–1<br />

0<br />

0<br />

–20<br />

Re<br />

–40<br />

( ) jω G 0<br />

ωř<br />

–120°<br />

–150°<br />

–180°<br />

–210°<br />

γ<br />

ωř<br />

δ<br />

ω<br />

f<br />

ω<br />

ω<br />

ϕ<br />

[°]<br />

Obr. 3.15 - a) frekvenční charakteristika v komplexní rovině; b) amplitudová a<br />

fázová charakteristika v logaritmických souřadnicích<br />

A<br />

[dB]<br />

–40 dB/dek<br />

–20 dB/dek<br />

∆L<br />

ω 1<br />

–20 dB/dek<br />

ω ř<br />

ω 2<br />

∆L<br />

ω<br />

–40 dB/dek<br />

–60 dB/dek<br />

Obr. 3.16 - Logaritmická frekvenční charakteristika ve formě asymptot<br />

Na obr. 3.16 je vyobrazena logaritmická frekvenční charakteristika ve formě<br />

asymptot. V této formě je to nejčastěji používaná logaritmická charakteristika, protože její<br />

konstrukce je jednoduchá a rychlá. Frekvence ω 1 a ω 2 se nazývají lomové frekvence, tvoří<br />

průsečíky asymptot a jsou na frekvencích, které jsou převrácenými hodnotami časových<br />

konstant regulátoru nebo regulované soustavy.


Použití frekvenčních charakteristik u spojitých systémů Strana 41<br />

3.3.2 Frekvenční metody syntézy<br />

U syntézy uzavřeného regulačního obvodu se upravuje frekvenční charakteristika<br />

rozpojeného obvodu tak, aby výsledná frekvenční charakteristika uzavřené smyčky měla<br />

požadovaný průběh, tedy aby byly splněny ukazatele kvality regulace. Frekvenční<br />

charakteristiku rozpojeného regulačního obvodu můžeme upravovat dvojím způsobem:<br />

• Volbou zesílení rozpojeného regulačního obvodu. Tím vlastně určím konstantu<br />

proporcionálního regulátoru (kap. 3.3.2.1).<br />

• Sériovými korekčními členy. Ideálním sériovým korekčním členem jsou vlastně<br />

regulátory P, PD, PI, popř. PID (kap. 3.3.2.2).<br />

3.3.2.1 Volba zesílení rozpojeného regulačního obvodu<br />

Pro zvolené rezonanční převýšení A r frekvenčního přenosu uzavřeného regulačního<br />

obvodu můžu určit odpovídající zesílení rozpojeného regulačního obvodu. V komplexní<br />

rovině provedu jednoduchou konstrukci (Brown - Campbellova konstrukce).<br />

Brown - Campbellova konstrukce je zobrazena na obr. 3.17. Nakreslím frekvenční<br />

charakteristiku rozpojeného regulačního obvodu se zvoleným zesílením k o . Z počátku vedu<br />

tečnu ke kružnici zvoleného rezonančního převýšení A r . Pak snadno odvodím, že pro úhel<br />

α se zápornou reálnou poloosou platí<br />

1<br />

α = arcsin , . ≥ 1<br />

(3.31)<br />

A r<br />

A r<br />

Sestrojím kružnici se středem na záporné reálné ose tak, aby se dotýkala v jednom bodě<br />

nakreslené frekvenční charakteristiky (bod A) a měla tečnu sestrojenou v bodě B. V bodě<br />

dotyku kružnice a tečny (bod B) vedu kolmici na reálnou osu, přečtu její vzdálenost od<br />

počátku, označím ji x a zesílení přenosu rozpojeného obvodu vyjádřím vztahem<br />

k<br />

oopt<br />

= k<br />

o<br />

1<br />

x<br />

. (3.32)<br />

Im<br />

S<br />

x<br />

A<br />

–1<br />

α<br />

0<br />

Re<br />

B<br />

( ) jω G 0<br />

Obr. 3.17 - Brown - Campbellova konstrukce


Strana 42<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

Ještě snadněji určím zesílení rozpojené smyčky v logaritmických souřadnicích.<br />

Sestrojím frekvenční charakteristiku rozpojeného regulačního obvodu pro zvolené zesílení<br />

k v logaritmických souřadnicích. Tuto charakteristiku přenesu do Nicholsova diagramu.<br />

Frekvenční charakteristiku kreslím na průsvitku položenou na Nicholsův diagram. Poté<br />

posunu průsvitku s nakreslenou charakteristikou ve směru osy, na nichž je vynesena<br />

amplituda v decibelech přenosu rozpojené smyčky tak, až se v Nicholsově diagramu<br />

dotkne izoplety s hodnotou A rdB . Stejně jako v předchozím případě velikost posunu<br />

označím x [dB]. Optimální zesílení přenosu rozpojené smyčky je<br />

nebo<br />

( k<br />

opt<br />

) = k<br />

dB<br />

+ xdB<br />

dB<br />

k opt<br />

= kx<br />

(3.33)<br />

(3.34)<br />

Zde se hodí uvést také zjednodušený postup. Vyhovující hodnotu zesílení rozpojené<br />

smyčky určím přibližně také tak, že zajistím, aby amplitudová logaritmická frekvenční<br />

charakteristika rozpojené smyčky procházela osou 0 dB při fázi ϕ = −135°<br />

. Tím zajistím<br />

fázovou bezpečnost 45° a hodnota rezonančního převýšení tak bude kolem 1,3. Při tomto<br />

zjednodušení nemusím ani používat Nicholsův diagram. Posun o x decibelů provedu přímo<br />

v logaritmických amplitudových souřadnicích rozpojeného regulačního obvodu.<br />

3.3.2.2 Sériové korekční členy<br />

Jak jsem již zmínil v kapitole 3.3.2, sériovým korekčním členem jsou regulátory P,<br />

PD, PI, popř. PID nebo jimi také mohou být obecné korekční členy nahrazující tyto<br />

regulátory. Jsou to členy, které jsou zapojeny v sérii s ostatními členy a tvoří přímou větev<br />

regulačního obvodu. Syntéza regulačního obvodu ve frekvenční oblasti využívá vlastností<br />

typizované logaritmických amplitudových frekvenční charakteristik (TLAFCH),<br />

především však skutečnosti, že výsledná (TLAFCH) je dána geometrickým součtem<br />

(TLAFCH) jednotlivých členů zařazených do série. Tato skutečnost nám umožňuje snadné<br />

určení korekčního členu.<br />

3.3.2.3 Metoda typizované logaritmické frekvenční charakteristiky<br />

Vlastnosti uzavřeného regulačního obvodu jsou plně určeny frekvenčními<br />

charakteristikami jeho rozpojeného obvodu [Davidová, 2004]. Při návrhu regulačního<br />

obvodu se bere často v úvahu metody využívající frekvenční charakteristiky v<br />

logaritmických souřadnicích, díky kterým je k dispozici lepší náhled na chování systému v<br />

případě, že se některý z jeho parametrů bude měnit. Nevýhodou této i ostatních<br />

frekvenčních metod je fakt, že na časové závislosti (přechodová charakteristika apod.) lze<br />

usuzovat pouze podle nepřímých ukazatelů a k získání odezev je nutné volit zdlouhavý<br />

postup [Suchna, 2005].<br />

Nyní zkusím uvést stručný popis metody. Základem je tzv. žádaná logaritmická<br />

amplitudová frekvenční charakteristika rozpojeného obvodu, která předem splňuje<br />

stanovené požadavky na přesnost a kvalitu regulačního obvodu.


Použití frekvenčních charakteristik u spojitých systémů Strana 43<br />

Podle pravidel, která uvedu dále, můžu sestrojit logaritmickou amplitudovou<br />

charakteristiku rozpojeného obvodu G 0 (jω) dB . Ta splňuje dané požadavky na regulaci (obr.<br />

3.18).<br />

Vezmu-li v úvahu přenos rozpojeného regulačního obvodu G 0 (jω), pro který platí<br />

vztah<br />

a z toho<br />

G<br />

0<br />

( jω) = G ( jω) G ( jω)<br />

G<br />

R<br />

( jω)<br />

R<br />

G<br />

=<br />

G<br />

0<br />

S<br />

S<br />

( )<br />

( )<br />

(3.35)<br />

jω<br />

. (3.36)<br />

jω<br />

Pak pro logaritmické souřadnice platí vztahy<br />

G<br />

R<br />

( jω) = G ( jω) − G ( jω) ,<br />

dB<br />

ϕ<br />

0 dB S dB<br />

ϕ<br />

R<br />

= 0<br />

−ϕ S<br />

(3.37)<br />

(3.38)<br />

Logaritmická amplitudová frekvenční charakteristika se tedy stanoví z rozdílů těchto<br />

charakteristik pro rozpojený regulační obvod a regulovanou soustavu.<br />

Metoda typizované logaritmické amplitudové frekvenční charakteristiky je<br />

vhodným příkladem toho, jak lze frekvenční charakteristiku využít pro návrh spojitého<br />

regulačního obvodu.<br />

A<br />

[dB]<br />

–20 dB/dek<br />

G 0 (jω)<br />

0 dB/dek<br />

–40 dB/dek<br />

ω 1<br />

–20 dB/dek<br />

–20 dB/dek<br />

G S (jω)<br />

ω2<br />

ω3<br />

ω<br />

–20 dB/dek<br />

G R (jω)<br />

0 dB/dek<br />

ω ř<br />

+20 dB/dek<br />

–20 dB/dek<br />

–40 dB/dek<br />

–60 dB/dek<br />

Jak je znázorněno na obr. 3.19, typizovanou logaritmickou amplitudovou<br />

frekvenční charakteristiku lze rozdělit na tři oblasti, a to nízkofrekvenční, středofrekvenční<br />

a vysokofrekvenční.<br />

Nízkofrekvenční oblast je úsek v rozsahu úhlových frekvencí ω ∈< 0,<br />

ω1<br />

> , kde<br />

ω = 1 T je první zlom frekvenční charakteristiky. Tato část je v logaritmických<br />

1<br />

/<br />

Obr. 3.18 - Typizovaná logaritmická amplitudová frekvenční charakteristika<br />

max


Strana 44<br />

Použití frekvenčních charakteristik u spojitých systémů<br />

souřadnicích kreslena se sklony 0 dB/dek, –20 dB/dek, –40 dB/dek. Asymptota<br />

nízkofrekvenční části charakteristiky musí procházet hodnotou 20 log k 0 při frekvenci<br />

ω = 1, kde k 0 je zesílení rozpojeného regulačního obvodu. Tato část charakteristiky určuje<br />

typ regulačního obvodu a mé vliv na přesnost regulačního obvodu v ustáleném stavu.<br />

Středofrekvenční oblast určuje kvalitu přechodového děje. Je to úsek v okolí<br />

průsečíku amplitudové frekvenční charakteristiky rozpojeného regulačního obvodu s osou<br />

0 dB při úhlové frekvenci ω<br />

ř<br />

. Na obr. 3.18 je šířka této části omezena úhlovými<br />

frekvencemi ω ∈< ω ,ω 2 3<br />

> . Z důvodu optimálního průběhu přechodové charakteristiky<br />

uzavřeného regulačního obvodu, je požadovaný skon charakteristiky v okolí úhlové<br />

frekvence ω –20 dB/dek. Úhlovou frekvenci ω stanovíme přibližně podle vztahu<br />

ř<br />

4π<br />

ω , (3.39)<br />

ř<br />

=<br />

tr<br />

kde t r je požadovaná doba regulace.<br />

Šířka středofrekvenční oblasti se volí v rozsahu 0,5 až 0,9 dekády. Čím je šířka oblasti<br />

větší, tím je regulační obvod více tlumen a zároveň je větší i jeho fázová bezpečnost γ.<br />

Vysokofrekvenční oblast začíná pro úhlové frekvence ω = (6 až 8)ω ř . Tato část<br />

nemá podstatný vliv na průběh přechodového děje uzavřeného regulačního obvodu. Její<br />

sklon bývá obvykle –40 dB/dek nebo –60 dB/dek.<br />

Spojení nízkofrekvenční a středofrekvenční oblasti bývá se sklonem asymptoty –40<br />

dB/dek až –60 dB/dek. Stejně tak spojení středofrekvenční a vysokofrekvenční části mívá<br />

sklon asymptot –40 dB/dek až –60 dB/dek.<br />

Metoda typizované amplitudové logaritmické frekvenční charakteristiky<br />

rozpojeného regulačního obvodu se spíše aplikuje při návrhu servomechanismů než u<br />

klasických regulačních obvodů složených z dvojice regulátor - regulovaná soustava. V<br />

tomto případě je totiž struktura obvodu plně dána a volba typu regulátoru a jeho parametrů<br />

je omezen.<br />

ř<br />

A<br />

[dB]<br />

–40 dB/dek<br />

–20 dB/dek<br />

0 dB/dek –40 dB/dek<br />

∆L<br />

ω 1<br />

ω 2<br />

–20 dB/dek<br />

ω ř<br />

ω 3<br />

ω 4<br />

∆L<br />

–40 dB/dek<br />

ω<br />

nízkofrekvenční<br />

oblast<br />

středofrekvenční<br />

oblast<br />

–60 dB/dek<br />

vysokofrekvenční<br />

oblast<br />

Obr. 3.19 - Typizovaná logaritmická amplitudová frekvenční charakteristika


Strana 45<br />

4 POUŽITÍ FREKVENČNÍCH CHARAKTERISTIK U<br />

DISKRÉTNÍCH SYSTÉMŮ<br />

Diskrétní regulační obvody jsou takové, v nichž alespoň jeden člen pracuje<br />

diskrétně, tj. informaci přijímá nebo vydává, eventuálně obojí, v diskrétních časových<br />

okamžicích (zpravidla rovnoměrných - ekvidistantních). Jinými slovy, alespoň jedna<br />

veličina obvodu má tvar posloupnosti diskrétních hodnot. [Švarc, 2002]<br />

Při analýze a syntéze diskrétních regulačních obvodů se zpravidla pracuje s<br />

diskrétními modely spojitých částí regulačních obvodů.<br />

4.1 Frekvenční popis diskrétních systémů<br />

V této části čtvrté kapitoly se budu zabývat frekvenčním popisem diskrétních<br />

regulačních obvodů. Vysvětlím, které pojmy a výpočtové vztahy využiji při analýze a<br />

syntéze regulačních obvodů. Důležité rovněž bude ukázat si postup při konstrukci<br />

frekvenčních charakteristik.<br />

4.1.1 Frekvenční přenos<br />

Nejprve je vhodné uvažovat lineární diskrétní systém, jehož Z - přenos má tvar s<br />

kladnými mocninami z<br />

G<br />

( z)<br />

bm<br />

z<br />

=<br />

a z<br />

+ ... + b1<br />

z + b<br />

+ ... + a z + a<br />

Frekvenční přenos diskrétního systému se potom definuje vztahem<br />

G<br />

n<br />

m<br />

n<br />

Y<br />

ω =<br />

U<br />

( j T )<br />

(4.1)<br />

(4.2)<br />

což je podíl Fourierova obrazu výstupní a Fourierova vstupní funkce. Tento frekvenční<br />

přenos je komplexní funkce bezrozměrné frekvence ωT. Získáme ho ze Z - přenosu G(z)<br />

dosazením<br />

jωT<br />

z = e<br />

(4.3)<br />

tedy<br />

G j T = G z<br />

(4.4)<br />

Frekvenční přenos diskrétního systému je (na rozdíl od frekvenčního přenosu<br />

spojitého systému) periodická funkce frekvence ωT s periodou 2π<br />

1<br />

( jωT<br />

)<br />

( jωT<br />

)<br />

( ) ( ) jωT<br />

ω<br />

z=<br />

e<br />

( jωT<br />

) = G[ j( ωT<br />

+ kπ<br />

)]<br />

G 2<br />

0<br />

0<br />

(4.5)


Strana 46<br />

Použití frekvenčních charakteristik u diskrétních systémů<br />

4.1.2 Frekvenční charakteristika v komplexní rovině<br />

Frekvenční charakteristika diskrétního systému v komplexní rovině je grafické<br />

znázornění frekvenčního přenosu G(jωT) v závislosti na bezrozměrné frekvenci ωT, která<br />

se mění obecně od 0 do 2π (z důvodu toho, že se jedná o periodickou funkci, její průběh by<br />

se opakoval) a znázorňuje se v komplexní rovině. Snadno lze zjistit, že frekvenční<br />

charakteristika je souměrná podle reálné osy [Švarc, 2002]. Vzhledem k této souměrnosti a<br />

k periodičnosti G(jωT) ji stačí znázorňovat a počítat v rozsahu<br />

0 ≤ ωT<br />

≤ π<br />

(4.6)<br />

Na obr. 4.1 je právě znázorněn příklad zobrazení frekvenční charakteristiky v<br />

komplexní rovině diskrétního regulačního obvodu.<br />

a) b)<br />

Im<br />

Im<br />

ωT=π/2<br />

ωT=2<br />

Re<br />

ωT=π/2<br />

ωT=0<br />

Re<br />

ωT=π<br />

Obr. 4.1 - Frekvenční charakteristika v komplexní rovině<br />

4.1.3 Logaritmická frekvenční charakteristika<br />

Jak je uvedeno v [Davidová, 2002], k sestrojení amplitudové a fázové frekvenční<br />

charakteristiky v logaritmických souřadnicích se využívá vztah<br />

pak<br />

G<br />

1+<br />

jΩ<br />

z =<br />

, (4.7)<br />

1−<br />

jΩ<br />

( jΩ) = G( z)<br />

1+<br />

jΩ<br />

z=<br />

1−<br />

jΩ<br />

Pro vyjádření amplitudové části frekvenční charakteristiky G( jΩ ) dB<br />

lze použít asymptot<br />

jako u spojitých systémů, protože<br />

(4.8)


Použití frekvenčních charakteristik u diskrétních systémů Strana 47<br />

Ω = tg ωT<br />

2<br />

(4.9)<br />

a frekvenčnímu rozsahu ω T ∈< 0,<br />

π > odpovídá frekvenční rozsah modifikované<br />

frekvence Ω ∈< 0 , ∞ > .<br />

4.2 Analýza diskrétních systémů frekvenčními metodami<br />

V této kapitole budou analyzovány vlastnosti diskrétních systémů. Stejně jako u<br />

spojitých sytému, analýza diskrétních systémů se orientuje na stabilitu regulačního obvodu.<br />

Budou analyzovány různé popisy řízených systémů a uvedu zde některé metody určení<br />

vlastností a stability systému.<br />

4.2.1 Stabilita regulačních obvodů<br />

Pro diskrétní systémy platí stejná definice stability jako pro každý lineární systém:<br />

Systém je stabilní, jestliže se po odeznění budicího signálu vrátí do rovnovážného<br />

stavu. Názorně zobrazeno je to na obr. 4.2 a matematicky zapsáno (už pro diskrétní<br />

systémy)<br />

( k ) 0<br />

lim y hom<br />

=<br />

k →∞<br />

(4.10)<br />

y hom (k)<br />

k<br />

Obr. 4.2<br />

Tato definice bývá pro diskrétní systémy často uváděna v trochu odlišné formě, a to:<br />

Obvod je stabilní, jestliže odezva na omezenou (konečnou) vstupní veličinu je opět<br />

omezená (konečná) výstupní veličina.<br />

Vezmu-li v úvahu diskrétní regulační obvod podle obr 4.3., jeho přenos řízení je dán<br />

vztahem


Strana 48<br />

Použití frekvenčních charakteristik u diskrétních systémů<br />

Gw<br />

( z)<br />

Y<br />

=<br />

U<br />

( z)<br />

( z)<br />

GR<br />

=<br />

1+<br />

G<br />

( z) GS<br />

( z)<br />

( z) G ( z)<br />

R<br />

S<br />

(4.11)<br />

T<br />

y(k)<br />

T<br />

T<br />

w(t) e(t) G G S (s)<br />

y(t)<br />

R (z)<br />

Obr. 4.3<br />

Tento vztah je možné si vyjádřit jako podíl polynomů s kladnými mocninami z<br />

G<br />

w<br />

( z)<br />

Y<br />

=<br />

U<br />

( z)<br />

( z)<br />

m<br />

( z) GS<br />

( z)<br />

bm<br />

z + ... + b1<br />

z +<br />

=<br />

n<br />

R<br />

( z) GS<br />

( z) anz<br />

+ ... a1z<br />

+<br />

0<br />

GR<br />

b<br />

=<br />

1+<br />

G<br />

a<br />

0<br />

(4.12)<br />

Pro odvození obecné podmínky stability, převedu Z - přenos řízení (4.12) na diferenční<br />

rovnici řízení<br />

a<br />

y<br />

( k + n) + ... + a y( k + 1) + a y( k ) = b w( k + m) + ... + b w( k + ) + b w( k )<br />

n 1 0<br />

m<br />

1<br />

1<br />

0<br />

(4.13)<br />

Na základě této rovnice mohu stanovit průběh regulované veličiny y(k) při různých<br />

změnách řídicí veličiny eventuálně poruchy. Přechodný děj při těchto změnách je dán<br />

řešením této rovnice a toto řešení má tvar<br />

( k ) y ( k ) + y ( k )<br />

kde y hom<br />

( k)<br />

je řešení homogenní rovnice a ( k)<br />

y<br />

= hom<br />

(4.14)<br />

y part<br />

je partikulární integrál daný pravou<br />

stranou rovnice (4.13). Partikulární integrál závisí čistě na vstupující řídicí nebo poruchové<br />

veličině a na stabilitu nemá vliv, protože stabilita je posuzována až po skončení působení<br />

vzruchu, který vyvedl regulační obvod z rovnovážného stavu. Důležitá je tedy pouze první<br />

část řešení diferenční rovnice a to je y hom<br />

( k)<br />

. Řeším tedy homogenní rovnici k rovnici<br />

(4.13)<br />

part<br />

a n<br />

0<br />

y<br />

( k + n) + + a y( k + 1) + a y( k) 0<br />

...<br />

0<br />

1<br />

=<br />

(4.15)<br />

K získání tohoto řešení nejprve sestavím charakteristickou rovnici k dané diferenční<br />

rovnici (budu zde používat symbol z)<br />

a<br />

n<br />

z<br />

n<br />

+ ... + a z + a0<br />

1<br />

=<br />

(4.16)


Použití frekvenčních charakteristik u diskrétních systémů Strana 49<br />

a určím kořeny z 1 , z 2 , … z n .<br />

Charakteristická rovnice (4.13) respektive její levá strana se objevuje ve<br />

jmenovateli přenosu řízení (4.12) a proto je charakteristická rovnice také<br />

( z) G ( ) 0<br />

1 + G z =<br />

R<br />

S<br />

(4.17)<br />

Jsou-li kořeny charakteristické rovnice všechny různé (násobné kořeny), je řešení<br />

homogenní diferenční rovnice (4.15) dáno rovnicí<br />

k k<br />

k<br />

( k ) = C z + C z + ... C z<br />

y +<br />

hom<br />

1<br />

1<br />

2<br />

2<br />

n<br />

n<br />

(4.18)<br />

Vztah pro limitu tohoto řešení je<br />

lim y hom<br />

k =<br />

Aby byla splněna podmínka stability ( ) 0<br />

rovnice<br />

k →∞<br />

k<br />

k<br />

k<br />

( k ) = C lim z + C lim z + ... + C lim z<br />

lim yhom<br />

1<br />

2<br />

1<br />

k →∞<br />

k→∞<br />

2<br />

k→∞<br />

n<br />

k→∞<br />

n<br />

(4.19)<br />

, musí být kořeny charakteristické<br />

z i<br />

< 1<br />

pro i = 1, 2, … n (4.20)<br />

Určením kořenů z1 , z2 , … zn. je možné<br />

vyslovit obecnou podmínku stability<br />

diskrétních regulačních obvodů, která se<br />

trochu liší od obecné podmínky stability<br />

spojitých obvodů a je uvedena v [Švarc,<br />

2003]: Diskrétní obvod je stabilní, ležíli<br />

kořeny charakteristické rovnice<br />

1 + G R (z)G S (z) = 0 uvnitř jednotkové<br />

kružnice. Tato definice je graficky<br />

znázorněna na obr. 4.4<br />

nestabilní<br />

oblast<br />

–1<br />

j<br />

Im<br />

stabilní<br />

oblast<br />

na hranici<br />

stability<br />

1<br />

Re<br />

–j<br />

Obr. 4.4<br />

4.2.2 Michajlov-Leonhardovo kritérium<br />

Michajlov - Leonhardovo kritérium pro diskrétní regulační obvod se také někdy<br />

označuje jako křivkové kritérium. Stejně jako u spojitých systému je z charakteristické<br />

rovnice (4.16) sestavena charakteristická funkce diskrétního regulačního obvodu<br />

H<br />

n<br />

( z) = a z + ... + a1z<br />

a0<br />

n<br />

+<br />

(4.21)


Strana 50<br />

Použití frekvenčních charakteristik u diskrétních systémů<br />

U Michajlov-Leonhardova kritéria pro spojité obvody se do charakteristické funkce<br />

n<br />

H ( s) = an s + ... + a0<br />

spojitého obvodu za proměnnou s dosadí s = jω, tj. mez stability.<br />

Podobně je tomu u charakteristické funkce H(z) diskrétního obvodu, kdy se za proměnnou<br />

z příslušné hranice stability dosazuje jednotková kružnice, jak je uvedeno v (4.22)<br />

z = e<br />

sT<br />

= e<br />

jωT<br />

= cos ωT<br />

+<br />

j sin ωT<br />

(4.22)<br />

Touto substitucí lze dostat funkce bezrozměrné proměnné ωT<br />

H<br />

+<br />

nωT<br />

nωT<br />

j<br />

j<br />

( jωT<br />

) = ane<br />

+ ... + a1e<br />

+ a0<br />

= an<br />

( cos nωT<br />

+ j sin nωT<br />

) + ... + a1( cosωT<br />

+ j sin ωT<br />

)<br />

a = a cos nωT<br />

+ ... + a cosωT<br />

+ a + j( a sin nωT<br />

+ ... + a sin ωT<br />

)<br />

0<br />

n<br />

1<br />

0<br />

1444442<br />

444443<br />

Re<br />

( ωT<br />

)<br />

n<br />

1<br />

14444244443<br />

(4.23)<br />

ke které se sestrojí křivka H(jωT), podobně jako se sestrojovala frekvenční charakteristika<br />

k jakémukoliv frekvenčnímu přenosu spojitého systému.<br />

Skutečnost, zda je regulační obvod stabilní, lze zjistit podle definice uvedené v<br />

[Švarc, 2003]: Diskrétní obvod je stabilní, když křivka H(jωT) začíná na kladné reálné<br />

poloose a průvodič H opíše v kladném smyslu úhel πn (n.180°) (n je stupeň<br />

charakteristické rovnice) - (obr. 4.5).<br />

Im<br />

( ωT<br />

)<br />

+<br />

Im<br />

H(jωT)<br />

π 0<br />

Re<br />

Obr. 4.5


Použití frekvenčních charakteristik u diskrétních systémů Strana 51<br />

4.2.3 Nyquistovo kritérium<br />

Z frekvenčních kritérií se u spojitých systémů nejčastěji používá kritérium<br />

Nyquistovo. Také diskrétní regulační obvody se mohou pyšnit tímto kritériem, nicméně se<br />

diskrétní verze Nyquistova kritéria stability používání spíše zřídka, protože je k výpočtu a<br />

konstrukci křivek zapotřebí vhodný software. Ten ovšem není podmínkou, ale bez něj se<br />

stává výpočet příliš pracný. Proto zde řešení Nyquistova kritéria jen nastíním.<br />

Jedná se pochopitelně o frekvenční charakteristiku rozpojeného obvodu,<br />

sestrojovanou pro bezrozměrnou frekvenci ωT od 0 do π, protože pak je od π do 2π<br />

symetrická podle reálné osy a dál se její průběh periodicky opakuje.<br />

Pokud je uvažován jednoduchý diskrétní regulační obvod uvedený na obr. 4.6, pak<br />

pro určení stability je potřeba stanovit přenos sériového zapojení regulátor - tvarovač -<br />

regulovaná soustava, a to jako z-přenos G(z).<br />

w<br />

G R (z)<br />

T<br />

G<br />

T<br />

tvarovač<br />

( s)<br />

e<br />

= 1−<br />

s<br />

−sT<br />

G S (s)<br />

y<br />

T<br />

Obr. 4.6 - Diskrétní regulační obvod<br />

Po odvození uvedeném v [Davidová, 2002] je výsledný vztah pro z-přenos<br />

rozpojeného diskrétního regulačního obvodu<br />

−1<br />

⎧GS<br />

( ) ( ) ( ) ( ) ( )<br />

( s)<br />

⎫<br />

G0 z = GR<br />

z ⋅GS<br />

z = GR<br />

z ⋅ 1−<br />

z ⋅ Z⎨<br />

⎬<br />

(4.24)<br />

⎩ s ⎭<br />

Aby bylo možné sestrojit frekvenční charakteristiku nutnou k určení stability, je<br />

třeba za hodnotu z dosadit výraz (4.22), čímž se stanoví frekvenční přenos diskrétního<br />

systému G(jωT) a z něho se vykreslí frekvenční charakteristika buď v komplexní rovině<br />

nebo v logaritmických souřadnicích.<br />

Při určováni stability diskrétního regulačního obvodu Nyquistovým kritériem platí<br />

stejná pravidla, jako pro spojité regulační obvody, tedy o stabilitě rozhoduje poloha<br />

kritického bohu –1 vzhledem k této charakteristice.<br />

Pro představu, jak diskrétní verzi Nyquistova kritéria řešit, pro názornost uvedu v kapitole<br />

5 příklad.<br />

4.3 Syntéza regulačních obvodů<br />

Na rozdíl od spojitých obvodů, v diskrétních regulačních obvodech je zapotřebí brát<br />

v úvahu další důležitý parametr. Je to vzorkovací perioda T a její velikost se obvykle volí<br />

na základě dynamických vlastností regulované soustavy. Pro určení hodnoty vzorkovací<br />

periody se využívá některý z těchto empirických vztahů [Balátě, 2003]


Strana 52<br />

Použití frekvenčních charakteristik u diskrétních systémů<br />

T<br />

T<br />

T<br />

T1<br />

≈ ,<br />

10<br />

⎛ 1 1 ⎞<br />

≈ ⎜ ÷ ⎟T<br />

⎝15<br />

6 ⎠<br />

⎛ 1 1 ⎞<br />

≈ ⎜ ÷<br />

⎝ 4 2 ⎠<br />

⎟∑<br />

95<br />

,<br />

τ ,<br />

j<br />

(4.25)<br />

(4.26)<br />

(4.27)<br />

kde T 1 je největší časová konstanta regulované soustavy, T 95 je doba, za kterou přechodová<br />

charakteristika regulované soustavy dosáhne 95% ustálené hodnoty a ∑τ j<br />

je součet<br />

časových konstant regulované soustavy. Vztah (4.26) se používá pro nekmitavou<br />

proporcionální regulovanou soustavu.<br />

V případě velkého dopravního zpoždění regulované soustavy se vzorkovací perioda<br />

stanovuje podle vzorce [Isermann, 2003]<br />

⎛ 1 1 ⎞<br />

T ≈ ⎜ ÷ ⎟ Td ,<br />

⎝ 8 3 ⎠<br />

nebo vztahu [Vítečková, 1998]<br />

T < 0,32Td<br />

,<br />

(4.28)<br />

(4.29)<br />

kde T d je hodnota dopravního zpoždění.<br />

Za vhodnou se považuje taková velikost vzorkovací periody, při které nedojde ke<br />

zhoršení kvality regulace o více než 15% ve srovnání s použitým analogickým spojitým<br />

regulátorem. Pokud je vzorkovací perioda příliš malá, zvyšuje se nárok na rychlost nejen<br />

číslicového regulátoru, ale rovněž převodníků a měřícího a akčního členu. Naopak<br />

zvětšováním vzorkovací periody dochází ke ztrátě informace v regulované veličině mezi<br />

jednotlivými okamžiky vzorkování, což vede k destabilizaci regulačního obvodu [Balátě,<br />

2003].<br />

Při návrhu diskrétního regulačního obvodu uvažujeme jednoduchý obvod, uvedený<br />

na obr. 4.6, který obsahuje spojitou regulovanou soustavu a před ní umístěný tvarovač.<br />

Spojitá regulovaná soustava je popsána z-přenosem, jež bude v dalším textu<br />

označován G S (z) a určí se podle vztahu uvedeném např. v [Balátě, 2003] a také v kapitole<br />

4.2.3 a po zjednodušení je sejný jako vztah (4.24).<br />

4.3.1 Typizovaná logaritmická amplitudová frekvenční charakteristika<br />

Metodu typizované logaritmické amplitudové frekvenční charakteristiky<br />

(TLAFCH) pro spojité regulační obvody jsem zpracoval v kapitole 3.3.2.3. Zde jsem uvedl<br />

návrh parametrů spojitých regulátorů. V této kapitole se pokusím zpracovat tuto metodu<br />

pro diskrétní regulační obvod. Tato metoda byla takto zpracována a publikována pouze v<br />

[Davidová, 2004].


Použití frekvenčních charakteristik u diskrétních systémů Strana 53<br />

Princip metody (TLAFCH) pro diskrétní regulační obvody je stejný jako v kapitole<br />

3.3.2.3, kde je tato metoda zpracována pro spojité regulační obvody.<br />

Je-li frekvenční přenos rozpojeného regulačního obvodu určen vztahem<br />

G<br />

0<br />

( jΩ) = G ( jΩ) G ( jΩ)<br />

R<br />

S<br />

, (4.30)<br />

pak přenos regulátoru je dán<br />

( jΩ)<br />

G<br />

=<br />

G<br />

( jΩ)<br />

( jΩ)<br />

což se v logaritmických souřadnicích vyjádří vztahem<br />

G<br />

R<br />

0<br />

S<br />

, (4.31)<br />

G<br />

R<br />

( jΩ) = G ( jΩ) − GS<br />

( jΩ) dB 0 dB<br />

dB<br />

. (4.32)<br />

Logaritmická amplitudová frekvenční charakteristika se tedy stanoví z rozdílů<br />

těchto charakteristik pro rozpojený regulační obvod a regulovanou soustavu. Typ a<br />

parametry regulátoru lze určit podle frekvencí zlomů a sklonů asymptot frekvenčního<br />

přenosu G R (jΩ), který se převede na z-přenos číslicového korekčního členu.<br />

Při sestrojování amplitudové frekvenční charakteristiky v logaritmických<br />

souřadnicích vycházíme z frekvenčního přenosu G S (jΩ). Pro určení frekvenčního přenosu<br />

je možné použít metodiku, která vychází z toho, že frekvenční charakteristiku lze rozdělit<br />

na dvě oblasti. První je charakteristika v oblast nízkých frekvencí a druhá v oblasti<br />

vysokých frekvencí.<br />

V oblasti nízkých frekvencí je dle odvozených vzorců uvedených v [Davidová,<br />

2004] zřejmé, že logaritmická amplitudová frekvenční charakteristika spojité regulované<br />

soustavy se v této oblasti shoduje s logaritmickou amplitudovou frekvenční charakteristiku<br />

spojité regulované soustavy s tvarovačem na vstupu. Proto je možné použít metodiku<br />

sestrojování frekvenční charakteristiky, která byla použita pro spojité systémy.<br />

Začátek logaritmické frekvenční charakteristiky v oblasti vysokých frekvencí<br />

splývá s koncem této charakteristiky v oblasti vysokých frekvencí.<br />

Nyní zde popíšu sestrojení typizované logaritmické frekvenční charakteristiky<br />

rozpojeného regulačního obvodu. Stejně jako u spojitých systémů platí, že průběh<br />

nízkofrekvenční části charakteristiky určuje přesnost regulace v ustáleném stavu (velikosti<br />

trvalé regulační odchylky) a středofrekvenční část charakteristiky určuje dynamické<br />

vlastnosti uzavřeného regulačního obvodu. Na obr. 4.7 je uveden příklad tvaru<br />

charakteristiky se zakreslenými významnými lomovými frekvencemi. Snahou je, aby<br />

zvolený postup konstrukce této charakteristiky byl co možná nejjednodušší. Nejprve tady<br />

vypočítám frekvenci řezu Ω<br />

ř<br />

, při které asymptota protíná osu 0 dB. Tu stanovím na<br />

základě předpokládané doby regulace t r podle vztahu<br />

π<br />

Ω<br />

. (4.33)<br />

ř<br />

= ( 1 ÷ 4)<br />

t<br />

r<br />

Tímto bodem vedu asymptotu se sklonem –20 dB/dek, která tvoří středofrekvenční oblast<br />

typizované charakteristiky. K určení frekvencí ohraničujících středofrekvenční oblast<br />

použiji vztahy pro spojité systémy aplikované na diskrétní systémy


Strana 54<br />

Použití frekvenčních charakteristik u diskrétních systémů<br />

Ω<br />

Ω<br />

2<br />

≈<br />

( 2 ÷ )<br />

3<br />

4<br />

Ω<br />

≈<br />

Ω<br />

2<br />

ř<br />

3<br />

Ω<br />

ř<br />

(4.34)<br />

(4.35)<br />

A<br />

[dB]<br />

0 dB/dek<br />

–40 dB/dek<br />

–20 dB/dek<br />

Ω 3<br />

2 /( T − 2T∑<br />

)<br />

Ω 1<br />

Ω 2<br />

Ω ř<br />

Ω<br />

–40 dB/dek<br />

–60 dB/dek<br />

Obr. 4.7 - Typizovaná logaritmická amplitudová frekvenční charakteristika<br />

Podle vztahu (4.34) určím lomovou frekvenci Ω 3 z rozmezí hodnot, což umožňuje volit<br />

fázovou bezpečnost, neboť ta je na šířce středofrekvenční oblasti závislá. Čím se volí větší<br />

hodnota koeficientu z intervalu od 2 do 4, tím je větší také fázová bezpečnost. Podle<br />

vztahu (4.38) se stanoví lomová frekvence Ω 2 . Jelikož výraz neobsahuje rovnítko, i zde je<br />

možná případná úprava hodnoty této frekvence. Pokud totiž některá z převrácených hodnot<br />

časových konstant regulované soustavy se této vypočítané hodnotě blíží, je lepší použít<br />

tuto. Totéž platí i pro lomovou frekvenci Ω 3 . Je to z důvodu toho, že čím více se<br />

typizovaná frekvenční charakteristika rozpojeného regulačního obvodu svým tvarem blíží<br />

k tvaru charakteristiky regulované soustavy, tím je navrhovaný korekční člen jednodušší.<br />

Z bodu tvořícího začátek středofrekvenční oblasti se sestrojí asymptota o sklonu –<br />

40 dB/dek, případně –60 dB/dek. Její průsečík s frekvenční charakteristikou regulované<br />

soustavy určí lomovou frekvenci Ω 1 , která tvoří konec nízkofrekvenční oblasti. Sklon<br />

asymptoty v nízkofrekvenční oblasti se volí na základě požadovaného řádu integrace<br />

rozpojeného regulačního obvodu. Jestliže je regulovaná soustava proporcionální a chceme<br />

dosáhnout regulace s nulovou regulační odchylkou, pak je třeba volit sklon charakteristiky<br />

v nízkofrekvenční oblasti –20 dB/dek. Je-li však regulovaná soustava integrační s řádem<br />

integrace 1, případně 2, pak je výhodné, aby nízkofrekvenční část typizované<br />

charakteristiky splývala s počátkem charakteristiky regulované soustavy a měla tedy sklon<br />

–20 dB/dek, resp. –40 dB/dek.


Použití frekvenčních charakteristik u diskrétních systémů Strana 55<br />

Sklon asymptot ve vysokofrekvenční oblasti se ukázalo vhodné volit stejný jako u<br />

asymptot frekvenčních charakteristik regulované soustavy.<br />

Při určování parametrů číslicového regulátoru, popř. korekčního členu, se provede<br />

grafické sestrojení logaritmické amplitudové frekvenční charakteristiky podle vztahu<br />

(4.32), tedy jako rozdíl amplitudové frekvenční charakteristiky rozpojeného obvodu<br />

(typizované) a regulované soustavy, jak je vyobrazeno na obr. 4.8. Na základě tvaru této<br />

charakteristiky, podle frekvencí zlomů a sklonů asymptot, se stanoví frekvenční přenos<br />

G R (jΩ).<br />

A<br />

[dB]<br />

G S (jω)<br />

Ω1<br />

Ω 2<br />

Ω3<br />

Ω ř<br />

Ω 4<br />

Ω5<br />

Ω6<br />

Ω<br />

G R (jω)<br />

G 0 (jω)<br />

Obr. 4.8 - Sestrojení logaritmické amplitudové frekvenční charakteristiky<br />

číslicového regulátoru<br />

4.3.2 Metoda transformovaných frekvenčních charakteristik<br />

Mezi metody syntézy diskrétních regulačních obvodů, které jsou vhodné i pro<br />

průmyslovou praxi, patří metoda transformovaných frekvenčních charakteristik uvedená v<br />

[Černý, 1984], [Zeman, 1987]. Pro regulační obvod na obr. 4.6 se na základě vztahu (4.32)<br />

určí z-přenos spojité regulované soustavy G S (z). Vzájemný vztah mezi operátorem s a<br />

operátorem z je dán rovnicí<br />

T<br />

1+<br />

s<br />

sT<br />

z = e ≈ 2<br />

(4.36)<br />

T<br />

1−<br />

s<br />

2


Strana 56<br />

Použití frekvenčních charakteristik u diskrétních systémů<br />

Vzhledem k tomu, že frekvenční charakteristika soustavy je transcendentní funkcí<br />

proměnné jω, je znesnadněno přímé sestrojení frekvenční charakteristiky diskrétního<br />

systému pouhým dosazením s=jω. Proto se zavádí operátor q, který je definován<br />

transformací<br />

T<br />

1+<br />

q<br />

z = 2<br />

T<br />

1−<br />

q<br />

2<br />

2<br />

q =<br />

T<br />

. (4.37a,b)<br />

Vztahy (4.36) a (4.37) představují vzájemné zobrazení komplexních rovin Z, S a Q.<br />

Pomocí těchto transformací se levá komplexní polorovina roviny S převede do vnitřní části<br />

jednotkové kružnice roviny Z. Body s=jω, které leží na imaginární ose roviny S, se<br />

transformují na body q=jΩ, které leží rovněž na imaginární ose, ovšem roviny Q.<br />

Transformovaný přenos v proměnné q se stanoví ze vztahu (4.32) dosazením<br />

transformace (4.37a) za hodnotu<br />

z<br />

z<br />

−1<br />

+ 1<br />

G<br />

S<br />

−1<br />

⎧GS<br />

( ) ( )<br />

( s)<br />

q = 1−<br />

z ⋅ Z<br />

⎨<br />

⎩<br />

s<br />

⎫<br />

⎬<br />

⎭<br />

T<br />

1+<br />

q<br />

z=<br />

2<br />

T<br />

1−<br />

q<br />

2<br />

(4.38)<br />

Při praktickém výpočtu se přenos G S (s) rozloží na parciální zlomky a potom lze pro<br />

výpočet transformovaného přenosu G S (q) použít tabulku 4.1 převzatou z [Černý, 1984]. V<br />

ní jsou uvedeny základní typy přenosů a jim příslušející jednotlivé transformované<br />

přenosy.<br />

Transformovaný frekvenční přenos se určí dosazením<br />

q = jΩ<br />

(4.39)<br />

do transformovaného přenosu G S (q). Z něho se sestrojí amplitudová, případně fázová<br />

transformovaná logaritmická frekvenční charakteristika stejným postupem jako v případě<br />

spojitých systémů.


Použití frekvenčních charakteristik u diskrétních systémů Strana 57<br />

Tabulka 4.1 - Transformované přenosy<br />

Č. G S (s) G S (q)<br />

1<br />

1 s<br />

T<br />

1− q<br />

2<br />

q<br />

1<br />

2 2<br />

s<br />

T<br />

1−<br />

q<br />

2<br />

2<br />

q<br />

1<br />

3 1+<br />

s<br />

T 1<br />

T<br />

1−<br />

q<br />

2<br />

*<br />

1+<br />

T q<br />

1<br />

1<br />

4 ( ) 2<br />

1+ T s<br />

1<br />

⎛ T ⎞⎡<br />

⎛ T<br />

⎜1<br />

− q⎟⎢1<br />

+<br />

2<br />

⎜<br />

⎝ ⎠⎢⎣<br />

⎝ 2T<br />

*<br />

( 1+<br />

T q) 2<br />

1<br />

1<br />

T<br />

2T<br />

*<br />

1<br />

⎞<br />

− β<br />

⎟T<br />

⎠<br />

*<br />

⎤<br />

⎥<br />

⎥⎦<br />

1<br />

5 ( 1+<br />

T s)<br />

s<br />

1<br />

1<br />

6 2 ( 1+<br />

T s)<br />

s<br />

1<br />

1<br />

7 ( 1+<br />

s)( + T s)<br />

T1 1<br />

2<br />

⎛ T ⎞<br />

⎜1<br />

− q⎟<br />

⎝ 2 ⎠<br />

q<br />

⎛ T ⎞<br />

⎜1<br />

− q⎟<br />

⎝ 2 ⎠<br />

q<br />

( 1+<br />

T β q)<br />

*<br />

( 1+<br />

T q)<br />

⎛ T ⎞⎛<br />

T1T<br />

2<br />

⎜1<br />

− q⎟<br />

1<br />

2<br />

⎜ +<br />

⎝ ⎠⎝<br />

T1<br />

− T<br />

1<br />

2 2<br />

( 1+<br />

T1<br />

β q − T1<br />

β q )<br />

2 *<br />

( 1+<br />

T q)<br />

*<br />

*<br />

( 1+<br />

T q)( 1+<br />

T q)<br />

1<br />

1<br />

2<br />

1<br />

( β − β )<br />

2<br />

2<br />

1<br />

⎞<br />

q<br />

⎟<br />

⎠<br />

1<br />

8 ( 1+ T s)( + T s)<br />

⎜1<br />

q⎟⎢1<br />

+ ( T1β<br />

1<br />

+ T2β<br />

2<br />

)<br />

s<br />

1<br />

1<br />

2<br />

Poznámka:<br />

T<br />

⎛ T −<br />

⎝ 2<br />

*<br />

i<br />

⎞⎡<br />

⎠⎣<br />

i<br />

q<br />

⎛ T1<br />

β1<br />

− T2β<br />

2<br />

⎞<br />

q +<br />

⎜ β1β<br />

2<br />

+<br />

T1T<br />

2q<br />

T1<br />

T<br />

⎟<br />

⎝<br />

−<br />

2 ⎠<br />

*<br />

*<br />

( 1+<br />

T q)( 1+<br />

T q)<br />

T T<br />

= h<br />

2 2T<br />

= T<br />

β<br />

i − 1<br />

Ti* <br />

1<br />

2<br />

⎤<br />

⎥<br />

⎦<br />

2


Strana 58<br />

Použití frekvenčních charakteristik u diskrétních systémů<br />

Pokud je splněna podmínka<br />

ωT<br />

2 ≤ 0,25<br />

pak je ω ≈ Ω a současně platí přibližný vztah<br />

, (4.40)<br />

G<br />

S<br />

( jΩ) ≈ G ( jω)<br />

S<br />

. (4.41)<br />

To znamená, že transformovaná logaritmická amplitudová frekvenční charakteristika<br />

(LAFCH)souhlasí přibližně s touto charakteristikou spojitého systému bez vzorkování při<br />

dané době vzorkování T při frekvencích<br />

2<br />

Ω ≤ 0,25 . (4.42)<br />

T<br />

Při určování parametrů číslicového regulátoru se vychází z jeho LAFCH, která se stanoví<br />

na základě LAFCH rozpojeného regulačního obvodu G 0 (jΩ), jež musí mít požadovaný<br />

průběh v nízkofrekvenční a středofrekvenční oblasti (podobně jako u spojitých systémů) a<br />

LAFCH regulované soustavy G S (jΩ) podle vztahu (4.32). Podle tvaru frekvenční<br />

charakteristiky regulátoru se určí transformovaný přenos regulátoru G R (q), do něhož se<br />

dosadí transformační vztah (4.40b) a získá se tak z-přenos regulátoru<br />

G<br />

R<br />

( z) = G ( z)<br />

R<br />

2 z−1<br />

q=<br />

T z+<br />

1<br />

(4.43)<br />

Z LAFCH rozpojeného regulačního obvodu lze odečíst frekvenci řezu Ω ř , z které je možné<br />

stanovit dobu regulace t r vztahem<br />

t r<br />

≈ , 5<br />

(4.44)<br />

Ω<br />

1 ř<br />

γ<br />

Velikost překmitu σ na přechodové charakteristice, udávaná v procentech, je dána fázovou<br />

bezpečností γ podle rovnice<br />

.<br />

σ = 70 −<br />

(4.45)


Strana 59<br />

5 POUŽITÍ VYBRANÝCH METOD NA PŘÍKLADECH<br />

V této kapitole uvedu příklady použití frekvenčních charakteristik při analýze a<br />

syntéze regulačních obvodů. Ukážu zde i příklady, na kterých budu demonstrovat diskrétní<br />

verze metod, které jsem v předchozí kapitole uvedl, kvůli jejich náročnosti, jen stručně.<br />

Uvedu zde tyto metody:<br />

• Michajlov-Leonhardovo kritérium pro spojité regulační obvody<br />

• Nyquistovo kritérium pro spojité regulační obvody<br />

• Michajlov-Leonhardovo kritérium pro diskrétní regulační obvody<br />

• Nyquistovo kritérium pro diskrétní regulační obvody<br />

• Metoda transformovaných frekvenčních charakteristik<br />

Z důvodu omezeného rozsahu diplomové práce, mohu uvést příklady pouze v<br />

omezeném množství. Rovnice a grafy byly zpracovány v programu Microsoft Office Excel<br />

2003, tyto poznámky jsou na přiloženém CD uvedeny v souboru (Výpočty a grafy.doc).<br />

5.1 Příklady pro spojité regulační obvody<br />

Zde uvedu dvě kritéria na analýzu regulačního obvodu, a to Michajlov-<br />

Leonhardovo a Nyquistovo.<br />

5.1.1 Michajlov-Leonhardovo kritérium<br />

Příklad 5.1: Podle obr. 5.1 určete stabilitu regulačního obvodu podle pomocí Michajlov-<br />

Leonhardova kritéria.<br />

v<br />

G S<br />

( s)<br />

=<br />

s<br />

1<br />

( 0,1s<br />

+ 1)( 0,5s<br />

+ 1)<br />

1<br />

G R<br />

1<br />

⎝ 0,5s<br />

⎛<br />

⎞<br />

( s) = 40 ⎜1+<br />

+ 0, s⎟ ⎠<br />

Obr. 5.1<br />

y<br />

w<br />

Řešení:<br />

Nejprve určím přenos rozpojeného obvodu, který ze známého vztahu<br />

( s) = G ( s) ⋅G<br />

( s)<br />

G0 S R<br />

je<br />

G<br />

0<br />

( s)<br />

=<br />

s<br />

1<br />

⎛<br />

40⎜1<br />

+<br />

⎝<br />

1<br />

0,5<br />

⎞<br />

+ 0,1s<br />

⎟ =<br />

⎠<br />

4s<br />

4 3 2<br />

( 0,1s<br />

+ 1)( 0,5s<br />

+ 1) s 0,05s<br />

+ 0,6s<br />

+ s<br />

2<br />

+ 40s<br />

+ 80


Strana 60<br />

Použití vybraných metod na příkladech<br />

Dále vyjádřím charakteristickou rovnici<br />

0,05s<br />

Michajlův-Leonhardův vektor je<br />

Tabulka 5.1<br />

H<br />

ω Re Im<br />

4<br />

3 2<br />

+ 0,6s<br />

+ 5s<br />

+ 40s<br />

+ 80 = 0<br />

4<br />

3<br />

2<br />

( jω) = 0,05( jω) + 0,6( jω) + 5( jω) + 40( jω)<br />

4 2<br />

= 0,05<br />

14<br />

ω<br />

44<br />

−<br />

2<br />

5<br />

4<br />

ω<br />

44<br />

+<br />

3<br />

80 +<br />

Re<br />

0 80 0<br />

0,5 78,75 19,93<br />

1 75,05 39,4<br />

2 60,8 75,2<br />

3 39,05 103,8<br />

5 -13,75 125<br />

7 -44,95 74,2<br />

8 -35,2 12,8<br />

9 3,05 -77,4<br />

9,5 36,00 -134,43<br />

10 80 -200<br />

10,2 101 -228,73<br />

3<br />

( 40ω<br />

− 0,6ω<br />

)<br />

j<br />

144<br />

2443<br />

Im<br />

+ 80 =<br />

K tomuto výrazu sestrojím na základě<br />

tabulky 5.1 Michajlov.Leonhardovu křivku H(jω)<br />

- obr. 5.2. Protože stupeň charakteristické rovnice<br />

je n=4 a křivka prochází v kladném smyslu<br />

čtyřmi za sebou jdoucími kvadranty, je regulační<br />

obvod stabilní.<br />

150<br />

Im<br />

100<br />

50<br />

0<br />

-60 -40 -20 0 20 40 60 80 100 120<br />

-50<br />

-100<br />

Re<br />

-150<br />

-200<br />

-250<br />

Obr. 5.2 - Michajlov-Leonhardova křivka H(jω)


Použití vybraných metod na příkladech Strana 61<br />

5.1.2 Nyquistovo kritérium<br />

Příklad 5.2: Určete stabilitu regulačního obvodu podle obr. 5.3 pomocí Nyquistova<br />

kritéria.<br />

v<br />

G S<br />

( s)<br />

=<br />

1<br />

G R<br />

( s)<br />

=<br />

s<br />

1<br />

( s + 1)( 12s<br />

+ 1)<br />

Obr. 5.3<br />

y<br />

w<br />

Řešení: Nejprve určím, stejně jako v 5.1.1, přenos rozpojeného obvodu, který ze<br />

známého vztahu G0 ( s) = GS<br />

( s) ⋅GR<br />

( s)<br />

je<br />

G<br />

0<br />

G<br />

0<br />

= 1<br />

1<br />

=<br />

3<br />

s( s + 1)( 10s<br />

+ 1) 10s<br />

+ 11s<br />

2 + s<br />

1<br />

1<br />

jω<br />

=<br />

=<br />

3<br />

2<br />

3 2<br />

10 jω<br />

+ 11 jω<br />

+ jω<br />

−10<br />

jω<br />

−11ω<br />

+<br />

( s)<br />

( )<br />

( ) ( )<br />

3 2<br />

−10<br />

jω<br />

+ 11ω<br />

+ jω<br />

=<br />

2 6<br />

5 2 4<br />

5<br />

4<br />

2 2 4<br />

2<br />

100 j ω −110<br />

jω<br />

−10<br />

j ω + 110 jω<br />

−121ω<br />

−11<br />

jω<br />

−10<br />

j ω + 11 jω<br />

+<br />

3 2<br />

−10<br />

jω<br />

+ 11ω<br />

+ jω<br />

=<br />

2 6 2 4<br />

4 2 4<br />

100 j ω −10<br />

j ω −121ω<br />

−10<br />

j ω +<br />

3<br />

( ω −10ω<br />

)<br />

3 2<br />

2<br />

−10<br />

jω<br />

+ 11ω<br />

+ jω<br />

11ω<br />

=<br />

=<br />

+ j<br />

6<br />

4 2<br />

6<br />

4 2<br />

6<br />

4 2<br />

−100ω<br />

−101ω<br />

− ω 14 −10044<br />

ω −2<br />

101 444<br />

ω − 3 ω 14 −100<br />

4ω<br />

44 −2<br />

101 4444<br />

ω − 3 ω<br />

Re<br />

2 2<br />

j ω<br />

3 2<br />

−10<br />

jω<br />

+ 11ω<br />

+<br />

⋅<br />

3 2<br />

jω<br />

−10<br />

jω<br />

+ 11ω<br />

+<br />

3 2<br />

−10<br />

jω<br />

+ 11ω<br />

+ jω<br />

=<br />

6 4<br />

4 4 2<br />

−100ω<br />

+ 10ω<br />

−121ω<br />

+ 10ω<br />

− ω<br />

Im<br />

jω<br />

=<br />

jω<br />

=<br />

2 2<br />

j ω<br />

=<br />

Tabulka 5.2<br />

ω Re Im<br />

0,2 -2,12 -0,58<br />

0,22 -1,80 -0,38<br />

0,24 -1,54 -0,25<br />

0,26 -1,33 -0,15<br />

0,28 -1,15 -0,08<br />

0,3 -1,01 -0,03<br />

0,4 -0,56 0,08<br />

0,5 -0,34 0,09<br />

0,6 -0,22 0,09<br />

1 -0,05 0,04<br />

Kořeny jmenovatele jsou 0, -0,1, -1. Žádný z<br />

nich není kladný (neleží v pravé komplexní<br />

polorovině), rozpojený obvod je tedy stabilní, a<br />

proto je možno Nyquistova kritéria použít.<br />

Frekvenční přenos G 0 (jω) rozdělím na reálnou a<br />

imaginární část a sestrojím frekvenční<br />

charakteristiku rozpojeného obvodu v<br />

komplexní rovině (tab. 5.2, obr. 5.4). Kritický<br />

bod [–1, 0] leží vlevo od frekvenční<br />

charakteristiky G 0 (jω)a proto je obvod stabilní.


Strana 62<br />

Použití vybraných metod na příkladech<br />

0,20<br />

Im<br />

Re<br />

0,10<br />

0,00<br />

-2,50 -2,00 -1,50 -1,00 -0,50 0,00<br />

-0,10<br />

-0,20<br />

-0,30<br />

-0,40<br />

-0,50<br />

Obr. 5.4<br />

-0,60<br />

-0,70<br />

5.2 Příklady pro diskrétní regulační obvody<br />

Zde opět uvedu dvě kritéria na analýzu regulačního obvodu, a to Michajlov-<br />

Leonhardovo a Nyquistovo a také příklad na syntézu regulačního obvodu, tedy metodu<br />

typizované logaritmické frekvenční charakteristiky.<br />

5.2.1 Michajlov-Leonhardovo kritérium<br />

Příklad 5.3: Na základě charakteristické rovnice určete stabilitu diskrétního regulačního<br />

obvodu:<br />

3 2<br />

z − 0,3z + 0,7z − 0,5 = 0<br />

Řešení: Do charakteristické funkce obvodu<br />

H<br />

( jωT)<br />

H z<br />

dosadím za z vztah<br />

sT<br />

z = e<br />

= e<br />

3jωT<br />

− 0,3e<br />

2 jωT<br />

= cos3ωT<br />

+ jsin 3ωT<br />

− 0,3<br />

3 2<br />

( ) = z − 0,3z + 0,7z − 0, 5<br />

= e<br />

jωT<br />

= cos3<br />

14444444<br />

ωT<br />

− 0,3cos 2<br />

2<br />

ωT<br />

4<br />

+<br />

444444 3<br />

Re<br />

+ 0,7e<br />

jωT<br />

= cos ωT<br />

+<br />

− 0,5 =<br />

j sinωT<br />

( cos 2ωT<br />

+ jsin 2ωT) + 0,7( cosωT<br />

+ jsinωT)<br />

− 0,5 =<br />

0,7 cosωT<br />

− 0,5 + j( sin 3ωT<br />

− 0,3sin 2ωT<br />

+ 0,7sinωT)<br />

1444444<br />

2444444<br />

3<br />

Im


Použití vybraných metod na příkladech Strana 63<br />

Tabulka 5.3<br />

ωΤ Re Im<br />

0,00 1,00 0,00<br />

0,52 0,06 1,09<br />

1,05 -0,90 0,35<br />

1,57 -0,10 -0,30<br />

1,83 0,39 0,12<br />

2,09 0,40 0,87<br />

2,36 -0,19 1,50<br />

2,62 -1,16 1,61<br />

3,14 -2,40 0,00<br />

Pokračuji výpočtem, který je uveden v<br />

tab. 5.3 a následně vykreslením grafu (obr. 5.5).<br />

Podle definice stability regulačního obvodu,<br />

uvedené v kapitole 4.2.2, se dle obr. 5.5 jedná o<br />

obvod stabilní, protože křivka H(jωT) začíná na<br />

kladné reálné poloose a průvodič H opsal v<br />

kladném smyslu úhel πn (n.180°). Vzhledem k<br />

tomu, že stupeň charakteristické rovnice n=3, pak<br />

3.180° = 540°, což podle grafu na obr. 5.5<br />

souhlasí.<br />

2,00<br />

Im<br />

1,50<br />

1,00<br />

0,50<br />

-3,00 -2,50 -2,00 -1,50 -1,00 -0,50 0,00 0,50 1,00 1,50<br />

0,00<br />

Re<br />

-0,50<br />

Obr. 5.5


Strana 64<br />

Použití vybraných metod na příkladech<br />

5.2.2 Nyquistovo kritérium<br />

Příklad 5.4: Pomocí Nyquistova kritéria vyšetřete stabilitu diskrétního regulačního<br />

obvodu uvedeném na obr. 5.6.<br />

w(t)<br />

T=1[s]<br />

1−<br />

e<br />

s<br />

−sT<br />

1<br />

( 2s + 1)( s + 1)<br />

T=1[s]<br />

y(k)<br />

y(t)<br />

Obr. 5.6<br />

Řešení:<br />

Základem je nejprve vyjádření přenosu řízení daného rovnicí<br />

G<br />

w<br />

( z)<br />

( )<br />

( )<br />

Y z<br />

=<br />

W z<br />

G<br />

R<br />

=<br />

1+<br />

G<br />

( z) G<br />

S<br />

( z)<br />

( z) G ( z)<br />

R<br />

S<br />

Vzhledem k tomu, že G R<br />

( z) = 1, stačí na základě vztahu (4.24) vyjádřit ( z)<br />

dosadit ho do G w<br />

( z)<br />

.<br />

G<br />

S<br />

−1<br />

( z) = ( 1−<br />

z )<br />

⎧<br />

Z⎨<br />

⎩s<br />

1<br />

( 2s + 1)( s + 1)<br />

⎫<br />

⎬ =<br />

⎭<br />

−1<br />

( 1−<br />

z )<br />

⎧1<br />

2 1 ⎫<br />

Z⎨<br />

− + ⎬ =<br />

⎩s<br />

s + 0,5 s + 1⎭<br />

G S<br />

a<br />

z −1<br />

⎡ z 2z z ⎤ ⎡ 1 2 1 ⎤<br />

=<br />

( z 1)<br />

z ⎢<br />

− + = −<br />

0,5T<br />

T<br />

z 1 z e z e<br />

⎢ − + =<br />

−<br />

− ⎥ z 1 z 0,607 z 0,368<br />

⎥<br />

⎣ − − − ⎦ ⎣ − − − ⎦<br />

0,154z + 0,094<br />

=<br />

2<br />

z − 0,975z + 0,223<br />

K získání Z-obrazů byl použit slovník Z-transformace uvedený v [Švarc, 2002]. Přenos<br />

řízení je tedy<br />

0,154z + 0,094<br />

2<br />

z − 0,975z + 0,223 0,154z + 0,094<br />

G<br />

w<br />

( z)<br />

=<br />

=<br />

2<br />

0,154z + 0,094 z − 0,821z + 0,317<br />

1+<br />

z2 − 0,975z + 0,223<br />

Teď využiji vztahu (4.22) pro vyjádření G( jω<br />

T)


Použití vybraných metod na příkladech Strana 65<br />

G<br />

=<br />

( jωT)<br />

=<br />

e<br />

2 jωT<br />

0,154( cosωT<br />

+ jsinωT)<br />

+ 0,094<br />

( cos 2ωT<br />

+ jsin 2ωT) − 0,975( cosωT<br />

+ jsinωT)<br />

= ....................... =<br />

jωT<br />

0,154e + 0,094<br />

=<br />

jωT<br />

− 0,975e + 0,223<br />

Re<br />

644444444444444<br />

744444444444444<br />

8<br />

⎧( cos 2ωT<br />

− 0,975cosωT<br />

+ 0,223) ⋅ ( 0,154cosωT<br />

+ 0,094) + ( 0,154sinωT)<br />

⎫<br />

2<br />

2<br />

⎪<br />

⋅<br />

( cos 2ωT<br />

− 0,975cosωT<br />

+ 0,223) + ( sin 2ωT<br />

− 0,975sinωT)<br />

⎪<br />

= ⎨<br />

⎬ +<br />

⎪ ( sin 2ωT<br />

− 0,975sinωT)<br />

⋅<br />

⎪<br />

⎪⎩<br />

1<br />

⎪⎭<br />

( sin 2ωT<br />

− 0,975sinωT) ⋅ ( − 0,154cosωT<br />

− 0,097) + ( 0,154sinωT)<br />

2<br />

2<br />

( cos 2ωT<br />

− 0,975cosωT<br />

+ 0,223) + ( sin 2ωT<br />

− 0,975sinωT)<br />

( cos 2ωT<br />

− 0,975cosωT<br />

+ 0,223)<br />

⎧<br />

⎫<br />

j<br />

⎪<br />

⋅<br />

⎪<br />

+ ⎨<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⋅<br />

⎩14 444444444444<br />

1<br />

⎪<br />

24444444444444<br />

3⎭<br />

Im<br />

=<br />

+ 0,223<br />

Tabulka 5.4<br />

ωΤ Re Im<br />

0 1,00 0,000<br />

0,1 0,94 -0,333<br />

0,2 0,79 -0,571<br />

0,3 0,60 -0,681<br />

0,4 0,42 -0,684<br />

0,7 0,09 -0,439<br />

1 -0,02 -0,207<br />

1,2 -0,04 -0,111<br />

1,5 -0,05 -0,030<br />

2 -0,04 0,017<br />

2,5 -0,03 0,018<br />

3 -0,03 0,005<br />

3,14 -0,03 0,000<br />

Vzhledem k velmi pracnému výpočtu<br />

zde uvádím pouze výsledek, do kterého jsou<br />

dosazeny hodnoty ωT (viz. tab. 5.4).<br />

Kritický bod [–1, 0] leží vlevo od této<br />

charakteristiky a proto je obvod stabilní.


Strana 66<br />

Použití vybraných metod na příkladech<br />

0,100<br />

Im<br />

0,000<br />

Re<br />

-0,20 0,00 0,20 0,40 0,60 0,80 1,00 1,20<br />

-0,100<br />

-0,200<br />

-0,300<br />

-0,400<br />

-0,500<br />

-0,600<br />

-0,700<br />

-0,800<br />

Obr. 5.7<br />

5.2.3 Metoda transformovaných frekvenčních charakteristik<br />

Příklad 5.5: Navrhněte parametry číslicového metodou transformovaných frekvenčních<br />

charakteristik pro regulovanou soustavu<br />

1<br />

G S<br />

( s)<br />

=<br />

1+<br />

4s 1+<br />

0,5s<br />

( )( )<br />

je-li vzorkovací perioda T = 1s.<br />

V tomto příkladu využiji tabulku 4.1 uvedenou v kapitole 4.3.2. Uvedený přenos<br />

odpovídá v této tabulce výrazu č. 7. Časové konstanty regulované soustavy jsou T 1 = 4s a<br />

T 2 = 0,5s a její zesílení K = 1. Nejprve určím jednotlivé dílčí hodnoty<br />

* T T 1 1 1<br />

T1 ≈ cot gh = cot gh = ⋅8,042<br />

≈ 4,021s<br />

2 2T1<br />

2 2 ⋅ 4 2<br />

* T T 1 1 1<br />

T2 ≈ cot gh = cot gh = ⋅1,313<br />

≈ 0,675s<br />

2 2T 2 2 ⋅ 0,5 2<br />

*<br />

T1<br />

4,021<br />

β<br />

1<br />

= −1<br />

= −1<br />

= 0,00525<br />

T 4<br />

1<br />

*<br />

T2<br />

0,657<br />

β<br />

2<br />

= −1<br />

= −1<br />

= 0,314<br />

T 0,5<br />

2<br />

2


Použití vybraných metod na příkladech Strana 67<br />

T1T<br />

2<br />

=<br />

2<br />

β<br />

T − T<br />

4 ⋅ 0,5<br />

4 − 0,5<br />

( β − ) ≈ ( 0,314 − 0,00525) 0,176s<br />

T3 1<br />

≈<br />

1 2<br />

Dosadím tyto číselné hodnoty a stanovím tím transformovaný přenos regulované soustavy<br />

( q)<br />

( 1−<br />

0,5q)( 1+<br />

0,176q)<br />

( 1+<br />

4,021q )( 1 0,657q)<br />

G S<br />

= .<br />

+<br />

Z něho potom za pomoci vztahu (4.42) určím transformovaný frekvenční přenos<br />

( jΩ<br />

)<br />

( 1−<br />

0,5jΩ<br />

)( 1+<br />

0,176jΩ<br />

)<br />

( 1+<br />

4,021jΩ<br />

)( 1 0,657jΩ<br />

)<br />

G S<br />

= .<br />

+<br />

Pro sestrojení transformované LAFCH regulované určím hodnotu amplitudy<br />

[ ]<br />

A dB = 20log1−<br />

20log<br />

2 2<br />

1+<br />

4,021 Ω<br />

− 20log<br />

2 2<br />

2 2<br />

+ 20log 1+<br />

0,5 Ω + 20log 1+<br />

0,176 Ω<br />

2 2<br />

1+<br />

0,657 Ω<br />

+<br />

Při sestrojení vlastní charakteristiky dále postupuji podle postupu pro konstrukci<br />

amplitudové logaritmické frekvenční charakteristiky. Na obr. 5.8 je zobrazena černou<br />

barvou.<br />

A<br />

[dB]<br />

20<br />

*<br />

Ω 1<br />

*<br />

Ω 2<br />

2 / T<br />

Ω 3<br />

–20<br />

Ω ř<br />

G 0 (jΩ)<br />

G R (jΩ)<br />

G S (jΩ)<br />

Ω<br />

Obr. 5.8 - Transformované frekvenční charakteristiky


Strana 68<br />

Použití vybraných metod na příkladech<br />

tvaru<br />

Jestliže použiji číslicový PI regulátor, jehož transformovaný přenos je v obecném<br />

G<br />

⎛ +<br />

⎛ 1+<br />

T q ⎞<br />

( ) =<br />

i<br />

q r<br />

⎜1<br />

⎟ = r<br />

⎜<br />

T q T q ⎟ 0<br />

0<br />

i<br />

i ⎠<br />

⎝<br />

1<br />

R<br />

,<br />

⎞<br />

⎠<br />

⎝<br />

pak při volbě časové konstanty regulátoru T i rovnající se transformované časové konstantě<br />

regulované soustavy T je transformovaný přenos rozpojeného obvodu<br />

( q)<br />

*<br />

1<br />

( 1−<br />

0,5q)( 1+<br />

0,176q)<br />

4,021q( 1 0,657q)<br />

G O<br />

= 0,8<br />

.<br />

+<br />

Amplituda odpovídající jeho logaritmické frekvenční charakteristiky je<br />

[ ]<br />

A dB = 20log0,8 − 20log 4,021Ω<br />

− 20log<br />

2 2<br />

2 2<br />

+ 20log 1+<br />

0,5 Ω + 20log 1+<br />

0,176 Ω<br />

1+<br />

0,657<br />

2<br />

2<br />

Ω<br />

+<br />

Na obr. 4.4 je průběh LAFCH rozpojeného regulačního obvodu uveden červenou<br />

čarou. Pro zvolené zesílení r 0<br />

= 0, 8 a integrační časovou konstantu T i<br />

= 4,021s<br />

protíná<br />

tato charakteristika osu 0dB pod sklonem –20dB/dek při frekvenci<br />

odpovídá doba regulace podle vztahu (4.47) hodnotě<br />

−1<br />

Ω<br />

ř<br />

= 0,209s<br />

a tomu<br />

t<br />

r<br />

1,5<br />

≈ ≈ 7,18s . Ω<br />

ř<br />

Na závěr do transformovaného přenosu regulátoru s využitím vztahů (4.40b) a (4.46) určím<br />

z-přenos regulátoru<br />

G R<br />

( z)<br />

( z − 0,779)<br />

9,042z − 7,042<br />

= 0,8 ⋅<br />

= 0,899 .<br />

8,042z − 8,042 z −1


Strana 69<br />

6 ZÁVĚR<br />

Tématem této diplomové práce bylo použití frekvenčních charakteristik při analýze<br />

a syntéze regulačních obvodů.<br />

Cílem úvodní a druhé kapitoly bylo získání uceleného obecného přehledu o<br />

pojmech analýza a syntéza. Tyto pojmy zde byly rozebrány podrobně, takže má o nich<br />

čtenář dokonalý přehled a informovanost.<br />

Třetí kapitola se zabývala použitím frekvenčních charakteristik u spojitých<br />

systémů. Kapitola byla rozdělena do tří důležitých podkapitol, kde se první věnuje<br />

obecnému popisu frekvenčních charakteristik u spojitých systémů, druhá popisuje analýzu<br />

regulačních obvodů, tedy určování jejich stability a třetí rozebírá jejich syntézu, tedy návrh<br />

parametrů regulátoru. To vše samozřejmě pro spojité regulační obvody. U každé z metod<br />

byly uvedeny buď kompletní metody, nebo zkrácené včetně výpočtových vztahů a grafů.<br />

Ve čtvrté části byly, obdobně jako u třetí, zdůrazněny tři podkapitoly, které se<br />

věnují diskrétním regulačním obvodům. Obecné pojetí frekvenčních charakteristik plynule<br />

navazuje na druhou a třetí podkapitolu věnující se analýze a syntéze regulačních obvodů.<br />

Metody stability regulačních obvodů jsou podobné jako pro spojité obvody, ale byly<br />

upraveny pro obvody diskrétní. U každé z metod byly uvedeny, stejně jako u spojitých<br />

regulačních obvodů, buď kompletní metody, nebo zkrácené včetně výpočtových vztahů a<br />

grafů.<br />

V poslední, páté kapitole, byly některé uvedené metody aplikovány na vhodných<br />

příkladech, díky kterým je možné snadno pochopit danou metodu. Grafy a postupy<br />

výpočtu u jednotlivých příkladů pomohou pochopit problém a jeho řešení.<br />

Po vypracování této diplomové práce jsem dospěl k závěru, že frekvenční<br />

charakteristiky jsou dobrým nástrojem pro analýzu i syntézu regulačních obvodů.<br />

Vzhledem k tomu, že se ale jedná o grafické vyjádření, doporučuji pro vykreslování grafů<br />

a výpočet rovnic použít počítač a software typu Microsoft Office Excel.


Strana 71<br />

SEZNAM POUŽITÉ LITERATURY<br />

1. BALÁTĚ, J. Automatické řízení. 2. přepracované vydání Praha : BEN, 2004. 664<br />

stran. ISBN 80-7300-148-9.<br />

2. ČERNÝ, M. Číslicová regulace elektrických pohonů. Praha : SNTL, 1984. 208s.<br />

3. DAVIDOVÁ, O. Využití frekvenčních metod při navrhování diskrétních systémů<br />

řízení. Brno, 2004. 106s. Disertační práce na FSI VUT v Brně na Ústavu<br />

<strong>automatizace</strong> a <strong>informatiky</strong>. Vedoucí disertační práce Doc. Ing. Ivan Švarc,<br />

CSc.<br />

4. FENCLOVÁ, M. Teorie automatického řízení - Návody ke cvičení. Praha :<br />

Vydavatelství ČVUT, 1996. 152s.<br />

5. ISERMANN, R. Digital Control Systems. Berlin : Springer-Verlag, 1981.<br />

6. KUBÍK, S. Teorie automatického řízení I. Praha : SNTL, 1980. 528 s. pod. č.j. 12<br />

037/80-30.<br />

7. SUCHNA, J. Metody syntézy spojitých regulačních obvodů. Brno, 2005. 74s.<br />

Diplomová práce na FSI VUT v Brně na Ústavu <strong>automatizace</strong> a<br />

<strong>informatiky</strong>. Vedoucí diplomové práce Ing. Olga Davidová, Ph.D.<br />

8. ŠVARC, I. Automatizace: Automatické řízení. 1. vydání Brno : CERM,<br />

9. 2002. 201 s. ISBN 80-214-2087-1.<br />

10. ŠVARC, I. Teorie automatického řízení I. 1. vydání Brno : Skriptum VUT,<br />

11. 1987. 210 s.<br />

12. VÍTEČKOVÁ, M. Seřízení regulátorů metodou inverze dynamiky. Ostrava : VŠB -<br />

Technická univerzita Ostrava, 1998. 56s. ISBN 80-7078-628-0.<br />

13. ZEMAN, K. Automatická regulace v elektrických pohonech I. Plzeň : Ediční<br />

středisko VŠSE Plzeň, 1987. 220s.<br />

14. ZÍTEK, P. Automatické řízení. Praha : Vydavatelství ČVUT, 2001. 148s.<br />

15. ZÍTEK, P. Matematické a simulační modely 1 - Modely v komplexním oboru.<br />

Praha: Vydavatelství ČVUT, 2001. 111s.<br />

16. HAVEL, P. Frekvenční metody syntézy. Praha, 2004. Dostupné na WWW:<br />

http://dce.felk.cvut.cz/sri2/ss/Synteza/Havel_frekv_met2.pdf.<br />

17. Internetová učebnice na Katedře řídicí techniky ČVUT-FEL Praha. Dostupná na<br />

WWW: http://dce.felk.cvut.cz/sri2/ss/<br />

18. MINÁR, K. Vyhledávací služba Gogole začala indexovat soubory ANALÝZA.<br />

VSB [online]. Dostupné na WWW: http://www.fs.vsb.cz/books/Analyza/prvni.html<br />

19. MODRLÁK, O. Analýza diskrétních regulačních obvodů. Liberec, 2004. Dostupné<br />

na WWW: http://www.fm.vslib.cz/~krtsub/fm/modrlak/pdf/cr_analyza.pdf<br />

20. MODRLÁK, O. Teorie automatického řízení I - Syntéza regulačních obvodů.<br />

Liberec, 2004. Dostupné na WWW:<br />

http://www.fm.vslib.cz/~krtsub/fm/tr1/tar1_syn.pdf<br />

21. MODRLÁK, O. Teorie automatického řízení I - Syntéza ve frekvenční oblasti a<br />

citlivost a robustnost regulačních obvodů. Liberec, 2004. Dostupné na WWW:<br />

http://www.fm.vslib.cz/~krtsub/fm/tr1/tar1_K5.9-10.pdf<br />

22. PSTRUŽINA, K. Atlas filosofie vědy. Dostupné na WWW:<br />

http://nb.vse.cz/kfil/win/atlas1/analyza.htm

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!