JP2003157090A

JP2003157090A - 残響音生成方法及び残響音シミュレータ

Info

Publication number: JP2003157090A
Application number: JP2001354527A
Authority: JP
Inventors: Kuniaki Osawa; 邦昭大澤
Original assignee: Matsushita Electric Industrial Co Ltd
Current assignee: Panasonic Holdings Corp
Priority date: 2001-11-20
Filing date: 2001-11-20
Publication date: 2003-05-30

Abstract

(57)【要約】【課題】実際に即した残響音を生成し違和感の無い残
響音を出力する。【解決手段】音場空間の直接音の初期エネルギの大き
さを周波数帯域毎に求め（ステップＳ１０、Ｓ１１）、
周波数帯域毎の直接音の初期エネルギの大きさと音場空
間の残響時間と残響音減衰式とから残響音を求めるに際
し、前記残響時間を前記周波数帯域毎に求めると共に前
記残響音減衰式から周波数帯域毎の残響音を求め、求め
た周波数帯域毎の残響音を合成して出力する。このよう
に、周波数帯域毎に残響音を求めることで、音場空間の
周波数特性に応じた残響音が生成される。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は残響音生成方法及び
残響音シミュレータに係り、特に、コンサートホールや
教室等の音場空間に応じた残響感、臨場感を実際に即し
て模擬するのに好適な残響音生成方法及び残響音シミュ
レータに関する。

【０００２】

【従来の技術】コンサートホール等の音場空間で発生さ
れた音（直接音）の響きをシミュレートする場合、例え
ば特開平１０―２５４４６９号公報に記載されている様
に、音場空間の音響環境条件に基づいて直接音の初期反
射音を求め、次に、初期反射音の後部残響音を求め、初
期反射音と後部残響音とを接続することで残響音を生成
している。

【０００３】また、例えば特開平９―２５８７６１号公
報記載の従来技術では、オーディオ信号の残響音を人工
的に作り出してオーディオ再生や楽器演奏の臨場感を高
めるために、初期反射音群に続く後部残響音の１つの反
射音から次の反射音までの時間間隔を、音源での音の発
生タイミングを基準とした反射音の発生時刻の自乗に或
る比例定数をもって反比例させ、この比例定数を初期反
射音群の各反射音から次の初期反射音までの時間間隔に
関して得られる比例定数のゆらぎを反映した漸化式写像
に基づいて定めている。

【０００４】

【発明が解決しようとする課題】上述した特開平１０―
２５４４６９号公報記載の従来技術では、後部残響音を
音線数の時間変化率から求めているが、後部残響音の周
波数特性について考慮していないため、残響音に不自然
さが残るという問題がある。

【０００５】また、特開平９―２５８７６１号公報記載
の従来技術も、残響音の周波数特性について考慮してお
らず、更に、この従来技術は、楽器演奏の臨場感などを
高めるために人工的に残響音を付加する技術に過ぎず、
実際の音場空間に即した残響音を生成するものではな
い。従って、生成した残響音をコンサートホール等の音
場空間の設計に利用することができない。

【０００６】本発明の目的は、コンサートホールや教室
等の音場空間に応じた残響感、臨場感を実際に即して模
擬する残響音生成方法及び残響音シミュレータを提供す
ることにある。

【０００７】

【課題を解決するための手段】上記目的を達成する残響
音生成方法は、音場空間の直接音及び初期反射音の初期
エネルギの大きさを周波数帯域毎に求め、前記直接音及
び初期反射音に応じた残響音を残響音減衰式により求め
るに際し前記周波数帯毎に前記残響音を求め、前記周波
数帯域毎の残響音を合成して出力することを特徴とす
る。この構成により、音場空間の周波数特性に応じた残
響音が生成でき、実際に即した音響効果を得ることがで
きる。

【０００８】上記目的を達成する残響音生成方法は、音
場空間の直接音及び初期反射音の初期エネルギの大きさ
と前記音場空間の残響時間と残響音減衰式とから残響音
を求めるに際し、前記残響時間を周波数帯域毎に求め、
前記残響音減衰式から前記周波数帯域毎の残響音を求
め、前記周波数帯域毎の残響音を合成して出力すること
を特徴とする。この構成により、音場空間の周波数特性
に応じた残響音が生成でき、実際に即した音響効果を得
ることができる。

【０００９】上記目的を達成する残響音生成方法は、音
場空間の直接音及び初期反射音の初期エネルギの大きさ
を周波数帯域毎に求め、前記周波数帯域毎の直接音及び
初期反射音の初期エネルギの大きさと前記音場空間の残
響時間と残響音減衰式とから残響音を求めるに際し、前
記残響時間を前記周波数帯域毎に求めると共に前記残響
音減衰式から前記周波数帯域毎の残響音を求め、前記周
波数帯域毎の残響音を合成して出力することを特徴とす
る。この構成により、この構成により、音場空間の周波
数特性に応じた残響音が生成でき、実際に即した音響効
果を得ることができる。

【００１０】上記目的を達成する残響音シミュレータ
は、直接音の音源データベースと、音場空間の音響条件
入力手段と、前記直接音及び初期反射音の周波数帯域毎
の初期エネルギを算出する手段と、前記直接音及び初期
反射音に応じた残響音を前記周波数帯域毎に残響音減衰
式から算出する手段と、前記周波数帯域毎の残響音を合
成する手段と、前記残響音に周波数位相情報を畳み込む
畳み込み演算手段とを備えたことを特徴とする。この構
成により、音場空間の周波数特性に応じた残響音が生成
でき、実際に即した音響効果を得ることができ、音場空
間の設計に利用することが可能となる。

【００１１】上記目的を達成する残響音シミュレータ
は、直接音の音源データベースと、音場空間の音響条件
入力手段と、前記音場空間の前記音響条件に応じた残響
時間を周波数帯域毎に算出する手段と、前記直接音及び
初期反射音の初期エネルギの大きさと前記周波数毎の前
記残響時間と残響音減衰式とから前記周波数帯域毎の残
響音を算出する手段と、前記周波数帯域毎の残響音を合
成する手段と、前記残響音に周波数位相情報を畳み込む
畳み込み演算手段とを備えたことを特徴とする。この構
成により、音場空間の周波数特性に応じた残響音が生成
でき、実際に即した音響効果を得ることができ、音場空
間の設計に利用することが可能となる。

【００１２】上記目的を達成する残響音シミュレータ
は、直接音の音源データベースと、音場空間の音響条件
入力手段と、前記直接音の周波数帯域毎の初期エネルギ
を算出する手段と、前記音場空間の前記音響条件に応じ
た残響時間を前記周波数帯域毎に算出する手段と、前記
周波数帯域毎の初期エネルギと前記周波数帯域毎の残響
時間と残響音減衰式とから前記周波数帯域毎の残響音を
算出する手段と、前記周波数帯域毎の残響音を合成する
手段と、前記残響音に周波数位相情報を畳み込む畳み込
み演算手段とを備えたことを特徴とする。この構成によ
り、音場空間の周波数特性に応じた残響音が生成でき、
実際に即した音響効果を得ることができ、音場空間の設
計に利用することが可能となる。

【００１３】好適には、上記において、前記残響音を算
出する手段は、算出した残響音の時間軸上の配置位置を
乱数によって決めたことを特徴とし、また、前記残響音
を算出する手段は、算出した残響音レベルの大きさを前
記残響音減衰式で定まるレベル以下の任意レベルにする
と共に前記任意レベルを乱数によって決定したことを特
徴とし、更に、前記残響音を算出する手段と、残響音の
到来方向を乱数によって決定したことを特徴とする。こ
の構成により、音場空間で完全に拡散状態となった残響
音が生成可能となる。

【００１４】更に好適には、前記直接音に対する前記音
場空間の初期反射音を幾何音響理論により算出する手段
と、前記合成した残響音と前記初期反射音と前記直接音
とを合成して出力する手段とを備えたことを特徴とす
る。これにより、残響音を可聴化したとき自然な音を出
力することができる。

【００１５】更に好適には、上記において、前記合成し
た残響音のうち、音の拡散状態が所定状態となる残響音
付加スタート時間以後の残響音を前記初期反射音及び前
記直接音に合成して出力することを特徴とする。これに
より、更に実際に即した音を出力することができる。

【００１６】更に好適には、上記において、前記拡散状
態が所定状態となる残響音付加スタート時間のデータを
音場空間毎に格納したデータベースを備えたことを特徴
とする。これにより、様々な音場空間に最適な残響音の
出力が可能となる。

【００１７】

【発明の実施の形態】以下、本発明の一実施形態につい
て、図面を参照して説明する。

【００１８】図１は、本発明の一実施形態に係る残響音
シミュレータの構成図である。この残響音シミュレータ
は、無響室などで録音した楽音ソース等の音源信号デー
タベース１と、この音源信号データベース１から任意に
取り出したアナログの音源信号（以下、直接音とい
う。）をデジタル信号に変換するＡ／Ｄ変換器２と、直
接音から残響音をインパルス応答として算出する（詳細
は後述する）インパルス応答算出部３と、デジタルの直
接音（Ａ／Ｄ変換器２の出力）と残響音（インパルス応
答算出部３の出力）とを畳み込み演算する畳み込み演算
手段４と、この畳み込み演算結果をアナログ信号に変換
するＤ／Ａ変換器５と、このＤ／Ａ変換器５から出力さ
れるアナログ信号にアナログの直接音を混合して増幅す
るアンプ６と、アンプ６の出力を音波として出力するス
ピーカ７とを備えて構成される。

【００１９】この図１に示す実施形態は、残響音をモノ
ラル信号として直接音に付加するものであるが、ステレ
オ信号として付加することも可能である。図２は、残響
音をステレオで付加する残響音シミュレータの構成図で
あり、基本構成は図１と同じである。この残響音シミュ
レータでは、右耳用と左耳用のインパルス応答算出部３
Ｒ、３Ｌを別々に設け、同様に、畳み込み演算手段４
Ｒ、４Ｌ、アンプ６Ｒ、６Ｌ、スピーカ７Ｒ、７Ｌを別
々に設けている。

【００２０】次に、本発明の一実施形態に係る残響音生
成方法の原理を説明する。図３は、幾何音響理論に基づ
く反射音を示す図である。直線音がスピーカから出力さ
れると、その直接音は、音場空間の様々な壁面で反射さ
れる。反射音は、直接音を１回反射した１回反射音と、
１回反射音を更に１回反射した２回反射音と、２回反射
音を更に１回反射した３回反射音と、…に分けることが
できる。何回までの反射音を採用するかは任意である
が、この例では３回反射音まで採用することとする。

【００２１】音場空間の体積をＶ、表面積をＳとした場
合、音波の平均自由行路は４Ｖ／Ｓとして算出される。
直接音発生時点から、１回反射の平均自由行路Ｌ1を音
波が速度ｃで通過する時刻Ｔ1は、４Ｖ／（Ｓｃ）とし
て算出される。１回反射音は、時刻Ｔ1を中心に分布
し、２回反射音は、２回反射の平均自由行路Ｌ2を音波
が通る時刻Ｔ2を中心に分布し、３回反射音は、３回反
射の平均自由行路Ｌ3を音波が通過する時刻Ｔ3を中心に
分布する。

【００２２】１回反射音は、音場空間の様々な壁面から
反射する音であるため、直接音１つに対して複数発生
し、また、１回反射音の各々に対して、夫々複数の２回
反射音が発生するため、反射音の音線数は、図４に示す
様に、指数関数的に増大することになる。

