JPH1019657A - 励振力推定装置 - Google Patents

Abstract

(57)【要約】 【課題】ノイズに対する安定性が強く、振幅、加振力持
続時間とも良好な推定結果が得られ、連続的な励振力の
推定にも適用できるようにする。 【解決手段】加振試験時にインパルスハンマ１２を用い
て被測定機器１１の励振力着力点Ａに既知の荷重を加
え、振動応答測定点Ｂにおける加速度応答を加速度検出
器１３により検出して、インパルスハンマ１２の励振力
及び加速度応答を第１のメモリ装置１４に記憶する。第
１の演算装置１５は、第１のメモリ装置１４に記憶した
加振試験時の励振力及び加速度応答に基づいて機器１１
の振動応答測定点Ｂと励振力着力点Ａとの間のインパル
ス応答関数を計算する。タイムウィンドウ処理ルーチン
１８は、振動応答を時間領域で短く切り、第２の演算装
置１７はそれぞれをパルス状の励振力に対する過渡応答
とみなして励振力を推定した後、重畳処理部１９で重ね
合わせて連続波形とする。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、例えばエンジン
等、稼働時に衝撃的励振力が作用する機器における励振
力を推定する励振力推定装置に関する。

【０００２】

【従来の技術】従来、例えばエンジン等、稼働時に衝撃
的励振力が作用する機器における励振力を推定する励振
力推定装置は、図１０に示すように構成されている。図
１０において、ａはインパルス応答関数計測部、ｂは機
器稼働時の振動応答計測部、ｃは励振力推定演算部で、
インパルス応答関数計測部ａで求めたインパルス応答関
数及び機器稼働時の振動応答計測部ｂで計測した機器稼
働時の振動応答を用いて、励振力推定演算部ｃで励振力
推定値を計算する。

【０００３】上記インパルス応答関数計測部ａにおい
て、１は構造系、つまり、被計測機器で、加振試験を行
なう際、励振力着力点Ａにインパルスハンマ２により既
知の荷重（インパルス）が加えられる。また、被測定機
器１には、振動応答測定点Ｂに加速度検出器３を設け、
上記インパルスハンマ２により荷重が加えられた際の加
速度応答（加速度信号）を検出している。上記インパル
スハンマ２の励振力及び加速度検出器３により検出した
加速度応答は、第１のメモリ装置４に記憶され、第１の
演算装置５に送られる。この第１の演算装置５は、上記
入力信号に基づいてインパルス応答関数計算を行ない、
この計算によって得られたインパルス応答関数を振動応
答計測部ｂの第２のメモリ装置６に出力する。

【０００４】また、この第２のメモリ装置６には、被測
定機器１の稼働時における加速度応答が加速度検出器３
により検出されて記憶される。そして、第２のメモリ装
置６に記憶された機器稼働時の加速度応答及びインパル
ス応答関数は、励振力推定演算部ｃの第２の演算装置７
に入力される。

【０００５】上記の構成において、励振力推定値は次の
ようにして求められる。 （１）まず、インパルスハンマ２により被計測機器１の
励振力着力点Ａを叩き、インパルスハンマ２の励振力と
加速度検出器３により検出した加速度応答とを第１のメ
モリ装置４に一時保存し、これを第１の演算装置５にて
インパルス応答関数を計算する。

【０００６】（２）次に機器稼働時の振動応答計測部ｂ
は、被計測機器１の稼働時の加速度応答を加速度検出器
３で検出し、その計測データを第２のメモリ装置６に保
存する。

【０００７】（３）励振力推定演算部ｃは、第１の演算
装置５にて計算されたインパルス応答関数と振動応答計
測部ｂで計測した加速度応答の計測データとを第２の演
算装置７に入力し、周波数領域においての励振力推定計
算を行なう。

【０００８】構造系の振動応答は、周波数領域では、励
振力と系の周波数応答関数の積で表されるため、機器稼
働時の振動応答スペクトルＸ（ω）（振動応答の周波数
領域における表現）を周波数応答関数Ｈ（ω）で除する
ことにより、構造系に負荷された励振力のスペクトルＦ
（ω）を求める。周波数領域で求められた励振力Ｆ
（ω）を時間領域に変換することにより、励振力の時系
列波形ｆ（ｔ）を得る。

【０００９】上記励振力推定演算部ｃは、図１１に示す
ように離散フーリエ変換器７ａ及び離散フーリエ逆変換
器７ｂにより構成されている。上記離散フーリエ変換器
７ａは、時系列データを周波数領域のデータに変換する
演算器であり、機器稼働時の振動応答及び予めインパル
ス応答関数計測部ａにて計測した振動応答測定点Ｂと励
振力着力点の間のインパルス応答関数を周波数領域のデ
ータに変換する。

【００１０】また、離散フーリエ逆変換器７ｂは、周波
数領域の時系列データに変換する機能を有するものであ
る。振動応答スペクトルを周波数応答関数で除すること
により、励振力推定値スペクトルが得られるが、この励
振力推定値スペクトルに離散フーリエ逆変換器７ｂによ
る演算を施すことにより、励振力推定値を時系列として
求める。

【００１１】

【発明が解決しようとする課題】上記従来装置で計測し
たデータの解析例を図１２及び図１３に示す。図１２及
び図１３はハンマリング時の機器１の周波数応答関数の
例で、図１２は振幅（加速度／力）を示し、図１３は位
相（ｒａｄ）を示しているが、機器の周波数応答関数は
多数の共振点（図１２で振幅がピークを与える部分）及
び反共振点（図１２では振幅が落ち込む部分）を有して
いる。従来技術を用いて周波数領域で励振力を正確に求
めるためには、これらの値を正確に計測する必要があ
る。特に反共振点では、振幅値がノイズレベルまで低下
する、あるいは計測機器のダイナミックレンジの制限を
受けるなどの理由で、正確な測定が困難であることが多
い。すなわち、反共振点においては、振幅がノイズある
いは計測器のダイナミックレンジの制限の影響で正しく
評価できない場合、その区間の位相値は不定となり、周
波数応答関数評価の上で大きな誤差要因となる。この誤
差を含んだ周波数応答関数を用いて励振力推定値を計算
すると、推定波形にはこの評価誤差に起因する誤差が生
ずる。

【００１２】本発明は上記の課題を解決するためになさ
れたもので、ノイズに対する安定性が強く、振幅、加振
力持続時間とも良好な推定結果が得られ、かつ、連続的
な励振力の推定にも適用できる励振力推定装置を提供す
ることを目的とする。

【００１３】

【課題を解決するための手段】本発明は、構造系の振動
応答測定点と励振力着力点との間のインパルス応答関数
を用いて構造系に負荷される励振力を推定する励振力推
定装置において、上記構造系の振動応答測定点の振動加
速度を検出する加速度検出器と、上記加速度検出器によ
り検出された加速度応答及び加振試験時に上記構造系に
加えられる励振力信号を記憶する第１の記憶手段と、上
記第１の記憶手段に記憶された加振試験時の励振力及び
加速度応答に基づいて上記構造系の振動応答測定点と励
振力着力点との間のインパルス応答関数を算出する第１
の演算手段と、上記第１の演算手段により計算されたイ
ンパルス応答関数と上記加速度検出器により検出される
機器稼働時の加速度応答を記憶する第２の記憶手段と、
上記第２の記憶手段に記憶された機器稼働時の加速度応
答を時分割して順次取り出すウタイムウィンドウ処理手
段と、上記タイムウィンドウ処理手段による各時分割出
力をそれぞれパルス状の励振力に対する過渡応答とみな
し、複素ケプストラム処理及びケフレンシフウィンドウ
処理により励振力を推定する第２の演算手段と、上記第
２の演算手段により推定された部分的な励振力を重ね合
わせて連続的な励振力推定値を得る重畳処理部とを具備
したことを特徴とする。

【００１４】上記第２の演算手段は、第２の記憶手段に
記憶されたインパルス応答関数及び機器稼働時の加速度
応答を周波数領域のデータに変換する離散フーリエ変換
器と、周波数応答関数と応答スペクトルとの複素対数を
とる複素対数変換器と、応答の対数スペクトルと対数周
波数応答関数に離散フーリエ逆変換を施し、複素ケプス
トラムに変換する離散フーリエ逆変換器と、上記複素ケ
プストラムに窓関数処理を施して短いケフレンシのデー
タだけを抽出する短ケフレンシ域抽出器と、上記窓関数
処理を施した応答とインパルス応答関数それぞれの複素
ケプストラムの演算により、励振力推定値の複素ケプス
トラムを求めた後、フーリエ変換により周波数領域に変
換する離散フーリエ変換器と、励振力スペクトルの複素
対数を真数に変換する複素指数変換器と、上記真数に変
換した励振力推定値スペクトルに離散フーリエ逆変換を
施して励振力推定値を時系列として求める離散フーリエ
逆変換器とからなることを特徴とする。

【００１５】（作用）加振試験時にインパルスハンマ等
を用いて構造系の励振力着力点に既知の荷重を加え、振
動応答測定点における加速度応答を加速度検出器により
検出して、加振試験における励振力及び加速度応答をメ
モリ装置に記憶する。第１の演算手段は、上記メモリ装
置に記憶した加振試験時の励振力及び加速度応答に基づ
いて構造系の振動応答測定点と励振力着力点との間のイ
ンパルス応答関数を計算する。第２の演算手段は、第１
の演算手段で計算したインパルス応答関数と上記加速度
検出器により検出される機器稼働時の加速度応答をケフ
レンシ領域のデータ（複素ケプストラム）に変換して構
造系に負荷される励振力を推定する。構造系の共振／反
共振の影響は、複素ケプストラムの長ケフレンシ領域に
現れるので、ショートパスケフレンシウィンドウ処理を
施してこれらの影響を除去する。

【００１６】連続的な励振力推定に対しては、タイムウ
ィンドウ処理手段により、連続波形を仮想的に衝撃的波
形が連続したものとして取り扱うことにより、上記と同
様の処理で励振力を推定することができる。すなわち、
振動応答を時間領域で細かく切り、それぞれをパルス状
の励振力に対する過渡応答とみなして第２の演算手段に
より励振力を推定した後、その推定した仮想的な衝撃波
形を重ね合わせて連続波形とする。

【００１７】

【発明の実施の形態】以下、図面を参照して本発明の一
実施形態を説明する。図１は、本発明の一実施形態に係
る励振力推定装置の構成を示すブロック図である。

【００１８】図１において、ａはインパルス応答関数計
測部、ｂは機器稼働時の振動応答計測部、ｃは励振力推
定演算部で、インパルス応答関数計測部ａで求めたイン
パルス応答関数及び機器稼働時の振動応答計測部ｂで計
測した機器稼働時の振動応答を用いて、励振力推定演算
部ｃで励振力推定値を計算する。

【００１９】上記インパルス応答関数計測部ａにおい
て、１１は構造系、つまり、被計測機器で、加振試験を
行なう際、励振力着力点Ａにインパルスハンマ１２によ
り既知の荷重（インパルス）が加えられる。また、被測
定機器１１には、振動応答測定点Ｂに加速度検出器１３
を設け、上記インパルスハンマ１２により荷重が加えら
れた際の加速度応答（加速度信号）を検出している。上
記インパルスハンマ１２の励振力及び加速度検出器１３
により検出した加速度応答は、第１のメモリ装置１４に
記憶され、第１の演算装置１５に送られる。この第１の
演算装置１５は、上記入力信号に基づいてインパルス応
答関数計算を行ない、この計算によって得られたインパ
ルス応答関数を第２のメモリ装置１６に出力する。

【００２０】また、この第２のメモリ装置１６には、被
測定機器１１の稼働時における加速度応答が加速度検出
器１３により検出されて記憶される。そして、第２のメ
モリ装置１６に記憶された機器稼働時の加速度応答及び
インパルス応答関数は、励振力推定演算部ｃに入力され
る。

【００２１】励振力推定演算部ｃは、第２の演算装置１
７、タイムウィンドウ処理ルーチン１８及び重畳処理部
１９からなり、第２のメモリ装置１６に記憶されたイン
パルス応答関数及び機器稼働時の加速度応答が入力され
る。

【００２２】第２の演算装置１７は、図２に示すように
離散フーリエ変換器２１、複素対数変換器２２、離散フ
ーリエ逆変換器２３、短ケフレンシ域抽出器２４、離散
フーリエ変換器２５、複素指数変換器２６、離散フーリ
エ逆変換器２７からなり、複素ケプストラム計算、ショ
ートパスケフレンシフウィンドウ処理により励振力推定
値を求める演算機能を備えている。

【００２３】上記第２の演算装置１７は、タイムウィン
ドウ処理ルーチン１８により振動応答を細かく切り、そ
れぞれをパルス状の励振力に対する過渡応答とみなして
複素ケプストラム計算、ショートパスケフレンシフウィ
ンドウ処理を行ない、励振力推定値を求めて重畳処理部
１９へ出力する。この重畳処理部１９は、第２の演算装
置１７の演算結果を重ね合わせ、連続した励振力推定値
として出力する。

【００２４】次に上記実施形態の動作を説明する。ま
ず、インパルス応答関数計測部ａにおいて、インパルス
ハンマ１２により被計測機器１１の励振力着力点Ａを叩
き、インパルスハンマ１２の励振力と加速度検出器３に
より検出した加速度応答とを第１のメモリ装置１４に一
時保存し、これを第１の演算装置１５にてインパルス応
答関数を計算する。このインパルス応答関数は、第２の
メモリ装置１６に記憶される。

【００２５】次に機器稼働時の振動応答計測部ｂは、被
計測機器１１の稼働時の加速度応答を加速度検出器１３
で検出し、その計測データを第２のメモリ装置１６に保
存する。

【００２６】そして、励振力推定演算部ｃは、第１の演
算装置１５にて計算されたインパルス応答関数と振動応
答計測部ｂで計測した加速度応答の計測データに基づい
て、第２の演算装置１７により、周波数領域においての
励振力推定計算を行なう。

【００２７】この場合、連続的な励振力推定に対して
は、タイムウィンドウ処理ルーチン１８を用いて、連続
波形を仮想的に衝撃的波形が連続したものとして取り扱
うことにより、励振力を推定することができる。

【００２８】上記第２の演算装置１７の演算動作は、次
のようにして行なわれる。第１の演算装置１５にて計算
されたインパルス応答関数と振動応答計測部ｂで計測し
た加速度応答の計測データは、離散フーリエ変換器２１
に入力されて処理され、続いて、複素対数変換器２２、
離散フーリエ逆変換器２３、短ケフレンシ域抽出器２
４、離散フーリエ変換器２５、複素指数変換器２６、離
散フーリエ逆変換器２７で順次処理され、この離散フー
リエ逆変換器２７から励振力推定値が出力される。この
場合、振動応答及びインパルス応答関数を複素ケプスト
ラム＾ｘ（τ）に変換して励振力の推定を行なうが、以
下にその詳細について説明する。上記複素ケプストラム
＾ｘ（τ）は、時系列ｘ（ｔ）のフーリエ変換の複素対
数の逆フーリエ変換で与えられる。定義式は、次に示す
ように

【００２９】

【数１】 Ｆ［…］：［…］のフーリエ変換で示される。

