JPH09292885A - Acoustic space impulse response estimating device - Google Patents

Acoustic space impulse response estimating device

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
Application number
JP8108828A
Other languages
Japanese (ja)
Inventor
Shinichi Kawada
眞一 川田
Masaharu Shimada
正治 島田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP8108828A priority Critical patent/JPH09292885A/en
Publication of JPH09292885A publication Critical patent/JPH09292885A/en
Withdrawn legal-status Critical Current

Links

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

PROBLEM TO BE SOLVED: To obtain an acoustic space impulse response estimating device with a small number of times of product-sum operation. SOLUTION: A digital filter 16 sets the transport time and level of direct and reflection waves against an impulse input from an input terminal 10. A digital filter 13 has as characteristics components eliminating the transport time from the impulse response of only direct waves, performs product-sum operation between the output of the digital filter 16 and its component, and obtains the direct and wave and plural reflection waves.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【発明の属する技術分野】本発明は、音声会議或いはテ
レビ会議(以下、TV会議という)の音声マンマシーン
インタフェースや、バーチャルリアリティの音声処理等
に使用される音響空間インパルス応答を推定する音響空
間インパルス応答推定装置に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an acoustic space impulse for estimating an acoustic space impulse response used for a voice man-machine interface of a voice conference or a video conference (hereinafter referred to as a TV conference), a voice processing of virtual reality, and the like. The present invention relates to a response estimation device.

【0002】[0002]

