(57)【要約】【課題】計算時間が少なく、且つ指向特性を柔軟に制御可能な音像定位装置を提供する。【解決手段】複数のスピーカＳＰ１〜ＳＰＱを直線上に配置させたスピーカアレー４０から放射させた音響信号を、反射板５０に反射させて音像を定位させる音像定位装置であって、所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出する展開係数算出部１０と、展開係数をスピーカＳＰ１〜ＳＰＱのそれぞれに対応させたフィルタ係数に変換するフィルタ係数生成部２０と、音声信号に、フィルタ係数を畳み込んでスピーカＳＰ１〜ＳＰＱをそれぞれ駆動するスピーカ駆動信号を生成するスピーカ駆動部３０とを備える。【選択図】図１
2020170961.tif 000002
【特許請求の範囲】
【請求項１】
複数のスピーカを直線上に配置させたスピーカアレーから放射させた音響信号を、反射板に反射させて音像を定位させる音像定位装置であって、
所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出する展開係数算出部と、
前記展開係数から、前記スピーカのそれぞれに対応させたフィルタ係数を生成するフィルタ係数生成部と、
音声信号に、前記フィルタ係数を畳み込んで前記スピーカをそれぞれ駆動するスピーカ駆動信号を生成するスピーカ駆動部と
を備えることを特徴とする音像定位装置。
【請求項２】
前記フィルタ係数生成部は、
前記展開係数のそれぞれに、球面調和関数を前記スピーカで再現する場合の該スピーカのそれぞれ対応する重みを乗じて前記フィルタ係数を生成する
ことを特徴とする請求項１に記載の音像定位装置。
【請求項３】
複数のスピーカを直線上に配置させたスピーカアレーから放射させた音響信号を、反射板に反射させて音像を定位させる音像定位装置が実行する音像定位方法であって、
所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出する展開係数算出ステップと、
前記展開係数から、前記スピーカのそれぞれに対応させたフィルタ係数を生成するフィルタ係数生成ステップと、
音声信号に、前記フィルタ係数を畳み込んで前記スピーカをそれぞれ駆動するスピーカ駆動信号を生成するスピーカ駆動ステップと
を行うことを特徴とする音像定位方法。
【請求項４】
請求項１又は２に記載の音像定位装置としてコンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音像定位装置、音像定位方法及びプログラムに関し、特に、スピーカ本体ではなく任意の位置に仮想的な音源を生成する演出効果を持つ音響再生技術に関する。
【背景技術】
【０００２】
近年、パブリックビューイングや家庭において、スピーカを複数配置した再生方式が広まっている。また、3D映像・ワイド映像などの映像技術が広まるとともに、音響に関してもスピーカ本体ではなく任意の位置に仮想的な音源を生成することで、より高い臨場感の感じられる再生を実現する取り組みが行われている。
【０００３】
音の反射を用いて仮想的なスピーカを創り出す音響再生技術として、指向性スピーカからの放射音と反射板からの反射音の総和が任意の点で最大になるように指向性を制御し、局所再生を実現する方法が例えば特許文献１に記載されている。また、正多面体スピーカの指向性再生により、天井に音を反射させることで上方への音像定位を実現する方法が例えば非特許文献１に開示されている。
【０００４】
非特許文献１では、天井からの反射音とスピーカからの直接音の差が5dBより大きい場合に、上方への音像の定位が可能であることが報告されている。複数人に、上方に音像が在ることを知覚させるためには、再生音の指向特性を任意の形状に制御する必要がある。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１２−８１５６号公報
【非特許文献】
【０００６】
【非特許文献１】H. Sakamoto, Y. Haneda,”Sound Localization of Beamforming-Controlled Reflected Sound from Ceiling in Presence of Direct Sound,” in 144th Audio Engineering Society Convension paper 9949, 2018, May.
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、従来の指向性制御により指向特性を柔軟に変化させるためには、多くの制御点を用いる必要があり、多くの計算時間を要するという課題がある。
【０００８】
本発明は、この課題に鑑みてなされたものであり、計算時間が少なく、且つ指向特性を柔軟に制御可能な音像定位装置、音像定位方法及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明の一態様に係る音像定位装置は、複数のスピーカを直線上に配置させたスピーカアレーから放射させた音響信号を、反射板に反射させて音像を定位させる音像定位装置であって、所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出する展開係数算出部と、前記展開係数を前記スピーカのそれぞれに対応させたフィルタ係数に変換するフィルタ係数生成部と、音声信号に、前記フィルタ係数を畳み込んで前記スピーカをそれぞれ駆動するスピーカ駆動信号を生成するスピーカ駆動部とを備えることを要旨とする。
【００１０】
また、本発明の一態様に係る音像定位方法は、上記の音像定位装置が実行する音像定位方法であって、複数のスピーカを直線上に配置させたスピーカアレーから放射させた音響信号を、反射板に反射させて音像を定位させる音像定位装置が実行する音像定位方法であって、所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出する展開係数算出ステップと、前記展開係数から、前記スピーカのそれぞれに対応させたフィルタ係数を生成するフィルタ係数生成ステップと、音声信号に、前記フィルタ係数を畳み込んで前記スピーカをそれぞれ駆動するスピーカ駆動信号を生成するスピーカ駆動ステップとを行うことを要旨とする。
【００１１】
また、本発明の一態様に係るプログラムは、上記の音像定位装置としてコンピュータを機能させるためのプログラムであることを要旨とする。
【発明の効果】
【００１２】
本発明によれば、計算時間が少なく、且つ指向特性を柔軟に制御可能な音像定位装置、音像定位方法及びプログラムを提供することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る音像定位装置の構成例を示すブロック図である。
【図２】エンドファイア方向への音響ビームを模式的に示す図である。
【図３】極座標系を示す図である。
【図４】次数n=3までの球面調和関数の例を模式的に示す図である。
【図５】図１に示す音像定位装置が実行する処理手順を示すフローチャートである。
【図６】図１に示す音像定位装置で提供される音像定位の様子を模式的に示す図である。
【図７】最小二乗法による指向性制御のフィルタを設計する場合の観測系を模式的に示す図である。
【発明を実施するための形態】
【００１４】
以下、本発明の実施形態について図面を用いて説明する。複数の図面中同一のものに
は同じ参照符号を付し、説明は繰り返さない。
【００１５】
図１は、本発明の実施形態に係る音像定位装置の構成例を示すブロック図である。
【００１６】
図１に示す音像定位装置１００は、仮想的なスピーカを、従来の指向性制御のように制御点を配置するのではなく任意の窓幅を有する窓関数を球面調和関数展開することで、解析的に展開係数を導出し、直線スピーカアレーによる多重極音源によって再現する。この方法によれば、計算時間が少なく、且つビーム幅を柔軟に制御可能なエンドファイア方向への音響ビームを生成して仮想的なスピーカを形成し、複数の受聴者に音像を提示することができる。エンドファイア方向とは、一次元アレーの軸に沿う方向である。
【００１７】
図２は、エンドファイア方向への音響ビームを模式的に示す図である。図２（ａ）及び（ｂ）は、音響ビームの幅の違いを模式的に示す。図２（ａ）は、音響ビームの幅が狭い。図２（ｂ）は、音響ビームの幅が広い。本実施形態に係る音像定位装置１００は、従来のように多くの制御点を設けることなく、図２に示す音響ビームの幅の制御を実現する。
【００１８】
図１に示すように、本実施形態に係る音像定位装置１００は、展開係数算出部１０、フィルタ係数生成部２０、スピーカ駆動部３０、スピーカアレー４０、及び反射板５０を備える。スピーカアレー４０と反射板５０を除く音像定位装置１００は、例えば、ＲＯＭ、ＲＡＭ、ＣＰＵ等からなるコンピュータで実現することができる。その場合、音像定位装置１００が有すべき機能の処理内容はプログラムによって記述される。なお、スピーカアレー４０は、複数のスピーカＳＰ１〜ＳＰＱが直線上に配置されて形成される例を示す。
【００１９】
展開係数算出部１０は、所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出する。所望の指向性は、外部からビーム幅θｗ（0＜θｗ≦π）で与えられる。
【００２０】
窓関数は、コサイン窓（式（１））を例に説明する。なお、他の窓関数としては矩形窓がある。
【００２１】
【数１】
2020170961.tif 000003
【００２２】
（球面調和関数）
ここで図３に示す極座標系を考える。この場合、球面上の任意の点で観測される音圧S(r,θ,φ,ω)は次式で表せる。
【００２３】
【数２】
2020170961.tif 000004
【００２４】
ここで、Ymn(θ,φ)は球面調和関数、Amn(ω)はその展開係数であり次式で表せる。
【００２５】
【数３】
2020170961.tif 000005
【００２６】
ここで、Pmn(・)はルジャンドル陪関数であり、式（４）は球面調和関数展開と称される。
【００２７】
図４は、次数n=3までの球面調和関数の例を模式的に示す図である。位数mが0以上のものは実部を示し、位数mが0より小さいものは虚部を示す。
【００２８】
式（１）でモデル化した所望特性d（θ）を式（２）のS(r,θ,φ,ω)に代入し、球面調和関数の位数mを0とし、球面調和展開を行うことで多重極音源に対応する展開係数A0nが得られる。
【００２９】
【数４】
2020170961.tif 000006
【００３０】
次数n=2までの展開係数を以下に示す。
【００３１】
【数５】
2020170961.tif 000007
【００３２】
n=2以降の次数に対しても同様に解析的に展開係数を導出することが可能である。
【００３３】
フィルタ係数生成部２０は、展開係数Amnから、次式によりスピーカアレー４０を構成するそれぞれのスピーカに対応させたフィルタ係数を生成する（ステップＳ２（図２））。
【００３４】
（多重極音源による指向性制御技術）
球面調和関数により所望の指向特性を展開し、得られた展開係数A0nを多重極音源に適用することで指向性を形成する手法が知られている（例えば参考文献：羽田陽一他２名「球面調和関数展開に基づく多重極音源を用いた指向性合成」、日本音響学会誌、69.11,2013,577-588）。
【００３５】
多重極音源とは、原点に限りなく近い位置に同じ振幅を持つ点音源を逆位相で分布している音源である。例えば、ｚ軸方向に微少間隔dで点音源を配置した場合、多重極音源の音圧分布M0n(r,θ,φ,ω)は次式で表せる。
【００３６】
【数６】
2020170961.tif 000008
【００３７】
近似は1≪krのときに成り立ちz=cosθである。Qは点音源の強度を表す。kは波数（k=ω/c）である。また、多重極音源は球面調和関数と非常に類似した指向特性を持ち、z軸方向に配置したスピーカアレー４０では、位数mが0の球面調和関数と同様の指向特性を再現することができる。
【００３８】
つまり、多重極音源への適用は次式で表せる。
【００３９】
【数７】
2020170961.tif 000009
【００４０】
フィルタ係数生成部２０は、展開係数Amnのそれぞれに、球面調和関数をスピーカＳＰ１〜ＳＰＱで再現する場合の該スピーカのそれぞれ対応する重みD0n(ω)を乗じてフィルタ係数w(ω)を生成する（式（１１））。
【００４１】
【数８】
2020170961.tif 000010
【００４２】
重みD0n(ω)は、例えば次数n=2までの球面調和関数に対応するスピーカの数を5個とした場合、次式で表せる。
【００４３】
【数９】
2020170961.tif 000011
【００４４】
ここでdは、スピーカＳＰ１〜ＳＰＱのそれぞれの間隔（上記の微少間隔）である。また、kは波数（k=ω/c）、cは光速である。
【００４５】
スピーカ駆動部は、外部から入力される音声信号に、フィルタ係数w(ω)を畳み込んでスピーカＳＰ１〜ＳＰＱをそれぞれ駆動するスピーカ駆動信号を生成する。式（１２）から明らかなように、次数n=0のスピーカ駆動信号はスピーカＳＰ３のみにA0n（1/4π）0.5が入力される。次数n=1のスピーカ駆動信号はスピーカＳＰ２とＳＰ４に入力される。次数n=2のスピーカ駆動信号はスピーカＳＰ２とＳＰ３とＳＰ４に入力される。
【００４６】
このようなスピーカ駆動信号をスピーカアレー４０に入力することで、所望の指向特性に対応する音響信号を再生することができる。
【００４７】
以上説明したように本実施形態に係る音像定位装置１００は、複数のスピーカを直線上に配置させたスピーカアレー４０から放射させた音響信号を、反射板５０に反射させて音像を定位させる音像定位装置であって、所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出する展開係数算出部１０と、展開係数Amnから、スピーカＳＰ１〜ＳＰＱのそれぞれに対応させたフィルタ係数w(ω)を生成するフィルタ係数生成部２０と、音声信号に、フィルタ係数w(ω)を畳み込んでスピーカＳＰ１〜ＳＰＱをそれぞれ駆動するスピーカ駆動信号を生成するスピーカ駆動部３０とを備える。
【００４８】
これにより計算時間が少なく、且つ指向特性を柔軟に制御可能な音像定位装置１００を提供することができる。
【００４９】
（音像定位方法）
次に、音像定位装置１００が実行する音像定位方法について説明する。
【００５０】
図５は、音像定位装置１００が実行する処理手順を示すフローチャートである。
【００５１】
先ず始めに、音像定位装置１００に所望の指向性を表すビーム幅を設定する（ステップＳ１）。ビーム幅θw（式（１））は、展開係数算出部１０に外部から入力される(ステップＳ１)。
【００５２】
次に、展開係数算出部１０は、所望の指向性d(θ)を表す窓関数を球面調和関数展開することで解析的に展開係数Amnを算出する（ステップＳ２）。
【００５３】
次に、フィルタ係数生成部２０は、展開係数Amnをスピーカアレー４０を構成するそれぞれのスピーカＳＰ１〜ＳＰＱに対応させたフィルタ係数w(ω)を生成する（ステップＳ３）。フィルタ係数生成部２０は、展開係数Amnのそれぞれに、球面調和関数をスピーカＳＰ１〜ＳＰＱで再現する場合の該スピーカＳＰ１〜ＳＰＱのそれぞれ対応する重みD0n(ω)を乗じてフィルタ係数w(ω)を生成する（式（１１））。
【００５４】
スピーカ駆動部３０は、外部から入力される音声信号に、フィルタ係数w(ω)を畳み込んでスピーカＳＰ１〜ＳＰＱをそれぞれ駆動するスピーカ駆動信号を生成する（ステップＳ４）。
【００５５】
以上説明したように本実施形態に係る音像定位方法は、複数のスピーカＳＰ１〜ＳＰＱを直線上に配置させたスピーカアレー４０から放射させた音響信号を、反射板５０に反射させて音像を定位させる音像定位装置１００が実行する音像定位方法であって、所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数Amnを算出する展開係数算出ステップＳ２と、展開係数AmnをスピーカＳＰ１〜ＳＰＱのそれぞれに対応させたフィルタ係数w(ω)を生成するフィルタ係数生成ステップＳ３と、音声信号に、フィルタ係数w(ω)を畳み込んでスピーカＳＰ１〜ＳＰＱをそれぞれ駆動するスピーカ駆動信号を生成するスピーカ駆動ステップＳ４とを行う。これにより、計算時間が少なく、且つ指向特性を柔軟に制御可能な音像定位方法を提供することができる。
【００５６】
図６は、本実施形態に係る音像定位装置１００及び音像定位方法で提供される音像定位の様子を模式的に示す図である。図６に示すように、音像定位装置１００は、反射板５０（例えば天井）に音響信号を放射し、上方への音像定位（仮想スピーカＫＳＰ）を実現する。
【００５７】
符号１０３は直接音、符号１０４は反射音、及び符号１０５は受聴点を示す。音像定位装置１００によれば、多くの制御点を用いることなく受聴点１０５に位置する受聴者に、上方への音像定位を知覚させることができる。
【００５８】
（比較例）
図７は、最小二乗法による指向性制御のフィルタを設計する場合の観測系を模式的に示す図である。図７に示すスピーカアレー４０の周囲を円環状に囲むのは制御点1〜Mである。
【００５９】
最小二乗法による指向性制御は、所望の指向特性と、制御点で観測される指向特性との誤差に二乗和を最小化するフィルタ係数を求める。よって、計算量が多くなる。最小二乗法による指向性制御は、周知であるので式を示した説明は省略する。
【００６０】
また、非特許文献１に基づく方法は、正多面体スピーカの指向性再生により、天井で音を反射させることで上方への音像定位を実現している。この方法は、正規化マッチドフィルタを用いて指向性を形成している。
【００６１】
正規化マッチドフィルタは、スピーカから放射された音響信号が観測点で観測されるとき、観測される音響信号とスピーカが発する音響信号とが一致するフィルタを与えることにより得られる。よって、目標とする観測点までの伝達関数を全てのスピーカについて求める必要があり計算量が多くなる。
【００６２】
これらの比較例に対して本実施形態に係る音像定位方法は、所望の指向性を表す窓関数を球面調和関数展開することで解析的に展開係数を算出し、展開係数から、スピーカのそれぞれに対応させたフィルタ係数を生成するので計算量を少なくすることができる。つまり、計算時間が少なく、且つ指向特性を柔軟に制御可能な音像定位方法を提供することができる。
【００６３】
本実施形態に係る音像定位装置１００の特徴的な機能部は、ＲＯＭ、ＲＡＭ、ＣＰＵ等からなるコンピュータで実現することができる。その場合、各機能部が有すべき機能の処理内容はプログラムによって記述される。そのようなプログラムは、ＣＤ−ＲＯＭ等の記録媒体やインターネット等の伝送媒体を介して配信することができる。
【００６４】
本発明はここでは記載していない様々な実施形態等を含むことは勿論である。したがって、本発明の技術的範囲は上記の説明から妥当な特許請求の範囲に係る発明特定事項によってのみ定められるものである。
【符号の説明】
【００６５】
１０：展開係数算出部
２０：フィルタ係数生成部
３０：スピーカ駆動部
４０：スピーカアレー
５０：反射板（天井）
１００：音像定位装置
１０３：直接音
１０４：反射音
１０５：受聴点
【図１】
2020170961.tif 000012
【図２】
2020170961.tif 000013
【図３】
2020170961.tif 000014
【図４】
2020170961.tif 000015
【図５】
2020170961.tif 000016
【図６】
2020170961.tif 000017
【図７】
2020170961.tif 000018