特開 2023-128179 パラメータ推定装置および方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023128179

(43)【公開日】2023-09-14

(54)【発明の名称】パラメータ推定装置および方法

(51)【国際特許分類】

   G05B 11/36 20060101AFI20230907BHJP

   G05B 13/02 20060101ALI20230907BHJP

【ＦＩ】

G05B11/36 501B

G05B13/02 C

【審査請求】未請求

【請求項の数】4

【出願形態】ＯＬ

(21)【出願番号】P 2022032348

(22)【出願日】2022-03-03

(71)【出願人】

【識別番号】000006666

【氏名又は名称】アズビル株式会社

(74)【代理人】

【識別番号】100098394

【弁理士】

【氏名又は名称】山川 茂樹

(72)【発明者】

【氏名】田中 雅人

【テーマコード（参考）】

5H004

【Ｆターム（参考）】

5H004GA03

5H004GA07

5H004HA01

5H004HB01

5H004KA69

5H004KB02

5H004KB04

5H004KB06

5H004KC53

5H004KC55

5H004LB05

(57)【要約】

【課題】オーバーシュート抑制機能のための適切なパラメータの決定・調整を効率化する。
【解決手段】パラメータ推定装置１は、制御装置２のオーバーシュート抑制機能が無効の状態において外乱印加前後の整定時における操作量ＭＶの変化分を検出する操作量差検出部１０と、オーバーシュート量の最大値を検出するオーバーシュート量検出部１１と、外乱偏差の最大値を検出する外乱偏差検出部１２と、制御装置２に設定されている比例帯を取得する比例帯取得部１３と、操作量ＭＶの変化分と比例帯に基づいて、オーバーシュート抑制のための制御装置２の偏差補正において確保されるべき補正後偏差を推定する補正後偏差推定部１４と、オーバーシュート量の最大値と外乱偏差の最大値との比率と、補正後偏差に基づいて、偏差補正のためのパラメータの適正値を算出する適正値算出部１５を備える。
【選択図】 図１

【特許請求の範囲】

【請求項1】

制御装置が設定値と制御量との偏差を補正して前記制御量のオーバーシュートを抑制する機能が無効の状態において、前記制御量に外乱が印加されたときに外乱印加前の整定時に前記制御装置から制御対象に出力されていた操作量に対して外乱印加後の整定時に出力された前記操作量の変化分を検出するように構成された操作量差検出部と、
前記外乱が印加された後に発生した前記制御量のオーバーシュート量の最大値を検出するように構成されたオーバーシュート量検出部と、
前記外乱が印加されたときの外乱偏差の最大値を検出するように構成された外乱偏差検出部と、
前記制御装置に設定されているＰＩＤパラメータ中の比例帯を取得するように構成された比例帯取得部と、
前記操作量の変化分と前記比例帯とに基づいて、オーバーシュート抑制のための前記制御装置の偏差補正において確保されるべき補正後偏差を推定するように構成された補正後偏差推定部と、
前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記補正後偏差とに基づいて、前記偏差補正のためのパラメータξの適正値を算出するように構成された適正値算出部と、
前記パラメータξの適正値を出力するように構成された適正値出力部とを備えることを特徴とするパラメータ推定装置。

【請求項2】

請求項１記載のパラメータ推定装置において、
前記適正値算出部は、前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記パラメータξとに基づいて、前記制御装置において前記偏差補正と等価な設定値補正を行った後の前記設定値に対する前記制御量のオーバーシュート量を推定し、このオーバーシュート量が、前記偏差に前記パラメータξを乗算した偏差補正量に対して規定された割合以下になるように前記パラメータξを修正することを特徴とするパラメータ推定装置。

【請求項3】

制御装置が設定値と制御量との偏差を補正して前記制御量のオーバーシュートを抑制する機能が無効の状態において、前記制御量に外乱が印加されたときに外乱印加前の整定時に前記制御装置から制御対象に出力されていた操作量に対して外乱印加後の整定時に出力された前記操作量の変化分を検出する第１のステップと、
前記外乱が印加された後に発生した前記制御量のオーバーシュート量の最大値を検出する第２のステップと、
前記外乱が印加されたときの外乱偏差の最大値を検出する第３のステップと、
前記制御装置に設定されているＰＩＤパラメータ中の比例帯を取得する第４のステップと、
前記操作量の変化分と前記比例帯とに基づいて、オーバーシュート抑制のための前記制御装置の偏差補正において確保されるべき補正後偏差を推定する第５のステップと、
前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記補正後偏差とに基づいて、前記偏差補正のためのパラメータの適正値を算出する第６のステップと、
前記パラメータの適正値を出力する第７のステップとを含むことを特徴とするパラメータ推定方法。

