WO2020246531A1

WO2020246531A1 - 情報処理装置、プログラム、及び算出方法

Info

Publication number: WO2020246531A1
Application number: PCT/JP2020/022049
Authority: WO
Inventors: 修一矢作
Original assignee: いすゞ自動車株式会社
Priority date: 2019-06-07
Filing date: 2020-06-04
Publication date: 2020-12-10
Also published as: JP2020201626A

Abstract

制御器Ｃと、制御器Ｃの出力を入力とする制御対象Ｐとを備える閉ループ系の制御システムＳにおいて、制御器Ｃのパラメータを算出する情報処理装置１において、時系列データ取得部３０は、制御器Ｃの出力の時系列データである第１時系列データと、制御対象Ｐの出力の時系列データである第２時系列データとを取得する。入力信号推定部３１は、制御器Ｃのパラメータ、第１時系列データ、及び第２時系列データから、制御器Ｃに入力する入力信号の時系列データの推定値である第３時系列データを推定する。相補感度関数算出部３２は、第２時系列データと第３時系列データとに基づいて、制御器Ｃに対する相補感度関数を算出する。プラント出力算出部３３は、制御システムＳの入力信号を相補感度関数に印加したときの出力である第４時系列データを算出する。

Description

情報処理装置、プログラム、及び算出方法

　本開示は情報処理装置、プログラム、及び算出方法に関し、特に、閉ループ系の制御器の制御パラメータを設定する技術に関する。

　閉ループ系における制御対象のモデルを用いない制御手法が種々提案されている。このようなモデルを用いない制御器パラメータの自動調整手法の一つとしてＦＲＩＴ（Fictitious Reference Iterative Tuning）が知られている（特許文献１参照）。

日本国特開２０１７－１８２６２４号公報

　ＦＲＩＴ等のデータ駆動制御により得られた制御パラメータを用いて閉ループ系を制御すると、系が不安定となる場合がある。特に、ＦＲＩＴで用いられる参照モデルの応答を実現する制御パラメータが存在しない場合に系が不安定となりやすいことが知られている。

　本開示はこれらの点に鑑みてなされたものであり、ＦＲＩＴにおいて制御器パラメータを用いた制御の安定性を評価する技術を提供することを目的とする。

　本開示のある態様は、制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出する情報処理装置である。この装置は、前記制御器の出力の時系列データである第１時系列データと、前記制御対象の出力の時系列データである第２時系列データとを取得する時系列データ取得部と、前記制御器のパラメータ、前記第１時系列データ、及び前記第２時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第３時系列データを推定する入力信号推定部と、前記第２時系列データと前記第３時系列データとに基づいて、前記制御器に対する相補感度関数を算出する相補感度関数算出部と、前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第４時系列データを算出するプラント出力算出部と、を備える。

　前記制御システムは、前記制御対象の入力に対する所定の応答を実現するモデルである参照モデルをさらに備えていてもよく、前記情報処理装置は、前記参照モデルに前記入力信号を入力した場合の出力の時系列データである第５時系列データを取得するモデル出力取得部と、前記第４時系列データと前記第５時系列データとの誤差に関する評価関数の評価値に基づいて前記制御器のパラメータを更新するパラメータ更新部と、をさらに備えてもよい。

　前記評価関数は、前記第４時系列データと前記第５時系列データとの誤差の二乗和であってもよく、前記パラメータ更新部は、前記評価関数の評価値が小さくなるように、反復処理によって前記制御器のパラメータを更新してもよい。

