• 検索結果がありません。

Millerの最適オブザーバ併合系の実験による考察 : 二重倒立振子の安定化制御において: University of the Ryukyus Repository

N/A
N/A
Protected

Academic year: 2021

シェア "Millerの最適オブザーバ併合系の実験による考察 : 二重倒立振子の安定化制御において: University of the Ryukyus Repository"

Copied!
11
0
0

読み込み中.... (全文を見る)

全文

(1)

Title

Millerの最適オブザーバ併合系の実験による考察 : 二重倒

立振子の安定化制御において

Author(s)

石田, 力; 長堂, 勤

Citation

琉球大学工学部紀要(34): 209-218

Issue Date

1987-10

URL

http://hdl.handle.net/20.500.12000/1482

Rights

(2)

琉球大学工学部紀婆館34秒.1987年

209

Millerの最適オブザーバ併合系の実験による考察

一二重倒立振子・の安定化ルリ御において-

石田力.長堂勘**

considerationfortheRegulatorwith

Miller'sObserverthroughtheExperiment

-ThroughStabilizingControlforaDoublelnvertedPendulum-TsutomulSHIDAcandTsutomuNAGADO…

AbStract

Whenwewerestabilizingourdoubleinvertedpendulumbyapersonalcom・

puter,weobtainedthedatathatthe妃gu]atorwithMiller,sobserverisrobust

compairedtotheomewiththeobserverhavingthesameobserverpolesas

Mi1ler,sbutthedifferentobservergainsfromMiI1er's・Thenwereportit,

Keywords:RegulatoT,ObserverDStabiIizing,Doubleinvertedpendulum,

