(57)【特許請求の範囲】
【請求項１】
制御対象の伝達関数を含むフィ−ドバック制御系に自励振動を発生させて、制御系の限界ゲインを同定する限界ゲイン同定方法であって、
前記フィードバック制御系は、可変ゲイン要素と飽和要素とを有しており、
前記可変ゲイン要素の入力で自励振動の振幅および振幅の変化率を計測しながら、下記の式によって前記可変ゲイン要素のゲインｋを十分に大きな値から漸減させ、
自励振動の計測から、自励振動が一定振幅に収束したときのゲインｋを限界ゲインとする
ことを特徴とするフィ−ドバック制御系の限界ゲイン同定方法。
【数１】
0003774376.tif 000036
ここで、
ｅ（ｔ）：可変ゲイン要素の入力
ｋ：可変ゲイン要素の値
バ−ｅ：ｅ（ｔ）の振幅
ｕ０：飽和要素の飽和値
β１＞０，β２≧０：定数
【請求項２】
請求項１記載の限界ゲイン同定方法において、前記フィードバック制御系は、さらに、位相を設定できる一次遅れ要素を有するとともに、自励振動の周波数ωおよび周波数ωの変化率を、前記可変ゲイン要素の入力で計測し、
前記可変ゲイン要素のゲインｋを十分に大きな値から漸減するとともに、前記一次遅れ要素に対して、下記の式によって十分小さな時定数Ｔ（Ｔ＞０）から適応的に調整し、
自励振動の計測により、自励振動が一定振幅に収束したときのゲインｋ，時定数Ｔ，周波数ωで、所定位相における前記フィードバック制御系における前記伝達関数を同定する
ことを特徴とする伝達関数同定方法。
【数２】
0003774376.tif 000037
ここで、
Ｔ：一次遅れ要素の時定数
ω：周波数
θｃ：所定位相値
α＞０：定数
【請求項３】
制御対象の伝達関数を含むフィ−ドバック制御系に自励振動を発生させて、制御系の限界ゲインを同定する限界ゲイン同定装置であって、
前記フィードバック制御系は、可変ゲイン要素と飽和要素とを有し、
さらに、前記ゲイン要素の入力で自励振動の振幅および振幅の変化率を計測する計測手段と、
下記の式によって前記可変ゲイン要素のゲインｋを十分大きな値から漸減させる調整手段とを有し、
前記自励振動が一定振幅に収束したことを前記計測手段で計測したとき、そのゲインｋを前記フィードバック制御系の限界ゲインとすることを特徴とするフィ−ドバック制御系の限界ゲイン同定装置。
【数３】
0003774376.tif 000038
ここで、
ｅ（ｔ）：可変ゲイン要素の入力
ｋ：可変ゲイン要素のゲイン
バ−ｅ：ｅ（ｔ）の振幅
ｕ０：飽和要素の飽和値
β１＞０，β２≧０：定数
【請求項４】
請求項３記載の限界ゲイン同定装置において、前記フィードバック制御系は、さらに位相を設定できる一次遅れ要素を有しており、
前記可変ゲイン要素の入力で計測する前記計測手段は、自励振動の周波数ω及び周波数ωの変化率も計測し、
前記調整手段は、前記一次遅れ要素を所定位相に設定することもでき、
前記調整手段は、前記可変ゲイン要素のゲインとともに、前記一次遅れ要素に対して位相を、下記の式によって十分小さな時定数Ｔ（Ｔ＞０）から所定位相に適応的に調整することにより、
前記計測手段により、自励振動が一定振幅に収束したことを計測したときのゲインｋ，時定数Ｔ，周波数ωで、所定位相における前記フィードバック制御系における前記伝達関数を同定することを特徴とする伝達関数同定装置。
【数４】
0003774376.tif 000039
ここで、
Ｔ：一次遅れ要素の時定数
ω：周波数
θｃ：所定位相値
α＞０：定数
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、制御系のゲインの同定方法および同定装置に関する。
【０００２】
【従来の技術】
ＰＩＤ制御器を用いたフィ−ドバック制御系の制御パラメ−タを適性値に設定するためにジ−グラ・ニコルスの限界感度法が知られている。この方法は、実稼動状態にある系で比例ゲインを少しずつ増加していくと系の振動が始まるので、このときのゲインと位相を元に制御系の制御パラメ−タを調整する。閉ル−プ系を安定限界にする比例ゲインが限界ゲインであり、その時の振動周期が限界周期である。
ジ−グラ・ニコルスの限界感度法を用いてパラメ−タの設定を自動的に行なう技術としては、特開平９−３４５０３号公報「ＰＩＤコントロ−ラの調整法」（明電舎）などがある。しかし、これらの方法は限界ゲインよりもかなり大きなゲインにならないと、振動が観測されないので、求められる限界ゲインの精度が悪い。また、特開平９−１０６３０３号公報「制御系のゲインの自動決定法」（ファナック）では、速度偏差によるゲイン調整系により自励振動を起こし、非線形領域での近似ゲインを求めているので、正確な限界ゲインが同定できない。
米国特許明細書第４，５４９，１２３号には、ル−プ内にリレ−要素を挿入して自励振動を起こし、記述関数法の考え方を導入して、精度のより高い限界ゲインを求める方法が示されている。しかし、この方法は、制御対象が１入力１出力系に限られる。また、外乱やノイズによる限界ゲイン同定精度の低下が見られる点で課題がある。
【０００３】
発明者らは、制御ル−プ内に飽和要素を導入し、その入出力信号の偏差を用いて適応的にル−プゲインを調整することで、精度良く限界ゲインを求めることが理論的に可能であることを提案した（「限界ゲインの高精度適応同定法の提案」計測自動制御学会第１８回適応制御シンポジウム１９９８年１月）。この方法は、記述関数による近似が不要である、振幅を測定せずに限界ゲインを陽に与えることができる、多変数系に適用できる、などの利点がある。しかし、これを実現する具体的な方法はいまだ提示されていないし、求められるゲイン値は細かな振動を示し、限界ゲインの一定値に収束しないという課題もある。また、限界ゲイン・限界周期が同じでも、特性が異なる制御系があり、このときは、限界ゲインを与える位相（−１８０°）以外の他の位相におけるゲインを知ることが、系の調整上で重要であるが、上記の文献にはその方法が示されていない。
【０００４】
【発明が解決しようとする課題】
本発明は、フィ−ドバック制御系における限界ゲインや限界周期を、プロセスの運転状態で精度良く同定する方法と装置を提供することを目的とする。また、位相交点以外の位相における伝達関数点を、同様に同定する方法・装置を提供することを目的とする。
【０００５】
【課題を解決するための手段】
上記目的を達成するために、本発明は、制御対象の伝達関数を含むフィ−ドバック制御系自励振動を発生させて、制御系の限界ゲインを同定する限界ゲイン同定方法であって、前記フィードバック制御系は、可変ゲイン要素と飽和要素とを有しており、前記可変ゲイン要素の入力で自励振動の振幅および振幅の変化率を計測しながら、下記の式によって前記可変ゲイン要素のゲインｋを十分に大きな値から漸減させ、自励振動の計測から、自励振動が一定振幅に収束したときのゲインｋを限界ゲインとするフィ−ドバック制御系の限界ゲインの同定方法である。
【数５】
0003774376.tif 000002
ここで、ｅ（ｔ）：可変ゲイン要素の入力、ｋ：可変ゲイン要素の値、バ−ｅ：ｅ（ｔ）の振幅、ｕ０：飽和要素の飽和値、β１＞０、β２≧０：定数
上述の限界ゲイン同定方法において、前記フィードバック制御系は、さらに位相を設定できる一次遅れ要素を有しており、また、自励振動の周波数ω及び周波数ωの変化率を、前記可変ゲイン要素の入力で計測し、計測された振幅と周波数から、前記可変ゲイン要素のゲインを漸減するとともに、前記一次遅れ要素に対して位相を下記の式によって所定位相に十分小さな時定数Ｔ（Ｔ＞０）から適応的に調整し、自励振動の計測により、自励振動が一定振幅に収束したときのゲインｋ，時定数Ｔ，周波数ωで、所定の位相における前記フィードバック制御系における前記伝達関数を同定することもできる。
【数６】
0003774376.tif 000003
ここで、Ｔ：一次遅れ要素の時定数、ω：周波数、θｃ：所定位相値、α＞０：定数
また、本発明は、制御対象の伝達関数を含むフィ−ドバック制御系に自励振動を発生させて、制御系の限界ゲインを同定する限界ゲイン同定装置であって、前記フィードバック制御系は、可変ゲイン要素と飽和要素とを有し、さらに、前記ゲイン要素の入力で自励振動の振幅および振幅の変化率を計測する計測手段と、下記の式によって前記可変ゲイン要素のゲインｋを十分大きな値から漸減させる調整手段とを有し、自励振動が一定振幅に収束したことを前記計測手段で計測したとき、そのゲインｋを前記フィードバック制御系の限界ゲインとするフィ−ドバック制御系の限界ゲインの同定装置である。
【数７】
0003774376.tif 000004
ここで、ｅ（ｔ）：可変ゲイン要素の入力、ｋ：可変ゲイン要素のゲイン、バ−ｅ：ｅ（ｔ）の振幅、ｕ０：飽和要素の飽和値、β１＞０、β２≧０：定数
上述の限界ゲインの同定装置において、前記フィードバック制御系は、さらに位相を設定できる一次遅れ要素を有しており、前記可変ゲイン要素の入力で振幅を計測する計測手段は、自励振動の周波数ω及び周波数ωの変化率も計測し、前記調整手段は、前記一次遅れ要素を所定位相に設定することもでき、前記調整手段は、計測された振幅と周波数から、前記可変ゲイン要素のゲインとともに、前記一次遅れ要素に対して位相を、下記の式によって十分小さな時定数Ｔ（Ｔ＞０）から所定位相に適応的に調整し、自励振動の計測により、自励振動が一定振幅に収束したときのゲインｋ，時定数Ｔ，周波数ωで、所定の位相における前記フィードバック制御系における前記伝達関数を同定することもできる。
【数８】
0003774376.tif 000005
ここで、Ｔ：一次遅れ要素の時定数、ω：周波数、θｃ：所定位相値、α＞０：定数
【０００６】
【発明の実施の形態】
本発明の実施形態を、図面を参照して詳細に説明する。
本発明は、系に自励信号を発生させて、適応的にゲインを調整して、限界ゲインを求めている。以下に、この方法により限界ゲインが求まることをまず説明し、また、このときのゲインの調整の方法を説明する。
さて、図１に示すゲイン要素１０および関数要素２０で構成されたフィ−ドバック系を例に、限界ゲインの求め方を説明する。
図１において、関数要素２０のＧ（ｓ）は周波数特性を同定したいプラントの伝達関数であり、ゲイン要素１０のゲインｋはｋ∈Ｒである。図１に示したこのフィ−ドバック系は、０＜ｋ＜Ｋｃｒで安定であり、ｋ＞Ｋｃｒで不安定になるとする。この場合、ｋ＝Ｋｃｒにおいては自励振動が持続する。このとき、Ｋｃｒを限界ゲインと呼び、自励振動の周期を限界周期ωｃｒと呼ぶ。Ｋｃｒとωｃｒは、次式のように、図１のフィ−ドバック系の特性方程式を満たす。
【数９】
１＋ＫｃｒＧ（ｊωｃｒ）＝０                  （１）
この特定方程式により、Ｋｃｒとωｃｒを同定することで、周波数特性が
【数１０】
Ｇ（ｊωｃｒ）＝−１／Ｋｃｒ            （２）
により同定できる。
【０００７】
（Ｋｃｒとωｃｒの同定の考え方）
図２の非線形フィ−ドバック系を考える。図２で示したフィードバック系は、図１のフィードバック系に飽和要素３０を挿入した構成を有している。ここに飽和要素３０のφは、図３に示すような、飽和値ｕ０で飽和する関数出力を有している。
ゲイン要素の値（ゲイン）ｋが無限大に近づくとき、飽和要素の出力関数φ（ｋσ）の極限はリレ−要素の特性をもつので、この系は、十分大きなｋに対して、リレ−制御系と同様に振舞う。すなわち、ｋが十分大きいとき、不安定モ−ドの振幅が急速に増大し、しばらくして振幅一定な自励振動状態に落ち着く。このとき、ゲインｋが大きいので、この定常状態におけるゲイン要素１０の出力ｖの振幅ｖＥ（バーｖと示す場合もある）は、飽和値ｕ０に比べて十分大きい。そこで、ゲインｋを漸減させることにより、振幅ｖＥを飽和値ｕ０に近づける。すなわち、自励振動を維持しながら飽和状態から非飽和状態に近づけ、最終的に、振幅ｖＥ＝ｕ０で振動を保つようにする。このとき、非飽和状態での振動が生じているので、そのときのゲインｋと振動周波数ωが、限界ゲインＫｃｒと限界周期ωｃｒを与える。これを達成するために、ｖＥ＞ｕ０であればｋを減少させ、ｖＥ＜ｕ０であればｋを増加させる調整則：
【数１１】
0003774376.tif 000006
を用いる。ただし、β１＞０、β２≧０であり、初期値は十分に大きな値とする。右辺第２項は飽和値の変化速度を微調整に利用するために用いるので、β２はゼロでも良い。また、（３）式では、ゲイン要素１０の入力ｅ（ｔ）の振幅ｅＥ（バーｅ（ｔ）とも記載）を推定し、ｖＥ＝ｋｅＥにより、ゲイン要素１０の出力ｖＥを求めることとし、その有用性は後述する。
【０００８】
＜ゲイン調整則の特性＞
（ｋ（ｔ）の過渡応答特性）
図４のグラフは、調整則を適用した場合のｅ（ｔ）の過渡応答例を示す。点線がゲイン要素１０の入力ｅ（ｔ）であり、実線がその振幅の周波数推定器（後述）による推定値ｅＥ（ｔ）である。このように、ｅＥ（ｔ）はゼロから急速に増加し一定値に落ち着くので、適当な正数ｅ１、ｅ２に対し、ｅ１≧ｅ（ｔ）≧ｅ２が満たされる。この性質から、式（３）を満たすゲインｋ（ｔ）は、β２＝０のとき、つぎの区間内に存在することが保証される。
【数１２】
ｋ２（ｔ）≧ｋ（ｔ）≧ｋ１（ｔ）                      （４）
ここに、
【数１３】
0003774376.tif 000007
このことの重要性を次に述べる。（３）式では、振幅ｅＥ（ｔ）を用いて、ゲイン要素の出力の振幅ｖＥ＝ｋｅＥ（ｔ）を求めたが、直接に、ゲイン要素の出力ｖ（ｔ）から周波数推定器を用いてｖＥを推定する方法：
【数１４】
0003774376.tif 000008
が考えられる。しかし、式（７）を用いた調整では、ｋ（ｔ）が負になることで非振動的な不安定現象が生じ、同定が失敗する場合があった。これに対し、（３）式を用いた調整では、（４）式から分かるようにｋ（ｔ）を正に保てるので、そのような欠点はない。このため、周波数推定器のパラメ−タの設定が、（３）式のほうが遥かに容易である。
【０００９】
（平衡状態の安定性）
上記の過渡特性に従い、ｋ（ｔ）を十分大きな初期値ｋ（０）から漸減させている。ｋ＝Ｋｃｒが同定できるには、（３）式により振幅ｖＥ＝ｕ０の自励振動が安定に保たれることが必要であり、以下ではその条件を求める。
図１のフィ−ドバック系の１＋ｋＧ（ｓ）＝０の主要極ｓはｋ＝Ｋｃｒのときｓ＝−ｊωｃｒであるので、ｋ＝Ｋｃｒの付近における主要極の近似式は
【数１５】
0003774376.tif 000009
で与えられる。これに対する振動モ−ドの振幅は
【数１６】
0003774376.tif 000010
で与えられるので、つぎのモデルを得る。
【数１７】
0003774376.tif 000011
ただし、ｋ＞ＫｃｒではｅＥが増加し、逆にｋ＜ＫｃｒではｅＥが減少するので、ａ＞０である。
このモデルと（３）式の調整則からなる系
【数１８】
0003774376.tif 000012
の平衡点は、ｋ＝Ｋｃｒ、ｅＥ＝ｅＥ０（＝ｕ０／Ｋｃｒ）で与えられる。そこで、新たに、変数をｘ１＝ｋ−Ｋｃｒ、ｘ２＝ｅＥ−ｅＥ０とおいて、平衡点まわりの線形近似式を求めると
【数１９】
0003774376.tif 000013
である。この特性方程式は
【数２０】
0003774376.tif 000014
であるので、β１、β２により任意に極配置でき、安定化のためには、β１＞０、β２＝０で十分である。以上より、この条件を満たすようにβ１、β２を選ぶことで、安定な自励振動が保たれ、ｋ＝Ｋｃｒが高精度に同定できる。
【００１０】
（自励振動の周期と振幅の同定法）
式（３）のアルゴリズムを実現するためには、自励振動の振幅の包絡線であるｅＥ（ｔ）を測定する必要がある。ここでは、適応ノッチ・フィルタを伴った周波数推定器を用いて、振動周期と共に、新たに、振幅の包絡線ｅＥ（ｔ）をオンラインで同定する方法を与える。
周波数推定器への入力信号が正弦波信号
【数２１】
0003774376.tif 000015
であるとする。このとき、周波数推定器が次式で与えられる。
【数２２】
0003774376.tif 000016
ここで、
【数２３】
0003774376.tif 000017
および、振動周波数ωは推定された周波数、ζ＞０はダンピング係数、γは可変パラメ−タである。また、ε＞０、Ｎ＞０、μ＞０、α２≧１である。
入力信号ｎ（ｔ）に対する定常応答は、次式で表される周期解となる。
【数２４】
0003774376.tif 000018
このとき、
【数２５】
0003774376.tif 000019
の関係が成り立つので、過渡時における振幅ｖ＞０と変化率ｄｖ／ｄｔの推定値が次式で与えられる。
【数２６】
0003774376.tif 000020
以上のように作用するので、系に自励振動を発生させて、適応的にゲインを調整して、限界ゲインを求めることができる。
【００１１】
＜任意の位相角における伝達関数の同定＞
次に、フィ−ドバック系に、ゲイン要素と一次遅れ要素（位相設定器と同じ意味である）を挿入し、限界ゲインを求める方法によって、任意の位相角における伝達関数を同定する方法について説明する。
図５のフィ−ドバック系を考える。図５は、図１のフィードバック系に一次遅れ要素４０を挿入した構成である。ここで、図１と同様に、Ｇ（ｓ）は周波数応答を求めたいプラントの伝達関数であり、ゲインｋはｋ∈Ｒである。図中の一次遅れ要素４０は、一次遅れ要素の時定数Ｔを用いて
【数２７】
0003774376.tif 000021
で表され、これによる位相遅れθは、一次遅れ要素の時定数Ｔにより、
【数２８】
0003774376.tif 000022
で与えられ、０＜θ＜π／２である。この系において、
【数２９】
0003774376.tif 000023
を満たすゲインｋ、一次遅れ要素の時定数Ｔ、ω＝ωｃが与えられたとする。このとき、
【数３０】
0003774376.tif 000024
であるから、
【数３１】
0003774376.tif 000025
である。したがって、Ｇ（ｊωｃ）は位相が−π＋θｃとなる周波数の伝達関数である。
そこで、θｃを与えて、
【数３２】
0003774376.tif 000026
を満たす一次遅れ要素の時定数Ｔを求める方法を検討する。
この一次遅れ要素の時定数Ｔに対し、系が安定限界にあれば式（２６）を満足するので、ω＝ωｃの自励振動を生じるように、上述した適応同定法によりゲインｋを調整する。これは、図６に示すように、θｃだけ回転させた座標系で見た限界ゲインを求めることと考えられる。
【００１２】
（調整則）
式（２９）を満たすＴを求める方法のための、一次遅れ要素の時定数Ｔの調整則を以下に説明する。まず、
【数３３】
0003774376.tif 000027
とおく。このとき、
【数３４】
0003774376.tif 000028
が成り立てば、ωＴとｔａｎθｃの偏差が減少する。この条件と
【数３５】
0003774376.tif 000029
より、
【数３６】
0003774376.tif 000030
を得る。したがって、一次遅れ要素の時定数Ｔの調整則として
【数３７】
0003774376.tif 000031
を用いる。
ところで、式（３１）より、
【数３８】
Ｖ（ｔ）＝Ｖ（０）ｅ−２αｔ                （35）
であるから
【数３９】
0003774376.tif 000032
ゆえに、初期偏差がｅ−αｔで減衰する。
ゲインｋの調整には、これまでの調整則を用いる。したがって、
【数４０】
0003774376.tif 000033
の両式を調整則とする。ここで、ｅＥとωは自励振動の振幅と周期であり、ｅＥ、ｄｅＥ／ｄｔ、ω、ｄω／ｄｔは周波数推定器により推定する。
初期値としてＴ＝Ｔ（０）を正で十分小さく、また、ｋ＝Ｋ（０）を十分大きく与える。この条件の下で実験し、ｋとＴが収束すれば、ωｃにおける伝達関数が
【数４１】
0003774376.tif 000034
で与えられる。以上のように作用するので、フィ−ドバック系に、図５に示すように、ゲイン要素１０と一次遅れ要素４０を挿入し、限界ゲインを求める方法によって、任意の位相角における伝達関数を同定することができる。
【００１３】
【実施例】
（数値シミュレ−ション）
本発明の限界ゲインを同定する上述のアルゴリズムを実施するために、図７に示すシミュレーション・システムを構成した。図７において、このシステムは、周波数推定器５０、上述の式（３）に示した調整を行う適応アルゴリズムを内蔵したゲイン調整器６０、可変ゲイン要素１０、飽和要素７０、対象制御系の一巡伝達関数器２０で構成されている。本発明を実施するために挿入される周波数推定器５０、調整器６０、可変ゲイン要素１０、飽和要素７０は、プロセスで一般的に使用されているデジタル制御システムの機能の一部として限界ゲインを同定するときに挿入されるものである。あるいは、別のデジタル制御装置として構成し、既存の制御系に限界ゲインを同定するときに挿入してもよい。
周波数推定器５０は、（１６）〜（１９）式で与えられるアルゴリズムが組み込まれており、（２２）式と（２３）式に従いｅＥ、ｄｅＥ／ｄｔの推定値を出力する（ただし、ｖがｅＥに対応する）。また、周波数推定器の機能により振動の周期ωを観測する。調整器６０は、（３）式で与えられる調整則を内臓して、可変ゲイン要素１０のゲインｋを調整する。飽和要素７０は、同定中に出力が異常とならないよう、安全のために挿入している。
限界ゲインの同定時は、自励振動の信号ｅ（ｔ）を周波数推定器５０に入力し、自励振動の振幅ｅＥ（ｔ）を推定する。この推定された振幅をもとに、調整器６０の調整則により可変ゲイン要素６０のゲインｋを適応的に変化させる。これにより、限界ゲインＫｃｒを同定する。このとき、限界周期ωｃｒも同時に同定される。
以下に、上記図７で示されたシステムに対して、具体的に一巡伝達関数Ｇ（ｓ）の限界ゲイン等を特定した場合を以下に説明する。
【００１４】
（数値例  １）
一巡伝達関数が（４０）式で表されるプラントを考える。
【数４２】
0003774376.tif 000035
パラメ−タは次のようにする。
調整則：  β１＝０．２、β２＝０、ｋ（０）＝１０、ｕ０＝１
適応ノッチフィルタ（ＡＮＦ）：ε＝０．６、ζ＝０．６、μ＝０．２、Ｎ＝１、α２＝２
そして、外乱として、大きさ０．００１のステップ入力を時刻０に加えたときの結果を図８〜図１０に示す。ナイキストの安定条件より計算される真値、ｋｃｒ＝０．１５０７とωｃｒ＝１．０３２に速やかに収束している。
これにより、本発明の方法によって、理論的に求められる限界ゲインが高精度に求められていることが示されている。
【００１５】
（第２の実施例）
本発明の第２の実施例を図１１に示す。
本発明の所定位相における伝達関数点を同定するアルゴリズムを実施するために、図１１のシステムを構成する。図１１に示すように、システムは、周波数推定器５２、調整器６２、可変ゲイン８０、一次遅れ要素９０、飽和要素７０、対象制御系の伝達関数２０で構成されている。本発明を実施するために挿入される周波数推定器５２、調整器６２、可変ゲイン８０、一次遅れ要素９０、飽和要素７０は、プロセスで一般的に使用されているデジタル制御システムの機能の一部として限界ゲインを同定するときに挿入されるものである。
周波数推定器５２は、振動の周期の推定値を出力する。また、（１６）〜（１９）式で与えられるアルゴリズムが組み込まれており、（２２）式と（２３）式に従い推定値ｅＥ、ｄｅＥ／ｄｔを出力し（ただし、ｖがｅＥに対応する）、（１８）式に従いｄω／ｄｔを出力する。調整器６２は、（３７）および（３８）式で与えられる調整則を内臓して、可変ゲイン８０、一次遅れ要素９０のゲインおよび一次遅れ要素部分（位相角）を調整する。また、飽和要素７０は、実施例１と同様に、同定中に出力が異常とならないよう、安全の目的で挿入したものである。
【００１６】
（数値例  ２）
次のプラントに対し、φｃ＝４０°となるときの伝達関数点を求める。
パラメ−タは次のようにする。
適応則：β１＝０．３、β２＝０、ｕ０＝１、α＝５．０
初期値：ｋ（０）＝１０、ω（０）＝１０、Ｔ（０）＝０．００１
外乱として、大きさ０．００１のステップ入力を時刻０に加えたときの結果を  図１２〜図１５に示す。
同様に、φｃ＝２０°、６０°、８０°においてもシミュレ−ションを行い、定常状態におけるゲインｋ、遅れ時間Ｔ、周波数ωの値から計算される伝達関数点を図１６に示す。実線がナイキスト線図であり、＊で示した点が本実施の形態により推定された伝達関数点である。図１６から、それぞれの位相において、高精度に同定されているのが確認できる。
上述で説明したように、限界ゲインの適応同定法を用いて指定した位相のゲインと周波数を求める方法をシミュレ−ションした。これにより、位相交点以外の周波数においても、ゲインと周波数を高精度に同定することができることが分かる。
（応用ができる技術・製品）
加熱・蒸留・合成・などの化学プロセスなどフィ−ドバック制御系を用いているプロセス一般に適用可能。
【００１７】
【発明の効果】
（１）運転中の実機において制御系の精度のよい限界ゲイン・限界周期の同定ができる。また、これは、１変数１出力系のみならず、多変数系多出力系においても適用できる。これにより、制御系のパラメ−タの最適な値が設定でき、性能のよいプロセスの制御が可能となる。
（２）位相が９０度〜１８０度遅れる時の系の伝達関数点を同定することができる。これによって、対象制御系のゲインを適切に調整することができる。
【図面の簡単な説明】
【図１】ゲイン要素を持つフィ−ドバック系の定義を説明する図である。
【図２】非線型フィ−ドバック系を説明する図である。
【図３】飽和要素の特性を示す図である。
【図４】振幅の包絡線の変化を示す図である。
【図５】位相遅れ要素を持つフィ−ドバック制御系を説明する図である。
【図６】位相回転の意味を説明する図である。
【図７】限界ゲインの同定装置を示す図である。
【図８】出力と包絡線の変化を示す図である。
【図９】ゲインの変化を示す図である。
【図１０】周期の変化を示す図である。
【図１１】所定位相における限界ゲインの同定装置を示す図である。
【図１２】ゲインの変化を示す図である。
【図１３】周期の変化を示す図である。
【図１４】出力と包絡線の変化を示す図である。
【図１５】遅れ時間の変化を示す図である。
【図１６】ナイキスト線図と同定点を示す図である。
【図１】
0003774376.tif 000040
【図２】
0003774376.tif 000041
【図３】
0003774376.tif 000042
【図４】
0003774376.tif 000043
【図５】
0003774376.tif 000044
【図６】
0003774376.tif 000045
【図７】
0003774376.tif 000046
【図８】
0003774376.tif 000047
【図９】
0003774376.tif 000048
【図１０】
0003774376.tif 000049
【図１１】
0003774376.tif 000050
【図１２】
0003774376.tif 000051
【図１３】
0003774376.tif 000052
【図１４】
0003774376.tif 000053
【図１５】
0003774376.tif 000054
【図１６】
0003774376.tif 000055