JP5741866B2 - 音場収音再生装置、方法及びプログラム - Google Patents

Description

この発明は、ある音場に設置されたマイクロホンで音信号を収音し、その音信号を用いてスピーカでその音場を再現する波面合成法（Wave Field Synthesis）の技術、アンビソニックス(Ambisonics)の技術に関する。

ある音場に設置されたマイクアレーで信号を収音し、その信号を用いてスピーカアレーでその音場を再現する波面合成法（Wave Field Synthesis）の技術、アンビソニックス(Ambisonics)の技術として、例えば非特許文献１に記載された技術が知られている。

小山翔一，外３名，「音場収音・再現のための時空間周波数領域信号変換法」，日本音響学会講演論文集，２０１１年９月，Ｐ．６３５−６３６

非特許文献１に記載された技術では、上下方向も含めて再現するためには、平面状アレーを用いる必要があった。また、上下方向に用いるアレーが少数であった場合、再現誤差が大きくなってしまい、上下感の再現をすることができない可能性があった。

この発明の目的は、従来よりも高精度に音場を再現することができる音場収音再生装置、方法及びプログラムを提供することである。

上記の課題を解決するために、この発明の一態様による音場収音再生装置は、半径R_bの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする半径R_mの２個以上の円のそれぞれに少なくとも２個のマイクロホンが配置されているとし、R_m>R_bとし、バッフルの軸方向をz軸方向とし、バッフルの周方向とφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、nをφ方向の次数とし、k_z,lをz軸方向の波数とし、lをそのインデックスとし、J_n(・)をn次の第一種ベッセル関数とし、H_n ⁽¹⁾(・)をn次の第一種ハンケル関数とし、J_n’(・)をn次の第一種ベッセル関数J_n(・)の微分とし、H_n ⁽¹⁾’(・)をn次の第一種ハンケル関数H_n ⁽¹⁾(・)の微分とし、A(ω)を所定の複素数とし、w_nlをn,lに基づいて定まる重みとし、少なくとも４個のスピーカが第一の空間と異なる第二の空間の半径R_sの仮想円筒の周面に配置されているとして、マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~_nl(ω)に対して次式により定義されるフィルタF~_nl(ω)を適用してフィルタ処理後信号D~_nl(ω)を生成する変換フィルタ部と、

空間の逆フーリエ変換により、フィルタ処理後信号D~_nl(ω)を周波数領域信号に変換する空間周波数逆変換部と、周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号をスピーカに出力する周波数逆変換部と、を備える。

この発明の他の一態様による音場収音再生装置は、半径R_bの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする半径R_mの２個以上の円のそれぞれに少なくとも２個のマイクロホンが配置されているとし、R_m>R_bとし、バッフルの円筒の軸方向をz軸方向とし、バッフルの円筒の周方向とφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、nをφ方向の次数とし、k_z,lをz軸方向の波数とし、lをそのインデックスとし、J_n(・)をn次の第一種ベッセル関数とし、H_n ⁽¹⁾(・)をn次の第一種ハンケル関数とし、J_n’(・)をn次の第一種ベッセル関数J_n(・)の微分とし、H_n ⁽¹⁾’(・)をn次の第一種ハンケル関数H_n ⁽¹⁾(・)の微分とし、A(ω)を所定の複素数とし、w_nlをn,lに基づいて定まる重みとし、少なくとも４個のスピーカが第一の空間と異なる第二の空間の半径R_sの仮想円筒の周面に配置されているとして、マイクアレーで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換部と、空間のフーリエ変換により、周波数領域信号を時空間周波数領域信号P~_nl(ω)に変換する空間周波数変換部と、時空間周波数領域信号P~_nl(ω)に対して次式により定義されるフィルタF~_nl(ω)を適用してフィルタ処理後信号D~_nl(ω)を生成する変換フィルタ部と、を備える。

マイクロホンアレー及びスピーカアレーをそれぞれ構成するマイクロホン及びスピーカの数が少ない場合であっても、マイクロホンアレー及びスピーカアレーを円筒状とし、左右方向に密な素子配置とし、上下方向は粗な素子配置として、上下感を再現することができる。したがって、従来よりも高精度に音場を再現することができる

音場収音再生装置の例を示す機能ブロック図。 マイクロホン及びスピーカの配置の例を説明するための図。 マイクロホン及びスピーカの配置の例を説明するための図。 音場収音再生方法の例を示す流れ図。

以下、図面を参照してこの発明の実施形態を説明する。以下の説明において、テキスト中で使用する記号「~」、「⁻」等は、本来直前の文字の真上に記載されるべきものであるが、テキスト記法の制限により、当該文字の直後に記載する。式中においてはこれらの記号は本来の位置に記述している。また、ベクトルや行列の各要素単位で行われる処理は、特に断りが無い限り、そのベクトルやその行列の全ての要素に対して適用されるものとする。

＜マイクロホンアレー及びスピーカアレーの配置＞
音場収音再生装置及び方法は、図１に示すように、第一の空間の半径R_bの円筒形状の剛体のバッフルの軸からR_m離れた位置に配置されているN_z×N_φ個のマイクロホンで構成されるマイクロホンアレーと、第一の空間と異なる第二の空間の半径R_sの仮想円筒の周面に配置されているN_z×N_φのスピーカで構成されるスピーカアレーとを用いて、第一の空間の音源Ｓｏで発生した音によって形成された第一の空間の音場を第二の空間で再現する。図１では、第二の空間で再現された音源Ｓｏを音源Ｓｏ’と表現している。円筒の軸方向をz軸方向とする。第一の空間及び第二の空間は、互い異なる空間である。

第一の空間に配置されたマイクロホンの数と第二の空間に配置されたスピーカの数は異なっていてもよい。マイクロホンの数が、第二の空間に配置されたスピーカの数よりも多い場合には、再生信号を間引けばよい。一方、マイクロホンの数が、第二の空間に配置されたスピーカの数よりも少ない場合には、再生信号をチャネル間で平均を取るなどして補間を行えばよい。補間を行う方法は、例えば、線形補間やsinc補間などを適用することができる。

図３に示すように、マイクロホンは半径R_bの円筒形状のバッフルＢの軸からR_m離れた位置に配置されることになる。言い替えると、R_m>R_bとして、マイクロホンはバッフルＢの周面の表面から(R_m-R_b)離れた位置に配置される。例えば、バッフルＢの周面から垂直に突き出す細い棒状の部材により支持することで、マイクロホンは配置される。

換言すれば、円筒形状の剛体のバッフルＢの軸を中心とし円筒形状の剛体のバッフルＢの周方向を円周方向とするN_z個の円のそれぞれにおいて、N_φ個のマイクロホンが等間隔に配置される。N_z,N_φは予め定められた２以上の整数である。すなわち、バッフルＢの周方向を円周方向とする２個の円のそれぞれにおいて２個のマイクロホンが配置されることにより、マイクロホンはバッフルＢの軸からR_m離れた位置に少なくとも４個配置される。

マイクロホンは、どのような間隔で配置されてもよい。すなわち、隣接するマイクロホンとの間隔であるz_c,φ_cのそれぞれは、任意の値を取ることができる。ただし、マイクロホンを等間隔に配置する、すなわち隣接するマイクロホンとの間隔であるz_c,φ_cのそれぞれを同じ値とすることで、高精度に音場を再現することができる。

マイクロホンは、円筒形状の剛体のバッフルの周面の外側方向に向けて配置される。

スピーカもマイクロホンと同様に配置される。すなわち、図３に示すように、仮想円筒の周面のN_z個の円のそれぞれにおいて、N_φ個のスピーカが等間隔に配置される。N_z,N_φは予め定められた２以上の整数である。スピーカは、仮想円筒の周面に少なくとも４個配置される。

仮想円筒の周面のN_z個の円は、例えばz_cを所定の距離として、z_c間隔で位置している。また、同一の円に配置されたN_φ個のスピーカは、φ_cを所定の角度として、φ_c度の間隔で位置している。

スピーカは、ほぼ等間隔に配置されていれば、厳密に等間隔に配置されている必要はない。すなわち、隣接するスピーカとの間隔であるz_c,φ_cのそれぞれは、厳密に同じ値である必要はなく、ほぼ同じ値であればよい。また、スピーカは、どのような間隔で配置されてもよい。すなわち、隣接するスピーカとの間隔であるφ_cは、任意の値を取ることができる。ただし、スピーカをほぼ等間隔に配置する、すなわち隣接するスピーカとの間隔であるφ_cをほぼ同じ値とすることで、高精度に音場を再現することができる。

マイクロホンが配置されるバッフルの軸からの距離R_mは、例えば2cm程度とする。バッフルの軸からの距離R_mは、値が大きいほど鋭い指向性を形成することができるが、より多くのマイクロホンが必要となる。また、R_m,R_sは、その値が大きいほど広い領域を再現することができるが、より多くのマイクロホン及びスピーカが必要となる。R_m,R_sは、収音する信号の周波数を考慮して実験的に設定することが望ましい。

また、スピーカが配置される円筒の半径R_sは、例えば1.5m程度とする。なお、R_m≦R_sでもR_s≦R_mでもよいが、R_m≦R_sのときに音場再生の精度がよくなる。

スピーカは、音響的に透明な状態で第二の空間の空中に配置されてもよいし、音響的に透明でない状態で第二の空間に配置されてもよい。音響的に透明な状態とは、スピーカが配置されていない第二の空間の伝達特性と同じ伝達特性を保った状態ということである。例えば、スピーカは、糸で吊るされるか、細い棒で固定されることにより、第二の空間の空中に配置される。また、スピーカは、マイクロホンと同様に、円筒形状の剛体のバッフルの周面から所定の距離だけ離れた位置に配置されてもよい。この場合、スピーカは、剛体のバッフルの周面の内側方向に向けて配置される。

第一の空間のマイクロホンＭｉ−ｊの位置を(R_m,φ_m,i,z_m,j)[i=1,2,…,N_φ,j=1,2,…,N_z]と表現する。第二の空間のスピーカＳｉ−ｊの位置を(R_s,φ_s,i,z_s,j)[i=1,2,…,N_φ,j=1,2,…,N_z]と表現する。

＜音場収音再生装置＞
音場収音再生装置は、図１に示すように周波数変換部１、空間周波数変換部２、変換フィルタ部３、空間周波数逆変換部４、周波数逆変換部５及び窓関数部６を例えば含み、図４に例示された各ステップの処理を行う。

第一の空間に配置されたマイクロホンＭ１−１，Ｍ２−１，…，ＭＮ_φ−Ｎ_ｚは、第一の空間の音源Ｓで発せられた音を収音して時間領域の信号を生成する。生成された信号は、周波数変換部１に送られる。(R_m,φ_m,i,z_m,j)のマイクロホンＭｉ−ｊで収音された時間領域の時刻ｔの信号をp_ij(t)と表記する。

＜周波数変換部１＞
周波数変換部１は、マイクロホンＭ１−１，Ｍ２−１，…，ＭＮ_φ−Ｎ_ｚで収音された信号p_ij(t)をフーリエ変換により周波数領域信号P_ij(ω)に変換する（ステップＳ１）。生成された周波数領域信号P_ij(ω)は、空間周波数変換部２に送られる。ωは周波数である。例えば、短時間離散フーリエ変換により周波数領域信号P_ij(ω)が生成される。もちろん、他の既存の方法により周波数領域信号P_ij(ω)を生成してもよい。また、オーバーラップアド等の方法を用いて周波数領域信号P_ij(ω)を生成してもよい。入力信号が長い場合や、リアルタイム処理のように連続して信号が入力される場合には、例えば１０ｍｓごとといったフレームごとに処理を行う。周波数領域信号P_ij(ω)は、例えば以下のように定義される。関数expの引数の中のjは虚数単位である。