【請求項4】

請求項３記載のパラメータ推定方法において、
前記第６のステップは、前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記パラメータξとに基づいて、前記制御装置において前記偏差補正と等価な設定値補正を行った後の前記設定値に対する前記制御量のオーバーシュート量を推定し、このオーバーシュート量が、前記偏差に前記パラメータξを乗算した偏差補正量に対して規定された割合以下になるように前記パラメータξを修正するステップを含むことを特徴とするパラメータ推定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、制御装置のパラメータを推定する技術に係り、特にオーバーシュート抑制機能のためのパラメータを推定するパラメータ推定装置および方法に関するものである。

【背景技術】

【0002】

代表的なフィードバック制御であるＰＩＤ制御の動作には、例えば温度制御系における外乱印加時のリカバリー制御（外乱リカバリー応答）があり、加熱装置などの保温性が高いことなどの要因により、制御量にオーバーシュートが発生することが多い。

【0003】

図１０の加熱装置は、処理対象のワークを加熱する熱処理炉１００と、電気ヒータ１０１と、熱処理炉１００内の温度を計測する温度センサ１０２と、熱処理炉１００内の温度を制御する温調計１０３と、電力調整器１０４と、電力供給回路１０５と、加熱装置全体を制御するＰＬＣ（Programmable Logic Controller）１０６とから構成される。温調計１０３は、温度センサ１０２が計測した温度ＰＶ（制御量）が温度設定値ＳＰと一致するように操作量ＭＶを算出する。電力調整器１０４は、操作量ＭＶに応じた電力を決定し、この決定した電力を電力供給回路１０５を通じて電気ヒータ１０１に供給する。

【0004】

図１１は、図１０に示したような加熱装置の温度制御中に発生したオーバーシュートの例を示す図である。図１１の例では、温度制御中に制御量ＰＶが下降する外乱が発生し、外乱に応じた操作量ＭＶが温調計１０３から出力され外乱リカバリー制御が行われたときに、過剰な操作量修正によって制御量ＰＶが温度設定値ＳＰ＝３００℃を超過するオーバーシュートが発生した例を示している。

【0005】

図１１のような現象に対応するため、オーバーシュート抑制機能が提案されている（特許文献１参照）。特許文献１に開示された技術では、外乱印加時の前半（図１１の降温時）はパラメータξにより設定値ＳＰと制御量ＰＶとの偏差Ｅｒ（＝ＳＰ－ＰＶ）の補正量Ｅｒｘを算出し、外乱印加時の後半（収束時）はパラメータλにより補正量Ｅｒｘを算出する。

【0006】

具体的には、偏差Ｅｒが正または負の値で、その絶対値｜Ｅｒ｜が増大するときに外乱が印加されたと判定し、偏差補正量Ｅｒｘを次式のように算出する。
Ｅｒｘ＝ξＥｒ ・・・（１）

【0007】

パラメータξは例えば０．８である。続いて、偏差補正量Ｅｒｘが一定規則で徐々に０に収束するように収束演算を行う。具体的には、収束演算後の偏差補正量Ｅｒｘ’を次式のように算出する。
Ｅｒｘ’＝λＥｒｘ ・・・（２）

【0008】

パラメータλ（０＜λ＜１）は例えばλ＝０．９５である。以上のような偏差補正量Ｅｒｘ’によって偏差Ｅｒを補正（Ｅｒ－Ｅｒｘ’）することにより、オーバーシュートを抑制する。特許文献１に開示された技術の効果を図１２に示す。なお、図１２では、偏差補正量Ｅｒｘ’の代わりに、補正後の設定値ＳＰｈ＝ＳＰ－Ｅｒｘ’を示している。特許文献１に開示されているように、Ｅｒ－Ｅｒｘ’という偏差補正は、（ＳＰ－Ｅｒｘ’）－ＰＶという設定値補正と等価である。

【0009】

特許文献１に開示された技術では、パラメータξを定数としている。対象とする制御系に特有の適切なパラメータξの数値を決定できれば、頻繁に数値変更をする必要はない。しかしながら、特許文献１に開示された従来では、適切なパラメータξを決定する技術が実現できていない、という課題があった。

【先行技術文献】

【特許文献】

【0010】

【特許文献1】特許第３４３７８０７号公報

【発明の概要】

【発明が解決しようとする課題】

【0011】

本発明は、上記課題を解決するためになされたもので、制御装置のオーバーシュート抑制機能のための適切なパラメータの決定・調整を効率化することができるパラメータ推定装置および方法を提供することを目的とする。

【課題を解決するための手段】

【0012】

本発明のパラメータ推定装置は、制御装置が設定値と制御量との偏差を補正して前記制御量のオーバーシュートを抑制する機能が無効の状態において、前記制御量に外乱が印加されたときに外乱印加前の整定時に前記制御装置から制御対象に出力されていた操作量に対して外乱印加後の整定時に出力された前記操作量の変化分を検出するように構成された操作量差検出部と、前記外乱が印加された後に発生した前記制御量のオーバーシュート量の最大値を検出するように構成されたオーバーシュート量検出部と、前記外乱が印加されたときの外乱偏差の最大値を検出するように構成された外乱偏差検出部と、前記制御装置に設定されているＰＩＤパラメータ中の比例帯を取得するように構成された比例帯取得部と、前記操作量の変化分と前記比例帯とに基づいて、オーバーシュート抑制のための前記制御装置の偏差補正において確保されるべき補正後偏差を推定するように構成された補正後偏差推定部と、前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記補正後偏差とに基づいて、前記偏差補正のためのパラメータξの適正値を算出するように構成された適正値算出部と、前記パラメータξの適正値を出力するように構成された適正値出力部とを備えることを特徴とするものである。

【0013】

また、本発明のパラメータ推定装置の１構成例において、前記適正値算出部は、前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記パラメータξとに基づいて、前記制御装置において前記偏差補正と等価な設定値補正を行った後の前記設定値に対する前記制御量のオーバーシュート量を推定し、このオーバーシュート量が、前記偏差に前記パラメータξを乗算した偏差補正量に対して規定された割合以下になるように前記パラメータξを修正することを特徴とするものである。

【0014】

また、本発明のパラメータ推定方法は、制御装置が設定値と制御量との偏差を補正して前記制御量のオーバーシュートを抑制する機能が無効の状態において、前記制御量に外乱が印加されたときに外乱印加前の整定時に前記制御装置から制御対象に出力されていた操作量に対して外乱印加後の整定時に出力された前記操作量の変化分を検出する第１のステップと、前記外乱が印加された後に発生した前記制御量のオーバーシュート量の最大値を検出する第２のステップと、前記外乱が印加されたときの外乱偏差の最大値を検出する第３のステップと、前記制御装置に設定されているＰＩＤパラメータ中の比例帯を取得する第４のステップと、前記操作量の変化分と前記比例帯とに基づいて、オーバーシュート抑制のための前記制御装置の偏差補正において確保されるべき補正後偏差を推定する第５のステップと、前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記補正後偏差とに基づいて、前記偏差補正のためのパラメータの適正値を算出する第６のステップと、前記パラメータの適正値を出力する第７のステップとを含むことを特徴とするものである。

【0015】

また、本発明のパラメータ推定方法の１構成例において、前記第６のステップは、前記オーバーシュート量の最大値と前記外乱偏差の最大値との比率と、前記パラメータξとに基づいて、前記制御装置において前記偏差補正と等価な設定値補正を行った後の前記設定値に対する前記制御量のオーバーシュート量を推定し、このオーバーシュート量が、前記偏差に前記パラメータξを乗算した偏差補正量に対して規定された割合以下になるように前記パラメータξを修正するステップを含むことを特徴とするものである。

【発明の効果】

【0016】

本発明によれば、操作量差検出部とオーバーシュート量検出部と外乱偏差検出部と比例帯取得部と補正後偏差推定部と適正値算出部と適正値出力部とを設けることにより、制御装置のオーバーシュート抑制機能のための適切なパラメータξの決定・調整を効率化することができる。

【図面の簡単な説明】

【0017】

【図1】図１は、本発明の実施例に係るパラメータ推定装置の構成を示すブロック図である。

【図2】図２は、本発明の実施例に係る制御装置の動作を説明するフローチャートである。

【図3】図３は、本発明の実施例に係るパラメータ推定装置の動作を説明するフローチャートである。

【図4】図４は、従来の技術を適用した制御装置による外乱リカバリー応答のシミュレーション結果を示す図である。

【図5】図５は、従来のオーバーシュート抑制機能を無効にした制御装置による外乱リカバリー応答のシミュレーション結果を示す図である。

【図6】図６は、パラメータを適正値に調整した後の制御装置による外乱リカバリー応答のシミュレーション結果を示す図である。

【図7】図７は、従来のオーバーシュート抑制機能を無効にした制御装置による外乱リカバリー応答のシミュレーション結果を示す図である。

【図8】図８は、パラメータを適正値に調整した後の制御装置による外乱リカバリー応答のシミュレーション結果を示す図である。

【図9】図９は、本発明の実施例に係るパラメータ推定装置と制御装置を実現するコンピュータの構成例を示すブロック図である。

【図10】図１０は、加熱装置の構成を示すブロック図である。

【図11】図１１は、加熱装置の温度制御中に発生したオーバーシュートの例を示す図である。

【図12】図１２は、従来のオーバーシュート抑制機能の効果を示す図である。

【発明を実施するための形態】

【0018】

［発明の原理］
以下の説明では、設定値ＳＰを補正するという解釈で発明の原理を説明する。設定値ＳＰを補正する動作として捉えると、特許文献１に開示された技術の動作は、制御量ＰＶの下降に同調して、パラメータξで規定される比率で補正された設定値ＳＰｇ（＝ＳＰ－Ｅｒｘ）も下降する。すなわち、実際の偏差ＥｒよりもＰＩＤ演算に入力される偏差が少なくなり、偏差の減少に対応して操作量ＭＶの変化分も少なくなる。このとき、オーバーシュートに寄与してしまう操作量ＭＶの分だけ少なくなるのが適切と考えると、外乱リカバリー応答の前後での整定時の操作量ＭＶの変化分ΔＭＶが残るような設定値ＳＰの補正が好ましい。

【0019】

外乱には、大きく分けてステップ外乱とインパルス外乱があり、ステップ外乱では操作量ＭＶの変化分ΔＭＶが発生し、インパルス外乱では発生しない。この両者の混合がパラメータξの決定を難しくしているのであるから、少なくとも１回の外乱リカバリー応答を試行し、その際の操作量ＭＶの変化分ΔＭＶに応じてパラメータξを決定するように規定すればよい。

【0020】

なお、本発明では、例えばパラメータλを適切に調整した後に、パラメータξに関わる制御対象や外乱の性質の変化があった場合などを想定しているが、パラメータξを調整した後に、パラメータλを試行錯誤で調整することも可能である。

【0021】

［実施例］
以下、本発明の実施例について図面を参照して説明する。図１は本発明の実施例に係るパラメータ推定装置の構成を示すブロック図である。パラメータ推定装置１は、制御装置２が設定値ＳＰと制御量ＰＶとの偏差Ｅｒを補正して制御量ＰＶのオーバーシュートを抑制する機能が無効の状態において、制御量ＰＶに外乱が印加されたときに外乱印加前の整定時に制御装置２から制御対象に出力されていた操作量ＭＶに対して外乱印加後の整定時に出力された操作量ＭＶの変化分ΔＭＶを検出する操作量差検出部１０と、外乱が印加された後に発生した制御量ＰＶのオーバーシュート量の最大値ＯＳ＿ｍａｘを検出するオーバーシュート量検出部１１と、外乱が印加されたときの外乱偏差の最大値Ｅｒ＿ｍａｘを検出する外乱偏差検出部１２と、制御装置２に設定されているＰＩＤパラメータ中の比例帯Ｐｂを取得する比例帯取得部１３と、操作量ＭＶの変化分ΔＭＶと比例帯Ｐｂとに基づいて、オーバーシュート抑制のための制御装置２の偏差補正において確保されるべき補正後偏差Ｅｒｈを推定する補正後偏差推定部１４と、オーバーシュート量の最大値ＯＳ＿ｍａｘと外乱偏差の最大値Ｅｒ＿ｍａｘとの比率と、補正後偏差Ｅｒｈとに基づいて、偏差補正のためのパラメータξの適正値を算出する適正値算出部１５と、パラメータξの適正値を出力する適正値出力部１６とを備えている。

【0022】

制御装置２は、制御量ＰＶを入力とする制御量入力部２０と、設定値ＳＰを入力とする設定値入力部２１と、設定値ＳＰから制御量ＰＶを減算して偏差Ｅｒを算出する減算部２２と、偏差Ｅｒに基づいて外乱が印加されたか否かを制御周期毎に判定し、偏差Ｅｒの絶対値が増大するとき外乱が印加されたと判定する外乱印加検出部２３と、外乱印加検出部２３によって外乱の印加が検出されたときに偏差Ｅｒの大きさに基づいて偏差補正量Ｅｒｘを算出する偏差補正量算出部２４と、偏差補正量算出部２４によって算出された偏差補正量Ｅｒｘが一定規則で徐々に０に収束するよう収束演算を行う偏差補正量収束算出部２５と、偏差Ｅｒと収束演算後の偏差補正量とに基づき操作量ＭＶをＰＩＤ制御演算によって算出する制御演算部２６と、操作量ＭＶを制御対象に出力する操作量出力部２７とを備えている。

【0023】

まず、パラメータ推定装置１について説明する前に、特許文献１に開示された制御装置２の動作について簡単に説明する。図２は制御装置２の動作を説明するフローチャートである。制御量入力部２０には、図示しないセンサによって計測された制御量ＰＶが入力される（図２ステップＳ１０１）。

【0024】