　本開示の別の態様は、プログラムである。このプログラムは、制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出するコンピュータに、前記制御器の出力の時系列データである第１時系列データを取得する機能と、前記制御対象の出力の時系列データである第２時系列データを取得する機能と、前記制御器のパラメータ、前記第１時系列データ、及び前記第２時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第３時系列データを推定する機能と、前記第２時系列データと前記第３時系列データとに基づいて、前記制御器に対する相補感度関数を算出する機能と、前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第４時系列データを算出する機能と、を実現させる。
本開示の別の態様は、算出方法である。この算出方法は、制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出する算出方法であって、前記制御器の出力の時系列データである第１時系列データを取得すること、前記制御対象の出力の時系列データである第２時系列データを取得すること、前記制御器のパラメータ、前記第１時系列データ、及び前記第２時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第３時系列データを推定すること、前記第２時系列データと前記第３時系列データとに基づいて、前記制御器に対する相補感度関数を算出すること、及び前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第４時系列データを算出すること、を含む。
本開示の別の態様は、記憶媒体である。この記憶媒体は、制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出するコンピュータによって読取可能なコンピュータプログラムを記憶する記憶媒体であって、前記コンピュータプログラムは、前記コンピュータにより実行されると、前記コンピュータに以下を実行させる：前記制御器の出力の時系列データである第１時系列データを取得すること、前記制御対象の出力の時系列データである第２時系列データを取得すること、前記制御器のパラメータ、前記第１時系列データ、及び前記第２時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第３時系列データを推定すること、前記第２時系列データと前記第３時系列データとに基づいて、前記制御器に対する相補感度関数を算出すること、及び前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第４時系列データを算出すること。

　なお、以上の構成要素の任意の組み合わせ、本開示の表現を方法、装置、システム、コンピュータプログラム、データ構造、記録媒体などの間で変換したものもまた、本開示の態様として有効である。

　本開示によれば、ＦＲＩＴにおいて制御器パラメータを用いた制御の安定性を評価することができる。

図１は、標準的なＦＲＩＴを説明するための図である。図２は、安定性を考慮したＦＲＩＴの構成を模式的に示す図である。図３Ａは、制御システムにステップ入力を印加したときのシミュレーションの結果を示す図であり、制御システムの出力を示している。図３Ｂは、制御システムにステップ入力を印加したときのシミュレーションの結果を示す図であり、図３Ａに対応する入力信号を示している。図４は、実施の形態に係る情報処理装置の機能構成を模式的に示す図である。図５は、実施の形態に係る情報処理装置が実行する情報処理の流れを説明するためのフローチャートである。

［１．概略］
　ＦＲＩＴや非反証制御等のデータ駆動制御により得られた制御パラメータを用いて制御を行うと、閉ループ系が不安定化する場合がある。特に、制御器の構造上、参照モデルの応答を実現する制御パラメータが存在しない場合に不安定化が生じやすいことが知られている。この問題はSafonovが提案した疑似参照入力を使用する場合に生じる。

　Engellらは、上記の問題が発生する原因として、プラント出力から疑似誤差（疑似参照信号とプラント出力の誤差）の伝達関数を求める際に不安定極が相殺されてしまうため、閉ループ系の不安定化を検知できないことを示している。そこで、Engellらは、疑似参照入力と疑似誤差に関する感度関数を入出力データから求めた後、得られた感度関数に目標値を印加し感度関数の出力である誤差を求めることを提案した。これにより、不安定極を相殺することなく感度関数の出力である誤差を求めることができ、閉ループ系の不安定化を検知することができる。Engellらの手法において、感度関数はＦＩＲ（Finite Impulse Response）モデルに基づいて同定される。このため、感度関数の同定にプラントモデルの構造は不要である。さらに、時間領域で計算していることからオンライン計算への展開が可能となる。

　Engellらは感度関数の出力である誤差を最小にする評価関数を設定した。これに対し、実施の形態に係る手法は、閉ループ系と設計者が設定した参照モデルとが一致するような評価関数を設定する。具体的には、本実施の形態では、まず、疑似参照入力を利用し、調整する制御器に対する相補感度関数を求める。すなわち相補感度関数の出力は制御パラメータの関数となる。次に、設計者が与えたい目標値を相補感度関数に印加しその出力を求める。相補感度関数の出力はプラントからの出力である。すなわち、相補感度関数の出力を見ることにより、ＦＲＩＴにおいて制御器パラメータを用いた制御の安定性を評価することができる。

　実施の形態に係る手法では、相補感度関数から得られた出力と設計者が設定した参照モデルの出力の二乗誤差が最小になる制御器パラメータを粒子群最適化等の最適化手法により求める。これにより、標準的なＦＲＩＴでは閉ループ系が不安定になる場合であっても、本実施の形態では安定な制御パラメータが得られる。これより、実施の形態は標準ＦＲＩＴの利点を活かしたまま不安定化を軽減することができる。

