(57)【要約】【課題】複数の信号が空間内で混合された混合信号から、源信号を正確に復元する。【解決手段】時間領域のセンサ信号を混合した混合信号ベクトルは周波数領域変換手段１０において短時間フーリエ変換により周波数毎の時系列を表現する観測信号ベクトルに変換され、次に分離信号算出手段２０において、観測信号ベクトルから周波数毎に分離信号ベクトルy(f,t)と基底ベクトルを列に持つ行列Aを算出し、パーミュテーション解決手段３０において、分離信号ベクトルy(f,t)の要素と行列Aの列を並べ替え、同じ信号源に対応するものが全ての周波数ｆで同じ添字iになるように並べ替え源信号を抽出する。【選択図】図２
2008227916.tif 000002
【特許請求の範囲】
【請求項１】
複数の信号源から発せられた信号が混合され、この混合された混合信号を各信号源からの信号に分離する信号分離装置において、
互いに異なる位置に設置された複数のセンサで観測された前記混合信号を周波数領域信号に変換する周波数領域変換手段と、
前記周波数領域信号に独立成分分析を適用し、周波数毎に分離信号ベクトル及び基底ベクトルを列に持つ行列を算出する分離信号算出手段と、
前記分離信号算出手段の算出結果を用いて周波数毎の観測信号における各分離信号の占有度を時間方向の系列として算出する占有度系列算出手段と、
前記占有度系列算出手段で算出された占有度系列に従って周波数毎の順列を算出するクラスタリング手段と、
前記クラスタリング手段で得られた順列に従って前記分離信号ベクトルの要素と基底ベクトルを並べ替える並べ替え手段と、
を備えることを特徴とする信号分離装置。
【請求項２】
請求項１記載の信号分離装置において、前記占有度系列算出手段ではｉ番目の分離信号のパワーとすべての分離信号のパワーの総和との割合を指標として占有度系列を算出することを特徴とする信号分離装置。
【請求項３】
請求項１記載の信号分離装置において、前記占有度系列算出手段は観測信号ベクトル白色化空間での基底ベクトルと観測信号ベクトルのコサイン距離を指標として、占有度系列を算出する構成されることを特徴とする信号分離装置。
【請求項４】
請求項１乃至３の何れかに記載の信号分離装置において、前記クラスタリング手段は、前記占有度系列と信号源に対応するセントロイドとの相関係数を全ての信号源および全ての周波数で足し合わせて得られる関数を最大化する条件で得られた順列を出力する大域的最適化手段を備えることを特徴とする信号分離装置。
【請求項５】
複数の信号源から発せられた信号が混合され、この混合された混合信号を各信号源からの信号に分離する信号分離方法において、
互いに異なる位置に設置された複数のセンサで観測された前記混合信号を周波数領域信号に変換する周波数領域変換処理と、
前記周波数領域信号に独立成分分析を適用し、周波数毎に分離信号ベクトル及び基底ベクトルを列に持つ行列を算出する分離信号算出処理と、
前記分離信号算出処理の算出結果を用いて周波数毎の観測信号における各分離信号の占有度を時間方向の系列として算出する占有度系列算出処理と、
前記占有度系列算出処理で算出された占有度系列に従って周波数毎の順列を算出するクラスタリング処理と、
前記クラスタリング処理で得られた順列に従って前記分離信号ベクトルの要素と基底ベクトルを並べ替える並べ替え処理と、
を含むことを特徴とする信号分離方法。
【請求項６】
請求項５記載の信号分離方法において、前記占有度系列算出処理はｉ番目の分離信号のパワーとすべての分離信号のパワーの総和との割合を指標として占有度系列を算出する処理とされることを特徴とする信号分離方法。
【請求項７】
請求項５記載の信号分離方法において、前記占有度系列算出処理は観測信号ベクトル白色化空間での基底ベクトルと観測信号ベクトルのコサイン距離を指標として、占有度系列を算出する処理されることを特徴とする信号分離方法。
【請求項８】
請求項５乃至７の何れかに記載の信号分離方法において、前記クラスタリング処理は、前記占有度系列と信号源に対応するセントロイドとの相関係数を全ての信号源および全ての周波数で足し合わせて得られる関数を最大化する条件で得られた順列を出力する大域的最適化処理を含むことを特徴とする信号分離方法。
【請求項９】
コンピュータが解読可能なプログラム言語によって記述され、コンピュータに請求項１乃至３の何れかに記載の信号分離装置として機能させる信号分離プログラム。
【請求項１０】
コンピュータが読み取り可能な記録媒体によって構成され、この記録媒体に請求項９記載の信号分離プログラムを記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号処理の技術分野に属し、特に複数の音波が空間内で混合されたものから、源音波に対応する源信号をできるだけ正確に復元する信号分離の技術に関する。
【背景技術】
【０００２】
［ブラインド信号分離］
以下では音波も信号と称して説明する。まず、ブラインド信号分離の定式化を行う。扱う信号は、あるサンプリング周波数fsでサンプリングされ、離散的に表現される。Ｎ個の信号が混合されてＭ個のセンサで観測されたとする。本発明では、信号の発生源からセンサまでの距離により信号が減衰・遅延し、また壁などにより信号が反射して伝送路歪みが発生しうる状況を扱う。このような状況の混合は、源信号sk(t)からセンサxj(t)へのインパルス応答hjk(r)による畳み込み混合
xj(t)＝Σk=1NΣr=0∞hjk(r)sk(t-r)          （１）
となる。ブラインド信号分離の目的は、源信号s1(t),...,sN(t)やインパルス応答h11(r
),...,h1N(r),...,hM1(r),...,hMN(r)を知らずに、観測信号x1(t),...,xM(t)のみから、源信号s1(t),...,sN(t)にそれぞれ対応する分離信号y1(t),...,yN(t)を求めることにある。
【０００３】
［周波数領域］
本発明では周波数領域において分離の操作を行う。そのために、センサでの観測信号xj(t)にL点の短時間フーリエ変換を適用して周波数毎の時間系列
【数１】
2008227916.tif 000003
を求める。ここでfは周波数であり、f=0,(1/2)fs,...,((L-1)/L)fsと離散化されている（fsはサンプリング周波数）。g(r)は窓関数である。ハニング窓g(r)=(1/2)(1+cos(2πr/L))などのg(0)にパワーの中心を持つ窓関数を用いることで、xj(f,t)は時刻tを中心とする観測信号xj(t)の周波数特性を表現する。xj(f,t)はLサンプルにわたる情報を含んでいるため、すべての時刻tに対してxj(f,t)を求める必要はなく、適当な感覚の時刻t毎にxj(f,t)を求める。また、複素共役の関係
【数２】
2008227916.tif 000004
から、周波数領域での処理においては、処理の対象とする周波数ビンの集合を
【数３】
2008227916.tif 000005
と絞り込むのが一般的である。
【０００４】
周波数領域で処理を行うと、式（１）で示される時間領域での畳み込み混合が、周波数領域では
【数４】
2008227916.tif 000006
と各周波数での単純混合に近似表現でき、分離の操作が単純になる。ここでhjk(f)は源信号skからセンサxjまでの周波数応答、sk(f,t)は式（２）と同様の式に従って源信号sk(t)に短時間フーリエ変換を施したものである。式（５）をベクトルを用いて表現すると、
【数５】
2008227916.tif 000007
となる。ここでx＝[x1,...,xM]Tは観測信号ベクトル、h＝[h1k,...,hMk]Tは源信号skから各センサへの周波数応答をまとめたベクトルである。
【０００５】
［独立成分分析による信号分離］
次に、独立成分分析(ICA:Independent Component Analysis)を用いて、周波数毎に信号分離を行う。周波数領域での処理により、畳み込み混合、式（１）を単純混合、式（５）およびそのベクトル表記、式（６）に近似できるため、ICAも単純混合のモデルに従って行う。ICAでは、観測信号ベクトルx(f,t)のみから、N行M列の分離行列W(f)および分離信号ベクトル
【数６】
2008227916.tif 000008
を算出する。ここで、分離信号ベクトルy＝[y1,...,yN]Tの各要素y1,...,yNが互いに独立になるように分離行列Wが算出される。このためのアルゴリズムは、［非特許文献１］などに様々なものが紹介されている。
【０００６】
次に、観測信号ベクトルxを基底ベクトルai=[a1i,...,aMi]Tの和で
【数７】
2008227916.tif 000009
を計算する。もし、分離行列Wの逆行列が存在すれば、行列AはA=W-1として計算できる。それ以外の場合は、二乗誤差の期待値E{‖x-Ay‖2}を最小にする
【数８】
2008227916.tif 000010
として算出される。ここで、[・H]はベクトル[・]の複素数共役転置、E{・}は行列[・]の要素毎の期待値から成る行列である。
【０００７】
［パーミュテーション問題］
ICAでは信号の独立性に着目して分離を行うため、得られる分離信号y1,...,yNには順序の任意性がある。その順序が入れ替わっても、独立性が保たれているからである。この任意性は、同じ源信号sk(t)に対応する分離信号成分yi(f,t)がすべての周波数fで同じ添え字iになるように解決されなければならない。これは、パーミュテーション問題と呼ばれ、周波数領域での信号分離において非常に重要な問題である。この問題を解決する従来技術は、周波数毎の分離信号のエンベロープの相関関数に基づくもの（［非特許文献２］、［非特許文献３］など）と、ICAの結果から信号源の方向や位置に相当する値を推定することに基づくもの（［特許文献１］）に大別できる。
【０００８】
周波数毎の分離信号のエンベロープを用いる従来技術では、分離信号yiのエンベロープを
【数９】
2008227916.tif 000011
などとして計算し、２つのエンベロープυiとυjのρ(υi,υj)を、相関、平均、標準偏差を、それぞれ、
【数１０】
2008227916.tif 000012
として計算する。相関係数ρは、−１から１までの間の値をとり、二つの系列が等しい時に１となる。二つのエンベロープυif,υkgの相関係数は、それらが同じ信号源に対応する場合に高くなると期待される。そのような傾向は、二つの周波数fとgが近いものであれば良く確認されるが、二つの周波数fとgが大きく異なる場合はあまり成り立たない。図９に二つの分離信号y1とy2のエンベロープを示す。図１０には、図９に示した分離信号エンベロープ間の相関係数を示す。図１０のＡとＤに示すように同じ音源の分離信号の相互でも、若干の数の周波数ペアで相関係数ρが０に近い値を取っており、
【数１１】
2008227916.tif 000013
となっていることが見受けられる。これは、音声信号のダイナミックレンジが広く、平均μi=0、分散σi2=1に正規化したとしても、信号がアクティブであるという現象が様々な値で表現されているからである。エンベロープの相関関数に基づくものは、局所的な周波数範囲でのパーミュテーション最適化には適しているが、周波数全体に渡って一貫性のある解を得るには不向きである。
【０００９】
一方、信号源の方向や位置に相当する値に基づくものは、周波数全体に渡って一貫性のある解をある程度の精度で容易に算出することができる。従って、これら２種類の手法を組み合わせた方法が［特許文献２］や［非特許文献４］などで提案されており、パーミュテーション問題への高精度な解を提供している。
【特許文献１】特許第３８８１３６７号明細書
【特許文献２】特開２００４−１４５１７２号公報
【非特許文献１】A. Hyvarinen and J. Karhunen and E Oja, “Independent Component Analysis,” John Wiley & Sons, 2001, ISBN 0-471-40540-X.
【非特許文献２】J. Anemuller and B. Kollmeier, “Amplitude modulation decorrelation for convolutive blind source separation,”in Proc. ICA 2000, June 2000, pp. 215-220.
【非特許文献３】N. Murata, S. Ikeda, and A. Ziehe, “An approach to blind source  separation based on temporal structure of speech signals,” Neurocomputing, vol. 41, pp. 1-24, Oct. 2001.
【非特許文献４】H. Sawada, R. Mukai, S. Araki, S. Makino, “A Robust and Precise Method for Solving the Permutation Problem of Frequency-Domain Blind Source Separation,” IEEE Trans. Speech and Audio Processing, vol.12, no. 5, pp. 530-538, Sep. 2004.
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、上述したような従来技術による信号分離方法が適切に動作するのは、信号源の方向や位置がある程度正確に推定できる状況に限られる。すなわち、反射の影響がそれほど大きくな状況で、しかも信号源からセンサまで信号が直接到達する経路が存在する場合に限られる。そのような条件を満たさない、例えば非常に残響時間の長い部屋で混合された音を分離するような状況では、信号源の位置や方向の推定値が非常に不正確になるため、従来のような方法では利用しにくい。
【００１１】
従来技術である分離信号エンベロープの相関係数に基づくパーミュテーション問題の解法は、信号源が全周波数に渡って同じ様な振幅のエンベロープを持たない限り、周波数全体に渡って一貫性のある解を得るには不向きである。また、ICAの結果から信号源の方向や位置に相当する値を推定する方法は、反射などの影響が強い場合には、信号源の方向や位置が正確に維持できずに性能が劣化する。
【００１２】
本発明の目的は、反射などの影響が強い場合でも、また、信号源が全周波数に渡って同じような振幅のエンベロープを持たない場合でも、正確にパーミュテーションを解決できる信号分離装置及び信号分離方法を提供することにある。
【課題を解決するための手段】
【００１３】
従来技術の課題を解決するために本発明では、ICAにより算出された分離信号が元の観測信号をどれぐらい占有しているかを示す値（以下、これを占有度と呼ぶ）を計算し、その値の時間方向の系列を従来技術のエンベロープに代わるものとして用いる。その新たな系列は、同じ信号源に対応する分離信号に関して、多くの周波数の組み合わせで高い相関係数を示すという特徴を持つ。本発明の原理的な詳細を以下に示す。
【００１４】
占有度を示す具体的な値を算出するために、あるJ番目のセンサの観測信号xJに着目する。式（９）のJ番目の要素のみを表現すると
【数１２】
2008227916.tif 000014
ある時間周波数スロット(f,t)において、ｉ番目の分離信号yiの占有度が高ければ式（１３）の値は１に近くなり、逆に、占有度が低ければ式（１３）の値は０に近くなる。
【００１５】
周波数毎の分離信号yi(f,t)の特徴を表現するものとして、従来技術で用いられている分離信号のエンベロープ式（１１）の代わりに式（１３）を用いてみる。
【数１３】
2008227916.tif 000015
すると、従来技術で説明した図９と同じ分離信号に対応する式（１４）の値は、図６に示すものとなり、それらの間の相関係数は図７に示すものとなる。図７を見ると、従来技術と比べて図７のＡとＤに示すように同一信号源間の相関係数が高くなる場合が多く、また、異なる信号源間の場合は図７のＢとＣに示すように相関係数が小さく（多くの場合負の値、白色部分）なっていることが分かる。従って、式（１４）によって算出される占有度の値は、従来のエンベロープよりも明確に、分離信号間の同一信号源依存性を示すものとなる。それは、以下の二つの理由による。
【００１６】
１点目は、式（１３）によって算出される占有度の値の範囲は０≦powRatio≦１と限定されており、ある分離信号が（他の分離信号と比べて）アクティブであるという事実は１に近い値によって一様に表現されるということである。従来技術のエンベロープ式（１１）では、それらを平均０，分散１に正規化したとしても、様々な値によって分離信号がアクティブであるという事実が表現されていたため、同一信号源得も相関係数が必ずしも高くならなかった。
【００１７】
２つ目は、式（１３）によって算出される占有度の値が、排他的であるという点である。図６を見ると、同じ時間周波数において、一方の分離信号の占有度の値が１（図では黒で表示される）に近ければ、他方の分離信号の占有度の値は０（図では白で表示される）に近くなっていることが分かる。これはある分離信号が観測信号をどれくらい占有しているか、ということを式（１３）が表現しているためであり、一方の占有度が高ければ、当然、他方の占有度は下がる。この排他的性質により、異なる信号源に対応する分離信号間での相関係数は、負の値になることが多い。
【００１８】
上述のように計算された占有度の系列が与えられた場合に、どのようにして最適なパーミュテーションを効率良く求めるかということも重要である。本発明では、以下の２段階から成る最適化手順も提供する。１段階目は、全ての周波数を同時に考慮する大域的な最適化である。占有度の系列は、多くの周波数の組み合わせで高い相関係数を示すため、そのような最適化方法は有効に働く。２段階目は局所的な最適化である。前段の大域的な最適化が、ある程度性質の良い解を生成した後に、局所的な最適化を行うことで、さらに性能を高めることが可能となる。
【００１９】
以上説明した本発明の特徴とする機能を実現する信号分離装置の構成は、複数の信号源から発せられた信号が混合され、この混合された混合信号を各信号源からの信号に分離する信号分離装置であって、互いに異なる位置に設置された複数のセンサで観測された混合信号を周波数領域信号に変換する周波数領域変換手段と、周波数領域信号に独立成分分析を適用し、周波数毎に分離信号ベクトル及び基底ベクトルを列に持つ行列を算出する分離信号算出手段と、分離信号算出手段の算出結果を用いて周波数毎の観測信号における各分離信号の占有度を時間方向の系列として算出する占有度系列算出手段と、占有度系列算出手段で算出された占有度系列に従って周波数毎の順列を算出するクラスタリング手段と、クラスタリング手段で得られた順列に従って分離信号ベクトルの要素と基底ベクトルを並べ替える並べ替え手段とを備えることを特徴とする。
【００２０】
更に、前記記載の信号分離装置において、占有度系列算出手段はｉ番目の分離信号のパワーとすべての分離信号のパワーの総和との割合を指標として占有度系列を算出することを特徴とする。
【００２１】
更に、前記記載の信号分離装置において、占有度系列算出手段は観測信号ベクトル白色化空間での基底ベクトルと観測信号ベクトルのコサイン距離を指標として、占有度系列を算出する構成されることを特徴とする。
【００２２】
更に、本発明では前記記載の信号分離装置において、クラスタリング手段は、占有度系列と信号源に対応するセントロイドとの相関係数を全ての信号源および全ての周波数で足し合わせて得られる関数を最大化する条件で得られた順列を出力する大域的最適化手段を備える構成を特徴とし、更に望ましくは大域的最適化手段で求めた順列を、それぞれの周波数について、それぞれの周波数と局所関係にある周波数との相関係数が最大化される条件で求められる順列に修正する局所的最適化手段とを備えることを特徴とする。
【００２３】
更に、本発明による信号分離方法は複数の信号源から発せられた信号が混合され、この混合された混合信号を各信号源からの信号に分離する信号分離方法にであって、互いに異なる位置に設置された複数のセンサで観測された混合信号を周波数領域信号に変換する周波数領域変換処理と、周波数領域信号に独立成分分析を適用し、周波数毎に分離信号ベクトル及び基底ベクトルを列に持つ行列を算出する分離信号算出処理と、分離信号算出処理の算出結果を用いて周波数毎の観測信号における各分離信号の占有度を時間方向の系列として算出する占有度系列算出処理と、占有度系列算出処理で算出された占有度系列に従って周波数毎の順列を算出するクラスタリング処理と、クラスタリング処理で得られた順列に従って分離信号ベクトルの要素と基底ベクトルを並べ替える並べ替え処理とを含むことを特徴とする。
【００２４】
更に、本発明による信号分離方法は前記記載の信号分離方法において、占有度系列算出処理はｉ番目の分離信号のパワーとすべての分離信号のパワーの総和との割合を指標として占有度系列を算出する処理とされることを特徴とする。
【００２５】
更に本発明による信号分離方法は、前記記載の信号分離方法において、占有度系列算出処理は観測信号ベクトル白色化空間での基底ベクトルと観測信号ベクトルのコサイン距離を指標として、占有度系列を算出する処理されることを特徴とする。
【００２６】
更に、本発明による信号分離方法は前記記載の信号分離方法において、クラスタリング処理は、占有度系列と信号源に対応するセントロイドとの相関係数を全ての信号源および全ての周波数で足し合わせて得られる関数を最大化する条件で得られた順列を出力する大域的最適化処理を含むことを特徴とし、更に望ましくは大域的最適化処理で求めた順列を、それぞれの周波数について、それぞれの周波数と局所関係にある周波数との相関係数が最大化される条件で求められる順列に修正する局所的最適化処理とを含むことを特徴とする。
【発明の効果】
【００２７】
本発明による信号分離装置及び信号分離方法によれば、様々な妨害信号が発生する実環境において、目的の信号を精度良く取り出すことが可能となる。音信号に対する応用例としては、音声認識器のフロントエンドとして働く音源分離システムなどが挙げられる。話者とマイクが離れた位置にあり、マイクが話者の音声以外を収音してしまうような状況でも、本発明による信号分離装置を使うことで、話者の音声のみを取り出して正しく音声を認識することができる。
【発明を実施するための最良の形態】
【００２８】
本発明による信号分離装置は全てをハードウェアによって構成することも可能であるが、最も簡素に実現するには、コンピュータに本発明による信号分離プログラムをインストールし、コンピュータに本発明による信号分離装置として機能させる実施形態が最良の実施形態である。
【００２９】
コンピュータに本発明による信号分離装置として機能させるには、コンピュータにインストールした信号分離プログラムにより、コンピュータに互いに異なる位置に設置された複数のセンサで観測された混合信号を周波数領域信号に変換する周波数領域変換手段と、周波数領域信号に独立成分分析を適用し、周波数毎に分離信号ベクトル及び基底ベクトルを列に持つ行列を算出する分離信号算出手段と、分離信号算出手段の算出結果を用いて周波数毎の観測信号における各分離信号の占有度を時間方向の系列として算出する占有度系列算出手段と、占有度系列算出手段で算出された占有度系列に従って周波数毎の順列を算出するクラスタリング手段と、クラスタリング手段で得られた順列に従って前記分離信号ベクトルの要素と基底ベクトルを並べ替える並べ替え手段とを構築し、信号分離装置として機能させる。
【実施例１】
【００３０】
図１に本発明による信号分離装置の実施例を示す。図中１００は本発明による信号分離装置の全体を示す。本発明による信号分離装置１００は、周波数領域変換手段１０と、分離信号算出手段２０と、パーミュテーション解決手段３０と、スケーリング調整手段４０と、時間領域変換手段５０とによって構成される。
パーミュテーション解決手段３０は図２に示すように、占有度算出手段３０−１と、クラスタリング部３０−２と、並べ替え手段３０−３とによって構成される。
【００３１】
更に、クラスタリング手段３０−２は図３に示すように、大域的最適化手段３０−２−１と、局所的最適化手段３０−２−２とによって構成される。尚、図３に示す実施例ではクラスタリング手段３０−２を大域的最適化手段３０−２−１と局所的最適化手段３０−２−２とによって構成した場合を示したが、後に説明する実験例で明らかなように、大域的最適化手段３０−２−１のみによるクラスタリング処理によっても優れた信号分離性能を得ることができるから、この発明では局所的最適化手段３２−２を必須要件として含めないこととする。
【００３２】
以下に各部の構成及び動作について説明する。
時間領域のセンサ信号をまとめたベクトル[x1(t),...,xM(t)]Tは、周波数領域変換手段１０において、短時間フーリエ変換により周波数毎の時系列を表現する観測信号ベクトルx(f,t)＝[x1(f,t),...,xM(f,t)]Tに変換される。次に、分離信号算出手段２０において、観測信号ベクトルx(f,t)から周波数f毎に分離信号ベクトルy(f,t)＝[y1(f,t),...,yN(f,t)]Tと基底ベクトルを列に持つ行列A(f)=[a1(f),...,aN(f)]を算出する。
【００３３】
次に、パーミュテーション解決手段３０において、分離信号ベクトルの要素y(f,t)と行列Aの列ai(f)を並べ替え、同じ信号源に対応するものがすべての周波数fで同じ添字iとなるようにする。そして、スケーリング調整手段４０で分離信号y1(f,t),...,yN(f,t)のスケーリングを調整し、最後に、時間領域変換手段５０において、それら分離信号yk(f,t)を短時間逆フーリエ変換して、時間領域での分離信号ベクトル[y1(t),...,yN(t)]Tを得る。
【００３４】
（周波数領域変換手段）
周波数領域変換手段１０では、従来技術の説明で記載したように、式（２）に従って短時間フーリエ変換を行う。
【００３５】
（分離信号算出手段）
分離信号算出手段２０では、周波数f毎に処理を行う。まず、従来技術の説明に記載したように、式（７）に示すICAを適用し、分離信号ベクトルyを算出する。次に、式（１０）を用いるか、あるいは分離行列Wの逆行列として、基底ベクトルを列に持つ行列Aを算出する。
【００３６】
（パーミュテーション解決手段）
図２にパーミュテーション解決手段３０の構成を示す。パーミュテーション解決手段３０はこの発明で特徴とする占有度算出手段３０−１と、クラスタリング部３０−２と、並べ替え手段３０−３とによって構成される。
この発明で特徴とする占有度算出手段３０−１では従来技術のエンベロープに代わるものとして分離信号の占有度を示す値を系列υif(t)として計算する。ここでのクラスタリングは、一般的なものとは異なり、クラスタリング結果は周波数毎の順列Πfで表現される。最後に並べ替え手段３０−３において、順列Πfに従って分離ベクトルy(f,t)の要素と基底ベクトルの並びを変更する。
【００３７】
（占有度算出手段）
占有度算出手段３０−１では、分離信号y(f,t)が観測信号x(f,t)をどれほど占有しているかを示す値を計算する。その方法の一つとして、この実施例では、課題を解決するための手段の項で説明した複数のセンサの中の特定したセンサにおけるｉ番目の分離信号のパワーとすべての分離信号のパワーの総和との割合を占有度として求める式（１３）或は全てのセンサに着目して占有度を求める式（１３’）、或は部分集合に該当するセンサに着目して占有度を求める方法を提案する。
【００３８】
（クラスタリング手段）
占有度算出手段３０−１で算出された占有度系列υif(t)は、次にクラスタリング手段３０−２に入力される。クラスタリング手段３０−２の構成を図３に示す。クラスタリング手段３０−２はこの実施例では大域的最適化手段３０−２−１と、局所的最適化手段３０−２−２とによって構成した場合を示す。
【００３９】
大域的最適化手段３０−２−１では信号源k毎にセントロイド（クラスタの中心ベクトル）ckを推定する形でクラスタリングを行う。最大化する目的関数として
【数１４】
2008227916.tif 000016
を考える。これは、順列Πfによって並べ替えられた占有度系列υif,i=Πf(k)と信号源kに対応するセントロイドckとの相関係数ρを全ての信号源および全ての周波数で足し合わせたものである。この目的関数は、良く知られたk-means法と同じように、セントロイドckと順列Πfを交互に最適化することで、最大化することができる。セントロイドckは、現在の順列Πfによって並べ替えられた占有度系列υifの値の周波数方向の平均値として算出される。
【数１５】
2008227916.tif 000017
式（１７）と式（１８）を、変化が起きなくなるまで繰り返すことで、目的関数を表わす式（１６）が最大化される。
【００４０】
尚、クラスタリング部３０−２を大域的最適化手段３０−２−１のみで構成する場合は上述した大域的最適化手段３０−２−１が算出した順列Πfを最終値として出力し、並べ替え手段３０−３に入力する。ただし、ここでは局所的最適化手段３０−２−２を用いる場合を例示しているから、大域的最適化手段３０−２−１で算出した順列Πfは暫定値として局所的最適化手段３０−２−２に入力される。
【００４１】
次に、局所的最適化手段３０−２−２では、それぞれの周波数fについて、各周波数の近隣に位置する局所的範囲内のいくつかの特定関係にある周波数との相関係数ρを、以下の式に従って順次最大化していく。
【数１６】
2008227916.tif 000018
特定関係にある周波数の集合R(f)は、同じ信号源である場合にυifとυi'gの相関係数が高くなるような周波数gを要素に持つようなものが望ましい。典型的には、近傍の集合Α(f)と倍音関係にある周波数の集合Η(f)を考慮し、R(f)＝Α(f)∪Η(f)として特定関係にある周波数gを定義する。近傍の周波数の集合Αは、例えば、
Α(f)=｛f−3Δf，f−2Δf，f−Δf，f＋Δf，f＋2Δf，f＋3Δf｝
として定義できる。ここで、Δf＝(1/L)fsは、隣り合う周波数ビン間の周波数の差である。また、倍音関係にある周波数の集合Ηは、例えば、
Η(f)＝｛round(f/2)−Δf，round(f/2)，round(f/2)＋Δf，2f−Δf，2f，2f＋Δf｝
として定義できる。ここで、round(・)は、周波数の集合Fから・に最も近い周波数を意味する。局所的最適化部では、式（１９）に従った順列Πfの最適化を、すべての周波数で改善が起こらなくなるまで繰り返す。
【００４２】
［並べ替え手段］
並べ替え手段３０−３では、クラスタリング手段３０−２で算出された順列Πfに基づいて、分離信号ベクトルy＝[y1,...,yN]Tの要素と、基底ベクトルから成る行列A＝[a1,...,aN]の列を以下のように並べ替える。
yk(f,t)←yΠf(k)(f,t)，ak(f)←aΠf(k)(f),∀k,f,t
以上でパーミュテーション解決手段３０の説明を終わる。
【００４３】
［スケーリング調整手段４０］
ICAでは信号の独立性に着目して分離を行うため、得られる分離信号には、振幅や移相、すなわちスケーリングの任意性がある。時間領域に戻した際に適切な分離信号となるためには、振幅や位相を適切に設定する必要がある。この実施例では、あるセンサJでの観測信号に合わせるという考え方を採用し、
yk(f,t)←ajk(f)yk(f,t),∀k,f,t
によりスケーリング調整を行う。
【００４４】
［時間領域変換手段］
最後に時間領域変換手段５０で、これまでに得られた周波数領域の分離信号に、短時間逆フーリエ変換などを施して、時間領域の分離信号を作成する。
【実施例２】
【００４５】
図４および図５を用いて本発明の実施例２を説明する。この実施例２ではパーミュテーション解決手段３０に観測信号ベクトルx(f,t)と、分離信号y(f,t)と、基底ベクトルを列に持つ行列A(f)とを入力し、主に占有度算出手段３０−１’は図５に示すように観測信号ベクトルx(f,t)と、基底ベクトルを列に持つ行列A(f,t)により占有度系列υif(t)を算出する構成とするものである。
つまり、実施例１では特定のセンサ番号Jを選択し、選択したセンサにおける分離信号のパワーと分離信号の総パワーの比により占有度系列υif(t)を算出する構成としたから、或るセンサJを指定しなければならないが、状況によっては不適切なセンサ（例えば感度の悪いセンサ）を選んでしまう恐れがある。
【００４６】
この課題を解決するために、占有度の算出方法としてこの実施例２では観測信号ベクトルxが白色化された空間での基底ベクトルaiと観測信号ベクトルxのコサイン距離
【数１７】
2008227916.tif 000019
を考えることもできる。ここで、
bi＝Vai，z＝Vx
は、観測信号ベクトル白色化空間での、それぞれ、biは白色化された基底ベクトル、zは白色化された観測信号ベクトルであり、白色化行列VはVHV＝(E{xxH})-1を満たすものとして、例えば、固有値分解E{xxH}＝EDEHを行った後
V＝D-1/2EH
として計算できる。式（１３）及び（１３’）では分離信号ベクトルyと基底ベクトルを列に持つ行列Aを用いるのに対し、この式（１５）では図５に記載されている観測信号ベクトルxと行列Aを用い、分離信号ベクトルyは占有度の算出には用いない。また、ここでは観測信号ベクトルｘを演算指標として用いるから、全てのセンサの観測信号を対象としている。従ってこの実施例によれば特定のセンサを選択する必要はない。
【００４７】
式（１３）で算出される値の性質と同様、式（１５）の値の範囲も０≦cosdw(ai,x)≦１と限定されている。そして、i番目の分離信号yiの占有度が高ければ式（１５）の値は１に近くなり、逆に、占有度が低ければ式（１５）の値は０に近くなる。式（１５）の計算過程に分離信号yiは直接関わっていないが、式（８）に示す通り、基底ベクトルaiと分離信号は密接に関連しているため、式（１５）も分離信号yiの性質を表現する値を算出していることになる。
【００４８】
尚、基準センサの影響を受けないようにする他の方法として、式（１３）を次式（１３’）に書き換えることにより、全てのセンサの分離信号を演算対象とするため特定のセンサの影響を除去することができる。また、複数のセンサの中の一部（部分集合）を演算対象とすることもできる。この場合も、特定のセンサの影響を受ける不都合を解消することができる。
【数１８】
2008227916.tif 000020
【実施例３】
【００４９】
図８に本発明の実施例３を示す。この実施例ではコンピュータに本発明による信号分離プログラムをインストールし、インストールしたプログラムによりコンピュータ内に周波数領域変換手段１０と、分離信号算出手段２０と、パーミュテーション解決手段３０と、スケーリング調整手段４０と、時間領域変換手段５０とを構築し、信号分離装置１００として機能させる実施例を示す。
【００５０】
つまり、コンピュータは周知のように中央演算処理装置（以下CPUと称す）１０１と、読み出し専用メモリ（以下ROMと称す）１０２と、プログラムを記録し、記録したプログラムをCPU１０１で解読する。書き替え可能なメモリ（以下RAM）１０３と、起動時にプログラム等をRAM１０３に展開する例えばハードディスクのような記録装置１０４と、入力ポート１０５、出力ポート１０６等により構成される。
【００５１】
RAM１０３に周波数領域変換手段１０、分離信号算出手段２０、パーミュテーション解決手段３０，スケーリング調整手段４０、時間領域変換手段５０と記載した部分はそれぞれ、これらの手段を構築するためのプログラムを記憶した領域を指し、各領域に記憶したプログラムをCPU１０１が解読し、コンピュータが各手段として機能する。
入力ポート１０５には特に図示していないが例えば磁気ディスク読取装置、CD-ROM読取装置、メモリカード読取装置、或いはモデム等を接続し、コンピュータが読取可能な記録媒体から入力ポート１０５を通じて記録装置１０４にプログラムをインストールする。
【００５２】
図示する例では入力ポート１０５にマウスのような入力手段２０１とマイク群２０２を接続し、このマイク群２０２で受音した観測信号を周波数領域変換手段１０に取り込むように構成し、更に、また出力ポート１０５にはスピーカ３０１を接続し、スピーカ３０１で分離した源音を再生する構成とした場合を示す。また、スピーカ３０１の他に、モニタ３０２、プリンタ３０３等を出力ポート１０６に接続した状態を示している。
以上により本発明による信号分離プログラムをコンピュータにインストールすることにより、コンピュータを信号分離装置として機能させることができることが理解できよう。
【００５３】
（実験結果）
本発明の効果を示すために、図１１に示す実験条件と図１２に示すマイクロホンＭ１〜Ｍ３とスピーカＳＰ１〜ＳＰ３の配置を用いて実験を行った。３つの音を同時に鳴らした時の混合音を３つのマイクロホンＭ１〜Ｍ３で観測し、その観測信号のみからそれぞれの音に対応する分離信号を算出するという問題設定である。様々な音声信号の組合せで評価できるように、スピーカＳＰ１〜ＳＰ３からマイクロホンＭ１〜Ｍ３までのインパルス応答を測定し、音声信号をインパルス応答に畳み込んで混合することで観測信号を生成した。図１２に示す配置では、二つの音源が同じ方向からマイクロホンＭ１〜Ｍ３に到来し、しかもマイクロホンＭ１〜Ｍ３の間隔は４ｃｍと狭いため、信号源の方向や位置に相当するICAの結果から推定するパーミュテーションの解法（［特許文献１］、［特許文献２］、［非特許文献４］に記載）は、適用しにくい。分離性能は、signal-to-interference ratio  (SIR)の改善量で評価した。これは、各出力i毎に、出力SIRと入力SIRの差OutputSIRi-InputSIRiとして計算される。入力SIRと出力SIRは、それぞれ以下の式で計算される。
【数１９】
2008227916.tif 000021
ここで、J∈｛1,...,M｝はある選択された基準センサの番号を示す。また、源信号skのうち出力yi(t)に出てきた成分をyik(t)と表記する。この定義により、yi(t)=Σk=1Nyik(t)が満たされる。
【００５４】
３つの音声の組合せを８通り用いて実験を行った。図１３に、種々のパーミュテーション解決法を用いた結果を、SIR改善量の全出力に関する平均値で示す。省略形として、“Env”と“Cos”は、それぞれ、式（１１）に示す従来のエンベロープ|yi|を用いた場合と、本発明で提案した式（１５）に示すcosdw(ai,x)の値を用いた場合を示す。また、省略形“Gl”と“Lo”は、それぞれ、大域的最適化と局所的最適化を示す。“Optimal”と書かれたものは、信号源に関する情報を用いて最適なパーミュテーションを算出したものである。現実的な状況では、そのような情報は得られないが、性能の上限を示す目的で掲載した。
【００５５】
図１３の結果から、以下のことが考察できる。大域的最適化“Gl”は、提案した“Cos”を用いた場合には良好に働くが、従来の“Env”では良い結果が得られていない。局所的最適化“Lo”のみを適用した場合は、双方“Env”、“Cos”の特徴系列について、良い結果は得られない。しかし、図１３に示すCos(Gl)が示す値から明らかなように大域的最適化“Gl”のみによっても、従来技術では得られない高い評価値を示す。望ましくは大域的最適化“Gl”によってある程度良い解を得た後に局所的最適化(Lo)を適用すれば、効果的に解を改善できることが分かる。提案手法としての“Cos(Gl+Lo)”は、ほぼ最適解に近い分離性能を得ている。
【産業上の利用可能性】
【００５６】
様々な妨害信号が発生する実環境において、本発明の信号分離装置によれば目的の信号を精度良く取り出すことが可能となる。応用例としては音声認識器の前処理装置として働く信号分離装置として活用できる。
【図面の簡単な説明】
【００５７】
【図１】本発明の信号分離装置の実施例１を説明するためのブロック図。
【図２】本発明の信号分離装置に用いるパーミュテーション解決手段３０の構成を説明するためのブロック図。
【図３】図２に示したパーミュテーション解決手段に用いるクラスタリング手段の構成例を説明するためのブロック図。
【図４】本発明の信号分離装置の実施例２を説明するためのブロック図。
【図５】本発明の実施例２に用いるパーミュテーション解決手段の構成を説明するためのブロック図。
【図６】本発明の作用効果を説明するためのグラフ。
【図７】図６と同様に本発明の作用効果を説明するためのグラフ。
【図８】本発明による信号分離装置をコンピュータで構成した場合の、コンピュータの内部の様子を説明するためのブロック図。
【図９】従来技術の動作を説明するためのグラフ。
【図１０】従来技術の分離性能を説明するためのグラフ。
【図１１】本発明による信号分離装置の作用効果を実証するために行なった実験の条件を説明するための図。
【図１２】図１１に示した実験条件を図で示した配置図。
【図１３】本発明と従来技術との分離性能を比較するためのグラフ。
【符号の説明】
【００５８】
１００    信号分離装置              ３０−１、３０−１’    占有度算出手段
１０    周波数領域変換手段                    ３０−２    クラスタリング手段
２０    分離信号算出手段                      ３０−３    並べ替え手段
３０    パーミュテーション解決手段        ３０−２−１    大域的最適化手段
４０    スケーリング調整手段              ３０−２−２    局所的最適化手段
５０    時間領域変換手段
【図１】
2008227916.tif 000022
【図２】
2008227916.tif 000023
【図３】
2008227916.tif 000024
【図４】
2008227916.tif 000025
【図５】
2008227916.tif 000026
【図８】
2008227916.tif 000027
【図１１】
2008227916.tif 000028
【図１２】
2008227916.tif 000029
【図１３】
2008227916.tif 000030
【図６】
2008227916.tif 000031
【図７】
2008227916.tif 000032
【図９】
2008227916.tif 000033
【図１０】
2008227916.tif 000034