設定値入力部２１には、オペレータによって設定された設定値ＳＰが入力される（図２ステップＳ１０２）。
減算部２２は、設定値ＳＰと制御量ＰＶとの偏差Ｅｒ＝ＳＰ－ＰＶを算出する（図２ステップＳ１０３）。

【0025】

次に、外乱印加検出部２３は、外乱が印加されたか否かを制御周期毎に判定し、外乱が印加されたと判定したときに偏差補正量算出部２４に起動信号ｓ１を出力する（図２ステップＳ１０４）。具体的には、外乱印加検出部２３は、偏差Ｅｒが正の値で、絶対値｜Ｅｒ｜が増大する場合（Ｅｒ＞０で、かつＥｒ＞Ｅｒ’（Ｅｒ’は１制御周期前の偏差）が成立する場合）、または偏差Ｅｒが負の値で、絶対値｜Ｅｒ｜が増大する場合（Ｅｒ＜０で、かつＥｒ＜Ｅｒ’が成立する場合）に、外乱が印加されたと判定する。

【0026】

偏差補正量算出部２４は、外乱印加検出部２３から起動信号ｓ１が出力されると（ステップＳ１０４においてＹＥＳ）、式（１）のように偏差Ｅｒにパラメータξを乗算した偏差補正量Ｅｒｘを算出し、これを偏差補正量収束算出部２５に出力する（図２ステップＳ１０５）。

【0027】

次に、偏差補正量収束算出部２５は、収束演算後の偏差補正量Ｅｒｘ’を式（２）のように算出する（図２ステップＳ１０６）。パラメータλは定数（０＜λ＜１）である。なお、偏差補正量算出部２４は、外乱が印加されたと外乱印加検出部２３が判定したときのみ偏差補正量Ｅｒｘを出力する。すなわち、偏差の絶対値｜Ｅｒ｜が増大する状況（制御量ＰＶが設定値ＳＰから離れる状況）が終了すると、起動信号ｓ１の出力が停止し、偏差補正量算出部２４による偏差補正量Ｅｒｘの出力が停止する。

【0028】

偏差補正量収束算出部２５は、偏差補正量算出部２４から偏差補正量Ｅｒｘが出力されない場合、１制御周期前に収束演算を行った偏差補正量Ｅｒｘ’を現制御周期の偏差補正量Ｅｒｘとして式（２）の収束演算を行う。偏差補正量Ｅｒｘ’（Ｅｒｘ）の初期値は０である。

【0029】

制御演算部２６は、偏差Ｅｒから収束演算後の偏差補正量Ｅｒｘ’を減算した値に基づいて次式のようなＰＩＤ制御演算を行い、操作量ＭＶを算出する（図２ステップＳ１０７）。
ＭＶ＝（１００／Ｐｂ）｛１＋（１／Ｔｉｓ）＋Ｔｄｓ｝（Ｅｒ－Ｅｒｘ’）
・・・（３）

【0030】

式（３）において、Ｐｂは比例帯、Ｔｉは積分時間、Ｔｄは微分時間、ｓはラプラス演算子である。
制御演算部２６によって算出された操作量ＭＶは、操作量出力部２７を介して制御対象（実際には例えば図１０の電力調整器１０４）へ出力される（図２ステップＳ１０８）。

【0031】

制御装置２は、ステップＳ１０１～Ｓ１０８の処理を例えばオペレータの指示によって制御が終了するまで（図２ステップＳ１０９においてＹＥＳ）、制御周期毎に実行する。

【0032】

次に、本実施例のパラメータ推定装置１の動作を図３を用いて説明する。パラメータ推定装置１は、パラメータξの調整時に制御装置２と接続される。
最初に、操作量差検出部１０は、制御装置２のオーバーシュート抑制機能を無効にする（図３ステップＳ２００）。オーバーシュート抑制機能が無効の場合、偏差補正量Ｅｒｘ’が０になるので、制御装置２の制御演算部２６は、式（３）の代わりに次の式（４）により操作量ＭＶを算出する。
ＭＶ＝（１００／Ｐｂ）｛１＋（１／Ｔｉｓ）＋Ｔｄｓ｝Ｅｒ ・・・（４）

【0033】

操作量差検出部１０は、制御装置２のオーバーシュート抑制機能が無効の状態で外乱が印加されたと制御装置２の外乱印加検出部２３が判定した後に（図３ステップＳ２０１においてＹＥＳ）、外乱リカバリー応答によって制御量ＰＶが整定したときに（図３ステップＳ２０２においてＹＥＳ）、外乱印加前の整定時の操作量ＭＶに対する外乱印加後の整定時の操作量ＭＶの変化分ΔＭＶを検出する（図３ステップＳ２０３）。

【0034】

