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

Description

この発明は、ある音場に設置されたマイクアレーで音信号を収音し、その音信号を用いてスピーカアレーでその音場を再生する波面合成法（Wave Field Synthesis）の技術に関する。

波面合成法は、複数のマイクとスピーカを用いて遠隔地の音場を仮想的に再現する技術である。波面合成法を実行する際には、音場を再現するための波面の位置を自在にコントロールできることが重要である。これにより、再生した音場に仮想的なイメージを与えたり、音響効果を加えることができるためである。

非特許文献１に記載された技術においては、図１６に示すマイクアレーＭで収音した音信号を用いて、マイクアレーＭよりも音源Ｓから離れた位置にある波面Ａを再現することができた。

Sascha Spors, Rudolf Rabenstein, and Jens Ahrens, "The Theory of Wave Field Synthesis Revisited", 124th Convention of the Audio Engineering Society Amsterdam, 2008 May 17-20

しかしながら、非特許文献１に記載された技術においては、マイクアレーＭよりも音源Ｓ側にある波面Ｂ（図１６参照）を再現することはできないという課題があった。

上記の課題を解決するために、マイクアレーで収音された音信号をフーリエ変換により周波数領域信号に変換し、周波数領域信号に窓関数を乗じて窓関数後周波数領域信号を生成し、マイクアレーよりも音源に近い位置での波面を模擬したフィルタを窓関数後周波数領域信号に適用し、フィルタを適用した後の信号を逆フーリエ変換により時間領域信号に変換する。

不均質波を無視することにより計算したフィルタを適用することにより、マイクアレーＭよりも音源Ｓ側にある波面Ｂ（図１６参照）を再現することができる。

第一実施形態の音場収音再生装置の例の機能ブロック図。 第一実施形態の音場収音再生装置のマイクアレー及びスピーカアレーの配置の例を説明するための図。 第二実施形態の音場収音再生装置の例の機能ブロック図。 第一実施形態の音場収音再生装置のマイクアレー及びスピーカアレーの配置の例を説明するための図。 音場収音再生方法の例を示す流れ図。 ダイポール特性のマイクの例を説明するための図。 シミュレーションの条件を示すための図。 従来のＷＦＳによってｚ＝０の波面を再現することにより音場を再現したシミュレーション結果を示す図。 図５に対応する理想的な音場を示す図。 この発明によりｚ＝−０．６ｍの波面を再現することにより音場を再現したシミュレーション結果を示す図。 図７に対応する理想的な音場を示す図。 この発明によりｚ＝−２ｍの波面を再現することにより音場を再現したシミュレーション結果を示す図。 図９に対応する理想的な音場を示す図。 この発明により音源が２つある音場を再現したシミュレーション結果を示す図。 図１１に対応する理想的な音場を示す図。 再現しようとする波面の位置を説明するための図。

以下、図面を参照してこの発明の一実施形態を説明する。

［第一実施形態］
第一実施形態の音場収音再生装置及び方法は、図２に示す第一の部屋のｚ＝ｚ_１の位置に配置されたＮ_ｘ×Ｎ_ｙ個のモノポール特性のマイクで構成される二次元マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙで収音した音信号を用いて以下に説明する方法でｚ＝ｚ_２の位置の波面Ｂを計算して、第二の部屋に配置されたＮ_ｘ’×Ｎ_ｙ’個のモノポール特性のスピーカで構成される二次元スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’で再現する。

Ｎ_ｘ，Ｎ_ｙ，Ｎ_ｘ’，Ｎ_ｙ’は任意の整数である。Ｎ_ｘとＮ_ｘ’は異なっていてもよい。また、Ｎ_ｙとＮ_ｙ’も異なっていてもよい。すなわち、マイクアレーを構成するマイクの数とスピーカアレーを構成するスピーカの数は同じでなくてもよく、マイクアレーを構成するマイクとスピーカアレーを構成するスピーカとは必ずしも一対一対応している必要はない。

第一の部屋のｚ＝ｚ_１の位置に配置された二次元マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙを構成する各マイクの位置をｒ_ｓ＝（ｘ_ｉ，ｙ_ｊ，ｚ_１）と表わすことにする。第一の部屋のｚ＝ｚ_２での波面Ｂを再現する第二の部屋に配置された二次元スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’を構成する各スピーカの位置をｒ_ｓ’＝（ｘ_ｎ，ｙ_ｎ，ｚ_２）と表わすことにする。

