JP3344864B2 - 残響時間自動測定方法および残響時間自動測定装置 - Google Patents
残響時間自動測定方法および残響時間自動測定装置Info
- Publication number
- JP3344864B2 JP3344864B2 JP05273895A JP5273895A JP3344864B2 JP 3344864 B2 JP3344864 B2 JP 3344864B2 JP 05273895 A JP05273895 A JP 05273895A JP 5273895 A JP5273895 A JP 5273895A JP 3344864 B2 JP3344864 B2 JP 3344864B2
- Authority
- JP
- Japan
- Prior art keywords
- time
- reverberation time
- impulse response
- noise
- curve
- 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 - Fee Related
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、室内音響特性を表す音
響エネルギー減衰曲線の傾きを精度良く推定することの
できる残響時間自動測定方法および残響時間自動測定装
置に関する。
響エネルギー減衰曲線の傾きを精度良く推定することの
できる残響時間自動測定方法および残響時間自動測定装
置に関する。
【０００２】［発明の概要］本発明は、室内の音響伝達
特性を表すインパルス応答から残響時間を算出するため
に必要な音響エネルギー減衰曲線を求める際に、低い周
波数帯域等で十分な信号対雑音比（以後、Ｓ／Ｎとい
う）が得られない場合でも、音響エネルギー減衰がノイ
ズの影響を受けてしまう以前の部分を観測区間として採
用することによって、音響エネルギーの減衰の傾きを推
定するという処理を自動的に行う方法に関するものであ
る。これまで、インパルス応答から残響時間を求める
際、測定現場で評価されるＳ／Ｎは、中域、高域におい
て十分であっても、特に低い周波数帯域において十分で
ない場合が多く、このような場合、減衰曲線がノイズの
影響を受けて折れ曲がりを生じ、本来の減衰から離れて
しまうため、正確な残響時間が求められないケースがあ
った。これに対して、得られたインパルス応答の無信号
の部分からノイズの自乗平均値を求めて減じることによ
り、信号成分を抜き出すという改善策が提案されている
が、十分にＳ／Ｎが改善されないため、得られた減衰曲
線のうちでノイズの影響が少ない部分を目で確認し、残
響時間を求め直す必要があった。本発明では、直接音到
来時刻以降の時々刻々のＳ／Ｎ（時間平均Ｓ／Ｎ）を観
測することによって、信号に対するノイズの影響をモニ
タし、時間平均Ｓ／Ｎが０ｄＢ以上という基準をもっ
て、ノイズの影響が少ない部分を観測区間とし、その部
分の最小自乗近似を行うことでエネルギー減衰曲線の傾
きを推定し、効率的に高精度な残響時間を求めるように
したものである。
特性を表すインパルス応答から残響時間を算出するため
に必要な音響エネルギー減衰曲線を求める際に、低い周
波数帯域等で十分な信号対雑音比（以後、Ｓ／Ｎとい
う）が得られない場合でも、音響エネルギー減衰がノイ
ズの影響を受けてしまう以前の部分を観測区間として採
用することによって、音響エネルギーの減衰の傾きを推
定するという処理を自動的に行う方法に関するものであ
る。これまで、インパルス応答から残響時間を求める
際、測定現場で評価されるＳ／Ｎは、中域、高域におい
て十分であっても、特に低い周波数帯域において十分で
ない場合が多く、このような場合、減衰曲線がノイズの
影響を受けて折れ曲がりを生じ、本来の減衰から離れて
しまうため、正確な残響時間が求められないケースがあ
った。これに対して、得られたインパルス応答の無信号
の部分からノイズの自乗平均値を求めて減じることによ
り、信号成分を抜き出すという改善策が提案されている
が、十分にＳ／Ｎが改善されないため、得られた減衰曲
線のうちでノイズの影響が少ない部分を目で確認し、残
響時間を求め直す必要があった。本発明では、直接音到
来時刻以降の時々刻々のＳ／Ｎ（時間平均Ｓ／Ｎ）を観
測することによって、信号に対するノイズの影響をモニ
タし、時間平均Ｓ／Ｎが０ｄＢ以上という基準をもっ
て、ノイズの影響が少ない部分を観測区間とし、その部
分の最小自乗近似を行うことでエネルギー減衰曲線の傾
きを推定し、効率的に高精度な残響時間を求めるように
したものである。
【０００３】
【従来の技術】従来より、室内の音響伝達特性を表すイ
ンパルス応答から音響エネルギー減衰曲線を求め、この
音響エネルギー減衰曲線に基づいて残響時間を算出する
ことが行われている。
ンパルス応答から音響エネルギー減衰曲線を求め、この
音響エネルギー減衰曲線に基づいて残響時間を算出する
ことが行われている。
【０００４】ところが、インパルス応答から残響時間を
求める際、測定現場で評価されるＳ／Ｎは、中域、高域
において十分であっても、特に低い周波数帯域において
十分でない場合が多い。このような場合、減衰曲線がノ
イズの影響を受けて折れ曲がりを生じ、本来の減衰から
離れてしまい、正確な残響時間が求められないケースが
あった。
求める際、測定現場で評価されるＳ／Ｎは、中域、高域
において十分であっても、特に低い周波数帯域において
十分でない場合が多い。このような場合、減衰曲線がノ
イズの影響を受けて折れ曲がりを生じ、本来の減衰から
離れてしまい、正確な残響時間が求められないケースが
あった。
【０００５】Ｓ／Ｎが十分でない音響エネルギー減衰曲
線から残響時間を求める際に、減衰曲線の観測区間が問
題になることがこれまでにも論じられており、ノイズの
エネルギーの自乗平均値を無信号部分から求め、得られ
た音響エネルギー減衰曲線から減じることによりＳ／Ｎ
の改善を図り、正確な残響時間が算出できるという報告
がある（Ｗ．Ｔ．Ｃｈｕ，Ｊ．Ａ．Ｓ．Ａ．Ｖｏｌ．６
３，Ｎｏ．５ １４４４−１４５９（１９７８））。
線から残響時間を求める際に、減衰曲線の観測区間が問
題になることがこれまでにも論じられており、ノイズの
エネルギーの自乗平均値を無信号部分から求め、得られ
た音響エネルギー減衰曲線から減じることによりＳ／Ｎ
の改善を図り、正確な残響時間が算出できるという報告
がある（Ｗ．Ｔ．Ｃｈｕ，Ｊ．Ａ．Ｓ．Ａ．Ｖｏｌ．６
３，Ｎｏ．５ １４４４−１４５９（１９７８））。
【０００６】しかし、上記改善策でも十分にＳ／Ｎが改
善されないため、得られた減衰曲線の観測区間のうちで
ノイズの影響が少ない部分を人間が目で確認し、残響時
間を求め直すという作業が必要であった。
善されないため、得られた減衰曲線の観測区間のうちで
ノイズの影響が少ない部分を人間が目で確認し、残響時
間を求め直すという作業が必要であった。
【０００７】
【発明が解決しようとする課題】上述したように、イン
パルス応答から求められる音響エネルギー減衰曲線は、
低い周波数帯域等においてＳ／Ｎが十分でない状態で得
られる場合が多く、音響エネルギー減衰曲線はノイズの
影響により折れ曲がりが生じ、本来の減衰から離れてし
まうため、正確な残響時間を求めるためには得られた音
響エネルギー減衰曲線からノイズの影響がない適切な部
分を評価する方法が課題とされ、最適な観測区間を定量
的に評価して自動的に決定する方法の開発が切望されて
いる。
パルス応答から求められる音響エネルギー減衰曲線は、
低い周波数帯域等においてＳ／Ｎが十分でない状態で得
られる場合が多く、音響エネルギー減衰曲線はノイズの
影響により折れ曲がりが生じ、本来の減衰から離れてし
まうため、正確な残響時間を求めるためには得られた音
響エネルギー減衰曲線からノイズの影響がない適切な部
分を評価する方法が課題とされ、最適な観測区間を定量
的に評価して自動的に決定する方法の開発が切望されて
いる。
【０００８】本発明は上記の事情に鑑みてなされたもの
であり、その目的は、高効率でかつ高精度に残響時間を
自動測定できる残響時間自動測定方法および残響時間自
動測定装置を提供することにある。
であり、その目的は、高効率でかつ高精度に残響時間を
自動測定できる残響時間自動測定方法および残響時間自
動測定装置を提供することにある。
【０００９】
【課題を解決するための手段】上記の目的を達成するた
めに請求項１記載の発明は、測定対象空間の音響測定に
おいて、インパルス応答から音響エネルギー減衰曲線を
求めて残響時間を算出するにあたり、前記インパルス応
答における直接音到来時刻から時々刻々のＳ／Ｎ（時間
平均Ｓ／Ｎ）が最初に０ｄＢになる時刻までの時間帯に
おける前記音響エネルギー減衰曲線の最小自乗近似線の
傾斜から残響時間を算出することを特徴としている。
めに請求項１記載の発明は、測定対象空間の音響測定に
おいて、インパルス応答から音響エネルギー減衰曲線を
求めて残響時間を算出するにあたり、前記インパルス応
答における直接音到来時刻から時々刻々のＳ／Ｎ（時間
平均Ｓ／Ｎ）が最初に０ｄＢになる時刻までの時間帯に
おける前記音響エネルギー減衰曲線の最小自乗近似線の
傾斜から残響時間を算出することを特徴としている。
【００１０】また、請求項２記載の発明は、測定対象と
なる空間におけるインパルス応答から音響エネルギー減
衰曲線を求めて残響時間を算出する装置において、前記
インパルス応答における時間平均Ｓ／Ｎを求める手段
と、インパルス応答における直接音到来時刻から前記時
間平均Ｓ／Ｎが最初に０ｄＢになる時刻までの時間帯を
検出する手段と、検出された時間帯における音響エネル
ギー減衰曲線の最小自乗近似線の傾斜を算出する手段
と、求められた傾斜から残響時間を算出する手段と、を
具備することを特徴としている。
なる空間におけるインパルス応答から音響エネルギー減
衰曲線を求めて残響時間を算出する装置において、前記
インパルス応答における時間平均Ｓ／Ｎを求める手段
と、インパルス応答における直接音到来時刻から前記時
間平均Ｓ／Ｎが最初に０ｄＢになる時刻までの時間帯を
検出する手段と、検出された時間帯における音響エネル
ギー減衰曲線の最小自乗近似線の傾斜を算出する手段
と、求められた傾斜から残響時間を算出する手段と、を
具備することを特徴としている。
【００１１】
【作用】上記の構成の本発明では、インパルス応答にお
ける直接音到来時刻から時々刻々のＳ／Ｎ（時間平均Ｓ
／Ｎ）を観測して信号成分の音響エネルギー減衰曲線に
対するノイズの影響をモニタする。そして、インパルス
応答における直接音到来時刻から前記時間平均Ｓ／Ｎが
最初に０ｄＢになる時刻までの時間帯を検出する。この
時間帯がノイズの影響が少ない部分であるから観測区間
とし、この観測区間（時間帯）における音響エネルギー
減衰曲線の最小自乗近似線の傾斜を算出（推定）し、求
められた傾斜から残響時間を算出するようにしている。
ける直接音到来時刻から時々刻々のＳ／Ｎ（時間平均Ｓ
／Ｎ）を観測して信号成分の音響エネルギー減衰曲線に
対するノイズの影響をモニタする。そして、インパルス
応答における直接音到来時刻から前記時間平均Ｓ／Ｎが
最初に０ｄＢになる時刻までの時間帯を検出する。この
時間帯がノイズの影響が少ない部分であるから観測区間
とし、この観測区間（時間帯）における音響エネルギー
減衰曲線の最小自乗近似線の傾斜を算出（推定）し、求
められた傾斜から残響時間を算出するようにしている。
【００１２】
【実施例】図１は本発明に係る残響時間自動測定装置の
一実施例構成を示すブロック図、図２はこの残響時間自
動測定装置によって測定される空間の測定環境を示す説
明図である。
一実施例構成を示すブロック図、図２はこの残響時間自
動測定装置によって測定される空間の測定環境を示す説
明図である。
【００１３】残響時間自動測定装置１は、ノイズ部分検
出部２と、ノイズレベル算出部３と、入力信号レベル算
出部４と、減算部５と、除算部６と、時間平均Ｓ／Ｎ算
出部７と、０ｄＢ検出部８と、ノイズ除去部９と、減衰
曲線算出部１０と、観測区間決定部（ゲート部）１１
と、傾斜算出部１２と、残響時間算出部１３とを備え、
図２に示すように測定対象空間３０のインパルス応答を
入力してその残響時間を自動測定する。この測定対象空
間３０には、図２に示すように、インパルス応答測定用
信号発生器２０が設置されている。この装置２０から発
生されたインパルス応答測定用信号は、スピーカーＳＰ
を介して測定対象空間３０中に出力され、その直接音波
および反射音波はマイクＭＣにより収音され、アンプ・
フィルタ部４０を介して前記残響時間測定装置１に供給
されるようになっている。
出部２と、ノイズレベル算出部３と、入力信号レベル算
出部４と、減算部５と、除算部６と、時間平均Ｓ／Ｎ算
出部７と、０ｄＢ検出部８と、ノイズ除去部９と、減衰
曲線算出部１０と、観測区間決定部（ゲート部）１１
と、傾斜算出部１２と、残響時間算出部１３とを備え、
図２に示すように測定対象空間３０のインパルス応答を
入力してその残響時間を自動測定する。この測定対象空
間３０には、図２に示すように、インパルス応答測定用
信号発生器２０が設置されている。この装置２０から発
生されたインパルス応答測定用信号は、スピーカーＳＰ
を介して測定対象空間３０中に出力され、その直接音波
および反射音波はマイクＭＣにより収音され、アンプ・
フィルタ部４０を介して前記残響時間測定装置１に供給
されるようになっている。
【００１４】前記ノイズ部分検出部２は、入力されたイ
ンパルス応答における無信号部分をノイズ部分として検
出する。具体的には、インパルス応答測定用信号発生器
２０からの信号発生と同時に出力される取り込み開始信
号を入力し、この信号入力から、受音点のマイクＭＣに
直接音が到来するまでの遅れ時間部分を無信号部分とす
る。
ンパルス応答における無信号部分をノイズ部分として検
出する。具体的には、インパルス応答測定用信号発生器
２０からの信号発生と同時に出力される取り込み開始信
号を入力し、この信号入力から、受音点のマイクＭＣに
直接音が到来するまでの遅れ時間部分を無信号部分とす
る。
【００１５】ノイズレベル算出部３は、検出されたノイ
ズ部分の自乗平均を算出してノイズレベルを検出する。
ズ部分の自乗平均を算出してノイズレベルを検出する。
【００１６】入力信号レベル算出部４は、入力された前
記インパルス応答における信号部分において十数ミリか
ら数ミリ秒（例えば１６ｍｓ）という短時間毎に入力信
号の自乗平均値を求めて入力信号レベルを検出する。
記インパルス応答における信号部分において十数ミリか
ら数ミリ秒（例えば１６ｍｓ）という短時間毎に入力信
号の自乗平均値を求めて入力信号レベルを検出する。
【００１７】減算部５は、入力信号レベル算出部４で検
出された入力信号レベルからノイズレベル算出部３で検
出されたノイズレベルを減算し、また、除算部６は、減
算部５の減算出力とノイズレベル算出部３で検出された
ノイズレベルとの比を求め、これを時間平均Ｓ／Ｎ算出
部７に供給する。
出された入力信号レベルからノイズレベル算出部３で検
出されたノイズレベルを減算し、また、除算部６は、減
算部５の減算出力とノイズレベル算出部３で検出された
ノイズレベルとの比を求め、これを時間平均Ｓ／Ｎ算出
部７に供給する。
【００１８】時間平均Ｓ／Ｎ算出部７は、十数ミリから
数ミリ秒（例えば１６ｍｓ）毎に供給される比データ
を、直接音到来から入力信号が十分に減衰するまで求め
続けて時間平均Ｓ／Ｎ曲線を算出する。
数ミリ秒（例えば１６ｍｓ）毎に供給される比データ
を、直接音到来から入力信号が十分に減衰するまで求め
続けて時間平均Ｓ／Ｎ曲線を算出する。
【００１９】０ｄＢ検出部８は、時間平均Ｓ／Ｎ算出部
７で算出された時間平均Ｓ／Ｎが０ｄＢとなった時点を
求めて、観測区間決定部１１に供給する。
７で算出された時間平均Ｓ／Ｎが０ｄＢとなった時点を
求めて、観測区間決定部１１に供給する。
【００２０】一方、ノイズ除去部９は、入力されたイン
パルス応答から前記無信号部分を検出してこれをノイズ
として除去し、ノイズが除去された信号を減衰曲線算出
部１０に供給する。
パルス応答から前記無信号部分を検出してこれをノイズ
として除去し、ノイズが除去された信号を減衰曲線算出
部１０に供給する。
【００２１】減衰曲線算出部１０は、入力信号からノイ
ズが除去された信号より音響エネルギーの減衰曲線を求
める操作を実行する。
ズが除去された信号より音響エネルギーの減衰曲線を求
める操作を実行する。
【００２２】観測区間決定部（ゲート部）１１は、求め
られた減衰曲線の観測区間を決定する部分であり、０ｄ
Ｂ検出部８から出力される、直接音到来時刻から時間平
均Ｓ／Ｎが最初に０ｄＢになるまでの時間を観測区間と
するゲート信号により減衰曲線算出部１０から出力され
る減衰曲線を切り出す処理をする。
られた減衰曲線の観測区間を決定する部分であり、０ｄ
Ｂ検出部８から出力される、直接音到来時刻から時間平
均Ｓ／Ｎが最初に０ｄＢになるまでの時間を観測区間と
するゲート信号により減衰曲線算出部１０から出力され
る減衰曲線を切り出す処理をする。
【００２３】傾斜算出部１２は、観測区間決定部１１で
切り出された減衰曲線における最小自乗近似線の傾斜を
求める。
切り出された減衰曲線における最小自乗近似線の傾斜を
求める。
【００２４】残響時間算出部１３は、傾斜算出部１２で
求められた傾斜から残響時間を求める処理を実行する。
求められた傾斜から残響時間を求める処理を実行する。
【００２５】次に本実施例例における残響時間測定処理
について図１の機能ブロック図および図３のフローチャ
ートを参照しつつ系統的に説明する。
について図１の機能ブロック図および図３のフローチャ
ートを参照しつつ系統的に説明する。
【００２６】また、図４には本実施例で使用された観測
波形例が示されている。
波形例が示されている。
【００２７】＜ノイズレベル検出処理（ステップＳＴ
１）＞先ず、ノイズ部分検出部２によって入力されたイ
ンパルス応答における無信号部分がノイズ部分として検
出され、ノイズレベル算出部３によってこの無信号部分
からノイズのエネルギー自乗平均値が求められ、ノイズ
レベルが算出される。無信号部分は、計測時に取り込み
時間を長くとることにより作成できる。例えば、前述し
たようにインパルス応答測定用信号発生器２０から供給
される取り込み開始信号と同時に取り込みを開始すれ
ば、受信音のマイクＭＣに直接音が到来するまでの遅れ
時間が無信号部分となる。
１）＞先ず、ノイズ部分検出部２によって入力されたイ
ンパルス応答における無信号部分がノイズ部分として検
出され、ノイズレベル算出部３によってこの無信号部分
からノイズのエネルギー自乗平均値が求められ、ノイズ
レベルが算出される。無信号部分は、計測時に取り込み
時間を長くとることにより作成できる。例えば、前述し
たようにインパルス応答測定用信号発生器２０から供給
される取り込み開始信号と同時に取り込みを開始すれ
ば、受信音のマイクＭＣに直接音が到来するまでの遅れ
時間が無信号部分となる。
【００２８】＜信号レベル検出処理（ステップＳＴ２）
＞入力信号レベル算出部４では、入力されたインパルス
応答について十数ミリ秒から数ミリ秒（例えば１６ｍ
ｓ）という短時間毎に入力信号の自乗平均値が求められ
る。
＞入力信号レベル算出部４では、入力されたインパルス
応答について十数ミリ秒から数ミリ秒（例えば１６ｍ
ｓ）という短時間毎に入力信号の自乗平均値が求められ
る。
【００２９】＜時間平均Ｓ／Ｎ算出処理（ステップＳＴ
３）＞上記処理により求めた入力信号とノイズとの自乗
平均値の比が数ミリ秒から十数ミリ秒毎に求められ、こ
れが直接音の到来から入力信号が十分に減衰するまで、
繰り返し求められて時間平均Ｓ／Ｎ曲線が求められる
（時間平均Ｓ／Ｎ算出部７の処理）。この時間平均Ｓ／
Ｎは次式により求めることができる。また、このように
して求められた時間平均Ｓ／Ｎ曲線を図６に示す。
３）＞上記処理により求めた入力信号とノイズとの自乗
平均値の比が数ミリ秒から十数ミリ秒毎に求められ、こ
れが直接音の到来から入力信号が十分に減衰するまで、
繰り返し求められて時間平均Ｓ／Ｎ曲線が求められる
（時間平均Ｓ／Ｎ算出部７の処理）。この時間平均Ｓ／
Ｎは次式により求めることができる。また、このように
して求められた時間平均Ｓ／Ｎ曲線を図６に示す。
【００３０】
【数１】 ＜ノイズ除去処理（ステップＳＴ４）＞上記の処理に並
行して、ノイズ部分検出部２で切り出されたノイズ部分
と、得られたインパルス応答とからノイズ除去部９にお
いて前記従来方法によりノイズ除去がされる。
行して、ノイズ部分検出部２で切り出されたノイズ部分
と、得られたインパルス応答とからノイズ除去部９にお
いて前記従来方法によりノイズ除去がされる。
【００３１】＜減衰曲線算出処理（ステップＳＴ５）＞
ノイズが除去されると、減衰曲線算出部１０では、その
信号に対して減衰曲線を求める操作が行われる。
ノイズが除去されると、減衰曲線算出部１０では、その
信号に対して減衰曲線を求める操作が行われる。
【００３２】このようにして得られた減衰曲線を図５に
示す。
示す。
【００３３】＜観測区間決定処理（ステップＳＴ６）＞
観測区間決定部１１では、０ｄＢ検出部８から時間平均
Ｓ／Ｎが最初に０ｄＢになる時刻が通知されると、直接
音到来時刻からこの０ｄＢになるまでの時間を観測区間
とするゲート信号により減衰曲線算出部１０から出力さ
れる減衰曲線を切り出す処理が行われる。
観測区間決定部１１では、０ｄＢ検出部８から時間平均
Ｓ／Ｎが最初に０ｄＢになる時刻が通知されると、直接
音到来時刻からこの０ｄＢになるまでの時間を観測区間
とするゲート信号により減衰曲線算出部１０から出力さ
れる減衰曲線を切り出す処理が行われる。
【００３４】＜傾斜算出処理（ステップＳＴ７）＞求め
られた観測区間における減衰曲線について最小自乗近似
線の傾斜が求られる。
られた観測区間における減衰曲線について最小自乗近似
線の傾斜が求られる。
【００３５】＜残響時間算出処理（ステップＳＴ８）＞
求められた傾斜に基づいて残響時間が求められる。この
ようにして求められた残響時間が結果出力とされる。
求められた傾斜に基づいて残響時間が求められる。この
ようにして求められた残響時間が結果出力とされる。
【００３６】次に、以上の処理による作用効果を説明す
る。
る。
【００３７】減衰曲線算出処理で求められた図５に示す
減衰曲線において、曲線Ａは何も処理を施さないときの
生データ曲線であるのに対し、曲線Ｂは従来方法により
ノイズ除去を行った場合の減衰を示している。また、図
６は、時間平均Ｓ／Ｎ曲線を示しており、この時間平均
Ｓ／Ｎが最初に０ｄＢになる時点（本実施例では、約
０．３５秒）まで、すなわち、０〜０．３５秒を観測区
間として推定した減衰曲線が図５の直線Ｃである。曲線
Ｄは、参考用にあらかじめ設定した残響時間（１秒）の
減衰を示しているが、直線Ｃは曲線Ｄに非常に近いこと
がわかる。
減衰曲線において、曲線Ａは何も処理を施さないときの
生データ曲線であるのに対し、曲線Ｂは従来方法により
ノイズ除去を行った場合の減衰を示している。また、図
６は、時間平均Ｓ／Ｎ曲線を示しており、この時間平均
Ｓ／Ｎが最初に０ｄＢになる時点（本実施例では、約
０．３５秒）まで、すなわち、０〜０．３５秒を観測区
間として推定した減衰曲線が図５の直線Ｃである。曲線
Ｄは、参考用にあらかじめ設定した残響時間（１秒）の
減衰を示しているが、直線Ｃは曲線Ｄに非常に近いこと
がわかる。
【００３８】また、算出された残響時間を表１に示す。
この表からも分かるように、設定した値（１秒）との誤
差は人間の検知限である５％以内に収まっている。
この表からも分かるように、設定した値（１秒）との誤
差は人間の検知限である５％以内に収まっている。
【００３９】
【表１】 このように本実施例によれば、高精度な減衰曲線の観測
区間決定を自動で行うことができ、残響時間の測定を効
率良く行うことが可能となる。また、測定対象空間３０
の本来の減衰曲線自体が折れ曲がっている場合、従来通
りに人間が目視判断すると、ノイズによる折れ曲がりと
誤って判断する虞があるが、本実施例では、常に時間平
均Ｓ／Ｎを確認することによりノイズの影響か否かを自
動的に正確に判断でき、減衰曲線の観測区間を正しく評
価することが可能となる。
区間決定を自動で行うことができ、残響時間の測定を効
率良く行うことが可能となる。また、測定対象空間３０
の本来の減衰曲線自体が折れ曲がっている場合、従来通
りに人間が目視判断すると、ノイズによる折れ曲がりと
誤って判断する虞があるが、本実施例では、常に時間平
均Ｓ／Ｎを確認することによりノイズの影響か否かを自
動的に正確に判断でき、減衰曲線の観測区間を正しく評
価することが可能となる。
【００４０】
【発明の効果】以上説明したように各請求項記載の発明
によれば、インパルス応答から音響エネルギー減衰曲線
を求めて残響時間を算出するにあたり、前記インパルス
応答における直接音到来時刻から時々刻々のＳ／Ｎ（時
間平均Ｓ／Ｎ）が最初に０ｄＢになる時刻までの時間帯
における前記音響エネルギー減衰曲線の最小自乗近似線
の傾斜から残響時間を算出するようにしたので、ノイズ
の影響を考慮でき、人手を介在させることなく高効率で
かつ高精度に残響時間を自動測定することが可能とな
る。
によれば、インパルス応答から音響エネルギー減衰曲線
を求めて残響時間を算出するにあたり、前記インパルス
応答における直接音到来時刻から時々刻々のＳ／Ｎ（時
間平均Ｓ／Ｎ）が最初に０ｄＢになる時刻までの時間帯
における前記音響エネルギー減衰曲線の最小自乗近似線
の傾斜から残響時間を算出するようにしたので、ノイズ
の影響を考慮でき、人手を介在させることなく高効率で
かつ高精度に残響時間を自動測定することが可能とな
る。
【図１】本発明に係る残響時間自動測定装置の一実施例
構成を示すブロック図である。
構成を示すブロック図である。
【図２】本発明に係る残響時間自動測定装置が適用され
る測定環境の構成を示す説明図である。
る測定環境の構成を示す説明図である。
【図３】本発明の一実施例における処理手順を示すフロ
ーチャートである。
ーチャートである。
【図４】本発明の一実施例で使用される観測波形の一例
を示す波形図である。
を示す波形図である。
【図５】本発明の一実施例における減衰曲線の一例を示
す波形図である。
す波形図である。
【図６】本発明の一実施例における時間平均Ｓ／Ｎ曲線
の一例を示す波形図である。
の一例を示す波形図である。
１ 残響時間自動測定装置 ２ ノイズ部分検出部 ３ ノイズレベル算出部 ４ 入力信号レベル算出部 ５ 減算部 ６ 除算部 ７ 時間平均Ｓ／Ｎ算出部 ８ ０ｄＢ検出部８ ９ ノイズ除去部 １０ 減衰曲線算出部１０ １１ 観測区間決定部 １２ 傾斜算出部 １３ 残響時間算出部 ２０ インパルス応答測定用信号発生器 ３０ 測定対象空間
───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7，ＤＢ名) G01H 7/00
Claims (2)
- 【請求項１】 測定対象空間の音響測定において、イン
パルス応答から音響エネルギー減衰曲線を求めて残響時
間を算出するにあたり、前記インパルス応答における直
接音到来時刻から時々刻々のＳ／Ｎ（時間平均Ｓ／Ｎ）
が最初に０ｄＢになる時刻までの時間帯における前記音
響エネルギー減衰曲線の最小自乗近似線の傾斜から残響
時間を算出することを特徴とする残響時間自動測定方
法。 - 【請求項２】 測定対象となる空間におけるインパルス
応答から音響エネルギー減衰曲線を求めて残響時間を算
出する装置において、 前記インパルス応答における時間平均Ｓ／Ｎを求める手
段と、 インパルス応答における直接音到来時刻から前記時間平
均Ｓ／Ｎが最初に０ｄＢになる時刻までの時間帯を検出
する手段と、 検出された時間帯における音響エネルギー減衰曲線の最
小自乗近似線の傾斜を算出する手段と、 求められた傾斜から残響時間を算出する手段と、 を具備することを特徴とする残響時間自動測定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05273895A JP3344864B2 (ja) | 1995-03-13 | 1995-03-13 | 残響時間自動測定方法および残響時間自動測定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05273895A JP3344864B2 (ja) | 1995-03-13 | 1995-03-13 | 残響時間自動測定方法および残響時間自動測定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08247837A JPH08247837A (ja) | 1996-09-27 |
| JP3344864B2 true JP3344864B2 (ja) | 2002-11-18 |
Family
ID=12923277
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP05273895A Expired - Fee Related JP3344864B2 (ja) | 1995-03-13 | 1995-03-13 | 残響時間自動測定方法および残響時間自動測定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3344864B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ATE532346T1 (de) * | 2003-09-24 | 2011-11-15 | Hewlett Packard Development Co | Verfahren und kommunikationsgerät mit mitteln zur entstörung von audiosignalen |
| JP2012015969A (ja) * | 2010-07-05 | 2012-01-19 | Tdk Corp | 無線通信システム |
| JP2012015968A (ja) * | 2010-07-05 | 2012-01-19 | Tdk Corp | 無線通信機器 |
| JP6197662B2 (ja) * | 2014-01-22 | 2017-09-20 | ヤマハ株式会社 | 残響時間解析装置 |
| JP6886890B2 (ja) * | 2017-08-08 | 2021-06-16 | 株式会社竹中工務店 | 減衰時間分析方法、装置、及びプログラム |
| CN114791360B (zh) * | 2021-01-25 | 2024-08-16 | 广州汽车集团股份有限公司 | 性能衰减分析处理方法、装置及存储介质 |
-
1995
- 1995-03-13 JP JP05273895A patent/JP3344864B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08247837A (ja) | 1996-09-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2020508168A5 (ja) | ||
| KR950002442B1 (ko) | 오디오 신호처리 시스템을 검사하기 위한 장치 | |
| JP3344864B2 (ja) | 残響時間自動測定方法および残響時間自動測定装置 | |
| JP2004193782A (ja) | スピーカとマイクロホン間の音波伝搬時間測定方法およびその装置 | |
| JPH0718746B2 (ja) | 回転体の異常検出装置 | |
| US20180024184A1 (en) | Method for characterising a soft fault in a cable | |
| JP3298318B2 (ja) | ガラス割れ検出装置 | |
| CN107785025B (zh) | 基于房间脉冲响应重复测量的噪声去除方法及装置 | |
| US9407992B2 (en) | Estimation of reverberation decay related applications | |
| JP5077847B2 (ja) | 残響時間推定装置及び残響時間推定方法 | |
| EP1806593B1 (en) | Determination of the adequate measurement window for sound source localization in echoic environments | |
| JP3630393B2 (ja) | 異常箇所検出装置 | |
| CN116519112A (zh) | 一种基于激光测振的超声刀头振动能量效率评价方法 | |
| JP3132388B2 (ja) | 回転機械の異音検査方法およびその装置 | |
| JP3387331B2 (ja) | 異常信号検出装置および異常信号源検出装置 | |
| JP2004333199A (ja) | 異音判定装置及び異音判定方法 | |
| JP2981480B2 (ja) | 配管の詰まり診断方法および診断装置 | |
| JP3430041B2 (ja) | 可搬型漏水位置探知装置による漏水位置の探知方法 | |
| JPH04318700A (ja) | 閾値決定手段を有する車間距離測定装置 | |
| JP3144878B2 (ja) | 回転機械の異常音診断方法、その装置及び回転機械の製造ライン | |
| JP2790909B2 (ja) | パッシブソーナの広帯域信号受信方法 | |
| JPS6255605B2 (ja) | ||
| WO2024018697A1 (ja) | 心拍数検出システム、心拍数検出方法及びプログラム | |
| CN112449278B (zh) | 一种自动校准延时输出声音的方法、装置、设备及存储介质 | |
| JPH08189921A (ja) | 割れ検出装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090830 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100830 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110830 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120830 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |