JP2000275098A - 波形信号解析装置 - Google Patents

波形信号解析装置

Info

Publication number
JP2000275098A
JP2000275098A JP11084781A JP8478199A JP2000275098A JP 2000275098 A JP2000275098 A JP 2000275098A JP 11084781 A JP11084781 A JP 11084781A JP 8478199 A JP8478199 A JP 8478199A JP 2000275098 A JP2000275098 A JP 2000275098A
Authority
JP
Japan
Prior art keywords
waveform signal
wavelet transform
cycle
unit
wavelet
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.)
Granted
Application number
JP11084781A
Other languages
English (en)
Other versions
JP2000275098A5 (ja
JP3730435B2 (ja
Inventor
Keiichi Sasaki
恵 一 佐々木
Akira Sawada
田 彰 澤
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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP08478199A priority Critical patent/JP3730435B2/ja
Priority to US09/533,937 priority patent/US6542836B1/en
Publication of JP2000275098A publication Critical patent/JP2000275098A/ja
Publication of JP2000275098A5 publication Critical patent/JP2000275098A5/ja
Application granted granted Critical
Publication of JP3730435B2 publication Critical patent/JP3730435B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

(57)【要約】 【課題】 波形信号の周期性を見出し、当該周期でウェ
ーブレット変換結果を分割し、分割された波形信号のウ
ェーブレット変換結果の時系列的な変動を容易に比較す
ることができる波形信号解析装置を提供すること。 【解決手段】 本発明の波形信号解析装置は、波形信号
を入力する入力手段10と、入力手段10により入力さ
れた波形信号の周期を検出する周期検出手段20と、周
期検出手段20で検出した周期によって波形信号を分割
する分割手段30とを備える。ウェーブレット変換手段
40は、分割された個々の波形信号をウェーブレット変
換する。表示手段50は、分割された個々の波形信号の
ウェーブレット変換結果を時系列に連続表示する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、波形信号にウェー
ブレット変換を施して時間−周波数解析を行う波形信号
解析装置に係り、とりわけ、波形信号の周期性を検出し
て当該周期で波形信号を分割し、分割された波形信号の
ウェーブレット変換結果を容易に比較することができる
波形信号解析装置に関する。
【0002】
【従来の技術】従来の有力な波形信号解析装置として、
波形信号にFFT(高速フーリエ変換)を施すFFTア
ナライザが知られている。これは、入力波形に含まれる
周波数成分をフーリエスペクトル強度として解析表示す
る装置であり、フーリエスペクトルの分布から、入力波
形に異常な周波数成分が含まれているかどうかをチェッ
クすることができる。
【0003】但し、FFTは、変換結果に時間(位相)
情報を持たない解析方法であるため、入力波形の解析対
象時間内に含まれる同じ周波数成分毎の合計値の分布を
示すに止まり、各周波数成分の発生時刻情報を得ること
はできない。
【0004】従来において、周波数成分の時系列の変化
を知りたい場合には、窓関数を用いた短時間フーリエ変
換によって、解析対象波形を微小時間だけ切り出し各時
刻におけるフーリエスペクトルを時間軸方向に時系列に
並べて表示させることで、時間−周波数解析を実現させ
ている。
【0005】また近年、短時間フーリエ変換と同様に、
時間情報と周波数情報の両方を取り扱う事ができる時間
−周波数解析の手法として、ウェーブレット変換が提案
され、普及が進んでいる。
【0006】ウェーブレット変換は、マザーウェーブレ
ットと呼ばれる局在関数のスケールパラメータ(周波
数)とシフトパラメータ(時間)とを独立に変えなが
ら、対象波形との内積をとり、対象波形を時間−周波数
スペクトルに変換する信号処理方法である。
【0007】ウェーブレット変換の利点は、対象波形内
に存在する様々な周波数成分について、その発生時刻と
強度との両方を求めることができることである。
【0008】
【発明が解決しようとする課題】ウェーブレット変換の
変換結果の表示方法については、時間−周波数分布を、
例えば時間−周波数平面に強度軸を立てた3次元表示
や、時間−周波数平面上で強度を段階的に色分けした等
高線図状の表示とすることができる。このような表示
は、対象波形の特徴成分の検出、判別のために非常に有
効である。
【0009】ここで、時間−周波数解析の対象とする入
力波形が、電源波形などの何らかのパイロット信号に同
期した特徴をもつ場合には、この信号の周期で入力波形
を分割してウェーブレット変換することで、パイロット
信号の位相と解析対象波形との相関を明確に調べること
ができる。従って、パイロット信号の周期を一区間(一
単位)の解析時間範囲として選択して時間−周波数解析
をすることが有効である。
【0010】そして、解析対象波形内の異常を検出した
い場合は、前記のように連続的に分割した解析対象波形
をそれぞれウェーブレット変換し、そのウェーブレット
変換結果を視覚化して見比べる方法が有効である。視覚
化の具体的な手法としては、変換結果(グラフ)のディ
スプレイ上での並列表記や、印刷物として並べること等
が考えられる。
【0011】しかし、このような視覚化の態様では、ウ
ェーブレット変換結果の特徴ピークの時系列の連続的変
動を捉えることが極めて困難である。また、長時間の時
系列変化を連続的に比較する場合には、連続的に結果を
並べるための極めて広いスペースが必要である。
【0012】また、予めパイロット信号の周期が分かっ
ていない場合には、ウェーブレット変換による時間−周
波数解析の対象波形を適切に分割することができず、ウ
ェーブレット変換結果における特徴ピークの時系列の挙
動を調べることが困難である。
【0013】本発明は、このような点を考慮してなされ
たものであり、入力波形信号の周期性を見出して当該周
期で波形信号を分割し、分割された波形信号のウェーブ
レット変換結果の時系列的な変動を容易に比較すること
ができる波形信号解析装置を提供することを目的とす
る。
【0014】あるいは本発明は、ウェーブレット変換結
果から波形信号の周期性を見出し、当該周期でウェーブ
レット変換結果を分割して、ウェーブレット変換結果の
時系列的な変動を容易に比較することができる波形信号
解析装置を提供することを目的とする。
【0015】
【課題を解決するための手段】本発明は、波形信号を入
力する入力手段と、入力手段により入力された波形信号
の周期を検出する周期検出手段と、周期検出手段で検出
した周期に基づいて、波形信号を分割する分割手段と、
分割された個々の波形信号をウェーブレット変換するた
めのウェーブレット変換手段と、分割された個々の波形
信号のウェーブレット変換結果を時系列に連続表示する
表示手段と、を備えたことを特徴とする波形信号解析装
置である。
【0016】本発明によれば、周期検出手段によって波
形信号の周期が検出され、その周期によって波形信号が
分割されるため、ウェーブレット変換結果における特徴
ピークの時系列な挙動を調べることが容易である。ま
た、本発明によれば、分割された個々の波形信号のウェ
ーブレット変換結果が時系列に連続表示されるため、ウ
ェーブレット変換結果についての比較、検討を容易に行
うことができる。
【0017】また、本発明は、波形信号を入力する入力
手段と、入力手段により入力された波形信号をウェーブ
レット変換するウェーブレット変換手段と、ウェーブレ
ット変換手段による変換結果の周期を検出する周期検出
手段と、周期検出手段で検出した周期に基づいて、変換
結果を分割する分割手段と、分割された変換結果を時系
列に連続表示する表示手段と、を備えたことを特徴とす
る波形信号解析装置である。
【0018】本発明によれば、周期検出手段によってウ
ェーブレット変換結果の周期が検出され、その周期によ
ってウェーブレット変換結果が分割されるため、ウェー
ブレット変換結果における特徴ピークの時系列な挙動を
調べることが容易である。また、本発明によれば、分割
されたウェーブレット変換結果が時系列に連続表示され
るため、分割された各ウェーブレット変換結果について
の比較、検討を容易に行うことができる。
【0019】
【発明の実施の形態】以下、図面を参照して本発明の実
施の形態について説明する。
【0020】図1は、本発明の第1の実施の形態による
波形信号解析装置を示す構成概略図である。図1に示す
ように、本発明の第1の実施の形態の波形信号解析装置
1は、波形信号を入力する入力手段10と、入力手段1
0により入力された波形信号の周期性を検出する周期検
出手段20と、周期検出手段20で検出した周期によっ
て波形信号を分割する分割手段30とを備えている。
【0021】本実施の形態の周期検出手段20は、入力
手段10により入力された波形信号にフーリエ変換を施
してフーリエスペクトルを得るフーリエ変換部21と、
フーリエ変換部21により得られたフーリエスペクトル
から、その強度の最も大きいピーク間隔を波形信号の周
期として識別する周期識別部22とを有している。
【0022】また、本実施の形態の分割手段30は、周
期検出手段20で検出した周期を整数倍して分割周期と
する周期設定部31を有している。この「整数倍」のた
めの「整数」は、操作者によって適宜に選択されるよう
になっている。
【0023】分割手段30には、分割された個々の波形
信号をウェーブレット変換するためのウェーブレット変
換手段40が接続され、ウェーブレット変換手段40に
は、分割された個々の波形信号のウェーブレット変換結
果を時系列に連続表示する表示手段50が接続されてい
る。
【0024】本実施の形態のウェーブレット変換手段4
0は、マザーウェーブレットの種類を指定するマザーウ
ェーブレット選択部41と、ウェーブレット変換の対象
とする時間領域を指定するウェーブレット変換範囲指定
部42と、ウェーブレット変換のスケールパラメータ及
びシフトパラメータを設定する変換パラメータ設定部4
3と、ウェーブレット係数を計算するウェーブレット演
算部44と、ウェーブレット変換結果を出力する出力部
45とを有している。
【0025】本実施の形態の表示手段50は、ウェーブ
レット変換結果の強度について、設定数値範囲毎に異な
る色で分けて表示するようになっている。設定数値範囲
は、操作者によって適宜に複数が設定され、結果として
等高線図状の表示が得られるようになっている。
【0026】また、本実施の形態の表示手段50は、複
数のウェーブレット変換結果を、指定した周波数領域毎
に統合して表示可能となっている。
【0027】また、本実施の形態の表示手段50は、表
示画面51と、表示画面中の任意の位置を指定するため
のマーカー操作部52とを有している。そして、マーカ
ー操作部52によって表示画面51中に表示されたウェ
ーブレット変換結果が指定された場合、その指定された
ウェーブレット変換結果についての時刻、周波数及び強
度が表示画面51中に表示されるようになっている。
【0028】また、本実施の形態の表示手段50は、時
系列の異なるウェーブレット変換結果の差分、変化率及
び全変換結果の平均を計算する演算部53を有してい
る。そして、演算部53による計算結果は、表示画面5
1中に表示されるようになっている。
【0029】さらに本実施の形態の表示手段50は、ウ
ェーブレット変換結果の周波数軸での断面である「指定
する周波数におけるピーク強度の時間変化」もしくは、
時間軸での断面である「指定する時刻での各周波数成分
のピーク強度分布」を表示可能となっている。
【0030】その他、本実施の形態の波形信号解析装置
1は、表示手段50に、診断手段60が接続されてい
る。診断手段60は。診断対象種毎に定められている診
断アルゴリズムが記憶された記憶部62と、記憶部62
に記憶された診断アルゴリズムを用いて表示手段50の
出力状態を判定する状態判定部61とを有している。状
態判定部61は、具体的には、例えば解析対象信号の発
信源の異常度等を判定することができる。また、判定結
果についても、記憶部62に記憶可能である。
【0031】次に、このような構成からなる本実施の形
態の波形信号解析装置の作用について、図2を用いて説
明する。図2は、本実施の形態の波形信号解析装置1の
作用の一例を説明する図である。
【0032】まず、入力手段10において、解析対象と
なる波形信号f(t)が入力される。
【0033】次に周期検出手段20において、フーリエ
変換部21が入力波形信号をフーリエ変換し、周期識別
部22がフーリエ変換部21により得られたフーリエス
ペクトルから、その強度の最も大きいピーク間隔を周期
として識別する。図2に示す例では、波形信号f(t)
が電源波形などのパイロット信号と同期する特徴を有し
ており、パイロット信号の周期t1 が波形信号f(t)
の周期として識別される。
【0034】次に、分割手段30の周期設定部31が、
周期識別部22によって設定された周期の整数倍で分割
周期を定義する。この「整数倍」のための「整数」は、
操作者によって適宜に選択される。分割周期が定義され
た後、その分割周期によって解析対象波形が分割され
る。図2に示す例では、識別された周期t1 の等倍(1
倍)に分割周期が定義され、波形信号f(t)が周期t
1 毎に分割される。
【0035】分割された個々の波形信号は、ウェーブレ
ット変換手段40によって、それぞれウェーブレット変
換が施される。ウェーブレット変換式を以下に示す。
【0036】
【数1】 ウェーブレット変換手段40においては、まずマザーウ
ェーブレット選択部41が、入力波形の特徴を検出し易
いマザーウェーブレットψ* ((t−b)/a)を選択
する。変換範囲指定部42は、分割された各々の入力波
形を指定する。変換パラメータ設定部43は、ウェーブ
レット変換のスケールパラメータ及びシフトパラメータ
を適宜に設定する。そして、ウェーブレット演算部44
が、ウェーブレット係数を計算し、出力部45からこれ
を表示手段50に出力する。これにより、それぞれの分
割領域毎の時間−周波数解析が実現される。
【0037】出力部45から出力される各ウェーブレッ
ト変換結果、すなわち、それぞれの分割領域毎の時間−
周波数解析結果は、図2に示すように、表示手段50に
よって時系列に連続に動画表示される。
【0038】本実施の形態の表示手段50は、図2に示
すように、3次元情報である時間−周波数解析結果を、
時間−周波数2次元平面において、強度について設定数
値範囲毎に異なる色で分けて等高線図状に表示する。
【0039】ここで、本実施の形態による表示手段50
の各種の機能について、図3乃至図6を用いて説明す
る。
【0040】本実施の形態の表示手段50は、図3に示
すように、複数のウェーブレット変換結果を、指定した
周波数領域毎に統合して表示することが可能である(複
数変換結果統合表示機能)。この複数変換結果統合表示
機能は、同時に測定した複数の入力波形信号をすべて等
しい周期で分割して、それぞれのウェーブレット変換結
果を統合して表示する場合に有効である。
【0041】図3を用いて、複数変換結果統合表示機能
を具体的に説明する。図3(a)に示すように、例えば
ある対象物にセンサ1〜4が取り付けられ、その対象物
が発生する波形信号が各センサ1〜4で同時に測定され
ると、図3(b)に示すように、それぞれのセンサ信号
のウェーブレット変換結果が各々別個に演算される。
【0042】表示手段50は、例えばセンサ1〜4の時
刻0〜t及び時刻t〜2tにおけるウェーブレット変換
結果を、200Hz〜400Hzの周波数領域と400
Hz〜800Hzの周波数領域とに分け、図3(c)に
示すように、それぞれを統合して表示する。この場合、
200Hz〜400Hzの周波数領域では特徴的な等高
線表示が得られ、400Hz〜800Hzの周波数領域
では何のピークも見られない。また、時刻0〜tの結果
と時刻t〜2tの結果は連続動画表示されるが、200
Hz〜400Hzの周波数領域についての等高線図表示
はその表示位置が変化する一方、400Hz〜800H
zの周波数領域について何の変化も現れない。
【0043】このように、複数信号のウェーブレット変
換結果を統合表示することにより、複数信号の表示情報
をあたかも一つの表示情報として認識し、対象毎の状態
比較や特徴抽出をより容易に行うことが可能となる。ま
た、対象物について空間的な情報をとらえることができ
たり、不必要な情報の効率的な省略化なども可能とな
る。さらには、複数の入力波形信号の時間−周波数分布
の相対的変動を調べることもでき、入力波形間の相関も
分かる。
【0044】また本実施の形態の表示手段50は、図4
に示すように、マーカー操作部52によって表示画面5
1中に表示されたウェーブレット変換結果が指定された
場合、そのウェーブレット変換結果の時刻、周波数及び
強度を表示画面51中に表示することができる(マーカ
ー機能)。
【0045】このマーカー機能は、表示画面51中のウ
ェーブレット変換結果やその他の数値情報を速やかに表
示させる機能である。図4に示す例では、数値情報を知
りたい表示画面51中の点をマーカー操作部52の操作
により矢印形のポインタによってマークすると、矢印形
のポインタが示す位置に表示されている情報についての
詳細(座標や強度など)が速やかに表示画面51中に表
示される。
【0046】このマーカー機能により、ウェーブレット
変換で求めた時間−周波数分布の各ピークの数値情報、
例えばピーク間の時間差(位相差)や強度差などを容易
に得ることができる。
【0047】さらに表示手段50は、図5に示すよう
に、時系列のウェーブレット変換結果の差分、変化率、
さらには全変換結果の平均等を演算部53により計算し
て、表示画面51中に表示させることができる(演算機
能)。
【0048】この演算機能は、ウェーブレット変換結果
の異なる2画像間の差分や変化率、各ウェーブレット変
換結果の絶対値、全変換結果の平均などを計算して表示
する機能である。図5に示す例では、2つの時間−周波
数分布画像のピーク強度の差をとって差分を表示した
り、2つの画像の差分を時系列的に早い画像の強度で割
ることで変化率を演算、表示している。
【0049】この演算機能により、ウェーブレット変換
結果の時間変化分等の特徴等を抽出して表示することが
できるので、時間−周波数分布の傾向解析がより容易に
なる。
【0050】また表示手段50は、図6に示すように、
ウェーブレット変換結果の周波数軸での断面である「指
定する周波数におけるピーク強度の時間変化」もしく
は、時間軸での断面である「指定する時刻での各周波数
成分のピーク強度分布」を表示することも可能である
(断面表示機能)。
【0051】この断面表示機能は、ウェーブレット変換
結果の周波数軸での断面や時間軸での断面を表示する機
能である。図6に示す例では、「指定周波数νHz成分
強度の時間変化」や「指定時刻t1 sec(位相θra
d)における周波数成分のピーク強度分布」を表示して
いる。
【0052】この断面表示機能も、時間−周波数分布の
傾向解析にとって効果的である。
【0053】表示手段50に表示された解析結果に基づ
いて、診断手段60が波形信号の発信源の異常等を判定
する。具体的には、記憶部62に記憶された診断対象種
毎に定められている診断アルゴリズムを用いて、状態判
定部61が表示手段50に出力された解析結果を判定す
る。判定結果は、記憶部62に記憶される。
【0054】以上のように本実施の形態によれば、周期
検出手段20によって波形信号の周期が検出され、その
周期によって波形信号f(t)が分割されるため、ウェ
ーブレット変換結果における特徴ピークの時系列な挙動
を調べることが容易である。
【0055】また、本実施の形態によれば、分割された
個々の波形信号のウェーブレット変換結果が表示手段5
0によって時系列に連続表示されるため、ウェーブレッ
ト変換結果についての比較、検討を容易に行うことがで
きる。
【0056】また本実施の形態によれば、周期設定部3
1によって、一度に表示する時間領域を周期性を保ちな
がら拡大できるので、解析時間周波数分布の時系列変化
を大まかに解析することも可能である。
【0057】また本実施の形態によれば、マザーウェー
ブレット変換手段40によって、解析対象波形の時間−
周波数解析にもっとも有効なマザーウェーブレット関数
を選択することができるので、時間−周波数分布内の特
徴成分の抽出感度を上げることができる。また、変換時
間範囲とスケールパラメータ及びシフトパラメータを自
由に変えられるので、ウェーブレット係数の演算におけ
る計算量の効率化を図ることができる。
【0058】また本実施の形態の表示手段50には、前
述のような種々の便利な機能が付加されており、ウェー
ブレット変換結果のピークの特徴抽出や傾向解析を容易
に行うことができ、解析波形に含まれる異常などを判断
しやすい表示が実現できる。
【0059】また本実施の形態の波形信号解析装置1
は、診断装置60により、波形信号の状態の正常あるい
は異常の判定、異常モードの特定、異常の度合い判定な
どが可能である。また、必要なデータベースを記憶部6
2内に取り揃えることにより、波形信号による異常診断
対象を拡張することができる。
【0060】なお、周期検出手段20は、フーリエ変換
部21を有する代わりに、入力波形信号のピーク値を検
出する等の機能を有していてもよい。この場合、ピーク
値の間隔から、ウェーブレット変換結果を動画表示する
のに最適な分割周期を設定することができる。
【0061】次に、本発明の第2の実施の形態の波形信
号解析装置について図7を用いて説明する。図7は、第
2の実施の形態の波形信号解析装置の概略構成図であ
る。
【0062】図7に示すように、本実施の形態の波形信
号解析装置1’は、ウェーブレット変換手段40’が直
接入力手段10に接続され、入力手段10により入力さ
れた波形信号をそのままウェーブレット変換するように
なっており、周期検出手段20’は、入力波形信号では
なくウェーブレット変換手段40’による変換結果の周
期性を検出するようになっており、分割手段30’は、
周期検出手段20’で検出した周期によって、入力波形
信号ではなく変換結果を分割するようになっている。
【0063】その他の構成は、図1に示す第1の実施の
形態と同様の構成である。第2の実施の形態において、
図1に示す第1の実施の形態と同一の部分には同一の符
号を付して詳細な説明は省略する。
【0064】周期検出手段20’は、例えば、ウェーブ
レット変換結果にフーリエ変換を施してフーリエスペク
トルを得るフーリエ変換部21’と、フーリエ変換部2
1’により求められたフーリエスペクトルから、その強
度の最も大きいピーク間隔を分割周期として識別する周
期識別部22’とを有する。具体的には、図8に示すよ
うに、フーリエ変換部21’がウェーブレット変換結果
に基づいて特徴的な指定周波数νHz(例えば、ピーク
が多く出現していると思われる周波数)を指定し、その
時間変化に対してフーリエ変換(FFT)を施し、その
結果から周期識別部22’が分割周期を識別する。
【0065】あるいは周期検出手段20’は、ウェーブ
レット変換結果において、図9に示すように、指定した
時間−周波数領域、例えばピーク強度が指定したしきい
値を超える時間周波数領域(危険領域)内の最大値tb
や最小値、重心tc 、指定領域開始時刻ta または指定
領域終了時刻td を検出し、これから分割周期を求めて
も良い。この場合、分割周期は不定周期となる場合があ
る。
【0066】本実施の形態によれば、周期検出手段2
0’によってウェーブレット変換結果の周期が検出さ
れ、その周期によってウェーブレット変換結果が分割さ
れるため、ウェーブレット変換結果における特徴ピーク
の時系列な挙動を調べることが容易である。
【0067】また、本実施の形態によれば、分割された
ウェーブレット変換結果が時系列に連続表示されるた
め、ウェーブレット変換結果についての検討を容易に行
うことができる。
【0068】
【発明の効果】本発明によれば、周期検出手段によって
波形信号の周期性が検出され、その周期によって波形信
号が分割されるため、ウェーブレット変換結果における
特徴ピークの時系列な挙動を調べることが容易である。
また、本発明によれば、分割された個々の波形信号のウ
ェーブレット変換結果が時系列に連続表示されるため、
ウェーブレット変換結果についての検討を容易に行うこ
とができる。
【0069】あるいは本発明によれば、周期検出手段に
よってウェーブレット変換結果の周期性が検出され、そ
の周期によってウェーブレット変換結果が分割されるた
め、ウェーブレット変換結果における特徴ピークの時系
列な挙動を調べることが容易である。また、本発明によ
れば、分割されたウェーブレット変換結果が時系列に連
続表示されるため、ウェーブレット変換結果についての
検討を容易に行うことができる。
【図面の簡単な説明】
【図1】本発明による波形信号解析装置の第1の実施の
形態を示す構成概略図。
【図2】図1の波形信号解析装置の作用を説明する図。
【図3】図1の表示手段の複数変換結果統合表示機能を
説明する図。
【図4】図1の表示手段のマーカー機能を説明する図。
【図5】図1の表示手段の演算機能を説明する図。
【図6】図1の表示手段の断面表示機能を説明する図。
【図7】本発明による波形信号解析装置の第2の実施の
形態を示す構成概略図。
【図8】図7の周期検出手段の一例の作用を説明する
図。
【図9】図7の周期検出手段の他の例の作用を説明する
図。
【符号の説明】
1、1’ 波形信号解析装置 10 入力手段 20、20’ 周期検出手段 21、21’ フーリエ変換部 22、22’ 周期識別部 30、30’ 分割手段 31、31’ 周期設定部 40、40’ ウェーブレット変換手段 41 マザーウェーブレット選択部 42 変換範囲指定部 43 変換パラメータ設定部 44 ウェーブレット演算部 45 出力部 50 表示手段 51 表示画面 52 マーカー操作部 53 演算部 60 診断手段 61 状態判定部 62 記憶部

Claims (12)

    【特許請求の範囲】
  1. 【請求項1】波形信号を入力する入力手段と、 入力手段により入力された波形信号の周期を検出する周
    期検出手段と、 周期検出手段で検出した周期に基づいて、波形信号を分
    割する分割手段と、 分割された個々の波形信号をウェーブレット変換するた
    めのウェーブレット変換手段と、 分割された個々の波形信号のウェーブレット変換結果を
    時系列に連続表示する表示手段と、を備えたことを特徴
    とする波形信号解析装置。
  2. 【請求項2】ウェーブレット変換手段は、 マザーウェーブレットの種類を指定するマザーウェーブ
    レット選択部と、 ウェーブレット変換の対象とする時間領域を指定する変
    換範囲指定部と、 ウェーブレット変換のスケールパラメータ及びシフトパ
    ラメータを設定する変換パラメータ設定部と、 ウェーブレット係数を計算するウェーブレット演算部
    と、 ウェーブレット変換結果を出力する出力部と、を有する
    ことを特徴とする請求項1に記載の波形信号解析装置。
  3. 【請求項3】周期検出手段は、 入力手段により入力された波形信号にフーリエ変換を施
    してフーリエスペクトルを得るフーリエ変換部と、 フーリエ変換部により得られたフーリエスペクトルから
    波形信号の周期を識別する周期識別部と、を有すること
    を特徴とする請求項1または2に記載の波形信号解析装
    置。
  4. 【請求項4】分割手段は、周期検出手段で検出した周期
    の整数倍に分割周期を設定する周期設定部を有すること
    を特徴とする請求項1乃至3のいずれかに記載の波形信
    号解析装置。
  5. 【請求項5】表示手段は、ウェーブレット変換結果の強
    度について、設定数値範囲毎に異なる色で分けて表示す
    るようになっていることを特徴とする請求項1乃至4の
    いずれかに記載の波形信号解析装置。
  6. 【請求項6】表示手段は、複数のウェーブレット変換結
    果を指定した周波数領域毎に統合して表示可能となって
    いることを特徴とする請求項1乃至5のいずれかに記載
    の波形信号解析装置。
  7. 【請求項7】表示手段は、表示画面と、表示画面中の任
    意の位置を指定するためのマーカー操作部とを有し、 マーカー操作部によって表示画面中のウェーブレット変
    換結果が指定された場合に、そのウェーブレット変換結
    果の時刻、周波数及び強度を表示画面中に表示するよう
    になっていることを特徴とする請求項1乃至6のいずれ
    かに記載の波形信号解析装置。
  8. 【請求項8】表示手段は、異なるウェーブレット変換結
    果の差分、変化率及び全変換結果の平均を計算する演算
    部を有していることを特徴とする請求項1乃至7のいず
    れかに記載の波形信号解析装置。
  9. 【請求項9】表示手段は、 ウェーブレット変換結果の周波数軸での断面である、指
    定する周波数におけるピーク強度の時間変化、もしく
    は、 時間軸での断面である、指定する時刻での各周波数成分
    のピーク強度分布を表示可能となっていることを特徴と
    する請求項1乃至8のいずれかに記載の波形信号解析装
    置。
  10. 【請求項10】波形信号を入力する入力手段と、 入力手段により入力された波形信号をウェーブレット変
    換するウェーブレット変換手段と、 ウェーブレット変換手段による変換結果の周期を検出す
    る周期検出手段と、 周期検出手段で検出した周期に基づいて、変換結果を分
    割する分割手段と、 分割された変換結果を時系列に連続表示する表示手段
    と、を備えたことを特徴とする波形信号解析装置。
  11. 【請求項11】周期検出手段は、 ウェーブレット変換結果にフーリエ変換を施してフーリ
    エスペクトルを得るフーリエ変換部と、 フーリエ変換部により得られたフーリエスペクトルから
    ウェーブレット変換結果の周期を識別する周期識別部
    と、を有することを特徴とする請求項10に記載の波形
    信号解析装置。
  12. 【請求項12】周期検出手段は、 ウェーブレット変換結果において、指定した時間周波数
    領域内の最大値や最小値、重心、指定領域開始時刻また
    は指定領域終了時刻を検出するようになっていることを
    特徴とする請求項10に記載の波形信号解析装置。
JP08478199A 1999-03-26 1999-03-26 波形信号解析装置 Expired - Lifetime JP3730435B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP08478199A JP3730435B2 (ja) 1999-03-26 1999-03-26 波形信号解析装置
US09/533,937 US6542836B1 (en) 1999-03-26 2000-03-23 Waveform signal analyzer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP08478199A JP3730435B2 (ja) 1999-03-26 1999-03-26 波形信号解析装置

Publications (3)

Publication Number Publication Date
JP2000275098A true JP2000275098A (ja) 2000-10-06
JP2000275098A5 JP2000275098A5 (ja) 2005-05-26
JP3730435B2 JP3730435B2 (ja) 2006-01-05

Family

ID=13840247

Family Applications (1)

Application Number Title Priority Date Filing Date
JP08478199A Expired - Lifetime JP3730435B2 (ja) 1999-03-26 1999-03-26 波形信号解析装置

Country Status (2)

Country Link
US (1) US6542836B1 (ja)
JP (1) JP3730435B2 (ja)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001099757A (ja) * 1999-09-29 2001-04-13 Takasago Thermal Eng Co Ltd ベアリングの異常検出方法
JP2003207389A (ja) * 2002-01-11 2003-07-25 Denso Corp 振動波判定装置
JP2004077279A (ja) * 2002-08-19 2004-03-11 Denso Corp 振動波判定装置
JP2005331300A (ja) * 2004-05-18 2005-12-02 Tektronix Japan Ltd 周波数分析装置の表示方法
JP2008076300A (ja) * 2006-09-22 2008-04-03 Shimadzu Corp 測定装置および材料試験機
JP2008199380A (ja) * 2007-02-14 2008-08-28 Mitsubishi Electric Corp 波形信号解析装置
JP2008250978A (ja) * 2007-03-05 2008-10-16 Nagaoka Univ Of Technology 触覚の視覚的表現方法および触覚の視覚的表現装置
JP2010266327A (ja) * 2009-05-14 2010-11-25 Denso Corp 設備診断装置及び設備診断方法
US7933742B2 (en) 2006-02-01 2011-04-26 Jtekt Corporation Abnormality diagnosing method for sound or vibration and abnormality diagnosing apparatus for sound or vibration
JP2011237414A (ja) * 2010-05-06 2011-11-24 Tektronix Inc 試験測定機器
JP5937727B1 (ja) * 2015-05-22 2016-06-22 Dmg森精機株式会社 工具摩耗評価装置
JPWO2016157912A1 (ja) * 2015-03-31 2018-01-25 株式会社東芝 電力機器の絶縁診断装置および絶縁診断方法
JP2020169916A (ja) * 2019-04-04 2020-10-15 株式会社日立製作所 診断装置
WO2020245970A1 (ja) * 2019-06-06 2020-12-10 三菱電機ビルテクノサービス株式会社 分析装置
CN112859011A (zh) * 2021-01-12 2021-05-28 武汉大学 一种单波长机载测深雷达波形信号提取方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7130749B2 (en) * 2003-07-11 2006-10-31 Texas Instruments Incorporated Wavelet analysis of signals to determine characteristics of anomalies in a wire
US7680208B2 (en) * 2004-02-25 2010-03-16 Nokia Corporation Multiscale wireless communication
GR1005713B (el) * 2006-05-11 2007-11-07 Ιωαννης Τσουμας Ταχεια ανιχνευση κομμενων μπαρων σε ηλεκτρικες ασυγχρονες μηχανες μεσω αναλυσης του ρευματος του στατη με τη χρηση κυματοειδων συναρτησεων μικρης διαρκειας
JP4763522B2 (ja) * 2006-06-14 2011-08-31 株式会社ブリヂストン タイヤ検査装置
US8412468B1 (en) 2009-04-08 2013-04-02 Western Kentucky University Method and apparatus for wavelet based elemental spectrum analysis
CN106443760B (zh) * 2016-09-14 2019-02-12 湖北锐世数字医学影像科技有限公司 一种闪烁脉冲简化模型、重建及能量获取的方法
US11475341B2 (en) * 2017-12-28 2022-10-18 Tata Consultancy Services Limited Systems and methods for obtaining optimal mother wavelets for facilitating machine learning tasks

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6462797A (en) * 1987-09-02 1989-03-09 Sharp Kk Periodical variance detector
JPH0883265A (ja) * 1994-09-13 1996-03-26 Toshiba Corp 振動信号解析装置
JPH10258974A (ja) * 1996-09-13 1998-09-29 Toshiba Corp 非定常信号解析装置及び非定常信号解析プログラムを記録した媒体
JPH10274558A (ja) * 1997-01-28 1998-10-13 Matsushita Electric Works Ltd 回転機器の異常診断方法およびその装置
JPH10281859A (ja) * 1997-04-07 1998-10-23 Mitsui Eng & Shipbuild Co Ltd 異常診断方法および装置
JPH10288639A (ja) * 1997-02-13 1998-10-27 Toshiba Corp 電気機器の異常監視装置および方法
JPH1153342A (ja) * 1997-07-31 1999-02-26 Yamatake Honeywell Co Ltd データ変換方法およびデータ変換装置およびプログラム記憶媒体

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR940002854B1 (ko) * 1991-11-06 1994-04-04 한국전기통신공사 음성 합성시스팀의 음성단편 코딩 및 그의 피치조절 방법과 그의 유성음 합성장치
JP2736584B2 (ja) * 1992-06-30 1998-04-02 花王株式会社 高速演算装置および方法,信号解析システム,ならびに計測装置および方法
JPH0895955A (ja) 1994-09-14 1996-04-12 Nippon Butsusei Kk ウェーブレット変換波形信号解析システムとウェーブレット変換波形信号解析ツールによる解析方法及びそれに用いるウェーブレット変換波形信号解析装置
US5737716A (en) * 1995-12-26 1998-04-07 Motorola Method and apparatus for encoding speech using neural network technology for speech classification
JP2735064B2 (ja) 1996-01-31 1998-04-02 日本電気株式会社 波形解析装置
US5841890A (en) * 1996-05-06 1998-11-24 Northrop Grumman Corporation Multi-dimensional wavelet tomography
JPH1065527A (ja) * 1996-08-21 1998-03-06 Nec Corp 位相同期発振回路
US6108609A (en) * 1996-09-12 2000-08-22 National Instruments Corporation Graphical system and method for designing a mother wavelet
US5870691A (en) * 1996-12-06 1999-02-09 Amoco Corporation Spectral decomposition for seismic interpretation
US6208951B1 (en) * 1998-05-15 2001-03-27 Council Of Scientific & Industrial Research Method and an apparatus for the identification and/or separation of complex composite signals into its deterministic and noisy components

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6462797A (en) * 1987-09-02 1989-03-09 Sharp Kk Periodical variance detector
JPH0883265A (ja) * 1994-09-13 1996-03-26 Toshiba Corp 振動信号解析装置
JPH10258974A (ja) * 1996-09-13 1998-09-29 Toshiba Corp 非定常信号解析装置及び非定常信号解析プログラムを記録した媒体
JPH10274558A (ja) * 1997-01-28 1998-10-13 Matsushita Electric Works Ltd 回転機器の異常診断方法およびその装置
JPH10288639A (ja) * 1997-02-13 1998-10-27 Toshiba Corp 電気機器の異常監視装置および方法
JPH10281859A (ja) * 1997-04-07 1998-10-23 Mitsui Eng & Shipbuild Co Ltd 異常診断方法および装置
JPH1153342A (ja) * 1997-07-31 1999-02-26 Yamatake Honeywell Co Ltd データ変換方法およびデータ変換装置およびプログラム記憶媒体

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001099757A (ja) * 1999-09-29 2001-04-13 Takasago Thermal Eng Co Ltd ベアリングの異常検出方法
JP2003207389A (ja) * 2002-01-11 2003-07-25 Denso Corp 振動波判定装置
JP2004077279A (ja) * 2002-08-19 2004-03-11 Denso Corp 振動波判定装置
JP2005331300A (ja) * 2004-05-18 2005-12-02 Tektronix Japan Ltd 周波数分析装置の表示方法
US7933742B2 (en) 2006-02-01 2011-04-26 Jtekt Corporation Abnormality diagnosing method for sound or vibration and abnormality diagnosing apparatus for sound or vibration
JP2008076300A (ja) * 2006-09-22 2008-04-03 Shimadzu Corp 測定装置および材料試験機
JP2008199380A (ja) * 2007-02-14 2008-08-28 Mitsubishi Electric Corp 波形信号解析装置
JP4684250B2 (ja) * 2007-02-14 2011-05-18 三菱電機株式会社 波形信号解析装置
JP2008250978A (ja) * 2007-03-05 2008-10-16 Nagaoka Univ Of Technology 触覚の視覚的表現方法および触覚の視覚的表現装置
JP2010266327A (ja) * 2009-05-14 2010-11-25 Denso Corp 設備診断装置及び設備診断方法
JP2011237414A (ja) * 2010-05-06 2011-11-24 Tektronix Inc 試験測定機器
JPWO2016157912A1 (ja) * 2015-03-31 2018-01-25 株式会社東芝 電力機器の絶縁診断装置および絶縁診断方法
JP5937727B1 (ja) * 2015-05-22 2016-06-22 Dmg森精機株式会社 工具摩耗評価装置
WO2016189911A1 (ja) * 2015-05-22 2016-12-01 Dmg森精機株式会社 工具摩耗評価装置
JP2016215333A (ja) * 2015-05-22 2016-12-22 Dmg森精機株式会社 工具摩耗評価装置
JP2020169916A (ja) * 2019-04-04 2020-10-15 株式会社日立製作所 診断装置
JP7352371B2 (ja) 2019-04-04 2023-09-28 株式会社日立製作所 診断装置
WO2020245970A1 (ja) * 2019-06-06 2020-12-10 三菱電機ビルテクノサービス株式会社 分析装置
JPWO2020245970A1 (ja) * 2019-06-06 2021-10-14 三菱電機ビルテクノサービス株式会社 分析装置
JP7001200B2 (ja) 2019-06-06 2022-01-19 三菱電機ビルテクノサービス株式会社 分析装置
CN112859011A (zh) * 2021-01-12 2021-05-28 武汉大学 一种单波长机载测深雷达波形信号提取方法
CN112859011B (zh) * 2021-01-12 2022-06-07 武汉大学 一种单波长机载测深雷达波形信号提取方法

Also Published As

Publication number Publication date
JP3730435B2 (ja) 2006-01-05
US6542836B1 (en) 2003-04-01

Similar Documents

Publication Publication Date Title
JP3730435B2 (ja) 波形信号解析装置
JP5783808B2 (ja) 異常音診断装置
US20070272599A1 (en) Moving body inspection apparatus and method of comparing phases between movement waveforms
US6269263B1 (en) Method for estimating heart rate variability and apparatus for embodying estimation
KR20160016460A (ko) 관심영역의 시각화 장치 및 방법
JP6231681B2 (ja) 呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体
US20020169365A1 (en) Apparatus for measuring activity of autonomic nerve
KR102539142B1 (ko) 스펙트럼 분석 장치와 방법, 및 혈당 측정 장치
KR20160000810A (ko) 파형종류 도수분포를 이용한 가속도맥파 분석 장치 및 분석 방법
JP2008061824A (ja) 医療用測定器、生体信号波形抽出方法、および生体信号波形抽出プログラム
KR20230003205A (ko) 광 파형 계측 장치 및 계측 방법
KR101543146B1 (ko) 진동 장치의 상태 판단 방법
US11109810B2 (en) Information display device, biological signal measurement system, and computer program product
US20220079498A1 (en) Data processing apparatus, data processing method, and storage medium storing program
JP2003065956A (ja) 蛍光ピーク検出方法及び分光蛍光光度計
US11762360B2 (en) Sensor device, data acquisition method, and current monitoring system
JP6403311B2 (ja) 心拍状態解析装置
JP3561151B2 (ja) 異常診断装置および異常診断方法
JP2008220556A (ja) 心拍測定システム及び方法
JP2009300192A (ja) 亀裂検出装置及び亀裂検出方法
JP2022041186A (ja) ラマン散乱光による温度測定方法及びラマン分光分析装置
JP6793077B2 (ja) 超音波診断装置及び波形処理方法
JPH08329046A (ja) ウェーブレット変換を用いた信号解析装置
JP7001806B2 (ja) 呼吸音解析装置、呼吸音解析方法、コンピュータプログラム及び記録媒体
JP3952176B2 (ja) 生体波形解析装置及び生体波形解析方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040730

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040730

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050527

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050628

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050829

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: 20050930

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20051006

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081014

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091014

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101014

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111014

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111014

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121014

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131014

Year of fee payment: 8

EXPY Cancellation because of completion of term