第一実施形態の音場収音再生装置は、図１に示すように周波数領域変換部１、窓関数部２、フィルタ部３、時間領域変換部４を例えば含み、図５の実線で示された処理を行う。

第一の部屋のｚ＝ｚ_１の位置に配置された二次元マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙは、第一の部屋の音源Ｓで発せられた音を収音して時間領域の音信号を生成する。生成された音信号は、周波数領域変換部１に送られる。ｒ_ｓ＝（ｘ_ｉ，ｙ_ｊ，ｚ_１）のマイクＭｉ−ｊで収音された時間領域の時刻ｔの音信号をｆ（ｉ，ｊ，ｔ）と表記する。

周波数領域変換部１は、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙで収音された音信号ｆ（ｉ，ｊ，ｔ）をフーリエ変換により周波数領域信号Ｆ（ｉ，ｊ，ω）に変換する（ステップＳ１）。生成された周波数領域信号Ｆ（ｉ，ｊ，ω）は、窓関数部２に送られる。ωは周波数である。例えば、短時間フーリエ変換により周波数領域信号Ｆ（ｉ，ｊ，ω）が生成される。もちろん、他の既存の方法により周波数領域信号Ｆ（ｉ，ｊ，ω）を生成してもよい。

窓関数部２は、周波数領域信号Ｆ（ｉ，ｊ，ω）に窓関数を乗じて窓関数後周波数領域信号Ｆ_ｗ（ｊ，ｊ，ω）を生成する（ステップＳ２）。窓関数後周波数領域信号Ｆ_ｗ（ｊ，ｊ，ω）は、フィルタ部３に送られる。窓関数として、以下の式より定義されるいわゆるターキー（Turkey）窓関数ｗ（ｉ，ｊ）を例えば用いる。Ｎ_ｔｐｒは、テーパーを適用する点数であり１以上Ｎ_ｘ，Ｎ_ｙ以下の整数である。

フィルタ部３は、マイクアレーよりも音源に近い位置での波面Ｂを不均質波を無視することにより模擬したフィルタを窓関数後周波数領域信号Ｆ_ｗ（ｊ，ｊ，ω）に適用する（ステップＳ３）。フィルタを適用した後の信号をＧ（ｎ，ｍ，ω）と表記すると、フィルタ部３の処理は次式のように表わされる。ｋは波数であり、ｃを音速とすると、ｋ＝ω／ｃである。ｅはネピア数、πは円周率である。フィルタが次式のように表わされる理由については後述する。Ｇ（ｎ，ｍ，ω）は、時間領域変換部４に送られる。

時間領域変換部４は、フィルタを適用した後の信号Ｇ（ｎ，ｍ，ω）を逆フーリエ変換により時間領域信号ｇ（ｎ，ｍ，ｔ）に変換する（ステップＳ４）。逆フーリエ変換は既存の方法を用いればよい。時間領域信号ｇ（ｎ，ｍ，ｔ）は、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’に送られる。逆フーリエ変換によりフレーム毎に得られた時間領域信号ｇ（ｎ，ｍ，ｔ）は適宜シフトされて線形和が取られて、連続した時間領域信号となる。

スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’は、時間領域信号ｇ（ｎ，ｍ，ｔ）に基づいて音を再生する。より詳細には、ｎ＝１，…，Ｎ_ｘ’，ｍ＝１，…，Ｎ_ｙ’として、スピーカＳｎ−ｍが時間領域信号ｇ（ｎ，ｍ，ｔ）に基づいて音を再生する。これにより、第一の部屋のｚ＝ｚ_２の位置の波面Ｂを第二の部屋のスピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’で再現して、第一の部屋の音場を第二の部屋に再現することができる。

このように、マイクアレーよりも音源に近い位置での波面を不均質波を無視することにより模擬したフィルタを適用することにより、マイクアレーよりも音源Ｓ側にある波面Ｂを再現することができる。図２のように、第一の部屋において音源Ｓが位置ｚ＝ｚ_１とｚ＝ｚ_２との間にある場合には、第二の部屋においてスピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’の前面に音源Ｓ’があるような音場を再現することができる。

