JP2967208B2 - Transient sound frequency analysis method and apparatus - Google Patents
Transient sound frequency analysis method and apparatusInfo
- Publication number
- JP2967208B2 JP2967208B2 JP7113914A JP11391495A JP2967208B2 JP 2967208 B2 JP2967208 B2 JP 2967208B2 JP 7113914 A JP7113914 A JP 7113914A JP 11391495 A JP11391495 A JP 11391495A JP 2967208 B2 JP2967208 B2 JP 2967208B2
- Authority
- JP
- Japan
- Prior art keywords
- power
- band
- frequency
- group
- expressed
- 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.)
- Expired - Lifetime
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION
【０００１】[0001]
【産業上の利用分野】本発明は、損傷発見のための打音
試験における打音の判定とか、生物の鳴き声の区別とか
に利用可能な過渡音の周波数とレベルとを解析する方法
及び同装置に関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method and apparatus for analyzing the frequency and level of a transient sound that can be used for determining a tapping sound in a tapping test for detecting damage and for distinguishing the sound of living things. It is about.
【０００２】[0002]
【従来の技術】従来、打音試験における打音や生物の鳴
き声の区別は、人の聴感によるか、連続した周波数分析
を行い、その都度得られる各周波数のレベルをスペクト
ログラフ（横軸に時間、縦軸に周波数を取り、各時刻の
周波数成分のレベルを濃淡または色相で表した表示図）
に表示することにより行われていた。2. Description of the Related Art Conventionally, a distinction between a tapping sound and a creature's sound in a tapping sound test is made based on human hearing or by performing continuous frequency analysis. , The vertical axis indicates the frequency, and the level of the frequency component at each time is represented by shading or hue.)
Was displayed.
【０００３】しかしながら、上記従来の音の特徴を判別
する手段では、周波数成分が広帯域に分布する過渡音の
場合は、表示されるデータ量が多大となり、人がそれを
見て音の特徴を判定するには表示が複雑過ぎるという問
題点がある。However, according to the above-described conventional means for determining the characteristics of sounds, in the case of transient sounds in which frequency components are distributed over a wide band, the amount of data to be displayed becomes large, and a person judges the characteristics of sounds by looking at it. There is a problem that the display is too complicated.
【０００４】[0004]
【発明が解決しようとする課題】解決しようとする課題
は、上記従来の音の特徴を判別する手段は、周波数成分
が広帯域に分布する過渡音の場合、表示されるデータ量
が多大となり、人がそれを見て音の特徴を判定するには
表示が複雑過ぎることである。この発明は、上記問題点
に着目したものであり、周波数成分が広帯域に分布する
過渡音に対して時刻毎にひとつの周波数とレベルの値の
組とを求める過渡音周波数解析方法及び同装置とを提供
することを目的としている。The problem to be solved is that the above-mentioned conventional means for discriminating the characteristics of sound requires a large amount of data to be displayed in the case of a transient sound in which frequency components are distributed over a wide band. However, the display is too complicated to determine the characteristics of the sound by looking at it. The present invention focuses on the above problem, and a transient sound frequency analysis method and a transient sound frequency analyzing method for finding one frequency and a set of level values at each time for transient sound in which frequency components are distributed over a wide band. It is intended to provide.
【０００５】[0005]
【課題を解決するための手段】請求項１の発明は、図１
に示すように、入力信号をそれぞれの通過周波数帯域が
互いに隣接しているバンドパスフィルタ群１０に加え、
各々バンドパスフィルタ出力を個別に二乗検波器群２０
により二乗検波し、さらにそれらの出力を平均器群３０
により個別に平滑化してパワーPi(i=1,2,..,N；N はバ
ンドの個数）を求めた後に、演算器５０により信号のパ
ワーPtをMeans for Solving the Problems The first aspect of the present invention is shown in FIG.
As shown in (1), the input signal is applied to a group of band-pass filters 10 whose respective pass frequency bands are adjacent to each other,
The output of each bandpass filter is individually set to a square detector group 20
, And their outputs are further averaged by the averager group 30.
Power P i and individually smoothed by (i = 1,2, .., N ; N is the band number of) after obtaining a power P t of the signal by the arithmetic unit 50
【０００６】[0006]
【数１】 (Equation 1)
【０００７】と求める。信号のレベルLtは音の基準パワ
ーPrとして Lt= 10 log(Pt/Pr) により求める。続いて、演算器４０により、対数で表し
た各バンドパスフィルタの中心周波数をFiとして、対数
で表した信号の周波数Ftを[0007] Level L t of the signal is determined by L t = 10 log based power P r of the sound (P t / P r). Subsequently, the arithmetic unit 40 sets the center frequency of each band-pass filter expressed by logarithm as F i , and calculates the frequency F t of the signal expressed by logarithm.
【０００８】[0008]
【数２】 (Equation 2)
【０００９】と求めることよりなる過渡音周波数解析方
法である。And a transient sound frequency analysis method.
【００１０】請求項２の発明は、図１に示すように、入
力信号からそれぞれの通過周波数帯域の成分を選択する
通過周波数帯域が互いに隣接しているバンドパスフィル
タ群１０、各々のバンドパスフィルタ出力を個別に二乗
検波する二乗検波器群２０、さらにそれらの出力を平滑
化してパワーPi(i=1,2,....N；N はバンドの個数）を求
める平均器群３０、それらのパワー信号から信号のパワ
ーPt及びレベルLtをA second aspect of the present invention is a band pass filter group 10 in which pass frequency bands for selecting components of each pass frequency band from an input signal are adjacent to each other, as shown in FIG. A group of square detectors 20 for individually square-detecting the outputs, and an averager group 30 for smoothing their outputs to obtain powers P i (i = 1, 2,... N; N is the number of bands); From those power signals, the signal power Pt and level Lt
【００１１】[0011]
【数３】 (Equation 3)
【００１２】及び Lt= 10 log(Pt/Pr) により求める演算器５０、続いて、対数で表した各バン
ドパスフィルタの中心周波数をFiとして、対数で表した
信号の周波数FtをAnd an arithmetic unit 50 obtained by L t = 10 log (P t / P r ), and then, assuming that the center frequency of each band-pass filter expressed by logarithm is Fi, the frequency Ft of the signal expressed by logarithm is obtained.
【００１３】[0013]
【数４】 (Equation 4)
【００１４】により求める演算器４０よりなる過渡音周
波数解析装置である。This is a transient sound frequency analysis device comprising an arithmetic unit 40 obtained by the above.
【００１５】[0015]
【実施例】実施例について図１によって説明すると、１
０はＮ個よりなるバンドバスフィルタ群であって、入力
信号からそれぞれの通過周波数帯域の成分を選択通過さ
せるものである。利用可能なフィルタとしては通過周波
数帯域が互いに隣接していればよく、例えば、定比フィ
ルタや高速フーリエ変換（ＦＦＴ）器が使用可能であ
る。２０は二乗検波器群であって前記バンドパスフィル
タ群１０のそれぞれの出力を二乗検波する。３０はそれ
らの出力を平滑化してパワーＰｉ（ｉ＝１，
２，．．．．，Ｎ；Ｎはバンドの個数）を求める平均器
群である。５０は演算器であって、それらのパワー信号
から入力信号のパワ−Ｐｔ及びレベルＬｔをDESCRIPTION OF THE PREFERRED EMBODIMENTS Referring to FIG.
Numeral 0 denotes a group of N band-pass filters for selectively passing the components of each pass frequency band from the input signal. It is sufficient that the pass frequency bands are adjacent to each other as an available filter. For example, a stoichiometric filter or a fast Fourier transform (FFT) device can be used. Reference numeral 20 denotes a square-law detector group, which square-detects the output of each of the band-pass filter groups 10. 30 smoothes their outputs and powers Pi (i = 1,
2,. . . . , N; N is the number of bands. 50 is a calculator, a power -P t and level L t input signals from their power signal
【００１６】[0016]
【数５】 (Equation 5)
【００１７】及び Lt=10 log(Pt/Pr) により求める。４０も演算器であって、対数で表した各
バンドパスフイルタの中心周波数をFiとして、対数で表
した入力信号の周波数FtをAnd L t = 10 log (P t / P r ). 40 is also an arithmetic unit, where F i is the center frequency of each bandpass filter expressed in logarithm, and the frequency F t of the input signal expressed in logarithm is
【００１８】[0018]
【数６】 (Equation 6)
【００１９】により求めるものである。[0019]
【００２０】[0020]
【作用】作用について、図２及び図３により説明する。
図２及び図３における破線及び実線は、それぞれ正常状
態及びその表面の一部に制振材を貼付した機械部品を叩
いて得られる打音に対して得られる本発明の解析結果を
横軸に時間を取って表したものである。両者を比較する
と、特に周波数の時間的変化が明確に異なっており、状
態の相違を容易に判定出来る結果が得られる。The operation will be described with reference to FIGS.
The dashed and solid lines in FIGS. 2 and 3 indicate the analysis results of the present invention obtained for hitting sound obtained by hitting a machine part having a normal state and a damping material attached to a part of its surface, respectively, on the horizontal axis. It is represented by taking time. When the two are compared, the temporal change of the frequency is clearly different, and the result that the difference in the state can be easily determined is obtained.
【００２１】[0021]
【発明の効果】本発明は以上のように構成されるため、
レベル及び周波数成分が時間的に変化する過渡音、特に
周波数成分が広帯域に分布するような過渡音に対して、
音の特徴を表す単純な解析結果を得ることができ、極め
て有用である。Since the present invention is configured as described above,
For transient sounds whose level and frequency components change over time, especially for transient sounds whose frequency components are distributed over a wide band,
A simple analysis result showing the characteristics of the sound can be obtained, which is extremely useful.
【図１】本発明の実施例を示す機器構成図である。FIG. 1 is a device configuration diagram showing an embodiment of the present invention.
【図２】本発明による機械部品の打音解析例を示す周波
数特性図である。FIG. 2 is a frequency characteristic diagram showing an example of a hitting sound analysis of a mechanical component according to the present invention.
【図３】本発明による機械部品の打音解析例を示すレベ
ル特性図である。FIG. 3 is a level characteristic diagram showing an example of hitting sound analysis of a mechanical component according to the present invention.
１０ バンドパスフィルタ群 ２０ 二乗検波器群 ３０ 平均器群 ４０ 演算器 ５０ 演算器 Reference Signs List 10 band-pass filter group 20 square-law detector group 30 averager group 40 arithmetic unit 50 arithmetic unit
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭63−65500（ＪＰ，Ａ) 特開 昭56−66899（ＪＰ，Ａ) (58)調査した分野(Int.Cl.6，ＤＢ名) G10L 3/00 - 9/20 G01H 17/00 G01R 23/165 ＪＩＣＳＴファイル（ＪＯＩＳ)────────────────────────────────────────────────── ─── Continuation of front page (56) References JP-A-63-65500 (JP, A) JP-A-56-66899 (JP, A) (58) Fields investigated (Int. Cl. 6 , DB name) G10L 3/00-9/20 G01H 17/00 G01R 23/165 JICST file (JOIS)
Claims (2)
数帯域が互いに隣接しているバンドパスフィルタ群（１
０）に加え、各々バンドパスフィルタ出力を個別に二乗
検波器群（２０）により二乗検波し、さらにそれらの出
力を平均器群（３０）により個別に平滑化して、各バン
ドについてパワーPi(i=1,2,..,N；Nはバンドの個数）を
求めた後に、演算器（５０）により、パワーPiの総和Pt
を求め、対数で表した、そのパワーPiの総和の基準パワ
ーPrに対する比で定義される信号のレベルLtを求めるこ
と、及び演算器（４０）により、対数で表した各バンド
パスフィルタ群（１０）の中心周波数をFiとして、対数
で表した、前記パワーPiと中心周波数Fiとの積の総和を
Ptで除した商で定義されるパワー信号の周波数Ftを求め
ることよりなる過渡音周波数解析方法。A band-pass filter group (1) having a pass frequency band adjacent to each other for an input signal of a transient sound.
0), the output of each band-pass filter is individually squared detected by a group of square detectors (20), and their outputs are individually smoothed by an averager group (30) to obtain a power P i ( i = 1,2, .., N; N is the sought after band number of), the computing unit (50), the power P i of total P t
The calculated and expressed in log, to determine the level L t of the signal is defined as the ratio relative to the reference power P r of the sum of the power P i, and the calculator by (40), each band pass filter expressed in logarithmic The center frequency of the group (10) is defined as F i, and the total sum of the product of the power P i and the center frequency F i expressed in logarithm is
Transient sound frequency analysis method consists in determining the frequency F t of the power signal that is defined by the quotient obtained by dividing the P t.
過周波数帯域の成分を選択通過させる、通過周波数帯域
が互いに隣接しているバンドパスフィルタ群（１０）、
各々バンドパスフィルタ出力を個別に二乗検波する二乗
検波器群（２０）、さらにそれらの出力を個別に平滑化
して、各バンドについてパワーPi(i=1,2,..,N；N はバ
ンドの個数）を求める平均器群（３０）、パワーPiの総
和Ptを求め、対数で表した、そのパワーPiの総和の基準
パワーPrに対する比で定義される信号のレベルLtを求め
る演算器（５０）、及び対数で表した各バンドパスフィ
ルタ群（１０）の中心周波数をFiとして、対数で表し
た、前記パワーPiと中心周波数Fiとの積の総和をPtで除
した商で定義されるパワー信号の周波数Ftを求める演算
器（４０）よりなる過渡音周波数解析装置。2. A band-pass filter group (10), which selectively passes respective pass frequency band components of the input signal of the added transient sound and whose pass frequency bands are adjacent to each other,
A group of square detectors (20) for individually square-detecting the output of each band-pass filter, and further smoothing their outputs individually, the power P i (i = 1, 2,..., N; the average unit group to obtain the band number of) (30), the power P i obtains the sum P t of, expressed in logarithmic level L t of the signal is defined as the ratio relative to the reference power P r of the sum of the power P i the determined arithmetic unit (50), and each band-pass filter group expressed in logarithmic center frequency of (10) as F i, expressed in logarithmic, the sum of product of the power P i and the center frequency F i P calculator for determining the frequency F t of the power signal that is defined by dividing the quotient by the t (40) transient sound frequency analyzer consisting.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7113914A JP2967208B2 (en) | 1995-04-14 | 1995-04-14 | Transient sound frequency analysis method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7113914A JP2967208B2 (en) | 1995-04-14 | 1995-04-14 | Transient sound frequency analysis method and apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08286699A JPH08286699A (en) | 1996-11-01 |
| JP2967208B2 true JP2967208B2 (en) | 1999-10-25 |
Family
ID=14624361
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7113914A Expired - Lifetime JP2967208B2 (en) | 1995-04-14 | 1995-04-14 | Transient sound frequency analysis method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2967208B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100339886C (en) * | 2003-04-10 | 2007-09-26 | 联发科技股份有限公司 | Coding device capable of detecting transient position of sound signal and its coding method |
| CN103743470B (en) * | 2013-12-23 | 2016-05-18 | 广西科技大学 | A kind of automobile noise frequency spectrum analysis method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5918718B2 (en) * | 1979-11-02 | 1984-04-28 | 株式会社東芝 | Speech analysis/synthesis system |
| JPS6365500A (en) * | 1986-09-08 | 1988-03-24 | 松下電器産業株式会社 | Friction property detector |
-
1995
- 1995-04-14 JP JP7113914A patent/JP2967208B2/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08286699A (en) | 1996-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Klein et al. | Spectral and modulation indices for annoyance-relevant features of urban road single-vehicle pass-by noises | |
| Karjalainen | A new auditory model for the evaluation of sound quality of audio systems | |
| US8065115B2 (en) | Method and system for identifying audible noise as wind noise in a hearing aid apparatus | |
| US7508948B2 (en) | Reverberation removal | |
| JPH11506970A (en) | Method and apparatus for characterizing gastrointestinal sounds | |
| DE69937789T2 (en) | METHOD AND DEVICE FOR DETERMINING THE RHYTHM FOR OUTCULTURE SIGNALS | |
| Durlach et al. | Informational masking for simultaneous nonspeech stimuli: Psychometric functions for fixed and randomly mixed maskers | |
| EP0186936A1 (en) | Method and apparatus for evaluating auditory distortions of an audio system | |
| Lee et al. | Development of a model to predict the likelihood of complaints due to assorted tone-in-noise combinations | |
| JP2967208B2 (en) | Transient sound frequency analysis method and apparatus | |
| Bruijn | Timbre-classification of complex tones | |
| JPH06201744A (en) | Device and method for measuring distortion of audio signal | |
| Olson | The measurement of loudness | |
| JP3350713B2 (en) | Method, apparatus and medium for identifying type of noise source | |
| Sakurai et al. | A computational software for noise measurement and toward its identification | |
| Torija et al. | Subjective dominance as a basis for selecting frequency weightings | |
| Blommer et al. | A new method for calculating Fluctuation Strength in electric motors | |
| JP4590545B2 (en) | Acoustic evaluation method and system | |
| Sottek et al. | Sound quality evaluation of fan noise based on hearing-related parameters | |
| JP2776307B2 (en) | Signal display device | |
| JP2784197B2 (en) | Frequency analysis method and apparatus | |
| JP3205109B2 (en) | Abnormal diagnosis device for rotating body | |
| JP2999868B2 (en) | Signal detection method | |
| JPH063474B2 (en) | Signal analyzer | |
| JP3148466B2 (en) | Device for discriminating between helicopter sound and vehicle sound |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |