(57)【要約】    【構成】  局所再生装置（１０）は、複数のスピーカ（１４１、１４２、…、１４Ｎ）を一定間隔でｘ軸方向の一直線上に配置した直線スピーカアレイ（１２）を含み、駆動装置（１６）は、ブライト領域が「１」でダーク領域が「０」の矩形窓で、直線スピーカアレイから所定距離ｙｂの聴取位置における音圧ｐをモデル化し、ｘ軸方向の空間フーリエ変換を行なう。必要なら、ブライト領域の位置ｘｂを変更して、フィルタ関数Ｆを求める。この空間フィルタ関数ｆと、音源ｓで各スピーカの駆動信号ｄを生成して、各スピーカを駆動する。    【効果】  空間フィルタ関数を解析的に導出できるので、煩瑣な計算なしに、安定した駆動信号を算出でき、精度のよい局所再生を実現できる。さらに、局所再生の位置や幅を変えた複数の局所再生を同時に再生することにより、マルチ局所音響再生を精度よく実現できる。【選択図】  図１
2015231087.tif 000002
【特許請求の範囲】
【請求項１】
複数のスピーカを互いに間隔を隔てて直線状に配置した直線スピーカアレイを用いる局所音響再生装置であって、
再生位置における音圧を任意幅の矩形窓でモデル化して音圧モデルを求めるモデル化手段、
前記音圧モデルに基づいて前記複数のスピーカの各々の空間フィルタ関数を導出する導出手段、および
音源および前記空間フィルタ関数に基づいて前記複数のスピーカのための駆動信号を生成して、各スピーカを駆動する駆動手段を備える、局所音響再生装置。
【請求項２】
前記モデル化手段は位置が異なる複数の矩形窓でモデル化した音圧モデルを求め、前記導出手段は当該複数の矩形窓でモデル化した音圧モデルに従って前記空間フィルタ関数を導出する、請求項１記載の局所音響再生装置。
【請求項３】
複数のスピーカを互いに間隔を隔てて直線状に配置した直線スピーカアレイを用いる局所音響再生装置のコンピュータによって実行される局所音響再生プログラムであって、前記局所音響再生プログラムは、前記コンピュータを、
再生位置における音圧を任意幅の矩形窓でモデル化して音圧モデルを求めるモデル化手段、
前記音圧モデルに基づいて前記複数のスピーカの各々の空間フィルタ関数を導出する導出手段、および
音源および前記空間フィルタ関数に基づいて前記複数のスピーカのための駆動信号を生成して、各スピーカを駆動する駆動手段
として機能させる、局所音響再生プログラム。
【請求項４】
前記モデル化手段は位置が異なる複数の矩形窓でモデル化した音圧モデルを求め、前記導出手段は当該複数の矩形窓でモデル化した音圧モデルに従って前記空間フィルタ関数を導出する、請求項３記載の局所音響再生プログラム。
【請求項５】
直線状に配置された多数のスピーカから成る直線スピーカアレイを用いて該直線スピーカアレイに対して任意の位置で任意の幅の局所に音響を再生する局所音響再生装置であって、
上記直線スピーカアレイのスピーカの並びに平行に音響が聞こえる幅を規定する幅規定手段と、
上記直線スピーカアレイのスピーカの並びに対して平行な位置で音響が聞こえる位置を規定する位置規定手段と、
上記幅規定手段および位置規定手段によって設定される空間フィルタ関数と音源とによって駆動信号を生成して上記直線状スピーカアレイの各スピーカを駆動するスピーカ駆動手段とを備える、局所音響再生装置。
【請求項６】
上記直線スピーカアレイからの距離を規定する距離規定手段をさらに備え、
上記空間フィルタ関数は、さらに上記距離規定手段によって設定される、請求項５記載の局所音響再生装置。
【請求項７】
直線状に配置された複数のスピーカから成る直線スピーカアレイを用いて該直線スピーカアレイに対して任意の位置で任意の幅の局所に音響を再生する局所音響再生装置のコンピュータによって実行される局所音響再生プログラムであって、前記局所音響再生プログラムは、前記コンピュータを、
上記直線スピーカアレイのスピーカの並びに平行に音響が聞こえる幅を規定する幅規定手段、
上記直線スピーカアレイのスピーカの並びに対して平行な位置で音響が聞こえる位置を規定する位置規定手段、および
上記幅規定手段および位置規定手段によって設定される空間フィルタ関数と音源とによって駆動信号を生成して上記直線状スピーカアレイの各スピーカを駆動するスピーカ駆動手段として機能させる、局所音響再生プログラム。
【請求項８】
さらに、前記コンピュータを、前記直線スピーカアレイからの距離を規定する距離規定手段として機能させ、
上記空間フィルタ関数は、さらに上記距離規定手段によって設定される、請求項７記載の局所音響再生プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は局所音響再生装置およびプログラムに関し、特に直線状に配置された多数のスピーカから成る直線スピーカアレイに対して任意の位置の任意の幅で局所的に音響を再生する、新規な局所音響再生装置およびプログラムに関する。
【背景技術】
【０００２】
多数のスピーカを用いた局所音響再生には大きく分けて２つの方式がある。１つは、スピーカアレイの指向性を制御するビームフォーミングに基づく方法であり、たとえば特許文献１などに開示されている。もう１つが、各スピーカで制御点の音圧を直接制御する多点制御に基づく方法であり、たとえば非特許文献１および非特許文献２などにおいて、近年いくつかの方法が提案されている。ビームフォーミングに基づく方法は、目的方向の指向性を最大にする方法であるのに対し、多点制御に基づく方法は制御点間の音響エネルギ差を最大にする方法となる。そのため、音が聞こえないエリアを形成するようなシステムの場合は後者の方法が性能が高い。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開2012-147414 [H04R 3/12, 3/00, 1/40, 1/34]
【非特許文献】
【０００４】
【非特許文献１】J.-W. Choi and Y.-H. Kim, “Generation of an acoustically bright zone with an illuminated region using multiple sources,” J. Acoust.Soc. Am., vol. 111, no. 4, pp. 1695- 1700, Apr. 2002.
【非特許文献２】M. Shin, S. Q. Lee, F. M. Fazi, P. A. Nelson, D. Kim, S. Wang, K. H. Park, and J. Seo, “Maximization of acoustic energy difference between two spaces,” J. Acoust. Soc. Am., vol. 128, no. 1, pp. 121-131, July 2010.
【発明の開示】
【発明が解決しようとする課題】
【０００５】
コントラスト最大化法を採用する非特許文献１では、各スピーカと制御点間の空間相関行列の逆行列演算を用いるが、音響逆問題は極めて悪条件であるため、計算により得られた解は不安定となる。この問題を解決するためには、空間相関行列の正則化が必要であるが、適切な正則化パラメータにも限界があり、さらにはパラメータ決定のためには各パラメータ毎の反復計算が必要となる。
【０００６】
それに対して、非特許文献２のようなエネルギ差最大化法(Acoustic energy difference maximization: ＥＤＭ)では、空間相関行列の逆行列ではなく、固有ベクトルの計算に基づくため、解の不安定という問題はなく、非特許文献１のものより高い精度で局所音響再生を実現している。
【０００７】
しかし、この方法ではチューニング係数を周波数毎に決める必要があるため、正則化問題と同様、適切な係数を決めるための反復計算が必要となる。
【０００８】
それゆえに、この発明の主たる目的は、新規な、局所音響再生装置およびプログラムを提供することである。
【０００９】
この発明の他の目的は、煩瑣な計算を不要として、精度の高い局所音響再生を可能にする、局所音響再生装置およびプログラムを提供することである。
【課題を解決するための手段】
【００１０】
本発明は、上記の課題を解決するために、以下の構成を採用した。なお、括弧内の参照符号および補足説明等は、本発明の理解を助けるために後述する実施の形態との対応関係を示したものであって、本発明を何ら限定するものではない。
【００１１】
第１の発明は、複数のスピーカを互いに間隔を隔てて直線状に配置した直線スピーカアレイを用いる局所音響再生装置であって、再生位置における音圧を任意幅の矩形窓でモデル化して音圧モデルを求めるモデル化手段、音圧モデルに基づいて複数のスピーカの各々の空間フィルタ関数を導出する導出手段、および音源および空間フィルタ関数に基づいて複数のスピーカのための駆動信号を生成して、各スピーカを駆動する駆動手段を備える、局所音響再生装置である。
【００１２】
第１の発明では、局所音響再生装置（１０：実施例において相当する部分を例示する参照符号。以下、同様。）は、複数のスピーカ（１４１、１４２、…、１４Ｎ）を互いに間隔（Δｘ）を隔てて直線状に配置した直線スピーカアレイ（１２）を用いる。モデル化手段（１６）は、たとえば、この直線スピーカアレイ（１２）に平行な再生位置（ｙｂ）における音圧（ｐ）を任意幅（ｌｂ）の矩形窓でモデル化して音圧モデル（たとえば、数［９］）を求める。導出手段（１６）は、この音圧モデル（ｐ）に基づいて、たとえば空間フーリエ変換等の手法を用いて、複数のスピーカの各々の空間フィルタ関数（ｆｉ：たとえば数［１２］）を解析的に導出する。そして、駆動手段（１６）は、音源（ｓ）および空間フィルタ関数（ｆｉ）に基づいて複数のスピーカ（１４ｉ）のための駆動信号（ｄｉ）を生成し、その駆動信号（ｄ１、ｄ２、…、ｄＮ）で各スピーカ（１４１、１４２、…、１４Ｎ）を駆動する。
【００１３】
第１の発明によれば、空間フィルタ関数を解析的に導出できるので、煩瑣な計算は不要で、しかも高い精度でブライト領域およびダーク領域を設定して、局所音響再生を行うことができる。
【００１４】
第２の発明は、第１の発明に従属し、モデル化手段は位置が異なる複数の矩形窓でモデル化した音圧モデルを求め、導出手段は当該複数の矩形窓でモデル化した音圧モデルに従って空間フィルタ関数を導出する、局所音響再生装置である。
【００１５】
第２の発明では、モデル化手段（１６）は、位置が異なる複数の矩形窓で音圧モデルを求め、導出手段（１６）は、各音圧モデルに従って空間フィルタ関数を解析的に導出する。
【００１６】
第２の発明によれば、ブライト領域が異なる、いわゆるマルチ局所音響再生が可能になる。
【００１７】
第３の発明は、複数のスピーカを互いに間隔を隔てて直線状に配置した直線スピーカアレイを用いる局所音響再生装置のコンピュータによって実行される局所音響再生プログラムであって、局所音響再生プログラムは、コンピュータを、再生位置における音圧を任意幅の矩形窓でモデル化して音圧モデルを求めるモデル化手段、音圧モデルに基づいて複数のスピーカの各々の空間フィルタ関数（ｆ）を導出する導出手段、および音源および空間フィルタ関数に基づいて複数のスピーカのための駆動信号を生成して、各スピーカを駆動する駆動手段として機能させる、局所音響再生プログラムである。
【００１８】
第３の発明においても、第１の発明と同様の効果が期待できる。
【００１９】
第４の発明は、第３の発明に従属し、モデル化手段は第１が異なる複数の矩形窓でモデル化した音圧モデルを求め、導出手段は当該複数の矩形窓でモデル化した音圧モデルに従って空間フィルタ関数を導出する、局所音響再生プログラムである。
【００２０】
第４の発明においても、第２の発明と同様の効果が期待できる。
【００２１】
第５の発明は、直線状に配置された多数のスピーカから成る直線スピーカアレイを用いて該直線スピーカアレイに対して任意の位置で任意の幅の局所に音響を再生する局所音響再生装置であって、直線スピーカアレイのスピーカの並びに平行に音響が聞こえる幅を規定する幅規定手段と、直線スピーカアレイのスピーカの並びに対して平行な位置で音響が聞こえる位置を規定する位置規定手段と、幅規定手段および位置規定手段によって設定される空間フィルタ関数と音源とによって駆動信号を生成して直線状スピーカアレイの各スピーカを駆動するスピーカ駆動手段と、を備える、局所音響再生装置である。
【００２２】
第５の発明では、局所音響再生装置（１０）は、複数のスピーカ（１４１、１４２、…、１４Ｎ）を互いに間隔（Δｘ）を隔てて直線状に配置した直線スピーカアレイ（１２）を用いる。幅規定手段（１６）は、たとえばブライト領域を設定すべき、直線スピーカアレイのスピーカの並び（ｘ方向）に平行に音響が聞こえる幅（ｌｂ）を規定する。そして、位置規定手段（１６）は、たとえばブライト領域を設定すべき、直線スピーカアレイのスピーカの並び（ｘ方向）に平行に音響が聞こえる位置（ｘｂ）を規定する。スピーカ駆動手段（１６）は、幅規定手段および位置規定手段によって設定される空間フィルタ関数（ｆｉ）およびたとえば音源ユニット（１８）から与えられる音源（ｓ）に基づいて駆動信号（ｄｉ）を生成し、その駆動信号で各スピーカを駆動する。
【００２３】
第５の発明によれば、幅および位置を規定することによって空間フィルタ関数を解析的に導出できるので、煩瑣な計算は不要で、しかも高い精度でブライト領域およびダーク領域を設定して、局所音響再生を行うことができる。
【００２４】
第６の発明は、第５の発明に従属し、直線スピーカアレイからの距離を規定する距離規定手段をさらに備え、空間フィルタ関数はさらに距離規定手段によって設定される、局所音響再生装置である。
【００２５】
第６の発明では、距離規定手段（１６）は、直線スピーカアレイに直交する方向での直線スピーカアレイからの距離（ｙｂ）を規定する。そして、空間フィルタ関数はさらにその距離規定手段によって設定され、駆動信号は、その空間フィルタ関数と音源とによって生成される。
【００２６】
第６の発明によれば、距離規定手段によって距離を規定するので、精度がより高くなる。
【００２７】
第７の発明は、直線状に配置された複数のスピーカから成る直線スピーカアレイを用いて該直線スピーカアレイに対して任意の位置で任意の幅の局所に音響を再生する局所音響再生装置のコンピュータによって実行される局所音響再生プログラムであって、局所音響再生プログラムは、コンピュータを、直線スピーカアレイのスピーカの並びに平行に音響が聞こえる幅を規定する幅規定手段、直線スピーカアレイのスピーカの並びに対して平行な位置で音響が聞こえる位置を規定する位置規定手段、および幅規定手段および位置規定手段によって設定される空間フィルタ関数と音源とによって駆動信号を生成して直線状スピーカアレイの各スピーカを駆動するスピーカ駆動手段として機能させる、局所音響再生プログラムである。
【００２８】
第７の発明においても、第５の発明と同様の効果が期待できる。
【００２９】
第８の発明は、第７の発明に従属し、さらに、コンピュータを、直線スピーカアレイからの距離を規定する距離規定手段として機能させ、駆動信号は、空間フィルタ関数に加えて、距離規定手段によって設定される、局所音響再生プログラムである。
【００３０】
第８の発明においても、第６の発明と同様の効果が期待できる。
【発明の効果】
【００３１】
この発明によれば、空間フィルタ関数を解析的に導出できるので、煩瑣な計算なしに、精度の高い局所音響再生が行える。
【００３２】
この発明の上述の目的，その他の目的，特徴および利点は、図面を参照して行う。以下の実施例の詳細な説明から一層明らかとなろう。
【図面の簡単な説明】
【００３３】
【図１】図１はこの発明の一実施例の局所音響再生装置の構成を示す概略図である。
【図２】図２はスペクトル除算法（ＳＤＭ：Spectrum Division Method）に基づく無限直線二次音源を用いた音場再現の状態を示す概略図である。
【図３】図３は波数領域空間フィルタリングによる無限直線二次音源を用いた局所音響再生を示す概略図である。
【図４】図４はこの発明の他の実施例の局所音響再生装置の構成を示す概略図である。
【図５】図５は有限個の直線スピーカアレイを用いる場合の駆動信号を計算するための矩形窓を示す概略図である。
【図６】図６はＥＤＭにおけるスピーカと局所音響再生制御(ブライト領域とダーク領域)の配置を示す概略図である。
【図７】図７は計算機シミュレーションでの時空間平均音圧レベルの結果をＥＤＭとの比較で示すグラフである。
【図８】図８は計算機シミュレーションでのＢＤＲの結果をＥＤＭとの比較で示すグラフであり、上から順に５００Ｈｚ、１ｋＨｚ、２ｋＨｚの結果を示し、いずれも左側がＥＤＭで、右側が実施例のものである。
【図９】図９は計算機シミュレーションでの時間平均音圧レベルの結果をＥＤＭとの比較で示すグラフであり、上から順に３ｋＨｚ、４ｋＨｚ、５ｋＨｚの結果を示し、いずれも左側がＥＤＭで、右側が実施例のものである。
【発明を実施するための最良の形態】
【００３４】
図１はこの発明の一実施例の局所音響再生装置を示す概略図である。この実施例の局所音響再生装置１０は、直線状にスピーカを多数配置した直線スピーカアレイ１２を用いる。直線スピーカアレイ１２は、複数（実施例では、Ｎ=６４個）のスピーカ１４１、１４２、…、１４Ｎを所定間隔、たとえば５ｃｍ毎に直線状に配置したものであり、各スピーカ１４は一例として１インチ（２．５ｃｍ）の口径を有する。ただし、これらの数値は適宜変更可能である。
【００３５】
たとえばコンピュータからなる駆動装置(コンピュータ)１６は、適宜の記憶媒体（メモリ、ＨＤＤ、光ディスクやネットワークなどの音源ユニット１８から与えられる音源ｓ（ω）に基づく駆動信号ｄ１、ｄ２、…、ｄＮを作成し、これらの駆動信号によって、各スピーカ１４１、１４２、…、１４Ｎが駆動される。
【００３６】
なお、この発明が向けられる局所音響再生において、「音響」とは、音声、音楽その他、音に関する音信号を意味しており、上述の音源１８から音信号を駆動装置１６に与える。
【００３７】
ただし、この実施例では、直線スピーカアレイ１２の前に存在する聴取者Ｐ１、Ｐ２、Ｐ３、…のうち、一例として、聴取者Ｐ２にのみ音源ｓを聞かせ、他の聴取者には聞かせない、という局所音響再生を行う。ここで、音源ｓを聞かせる領域を図１では「Bright（ブライト）」と表現し、聞かせない領域を「Dark（ダーク）」と表現する。ブライト領域では、その領域に存在する聴取者は、１人でも複数人でも、すべての聴取者が音源を聴取できる。
【００３８】
このようなブライト領域およびダーク領域を形成するために、この実施例では、駆動装置すなわちコンピュータ１６は、後に説明するように、スピーカ１４毎の空間フィルタ関数Ｆｉを計算し、その空間フィルタ関数と音源ｓ（ω）とを乗算することによって、各スピーカのための駆動信号ｄｉを生成する。そして、生成した各スピーカ用の駆動信号ｄ１、ｄ２、…、ｄＮで各スピーカ１４１、１４２、…、１４Ｎを駆動する。
【００３９】
図２は、図１実施例の前提としてのＳＤＭによる音場再現を示す概略図である。ただし、図１実施例は有限個のスピーカ１４１、１４２、…、１４Ｎからなる直線スピーカアレイ１２を用いるが、以下の説明では、計算のために、無限に連続した直線スピーカ（無限直線二次音源）を想定して音場再現を考える。
【００４０】
ただし、ＳＤＭについては、参考文献１（J. Ahrens, and S. Spors,“Sound field reproduction using planar and linear arrays of loudspeakers,”IEEE Trans. Audio, Speech, Lang. Process., vol. 18, no. 8, pp. 2038-2050, Nov. 2010.）に詳しい。
【００４１】
無限直線二次音源はｘ軸上（ｙ＝０）にあると仮定し、ｚ軸は省略する。周波数をｆとし、角周波数をω＝２πｆとし、無限直線二次音源の位置ｘ０＝［ｘ，０］Ｔにおける駆動信号をｄ（ｘ０，ω）とすると，位置ｘ＝［ｘ，ｙ］Ｔでの音圧ｐ（ｘ，ω）は、数［１］で与えられる。
【００４２】
なお、以下の各式において、フーリエ変換前の関数をアルファベットの小文字で表現し、フーリエ変換後の関数を、「^(ハット)」を用いる代わりに、大文字で表現する。
【００４３】
【数１】
2015231087.tif 000003
【００４４】
ここで、
【００４５】
【数２】
2015231087.tif 000004
【００４６】
は、無限直線二次音源からの音声を伝達する特性を示す３次元グリーン関数であり、ｊ＝√−１であり、ｋ＝ω／ｃは波数を示し、ｃは音速を示す。ただし、ここでは、反射などは考慮していない。
【００４７】
数［１］は、無限直線二次音源の前方の受音位置ｙｒｅｆにおける各スピーカからの音圧を合計した音圧ｐを示し、この数［１］をｘ軸方向に空間フーリエ変換すると、畳み込みの定理より、数３が得られる。
【００４８】
【数４】
2015231087.tif 000005
【００４９】
ここで、ｋｘはｘ軸方向の空間周波数、Ｄ(ｋｘ,ω)は駆動信号のフーリエ変換、Ｇ（ｋｘ，ｙ，ω）は３次元グリーン関数ｇ３Ｄ（ｘ−ｘ０，ω）の空間フーリエ変換であり、数［４］で表わされる。
【００５０】
ただし、空間フーリエ変換については、たとえば、参考文献２（E. G. Williams, Fourier Acoustics: Sound Radiation and Nearfield Acoustic Holography, London: Academic Press, 1999.）を参照されたい。
【００５１】
空間フーリエ変換することによって、受音位置における各方向からの波の成分量が得られる。
【００５２】
【数４】
2015231087.tif 000006
【００５３】
Ｈ０（２）およびｋ０はそれぞれ０次の第２種ハンケル関数および０次の変形ベッセル関数を示す。
【００５４】
したがって、図２に示すようにｘ軸と平行に連続した受音位置（ｙ＝ｙｒｅｆ）における音圧の空間フーリエ変換をＰ（ｋｘ，ｙｒｅｆ，ω）とすると、ＳＤＭでは、波数領域における無限直線二次音源の駆動信号Ｄ（ｋｘ，ω）は、数［３］より、数［５］として解析的に求めることができる。つまり、受音位置ｙｒｅｆにおけるＳＤＭに基づく無限直線二次音源からの音場が数［５］で求められる。つまり、ＳＤＭにおいて、数［５］のように、音圧ｐおよび伝達関数ｇをフーリエ変換することによって、駆動信号Ｄが求まることは既に、上記した参考文献１において公知である。
【００５５】
【数５】
2015231087.tif 000007
【００５６】
以上を前提にして、次に、図１実施例のように、駆動装置すなわちコンピュータ１６によって、ブライト領域およびダーク領域を形成するための手法について、説明する。
【００５７】
図２に示す無限直線二次音源を用いて、或る領域にのみ音源ｓ（ω）を再現するための各スピーカのフィルタをｆ（ｘ０，ω）とすると、無限直線二次音源の駆動信号ｄ（ｘ０，ω）は、数［６］として記述できる。
【００５８】
【数６】
2015231087.tif 000008
【００５９】
ただし、図１に示す音源ユニット１８から供給される音源ｓは任意の音源であってよく、モノラル音源を想定している。
【００６０】
この実施例では、数［６］において最適なフィルタf（x０,ω）を設計することにより、音源ｓ（ω）の図１実施例のような局所音響再生を実現する。
【００６１】
フィルタは音源に依存しないため，ｓ（ω）＝１とする。よって、フィルタｆ（ｘ０，ω）による受音位置ｙ＝ｙｒｅｆにおける音場の空間フーリエ変換は，数［３］により、数［７］となる。
【００６２】
【数７】
2015231087.tif 000009
【００６３】
したがって、数［５］と同様に、波数領域の空間フィルタＦ（ｋｘ，ω）は、数［８］となる。この数［８］で示す空間フィルタＦは、スピーカ毎のフィルタ関数ｆを空間フーリエ変換した結果を示す。
【００６４】
【数８】
2015231087.tif 000010
【００６５】
多点制御法において局所音響再生するためには、音の聞こえる場所の音圧を「１」とし、聞こえない場所を「０」としている。そこで、この実施例では、図３に示すように、ｘ＝０を中心にｙ＝ｙｂの距離に幅ｌｂの局所音響再生を実現するために、ｙ＝ｙｂにおける音圧ｐ（ｘ，ｙｙ）を聞こえるところを「１」、聞こえないところを「０」とする音圧分布の矩形窓で以下のようにモデル化し、コンピュータ１６は、数［９］に示す音圧モデルを求める。その意味で、数［９］に応じて音圧モデルを計算するコンピュータ１６は、モデル化手段として機能しているということができる。
【００６６】
【数９】
2015231087.tif 000011
【００６７】
この矩形窓は周波数に依存しないためωは省略できる。
【００６８】
そして、矩形窓ｐ（ｘ，ｙｂ）のｘ軸方向の空間フーリエ変換（参考文献２参照）は、数［１０］の通りである。
【００６９】
【数１０】
2015231087.tif 000012
【００７０】
さらに、局所音響再生を任意の位置ｘｂで実現するためには、シフトの定理（参考文献２参照）により、数［１０］は、数［１１］となる。
【００７１】
【数１１】
2015231087.tif 000013
2015231087.tif 000014
【００７２】
したがって、数［１１］を数［８］に代入すると、波数領域における無限直線二次音源による局所再生フィルタＦ（ｋｘ，ω）は、数［１２］として解析的に導出できる。
【００７３】
【数１２】
2015231087.tif 000015
【００７４】
この空間フィルタリングにより、無限直線二次音源を用いて任意の位置［ｘｂ，ｙｂ］Ｔに任意の幅ｌｂの領域の局所音響再生を実現できる。このように、駆動装置すなわちコンピュータ１６は、無限直線二次音源の方向（実施例では、直線スピーカアレイ１２のスピーカ１４の並び方向）における位置ｘｂと窓幅ｌｂを設定することによって、空間フィルタを導出する。したがって、コンピュータ１６は、位置規定手段および幅規定手段として機能しているということができる。さらに、駆動装置すなわちコンピュータ１６は、直線スピーカアレイ１２からの距離（ｙｂ）を用いて、空間フィルタを導出する。したがって、コンピュータ１６は、距離規定手段としても機能する。
【００７５】
上記数［７］の前提として音源ｓ＝１としたが、実際に図１の各スピーカ１４１、１４２、…、１４Ｎを駆動するためには、駆動装置すなわちコンピュータ１６は、数［１２］で導出した空間フィルタ関数ｆと音源ｓとによって数［１３］に従って駆動信号ｄを生成する。
［数１３］
ｄｉ＝ｓｉｆｉ
そして、駆動装置すなわちコンピュータ１６は、数［１３］で各スピーカ毎に生成した駆動信号ｄ１、ｄ２、…、ｄＮによって、スピーカ１４１、１４２、…、１４Ｎを駆動する。したがって、駆動装置すなわちコンピュータ１６は、駆動手段として機能している。
【００７６】
数［１２］までは計算のために無限直線二次音源を想定したが、数［１２］を用いた局所音響再生を図１に示すように、有限個Ｎ（実施例では６４個）のスピーカ１４１、１４２、…、１４Ｎを直線上に配置した直線スピーカアレイ１２を用いて実現する。そのためには数［１２］を離散化し、有限長で打ち切る必要がある。
【００７７】
ここで、直線スピーカアレイ１２のチャネル数（スピーカの個数）をＮとし、素子間隔をΔｘとすると、アレイ長はＬ＝（Ｎ-1）Δｘとなる．数［１２］の離散逆フーリエ変換（ＩＤＦＴ）は、数［１４］で得られる。
【００７８】
【数１４】
2015231087.tif 000016
【００７９】
ここで、ｘ＝ｎΔｘであり、ｋｘ＝２πｍ=／ＮΔｘであり、−Ｎ／２≦ｎ≦Ｎ／２−１である。
【００８０】
近接音響ホログラフィ（参考文献２）では、アレイの開口長の打ち切りによる誤差を低減するために、音圧収録面の開口長の両端に仮想的に「０」を挿入し、開口長の２倍以上の長さの面を用いて２Ｎ以上の点を用いてＤＦＴを計算する。この方法と同様に、図１実施例においても、図５に示すように、音圧分布の矩形窓の長さをアレイ長の２倍である２Ｌとし、２Ｎ点を用いたＩＤＦＴにより数［１４］を計算する。
【００８１】
そして、得られた２Ｎチャネルの駆動信号のうち、中心のＮチャネル分を実際の駆動信号として用いる。図１実施例でいえば、１２８チャネルの駆動信号のうちの中心の６４チャネル分を駆動信号ｄとする。このようにして、数［１２］を離散逆フーリエ変換することによって、実際の駆動信号ｄを計算する。
【００８２】
図１の実施例では、直線スピーカアレイ１２からの音を聴取者Ｐ２（の領域）にのみ再生し、他の聴取者では聞き取れないようにした。つまり、聴取者Ｐ２のための１つのブライト領域を形成した。これによって、１つの音源を聴取者Ｐ２にのみ再生できる。
【００８３】
これに対して、図４の実施例では、２つの音源ユニット１８１および１８２を設け、各音源ユニット１８１および１８２から２つの音源ｓ１（ω）および音源ｓ２（ω）を駆動装置すなわちコンピュータ１６に与える。他方、以下に説明する方法に従って、聴取者Ｐ２に対する第１のブライト領域と、聴取者Ｐ１に対する第２のブライト領域を設定する。そして、第１のブライト領域ではたとえば英語の音源ｓ１（ω）を再生し、第２のブライト領域ではたとえば日本語の音源ｓ２（ω）を再生する。したがって、この実施例は、先の図１の実施例と同様に、たとえばパーソナルオーディオシステム、美術館やオリンピックなどでの多言語同時ガイド、その他バーチャルリアリティ技術などに応用可能である。
【００８４】
また、特急列車や普通列車が共通に停車する駅のプラットホームにおいて、特急列車乗車予定者が並ぶ位置と普通列車のそれとが異なる場合が多々あるが、このような場合に、図４に示す実施例を応用して、たとえば特急列車に乗る人Ｐ１に対して第１ブライト領域を設定して特急列車向けのアナウンス（音源ｓ１）を流し、それと同時に普通列車に乗る人Ｐ２に対して第２のブライト領域を設定して普通列車向けのアナウンス（音源ｓ２）を流すことによって、それぞれの乗車予定者が誤乗車することなく正しい列車に乗車することが可能となる。
【００８５】
このように、Ｍ個の音源ｓｉ（ω）をＭ箇所別々のブライト領域で再生するマルチ局所音響再生を実現するためには、図２のように、数［１２］で導出された各局所音響再生位置でのフィルタｆｉと、対応する音源ｓｉの組み合わせｓｉ（ω）ｆｉ（ｘ０，ω）の重ね合わせにより、駆動信号を数［１５］で得る。つまり、図４の実施例では、駆動装置すなわちコンピュータ１６は、受音位置ｙ＝ｙｂにおいて、スピーカ１４の配列方向すなわちｘ方向において異なる２つの位置において、矩形窓を設定し、それぞれの矩形窓について数［１２］で導出した空間フィルタ関数と、音源ユニット１８１および１８２からの音源ｓ１およびｓ２とで駆動信号ｄｉを生成する。
【００８６】
【数１５】
2015231087.tif 000017
【００８７】
つまり、数［１５］に従って、図４に示す実施例において、駆動装置すなわちコンピュータ１６は、音源ｓｉと各スピーカ１４ｉの空間フィルタｆｉによって各スピーカ１４１、１４２、…、１４Ｎのための駆動信号ｄ１、ｄ２、…、ｄＮを生成して、各スピーカを駆動する。
【００８８】
ただし、先の説明では、複数の音源を複数のブライト領域で局所的に再生するようにしたが、１つの音源を複数のブライト領域で局所的に再生する場合でも、図４に示す実施例を利用可能であることは勿論である。つまり、図４に示す実施例は、複数のブライト領域を個別に設定する実施例である。
【００８９】
また、図４では、２つのブライト領域が直線スピーカアレイ１２から異なる聴取位置に設定されているように描かれているが、これは単に図解のためであり、実際には２つのブライト領域は同じ距離ｙ＝ｙｂ上に設定されるものであることに留意されたい。
【００９０】
ここで、コンピュータシミュレーションによって、この発明の有効性を検証する。局所音響再生精度を従来法の１つであるＥＤＭと比較する。
【００９１】
ＥＤＭではＮチャネルのスピーカを用いてＫチャネルの制御点を制御する。各スピーカ位置をｘｓｐとし、制御点位置をｘｃｏとすると、各スピーカと制御点間の空間相関行列は、
【００９２】
【数１６】
2015231087.tif 000018
【００９３】
として計算される。ここで、
【００９４】
【数１７】
2015231087.tif 000019
【００９５】
である。
【００９６】
ＥＤＭでは、各スピーカと局所音響再生領域の制御点との空間相関行列をｒｂ（ω）、各スピーカと無音領域の制御点との空間相関行列をｒｄ（ω）とし、各スピーカのフィルタｆ（ω）は、
【００９７】
【数１８】
2015231087.tif 000020
【００９８】
の最大固有値に対応する固有ベクトルとして計算される。ここで、αはチューニング係数であり、角周波数ω毎に設定する。
【００９９】
コンピュータシミュレーションによって、３次元自由音場を仮定した局所音響再生の精度を評価する。各スピーカは無指向性の点音源を仮定する。シミュレーション条件およびスピーカと局所音響再生制御ライン配置をそれぞれ表１および図６に示す。スピーカアレイはｘ軸上にｘ＝０を中心に配置される。ｘ＜０をブライト領域(ｘｂ， p＝１)とし、ｘ＞０をダーク領域(ｘｄ，p＝０)とする。ＥＤＭでは、制御点をｙ＝ｙｂ（＝２ｍ）上にスピーカと平行して設置し、ｘ＜０の領域の３２チャネルを局所音響再生の制御点ｘｂとし、残り半分の３２チャネルを無音制御領域の制御点ｘｄとして数［１７］を計算する。
【０１００】
チューニング係数は予備検討により、ｘｂ領域とｘｄ領域間でエネルギ差が最大となるように設定した。直線スピーカアレイを用いた場合、周波数に依存せず、α＝０.９９９９で最大となった。
【０１０１】
【表１】
2015231087.tif 000021
【０１０２】
たとえば図１実施例の方法では、数［１２］の窓関数の長さを２Ｌとし、２Ｎ＝１２８点のＩＤＦＴにより数［１４］を計算する。得られた１２８チャネルのｄ（ｘ０，ω）のうち、中心の６４チャネルを実際の駆動信号として用いる。
【０１０３】
局所音響再生の位置による精度を評価するため、位置ｘにおける音圧レベルＰＳＰＬ（ｘ）を、数［１９］のように定義する。
【０１０４】
【数１９】
2015231087.tif 000022
【０１０５】
ここで、ｐ（ｘ，ω）は位置ｘ、角周波数ωの音圧である。
【０１０６】
また、局所音響再生領域ｘｂと無音制御領域ｘｄ間の音圧レベルを評価するために、ブライト領域対ダーク領域比（Bright to dark ratio：ＢＤＲ)を、数［２０］のように定義する。
【０１０７】
【数２０】
2015231087.tif 000023
【０１０８】
それぞれの領域は図５に示すようにｙ＝ｙｂを中心とし、幅０．４ｍ、長さはＬ／２＝１.６ｍとして計算した。
【０１０９】
シミュレーション結果である音圧レベルの結果を図７に示し、ＢＤＲの結果を図８および図９に示す。これらの結果より、特に１０００Ｈｚ以上では実施例の方がＥＤＭよりも精度よく局所音響再生を実現できることがわかる。また、図８および図９の結果より、ＥＤＭでは周波数が高くなるにつれて局所音響再生の幅が狭くなるのに対して、実施例では周波数に依存せず、設定した幅の局所音響再生を実現できることが確認できる。
【０１１０】
以上より、実施例に従った空間フィルタリングにより、設定した局所音響再生位置を効果的に制御することができる。
【０１１１】
なお、上で挙げた寸法などの具体的数値は、いずれも単なる一例であり、製品の仕様などの必要に応じて適宜変更可能である。
【符号の説明】
【０１１２】
１０  …局所音響再生装置
１２  …直線スピーカアレイ
１４１、１４２、…、１４Ｎ  …スピーカ
１６  …駆動装置（コンピュータ）
【図１】
2015231087.tif 000024
【図２】
2015231087.tif 000025
【図３】
2015231087.tif 000026
【図４】
2015231087.tif 000027
【図５】
2015231087.tif 000028
【図６】
2015231087.tif 000029
【図７】
2015231087.tif 000030
【図８】
2015231087.tif 000031
【図９】
2015231087.tif 000032