Robustness 1.まえがき 2.オブザーバ 状態フィードバック制御をオブザーバを併用して 実現する際に,オブザーバの極およびオブザーバゲ インをどのように決定するかということが問題とな る。Millerはシステムの初期{ihの平均値と分散が既 知のとき,オブザーバを用いることによる評価関数 の劣化魁△Jを最小にするという意味において破適 なオブザーバゲインの決定法を提案した:〕 本摘では当研究室で作製した二愈倒立振子の安 定化実験において,蝋適レギュレータと併合した MillerのIF上適オブザーバが同じオブザーバの極を持 つ他のオブザーバ併合系と比較してロバストである という実験データが得られたのでそれを報告する。 2.1最'】吹元オブザーバ 状態フィードバック制御を行うには, 状態フィードバック制御を行うには,システムの すべての状態を入手する必要があるが,入手できな い状態が存在する場合,(2.1),(2.2)式のよう なダイナミックシステムにより状態を推定する。 へ 。)(IノーAmrJノ+KWIノ+庇、ノ (2.1) 昂イノノーDCWYj+〃9W (2.2) ただし,J仰1Wはそれぞれプラントの出力変数, 入力変数であり,畠のは状態変数の推定値である。 受付:1987年5月9日 ・琉球大学工学部電子・禰報工学科 Dept・ofEIectronicsandlnformationEngmeering,Fac・ofEng. ・・琉球大学院研究科 GraduateStudent,ElectricalandlnformationEngineeringb

(3)

MiIlerの峻迩オブサーバイル介系グ)リミ賊による必察:イjIIl.Lと堂 210 ただし,ヱィノノは状態変数,q、(は'、みイj列,ノj(よ次 のりか,チノノイM式のjIf遮解である” ilII定不1i丁能な状態変数のみを推定する岐小次)Gオブ ザーバでは,ゴピナスの正f仏形式を11Iいて升係'数イ↑ 列を決定する。特にオブザーバの拳jlUlに人きな膨料i を与えるAは パ1.1jトルルI(J--lW/111ノィ・I1J-U (2.7) P /I=/l22--MIHR オブザーバを川いた(ル介ンステムでは.状態蛮数 エイ"の代わりにその批定IlIl【であるjMを使うので人ソj 1Uイ"は, Mノーィ(’/jrノjf「リノ (2.8) となり,lItj1m{IiI1伽人ノノで【よむいのでlI1然,ilillⅡillL1激./ のllIIより7化-1る”MiⅡc『‘〕腿Ⅱハした〃UiIjLロブ・ル ト〃》IiUlUlIIIlの銃,ilrl1な'rl1iMMベルか.’ているときに, こび)0岬li1ilM散の坊化111ムノセ11t'|、にWIDようにオブ リ`-バゲインLを決定1.る(」のでfDる.、こぬノノi」:Iよ 〃)慨を柵定する必狸がなく.]没I1llfのIilⅡがかiAF I)擁iLIビされる緬以下にそび)繊パlTLIIlIiを,l《す31 (1)状態内初101仙の、14均伽〔と此分1M(イ「ダ11が次のよう に!』えられるとする。 (2.3) のように疋まる“(23)式においてjが安定行ダリ となるように,没,汁パラメータであるオブザーバゲイ ンLを決定しなければならない。しかし多111力系の 」ルイfには/iの櫛をlfiiIとしてもノ`には111111災が'’1) llIlM-には定f',ない〃さらにまた〃)樋をどこに111 》とすべきかというI1IIlu6ダ“。そこでMiIlerLMIt jlmレギ百ルークのIlliIIIllIU数にil;|イし,オブザーバゲ インLcノ)決巡にイ「効とノノijjを腿樅している“ 2.2Millerの蛾適オブザーバⅡ (2.4) Jfr〃=ハェィルトノイ『((〃 E(ェ⑩'〕=畑。 なるシステムに対し,(25)式の二次形式(`liliHillU 数をもつ股適レギュレータ|H1題の肢適入力は(2.6) 式のようになる。 (2.9) E〔化(01-'1m)(工(O)-地,)T〕=塾

[:M糞1

J=」1M"…+腱冊JR“('′M1

(2.5) (2.10) uイノノーーノjrr〃 =-lt-llfr「〕エィノノ (2.6) ただし,sはj911桃([を表わす演算子である。 ,り卍・F1

蔀Iii…

E最Qi:ii竺幽

lmuqn POtGntio.’ e 図1二1ft倒立振子系の構成図

(4)

li酢k大学l:学麟紀撰堀34.ザ・’1987年

211 (H)(2.11)式の仮想システムに対し、(2.12) 式②皇次l形式iⅡ;IiniilLI数をjlt小化するjiとj凶iIiI1伽Ⅱリ (2.」3)式を派ぬる!, 3.シミュレーション 本節では当1ili究鑓で作製した二JK倒立搬子系を MillcI・の妓適オブザーバ併合系により安定化すると きのシミュレーション結果を示す。二m倒立限子シ ステムは図1のような榊成になっており,その状態 空lillモデルは以下のようになる。

りi"=('ll2-AHZI1Z11)〃i"十J1,2W'

・I.、. (2.11)

ノミ1W.'ヅ'カバど蕊-ど即画Bi遵聰)w1jやざhldl

(2.12) ・r lrl'("=-.-/`2リイノノ (2.13) f(n=ハエru)+Bmrfj W"=CMノ ただし状態変数,人力変数および杵係数行列は エイ'ノー〔識、LlI,仇-01,A,‘,、A-j,〕T 町〃=(P (3)オブザーバゲインLと(`リィィノグ〕初lU11i(〔を次のよう にik達する。 1.=Z2lZ7ll+L2 ⑪((1)=〔-l、1m"`,

帷ili

l-肌

jM

1 0 0 -18.2 64.5 -77.4 0 1 0 0.025 -0.32 0.269

0 1 0 001 000 である。 レギュレータについては,(2.6)式の舷適レギ ニルータを榊成する。爪み行列Q,ノヒを以下のように 股通してシミュレーションを行う。 Q=diaR〔0.2220000〕 R=l このとき職適レギュレータの極は -24.2±j18.3 -2.79±j219 -L24±j0.41 となり,レギュレータゲインは F=〔0.4454.21811.350.4761.761 1.548〕 となる。 次にMillerの般適オブザーバを股計する。繩,2. を次のよう膣与えてシミュレーションを行う。 湖⑨=0

ユーdiag〔0.10.10.1000〕

このとき』とLは

(5)

Miuerの賦適オブザーバ併合系の実験による考察:r1ul・催堂 212 1- 一一一一 〈八L

,liliil

二|劃

I-L;蕊

1劃

-19.14 1.47 -1.75 -1.44 12.58 -7.63 -14.83 14.54 0.269 0.934 1.465 -1.749 (3.2) 1.465 12.58 -7.627 以上の二つのオブザーバについてそれぞれWLIルー プ系をlilIi成し,シミュレーションを行った。レギュ レータについてはI1liシミュレーションにおいて全く liil-のものを仙川Lている。状態変数およびオブ ザーバの例JIM、を次のようにしたときのこ疏倒jxlli f燐定化のシミュレーション紬来をlxi2~Ⅸ17に,パ リ-. 災'01-〔I).l(】」()10()()J「 (3.1) となる。このときのオブザーバの蝿は -烈入7土ノ14.9 -6.34 MⅡ)=〔(〕()(〕〕.『 となる。比絞のためにMi化「。')舷jliオブザーバと |可ヒォプザーパの櫛を抄つが’別のオブザーバゲイ ンを有するオブザーバを櫛砿する。その一例として 次のような(,Lを採用する。 は0M上越をもつ別のオブザーバ併合系に比べ,娠これらの図より.Millerの峨適オブザーバ併合系 幡が小さく収束も速いことがわかる。これはMiller の峨遮オブザーバのオブザーバゲインLのノルムが 比鮫的小さく、そのために初期再現誤差が小さくな り,応答が良くなったものと魁われる。(付録参照)

に;;

14.9

41

-20.7 0 押し バー 9 - & [Eu] ・ ・⑭け い晒垣e鶴鞄 ] 9. ̄ 図2台Ⅱiの位置の振子舞い

(6)

琉球大学工学部紀要第34号,1987年

213 ○毎 。 [山⑩ロ]|の圏梗e什曠1灘

1.

] ○ぬI 図3第一振子(下の繊子)の角度の娘子郷い C頃 ・ 肩・ [則①□]竜‐ご鉦や科e伴暖I漣却伴喋Ⅱ灘 ]

図4第一披子(下の振子)と第二搬子(上の磯子)のなす角度の扱子舞い

。●西 C [い、目]畑囲倒e鶴釦

企岑--房二=

□・面I 台叩の迎皮の振子難い 図5

(7)

MiⅡerの岐迩オブサー′{イノド令系〃)Pノミ験による稀祭:イI.Ⅱ1.1違蝋 214 涯鈴、つ一 つ 〔⑪、已甸』]|囮髄掴極S膀堕1慨 --- ̄ 白T[$】 a 3 01 M 11 揮碗へC-b MIII『・I・のオーブザーバ lj1化伸をもつMIIのオブザーバ 図6軸・ル「‘)/「IjQkl哩`/)1hK''1W:い 侭急へつ一 [⑫、口輸』]己 C R碗へ。〒 凹色遡燭極e虹や科e卜瞳1麗却峠鴨川雛 ST[s] 3 4 Millerのオブザーバ 同じ極老もつ別のオブザーバ 図7輔一縦「・と箙二縦「・のなす『(lの角速度‘)賑7.1Niい 4.実験 4.2制御系股貯 12.5)JICのlld 12.5)JICの11t過しキー‐ルータに」jける」[ぬii ダ11(J,(の選び〃については現、;のところ体系的なノノ 法はなく,小火験ではノ(=lと1M定し.Qをいろいろ 変えてシュミレーション結采を見ながら決定した。 そのときのQは Q=diag〔0.181.818000〕 となり,レギュレータゲインFは F=〔1.1815.92912.891.0642.227 1974〕 災際c'xニノMIl1>:11,t1..糸にMi化「の雌jUhオブザーバ IjI介系を適川し、そのイ「効lfliを検jけする《, 4.1実験叢歴 実験装置の榊成図を図1に示す。制御装謎として パーソナルコンピュータを使用しており.状態フ ィードバックを行うレギュレータと状態推定を行う オブザーバがプログラムされている。数値演算プロ セッサを使用して潤迷化をはかっている。

(8)

flit球大半:T:IF族IHI紀要jOMり・’1987年 となった。そのときレギエレータの極は以下のようで,このときの/i及びLは IこlHLiiHされた。

‘{篝二,

-24.2±ノ16.93 -2.84±ノ2.14 -1.73±j1.28

巖ツツ蝋為鰯鯨蝋:ルーL菫1-1

〃法をとった。実験に用いた、..Z・は 215

←に篝二嚢靭

←L霊_蕊苣讓|

のようになった。このときのオブザーバの極は

?"⑥=O Zm=diag〔O」O」0.1555〕 -22.3±ノ16.3 -7.09 R E:lIU U  ̄:20 l(b no -lD -2n -3⑥ ワ】(皿八壱 ロー [ぬ①ロ〕 onB2 』I ⑱》 一 2⑪40}4宮 [四じつ]』へUl向へ四 -fF -lZ 図8Mille「のオブザーバを用いた時の応答

(9)

Millerのl戯趣オブザーバ併合系の災験にょあみ箙:ィiⅡ1.1と堂 216 となった。比較のために上述のMillerの岐適オブ ザーバとliilヒオプザーバの極を特つか別のオブザー バゲインを有するオブザーバを柵旗する。そのオブ

ザーバの』とLは以下のようである。

4.3実験結果 Millerの收適オブザーバ併イチ系をli脱た時の応騨 例をlXl8に.Iiiルオプサーバの極をもつ別のオブ ザーバイル介系をⅡ川左IIlrのL心縛を脚!)に,I《す。二れ らは外乱として二jIMjlM:雌/・@')l:のIAtrを一定の強 さでulIいた時の応禅である。AIj:線はシミュレーショ ン結果をjP戒す。11,図を比核すると.Millerの岐迩オ ブザーバ併合系を用いた方が制御性能が良いことが わかる.実験は数回行なわれたが,系の振灘いはそ オLぞオLtXI8,偶19とほばliil様であった、爽嶮胤蛾を

に蕊

lZYili

」鬘l

-1.89 -14.0 11.7 バー 1.90 13.8 -11.7 L= tXllOに71《しておく,, 1J ■q。 [E⑥] -】 -2 -忽 I [四ロマ] 。 ロミーハーハハムヒハワ4 0 -1 l [凶①ロ] へU1円へ四 qヂ -1 図9Millerのオブザーバと同じオブザーバの極を持つが別のオ ブザーバゲインで構成されたオブザーバを用いた時の応答

州⑪⑪、、ひq 2840 I 8 ■ 句』 ■ ⑩} ■ 8 ■ 1 『 ⑪ ■ 4 ■ 、 ■ 2

(10)

琉球大学工学部紀要節34号,1987年 217 プザーバ併合系に比べてロバストな制御系である といえる。今後はこのMillerのオブザーバ併合系 のロバスト性について理論的検証を行っていきた い、 8.あとがき 本稿においては、当研究室で作製した二通倒立振 子の安定化、11御実験において,Millerの股適オブ ザーバ(沖合系がMillerの職適オブザーバと同一のオ ブザーバの極を持つ別のオブザーバ併合系に比べて 線形系から非線形系へのシステム変動に対し,ロバ ストであるという実験データが得られたのでこれ を報告した。今後は,このMillerの熾適オブザー バ併合系のロバスト性の理絵的検証を行っていき たい。

図10二、倒立振子の安定化制御実験風景 5.考察 二m倒立搬子の制御装置の設計(レギュレータ, オブザーバの係数の決定)過程を考える。二玖倒立 振子系は非線形であるが制御装置の股叶は非線形系 を線形化した線形系で行う。この線形系で殻31.した 制御装髄を現実の非線形系に適用する。このことは 股計時のシステムに一種のシステム変動が起こった と考えることができる。線形系で設計したものをも との非線形系に適用した時制御性能が良いというこ とは,上述の譲論を考慮すれば,Millerのオブザー バ併合系はこの実験においてはロバストな制御系で あるといえる。ちなみにシステム変動が起こらない 場合Iすなわち,非線形系を線形化した線形系で制 御装湿を汲:|してそれを線形系の数式モデルに適用 してシミュレーションした結采(図8,図9の点線の 波形)を見ると,Millerのオブザーバ併合系を用い た方が,同じオブザーバの極を持つ他のオブザーバ 併合系を用いた系よりもわずかに優れてはいるもの (1110》 のほとんど変わらなし、ことがわかる。したがって, システム変動がない場合はMillerのオブザーバ併合 系も同じオブザーバの極を持つ他のオブザーバ併合 系もたいした差はないが,線形系から非jMl形系にシ ステム変動があったときにMillerのオブザーバ併合 系は同じオブザーバの極を持つ他のオブザーバ併合 糸よりMlll御性能が優れているわけであるから. Millerのオブザーバ併合系は議鎗している他のオ 鮒 辞 本研究を進めるにあたり賀璽な示唆を与えて下さ った本大学工学部堀子・傭報工学科の山下勝巳氏に 心より感DMの窓を表わします。また二重倒立撮子の 安定化制御システムを作製してくれた牧野良二君 (現在,九州窟士通システムエンジニアリング),な らびに,シュミレーションを行ってくれた照屋秦一 君(現在,沖繩箇士通システムエンジニアリング) に懸鮒いたします。 参考文献 (1)R、AMiller:SpecificoptimaIcontTolof theIinearregulaterusimgaminimalorder obsewer,1,t.』・ComtroI,Vol、18,11,1,pp、 139-15901973. (2)古田I美多,川脇,原:メカニカルシステム制 御,オーム社,1984. (3)牧野良二:琉球大学工学部電気エ学科卒業研究 稔文昭和60年度 (4)長堂助:XipMt大学工学部冠子傭報エ学科卒業 研究鎗文昭和60年度 (注1)実際の実験においては,Millerのオブザーバと同じオブザーバの極を持つ他のオブザーバにおいて Millerのオブザーバゲインと極蝋に異なるオブザーバゲインを用いると二重倒立振子が立たなくな るのでMillerのオブザーバゲインと若干異なるオブザーバゲインを用いた。そのためにシュミレー ションにおいてはたいして麓が出なかったものと思われる。

(11)

218 Millerの11t鰹オブザーバIjr行系のリビ験による捗察:ィiIII.【之`;Ir

付録肘をⅡi位行タリとするとHW'14f次グ)よJMIこ炎わせる.

ゴピナスの、':'w…'岬…パビュ'パI』鳥"'繧し(`!〕

次のような形になる③ jlDIli、'"ノ勿初IUlMi⑰JII''は群とする〃)で.オブ.Ⅲ`-ハ

…'[ルルトs-'[小,〃,M,錐榊M`…,"`,胸,,,,,M1瞳,…`…

かってj、のノルムが大きければ,初lUliV現談艦か火

ただし5は適当な行ダリWにより次のように柵I皮きれきくなる..MiⅡerのllkjISiオブザ.-(の,役I;1.鱗によ

ぁ(《l」)火の1.は他刎没iil・洪による(3.2〕式の

r・に比べてノ〃》ノルムカ『小さくな知ていゐことかオ’ た正Ⅱl行列である..

鰍-M

がZJ..

参照

関連したドキュメント

算処理の効率化のliM点において従来よりも優れたモデリング手法について提案した.lMil9f

スキルに国境がないIT系の職種にお いては、英語力のある人材とない人 材の差が大きいので、一定レベル以

2 E-LOCA を仮定した場合でも,ECCS 系による注水流量では足りないほどの原子炉冷却材の流出が考

上位系の対策が必要となる 場合は早期連系は困難 上位系及び配電用変電所の 逆潮流対策等が必要となる

(参考)系統連系希望者がすべて旧費用負担ルール ※4 適用者 ※5 の場合における工事費用 特定負担 約1,310百万円.. ※1

上位系の対策が必要となる 場合は早期連系は困難 上位系及び配電用変電所の 逆潮流対策等が必要となる

 自然科学の場合、実験や観測などによって「防御帯」の

(参考)系統連系希望者がすべて旧費用負担ルール ※4 適用者 ※5 の場合における工事費用 特定負担 約6,740百万円.. ※2