(57)【特許請求の範囲】
【請求項１】
同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出されたインパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源それぞれの到来時刻及び到来方向を推定する仮想音源情報推定部と、
前記仮想音源情報推定部により推定された前記仮想音源ごとに、複数の前記マイクロホンのうち前記任意の位置からの方向が前記仮想音源の前記到来方向に対応したマイクロホンを選択し、前記仮想音源の前記到来時刻と、当該仮想音源について選択された前記マイクロホンにより収音した前記信号から算出された前記インパルス応答における局所的なピークそれぞれが出現する時刻との差分に基づいて、当該インパルス応答における前記局所的なピークのうち、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音のピークを特定する特定処理を行い、前記特定処理により特定された前記インパルス応答における前記ピークを、当該インパルス応答の算出の元となった前記信号が得られた前記マイクロホンに対応する方向から到来する反射音と判定する比較判定部と、
前記インパルス応答から、当該インパルス応答における前記ピークのうち前記比較判定部により反射音と判定された前記ピーク以外のピークを除去して方向別のインパルス応答を生成する反射音除去部と、
を備えることを特徴とするインパルス応答生成装置。
【請求項２】
前記比較判定部は、前記仮想音源ごとの前記特定処理において、複数の前記マイクロホンのうち前記任意の位置からの方向が前記仮想音源の前記到来方向に対応したマイクロホンを選択し、選択した前記マイクロホンにより収音した信号から算出された前記インパルス応答における局所的なピークのうち、選択した前記マイクロホンと前記任意の位置との距離に応じた遅延を前記仮想音源の前記到来時刻から減算した時刻を含んだ所定の時刻範囲内で生じた前記ピークを、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音のピークと特定する、
ことを特徴とする請求項１に記載のインパルス応答生成装置。
【請求項３】
前記仮想音源情報推定部は、前記任意の位置からの方向が異なる複数の前記マイクロホンそれぞれにより収音された前記信号に基づき算出された前記インパルス応答と、前記任意の位置に設置された前記マイクロホンにより前記音響信号を収音して得られた信号に基づき算出されたインパルス応答とを用いて算出した瞬時インテンシティに基づいて前記仮想音源の前記到来時刻及び前記到来方向を得る、又は、前記任意の位置からの方向が異なる複数の前記マイクロホンそれぞれにより収音された前記信号に基づき算出された前記インパルス応答と、前記任意の位置に設置された前記マイクロホンにより前記音響信号を収音して得られた信号に基づき算出されたインパルス応答との相関に基づいて前記仮想音源の前記到来時刻及び前記到来方向を得る、
ことを特徴とする請求項１又は請求項２に記載のインパルス応答生成装置。
【請求項４】
コンピュータを、
同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出されたインパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源それぞれの到来時刻及び到来方向を推定する仮想音源情報推定手段と、
前記仮想音源情報推定手段により推定された前記仮想音源ごとに、複数の前記マイクロホンのうち前記任意の位置からの方向が前記仮想音源の前記到来方向に対応したマイクロホンを選択し、前記仮想音源の前記到来時刻と、当該仮想音源について選択された前記マイクロホンにより収音した前記信号から算出された前記インパルス応答における局所的なピークそれぞれが出現する時刻との差分に基づいて、当該インパルス応答における前記局所的なピークのうち、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音のピークを特定する特定処理を行い、前記特定処理により特定された前記インパルス応答における前記ピークを、当該インパルス応答の算出の元となった前記信号が得られた前記マイクロホンに対応する方向から到来する反射音と判定する比較判定手段と、
前記インパルス応答から、前記比較判定手段により当該インパルス応答において反射音と判定された前記ピーク以外の前記ピークを除去して方向別のインパルス応答を生成する反射音除去手段と、
を具備するインパルス応答生成装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、インパルス応答生成装置及びプログラムに関する。
【背景技術】
【０００２】
テレビ番組や音楽、映画などの音響制作において、音の広がり感や臨場感を与えるために、残響を付加する技術が用いられている。空間の残響をシミュレーションで模擬したＩＩＲ（Infinite impulse response、無限インパルス応答）フィルタ型のものや、実空間で測定したインパルス応答を用いて残響を生成するＦＩＲ（Finite impulse response、有限インパルス応答）フィルタ型のものなど、様々な残響付加装置が開発されている（例えば、特許文献１参照）。特に、ＦＩＲフィルタ型の残響付加装置では、インパルス応答を測定した空間の残響特性が忠実に再現できるという特長がある。
【０００３】
ところで、５．１チャネルサラウンド方式やスーパーハイビジョン用の２２．２マルチチャンネル方式といったマルチチャンネル音響方式は、平面上もしくは空間内に配置した複数のスピーカを用いて、様々な方向から到来する音を表現することができるという特長を持っている（例えば、非特許文献１参照）。一般的にこれらのコンテンツを制作する過程でも、残響付加装置を用いて豊かな臨場感と空間的な広がり感を付与することが多い。実空間において壁や床、天井など様々な方向から反射音が到来するように、マルチチャンネル音響方式における残響付加でも、ある任意の方向から到来する特徴的な反射音（方向別の残響）を、その方向に設置したスピーカから再生することで、インパルス応答を測定した空間の残響特性をできるだけ忠実に再現するのが理想である。特に直接音の後、１００ｍｓ（ミリ秒）くらいまでの間に到来する初期反射音は、音像の幅や方向感などの知覚に寄与することが分かっており、残響付加において初期反射音を方向別に再生することは、インパルス応答を測定した空間の印象を再現する上で重要である。そこで、指向性マイクロホンで測定したインパルス応答を用いて、方向別のインパルス応答を生成する技術が提案されてきた（例えば、特許文献２、３参照）。
【０００４】
一方、近接四点法（例えば、非特許文献２参照）やそれを拡張した方法（例えば、非特許文献３、４参照）など、複数のマイクロホンを用いて、音源や仮想音源（反射音）の位置を推定する様々な技術が提案されてきた。これらの手法では、各マイクロホンの相互相関やインテンシティ算出することにより、受音点における反射音の到来時刻、到来方向を予測することができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１４−４５２８２号公報
【特許文献２】特開２０１３−２３８６４３号公報
【特許文献３】特開２０１４−１１２７６７号公報
【非特許文献】
【０００６】
【非特許文献１】W.woszczyk，et al.，“Space Builder: An impulse response-based tool for immersive 22.2 channel ambiance design”,日本，AES 40th international conference，2010年
【非特許文献２】山崎芳男，「近接4点法による音場の空間情報の可視化」，テレビジョン学会誌，1990年，Vol.44，No.3，p.253-258
【非特許文献３】西村昌浩,外2名,「BoSCマイクロホンによる仮想音源推定のための近接４点法の拡張アルゴリズムの提案」，2012年3月，日本音響学会講演論文集、p.789-790
【非特許文献４】小林万理恵他，外２名，「分散配置した近接４点法マイクロホンによる音源位置推定と音源分離」，日本音響学会講演論文集，2012年3月，p.685-686
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上述した技術によって生成したインパルス応答には、厳密にはある任意の方向以外から到来する反射音も含まれている。指向性マイクロホンとは言っても、その指向性には幅があるため、近接する複数のマイクロホンで同一の反射音を拾ってしまうからである。このため、これらのインパルス応答を用いて残響付加を行う場合、ある反射音が複数のスピーカから再生され、リスニングポイントでの反射音の密度が高くなってしまい、測定音場の空間印象が厳密には再現されないという問題点がある。
【０００８】
本発明は、このような事情を考慮してなされたもので、残響付加において、測定音場の空間印象をより厳密に再現できる方向別のインパルス応答を生成するインパルス応答生成装置及びプログラムを提供する。
【課題を解決するための手段】
【０００９】
本発明の一態様は、同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出された各インパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源の情報を推定する仮想音源情報推定部と、前記仮想音源情報推定部により推定された前記仮想音源の情報と、前記インパルス応答との比較に基づいて、前記インパルス応答において、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音を判定する比較判定部と、前記インパルス応答から、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来すると前記比較判定部により判定された反射音以外の反射音を間引いて方向別のインパルス応答を生成する反射音除去部と、を備えることを特徴とするインパルス応答生成装置である。
【００１０】
本発明の一態様は、上述するインパルス応答生成装置であって、前記仮想音源の情報は、前記任意の位置における前記仮想音源の到来時刻及び到来方向を示し、前記比較判定部は、前記仮想音源の到来時刻と、前記仮想音源の到来方向に対応したマイクロホンにより収音した信号から算出された前記インパルス応答の波形構造が表す経時的な音の大きさの変化との比較に基づいて、前記インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音を判定する、ことを特徴とする。
【００１１】
本発明の一態様は、上述するインパルス応答生成装置であって、前記比較判定部は、複数の前記マイクロホンのうち前記仮想音源の到来方向に対応したマイクロホンを選択し、選択した前記マイクロホンにより収音した信号から算出された前記インパルス応答において、選択した前記マイクロホンと前記任意の位置との距離に応じた遅延を前記仮想音源の到来時刻から減算した時刻を含んだ所定の時刻範囲内で生じた局所的なピークを特定する処理を前記仮想音源のそれぞれについて行い、前記インパルス応答において前記処理により特定された局所的なピークを、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音と判定する、ことを特徴とする。
【００１２】
本発明の一態様は、コンピュータを、同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出された各インパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源の情報を推定する仮想音源情報推定手段と、前記仮想音源情報推定手段により推定された前記仮想音源の情報と、前記インパルス応答との比較に基づいて、前記インパルス応答において、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音を判定する比較判定手段と、前記インパルス応答から、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来すると前記比較判定手段により判定された反射音以外の反射音を間引いて方向別のインパルス応答を生成する反射音除去手段と、を具備するインパルス応答生成装置として機能させるためのプログラムである。
【発明の効果】
【００１３】
本発明によれば、残響付加において、測定音場の空間印象をより厳密に再現できる方向別のインパルス応答を生成することが可能となる。
【図面の簡単な説明】
【００１４】
【図１】単一指向性マイクロホンの指向特性を示す図である。
【図２】単一指向性マイクロホンにより測定されたインパルス応答を示す図である。
【図３】本発明の一実施形態によるインパルス応答生成装置の構成を示す機能ブロック図である。
【図４】同実施形態によるインパルス応答測定方法の例を示す図である。
【図５】同実施形態によるインパルス応答生成装置における方向別インパルス応答生成処理のフロー図である。
【図６】同実施形態によるインパルス応答測定方法の例を示す図である。
【図７】同実施形態によるインテンシティを用いた仮想音源情報の推定処理を説明するための図である。
【図８】同実施形態による相互相関係数を用いた仮想音源情報推定処理のフロー図である。
【図９】同実施形態による比較判定処理及び反射音除去処理の例を示す図である。
【発明を実施するための形態】
【００１５】
以下、図面を参照しながら本発明の実施形態を詳細に説明する。
まず、図１及び図２を用いて、複数の単一指向性マイクロホンによる収音を説明する。
図１は、単一指向性マイクロホンの指向特性を示す図である。マイクロホンＭａ及びＭｂは、単一指向性マイクロホンである。同図では、方向別の反射音を収音するため、マイクロホンＭａをＤ１方向に向けて、マイクロホンＭｂをＤ２方向に向けて配置している。マイクロホンＭａの指向特性はＲａであり、マイクロホンＭｂの指向特性はＲｂである。点線は、マイクロホンＭａ、Ｍｂが全指向性マイクロホンであった場合の収音範囲を示す。Ｄ１方向からは反射音Ａが、Ｄ２方向からは反射音Ｂが到来する。
【００１６】
図２は、図１に示す単一指向性マイクロホンにより測定されたインパルス応答を示す図である。インパルス応答は、ＴＳＰ（Time Stretched Pulse）信号などの音響測定用の信号を放射し、マイクロホンにより収音（測定）された信号を処理することによって算出されるものであり、音の大きさの経時的な変化を表す。本明細書では、あるマイクロホンにより収音された信号を処理することによって算出されたインパルス応答を、そのマイクロホンにより測定されたインパルス応答と記載する。具体的には、ある任意の位置のマイクロホンにより測定されたインパルス応答は、時刻ｔと、時刻ｔにおいてそのマイクロホンが収音した音の振幅とを対応付けた情報として得られる。図２（ａ）は、図１に示すマイクロホンＭａにより測定されたインパルス応答ｈａ（ｔ）を示し、図２（ｂ）は、図１に示すマイクロホンＭｂにより測定されたインパルス応答ｈｂ（ｔ）を示す。
【００１７】
本来、Ｄ１方向に向けて設置したマイクロホンＭａのみでＤ１方向からの反射音Ａを収音し、Ｄ２方向に向けて設置したマイクロホンＭｂのみでＤ２方向からの反射音Ｂを収音するのが理想である。しかし、図１に示すように、実際には指向性マイクロホンの指向特性には幅がある。そのため、反射音Ａ、ＢともマイクロホンＭａ、Ｍｂの両方で収音され、図２に示すインパルス応答ｈａ（ｔ）、ｈｂ（ｔ）のような波形となってしまう。つまり、インパルス応答ｈａ（ｔ）では、反射音Ａより反射音Ｂの方がレベルは小さく、インパルス応答ｈｂ（ｔ）では、反射音Ｂより反射音Ａの方がレベルは小さいものの、それぞれのマイクロホンに対応する方向（測定したい方向）以外から到来する反射音も含まれてしまう。
【００１８】
そのため、複数の指向性マイクロホンのそれぞれにより測定されたインパルス応答を使って残響付加を行うと、各マイクロホンに対応する方向（測定したい方向）以外から到来する反射音が複数のスピーカから再生されてしまい、測定空間の印象が変わってしまう可能性がある。そこで、本実施形態のインパルス応答生成装置は、マイクロホンＭａにより収音された反射音Ｂや、マイクロホンＭｂにより収音された反射音Ａのような、マイクロホンに対応する方向（測定したい方向）以外から到来する反射音を、測定されたインパルス応答から除去し、より正確な方向別のインパルス応答を生成する。
【００１９】
方向別のインパルス応答を生成するため、本実施形態では、まず、複数のマイクロホンを空間内の任意の位置に配置し、インパルス応答を測定する。インパルス応答生成装置は、各マイクロホンにより測定されたインパルス応答を用いて、マイクロホン群に囲まれた任意の位置における、空間内の仮想音源（反射音）の到来方向及び到来時刻を推定する。インパルス応答生成装置は、推定したそれらの情報と、各マイクロホンにより測定されたインパルス応答に含まれる反射音の到来時刻とを照らし合わせる。インパルス応答生成装置は、照らし合わせの結果に基づき、各マイクロホンに対応する方向（測定したい方向）以外から到来する反射音を、そのマイクロホンにより測定されたインパルス応答からの除去対象として特定する。インパルス応答生成装置は、各マイクロホンにより測定されたインパルス応答から、除去対象の任意の反射音を除去し、方向別のインパルス応答を生成する。なお、各マイクロホンにより収音されるべき到来方向は、そのマイクロホンの設置方向に加え、隣接するマイクロホンの設置方向等によっても異なり得る。
【００２０】
図３は、本発明の一実施形態によるインパルス応答生成装置１の構成を示す機能ブロック図であり、本実施形態と関係する機能ブロックのみを抽出して示してある。インパルス応答生成装置１は、例えば、コンピュータ装置により実現することができる。同図に示すように、インパルス応答生成装置１は、記憶部１１と、測定情報取得部１２と、仮想音源情報推定部１３と、比較判定部１４と、反射音除去部１５とを備えて構成される。
【００２１】
記憶部１１は、各マイクロホンにより測定されたインパルス応答など各種情報を記憶する。
測定情報取得部１２は、空間内の任意の位置に配置した複数のマイクロホンのそれぞれにより測定されたインパルス応答を取得し、記憶部１１に書き込む。以下では、ｎ本のマイクロホンのうちｉ番目のマイクロホンをＭｉ（ｉは１以上ｎ以下の整数）と記載し、マイクロホンＭｉにより測定されたインパルス応答をｈｉ（ｔ）と記載する。また、インパルス応答ｈｉ（ｔ）に含まれる局所的なピーク（包絡の極大点）における時刻を、時刻順にｔｉ１，ｔｉ２，…とする。
【００２２】
仮想音源情報推定部１３は、各マイクロホンにより測定されたインパルス応答からインテンシティや相互相関係数などを算出することにより、マイクロホン群（マイクロホンＭ１〜Ｍｎ）で囲まれた任意の点における、空間内の仮想音源の情報を推定する。推定の対象は、方向知覚への寄与が大きいと考えられるレベルの仮想音源である。方向知覚への寄与が大きいと考えられるレベルの仮想音源は、主に初期反射音部分、直接音から１００ｍｓ（ミリ秒）くらいまでの間に受音点に到来する仮想音源である。
比較判定部１４は、仮想音源情報推定部１３が推定した仮想音源の情報と元のインパルス応答の波形構造とを比較する。比較判定部１４は、残響付加の際に、ある反射音が複数のスピーカから重複して再生されないよう、比較結果に基づき、元のインパルス応答から除去すべき反射音（局所的なピーク）を判定する。
反射音除去部１５は、元のインパルス応答から、比較判定部１４が判定した反射音（局所的なピーク）を除去し、方向別のインパルス応答を生成する。この生成された方向別のインパルス応答を用いて、従来技術により残響付加が行われる。
【００２３】
図４は、インパルス応答測定方法の例を示す図である。インパルス応答の測定空間内に、インパルス応答測定用信号を再生するスピーカと複数のマイクロホンＭ０〜Ｍｎとを配置する。ここでは、マイクロホンＭ０〜Ｍｎは、単一指向性マイクロホンであるが、全指向性マイクロホンを使用してもよい。スピーカの位置や、マイクロホンの個数及び位置などは、測定空間の広さや残響の使用用途、再生方法に応じて設定する。例えば、２２．２マルチチャンネル方式の場合は、２４本もしくはＬＦＥ（Low Frequency Effect）チャンネルを除いた２２本のマイクロホンを使用することが考えられる。そして、適当な位置に設置したスピーカから、ＴＳＰ信号などのインパルス応答測定用の音響信号を放射し、マイクロホンＭ０〜Ｍｎからなるマイクロホン群で録音する。これにより各マイクロホンＭ０、Ｍ１、Ｍ２、…、Ｍｎのそれぞれにより測定されたインパルス応答ｈ０（ｔ）、ｈ１（ｔ）、ｈ２（ｔ）、…、ｈｎ（ｔ）が得られる。そして、マイクロホン群で囲まれた空間内の任意の点をＰとする。図４（ａ）では、マイクロホンＭ０の設置位置が点Ｐであり、図４（ｂ）では、マイクロホンＭ０〜Ｍｎのいずれも設置されていない位置が点Ｐである。
【００２４】
図５は、インパルス応答生成装置１における方向別インパルス応答生成処理のフロー図である。
まず、インパルス応答生成装置１の測定情報取得部１２は、各マイクロホンＭ０、Ｍ１、Ｍ２、…、Ｍｎにより測定されたインパルス応答ｈ０（ｔ）、ｈ１（ｔ）、ｈ２（ｔ）、…、ｈｎ（ｔ）を取得する（ステップＳ１０）。例えば、測定情報取得部１２は、インパルス応答を各マイクロホンＭ０〜Ｍｎから受信してもよく、他のコンピュータ装置から受信してもよく、記録媒体から読み出してもよい。測定情報取得部１２は、取得したインパルス応答を記憶部１１に書き込む。
【００２５】
仮想音源情報推定部１３は、記憶部１１から読み出したインパルス応答ｈ０（ｔ）、ｈ１（ｔ）、ｈ２（ｔ）、…、ｈｎ（ｔ）を元に、マイクロホン群で囲まれた任意の点Ｐにおける、空間内の仮想音源の情報ｅ１（ｒ１,θ１，φ１，ｔ１）、ｅ２（ｒ２,θ２，φ２，ｔ２）、…、ｅｍ（ｒｍ,θｍ，φｍ，ｔｍ）を予測する（ステップＳ２０）。ここで、ｒは振幅、θは方位角、φは仰角、ｔは到来時刻、ｍは予測された仮想音源の個数であり、ｒ、θ、φ及びｔの添え字は、ｍ個の仮想音源のうち何番目の仮想音源であるかを示す。振幅ｒは仮想音源のレベルを表し、方位角θ及び仰角φは点Ｐにおける仮想音源の到来方向を表す。方位角θ及び仰角φを表すために用いられる３次元座標空間は任意に設定し得る。ここでは、点Ｐの位置を原点０とし、ｘ軸をスピーカと点Ｐとを結ぶ直線、ｙ軸を地面と平行な平面上でｘ軸と垂直な方向、ｚ軸をｘｙ平面と垂直な方向とする３次元座標空間を用いる。また、ここでは、点ＰにマイクロホンＭ０を設置した場合を考える。仮想音源情報推定部１３は、近接四点法やそれを拡張した方法など任意の既存の手法を用いてインテンシティや各チャンネル間の相互相関を算出することにより、仮想音源情報を予測する。算出手順の詳細については後述する。仮想音源情報推定部１３は、予測した仮想音源情報ｅ１〜ｅｍを記憶部１１に書き込む。
【００２６】
次に、比較判定部１４は、各インパルス応答ｈ１（ｔ）〜ｈｎ（ｔ）から除去すべき反射音を判定する比較判定処理を行う（ステップＳ３０）。まず、比較判定部１４は、記憶部１１からインパルス応答ｈ０（ｔ）、ｈ１（ｔ）、ｈ２（ｔ）、…、ｈｎ（ｔ）及び仮想音源情報ｅ１（ｒ１,θ１，φ１，ｔ１）、ｅ２（ｒ２,θ２，φ２，ｔ２）、…、ｅｍ（ｒｍ,θｍ，φｍ，ｔｍ）を読み出す。比較判定部１４は、仮想音源情報推定部１３が推定した、点Ｐにおける仮想音源の情報ｅ１〜ｅｍから、各仮想音源の到来時刻（ｔ１，ｔ２，…，ｔｍ）を抽出する。比較判定部１４は、各仮想音源の到来時刻（ｔ１，ｔ２，…，ｔｍ）と、測定されたインパルス応答ｈ１（ｔ）、ｈ２（ｔ）、…、ｈｎ（ｔ）のそれぞれに含まれる局所的なピークの到来時刻（ｔ１１，ｔ１２，…）、（ｔ２１，ｔ２２，…）、…、（ｔｎ１，ｔｎ２，…）とを比較する。比較判定部１４は、比較結果に基づいて、各インパルス応答から取り除くべき反射音を判定する。取り除くべき反射音とは、あるマイクロホンにより測定されたインパルス応答に含まれる反射音のうち、そのマイクロホンに対応する方向（マイクロホンを向けた方向）とは異なる方向から到来する反射音である。つまり、この取り除くべき反射音とは、その到来方向に向けて配置したマイクロホンだけではなく、その到来方向以外に向けて配置した他のマイクロホンにより測定されたインパルス応答において、重複して測定された反射音である。比較判定部１４は、取り除くべき反射音の判定にあたって、点Ｐから各マイクロホンＭ１、Ｍ２、…、Ｍｎまでの距離に応じた遅延τ１、τ２、…、τｎも考慮する。さらに、判定には、測定誤差等も加味する必要がある。そこで、比較判定部１４は、ある所定の窓幅ｗを設定する。比較判定部１４は、遅延を考慮した仮想音源の到来時刻と、インパルス応答においてその到来時刻を中心とした窓幅ｗの時刻範囲内で局所的なピークが生じた時刻とを比較して判定を行う。
【００２７】
反射音除去部１５は、各インパルス応答ｈ１（ｔ）〜ｈｎ（ｔ）から比較判定部１４が除去すべきと判定した反射音を除去する反射音除去処理を行い、方向別のインパルス応答を生成する（ステップＳ４０）。まず、反射音除去部１５は、記憶部１１から読み出したインパルス応答ｈ１（ｔ）、ｈ２（ｔ）、…、ｈｎ（ｔ）のそれぞれから、比較判定部１４によりそのインパルス応答から除去すべきと判定された反射音を除去し、方向別のインパルス応答ｈ’１（ｔ）、ｈ’２（ｔ）、…、ｈ’ｎ（ｔ）を生成する。
【００２８】
以下に、具体例を交えながら、図５のステップＳ２０〜ステップＳ４０の詳細な処理を説明する。
本来、三次元空間内の仮想音源情報を得るためには、同一平面内にない４点以上にマイクロホンを設置する必要があるが、ここでは原理を分かりやすく説明するため、同一平面内に複数のマイクロホンを置いて測定した場合を考える。
【００２９】
図６は、インパルス応答測定方法の例を示す図である。マイクロホンＭ０を中心とし、それを取り囲むようにマイクロホンＭ１〜Ｍ４を、マイクロホンＭ０から等距離かつ等間隔に設置する。隣り合うマイクロホンＭ１〜Ｍ４とマイクロホンＭ０とがなす角度は９０度であり、マイクロホンＭ０と各マイクロホンＭ１〜Ｍ４との距離はＬである。具体的には、マイクロホンＭ１、Ｍ２、Ｍ３、Ｍ４は、マイクロホンＭ０を原点とした２次元空間のｘｙ平面上に設置される。マイクロホンＭ０、Ｍ１、Ｍ２、Ｍ３、Ｍ４それぞれの設置方向は、ｘ軸に対してそれぞれ０度、０度、９０度、２７０度、１８０度である。
【００３０】
そして、スピーカから測定信号を再生し、各マイクロホンＭ０〜Ｍ４により収音された信号からインパルス応答ｈ０（ｔ）〜ｈ４（ｔ）を算出する。同図に示すように、各インパルス応答ｈｉ（ｔ）には（ｉは１以上４以下の整数）、マイクロホンＭｉに対応する方向（マイクロホンを向けた方向）から到来する反射音のレベル（実線）だけでなく、それ以外の方向から到来する反射音のレベル（破線）も含まれてしまう。指向性マイクロホンで測定すれば、全指向性マイクロホンで測定するよりも、マイクロホンに対応する方向（マイクロホンを向けた方向）以外から到来する反射音を抑制することができる。しかしながら、指向性マイクロホンとは言っても、その指向性には幅があるため、隣接する複数のマイクロホンで同じ反射音を拾ってしまう。そこで、インパルス応答生成装置１は、後述する処理によって、マイクロホン群の中心点Ｍ０（マイクロホンＭ０の位置）における仮想音源情報を推定し、その推定した仮想音源情報を元に、各マイクロホンＭ１〜Ｍ４のインパルス応答ｈ１（ｔ）〜ｈ４（ｔ）から重複する反射音を間引く。
【００３１】
次に、ステップＳ２０において仮想音源情報推定部１３が、測定されたインパルス応答から仮想音源情報を推定する処理を説明する。推定には様々な手法が用いられるが、ここでは例として、インテンシティを用いる手法と、相互相関係数を用いる手法とを示す。
【００３２】
まず、インテンシティを用いた仮想音源情報の推定処理を説明する。
図７は、インテンシティを用いた仮想音源情報の推定処理を説明するための図である。
最初に、仮想音源情報推定部１３は、インパルス応答ｈ０（ｔ）及びインパルス応答ｈ１（ｔ）から０−Ｘ方向（Ｍ０−Ｍ１の方向）の瞬時インテンシティＩｘ（ｔ）を求める。さらに、仮想音源情報推定部１３は、インパルス応答ｈ０（ｔ）及びインパルス応答ｈ２（ｔ）とから０−Ｙ方向（Ｍ０−Ｍ２の方向）の瞬時インテンシティＩｙ（ｔ）を求める。
【００３３】
瞬時インテンシティＩｘ（ｔ）は、以下の式（１）により求められる。
【００３４】
【数１】
0006650245.tif 000002
【００３５】
式（１）において、ｐ１（ｔ）、ｐ２（ｔ）は、マイクロホンで測定された時刻ｔにおける音圧であり、それぞれ、インパルス応答ｈ０（ｔ）、ｈ１（ｔ）を用いる。また、ｄは、マイクロホン間の間隔であり、マイクロホンＭ０とマイクロホンＭ１の距離Ｌを用いる。ρ０は、媒質の密度である。瞬時インテンシティＩｘ（ｔ）の値の正負は、ｘ軸方向の到来方向を表す。
瞬時インテンシティＩｙ（ｔ）も同様に求められる。ｐ１（ｔ）、ｐ２（ｔ）には、インパルス応答ｈ０（ｔ）、ｈ２（ｔ）を用い、ｄには、マイクロホンＭ０とマイクロホンＭ２の距離Ｌを用いる。瞬時インテンシティＩｙ（ｔ）の値の正負は、ｙ軸方向の到来方向を表す。
【００３６】
なお、瞬時インテンシティＩｘ（ｔ）は、ｘ軸上に設置した任意のマイクロホンにより測定されたインパルス応答の組合せ（例えば、インパルス応答ｈ０（ｔ）及びｈ４（ｔ））を用いて算出することができる。同様に、瞬時インテンシティＩｙ（ｔ）は、ｙ軸上に設置した任意のマイクロホンにより測定されたインパルス応答の組合せ（例えば、インパルス応答ｈ０（ｔ）及びｈ３（ｔ））を用いて算出することができる。
【００３７】
続いて、仮想音源情報推定部１３は、ヒルベルト変換により瞬時インテンシティＩｘ（ｔ）、瞬時インテンシティＩｙ（ｔ）それぞれの包絡を取ったエンベロープインテンシティＩｈｘ（ｔ）、Ｉｈｙ（ｔ）を求める。仮想音源情報推定部１３は、エンベロープインテンシティＩｈｘ（ｔ）、Ｉｈｙ（ｔ）のレベルをベクトル合成し、合成インテンシティのレベルＬｅｖ＿Ｉｈ（ｔ）を求める。エンベロープインテンシティＩｈｘ（ｔ）、Ｉｈｙ（ｔ）のレベルをベクトル合成することにより、時刻ｔにおける仮想音源の到来方向も得られるが、同図では、レベルの経時的な変化のみを示している。仮想音源情報推定部１３は、合成インテンシティレベルＬｅｖ＿Ｉｈ（ｔ）の最初のピークレベルＬｅｖ＿Ｐｓ（直接音）から任意のレベルＬｅｖ＿Ｐｅまでの範囲に表れる局所的なピークＰ１、Ｐ２、…、Ｐｍを検出する。仮想音源情報推定部１３は、検出したピークＰ１、Ｐ２、…、Ｐｍにおける合成インテンシティから、室内（測定空間）における仮想音源の情報ｅ１（ｒ１,θ１，φ１，ｔ１）、ｅ２（ｒ２,θ２，φ２，ｔ２）、…、ｅｍ（ｒｍ,θｍ，φｍ，ｔｍ）を得る。なお、到来時刻ｔ１＜ｔ２＜ｔ３＜…＜ｔｍとする。仮想音源の情報ｅ１〜ｅｍを到来時刻によって並べた情報をｅ（ｔ）とする。
【００３８】
上記は、マイクロホンＭ０〜Ｍ４が同一平面（ｘｙ平面）上に設置されている例であるため、仮想音源のＺ座標は０であり、仰角φ１、φ２、…、φｍも０である。ｚ軸上にもマイクロホンを設置した場合、つまり三次元的にマイクロホンを設置した場合は、その設置したマイクロホンにより測定されたインパルス応答の組合せをさらに用いて上記と同様に合成インテンシティを算出することにより、各仮想音源の仰角を算出することができる。つまり、仮想音源情報推定部１３は、ｚ軸上に設置した２つのマイクロホンにより測定されたインパルス応答の組合せを用いて、0−Ｚ方向の瞬時インテンシティＩｚ（ｔ）の包絡を取ったエンベロープインテンシティＩｈｚ（ｔ）を算出する。仮想音源情報推定部１３は、エンベロープインテンシティＩｈｘ（ｔ）、Ｉｈｙ（ｔ）、Ｉｈｚ（ｔ）のレベルをベクトル合成し、Ｌｅｖ＿ＰｓとＬｅｖ＿Ｐｅの範囲に表れる局所的なピークＰ１、Ｐ２、…、Ｐｍを検出する。
【００３９】
次に、相互相関係数を用いた仮想音源情報の推定処理を説明する。
時刻ｔａから時刻ｔａ＋Δにおける２つの各インパルス応答ｈｉ（ｔ）、ｈｊ（ｔ）間の短時間相関係数Ｓｉｊ（τ）は、次式（２）で与えられる。Δは分析窓長であり、τは時間差である。
【００４０】
【数２】
0006650245.tif 000003
【００４１】
異なる位置に設置された各マイクロホンに同一の反射音が到達する時刻にはずれがあるため、各マイクロホンにより測定されたインパルス応答には同一の反射音が異なる時刻に表れる。短時間相関係数Ｓｉｊの値が大きいということは、同一の反射音が時間差τだけずれてｈｉ（ｔ）とｈｊ（ｔ）に出現することを表す。つまり、同一の反射音がτだけずれてマイクロホンＭｉ、マイクロホンＭｊに到達したことを表す。本実施形態では、ｘ軸上、ｙ軸上のそれぞれに設置した任意のマイクロホンにより測定されたインパルス応答の組合せを用いるため、ここでは、ｉ＝０，１，２、かつ、ｊ＝１，２（ｉ≠ｊ）とする。つまり、（ｉ，ｊ）の組合せには、（０，１）、（０，２）、（１，２）又は（２，１）の３つがある。以下では、ｉ＝０，１，２かつｊ＝１，２である場合を例に説明するが、ｘ軸上、ｙ軸上のそれぞれに設置した任意のマイクロホンにより測定されたインパルス応答の組合せ（例えば、ｉ＝０，３，４かつｊ＝３，４など）を用いることができる。
【００４２】
図８は、相互相関係数を用いた仮想音源情報推定処理のフロー図である。
まず、仮想音源情報推定部１３は、時刻ｔａに初期値を設定する（ステップＳ１１０）。初期値は、例えば、インパルス応答ｈ０（ｔ）において直接音が検出された時刻である。仮想音源情報推定部１３は、各（ｉ，ｊ）の組合せについて、相関係数Ｓｉｊ（τ）に極大値を与える時間差τ＝τｉｊを見つける（ステップＳ１２０）。仮想音源情報推定部１３は、各（ｉ，ｊ）の組合せごとに、相関係数Ｓｉｊ（τ）がある一定の値（例えば、０．８）以上となる時間長（窓幅）Δｉｊを求め、求めた中で最小の値の時間長Δｉｊ＿ｍｉｎを選択する（ステップＳ１３０）。時間長Δｉｊ＿ｍｉｎは、ある反射音が到達してから、次の反射音が到達するまでの最も短い時間に相当する。時間長Δｉｊの最小値を選択することにより、窓幅の中に次の反射音が含まれないようにする。
【００４３】
次に、仮想音源情報推定部１３は、この時間長Δｉｊ＿ｍｉｎ付近で補間処理を行う。ここでは、補間処理として、サンプリング周波数ｆｓを数十倍にあげる。本実施形態では、サンプリング周波数ｆｓを１６倍に上げる。仮想音源情報推定部１３は、補間処理により得られたインパルス応答ｈｉ（ｔ）、ｈｊ（ｔ）を用い、ステップＳ１２０と同様の処理を行い、相関係数Ｓｉｊ（τ）に極大値を与え、かつ、より精度の高い時間差τ＝ｈτｉｊを求める（ステップＳ１４０）。仮想音源情報推定部１３は、求めた時間差ｈτｉｊを用いて仮想音源の座標及び角度を算出する（ステップＳ１５０）。
【００４４】
仮想音源とマイクロホンＭ０の距離は、ｃを音速、仮想音源が放射された時刻をｔ０とすると、（ｔａ−ｔ０）×ｃにより得られる。また、仮想音源とマイクロホンＭ１の距離は、（ｔａ＋ｈτ０１−ｔ０）×ｃにより、仮想音源とマイクロホンＭ２の距離は、（ｔａ＋ｈτ０２−ｔ０）×ｃにより得られる。マイクロホンＭ０とｘ軸上のマイクロホンＭ１との距離、及び、マイクロホンＭ０とｙ軸上のマイクロホンＭ２との距離は既知（距離Ｌ）である。よって、これらのマイクロホン間の距離、及び、仮想音源と各マイクロホンＭ０、Ｍ１、Ｍ２との距離とに基づいて、時刻ｔａにおいて中心点Ｍ０に到達した仮想音源の座標（Ｘ，Ｙ，Ｚ）が得られ、この座標から方位角及び仰角が得られる。
【００４５】
仮想音源情報推定部１３は、時刻ｔａが分析終了の時刻に達していなければ（ステップＳ１６０：ＮＯ）、分析する窓のスタート地点である時刻ｔａをΔｉｊ＿ｍｉｎだけずらして（ステップＳ１７０）、ステップＳ１２０からの処理を行う。仮想音源情報推定部１３は、時刻ｔａが分析終了の時刻に達するまで（ステップＳ１６０：ＹＥＳ）、次々に仮想音源を検出していく。このような処理により、仮想音源情報推定部１３は、マイクロホン群の中心点Ｍ０における仮想音源情報ｅ１（ｒ１,θ１，φ１，ｔ１）、ｅ２（ｒ２,θ２，φ２，ｔ２）、…、ｅｍ（ｒｍ,θｍ，φｍ，ｔｍ）を推定することができる。
【００４６】
上記では、マイクロホンＭ０〜Ｍ４が同一平面（ｘｙ平面）上に設置されているため、仮想音源のＺ座標は０となり、従って、仰角φ１、φ２、…、φｍも０となる。ｘｙ平面上にないｚ軸上のマイクロホンなど、同一平面上にないマイクロホンにより測定されたインパルス応答をさらに用いて上記と同様の処理を行うことで、各仮想音源のＺ座標が求められ、仰角も算出できる。
【００４７】
次に、図５のステップＳ３０における比較判定処理及び図５のステップＳ４０における反射音除去処理を説明する。
図９は、比較判定処理及び反射音除去処理の例を示す図である。ここでは、仮想音源情報推定の結果、４つの仮想音源情報ｅ１（ｒ１,θ１，φ１，ｔ１）、ｅ２（ｒ２,θ２，φ２，ｔ２）、ｅ３（ｒ３,θ３，φ３，ｔ３）、ｅ４（ｒ４,θ４，φ４，ｔ４）が推定された場合を例に説明する。比較判定部１４は、中心点Ｍ０における仮想音源の情報ｅ１〜ｅ４から各仮想音源の到来時刻（ｔ１，ｔ２，ｔ３，ｔ４）を抽出する。比較判定部１４は、各仮想音源の到来時刻（ｔ１，ｔ２，ｔ３，ｔ４）と、測定されたインパルス応答ｈ１（ｔ）、ｈ２（ｔ）、ｈ３（ｔ）、ｈ４（ｔ）のそれぞれに含まれる局所的なピークの到来時刻（ｔ１１，ｔ１２，…），（ｔ２１，ｔ２２，…），（ｔ３１，ｔ３２，…），（ｔ４１，ｔ４２，…）とを比較し、取り除くべき反射音を判定する。取り除くべき反射音とは、あるインパルス応答において、そのインパルス応答が測定されたマイクロホンに対応する方向（マイクロホンを向けた方向）以外から到来する反射音である。
【００４８】
判定の流れを下記に示す。なお、ここで、ｍは予測した仮想音源の数（図９ではｍ＝４）であり、ｎはマイクロホンの数（図９ではｎ＝４）である。
【００４９】
（手順１）比較判定部１４は、仮想音源情報ｅｋ（ｋは１以上ｍ以下の整数）ごとに、方位角θｋに基づいて仮想音源の到来方向に対応するマイクロホンＭｉ（ｉは１以上ｎ以下の整数）を特定し、特定したマイクロホンＭｉにより測定されたインパルス応答ｈｉ（ｔ）を選択する。具体的には、比較判定部１４は、設置方向が方位角θｋであるマイクロホンにより測定されたインパルス応答を選択する。方位角θｋが、いずれのマイクロホンＭ１〜Ｍｎの設置方向とも一致しない場合、比較判定部１４は、下記の（１）又は（２）のいずれかによりインパルス応答を選択する。
【００５０】
（１）比較判定部１４は、θｋ方向の両隣に存在するマイクロホンのうち、設置方向がθｋに近い方のマイクロホンのインパルス応答を選択する。この場合、図６に示すマイクロホンＭ１、Ｍ２、Ｍ４、Ｍ３で測定すべき反射音の方向は、例えば、方位角がｘ軸に対してそれぞれ０〜４５及び３１５〜３６０度、４５〜１３５度、１３５〜２２５度、２２５〜３１５度である。
【００５１】
（２）比較判定部１４は、θｋ方向の両隣に存在するマイクロホンのそれぞれにより測定された２つのインパルス応答を選択する。この場合、マイクロホンＭ１、Ｍ２、Ｍ４、Ｍ３で測定すべき反射音の方向は、０〜９０及び２７０〜３６０度、０〜１８０度、９０〜２７０度、１８０〜３６０度である。
【００５２】
（手順２）比較判定部１４は、手順１により仮想音源情報ｅｋについて選択したインパルス応答ｈｉ（ｔ）において、ｔｋ−τｉ−ｗ／２≦ｔ≦ｔｋ−τｉ＋ｗ／２の範囲内で、局所的なピークが生ずる時刻ｔ＝ｔｉｓを選択する。ここで、ｔｋは、中心点Ｍ０における仮想音源の到来時刻、τｉは手順１において選択されたマイクロホンＭｉから中心点Ｍ０までの距離Ｌに応じた遅延、ｗは窓幅である。なお、窓幅ｗは予め設定されるが、隣接するマイクロホンとの距離が長い程大きな値とすることが望ましい。また、ｓは各インパルス応答における局所的なピークのインデックスであり、ｔｉｓはインパルス応答ｈｉ（ｔ）の波形構造においてｓ番目の局所的なピークが生ずる時刻である。
【００５３】
（手順３）比較判定部１４は、推定した仮想音源情報ごとに（手順１）及び（手順２）を繰り返し、各インパルス応答において選ばれなかった局所的なピークを、そのインパルス応答から除去すべき反射音であると判定する。反射音除去部１５は、各インパルス応答から、そのインパルス応答について比較判定部１４が除去すべきと判定した局所的なピークを除去する。
【００５４】
図９に示した例では、まず、比較判定部１４は、仮想音源情報ｅ１（ｒ１,θ１，φ１，ｔ１）の到来方向を表すθ１がマイクロホンＭ１の設置方向であるため、マイクロホンＭ１により測定されたインパルス応答ｈ１（ｔ）を選択する。比較判定部１４は、インパルス応答ｈ１（ｔ）において、ｔ１−τ１−ｗ／２≦ｔ≦ｔ１−τ１＋ｗ／２の範囲内で局所的なピークが生ずる時刻ｔ＝ｔ１１を選択する。
【００５５】
次に、仮想音源情報ｅ２（ｒ２，θ２，φ２，ｔ２）の到来方向を表すθ２は、マイクロホンＭ１の設置方向とマイクロホンＭ２の設置方向の間の方位角であり、仮想音源はマイクロホンＭ１とマイクロホンＭ２の間から到来する。比較判定部１４は、マイクロホンＭ１の設置方向よりもマイクロホンＭ２の設置方向のほうがθ２に近いため、マイクロホンＭ２により測定されたインパルス応答ｈ２（ｔ）を選択する（上記の手順１の（１）を適用）。比較判定部１４は、インパルス応答ｈ２（ｔ）において、ｔ２−τ２−ｗ／２≦ｔ≦ｔ２−τ２＋ｗ／２の範囲内で局所的なピークが生ずる時刻ｔ＝ｔ２２を選択する。
【００５６】
比較判定部１４は、仮想音源情報ｅ３（ｒ３，θ３，φ３，ｔ３）及びｅ４（ｒ４，θ４，φ４，ｔ４）についても同様に処理し、インパルス応答ｈ３（ｔ）において時刻ｔ＝ｔ３２、インパルス応答ｈ４（ｔ）において時刻ｔ＝ｔ４４を選択する。
【００５７】
反射音除去部１５は、インパルス応答ｈ１（ｔ）の時刻ｔ１１におけるピークＰ（ｔ１１）、インパルス応答ｈ２（ｔ）の時刻ｔ２２におけるピークＰ（ｔ２２）、インパルス応答ｈ３（ｔ）の時刻ｔ３２におけるピークＰ（ｔ３２）、及び、インパルス応答ｈ４（ｔ）の時刻ｔ４４におけるピークＰ（ｔ４４）以外の局所的なピークを除去する。これより、反射音除去部１５は、新たな方向別インパルス応答ｈ１’（ｔ），ｈ２’（ｔ），ｈ３’（ｔ），ｈ４’（ｔ）を生成する。
【００５８】
以上説明した実施形態によって、従来の技術よりも、測定音場に近似した方向別のインパルス応答を生成することが可能となる。
なお、平面上にない複数のマイクロホンを使用する場合、比較判定部１４は、方位角θｋだけでなく、仰角φｋをさらに加えて表される仮想音源の到来方向に基づいて、インパルス応答を選択してもよい。つまり、比較判定部１４は、方位角θｋ及び仰角φｋにより表される到来方向に主軸を向けたマイクロホン、あるいは、その到来方向の両隣に存在するマイクロホン又はそれらマイクロホンのうち主軸を向けた方向が到来方向に近いマイクロホンを特定し、特定したマイクロホンにより測定されたインパルス応答を選択する。
また、上記で示した例では中心点Ｍ０と、それを取り囲むようにマイクロホン群を等間隔に設置しているが、必ずしもその必要はなく、任意のマイクロホン配置から、マイクロホン群で囲まれた空間内の任意の点における仮想音源情報を推定することができる。
【００５９】
上述した実施形態によれば、インパルス応答生成装置１は、空間内の異なる位置に設置された複数のマイクロホンが同一の音響信号を収音して得られた信号それぞれに基づき算出された各インパルス応答を用いて、複数のマイクロホンに囲まれた任意の位置における、空間内の仮想音源の情報を推定する。インパルス応答生成装置１は、仮想音源の情報とインパルス応答との比較に基づいて、任意のインパルス応答において、そのインパルス応答の算出の元となった信号が得られたマイクロホンに対応する方向から到来する反射音を判定する。インパルス応答生成装置１は、インパルス応答から、そのインパルス応答の算出の元となった信号が得られたマイクロホンに対応する方向以外から到来する余計な反射音を間引いて方向別のインパルス応答を生成する。これにより、インパルス応答生成装置１は、従来の手法よりも、測定音場に近い方向別のインパルス応答を生成することができる。
【００６０】
なお、上述のインパルス応答生成装置１は、内部にコンピュータシステムを有している。そして、インパルス応答生成装置１の動作の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータシステムが読み出して実行することによって、上記処理が行われる。ここでいうコンピュータシステムとは、ＣＰＵ及び各種メモリやＯＳ、周辺機器等のハードウェアを含むものである。
【００６１】
また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【符号の説明】
【００６２】
１  インパルス応答生成装置
１１  記憶部
１２  測定情報取得部
１３  仮想音源情報推定部
１４  比較判定部
１５  反射音除去部
【図１】
0006650245.tif 000004
【図２】
0006650245.tif 000005
【図３】
0006650245.tif 000006
【図４】
0006650245.tif 000007
【図５】
0006650245.tif 000008
【図６】
0006650245.tif 000009
【図７】
0006650245.tif 000010
【図８】
0006650245.tif 000011
【図９】
0006650245.tif 000012