JPH07209068A - 音源探査装置 - Google Patents
音源探査装置Info
- Publication number
- JPH07209068A JPH07209068A JP450194A JP450194A JPH07209068A JP H07209068 A JPH07209068 A JP H07209068A JP 450194 A JP450194 A JP 450194A JP 450194 A JP450194 A JP 450194A JP H07209068 A JPH07209068 A JP H07209068A
- Authority
- JP
- Japan
- Prior art keywords
- sound source
- microphones
- time
- wavelet
- sound
- 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.)
- Withdrawn
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】
【目的】 種々の機械装置から発生する騒音をマイクロ
フォンで検出して演算により音源を探査す装置に関し、
時間的に変動する音源の探査が精度良くできる。 【構成】 ３個のマイクロフォン１ａ，１ｂ，１ｃが音
源２の音を検出し、Ａ／Ｄ変換器４を介してデジタル信
号としてウェーブレット変換器５へ入力する。同変換器
５ではこれらの信号にウェーブレット変換を施し、周波
数特性成分と時刻歴特性成分を抽出して音源位置演算器
７へ入力すると共に表示装置６へ６ａ，６ｂ，６ｃの如
く表示する。演算器７は、３個のマイクロフォン配置関
係より音源の方向角を求めると共にウェーブレット変換
した信号の時刻成分から音源までの距離を算出して音源
の方向を探査するので時間的に変動する音源に対しても
精度の良い探査ができる。
フォンで検出して演算により音源を探査す装置に関し、
時間的に変動する音源の探査が精度良くできる。 【構成】 ３個のマイクロフォン１ａ，１ｂ，１ｃが音
源２の音を検出し、Ａ／Ｄ変換器４を介してデジタル信
号としてウェーブレット変換器５へ入力する。同変換器
５ではこれらの信号にウェーブレット変換を施し、周波
数特性成分と時刻歴特性成分を抽出して音源位置演算器
７へ入力すると共に表示装置６へ６ａ，６ｂ，６ｃの如
く表示する。演算器７は、３個のマイクロフォン配置関
係より音源の方向角を求めると共にウェーブレット変換
した信号の時刻成分から音源までの距離を算出して音源
の方向を探査するので時間的に変動する音源に対しても
精度の良い探査ができる。
Description
【０００１】
【産業上の利用分野】本発明は種々の機械装置から発生
する音を３個以上のマイクロフォンで検出して演算によ
り音源を探査する音源探査装置に関する。
する音を３個以上のマイクロフォンで検出して演算によ
り音源を探査する音源探査装置に関する。
【０００２】
【従来の技術】機械装置や工場等における騒音低減対策
は、まずその騒音発制限を正確に検出することが糸口と
なるものであって、この騒音源の位置を探査する方式と
して種々のものが開発され、実用化されているが、この
うち最も代表的なものとして音響ホログラフィ法があ
る。
は、まずその騒音発制限を正確に検出することが糸口と
なるものであって、この騒音源の位置を探査する方式と
して種々のものが開発され、実用化されているが、この
うち最も代表的なものとして音響ホログラフィ法があ
る。
【０００３】音響ホログラフィ法による音源探査の過程
を図５および図６により具体的に説明する。図５におい
て、騒音を発している供試体２１はほぼｕ−ｖ面内にあ
り、供試体２１から距離Ｚ隔てた位置の平面ｘ−ｙ内に
記録面２２がある。記録面２２内の走査用マイクロフォ
ン２３は縦横に移動することが可能である。さらに、記
録面２２内又はそとの任意点に固定マイクロフォン２４
が設置されている。図６において、走査用マイクロフォ
ン２３、固定マイクロフォン２４の信号は、それぞれロ
ーパスフィルタ２５，２５′を介してＦＦＴ装置（高速
フーリエ変換機）２６，２６′によりフーリエスペクト
ルに変換され、両者のクロススペクトル２７を演算する
ことにより記録面２２での周波数別の音圧の複素振幅分
布が得られる。これをホログラム２８と呼び、いま、次
のように表す。
を図５および図６により具体的に説明する。図５におい
て、騒音を発している供試体２１はほぼｕ−ｖ面内にあ
り、供試体２１から距離Ｚ隔てた位置の平面ｘ−ｙ内に
記録面２２がある。記録面２２内の走査用マイクロフォ
ン２３は縦横に移動することが可能である。さらに、記
録面２２内又はそとの任意点に固定マイクロフォン２４
が設置されている。図６において、走査用マイクロフォ
ン２３、固定マイクロフォン２４の信号は、それぞれロ
ーパスフィルタ２５，２５′を介してＦＦＴ装置（高速
フーリエ変換機）２６，２６′によりフーリエスペクト
ルに変換され、両者のクロススペクトル２７を演算する
ことにより記録面２２での周波数別の音圧の複素振幅分
布が得られる。これをホログラム２８と呼び、いま、次
のように表す。
【０００４】Ｑ（xm,yn)：但しｍ＝１、２、〜、Ｍ、ｎ
＝１、２、〜、Ｎとする。ここで、（xm,yn)は記録面で
のサンプル点での座標である。ホログラム情報Ｑ(xm,y
n) を用いて、記録面２２から任意の距離Ｚにある面
（u,v)の音波情報Ｋ（u,y)３０を得るには次の式（１）
のヘルムホルツ／キルヒホフの積分定理を用い計算機２
９により演算する。
＝１、２、〜、Ｎとする。ここで、（xm,yn)は記録面で
のサンプル点での座標である。ホログラム情報Ｑ(xm,y
n) を用いて、記録面２２から任意の距離Ｚにある面
（u,v)の音波情報Ｋ（u,y)３０を得るには次の式（１）
のヘルムホルツ／キルヒホフの積分定理を用い計算機２
９により演算する。
【０００５】
【数１】
【０００６】式（１）は畳み込み積分であるので次の式
となる。
となる。
【０００７】 Ｑ（x,y)＝Ｋ（u,y)＊Ｔ（u,v) （２） と表わせる。
【０００８】式（１）において、右辺Ｋ（u,v)を除いた
項を伝達関数Ｔとする。
項を伝達関数Ｔとする。
【０００９】
【数２】
【００１０】従って、フーリエ変換により式（２）は次
の様になる。
の様になる。
【００１１】 Ｆ〔Ｑ（x,y)〕＝Ｆ〔Ｋ（u,v)〕・Ｆ〔Ｔ（u,v)〕 （３） 式（３）により、Ｋ（u,v)は次の式（４）により求めら
れる。
れる。
【００１２】
【数３】
【００１３】この音響ホログラフィ法は高精度で音源位
置の探査・測定を行うことができるが、大規模マイクロ
フォン掃引装置や複雑な演算処理機器を必要とし、しか
もその演算処理に多大の時間をようするため、瞬間的騒
音、変動する騒音の音源位置の探査・測定が困難になる
等の問題があった。
置の探査・測定を行うことができるが、大規模マイクロ
フォン掃引装置や複雑な演算処理機器を必要とし、しか
もその演算処理に多大の時間をようするため、瞬間的騒
音、変動する騒音の音源位置の探査・測定が困難になる
等の問題があった。
【００１４】
【発明が解決しようとする課題】前述の様に従来の方法
では音源から発した信号をフーリエ変換により周波数領
域に変換するため、時間的に平均化された情報しか入手
することができず、しかも、その演算処理には多くの時
間を必要とし、このため、時間的に変動する音源の探索
はほとんど不可能であった。また音源が移動する様な場
合についても同様である。
では音源から発した信号をフーリエ変換により周波数領
域に変換するため、時間的に平均化された情報しか入手
することができず、しかも、その演算処理には多くの時
間を必要とし、このため、時間的に変動する音源の探索
はほとんど不可能であった。また音源が移動する様な場
合についても同様である。
【００１５】本発明は上記の時間的に変動し、しかも分
散性のある音源に対して３個のマイクロフォンを用いて
これらの検出した信号にウェーブレット変換を施して周
波数と時間を把握して音源の方向，距離の精度のよい探
査が可能な装置提供するものである。
散性のある音源に対して３個のマイクロフォンを用いて
これらの検出した信号にウェーブレット変換を施して周
波数と時間を把握して音源の方向，距離の精度のよい探
査が可能な装置提供するものである。
【００１６】
【課題を解決するための手段】本発明はそのために正三
角形に配置された３個のマイクロフォンの時間および空
間的に分散性のある騒音信号に対してウェーブレット変
換を行い、これらの信号から音源の方向角を演算すると
共に時間特性および周波数特性を分離した形で詳細に把
握することにより、各マイクロフォンから音源までの距
離を演算し、それぞれの距離から音源の方向を演算する
ことにより音源を正確に探索する装置を提供する。
角形に配置された３個のマイクロフォンの時間および空
間的に分散性のある騒音信号に対してウェーブレット変
換を行い、これらの信号から音源の方向角を演算すると
共に時間特性および周波数特性を分離した形で詳細に把
握することにより、各マイクロフォンから音源までの距
離を演算し、それぞれの距離から音源の方向を演算する
ことにより音源を正確に探索する装置を提供する。
【００１７】即ち、本発明は、正三角形に配置され、音
源の騒音波形を電気信号に変換するための３個のマイク
ロフォンと、これらマイクロフォンの出力信号をＡ／Ｄ
変換器を介して取り込み、該マイクロフォンの出力信号
にそれぞれウェーブレット変換を施すことにより周波数
特性成分と時刻歴特性成分を抽出するための処理を行う
ウェーブレット変換器と、その演算結果から得られる前
記マイクロフォンへの各周波数成分の到達時間、前記３
個のマイクロフォン間の距離及び配置角度との関係から
前記音源の入射角を含むその音源の距離、方向を演算す
る演算処理装置と、該演算処理装置の演算結果を表示す
る表示装置とを具備してなることを特徴とする音源探査
装置を提供する。
源の騒音波形を電気信号に変換するための３個のマイク
ロフォンと、これらマイクロフォンの出力信号をＡ／Ｄ
変換器を介して取り込み、該マイクロフォンの出力信号
にそれぞれウェーブレット変換を施すことにより周波数
特性成分と時刻歴特性成分を抽出するための処理を行う
ウェーブレット変換器と、その演算結果から得られる前
記マイクロフォンへの各周波数成分の到達時間、前記３
個のマイクロフォン間の距離及び配置角度との関係から
前記音源の入射角を含むその音源の距離、方向を演算す
る演算処理装置と、該演算処理装置の演算結果を表示す
る表示装置とを具備してなることを特徴とする音源探査
装置を提供する。
【００１８】
【作用】本発明は前述の手段により、３個のマイクロフ
ォンで検出した騒音信号はＡ／Ｄ変換器でデジタル信号
としてウェーブレット演算装置に入力される。これらの
入力された信号はそれぞれウェーブレット変換を施すこ
とにより周波数特性成分と時刻歴特性成分で把握される
信号となり演算処理装置に入力される。これらの入力さ
れた信号は３個のマイクロフォン信号から演算されたも
のであるため、これらのマイクロフォンの配置角度、マ
イクロフォン間の距離、等の関係より、各周波数ごとの
マイクロフォン間の到達時間差が求められ、これらの到
達時間差の正・負より音源からの方位角、が求められ
る。これと共に、ウェーブレット変換された信号には時
刻特性より音がマイクロフォンに到達した時刻が把握さ
れているので音速とで音源までの距離が演算され、音源
の正確な方向探査が可能となる。
ォンで検出した騒音信号はＡ／Ｄ変換器でデジタル信号
としてウェーブレット演算装置に入力される。これらの
入力された信号はそれぞれウェーブレット変換を施すこ
とにより周波数特性成分と時刻歴特性成分で把握される
信号となり演算処理装置に入力される。これらの入力さ
れた信号は３個のマイクロフォン信号から演算されたも
のであるため、これらのマイクロフォンの配置角度、マ
イクロフォン間の距離、等の関係より、各周波数ごとの
マイクロフォン間の到達時間差が求められ、これらの到
達時間差の正・負より音源からの方位角、が求められ
る。これと共に、ウェーブレット変換された信号には時
刻特性より音がマイクロフォンに到達した時刻が把握さ
れているので音速とで音源までの距離が演算され、音源
の正確な方向探査が可能となる。
【００１９】図２は測定地点に正三角形に配置された３
個のマイクロフォンＭ1 , Ｍ2 ，Ｍ 3 の配置図であり、
この例で探査の原理を説明する。いま、Ｌ：マイクロフ
ォン間隔、ｖ：音速、Ｍ1 〜Ｍ3 ：マイクロフォンＩ
ａ，Ｉｂ，Ｉｃ、θ：音源方向角（時計方向を正とす
る）、α：マイクロフォンの開き角、Ｄij：マイクロフ
ォンＭｉに対するマイクロフォンＭｊの遅れ、とすれ
ば、各マイクロフォンへの到達時間差Ｄijが判れば音源
方向角は従来の方法と同様に次のように求められる。
個のマイクロフォンＭ1 , Ｍ2 ，Ｍ 3 の配置図であり、
この例で探査の原理を説明する。いま、Ｌ：マイクロフ
ォン間隔、ｖ：音速、Ｍ1 〜Ｍ3 ：マイクロフォンＩ
ａ，Ｉｂ，Ｉｃ、θ：音源方向角（時計方向を正とす
る）、α：マイクロフォンの開き角、Ｄij：マイクロフ
ォンＭｉに対するマイクロフォンＭｊの遅れ、とすれ
ば、各マイクロフォンへの到達時間差Ｄijが判れば音源
方向角は従来の方法と同様に次のように求められる。
【００２０】 Ｄ12＝（Ｌ／ｖ）cos （α−θ） cosφ （５） Ｄ13＝（Ｌ／ｖ）cos （α＋θ） cosφ （６） 式（５）（６）を展開して差および和をとれば、 Ｄ12−Ｄ13＝２（Ｌ／ｖ） sinα sinθ cosφ （７） Ｄ12＋Ｄ13＝２（Ｌ／ｖ） cosα cosθ cosφ （８） 式（７）を式（８）で割ると、
【００２１】
【数４】
【００２２】（１０）あるいは（１０′）よりＤ12，Ｄ
13の正負を考慮して音源方向角θが求められる。音源迎
角φはθを（５）あるいは（６）式に用いることにより
求められる。
13の正負を考慮して音源方向角θが求められる。音源迎
角φはθを（５）あるいは（６）式に用いることにより
求められる。
【００２３】次に本発明のポイントである到達時間差Ｄ
12，Ｄ13の評価法について説明する。ウェーブレット変
換は、時間的にも周波数的にも局在した関数Ψ（ｔ）に
スケール変換ａ（＞０）とシフト変換ｂを作用させて得
られる関数Ψ((ｔ−ｂ）／a)を積分核とする積分変換で
あり、入力信号をｆ（ｔ）とすると次式で表わされる。
12，Ｄ13の評価法について説明する。ウェーブレット変
換は、時間的にも周波数的にも局在した関数Ψ（ｔ）に
スケール変換ａ（＞０）とシフト変換ｂを作用させて得
られる関数Ψ((ｔ−ｂ）／a)を積分核とする積分変換で
あり、入力信号をｆ（ｔ）とすると次式で表わされる。
【００２４】
【数５】
【００２５】ここで上つきバーは複素共役を表わし、Ψ
（ｔ）は次式の条件を満たす関数である。
（ｔ）は次式の条件を満たす関数である。
【００２６】
【数６】
【００２７】ここで、Ψ（ω）はΨ（ｔ）のフーリエ変
換を表わす。
換を表わす。
【００２８】Ψ（ｔ）としては次式の Morlet wavelet
(Gabor関数：Ψ１（ｔ））およびMexican Hat (Lapac
ian Gaussian関数：Ψ (ｔ) ）等を用いることができ
る。
(Gabor関数：Ψ１（ｔ））およびMexican Hat (Lapac
ian Gaussian関数：Ψ (ｔ) ）等を用いることができ
る。
【００２９】
【数７】
【００３０】図３に示すような分散波形について上記の
ウェーブレット変換を行った場合の結果を図４に示す。
図４においてｂ軸は時刻を、ａ軸は周波数を示してい
る。ａ−ｂ平面において変換値がピークを示す線ｇは音
がマイクロフォンに到達した時刻を表わしており、この
変換結果と音速により音源までの距離を推定することが
可能となる。
ウェーブレット変換を行った場合の結果を図４に示す。
図４においてｂ軸は時刻を、ａ軸は周波数を示してい
る。ａ−ｂ平面において変換値がピークを示す線ｇは音
がマイクロフォンに到達した時刻を表わしており、この
変換結果と音速により音源までの距離を推定することが
可能となる。
【００３１】
【実施例】以下、本発明を図面に示す実施例に基づいて
具体的に説明する。図１は本発明の一実施例を示す音源
探査装置のシステムの構成図で、１ａ，１ｂ，１ｃ，は
前述のように正三角形に配置された音源探査用のマイク
ロフォン、２は音源、３ａ，３ｂ，３ｃはそれぞれ３個
のマイクロフォンに対応するアンプ、４はマイクロフォ
ン３ａ，３ｂ，３ｃの信号をデジタル信号に変換するた
めのＡ／Ｄ変換器、５はそれらのデジタル信号をウェー
ブレット変換するために演算するウェーブレット変換
器、６は３個のマイクロフォン１ａ，１ｂ，１ｃからの
信号のウェーブレット変換結果を表示するための表示装
置で６ａ，６ｂ，６ｃのように表示する。７はウェーブ
レット変換結果を基にマイクロフォン１ａ，１ｂ，１ｃ
への音の到達時間差から音源の方向位置を演算するため
の演算器、この演算結果は６の表示装置で表示される。
具体的に説明する。図１は本発明の一実施例を示す音源
探査装置のシステムの構成図で、１ａ，１ｂ，１ｃ，は
前述のように正三角形に配置された音源探査用のマイク
ロフォン、２は音源、３ａ，３ｂ，３ｃはそれぞれ３個
のマイクロフォンに対応するアンプ、４はマイクロフォ
ン３ａ，３ｂ，３ｃの信号をデジタル信号に変換するた
めのＡ／Ｄ変換器、５はそれらのデジタル信号をウェー
ブレット変換するために演算するウェーブレット変換
器、６は３個のマイクロフォン１ａ，１ｂ，１ｃからの
信号のウェーブレット変換結果を表示するための表示装
置で６ａ，６ｂ，６ｃのように表示する。７はウェーブ
レット変換結果を基にマイクロフォン１ａ，１ｂ，１ｃ
への音の到達時間差から音源の方向位置を演算するため
の演算器、この演算結果は６の表示装置で表示される。
【００３２】ウェーブレット変換は、前述のようり時間
的にも周波数的にも局在した関数Ψ（ｔ）にスケール変
換ａ（＞０）とシフト変換ｂを作用させて得られる関数
Ψ（（ｔ−ｂ）／ａ）を積分核とする積分変換であり、
入力信号をｆ（ｔ）とすると前述の式（１１）で表わさ
れる。
的にも周波数的にも局在した関数Ψ（ｔ）にスケール変
換ａ（＞０）とシフト変換ｂを作用させて得られる関数
Ψ（（ｔ−ｂ）／ａ）を積分核とする積分変換であり、
入力信号をｆ（ｔ）とすると前述の式（１１）で表わさ
れる。
【００３３】図３に示すような分散波について上記のウ
ェーブレット変換を行った場合の結果は前述の通り、図
４のようになり、図４においてｂ軸は時刻を、ａ軸は周
波数を示している、ａ−ｂ平面において変換値がピーク
を示す線ｇは音がマイクロフォンに到達した時刻を表わ
しており、音源位置演算器７でこの変換結果と音速によ
り音源までの距離を推定することが可能となる。
ェーブレット変換を行った場合の結果は前述の通り、図
４のようになり、図４においてｂ軸は時刻を、ａ軸は周
波数を示している、ａ−ｂ平面において変換値がピーク
を示す線ｇは音がマイクロフォンに到達した時刻を表わ
しており、音源位置演算器７でこの変換結果と音速によ
り音源までの距離を推定することが可能となる。
【００３４】
【発明の効果】以上、具体的に説明したように従来音源
の探査はマイクロフォンでとらえた信号をＦＦＴ等の手
法を用いて周波数領域に変換して時間的に平均化された
探査情報を入手していた。この方法では時間的に変動す
る音源の場合、移動する音源の場合その探査精度が極端
に低下するという問題があった。そこで本発明では、マ
イクロフォンでとらえた信号を周波数のみでなく時間的
にも信号として把握可能なウェーブレット変換により変
換しその結果をもとに音源の方向距離を同定する装置と
したために時間的に変動する音源、移動する音源の探査
が精度良く求めることができる。
の探査はマイクロフォンでとらえた信号をＦＦＴ等の手
法を用いて周波数領域に変換して時間的に平均化された
探査情報を入手していた。この方法では時間的に変動す
る音源の場合、移動する音源の場合その探査精度が極端
に低下するという問題があった。そこで本発明では、マ
イクロフォンでとらえた信号を周波数のみでなく時間的
にも信号として把握可能なウェーブレット変換により変
換しその結果をもとに音源の方向距離を同定する装置と
したために時間的に変動する音源、移動する音源の探査
が精度良く求めることができる。
【図１】本発明の一実施例に係る音源探査装置の構成図
である。
である。
【図２】本発明のマイクロフォンの配置を示す図であ
る。
る。
【図３】本発明の音源探査装置に入力してウェーブレッ
ト変換を行う分散波形（ひずみ応答）の一例を示す波形
図である。
ト変換を行う分散波形（ひずみ応答）の一例を示す波形
図である。
【図４】本発明の音源探査装置でウェーブレット変換を
施した結果を示す斜視図である。
施した結果を示す斜視図である。
【図５】従来の音源探査の方法を示す説明図である。
【図６】従来の音源探査装置の構成を示すブロック図で
ある。
ある。
１ａ マイクロフォン １ｂ マイクロフォン １ｃ マイクロフォン ２ 音源 ４ Ａ／Ｄ変換器 ５ ウェーブレット変換器 ６ 表示装置 ７ 音源位置演算装置
Claims (1)
- 【請求項１】 正三角形に配置され、音源の騒音波形を
電気信号に変換するための３個のマイクロフォンと、こ
れらマイクロフォンの出力信号をＡ／Ｄ変換器を介して
取り込み該マイクロフォンの出力信号にそれぞれウェー
ブレット変換を施すことにより周波数特性成分と時刻歴
特性成分を抽出するための処理を行うウェーブレット変
換器と、その演算結果から得られる前記マイクロフォン
への各周波数成分の到達時間、前記３個のマイクロフォ
ン間の距離及び配置角度との関係から前記音源の入射角
を含むその音源の距離、方向を演算する演算処理装置
と、該演算処理装置の演算結果を表示する表示装置とを
具備してなることを特徴とする音源探査装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP450194A JPH07209068A (ja) | 1994-01-20 | 1994-01-20 | 音源探査装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP450194A JPH07209068A (ja) | 1994-01-20 | 1994-01-20 | 音源探査装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH07209068A true JPH07209068A (ja) | 1995-08-11 |
Family
ID=11585813
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP450194A Withdrawn JPH07209068A (ja) | 1994-01-20 | 1994-01-20 | 音源探査装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07209068A (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998011417A1 (fr) * | 1996-09-13 | 1998-03-19 | Kabushiki Kaisha Toshiba | Analyseur de signal instable et support d'enregistrement de programme d'analyse de signal instable |
| JP2000214052A (ja) * | 1999-01-28 | 2000-08-04 | Nichiha Corp | 異常音検出システム及び記録媒体 |
| JP2003083803A (ja) * | 2001-09-17 | 2003-03-19 | Sato Kogyo Co Ltd | 騒音・振動の監視低減方法およびそのための騒音・振動モニタリング装置 |
| KR20030089528A (ko) * | 2002-05-15 | 2003-11-22 | 현대자동차주식회사 | 차량의 소음 발생위치 확인장치 및 그 방법 |
| JP2010169430A (ja) * | 2009-01-20 | 2010-08-05 | Nec System Technologies Ltd | 電子機器、および電子機器における機構部品の診断方法、ならびに診断プログラム |
| TWI450602B (zh) * | 2012-06-06 | 2014-08-21 | Nat Univ Tsing Hua | 微型電子式獵槍麥克風 |
-
1994
- 1994-01-20 JP JP450194A patent/JPH07209068A/ja not_active Withdrawn
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998011417A1 (fr) * | 1996-09-13 | 1998-03-19 | Kabushiki Kaisha Toshiba | Analyseur de signal instable et support d'enregistrement de programme d'analyse de signal instable |
| JP2000214052A (ja) * | 1999-01-28 | 2000-08-04 | Nichiha Corp | 異常音検出システム及び記録媒体 |
| JP2003083803A (ja) * | 2001-09-17 | 2003-03-19 | Sato Kogyo Co Ltd | 騒音・振動の監視低減方法およびそのための騒音・振動モニタリング装置 |
| KR20030089528A (ko) * | 2002-05-15 | 2003-11-22 | 현대자동차주식회사 | 차량의 소음 발생위치 확인장치 및 그 방법 |
| JP2010169430A (ja) * | 2009-01-20 | 2010-08-05 | Nec System Technologies Ltd | 電子機器、および電子機器における機構部品の診断方法、ならびに診断プログラム |
| TWI450602B (zh) * | 2012-06-06 | 2014-08-21 | Nat Univ Tsing Hua | 微型電子式獵槍麥克風 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ding et al. | A new method for waveform analysis for estimating AE wave arrival times using wavelet decomposition | |
| Gaul et al. | Identification of the impact location on a plate using wavelets | |
| Jeong et al. | Fracture source location in thin plates using the wavelet transform of dispersive waves | |
| US6587798B2 (en) | Method and system for determining the speed of sound in a fluid within a conduit | |
| JP2958362B2 (ja) | 地盤構造の計測解析判定方法 | |
| US6360609B1 (en) | Method and system for interpreting and utilizing multimode dispersive acoustic guided waves | |
| CN101893698B (zh) | 噪声源测试分析方法及其装置 | |
| CN113358743B (zh) | 一种基于时频分布相似度分析的兰姆波模态分离方法 | |
| Park et al. | Localization of a sound source in a noisy environment by hyperbolic curves in quefrency domain | |
| JPH07209068A (ja) | 音源探査装置 | |
| JP2006058278A (ja) | 弾性波発生位置算出装置及び方法 | |
| Harley et al. | Broadband localization in a dispersive medium through sparse wavenumber analysis | |
| Wu et al. | Defect detection in pipes using Van der Pol systems based on ultrasonic guided wave | |
| CN110045014B (zh) | 基于贝叶斯学习的Lamb波频散消除方法及其系统 | |
| Lee et al. | Identification of impact force on a thick plate based on the elastodynamic and higher-order time-frequency analysis | |
| Apostoloudia et al. | Time–frequency analysis of transient dispersive waves: A comparative study | |
| JP3639160B2 (ja) | 異常箇所検出装置 | |
| JP4519852B2 (ja) | 共振現象を利用した超音波探査方法およびその装置 | |
| CN107741263A (zh) | 基于频域信息融合的液位估计方法 | |
| JP2003114219A (ja) | 弾性波伝播速度測定演算方法及び該方法を用いた非破壊圧縮強度試験装置 | |
| JPH09178863A (ja) | 地盤強度解析方法 | |
| JPH07209069A (ja) | 音源探査装置 | |
| US9441993B2 (en) | Flow measurement systems and methods for gas and liquid applications | |
| Szwoch et al. | Detection of the incoming sound direction employing MEMS microphones and the DSP | |
| JP3639145B2 (ja) | 異常箇所検出装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20010403 |