［２．安定性を考慮したＦＲＩＴの導出］［２．１．標準ＦＲＩＴ］
　図１は、標準的なＦＲＩＴを説明するための図である。図１に示す閉ループ系の制御システムＳにおいて、制御器Ｃは制御に用いるパラメータθを引数とする関数Ｃ（θ）で表現されている。図１に示す制御システムＳは、制御対象Ｐの出力ｙと後述する参照モデルＭの出力とを一致させることが目的である。具体的には、この目的を達成するために制御対象Ｐに入力すべき制御量ｕを出力するようなパラメータθを特定することが目的である。なお、図１において、θは自由に調整できるパラメータであり、ｕ及びｙは観測により取得可能なデータである。また、ｄは制御システムＳの目標値である。

　ＦＲＩＴは１組の入出力データと参照モデルＭとから、制御システムＳの制御器Ｃのパラメータθを自動調整する手法である。初期のパラメータθを用いて１組の閉ループ実験を行い、そのときの入出力データｕ及びｙをサンプリングして計測する。このとき、制御システムＳは安定であるとする。図１より、Ｃ（θ）（ｒ－ｙ）＝ｕであるから、制御システムＳの目標値ｄの推定値である疑似参照信号ｒはｕ及びｙを用いて以下の式（１）で表せる。
　　ｒ（θ，ｋ）＝Ｃ^-１（θ）ｕ（ｋ）＋ｙ（ｋ）　　（１）

　図１に示す一般的なフィードバック制御の応答と、参照モデルＭ（ｚ）及び擬似参照信号ｒ（θ，ｋ）から得られる目標応答との誤差に関する評価関数Ｊ_θは次式（２）で表される。

　この評価関数Ｊ_θを最小にする制御器Ｃのパラメータθは、制御対象Ｐの出力であるプラント出力ｙ（ｋ）と参照モデルＭの出力Ｍ（ｚ）ｒ（θ，ｋ）との二乗誤差を最小化するという意味において、制御器Ｃの最適なパラメータθである。一般的なＦＲＩＴでは、オフライン計算で最適なパラメータθが計算される。なお、評価関数は式（２）に示す形に限られず、制御入力等の制約を考慮したものであってもよい。

　式（１）から明らかなように、ＦＲＩＴは閉ループ系の制御システムＳの伝達関数と参照モデルがマッチングする最適な制御器パラメータを求めることを目的としている。すなわち、ＦＲＩＴは、以下の式（２）で表される評価関数を最小にする最適パラメータを求める。ＦＲＩＴは擬似参照信号ｒ（θ，ｋ）を利用して、閉ループ試験を繰り返すことなく実験により取得した１組の入出力データを用いてオフラインで最適制御パラメータを求めるデータ駆動制御の１つであるともいえる。

　ＦＲＩＴの評価関数Ｊ_θ、疑似参照信号ｒ（θ，ｋ）、及び疑似誤差ｅ（θ，ｋ）は次式で表される。

［２．２．安定性を考慮したＦＲＩＴ］
　疑似参照信号ｒ（θ，ｋ）を用いる場合にはＦＲＩＴの閉ループ系の不安定化を検知できない。そこで、まず、疑似参照信号ｒ（θ，ｋ）とプラント出力ｙ_０を用いて時間領域における相補感度関数を求める。求めた相補感度関数に目標値ｄを印加し、その応答ｙ^＊を求める。この応答ｙ^＊が参照モデルＭの応答Ｍ（ｚ）ｒ（θ，ｋ）と一致するような制御器パラメータを求める。なお、相補感度関数の同定はＦＩＲモデルを用いる。これにより、制御対象Ｐの構造を知る必要はなく、取得したデータのみを用いて相補感度関数を同定することができる。

　まず、ラプラス領域における感度関数と相補感度関数とについて説明する。
　パラメータθの調整対象とする制御器Ｃの感度関数Ｓは次式（６）で表される。

　調整対象とする制御器Ｃの相補感度関数Ｔは次式（７）で表される。

　なお、感度関数Ｓと相補感度関数Ｔとの関係は次式（８）となる。
　　Ｓ（ｓ）＋Ｔ（ｓ）＝１　　（８）

　続いて、時間領域における疑似参照信号ｒ（θ，ｋ）とプラント出力ｙ_０との関係について説明する。時間領域における疑似参照信号ｒ（θ，ｋ）とプラント出力ｙ_０との関係は次式（９）となる。

　式（９）において、記号＊は畳み込みを表し、ｔ（ｋ）は相補感度関数Ｔのインパルス応答を表す。疑似参照信号ｒ（θ，ｋ）とプラント出力ｙ_０とは観測可能であるが、相補感度関数Ｔのインパルス応答ｔ（ｋ）は未知である。

　式（９）は行列を用いて表すと以下の式（１０）となる。

　式（１０）の左辺をベクトルｙ_０、右辺第１項を行列Ｒ_θ、右辺第２項をベクトルｔとすると、相補感度関数Ｔのインパルス応答ｔ（ｋ）は、以下の式（１１）で表される。

　式（１１）は、式（９）のデコンボリューションともいえる。ｔは、パラメータθに依存する。

　制御システムＳの入力信号ｄ（目標値）を相補感度関数Ｔに印加した場合の出力ｙ^＊を時間領域で表すと、次式（１２）となる。

　式（１２）は行列を用いて表すと以下の式（１３）となる。

　式（１３）の左辺をベクトルｙ^＊、右辺第１項を行列Ｄ、右辺第２項をベクトルｔとし、式（１１）を用いてｔを消去すると、式（１３）は式（１４）に変形できる。

　式（１４）の右辺は全て観測により取得可能であるから、制御器Ｃを用いたときの目標値ｄを閉ループ系の制御システムＳに印加したときの出力ベクトルｙ^＊を演算により取得することができる。この出力ベクトルｙ^＊と参照モデルＭの出力ｙ_ｄ（ｋ）＝Ｍ（ｚ）ｄ（ｋ）との誤差ｅ^＊の二乗和を、次式（１５）で表される評価関数Ｊ^＊（θ）とする。

　式（１５）において、
　ｅ^＊（θ，ｋ）＝ｙ^＊（θ，ｋ）－ｙ_ｄ（ｋ）　　（１６）である。

　標準的なＦＲＩＴでは、取得したプラント出力ｙ_０に一致するようにパラメータθの関数である疑似参照信号ｒ（θ，ｋ）を調整する。これに対し、実施の形態に係る手法では、参照モデルＭの出力Ｍ（ｚ）ｄ（ｋ）にプラント出力ｙ_０が一致するように、パラメータθを調整してプラント出力ｙ^＊を変更する。すなわち、標準的なＦＲＩＴは実験によりあらかじめ取得したプラント出力に合うようにパラメータθを求めているのに対し、実施の形態に係る手法は参照モデルの出力Ｍ（ｚ）ｄ（ｋ）に合うように、パラメータθを求める。

［２．３．安定性を考慮したＦＲＩＴの構成］
　図２は、安定性を考慮したＦＲＩＴの構成を模式的に示す図である。まず、パラメータθが初期値θ_０のときの制御システムＳの入出力データから疑似参照信号ｒ（θ）が算出される。算出された疑似参照入力とプラント出力ｙとを用いて相補感度関数ｔを上記式（１１）を用いて求められ、制御システムＳの入力信号である目標値ｄを相補感度関数ｔに入力される。相補感度関数ｔの出力である出力ｙ^＊と、参照モデルの出力Ｍ（ｚ）ｄとの誤差を最小にする制御器Ｃのパラメータθが、最適化手法により求められる。

［３．シミュレーションによる検証］［３．１．システムの定式化］
　本願の発明者は、安定性を考慮したＦＲＩＴについてシミュレーションを用いて検証した。そこで、まず、シミュレーションで用いる制御対象Ｐ、参照モデルＭ、制御器Ｃの定式化を行う。

　制御対象Ｐは次式（１７）に示すように、ばね質量系とする。ここで、ｍ、ｃ、ｋ、及びＬは、それぞれ質量、粘性係数、ばね剛性、及びむだ時間を表す。

　参照モデルＭは、次式（１８）に示すように１次遅れ系とする。ここで、τ_Ｍはシステムの応答に関するパラメータであり、設計者が設定する。

　制御器Ｃは次式（１９）に示すＰＩＤ（Proportional-Integral-Differential）制御器を用いる。ここで、Ｋ_ｐ、Ｋ_ｉ、及びＫ_ｄはそれぞれ、Ｐゲイン、Ｉゲイン、及びＤゲインである。

