JP2003005770A - 残響生成付加方法とその装置 - Google Patents

Abstract

(57)【要約】 【課題】より自然な残響音をより簡単な構成により生成
しソース音源に付加することのできる残響生成付加装置
とその方法を提供する。 【解決手段】残響音を表すインパルス応答を、後部と初
期部に分割して生成する。また、後部インパルス応答
は、固有周波数列の生成とそれらの固有周波数を振動数
とする減衰自由振動の生成の２段階の処理により行う。
したがって、まず隣接固有周波数間隔統計により室の固
有周波数数列を生成し（ステップＳ１）、これに基づい
て固有周波数減衰振動を生成して後部残響音を生成する
（ステップＳ２）。一方で、初期反射音を生成し（ステ
ップＳ３）、これら生成した後部残響音を初期反射音と
結合し（ステップＳ４）、さらに両耳間時間差を用いて
２チャンネル化を行う（ステップＳ５）。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、たとえばカーオー
ディオ機器などに適用して好適な、残響生成付加装置お
よびその方法に関する。

【０００２】

【従来の技術および発明が解決しようとする課題】従来
の代表的な残響付加装置は、図２３に示すようなＦＩＲ
ディジタルフィルタ形、あるいは図２４に示すような帰
還路を有するＩＩＲディジタルフィルタ形のいずれかに
大別される。図２３に示したＦＩＲ形フィルタの係数
は、コンサートホールで予め収録した残響特性から決定
されることが多い。したがって、この残響生成付加方式
では、自然な残響音が生成される代わりに、コンサート
ホールの残響特性実測データベースの利用に起因する制
限がある。

【０００３】また、図２４に示す方式では、実測データ
を不要とする代わりに、残響音生成アルゴリズムが未だ
研究開発対象であるため、自然な残響音を生成しにくい
という困難を有していた。特に、室容積、室形、広がり
感を含む室内音場パラメータによる残響音制御、自然な
初期反射音群の生成、原音の音色を損なわない自然な残
響音減衰特性と周波数特性の実現、規則的な繰り返しが
なく密度の高い残響反射音の生成、さらに、商用ステレ
オ音楽ソースの定位感を損なわずに豊かな広がり感を付
加した残響音の生成、というような未解決課題があり改
善の余地があった。

【０００４】したがって本発明の目的は、より自然な残
響音をより簡単な構成により生成しソース音源に付加す
ることのできる残響生成付加装置とその方法を生成する
ことにある。

【０００５】

【課題を解決するための手段】前記課題を解決するため
に、本発明の残響付加方法は、室内音場の固有周波数を
振動数とする減衰自由振動の重畳によって室内インパル
ス応答を表現、生成し、さらに、残響音の特徴に応じて
複数の初期反射音を加えることによって、室内インパル
ス応答を生成するようにした。そして、この生成したイ
ンパルス応答を音源信号に畳み込むことにより、残響音
を付加するようにした。また、特に本発明においては、
新たに開拓された 室内音場理論に基づいて、（ａ）減
衰自由振動の周波数選択、（ｂ）減衰自由振動の振幅と
位相の配列、（ｃ）減衰自由振動の減衰特性の設定、
（ｄ）初期反射音の時間間隔、減衰特性と周波数特性の
制御、（ｅ）初期残響音における微小な乱反射音の生
成、（ｆ）反射音の到来方向と両耳効果に基づくステレ
オ化処理、を行うことにより、従来の課題を解決し好ま
しい残響音を付加するようにした。

【０００６】したがって、本発明の残響付加方法は、室
内音場の固有周波数を振動数とし、振幅、位相および減
衰の各特性が各々所定の特性に規定された減衰自由振動
を後部残響音として生成し、振幅、減衰、周波数および
遅れ時間間隔分布の各特性が各々所望の特性に規定され
た初期反射音を生成し、前記生成された後部残響音およ
び初期反射音、および直接音を合成して室内インパルス
応答を生成し、前記生成された室内インパルス応答を音
源信号に付加する。

【０００７】好適には、前記後部残響音は、前記初期反
射音と前記後部残響音のエネルギー比、残響時間周波数
特性、共振（固有）周波数分布、減衰自由振動振幅およ
び位相分布、反射音到来方向分布に基づいて確率的に抽
出された共振周波数を振動数とし、振幅がレーレー分
布、初期位相が一様分布として生成される指数減衰自由
振動である。

【０００８】また好適には、前記初期反射音は、要求さ
れる室内広さ感に応じて設定される前記遅れ時間間隔分
布、要求される残響付加音の明瞭度に応じて設定される
前記初期反射音の振幅および減衰特性および所望の特性
に設定される前記周波数特性に基づいて生成し、さらに
微小な乱反射音が付加され生成される。

【０００９】さらに好適には、周波数帯域別に、１つの
反射音群とみなせるサンプル個数の遅れ時間を、反射音
の到来方向確率分布と反射音の両耳間時間差に基づいて
変化させることにより、前記インパルス応答を確率的に
複数チャンネル化する

【００１０】また好適には、前記室内インパルス応答の
付加は、帯域合成後部残響波形、帯域合成後部残響波形
の包絡線と正弦波搬送波、あるいは、後部帯域残響波形
の包絡線のいずれかを保存し、かつ波形を二乗誤差尺度
で最良近似する正弦波搬送波を用いて行う。

【００１１】また、本発明の残響生成付加装置は、室内
音場の固有周波数を振動数とし、振幅、位相および減衰
の各特性が各々所定の特性に規定された後部残響音を生
成する後部残響音生成手段と、振幅、減衰、周波数、遅
れ時間間隔分布の各特性が各々所望の特性に規定された
初期反射音を生成する初期反射音生成手段と、前記生成
された後部残響音および初期反射音、および、直接音を
合成して室内インパルス応答を生成する室内インパルス
応答生成手段と、前記生成された室内インパルス応答を
音源信号に付加する付加手段とを有する。

【００１２】

【発明の実施の形態】本発明の一実施の形態の残響生成
付加方法について、図１〜図２２を参照して説明する。
本実施の形態の残響生成付加方法は、商用音楽ソースを
含む２チャンネルの音源信号に、所望の残響・広がり感
を付加するものであって、まず、室内音場特性と両耳受
聴効果に基づいて残響音を室内のインパルス応答として
求め、生成されるインパルス応答を、音楽あるいは音声
信号を含む音源信号に畳み込むことによって所望の残響
を付加するものである。

【００１３】以下、その残響生成付加方法の基本的な処
理手順、すなわち残響音場のインパルス応答を生成する
処理の流れについて説明する。図１は、残響音場のイン
パルス応答を生成する処理の概略の流れを示すフローチ
ャートである。残響音場のインパルス応答は、直接音、
いくつかの独立した反射音である初期反射音、および、
それに後続する残響音群である後部残響音に分割して考
えることができ、この後部残響音は、主として音場の統
計的性質に基づいて特徴つけられる。したがって、本方
式では、まず、残響音を表すインパルス応答を、後部と
初期部に分割して生成する。また、後部インパルス応答
は、固有周波数列の生成とそれらの固有周波数を振動数
とする減衰自由振動の生成の２段階の処理により行う。

【００１４】したがって、本実施の形態の方法において
は、図1 に示すように、まず隣接固有周波数間隔統計に
より室の固有周波数数列を生成し（ステップＳ１）、こ
れに基づいて固有周波数減衰振動を生成して後部残響音
を生成する（ステップＳ２）。一方で、初期反射音を生
成し（ステップＳ３）、これら生成した後部残響音を初
期反射音と結合し（ステップＳ４）、さらに両耳間時間
差を用いて２チャンネル化を行う（ステップＳ５）。

【００１５】以下、各ステップの処理について詳細に説
明する。

【００１６】固有周波数列の生成 残響音場解析の手法としては、固有周波数感覚分布に関
する研究が古くから行われているが（M. R. Schroeder,
J. Acoust Soc. Am., 45, 545-565 (1969)）、近年、
隣接固有周波数の間隔分布と音線伝搬のカオス性との関
連が指摘され（M. R. Schroeder, J. Audio Eng. Soc.,
37, 795-808 (1989) ）、音線伝搬のカオス性の増大に
伴って隣接固有周波数間隔分布が自由度１から非整数自
由度を経て自由度２のガンマ分布へ推移する予測が述べ
られている。ここで、式（１）に示すような、平均値を
１とする自由度２のガンマ分布（Ｗｉｇｎｅｒ分布）に
着目する。

【００１７】

【数１】

【００１８】但し、ｘは式（２）で表される。

【００１９】

【数２】

【００２０】ここでｃは音速（ｍ／ｓ）、Ｖは室容積
（ｍ³ ）、ｆは周波数（Ｈｚ）を各々表す。また、ｄ
（ｆ）は平均固有周波数を表し、固有周波数の密度ｎ
（ｆ）の逆数として、式（３）により表される。

【００２１】

【数３】

【００２２】式（１）に基づく、不整形室における固有
周波数間隔分布を図２（Ａ）に示す。図２（Ａ）に示す
ように、不整形室における固有周波数間隔分布において
は、固有周波数の縮退が起こりにくいことがわかる。な
お、直方体室における固有周波数間隔分布は、自由度１
のガンマ分布、すなわちポアソン分布として与えられ、
その固有周波数間隔分布は、図２（Ｂ）のようになる。
この場合、図２（Ａ）に示した不整形音場と比較して、
固有周波数の縮退が起こりやすいことがわかる。本実施
の形態による残響生成方式においては、コンサートホー
ルのような複雑形状で囲まれた音場、すなわち図２
（Ａ）に示すような不整形室の固有周波数分布に着目す
る。

【００２３】式（３）は、固有周波数の密度が周波数の
二乗に比例して上昇することを示している。この周波数
の上昇に伴う密度の増大は音場の波動理論に基づく解析
を困難にする要因の１つとなっているが、音場統計論で
は多数の隣接する固有周波数による自由振動の重畳が基
本的原理になっている。固有周波数による自由振動の励
振特性が有する周波数帯域幅（Modal Bandwidth ）は、
式（４）で表される。ここで、Ｔ_R は残響時間（ｓ）、
Ｈ（ｆ）は減衰自由振動の励振特性である。

【００２４】

【数４】

【００２５】本実施の形態においては、図２（Ａ）に示
した隣接固有周波数間隔分布と式（２）に基づいて、式
（５）より固有周波数間隔δｆを抽出し、固有周波数列
ｐ（ｆ）を発生させる。

【００２６】

【数５】

【００２７】なお、周波数の増大とともに固有周波数間
隔が小さくなるため、ここでは固有周波数間隔が周波数
帯域幅（Modal Bandwidth ）より大きいものだけを抽出
するものとする。具体的には、フローチャートを図３に
示すように、まず図２（Ａ）に示したような隣接固有周
波数間隔のヒストグラムより隣接固有周波数間隔ｘを抽
出し（ステップＳ１１），前述した式（２）および式
（４）に基づいて固有周波数間隔δｆと周波数帯域幅Ｂ
を求める（ステップＳ１２）。固有周波数間隔δｆが周
波数帯域幅Ｂより大きい場合には（ステップＳ１３）、
周波数ｆ＋δｆに対応する固有周波数列Ｐの値を１とす
ることにより、これを抽出する（ステップＳ１４）。

【００２８】そして、その周波数ｆ＋δｆが、標本化周
波数ｆ_s の１／２よりも大きい場合には（ステップＳ１
５）、それまで抽出した固有周波数をもって、固有周波
数列収集の処理を終了する（ステップＳ１７）。また、
ステップＳ１３において固有周波数間隔δｆが周波数帯
域幅Ｂ以下の場合、または、ステップＳ１５において周
波数ｆ＋δｆが標本化周波数ｆ_s の１／２以下の場合に
は、周波数ｆを固有周波数間隔δｆだけ増加させて（ス
テップＳ１６）、再びステップＳ１１以下の処理を繰り
返す。以上の処理により、固有周波数列ｐ（ｆ）を発生
させる。

【００２９】図４に、生成された固有周波数列の例を示
す。周波数の増大とともに固有周波数の間隔が狭まる状
態が理解できる。この固有周波数の密度の周波数に伴う
増加は、自然な残響音を生成する上で非常に重要であ
る。

【００３０】減衰自由振動の生成による後部残響音の生
成 前述したように生成された固有周波数列の各固有周波数
を振動数ｆ_i とし、式（６）に示すような減衰自由振動
を考え、これを重畳することにより、後部残響音のイン
パルス応答を生成する。

【００３１】

【数６】

【００３２】ここでＡ_i 及びθ_i はそれぞれ減衰振動の
初期振幅と初期位相を表す。また、αは振動の減衰係数
で、予め設定される残響時間に基づいて決定される。換
言すれば、着目する周波数帯域別に固有インパルス応答
の減衰特性を制御することによって、所望の残響特性が
実現される。また、初期振幅は確率変数Ａと表すことに
し、標準偏差σ、平均値０のガウス雑音の振幅分布を表
す式（７）に示すレーレー分布に従う母集団から抽出さ
れる。

【００３３】

【数７】

【００３４】また、初期位相角も確率変数θとして一様
分布する母集団から選び出す。全帯域にわたって自然な
残響減衰特性を実現するためには、それぞれの固有イン
パルス応答の振幅と位相特性を制御する必要がある。そ
のため本方式では、不規則音場の性質を鑑み、このよう
にレーレー分布および一様分布を用いる。

【００３５】このようにして得られた後部インパルス応
答、すなわち後部残響音の例を、図５に示す。また、こ
の後部残響音の振幅周波数特性を図６に、残響減衰曲線
を図７に、帯域別残響特性を図８に、残響時間周波数特
性を図９に各々示す。

【００３６】初期反射音の生成 初期反射音群は、反射音の到来時間間隔を確率変数とし
て、直接音から５０ｍｓの間に到来する反射音間隔を正
規乱数として生成する。本実施の形態においては、個々
の初期反射音に対して、初期反射音の遅れ時間間隔分
布、初期反射音の大きさと減衰特性、および、初期反射
音の周波数特性の３つ条件を設定することにより、生成
する初期反射音を規定する。

【００３７】具体的には、初期反射音の遅れ時間間隔分
布は、残響を生成する室内の広さ感に応じて、たとえば
図１０に示すような分布のいずれかを選択することによ
り設定する。また、初期反射音の大きさと減衰特性は、
音源の明瞭度に応じて、図１１に示すような初期減衰特
性のいずれかを選択することにより設定する。また、初
期反射音の周波数特性としては、たとえば図１２に示す
ような特性を設定する。このようにして得られた初期反
射音の例を、図１３に示す。また、この初期反射音の振
幅周波数特性を図１４示す。

【００３８】後部残響音と初期反射音の結合 図５に示したような後部残響音と、図１３に示したよう
な初期反射音を所望のエネルギー比をもって合成するこ
とにより、残響音を表すインパルス応答を生成する。後
部と初期部を合成したインパルス応答の例を図１５に、
また、この振幅周波数特性を図１６に、残響曲線を図１
７に示す。

【００３９】両耳間時間差による２チャンネル化 残響音は直接音に響きを付加する以外に、反射音の到来
方向による両耳間時間差の変化に基づく２チャンネル化
を施すことにより、音場の広がり感の制御が期待でき
る。本実施の形態の方法においては、両耳を２つのマイ
クロホンに置き換えるときに必要なマイクロホン距離、
すなわち等価両耳間距離を用いて、図１８に示すように
音波の到来方向によって決定される両耳間時間差τに基
づいて、周波数帯域別に２チャンネルの残響音をを生成
する。

