JPH09292885A - 音響空間インパルス応答推定装置 - Google Patents
音響空間インパルス応答推定装置Info
- Publication number
- JPH09292885A JPH09292885A JP8108828A JP10882896A JPH09292885A JP H09292885 A JPH09292885 A JP H09292885A JP 8108828 A JP8108828 A JP 8108828A JP 10882896 A JP10882896 A JP 10882896A JP H09292885 A JPH09292885 A JP H09292885A
- Authority
- JP
- Japan
- Prior art keywords
- digital filter
- impulse response
- sound
- wave
- direct wave
- 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
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 8
- 230000001902 propagating effect Effects 0.000 description 6
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 101150087426 Gnal gene Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000000034 method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Two-Way Televisions, Distribution Of Moving Picture Or The Like (AREA)
- Reverberation, Karaoke And Other Acoustics (AREA)
- Telephonic Communication Services (AREA)
Abstract
(57)【要約】
【課題】 積和演算回数の少ない音響空間インパルス応
答推定装置を実現する。 【解決手段】 入力端子１０から入力されたインパルス
に対し、ディジタルフィルタ１６は、直接波と反射波の
伝搬時間とレベルを設定する。ディジタルフィルタ１３
は、直接波のみのインパルス応答から伝搬時間を除いた
成分を特性として持ち、ディジタルフィルタ１６の出力
信号とその成分との積和演算を行うことで、直接波と複
数の反射波のインパルス応答を求める。
答推定装置を実現する。 【解決手段】 入力端子１０から入力されたインパルス
に対し、ディジタルフィルタ１６は、直接波と反射波の
伝搬時間とレベルを設定する。ディジタルフィルタ１３
は、直接波のみのインパルス応答から伝搬時間を除いた
成分を特性として持ち、ディジタルフィルタ１６の出力
信号とその成分との積和演算を行うことで、直接波と複
数の反射波のインパルス応答を求める。
Description
【０００１】
【発明の属する技術分野】本発明は、音声会議或いはテ
レビ会議（以下、ＴＶ会議という）の音声マンマシーン
インタフェースや、バーチャルリアリティの音声処理等
に使用される音響空間インパルス応答を推定する音響空
間インパルス応答推定装置に関するものである。
レビ会議（以下、ＴＶ会議という）の音声マンマシーン
インタフェースや、バーチャルリアリティの音声処理等
に使用される音響空間インパルス応答を推定する音響空
間インパルス応答推定装置に関するものである。
【０００２】
【従来の技術】音声会議やＴＶ会議における音響のステ
レオ再生に際し、まるでコンサートホールで音楽を聞い
ているような臨場感を出すために、電気的に反射音や残
響音を操作することが行われている。この操作を実現す
るためには、ホールや部屋での音響空間インパルス応答
を求める必要がある。つまり、音源のスピーカから放射
したインパルスが、マイクにどのような波形で到達する
かを求める必要がある。従来では、ＤＳＰ（Digital Si
gnal Processor；ディジタル信号処理プロセッサ）等を
用いた音響空間インパルス応答推定装置で、そのインパ
ルス応答を推定していた。
レオ再生に際し、まるでコンサートホールで音楽を聞い
ているような臨場感を出すために、電気的に反射音や残
響音を操作することが行われている。この操作を実現す
るためには、ホールや部屋での音響空間インパルス応答
を求める必要がある。つまり、音源のスピーカから放射
したインパルスが、マイクにどのような波形で到達する
かを求める必要がある。従来では、ＤＳＰ（Digital Si
gnal Processor；ディジタル信号処理プロセッサ）等を
用いた音響空間インパルス応答推定装置で、そのインパ
ルス応答を推定していた。
【０００３】
【発明が解決しようとする課題】しかしながら、従来の
音響空間インパルス応答推定装置では、次のような課題
があった。通常の部屋でのインパルス応答の長さは、高
次の反射波の影響で０．２〜０．６秒に達する。このよ
うな高次の反射を処理しようとすると、演算量が増加す
るという課題があった。以下、その理由を説明する。図
２は、音の反射波と直接波を幾何学的に説明する図であ
り、部屋における直接波と反射波を示している。部屋１
内に配置されたスピーカ２から放射した音波は、種々の
経路ｋｄ，ｋ１，ｋ２，ｋ３，ｋ４，ｋ５，ｋ６を通過
してマイク３に達する。経路ｋｄは、スピーカ２から直
接マイク３に入力される直接波の経路である。各経路ｋ
１，ｋ２を伝搬する音波は、部屋１の壁で２度反射して
マイク３にそれぞれ到達する。各経路ｋ３，ｋ４を伝搬
する音波は、１度反射してマイク３に到達する。各経路
ｋ５，ｋ６を伝搬する音波も、壁で２度反射してマイク
３に到達する。つまり、各経路ｋ１〜ｋ６を伝搬してマ
イク３に達した音波は、反射波である。経路ｋ３，ｋ４
を伝搬してマイク３に達した音波は、一次反射波であ
り、経路ｋ１，ｋ２，ｋ５，ｋ６を伝搬してマイク３に
達した音波は二次反射波である。図２中のｚ１，ｚ２，
ｚ３，ｚ４，ｚ５，ｚ６は、音波をマイク３側から見た
スピーカ２の虚像である。つまり、マイク３には、虚像
ｚ１〜ｚ６の位置にスピーカ２が配置してあるように、
音波が入力されている。
音響空間インパルス応答推定装置では、次のような課題
があった。通常の部屋でのインパルス応答の長さは、高
次の反射波の影響で０．２〜０．６秒に達する。このよ
うな高次の反射を処理しようとすると、演算量が増加す
るという課題があった。以下、その理由を説明する。図
２は、音の反射波と直接波を幾何学的に説明する図であ
り、部屋における直接波と反射波を示している。部屋１
内に配置されたスピーカ２から放射した音波は、種々の
経路ｋｄ，ｋ１，ｋ２，ｋ３，ｋ４，ｋ５，ｋ６を通過
してマイク３に達する。経路ｋｄは、スピーカ２から直
接マイク３に入力される直接波の経路である。各経路ｋ
１，ｋ２を伝搬する音波は、部屋１の壁で２度反射して
マイク３にそれぞれ到達する。各経路ｋ３，ｋ４を伝搬
する音波は、１度反射してマイク３に到達する。各経路
ｋ５，ｋ６を伝搬する音波も、壁で２度反射してマイク
３に到達する。つまり、各経路ｋ１〜ｋ６を伝搬してマ
イク３に達した音波は、反射波である。経路ｋ３，ｋ４
を伝搬してマイク３に達した音波は、一次反射波であ
り、経路ｋ１，ｋ２，ｋ５，ｋ６を伝搬してマイク３に
達した音波は二次反射波である。図２中のｚ１，ｚ２，
ｚ３，ｚ４，ｚ５，ｚ６は、音波をマイク３側から見た
スピーカ２の虚像である。つまり、マイク３には、虚像
ｚ１〜ｚ６の位置にスピーカ２が配置してあるように、
音波が入力されている。
【０００４】図３は、図２中のマイク３に入力される音
波とインパルス波形を示す図である。この図３では、経
路ｋｄを伝搬した直接波Ｓｄと、経路ｋ３，ｋ４を伝搬
した一次反射波Ｓ３，Ｓ４と、経路ｋ１，ｋ２を伝搬し
た二次反射波Ｓ１，Ｓ２と、経路ｋ５，ｋ６を伝搬した
二次反射波Ｓ５，Ｓ６とが、示されている。なお、簡単
のため、経路ｋ１と経路ｋ２、経路ｋ３と経路ｋ４、及
び経路ｋ５と経路ｋ６は対称であり、音波の伝搬距離は
等しいものとしている。直接波Ｓｄは、時間ｔ0 後にマ
イク３に到達する。各一次反射波Ｓ３，Ｓ４は、時間ｔ
１後にマイク３にそれぞれ到達する。各二次反射波Ｓ
１，Ｓ２は、時間ｔ２後にマイク３にそれぞれ到達し、
各二次反射波Ｓ５，Ｓ６は、時間ｔ３後にマイク３にそ
れぞれ到達する。これらの反射波Ｓ１，Ｓ２と、反射波
Ｓ３，Ｓ４と反射波Ｓ５，Ｓ６は、スピーカ２からマイ
ク３までの経路上の距離が異なるので、該マイク３への
到着時間が異なる。また、これら反射波Ｓ１〜Ｓ６は、
壁で反射する度に、音のエネルギの一部が吸収されて減
衰する。これらの直接波Ｓｄ、反射波Ｓ１〜Ｓ６が合成
されて、インパルス波形が得られる。これがインパルス
応答である。図２及び図３では二次反射までしか示して
いないが、実際にはもっと高次の反射まで存在するの
で、インパルス波形の継続時間はもっと長時間になり、
通常の部屋では０．２〜０．６秒になる。
波とインパルス波形を示す図である。この図３では、経
路ｋｄを伝搬した直接波Ｓｄと、経路ｋ３，ｋ４を伝搬
した一次反射波Ｓ３，Ｓ４と、経路ｋ１，ｋ２を伝搬し
た二次反射波Ｓ１，Ｓ２と、経路ｋ５，ｋ６を伝搬した
二次反射波Ｓ５，Ｓ６とが、示されている。なお、簡単
のため、経路ｋ１と経路ｋ２、経路ｋ３と経路ｋ４、及
び経路ｋ５と経路ｋ６は対称であり、音波の伝搬距離は
等しいものとしている。直接波Ｓｄは、時間ｔ0 後にマ
イク３に到達する。各一次反射波Ｓ３，Ｓ４は、時間ｔ
１後にマイク３にそれぞれ到達する。各二次反射波Ｓ
１，Ｓ２は、時間ｔ２後にマイク３にそれぞれ到達し、
各二次反射波Ｓ５，Ｓ６は、時間ｔ３後にマイク３にそ
れぞれ到達する。これらの反射波Ｓ１，Ｓ２と、反射波
Ｓ３，Ｓ４と反射波Ｓ５，Ｓ６は、スピーカ２からマイ
ク３までの経路上の距離が異なるので、該マイク３への
到着時間が異なる。また、これら反射波Ｓ１〜Ｓ６は、
壁で反射する度に、音のエネルギの一部が吸収されて減
衰する。これらの直接波Ｓｄ、反射波Ｓ１〜Ｓ６が合成
されて、インパルス波形が得られる。これがインパルス
応答である。図２及び図３では二次反射までしか示して
いないが、実際にはもっと高次の反射まで存在するの
で、インパルス波形の継続時間はもっと長時間になり、
通常の部屋では０．２〜０．６秒になる。
【０００５】図４は、従来の音響空間インパルス応答推
定装置を示す構成ブロック図である。入力端子１０は、
スピーカの出力するインパルスを入力する端子であり、
その入力端子１０には、複数の遅延素子１１−１，１１
−２，…，１１−ｎが接続されている。各遅延素子１１
−１〜１１−ｎは、直接波及び反射波の伝搬距離に対応
する遅延時間を設定する機能をそれぞれ有し、これら遅
延素子１１−１〜１１−ｎの数ｎは、直接波と反射波
（高次の反射波も含む）の数である。設定される遅延時
間は、図２の時間ｔ０，ｔ１，ｔ２，ｔ３，…に相当す
る。各遅延素子１１−１，１１−２，…１１−ｎの出力
側が、直接波及び各反射波のレベルを設定する乗算器１
２−１，１２−２，…，１２−ｎに、それぞれ接続され
ている。各乗算器１２−１，１２−２，…，１２−ｎの
出力側が、ディジタルフィルタ１３−１，１３−２，
…，１３−ｎにそれぞれ接続されている。ディジタルフ
ィルタ１３−１，１３−２，…，１３−ｎの出力側が累
算器１４に接続され、該累算器１４の出力側が出力端子
１５に接続されている。出力端子１５からの出力信号が
マイクに入力される構成である。各ディジタルフィルタ
１３−１〜１３−ｎの有する特性は同じであり、マイク
３へ到達する直接波の特性のうち、遅延時間を除いたも
のである。各ディジタルフィルタ１３−１〜１３−ｎで
は、入力された信号に対して、その特性で積和演算を行
う機能をそれぞれ有している。
定装置を示す構成ブロック図である。入力端子１０は、
スピーカの出力するインパルスを入力する端子であり、
その入力端子１０には、複数の遅延素子１１−１，１１
−２，…，１１−ｎが接続されている。各遅延素子１１
−１〜１１−ｎは、直接波及び反射波の伝搬距離に対応
する遅延時間を設定する機能をそれぞれ有し、これら遅
延素子１１−１〜１１−ｎの数ｎは、直接波と反射波
（高次の反射波も含む）の数である。設定される遅延時
間は、図２の時間ｔ０，ｔ１，ｔ２，ｔ３，…に相当す
る。各遅延素子１１−１，１１−２，…１１−ｎの出力
側が、直接波及び各反射波のレベルを設定する乗算器１
２−１，１２−２，…，１２−ｎに、それぞれ接続され
ている。各乗算器１２−１，１２−２，…，１２−ｎの
出力側が、ディジタルフィルタ１３−１，１３−２，
…，１３−ｎにそれぞれ接続されている。ディジタルフ
ィルタ１３−１，１３−２，…，１３−ｎの出力側が累
算器１４に接続され、該累算器１４の出力側が出力端子
１５に接続されている。出力端子１５からの出力信号が
マイクに入力される構成である。各ディジタルフィルタ
１３−１〜１３−ｎの有する特性は同じであり、マイク
３へ到達する直接波の特性のうち、遅延時間を除いたも
のである。各ディジタルフィルタ１３−１〜１３−ｎで
は、入力された信号に対して、その特性で積和演算を行
う機能をそれぞれ有している。
【０００６】各反射波は、直接波に類似な波形を遅延さ
せて減衰させた波形として表すことができる。反射波の
遅延時間は、反射経路を音速で割ったものであり、部屋
の幾何学的配置から推定できる。各遅延素子１１−１〜
１１−ｎは、入力端子１０から入力されたインパルスに
対して、その推定された遅延時間分だけ、それぞれ遅延
させる。各乗算器１２−１〜１２−ｎは、各遅延素子１
１−１〜１１−ｎの出力信号に対して、係数α0 〜αn
を掛けて直接波と各反射波のレベルをそれぞれ設定す
る。各ディジタルフィルタ１３−１〜１３−ｎは、前記
特性と乗算器１２−１〜１２−ｎの出力信号の積和を求
め、直接波と各反射波の波形を形成する。累算器１４の
累算によって、合成したインパルス応答のインパルス波
形が得られる。直接波の遅延時間を除いた信号の継続時
間は、１０〜２０ｍｓ程度である。例えば、音声を１６
ｋＨｚでサンプリングする場合、約１６０〜３２０サン
プルの積和演算が必要である。これ自身はそれほどの演
算量ではないが、反射波も含めると、ｎ倍の演算量が必
要である。例えば、図２の場合には、約６４０〜１２８
０サンプルの積和演算が必要になる。
せて減衰させた波形として表すことができる。反射波の
遅延時間は、反射経路を音速で割ったものであり、部屋
の幾何学的配置から推定できる。各遅延素子１１−１〜
１１−ｎは、入力端子１０から入力されたインパルスに
対して、その推定された遅延時間分だけ、それぞれ遅延
させる。各乗算器１２−１〜１２−ｎは、各遅延素子１
１−１〜１１−ｎの出力信号に対して、係数α0 〜αn
を掛けて直接波と各反射波のレベルをそれぞれ設定す
る。各ディジタルフィルタ１３−１〜１３−ｎは、前記
特性と乗算器１２−１〜１２−ｎの出力信号の積和を求
め、直接波と各反射波の波形を形成する。累算器１４の
累算によって、合成したインパルス応答のインパルス波
形が得られる。直接波の遅延時間を除いた信号の継続時
間は、１０〜２０ｍｓ程度である。例えば、音声を１６
ｋＨｚでサンプリングする場合、約１６０〜３２０サン
プルの積和演算が必要である。これ自身はそれほどの演
算量ではないが、反射波も含めると、ｎ倍の演算量が必
要である。例えば、図２の場合には、約６４０〜１２８
０サンプルの積和演算が必要になる。
【０００７】
【課題を解決するための手段】第１の発明は、前記課題
を解決するために、音源から受音点に到達する音波の経
路に、直接波の経路と複数の反射波の経路とが存在する
音響空間におけるインパルス応答を求める音響空間イン
パルス応答推定装置において、次のような構成にしてい
る。即ち、第１の発明の音響空間インパルス応答推定装
置は、入力されたインパルスに対し、前記音源から前記
受音点に到達する前記直接波と前記複数の反射波の伝搬
時間とそれらのレベルとを推定して設定する第１のディ
ジタルフィルタと、前記直接波のみのインパルス応答か
ら前記伝搬時間を除いた成分を特性として持ち、前記第
１のディジタルフィルタの出力信号と該特性との積和演
算を行うことで、前記直接波と前記複数の反射波のイン
パルス応答を求める第２のディジタルフィルタとを備え
ている。第２の発明は、第１の発明における第１のディ
ジタルフィルタは、主観受音として有効な分だけの前記
直接波と前記複数の反射波の伝搬時間とそれらのレベル
を推定して設定し、前記第２のディジタルフィルタにお
ける積和演算回数を制限構成にしている。第１，第２の
発明によれば、以上のように音響空間インパルス応答推
定装置を構成したので、第１のディジタルフィルタによ
り、前記直接波と前記複数の反射波の伝搬時間とそれら
のレベルが、１系統の信号として設定される。その第１
のディジタルフィルタの出力信号に対してのみ、第２の
ディジタルフィルタの積和演算が行われる。従って、前
記課題を解決できるのである。
を解決するために、音源から受音点に到達する音波の経
路に、直接波の経路と複数の反射波の経路とが存在する
音響空間におけるインパルス応答を求める音響空間イン
パルス応答推定装置において、次のような構成にしてい
る。即ち、第１の発明の音響空間インパルス応答推定装
置は、入力されたインパルスに対し、前記音源から前記
受音点に到達する前記直接波と前記複数の反射波の伝搬
時間とそれらのレベルとを推定して設定する第１のディ
ジタルフィルタと、前記直接波のみのインパルス応答か
ら前記伝搬時間を除いた成分を特性として持ち、前記第
１のディジタルフィルタの出力信号と該特性との積和演
算を行うことで、前記直接波と前記複数の反射波のイン
パルス応答を求める第２のディジタルフィルタとを備え
ている。第２の発明は、第１の発明における第１のディ
ジタルフィルタは、主観受音として有効な分だけの前記
直接波と前記複数の反射波の伝搬時間とそれらのレベル
を推定して設定し、前記第２のディジタルフィルタにお
ける積和演算回数を制限構成にしている。第１，第２の
発明によれば、以上のように音響空間インパルス応答推
定装置を構成したので、第１のディジタルフィルタによ
り、前記直接波と前記複数の反射波の伝搬時間とそれら
のレベルが、１系統の信号として設定される。その第１
のディジタルフィルタの出力信号に対してのみ、第２の
ディジタルフィルタの積和演算が行われる。従って、前
記課題を解決できるのである。
【０００８】
【発明の実施の形態】図１は、本発明の実施形態を示す
音響空間インパルス応答推定装置の構成ブロック図であ
る。このインパルス応答推定装置では、スピーカ等の音
源から出力するインパルスを入力する入力端子１０に、
第１のディジタルフィルタ１６が接続されている。ディ
ジタルフィルタ１６の出力側には、従来のディジタルフ
ィルタ１３−１〜１３−ｎと同等の機能を有する第２の
ディジタルフィルタ１３が接続され、ディジタルフィル
タ１３の出力側が、出力端子１５に接続されている。２
つのディジタルフィルタ１３，１６は、例えばＤＳＰ等
で構成され、出力端子１５から出力された信号が受音点
に相当するマイクに入力されることになる。このインパ
ルス応答推定装置の動作を説明する前に、図５を参照し
つつ、ディジタルフィルタ１６の機能を説明する。
音響空間インパルス応答推定装置の構成ブロック図であ
る。このインパルス応答推定装置では、スピーカ等の音
源から出力するインパルスを入力する入力端子１０に、
第１のディジタルフィルタ１６が接続されている。ディ
ジタルフィルタ１６の出力側には、従来のディジタルフ
ィルタ１３−１〜１３−ｎと同等の機能を有する第２の
ディジタルフィルタ１３が接続され、ディジタルフィル
タ１３の出力側が、出力端子１５に接続されている。２
つのディジタルフィルタ１３，１６は、例えばＤＳＰ等
で構成され、出力端子１５から出力された信号が受音点
に相当するマイクに入力されることになる。このインパ
ルス応答推定装置の動作を説明する前に、図５を参照し
つつ、ディジタルフィルタ１６の機能を説明する。
【０００９】図５は、図４の推定装置の構成の変形を示
す図であり、図４及び図１中と同じ要素には共通の符号
が付されている。従来のインパルス応答推定装置と同等
の機能を得るために、図４の要素を変形すると、図５の
構成が考えられる。これによると、入力端子１０には複
数の乗算器１２−１，１２−２，…，１２−ｎが接続さ
れ、該各乗算器１２−１，１２−２，…，１２−ｎの出
力側に、遅延素子１１−１，１１−２，…，１１−ｎ
が、それぞれ接続されている。遅延素子１１−１，１１
−２，…，１１−ｎの出力側が、累算器１４が接続さ
れ、累算器１４にディジタルフィルタ１３が、接続され
ている。乗算器１２−１〜１２−ｎと遅延素子１１−１
〜１１−ｎと累算器１４とは、一つディジタルフィルタ
で構成でき、それが図１のディジタルフィルタ１６であ
る。ディジタルフィルタ１６は、入力されたインパルス
に対して複数の係数を乗算し、各乗算結果をそれぞれ遅
延し、それらを累算する機能を有している。但し、これ
らの処理は、主観受音に影響を与える有効な直接波と反
射波の成分だけとする。このように個数ｎを設定するこ
とで、ディジタルフィルタ１３での演算量を制限する。
一方、ディジタルフィルタ１３は、従来のデジタルフィ
ルタ１３−１〜１３−ｎと同様に、図３の直接波のイン
パルス応答から遅延時間を除いた成分を特性として持
つ。
す図であり、図４及び図１中と同じ要素には共通の符号
が付されている。従来のインパルス応答推定装置と同等
の機能を得るために、図４の要素を変形すると、図５の
構成が考えられる。これによると、入力端子１０には複
数の乗算器１２−１，１２−２，…，１２−ｎが接続さ
れ、該各乗算器１２−１，１２−２，…，１２−ｎの出
力側に、遅延素子１１−１，１１−２，…，１１−ｎ
が、それぞれ接続されている。遅延素子１１−１，１１
−２，…，１１−ｎの出力側が、累算器１４が接続さ
れ、累算器１４にディジタルフィルタ１３が、接続され
ている。乗算器１２−１〜１２−ｎと遅延素子１１−１
〜１１−ｎと累算器１４とは、一つディジタルフィルタ
で構成でき、それが図１のディジタルフィルタ１６であ
る。ディジタルフィルタ１６は、入力されたインパルス
に対して複数の係数を乗算し、各乗算結果をそれぞれ遅
延し、それらを累算する機能を有している。但し、これ
らの処理は、主観受音に影響を与える有効な直接波と反
射波の成分だけとする。このように個数ｎを設定するこ
とで、ディジタルフィルタ１３での演算量を制限する。
一方、ディジタルフィルタ１３は、従来のデジタルフィ
ルタ１３−１〜１３−ｎと同様に、図３の直接波のイン
パルス応答から遅延時間を除いた成分を特性として持
つ。
【００１０】図６は、図１中のディジタルフィルタ１６
の出力信号を示す図であり、この図６を参照しつつ、図
１の音響空間インパルス応答推定装置の動作を説明す
る。入力端子１０から入力されたインパルスに対し、デ
ィジタルフィルタ１６は、直接波と各反射波の伝搬時間
とそれらのレベルを設定する。つまり、ディジタルフィ
ルタ１６はその機能により、伝搬による減衰を考慮した
係数α0 〜αn を掛け、掛けた結果を、対象音響空間で
ある部屋の幾何学的配置から推定された直接波と各反射
波の遅延時間ｔ０，ｔ１，…，ｔｎ分だけずらし、図６
のような１系統の出力信号を生成する。ディジタルフィ
ルタ１３は、ディジタルフィルタ１６の出力信号と所有
する特性の成分とを、積和演算する。これにより、直接
波と反射波のインパルス応答が求められ、そのインパル
ス応答が出力端子１５を介して出力される。
の出力信号を示す図であり、この図６を参照しつつ、図
１の音響空間インパルス応答推定装置の動作を説明す
る。入力端子１０から入力されたインパルスに対し、デ
ィジタルフィルタ１６は、直接波と各反射波の伝搬時間
とそれらのレベルを設定する。つまり、ディジタルフィ
ルタ１６はその機能により、伝搬による減衰を考慮した
係数α0 〜αn を掛け、掛けた結果を、対象音響空間で
ある部屋の幾何学的配置から推定された直接波と各反射
波の遅延時間ｔ０，ｔ１，…，ｔｎ分だけずらし、図６
のような１系統の出力信号を生成する。ディジタルフィ
ルタ１３は、ディジタルフィルタ１６の出力信号と所有
する特性の成分とを、積和演算する。これにより、直接
波と反射波のインパルス応答が求められ、そのインパル
ス応答が出力端子１５を介して出力される。
【００１１】ここで、積和演算の回数を説明する。従来
では、考慮に入れる直接波及び反射波の数ｎと、各ディ
ジタルフィルタ１３−１〜１３−ｎのサンプル数つまり
タップ数（Ｎｔ）とを乗じた回数（ｎ×Ｎｔ）の積和演
算が必要であった。これに対し、図１の装置の積和演算
回数は、ディジタルフィルタ１６のタップ数（ｎ）と、
ディジタルフィルタ１３のタップ数（Ｎｔ）の和（ｎ＋
Ｎｔ）になる。以上のように、本実施形態では、ディジ
タルフィルタ１６と、該ディジタルフィルタ１６の出力
信号に積和演算を行うディジタルフィルタ１３とで、音
響空間インパルス応答推定装置を構成し、図５の機能を
展開するようにしているので、積和演算の回数を、従来
に比べて大幅に減じることができる。なお、本発明は、
ディジタルフィルタ１３，１６を、上述のようにＤＳＰ
をはじめとする各種プロセッサで構成することができる
が、計算機でのソフトウエア処理で実現することも可能
である。
では、考慮に入れる直接波及び反射波の数ｎと、各ディ
ジタルフィルタ１３−１〜１３−ｎのサンプル数つまり
タップ数（Ｎｔ）とを乗じた回数（ｎ×Ｎｔ）の積和演
算が必要であった。これに対し、図１の装置の積和演算
回数は、ディジタルフィルタ１６のタップ数（ｎ）と、
ディジタルフィルタ１３のタップ数（Ｎｔ）の和（ｎ＋
Ｎｔ）になる。以上のように、本実施形態では、ディジ
タルフィルタ１６と、該ディジタルフィルタ１６の出力
信号に積和演算を行うディジタルフィルタ１３とで、音
響空間インパルス応答推定装置を構成し、図５の機能を
展開するようにしているので、積和演算の回数を、従来
に比べて大幅に減じることができる。なお、本発明は、
ディジタルフィルタ１３，１６を、上述のようにＤＳＰ
をはじめとする各種プロセッサで構成することができる
が、計算機でのソフトウエア処理で実現することも可能
である。
【００１２】
【発明の効果】以上詳細に説明したように、第１及び第
２の発明によれば、入力されたインパルスに対し、直接
波と複数の反射波の伝搬時間とそれらのレベルを設定す
る第１のディジタルフィルタと、第１のディジタルフィ
ルタの出力信号と積和演算を行う第２のディジタルフィ
ルタとを備えているので、音響空間インパルス応答推定
装置における積和演算の回数を、従来に比べて大幅に減
じることができる。
２の発明によれば、入力されたインパルスに対し、直接
波と複数の反射波の伝搬時間とそれらのレベルを設定す
る第１のディジタルフィルタと、第１のディジタルフィ
ルタの出力信号と積和演算を行う第２のディジタルフィ
ルタとを備えているので、音響空間インパルス応答推定
装置における積和演算の回数を、従来に比べて大幅に減
じることができる。
【図１】本発明の実施形態を示す音響空間インパルス応
答推定装置の構成ブロック図である。
答推定装置の構成ブロック図である。
【図２】音の反射波と直接波を幾何学的に説明する図で
ある。
ある。
【図３】図２中のマイク３に入力される音波とインパル
ス波形を示す図である。
ス波形を示す図である。
【図４】従来の音響空間インパルス応答推定装置を示す
構成ブロック図である。
構成ブロック図である。
【図５】図４の推定装置の構成の変形を示す図である。
【図６】図１中のディジタルフィルタ１６の出力信号を
示す図である。
示す図である。
１ 部屋 ２ スピーカ ３ マイク １３ 第１のディジタルフィルタ １６ 第２のディジタルフィルタ Ｓｄ 直接波 Ｓ１〜Ｓ６ 反射波 ｔ０〜ｔｎ 遅延時間
Claims (2)
- 【請求項１】 音源から受音点に到達する音波の経路
に、直接波の経路と複数の反射波の経路とが存在する音
響空間におけるインパルス応答を求める音響空間インパ
ルス応答推定装置において、 入力されたインパルスに対し、前記音源から前記受音点
に到達する前記直接波と前記複数の反射波の伝搬時間と
それらのレベルとを推定して設定する第１のディジタル
フィルタと、 前記直接波のみのインパルス応答から前記伝搬時間を除
いた成分を特性として持ち、前記第１のディジタルフィ
ルタの出力信号と該特性との積和演算を行うことで、前
記直接波と前記複数の反射波のインパルス応答を求める
第２のディジタルフィルタとを、 備えたことを特徴とする音響空間インパルス応答推定装
置。 - 【請求項２】 前記第１のディジタルフィルタは、主観
受音として有効な分だけの前記直接波と前記複数の反射
波の伝搬時間とそれらのレベルとを推定して設定し、前
記第２のディジタルフィルタにおける積和演算回数を制
限する構成にしたことを特徴とする請求項１記載の音響
空間インパルス応答推定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP8108828A JPH09292885A (ja) | 1996-04-30 | 1996-04-30 | 音響空間インパルス応答推定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP8108828A JPH09292885A (ja) | 1996-04-30 | 1996-04-30 | 音響空間インパルス応答推定装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH09292885A true JPH09292885A (ja) | 1997-11-11 |
Family
ID=14494583
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP8108828A Withdrawn JPH09292885A (ja) | 1996-04-30 | 1996-04-30 | 音響空間インパルス応答推定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH09292885A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007333545A (ja) * | 2006-06-15 | 2007-12-27 | Clarion Co Ltd | 吸音特性測定方法および吸音特性測定装置 |
| JP2015535962A (ja) * | 2012-10-09 | 2015-12-17 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | オーディオ干渉推定のための方法及び装置 |
-
1996
- 1996-04-30 JP JP8108828A patent/JPH09292885A/ja not_active Withdrawn
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007333545A (ja) * | 2006-06-15 | 2007-12-27 | Clarion Co Ltd | 吸音特性測定方法および吸音特性測定装置 |
| JP2015535962A (ja) * | 2012-10-09 | 2015-12-17 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | オーディオ干渉推定のための方法及び装置 |
| RU2651616C2 (ru) * | 2012-10-09 | 2018-04-23 | Конинклейке Филипс Н.В. | Способ и устройство оценки аудиопомехи |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5371799A (en) | Stereo headphone sound source localization system | |
| US6002776A (en) | Directional acoustic signal processor and method therefor | |
| JP3397269B2 (ja) | 多チャネル反響消去方法 | |
| CN1762178B (zh) | 用于处理电声音信号的方法 | |
| Buchner et al. | Wave-domain adaptive filtering: Acoustic echo cancellation for full-duplex systems based on wave-field synthesis | |
| US8306243B2 (en) | Audio device | |
| EP3026666B1 (en) | Reverberant sound adding apparatus, reverberant sound adding method, and reverberant sound adding program | |
| JP3359460B2 (ja) | 適応フィルタ及びエコーキャンセラ | |
| JP2007129736A (ja) | 音声信号中の暗騒音を抑制する方法及び装置、並びに反響消去を伴う対応する装置 | |
| JP2021527360A5 (ja) | ||
| JPH09261792A (ja) | 受音方法及びその装置 | |
| TWI245258B (en) | Method and related apparatus for generating audio reverberation effect | |
| JPH09292885A (ja) | 音響空間インパルス応答推定装置 | |
| JP7403436B2 (ja) | 異なる音場の複数の録音音響信号を合成する音響信号合成装置、プログラム及び方法 | |
| JP2004228621A (ja) | エコー抑圧装置 | |
| JP4928376B2 (ja) | 収音装置、収音方法、その方法を用いた収音プログラム、および記録媒体 | |
| WO2022048376A1 (zh) | 一种声场扩展方法、电路、电子设备及存储介质 | |
| Menzer | Binaural reverberation using two parallel feedback delay networks | |
| JP3583980B2 (ja) | 収音装置及び受信装置 | |
| JP2009005261A (ja) | 収音装置、収音方法、その方法を用いた収音プログラム、および記録媒体 | |
| JPH0126598B2 (ja) | ||
| JP3579508B2 (ja) | オーディオ装置 | |
| US20080285768A1 (en) | Method and System for Modifying and Audio Signal, and Filter System for Modifying an Electrical Signal | |
| CN110719564A (zh) | 音效处理方法和装置 | |
| JP3941351B2 (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: 20030701 |