［３．２．シミュレーション］
　シミュレーションにおいてサンプリング周期は４ミリ秒とし、制御器Ｃの離散化はゼロ次ホールドを用いた。なお、最適化計算には粒子群最適化を用いた。

　図３Ａと図３Ｂは、従来のＦＲＩＴと実施の形態に係る手法とのそれぞれにより調整したパラメータθによる制御システムＳの入出力データの時間履歴を示す図である。具体的には、図３Ａにおいて、破線は参照応答、実線は実施の形態に係る手法における制御システムＳの出力、一点鎖線は従来のＦＲＩＴにおける制御システムＳの出力を示している。図３Ａより、標準ＦＲＩＴでは制御システムＳの出力が発散してグラフの枠内に収まっていないのに対し、実施の形態に係る手法は、制御の構造上モデルマッチングが実現できない場合においても閉ループ系である制御システムＳが安定となる制御器Ｃのパラメータθが得られていることがわかる。

　図３Ｂは、図３Ａに対応する入力信号を示す図である。図３Ａと同様に、実線は実施の形態に係る手法における入力、一点鎖線は従来のＦＲＩＴにおける入力を示している。図３Ｂより、標準ＦＲＩＴでは入力も不安定となっている。一方、実施の形態に係る手法では、入力も安定している。

［４．小括］
　従来のＦＲＩＴでは、制御器Ｃのパラメータθの最適解を用いて制御しても、閉ループ系の挙動が不安定になるという問題があった。実施の形態に係る手法では、時間領域で擬似参照信号ｒ（θ，ｋ）とプラント出力ｙとの伝達関数である相補感度関数Ｔを同定することで、誤差を正しく検知することを可能にした。実施の形態に係る手法により、閉ループ系が安定である制御器Ｃのパラメータθを求められる。

［５．参考文献］
（ＦＲＩＴ）・相馬将太郎, 金子修, 藤井隆雄, 一回の実験データに基づく制御器パラメータチューニングの新しいアプローチ Fictitious Reference Iterative Tuning の提案, システム制御情報学会論文誌, Vol. 17, No.12 (2004), pp. 528-536・奥谷明大, 金子修, 山本茂, ＦＲＩＴを用いた多入出力むだ時間系に対するスミス補償器のチューニング, システム制御情報学会論文誌, Vol. 28, No 2 (2015), pp. 58-65
・データを直接用いた制御器のパラメータチューニング, 金子修, 計測と制御, Vol.43, No.11 (2008), pp903-908

（非反証制御）・M. G. Safonov and T. C. Tsao, The unfalsified control, concept and learning, IEEE Trans. on Automat.Contr., Vol. 42, No. 6, pp. 843-847 (1997)

（安定性の考慮）・弓場井一裕, 藤井宏樹, 平井淳之, パラメータ更新時の閉ループシステムの安定性を考慮したFCbTの提案, 電気学会論文誌Ｄ（産業応用部門誌）, Vol.132, No.6 (2011),pp. 607-615・Kazuhiro Yubai, Hiroki Fujii, Junji Hirai, Fictitious Correlation-based Tuning Integrating the Data-Based Stability Test at Each Parameter Update, Electrical Power Systems and Computers, LNEE 99, pp. 511-518.

（疑似参照信号の問題と安定性）・S. Engell, T. Tometzki and T. Wonghong, A New Approach to Adaptive Unfalsified Control. In Proc. European Control Conf., Kos, 2007, 1328-1333.・T. Wonghong and S. Engell, Application of a New Scheme for Adaptive Unfalsified Control to a CSTR. Proc. IFAC World Congress, Korea, 13247-13252, 2008.

［６．実施の形態］
　以上を踏まえ、本開示の実施の形態について説明する。

　図４は、実施の形態に係る情報処理装置１の機能構成を模式的に示す図である。情報処理装置１は、記憶部２と制御部３とを備える。図４において、矢印は主なデータの流れを示しており、図４に示していないデータの流れがあってもよい。図４において、各機能ブロックはハードウェア（装置）単位の構成ではなく、機能単位の構成を示している。そのため、図４に示す機能ブロックは単一の装置内に実装されてもよく、あるいは複数の装置内に分かれて実装されてもよい。機能ブロック間のデータの授受は、データバス、ネットワーク、可搬記憶媒体等、任意の手段を介して行われてもよい。

