JPH0550002B2 - - Google Patents
Info
- Publication number
- JPH0550002B2 JPH0550002B2 JP62262031A JP26203187A JPH0550002B2 JP H0550002 B2 JPH0550002 B2 JP H0550002B2 JP 62262031 A JP62262031 A JP 62262031A JP 26203187 A JP26203187 A JP 26203187A JP H0550002 B2 JPH0550002 B2 JP H0550002B2
- Authority
- JP
- Japan
- Prior art keywords
- membership function
- value
- control
- inference
- composite
- 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.)
- Expired - Lifetime
Links
Landscapes
- Feedback Control In General (AREA)
- Devices For Executing Special Programs (AREA)
Description
【発明の詳細な説明】
〔産業上の利用分野〕
この発明は、各種工業プロセスを監視して、当
該工業プロセスに対してふさわしいパラメータの
値を推論する、再帰型のフアジイ推論装置に関す
るものである。
該工業プロセスに対してふさわしいパラメータの
値を推論する、再帰型のフアジイ推論装置に関す
るものである。
第5図は例えば「省力と自動化」の1986年11月
号第61〜66頁の『フアジイ・システム理論とフア
ジイ制御』(浅居 喜代治)に示された従来のフ
アジイ推論装置の動作原理を示す説明図であり、
図において、1及び2は推論のルール、3及び4
は当該フアジイ推論装置に入力される特徴量であ
り、ここではそれぞれ制御系の制御偏差eと制御
偏差の変化率Δeである。また、5及び6は前記
ルール1の前件部のメンバーシツプ関数、7は当
該ルール1の後件部のメンバーシツプ関数、8及
び9は前記ルール2の前件部のメンバーシツプ関
数、10は当該ルール2の後件部のメンバーシツ
プ関数である。さらに、11は前記メンバーシツ
プ関数7及び10を合成したメンバーシツプ関
数、12はこのメンバーシツプ関数11から重心
をとつて得られる推論値であり、この例では操作
量の変化率Δuとして当該フアジイ推論装置より
出力される。
号第61〜66頁の『フアジイ・システム理論とフア
ジイ制御』(浅居 喜代治)に示された従来のフ
アジイ推論装置の動作原理を示す説明図であり、
図において、1及び2は推論のルール、3及び4
は当該フアジイ推論装置に入力される特徴量であ
り、ここではそれぞれ制御系の制御偏差eと制御
偏差の変化率Δeである。また、5及び6は前記
ルール1の前件部のメンバーシツプ関数、7は当
該ルール1の後件部のメンバーシツプ関数、8及
び9は前記ルール2の前件部のメンバーシツプ関
数、10は当該ルール2の後件部のメンバーシツ
プ関数である。さらに、11は前記メンバーシツ
プ関数7及び10を合成したメンバーシツプ関
数、12はこのメンバーシツプ関数11から重心
をとつて得られる推論値であり、この例では操作
量の変化率Δuとして当該フアジイ推論装置より
出力される。
また、第6図はこのような動作原理に基づく従
来のフアジイ推論装置の一例を示すブロツク図で
あり、図において、13は前記各ルール1及び2
に対し、入力された特徴量3,4から前件部の成
立度を評価し、その成立度に基づいて後件部のメ
ンバーシツプ関数の重み付けを行う加重手段、1
4はこの加重手段13によつて重み付けされたメ
ンバーシツプ関数を合成する合成手段、15はこ
の合成手段14によつて合成されたメンバーシツ
プ関数から推論値12を決定して、操作量の変化
率Δuを出力する推論値決定手段である。
来のフアジイ推論装置の一例を示すブロツク図で
あり、図において、13は前記各ルール1及び2
に対し、入力された特徴量3,4から前件部の成
立度を評価し、その成立度に基づいて後件部のメ
ンバーシツプ関数の重み付けを行う加重手段、1
4はこの加重手段13によつて重み付けされたメ
ンバーシツプ関数を合成する合成手段、15はこ
の合成手段14によつて合成されたメンバーシツ
プ関数から推論値12を決定して、操作量の変化
率Δuを出力する推論値決定手段である。
次に動作について説明する。ここで、ルール1
は、「特徴量3(制御偏差e)が負に少しずれ、
かつ特徴量4(制御偏差の変化率Δe)が正に少
しずれているならば、、推論値12(操作量Δu)
を正に少しずらす」という意味であり、このう
ち、「…ならば」の部分を前述の前件部と呼び、
それより後の部分を前述の後件部と呼んでいる。
従つて、このルール1の前件部のメンバーシツプ
関数5は「負に少しずれた制御偏差の集合」を規
定し、メンバーシツプ関数6は「正に少しずれた
制御偏差の変化率の集合」を規定している。
は、「特徴量3(制御偏差e)が負に少しずれ、
かつ特徴量4(制御偏差の変化率Δe)が正に少
しずれているならば、、推論値12(操作量Δu)
を正に少しずらす」という意味であり、このう
ち、「…ならば」の部分を前述の前件部と呼び、
それより後の部分を前述の後件部と呼んでいる。
従つて、このルール1の前件部のメンバーシツプ
関数5は「負に少しずれた制御偏差の集合」を規
定し、メンバーシツプ関数6は「正に少しずれた
制御偏差の変化率の集合」を規定している。
今、加重手段13に入力された特徴量3として
の制御偏差の実際の値がe0である場合、第5図に
示すように、前記値e0が「負に少しずれた制御偏
差」である度合いは、メンバーシツプ関数5によ
つて“0.8”と評価され、前記値Δe0が「正に少し
ずれた制御偏差の変化率」である度合いは、メン
バーシツプ関数6によつて“0.7”と評価される。
これら両評価値の内、低い方の値“0.7”が採用
されて当該ルール1の前件部の成立度となる。さ
らに、このルール1の後件部のメンバーシツプ関
数7は「操作量を正に少しずらす」ということを
意味いており、当該メンバーシツプ関数7は、前
記前件部の成立度の値に従つて0.7倍に重み付け
される。
の制御偏差の実際の値がe0である場合、第5図に
示すように、前記値e0が「負に少しずれた制御偏
差」である度合いは、メンバーシツプ関数5によ
つて“0.8”と評価され、前記値Δe0が「正に少し
ずれた制御偏差の変化率」である度合いは、メン
バーシツプ関数6によつて“0.7”と評価される。
これら両評価値の内、低い方の値“0.7”が採用
されて当該ルール1の前件部の成立度となる。さ
らに、このルール1の後件部のメンバーシツプ関
数7は「操作量を正に少しずらす」ということを
意味いており、当該メンバーシツプ関数7は、前
記前件部の成立度の値に従つて0.7倍に重み付け
される。
このことはルール2についても全く同様であ
り、入力された特徴量3の制御偏差の実際の値e0
と、特徴量4の制御偏差の変化率の実際の値Δe0
に基づいて前件部の成立度を評価し、その成立度
の値“0.5”に基づいて、後件部のメンバーシツ
プ関数10に0.5倍の重み付けを行う。このよう
にして重み付けされたメンバーシツプ関数7及び
10は合成手段14へ入力されて合成され、合成
メンバーシツプ関数11が得られる。さらに、こ
の合成メンバーシツプ関数11は推論値決定手段
25に入力されて重心が計算され、その結果、操
作量の変化率Δu0が推論値12として当該フアジ
イ推論装置より出力される。
り、入力された特徴量3の制御偏差の実際の値e0
と、特徴量4の制御偏差の変化率の実際の値Δe0
に基づいて前件部の成立度を評価し、その成立度
の値“0.5”に基づいて、後件部のメンバーシツ
プ関数10に0.5倍の重み付けを行う。このよう
にして重み付けされたメンバーシツプ関数7及び
10は合成手段14へ入力されて合成され、合成
メンバーシツプ関数11が得られる。さらに、こ
の合成メンバーシツプ関数11は推論値決定手段
25に入力されて重心が計算され、その結果、操
作量の変化率Δu0が推論値12として当該フアジ
イ推論装置より出力される。
以上のように3フアジイ推論装置では複数のル
ールが同時に働き、各々の前件部の成立度に応じ
た後件部の重み付けがなされ、全体として釣り合
いのとれた値を推論値として出力している。
ールが同時に働き、各々の前件部の成立度に応じ
た後件部の重み付けがなされ、全体として釣り合
いのとれた値を推論値として出力している。
第7図は上述した従来のフアジイ推論装置を制
御装置として適用した制御系の構成を示すブロツ
ク線図である。図において、40は制御系の外部
より与えられる目標値(r)、41はプロセス36よ
り出力される制御量(y)、3は目標値40と制御量
41との制御偏差(e)、100は制御偏差(e)3を入
力して制御偏差の変化率(Δe)4を出力する微
分要素、101はフアジイ推論装置であり、前記
第6図の全体がこれに相当する。
御装置として適用した制御系の構成を示すブロツ
ク線図である。図において、40は制御系の外部
より与えられる目標値(r)、41はプロセス36よ
り出力される制御量(y)、3は目標値40と制御量
41との制御偏差(e)、100は制御偏差(e)3を入
力して制御偏差の変化率(Δe)4を出力する微
分要素、101はフアジイ推論装置であり、前記
第6図の全体がこれに相当する。
12はフアジイ推論装置101より出力される
推論値であり、この適用例では操作量の変化率
(Δu)となる。102は操作量の変化率Δuを入
力して操作量x43を上記プロセス36に出力す
る積分要素である。
推論値であり、この適用例では操作量の変化率
(Δu)となる。102は操作量の変化率Δuを入
力して操作量x43を上記プロセス36に出力す
る積分要素である。
第8図は例えば長田 正著「フイードバツク制
御」(オーム社、1971発行)の第16頁に示された
PID制御装置のブロツク線図であり、図におい
て、103は制御偏差(e)3を入力してそのまま出
力する形式要素、104は制御偏差(e)3を入力し
て積分値の1/TI倍を出力する積分要素、10
5は制御偏差(e)を入力して該制御偏差の変化率の
TD倍を出力する微分要素、106は形式要素1
03、積分要素104、微分要素105の各出力
の和を入力し、この和をKC倍した信号を操作量
(y)41として出力するゲイン要素である。
御」(オーム社、1971発行)の第16頁に示された
PID制御装置のブロツク線図であり、図におい
て、103は制御偏差(e)3を入力してそのまま出
力する形式要素、104は制御偏差(e)3を入力し
て積分値の1/TI倍を出力する積分要素、10
5は制御偏差(e)を入力して該制御偏差の変化率の
TD倍を出力する微分要素、106は形式要素1
03、積分要素104、微分要素105の各出力
の和を入力し、この和をKC倍した信号を操作量
(y)41として出力するゲイン要素である。
第9図は上記第8図に示すPID制御装置107
を制御装置として適用した場合の制御系の構成を
示すブロツク線図である。この場合、PID制御装
置107は制御偏差(e)3を入力して操作量(y)41
をプロセス36へ出力する。
を制御装置として適用した場合の制御系の構成を
示すブロツク線図である。この場合、PID制御装
置107は制御偏差(e)3を入力して操作量(y)41
をプロセス36へ出力する。
以上のように、従来の制御系は制御装置として
フアジイ推論装置あるいはPID制御装置を用いて
構成されている。
フアジイ推論装置あるいはPID制御装置を用いて
構成されている。
しかし、いずれの場合にも、制御対象であるプ
ロセス36の性質により、制御装置のチユーニン
グが必要である。従つて、制御装置としてフアジ
イ推論装置を用いた場合は、第5図の各メンバー
シツプ関数5〜10のチユーニングが必要であ
り、制御装置としてPID制御装置を用いた場合
は、第8図の各係数KC,TI,TDのチユーニング
が必要である。
ロセス36の性質により、制御装置のチユーニン
グが必要である。従つて、制御装置としてフアジ
イ推論装置を用いた場合は、第5図の各メンバー
シツプ関数5〜10のチユーニングが必要であ
り、制御装置としてPID制御装置を用いた場合
は、第8図の各係数KC,TI,TDのチユーニング
が必要である。
そこで、上記チユーニングを自動的に行うしく
みを付加したオートチユーニング制御系を構成し
ようとするならば、制御装置としては実績の高い
PID制御装置を適用し、その各係数KC,TI,TD
のチユーニングには、オペレータのノウハウを自
動化するのに適したフアジイ推論装置を利用する
ことが考えられる。
みを付加したオートチユーニング制御系を構成し
ようとするならば、制御装置としては実績の高い
PID制御装置を適用し、その各係数KC,TI,TD
のチユーニングには、オペレータのノウハウを自
動化するのに適したフアジイ推論装置を利用する
ことが考えられる。
しかし、チユーニング時にフアジイ推論装置の
入力となる量は、信号の瞬時値であつて常に入力
できる制御偏先(e)とその変化率(Δe)ではなく、
制御偏差の発散傾向、制御偏差の大きさ、目標値
の変化に対する制御量の追従度など、信号の傾向
をある時間観測して間欠的に抽出される制御系の
状態を示す特徴量(Si)である。
入力となる量は、信号の瞬時値であつて常に入力
できる制御偏先(e)とその変化率(Δe)ではなく、
制御偏差の発散傾向、制御偏差の大きさ、目標値
の変化に対する制御量の追従度など、信号の傾向
をある時間観測して間欠的に抽出される制御系の
状態を示す特徴量(Si)である。
第10図aは制御偏差eまたはその変化率Δe
の値の時間的経緯を示し、第10図bは制御系の
状態を示す特徴量(Si)の値の時間的経緯を示し
ている。
の値の時間的経緯を示し、第10図bは制御系の
状態を示す特徴量(Si)の値の時間的経緯を示し
ている。
上記特徴量(Si)は普段はSi=0であつて、信
号の傾向から制御系の状態が抽出された時のみ0
<Si≦1となる。この0という値は、真に0とし
ての意味がある0ではなく、特徴量(Si)が次に
抽出されるまでの間の便宜の値である。
号の傾向から制御系の状態が抽出された時のみ0
<Si≦1となる。この0という値は、真に0とし
ての意味がある0ではなく、特徴量(Si)が次に
抽出されるまでの間の便宜の値である。
したがつて、従来のフアジイ推論装置では、特
徴量(Si)を間欠的に入力したのでは、前件部
5・6,8・9での評価が行えないために推論は
不可能である。また、特徴量(Si)の全てが同時
に入力された場合に推論ができても、その推論は
その瞬間のみのことであり、特徴量(Si)の一部
のみが入力されても、それらを無駄にしてしま
う。
徴量(Si)を間欠的に入力したのでは、前件部
5・6,8・9での評価が行えないために推論は
不可能である。また、特徴量(Si)の全てが同時
に入力された場合に推論ができても、その推論は
その瞬間のみのことであり、特徴量(Si)の一部
のみが入力されても、それらを無駄にしてしま
う。
さらに、PID制御装置の係数Kc,TI,TDの値
は一定、変化したとしても緩慢な変化のみである
から、当然、推論される値は収束すべきである。
しかしながら、従来のフアジイ推論装置では、瞬
時値のみに反応し、過去の推論の蓄積は行われな
いので、この収束性は得られなかつた。
は一定、変化したとしても緩慢な変化のみである
から、当然、推論される値は収束すべきである。
しかしながら、従来のフアジイ推論装置では、瞬
時値のみに反応し、過去の推論の蓄積は行われな
いので、この収束性は得られなかつた。
すなわち、第1の従来技術である従来のフアジ
イ推論装置と、第2の従来技術であるPID制御装
置とを単純に組み合わせても、オートチユーニン
グ制御系を構成することは不可能であつた。
イ推論装置と、第2の従来技術であるPID制御装
置とを単純に組み合わせても、オートチユーニン
グ制御系を構成することは不可能であつた。
この発明は、上記のような問題点を解決するた
めになされたもので、オートチユーニング制御系
への適用は勿論、それ以外にも広く、間欠的にし
か入力値が意味を持たない場合、すべての入力値
がそろわない場合にも推論に利用することが可能
で、もし推論すべき値が一定または緩慢な変化し
かしない場合には、収束性がある推論を行えるよ
うなフアジイ推論装置を得ることを目的としてい
る。
めになされたもので、オートチユーニング制御系
への適用は勿論、それ以外にも広く、間欠的にし
か入力値が意味を持たない場合、すべての入力値
がそろわない場合にも推論に利用することが可能
で、もし推論すべき値が一定または緩慢な変化し
かしない場合には、収束性がある推論を行えるよ
うなフアジイ推論装置を得ることを目的としてい
る。
この発明に係るフアジイ推論装置は、合成メン
バーシツプ関数の合成に際して、各ルールの後件
部のメンバーシツプ関数だけでなく、前回の合成
メンバーシツプ関数も合わせて合成するようにし
たものである。
バーシツプ関数の合成に際して、各ルールの後件
部のメンバーシツプ関数だけでなく、前回の合成
メンバーシツプ関数も合わせて合成するようにし
たものである。
この発明におけるフアジイ推論装置は、現在だ
けでなく、過去の各ルールの働きも加えたメンバ
ーシツプ関数の合成を行うため、一時に多数のル
ールが働いた場合と同様な推論を可能とし、しか
も、その推論値は連続的な値をとり得るととも
に、合成メンバーシツプ関数に学習効果が生じ、
収束性のある推論をも可能とする。
けでなく、過去の各ルールの働きも加えたメンバ
ーシツプ関数の合成を行うため、一時に多数のル
ールが働いた場合と同様な推論を可能とし、しか
も、その推論値は連続的な値をとり得るととも
に、合成メンバーシツプ関数に学習効果が生じ、
収束性のある推論をも可能とする。
以下、この発明の一実施例を図について説明す
る。第1図はこの発明による再帰型のフアジイ推
論装置の一実施例を示すブロツク図、第2図はそ
のプロセス制御用制御装置の制御ゲインのチユー
ニングへの適用例を示すブロツク図、第3図はそ
の動作原理を示す説明図である。これら各図にお
いて、20及び21は推論のルール、22〜25
は当該再帰型のフアジイ推論装置に入力される特
徴量、26及び27は前記ルール20の前件部の
メンバーシツプ関数、28はその後件部のメンバ
ーシツプ関数、29及び30は前記ルール21の
前件部のメンバーシツプ関数、31はその後件部
のメンバーシツプ関数である。さらに、32は前
回に合成された不満足さを表す合成メンバーシツ
プ関数、33はこの不満足さを表す前回の合成メ
ンバーシツプ関数32に、プロセスの特性変化に
応じて重み付けした合成メンバーシツプ関数、3
4は前記ルール20の後件部のメンバーシツプ関
数28及びルール21の後件部のメンバーシツプ
関数31と、前記重み付けされた前回のメンバー
シツプ関数33を合成して得られる合成メンバー
シツプ関数で、「不満足な制御ゲインKc」のフア
ジイ集合を表しており、また、35はこの合成メ
ンバーシツプ関数34から得られる推論値で、こ
の例では制御ゲインKcとして当該再帰型のフア
ジイ推論装置から出力される。
る。第1図はこの発明による再帰型のフアジイ推
論装置の一実施例を示すブロツク図、第2図はそ
のプロセス制御用制御装置の制御ゲインのチユー
ニングへの適用例を示すブロツク図、第3図はそ
の動作原理を示す説明図である。これら各図にお
いて、20及び21は推論のルール、22〜25
は当該再帰型のフアジイ推論装置に入力される特
徴量、26及び27は前記ルール20の前件部の
メンバーシツプ関数、28はその後件部のメンバ
ーシツプ関数、29及び30は前記ルール21の
前件部のメンバーシツプ関数、31はその後件部
のメンバーシツプ関数である。さらに、32は前
回に合成された不満足さを表す合成メンバーシツ
プ関数、33はこの不満足さを表す前回の合成メ
ンバーシツプ関数32に、プロセスの特性変化に
応じて重み付けした合成メンバーシツプ関数、3
4は前記ルール20の後件部のメンバーシツプ関
数28及びルール21の後件部のメンバーシツプ
関数31と、前記重み付けされた前回のメンバー
シツプ関数33を合成して得られる合成メンバー
シツプ関数で、「不満足な制御ゲインKc」のフア
ジイ集合を表しており、また、35はこの合成メ
ンバーシツプ関数34から得られる推論値で、こ
の例では制御ゲインKcとして当該再帰型のフア
ジイ推論装置から出力される。
また、36は制御の対象となるプロセス、37
はこのプロセス36を制御する、例えばPID制御
装置のような制御装置、38はこの制御装置37
に推論値(制御ゲインKc)35を供給する、こ
の発明に係る再帰型のフアジイ推論装置、39は
このフアジイ推論装置38に特徴量22〜25を
供給する特徴量抽出器であり、40は制御系の外
部より与えられる目標値(r)、41はプロセス36
から出力される制御量(y)、42は制御装置37へ
入力される、前記目標値40と制御量41との変
化(e)、43は制御装置37よりプロセス36へ与
えられる操作量(x)、44はプロセス36からフア
ジイ推論装置38へ送られるプロセス特性変化量
である。
はこのプロセス36を制御する、例えばPID制御
装置のような制御装置、38はこの制御装置37
に推論値(制御ゲインKc)35を供給する、こ
の発明に係る再帰型のフアジイ推論装置、39は
このフアジイ推論装置38に特徴量22〜25を
供給する特徴量抽出器であり、40は制御系の外
部より与えられる目標値(r)、41はプロセス36
から出力される制御量(y)、42は制御装置37へ
入力される、前記目標値40と制御量41との変
化(e)、43は制御装置37よりプロセス36へ与
えられる操作量(x)、44はプロセス36からフア
ジイ推論装置38へ送られるプロセス特性変化量
である。
さらに、45は前記各ルール20及び21に対
して入力された特徴量22〜25から前件部の成
立度を評価し、その成立度に基づいて後件部のメ
ンバーシツプ関数の重み付けを行う加重手段、4
6はこの加重手段45によつて重み付けされた各
メンバーシツプ関数28,31と、前回の合成メ
ンバーシツプ関数33とを合成して、各推論値3
5(例えばゲインK)の不満足さの度合いを表す
新たな合成メンバーシツプ関数34を得る合成手
段、47はこの合成手段46によつて合成された
新たな合成メンバーシツプ関数34から推論値3
5を決定して出力する推論値決定手段、48は前
記合成手段46にて合成された新たな合成メンバ
ーシツプ関数34を遅延させる遅延手段、49は
非線形性が強いために、動作点により制御特性が
変化する場合の目標値信号、たとえば弁の開閉の
みが制御できる流量制御において、流体の粘性に
影響を与えて制御特性を変化させてしまう流体温
度のように直接に制御できない信号を、プロセス
特性変化量44として、前記プロセス36から入
力し、その制御特性の変化の度合いに応じて過去
の忘却のための加重係数を出力する特性変化評価
手段、50はこの特性変化評価手段49とともに
評価加重手段51を形成し、遅延手段48で遅延
された前回の合成メンバーシツプ関数32に、こ
の特性変化評価手段49から出力される荷重係数
値を乗算して重み付けを行い、合成メンバーシツ
プ関数33として前記合成手段46へフイードバ
ツクする乗算手段である。なお、第1図におい
て、二重線はフアジイ量、単一線は非フアジイ量
を表わしている。
して入力された特徴量22〜25から前件部の成
立度を評価し、その成立度に基づいて後件部のメ
ンバーシツプ関数の重み付けを行う加重手段、4
6はこの加重手段45によつて重み付けされた各
メンバーシツプ関数28,31と、前回の合成メ
ンバーシツプ関数33とを合成して、各推論値3
5(例えばゲインK)の不満足さの度合いを表す
新たな合成メンバーシツプ関数34を得る合成手
段、47はこの合成手段46によつて合成された
新たな合成メンバーシツプ関数34から推論値3
5を決定して出力する推論値決定手段、48は前
記合成手段46にて合成された新たな合成メンバ
ーシツプ関数34を遅延させる遅延手段、49は
非線形性が強いために、動作点により制御特性が
変化する場合の目標値信号、たとえば弁の開閉の
みが制御できる流量制御において、流体の粘性に
影響を与えて制御特性を変化させてしまう流体温
度のように直接に制御できない信号を、プロセス
特性変化量44として、前記プロセス36から入
力し、その制御特性の変化の度合いに応じて過去
の忘却のための加重係数を出力する特性変化評価
手段、50はこの特性変化評価手段49とともに
評価加重手段51を形成し、遅延手段48で遅延
された前回の合成メンバーシツプ関数32に、こ
の特性変化評価手段49から出力される荷重係数
値を乗算して重み付けを行い、合成メンバーシツ
プ関数33として前記合成手段46へフイードバ
ツクする乗算手段である。なお、第1図におい
て、二重線はフアジイ量、単一線は非フアジイ量
を表わしている。
次に動作について説明する。この再帰型のフア
ジイ推論装置における推論の目的は、プロセス3
6の特徴量22〜25を監視することで制御ゲイ
ンKcのチユーニングを行うことである。そこで、
まず上記特徴量22〜25を具体的に示す。即
ち、特徴量22は偏差(e)42の発散傾向Sa、特
徴量23は偏差42の大きさSb、特徴量24は
目標値(r)40の変化に対する制御量(y)の追従度
Sc、特徴量25は偏差42の大きさSd(=Sb)で
ある。このとき、ルール20は「もし偏差(e)42
の発散傾向が大きく、またその絶対値も大きけれ
ば、現在の制御ゲインKcは大きすぎると判断で
きる」という意味を持つている。ここで、加重手
段45に入力された特徴量22〜25の実際の値
がSa0,Sb0,Sc0,Sd0である場合、前記ルール
20において、この値Sa0が「大きい」かどう
か、値Sb0が「大きい」かどうかは、それぞれ当
該ルール20の前件部のメンバーシツプ関数26
及び27によつて評価される。第3図に図示の例
では、それぞれの評価値は“0.4”と“1.0”であ
り、この両評価の内、低い方の値“0.4”が当該
ルール20の前件部の成立度として採用される。
ジイ推論装置における推論の目的は、プロセス3
6の特徴量22〜25を監視することで制御ゲイ
ンKcのチユーニングを行うことである。そこで、
まず上記特徴量22〜25を具体的に示す。即
ち、特徴量22は偏差(e)42の発散傾向Sa、特
徴量23は偏差42の大きさSb、特徴量24は
目標値(r)40の変化に対する制御量(y)の追従度
Sc、特徴量25は偏差42の大きさSd(=Sb)で
ある。このとき、ルール20は「もし偏差(e)42
の発散傾向が大きく、またその絶対値も大きけれ
ば、現在の制御ゲインKcは大きすぎると判断で
きる」という意味を持つている。ここで、加重手
段45に入力された特徴量22〜25の実際の値
がSa0,Sb0,Sc0,Sd0である場合、前記ルール
20において、この値Sa0が「大きい」かどう
か、値Sb0が「大きい」かどうかは、それぞれ当
該ルール20の前件部のメンバーシツプ関数26
及び27によつて評価される。第3図に図示の例
では、それぞれの評価値は“0.4”と“1.0”であ
り、この両評価の内、低い方の値“0.4”が当該
ルール20の前件部の成立度として採用される。
さらに、このルール20の後件部は「大きすぎ
る制御ゲインKc」のフアジイ集合を定義してお
り、現在の制御ゲインKc≧Kc0は、少なくとも
前件部の成立度“0.4”の度合いで「大きすぎる」
といえる。そこで、この前件部の成立度“0.4”
に応じた重み付けをして、メンバーシツプ関数2
8を作成する。このことはルール21についても
全く同様であり、入力された特徴量24,25の
実際の値Sc0とSd0を前件部のメンバーシツプ関
数29,30で評価した成立度に基づいて後件部
のメンバーシツプ関数31を作成する。この例で
はメンバーシツプ関数31は全ての値が“0”で
ある関数となつている。
る制御ゲインKc」のフアジイ集合を定義してお
り、現在の制御ゲインKc≧Kc0は、少なくとも
前件部の成立度“0.4”の度合いで「大きすぎる」
といえる。そこで、この前件部の成立度“0.4”
に応じた重み付けをして、メンバーシツプ関数2
8を作成する。このことはルール21についても
全く同様であり、入力された特徴量24,25の
実際の値Sc0とSd0を前件部のメンバーシツプ関
数29,30で評価した成立度に基づいて後件部
のメンバーシツプ関数31を作成する。この例で
はメンバーシツプ関数31は全ての値が“0”で
ある関数となつている。
このように重み付けされたメンバーシツプ関数
28及び31と遅延手段48で1イテレーシヨン
分の遅延が与えられた前回の合成メンバーシツプ
関数32に評価加重手段51で重み付けした合成
メンバーシツプ関数33とを合成手段46で合成
し、新たな合成メンバーシツプ関数34を推論、
結果として得られる。
28及び31と遅延手段48で1イテレーシヨン
分の遅延が与えられた前回の合成メンバーシツプ
関数32に評価加重手段51で重み付けした合成
メンバーシツプ関数33とを合成手段46で合成
し、新たな合成メンバーシツプ関数34を推論、
結果として得られる。
この推論結果は不満の最も低いゲインとして、
新たな合成メンバーシツプ関数34の最も低いゲ
インに設定され、これが制御ゲインとして用いら
れるものである。
新たな合成メンバーシツプ関数34の最も低いゲ
インに設定され、これが制御ゲインとして用いら
れるものである。
ここで、第4図はその動作の流れを示すフロー
チヤートである。前回のイテレーシヨンで合成手
段46にて合成され、次回のイテレーシヨンの入
力として(ステツプST8)送出された合成メンバ
ーシツプ関数34は、遅延手段48に送られて1
イテレーシヨン分の遅延が与えられ、メンバーシ
ツプ関数32となる(ステツプST1)。また、こ
れとは別に、特徴量22〜25が加重手段45
に、プロセス36からのプロセス特性変化量44
が特性変化評価手段49にそれぞれ入力され(ス
テツプST2)、加重手段45は入力された特徴量
22〜25をルール20及び21の前件部のメン
バーシツプ関数26,27及び29,30で評価
し、求めた成立度に基づいて後件部のメンバーシ
ツプ関数28及び31を作成する(ステツプ
ST3)。
チヤートである。前回のイテレーシヨンで合成手
段46にて合成され、次回のイテレーシヨンの入
力として(ステツプST8)送出された合成メンバ
ーシツプ関数34は、遅延手段48に送られて1
イテレーシヨン分の遅延が与えられ、メンバーシ
ツプ関数32となる(ステツプST1)。また、こ
れとは別に、特徴量22〜25が加重手段45
に、プロセス36からのプロセス特性変化量44
が特性変化評価手段49にそれぞれ入力され(ス
テツプST2)、加重手段45は入力された特徴量
22〜25をルール20及び21の前件部のメン
バーシツプ関数26,27及び29,30で評価
し、求めた成立度に基づいて後件部のメンバーシ
ツプ関数28及び31を作成する(ステツプ
ST3)。
また、特性変化評価手段49は入力されたプロ
セス特性変化量44からプロセス36の特性変化
の度合いを評価してその評価値を乗算手段50へ
送り、遅延手段48で遅延された前回の合成メン
バーシツプ関数32にこの評価値を乗算して重み
付けを行い、重み付けされた合成メンバーシツプ
関数33を得る(ステツプST4)。
セス特性変化量44からプロセス36の特性変化
の度合いを評価してその評価値を乗算手段50へ
送り、遅延手段48で遅延された前回の合成メン
バーシツプ関数32にこの評価値を乗算して重み
付けを行い、重み付けされた合成メンバーシツプ
関数33を得る(ステツプST4)。
前記各ルール20及び21の後件部のメンバー
シツプ関数28及び31と、この重み付けされた
合成メンバーシツプ関数33とは合成手段46に
入力された合成され、新たな合成メンバーシツプ
関数34が生成される(ステツプST5)。
シツプ関数28及び31と、この重み付けされた
合成メンバーシツプ関数33とは合成手段46に
入力された合成され、新たな合成メンバーシツプ
関数34が生成される(ステツプST5)。
ここで、この合成演算には、和集合の演算が用
いられる。従つて、この合成メンバーシツプ関数
34は「大きすぎる制御ゲインKc」のフアジイ
集合と「小さすぎる制御ゲインKc」のフアジイ
集合との和集合であるから、結局「不満足な制御
ゲインKc」のフアジイ集合と解することができ
る。
いられる。従つて、この合成メンバーシツプ関数
34は「大きすぎる制御ゲインKc」のフアジイ
集合と「小さすぎる制御ゲインKc」のフアジイ
集合との和集合であるから、結局「不満足な制御
ゲインKc」のフアジイ集合と解することができ
る。
特にオートチユーニングの場合には、現在の値
をどう変えればよいのかがわからないし、現在の
値よりも良い値があるか否かもわからないことが
多い。このような場合には、今の値が大きすぎる
ために生ずる不満と、今の値が小さすぎるために
生ずる不満のみが評価である。一般に、このよう
な場合には、不満足さを後件部で表現することが
適切であり、この発明によるフアジイ推論装置
は、この後件部を用いたフアジイ推論を実行でき
る。
をどう変えればよいのかがわからないし、現在の
値よりも良い値があるか否かもわからないことが
多い。このような場合には、今の値が大きすぎる
ために生ずる不満と、今の値が小さすぎるために
生ずる不満のみが評価である。一般に、このよう
な場合には、不満足さを後件部で表現することが
適切であり、この発明によるフアジイ推論装置
は、この後件部を用いたフアジイ推論を実行でき
る。
生成された合成メンバーシツプ関数34は推論
値決定手段47に入力され、推論値決定手段47
は不満足さの度合が最も低い値」として、フアジ
イ集合のメンバーシツプ関数34が最低となる値
を選択することにより、推論値35としての制御
ゲインKc0を決定し、当該再帰型のフアジイ推論
装置38より制御装置37へ出力する(ステツプ
ST6)。具体的には、当該合成メンバーシツプ関
数34の値が一番低くなる制御ゲインを選べばよ
い。
値決定手段47に入力され、推論値決定手段47
は不満足さの度合が最も低い値」として、フアジ
イ集合のメンバーシツプ関数34が最低となる値
を選択することにより、推論値35としての制御
ゲインKc0を決定し、当該再帰型のフアジイ推論
装置38より制御装置37へ出力する(ステツプ
ST6)。具体的には、当該合成メンバーシツプ関
数34の値が一番低くなる制御ゲインを選べばよ
い。
不満足さを後件部で表すこの発明によるフアジ
イ推論装置では、合成メンバーシツプ関数は必ず
凹型となり、その値の最も低いところが最も不満
の少ない制御ゲインであるからである。
イ推論装置では、合成メンバーシツプ関数は必ず
凹型となり、その値の最も低いところが最も不満
の少ない制御ゲインであるからである。
次に、動作の打ち切りの判断を行い(ステツプ
ST7)、動作を継続する場合には処理をステツプ
ST8へ戻し、ステツプST5で得られた合成メンバ
ーシツプ関数を次回のイテレーシヨンの入力とす
る。
ST7)、動作を継続する場合には処理をステツプ
ST8へ戻し、ステツプST5で得られた合成メンバ
ーシツプ関数を次回のイテレーシヨンの入力とす
る。
なお、上記実施例では推論のルールの数が2つ
のものを示したが、その数は3つ以上であつても
よく、入力と出力の数、各ルールの前件部の条件
の段数も任意に設定できる。
のものを示したが、その数は3つ以上であつても
よく、入力と出力の数、各ルールの前件部の条件
の段数も任意に設定できる。
また、上記実施例ではプロセス制御用制御装置
における制御ゲインのチユーニングに適用した場
合について説明したが、その他のパラメータの推
定などに応用してもよく、上記実施例と同様の効
果を奏する。
における制御ゲインのチユーニングに適用した場
合について説明したが、その他のパラメータの推
定などに応用してもよく、上記実施例と同様の効
果を奏する。
以上のように、この発明によれば前回の合成演
算で得られた合成メンバーシツプ関数をフイード
バツクして、今回のメンバーシツプ関数の合成に
再度使用するように構成したので、特徴量{Si}
として、通常はSi=0で、特定の現象が起きた時
のみ0<Si≦1となるようなものを利用した場合
でも推論が可能となり、しかも、推論値が連続的
な値をとり得るばかりか、推定すべきパラメータ
が、一定もしくは緩慢な変化しかしない場合に
は、収束性のある推論値が得られる効果がある。
また目標値がよく変化する場合、次々に特徴量が
得られ、常に新しい情報を優先できるという効果
がある。
算で得られた合成メンバーシツプ関数をフイード
バツクして、今回のメンバーシツプ関数の合成に
再度使用するように構成したので、特徴量{Si}
として、通常はSi=0で、特定の現象が起きた時
のみ0<Si≦1となるようなものを利用した場合
でも推論が可能となり、しかも、推論値が連続的
な値をとり得るばかりか、推定すべきパラメータ
が、一定もしくは緩慢な変化しかしない場合に
は、収束性のある推論値が得られる効果がある。
また目標値がよく変化する場合、次々に特徴量が
得られ、常に新しい情報を優先できるという効果
がある。
第1図はこの発明の一実施例による再帰型のフ
アジイ推論装置の一実施例を示すブロツク図、第
2図はそのプロセス制御用制御装置の制御ゲイン
のチユーニングへの適用例を示すブロツク図、第
3図はその動作原理を示す説明図、第4図はその
動作の流れを示すフローチヤート、第5図は従来
のフアジイ推論装置の動作原理を示す説明図、第
6図はそのブロツク図、第7図は従来のフアジイ
推論装置を制御装置として適用した場合の制御系
の構成を示すブロツク図、第8図はPID制御装置
の構成を示すブロツク図、第9図はこのPID制御
装置を制御装置として適用した場合の制御系の構
成を示すブロツク図、第10図は制御偏差または
その変化率の時間的経緯および特徴量の時間的経
緯を示す図である。 20,21はルール、22〜25は特徴量、3
5は推論値、36はプロセス、38はフアジイ推
論装置、45は加重手段、46は合成手段、47
は推論値決定手段、48は遅延手段、51は評価
加重手段。なお、図中、同一符号は同一、又は相
当部分を示す。
アジイ推論装置の一実施例を示すブロツク図、第
2図はそのプロセス制御用制御装置の制御ゲイン
のチユーニングへの適用例を示すブロツク図、第
3図はその動作原理を示す説明図、第4図はその
動作の流れを示すフローチヤート、第5図は従来
のフアジイ推論装置の動作原理を示す説明図、第
6図はそのブロツク図、第7図は従来のフアジイ
推論装置を制御装置として適用した場合の制御系
の構成を示すブロツク図、第8図はPID制御装置
の構成を示すブロツク図、第9図はこのPID制御
装置を制御装置として適用した場合の制御系の構
成を示すブロツク図、第10図は制御偏差または
その変化率の時間的経緯および特徴量の時間的経
緯を示す図である。 20,21はルール、22〜25は特徴量、3
5は推論値、36はプロセス、38はフアジイ推
論装置、45は加重手段、46は合成手段、47
は推論値決定手段、48は遅延手段、51は評価
加重手段。なお、図中、同一符号は同一、又は相
当部分を示す。
Claims (1)
- 1 第1の値から第2の値の間の値をとるメンバ
ーシツプ関数を用いて前件部を形成し、それぞれ
プロセスの特徴量の状態による制御パラメータ値
の不満足さについて記述した複数のルールを有
し、各ルールに対し、プロセスにおいて時折現れ
る特徴量から前記前件部の成立度を評価し該成立
度に応じた重み付けをして前記後件部のメンバー
シツプ関数を作成する加重手段と、前回の合成メ
ンバーシツプ関数を出力する遅延手段と、この前
回の合成メンバーシツプ関数に前記プロセスの特
性変化の度合いに応じた重み付けを行う評価加重
手段と、この評価加重手段から出力された合成メ
ンバーシツプ関数と前記各ルールの後件部のメン
バーシツプ関数とを合成して新たな合成メンバー
シツプ関数を得て前記遅延手段に供給する合成手
段と、前記新たな合成メンバーシツプ関数の値が
最小となる値を選択して推論値としての制御パラ
メータを決定し出力する推論値決定手段とを備え
たフアジイ推論装置。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP62262031A JPH01103704A (ja) | 1987-10-16 | 1987-10-16 | フアジイ推論装置 |
| GB8819360A GB2211324B (en) | 1987-10-16 | 1988-08-15 | Fuzzy inference apparatus |
| CA000575019A CA1288168C (en) | 1987-10-16 | 1988-08-17 | Fuzzy inference apparatus |
| DE3832789A DE3832789A1 (de) | 1987-10-16 | 1988-09-27 | Stoerungsfolgevorrichtung |
| US07/950,166 US5432885A (en) | 1987-10-16 | 1992-09-24 | Recurrent fuzzy inference apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP62262031A JPH01103704A (ja) | 1987-10-16 | 1987-10-16 | フアジイ推論装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH01103704A JPH01103704A (ja) | 1989-04-20 |
| JPH0550002B2 true JPH0550002B2 (ja) | 1993-07-27 |
Family
ID=17370066
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP62262031A Granted JPH01103704A (ja) | 1987-10-16 | 1987-10-16 | フアジイ推論装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH01103704A (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5193144A (en) * | 1988-12-14 | 1993-03-09 | Shimano, Inc. | Fuzzy system |
| US5179634A (en) * | 1989-04-14 | 1993-01-12 | Omron Corporation | System for synthesizing new fuzzy rule sets from existing fuzzy rule sets |
-
1987
- 1987-10-16 JP JP62262031A patent/JPH01103704A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPH01103704A (ja) | 1989-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0617806B1 (en) | A closed loop neural network automatic tuner | |
| US5268835A (en) | Process controller for controlling a process to a target state | |
| EP2175332B1 (en) | Two-stage model predictive control technique | |
| CN1940780B (zh) | 过程控制系统中的在线自适应模型预测控制 | |
| US20040249483A1 (en) | Multiple-input/multiple-output control blocks with non-linear predictive capabilities | |
| Joseph | A tutorial on inferential control and its applications | |
| US6847954B1 (en) | Control-loop auto-tuner with nonlinear tuning rules estimators | |
| Haber et al. | Nonlinear internal model control using neural networks: An application for machining processes | |
| Alshammari et al. | A neural network-based adaptive backstepping control law with covariance resetting for asymptotic output tracking of a CSTR plant | |
| CN113625549A (zh) | 一种基于优化路径的pid参数设置方法与系统 | |
| JPH0550002B2 (ja) | ||
| GB2211324A (en) | Fuzzy inference apparatus | |
| JPH0550761B2 (ja) | ||
| JPH0551923B2 (ja) | ||
| Campos et al. | PD+ I Fuzzy Controller optimized by PSO applied to a variable dead time process | |
| JPH0235501A (ja) | ファジィ制御方式 | |
| JPH0551922B2 (ja) | ||
| US5649062A (en) | Auto-tuning controller and method of use therefor | |
| JPH04256102A (ja) | モデル予測制御装置 | |
| JPH06348305A (ja) | Pid調節器のパラメータ調整装置 | |
| Vega et al. | Model-free fuzzy controller for time delay systems | |
| JPS62241002A (ja) | オ−ト・チユ−ニング・コントロ−ラ | |
| JPH03268103A (ja) | オートチューニングコントローラ | |
| JPH0628005A (ja) | スライディングモード制御系を用いた制御方法 | |
| Al-Olimat et al. | Intelligent parallel controller based on fuzzy linguistic rules assessment |