なお、再現できる波面の大きさは、マイクアレーで収音した平面の面積と同程度の大きさである。

ダイポール収音−モノポール再生
第一実施形態の音場収音再生装置において、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙをダイポール特性のマイクで構成し、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’をモノポール特性のスピーカで構成した場合には、フィルタ部３は以下の式により表わされる処理を行えばよい。

ダイポール特性のマイクｍ１は、図６に示すように例えば２個のモノポール特性のマイクｍ２により構成することができる。この場合、２個のモノポール特性のマイクｍ２で収音した音信号の差分を、ダイポール特性のマイクｍ１で収音した音信号とすればよい。図６において、２個のモノポール特性のマイクｍ２の中点がダイポール特性のマイクｍ１の位置となり、この位置にマイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙの面Ｍが形成される。２個のモノポール特性のマイクｍ２を結ぶ直線は、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙの面Ｍと直交する。

モノポール収音−ダイポール再生
第一実施形態の音場収音再生装置において、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙをモノポール特性のマイクで構成し、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’をダイポール特性のスピーカで構成した場合には、フィルタ部３は以下の式により表わされる処理を行えばよい。

ダイポール特性のスピーカは、図６に示したダイポール特性のマイクｍ１と同様に、例えば２個のモノポール特性のスピーカにより構成することができる。この場合、２個のモノポール特性のスピーカを逆相で駆動すれば、ダイポール特性を実現することができる。２個のモノポール特性のスピーカの中点にスピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’の面が位置し、かつ、２個のモノポール特性のスピーカを結ぶ直線は、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’の面と直交する。

ダイポール収音−ダイポール再生
第一実施形態の音場収音再生装置において、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−Ｎ_ｙをダイポール特性のマイクで構成し、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−Ｎ_ｙ’をダイポール特性のスピーカで構成した場合には、フィルタ部３は以下の式により表わされる処理を行えばよい。

［第二実施形態］
第一実施形態の音場収音再生装置及び方法は二次元マイクアレー及び二次元スピーカアレーを用いるのに対して、第二実施形態の音場収音再生装置及び方法は一次元マイクアレー及び二次元スピーカアレーを用いる。これにより、マイク数、スピーカ数及びチャネル数を少なくすることができるため、実装が比較的容易となる。

第二実施形態の音場収音再生装置及び方法は、図４に示す第一の部屋のｙ＝ｙ_１，ｚ＝ｚ_１の位置に配置されたＮ_ｘ個のモノポール特性のマイクで構成される一次元マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−１で収音した音信号を用いて以下に説明する方法でｙ＝ｙ_２，ｚ＝ｚ_２の位置の波面Ｂを計算して、第二の部屋に配置されたＮ_ｘ’個のモノポール特性のスピーカで構成される一次元スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１’で再現する。

Ｎ_ｘ，Ｎ_ｘ’は任意の整数である。Ｎ_ｘとＮ_ｘ’は異なっていてもよい。すなわち、マイクアレーを構成するマイクの数とスピーカアレーを構成するスピーカの数は同じでなくてもよく、マイクアレーを構成するマイクとスピーカアレーを構成するスピーカとは必ずしも一対一対応している必要はない。

第一の部屋のｙ＝ｙ_１，ｚ＝ｚ_１の位置に配置された一次元マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−１を構成する各マイクの位置をｒ_ｓ＝（ｘ_ｉ，ｙ_１，ｚ_１）と表わすことにする。第一の部屋のｙ＝ｙ_２，ｚ＝ｚ_２での波面Ｂを再現する第二の部屋に配置された一次元スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１を構成する各スピーカの位置をｒ_ｓ’＝（ｘ_ｎ，ｙ_２，ｚ_２）と表わすことにする。

第二実施形態の音場収音再生装置は、図３に示すように周波数領域変換部１、窓関数部２、フィルタ部３、時間領域変換部４、補正フィルタ部５を例えば含み、図５の実線及び破線で示された処理を行う。

第一の部屋のｙ＝ｙ_１，ｚ＝ｚ_１の位置に配置された一次元マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−１は、第一の部屋の音源Ｓで発せられた音を収音して時間領域の音信号を生成する。生成された音信号は、周波数領域変換部１に送られる。ｒ_ｓ＝（ｘ_ｉ，ｙ_１，ｚ_１）のマイクＭｉ−１で収音された時間領域の時刻ｔの音信号をｆ（ｉ，ｔ）と表記する。

