(57)【特許請求の範囲】
【請求項１】
（ａ）インパルス形態の音源が加振されたときにマイクロフォンによって得られる音圧信号の入力を受けて室内インパルス応答を検出するステップと、
（ｂ）前記室内インパルス応答が時間領域で複数の区間に分離されるように前記区間の境界を決定するステップと、
（ｃ）前記室内インパルス応答を前記分離された区間別に互いに異なるモデリング技法を適用してモデリングするステップとを含み、
前記複数の区間は、前記室内インパルス応答の時間経過による音響特性を考慮して分離されると共に、
前記複数の区間は、前記音源と前記マイクロフォンとの距離による前記室内インパルス応答の初期遅延を表す第１区間、前記第１区間を除いた区間で前記音源又は前記マイクロフォンの位置変化に対する前記室内インパルス応答の敏感な程度によって区分される第２区間及び第３区間を含む室内インパルス応答モデリング方法。
【請求項２】
前記第２区間は、前記音源又は前記マイクロフォンの位置変化による前記室内インパルス応答の敏感な程度を表す指標が所定値以上の区間であり、前記第３区間は、前記敏感度が所定値未満の区間である請求項１に記載の室内インパルス応答モデリング方法。
【請求項３】
前記指標は、前記音源又は前記マイクロフォンの位置変化による前記室内インパルス応答の変化量と前記位置変化量との比率である請求項２に記載の室内インパルス応答モデリング方法。
【請求項４】
前記（ｃ）ステップは、
前記第１区間を遅延モデルを用いてモデリングするステップと、
前記第２区間を移動平均モデルを用いてモデリングするステップと、
前記第３区間を自己回帰移動平均モデル又は自己回帰モデルを用いてモデリングするステップとを含む請求項２に記載の室内インパルス応答モデリング方法。
【請求項５】
前記（ａ）ステップで、数回にわたって前記室内インパルス応答を検出し、
前記（ｂ）ステップは、現在検出された第１室内インパルス応答に対する前記区間の境界を、前記第１室内インパルス応答以前に検出された第２室内インパルス応答に対して決定された前記区間の境界に基づいて決定する請求項３に記載の室内インパルス応答モデリング方法。
【請求項６】
前記（ｂ）ステップは、
（ｂ１）前記第１室内インパルス応答に対して前記第１区間と前記第２区間との境界である第１境界を決定するステップと、
（ｂ２）前記第１室内インパルス応答に対して前記第２区間と前記第３区間との境界である第２境界を決定するステップとを含む請求項５に記載の室内インパルス応答モデリング方法。
【請求項７】
前記（ｂ１）ステップは、前記第２室内インパルス応答に対して既に決定された前記第１境界に該当する時間における前記第１室内インパルス応答のレベルを所定の閾値と比較し、前記比較結果を用いて前記第１境界を決定する請求項６に記載の室内インパルス応答モデリング方法。
【請求項８】
前記（ｂ２）ステップは、前記第２室内インパルス応答に対してすでに決定された前記第２境界に該当する時間における前記音源又は前記マイクロフォンの位置変化による前記第１室内インパルス応答の敏感な程度を表す指標を所定の閾値と比較し、前記比較結果を用いて前記第２境界を決定する請求項６に記載の室内インパルス応答モデリング方法。
【請求項９】
請求項１乃至８の何れか一項に記載の室内インパルス応答モデリング方法を実行させるためのプログラムを記録したコンピュータで読み取り可能な記録媒体。
【請求項１０】
インパルス形態の音源が加振されたときにマイクロフォンによって得られる音圧信号の入力を受けて室内インパルス応答を検出する室内インパルス応答検出部と、
前記室内インパルス応答が時間領域で複数の区間に分離されるように前記区間の境界を決定する境界決定部と、
前記室内インパルス応答を前記分離された区間別に互いに異なるモデリング技法を適用してモデリングするモデリング部とを備え、
前記複数の区間は、前記室内インパルス応答の時間経過による音響特性を考慮して分離されると共に、
前記複数の区間は、前記音源と前記マイクロフォンとの距離による前記室内インパルス応答の初期遅延を表す第１区間、前記第１区間を除いた区間で前記音源又は前記マイクロフォンの位置変化に対する前記室内インパルス応答の敏感な程度によって区分される第２区間及び第３区間を含む室内インパルス応答モデリング装置。
【請求項１１】
前記第２区間は、前記音源又は前記マイクロフォンの位置変化による前記室内インパルス応答の敏感な程度を表す指標が所定値以上の区間であり、前記第３区間は、前記敏感度が所定値未満の区間である請求項１０に記載の室内インパルス応答モデリング装置。
【請求項１２】
前記指標は、前記音源又は前記マイクロフォンの位置変化による前記室内インパルス応答の変化量と前記位置変化量との比率である請求項１１に記載の室内インパルス応答モデリング装置。
【請求項１３】
前記モデリング部は、
前記第１区間を遅延モデルを用いてモデリングする第１モデリング部と、
前記第２区間を移動平均モデルを用いてモデリングする第２モデリング部と、
前記第３区間を自己回帰移動平均モデル又は自己回帰モデルを用いてモデリングする第３モデリング部とを備える請求項１１に記載の室内インパルス応答モデリング装置。
【請求項１４】
前記室内インパルス応答検出部は、数回にわたって前記室内インパルス応答を検出し、
前記室内インパルス応答検出部は、現在検出された第１室内インパルス応答に対する前記区間の境界を、前記第１室内インパルス応答以前に検出された第２室内インパルス応答に対して決定された前記区間の境界に基づいて決定する請求項１１に記載の室内インパルス応答モデリング装置。
【請求項１５】
前記境界決定部は、前記第１室内インパルス応答に対して前記第１区間と前記第２区間との境界である第１境界を決定しようとするときに、前記第２室内インパルス応答に対して決定された前記第１境界に該当する時間における前記第１室内インパルス応答のレベルを所定の閾値と比較し、前記比較結果を用いて前記第１境界を決定する請求項１４に記載の室内インパルス応答モデリング装置。
【請求項１６】
前記境界決定部は、前記第１室内インパルス応答に対して前記第２区間と前記第３区間との境界である第２境界を決定しようとするときに、前記第２室内インパルス応答に対して決定された前記第２境界に該当する時間における前記音源又は前記マイクロフォンの位置変化による前記第１室内インパルス応答の敏感な程度を表す指標を所定の閾値と比較し、前記比較結果を用いて前記第２境界を決定する請求項１４に記載の室内インパルス応答モデリング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、室内インパルス応答モデリング方法および装置に関し、より詳しくは、正確度を保持しながらも、パラメータの数を減少させることができる室内インパルス応答モデリング方法および装置に関する。
【背景技術】
【０００２】
日常生活で聞く音は、ほとんどすべての場合が残響、すなわち反射音が混ざった音である。自由空間でのように、音源から放射される音だけを聞く場合は、極めて特別な場合を除いては想像することが難しい。室内で音を聴取する場合には、残響の程度に応じて室内空間の大きさ、壁の材質などの空間感を感じるようになる。無響室で録音した原音に人工的に合成された残響を付加して自然かつリアルな音を再生しようとする試みは、以前から放送や音楽産業などに適用されており、最近では３次元音場システム（３Ｄ  ｓｏｕｎｄ  ｆｉｅｌｄ  ｓｙｓｔｅｍ）開発にも適用され、自然な空間的聴取（ｓｐａｔｉａｌ  ｈｅａｒｉｎｇ）を再現しようとする研究が行われている。
【０００３】
残響効果を得るために最も広く用いられる方法は、所望する残響を有する室内で室内インパルス応答（Ｒｏｏｍ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ）を測定し、無響室で録音した原音とコンボリューションする方法である。室内インパルス応答とは、音源がインパルス関数で加振されたときに、聴取者の位置で測定される音圧である。室内インパルス応答をモデリングする技法としては、ＦＩＲ（Ｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ）基盤のオールゼロ（ａｌｌ−ｚｅｒｏ）モデルとＩＩＲ（Ｉｎｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ）基盤のポールゼロ（ｐｏｌｅ−ｚｅｒｏ）モデルなどの多様なモデルがある。一般的に、室内インパルス応答をより正確にモデリングしようとすれば、より多数のパラメータが必要となる。従来には、測定された室内インパルス応答の音響特性を考慮せず、時間領域または周波数領域で１つのモデリング技法を用いて一括的にモデリングしていた。室内インパルス応答をモデリングするパラメータの数が多くなるほど、データ処理に要される時間と必要なメモリやプロセッサなどの容量が増え、より高性能のプロセッサが要求される上に、音響再生をリアルタイムで処理しようとする場合には、多くの制約が生じるようになる。
【発明の開示】
【発明が解決しようとする課題】
【０００４】
本発明が解決しようとする技術的課題は、室内インパルス応答をモデリングすることにおいて、正確度を保持しながらも、パラメータの数を減少させることができる室内インパルス応答モデリング方法および装置を提供することを目的とする。
【課題を解決するための手段】
【０００５】
前記技術的課題を解決するために、本発明に係る室内インパルス応答モデリング方法は、（ａ）インパルス形態の音源が加振されたときにマイクロフォンによって得られる音圧信号の入力を受けて室内インパルス応答を検出するステップと、（ｂ）前記室内インパルス応答が時間領域で複数の区間に分離されるように前記区間の境界を決定するステップと、（ｃ）前記室内インパルス応答を前記分離された区間別に互いに異なるモデリング技法を適用してモデリングするステップとを含むことを特徴とする。
【０００６】
このとき、前記複数の区間は、前記室内インパルス応答の時間経過による音響特性を考慮して分離されることが好ましく、さらに、前記複数の区間は、前記音源と前記マイクロフォンとの距離による前記室内インパルス応答の初期遅延を示す第１区間、前記第１区間を除いた区間で前記音源または前記マイクロフォンの位置変化に対する前記室内インパルス応答の敏感な程度によって区分される第２区間および第３区間を含むことが好ましい。
【０００７】
前記他の技術的課題を解決するために、本発明に係る室内インパルス応答モデリング装置は、インパルス形態の音源が加振されなかったときにマイクロフォンによって得られる音圧信号の入力を受けて室内インパルス応答を検出する室内インパルス応答検出部と、前記室内インパルス応答が時間領域で複数の区間に分離されるように前記区間の境界を決定する境界決定部と、前記室内インパルス応答を前記分離された区間別に互いに異なるモデリング技法を適用してモデリングするモデリング部とを備えることを特徴とする。
【０００８】
前記さらに他の技術的課題を解決するために、前記した室内インパルス応答モデリング方法を実行させるためのプログラムを記録したコンピュータで読み取り可能な記録媒体を提供する。
【発明の効果】
【０００９】
本発明によれば、室内インパルス応答をモデリングすることにおいて、正確度を保持しながらも、パラメータの数を著しく減少させることができる。したがって、音響再生時に室内音響特性を正確に表現しながらも、データ処理に要される時間または必要なメモリやプロセッサなどの容量を減少させることができる。
【発明を実施するための最良の形態】
【００１０】
以下、図面を参照しながら、本発明の好ましい実施形態について詳細に説明する。本発明を説明することにおいて、関連した公知機能あるいは構成に対する具体的な説明が本発明の要旨を不必要かつ不明確にすると判断される場合には、これに対する詳細な説明は省略する。
【００１１】
まず、本発明の好ましい実施形態についての説明に先立ち、室内インパルス応答の時間経過による音響特性を調べることにする。図１は、任意の室内でインパルス音源に対する室内インパルス応答の一例を示す図である。図１に示すように、室内インパルス応答は、音源とマイクロフォン（聴取者）との距離による初期遅延区間１０を有しており、初期反射音が主に表れる区間２０と後期残響音が主に表れる区間３０とに分けることができる。
【００１２】
まず、初期遅延区間１０で、室内インパルス応答は、遅延の大きさだけを表すだけで、音圧の大きさは表さない。したがって、初期遅延区間１０では、室内インパルス応答を簡単な遅延モデルを用いてモデリングすれば、モデリング結果によるパラメータの数を減らすことができるはずである。
【００１３】
次に、初期遅延区間１０以後の室内インパルス応答の音響特性を詳察することにする。図２は、任意の室内で音源を固定してマイクロフォンの位置を変化させながら、室内インパルス応答を測定する様子を説明するための図である。図３は、図２の（ａ）位置と（ｂ）位置それぞれで測定された室内インパルス応答を示している。図３を参照すれば、互いに異なる位置で測定された室内インパルス応答は、異なる形態であることが分かる。これは勿論、音源とマイクロフォンとの距離、音の反射経路などが異なるためである。一方、図２では、音源を固定してマイクロフォンの位置を変化させたが、マイクロフォンを固定して音源の位置を変化させても、互いに異なる形態の室内インパルス応答を得ることができるはずである。
【００１４】
図３に示すように、互いに異なる位置で測定された室内インパルス応答が異なる形態であるため、位置変化によって室内インパルス応答の差がどれくらい表れるかを考慮して見ることができる。したがって、音源またはマイクロフォンの位置変化によって室内インパルス応答の敏感な程度を表す指標を、下記の数式のように定義することにし、これを便宜上、室内インパルス応答の位置敏感度とする。下記の数式は、マイクロフォンの位置変化による室内インパルス応答の変化量と位置変化量との比率を表す。
【００１５】
【数１】
0004829278.tif 000002
ここで、ｎはサンプリング時間、ｈ’［ｎ］は室内インパルス応答の位置敏感度、下付き添字ｌ０は基準位置、Δｌは位置変化量を示す。ｈＮ’［ｎ］は正規化された位置敏感度を示すものであって、％単位で表示される値である。
【００１６】
図４は、前記数式１によって計算された室内インパルス応答の正規化された位置敏感度をサンプリング時間によって示すグラフである。図４を参照すれば、室内インパルス応答の位置敏感度は、サンプリング時間によって、すなわち時間が経過するほどその値が減少することが分かる。図１と結び付けて説明すれば、後期残響音が主に表れる区間３０における位置敏感度は、初期反射音が主に表れる区間２０における位置敏感度よりも相対的に小さい値を有するようになる。
【００１７】
初期反射音とは、インパルス音源から発生した音が壁面や天井のような反射面に反射し、最初に聴取者に伝達された音よりも相対的に長くなった経路を有して伝達される音であって、上述したように、音源またはマイクロフォンの位置変化に相対的に敏感である。後期残響音とは、様々な経路を介して伝達される反射音が重畳したものであって、反射面でエネルギの一部が吸収されて音圧が小さく表れる音であって、上述したように、音源またはマイクロフォンの位置変化に相対的に敏感ではない。
【００１８】
このように、室内インパルス応答の位置敏感度が時間経過によって変化するため、初期遅延区間１０以後の区間を室内インパルス応答の位置敏感度によって区分して適応的にモデリングすれば、モデリング結果によるパラメータの数を減少することができるはずである。したがって、本発明では、室内インパルス応答を音源またはマイクロフォンの位置変化に敏感な区間、例えば図１に示すような初期反射音が主に表れる区間２０と、音源またはマイクロフォンの位置変化にあまり敏感ではない区間、例えば図１に示すような後期残響音が主に表れる区間３０とに分け、前者の区間に対しては、パラメータの数は多いが正確度が高いモデリング方法を用い、後者の区間に対しては、正確度が大きく問題にならないため、パラメータの数が少ないモデリング方法を用いる。
【００１９】
図５は、室内インパルス応答を時間領域で区間に分け、それぞれ周波数特性を分析した結果を示している。図５において、（ａ）は、任意の室内における室内インパルス応答を示しており、（ｂ）は、室内インパルス応答の周波数領域における特性を表した室内伝達関数（Ｒｏｏｍ  Ｔｒａｎｓｆｅｒ  Ｆｕｎｃｔｉｏｎ）を示している。図５の（ｂ）、（ｃ）、および（ｄ）はそれぞれ、（ａ）に示す室内インパルス応答を初期遅延区間、初期反射音区間、および後期残響音区間に分離し、それぞれの周波数領域特性を示すグラフである。このように、室内インパルス応答が時間領域において区間別にそれぞれ周波数特性が分析されるという点を勘案すれば、室内インパルス応答を時間領域で区間に分けてモデリングできることが分かる。
【００２０】
図６は、本発明の一実施形態に係る室内インパルス応答モデリング装置を示すブロック図である。図６を参照すれば、室内インパルス応答モデリング装置は、マイクロフォン６０２と、室内インパルス応答検出部６１０と、境界決定部６２０と、モデリング部６３０と、合成部６４０とを備えてなされる。本実施例に係る室内インパルス応答モデリング装置の動作を、本発明の一実施形態に係る室内インパルス応答モデリング方法を示すフローチャートである図７と結び付けて説明する。
【００２１】
マイクロフォン６０２は、任意の室内６００でインパルス形態の音源６０１が加振されたときに時音圧信号を得て、室内インパルス応答検出部６１０は、マイクロフォン６０２によって得られる音圧信号の入力を受けて室内インパルス応答を検出する（Ｓ７１０）。
【００２２】
境界決定部６２０は、室内インパルス応答検出部６１０で検出された室内インパルス応答が時間領域で時間経過による音響特性が考慮された複数の区間に分離されるように区間の境界を決定する（Ｓ７２０）。一実施形態において、境界決定部６２０は、音源６０１とマイクロフォン６０２との距離による室内インパルス応答の初期遅延を表す第１区間、または第１区間を除いた残りの区間ですでに説明したような室内インパルス応答の位置敏感度によって区分される第２区間および第３区間に分離されるように区間の境界を決定する。ここで、第１区間は、図１に示す初期遅延区間１０と同じである。
【００２３】
以下、便宜上、第１区間と第２区間との境界を第１境界、第２区間と第３区間との境界を第２境界として表現する。
境界決定部６２０は、第１境界を決定するときに、例えば室内インパルス応答を時間によってスキャンしながら、レベルが所定の大きさ以上になるか否かを検査して決定することができる。また、境界決定部６２０は、第２境界を決定するときに、室内インパルス応答の位置敏感度が所定の値以上である区間が第２区間となり、位置敏感度が前記所定の値以下である区間が第３区間となるように第２区間と第３区間との境界を決定することができる。再び図４を参照すれば、室内インパルス応答の位置敏感度が所定の閾値、例えば１０％となる地点に該当する時間は、約０．４２ｓｅｃであることが分かる。したがって、前記所定の値を１０％としたときに、境界決定部６２０は第２境界を０．４２ｓｅｃ地点として決定する。
【００２４】
境界決定部６２０によって室内インパルス応答の区間の境界が決定されれば、モデリング部６３０は、各区間のサンプルに対して区間別に互いに異なるモデリング技法を適用してモデリングする（Ｓ７３０）。
【００２５】
一実施形態において、モデリング部６３０は、図６に示すように、第１モデリング部６３１と、第２モデリング部６３２と、第３モデリング部６３３とを備えてなされる。この場合に、モデリング部６３０の動作を、前記ステップ７３０のより具体的な一実施形態を示すフローチャートである図８と結び付けて説明する。
【００２６】
第１モデリング部６３１は、境界決定部６２０から室内インパルス応答の第１区間のサンプルの入力を受け、例えば下記の数式のような形態の遅延モデルを用いてモデリングする（Ｓ８１０）。
【００２７】
【数２】
0004829278.tif 000003
ここで、Ｄは、第１区間に含まれたサンプリング数を意味する。
【００２８】
第２モデリング部６３２は、境界決定部６２０から室内インパルス応答の第２区間のサンプルの入力を受け、例えば下記の数式のような形態のＭＡ（Ｍｏｖｉｎｇ  Ａｖｅｒａｇｅ）モデル（Ａｌｌ−ｚｅｒｏモデル）を用いてモデリングする（Ｓ８２０）。
【００２９】
【数３】
0004829278.tif 000004
【００３０】
第３モデリング部６３３は、境界決定部６２０から室内インパルス応答の第３区間のサンプルの入力を受け、例えば下記の数式のような形態のＡＲＭＡ（Ａｕｔｏ−ｒｅｇｒｅｓｓｉｖｅ  Ｍｏｖｉｎｇ  Ａｖｅｒａｇｅ）モデル（ｐｏｌｅ−ｚｅｒｏモデル）またはＡＲ（Ａｕｔｏ−ｒｅｇｒｅｓｓｉｖｅ）モデル（Ａｌｌ−ｐｏｌｅモデル）を用いてモデリングする（Ｓ８３０）。
【００３１】
【数４】
0004829278.tif 000005
【００３２】
【数５】
0004829278.tif 000006
【００３３】
再び図７を参照すれば、最後に、合成部６４０は、第１〜第３モデリング部６３１、６３２、６３３でモデリング結果、生成されたパラメータを合成して最終モデリングとなった室内インパルス応答のパラメータを出力する（Ｓ７４０）。
【００３４】
図９は、本実施形態によって室内インパルス応答が区間別にモデリングされる形態を示す図である。このとき、第１区間の伝達関数は時間遅延フィルタで、第２区間の伝達関数はＦＩＲ（Ｆｉｎｉｔｅ  ｉｍｐｌｕｓｅ  ｒｅｓｐｏｎｓｅ）フィルタで、第３区間の伝達関数はＩＩＲ（Ｉｎｆｉｎｉｔｅ  ｉｍｐｕｌｓｅ  ｒｅｓｐｏｎｓｅ）フィルタで実現される。
【００３５】
Ａｌｌ−ｚｅｒｏモデルは、室内インパルス応答のモデリングに主に用いられる技法であって、ｐｏｌｅ−ｚｅｒｏモデルまたはａｌｌ−ｐｏｌｅモデルよりも正確にモデリングすることができるが、必要なパラメータの数が多い。本実施形態では、従来のように、室内インパルス応答全体をＡｌｌ−ｚｅｒｏモデルを用いてモデリングせずに、室内インパルス応答の位置敏感度が相対的に敏感な第２区間に対してのみＡｌｌ−ｚｅｒｏモデルを用いてモデリングする。さらに、室内インパルス応答の位置敏感度が相対的にあまり敏感でない第３区間は、正確度が大きく問題にならないため、必要なパラメータの数が少ないｐｏｌｅ−ｚｅｒｏモデルまたはａｌｌ−ｐｏｌｅモデルを用いてモデリングすることによって、パラメータの数を著しく減らすことができる。
【００３６】
一方、本実施形態によってモデリングされた室内インパルス応答は、ブロックコンボリューションを用いて原音と合成される。この場合に、室内インパルス応答を長さが異なる３つのブロック、すなわち第１区間に該当するブロック、第２区間に該当するブロック、および第３区間に該当するブロックに分けて原音とブロックコンボリューションを実効し、音響を再生することができる。
【００３７】
室内インパルス応答検出部６１０は、室内構造、音源、または聴取者の位置が変化する環境で連続的に数回にわたって室内インパルス応答を検出するように実現される。以下、このような場合に、境界決定部６２０が区間の境界をより効率的に決定することができる方法について説明する。図１０は、室内インパルス応答検出部６１０が連続的に数回にわたって室内インパルス応答を検出するときに、図７に示すステップ７２０の具体的な一実施形態を示すフローチャートである。本実施例によれば、境界決定部６２０は、現在検出された室内インパルス応答に対する区間の境界を、以前に検出された室内インパルス応答に対してすでに決定された区間の境界に基づいて決定する。図１０を参照すれば、境界決定部６２０は、以前に検出された室内インパルス応答の第１境界に基づいて現在検出された室内インパルス応答の第１境界を決定し（Ｓ１０１０）、以前に検出された室内インパルス応答の第２境界に基づいて現在検出された室内インパルス応答の第２境界を決定する（Ｓ１０２０）。
【００３８】
図１１は、図１０に示すステップ１０１０のより具体的な一実施形態を示すフローチャートである。
境界決定部６２０は、以前に検出された室内インパルス応答に対してすでに決定された第１境界に該当する時間ｎの入力を受け（Ｓ１１１０）、現在検出された室内インパルス応答の入力を受ける（Ｓ１１２０）。
【００３９】
次に、時間ｎにおける現在検出された室内インパルス応答のレベルｈ［ｎ］が所定の閾値である第１閾値を越えるか否かを判断し（Ｓ１１２５）、超える場合にはステップ１１３０に進み、超えない場合にはステップ１１４０に進む。ここで、第１境界は、信号が存在しない部分と存在する部分との間の境界であるため、第１境界を正確に判断するためには、室内インパルス応答の最高レベルよりも相当に小さい値を用いることが好ましい。例えば、最高レベルが１である場合には、０．００５を第１閾値として用いることができる。
【００４０】
ｈ［ｎ］が前記第１閾値を超える場合は、ｎが現在検出された室内インパルス応答の第１境界よりも右側に位置するという意味である。したがって、ステップ１１３０でｎを１だけ減少させ、ステップ１１３５で再びｈ［ｎ］が第１閾値を超えるか否かを判断し、超える場合にはステップ１１３０に再び戻り、超えなければ、そのときのｎを現在検出された室内インパルス応答の第１境界として決定する（Ｓ１１５０）。
【００４１】
ｈ［ｎ］が前記第１閾値を超えない場合は、ｎが現在検出された室内インパルス応答の第１境界よりも左側に位置するという意味である。したがって、ステップ１１４０でｎを１だけ増加させ、ステップ１１４５でｈ［ｎ］が第１閾値よりも小さいか否かを判断し、小さい場合にはステップ１１４０に再び戻り、小さくなければ、そのときのｎを現在検出された室内インパルス応答の第１境界として決定する（Ｓ１１５０）。
【００４２】
図１２は、図１０に示すステップ１０２０のより具体的な一実施形態を示すフローチャートである。
境界決定部６２０は、以前に検出された室内インパルス応答に対してすでに決定された第２境界に該当する時間ｎの入力を受け（Ｓ１２１０）、現在検出された室内インパルス応答の入力を受ける（Ｓ１２２０）。
【００４３】
次に、時間ｎにおける現在検出された室内インパルス応答の位置敏感度ｈ’［ｎ］が所定の閾値である第２閾値を超えるか否かを判断し（Ｓ１２２５）、超える場合にはステップ１２４０に進み、超えない場合にはステップ１２３０に進む。ここで、第２閾値は勿論、第２区間と第３区間とを区分するための基準になる室内インパルス応答の位置敏感度である。
【００４４】
ｈ’［ｎ］が前記第２閾値を超える場合は、ｎが現在検出された室内インパルス応答の第２境界よりも左側に位置するという意味である。したがって、ステップ１２４０でｎを１だけ増加させ、ステップ１２４５で再びｈ’［ｎ］が第２閾値を超えるか否かを判断し、超える場合にはステップ１２４０に再び戻り、超えなければ、そのときのｎを現在検出された室内インパルス応答の第２境界として決定する（Ｓ１２５０）。
【００４５】
ｈ’［ｎ］が前記第２閾値を超えない場合には、ｎが現在検出された室内インパルス応答の第２境界よりも右側に位置するという意味である。したがって、ステップ１２３０でｎを１だけ減少させ、ステップ１２３５でｈ’［ｎ］が第２閾値よりも小さいか否かを判断し、小さい場合にはステップ１２３０に再び戻り、小さくなければ、そのときのｎを現在検出された室内インパルス応答の第２境界として決定する（Ｓ１２５０）。
【００４６】
このように、連続的に数回にわたって室内インパルス応答が検出される場合に、区間の境界を決定するときに、以前に検出された室内インパルス応答に対してすでに決定された区間の境界を用いることによって、区間の境界をより迅速かつ効率的に決定することができる。
【００４７】
上述した本発明によれば、室内インパルス応答のモデリング性能において、相当な正確度を保持しながらも、モデリングするためのパラメータの数を著しく減少させることができる。したがって、音響再生時に室内音響特性を正確に表現しながらも、データ処理に要される時間または必要なメモリやプロセッサなどの容量を減少できるようになる。さらに、上述した第２区間と第３区間との境界を相違して設定することによって、モデリング性能とパラメータ数とのトレードオフ関係を適切に調節することができる。例えば、第２区間と第３区間との境界の基準となる位置敏感度の閾値を低く設定することによって、パラメータの数を増やす代わりに、より正確にモデリングすることができる。
【００４８】
一方、上述した本発明の実施形態は、コンピュータで実行されるプログラムで作成が可能であり、コンピュータで読み取り可能な記録媒体を用いて前記プログラムを動作させる汎用デジタルコンピュータで実現することができる。前記コンピュータで読み取り可能な記録媒体は、磁気格納媒体（例えば、ＲＯＭ、フロッピーディスク、ハードディスクなど）、光学的読み取り媒体（例えば、ＣＤ−ＲＯＭ、ＤＶＤなど）、および搬送波（例えば、インターネットを介した送信）のような格納媒体を含む。
【００４９】
これまで、本発明に対して好ましい実施形態を中心として詳察した。本発明が属する技術分野において通常の知識を有する者は、本発明が本発明の本質的な特性から逸脱しない範囲内で変形した形態で具現されることを理解することができるであろう。したがって、開示された実施形態は、限定的な観点ではなく説明的な観点で考慮されなければならない。本発明の範囲は、上述した説明ではなく特許請求の範囲に示されており、これと同等な範囲内に存在するすべての差異点は、本発明に含まれたものとして解釈されるべきである。
【図面の簡単な説明】
【００５０】
【図１】任意の室内でインパルス音源に対する室内インパルス応答の一例を示す図である。
【図２】任意の室内で音源を固定してマイクロフォンの位置を変化させながら、室内インパルス応答を測定する様子を説明するための図である。
【図３】図２の（ａ）位置と（ｂ）位置それぞれで測定された室内インパルス応答を示す図である。
【図４】前記数式１によって計算された室内インパルス応答の位置敏感度をサンプリング時間に応じて示すグラフである。
【図５】室内インパルス応答を時間領域で区間に分け、それぞれ周波数特性を分析した結果を示す図である。
【図６】本発明の一実施形態に係る室内インパルス応答モデリング装置を示すブロック図である。
【図７】本発明の一実施形態に係る室内インパルス応答モデリング方法を示すフローチャートである。
【図８】図７に示すステップ７３０のより具体的な一実施形態を示すフローチャートである。
【図９】本発明の一実施形態によって室内インパルス応答が区間別にモデリングされる形態を表現した図である。
【図１０】連続的に数回にわたって室内インパルス応答が検出されるときに、図７に示すステップ７２０の具体的な一実施形態を示すフローチャートである。
【図１１】図１０に示すステップ１０１０のより具体的な一実施形態を示すフローチャートである。
【図１２】図１０に示すステップ１０２０のより具体的な一実施形態を示すフローチャートである。
【符号の説明】
【００５１】
６００：任意の室内
６０１：音源
６０２：マイクロフォン
６１０：室内インパルス応答検出部
６２０：境界決定部
６３０：モデリング部
６３１：第１モデリング部
６３２：第２モデリング部
６３３：第３モデリング部
６４０：合成部
【図１】
0004829278.tif 000007
【図２】
0004829278.tif 000008
【図３】
0004829278.tif 000009
【図４】
0004829278.tif 000010
【図５】
0004829278.tif 000011
【図６】
0004829278.tif 000012
【図７】
0004829278.tif 000013
【図８】
0004829278.tif 000014
【図９】
0004829278.tif 000015
【図１０】
0004829278.tif 000016
【図１１】
0004829278.tif 000017
【図１２】
0004829278.tif 000018