　記憶部２は、情報処理装置１を実現するコンピュータのＢＩＯＳ（Basic Input Output System）等を格納するＲＯＭ（Read Only Memory）や情報処理装置１の作業領域となるＲＡＭ（Random Access Memory）、ＯＳ（Operating System）やアプリケーションプログラム、当該アプリケーションプログラムの実行時に参照される種々の情報を格納するＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）等の大容量記憶装置である。

　制御部３は、情報処理装置１のＣＰＵ（Central Processing Unit）やＧＰＵ（Graphics Processing Unit）等のプロセッサであり、記憶部２に記憶されたプログラムを実行することによって時系列データ取得部３０、入力信号推定部３１、相補感度関数算出部３２、プラント出力算出部３３、モデル出力取得部３４、及びパラメータ更新部３５として機能する。

　なお、図４は、情報処理装置１が単一の装置で構成されている場合の例を示している。しかしながら、情報処理装置１は、例えばクラウドコンピューティングシステムのように複数のプロセッサやメモリ等の計算リソースによって実現されてもよい。この場合、制御部３を構成する各部は、複数の異なるプロセッサの中の少なくともいずれかのプロセッサがプログラムを実行することによって実現される。

　情報処理装置１は、制御器Ｃと制御器Ｃの出力を入力とする制御対象Ｐとを備える制御システムＳにおいて、制御器Ｃのパラメータθを算出するための装置である。図４に示すように、制御システムＳは、制御対象Ｐの出力が制御器Ｃの入力にフィードバックされる閉ループ系である。

　時系列データ取得部３０は、制御器Ｃの出力の時系列データである第１時系列データと、制御対象Ｐの出力の時系列データである第２時系列データとを取得する。ここで、第１時系列データは上述した制御対象Ｐに入力すべき制御量ｕに対応し、第２時系列データは上述した制御対象Ｐの出力ｙに対応する。したがって、以下本明細書において、「第１時系列データｕ」、「第２時系列データｙ」と記載することがある。

　入力信号推定部３１は、制御器Ｃのパラメータθ、第１時系列データｕ、及び第２時系列データｙから、制御器Ｃに入力する入力信号の時系列データの推定値である第３時系列データを推定する。第３時系列データは、上述した擬似参照信号ｒ（θ，ｋ）に対応する。したがって、以下、「第３時系列データｒ（θ）」と記載することがある。

　相補感度関数算出部３２は、第２時系列データｕと第３時系列データｒ（θ）とに基づいて、上述した式（１１）を用いて制御器Ｃに対する相補感度関数ｔを算出する。プラント出力算出部３３は、上述した式（１２）又は式（１３）を用いて、制御システムＳの入力信号ｄを相補感度関数ｔに印加したときの出力である第４時系列データを算出する。第４時系列データは、制御システムＳの入力信号ｄ（目標値）を相補感度関数Ｔに印加した場合の出力ｙ^＊に対応する。したがって、以下、「第４時系列データｙ^＊」と記載することがある。

　制御対象Ｐに入力すべき制御量である第１時系列データｕ、制御対象Ｐの出力である第２時系列データｙ、及び制御システムＳの入力信号である第４時系列データｄは、いずれも観測により取得可能な量である。情報処理装置１は、観測により取得した第１時系列データｕ、第２時系列データｙ、及び入力信号ｄを用いることで、制御システムＳに入力信号ｄを入力したときの出力ｙ^＊を算出することができる。出力ｙ^＊は、プラント（制御対象Ｐ）からの出力である。情報処理装置１は、この出力ｙ^＊の振る舞い（例えば、発散するか否か、振動するか否か、収束するか否か等）を解析することで、ＦＲＩＴにおいて制御器Ｃのパラメータθを用いたプラントの制御の安定性を評価することができる。