周波数領域変換部１は、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−１で収音された音信号ｆ（ｉ，ｔ）をフーリエ変換により周波数領域信号Ｆ（ｉ，ω）に変換する（ステップＳ１）。生成された周波数領域信号Ｆ（ｉ，ω）は、窓関数部２に送られる。ωは周波数である。例えば、短時間フーリエ変換により周波数領域信号Ｆ（ｉ，ω）が生成される。もちろん、他の既存の方法により周波数領域信号Ｆ（ｉ，ω）を生成してもよい。

窓関数部２は、周波数領域信号Ｆ（ｉ，ω）に窓関数を乗じて窓関数後周波数領域信号Ｆ_ｗ（ｉ，ω）を生成する（ステップＳ２）。窓関数後周波数領域信号Ｆ_ｗ（ｉ，ω）は、フィルタ部３に送られる。窓関数として、以下の式より定義されるいわゆるターキー（Turkey）窓関数ｗ_ｘ（ｉ）を例えば用いる。Ｎ_ｔｐｒは、テーパーを適用する点数であり１以上Ｎ_ｘ以下の整数である。

フィルタ部３は、マイクアレーよりも音源に近い位置での波面Ｂを不均質波を無視することにより模擬したフィルタを窓関数後周波数領域信号Ｆ_ｗ（ｉ，ω）に適用する（ステップＳ３）。フィルタを適用した後の信号をＧ（ｎ，ω）と表記すると、フィルタ部３の処理は次式のように表わされる。ｋは波数であり、ｃを音速とすると、ｋ＝ω／ｃである。ｅはネピア数、πは円周率である。フィルタが次式のように表わされる理由については後述する。Ｇ（ｎ，ω）は、時間領域変換部４に送られる。

Ｈ_ｎ ^（１）（ｘ）は第一種ハンケル関数であり、Ｈ_ｎ ^（２）（ｘ）は第二種ハンケル関数である。第一種ハンケル関数Ｈ_ｎ ^（１）（ｘ）及び第二種ハンケル関数Ｈ_ｎ ^（２）（ｘ）は、第一種ベッセル関数Ｊ_ｎ（ｘ）及び第二種ベッセル関数Ｙ_ｎ（ｘ）を用いて以下のように定義される。

補正フィルタ部５は、１次アレーで近似することによる誤差を補正するために、以下の補正フィルタを用いてＧ（ｎ，ω）を補正して、補正後の信号Ｇ’（ｎ，ω）を得る（ステップＳ５）。

端数ｋ＝ω／ｃであり、ω＝０のＤＣ成分の場合には上記式の右辺の平方根の分母ｊｋが０となるため、上記の補正フィルタを適用することはできない。このため、ω＝０の場合には、Ｇ’（ｎ，ω）＝Ｇ’（ｎ，０）＝０とする。

時間領域変換部４は、補正フィルタを適用した後の信号Ｇ’（ｎ，ω）を逆フーリエ変換により時間領域信号ｇ（ｎ，ｔ）に変換する（ステップＳ４）。逆フーリエ変換は既存の方法を用いればよい。時間領域信号ｇ（ｎ，ｔ）は、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１に送られる。逆フーリエ変換によりフレーム毎に得られた時間領域信号ｇ（ｎ，ｔ）は適宜シフトされて線形和が取られて、連続した時間領域信号となる。

スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１は、時間領域信号ｇ（ｎ，ｔ）に基づいて音を再生する。より詳細には、ｎ＝１，…，Ｎ_ｘ’として、スピーカＳｎ−１が時間領域信号ｇ（ｎ，ｔ）に基づいて音を再生する。これにより、第一の部屋のｙ＝ｙ_２，ｚ＝ｚ_２の位置の波面Ｂを第二の部屋のスピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１で再現して、第一の部屋の音場を第二の部屋に再現することができる。

