(57)【特許請求の範囲】
【請求項１】
状態診断やパターン認識のために対象物から測定した信号の波形データについてオーバーフロー有無ならびに下記の「外れ値程度算出法」と下記の「波形非安定程度算出法」を用いて総合的に評価することにより前記波形データが状態診断やパターン認識に使用できるか否かを検査して決定することを特徴とする信号検査モジュール。
「外れ値程度算出法」：
状態診断やパターン認識のために得られた信号の波形データをｘ（ｉ）（ｉ＝１〜Ｎ）とし、ｘ（ｉ）の絶対値を｜ｘ（ｉ）｜とし、｜ｘ（ｉ）｜の平均値と標準偏差をそれぞれμとＳとし、係数ｋの範囲をｋｍｉｎ〜ｋｍａｘと設定した後、｜ｘ（ｉ）｜の中にある値｜ｘ（ｊ）｜（ｊ＝１〜Ｎ）＞μ＋ｋｍｉｎＳならば、ｘ（ｊ）を外れ値と判定し、
または、ｘ（ｉ）のＬ個の絶対最大値を｜ｘＭｊ｜（ｊ＝１〜Ｌ）とし、ｘＭｊ除いた後の平均値と標準偏差をそれぞれμａとＳａとすると、係数ｋａの範囲をｋａｍｉｎ〜ｋａｍａｘと設定した後、｜ｘＭｊ｜＞μａ＋ｋａｍｉｎＳａならば、ｘＭｊを外れ値と判定し、
外れ値の発生程度（以下「外れ程度」と呼ぶ）をパーセンテージで表す場合、係数ｋ（または、ｋａ）の上限値と下限値をそれぞれｋｍａｘとｋｍｉｎ（または、ｋａｍａｘとｋａｍｉｎ）とし、｜ｘ（ｊ）｜≧μ＋ｋｍａｘＳ（または、｜ｘＭｊ｜≧μａ＋ｋａｍａｘＳａ）のとき、外れ程度を１００％とし、｜ｘ（ｊ）｜≦μ＋ｋｍｉｎＳ（または、｜ｘＭｊ｜≦μａ＋ｋａｍｉｎＳａ）のとき、外れ程度を０％とし、｜ｘ（ｊ）｜＝μ＋ｋｍｉｎＳ〜μ＋ｋｍａｘＳ（または、｜ｘＭｊ｜＝μａ＋ｋａｍｉｎＳａ〜μａ＋ｋａｍａｘＳａ）のとき、外れ程度を０％〜１００％とすると、外れ程度が予め設定した閾値より大きい時に前記信号を不備な信号と判定する。
「波形非安定程度算出法」：
「波形非安定程度算出法」は下記の「（１）有次元特徴パラメータの安定性による波形非安定程度算出法」と下記の「（２）無次元特徴パラメータの安定性による波形非安定程度算出法」とから構成され、状態診断やパターン認識などのために得られた信号に対して、下記の「（１）有次元特徴パラメータの安定性による波形非安定程度算出法」を用いて算出した有次元特徴パラメータの非定常程度が予め設定した閾値より大きい場合、または、下記の「（２）無次元特徴パラメータの安定性による波形非安定程度算出法」を用いて算出した無次元特徴パラメータの非定常程度が予め設定した閾値より大きい場合、前記信号を非安定な信号と判定する。
（１）有次元特徴パラメータの安定性による波形非安定程度算出法：
状態診断やパターン認識などのために得られた信号の波形データｘ（ｉ）をＭ区間に分割し、第ｊ区間と第ｋ区間の波形データをそれぞれｘｊ（ｉ）とｘｋ（ｉ）で表し、第ｊ区間と第ｋ区間の波形データ数をそれぞれＮｊとＮｋとし、第ｊ区間と第ｋ区間における有次元特徴パラメータの平均値をそれぞれμｊとμｋとし、第ｊ区間と第ｋ区間における有次元特徴パラメータの標準偏差をそれぞれＳｊとＳｋとし、Ｍ個の区間における平均値μｉ（ｉ＝１〜Ｍ）が全て等しい（すなわちμ１＝μ２＝・・・μｊ＝μｋ＝・・・＝μＭ）またはＭ個の区間における標準偏差Ｓｉ（ｉ＝１〜Ｍ）が全て等しい（すなわちＳ１＝Ｓ２＝・・・Ｓｊ＝Ｓｋ＝・・・＝ＳＭ）という帰無仮説を統計理論により検定するために与えられた有意水準αの最大値と最小値をそれぞれαｍａｘとαｍｉｎとを設定した後、αを０から１まで変化させた時に前記の帰無仮説が丁度棄却された時点でのαをαｄとすると、αｄ＜αｍｉｎのときに有次元特徴パラメータが非定常と判定し、
有次元特徴パラメータの非安定程度（以下、「有次元非安定程度」とよぶ）をパーセンテージで表す場合、αｄ≦αｍｉｎのとき、有次元非安定程度を１００％とし、αｄ≧αｍａｘのとき、有次元非安定程度を０％とし、αｄ＝αｍｉｎ〜αｍａｘのとき、有次元非安定程度を１００％〜０％とすると、有次元非安定程度が予め設定した閾値より大きい時に前記信号を非安定な信号と判定する。
（２）無次元特徴パラメータの安定性による波形非安定程度算出法：
状態診断やパターン認識などのために得られた信号の波形データｘ（ｉ）から算出された無次元特徴パラメータをｐとし、ｘ（ｉ）における無次元特徴パラメータｐの平均値と標準偏差をそれぞれμｐとＳｐとし、ｘ（ｉ）をＭ区間に分割し、第ｊ区間の波形データ数をＮｊとし、第ｊ区間の波形データをｘｊ（ｉ）（ｊ＝１〜Ｍ）で表し、ｘｊ（ｉ）における無次元特徴パラメータをｐｊで表すと、
「無次元特徴パラメータの安定性による波形非安定程度算出法」は下記の「無次元非定常程度算出法１」と下記の「無次元非定常程度算出法２」と下記の「無次元非定常程度算出法３」と下記の「無次元非定常程度算出法４」から構成され、前記信号に対して、下記の「無次元非定常程度算出法１」と下記の「無次元非定常程度算出法２」と下記の「無次元非定常程度算出法３」と下記の「無次元非定常程度算出法４」とのうちの少なくとも一方を用いて算出した無次元特徴パラメータの非安定程度が予め設定した閾値より大きい場合、前記信号を非安定な信号と判定する。
無次元非定常程度算出法１：
予め係数ｋｐの範囲をｋｐｎｉｎ〜ｋｐｍａｘと設定した後、｜ｐｊ｜＞μｐ＋ｋｐｍｉｎＳｐならばｐｊを特異な無次元特徴パラメータと判定し、
無次元特徴パラメータの非安定程度（以下、「無次元非安定程度」とよぶ）をパーセンテージで表す場合、｜ｐｊ｜≧μｐ＋ｋｐｍａｘＳｐのとき、無次元非安定程度を１００％とし、｜ｐｊ｜≦μｐ＋ｋｐｍｉｎＳｐのとき、無次元非安定程度を０％とし、｜ｐｊ｜＝μｐ＋ｋｐｍｉｎＳｐ〜μｐ＋ｋｐｍａｘＳｐのとき、無次元非安定程度を０％〜１００％とすると、無次元非安定程度が予め設定した閾値より大きい時に前記信号を非安定な信号と判定する。
無次元非定常程度算出法２：
予め係数ｋＬの範囲をｋＬｎｉｎ〜ｋＬｍａｘと設定した後、｜Ｓｐ／μｐ｜＞ｋＬｍｉｎならば無次元特徴パラメータｐを非安定と判定し、
無次元特徴パラメータの非安定程度（以下、「無次元非安定程度」とよぶ）をパーセンテージで表す場合、｜Ｓｐ／μｐ｜≧ｋＬｍａｘのとき、無次元非安定程度を１００％とし、｜Ｓｐ／μｐ｜≦ｋＬｍｉｎのとき、無次元非安定程度を０％とし、｜Ｓｐ／μｐ｜＝ｋＬｍｉｎ〜ｋＬｍａｘのとき、無次元非安定程度を０％〜１００％とすると、無次元非安定程度が予め設定した閾値より大きい時に前記信号を非安定な信号と判定する。
無次元非定常程度算出法３：
ｘｊ（ｉ）における無次元特徴パラメータをｐｊのＵ個の最大値をｐＭｊ（ｊ＝１〜Ｕ）とし、ｐＭｊ除いた後のｘ（ｉ）の平均値と標準偏差をそれぞれμｐｎとＳｐｎとすると、予め係数ｋＭの範囲をｋＭｎｉｎ〜ｋＭｍａｘと設定した後、｜ｐＭｊ｜＞μｐｎ＋ｋＭｍｉｎＳｐｎならば無次元特徴パラメータｐを非安定と判定し、
無次元特徴パラメータの非安定程度（以下、「無次元非安定程度」とよぶ）をパーセンテージで表す場合、｜ｐＭｊ｜≧μｐｎ＋ｋＭａｘＳｐｎのとき、無次元非安定程度を１００％とし、｜ｐＭｊ｜≦μｐｎ＋ｋＭｉｎＳｐｎのとき、無次元非安定程度を０％とし、｜ｐＭｊ｜＝μｐｎ＋ｋＭｍｉｎＳｐｎ〜μｐｎ＋ｋＭａｘＳｐｎのとき、無次元非安定程度を０％〜１００％とすると、無次元非安定程度が予め設定した閾値より大きい時に前記信号を非安定な信号と判定する。
無次元非定常程度算出法４：
予め係数ｋｉの範囲をｋｉｎｉｎ〜ｋｉｍａｘと設定した後、｜ｐｊ×ｐｊ−１｜（ｊ＝２〜Ｍ）の最大値＞ｋｉｍｉｎならば無次元特徴パラメータｐを非安定と判定し、
無次元特徴パラメータの非安定程度（以下、「無次元非安定程度」とよぶ）をパーセンテージで表す場合、｜ｐｊ×ｐｊ−１｜（ｊ＝２〜Ｍ）の最大値≧ｋｉｍａｘのとき、無次元非安定程度を１００％とし、｜ｐｊ×ｐｊ−１｜（ｊ＝２〜Ｍ）の最大値≦ｋｉｍｉｎのとき、無次元非安定程度を０％とし、｜ｐｊ×ｐｊ−１｜（ｊ＝２〜Ｍ）の最大値＝ｋｉｍｉｎ〜ｋｉｍａｘのとき、無次元非安定程度を０％〜１００％とすると、無次元非安定程度が予め設定した閾値より大きい時に前記信号を非安定な信号と判定する。
【請求項２】
状態診断やパターン認識のために信号計測部により設定のサンプリング時間で信号を測定して波形データを得る第１工程と、請求項１に記載の信号検査モジュールを用いて信号検査を行う第２工程と、信号検査の結果を表示部で表示する第３工程と、前記信号検査の結果により前記波形データが状態診断やパターン認識に用いられるか否かの最終決定を行う第４工程と、前記最終決定の結果に従って信号の再測定または状態診断やパターン認識などを開始させる第５工程と、を有することを特徴とする信号検査方法。
【請求項３】
請求項２に記載の信号検査方法に基づいた信号検査機能を有することを特徴とする、状態診断やパターン認識のための装置。
【発明の詳細な説明】
【技術分野】
【０００１】本発明は、対象物の状態診断やパターン認識などの分野において、計測により得られた対象物の信号はオーバーフロー、外れ値および非定常性があるか否かを総合的に評価し、計測された信号の良否を検査する方法であり、診断装置やパターン認識装置などにおいて誤診断や誤判定を防ぐために提供できる信号検査用のモジュールに関するものである。
【背景技術】
【０００２】対象物の状態診断やパターン認識などのために信号を計測するときに、測定された信号が定常な信号か、非定常な信号かによって処理方法が違う。定常な信号を処理するための診断装置やパターン認識装置などに非定常な信号を入力すると、間違った結果が出力される。【０００３】また、診断やパターン認識などの対象物から発する真の信号が定常な信号であるが、計測ミスや外乱などの原因により計測された信号は非定常な信号となり、更にオーバーフローなどの発生により計測された信号は対象物の真の信号でない不備な信号となってしまう。このような非定常な信号や不備な信号を診断装置やパターン認識装置などに入力しても正確な結果が得られない。【０００４】従来、様々な分野で色々な診断装置やパターン認識装置などが提案されているが、診断やパターン認識のための処理を行う前に、測定された信号の良否（測定ミスや処理上の不備などの有無）についての検査手法およびモジュールが示されていない。この点については特に関連性のある過去の文献は先願［特許文献１］であり、その他に関連のある文献は、たとえば、［特許文献２，３，４，５、６］である。
【特許文献１】特願２００４−１６１８０２
【特許文献２】特許公開平７−３２４９７６
【特許文献３】特許公開２０００−１７１２９１
【特許文献４】特許公開平１１−２９０３２３
【特許文献５】特許公開２００５−３３５５９
【特許文献６】国際公開番号Ｗ００１／０７０１１３
【発明の開示】
【発明が解決しようとする課題】
【０００５】従来の診断装置やパターン認識装置には、測定された信号の測定ミスがあるか否かを判断する手段は、信号のオーバーフローの有無だけをチェックするか、または計測者の目視や勘などにより不備（測定ミス）の有無を判断することである。しかし、前記のように、定常な信号を処理するための診断装置やパターン認識装置は、測定された信号のオーバーフローの有無だけの確認では不十分であり、信号の定常性を検査する必要がある。また、計測者の目視や勘などにより不備の有無を判断するのは限界があり、特に測定された生信号が直接に確認できない計測装置の場合は判定の精度が保証できない。【０００６】本発明は、様々な分野で色々な診断装置やパターン認識装置などにおいて、診断やパターン認識のための処理を行う前に、測定された信号が定常な信号か、非定常な信号か、また測定された信号の測定ミスがあるか否かについての信号検査手法および信号検査モジュールを提供することを目的とする。
【課題を解決するための手段】
【０００７】上記の問題点を解決するために、本発明においては、次のような手段を採る。測定された信号が定常な信号か、非定常な信号かについては、有次元特徴パラメータおよび無次元特徴パラメータの定常性を確認することにより検査する。また、測定された信号に不備（測定ミス）があるか否かについては、オーバーフローの有無および外れ値の有無を確認することにより検査する。
【発明の効果】
【０００８】本発明においては、対象物の状態診断やパターン認識などのために測定された信号が定常な信号か、非定常な信号か、また不備があるか否かを自動的に検査する方法およびモジュールを提供できる。よって、本発明の方法およびモジュールは、様々な分野で色々な診断装置やパターン認識装置に適用でき、状態診断やパターン認識などの精度を高くすることができ、信号検査の自動化に役立つ。
【発明を実施するための最良の形態】
【０００９】本発明の処理の流れは図１に示す。ここで、この流れに沿って本発明の最良の形態について説明する。【００１０】信号の測定（図１のＡ）
ここで、対象とする信号は、振動信号、電圧信号、電流信号、音響信号および生体信号などであり、対象物の状態診断やパターン認識に用いられるものである。【００１１】時系列信号（図１のＢ）
本発明の検査対象は時系列信号であり、図５，図６，図７および図８の例に示すように、計測部で測定されて状態診断やパターン認識などの処理をしようとする対象信号の波形データ全体である。すなわち、測定された信号は不備などがあるか否か、また、時系列信号の統計的な性質が時間と共に変化するか否かを検査することにより信号の定常性を判定する。なお、計測された信号はサンプリングによりＮ個の離散データ波形となり、ｘ（ｉ）（ｉ＝１〜Ｎ）で表す。なお、ｘ（ｉ）の平均値μｒと標準偏差Ｓｒを求めたら、次のようにｘ（ｉ）の正規化を行えば、
0005017678.tif 000002
【００１２】
【数１】
0005017678.tif 000003
【００１３】オーバーフロー検出法（図１のＣ，Ｄ）
オーバーフローとは、信号の最大値（プラス側）が計測装置の最大計測レンジより大きいか、最小値（マイナス側）が計測装置の最大計測レンジより小さいことである。オーバーフローが発生したか否かは次のように判定する。波形データｘ（ｉ）の最大値と最小値をそれぞれｘｍａｘ（ｉ）とｘｍｉｎ（ｉ）とし、計測装置の最大計測レンジと最小計測レンジをそれぞれＲｍａｘとＲｍｉｎとすると、ｘｍａｘ（ｉ）≧Ｒｍａｘはるいはｘｍｉｎ（ｉ）≦Ｒｍｉｎならば、オーバーフローが発生したと判定する。【００１４】外れ値の検出法（図１のＥ，Ｆ，Ｇ，Ｈ）
外れ値は、「統計においてデータの集団から外れた極端に大きな値あるいは極端に小さな値である」と定義されている。外れ値の検定方法は、様々あるが、本発明では、処理の迅速化のために波形データの中に「平均値＋４×標準偏差」以上、あるいは「平均値−４×標準偏差」以下の値を外れ値と見做し、外れ値の程度は「平均値＋４×標準偏差」より大きいほど、あるいは「平均値−４×標準偏差」より小さいほど、外れ値の発生程度が大きい。外れ値の検出は次の２つの方法がある。
（１）波形データの平均値と標準偏差を用いた外れ値検出法
状態診断やパターン認識などのために得られた信号の離散データをｘ（ｉ）（ｉ＝１〜Ｎ）とし、ｘ（ｉ）の絶対値を｜ｘ（ｉ）｜とし、｜ｘ（ｉ）｜の平均値と標準偏差をそれぞれμとＳとし、係数ｋの範囲をｋｍｉｎ〜ｋｍａｘと設定した後、｜ｘ（ｉ）｜の中にある値｜ｘ（ｊ）｜（ｊ＝１〜Ｎ）＞μ＋ｋｍｉｎＳならば、ｘ（ｊ）を外れ値と判定される。一般にｋｍｉｎ〜ｋｍａｘ＝４〜５０である。望ましいｋｍｉｎ〜ｋｍａｘは５〜２０である。【００１５】（２）波形データの絶対最大値を用いた外れ値検出法
前記のｘ（ｉ）のＬ個の絶対最大値を｜ｘＭｊ｜（ｊ＝１〜Ｌ）とし、ｘＭｊ除いた後の平均値と標準偏差をそれぞれμａとＳａとすると、係数ｋａの範囲をｋａｍｉｎ〜ｋａｍａｘと設定した後、｜ｘＭｊ｜＞μａ＋ｋｍｉｎＳａならば、ｘＭｊを外れ値と判定される。ここで、一般にｋａｍｉｎ〜ｋａｍａｘ＝４〜３０である。望ましいｋａｍｉｎ〜ｋａｍａｘは５〜２０である。また、一般にＬ／Ｎ＝０．００００１〜０．１であり、望ましいＬ／Ｎの範囲は０．０００１〜０．００１である。【００１６】（３）外れ値の発生程度の計算法
外れ値の発生程度（以下「外れ程度」と呼ぶ）をパーセンテージで表す場合、係数ｋ（または、ｋａ）の上限値と下限値をそれぞれｋｍａｘとｋｍｉｎ（または、ｋａｍａｘとｋａｍｉｎ）とし、｜ｘ（ｊ）｜≧μ＋ｋｍａｘＳ（または、｜ｘＭｊ｜≧μａ＋ｋａｍａｘＳａ）のとき、外れ程度を１００％とし、｜ｘ（ｊ）｜≦μ＋ｋｍｉｎＳ（または、｜ｘＭｊ｜≦μａ＋ｋａｍｉｎＳａ）のとき、外れ程度を０％とし、｜ｘ（ｊ）｜＝μ＋ｋｍｉｎＳ〜μ＋ｋｍａｘ（または、｜ｘＭｊ｜＝μａ＋ｋａｍｉｎＳａ〜μａ＋ｋａｍａｘＳａ）のとき、外れ程度を０％〜１００％とすると、外れ程度が予め設定した閾値より大きい時に前記信号を不備な信号と判定して適宜な対処を行う。
たとえば、ｋｍａｘ＝１５とｋｍｉｎ＝５のとき、ｘｍａｘ＝μ＋５Ｓであれば０％、ｘｍａｘ＝μ＋１０Ｓであれば５０％、ｘｍａｘ＞μ＋１５Ｓであれば１００％。
外れ値の発生程度はｘ％以上の場合、信号の不備や測定ミスがあると判定されるが、ｘ％は「外れパーセンテージ閾値」といって、測定対象の性質によって測定者が決める。また、外れ値の発生程度が大きいと判明した後、その対処として、信号を測定し直すか、外れ値を取り除いて処理するか測定者が決めることができる。【００１７】有次元特徴パラメータの定常性について（図１のＩ，Ｊ，Ｋ）
有次元特徴パラメータとは、波形データから求めた単位付きのパラメータのことである。波形データｘ（ｉ）の有次元特徴パラメータの例としては、次のものがある。
0005017678.tif 000004
なお、上記の有次元特徴パラメータ以外に実効値もあるが、実効値の定常性については標準偏差と同様に考えることができる。【００１８】有次元特徴パラメータの定常性を検査するとき、波形データｘ（ｉ）（あるいは、
0005017678.tif 000005
いは、ｘｊ（ｉ）とｘｋ（ｉ））で表す。一般にＭ／Ｎ＝０．０００１＝０．５であり、望ましいＭ／Ｎの範囲は０．０１〜０．１である。第ｊ区間と第ｋ区間の波形データ数をＮｊとＮｋとする。
0005017678.tif 000006
【００１９】（１）有次元特徴パラメータの定常性評価法Ｉ
μの定常性を検定するときに、ｘｊ（ｉ）とｘｋ（ｉ）における平均値をμｊとμｋとし、ｘｊ（ｉ）とｘｋ（ｉ）における標準偏差をＳｊとＳｋとすると、μｊ＝μｋならば、第ｊ区間と第ｋ区間はμに関する定常区間と定義する。Ｍ個の区間における平均値μが全て等しければ、すなわちμ１＝μ２＝・・・μｊ＝μｋ＝・・・＝μＭならば、波形データ（信号）の平均値が定常であると判定できる。しかし、実際にはデータ波形のばらつきを考慮して、確率の有意水準αを与えて、仮説μ１＝μ２＝・・・μｊ＝μｋ＝・・・＝μＭについての統計検定を行う。統計検定方法［非特許文献１］は多く提案されているが、ここでその１例を示す。たとえば、有意水準αが与えられた場合、【００２０】
【数２】
0005017678.tif 000007
であるば、上記の仮説が成り立たなく、信号が非定常と判定される。ここで、ｊ≠ｋ，ｊ，ｋ＝１〜Ｍで、ｔ（α，∞）は自由度∞のｔ分布の確率密度関数が上側確率αに対するパーセント点である。αの最小値と最大値をそれぞれαｍｉｎとαｍａｘとすると、一般にαｍｉｎ〜αｍａｘ＝０．００００１〜０．９であるが、望ましい値としてαｍｉｎ〜αｍａｘ＝０．００１〜０．５である。
0005017678.tif 000008
る。
【非特許文献１】Ｋ．Ａ．Ｂｒｏｗｎｌｅｅ：Ｓｔａｔｉｓｔｉｃａｌ  Ｔｈｅｏｒｙ  ａｎｄ  Ｍｅｔｈｏｄｏｌｏｇｙ  ｉｎ  Ｓｃｉｅｎｃｅ  ａｎｄ  Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｅｃｏｎｄ  Ｅｄｉｔｉｏｎ，Ｔｈｅ  Ｕｎｉｖｅｒｓｉｔｙ  ｏｆ  Ｃｈｉｃａｇｏ，１９６５．
【００２１】（２）有次元特徴パラメータの定常性評価法ＩＩ
0005017678.tif 000009
【００２２】
【数３】
0005017678.tif 000010
であるば、信号が非定常と判定される。ここで、ｔ（α／２，Ｎｊ−１）は自由度（Ｎｊ−１）のｔ分布の確率密度関数が上側確率α／２に対するパーセント点である。αの最小値と最大値をそれぞれαｍｉｎとαｍａｘとすると、一般にαｍｉｎ〜αｍａｘ＝０．００００１〜０．９であるが、望ましい値としてαｍｉｎ〜αｍａｘ＝０．００１〜０．５である。【００２３】（３）有次元特徴パラメータの定常性評価法ＩＩＩ
Ｓの定常性を検定するとき、ｘｊ（ｉ）とｘｋ（ｉ）における標準偏差をＳｊとＳｋとすると、【００２４】
【数４】
0005017678.tif 000011
であるば、Ｓｊ≠Ｓｋであり、信号が非定常と判定される。ここで、ここで、Ｆ（α／２，Ｎｊ−１，Ｎｋ−１）は自由度（Ｎｊ−１，Ｎｋ−１）のＦ分布の確率密度関数が上側確率α／２に対するパーセント点である。αの最小値と最大値をそれぞれαｍｉｎとαｍａｘとすると、一般にαｍｉｎ〜αｍａｘ＝０．００００１〜０．９であるが、望ましい値としてαｍｉｎ〜αｍａｘ＝０．００１〜０．５である。
0005017678.tif 000012
る。【００２５】有次元特徴パラメータの定常性評価法ＩＶ
（４）有次元特徴パラメータの定常性評価法ＩＶ
0005017678.tif 000013
と、【００２６】
【数５】
0005017678.tif 000014
であるば、信号が非定常と判定される。ここで、Ｘ２（α／２，Ｎｊ−１）は自由度（Ｎｊ−１）のカイ２乗分布の確率密度関数が上側確率α／２に対するパーセント点である。αの最小値と最大値をそれぞれαｍｉｎとαｍａｘとすると、一般にαｍｉｎ〜αｍａｘ＝０．００００１〜０．９であるが、望ましい値としてαｍｉｎ〜αｍａｘ＝０．００１〜０．５である。【００２７】（５）有次元特徴パラメータの非定常程度の計算法
有次元特徴パラメータの非安定程度（以下、「有次元非安定程度」とよぶ）をパーセンテージで表す場合、、有意水準αｄを与えたときに、［数２］、［数３］、［数４］および［数５］が左側と右側が等しく、しかもαｄ≧αｍａｘのときに、有次元非安定程度を０％とし、また、αｄ≦αｍｉｎのときに、有次元非安定程度を１００％とし、αｄ＝αｍｉｎ〜αｍａｘのとき、有次元非安定程度を１００％〜０％とすると、有次元非安定程度が予め設定した閾値より大きい時に前記信号を非安定な信号と判定して適宜な対処を行う。
たとえば、［数２］において、αｍａｘ＝０．２、αｍｉｎ＝０．００１とする。α＝０．３のときに、［数２］の左側と右側が等しければ０％、α＝０．０００５のとき、［数２］の左側と右側が等しければ１００％、α＝０．１のとき、左側と右側が等しければ５０．３％、α＝０．０５のとき、左側と右側が等しければ２５．１％である。
有次元非安定程度はｘ％以上の場合、信号が「非定常」と判定すれば、ｘ％は有次元特徴パラメータの「非定常パーセンテージ閾値」といって、測定対象の性質によって測定者が決める。【００２８】信号が非定常と判定されたとき、信号を測定しなおすか否かは対象物の性質によって決められる。場合によって検定の結果を測定者に示し、測定者が信号を測定しなおすか否かを決める。
なお、定常な信号しか処理できない診断装置やパターン認識装置などには、非定常な信号を入力することは明らかに不適切である。【００２９】無次元特徴パラメータの定常性評価法について（図１のＬ，Ｍ，Ｎ）
無次元特徴パラメータとは、波形データから求めた単位無しのパラメータのことである。波形データｘ（ｉ）（＝ｘｉ）の無次元特徴パラメータの例としては、次のものがある。【００３０】
【数６】
0005017678.tif 000015
ここで、【００３１】
【数７】
0005017678.tif 000016
【００３２】
【数８】
0005017678.tif 000017
【００３３】
【数９】
0005017678.tif 000018
【００３４】ここで、μはｘｉの平均値で、σはｘｉの標準偏差である。【００３５】
【数１０】
0005017678.tif 000019
【００３６】
【数１１】
0005017678.tif 000020
ここで，μｐは波形の極大値（ピーク値）の平均値である。【００３７】
【数１２】
0005017678.tif 000021
ここで、μｍａｘはの１０個最大値の平均値である。【００３８】
【数１３】
0005017678.tif 000022
ここで、σｐは極大値の標準偏差値である。【００３９】
【数１４】
0005017678.tif 000023
ここで、μＬとσＬはそれぞれ極小値（谷値）の平均値と標準偏差値である。【００４０】
【数１５】
0005017678.tif 000024
【００４１】
【数１６】
0005017678.tif 000025
【００４２】
【数１７】
0005017678.tif 000026
【００４３】
【数１８】
0005017678.tif 000027
但し、ｘｋｉ＞ｋσ（例えば、ｋ＝１，２）。【００４４】
【数１９】
0005017678.tif 000028
但し、ｘｈｉ＞ｋσ。（例えば、Ｈ＝１，２）。
なお、以下で述べる無次元特徴パラメータの定常性の判定法は上記以外に他の無次元特徴パラメータにも適用できる。【００４５】無次元特徴パラメータの定常性を検査するとき、波形データｘ（ｉ）をＭ区間に分割し、第ｊ区間の波形データをｘｊ（ｉ）で表す。一般にＭ／Ｎ＝０．０００１〜０．１であり、望ましいＭ／Ｎの範囲は０．００５〜０．０１である。第ｊ区間の波形データ数をＮｊとする。舞次元特徴パラメータを統一にｐで表し、ｘｊ（ｉ）における無次元特徴パラメータｐをｐｊ（ｊ＝１〜Ｍ）で表す。
無次元特徴パラメータｐの定常性を判定するとき、次のような方法がある。【００４６】（１）無次元特徴パラメータの定常性評価法Ｉ
予め係数ｋｐの範囲をｋｐｎｉｎ〜ｋｐｍａｘと設定した後、｜ｐｊ｜＞μｐ＋ｋｐｍｉｎＳｐならばｐｊを特異な無次元特徴パラメータと判定され、信号が非定常と判定される。一般にｋｐｎｉｎ〜ｋｐｍａｘ＝１〜２０である。望ましいｋｐｎｉｎ〜ｋｐｍａｘは２〜３である。【００４７】（２）無次元特徴パラメータの定常性評価法ＩＩ
予め係数ｋＬの範囲をｋＬｎｉｎ〜ｋＬｍａｘと設定した後、｜Ｓｐ／μｐ｜＞ｋＬｍｉｎならば無次元特徴パラメータｐを非安定と判定される。一般にｋＬｎｉｎ〜ｋＬｍａｘ＝０．１〜６である。望ましいｋＬｎｉｎ〜ｋＬｍａｘは０．５〜３である。【００４８】（３）無次元特徴パラメータの定常性評価法ＩＩＩ
ｘｊ（ｉ）における無次元特徴パラメータをｐｊのＵ個の最大値をｐＭｊ（ｉ＝１〜Ｕ）とし、ｐＭｊ除いた後のｘ（ｉ）の平均値と標準偏差をそれぞれμｐｎとＳｐｎとすると、予め係数ｋＭの範囲をｋＭｎｉｎ〜ｋＭｍａｘと設定した後、｜ｐＭｊ｜＞μｐｎ＋ｋＭｍｉｎＳｐｎならば無次元特徴パラメータｐを非安定と判定され、信号が非定常と判定される。一般にｋＭｎｉｎ〜ｋＭｍａｘ＝２〜１０である。望ましいｋＭｎｉｎ〜ｋＭｍａｘは３〜６である。また、一般にＵ／Ｍ＝０．００１〜０．５であり、望ましいＵ／Ｍの範囲は０．００１〜０．０１である。【００４９】（４）無次元特徴パラメータの定常性評価法ＩＶ
予め係数ｋｉの範囲をｋｉｎｉｎ〜ｋｉｍａｘと設定した後、｜ｐｊ×ｐｊ−１｜（ｊ＝２〜Ｍ）の最大値＞ｋｉｍｉｎならば無次元特徴パラメータｐを非安定と判定される。ｋｉの決め方は経験
0005017678.tif 000029
【００５０】（５）無次元特徴パラメータの非定常程度の計算法
無次元特徴パラメータの非定常程度（以下、「無次元非安定程度」とよぶ）をパーセンテージで表す場合、係数ｋｐ（ｋＬ、ｋＭ、ｋｉも同様）の上限値と下限値をそれぞれｋｍａｘとｋｍｉｎとし、｜ｐｊ｜≧μｐ＋ｋｐｍａｘＳｐのとき、無次元非安定程度を１００％とし、、｜ｐｊ｜≦μｐ＋ｋｐｍｉｎＳｐのとき、無次元非安定程度を０％とし、ｋｐ＝μｐ＋ｋｐｍｉｎＳｐ〜ｋｐｍａｘＳｐのとき、無次元非安定程度を０％〜１００％とすると、無次元非安定程度が予め設定した閾値より大きい時に前記信号を非安定な信号と判定して適宜な対処を行う。
たとえば、ｋｐｍａｘ＝５とｋｐｍｉｎ＝１のとき、｜ｐｊ｜＝μｐ＋０．８Ｓｐであれば０％、｜ｐｊ｜＝ｕｐ＋３Ｓｐであれば５０％、｜ｐｊ｜＝μｐ＋６Ｓｐであれば１００％。
また、係数ｋｉ（なお、以下の方法はｋＬにも適用できる。）の上限値と下限値をそれぞれｋｉｍａｘとｋｉｍｉｎとし、｜ｐｊｐｊ−１｜（ｊ＝２〜Ｍ）の最大値をΔｐｊｍａｘとすると、Δｐｊｍａｘ≧ｋｉｍａｘのとき、１００％とし、Δｐｊｍａｘ＜ｋｉｍｉｎのとき、０％とする。
たとえば、ｋｉｍａｘ＝１１とｋｉｍａｘ＝１とのとき、Δｐｊｍａｘ＝０．９であれば０％、Δｐｊｍａｘ＝６であれば５０％、Δｐｊｍａｘ＝１２であれば１００％である。
無次元特徴パラメータの非定常程度はｘ％以上の場合、信号が「非定常」と判定すれば、ｘ％は無次元特徴パラメータの「非定常パーセンテージ閾値」といって、測定対象の性質によって測定者が決める。【００５１】信号が非定常と判定されたとき、信号を測定しなおすか否かは対象物の性質によって決められる。場合によって検定の結果を測定者に示し、測定者が信号を測定しなおすか否かを決める。
なお、定常な信号しか処理できない診断装置やパターン認識装置などには、非定常な信号を入力することは明らかに不適切である。【００５２】以上、各項目による信号検査の結果を信号の測定者に示し、信号を計測しなおすか否かは測定者によって決めることができる。複数の有・無次元特徴パラメータで評価する場合、得られた複数の有・無次元特徴パラメータの非定常程度（％）の最大値を最終結果として表示する。
図２は棒グラフによる信号検査結果の表示法（Ｉ、ＩＩ）を示す。
図３はランプによる信号検査結果の表示法（Ｉ、ＩＩ、ＩＩＩ、ＩＶ）を示す。ランプの点灯か消灯かは対象物の特性により決められる。望ましい決め方は次の通りである。
ランプによる信号検査結果の表示法ＩとＩＩＩにおいては、８０％以上は「赤」、３０％以上〜８０％以下は「黄」、３０％未満は「青」とする。ランプによる信号検査結果の表示法ＩＩとＩＶにおいては、５０％以上は「点灯」、５０％未満は「消灯」とする。【００５３】更に、ブザーなどの音による表示も出来る。信号検査結果（％）が設定されたパーセンテージ以上になった場合、ブザーを鳴らす。たとえば、５０％以上は「ブザー音」とする。【００５４】図４は本発明の方法により構築される信号検査モジュールの構成図を示す。すなわち、対象物の状態診断やパターン認識などのために、信号計測部で測定された信号が定常な信号か、非定常な信号か、また測定された信号の測定ミスがあるか否かを信号検査部で自動的に検査し、検査結果を検査結果の表示部で表示し、診断部やパターン認識部で診断やパターン認識のための処理を行う前に、信号の測定者が検査結果に基づき信号を測定しなおすか否かをその場で決定するか、事前に設定したルールにより信号検査部で自動的に対処することができ、状態診断やパターン認識などの精度を確保することができる。
【実施例１】
【００５５】図５は設備診断装置によって測定された回転機械の振動加速度信号である。この設備診断装置は、定常回転のモータ、ポンプ、発電機および歯車装置などの回転機械の状態診断に用いられる。この診断装置は診断対象の回転機械の定常な振動信号を測定して、定常な振動信号を処理することにより状態診断を行う。
この信号を検査した結果を次に示す。
データの測定条件：データ数：Ｎ＝２０４８、サンプリング周波数：１０００Ｈｚ。
オーバーフローの検査結果：オーバーフロー検出のレンジはＲｍａｘ＝４．０（ｖ）、Ｒｍｉｎ＝−４．０（ｖ）である。この信号のｘｍａｘ（ｉ）＝２．７９（ｖ）、ｘｍｉｎ（ｉ）＝−２．９５（ｖ）である。よって、ｘｍａｘ（ｉ）≦Ｒｍａｘ、ｘｍｉｎ（ｉ）≧Ｒｍｉｎから、オーバーフローは無し。
外れ値の検査結果：外れ値検出法Ｉの係数ｋをｋｍａｘ＝２０とｋｍｉｎ＝５とする。この信号のμ＝０．０１３７８４、Ｓ：０．７９９である。よって、｜ｘ（ｉ）｜の最大値２．７９＜０．０１３７８４＋５×０．７９９から、外れ値の程度は０％である。
有次元特徴パラメータ（μｒとμａ）の定常性の判定結果：
波形データの分割数Ｍ＝８、
有次元特徴パラメータの定常性評価法Ｉにおける有意水準α＝０．００１、
有意水準の最小値αｍｉｎ＝０．００１、
有意水準の最大値αｍｉｎ＝０．０１、
ｔ（０．００１，２０４８／８）＝３．３２９８６７、
ｔ（０．０１，２０４８／８）＝２．５９６、
第１分割区間のｘ１（ｉ）の平均値μ１および｜ｘ１（ｉ）｜の平均値μａ１を基準として、有次元特徴パラメータの定常性評価法Ｉにより求めた分割区間のｔｉおよびｔａｉによる判定結果は次の通りである。
第２分割区間のｔ２＝１．４５、非定常程度＝０％
第２分割区間のｔａ２＝０．８０、非定常程度＝０％
第３分割区間のｔ３＝０．２３、非定常程度＝０％
第３分割区間のｔａ３＝０．９７、非定常程度＝０％
第４分割区間のｔ４＝１．３５、非定常程度＝０％
第４分割区間のｔａ４＝０．８７、非定常程度＝０％
第５分割区間のｔ５＝０．７９、非定常程度＝０％
第５分割区間のｔａ５＝０．１０、非定常程度＝０％
第６分割区間のｔ６＝２．３６、非定常程度＝０％
第６分割区間のｔａ６＝１．０８、非定常程度＝０％
第７分割区間のｔ７＝０．４０、非定常程度＝０％
第７分割区間のｔａ７＝０．２５、非定常程度＝０％
第８分割区間のｔ８＝１．９５、非定常程度＝０％
第８分割区間のｔａ８＝１．０７、非定常程度＝０％
以上の結果により、この波形は有次元特徴パラメータの非定常の程度は０％である。
無次元特徴パラメータ（ｐ２とｐ３）の定常性の判定結果：
波形データの分割数Ｍ＝８、
無次元特徴パラメータの定常性評価法ＩＶにおけるｋｊは
ｐ２の係数ｋ２ｍａｘ＝２、ｋ２ｍｉｎ＝１
ｐ３の係数ｋ３ｍａｘ＝３、ｋ３ｍｉｎ＝１
ｐ２のΔｐｊｍａｘ＝０．８９、非定常程度＝０％
ｐ３のΔｐｊｍａｘ＝１．２３、非定常程度＝１１．５％
以上の結果により、この波形は無次元特徴パラメータの非定常の程度が１１．５％である。
よって、この信号は正常に測定された信号として、状態診断に用いられた。
【実施例２】
【００５６】図６は図５の波形と同じ診断装置と測定条件で測定された波形例である。この信号を検査した結果を次に示す。
データの測定条件：データ数：Ｎ＝２０４８、サンプリング周波数：１０００Ｈｚ。
オーバーフローの検査結果：オーバーフロー検出のレンジはＲｍａｘ＝４．０（ｖ）、Ｒｍｉｎ＝−４．０（ｖ）である。この信号のｘｍａｘ（ｉ）＝３．１４（ｖ）、ｘｍｉｎ（ｉ）＝−３．１９（ｖ）である。よって、ｘｍａｘ（ｉ）≦Ｒｍａｘ、ｘｍｉｎ（ｉ）≧Ｒｍｉｎから、オーバーフローは無し。
外れ値の検査結果：外れ値検出法Ｉの係数ｋをｋｍａｘ＝２０とｋｍｉｎ＝５とする。この信号のμ＝０．２６、Ｓ＝０．７０１である。よって、｜ｘ（ｉ）｜の最大値３．１９＜０．２６＋５×０．７０１から、外れ値の程度は０％である。
有次元特徴パラメータ（μｒとμａ）の定常性の判定結果：
波形データの分割数Ｍ＝８、
有次元特徴パラメータの定常性評価法Ｉにおける有意水準α＝０．００１、
有意水準の最小値αｍｉｎ＝０．００１、
有意水準の最大値αｍｉｎ＝０．０１、
ｔ（０．００１，２０４８／８）＝３．３２９８６７、
ｔ（０．０１，２０４８／８）＝２．５９６、
第１分割区間のｘ１（ｉ）の平均値μ１および｜ｘ１（ｉ）｜の平均値μａ１を基準として、有次元特徴パラメータの定常性評価法Ｉにより求めた分割区間のｔｉおよびｔａｉによる判定結果は次の通りである。
第２分割区間のｔ２＝４．９８、非定常程度＝１００％
第２分割区間のｔａ２＝０．７８、非定常程度＝０％
第３分割区間のｔ３＝５．３２、非定常程度＝１００％
第３分割区間のｔａ３＝１．９１、非定常程度＝０％
第４分割区間のｔ４＝４．４８、非定常程度＝１００％
第４分割区間のｔａ４＝０．８２、非定常程度＝０％
第５分割区間のｔ５＝４．９２、非定常程度＝１００％
第５分割区間のｔａ５＝０．１０、非定常程度＝０％
第６分割区間のｔ６＝０．４４、非定常程度＝０％
第６分割区間のｔａ６＝０．７６、非定常程度＝０％
第７分割区間のｔ７＝２．７５、非定常程度＝５％
第７分割区間のｔａ７＝０．３３、非定常程度＝０％
第８分割区間のｔ８＝６．７２、非定常程度＝１００％
第８分割区間のｔａ８＝０．９９、非定常程度＝０％
以上の結果により、この波形は有次元特徴パラメータの非定常の程度が１００％である。
無次元特徴パラメータ（ｐ２とｐ３）の定常性の判定結果：
波形データの分割数Ｍ＝８、
無次元特徴パラメータの定常性評価法ＩＶにおけるｋｊは
ｐ２の係数ｋ２ｍａｘ＝２、ｋ２ｍｉｎ＝１
ｐ３の係数ｋ３ｍａｘ＝３、ｋ３ｍｉｎ＝１
ｐ２のΔｐｊｍａｘ＝１．６８、非定常程度＝６８％
ｐ３のΔｐｊｍａｘ＝１．７５、非定常程度＝３７．５％
以上の結果により、この波形は有次元特徴パラメータの非定常の程度が６８％である。
この例では、有次元特徴パラメータの非定常の程度が１００％であるから、定常信号しか処理できない上記の設備診断装置は、この信号を処理できないと判定された。
【実施例３】
【００５７】図７も図５の波形と同じ診断装置と測定条件で測定された波形例である。この信号を検査した結果を次に示す。
データの測定条件：データ数：Ｎ＝２０４８、サンプリング周波数：１０００Ｈｚ。
オーバーフローの検査結果：オーバーフロー検出のレンジはＲｍａｘ＝２．０（ｖ）、Ｒｍｉｎ＝−２．０（ｖ）である。この信号のｘｍａｘ（ｉ）＝３．１７（ｖ）、ｘｍｉｎ（ｉ）＝−２．７９（ｖ）である。よって、ｘｍａｘ（ｉ）＞Ｒｍａｘ、ｘｍｉｎ（ｉ）＜Ｒｍｉｎから、オーバーフローが発生したと判定され、信号を測定しなおした。
【実施例４】
【００５８】図８も図５の波形と同じ診断装置と測定条件で測定された波形例である。この信号を検査した結果を次に示す。
データの測定条件：データ数：Ｎ＝２０４８、サンプリング周波数：１０００Ｈｚ。
オーバーフローの検査結果：オーバーフロー検出のレンジはＲｍａｘ＝１５．０（ｖ）、Ｒｍｉｎ＝−１５．０（ｖ）である。この信号のｘｍａｘ（ｉ）＝９．７（ｖ）、ｘｍｉｎ（ｉ）＝−１４．０（ｖ）である。よって、ｘｍａｘ（ｉ）≦Ｒｍａｘ、ｘｍｉｎ（ｉ）≧Ｒｍｉｎから、オーバーフローは無し。
外れ値の検査結果：外れ値検出法ＩＩのＬをＬ＝２、係数ｋをｋｍａｘ＝１２とｋｍｉｎ＝６とする。この信号のμ＝−０．００５４、Ｓ＝０．８２２である。よって、｜ｘ（ｉ）｜の最大値１２．３＞−０．００５４＋１２×０．８２２から、外れ値の程度は１００％であり、信号を測定しなおした。
【図面の簡単な説明】
【００５９】
【図１】本発明の処理の流れを示すグラフである。
【図２】棒グラフによる信号検査結果の表示法を示すグラフである。
【図３】ランプによる信号検査結果の表示法を示すグラフである。
【図４】本発明の信号検査モジュールの構成を示すグラフである。
【図５】定常な信号の例を示すグラフである。
【図６】非定常な信号の例を示すグラフである。
【図７】オーバーフローが発生した信号の例を示すグラフである。
【図８】外れ値が存在する信号の例を示すグラフである。
【００６０】図４中の符号について、
１  状態診断装置やパターン認識装置、２  本発明の信号検査モジュール、３  信号計測部、４  信号検査部、５  検査結果表示部、６  診断部やパターン認識部。
【図１】
0005017678.tif 000030
【図２】
0005017678.tif 000031
【図３】
0005017678.tif 000032
【図４】
0005017678.tif 000033
【図５】
0005017678.tif 000034
【図６】
0005017678.tif 000035
【図７】
0005017678.tif 000036
【図８】
0005017678.tif 000037