(57)【特許請求の範囲】
【請求項１】
音源の方向を示す方位情報を取得する方位情報取得部と、
前記方位情報に基づいて、空間フィルタの特性として中心周波数およびバンド幅を決定し、複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記空間フィルタを適用する空間フィルタ適用部と
を備える音声処理装置。
【請求項２】
前記空間フィルタは、前記中心周波数および前記バンド幅により定まる空間周波数の帯域を透過周波数帯域として、前記収音信号の透過周波数帯域の成分を透過させるフィルタである
請求項１に記載の音声処理装置。
【請求項３】
前記空間フィルタは、前記中心周波数および前記バンド幅により定まる時間周波数の帯域を透過周波数帯域として、前記収音信号の透過周波数帯域の成分を透過させるフィルタである
請求項１または請求項２に記載の音声処理装置。
【請求項４】
前記空間フィルタ適用部は、前記方位情報により示される前記音源の方向と、前記マイクアレイとのなす角度がπ/2に近くなるほど前記バンド幅が広くなるように前記空間フィルタの特性を決定する
請求項１乃至請求項３の何れか一項に記載の音声処理装置。
【請求項５】
前記マイクアレイは直線マイクアレイである
請求項１乃至請求項４の何れか一項に記載の音声処理装置。
【請求項６】
音源の方向を示す方位情報を取得し、
前記方位情報に基づいて、空間フィルタの特性として中心周波数およびバンド幅を決定し、複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記空間フィルタを適用する
ステップを含む音声処理方法。
【請求項７】
音源の方向を示す方位情報を取得し、
前記方位情報に基づいて、空間フィルタの特性として中心周波数およびバンド幅を決定し、複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記空間フィルタを適用する
ステップを含む処理をコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は音声処理装置および方法、並びにプログラムに関し、特に、より低コストで音像の定位を改善することができるようにした音声処理装置および方法、並びにプログラムに関する。
【背景技術】
【０００２】
従来、平面スピーカアレイや直線スピーカアレイを用いて音場を再現する波面合成技術が知られている。このような波面合成技術は、例えば図１に示すように次世代双方向通信などに利用することができる。【０００３】
図１では、通話者W11がいる空間P11と、通話者W12がいる空間P12との間で次世代双方向通信が行われる。【０００４】
具体的には、空間P11では、主に通話者W11が発した音声からなる音場Aが図中、縦方向に並べられた複数のマイクロフォンからなる直線マイクアレイMCA11により収音され、その結果得られた音源信号が空間P12へと伝送される。【０００５】
この例では、図中の矢印は音源である通話者W11からの音声の伝搬方向を示しており、通話者W11の音声は直線マイクアレイMCA11から見て角度θの方向から到来し、収音される。以下では、この角度θ、つまり音源からの音声の伝搬方向と、マイクアレイを構成するマイクロフォンが並ぶ方向とのなす角度を到来角度θと称することとする。【０００６】
空間P12では、空間P11から伝送されてきた音源信号から、音場Aを再現するためのスピーカ駆動信号が生成される。そして、空間P12内に図中、縦方向に並べられた複数のスピーカからなる直線スピーカアレイSPA11により、生成されたスピーカ駆動信号に基づいて音場Aが再現される。【０００７】
この例では、図中の矢印は直線スピーカアレイSPA11から出力され、通話者W12へと伝搬される音声の伝搬方向を示している。この伝搬方向と直線スピーカアレイSPA11とのなす角度は、到来角度θと同じ角度となる。【０００８】
なお、ここでは図示が省略されているが、空間P12にも直線マイクアレイが設けられており、この直線マイクアレイにより主に通話者W12が発した音声からなる音場Bが収音され、得られた音源信号が空間P11へと伝送される。また、空間P11では、空間P12から伝送されてきた音源信号からスピーカ駆動信号が生成され、得られたスピーカ駆動信号に基づいて図示せぬ直線スピーカアレイにより音場Bが再現される。【０００９】
ところで、このようにマイクアレイやスピーカアレイを利用して音場を再現する場合、物理的に正確な音場再現を行うためにはスピーカおよびマイクロフォンを無限個数配置しなければならない。例えば、図１に示した例のようにスピーカやマイクロフォンを離散的に配置すると、空間エリアシングが生じてしまう。【００１０】
空間エリアシングに侵されない最も高い空間周波数（以下、上限空間周波数と称する）klimは、スピーカアレイを構成するスピーカの間隔、またはマイクアレイを構成するマイクロフォンの間隔で計算される空間ナイキスト周波数のより低い方によって決定される。【００１１】
すなわち、マイクロフォンの間隔をdmicとし、スピーカの間隔をdspkとすると、上限空間周波数klimは次式（１）により求まる。【００１２】
【数１】
0006604331.tif 000002
【００１３】
このようにして得られる上限空間周波数klimは音像の定位に影響を与え、一般的に値が高い方が好ましい。【００１４】
また、音源信号の周波数（以下、時間周波数と称する）fと空間周波数kの関係は次式（２）に示す通りである。なお、式（２）において、cは音速を示している。【００１５】
【数２】
0006604331.tif 000003
【００１６】
したがって、特に工夫を行わない場合には、空間エリアシングに侵されない最も高い時間周波数（以下、上限時間周波数と称する）flimを、式（２）によって求めることが可能となる。上限時間周波数flimは音質に影響を与え、一般的には値が高い方が再現性が高い、つまりHiFi（High Fidelity）とされる。【００１７】
ここで、空間エリアシングについて説明する。図２は、音源からの音声の平面波の到来角度の違いによる空間スペクトルを示したものであり、空間スペクトルは平面波の到来角度によってスペクトルピークの位置が変化することから角度スペクトルとも呼ばれている。なお、図２において、縦軸は時間周波数fを示しており、横軸は空間周波数kを示している。また、直線L11乃至直線L13は、それぞれスペクトルピークを示している。【００１８】
図中、左側には、空間サンプリングされる前、すなわちマイクアレイにより収音される前の元の平面波の到来角度θ=0のときの角度スペクトルの様子が示されている。この例では、直線L11に示されるように、スペクトルピークは空間周波数kの正方向に現れる。【００１９】
これに対して図中、右側には、離散配置されたマイクロフォンからなるマイクアレイにより到来角度θ=0の平面波が空間サンプリング、つまり収音されて得られた音源信号の角度スペクトルの様子が示されている。【００２０】
この例では、直線L12は直線L11に対応し、本来現れるべきスペクトルピークを示している。また、直線L13は空間エリアシングにより出現したスペクトルピークを示しており、この例では時間周波数fが上限時間周波数flimより高く、かつ空間周波数kが負である領域に顕著に空間エリアシングが生じている。【００２１】
空間エリアシングがなければ、本来、空間周波数kが負である領域にスペクトルピークが現れるのは平面波の到来角度θがπ/2≦θ≦πであるときのはずである。【００２２】
したがって、図中、右側に示す例では、収音された平面波（音場）を再現すると、空間エリアシングにより生じたスペクトルピークの影響で、様々な異なる角度からの平面波信号が混在しているように音像定位してしまう。【００２３】
そのため、収音により得られた音源信号から波面合成のスピーカ駆動信号を生成し、そのスピーカ駆動信号に基づいてスピーカアレイで平面波を再現しても、聴者に正しい平面波を感じさせることができない。なお、上述した次世代双方向通信の例等で通話者同士が近づいた場合には平面波ではなく球面波になるが、平面波の場合と同様のことがいえる。【００２４】
以上のように、空間エリアシングが生じると音像を正確に定位させることができなくなってしまう。そこで、空間エリアシングの低減に関して、互いにスピーカ間隔が異なる高音用スピーカユニットと低音用スピーカユニットとの二種類のスピーカアレイを用いることで、空間エリアシングに侵されない上限時間周波数flimをより高くする技術が提案されている（例えば、特許文献１参照）。この技術によれば、より高い時間周波数を持つ信号の正確な再現が可能となる。
【先行技術文献】
【特許文献】
【００２５】
【特許文献１】特許第４１２４１８２号公報
【発明の概要】
【発明が解決しようとする課題】
【００２６】
しかしながら、上述した技術では、音像の定位を改善することはできるが、上限時間周波数flimをより高くするために追加のスピーカアレイ、つまり二種類のスピーカアレイを用意しなければならず、コストがかかってしまう。音像の定位を改善するにあたっては、既存のスピーカアレイのスピーカ間隔を保ちつつ、上限時間周波数flimをより高くすることが望ましい。【００２７】
本技術は、このような状況に鑑みてなされたものであり、より低コストで音像の定位を改善することができるようにするものである。
【課題を解決するための手段】
【００２８】
本技術の一側面の音声処理装置は、音源の方向を示す方位情報を取得する方位情報取得部と、前記方位情報に基づいて、空間フィルタの特性として中心周波数およびバンド幅を決定し、複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記空間フィルタを適用する空間フィルタ適用部とを備える。
【００３０】
前記空間フィルタを、前記中心周波数および前記バンド幅により定まる空間周波数の帯域を透過周波数帯域として、前記収音信号の透過周波数帯域の成分を透過させるフィルタとすることができる。【００３１】
前記空間フィルタを、前記中心周波数および前記バンド幅により定まる時間周波数の帯域を透過周波数帯域として、前記収音信号の透過周波数帯域の成分を透過させるフィルタとすることができる。【００３２】
前記空間フィルタ適用部には、前記方位情報により示される前記音源の方向と、前記マイクアレイとのなす角度がπ/2に近くなるほど前記バンド幅が広くなるように前記空間フィルタの特性を決定させることができる。【００３３】
前記マイクアレイを直線マイクアレイとすることができる。【００３４】
本技術の一側面の音声処理方法またはプログラムは、音源の方向を示す方位情報を取得し、前記方位情報に基づいて、空間フィルタの特性として中心周波数およびバンド幅を決定し、複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記空間フィルタを適用するステップを含む。
【００３５】
本技術の一側面においては、音源の方向を示す方位情報が取得され、前記方位情報に基づいて、空間フィルタの特性として中心周波数およびバンド幅が決定され、複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記空間フィルタが適用される。
【発明の効果】
【００３６】
本技術の一側面によれば、より低コストで音像の定位を改善することができる。【００３７】
なお、ここに記載された効果は必ずしも限定されるものではなく、本開示中に記載された何れかの効果であってもよい。
【図面の簡単な説明】
【００３８】
【図１】次世代双方向通信について説明する図である。
【図２】空間エリアシングについて説明する図である。
【図３】空間エリアシング制御器の構成例を示す図である。
【図４】空間バンドパスフィルタの特性について説明する図である。
【図５】空間バンドパスフィルタの特性について説明する図である。
【図６】音場再現処理を説明するフローチャートである。
【図７】コンピュータの構成例を示す図である。
【発明を実施するための形態】
【００３９】
以下、図面を参照して、本技術を適用した実施の形態について説明する。
【００４０】
〈第１の実施の形態〉
〈空間エリアシング制御器の構成例〉
本技術は、平面スピーカアレイや直線スピーカアレイを用いて音場を再現する場合に、波面合成を行うためのスピーカ駆動信号を生成する際、スピーカ駆動信号に適切な空間フィルタを適用することで、スピーカの離散配置によって生じる空間エリアシングを低減させるものである。
【００４１】
具体的には、本技術では、目的外方向に伝搬する波面を犠牲にして空間エリアシングを低減させることで上限時間周波数flimの上昇が実現される。
【００４２】
例えば、図１を参照して説明した次世代双方向通信のように通話者同士の位置関係が明らかな場合には、再現が優先される波面の伝搬方向とその他の目的外方向を特定することができる。そこで、特定された目的外方向の空間周波数を遮断することで、上限時間周波数flimを上昇させることができる。
【００４３】
次に、本技術を空間エリアシング制御器に適用した場合を例として、本技術を適用した具体的な実施の形態について説明する。
【００４４】
図３は、本技術を適用した空間エリアシング制御器の一実施の形態の構成例を示す図である。
【００４５】
空間エリアシング制御器１１は、送信器２１および受信器２２を有している。この例では、例えば送信器２１が音場を収音する収音空間に配置され、受信器２２が収音空間で収音された音場を再現する再現空間に配置される。
【００４６】
送信器２１は、音場を収音し、収音により得られた収音信号から空間周波数スペクトルを生成して受信器２２に送信する。受信器２２は、送信器２１から送信されてきた空間周波数スペクトルを受信してスピーカ駆動信号を生成するとともに、得られたスピーカ駆動信号に基づいて音場を再現する。
【００４７】
送信器２１は、マイクアレイ３１、時間周波数分析部３２、空間周波数分析部３３、および通信部３４を有している。また、受信器２２は、方位情報取得部３５、通信部３６、駆動信号生成部３７、空間フィルタ適用部３８、空間周波数合成部３９、時間周波数合成部４０、およびスピーカアレイ４１を有している。
【００４８】
マイクアレイ３１は、例えば直線状に配置された複数のマイクロフォンからなる直線マイクアレイなどからなり、到来した音声の平面波を収音し、その結果各マイクロフォンで得られた収音信号を時間周波数分析部３２に供給する。
【００４９】
時間周波数分析部３２は、マイクアレイ３１から供給された収音信号に対して時間周波数変換を行い、その結果得られた時間周波数スペクトルを空間周波数分析部３３に供給する。空間周波数分析部３３は、時間周波数分析部３２から供給された時間周波数スペクトルに対して空間周波数変換を行い、その結果得られた空間周波数スペクトルを通信部３４に供給する。
【００５０】
通信部３４は空間周波数分析部３３から供給された空間周波数スペクトルを、有線または無線により受信器２２の通信部３６に送信する。
【００５１】
また、受信器２２の方位情報取得部３５は、マイクアレイ３１により収音された音声の音源である通話者の方位（方向）を示す話者方位情報を取得して空間フィルタ適用部３８に供給する。
【００５２】
ここでは、収音される音場の音源が通話者である例について説明するが、音源は通話者に限らず、車等のオブジェクトや環境音の音源など、どのようなものであってもよい。また、話者方位情報は、音源である通話者からの音声を聴取する受話者に対する通話者の方向など、主たる音源と受話者の相対的な位置関係を示すものであればどのようなものであってもよいが、ここでは話者方位情報は、上述した到来角度θであるものとして説明を続ける。この場合、例えば図１に示した次世代双方向通信の例において、通話者W11が音源とされ、通話者W12が受話者とされるときには、通話者W11から通話者W12へと向かう音声の波面の伝搬方向を示す到来角度θが話者方位情報とされる。また、以下では、到来角度θを示す話者方位情報を話者方位情報θとも記すこととする。
【００５３】
通信部３６は、通信部３４から送信されてきた空間周波数スペクトルを受信して駆動信号生成部３７に供給する。駆動信号生成部３７は、通信部３６から供給された空間周波数スペクトルに基づいて、収音された音場を再現するための空間領域のスピーカ駆動信号を生成し、空間フィルタ適用部３８に供給する。
【００５４】
空間フィルタ適用部３８は、駆動信号生成部３７から供給されたスピーカ駆動信号に対して、方位情報取得部３５から供給された話者方位情報により定まる特性の空間フィルタを用いたフィルタ処理を施し、その結果得られた空間フィルタスペクトルを空間周波数合成部３９に供給する。
【００５５】
空間周波数合成部３９は、空間フィルタ適用部３８から供給された空間フィルタスペクトルの空間周波数合成を行い、その結果得られた時間周波数スペクトルを時間周波数合成部４０に供給する。
【００５６】
時間周波数合成部４０は、空間周波数合成部３９から供給された時間周波数スペクトルの時間周波数合成を行い、その結果得られたスピーカ駆動信号をスピーカアレイ４１に供給する。スピーカアレイ４１は、例えば直線状に配置された複数のスピーカからなる直線スピーカアレイなどからなり、時間周波数合成部４０から供給されたスピーカ駆動信号に基づいて音声を再生する。これにより、収音空間における音場が再現される。
【００５７】
ここで、空間エリアシング制御器１１を構成する各部についてより詳細に説明する。
【００５８】
（時間周波数分析部）
時間周波数分析部３２は、マイクアレイ３１を構成する各マイクロフォンで得られた収音信号s(nmic,t)の時間周波数情報を分析する。
【００５９】
但し、収音信号s(nmic,t)においてnmicはマイクアレイ３１を構成するマイクロフォンを示すマイクインデクスであり、マイクインデクスnmic＝0,…,Nmic-1である。ここで、Nmicはマイクアレイ３１を構成するマイクロフォンの数である。また、収音信号s(nmic,t)においてtは時間を示している。
【００６０】
時間周波数分析部３２は、収音信号s(nmic,t)に対して固定サイズの時間フレーム分割を行って入力フレーム信号sfr(nmic,nfr,l)を得る。そして、時間周波数分析部３２は、次式（３）に示す窓関数wT(nfr)を入力フレーム信号sfr(nmic,nfr,l)に乗算し、窓関数適用信号sw(nmic,nfr,l)を得る。すなわち、以下の式（４）の計算が行われて窓関数適用信号sw(nmic,nfr,l)が算出される。
【００６１】
【数３】
0006604331.tif 000004
【００６２】
【数４】
0006604331.tif 000005
【００６３】
ここで、式（３）および式（４）において、nfrは時間フレーム内のサンプルを示す時間インデクスを示しており、時間インデクスnfr＝0,…,Nfr-1である。また、lは時間フレームインデクスを示しており、時間フレームインデクスl＝0,…,L-1である。なお、Nfrはフレームサイズ（時間フレームのサンプル数）であり、Lは総フレーム数である。
【００６４】
また、フレームサイズNfrは、時間サンプリング周波数fsT[Hz]における一フレームの時間Tfr[s]相当のサンプル数Nfr（＝R(fsT×Tfr)、但しR()は任意の丸め関数）である。この実施の形態では、例えば一フレームの時間Tfr＝1.0[s]であり、丸め関数R()は四捨五入であるが、それ以外でも構わない。さらに、フレームのシフト量はフレームサイズNfrの50％としているが、それ以外でも構わない。
【００６５】
さらに、ここでは窓関数としてハニング窓の平方根を用いているが、ハミング窓やブラックマンハリス窓などのその他の窓を用いるようにしてもよい。
【００６６】
このようにして窓関数適用信号sw(nmic,nfr,l)が得られると、時間周波数分析部３２は、以下の式（５）および式（６）を計算することで、窓関数適用信号sw(nmic,nfr,l)に対して時間周波数変換を行い、時間周波数スペクトルS(nmic,nT,l)を算出する。
【００６７】
【数５】
0006604331.tif 000006
【００６８】
【数６】
0006604331.tif 000007
【００６９】
すなわち、式（５）の計算によりゼロ詰め信号sw’(nmic,mT,l)が求められ、得られたゼロ詰め信号sw’(nmic,mT,l)に基づいて式（６）が計算され、時間周波数スペクトルS(nmic,nT,l)が算出される。
【００７０】
なお、式（５）および式（６）において、MTは時間周波数変換に用いるポイント数を示している。また、nTは時間周波数スペクトルインデクスを示している。ここで、nT＝0,…,NT-1であり、NT＝MT/2+1である。さらに、式（６）においてiは純虚数を示している。
【００７１】
また、この実施の形態では、STFT（Short Time Fourier Transform）（短時間フーリエ変換）による時間周波数変換を行っているが、DCT（Discrete Cosine Transform）（離散コサイン変換）やMDCT（Modified Discrete Cosine Transform）（修正離散コサイン変換）などの他の時間周波数変換を用いてもよい。
【００７２】
さらに、STFTのポイント数MTは、Nfr以上である、Nfrに最も近い2のべき乗の値としているが、それ以外のポイント数MTでも構わない。
【００７３】
時間周波数分析部３２は、以上において説明した処理で得られた時間周波数スペクトルS(nmic,nT,l)を、空間周波数分析部３３に供給する。
【００７４】
（空間周波数分析部）
続いて空間周波数分析部３３は、時間周波数分析部３２から供給された時間周波数スペクトルS(nmic,nT,l)に対して、次式（７）を計算することで空間周波数変換を行い、空間周波数スペクトルSSP(nS,nT,l)を算出する。
【００７５】
【数７】
0006604331.tif 000008
【００７６】
なお、式（７）においてMSは空間周波数変換に用いるポイント数を示しており、mS＝0,…,MS-1である。また、S’(mS,nT,l)は時間周波数スペクトルS(nmic,nT,l)に対してゼロ詰めを行うことにより得られるゼロ詰め時間周波数スペクトルを示しており、iは純虚数を示している。さらに、nSは空間周波数スペクトルインデクスを示している。
【００７７】
この実施の形態では、式（７）の計算によってIDFT（Inverse Discrete Fourier Transform）（逆離散フーリエ変換）による空間周波数変換が行われている。
【００７８】
また、必要であればIDFTのポイント数MSに合わせて適切にゼロ詰めを行ってもよい。この例では、0≦mS≦Nmic-1であるポイントmSについては、ゼロ詰め時間周波数スペクトルS’(mS,nT,l)＝時間周波数スペクトルS(nmic,nT,l)とされ、Nmic≦mS≦MS-1であるポイントmSについては、ゼロ詰め時間周波数スペクトルS’(mS,nT,l)＝0とされる。
【００７９】
以上において説明した処理により得られる空間周波数スペクトルSSP(nS,nT,l)は、時間フレームlに含まれている時間周波数nTの信号が空間上においてどのような波形となっているかを示している。空間周波数分析部３３は、空間周波数スペクトルSSP(nS,nT,l)を通信部３４に供給する。
【００８０】
（駆動信号生成部）
駆動信号生成部３７には、通信部３６および通信部３４を介して、空間周波数分析部３３から空間周波数スペクトルSSP(nS,nT,l)が供給される。
【００８１】
駆動信号生成部３７は、空間周波数スペクトルSSP(nS,nT,l)に基づいて、以下の式（８）を計算し、スピーカアレイ４１で音場（波面）を再現するための空間領域のスピーカ駆動信号DSP(mS,nT,l)を求める。すなわち、SDM(Spectral Division Method)により、空間周波数スペクトルであるスピーカ駆動信号DSP(mS,nT,l)が算出される。
【００８２】
【数８】
0006604331.tif 000009
【００８３】
なお、式（８）においてyrefはSDMの基準距離を示しており、基準距離yrefは波面が正確に再現される位置である。この基準距離yrefはマイクアレイ３１を構成するマイクロフォンが並ぶ方向と垂直な方向の距離となる。例えば、ここでは基準距離yref＝1[m]とされているが、その他の値でもよい。
【００８４】
さらに、式（８）においてH0(2)は第二種ハンケル関数を示しており、K0はベッセル関数を示している。また、式（８）においてiは純虚数を示しており、cは音速を示しており、ωは時間角周波数を示している。
【００８５】
さらに、式（８）においてkは空間周波数を示しており、mS、nT、lは、それぞれ空間周波数スペクトルインデクス、時間周波数スペクトルインデクス、時間フレームインデクスを示している。
【００８６】
なお、ここではSDMによりスピーカ駆動信号DSP(mS,nT,l)を算出する手法を例として説明したが、他の手法によりスピーカ駆動信号が算出されてもよい。また、SDMについては、特に「Jens Adrens, Sascha Spors, “Applying the Ambisonics Approach on Planar and Linear Arrays of Loudspeakers”, in 2nd International Symposium on Ambisonics and Spherical Acoustics」に詳細に記載されている。
【００８７】
駆動信号生成部３７は、以上のようにして得られたスピーカ駆動信号DSP(mS,nT,l)を空間フィルタ適用部３８に供給する。
【００８８】
（空間フィルタ適用部）
空間フィルタ適用部３８は、駆動信号生成部３７から供給されたスピーカ駆動信号DSP(mS,nT,l)と、方位情報取得部３５から供給された話者方位情報θにより特性が決定される空間バンドパスフィルタBθ(mS,nT)とを用いて、空間フィルタスペクトルF(mS,nT,l)を求める。なお、ここでは空間バンドパスフィルタBθ(mS,nT)の形状が矩形状であるものとするが、空間バンドパスフィルタBθ(mS,nT)の形状は他のどのような形状であってもよい。
【００８９】
具体的には、空間フィルタ適用部３８は話者方位情報θに基づいて、空間バンドパスフィルタBθ(mS,nT)の中心周波数kcenおよびバンド幅klenを決定することで、空間バンドパスフィルタBθ(mS,nT)の特性を決定する。すなわち、目的とする主たる音源からの音声の平面波の到来角度θに応じて、空間バンドパスフィルタBθ(mS,nT)の特性が決定される。
【００９０】
例えば、空間フィルタ適用部３８は、以下の式（９）を計算することで中心周波数kcenを算出し、以下の式（１０）を計算することでバンド幅klenを算出する。
【００９１】
【数９】
0006604331.tif 000010
【００９２】
【数１０】
0006604331.tif 000011
【００９３】
なお、式（９）および式（１０）において、θは話者方位情報、すなわち音源から出力され、受話者へと向かう平面波（音声）のマイクアレイ３１への到来角度を示している。また、式（９）および式（１０）において、klimはマイクアレイ３１のマイクロフォン間隔と、スピーカアレイ４１のスピーカ間隔とから定まる上限空間周波数を示している。
【００９４】
空間バンドパスフィルタBθ(mS,nT)は、中心周波数kcenを中心とするバンド幅klenの空間周波数帯域を透過周波数帯域（パスバンド）とし、それ以外の空間周波数帯域を遮断周波数帯域（ストップバンド）とするバンドパスフィルタである。
【００９５】
したがって、空間バンドパスフィルタBθ(mS,nT)の値は、空間周波数スペクトルインデクスmSにより示される空間周波数が透過周波数帯域内の周波数であれば1となり、空間周波数スペクトルインデクスmSにより示される空間周波数が遮断周波数帯域内の周波数であれば0となる。これにより、空間バンドパスフィルタBθ(mS,nT)は、透過周波数帯域の成分のみを透過させる空間フィルタとなる。
【００９６】
このようにして空間バンドパスフィルタが決定されると、空間フィルタ適用部３８は、次式（１１）に示すように空間バンドパスフィルタBθ(mS,nT)をスピーカ駆動信号DSP(mS,nT,l)に乗算し、空間フィルタスペクトルF(mS,nT,l)を得る。
【００９７】
【数１１】
0006604331.tif 000012
【００９８】
空間フィルタ適用部３８は、式（１１）の計算により得られた空間フィルタスペクトルF(mS,nT,l)を、空間周波数合成部３９に供給する。
【００９９】
ここで、空間バンドパスフィルタBθ(mS,nT)について、さらに説明する。
【０１００】
式（９）により示される中心周波数kcenは、話者方位情報θ、つまり到来角度θがθ＝0であるときにはklim/2となり、到来角度θがπ/2に近づくにしたがって小さくなり、到来角度θ＝π/2のときには0となる。また、中心周波数kcenは、到来角度θがθ＝π/2からθ＝πに近づくにしたがって小さくなり、到来角度θ＝πのときには-klim/2となる。
【０１０１】
さらに式（１０）により示されるバンド幅klenは、到来角度θがθ＝0であるときにはklimとなり、到来角度θがπ/2に近づくにしたがって大きくなり、到来角度θ＝π/2のときには2klimとなる。また、バンド幅klenは、到来角度θがθ＝π/2からθ＝πに近づくにしたがって小さくなり、到来角度θ＝πのときにはklimとなる。
【０１０２】
したがって、例えば到来角度θ＝0であるときには、図４に示すように空間周波数kが0から上限空間周波数klimまでの間の帯域が透過周波数帯域となる。なお、図４において縦軸は時間周波数fを示しており、横軸は空間周波数kを示している。
【０１０３】
到来角度θ＝0である場合、矢印A11に示すようにマイクアレイ３１で収音された平面波の空間スペクトル（角度スペクトル）では、直線L21に示されるスペクトルピークと、直線L22に示されるスペクトルピークとが観測される。
【０１０４】
ここで、空間周波数kがk≧0である領域に現れている直線L21に示されるスペクトルピークは、本来現れるべきスペクトルピークである。これに対して、直線L22に示されるスペクトルピークは、空間エリアシングにより現れたものであり、空間周波数kが負である領域で空間エリアシングが顕著であることが分かる。
【０１０５】
この例では、空間エリアシングに起因する、直線L22に示されるスペクトルピークのない時間周波数fの領域が非エリアシング帯域R11となっている。これに対して、非エリアシング帯域R11の上限の時間周波数、つまり上述した上限時間周波数flimより時間周波数が高い領域が、空間エリアシングに侵されたエリアシング帯域R12となっている。
【０１０６】
また、到来角度θ＝0である場合、空間バンドパスフィルタBθ(mS,nT)の特性は、上述した式（９）および式（１０）から、折れ線L23に示される特性となる。
【０１０７】
すなわち、矢印A12に示すように空間バンドパスフィルタBθ(mS,nT)は、中心周波数kcen＝klim/2であり、バンド幅klen＝klimであるバンドパスフィルタとなる。したがって、この例では空間周波数kが開始周波数sb＝0から終了周波数eb＝klimまでの間の成分がそのまま透過とされ、それ以外の遮断周波数帯域の成分は空間バンドパスフィルタBθ(mS,nT)により遮断（除去）される。
【０１０８】
図４では、図中の斜線が施された領域が空間バンドパスフィルタBθ(mS,nT)により遮断される領域を示しており、この例では直線L22に示される空間エリアシングのスペクトルピークのうちの空間周波数kが負である領域の部分が除去されることになる。
【０１０９】
その結果、空間エリアシングのスペクトルピークのない時間周波数fの領域である非エリアシング帯域R13は非エリアシング帯域R11よりも広くなり、その分だけ空間エリアシングに侵されたエリアシング帯域R14が狭くなる。換言すれば、空間バンドパスフィルタBθ(mS,nT)によるフィルタ処理によって、上限時間周波数flimをより高くすることができる。この例では、負の空間周波数成分を低減させることにより、空間エリアシングに侵されない上限時間周波数flimが二倍に引き上げられている。
【０１１０】
このように空間バンドパスフィルタBθ(mS,nT)によるフィルタ処理を行なえば、上限時間周波数flimを上昇させることができるので、特に到来角度θ、つまり角度θで伝搬する平面波の音質を向上させることができる。また、空間エリアシングを低減させることができるので、本来ないはずの他の角度から伝搬してきた平面波が混在しているような音像の定位を改善することができる。つまり、より正確な音像定位を実現することができる。
【０１１１】
同様に式（９）および式（１０）から、例えば到来角度θ＝π/2であるときには図５に示すように、空間周波数kが-klimから上限空間周波数klimまでの間の帯域が透過周波数帯域となる。なお、図５において縦軸は時間周波数fを示しており、横軸は空間周波数kを示している。
【０１１２】
到来角度θ＝π/2である場合、矢印A21に示すようにマイクアレイ３１で収音された平面波の空間スペクトル（角度スペクトル）では、直線L31に示されるスペクトルピークが観測される。
【０１１３】
ここで、空間周波数kがk＝0である領域に現れている直線L31に示されるスペクトルピークは、本来現れるべきスペクトルピークである。また、到来角度θ＝π/2である場合には空間エリアシングは生じないため、空間エリアシングに起因するスペクトルピークは出現しない。この例では、空間エリアシングは発生しないため、時間周波数fの全領域が非エリアシング帯域となっている。
【０１１４】
このように到来角度θ＝π/2である場合、空間バンドパスフィルタBθ(mS,nT)の特性は、上述した式（９）および式（１０）から、折れ線L32に示される特性となる。
【０１１５】
すなわち、矢印A22に示すように空間バンドパスフィルタBθ(mS,nT)は、中心周波数kcen＝0であり、バンド幅klen＝2klimであるバンドパスフィルタとなる。したがって、この例では空間周波数kが開始周波数sb＝-klimから終了周波数eb＝klimまでの間の成分がそのまま透過とされ、それ以外の遮断周波数帯域の成分は空間バンドパスフィルタBθ(mS,nT)により遮断される。
【０１１６】
なお、図５では、図中の斜線が施された領域が空間バンドパスフィルタBθ(mS,nT)により遮断される領域を示している。この例では、開始周波数sbが-klimであり、終了周波数ebがklimであるため、特に正および負の空間周波数成分の低減は行われない。
【０１１７】
以上のように空間バンドパスフィルタBθ(mS,nT)によるフィルタ処理を行なえば、上限時間周波数flimを上昇させ、特に目的とする到来角度θの方向に伝搬する平面波について音質を向上させるとともに、音像の定位を改善することができる。
【０１１８】
但し、空間バンドパスフィルタBθ(mS,nT)により遮断周波数帯域の成分を除去すると、角度θとは異なる他の角度で伝搬する平面波の音質は、除去された成分に応じて劣化してしまう。そのため、その分だけ再現空間において良好な音質で音声を聴取できる領域の範囲が狭くなってしまう。
【０１１９】
しかし、空間エリアシング制御器１１では、到来角度θがπ/2に近くなるにしたがって、つまり空間エリアシングが少なくなるにしたがってバンド幅klenがより広くなり、良好な音質で音声を聴取できる領域が広くなるようになされているので、フィルタ処理により生じる影響を小さく抑えることができる。
【０１２０】
なお、以上においては空間バンドパスフィルタBθ(mS,nT)の特性として、空間周波数について話者方位情報θに応じた透過周波数帯域を設定する例について説明したが、時間周波数についてのみ、話者方位情報θに応じた透過周波数帯域を設定するようにしてもよい。
【０１２１】
また、空間バンドパスフィルタBθ(mS,nT)の特性として、空間周波数と時間周波数の両方について、それぞれ話者方位情報θに応じた透過周波数帯域を設定するようにしてもよい。そのような場合、空間周波数だけでなく、時間周波数についても話者方位情報θに応じた中心周波数およびバンド幅、つまり透過周波数帯域が決定される。そして、空間周波数スペクトルインデクスmSにより示される空間周波数が透過周波数帯域内の周波数であり、かつ時間周波数スペクトルインデクスnTにより示される時間周波数が透過周波数帯域内の周波数であるときに、空間バンドパスフィルタBθ(mS,nT)の値は1となる。つまり、空間バンドパスフィルタBθ(mS,nT)は、空間周波数の透過周波数帯域、かつ時間周波数の透過周波数帯域の成分のみを透過させる空間フィルタとなる。
【０１２２】
（空間周波数合成部）
続いて、空間周波数合成部３９について説明する。
【０１２３】
空間周波数合成部３９は、次式（１２）を計算することで、空間フィルタ適用部３８から供給された空間フィルタスペクトルF(mS,nT,l)の空間周波数合成、すなわち空間フィルタスペクトルF(mS,nT,l)に対する逆空間周波数変換を行い、時間周波数スペクトルD(nspk,nT,l)を算出する。式（１２）では、逆空間周波数変換として、DFT（Discrete Fourier Transform）（離散フーリエ変換）が行われる。
【０１２４】
【数１２】
0006604331.tif 000013
【０１２５】
なお、式（１２）において、nspkはスピーカアレイ４１を構成するスピーカを特定するスピーカインデクスを示している。また、MSはDFTのポイント数を示しており、iは純虚数を示している。
【０１２６】
空間周波数合成部３９は、このようにして得られた時間周波数スペクトルD(nspk,nT,l)を時間周波数合成部４０に供給する。
【０１２７】
（時間周波数合成部）
時間周波数合成部４０は、次式（１３）の計算を行うことで、空間周波数合成部３９から供給された時間周波数スペクトルD(nspk,nT,l)の時間周波数合成を行い、出力フレーム信号dfr(nspk,nfr,l)を得る。ここでは、時間周波数合成として、ISTFT（Inverse Short Time Fourier Transform）（短時間逆フーリエ変換）が用いられているが、時間周波数分析部３２で行われる時間周波数変換（順変換）の逆変換に相当するものを用いればよい。
【０１２８】
【数１３】
0006604331.tif 000014
【０１２９】
なお、式（１３）におけるD’(nspk,mT,l)は、次式（１４）により得られるものである。
【０１３０】
【数１４】
0006604331.tif 000015
【０１３１】
式（１３）においてiは純虚数を示しており、nfrは時間インデクスを示している。また、式（１３）および式（１４）において、MTはISTFTのポイント数を示しており、nspkはスピーカインデクスを示している。
【０１３２】
さらに、時間周波数合成部４０は、得られた出力フレーム信号dfr(nspk,nfr,l)に、窓関数wT(nfr)を乗算し、オーバーラップ加算を行うことでフレーム合成を行う。例えば、次式（１５）の計算によりフレーム合成が行われて、出力信号d(nspk,t)が求められる。
【０１３３】
【数１５】
0006604331.tif 000016
【０１３４】
なお、ここでは、出力フレーム信号dfr(nspk,nfr,l)に乗算する窓関数wT(nfr)として、時間周波数分析部３２で用いた窓関数と同じものを用いているが、ハミング窓などのその他の窓の場合は矩形窓で構わない。
【０１３５】
また、式（１５）において、dprev(nspk,nfr+lNfr)およびdcurr(nspk,nfr+lNfr)は、どちらも出力信号d(nspk,t)を示しているが、dprev(nspk,nfr+lNfr)は更新前の値を示し、dcurr(nspk,nfr+lNfr)は更新後の値を示している。
【０１３６】
時間周波数合成部４０は、このようにして得られた出力信号d(nspk,t)を、スピーカ駆動信号としてスピーカアレイ４１に供給する。
【０１３７】
〈音場再現処理の説明〉
次に、以上において説明した空間エリアシング制御器１１により行われる処理の流れについて説明する。空間エリアシング制御器１１は、収音空間における音声の平面波の収音が指示されると、その平面波の収音を行って音場を再現する音場再現処理を行う。
【０１３８】
以下、図６のフローチャートを参照して空間エリアシング制御器１１による音場再現処理について説明する。
【０１３９】
ステップＳ１１において、マイクアレイ３１は収音空間において音声の平面波を収音し、その結果得られた収音信号s(nmic,t)を時間周波数分析部３２に供給する。
【０１４０】
ステップＳ１２において、時間周波数分析部３２はマイクアレイ３１から供給された収音信号s(nmic,t)の時間周波数情報を分析する。
【０１４１】
具体的には、時間周波数分析部３２は収音信号s(nmic,t)に対して時間フレーム分割を行い、その結果得られた入力フレーム信号sfr(nmic,nfr,l)に窓関数wT(nfr)を乗算し、窓関数適用信号sw(nmic,nfr,l)を算出する。
【０１４２】
また、時間周波数分析部３２は、窓関数適用信号sw(nmic,nfr,l)に対して時間周波数変換を行い、その結果得られた時間周波数スペクトルS(nmic,nT,l)を空間周波数分析部３３に供給する。すなわち、式（６）の計算が行われて時間周波数スペクトルS(nmic,nT,l)が算出される。
【０１４３】
ステップＳ１３において、空間周波数分析部３３は、時間周波数分析部３２から供給された時間周波数スペクトルS(nmic,nT,l)に対して空間周波数変換を行い、その結果得られた空間周波数スペクトルSSP(nS,nT,l)を通信部３４に供給する。
【０１４４】
具体的には、空間周波数分析部３３は式（７）を計算することで、時間周波数スペクトルS(nmic,nT,l)を空間周波数スペクトルSSP(nS,nT,l)に変換する。
【０１４５】
ステップＳ１４において、通信部３４は、空間周波数分析部３３から供給された空間周波数スペクトルSSP(nS,nT,l)を、無線通信により再現空間に配置された受信器２２に送信する。そして、ステップＳ１５において、受信器２２の通信部３６は、無線通信により送信されてきた空間周波数スペクトルSSP(nS,nT,l)を受信し、駆動信号生成部３７に供給する。
【０１４６】
ステップＳ１６において、方位情報取得部３５は話者方位情報θを取得し、空間フィルタ適用部３８に供給する。例えば話者方位情報θは、予め定められていてもよいし、送信器２１等から取得されるようにしてもよい。
【０１４７】
ステップＳ１７において、駆動信号生成部３７は、通信部３６から供給された空間周波数スペクトルSSP(nS,nT,l)に基づいて空間領域のスピーカ駆動信号DSP(mS,nT,l)を算出し、空間フィルタ適用部３８に供給する。例えば駆動信号生成部３７は、式（８）を計算することで、空間領域のスピーカ駆動信号DSP(mS,nT,l)を算出する。
【０１４８】
ステップＳ１８において、空間フィルタ適用部３８は、方位情報取得部３５から供給された話者方位情報θに基づいて、空間バンドパスフィルタBθ(mS,nT)の特性を決定する。
【０１４９】
例えば、空間フィルタ適用部３８は、上述した式（９）および式（１０）を計算して、空間バンドパスフィルタBθ(mS,nT)の中心周波数kcenおよびバンド幅klenを算出することで空間バンドパスフィルタBθ(mS,nT)の特性、つまり透過周波数帯域を決定する。
【０１５０】
ステップＳ１９において、空間フィルタ適用部３８は、駆動信号生成部３７から供給されたスピーカ駆動信号DSP(mS,nT,l)に対して、決定された特性の空間バンドパスフィルタBθ(mS,nT)を適用する。
【０１５１】
すなわち、空間フィルタ適用部３８は式（１１）の計算を行うことで、スピーカ駆動信号DSP(mS,nT,l)に対して、空間バンドパスフィルタBθ(mS,nT)を用いたフィルタ処理を施し、空間フィルタスペクトルF(mS,nT,l)を得る。
【０１５２】
空間フィルタ適用部３８は、フィルタ処理により得られた空間フィルタスペクトルF(mS,nT,l)を、空間周波数合成部３９に供給する。
【０１５３】
ステップＳ２０において、空間周波数合成部３９は、空間フィルタ適用部３８から供給された空間フィルタスペクトルF(mS,nT,l)に対して逆空間周波数変換を行い、その結果得られた時間周波数スペクトルD(nspk,nT,l)を時間周波数合成部４０に供給する。例えば空間周波数合成部３９は、式（１２）を計算することで逆空間周波数変換を行う。
【０１５４】
ステップＳ２１において、時間周波数合成部４０は、空間周波数合成部３９から供給された時間周波数スペクトルD(nspk,nT,l)の時間周波数合成を行う。
【０１５５】
具体的には、時間周波数合成部４０は式（１３）の計算を行って、時間周波数スペクトルD(nspk,nT,l)から出力フレーム信号dfr(nspk,nfr,l)を算出する。さらに時間周波数合成部４０は、出力フレーム信号dfr(nspk,nfr,l)に窓関数wT(nfr)を乗算して式（１５）の計算を行い、フレーム合成により出力信号d(nspk,t)を算出する。
【０１５６】
時間周波数合成部４０は、このようにして得られた出力信号d(nspk,t)を、スピーカ駆動信号としてスピーカアレイ４１に供給する。
【０１５７】
ステップＳ２２において、スピーカアレイ４１は、時間周波数合成部４０から供給されたスピーカ駆動信号に基づいて音声を再生し、音場再現処理は終了する。このようにしてスピーカ駆動信号に基づいて音声が再生されると、再現空間において収音空間の音場が再現される。
【０１５８】
以上のようにして空間エリアシング制御器１１は、話者方位情報θに基づいて空間バンドパスフィルタBθ(mS,nT)の特性を決定するとともに、空間バンドパスフィルタBθ(mS,nT)をスピーカ駆動信号DSP(mS,nT,l)に適用し、空間エリアシングを低減させる。
【０１５９】
このように話者方位情報θに応じた特性の空間バンドパスフィルタBθ(mS,nT)を用いて空間エリアシングを低減させることで上限時間周波数flimを上昇させ、音質を向上させるとともに音像の定位を改善することができる。しかも、空間エリアシング制御器１１では、特別なスピーカアレイを用意する必要もなく、フィルタ処理という簡単な処理で空間エリアシングを低減させることができるので、より低コストで上限時間周波数を上昇させることができる。
【０１６０】
ところで、上述した一連の処理は、ハードウェアにより実行することもできるし、ソフトウェアにより実行することもできる。一連の処理をソフトウェアにより実行する場合には、そのソフトウェアを構成するプログラムが、コンピュータにインストールされる。ここで、コンピュータには、専用のハードウェアに組み込まれているコンピュータや、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のパーソナルコンピュータなどが含まれる。
【０１６１】
図７は、上述した一連の処理をプログラムにより実行するコンピュータのハードウェアの構成例を示すブロック図である。
【０１６２】
コンピュータにおいて、CPU（Central Processing Unit）５０１，ROM（Read Only Memory）５０２，RAM（Random Access Memory）５０３は、バス５０４により相互に接続されている。
【０１６３】
バス５０４には、さらに、入出力インターフェース５０５が接続されている。入出力インターフェース５０５には、入力部５０６、出力部５０７、記録部５０８、通信部５０９、及びドライブ５１０が接続されている。
【０１６４】
入力部５０６は、キーボード、マウス、マイクロフォン、撮像素子などよりなる。出力部５０７は、ディスプレイ、スピーカなどよりなる。記録部５０８は、ハードディスクや不揮発性のメモリなどよりなる。通信部５０９は、ネットワークインターフェースなどよりなる。ドライブ５１０は、磁気ディスク、光ディスク、光磁気ディスク、又は半導体メモリなどのリムーバブルメディア５１１を駆動する。
【０１６５】
以上のように構成されるコンピュータでは、CPU５０１が、例えば、記録部５０８に記録されているプログラムを、入出力インターフェース５０５及びバス５０４を介して、RAM５０３にロードして実行することにより、上述した一連の処理が行われる。
【０１６６】
コンピュータ（CPU５０１）が実行するプログラムは、例えば、パッケージメディア等としてのリムーバブルメディア５１１に記録して提供することができる。また、プログラムは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の伝送媒体を介して提供することができる。
【０１６７】
コンピュータでは、プログラムは、リムーバブルメディア５１１をドライブ５１０に装着することにより、入出力インターフェース５０５を介して、記録部５０８にインストールすることができる。また、プログラムは、有線または無線の伝送媒体を介して、通信部５０９で受信し、記録部５０８にインストールすることができる。その他、プログラムは、ROM５０２や記録部５０８に、あらかじめインストールしておくことができる。
【０１６８】
なお、コンピュータが実行するプログラムは、本明細書で説明する順序に沿って時系列に処理が行われるプログラムであっても良いし、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで処理が行われるプログラムであっても良い。
【０１６９】
また、本技術の実施の形態は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。
【０１７０】
例えば、本技術は、１つの機能をネットワークを介して複数の装置で分担、共同して処理するクラウドコンピューティングの構成をとることができる。
【０１７１】
また、上述のフローチャートで説明した各ステップは、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０１７２】
さらに、１つのステップに複数の処理が含まれる場合には、その１つのステップに含まれる複数の処理は、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０１７３】
さらに、本技術は、以下の構成とすることも可能である。
【０１７４】
［１］
音源の方向を示す方位情報を取得する方位情報取得部と、
複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記方位情報により定まる特性の空間フィルタを適用する空間フィルタ適用部と
を備える音声処理装置。
［２］
前記空間フィルタ適用部は、前記方位情報に基づいて、前記空間フィルタの特性として中心周波数およびバンド幅を決定する
［１］に記載の音声処理装置。
［３］
前記空間フィルタは、前記中心周波数および前記バンド幅により定まる空間周波数の帯域を透過周波数帯域として、前記収音信号の透過周波数帯域の成分を透過させるフィルタである
［２］に記載の音声処理装置。
［４］
前記空間フィルタは、前記中心周波数および前記バンド幅により定まる時間周波数の帯域を透過周波数帯域として、前記収音信号の透過周波数帯域の成分を透過させるフィルタである
［２］または［３］に記載の音声処理装置。
［５］
前記空間フィルタ適用部は、前記方位情報により示される前記音源の方向と、前記マイクアレイとのなす角度がπ/2に近くなるほど前記バンド幅が広くなるように前記空間フィルタの特性を決定する
［２］乃至［４］の何れか一項に記載の音声処理装置。
［６］
前記マイクアレイは直線マイクアレイである
［１］乃至［５］の何れか一項に記載の音声処理装置。
［７］
音源の方向を示す方位情報を取得し、
複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記方位情報により定まる特性の空間フィルタを適用する
ステップを含む音声処理方法。
［８］
音源の方向を示す方位情報を取得し、
複数のマイクロフォンからなるマイクアレイにより前記音源からの音声を収音して得られた収音信号に対して、前記方位情報により定まる特性の空間フィルタを適用する
ステップを含む処理をコンピュータに実行させるプログラム。
【符号の説明】
【０１７５】
１１  空間エリアシング制御器，  ３１  マイクアレイ，  ３２  時間周波数分析部，  ３３  空間周波数分析部，  ３５  方位情報取得部，  ３７  駆動信号生成部，  ３８  空間フィルタ適用部，  ３９  空間周波数合成部，  ４０  時間周波数合成部，  ４１  スピーカアレイ
【図１】
0006604331.tif 000017
【図２】
0006604331.tif 000018
【図３】
0006604331.tif 000019
【図４】
0006604331.tif 000020
【図５】
0006604331.tif 000021
【図６】
0006604331.tif 000022
【図７】
0006604331.tif 000023