(57)【要約】 【課題】既存の楽曲の音波形から所望の態様で音が変化している区間の音波形を素材として抽出し、新たな楽曲を作成することを可能にする。【解決手段】既存の楽曲の音波形をフレーム単位で解析し、各フレームの音波形からそのフレームの音響的な特徴を示す特徴量を算出する特徴量算出手段と、連続する所定数のフレームからなる区間毎に当該区間内での前記特徴量の変化の大きさおよび前記特徴量の変化の滑らかさに応じた指標値を算出する指標値算出手段と、前記指標値が利用者により指定された値となっている区間の音波形を前記楽曲の音波形から抽出する音波形抽出手段とを有することを特徴とする音波形抽出装置を提供する。【選択図】図２
2010217408.tif 000002
【特許請求の範囲】
【請求項１】
楽曲の音波形をフレーム単位で解析し、各フレームの音波形からそのフレームの音響的な特徴を示す特徴量を算出する特徴量算出手段と、
連続する所定数のフレームからなる区間毎に当該区間内での前記特徴量の変化の大きさおよび前記特徴量の変化の滑らかさに応じた指標値を算出する指標値算出手段と、
前記指標値が予め設定された値となっている区間の音波形を前記楽曲の音波形から抽出する音波形抽出手段と
を有することを特徴とする音波形抽出装置。
【請求項２】
前記指標値算出手段は、
前記連続する所定数のフレームからなる区間の末尾のフレームにおける前記特徴量と当該区間の先頭のフレームにおける前記特徴量の差を、当該区間における前記特徴量の最大値と最小値の差、または当該差をより大きな値に補正した補正値で規格化して当該区間についての前記指標値を算出する
ことを特徴とする請求項１に記載の音波形抽出装置。
【請求項３】
前記指標値算出手段は、
前記連続する所定数のフレームからなる区間における前記特徴量のばらつきが予め定められた第１の閾値よりも大きい場合、または当該区間における前記特徴量の最大値と最小値の差が予め定められた第２の閾値よりも小さい場合に、前記最大値と最小値の差をより大きな値に補正する
ことを特徴とする請求項２に記載の音波形抽出装置。
【請求項４】
コンピュータを、
楽曲の音波形をフレーム単位で解析し、各フレームの音波形からそのフレームの音響的な特徴を示す特徴量を算出する特徴量算出手段と、
連続する所定数のフレームからなる区間毎に当該区間内での前記特徴量の変化の大きさおよび前記特徴量の変化の滑らかさに応じた指標値を算出する指標値算出手段と、
前記指標値が予め設定された値となっている区間の音波形を前記楽曲の音波形から抽出する音波形抽出手段
として機能させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、既存の楽曲の音波形から新たな楽曲の素材となる音波形を抽出する技術に関する。
【背景技術】
【０００２】
既存の楽曲から抽出した音波形をパッチワークのように繋げて新たな楽曲を作成することが一般に行われている。例えば、特許文献１には、複数の楽曲の各々を拍に同期した時点で区切って得られる各音波形を素材として新たな楽曲を生成する技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００８−１２９１３５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
このようにパッチワークのように素材（音波形）を繋げて新たな楽曲を作成する場合、例えば連続的に音が変化しているなど所望の態様で音が変化している素材を選んで使用したいといったニーズがあった。しかし、従来は各素材の音を実際に聴いてみなければ音がどのように変化しているのかを把握することができず、膨大な数の素材のうちから所望の態様で音が変化しているものを探し出すことは略不可能であった。
本発明は上記課題に鑑みて為されたものであり、既存の楽曲の音波形から所望の態様で音が変化している区間の音波形を素材として抽出し、新たな楽曲の作成に利用することを可能にする技術を提供することを目的する。
【課題を解決するための手段】
【０００５】
上記課題を解決するため本発明は、楽曲の音波形をフレーム単位で解析し、各フレームの音波形からそのフレームの音響的な特徴を示す特徴量を算出する特徴量算出手段と、連続する所定数のフレームからなる区間毎に当該区間内での前記特徴量の変化の大きさおよび前記特徴量の変化の滑らかさに応じた指標値を算出する指標値算出手段と、前記指標値が予め設定された値（目標値）となっている区間の音波形を前記楽曲の音波形から抽出する音波形抽出手段とを有することを特徴とする音波形抽出装置、およびコンピュータを上記各手段として機能させることを特徴とするプログラムを提供する。
【０００６】
本発明によれば、既存の楽曲の音波形から上記指標値算出手段により算出される指標値が予め設定された目標値となっている区間の音波形が新たな楽曲の素材として抽出される。ここで、上記指標値は、既存の楽曲の音波形についてフレーム毎に算出される特徴量の変化の大きさや変化の滑らかさを示す値である。例えば、上記特徴量としてピッチを用いるようにすれば上記指標値は所定数のフレーム内でのピッチの変化の大きさおよび滑らかさを示すため、上記目標値に応じてピッチが滑らかに上昇している区間の音波形や、逆に、ピッチが滑らかに下降している区間の音波形が新たな楽曲の素材として抽出される。また、音量を上記特徴量として用いると、上記指標値は所定数のフレーム内での音量の変化の大きさおよび変化の滑らかさを示すため、上記目標値に応じて音量が滑らかに増加してゆく区間の音波形や、逆に、音量が滑らかに減少してゆく区間の音波形が新たな楽曲の素材として抽出される。
【０００７】
上記指標値算出手段による指標値の算出態様としては種々のものが考えられる。例えば、連続する所定数のフレームからなる区間の末尾のフレームにおける特徴量と当該区間の先頭のフレームにおける特徴量の差を当該区間における特徴量の最大値と最小値の差で規格化（除算）して当該区間についての指標値を算出する態様や、当該区間における特徴量の最大値と最小値の差をより大きな値に補正し、その補正後の値で当該区間の末尾のフレームにおける特徴量と当該区間の先頭のフレームにおける特徴量の差を規格化して当該区間についての指標値を算出する態様が考えられる。また、後者の態様においては、連続する所定数のフレームからなる区間における特徴量のばらつきの大きさ（例えば、上記所定数のフレーム内での上記変化量の分散や標準偏差）が予め定められた第１の閾値よりも大きい場合、または当該区間における特徴量の最大値と最小値の差が予め定められた第２の閾値よりも小さい場合に、当該区間における特徴量の最大値と最小値の差をより大きな値に補正して上記指標値を算出するとしても良い。また、これらの態様の他に、連続する所定数のフレームの先頭のものから順に、互いに隣接するフレーム間の特徴量の変化量を算出し、その変化量の平均値をその変化量のばらつきの大きさで補正した後に、その補正値を所定の範囲（例えば、−１から１の範囲など）に収まる値に規格化して上記指標値とする態様も考えられる。
【図面の簡単な説明】
【０００８】
【図１】この発明の一実施形態である音波形抽出装置１００の構成例を示す図である。
【図２】同音波形抽出装置１００の制御部１１０が音波形抽出プログラムにしたがって実行する処理の流れを示すフローチャートである。
【図３】同制御部１１０が実行する特徴量算出処理ＳＡ１００の一例を説明するための図である。
【図４】同制御部１１０が実行するおよび指標値算出処理ＳＡ１１０の一例を説明するための図である。
【図５】複数のフレームに亘る特徴量の変化態様と指標値の時間変化の関係を示す図である。
【図６】音波形抽出処理ＳＡ１２０にて制御部１１０が表示部１４０に表示させるユーザインタフェース画面の一例を示す図である。
【発明を実施するための形態】
【０００９】
以下、本発明の実施形態について図面を参照しつつ説明する。
（Ａ：構成）
図１は、本発明の一実施形態に係る音波形抽出装置１００の構成例を示すブロック図である。この音波形抽出装置１００は、既存の楽曲の音波形から新たな楽曲の素材となる音波形を抽出する装置である。この音波形抽出装置１００の利用者は、既存の楽曲中で音響的な特徴量（例えば、ピッチや音量など）が滑らかに変化している区間の音波形を新たな楽曲の素材として抽出することができる。図１に示すように、音波形抽出装置１００は、制御部１１０、インタフェース群１２０、操作部１３０、表示部１４０、記憶部１５０、および、これら構成要素間のデータ授受を仲介するバス１６０を有している。
【００１０】
制御部１１０は、例えばＣＰＵ（Central Processing Unit）である。制御部１１０は、記憶部１５０に格納されている音波形抽出プログラムを実行することにより、音波形抽出装置１００の制御中枢の役割を果たす。この音波形抽出プログラムにしたがって制御部１１０が実行する処理については後に詳細に説明する。
【００１１】
インタフェース群１２０は、ネットワークを介して他の装置との間でデータ通信を行うためのネットワークインタフェースや、ＣＤ−ＲＯＭ（Compact Disk-Read Only Memory）などの外部記録媒体との間でデータの授受を行うためのドライバを含んでいる。本実施形態では、インタフェース群１２０のうちの適切なものを介して楽曲の音波形を示す楽曲データ（楽曲の波形を所定のサンプリング周期でサンプリングしたサンプル列）が入力される。本実施形態では、インタフェース群１２０を介して入力される楽曲データの表す音波形が新たな楽曲の素材となる音波形の抽出元となるのである。なお、上記楽曲データとしてアナログ形式のオーディオ信号が与えられる場合には、Ａ／Ｄ変換器を介して入力するようにすれば良い。
【００１２】
操作部１３０は、例えばマウスなどのポンティングデバイスや複数のキーを備えたキーボードなどであり、各種の入力操作を利用者に行わせるためのものである。操作部１３０は、利用者により為された操作を示すデータを制御部１１０に与える。これにより操作部１３０に対して利用者が行った操作の内容が制御部１１０に伝達される。表示部１４０は、例えば液晶ディスプレイとその駆動回路である。この表示部１４０には、音波形抽出装置１００の利用を促すためのユーザインタフェース画面（図６参照）が表示される。
【００１３】
記憶部１５０は、図１に示すように、揮発性記憶部１５０ａと不揮発性記憶部１５０ｂを含んでいる。揮発性記憶部１５０ａは、例えばＲＡＭ（Random Access Memory）であり、各種プログラムを実行するためのワークエリアとして利用される。一方、不揮発性記憶部１５０ｂは、例えばハードディスクである。不揮発性記憶部１５０ｂには、前述した音波形抽出プログラムが格納されている他、新たな楽曲の素材として抽出した音波形の波形データの集合体である素材データベース（図示略）が格納される。図２は、音波形抽出プログラムにしたがって制御部１１０が実行する処理の流れを示すフローチャートである。図２に示すように、音波形抽出プログラムにしたがって制御部１１０が実行する処理には、特徴量算出処理ＳＡ１００、指標値算出処理ＳＡ１１０および音波形抽出処理ＳＡ１２０の３つの処理が含まれる。これら３つの処理の詳細については後に明らかにするが、その概要は以下の通りである。
【００１４】
図２の特徴量算出処理ＳＡ１００は、インタフェース群１２０を介して入力される楽曲データの表す音波形を所定時間長（或いは、所定サンプル数）のフレーム単位で解析し、各フレームの音波形からそのフレームの音響的な特徴を示す特徴量を算出する処理である。前述したように、この特徴量算出処理ＳＡ１００にてフレーム毎に算出する特徴量の一例としてはピッチや音量が挙げられる。この特徴量算出処理ＳＡ１００にてどのような特徴量を算出するのかについては、音波形抽出装置１００の工場出荷時点で予め定めておいても良く、また、予め定められた複数種の特徴量のうちから音波形抽出装置１００の使用のたびに利用者に選択させるようにしても良い。
【００１５】
指標値算出処理ＳＡ１１０は、連続する所定数のフレーム（例えば、５フレームや８フレームなど）からなる区間毎に、当該区間内での特徴量の変化の大きさおよび変化の滑らかさに応じた指標値を算出する処理である。そして、音波形抽出処理ＳＡ１２０は、指標値算出処理ＳＡ１１０にて算出された指標値が利用者により設定された値となっている区間の音波形を上記楽曲データの表す音波形から抽出する処理である。
以上が音波形抽出装置１００の構成である。
【００１６】
（Ｂ：動作）
以下、特徴量算出処理ＳＡ１００にて算出する特徴量として、新たな楽曲の素材となる音波形の抽出元の楽曲を構成する複数のパートのうちの特定のパートのピッチを採用した場合を例にとって音波形抽出装置１００が実行する動作を説明する。ここで、複数のパートのうちから特定のパートを指定させる態様の具体例としては、例えば、そのパートの担当楽器の演奏音の周波数帯域を示す値（例えば、上限周波数および下限周波数）で指定させる態様が考えられる。なお、インタフェース群１２０を介して入力される楽曲データの表す音には、各パートの演奏音が含まれているため、上記のように周波数範囲を指定するだけでは上記特定のパートの音のみを抽出することは一般にはできないのであるが、本実施形態では上記特定のパートのピッチ変化を高い確度で追跡するための工夫が為されている。
【００１７】
図３は、特徴量算出処理ＳＡ１００の流れを説明するための図である。この特徴量算出処理ＳＡ１００では、まず、インタフェース群１２０を介して入力された楽曲データを所定の時間長（サンプル数）ずつのフレームに区切り、フレーム毎に高速フーリエ変換（Fast Fourier Transform：以下、ＦＦＴ）を施して振幅スペクトラムを算出する処理が実行される。例えば、図３では、楽曲の先頭に位置するフレームを「０番目のフレーム」とした場合における、ｎ番目のフレームについてＦＦＴを施し振幅スペクトラムを算出する場合について例示されている。
【００１８】
次いで、特徴量算出処理ＳＡ１００では、各フレームの振幅スペクトラムを解析して特徴量を算出し、フレーム毎に算出された特徴量を書き込んだ特徴量テーブルを生成する処理が実行される。前述したように本実施形態では、特徴量として特定のパートのピッチが採用されるため、各フレームの振幅スペクトラムから当該パートに属すると推定される音のピッチが抽出され、そのピッチを示すデータを格納した特徴量テーブルが生成される。
【００１９】
図３には、特徴量算出処理ＳＡ１００で生成される特徴量テーブルの一例が示されている。この特徴量テーブルには、フレーム番号が若いものから順に、そのフレームの振幅スペクトラムの特定の周波数帯域（上記特定のパートに対応する周波数帯域）に現れているローカルピークのうちで振幅値が上位Ｎ（本実施形態では、５）番目までのＮ個のローカルピークの中心周波数が格納される。図３では、ｎ番目のフレームの振幅スペクトラムの特定の周波数帯域でｉ番目に振幅値が大きいローカルピークの中心周波数がｆｎｉで表されている。ここで、ローカルピークとは振幅スペクトラムを表すグラフにて振幅値が増加から減少に転じている箇所のことである。
【００２０】
ここで、各フレームの特定の周波数帯域に現れているローカルピークのうち、上位５番目までの振幅値を有するもののみに着目するのは、他のパートの演奏音や雑音により現れているローカルピークを排除するためである。これは、上記特定のパートの担当楽器ほど、そのパートに対応する周波数帯域に属する音を高い音量で放音していると考えられるからである。なお、図３に示す特徴量テーブルを生成するために、何番目までのローカルピークに着目するのかについては他のパートの演奏音や雑音の大きさとの兼ね合いで定めるようにすれば良い。また、本実施形態では、各フレームの特定の周波数帯域に現れているローカルピークのみに着目したが、振幅スペクトラムに現れているすべてローカルピークのうちで上位５番目までの振幅値を有するものに着目しても勿論良い（この場合は、楽曲全体でピッチが滑らかに変化している区間を新たな楽曲の素材として切り出すこととなる）。
以上が特徴量算出処理ＳＡ１００の内容である。
【００２１】
次いで、指標値算出処理ＳＡ１１０について説明する。
指標値算出処理ＳＡ１１０は、特徴量算出処理ＳＡ１００にて算出した特徴量の時間変化を追跡する処理と、その追跡結果から特徴量の時間変化の態様を示す指標値を算出する処理とに大別される。図４は前者の処理内容を説明するための図であり、図５は後者の処理にて算出される指標値について説明するための図である。
【００２２】
指標値算出処理ＳＡ１１０では、まず、特徴量算出処理ＳＡ１００にて生成した特徴量テーブル（図３参照）の格納内容から、特徴量の時間変化を表す特徴量変化追跡テーブル（図４参照）が以下の要領で生成される。すなわち、制御部１１０は、まず、フレーム番号が若いものから順に、特徴量テーブルに格納されている５つの中心周波数のうちで振幅値が最も大きいローカルピークの中心周波数を予め定められた閾値ｆｔｈ１と比較し、当該中心周波数がはじめて当該閾値ｆｔｈ１を下回るフレームのフレーム番号ＮＡを特定する。ここで、特徴量テーブルに格納されている５つの中心周波数のうちで振幅値が最も大きいローカルピークの中心周波数に着目するのは、振幅値が最も大きいローカルピークは、利用者に指定されたパートの担当楽器の演奏音に対応している可能性が高いと考えられるからである。また、閾値ｆｔｈ１の具体的な値については、利用者により指定されたパートの担当楽器の演奏音の周波数帯域との関係で定めれば良い（例えば、基本周波数成分および倍音成分が適宜含まれるように定めれば良い）。
【００２３】
上記のようにしてフレーム番号ＮＡが特定されると、制御部１１０は、そのフレームから予め定められた範囲（例えば、１０フレームの範囲）で、振幅値の最も大きいローカルピークの中心周波数が最も低くなっているもの（但し、当該中心周波数が閾値ｆｔｈ２（ｆｔｈ２＜ｆｔｈ１：本実施形態では、ｆｔｈ２＝０）を上回っていることが条件）のフレーム番号ＮＢを特定する。前述したように、振幅値が最も大きいローカルピークは、利用者により指定されたパートの担当楽器の演奏音に対応している可能性が高いと考えられる。したがって、上記の要領でフレーム番号ＮＢを特定することで、上記担当楽器が基本周波数の演奏音を奏でている区間が特定されると考えられる。なお、本実施形態では閾値ｆｔｈ２＝０としたが、この閾値ｆｔｈ２についても利用者により指定されたパートの担当楽器の演奏音の周波数帯域との関係で定めるようにしても良い。また、本実施形態のように閾値ｆｔｈ２を０とする態様においては、特徴量テーブル（図３参照）に格納されているｆｋ１（ｋ：フレーム番号）のうちで最小のもののフレーム番号ｋを上記ＮＢとしても勿論良い。
【００２４】
次いで、制御部１１０は、フレーム番号ｋが０〜ＮＡ−１までのフレームについては、そのフレーム番号ｋと、特徴量テーブル（図３参照）に格納されている中心周波数ｆｋ１とを対応付けて特徴量変化追跡テーブルに書き込む。フレーム番号ＮＡのフレームについては、制御部１１０は、そのフレーム番号ＮＡと周波数ｆＮＢ１（すなわち、フレーム番号ＮＢに関して特徴量テーブルに格納されているローカルピークの中心周波数のうち、振幅が最大のローカルピークの中心周波数）とを対応付けて特徴量変化追跡テーブルに書き込む。フレーム番号ｋがＮＡ＋１以降のフレームについては、制御部１１０は以下の要領で特徴量変化追跡テーブルに書き込む周波数を選択する。すなわち、制御部１１０は、まず、フレーム番号ｋ（ｋ≧ＮＡ＋１）のフレームに関して特徴量テーブルに格納されている５種類の周波数ｆｋｉ（ｉ＝１〜５）のうち、一つ手前のフレーム（すなわち、フレーム番号がｋ−１のフレーム）について特徴量変化追跡テーブルに書き込んだ周波数に最も近いものを選択する。そして、制御部１１０は、このようにして選択した周波数をフレーム番号ｋに対応付けて特徴量変化追跡テーブルに書き込む。これにより、図４に示す特徴量変化追跡テーブルが得られるのである。前述したように、フレーム番号ＮＡに対応付けて特徴量変化追跡テーブルに格納されている周波数ｆＮＢ１は、利用者により指定されたパートの担当楽器の基本周波数と考えられる。したがって、上記のようにして特徴量変化追跡テーブルを生成することによって、フレーム番号がＮＡ＋１以降のフレームについては、利用者により指定されたパートの基本周波数の時間変化が追跡されていると期待される。
【００２５】
次いで、制御部１１０は、５フレームや８フレームなど連続する所定数のフレームからなる区間毎に、当該区間内での特徴量の変化態様（増加しているのか、減少しているのか、また、その変化は滑らかであるのか否かなど）を示す指標値Ｘを上記特徴量変化追跡テーブルの格納内容を利用して算出する。ここで指標値Ｘの算出態様としては、種々の態様が考えられる。例えば、上記区間を構成する所定数のフレームの先頭のものから順に、互いに隣接するフレーム間の特徴量の変化量（すなわち、ｋ＋１番目のフレームの特徴量とｋ番目のフレームの特徴量との差）を算出し、その変化量の平均値をその変化量のばらつきの大きさ（例えば、分散や標準偏差）で補正した後に、その補正後の値を−１から１などの所定の範囲に収まるように規格化して指標値Ｘとする態様が考えられる。本実施形態では、以下の式（１）にしたがって上記指標値Ｘを算出する。
Ｘ＝（Ｅｎｄ−Ｂｅｇｉｎ）／（ＭＡＸ−ＭＩＮ＋ｆａｓｔＰ×ＶＡＲ＋ｓｌｏｗＰ）・・・（１）
【００２６】
式（１）の右辺の構成要素の意味は夫々以下の通りである。
Ｅｎｄは上記所定数フレームからなる区間の末尾のフレームにおける特徴量（特徴量変化追跡テーブルに書き込まれている中心周波数）であり、Ｂｅｇｉｎは同区間の先頭のフレームにおける特徴量である。ＭＡＸは上記区間における特徴量の最大値であり、ＭＩＮは同区間における特徴量の最小値である。ＶＡＲは上記区間において互いに隣接するフレーム間の特徴量の変化量のばらつきの大きさを表す値（本実施形態では、互いに隣接するフレーム間の特徴量の差分についての分散）である。そして、ｆａｓｔＰは、上記区間内での特徴量の変化のばらつきが大きすぎる場合に課されるペナルティを表す値であり、ｓｌｏｗＰは上記区間内での特徴量の変化が緩やかすぎる場合に課されるペナルティを表す値である。これらｆａｓｔＰおよびｓｌｏｗＰは、何れも正の固定値である。これらペナルティに応じて、上記区間における特徴量の最大値と最小値の差はより大きな値に補正される。つまり、式（１）は、上記区間における特徴量の最大値と最小値の差をより大きな値に補正し、その補正後の値で、同区間の末尾のフレームにおける特徴量と先頭のフレームにおける特徴量の差を規格化（除算）して指標値Ｘを算出することを意味するのである。このように特徴量の変化のばらつきや変化のゆるやかさに応じたペナルティを課す理由については後に明らかにする。
【００２７】
例えば、図５（Ａ）に示すように、特徴量が単調に増加している場合、ＭＡＸ＝Ｅｎｄ、ＭＩＮ＝Ｂｅｇｉｎであり、かつＶＡＲは略ゼロとなるから、式（１）にしたがって算出される指標値Ｘは、正の値（ＭＡＸとＭＩＮの差がｓｌｏｗＰに比べて大きければ大きいほど１に近い値）となる。逆に、図５（Ｂ）に示すように、特徴量が単調に減少している場合、式（１）にしたがって算出される指標値Ｘは、負の値（ＭＡＸとＭＩＮの差がｓｌｏｗＰに比べて大きければ大きいほど−１に近い値）となる。これに対して、図５（Ｃ）に示すように、上記所定数フレームからなる区間内の前半では特徴量が次第に減少し、その後増加に転じてＥｎｄとＢｅｇｉｎが略等しくなっている場合（すなわち、ＥｎｄとＢｅｇｉｎの差が略ゼロの場合）には、式（１）にしたがって算出される指標値Ｘは略ゼロになる。つまり、式（１）にしたがって算出される指標値Ｘは、その値が１に近いほど特徴量が滑らかに増加していることを表し、その値が−１に近いほど特徴量が滑らかに減少していることを表し、その値が０に近い場合には上記所定数のフレームに亘って特徴量があまり変化していないか、または特徴量の変化が滑らかではないことを表すのである。
【００２８】
ここで、特徴量の変化のばらつきが大きすぎる場合や特徴量の変化が緩やかすぎる場合にペナルティを課す理由は以下の通りである。特徴量の変化のばらつきの大きさや特徴量の変化の緩やかさを考慮せずに常に（Ｅｎｄ−Ｂｅｇｉｎ）を（ＭＡＸ−ＭＩＮ）で除算（規格化）して指標値Ｘを算出するとしたならば、図５（Ｅ）に示すように所定数のフレームの末尾の部分で特徴量が急激に変化している場合や図５（Ｆ）に示すように特徴量の変化が小さい場合についても指標値Ｘの値は１になり、図５（Ａ）に示す変化態様と区別がつかなくなる。そこで、本実施形態では、図５（Ｅ）に示す変化態様や図５（Ｆ）に示す変化態様と図５（Ａ）に示す変化態様とを指標値Ｘの値で区別することができるようにするために特徴量の変化のばらつきの大きさや特徴量の変化の緩やかさに応じてぺナルティを課すこととしたのである。また、図５（Ｆ）に示す変化態様では、（ＭＡＸ−ＭＩＮ）は０に近い値となり、ＶＡＲも略０となるから、ｓｌｏｗＰに正の固定値を設定しないとゼロ除算の危険もあるが、ｓｌｏｗＰに正の固定値を設定することでゼロ除算の危険を回避することもできる。
【００２９】
なお、本実施形態では、ｆａｓｔＰおよびｓｌｏｗＰとして正の固定値を用いたが、所定数フレームからなる区間における特徴量のばらつきの大きさが所定の閾値ｔｈ１を超えた場合にｆａｓｔＰに正の固定値をセットし、その他の場合（すなわち、上記ばらつきの大きさが閾値ｔｈ１以下である場合）には、ｆａｓｔＰにゼロをセットするとしても良い。同様に、所定数フレームからなる区間における特徴量の最大値と最小値との差が所定の閾値ｔｈ２より小さい場合にｓｌｏｗＰに正の固定値をセットし、その他の場合（すなわち、上記差が閾値ｔｈ２以上である場合）にはｓｌｏｗＰにゼロをセットするとしても良い。また、本実施形態では、ＭＡＸとＭＩＮの差に特徴量の変化のばらつきの大きさや特徴量の変化の緩やかさに応じた値を加算してその補正を行ったが、特徴量の変化のばらつきの大きさや特徴量の変化の緩やかさに応じた値を乗算することで補正を行っても勿論良い。
【００３０】
このように、本実施形態では、特徴量の変化のばらつきの大きさと特徴量の変化の緩やかさの各々に応じてペナルティを課したが、何れか一方についてのみペナルティを課しても良い。例えば、特徴量の変化のばらつきの大きさに関してだけペナルティを課す場合には、式（１）において常にｓｌｏｗＰ＝０とすれば良く、逆に、特徴量の変化の緩やかさに関してだけペナルティを課す場合には常にｆａｓｔＰ＝０とすれば良い。また、上記所定数フレームの区間内で特徴量が増加傾向にあるのか、減少傾向にあるのか、それともどうちらでもないのかを大まかに判別することができれば良い場合には、特徴量の変化のばらつきの大きさや特徴量の変化の緩やかさに応じたペナルティを課す必要はなく、（Ｅｎｄ−Ｂｅｇｉｎ）を（ＭＡＸ−ＭＩＮ）で除算（規格化）して指標値Ｘを算出すれば良い。
以上が指標値算出処理ＳＡ１１０の内容である。
【００３１】
次いで、音波形抽出処理ＳＡ１２０の内容について説明する。
この音波形抽出処理ＳＡ１２０では、制御部１１０は、まず、新たな楽曲の素材として切り出す区間を利用者に指定させるためのユーザインタフェース画面を表示部１４０に表示させる。図６は、表示部１４０に表示されるユーザインタフェース画面の一例を示す図である。図６に示すように、このユーザインタフェース画面には、２つの表示領域（表示領域Ａ０１およびＡ０２）と入力領域Ｉ０１と、仮想ボタンＢ０１とが配置されている。
【００３２】
図６の表示領域Ａ０１には、インタフェース群１２０を介して入力された楽曲データの表す音波形が表示され、表示領域Ａ０２には、指標値算出処理ＳＡ１１０にて所定数のフレーム毎に算出された指標値の時間変化を表すグラフが上記音波形と時間軸Ｔを共通にして表示される。音波形抽出装置１００の利用者は、表示領域Ａ０２に表示されるグラフを参考に上記楽曲データの表す楽曲の音波形から新たな楽曲の素材として抽出すべき音波形の区間を指標値で指定することができる。例えば、ピッチが滑らかに上昇している区間の抽出を所望する場合には、指標値が略１になっている区間を表示領域Ａ０２に表示されるグラフにて指定すれば良く、逆に、ピッチが滑らかに下降している区間の抽出を所望する場合には、同指標値が略−１になっている区間を指定すれば良い。なお、新たな楽曲の素材として抽出すべき音波形の区間を指標値で指定する際には、抽出を所望する区間の指標値についての目標値を操作部１３０の操作により数値で入力させるとしても勿論良い。
【００３３】
図６の入力領域Ｉ０１には、新たな楽曲の素材として切り出す音波形を一意に識別するための識別子（例えば、当該音波形に付与する名称を表す文字列）が操作部１３０に対する操作により入力される。このような識別子を入力させるようにしたのは、新たな楽曲の素材として複数の音波形を切り出す場合に、それらを容易に区別することができるようにするためである。本実施形態では、上記識別子を利用者に入力させるようにしたが、素材の抽出元の楽曲の名称を表す文字列とその音波形の開始時刻および終了時刻を表す文字列とを連結して上記識別子を自動生成するようにしても良い。
【００３４】
そして、図６の仮想ボタンＢ０１は、利用者により設定された指標値を有する区間の音波形の波形データと、その音波形を一意に識別する識別子とを対応付けて記憶する処理の実行を制御部１１０に対して指示するための仮想操作子である。仮想ボタンＢ０１がクリックされると、制御部１１０は、インタフェース群１２０を介して入力された楽曲データ（サンプル列）のうち、指標値算出処理ＳＡ１１０にて算出された指標値Ｘが利用者により設定された値となっている区間のサンプル列（すなわち、当該区間の音波形を示す波形データ）のコピーを生成し、当該コピーサンプル列に上記識別子を対応付けて不揮発性記憶部１５０ｂに書き込む。なお、本実施形態では、指標値Ｘが利用者により設定された値となっている区間の音波形を表す波形データとその音波形を一意に識別する識別子とを対応付けて不揮発性記憶部１５０ｂに格納することで新たな楽曲の素材となる音波形の抽出を実現した。しかし、上記楽曲データと時間軸を共通にするデータトラックであって、指標値Ｘが利用者により設定された値となっている区間の開始時刻および終了時刻を示す情報と上記識別子とを格納したデータトラックを生成し、このデータトラックを上記楽曲データと対応付けて不揮発性記憶部１５０ｂに格納することで新たな楽曲の素材となる音波形の抽出を実現しても良い。
以上が音波形抽出処理ＳＡ１２０の内容である。
【００３５】
以上説明したように本実施形態に係る音波形抽出装置１００によれば、ピッチが滑らかに変化している区間や音量が滑らかに変化している区間の音波形を既存の楽曲から抽出し、新たな楽曲の素材とすることが容易になる。また、上記のようにして切り出される波形データの集合（素材データベース）をインタフェース群１２０のうちの適切なものを介して通信回線経由で他の装置へ送信しても良く、また、インタフェース群１２０のうちの適切なものを介してＣＤ−ＲＯＭなどの記録媒体に書き込み、その記録媒体を配布しても良い。
【００３６】
（Ｃ：他の実施形態）
以上、本発明の一実施形態について説明したが、かかる実施形態に以下に述べる変形を加えても勿論良い。
（１）上述した実施形態では、新たな楽曲の素材として抽出する音波形を示す値（指標値についての目標値）を、ユーザインタフェース画面（図６参照）の表示内容を参考に利用者に指定させた。しかしながら、上記目標値を予め設定した状態で音波形抽出プログラムの実行を開始させるようにしても良く、このような態様においては、図６に示すようなユーザインタフェース画面の表示を制御部１１０に実行させる必要はない。また、上述した実施形態では、指標値Ｘが利用者により設定された値となっている区間の音波形を表す波形データとその音波形を一意に識別する識別子とを対応付けて不揮発性記憶部１５０ｂに格納することで素材データベースを生成した。しかし、指標値Ｘが利用者により設定された値となっている区間毎にその波形データとその指標値Ｘとを対応付けて不揮発性記憶部１５０ｂに格納して素材データベースを構成し、指標値Ｘをキーとしてその素材データベースから所望の態様で音が変化している音波形の波形データを読み出すようにしても勿論良い。
【００３７】
（２）上述した実施形態では、特徴量算出処理ＳＡ１００にて算出する特徴量としてピッチを採用したが、音量や他の種類の特徴量を採用しても勿論良い。例えば、音量を特徴量として採用する場合には、フレーム毎に音量（ラウドネス）を算出し、その算出値をフレーム番号と対応付けて揮発性記憶部１５０ａに書き込むことで特徴量変化追跡テーブルを生成すれば良い。このように音量を特徴量として選択すると、音量が滑らかに増加してゆく区間（すなわち、音がフェードインする区間）や音量が滑らかに減少してゆく区間（すなわち、音がフェードアウトする区間）を新たな楽曲の素材として切り出すことが可能になる。
【００３８】
また、ピッチや音量の他の特徴量としては、スペクトラルクレストやスペクトラルフラットネス、コードストレングスなどを用いることが考えられる。スペクトラルクレストとスペクトラルフラットネスは、スペクトラルエナジの平坦さを示す特徴量である。スペクトラルクレストは、フレーム内でのスペクトラルエナジの最大値を同フレーム内での平均値で除算して算出され、スペクトラルフラットネスは、フレーム内でのスペクトラルエナジの幾何平均（相乗平均）を算術平均（相加平均）で除算して算出される。スペクトラルエナジが平坦であるほどノイズのような音を表すと考えられる。したがって、特徴量算出処理ＳＡ１００にて算出する特徴量としてスペクトラルクレストやスペクトラルフラットネスを採用すると、次第にノイズのようになってゆく音などを新たな楽曲の素材として切り出すことが可能になる。また、コードストレングスは、コード感の有無（すなわち、調和的であるか否か）を示す特徴量である。特徴量算出処理ＳＡ１００にて算出する特徴量としてコードスレングスを採用すると、次第にコード感が薄れてゆく音や、逆に、次第にコード感が強くなってゆく音を新たな楽曲の素材として切り出すことが可能になる。さらに、複数種の特徴量から算術演算により求まる新たな特徴量、または１または複数種の特徴量から関数演算により求まる新たな特徴量を定義し、当該新たな特徴量の時間変化を追跡しても良いことは勿論である。
【００３９】
（３）上述した実施形態では、式（１）にしたがって算出される指標値Ｘの時間変化を示すグラフを楽曲の音波形と共通の時間軸で表示することで、そのグラフから特徴量が滑らかに増加（或いは減少）している区間を把握することができるようにした。しかし、所定数のフレーム内で特徴量がどの程度変化しているのかを示す数値（すなわち、特徴量の変化量）を表示させても良い。例えば、特徴量算出処理ＳＡ１００にて算出する特徴量が特定のパートのピッチである場合、上記変化量をセント値に換算して表示させるようにすれば良い。このようにピッチの変化量をセント値に換算して表示させるようにすると、指標値算出処理ＳＡ１１０にて基本周波数ではなく倍音の周波数の時間変化が追跡されている場合であっても、そのパートのピッチが概ね何オクターブ変化しているのかを把握することができる。
【００４０】
（４）上述した実施形態では、操作部１３０および表示部１４０を音波形抽出装置１００の構成要素としたが、これらは必ずしも必須ではない。例えばキーボードやマウスなどの入力装置をインタフェース群１２０のうちの適当なものを介して接続し、この入力装置に操作部１３０の役割を担わせても良い。同様に、液晶ディスプレイとその駆動回路とを含む表示装置をインタフェース群１２０のうちの適当なものを介して接続し、この表示装置に表示部１４０の役割を担わせても良い。
【００４１】
（５）上述した実施形態では、制御部１１０に、特徴量算出処理ＳＡ１００、指標値算出処理ＳＡ１１０および音波形抽出処理ＳＡ１２０を実行させる音波形抽出プログラムが不揮発性記憶部１５０ｂに予め格納されていた。しかし、ＣＤ−ＲＯＭなどのコンピュータ読み取り可能な記録媒体に上記音波形抽出プログラムを書き込んで配布しても良く、また、インターネットなどの電気通信回線経由のダウンロードにより上記音波形抽出プログラムを配布しても良い。このようにして配布される音波形抽出プログラムにしたがって、パーソナルコンピュータなどの一般的なコンピュータを作動させることで、そのコンピュータを音波形抽出装置１００として機能させることが可能になるからである。
【００４２】
（６）上述した実施形態では、本発明に係る音波形抽出装置１００の特徴である３つの処理（すなわち、特徴量算出処理ＳＡ１００、指標値算出処理ＳＡ１１０および音波形抽出処理ＳＡ１２０）をソフトウェアで実現した。しかし、特徴量算出処理ＳＡ１００を実行する特徴量算出手段、指標値算出処理ＳＡ１１０を実行する指標値算出手段、および音波形抽出処理ＳＡ１２０を実行する音波形抽出手段の各々を電子回路で構成し、これら電子回路を組み合わせて音波形抽出装置を構成しても良い。
【符号の説明】
【００４３】
１００…音波形抽出装置、１１０…制御部、１２０…インタフェース群、１３０…操作部、１４０…表示部、１５０…記憶部、１５０ａ…揮発性記憶部、１５０ｂ…不揮発性記憶部、１６０…バス。
【図１】
2010217408.tif 000003
【図２】
2010217408.tif 000004
【図３】
2010217408.tif 000005
【図４】
2010217408.tif 000006
【図５】
2010217408.tif 000007
【図６】
2010217408.tif 000008