このように、マイクアレーよりも音源に近い位置での波面を不均質波を無視することにより模擬したフィルタを適用することにより、マイクアレーよりも音源Ｓ側にある波面Ｂを再現することができる。図４のように、第一の部屋において音源Ｓが位置ｚ＝ｚ_１とｚ＝ｚ_２との間にある場合には、第二の部屋においてスピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１の前面に音源Ｓ’があるような音場を再現することができる。

ダイポール収音−モノポール再生
第二実施形態の音場収音再生装置において、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−１をダイポール特性のマイクで構成し、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１をモノポール特性のスピーカで構成した場合には、フィルタ部３は以下の式により表わされる処理を行えばよい。

モノポール収音−ダイポール再生
第二実施形態の音場収音再生装置において、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−１をモノポール特性のマイクで構成し、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１をダイポール特性のスピーカで構成した場合には、フィルタ部３は以下の式により表わされる処理を行えばよい。

ダイポール収音−ダイポール再生
第一実施形態の音場収音再生装置において、マイクアレーＭ１−１，Ｍ２−１，…，ＭＮ_ｘ−１をダイポール特性のマイクで構成し、スピーカアレーＳ１−１，Ｓ２−１，…，ＳＮ_ｘ’−１をダイポール特性のスピーカで構成した場合には、フィルタ部３は以下の式により表わされる処理を行えばよい。

［理論的背景］
この発明の理論的背景を説明する。ここでは、第一実施形態においてダイポール収音−モノポール再生を行う場合にフィルタ部３の処理が式（２）となる理由を例に挙げて説明する。以下では、表記の簡略化のためにωを記載しない場合がある。

レイリー方程式に基づいて、図２においてｚ＝ｚ_１の位置ｒ_ｓの音圧勾配∂Ｐ（ｒ_ｓ）／∂ｚを用いてｚ＝ｚ_２の位置ｒ_ｓ’の音圧Ｐ（ｒ_ｓ’）を推定すると以下のようになる。

ここで、κ^−１は、核関数である。

ｋ_ｘ、ｋ_ｙは、それぞれｘ方向の波数、ｙ方向の波数である。核関数κ^−１は、均質波に対応する成分κ_１ ^−１と不均質波に対応する成分κ_２ ^−２との和に分解することが可能である。

ここで、ＷＦＳにおいて、不均質波の成分は音場の再生に寄与しない。したがって、核関数κ^−１は、均質波に対応する成分κ_１ ^−１のみで近似可能である。すなわち、以下の式が成立する。

このκ^−１を用いることにより、位置ｒ_ｓ’の音圧Ｐ（ｒ_ｓ’）は次のように表される。

再生するスピーカがモノポール特性である場合には、音圧をｚで偏微分した音圧勾配∂Ｐ（ｒ_ｓ’）／∂ｚを２倍した信号を再生することにより音場を再生することができる。

ダイポール特性のマイクで収音した場合には、Ｆ_ｗ（ｉ，ｊ，ω）＝∂Ｐ（ｒ_ｓ，ω）／∂ｚとなることを考慮して、上記式を離散的に表現したものが式（２）である。

ダイポールかモノポールか、アレーが二次元か一次元かの違いはあるが、他の場合についても同様にして式（１），式（３）から式（８）を導出することができるが、ここではそれらの導出を省略する。

［シミュレーション結果］
図７のように、ｚ＝０の位置に直線状に４ｃｍ間隔で配置されたマイクアレーでｚ＜０の領域の音場の音信号を収音して、ｚ＝０の位置に直線状に４ｃｍ間隔で配置されたスピーカアレーでｚ＞０の領域において収音された音信号を用いてその音場を再現するシミュレーションを行った。マイクアレー及びスピーカアレーの全体の長さは３．８４ｍである。音源は、（ｘ，ｙ，ｚ）＝（−０．４ｍ，０ｍ，−１．０ｍ）の位置にある点音源である。音源信号は、１ｋＨｚの正弦波である。

図８は従来のＷＦＳによってｚ＝０の波面を再現することにより音場を再現したシミュレーション結果であり、図９はこの場合の理想的な音場である。図１０はこの発明によりｚ＝−０．６ｍの波面を再現することにより音場を再現したシミュレーション結果であり、図１１はこの場合の理想的な音場である。図１２はこの発明によりｚ＝−２ｍの波面を再現することにより音場を再現したシミュレーション結果であり、図１３はこの場合の理想的な音場である。図８と図１２とを比較すると、図１２では波面がｚ軸の正方向に移動したように合成されていることがわかる。また、図８と図１３とを比較すると、スピーカアレーの前面に音源が形成されていることがわかる。

図１４は音源が（−０．４ｍ，０ｍ，−１．０ｍ）と（０．６ｍ，０ｍ，−１．２ｍ）の２つの位置にある場合にｚ＝−２．０ｍの波面を再現することにより音場を再現したシミュレーション結果であり、図１２この場合の理想的な音場である。図１４をみると、音源が複数ある場合にも音場を移動させることが可能であることがわかる。

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

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

窓関数の処理は、周波数領域で行っても、時間領域で行ってもよい。

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

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

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

１ 周波数領域変換部
２ 窓関数部
３ フィルタ部
４ 時間領域変換部
５ 補正フィルタ部

Claims (14)

1. マイクアレーで収音された音信号をフーリエ変換により周波数領域信号に変換する周波数領域変換部と、
   上記周波数領域信号に窓関数を乗じて窓関数後周波数領域信号を生成する窓関数部と、
   上記マイクアレーよりも音源に近い位置での波面を模擬したフィルタを上記窓関数後周波数領域信号に適用するフィルタ部と、
   上記フィルタを適用した後の信号を逆フーリエ変換により時間領域信号に変換する時間領域変換部と、
   を含む音場収音再生装置。
2. 請求項１に記載の音場収音再生装置において、
   上記音場収音再生装置は、収音装置及びと再生装置を備え、
   上記収音装置は、上記周波数領域変換部と上記窓関数部と上記フィルタ部とを含み、
   上記再生装置は、上記時間領域変換部を含む、
   ことを特徴とする音場収音再生装置。
3. 請求項１に記載の音場収音再生装置において、
   上記音場収音再生装置は、収音装置及びと再生装置を備え、
   上記収音装置は、上記周波数領域変換部を含み、
   上記再生装置は、上記窓関数部と上記フィルタ部と上記時間領域変換部とを含む、
   ことを特徴とする音場収音再生装置。
4. 請求項１に記載の音場収音再生装置において、
   上記音場収音再生装置は、収音装置及びと再生装置を備え、
   上記収音装置は、上記周波数領域変換部と上記窓関数部とを含み、
   上記再生装置は、上記フィルタ部と上記時間領域変換部を含む、
   ことを特徴とする音場収音再生装置。
5. 請求項１から４の何れかに記載の音場収音再生装置において、
   Ｎ_ｘ，Ｎ_ｙ，Ｎ_ｘ’，Ｎ_ｙ’を正の整数とし、上記マイクアレーは二次元平面に配置されたＮ_ｘ×Ｎ_ｙ個のモノポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは二次元平面に配置されたＮ_ｘ’×Ｎ_ｙ’個のモノポール特性のスピーカで構成されており、ｘ方向にｉ番目ｙ方向にｊ番目のマイクの位置を（ｘ_ｉ，ｙ_ｊ，ｚ_１）と表記し、ｘ方向にｎ番目ｙ方向のｍ番目のスピーカの位置を（ｘ_ｎ，ｙ_ｍ，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎｍ（ｉ，ｊ，ω）をフィルタ、Ｆ_ｗ（ｉ，ｊ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ｍ，ω）を上記フィルタを適用した後の信号として、
   上記フィルタ部は次式で定義される処理を行う、
   
   ことを特徴とする音場収音再生装置。
6. 請求項１から４の何れかに記載の音場収音再生装置において、
   Ｎ_ｘ，Ｎ_ｙ，Ｎ_ｘ’，Ｎ_ｙ’を正の整数とし、上記マイクアレーは二次元平面に配置されたＮ_ｘ×Ｎ_ｙ個のダイポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは二次元平面に配置されたＮ_ｘ’×Ｎ_ｙ’個のモノポール特性のスピーカで構成されており、ｘ方向にｉ番目ｙ方向にｊ番目のマイクの位置を（ｘ_ｉ，ｙ_ｊ，ｚ_１）と表記し、ｘ方向にｎ番目ｙ方向のｍ番目のスピーカの位置を（ｘ_ｎ，ｙ_ｍ，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎｍ（ｉ，ｊ，ω）をフィルタ、Ｆ_ｗ（ｉ，ｊ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ｍ，ω）を上記フィルタを適用した後の信号として、
   上記フィルタ部は次式で定義される処理を行う、
   
   ことを特徴とする音場収音再生装置。
7. 請求項１から４の何れかに記載の音場収音再生装置において、
   Ｎ_ｘ，Ｎ_ｙ，Ｎ_ｘ’，Ｎ_ｙ’を正の整数とし、上記マイクアレーは二次元平面に配置されたＮ_ｘ×Ｎ_ｙ個のモノポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは二次元平面に配置されたＮ_ｘ’×Ｎ_ｙ’個のダイポール特性のスピーカで構成されており、ｘ方向にｉ番目ｙ方向にｊ番目のマイクの位置を（ｘ_ｉ，ｙ_ｊ，ｚ_１）と表記し、ｘ方向にｎ番目ｙ方向のｍ番目のスピーカの位置を（ｘ_ｎ，ｙ_ｍ，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎｍ（ｉ，ｊ，ω）をフィルタ、Ｆ_ｗ（ｉ，ｊ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ｍ，ω）を上記フィルタを適用した後の信号として、
   上記フィルタ部は次式で定義される処理を行う、
   
   ことを特徴とする音場収音再生装置。
8. 請求項１から４の何れかに記載の音場収音再生装置において、
   Ｎ_ｘ，Ｎ_ｙ，Ｎ_ｘ’，Ｎ_ｙ’を正の整数とし、上記マイクアレーは二次元平面に配置されたＮ_ｘ×Ｎ_ｙ個のダイポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは二次元平面に配置されたＮ_ｘ’×Ｎ_ｙ’個のダイポール特性のスピーカで構成されており、ｘ方向にｉ番目ｙ方向にｊ番目のマイクの位置を（ｘ_ｉ，ｙ_ｊ，ｚ_１）と表記し、ｘ方向にｎ番目ｙ方向のｍ番目のスピーカの位置を（ｘ_ｎ，ｙ_ｍ，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎｍ（ｉ，ｊ，ω）をフィルタ、Ｆ_ｗ（ｉ，ｊ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ｍ，ω）を上記フィルタを適用した後の信号として、
   上記フィルタ部は次式で定義される処理を行う、
   
   ことを特徴とする音場収音再生装置。
9. 請求項１から４の何れかに記載の音場収音再生装置において、
   Ｎ_ｘ，Ｎ_ｘ’を正の整数とし、上記マイクアレーは一次元平面に配置されたＮ_ｘ個のモノポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは一次元平面に配置されたＮ_ｘ’個のモノポール特性のスピーカで構成されており、ｘ方向にｉ番目のマイクの位置を（ｘ_ｉ，ｙ_１，ｚ_１）と表記し、ｘ方向にｎ番目のスピーカの位置を（ｘ_ｎ，ｙ_２，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎ（ｉ，ω）をフィルタ、Ｆ_ｗ（ｉ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ω）を上記フィルタを適用した後の信号、Ｈ_ｎ ^（１）（・）を第一種ハンケル関数、Ｈ_ｎ ^（２）（・）を第二種ハンケル関数として、
   上記フィルタ部は次式で定義される処理を行い、
   
   上記フィルタを適用した後の信号Ｇ（ｎ，ω）に対して次式により定義されるフィルタ処理を行う補正フィルタ部を更に含み、
   
   上記時間領域変換部は、上記フィルタを適用した後の信号に代えて、上記補正フィルタによる補正後の信号Ｇ’（ｎ，ω）を逆フーリエ変換により時間領域信号に変換する、
   ことを特徴とする音場収音再生装置。
