JP6558861B2 - 生存確率推定装置、方法、及びプログラム - Google Patents
生存確率推定装置、方法、及びプログラム Download PDFInfo
- Publication number
- JP6558861B2 JP6558861B2 JP2016132060A JP2016132060A JP6558861B2 JP 6558861 B2 JP6558861 B2 JP 6558861B2 JP 2016132060 A JP2016132060 A JP 2016132060A JP 2016132060 A JP2016132060 A JP 2016132060A JP 6558861 B2 JP6558861 B2 JP 6558861B2
- Authority
- JP
- Japan
- Prior art keywords
- probability
- event
- occurrence
- time
- elapsed time
- 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.)
- Active
Links
Landscapes
- Testing And Monitoring For Control Systems (AREA)
- Complex Calculations (AREA)
Description
・故障発生が確認されているか否か(ci)
・機械が設置(製造)されてから故障していないことが確認された最後の時刻までのイベント未発生経過時間(Ti,N)
・機械が設置(製造)されてから故障が実際に発生した時刻までのイベント発生経過時間(Ti,R)
・機械が設置(製造)されてから故障か否かの情報が最後に観測された時刻までの最終観測経過時間(Ti)
・(1,3,5,7)
・(1,6,7,8)
・(1,2,3,5)
・(1,10,12,13)
・(0,13,−,15)
・(0,4,−,8)
20 演算部
30 初期化処理部
32 発生時刻推定部
34 パラメータ推定部
36 収束判定部
50 出力部
100 生存確率推定装置
Claims (7)
- イベントが発生するまで生存する対象物の生存確率を推定する生存確率推定装置であって、
前記対象物の各々について観測されたイベントに関するデータに基づいて、i番目のデータの各々に関して、観測が開始されてからイベントが実際に発生した時刻までのイベント発生経過時間Ti,Rを初期化し、経過時間tにおいてイベントが発生する確率f(t)及び経過時間tまでにイベントが発生していない確率S(t)を推定するためのモデルに関するパラメータを初期化する初期化処理部と、
前記イベントの発生が観測されたi番目のデータの各々に関して、観測が開始されてからイベントが発生していないことが最後に確認された時刻までのイベント未発生観測経過時間Ti,N における、故障が発生していない確率S(T i,N )と、観測が開始されてからイベントが発生したか否かが最後に観測された時刻までの最終観測経過時間Ti における、故障が発生していない確率S(T i )とを、初期化された前記パラメータに基づいて計算し、計算した当該確率S(T i,N )から当該確率S(T i )までの範囲から、一様乱数を用いて確率S(T i,R )を求め、逆関数を用いることにより、前記イベント発生経過時間Ti,Rを推定する発生時刻推定部と、
i番目のデータの各々についての前記イベント発生経過時間Ti,R 及び前記最終観測経過時間T i により表される尤度式を最適化するように、前記パラメータを推定するパラメータ推定部と、
前記発生時刻推定部による推定、及び前記パラメータ推定部による推定を、予め定めた条件を満たすまで繰り返す収束判定処理部と、
を含み、
前記繰り返しにおいて、前記発生時刻推定部は、前記イベントの発生が観測されたi番目のデータの各々に関して、前記確率S(T i,N )と、前記確率S(T i )とを、前記収束判定処理部による繰り返しの一つ前の処理の前記パラメータ推定部によって推定された前記パラメータに基づいて計算し、計算した当該確率S(T i,N )から当該確率S(T i )までの範囲から、一様乱数を用いて確率S(T i,R )を求め、逆関数を用いることにより、前記イベント発生経過時間T i,R を推定する生存確率推定装置。 - 前記パラメータ推定部は、前記モデルとして、指数分布のモデルを用いる請求項1に記載の生存確率推定装置。
- 前記発生時刻推定部により推定された前記イベント発生経過時間Ti,Rを出力する出力部を更に含む請求項1又は請求項2に記載の生存確率推定装置。
- 前記パラメータ推定部は、前記イベントの発生が観測されたi番目のデータの各々に関する、前記イベント発生経過時間Ti,Rにおける、イベントが発生する確率f(Ti,R)と、前記イベントの発生が観測されていないi番目のデータの各々に関する前記確率S(Ti)とを用いて表される前記尤度式を最適化するように、前記パラメータを推定する請求項1〜請求項3のいずれか1項に記載の生存確率推定装置。
- イベントが発生するまで生存する対象物の生存確率を推定する生存確率推定装置における生存確率推定方法であって、
初期化処理部が、前記対象物の各々について観測されたイベントに関するデータに基づいて、i番目のデータの各々に関して、観測が開始されてからイベントが実際に発生した時刻までのイベント発生経過時間Ti,Rを初期化し、経過時間tにおいてイベントが発生する確率f(t)及び経過時間tまでにイベントが発生していない確率S(t)を推定するためのモデルに関するパラメータを初期化するステップと、
発生時刻推定部が、前記イベントの発生が観測されたi番目のデータの各々に関して、観測が開始されてからイベントが発生していないことが最後に確認された時刻までのイベント未発生観測経過時間Ti,N における、故障が発生していない確率S(T i,N )と、観測が開始されてからイベントが発生したか否かが最後に観測された時刻までの最終観測経過時間Ti における、故障が発生していない確率S(T i )とを、初期化された前記パラメータに基づいて計算し、計算した当該確率S(T i,N )から当該確率S(T i )までの範囲から、一様乱数を用いて確率S(T i,R )を求め、逆関数を用いることにより、前記イベント発生経過時間Ti,Rを推定するステップと、
パラメータ推定部が、i番目のデータの各々についての前記イベント発生経過時間Ti,R 及び前記最終観測経過時間T i により表される尤度式を最適化するように、前記パラメータを推定するステップと、
収束判定処理部が、前記発生時刻推定部による推定、及び前記パラメータ推定部による推定を、予め定めた条件を満たすまで繰り返すステップと、を含み、
前記繰り返しにおいて、前記発生時刻推定部が推定するステップは、前記イベントの発生が観測されたi番目のデータの各々に関して、前記確率S(T i,N )と、前記確率S(T i )とを、前記収束判定処理部による繰り返しの一つ前の処理の前記パラメータ推定部によって推定された前記パラメータに基づいて計算し、計算した当該確率S(T i,N )から当該確率S(T i )までの範囲から、一様乱数を用いて確率S(T i,R )を求め、逆関数を用いることにより、前記イベント発生経過時間T i,R を推定する生存確率推定方法。 - 出力部が、前記発生時刻推定部により推定された前記イベント発生経過時間Ti,Rを出力するステップを更に含む請求項5に記載の生存確率推定方法。
- コンピュータを、請求項1〜請求項4のいずれか1項に記載の生存確率推定装置の各部として機能させるためのプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016132060A JP6558861B2 (ja) | 2016-07-01 | 2016-07-01 | 生存確率推定装置、方法、及びプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016132060A JP6558861B2 (ja) | 2016-07-01 | 2016-07-01 | 生存確率推定装置、方法、及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018005580A JP2018005580A (ja) | 2018-01-11 |
JP6558861B2 true JP6558861B2 (ja) | 2019-08-14 |
Family
ID=60949305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016132060A Active JP6558861B2 (ja) | 2016-07-01 | 2016-07-01 | 生存確率推定装置、方法、及びプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6558861B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6903246B2 (ja) * | 2019-03-04 | 2021-07-14 | 三菱電機株式会社 | 情報処理装置、情報処理方法及び情報処理プログラム |
-
2016
- 2016-07-01 JP JP2016132060A patent/JP6558861B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2018005580A (ja) | 2018-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108700852B (zh) | 模型参数值推定装置及推定方法、记录介质、模型参数值推定*** | |
Echard et al. | A combined importance sampling and kriging reliability method for small failure probabilities with time-demanding numerical models | |
CN104823191B (zh) | 用于概率性疲劳裂纹寿命估计的方法和*** | |
Wang et al. | Bayesian dynamic forecasting of structural strain response using structural health monitoring data | |
JP5261740B2 (ja) | コンクリート構造物の中性化深さ予測装置および中性化深さをコンピュータに計算させるためのプログラム | |
Gobbato et al. | A recursive Bayesian approach for fatigue damage prognosis: An experimental validation at the reliability component level | |
Roy et al. | An NHPP software reliability growth model with imperfect debugging and error generation | |
JP6610987B2 (ja) | 異常診断方法及び異常診断システム | |
EP3979120A2 (en) | Probabilistic fatigue and blend limit assessment and visualization methods for airfoils | |
KR102240646B1 (ko) | 수명 예측 방법 및 시스템 | |
JP6845126B2 (ja) | 故障確率算出装置、故障確率算出方法及びプログラム | |
KR20140139622A (ko) | Eifs 불확실성을 고려하여 초음파 검사 데이터를 사용한 확률적 피로 수명 예측 | |
US11361233B2 (en) | Estimating fatigue life of technical systems | |
JP2015529813A (ja) | 遺伝的プログラミングを用いて発見された前兆的特徴からの、残存耐用寿命の推定 | |
Blesa et al. | Set‐membership parity space approach for fault detection in linear uncertain dynamic systems | |
JP2020051385A (ja) | 熱機器の内部状態推定方法及び装置 | |
Bansal | A new Gibbs sampling based Bayesian model updating approach using modal data from multiple setups | |
CN111160666B (zh) | 强噪声与非周期状态监测的健康状态与可靠性评估方法 | |
Proverbio et al. | Adaptive sampling methodology for structural identification using radial basis functions | |
WO2020105468A1 (ja) | 情報処理装置、情報処理システム、情報処理方法及びプログラムが格納された非一時的なコンピュータ可読媒体 | |
JP6558861B2 (ja) | 生存確率推定装置、方法、及びプログラム | |
JP6558862B2 (ja) | 生存確率推定装置、方法、及びプログラム | |
Sun et al. | Reliability Analysis with Multiple Dependent Features from a Vibration‐Based Accelerated Degradation Test | |
Li et al. | Global Sensitivity Analysis Using an Auxiliary Variable Approach | |
JP2008165412A (ja) | 性能演算装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180817 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190507 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20190426 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190702 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20190709 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190715 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6558861 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |