状態判定方法と状態予測方法及び装置
対象物の状態監視のために測定された波形データや特徴パラメータが正規分布に従わなくても、対象物の状態診断と状態予測は統計検定などの手法により高精度に行わ、また、測定された対象物の信号が脈動信号である場合、信号の局所に存在する特異成分はリアルタイムに検出され、対象物の状態は効率的に判定される必要がある。本発明は、測定した波形データ、あるいは、波形データから算出された特徴パラメータを既知確率分布（例えば、正規分布）に従うように変換した後、統計検定などの手法により対象物の状態判定と状態予測を行う方法と装置、また、測定した対象物の信号が脈動信号である場合は、雑音を除去された信号の包絡線波形データなどを求めることにより信号中の特異成分を検出して対象物の状態判定を行う方法と装置を提供する。
【発明の詳細な説明】
一．技術分野
本発明は設備診断、医療診断などにおいて、対象物の状態変化の有無を判定するための状態判定装置及びオンライン状態監視・診断システムに関するものである。
二．背景技術
（１）従来、状態監視の対象物から測定した波形データの確率密度関数は、正常状態時に「正規分布」に従うと仮定して状態判定を行う（［１］特許公開２０００−１７１２９１、［２］豊田利夫、二保知也：「正常時の振動波形のみを用いた回転機械の異常診断」、日本設備管理学会誌、Ｖｏｌ．１１，Ｎｏ．１，１９９９年、ｐ．４−１１．［３］劉 信芳、豊田利夫、陳 鵬、馮 芳、二保知也：「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅによる回転機械の異常診断」、精密工学会誌、Ｖｏｌ．６６、Ｎｏ．１、２０００年、ｐ．１５７−１６２．）
（２）対象物の状態変化の有無を判定するために用いる特徴パラメータを正規分布の確率変数と仮定して統計検定を行う（［４］河部佳樹、豊田利夫、江口透、福田和久：「無次元兆候パラメータを用いた回転機械の劣化傾向管理（ＩＩＩ）」、平成７年度日本設備管理学会秋季研究発表大会論文集、Ｖｏｌ．２、１９９５年、ｐ．３２．）
状態監視の対象物から測定した信号が脈動信号の場合は、特徴パラメータや確率密度分布関数などにより状態判定が困難である。従来、脈動信号の中に含まれる特異成分（あるいは、異常成分）を検出する方法として、
（３）時間−周波数解析手法（ウェーブレット解析、ウィグナ分布解析、短時間ＦＦＴなど）を用いた局所異常の検出法がある（［５］章 忠、中堀智之、川畑洋昭：「高速ウェーブレット変換およびその脳波解析への応用」、日本機械学会論文集（Ｃ編）、Ｖｏｌ．６５，Ｎｏ．６３２，１９９９年、ｐ．１９１５−１９２１．）
（４）脈動信号のピーク値などに合わせた所定の間隔でサンプリングした値は所定の値以下、または、前回サンプリングした値との差が所定の値以上の場合に異常と判定する。（［６］特開平６−３１７２１５）
（５）脈動信号のピーク値などに合わせた所定の間隔でサンプリング毎の代表値のばらつき（標準偏差）が所定のしきい値以上の場合に異常と判定する。（［７］特開平７−２９３３１１）
三．発明が解決しようとする課題
しかしながら、前記従来の諸手法には次のような問題があった。
第１と第２の方法では、波形データと特徴パラメータの値が正規分布に従うと仮定しているが、実際に正常状態でも必ずしも正規分布に従うとは限らないから、波形データと特徴パラメータの値が正規分布に従うと仮定して対象物の状態変化を統計検定により判定すると、誤判定をもたらす。
第３の方法では、処理に必要な時間が長いため、アルタイムの特異点の検出が困難である。
第４と第５の方法では、波形データを所定の間隔でサンプリングする場合、回転数の変化や負荷の変化などによるピーク位置が変化した時、サンプリング値は正常でも大きく変化するから、誤判定をもたらす。また、サンプリング毎の代表値のばらつき（標準偏差）を計算すると、リアルタイムの特異検出が困難である。
四．課題を解決するための手段
上記に述べたような問題点を解決するために、本発明においては、測定した波形データを既知確率分布（例えば、正規分布）の波形データに変換した後、あるいは、波形データから算出された特徴パラメータを既知確率分布（例えば、正規分布）の特徴パラメータに変換した後、統計検定や可能性理論や情報理論などにより対象物の状態判定を行う。
また、測定した対象物の信号が脈動信号である場合は、特徴パラメータやスペクトルなどによる特異成分の抽出及び状態判定が困難であるので、雑音を除去した信号の包絡線波形データを求め、正規化された包絡線波形データおよび周期パルスにより信号中の特異成分を検出して状態判定を行う。または、上記に雑音を除去した脈動信号のピーク波形データを求め、ピーク波形データの正規化処理を行い、正規化されたピーク波形データおよび周期パルス波形データにより脈動信号中の特異成分を検出して状態判定を行う。
五．発明の効果
対象物の状態監視のために測定された波形データや特徴パラメータは必ずしも正規分布に従うと限らないため、これらの波形データあるいは特徴パラメータを正規分布に従うと仮定して対象物の状態変化の有無を統計検定により判定し、或いは、状態予測を行うと、大きな誤差をもたらす。本発明においては、測定した波形データを既知確率分布（例えば、正規分布）の波形データに変換した後、あるいは、波形データから算出された特徴パラメータを既知確率分布（例えば、正規分布）の特徴パラメータに変換した後、統計検定や可能性理論や情報理論などにより対象物の状態判定を行う。従って、本発明の手法は状態判定、或いは、状態予測の精度が従来正規分布に従うと仮定した場合より高い。
脈動信号の場合、測定した脈動信号の波形データには一部のピークが欠落し、または波形データの局所に微小な特異成分が存在する時、従来の判定法（特徴パラメータやスペクトルや確率密度関数など）による検出が困難である。本発明においては、雑音を除去した脈動信号の包絡線波形データを求め、正規化された包絡線波形データおよび周期パルスにより信号中の特異成分を検出して状態判定を行う。または、上記に雑音を除去した脈動信号のピーク波形データを求め、ピーク波形データの正規化処理を行い、正規化されたピーク波形データおよび周期パルス波形データにより脈動信号中の特異成分を検出して状態判定を行う。本発明の方法では、主な信号処理はハードウェアでも実現でき、数値処理装置（計算機）にかかる負担が小さいから、リアルタイムに異常の検出が実現できる。また、正常状態の波形データを基準としないため、脈動周期が変化しても検出と判定の結果に与える影響が少ない。
六．発明を実施するための最良の形態
１．特徴パラメータについて
状態判定のために用いる特徴パラメータは時間領域と周波数領域の特徴パラメータがある。周波数領域の特徴パラメータは（参考文献１）で定義されているものがある。（［８］陳 鵬，豊田利夫：遺伝的プログラミングによる周波数領域の特徴パラメータの自己再組織化，日本機械論文集（Ｃ編），Ｖｏｌ．６５ Ｎｏ．６３３，ｐｐ．１９４６−１９５３，１９９８．）ここで時間領域の特徴パラメータについて詳述する。
対象物の状態変化の有無を判定するために用いる時間領域の特徴パラメータは次の通りである。
１）無次元特徴パラメータ
測定した時系列波形データから、フィルタを用いて、低、中、高周波数領域の波形データを抽出する。抽出した波形データｘ（ｔ）を次式で正規化する。
ここで、ｘ’ｉはＡ／Ｄ変換後のｘ（ｔ）の離散値であり、μとＳはそれぞれｘ’ｉの平均値と標準偏差である。
従来用いられる無次元特徴パラメータは式（２）〜式（１３）に示す（［９］Ｐｅｎｇ ＣＨＥＮ，Ｔｏｓｈｉｏ ＴＯＹＯＴＡ，Ｙｕｅｔｏｎ ＬＩＮ，Ｆｅｉｙｕｅ ＷＡＮＧ：ＦＡＩＬＵＲＥ ＤＩＡＧＮＯＳＩＳ ＯＦ ＭＡＣＨＩＮＥＲＹ ＢＹ ＳＥＬＦ−ＲＥＯＲＧＡＮＩＺＡＴＩＯＮ ＯＦ ＳＹＭＰＴＯＭ ＰＡＲＡＭＥＴＥＲＳ ＩＮ ＴＩＭＥ ＤＯＭＡＩＮ ＵＳＩＮＧ ＧＥＮＥＴＩＣ ＡＬＧＯＲＩＴＨＭＳ，Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｉｎｔｅｌｌｉｇｅｎｔ Ｃｏｎｔｒｏｌ ａｎｄ Ｓｙｓｔｅｍ，Ｖｏｌ．３，Ｎｏ．４，ｐｐ．５７１−５８５，１９９９．）
ここで、
は絶対平均値で，Ｎはデータの総数である。
は標準偏差である。
ここで、μｐは波形の極大値（ピーク値）の平均値である。
ここで、｜μｍａｘ｜は波形の１０個の最大値の平均値である。
ここで、σｐは極大値の標準偏差値である。
ここで、μＬとσＬはそれぞれ極小値（谷値）の平均値と標準偏差値である。
式（２）〜式（１３）は従来の特徴パラメータであるが、数値計算で容易に高速計算するために、「区間特徴パラメータ」は式（１４）〜式（２１）のように新たに提案する。
但し、ｘｉ≧ｋσ、ｋは任意に設定できるが、例えば、ｋ＝０．５、１、２。μｋ１はｘｉの平均値である。ｔは任意に設定できるが、例えば、ｔ＝０．５、１、２、３、４。
但し、ｘｈ≦−ｈσ、ｈは任意に設定できるが、例えば、ｋ＝０．５、１、２。μｈ１はｘｈの平均値である。ｔは任意に設定できるが、例えば、ｔ＝０．５、１、２、３、４。
但し、ｈ０は時系列波形が単位時間あたり０レベルをクロースする頻度、ｈｐは単位時間あたり時系列波形のピークの数である。
但し、ｈｎσは波形が単位時間あたりｎσレベルをクロースする頻度である。ｎは任意に設定できるが、例えば、ｎ＝０．５、１、２。
但し、ｈ−ｎσは波形が単位時間あたり−ｎσレベルをクロースする頻度。である。ｎは任意に設定できるが、例えば、ｎ＝０．５、１、２。
２）有次元特徴パラメータ
有次元特徴パラメータを計算するとき、測定した波形データに対して、式（１）のような正規化をしない。
ここで、｜ｘｉ｜ｐは波形データの絶対値のピーク値（極大値）、Ｎｐはピーク値の総数である。
なお、上記の特徴パラメータ以外にも多くの特徴パラメータが定義できるが、本方法を応用するとき、まず上記の特徴パラメータで試し、もし状態識別の効果が良くなければ、更に他の特徴パラメータを追加定義すればよい。
２．波形データと特徴パラメータを指定の確率分布への変換
測定した波形データをｘ＊ｉで、波形データから算出した特徴パラメータをｐ＊ｉで表す。ｘ＊ｉあるいはｐ＊ｉを用いて、統計理論により状態判定や状態予測を行う場合、ｘ＊ｉあるいはｐ＊ｉがどのような確率分布に従うか事前に知る必要がある。しかし、ｘ＊ｉあるいはｐ＊ｉがどのような確率分布従うかは事前に知らない場合が多い。そこで、指定の既知確率分布関数をΣとすると、次式を用いてｘ＊ｉあるいはｐ＊ｉをΣに従う確率変数ｘｉあるいはｐｉに変換することができる。
ここで、Ｆｘｉ（ｘ＊ｉ）とＦｐｉ（ｐ＊ｉ）はそれぞれｘ＊ｉとｐ＊ｉとの累積確率分布（或いは、累積度数分布）であり、Σ−１はΣの逆関数である。例えば、Σは正規分布、ワイブル分布、指数分布、ガンマ分布などである。ｘ＊ｉあるいはｐ＊ｉをｘｉあるいはｐｉに変換した後、統計検定などを用いて状態判定や状態予測を行う。
なお、元の波形データｘ＊ｉは図１のように４種類に分けられる。つまり、平均値より大きいデータｘｉ＋、平均値より小さいデータｘｉ−、式（１）による正規化後の絶対値データ｜ｘｉ｜及び全体波形データｘ＊ｉＡであり、それぞれの累積確率分布（或いは、累積度数分布）はＦｘｉ＋（ｘｉ＋）、Ｆｘｉ−（ｘｉ−）、Ｆ｜ｘｉ｜（｜ｘｉ｜）及びＦｘｉＡ（ｘ＊ｉＡ）で表すが、以下、特に指定しなければ、統一にＦｘｉ（ｘ＊ｉ）で表す。また、特徴パラメータｐ＊ｉは正規化後のｘｉ＋とｘｉ−、｜ｘｉ｜及びｘ＊ｉＡのいずれかで計算される。
ここで、例として、Σを正規分布に指定した場合、ｘ＊ｉあるいはｐ＊ｉを正規分布に従うように変換する方法について詳述する。
正規分布の確率密度関数ｆ（ｔ）は次式で表す。
ここで、μは確率変数ｔの平均値で、σは標準偏差である。
（１）基準状態を参照する場合
測定対象の基準状態、例えば、１回目測定した時の状態、を決めて、この時の波形データｘ＊ｉ０あるいは特徴パラメータｐ＊ｉ０の確率密度関数と累積確率分布関数をそれぞれｆｘ０（ｘ＊ｉ０）とＦｘ０（ｘ＊ｉ０）あるいはｆｐ０（ｐ＊ｉ０）とＦｐ０（ｐ＊ｉ０）とする。平均値μが０で標準偏差σが１である標準正規分布の確率密度関数をφ（ｘｉ）、標準正規分布の確率分布関数をΦ（ｘｉ）とする。なお、離散データであるｘ＊ｉ０あるいはｐ＊ｉ０の確率密度関数の替わりに「頻度分布関数」や「ヒスとグラム」を用いてもよいが、以下には「確率密度関数」を用いて説明する。
１）正規分布に従う平均値に変換する方法
状態ｋでのｘ＊ｉｋとｐ＊ｉｋを次式でそれぞれ正規分布に従う平均値μｘｉｋｏとμｐｉｋ０に変換される。なお、状態ｋは任意の状態であり、基準状態も含む。
ここで、Φ−１はΦの逆関数で、σｘｉ０とσｐｉ０はそれぞれ正規分布に変換されたｘ＊ｉｋとｐ＊ｉｋとの標準偏差であり、次式で求められる。
平均値μｘｉｋ０とμｐｉｋ０を用いて状態判定や状態予測を行う。
２）直接変換法
状態ｋでのｘ＊ｉｋとｐ＊ｉｋを次式で正規分布の確率変数に変換する。
ここで、ＳｘｋとＳｐｋはそれぞれｘ＊ｉｋとｐ＊ｉｋの標準偏差で、μｘｋとμｐｋはそれぞれｘ＊ｉｋとｐ＊ｉｋの平均値である。ｘ’ｉｋ０あるいはｐ’ｉｋ０を用いて状態判定や状態予測を行う。
（２）基準状態を参照しない場合
まず、状態ｋでの波形データｘ＊ｉｋと特徴パラメータｐ＊ｉｋとの確率密度関数（あるいは、頻度分布）をそれぞれｆｘｋ（ｘ＊ｉｋ）とｆｐｋ（ｐ＊ｉｋ）とし、確率分布関数（あるいは、累積度数分布）をＦｘｋ（ｘ＊ｉｋ）とＦｐｋ（ｐ＊ｉｋ）とする。
１）直接変数変換法
波形データｘ＊ｉｋと特徴パラメータｐ＊ｉｋを次式で正規分布の確率変数に変換する。
ここで、ＳｘｋとＳｐｋはそれぞれｘ＊ｉｋとｐ＊ｉｋとの標準偏差で、μｘｋとμｐｋはｘ＊ｉｋとｐ＊ｉｋとの平均値である。ｘ’ｉｋあるいはｐ’ｉｋを用いて状態判定や状態予測を行う。
２）正規分布に従う平均値に変換する方法
正規分布に従う平均値は次式で求められる。
ここで、σｘｉｋとσｐｉｋは次式で求められる。
μｘｉｋとμｐｉｋは正規分布に従うので、μｘｉｋあるいはμｐｉｋを用いて状態判定や状態予測を行う。
３）間接変換法
Ｍ個の特徴パラメータｐ＊ｉｋの最小値と最大値をそれぞれ（ｐ＊ｉｋ）ｍｉｎと（ｐ＊ｉｋ）ｍａｘとする。（ｐ＊ｉｋ）ｍｉｎから（ｐ＊ｉｋ）ｍａｘまでＮ個の等間隔区間に分割する。各区間の代表値（例えば、中央値）をｐ＊ｉｋｊとする。ここで、ｊ＝１〜Ｎ。ｐ＊ｉｋｊをｐ＊ｉｋの替わりに式（３６）または式（３８）に代入すると、Ｎ個のｐ’’ｉｋまたはμ’ｐｉｋが得られる。ｐ’’ｉｋまたはμ’ｐｉｋを用いて状態判定や状態予測を行う。
同様に、Ｍ個の波形データｘ＊ｉｋの最小値と最大値をそれぞれ（ｘ＊ｉｋ）ｍｉｎと（ｘ＊ｉｋ）ｍａｘとする。（ｘ＊ｉｋ）ｍｉｎから（ｘ＊ｉｋ）ｍａｘまでＮ個の等間隔区間に分割する。各区間の代表値（例えば、中央値）をｘ＊ｉｋｊとする。ここで、ｊ＝１〜Ｋ。ｘ＊ｉｋｊをｘ＊ｉｋの替わりに式（３５）または式（３４）に代入すると、Ｋ個のｘ’’ｉｋまたはμ’ｘｉｋが得られる。ｘ’’ｉｋまたはμ’ｘｉｋを用いて状態判定や状態予測を行う。
４）特徴パラメータｐ＊ｉｋの区間平均値を求める方法
Ｎ個特徴パラメータｐ＊ｉｋを求めた後、Ｍ組に分割して第ｊ組の平均値は次のように求める。
ここで、Ｎｊは第ｊ組にあるｐ＊ｉｋの数である。μｋ（ｊ）は近似的に正規分布に従うので、μｋ（ｊ）を用いて状態判定や状態予測を行う。
波形データｘ＊ｉｋを式（２９）、式（３３）、式（３５）、式（３７）で正規確率分布の波形データに変換したμｘｉｋ０、ｘ’ｉｋ０、ｘ’ｉｋ、μｘｉｋ、ｘ’’ｉｋ、μ’ｘｉｋは「正規分布の波形データ」と呼ぶ。また、特徴パラメータｐ＊ｉを式（３０）、式（３４）、式（３６）、式（３８）で正規確率変数に変換したμｐｉｋ０、ｐ’ｉｋ０、ｐ’ｉｋ、μｐｉｋ、ｐ’’ｉｋ、μ’ｐｉｋ、μｋ（ｊ）は「正規分布の特徴パラメータ」と呼ぶ。
３．正規分布の特徴パラメータによる状態判別方法
ここで、正規分布の特徴パラメータを用いて、対象物の状態を判別する方法を述べる。
（１）統計理論による判別
１）正規分布の特徴パラメータの平均値の検定
状態ｋと状態ｙで求めた正規分布の特徴パラメータをそれぞれｐｉｋとｐｉｙとする。ここで、ｉ＝１〜Ｍ、Ｍは使用する正規分布の特徴パラメータの総数を表す。ｐｉｋとｐｉｙとの平均値をそれぞれμｉｋとμｉｙ、ｐｉｋとｐｉｙとの標準偏差をそれぞれＳｉｋとＳｉｙとする。一般にＪ個のｐｊの平均値μと標準偏差Ｓは次の式で計算する。
μｉｋとμｉｙが等しいか否かの検定は次のように行う（参考文献３参照）
（参考文献３）Ｋ．Ａ．Ｂｒｏｗｎｌｅｅ．Ｓｔａｔｉｓｔｉｃａｌ Ｔｈｅｏｒｙ ａｎｄ Ｍｅｔｈｏｄｏｌｏｇｙ ｉｎ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｅｃｏｎｄ Ｅｄｉｔｉｏｎ，Ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ Ｃｈｉｃａｇｏ，１９６５
が成立すれば、有意水準αで「μｉｋとμｉｙとが等しくない」と判定する。ここで、ｔα／２（Ｊ−１）は自由度Ｊ−１のｔ分布の確率密度関数が下側確率α／２に対するパーセント点である。
２）正規分布の特徴パラメータの分散の検定
ＳｉｋとＳｉｙが同じか否かの検定は次のように行う（参考文献４参照）。
（参考文献４）Ｋ．Ａ．Ｂｒｏｗｎｌｅｅ．Ｓｔａｔｉｓｔｉｃａｌ Ｔｈｅｏｒｙ ａｎｄ Ｍｅｔｈｏｄｏｌｏｇｙ ｉｎ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｅｃｏｎｄ Ｅｄｉｔｉｏｎ，Ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ Ｃｈｉｃａｇｏ，１９６５
が成立すれば、有意水準αで「ＳｉｋとＳｉｙとが等しくない」と判定する。ここで、Ｆα／２（Ｊ−１，Ｊ−１）は自由度Ｊ−１のＦ分布の確率密度関数が下側確率α／２に対するパーセント点である。
有意水準αを変えたとき、式（４４）または式（４５）を満足するか否かを確認することにより、状態ｙが状態ｋに対する状態変化の程度を決める。状態変化の程度を有意水準αにより決定する例は表１に示す。なお、設備診断の場合、状態ｋを正常状態とすれば、状態ｙは正常状態か、注意状態か、危険状態かの判別については、表１のように「正常」（α１）、「注意」（α２）、及び「危険」（α３）のように設定し検定することができる。つまり、式（４４）または式（４５）はα１の時に成立しなければ「正常」と判定する。また、式（４４）または式（４５）はα２の時に成立すれば「注意」と、α３の時に成立すれば「危険」と判定する。なお、表１中のαの数値範囲は例であり、設備の重要度などによって決められる。
数個の特徴パラメータを用いて状態判定を行う場合、判定の結果は最も状態変化が大きいと示した特徴パラメータの判定結果に従われる。例えば、３つの特徴パラメータｐ１、ｐ２、ｐ３を用いて判定するとき、ｐ１の判定結果が「注意」、ｐ２の判定結果が「正常」、ｐ３の判定結果が「危険」である場合、最終の判定結果は「危険」とする。
３）信頼区間による判定
基準時点で測定した波形データから求めた正規分布の特徴パラメータの平均値をμｉ０とし、他の時点で測定した波形データから求めた正規分布の特徴パラメータの平均値をμｉｋとすると、μｉ０の信頼区間は次式で与えられる。
ここで、ｔα／２（Ｊ−１）は自由度Ｊ−１のｔ分布の確率密度関数が下側確率α／２に対するパーセント点である。Ｓｉ０は波形データから求めたｐｉ０の標準偏差である。μｉｋは式（４６）に示す区間内にあれば、１−αの確率でμｉ０との差がないという。μｉ０の９９％信頼区間は、Ｊ＞１０のとき、近似的に次のようになる。
従って、μｉｋは式（４７）の範囲を超えれば、９９％の確率でμｉ０と違うと判定する。更に、測定した波形データで求めたμｉｋの信頼区間は、次式で求められる。
ここで、Ｓｉｋは波形データから求めたｐｉｋの標準偏差である。
表１のαを式（４６）に代入すれば、次のような信頼区間が得られる。
μｉｋがこれらの区間内にあるか否かによって状態の判定を行う。
数個の特徴パラメータを用いて状態判定を行う場合、最終の判定結果は式（４９），（５０），（５１）に記述される内容と同じである。
（２）可能性理論による判別
１）可能性分布関数の作成
状態ｋにおける波形データで正規分布の特徴パラメータｐｉの値を算出した後、ｐｉの確率密度関数ｆｋ（ｐｉ）から可能性分布関数Ｐｋ（ｐｉ）を式（５２）で求める。可能性理論によれば、ｐｉがどのような確率分布に従っても、その可能性分布関数が求められる。ｐｉが正規分布に従う場合、Ｎ段の可能性分布関数ｐｋ（ｐｉ）は次のように求める（参考文献５参照）。
ここで、
但し、上式中、ｐｉｘ＝ｍｉｎ｛ｐｉ｝＋ｘ×（ｍａｘ｛ｐｉ｝−ｍｉｎ｛ｐｉ｝）／Ｎ、ｘ＝１〜Ｎ，Ｓｉはｐｉの標準偏差，μｉはのｐｉ平均値である。
（参考文献５）Ｌ．Ｄａｖｉｓ：ＨＡＮＤＢＯＯＫ ＯＦ ＧＥＮＥＴＩＣ ＡＬＧＯＲＩＴＨＭＳ，Ｖａｎ Ｎｏｓｔｒａｎｄ Ｒｅｉｎｈｏｌｄ，Ａ Ｄｉｖｉｓｉｏｎ ｏｆ Ｗａｄｓｗｏｒｔｈ，Ｉｎｃ（１９９０）
２）可能性の求め方
図２のように、状態ｋと状態ｙで求めた正規分布の特徴パラメータｐｉの可能性分布関数をＰｋ（ｐｉ）とＰｙ（ｐｉ）とし、状態ｙで求めた正規分布の特徴パラメータの値をｐ’ｉとすると、「状態ｙが状態ｋと同じである」という可能性ｗは次のように求められる。
ａ）ｐ’ｉの平均値ｐ’ｉｍｅａｎとＰｋ（ｐｉ）とのマッチングによるｗの決定、
ｂ）Ｐｙ（ｐｉ）とＰｋ（ｐｉ）とのマッチングによるｗの決定。
なお、Ｐｙ（ｐｉ）とＰｋ（ｐｉ）とのマッチングによるｗを求める式は次に示す。
３）状態変化の判別
状態ｋで求めた正規分布の特徴パラメータｐｉの可能性分布関数ｐｋ（ｐｉ）が得られた後、左右両側の「状態変化が小さい」の可能性分布関数（ｐｃ１（ｐｉ）とｐｃ２（ｐｉ））、および「状態変化が大きい」の可能性分布関数（ｐｄ１（ｐｉ）とｐｄ２（ｐｉ））は図２に示すように決定する。境界値の
におけるｉ，ｊはユーザ入力により決定するが、標準値としてｉ＝３、ｊ＝６とする。
設備診断の場合、正常状態の可能性分布関数をｐｋ（ｐｉ）、注意状態の可能性分布関数をｐｃ１（ｐｉ）とｐｃ２（ｐｉ）、危険状態の可能性分布関数をｐｄ１（ｐｉ）とｐｄ２（ｐｉ）とする。実際の識別時に得られた「正常」と「注意」と「危険」との可能性は図２のように表示する。また、「危険」と判定した場合、警報を出すことも可能である。
（３）情報理論による状態判別法
対象物の基準状態における正規分布の特徴パラメータの確率密度関数をｆｐ０（ｐｉ）し、基準状態以外における正規分布の特徴パラメータの確率密度関数をｆｐｋ（ｐｉ）とする。なお、基準状態以外の時点での状態を「テスト状態」と呼ぶ。テスト状態が基準状態と同じ状態であるか否かは次の「Ｋｕｌｌｂａｃｋ−Ｌｅｉｂｌｅｒ Ｉｎｆｏｒｍａｔｉｏｎ（カルバック情報量、ＫＩ）」と「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅ（情報ダイバージェンス、ＩＤ）」で判定できる。
ＫＩｐとＩＤｐによる状態判別法は［１０］に詳細に記述しているので、ここで省略する。（［１０］劉 信芳、豊田利夫、陳 鵬、馮 芳、二保知也：「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅによる回転機械の異常診断」、精密工学会誌、Ｖｏｌ．６６、Ｎｏ．１、２０００年、ｐ．１５７−１６２．）
（４）複数の特徴パラメータを統合することによる状態判別法
数個の特徴パラメータを統合して、状態変化の有無や状態予測を行うこともできる。特徴パラメータの統合法は、主成分分析法やＫＬ展開法などがあるが、特徴パラメータの統合法により求めた新たな特徴パラメータは「統合特徴パラメータ」とよぶ。ここで、主成分分析法の例を示す。（［１１］大津，栗田，関田著：パターン認識，朝倉書店、１９９６．）
例えば、設備診断の場合、設備の正常状態下で求めた無次元特徴パラメータ（ｐ１，ｐ２，・・・，ｐｍ）に対して、ｍ個の主成分は次のように表されます。
各主成分ｚ１〜ｚｍは「統合特徴パラメータ」とも呼ぶ。
相関行列は次のように求められます。
ここで、正常状態で収録したｎ組のデータを｛ｐ１ｋ，ｐ２ｋ，・・・，ｐｍｋ｝、ｋ＝１，２，・・・，ｎとすると、
相関行列Ｒの固有ベクトル、λ＝（λ１，λ２，・・・，λｍ）、λ１，≧λ２≧・・・≧λｍとすると、
固有値λｉに対応する式（５９）の係数が求まり、第ｉ番目の主成分は次のように求められます。
主成分を用いて次式のように状態判定を行う。
ここで、Ｋは使用する主成分の数、αは有意水準、χ２（Ｋ，α）は自由度Ｋのカイ２乗分布の確率密度関数が上側確率αに対するパーセント点である。αの決定法は表１同じである。なお、Ｋ＝３に対してα＝０．０５のとき、χ２（３，０．０５）＝７．８１５である。
なお、式（６４）に示す判定方法以外に、特徴パラメータの統合法により求めた統合特徴パラメータ、例えば、式（６３）に示す主成分、を正規分布の確率変数に変換した後、統計検定や可能性理論により状態判定を行うこともできる。
ここで、実例を示す。図３はある回転機械の正常状態時（図３（ａ））と回転軸ミスアライメント状態時（図３（ｂ））に測定した振動の加速度波形データである。図４（ａ）及び図５（ａ）は式（２）〜式（１３）及び式（１８）〜式（２１）に示す１４個の無次元特徴パラメータを６０回求め、正規分布変換を行う前に、Ｋ＝３のときに求めた式（６４）右側の値（６０個）を示す。図４（ｂ）及び図５（ｂ）は式（３１）に基づいて正規分布変換を行った後、Ｋ＝３のときに求めた式（６４）右側の値（９個）を示す。これらの図により、正規分布変換後の結果は正規分布変換前の結果より良好であることがわかる。なお、「良好である」とは、式（６４）左側の値は正常状態の時に小さく、異常状態の時に大きいということである。
４．正規分布の波形データによる状態判別方法
（１）特徴パラメータによる状態判別法
対象物から計測された波形データｘ＊ｉ（例えば、図６（ａ）、（ｃ））を式（３３）により正規分布の波形データｘ’ｉｋ０（例えば、図６（ｂ）、（ｄ））に変換した後、正規分布の波形データｘ’ｉｋ０を用いて前述の正規分布の特徴パラメータによる状態判別方法で状態判定を行うことができる。なお、この方法は正規分布の波形データμｘｉｋ０とｘ’ｉｋ０に適用できる。
（２）情報理論による状態判別法
基準状態の正規分布の波形データの確率密度関数をｆｘ０（ｘｉ）とし、基準状態以外の時点での正規分布の波形データの確率密度関数をｆｘｋ（ｘｉ）とする。なお、基準状態以外の時点での状態を「テスト状態」と呼ぶ。テスト状態が基準状態と同じ状態であるか否かは次の「Ｋｕｌｌｂａｃｋ−Ｌｅｉｂｌｅｒ Ｉｎｆｏｒｍａｔｉｏｎ（カルバック情報量、ＫＩ）」と「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅ（情報ダイバージェンス、ＩＤ）」で判定できる。
ＫＩｘとＩＤｘによる状態判別法はに詳細に記述しているので、ここで省略する。（［１２］劉 信芳、豊田利夫、陳 鵬、馮 芳、二保知也：「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅによる回転機械の異常診断」、精密工学会誌、Ｖｏｌ．６６、Ｎｏ．１、２０００年、ｐ．１５７−１６２．）
ここで、実例を示す。図６（ａ）はある回転機械の正常状態時に測定した振動加速度の波形データである。図６（ｃ）は同じ回転機械のアンバランス時に測定した波形データである。
平均値以上の波形データｘｉ＋、平均値以下の波形データｘｉ−、及び絶対値データ｜ｘｉ｜を正規分布の波形データに変換した後、それぞれｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’で表す。ｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’の平均値と標準偏差は次の通りである。
ｘ’ｉ＋：正常時の平均値＝０．６９、アンバランス時の平均値＝１．９６
正常時の標準偏差＝１．０３、アンバランス時の標準偏差＝２．３７
ｘ’ｉ−：正常時の平均値＝−０．３７、アンバランス時の平均値＝−０．６６
正常時の標準偏差＝０．４７、アンバランス時の標準偏差＝０．６９
｜ｘｉ｜’：正常時の平均値＝１．４５、アンバランス時の平均値＝３．２９
正常時の標準偏差＝１．５２、アンバランス時の標準偏差＝３．５３
（参考文献８）によると、過検出率α＝０．１５、見逃率β＝０．１５とすると、ＫＩ＞１．２１或いはＩＤ＞２．４３ならば、８５％ｎの確率で「テスト状態は基準状態と違う」と判定できる。上記のｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’のＫＩとＩＤは次の通りである。
ｘ’ｉ＋：ＫＩ＝０．５１９， ＩＤ＝２．２４
ｘ’ｉ−：ＫＩ＝０．５３， ＩＤ＝０．４０
｜ｘｉ｜’：ＫＩ＝０．５７， ＩＤ＝２．６７
｜ｘｉ｜’のＩＤ＞２．４３から、８５％ｎの確率で「図６（ａ）の状態は図６（ｃ）の状態と違う」と判定できる。
また、式（４４）と式（４５）に示す平均値と分散の検定により状態判定を行ってもよい。
更に、上記のｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’を用いて、式（２）〜式（２５）で特徴パラメータを求めて、統計検定や可能性理論及び特徴パラメータの統合法により状態判定を行うこともできる。
５．状態予測
正規分布の特徴パラメータｐｉが求められた後、従来の状態予測手法を用いて、測定対象の状態を予測することが行える（［１３］石川真澄、武藤博道：予測手法、計測と制御、１９８２．３．［１４］Ｏｇａｗａ，Ｍ．：Ｔｉｍｅ ｓｅｒｉｅｓ ａｎａｌｙｓｉｓ ａｎｄ ｓｔｏｃｈａｓｔｉｃ ｐｒｅｄｉｃｔｉｏｎ，Ｂｕｌｌ．Ｍａｔｈ．Ｓｔａｔ．，８，８−７２，１９５８．［１５］Ｂ．ピカニオル著、小村賢二、柴山宮恵子訳：予測のための統計学、晃洋書房、１９８７．）
図７は状態予測の方法を示す。各測定時点（ｘ１〜ｘ７）で測定した特徴パラメータ値、或いは、主成分値を正規確率分布に従うように変換した後、回帰解析により予測曲線とその信頼区間を求め、「寿命限界」との交点で「最短寿命」、「平均寿命」及び「最長寿命」を求める。
図８はある回転機械が正常状態からアンバランス状態へ変化していく間に８回測定した波形例である。測定１は正常状態であり、測定８は最も程度の重いアンバランス状態である。図９は式（２）〜式（１３）で求めた、各測定における特徴パラメータの平均値を示す。この例でわかるように、特徴パラメータは異常状態の程度が重くなっていくにつれて値が単調増加の特徴パラメータ（例えば、ｐ４，ｐ５）があれば、値が単調減少の特徴パラメータ（例えば、ｐ２，ｐ８，ｐ９，ｐ１０）もある。また、異常状態の程度が重くなっても値がほぼ不変な特徴パラメータ（例えば、ｐ１，ｐ７）もある。従って、値が単調増加（或いは、単調減少）の特徴パラメータを選んで、寿命予測を行うべきである。また、主成分により寿命予測を行う時、値が単調増加の特徴パラメータのみ（或いは、単調減少の特徴パラメータのみ）を選んで、主成分を求めて寿命予測を行うべきである。
上記に述べた波形データと特徴パラメータを正規分布の確率変数に変換する手法、状態判定手法及び状態予測手法を実現するための計測と処理の流れは図１６（ａ）に示す。また、図１６（ａ）を実現するための波形データ計測及び状態判定装置の回路は図１７に示す。
６．脈動信号から特異成分の検出法及び状態判定法
ここで、ガスエンジン失火診断の例を用いて説明する。
図１０（ａ）、図１１（ａ）、図１２（ａ）、図１３（ａ）はガスエンジン失火（図中のピーク値の異常個所）が発生したときのシリンダ圧力波形データである。なお、図１２（ａ）、図１３（ａ）は多数の異常ピークが連続して発生する例を示す。このような脈動信号に発生する局所異常（特異成分）に対して、オンライン特異点の検出及び状態判定の流れは図１５に示す。以下、その手順を説明する。
１）．診断ための準備（図１５（ａ））
（１）診断対象の脈動信号と回転パルス信号を同時に計測する。なお、回転パルス信号は周期パルス信号とも言い、脈動信号の各ピーク値のタイミングを決めるときに用いられる。
（２）雑音を除去するために，ローパスフィルタリングを行う。ローパスフィルタのカットオフ周波数ｆＬは次のように決定する。
ここで、ｎは軸の回転数（ｒｐｍ）、ｚは１回転毎のピーク数（ピーク／１回転）、ｆ０は余裕周波数（＞ｎ／６０、フィルタリング後のノイズ除去効果を観察することにより決定する）である。
（３）ローパスフィルタリング後の包絡線波形データ、または、ピーク波形データ、または、移動平均波形データを求める。
（４）上記の包絡線波形データ、または、ピーク波形データ、または、移動平均波形データの正規化を次式のように行う。
ここで、ｘ（ｔ）は正規化後の波形データ、ｘ’（ｔ）は元の波形データ、μ’（ｔ）はｘ’（ｔ）の平均値、ｓはｘ’（ｔ）の標準偏差である。なお、μ’（ｔ）とｓを求めるとき、正常状態の波形データが望ましいが、正常状態の波形データでなくてもよい。また、μ’（ｔ）とｓを求めるとき、運転条件（回転数や負荷）をできるだけ判定のときと同じように設定する。但し、運転条件が多少変化しても判定結果に大きく影響しない。
（５）回転パルス信号のピークと脈動信号のピークとの相対位置関係を調べて置き、異常を識別するための閾値を決定する。閾値の決定は次のように考える。
のとき、異常が発生したと判定する。ここで、ｋ（＝２〜４）は診断対象物によって決定する。
｜ｘ（ｔ）｜がｋσより大きい間、回転パルス信号に対応する脈動信号のピークの所に異常があると判定する。
２）．オンライン診断（図１５（ｂ））
（１）判定対象物の脈動信号と回転パルス信号を同時に計測する。
（２）ローパスフィルタリングを行う。ローパスフィルタのカットオフ周波数ｆＬは式（６６）のように決定する。
（３）ローパスフィルタリング後の包絡線波形データ、または、ピーク波形データを求める。
（４）上記の包絡線波形データ、または、ピーク波形データの正規化は式（６８）により求める。
（５）正規化された包絡線波形データ、または、ピーク波形データの絶対値は閾値（ｋσ）より大きいか否かを監視する。大きくなければ、異常がないと判定する。大きければ、異常が発生したと判定し、｜ｘ（ｔ）｜がｋσより大きい間、前もって測って置いた回転パルス信号と脈動信号のピークとの関係により、異常の箇所を判定する。
上記の手順に基づく信号処理の例を図１０、図１１、図１２、図１３に示す。
図１０（ａ）（ｂ）、図１１（ａ）（ｂ）、図１２（ａ）（ｂ）、図１３（ａ）（ｂ）ガスエンジン失火（図中の異常箇所）が発生したときのシリンダ圧力の時系列波形データとスペクトルである。
図１０（ｃ）（ｄ）、図１１（ｃ）（ｄ）、図１２（ｃ）（ｄ）、図１３（ｃ）（ｄ）はローパスフィルタにより雑音を除去した時系列波形データとスペクトルである。
図１０（ｅ）、図１１（ｅ）、図１２（ｅ）、図１３（ｅ）は回転パルス波形データである。
図１０（ｆ）、図１２（ｆ）は包絡線波形データである。
図１１（ｆ）、図１３（ｆ）はピーク波形データである。
図１０，１１，１２，１３の（ｆ）で分かるように、異常（ピークの欠落、または、小さくなる）が発生した場合、その絶対値｜ｘ（ｔ）｜が２σより大きくなることが判明できる。
また、絶対値｜ｘ（ｔ）｜が２σより大きくなる箇所と回転パルス波形データとの対応関係により異常箇所を判明できる。
上記に包絡線波形データとピーク波形データの例を示したが、図１４に示すように、移動平均波形データ（ｕｉ）を求めたら、ｘ（ｔ）の替わりにｕ（ｔ）を用いて上記の手順と同様に脈動信号の特異点検出ができる。なお、カットオフ周波数（ｆｃ）は式（６７）により求めたら、移動平均の点数（Ｍ）が求められる。
上記に述べた脈動信号の特異成分の検出と判定方法を実現するための計測と処理の流れは図１６（ｂ）に示す。また図１６（ｂ）を実現するための脈動信号計測及び状態判定装置の回路は図１８に示す。
７．状態判定装置、または、オンライン状態判定システムの処理流れ
状態判定装置装置、または、オンライン状態判定システムの処理流れは図１６に示す。図１６（ａ）に示すように、測定波形データからノイズを除去し、特徴パラメータを求め、その特徴パラメータを正規分布の確率変数に変換した後、統計検定や可能性理論や特徴パラメータの統合により状態変化の有無判定を行う。これらの処理は計算機または専用装置で実現できる。また、図１６（ｂ）に示すように、脈動信号からの特異成分検出及び状態判定は、ローパスフィルタ、包絡線（または、ピーク値）はハードウェアで実現できる。数値演算装置（または、計算機）は正規化したｘ（ｔ）、｜ｘ（ｔ）＞ｋσ｜の判定、特異箇所の検出と判定および判定結果の表示を行えばよいから、リアルタイムに行える。
【図面の簡単な説明】
図１は、測定した波形データを４種類、つまり、平均値より大きいデータｘｉ＋、平均値より小さいデータｘｉ−、式（１）による正規化後の絶対値データ｜ｘｉ｜及び全体波形データｘ＊ｉＡに分ける例を示すグラフである。
図２は可能性分布関数の例を示すグラフである。
図３は正常状態と回転軸ミスアライメント状態の振動波形データの例を示すグラフである。
図４は正常状態の振動波形データに対する主成分分析結果を示すグラフである。
図５は回転軸ミスアライメント状態の振動波形データに対する主成分分析結果を示すグラフである。
図６は正常状態とアンバランス状態の振動波形データの例を示すグラフである。
図７は状態予測方法を示すグラフである。
図８は正常状態からアンバランス状態へ変化していく時の振動波形データの例を示すグラフである。
図９は各測定の波形の特徴パラメータ値を示すグラフである。
図１０は包絡線波形データによる１個のピークが欠落した場合の信号処理例を示すグラフである。
図１１はピーク波形データによる１個のピークが欠落した場合の信号処理例を示すグラフである。
図１２は包絡線波形データによる多数のピーク異常の場合の信号処理例を示すグラフである。
図１３はピーク波形データによる多数のピーク異常の場合の信号処理例を示すグラフである。
図１４は移動平均波形データの求め方を示すグラフである。
図１５は脈動信号の特異成分検出及び異常診断処理の流れを示すフローチャートである。
図１６は状態判定装置、または、オンライン状態判定システムの処理流れを示すフローチャートである。
図１７は信号計測及び状態判定装置の回路の一例を示す回路図であり、図中の符号は次の通りである。
１ センサ、２ アンプ、３ フィルタ、４ 信号処理・計算装置、５ 表示出力装置、６ データ用ＲＡＭ、７ ＡＤ変換器、８ ＤＣポート、９ ＳＣＩ（シリアル・コミュニケーション・インタフェース）、１０ ＣＰＵ、１１ フラッシュＲＯＭ、１２ 外部コンピュータ．
図１８は脈動信号計測及び状態判定装置の回路の一例を示す回路図であり、図中の符号は次の通りである。
１ センサ、２ アンプ、３ フィルタ、４ 包絡線処理部、５ 周期パルスセンサ、６ 表示出力装置、７ ＡＤ変換器、８ ＤＣポート、９ ＳＣＩ（シリアル・コミュニケーション・インタフェース）、１０ ＣＰＵ、１１ フラッシュＲＯＭ、１２ 外部コンピュータ、１３ データ用ＲＡＭ、１４ 信号処理・計算装置．
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
一．技術分野
本発明は設備診断、医療診断などにおいて、対象物の状態変化の有無を判定するための状態判定装置及びオンライン状態監視・診断システムに関するものである。
二．背景技術
（１）従来、状態監視の対象物から測定した波形データの確率密度関数は、正常状態時に「正規分布」に従うと仮定して状態判定を行う（［１］特許公開２０００−１７１２９１、［２］豊田利夫、二保知也：「正常時の振動波形のみを用いた回転機械の異常診断」、日本設備管理学会誌、Ｖｏｌ．１１，Ｎｏ．１，１９９９年、ｐ．４−１１．［３］劉 信芳、豊田利夫、陳 鵬、馮 芳、二保知也：「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅによる回転機械の異常診断」、精密工学会誌、Ｖｏｌ．６６、Ｎｏ．１、２０００年、ｐ．１５７−１６２．）
（２）対象物の状態変化の有無を判定するために用いる特徴パラメータを正規分布の確率変数と仮定して統計検定を行う（［４］河部佳樹、豊田利夫、江口透、福田和久：「無次元兆候パラメータを用いた回転機械の劣化傾向管理（ＩＩＩ）」、平成７年度日本設備管理学会秋季研究発表大会論文集、Ｖｏｌ．２、１９９５年、ｐ．３２．）
状態監視の対象物から測定した信号が脈動信号の場合は、特徴パラメータや確率密度分布関数などにより状態判定が困難である。従来、脈動信号の中に含まれる特異成分（あるいは、異常成分）を検出する方法として、
（３）時間−周波数解析手法（ウェーブレット解析、ウィグナ分布解析、短時間ＦＦＴなど）を用いた局所異常の検出法がある（［５］章 忠、中堀智之、川畑洋昭：「高速ウェーブレット変換およびその脳波解析への応用」、日本機械学会論文集（Ｃ編）、Ｖｏｌ．６５，Ｎｏ．６３２，１９９９年、ｐ．１９１５−１９２１．）
（４）脈動信号のピーク値などに合わせた所定の間隔でサンプリングした値は所定の値以下、または、前回サンプリングした値との差が所定の値以上の場合に異常と判定する。（［６］特開平６−３１７２１５）
（５）脈動信号のピーク値などに合わせた所定の間隔でサンプリング毎の代表値のばらつき（標準偏差）が所定のしきい値以上の場合に異常と判定する。（［７］特開平７−２９３３１１）
三．発明が解決しようとする課題
しかしながら、前記従来の諸手法には次のような問題があった。
第１と第２の方法では、波形データと特徴パラメータの値が正規分布に従うと仮定しているが、実際に正常状態でも必ずしも正規分布に従うとは限らないから、波形データと特徴パラメータの値が正規分布に従うと仮定して対象物の状態変化を統計検定により判定すると、誤判定をもたらす。
第３の方法では、処理に必要な時間が長いため、アルタイムの特異点の検出が困難である。
第４と第５の方法では、波形データを所定の間隔でサンプリングする場合、回転数の変化や負荷の変化などによるピーク位置が変化した時、サンプリング値は正常でも大きく変化するから、誤判定をもたらす。また、サンプリング毎の代表値のばらつき（標準偏差）を計算すると、リアルタイムの特異検出が困難である。
四．課題を解決するための手段
上記に述べたような問題点を解決するために、本発明においては、測定した波形データを既知確率分布（例えば、正規分布）の波形データに変換した後、あるいは、波形データから算出された特徴パラメータを既知確率分布（例えば、正規分布）の特徴パラメータに変換した後、統計検定や可能性理論や情報理論などにより対象物の状態判定を行う。
また、測定した対象物の信号が脈動信号である場合は、特徴パラメータやスペクトルなどによる特異成分の抽出及び状態判定が困難であるので、雑音を除去した信号の包絡線波形データを求め、正規化された包絡線波形データおよび周期パルスにより信号中の特異成分を検出して状態判定を行う。または、上記に雑音を除去した脈動信号のピーク波形データを求め、ピーク波形データの正規化処理を行い、正規化されたピーク波形データおよび周期パルス波形データにより脈動信号中の特異成分を検出して状態判定を行う。
五．発明の効果
対象物の状態監視のために測定された波形データや特徴パラメータは必ずしも正規分布に従うと限らないため、これらの波形データあるいは特徴パラメータを正規分布に従うと仮定して対象物の状態変化の有無を統計検定により判定し、或いは、状態予測を行うと、大きな誤差をもたらす。本発明においては、測定した波形データを既知確率分布（例えば、正規分布）の波形データに変換した後、あるいは、波形データから算出された特徴パラメータを既知確率分布（例えば、正規分布）の特徴パラメータに変換した後、統計検定や可能性理論や情報理論などにより対象物の状態判定を行う。従って、本発明の手法は状態判定、或いは、状態予測の精度が従来正規分布に従うと仮定した場合より高い。
脈動信号の場合、測定した脈動信号の波形データには一部のピークが欠落し、または波形データの局所に微小な特異成分が存在する時、従来の判定法（特徴パラメータやスペクトルや確率密度関数など）による検出が困難である。本発明においては、雑音を除去した脈動信号の包絡線波形データを求め、正規化された包絡線波形データおよび周期パルスにより信号中の特異成分を検出して状態判定を行う。または、上記に雑音を除去した脈動信号のピーク波形データを求め、ピーク波形データの正規化処理を行い、正規化されたピーク波形データおよび周期パルス波形データにより脈動信号中の特異成分を検出して状態判定を行う。本発明の方法では、主な信号処理はハードウェアでも実現でき、数値処理装置（計算機）にかかる負担が小さいから、リアルタイムに異常の検出が実現できる。また、正常状態の波形データを基準としないため、脈動周期が変化しても検出と判定の結果に与える影響が少ない。
六．発明を実施するための最良の形態
１．特徴パラメータについて
状態判定のために用いる特徴パラメータは時間領域と周波数領域の特徴パラメータがある。周波数領域の特徴パラメータは（参考文献１）で定義されているものがある。（［８］陳 鵬，豊田利夫：遺伝的プログラミングによる周波数領域の特徴パラメータの自己再組織化，日本機械論文集（Ｃ編），Ｖｏｌ．６５ Ｎｏ．６３３，ｐｐ．１９４６−１９５３，１９９８．）ここで時間領域の特徴パラメータについて詳述する。
対象物の状態変化の有無を判定するために用いる時間領域の特徴パラメータは次の通りである。
１）無次元特徴パラメータ
測定した時系列波形データから、フィルタを用いて、低、中、高周波数領域の波形データを抽出する。抽出した波形データｘ（ｔ）を次式で正規化する。
ここで、ｘ’ｉはＡ／Ｄ変換後のｘ（ｔ）の離散値であり、μとＳはそれぞれｘ’ｉの平均値と標準偏差である。
従来用いられる無次元特徴パラメータは式（２）〜式（１３）に示す（［９］Ｐｅｎｇ ＣＨＥＮ，Ｔｏｓｈｉｏ ＴＯＹＯＴＡ，Ｙｕｅｔｏｎ ＬＩＮ，Ｆｅｉｙｕｅ ＷＡＮＧ：ＦＡＩＬＵＲＥ ＤＩＡＧＮＯＳＩＳ ＯＦ ＭＡＣＨＩＮＥＲＹ ＢＹ ＳＥＬＦ−ＲＥＯＲＧＡＮＩＺＡＴＩＯＮ ＯＦ ＳＹＭＰＴＯＭ ＰＡＲＡＭＥＴＥＲＳ ＩＮ ＴＩＭＥ ＤＯＭＡＩＮ ＵＳＩＮＧ ＧＥＮＥＴＩＣ ＡＬＧＯＲＩＴＨＭＳ，Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｉｎｔｅｌｌｉｇｅｎｔ Ｃｏｎｔｒｏｌ ａｎｄ Ｓｙｓｔｅｍ，Ｖｏｌ．３，Ｎｏ．４，ｐｐ．５７１−５８５，１９９９．）
ここで、
は絶対平均値で，Ｎはデータの総数である。
は標準偏差である。
ここで、μｐは波形の極大値（ピーク値）の平均値である。
ここで、｜μｍａｘ｜は波形の１０個の最大値の平均値である。
ここで、σｐは極大値の標準偏差値である。
ここで、μＬとσＬはそれぞれ極小値（谷値）の平均値と標準偏差値である。
式（２）〜式（１３）は従来の特徴パラメータであるが、数値計算で容易に高速計算するために、「区間特徴パラメータ」は式（１４）〜式（２１）のように新たに提案する。
但し、ｘｉ≧ｋσ、ｋは任意に設定できるが、例えば、ｋ＝０．５、１、２。μｋ１はｘｉの平均値である。ｔは任意に設定できるが、例えば、ｔ＝０．５、１、２、３、４。
但し、ｘｈ≦−ｈσ、ｈは任意に設定できるが、例えば、ｋ＝０．５、１、２。μｈ１はｘｈの平均値である。ｔは任意に設定できるが、例えば、ｔ＝０．５、１、２、３、４。
但し、ｈ０は時系列波形が単位時間あたり０レベルをクロースする頻度、ｈｐは単位時間あたり時系列波形のピークの数である。
但し、ｈｎσは波形が単位時間あたりｎσレベルをクロースする頻度である。ｎは任意に設定できるが、例えば、ｎ＝０．５、１、２。
但し、ｈ−ｎσは波形が単位時間あたり−ｎσレベルをクロースする頻度。である。ｎは任意に設定できるが、例えば、ｎ＝０．５、１、２。
２）有次元特徴パラメータ
有次元特徴パラメータを計算するとき、測定した波形データに対して、式（１）のような正規化をしない。
ここで、｜ｘｉ｜ｐは波形データの絶対値のピーク値（極大値）、Ｎｐはピーク値の総数である。
なお、上記の特徴パラメータ以外にも多くの特徴パラメータが定義できるが、本方法を応用するとき、まず上記の特徴パラメータで試し、もし状態識別の効果が良くなければ、更に他の特徴パラメータを追加定義すればよい。
２．波形データと特徴パラメータを指定の確率分布への変換
測定した波形データをｘ＊ｉで、波形データから算出した特徴パラメータをｐ＊ｉで表す。ｘ＊ｉあるいはｐ＊ｉを用いて、統計理論により状態判定や状態予測を行う場合、ｘ＊ｉあるいはｐ＊ｉがどのような確率分布に従うか事前に知る必要がある。しかし、ｘ＊ｉあるいはｐ＊ｉがどのような確率分布従うかは事前に知らない場合が多い。そこで、指定の既知確率分布関数をΣとすると、次式を用いてｘ＊ｉあるいはｐ＊ｉをΣに従う確率変数ｘｉあるいはｐｉに変換することができる。
ここで、Ｆｘｉ（ｘ＊ｉ）とＦｐｉ（ｐ＊ｉ）はそれぞれｘ＊ｉとｐ＊ｉとの累積確率分布（或いは、累積度数分布）であり、Σ−１はΣの逆関数である。例えば、Σは正規分布、ワイブル分布、指数分布、ガンマ分布などである。ｘ＊ｉあるいはｐ＊ｉをｘｉあるいはｐｉに変換した後、統計検定などを用いて状態判定や状態予測を行う。
なお、元の波形データｘ＊ｉは図１のように４種類に分けられる。つまり、平均値より大きいデータｘｉ＋、平均値より小さいデータｘｉ−、式（１）による正規化後の絶対値データ｜ｘｉ｜及び全体波形データｘ＊ｉＡであり、それぞれの累積確率分布（或いは、累積度数分布）はＦｘｉ＋（ｘｉ＋）、Ｆｘｉ−（ｘｉ−）、Ｆ｜ｘｉ｜（｜ｘｉ｜）及びＦｘｉＡ（ｘ＊ｉＡ）で表すが、以下、特に指定しなければ、統一にＦｘｉ（ｘ＊ｉ）で表す。また、特徴パラメータｐ＊ｉは正規化後のｘｉ＋とｘｉ−、｜ｘｉ｜及びｘ＊ｉＡのいずれかで計算される。
ここで、例として、Σを正規分布に指定した場合、ｘ＊ｉあるいはｐ＊ｉを正規分布に従うように変換する方法について詳述する。
正規分布の確率密度関数ｆ（ｔ）は次式で表す。
ここで、μは確率変数ｔの平均値で、σは標準偏差である。
（１）基準状態を参照する場合
測定対象の基準状態、例えば、１回目測定した時の状態、を決めて、この時の波形データｘ＊ｉ０あるいは特徴パラメータｐ＊ｉ０の確率密度関数と累積確率分布関数をそれぞれｆｘ０（ｘ＊ｉ０）とＦｘ０（ｘ＊ｉ０）あるいはｆｐ０（ｐ＊ｉ０）とＦｐ０（ｐ＊ｉ０）とする。平均値μが０で標準偏差σが１である標準正規分布の確率密度関数をφ（ｘｉ）、標準正規分布の確率分布関数をΦ（ｘｉ）とする。なお、離散データであるｘ＊ｉ０あるいはｐ＊ｉ０の確率密度関数の替わりに「頻度分布関数」や「ヒスとグラム」を用いてもよいが、以下には「確率密度関数」を用いて説明する。
１）正規分布に従う平均値に変換する方法
状態ｋでのｘ＊ｉｋとｐ＊ｉｋを次式でそれぞれ正規分布に従う平均値μｘｉｋｏとμｐｉｋ０に変換される。なお、状態ｋは任意の状態であり、基準状態も含む。
ここで、Φ−１はΦの逆関数で、σｘｉ０とσｐｉ０はそれぞれ正規分布に変換されたｘ＊ｉｋとｐ＊ｉｋとの標準偏差であり、次式で求められる。
平均値μｘｉｋ０とμｐｉｋ０を用いて状態判定や状態予測を行う。
２）直接変換法
状態ｋでのｘ＊ｉｋとｐ＊ｉｋを次式で正規分布の確率変数に変換する。
ここで、ＳｘｋとＳｐｋはそれぞれｘ＊ｉｋとｐ＊ｉｋの標準偏差で、μｘｋとμｐｋはそれぞれｘ＊ｉｋとｐ＊ｉｋの平均値である。ｘ’ｉｋ０あるいはｐ’ｉｋ０を用いて状態判定や状態予測を行う。
（２）基準状態を参照しない場合
まず、状態ｋでの波形データｘ＊ｉｋと特徴パラメータｐ＊ｉｋとの確率密度関数（あるいは、頻度分布）をそれぞれｆｘｋ（ｘ＊ｉｋ）とｆｐｋ（ｐ＊ｉｋ）とし、確率分布関数（あるいは、累積度数分布）をＦｘｋ（ｘ＊ｉｋ）とＦｐｋ（ｐ＊ｉｋ）とする。
１）直接変数変換法
波形データｘ＊ｉｋと特徴パラメータｐ＊ｉｋを次式で正規分布の確率変数に変換する。
ここで、ＳｘｋとＳｐｋはそれぞれｘ＊ｉｋとｐ＊ｉｋとの標準偏差で、μｘｋとμｐｋはｘ＊ｉｋとｐ＊ｉｋとの平均値である。ｘ’ｉｋあるいはｐ’ｉｋを用いて状態判定や状態予測を行う。
２）正規分布に従う平均値に変換する方法
正規分布に従う平均値は次式で求められる。
ここで、σｘｉｋとσｐｉｋは次式で求められる。
μｘｉｋとμｐｉｋは正規分布に従うので、μｘｉｋあるいはμｐｉｋを用いて状態判定や状態予測を行う。
３）間接変換法
Ｍ個の特徴パラメータｐ＊ｉｋの最小値と最大値をそれぞれ（ｐ＊ｉｋ）ｍｉｎと（ｐ＊ｉｋ）ｍａｘとする。（ｐ＊ｉｋ）ｍｉｎから（ｐ＊ｉｋ）ｍａｘまでＮ個の等間隔区間に分割する。各区間の代表値（例えば、中央値）をｐ＊ｉｋｊとする。ここで、ｊ＝１〜Ｎ。ｐ＊ｉｋｊをｐ＊ｉｋの替わりに式（３６）または式（３８）に代入すると、Ｎ個のｐ’’ｉｋまたはμ’ｐｉｋが得られる。ｐ’’ｉｋまたはμ’ｐｉｋを用いて状態判定や状態予測を行う。
同様に、Ｍ個の波形データｘ＊ｉｋの最小値と最大値をそれぞれ（ｘ＊ｉｋ）ｍｉｎと（ｘ＊ｉｋ）ｍａｘとする。（ｘ＊ｉｋ）ｍｉｎから（ｘ＊ｉｋ）ｍａｘまでＮ個の等間隔区間に分割する。各区間の代表値（例えば、中央値）をｘ＊ｉｋｊとする。ここで、ｊ＝１〜Ｋ。ｘ＊ｉｋｊをｘ＊ｉｋの替わりに式（３５）または式（３４）に代入すると、Ｋ個のｘ’’ｉｋまたはμ’ｘｉｋが得られる。ｘ’’ｉｋまたはμ’ｘｉｋを用いて状態判定や状態予測を行う。
４）特徴パラメータｐ＊ｉｋの区間平均値を求める方法
Ｎ個特徴パラメータｐ＊ｉｋを求めた後、Ｍ組に分割して第ｊ組の平均値は次のように求める。
ここで、Ｎｊは第ｊ組にあるｐ＊ｉｋの数である。μｋ（ｊ）は近似的に正規分布に従うので、μｋ（ｊ）を用いて状態判定や状態予測を行う。
波形データｘ＊ｉｋを式（２９）、式（３３）、式（３５）、式（３７）で正規確率分布の波形データに変換したμｘｉｋ０、ｘ’ｉｋ０、ｘ’ｉｋ、μｘｉｋ、ｘ’’ｉｋ、μ’ｘｉｋは「正規分布の波形データ」と呼ぶ。また、特徴パラメータｐ＊ｉを式（３０）、式（３４）、式（３６）、式（３８）で正規確率変数に変換したμｐｉｋ０、ｐ’ｉｋ０、ｐ’ｉｋ、μｐｉｋ、ｐ’’ｉｋ、μ’ｐｉｋ、μｋ（ｊ）は「正規分布の特徴パラメータ」と呼ぶ。
３．正規分布の特徴パラメータによる状態判別方法
ここで、正規分布の特徴パラメータを用いて、対象物の状態を判別する方法を述べる。
（１）統計理論による判別
１）正規分布の特徴パラメータの平均値の検定
状態ｋと状態ｙで求めた正規分布の特徴パラメータをそれぞれｐｉｋとｐｉｙとする。ここで、ｉ＝１〜Ｍ、Ｍは使用する正規分布の特徴パラメータの総数を表す。ｐｉｋとｐｉｙとの平均値をそれぞれμｉｋとμｉｙ、ｐｉｋとｐｉｙとの標準偏差をそれぞれＳｉｋとＳｉｙとする。一般にＪ個のｐｊの平均値μと標準偏差Ｓは次の式で計算する。
μｉｋとμｉｙが等しいか否かの検定は次のように行う（参考文献３参照）
（参考文献３）Ｋ．Ａ．Ｂｒｏｗｎｌｅｅ．Ｓｔａｔｉｓｔｉｃａｌ Ｔｈｅｏｒｙ ａｎｄ Ｍｅｔｈｏｄｏｌｏｇｙ ｉｎ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｅｃｏｎｄ Ｅｄｉｔｉｏｎ，Ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ Ｃｈｉｃａｇｏ，１９６５
が成立すれば、有意水準αで「μｉｋとμｉｙとが等しくない」と判定する。ここで、ｔα／２（Ｊ−１）は自由度Ｊ−１のｔ分布の確率密度関数が下側確率α／２に対するパーセント点である。
２）正規分布の特徴パラメータの分散の検定
ＳｉｋとＳｉｙが同じか否かの検定は次のように行う（参考文献４参照）。
（参考文献４）Ｋ．Ａ．Ｂｒｏｗｎｌｅｅ．Ｓｔａｔｉｓｔｉｃａｌ Ｔｈｅｏｒｙ ａｎｄ Ｍｅｔｈｏｄｏｌｏｇｙ ｉｎ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｅｃｏｎｄ Ｅｄｉｔｉｏｎ，Ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ Ｃｈｉｃａｇｏ，１９６５
が成立すれば、有意水準αで「ＳｉｋとＳｉｙとが等しくない」と判定する。ここで、Ｆα／２（Ｊ−１，Ｊ−１）は自由度Ｊ−１のＦ分布の確率密度関数が下側確率α／２に対するパーセント点である。
有意水準αを変えたとき、式（４４）または式（４５）を満足するか否かを確認することにより、状態ｙが状態ｋに対する状態変化の程度を決める。状態変化の程度を有意水準αにより決定する例は表１に示す。なお、設備診断の場合、状態ｋを正常状態とすれば、状態ｙは正常状態か、注意状態か、危険状態かの判別については、表１のように「正常」（α１）、「注意」（α２）、及び「危険」（α３）のように設定し検定することができる。つまり、式（４４）または式（４５）はα１の時に成立しなければ「正常」と判定する。また、式（４４）または式（４５）はα２の時に成立すれば「注意」と、α３の時に成立すれば「危険」と判定する。なお、表１中のαの数値範囲は例であり、設備の重要度などによって決められる。
数個の特徴パラメータを用いて状態判定を行う場合、判定の結果は最も状態変化が大きいと示した特徴パラメータの判定結果に従われる。例えば、３つの特徴パラメータｐ１、ｐ２、ｐ３を用いて判定するとき、ｐ１の判定結果が「注意」、ｐ２の判定結果が「正常」、ｐ３の判定結果が「危険」である場合、最終の判定結果は「危険」とする。
３）信頼区間による判定
基準時点で測定した波形データから求めた正規分布の特徴パラメータの平均値をμｉ０とし、他の時点で測定した波形データから求めた正規分布の特徴パラメータの平均値をμｉｋとすると、μｉ０の信頼区間は次式で与えられる。
ここで、ｔα／２（Ｊ−１）は自由度Ｊ−１のｔ分布の確率密度関数が下側確率α／２に対するパーセント点である。Ｓｉ０は波形データから求めたｐｉ０の標準偏差である。μｉｋは式（４６）に示す区間内にあれば、１−αの確率でμｉ０との差がないという。μｉ０の９９％信頼区間は、Ｊ＞１０のとき、近似的に次のようになる。
従って、μｉｋは式（４７）の範囲を超えれば、９９％の確率でμｉ０と違うと判定する。更に、測定した波形データで求めたμｉｋの信頼区間は、次式で求められる。
ここで、Ｓｉｋは波形データから求めたｐｉｋの標準偏差である。
表１のαを式（４６）に代入すれば、次のような信頼区間が得られる。
μｉｋがこれらの区間内にあるか否かによって状態の判定を行う。
数個の特徴パラメータを用いて状態判定を行う場合、最終の判定結果は式（４９），（５０），（５１）に記述される内容と同じである。
（２）可能性理論による判別
１）可能性分布関数の作成
状態ｋにおける波形データで正規分布の特徴パラメータｐｉの値を算出した後、ｐｉの確率密度関数ｆｋ（ｐｉ）から可能性分布関数Ｐｋ（ｐｉ）を式（５２）で求める。可能性理論によれば、ｐｉがどのような確率分布に従っても、その可能性分布関数が求められる。ｐｉが正規分布に従う場合、Ｎ段の可能性分布関数ｐｋ（ｐｉ）は次のように求める（参考文献５参照）。
ここで、
但し、上式中、ｐｉｘ＝ｍｉｎ｛ｐｉ｝＋ｘ×（ｍａｘ｛ｐｉ｝−ｍｉｎ｛ｐｉ｝）／Ｎ、ｘ＝１〜Ｎ，Ｓｉはｐｉの標準偏差，μｉはのｐｉ平均値である。
（参考文献５）Ｌ．Ｄａｖｉｓ：ＨＡＮＤＢＯＯＫ ＯＦ ＧＥＮＥＴＩＣ ＡＬＧＯＲＩＴＨＭＳ，Ｖａｎ Ｎｏｓｔｒａｎｄ Ｒｅｉｎｈｏｌｄ，Ａ Ｄｉｖｉｓｉｏｎ ｏｆ Ｗａｄｓｗｏｒｔｈ，Ｉｎｃ（１９９０）
２）可能性の求め方
図２のように、状態ｋと状態ｙで求めた正規分布の特徴パラメータｐｉの可能性分布関数をＰｋ（ｐｉ）とＰｙ（ｐｉ）とし、状態ｙで求めた正規分布の特徴パラメータの値をｐ’ｉとすると、「状態ｙが状態ｋと同じである」という可能性ｗは次のように求められる。
ａ）ｐ’ｉの平均値ｐ’ｉｍｅａｎとＰｋ（ｐｉ）とのマッチングによるｗの決定、
ｂ）Ｐｙ（ｐｉ）とＰｋ（ｐｉ）とのマッチングによるｗの決定。
なお、Ｐｙ（ｐｉ）とＰｋ（ｐｉ）とのマッチングによるｗを求める式は次に示す。
３）状態変化の判別
状態ｋで求めた正規分布の特徴パラメータｐｉの可能性分布関数ｐｋ（ｐｉ）が得られた後、左右両側の「状態変化が小さい」の可能性分布関数（ｐｃ１（ｐｉ）とｐｃ２（ｐｉ））、および「状態変化が大きい」の可能性分布関数（ｐｄ１（ｐｉ）とｐｄ２（ｐｉ））は図２に示すように決定する。境界値の
におけるｉ，ｊはユーザ入力により決定するが、標準値としてｉ＝３、ｊ＝６とする。
設備診断の場合、正常状態の可能性分布関数をｐｋ（ｐｉ）、注意状態の可能性分布関数をｐｃ１（ｐｉ）とｐｃ２（ｐｉ）、危険状態の可能性分布関数をｐｄ１（ｐｉ）とｐｄ２（ｐｉ）とする。実際の識別時に得られた「正常」と「注意」と「危険」との可能性は図２のように表示する。また、「危険」と判定した場合、警報を出すことも可能である。
（３）情報理論による状態判別法
対象物の基準状態における正規分布の特徴パラメータの確率密度関数をｆｐ０（ｐｉ）し、基準状態以外における正規分布の特徴パラメータの確率密度関数をｆｐｋ（ｐｉ）とする。なお、基準状態以外の時点での状態を「テスト状態」と呼ぶ。テスト状態が基準状態と同じ状態であるか否かは次の「Ｋｕｌｌｂａｃｋ−Ｌｅｉｂｌｅｒ Ｉｎｆｏｒｍａｔｉｏｎ（カルバック情報量、ＫＩ）」と「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅ（情報ダイバージェンス、ＩＤ）」で判定できる。
ＫＩｐとＩＤｐによる状態判別法は［１０］に詳細に記述しているので、ここで省略する。（［１０］劉 信芳、豊田利夫、陳 鵬、馮 芳、二保知也：「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅによる回転機械の異常診断」、精密工学会誌、Ｖｏｌ．６６、Ｎｏ．１、２０００年、ｐ．１５７−１６２．）
（４）複数の特徴パラメータを統合することによる状態判別法
数個の特徴パラメータを統合して、状態変化の有無や状態予測を行うこともできる。特徴パラメータの統合法は、主成分分析法やＫＬ展開法などがあるが、特徴パラメータの統合法により求めた新たな特徴パラメータは「統合特徴パラメータ」とよぶ。ここで、主成分分析法の例を示す。（［１１］大津，栗田，関田著：パターン認識，朝倉書店、１９９６．）
例えば、設備診断の場合、設備の正常状態下で求めた無次元特徴パラメータ（ｐ１，ｐ２，・・・，ｐｍ）に対して、ｍ個の主成分は次のように表されます。
各主成分ｚ１〜ｚｍは「統合特徴パラメータ」とも呼ぶ。
相関行列は次のように求められます。
ここで、正常状態で収録したｎ組のデータを｛ｐ１ｋ，ｐ２ｋ，・・・，ｐｍｋ｝、ｋ＝１，２，・・・，ｎとすると、
相関行列Ｒの固有ベクトル、λ＝（λ１，λ２，・・・，λｍ）、λ１，≧λ２≧・・・≧λｍとすると、
固有値λｉに対応する式（５９）の係数が求まり、第ｉ番目の主成分は次のように求められます。
主成分を用いて次式のように状態判定を行う。
ここで、Ｋは使用する主成分の数、αは有意水準、χ２（Ｋ，α）は自由度Ｋのカイ２乗分布の確率密度関数が上側確率αに対するパーセント点である。αの決定法は表１同じである。なお、Ｋ＝３に対してα＝０．０５のとき、χ２（３，０．０５）＝７．８１５である。
なお、式（６４）に示す判定方法以外に、特徴パラメータの統合法により求めた統合特徴パラメータ、例えば、式（６３）に示す主成分、を正規分布の確率変数に変換した後、統計検定や可能性理論により状態判定を行うこともできる。
ここで、実例を示す。図３はある回転機械の正常状態時（図３（ａ））と回転軸ミスアライメント状態時（図３（ｂ））に測定した振動の加速度波形データである。図４（ａ）及び図５（ａ）は式（２）〜式（１３）及び式（１８）〜式（２１）に示す１４個の無次元特徴パラメータを６０回求め、正規分布変換を行う前に、Ｋ＝３のときに求めた式（６４）右側の値（６０個）を示す。図４（ｂ）及び図５（ｂ）は式（３１）に基づいて正規分布変換を行った後、Ｋ＝３のときに求めた式（６４）右側の値（９個）を示す。これらの図により、正規分布変換後の結果は正規分布変換前の結果より良好であることがわかる。なお、「良好である」とは、式（６４）左側の値は正常状態の時に小さく、異常状態の時に大きいということである。
４．正規分布の波形データによる状態判別方法
（１）特徴パラメータによる状態判別法
対象物から計測された波形データｘ＊ｉ（例えば、図６（ａ）、（ｃ））を式（３３）により正規分布の波形データｘ’ｉｋ０（例えば、図６（ｂ）、（ｄ））に変換した後、正規分布の波形データｘ’ｉｋ０を用いて前述の正規分布の特徴パラメータによる状態判別方法で状態判定を行うことができる。なお、この方法は正規分布の波形データμｘｉｋ０とｘ’ｉｋ０に適用できる。
（２）情報理論による状態判別法
基準状態の正規分布の波形データの確率密度関数をｆｘ０（ｘｉ）とし、基準状態以外の時点での正規分布の波形データの確率密度関数をｆｘｋ（ｘｉ）とする。なお、基準状態以外の時点での状態を「テスト状態」と呼ぶ。テスト状態が基準状態と同じ状態であるか否かは次の「Ｋｕｌｌｂａｃｋ−Ｌｅｉｂｌｅｒ Ｉｎｆｏｒｍａｔｉｏｎ（カルバック情報量、ＫＩ）」と「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅ（情報ダイバージェンス、ＩＤ）」で判定できる。
ＫＩｘとＩＤｘによる状態判別法はに詳細に記述しているので、ここで省略する。（［１２］劉 信芳、豊田利夫、陳 鵬、馮 芳、二保知也：「Ｉｎｆｏｒｍａｔｉｏｎ Ｄｉｖｅｒｇｅｎｃｅによる回転機械の異常診断」、精密工学会誌、Ｖｏｌ．６６、Ｎｏ．１、２０００年、ｐ．１５７−１６２．）
ここで、実例を示す。図６（ａ）はある回転機械の正常状態時に測定した振動加速度の波形データである。図６（ｃ）は同じ回転機械のアンバランス時に測定した波形データである。
平均値以上の波形データｘｉ＋、平均値以下の波形データｘｉ−、及び絶対値データ｜ｘｉ｜を正規分布の波形データに変換した後、それぞれｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’で表す。ｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’の平均値と標準偏差は次の通りである。
ｘ’ｉ＋：正常時の平均値＝０．６９、アンバランス時の平均値＝１．９６
正常時の標準偏差＝１．０３、アンバランス時の標準偏差＝２．３７
ｘ’ｉ−：正常時の平均値＝−０．３７、アンバランス時の平均値＝−０．６６
正常時の標準偏差＝０．４７、アンバランス時の標準偏差＝０．６９
｜ｘｉ｜’：正常時の平均値＝１．４５、アンバランス時の平均値＝３．２９
正常時の標準偏差＝１．５２、アンバランス時の標準偏差＝３．５３
（参考文献８）によると、過検出率α＝０．１５、見逃率β＝０．１５とすると、ＫＩ＞１．２１或いはＩＤ＞２．４３ならば、８５％ｎの確率で「テスト状態は基準状態と違う」と判定できる。上記のｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’のＫＩとＩＤは次の通りである。
ｘ’ｉ＋：ＫＩ＝０．５１９， ＩＤ＝２．２４
ｘ’ｉ−：ＫＩ＝０．５３， ＩＤ＝０．４０
｜ｘｉ｜’：ＫＩ＝０．５７， ＩＤ＝２．６７
｜ｘｉ｜’のＩＤ＞２．４３から、８５％ｎの確率で「図６（ａ）の状態は図６（ｃ）の状態と違う」と判定できる。
また、式（４４）と式（４５）に示す平均値と分散の検定により状態判定を行ってもよい。
更に、上記のｘ’ｉ＋、ｘ’ｉ−、及び｜ｘｉ｜’を用いて、式（２）〜式（２５）で特徴パラメータを求めて、統計検定や可能性理論及び特徴パラメータの統合法により状態判定を行うこともできる。
５．状態予測
正規分布の特徴パラメータｐｉが求められた後、従来の状態予測手法を用いて、測定対象の状態を予測することが行える（［１３］石川真澄、武藤博道：予測手法、計測と制御、１９８２．３．［１４］Ｏｇａｗａ，Ｍ．：Ｔｉｍｅ ｓｅｒｉｅｓ ａｎａｌｙｓｉｓ ａｎｄ ｓｔｏｃｈａｓｔｉｃ ｐｒｅｄｉｃｔｉｏｎ，Ｂｕｌｌ．Ｍａｔｈ．Ｓｔａｔ．，８，８−７２，１９５８．［１５］Ｂ．ピカニオル著、小村賢二、柴山宮恵子訳：予測のための統計学、晃洋書房、１９８７．）
図７は状態予測の方法を示す。各測定時点（ｘ１〜ｘ７）で測定した特徴パラメータ値、或いは、主成分値を正規確率分布に従うように変換した後、回帰解析により予測曲線とその信頼区間を求め、「寿命限界」との交点で「最短寿命」、「平均寿命」及び「最長寿命」を求める。
図８はある回転機械が正常状態からアンバランス状態へ変化していく間に８回測定した波形例である。測定１は正常状態であり、測定８は最も程度の重いアンバランス状態である。図９は式（２）〜式（１３）で求めた、各測定における特徴パラメータの平均値を示す。この例でわかるように、特徴パラメータは異常状態の程度が重くなっていくにつれて値が単調増加の特徴パラメータ（例えば、ｐ４，ｐ５）があれば、値が単調減少の特徴パラメータ（例えば、ｐ２，ｐ８，ｐ９，ｐ１０）もある。また、異常状態の程度が重くなっても値がほぼ不変な特徴パラメータ（例えば、ｐ１，ｐ７）もある。従って、値が単調増加（或いは、単調減少）の特徴パラメータを選んで、寿命予測を行うべきである。また、主成分により寿命予測を行う時、値が単調増加の特徴パラメータのみ（或いは、単調減少の特徴パラメータのみ）を選んで、主成分を求めて寿命予測を行うべきである。
上記に述べた波形データと特徴パラメータを正規分布の確率変数に変換する手法、状態判定手法及び状態予測手法を実現するための計測と処理の流れは図１６（ａ）に示す。また、図１６（ａ）を実現するための波形データ計測及び状態判定装置の回路は図１７に示す。
６．脈動信号から特異成分の検出法及び状態判定法
ここで、ガスエンジン失火診断の例を用いて説明する。
図１０（ａ）、図１１（ａ）、図１２（ａ）、図１３（ａ）はガスエンジン失火（図中のピーク値の異常個所）が発生したときのシリンダ圧力波形データである。なお、図１２（ａ）、図１３（ａ）は多数の異常ピークが連続して発生する例を示す。このような脈動信号に発生する局所異常（特異成分）に対して、オンライン特異点の検出及び状態判定の流れは図１５に示す。以下、その手順を説明する。
１）．診断ための準備（図１５（ａ））
（１）診断対象の脈動信号と回転パルス信号を同時に計測する。なお、回転パルス信号は周期パルス信号とも言い、脈動信号の各ピーク値のタイミングを決めるときに用いられる。
（２）雑音を除去するために，ローパスフィルタリングを行う。ローパスフィルタのカットオフ周波数ｆＬは次のように決定する。
ここで、ｎは軸の回転数（ｒｐｍ）、ｚは１回転毎のピーク数（ピーク／１回転）、ｆ０は余裕周波数（＞ｎ／６０、フィルタリング後のノイズ除去効果を観察することにより決定する）である。
（３）ローパスフィルタリング後の包絡線波形データ、または、ピーク波形データ、または、移動平均波形データを求める。
（４）上記の包絡線波形データ、または、ピーク波形データ、または、移動平均波形データの正規化を次式のように行う。
ここで、ｘ（ｔ）は正規化後の波形データ、ｘ’（ｔ）は元の波形データ、μ’（ｔ）はｘ’（ｔ）の平均値、ｓはｘ’（ｔ）の標準偏差である。なお、μ’（ｔ）とｓを求めるとき、正常状態の波形データが望ましいが、正常状態の波形データでなくてもよい。また、μ’（ｔ）とｓを求めるとき、運転条件（回転数や負荷）をできるだけ判定のときと同じように設定する。但し、運転条件が多少変化しても判定結果に大きく影響しない。
（５）回転パルス信号のピークと脈動信号のピークとの相対位置関係を調べて置き、異常を識別するための閾値を決定する。閾値の決定は次のように考える。
のとき、異常が発生したと判定する。ここで、ｋ（＝２〜４）は診断対象物によって決定する。
｜ｘ（ｔ）｜がｋσより大きい間、回転パルス信号に対応する脈動信号のピークの所に異常があると判定する。
２）．オンライン診断（図１５（ｂ））
（１）判定対象物の脈動信号と回転パルス信号を同時に計測する。
（２）ローパスフィルタリングを行う。ローパスフィルタのカットオフ周波数ｆＬは式（６６）のように決定する。
（３）ローパスフィルタリング後の包絡線波形データ、または、ピーク波形データを求める。
（４）上記の包絡線波形データ、または、ピーク波形データの正規化は式（６８）により求める。
（５）正規化された包絡線波形データ、または、ピーク波形データの絶対値は閾値（ｋσ）より大きいか否かを監視する。大きくなければ、異常がないと判定する。大きければ、異常が発生したと判定し、｜ｘ（ｔ）｜がｋσより大きい間、前もって測って置いた回転パルス信号と脈動信号のピークとの関係により、異常の箇所を判定する。
上記の手順に基づく信号処理の例を図１０、図１１、図１２、図１３に示す。
図１０（ａ）（ｂ）、図１１（ａ）（ｂ）、図１２（ａ）（ｂ）、図１３（ａ）（ｂ）ガスエンジン失火（図中の異常箇所）が発生したときのシリンダ圧力の時系列波形データとスペクトルである。
図１０（ｃ）（ｄ）、図１１（ｃ）（ｄ）、図１２（ｃ）（ｄ）、図１３（ｃ）（ｄ）はローパスフィルタにより雑音を除去した時系列波形データとスペクトルである。
図１０（ｅ）、図１１（ｅ）、図１２（ｅ）、図１３（ｅ）は回転パルス波形データである。
図１０（ｆ）、図１２（ｆ）は包絡線波形データである。
図１１（ｆ）、図１３（ｆ）はピーク波形データである。
図１０，１１，１２，１３の（ｆ）で分かるように、異常（ピークの欠落、または、小さくなる）が発生した場合、その絶対値｜ｘ（ｔ）｜が２σより大きくなることが判明できる。
また、絶対値｜ｘ（ｔ）｜が２σより大きくなる箇所と回転パルス波形データとの対応関係により異常箇所を判明できる。
上記に包絡線波形データとピーク波形データの例を示したが、図１４に示すように、移動平均波形データ（ｕｉ）を求めたら、ｘ（ｔ）の替わりにｕ（ｔ）を用いて上記の手順と同様に脈動信号の特異点検出ができる。なお、カットオフ周波数（ｆｃ）は式（６７）により求めたら、移動平均の点数（Ｍ）が求められる。
上記に述べた脈動信号の特異成分の検出と判定方法を実現するための計測と処理の流れは図１６（ｂ）に示す。また図１６（ｂ）を実現するための脈動信号計測及び状態判定装置の回路は図１８に示す。
７．状態判定装置、または、オンライン状態判定システムの処理流れ
状態判定装置装置、または、オンライン状態判定システムの処理流れは図１６に示す。図１６（ａ）に示すように、測定波形データからノイズを除去し、特徴パラメータを求め、その特徴パラメータを正規分布の確率変数に変換した後、統計検定や可能性理論や特徴パラメータの統合により状態変化の有無判定を行う。これらの処理は計算機または専用装置で実現できる。また、図１６（ｂ）に示すように、脈動信号からの特異成分検出及び状態判定は、ローパスフィルタ、包絡線（または、ピーク値）はハードウェアで実現できる。数値演算装置（または、計算機）は正規化したｘ（ｔ）、｜ｘ（ｔ）＞ｋσ｜の判定、特異箇所の検出と判定および判定結果の表示を行えばよいから、リアルタイムに行える。
【図面の簡単な説明】
図１は、測定した波形データを４種類、つまり、平均値より大きいデータｘｉ＋、平均値より小さいデータｘｉ−、式（１）による正規化後の絶対値データ｜ｘｉ｜及び全体波形データｘ＊ｉＡに分ける例を示すグラフである。
図２は可能性分布関数の例を示すグラフである。
図３は正常状態と回転軸ミスアライメント状態の振動波形データの例を示すグラフである。
図４は正常状態の振動波形データに対する主成分分析結果を示すグラフである。
図５は回転軸ミスアライメント状態の振動波形データに対する主成分分析結果を示すグラフである。
図６は正常状態とアンバランス状態の振動波形データの例を示すグラフである。
図７は状態予測方法を示すグラフである。
図８は正常状態からアンバランス状態へ変化していく時の振動波形データの例を示すグラフである。
図９は各測定の波形の特徴パラメータ値を示すグラフである。
図１０は包絡線波形データによる１個のピークが欠落した場合の信号処理例を示すグラフである。
図１１はピーク波形データによる１個のピークが欠落した場合の信号処理例を示すグラフである。
図１２は包絡線波形データによる多数のピーク異常の場合の信号処理例を示すグラフである。
図１３はピーク波形データによる多数のピーク異常の場合の信号処理例を示すグラフである。
図１４は移動平均波形データの求め方を示すグラフである。
図１５は脈動信号の特異成分検出及び異常診断処理の流れを示すフローチャートである。
図１６は状態判定装置、または、オンライン状態判定システムの処理流れを示すフローチャートである。
図１７は信号計測及び状態判定装置の回路の一例を示す回路図であり、図中の符号は次の通りである。
１ センサ、２ アンプ、３ フィルタ、４ 信号処理・計算装置、５ 表示出力装置、６ データ用ＲＡＭ、７ ＡＤ変換器、８ ＤＣポート、９ ＳＣＩ（シリアル・コミュニケーション・インタフェース）、１０ ＣＰＵ、１１ フラッシュＲＯＭ、１２ 外部コンピュータ．
図１８は脈動信号計測及び状態判定装置の回路の一例を示す回路図であり、図中の符号は次の通りである。
１ センサ、２ アンプ、３ フィルタ、４ 包絡線処理部、５ 周期パルスセンサ、６ 表示出力装置、７ ＡＤ変換器、８ ＤＣポート、９ ＳＣＩ（シリアル・コミュニケーション・インタフェース）、１０ ＣＰＵ、１１ フラッシュＲＯＭ、１２ 外部コンピュータ、１３ データ用ＲＡＭ、１４ 信号処理・計算装置．
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
【特許請求の範囲】
【請求項１】
状態監視の対象物の波形データを複数の周波数帯域において測定するステップと、前記測定された波形データから雑音を除去するステップと、前記雑音を除去された波形データを用いて時間領域及び周波数領域の特徴パラメータを原始特徴パラメータとして定義し計算するステップと、前記計算された原始特徴パラメータを、予め指定した確率分布に従う、既知分布の特徴パラメータに変換するステップと、前記既知分布の特徴パラメータを用いて、対象物の状態判定及び状態予測を行うステップと、を有し、対象物の状態判定と状態予測を行う方法であって、前記時間領域及び周波数領域の特徴パラメータを原始特徴パラメータとして定義し計算するステップにおいて、前記原始特徴パラメータは無次元特徴パラメータであれば、前記雑音を除去された波形データを正規化した後、算出されることと、前記原始特徴パラメータを既知分布の特徴パラメータに変換するステップにおいて、対象物の基準状態を参照する場合と、対象物の基準状態を参照しない場合と、に分けて変換することと、を特徴とする方法であり、この方法を実行するために、対象物の波形データを取得するためのセンサ、アンプ、フィルタ、信号処理・計算装置、及び表示出力装置を備えた、対象物の状態判定と状態予測を行う装置。
【請求項２】
前記原始特徴パラメータは、前記正規化された波形データを平均値より大きい波形データ、平均値より小さい波形データ、絶対値波形データ及び全体波形データに分けて正規分布に変換した後、これらの波形データにより算出される方法。
【請求項３】
前記原始特徴パラメータのうち、区間特徴パラメータの定義と計算を行う方法。
【請求項４】
前記原始特徴パラメータを既知分布の特徴パラメータに変換するステップにおいて、前記原始特徴パラメータの確率分布と平均値と分散を求めるステップと、前記原始特徴パラメータの確率分布と平均値と分散を用いて既知分布の特徴パラメータの値を求めるステップと、を有する。
【請求項５】
前記既知分布の特徴パラメータを用いて、統計理論と可能性理論と特徴パラメータの統合法と情報理論とにより対象物の状態判定及び状態予測を行う方法。
【請求項６】
前記既知分布の特徴パラメータを用いて特徴パラメータの統合法により求めた統合特徴パラメータを、予め指定した確率分布に従う、既知分布の統合特徴パラメータに変換し、既知分布の統合特徴パラメータを用いて統計理論と可能性理論により対象物の状態判定及び状態予測を行う方法。
【請求項７】
前記の諸方法において、状態変化につれて値が単調増加する特徴パラメータあるいは統合特徴パラメータのみ、あるいは、値が単調減少する特徴パラメータあるいは統合特徴パラメータのみ、を選出して、状態判定と状態予測を行う方法。
【請求項８】
前記予め指定した確率分布が正規分布である場合の状態判定と状態予測を行う方法。
【請求項９】
対象物の状態監視のために、波形データを複数の周波数帯域において測定するステップと、前記測定された波形データから雑音を除去するステップと、前記雑音を除去された波形データを、予め指定した確率分布に従う、既知分布の波形データに変換するステップと、前記既知分布の波形データを用いて対象物の状態判定及び状態予測を行うステップと、を有し、対象物の状態判定と状態予測を行う方法であって、前記雑音を除去された波形データを既知分布の波形データに変換するステップにおいて、波形データを正規化し、対象物の基準状態を参照する場合と、対象物の基準状態を参照しない場合と、に分けて変換することを特徴とする方法であり、この方法を実行するために、対象物の波形データを取得するためのセンサ、アンプ、フィルタ、信号処理・計算装置、及び表示出力装置を備えた、対象物の状態判定と状態予測を行うコンピュータ装置。
【請求項１０】
前記雑音を除去された波形データを平均値より大きい波形データ、平均値より小さい波形データ、絶対値波形データ及び全体波形データに分けて、予め指定した確率分布に従う、既知分布の波形データ変換する方法。
【請求項１１】
前記既知分布の波形データを用いて特徴パラメータと情報理論により対象物の状態判定及び状態予測を行う方法。
【請求項１２】
前記予め指定した確率分布が正規分布である場合、波形データを正規分布の波形データに変換する方法。
【請求項１３】
センサーで測定した対象物の信号が脈動信号である場合は、前記測定された脈動信号から雑音を除去するステップと、前記雑音を除去された脈動信号の特徴波形データを求めるステップと、前記特徴波形データの正規化を行うステップと、前記正規化された特徴波形データにより特異成分を検出し、特異成分の発生位置を特定するステップと、を有し、対象物の状態判定を行う方法であって、前記特徴波形データは、前記雑音を除去された脈動信号の包絡線波形データである場合と、前記雑音を除去された脈動信号のピーク波形データである場合と、前記雑音を除去された脈動信号の移動平均波形データと、に分けて求められることを特徴とする方法であり、この方法を実行するために、対象物の脈動信号を取得するためのセンサ、周期パルス信号を取得するための周期パルスセンサ、アンプ、フィルタ、包絡線処理部、信号処理・計算装置、及び表示出力装置を備えた、対象物の脈動信号の特異成分の検出と状態判定を行う装置。
【請求項１４】
前記測定された脈動信号からローパスフィルタを用いて雑音を除去する方法。
【請求項１５】
前記正規化された特徴波形データと、同時に測定した周期パルス波形データにより特異成分を検出し、特異成分の発生位置を特定する方法。
【請求項１】
状態監視の対象物の波形データを複数の周波数帯域において測定するステップと、前記測定された波形データから雑音を除去するステップと、前記雑音を除去された波形データを用いて時間領域及び周波数領域の特徴パラメータを原始特徴パラメータとして定義し計算するステップと、前記計算された原始特徴パラメータを、予め指定した確率分布に従う、既知分布の特徴パラメータに変換するステップと、前記既知分布の特徴パラメータを用いて、対象物の状態判定及び状態予測を行うステップと、を有し、対象物の状態判定と状態予測を行う方法であって、前記時間領域及び周波数領域の特徴パラメータを原始特徴パラメータとして定義し計算するステップにおいて、前記原始特徴パラメータは無次元特徴パラメータであれば、前記雑音を除去された波形データを正規化した後、算出されることと、前記原始特徴パラメータを既知分布の特徴パラメータに変換するステップにおいて、対象物の基準状態を参照する場合と、対象物の基準状態を参照しない場合と、に分けて変換することと、を特徴とする方法であり、この方法を実行するために、対象物の波形データを取得するためのセンサ、アンプ、フィルタ、信号処理・計算装置、及び表示出力装置を備えた、対象物の状態判定と状態予測を行う装置。
【請求項２】
前記原始特徴パラメータは、前記正規化された波形データを平均値より大きい波形データ、平均値より小さい波形データ、絶対値波形データ及び全体波形データに分けて正規分布に変換した後、これらの波形データにより算出される方法。
【請求項３】
前記原始特徴パラメータのうち、区間特徴パラメータの定義と計算を行う方法。
【請求項４】
前記原始特徴パラメータを既知分布の特徴パラメータに変換するステップにおいて、前記原始特徴パラメータの確率分布と平均値と分散を求めるステップと、前記原始特徴パラメータの確率分布と平均値と分散を用いて既知分布の特徴パラメータの値を求めるステップと、を有する。
【請求項５】
前記既知分布の特徴パラメータを用いて、統計理論と可能性理論と特徴パラメータの統合法と情報理論とにより対象物の状態判定及び状態予測を行う方法。
【請求項６】
前記既知分布の特徴パラメータを用いて特徴パラメータの統合法により求めた統合特徴パラメータを、予め指定した確率分布に従う、既知分布の統合特徴パラメータに変換し、既知分布の統合特徴パラメータを用いて統計理論と可能性理論により対象物の状態判定及び状態予測を行う方法。
【請求項７】
前記の諸方法において、状態変化につれて値が単調増加する特徴パラメータあるいは統合特徴パラメータのみ、あるいは、値が単調減少する特徴パラメータあるいは統合特徴パラメータのみ、を選出して、状態判定と状態予測を行う方法。
【請求項８】
前記予め指定した確率分布が正規分布である場合の状態判定と状態予測を行う方法。
【請求項９】
対象物の状態監視のために、波形データを複数の周波数帯域において測定するステップと、前記測定された波形データから雑音を除去するステップと、前記雑音を除去された波形データを、予め指定した確率分布に従う、既知分布の波形データに変換するステップと、前記既知分布の波形データを用いて対象物の状態判定及び状態予測を行うステップと、を有し、対象物の状態判定と状態予測を行う方法であって、前記雑音を除去された波形データを既知分布の波形データに変換するステップにおいて、波形データを正規化し、対象物の基準状態を参照する場合と、対象物の基準状態を参照しない場合と、に分けて変換することを特徴とする方法であり、この方法を実行するために、対象物の波形データを取得するためのセンサ、アンプ、フィルタ、信号処理・計算装置、及び表示出力装置を備えた、対象物の状態判定と状態予測を行うコンピュータ装置。
【請求項１０】
前記雑音を除去された波形データを平均値より大きい波形データ、平均値より小さい波形データ、絶対値波形データ及び全体波形データに分けて、予め指定した確率分布に従う、既知分布の波形データ変換する方法。
【請求項１１】
前記既知分布の波形データを用いて特徴パラメータと情報理論により対象物の状態判定及び状態予測を行う方法。
【請求項１２】
前記予め指定した確率分布が正規分布である場合、波形データを正規分布の波形データに変換する方法。
【請求項１３】
センサーで測定した対象物の信号が脈動信号である場合は、前記測定された脈動信号から雑音を除去するステップと、前記雑音を除去された脈動信号の特徴波形データを求めるステップと、前記特徴波形データの正規化を行うステップと、前記正規化された特徴波形データにより特異成分を検出し、特異成分の発生位置を特定するステップと、を有し、対象物の状態判定を行う方法であって、前記特徴波形データは、前記雑音を除去された脈動信号の包絡線波形データである場合と、前記雑音を除去された脈動信号のピーク波形データである場合と、前記雑音を除去された脈動信号の移動平均波形データと、に分けて求められることを特徴とする方法であり、この方法を実行するために、対象物の脈動信号を取得するためのセンサ、周期パルス信号を取得するための周期パルスセンサ、アンプ、フィルタ、包絡線処理部、信号処理・計算装置、及び表示出力装置を備えた、対象物の脈動信号の特異成分の検出と状態判定を行う装置。
【請求項１４】
前記測定された脈動信号からローパスフィルタを用いて雑音を除去する方法。
【請求項１５】
前記正規化された特徴波形データと、同時に測定した周期パルス波形データにより特異成分を検出し、特異成分の発生位置を特定する方法。
【国際公開番号】ＷＯ２００４／０６８０７８
【国際公開日】平成１６年８月１２日（２００４．８．１２）
【発行日】平成１８年５月２５日（２００６．５．２５）
【国際特許分類】
【出願番号】特願２００５−５０４６６５（Ｐ２００５−５０４６６５）
【国際出願番号】ＰＣＴ／ＪＰ２００４／０００１６３
【国際出願日】平成１６年１月１３日（２００４．１．１３）
【出願人】（５９９０６４８７８）
【Ｆターム（参考）】
【国際公開日】平成１６年８月１２日（２００４．８．１２）
【発行日】平成１８年５月２５日（２００６．５．２５）
【国際特許分類】
【国際出願番号】ＰＣＴ／ＪＰ２００４／０００１６３
【国際出願日】平成１６年１月１３日（２００４．１．１３）
【出願人】（５９９０６４８７８）
【Ｆターム（参考）】
[ Back to top ]