JP5644253B2 - Sound field support apparatus and program - Google Patents
Sound field support apparatus and program Download PDFInfo
- Publication number
- JP5644253B2 JP5644253B2 JP2010183078A JP2010183078A JP5644253B2 JP 5644253 B2 JP5644253 B2 JP 5644253B2 JP 2010183078 A JP2010183078 A JP 2010183078A JP 2010183078 A JP2010183078 A JP 2010183078A JP 5644253 B2 JP5644253 B2 JP 5644253B2
- Authority
- JP
- Japan
- Prior art keywords
- impulse response
- sound
- upsideir
- filter coefficient
- convolution operation
- 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
- 230000005236 sound signal Effects 0.000 claims description 56
- 230000002238 attenuated effect Effects 0.000 claims description 4
- 238000000034 method Methods 0.000 description 20
- 238000012545 processing Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000005192 partition Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000003014 reinforcing effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Landscapes
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Description
本発明は、音響空間内における音響効果を制御する技術に関する。 The present invention relates to a technique for controlling an acoustic effect in an acoustic space.
クラシック音楽等の演奏を催す専用のホールには、ステージの前方に複数列の客席を平行に並べて配したジュークボックス型ホールや、ステージを取り囲むようにして客席のブロックを設けたワインヤード型ホール等の種類がある。この種の音響空間のステージ上で演奏が行われている間に客席にいる者が感じる聴感的な印象は、空間内におけるステージと客席との位置関係、天井の高さ、壁間の距離や向き、天井や壁の材質などといった諸条件に依存して決まる(例えば、特許文献1を参照)。 Special halls for performances such as classical music include a jukebox type hall with multiple rows of seats arranged in parallel in front of the stage, and a wineyard type hall with a block of seats surrounding the stage. There are different types. The auditory impression that a person in the audience feels while performing on a stage in this type of acoustic space is the positional relationship between the stage and the audience in the space, the height of the ceiling, the distance between the walls, It is determined depending on various conditions such as the orientation, the material of the ceiling and the wall (for example, see Patent Document 1).
ところで、この種の音響空間の中には、音響空間における後壁の上方から両側壁に沿ってステージ側に突出したバルコニーを有するものがある。このような設計になっている音響空間では、音響空間の天井で反射した音がバルコニーの下の客席に至る途中でバルコニーに遮られてしまう。このため、バルコニーの下の客席の受聴者が感じる残響感や拡がり感などの聴感的な印象である空間感が、上方にバルコニーがない客席の受聴者が感じる空間感に比べて貧弱で違和感のあるものになってしまうという問題があった。
本発明は、このような背景の下に案出されたものであり、バルコニーが設けられた音響空間におけるバルコニーの下にいる者が感じる空間感を改善することを目的とする。
By the way, in this kind of acoustic space, there is one having a balcony protruding to the stage side along both side walls from above the rear wall in the acoustic space. In the acoustic space having such a design, the sound reflected from the ceiling of the acoustic space is blocked by the balcony on the way to the audience seat under the balcony. For this reason, the sense of space, such as the feeling of reverberation and spread, felt by the listeners in the audience seats under the balcony, is weaker and uncomfortable than the sense of space felt by the listeners in the seats without the balcony on the upper side. There was a problem of becoming something.
The present invention has been devised under such a background, and an object thereof is to improve a feeling of space felt by a person under a balcony in an acoustic space provided with the balcony.
本発明は、音響空間内の音源の近傍の収音手段により収音された音の音信号とフィルタ係数列とを畳込む畳込み演算を行い、この畳込み演算の演算結果を前記音源から前記音響空間内の受聴エリアに至る反射音の伝搬経路上にある反射音遮蔽体から前記受聴エリアに音として出力させるディジタルフィルタ手段と、前記音源と前記反射音遮蔽体との間の区間のインパルス応答である第1のインパルス応答と、前記音源と前記収音手段との間の区間のインパルス応答である第2のインパルス応答とを求めるインパルス応答算出手段と、前記ディジタルフィルタ手段の畳込み演算に用いるフィルタ係数列を決定する手段であって、当該フィルタ係数列と前記第2のインパルス応答との畳込み演算結果であるインパルス応答の実効値波形の傾きが前記第1のインパルス応答の実効値波形の傾きに近似するフィルタ係数列を求め、このフィルタ係数列を前記ディジタルフィルタ手段の畳込み演算に用いるフィルタ係数列とするフィルタ係数列決定手段とを具備する音場支援装置を提供する。 The present invention performs a convolution operation that convolves the sound signal of the sound collected by the sound collection means near the sound source in the acoustic space and the filter coefficient sequence, and the operation result of the convolution operation is calculated from the sound source. Digital filter means for outputting a sound from the reflected sound shield on the propagation path of the reflected sound to the listening area in the acoustic space to the listening area, and an impulse response in a section between the sound source and the reflected sound shield And an impulse response calculation means for obtaining a second impulse response which is an impulse response in a section between the sound source and the sound collection means, and is used for a convolution operation of the digital filter means. Means for determining a filter coefficient sequence, wherein the slope of the RMS waveform of the impulse response, which is the result of convolution of the filter coefficient sequence and the second impulse response, is A sound comprising: a filter coefficient string determining unit that obtains a filter coefficient string that approximates the slope of the effective value waveform of the first impulse response, and uses the filter coefficient string as a filter coefficient string used for a convolution operation of the digital filter means. A field support device is provided.
この発明によると、バルコニーのような反射音遮蔽体を有する音響空間における反射音遮蔽体の背後にいる者に対し、反射音遮蔽体のない音響空間において感じる空間感と同じ空間感を感じさせることができる。 According to the present invention, a person behind a reflection sound shield in an acoustic space having a reflection sound shield such as a balcony is caused to feel the same spatial feeling as that felt in an acoustic space without the reflection sound shield. Can do.
なお、特許文献2には、音響空間におけるバルコニーの背後の客席の残響感を改善するための技術の開示がある。しかし、この技術は、ステージからバルコニーの背後の客席に直接伝搬されるはずの直接音を生成してその客席に放音するものであり、本発明とは内容を異にするものである。
Note that
以下、図面を参照しつつ本発明の実施形態について説明する。図1は、本発明の一実施形態である音場支援装置10の構成を示す図である。図2(A)は、音場支援装置10が使用される音響空間80を同空間80の天井82側から見た平面図である。図2(B)は、図2(A)のA−A線断面図である。図2(A)及び図2(B)に示す音響空間80において、仮にバルコニー86がないとすると、ステージ81上の音源90(例えば、ステージ81で奏者により演奏される楽器)から発音された音は、音響空間80の内郭面(天井82や天井82に交差する壁83F,83B,83L,83R)において反射され、各反射面からの反射音がバルコニー86に遮られることなく1階席85の後方の受聴エリア91に到達する。以下、この反射音の伝搬経路を反射音伝搬経路REFという。ここで、音響空間80には、反射音伝搬経路REFの途中に反射音遮蔽体たるバルコニー86がある。このため、何ら策を講じないとすると、音響空間80の上記内郭面で反射された反射音の一部が反射音伝搬経路REFを介して伝搬される過程においてバルコニー86によって遮蔽され、受聴エリア91に到達しなくなる。よって、受聴エリア91の客席にいる者が感じる空間感は、上方にバルコニー86がない客席にいる者が感じる空間感に比べて貧弱で違和感のあるものになってしまう。本実施形態の目的はこの問題を解決することにある。
Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a diagram illustrating a configuration of a sound
図2(A)及び図2(B)に示す音響空間80では、天井82からマイクロホン1L−1,1L−2,1R−1,及び1R−2が吊り下げられている。マイクロホン1L−1,1L−2,1R−1,及び1R−2は、ステージ81上の音源90の近傍に各々位置して、音源90の直接音を独立して収音可能なように配置されている。具体的には、少なくとも臨界距離(直接音と間接音との音響エネルギー比が等しくなる距離)よりも音源90に近い距離であればよく、音源90直近で間接音の成分が少ないほどなおよい。マイクロホン1L−1,1L−2,1R−1,及び1R−2は、音源90の方向から到来する音波に対して高い受音感度を持った単一指向性マイクロホンである。バルコニー86における受聴エリア91側の面78にはスピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)が設けられている。スピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)の各々の放音面は受聴エリア91に向けられている。
In the
図1に示す音場支援装置10は、マイクロホン1L−1,1L−2,1R−1,及び1R−2からステージ81上の音源90が発生した楽音の音信号Sを取得し、この音信号Sから受聴エリア91にいる者が聴感的に感じる空間感を補強・補正するための音信号Z”を生成し、この音信号Z”をスピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)からバルコニー86下の受聴エリア91に向けて出力させる装置である。本実施形態では、音場支援装置10による音信号Z”の生成を開始する前に、音信号Z”を生成するための信号処理の内容を決定するための設定作業が行われる。音場支援装置10の構成とその設定作業の詳細については、後述する。
The sound
また、図2(A)及び図2(B)に示す音響空間80では、ステージ81における音源90の位置と同等とみなせる位置とバルコニー86における天井82側の面79に音信号発生装置99と指向性マイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)が各々配置されている。音信号発生装置99とマイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)は、音場支援装置10の設定作業の際に利用され、設定作業の終了後に撤去されるものである。設定作業の際は、音信号発生装置99は、音場支援装置10による制御の下、各種音信号を発生する。マイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)は、音信号発生装置99により発生された音を収音し、収音した音を示す音信号を音場支援装置10に供給する。詳しくは、後述する。
Also, in the
次に、音場支援装置10の構成について説明する。音場支援装置10は、マイクロホン1L−1,1L−2,1R−1,1R−2及びスピーカ2L−k(k=1〜8),2R−k(k=1〜8)と接続されている。音場支援装置10は、当該音場支援装置10の設定作業の際は、ステージ81上及びバルコニー86上に配置される音信号発生装置99及びマイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)とも接続される。音場支援装置10において、マイクロホン1L−1及び1L−2から入力される各音信号Sは、A/D変換器(不図示)にてディジタル形式に変換され、加算部11Lによって加算される。マイクロホン1R−1及び1R−2から入力される各音信号Sは、A/D変換器(不図示)にてディジタル形式に変換され、加算部11Rによって加算される。加算部11Lの出力信号SはPEQ(Parametric Equalizer)12Lに入力され、加算部11Rの出力信号SはPEQ12Rに入力される。
Next, the configuration of the sound
PEQ12L及び12Rは、マイクロホン1L−1,1L−2,1R−1,及び1R−2の音信号Sに対して、音響空間80→マイクロホン1L−1,1L−2,1R−1及び1R−2→音場支援装置10→スピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)→音響空間80を一巡する音響帰還系の伝達関数(周波数応答)の振幅特性に急峻なピークが発生するのを防止するための補正を施す装置である。
The
PEQ12L及び12Rの伝達関数(周波数応答)の振幅特性は、PEQパラメータps(s=1〜8)によって決まる。PEQパラメータps(s=1〜8)は、音信号Sの帯域を8つのバンドに分割した場合における各バンドの中心周波数cfs(s=1〜8)、ゲインgs(s=1〜8)(各バンドの中心周波数cfs(s=1〜8)におけるレベルの増減量)、及びQ値qs(s=1〜8)(各バンドの周波数特性のピーク/ディップの鋭さ)を指定するパラメータである。PEQパラメータps(s=1〜8)はCPU25によって決定される。詳しくは、後述する。
The amplitude characteristics of the transfer functions (frequency response) of the
PEQ12Lによって補正された音信号S’は、スイッチ14L及び加算部15Lを経由してから8チャネルの音信号S’に分岐され、分岐された8チャネルの音信号S’の各々はFIR(Finite Impulse Response)フィルタ16L−k(k=1〜8)に入力される。PEQ12Rによって補正された音信号S’は、スイッチ14R及び加算部15Rを経由してから8チャネルの音信号S’に分岐され、分岐された8チャネルの音信号S’の各々はFIRフィルタ16R−k(k=1〜8)に入力される。加算部15L及び15Rは、ノイズジェネレータ30からテスト音信号が出力されている場合に、その信号をFIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)に供給する役割を果たす。このノイズジェネレータ30の出力信号については、後述する。
The sound signal S ′ corrected by the
FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)は、PEQ12L及び12Rの出力信号S’とフィルタ係数hm(m=0,1,2…M)を畳込む畳込み演算を行い、この演算結果を示す畳込み信号ZをPEQ17L−k(k=1〜),17R−k(k=1〜8)及び遅延部18L−k(k=1〜),18R−k(k=1〜8)とパワーアンプ部19L,19Rとを介してスピーカ2L−k(k=1〜8),2R−k(k=1〜8)へ供給する。これにより、スピーカ2L−k(k=1〜8),2R−k(k=1〜8)からは、受聴エリア91における空間感を補強・補正する音が、上述した反射音伝搬経路REFに沿って受聴エリア91に向けて出力される。ここで、FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)による畳込み演算に使用されるフィルタ係数hm(m=0,1,2…M)は、CPU25によって決定される。詳しくは、後述する。
The FIR filters 16L-k (k = 1 to 8) and 16R-k (k = 1 to 8) output the output signals S ′ and filter coefficients h m (m = 0, 1, 2,... M) of the
FIRフィルタ16L−k(k=1〜8)における各FIRフィルタ16L−kが出力した畳込み音信号ZはPEQ17L−kに入力される。また、FIRフィルタ16R−k(k=1〜8)における各FIRフィルタ16R−kが出力した畳込み音信号ZはPEQ17R−kに入力される。PEQ17L−k(k=1〜8)は、各FIRフィルタ16L−kが出力した音信号Zに対してスピーカ2L−k(k=1〜8)の設置環境を考慮した補正を施す装置である。PEQ17R−k(k=1〜8)は、各FIRフィルタ16R−kが出力した音信号Zに対してスピーカ2R−k(k=1〜8)の設置環境を考慮した補正を施す装置である。PEQ17L−k(k=1〜8)及び17R−k(k=1〜8)のPEQパラメータpsは、操作部31の操作によって書き換えられる。
The convolution sound signal Z output from each
PEQ17L−k(k=1〜8)における各PEQ17L−kによって補正された音信号Z’は遅延部18L−kに入力される。また、PEQ17R−k(k=1〜8)における各PEQ17R−kによって補正された音信号Z’は遅延部18R−kに入力される。遅延部18L−k(k=1〜8)及び18R−k(k=1〜8)は、PEQ17L−k(k=1〜8)及び17R−k(k=1〜8)の出力信号Z’に、マイクロホン1L−1,1L−2,1R−1,及び1R−2からスピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)までの音の伝搬時間に相当する遅延を与える装置である。遅延部18L−k(k=1〜8)及び18R−k(k=1〜8)の遅延量は、CPU25によって設定される。詳しくは、後述する。
The sound signal Z ′ corrected by each
遅延部18L−k(k=1〜8)における各遅延部18L−kによって遅延が与えられた音信号Z”は、パワーアンプ部19Lにより増幅された後、D/A変換器(不図示)にてアナログ信号に変換され、変換されたアナログ信号がスピーカ2L−kに供給される。また、遅延部18R−k(k=1〜8)における各遅延部18R−kによって遅延が与えられた音信号Z”は、パワーアンプ部19Rより増幅された後、D/A変換器(不図示)にてアナログ信号に変換され、変換されたアナログ信号がスピーカ2L−kに供給される。そして、スピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)は、各々に供給されたアナログ信号を音として出力する。
The sound signal Z ″ delayed by each
CPU25は、当該音場支援装置10の制御中枢である。CPU25は、RAM26をワークエリアとして利用しつつROM27に記憶された音場支援プログラム28を実行する。音場支援プログラム28は、次の2つの機能を有する。
a1.インパルス応答算出機能
これは、反射音伝搬経路REFにおける音源90とバルコニー86との間の区間のインパルス応答UpsideIR、及び音源90とマイクロホン1L−1,1L−2,1R−1,及び1R−2との間の区間のインパルス応答IRssを求める機能である。
b1.フィルタ係数列決定機能
これは、当該フィルタ係数hm(m=0,1,2…M)とインパルス応答IRssとの畳込み演算結果であるインパルス応答の実効値波形の傾きがインパルス応答UpsideIRの実効値波形の傾きに近似するフィルタ係数hm(m=0,1,2…M)を決定し、このフィルタ係数hm(m=0,1,2…M)をFIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)に設定する機能である。
The
a1. Impulse response calculation function This is the impulse response UpsideIR of the section between the
b1. Filter coefficient sequence determination function This is because the slope of the effective value waveform of the impulse response, which is the result of convolution of the filter coefficient h m (m = 0, 1, 2... M) and the impulse response IRss, is the effective of the impulse response UpsideIR. A filter coefficient h m (m = 0, 1, 2,... M) that approximates the slope of the value waveform is determined, and this filter coefficient h m (m = 0, 1, 2,... M) is determined as the
次に、本実施形態の動作を説明する。本実施形態の動作の特徴は、当該音場支援装置の設定作業の際にCPU25によって行われる処理にある。図3は、設定作業の際にCPU25が行う一連の処理を示すフローチャートである。ここで、図3に示す一連の処理のうちステップS100及びステップS120は、インパルス応答算出機能の働きによって実行される処理である。また、ステップS130(S131〜S136)は、フィルタ係数決定機能の働きによって実行される処理である。CPU25は、操作部31によって音場支援装置10の設定の指示が下されると、図3に示す処理を開始する。
Next, the operation of this embodiment will be described. The feature of the operation of the present embodiment is the processing performed by the
図3において、CPU25は、反射音伝搬経路REFにおける音源90とバルコニー86上の各マイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)との間の区間のインパルス応答UpsideIRを測定する(S100)。このステップS100では、ステージ81上における音源90の位置に置かれた音信号発生装置99を音源90とみなし、音信号発生装置99とバルコニー86の間の区間のインパルス応答をインパルス応答UpsideIRとする。より具体的に説明すると、CPU25は、スイッチ14L及び14Rをオフにした上で、音信号発生装置99から単位インパルスを発生させる。図4に示すように、信号発生装置99から単位インパルスを発生させると、天井82や壁83F,83B,83L及び83Rを経由することなくマイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)に到達した直接音D−SNDと天井82や壁83F,83B,83L及び83Rを経由してからマイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)に到達した反射音R−SNDとがマイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)によって収音される。CPU25は、単位インパルスを発生させた後にマイクロホン3L−k(k=1〜8)及び3R−k(k=1〜8)の各々から出力された音信号(直接音D−SNDと反射音R−SNDとを含む音信号)をインパルス応答UpsideIRとする。
In FIG. 3, the
次に、CPU25は、遅延量設定処理を行う(S110)。遅延量設定処理では、CPU25は、音源90からスピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)の各々までの音波の伝搬時間を個別に計測し、計測した伝搬時間を遅延部18L−k(k=1〜8)及び18R−k(k=1〜8)の遅延量として設定する。
Next, the
次に、CPU25は、インパルス応答IRssを測定する(S120)。このステップS120では、ステージ81上における音源90の位置に置かれた音信号発生装置99を音源90とみなし、音信号発生装置99とマイクロホン1L−1,1L−2,1R−1,及び1R−2の各々の間の区間のインパルス応答をインパルス応答IRssとする。より具体的に説明すると、CPU25は、スイッチ14L及び14Rをオフにした上で、音信号発生装置99から単位インパルスを発生させる。図5に示すように、音信号発生装置99から単位インパルスを発生させると、天井82や壁83F,83B,83L及び83Rを経由することなくマイクロホン1L−1,1L−2,1R−1,及び1R−2に到達した直接音D−SNDと天井82や壁83F,83B,83L及び83Rを経由してからマイクロホン1L−1,1L−2,1R−1,及び1R−2に到達した反射音R−SNDとがマイクロホン1L−1,1L−2,1R−1,及び1R−2によって収音される。CPU25は、単位インパルスを発生させた後、PEQ12Lの後段の信号採取点13LまたはPEQ12Rの後段の信号採取点13R(例えば、信号採取点13R)から音信号S(直接音D−SNDと反射音R−SNDとを含む音信号S)を取り込み、この取り込んだ音信号Sをインパルス応答IRssとする。図6(A)は、インパルス応答IRssの波形を示す図である。
Next, the
CPU25は、フィルタ係数列決定処理を行う(S130)。フィルタ係数列決定処理は、FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)による畳込み演算に使用するフィルタ係数hm(m=0,1,2…M)を決定する処理である。ここで、ステップS100で測定したインパルス応答UpsideIRは、音響空間80における音源90とバルコニー86の間の区間のインパルス応答である。よって、原理的には、音源90が発生した音とフィルタ係数hm(m=0,1,2…M)の畳込み演算結果である畳込み音信号Zをスピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)から受聴エリア91に音として出力させることにより、バルコニー86を無くした場合と同じ空間感が創出されるはずである。しかし、本実施形態では、音源90が発生した音を収音する収音手段としてマイクロホン1L−1,1L−2,1R−1及び1R−2を利用しているので、FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)による畳込み演算結果である畳込み音信号Zは音源90が発生した音の直接音D−SNDだけでなくその反射音R−SNDまでも冗長に畳み込まれたものとなる。このため、この畳込み音信号Zを受聴エリア91に音として出力した場合に同エリア91において創出される空間感は、依然として違和感を持ったものとなってしまう。
The
そこで、フィルタ係数列決定処理では、当該フィルタ係数hm(m=0,1,2…M)とインパルス応答IRssとの畳込み演算結果であるインパルス応答の実効値波形の傾きがインパルス応答UpsideIRの実効値波形の傾きに近似するフィルタ係数hm(m=0,1,2…M)を求め、このフィルタ係数hm(m=0,1,2…M)をFIRフィルタ16L−k及び16R−kの畳込み演算に使用するフィルタ係数hm(m=0,1,2…M)とする。 Therefore, in the filter coefficient sequence determination process, the slope of the effective value waveform of the impulse response that is the result of convolution of the filter coefficient h m (m = 0, 1, 2,... M) and the impulse response IRss is the impulse response UpsideIR. A filter coefficient h m (m = 0, 1, 2,... M) approximating the slope of the effective value waveform is obtained, and this filter coefficient h m (m = 0, 1, 2,... M) is obtained from the FIR filters 16L-k and 16R. The filter coefficient h m (m = 0, 1, 2,... M) used for the −k convolution operation is used.
具体的には、以下のような処理を行う。まず、CPU25は、ステップS100で求めたインパルス応答UpsideIRの実効値波形RMS(UpsideIR)を算出する(S131)。このステップS131では、実効値波形RMS(UpsideIR)を次のようにして算出する。まず、CPU25は、インパルス応答UpsideIRにおける単位インパルスの発生時刻t=0以降の各時刻tの振幅を各々2乗してエネルギーに変換する。次に、これらのエネルギーをτ(例えば、τ=25ミリ秒とする)秒分ずつ移動平均した移動平均エネルギーを求める。そして、各時刻tにおける移動平均エネルギーからなる波形を実効値波形RMS(UpsideIR)とする。
Specifically, the following processing is performed. First, the
次に、CPU25は、ステップS100で求めた16種類のインパルス応答UpsideIRのうち1つ(例えば、音源90及びマイクロホン3L−1の間の区間のインパルス応答UpsideIR)を代表値として選択し、このインパルス応答UpsideIRにおける直接音に相当する波形の部分を除いたインパルス応答UpsideIR’を算出する(S132)。より具体的には、インパルス応答UpsideIRにおける各時刻t=0,t=1,t=2…t=Mにオフセット時間TOFFSETを加算し、インパルス応答UpsideIRにおける各時刻t=0+TOFFSET,t=1+TOFFSET,t=2+TOFFSET,t=M+TOFFSET…の振幅をインパルス応答UpsideIR’におけるt=0,t=1,t=2…t=Mの振幅とする。図6(B)は、インパルス応答UpsideIR’の波形を示す図である。
Next, the
CPU25は、FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)のうちの1つを設定対象として選択する(S133)。
The
ここで、上述したように、本実施形態では収音手段としてマイクロホン1L−1,1L−2,1R−1及び1R−2を使用する。図6(C)に示すように、音源90とこのマイクロホン1L−1,1L−2,1R−1及び1R−2の間のインパルス応答IRssとインパルス応答UpsideIR’の畳込み演算結果であるインパルス応答IRss*UpsideIR’の傾きは、インパルス応答UpsideIR’の傾きよりも緩いものとなる。そして、このインパルス応答IRss*UpsideIR’の傾きの緩さが、受聴エリア91における空間感を違和感を持ったものとしてしまう原因となる。
Here, as described above, in this embodiment, the
そこで、CPU25は、ステップS134では、インパルス応答UpsideIR’を減衰率Di(i=1,2,3,4:ただし、D1<D2<D3<D4)によって減衰させた4種類のインパルス応答UpsideIR’×Di(i=1,2,3,4:ただし、D1<D2<D3<D4)を求め、それらの4種類の中から当該インパルス応答UpsideIR’×Diとインパルス応答IRssとを畳込んだ演算結果であるインパルス応答(UpsideIR’×Di)*IRssの実効値波形RMS((UpsideIR’×Di)*IRss)の傾きが実効値波形RMS(UpsideIR)の傾きに最も近くなるインパルス応答UpsideIR’×Diを選択し、この選択したインパルス応答UpsideIR’×Diを模擬したフィルタ係数hm(m=0,1,2…M)を求める(S134)。
Therefore, in step S134, the
ここで、本実施形態における減衰率Diは、実効値波形RMSの終端での値を減衰させる定数値である。本実施形態において、実効値波形RMSを減衰値Diだけ減衰させるとは、当該実効値波形RMSの終端の振幅をその振幅に減衰率Diを乗算した値とし、当該実効値波形RMSの始点と終端の間の各振幅を、終端における減衰率Diを乗算した値に応じた内挿(比例配分)により求まる値にすることを意味する。 Here, the attenuation factor D i in this embodiment is a constant value for attenuating the value at the end of the effective value waveform RMS. In the present embodiment, to attenuate the effective value waveform RMS by the attenuation value D i means that the amplitude at the end of the effective value waveform RMS is a value obtained by multiplying the amplitude by the attenuation factor D i, and the starting point of the effective value waveform RMS and each amplitude between the end, which means that the value obtained by interpolation in accordance with the value obtained by multiplying the attenuation factor D i at the end (prorated).
ステップS134についてより詳細に説明すると、CPU25は、ステップS131における実効値波形RMS(UpsideIR)の算出の手順と同様の手順に従って、実効値波形RMS((UpsideIR’×D1)*IRss),RMS((UpsideIR’×D2)*IRss),RMS((UpsideIR’×D3)*IRss),及びRMS((UpsideIR’×D4)*IRss)を各々算出する。次に、CPU25は、これらの各々について、当該実効値波形RMS((UpsideIR’×Di)*IRss)の各時刻tのエネルギーを各々の平方根であるレベルXtに変換する。また、実効値波形RMS(UpsideIR)についても同様に、各時刻tのエネルギーを各々の平方根であるレベルX0tに変換する。そして、CPU25は、実効値波形RMS((UpsideIR’×D1)*IRss),RMS((UpsideIR’×D2)*IRss),RMS((UpsideIR’×D3)*IRss),及びRMS((UpsideIR’×D4)*IRss)の各々について、当該実効値波形RMS((UpsideIR’×Di)*IRss)の始点の時刻t=0から終端の時刻t=tmaxまでの各時刻tのレベルXtと、実効値波形RMS(UpsideIR)の各時刻tのレベルX0tとを次式に代入することにより、当該実効値波形RMS((UpsideIR’×Di)*IRss)におけるレベルXtと実効値波形RMS(UpsideIR)におけるレベルX0tの誤差の2乗和MinValueを算出する。
CPU25は、実効値波形RMS((UpsideIR’×D1)*IRss),RMS((UpsideIR’×D2)*IRss),RMS((UpsideIR’×D3)*IRss),及びRMS((UpsideIR’×D4)*IRss)のうち式(1)により求めた誤差の2乗和MinValueが最小であった実効値波形RMS((UpsideIR’×Di)*IRss)を選択する。CPU25は、この選択した実効値波形RMS((UpsideIR’×Di)*IRss)における減衰率Diによってインパルス応答UpsideIR’を減衰させたインパルス応答UpsideIR’×Diを取得し、このインパルス応答UpsideIR’×Diにおける時刻t=0,t=1,t=2…t=Mの振幅を、設定対象のFIRフィルタ16L−k(または16R−k)の畳込み演算に用いるフィルタ係数hm(m=0,1,2…M)とする。
The
図7は、減衰率D1をD1=30(dB),減衰率D2をD2=40(dB),減衰率D3をD3=50(dB),減衰率D4をD4=60(dB)とした場合における実効値波形RMS((UpsideIR’×D1)*IRss),RMS((UpsideIR’×D2)*IRss),RMS((UpsideIR’×D3)*IRss),及びRMS((UpsideIR’×D4)*IRss)と、インパルス応答UpsideIR’とインパルス応答IRssとの畳込み演算結果であるインパルス応答UpsideIR’*IRssについて求めた実効値波形RMS(UpsideIR’*IRss)と、インパルス応答UpsideIRについて求めた実効値波形RMS(UpsideIR)とを示す図である。この図7の例では、インパルス応答UpsideIR’を減衰率D4により減衰させたものの実効値波形RMS((UpsideIR’×D4)*IRss)の傾きが、実効値波形RMS(UpsideIR)の傾きに最も近くなっている。よって、この例の場合、ステップS134では、インパルス応答UpsideIR’×D4を模擬したフィルタ係数hm(m=0,1,2…M)が算出される。 FIG. 7 shows that the attenuation rate D 1 is D 1 = 30 (dB), the attenuation rate D 2 is D 2 = 40 (dB), the attenuation rate D 3 is D 3 = 50 (dB), and the attenuation rate D 4 is D 4. = RMS ((UpsideIR ′ × D 1 ) * IRss), RMS ((UpsideIR ′ × D 2 ) * IRss), RMS ((UpsideIR ′ × D 3 ) * IRss) , And RMS ((UpsideIR ′ × D 4 ) * IRss) and the impulse response UpsideIR ′ * IRss which is the result of convolution of the impulse response UpsideIR ′ and the impulse response IRss, the RMS waveform RMS (UpsideIR ′ * IRss ) And the RMS waveform RMS (UpsideIR) obtained for the impulse response UpsideIR. In the example of FIG. 7, 'although the attenuated by the attenuation factor D 4 rms waveform RMS ((UpsideIR' impulse response UpsideIR the inclination of × D 4) * IRss), the inclination of the effective value waveform RMS (UpsideIR) Is closest. Thus, in this example, in step S134, the filter coefficient simulating the impulse response UpsideIR '× D 4 h m ( m = 0,1,2 ... M) is calculated.
CPU25は、ステップS134においてフィルタ係数hm(m=0,1,2…M)を算出した後、算出したフィルタ係数hm(m=0,1,2…M)を設定対象のFIRフィルタ16L−k(または16R−k)に設定する(S135)。次に、CPU25は、FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)のうち未だ設定対象となっていないものがあるかを判断する(S136)。そして、CPU25は、未だ設定対象となっていない一又は複数のFIRフィルタ16L−k及び16R−kがある場合(S135:Yes)、ステップS133に戻って、それらのうち一つを新たな設定対象として選択し、以降の処理を繰り返す。
CPU25, after calculating the filter coefficients h m (m = 0,1,2 ... M ) at step S134, the calculated filter coefficients h m (m = 0,1,2 ... M ) the setting target of the
CPU25は、ステップS136において、未だ設定対象となっていないものがないと判断した場合(S136:No)、PEQパラメータ決定処理を行う(S140)。PEQパラメータ決定処理は、PEQパラメータps(s=1〜8)を決定し、決定したPEQパラメータps(s=1〜8)をPEQ12L及び12Rに設定する処理である。
If the
PEQパラメータ決定処理では、CPU25は、次のような処理を行う。まず、CPU25は、スイッチ14L及び14Rをオフにした上で、ノイズジェネレータ30から加算部15L及び15Rにテスト音信号を供給させる。そして、加算部15L及び15R→FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)→PEQ17L−k(k=1〜8)及びPEQ17R−k(k=1〜8)→遅延部18L−k(k=1〜8)及び遅延部18R−k(k=1〜8)→パワーアンプ部19L及び19R→スピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)→音響空間80→マイクロホン1L−1,1L−2,1R−1,及び1R−2→加算部11L及び11R→PEQ12L及び12Rを経由して信号採取点13L及び13Rに伝搬された応答音信号を取り込む。
In the PEQ parameter determination process, the
CPU25は、テスト音信号と応答音信号の各々のパワースペクトルの差を、加算部15L及び15R→FIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)→PEQ17L−k(k=1〜8)及びPEQ17R−k(k=1〜8)→遅延部18L−k(k=1〜8)及び遅延部18R−k(k=1〜8)→パワーアンプ部19L及び19R→スピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)→音響空間80→マイクロホン1L−1,1L−2,1R−1,及び1R−2→加算部11L及び11R→PEQ12L及び12Rの系の伝達関数(周波数応答)の振幅特性R(ω)とする。
The
CPU25は、振幅特性R(ω)における全帯域のレベルを目標レベルLEV(例えば、LEV=6(dB))以下にするために必要な目標振幅特性C(ω)を求め、PEQ12L及び12Rの伝達関数(周波数応答)の振幅特性G(ω)を目標振幅特性C(ω)に最も近づけるために必要なPEQパラメータps(s=1〜8)を算出し、このPEQパラメータps(s=1〜8)をPEQ12L及び12Rに設定する。振幅特性G(ω)を目標振幅特性C(ω)に最も近づけるために必要なPEQパラメータps(s=1〜8)の算出の具体的な手順については、本発明との関連性が薄いため、この説明を割愛する。詳しくは、特許文献3を参照されたい。
The
以上説明した本実施形態では、当該フィルタ係数hm(m=0,1,2…M)とインパルス応答IRssとの畳込み演算結果であるインパルス応答の実効値波形の傾きがインパルス応答UpsideIRの実効値波形の傾きに近似するフィルタ係数hm(m=0,1,2…M)を求め、このフィルタ係数hm(m=0,1,2…M)をFIRフィルタ16L−k及び16R−kの畳込み演算に使用するフィルタ係数hm(m=0,1,2…M)とする。これにより、音響空間80におけるバルコニー86の下の受聴エリア91においてバルコニー86がない場合と同様の空間感が模擬される。すなわち、本実施形態では、音源90とバルコニー86上のスピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)との間のインパルス応答UpsideIRから直接音の波形部分を除去したインパルス応答UpsideIR’を求め、インパルス応答UpsideIR’を減衰率Diにより減衰させたインパルス応答UpsideIR’×Diを、FIRフィルタ16L−k及び16R−kの畳込み演算に用いるフィルタ係数hm(m=0,1,2…M)とする。これにより、元々のインパルス応答UpsideIRの反射音群の微視的構造(反射音群の相対的な関係)を保存しつつ、インパルス応答UpsideIRの反射音群の巨視的構造(包絡線構造)を修正したフィルタ係数hm(m=0,1,2…M)が得られる。そして、このようにして生成したフィルタ係数hm(m=0,1,2…M)を、FIRフィルタ16L−k及び16R−kの畳込み演算に使用するフィルタ係数hm(m=0,1,2…M)とする。よって、本実施形態によると、単なる残響付加装置と異なり、受聴エリア91にいる者にバルコニー86がない場合と同様の空間感を感じさせることができる。
In the present embodiment described above, the slope of the effective value waveform of the impulse response, which is the result of convolution of the filter coefficient h m (m = 0, 1, 2... M) and the impulse response IRss, is the effective of the impulse response UpsideIR. A filter coefficient h m (m = 0, 1, 2,... M) approximating the slope of the value waveform is obtained, and this filter coefficient h m (m = 0, 1, 2,... M) is obtained from the FIR filters 16L-k and 16R-. to use the convolution operation of k and filter coefficient h m (m = 0,1,2 ... M ). As a result, a spatial feeling similar to that in the case where there is no
以上、この発明の一実施形態について説明したが、この発明には他にも実施形態があり得る。例えば、以下の通りである。
(1)上記実施形態において、音源90の近傍のマイクロホン1L−1,1L−2,1R−1,1R−2の数を1個〜3個にしてもよいし、5個以上にしてもよい。また、バルコニー86における各面79及び78に各々設けられたマイクロホン3L−k及び3R−kとスピーカ2L−k及び2R−kを1〜7対にしてもよいし、9対以上にしてもよい。また、音場支援装置10内におけるFIRフィルタ16L−k(k=1〜8)及び16R−k(k=1〜8)の個数を17個以上にしてもよいし、15個以下にしてもよい。また、スピーカ2L−k(k=1〜8)及び2R−k(k=1〜8)の個数を17個以上にしてもよいし、15個以下にしてもよい。また、PEQ12L及び12Rにおけるバンド数は9バンド以上であってもよいし、7バンド以下であってもよい。
Although one embodiment of the present invention has been described above, the present invention may have other embodiments. For example, it is as follows.
(1) In the above embodiment, the number of
(2)上記実施形態は、本発明を、バルコニー86の下の受聴エリア91にいる者の空間感の補強に適用したものであった。しかし、音源90から受聴エリア91に至る直接音の伝搬を妨げず、音源90から音響空間80の内郭面を介して受聴エリア91に至る反射音の伝搬を妨げる位置にある別の反射音遮蔽体がある場合に、その反射音遮蔽体の背後にいる者の空間感の補強に本発明を適用してもよい。図8は、本発明を適用する音響空間80”を天井側から見た平面図である。この音響空間80”では、ステージ81”と対向する壁83B”からステージ81”に向かって仕切板87が突出している。このため、ステージ81”上の音源90から見て仕切板87の背後の受聴エリア91の者が感じる空間感が貧弱で違和感を有するものとなる。これに対し、上記実施形態の音場支援装置10によって生成した音信号Zを仕切板87における受聴エリア91側に設けたスピーカ(不図示)から出力することにより、この受聴エリア91の者が感じる空間感を仕切板87がない場合と同じ程度にまで補強することが可能である。この実施形態は、ホール後方から前方客席に上っていく「階段エリア(通称ダグアウト)」と呼ばれるエリアにおける空間感の補強にも好適である。
(2) In the above embodiment, the present invention is applied to the reinforcement of the spatial feeling of the person in the listening
(3)上記実施形態におけるフィルタ係数列決定処理では、当該インパルス応答UpsideIR’を減衰させる減衰率Diの候補を4種類とし、この4種類の減衰率Diのうちの最適値を選択した。しかし、インパルス応答UpsideIR’を減衰させる減衰率Diの候補を2種類または3種類としてもよいし、5種類以上としてもよい。要するに、減衰率Diの候補の数は任意の有限値にするとよい。 (3) in the filter coefficient sequence determination process in the above embodiment, the four candidates for attenuation rate D i attenuate the impulse response UpsideIR ', and select the optimal value of the four damping factor D i. However, to the candidate of the attenuation factor D i for attenuating impulse response UpsideIR 'may be two or three, may be 5 or more. In short, the number of attenuation rate D i candidates may be any finite value.
10…音場支援装置、1L,1R,3L,3R…マイクロホン、2L,2R…スピーカ、11L,11R…加算部、12L,12R,17L,17R…PEQ、13L,13R…信号採取点、14L,14R…スイッチ、15L,15R…加算部、16L,16R…FIRフィルタ、18L,18R…遅延部、19L,19R…パワーアンプ部、30…ノイズジェネレータ、31…操作部、99…音信号発生装置。
DESCRIPTION OF
Claims (2)
前記音源と前記反射音遮蔽体との間の区間のインパルス応答である第1のインパルス応答と、前記音源と前記収音手段との間の区間のインパルス応答である第2のインパルス応答とを求めるインパルス応答算出手段と、
前記ディジタルフィルタ手段の畳込み演算に用いるフィルタ係数列を決定する手段であって、前記第1のインパルス応答から直接音に相当する波形部分を除いたインパルス応答である第3のインパルス応答を求める手段を含み、前記第3のインパルス応答を複数種類の減衰率によって減衰させた複数種類のインパルス応答のうち、当該インパルス応答と前記第2のインパルス応答との畳込み演算結果であるインパルス応答の実効値波形の傾きが前記第1のインパルス応答の実効値波形の傾きに最も近くなるインパルス応答を選択し、選択したインパルス応答を模擬したフィルタ係数列を前記ディジタルフィルタ手段の畳込み演算に用いるフィルタ係数列とするフィルタ係数列決定手段と
を具備することを特徴とする音場支援装置。 A convolution operation is performed to convolve the sound signal of the sound collected by the sound collection means in the vicinity of the sound source in the acoustic space and the filter coefficient sequence, and the result of the convolution operation is calculated from the sound source in the acoustic space. Digital filter means for outputting a sound from the reflected sound shield on the propagation path of the reflected sound to the listening area to the listening area;
A first impulse response that is an impulse response in a section between the sound source and the reflection sound shielding body and a second impulse response that is an impulse response in a section between the sound source and the sound collecting means are obtained. Impulse response calculating means;
Means for determining a filter coefficient sequence used for the convolution operation of the digital filter means, and means for obtaining a third impulse response which is an impulse response obtained by removing a waveform portion corresponding to a direct sound from the first impulse response. An effective value of an impulse response that is a convolution operation result of the impulse response and the second impulse response among a plurality of types of impulse responses attenuated by a plurality of types of attenuation rates A filter coefficient sequence that selects an impulse response whose waveform slope is closest to the slope of the effective value waveform of the first impulse response, and uses a filter coefficient sequence that simulates the selected impulse response for the convolution operation of the digital filter means And a filter coefficient sequence determination means.
音響空間内の音源の近傍の収音手段により収音された音の音信号とフィルタ係数列とを畳込む畳込み演算を行い、この畳込み演算の演算結果を前記音源から前記音響空間内の受聴エリアに至る反射音の伝搬経路上にある反射音遮蔽体から前記受聴エリアに音として出力させるディジタルフィルタ手段と、A convolution operation is performed to convolve the sound signal of the sound collected by the sound collection means in the vicinity of the sound source in the acoustic space and the filter coefficient sequence, and the result of the convolution operation is calculated from the sound source in the acoustic space. Digital filter means for outputting a sound from the reflected sound shield on the propagation path of the reflected sound to the listening area to the listening area;
前記音源と前記反射音遮蔽体との間の区間のインパルス応答である第1のインパルス応答と、前記音源と前記収音手段との間の区間のインパルス応答である第2のインパルス応答とを求めるインパルス応答算出手段と、A first impulse response that is an impulse response in a section between the sound source and the reflection sound shielding body and a second impulse response that is an impulse response in a section between the sound source and the sound collecting means are obtained. Impulse response calculating means;
前記ディジタルフィルタ手段の畳込み演算に用いるフィルタ係数列を決定する手段であって、前記第1のインパルス応答から直接音に相当する波形部分を除いたインパルス応答である第3のインパルス応答を求める手段を含み、前記第3のインパルス応答を複数種類の減衰率によって減衰させた複数種類のインパルス応答のうち、当該インパルス応答と前記第2のインパルス応答との畳込み演算結果であるインパルス応答の実効値波形の傾きが前記第1のインパルス応答の実効値波形の傾きに最も近くなるインパルス応答を選択し、選択したインパルス応答を模擬したフィルタ係数列を前記ディジタルフィルタ手段の畳込み演算に用いるフィルタ係数列とするフィルタ係数列決定手段とMeans for determining a filter coefficient sequence used for the convolution operation of the digital filter means, and means for obtaining a third impulse response which is an impulse response obtained by removing a waveform portion corresponding to a direct sound from the first impulse response. An effective value of an impulse response that is a convolution operation result of the impulse response and the second impulse response among a plurality of types of impulse responses attenuated by a plurality of types of attenuation rates A filter coefficient sequence that selects an impulse response whose waveform slope is closest to the slope of the effective value waveform of the first impulse response, and uses a filter coefficient sequence that simulates the selected impulse response for the convolution operation of the digital filter means And a filter coefficient sequence determining means
を実現させるプログラム。A program that realizes
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010183078A JP5644253B2 (en) | 2010-08-18 | 2010-08-18 | Sound field support apparatus and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010183078A JP5644253B2 (en) | 2010-08-18 | 2010-08-18 | Sound field support apparatus and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012042661A JP2012042661A (en) | 2012-03-01 |
JP5644253B2 true JP5644253B2 (en) | 2014-12-24 |
Family
ID=45899074
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010183078A Expired - Fee Related JP5644253B2 (en) | 2010-08-18 | 2010-08-18 | Sound field support apparatus and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5644253B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106714040A (en) * | 2017-01-17 | 2017-05-24 | 宁波音王电声股份有限公司 | Digital variable reverberation audio system and realization method thereof |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7447533B2 (en) * | 2020-02-19 | 2024-03-12 | ヤマハ株式会社 | Sound signal processing method and sound signal processing device |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0583786A (en) * | 1991-09-20 | 1993-04-02 | Matsushita Electric Ind Co Ltd | Reflection sound extracting method and sound field correcting method |
JP3147618B2 (en) * | 1993-09-21 | 2001-03-19 | ヤマハ株式会社 | Acoustic characteristic correction device |
JP3957819B2 (en) * | 1997-06-26 | 2007-08-15 | 松下電器産業株式会社 | Sound field correction method and sound field correction apparatus |
-
2010
- 2010-08-18 JP JP2010183078A patent/JP5644253B2/en not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106714040A (en) * | 2017-01-17 | 2017-05-24 | 宁波音王电声股份有限公司 | Digital variable reverberation audio system and realization method thereof |
CN106714040B (en) * | 2017-01-17 | 2019-05-03 | 音王电声股份有限公司 | A kind of digital variable reverberant audio system and its implementation |
Also Published As
Publication number | Publication date |
---|---|
JP2012042661A (en) | 2012-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9432793B2 (en) | Head-related transfer function convolution method and head-related transfer function convolution device | |
JP5444472B2 (en) | Sound source separation apparatus, sound source separation method, and program | |
US8199921B2 (en) | Sound field controlling device | |
EP2320675B1 (en) | Audio processing device | |
US10115412B2 (en) | Signal processor with side-tone noise reduction for a headset | |
JP4984683B2 (en) | Sound emission and collection device | |
EP2268065A2 (en) | Audio signal processing device and audio signal processing method | |
JP6251054B2 (en) | Sound field correction apparatus, control method therefor, and program | |
CN1440628A (en) | Interference Suppression Technology | |
KR20180071395A (en) | Generating binaural audio in response to multi-channel audio using at least one feedback delay network | |
JP2016144112A (en) | Microphone selection device, microphone system and microphone selection method | |
JP5644253B2 (en) | Sound field support apparatus and program | |
WO2024169940A1 (en) | Apparatus and method for echo cancellation without reference loop | |
JP4552876B2 (en) | Audio signal transmitting / receiving apparatus and audio conference apparatus | |
US11895485B2 (en) | Sound signal processing method and sound signal processing device | |
JP6631657B2 (en) | Sound emission and collection device | |
JP5633145B2 (en) | Sound signal processing device | |
KR970004178B1 (en) | Audio echo sound adding device | |
EP3869500B1 (en) | Sound signal processing method and sound signal processing device | |
EP3869501B1 (en) | Sound signal processing method and sound signal processing device | |
JP5224613B2 (en) | Sound field correction system and sound field correction method | |
JP2013247524A (en) | Sound-reproducing device and method for controlling the same | |
JP2007318521A (en) | Sound emission/pickup apparatus | |
JPH06266379A (en) | Reverberation method | |
JP5633144B2 (en) | Sound signal processing device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130620 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140522 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140527 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140624 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20141007 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141020 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5644253 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |