JP2019054344A

JP2019054344A - フィルタ係数算出装置、収音装置、その方法、及びプログラム

Info

Publication number: JP2019054344A
Application number: JP2017175898A
Authority: JP
Inventors: 江村　暁; Akira Emura; 暁江村
Original assignee: Nippon Telegraph and Telephone Corp
Current assignee: Nippon Telegraph and Telephone Corp
Priority date: 2017-09-13
Filing date: 2017-09-13
Publication date: 2019-04-04
Anticipated expiration: 2037-09-13
Also published as: JP6815956B2

Abstract

【課題】マイクロホン信号からターゲット音をよりクリアに抽出するためのフィルタ係数ベクトルを求める。【解決手段】フィルタ係数算出装置は、Nを2以上の整数の何れかとし、N個のマイクロホンからなるマイクロホンアレーにおける音源方向に対するステアリングベクトルをもとに、ターゲット音を無歪みで出力するという拘束条件のもとで、マイクロホンアレーの周波数領域のNチャネルマイクロホン信号にフィルタ係数ベクトルを適用して得られる値が疎になるようにフィルタ係数ベクトルを求める。【選択図】図３

Description

本発明は、複数のマイクロホンを用いてビームを形成するビームフォーミング技術を用いた収音装置、収音装置において用いるフィルタ係数を算出するフィルタ係数算出装置、その方法、及びプログラムに関する。

複数のマイクロホンを音場に設置してマルチチャネルのマイクロホン信号を取得し、そこからターゲットとする音声や音(以下、ターゲット音ともいう)をクリアに、ノイズやその他音声(以下、非ターゲット音ともいう)をできる限り取り除いて、取り出す技術のニーズが近年高まっている。そのために、複数のマイクロホンをもちいてビームを形成するビームフォーミング技術が近年さかんに研究開発されている。

ビームフォーミング技術では、図１のようにN個のマイクロホン９１−ｎ(ただし、n=1,2,…,N)で収音した各マイクロホン信号y_n(t)にフィルタリング部９２−ｎにおいてフィルタを適用する。なお、tは時刻を示すインデックスである。次に、加算部９３においてフィルタリング部９２−ｎの出力値の総和をとる。求めた総和を収音装置の出力信号z(t)として出力する。このような構成により雑音を大幅に減らし、ターゲット音をより明瞭に取り出すことができる。このようなビームフォーミングのフィルタを求める方法として、minimum variance distortionless response法（MVDR法）がよく使われる（非特許文献１参照）。

図２を用いて、MVDR法を説明する。なお、以下では、周波数fにおいて目的音源に対する各マイクロホンの応答a(f)(ステアリングベクトル)が既知であると想定している。非特許文献２では、マイクロホン信号からa(f)を相対的に推定する方法が示されており、この想定は妥当である。

マイクロホンアレイからのNチャネルマイクロホン信号y_n(t)(1≦n≦N)はフレームごとに短時間フーリエ変換部１０７において短時間フーリエ変換される。その周波数f、フレームlでの変換結果を

のようにベクトル化して扱う。このNチャネルマイクロホン信号y(f,l)は、
y(f,l)=x(f,l)+v(f,l) (2)
のようにターゲット音の直接波のマルチチャネル信号x(f,l）と、反射残響成分および雑音のマルチチャネル信号v(f,l）からなる。

相関行列算出部１０１では、Nチャネルマイクロホン信号y(f,l)の周波数fでの空間相関行列を
R(f,l)=E[y(f,l)y^H(f,l)] (3)
で算出する。ただしE[ ]は期待値をとることを意味する。またy^H(f,l)は、y(f,l)を転置し複素共役をとったベクトルである。実際の処理では通常E[ ]の代わりに短時間平均を用いる。

ステアリングベクトル決定部１０２では、非特許文献２等の方法により、ステアリングベクトルa(f)を決める。非特許文献２で推定したa(f)を用いる場合、参照マイクロホンが収音する目的音源の音がターゲット音となる。

アレーフィルタ推定部１０４では、次の拘束条件つき最適化問題を解いて、フィルタ係数ベクトル(すなわちN次元複素数ベクトル)h(f,l)を求める。
h(f,l)=arg min h^H(f,l)R(f,l)h(f,l) (4)
拘束条件
h^H(f,l)a(f)=1
上記の最適化問題は、周波数fにおいてターゲット音を無歪みで出力するという拘束条件のもとで、フィルタ係数ベクトルh(f,l)をNチャネルマイクロホン信号y(f,l)に適用した際に得られる値のパワーが最小になるようにフィルタ係数ベクトルh(f,l)を求めている。

アレーフィルタリング部１０５では、次式により、推定されたフィルタ係数ベクトルh(f,l)をマイクロホン信号y(f,l)に適用する。
z(f,l)=h^H(f,l)y(f,l) (5)
これにより、ターゲット音以外の成分を極力抑えてターゲット音を取り出すことができる。全周波数での処理結果を短時間逆フーリエ変換部１０８において短時間逆フーリエ変換することで、ターゲット音を取り出すことができる。

D. H. Johnson, D. E. Dudgeon, "Array Signal Processing", Prentice Hall, 1993. S. Araki, H. Sawada, and S. Makino," Blind speech separation in a meeting situation with maximum SNR beamformer", in proc. IEEE Int. Conf. Acoust. Speech Signal Process. (ICASSP2007), 2007, pp. 41-44.

上記のMVDR法によりターゲット音がクリアに抽出されるためには、相関行列が高精度に推定されている必要がある。しかしながら、実際に処理に使われる相関行列は、有限のサンプルから短時間平均をもちいて算出されている。そのために相関行列の各要素の値は真値と完全に一致するわけではなく、一定のばらつき（分散）をもってしまう。統計学上、サンプル数を増加させることで、この分散を小さくすることは可能だが0にはできない。

真の相関行列からフィルタ係数ベクトルを求めることができれば、理想的な収音ビームを形成できる。しかしながら、実際には推定誤差のある相関行列からマイクロホンアレイのフィルタ係数ベクトルが求められる。このフィルタ係数が形成する収音ビームは理想の収音ビームからズレが生じ、雑音や残響の抑え込みが本来より弱くなる。そのため、理想的な収音ビームのようにはターゲット音をクリアに抽出できない。この場合、スペクトルグラムでフィルタの出力信号をみると、ターゲット音成分のないところでも、残響音の成分や雑音の成分が重畳してしまう。

そこで本発明では、マイクロホン信号からターゲット音を従来よりクリアに抽出するためのフィルタ係数ベクトルを求めるフィルタ係数算出装置、求めたフィルタ係数ベクトルを用いる収音装置、その方法、及びプログラムを提供することを目的とする。

上記の課題を解決するために、本発明の一態様によれば、フィルタ係数算出装置は、Nを2以上の整数の何れかとし、N個のマイクロホンからなるマイクロホンアレーにおける音源方向に対するステアリングベクトルをもとに、ターゲット音を無歪みで出力するという拘束条件のもとで、マイクロホンアレーの周波数領域のNチャネルマイクロホン信号にフィルタ係数ベクトルを適用して得られる値が疎になるようにフィルタ係数ベクトルを求める。

上記の課題を解決するために、本発明の他の態様によれば、フィルタ係数算出方法は、Nを2以上の整数の何れかとし、N個のマイクロホンからなるマイクロホンアレーにおける音源方向に対するステアリングベクトルをもとに、ターゲット音を無歪みで出力するという拘束条件のもとで、マイクロホンアレーの周波数領域のNチャネルマイクロホン信号にフィルタ係数ベクトルを適用して得られる値が疎になるようにフィルタ係数ベクトルを求める。

本発明によれば、マイクロホン信号からターゲット音を従来よりクリアに抽出することができるという効果を奏する。

従来技術に係る収音装置の機能ブロック図。従来技術に係るMVDR法を説明するための図。第一、第二及び第三実施形態に係る収音装置の機能ブロック図。第一、第二及び第三実施形態に係る収音装置の処理フローの例を示す図。第二実施形態に係るステアリングベクトル決定部の機能ブロック図。第二実施形態に係るステアリングベクトル決定部の処理フローの例を示す図。

以下、本発明の実施形態について、説明する。なお、以下の説明に用いる図面では、同じ機能を持つ構成部や同じ処理を行うステップには同一の符号を記し、重複説明を省略する。以下の説明において、ベクトルや行列の各要素単位で行われる処理は、特に断りが無い限り、そのベクトルやその行列の全ての要素に対して適用されるものとする。

＜第一実施形態＞
図３は第一実施形態に係る収音装置の機能ブロック図を、図４はその処理フローを示す。

以下、N個のマイクロホン９１−ｎ(n=1,2,…,N)からなるマイクロホンアレイの出力信号(以下、Nチャネルマイクロホン信号ともいう)y_n(t)を対象とし、そのNチャネルマイクロホン信号y_n(t)から、ターゲット音をとりだす収音装置について説明する。この例では、マイクロホン９１−ｎは無指向性のマイクロホン素子からなる。

収音装置２００は、Nチャネルマイクロホン信号y_n(t)を入力とし、ターゲット音を取り出し、出力信号z(t)とし、出力する。

収音装置２００は、短時間フーリエ変換部１０７と、ステアリングベクトル決定部１０２と、アレーフィルタ推定部２０４と、アレーフィルタリング部１０５と、短時間逆フーリエ変換部１０８とを含む。

収音装置２００は、例えば、中央演算処理装置（CPU: Central Processing Unit）、主記憶装置（RAM: Random Access Memory）などを有する公知又は専用のコンピュータに特別なプログラムが読み込まれて構成された特別な装置である。収音装置２００は、例えば、中央演算処理装置の制御のもとで各処理を実行する。収音装置２００に入力されたデータや各処理で得られたデータは、例えば、主記憶装置に格納され、主記憶装置に格納されたデータは必要に応じて中央演算処理装置へ読み出されて他の処理に利用される。収音装置２００の各処理部は、少なくとも一部が集積回路等のハードウェアによって構成されていてもよい。収音装置２００が備える各記憶部は、例えば、RAM（Random Access Memory）などの主記憶装置、ハードディスクや光ディスクもしくはフラッシュメモリ（Flash Memory）のような半導体メモリ素子により構成される補助記憶装置、またはリレーショナルデータベースやキーバリューストアなどのミドルウェアにより構成することができる。

＜短時間フーリエ変換部１０７＞
短時間フーリエ変換部１０７は、Nチャネルの時間領域のマイクロホン信号y_n(t)を入力とし、フレームl(エル)毎に周波数領域のマイクロホン信号Y_n(f,l)に短時間フーリエ変換し(Ｓ１０７)、出力する。その周波数f、フレームlでの変換結果を

のようにベクトル化して扱う。

＜ステアリングベクトル決定部１０２＞
ステアリングベクトル決定部１０２は、ステアリングベクトルa(f)を求め（Ｓ１０２）、出力する。ステアリングベクトルの決定方法には様々な公知の技術を利用できる。例えば、ステアリングベクトル決定部１０２は、Nチャネルの周波数領域のマイクロホン信号y(f,l)を入力とし、非特許文献２の方法により、ステアリングベクトルa(f)を求める。

＜アレーフィルタ推定部２０４＞
アレーフィルタ推定部２０４は、Nチャネルの周波数領域のマイクロホン信号y(f,l)とステアリングベクトルa(f)とを入力とし、ステアリングベクトルa(f)をもとに、以下の拘束条件付き最適化問題を解いて、フィルタ係数ベクトルh(f,l)を求め（Ｓ２０４）、出力する。

拘束条件
h^H(f,l)a(f)=1
ただし、|・|₁は、ベクトル各成分の絶対値の総和すなわちＬ１ノルムを取ることを意味する。従来の最適化では、2乗ノルムすなわちL2ノルムを用いるコスト関数が良く使われてきた。本実施形態では、L2ノルムの代わりにL1ノルムを用いてコスト関数を最適化することで、スパースなベクトル、すなわち、0を多く含むベクトルが得られる。L1ノルムを用いてコスト関数を最適化する方法は、近年、圧縮センシングの分野で知られている（参考文献１参照）。
（参考文献１）田中利幸、「圧縮センシングの数理」、IEICE Fundamental Review, vol. 4, no. 1, pp. 39-47, 2010

上記の最適化問題では、周波数fごとに、ターゲット音（参照マイクロホンが収音する目的音源の音）を無歪みで出力するという拘束条件のもとで、Lフレーム分のNチャネルの周波数領域のマイクロホン信号y(f,l-L+1),y(f,l-L+2),…,y(f,l)にフィルタ係数ベクトルh(f,l)を適用して得られる値(例えば、h^H(f,l)y(f,l-L+1),h^H(f,l)y(f,l-L+2),…,h^H(f,l)y(f,l))が疎になるようにフィルタ係数ベクトルh(f,l)を求めている。スペクトルグラム上で信号成分のスパース化をかけることで、雑音成分や残響成分を抑えることが可能となる。なお、フレーム数Lは、1以上の整数の何れかとし、例えば、100前後の値に設定することがのぞましい。また、フィルタ係数ベクトルh(f,l)は、適切なタイミングで推定し、出力すればよく、例えば、Lフレームの分のNチャネルの周波数領域のマイクロホン信号y(f,l)を蓄積しておき0.5〜1秒ごとにフィルタ係数ベクトルh(f,l)を求める構成としてもよい。

＜アレーフィルタリング部１０５＞
アレーフィルタリング部１０５は、Nチャネルの周波数領域のマイクロホン信号y(f,l)とフィルタ係数ベクトルh(f,l)とを入力とし、次式のようにNチャネルマイクロホン信号y(f,l)にフィルタ係数ベクトルh(f,l)を適用し、周波数領域の出力信号z(f,l)を求め（Ｓ１０５）、出力する。
z(f,l)=h^H(f,l)y(f,l) (5)

＜短時間逆フーリエ変換部１０８＞
短時間逆フーリエ変換部１０８は、周波数領域の出力信号z(f,l)を入力とし、全周波数での処理結果を短時間逆フーリエ変換し（Ｓ１０８）、時間領域の出力信号z(t)を得、出力する。

＜効果＞
以上の構成により、Nチャネルマイクロホン信号からターゲット音を従来よりクリアに抽出することができる。

＜変形例＞
本実施形態では、収音装置として説明したが、アレーフィルタ推定部のみからなるフィルタ係数算出装置であっても、マイクロホン信号からターゲット音を従来よりクリアに抽出するためのフィルタ係数を求めることができる。また、収音装置は、アレーフィルタ推定部とアレーフィルタリング部のみを含む構成としてもよい。別装置からNチャネルの周波数領域のマイクロホン信号y(f,l)とステアリングベクトルa(f)とを受け取り、周波数領域の出力信号z(f,l)を求め、別装置に出力する構成としてもよい。

なお、本実施形態では、N個の無指向性のマイクロホンからなるマイクロホンアレイのマイクロホン信号y_n(t)を入力としているが、必ずしも無指向性である必要はない。また、あるマイクロホン素子の指向性と他のマイクロホン素子の指向性とは、同一であってもよいし、異なってもよい。

＜第二実施形態＞
第一実施形態と異なる部分を中心に説明する。

図３は第二実施形態に係る収音装置の機能ブロック図を、図４はその処理フローを示す。

収音装置３００は、ステアリングベクトル決定部１０２に代えてステアリングベクトル決定部３０２を含む。

＜ステアリングベクトル決定部３０２＞
図５はステアリングベクトル決定部３０２の機能ブロック図を、図６はその処理フローの例を示す。

ステアリングベクトル決定部３０２は、Nチャネルの周波数領域のマイクロホン信号y(f,l)を入力とし、ステアリングベクトルa(f)を求め（Ｓ３０２）、出力する。本実施形態では、ステアリングベクトル決定部３０２は、ノイズ・到来波分解部３０２１及びターゲット方向判定部３０２２を含む。

＜ノイズ・到来波分解部３０２１＞
ノイズ・到来波分解部３０２１は、Nチャネルマイクロホン信号y(f,l)の空間相関行列R(f,l)から、周波数fにおける複数到来波の強度および各マイクロホンのノイズパワーを推定し（Ｓ３０２１）、出力する。

例えば、ノイズ・到来波分解部３０２１は、周波数領域のマイクロホン信号y(f,l)を入力とし、周波数f、フレームlにおけるマイクロホン信号y(f,l)を用いて、その空間相関行列R(f,l)を算出する。例えば、次式により算出する。
R(f,l)=E[y(f,l)y(f,l)^H] (22)

ただしE[ ]は期待値をとることを意味する。また、y(f,l)^Hは、y(f,l)を転置し複素共役をとったベクトルである。実際の処理では通常E[ ]の代わりに短時間平均を用いる。

そして、ノイズ・到来波分解部３０２１は、空間相関行列R(f,l)からK個の方向からの到来波の強度の推定値p_k(f,l)及び各マイクロホン信号Y_n(f,l)に含まれるノイズパワーの推定値q_n(f,l)を求め（Ｓ３０２１）、p_k(f,l)及びq_n(f,l)を対角成分とする対角行列V(f,l)を出力する。ただし、kを到来方向のインデックスとし、平面波の到来可能方向としてK方向を想定し、k=1,2,…,Kとする。よって、対角行列V(f,l)は、以下のように表される。

なおK>Nである。強度の推定値p_k(f,l)及びノイズパワーの推定値q_n(f,l)の推定方法として、例えば、参考文献２の方法をもちいることができる。
（参考文献２）P. Stoica, P. Babu, and J. Li, "SPICE A sparse covariance-based estimation method for array processing", IEEE Transactions on signal processing, vol. 59, no. 2, 2011, 629-638.

この方法では、あらかじめ平面波の到来可能方向としてK方向(>N)を想定する。周波数ｆにおいて、k番目の方向から振幅１の平面波がマイクロホンアレイに到達したとき、その各マイクロホンのレスポンス(出力信号)をa_k(f)=[a_k,1(f) a_k,2(f) … a_k,N(f)]^Tとする。a_k(f)をk番目の方向に対するステアリングベクトルともいう。a_k,n(f)は、周波数ｆにおいて、k番目の方向からの到来する振幅１の平面波に対するn番目のマイクロホンのレスポンス(出力信号)を表す。なお、a_k(f)は、収音に先立ち予め求めておく。だだし、a_k(f)は、実験(実測)やシミュレーションにより予め求めてもよいし、計算による理論値を用いてもよい。K個のレスポンスベクトルa_k(f)とN×N単位行列I_Nからなる行列
A(f)^H=[a₁(f) a₂(f) … a_K(f) I_N] (23)
を用いて、参考文献２では
R(f,l)=A(f)^HV(f,l)A(f) (24)
の形に行列R(f,l)を行列A(f)^H、対角行列V(f,l)と行列A(f)の積に分解する。この分解により、対角行列V(f,l)に含まれるk番目の方向からの平面波の強度の推定値p_k(f,l)と、n番目のマイクロホン９１−ｎのノイズパワーの推定値q_n(f,l)とが得られる。なお実際には、上記の分解は、
||(A(f)^HV(f,l)A(f))^-1/2(R(f,l)-A(f)^HV(f,l)A(f))R(f,l)^-1/2||² (25)
を最小にする対角行列V(f,l)を求めることに対応する。なお、この式(25)で||x||は行列xのフロベニウスノルムをとることを意味する。

＜ターゲット方向判定部３０２２＞
ターゲット方向判定部３０２２は、ターゲット音の到来方向の推定値k_tを求め（Ｓ３０２２）、推定値k_tに対応するステアリングベクトルa(f)=a_kt(f)を出力する。例えば、ターゲット方向判定部３０２２は、対角行列V(f,l)を入力とし、対角行列V(f,l)に含まれる各到来方向kの強度の推定値p_k(f,l)を用いて、強度が所定の値より大きい方向をターゲット音の到来方向と判定し、判定結果(到来方向の推定値)k_tを求める。この例では、ターゲット方向判定部３０２２は、音声パワーが集中している帯域100〜500Hzの強度の推定値p_k(f,l)を用いてターゲット音の到来方向の推定値k_tを求める。この帯域で各到来方向kの強度は

になる。この例では、f₀は100Hz、f₁は500Hzに相当する。b(k,l)がピークを取る方向k_peakを、フレームlでのターゲット音の到来方向(以下、ターゲット音方向ともいう)の候補とする。そして、b(k,l)の最大値をb_maxとし、ピーク値がb_max×α以上の値を持つピークの位置をターゲット音方向k_tとして抽出する。ターゲット音方向k_tは複数の場合もあるし1つの場合もある。αとしては、例えば-12dB〜-6dBの範囲の値をもちいればよい。ターゲット音の到来方向の推定値k_tに対応するステアリングベクトルa(f)=a_kt(f)、つまり、行列A(f)^H=[a₁(f) a₂(f) … a_K(f) I_N]の第k_t列のベクトルを出力する。

アレーフィルタ推定部２０４は、第一実施形態と同様の処理を行う。ただし、ターゲット音方向k_tが複数の場合には、各ターゲット音方向k_tに対してフィルタ係数ベクトルh_kt(f,l)を求めればよい。

＜効果＞
このような構成とすることで、第一実施形態と同様の効果を得ることができる。なお、本実施形態では、空間相関行列を算出し、利用しているが、空間相関行列はターゲット音の到来方向を推定するために利用しているのであって、フィルタ係数ベクトルを求める最適化問題に利用しているのではない。そのため、推定誤差のある相関行列からマイクロホンアレイのフィルタ係数ベクトルが求め、理想の収音ビームからズレが生じ、雑音や残響の抑え込みが本来より弱くなるという問題は生じない。

＜第三実施形態＞
第一実施形態と異なる部分を中心に説明する。

第一実施形態のアレーフィルタ推定部２０４は、ノルムの大きいフィルタ係数ベクトルh(f,l)を推定することがある。このようなフィルタ係数ベクトルh(f,l)は、ターゲット音やノイズの特性変動への感度が高くなり、現実の信号に対して雑音抑圧性能や非ターゲット音抑圧性能が劣化してしまう場合がある。

この劣化を防止するために、フィルタ係数ベクトルh(f,l)のノルムを一定量以下に設定する仕組みを導入する。

図３は第三実施形態に係る収音装置の機能ブロック図を、図４はその処理フローを示す。

収音装置４００は、アレーフィルタ推定部２０４に代えてアレーフィルタ推定部４０４を含む。

＜アレーフィルタ推定部４０４＞
アレーフィルタ推定部４０４は、Nチャネルの周波数領域のマイクロホン信号y(f,l)とステアリングベクトルa(f)とを入力とし、ステアリングベクトルa(f)をもとに、以下の拘束条件付き最適化問題を解いて、フィルタ係数ベクトルh(f,l)を求め（Ｓ４０４）、出力する。

拘束条件
h^H(f,l)a(f)=1

ただし||h(f,l)||はベクトルh(f,l)の２乗ノルムをとることを意味する。またλは事前に指定するパラメータであり、10前後の値が設定される。つまり、上記の最適化問題では、周波数fごとに、ターゲット音を無歪みで出力し、かつ、フィルタ係数ベクトルh(f,l)のノルムが所定の値以下であるというという拘束条件のもとで、Lフレーム分のNチャネルの周波数領域のマイクロホン信号y(f,l-L+1),y(f,l-L+2),…,y(f,l)にフィルタ係数ベクトルh(f,l)を適用して得られる値が疎になるようにフィルタ係数ベクトルh(f,l)を求めている。

＜効果＞
このような構成とすることで、第一実施形態と同様の効果を得ることができる。さらに、ノルムの大きいフィルタ係数ベクトルh(f,l)を推定することを防ぎ、雑音抑圧性能や非ターゲット音抑圧性能の劣化を抑えることができる。なお、本実施形態と第二実施形態を組合せてもよい。

＜その他の変形例＞
本発明は上記の実施形態及び変形例に限定されるものではない。例えば、上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。

＜プログラム及び記録媒体＞
また、上記の実施形態及び変形例で説明した各装置における各種の処理機能をコンピュータによって実現してもよい。その場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記各装置における各種の処理機能がコンピュータ上で実現される。

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。

また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させてもよい。

このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶部に格納する。そして、処理の実行時、このコンピュータは、自己の記憶部に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実施形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよい。さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、プログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。

また、コンピュータ上で所定のプログラムを実行させることにより、各装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。

Claims

Nを2以上の整数の何れかとし、N個のマイクロホンからなるマイクロホンアレーにおける音源方向に対するステアリングベクトルをもとに、ターゲット音を無歪みで出力するという拘束条件のもとで、マイクロホンアレーの周波数領域のNチャネルマイクロホン信号にフィルタ係数ベクトルを適用して得られる値が疎になるように前記フィルタ係数ベクトルを求める、
フィルタ係数算出装置。
請求項１のフィルタ係数算出装置であって、
fを周波数番号とし、lをフレーム番号とし、Lを1以上の整数の何れかとし、周波数領域のNチャネルマイクロホン信号をy(f,l)とし、フィルタ係数ベクトルをh(f,l)とし、h^H(f,l)をh(f,l)を転置し複素共役をとったベクトルとし、ステアリングベクトルをa(f)とし、ベクトルAのL1ノルムを|A|₁とし、

拘束条件
h^H(f,l)a(f)=1
により、フィルタ係数ベクトルh(f,l)を求める、
フィルタ係数算出装置。
請求項１または請求項２のフィルタ係数算出装置であって、
前記拘束条件は、
ターゲット音を無歪みで出力し、かつ、フィルタ係数ベクトルのノルムが所定の値以下であるという条件を含む、
フィルタ係数算出装置。
請求項３のフィルタ係数算出装置であって、
fを周波数番号とし、lをフレーム番号とし、Lを1以上の整数の何れかとし、周波数領域のNチャネルマイクロホン信号をy(f,l)とし、フィルタ係数ベクトルをh(f,l)とし、h^H(f,l)をh(f,l)を転置し複素共役をとったベクトルとし、ステアリングベクトルをa(f)とし、ベクトルAのL1ノルムを|A|₁とし、ベクトルAのL2ノルムを||A||とし、λを所定のパラメータとし、

拘束条件
h^H(f,l)a(f)=1

により、フィルタ係数ベクトルh(f,l)を求める、
フィルタ係数算出装置。
請求項１から請求項４の何れかのフィルタ係数算出装置を含む収音装置であって、
マイクロホンアレーの周波数領域のNチャネルマイクロホン信号に前記フィルタ係数ベクトルを適用して周波数領域の出力信号z(f,l)を求めるアレーフィルタリング部を含む、
収音装置。
Nを2以上の整数の何れかとし、N個のマイクロホンからなるマイクロホンアレーにおける音源方向に対するステアリングベクトルをもとに、ターゲット音を無歪みで出力するという拘束条件のもとで、マイクロホンアレーの周波数領域のNチャネルマイクロホン信号にフィルタ係数ベクトルを適用して得られる値が疎になるように前記フィルタ係数ベクトルを求める、
フィルタ係数算出方法。
請求項６のフィルタ係数算出方法を含む収音方法であって、
マイクロホンアレーの周波数領域のNチャネルマイクロホン信号に前記フィルタ係数ベクトルを適用して周波数領域の出力信号z(f,l)を求めるアレーフィルタリングステップを含む、
収音方法。
請求項１から請求項４の何れかのフィルタ係数算出装置、または、請求項５の収音装置としてコンピュータを機能させるためのプログラム。