(57)【要約】【課題】制御器の挙動を容易に知ることができる技術を提供する。【解決手段】本開示の一側面によれば、制御器に目標指令ｒを入力したときの制御入力ｕ及び制御出力ｙを予測するための方法が提供される。本方法では、制御入力ｕｍと制御出力ｙｍとの組合せの時系列データが取得される（Ｓ１１０）。この時系列データから、目標応答ｙｒと制御出力ｙｍとの間の誤差を評価するための評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、設計変数θを固定値θ０に設定した状態で、評価関数の値を最小にする設計変数ρの値ρ＊が算出される（Ｓ１２０）。評価関数は、制御器の伝達関数Ｃ（ρ）と、目標指令ｒに対する目標応答ｙｒを規定する目標応答伝達関数Ｔｄ（θ）と、を要素に含む。更に、設計変数ρを値ρ＊に設定した状態で、評価関数の値を最小にする設計変数θの値θ＊が算出される（Ｓ１３０）。制御入力ｕ及び制御出力ｙの予測は、値ρ＊及び値θ＊に基づく。【選択図】図３
2019175319.tif 000002
【特許請求の範囲】
【請求項１】
目標指令ｒ及び制御出力ｙに基づき制御対象への制御入力ｕを決定する制御器に前記目標指令ｒを入力したときの前記制御入力ｕ及び前記制御出力ｙを予測するための予測方法であって、
前記制御入力ｕ及び前記制御出力ｙのサンプルとして、制御入力ｕｍと制御出力ｙｍとの組合せの時系列データを取得することと、
設計変数ρを有する前記制御器の伝達関数Ｃ（ρ）と、前記目標指令ｒに対する目標応答ｙｒを規定する設計変数θを有する目標応答伝達関数Ｔｄ（θ）と、を要素に含み、前記時系列データから、前記目標応答ｙｒと前記制御出力ｙｍとの間の誤差を評価するための評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、前記設計変数θを固定値θ０に設定した状態で、前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする前記設計変数ρの値ρ＊を算出することと、
前記設計変数ρを前記値ρ＊に設定した状態で、前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする前記設計変数θの値θ＊を算出することと、
前記設計変数ρを前記値ρ＊に設定した伝達関数Ｃ（ρ＊）及び前記設計変数θを前記値θ＊に設定した目標応答伝達関数Ｔｄ（θ＊）に基づいて、前記目標指令ｒに対する前記制御入力ｕの予測値ｕｐ及び前記制御出力ｙの予測値ｙｐを算出することと、
を含む予測方法。
【請求項２】
前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）が、ＦＲＩＴ（Ｆｉｃｔｉｔｉｏｕｓ  Ｒｅｆｅｒｅｎｃｅ  Ｉｔｅｒａｔｉｖｅ  Ｔｕｎｉｎｇ）技術に基づく関数である請求項１記載の予測方法。
【請求項３】
前記制御器が前記目標指令ｒと前記制御出力ｙとの偏差（ｒ−ｙ）に基づき前記制御入力ｕ＝Ｃ（ρ）（ｒ−ｙ）を決定するように構成され、
前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）が、式
【数１】
2019175319.tif 000015
に従う関数ｇ（ρ，θ，ｕｍ，ｙｍ）のノルムに対応する請求項１又は請求項２記載の予測方法。
【請求項４】
前記予測値ｕｐ及び前記予測値ｙｐが、式
【数２】
2019175319.tif 000016
に従って算出される請求項３記載の予測方法。
【請求項５】
前記予測値ｕｐ及び前記予測値ｙｐが、外乱を加味して算出される請求項１〜請求項３のいずれか一項記載の予測方法。
【請求項６】
前記予測値ｕｐ及び前記予測値ｙｐが、前記制御出力ｙに加算される外乱ｗを加味した式
【数３】
2019175319.tif 000017
に従って算出される請求項３記載の予測方法。
【請求項７】
前記予測値ｕｐ及び前記予測値ｙｐが、前記制御入力ｕに加算される外乱ｆを加味した式
【数４】
2019175319.tif 000018
に従って算出される請求項３記載の予測方法。
【請求項８】
前記予測値ｕｐ及び前記予測値ｙｐが、前記制御出力ｙに乗算される外乱Ｇを加味した式
【数５】
2019175319.tif 000019
に従って算出される請求項３記載の予測方法。
【請求項９】
前記予測値ｕｐ及び前記予測値ｙｐが、前記制御対象の変化を加味して算出される請求項１〜請求項３のいずれか一項記載の予測方法。
【請求項１０】
前記予測値ｕｐ及び前記予測値ｙｐが、前記制御対象の変化δＱを加味した式
【数６】
2019175319.tif 000020
に従って算出される請求項３記載の予測方法。
【請求項１１】
目標指令ｒ及び制御出力ｙに基づき制御対象への制御入力ｕを決定する制御器に前記目標指令ｒを入力したときの前記制御入力ｕ及び前記制御出力ｙを予測する処理をコンピュータに実行させるためのコンピュータプログラムであって、前記コンピュータに、
前記制御入力ｕ及び前記制御出力ｙのサンプルとして、制御入力ｕｍと制御出力ｙｍとの組合せの時系列データを取得することと、
設計変数ρを有する前記制御器の伝達関数Ｃ（ρ）と、前記目標指令ｒに対する目標応答ｙｒを規定する設計変数θを有する目標応答伝達関数Ｔｄ（θ）と、を要素に含み、前記時系列データから、前記目標応答ｙｒと前記制御出力ｙｍとの間の誤差を評価するための評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、前記設計変数θを固定値θ０に設定した状態で、前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする前記設計変数ρの値ρ＊を算出することと、
前記設計変数ρを前記値ρ＊に設定した状態で、前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする前記設計変数θの値θ＊を算出することと、
前記設計変数ρを前記値ρ＊に設定した伝達関数Ｃ（ρ＊）及び前記設計変数θを前記値θ＊に設定した目標応答伝達関数Ｔｄ（θ＊）に基づいて、前記目標指令ｒに対する前記制御入力ｕの予測値ｕｐ及び前記制御出力ｙの予測値ｙｐを算出することと、
を実行させるためのコンピュータプログラム。
【請求項１２】
目標指令ｒ及び制御出力ｙに基づき制御対象への制御入力ｕを決定する制御器に前記目標指令ｒを入力したときの前記制御入力ｕ及び前記制御出力ｙを予測するための予測装置であって、
前記制御入力ｕ及び前記制御出力ｙのサンプルとして、制御入力ｕｍと制御出力ｙｍとの組合せの時系列データを取得するように構成される取得ユニットと、
設計変数ρを有する前記制御器の伝達関数Ｃ（ρ）と、前記目標指令ｒに対する目標応答ｙｒを規定する設計変数θを有する目標応答伝達関数Ｔｄ（θ）と、を要素に含み、前記時系列データから、前記目標応答ｙｒと前記制御出力ｙｍとの間の誤差を評価するための評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、前記設計変数θを固定値θ０に設定した状態で、前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする前記設計変数ρの値ρ＊を算出するように構成される第一算出ユニットと、
前記設計変数ρを前記値ρ＊に設定した状態で、前記評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする前記設計変数θの値θ＊を算出するように構成される第二算出ユニットと、
前記設計変数ρを前記値ρ＊に設定した伝達関数Ｃ（ρ＊）及び前記設計変数θを前記値θ＊に設定した目標応答伝達関数Ｔｄ（θ＊）に基づいて、前記目標指令ｒに対する前記制御入力ｕの予測値ｕｐ及び前記制御出力ｙの予測値ｙｐを算出するように構成される予測ユニットと、
を含む予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、制御入力及び制御出力の予測方法、予測装置、及びコンピュータプログラムに関する。
【背景技術】
【０００２】
制御器のパラメータを調整するための技術として、ＦＲＩＴ（Ｆｉｃｔｉｔｉｏｕｓ  Ｒｅｆｅｒｅｎｃｅ  Ｉｔｅｒａｔｉｖｅ  Ｔｕｎｉｎｇ）技術が既に知られている。この技術では、制御出力が目標応答に近づくような制御器のパラメータが、評価関数に基づき算出される。この算出値が制御器に設定されることにより、制御器は、目標応答に対応した制御を実現する。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１７−６８６５８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記評価関数に基づき算出される制御器のパラメータは、制御出力を目標応答に近づけるものである。しかしながら、上記技術に従って単に評価関数に基づく算出値を制御器のパラメータに設定しても、制御出力と目標応答との間には誤差が依然として残る。即ち、従来技術に従って制御器のパラメータを算出するだけでは、制御器の挙動を知ることはできない。この場合には制御器の挙動を知るために、制御器を組み込んだ制御系を実際に作製して動作させる必要があり、設計者の作業負荷が大きい。
【０００５】
そこで、本開示の一側面によれば、制御器の挙動を容易に知ることができる技術を提供できることが望ましい。
【課題を解決するための手段】
【０００６】
本開示の一側面によれば、制御器に目標指令ｒを入力したときの制御入力ｕ及び制御出力ｙを予測するための予測方法が提供される。制御器は、目標指令ｒ及び制御出力ｙに基づき制御対象への制御入力ｕを決定するように構成される。
【０００７】
この予測方法は、制御入力ｕ及び制御出力ｙのサンプルとして、制御入力ｕｍと制御出力ｙｍとの組合せの時系列データを取得することを含む。予測方法は更に、時系列データから、目標応答ｙｒと制御出力ｙｍとの間の誤差を評価するための評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、設計変数θを固定値θ０に設定した状態で評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする設計変数ρの値ρ＊を算出することを含む。
【０００８】
評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）は、設計変数ρを有する制御器の伝達関数Ｃ（ρ）と、目標指令ｒに対する目標応答ｙｒを規定する目標応答伝達関数Ｔｄ（θ）と、を要素に含む。目標応答伝達関数Ｔｄ（θ）は、設計変数θを有する。
【０００９】
予測方法は更に、設計変数ρを値ρ＊に設定した状態で、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする設計変数θの値θ＊を算出することを含む。予測方法は更に、設計変数ρを値ρ＊に設定した伝達関数Ｃ（ρ＊）及び設計変数θを値θ＊に設定した目標応答伝達関数Ｔｄ（θ＊）に基づいて、目標指令ｒに対する制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出することを含む。
【００１０】
この予測方法によれば、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、設計変数ρ及び設計変数θの適値ρ＊及びθ＊が求められる。従って、伝達関数Ｃ（ρ＊）及び目標応答伝達関数Ｔｄ（θ＊）に基づけば、目標指令ｒに対する精度のよい制御入力ｕ及び制御出力ｙの予測値ｕｐ及び予測値ｙｐを算出することができる。従って、本開示によれば、目標指令ｒに対する制御器の挙動を、実機を用いた繰返しの実験なしで、容易に知ることができる。
【００１１】
本開示の一側面によれば、目標指令ｒ及び制御出力ｙに基づき制御対象への制御入力ｕを決定する制御器に目標指令ｒを入力したときの制御入力ｕ及び制御出力ｙを予測する処理をコンピュータに実行させるためのコンピュータプログラムが提供されてもよい。
【００１２】
コンピュータプログラムは、制御入力ｕ及び制御出力ｙのサンプルとして、制御入力ｕｍと制御出力ｙｍとの組合せの時系列データを取得することと、設計変数ρを有する制御器の伝達関数Ｃ（ρ）と、目標指令ｒに対する目標応答ｙｒを規定する設計変数θを有する目標応答伝達関数Ｔｄ（θ）と、を要素に含み、時系列データから、目標応答ｙｒと制御出力ｙｍとの間の誤差を評価するための評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、設計変数θを固定値θ０に設定した状態で、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする設計変数ρの値ρ＊を算出することと、設計変数ρを値ρ＊に設定した状態で、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする設計変数θの値θ＊を算出することと、設計変数ρを値ρ＊に設定した伝達関数Ｃ（ρ＊）及び設計変数θを値θ＊に設定した目標応答伝達関数Ｔｄ（θ＊）に基づいて、目標指令ｒに対する制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出することと、をコンピュータに実行させるように構成されてもよい。本開示の一側面では、上述した予測方法をコンピュータに実行させるためのコンピュータプログラムが提供されてもよい。本開示の一側面によれば、コンピュータプログラムを記録したコンピュータ読取可能な一時的でない記録媒体が提供されてもよい。
【００１３】
本開示の一側面によれば、目標指令ｒ及び制御出力ｙに基づき制御対象への制御入力ｕを決定する制御器に目標指令ｒを入力したときの制御入力ｕ及び制御出力ｙを予測するための予測装置が提供されてもよい。
【００１４】
予測装置は、取得ユニットと、第一算出ユニットと、第二算出ユニットと、予測ユニットとを備えていてもよい。取得ユニットは、制御入力ｕ及び制御出力ｙのサンプルとして、制御入力ｕｍと制御出力ｙｍとの組合せの時系列データを取得するように構成されてもよい。第一算出ユニットは、制御器の伝達関数Ｃ（ρ）と、目標指令ｒに対する目標応答ｙｒを規定する目標応答伝達関数Ｔｄ（θ）と、を要素に含み、時系列データから、目標応答ｙｒと制御出力ｙｍとの間の誤差を評価するための評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）に従って、設計変数θを固定値θ０に設定した状態で、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする設計変数ρの値ρ＊を算出するように構成されてもよい。
【００１５】
第二算出ユニットは、設計変数ρを値ρ＊に設定した状態で、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）の値を最小にする設計変数θの値θ＊を算出するように構成されてもよい。予測ユニットは、設計変数ρを値ρ＊に設定した伝達関数Ｃ（ρ＊）及び設計変数θを値θ＊に設定した目標応答伝達関数Ｔｄ（θ＊）に基づいて、目標指令ｒに対する制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出するように構成されてもよい。
【図面の簡単な説明】
【００１６】
【図１】予測装置の構成を表すブロック図である。
【図２】フィードバック制御系の構成を表すブロック図である。
【図３】プロセッサが実行する予測前処理を表すフローチャートである。
【図４】プロセッサが実行する予測処理を表すフローチャートである。
【図５】図５Ａは、制御系の具体例を表すブロック図であり、図５Ｂは、外乱入力を含む制御系の具体例を表すブロック図である。
【図６】図６Ａ−６Ｂは、制御入力の予測値を表すグラフである。
【図７】図７Ａ−７Ｃは、制御出力の予測値を表すグラフである。
【図８】図８Ａは、制御出力に加算される外乱を含む制御系を表すブロック図であり、図８Ｂは、図８Ａに示す制御系と等価な制御系を表すブロック図である。
【図９】図９Ａは、正弦波外乱を含む制御系における制御入力の予測値を表すグラフであり、図９Ｂは、図９Ａに対応する制御出力の予測値を表すグラフである。
【図１０】図１０Ａは、制御出力に乗算される外乱を含む制御系を表すブロック図であり、図１０Ｂは、図１０Ａに示す制御系と等価な制御系を表すブロック図である。
【図１１】図１１Ａは、制御入力に加算される外乱を含む制御系を表すブロック図であり、図１１Ｂは、図１１Ａに示す制御系と等価な制御系を表すブロック図である。
【図１２】図１２Ａは、プラントの伝達特性が変化する制御系の第一例を表すブロック図であり、図１２Ｂは、図１２Ａに示す制御系と等価な制御系を表すブロック図であり、図１２Ｃは、プラントの伝達特性の具体例を表すブロック図である。
【図１３】図１３Ａは、プラントの伝達特性が変化する制御系の第二例を表すブロック図であり、図１３Ｂは、図１３Ａに示す制御系と等価な制御系を表すブロック図である。
【発明を実施するための形態】
【００１７】
以下に本開示の例示的実施形態を、図面と共に説明する。
本実施形態の予測装置１は、汎用の情報処理装置に、専用のコンピュータプログラムがインストールされて構成される。この予測装置１は、プロセッサ１１と、メモリ１２と、ストレージデバイス１３と、操作デバイス１５と、表示デバイス１７と、入出力インタフェース１９と、を備える。プロセッサ１１は、ＣＰＵを含み、メモリ１２は、ＲＡＭを含む。ストレージデバイス１３は、ハードディスクドライブ及びソリッドステートドライブの少なくとも一方を含む。ストレージデバイス１３は、コンピュータプログラム及びデータを備える。
【００１８】
プロセッサ１１は、ストレージデバイス１３に記録されたコンピュータプログラムに従う処理を実行するように構成される。メモリ１２は、プロセッサ１１による処理実行時に作業領域として使用される。
【００１９】
操作デバイス１５は、ユーザからの操作信号をプロセッサ１１に入力するように構成される。操作デバイス１５の例には、キーボード及びポインティングデバイスが含まれる。表示デバイス１７は、プロセッサ１１により制御されて、ユーザ向けの情報を表示するように構成される。表示デバイス１７の例には、液晶ディスプレイ及び有機ＥＬディスプレイが含まれる。
【００２０】
入出力インタフェース１９は、予測装置１の内部にデータを入力したり、予測装置１の外部にデータを出力したりするためのインタフェースである。入出力インタフェース１９の例には、ＵＳＢインタフェースが含まれる。
【００２１】
本実施形態によれば、この予測装置１を用いて、図２に示すフィードバック制御系１００における制御入力ｕ及び制御出力ｙが予測される。この予測は、制御器３０の挙動を、制御器３０の設計者が分析及び把握するために行われる。このフィードバック制御系１００は、減算器３１と制御器本体３５とを含む制御器３０と、制御対象のプラント５０とを備える。
【００２２】
減算器３１は、目標指令ｒと制御出力ｙとの偏差Ｅ＝ｒ−ｙを、制御器本体３５に入力するように構成される。制御器本体３５は、設定された伝達関数Ｃ（ρ）に従って、偏差Ｅに対応するプラント５０への制御入力ｕ＝Ｃ（ρ）Ｅを決定するように構成される。プラント５０からの制御出力ｙは、図示しない計測器により計測されて減算器３１に入力される。
【００２３】
設計者は、制御出力ｙが目標応答ｙｒに一致するように制御器３０の設計変数ρの値を決定する。しかしながら、従来技術では、設計変数ρの値を決定しても、その決定値に基づく制御器３０を実機で駆動させなければ、その挙動を正確に知ることができなかった。予測装置１は、そのような実機の駆動なしに、目標指令ｒに対する制御入力ｕ及び制御出力ｙを高精度に予測可能に構成される。
【００２４】
設計者は、予測装置１を利用して、制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを得て、制御器３０の挙動を分析及び把握することができる。設計者は、予測装置１を利用するために、制御入力ｕ及び制御出力ｙのサンプルデータを用意することができる。
【００２５】
具体的には、設計者は、設計変数ρを初期値ρ０に設定した制御器３０を用いて、フィードバック制御系１００を試験的に動作させると共にフィードバック制御系１００を観測し、サンプルデータとして、目標指令ｒに対する制御入力ｕの観測値ｕｍ及び制御出力ｙの観測値ｙｍの時系列データを作成することができる。初期値ρ０は、設計者が適宜に定めることができる。以下で述べる制御入力ｕｍは、制御入力ｕの観測値ｕｍを意味し、制御出力ｙｍは、制御出力ｙの観測値ｙｍを意味する。
【００２６】
設計者は、上記作成した時系列データを、入出力インタフェース１９を通じて予測装置１に入力することができる。予測装置１は、この時系列データを、ストレージデバイス１３に記憶することができる。更に、設計者は、予測装置１に、制御器３０の伝達関数Ｃ（ρ）及び、目標応答伝達関数Ｔｄ（θ）を、操作デバイス１５又は入出力インタフェース１９を通じて予測装置１に登録することができる。
【００２７】
目標応答伝達関数Ｔｄ（θ）は、目標指令ｒに対する制御出力ｙの規範値である目標応答ｙｒを規定する関数である。目標応答伝達関数Ｔｄ（θ）は、式ｙｒ＝Ｔｄ（θ）ｒに従う。登録される伝達関数Ｃ（ρ）及び目標応答伝達関数Ｔｄ（θ）は、それぞれ、可変パラメータρ及び可変パラメータθを含み、可変パラメータρは、制御器３０の設計変数ρに対応し、可変パラメータθは、目標応答ｙｒの設計変数θに対応する。設計者は更に、所望の目標応答ｙｒを規定する設計変数θの値θ０を、初期値θ０として予測装置１に登録することができる。
【００２８】
その後、設計者は、操作デバイス１５を通じて予測前処理の実行指令を予測装置１に入力することができる。予測前処理の実行指令が入力されると、予測装置１のプロセッサ１１は、図３に示す予測前処理の実行を開始し、予めストレージデバイス１３に登録された制御入力ｕｍ及び制御出力ｙｍの時系列データ、伝達関数Ｃ（ρ）、目標応答伝達関数Ｔｄ（θ）、及び、設計変数θの初期値θ０を読み出す（Ｓ１１０）。
【００２９】
読出後、プロセッサ１１は、設計変数θを初期値θ０に固定して、評価関数の値を最小にする設計変数ρの値ρ＊を算出する（Ｓ１２０）。ここで用いる評価関数は、ＦＲＩＴ技術に従う次の評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）である。
【００３０】
【数１】
2019175319.tif 000003
評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）は、次の関数ｇのノルム（ユークリッドノルム）の二乗に対応する。
【００３１】
【数２】
2019175319.tif 000004
上式から理解できるように、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）は、制御出力ｙと目標応答ｙｒとの誤差を評価可能に構成される。具体的に、評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）は、誤差が大きいほど大きな値を示す。
【００３２】
プロセッサ１１は、評価関数Ｊ（ρ，θ＝θ０，ｕｍ，ｙｍ）の値を最小にする値ρ＊を算出し終えると、次に設計変数ρを当該算出値ρ＊に固定して、評価関数Ｊ（ρ＝ρ＊，θ，ｕｍ，ｙｍ）の値を最小にする設計変数θの値θ＊を算出する（Ｓ１３０）。
【００３３】
その後、プロセッサ１１は、Ｓ１２０での算出値ρ＊及びＳ１３０での算出値θ＊を、予測式の定義値としてストレージデバイス１３に登録し（Ｓ１４０）、予測前処理を終了する。
【００３４】
制御器３０の設計者は、予測前処理の終了後、操作デバイス１５を通じて、図４に示す予測処理の実行指令を、予測装置１に入力することができる。予測処理の実行指令が入力されると、予測装置１のプロセッサ１１は、図４に示すように、予測式の定義情報を読み出す（Ｓ２１０）。具体的には、Ｓ１４０で登録された予測式の定義値ρ＊及びθ＊を読み出す。更には、予測式を定義する伝達関数Ｃ（ρ）及び目標応答伝達関数Ｔｄ（θ）を読み出す。
【００３５】
その後、プロセッサ１１は、読み出した定義情報に基づき、制御入力ｕの予測式及び制御出力ｙの予測式を組み立てる。組み立てられる制御入力ｕの予測式及び制御出力ｙの予測式、即ち、制御入力ｕの予測値ｕｐの算出式及び制御出力ｙの予測値ｙｐの算出式は、次の通りであり、Ｓ１４０で登録された値ρ＊及び値θ＊を含む。
【００３６】
【数３】
2019175319.tif 000005
【００３７】
プロセッサ１１は更に、予測に用いる目標指令ｒを取得する（Ｓ２３０）。プロセッサ１１は、この目標指令ｒを、操作デバイス１５を通じて操作者から取得してもよいし、ストレージデバイス１３から読み出すことにより取得してもよい。
【００３８】
その後、プロセッサ１１は、取得した目標指令ｒをＳ２２０で組み立てた予測式に代入して、制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出する（Ｓ２４０）。これにより、Ｓ２４０では、目標指令ｒに対応した予測値ｕｐ及び予測値ｙｐの時系列データが生成される。
【００３９】
プロセッサ１１は更に、予測値ｕｐ及び予測値ｙｐの情報を出力する（Ｓ２５０）。具体的に、プロセッサ１１は、Ｓ２４０で生成した予測値ｕｐ及び予測値ｙｐの時系列データをストレージデバイス１３に記録することができる。プロセッサ１１は、表示デバイス１７を制御して、予測値ｕｐ及び予測値ｙｐのグラフを表示デバイス１７に表示させてもよい。その後、プロセッサ１１は、予測処理を終了する。
【００４０】
制御器３０の設計者は、このようにして予測装置１から提供される予測値ｕｐ及び予測値ｙｐの時系列データ又は表示デバイス１７に表示されるグラフに基づいて、制御器３０の挙動を分析及び把握することができる。
【００４１】
ここで、図５Ａに示す具体的なフィードバック制御系１１０に基づく制御入力ｕ及び制御出力ｙの予測精度に関する実験結果を図６Ａ−６Ｂ及び図７Ａ−７Ｃに示す。このフィードバック制御系１１０は、減算器６１，６２と、ゲインＫｐを有する増幅器６３と、ゲインＫｖを有する増幅器６５と、微分器６７とを含む制御器６０と、プラント７０とを備える。実験は、このフィードバック制御系１１０に対し、次式に従う目標応答伝達関数Ｔｄ（θ）を設定して行われた。
【００４２】
【数４】
2019175319.tif 000006
【００４３】
このフィードバック制御系１１０に対する評価関数Ｊ（ρ，θ，ｕｍ，ｙｍ）は、上式（１）に、次のｒ（ρ）を設定して定義される。ここでＦは、微分作用素を表す。
【００４４】
【数５】
2019175319.tif 000007
【００４５】
図６Ａ−６Ｂ及び図７Ａ−７Ｃは、更に、設計変数ρの初期値ρ０を、Ｋｐ＝１０及びＫｖ＝０．１で定義し、設計変数θの初期値θ０を、ζ１＝ζ２＝１及びω１＝ω２＝２０で定義したときの実験結果を示す。
【００４６】
実験は、初期値ρ０の制御器６０を用いてフィードバック制御系１１０における制御入力ｕｍ及び制御出力ｙｍの時系列データを取得し、この時系列データ及び初期値θ０を用いて図３及び図４に示す処理を予測装置１に実行させ、それにより予測値ｕｐ及び予測値ｙｐの時系列データを取得することにより行われた。
【００４７】
この実験では、ρ＊に対応するＫｐ及びＫｖの値として、Ｋｐ＝１６．２及びＫｖ＝４．３９が得られ、θ＊に対応するζ１、ζ２、ω１、及びω２の値として、ζ１＝２．４７×１０５、ζ２＝８．０２×１０−１、ω１＝５．９４×１０７、及びω２＝８．７２が得られた。
【００４８】
図６Ａのグラフは、この値ρ＊及びθ＊を用いて算出した制御入力ｕの予測値ｕｐを破線で表し、予測前の観測値ｕｍを実線で表す、横軸が時間、縦軸が制御入力のグラフである。図６Ｂのグラフは、値ρ＊及びθ＊を用いた別の実験で得られた制御入力ｕの観測値ｕｂを実線で表すグラフである。このグラフには、図６Ａと同じ予測値ｕｐが示される。
【００４９】
図７Ａ−７Ｃのグラフは、図６Ａ−６Ｂと同じ値ρ＊及びθ＊を用いて算出した制御出力ｙの予測値ｙｐを破線で表す、横軸が時間、縦軸が制御出力のグラフである。図７Ａのグラフは更に、予測前の観測値ｙｍを実線で表し、図７Ｂのグラフは更に、初期値θ０に従う目標応答ｙｒ＝Ｔｄ（θ０）ｒを実線で表し、図７Ｃのグラフは更に、値ρ＊及びθ＊を用いた別の実験で得られた制御出力ｙの観測値ｙｂを実線で表す。
【００５０】
予測装置１から得られる予測値ｕｐ及び予測値ｙｐが高精度であることは、図６Ｂ及び図７Ｃのグラフから良く理解できる。制御出力ｙと目標応答ｙｒとの間の誤差は、ρ及びθが初期値ρ０及びθ０であるときと比較して、ρ及びθが調整後の値ρ＊及びθ＊であるとき、およそ１／２０倍まで改善した。
【００５１】
以上には、外乱やプラントの変化を加味せずに制御入力ｕ及び制御出力ｙを予測する例を説明したが、予測値ｕｐ及び予測値ｙｐは外乱やプラントの変化を加味して算出されてもよい。
【００５２】
以下で説明する変形例は、予測式が上記実施形態と異なる例である。一方、これらの変形例は、他の構成に関して上記実施形態と同じ構成を有する。従って、以下の変形例の説明では、上記実施形態と同一構成に関する説明を適宜省略する。
【００５３】
［第一変形例］
第一変形例の予測装置１は、制御出力ｙに加算される外乱ｗを加味して、制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出するように構成される。対応するフィードバック制御系１２０は、図８Ａに示され、このフィードバック制御系１２０と等価な制御系１２５が図８Ｂに示される。
【００５４】
この外乱ｗを加味した予測は、Ｓ２２０で次式に従う予測式を組み立てることにより実現することができる。プロセッサ１１は、外乱ｗの情報を、Ｓ２１０において予測式の定義情報の一部として取得する、又は、Ｓ２３０において目標指令ｒと共に取得することができる。
【００５５】
【数６】
2019175319.tif 000008
【００５６】
この予測式によれば、制御出力ｙに加算される外乱ｗを加味した制御入力ｕ及び制御出力ｙの予測が可能である。外乱ｗの例には、制御出力ｙを計測する計測器において制御出力ｙに加算されるノイズが含まれる。設計者は、予測式に与える外乱ｗを変更して、様々な予測値ｕｐ及び予測値ｙｐを得ることで、外乱ｗの影響を受ける制御器３０の挙動を詳しく分析及び把握することができる。
【００５７】
ここで図５Ｂに示す具体的なフィードバック制御系１３０に基づく制御入力ｕ及び制御出力ｙの予測精度に関する実験結果を、図９Ａ及び９Ｂに示す。このフィードバック制御系１３０は、図５Ａに示すフィードバック制御系１１０において、制御出力ｙに外乱ｗが加算される例に対応する。実験は、外乱ｗとして、正弦波外乱をフィードバック制御系１３０に入力したときの実験結果に対応する。
【００５８】
図９Ａのグラフは、制御入力ｕの予測値ｕｐを破線で表し、別の実験で得られた制御入力ｕの観測値ｕｂを実線で表す、横軸が時間、縦軸が制御入力のグラフである。図９Ｂのグラフは、制御出力ｙの予測値ｙｐを破線で表し、別の実験で得られた制御出力ｙの観測値ｙｂを実線で表す、横軸が時間、縦軸が制御出力のグラフである。図９Ｂには更に、初期値θ０の目標応答ｙｒ＝Ｔｄ（θ０）ｒが一点鎖線で表される。本変形例によれば、外乱ｗを含むフィードバック制御系１３０における制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを高精度に得られることが、図９Ａ及び図９Ｂから良く理解できる。
【００５９】
外乱ｗを加味した制御入力ｕ及び制御出力ｙの予測機能は、動作モードの一つとして、上記実施形態の予測装置１に組み込まれてもよい。即ち、予測装置１は、予測処理の実行指令を動作モードの指定と共に受けるように構成されてもよい。プロセッサ１１は、図４に示す予測処理として、指定された動作モードに対応した処理を実行するように構成されてもよい。
【００６０】
プロセッサ１１は、外乱を加味しない動作モードでは、最初に説明した実施形態に従う予測式、即ち上式（４）及び（５）に従う予測式に基づいた予測値ｕｐ及び予測値ｙｐを算出し、外乱ｗを加味した動作モードでは、本変形例に従う予測式、即ち上式（１２）及び（１３）に従う予測式に基づいた予測値ｕｐ及び予測値ｙｐを算出するように動作してもよい。指定された動作モードに従って予測式を切り替える予測装置１によれば、設計者は、制御器３０の挙動をより詳しく分析及び把握することができる。
【００６１】
［第二変形例］
第二変形例の予測装置１は、制御出力ｙに乗算される外乱Ｇを加味して、制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出するように構成される。対応するフィードバック制御系１４０は、図１０Ａに示され、このフィードバック制御系１４０と等価な制御系１４５が図１０Ｂに示される。
【００６２】
この外乱Ｇを加味した予測は、Ｓ２２０で次式に従う予測式を組み立てることにより実現することができる。プロセッサ１１は、外乱Ｇの情報を、Ｓ２１０において予測式の定義情報の一部として取得する、又は、Ｓ２３０において目標指令ｒと共に取得することができる。
【００６３】
【数７】
2019175319.tif 000009
【００６４】
この予測式によれば、制御出力ｙに乗算される外乱Ｇを加味した制御入力ｕ及び制御出力ｙの予測が可能である。外乱Ｇの例には、制御出力ｙを計測する計測器の特性により制御出力ｙに乗算される外乱が含まれる。本変形例に関する外乱Ｇを加味した制御入力ｕ及び制御出力ｙの予測機能は、第一変形例と同様、動作モードの一つとして、上記実施形態の予測装置１に組み込まれてもよい。
【００６５】
［第三変形例］
第三変形例の予測装置１は、制御入力ｕに加算される外乱ｆを加味して、制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出するように構成される。対応するフィードバック制御系１５０は、図１１Ａに示され、このフィードバック制御系１５０と等価な制御系１５５が図１１Ｂに示される。
【００６６】
この外乱ｆを加味した予測は、Ｓ２２０で次式に従う予測式を組み立てることにより実現することができる。プロセッサ１１は、外乱ｆの情報を、Ｓ２１０において予測式の定義情報の一部として取得する、又は、Ｓ２３０において目標指令ｒと共に取得することができる。
【００６７】
【数８】
2019175319.tif 000010
【００６８】
この予測式によれば、制御入力ｕに加算される外乱ｆを加味した制御入力ｕ及び制御出力ｙの予測が可能である。外乱ｆの例には、プラント５０に含まれる機構で生じる摩擦、負荷、及び外力が含まれる。本変形例に関する外乱ｆを加味した制御入力ｕ及び制御出力ｙの予測機能は、先に説明した変形例と同様、動作モードの一つとして、上記実施形態の予測装置１に組み込まれてもよい。
【００６９】
［第四変形例］
第四変形例の予測装置１は、プラントにおける伝達特性の変化δＱを加味して、制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出するように構成される。対応するフィードバック制御系１６０は、図１２Ａに示され、このフィードバック制御系１６０と等価な制御系１６５が図１２Ｂに示される。このフィードバック制御系１６０が備える、制御入力ｕと制御出力ｙとの関係を示すプラント８０の伝達関数Ｐは、次式で表される。
【００７０】
【数９】
2019175319.tif 000011
【００７１】
プラント８０の伝達特性の変化δＱを加味した予測は、Ｓ２２０で次式に従う予測式を組み立てることにより実現することができる。プロセッサ１１は、変化δＱの情報を、Ｓ２１０において予測式の定義情報の一部として取得する、又は、Ｓ２３０において目標指令ｒと共に取得することができる。
【００７２】
【数１０】
2019175319.tif 000012
【００７３】
この予測式によれば、プラント８０の伝達特性の変化δＱを加味した制御入力ｕ及び制御出力ｙの予測が可能である。プラント８０の伝達特性の変化δＱの例には、プラント８０内の粘性摩擦Ｄの変化が含まれる。モータの粘性摩擦が変化するプラント８０の伝達特性は、図１２Ｃに示される。この例によれば、プラント８０の伝達関数Ｐは、次式で表される。
【００７４】
【数１１】
2019175319.tif 000013
ここでＤは、粘性摩擦係数に対応し、Ｊは、慣性モーメントに対応し、ｓは、ラプラス演算子に対応する。この例によれば、設計者は、粘性摩擦の変化を加味した制御器３０の挙動を詳細に分析及び把握可能である。本変形例に関するプラントの伝達特性の変化δＱを加味した制御入力ｕ及び制御出力ｙの予測機能は、先に説明した変形例と同様、動作モードの一つとして、上記実施形態の予測装置１に組み込まれてもよい。
【００７５】
［第五変形例］
第五変形例の予測装置１は、プラントにおける伝達特性の変化（ｘ−１）Ｑを加味して、制御入力ｕの予測値ｕｐ及び制御出力ｙの予測値ｙｐを算出するように構成される。変数ｘが変動成分に対応する。対応するフィードバック制御系１７０は、図１３Ａに示され、このフィードバック制御系１７０と等価な制御系１７５が図１３Ｂに示される。
【００７６】
フィードバック制御系１７０におけるプラント９０の伝達特性の変化（ｘ−１）Ｑを加味した予測は、Ｓ２２０で次式に従う予測式を組み立てることにより実現することができる。プロセッサ１１は、変化（ｘ−１）Ｑの情報を、Ｓ２１０において予測式の定義情報の一部として取得する、又は、Ｓ２３０において目標指令ｒと共に取得することができる。
【００７７】
【数１２】
2019175319.tif 000014
【００７８】
この予測式によれば、プラント９０の伝達特性の変化（ｘ−１）Ｑを加味した制御入力ｕ及び制御出力ｙの予測が可能である。本変形例の上記変化（ｘ−１）Ｑは、第四変形例の変化δＱに対応する。本変形例に関するプラント９０の伝達特性の変化（ｘ−１）Ｑを考慮した制御入力ｕ及び制御出力ｙの予測機能は、先に説明した変形例と同様、動作モードの一つとして、上記実施形態の予測装置１に組み込まれてもよい。
【００７９】
［その他］
本開示は、上述した変形例を含む予測装置１の実施形態に限定されるものではなく、種々の態様を採ることができる。
【００８０】
予測装置は、コンピュータプログラムに基づく処理を、汎用のコンピュータに実行させることにより実現されなくてもよい。即ち、予測装置は、専用のハードウェア回路で対応する処理を実行するように構成されてもよい。特許請求の範囲に記載の文言から特定される技術思想に含まれるあらゆる態様が本開示の実施形態である。
【００８１】
最後に用語間の対応関係を説明する。プロセッサ１１が実行するＳ１１０の処理が、取得ユニットにより実現される処理の一例に対応し、プロセッサ１１が実行するＳ１２０の処理が、第一算出ユニットにより実現される処理の一例に対応し、プロセッサ１１が実行するＳ１３０の処理が、第二算出ユニットにより実現される処理の一例に対応し、プロセッサ１１が実行する予測処理が、予測ユニットにより実現される処理の一例に対応する。
【符号の説明】
【００８２】
１…予測装置、１１…プロセッサ、１２…メモリ、１３…ストレージデバイス、１５…操作デバイス、１７…表示デバイス、１９…入出力インタフェース、３０，６０…制御器、５０，７０，８０，９０…プラント、１００，１１０，１２０，１３０，１４０，１５０，１６０，１７０…フィードバック制御系。
【図１】
2019175319.tif 000021
【図２】
2019175319.tif 000022
【図３】
2019175319.tif 000023
【図４】
2019175319.tif 000024
【図５】
2019175319.tif 000025
【図６】
2019175319.tif 000026
【図７】
2019175319.tif 000027
【図８】
2019175319.tif 000028
【図９】
2019175319.tif 000029
【図１０】
2019175319.tif 000030
【図１１】
2019175319.tif 000031
【図１２】
2019175319.tif 000032
【図１３】
2019175319.tif 000033