JPH07219605A - ファジィ論理出力の非ファジィ化方法 - Google Patents

ファジィ論理出力の非ファジィ化方法

Info

Publication number
JPH07219605A
JPH07219605A JP6280711A JP28071194A JPH07219605A JP H07219605 A JPH07219605 A JP H07219605A JP 6280711 A JP6280711 A JP 6280711A JP 28071194 A JP28071194 A JP 28071194A JP H07219605 A JPH07219605 A JP H07219605A
Authority
JP
Japan
Prior art keywords
value
rule
output
membership function
area
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.)
Pending
Application number
JP6280711A
Other languages
English (en)
Inventor
Dinu Petre Madau
ペトレ マダウ ディヌ
Lee Albert Feldkamp
アルバート フェルドカンプ リー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ford Motor Co
Original Assignee
Ford Motor Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ford Motor Co filed Critical Ford Motor Co
Publication of JPH07219605A publication Critical patent/JPH07219605A/ja
Pending legal-status Critical Current

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16HGEARING
    • F16H61/00Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing
    • F16H61/02Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing characterised by the signals used
    • F16H61/0202Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing characterised by the signals used the signals being electric
    • F16H61/0204Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing characterised by the signals used the signals being electric for gearshift control, e.g. control functions for performing shifting or generation of shift signal
    • F16H61/0213Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing characterised by the signals used the signals being electric for gearshift control, e.g. control functions for performing shifting or generation of shift signal characterised by the method for generating shift signals
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0265Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
    • G05B13/0275Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion using fuzzy logic only
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/02Computing arrangements based on specific mathematical models using fuzzy logic
    • G06N7/04Physical realisation
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W50/00Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
    • B60W2050/0001Details of the control system
    • B60W2050/0043Signal treatments, identification of variables or parameters, parameter estimation or state estimation
    • B60W2050/0057Frequency analysis, spectral techniques or transforms
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16HGEARING
    • F16H61/00Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing
    • F16H2061/0075Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing characterised by a particular control method
    • F16H2061/0081Fuzzy logic
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S706/00Data processing: artificial intelligence
    • Y10S706/90Fuzzy logic

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Automation & Control Theory (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Molecular Biology (AREA)
  • Computational Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Algebra (AREA)
  • Mechanical Engineering (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computing Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Feedback Control In General (AREA)

Abstract

(57)【要約】 (修正有) 【目的】 単純化により高速化した非ファジィ化方法を
使用することにより、比較的廉価なファジィ論理コント
ローラで実現可能なファジィ論理処理モジュール101
を提供する。 【構成】 本ファジィ論理モジュールによる処理は、入
力信号の現在値が入力メンバーシップ関数を満足する程
度を表す「ルール強度値」を決定し、つぎに出力メンバ
ーシップ関数に割当られたルール強度値に基づき重心を
決定する。ルール強度と、それがコマンド出力値に及ぼ
す影響との間の非線形関係を反映するように、大きさが
調整される「影響値」で各出力メンバーシップ関数を表
すことにより、重心を決定する。事前に計算した重心に
より重み付けされた影響値を使用し、区分求積を使用せ
ずに重心530を正確に近似する。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は一般に電子制御システム
に関し、より詳細には、ファジィ論理コントローラにお
いて条件を満たす複数のファジィ論理ルール(satistie
d fuzzy logic rules )を非ファジィ化(defuzzifyin
g)して1つの出力コマンドの値を発生する方法に関す
る。
【0002】
【従来の技術および発明が解決しようとする課題】普
通、制御システムは、感知したシステムの状態の変化に
応答して、1つまたはそれ以上の出力コマンド信号の値
を変更するフィードバック機構を使用している。感知し
た信号の積分あるいは導関数に関する信号はもちろん、
感知したシステムの状態に直接比例する信号を処理(ma
nipulate)することにより、多くの動的システムを正確
に制御することができるのである。かかるフィードバッ
ク制御システムは、希望するシステムの性能のモデルを
つくり、実際のシステムの状態を感知して希望するモデ
ルと比較し、実際のシステムが、希望するモデルに従っ
て作動するようにエラー訂正フィードバック信号を発生
させることができるため、フィードバック制御システム
はうまく機能するのである。
【0003】しかし、精妙な比例・積分・微分(PI
D)コントローラであっても、より複雑な実システム
(physical system ) に対しては正確なモデルをつれ
くないことがしばしばあるのである。人間ならば容易に
実行できる多数の制御機能も、従来のフィードバックコ
ントローラによって自動化することができないのであ
る。たとえば、経験豊かな自動車ドライバが実行する
(ハンドルを回転する、ブレーキやアクセルペダルを押
す、ギヤをシフトする)機能は従来のコントローラによ
って処理することができない。その理由は、ドライバの
動きは、過去の経験をベースにして処理する複雑な入力
に基づいているからである。
【0004】Dr. Lotfi Zadehによって
開発され、「情報と制御」、第8巻、第3号、1965
年6月、頁338〜353の「ファジィ集合」に著され
ている「ファジィ論理」と呼ぶ手法は、複雑なシステム
を制御するときに非常に有効なことが証明された。ファ
ジィ論理を使用すると、入力変数はルールの集りと照ら
し合わせて処理されるが、各ルールは、たとえば「速度
が遅くかつ回転数が低いならば第1段ギヤにシフトせ
よ」という具合に、命題の形でシステムの目的を表現し
ている。条件が満足されているか、満足されていないか
の何れかであるという従来の論理とは異なり、「速度が
遅い」および「回転数が小さい」の2つの条件は部分的
に満足されているだけであるから、このルールも部分的
に満足されているだけである。結論として、入力条件が
部分的に満足されているだけである場合、「ロウギヤに
シフトせよ」に対するルールとその出力指令には、より
弱い「ルール強度(rule strength )」が割当られる。
ファジィ理論を使用するコントローラでは、多数のこの
ようなルールが同時に評価される。実際のギヤ選択は、
各ルールが満足した程度に従って満足した全ルールから
の全指令を組合わせることにより実行され、これによっ
て一致した指令(consensus directive )に到達する。
【0005】ファジィ論理制御システムにより、システ
ムによって想定できる起こりうる状態あるいは起こりう
る信号の値を、それぞれがメンバーシップ関数によって
定義された複数の「ファジィ集合」に分類することが可
能になる。したがって、与えられた信号に関連するメン
バーシップ関数は、ファジィ集合に関してその信号の現
在値が持っているメンバーシップの程度(degree-of-me
mbership:以下「メンバーシップ度」とする)を表示す
る。ファジィ集合では、ルールを、ルールの条件とルー
ルの指令とで表現する。たとえば、「速度が遅い」とい
う条件は、以下に述べる三角形型メンバーシップ関数
(triangular membership functions )によって表され
る。すなわち、4mphのとき最大のメンバーシップ度
の値が生じ、さらに−3mphと+11mphではゼロ
になるという三角形型メンバーシップ関数である。速度
の現在値が与えられると、第1のルールの条件の中で指
定されたメンバーシップ関数からメンバーシップ度の値
が生じ、つぎにこの値は第2の条件、「回転数が低い」
からのメンバーシップ度の値と組合わされて現在のルー
ル強度を決定する。また、ルールの指令、「ロウギヤに
シフトせよ」も、起こりうる各ギヤレシオに対するメン
バーシップ度の値を生じるメンバーシップ関数によって
表わされる「ロウギヤ」とともに、ファジィ集合として
表現される。
【0006】適用可能な全てのルールの中で「一致」を
決定することは、「非ファジィ化」(defuzzification
)と呼ぶ処理によって達成されるのが普通である。非
ファジィ化においては、各ルールの(複数の)条件を満
足した程度によって決定されるようなルール強度と、ル
ールの出力指令のメンバーシップ関数のタイプとを考慮
にいれて、全ての満足されたルールの重心(center of
gravity )が決定される。非常に多数のルールを使用す
る複雑なファジィ論理システムでは、希望する一致出力
の値に到達するために要求される計算は大量になる。制
御システムが実時間で動作しなければならないとき、急
速に変化する入力信号についていけるように、希望する
出力制御値を計算しなければならないという速度は、シ
ステムの計算用資源にかなりの負荷をかけることにな
る。
【0007】普通の非ファジィ化方法は制御プロセッサ
にかなりの負荷をかけており、入れ子型ループを多数回
繰り返すことが必要になっている。最も一般的な方法
は、各増分面積掛けるそのモーメントアーム(moment a
rm)の乗算を実行する区分求積法(piece-wise integra
tion)を使用して、出力重心(output centroid )用に
面積の中心を決定して最終的に重心の値を出す。出力コ
マンド信号の領域がη個の値に分割されていれば、積分
の過程でη個の増分面積が決定され、モーメントアーム
を計算するためにη回の乗算が要求される。この大量の
乗算に対する要求条件のため、ハードウエア乗算器を備
えた高価なマイクロコントローラ、あるいはカスタム集
積回路を使って、ファジィ論理コントローラを実現しな
ければならないことが多い。
【0008】本発明の目的は、ファジィ論理コントロー
ラによる非ファジィ化を単純化し、その結果として非フ
ァジィ化の速度を速めることである。
【0009】本発明の関連する目的は、費用があまりか
からない処理装置を使用して非ファジィ化を実行するこ
とである。
【0010】
【課題を解決する手段】本発明は、異なる信号源(sour
ce)からの複数の入力信号を1つの出力コマンド信号値
に変換する方法及び/或いはメカニズムの形式をとって
いる。制御メカニズムの目的を定める複数のルールは参
照するために格納されており、各ルールには、1つまた
はそれ以上の入力メンバーシップ関数および1つの出力
メンバーシップ関数から成るデータが含まれている。
【0011】本発明によれば、入力信号の処理が始まる
前に、各ルールの出力メンバーシップ関数は、事前に計
算された個別の面積中心(area-center )の値、すなわ
ち出力メンバーシップ関数の面積中心における出力信号
値を指定する値によって示される。さらに、与えられた
メンバーシップ関数のタイプに対する出力ルール強度を
影響値(influence value )に変換する影響伝達関数
(ITF:influence transfer function )が与えられ
ており、この影響値は、そのルール強度が出力コマンド
信号の値に影響を及ぼす程度を示すものである。
【0012】信号を処理している間に、入力信号の現在
値は各ルールのメンバーシップ関数によって評価され、
出力ルール強度の値が生じる。各ルールの出力ルール強
度は、そのルールの入力メンバーシップ関数が入力信号
の現在値によって満足されている程度を示している。本
発明によれば、つぎに各ルールに関連する影響伝達関数
は、かかる各ルール強度値を、後で行う処理に使用され
る影響値の現在値に変換する。
【0013】各メンバーシップの面積の形状に対して使
用される影響伝達関数(ITF)は、ルール強度値を影
響値に変換(transform )するが、この影響値は、形状
の相対面積と、ルール強度の大きさが希望する出力信号
に及ぼす影響とルール強度の大きさとの間の非線形関係
とを反映している。ITFの変換はアルゴリズムによる
処理によって実行することができるが、1ステップのテ
ーブルルックアップとして実行されることが望ましい。
この1ステップのテーブルルックアップは、(選択した
メンバーシップ関数のサイズと形状を識別するため)ル
ール番号によって索引がつけられた2次元のテーブルを
使用して実行される。
【0014】本発明のさらに優れた特徴によれば、隣接
するメンバーシップ関数の間のオーバーラップを表す負
の影響値も算出される。共通部分の影響値(intersecti
on influence value)は、隣接する出力メンバーシップ
関数を持つ2つのルールのうち小さいほうのルール強度
値を使用してつくられる。この小さいほうのルール強度
はオーバーラップ面積に対するITFによって変換さ
れ、得られた負の影響値はその後の処理の中で使用され
る。
【0015】本発明のさらに別の特徴によれば、全面積
値(total area value)は、隣接する出力メンバーシッ
プ関数の間でオーバーラップする共通部分の面積を示す
負の影響値を加算することにより減少する全面積値とと
もに、個々の影響値を累積することによりつくられる。
【0016】本発明によれば、全モーメントアームの値
は、各ルールに対するモーメントアームの値を累積し、
共通部分の各面積に対するモーメントアームの値を減算
することによりつくられる。ルールのモーメントアーム
の値は、そのルールの影響値と事前に計算された面積中
心の値との積によりつくられる。 共通部分の面積の中
心は、共通部分の面積の底線の中点と、2つの隣接する
メンバーシップ関数が交差する点の中間点を決定するこ
とにより求められる。
【0017】現在の入力信号の集合に対する出力コマン
ド信号の最終値は、全モーメントアームの値を全面積の
値で割ることにより求められる。
【0018】本発明によって考案された非ファジィ化方
法は、従来のファジィ論理コントローラで使用されてい
る重負荷の区分求積法によって生じる出力コマンド信号
値に非常に近似した出力コマンド信号値を発生するが、
必要とする処理は大幅に少ない。
【0019】本発明のその他の目的、特徴および利点も
含めて、これらの目的、特徴および利点は、以下に述べ
る詳細な説明を検討することによりより良く理解できる
であろう。
【0020】
【実施例】以下に添付の図面を参照して、本発明の好適
実施例を詳細に説明する。
【0021】概要 ここに説明する本発明の実施例は、添付図面の図1に一
般的に示す通り、汎用ファジィ論理制御機能を実行する
ようにプログラムされたマイクロコントローラを利用し
ている。ファジィ論理モジュール101はいくつかの信
号源103、105、107から入力信号を受信する
が、各信号は109で前処理されて正規化されたディジ
タル信号になる。このディジタル信号は整数値をとり、
ファジィ論理入力110によりモジュール101に供給
される。モジュール101は、112に一連のディジタ
ル整数出力値を発生するためにディジタル入力信号11
0の値を繰り返し処理する。この出力値のシーケンス
は、113で適切な後処理によって、図1に見られるア
クチュエータ115などの出力デバイス(utilizationd
evices )を駆動する制御信号に変換される。
【0022】論理モジュール101は、117のメモリ
に格納されているルール群に定義されている条件と11
0における入力信号を比較することにより、112で出
力信号の整数値を生成する。ルールに表現されている条
件は、119に格納されている入力メンバーシップ関数
を識別する。つぎに、119からの識別されたメンバー
シップ関数は120で使用され、ルールの入力条件が、
110における1つまたはそれ以上の入力信号値によっ
て満足されている程度を示すルール強度値を決定する。
得られたルール強度値は122で一時的に記憶される。
【0023】これらの記憶されたルール強度値と128
に格納されている出力メンバーシップ関数とを使用し
て、モジュールはつぎに「重心」の値を124で決定
し、出力整数値を112に発生する。
【0024】たとえば1993年6月6日出願の米国特
許出願出願番号第08/085799に説明されている
汎用ファジィ論理コントローラなど、以前の装置では、
重心の決定は、122で記憶されたルール強度値によっ
て制約をうける各メンバーシップ度の値とともに、11
8に格納されている出力メンバーシップ関数によって発
生するメンバーシップ度の値を、起こりうる出力整数値
の範囲にわたって積分することにより一部分は実行され
ている。つぎに積分の結果が処理されて、非ファジィ化
された出力整数値として112に送られる重心値がつく
られる。重心値を発生するこの以前の区分求積法につい
て、添付図面の図11に関連して以下に詳細に考察す
る。
【0025】本発明によれば、図1の124に示す重心
の決定は改良されたメカニズムによって達成されるが、
このメカニズムはファジィ論理コントローラに課せられ
る計算量の負荷を大幅に減少し、小容量のマイクロプロ
セッサが同等の結果を得ることを可能にしている。添付
図面の図13〜図17に関連して以下に詳細に説明する
通り、改良された非ファジィ化メカニズムは、事前に計
算された値を使用して、図1の118で見られる出力メ
ンバーシップ関数を表し、マイクロコントローラの読出
し専用メモリ(ROM)の中の、1つまたはそれ以上の
ルックアップテーブルの形式をとる影響伝達関数により
ルール強度122を処理し、さらに出力信号の全領域に
わたって区分求積法を実行することを必要とせずに、希
望する出力コマンド値を決定する。
【0026】非ファジィ化処理によって生じた出力整数
112は、程度は変わるとしても、110の入力信号値
が満足する条件を持つ全てのルールが一致したことを示
しており、113の後処理によって、図1に示すアクチ
ュエータ115などの出力デバイスを制御するのに適切
な信号の大きさに変換される。
【0027】標準メンバーシップ関数 本発明の原理は、たとえば図1に示す汎用ファジィ論理
システムのような、従来のマイクロプロセッサを使用し
たファジィ論理コントローラのコストの低減に便利に使
用することができる。ファジィ論理モジュールは、メン
バーシップ関数によって定義される「ファジィ集合」を
使用するが、所定の様式を備えたメンバーシップ関数が
メモリに格納されており、入力信号および希望する出力
信号によって想定されうる値を特性づけている。標準的
な処理方法が使用できるように、各入力信号は、ファジ
ィ論理モジュールに供給される前に109で前処理され
るので、所定の範囲の、たとえば0〜40の起こりうる
値を持つディジタル整数の形をした正規化された信号レ
ベルが110に供給される。
【0028】その上、各入力信号レベルによって想定さ
れうる起こりそうな値を特徴づけるために、所定の数の
ファジィ集合が使用される。好適実施例においては、各
信号値は、7つ集合のうちの1つまたはそれ以上のメン
バーであり、この7つの集合は意味のある指定(semant
ic designation)、「負の大」、「負の中位」、「負の
小」、「ゼロに近い」、「正の小」、「正の中位」、
「正の大」によって示されている。与えられた集合の内
部にある起こりうる各信号値は、所定の範囲を持った整
数として表現されるメンバーシップ度の値に関連してい
るのに対して、与えられた集合の外部にある信号値はメ
ンバーシップ度の値ゼロが割当てられる。ここに説明し
ている好適実施例においては、集合の最大のメンバーシ
ップ度は整数32によって表されている。メンバーシッ
プ関数は、各整数信号値を、その信号に関連するファジ
ィ集合の指定されたどれかに対するメンバーシップ度の
値に翻訳するため、メモリに格納されている。
【0029】三角形型メンバーシップ関数 いくつかの入力信号に必要なメンバーシップ関数を定義
するために必要な膨大なデータ点(data points )を効
率よく格納するため、本好適実施例は三角形の3つの頂
点を定義することによって各メンバーシップ関数を格納
している。三角形型メンバーシップ関数は、定義されて
いるファジィ集合の内部にある信号値の外側の限界を示
す指定された信号値により左と右の境界が決められてい
る。第3の指定された信号値は、その集合のメンバーシ
ップ度が最大になる信号値を示している。3個の整数値
L 、XC 、XR によって決まる1つの代表的な標準三
角形型メンバーシップ関数が図2に示されている。図3
は格納される整数のアレーを示す表であり、これらの整
数は各入力信号を評価するために使用される7つの三角
形型ファジィ集合を定義している。
【0030】ファジィ集合の三角形型メンバーシップ関
数のメンバーシップ度の値は、信号値を表す垂直の線
と、三角形との交点の位置を決めることにより決定され
る。処理ルーチンは、与えられた入力値が三角形の境界
内にあるか否かを第1に決定する。入力値が三角形の外
部にあれば、入力値がその集合に属していないのである
からメンバーシップ度の出力値はゼロであることがわか
る。入力が三角形と交差していると決定されると、ルー
チンは入力が最大値あるいは「中心値」XC の左にある
か右にあるかを判断する。つぎにルーチンは、傾斜線方
程式を使用してX座標が与える交点のY(メンバーシッ
プ度)座標を決定する。中心の右にある入力(すなわち
C より大きい入力)に対しては、
【数1】 中心の左にある入力(すなわちXC より小さい入力)に
対しては、
【数2】 スケール因子32は、動的範囲内で与えられたルール強
度が適切な分解能になるように選ばれる。「中心値」X
C における三角形型メンバーシップ関数の最大メンバー
シップ度は、2つの外側限界XR 、XL の中間にある必
要はない。たとえば、ファジィ集合「負の大」は直角三
角形で表されるかもしれないが、その場合XC もXL
共に、整数変数によって想定されうる最小値、0に等し
い。
【0031】1グループの7つの三角形型メンバーシッ
プ関数に対する各頂点の指定は各入力関数に関連してい
ることが望ましい。7つの三角形型集合を定義する整数
のアレーは添付図面の図3に示されている。
【0032】さらに各メンバーシップ関数は、タイプイ
ンジケータによってメモリの中で識別される。たとえ
ば、値1は三角形型メンバーシップ関数を指定し、値2
はテーブルルックアップメンバーシップ関数を指定し、
値3はベル型関数を指定する。
【0033】テーブルルックアップメンバーシップ関数 計算を速くするため、118に記憶されており頻繁にア
クセスされる出力メンバーシップ関数は、メモリの中で
ルックアップ値により表されていることが望ましい。各
テーブルルックアップメンバーシップ関数は、数列とし
てメモリに格納されているが、この数は、範囲の左およ
び右の外側の限界の組で、メンバーシップ度の値をゼロ
とするXR 、XL の組に始まり、N=(XR −XL
1)がその後に続く。このNは、それぞれ出力信号値
(XL +1)から(XR −1)に関するメンバーシップ
度の値である。
【0034】テーブルルックアップメンバーシップ関数
は、(どのような適切な形もとることができるが)その
代表的な形が図4に示されている。ルックアップ関数の
データ値は、図5のメモリマップが示すようなメモリの
中に格納されている。各関数の第1の2つの数XR 、X
L はファジィ集合の外側限界を指定し、各集合の残りの
値M1〜MNは中間にある点のメンバーシップ度の値を
表している。図4にグラフで示しているルックアップの
値は、格納されているデータポイント、0、12、2
9、27、25、21、19、16、13、11、8、
5および3によって定義される「負の大」メンバーシッ
プ関数に対して適切であろう。格納されているテーブル
ルックアップメモリに関しては、メンバーシップ関数の
形状に制約がない。しかし、通常は、連続するメンバー
シップ関数を備えていることが望ましいので、入力にお
ける変化により出力に予期しない変化が生じることはな
い。
【0035】ベル型メンバーシップ関数 また、所定の形状をしたメンバーシップ関数は、格納さ
れているテンプレート情報から実行時に動的につくられ
る。たとえば、図6に示す形をした標準のベル型メンバ
ーシップ関数を、図7(a)に見られるような標準「高
分解能」のルックアップテーブル構造体の中に格納する
ことができる。図7(b)に見られる通り、上部限界X
R および下部限界XL は、出力メンバーシップ関数の意
味のある指定に関連する7つのファジィ集合のそれぞれ
に対して格納されている。実際のルックアップ関数は高
分解能の関数から内挿法によって引き出され、その関数
を7つのファジィ集合のそれぞれに対する適切な限界に
スケーリングする。高速の処理が要求されていなけれ
ば、この内挿法は、各入力信号値の処理の一部として達
成される。もう1つの方法として、内挿法を処理の初期
化の間に実行して、与えられた信号変数に関する7つの
ファジィ集合のそれぞれに対して、添付図面の図5に示
すように構成された完全なルックアップデータポイント
アレーを生成してもよい。
【0036】ルールの格納 ファジィ論理モジュールの挙動は、図1の117に見ら
れるようなメモリに格納されているルールの集合によっ
て決定される。ルールを定義するデータ構造は次に示す
所定の記憶場所、すなわち、 1.与えられた非ファジィ化された出力信号を発生させ
るために使用されるルールの数を示す値、NoOfRu
lesを含む記憶場所 2.ルール番号によって索引がつけられたルールアドレ
スのアレー 3.1つまたはそれ以上の入力信号状態および出力指令
を定義する各ルールの可変長レコード構造 から構成されている。
【0037】ルールによって定義された各入力信号の状
態は、番号によって特定の入力信号を指定する入力識別
子と、状態を満足するため識別された入力信号の現在値
が所属しなければならない特定のファジィ集合を指定す
る入力メンバーシップ関数識別子(たとえば、nl、n
m、ns、ze、ps、pmあるいはpl)から構成さ
れている。したがって、入力条件は、入力3の現在値が
いつでも集合ps(正の小)のメンバーであることが満
足する条件「3、ps」としてルールに表現されるであ
ろう。
【0038】同様に、出力指令は出力識別子の値とメン
バーシップ関数との組合せで構成される。たとえば、
「output、pm」という出力指令は、ルールの入
力条件が満足される場合、出力信号が集合pm(正の中
位)のメンバーでなければならないことを示している。
【0039】ルールには1つ以上の入力条件が含まれて
よい。したがって、「1、nl、2、ze、outpu
t、ns」というルールを、「もし入力が負の大で、か
つ入力2がゼロに近いであれば、出力は負の小でなけれ
ばならない」という言語的命題として記述することがで
きる。
【0040】添付図面の図8は、格納されているルール
の集まりが与えられると、ファジィ論理モジュールは
「ルール強度」を計算するために要求される情報構造を
配置することを示している。
【0041】ルール自体はそれぞれのアドレスでメモリ
に中に格納されており、アドレスは201のルールアド
レスのアレー中に格納されている。このように、ルール
は、1から値NoOfRulesを含むNoOfRul
esまでのルール番号によってアクセスされるが、記憶
されている最終ルールを指定するため、NoOfRul
esはメモリに格納されている。図8に示す例では、ア
レー201は、第2のルール202のアドレスを生じる
ために、ルール2によって索引がつけられる。第2のル
ール202には1対の入力条件と1つの出力指令とが含
まれている。第1の条件は、アレー203の行4を指定
するため、入力1=4がテストされる特定の入力信号を
識別するが、110の中で定義された通り、アレー20
3の行4には信号番号4に対するメンバーシップ関数識
別情報が含まれている。アレー203の列は、ルール2
02の中の第1のメンバーシップ関数識別子Memb1
によって指定される。図8で、Memb1の値は「p
s」(正の小)であるから、セル204を指定する。セ
ル204には2つの値が含まれている。すなわち、第1
の値「1」はメンバーシップ関数のタイプ(三角形型)
を指定するが、このメンバーシップ関数は第4の入力信
号の「ps」(正の小)集合を定義するために既に選択
されている。セル204の第2の部分には、三角形型関
数を定義するデータが含まれているが、この場合は、2
05に格納されている3つの頂点の値、XL 、XC 、X
R によって定義される三角形型関数である。同じよう
に、ファジィ論理モジュールは、感知した入力信号を希
望する出力信号に変換するファジィ論理変換を支配する
すべてのルールを評価するため、各メンバーシップ関数
のタイプと場所およびそれに関連する入力信号を識別す
ることができる。
【0042】ルール強度の決定 図1のブロック図に示すように、ファジィ論理モジュー
ルは、第1にルールの強度を決定することにより、入力
の現在値を与える110の入力信号を処理し、決定した
ルール強度を122に記憶する。ルール強度が決定され
ると、プロセスは124で重心の値を決定する。この処
理は、図9の流れ図でより詳細に示されている。
【0043】図9に見られる通り、ルールの処理は22
0で開始する。各ルールが評価される前に、221でル
ールのルール強度変数が1に初期化される。つぎに各ル
ールの条件が222で識別され、224でそのメンバー
シップ関数が配置され、226でその条件に関連する入
力信号の現在値が求められる。つぎに配置されたメンバ
ーシップ関数が与える関連入力信号値によって指定され
たメンバーシップ関数の値が228で決定される。
【0044】図9に示す通り、1つ以上の条件を含むそ
れらのルールは、論理的AND演算子を使用して、それ
らの条件を組合せる。論理的OR演算子を使用して条件
を組合せるルールは使用されず、その代わりに、別々の
条件を表現する個別のルールを使用した等価関係が確立
される。このようにルールの中の条件の間の関係を制約
することによって処理が単純化され、ルールに条件と条
件との間の論理関係を含める必要がないので、ルールに
よる記憶スペースは少ししか使用されない。
【0045】AND演算と等価なファジィ論理は、ルー
ルの中の複数の条件から最小状態(MIN)のメンバー
シップの値を選ぶことにより実行される。したがって、
全体としてどれかの条件を満足しなければ、(つまり、
ゼロ値に評価されれば)結果はゼロである。単一条件に
対しては、各ルールの「ルール強度」は、単純に云って
識別された1つのメンバーシップ関数が与える入力の現
在値であるに過ぎない。ルールが1つ以上の条件を持っ
ていれば、後の処理に使用するためセーブされるルール
強度は、236でテストされ237でセーブされること
により配置される最小メンバーシップの値である。与え
られたルールの中の諸条件が論理的OR演算子によって
結合されたならば、妥当な結果の論理の組合せは、最大
状態(MAX)のメンバーシップの値によってつくられ
るであろう。
【0046】各ルール強度が格納されると、現在のルー
ル番号は、もっと多くのルールを処理すべきか否かを決
定するために、238でNoOfRulesと比較され
る。そうでなければ、次のステップ240で、重心の決
定の処理を続ける。
【0047】上に説明した処理ステップは図10の
(a)から(f)にグラフで示されている。これは、2
つのルールに基づいて、位置信号Pと偏差信号Dの2つ
の入力信号から、求める速度Vを示す出力信号を発生さ
せる制御システムを示している。図10(a)〜図10
(c)は、「P、ps、D、ps、output、p
s」すなわち、「もし入力位置Pが正の小で、かつ偏差
Dが正の小ならば、速度Vは正の小でなければならな
い」と書くことができる第1のルールを図示している。
図10(d)〜図10(f)は、「P、ze、D、p
s、output、ze」あるいは「もしPがゼロに近
いで、かつDが正の小ならば、Vはゼロに近いである」
と書くことができる第2のルールを図示している。図1
0では、例として入力信号の値P、Dが、縦の点線P、
Dで示されている。
【0048】図10(a)に見られる第1の条件は、位
置信号Pが正の小であれば満足している。Pは三角形型
メンバーシップ関数と交差し、水平の線301で示され
るメンバーシップ度の値をつくり出す。第1のルールに
対する第2の条件に対して、少し小さい振幅の入力信号
値D(偏差:deviation ) も、図10(b)で見られ
る三角形型メンバーシップ関数を満足しており、メンバ
ーシップ度の値303をつくり出す。値303は値30
1より小さいので、値303は第1のルールに対して記
憶された「ルール強度」となる。
【0049】また、図10(d)および図10(f)に
図示されている通り、PおよびDは共に第2のルールの
入力条件を満足している。図10(d)において、第1
の条件からのより小さいメンバーシップ度の値305
は、図10(e)のメンバーシップ関数によって生じた
値307より小さいので、305によって示された値は
第2のルールに対するルール強度として記憶される。記
憶されたこれら2つのルール強度の値は、次に考察する
通り、重心を決定するために使用される。
【0050】重心の決定 区分求積法の使用 図1の122および124で示した通り、入力信号の現
在値に対する各ルールの強度が決定されると、モジュー
ルは「重心」、すなわち、ゼロでないルール強度の値に
より示されるように、満足されているルールの全てが
「非ファジィ化された」一致を表す重心の値を決定す
る。
【0051】非ファジィ化された出力値は、図10
(g)によって示され、図12で拡大して示されている
数値積分の処理によって決定される。出力値は、そのル
ールの強度によって制約されている通り、満足された各
ルールの出力メンバーシップ度の最大値の点から点をつ
ないで形成される曲線の下の面積を積分することによっ
て決定される。この曲線は図10(g)および図12に
見られるとおりであり、2つの台形330、340の最
大の形を表している。台形330は、第1のルールの強
度値303によって制約された通り、図10(c)に見
られる第1のルールの三角形型出力メンバーシップ関数
によって形成された曲線である。小さいほうの台形34
0は、第2のルールの強度値305によって制約された
通り、図10(f)に見られる第2のルールの三角形型
出力メンバーシップ関数によって形成されている。第1
のルールは第2のルールよりも満足の程度が高かったの
で、出力値Vは集合「正の小」の内部になければならな
いことを示しており、第2のルールは満足の程度が低か
ったので、出力値Vは集合「ゼロに近い」の内部になけ
ればならないことを示している。非ファジィ化のプロセ
スにより、一致した値を示す非ファジィ化された結果の
信号がを与えられるが、この信号は、図10(g)およ
び図12を見て判るように、2つの台形330、340
の最大の形によって形成された曲線の下の面積の重心3
50である。
【0052】以前に使用した重心の値を決定する区分求
積法は、添付図面の図11の流れ図でより詳細に説明さ
れている。この手続きは数値積分によって2つの値を発
生する。すなわち、全モーメントアームの値MArmと
全面積の値Areaであり、図11の400に示す開始
点でルーチンが開始されると直ちに、両者ともに401
でゼロに初期化される。積分は、出力変数Xを0からM
axXまで変化させることにより繰り返される。変数R
uleNo、MaxStrengthの値は、403で
それぞれ1、0に初期化され(かつXが新しい値になる
都度、403でこれらの値に再初期化され)る。
【0053】(図9とともに上で考察した通り)Rul
eNoにより識別されたルールに対して前に決定された
ルール強度がつぎに405でチェックされる。ルール強
度がゼロであれば、そのルールは出力値に寄与しないの
で、直ちに408にジャンプし、ここでそのルールがR
uleNoに等しいか否かを決定するため、RuleN
oの現在値がチェックされ、処理すべき最後のルールを
識別する。RuleNoがNoOfRulesより小さ
ければ、410でRuleNoは次のルール番号に増分
されて405に戻り、そのルールの処理を続ける。
【0054】RuleStrengthがゼロでなけれ
ば、制御は412に渡され、ここでメンバーシップ度の
値LocalStrengthが与えるXの値を決定す
るため、RuleNoにより指定されたルールの出力メ
ンバーシップ関数が使用される。つぎに414でLoc
alStrengthが、RuleNoにより指定され
たルールに対して前に決定されたRuleStreng
thの値と比較される。そのルールのルール強度がLo
calStrengthよりも小さければ、Local
Strengthは416で小さいほうの値に交換され
る。この時、LocalStrengthは、つぎのこ
とを考慮して、そのルールの指令の値を表す整数にな
る。すなわち、入力信号の現在値、それらの入力信号値
がルールの入力条件を満足している程度、そしてXの現
在値に対するそのルールの出力メンバーシップ関数によ
り指定された値を考慮するのである。
【0055】414でテストされた後、得られたルール
指令の値LocalStrengthがMaxStre
ngthより大きければ、422でMaxStreng
thの前の値をLocalStrengthに交換する
ので、MaxStrengthは、与えられたXの値に
対して全てのルールが処理された後、(408でテスト
されたとき、もはやRuleNoがNoOfRules
よりも小さくない場合を示すように)ルールの1つによ
って発生した指令の最大値を保持する。
【0056】したがって、図12のグラフに見られる通
り、上に説明したルーチンはXの各値に対して値Max
Strengthを決定し、MaxStrength
は、図12の影をつけた四角形360によって示される
積分の区分面積(subarea )の高さに対応する。Xの変
化に伴ってつぎつぎにできるMaxStrengthに
よって定義される曲線の下の全面積は、図11の440
で示されているように変数AreaにMaxStren
gthの値を累積することにより決定される。450で
値MArmは区分面積のモーメントアームにより増分さ
れるが、区分面積のモーメントアームは積(X*Max
Strength)に等しい。つぎにXの値が積分区間
の最終MaxXに到達したか否か460でテストされ
る。到達していなければ、470でXは値Stepsi
zeだけ増分され、制御は403に戻り、処理すべき次
のXの値に対するRuleNoとMaxStrengt
hを再初期化する。
【0057】積分区間が終結すると、480でMArm
をAreaで割算することにより、重心の値COAが決
定される。つぎにモジュールは490で、非ファジィ化
された希望する出力信号を表すこの値を、適切な後処理
用ユーティリティルーチンに渡す。
【0058】つぎに説明するように、本発明によって考
案された非ファジィ化メカニズムは、図11に示した区
分求積法を大幅に改良していることを表している。つぎ
にこの改良案がベースにしている原理をより良く理解す
るためには、オーバラッピング出力メンバーシップ関数
と、かかるオーバラッピング出力メンバーシップ関数が
改良案で表されている方法との間の関係によく注意しな
ければならない。
【0059】三角形型メンバーシップ対単集合(Single
tons) これまでに説明した汎用ファジィ論理システムと同様、
ほとんどのファジィ論理システムは、出力重心を定義す
るために三角形型および台形型出力メンバーシップ関数
を使用している。出力メンバーシップはオーバラップさ
れて重心(centroid)を発生させるが、この面積の中心
がファジィ論理ブロックの出力である。メンバーシップ
関数は代わりに「単集合」(singletons)によって表さ
れるが、単集合は、出力メンバーシップを離散的出力領
域の値およびそのメンバーシップのルール強度とによっ
て定義される。面積あるいは面積のモーメントアームを
決定する必要がないので、単集合を使用することはファ
ジィ論理プロセッサが負担する計算量を大幅に減少する
ことができる。しかし、単集合だけでは希望する出力量
の貧弱な近似を与えるに過ぎず、単集合はほとんどのア
プリケーションに不適切である。
【0060】単集合を使用することから來る正確性の欠
如は、方法上の2つの主要な相違点が原因になってい
る。第1に、三角形型あるいは他の形状のメンバーシッ
プ関数は、その領域の重み付き面積を説明しているのに
対しては、単集合はその領域に1つの離散的出力を定義
するだけである。ファジィ論理出力を引き出すために複
数の三角形型メンバーシップ出力を組合せる場合、オー
バーラップする面積は重み付きの出力に2度カウントさ
れない。図12の例に見られる通り、メンバーシップ関
数の形状330、340のオーバーラップは、重心の値
350を計算するのに1回カウントされるだけである。
しかし、単集合は、面積の無い出力メンバーシップの重
みを加算することによって累積されるため、隣接メンバ
ーシップ関数の間のどんなオーバーラップも考慮されて
いない。
【0061】第2の相違点は、個々のメンバーシップ出
力が全出力に対して持っている影響に由来する。三角形
型/台形型メンバーシップに対して、メンバーシップの
強度が増加すると、出力の影響(あるいは面積)は、単
集合の出力が線形で増加するようには増加しない。メン
バーシップのルール強度が小さいところでは三角形型メ
ンバーシップの影響は急激にに増加するが、メンバーシ
ップのルール強度が1に近づくのに伴い減少する。
【0062】影響値による方法(Influence Value Meth
od) 面積の中心の出力値を決定するため本発明により考案さ
れた方法は、一方で単集合による方法の簡潔さを利用し
ながら、前に説明した三角形型メンバーシップによる方
法の特性を考慮している。この方法は、各出力メンバー
シップ関数および隣接メンバーシップ関数間のオーバラ
ップした面積に対して、記号をつけかつスケーリングを
行った代表的影響値(IV)を計算することにより、重
心用に面積中心を決定することを包含している。影響値
による方法は、(1)2つの満足している出力ルール間
の潜在的なオーバラップ、(2)ルール強度とその相対
面積との間の非線形関係および(3)ルール強度の関数
としてのルールの面積中心における変動の可能性を考慮
に入れることにより、区分求積法を使用せずに正確な結
果をつくり出している。さらに、事前に計算されたルッ
クアップテーブルを使った影響値による方法を実現する
ことにより、信号処理中の乗算の必要性を大幅に減少も
しくは除去することができ、さらにプロセッサに対して
重負荷の計算量を軽減している。
【0063】添付図面の図13は2つの交差する出力関
数の間のオーバーラップする面積の影響を示している。
図13に見られる通り、2つの三角形型出力メンバーシ
ップ関数501、502がオーバーラップしている場
合、出力重心用の有効面積は(図12により前に説明し
た通り)従来の区分求積法により決定できるが、この方
法は(505に見られる)交差する面積を2度カウント
しない。単集合による方法が使用される場合、交差面積
は、出力コマンド信号を決定するときに有効かつ不適切
に2度カウントされる。本発明によれば、添付図面の図
14に示す通り、交差する面積は負の影響値510によ
って表され、この負の影響値は、図13に見られる交差
する面積505を表している。図13に見られる通り、
負の値510はーバーラップ部分505の面積の中心に
位置にあり、それぞれ出力メンバーシップ関数の面積5
01、502に対応しかつ面積501、502を表して
いる正の影響値512、514の間にある。
【0064】出力メンバーシップは、初めから形状を持
ったメンバーシップ関数として定義されている。本発明
の原理は予め定められた如何なる図形にも応用が可能で
あるが、以下の例における考察を単純化するために、こ
れらの初めに定義されたメンバーシップ関数は、標準台
形型およびベル型メンバーシップ関数を含む三角形をし
ていると想定する。メンバーシップ関数の形状が定義さ
れると、その形状の面積中心が決定できる。図形が面積
中心に関して対象であれば、COAの値はルール強度に
独立であり、さもなければ、以下に考察する通り、CO
Aの値の決定はルール強度の関数としてアルゴリズムに
よって決定されるか、以下に考察するように現在のルー
ル強度値に基づくテーブルルックアップ動作の一部とし
て決定される。
【0065】図15は、2つの交差する三角形型メンバ
ーシップ関数のオーバーラップする部分に対する面積中
心(COA)の位置をアルゴリズムによって決定する方
法を示している。決定されたオーバーラップする部分の
COAは、交差三角形375の底辺の中点と、2つのメ
ンバーシップ関数の面積501、502の交差点、52
1の値との間の中間にある。したがって、得られた値5
30は、関数IntCOA=(Midbase+Int
point)/2によって決定できる。ここにIntC
OAは交差面積の面積中心530、Midbaseは交
差三角形の中点520、Intpointは2つのメン
バーシップ関数の交差点の出力領域の値521である。
このオーバーラップする面積中心値を決定する方法は、
交差する出力メンバーシップ関数が三角形型であり、面
積中心はルール強度に独立であることを想定している。
多くの場合、非三角形型関数が使用される場合でも、出
力メンバーシップ関数の相対面積に比較してオーバーラ
ップの相対面積が小さい場合、この近似により十分な精
度になる。
【0066】交差面積を表す負のIVの大きさを計算す
るためには、オーバーラップする2つのメンバーシップ
関数のルール強度が比較され、ITFを使用して小さい
ほうがIVを計算するために使用され、その小さいほう
のルール強度は、オーバーラップ面積を表す負の影響値
に変換される。以下に考察する通り、この負の影響値の
大きさと出力信号領域における位置は、後続する処理の
中で使用され、満足する出力メンバーシップ関数の全て
に対して希望する面積中心(COA)をより急速に発生
させる。そしてこのCOAが非ファジィ化された最終の
出力コマンド値として使用される。
【0067】図12に示す通り、ルール強度値は出力メ
ンバーシップを上限値まで切り取る効果を持っている。
したがって、三角形型出力メンバーシップ関数は台形型
に切り取られる。この結果、出力メンバーシップ関数の
相対面積は、その関数に適用可能な現在のルール強度の
非線形関数になる。台形になった(shaped)メンバーシ
ップのルール強度と, 最終コマンドに及ばす影響を持つ
ルール強度との間の非線形関係を補償するため、この効
果を説明する正と負のIVの大きさをさらに修正するこ
とが望ましい。添付図面の図16は、(1〜32の整数
として表され、0〜1.0のルール強度の領域をディジ
タル表示するために使用される)スケーリングされたル
ール強度と実効面積との間の関係を示している。図16
に見られる通り、ルール強度が増加するのに伴って、出
力関数の影響(面積)は小さなルール強度に対して急速
に増加するが、ルール強度が最大に近づくのに伴いこの
増加は小さくなる。ルール強度に正比例する大きさを持
つ単集合はエラーを発生するが、このエラーは影響伝達
関数(IFT)を使用することにより回避することがで
きる。
【0068】影響伝達関数(IFT)は、各出力メンバ
ーシップ関数と各オーバーラップ面積とに対して、相対
面積を示す影響値と、影響値掛けるその面積中心値の積
を示すモーメントアーム値とを発生するため、ルックア
ップテーブルにより部分的にあるいは全体的に便利に実
現される。
【0069】希望するならば、ルール強度値を、ルール
強度とその出力コマンド値に関する影響との間の非線形
関係を反映する「影響強度値」を格納するルックアップ
テーブルに索引をつけるために使用しても良い。つぎに
そのルール強度は、1つの記号付き影響値(IV)を求
めるため、面積の形状とサイズとによって値が決定され
る回数だけ乗算され、この記号付き影響値は、区分求積
を必要とせずに全面積(total-area value)の値に直接
加算される。これに代わる方法として、IV値自体がル
ックアップテーブルに格納され、従来の区分求積法を使
用する場合に実行されるような、各相対面積の値を求め
る乗算の必要性が除かれている。
【0070】同様に、モーメントアームも、(上に示す
ように求められた)IVを面積中心の値と掛ける乗算を
行わずに、テーブルルックアップにより直接つくり出す
ことができる。ルール強度値により索引がつけられたモ
ーメントアームの値を直接格納するテーブルを使用する
ことにより、面積中心がルール強度の関数である場合で
も、1回のテーブルルックアップを行うだけで、正確な
モーメントアームの値を得ることができる。
【0071】これまでに考察したとおり、影響伝達関数
(ITV)は、各々独自のメンバーシップ関数の形状と
サイズに対して定義されることが望ましい。つぎに影響
伝達関数は、ゼロでない各出力ルール強度を、面積に比
例しかつ形状が変換された影響値(IV)とモーメント
アームの値に写像するために使用される。ルールの影響
値は、出力信号値に及ばす個々のルールの影響を決定す
る。必要とする変換は、個別のルックアップテーブルを
使用して実現されることが望ましいが、この個別のルッ
クアップテーブルは、各メンバーシップ関数の形状とサ
イズに対して事前に計算され、(たとえばROMの中
に)格納されている。メンバーシップ関数が同じ様な形
状と同じ様な規模を持っている限り、これらの関数は同
一ITFテーブルを共用することができる。
【0072】本発明を使用した好適な非ファジィ化メカ
ニズムは、図17の流れ図に示されている。簡潔に云う
と、この流れ図に示されている方法は次のように要約で
きるであろう。すなわち、 1.出力メンバーシップは、各メンバーシップの面積の
中心を表す位置のスカラー値に変換される。 2.(η個の起こりうる離散的なルール強度値の領域を
表すη個の値を含む)ITFルックアップテーブルを形
成するデータは、各々異なるメンバーシップ関数の形状
とサイズに対する影響伝達関数(ITF)を参照するた
めに格納されている。 3.図9に関連して説明されたように計算される出力メ
ンバーシップのルール強度は、得られたITF変換メカ
ニズムを使用してIVの大きさに変換される。各ルール
強度の大きさはITFによって変換され、そのルールの
影響値IVになる。2つの隣接ルール強度の小さいほう
が、交差面積のルール強度として使用される。正のIV
のCOAの位置は、メンバーシップ関数の形状に関して
事前に計算された面積中心であるから、交差メンバーシ
ップに対するCOA位置の値も同様に事前に計算できる
(たとえば、図15に示す通り、三角形型オーバーラッ
プ面積に対して、COAの位置はオーバーラップ三角形
の底辺の中点と頂点との間の中間の点で極めて良く近似
させることができる)。IVの値の決定は、それぞれ独
自な面積のサイズと形状に対して使用される個別のテー
ブルを備え、ルール強度値によって索引がつけられたI
V値を格納するルックアップテーブルにより実現される
ことが望ましい。 5.全面積の値は、このように決定された影響値を加算
することによりつくられ、全モーメントアームの値は増
分モーメントアーム値を加算することによりつくられる
が、各面積は、引き出されたルールの影響値にそのCO
A位置の値を乗算することによりつくられるか、あるい
はルール強度の値によって索引がつけられたモーメント
アームの値を直接格納するテーブルを使用して、直接テ
ーブルルックアップすることによりつくり出される。交
差部分を表す負のIVは、全面積の値に加算(すなわち
減算)され、交差部分に対する負のモーメントアーム値
も同様に全モーメントアーム値に加算(すなわち減算)
され、これによって隣接メンバーシップ関数の形状の間
のオーバーラップを説明している。 6.非ファジィ化された最終の面積中心(COA)の値
は全モーメントアーム値を全面積の値で除算することに
より求められる。
【0073】図17の流れ図は、好適な非ファジィ化メ
カニズムの動作をより詳細に示している。処理は、処理
変数、MArm、Area、RuleNo、全モーメン
トアーム、全面積、現在のルール番号を表す整数をゼロ
にすることにより、601で開始する。RuleNoは
現在処理されているルールを識別し、その初期値の1か
ら最後のルールの値NoOfRules(定数)迄の範
囲の値を想定する。したがって、603でRuleNo
の初期値ゼロが1に増分された後、処理はRuleNo
=1により開始する。
【0074】RuleNoによって指定されたRule
Strengthは、605でテストされ、ゼロである
か否かが決定される。RuleStrength=0で
あれば、そのルールは出力信号の大きさに何も寄与しな
いので、607で示されるように、そのルールのこれ以
上の処理はバイパスされる。
【0075】613におけるテストで示されるように、
隣接ルールの強度Adj_RuleStrengthが
ゼロであれば、616で示されるように、オーバーラッ
プする交差した面積の処理はバイパスされる。ルールの
出力メンバーシップ関数のCOA値の順番になるように
ルールが配列されているため、「隣接ルール」が、量
(RuleNo−1)によって識別されることが望まし
い。
【0076】625でRuleStrengthとAd
j_RuleStrengthの小さいほうがIntR
uleStrengthとして選択されると、このルー
ル強度が交差部分に割当てられる。つぎにこの値は63
0で、ITFを使用して変換され、影響値IVがつくら
れる。妥当な近似として、オーバーラップ面積の形状の
全てを変換するために三角形型をベースにした1つのI
TFルックアップテーブルが使用される。正のメンバー
シップ関数に対しては、使用されるITFは、多分に変
化する形状と出力メンバーシップ関数の相対的な大きさ
とをベースにして、ルールごとに変わる。したがって、
635では、2次元のITFルックアップテーブルが使
用されるので、この変換の出力IVは、RuleNoと
RuleStrengthの値に依存して決まる。
【0077】640で示す通り、オーバーラップ面積の
負のIV値は、全面積変数Areaの大きさを減少する
ために使用され、全モーメントアーム変数MArmの大
きさを減少するために使用される負の積の値をつくるた
め、オーバーラップ面積の面積中心の値IntCOAに
乗算される。面積中心の値IntCOAは、図13に関
連して説明された関係を使用することにより事前に計算
され、各オーバーラップ面積に対して格納されている。
これに代わる方法として、MArmの値が、ルール強度
により索引がつけられたルックアップテーブルに格納さ
れている(IV*IntCOA[RuleNo])と等
価な交差部分のモーメントアーム値だけ増分される。こ
の場合IntCOAの値は非対象なオーバーラップ面積
のルール強度を変数とする関数(variable function )
でもよい。
【0078】交差部分の処理が完了すると、RuleN
oに対する出力メンバーシップ関数の影響値IVが63
5で決定され、この結果は650で、全面積の値Are
aを増分するために使用される。また、650で、この
IVは、そのRuleNoに対して事前に計算され格納
されている面積中心の値COAに乗算され、得られた積
は、そのルールのモーメントアームの値を示しており、
つぎに全モーメントアームの値MArmに加算される。
これに代わる方法として、そのルールのモーメントアー
ムの値がルール強度によって索引がつけられたモーメン
トアーム値ルックアップテーブルにより直接決定できる
のであり、このことは乗算の必要性を除くとともに、非
対象出力メンバーシップ関数を収容することができる。
【0079】660におけるテストによって示される通
り、評価されたルールが最終ルールでなければ、制御は
603に戻り、次のルールを評価するためにRuleN
oを1だけ増分する。RuleNoがNoOfRule
sに等しい場合が示すように全てのルールが処理される
と、690で示す通りMArmをAreaで除算するこ
とにより、出力重心の最終値COAが決定される。
【0080】図17に示すように、本発明により考案さ
れた非ファジィ化方法は、たとえば図1に見られるファ
ジィ論理制御モジュールを実現するマイクロコントロー
ラによって極めて短時間で実行することができる。図1
7に示すような7つのルールと7つのメンバーシップ関
数を備えているシステムにおいては、この方法は(多く
ても)13の乗算を必要とする。すなわち、(650で
の)7つの乗算は7つの出力メンバーシップの面積のモ
ーメントアームの値を決定するために、そして(640
での)6つの乗算は6つの交差部分のモーメントアーム
の値を決定するために必要になる。最終の重心の計算の
ために(690で)1つの除算が実行される。したがっ
て、1乗算当たり3マイクロセコンドを必要とし、1除
算当たり4マイクロセコンドを必要とするマイクロコン
トローラを使用すると、非ファジィ化に要する時間はほ
ぼ43マイクロセコンドになる。モーメントアームの値
をつくるためにルックアップテーブルが使用されるとす
れば、必要なテーブルのために大きくしたROM容量の
コストで、処理がより速くなるために、この13の乗算
も無くなることになる。
【0081】IV法を実現する人は、ここに説明したI
V法についてのコスト、有用性および正確さに関連する
いくつかの要因が存在することを、念頭に入れておかな
ければならない。第1の要因は、整数を見つける変換を
実行するためルックアップテーブルを使用することに関
係しており、このテーブルルックアップはルール強度の
値がより精緻に定義されている場合量子化エラーにな
る。第2の要因は、実際の面積の値と面積中心の値と
が、想定されている形状に実際に一致していない場合
に、特定の関数の形状を想定したITFアルゴリズムあ
るいはルックアップテーブルを使用することに関係して
いる。想定した形状と、きめの粗い近似は出力信号の正
確さに影響するが、説明済みの実施例に対して実行され
た数多くの評価テストは、IV法を使用して発生した出
力値と、より精密な区分求積法によって発生した値との
間で、最大偏差が2パーセント以下であることが期待で
きることを示している。したがって、この偏差は、ルー
ル強度を32の整数領域で表しているために予想される
量子化エラーに比例した精度レベルになる。
【0082】ここに説明した本発明の実施例は、本発明
の諸原理の一応用の説明にすぎない。本発明の真の精神
と範囲から逸脱せずに、説明した方法および装置に対し
て多数の修正を施すことができるであろう。
【図面の簡単な説明】
【図1】本発明を使用するのに適したファジィ論理コン
トローラの機能構成を示すブロック図。
【図2】信号値をファジィ集合のメンバーシップ度の値
に変換するために使用される三角形型メンバーシップ関
数の図。
【図3】与えられた1つの信号に対する7個の三角形型
メンバーシップ関数を定義するために必要な値を指定す
るために必要な値のメモリマップを示す図。
【図4】代表的テーブルルックアップメンバーシップ関
数用のメモリに格納される値の相対位置を示す図。
【図5】与えられた信号に関連する7つの集合のそれぞ
れに対するルックアップメンバーシップ関数を定義する
値のメモリマップを示す図。
【図6】形状関数用のメモリに格納される値の位置を示
す図。
【図7】与えられた信号に対する7つの形状のメンバー
シップ関数を指定するために必要な値のメモリマップを
示す図。
【図8】指定されたグループの入力値から非ファジィ化
された出力値を決定するルールの集まりを処理中に、メ
ンバーシップ関数の集合を指定するために使用されるア
ドレス指定メカニズムのメモリマップを示す図。
【図9】出力制御信号の大きさをつくる、ファジィ論理
のルールの集つまりの適用可能性を評価するために使用
される、処理ステップのフローチャートを示す図。
【図10】2つのルールによって定義されかつこれらの
2つのルールは、2つの入力信号について処理され1つ
の出力制御値を生じるように定義された、三角形型の入
力メンバーシップ関数および出力メンバーシップ関数を
示す図。
【図11】希望する出力信号の大きさをつくるため、区
分求積法を使用して重心の決定を実行するときに使用さ
れる処理ステップのフローチャートを示す図。
【図12】重心の値を決定する区分求積法を示す図9
(g)の拡大図。
【図13】隣接しかつ共通部分のある2つの出力メンバ
ーシップ関数を示す図。
【図14】図13で見る2つのメンバーシップ関数とそ
れらの共通部分を示すために使用される3つの離散的な
影響値を示す図。
【図15】隣接するメンバーシップ関数の間の共通部部
の面積中心の値を決定する方法を示す図。
【図16】ルール強度を影響値に翻訳するために使用さ
れる型の影響伝達関数を示す図。
【図17】本発明の原理を使用して重心の決定を実行す
るときに使用される処理ステップのフローチャートを示
す図。
【符号の説明】
101 ファジィ論理モジュール 103、105、107 入力信号源 109 前処理(部分) 110 ファジィ論理入力(部分) 112 非ファジィ化されたデータの出力(部分) 113 後処理(部分) 115 アクチュエータ(その他の出力デバイス) 117 ルール格納(用メモリ) 118 出力メンバーシップ関数(用メモリ) 119 入力メンバーシップ関数(用メモリ) 120 ルール強度の決定(部分) 122 ルール強度(用一時的メモリ) 124 重心の決定(部分) 201 ルールのアドレスアレー 202 第2のルールのアドレスアレー 203 メモリデータアレー 204 説明中のデータがあるセル 205 三角形型メンバーシップ関数の形を表すデータ
アレー 301、303、305、307 メンバーシップ度の
値 330、340 メンバーシップ関数の形 350 重心の値 360 区分求積法における区分面積 501、502 メンバーシップ関数の面積 505 2つのメンバーシップ関数の交差面積 510 501、502に対応する負の影響値 512、514 501、502を表す正の影響値 520 交差三角形の中点 521 2つのメンバーシップ関数の交差点の出力領域
の値 530 交差面積の面積中心/重心

Claims (3)

    【特許請求の範囲】
  1. 【請求項1】 異なる信号源からの複数の入力信号の変
    換に適したファジィ論理コントローラであって、 それぞれのルールが、1つの出力メンバーシップ関数
    と、前記入力信号に関連している1つまたはそれ以上の
    入力メンバーシップ関数との組合わせを示すデータを含
    む、複数のルールを格納する手段と、 当該ルールの入力メンバーシップ関数が満足している場
    合、前記ルールのそれぞれの出力メンバーシップ関数か
    ら、前記出力信号に対する推奨値を表す面積中心の値を
    引き出す手段と、 前記ルールのそれぞれの出力メンバーシップ関数から、
    少なくとも1つの入力メンバーシップ関数が満足する程
    度を示すルール強度の値を、そのルール強度の値が前記
    出力コマンド信号の値を変える程度を表す影響値に変換
    する、影響伝達関数を引き出す手段と、 それぞれのルール強度値が、与えられたルールの入力メ
    ンバーシップ関数が1つまたはそれ以上の前記入力信号
    の現在値によって満足する程度を示している、複数のル
    ール強度値の現在値をつくり出すために、前記ルールの
    それぞれの入力メンバーシップ関数に関連する前記入力
    信号の現在値を評価する手段と、 そのルールの影響値の現在値を引き出すため、前記ルー
    ルのそれぞれに関連する影響伝達関数をそのルールの前
    記ルール強度の現在値に適用する手段と、 隣接出力メンバーシップ関数間でオーバーラップする領
    域を表す交差部分の面積中心の値をつくる手段と、 前記オーバーラップする領域のそれぞれを代表する負の
    影響値であって、少なくとも1つの前記隣接出力メンバ
    ーシップ関数のルール強度から引き出される前記負の影
    響値をつくる手段と、 前記負の影響値の和だけ減算された、前記ルールの全影
    響値の現在値を加算することによって全面積の値をつく
    る手段と、 各ルールの影響値の現在値に、そのルールの面積中心の
    値を乗算することによりつくられた積を加算することに
    より、全モーメントアームの値をつくり、そして前記全
    モーメントアームの値を、前記負の影響値に、その関連
    する交差部分の面積中心の値を乗算することによりつく
    られた積の和だけ減少する手段と、 前記全モーメントアームの値を前記全面積の値で割算す
    ることにより前記出力コマンド信号をつくる手段と、を
    含むことを特徴とするファジィ論理コントローラ。
  2. 【請求項2】 請求項1記載のファジィ論理コントロー
    ラであって、前記影響伝達関数を前記ルール強度値に適
    用する手段は、前記ルールの前記影響値をつくり出すた
    めに前記ルール強度値によって索引がつけられた第1の
    テーブルルックアップ関数を含むことを特徴とするファ
    ジィ論理コントローラ。
  3. 【請求項3】 請求項2記載のファジィ論理コントロー
    ラであって、各ルール強度値によってつくられた影響値
    に前記ルール強度値に対応する面積中心の値を乗算する
    ことによりつくられた前記積を格納するため、前記ルー
    ル強度値によって索引がつけられた第2のルックアップ
    テーブルを、さらに含むことを特徴とするファジィ論理
    コントローラ。
JP6280711A 1993-12-27 1994-11-15 ファジィ論理出力の非ファジィ化方法 Pending JPH07219605A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US17299193A 1993-12-27 1993-12-27
US172991 1993-12-27

Publications (1)

Publication Number Publication Date
JPH07219605A true JPH07219605A (ja) 1995-08-18

Family

ID=22630044

Family Applications (1)

Application Number Title Priority Date Filing Date
JP6280711A Pending JPH07219605A (ja) 1993-12-27 1994-11-15 ファジィ論理出力の非ファジィ化方法

Country Status (4)

Country Link
US (1) US5751908A (ja)
JP (1) JPH07219605A (ja)
DE (1) DE4439225B4 (ja)
GB (1) GB2285146B (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19544517C2 (de) * 1995-11-29 1999-08-05 Siemens Ag Steuereinrichtung für ein automatisches Kraftfahrzeuggetriebe
IT1283672B1 (it) * 1996-08-02 1998-04-23 Sgs Thomson Microelectronics Processore elettronico operante in logica fuzzy
DE19640635C2 (de) * 1996-10-02 1999-02-11 Harro Prof Dr Kiendl Verfahren zur Defuzzifizierung für signalverarbeitende Fuzzy-Baueinheiten
US6094646A (en) * 1997-08-18 2000-07-25 Siemens Aktiengesellschaft Method and apparatus for computerized design of fuzzy logic rules from training data vectors of a training dataset
US6205439B1 (en) * 1998-07-15 2001-03-20 Micron Technology, Inc. Optimization of simulation run-times based on fuzzy-controlled input values
US6763337B1 (en) * 1999-03-05 2004-07-13 Lockheed Martin Corporation Weighted wedge defuzzification for conceptual system design evaluation
US6662140B2 (en) 2001-06-12 2003-12-09 Rolls-Royce Canada Limited Fuzzy logic estimator for minimizing signal measurement inaccuracy effects in a signal processing system
GB0321305D0 (en) * 2003-09-11 2003-10-15 Univ Reading The Controlling an environment's characteristics using occupant feedback
DE102004021801A1 (de) * 2004-05-03 2005-12-01 Zf Friedrichshafen Ag Verfahren zur Steuerung eines automatisierten Schaltgetriebes
WO2012061669A2 (en) 2010-11-05 2012-05-10 Cytognomix,Inc. Centromere detector and method for determining radiation exposure from chromosome abnormalities
CN105867114B (zh) * 2016-04-21 2019-02-01 北京航空航天大学 一种新型的电动加载系统控制方法
US20200065692A1 (en) * 2016-12-05 2020-02-27 British Telecommunications Public Limited Company Defuzzification apparatus and method
WO2023101574A1 (en) * 2021-12-03 2023-06-08 Limited Liability Company Solar Security Method and system for static analysis of binary executable code
DE102022120513B4 (de) 2022-08-15 2024-04-25 Krohne S.A.S. Messgerät und Verfahren zum Betreiben eines Messgeräts
CN117648216B (zh) * 2024-01-30 2024-04-23 苏州元脑智能科技有限公司 调整内存可纠正错误阈值的方法、装置、服务器及介质
CN119211317A (zh) * 2024-11-26 2024-12-27 深圳深安阳光电子有限公司 基于物联网的水阀运行控制方法、装置和计算机设备

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5249258A (en) * 1988-09-30 1993-09-28 Omron Tateisi Electronics Co. Reasoning computer system
US5247472A (en) * 1989-02-09 1993-09-21 G.D.S. Co., Ltd. Fuzzy logic operation circuit utilizing charge coupled devices
JPH02208787A (ja) * 1989-02-09 1990-08-20 Yasuo Nagazumi フアジイ演算回路および該回路を用いたファジイ計算機
JPH0771094B2 (ja) * 1989-05-19 1995-07-31 オムロン株式会社 通信ネットワークシステム
US5185849A (en) * 1989-06-27 1993-02-09 Olympus Optical Co., Ltd. Digital fuzzy inference apparatus having time divisional control function
US5148089A (en) * 1989-08-31 1992-09-15 Canon Kabushiki Kaisha Method and apparatus for controlling motors
US5261036A (en) * 1989-10-24 1993-11-09 Mitsubishi Denki K.K. Programmable controller with fuzzy control function, fuzzy control process and fuzzy control monitoring process
JP2656637B2 (ja) * 1989-11-22 1997-09-24 株式会社日立製作所 プロセス制御システム及び発電プラントプロセス制御システム
JPH048964A (ja) * 1990-04-27 1992-01-13 Honda Motor Co Ltd 自動変速機の制御装置
US5267144A (en) * 1990-06-18 1993-11-30 Matsushita Electric Industrial Co., Inc. Access method of an actuator and control apparatus therefor
US5401949A (en) * 1991-06-12 1995-03-28 American Neurologix, Inc. Fuzzy logic barcode reader
US5245695A (en) * 1991-06-12 1993-09-14 American Neuralogix Inc. Fuzzy microcontroller
JPH07503803A (ja) * 1992-02-12 1995-04-20 シーメンス アクチエンゲゼルシヤフト ファジイコントローラ
DE4215406A1 (de) * 1992-05-11 1993-11-18 Zahnradfabrik Friedrichshafen Steuersystem zum Schalten eines automatischen Getriebes
US5371832A (en) * 1992-06-12 1994-12-06 Siemens Aktiengesellschaft Fuzzy logic controller having high processing speed
US5309485A (en) * 1992-07-06 1994-05-03 General Electric Company Core automated monitoring system
DE4230756A1 (de) * 1992-09-15 1993-11-04 Daimler Benz Ag Durch speicherprogrammierbare steuerungen ausfuehrbares steuerungsverfahren nach den regeln der fuzzy-logik
US5261247A (en) * 1993-02-09 1993-11-16 Whirlpool Corporation Fuzzy logic apparatus control
DE69415733T2 (de) * 1993-07-06 1999-05-20 Ford Motor Co., Dearborn, Mich. Elektronisches Steuerungssystem

Also Published As

Publication number Publication date
DE4439225B4 (de) 2005-01-27
GB2285146A (en) 1995-06-28
GB9425155D0 (en) 1995-02-08
DE4439225A1 (de) 1995-06-29
GB2285146B (en) 1998-01-28
US5751908A (en) 1998-05-12

Similar Documents

Publication Publication Date Title
JPH07219605A (ja) ファジィ論理出力の非ファジィ化方法
US5025499A (en) Process control method and control system
Smith et al. Automated calibration of a fuzzy logic controller using a cell state space algorithm
US5491775A (en) Microcontroller fuzzy logic processing module using selectable membership function shapes
JP3076399B2 (ja) ファジィ推論ルールの自動生成装置
Hoffmann Incremental tuning of fuzzy controllers by means of evolution strategy
CA2102944C (en) Method of evaluation of a set of linguistic rules
US5526467A (en) Membership function data preparation method and an apparatus thereof and an adaption degree operation method and an apparatus thereof
Fathi-Torbaghan et al. Model-free optimization of fuzzy rule-based systems using evolution strategies
JP2696972B2 (ja) フアジイ・ルール発生装置および方法,ならびに確認装置および方法
US5988844A (en) Method for arranging a fuzzy controller or a system of interlinked fuzzy controllers
Cordón et al. A practical study on the implementation of fuzzy logic controllers
Grass et al. Programming with anytime algorithms
Barriga et al. Automatic synthesis of fuzzy logic controllers.
Gaines Training, stability and control
CN121447652B (zh) 基于强化学习的人形机器人运动控制模型训练方法、人形机器人运动控制方法、装置、设备、人形机器人及存储介质
JP2532610B2 (ja) プログラム制御の目標値生成装置
JPH02140802A (ja) 自己修正型ファジィ制御方法
JPH02308303A (ja) 自動調整制御装置
Ilg et al. Continuous Q-Learning Resource Allocation Network
JPH04302328A (ja) ファジィ処理装置および方法
JPH04293078A (ja) 地形モデル用3次元データの作成方法
JPH06332708A (ja) ファジィ推論における重心演算装置
JPH0740763A (ja) 制御システム及び定速走行装置
JPH06348496A (ja) メンバシップ関数発生装置及び該メンバシップ関数発生装置を用いたファジィ推論装置