Uporaba metode QFT »Quantitative Feedback Theory« pri načrtovanju robustnega vodenja Dalibor Igrec, Amor Chowdhury, Rajko Svečko FERIMaribor, Smetanova ulica 17, 2000 Maribor E-pošta: dalibor. igrec@uni-mb.si Povzetek. Namen članka je predstaviti načrtovanje regulatorja z metodo QFT, ki je verjetno edina znana tehnika načrtovanja vodenja, kjer je zajeto hkratno upoštevanje faze in negotovosti parametrov objekta. Prednost metode je možnost doseganja robustne stabilnosti in robustnega učinka z minimalnim učinkom povratne vezave [7]. QFT je grafično-analitični postopek načrtovanja vodenja, ki zahteva precej predpriprav pri oblikovanju vzorcev objekta in empiričnih izkušenj, obenem pa daje načrtovalcu precej manevrskega prostora in neposrednega vpogleda v spremembe regulatorja pri načrtovanju. Temeljna hiba metode je določanje mej objekta v Nicholsovem grafu, saj izračun meja metode QFT eksponentno narašča z natančnostjo vzorca objekta. Za nazornejšo predstavitev metode smo predstavili poskus na realnem objektu, kjer smo izvedli vodenje sistema z regulatorjem načrtovanim z metodo QFT ter preverili lastnosti sistema in primerjali dobljene rezultate z rezultati, ko je bil sistem voden z regulatorjem, načrtovanim po metodi H(. Ključne besede: metoda QFT, Nicholsov graf, negotovost parametrov, vzorec objekta, robustna stabilnost Robust control design with QFT Abstract. The quantitative feedback theory (QFT) design is a graphic technique for designing feedback controllers which is probably the only known technique that simultaneously considers large parametric uncertainty and phase information. The ability to satisfy robust stability and different performance constraints with the minimum possible cost of feedback is the biggest advantage of the method. The downside is that the method, though systematic and powerful in hands of an experienced control engineer, has only recently lent itself to a formal mathematical form as is the case with the more recent paradigms such as H control and /-synthesis. A major advantage of QFT is that the design is performed in the frequency domain. This enables a good insight into the plant operation and difficulties that may arise during the controller design. Uncertainties can be caused either by changing the plant characteristics or ambient conditions or by unknown external disturbances. QFT starts by defining the plant and then specifying its uncertainties. The defined uncertainties are then used to determine the differential gain and phase from the nominal ones, over the range of frequencies through which the plant operates. At each distinct frequency, differential gains and phases are used to generate the Plant template. The given example illustrates the steps taken in the QFT controller design. The aim of the present work is to present the usage of the QFT method for the controller design. To allow for a more illustrative presentation we made an experiment with a real object with the controller designed according to the QFT method. In the paper we show the complete procedure of the QFT design from the model analysis to the controller design. At the end we also compared the system performances of the QFT and H( controller design. Keywords: QFT design, Nichols chart, parameter uncertainty, plant template, robust stability 1 Uvod Načrtovanje vodenja sistemov je pogosto povezano z negotovostjo sistemov in težavami pri uporabi ustreznih regulatorjev, vzdrževanju stabilnosti sistema in preprečevanju neželenih učinkov, kot npr. eksternih motenj in perturbacij sistema. Za odpravljanje učinka negotovosti se uporabljajo adaptivni regulatorji, kjer se parametri objekta (proge) sproti identificirajo, kar ustrezno vpliva na dinamično sintezo regulatorja in robustni regulatorji, ki upoštevajo najslabše mogoče razmere celotne družine modelov objekta s stališča negotovosti ob nespremenjenem regulatorju. Delovanja sistemov v realnosti nikoli ne poznamo v celoti, saj se njihovo delovanje spreminja s časom zaradi staranja komponent, spreminjanja parametrov ali s spremembo delovnih razmer (obremenitve, motnje). Cilj načrtovanja je izdelati takšen regulator, ki bo zagotavljal stabilno in performančno zadovoljivo delovanje vodenega sistema v realnem okolju. Metoda QFT načrtovanja robustnih regulatorjev [2] se je razvila v zadnjih 25 letih sistematičnega ukvarjanja z negotovostjo modelov. Začetki reševanja problematike segajo v leto 1980, ko sta Gera in Horowitz objavila članek o uporabi Bodejevega ojačitveno-faznega integrala za ugotavljanje karakteristike nominalne zanke z iteracijskim postopkom [12]. Postopek ni bil vedno konvergenčen, prav tako je bila za rešitev potrebna aproksimacija. Postopek je bil avtomatiziran z uporabo orodja QFT Toolbox [1], ki je poenostavil postopek iteriranja in uporabo aproksimacij višjega reda. Thompson in Nwokah [7] sta za izračunavanje uporabila nelinearne tehnike programiranja, kjer sta meje funkcij QFT določila s pretvorbo funkcij QFT v funkcije H(. Bryant in Halikias [3] sta uporabila linearne tehnike programiranja, vendar so bili njuni rezultati močno poenostavljeni in zaradi neupoštevanja polov in ničel funkcij nesposobni zagotavljati stabilnost sistemov. Zhao in Jayasuriya [9] sta uporabila Youlovo parametrizacijo za transformiranje funkcije QFT v enodimenzionalen problem, vendar se s tem lahko avtomatično spreminja le en parameter regulatorja. V primerjavi z drugimi optimizacij skimi metodami načrtovanja robustnega vodenja ima metoda QFT nekaj prednosti, kot na primer možnost kvantitativnega ugotavljanja "izgube povratne zanke", možnost upoštevanja faznih zamikov v procesu načrtovanja in možnost izbire med zahtevnostjo načrtovanja procesa in kompleksnostjo regulatorja. Omenjena prednost je še posebno pomembna, saj omogoča izdelavo preprostih regulatorjev, ki jih je mogoče tudi preprosto implementirati. 2 Princip načrtovanja Cilj robustnega vodenja je zadovoljivo performančno obnašanje vodenega sistema kljub nepopolnemu poznavanju objekta. Zadovoljivo performančno obnašanje lahko opišemo kot ojačenje |H(j®)| prenosne funkcije zaprte zanke, ki leži v dovoljenem območju Bodejevega diagrama na sliki 1: 0 <|a( jm)\ <¡H(jm)\ <\p{ jm)\ (1) Amplituda želene prenosne funkcije zaprte zanke je do določene frekvence enaka 0dB, nato pa hitro pada. Ponavadi so meje a (m) in p(a) določene tako, da je sprememba \p{m) -a(m)| dovolj majhna, da ne povzroča nestabilnosti in zadovolji performančne kriterije, v frekvenčnem območju m< mh. mh pomeni frekvenco pri absolutni vrednosti razlike mej \p{m)-a(m)| .Pri izbiri mej a(m) in p(m) je nerealno določiti ozko spremembo nad frekvenco mh, saj je naraščanje občutljivosti sistema v visokofrekvenčnem območju neizogibno. Predvidevamo, da so vrednosti |P(a>)ja(m)| pri mh zelo velike, medtem ko so dejanske vrednosti |a(m)| in \p(m)\ majhne. 10 5 o „ -5 I S -10 -20 -25 -30 Slika 1: Meje tolerance za ojačenje zaprte zanke Figure 1. Tolerance band for the closed-loop gain Za doseganje robustnosti je treba zadostiti vsem pogojem, prikazanim na sliki 1. Pri objektih z ne-minimalno fazo moramo za fazo prenosne funkcije zaprte zanke dodatno definirati tolerančne meje [6]. Elektromehanski objekti so lahko predstavljeni kot modeli, ki vsebujejo strukturirano in nestrukturirano odstopanje [11]. Nestrukturirana odstopanja so zajeta z metodo QFT prek omejitve amplitude zaprte zanke, ki zagotavlja stabilnost za vse objekte iz družine objektov z nestrukturiranim odstopanjem. Cilj metode QFT je, da dosežemo robustno stabilnost in hkratni robustni učinek za družino modelov objekta. Med vsemi regulatorji, ki zadostujejo pogojem, imenujemo tistega z minimalnim visokofrekvenčnim ojačenjem optimalni regulator [7]. Postopek načrtovanja regulatorja po metodi QFT je sestavljen iz več korakov. V prvem koraku izberemo prenosno funkcijo nominalnega objekta P0(s), ki ga dobimo iz lineariziranega modela objekta. Koeficiente modela izberemo iz njihovih intervalov negotovosti. Edina omejitev pri tej izbiri je, da mora nominalna zanka L0(s) zadostiti pogoju ničelne izključitve [8]. V drugem koraku izberemo množico kritičnih frekvenc tako, da bodo meje pri izbranih frekvencah zadostovale pogojem za vse frekvence v opazovanem območju [1]. Frekvence so empirično izbrane po priporočilu v [14]. Z večanjem števila izbranih frekvenc se povečuje zapletenost načrtovanja. S tretjim korakom določimo vzorce, ki prikazujejo učinek negotovosti parametrov na ojačenje in fazo nominalne zanke pri kritičnih frekvencah. Določimo meje stabilnosti, ki so vodilo pri načrtovanju regulatorja C(s) za kritične frekvence. Mejo stabilnosti vseh kritičnih frekvenc sestavljajo nominalne točke, vzorec se dotika kroga območja stabilnosti v Nicholsovem grafu (slika 2). Iz tega sledi, da je maksimalna dovoljena amplituda odvisna od omejitvenega kriterija p(a>), ki je supremum zgornje meje sup|. Tako dobimo mejo stabilnosti za nizke frekvence. Postavitev meje stabilnosti za visoke frekvence si zamislimo tako, da postavimo pisalo na nominalno točko vzorca za določeno kritično frekvenco. Nato celoten vzorec premikamo okoli kroga območja stabilnosti, medtem ko še vedno držimo pisalo v nominalni točki vzorca. Slika 2 prikazuje, kako v korakih pri določeni kritični frekvenci premikamo vzorec v smeri urnega kazalca okoli kroga območja stabilnosti. Obris na sliki 2 je meja stabilnosti za pripadajočo kritično frekvenco. Dobljen obris je na zgornji strani fiksno določen z območjem stabilnosti (polna črta), na spodnji pa odvisen od oblike vzorca (črtkana črta). Potem določimo mejno črto pri kritični frekvenci. Točko mejne črte dobimo s postavitvijo nominalne točke vzorca na določeno fazo odprte zanke in s premikanjem vzorca navzgor in navzdol tako dolgo, dokler ne najdemo najnižjega položaja, kjer je sprememba ojačenja zaprte zanke enaka dovoljeni spremembi, določeni z mejo tolerance. Amplitudo zaprte zanke razberemo iz zaprtozančne mreže v Nicholsovem grafu (M-linije). Načrtovanje regulatorja z oblikovanjem odprte zanke izvajamo tako, da dosežemo: - robustno stabilnost: to pomeni, da je pri vsaki kritični frekvenci pripadajoči odziv odprte zanke desno od pripadajoče meje stabilnosti ali na njej, - robustni učinek: to pomeni, da je pri vsaki kritični frekvenci pripadajoči odziv odprte zanke nad pripadajočo mejno črto ali na njej. Spreminjanje oblike pomeni spreminjanje ojačenja ter dodajanje ničel ali polov v odprto zanko. Najpomembnejša zahteva je, da mora L (s) = C (s)P (s) zadovoljiti Nyquistov stabilnostni kriterij [14]. Pri nestabilnih modelih in modelih z neminimalno fazo obstajajo določene omejitve pri izvajanju metode zaradi povezave med amplitudo in fazo, ki sta podana s Hilbertovim transformom in Bodejevim integralom [14]. Omejitev zmanjšuje svobodo pri izbiri strukture regulatorja in v določenih primerih vsiljuje nasprotujoče si zahteve, posledica pa je, da regulatorja sploh ni mogoče načrtati. Osnovni elementi, ki jih dodajamo v odprto zanko objekta L(s) med načrtovanjem regulatorja, so: - ojačenje: k, p s + p Priporočila za uporabo predstavljenih elementov pri načrtovanje po metodi QFT so podana v [14]. pol ali ničla: element lead ali lag: s + b 2 . O e „„ I „2 jri.L t. t , i s + 2£,ms + m filter Notch: —-21- s + 2^2ms + m (2) (3) (4) meja stabilnosti za določeno kritično frekvenco -350 -300 -250 -200 -150 -100 -50 0 Odprto-zančna faza arg(L(j (s) = 0, W3(s) = VÔÎ5s + 250 (18) regulatorja H(: _ 1.105 • 107 s3 + 3.823 • 108 s2 + 2.923 • 109 s + 3.934 • 108 (s) _ s3 + 2.360 • 103 s2 +1.251 • 106s (19) Slika 11 predstavlja stopnične odzive vodenega sistema z regulatorjem QFT pri treh različnih razdaljah med kroglo in elektromagnetom (5, 12 in 20 mm). a in / sta stopnična odziva performančnih mej a(o>) in P(o) omenjenih v drugem poglavju. Občutljivost in komplementarno občutljivost sistema smo preverili v treh delovnih točkah (slika 10) in ugotovili, da regulator QFT zagotavlja sistemu zahtevano občutljivost in komplementarno občutljivost v celotnem delovnem območju. Slika 11: Stopnični odziv zaprte zanke sistema Figure 11. Closed-loop system step response Slika 12 prikazuje odsekoma stopnični odziv sistemov, vodenih z regulatorjem QFT in z reduciranim suboptimalnim regulatorjem H( na celotnem delovnem območju. Razvidno je, da performančno sistem, voden z regulatorjem QFT, deluje bolje kot sistem, voden z reduciranim suboptimalnim regulatorjem H(. Poleg direktne primerjave odzivov je pri primerjavi metod treba upoštevati tudi, da je prenosna funkcija suboptimalnega regulatorja H( pred redukcijo sedmega reda. Slednje je vsekakor pomanjkljivost pri izvedbi suboptimalnega regulatorja H(, ki pa se jo da deloma kompenzirati z redukcijo. Slika 10: Občutljivosti sistema za tri različne odmike krogle Figure 10. System sensitivity at three varied deviations Sistem z regulatorjem QFT smo preizkusili tudi performančno z odzivom na stopnično vzbujanje in ugotovili, da regulator zagotavlja zahtevano obnašanje sistema v celotnem delovnem območju. Delovanje sistema, vodenega z regulatorjem QFT, smo primerjali tudi z delovanjem sistema, vodenega z regulatorjem H(. Reducirani suboptimalni regulator H( smo zasnovali po postopku [10]. Utežnostne funkcije za določitev performančnih zahtev in modeliranje dinamike nestrukturiranih odstopanj pa smo določili po priporočilih [11]. Utežnostne funkcije metode H(: ( 250 V Prenosna funkcija reduciranega suboptimalnega Slika 12: Primerjava odzivov v delovnem območju Figure 12. Response comparison in the working range 6 Sklep Metoda QFT je grafična metoda načrtovanja regulatorjev za SISO in deloma MIMO negotove objekte. Pri načrtovanju regulatorjev za MIMO sisteme je metoda QFT uporabna, če se problem lahko razgradi v preproste nepovezane probleme SISO. Osnovna ideja metode je razdeliti proces načrtovanja na več stopenj, od katerih vsaka pomeni poenostavljen problem povratne zanke SISO ali MISO. Ker proces zahteva veliko računanja in grafične podpore, so računalniški programi še posebno primerni za načrtovanje [1]. Pri enostavnih objektih SISO metoda QFT dosega podobne rezultate kot druge metode načrtovanja robustnega vodenja. Prednost metode QFT je, da ima načrtovalec neposreden vpogled v potek načrtovanja, saj sam postavi strukturo regulatorja. Prav tako ima možnost upoštevanja faznih zamikov v procesu načrtovanja, rezultat pa je znižana struktura regulatorja, kar npr. pri metodi LQG, H( ali H2 nima [10]. Kljub določenim prednostim metode QFT še vedno ostajajo pomanjkljivosti. Načrtovanje po metodi QFT zahteva nabor značilnih dinamik objekta za celotno območje variacij negotovosti. Žal pa preslikava negotovosti modelov v frekvenčni prostor ni enolično določena. Zato je včasih težko določiti obliko vzorcev objekta, ker ni transparentno, kako fizikalni parametri vplivajo na obliko. Rezultat tega načina ni optimalen zaradi svoje odvisnosti od izkušenosti načrtovalca. Metoda QFT se izkaže za uporabno pri načrtovanju vodenja kompleksnih sistemov SISO. Pri metodi QFT načrtovanje ni avtomatizirano in je načrtovalec potreben na vsaki stopnji načrtovanja. Pozitivna lastnost tega je, da ima načrtovalec nadzor nad celotnim postopkom načrtovanja in se lahko odloči za zanemarjanje določenih specifikacij oziroma omejitev. Načrtovanje z metodo QFT je postopek poskusov in napak, ki je močno odvisen od izkušenj načrtovalca. Rezultat načrtovanja z metodo QFT so regulatorji nižjega reda, ki so uporabni zaradi obvladovanja strukture regulatorja. Vsi poli in ničle so rezultat metode in ne samega analitičnega postopka kot pri metodi LQG, H« ali H2. 7 Literatura [1] C. Borghesani, Y. Chait, O. Yaniv, The QFT Frequency Domain Control Design Toolbox, 2001. [2] Isac M. Horowitz, Quantitative Feedback Design Theory (QFT), QFT Publications, Boulder, Colorado, 1993. [3] R. Nordgren, O.D.I. Nwokah and M.A. Franchek, A New Perspective on the Formulation of QFT, Proceedings of the American Control Conference, 2, 1716-1720, San Francisco, California, 1993. [4] K. Zhou, John C. Doyle, K. Glover, J. Doyle Robust and Optimal Control, Prentice Hall, 1st edition, 1995. [5] C. Houpis, John D Azzo, Linear Control System Analysis and Design, The McGraw-Hill Companies, 1995. [6] I. M. Horowitz, M. Sidi, Optimum synthesis of non-minimum phase feedback system with plant uncertainty, Int. J. Control, 1978. [7] Thompson, D.F., Gain-Bandwidth Optimal Design for the New Formulation Quantitative Feedback Theory, J. of Dynamic Systems, Measurement, and Control, V. 120, No. 3, 1998. [8] J. Ackermann, A. Bartlett, D. Kaesbauer, W. Sienel, and R. Steinhauser, Robust Control: Systems with Uncertain Physical Parameters, London: Springer-Verlag, 1993. [9] Zhao, Y., and Jayasuriya, S., Robust Stabilization of Uncertain Systems with Parametric Uncertainties, Procs. 12th IFAC Conf., Sydney, Australia, Vol. 6, 1993. [10] A. Chowdhury, Robustna sinteza regulacijskih sistemov z upoštevanjem performancnih kriterijev : doktorska disertacija. Maribor, 2001 [11] K. Zhou, John C. Doyle, Essentials of Robust Control, Prentice Hall, 1998 [12] Wenhua Chen, Donald J. Ballanc, Plant Template Generation in Quantitative Feedback Theory, University of Glasgow, 1998 [13] Wenhua Chen, Donald J. Ballanc, On Choice of the Nominal Plant in Quantitative Feedback Theory, University of Glasgow, 1997 [14] O. Yaniv, Quantitative Feedback Design of Linear and Nonlinear Control Systems, Kluwer Academic Publishers, 1999 Dalibor Igrec graduated in 2000 from the Faculty of Electrical Engineering and Computer Sciences of the University of Maribor, Slovenia. Since 2001, he has been employed as a researcher with Ultra d.o.o., R&D center Maribor. His work is focused on mobile communications networks and solutions. Amor Chowdhury received his B.Sc., M.Sc. and Ph.D. degrees in electrical engineering from the Faculty of Electrical Engineering and Computer Sciences of the University of Maribor, Slovenia, in 1994, 1997 and 2001, respectively. Since 1994, he has been employed with the Institute of Automation, Maribor as a researcher. His main research interests include control systems design and H2 and Hw theory. Rajko Svečko received his B.Sc., M.Sc. and Ph.D. degrees in electrical engineering from the Faculty of Electrical Engineering and Computer Sciences of the University of Maribor, Slovenia, in 1981, 1984 and 1989, respectively. Since 1981, he has been employed with the Institute of Automation, Maribor. He is the Head of the Laboratory for Control Systems. From 1995 he has been Associate Professor. His main research interests include system theory, intelligent control systems, classical and modern control design.