(57)【特許請求の範囲】
【請求項１】
第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の移動しうる音源からの音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力する混合音時間周波数展開部と、
前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、
前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分を要素に持つテンソルと、
前記複数の音源の各々から各マイクまでの伝達特性を表すステアリングベクトルの絶対値成分、前記複数の音源の各々についての前記基底パワースペクトル、及び前記複数の音源の各々についての、各時刻及び各基底におけるパワーを表すアクティベーションパラメータから求められるテンソルとの距離を表す目的関数を最小化するように、
前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータとの更新の各々を行うパラメータ更新部と、
予め定められた収束条件を満たしたか否かを判定し、前記収束条件を満たすまで、前記パラメータ更新部による更新を繰り返し行う収束判定部と、
を含む信号解析装置。
【請求項２】
前記ステアリングベクトルは、前記絶対値成分と偏角成分に式(10)の通り分解され、
0006519801.tif 000116
前記式(10)に基づき、前記観測時間周波数成分の瞬時混合モデルは、式(11)と式(12)の積で表現される
0006519801.tif 000117
請求項１記載の信号解析装置。
ただし、ｉは音源インデックスを表し、ωは周波数インデックスを表し、
0006519801.tif 000118
は虚数単位を表し、
0006519801.tif 000119
は、各音源の複素スペクトログラムを表し、τは時刻フレームインデックスを表す。
【請求項３】
前記偏角成分を式(13)の通りとすることで、
0006519801.tif 000120
前記瞬時混合モデルを式(14)と変形し、
0006519801.tif 000121
上記式(14)に基づき前記観測信号を解析する
請求項２記載の信号解析装置。
ただし、ｐｉ，ω，τは、各音源のパワースペクトログラムの期待値を表すパラメータであり、
0006519801.tif 000122
は、（２×２）の単位行列である。
【請求項４】
前記パラメータ更新部は、前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、前記目的関数を最小化するように、前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータと、前記基底パワースペクトルとの更新の各々を行う請求項１〜３の何れか１項記載の信号解析装置。
【請求項５】
前記第１のマイクを、音声を体表から取得する非可聴つぶやきマイクとした請求項１〜４の何れか１項記載の信号解析装置。
【請求項６】
前記距離の距離尺度を板倉齋藤擬距離基準とする請求項１〜５の何れか１項記載の信号解析装置。
【請求項７】
混合音時間周波数展開部と、パラメータ更新部と、収束判定部とを含む信号解析装置における、信号解析方法であって、
前記混合音時間周波数展開部は、第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の移動しうる音源からの音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、
前記パラメータ更新部は、前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、
前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分を要素に持つテンソルと、
前記複数の音源の各々から各マイクまでの伝達特性を表すステアリングベクトルの絶対値成分、前記複数の音源の各々についての前記基底パワースペクトル、及び前記複数の音源の各々についての、各時刻及び各基底におけるパワーを表すアクティベーションパラメータから求められるテンソルとの距離を表す目的関数を最小化するように、
前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータとの更新の各々を行い、
前記収束判定部は、予め定められた収束条件を満たしたか否かを判定し、前記収束条件を満たすまで、前記パラメータ更新部による更新を繰り返し行う
信号解析方法。
【請求項８】
前記パラメータ更新部により更新の各々を行うことは、前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、前記目的関数を最小化するように、前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータと、前記基底パワースペクトルとの更新の各々を行う請求項７記載の信号解析方法。
【請求項９】
前記第１のマイクを、音声を体表から取得する非可聴つぶやきマイクとした請求項７又は８記載の信号解析方法。
【請求項１０】
コンピュータを、請求項１〜６の何れか１項記載の信号解析装置を構成する各部として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号解析装置、方法、及びプログラムに係り、特に、パラメータを推定する信号解析装置、方法、及びプログラムに関する。
【背景技術】
【０００２】
音声コミュニケーションは、我々の生活における最も効率的なコミュニケーション手段の一つであり、携帯電話などの通信端末の普及により、異なる環境下にいる人物との会話も容易となっている。
【０００３】
しかしながら、聞き手が周囲にいない状況での発声行為は周囲の人物にとって迷惑になりやすい、また、会話内容が漏えいしてしまう、といった理由から発声を躊躇するような状況が数多く存在する。
【０００４】
このような問題に対し、周囲が聴取可能な声を発することなく音声コミュニケーションを行う技術として、サイレント音声インタフェースが注目を浴びている。
【０００５】
非可聴つぶやき(Non-Audible Murmur:NAM)マイク（以後、ＮＡＭマイクとする）は、そのようなサイレント音声インタフェースの一つであり、周囲が聴取困難なほど微弱なささやき声を体表から収録することを可能にする。
【０００６】
また、ＮＡＭマイクは、微弱な音声だけでなく、通常音声やささやき声など様々な音声を収録することができ、次世代の音声インタフェースとして期待されている。
【０００７】
ＮＡＭマイクによる音声の体内伝導収録は、マイクの構造上、外部雑音に対して比較的頑健である。通常音声を収録する場合、一般的な空気伝導マイクによる収録と比較して、外部雑音に対して高いＳＮ比を確保できる。一方で、微弱な音声を収録する場合、十分なＳＮ比を得ることは困難となり、外部雑音によって収録音声の品質が容易に劣化する。
【０００８】
雑音環境下では、マイクを口唇付近に設置しない限り、非常に微弱な音声を空気伝導音声として収録することは困難である。このような発声音の微弱性に着目し、図７に示すような位置に空気伝導マイクを配置することで、口唇から放射される音声が十分に抑えられた、かつ、ＮＡＭマイクに混入するものと近い音響特性を有する外部雑音を収録できる。
【０００９】
このとき、時刻ｔにおける体内伝導微弱音声をｓ１，ｔ、空気伝導外部雑音をｓ２，ｔとすると、ＮＡＭマイクおよび空気伝導マイクにおける観測信号は、下記（１）式、及び（２）式のように表される。なお、下記（１）式がＮＡＭマイクにおける観測信号であり、下記（２）式が空気伝導マイクにおける観測信号である。
【００１０】
【数１】
0006519801.tif 000002
【００１１】
したがって、上記（２）式が成立し、空気伝導マイクで観測した信号が外部雑音の参照信号として利用できれば、古典的な適応アルゴリズムやセミブラインド音源分離（Semi-Blind Source Separation: Semi-BSS）の適用により、マイク間特性ατを推定することで体内伝導信号に含まれる外部雑音を抑圧することができる。
【００１２】
空気伝導信号を外部雑音の参照信号として扱い、ＮＡＭマイクに混入した外部雑音を抑圧することは、典型的なエコーキャンセリングと類似した問題である。本枠組みでは、微弱音声の発生時がエコーキャンセリングにおけるダブルトーク状態に相当する。したがって、Least Mean Square (LMS)アルゴリズムなどの古典的な適応アルゴリズムを用いる場合、外部雑音が混入した体内伝導信号に対して音声区間検出(Voice Activity Detection: VAD)を行い、発話時には更新処理を停止させるといった対処が必要になる。しかしながら、雑音環境下において、高精度なＶＡＤを実現することは容易ではない。近年、ダブルトーク対策が不要な手法としてSemi-BSSを適用した手法が提案されている（非特許文献１、及び非特許文献２）。
【００１３】
音源複素スペクトログラムを
【００１４】
【数２】
0006519801.tif 000003
【００１５】
観測複素スペクトログラムを
【００１６】
【数３】
0006519801.tif 000004
【００１７】
とする。このとき、分析フレーム超に対してマイク間特性のフィルタ長が短いと仮定すると、観測スペクトログラム
【００１８】
【数４】
0006519801.tif 000005
【００１９】
は、下記（３）式のように表すことができる。
【００２０】
【数５】
0006519801.tif 000006
【００２１】
ここで、ωは周波数インデックス、τは時刻フレームインデックスを表す。また、
【００２２】
【数６】
0006519801.tif 000007
【００２３】
は音源からマイクまでの伝達特性を表す（２×２）の混合行列であり、時間的に変動しないものとみなす。ＢＳＳの目的は、
【００２４】
【数７】
0006519801.tif 000008
【００２５】
の逆行列である分離行列
【００２６】
【数８】
0006519801.tif 000009
【００２７】
を観測信号
【００２８】
【数９】
0006519801.tif 000010
【００２９】
のみから推定することであり、音源の統計的独立性に基づく手法である独立成分分析（Independent Component Analysis:ICA）が広く用いられている。ここでは、ＩＣＡのアルゴリズムとして、自然勾配法を使用し、下記（４）式、及び（５）式による更新を繰り返すことで、分離信号間の独立性を最大にする分離行列
【００３０】
【数１０】
0006519801.tif 000011
【００３１】
を得る。
【００３２】
【数１１】
0006519801.tif 000012
【００３３】
ここで、
【００３４】
【数１２】
0006519801.tif 000013
【００３５】
は、（２×２）の単位行列、
【００３６】
【数１３】
0006519801.tif 000014
【００３７】
は分離信号の複素スペクトログラム、
【００３８】
【数１４】
0006519801.tif 000015
【００３９】
はステップ幅、
【００４０】
【数１５】
0006519801.tif 000016
【００４１】
はスコア関数と呼ばれる非２次的関数、
【００４２】
【数１６】
0006519801.tif 000017
【００４３】
は時間的平均演算を表す。
【００４４】
スコア関数には、目的信号の分布がスーパーガウシアンであることを仮定し、下記（６）式を用いる。
【００４５】
【数１７】
0006519801.tif 000018
ただし、
【数１８】
0006519801.tif 000019
は虚数単位を表す。
【００４６】
本枠組みでは、微弱音声がＮＡＭマイクでのみ観測されると仮定しているため、分離行列の要素の一部を下記（７）式のように固定でき、セミブラインドな問題となる。
【００４７】
【数１９】
0006519801.tif 000020
【００４８】
さらに、上記（１）式、及び（２）式に示した混合過程においては、ｗ１，１，ω＝１と固定できる。これにより、ｗ１，１，ωを更新する場合と比較して、周波数ＩＣＡにおけるスケールの不定性が解消され、逆投影法（Projection Back: PB）の適用が不要になる。さらに、ｗ２，２，ωについても同様に、ｗ２，２，ω＝１と同定できる。その結果、分離行列を下記（８）式のように設定することができ、この場合、推定すべき要素はｗ１，２，ωのみとなる。なお、雑音源から空気伝導マイクへのパスの逆フィルタが作れない場合を考慮すると、ｗ２，２，ωも更新する方が無難である。
【００４９】
【数２０】
0006519801.tif 000021
【先行技術文献】
【非特許文献】
【００５０】
【非特許文献１】田尻祐介, 戸田智基, Graham Neubig, Sakriani Sakti, 中村哲, “外部雑音モニタリングを用いた非可聴つぶやきに対する雑音抑圧法,” 日本音響学会秋季研究発表会講演論文集, pp. 591-592, Sep. 2015.
【非特許文献２】田尻祐介, 戸田智基, 中村哲, “外部雑音モニタリングに基づく体内伝導微弱音声に対する雑音抑圧法,”信学技報, vol. 115, no. 302, EA2015-31, pp. 41-46, Nov. 2015.
【発明の概要】
【発明が解決しようとする課題】
【００５１】
上記非特許文献１及び非特許文献２に記載の技術では、雑音源が一つの点音源である条件下でのみ、外部雑音モニタリングに基づく雑音抑圧法の有効性が確認されている。しかしながら、実環境においてそのような条件が成立することは希であり、複数音源の存在を考慮する必要があるという問題がある。
【００５２】
また、複数の雑音源がユーザに対して異なる方位に存在する場合、マイクごとに異なる振幅・位相で入力された雑音が加算されるため、空気伝導信号に対するフィルタリング処理によって、ＮＡＭマイクに混入した外部雑音を表現することは困難となり、従来法による雑音抑圧効果が期待できないという問題がある。
【００５３】
また、雑音源が一つであっても、その音源が移動する場合、上記と同様の問題が生じるという問題がある。
【００５４】
本発明は、上記の事情に鑑みてなされたもので、雑音源が複数、又は移動する場合であっても雑音を抑圧することができる信号解析装置、方法、及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【００５５】
上記目的を達成するために、第１の発明に係る信号解析装置は、第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の音源からの音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力する混合音時間周波数展開部と、前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分を要素に持つテンソルと、前記複数の音源の各々から各マイクまでの伝達特性を表すステアリングベクトルの絶対値成分、前記複数の音源の各々についての前記基底パワースペクトル、及び前記複数の音源の各々についての、各時刻及び各基底におけるパワーを表すアクティベーションパラメータから求められるテンソルとの距離を表す目的関数を最小化するように、前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータとの更新の各々を行うパラメータ更新部と、予め定められた収束条件を満たしか否かを判定し、前記収束条件を満たすまで、前記パラメータ更新部による更新を繰り返し行う収束判定部と、を含んで構成されている。
【００５６】
第２の発明に係る信号解析方法は、混合音時間周波数展開部と、パラメータ更新部と、収束判定部とを含む信号解析装置における、信号解析方法であって、前記混合音時間周波数展開部は、第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の音源からの音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、及び第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、前記パラメータ更新部は、前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分を要素に持つテンソルと、前記複数の音源の各々から各マイクまでの伝達特性を表すステアリングベクトルの絶対値成分、前記複数の音源の各々についての前記基底パワースペクトル、及び前記複数の音源の各々についての、各時刻及び各基底におけるパワーを表すアクティベーションパラメータから求められるテンソルとの距離を表す目的関数を最小化するように、前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータとの更新の各々を行い、前記収束判定部は、予め定められた収束条件を満たしたか否かを判定し、前記収束条件を満たすまで、前記パラメータ更新部による更新を繰り返し行う。
【００５７】
第１及び第２の発明によれば、前記混合音時間周波数展開部により、第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の音源からの音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、各時刻における各周波数の観測時間周波数成分を表す観測スペクトログラムを出力し、前記パラメータ更新部により、前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分を要素に持つテンソルと、前記複数の音源の各々から各マイクまでの伝達特性を表すステアリングベクトルの絶対値成分、前記複数の音源の各々についての前記基底パワースペクトル、前記複数の音源の各々についての、各時刻及び各基底におけるパワーを表すアクティベーションパラメータから求められるテンソルとの距離を表す目的関数を最小化するように、前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータとの更新の各々を行い、前記収束判定部により、予め定められた収束条件を満たしたか否かを判定し、収束条件を満たすまで、前記パラメータ更新部による更新を繰り返し行う。
【００５８】
このように、第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の音源からの音源信号が混合された観測信号の時系列データを入力として、観測スペクトログラムを出力し、第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、観測スペクトログラムを出力し、取得した第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の基底パワースペクトルと、に基づいて、目的関数を最小化するように、複数の音源の各々のステアリングベクトルの絶対値成分と、複数の音源の各々のアクティベーションパラメータとの更新の各々を行うことを繰り返すことにより、雑音源が複数、又は移動する場合であっても雑音を抑圧することができる。
【００５９】
また、第１及び第２の発明において、前記パラメータ更新部は、前記混合音時間周波数展開部により取得した前記第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び前記第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の、各周波数及び各基底におけるパワースペクトルを表す基底パワースペクトルと、に基づいて、前記目的関数を最小化するように、前記複数の音源の各々の前記ステアリングベクトルの絶対値成分と、前記複数の音源の各々の前記アクティベーションパラメータと、前記基底パワースペクトルとの更新の各々を行ってもよい。
【００６０】
また、第１及び第２の発明において、前記第１のマイクを、音声を体表から取得する非可聴つぶやきマイクとしてもよい。
【００６１】
また、第１の発明において、前記距離の距離尺度を板倉齋藤擬距離基準としてもよい。
【００６２】
また、本発明のプログラムは、上記信号解析装置を構成する各部として機能させるためのプログラムである。
【発明の効果】
【００６３】
以上説明したように、本発明の信号解析方法、装置、及びプログラムによれば、第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の音源からの音源信号が混合された観測信号の時系列データを入力として、観測スペクトログラムを出力し、第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、観測スペクトログラムを出力し、取得した第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の基底パワースペクトルと、に基づいて、目的関数を最小化するように、複数の音源の各々のステアリングベクトルの絶対値成分と、複数の音源の各々のアクティベーションパラメータとの更新の各々を行うことを繰り返すことにより、雑音源が複数、又は移動する場合であっても雑音を抑圧することができる、という効果が得られる。
【図面の簡単な説明】
【００６４】
【図１】本実施の形態に係る信号解析装置の機能的な構成例を示すブロック図である。
【図２】本実施の形態に係る２つのマイクの位置関係の一例を示す図である。
【図３】本実施の形態に係る信号解析装置のパラメータ推定処理ルーチンを示すフローチャートである。
【図４】クリーンな体内伝導微弱音声のスペクトログラムの一例を示す図である。
【図５】駅構内の雑音に対する結果の一例を示す図である。
【図６】飲食店内の雑音に対する結果の一例を示す図である。
【図７】２つのマイクの位置関係の一例を示す図である。
【発明を実施するための形態】
【００６５】
以下、図面を参照して、本発明の実施の形態を詳細に説明する。
【００６６】
＜本発明の実施の形態の概要＞
【００６７】
まず、本発明の実施の形態の概要について説明する。本発明の実施の形態は、信号処理の技術分野に属し、体内伝導微音声信号に対し雑音抑圧することを目的とする。
【００６８】
また、本発明の実施の形態は、各音源の時間周波数成分の絶対値成分とステアリングベクトルの絶対値成分とで表される非負値瞬時混合モデルに基づき、空気伝導マイクでモニタリングした外部雑音の信号とＮＡＭマイクで取得した信号の時間周波数成分の絶対値成分から、学習サンプルから事前学習した音声の非負値基底パワースペクトルを用いて、モニタリングした外部雑音の時間周波数成分の絶対値成分を手がかりにステアリングベクトルの絶対値成分と音声の基底パワースペクトルのアクティベーションパラメータを推定する。
【００６９】
＜本実施の形態の原理＞
次に、本実施形態の原理について説明する。
【００７０】
複数の雑音源が存在し、それらが移動する場合、空気伝導マイクでモニタリングした外部雑音の位相情報は参照信号として意味を持たない。そこで、観測信号の振幅情報を重視したモデルについて検討する。
【００７１】
観測複素スペクトログラムを
【００７２】
【数２１】
0006519801.tif 000022
【００７３】
とする。ここで、観測複素スペクトログラムの各要素はｙｊ，ω，τであるものとする。また、各音源の複素スペクトログラムを
【００７４】
【数２２】
0006519801.tif 000023
【００７５】
とすると、時間周波数領域における瞬時混合モデルは、下記（９）式で表される。
【００７６】
【数２３】
0006519801.tif 000024
【００７７】
ここで、ｉは音源インデックスを表す。
【００７８】
【数２４】
0006519801.tif 000025
【００７９】
は音源から各マイクまでの伝達特性を表わす時不変なステアリングベクトルで、下記（１０）式に示すように絶対値と偏角の成分に分解できる。
【００８０】
【数２５】
0006519801.tif 000026
【００８１】
ここで、新たに下記（１１）式、及び（１２）式を定義する。
【００８２】
【数２６】
0006519801.tif 000027
0006519801.tif 000028
【００８３】
このとき、下記の１）〜３）を仮定する。
【００８４】
１）
【数２７】
0006519801.tif 000029
【００８５】
は区間［０，２π）で一様分布に従う、
【００８６】
２）
【数２８】
0006519801.tif 000030
【００８７】
は互いに独立である、
【００８８】
３）
【数２９】
0006519801.tif 000031
【００８９】
は複素正規分布
【００９０】
【数３０】
0006519801.tif 000032
【００９１】
に従う、上記１）〜３）を仮定し、本枠組みにおいて信頼度の低い偏角成分を確率変数とみなして下記（１３）式のように周辺化する。
【００９２】
【数３１】
0006519801.tif 000033
【００９３】
ここで、ｐｉ，ω，τは、各音源のパワースペクトログラムの期待値を表すパラメータである。結果、上記（９）〜（１３）の関係式より、上記（９）式の観測スペクトログラムも同様に、下記（１４）式のような複素正規分布でモデル化される。
【００９４】
【数３２】
0006519801.tif 000034
【００９５】
また、
【００９６】
【数３３】
0006519801.tif 000035
【００９７】
は、｜αｊ，ｉ，ω｜２を要素に持つ対角行列であり、下記（１５）式のように表される。
【００９８】
【数３４】
0006519801.tif 000036
【００９９】
よって、
【０１００】
【数３５】
0006519801.tif 000037
【０１０１】
とおくと、ｙｊ，ω，τが与えられたもとでのｘｊ，ω，τの対数尤度は、下記（１６）式となる。
【０１０２】
【数３６】
0006519801.tif 000038
【０１０３】
この対数尤度は、
【０１０４】
【数３７】
0006519801.tif 000039
【０１０５】
のとき最大となるので、
【０１０６】
【数３８】
0006519801.tif 000040
【０１０７】
の最大化問題は、対数尤度差
【０１０８】
【数３９】
0006519801.tif 000041
【０１０９】
の最小化問題に置き換えられる。また、この対数尤度差は下記（１７）式に示すとおりであり、板倉齋藤擬距離
【０１１０】
【数４０】
0006519801.tif 000042
【０１１１】
と等しいことがわかる。
【０１１２】
【数４１】
0006519801.tif 000043
【０１１３】
したがって、ｐｉ，ω，τに対し、非負値行列因子分解（Non-negative Matrix Factorization: NMF）で用いられる構造
【０１１４】
【数４２】
0006519801.tif 000044
【０１１５】
を仮定すると、パラメータ
【０１１６】
【数４３】
0006519801.tif 000045
、
0006519801.tif 000046
、及び
0006519801.tif 000047
【０１１７】
の最大尤度推定問題は、観測パワースペクトログラム｜ｙｊ，ω，τ｜２を要素に持つテンソル
【０１１８】
【数４４】
0006519801.tif 000048
【０１１９】
に対し、板倉齋藤擬距離基準の非負値テンソル分解（Non-negative Tensor Factorization: NTF）を行うことと等価となる。
【０１２０】
上記（１７）式より、板倉齋藤擬距離基準のＮＴＦの目的関数は、パラメータ
【０１２１】
【数４５】
0006519801.tif 000049
【０１２２】
に依らない項を省略すると下記（１８）式で表される。
【０１２３】
【数４６】
0006519801.tif 000050
【０１２４】
次に、補助関数法により、上記（１８）式の上限関数
【０１２５】
【数４７】
0006519801.tif 000051
【０１２６】
を、下記（１９）式のように設計する。なお、下記（１９）式も上記（１８）式と同様に、パラメータ
【０１２７】
【数４８】
0006519801.tif 000052
【０１２８】
に依らない項を省略している。
【０１２９】
【数４９】
0006519801.tif 000053
【０１３０】
ここで、
【０１３１】
【数５０】
0006519801.tif 000054
【０１３２】
は補助変数であり、それぞれ、下記（２０）式に示されるときに、不等式
【０１３３】
【数５１】
0006519801.tif 000055
【０１３４】
の等号が成立する。
【０１３５】
【数５２】
0006519801.tif 000056
【０１３６】
上記（１９）式の補助関数を用いることで、下記（２１）式、（２２）式、及び（２３）式が導かれる。
【０１３７】
【数５３】
0006519801.tif 000057
【０１３８】
雑音環境下では、空気伝導マイクにわずかに漏れ込んだ（空気伝導）微弱音声を無視できる。したがって、
【０１３９】
【数５４】
0006519801.tif 000058
【０１４０】
を体内伝導微弱音声のスペクトログラムとすると、
【０１４１】
【数５５】
0006519801.tif 000059
【０１４２】
より、下記（２４）式のように固定することができる。
【０１４３】
【数５６】
0006519801.tif 000060
【０１４４】
また、マイク間距離が短いため、雑音源がユーザに極めて近い場合を除き、マイク音源間の距離差によって生じる減衰量の差は無視できる。このことから、各マイクの特性やアンプの特性などが時不変であると仮定すれば、空気伝導マイクで観測した複数の雑音をＮＡＭマイクに混入する外部雑音の参照信号として利用することができる。結果、雑音源の数によらず、音源インデックスは
【０１４５】
【数５７】
0006519801.tif 000061
【０１４６】
と固定できる。なお、マイク間距離は例えば５センチ程度としてもよいが、さらに短くしてもよい。
【０１４７】
＜本発明の実施形態に係る信号解析装置の構成＞
次に、本発明の実施形態に係る信号解析装置の構成について説明する。図１に示すように、本発明の実施形態に係る信号解析装置１００は、ＣＰＵと、ＲＡＭと、後述するパラメータ推定処理ルーチンを実行するためのプログラムや各種データを記憶したＲＯＭと、を含むコンピュータで構成することができる。この信号解析装置１００は、機能的には図１に示すようにＮＡＭマイク１０と、空気伝導マイク１２と、演算部２０と、出力部９０と、を含んで構成されている。なお、第１のマイクの音声を体表から取得する非可聴つぶやきマイクであるＮＡＭマイク１０と、第２のマイクである空気伝導マイク１２の位置関係を図２に示す。図２の場合、ＮＡＭマイク１０は、耳介後方部（乳様突起直下の胸鎖乳突筋上）に設置され、空気伝導マイク１２は、ＮＡＭマイク１０の周辺に設置され、両マイクの間の距離は５センチ程度に設定されている。また、本実施形態においてｊ＝１はＮＡＭマイク１０を表し、ｊ＝２は空気伝導マイク１２を表すものとする。また、ｉ＝１は体内伝導微弱音声の音源を表し、ｉ＝２は雑音源を表す。
【０１４８】
ＮＡＭマイク１０は、空気伝導外部雑音を含まない、クリーンな体内伝導微弱音声の音源信号（以後、クリーン音源信号）の時系列データを出力し、データベース２２に記憶する。また、ＮＡＭマイク１０は、空気伝導外部雑音の音源信号と体内伝導微弱音声の音源信号とが混じっている混合音の音源信号（以後、混合音音源信号）の時系列データを、混合音時間周波数展開部に出力する。
【０１４９】
空気伝導マイク１２は、空気伝導外部雑音の音源信号の時系列データを、混合音時間周波数展開部に出力する。
【０１５０】
演算部２０は、データベース２２と、時間周波数展開部２４と、パラメータ推定部２６と、パラメータ更新部３２と、収束判定部３４と、音声信号生成部４０と、を含んで構成されている。
【０１５１】
データベース２２には、ＮＡＭマイク１０から出力されたクリーン音源信号の時系列データが記憶されている。また、データベース２２には、後述するパラメータ推定部２６により推定された、処理対象のクリーン音源信号に対応したパワースペクトルｗｉ，ω，ｋの各々が記憶されている。また、データベース２２には、後述するパラメータ更新部３２により更新されたパラメータの各々が記憶されている。また、データベース２２には、収束判定部３４から入力される目的関数の値が記憶されている。
【０１５２】
時間周波数展開部２４は、データベース２２に記憶されているクリーン音源信号の時系列データに基づいて、各時刻τにおける各周波数ωの観測時間周波数成分を表す観測複素スペクトログラム
【０１５３】
【数５８】
0006519801.tif 000062
【０１５４】
を計算する。なお、本実施形態においては、短時間フーリエ変換やウェーブレット変換などの時間周波数展開を行う。
【０１５５】
パラメータ推定部２６は、時間周波数展開部２４において取得した観測複素スペクトログラム
【０１５６】
【数５９】
0006519801.tif 000063
【０１５７】
に対して従来技術である、多チャンネルＮＭＦを用いて各音源ｉ、各周波数ω、及び各基底ｋにおけるパワースペクトルｗｉ，ω，ｋの推定を行い、基底パワースペクトル
【０１５８】
【数６０】
0006519801.tif 000064
【０１５９】
を取得する。具体的には、パワースペクトルｗｉ，ω，ｋの推定を行い、結果として得られたパワースペクトルを、処理対象としたクリーン音源信号に対応したパワースペクトルｗｉ，ω，ｋとして用いる。また、処理対象のクリーン音源信号に対応したパワースペクトルｗｉ，ω，ｋを保持する。
【０１６０】
混合音時間周波数展開部２８は、時間周波数展開部２４と同様に、ＮＡＭマイク１０から入力された混合音音源信号の時系列データに基づいて、各時刻τにおける各周波数ωの観測時間周波数成分を表す観測複素スペクトログラム
【０１６１】
【数６１】
0006519801.tif 000065
【０１６２】
を計算する。
【０１６３】
また、混合音時間周波数展開部２８は、時間周波数展開部２４と同様に、空気伝導マイク１２から入力された空気伝導外部雑音の音源信号の時系列データに基づいて、各時刻τにおける各周波数ωの観測時間周波数成分を表す観測複素スペクトログラム
【０１６４】
【数６２】
0006519801.tif 000066
【０１６５】
を計算する。
【０１６６】
パラメータ更新部３２は、上記（１８）式に示す距離尺度を板倉齋藤擬距離基準とするＮＴＦの目的関数を最小化するように、混合音時間周波数展開部２８において取得した混合音音源信号の時系列データに基づく観測複素スペクトログラム
【０１６７】
【数６３】
0006519801.tif 000067
【０１６８】
、及び空気伝導外部雑音の音源信号の時系列データに基づく観測複素スペクトログラム
【０１６９】
【数６４】
0006519801.tif 000068
【０１７０】
と、データベース２２に記憶されている基底パワースペクトル
【０１７１】
【数６５】
0006519801.tif 000069
【０１７２】
に含まれるパワースペクトルｗｉ，ω，ｋと初期値であるパワースペクトルｗｉ，ω，ｋとを含む基底パワースペクトル
【０１７３】
【数６６】
0006519801.tif 000070
【０１７４】
、又は前回更新した基底パワースペクトル
【０１７５】
【数６７】
0006519801.tif 000071
【０１７６】
と、初期値である、又は前回更新した、各時刻τ及び各基底ｋにおける各音源ｉの音源信号のパワーｈｉ，ｋ，τを表すアクティベーションパラメータ
【０１７７】
【数６８】
0006519801.tif 000072
【０１７８】
と、初期値である、又は前回更新した、各マイクｊの各音源iの各周波数ωにおける音源iからマイクｊまでの伝達特性に関するステアリングベクトルの絶対値成分｜αｊ,i,ω｜２と、に基づいて、基底パワースペクトル
【０１７９】
【数６９】
0006519801.tif 000073
【０１８０】
、アクティベーションパラメータ
【０１８１】
【数７０】
0006519801.tif 000074
【０１８２】
、及びステアリングベクトルの絶対値成分｜αｊ,i,ω｜２の更新を行う。
【０１８３】
なお、パラメータ更新部３２は、更新した各々の値を更新した値として、データベース２２に記憶する。また、パラメータ更新部３２は、前回更新した値をデータベース２２から取得するものとする。なお、処理に必要なパラメータの各々についてデータベース２２に値が存在しない場合については、事前にパラメータの各々について初期値を設定し、処理を行うものとする。
【０１８４】
具体的には、混合音時間周波数展開部２８において取得した混合音音源信号の時系列データに基づく観測複素スペクトログラム
【０１８５】
【数７１】
0006519801.tif 000075
【０１８６】
、及び空気伝導外部雑音の音源信号の時系列データに基づく観測複素スペクトログラム
【０１８７】
【数７２】
0006519801.tif 000076
【０１８８】
と、データベース２２に記憶されている基底パワースペクトル
【０１８９】
【数７３】
0006519801.tif 000077
【０１９０】
に含まれるパワースペクトルｗｉ，ω，ｋと初期値であるパワースペクトルｗｉ，ω，ｋとを含む基底パワースペクトル
【０１９１】
【数７４】
0006519801.tif 000078
【０１９２】
、又は前回更新した基底パワースペクトル
【０１９３】
【数７５】
0006519801.tif 000079
【０１９４】
と、初期値である、又は前回更新したアクティベーションパラメータ
【０１９５】
【数７６】
0006519801.tif 000080
【０１９６】
と、初期値である、又は前回更新したステアリングベクトルの絶対値成分｜αｊ,i,ω｜２とに基づいて、上記（２１）式、（２２）式、及び（２３）式の各々に従って、基底パワースペクトル
【０１９７】
【数７７】
0006519801.tif 000081
【０１９８】
、アクティベーションパラメータ
【０１９９】
【数７８】
0006519801.tif 000082
【０２００】
、及びステアリングベクトルの絶対値成分｜αｊ,i,ω｜２の各々を更新し、各パラメータの更新処理毎に更新した各値をデータベースに記憶する。
【０２０１】
また、パラメータ更新部３２は、混合音時間周波数展開部２８において取得した混合音音源信号の時系列データに基づく観測複素スペクトログラム
【０２０２】
【数７９】
0006519801.tif 000083
【０２０３】
、及び空気伝導外部雑音の音源信号の時系列データに基づく観測複素スペクトログラム
【０２０４】
【数８０】
0006519801.tif 000084
【０２０５】
と、更新した基底パワースペクトル
【０２０６】
【数８１】
0006519801.tif 000085
【０２０７】
、アクティベーションパラメータ
【０２０８】
【数８２】
0006519801.tif 000086
【０２０９】
、及びステアリングベクトルの絶対値成分｜αｊ,i,ω｜２の各々の値とに基づいて、上記（１８）式に従って、目的関数の値を算出し、収束判定部３４に出力する。
【０２１０】
なお、上記各パラメータの更新処理の順番は任意であってもよい。
【０２１１】
【数８３】
0006519801.tif 000087
【０２１２】
であり、かつ
【０２１３】
【数８４】
0006519801.tif 000088
【０２１４】
として各パラメータの更新処理を行うものとする。
【０２１５】
収束判定部３４は、パラメータ更新部３２において取得した目的関数の値とデータベース２２に記憶されている前回の目的関数の値との差分が、予め定められた閾値以下である場合に、収束条件を満たすと判定する。収束条件を満たすまで、パラメータ更新部３２における更新処理と、収束判定部３４における収束条件の判定処理を繰り返す。なお、収束条件を満たさない場合には、パラメータ更新部３２において取得した目的関数の値を前回の目的関数の値としてデータベース２２に記憶する。
【０２１６】
音声信号生成部４０は、パラメータ更新部３２において更新した基底パワースペクトル
【０２１７】
【数８５】
0006519801.tif 000089
【０２１８】
と、アクティベーションパラメータ
【０２１９】
【数８６】
0006519801.tif 000090
【０２２０】
とに基づいて、空気伝導外部雑音が除去された分離音を算出し、出力部９０から出力する。
【０２２１】
＜本発明の実施形態に係る信号解析装置の作用＞
次に、本発明の実施形態に係る信号解析装置１００の作用について説明する。まず、ＮＡＭマイク１０から出力された、クリーン音源信号の時系列データを受け付け、データベース２２に記憶する。次に、ＮＡＭマイク１０から出力された、混合音音源信号の時系列データを受け付け、かつ空気伝導マイク１２から出力された、空気伝導外部雑音の時系列データを受け付けると、信号解析装置１００は、図３に示すパラメータ推定処理ルーチンを実行する。
【０２２２】
まず、ステップＳ１００では、データベース２２に記憶されているクリーン音源信号の時系列データを読み込む。
【０２２３】
次に、ステップＳ１０２では、ステップＳ１００において取得したクリーン音源信号の時系列データに基づいて、観測複素スペクトログラム
【０２２４】
【数８７】
0006519801.tif 000091
【０２２５】
を取得する。
【０２２６】
次に、ステップＳ１０４では、ステップＳ１０２において取得した観測複素スペクトログラム
【０２２７】
【数８８】
0006519801.tif 000092
【０２２８】
に基づいて、基底パワースペクトル
【０２２９】
【数８９】
0006519801.tif 000093
【０２３０】
を取得する。
【０２３１】
次に、ステップＳ１０５では、ＮＡＭマイク１０から受け付けた混合音音源信号の時系列データに基づいて、観測複素スペクトログラム
【０２３２】
【数９０】
0006519801.tif 000094
【０２３３】
を取得する。
【０２３４】
次に、ステップＳ１０６では、空気伝導マイク１２から受け付けた空気伝導外部雑音の時系列データに基づいて、観測複素スペクトログラム
【０２３５】
【数９１】
0006519801.tif 000095
【０２３６】
を取得する。
【０２３７】
次に、ステップＳ１０７では、データベース２２に記憶されていない各パラメータの値について初期値を設定する。
【０２３８】
次に、ステップＳ１０８では、ステップＳ１０４及びステップＳ１０７において取得した、又はデータベース２２に記憶されている前回ステップＳ１１０において更新した基底パワースペクトル
【０２３９】
【数９２】
0006519801.tif 000096
【０２４０】
と、ステップＳ１０５において取得した混合音音源信号の観測複素スペクトログラム
【０２４１】
【数９３】
0006519801.tif 000097
【０２４２】
と、ステップＳ１０６において取得した空気伝導外部雑音についての観測複素スペクトログラム
【０２４３】
【数９４】
0006519801.tif 000098
【０２４４】
と、ステップＳ１０７において初期値が設定された、又はデータベース２２に記憶されているステップＳ１１２において前回更新したアクティベーションパラメータ
【０２４５】
【数９５】
0006519801.tif 000099
【０２４６】
と、ステップＳ１０７において初期値が設定された、又はデータベース２２に記憶されているステップＳ１０８において前回更新したステアリングベクトルの絶対値成分｜αｊ,i,ω｜２とに基づいて、上記（２１）式に従って、ステアリングベクトルの絶対値成分｜αｊ,i,ω｜２を更新し、データベース２２に記憶する。
【０２４７】
次に、ステップＳ１１０では、ステップＳ１０４及びステップＳ１０７において取得した、又はデータベース２２に記憶されているステップＳ１１０において前回更新した基底パワースペクトル
【０２４８】
【数９６】
0006519801.tif 000100
【０２４９】
と、ステップＳ１０５において取得した混合音音源信号の観測複素スペクトログラム
【０２５０】
【数９７】
0006519801.tif 000101
【０２５１】
と、ステップＳ１０６において取得した空気伝導外部雑音についての観測複素スペクトログラム
【０２５２】
【数９８】
0006519801.tif 000102
【０２５３】
と、ステップＳ１０７において初期値が設定された、又はデータベース２２に記憶されているステップＳ１１２において前回更新したアクティベーションパラメータ
【０２５４】
【数９９】
0006519801.tif 000103
【０２５５】
と、ステップＳ１０７において初期値が設定された、又はデータベース２２に記憶されているステップＳ１０８において前回更新したステアリングベクトルの絶対値成分｜αｊ,i,ω｜２とに基づいて、上記（２２）式に従って、基底パワースペクトル
【０２５６】
【数１００】
0006519801.tif 000104
【０２５７】
を更新し、データベース２２に記憶する。
【０２５８】
次に、ステップＳ１１２では、ステップＳ１０４及びステップＳ１０７において取得した、又はデータベース２２に記憶されているステップＳ１１０において前回更新した基底パワースペクトル
【０２５９】
【数１０１】
0006519801.tif 000105
【０２６０】
と、ステップＳ１０５において取得した混合音音源信号の観測複素スペクトログラム
【０２６１】
【数１０２】
0006519801.tif 000106
【０２６２】
と、ステップＳ１０６において取得した空気伝導外部雑音についての観測複素スペクトログラム
【０２６３】
【数１０３】
0006519801.tif 000107
【０２６４】
と、ステップＳ１０７において初期値が設定された、又はデータベース２２に記憶されているステップＳ１１２において前回更新したアクティベーションパラメータ
【０２６５】
【数１０４】
0006519801.tif 000108
【０２６６】
と、ステップＳ１０７において初期値が設定された、又はデータベース２２に記憶されているステップＳ１０８において前回更新したステアリングベクトルの絶対値成分｜αｊ,i,ω｜２とに基づいて、上記（２３）式に従って、アクティベーションパラメータ
【０２６７】
【数１０５】
0006519801.tif 000109
【０２６８】
を更新し、データベース２２に記憶する。
【０２６９】
次に、ステップＳ１１３では、ステップＳ１０５において取得した混合音音源信号の観測複素スペクトログラム
【０２７０】
【数１０６】
0006519801.tif 000110
【０２７１】
と、ステップＳ１０６において取得した空気伝導外部雑音についての観測複素スペクトログラム
【０２７２】
【数１０７】
0006519801.tif 000111
【０２７３】
と、ステップＳ１０８において取得したステアリングベクトルの絶対値成分｜αｊ,i,ω｜２と、ステップＳ１１０において取得した基底パワースペクトル
【０２７４】
【数１０８】
0006519801.tif 000112
【０２７５】
と、ステップＳ１１２において取得したアクティベーションパラメータ
【０２７６】
【数１０９】
0006519801.tif 000113
【０２７７】
とに基づいて、上記（１８）式に従って、目的関数の値を算出する。
【０２７８】
次に、ステップＳ１１４では、ステップＳ１１３において取得した目的関数の値と、データベース２２に記憶されている前回の目的関数の値との差分とに基づいて、収束条件を満たすか否かを判定する。収束条件を満たした場合には、パラメータ推定処理ルーチンを終了する。一方、収束条件を満たしていない場合には、ステップＳ１１３において取得した目的関数の値をデータベース２２に記憶し、ステップＳ１０８へ移行し、ステップＳ１０８〜ステップＳ１１４までの処理を繰り返す。
【０２７９】
また、本実施形態に係る信号解析装置１００は、上述のパラメータ推定処理ルーチンにより取得した基底パワースペクトル
【０２８０】
【数１１０】
0006519801.tif 000114
【０２８１】
と、アクティベーションパラメータ
【０２８２】
【数１１１】
0006519801.tif 000115
【０２８３】
とに基づいて、空気伝導外部雑音が除去された分離音を算出し、出力部９０から出力する。
【０２８４】
なお、上述のステップＳ１００〜ステップＳ１０４までの処理を１つの処理ルーチンとし、ステップＳ１０５〜ステップＳ１１４までの処理を１つの処理ルーチンとして別々に行ってもよい。この場合、ＮＡＭマイク１０からクリーン音源信号の時系列データを受け付けると、信号解析装置１００は、ステップＳ１００〜ステップＳ１０４までの処理の処理ルーチンを実行する。一方、ＮＡＭマイク１０から、混合音音源信号の時系列データを受け付け、かつ空気伝導マイク１２から、空気伝導外部雑音の時系列データを受け付けると、信号解析装置１００は、ステップＳ１０５〜ステップＳ１１４までの処理の処理ルーチンを実行する。
【０２８５】
＜実験例＞
本実施形態に係る信号解析装置１００を用いて、実環境下で収録した雑音に対する雑音抑圧性能を検証した。雑音は駅構内の雑音と、飲食店内の雑音の２種類とし、それぞれＳＮ比０ｄＢでクリーンな体内伝導微弱音声および空気伝導（漏れ込み）微弱音声に重畳する。
【０２８６】
なお、収録信号のサンプリング周波数は１６ｋＨｚ、ＦＦＴ分析のフレーム長は６４ｍｓ、シフト長は３２ｍｓとする。ＢＳＳの更新回数は２００回、ＮＴＦおよびＭＮＭＦＤの更新回数は５０回、各音源の基底数は２０、ＮＭＦＤ基底のタップ数は１０＊４とする。
【０２８７】
図４に目的信号であるクリーンな体内伝導微弱音声のスペクトログラムを示す。また、図５に駅構内の雑音に対する雑音抑圧結果の例を示す。また、図６に飲食店内の雑音に対する雑音抑圧結果の例を示す。
【０２８８】
以上説明したように、本発明の実施形態に係る信号解析装置によれば、第１のマイクで取得した体内伝導微弱音声及び外部雑音を含む複数の音源からの音源信号が混合された観測信号の時系列データを入力として、観測スペクトログラムを出力し、第２のマイクで取得した前記外部雑音の音源信号が混合された観測信号の時系列データを入力として、観測スペクトログラムを出力し、取得した第１のマイクで取得した音源信号の観測スペクトログラムの絶対値成分、及び第２のマイクで取得した音源信号の観測スペクトログラムの絶対値成分と、予め学習された、体内伝導微弱音声の基底パワースペクトルと、に基づいて、目的関数を最小化するように、複数の音源の各々のステアリングベクトルの絶対値成分と、複数の音源の各々のアクティベーションパラメータとの更新の各々を行い、予め定められた収束条件を満たすまで、各パラメータの更新を繰り返し行うことにより、雑音源が複数、又は移動する場合であっても雑音を抑圧することができる。
【０２８９】
なお、本発明は、上述した実施形態に限定されるものではなく、この発明の要旨を逸脱しない範囲内で様々な変形や応用が可能である。
【０２９０】
また、上述の信号解析装置１００は、内部にコンピュータシステムを有しているが、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
【０２９１】
また、本願明細書中において、プログラムが予めインストールされている実施形態として説明したが、当該プログラムを、コンピュータ読み取り可能な記録媒体、例えばＣＤ−ＲＯＭやメモリーカード等に格納して提供することも可能である。
【符号の説明】
【０２９２】
１０   ＮＡＭマイク
１２   空気伝導マイク
２０   演算部
２２   データベース
２４   時間周波数展開部
２６   パラメータ推定部
２８   混合音時間周波数展開部
３２   パラメータ更新部
３４   収束判定部
４０   音声信号生成部
９０   出力部
１００ 信号解析装置
【図１】
0006519801.tif 000123
【図２】
0006519801.tif 000124
【図３】
0006519801.tif 000125
【図４】
0006519801.tif 000126
【図５】
0006519801.tif 000127
【図６】
0006519801.tif 000128
【図７】
0006519801.tif 000129