【００３０】複素ケプストラム（cepstrum）の名称は、
スペクトラム（spectrum）に逆変換を施したと言う観点
から付けられたものである。また、複素ケプストラム
は、時間と同じ次元の領域の関数であるが、特に周波数
（frequency ）領域の関数に逆フーリエ変換を施したも
のという観点から、ケフレンシ（quefrency ）領域の関
数として定義される。ケフレンシが“０”に近い領域は
時間と対応付けて短ケフレンシ領域、その逆の領域は長
ケフレンシ領域と称する。

【００３１】そして、振動応答及びインパルス応答関数
の複素ケプストラム＾ｘ（τ）を計算し、短ケフレンシ
領域の成分だけを抽出する窓関数処理を施す。周波数ス
ペクトルの共振／反共振の影響は長ケフレンシ領域に現
れるため、この窓関数処理（ケフレンシ領域における窓
関数をケフレンシウィンドウ、短ケフレンシ領域の成分
だけを抽出するケフレンシウィンドウをショートパスケ
フレンシウィンドウと称する）により、振動応答／イン
パルス応答関数双方のスペクトルへの共振／反共振の影
響を低減することができる。

【００３２】構造系の共振／反共振は、構造中を伝搬す
る振動の反射に起因するものであるから、この影響を低
減させることは励振力に対する直接応答のみを検出する
ことを意味し、従って、衝撃的励振力の推定に対して
は、精度の向上が期待できる。

【００３３】インパルス応答関数の複素ケプストラムを
計算する際、周波数応答関数の反共振点近傍ではノイズ
の影響が大きく位相値が不定となるため、複素ケプスト
ラム計算値にこれに起因する誤差が含まれる恐れがあ
る。この現象を回避するためインパルス応答関数は、最
小位相条件を満足すると仮定し、周波数応答関数の位相
は振幅値から定める。

【００３４】上記振動応答の複素ケプストラムは、入力
の複素ケプストラムとインパルス応答関数の複素ケプス
トラムとの和で与えられるので、ケフレンシウィンドウ
処理を施した振動応答の複素ケプストラムとの差を取る
ことにより、励振力波形推定値の複素ケプストラムを得
ることができる。そして、これを時系列に変換して励振
力波形推定値を得る。

【００３５】また、連続的な励振力推定に対しては、タ
イムウィンドウ処理ルーチン１８を用いて、連続波形を
衝撃的波形が連続したものとして取り扱うことにより、
同様の処理で励振力を推定することができる。図３は、
タイムウィンドウ処理の概念図を示したものである。

【００３６】図３（ａ）に示す振動応答を同図（ｂ）に
示すように時間領域で短く切り、それぞれを同図（ｃ）
に示すようにパルス状の励振力に対する過渡応答とみな
し、複素ケプストラム＋ケフレンシフウィンドウ処理に
より同図（ｄ）に示すように励振力を推定する。その
後、推定した仮想的な衝撃的波形を重畳処理部１９で重
ね合わせて同図（ｅ）に示す連続波形とする。

【００３７】次に上記の演算を行なう各演算器の処理機
能について説明する。離散フーリエ変換器２１は、機器
稼働時の加速度応答及び予めインパルス応答関数計測部
ａにて計測した応答測定点Ｂと励振力着力点Ａとの間の
インパルス応答関数ｈ′（ｔ）に離散フーリエ変換を施
して周波数領域のデータに変換する。応答は、時間領域
では励振力と系のインパルス応答の畳み込みで与えられ
るが（i.e. ｘ（ｔ）＝ｆ（ｔ）＊ｈ（ｔ））、周波数
領域に変換することにより、励振力スペクトルと周波数
応答関数の積で与えられるようになる（i.e. Ｘ（ω）
＝Ｆ（ω）×Ｈ（ω））。

