多和音の音名と音高推定手法
【課題】
４和音以上の多和音について各音高を精度良く推定するために、縦続接続くし形フィルタシステムを用い、フィルタ出力の時間波形における振幅の変化や周波数の変化から多和音の音高を推定する手法を提供することである。
【解決手段】
実楽器音の多和音を縦続接続くし形フィルタで処理する時、入出力振幅比だけでなく、入出力周波数比にも注目する。その周波数比の大きな変化点を示すくし形フィルタから、対応する音名が推定され、さらに、そのくし形フィルタの次数を各オクターブに対応するように変化させ、その出力波形の振幅変化か周波数変化からその音のオクターブが推定され、多和音のうちの１音の音高が推定できる。次に、推定した１音に対応するくし形フィルタをフィルタシステムの先頭に移動させる。この操作を繰り返して、すべての音高を推定する。
４和音以上の多和音について各音高を精度良く推定するために、縦続接続くし形フィルタシステムを用い、フィルタ出力の時間波形における振幅の変化や周波数の変化から多和音の音高を推定する手法を提供することである。
【解決手段】
実楽器音の多和音を縦続接続くし形フィルタで処理する時、入出力振幅比だけでなく、入出力周波数比にも注目する。その周波数比の大きな変化点を示すくし形フィルタから、対応する音名が推定され、さらに、そのくし形フィルタの次数を各オクターブに対応するように変化させ、その出力波形の振幅変化か周波数変化からその音のオクターブが推定され、多和音のうちの１音の音高が推定できる。次に、推定した１音に対応するくし形フィルタをフィルタシステムの先頭に移動させる。この操作を繰り返して、すべての音高を推定する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、楽音を解析し楽譜として出力する採譜システムの実現において、出力波形における振幅の変化や周波数の変化に注目して多和音を対象とした音高推定手法に関するものである。
【背景技術】
【０００２】
音楽信号を楽譜に書き起こす処理のことを採譜と呼ぶ。従来、採譜は、音楽に関する知識や経験を持つ専門家にしか行えないものであった。この作業を計算機によって自動的に行う技術が自動採譜である。ここで、自動採譜システムを計算機上に実現するためには、楽音の音の高さ（音高）を推定する処理が必要となる。音楽信号は、基本周波数成分とその整数倍の周波数を有する多数の高調波成分から構成されており、音高推定とは主に基本周波数（ピッチ）を推定する処理のことを指す。
【０００３】
従来技術として、マイクなどの入力装置、あるいはＣＤ再生などによって入力された音響信号を楽譜に直す採譜システム、例えば、ヤマハ社のＸＧワークスＶ２.０、ＵＳＡワイルドキャットキャニオン社のオートスコアなどがある。これらのほとんどは独唱や単独の楽器によって演奏される単音を対象にしたものであり、周波数特性の異なる複数の音源によって構成された和音までを扱うことが可能なものは、いまだ実用化されていない。単音に関する特許文献として、特開平０５−１０８０６２号公報、特開平０５−１２７６６８号公報、特開平０５−１６５４６４号公報、特開平０９−０７２７７９号公報がある。
【０００４】
一方、これまでの和音を対象とした音高推定のための技術として、和音の採譜に関する特許文献では、特開平０４−２７８５４４号公報、特開平０５−１００６６０２号公報、特開平０５−１７３５５７号公報、特開２００５−２０２３５４Ａ号公報がある。しかし、これらに開示されている技術では、１０和音などの多和音に対応することは難しい。
【０００５】
また、非特許文献１では、楽音のすべての成分を消去することができるくし形フィルタを使用し、１オクターブにある１２音の各音に対応する１２個のくし形フィルタを用いることで、何和音でもその各音の音高を推定する方法が提案されている。しかし、使用されている楽音はＭＩＤＩ音源（電子楽器）に対する多和音の音高推定手法に関するものである。この手法を実楽器音源データ（ＲＷＣ音楽データベースに収録、非特許文献３を参照）に使用した場合、実楽器のもつ非定常性のため、フィルタの振幅に注目した音高推定は困難になる。
【０００６】
最後に、非特許文献２では、非特許文献1で明らかになった実楽器の非定常性に対応するための音高推定法が提案されている。並列構成くし形フィルタの入出力振幅比を基本にして音高を推定する方法である。しかしながら、和音の数だけのくし形フィルタが接続されることになり、和音数の増加とともに音高を推定する能力は低下する。
【非特許文献１】三輪、田所、斉藤、「くし形フィルタを利用した採譜のための異楽器音中のピッチ推定」、信学論、Ｖｏｌ．Ｊ８１−Ｄ−ＩＩ、Ｎｏ．９、ｐｐ．１９６５−１９７４、１９９８
【非特許文献２】森田、山口、田所、「並列構成くし形フィルタの出力値に注目した採譜のための音高推定法」、信学論、Ｖｏｌ．Ｊ８７Ｄ−ＩＩ、Ｎｏ．１２、ｐｐ．２２７１−２２７９、２００４
【非特許文献３】後藤、橋口、西村、「ＲＷＣ研究用音楽データベース：音楽ジャンルデータベースと楽器音データベース」、情処研報、２００２−ＭＵＳ−４５、ｐｐ．１９−２６、２００２
【発明の開示】
【発明が解決しようとする課題】
【０００７】
楽音は、基本周波数とその倍音成分から構成されるため、和音の場合、各音の成分が非常に複雑な関係になり、信号処理の一般的手法である離散フーリエ変換を基本にしたアルゴリズムでは、一般に高精度の推定は３和音程度までの音高推定である。しかし、実際の演奏では４和音以上も存在することが多い。例えば、ピアノ演奏では、１０本の指を使えば１０和音も生成できる。そのため４和音以上の多和音についても、その各音高を精度良く推定する手法を開発することは、実用的な採譜システムを実現するために是非必要とされる技術である。
【０００８】
そこで、本発明で解決する課題は、縦続接続くし形フィルタシステムを用い、フィルタ出力における振幅の変化や周波数の変化から多和音の音高を推定する手法を提供することである。
【課題を解決するための手段】
【０００９】
実楽器音の多和音を縦続接続くし形フィルタで処理する時、くし形フィルタの入出力振幅比だけでなく、入出力周波数比にも注目する点が本発明のポイントである。楽音の非定常性が縦続接続くし形フィルタの周波数特性のために増幅され、振幅の大小だけでは音高推定が不可能である。しかし、楽音の基本となる各周波数成分がくし形フィルタでほぼ除かれた後に、非定常成分の低周波成分が出力される。このため、くし形フィルタの入出力周波数比に注目し、その周波数比の大きな変化点を示すくし形フィルタから、対応する音名（Ｃ、Ｃ＃、・・・、Ｂなど）が推定され、さらに、そのくし形フィルタの次数を各オクターブに対応するように変化させ、その出力波形の変化（振幅変化か周波数変化）からその音のオクターブが推定され、多和音のうちの１音の音高が推定できる。
【００１０】
このようにして推定した１音に対応したくし形フィルタを縦続接続くし形フィルタの先頭に移動させ、再度、波形の変化、すなわち振幅比や周波数比に注目して残りの音名（オクターブ推定の処理を含めると音高）の一つを推定し、再度このくし形フィルタを縦続接続くし形フィルタの先頭に移動させる。この操作を繰り返して、最終的に最初に移動させたくし形フィルタでの波形の変化に大きな変化が生じた時点で音高推定操作が終了する。そして、移動させたくし形フィルタから多和音の各音高が推定される。
【００１１】
対象音域の最低オクターブ中に存在する１２音に対応するくし形フィルタを縦続接続し、その各くし形フィルタの入出力振幅比、周波数比の変化を測定し、その変化の大きいくし形フィルタを順次、縦続接続くし形フィルタの先頭に移動させる。なお、縦続接続するくし形フィルタの数は、対象音域の最低オクターブよりひとつ低いオクターブに対応するくし形フィルタを使用すれば、等価的に７個で済まされることも可能である。
【発明の効果】
【００１２】
くし形フィルタは、周波数軸上に等間隔の零点を有しており、単音による音楽信号に含まれる全調波成分を除去する性質を有する。この性質を利用し、各音高に対応するくし形フィルタの出力から零出力を検出することにより、音高の推定が可能となる。また、くし形フィルタは遅延器のみを用いて構成されるため、時間領域においてはサンプル点間の差分のみで処理が可能であり、演算量の削減および処理系の簡単化を図ることができる。
【発明を実施するための最良の形態】
【００１３】
１オクターブ１２音に対応するくし形フィルタを図１のように縦続接続したシステムについて考える。このシステムは、原理的には、入力ｘ（ｎ）に含まれる音を１音ずつ除去し、入力中の全ての調波成分を除去した段階で、零値を出力する。ここで、このフィルタの最終出力の振幅特性は図２のようになる。原理的には、入力和音中の全ての音が消去された段階でくし形フィルタの出力は零となるが、実際の楽音では非定常性、周波数変動のため、完全な零出力とならず、むしろ図２の特性により、低域の周波数成分が増大し、振幅の大小での音高推定は困難になる。
【００１４】
そこで、振幅値に注目するのではなく、時間波形（周波数）の変化に注目した方法を提案する。縦続接続くし形フィルタに信号を入力し、和音の周波数成分が全て除去されると、くし形フィルタの出力は零出力にはならないが、その時点で時間波形は楽音の非定常性により、周波数変動分の低い周波数の出力波形になる。この時間波形の変化に注目し、和音の音高推定を行う。
【００１５】
例として、図３に５和音（Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）を入力した時の縦続接続くし形フィルタの各フィルタの出力を示す。各フィルタの出力波形を見ると、ＢからＤ＃までのフィルタの出力波形に対し、音高Ｄのフィルタまで通した時点で時間波形（周波数）は、低い周波数の出力波形になっているのがわかる。この変化により、入力和音にＤが含まれていると推定でき、ＢからＤまでの間の音高の組み合わせで入力和音が構成されているとも考えられる。
【００１６】
以上のことより、提案する音高推定アルゴリズムは図４のようになる。
【００１７】
（くし形フィルタによる音高推定の基本原理）
音高ｐの基本周波数をｆｐとした場合、音楽信号は基本周波数成分と、その整数倍の高調波成分から構成される。ここで、ｐ＝１（Ｃ）、２（Ｃ＃）、・・・、１２（Ｂ）である。この音楽信号を消去するくし形フィルタの伝達関数Ｈｐ（ｚ）は、数１で表され、その出力ｙｐ（ｎ）は、数２に示す減算で実現される。
【００１８】
【数１】
【００１９】
【数２】
【００２０】
ここで、Ｎｐは遅延数であり、サンプリング周波数をｆｓとするとＮｐは、数３で求められる。
【００２１】
【数３】
【００２２】
図５にオクターブ４の音名Ｅの音楽信号（基本周波数は３２９．６３Ｈｚ）の周波数特性と、対応するくし形フィルタの振幅特性を示す。くし形フィルタは、周波数領域において等間隔に零点が存在するような特性を有しており、音楽信号を構成する複数の周波数成分を同時に除去することができる。そのため、くし形フィルタの出力値から零出力を検出することで、音高の推定が可能となる。
【実施例】
【００２３】
対象音域をオクターブ３以上としたときは、オクターブ３の各音名に対応するくし形フィルタを図１のように１２個縦続接続し、このシステムに多和音を入力させる。この縦続接続くし形フィルタの振幅特性は、図２のようになり、オクターブ３以下の音域の振幅ゲインが大きくなっている。そのため、多和音の各周波数成分が除かれた後には、理想的にはくし形フィルタの出力振幅は零になるはずであるが、実楽器の多和音の非定常性成分（低周波成分）が出力されることになる。その出力振幅は後段のくし形フィルタになるほど大きくなり、入出力振幅比での音名推定は困難になる。しかし、くし形フィルタの入出力周波数比に注目することで、多和音がすべて除かれたかどうかを判断できる。
【００２４】
入出力周波数比に注目して多和音の各音名を推定する例を図６に示す。図６では、左側が縦続接続くし形フィルタの移動を示し、対応する右側は和音の全周波数成分が消去される前後の時間波形をｉｎｐｕｔとｏｕｔｐｕｔとしてそれぞれ示している。多和音の例として、５和音（オクターブ４：Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）を考える。縦続接続くし形フィルタ（オクターブ４の１２音に対応するくし形フィルタから構成）は、音名の高い順に並べている。より確実な推定をするためには、並べ方を逆にした二つの縦続接続くし形フィルタを使用するとよい。
【００２５】
最初に入出力周波数比が大きく変わるのは、音名Ｄを除くくし形フィルタである。なお、入出力振幅比(出力/入力)も１以下になっている。このＤに対応するくし形フィルタを縦続接続くし形フィルタの先頭に移動させ、２回目の入出力周波数比が大きく変わるくし形フィルタを見つける。今度は音名Ｆを除くくし形フィルタで検出される。このＦのくし形フィルタを先頭に移動させる。３回目の入出力周波数比の大きな変化は、音名Ｇに対応するくし形フィルタである。このＧのくし形フィルタを先頭に移動させる。４回目の入出力周波数比の大きな変化点は、音名Ａのくし形フィルタである。Ａのくし形フィルタを先頭に移す。そして、５回目の入出力周波数比の大きな変化するくし形フィルタは音名Ｂのくし形フィルタである。このＢに対応するくし形フィルタを先頭に移動し、6回目の入出力周波数比の大きな変化するくし形フィルタを検出する。今度は音名Ｄに相当するくし形フィルタであり、これは最初に先頭に移動されたくし形フィルタであるので、ここで音名推定の操作は終了する。
【００２６】
このようにして、各操作で先頭に移動されたくし形フィルタに相当する音名が、入力多和音の各音名に対応する。すなわち、入力されたＤ、Ｆ、Ｇ、Ａ、Ｂが推定されている。
【００２７】
次に、これらの各音名のオクターブンの推定は、各音名が推定された時点で、そのくし形フィルタの次数を２分の１、４分の１、・・・と変化させて、入出力周波数比の大きな変化が生じなく、または入出力振幅比(出力/入力)が１以上になる、次数の手前の次数が対応するオクターブとして決定される。そして、音名とオクターブから音高が決定される。なお、高調波成分のため、入出力周波数比（ゼロクロス数、ピーク数）の判定が入出力波形そのものから難しい場合は、平均化フィルタや自己相関関数の処理で明確にできる。この例の場合は、オクターブ推定でくし形フィルタの次数を２分の１にしたときに入出力周波数比の大きな変化が生じなくなる。よって、これらのオクターブは４と決定される。
【産業上の利用可能性】
【００２８】
本発明により、和音の数に制限されずに多和音の各音高を推定することができ、本格的な採譜システムを実現できる。これにより、音楽配信ビジネス、カラオケ産業での楽曲検索、民族音楽の記録の効率化、即興音楽の記録化、作曲作業の能率化および音楽教育への利用など、その応用範囲は多岐にわたる。音楽配信ビジネス（通信カラオケ、携帯電話やＰＨＳの着信メロディ）では、新曲をいち早く楽譜化することが必要であり、これを人手ではなく計算機による自動採譜システムで行なうことが熱望されている。この自動採譜システムとして４和音以上の多和音の精度良い音高推定が可能な本格的な採譜システムが実現されれば、この業界に多大な貢献をすることが可能である。
【００２９】
また、楽音のリアルタイム認識によるセッションシステム（相手の演奏の理解に基づいた演奏応答システム）のようなインタラクティブ音楽コンテンツの実現、さらには和音の組み合わせの多さを利用した個人認証システム（例えば、特許文献１を参照）への応用など考えられる。携帯電話の音源を利用した個人認証システムでは、採譜システムによる音高推定能力により、その性能（秘話性）が決定され、本発明による本格的な採譜システムは、そのための有望な技術になり、これに関係する産業界、ＩＴ社会での本発明の使用が期待される。
【特許文献１】特開２００１−３１８８９４号公報
【図面の簡単な説明】
【００３０】
【図１】くし形フィルタ（オクターブ３）の１２段縦続接続によるシステム構成
【図２】Ｈｐ（ｚ）（ｐは１から１２）を１２段接続した縦続接続くし形フィルタの振幅特性
【図３】５和音（Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）を入力したときの縦続接続くし形フィルタの各出力
【図４】図１のように構成されるくし形フィルタシステムの処理アルゴリズム
【図５】オクターブ４の音名Ｅの周波数特性とくし形フィルタの振幅特性
【図６】本発明の実施例に係る多和音（オクターブ４：Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）の音名推定過程
【符号の説明】
【００３１】
１ Ｈｐ（ｚ）＝１−ｚ−Ｎｐの特性を持つくし形フィルタ
２ くし形フィルタを１２段縦続接続構成にしたくし形フィルタシステム
【技術分野】
【０００１】
本発明は、楽音を解析し楽譜として出力する採譜システムの実現において、出力波形における振幅の変化や周波数の変化に注目して多和音を対象とした音高推定手法に関するものである。
【背景技術】
【０００２】
音楽信号を楽譜に書き起こす処理のことを採譜と呼ぶ。従来、採譜は、音楽に関する知識や経験を持つ専門家にしか行えないものであった。この作業を計算機によって自動的に行う技術が自動採譜である。ここで、自動採譜システムを計算機上に実現するためには、楽音の音の高さ（音高）を推定する処理が必要となる。音楽信号は、基本周波数成分とその整数倍の周波数を有する多数の高調波成分から構成されており、音高推定とは主に基本周波数（ピッチ）を推定する処理のことを指す。
【０００３】
従来技術として、マイクなどの入力装置、あるいはＣＤ再生などによって入力された音響信号を楽譜に直す採譜システム、例えば、ヤマハ社のＸＧワークスＶ２.０、ＵＳＡワイルドキャットキャニオン社のオートスコアなどがある。これらのほとんどは独唱や単独の楽器によって演奏される単音を対象にしたものであり、周波数特性の異なる複数の音源によって構成された和音までを扱うことが可能なものは、いまだ実用化されていない。単音に関する特許文献として、特開平０５−１０８０６２号公報、特開平０５−１２７６６８号公報、特開平０５−１６５４６４号公報、特開平０９−０７２７７９号公報がある。
【０００４】
一方、これまでの和音を対象とした音高推定のための技術として、和音の採譜に関する特許文献では、特開平０４−２７８５４４号公報、特開平０５−１００６６０２号公報、特開平０５−１７３５５７号公報、特開２００５−２０２３５４Ａ号公報がある。しかし、これらに開示されている技術では、１０和音などの多和音に対応することは難しい。
【０００５】
また、非特許文献１では、楽音のすべての成分を消去することができるくし形フィルタを使用し、１オクターブにある１２音の各音に対応する１２個のくし形フィルタを用いることで、何和音でもその各音の音高を推定する方法が提案されている。しかし、使用されている楽音はＭＩＤＩ音源（電子楽器）に対する多和音の音高推定手法に関するものである。この手法を実楽器音源データ（ＲＷＣ音楽データベースに収録、非特許文献３を参照）に使用した場合、実楽器のもつ非定常性のため、フィルタの振幅に注目した音高推定は困難になる。
【０００６】
最後に、非特許文献２では、非特許文献1で明らかになった実楽器の非定常性に対応するための音高推定法が提案されている。並列構成くし形フィルタの入出力振幅比を基本にして音高を推定する方法である。しかしながら、和音の数だけのくし形フィルタが接続されることになり、和音数の増加とともに音高を推定する能力は低下する。
【非特許文献１】三輪、田所、斉藤、「くし形フィルタを利用した採譜のための異楽器音中のピッチ推定」、信学論、Ｖｏｌ．Ｊ８１−Ｄ−ＩＩ、Ｎｏ．９、ｐｐ．１９６５−１９７４、１９９８
【非特許文献２】森田、山口、田所、「並列構成くし形フィルタの出力値に注目した採譜のための音高推定法」、信学論、Ｖｏｌ．Ｊ８７Ｄ−ＩＩ、Ｎｏ．１２、ｐｐ．２２７１−２２７９、２００４
【非特許文献３】後藤、橋口、西村、「ＲＷＣ研究用音楽データベース：音楽ジャンルデータベースと楽器音データベース」、情処研報、２００２−ＭＵＳ−４５、ｐｐ．１９−２６、２００２
【発明の開示】
【発明が解決しようとする課題】
【０００７】
楽音は、基本周波数とその倍音成分から構成されるため、和音の場合、各音の成分が非常に複雑な関係になり、信号処理の一般的手法である離散フーリエ変換を基本にしたアルゴリズムでは、一般に高精度の推定は３和音程度までの音高推定である。しかし、実際の演奏では４和音以上も存在することが多い。例えば、ピアノ演奏では、１０本の指を使えば１０和音も生成できる。そのため４和音以上の多和音についても、その各音高を精度良く推定する手法を開発することは、実用的な採譜システムを実現するために是非必要とされる技術である。
【０００８】
そこで、本発明で解決する課題は、縦続接続くし形フィルタシステムを用い、フィルタ出力における振幅の変化や周波数の変化から多和音の音高を推定する手法を提供することである。
【課題を解決するための手段】
【０００９】
実楽器音の多和音を縦続接続くし形フィルタで処理する時、くし形フィルタの入出力振幅比だけでなく、入出力周波数比にも注目する点が本発明のポイントである。楽音の非定常性が縦続接続くし形フィルタの周波数特性のために増幅され、振幅の大小だけでは音高推定が不可能である。しかし、楽音の基本となる各周波数成分がくし形フィルタでほぼ除かれた後に、非定常成分の低周波成分が出力される。このため、くし形フィルタの入出力周波数比に注目し、その周波数比の大きな変化点を示すくし形フィルタから、対応する音名（Ｃ、Ｃ＃、・・・、Ｂなど）が推定され、さらに、そのくし形フィルタの次数を各オクターブに対応するように変化させ、その出力波形の変化（振幅変化か周波数変化）からその音のオクターブが推定され、多和音のうちの１音の音高が推定できる。
【００１０】
このようにして推定した１音に対応したくし形フィルタを縦続接続くし形フィルタの先頭に移動させ、再度、波形の変化、すなわち振幅比や周波数比に注目して残りの音名（オクターブ推定の処理を含めると音高）の一つを推定し、再度このくし形フィルタを縦続接続くし形フィルタの先頭に移動させる。この操作を繰り返して、最終的に最初に移動させたくし形フィルタでの波形の変化に大きな変化が生じた時点で音高推定操作が終了する。そして、移動させたくし形フィルタから多和音の各音高が推定される。
【００１１】
対象音域の最低オクターブ中に存在する１２音に対応するくし形フィルタを縦続接続し、その各くし形フィルタの入出力振幅比、周波数比の変化を測定し、その変化の大きいくし形フィルタを順次、縦続接続くし形フィルタの先頭に移動させる。なお、縦続接続するくし形フィルタの数は、対象音域の最低オクターブよりひとつ低いオクターブに対応するくし形フィルタを使用すれば、等価的に７個で済まされることも可能である。
【発明の効果】
【００１２】
くし形フィルタは、周波数軸上に等間隔の零点を有しており、単音による音楽信号に含まれる全調波成分を除去する性質を有する。この性質を利用し、各音高に対応するくし形フィルタの出力から零出力を検出することにより、音高の推定が可能となる。また、くし形フィルタは遅延器のみを用いて構成されるため、時間領域においてはサンプル点間の差分のみで処理が可能であり、演算量の削減および処理系の簡単化を図ることができる。
【発明を実施するための最良の形態】
【００１３】
１オクターブ１２音に対応するくし形フィルタを図１のように縦続接続したシステムについて考える。このシステムは、原理的には、入力ｘ（ｎ）に含まれる音を１音ずつ除去し、入力中の全ての調波成分を除去した段階で、零値を出力する。ここで、このフィルタの最終出力の振幅特性は図２のようになる。原理的には、入力和音中の全ての音が消去された段階でくし形フィルタの出力は零となるが、実際の楽音では非定常性、周波数変動のため、完全な零出力とならず、むしろ図２の特性により、低域の周波数成分が増大し、振幅の大小での音高推定は困難になる。
【００１４】
そこで、振幅値に注目するのではなく、時間波形（周波数）の変化に注目した方法を提案する。縦続接続くし形フィルタに信号を入力し、和音の周波数成分が全て除去されると、くし形フィルタの出力は零出力にはならないが、その時点で時間波形は楽音の非定常性により、周波数変動分の低い周波数の出力波形になる。この時間波形の変化に注目し、和音の音高推定を行う。
【００１５】
例として、図３に５和音（Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）を入力した時の縦続接続くし形フィルタの各フィルタの出力を示す。各フィルタの出力波形を見ると、ＢからＤ＃までのフィルタの出力波形に対し、音高Ｄのフィルタまで通した時点で時間波形（周波数）は、低い周波数の出力波形になっているのがわかる。この変化により、入力和音にＤが含まれていると推定でき、ＢからＤまでの間の音高の組み合わせで入力和音が構成されているとも考えられる。
【００１６】
以上のことより、提案する音高推定アルゴリズムは図４のようになる。
【００１７】
（くし形フィルタによる音高推定の基本原理）
音高ｐの基本周波数をｆｐとした場合、音楽信号は基本周波数成分と、その整数倍の高調波成分から構成される。ここで、ｐ＝１（Ｃ）、２（Ｃ＃）、・・・、１２（Ｂ）である。この音楽信号を消去するくし形フィルタの伝達関数Ｈｐ（ｚ）は、数１で表され、その出力ｙｐ（ｎ）は、数２に示す減算で実現される。
【００１８】
【数１】
【００１９】
【数２】
【００２０】
ここで、Ｎｐは遅延数であり、サンプリング周波数をｆｓとするとＮｐは、数３で求められる。
【００２１】
【数３】
【００２２】
図５にオクターブ４の音名Ｅの音楽信号（基本周波数は３２９．６３Ｈｚ）の周波数特性と、対応するくし形フィルタの振幅特性を示す。くし形フィルタは、周波数領域において等間隔に零点が存在するような特性を有しており、音楽信号を構成する複数の周波数成分を同時に除去することができる。そのため、くし形フィルタの出力値から零出力を検出することで、音高の推定が可能となる。
【実施例】
【００２３】
対象音域をオクターブ３以上としたときは、オクターブ３の各音名に対応するくし形フィルタを図１のように１２個縦続接続し、このシステムに多和音を入力させる。この縦続接続くし形フィルタの振幅特性は、図２のようになり、オクターブ３以下の音域の振幅ゲインが大きくなっている。そのため、多和音の各周波数成分が除かれた後には、理想的にはくし形フィルタの出力振幅は零になるはずであるが、実楽器の多和音の非定常性成分（低周波成分）が出力されることになる。その出力振幅は後段のくし形フィルタになるほど大きくなり、入出力振幅比での音名推定は困難になる。しかし、くし形フィルタの入出力周波数比に注目することで、多和音がすべて除かれたかどうかを判断できる。
【００２４】
入出力周波数比に注目して多和音の各音名を推定する例を図６に示す。図６では、左側が縦続接続くし形フィルタの移動を示し、対応する右側は和音の全周波数成分が消去される前後の時間波形をｉｎｐｕｔとｏｕｔｐｕｔとしてそれぞれ示している。多和音の例として、５和音（オクターブ４：Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）を考える。縦続接続くし形フィルタ（オクターブ４の１２音に対応するくし形フィルタから構成）は、音名の高い順に並べている。より確実な推定をするためには、並べ方を逆にした二つの縦続接続くし形フィルタを使用するとよい。
【００２５】
最初に入出力周波数比が大きく変わるのは、音名Ｄを除くくし形フィルタである。なお、入出力振幅比(出力/入力)も１以下になっている。このＤに対応するくし形フィルタを縦続接続くし形フィルタの先頭に移動させ、２回目の入出力周波数比が大きく変わるくし形フィルタを見つける。今度は音名Ｆを除くくし形フィルタで検出される。このＦのくし形フィルタを先頭に移動させる。３回目の入出力周波数比の大きな変化は、音名Ｇに対応するくし形フィルタである。このＧのくし形フィルタを先頭に移動させる。４回目の入出力周波数比の大きな変化点は、音名Ａのくし形フィルタである。Ａのくし形フィルタを先頭に移す。そして、５回目の入出力周波数比の大きな変化するくし形フィルタは音名Ｂのくし形フィルタである。このＢに対応するくし形フィルタを先頭に移動し、6回目の入出力周波数比の大きな変化するくし形フィルタを検出する。今度は音名Ｄに相当するくし形フィルタであり、これは最初に先頭に移動されたくし形フィルタであるので、ここで音名推定の操作は終了する。
【００２６】
このようにして、各操作で先頭に移動されたくし形フィルタに相当する音名が、入力多和音の各音名に対応する。すなわち、入力されたＤ、Ｆ、Ｇ、Ａ、Ｂが推定されている。
【００２７】
次に、これらの各音名のオクターブンの推定は、各音名が推定された時点で、そのくし形フィルタの次数を２分の１、４分の１、・・・と変化させて、入出力周波数比の大きな変化が生じなく、または入出力振幅比(出力/入力)が１以上になる、次数の手前の次数が対応するオクターブとして決定される。そして、音名とオクターブから音高が決定される。なお、高調波成分のため、入出力周波数比（ゼロクロス数、ピーク数）の判定が入出力波形そのものから難しい場合は、平均化フィルタや自己相関関数の処理で明確にできる。この例の場合は、オクターブ推定でくし形フィルタの次数を２分の１にしたときに入出力周波数比の大きな変化が生じなくなる。よって、これらのオクターブは４と決定される。
【産業上の利用可能性】
【００２８】
本発明により、和音の数に制限されずに多和音の各音高を推定することができ、本格的な採譜システムを実現できる。これにより、音楽配信ビジネス、カラオケ産業での楽曲検索、民族音楽の記録の効率化、即興音楽の記録化、作曲作業の能率化および音楽教育への利用など、その応用範囲は多岐にわたる。音楽配信ビジネス（通信カラオケ、携帯電話やＰＨＳの着信メロディ）では、新曲をいち早く楽譜化することが必要であり、これを人手ではなく計算機による自動採譜システムで行なうことが熱望されている。この自動採譜システムとして４和音以上の多和音の精度良い音高推定が可能な本格的な採譜システムが実現されれば、この業界に多大な貢献をすることが可能である。
【００２９】
また、楽音のリアルタイム認識によるセッションシステム（相手の演奏の理解に基づいた演奏応答システム）のようなインタラクティブ音楽コンテンツの実現、さらには和音の組み合わせの多さを利用した個人認証システム（例えば、特許文献１を参照）への応用など考えられる。携帯電話の音源を利用した個人認証システムでは、採譜システムによる音高推定能力により、その性能（秘話性）が決定され、本発明による本格的な採譜システムは、そのための有望な技術になり、これに関係する産業界、ＩＴ社会での本発明の使用が期待される。
【特許文献１】特開２００１−３１８８９４号公報
【図面の簡単な説明】
【００３０】
【図１】くし形フィルタ（オクターブ３）の１２段縦続接続によるシステム構成
【図２】Ｈｐ（ｚ）（ｐは１から１２）を１２段接続した縦続接続くし形フィルタの振幅特性
【図３】５和音（Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）を入力したときの縦続接続くし形フィルタの各出力
【図４】図１のように構成されるくし形フィルタシステムの処理アルゴリズム
【図５】オクターブ４の音名Ｅの周波数特性とくし形フィルタの振幅特性
【図６】本発明の実施例に係る多和音（オクターブ４：Ｄ＋Ｆ＋Ｇ＋Ａ＋Ｂ）の音名推定過程
【符号の説明】
【００３１】
１ Ｈｐ（ｚ）＝１−ｚ−Ｎｐの特性を持つくし形フィルタ
２ くし形フィルタを１２段縦続接続構成にしたくし形フィルタシステム
【特許請求の範囲】
【請求項１】
対象音域の最低オクターブにある１２音の各音を消去するくし形フィルタを縦続接続し、この縦続接続くし形フィルタに入力音として多和音を入力させ、各くし形フィルタの出力の波形の変化、つまり入出力振幅比か入出力周波数比の変化に注目し、多和音の一つの音名を検出する。その変化したくし形フィルタを縦続接続くし形フィルタの先頭に移動させ、再度各くし形フィルタの各出力の波形変化を検出して多和音の残りの音名の一つを推定する。この動作を繰り返し、先頭に移動させた最初のくし形フィルタでの出力波形変化が生じるまでこの操作を繰り返し、多和音のすべての音名を推定する方法。
【請求項２】
請求項１に記述された方法により推定された音名に対応するくし形フィルタの次数（遅延数）を一つ上のオクターブに相当するように２分の１、更に一つ上に相当するように４分の１と順次変化させ、その出力波形の変化の有無により、その音のオクターブを推定し、音名とオクターブより音高を決定し、多和音の音高を推定する方法。
【請求項１】
対象音域の最低オクターブにある１２音の各音を消去するくし形フィルタを縦続接続し、この縦続接続くし形フィルタに入力音として多和音を入力させ、各くし形フィルタの出力の波形の変化、つまり入出力振幅比か入出力周波数比の変化に注目し、多和音の一つの音名を検出する。その変化したくし形フィルタを縦続接続くし形フィルタの先頭に移動させ、再度各くし形フィルタの各出力の波形変化を検出して多和音の残りの音名の一つを推定する。この動作を繰り返し、先頭に移動させた最初のくし形フィルタでの出力波形変化が生じるまでこの操作を繰り返し、多和音のすべての音名を推定する方法。
【請求項２】
請求項１に記述された方法により推定された音名に対応するくし形フィルタの次数（遅延数）を一つ上のオクターブに相当するように２分の１、更に一つ上に相当するように４分の１と順次変化させ、その出力波形の変化の有無により、その音のオクターブを推定し、音名とオクターブより音高を決定し、多和音の音高を推定する方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図２】
【図３】
【図４】
【図５】
【図６】
【公開番号】特開２００８−５２０２３（Ｐ２００８−５２０２３Ａ）
【公開日】平成２０年３月６日（２００８．３．６）
【国際特許分類】
【出願番号】特願２００６−２２７７７３（Ｐ２００６−２２７７７３）
【出願日】平成１８年８月２４日（２００６．８．２４）
【出願人】（３０４０２７３４９）国立大学法人豊橋技術科学大学 (391)
【Ｆターム（参考）】
【公開日】平成２０年３月６日（２００８．３．６）
【国際特許分類】
【出願日】平成１８年８月２４日（２００６．８．２４）
【出願人】（３０４０２７３４９）国立大学法人豊橋技術科学大学 (391)
【Ｆターム（参考）】
[ Back to top ]