＜空間周波数変換部２＞
空間周波数変換部２は、空間のフーリエ変換により周波数領域信号P_ij(ω)を時空間周波数領域信号P~_nl(ω)に変換する（ステップＳ２）。時空間周波数領域信号P~_nl(ω)は、各ωごとに計算される。変換された時空間周波数領域信号P~_nl(ω)は、変換フィルタ部３に送られる。空間周波数変換部２は、具体的には式（１）により定義されるP~_nl(ω)を計算する。

k_z,lはz軸方向の波数であり、lはそのインデックスである。波数とは、いわゆる空間周波数又は角度スペクトルのことである。式（１）は、時空間周波数領域への変換の一例であり、他の方法により空間のフーリエ変換を行ってもよい。また、式（０）、（１）を合わせて二次元のDFTを行うような方法でもよい。

＜変換フィルタ部３＞
変換フィルタ部３は、時空間周波数領域信号P~_nl(ω)に対して式（２）により定義されるフィルタF~_nl(ω)を適用してフィルタ処理後信号D~_nl(ω)を生成する（ステップＳ３）。フィルタ処理後信号D~_nl(ω)は、空間周波数逆変換部４に送信される。

式（２）において、cを音速として、k=ω/cは波数である。A(ω)は、所定の複素数である。例えば、A(ω)=1+0×i=1である。また、w_nlはn,lに基づいて例えば以下のように定まる重みである。以下の式において、n_cは、予め定められた値でありnのカットオフ値である。k_cは、予め定められた値でありk_zのカットオフ値である。α_n,α_zは、予め定められた値であり例えば0.05である。もちろん、w_nlとして、他の重み関数を用いてもよい。

H_n ⁽¹⁾(・)はn次の第一種ハンケル関数である。n次の第一種ハンケル関数H_n ⁽¹⁾(x)は、n次の第一種ベッセル関数J_n(x)及び第二種ベッセル関数Y_n(x)を用いて以下のように定義される。

J_n’(・)はn次の第一種ベッセル関数J_n(・)の微分であり、H_n ⁽¹⁾’(・)はn次の第一種ハンケル関数H_n ⁽¹⁾(・)の微分である。J_n’(ψ)及びH_n ⁽¹⁾’(ψ)は以下のように定義される。

＜空間周波数逆変換部４＞
空間周波数逆変換部４は、フィルタ処理後信号D~_nl(ω)を空間の逆フーリエ変換により周波数領域信号D_ij(ω)に変換する（ステップＳ４）。変換された周波数領域信号D_ij(ω)は、周波数逆変換部５に送られる。空間周波数逆変換部４は、具体的には式（３）により定義される周波数領域信号D_ij(ω)を計算する。

＜周波数逆変換部５＞
周波数逆変換部５は、周波数領域信号D_ij(ω)を逆フーリエ変換により時間領域信号P^d _ij(t)に変換する（ステップＳ５）。逆フーリエ変換によりフレーム毎に得られた時間領域信号P^d _ij(t)は適宜シフトされて線形和が取られて、連続した時間領域信号となる。逆フーリエ変換は短時間離散逆フーリエ変換等の既存の方法を用いればよい。時間領域信号P^d _ij(t)は、窓関数部６に送られる。

＜窓関数部６＞
窓関数部６は、時間領域信号P^d _ij(t)に窓関数を乗じて窓関数後時間領域信号d_ij（t）を生成する（ステップＳ６）。窓関数後時間領域信号d_ij（t）は、スピーカＳｉ−ｊ，Ｓ２−１，…，ＳＮ_φ−Ｎ_ｚに送られる。

窓関数として、以下の式より定義されるいわゆるターキー（Tukey）窓関数γ_jを例えば用いる。N_tprは、テーパーを適用する点数であり１以上N_φ,N_z以下の整数である。もちろん、他の窓関数を用いてもよい。

スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_φ−Ｎ_ｚは、窓関数後時間領域信号d_ij（t）に基づいて音を再生する。具体的には、ｉ＝１，…，Ｎ_φ，ｊ＝１，…，Ｎ_ｚとして、スピーカＳｉ−ｊが窓関数後時間領域信号d_ij（t）に基づいて音を再生する。これにより、第一の空間の音場を第二の空間に再現することができる。

このように、マイクロホンアレー及びスピーカアレーを円筒状とすることにより、マイクロホンアレー及びスピーカアレーをそれぞれ構成するマイクロホン及びスピーカの数が少ない場合であっても、左右方向に密な素子配置とし、上下方向は粗な素子配置として、上下感を再現することができる。したがって、従来よりも高精度に音場を再現することができる。

＜理論的背景＞
以下、フィルタF~_nl(ω)が式（２）のように表される理由について説明する。

再現領域の位置ベクトルをr’=(r,φ,z)とし、二次音源が配置されている円筒面をSとし、S上の位置ベクトルをr_s’=(R_s,φ_s,z_s)とする。周波数ωにおいて、再現領域における位置r’の音圧をP(r’,ω)とし、位置r_s’の二次音源の信号をD(r_s’,ω)、r_s’からr’までの伝達関数をG(r’-r_s’,ω)とする。このとき、二次音源によって再現領域内に合成される音場は、以下のように書ける。

この式は関数D(・)とG(・) との変数φ，zに関する畳み込み演算に他ならない。したがって、畳み込みの定理より、式（４）はヘリカル波スペクトル領域において、以下のように表せる。

ここで、二次音源がモノポール特性として近似できる場合について、式（５）の展開を行う。このとき、

であるから、

となる。ここで、

である。Hankel関数の加法定理より、

であるから、

となる。式（６）を式（５）に代入すると、

となる。

理想音場は、無指向性マイクロホンアレーなどを用いて、円筒形状の剛体バッフル上の音圧から与えることができる。収音側の円筒面状の音圧分布をP(r_m’,ω)とすると、P(r’,ω)は入射音場P_i(r’,ω)と散乱音場P_s(r’,ω)との和で書ける。

剛体バッフル上で音圧勾配が０となるという境界条件より、

となる。P_i(r’,ω)とP_s(r’,ω)をヘリカル波スペクトル領域で書くと、

となり、式（８）の境界条件より、以下の関係式が成り立つ。

したがって、

となる。再現すべき音場は入射音場P_i(・)であるから、理想音場はヘリカル波スペクトル領域で以下のように書ける。

ここで求めたいのは、円筒面状の音圧分布P(r_m’,ω)から、二次音源の駆動信号D(r_s’,ω) への変換式である。二次音源をモノポール特性と仮定して変換式を求める場合、式（７）と式（９）とを連立して解けば良い。

フィルタ形式で書くと、以下のようになる。

ただし、

である。

なお、フィルタの特性として本質的には変わらないため、この式を４倍し、周波数特性を調整するための所定の定数であるA(ω)、及び、エバネッセント波を減衰させるための重みw_nlを乗算してもよい。そうすると、この式は、式（２）と一致する。

［変形例等］
音場収音再生装置を構成する各部は、第一の空間に配置された収音装置と第二の空間に配置された再生装置の何れに備えられていてもよい。換言すれば、周波数変換部１、空間周波数変換部２、変換フィルタ部３、空間周波数逆変換部４、周波数逆変換部５、窓関数部６のそれぞれの処理は、第一の空間に配置された収音装置で実行されてもよいし、第二の空間に配置された再生装置で実行されてもよい。収音装置で生成された信号は、再生装置に送信される。

第一の空間と第二の空間の位置は、図２に示したものに限定されない。第一の空間と第二の空間は、隣接していても互いに離れた位置にあってもよい。また、第一の空間と第二の空間の向きもどのようなものであってもよい。

窓関数部６による窓関数の処理は、どの段階で行ってもよいし、多段で行ってもよい。すなわち、窓関数部６は、マイクアレーと周波数変換部１との間、周波数変換部１と空間周波数変換部２との間、空間周波数変換部２と変換フィルタ部３との間、変換フィルタ部３と空間周波数逆変換部４との間、空間周波数逆変換部４と周波数逆変換部５との間の少なくとも１つの間に備えられていてもよい。音場収音再生装置の各部は、その各部に入力される信号について窓関数の処理が行われた場合には、その入力される信号に代えて上記と同様にしてその窓関数の処理がされた後の信号に対して処理を行う。

また、窓関数部６はなくてもよい。この場合、ｉ＝１，…，Ｎ_φ，ｊ＝１，…，Ｎ_ｚとして、スピーカＳｉ−ｊが時間領域信号P^d _ij(t)に基づいて音を再生する。

音場収音再生装置は、変換フィルタ部３を含みさえすれば、他の部を備えていなくてもよい。例えば、音場収音再生装置は、変換フィルタ部３、空間周波数逆変換部４及び周波数逆変換部５から構成されていてもよい。また、音場収音再生装置は、周波数変換部１、空間周波数変換部２及び変換フィルタ部３から構成されていてもよい。

周波数変換部１の処理と空間周波数変換部２の処理とを同時に行ってもよい。同様に、空間周波数逆変換部４の処理と周波数逆変換部５の処理とを同時に行ってもよい。また、空間周波数変換部２と空間周波数逆変換部４とを入れ替えてもよい。

音場収音再生装置は、コンピュータによって実現することができる。この場合、この装置の各部の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、この装置における各部がコンピュータ上で実現される。

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、これらの装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。

この発明は、上述の実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。

１ 周波数変換部
２ 空間周波数変換部
３ 変換フィルタ部
４ 空間周波数逆変換部
５ 周波数逆変換部
６ 窓関数部

Claims (6)

1. 半径R_bの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする半径R_mの２個以上の円のそれぞれに少なくとも２個のマイクロホンが配置されているとし、R_m>R_bとし、上記バッフルの軸方向をz軸方向とし、上記バッフルの周方向とφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、nをφ方向の次数とし、k_z,lをz軸方向の波数とし、lをそのインデックスとし、J_n(・)をn次の第一種ベッセル関数とし、H_n ⁽¹⁾(・)をn次の第一種ハンケル関数とし、J_n’(・)をn次の第一種ベッセル関数J_n(・)の微分とし、H_n ⁽¹⁾’(・)をn次の第一種ハンケル関数H_n ⁽¹⁾(・)の微分とし、A(ω)を所定の複素数とし、w_nlをn,lに基づいて定まる重みとし、少なくとも４個のスピーカが上記第一の空間と異なる第二の空間の半径R_sの仮想円筒の周面に配置されているとして、
   上記マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~_nl(ω)に対して次式により定義されるフィルタF~_nl(ω)を適用してフィルタ処理後信号D~_nl(ω)を生成する変換フィルタ部と、
   

   

   空間の逆フーリエ変換により、上記フィルタ処理後信号D~_nl(ω)を周波数領域信号に変換する空間周波数逆変換部と、
   上記周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号を上記スピーカに出力する周波数逆変換部と、
   を含む音場収音再生装置。
2. 半径R_bの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする半径R_mの２個以上の円のそれぞれに少なくとも２個のマイクロホンが配置されているとし、R_m>R_bとし、上記バッフルの軸方向をz軸方向とし、上記バッフルの円筒の周方向とφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、nをφ方向の次数とし、k_z,lをz軸方向の波数とし、lをそのインデックスとし、J_n(・)をn次の第一種ベッセル関数とし、H_n ⁽¹⁾(・)をn次の第一種ハンケル関数とし、J_n’(・)をn次の第一種ベッセル関数J_n(・)の微分とし、H_n ⁽¹⁾’(・)をn次の第一種ハンケル関数H_n ⁽¹⁾(・)の微分とし、A(ω)を所定の複素数とし、w_nlをn,lに基づいて定まる重みとし、少なくとも４個のスピーカが上記第一の空間と異なる第二の空間の半径R_sの仮想円筒の周面に配置されているとして、
   上記マイクアレーで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換部と、
   空間のフーリエ変換により、上記周波数領域信号を時空間周波数領域信号P~_nl(ω)に変換する空間周波数変換部と、
   上記時空間周波数領域信号P~_nl(ω)に対して次式により定義されるフィルタF~_nl(ω)を適用してフィルタ処理後信号D~_nl(ω)を生成する変換フィルタ部と、
   

   

   を含む音場収音再生装置。
3. 請求項１又は２に記載された音場収音再生装置において、
   上記時空間周波数領域信号P~_nl(ω), D~_nl(ω)と、上記周波数領域信号と、周波数逆変換部により変換された時間領域信号との少なくともひとつは、所定の窓関数により窓関数処理が行われた信号である、
   音場収音再生装置。
4. 半径R_bの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする半径R_mの２個以上の円のそれぞれに少なくとも２個のマイクロホンが配置されているとし、R_m>R_bとし、上記バッフルの軸方向をz軸方向とし、上記バッフルの周方向とφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、nをφ方向の次数とし、k_z,lをz軸方向の波数とし、lをそのインデックスとし、J_n(・)をn次の第一種ベッセル関数とし、H_n ⁽¹⁾(・)をn次の第一種ハンケル関数とし、J_n’(・)をn次の第一種ベッセル関数J_n(・)の微分とし、H_n ⁽¹⁾’(・)をn次の第一種ハンケル関数H_n ⁽¹⁾(・)の微分とし、A(ω)を所定の複素数とし、w_nlをn,lに基づいて定まる重みとし、少なくとも４個のスピーカが上記第一の空間と異なる第二の空間の半径R_sの仮想円筒の周面に配置されているとして、
   変換フィルタ部が、上記マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~_nl(ω)に対して次式により定義されるフィルタF~_nl(ω)を適用してフィルタ処理後信号D~_nl(ω)を生成する変換フィルタステップと、
   

   

   空間周波数逆変換部が、空間の逆フーリエ変換により、上記フィルタ処理後信号D~_nl(ω)を周波数領域信号に変換する空間周波数逆変換ステップと、
   周波数逆変換部が、上記周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号を上記スピーカに出力する周波数逆変換ステップと、
   を含む音場収音再生方法。
5. 半径R_bの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする半径R_mの２個以上の円のそれぞれに少なくとも２個のマイクロホンが配置されているとし、R_m>R_bとし、上記バッフルの軸方向をz軸方向とし、上記バッフルの周方向とφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、nをφ方向の次数とし、k_z,lをz軸方向の波数とし、lをそのインデックスとし、J_n(・)をn次の第一種ベッセル関数とし、H_n ⁽¹⁾(・)をn次の第一種ハンケル関数とし、J_n’(・)をn次の第一種ベッセル関数J_n(・)の微分とし、H_n ⁽¹⁾’(・)をn次の第一種ハンケル関数H_n ⁽¹⁾(・)の微分とし、A(ω)を所定の複素数とし、w_nlをn,lに基づいて定まる重みとし、少なくとも４個のスピーカが上記第一の空間と異なる第二の空間の半径R_sの仮想円筒の周面に配置されているとして、
   周波数変換部が、上記マイクアレーで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換ステップと、
   空間周波数変換部が、空間のフーリエ変換により、上記周波数領域信号を時空間周波数領域信号P~_nl(ω)に変換する空間周波数変換ステップと、
   変換フィルタ部が、上記時空間周波数領域信号P~_nl(ω)に対して次式により定義されるフィルタF~_nl(ω)を適用してフィルタ処理後信号D~_nl(ω)を生成する変換フィルタステップと、
   

   

   を含む音場収音再生方法。
6. 請求項１から３の何れかに記載された音場収音再生装置の各部としてコンピュータを機能させるための音場収音再生プログラム。

Images