【００３８】複素対数変換器２２は、応答スペクトルＸ
（ω）と、周波数応答関数スペクトルＨ′（ω）の複素
対数をとる。応答スペクトルの複素対数「＾Ｘ（ω）＝
logＸ（ω）＝log ｜Ｘ（ω）｜＋ｊφ」は、励振力ス
ペクトルの複素対数と周波数応答関数の複素対数の和と
なる（i.e. ＾Ｘ（ω）＝＾Ｆ（ω）＋＾Ｈ（ω））。

【００３９】離散フーリエ逆変換器２３は、応答の対数
スペクトル＾Ｘ（ω）と対数周波数応答関数＾Ｈ′
（ω）とに離散フーリエ逆変換を施し、ケフレンシ領域
のデータ（複素ケプストラム）に変換する。図４は、こ
の構造系のインパルス応答関数の複素ケプストラムの例
を示したものである。構造系の共振の影響は、長ケフレ
ンシ領域に現れる。

【００４０】短ケフレンシ域抽出器２４は、複素ケプス
トラムに短いケフレンシのデータだけ抽出する窓関数処
理を施す。窓関数の幅は、対象とする機器により変化さ
せる必要があるため、外部から入力可能とする。図５
は、この窓関数処理、つまり、ケフレンシウィンドウ処
理を施した複素ケプストラムの例を示したもので、図４
の短ケフレンシ領域のデータのみを抽出したものであ
る。

【００４１】離散フーリエ変換器２５は、窓関数処理を
施した応答とインパルス応答関数それぞれの複素ケプス
トラムの演算により、励振力推定値の複素ケプストラム
を求めた後、これにフーリエ変換を施し、周波数領域に
変換する。図６及び図７は、ケフレンシウィンドウ処理
を施した周波数応答関数の例、すなわち、図５の複素ケ
プストラムを周波数領域に変換したものであり、図６は
振幅、図７は位相を表している。

【００４２】複素指数変換器２６は、励振力推定値スペ
クトルの複素対数を真数に変換する。離散フーリエ逆変
換器２７は、励振力推定値スペクトルに離散フーリエ逆
変を施すことにより、励振力推定値を時系列として求め
る。

【００４３】上記のように第２の演算装置１７では、振
動応答及びインパルス応答関数を複素ケプストラムに変
換し、ケフレンシ領域で励振力の推定を行なう。また、
上記実施形態で示したように複素ケプストラム＋ショー
トパスケフレンシフウィンドウ処理を用いることによ
り、共振／反共振の影響を除去して直接応答のみを抽出
でき、衝撃的波形を同定する場合に有効である。更に、
タイムウインドウ処理により仮想的な衝撃力の推定を行
なうことによって、図４に示すような連続的な励振力波
形も図５に示す解析値と比較して良好な一致が見られ、
高い推定精度を得ることができる。なお、上記図４及び
図５は、ピストンスラップ力の波形を例として示したも
のである。

【００４４】

【発明の効果】以上詳記したように本発明によれば、周
波数応答関数の振幅をノイズレベルあるいは計測機器の
ダイナミックレンジの制限により正確に測定できない場
合においても、共振／反共振の影響をショートパスケフ
レンシウィンドウ処理により除去することにより、ノイ
ズに対する安定性が強く、振幅、加振力持続時間ともに
妥当な推定結果を得ることができる。従って、本発明に
よる励振力推定装置は、従来装置に比較して優れた機能
を有し、特に衝撃的加振力波形の同定に有効である。更
に、タイムウィンドウ処理により仮想的な衝撃力の推定
を行なうことによって、連続的な励振力の推定にも適用
することができる。

【図面の簡単な説明】

【図１】本発明の一実施形態に係る励振力推定装置の構
成を示すブロック図。

【図２】同実施形態における第２の演算装置の詳細を示
す構成図。

【図３】同実施形態におけるタイムウィンドウ処理の概
念図。

【図４】構造系のインパルス応答関数の複素ケプストラ
ムの例を示す図。

【図５】ケフレンシウィンドウ処理を施した複素ケプス
トラムの例を示す図。

