JP2020201626A - 情報処理装置及びプログラム - Google Patents
情報処理装置及びプログラム Download PDFInfo
- Publication number
- JP2020201626A JP2020201626A JP2019106908A JP2019106908A JP2020201626A JP 2020201626 A JP2020201626 A JP 2020201626A JP 2019106908 A JP2019106908 A JP 2019106908A JP 2019106908 A JP2019106908 A JP 2019106908A JP 2020201626 A JP2020201626 A JP 2020201626A
- Authority
- JP
- Japan
- Prior art keywords
- series data
- time series
- output
- controller
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
【課題】FRITにおいて制御器パラメータを用いた制御の安定性を評価する。【解決手段】制御器Cと、制御器Cの出力を入力とする制御対象Pとを備える閉ループ系の制御システムSにおいて、制御器Cのパラメータを算出する情報処理装置1において、時系列データ取得部30は、制御器Cの出力の時系列データである第1時系列データと、制御対象Pの出力の時系列データである第2時系列データとを取得する。入力信号推定部31は、制御器Cのパラメータ、第1時系列データ、及び第2時系列データから、制御器Cに入力する入力信号の時系列データの推定値である第3時系列データを推定する。相補感度関数算出部32は、第2時系列データと第3時系列データとに基づいて、制御器Cに対する相補感度関数を算出する。プラント出力算出部33は、制御システムSの入力信号を相補感度関数に印加したときの出力である第4時系列データを算出する。【選択図】図4
Description
本発明は情報処理装置及びプログラムに関し、特に、閉ループ系の制御器の制御パラメータを設定する技術に関する。
閉ループ系における制御対象のモデルを用いない制御手法が種々提案されている。このようなモデルを用いない制御器パラメータの自動調整手法の一つとしてFRIT(Fictitious Reference Iterative Tuning)が知られている(特許文献1参照)。
FRIT等のデータ駆動制御により得られた制御パラメータを用いて閉ループ系を制御すると、系が不安定となる場合がある。特に、FRITで用いられる参照モデルの応答を実現する制御パラメータが存在しない場合に系が不安定となりやすいことが知られている。
本発明はこれらの点に鑑みてなされたものであり、FRITにおいて制御器パラメータを用いた制御の安定性を評価する技術を提供することを目的とする。
本発明のある態様は、制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出する情報処理装置である。この装置は、前記制御器の出力の時系列データである第1時系列データと、前記制御対象の出力の時系列データである第2時系列データとを取得する時系列データ取得部と、前記制御器のパラメータ、前記第1時系列データ、及び前記第2時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第3時系列データを推定する入力信号推定部と、前記第2時系列データと前記第3時系列データとに基づいて、前記制御器に対する相補感度関数を算出する相補感度関数算出部と、前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第4時系列データを算出するプラント出力算出部と、を備える。
前記制御システムは、前記制御対象の入力に対する所定の応答を実現するモデルである参照モデルをさらに備えていてもよく、前記情報処理装置は、前記参照モデルに前記入力信号を入力した場合の出力の時系列データである第5時系列データを取得するモデル出力取得部と、前記第4時系列データと前記第5時系列データとの誤差に関する評価関数の評価値に基づいて前記制御器のパラメータを更新するパラメータ更新部と、をさらに備えてもよい。
前記評価関数は、前記第4時系列データと前記第5時系列データとの誤差の二乗和であってもよく、前記パラメータ更新部は、前記評価関数の評価値が小さくなるように、反復処理によって前記制御器のパラメータを更新してもよい。
本発明の別の態様は、プログラムである。このプログラムは、制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出するコンピュータに、前記制御器の出力の時系列データである第1時系列データを取得する機能と、前記制御対象の出力の時系列データである第2時系列データを取得する機能と、前記制御器のパラメータ、前記第1時系列データ、及び前記第2時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第3時系列データを推定する機能と、前記第2時系列データと前記第3時系列データとに基づいて、前記制御器に対する相補感度関数を算出する機能と、前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第4時系列データを算出する機能と、を実現させる。
なお、以上の構成要素の任意の組み合わせ、本発明の表現を方法、装置、システム、コンピュータプログラム、データ構造、記録媒体などの間で変換したものもまた、本発明の態様として有効である。
本発明によれば、FRITにおいて制御器パラメータを用いた制御の安定性を評価することができる。
[1.概略]
FRITや非反証制御等のデータ駆動制御により得られた制御パラメータを用いて制御を行うと、閉ループ系が不安定化する場合がある。特に、制御器の構造上、参照モデルの応答を実現する制御パラメータが存在しない場合に不安定化が生じやすいことが知られている。この問題はSafonovが提案した疑似参照入力を使用する場合に生じる。
FRITや非反証制御等のデータ駆動制御により得られた制御パラメータを用いて制御を行うと、閉ループ系が不安定化する場合がある。特に、制御器の構造上、参照モデルの応答を実現する制御パラメータが存在しない場合に不安定化が生じやすいことが知られている。この問題はSafonovが提案した疑似参照入力を使用する場合に生じる。
Engellらは、上記の問題が発生する原因として、プラント出力から疑似誤差(疑似参照信号とプラント出力の誤差)の伝達関数を求める際に不安定極が相殺されてしまうため、閉ループ系の不安定化を検知できないことを示している。そこで、Engellらは、疑似参照入力と疑似誤差に関する感度関数を入出力データから求めた後、得られた感度関数に目標値を印加し感度関数の出力である誤差を求めることを提案した。これにより、不安定極を相殺することなく感度関数の出力である誤差を求めることができ、閉ループ系の不安定化を検知することができる。Engellらの手法において、感度関数はFIR(Finite Impulse Response)モデルに基づいて同定される。このため、感度関数の同定にプラントモデルの構造は不要である。さらに、時間領域で計算していることからオンライン計算への展開が可能となる。
Engellらは感度関数の出力である誤差を最小にする評価関数を設定した。これに対し、実施の形態に係る手法は、閉ループ系と設計者が設定した参照モデルとが一致するような評価関数を設定する。具体的には、本実施の形態では、まず、疑似参照入力を利用し、調整する制御器に対する相補感度関数を求める。すなわち相補感度関数の出力は制御パラメータの関数となる。次に、設計者が与えたい目標値を相補感度関数に印加しその出力を求める。相補感度関数の出力はプラントからの出力である。すなわち、相補感度関数の出力を見ることにより、FRITにおいて制御器パラメータを用いた制御の安定性を評価することができる。
実施の形態に係る手法では、相補感度関数から得られた出力と設計者が設定した参照モデルの出力の二乗誤差が最小になる制御器パラメータを粒子群最適化等の最適化手法により求める。これにより、標準的なFRITでは閉ループ系が不安定になる場合であっても、本実施の形態では安定な制御パラメータが得られる。これより、実施の形態は標準FRITの利点を活かしたまま不安定化を軽減することができる。
[2.安定性を考慮したFRITの導出]
[2.1.標準FRIT]
図1は、標準的なFRITを説明するための図である。図1に示す閉ループ系の制御システムSにおいて、制御器Cは制御に用いるパラメータθを引数とする関数C(θ)で表現されている。図1に示す制御システムSは、制御対象Pの出力yと後述する参照モデルMの出力とを一致させることが目的である。具体的には、この目的を達成するために制御対象Pに入力すべき制御量uを出力するようなパラメータθを特定することが目的である。なお、図1において、θは自由に調整できるパラメータであり、u及びyは観測により取得可能なデータである。また、dは制御システムSの目標値である。
[2.1.標準FRIT]
図1は、標準的なFRITを説明するための図である。図1に示す閉ループ系の制御システムSにおいて、制御器Cは制御に用いるパラメータθを引数とする関数C(θ)で表現されている。図1に示す制御システムSは、制御対象Pの出力yと後述する参照モデルMの出力とを一致させることが目的である。具体的には、この目的を達成するために制御対象Pに入力すべき制御量uを出力するようなパラメータθを特定することが目的である。なお、図1において、θは自由に調整できるパラメータであり、u及びyは観測により取得可能なデータである。また、dは制御システムSの目標値である。
FRITは1組の入出力データと参照モデルMとから、制御システムSの制御器Cのパラメータθを自動調整する手法である。初期のパラメータθを用いて1組の閉ループ実験を行い、そのときの入出力データu及びyをサンプリングして計測する。このとき、制御システムSは安定であるとする。図1より、C(θ)(r−y)=uであるから、制御システムSの目標値dの推定値である疑似参照信号rはu及びyを用いて以下の式(1)で表せる。
r(θ,k)=C-1(θ)u(k)+y(k) (1)
r(θ,k)=C-1(θ)u(k)+y(k) (1)
図1に示す一般的なフィードバック制御の応答と、参照モデルM(z)及び擬似参照信号r(θ,k)から得られる目標応答との誤差に関する評価関数Jθは次式(2)で表される。
この評価関数Jθを最小にする制御器Cのパラメータθは、制御対象Pの出力であるプラント出力y(k)と参照モデルMの出力M(z)r(θ,k)との二乗誤差を最小化するという意味において、制御器Cの最適なパラメータθである。一般的なFRITでは、オフライン計算で最適なパラメータθが計算される。なお、評価関数は式(2)に示す形に限られず、制御入力等の制約を考慮したものであってもよい。
式(1)から明らかなように、FRITは閉ループ系の制御システムSの伝達関数と参照モデルがマッチングする最適な制御器パラメータを求めることを目的としている。すなわち、FRITは、以下の式(2)で表される評価関数を最小にする最適パラメータを求める。FRITは擬似参照信号r(θ,k)を利用して、閉ループ試験を繰り返すことなく実験により取得した1組の入出力データを用いてオフラインで最適制御パラメータを求めるデータ駆動制御の1つであるともいえる。
FRITの評価関数Jθ、疑似参照信号r(θ,k)、及び疑似誤差e(θ,k)は次式で表される。
[2.2.安定性を考慮したFRIT]
疑似参照信号r(θ,k)を用いる場合にはFRITの閉ループ系の不安定化を検知できない。そこで、まず、疑似参照信号r(θ,k)とプラント出力y0を用いて時間領域における相補感度関数を求める。求めた相補感度関数に目標値dを印加し、その応答y*を求める。この応答y*が参照モデルMの応答M(z)r(θ,k)と一致するような制御器パラメータを求める。なお、相補感度関数の同定はFIRモデルを用いる。これにより、制御対象Pの構造を知る必要はなく、取得したデータのみを用いて相補感度関数を同定することができる。
疑似参照信号r(θ,k)を用いる場合にはFRITの閉ループ系の不安定化を検知できない。そこで、まず、疑似参照信号r(θ,k)とプラント出力y0を用いて時間領域における相補感度関数を求める。求めた相補感度関数に目標値dを印加し、その応答y*を求める。この応答y*が参照モデルMの応答M(z)r(θ,k)と一致するような制御器パラメータを求める。なお、相補感度関数の同定はFIRモデルを用いる。これにより、制御対象Pの構造を知る必要はなく、取得したデータのみを用いて相補感度関数を同定することができる。
まず、ラプラス領域における感度関数と相補感度関数とについて説明する。
パラメータθの調整対象とする制御器Cの感度関数Sは次式(6)で表される。
パラメータθの調整対象とする制御器Cの感度関数Sは次式(6)で表される。
調整対象とする制御器Cの相補感度関数Tは次式(7)で表される。
なお、感度関数Sと相補感度関数Tとの関係は次式(8)となる。
S(s)+T(s)=1 (8)
S(s)+T(s)=1 (8)
続いて、時間領域における疑似参照信号r(θ,k)とプラント出力y0との関係について説明する。時間領域における疑似参照信号r(θ,k)とプラント出力y0との関係は次式(9)となる。
式(9)は行列を用いて表すと以下の式(10)となる。
式(10)の左辺をベクトルy0、右辺第1項を行列Rθ、右辺第2項をベクトルtとすると、相補感度関数Tのインパルス応答t(k)は、以下の式(11)で表される。
制御システムSの入力信号d(目標値)を相補感度関数Tに印加した場合の出力y*を時間領域で表すと、次式(12)となる。
式(12)は行列を用いて表すと以下の式(13)となる。
式(13)の左辺をベクトルy*、右辺第1項を行列D、右辺第2項をベクトルtとし、式(11)を用いてtを消去すると、式(13)は式(14)に変形できる。
式(14)の右辺は全て観測により取得可能であるから、制御器Cを用いたときの目標値dを閉ループ系の制御システムSに印加したときの出力ベクトルy*を演算により取得することができる。この出力ベクトルy*と参照モデルMの出力yd(k)=M(z)d(k)との誤差e*の二乗和を、次式(15)で表される評価関数J*(θ)とする。
標準的なFRITでは、取得したプラント出力y0に一致するようにパラメータθの関数である疑似参照信号r(θ,k)を調整する。これに対し、実施の形態に係る手法では、参照モデルMの出力M(z)d(k)にプラント出力y0が一致するように、パラメータθを調整してプラント出力y*を変更する。すなわち、標準的なFRITは実験によりあらかじめ取得したプラント出力に合うようにパラメータθを求めているのに対し、実施の形態に係る手法は参照モデルの出力M(z)d(k)に合うように、パラメータθを求める。
[2.3.安定性を考慮したFRITの構成]
図2は、安定性を考慮したFRITの構成を模式的に示す図である。まず、パラメータθが初期値θ0のときの制御システムSの入出力データから疑似参照信号r(θ)が算出される。算出された疑似参照入力とプラント出力yとを用いて相補感度関数tを上記式(11)を用いて求められ、制御システムSの入力信号である目標値dを相補感度関数tに入力される。相補感度関数tの出力である出力y*と、参照モデルの出力M(z)dとの誤差を最小にする制御器Cのパラメータθが、最適化手法により求められる。
図2は、安定性を考慮したFRITの構成を模式的に示す図である。まず、パラメータθが初期値θ0のときの制御システムSの入出力データから疑似参照信号r(θ)が算出される。算出された疑似参照入力とプラント出力yとを用いて相補感度関数tを上記式(11)を用いて求められ、制御システムSの入力信号である目標値dを相補感度関数tに入力される。相補感度関数tの出力である出力y*と、参照モデルの出力M(z)dとの誤差を最小にする制御器Cのパラメータθが、最適化手法により求められる。
[3.シミュレーションによる検証]
[3.1.システムの定式化]
本願の発明者は、安定性を考慮したFRITについてシミュレーションを用いて検証した。そこで、まず、シミュレーションで用いる制御対象P、参照モデルM、制御器Cの定式化を行う。
[3.1.システムの定式化]
本願の発明者は、安定性を考慮したFRITについてシミュレーションを用いて検証した。そこで、まず、シミュレーションで用いる制御対象P、参照モデルM、制御器Cの定式化を行う。
制御対象Pは次式(17)に示すように、ばね質量系とする。ここで、m、c、k、及びLは、それぞれ質量、粘性係数、ばね剛性、及びむだ時間を表す。
参照モデルMは、次式(18)に示すように1次遅れ系とする。ここで、τMはシステムの応答に関するパラメータであり、設計者が設定する。
制御器Cは次式(19)に示すPID(Proportional-Integral-Differential)制御器を用いる。ここで、Kp、Ki、及びKdはそれぞれ、Pゲイン、Iゲイン、及びDゲインである。
[3.2.シミュレーション]
シミュレーションにおいてサンプリング周期は4ミリ秒とし、制御器Cの離散化はゼロ次ホールドを用いた。なお、最適化計算には粒子群最適化を用いた。
シミュレーションにおいてサンプリング周期は4ミリ秒とし、制御器Cの離散化はゼロ次ホールドを用いた。なお、最適化計算には粒子群最適化を用いた。
図3(a)−(b)は、従来のFRITと実施の形態に係る手法とのそれぞれにより調整したパラメータθによる制御システムSの入出力データの時間履歴を示す図である。具体的には、図3(a)において、破線は参照応答、実線は実施の形態に係る手法における制御システムSの出力、一点鎖線は従来のFRITにおける制御システムSの出力を示している。図3(a)より、標準FRITでは制御システムSの出力が発散してグラフの枠内に収まっていないのに対し、実施の形態に係る手法は、制御の構造上モデルマッチングが実現できない場合においても閉ループ系である制御システムSが安定となる制御器Cのパラメータθが得られていることがわかる。
図3(b)は、図3(a)に対応する入力信号を示す図である。図3(a)と同様に、実線は実施の形態に係る手法における入力、一点鎖線は従来のFRITにおける入力を示している。図3(b)より、標準FRITでは入力も不安定となっている。一方、実施の形態に係る手法では、入力も安定している。
[4.小括]
従来のFRITでは、制御器Cのパラメータθの最適解を用いて制御しても、閉ループ系の挙動が不安定になるという問題があった。実施の形態に係る手法では、時間領域で擬似参照信号r(θ,k)とプラント出力yとの伝達関数である相補感度関数Tを同定することで、誤差を正しく検知することを可能にした。実施の形態に係る手法により、閉ループ系が安定である制御器Cのパラメータθを求められる。
従来のFRITでは、制御器Cのパラメータθの最適解を用いて制御しても、閉ループ系の挙動が不安定になるという問題があった。実施の形態に係る手法では、時間領域で擬似参照信号r(θ,k)とプラント出力yとの伝達関数である相補感度関数Tを同定することで、誤差を正しく検知することを可能にした。実施の形態に係る手法により、閉ループ系が安定である制御器Cのパラメータθを求められる。
[5.参考文献]
(FRIT)
・相馬 将太郎, 金子 修, 藤井 隆雄, 一回の実験データに基づく制御器パラメータチューニングの新しいアプローチ Fictitious Reference Iterative Tuning の提案, システム制御情報学会論文誌, Vol. 17, No.12 (2004), pp. 528-536
・奥谷 明大, 金子 修, 山本 茂, FRITを用いた多入出力むだ時間系に対するスミス補償器のチューニング, システム制御情報学会論文誌, Vol. 28, No 2 (2015), pp. 58-65
・データを直接用いた制御器のパラメータチューニング, 金子修, 計測と制御, Vol.43, No.11 (2008), pp903-908
(FRIT)
・相馬 将太郎, 金子 修, 藤井 隆雄, 一回の実験データに基づく制御器パラメータチューニングの新しいアプローチ Fictitious Reference Iterative Tuning の提案, システム制御情報学会論文誌, Vol. 17, No.12 (2004), pp. 528-536
・奥谷 明大, 金子 修, 山本 茂, FRITを用いた多入出力むだ時間系に対するスミス補償器のチューニング, システム制御情報学会論文誌, 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)
・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の提案, 電気学会論文誌D(産業応用部門誌), 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.
・弓場井 一裕, 藤井 宏樹, 平井 淳之, パラメータ更新時の閉ループシステムの安定性を考慮したFCbTの提案, 電気学会論文誌D(産業応用部門誌), 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.
・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.
[6.実施の形態]
以上を踏まえ、本発明の実施の形態について説明する。
以上を踏まえ、本発明の実施の形態について説明する。
図4は、実施の形態に係る情報処理装置1の機能構成を模式的に示す図である。情報処理装置1は、記憶部2と制御部3とを備える。図4において、矢印は主なデータの流れを示しており、図4に示していないデータの流れがあってもよい。図4において、各機能ブロックはハードウェア(装置)単位の構成ではなく、機能単位の構成を示している。そのため、図4に示す機能ブロックは単一の装置内に実装されてもよく、あるいは複数の装置内に分かれて実装されてもよい。機能ブロック間のデータの授受は、データバス、ネットワーク、可搬記憶媒体等、任意の手段を介して行われてもよい。
記憶部2は、情報処理装置1を実現するコンピュータのBIOS(Basic Input Output System)等を格納するROM(Read Only Memory)や情報処理装置1の作業領域となるRAM(Random Access Memory)、OS(Operating System)やアプリケーションプログラム、当該アプリケーションプログラムの実行時に参照される種々の情報を格納するHDD(Hard Disk Drive)やSSD(Solid State Drive)等の大容量記憶装置である。
制御部3は、情報処理装置1のCPU(Central Processing Unit)やGPU(Graphics Processing Unit)等のプロセッサであり、記憶部2に記憶されたプログラムを実行することによって時系列データ取得部30、入力信号推定部31、相補感度関数算出部32、プラント出力算出部33、モデル出力取得部34、及びパラメータ更新部35として機能する。
なお、図4は、情報処理装置1が単一の装置で構成されている場合の例を示している。しかしながら、情報処理装置1は、例えばクラウドコンピューティングシステムのように複数のプロセッサやメモリ等の計算リソースによって実現されてもよい。この場合、制御部3を構成する各部は、複数の異なるプロセッサの中の少なくともいずれかのプロセッサがプログラムを実行することによって実現される。
情報処理装置1は、制御器Cと制御器Cの出力を入力とする制御対象Pとを備える制御システムSにおいて、制御器Cのパラメータθを算出するための装置である。図4に示すように、制御システムSは、制御対象Pの出力が制御器Cの入力にフィードバックされる閉ループ系である。
時系列データ取得部30は、制御器Cの出力の時系列データである第1時系列データと、制御対象Pの出力の時系列データである第2時系列データとを取得する。ここで、第1時系列データは上述した制御対象Pに入力すべき制御量uに対応し、第2時系列データは上述した制御対象Pの出力yに対応する。したがって、以下本明細書において、「第1時系列データu」、「第2時系列データy」と記載することがある。
入力信号推定部31は、制御器Cのパラメータθ、第1時系列データu、及び第2時系列データyから、制御器Cに入力する入力信号の時系列データの推定値である第3時系列データを推定する。第3時系列データは、上述した擬似参照信号r(θ,k)に対応する。したがって、以下、「第3時系列データr(θ)」と記載することがある。
相補感度関数算出部32は、第2時系列データuと第3時系列データr(θ)とに基づいて、上述した式(11)を用いて制御器Cに対する相補感度関数tを算出する。プラント出力算出部33は、上述した式(12)又は式(13)を用いて、制御システムSの入力信号dを相補感度関数tに印加したときの出力である第4時系列データを算出する。第4時系列データは、制御システムSの入力信号d(目標値)を相補感度関数Tに印加した場合の出力y*に対応する。したがって、以下、「第4時系列データy*」と記載することがある。
制御対象Pに入力すべき制御量である第1時系列データu、制御対象Pの出力である第2時系列データy、及び制御システムSの入力信号である第4時系列データdは、いずれも観測により取得可能な量である。情報処理装置1は、観測により取得した第1時系列データu、第2時系列データy、及び入力信号dを用いることで、制御システムSに入力信号dを入力したときの出力y*を算出することができる。出力y*は、プラント(制御対象P)からの出力である。情報処理装置1は、この出力y*の振る舞い(例えば、発散するか否か、振動するか否か、収束するか否か等)を解析することで、FRITにおいて制御器Cのパラメータθを用いたプラントの制御の安定性を評価することができる。
図2及び図4に示すように、制御システムSは制御対象Pの入力信号dに対する出力を所定の出力とすることを実現する参照モデルMを備えている。なお、参照モデルMは、制御対象Pの出力が設計者が希望する出力となるように、設計者が決定する。モデル出力取得部34は、参照モデルMに入力信号dを入力した場合の出力の時系列データである第5時系列データを取得する。第5時系列データは、上述した参照モデルの出力M(z)dに対応する。以下、第5時系列データを「第5時系列データyd」と記載することがある。
ここで、情報処理装置1がFRITにおける制御器Cのパラメータθを用いた制御の安定性を評価することができるのであれば、情報処理装置1は、制御が安定化するようにパラメータθを最適化することもできる。これを実現するために、パラメータ更新部35は、第4時系列データy*と第5時系列データydとの誤差に関する評価関数J*(θ)の評価値に基づいて、制御器Cのパラメータθを更新する。
より具体的には、上述した式(15)及び式(16)に示すように、パラメータ更新部35が用いる評価関数J*(θ)は、第4時系列データy*と第5時系列データydとの誤差の二乗和である。パラメータ更新部35は、評価関数J*(θ)の評価値が小さくなるように、反復処理によって制御器Cのパラメータθを更新する。すなわち、パラメータ更新部35は、第4時系列データy*と第5時系列データydとの誤差の二乗和が小さくなるという意味において最適なパラメータθを反復によって求める。
一般に、第4時系列データy*が発散したり振動したりすると、第4時系列データy*と第5時系列データydとの誤差の二乗和は大きくなる。パラメータ更新部35が第4時系列データy*と第5時系列データydとの誤差の二乗和が小さくなるように制御器Cのパラメータθを更新することにより、情報処理装置1は、FRITにおいて制御器Cのパラメータθを用いたプラントの制御を安定化させることができる。
なお、パラメータ更新部35は、評価関数の評価値が小さくなるように制御器Cのパラメータθを更新できるのであればどのような最適化手法を用いてもよい。一例としてパラメータ更新部35は、あらかじめ定められた所定の回数を反復回数の上限として、粒子群最適化の手法を用いて制御器Cのパラメータθを更新してもよい。
<情報処理装置1が実行する情報処理方法の処理フロー>
図5は、実施の形態に係る情報処理装置1が実行する情報処理の流れを説明するためのフローチャートである。本フローチャートにおける処理は、例えば情報処理装置1が起動したときに開始する。
図5は、実施の形態に係る情報処理装置1が実行する情報処理の流れを説明するためのフローチャートである。本フローチャートにおける処理は、例えば情報処理装置1が起動したときに開始する。
時系列データ取得部30は、制御器Cの出力の時系列データである第1時系列データuを取得する(S2)。また、時系列データ取得部30は、制御対象Pの出力の時系列データである第2時系列データyを取得する(S4)。
入力信号推定部31は、制御器Cのパラメータθ、第1時系列データu、及び第2時系列データyから、式(1)を用いて、制御器Cに入力する入力信号dの時系列データの推定値である第3時系列データr(θ)を推定する(S6)。
相補感度関数算出部32は、第2時系列データyと第3時系列データr(θ)とに基づいて、式(11)を用いて制御器Cに対する相補感度関数t(すなわち、擬似参照信号r(θ,k)に対する相補感度関数t)を算出する(S8)。プラント出力算出部33は、式(14)を用いて、制御システムSの入力信号dを相補感度関数tに印加したときの出力である第4時系列データy*(すなわち、入力信号dに対する制御対象Pの出力)を算出する(S10)。
モデル出力取得部34は、参照モデルMに入力信号dを入力した場合の出力の時系列データである第5時系列データydを取得する(S12)。パラメータ更新部35は、式(15)及び式(16)を用いて、第4時系列データy*と第5時系列データydとの誤差に関する評価関数J(θ)の評価値を算出する(S14)。パラメータ更新部35は、評価関数J(θ)の評価値が小さくなるように、反復処理によって制御器Cのパラメータθを更新する(S16)。
パラメータ更新部35が制御器Cのパラメータθを更新すると、本フローチャートにおける処理は終了する。情報処理装置1は、上記の処理をオンラインで繰り返すことにより、制御器Cのパラメータθの更新を継続する。
<実施の形態に係る情報処理装置1が奏する効果>
以上説明したように、実施の形態に係る情報処理装置1によれば、FRITにおいて制御器Cのパラメータθを用いたプラントの制御の安定性を評価することができる。また、情報処理装置1は、閉ループ系の制御が安定化するようにFRITにおける制御器Cパラメータθを最適化することもできる。
以上説明したように、実施の形態に係る情報処理装置1によれば、FRITにおいて制御器Cのパラメータθを用いたプラントの制御の安定性を評価することができる。また、情報処理装置1は、閉ループ系の制御が安定化するようにFRITにおける制御器Cパラメータθを最適化することもできる。
以上、本発明を実施の形態を用いて説明したが、本発明の技術的範囲は上記実施の形態に記載の範囲には限定されず、その要旨の範囲内で種々の変形及び変更が可能である。例えば、装置の全部又は一部は、任意の単位で機能的又は物理的に分散・統合して構成することができる。また、複数の実施の形態の任意の組み合わせによって生じる新たな実施の形態も、本発明の実施の形態に含まれる。組み合わせによって生じる新たな実施の形態の効果は、もとの実施の形態の効果を併せ持つ。
1・・・情報処理装置
2・・・記憶部
3・・・制御部
30・・・時系列データ取得部
31・・・入力信号推定部
32・・・相補感度関数算出部
33・・・プラント出力算出部
34・・・モデル出力取得部
35・・・パラメータ更新部
C・・・制御器
M・・・参照モデル
P・・・制御対象
S・・・制御システム
2・・・記憶部
3・・・制御部
30・・・時系列データ取得部
31・・・入力信号推定部
32・・・相補感度関数算出部
33・・・プラント出力算出部
34・・・モデル出力取得部
35・・・パラメータ更新部
C・・・制御器
M・・・参照モデル
P・・・制御対象
S・・・制御システム
Claims (4)
- 制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出する情報処理装置であって、
前記制御器の出力の時系列データである第1時系列データと、前記制御対象の出力の時系列データである第2時系列データとを取得する時系列データ取得部と、
前記制御器のパラメータ、前記第1時系列データ、及び前記第2時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第3時系列データを推定する入力信号推定部と、
前記第2時系列データと前記第3時系列データとに基づいて、前記制御器に対する相補感度関数を算出する相補感度関数算出部と、
前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第4時系列データを算出するプラント出力算出部と、
を備える情報処理装置。 - 前記制御システムは、前記制御対象の入力に対する所定の応答を実現するモデルである参照モデルをさらに備えており、
前記情報処理装置は、
前記参照モデルに前記入力信号を入力した場合の出力の時系列データである第5時系列データを取得するモデル出力取得部と、
前記第4時系列データと前記第5時系列データとの誤差に関する評価関数の評価値に基づいて前記制御器のパラメータを更新するパラメータ更新部と、
をさらに備える請求項1に記載の情報処理装置。 - 前記評価関数は、前記第4時系列データと前記第5時系列データとの誤差の二乗和であり、
前記パラメータ更新部は、前記評価関数の評価値が小さくなるように、反復処理によって前記制御器のパラメータを更新する、
請求項2に記載の情報処理装置。 - 制御器と、前記制御器の出力を入力とする制御対象と、を備え、前記制御対象の出力が前記制御器の入力にフィードバックされる制御システムにおいて、前記制御器のパラメータを算出するコンピュータに、
前記制御器の出力の時系列データである第1時系列データを取得する機能と、
前記制御対象の出力の時系列データである第2時系列データを取得する機能と、
前記制御器のパラメータ、前記第1時系列データ、及び前記第2時系列データから、前記制御器に入力する入力信号の時系列データの推定値である第3時系列データを推定する機能と、
前記第2時系列データと前記第3時系列データとに基づいて、前記制御器に対する相補感度関数を算出する機能と、
前記制御システムの入力信号を前記相補感度関数に印加したときの出力である第4時系列データを算出する機能と、
を実現させるプログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019106908A JP2020201626A (ja) | 2019-06-07 | 2019-06-07 | 情報処理装置及びプログラム |
PCT/JP2020/022049 WO2020246531A1 (ja) | 2019-06-07 | 2020-06-04 | 情報処理装置、プログラム、及び算出方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019106908A JP2020201626A (ja) | 2019-06-07 | 2019-06-07 | 情報処理装置及びプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2020201626A true JP2020201626A (ja) | 2020-12-17 |
Family
ID=73652221
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019106908A Pending JP2020201626A (ja) | 2019-06-07 | 2019-06-07 | 情報処理装置及びプログラム |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP2020201626A (ja) |
WO (1) | WO2020246531A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116577995B (zh) * | 2023-07-05 | 2023-11-21 | 广州惠特节能科技有限公司 | 热泵烘干机组参数智能控制方法及系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5362339B2 (ja) * | 2008-12-03 | 2013-12-11 | 住友重機械工業株式会社 | モータ制御装置 |
WO2017088054A1 (en) * | 2015-11-23 | 2017-06-01 | The University Of British Columbia | Methods and systems for closed-loop control of drug administration |
JP6812653B2 (ja) * | 2016-03-31 | 2021-01-13 | ブラザー工業株式会社 | パラメータ更新方法、パラメータ更新装置、及びプログラム |
JP6870330B2 (ja) * | 2017-01-11 | 2021-05-12 | 横河電機株式会社 | データ処理装置及びデータ処理方法 |
-
2019
- 2019-06-07 JP JP2019106908A patent/JP2020201626A/ja active Pending
-
2020
- 2020-06-04 WO PCT/JP2020/022049 patent/WO2020246531A1/ja active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2020246531A1 (ja) | 2020-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6143943B2 (ja) | 二つのサブシステムの協調シミュレーション方法及び装置 | |
JP2021051462A (ja) | 情報処理装置及びプログラム | |
JP5723582B2 (ja) | スライディング・モード制御方式において、制御器を構成する方法、プログラム及びシステム | |
JP2013525910A (ja) | 制御装置及び制御のための出力変数を計算する方法 | |
JP6718500B2 (ja) | 生産システムにおける出力効率の最適化 | |
KR102038703B1 (ko) | 동적 전이 앙상블 모형을 통한 실시간 다변량 시계열 예측방법 및 그 시스템 | |
KR102577188B1 (ko) | 목표 시스템에 대한 제어 시스템 생성 | |
JP7497516B2 (ja) | 等式制約を代数モデルに課すための射影方法 | |
Miller et al. | 23 Real-Time Simulation of Physical Systems Using Simscape™ | |
Chen et al. | Online identification of time‐delay jump Markov autoregressive exogenous systems with recursive expectation‐maximization algorithm | |
WO2020246531A1 (ja) | 情報処理装置、プログラム、及び算出方法 | |
JP7402482B2 (ja) | 前段階コシミュレーション方法、デバイス、コンピュータ可読媒体、及びプログラム | |
EP3582153A1 (en) | Generating hybrid models of physical systems | |
CN117290965B (zh) | 车辆仿真软件的车辆模型仿真试验仿真方法、设备及介质 | |
WO2021006332A1 (ja) | 情報処理装置、プログラム、及び算出方法 | |
JP2010049392A (ja) | Pidコントローラのチューニング装置、pidコントローラのチューニング用プログラムおよびpidコントローラのチューニング方法 | |
JP7115656B1 (ja) | 制御装置、制御方法及びプログラム | |
JP6919856B2 (ja) | 強化学習プログラム、強化学習方法、および強化学習装置 | |
JP2021012605A (ja) | 伝達関数の予測方法 | |
JP6984597B2 (ja) | 線形パラメータ変動モデル推定システム、方法およびプログラム | |
Singh et al. | Randomized probabilistic approach for parametric uncertainties in unmanned helicopters | |
CN108920842B (zh) | 一种潜艇动力学模型参数在线估计方法及装置 | |
JP7207473B1 (ja) | 情報処理装置 | |
JP7207474B1 (ja) | 情報処理装置 | |
Gogilan et al. | Detection and characterization of local nonlinearity in a clamped‐clamped beam with state space realization |