【００４０】また、音像の定位感を安定させるため、周
波数帯域別に反射音群を定めて、反射音群ごとに図１９
に示すような反射音到来方向分布を生成する。図２０
は、音の到来方向と両耳間時間差に基づいて２チャンネ
ル化されたインパルス応答の例である。

【００４１】以上が、本実施の形態の残響生成付加方法
の基本的な処理手順の説明であるが、次に、このような
残響生成付加方法をコンピュータシステムにおいて、あ
るいは専用装置として適用する場合に好適な構成の一例
をを図２１および図２２に示す。図２１は、本実施の形
態の残響生成付加方法を具体的な装置に適用する場合の
処理構成の一例を模式的に示す図である。

【００４２】この例においては、まず、初期設定とし
て、室容積、室形、初期反射音到来間隔、初期減衰、初
期周波数特性、残響時間、残響周波数特性、反射音到来
方向、初期反射音と残響音のエネルギー比を設定する
(ステップＳ２１）。次に、後部残響音を設定するため
に、設定した室容積、室形および残響時間に基づいて共
振周波数分布を求め、これより減衰振動数を設定する
（ステップＳ２２）。

【００４３】次に、これを帯域分割し (ステップＳ２
３）、設定された残響時間および残響周波数特性に基づ
いて、特に本実施の形態においては振幅特性はレーレー
分布より、位相特性は一様乱数により求めるものとして
帯域ごとの減衰波形を生成し（ステップＳ２４）、ここ
では特に各帯域ごとに、正規分布と設定した反射音到来
方向および両耳間時間差に基づいて減衰波形の２ｃｈス
テレオ化を行う（ステップＳ２５）。そして、各帯域ご
との減衰波形を合成することにより、全域の後部残響波
形を生成する (ステップＳ２６）。

【００４４】一方、初期反射音を生成する処理として
は、設定された初期反射音到来間隔、すなわち本実施の
形態においては正規分布とした初期反射音到来間隔に基
づいて反射音を設定し (ステップＳ２７）、設定した初
期減衰、すなわち本実施の形態においては指数減衰とし
た初期減衰にもとづいて初期減衰波形を生成する（ステ
ップＳ２８）。次に、低域および高域通過フィルタを通
すことにより、同じく設定された周波数特性となるよう
に制御するとともに微小な乱反射音を生成し（ステップ
Ｓ２９）、さらに、後部残響波形と同様に、正規分布と
設定した反射音到来方向および両耳間時間差に基づいて
ステレオ化を行う（ステップＳ３０）。

【００４５】そして、２ｃｈ各々に初期反射音と後部残
響音を合成し、最後に、元の音源信号と合成して、残響
音を付加したインパルス応答を生成する (ステップＳ３
１）。

【００４６】また、図２２は、本実施の形態の残響生成
付加方法を適用した処理装置の具体的構成の一例を示す
図である。図２２に示す例においては、周波数帯域を１
／３２ｏｃｔ分析して２チャンネル化された帯域合成残
響波形を生成した後、全帯域に渡る重畳によって２チャ
ンネル全帯域合成残響波形を生成する。１／３２ｏｃｔ
分析されたそれぞれの帯域合成残響波形は、遅延加算装
置による初期反射音生成部と、振幅・位相・周波数・減
衰係数の異なる減衰自由振動発生器から構成される後部
残響音生成部と、２チャンネル可変遅延器を有する２チ
ャンネル残響音声部から構成される。

【００４７】それぞれの１／３２ｏｃｔ帯域内の処理は
全て同様である。帯域内の各処理において、初期反射音
生成部の初期反射音を生成する遅延加算装置は、図１０
〜図１２に示した初期反射音の遅れ時間間隔分布、初期
減衰特性ならびに初期反射音の周波数特性に基づいて、
遅延器、乗算器を設定する。後部残響音生成部は、減衰
自由振動発生器によって構成され、室容積を５００ｍ³
としたとき、図２（Ａ）に示した頻度分布に従って、図
３に示したような固有周波数による指数減衰自由振動を
生成する。指数減衰特性は、図１８の残響時間周波数特
性から設定する。

【００４８】次に、指数減衰自由振動の初期振幅を、レ
ーレー分布、位相を一様分布として、それぞれの指数減
衰自由振動を帯域別に（図２２の実施例では、１／３２
ｏｃｔ）合成したそれぞれの帯域合成残響波形を、最小
二乗誤差の意味で最良近似する帯域別減衰自由振動波形
を生成する。２チャンネル残響音生成部は、図１８に示
したような反射音帯域別到来方向分布と、同一方向から
到来する音波の帯域別継続時間に基づいて、両耳間時間
差を帯域別に設定する。全帯域合成部は、帯域別残響波
形の線形結合によって、全帯域残響波形を生成する。

【００４９】なお、本発明は本実施の形態に限られるも
のではなく、任意好適な種々の改変が可能である。たと
えば、固有周波数列を生成する際に固有周波数間隔が小
さくなった場合、本実施の形態においては、固有周波数
間隔が周波数帯域幅（Modal Bandwidth ）より大きいも
のだけを抽出するようにした。しかし、たとえば、生成
される固有周波数間隔の平均値が周波数帯域幅（Modal
Bandwidth ）より小さくなった場合には、周波数帯域幅
おきにランダムに固有周波数を抽出するようにしてもよ
い。得ようとする残響音の特徴や、実際に処理する際の
処理能力、処理精度などに応じて任意の方法で固有周波
数の抽出方法、さらには、図２１に示したような各処理
条件などを設定するようにしてよい。

【００５０】

【発明の効果】このように、本発明によれば、より自然
な残響音をより簡単な構成により生成しソース音源に付
加することのできる残響生成付加装置とその方法を生成
することができる。

【図面の簡単な説明】

【図１】図１は、残響音場のインパルス応答を生成する
処理の概略の流れを示すフローチャートである。

【図２】図２（Ａ）は、不整形室における固有周波数間
隔分布を示すヒストグラムであり、図２（Ｂ）は直方体
室における固有周波数間隔分布を示すヒストグラムであ
る。

【図３】図３は、固有周波数列を生成する処理の流れを
示すフローチャートである。

【図４】図４は、生成された固有周波数列の例を示す図
である。

【図５】図５は、生成された後部残響音を示す波形図で
ある。

【図６】図６は、図５に示した後部残響音の振幅周波数
特性を示す図である。

【図７】図７は、図５に示した後部残響音の残響減衰曲
線を示す図である。

【図８】図８は、図５に示した後部残響音の帯域別残響
特性を示す図である。

【図９】図９は、図５に示した後部残響音の残響時間周
波数特性を示す図である。

【図１０】図１０は、設定する初期反射音の遅れ時間間
隔分布を示す図である。

【図１１】図１１は、設定する初期反射音の減衰特性を
示す図である。

【図１２】図１２は、設定される初期反射音の周波数特
性の一例を示す図である。

【図１３】図１３は、生成された初期反射音を示す波形
図である。

【図１４】図１４は、図１３に示した初期反射音の振幅
周波数特性を示す図である。

【図１５】図１５は、後部と初期部を合成したインパル
ス応答を示す波形図である。

【図１６】図１６は、図１５に示した後部と初期部を合
成したインパルス応答の振幅周波数特性を示す図であ
る。

【図１７】図１７は、図１５に示した後部と初期部を合
成したインパルス応答の残響減衰曲線を示す図である。

【図１８】図１８は、等価両耳間距離を説明するための
図である。

【図１９】図１９は、帯域別に設定される反射音到来方
向分布を示す図である。

【図２０】図２０は、音の到来方向と両耳間時間差に基
づいて２チャンネル化されたインパルス応答の例であ
る。

【図２１】図２１は、本実施の形態の残響生成付加方法
を具体的な装置に適用する場合の処理手順の例を模式的
に示す図である。

【図２２】図２２は、本実施の形態の残響生成付加方法
を適用した処理装置の具体的構成の例を示す図である。

【図２３】図２３は、従来のＦＩＲディジタルフィルタ
形の残響付加装置の処理を説明するための図である。

【図２４】図２４は、従来のＩＩＲディジタルフィルタ
形の残響付加装置の処理を説明するための図である。

─────────────────────────────────────────────────────

【手続補正書】

【提出日】平成１３年７月１０日（２００１．７．１
０）

【手続補正１】

【補正対象書類名】図面

【補正対象項目名】図１

【補正方法】変更

【補正内容】

【図１】

【手続補正２】

【補正対象書類名】図面

【補正対象項目名】図３

【補正方法】変更

【補正内容】

【図３】

【手続補正３】

【補正対象書類名】図面

【補正対象項目名】図２１

【補正方法】変更

【補正内容】

【図２１】

───────────────────────────────────────────────────── フロントページの続き (72)発明者 風間 道子 東京都新宿区下落合３−４−７ (72)発明者 神谷 裕一 神奈川県川崎市宮前区鷺沼１−11−１ Ｄ ＩＫマンション308号 Ｆターム(参考） 5D108 AB07 AB08 AC02 AD04 AD05 AD07

Claims (6)

【特許請求の範囲】
1. 【請求項１】室内音場の固有周波数を振動数とし、振
   幅、位相および減衰の各特性が各々所定の特性に規定さ
   れた減衰自由振動を後部残響音として生成し、 振幅、減衰、周波数および遅れ時間間隔分布の各特性が
   各々所望の特性に規定された初期反射音を生成し、 前記生成された後部残響音および初期反射音、および直
   接音を合成して室内インパルス応答を生成し、 前記生成された室内インパルス応答を音源信号に付加す
   る残響生成付加方法。
2. 【請求項２】前記後部残響音は、前記初期反射音と前記
   後部残響音のエネルギー比、残響時間周波数特性、共振
   （固有）周波数分布、減衰自由振動振幅および位相分
   布、反射音到来方向分布に基づいて確率的に抽出された
   共振周波数を振動数とし、振幅がレーレー（Ｒａｙｌｅ
   ｉｇｈ）分布、初期位相が一様分布として生成される指
   数減衰自由振動である請求項１に記載の残響生成付加方
   法。
3. 【請求項３】前記初期反射音は、要求される室内広さ感
   に応じて設定される前記遅れ時間間隔分布、要求される
   残響付加音の明瞭度に応じて設定される前記初期反射音
   の振幅および減衰特性および所望の特性に設定される前
   記周波数特性に基づいて生成し、さらに微小な乱反射音
   が付加されて生成される請求項１または２に記載の残響
   生成付加方法。
4. 【請求項４】さらに、周波数帯域別に、１つの反射音群
   とみなせるサンプル個数の遅れ時間を、反射音の到来方
   向確率分布と反射音の両耳間時間差に基づいて変化させ
   ることにより、前記インパルス応答を確率的に複数チャ
   ンネル化する請求項１〜３のいずれかに記載の残響生成
   付加方法。
5. 【請求項５】前記室内インパルス応答の付加は、帯域合
   成後部残響波形、帯域合成後部残響波形の包絡線と正弦
   波搬送波、あるいは、後部帯域残響波形の包絡線のいず
   れかを保存し、かつ波形を二乗誤差尺度で最良近似する
   正弦波搬送波を用いて行う請求項１〜４のいずれかに記
   載の残響生成付加方法。
6. 【請求項６】室内音場の固有周波数を振動数とし、振
   幅、位相および減衰の各特性が各々所定の特性に規定さ
   れた後部残響音を生成する後部残響音生成手段と、 振幅、減衰、周波数、遅れ時間間隔分布の各特性が各々
   所望の特性に規定された初期反射音を生成する初期反射
   音生成手段と、 前記生成された後部残響音および初期反射音、および、
   直接音を合成して室内インパルス応答を生成する室内イ
   ンパルス応答生成手段と、 前記生成された室内インパルス応答を音源信号に付加す
   る付加手段とを有する残響生成付加装置。