【図６】ケフレンシウィンドウ処理を施した周波数応答
関数の振幅を示す図。

【図７】ケフレンシウィンドウ処理を施した周波数応答
関数の位相を示す図。

【図８】本発明の手法によって求めた励振力波形の例
（ピストンスラップ力）を示す図。

【図９】解析により求めた励振力波形の例（ピストンス
ラップ力）を示す図。

【図１０】従来の励振力推定装置の構成を示すブロック
図。

【図１１】従来の励振力推定装置の第２の演算装置の構
成を示すブロック図。

【図１２】最小位相条件を導入した構造系の周波数応答
関数の振幅を示す図。

【図１３】最小位相条件を導入した構造系の周波数応答
関数の位相を示す図。

【符号の説明】

１１ 被測定機器 １２ インパルスハンマ １３ 加速度検出器 １４ 第１のメモリ装置 １５ 第１の演算装置 １６ 第２のメモリ装置 １７ 第２の演算装置 １８ タイムウィンドウ処理ルーチン １９ 重畳処理部 ２１ 離散フーリエ変換器 ２２ 複素対数変換器 ２３ 離散フーリエ逆変換器 ２４ 短ケフレンシ域抽出器 ２５ 離散フーリエ変換器 ２６ 複素指数変換器 ２７ 離散フーリエ逆変換器

Claims (2)

【特許請求の範囲】
1. 【請求項１】 構造系の振動応答測定点と励振力着力点
   との間のインパルス応答関数を用いて構造系に負荷され
   る励振力を推定する励振力推定装置において、 上記構造系の振動応答測定点の振動加速度を検出する加
   速度検出器と、 上記加速度検出器により検出された加速度応答及び加振
   試験時に上記構造系に加えられる励振力信号を記憶する
   第１の記憶手段と、 上記第１の記憶手段に記憶された加振試験時の励振力及
   び加速度応答に基づいて上記構造系の振動応答測定点と
   励振力着力点との間のインパルス応答関数を算出する第
   １の演算手段と、 上記第１の演算手段により計算されたインパルス応答関
   数と上記加速度検出器により検出される機器稼働時の加
   速度応答を記憶する第２の記憶手段と、 上記第２の記憶手段に記憶された機器稼働時の加速度応
   答を時分割して順次取り出すウタイムウィンドウ処理手
   段と、 上記タイムウィンドウ処理手段による各時分割出力をそ
   れぞれパルス状の励振力に対する過渡応答とみなし、複
   素ケプストラム処理及びケフレンシフウィンドウ処理に
   より励振力を推定する第２の演算手段と、 上記第２の演算手段により推定された部分的な励振力を
   重ね合わせて連続的な励振力推定値を得る重畳処理部と
   を具備したことを特徴とする励振力推定装置。
2. 【請求項２】 第２の演算手段は、第２の記憶手段に記
   憶されたインパルス応答関数及び機器稼働時の加速度応
   答を周波数領域のデータに変換する離散フーリエ変換器
   と、周波数応答関数と応答スペクトルとの複素対数をと
   る複素対数変換器と、応答の対数スペクトルと対数周波
   数応答関数に離散フーリエ逆変換を施し、複素ケプスト
   ラムに変換する離散フーリエ逆変換器と、上記複素ケプ
   ストラムに窓関数処理を施して短いケフレンシのデータ
   だけを抽出する短ケフレンシ域抽出器と、上記窓関数処
   理を施した応答とインパルス応答関数それぞれの複素ケ
   プストラムの演算により、励振力推定値の複素ケプスト
   ラムを求めた後、フーリエ変換により周波数領域に変換
   する離散フーリエ変換器と、励振力スペクトルの複素対
   数を真数に変換する複素指数変換器と、上記真数に変換
   した励振力推定値スペクトルに離散フーリエ逆変換を施
   して励振力推定値を時系列として求める離散フーリエ逆
   変換器とからなることを特徴とする請求項１記載の励振
   力推定装置。