(57)【要約】【課題】任意の指向特性を持つ音源によって作り出される音空間について、マイクロホンの近傍以外の任意の位置の音空間情報を得る。【解決手段】指向特性表現部１２１が、指向性を持つ音源を複数の点音源で表現し、第１行列算出部１２２が、各点音源によって生じる音空間情報の球面調和係数行列Ｂ(enc)を算出し、第２行列算出部１２３が、各点音源から任意の位置ｒtgtまでの伝搬作用素Ｂ(tgt)を算出し、移動作用素算出部１２４が、球面調和係数行列Ｂ(enc)と伝搬作用素Ｂ(tgt)の線形演算により移動作用素Ｔrec→tgtを算出し、合成部１４が、収録位置ｒrecでの音空間情報の展開係数Ｂ(rec)に移動作用素Ｔrec→tgtを適用し、任意の位置の音空間情報の展開係数Ｂ(trans)を求める。【選択図】図２
2019075616.tif 000002
【特許請求の範囲】
【請求項１】
指向性を持つ音源を複数の点音源で表現する音源表現手段と、
前記複数の点音源のそれぞれによって生じる音空間情報の球面調和係数行列を算出し、前記複数の点音源のそれぞれから任意の位置までの伝搬作用素を算出し、前記球面調和係数行列と前記伝搬作用素の線形演算により移動作用素を算出する算出手段と、
マイクロホンアレイで得られた音空間情報に前記移動作用素を適用して前記任意の位置での音空間情報を合成する合成手段と、
を有することを特徴とする音場収録装置。
【請求項２】
前記球面調和係数行列及び前記伝搬作用素は、行列
【数１５】
2019075616.tif 000017
で表され、上記行列の各要素は、各点音源から出力される音の球面調和係数であり、次式
【数１６】
2019075616.tif 000018
（ただし、ｋは波数、ｒs，θs，φsは点音源の位置、ｈn(1)は第一種球ハンケル関数、Ｙnmは次数ｎ，陪数ｍの球面調和関数であり、＊は複素共役を表す）
で求めることを特徴とする請求項１に記載の音場収録装置。
【請求項３】
コンピュータが実行する音場収録方法であって、
指向性を持つ音源を複数の点音源で表現するステップと、
前記複数の点音源のそれぞれによって生じる音空間情報の球面調和係数行列を算出するステップと、
前記複数の点音源のそれぞれから任意の位置までの伝搬作用素を算出するステップと、
前記球面調和係数行列と前記伝搬作用素の線形演算により移動作用素を算出するステップと、
マイクロホンアレイで得られた音空間情報に前記移動作用素を適用して前記任意の位置での音空間情報を合成するステップと、
を有することを特徴とする音場収録方法。
【請求項４】
前記球面調和係数行列及び前記伝搬作用素は、行列
【数１７】
2019075616.tif 000019
で表され、上記行列の各要素は、各点音源から出力される音の球面調和係数であり、次式
【数１８】
2019075616.tif 000020
（ただし、ｋは波数、ｒs，θs，φsは点音源の位置、ｈn(1)は第一種球ハンケル関数、Ｙnmは次数ｎ，陪数ｍの球面調和関数であり、＊は複素共役を表す）
で求めることを特徴とする請求項３に記載の音場収録方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音場再現のための音空間情報を取得する技術に関する。
【背景技術】
【０００２】
音空間情報の収録には、通常、複数のマイクロホンを有するマイクロホンアレイを用いる。特に、剛球上にマイクロホンをできるだけ等密度に配置した球状マイクロホンアレイは、形状がコンパクト、かつ、対称形であるため、比較的多くの場面で使用されている。ただし、マイクロホンの数に応じて収録される音空間の距離が規定されるため、マイクロホンアレイの近傍の音空間しか正確に収録することができない。この範囲を広げるためには、極めて多数のマイクロホンを稠密に配置する必要がある。
【０００３】
この問題を解決するために、発明者らは、非特許文献１において、収録する音源の位置を既知とすることで、マイクロホンの近傍以外のあらゆる位置での音空間情報を算出することができる移動作用素の算出方法を提案した。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】J. Trevino, K. Wakayama, S. Sakamoto, Y. Suzuki, H. Takada, and M. Okamoto, “Recording of extended sound fields using spherical microphone arrays and a-priori knowledge of the sound source positions,”The Journal of the Acoustical Society of America, 2016, Volume 139, Issue 4, http://dx.doi.org/10.1121/1.4950265
【非特許文献２】羽田陽一、古家賢一、島内末廣、“球調和関数展開に基づく多重極音源を用いた指向性合成”、日本音響学会誌、一般社団法人  日本音響学会、２０１３年、第６９巻、第１１号、pp. 577-588
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、非特許文献１の手法では、音源の特性を無指向性の点音源と限定して移動作用素を算出していた。
【０００６】
本発明は、上記に鑑みてなされたものであり、任意の指向特性を持つ音源によって作り出される音空間について、マイクロホンの近傍以外の任意の位置の音空間情報を得ることを目的とする。
【課題を解決するための手段】
【０００７】
第１の本発明に係る音場収録装置は、指向性を持つ音源を複数の点音源で表現する音源表現手段と、前記複数の点音源のそれぞれによって生じる音空間情報の球面調和係数行列を算出し、前記複数の点音源のそれぞれから任意の位置までの伝搬作用素を算出し、前記球面調和係数行列と前記伝搬作用素の線形演算により移動作用素を算出する算出手段と、マイクロホンアレイで得られた音空間情報に前記移動作用素を適用して前記任意の位置での音空間情報を合成する合成手段と、を有することを特徴とする。
【０００８】
第２の本発明に係る音場収録方法は、コンピュータが実行する音場収録方法であって、指向性を持つ音源を複数の点音源で表現するステップと、前記複数の点音源のそれぞれによって生じる音空間情報の球面調和係数行列を算出するステップと、前記複数の点音源のそれぞれから任意の位置までの伝搬作用素を算出するステップと、前記球面調和係数行列と前記伝搬作用素の線形演算により移動作用素を算出するステップと、マイクロホンアレイで得られた音空間情報に前記移動作用素を適用して前記任意の位置での音空間情報を合成するステップと、を有することを特徴とする。
【発明の効果】
【０００９】
本発明によれば、任意の指向特性を持つ音源によって作り出される音空間について、マイクロホンの近傍以外の任意の位置の音空間情報を得ることができる。
【図面の簡単な説明】
【００１０】
【図１】本実施形態の音場収録装置を用いて任意の位置での音空間情報を得る様子を示す図である。
【図２】本実施形態の音場収録装置の構成を示す機能ブロック図である。
【図３】球面調和関数と９つの点音源の位置関係を示す図である。
【発明を実施するための形態】
【００１１】
以下、本発明の実施の形態について図面を用いて説明する。
【００１２】
図１は、本実施形態の音場収録装置１を用いて任意の位置での音空間情報を得る様子を示す図である。音場収録装置１は、音源位置が既知の指向性を持つ音源によって作り出される音空間の、球状マイクロホンアレイ３の設置位置以外の任意の位置における音空間情報を合成する。
【００１３】
（音場収録装置の構成）
図２は、本実施形態の音場収録装置１の構成を示す機能ブロック図である。同図に示す音場収録装置１は、音源位置入力部１１、算出部１２、音空間情報入力部１３、合成部１４、及び音空間情報出力部１５を有する。音場収録装置１が備える各部は、演算処理装置、記憶装置等を備えたコンピュータにより構成して、各部の処理がプログラムによって実行されるものとしてもよい。このプログラムは音場収録装置１が備える記憶装置に記憶されており、磁気ディスク、光ディスク、半導体メモリ等の記録媒体に記録することも、ネットワークを通して提供することも可能である。以下、各部について説明する。
【００１４】
音源位置入力部１１は、指向性を持つ音源の音源位置など算出部１２が移動作用素を算出するために必要な情報を入力する。
【００１５】
算出部１２は、入力した音源位置に基づいて指向性を持つ音源を複数の点音源で表現し、指向性を持つ音源により任意の位置で作り出される音空間情報を算出するための移動作用素を求める。算出部１２の詳細については後述する。
【００１６】
音空間情報入力部１３は、球状マイクロホンアレイ３で収録されたマルチチャネル信号を球面調和関数領域で表現して得られた、収録位置での音空間情報を入力する。
【００１７】
合成部１４は、収録位置での音空間情報に算出部１２の求めた移動作用素を適用し、任意の位置の音空間情報を得る。
【００１８】
音空間情報出力部１５は、合成部１４で得られた任意の位置の音空間情報を出力する。
【００１９】
（算出部の処理）
続いて、算出部１２について説明する。算出部１２は、指向特性表現部１２１、第１行列算出部１２２、第２行列算出部１２３、及び移動作用素算出部１２４を有する。
【００２０】
指向特性表現部１２１は、音源位置入力部１１が入力した音源位置の周囲に複数の点音源を近接して配置し、任意の指向性を持つ音源を複数の点音源で表現する。
【００２１】
第１行列算出部１２２は、複数の点音源のそれぞれによって生じる音空間情報の球面調和係数行列を算出する。
【００２２】
第２行列算出部１２３は、複数の点音源のそれぞれから任意の位置までの伝搬作用素を算出する。
【００２３】
移動作用素算出部１２４は、第１行列算出部１２２の算出した球面調和係数行列と第２行列算出部１２３の算出した伝搬作用素の線形演算により移動作用素を算出する。
【００２４】
以下、算出部１２の処理に関して、複数の点音源による指向性合成、指向性合成の球面調和関数表現、球面調和係数の推定、及び移動作用素の導出について説明する。
【００２５】
（複数の点音源による指向性合成）
音源のもつ指向性を球面調和関数を用いて表すと次式（１）のようになる。
【００２６】
【数１】
2019075616.tif 000003
【００２７】
ここで、Ｄ（θ，φ）は、とある周波数における音源の指向性を表し、θとφはそれぞれ水平角と仰角である。Ｙnmは、字数ｎ、陪数ｍの球面調和関数を表す。式（１）から、音源の指向性は、球面調和関数の展開係数Ａnmで表されることがわかる。
【００２８】
このように算出された指向性は、近接した極性の異なる複数の点音源の重ね合わせで表現することが可能である（非特許文献２参照）。例えば、非特許文献２では、例として、２次の球面調和関数で表現できる任意の指向性を、以下のように配置した９つの点音源の重ね合わせにより表現できることを示している。
【００２９】
【数２】
2019075616.tif 000004
【００３０】
なお、ｄはそれぞれの点音源の間隔であり、表現すべき指向性を持つ音源の波長より小さくなる必要がある。また、ベクトルｒrecは、音空間を収録したマイクロホンの位置を表す。
【００３１】
元々の音源の指向性は式（１）のように展開係数［Ａ］nm＝Ａnmで表現できることから、球面調和関数領域でこの指向性を９つの点音源の重ね合わせとして表現する際には、各点音源に掛け合わせる重み係数ベクトルＨ＝［Ｈ1 Ｈ2・・・Ｈ9］Tを、行列Ｃ，Ｐを用いて次式（３）のように算出することとなる。
【００３２】
【数３】
2019075616.tif 000005
【００３３】
例えば、先に例を示した２次の球面調和関数で表される指向性は、次式（４），（５）として表現される。
【００３４】
【数４】
2019075616.tif 000006
【００３５】
【数５】
2019075616.tif 000007
【００３６】
ここで、ｋは波数である。なお、式（４），（５）は、非特許文献２の式と異なっているが、これは、単に両者で使用している球面調和関数の座標系の違い、および、正規化の違いであり、本質的には同じものである。図３に、使用した球面調和関数と９つの点音源の位置関係を示す。
【００３７】
（指向性合成の球面調和関数表現）
位置ｒs＝（ｒs，θs，φs）にある点音源により任意の位置ｒで生じる音圧は、球面調和関数を用いて以下の式（６）で表される。
【００３８】
【数６】
2019075616.tif 000008
【００３９】
ここで、ｊnとｈn(1)はそれぞれ、球ベッセル関数、第一種球ハンケル関数であり、＊は複素共役を表す。
【００４０】
式（６）から、位置ｒsにおける点音源により生じる音空間情報を空間フーリエ変換で表すと、その係数Ｂnm(s)は、次式（７）のように導かれる。
【００４１】
【数７】
2019075616.tif 000009
【００４２】
したがって、非特許文献２で示した複数の点音源による指向性合成手法は、近接するそれぞれの点音源によって生じる音空間情報の空間フーリエ変換係数Ｂnm(s)を用いて、球面調和関数領域で表現し直すことが可能である。すなわち、マイクロホンにより収録された展開係数Ｂnm(rec)は、式（７）に示す各点音源の係数の和として次式（８）で表されることになる。
【００４３】
【数８】
2019075616.tif 000010
【００４４】
（球面調和係数の推定）
式（８）を各点音源における係数を用いた行列により表現し直すと、各点音源によって生じる音空間情報の球面調和係数行列Ｂ(enc)は次式（９）のように表される。
【００４５】
【数９】
2019075616.tif 000011
【００４６】
この球面調和係数行列Ｂ(enc)の各要素は各点音源から出力される音の球面調和係数であり、式（２）により算出された距離を用いて式（７）から求める。したがって、任意の指向性を持つ音源によって作り出される音空間の、収録位置における球面調和係数Ｂ(rec)は、音源の指向性を表現するために各点音源に重み付けする重み係数ベクトルＨを用いて、式（３）と同じように次式（１０）のように表される。
【００４７】
【数１０】
2019075616.tif 000012
【００４８】
球面調和係数Ｂ(rec)はマイクロホンアレイを用いて収録されるため、算出すべき元の音源の展開係数Ａは、式（１０）から次式（１１）のように推定される。
【００４９】
【数１１】
2019075616.tif 000013
【００５０】
ここで、［・］+は一般逆行列を表し、多くの場合、Ｍｏｏｒｅ−Ｐｅｎｒｏｓｅの疑似逆行列として算出する。さらに、各音源に対する重み係数ベクトルＨは次式（１２）のように表現される。
【００５１】
【数１２】
2019075616.tif 000014
【００５２】
（移動作用素の導出）
式（１２）は、指向性を持つ音源を複数の点音源で表した際の各点音源に掛け合わせる重み係数ベクトルＨを算出する。重み係数ベクトルＨを用いて、指向性を持つ音源が作り出す音空間を任意の位置で合成する。
【００５３】
式（９）で表される各点音源によって生じる音空間情報の球面調和係数行列Ｂ(enc)を、各点音源から該当する点までの音の伝搬行列として解釈する。すなわち、式（９）を、各点音源から収録位置ｒrecまでの音空間情報の伝搬特性として考える。このように考えると、各点音源から任意の位置ｒtgtまでの伝搬作用素Ｂ(tgt)は、式（２）におけるｒrecをｒtgtに入れ替えた上で、式（９），（７）を適用して求めることができる。
【００５４】
この伝搬作用素Ｂ(tgt)と式（１２）から、収録位置ｒrecからターゲットとなる任意の位置ｒtgtへ移動した際の展開係数を次式（１３）で求めることができる。
【００５５】
【数１３】
2019075616.tif 000015
【００５６】
収録位置ｒrecから任意の位置ｒtgtへの展開係数の移動作用素Ｔrec→tgtは次式（１４）で表される。
【００５７】
【数１４】
2019075616.tif 000016
【００５８】
本実施形態では、指向特性表現部１２１が、指向性を持つ音源を複数の点音源で表現し、第１行列算出部１２２が、各点音源によって生じる音空間情報の球面調和係数行列Ｂ(enc)を算出し、第２行列算出部１２３が、各点音源から任意の位置ｒtgtまでの伝搬作用素Ｂ(tgt)を算出し、移動作用素算出部１２４が、球面調和係数行列Ｂ(enc)と伝搬作用素Ｂ(tgt)の線形演算により移動作用素Ｔrec→tgtを算出する。そして、合成部１４が、収録位置ｒrecでの音空間情報の展開係数Ｂ(rec)に移動作用素Ｔrec→tgtを適用し、任意の位置の音空間情報の展開係数Ｂ(trans)を求めることにより、ある収録位置で球状マイクロホンアレイ３で収録しさえすれば、任意の指向性を持つ音源によって作り出される音空間について、収録位置以外の任意の位置の音空間情報を得ることができる。
【符号の説明】
【００５９】
１…音場収録装置
１１…音源位置入力部
１２…算出部
１２１…指向特性表現部
１２２…第１行列算出部
１２３…第２行列算出部
１２４…移動作用素算出部
１３…音空間情報入力部
１４…合成部
１５…音空間情報出力部
３…球状マイクロホンアレイ
【図１】
2019075616.tif 000021
【図２】
2019075616.tif 000022
【図３】
2019075616.tif 000023