【００２３】１回反射音、２回反射音、３回反射音は、
幾何音響理論に基づいて算出された音であり、各音の全
てが反射の方向を持っている。このため、これらをスピ
ーカから出力しても、真の意味での残響音すなわち後部
残響音としての効果はない。そこで、本実施形態では、
図３に示す時刻Ｔ1から時刻Ｔ2までに付加する残響音を
１回反射音から生成し、時刻Ｔ2から時刻Ｔ3に付加する
残響音を２回反射音から生成し、時刻Ｔ3から時刻Ｔ4に
付加する残響音を３回反射音から生成する。

【００２４】図５は、残響音のエンベロープを示すグラ
フである。残響音の大きさＥは、直接音及び初期反射音
の大きさ（初期エネルギ）をＥ0とした場合、次の残響
音減衰式Ｅ＝Ｅ0・ｅｘｐ（−ｃＡｔ／４Ｖ） …（式１）ここで、ｃ：音速Ａ：音場空間の吸音力（音の周波数によって異なる）Ｖ：音場空間の体積ｔ：残響時間で算出される。

【００２５】残響時刻ｔは、音場空間の表面積Ｓ、体積
Ｖ、壁面の吸音率（壁面の吸音率は音の周波数によって
異なる。）によって周知の算出式から求められる。この
残響時間ｔを式１に代入し、直接音の大きさＥ0を代入
することにより、残響音エンベロープの形状が求められ
る。本実施形態では、残響時間ｔや室内（音場空間）の
吸音力Ａを周波数帯域毎（帯域をどの様に分けるかは任
意である）に求め、また、直接音の大きさＥ0も周波数
帯毎の平均音圧レベルとして算出するため、残響音を求
めるエンベロープの形状は周波数帯域毎に少しずつ変化
することになる。

【００２６】残響音は、式１の残響音エンベロープを満
たす他、その密度（残響音を構成する音場密度）の変化
も考慮する必要がある。残響音も直接音に対するインパ
ルス応答として求めるものであるが、反射音線数と同様
に、時間的に遅れてくる残響音ほど、その残響音も間隔
が密になる。このため、本実施形態では、時刻Ｔ1から
時刻Ｔ2までの間に付加する残響音の音線数は１回反射
音の音線数と同数とし、時刻Ｔ2から時刻Ｔ3までの間に
付加する残響音の音線数は２回反射音の音線数と同数と
し、時刻Ｔ3から時刻Ｔ4までの間に付加する残響音の音
線数は３回反射音の音線数と同数とする。即ち、残響音
の音線数は、本実施形態では、図４のグラフと同様とな
る。

【００２７】更に、個々の残響音が反射音と異なるの
は、発生時刻が不定であることと、レベルの大きさが不
定（勿論、式１のエンベロープは満たす必要がある。）
であることと、到来方向が不定であることにある。そこ
で、本実施形態では、例えば、時刻Ｔ1から時刻Ｔ2まで
の間に付加する残響音の各々の発生時刻を乱数で決める
こととする。また、各残響音の大きさも、式１を満たす
エンベロープ関数以下の乱数で決めた任意レベルとし、
到来方向も乱数で決めた方向とする。

【００２８】到来方向を乱数制御する場合、具体的に
は、ダミーヘッドによって取得した方向別の頭部伝達関
数をデータベース化しておき、乱数によって得られた方
向に対する頭部伝達関数を畳み込み演算することで、方
向の情報を各残響音に付加することができる。

【００２９】時刻Ｔ2から時刻Ｔ3、時刻Ｔ3から時刻Ｔ4
までの残響音についても同様に乱数を用いて発生時刻、
発生レベル、到来方向を決める。このようにして求めた
残響音を時間軸上に並べた音列を図６に示す。

【００３０】図６に示す残響音は、発生時刻、発生レベ
ル、到来方向を乱数で決め、エンベロープを式１で決め
た残響音であるため、完全拡散状態に近い残響音という
ことができる。この残響音を、音場空間に関わらず、常
に時刻Ｔ1から実際に出力すると、音場空間によっては
違和感が生じてしまう。音場空間で音源を設定し過渡応
答を測定した場合、音が完全拡散状態になるまでの時間
は、音場空間の容積や平均吸音力、平均自由行路等によ
って異なってくる。

【００３１】このため、本実施形態では、音場空間毎
に、完全拡散状態を“１．０”としたとき、直接音発生
時点から拡散状態が“０．８”（８０％）になるまでの
時間を予め計測してデータベース化しておき、図７に示
す様に、設計対象とする音場空間の拡散度が８０％とな
る時間（付加スタート時間）をデータベースから読み出
し、この付加スタート時間以後の残響音を、実際に付加
することとする。即ち、この付加スタート時間以後の残
響音を、図３に示す直接音及び幾何音響による初期反射
音の後部残響音としてスピーカから出力する。

【００３２】図８から図１０は、上述した原理に基づい
て本実施形態の残響音シミュレータが実行する処理手順
を示すフローチャートである。

【００３３】先ず、音場空間の音響条件を入力する。即
ち、ステップ１で、室形すなわち音場空間の形状を入力
する。例えば、設計中の音場空間のＣＡＤ情報を入力す
る。次に、ステップＳ２で、壁面材料を入力する。そし
て、周知の残響時間算出式に基づき、周波数帯域毎の残
響時間ＲＴiを算出する（ステップＳ３）。次のステッ
プＳ４では、ステップＳ３で求めた周波数帯域毎の残響
時間ＲＴiや、吸音力Ａi、音場空間の体積Ｖ、表面積Ｓ
などを記憶装置に格納する。

【００３４】次のステップＳ５では、音源信号データベ
ースの中から音源（直接音）を選び、設定する。そし
て、直接音に対する初期反射音（この例では、上述した
１回反射音、２回反射音、３回反射音）のインパルス応
答を、幾何音響理論に基づいて算出し（ステップＳ
６）、このインパルス応答の結果を記憶装置に保存する
（ステップＳ７）。

【００３５】次のステップＳ８では、反射回数に対応し
た音線数の予測関数を算出する。即ち、図４で説明した
反射音線数を求め、次のステップＳ９で、反射回数毎の
音線数Ｎiを記憶装置に保存する。

【００３６】ステップＳ１０（図９）では、周波数帯域
ループ数を設定し、以下のステップＳ２０までの処理を
周波数帯域数ｎだけ繰り返し実行する。この周波数帯域
ループでは、先ず、ステップＳ１１で、初期エネルギＥ
0を算出し、次のステップＳ１２で、図５に示す残響音
エンベロープ関数を生成する。即ち、このステップＳ１
１とステップＳ１２により、残響エンベロープの形状が
決定される。そして、ステップＳ１３で、残響時間と平
均自由行路の関係より反射回数を決定する。

【００３７】次のステップＳ１４では、反射回数のルー
プ処理を開始する。即ち、以下のステップＳ１７までの
処理を、ステップＳ１３で求めた反射回数だけ繰り返
す。この反射回数ループ処理では、ステップＳ１５で、
各反射回数に対応した時間帯域内（図３の時刻Ｔ1から
時刻Ｔ2、時刻Ｔ2から時刻Ｔ3、時刻Ｔ3から時刻Ｔ4）
に残響音が入射するように、乱数制御により、発生時
刻、発生レベル等を決定する。そして、発生時刻、発生
レベル等が決定された残響音を、時間軸上に配置する
（ステップＳ１６）。

【００３８】このステップＳ１５とステップＳ１６の処
理を反射回数だけ繰り返し、ステップＳ１７で反射回数
ループが終了したと判断された場合には、この周波数帯
域ループ処理における残響音の音列（図６）が完成す
る。

【００３９】次のステップＳ１８では、周波数帯域毎に
周波数位相情報を残響音に付加するため、バンドパスフ
ィルタやトーンバースト波形などを畳み込み演算し、演
算結果を記憶装置に保存する（ステップＳ１９）。そし
て、次のステップＳ２０で周波数帯域ループ処理が全帯
域で終了したか否かを判断し、終了していない場合には
ステップＳ１０に戻り、終了したと判断された場合には
ステップＳ２１に進む。

【００４０】ステップＳ２１に進む場合には、図６に示
す残響音の音列が周波数帯域毎に完成されているため、
このステップＳ２１では、周波数帯域毎の図６に示す残
響音の音列を１つの時間軸上に合成する。そして、次の
ステップＳ２２で、図７で説明した付加スタート時間Ｔ
sを算出し、ステップＳ２３で、付加スタート時間Ｔs以
降の残響音の音列を、図３に示す直接音及び初期反射音
のインパルス応答に時間軸上で合成し、この処理を終了
する。

【００４１】以上述べたように、本実施形態による残響
音シミュレータは、音場空間の周波数特性に応じた残響
音を生成するため、実際に即した残響感、臨場感を模擬
することができ、音場設計に利用することが可能とな
る。

【００４２】また、上述した残響音シミュレータは、パ
ーソナルコンピュータ上で実現するが、音場の予測に関
しては高い精度の予測が可能となり、また、可聴化を実
施するためにスピーカから直接音、初期反射音、残響音
を出力した場合、より自然な音が出力される。このた
め、コンサートホール等の音場空間の設計段階でスピー
カの機種選定や、材料の選定などを含めた音響設計ツー
ルとして使用できる。更に、各種効果音用（風呂場、教
室など様々な音場の残響感、臨場感）の生成装置として
適用することもできる。例えば、スタジオ等の何らかの
制作現場において様々なシーンに最適な残響感を持った
音を制作可能となる。

【００４３】

【発明の効果】本発明によれば、音場空間に応じた残響
感、臨場感を実際に即して模擬することができる残響音
シミュレータを提供することができる。

【図面の簡単な説明】

【図１】本発明の一実施形態に係る残響音シミュレータ
の構成図

【図２】本発明の別実施形態に係る残響音シミュレータ
の構成図

【図３】本発明の一実施形態に係る残響音生成方法の原
理説明で用いる初期反射音を示す図

【図４】本発明の一実施形態に係る残響音生成方法の原
理説明で用いる反射音線数を示す図

【図５】本発明の一実施形態に係る残響音生成方法の原
理説明で用いる残響音エンベロープを示す図

【図６】本発明の一実施形態に係る残響音生成方法の原
理説明で用いる残響音の音列を示す図

【図７】本発明の一実施形態に係る残響音生成方法の原
理説明で用いる残響音の付加スタート時間の説明図

【図８】本発明の一実施形態に係る残響音シミュレータ
の処理手順を示すフローチャート

【図９】本発明の一実施形態に係る残響音シミュレータ
の処理手順を示すフローチャート

【図１０】本発明の一実施形態に係る残響音シミュレー
タの処理手順を示すフローチャート

【符号の説明】

１音源信号（直接音）データベース３、３Ｒ、３Ｌインパルス応答算出部４、４Ｒ、４Ｌ畳み込み演算手段

Claims

【特許請求の範囲】

