(57)【要約】      （修正有）【課題】円形のマイクロホンアレー・スピーカアレーを用いる場合に、従来よりも高精度に音場を再現することができる音場収音再生技術を提供する。【解決手段】半径Rbの円筒形状の剛体のバッフルの軸を中心としバッフルの周方向を円周方向とする第一空間の半径Rmの円周に少なくとも２個のマイクロホンが配置されているとし、少なくとも２個のスピーカが第一の空間と異なる第二の空間の半径Rsの仮想円筒の円周に配置されているとし、スピーカが配置される円を二次音源円とし、Rrefを二次音源円の中心から振幅を一致させる円周上の点までの距離として、マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~m(ω)に対して特定の式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する。【選択図】図１
2014150416.tif 000002
【特許請求の範囲】
【請求項１】
半径Rbの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする第一空間の半径Rmの円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記バッフルの周方向をφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、H(1)m’(・)をm次の第一種ハンケル関数H(1)m(・)の微分とし、Jm’(・)をm次のベッセル関数Jm(・)の微分とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想円筒の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
上記マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【数４４】
2014150416.tif 000046
空間の逆フーリエ変換により、上記フィルタ処理後信号D~m(ω)を周波数領域信号に変換する空間周波数逆変換部と、
上記周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号を上記スピーカに出力する周波数逆変換部と、
を含む音場収音再生装置。
【請求項２】
半径Rbの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする第一空間の半径Rmの円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記バッフルの周方向をφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、H(1)m’(・)をm次の第一種ハンケル関数H(1)m(・)の微分とし、Jm’(・)をm次のベッセル関数Jm(・)の微分とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想円筒の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
上記マイクロホンで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換部と、
空間のフーリエ変換により、上記周波数領域信号を時空間周波数領域信号P~m(ω)に変換する空間周波数変換部と、
上記時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【数４５】
2014150416.tif 000047
を含む音場収音再生装置。
【請求項３】
半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円の円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記マイクロホンが配置される円を受音円とし、上記受音円の周方向をφ方向とし、上記受音円に対して垂直な平面と上記バッフルとの大円の周方向をθ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、nをθ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、h(1)n(・)をn次の第一種球ハンケル関数とし、jn(・)をn次の球ベッセル関数とし、h(1)n’(・)をn次の第一種球ハンケル関数h(1)n(・)の微分とし、jn’(・)をn次の球ベッセル関数jn(・)の微分とし、Pmn(・)をルジャンドル陪関数とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想球の大円の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
上記マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【数４６】
2014150416.tif 000048
空間の逆フーリエ変換により、上記フィルタ処理後信号D~m(ω)を周波数領域信号に変換する空間周波数逆変換部と、
上記周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号を上記スピーカに出力する周波数逆変換部と、
を含む音場収音再生装置。
【請求項４】
半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円の円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記マイクロホンが配置される円を受音円とし、上記受音円の周方向をφ方向とし、上記受音円に対して垂直な平面と上記バッフルとの大円の周方向をθ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、nをθ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、h(1)n(・)をn次の第一種球ハンケル関数とし、jn(・)をn次の球ベッセル関数とし、h(1)n’(・)をn次の第一種球ハンケル関数h(1)n(・)の微分とし、jn’(・)をn次の球ベッセル関数jn(・)の微分とし、Pmn(・)をルジャンドル陪関数とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想球の大円の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
上記マイクロホンで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換部と、
空間のフーリエ変換により、上記周波数領域信号を時空間周波数領域信号P~m(ω)に変換する空間周波数変換部と、
上記時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【数４７】
2014150416.tif 000049
を含む音場収音再生装置。
【請求項５】
半径Rbの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする第一空間の半径Rmの円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記バッフルの周方向をφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、H(1)m’(・)をm次の第一種ハンケル関数H(1)m(・)の微分とし、Jm’(・)をm次のベッセル関数Jm(・)の微分とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想円筒の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
変換フィルタ部が、上記マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタステップと、
【数４８】
2014150416.tif 000050
空間周波数逆変換部が、空間の逆フーリエ変換により、上記フィルタ処理後信号D~m(ω)を周波数領域信号に変換する空間周波数逆変換ステップと、
周波数逆変換部が、上記周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号を上記スピーカに出力する周波数逆変換ステップと、
を含む音場収音再生方法。
【請求項６】
半径Rbの円筒形状の剛体のバッフルの軸を中心とし上記バッフルの周方向を円周方向とする第一空間の半径Rmの円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記バッフルの周方向をφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、H(1)m’(・)をm次の第一種ハンケル関数H(1)m(・)の微分とし、Jm’(・)をm次のベッセル関数Jm(・)の微分とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想円筒の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
周波数変換部が、上記マイクロホンで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換ステップと、
空間周波数変換部が、空間のフーリエ変換により、上記周波数領域信号を時空間周波数領域信号P~m(ω)に変換する空間周波数変換ステップと、
変換フィルタ部が、上記時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタステップと、
【数４９】
2014150416.tif 000051
を含む音場収音再生方法。
【請求項７】
半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円の円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記マイクロホンが配置される円を受音円とし、上記受音円の周方向をφ方向とし、上記受音円に対して垂直な平面と上記バッフルとの大円の周方向をθ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、nをθ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、h(1)n(・)をn次の第一種球ハンケル関数とし、jn(・)をn次の球ベッセル関数とし、h(1)n’(・)をn次の第一種球ハンケル関数h(1)n(・)の微分とし、jn’(・)をn次の球ベッセル関数jn(・)の微分とし、Pmn(・)をルジャンドル陪関数とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想球の大円の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
変換フィルタ部が、上記マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタステップと、
【数５０】
2014150416.tif 000052
空間周波数逆変換部が、空間の逆フーリエ変換により、上記フィルタ処理後信号D~m(ω)を周波数領域信号に変換する空間周波数逆変換ステップと、
周波数逆変換部が、上記周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号を上記スピーカに出力する周波数逆変換ステップと、
を含む音場収音再生方法。
【請求項８】
半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円の円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、上記マイクロホンが配置される円を受音円とし、上記受音円の周方向をφ方向とし、上記受音円に対して垂直な平面と上記バッフルとの大円の周方向をθ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、nをθ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、h(1)n(・)をn次の第一種球ハンケル関数とし、jn(・)をn次の球ベッセル関数とし、h(1)n’(・)をn次の第一種球ハンケル関数h(1)n(・)の微分とし、jn’(・)をn次の球ベッセル関数jn(・)の微分とし、Pmn(・)をルジャンドル陪関数とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが上記第一の空間と異なる第二の空間の半径Rsの仮想球の大円の円周に配置されているとし、上記スピーカが配置される円を二次音源円とし、Rrefを上記二次音源円の中心から振幅を一致させる円周上の点までの距離として、
周波数変換部が、上記マイクロホンで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換ステップと、
空間周波数変換部が、空間のフーリエ変換により、上記周波数領域信号を時空間周波数領域信号P~m(ω)に変換する空間周波数変換ステップと、
変換フィルタ部が、上記時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタステップと、
【数５１】
2014150416.tif 000053
を含む音場収音再生方法。
【請求項９】
請求項１から４の何れかに記載された音場収音再生装置の各部としてコンピュータを機能させるための音場収音再生プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、ある音場に設置されたマイクロホンで音信号を収音し、その音信号を用いてスピーカでその音場を再現する波面合成法（Wave Field Synthesis），アンビソニックス(Ambisonics)の技術に関する。
【背景技術】
【０００２】
波面合成法及びアンビソニックスは、複数のマイクロホンとスピーカを用いて、遠隔地の音場を仮想的に再現する技術である。そのような技術として例えば非特許文献１に記載された技術が知られている。遠隔コミュニケーションシステムなどの応用では、リアルタイムの収音・再現が必要になるため、一般的なマイクロホンアレーで収音した音圧を、一般的なスピーカアレーで出力するための音場再現信号へと一意に変換可能であることが必要となる。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】小山,古家,日和崎,羽田,鈴木,“円筒状マイクロホン・スピーカアレーのための波面再構成フィルタ”,2012年9月,音響学会秋季研究発表会講演論文集,pp.605-608
【発明の概要】
【発明が解決しようとする課題】
【０００４】
非特許文献１に記載された技術では、円筒状のマイクロホンアレー・スピーカアレーを用いることを仮定して、変換のためのフィルタが導出されていた。そのため、このフィルタを円形のマイクロホンアレー・スピーカアレーに適用すると高精度に音場を再現することができず、聴取者の周囲全方位から音が到来するような音場を高精度に再現することができない可能性があった。
【０００５】
この発明の目的は、円形のマイクロホンアレー・スピーカアレーを用いる場合に、従来よりも高精度に音場を再現することができ、聴取者の周囲全方位から音が到来するような音場を高精度に再現することができる音場収音再生装置、方法及びプログラムを提供することである。
【課題を解決するための手段】
【０００６】
上記の課題を解決するために、この発明の一態様による音場収音再生装置は、半径Rbの円筒形状の剛体のバッフルの軸を中心としバッフルの周方向を円周方向とする第一空間の半径Rmの円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、バッフルの周方向をφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、H(1)m’(・)をm次の第一種ハンケル関数H(1)m(・)の微分とし、Jm’(・)をm次のベッセル関数Jm(・)の微分とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが第一の空間と異なる第二の空間の半径Rsの仮想円筒の円周に配置されているとし、スピーカが配置される円を二次音源円とし、Rrefを二次音源円の中心から振幅を一致させる円周上の点までの距離として、マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【０００７】
【数１】
2014150416.tif 000003
【０００８】
空間の逆フーリエ変換により、フィルタ処理後信号D~m(ω)を周波数領域信号に変換する空間周波数逆変換部と、周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号をスピーカに出力する周波数逆変換部と、を含む。
【０００９】
この発明の他の一態様による音場収音再生装置は、半径Rbの円筒形状の剛体のバッフルの軸を中心としバッフルの周方向を円周方向とする第一空間の半径Rmの円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、バッフルの周方向をφ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、H(1)m’(・)をm次の第一種ハンケル関数H(1)m(・)の微分とし、Jm’(・)をm次のベッセル関数Jm(・)の微分とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが第一の空間と異なる第二の空間の半径Rsの仮想円筒の円周に配置されているとし、スピーカが配置される円を二次音源円とし、Rrefを二次音源円の中心から振幅を一致させる円周上の点までの距離として、マイクロホンで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換部と、空間のフーリエ変換により、周波数領域信号を時空間周波数領域信号P~m(ω)に変換する空間周波数変換部と、時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【００１０】
【数２】
2014150416.tif 000004
【００１１】
を含む。
【００１２】
この発明の他の一態様による音場収音再生装置は、半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円の円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、マイクロホンが配置される円を受音円とし、受音円の周方向をφ方向とし、受音円に対して垂直な平面とバッフルとの大円の周方向をθ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、nをθ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、h(1)n(・)をn次の第一種球ハンケル関数とし、jn(・)をn次の球ベッセル関数とし、h(1)n’(・)をn次の第一種球ハンケル関数h(1)n(・)の微分とし、jn’(・)をn次の球ベッセル関数jn(・)の微分とし、Pmn(・)をルジャンドル陪関数とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが第一の空間と異なる第二の空間の半径Rsの仮想球の大円の円周に配置されているとし、スピーカが配置される円を二次音源円とし、Rrefを二次音源円の中心から振幅を一致させる円周上の点までの距離として、マイクロホンで収音された信号に基づいて生成された時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【００１３】
【数３】
2014150416.tif 000005
【００１４】
空間の逆フーリエ変換により、フィルタ処理後信号D~m(ω)を周波数領域信号に変換する空間周波数逆変換部と、周波数領域信号を逆フーリエ変換により時間領域信号に変換し、変換された時間領域信号をスピーカに出力する周波数逆変換部と、を含む。
【００１５】
この発明の他の一態様による音場収音再生装置は、半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円の円周に少なくとも２個のマイクロホンが配置されているとし、Rm≧Rbとし、マイクロホンが配置される円を受音円とし、受音円の周方向をφ方向とし、受音円に対して垂直な平面とバッフルとの大円の周方向をθ方向とし、jを虚数単位とし、ωを周波数とし、cを音速とし、k=ω/cとし、mをφ方向の次数とし、nをθ方向の次数とし、H(1)m(・)をm次の第一種ハンケル関数とし、Jm(・)をm次のベッセル関数とし、h(1)n(・)をn次の第一種球ハンケル関数とし、jn(・)をn次の球ベッセル関数とし、h(1)n’(・)をn次の第一種球ハンケル関数h(1)n(・)の微分とし、jn’(・)をn次の球ベッセル関数jn(・)の微分とし、Pmn(・)をルジャンドル陪関数とし、wmをmに基づいて定まる重みとし、A(ω)をωに基づいて定まる複素数とし、少なくとも２個のスピーカが第一の空間と異なる第二の空間の半径Rsの仮想球の大円の円周に配置されているとし、スピーカが配置される円を二次音源円とし、Rrefを二次音源円の中心から振幅を一致させる円周上の点までの距離として、マイクロホンで収音された信号をフーリエ変換により周波数領域信号に変換する周波数変換部と、空間のフーリエ変換により、周波数領域信号を時空間周波数領域信号P~m(ω)に変換する空間周波数変換部と、時空間周波数領域信号P~m(ω)に対して次式により定義されるフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する変換フィルタ部と、
【００１６】
【数４】
2014150416.tif 000006
【００１７】
を含む。
【発明の効果】
【００１８】
円形のマイクロホンアレー・スピーカアレー用のフィルタを用いて、円形のマイクロホンアレーの収音信号を円形のスピーカアレーの駆動信号へ変換し、音場を再現することが可能であることから、聴取者の周囲全方位から音が到来するような音場を高精度に再現することが可能となる。
【図面の簡単な説明】
【００１９】
【図１】マイクロホン及びスピーカの配置の例を説明するための図。
【図２】音場収音再生装置の例を示す機能ブロック図。
【図３】音場収音再生方法の例を示す流れ図。
【図４】マイクロホン及びスピーカの配置の例を説明するための図。
【発明を実施するための形態】
【００２０】
［第一実施形態］
以下、図面を参照してこの発明の実施形態を説明する。以下の説明において、テキスト中で使用する記号「~」、「−」等は、本来直前の文字の真上に記載されるべきものであるが、テキスト記法の制限により、当該文字の直後に記載する。式中においてはこれらの記号は本来の位置に記述している。また、ベクトルや行列の各要素単位で行われる処理は、特に断りが無い限り、そのベクトルやその行列の全ての要素に対して適用されるものとする。
【００２１】
＜マイクロホンアレー及びスピーカアレーの配置について＞
第一実施形態の音場収音再生装置及び方法は、図１に示すように、第一の空間の半径Rbの円筒形状の剛体のバッフルの軸を中心とし、バッフルの周方向を円周方向とする半径Rmの円の円周に配置されているNch個のマイクロホンＭ１，Ｍ２，…，ＭＮｃｈで構成されるマイクロホンアレーと、第二の空間の半径Rsの仮想円筒の円周に配置されているNch個のスピーカＳ１，Ｓ２，…，ＳＮｃｈで構成されるスピーカアレーとを用いて、第一の空間の音源Ｓで発生した音によって形成された第一の空間の音場を第二の空間で再現する。なお、仮想円筒の円周とは、仮想円筒の軸に対して垂直な平面と仮想円筒との共通部分である円の円周を意味する。第一の空間及び第二の空間は互いに異なる空間である。図１では、第二の空間で再現された音源Ｓを音源Ｓ'と表現している。マイクロホンアレーが形成する円を受音円ともいい、スピーカアレーが形成する円を二次音源円ともいう。言い換えると、マイクロホンアレーを配置する円を受音円といい、スピーカアレーを配置する円を二次音源円という。バッフルの軸方向をz方向とし、半径Rbのバッフルの周方向をφ方向とする。言い換えると、φ方向とは、受音円の周方向である。第一の空間に配置されたマイクロホンの数と第二の空間に配置されたスピーカの数は異なっていてもよい。マイクロホンの数が、第二の空間に配置されたスピーカの数よりも多い場合には、再生信号を間引けばよい。一方、マイクロホンの数が、第二の空間に配置されたスピーカの数よりも少ない場合には、再生信号をチャネル間で平均を取るなどして補間を行えばよい。補間を行う方法は、例えば、線形補間やsinc補間などを適用することができる。
【００２２】
マイクロホンアレーは、２個以上のマイクロホンを使って円形を構成する。例えば、図１に示すように、円筒形状の剛体のバッフルの軸を中心とし、バッフルの周方向を円周方向とする半径Rmの円の円周にNch個のマイクロホンが等間隔に配置される。Nchは予め定められた２以上の整数である。
【００２３】
円筒形状の剛体のバッフルの半径Rbとバッフルの軸とマイクロホンの距離RmとはRm≧Rbであればどのような値でもよい。Rm=Rbの場合には、マイクロホンは半径Rb(=Rm)の円筒形状のバッフルの周面の表面に配置されることになる。Rm>Rbの場合には、図１に示すように、マイクロホンは半径Rbの円筒形状のバッフルの軸からRm離れた位置に配置されることになる。言い替えると、マイクロホンはバッフルの周面の表面から(Rm-Rb)離れた位置に配置される。例えば、バッフルの周面から垂直に突き出す細い棒状の部材により支持することで、マイクロホンは配置される。
【００２４】
なお、変換フィルタ部４において用いるフィルタF~m(ω)は、バッフルの円筒が無限長であることを想定して、導出している（詳細は後述する＜フィルタの導出について＞を参照）。そのため、（１）バッフルの円筒は長いものの方が好ましく、（２）マイクロホンの設置位置はバッフルの円筒の長さの中央が好ましい。なお、本実施形態は、バッフルの長さ、マイクロホンの設置位置を、（１）（２）に限定するものではないが、（１）（２）の場合に最も収音場に近い波面を合成することができる。
【００２５】
半径Rmの円の円周に配置されたNch個のマイクロホンは、φcを所定の角度として、それぞれφcradianの間隔で位置している。なお、マイクロホンは、ほぼ等間隔に配置されていれば、厳密に等間隔に配置されている必要はない。すなわち、隣接するマイクロホンとの間隔であるφcのそれぞれは、厳密に同じ値である必要はなく、ほぼ同じ値であればよい。また、マイクロホンは、どのような間隔で配置されてもよい。すなわち、隣接するマイクロホンとの間隔であるφcは、任意の値を取ることができる。ただし、マイクロホンをほぼ等間隔に配置する、すなわち隣接するマイクロホンとの間隔であるφcをほぼ同じ値とすることで、高精度に音場を再現することができる。
【００２６】
スピーカもマイクロホンと同様に配置される。すなわち、スピーカアレーは、２個以上のスピーカを使って円形を構成する。例えば、図１に示すように、仮想円筒の円周にNch個のスピーカが等間隔に配置される。Nchは予め定められた２以上の整数である。
【００２７】
仮想円筒の円周に配置されたNch個のスピーカは、φcを所定の角度として、それぞれφcradianの間隔で位置している。なお、スピーカは、ほぼ等間隔に配置されていれば、厳密に等間隔に配置されている必要はない。すなわち、隣接するスピーカとの間隔であるφcのそれぞれは、厳密に同じ値である必要はなく、ほぼ同じ値であればよい。また、スピーカは、どのような間隔で配置されてもよい。すなわち、隣接するスピーカとの間隔であるφcは、任意の値を取ることができる。ただし、スピーカをほぼ等間隔に配置する、すなわち隣接するスピーカとの間隔であるφcをほぼ同じ値とすることで、高精度に音場を再現することができる。
【００２８】
受音円の円周のNch個のマイクロホンの位置関係は、二次音源円の対応するNch個のスピーカの位置と同じであることが望ましいが、異なっていても良い。この位置が同じであれば、より忠実に音場の再生を行うことができる。
【００２９】
受音円の半径Rmは、例えば0.2m程度とする。また、二次音源円の半径Rsは、例えば、1.5m程度とする。なお、本実施形態では、Rs≧Rmであると仮定するが、これを満たさなくても（つまりRs＜Rmの場合にも）同様の処理を適用可能である。ただし、Rs≧Rmのときに音場再生の精度が良くなる。マイクロホンアレーが形成する円の半径Rm及び仮想円筒の半径Rsは、それぞれ値が大きいほど広い領域を再現することができるが、より多くのマイクロホン及びスピーカが必要となる。マイクロホンアレーが形成する円の半径Rm及び仮想円筒の半径Rsは、それぞれ収音及び再生する信号の周波数を考慮して実験的に設定することが望ましい。また、マイクロホンは受音円の外側に、スピーカは二次音源円の内側に向けて配置される。マイクロホンアレーが対象とする収音範囲は受音円の外側であり、スピーカアレーが対象とする再生範囲は二次音源円の内側である。
【００３０】
スピーカは、音響的に透明な状態で第二の空間の空中に配置されてもよいし、音響的に透明でない状態で第二の空間に配置されてもよい。音響的に透明な状態とは、スピーカが配置されていない第二の空間の伝達特性と同じ伝達特性を保った状態ということである。例えば、スピーカを音響的に透明な状態で第二の空間の空中に配置したい場合には、スピーカは、糸で吊るされるか、細い棒で固定されることにより、第二の空間の空中に配置される。
【００３１】
第一の空間のマイクロホンＭｉ(i=1,2,…,Nch)の位置を円筒座標系でr-m,i=(Rm,φm,i,0)と表現する。第二の空間のスピーカＳｉ(i=1,2,…,Nch)の位置を円筒座標系でr-s,i=(Rs,φs,i,0)と表現する。
【００３２】
＜音場収音再生装置＞
第一実施形態の音場収音再生装置は、図２に示すように周波数変換部１、空間周波数変換部３、変換フィルタ部４、空間周波数逆変換部５及び周波数逆変換部６を例えば含み、図３に例示された各ステップの処理を行う。
【００３３】
第一の空間に配置されたマイクロホンＭ１，Ｍ２，…，ＭＮchは、第一の空間の音源Ｓで発せられた音を収音して時間領域の信号を生成する。生成された信号は、周波数変換部１に送られる。r-m,i=(Rm,φm,i,0)に位置するマイクロホンＭｉで収音された時間領域の時刻ｔの信号をpi(t)と表記する。
【００３４】
＜周波数変換部１＞
周波数変換部１は、マイクロホンＭiで収音された信号pi(t)をフーリエ変換により周波数領域信号Pi(ω)に変換する（ステップＳ１）。生成された周波数領域信号Pi(ω)は、空間周波数変換部３に送られる。ωは周波数である。なお、ωを音速cで割った値k=ω/cを波数として定義する。波数とは、いわゆる空間周波数又は角度スペクトルのことである。例えば、短時間離散フーリエ変換により周波数領域信号Pi(ω)が生成される。もちろん、他の既存の方法により周波数領域信号Pi(ω)を生成してもよい。また、オーバーラップアド等の方法を用いて周波数領域信号Pi(ω)を生成してもよい。入力信号が長い場合や、リアルタイム処理のように連続して信号が入力される場合には、例えば10msごとといったフレームごとに処理を行う。周波数領域信号Pi(ω)は、例えば式（１）のように定義される。関数expの引数の中のjは虚数単位である。
【００３５】
【数５】
2014150416.tif 000007
【００３６】
＜空間周波数変換部３＞
空間周波数変換部３は、空間のフーリエ変換により周波数領域信号Pi(ω)を時空間周波数領域信号P~m(ω)に変換する（ステップＳ３）。時空間周波数領域信号P~m(ω)は、周波数ωごとに計算される。変換された時空間周波数領域信号P~m(ω)は、変換フィルタ部４に送られる。空間周波数変換部３は、例えば、式（２）により定義されるP~m(ω)を計算する。
【００３７】
【数６】
2014150416.tif 000008
【００３８】
mはφ方向の次数である。ここで、-M≦m≦Mであり、mは整数とする。式（２）は、時空間周波数領域への変換の一例であり、他の方法により空間のフーリエ変換を行ってもよい。また、式（１）、（２）を合わせて二次元のDFTを行うような方法でもよい。
【００３９】
＜変換フィルタ部４＞
変換フィルタ部４は、時空間周波数領域信号P~m(ω)に対して次式により所定のフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する（ステップＳ４）。フィルタ処理後信号D~m(ω)は、空間周波数逆変換部５に送信される。
【００４０】
【数７】
2014150416.tif 000009
【００４１】
Rrefは、二次音源円の中心から振幅を一致させる（二次音源円と同一平面上に存在する）円周上の点までの距離である。従って、Rref=0のとき、二次音源円の中心において振幅が一致し、Rref≠0のとき、半径Rrefの二次音源円の同心円の円周上において振幅が一致する。なお、Rref=0の場合、二次音源円の中心から振幅を一致させる点の集合は円を形成せず、点（二次音源円の中心）となるが、この点も「二次音源円の中心から振幅を一致させる円周」との表現に含まれるものとする。また、Rrefは、二次音源円の中心から振幅を一致させる、二次音源円と同一平面上に存在する点までの距離と言い換えてもよい。式（４）の場合、振幅がRrefで決まる二次音源円の同心円の円周上で目的音場と一致するのに対し、式（５）の場合にはある一点（二次音源円の中心）で振幅が一致する。A(ω)はωに基づいて定まる複素数（ωに依存する任意の複素数）であり、例えば高域のゲインを下げるなどのイコライジング操作を行う値である。式（５）の場合には、具体的には、スピーカアレーと同じ高さであり、二次音源円の中心、または、半径Rrefの二次音源円の同心円の円周上の位置で振幅が一致する。wmは重み関数であり、例えば以下のように定義される。
【００４２】
【数８】
2014150416.tif 000010
【００４３】
【数９】
2014150416.tif 000011
【００４４】
Γ(z),Yv(z)はそれぞれガンマ関数、ノイマン関数を表す。変換フィルタは上式を数値計算することによって得られる。式（４）の積分値などは、次式のように例えば離散化して近似的に計算してもよい。
【００４５】
【数１０】
2014150416.tif 000012
【００４６】
ここで、iは離散化のインデックスであり、Δkzは離散化間隔、Iは総和の範囲を表わす。式（１１）は離散化して近似的に積分値を求める方法の一例であり、他の方法により離散化してもよい。区分求積法の定義で用いられる様々な近似公式を用いて離散化することができる。
【００４７】
また、Rb=Rmのとき、つまり、マイクロホンが半径Rb(=Rm)の円筒形状のバッフルの周面の表面に配置されるとき、式（４）及び式（５）は、それぞれ以下のように簡略化できる。
【００４８】
【数１１】
2014150416.tif 000013
【００４９】
式（４）もしくは式（５）、または、式（４−２）もしくは式（５−２）で定義されるフィルタF~m(ω)を適用することにより、第二の空間で再現される信号の振幅を所定の円周上で一致させることができるため、従来よりも広い範囲で再現される信号の振幅が一致する。
【００５０】
＜空間周波数逆変換部５＞
空間周波数逆変換部５は、フィルタ処理後信号D~m(ω)を空間の逆フーリエ変換により周波数領域信号Di(ω)に変換する（ステップＳ５）。変換された周波数領域信号Di(ω)は、周波数逆変換部６に送られる。空間周波数逆変換部５は、例えば式（１２）により定義される周波数領域信号Di(ω)を計算する。
【００５１】
【数１２】
2014150416.tif 000014
【００５２】
＜周波数逆変換部６＞
周波数逆変換部６は、周波数領域信号Di(ω)を逆フーリエ変換により時間領域信号di(t)に変換する（ステップＳ６）。逆フーリエ変換によりフレーム毎に得られた時間領域信号di(t)は適宜シフトされて線形和が取られて、連続した時間領域信号となる。逆フーリエ変換は短時間離散逆フーリエ変換等の既存の方法を用いればよい。時間領域信号di(t)は、スピーカＳ１，Ｓ２，…，ＳＮchに送られる。
【００５３】
スピーカＳｉは、時間領域信号di(t)に基づいて音を再生する。例えば、ｉ＝１，…，Ｎchとして、r-s,i=(Rs,φs,i,0)に位置するスピーカＳｉが時間領域信号di(t)に基づいて音を再生する。これにより、第一の空間の音場を第二の空間に再現することができる。
【００５４】
上述の通り、マイクロホンの数が、スピーカの数よりも多い場合には、時間領域信号di(t)を間引いてもよい。一方、マイクロホンの数が、スピーカの数よりも少ない場合には、時間領域信号di(t)の平均を取るなどして補間を行ってもよい。補間を行う方法は、例えば、線形補間やsinc補間などを適用することができる。
【００５５】
＜効果＞
円形のマイクロホンアレー・スピーカアレー用のフィルタを用いて、円形のマイクロホンアレーの収音信号を円形のスピーカアレーの駆動信号へ変換し、音場を再現することが可能であることから、聴取者の周囲全方位から音が到来するような音場を高精度に再現することが可能となる。したがって、従来よりも高精度に音場を再現することができる。
【００５６】
＜フィルタの導出について＞
以下、フィルタF~m(ω)が式（４）及び式（５）のように表される理由について説明する。
【００５７】
ここでは円筒座標系r-=(r,φ,z)で考える。収音領域に配置した受音円で取得した音圧分布から、再現領域に配置した二次音源円を用いて、x-y平面上の音場を再現する。
【００５８】
円周上に連続的に配置した二次音源による合成音場Psyn(r-,ω)は、二次音源駆動信号D(r-s,ω)と、二次音源と目的音場との間の伝達関数G(r--r-s,ω)を用いて、以下のように表わされる。
【００５９】
【数１３】
2014150416.tif 000015
【００６０】
これはD(・)とG(・)とのφに関する畳みこみと見なせることから、円調和スペクトル領域で、以下のように積の形で表わせる。
【００６１】
【数１４】
2014150416.tif 000016
【００６２】
一方、目的音場Pdes(r-，ω)を、収音場Prcv(r-,ω)の半径Rmの円周上で取得した音圧分布Prcv(Rm,0,φ,ω)を用いて表わすことが必要となる。受音円と同じ位置を中心として、半径Rb(≦Rm)の無限長の剛体円筒が設置されているとすると、収音場Prcv(r-,ω)は次式のように入射音場Pinc(r-,ω)と散乱音場Psct(r-,ω)との和で表わせる。
Prcv(r-,ω)=Pinc(r-,ω)+Psct(r-,ω)   …(15)
また、剛体円筒上の境界条件より、
【００６３】
【数１５】
2014150416.tif 000017
【００６４】
が成り立つ。また、Pinc(r-,ω)とPsct(r-,ω)は、ヘリカル波スペクトル領域で以下のように書ける。
【００６５】
【数１６】
2014150416.tif 000018
【００６６】
したがって、
【００６７】
【数１７】
2014150416.tif 000019
【００６８】
となる。以上により、収音場がz軸方向に不変であると仮定すると、入射音場と収音場との関係式が以下のように導ける。
【００６９】
【数１８】
2014150416.tif 000020
【００７０】
目的音場は入射音場に等しいことから、
【００７１】
【数１９】
2014150416.tif 000021
【００７２】
となる。
合成音場Psyn(r-,ω)と目的音場Pdes(r-,ω)が一致すればよいことから、
【００７３】
【数２０】
2014150416.tif 000022
【００７４】
ここで簡単のため、G(r--r-s,ω)はモノポール特性であると仮定した。
【００７５】
【数２１】
2014150416.tif 000023
【００７６】
以上で受音円上の音圧分布から二次音源駆動信号の変換式が得られたが、r=Rrefとして、振幅が目的音場と一致する円を予め与える必要がある。
【００７７】
二次音源を線音源として仮定すると、変換式はより簡単化できる。
【００７８】
【数２２】
2014150416.tif 000024
【００７９】
であることより、
【００８０】
【数２３】
2014150416.tif 000025
【００８１】
となる。ただし、実際のスピーカはモノポール特性に近く、線音源近似のための補正が必要となる。
【００８２】
【数２４】
2014150416.tif 000026
【００８３】
であることより、各出力信号に対して、以下の周波数特性の補正を行う。
【００８４】
【数２５】
2014150416.tif 000027
【００８５】
Rb=Rmの場合は以下の関係式を用いて簡単化できる。
【００８６】
【数２６】
2014150416.tif 000028
【００８７】
［第二実施形態］
第一実施形態と異なる部分を中心に説明する。
【００８８】
＜マイクロホンアレー及びスピーカアレーの配置について＞
第二実施形態の音場収音再生装置及び方法は、図４に示すように、半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円の円周に配置されているNch個のマイクロホンＭ１，Ｍ２，…，ＭＮｃｈで構成されるマイクロホンアレーと、第二の空間の半径Rsの仮想球の大円の円周に配置されているNch個のスピーカＳ１，Ｓ２，…，ＳＮｃｈで構成されるスピーカアレーとを用いて、第一の空間の音源Ｓで発生した音によって形成された第一の空間の音場を第二の空間で再現する。マイクロホンアレーが形成する円を受音円ともいい、スピーカアレーが形成する仮想球の大円を二次音源円ともいう。受音円及び二次音源円に対して垂直な方向をz方向とし、受音円及び二次音源円の周方向をφ方向とし、受音円に対して垂直な平面とバッフルとの大円、及び、二次音源円に対して垂直な平面と仮想球との大円の周方向をθ方向とする。
【００８９】
第一実施形態では、半径Rbの円筒形状の剛体のバッフルの軸を中心とし、バッフルの周方向を円周方向とする半径Rmの円及び半径Rsの仮想円筒の円周にそれぞれマイクロホン及びスピーカが配置されている。一方、第二実施形態では、半径Rbの球状の剛体のバッフルの中心を中心とする、第一空間の半径Rmの円及び半径Rsの仮想球の大円の円周にそれぞれマイクロホン及びスピーカが配置されている。
【００９０】
第一の空間のマイクロホンＭｉ(i=1,2,…,Nch)の位置を球座標系でr-m,i=(Rm,π/2,φm,i)と表現する。第二の空間のスピーカＳｉ(i=1,2,…,Nch)の位置を球座標系でr-s,i=(Rs,π/2,φs,i)と表現する。
【００９１】
＜変換フィルタ部４＞
変換フィルタ部４は、時空間周波数領域信号P~m(ω)に対して次式により所定のフィルタF~m(ω)を適用してフィルタ処理後信号D~m(ω)を生成する（ステップＳ４）。フィルタ処理後信号D~m(ω)は、空間周波数逆変換部５に送信される。
【００９２】
【数２７】
2014150416.tif 000029
【００９３】
Rrefは、二次音源円の中心から振幅を一致させる（二次音源円と同一平面上に存在する）円周上の点までの距離である。nはθ方向の次数である。h(1)n(・)、jn(・)、h(1)n’(・)、jn’(・)、Pmn(・)は、それぞれn次の第一種球ハンケル関数、n次の球ベッセル関数、n次の第一種球ハンケル関数h(1)n(・)の微分、n次の球ベッセル関数jn(・)の微分、ルジャンドル陪関数であり、以下のように定義される。
【００９４】
【数２８】
2014150416.tif 000030
【００９５】
Pn(・)はルジャンドル多項式を表す。変換フィルタは上式を数値計算することによって得られる。
【００９６】
また、Rb=Rmのとき、つまり、マイクロホンが半径Rb(=Rm)の円筒形状のバッフルの周面の表面に配置されるとき、式（３１）及び式（３２）は、それぞれ以下のように簡略化できる。
【００９７】
【数２９】
2014150416.tif 000031
【００９８】
式（３１）もしくは式（３２）、または、式（３１−２）もしくは式（３２−２）で定義されるフィルタF~m(ω)を適用することにより、第二の空間で再現される信号の振幅を所定の円周上で一致させることができるため、従来よりも広い範囲で再現される信号の振幅が一致する。
【００９９】
このような構成により、第一実施形態と同様の効果を得ることができる。
【０１００】
＜フィルタの導出について＞
以下、フィルタF~m(ω)が式（３１）及び式（３２）のように表される理由について説明する。
【０１０１】
ここでは球座標系r-=(r,θ,φ)で考える。収音領域に配置した受音円で取得した音圧分布から、再現領域に配置した二次音源円を用いて、x-y平面上の音場を再現する。
【０１０２】
円周上に連続的に配置した二次音源による合成音場Psyn(r-,ω)は、第一実施形態の場合と同じく、円調和スペクトル領域で以下のように書ける。
【０１０３】
【数３０】
2014150416.tif 000032
【０１０４】
一方、目的音場Pdes(r-,ω)を、収音場Prcv(r-,ω)の半径Rmの円周上で取得した音圧分布P(Rm,π/2,φ,ω)を用いて表わすことが必要となる。受音円と同じ位置と中心として、半径Rb(≦Rm)の剛体球が設置されているとすると、収音場Prcv(r-,ω)は入射音場Pinc(r-,ω)と散乱音場Psct(r-,ω)との和で表わせる。
Prcv(r-,ω)=Pinc(r-,ω)+Psct(r-,ω)    (38)
また、剛体球上の境界条件より、
【０１０５】
【数３１】
2014150416.tif 000033
【０１０６】
が成り立つ。また、Pinc(r-,ω)とPsct(r-,ω)は、球調和スペクトル領域で以下のように書ける。
【０１０７】
【数３２】
2014150416.tif 000034
【０１０８】
ここで、Ymn(θ,φ)は、球調和関数を表す。したがって、
【０１０９】
【数３３】
2014150416.tif 000035
【０１１０】
となる。以上により、収音場がz軸方向に普遍であると仮定すると、入射音場と収音場との関係式が以下のように導ける。
【０１１１】
【数３４】
2014150416.tif 000036
【０１１２】
更に、
【０１１３】
【数３５】
2014150416.tif 000037
【０１１４】
であり、
【０１１５】
【数３６】
2014150416.tif 000038
【０１１６】
であることから、
【０１１７】
【数３７】
2014150416.tif 000039
【０１１８】
となる。
合成音場Psys(r-,ω)と目的音場Pdes(r-,ω)が一致すればよいことから、
【０１１９】
【数３８】
2014150416.tif 000040
【０１２０】
ここで簡単のため、G(r--r-s,ω)はモノポール特性であると仮定した。
【０１２１】
【数３９】
2014150416.tif 000041
【０１２２】
以上で、受音円上の音圧分布から二次音源駆動信号の変換式が得られたが、r=Rrefとして、振幅が目的音場と一致する円を予め与える必要がある。
【０１２３】
二次音源を線音源として仮定すると、変換式はより簡単化できる。
【０１２４】
【数４０】
2014150416.tif 000042
【０１２５】
であることより、
【０１２６】
【数４１】
2014150416.tif 000043
【０１２７】
となる。ただし、実際のスピーカはモノポール特性に近く、線音源近似のための補正が必要となる。式(26)より、各出力信号に対して、以下の周波数特性の補正を行う。
【０１２８】
【数４２】
2014150416.tif 000044
【０１２９】
Rb=Rmの場合は、以下の関係式を用いて導出できる。
【０１３０】
【数４３】
2014150416.tif 000045
【０１３１】
［変形例等］
音場収音再生装置を構成する各部は、第一の空間に配置された収音装置と第二の空間に配置された再生装置の何れに備えられていてもよい。換言すれば、周波数変換部１、空間周波数変換部３、変換フィルタ部４、空間周波数逆変換部５、周波数逆変換部６のそれぞれの処理は、第一の空間に配置された収音装置で実行されてもよいし、第二の空間に配置された再生装置で実行されてもよい。収音装置で生成された信号は、再生装置に送信される。
【０１３２】
第一の空間と第二の空間の位置は、図１や図４に示したものに限定されない。第一の空間と第二の空間は、隣接していても互いに離れた位置にあってもよい。また、第一の空間と第二の空間の向きもどのようなものであってもよい。
【０１３３】
音場収音再生装置は、変換フィルタ部４を含みさえすれば、他の部を備えていなくてもよい。例えば、音場収音再生装置は、変換フィルタ部４、空間周波数逆変換部５及び周波数逆変換部６から構成されていてもよい。また、音場収音再生装置は、周波数変換部１、空間周波数変換部３、変換フィルタ部４から構成されていてもよい。
【０１３４】
周波数変換部１の処理と空間周波数変換部３の処理とを同時に行ってもよい。同様に、空間周波数逆変換部５の処理と周波数逆変換部６の処理とを同時に行ってもよい。また、空間周波数変換部３と空間周波数逆変換部５とを入れ替えてもよい。
【０１３５】
音場収音再生装置は、コンピュータによって実現することができる。この場合、この装置の各部の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、この装置における各部がコンピュータ上で実現される。
【０１３６】
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、これらの装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
【０１３７】
この発明は、上述の実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。
【符号の説明】
【０１３８】
１  周波数変換部
３  空間周波数変換部
４  変換フィルタ部
５  空間周波数逆変換部
６  周波数逆変換部
【図１】
2014150416.tif 000054
【図２】
2014150416.tif 000055
【図３】
2014150416.tif 000056
【図４】
2014150416.tif 000057