(57)【要約】【課題】オーバシュートを抑制する。【解決手段】制御装置１は、ステップ応答制御手段４ａが２回のステップ応答の制御を行い、状態判定手段４ｂが２回のステップ応答による偏差と偏差の変化量が基準以下のときに定常状態と判定する。平均値算出手段４ｃが定常状態時の制御量と操作量の平均値を算出し、むだ時間計測手段４ｄが定常状態時のむだ時間を計測し、実測波形取得手段４ｅが定常状態時の実測波形を取得する。モデル制御手段４ｆが伝達関数Ｇ（ｓ）をモデルとしてシミュレーションし、ゲイン算出手段４ｇが制御量と操作量の平均値からゲインを算出し、予測波形取得手段４ｈがモデル制御時の予測波形を取得する。評価値が最も小さくなる固有角周波数と減衰比の組み合わせをパラメータ算出手段４ｊが算出し、制御量が目標値以上かつ制御量と目標値との差の最大値が最も小さくなるＰ、Ｉ、Ｄの組み合わせをＰＩＤ定数算出手段４ｋが算出する。【選択図】図１
2020160550.tif 000002
【特許請求の範囲】
【請求項１】
制御対象の制御量を変化させるための操作量を算出して前記制御量を目標値に一致させる制御装置であって、
前記制御量が前記目標値よりも小さい中途目標値、前記目標値の順に到達するように２回のステップ応答の制御を行うステップ応答制御手段と、
前記２回のステップ応答の制御によって得られる前記目標値と前記中途目標値それぞれと前記制御量との差を示す偏差と該偏差の変化量が基準以下のときに定常状態と判定する状態判定手段と、
前記２回のステップ応答の制御で前記定常状態と判定されたときの前記制御量と前記操作量の平均値を算出する平均値算出手段と、
前記制御量が前記中途目標値に到達するようにステップ応答の制御をして前記定常状態と判定されたときに該定常状態から前記制御量が所定量変化するのに要するむだ時間を計測するむだ時間計測手段と、
前記２回のステップ応答の制御で前記定常状態と判定されたときの波形を実測波形として取得する実測波形取得手段と、
ラプラス演算子ｓ、ゲインＫ、むだ時間Ｔ、固有角周波数ωｎ、減衰比ζとしたときの伝達関数Ｇ（ｓ）＝（Ｋωｎ2 ／ｓ2 ＋２ζωｎｓ＋ωｎ2 ）・（１−Ｔｓ／１＋Ｔｓ）をモデルとし、前記固有角周波数と前記減衰比を正規分布乱数により変えて前記制御量のシミュレーションを行うモデル制御手段と、
前記平均値算出手段にて算出した前記制御量と前記操作量の平均値に基づいてゲインを算出するゲイン算出手段と、
前記モデルにて制御を行ったときの波形を予測波形として取得する予測波形取得手段と、
前記実測波形と前記予測波形との波形の誤差の絶対値を加算した評価値が最も小さくなる前記固有角周波数と前記減衰比の組み合わせを算出するパラメータ算出手段と、
比例帯Ｐ、積分時間Ｉ、微分時間Ｄを正規分布乱数により変えて前記制御量のシミュレーションを行ったときの前記制御量が前記目標値以上かつ前記制御量と前記目標値との差の最大値が最も小さくなる比例帯Ｐ、積分時間Ｉ、微分時間Ｄの組み合わせを算出するＰＩＤ定数算出手段と、
を具備することを特徴とする制御装置。
【請求項２】
前記制御対象の特性に応じて初期時のＰＩＤ定数を設定することを特徴とする請求項１に記載の制御装置。
【請求項３】
制御対象の制御量を変化させるための操作量を算出して前記制御量を目標値に一致させる制御方法であって、
前記制御量が前記目標値よりも小さい中途目標値、前記目標値の順に到達するように２回のステップ応答の制御を行うステップと、
前記２回のステップ応答の制御によって得られる前記目標値と前記中途目標値それぞれと前記制御量との差を示す偏差と該偏差の変化量が基準以下のときに定常状態と判定するステップと、
前記２回のステップ応答の制御で前記定常状態と判定されたときの前記制御量と前記操作量の平均値を算出するステップと、
前記制御量が前記中途目標値に到達するようにステップ応答の制御をして前記定常状態と判定されたときに該定常状態から前記制御量が所定量変化するのに要するむだ時間を計測するステップと、
前記２回のステップ応答の制御で前記定常状態と判定されたときの波形を実測波形として取得するステップと、
ラプラス演算子ｓ、ゲインＫ、むだ時間Ｔ、固有角周波数ωｎ、減衰比ζとしたときの伝達関数Ｇ（ｓ）＝（Ｋωｎ2 ／ｓ2 ＋２ζωｎｓ＋ωｎ2 ）・（１−Ｔｓ／１＋Ｔｓ）をモデルとし、前記固有角周波数と前記減衰比を正規分布乱数により変えて前記制御量のシミュレーションを行うステップと、
前記制御量と前記操作量の平均値に基づいてゲインを算出するステップと、
前記モデルにて制御を行ったときの波形を予測波形として取得するステップと、
前記実測波形と前記予測波形との波形の誤差の絶対値を加算した評価値が最も小さくなる前記固有角周波数と前記減衰比の組み合わせを算出するステップと、
比例帯Ｐ、積分時間Ｉ、微分時間Ｄを正規分布乱数により変えて前記制御量のシミュレーションを行ったときの前記制御量が前記目標値以上かつ前記制御量と前記目標値との差の最大値が最も小さくなる比例帯Ｐ、積分時間Ｉ、微分時間Ｄの組み合わせを算出するステップと、
を含むことを特徴とする制御方法。
【請求項４】
前記制御対象の特性に応じて初期時のＰＩＤ定数を設定するステップを含むことを特徴とする請求項３に記載の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば電気炉などの制御対象の温度などの物理量を制御する制御装置および制御方法に関するものである。
【背景技術】
【０００２】
例えば電気炉では、炉内温度と目標値を比較し、その偏差に応じて演算を行ってヒータなどの加熱手段を制御する制御装置として調節計が一般的に知られている。
【０００３】
ところで、ディジタル式の調節計では、制御対象に応じたＰＩＤ定数の自動算出と設定を行うオートチューニング機能を搭載している。このオートチューニング機能では、調節計を一時的にオンオフ調節計として用い、オン／オフ動作により発生するハンチングの周期と振幅の値からＰＩＤ定数を算出するリミットサイクル法を採用している。
【０００４】
なお、本件出願人は、上述したリミットサイクル法として、下記特許文献１に開示されるように、設定値付近でのサイクリング（リミットサイクル）波形からＰＩＤ定数を求め、操作量を決定する制御装置を提案している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平７−１３６０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上述したリミットサイクル法では、ＰＩＤ定数が適切に決定されていない場合、制御量が目標値を超過するオーバシュートという現象が発生する。このオーバシュートが発生すると、制御対象（例えば、製造品）に過大なストレスをかけてしまい、歩留りが悪化してしまうという問題があり、オーバシュートの抑制が望まれていた。
【０００７】
そこで、本発明は上記問題点に鑑みてなされたものであって、オーバシュートを抑制することができる制御装置および制御方法を提供することを目的としている。
【課題を解決するための手段】
【０００８】
上記目的を達成するため、本発明の請求項１に記載された制御装置は、制御対象の制御量を変化させるための操作量を算出して前記制御量を目標値に一致させる制御装置であって、
前記制御量が前記目標値よりも小さい中途目標値、前記目標値の順に到達するように２回のステップ応答の制御を行うステップ応答制御手段と、
前記２回のステップ応答の制御によって得られる前記目標値と前記中途目標値それぞれと前記制御量との差を示す偏差と該偏差の変化量が基準以下のときに定常状態と判定する状態判定手段と、
前記２回のステップ応答の制御で前記定常状態と判定されたときの前記制御量と前記操作量の平均値を算出する平均値算出手段と、
前記制御量が前記中途目標値に到達するようにステップ応答の制御をして前記定常状態と判定されたときに該定常状態から前記制御量が所定量変化するのに要するむだ時間を計測するむだ時間計測手段と、
前記２回のステップ応答の制御で前記定常状態と判定されたときの波形を実測波形として取得する実測波形取得手段と、
ラプラス演算子ｓ、ゲインＫ、むだ時間Ｔ、固有角周波数ωｎ、減衰比ζとしたときの伝達関数Ｇ（ｓ）＝（Ｋωｎ2 ／ｓ2 ＋２ζωｎｓ＋ωｎ2 ）・（１−Ｔｓ／１＋Ｔｓ）をモデルとし、前記固有角周波数と前記減衰比を正規分布乱数により変えて前記制御量のシミュレーションを行うモデル制御手段と、
前記平均値算出手段にて算出した前記制御量と前記操作量の平均値に基づいてゲインを算出するゲイン算出手段と、
前記モデルにて制御を行ったときの波形を予測波形として取得する予測波形取得手段と、
前記実測波形と前記予測波形との波形の誤差の絶対値を加算した評価値が最も小さくなる前記固有角周波数と前記減衰比の組み合わせを算出するパラメータ算出手段と、
比例帯Ｐ、積分時間Ｉ、微分時間Ｄを正規分布乱数により変えて前記制御量のシミュレーションを行ったときの前記制御量が前記目標値以上かつ前記制御量と前記目標値との差の最大値が最も小さくなる比例帯Ｐ、積分時間Ｉ、微分時間Ｄの組み合わせを算出するＰＩＤ定数算出手段と、
を具備することを特徴とする。
【０００９】
請求項２に記載された制御装置は、請求項１の制御装置において、
前記制御対象の特性に応じて初期時のＰＩＤ定数を設定することを特徴とする。
【００１０】
請求項３に記載された制御方法は、制御対象の制御量を変化させるための操作量を算出して前記制御量を目標値に一致させる制御方法であって、
前記制御量が前記目標値よりも小さい中途目標値、前記目標値の順に到達するように２回のステップ応答の制御を行うステップと、
前記２回のステップ応答の制御によって得られる前記目標値と前記中途目標値それぞれと前記制御量との差を示す偏差と該偏差の変化量が基準以下のときに定常状態と判定するステップと、
前記２回のステップ応答の制御で前記定常状態と判定されたときの前記制御量と前記操作量の平均値を算出するステップと、
前記制御量が前記中途目標値に到達するようにステップ応答の制御をして前記定常状態と判定されたときに該定常状態から前記制御量が所定量変化するのに要するむだ時間を計測するステップと、
前記２回のステップ応答の制御で前記定常状態と判定されたときの波形を実測波形として取得するステップと、
ラプラス演算子ｓ、ゲインＫ、むだ時間Ｔ、固有角周波数ωｎ、減衰比ζとしたときの伝達関数Ｇ（ｓ）＝（Ｋωｎ2 ／ｓ2 ＋２ζωｎｓ＋ωｎ2 ）・（１−Ｔｓ／１＋Ｔｓ）をモデルとし、前記固有角周波数と前記減衰比を正規分布乱数により変えて前記制御量のシミュレーションを行うステップと、
前記制御量と前記操作量の平均値に基づいてゲインを算出するステップと、
前記モデルにて制御を行ったときの波形を予測波形として取得するステップと、
前記実測波形と前記予測波形との波形の誤差の絶対値を加算した評価値が最も小さくなる前記固有角周波数と前記減衰比の組み合わせを算出するステップと、
比例帯Ｐ、積分時間Ｉ、微分時間Ｄを正規分布乱数により変えて前記制御量のシミュレーションを行ったときの前記制御量が前記目標値以上かつ前記制御量と前記目標値との差の最大値が最も小さくなる比例帯Ｐ、積分時間Ｉ、微分時間Ｄの組み合わせを算出するステップと、
を含むことを特徴とする。
【００１１】
請求項４に記載された制御方法は、請求項３の制御方法において、
前記制御対象の特性に応じて初期時のＰＩＤ定数を設定するステップを含むことを特徴とする。
【発明の効果】
【００１２】
本発明によれば、２回のステップ応答と伝達関数Ｇ（ｓ）のモデルによる制御量のシミュレーションによりＰＩＤ定数を適切に決定して設定することができ、制御量のオーバシュートを低減することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明に係る制御装置の概略構成を示すブロック図である。
【図２】本発明に係る制御装置および制御方法によるチューニング時のタイミングチャートの一例を示す図である。
【発明を実施するための形態】
【００１４】
以下、本発明を実施するための形態について、添付した図面を参照しながら詳細に説明する。
【００１５】
図１に示すように、本実施の形態の制御装置１は、例えば電気炉、流量系や化学反応系などの制御対象の温度、流量、湿度、圧力などの物理量を制御量とし、この制御量を操作量により変化させて制御量を目標値に一致させる調節計であり、ＰＩＤ定数を決定するための構成として操作表示部２、記憶部３、制御部４を備えて概略構成される。
【００１６】
ここでは、図１に示すように、例えばヒータなどの加熱手段１１や温度センサ１２が炉内に設置された電気炉１３を制御対象とし、温度センサ１２が検出する炉内温度（制御量ＰＶ）が目標温度（目標値ＳＶ）に一致するように操作量ＭＶでアクチュエータ（不図示）を駆動して加熱手段１１を制御する場合を例にとって説明する。
【００１７】
操作表示部２は、装置前面に設けられる各種キー、液晶やＬＥＤなどの表示器を備えて構成される。操作表示部２は、オートチューニングの開始の指示、オートチューニング開始時のＰＩＤ定数、中途目標値ＳＶ１、目標値ＳＶ２の設定を行う。また、操作表示部２は、設定画面の表示の他、目標値や制御量の表示、プログラムパターンなどの各種データ、バーグラフ表示やトレンド表示などの各種モニタ表示を行う。
【００１８】
なお、中途目標値ＳＶ１は、目標値ＳＶ２に対して例えば負荷率１０〜９０％で設定される。例えば目標値ＳＶ２の半分である負荷率５０％を中途目標値ＳＶ１として設定する。
【００１９】
記憶部３は、操作表示部２の操作により設定された設定値（ＰＩＤ定数、中途目標値ＳＶ１、目標値ＳＶ２）などを含め、電気炉１３の炉内温度の制御に関わる各種情報を記憶する。
【００２０】
制御部４は、操作表示部２の操作により設定された設定値、記憶部３に記憶された各種情報、電気炉１３の温度センサ１２が検出する炉内温度に基づいて制御装置１を統括制御するものである。制御部４は、ＰＩＤ定数を決定するための構成として、ステップ応答制御手段４ａ、状態判定手段４ｂ、平均値算出手段４ｃ、むだ時間計測手段４ｄ、実測波形取得手段４ｅ、シミュレーションモデル制御手段４ｆ、ゲイン算出手段４ｇ、予測波形取得手段４ｈ、評価値算出手段４ｉ、パラメータ算出手段４ｊ、ＰＩＤ定数算出手段４ｋを備える。
【００２１】
ステップ応答制御手段４ａは、操作表示部２にて設定された中途目標値ＳＶ１によるステップ応答１と目標値ＳＶ２によるステップ応答２の２回の制御を行う。
【００２２】
具体的に、ステップ応答制御手段４ａは、ステップ応答１として、操作表示部２にて設定された中途目標値ＳＶ１に制御量ＰＶが到達するように操作量ＭＶによりアクチュエータを駆動して加熱手段１１を制御（例えばＰＩＤ制御）する。
【００２３】
また、ステップ応答制御手段４ｂは、ステップ応答２として、操作表示部２にて設定された目標値ＳＶ２に制御量ＰＶが到達するように操作量ＭＶによりアクチュエータを駆動して加熱手段１１を制御（例えばＰＩＤ制御）する。
【００２４】
状態判定手段４ｂは、ステップ応答制御手段４ａによるステップ応答１，２の２回の制御で得られる中途目標値ＳＶ１と目標値ＳＶ２それぞれと制御量ＰＶとの差を示す偏差Ｅと偏差Ｅの変化量ΔＥに基づいて定常状態（制御量ＰＶが時間とともに変化しない状態）か否かを判定する。
【００２５】
具体的に、状態判定手段４ｂは、ステップ応答１において、中途目標値ＳＶ１と制御量ＰＶとの差を示す偏差Ｅと偏差Ｅの変化量ΔＥが基準以下、すなわち、偏差Ｅと偏差Ｅの変化量ΔＥがほぼ０に近ければ定常状態と判定する。
【００２６】
また、状態判定手段４ｂは、ステップ応答２において、目標値ＳＶ２と制御量ＰＶとの差を示す偏差Ｅと偏差Ｅの変化量ΔＥが基準以下、すなわち、偏差Ｅと偏差Ｅの変化量ΔＥがほぼ０に近ければ定常状態と判定する。
【００２７】
平均値算出手段４ｃは、ステップ応答１，２による制御で定常状態と判定されたときの制御量ＰＶと操作量ＭＶそれぞれの平均値を算出する。
【００２８】
具体的に、平均値算出手段４ｃは、ステップ応答１で定常状態と判定すると、定常状態における制御量ＰＶと操作量ＭＶそれぞれの平均値ＰＶ１、ＭＶ１を算出する。
【００２９】
また、平均値算出手段４ｃは、ステップ応答２で定常状態と判定すると、定常状態における制御量ＰＶと操作量ＭＶそれぞれの平均値ＰＶ２、ＭＶ２を算出する。
【００３０】
むだ時間計測手段４ｄは、ステップ応答１で定常状態と判定すると、この定常状態から制御量ＰＶが所定量変化するのに要するむだ時間Ｔを計測する。
【００３１】
具体的に、むだ時間計測手段４ｄは、定常状態から制御量ＰＶが温度センサ１２のスケール（検出し得る温度範囲）の例えば０．１％変化するまでの時間をむだ時間Ｔとして計測する。
【００３２】
実測波形取得手段４ｅは、ステップ応答１，２による制御で定常状態と判定されたときの波形を実測波形Ｗ１として取得する。
【００３３】
モデル制御手段４ｆは、制御対象１３との間の入出力をオフし、ラプラス演算子ｓ、ゲインＫ、むだ時間Ｔ、固有角周波数ωｎ、減衰比ζとしたときの伝達関数Ｇ（ｓ）＝（Ｋωｎ2 ／ｓ2 ＋２ζωｎｓ＋ωｎ2 ）・（１−Ｔｓ／１＋Ｔｓ）をモデルとして、固有角周波数ωｎと減衰比ζを正規分布乱数により変え、上述したステップ応答２と同条件で制御量ＰＶのシミュレーションを行う。
【００３４】
また、モデル制御手段４ｆは、ＰＶｓ＝Ｇ（ｓ）・Ｃｐｉｄを制御系として、比例帯Ｐ、積分時間Ｉ、微分時間Ｄを正規分布乱数により変えて制御量ＰＶのシミュレーションを行う。なお、ＰＶｓは制御量ＰＶのシミュレーション結果、ＣｐｉｄはＰＩＤ制御器である。
【００３５】
ゲイン算出手段４ｇは、平均値算出手段４ｃにて算出した制御量ＰＶと操作量ＭＶの平均値ＰＶ１，ＭＶ１，ＰＶ２，ＭＶ２に基づき、制御量ＰＶの平均値ＰＶ１，ＰＶ２と操作量ＭＶの平均値ＭＶ１，ＭＶ２の傾きから伝達関数Ｇ（ｓ）のモデルにおけるゲインＫを算出する。
【００３６】
具体的に、ゲイン算出手段４ｇは、ゲインＫ＝（ＰＶ２−ＰＶ１）／（ＭＶ２−ＭＶ１）の式に対し、平均値算出手段４ｃにて算出した制御量ＰＶと操作量ＭＶの平均値ＰＶ１，ＭＶ１，ＰＶ２，ＭＶ２を代入して伝達関数Ｇ（ｓ）のモデルにおけるゲインＫを算出する。
【００３７】
予測波形取得手段４ｈは、伝達関数Ｇ（ｓ）のモデルにて制御を行ったときの波形を予測波形Ｗ２として取得する。
【００３８】
評価値算出手段４ｉは、実測波形取得手段４ｅにて取得した実測波形Ｗ１と予測波形取得手段４ｈにて取得した予測波形Ｗ２との波形の誤差の絶対値を加算した値を評価値Ｖａｌ（＝Σ｜Ｗ２−Ｗ１｜）として算出する。
【００３９】
パラメータ算出手段４ｊは、評価値算出手段４ｉにて算出した評価値Ｖａｌが最も小さくなる固有角周波数ωｎと減衰比ζの組み合わせを算出する。
【００４０】
ＰＩＤ定数算出手段４ｋは、比例帯Ｐ、積分時間Ｉ、微分時間Ｄを正規分布乱数により変えてモデル制御手段４ｆにより制御量ＰＶのシミュレーションを行ったときの制御量ＰＶが目標値ＳＶ２以上かつ制御量ＰＶと目標値ＳＶ２との差の最大値が最も小さくなる比例帯Ｐ、積分時間Ｉ、微分時間Ｄの組み合わせをＰＩＤ定数として算出する。
【００４１】
次に、上記のように構成される制御装置１によるオートチューニング時の制御方法について図２を参照しながら説明する。
【００４２】
まず、操作表示部２の操作によりＰＩＤ定数を任意に設定するとともに、中途目標値ＳＶ１、目標値ＳＶ２を設定する。ここでは、目標値ＳＶ２の負荷率５０％が中途目標値ＳＶ１に設定されているものとする。なお、初期時のＰＩＤ定数は、制御対象１３の特性に応じて設定するのが好ましい。
【００４３】
上記設定を終え、操作表示部２の操作によりオートチューニングを開始すると、制御部４のステップ応答制御手段４ａは、中途目標値ＳＶ１に制御量ＰＶが到達するように操作量ＭＶによりアクチュエータを駆動して加熱手段１１をＰＩＤ制御する（図２のステップ応答１）。
【００４４】
続いて、ステップ応答制御手段４ｂは、目標値ＳＶ２に制御量ＰＶが到達するように操作量ＭＶによりアクチュエータを駆動して加熱手段１１をＰＩＤ制御する（図２のステップ応答２）。
【００４５】
そして、状態判定手段４ｂは、図２のステップ応答１において、中途目標値ＳＶ１と制御量ＰＶとの差を示す偏差Ｅと偏差Ｅの変化量ΔＥが基準以下であれば定常状態（図２の点線で囲むＳ１，Ｓ２の状態）と判定する。
【００４６】
また、状態判定手段４ｂは、図２のステップ応答２において、目標値ＳＶ２と制御量ＰＶとの差を示す偏差Ｅと偏差Ｅの変化量ΔＥが基準以下であれば定常状態（図２の点線で囲むＳ３，Ｓ４の状態）と判定する。
【００４７】
そして、平均値算出手段４ｃは、図２のステップ応答１で定常状態（図２の点線で囲むＳ１，Ｓ２の状態）と判定すると、この定常状態における制御量ＰＶと操作量ＭＶそれぞれの平均値ＰＶ１、ＭＶ１を算出する。
【００４８】
また、平均値算出手段４ｃは、図２のステップ応答２で定常状態（図２の点線で囲むＳ３，Ｓ４の状態）と判定すると、この定常状態における制御量ＰＶと操作量ＭＶそれぞれの平均値ＰＶ２、ＭＶ２を算出する。
【００４９】
さらに、むだ時間計測手段４ｄは、図２のステップ応答１で定常状態と判定すると、この定常状態から制御量ＰＶが所定量変化するのに要するむだ時間Ｔを計測する。
【００５０】
そして、実測波形取得手段４ｅは、ステップ応答１，２による制御で定常状態と判定されたときの波形を実測波形Ｗ１として取得する。
【００５１】
次に、モデル制御手段４ｆは、制御対象１３との間の入出力をオフした状態で、ラプラス演算子ｓ、ゲインＫ、むだ時間Ｔ、固有角周波数ωｎ、減衰比ζとしたときの伝達関数Ｇ（ｓ）＝（Ｋωｎ2 ／ｓ2 ＋２ζωｎｓ＋ωｎ2 ）・（１−Ｔｓ／１＋Ｔｓ）をモデルとして、固有角周波数ωｎと減衰比ζを正規分布乱数により変えて上述したステップ応答２と同条件で制御量ＰＶのシミュレーションを行う。
【００５２】
そして、予測波形取得手段４ｈは、伝達関数Ｇ（ｓ）のモデルにて制御を行ったときの波形を予測波形Ｗ２として取得する。
【００５３】
ここで、ゲイン算出手段４ｇは、ゲインＫ＝（ＰＶ２−ＰＶ１）／（ＭＶ２−ＭＶ１）の式に対し、平均値算出手段４ｃにて算出した制御量ＰＶと操作量ＭＶの平均値ＰＶ１，ＭＶ１，ＰＶ２，ＭＶ２を代入して伝達関数Ｇ（ｓ）のモデルにおけるゲインＫを算出する。
【００５４】
また、評価値算出手段４ｉは、実測波形取得手段４ｅにて取得した実測波形Ｗ１と予測波形取得手段４ｈにて取得した予測波形Ｗ２との波形の誤差の絶対値を加算した値を評価値Ｖａｌ（＝Σ｜Ｗ２−Ｗ１｜）として算出する。
【００５５】
そして、パラメータ算出手段４ｊは、評価値算出手段４ｉにて算出した評価値Ｖａｌが最も小さくなる固有角周波数ωｎと減衰比ζの組み合わせを算出する。
【００５６】
次に、ＰＩＤ定数算出手段４ｋは、比例帯Ｐ、積分時間Ｉ、微分時間Ｄを正規分布乱数により変えてモデル制御手段４ｆの制御による制御量ＰＶのシミュレーションを行う。そして、この制御量ＰＶのシミュレーションを行ったときの制御量ＰＶが目標値ＳＶ２以上かつ制御量ＰＶと目標値ＳＶ２との差の最大値が最も小さくなる比例帯Ｐ、積分時間Ｉ、微分時間Ｄの組み合わせをＰＩＤ定数として算出する。これにより、最適な比例帯Ｐ、積分時間Ｉ、微分時間Ｄを求めてＰＩＤ定数を設定することができる。
【００５７】
ところで、上述した実施の形態では、制御対象として電気炉１３の炉内温度を目標値ＳＶに制御する制御装置を例にとって説明したが、これに限定されるものではない。例えば流量系や化学反応系などにおける流量、湿度、圧力などの物理量を目標値に制御する制御装置として用いることもできる。
【００５８】
このように、本実施の形態によれば、２回のステップ応答と伝達関数Ｇ（ｓ）のモデルによる制御量のシミュレーションによりＰＩＤ定数を適切に決定して設定することができ、制御量のオーバシュートを低減することができる。
【００５９】
以上、本発明に係る制御装置および制御方法の最良の形態について説明したが、この形態による記述及び図面により本発明が限定されることはない。すなわち、この形態に基づいて当業者等によりなされる他の形態、実施例及び運用技術などはすべて本発明の範疇に含まれることは勿論である。
【符号の説明】
【００６０】
１  制御装置
２  操作表示部
３  記憶部
４  制御部
１１  加熱手段
１２  温度センサ
１３  電気炉（制御対象）
【図１】
2020160550.tif 000003
【図２】
2020160550.tif 000004