(57)【特許請求の範囲】
【請求項１】
音響信号データが記憶された信号データ記憶手段と、
前記信号データ記憶手段から読み込んだ前記音響信号データに対する時間周波数解析により、前記音響信号データの時間周波数成分を格納したデータ行列Ｙを求める時間周波数解析手段と、
非負値行列因子分解手法を用いて、前記データ行列Ｙを、基底行列Ｈと係数行列Ｕの積で近似した場合に、求めるべきスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値それぞれの初期値を設定する初期値設定手段と、
前記初期値が設定されたスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値を用いて、前記係数行列Ｕの各行に対応するスペクトル基底アクティビティ時系列を畳み込み混合の形式で表現したスペクトログラムモデルの値を算出するモデル算出手段と、
前記スペクトログラムモデルの値と、前記データ行列Ｙと、前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値それぞれの値を用いて、前記初期値を設定した前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータと、前記パワーエンベロープ基底アクティビティ値とを更新する更新手段と、
前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値が収束するまで前記更新手段により値の更新を続け、収束した時点で、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ、前記パワーエンベロープ基底アクティビティ値及び前記スペクトログラムモデルの値を出力する出力手段と、
前記パワーエンベロープ基底アクティビティ値の時系列から各ノートのオンセットに対応するピークを検出し、各ピーク近傍の値と前記パワーエンベロープ基底パラメータと前記スペクトル基底パラメータとを用いて、楽譜における音符に相当するノート毎の音響信号のスペクトログラムを算出するスペクトログラム算出手段と
を備えたことを特徴とする信号分析装置。
【請求項２】
音響信号データが記憶された信号データ記憶手段を備えた信号分析装置のコンピュータに、信号分析処理を行わせる信号分析方法であって、
前記信号データ記憶手段から読み込んだ前記音響信号データに対する時間周波数解析により、前記音響信号データの時間周波数成分を格納したデータ行列Ｙを求める時間周波数解析ステップと、
非負値行列因子分解手法を用いて、前記データ行列Ｙを、基底行列Ｈと係数行列Ｕの積で近似した場合に、求めるべきスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値それぞれの初期値を設定する初期値設定ステップと、
前記初期値が設定されたスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値を用いて、前記係数行列Ｕの各行に対応するスペクトル基底アクティビティ時系列を畳み込み混合の形式で表現したスペクトログラムモデルの値を算出するモデル算出ステップと、
前記スペクトログラムモデルの値と、前記データ行列Ｙと、前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値それぞれの値を用いて、前記初期値を設定した前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータと、前記パワーエンベロープ基底アクティビティ値とを更新する更新ステップと、
前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値が収束するまで前記更新ステップにより値の更新を続け、収束した時点で、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ、前記パワーエンベロープ基底アクティビティ値及び前記スペクトログラムモデルの値を出力する出力ステップと、
前記パワーエンベロープ基底アクティビティ値の時系列から各ノートのオンセットに対応するピークを検出し、各ピーク近傍の値と前記パワーエンベロープ基底パラメータと前記スペクトル基底パラメータとを用いて、楽譜における音符に相当するノート毎の音響信号のスペクトログラムを算出するスペクトログラム算出ステップと
を有することを特徴とする信号分析方法。
【請求項３】
コンピュータを、請求項１に記載の信号分析装置として機能させるための信号分析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音響信号をノート（楽譜における音符）ごとの信号に分解するために音響信号を分析する信号分析装置、信号分析方法及び信号分析プログラムに関する。
【背景技術】
【０００２】
複数の音響信号が重畳する混合信号から個々の音響信号を分離抽出することは容易ではない。このような問題は音源分離と呼ばれ、特にモノラル信号を対象とした信号分離は典型的な不良設定問題であり、何の仮定も置かずに解くことは困難である。モノラル信号分離の研究は、これまで多くのアプローチが検討されているが、近年有効なアプローチとして脚光を浴びているのは、非負値行列因子分解（Non-negative Matrix Factorization；ＮＭＦ）の原理を応用した手法である（例えば、特許文献１参照）。このアプローチでは、観測信号の各時刻でのスペクトル（周波数成分の大きさ）を列ベクトルとして並べた非負値のデータ行列Ｙを非負値の基底行列Ｈと非負値の係数行列Ｕの積の形で近似する。
【数１】
0005580585.tif 000002
【０００３】
この結果、観測される全スペクトルを構成しているスペクトル基底関数が基底行列Ｈの各列に格納され、特定のスペクトル基底関数が各時刻でどの程度の大きさでアクティベートしているかを表すスペクトル基底アクティビティ値の時系列が係数行列Ｕの一つの行に格納される。以上のようにして、信号の分解表現を得ることができる。この手法は、音源分離の問題に対して、「観測信号は限られた種類のスペクトルをもった音だけで構成される」という仮定を活用している点が特徴的であり、この仮定に適合する信号に対しては有効な解法となる。
【０００４】
非負値行列因子分解によるスペクトログラムの分解表現は、スペクトログラムを
【数２】
0005580585.tif 000003
と見立て、
【数３】
0005580585.tif 000004
となるような
【数４】
0005580585.tif 000005
【数５】
0005580585.tif 000006
を決めてやることによって得られる。ただし、ω、ｔはそれぞれ周波数と時刻に対応するインデックスである。
【０００５】
ｙｔ：＝（Ｙ１，，・・・，ＹΩ，ｔ）Ｔ，ｈｉ：＝（Ｈ１，ｉ，・・・，ＨΩ，ｉ）Ｔとすると式（２）は、
【数６】
0005580585.tif 000007
と書き直せることから分かるとおり、すべてのｔにおける観測データｙｔを、高々Ｉ種類の「パーツ」ｈ１，・・・，ｈＩだけで構成されたものと見なそうとしていることになり、そのためにそれぞれのパーツをどのように置くのが最も妥当かを判断する問題と理解される。これにより得られるそれぞれの行列が表すものは、図２を見るとより分かりやすい。Ｈの各列ベクトルには楽曲中に繰り返し現れるスペクトルが典型的な構成パーツと見なされて表出される。従って、音楽スペクトログラムが楽器の種類や音階で決まる限られたパターンのスペクトルだけで構成されているとするなら、Ｈの各列ベクトルはおおよそ特定楽器の特定音階に対応したスペクトルとなる。一方、Ｕの各行ベクトルはそれぞれのスペクトルパーツがどの時刻にどの程度の強さでアクティベートしているかを表している。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】P. Smaragdis and J. C. Brown, “Non-negative matrix factorization for music transcription," in Proceedings of the IEEEWorkshop on Applications of Signal Processing to Audio and Acoustics (WASPAA), 2003, pp. 177.180.
【発明の概要】
【発明が解決しようとする課題】
【０００７】
前述のとおり、非負値行列因子分解は、観測スペクトルを格納したデータ行列を、スペクトル基底を各列に格納した基底行列とスペクトル基底アクティビティの時系列を各行に格納した係数行列の積で近似することで、観測スペクトルの集合からスペクトル基底関数を自動獲得し、スペクトル基底関数ごとに観測スペクトルを分解する機能をもつものである。
【０００８】
しかしながら、音楽信号をノート（楽譜における音符）ごとの信号に分解するためには、スペクトル基底アクティビティ時系列自体をも各ノートイベントに該当するスペクトル基底アクティビティ時系列に分解する必要がある。明らかに従来の非負値行列因子分解にはこのような階層的な分解表現を取得する機能がないという問題がある。
【０００９】
本発明は、このような事情に鑑みてなされたもので、前述したような階層的な分解表現を得ることができる信号分析装置、信号分析方法及び信号分析プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明は、音響信号データが記憶された信号データ記憶手段と、前記信号データ記憶手段から読み込んだ前記音響信号データに対する時間周波数解析により、前記音響信号データの時間周波数成分を格納したデータ行列Ｙを求める時間周波数解析手段と、非負値行列因子分解手法を用いて、前記データ行列Ｙを、基底行列Ｈと係数行列Ｕの積で近似した場合に、求めるべきスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値それぞれの初期値を設定する初期値設定手段と、前記初期値が設定されたスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値を用いて、前記係数行列Ｕの各行に対応するスペクトル基底アクティビティ時系列を畳み込み混合の形式で表現したスペクトログラムモデルの値を算出するモデル算出手段と、前記スペクトログラムモデルの値と、前記データ行列と、前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値それぞれの値を用いて、前記初期値を設定した前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータと、前記パワーエンベロープ基底アクティビティ値とを更新する更新手段と、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値が収束するまで前記更新手段により値の更新を続け、収束した時点で、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ、前記パワーエンベロープ基底アクティビティ値及び前記スペクトログラムモデルの値を出力する出力手段とを備えたことを特徴とする。
【００１１】
本発明は、音響信号データが記憶された信号データ記憶手段を備えた信号分析装置のコンピュータに、信号分析処理を行わせる信号分析方法であって、前記信号データ記憶手段から読み込んだ前記音響信号データに対する時間周波数解析により、前記音響信号データの時間周波数成分を格納したデータ行列Ｙを求める時間周波数解析ステップと、非負値行列因子分解手法を用いて、前記データ行列Ｙを、基底行列Ｈと係数行列Ｕの積で近似した場合に、求めるべきスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値それぞれの初期値を設定する初期値設定ステップと、前記初期値が設定されたスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値を用いて、前記係数行列Ｕの各行に対応するスペクトル基底アクティビティ時系列を畳み込み混合の形式で表現したスペクトログラムモデルの値を算出するモデル算出ステップと、前記スペクトログラムモデルの値と、前記データ行列と、前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値それぞれの値を用いて、前記初期値を設定した前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータと、前記パワーエンベロープ基底アクティビティ値とを更新する更新ステップと、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値が収束するまで前記更新ステップにより値の更新を続け、収束した時点で、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ、前記パワーエンベロープ基底アクティビティ値及び前記スペクトログラムモデルの値を出力する出力ステップとを有することを特徴とする。
【００１２】
本発明は、音響信号データが記憶された信号データ記憶手段を備えた信号分析装置のコンピュータに、信号分析処理を行わせる信号分析プログラムであって、前記信号データ記憶手段から読み込んだ前記音響信号データに対する時間周波数解析により、前記音響信号データの時間周波数成分を格納したデータ行列Ｙを求める時間周波数解析ステップと、非負値行列因子分解手法を用いて、前記データ行列Ｙを、基底行列Ｈと係数行列Ｕの積で近似した場合に、求めるべきスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値それぞれの初期値を設定する初期値設定ステップと、前記初期値が設定されたスペクトル基底パラメータ、パワーエンベロープ基底パラメータ及びパワーエンベロープ基底アクティビティ値を用いて、前記係数行列Ｕの各行に対応するスペクトル基底アクティビティ時系列を畳み込み混合の形式で表現したスペクトログラムモデルの値を算出するモデル算出ステップと、前記スペクトログラムモデルの値と、前記データ行列と、前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値それぞれの値を用いて、前記初期値を設定した前記スペクトル基底パラメータと、前記パワーエンベロープ基底パラメータと、前記パワーエンベロープ基底アクティビティ値とを更新する更新ステップと、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ及び前記パワーエンベロープ基底アクティビティ値が収束するまで前記更新ステップにより値の更新を続け、収束した時点で、前記スペクトル基底パラメータ、前記パワーエンベロープ基底パラメータ、前記パワーエンベロープ基底アクティビティ値及び前記スペクトログラムモデルの値を出力する出力ステップとを前記コンピュータに行わせることを特徴とする。
【発明の効果】
【００１３】
本発明によれば、複数の音が混在している音響信号からの特定音の検出、複数の音が混在している音響信号からの特定音の抽出、および複数の音が混在している音響信号からの特定音の加工などに信号分析結果のパラメータを利用することができるという効果が得られる。
【図面の簡単な説明】
【００１４】
【図１】本発明の一実施形態の構成を示すブロック図である。
【図２】非負値行列因子分解（ＮＭＦ）による音楽スペクトルグラムの分解表現を示す説明図である。
【発明を実施するための形態】
【００１５】
以下、図面を参照して、本発明の一実施形態による信号分析装置を説明する。初めに、本発明による信号分析装置の基本原理について説明する。音楽信号からデータ行列を取得した場合、非負値行列因子分解によって得られる係数行列の各行に格納されたスペクトル基底アクティビティ時系列は、限られた種類の局所パターンの組み合わせで表現されることが多い。これは、音楽においては、音価と呼ぶノートの発音時間長に関する単位があり、同じ音価のノートは音の立ち上がり方と減衰の仕方が類似することが多いからである。
【００１６】
本発明では、前述した音楽の性質を利用し、観測スペクトルの集合からスペクトル基底関数を自動獲得しつつスペクトル基底関数ごとに観測スペクトルを分解する非負値行列因子分解の従来の機能に加え、各スペクトル基底アクティビティ時系列の中に混在する局所パターンの基底関数（以後、これをパワーエンベロープ基底という）を自動獲得しつつパワーエンベロープ基底ごとにスペクトル基底アクティビティ時系列を分解する機能を実現している。
【００１７】
具体的に説明すると、係数行列Ｕの各行に対応するスペクトル基底アクティビティ時系列を畳み込み混合の形で表現する。すなわち、Ｕのｉ行目の要素をＵｉ，０，Ｕｉ，１，・・・とすると、
【数７】
0005580585.tif 000008
と表す。Ｇｊ，０，Ｇｊ，１，・・・はｊ番目のパワーエンベロープ基底であり、Ｏｉ，ｊ，０，Ｏｉ，ｊ，１，・・・はそのアクティビティを表すので以後、Ｏｉ，ｊ，ｔをパワーエンベロープ基底アクティビティと呼ぶ。ところで、畳み込み混合とは、複数の信号源から到来する信号を複数のマイクロホンで観測する際の、信号源の混合過程のモデルとして広く認知されているものであり、前述の場合、Ｕのｉ行目のスペクトル基底アクティビティ時系列がｉ番目のマイクロホンにおける観測信号に対応し、Ｇｊ，１，Ｇｊ，２，・・・がｊ番目の信号源の信号に対応し、Ｏｉ，ｊ，１，Ｏｉ，ｊ，２・・・がｊ番目の信号源からｉ番目のマイクロホンまでのインパルス応答に対応していることになる。
【００１８】
次に、スペクトログラムモデルについて説明する。図２に示すＵの各行ベクトルを見ると、ｈｉのスペクトル基底アクティベーションの時間エンベロープは、限られた種類の局所パターンだけで構成されていることに気づく。これは、個々のノートは限られた種類の音長で弾かれるため、立ち上がり・減衰パターンの種類もまた限られるからである。そこで、Ｕに関して、さらに、
【数８】
0005580585.tif 000009
のような畳み込み混合の形で表される分解表現を考える。Ｇｊ，τはｊ番目の時間エンベロープの局所パターンを表す。一方Ｏｉ，ｊ，ｔはそのアクティベーションを表し、理想的には各ノートのオンセット時刻にピークが立つようなイメージとなる。式（２）の右辺に式（５）を代入すると、
【数９】
0005580585.tif 000010
のように、式（２）で与えられるＮＭＦ型のモデルを拡張したスペクトログラムモデルが立てられる。ここで、分解のスケールの任意性を除く目的で、
【数１０】
0005580585.tif 000011
を仮定しておく。
【００１９】
次に、最適化アルゴリズムについて説明する。非負値行列因子分解でＹから図２のようにＨとＵを求めたように、提案するモデルのもとでＹからＨとＧとＯを求めたい。以下では、提案モデルによる音楽の階層的スパース表現を得るための最適化アルゴリズムについて説明する。まず、二乗誤差規準について説明する。ここでは、観測スペクトログラムＹ：＝（Ｙω，ｔ）Ω×Ｔのもとで設定される次の最適化問題
【数１１】
0005580585.tif 000012
を検討する。Ｓ（Ｇ，Ｏ）はＧ，Ｏをスパースな解へと誘導する正則化項であり、ここでは
【数１２】
0005580585.tif 000013
と定義する。ただし、０＜ｐｇ≦２，０＜ｐｏ≦２とする。
【００２０】
まず、Ｆを降下させるＨの更新式を導出する。Ｈ，Ｇ，Ｏの１ステップ前での更新値をそれぞれＨ’，Ｇ’，Ｏ’とすると、
【数１３】
0005580585.tif 000014
となる。ただし、
【数１４】
0005580585.tif 000015
である。詳細な説明は省略するが、ＦＨ（Ｈ，Ｇ’，Ｏ’）を最小化するＨは解析的に
【数１５】
0005580585.tif 000016
と求まり、このように更新すればＦ（Ｈ，Ｇ’，Ｏ’）は増加しないことが保証される。また、Ｈ’およびＵ’がともに非負値であればＨも必ず非負値となる。
【００２１】
次に、Ｆを降下させるＧの更新式を導出する。先と同様に、
【数１６】
0005580585.tif 000017
となる。ただし、
【数１７】
0005580585.tif 000018
である。ＦＧ（Ｈ’，Ｇ，Ｏ’）の第２項は
【数１８】
0005580585.tif 000019
の不等式による（右辺は接点±ｘ’で｜ｘ｜ｐに接する放物線であるため明らか）。以上より、ＦＧ（Ｈ’，Ｇ，Ｏ’）を用いてＧの更新式が
【数１９】
0005580585.tif 000020
と導ける。
【００２２】
最後に、Ｆを降下させるＯの更新式を導出する。先と同様に、
【数２０】
0005580585.tif 000021
が言え、これを用いてＯの更新式が
【数２１】
0005580585.tif 000022
と導ける。
【００２３】
次に、Ｉダイバージェンス規準について説明する。モデル化誤差をＩダイバージェンスで測った場合の最適化問題
【数２２】
0005580585.tif 000023
についても検討する。
【００２４】
まず、Ｈの更新式については、
【数２３】
0005580585.tif 000024
の不等式より導かれる
【数２４】
0005580585.tif 000025
を最小化するＨを求めれば良く、
【数２５】
0005580585.tif 000026
のように解析的に得られる。以下同様にして、Ｇの更新式については、
【数２６】
0005580585.tif 000027
の不等式、Ｏの更新式については、
【数２７】
0005580585.tif 000028
の不等式より導かれる
【数２８】
0005580585.tif 000029
をそれぞれ最小化するＧ，Ｏを求めれば良く、
【数２９】
0005580585.tif 000030
のように解析的に得られる。式（２１）の第２項および式（２２）の第３項は
∀ｘ＞０，｜ｘ｜ｐ≦ｐ｜ｘ’｜ｐ−１（ｘ−ｘ’）＋｜ｘ’｜ｐ（０＜ｐ≦２）・・・（２５）による（右辺は接点ｘ’における｜ｘ｜ｐの接線であるため明らか）。
【００２５】
次に、前述した基本原理を使用した信号分析装置について説明する。図１は第１の実施形態における信号分析装置の構成を示すブロック図である。信号分析装置は、コンピュータ装置で構成する。この図において、符号１は、音響信号を標本化・量子化することにより得られる音響信号データを入力して記憶した信号データ記憶部である。符号２は、時間周波数解析を行う時間周波数解析部である。符号３は、値の初期設定を行う初期設定部である。符号４は、スペクトログラムモデルを算出するスペクトログラムモデル算出部である。符号５は、スペクトル基底を更新するスペクトル基底更新部である。符号６は、パワーエンベロープ基底を更新するパワーエンベロープ基底更新部である。符号７は、パワーエンベロープ基底アクティビティを更新するパワーエンベロープ基底アクティビティ更新部である。符号８は、パラメータを規格化するパラメータ規格化部である。符号９は、処理が収束したか否かを判定する収束判定部である。符号１０は、パラメータを出力するパラメータ出力部である。符号１１は、出力されたパラメータを記憶するパラメータ記憶部である。
【００２６】
次に、図１を参照して、図１に示す信号分析装置の動作を説明する。まず、時間周波数解析部２は、信号データ記憶部１に記憶されている分析対象の信号データを読み込み、短時間フーリエ変換（Short-Time Fourier Transform；ＳＴＦＴ）やウェーブレット変換などを用いて、時間周波数解析を行って非負値で与えられる時間周波数成分｛Ｙω，ｔ｝０≦Ω≦Ω−１，０≦ｔ≦Ｔ−１を計算する。ただし、ω＝０，・・・，Ω−１，ｔ＝０，・・・，Ｔ−１はそれぞれ周波数、時刻に対応するインデックスとする。時間周波数解析部２は、時間周波数成分Ｙω，ｔを格納した行列Ｙ＝（Ｙω，ｔ）Ω×Ｔを出力する。
【００２７】
次に、初期設定部３は、スペクトル基底数Ｉ、パワーエンベロープ基底数Ｊおよび正則化パラメータλｇ，λｏ，ｐｇ，ｐｏを決定する。そして、初期設定部３は、非負値行列因子分解（ＮＭＦ）によりＹに対し、
【数３０】
0005580585.tif 000031
となる
【数３１】
0005580585.tif 000032
を出力する。これを用いて、スペクトル基底パラメータＨω，ｉ、パワーエンベロープ基底パラメータＧｊ，ｔ、パワーエンベロープ基底アクティビティ値Ｏｉ，ｊ，ｔの初期値をそれぞれ
【数３２】
0005580585.tif 000033
として出力する。ただし、ｉ＝１，・・・，Ｉはスペクトル基底のインデックス、ｊ＝１，・・・，Ｊはパワーエンベロープ基底のインデックスとする。また、［・］ａ，ｂは行列のａ行ｂ列の成分を表す。
【００２８】
次に、スペクトログラムモデル算出部４は、前段で得られたＨω，ｉ、Ｇｊ，ｔ、Ｏｉ，ｊ，ｔを用いてスペクトログラムモデルＸω，ｔを以下の手順で算出して出力する。まず、スペクトログラムモデル算出部４は、Ｇｊ，ｔとＯｉ，ｊ，ｔを用いてスペクトル基底アクティビティ値Ｕｉ，ｔを畳み込み混合演算
【数３３】
0005580585.tif 000034
により算出する。この畳み込み混合演算は高速フーリエ変換を用いて高速計算する。
【００２９】
次に、スペクトログラムモデル算出部４は、Ｈω，ｉおよび先に求まったＵｉ，ｔを用いて、Ｘω，ｔを積和演算
【数３４】
0005580585.tif 000035
により算出する。
【００３０】
次に、スペクトル基底更新部５は、Ｙω，ｔおよび前段で得られたＸω，ｔとＵｉ，ｔとＨω，ｉを用いて、Ｈω，ｉを
【数３５】
0005580585.tif 000036
により更新して出力する。
【００３１】
次に、パワーエンベロープ基底更新部６は、Ｙω，ｔおよび前段で得られたＸω，ｔとＨω，ｉとＯｉ，ｊ，ｔを用いて、Ｇｊ，ｔを以下の手順により算出する。まず、パワーエンベロープ基底更新部６は、ｉ番目のスペクトル基底Ｈ０，ｉ，・・・，ＨΩ−１，ｉと時刻ｔにおける観測スペクトルＹ０，ｔ，・・・，ＹΩ−１，ｔとの
【数３６】
0005580585.tif 000037
を
【数３７】
0005580585.tif 000038
により算出する。同様に、パワーエンベロープ基底更新部６は、Ｈ０，ｉ，・・・，ＨΩ−１，ｉとＸ０，ｔ，・・・，ＸΩ−１，ｔの
【数３８】
0005580585.tif 000039
を
【数３９】
0005580585.tif 000040
により算出する。
【００３２】
次に、パワーエンベロープ基底更新部６は、
【数４０】
0005580585.tif 000041
を
【数４１】
0005580585.tif 000042
により算出する。この相互相関の演算は、高速フーリエ変換を用いることで高速計算する。同様に、パワーエンベロープ基底更新部６は、
【数４２】
0005580585.tif 000043
を
【数４３】
0005580585.tif 000044
により算出する。この相互相関の演算も高速フーリエ変換を用いることで高速計算する。
【００３３】
最後に、パワーエンベロープ基底更新部６は、
【数４４】
0005580585.tif 000045
および前段で得られたＧｊ，ｔを用いて、Ｇｊ，ｔを
【数４５】
0005580585.tif 000046
により更新する。このλｇｐｇ｜Ｇｊ，τ｜ｐｇ−１はスパース正則化項に関係する項であり、Ｇｊ，ｔの要素をスパースにするように誘導する効果をもつものであればこれ以外の形のもので代替してもよい。
【００３４】
次に、パワーエンベロープ基底アクティビティ更新部７は、Ｙω，ｔおよび前段で得られたＸω，ｔとＨω，ｉとＯｉ，ｊ，ｔを用いて、Ｇｊ，ｔを以下の手順により算出して出力する。まず、パワーエンベロープ基底アクティビティ更新部７は、ｉ番目のスペクトル基底Ｈ０，ｉ，・・・，ＨΩ−１，ｉと時刻ｔにおける観測スペクトルＹ０，ｔ，・・・，ＹΩ−１，ｔとの
【数４６】
0005580585.tif 000047
を
【数４７】
0005580585.tif 000048
により算出する。同様に、パワーエンベロープ基底アクティビティ更新部７は、Ｈ０，ｉ，・・・，ＨΩ−１，ｉとＸ０，ｔ，・・・，ＸΩ−１，ｔの
【数４８】
0005580585.tif 000049
を
【数４９】
0005580585.tif 000050
により算出する。
【００３５】
次に、パワーエンベロープ基底アクティビティ更新部７は、
【数５０】
0005580585.tif 000051
を
【数５１】
0005580585.tif 000052
により算出する。この相互相関の演算は，高速フーリエ変換を用いることで高速計算する。同様に、パワーエンベロープ基底アクティビティ更新部７は、
【数５２】
0005580585.tif 000053
を
【数５３】
0005580585.tif 000054
により算出する。この相互相関の演算も高速フーリエ変換を用いることで高速計算する。
【００３６】
最後に、パワーエンベロープ基底アクティビティ更新部７は、
【数５４】
0005580585.tif 000055
を
【数５５】
0005580585.tif 000056
により更新する。この＋λｏｐｏ｜Ｏ’ｉ，ｊ，τ｜ｐｏ−１はスパース正則化項に関係する項であり、Ｏｉ，ｊ，ｔの要素をスパースにするように誘導する効果をもつものであればこれ以外の形のもので代替してもよい。
【００３７】
次に、パラメータ規格化部８は、前段で得られたＨω，ｉおよびＧｊ，ｔを規格化して出力する。例えば、足して１になるように両者を規格化する場合は、
【数５６】
0005580585.tif 000057
によりＨω，ｉおよびＧｊ，ｔをそれぞれ更新する。
【００３８】
次に、収束判定部９は、前段の反復計算が所定の回数を満たしたか否か、あるいは、反復計算においてパラメータの更新の変化率が所定値以下になったか否か、あるいは、目的関数値の変化率が所定値以下になったか否かを判定する。例えば、目的関数は、
【数５７】
0005580585.tif 000058
により計算する。ただし、Ｓ（Ｇ，Ｏ）はＧ，Ｏをスパースな解へと誘導する正則化項であり、
【数５８】
0005580585.tif 000059
と定義する。収束判定部９は、反復計算が収束していなければ、スペクトログラムモデル算出部４に対して、再度スペクトログラムモデルを算出し直す指示を出力し、これを受けて、スペクトログラムモデル算出部４、スペクトル基底更新部５、パワーエンベロープ基底更新部６、パワーエンベロープ基底アクティビティ更新部７及びパラメータ規格化部８は、反復計算が収束するまで前述した処理動作を繰り返す。
【００３９】
次に、パラメータ出力部９は、反復計算が収束したと見なされたＨω，ｉ，Ｇｊ，ｔ，Ｏｉ，ｊ，ｔ，Ｘω，ｔなどのパラメータをパラメータ記憶部１１に記憶する。
【００４０】
次に、第２の実施形態における信号分析装置を説明する。まず、第２の実施形態におけるスペクトル基底更新部５の処理動作を説明する。スペクトル基底更新部５は、Ｙω，ｔおよび前段で得られたＸω，ｔとＨω，ｉとＵｉ，ｔを用いて、以下の手順によりＨω，ｉを更新する。まず、観測スペクトログラムＹω，ｔとスペクトログラムモデルＸω，ｔとのスペクトログラム比Ｒω，ｔを
【数５９】
0005580585.tif 000060
により算出する。
【００４１】
次に、スペクトル基底更新部５は、先に求まったＲω，ｔおよびＨω，ｉとＵｉ，ｔを用いて、
【数６０】
0005580585.tif 000061
により算出する。
【００４２】
次に、第２の実施形態におけるパワーエンベロープ基底更新部６の処理動作を説明する。パワーエンベロープ基底更新部６、Ｙω，ｔおよび前段で得られたＸω，ｔとＨω，ｉとＧｊ，ｔとＯｉ，ｊ，ｔを用いて、以下の手順によりＧｊ，ｔを更新する。まず、パワーエンベロープ基底更新部６は、観測スペクトログラムＹω，ｔとスペクトログラムモデルＸω，ｔとのスペクトログラム比Ｒω，ｔを
【数６１】
0005580585.tif 000062
により算出する。次に、パワーエンベロープ基底更新部６は、Ｈ０，ｉ，・・・，ＨΩ−１，ｉとＲ０，ｔ，・・・，ＲΩ−１，ｔとの
【数６２】
0005580585.tif 000063
を
【数６３】
0005580585.tif 000064
により計算する。次に、パワーエンベロープ基底更新部６は、
【数６４】
0005580585.tif 000065
【数６５】
0005580585.tif 000066
により計算する。この相互相関の演算は、高速フーリエ変換を用いることで高速計算する。
【００４３】
最後に、パワーエンベロープ基底更新部６は、
【数６６】
0005580585.tif 000067
によりＧｊ，ｔを更新する。この２λｇｐｇ｜Ｇｊ，τ｜ｐｇ−１はスパース正則化項に関係する項であり、Ｇｊ，ｔの要素をスパースにするように誘導する効果をもつものであればこれ以外の形のもので代替してもよい。
【００４４】
次に、第２の実施形態におけるパワーエンベロープ基底アクティビティ更新部７の処理動作を説明する。パワーエンベロープ基底アクティビティ更新部７は、Ｙω，ｔおよび前段で得られたＸω，ｔとＨω，ｉとＧｊ，ｔとＯｉ，ｊ，ｔを用いて、以下の手順によりＯｉ，ｊ，ｔを更新する。まず、パワーエンベロープ基底アクティビティ更新部７は、観測スペクトログラムＹω，ｔとスペクトログラムモデルＸω，ｔとのスペクトログラム比Ｒω，ｔを
【数６７】
0005580585.tif 000068
により算出する。次に、パワーエンベロープ基底アクティビティ更新部７は、Ｈ０，ｉ，・・・，ＨΩ−１，ｉとＲ０，ｔ，・・・，ＲΩ−１，ｔとの
【数６８】
0005580585.tif 000069
を
【数６９】
0005580585.tif 000070
により計算する。次に、パワーエンベロープ基底アクティビティ更新部７は、
【数７０】
0005580585.tif 000071
を
【数７１】
0005580585.tif 000072
により計算する。この相互相関の演算は、高速フーリエ変換を用いることで高速計算する。
【００４５】
最後に、パワーエンベロープ基底アクティビティ更新部７は、
【数７２】
0005580585.tif 000073
によりＯｉ，ｊ，ｔを更新する。この２λｏｐｏ｜Ｏｉ，ｊ，τ｜ｐｏ−１はスパース正則化項に関係する項であり、Ｏｉ，ｊ，ｔの要素をスパースにするように誘導する効果をもつものであればこれ以外の形のもので代替してもよい。
【００４６】
次に、第２の実施形態における収束判定部９の処理動作を説明する。収束判定部９は、反復計算が所定の回数を満たしたか否か、あるいは、反復計算においてパラメータの更新の変化率が所定値以下になったか否か、あるいは、目的関数値の変化率が所定値以下になったか否かを判定する。例えば、目的関数は
【数７３】
0005580585.tif 000074
により計算する。ただし、Ｓ（Ｇ，Ｏ）はＧ，Ｏをスパースな解へと誘導する正則化項であり、ここでは
【数７４】
0005580585.tif 000075
と定義する。
【００４７】
以上説明したように、音響信号をノートごとの信号に分解するために、従来の非負値行列分解（ＮＭＦ）で用いられる分解要素（行列Ｕ）について、さらに音の立ち上がりや減衰のパターンの情報を表現できる分解（式（４）及び式（５））を導入した新たなモデルを用い、モデルの各パラメータを推定することにより、複数の音が混在している音響信号からの特定音の検出、複数の音が混在している音響信号からの特定音の抽出、および複数の音が混在している音響信号からの特定音の加工などに利用することができる。
【００４８】
なお、図１に示す時間周波数解析部２、初期設定部３、スペクトログラムモデル算出部４、スペクトル基底更新部５、パワーエンベロープ基底更新部６、パワーエンベロープ基底アクティビティ更新部７、パラメータ規格化部８、収束判定部９及びパラメータ出力部１０の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより信号分析処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、ホームページ提供環境（あるいは表示環境）を備えたＷＷＷシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
【００４９】
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。
【産業上の利用可能性】
【００５０】
複数の音が混在している音響信号からの特定音の検出、複数の音が混在している音響信号からの特定音の抽出、および複数の音が混在している音響信号からの特定音の加工などを行うことが不可欠な用途に適用できる。
【符号の説明】
【００５１】
１・・・信号データ記憶部、２・・・時間周波数解析部、３・・・初期設定部、４・・・スペクトログラムモデル算出部、５・・・スペクトル基底更新部、６・・・パワーエンベロープ基底更新部、７・・・パワーエンベロープ基底アクティビティ更新部、８・・・パラメータ規格化部、９・・・収束判定部、１０・・・パラメータ出力部、１１・・・パラメータ記憶部
【図１】
0005580585.tif 000076
【図２】
0005580585.tif 000077