(57)【要約】【目的】倍音構造や雑音の影響を軽減可能とし、かつ推定精度が良好な「メロディーの類似度算出方法及び類似度算出装置」を提供することである。【構成】２つの楽曲におけるメロディーの類似度を算出するメロディーの類似度算出に際して、比較する２つの楽曲信号のそれぞれに含まれるベース音の時間的変化を示す音高推移特性と各楽曲信号に含まれるベース音以外の楽器音のエネルギーの時間的変化を示すエネルギー特性とを取得する。ついで、各楽曲の音高推移特性の類似度を計算すると共に、各楽曲のベース音以外の楽器音のエネルギー特性の類似度を計算し、各類似度を用いてメロディーの類似度を算出する。【選択図】 図１０
2008257020.tif 000002
【特許請求の範囲】
【請求項１】
２つの楽曲におけるメロディーの類似度を算出するメロディーの類似度算出方法において、
比較する２つの楽曲信号のそれぞれに含まれるベース音の時間的変化を示す音高推移特性を取得する第１ステップ、
前記各楽曲信号に含まれるベース音以外の楽器音のエネルギーの時間的変化を示すエネルギー特性を取得する第２ステップ、
各楽曲の音高推移特性の類似度を計算する第３ステップ、
各楽曲のベース音以外の楽器音のエネルギー特性の類似度を計算する第４ステップ、
前記各類似度を用いてメロディーの類似度を算出する第５ステップ、
を有することを特徴とするメロディーの類似度算出方法。
【請求項２】
前記第２ステップにおいて、各楽曲信号の音名毎に、該音名の周波数成分のエネルギー及びその２k倍（ｋ＝１，２、・・・）の周波数成分のエネルギーを合計したエネルギーの時間変化特性を求め、
前記第４ステップにおいて、２つの楽曲信号の音名毎のエネルギー特性の相違度を計算し、該音名毎の相違度を合計し、該合計値の逆数により前記エネルギー特性の類似度を計算する、
ことを特徴とする請求項１記載のメロディーの類似度算出方法。
【請求項３】
前記第５ステップにおいて、前記音高推移特性の類似度とエネルギー特性の類似度を乗算してメロディーの類似度を算出する、
ことを特徴とする請求項１記載のメロディーの類似度算出方法。
【請求項４】
２つの楽曲におけるメロディーの類似度を算出するメロディーの類似度算出装置において、
比較する２つの楽曲信号のそれぞれに含まれるベース音の時間的変化を示す音高推移特性を取得する音高推移取得部、
前記各楽曲信号に含まれるベース音以外の楽器音のエネルギーの時間的変化を示すエネルギー特性を取得するエネルギー特性取得部、
各楽曲の音高推移特性の類似度を計算する第１の類似度計算部、
各楽曲のベース音以外の楽器音のエネルギー特性の類似度を計算する第２の類似度計算部、
前記各類似度を用いてメロディーの類似度を算出するメロディー類似度算出部、
を有することを特徴とするメロディーの類似度算出装置。
【請求項５】
前記エネルギー特性取得部は、各楽曲信号の音名毎に、該音名の周波数成分のエネルギー及びその２k倍（ｋ＝１，２、・・・）の周波数成分のエネルギーを合計したエネルギーの時間変化特性を取得する複数の音名エネルギー取得部を備え、
前記第２の類似度計算部は、２つの楽曲信号の音名毎のエネルギー特性の相違度を計算する相違度計算部、音名毎の相違度を合計し、該合計値の逆数により前記エネルギー特性の類似度を計算する計算部、
を有することを特徴とする請求項４記載のメロディーの類似度算出装置。
【請求項６】
前記メロディー類似度算出部は、前記音高推移特性の類似度とエネルギー特性の類似度を乗算してメロディーの類似度を算出する乗算部、
を備えることを特徴とする請求項４記載のメロディーの類似度算出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、メロディーの類似度算出方法及び類似度算出装置に係わり、特に、２つの楽曲におけるメロディーの類似度を算出するメロディーの類似度算出方法及び類似度算出装置に関する。
【背景技術】
【０００２】
楽曲のディジタル化に伴い、楽曲名、作曲者名等のキーワードによる検索だけでなく、音楽信号を直接用いた検索が可能となった。そこで、同一ジャンルの楽曲の検索や類似する楽曲の検索等、音楽信号の内容に基づく様々な検索技術の実現が期待されている。このような検索技術の実現のために、楽曲間の類似度を算出する手法が提案されている(非特許文献１，２参照)。
【非特許文献１】北川高嗣、中西崇文、清水康、”楽曲メディアデータを対象としたメタデータ自動抽出方式の実現とその意味的楽曲検索への適用、”電子情報通信学会論文誌、vol.J85-D-I, no. 6, pp. 512−526、2002
【非特許文献２】蔭山哲也、高島洋典、”ハミング歌唱を手がかりとするメロディー検索、“”電子情報通信学会論文誌、vol.J77-DII, no. 8, pp.1543-1551, 1994
【発明の開示】
【発明が解決しようとする課題】
【０００３】
非特許文献１の手法では、音楽信号より算出された複数の特徴量に基づき、楽曲にメタデータとして単語を付与する。さらに、人間の感性に基づく単語間の関係を用いて楽曲間の類似度を与え
る。従って、特徴量と単語を適切に設定することにより、楽曲の内容に基づく類似度算出を可能とする。しかしながら、メタデータが人間の感性に基づく単語から構成されるため、設定が困難であり、類似度の定量化も困難となる問題がある。また、類似度算出において楽曲の構造やメロディー等、音楽の構成要素を反映することは困難である。
一方、非特許文献２の手法では、歌唱から得られる音声と楽曲の類似度算出を可能とする。しかしながら、同一楽曲の検索を目的とするため、異なる楽曲間での類似度算出には十分ではない。
以上から本発明の目的は、上記従来の問題を解決する新たなメロディーの類似度算出手法及び算出装置を提供することである。
本発明の別の目的は、倍音構造や雑音の影響を軽減可能とし、推定精度が良好なメロディーの類似度算出手法及び算出装置を提供することである。
また、本発明の別の目的は、メロディーの伸縮や欠落が生じた場合にも類似度算出が可能なメロディーの類似度算出手法及び算出装置を提供することである。
【課題を解決するための手段】
【０００４】
本発明の第１の態様は、２つの楽曲におけるメロディーの類似度を算出するメロディーの類似度算出方法であり、比較する２つの楽曲信号のそれぞれに含まれるベース音の時間的変化を示す音高推移特性を取得する第１ステップ、前記各楽曲信号に含まれるベース音以外の楽器音のエネルギーの時間的変化を示すエネルギー特性を取得する第２ステップ、各楽曲の音高推移特性の類似度を計算する第３ステップ、各楽曲のベース音以外の楽器音のエネルギー特性の類似度を計算する第４ステップ、
前記各類似度を用いてメロディーの類似度を算出する第５ステップを有している。
前記第２ステップにおいて、各楽曲信号の音名毎に、該音名の周波数成分のエネルギー及びその２k倍（ｋ＝１，２、・・・）の周波数成分のエネルギーを合計したエネルギーの時間変化特性を求め、前記第４ステップにおいて、２つの楽曲信号の音名毎のエネルギー特性の相違度を計算し、該音名毎の相違度を合計し、該合計値の逆数により前記エネルギー特性の類似度を計算する。また、前記第５ステップにおいて、前記音高推移特性の類似度とエネルギー特性の類似度を乗算してメロディーの類似度を算出する。
【０００５】
本発明の第２の態様は、２つの楽曲におけるメロディーの類似度を算出するメロディーの類似度算出装置であり、比較する２つの楽曲信号のそれぞれに含まれるベース音の時間的変化を示す音高推移特性を取得する音高推移取得部、前記各楽曲信号に含まれるベース音以外の楽器音のエネルギーの時間的変化を示すエネルギー特性を取得するエネルギー特性取得部、各楽曲の音高推移特性の類似度を計算する第１の類似度計算部、各楽曲のベース音以外の楽器音のエネルギー特性の類似度を計算する第２の類似度計算部、前記各類似度を用いてメロディーの類似度を算出するメロディー類似度算出部を有している。
前記エネルギー特性取得部は、各楽曲信号の音名毎に、該音名の周波数成分のエネルギー及びその２k倍（ｋ＝１，２、・・・）の周波数成分のエネルギーを合計したエネルギーの時間変化特性を取得する複数の音名エネルギー取得部を備え、前記第２の類似度計算部は、２つの楽曲信号の音名毎のエネルギー特性の相違度を計算する相違度計算部、音名毎の相違度を合計し、該合計値の逆数により前記エネルギー特性の類似度を計算する計算部を有している。また、前記メロディー類似度算出部は、前記音高推移特性の類似度とエネルギー特性の類似度を乗算してメロディーの類似度を算出する乗算部を備えている。
【発明の効果】
【０００６】
本発明によれば、ベース音以外の楽器音のエネルギーとして、各音名が示す周波数のエネルギーを用いるため、倍音構造、雑音の影響を軽減可能であり、推定精度が良好なメロディーの類似度算出手法及び算出装置を提供することができる。
また、本発明によれば、低周波数域に基本周波数を持つベース音を併せて用いることで、倍音構造の影響をより軽減した類似度算出を可能とする。
さらに、本発明によれば、類似度の算出にはDTW を用いるため、メロディーの伸縮や欠落が生じた場合にも類似度算出が可能となる。
【発明を実施するための最良の形態】
【０００７】
（Ａ）本発明の概略
本発明は、音楽の構成要素であるメロディーに着眼する。音楽におけるメロディーとは、複数の音源により構成される基本周波数の時間推移である(石桁真礼生、“楽典−理論と実習、”音楽之友社発行、1980を参照)。本発明は、このメロディーの定義に従い、メロディーがベース音と、それ以外の楽器音から構成されると仮定する。さらに、この仮定に基づき、ベース音が示すエネルギーの推移、及びベース音以外の楽器音が示すエネルギーの推移についてマッチング処理を施すことで類似度を得る。ベース音が示すエネルギーとしてはベース音が存在する周波数域のパワースペクトルを用い、その他の楽器音が示すエネルギーとしてはＣ，Ｄ，Ｅ，・・・・等の音名が示す周波数のエネルギーを用いる。図１は音名説明図であり、１は日本語の音名(ハニホヘトイロハ)、２は英語標記の音名（CDEFGAB）、３は音符、４はハ長調での階名である。
ところで、音楽信号には以下の２つの特徴がある。第１の特徴は、楽器音に基本周波数の倍音が多く含まれる(以降、倍音構造という) 点であり、このため周波数域が高くなるに従い、基本周波数の特定が困難となる。第２の特徴は、楽曲中には発音の際に発生する擦弦音等の雑音が含まれる点であり、このため、音階上に存在しない周波数が楽器音の基本周波数として推定される可能性がある。
本発明はベース音が示すエネルギーの推移およびベース音以外の楽器音が示すエネルギーの推移、例えば、音名のエネルギーの推移を用いることにより、上記２つの特徴を有する音楽信号に有効に対応することができる。
本発明の実施例を説明する前に、本発明で使用するDTW(Dynamic Time Warping)技術と、音高推移特性の取得方法と、ベース音以外の楽器音が示すエネルギーの算出法を説明する。
【０００８】
（B）DTW技術
最初に、一般的なDTW について説明し、しかる後、本発明の類似度算出に用いるDTWの構成について説明する。
DTWとは、二つの一次元信号に伸縮を施し、信号間の相違度を算出する手法である。このため、時間軸方向に伸縮の生じる信号の比較に有効である。特に音楽信号では、演奏速度の変化が頻繁に発生することから、相違度より求められる類似度の算出にDTWを用いることは有効である。以降、相違度算出において、参照する信号を参照パターン、参照パターンとの相違度を求める信号を被参照パターンと呼ぶ。
まず、DTW によるパターン間の相違度の算出について説明する。長さI の一次元の参照パターンに含まれる各要素を順にa1，a2，…，aI とし、長さJ の被参照パターンに含まれる各要素を順にb1，b2，…，bJ と表現する。さらに、各パターンの位置集合を｛a1，a2，…，aI｝，[ b1，b2，…，bJ ] で表現すると、パターンの各要素間の対応を決定する伸縮写像ｗ：[1,2,…,I]→｛1,2,…,J｝は以下の性質を満たす。すなわち、
「wはパターンの始点、終点を一致させるもので、
【数１】
2008257020.tif 000003
である。また、wは単調写像であり、次式
【数２】
2008257020.tif 000004
が成立する。
【０００９】
このような写像w を用いたとき、パターン間の相違度の算出は図２における格子点(b1,a1) から格子点(bJ，aI) までの最短経路の探索問題に置換することができる。そこで、DTW では、
「初期状態の最初の決定が何であろうとも、以後の決定は最初の遷移から生じた状態に関して適切でなければならない」
という最適性の原理に基づいて上記の経路探索問題を解く。すなわち、全体の経路長を部分の経路長の和で求める。部分の経路長は、経路上の格子点(aj，bi) におけるコストd (j，i)、及び2 つの格子点(aj，bi)，(bj-1，ai-1)間の移動コストcj,i(j-1，i-1)に加え、cj,i(j，i-1)，cj,i(j-1，i)の３つを用いて算出する。部分の経路長の算出法を図２に示す。ここで、格子点上のコストd (j，i)は参照パターンと被参照パターンの間で対応する要素が異なる場合のペナルティであり、たとえば、要素が同じであれば０、異なればαである。また、移動コストcj,i(j-1，i-1) は参照パターンと被参照パターンの間で伸縮が生じた場合、格子点(bj-1，ai-1) から格子点(bj，ai) に移動するペナルティであり、たとえばi=jであれば０、パターンのずれによりi≠jであればβである。
上記のコストに基づいて部分の経路長を算出し、経路全体のコストが最小となる部分経路を選択する。最後に、選択された部分経路毎のコストの和を算出することで、全体の経路長が得られる。以上より、パターンの部分毎の相違度からパターン全体の相違度を得ることが可能となる。
【００１０】
本発明では、DTW を音楽信号に適用することから、音楽信号の類似度算出における特徴を考慮し、さらに詳細な相違度の算出法を提案する。本発明は、音楽の特徴として、同一楽曲の演奏速度が異なる場合にも、楽譜上の音符が欠落することがない点に着眼する。この特徴を、格子点間の移動による相違度算出に適用すると、参照パターンに含まれる全ての要素が被参照パターンに含まれ、要素間の対応が決定されることを意味する。従って、伸縮写像w は次式に示す傾斜制限を加えることが可能となる。
【数３】
2008257020.tif 000005
本発明上記の条件に従ってDTW による相違度の算出を行う。これより相違度は、(4)式を用いて経路長を漸化的に求めることで算出可能となる。
【数４】
2008257020.tif 000006
例えば、(４)式より、図４に示す格子点(1,1)の相違度D(１,１)は
D(１,１)=d(1,1)
であり、格子点（1,2)の相違度D(1,2）は
D(1,2）=d(1,2)+c1,2(1,1,)+D(1,1)
であり、格子点（2,1)の相違度D(2,1）は
D(2,1）=d(2,1)+c2,1(1,1,)+D(1,1)
である。同様にして第1行目の格子点の相違度及び第1列目の格子点の相違度が求まる。そして、以後、順次(4)式により他の格子点の相違度が求まってゆき、最終的に格子点（5,6）の相違度D（5,6）が求まる。この相違度D（5,6）が図3の2つのパターンa:｛a1，a2，…，a6｝，b:[ b1，b2，…，b5 ]の相違度となる。
【００１１】
（Ｃ）音高推移特性の取得
本発明では、楽曲中のベース音の推移として、ベース音が示す音高の推移を用いる。音高とは、楽譜上に記載される各音符が示す基本周波数である。従って、音高の推移はベース音に含まれる主要な周波数におけるエネルギーの推移を意味する。ベース音の音高は、下記文献１に記載されている音高推定手法に従って推定する。
文献１：今野聡司、他“音楽信号の低周波数域における楽曲の構造に着眼した時間周波数解析に関する検討"、映像情報メディア学会技術報告、vol.29, no.46, pp.13-16, 2005
図５は音高推定装置の構成図である。バンドパスフィルタ１はベース音が存在する周波数帯域に基づいて処理対象の楽曲信号に含まれる40〜250Hzの信号成分を通過する。パワースペクトル算出部２はバンドパスフィルタ出力信号からパワースペクトルを算出する。評価関数計算部３は該パワースペクトルの時間軸方向、及び周波数軸方向にガウス関数に基づく重みを付加し、最後に、音高推定部４は重み付けされた各時刻のパワースペクトルにおいて最大のエネルギーを与える周波数を音高として推定する。
【００１２】
すなわち、時刻t(0 ≦t ≦ T)、周波数f においてパワースペクトル算出部２が算出するエネルギーをP(ｔ，f)とすれば、評価関数計算部３は(5)式により重み付けされたパワースペクトルR（t，f）を評価関数として出力する。
【数５】
2008257020.tif 000007
なお、vt(s)は時間軸の重みであり、
【数６】
2008257020.tif 000008
ただし、σは音の持続時間の指標となる定数である。また、（5）式においてw(f)は周波数軸の重みであり、
【数７】
2008257020.tif 000009
ただし、mを自然数として、
【数８】
2008257020.tif 000010
で示されるFm は、DI（Musical Instrument Digital Interface）のm 番目のノートにおける周波数を表す。
(５)式の評価関数R（t，f）は、(6) 式の時間軸方向の重みにより、一定時間持続する基本周波数を音高として推定可能とする。また、(7) 式に示す周波数軸方向の重みにより、音階上に存在する周波数のみを音高として推定可能とする。以降、評価関数R（t，f）の各時刻t において最大値を与える周波数f をベースの音高とし、B(t) と表す。
【００１３】
（Ｄ）その他の楽器音によるエネルギーの推移
次に、ベース音以外の楽器音が示すエネルギーの算出法について説明する。一般的な音楽の構成では、主にベース音が楽曲の最低音となるため、その他の楽器音はベース音の周波数域より高い周波数を示す。また、ベース音より高い周波数域で、各音名は図６に示す周波数を持ち、各周波数の2k(k = 1，2，．．．) 倍の周波数も同一の音名として扱われる。そこで、本発明では、ベース音以外の楽器音が示すエネルギーを、ベース音の周波数より高く、かつ音名が持つ周波数のエネルギーとする。さらに、各音名が示す周波数のエネルギーには、図６の2k 倍の周波数が示すエネルギーの和を用いる。これにより、本発明は複数の楽器による倍音構造を軽減し、音高の推定が困難な周波数域に存在する楽器音についても類似度算出に用いることを可能とする。以上の手順概要を図７ に示し、以下で詳細を説明する。なお、音楽信号の信号長はT 秒、サンプリングレートはfs とし、また、時刻t(0 ≦t≦ T)および周波数f におけるエネルギーはパワースペクトルより算出されてP（t，f）と表す。
【００１４】
1 ）音名が示す周波数のエネルギーの算出：
パワースペクトルから、各音名が示す周波数のエネルギーを算出する。図６ において音名X に対応する周波数をfX として、音名X が示す周波数のエネルギーPX(t) を次式で定義する。
【数９】
2008257020.tif 000011
ただし、K はを越えない任意の整数とする。(9)式により各音名が示す周波数のエネルギーを定義することで、倍音の影響が軽減可能となる。
2 ）エネルギーの割合の算出：
1）で得られた各音名が示す周波数のエネルギーを全周波数域に対するエネルギーの割合で表現する。これにより、音名毎に時間軸方向での比較が可能となり、推移を得ることが可能となる。音名X が示す周波数のエネルギーの割合px(t) は次式で示される。
【数１０】
2008257020.tif 000012
以上を全てのt、X について施し、得られたpx(t) をベース音以外の楽器音におけるエネルギーの推移として用いる。
【００１５】
図８は各音名のエネルギー特性を取得する音名エネルギー取得部の構成図であり、音名の周波数成分、該周波数の2k(k = 1，2，．．．) 倍の周波数成分をそれぞれ通過するバンドパスフィルタ１１1〜１１n、各バンドパスフィルタ出力を合成することにより音名エネルギーPX(t)を出力する合成部１２、（10）式によりエネルギーの割合px(t)を出力するエネルギー割合計算部13を備えている。
【００１６】
(E)メロディー類似度算出法
図９は２つの楽曲におけるメロディーの類似度を算出する本発明のメロディーの類似度算出方法の処理フローである。
本発明は、まずメロディーがベース音、及びベース音以外の楽器音から構成されると仮定する。これは、ベース音と他楽器音により同時に発音される音が、メロディーの特徴を決定する和音や調の指標となるためである。本発明は上記の仮定に基づき、それぞれの楽器音のエネルギーにDTW を適用することで類似度を算出する。
すなわち、最初に比較する２つの楽曲信号のそれぞれに含まれるベース音の時間的変化を示す音高B(t)を取得し(ステップ１０１)、また、各楽曲信号に含まれるベース音以外の楽器音のエネルギーの時間的変化を示すエネルギー特性px(t)を取得する(ステップ１０２)。次に(4)式に基づくDTW を用いて、各楽曲の音高の相違度を求め、その逆数を演算することで類似度Ｓａを計算する(ステップ１０３)。しかる後、同様に、DTW を用いて、各楽曲のベース音以外の楽器音のエネルギー特性の相違度を求め、その逆数を演算することで類似度Sbを計算する(ステップ１０４)。最後に、音高B(t)、ベース音以外の楽器音のエネルギーpx(t)のそれぞれについて得られた類似度の積（Sa×Sb）を算出することでメロディーの類似度とする(ステップ１０５)。
【００１７】
以下で類似度算出の詳細について説明する。なお、DTW による相違度算出では、各楽曲信号を十分に小さい時間で分割し、一方を参照パターン、他方を被参照パターンとして用いる。
1）音高B(t) の相違度算出：
図5で説明したように、参照側のメロディーと被参照側のメロディーそれぞれについてベースの音高を推定し、DTWによる相違度算出を行う。ここで、DTW の(4)式で用いる各コストは以下のように設定する。
【数１１】
2008257020.tif 000013
【数１２】
2008257020.tif 000014
ただし、α＞βとする。これにより、メロディーの不一致によるコストと比較して、演奏速度の変化等に伴うメロディーのずれに対するコストが小さくなる。以上により得られた相違度をDb と表す。
【００１８】
2）px(t) の相違度算出：
参照側のメロディーと被参照側のメロディーそれぞれについて、ベース音以外の楽器音の音名X毎のネルギーpxr (t), pxi (t)を算出し、これらを用いて音名毎にDTW による相違度Daxの算出を行う。従って、相違度は音名の数である12個得られる。そして、ベース以外の楽器音の相違度は音名毎に得られた相違度の和により定義する。すなわち、音名X について得られる相違度をDax とすると、ベース以外の楽器による音の相違度Da は次式で表される。
Da = DaC+DaCis+DaD+DaDis+DaE+DaF+DaFis +DaG +DaGis +DaA +DaB +DaH   （13）
なお、DTWによる相違度算出に用いるコストは以下のように設定する。
【数１３】
2008257020.tif 000015
【数１４】
2008257020.tif 000016
以上、（13）式により、全ての音名が示す周波数のエネルギーの推移を用いた相違度算出が可能となる。また、（14）式に示すコストを設定することで、エネルギーの大きな周波数に対応する音名が、相違度全体に与える影響を増加する。これにより、メロディーを構成する主要な周波数成分を反映した相違度算出が可能となる。
【００１９】
3）各類似度の算出：
上記1）、2）で得られた相違度Db，Daの逆数を算出することで類似度とし、Sb， Sa と表す。
4）メロディーの類似度算出：
上記3）で得られた類似度Sb， Saから次式
S = Sb×Sa     （16）
によりメロディーの類似度S を算出する。
以上の処理により、メロディー間の類似度を算出する。メロディー間の類似度は、ベース音とその他の楽器による音のエネルギーそれぞれの類似度から算出する。このため、倍音構造や雑音の影響を軽減可能とし、メロディーを構成する複数の楽器音を反映した類似度算出が可能となる。また、DTW を用いるため、メロディー間に伸縮や欠落が生じた場合にも類似度算出が可能となる。
【００２０】
(F)メロディー類似度算出装置
図１０は２つの楽曲におけるメロディーの類似度を算出する本発明のメロディーの類似度算出装置の構成図である。参照側及び被参照側の音高推移取得部５１，７１はそれぞれ比較する２つの楽曲信号に含まれるベース音の時間的変化を示す音高推移特性を取得し、参照側及び被参照側のエネルギー特性取得部５２，７２は楽曲信号に含まれるベース音以外の楽器音のエネルギーの時間的変化を示すエネルギー特性を取得する。エネルギー特性取得部５２は、12個の音名のエネルギー特性をそれぞれ取得する１２個のエネルギー取得部５２1〜５２12を備えている。図示しないがエネルギー特性取得部７２も同様の構成を備えている。相違度計算部５３はDTW を用いて、各楽曲の音高の相違度Dbを求め、演算部５４は相違度の逆数を演算することで類似度Sbを計算する。
各音名の相違度計算部５５1〜５５12はDTW を用いて音名毎にエネルギー特性の相違度を計算し、演算部５６は音名毎のエネルギー特性の相違度を合計してベース音以外の楽器音のエネルギー特性の相違度Daを求め（（13）式参照）、演算部５７は相違度Daの逆数を演算することで類似度Saを計算する。メロディー類似度計算部５８は、音高B(t)、ベース音以外の楽器音のエネルギーpx(t)のそれぞれについて得られた類似度Sb，Saの積Sを算出し、該積をメロディーの類似度Sとして出力する。
【００２１】
（Ｇ）実験
本発明の有効性を確認するために実験を行った。実験にはCD から直接入手した44.1kHz、10秒の音楽信号を3 つ用いた。音楽信号は同一の楽曲から２ 箇所、異なる楽曲から１ 箇所切り出した信号であり、以降それぞれを順に信号1、 信号2、 信号3 とする。なお、本発明で用いるパラメータはそれぞれσ = 0:1、α＝3、β＝1、γ＝0.1、ｋ＝10とした。
実験結果を図１１、図１２に示す。図１１は各信号間で最終的に得られた類似度を示している。ただし、図における類似度は、左端に示す信号を参照パターン、右上端に示す信号を被参照パターンとして用いた場合の結果を表す。一方、図１２は、各時刻までで算出される類似度を示している。ただし、図１２（Ａ），（Ｂ），（Ｃ）は参照パターンを信号１、２、３とし、他の２つの信号を被参照パターンとして得られた類似度を示す。ここで、被参照パターンに用いた信号１、２、３は図中でそれぞれ実線、破線、点線で示してある。
【００２２】
まず、信号1 について考察を行う。図１１ より、信号3と比較して、信号2 との類似度が高いことが確認できる。また、図１２（Ａ） より、各時刻における類似度についても、同様に信号2 との類似度が高いことが確認できる。これは、信号１、２ が同一の楽曲の異なる時間を切り出した信号であるため、メロディー全体に渡り、同時に発音される音が類似するためと考えられる。また、図１２（Ａ）の３．５ 秒以降の類似度が特に高くなることが確認できる。これは、信号１ の演奏において、メロディーが強調される演奏がなされたことに起因すると考えられる。このため、 各音名が示す周波数のエネルギーが強くなり、楽曲２ との類似度が高くなったと考えられる。
次に、信号２ について考察を行う。信号１ と同様に、図１１より、信号３と比較して信号１との類似度が高いことが確認できる。しかしながら、図１２（Ｂ）より、信号１、３との類似度は各時刻で比較すると差が小さいことが確認できる。この点は、信号２が楽曲のサビであるため、同時に発音される楽器数が多く、信号のパワーが極めて強いことが原因と考えられる。このため、信号２と比較してパワーの小さい信号１、３ が類似すると推定されたと考えられる。
最後に、信号３について考察を行う。図１１より、信号３と異なる楽曲である信号１、２との類似度は低いことが確認できる。また、図１２より、各時刻においても信号1、２との類似度は低く、時間による変化が少ないことが確認できる。以上より、提案手法は楽曲の類似度を算出可能であることが確認できる。図１２では、信号１、２、３ いづれの類似度算出においても、時間軸に対して類似度が増加することが確認できる。また、時間の経過に従って類似度の変化が小さくなることが確認できる。これらは、類似度がDTW を用いた経路探索により得られるため、短い経路では類似度算出の精度が劣化することを意味する。このため、類似度の算出では、十分な信号長を持つ信号を用いる必要があると考えられる。
【図面の簡単な説明】
【００２３】
【図１】音名説明図である。
【図２】DTW(Dynamic Time Warping)技術の第１の説明図である。
【図３】DTW技術の第２の説明図である。
【図４】本発明のDTWの説明図である。
【図５】音高推定装置の構成図である。
【図６】各音名の周波数説明図である。
【図７】本発明の手順概要説明図である。
【図８】音名のエネルギー特性を取得する音名エネルギー取得部の構成図である。
【図９】２つの楽曲におけるメロディーの類似度を算出する本発明のメロディーの類似度算出方法の処理フローである。
【図１０】２つの楽曲におけるメロディーの類似度を算出する本発明のメロディー類似度算出装置の構成図である。
【図１１】第１の実験結果である。
【図１２】第２の実験結果である。
【符号の説明】
【００２４】
５１，７１  参照側及び被参照側の音高推移取得部
５２，７２  参照側及び被参照側のエネルギー特性取得部
５２1〜５２12  音名エネルギー取得部
５３  音高の相違度計算部
５４  演算部
５５1〜５５12  各音名の相違度計算部
５６、５７  演算部
５８  メロディー類似度計算部
【図１】
2008257020.tif 000017
【図２】
2008257020.tif 000018
【図３】
2008257020.tif 000019
【図４】
2008257020.tif 000020
【図５】
2008257020.tif 000021
【図６】
2008257020.tif 000022
【図７】
2008257020.tif 000023
【図８】
2008257020.tif 000024
【図９】
2008257020.tif 000025
【図１０】
2008257020.tif 000026
【図１１】
2008257020.tif 000027
【図１２】
2008257020.tif 000028