(57)【特許請求の範囲】
【請求項１】
球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成部と、
前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成部と
を備える音場再現装置。
【請求項２】
前記第１の駆動信号生成部は、前記収音信号から得られた空間周波数スペクトルに対して空間フィルタを用いたフィルタ処理を施すことで、前記収音信号を前記仮想スピーカアレイの駆動信号に変換する
請求項１に記載の音場再現装置。
【請求項３】
前記収音信号から得られた時間周波数スペクトルを前記空間周波数スペクトルに変換する空間周波数分析部をさらに備える
請求項２に記載の音場再現装置。
【請求項４】
前記第２の駆動信号生成部は、前記実スピーカアレイから前記仮想スピーカアレイまでの伝達関数に基づく逆フィルタを用いて、前記仮想スピーカアレイの駆動信号に対してフィルタ処理を施すことで、前記仮想スピーカアレイの駆動信号を前記実スピーカアレイの駆動信号に変換する
請求項１乃至請求項３の何れか一項に記載の音場再現装置。
【請求項５】
前記仮想スピーカアレイは球状または環状のスピーカアレイである
請求項１乃至請求項４の何れか一項に記載の音場再現装置。
【請求項６】
球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成ステップと、
前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成ステップと
を含む音場再現方法。
【請求項７】
球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成ステップと、
前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成ステップと
を含む処理をコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は音場再現装置および方法、並びにプログラムに関し、特に、より正確に音場を再現することができるようにした音場再現装置および方法、並びにプログラムに関する。
【背景技術】
【０００２】
従来、実空間にて球状または環状のマイクアレイで収音した信号を用いて、再現空間で実空間と同様の音場を再現する技術が提案されている。【０００３】
例えば、そのような技術として、コンパクトな球状マイクアレイによる収音およびスピーカアレイによる再生を可能とするものが提案されている（例えば、非特許文献１参照）。【０００４】
また、例えば任意のアレイ形状のスピーカアレイで再生可能であり、かつ予めスピーカからマイクロホンまでの伝達関数を収録し、逆フィルタを生成しておくことで個々のスピーカの特性の差を吸収することを可能とするものも提案されている（例えば、非特許文献２参照）。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Zhiyun Li et al,“Capture and Recreation of Higher Order 3D Sound Fields via Reciprocity,”Proceedings of ICAD 04-Tenth Meeting of the International Conference on Auditory Display, Sydney, 2004
【非特許文献２】伊勢史郎,“境界音場制御,”日本音響学会誌, 67巻, 11号, 2011
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところが、非特許文献１に記載の技術では、コンパクトな球状マイクアレイによる収音およびスピーカアレイによる再生が可能であるが、厳密な音場再現のためにはスピーカアレイの形状が球状または環状であり、かつスピーカが等密度の配置でなければならないという制約が求められる。【０００７】
例えば図１の左側に示すように、スピーカアレイSPA11を構成する各スピーカが環状に配置されており、図中、点線で表される基準点に対して、各スピーカ同士が等密度（簡略化のため図では等角度）の配置となる場合には、厳密な音場再現が可能である。この例では、互いに隣接する任意の２つのスピーカについて、一方のスピーカおよび基準点を結ぶ直線と、他方のスピーカおよび基準点を結ぶ直線とがなす角度が一定の角度となっている。【０００８】
これに対して、図中、右側に示すように正方形で等間隔に並んだスピーカからなるスピーカアレイSPA12の場合、スピーカ同士が、図中、点線で表される基準点から等密度とならないため、厳密に音場再現することができない。この例では、互いに隣接する２つのスピーカの一方および基準点を結ぶ直線と、他方のスピーカおよび基準点を結ぶ直線とがなす角度が、隣接する２つのスピーカの組ごとに異なる角度となっている。【０００９】
また、モノポール音源を発するような理想スピーカアレイを想定した駆動信号が生成されるため、実際のスピーカの特性の影響で実空間の音場を正確に再現することができなかった。【００１０】
さらに、非特許文献２に記載の技術では、任意のアレイ形状で再生可能であり、かつ予めスピーカからマイクロホンまでの伝達関数を収録し逆フィルタを生成しておけば、個々のスピーカの特性の差を吸収することが可能であった。一方で、予め収録された各スピーカから各マイクロホンへの伝達関数群がそれぞれ似通った性質を保つ場合、伝達関数から駆動信号を生成するための、安定した逆フィルタを求めることが困難であった。【００１１】
特に図２の右側に示す、球状マイクアレイMKA11を用いた例のように、球状マイクアレイMKA11を構成するマイクロホン同士が接近している場合、正方形で等間隔に並んだスピーカからなるスピーカアレイSPA21の特定のスピーカからの全てのマイクロホンへの距離がほぼ等距離となる。そのため、逆フィルタの安定解を求めることが困難であった。【００１２】
なお、図２中、左側には、スピーカアレイSPA21のスピーカからの、球状マイクアレイMKA21を構成する各マイクロホンへの距離が等距離とならず、伝達関数のばらつきが大きくなる例について示されている。この例では、スピーカアレイSPA21のスピーカからの各マイクロホンへの距離が異なるので、逆フィルタの安定解を求めることができる。しかしながら、逆フィルタの安定解を求めることができる程度に球状マイクアレイMKA21の半径を大きくすることは現実的ではない。【００１３】
本技術は、このような状況に鑑みてなされたものであり、より正確に音場を再現することができるようにするものである。
【課題を解決するための手段】
【００１４】
本技術の一側面の音場再現装置は、球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成部と、前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成部とを備える。【００１５】
前記第１の駆動信号生成部には、前記収音信号から得られた空間周波数スペクトルに対して空間フィルタを用いたフィルタ処理を施すことで、前記収音信号を前記仮想スピーカアレイの駆動信号に変換させることができる。【００１６】
音場再現装置には、前記収音信号から得られた時間周波数スペクトルを前記空間周波数スペクトルに変換する空間周波数分析部をさらに設けることができる。【００１７】
前記第２の駆動信号生成部には、前記実スピーカアレイから前記仮想スピーカアレイまでの伝達関数に基づく逆フィルタを用いて、前記仮想スピーカアレイの駆動信号に対してフィルタ処理を施すことで、前記仮想スピーカアレイの駆動信号を前記実スピーカアレイの駆動信号に変換させることができる。【００１８】
前記仮想スピーカアレイを球状または環状のスピーカアレイとすることができる。【００１９】
本技術の一側面の音場再現方法またはプログラムは、球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成ステップと、前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成ステップとを含む。【００２０】
本技術の一側面においては、球状または環状のマイクアレイが収音することで得られた収音信号が、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換され、前記仮想スピーカアレイの駆動信号が、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換される。
【発明の効果】
【００２１】
本技術の一側面によれば、より正確に音場を再現することができる。【００２２】
なお、ここに記載された効果は必ずしも限定されるものではなく、本開示中に記載された何れかの効果であってもよい。
【図面の簡単な説明】
【００２３】
【図１】従来の音場再現について説明する図である。
【図２】従来の音場再現について説明する図である。
【図３】本技術の音場再現について説明する図である。
【図４】本技術の音場再現の他の例について説明する図である。
【図５】音場再現器の構成例を示す図である。
【図６】実スピーカアレイ駆動信号生成処理を説明するフローチャートである。
【図７】音場再現システムの構成例を示す図である。
【図８】音場再現処理を説明するフローチャートである。
【図９】コンピュータの構成例を示す図である。
【発明を実施するための形態】
【００２４】
以下、図面を参照して、本技術を適用した実施の形態について説明する。
【００２５】
〈第１の実施の形態〉
〈本技術について〉
本技術では、実空間にて球状または環状のマイクアレイで収音した信号が用いられて、再現空間で実空間と同様の音場が再現されるように、実スピーカアレイの駆動信号が生成される。その際、マイクアレイは十分に小さくコンパクトなものが想定される。
【００２６】
また、実スピーカアレイの内側または外側に、球状または環状の仮想スピーカアレイが配置される。そして、第一の信号処理により、マイクアレイ収音信号から仮想スピーカアレイ駆動信号が生成される。また、第二の信号処理により、仮想スピーカアレイ駆動信号から実スピーカアレイ駆動信号が生成される。
【００２７】
例えば、図３に示す例では、実空間の球面波が、球状マイクアレイ１１で収音され、再現空間に正方形に配置された実スピーカアレイ１２に対して、その内側に配置された仮想スピーカアレイ１３の駆動信号から求めた駆動信号を供給することにより、実空間の音場が再現されている。
【００２８】
図３では、球状マイクアレイ１１は複数のマイクロホン（マイクセンサ）からなり、各マイクロホンは、所定の基準点を中心とする球の表面に配置されている。以下では、球状マイクアレイ１１を構成するスピーカが配置されている球の中心を球状マイクアレイ１１の中心とも称し、その球の半径を球状マイクアレイ１１の半径、またはセンサ半径とも称することとする。
【００２９】
また、実スピーカアレイ１２は複数のスピーカからなり、それらのスピーカが正方形状に並べられて配置されている。この例では、所定の基準点にいるユーザを囲むように水平面上に実スピーカアレイ１２を構成するスピーカが並べられている。
【００３０】
なお、実スピーカアレイ１２を構成するスピーカの配置は図３に示す例に限らず、所定の基準点を囲むように各スピーカが配置されていればよい。したがって、例えば実スピーカアレイを構成する各スピーカが部屋の天井や壁に設けられているようにしてもよい。
【００３１】
さらに、この例では、実スピーカアレイ１２の内側に、複数の仮想のスピーカが並べられて得られる仮想スピーカアレイ１３が配置されている。すなわち、仮想スピーカアレイ１３を構成するスピーカにより囲まれる空間の外側に実スピーカアレイ１２が配置されている。この例では、仮想スピーカアレイ１３を構成する各スピーカは、所定の基準点を中心として円形状（環状）に並べられており、それらのスピーカは図１に示したスピーカアレイSPA11と同様に、基準点に対して等密度で並ぶように配置されている。
【００３２】
以下では、仮想スピーカアレイ１３を構成するスピーカが配置されている円の中心を、仮想スピーカアレイ１３の中心とも称し、その円の半径を仮想スピーカアレイ１３の半径とも称することとする。
【００３３】
ここで、再現空間では、仮想スピーカアレイ１３の中心位置、つまり基準点は、再現空間において想定される球状マイクアレイ１１の中心位置（基準点）と同じ位置とされる必要がある。なお、仮想スピーカアレイ１３の中心位置と実スピーカアレイ１２の中心位置とは必ずしも同じ位置である必要はない。
【００３４】
本技術では、まず球状マイクアレイ１１で得られた収音信号から、仮想スピーカアレイ１３で実空間の音場を再現するための仮想スピーカアレイ駆動信号が生成される。仮想スピーカアレイ１３は円形状（環状）であり、その中心からみて各スピーカが等密度（等間隔）で配置されているので、正確に実空間の音場を再現することのできる仮想スピーカアレイ駆動信号が生成される。
【００３５】
さらに、このようにして得られた仮想スピーカアレイ駆動信号から、実スピーカアレイ１２で実空間の音場を再現するための実スピーカアレイ駆動信号が生成される。
【００３６】
このとき、実スピーカアレイ１２の各スピーカから、仮想スピーカアレイ１３の各スピーカまでの伝達関数から得られる逆フィルタが用いられて実スピーカアレイ駆動信号が生成される。したがって、実スピーカアレイ１２の形状を任意の形状とすることができる。
【００３７】
このように、本技術では、収音信号から一旦、環状または球状の仮想スピーカアレイ１３の仮想スピーカアレイ駆動信号を生成し、さらにその仮想スピーカアレイ駆動信号を実スピーカアレイ駆動信号に変換することで、実スピーカアレイ１２の形状によらず、正確に音場を再現することができる。
【００３８】
なお、以下では、図３に示したように実スピーカアレイ１２の内側に仮想スピーカアレイ１３が配置されている場合を例として説明するが、例えば図４に示すように実スピーカアレイ２１が、仮想スピーカアレイ２２を構成するスピーカにより囲まれる空間の内側に配置されているようにしてもよい。なお、図４において図３における場合と対応する部分には同一の符号を付してあり、その説明は適宜省略する。
【００３９】
図４の例では、実スピーカアレイ２１を構成する各スピーカは、所定の基準点を中心とする円上に配置されている。また、仮想スピーカアレイ２２を構成する各スピーカも所定の基準点を中心とする円上に等間隔で配置されている。
【００４０】
したがって、この例では上述した第一の信号処理によって、仮想スピーカアレイ２２で音場を再現するための仮想スピーカアレイ駆動信号が収音信号から生成される。また、第二の信号処理によって、仮想スピーカアレイ２２の半径よりも、より小さい半径の円上に配置されたスピーカからなる実スピーカアレイ２１で音場を再現するための実スピーカアレイ駆動信号が仮想スピーカアレイ駆動信号から生成される。
【００４１】
例えば、図３に示した実スピーカアレイ１２として家などの部屋の壁に設けられたスピーカアレイを想定しており、図４に示した実スピーカアレイ２１としてユーザの頭部を囲むポータブルなスピーカアレイを想定している。これらの図３および図４に示した例では、上述した第一の信号処理によって得られた仮想スピーカアレイ駆動信号を共通で使用することができる。
【００４２】
本技術によれば、例えば実空間において、直径が人の頭部程度の球状または環状のマイクアレイで音場を保存する収音部を備え、再現空間において、実空間と同様の音場となるよう、上記マイクアレイより直径が大きな球状または環状の仮想スピーカアレイへの駆動信号を生成する第一の駆動信号生成部を備え、上記駆動信号を、上記仮想スピーカアレイが囲む空間より内側または外側に配置された任意の形状の実スピーカアレイへ信号変換する第二の駆動信号生成部を備えるような音場再現装置を実現することができる。
【００４３】
そして本技術によれば、以下の効果（１）乃至効果（３）を得ることができる。
【００４４】
効果（１）
コンパクトな球状または環状のマイクアレイで収音した信号を、任意のアレイ形状から音場再現が可能である。
効果（２）
逆フィルタの計算の際、実録した伝達関数を使用することにより、スピーカ特性のばらつきや再現空間の反射特性を吸収した駆動信号を生成することが可能である。
効果（３）
球状または環状の仮想スピーカアレイの半径を広げることにより、伝達関数の逆フィルタを安定的に解くことが可能である。
【００４５】
〈音場再現器の構成例〉
次に、本技術を音場再現器に適用した場合を例として、本技術を適用した具体的な実施の形態について説明する。
【００４６】
図５は、本技術を適用した音場再現器の一実施の形態の構成例を示す図である。
【００４７】
音場再現器４１は、駆動信号生成器５１および逆フィルタ生成器５２を有している。
【００４８】
駆動信号生成器５１は、球状マイクアレイ１１を構成する各マイクロホン、すなわちマイクセンサが収音することで得られた収音信号に対して、逆フィルタ生成器５２で得られる逆フィルタを用いたフィルタ処理を施し、その結果得られた実スピーカアレイ駆動信号を実スピーカアレイ１２に供給し、音声を出力させる。すなわち、逆フィルタ生成器５２で生成された逆フィルタが使用されて、実際に音場再現を行うための実スピーカアレイ駆動信号が生成される。
【００４９】
逆フィルタ生成器５２は、入力された伝達関数に基づいて逆フィルタを生成し、駆動信号生成器５１に供給する。
【００５０】
ここで、逆フィルタ生成器５２に入力される伝達関数は、例えば図３に示した実スピーカアレイ１２を構成する各スピーカから仮想スピーカアレイ１３を構成する各スピーカ位置までのインパルスレスポンスとされる。
【００５１】
駆動信号生成器５１は、時間周波数分析部６１、空間周波数分析部６２、空間フィルタ適用部６３、空間周波数合成部６４、逆フィルタ適用部６５、および時間周波数合成部６６を有している。
【００５２】
また、逆フィルタ生成器５２は、時間周波数分析部７１および逆フィルタ生成部７２を有している。
【００５３】
以下、駆動信号生成器５１および逆フィルタ生成器５２を構成する各部について詳細に説明する。
【００５４】
（時間周波数分析部）
時間周波数分析部６１は、実空間の基準点に中心が合うように設置された球状マイクアレイ１１の各マイクセンサの位置Omic(p)＝［apcosθpcosφp，apsinθpcosφp，apsinφp］における収音信号s(p,t)の時間周波数情報を分析する。
【００５５】
但し、位置Omic(p)において、ａpはセンサ半径、つまり球状マイクアレイ１１の中心位置から、その球状マイクアレイ１１を構成する各マイクセンサ（マイクロホン）までの距離を示しており、θpはセンサ方位角を示しており、φpはセンサ仰角を示している。センサ方位角θpおよびセンサ仰角φpは、球状マイクアレイ１１の中心から見た各マイクセンサの方位角および仰角である。したがって、位置p（位置Omic(p)）は極座標で表現された球状マイクアレイ１１の各マイクセンサの位置を示している。
【００５６】
なお、以下では、センサ半径ａpを単にセンサ半径ａとも記すこととする。また、この実施の形態では、球状マイクアレイ１１を用いているが、水平面の音場のみが収録できる環状マイクアレイを用いても構わない。
【００５７】
初めに、時間周波数分析部６１は、収音信号s(p,t)から固定サイズの時間フレーム分割を行った入力フレーム信号sfr(p,n,l)を得る。そして、時間周波数分析部６１は、次式（１）に示す窓関数wana(n)を入力フレーム信号sfr(p,n,l)に乗算し、窓関数適用信号sw(p,n,l)を得る。すなわち、以下の式（２）の計算が行われて、窓関数適用信号sw(p,n,l)が算出される。
【００５８】
【数１】
0006458738.tif 000002
【００５９】
【数２】
0006458738.tif 000003
【００６０】
ここで、式（１）および式（２）において、ｎは時間インデクスを示しており、時間インデクスｎ＝0,…,Nfr-1である。また、ｌは時間フレームインデクスを示しており、時間フレームインデクスｌ＝0,…,L-1である。なお、Nfrはフレームサイズ（時間フレームのサンプル数）であり、Ｌは総フレーム数である。
【００６１】
また、フレームサイズNfrは、サンプリング周波数fsにおける一フレームの時間fsec相当のサンプル数Nfr（＝R(fs×fsec)、但しR()は任意の丸め関数）である。この実施の形態では、例えば一フレームの時間fsec＝0.02[s]であり、丸め関数R()は四捨五入であるが、それ以外でも構わない。さらに、フレームのシフト量はフレームサイズNfrの50％としているが、それ以外でも構わない。
【００６２】
さらに、ここでは窓関数としてハニング窓の平方根を用いているが、ハミング窓やブラックマンハリス窓などのその他の窓を用いるようにしてもよい。
【００６３】
このようにして窓関数適用信号sw(p,n,l)が得られると、時間周波数分析部６１は、以下の式（３）および式（４）を計算することで、窓関数適用信号sw(p,n,l)に対して時間周波数変換を行い、時間周波数スペクトルS(p,ω,l)を得る。
【００６４】
【数３】
0006458738.tif 000004
【００６５】
【数４】
0006458738.tif 000005
【００６６】
すなわち、式（３）の計算によりゼロ詰め信号sw’(p,q,l)が求められ、得られたゼロ詰め信号sw’(p,q,l)に基づいて式（４）が計算され、時間周波数スペクトルS(p,ω,l)が算出される。
【００６７】
なお、式（３）および式（４）において、Ｑは時間周波数変換に用いるポイント数を示しており、式（４）においてｉは純虚数を示している。また、ωは時間周波数インデクスを示している。ここで、Ω＝Ｑ／２＋１とすると、ω＝0,…,Ω-1である。
【００６８】
したがって、球状マイクアレイ１１の各マイクロホンから出力された収音信号ごとに、L×Ω個の時間周波数スペクトルS(p,ω,l)が得られる。
【００６９】
また、この実施の形態では、DFT（Discrete Fourier Transform）（離散フーリエ変換）による時間周波数変換を行っているが、DCT（Discrete Cosine Transform）（離散コサイン変換）やMDCT（Modified Discrete Cosine Transform）（修正離散コサイン変換）などの他の時間周波数変換を用いてもよい。
【００７０】
さらに、DFTのポイント数Ｑは、Nfr以上である、Nfrに最も近い２のべき乗の値としているが、それ以外のポイント数Ｑでも構わない。
【００７１】
時間周波数分析部６１は、以上において説明した処理で得られた時間周波数スペクトルS(p,ω,l)を、空間周波数分析部６２に供給する。
【００７２】
また、逆フィルタ生成器５２の時間周波数分析部７１も、実スピーカアレイ１２のスピーカから仮想スピーカアレイ１３のスピーカまでの伝達関数に対して時間周波数分析部６１と同様の処理を行って、得られた時間周波数スペクトルを逆フィルタ生成部７２に供給する。
【００７３】
（空間周波数分析部）
続いて空間周波数分析部６２は、時間周波数分析部６１から供給された時間周波数スペクトルS(p,ω,l)の空間周波数情報を分析する。
【００７４】
例えば、空間周波数分析部６２は、次式（５）を計算することで、球面調和関数Yn-m(θ,φ)による空間周波数変換を行い、空間周波数スペクトルSnm(a,ω,l)を得る。但し、Ｎは球面調和関数の次数であり、ｎ＝0,…,Nである。
【００７５】
【数５】
0006458738.tif 000006
【００７６】
なお、式（５）において、Pは球状マイクアレイ１１のセンサ数、つまりマイクセンサの数を示しており、ｎは次数を示している。また、θpはセンサ方位角を示しており、φpはセンサ仰角を示しており、aは球状マイクアレイ１１のセンサ半径を示している。ωは時間周波数インデクスを示しており、ｌは時間フレームインデクスを示している。
【００７７】
さらに、球面調和関数Ynm(θ,φ)は、次式（６）に示すように、ルジャンドル随伴多項式Pnm(z)によって与えられる。球面調和関数の最大次数Ｎはセンサ数Pによって制限され、Ｎ＝（Ｐ＋１）２である。
【００７８】
【数６】
0006458738.tif 000007
【００７９】
このようにして得られる空間周波数スペクトルSnm(a,ω,l)は、時間フレームｌに含まれている時間周波数ωの信号が空間上においてどのような波形となっているかを示しており、時間フレームｌごとにΩ×P個の空間周波数スペクトルが得られる。
【００８０】
空間周波数分析部６２は、以上で説明した処理により得られた空間周波数スペクトルSnm(a,ω,l)を空間フィルタ適用部６３に供給する。
【００８１】
（空間フィルタ適用部）
空間フィルタ適用部６３は、空間周波数分析部６２から供給された空間周波数スペクトルSnm(a,ω,l)に空間フィルタwn(a,r,ω)を適用することで、空間周波数スペクトルを、球状マイクアレイ１１のセンサ半径aより大きな半径rの環状の仮想スピーカアレイ１３の仮想スピーカアレイ駆動信号に変換する。すなわち、次式（７）が計算されて、空間周波数スペクトルSnm(a,ω,l)が、仮想スピーカアレイ駆動信号、すなわち空間周波数スペクトルDnm(r,ω,l)に変換される。
【００８２】
【数７】
0006458738.tif 000008
【００８３】
なお、式（７）における空間フィルタwn(a,r,ω)は、例えば次式（８）に示されるフィルタとされる。
【００８４】
【数８】
0006458738.tif 000009
【００８５】
さらに、式（８）におけるBn(ka)およびRn(kr)は、それぞれ次式（９）および式（１０）に示す関数とされる。
【００８６】
【数９】
0006458738.tif 000010
【００８７】
【数１０】
0006458738.tif 000011
【００８８】
なお、式（９）および式（１０）において、JnおよびHnは、それぞれ球面ベッセル関数および第一種球面ハンケル関数を示す。また、Jn’およびHn’は、それぞれJnおよびHnの微分値を示している。
【００８９】
このように空間周波数スペクトルに空間フィルタを用いたフィルタ処理を施すことで、球状マイクアレイ１１により収音して得られた収音信号を、仮想スピーカアレイ１３で再生したときに音場が再現される仮想スピーカアレイ駆動信号へと変換することができる。
【００９０】
このようにして収音信号を、仮想スピーカアレイ駆動信号へと変換する処理は、時間周波数領域では行うことができないので、音場再現器４１は収音信号を空間周波数スペクトルに変換し、空間フィルタを適用する。
【００９１】
空間フィルタ適用部６３は、このようにして得られた空間周波数スペクトルDnm(r,ω,l)を、空間周波数合成部６４に供給する。
【００９２】
（空間周波数合成部）
空間周波数合成部６４は、次式（１１）の計算を行うことで、空間フィルタ適用部６３から供給された空間周波数スペクトルDnm(r,ω,l)の空間周波数合成を行い、時間周波数スペクトルDt(xvspk,ω,l)を得る。
【００９３】
【数１１】
0006458738.tif 000012
【００９４】
なお、式（１１）において、Ｎは球面調和関数Ynm(θp,φp)の次数を示しており、ｎは次数を示している。また、θpはセンサ方位角を示しており、φpはセンサ仰角を示しており、rは仮想スピーカアレイ１３の半径を示している。ωは時間周波数インデクスを示しており、xvspkは仮想スピーカアレイ１３を構成するスピーカを示すインデクスである。
【００９５】
空間周波数合成部６４では、仮想スピーカアレイ１３を構成する各スピーカについて、時間フレームｌごとに時間周波数の数であるΩ個の時間周波数スペクトルDt(xvspk,ω,l)が得られる。
【００９６】
空間周波数合成部６４は、このようにして得られた時間周波数スペクトルDt(xvspk,ω,l)を、逆フィルタ適用部６５に供給する。
【００９７】
（逆フィルタ生成部）
また、逆フィルタ生成器５２の逆フィルタ生成部７２は、時間周波数分析部７１から供給された時間周波数スペクトルS(x,ω,l)に基づいて逆フィルタH(xvspk,xrspk,ω)を求める。
【００９８】
時間周波数スペクトルS(x,ω,l)は、実スピーカアレイ１２から仮想スピーカアレイ１３までの伝達関数g(xvspk,xrspk,n)を時間周波数分析した結果であり、ここでは図５下段の時間周波数分析部６１で得られる時間周波数スペクトルS(p,ω,l)と区別するためにG(xvspk,xrspk,ω)と記すこととする。
【００９９】
なお、伝達関数g(xvspk,xrspk,n)、時間周波数スペクトルG(xvspk,xrspk,ω)、および逆フィルタH(xvspk,xrspk,ω)におけるxvspkは仮想スピーカアレイ１３を構成するスピーカを示すインデクスであり、xrspkは実スピーカアレイ１２を構成するスピーカを示すインデクスである。また、ｎは時間インデクスを示しており、ωは時間周波数インデクスを示している。なお、時間周波数スペクトルG(xvspk,xrspk,ω)では、時間フレームインデクスｌは省略されている。
【０１００】
伝達関数g(xvspk,xrspk,n)は、仮想スピーカアレイ１３の各スピーカの位置にマイクロホン（マイクセンサ）を置くことで予め測定されている。
【０１０１】
例えば逆フィルタ生成部７２は、測定結果から逆フィルタを求めることで仮想スピーカアレイ１３から実スピーカアレイ１２までの逆フィルタH(xvspk,xrspk,ω)を求める。すなわち、次式（１２）の計算により、逆フィルタH(xvspk,xrspk,ω)が算出される。
【０１０２】
【数１２】
0006458738.tif 000013
【０１０３】
なお、式（１２）において、ＨおよびＧは、それぞれ逆フィルタH(xvspk,xrspk,ω)および時間周波数スペクトルG(xvspk,xrspk,ω)（伝達関数g(xvspk,xrspk,n)）を行列で表したものであり、（・）−１は擬似逆行列を示す。一般的に、行列のランクが低い場合には安定した解を求めることができない。
【０１０４】
すなわち、仮想スピーカアレイ１３の半径ｒが小さいと、つまり仮想スピーカアレイ１３の中心位置（基準位置）から、仮想スピーカアレイ１３のスピーカまでの距離が短いと、各伝達関数g(xvspk,xrspk,n)の特性のばらつきが小さくなる。そうすると、行列のランクが低くなり、安定した解を求めることができなくなる。そこで、安定解を求めることが可能な球状または環状の仮想スピーカの半径ｒを予め求めておく。
【０１０５】
このとき、安定解を求めることができるように、すなわち正確な逆フィルタH(xvspk,xrspk,ω)を得ることができるように、少なくとも仮想スピーカアレイ１３の半径ｒは、球状マイクアレイ１１のセンサ半径ａよりも大きい値となるように定められるものとする。
【０１０６】
伝達関数g(xvspk,xrspk,n)から逆フィルタH(xvspk,xrspk,ω)を求めておけば、逆フィルタを用いたフィルタ処理により、仮想スピーカアレイ１３により音場を再現するための仮想スピーカアレイ駆動信号を、任意の形状の実スピーカアレイ１２の実スピーカアレイ駆動信号へと変換することができる。
【０１０７】
逆フィルタ生成部７２は、このようにして得られた逆フィルタH(xvspk,xrspk,ω)を逆フィルタ適用部６５に供給する。
【０１０８】
（逆フィルタ適用部）
逆フィルタ適用部６５は、空間周波数合成部６４から供給された時間周波数スペクトルDt(xvspk,ω,l)に、逆フィルタ生成部７２から供給された逆フィルタH(xvspk,xrspk,ω)を適用し、逆フィルタ信号Di(xrspk,ω,l)を得る。すなわち、逆フィルタ適用部６５は、次式（１３）の計算を行って、フィルタ処理により逆フィルタ信号Di(xrspk,ω,l)を算出する。この逆フィルタ信号は、音場を再現するための実スピーカアレイ駆動信号の時間周波数スペクトルである。逆フィルタ適用部６５では、実スピーカアレイ１２を構成する各スピーカについて、時間フレームｌごとに時間周波数の数であるΩ個の逆フィルタ信号Di(xrspk,ω,l)が得られる。
【０１０９】
【数１３】
0006458738.tif 000014
【０１１０】
逆フィルタ適用部６５は、このようにして得られた逆フィルタ信号Di(xrspk,ω,l)を時間周波数合成部６６に供給する。
【０１１１】
（時間周波数合成部）
時間周波数合成部６６は、次式（１４）の計算を行うことで、逆フィルタ適用部６５から供給された逆フィルタ信号Di(xrspk,ω,l)、すなわち時間周波数スペクトルの時間周波数合成を行い、出力フレーム信号d’(xrspk,n,l)を得る。
【０１１２】
【数１４】
0006458738.tif 000015
【０１１３】
なお、式（１４）におけるD’(xrspk,ω,l)は、次式（１５）により得られるものである。
【０１１４】
【数１５】
0006458738.tif 000016
【０１１５】
また、ここではIDFT（Inverse Discrete Fourier Transform）（逆離散フーリエ変換）を用いる例について説明しているが、時間周波数分析部６１で使用した変換の逆変換に相当するものを用いればよい。
【０１１６】
さらに、時間周波数合成部６６は、得られた出力フレーム信号d’(xrspk,n,l)に、窓関数wsyn(n)を乗算し、オーバーラップ加算を行うことでフレーム合成を行う。例えば、次式（１６）に示される窓関数wsyn(n)が用いられて、式（１７）の計算によりフレーム合成が行われて、出力信号d(xrspk,t)が求められる。
【０１１７】
【数１６】
0006458738.tif 000017
【０１１８】
【数１７】
0006458738.tif 000018
【０１１９】
なお、ここでは、時間周波数分析部６１で用いた窓関数と同じものを用いているが、ハミング窓などのその他の窓の場合は矩形窓で構わない。
【０１２０】
また、式（１７）において、dprev(xrspk,n+lN)およびdcurr(xrspk,n+lN)は、どちらも出力信号d(xrspk,t)を示しているが、dprev(xrspk,n+lN)は更新前の値を示し、dcurr(xrspk,n+lN)は更新後の値を示している。
【０１２１】
時間周波数合成部６６は、このようにして得られた出力信号d(xrspk,t)を、実スピーカアレイ駆動信号として音場再現器４１の出力とする。
【０１２２】
以上のように、音場再現器４１によれば、より正確に音場を再現することができる。
【０１２３】
〈実スピーカアレイ駆動信号生成処理の説明〉
次に、以上において説明した音場再現器４１により行われる処理の流れについて説明する。音場再現器４１は、伝達関数と収音信号が供給されると、収音信号を実スピーカアレイ駆動信号へと変換して出力する実スピーカアレイ駆動信号生成処理を行う。
【０１２４】
以下、図６のフローチャートを参照して音場再現器４１による実スピーカアレイ駆動信号生成処理について説明する。なお、逆フィルタ生成器５２による逆フィルタの生成は予め行われてもよいが、ここでは実スピーカアレイ駆動信号の生成時に逆フィルタが生成されるものとして説明を続ける。
【０１２５】
ステップＳ１１において、時間周波数分析部６１は球状マイクアレイ１１から供給された収音信号s(p,t)の時間周波数情報を分析する。
【０１２６】
具体的には、時間周波数分析部６１は収音信号s(p,t)に対して時間フレーム分割を行い、その結果得られた入力フレーム信号sfr(p,n,l)に窓関数wana(n)を乗算し、窓関数適用信号sw(p,n,l)を算出する。
【０１２７】
また、時間周波数分析部６１は、窓関数適用信号sw(p,n,l)に対して時間周波数変換を行い、その結果得られた時間周波数スペクトルS(p,ω,l)を空間周波数分析部６２に供給する。すなわち、式（４）の計算が行われて時間周波数スペクトルS(p,ω,l)が算出される。
【０１２８】
ステップＳ１２において、空間周波数分析部６２は、時間周波数分析部６１から供給された時間周波数スペクトルS(p,ω,l)に対して空間周波数変換を行い、その結果得られた空間周波数スペクトルSnm(a,ω,l)を空間フィルタ適用部６３に供給する。
【０１２９】
具体的には、空間周波数分析部６２は式（５）を計算することで、時間周波数スペクトルS(p,ω,l)を空間周波数スペクトルSnm(a,ω,l)に変換する。
【０１３０】
ステップＳ１３において、空間フィルタ適用部６３は、空間周波数分析部６２から供給された空間周波数スペクトルSnm(a,ω,l)に空間フィルタwn(a,r,ω)を適用する。
【０１３１】
すなわち、空間フィルタ適用部６３は式（７）を計算することで、空間周波数スペクトルSnm(a,ω,l)に対して、空間フィルタwn(a,r,ω)を用いたフィルタ処理を施し、その結果得られた空間周波数スペクトルDnm(r,ω,l)を空間周波数合成部６４に供給する。
【０１３２】
ステップＳ１４において、空間周波数合成部６４は、空間フィルタ適用部６３から供給された空間周波数スペクトルDnm(r,ω,l)の空間周波数合成を行い、その結果得られた時間周波数スペクトルDt(xvspk,ω,l)を逆フィルタ適用部６５に供給する。すなわち、ステップＳ１４では、式（１１）の計算が行われて、時間周波数スペクトルDt(xvspk,ω,l)が求められる。
【０１３３】
ステップＳ１５において、時間周波数分析部７１は、供給された伝達関数g(xvspk,xrspk,n)の時間周波数情報を分析する。具体的には、時間周波数分析部７１は、伝達関数g(xvspk,xrspk,n)に対してステップＳ１１における処理と同様の処理を行い、その結果得られた時間周波数スペクトルG(xvspk,xrspk,ω)を逆フィルタ生成部７２に供給する。
【０１３４】
ステップＳ１６において、逆フィルタ生成部７２は、時間周波数分析部７１から供給された時間周波数スペクトルG(xvspk,xrspk,ω)に基づいて逆フィルタH(xvspk,xrspk,ω)を算出し、逆フィルタ適用部６５に供給する。例えばステップＳ１６では、式（１２）の計算が行われ、逆フィルタH(xvspk,xrspk,ω)が算出される。
【０１３５】
ステップＳ１７において、逆フィルタ適用部６５は、空間周波数合成部６４から供給された時間周波数スペクトルDt(xvspk,ω,l)に対して、逆フィルタ生成部７２から供給された逆フィルタH(xvspk,xrspk,ω)を適用し、その結果得られた逆フィルタ信号Di(xrspk,ω,l)を時間周波数合成部６６に供給する。例えば、ステップＳ１７では式（１３）の計算が行われ、フィルタ処理により逆フィルタ信号Di(xrspk,ω,l)が算出される。
【０１３６】
ステップＳ１８において、時間周波数合成部６６は、逆フィルタ適用部６５から供給された逆フィルタ信号Di(xrspk,ω,l)の時間周波数合成を行う。
【０１３７】
具体的には、時間周波数合成部６６は式（１４）の計算を行って、逆フィルタ信号Di(xrspk,ω,l)から出力フレーム信号d’(xrspk,n,l)を算出する。さらに時間周波数合成部６６は、出力フレーム信号d’(xrspk,n,l)に窓関数wsyn(n)を乗算して式（１７）の計算を行い、フレーム合成により出力信号d(xrspk,t)を算出する。時間周波数合成部６６は、このようにして得られた出力信号d(xrspk,t)を、実スピーカアレイ駆動信号として実スピーカアレイ１２に出力し、実スピーカアレイ駆動信号生成処理は終了する。
【０１３８】
以上のようにして音場再現器４１は、空間フィルタを用いたフィルタ処理により、収音信号から仮想スピーカアレイ駆動信号を生成し、さらに仮想スピーカアレイ駆動信号に対して逆フィルタを用いたフィルタ処理により実スピーカアレイ駆動信号を生成する。
【０１３９】
音場再現器４１では、球状マイクアレイ１１のセンサ半径ａよりも大きい半径ｒの仮想スピーカアレイ１３の仮想スピーカアレイ駆動信号を生成し、得られた仮想スピーカアレイ駆動信号を逆フィルタを用いて実スピーカアレイ駆動信号に変換することで、実スピーカアレイ１２の形状がどのような形状であっても、より正確に音場を再現することができる。
【０１４０】
〈第２の実施の形態〉
〈音場再現システムの構成例〉
なお、以上においては、収音信号を実スピーカアレイ駆動信号へと変換する処理を１つの装置が実行する例について説明したが、いくつかの装置からなる音場再現システムにより、収音信号を実スピーカアレイ駆動信号へと変換する処理が行われるようにしてもよい。
【０１４１】
そのような音場再現システムは、例えば図７に示すように構成される。なお、図７において、図３または図５における場合と対応する部分には同一の符号を付してあり、その説明は省略する。
【０１４２】
図７に示す音場再現システム１０１は、駆動信号生成器１１１および逆フィルタ生成器５２から構成される。逆フィルタ生成器５２には、図５における場合と同様に、時間周波数分析部７１および逆フィルタ生成部７２が設けられている。
【０１４３】
また、駆動信号生成器１１１は、無線により相互に通信を行って各種の情報等の授受を行う送信器１２１および受信器１２２から構成される。特に送信器１２１は、球面波（音声）の収音が行われる実空間に配置されており、受信器１２２は、収音された音声を再生する再現空間に配置されている。
【０１４４】
送信器１２１は、球状マイクアレイ１１、時間周波数分析部６１、空間周波数分析部６２、および通信部１３１を有している。通信部１３１は、アンテナなどからなり、空間周波数分析部６２から供給された空間周波数スペクトルSnm(a,ω,l)を、無線通信により受信器１２２に送信する。
【０１４５】
また、受信器１２２は、通信部１３２、空間フィルタ適用部６３、空間周波数合成部６４、逆フィルタ適用部６５、時間周波数合成部６６、および実スピーカアレイ１２を有している。通信部１３２は、アンテナなどからなり、無線通信により通信部１３１から送信されてきた空間周波数スペクトルSnm(a,ω,l)を受信して、空間フィルタ適用部６３に供給する。
【０１４６】
〈音場再現処理の説明〉
次に、図８のフローチャートを参照して、図７に示した音場再現システム１０１により行われる音場再現処理について説明する。
【０１４７】
ステップＳ４１において、球状マイクアレイ１１は、実空間において音声を収音し、その結果得られた収音信号を時間周波数分析部６１に供給する。
【０１４８】
収音信号が得られると、その後、ステップＳ４２およびステップＳ４３の処理が行われるが、これらの処理は図６のステップＳ１１およびステップＳ１２の処理と同様であるので、その説明は省略する。但し、ステップＳ４３では、空間周波数分析部６２は、得られた空間周波数スペクトルSnm(a,ω,l)を通信部１３１に供給する。
【０１４９】
ステップＳ４４において、通信部１３１は、空間周波数分析部６２から供給された空間周波数スペクトルSnm(a,ω,l)を、無線通信により受信器１２２に送信する。
【０１５０】
ステップＳ４５において、通信部１３２は、無線通信により通信部１３１から送信されてきた空間周波数スペクトルSnm(a,ω,l)を受信して、空間フィルタ適用部６３に供給する。
【０１５１】
空間周波数スペクトルが受信されると、その後、ステップＳ４６乃至ステップＳ５１の処理が行われるが、これらの処理は図６のステップＳ１３乃至ステップＳ１８の処理と同様であるので、その説明は省略する。但し、ステップＳ５１では、時間周波数合成部６６は、得られた実スピーカアレイ駆動信号を実スピーカアレイ１２に供給する。
【０１５２】
ステップＳ５２において、実スピーカアレイ１２は、時間周波数合成部６６から供給された実スピーカアレイ駆動信号に基づいて音声を再生し、音場再現処理は終了する。このようにして実スピーカアレイ駆動信号に基づいて音声が再生されると、再現空間において実空間の音場が再現される。
【０１５３】
以上のようにして音場再現システム１０１は、空間フィルタを用いたフィルタ処理により、収音信号から仮想スピーカアレイ駆動信号を生成し、さらに仮想スピーカアレイ駆動信号に対して逆フィルタを用いたフィルタ処理により実スピーカアレイ駆動信号を生成する。
【０１５４】
このとき、球状マイクアレイ１１のセンサ半径ａよりも大きい半径ｒの仮想スピーカアレイ１３の仮想スピーカアレイ駆動信号を生成し、得られた仮想スピーカアレイ駆動信号を逆フィルタを用いて実スピーカアレイ駆動信号に変換することで、実スピーカアレイ１２の形状がどのような形状であっても、より正確に音場を再現することができる。
【０１５５】
ところで、上述した一連の処理は、ハードウェアにより実行することもできるし、ソフトウェアにより実行することもできる。一連の処理をソフトウェアにより実行する場合には、そのソフトウェアを構成するプログラムが、コンピュータにインストールされる。ここで、コンピュータには、専用のハードウェアに組み込まれているコンピュータや、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のコンピュータなどが含まれる。
【０１５６】
図９は、上述した一連の処理をプログラムにより実行するコンピュータのハードウェアの構成例を示すブロック図である。
【０１５７】
コンピュータにおいて、ＣＰＵ（Central Processing Unit）５０１，ＲＯＭ（Read Only Memory）５０２，ＲＡＭ（Random Access Memory）５０３は、バス５０４により相互に接続されている。
【０１５８】
バス５０４には、さらに、入出力インターフェース５０５が接続されている。入出力インターフェース５０５には、入力部５０６、出力部５０７、記録部５０８、通信部５０９、及びドライブ５１０が接続されている。
【０１５９】
入力部５０６は、キーボード、マウス、マイクロホン、撮像素子などよりなる。出力部５０７は、ディスプレイ、スピーカなどよりなる。記録部５０８は、ハードディスクや不揮発性のメモリなどよりなる。通信部５０９は、ネットワークインターフェースなどよりなる。ドライブ５１０は、磁気ディスク、光ディスク、光磁気ディスク、又は半導体メモリなどのリムーバブルメディア５１１を駆動する。
【０１６０】
以上のように構成されるコンピュータでは、ＣＰＵ５０１が、例えば、記録部５０８に記録されているプログラムを、入出力インターフェース５０５及びバス５０４を介して、ＲＡＭ５０３にロードして実行することにより、上述した一連の処理が行われる。
【０１６１】
コンピュータ（ＣＰＵ５０１）が実行するプログラムは、例えば、パッケージメディア等としてのリムーバブルメディア５１１に記録して提供することができる。また、プログラムは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の伝送媒体を介して提供することができる。
【０１６２】
コンピュータでは、プログラムは、リムーバブルメディア５１１をドライブ５１０に装着することにより、入出力インターフェース５０５を介して、記録部５０８にインストールすることができる。また、プログラムは、有線または無線の伝送媒体を介して、通信部５０９で受信し、記録部５０８にインストールすることができる。その他、プログラムは、ＲＯＭ５０２や記録部５０８に、あらかじめインストールしておくことができる。
【０１６３】
なお、コンピュータが実行するプログラムは、本明細書で説明する順序に沿って時系列に処理が行われるプログラムであっても良いし、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで処理が行われるプログラムであっても良い。
【０１６４】
また、本技術の実施の形態は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。
【０１６５】
例えば、本技術は、１つの機能をネットワークを介して複数の装置で分担、共同して処理するクラウドコンピューティングの構成をとることができる。
【０１６６】
また、上述のフローチャートで説明した各ステップは、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０１６７】
さらに、１つのステップに複数の処理が含まれる場合には、その１つのステップに含まれる複数の処理は、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０１６８】
また、本明細書中に記載された効果はあくまで例示であって限定されるものではなく、他の効果があってもよい。
【０１６９】
さらに、本技術は、以下の構成とすることも可能である。
【０１７０】
（１）
球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成部と、
前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成部と
を備える音場再現装置。
（２）
前記第１の駆動信号生成部は、前記収音信号から得られた空間周波数スペクトルに対して空間フィルタを用いたフィルタ処理を施すことで、前記収音信号を前記仮想スピーカアレイの駆動信号に変換する
（１）に記載の音場再現装置。
（３）
前記収音信号から得られた時間周波数スペクトルを前記空間周波数スペクトルに変換する空間周波数分析部をさらに備える
（２）に記載の音場再現装置。
（４）
前記第２の駆動信号生成部は、前記実スピーカアレイから前記仮想スピーカアレイまでの伝達関数に基づく逆フィルタを用いて、前記仮想スピーカアレイの駆動信号に対してフィルタ処理を施すことで、前記仮想スピーカアレイの駆動信号を前記実スピーカアレイの駆動信号に変換する
（１）乃至（３）の何れか一項に記載の音場再現装置。
（５）
前記仮想スピーカアレイは球状または環状のスピーカアレイである
（１）乃至（４）の何れか一項に記載の音場再現装置。
（６）
球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成ステップと、
前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成ステップと
を含む音場再現方法。
（７）
球状または環状のマイクアレイが収音することで得られた収音信号を、前記マイクアレイの第１の半径よりも大きい第２の半径を有する仮想スピーカアレイの駆動信号に変換する第１の駆動信号生成ステップと、
前記仮想スピーカアレイの駆動信号を、前記仮想スピーカアレイが囲む空間の内側または外側に配置された実スピーカアレイの駆動信号に変換する第２の駆動信号生成ステップと
を含む処理をコンピュータに実行させるプログラム。
【符号の説明】
【０１７１】
１１  球状マイクアレイ，  １２  実スピーカアレイ，  １３  仮想スピーカアレイ，  ４１  音場再現器，  ５１  駆動信号生成器，  ５２  逆フィルタ生成器，  ６１  時間周波数分析部，  ６２  空間周波数分析部，  ６３  空間フィルタ適用部，  ６４  空間周波数合成部，  ６５  逆フィルタ適用部，  ６６  時間周波数合成部，  ７１  時間周波数分析部，  ７２  逆フィルタ生成部，  １３１  通信部，  １３２  通信部
【図１】
0006458738.tif 000019
【図２】
0006458738.tif 000020
【図３】
0006458738.tif 000021
【図４】
0006458738.tif 000022
【図５】
0006458738.tif 000023
【図６】
0006458738.tif 000024
【図７】
0006458738.tif 000025
【図８】
0006458738.tif 000026
【図９】
0006458738.tif 000027