【従来の技術】音声会議やTV会議における音響のステ
レオ再生に際し、まるでコンサートホールで音楽を聞い
ているような臨場感を出すために、電気的に反射音や残
響音を操作することが行われている。この操作を実現す
るためには、ホールや部屋での音響空間インパルス応答
を求める必要がある。つまり、音源のスピーカから放射
したインパルスが、マイクにどのような波形で到達する
かを求める必要がある。従来では、DSP(Digital Si
gnal Processor;ディジタル信号処理プロセッサ)等を
用いた音響空間インパルス応答推定装置で、そのインパ
ルス応答を推定していた。
2. Description of the Related Art In stereophonic reproduction of sound in a voice conference or a TV conference, the reflected sound and the reverberant sound are electrically operated in order to give a realistic sensation of listening to music in a concert hall. ing. In order to realize this operation, it is necessary to find the acoustic spatial impulse response in the hall or room. That is, it is necessary to find what waveform the impulse emitted from the speaker of the sound source reaches the microphone. Conventionally, DSP (Digital Si
The impulse response was estimated by an acoustic space impulse response estimation device using a gnal Processor (digital signal processor) or the like.

【0003】[0003]

【発明が解決しようとする課題】しかしながら、従来の
音響空間インパルス応答推定装置では、次のような課題
があった。通常の部屋でのインパルス応答の長さは、高
次の反射波の影響で0.2〜0.6秒に達する。このよ
うな高次の反射を処理しようとすると、演算量が増加す
るという課題があった。以下、その理由を説明する。図
2は、音の反射波と直接波を幾何学的に説明する図であ
り、部屋における直接波と反射波を示している。部屋1
内に配置されたスピーカ2から放射した音波は、種々の
経路kd,k1,k2,k3,k4,k5,k6を通過
してマイク3に達する。経路kdは、スピーカ2から直
接マイク3に入力される直接波の経路である。各経路k
1,k2を伝搬する音波は、部屋1の壁で2度反射して
マイク3にそれぞれ到達する。各経路k3,k4を伝搬
する音波は、1度反射してマイク3に到達する。各経路
k5,k6を伝搬する音波も、壁で2度反射してマイク
3に到達する。つまり、各経路k1〜k6を伝搬してマ
イク3に達した音波は、反射波である。経路k3,k4
を伝搬してマイク3に達した音波は、一次反射波であ
り、経路k1,k2,k5,k6を伝搬してマイク3に
達した音波は二次反射波である。図2中のz1,z2,
z3,z4,z5,z6は、音波をマイク3側から見た
スピーカ2の虚像である。つまり、マイク3には、虚像
z1〜z6の位置にスピーカ2が配置してあるように、
音波が入力されている。
However, the conventional acoustic space impulse response estimation device has the following problems. The length of the impulse response in a normal room reaches 0.2 to 0.6 seconds due to the influence of higher-order reflected waves. There is a problem in that the amount of calculation increases when processing such high-order reflections. Hereinafter, the reason will be described. FIG. 2 is a diagram for geometrically explaining the reflected wave and the direct wave of the sound, and shows the direct wave and the reflected wave in the room. Room 1
The sound waves emitted from the speaker 2 arranged inside reach the microphone 3 through various paths kd, k1, k2, k3, k4, k5 and k6. The route kd is a route of a direct wave directly input from the speaker 2 to the microphone 3. Each route k
The sound waves propagating through 1 and k2 are reflected twice by the wall of the room 1 and reach the microphones 3, respectively. The sound waves propagating through the paths k3 and k4 are reflected once and reach the microphone 3. The sound waves propagating through the paths k5 and k6 also reach the microphone 3 after being reflected twice by the wall. That is, the sound wave that propagates through each of the paths k1 to k6 and reaches the microphone 3 is a reflected wave. Route k3, k4
The sound wave that propagates through and reaches the microphone 3 is a primary reflected wave, and the sound wave that propagates through the paths k1, k2, k5, k6 and reaches the microphone 3 is a secondary reflected wave. Z1, z2 in FIG.
z3, z4, z5, and z6 are virtual images of the speaker 2 when sound waves are viewed from the microphone 3 side. That is, in the microphone 3, the speaker 2 is arranged at the positions of the virtual images z1 to z6.
Sound waves are being input.

【0004】図3は、図2中のマイク3に入力される音
波とインパルス波形を示す図である。この図3では、経
路kdを伝搬した直接波Sdと、経路k3,k4を伝搬
した一次反射波S3,S4と、経路k1,k2を伝搬し
た二次反射波S1,S2と、経路k5,k6を伝搬した
二次反射波S5,S6とが、示されている。なお、簡単
のため、経路k1と経路k2、経路k3と経路k4、及
び経路k5と経路k6は対称であり、音波の伝搬距離は
等しいものとしている。直接波Sdは、時間t0 後にマ
イク3に到達する。各一次反射波S3,S4は、時間t
1後にマイク3にそれぞれ到達する。各二次反射波S
1,S2は、時間t2後にマイク3にそれぞれ到達し、
各二次反射波S5,S6は、時間t3後にマイク3にそ
れぞれ到達する。これらの反射波S1,S2と、反射波
S3,S4と反射波S5,S6は、スピーカ2からマイ
ク3までの経路上の距離が異なるので、該マイク3への
到着時間が異なる。また、これら反射波S1〜S6は、
壁で反射する度に、音のエネルギの一部が吸収されて減
衰する。これらの直接波Sd、反射波S1〜S6が合成
されて、インパルス波形が得られる。これがインパルス
応答である。図2及び図3では二次反射までしか示して
いないが、実際にはもっと高次の反射まで存在するの
で、インパルス波形の継続時間はもっと長時間になり、
通常の部屋では0.2〜0.6秒になる。
FIG. 3 is a diagram showing a sound wave and an impulse waveform input to the microphone 3 in FIG. In FIG. 3, the direct wave Sd propagating on the route kd, the primary reflected waves S3 and S4 propagating on the routes k3 and k4, the secondary reflected waves S1 and S2 propagating on the routes k1 and k2, and the routes k5 and k6. The secondary reflected waves S5 and S6 propagated through are shown. Note that, for simplification, the route k1 and the route k2, the route k3 and the route k4, and the route k5 and the route k6 are symmetrical, and the sound wave propagation distances are equal. The direct wave Sd reaches the microphone 3 after the time t0. Each of the primary reflected waves S3 and S4 has time t
The microphones 3 arrive after 1 each. Each secondary reflected wave S
1, S2 reach the microphone 3 after time t2,
The secondary reflected waves S5 and S6 reach the microphone 3 after time t3. The reflected waves S1 and S2, the reflected waves S3 and S4, and the reflected waves S5 and S6 have different distances on the path from the speaker 2 to the microphone 3, and therefore have different arrival times at the microphone 3. Further, these reflected waves S1 to S6 are
Each time it is reflected by the wall, part of the sound energy is absorbed and attenuated. The direct wave Sd and the reflected waves S1 to S6 are combined to obtain an impulse waveform. This is the impulse response. Although only the secondary reflection is shown in FIGS. 2 and 3, since the higher-order reflection is actually present, the duration of the impulse waveform becomes longer,
It takes 0.2 to 0.6 seconds in a normal room.

【0005】図4は、従来の音響空間インパルス応答推
定装置を示す構成ブロック図である。入力端子10は、
スピーカの出力するインパルスを入力する端子であり、
その入力端子10には、複数の遅延素子11−1,11
−2,…,11−nが接続されている。各遅延素子11
−1〜11−nは、直接波及び反射波の伝搬距離に対応
する遅延時間を設定する機能をそれぞれ有し、これら遅
延素子11−1〜11−nの数nは、直接波と反射波
(高次の反射波も含む)の数である。設定される遅延時
間は、図2の時間t0,t1,t2,t3,…に相当す
る。各遅延素子11−1,11−2,…11−nの出力
側が、直接波及び各反射波のレベルを設定する乗算器1
2−1,12−2,…,12−nに、それぞれ接続され
ている。各乗算器12−1,12−2,…,12−nの
出力側が、ディジタルフィルタ13−1,13−2,
…,13−nにそれぞれ接続されている。ディジタルフ
ィルタ13−1,13−2,…,13−nの出力側が累
算器14に接続され、該累算器14の出力側が出力端子
15に接続されている。出力端子15からの出力信号が
マイクに入力される構成である。各ディジタルフィルタ
13−1〜13−nの有する特性は同じであり、マイク
3へ到達する直接波の特性のうち、遅延時間を除いたも
のである。各ディジタルフィルタ13−1〜13−nで
は、入力された信号に対して、その特性で積和演算を行
う機能をそれぞれ有している。
FIG. 4 is a block diagram showing the configuration of a conventional acoustic spatial impulse response estimation device. The input terminal 10 is
It is a terminal to input the impulse output from the speaker,
The input terminal 10 has a plurality of delay elements 11-1, 11
-2, ..., 11-n are connected. Each delay element 11
-1 to 11-n each have a function of setting a delay time corresponding to the propagation distance of the direct wave and the reflected wave, and the number n of these delay elements 11-1 to 11-n is the direct wave and the reflected wave. (Including higher-order reflected waves). The set delay time corresponds to the times t0, t1, t2, t3, ... Of FIG. The output side of each delay element 11-1, 11-2, ... 11-n sets a multiplier 1 for setting the level of the direct wave and each reflected wave.
2-1, 12-2, ..., 12-n are respectively connected. The output sides of the multipliers 12-1, 12-2, ..., 12-n are digital filters 13-1, 13-2,
..., 13-n, respectively. The output sides of the digital filters 13-1, 13-2, ..., 13-n are connected to the accumulator 14, and the output side of the accumulator 14 is connected to the output terminal 15. The output signal from the output terminal 15 is input to the microphone. The characteristics of the digital filters 13-1 to 13-n are the same, and the delay time is excluded from the characteristics of the direct wave reaching the microphone 3. Each of the digital filters 13-1 to 13-n has a function of performing a product-sum operation on the input signal with its characteristic.

【0006】各反射波は、直接波に類似な波形を遅延さ
せて減衰させた波形として表すことができる。反射波の
遅延時間は、反射経路を音速で割ったものであり、部屋
の幾何学的配置から推定できる。各遅延素子11−1〜
11−nは、入力端子10から入力されたインパルスに
対して、その推定された遅延時間分だけ、それぞれ遅延
させる。各乗算器12−1〜12−nは、各遅延素子1
1−1〜11−nの出力信号に対して、係数α0 〜αn
を掛けて直接波と各反射波のレベルをそれぞれ設定す
る。各ディジタルフィルタ13−1〜13−nは、前記
特性と乗算器12−1〜12−nの出力信号の積和を求
め、直接波と各反射波の波形を形成する。累算器14の
累算によって、合成したインパルス応答のインパルス波
形が得られる。直接波の遅延時間を除いた信号の継続時
間は、10〜20ms程度である。例えば、音声を16
kHzでサンプリングする場合、約160〜320サン
プルの積和演算が必要である。これ自身はそれほどの演
算量ではないが、反射波も含めると、n倍の演算量が必
要である。例えば、図2の場合には、約640〜128
0サンプルの積和演算が必要になる。
Each reflected wave can be represented as a delayed and attenuated waveform similar to the direct wave. The delay time of the reflected wave is the reflection path divided by the speed of sound and can be estimated from the geometrical arrangement of the room. Each delay element 11-1 to 11-1
11-n delays the impulse input from the input terminal 10 by the estimated delay time. Each of the multipliers 12-1 to 12-n includes each delay element 1
Coefficients α 0 to α n for the output signals of 1-1 to 11- n
To set the levels of the direct wave and each reflected wave. Each of the digital filters 13-1 to 13-n obtains the product sum of the above characteristics and the output signals of the multipliers 12-1 to 12-n and forms the waveforms of the direct wave and the reflected waves. By the accumulation of the accumulator 14, the impulse waveform of the combined impulse response is obtained. The signal duration excluding the delay time of the direct wave is about 10 to 20 ms. For example, sound 16
When sampling at kHz, a product-sum operation of about 160 to 320 samples is required. Although this is not so much calculation amount itself, when the reflected wave is also included, the calculation amount of n times is required. For example, in the case of FIG. 2, about 640 to 128.
A product-sum operation of 0 samples is required.

【0007】[0007]

【課題を解決するための手段】第1の発明は、前記課題
を解決するために、音源から受音点に到達する音波の経
路に、直接波の経路と複数の反射波の経路とが存在する
音響空間におけるインパルス応答を求める音響空間イン
パルス応答推定装置において、次のような構成にしてい
る。即ち、第1の発明の音響空間インパルス応答推定装
置は、入力されたインパルスに対し、前記音源から前記
受音点に到達する前記直接波と前記複数の反射波の伝搬
時間とそれらのレベルとを推定して設定する第1のディ
ジタルフィルタと、前記直接波のみのインパルス応答か
ら前記伝搬時間を除いた成分を特性として持ち、前記第
1のディジタルフィルタの出力信号と該特性との積和演
算を行うことで、前記直接波と前記複数の反射波のイン
パルス応答を求める第2のディジタルフィルタとを備え
ている。第2の発明は、第1の発明における第1のディ
ジタルフィルタは、主観受音として有効な分だけの前記
直接波と前記複数の反射波の伝搬時間とそれらのレベル
を推定して設定し、前記第2のディジタルフィルタにお
ける積和演算回数を制限構成にしている。第1,第2の
発明によれば、以上のように音響空間インパルス応答推
定装置を構成したので、第1のディジタルフィルタによ
り、前記直接波と前記複数の反射波の伝搬時間とそれら
のレベルが、1系統の信号として設定される。その第1
のディジタルフィルタの出力信号に対してのみ、第2の
ディジタルフィルタの積和演算が行われる。従って、前
記課題を解決できるのである。
According to a first aspect of the present invention, in order to solve the above-mentioned problems, there are a direct wave path and a plurality of reflected wave paths in a sound wave path from a sound source to a sound receiving point. In the acoustic space impulse response estimation device that obtains the impulse response in the acoustic space, the following configuration is used. That is, the acoustic spatial impulse response estimation apparatus of the first invention determines the propagation times of the direct wave and the plurality of reflected waves that reach the sound receiving point from the sound source and their levels with respect to the input impulse. It has a first digital filter that is estimated and set, and a component obtained by removing the propagation time from the impulse response of only the direct wave as a characteristic, and a product-sum operation of the output signal of the first digital filter and the characteristic is performed. And a second digital filter that obtains impulse responses of the direct wave and the plurality of reflected waves by performing. According to a second invention, the first digital filter in the first invention estimates and sets the propagation times and the levels of the direct wave and the plurality of reflected waves that are effective as subjective sound reception, The number of product-sum operations in the second digital filter is limited. According to the first and second aspects of the invention, since the acoustic spatial impulse response estimating apparatus is configured as described above, the propagation time of the direct wave and the plurality of reflected waves and their levels are determined by the first digital filter. It is set as a single-system signal. The first
The product-sum operation of the second digital filter is performed only on the output signal of the second digital filter. Therefore, the above problem can be solved.

【0008】[0008]

【発明の実施の形態】図1は、本発明の実施形態を示す
音響空間インパルス応答推定装置の構成ブロック図であ
る。このインパルス応答推定装置では、スピーカ等の音
源から出力するインパルスを入力する入力端子10に、
第1のディジタルフィルタ16が接続されている。ディ
ジタルフィルタ16の出力側には、従来のディジタルフ
ィルタ13−1〜13−nと同等の機能を有する第2の
ディジタルフィルタ13が接続され、ディジタルフィル
タ13の出力側が、出力端子15に接続されている。2
つのディジタルフィルタ13,16は、例えばDSP等
で構成され、出力端子15から出力された信号が受音点
に相当するマイクに入力されることになる。このインパ
ルス応答推定装置の動作を説明する前に、図5を参照し
つつ、ディジタルフィルタ16の機能を説明する。
1 is a block diagram showing the configuration of an acoustic spatial impulse response estimation apparatus according to an embodiment of the present invention. In this impulse response estimation device, an input terminal 10 for inputting an impulse output from a sound source such as a speaker,
The first digital filter 16 is connected. A second digital filter 13 having a function equivalent to that of the conventional digital filters 13-1 to 13-n is connected to the output side of the digital filter 16, and the output side of the digital filter 13 is connected to the output terminal 15. There is. Two
Each of the digital filters 13 and 16 is composed of, for example, a DSP, and the signal output from the output terminal 15 is input to the microphone corresponding to the sound receiving point. Before explaining the operation of the impulse response estimation device, the function of the digital filter 16 will be described with reference to FIG.

【0009】図5は、図4の推定装置の構成の変形を示
す図であり、図4及び図1中と同じ要素には共通の符号
が付されている。従来のインパルス応答推定装置と同等
の機能を得るために、図4の要素を変形すると、図5の
構成が考えられる。これによると、入力端子10には複
数の乗算器12−1,12−2,…,12−nが接続さ
れ、該各乗算器12−1,12−2,…,12−nの出
力側に、遅延素子11−1,11−2,…,11−n
が、それぞれ接続されている。遅延素子11−1,11
−2,…,11−nの出力側が、累算器14が接続さ
れ、累算器14にディジタルフィルタ13が、接続され
ている。乗算器12−1〜12−nと遅延素子11−1
〜11−nと累算器14とは、一つディジタルフィルタ
で構成でき、それが図1のディジタルフィルタ16であ
る。ディジタルフィルタ16は、入力されたインパルス
に対して複数の係数を乗算し、各乗算結果をそれぞれ遅
延し、それらを累算する機能を有している。但し、これ
らの処理は、主観受音に影響を与える有効な直接波と反
射波の成分だけとする。このように個数nを設定するこ
とで、ディジタルフィルタ13での演算量を制限する。
一方、ディジタルフィルタ13は、従来のデジタルフィ
ルタ13−1〜13−nと同様に、図3の直接波のイン
パルス応答から遅延時間を除いた成分を特性として持
つ。
FIG. 5 is a diagram showing a modification of the configuration of the estimation apparatus of FIG. 4, and the same elements as those in FIGS. 4 and 1 are designated by the same reference numerals. If the elements of FIG. 4 are modified in order to obtain a function equivalent to that of the conventional impulse response estimation device, the configuration of FIG. 5 can be considered. According to this, a plurality of multipliers 12-1, 12-2, ..., 12-n are connected to the input terminal 10, and the output side of each of the multipliers 12-1, 12-2 ,. , Delay elements 11-1, 11-2, ..., 11-n
Are connected respectively. Delay elements 11-1, 11
The output side of −2, ..., 11-n is connected to the accumulator 14, and the digital filter 13 is connected to the accumulator 14. Multipliers 12-1 to 12-n and delay element 11-1
.About.11-n and the accumulator 14 can be constituted by one digital filter, which is the digital filter 16 in FIG. The digital filter 16 has a function of multiplying the input impulse by a plurality of coefficients, delaying each multiplication result, and accumulating them. However, these processes are limited to effective direct wave and reflected wave components that affect subjective sound reception. By setting the number n in this way, the amount of calculation in the digital filter 13 is limited.
On the other hand, the digital filter 13 has, as the characteristics of the conventional digital filters 13-1 to 13-n, a component obtained by removing the delay time from the impulse response of the direct wave in FIG. 3 as a characteristic.

【0010】図6は、図1中のディジタルフィルタ16
の出力信号を示す図であり、この図6を参照しつつ、図
1の音響空間インパルス応答推定装置の動作を説明す
る。入力端子10から入力されたインパルスに対し、デ
ィジタルフィルタ16は、直接波と各反射波の伝搬時間
とそれらのレベルを設定する。つまり、ディジタルフィ
ルタ16はその機能により、伝搬による減衰を考慮した
係数α0 〜αn を掛け、掛けた結果を、対象音響空間で
ある部屋の幾何学的配置から推定された直接波と各反射
波の遅延時間t0,t1,…,tn分だけずらし、図6
のような1系統の出力信号を生成する。ディジタルフィ
ルタ13は、ディジタルフィルタ16の出力信号と所有
する特性の成分とを、積和演算する。これにより、直接
波と反射波のインパルス応答が求められ、そのインパル
ス応答が出力端子15を介して出力される。
FIG. 6 shows the digital filter 16 shown in FIG.
FIG. 7 is a diagram showing the output signal of FIG. 1, and the operation of the acoustic spatial impulse response estimation device of FIG. 1 will be described with reference to FIG. 6. The digital filter 16 sets the propagation time of the direct wave and each reflected wave and their levels with respect to the impulse input from the input terminal 10. That is, by its function, the digital filter 16 multiplies the coefficients α 0 to α n in consideration of attenuation due to propagation, and multiplies the result by the direct wave and each reflection estimated from the geometrical arrangement of the room as the target acoustic space. The wave delay time t0, t1, ...
One-system output signal such as The digital filter 13 multiplies and sums the output signal of the digital filter 16 and the characteristic component possessed by the digital filter 13. As a result, the impulse response of the direct wave and the reflected wave is obtained, and the impulse response is output via the output terminal 15.

【0011】ここで、積和演算の回数を説明する。従来
では、考慮に入れる直接波及び反射波の数nと、各ディ
ジタルフィルタ13−1〜13−nのサンプル数つまり
タップ数(Nt)とを乗じた回数(n×Nt)の積和演
算が必要であった。これに対し、図1の装置の積和演算
回数は、ディジタルフィルタ16のタップ数(n)と、
ディジタルフィルタ13のタップ数(Nt)の和(n+
Nt)になる。以上のように、本実施形態では、ディジ
タルフィルタ16と、該ディジタルフィルタ16の出力
信号に積和演算を行うディジタルフィルタ13とで、音
響空間インパルス応答推定装置を構成し、図5の機能を
展開するようにしているので、積和演算の回数を、従来
に比べて大幅に減じることができる。なお、本発明は、
ディジタルフィルタ13,16を、上述のようにDSP
をはじめとする各種プロセッサで構成することができる
が、計算機でのソフトウエア処理で実現することも可能
である。
Here, the number of product-sum operations will be described. Conventionally, the product-sum operation of the number of times (n × Nt) obtained by multiplying the number n of direct waves and reflected waves to be taken into consideration by the number of samples of each digital filter 13-1 to 13-n, that is, the number of taps (Nt) is performed. Was needed. On the other hand, the number of product-sum operations of the device of FIG.
Sum of tap numbers (Nt) of the digital filter 13 (n +
Nt). As described above, in the present embodiment, the digital filter 16 and the digital filter 13 that performs the sum of products operation on the output signals of the digital filter 16 constitute an acoustic spatial impulse response estimation device and expand the function of FIG. Since this is done, the number of product-sum operations can be greatly reduced compared to the conventional case. In addition, the present invention
The digital filters 13 and 16 are connected to the DSP as described above.
It can be configured by various processors such as, but can also be realized by software processing in a computer.

【0012】[0012]

【発明の効果】以上詳細に説明したように、第1及び第
2の発明によれば、入力されたインパルスに対し、直接
波と複数の反射波の伝搬時間とそれらのレベルを設定す
る第1のディジタルフィルタと、第1のディジタルフィ
ルタの出力信号と積和演算を行う第2のディジタルフィ
ルタとを備えているので、音響空間インパルス応答推定
装置における積和演算の回数を、従来に比べて大幅に減
じることができる。
As described in detail above, according to the first and second aspects of the present invention, the propagation times of the direct wave and the plurality of reflected waves and their levels are set for the input impulse. , And the second digital filter for performing the product sum calculation with the output signal of the first digital filter, the number of product sum calculations in the acoustic space impulse response estimation device is significantly larger than that of the conventional one. Can be reduced to

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明の実施形態を示す音響空間インパルス応
答推定装置の構成ブロック図である。
FIG. 1 is a configuration block diagram of an acoustic spatial impulse response estimation device showing an embodiment of the present invention.

【図2】音の反射波と直接波を幾何学的に説明する図で
ある。
FIG. 2 is a diagram for geometrically explaining a reflected wave and a direct wave of sound.

【図3】図2中のマイク3に入力される音波とインパル
ス波形を示す図である。
FIG. 3 is a diagram showing a sound wave and an impulse waveform input to a microphone 3 in FIG.

【図4】従来の音響空間インパルス応答推定装置を示す
構成ブロック図である。
FIG. 4 is a configuration block diagram showing a conventional acoustic space impulse response estimation device.

【図5】図4の推定装置の構成の変形を示す図である。5 is a diagram showing a modification of the configuration of the estimation device in FIG.

【図6】図1中のディジタルフィルタ16の出力信号を
示す図である。
FIG. 6 is a diagram showing an output signal of a digital filter 16 in FIG.

【符号の説明】[Explanation of symbols]

1 部屋 2 スピーカ 3 マイク 13 第1のディジタルフィルタ 16 第2のディジタルフィルタ Sd 直接波 S1〜S6 反射波 t0〜tn 遅延時間 1 room 2 speaker 3 microphone 13 1st digital filter 16 2nd digital filter Sd direct wave S1-S6 reflected wave t0-tn delay time

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 音源から受音点に到達する音波の経路
に、直接波の経路と複数の反射波の経路とが存在する音
響空間におけるインパルス応答を求める音響空間インパ
ルス応答推定装置において、 入力されたインパルスに対し、前記音源から前記受音点
に到達する前記直接波と前記複数の反射波の伝搬時間と
それらのレベルとを推定して設定する第1のディジタル
フィルタと、 前記直接波のみのインパルス応答から前記伝搬時間を除
いた成分を特性として持ち、前記第1のディジタルフィ
ルタの出力信号と該特性との積和演算を行うことで、前
記直接波と前記複数の反射波のインパルス応答を求める
第2のディジタルフィルタとを、 備えたことを特徴とする音響空間インパルス応答推定装
置。
1. An acoustic space impulse response estimation device for obtaining an impulse response in an acoustic space in which a direct wave path and a plurality of reflected wave paths exist in a sound wave path from a sound source to a sound receiving point. A first digital filter that estimates and sets the propagation times of the direct wave and the plurality of reflected waves that reach the sound receiving point from the sound source and their levels with respect to the impulse; The impulse response of the direct wave and the plurality of reflected waves is obtained by performing a product-sum operation of the output signal of the first digital filter and the characteristic, which has a characteristic that is a component excluding the propagation time from the impulse response. An acoustic space impulse response estimation apparatus comprising: a second digital filter to be obtained.
【請求項2】 前記第1のディジタルフィルタは、主観
受音として有効な分だけの前記直接波と前記複数の反射
波の伝搬時間とそれらのレベルとを推定して設定し、前
記第2のディジタルフィルタにおける積和演算回数を制
限する構成にしたことを特徴とする請求項1記載の音響
空間インパルス応答推定装置。
2. The first digital filter estimates and sets the propagation times and levels of the direct wave and the plurality of reflected waves, which are effective as subjective sound reception, and sets the second digital filter. The acoustic spatial impulse response estimating apparatus according to claim 1, wherein the number of multiplication and addition operations in the digital filter is limited.
JP8108828A 1996-04-30 1996-04-30 Acoustic space impulse response estimating device Withdrawn JPH09292885A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP8108828A JPH09292885A (en) 1996-04-30 1996-04-30 Acoustic space impulse response estimating device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP8108828A JPH09292885A (en) 1996-04-30 1996-04-30 Acoustic space impulse response estimating device

Publications (1)

Publication Number Publication Date
JPH09292885A true JPH09292885A (en) 1997-11-11

Family

ID=14494583

Family Applications (1)

Application Number Title Priority Date Filing Date
JP8108828A Withdrawn JPH09292885A (en) 1996-04-30 1996-04-30 Acoustic space impulse response estimating device

Country Status (1)

Country Link
JP (1) JPH09292885A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007333545A (en) * 2006-06-15 2007-12-27 Clarion Co Ltd Sound absorption characteristic measuring method and sound absorption characteristic measuring apparatus
JP2015535962A (en) * 2012-10-09 2015-12-17 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method and apparatus for audio interference estimation

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007333545A (en) * 2006-06-15 2007-12-27 Clarion Co Ltd Sound absorption characteristic measuring method and sound absorption characteristic measuring apparatus
JP2015535962A (en) * 2012-10-09 2015-12-17 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method and apparatus for audio interference estimation
RU2651616C2 (en) * 2012-10-09 2018-04-23 Конинклейке Филипс Н.В. Method and apparatus for audio interference estimation

Similar Documents

Publication Publication Date Title
US6002776A (en) Directional acoustic signal processor and method therefor
JP3397269B2 (en) Multi-channel echo cancellation method
Buchner et al. Wave-domain adaptive filtering: Acoustic echo cancellation for full-duplex systems based on wave-field synthesis
US8306243B2 (en) Audio device
JP2007129736A (en) Method and device for suppressing background noise in voice signal, and corresponding device with echo cancellation
JP3359460B2 (en) Adaptive filter and echo canceller
EP3026666B1 (en) Reverberant sound adding apparatus, reverberant sound adding method, and reverberant sound adding program
JPH0157880B2 (en)
US20080285768A1 (en) Method and System for Modifying and Audio Signal, and Filter System for Modifying an Electrical Signal
TWI245258B (en) Method and related apparatus for generating audio reverberation effect
JPH09292885A (en) Acoustic space impulse response estimating device
JPH09261792A (en) Sound receiving method and its device
JP2003309493A (en) Method, device and program for reducing echo
JP4928376B2 (en) Sound collection device, sound collection method, sound collection program using the method, and recording medium
WO2022048376A1 (en) Sound field expansion method and circuit, electronic device, and storage medium
JP3583980B2 (en) Sound collecting device and receiving device
JP2004228621A (en) Echo suppression apparatus
JPH0126598B2 (en)
US9641953B2 (en) Sound spatialization with room effect, optimized in terms of complexity
Adapa et al. Performance analysis of different adapative algorithms based on acoustic echo cancellation
JPH0922293A (en) Audio device
JPH0965497A (en) Sound image localization device
JP3941351B2 (en) Sound field processing apparatus and sound field processing method
JP3285070B2 (en) Sound receiving method and device
JP2002261659A (en) Multi-channel echo cancellation method, its apparatus, its program, and its storage medium

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