　図２及び図４に示すように、制御システムＳは制御対象Ｐの入力信号ｄに対する出力を所定の出力とすることを実現する参照モデルＭを備えている。なお、参照モデルＭは、制御対象Ｐの出力が設計者が希望する出力となるように、設計者が決定する。モデル出力取得部３４は、参照モデルＭに入力信号ｄを入力した場合の出力の時系列データである第５時系列データを取得する。第５時系列データは、上述した参照モデルの出力Ｍ（ｚ）ｄに対応する。以下、第５時系列データを「第５時系列データｙ_ｄ」と記載することがある。

　ここで、情報処理装置１がＦＲＩＴにおける制御器Ｃのパラメータθを用いた制御の安定性を評価することができるのであれば、情報処理装置１は、制御が安定化するようにパラメータθを最適化することもできる。これを実現するために、パラメータ更新部３５は、第４時系列データｙ^＊と第５時系列データｙ_ｄとの誤差に関する評価関数Ｊ^＊（θ）の評価値に基づいて、制御器Ｃのパラメータθを更新する。

　より具体的には、上述した式（１５）及び式（１６）に示すように、パラメータ更新部３５が用いる評価関数Ｊ^＊（θ）は、第４時系列データｙ^＊と第５時系列データｙ_ｄとの誤差の二乗和である。パラメータ更新部３５は、評価関数Ｊ^＊（θ）の評価値が小さくなるように、反復処理によって制御器Ｃのパラメータθを更新する。すなわち、パラメータ更新部３５は、第４時系列データｙ^＊と第５時系列データｙ_ｄとの誤差の二乗和が小さくなるという意味において最適なパラメータθを反復によって求める。

　一般に、第４時系列データｙ^＊が発散したり振動したりすると、第４時系列データｙ^＊と第５時系列データｙ_ｄとの誤差の二乗和は大きくなる。パラメータ更新部３５が第４時系列データｙ^＊と第５時系列データｙ_ｄとの誤差の二乗和が小さくなるように制御器Ｃのパラメータθを更新することにより、情報処理装置１は、ＦＲＩＴにおいて制御器Ｃのパラメータθを用いたプラントの制御を安定化させることができる。

　なお、パラメータ更新部３５は、評価関数の評価値が小さくなるように制御器Ｃのパラメータθを更新できるのであればどのような最適化手法を用いてもよい。一例としてパラメータ更新部３５は、あらかじめ定められた所定の回数を反復回数の上限として、粒子群最適化の手法を用いて制御器Ｃのパラメータθを更新してもよい。

＜情報処理装置１が実行する情報処理方法の処理フロー＞
　図５は、実施の形態に係る情報処理装置１が実行する情報処理の流れを説明するためのフローチャートである。本フローチャートにおける処理は、例えば情報処理装置１が起動したときに開始する。

　時系列データ取得部３０は、制御器Ｃの出力の時系列データである第１時系列データｕを取得する（Ｓ２）。また、時系列データ取得部３０は、制御対象Ｐの出力の時系列データである第２時系列データｙを取得する（Ｓ４）。

　入力信号推定部３１は、制御器Ｃのパラメータθ、第１時系列データｕ、及び第２時系列データｙから、式（１）を用いて、制御器Ｃに入力する入力信号ｄの時系列データの推定値である第３時系列データｒ（θ）を推定する（Ｓ６）。

　相補感度関数算出部３２は、第２時系列データｙと第３時系列データｒ（θ）とに基づいて、式（１１）を用いて制御器Ｃに対する相補感度関数ｔ（すなわち、擬似参照信号ｒ（θ，ｋ）に対する相補感度関数ｔ）を算出する（Ｓ８）。プラント出力算出部３３は、式（１４）を用いて、制御システムＳの入力信号ｄを相補感度関数ｔに印加したときの出力である第４時系列データｙ^＊（すなわち、入力信号ｄに対する制御対象Ｐの出力）を算出する（Ｓ１０）。

　モデル出力取得部３４は、参照モデルＭに入力信号ｄを入力した場合の出力の時系列データである第５時系列データｙ_ｄを取得する（Ｓ１２）。パラメータ更新部３５は、式（１５）及び式（１６）を用いて、第４時系列データｙ^＊と第５時系列データｙ_ｄとの誤差に関する評価関数Ｊ（θ）の評価値を算出する（Ｓ１４）。パラメータ更新部３５は、評価関数Ｊ（θ）の評価値が小さくなるように、反復処理によって制御器Ｃのパラメータθを更新する（Ｓ１６）。

