[go: up one dir, main page]

JPH0721024A - 推論装置 - Google Patents

推論装置

Info

Publication number
JPH0721024A
JPH0721024A JP5165654A JP16565493A JPH0721024A JP H0721024 A JPH0721024 A JP H0721024A JP 5165654 A JP5165654 A JP 5165654A JP 16565493 A JP16565493 A JP 16565493A JP H0721024 A JPH0721024 A JP H0721024A
Authority
JP
Japan
Prior art keywords
degree
cause
certainty
certainty factor
calculated
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
JP5165654A
Other languages
English (en)
Inventor
Hiroyoshi Yamaguchi
博義 山口
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.)
Komatsu Ltd
Original Assignee
Komatsu Ltd
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 Komatsu Ltd filed Critical Komatsu Ltd
Priority to JP5165654A priority Critical patent/JPH0721024A/ja
Publication of JPH0721024A publication Critical patent/JPH0721024A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Testing And Monitoring For Control Systems (AREA)
  • Feedback Control In General (AREA)

Abstract

(57)【要約】 【目的】駆動機械の故障等の原因の推論を精度よく行
う。 【構成】予測発生度合いデータに基づいて各種現象の予
測発生度合いが、各種原因の各確信度ごとに演算され
る。そして、この演算された各種原因の各確信度ごとの
予測発生度合いと発生度合いデータとが比較され、各種
原因の各確信度ごとに誤差が演算される。そしてこの演
算された誤差が各種原因ごとに比較され、各種原因にお
いてそれぞれ最も小さい誤差が得られる確信度が各種原
因における確信度であると推論される。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、いわゆる逆診断の診断
結果に基づいて駆動機械の故障等の原因を推論するよう
にした推論装置に関する。
【0002】
【従来の技術】従来より建設機械等の駆動機械の故障原
因を推論するために図2に示すような逆診断知識マトリ
ックスを用いた故障診断装置がある。
【0003】逆診断は、各種原因と各種現象間との因果
関係、つまり原因から予測した各種現象の標準発生パタ
ーンを示す逆診断知識と実際に生じた各種現象の発生パ
ターンとを比較することで故障原因を推定するものであ
る。この逆診断は、標準の発生パターンと実際の発生パ
ターンとの比較の結果、類似しているパターンとなる原
因ほど確信度が高く推定される傾向がある。
【0004】
【発明が解決しようとする課題】しかし、通常、逆診断
に用いられる標準発生パターンは、原因の確信度が10
0%であると仮定して得られたものであり、原因の確信
度が0〜100%で変化するような問題を推論するため
には、そのまま適用することができない。
【0005】また、原因の確信度が変化する場合を考慮
した逆診断としては、法政大学寺野教授が提案したファ
ジイ逆診断があるが(システムと制御vol.24、N
o.11「あいまい論理を用いた異常診断」)、この逆
診断は、複雑なマトリックス演算を実施する必要があ
り、多大な計算時間(2時間を要する場合もある)を要
し、迅速な診断が要求される故障診断には利用すること
ができないこととなっている。
【0006】本発明はこうした実状に鑑みてなされたも
のであり、原因の確信度が0〜100%で変化するよう
な問題を処理する場合に、より誤差の少ない推論結果が
得られるようにするとともに、複雑な演算処理を要する
ことなく短時間に推論を行うことができる推論装置を提
供することを目的としている。
【0007】
【課題を解決するための手段】そこで、この発明では、
各種現象の現在の発生度合いを示す発生度合いデータと
原因から予測した各種現象の発生度合いを示す予測発生
度合いデータとを用意し、前記予測発生度合いデータと
前記発生度合いデータとを各種原因ごとに比較し誤差の
小さい原因ほど確信度が大きくなるよう各種原因ごとの
確信度を推論する推論装置において、前記予測発生度合
いデータに基づいて各種現象の予測発生度合いを、各種
原因の各確信度ごとに演算する第1の演算手段と、前記
第1の演算手段で演算された各種原因の各確信度ごとの
予測発生度合いと前記発生度合いデータとを比較し、各
種原因の各確信度ごとに誤差を演算する第2の演算手段
と、前記第2の演算手段で演算された誤差について各種
原因ごとに比較を行い、各種原因においてそれぞれ最も
小さい誤差が得られる確信度を各種原因におけるそれぞ
れの確信度であると推論する手段とを具えている。
【0008】
【作用】かかる構成によれば、予測発生度合いデータに
基づいて各種現象の予測発生度合いが、各種原因の各確
信度ごとに演算される。そして、この演算された各種原
因の各確信度ごとの予測発生度合いと前記発生度合いデ
ータとが比較され、各種原因の各確信度ごとに誤差が演
算される。そしてこの演算された誤差が各種原因ごとに
比較され、各種原因においてそれぞれ最も小さい誤差が
得られる確信度が各種原因における確信度であると推論
される。
【0009】
【実施例】以下、図面を参照して本発明に係る推論装置
の実施例について説明する。
【0010】実施例では、建設機械の故障診断装置を想
定している。具体的にはパーソナルコンピュータ(以下
パソコンという)が図1に示す手順にしたがい建設機械
の故障診断を行うものである。図2は、逆診断処理を行
うための逆診断知識の構造をそれぞれマトリックス状で
示している。このマトリックスに示される現象の種類は
現象1、現象2、現象3と3種類であり、たとえば図4
に示すように現象1は「排気色が黒い」を、現象2は
「力がない」を、現象3は「徐々に悪化した」をそれぞ
れその内容とする。またマトリックスに示される原因の
種類も原因1、原因2、原因3と3種類であり、たとえ
ば原因1は「エアクリーナのつまり」を、原因2は「噴
射ノズルの不良」を、原因3は「バッテリーの不良」を
それぞれその内容とする。まず、最初にオペレータは、
上記パソコンのキーボードを操作する等して、上記各種
現象1、2、3の現在の発生度合い、つまり異常度Ii
(iは上記各種現象1、2、3に対応する数字である)
を示すデータを入力する。
【0011】ただし、異常度Iiは、0≦Ii≦1の範
囲内の値であり、オペレータの感覚的な数値として与え
る。この場合、異常度Ii=0は「異常が全く発生しな
い」という感覚に対応し、異常度Ii=1は「異常が大
きく発生している」という感覚に対応している。また、
異常が不明(あるいは未チェック)の場合は、上記数値
外を示す’−’をデータとして入力する。ただし、異常
度Ii=’−’の場合には、以下の演算は行われない。
【0012】図5中に異常度Iiを示すデータの具体例
を示しており、説明の便宜のためa、b、cで上記異常
度の数値範囲0〜1を区分している(a>b>c)。た
とえば、No1の例は、現象1「排気色が黒い」の異常
度が大きく(I1=a)、現象2「力がない」の異常度
が中くらいで(I2=b)、現象3「徐々に悪化した」
の異常度が小さい(I3=c)という場合である(ステ
ップ101)。
【0013】つぎに上記入力された異常度Iiに基づい
て逆診断処理が実行される。
【0014】逆診断処理を行うには、各種現象と各種原
因との間の因果関係を示す逆診断知識が必要である。こ
こに因果関係とは、所定の原因jを真の原因と仮定した
ときの所定の現象iの発生確率あるいは発生程度であ
る。この逆診断知識は予め所定のメモリに記憶されてい
るものとする。
【0015】逆診断知識は現象iと原因jの関連度Wi
jとして、0≦Wij≦1の範囲内のあいまいな数値と
して与えられる。この場合、関連度Wij=0は「原因
jを真の原因と仮定したとき現象iは全く発生しない」
に対応し、関連度Wij=1は「原因jを真の原因と仮
定したとき現象iは顕著に発生する」に対応している。
また、現象iと原因jの因果関係が不明の場合は、’
−’により上記数値範囲0〜1外であることを表すよう
にする。ただし、関連度Wij=’−’の場合には、以
下の演算は行われない。
【0016】図4は関連度Wijを示すデータの具体例
を示しており、説明の便宜のためa、b、cで上記関連
度Wijの数値範囲0〜1を区分している(a>b>
c)。たとえば、原因1「エアクリーナのつまり」と現
象1「排気色が黒い」の関連度は中程度であり(W11=
b)、原因1「エアクリーナのつまり」と現象2「力が
ない」の関連度は小さく(W21=c)、原因2「噴射ノ
ズルの不良」と現象1「排気色が黒い」の関連度は大き
い(W12=a)といったごとくである。つぎに、上記関
連度Wijに基づいて予測異常度Pij(k)が演算さ
れる。
【0017】ここに、予測異常度Pij(k)とは、原
因jの確信度がk(0≦k≦1)であると仮定したとき
に予測される各種現象iの発生度合いのことである。確
信度kは0〜1の範囲で連続的な値をとりうるが、演算
処理の便宜上、数値範囲0〜1を10分割あるいは10
0分割して演算される。分割数は、要求される演算精
度、演算時間を勘案して決定される。
【0018】予測異常度Pij(k)はたとえば下記
(1)式より求めることができる。この(1)式の演算
は図3(a)のように図示することができる。
【0019】 Pij(k)={min(Wij、k)} …(1) (ステップ102) ついで、実際の異常度Iiと上記予測異常度Pij
(k)との誤差Eij(k)、つまり原因jの確信度が
kである場合の異常度の誤差Eij(k)が下記(2)
式により演算される。
【0020】 Eij(k)=Ii−Pij(k) …(2) ただし、異常度Ii=’−’または予測異常度Pij
(k)=’−’(Wij=’−’)の場合は、誤差Ei
j(k)=0として扱われる。図3(b)は確信度kと
誤差Eij(k)との関係を例示したグラフである(ス
テップ103)。つぎに誤差Eij(k)の大きさを各
種原因jごとに示すべく、ERRORjが下記(3)式
のごとく演算される。
【0021】 ERRORj(k)=Σ(Eij(k))*2/Nj …(3) ここで、Njは、Iiが’−’でなく、かつWijが’
−’ではないという条件を満たす原因jに対応する現象
iの数である。なお、(X)*2はXの2乗を表すものと
定義する(ステップ104)。
【0022】つぎにERRORj(k)を最小にする確
信度kが原因jごとに求められ、この最小の確信度kが
原因jの最終的な確信度CFjとされる。
【0023】この実施例では、図3(c)に示すように
ERRORj(k)の分布から最小のERRORj
(k)、つまりminERRORj(k)が得られる確
信度kの範囲を求め、この範囲の中で最小の確信度を下
限値とし、また最大の確信度を上限値とし、これら下限
値と上限値の中心値を最小の確信度kとしている。
【0024】これは、中心値近傍の確信度で原因jが発
生していると仮定するのが、最も誤差を小さくし、判断
が妥当と考えられるからである。
【0025】ただし、中心値は、上限値と下限値の差幅
分だけあいまいさを含んでいるため、下記(4)式のよ
うな補正演算を行い、最終確信度CFjを求めるように
する。
【0026】 CFj=中心値−c1・minERRORj−c2・(上限値−下限値) …(4) ただし、c1、c2は定数であり、CFjは、常に0≦C
Fj≦1の範囲内に入るように(4)式の演算後、補正
される(ステップ105)。そして、確信度CFjがパ
ソコンの表示画面に表示される(ステップ106)。
【0027】図5中に逆診断処理による各種原因jの確
信度CFjを百分率で示している。たとえば、No1の
場合、原因1「エアクリーナのつまり」の確信度CF1
は66%となり、原因2「噴射ノズルの不良」の確信度
CF2は90%となり、原因3「バッテリーの不良」の
確信度は30%になることを示している。
【0028】以上説明したようにこの実施例によれば、
実際の異常度Iiと予測異常度Pij(k)との誤差E
ij(k)が最も小さくなるkが確信度CFjとして求
められる。このため、原因の確信度が0〜100%で変
化するような問題を推論する場合に、推論を精度よく行
うことができるようになる。しかも、上記(1)〜
(4)式に示されるように複雑な演算処理を要すること
がないので、短時間に推論を行うことができる。
【0029】
【発明の効果】以上説明したように本発明によれば、原
因の確信度が0〜100%で変化するような問題を推論
する場合に、推論を精度よく行うことができるようにな
る。しかも複雑な演算処理を要することがないので、短
時間に推論を行うことができるようになる。
【図面の簡単な説明】
【図1】図1は本発明に係る推論装置の実施例における
処理手順を示すフローチャートである。
【図2】図2は逆診断知識の構造を概念的に示す図であ
る。
【図3】図3は実施例の演算処理を説明するために用い
たグラフである。
【図4】図4は逆診断知識の具体例を示す図である。
【図5】図5は実施例における診断結果を示す図であ
る。

Claims (1)

    【特許請求の範囲】
  1. 【請求項1】 各種現象の現在の発生度合いを示す
    発生度合いデータと原因から予測した各種現象の発生度
    合いを示す予測発生度合いデータとを用意し、前記予測
    発生度合いデータと前記発生度合いデータとを各種原因
    ごとに比較し誤差の小さい原因ほど確信度が大きくなる
    よう各種原因ごとの確信度を推論する推論装置におい
    て、 前記予測発生度合いデータに基づいて各種現象の予測発
    生度合いを、各種原因の各確信度ごとに演算する第1の
    演算手段と、 前記第1の演算手段で演算された各種原因の各確信度ご
    との予測発生度合いと前記発生度合いデータとを比較
    し、各種原因の各確信度ごとに誤差を演算する第2の演
    算手段と、 前記第2の演算手段で演算された誤差について各種原因
    ごとに比較を行い、各種原因においてそれぞれ最も小さ
    い誤差が得られる確信度を各種原因におけるそれぞれの
    確信度であると推論する手段とを具えた推論装置。
JP5165654A 1993-07-05 1993-07-05 推論装置 Pending JPH0721024A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP5165654A JPH0721024A (ja) 1993-07-05 1993-07-05 推論装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP5165654A JPH0721024A (ja) 1993-07-05 1993-07-05 推論装置

Publications (1)

Publication Number Publication Date
JPH0721024A true JPH0721024A (ja) 1995-01-24

Family

ID=15816466

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5165654A Pending JPH0721024A (ja) 1993-07-05 1993-07-05 推論装置

Country Status (1)

Country Link
JP (1) JPH0721024A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007328645A (ja) * 2006-06-09 2007-12-20 Fuji Xerox Co Ltd 故障診断システム、画像形成装置および故障診断プログラム
WO2011132524A1 (ja) * 2010-04-20 2011-10-27 株式会社日立製作所 異常検知・診断方法、異常検知・診断システム、及び異常検知・診断プログラム

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007328645A (ja) * 2006-06-09 2007-12-20 Fuji Xerox Co Ltd 故障診断システム、画像形成装置および故障診断プログラム
WO2011132524A1 (ja) * 2010-04-20 2011-10-27 株式会社日立製作所 異常検知・診断方法、異常検知・診断システム、及び異常検知・診断プログラム
JP2011227706A (ja) * 2010-04-20 2011-11-10 Hitachi Ltd 異常検知・診断方法、異常検知・診断システム、及び異常検知・診断プログラム

Similar Documents

Publication Publication Date Title
US11593648B2 (en) Methods and systems for detection and isolation of bias in predictive models
Blockley Predicting the likelihood of structural accidents.
US10991053B2 (en) Long-term healthcare cost predictions using future trajectories and machine learning
EP2410312A1 (en) A method for computer-assisted analyzing of a technical system
JP2023052477A (ja) データ処理方法およびデータ処理プログラム
WO2023065584A1 (zh) 设备的剩余使用寿命的确定方法、装置和电子设备
EP4068031A1 (en) Health assessment method for an aircraft system and apparatus therefor
JP7180692B2 (ja) 推定装置、推定方法、及びプログラム
CN104598984A (zh) 一种基于模糊神经网络的故障预测方法
JPH0721024A (ja) 推論装置
JPH0721029A (ja) 推論装置
CN111816291A (zh) 一种设备维护方法及装置
CN119962654B (zh) 一种因果置信规则库的复杂装备故障诊断方法
Han et al. Cascaded LSTM-Based State Prediction of Equipment in Wastewater Treatment Process
WO2023042301A1 (ja) 情報処理装置、情報処理方法、及び、記録媒体
US20090037155A1 (en) Machine condition monitoring using a flexible monitoring framework
Yin et al. A new model based on belief rule base and membership function (BRB-MF) for health state prediction in sensor
JPH0721023A (ja) 推論装置
WO1996019768A1 (fr) Appareil de deduction
Shenoy Valuation network representation and solution of asymmetric decision problems
JP2019028619A (ja) アラーム品質評価装置、アラーム品質評価方法及びアラーム品質評価プログラム
US11669082B2 (en) Online fault localization in industrial processes without utilizing a dynamic system model
CN113408858A (zh) 操作班组风险预警方法、装置、计算机设备和存储介质
JPH02161567A (ja) Ft図利用設備故障診断装置
US5774629A (en) Inference apparatus using occurence rate data and degree of relevance data