【請求項１】音場空間の直接音及び初期反射音の初期
エネルギの大きさを周波数帯域毎に求め、前記直接音及
び初期反射音に応じた残響音を残響音減衰式により求め
るに際し前記周波数帯毎に前記残響音を求め、前記周波
数帯域毎の残響音を合成して出力することを特徴とする
残響音生成方法。
【請求項２】音場空間の直接音及び初期反射音の初期
エネルギの大きさと前記音場空間の残響時間と残響音減
衰式とから残響音を求めるに際し、前記残響時間を周波
数帯域毎に求め、前記残響音減衰式から前記周波数帯域
毎の残響音を求め、前記周波数帯域毎の残響音を合成し
て出力することを特徴とする残響音生成方法。
【請求項３】音場空間の直接音及び初期反射音の初期
エネルギの大きさを周波数帯域毎に求め、前記周波数帯
域毎の直接音及び初期反射音の初期エネルギの大きさと
前記音場空間の残響時間と残響音減衰式とから残響音を
求めるに際し、前記残響時間を前記周波数帯域毎に求め
ると共に前記残響音減衰式から前記周波数帯域毎の残響
音を求め、前記周波数帯域毎の残響音を合成して出力す
ることを特徴とする残響音生成方法。
【請求項４】直接音の音源データベースと、音場空間
の音響条件入力手段と、前記直接音及び初期反射音の周
波数帯域毎の初期エネルギを算出する手段と、前記直接
音及び初期反射音に応じた残響音を前記周波数帯域毎に
残響音減衰式から算出する手段と、前記周波数帯域毎の
残響音を合成する手段と、前記残響音に周波数位相情報
を畳み込む畳み込み演算手段とを備えたことを特徴とす
る残響音シミュレータ。
【請求項５】直接音の音源データベースと、音場空間
の音響条件入力手段と、前記音場空間の前記音響条件に
応じた残響時間を周波数帯域毎に算出する手段と、前記
直接音及び初期反射音の初期エネルギの大きさと前記周
波数毎の前記残響時間と残響音減衰式とから前記周波数
帯域毎の残響音を算出する手段と、前記周波数帯域毎の
残響音を合成する手段と、前記残響音に周波数位相情報
を畳み込む畳み込み演算手段とを備えたことを特徴とす
る残響音シミュレータ。
【請求項６】直接音の音源データベースと、音場空間
の音響条件入力手段と、前記直接音の周波数帯域毎の初
期エネルギを算出する手段と、前記音場空間の前記音響
条件に応じた残響時間を前記周波数帯域毎に算出する手
段と、前記周波数帯域毎の初期エネルギと前記周波数帯
域毎の残響時間と残響音減衰式とから前記周波数帯域毎
の残響音を算出する手段と、前記周波数帯域毎の残響音
を合成する手段と、前記残響音に周波数位相情報を畳み
込む畳み込み演算手段とを備えたことを特徴とする残響
音シミュレータ。
【請求項７】前記残響音を算出する手段は、算出した
残響音の時間軸上の配置位置を乱数によって決めたこと
を特徴とする請求項４乃至請求項６のいずれかに記載の
残響音シミュレータ。
【請求項８】前記残響音を算出する手段は、算出した
残響音レベルの大きさを前記残響音減衰式で定まるレベ
ル以下の任意レベルにすると共に前記任意レベルを乱数
によって決定したことを特徴とする請求項４乃至請求項
７のいずれかに記載の残響音シミュレータ。
【請求項９】前記残響音を算出する手段と、残響音の
到来方向を乱数によって決定したことを特徴とする請求
項４乃至請求項８のいずれかに記載の残響音シミュレー
タ。
【請求項１０】前記直接音に対する前記音場空間の初
期反射音を幾何音響理論により算出する手段と、前記合
成した残響音と前記初期反射音と前記直接音とを合成し
て出力する手段とを備えたことを特徴とする請求項４乃
至請求項９のいずれかに記載の残響音シミュレータ。
【請求項１１】前記合成した残響音のうち、音の拡散
状態が所定状態となる残響音付加スタート時間以後の残
響音を前記初期反射音及び前記直接音に合成して出力す
ることを特徴とする請求項１０に記載の残響音シミュレ
ータ。
【請求項１２】前記拡散状態が所定状態となる残響音
付加スタート時間のデータを音場空間毎に格納したデー
タベースを備えたことを特徴とする請求項１１記載の残
響音シミュレータ。