　パラメータ更新部３５が制御器Ｃのパラメータθを更新すると、本フローチャートにおける処理は終了する。情報処理装置１は、上記の処理をオンラインで繰り返すことにより、制御器Ｃのパラメータθの更新を継続する。

＜実施の形態に係る情報処理装置１が奏する効果＞
　以上説明したように、実施の形態に係る情報処理装置１によれば、ＦＲＩＴにおいて制御器Ｃのパラメータθを用いたプラントの制御の安定性を評価することができる。また、情報処理装置１は、閉ループ系の制御が安定化するようにＦＲＩＴにおける制御器Ｃパラメータθを最適化することもできる。

　以上、本開示を実施の形態を用いて説明したが、本開示の技術的範囲は上記実施の形態に記載の範囲には限定されず、その要旨の範囲内で種々の変形及び変更が可能である。例えば、装置の全部又は一部は、任意の単位で機能的又は物理的に分散・統合して構成することができる。また、複数の実施の形態の任意の組み合わせによって生じる新たな実施の形態も、本開示の実施の形態に含まれる。組み合わせによって生じる新たな実施の形態の効果は、もとの実施の形態の効果を併せ持つ。

　本出願は、2019年6月7日付で出願された日本国特許出願（特願2019-106908）に基づくものであり、その内容はここに参照として取り込まれる。

　本開示の情報処理装置、及びプログラム、及び算出方法は、ＦＲＩＴにおいて制御器パラメータを用いた制御の安定性を評価する技術を提供することができる、という点において有用である。

１・・・情報処理装置
２・・・記憶部
３・・・制御部
３０・・・時系列データ取得部
３１・・・入力信号推定部
３２・・・相補感度関数算出部
３３・・・プラント出力算出部
３４・・・モデル出力取得部
３５・・・パラメータ更新部
Ｃ・・・制御器
Ｍ・・・参照モデル
Ｐ・・・制御対象
Ｓ・・・制御システム

Claims

　制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出する情報処理装置であって、
　前記制御器の出力の時系列データである第１時系列データと、前記制御対象の出力の時系列データである第２時系列データとを取得する時系列データ取得部と、
　前記制御器のパラメータ、前記第１時系列データ、及び前記第２時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第３時系列データを推定する入力信号推定部と、
　前記第２時系列データと前記第３時系列データとに基づいて、前記制御器に対する相補感度関数を算出する相補感度関数算出部と、
　前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第４時系列データを算出するプラント出力算出部と、
　を備える情報処理装置。
　前記制御システムは、前記制御対象の入力に対する所定の応答を実現するモデルである参照モデルをさらに備えており、
　前記情報処理装置は、
　　前記参照モデルに前記入力信号を入力した場合の出力の時系列データである第５時系列データを取得するモデル出力取得部と、
　　前記第４時系列データと前記第５時系列データとの誤差に関する評価関数の評価値に基づいて前記制御器のパラメータを更新するパラメータ更新部と、
　をさらに備える請求項１に記載の情報処理装置。
　前記評価関数は、前記第４時系列データと前記第５時系列データとの誤差の二乗和であり、
　前記パラメータ更新部は、前記評価関数の評価値が小さくなるように、反復処理によって前記制御器のパラメータを更新する、
　請求項２に記載の情報処理装置。
　制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出するコンピュータに、
　前記制御器の出力の時系列データである第１時系列データを取得する機能と、
　前記制御対象の出力の時系列データである第２時系列データを取得する機能と、
　前記制御器のパラメータ、前記第１時系列データ、及び前記第２時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第３時系列データを推定する機能と、
　前記第２時系列データと前記第３時系列データとに基づいて、前記制御器に対する相補感度関数を算出する機能と、
　前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第４時系列データを算出する機能と、
　を実現させるプログラム。
　制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出する算出方法であって、
　前記制御器の出力の時系列データである第１時系列データを取得すること、
　前記制御対象の出力の時系列データである第２時系列データを取得すること、
　前記制御器のパラメータ、前記第１時系列データ、及び前記第２時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第３時系列データを推定すること、
　前記第２時系列データと前記第３時系列データとに基づいて、前記制御器に対する相補感度関数を算出すること、及び
　前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第４時系列データを算出すること、
　を含む算出方法。