10. 請求項１から４の何れかに記載の音場収音再生装置において、
    Ｎ_ｘ，Ｎ_ｘ’を正の整数とし、上記マイクアレーは一次元平面に配置されたＮ_ｘ個のダイポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは一次元平面に配置されたＮ_ｘ’個のモノポール特性のスピーカで構成されており、ｘ方向にｉ番目のマイクの位置を（ｘ_ｉ，ｙ_１，ｚ_１）と表記し、ｘ方向にｎ番目のスピーカの位置を（ｘ_ｎ，ｙ_２，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎ（ｉ，ω）をフィルタ、Ｆ_ｗ（ｉ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ω）を上記フィルタを適用した後の信号、Ｈ_ｎ ^（１）（・）を第一種ハンケル関数、Ｈ_ｎ ^（２）（・）を第二種ハンケル関数として、
    上記フィルタ部は次式で定義される処理を行い、
    
    上記フィルタを適用した後の信号Ｇ（ｎ，ω）に対して次式により定義されるフィルタ処理を行う補正フィルタ部を更に含み、
    
    上記時間領域変換部は、上記フィルタを適用した後の信号に代えて、上記補正フィルタによる補正後の信号Ｇ’（ｎ，ω）を逆フーリエ変換により時間領域信号に変換する、
    ことを特徴とする音場収音再生装置。
11. 請求項１から４の何れかに記載の音場収音再生装置において、
    Ｎ_ｘ，Ｎ_ｘ’を正の整数とし、上記マイクアレーは一次元平面に配置されたＮ_ｘ個のモノポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは一次元平面に配置されたＮ_ｘ’個のダイポール特性のスピーカで構成されており、ｘ方向にｉ番目のマイクの位置を（ｘ_ｉ，ｙ_１，ｚ_１）と表記し、ｘ方向にｎ番目のスピーカの位置を（ｘ_ｎ，ｙ_２，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎ（ｉ，ω）をフィルタ、Ｆ_ｗ（ｉ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ω）を上記フィルタを適用した後の信号、Ｈ_ｎ ^（１）（・）を第一種ハンケル関数、Ｈ_ｎ ^（２）（・）を第二種ハンケル関数として、
    上記フィルタ部は次式で定義される処理を行い、
    
    上記フィルタを適用した後の信号Ｇ（ｎ，ω）に対して次式により定義されるフィルタ処理を行う補正フィルタ部を更に含み、
    
    上記時間領域変換部は、上記フィルタを適用した後の信号に代えて、上記補正フィルタによる補正後の信号Ｇ’（ｎ，ω）を逆フーリエ変換により時間領域信号に変換する、
    ことを特徴とする音場収音再生装置。
12. 請求項１から４の何れかに記載の音場収音再生装置において、
    Ｎ_ｘ，Ｎ_ｘ’を正の整数とし、上記マイクアレーは一次元平面に配置されたＮ_ｘ個のダイポール特性のマイクで構成されており、上記時間領域信号を再生するスピーカアレーは一次元平面に配置されたＮ_ｘ’個のダイポール特性のスピーカで構成されており、ｘ方向にｉ番目のマイクの位置を（ｘ_ｉ，ｙ_１，ｚ_１）と表記し、ｘ方向にｎ番目のスピーカの位置を（ｘ_ｎ，ｙ_２，ｚ_２）と表記し、ｋを波数とし、Ｈ_ｎ（ｉ，ω）をフィルタ、Ｆ_ｗ（ｉ，ω）を窓関数後周波数領域信号、Ｇ（ｎ，ω）を上記フィルタを適用した後の信号、Ｈ_ｎ ^（１）（・）を第一種ハンケル関数、Ｈ_ｎ ^（２）（・）を第二種ハンケル関数として、
    上記フィルタ部は次式で定義される処理を行い、
    
    上記フィルタを適用した後の信号Ｇ（ｎ，ω）に対して次式により定義されるフィルタ処理を行う補正フィルタ部を更に含み、
    
    上記時間領域変換部は、上記フィルタを適用した後の信号に代えて、上記補正フィルタによる補正後の信号Ｇ’（ｎ，ω）を逆フーリエ変換により時間領域信号に変換する、
    ことを特徴とする音場収音再生装置。
13. マイクアレーで収音された音信号をフーリエ変換により周波数領域信号に変換する周波数領域変換ステップと、
    上記周波数領域信号に窓関数を乗じて窓関数後周波数領域信号を生成する窓関数ステップと、
    上記マイクアレーよりも音源に近い位置での波面を模擬したフィルタを上記窓関数後周波数領域信号に適用するフィルタステップと、
    上記フィルタを適用した後の信号を逆フーリエ変換により時間領域信号に変換する時間領域変換ステップと、
    を含む音場収音再生方法。
14. 請求項１に記載された音場収音再生装置の各部としてコンピュータを機能させるための音場収音再生プログラム。