なお、例えば図１０に示した加熱装置の場合であれば、パラメータξを決定するための試行時に、熱処理炉１００の扉を一時的に開くことで、熱処理炉１００内の温度（制御量ＰＶ）が降温する外乱を意図的に発生させることが可能である。

【0035】

オーバーシュート量検出部１１は、制御装置２のオーバーシュート抑制機能が無効の状態で、外乱が印加された後に外乱リカバリー応答によって発生した制御量ＰＶのオーバーシュート量の最大値ＯＳ＿ｍａｘ（上昇した制御量ＰＶの最大値ＰＶ＿ｍａｘと設定値ＳＰとの差ＰＶ＿ｍａｘ－ＳＰ）を検出する（図３ステップＳ２０４）。

【0036】

外乱偏差検出部１２は、制御装置２のオーバーシュート抑制機能が無効の状態で、外乱が印加されたときの外乱偏差の最大値Ｅｒ＿ｍａｘ（設定値ＳＰと下降した制御量ＰＶの最小値ＰＶ＿ｍｉｎとの差ＳＰ－ＰＶ＿ｍｉｎ）を検出する（図３ステップＳ２０５）。
比例帯取得部１３は、制御装置２の制御演算部２６に設定されているＰＩＤパラメータ中の比例帯Ｐｂを取得する（図３ステップＳ２０６）。

【0037】

補正後偏差推定部１４は、操作量差検出部１０によって検出された操作量ＭＶの変化分ΔＭＶと比例帯取得部１３によって取得された比例帯Ｐｂとに基づいて、オーバーシュート抑制のための偏差Ｅｒ補正において確保されるべき補正後偏差Ｅｒｈを推定する（図３ステップＳ２０７）。具体的には、補正後偏差Ｅｒｈは、制御装置２の、比例帯Ｐｂで規定される比例動作のみによって操作量変化分ΔＭＶを発生させるために必要な偏差である。補正後偏差推定部１４は、補正後偏差Ｅｒｈを次式により算出する。
Ｅｒｈ＝ΔＭＶ（Ｐｂ／１００） ・・・（５）

【0038】

適正値算出部１５は、外乱偏差検出部１２によって検出された外乱偏差の最大値Ｅｒ＿ｍａｘと補正後偏差推定部１４によって算出された補正後偏差Ｅｒｈとに基づいて、パラメータξの適正値候補（暫定値）ξａを次式により算出する（図３ステップＳ２０８）。
ξａ＝（Ｅｒ＿ｍａｘ－Ｅｒｈ）／Ｅｒ＿ｍａｘ ・・・（６）

【0039】

式（６）は、偏差補正量Ｅｒｘ＝ξＥｒ＿ｍａｘ＝Ｅｒ＿ｍａｘ－Ｅｒｈの関係から得られたものである。
続いて、適正値算出部１５は、パラメータξの適正値候補ξａにより補正後の設定値ＳＰｇが想定できるので、適正値候補ξａとオーバーシュート量検出部１１によって検出されたオーバーシュート量の最大値ＯＳ＿ｍａｘと外乱偏差検出部１２によって検出された外乱偏差の最大値Ｅｒ＿ｍａｘとに基づいて、偏差補正と等価な設定値補正を行った後の設定値ＳＰｇに対する制御量ＰＶのオーバーシュート量ＯＳｘを次式により推定する（図３ステップＳ２０９）。
ＯＳｘ＝ＥｒｈＯＳ＿ｍａｘ／Ｅｒ＿ｍａｘ
＝（Ｅｒ＿ｍａｘ－ξａＥｒ＿ｍａｘ）ＯＳ＿ｍａｘ／Ｅｒ＿ｍａｘ ・・（７）

【0040】

式（７）は、制御装置２のＰＩＤ制御系の特性に変化がない場合に、外乱偏差とオーバーシュート量との比率（オーバーシュート率）が維持される（Ｅｒ＿ｍａｘ：ＯＳ＿ｍａｘ＝Ｅｒｈ：ＯＳｘ）という、性質から得られたものである。

【0041】

なお、式（７）の計算で想定する補正後の設定値ＳＰｇとは、制御装置２が操作量ＭＶの算出の際に用いる補正後偏差に相当する補正後の設定値ＳＰｈ＝ＳＰ－Ｅｒｘ’そのものではなく、偏差補正量Ｅｒｘによって設定値ＳＰを補正した値ＳＰｇ＝ＳＰ－Ｅｒｘである。したがって、上記の「偏差補正と等価な設定値補正」の偏差補正とは、Ｅｒ－Ｅｒｘ’ではなく、Ｅｒ－Ｅｒｘを意味する。

【0042】

算出されたオーバーシュート量ＯＳｘが偏差補正量Ｅｒｘ＝ξａＥｒ＿ｍａｘよりも大きい場合（図１２のケースで言えば、オーバーシュート量が補正後の設定値ＳＰｈの降下量よりも大きいことで補正前の設定値ＳＰを制御量ＰＶが超える場合）、オーバーシュート抑制が不十分である。さらに厳密に言えば、偏差補正量Ｅｒｘが補正前の設定値ＳＰに収束する途中でオーバーシュート量ＯＳｘが発生するようになる。設定値ＳＰに収束する途中とは、パラメータλを時定数Ｔに等価変換した場合の時定数Ｔの概ね１．５倍くらいの時間帯になる。パラメータλと時定数Ｔとの関係は式（８）のようになる。ｄｔは制御装置２の制御周期である。
λ＝Ｔ／（Ｔ＋ｄｔ） ・・・（８）

【0043】

したがって、適正値算出部１５は、オーバーシュート量ＯＳｘが偏差補正量Ｅｒｘに対して規定された割合を超える場合、具体的にはＯＳｘ＞Ｅｒｘ／（２．７１８×１．５）が成立する場合（図３ステップＳ２１０においてＹＥＳ）、次式によりパラメータξの適正値候補ξａを修正する（図３ステップＳ２１１）。
ξａ＿ｎｅｗ＝ξａ＋０．０１ ・・・（９）

【0044】

上記の２．７１８は、偏差補正量Ｅｒｘの収束量に対応する値であり、指数関数で用いられるネイピア数である。
適正値算出部１５は、算出したξａ＿ｎｅｗを新たな適正値候補ξａとして、ステップＳ２０９のオーバーシュート量ＯＳｘの算出処理に戻る。こうして、適正値算出部１５は、オーバーシュート量ＯＳｘが偏差補正量Ｅｒｘに対して規定された割合以下になるまで、すなわちＯＳｘ≦Ｅｒｘ／（２．７１８×１．５）が成立するまで（ステップＳ２１０においてＮＯ）、ステップＳ２０９～Ｓ２１１の処理を繰り返す。

【0045】

適正値算出部１５は、ＯＳｘ≦Ｅｒｘ／（２．７１８×１．５）が成立した時点の適正値候補ξａをパラメータξの適正値として確定する。
適正値出力部１６は、パラメータξの適正値を出力する（図３ステップＳ２１２）。出力の例としては、適正値の表示などがある。オペレータは、表示された適正値を見て、制御装置２に設定されているパラメータξを変更することになる。あるいは、適正値出力部１６は、制御装置２の偏差補正量算出部２４に設定されているパラメータξを適正値に自動的に更新するようにしてもよい。

【0046】

なお、パラメータξの上限値をξ＿ｍａｘ＝０．９９のように予め規定しておくことが好ましい。この場合、適正値算出部１５は、ＯＳｘ≦Ｅｒｘ／（２．７１８×１．５）が成立する前に式（９）のξａ＿ｎｅｗが上限値ξ＿ｍａｘに達した時点でステップＳ２０９～Ｓ２１１の探索処理を終了し、上限値ξ＿ｍａｘをパラメータξの適正値として確定する。

【0047】

以上で、パラメータξの調整が終了する。調整終了後、オペレータは、パラメータ推定装置１を操作して、制御装置２のオーバーシュート抑制機能を有効に戻し、パラメータ推定装置１と制御装置２との接続を解除する（図３ステップＳ２１３）。

【0048】

制御量ＰＶにインパルス外乱の傾向が強い変化（操作量ＭＶの変化分ΔＭＶ＝１．０％）が発生した場合の本実施例の検証例を図４～図６に示す。なお、この検証例では、制御装置２のパラメータλに対応する時定数ＴがＴ＝１２０秒で概ね妥当な値に調整されていることを前提とする。ただし、パラメータ推定装置１によってパラメータξを調整した後に、時定数Ｔ（パラメータλ）を試行錯誤で微調整しても構わない。

【0049】

図４は、特許文献１に開示された技術を適用した制御装置２による外乱リカバリー応答のシミュレーション結果を示す図である。図４の例では、パラメータλに対応する時定数ＴがＴ＝１２０秒で概ね妥当な値に設定されているが、パラメータξが不適切な値ξ＝０．５０に設定されていることにより、制御量ＰＶが下降する外乱が発生したときに、過剰な操作量修正によって制御量ＰＶが設定値ＳＰ＝３００℃を超過するオーバーシュートが発生している。

【0050】

ここでは、パラメータξと時定数Ｔがもともと適切に調整されていたものとし、その後、インパルス外乱の傾向とステップ外乱の傾向のバランスが変化したことにより、パラメータξが不適切になったものと仮定する。すなわち、調整すべき対象はパラメータξであるとする。

【0051】

図５は、図４と同じ条件で、特許文献１に開示されたオーバーシュート抑制機能を無効にした制御装置２による外乱リカバリー応答のシミュレーション結果を示す図である。すなわち、図５の例は、ステップＳ２００～ステップＳ２０５で実施された外乱リカバリー応答に相当する。図５の例は操作量ＭＶの変化分がΔＭＶ＝１．０％と小さく、インパルス外乱の傾向がある。

【0052】

本実施例のオーバーシュート量検出部１１はオーバーシュート量の最大値ＯＳ＿ｍａｘ＝８℃を検出し、外乱偏差検出部１２は外乱偏差の最大値Ｅｒ＿ｍａｘ＝２０℃を検出する。制御装置２の制御演算部２６に設定されている比例帯Ｐｂは１００％である。本実施例の適正値算出部１５は、パラメータξの適正値をξ＝０．９５と算出する。

【0053】

図６は、パラメータξを適正値に調整した後の制御装置２による外乱リカバリー応答のシミュレーション結果を示す図である。図６の例では、パラメータλに対応する時定数ＴがＴ＝１２０秒で概ね妥当な値に設定され、さらにパラメータξを適切な値ξ＝０．９５に設定したことにより、制御量ＰＶのオーバーシュートが適切に抑制されていることが分かる。

【0054】

制御量ＰＶにステップ外乱の傾向が強い変化（操作量ＭＶの変化分ΔＭＶ＝３０．０％）が発生した場合の本実施例の検証例を図７、図８に示す。なお、この検証例では、制御装置２のパラメータλに対応する時定数ＴがＴ＝２３０秒で概ね妥当な値に調整されていることを前提とする。

【0055】

図７は、特許文献１に開示されたオーバーシュート抑制機能を無効にした制御装置２による外乱リカバリー応答のシミュレーション結果を示す図である。すなわち、図７の例は、ステップＳ２００～ステップＳ２０５で実施された外乱リカバリー応答に相当する。図７の例は操作量ＭＶの変化分がΔＭＶ＝３０．０％と大きく、ステップ外乱の傾向がある。

【0056】

本実施例のオーバーシュート量検出部１１はオーバーシュート量の最大値ＯＳ＿ｍａｘ＝５．５℃を検出し、外乱偏差検出部１２は外乱偏差の最大値Ｅｒ＿ｍａｘ＝２３℃を検出する。制御装置２の制御演算部２６に設定されている比例帯Ｐｂは１００％である。本実施例の適正値算出部１５は、パラメータξの適正値をξ＝０．５０と算出する。

【0057】

図８は、パラメータξを適正値に調整した後の制御装置２による外乱リカバリー応答のシミュレーション結果を示す図である。図８の例では、パラメータλに対応する時定数ＴがＴ＝２３０秒で概ね妥当な値に設定され、さらにパラメータξを適切な値ξ＝０．５０に設定したことにより、制御量ＰＶのオーバーシュートが適切に抑制されていることが分かる。

【0058】

なお、本実施例では、パラメータ推定装置１と制御装置２とを別々に設けているが、制御装置２の内部にパラメータ推定装置１を設けるようにしてもよい。

【0059】

本実施例のパラメータ推定装置１と制御装置２の各々は、ＣＰＵ（Central Processing Unit）、記憶装置および外部とのインタフェースを備えたコンピュータと、これらのハードウェア資源を制御するプログラムによって実現することができる。このコンピュータの構成例を図９に示す。コンピュータは、ＣＰＵ２００と、記憶装置２０１と、インタフェース装置（Ｉ／Ｆ）２０２とを備えている。

【0060】

パラメータ推定装置１の場合、Ｉ／Ｆ２０２には、制御装置２とディスプレイ装置等が接続される。制御装置２の場合、Ｉ／Ｆ２０２には、パラメータ推定装置１とセンサと電力調整器等が接続される。本発明のパラメータ推定方法を実現させるためのプログラムは記憶装置２０１に格納される。各装置のＣＰＵ２００は、記憶装置２０１に格納されたプログラムに従って本実施例で説明した処理を実行する。

【産業上の利用可能性】

【0061】

本発明は、制御装置のパラメータを推定する技術に適用することができる。

【符号の説明】

【0062】

１…パラメータ推定装置、２…制御装置、１０…操作量差検出部、１１…オーバーシュート量検出部、１２…外乱偏差検出部、１３…比例帯取得部、１４…補正後偏差推定部、１５…適正値算出部、１６…適正値出力部、２０…制御量入力部、２１…設定値入力部、２２…減算部、２３…外乱印加検出部、２４…偏差補正量算出部、２５…偏差補正量収束算出部、２６…制御演算部、２７…操作量出力部。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】