ホーム > 特許ランキング > カシオ計算機株式会社
知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-01-05
(45)【発行日】2024-01-16
(54)【発明の名称】拍子決定方法、拍子決定装置、及びプログラム
(51)【国際特許分類】
G10G 3/04 20060101AFI20240109BHJP
G10L 25/51 20130101ALI20240109BHJP
【ＦＩ】
G10G3/04
G10L25/51 300
【請求項の数】
8
(21)【出願番号】P 2021155215
(22)【出願日】2021-09-24
(65)【公開番号】P2023046560
(43)【公開日】2023-04-05
【審査請求日】2022-09-26
(73)【特許権者】
【識別番号】000001443
【氏名又は名称】カシオ計算機株式会社
(74)【代理人】
【識別番号】110004185
【氏名又は名称】インフォート弁理士法人
(74)【代理人】
【識別番号】100121083
【弁理士】
【氏名又は名称】青木 宏義
(74)【代理人】
【識別番号】100138391
【弁理士】
【氏名又は名称】天田 昌行
(74)【代理人】
【識別番号】100074099
【弁理士】
【氏名又は名称】大菅 義之
(74)【代理人】
【識別番号】100182936
【弁理士】
【氏名又は名称】矢野 直樹
(72)【発明者】
【氏名】南高 純一
【審査官】大野 弘
(56)【参考文献】
【文献】特開２００９－００３２２５（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ１０Ｇ ３／０４
Ｇ１０Ｌ ２５／５１
(57)【特許請求の範囲】
【請求項1】
コンピュータが、
楽曲データの或るタイミングにおける全周波数帯域である第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる１つ以上の或る周波数帯域を含む第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、前記第１周波数帯域と前記第２周波数帯域毎とで夫々ビートレベル変動波形を取得し、
前記ビートレベル変動波形が各周波数帯域に応じた重み付けで重み付き平均されることにより重み付き平均ビートレベル変動波形を取得し、
取得された前記重み付き平均ビートレベル変動波形の自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
拍子決定方法。
【請求項2】
前記コンピュータが、
前記重み付き平均ビートレベル変動波形に基づいて、前記楽曲データの或るタイミングから設定時間長の比較元と、前記或るタイミングから離れたタイミングから前記設定時間長の比較先を取得し、
取得された前記比較元と前記比較先に基づいて、前記自己相関値を計算する、
請求項１に記載の拍子決定方法。
【請求項3】
前記第１ビートレベルは、前記第１パワーレベルの変動がマイナスからプラスに転じた第１タイミングのパワーレベルと、前記変動が次にプラスからマイナスに転じた第２タイミングのパワーレベルとのレベル差を含み、前記コンピュータが、
前記レベル差を前記第２タイミングにおける前記第１ビートレベルとして取得する、
請求項１又は２に記載の拍子決定方法。
【請求項4】
前記コンピュータが、
前記重み付き平均ビートレベル変動波形に基づいて、前記自己相関値のヒストグラムを生成し、
生成された前記ヒストグラムに基づいて、前記拍子を決定する、
請求項１乃至３の何れかに記載の拍子決定方法。
【請求項5】
前記コンピュータが、
前記ヒストグラムから抽出した、値が上位所定数個までの相関値のヒストグラムのピークの夫々が、３拍子、４拍子又は５拍子のいずれの拍子の倍数関係にあるかを判定し、
判定結果に基づいて、前記拍子を決定する、
請求項４に記載の拍子決定方法。
【請求項6】
前記第２周波数帯域は、バスドラム帯域、スネアドラム帯域又はコード帯域の何れかを含む、請求項１乃至５の何れかに記載の拍子決定方法。
【請求項7】
楽曲データの或るタイミングにおける全周波数帯域である第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる１つ以上の或る周波数帯域を含む第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、前記第１周波数帯域と前記第２周波数帯域毎とで夫々ビートレベル変動波形を取得し、
前記ビートレベル変動波形が各周波数帯域に応じた重み付けで重み付き平均されることにより重み付き平均ビートレベル変動波形を取得し、
取得された前記重み付き平均ビートレベル変動波形の自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
処理を実行する拍子決定装置。
【請求項8】
コンピュータが、楽曲データの或るタイミングにおける全周波数帯域である第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる１つ以上の或る周波数帯域を含む第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、前記第１周波数帯域と前記第２周波数帯域毎とで夫々ビートレベル変動波形を取得し、
前記ビートレベル変動波形が各周波数帯域に応じた重み付けで重み付き平均されることにより重み付き平均ビートレベル変動波形を取得し、
取得された前記重み付き平均ビートレベル変動波形の自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
処理を実行するためのプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、拍子決定方法、拍子決定装置、及びプログラムに関する。
【背景技術】
【0002】
従来、楽曲音を示す楽曲音データのテンポを解析する技術が知られている（例えば、特許文献１）。楽曲音からテンポを抽出できれば、例えばテンポを変えたオーディオデータの再生や、他のＭＩＤＩ（Ｍｕｓｉｃａｌ Ｉｎｓｔｒｕｍｅｎｔ Ｄｉｇｉｔａｌ Ｉｎｔｅｒｆａｃｅ）データと重ねて、テンポを合わせた再生などが可能となる。
【先行技術文献】
【特許文献】
【0003】
【文献】特開２００７－２７２１１８号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、楽曲音は、テンポだけでなく、拍子によっても特徴付けられているため、拍子を認識したいという要請もある。しかしながら、従来、楽曲データからの拍子の認識は、演奏者が実際に聞いて認識しており、拍子を自動認識する技術は知られていなかった。
【0005】
そこで、本発明は、楽音データから良好に拍子を決定することを目的とする。
【課題を解決するための手段】
【0006】
態様の一例の拍子決定方法では、コンピュータが、楽曲データの或るタイミングにおける全周波数帯域である第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる１つ以上の或る周波数帯域を含む第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、前記第１周波数帯域と前記第２周波数帯域毎とで夫々ビートレベル変動波形を取得し、前記ビートレベル変動波形が各周波数帯域に応じた重み付けで重み付き平均されることにより重み付き平均ビートレベル変動波形を取得し、取得された前記重み付き平均ビートレベル変動波形の自己相関値が高い複数のタイミングを取得し、取得された前記複数のタイミングに基づいて、拍子を決定する。
【発明の効果】
【0007】
本発明によれば、楽音データから良好に拍子を決定することが可能となる。
【図面の簡単な説明】
【0008】
【図1】拍子決定装置の実施形態のブロック図である。
【図2】周波数毎のパワー変動の模式図である。
【図3】パワー変動波形からビートレベル変動波形の導出の説明図である。
【図4】重み付き平均ビートレベル変動波形の自己相関算出処理の説明図である。
【図5】重み付き平均ビートレベル変動波形に対する自己相関波形例（４拍子の場合）を示す図である。
【図6】重み付き平均ビートレベル変動波形に対する自己相関波形例（３拍子の場合）を示す図である。
【図7】自己相関ヒストグラムの例を示す図である。
【図8】拍子決定のメイン処理の例を示すフローチャートである。
【図9】拍解析処理の例を示すフローチャート（その１）である。
【図10】拍解析処理の例を示すフローチャート（その２）である。
【図11】自己相関計算処理の詳細例を示すフローチャートである。
【図12】比較元の検証処理の詳細例を示すフローチャートである。
【図13】３拍子の検証処理の詳細例を示すフローチャートである。
【図14】５拍子の検証処理の詳細例を示すフローチャートである。
【図15】小節線位置特定処理の例を示すフローチャートである。
【発明を実施するための形態】
【0009】
図１は、拍子決定装置の実施形態１００のブロック図である。拍子決定装置１００は、ＣＰＵ１０１、ＲＯＭ（リードオンリーメモリ）１０２、ＲＡＭ（ランダムアクセスメモリ）１０３、入力部１０４、表示部１０５、及び出力部１０６が、システムバス１０７によって相互に接続される構成を有する。
【0010】
ＣＰＵ１０１は、拍子決定装置１００の全体の制御と、拍解析処理を実行する。
【0011】
ＲＯＭ１０２は、制御プログラムと、データベースを格納する。
【0012】
ＲＡＭ１０３は、制御プログラム実行時の変数などを記憶する。
【0013】
入力部１０４は、音楽オーディオデータを入力する部分であり、オーディオファイル形式のデータを入力する。
【0014】
表示部１０５は、処理結果を表示する。
【0015】
出力部１０６は、音楽オーディオを再生する。
【0016】
図１の拍子決定装置１００の実施形態の動作概要について、以下に説明する。図２は、楽曲データの再生に対する周波数毎のパワー変動の模式図である。３次元奥行き方向の斜め軸は周波数［Ｈｚ（ヘルツ）］を示しており、横軸は経過時刻［秒］、縦軸はパワーレベル［ｄＢ］を示している。
【0017】
一般に音楽のテンポはその音楽で演奏される楽器や歌われる歌のリズム構造によって具現化されるものである。ドラムやベース、ギター、鍵盤楽器や歌など様々な楽器と歌で音楽が構成され、それぞれのパートがテンポやリズム構造に影響を与える。また、一般にテンポやリズムをキープするのはドラムやギター、キーボードなどの楽器であることが多く、歌はそのリズム構造の上で或る程度自由に揺れがあったり自由な表現がなされる。また、リズム構造は、拍、小節、ビートなどの夫々の階層での周期性を持つことで音楽の秩序が生まれる。
【0018】
図２に例示される周波数スペクトルの時間変化では、周波数帯域別に周期性の特徴が表れていることがわかる。例えば帯域Ａに着目すると、その帯域での周波数成分のパワーは、経過時刻に沿って周期的に大きく変動している。このパワー変動２０１は、楽器演奏のリズムに対応している。ここで、帯域Ａは低い周波数帯域である。従って、これらの４個の大きなパワーのピークは、例えば低い周波数成分が多く含まれる楽音を発するバスドラムが、例えば４拍子でリズム演奏されることに起因していると考えられる。
【0019】
次に、中間周波数帯域である帯域Ｂに着目すると、やはり経過時刻に沿ったパワー変動２０２が見られる。ただし、この場合は、大きなピークの数は２個になっている。従って、これらの２個の大きなパワーのピークは、例えば中間帯域の周波数成分が多く含まれる楽音を発するスネアドラムが、例えば４拍子における強拍又は弱拍の２つの発音タイミングでリズム演奏されることに起因していると考えられる。
【0020】
更に、高い周波数帯域である帯域Ｃに着目すると、やはり経過時刻に沿ったパワー変動２０３が見られる。ただし、この場合は、大きなピークの数は８個になっている。従って、これらの８個の大きなパワーのピークは、例えば高い周波数帯域の周波数成分が多く含まれる楽音を発するギターによるコード演奏により、例えば４拍子において８分音符の発音タイミングでリズム演奏されることに起因していると考えられる。
【0021】
以上の考察により、以下に説明する実施形態（以下「本実施形態」と記載）では、各楽器や歌の特徴を捉えやすいように周波数帯域別にスペクトルパワーの立ち上がり部分のパワーをビートレベルと定義し、このビートレベルを周波数解析結果から求める。これを、リズム構造の特徴として表れやすい周波数帯域毎に求めるようにする。
【0022】
図３は、パワー変動波形からビートレベル変動波形の導出の説明図である。フレーム単位の例えば短時間フーリエ変換で算出される或る帯域のパワー変動波形３０１に着目した場合、その波形変化がマイナス変動からプラス変動に転じたフレームｆｂから、波形変化が次にプラス変動からマイナス変動に転じるフレームｆｐまでをパワー変動波形３０１の立ち上がり部分と判断する。そして、フレームｆｂのパワーレベルからフレームｆｐのパワーレベルまでのレベル差を、フレームｆｐにおけるビートレベルｂｅａｔと定義する。図３（ａ）のパワー変動波形３０１からは、＃１、＃２、及び＃３の３個の大きなピークが抽出できる。
【0023】
従って、図３（ｂ）のビートレベル変動波形３０２は、＃１、＃２、及び＃３の各ピークフレームｆｐにおいて有意なビートレベル値を有し、その他の時刻において値０を有する波形となる。
【0024】
ところで、１つの帯域に対応するビートレベル変動波形３０２のみを拍子の検出に用いた場合、その帯域に対応した楽器の演奏形態によって、正確なビートが刻まれない可能性がある。そこで、本実施形態では、任意の周波数帯域のパワーレベルから図３に示されるようにして算出したビートレベルを、まず全周波数帯域である第１周波数帯域のパワーレベルである第１パワーレベルに応じた第１ビートレベルとして集計する。更に、上記算出したビートレベルを、上記任意の周波数帯域が含まれる第２周波数帯域（図２の帯域Ａ、Ｂ、Ｃの何れか）のパワーレベルである第２パワーレベルに応じた第２ビートレベルとして集計してゆく。この結果、第１周波数帯域、及び１つ以上の第２周波数帯域毎に、図３に示されるようなビートレベル変動波形３０２が得られる。本実施形態では、これらの帯域毎のビートレベル変動波形３０２の時間毎のビートレベル値が各帯域に応じた重み付けで重み付き平均されることにより、重み付き平均ビートレベル変動波形が算出される。そして、この重み付き平均ビートレベル変動波形に基づいて、拍子が決定される。
【0025】
このように、（１）バスドラムの帯域、（２）スネアドラムの帯域、（３）コード楽器の帯域、及び（４）帯域全体について算出された各ビートレベル変動波形３０２が重ね合わされて重み付き平均される。これにより、拍や小節に基づく規則的な音の特徴をより際立たせることができ、拍子の抽出がし易くなる。楽曲に含まれるメロディ等の規則的でない音は重ね合わせによってより強調されなくなるので、結果的に拍にかかわる音がより強調されることになる。上記（１）から（４）の重ね合わせにより、より多くの種類の音楽に対応して拍子を決定することが可能となる。
【0026】
次に、本実施形態では、上述のように算出された重み付き平均ビートレベル変動波形に基づいて、次のような比較元データと各比較先データとの各自己相関が計算される。比較元データは、楽曲データの各経過時刻から設定時間分のデータである。各比較先データは、実際に設定され得る各テンポに対応する同じ楽曲データの各経過時刻だけ比較元データから離れ（ずらされ）、各経過時刻から設定時間分の各データである。そして、自己相関計算により得られる各相関値のうち、値が高い複数（例えば５個）のタイミング(ピーク位置)及び各タイミングの相関値が取得される。そして、その取得された複数のタイミング及び各タイミングの相関値に基づいて、拍子が決定される。
【0027】
図４は、重み付き平均ビートレベル変動波形の自己相関算出処理の説明図である。４０１は、上述した重み付き平均ビートレベル変動波形である。この重み付き平均ビートレベル変動波形４０１では、重み付き平均処理により、拍子に対応する規則性が際立っていることがわかる。
【0028】
本実施形態では、この重み付き平均ビートレベル変動波形４０１において、各経過時刻が楽曲の先頭である経過時刻０秒から例えば２秒ずつ順次進められてゆきながら（図４の時刻ｔｎ１、ｔｎ２等）、楽曲データの各経過時刻の時点から設定時間＝Ｔ秒分のデータ、例えば０秒時点～Ｔ秒、２秒時点～Ｔ秒、４秒時点～Ｔ秒、つまり例えば２秒毎時点からＴ秒分のデータが、比較元データ４０２として設定される。また、実際に設定される可能性のある各テンポ、例えばテンポ＝７０～１８０ｂｐｍ（ビート／秒）に対応する各時間、例えば１．３３秒～３．４３秒だけ比較元データから離され（ずらされ）、各経過時刻の時点から設定時間分のデータが比較先データ４０３として設定される。そして、比較元データ４０２と、各比較先データ４０３との間で自己相関が計算される。
【0029】
図５は、図４のようにして算出された重み付き平均ビートレベル変動波形に対する自己相関波形例（４拍子の場合）を示す図である。
図５において、ｘ軸は、自己相関における比較元に対する比較先のずらし時間であり、テンポが７０ｂｐｍから１８０ｂｐｍまでの１小節あたりの時間（１．３３秒～３．４３秒）である。
ｙ軸は、楽曲の先頭から終端までの経過時刻である。
ｚ軸は、自己相関計算の結果である相関値である。
図５において、ｘ軸は、自己相関における比較元に対する比較先のずらし時間であり、テンポが７０ｂｐｍから１８０ｂｐｍまでの１小節あたりの時間（１．３３秒～３．４３秒）である。
ｙ軸は、楽曲の先頭から終端までの経過時刻である。
ｚ軸は、自己相関計算の結果である相関値である。
【0030】
図５に示される自己相関波形としては、ｙ軸方向の各経過時刻毎に、その経過時刻からＴ秒分の比較元データ４０２（図４参照）について、１．３３秒～３．４３秒だけ比較元データから離され（ずらされ）、各経過時刻からＴ秒分の各比較先データ４０３（図４参照）に対応するｚ軸方向の相関値が、ｘ軸方向に２次元波形としてプロットされる。更に、ｙ軸方向の各経過時刻毎に比較元データ４０２がシフトされることにより、全体として３次元の自己相関波形がプロットされている。
ｙ軸方向に一番手前の、ｘ軸方向及びｚ軸方向からなる２次元の自己相関波形は、楽曲の先頭から例えば４０秒時点における自己相関波形を示している。
ｙ軸方向に一番手前の、ｘ軸方向及びｚ軸方向からなる２次元の自己相関波形は、楽曲の先頭から例えば４０秒時点における自己相関波形を示している。
【0031】
図５に例示される３次元波形全体として、まず、比較先データ４０３が比較元データ４０２から夫々、１．４４秒、１．９２秒、２．４０秒、２．８８秒、及び３．３６秒ずれた各時点付近に、＃１～＃５の相関値のピーク５０１が現れていることがわかる。前述したように、自己相関の計算における比較元に対する比較先のずらし時間としては、楽曲として想定され得るテンポが７０ｂｐｍから１８０ｂｐｍまでの１小節あたりの時間（１．３３秒～３．４３秒）が設定されている。従って、図５に例示される自己相関波形において現れる相関値のピーク５０１は、楽曲演奏により発生するテンポ及び拍に同期したビート成分が現れていると考えられる。
【0032】
従って、まず、４拍子の楽曲を仮定した場合には、小節内の４拍に相当する拍間隔及び小節間隔で相関値のピーク５０１が並ぶはずである。
また、楽曲のテンポは不明であるため、図５に例として現れる＃１から＃５の相関値のピーク５０１のうち、どの相関値のピーク５０１に対応するずらし時間が小節長に対応するかはわからない。
また、楽曲のテンポは不明であるため、図５に例として現れる＃１から＃５の相関値のピーク５０１のうち、どの相関値のピーク５０１に対応するずらし時間が小節長に対応するかはわからない。
【0033】
しかしながら、楽曲が４拍子であると仮定され、ある相関値のピーク５０１に対応するずらし時間がもし小節長に対応していると仮定するならば、他の相関値のピーク５０１に対応するずらし時間は、その小節長に対して４拍子の各拍子タイミングに対応する倍数関係になるはずである。具体的には、４拍子の楽曲を仮定した場合、他の相関値のピーク５０１に対応するずらし時間長は、小節長に対して、３／４倍、４／４倍、５／４倍、６／４倍、７／４倍等の倍数関係になるはずである。
【0034】
本実施形態では、上記の４拍子に関する関係性が成立しない場合、次に、楽曲が３拍子であると仮定される。図６は、図４のようにして算出された重み付き平均ビートレベル変動波形に対する自己相関波形例（３拍子の場合）を示す図である。図６において、ｘ軸、ｙ軸、ｚ軸の意味は、夫々図５の場合と同じである。
【0035】
図６に示される自己相関波形の３次元形状の意味も、図５の場合と同様である。図６に例示される３次元波形全体として、まず、比較先データ４０３が比較元データ４０２から夫々、１．３７秒、１．８３秒、２．２９秒、及び２．７４秒ずれた各時点付近に、＃１～＃４の相関値のピーク６０１が現れていることがわかる。
【0036】
ある相関値のピーク６０１に対応するずらし時間がもし小節長に対応していると仮定するならば、他の相関値のピーク６０１に対応するずらし時間は、その小節長に対して３拍子の各拍子タイミングに対応する倍数関係になるはずである。具体的には、３拍子の楽曲を仮定した場合、他の相関値のピーク６０１に対応するずらし時間長は、小節長に対して、３／３倍、４／３倍、５／３倍、６／３倍等の倍数関係になるはずである。
【0037】
本実施形態では、上記の３拍子に関する関係性も成立しない場合、特には図示しないが、楽曲が５拍子であると仮定される。ある相関値のピークに対応するずらし時間がもし小節長に対応していると仮定するならば、他の相関値のピークに対応するずらし時間は、その小節長に対して５拍子の各拍子タイミングに対応する倍数関係になるはずである。具体的には、５拍子の楽曲を仮定した場合、他の相関値のピークに対応するずらし時間長は、小節長に対して、３／５倍、４／５倍、５／５倍、６／５倍、７／５倍等の倍数関係になるはずである。
【0038】
楽曲が他の拍子である場合も考えられるが、通常は３、４、５拍子を仮定すれば十分である。
【0039】
本実施形態では、上述のアルゴリズムを実現するために、次のような処理が実行される。まず、例えば図４のようにして算出された図５又は図６のような３次元の相関値のピーク５０１又は６０１について、ｙ軸方向の或る経過時刻における自己相関波形中の、例えば相関値が高い上位５個のピークの相関値が抽出される。そして、それらの相関値が、そのピーク位置（ずらし時間）に対応するヒストグラムのビン位置に累算される。この動作が、ｙ軸方向の各経過時刻毎の自己相関波形に対して実行され、同じヒストグラムに累算されてゆく。
【0040】
この結果、４拍子の楽曲に対する相関値のヒストグラムとして、例えば図７（ａ）に例示される相関値のヒストグラム７００（ａ）が得られる。この相関値のヒストグラム７００（ａ）においては、図５の＃１～＃５の相関値のピーク５０１に対応して、それらのピーク位置（ずらし時間）に対応するビン位置に、＃１～＃５の相関値のヒストグラムのピーク７０１が得られる。
同様に、３拍子の楽曲に対する相関値のヒストグラムとして、例えば図７（ｂ）に示される相関値のヒストグラム７００（ｂ）が得られる。この相関値のヒストグラム７００（ｂ）においては、図６の＃１～＃４の相関値のピーク６０１に対応して、それらのピーク位置（ずらし時間）に対応するビン位置に、＃１～＃４の相関値のヒストグラムのピーク７０２が得られる。
同様に、３拍子の楽曲に対する相関値のヒストグラムとして、例えば図７（ｂ）に示される相関値のヒストグラム７００（ｂ）が得られる。この相関値のヒストグラム７００（ｂ）においては、図６の＃１～＃４の相関値のピーク６０１に対応して、それらのピーク位置（ずらし時間）に対応するビン位置に、＃１～＃４の相関値のヒストグラムのピーク７０２が得られる。
【0041】
続いて、本実施形態では、後述する図８、図９、及び図１２の拍解析処理及び比較元の検証処理が実行される。これらの処理において、図７（ａ）のように得られた相関値のヒストグラム７００（ａ）から抽出した、値が例えば上位７個までの相関値のヒストグラムのピーク７０１の夫々について、そのピーク位置（ずらし時間）に対応するビン位置が小節時間長であると仮定される。そして、他の相関値のヒストグラムのピーク７０１に対応するビン位置（＝ずらし時間）の２つ以上が前述した４拍子の倍数関係の何れかにあるか否かを判定する処理が実行される。この判定結果が肯定的であれば、上記仮定した小節時間長は正しいと判定して、その４拍子の小節時間長を確定させ、同時にテンポも確定させる。
【0042】
相関値のヒストグラム７００が、図７（ａ）の７００（ａ）ではなく、図７（ｂ）の７００（ｂ）であったような場合であって、抽出された相関値のヒストグラムのピーク７０２の何れにおいても、４拍子に関する上記判定が肯定的にならなかった場合には、後述する図１３の３拍子の検証処理が実行される。この３拍子の検証処理において、図７（ｂ）のように得られた相関値のヒストグラム７００（ｂ）から抽出した、値が例えば上位７個までの抽出された相関値のヒストグラムのピーク７０２の夫々について、そのピーク位置（ずらし時間）に対応するビン位置が小節時間長であると仮定される。そして、他の相関値のヒストグラムのピーク７０２に対応するビン位置（＝ずらし時間）の２つ以上が、前述した３拍子の倍数関係の何れかにあるか否かを判定する処理が実行される。この判定結果が肯定的であれば、上記仮定した小節長は正しいと判定して、その３拍子の小節時間長を確定させ、同時にテンポも確定させる。
【0043】
相関値のヒストグラム７００が、図７（ａ）の７００（ａ）でも図７（ｂ）の７００（ｂ）でもなく、抽出された相関値のヒストグラムのピークの何れにおいても、４拍子に関する上記判定及び３拍子に関する上記判定が肯定的にならなかった場合には、後述する図１４の５拍子の検証処理が実行される。この５拍子の検証処理において、特には図示しないが、抽出された相関値のヒストグラムのピークの夫々について、そのピーク位置（ずらし時間）に対応するビン位置が小節時間長であると仮定される。そして、他の相関値のヒストグラムのピークに対応するビン位置（＝ずらし時間）の２つ以上が前述した５拍子の倍数関係の何れかにあるか否かを判定する処理が実行される。この判定結果が肯定的であれば、上記仮定した小節長は正しいと判定して、その５拍子の小節時間長を確定させ、同時にテンポも確定させる。
【0044】
上述のようにして本実施形態では、楽音データから良好に拍子を決定することが可能となる。
【0045】
図８は、上述した実施形態の動作を実現する拍子決定のメイン処理の例を示すフローチャートである。このメイン処理は、図１のＣＰＵ１０１がＲＯＭ１０２に記憶されているメイン処理プログラムをＲＡＭ１０３に読み出して実行する処理である。このメイン処理では、各楽器や歌の特徴を捉えやすいように周波数解析（短時間フーリエ変換）により図２に例示されるように周波数帯域別に算出したスペクトル成分のパワーの立上り部分を、図３で説明したビートレベルとして算出する。メイン処理では、この処理を、リズム構造の特徴として表れやすい周波数帯域毎に算出する。更に、メイン処理では、これらの処理を、楽曲データの全体に対して、時間（フレーム）を少しずつずらしながら、実行していく。
【0046】
ＣＰＵ１０１はまず、初期化処理として、楽曲の経過時刻を一定時間（例えば２５６～１０２４ミリ秒）毎のフレームとして指定するためのＲＡＭ１０３に記憶されている変数であるフレームカウンタ変数ｆｒｍの値を０（ゼロ）にリセットする（ステップＳ８０１）。
【0047】
次に、ＣＰＵ１０１は、ステップＳ８０２で楽曲データの全てのフレームに対して処理が終了したと判定するまで、ステップＳ８２０でフレームカウンタ変数ｆｒｍの値を＋１ずつインクリメントしながら、ステップＳ８０２からステップＳ８２０までの一連の処理を繰り返し実行する。
【0048】
この繰返し処理において、ＣＰＵ１０１はまず、図１の入力部１０４からＲＡＭ１０３に読み込んだフレームカウンタ変数ｆｒｍが示す現在のフレームの楽曲データに対して、周波数解析処理である短時間フーリエ変換の演算を実行する（ステップＳ８０３）。
【0049】
次に、ＣＰＵ１０１は、ステップＳ８０３の短時間フーリエ変換の演算により算出された各周波数成分から、各周波数成分のパワーを計算する（ステップＳ８０４）。周波数成分毎のパワー値は、周波数成分の周波数位置であるｂｉｎ値をキーとして、ＲＡＭ１０３上の配列変数であるパワー配列変数ｄｏＤａｔａ［ｂｉｎ］に得られる。
【0050】
続いて、ＣＰＵ１０１は、ステップＳ８０５で、前述したｂｉｎ値を指定するためのＲＡＭ１０３に記憶される変数であるｂｉｎ変数の値を１にリセットした後、ステップＳ８１８でそのｂｉｎ変数の値を＋１ずつインクリメントしながら、ステップＳ８０６で全てのｂｉｎ値について処理が完了したと判定するまで、ステップＳ８０６からステップＳ８１８までの一連の処理をｂｉｎ値毎に繰返し実行する。
【0051】
このｂｉｎ値の繰返し処理において、ＣＰＵ１０１はまず、ｂｉｎ変数が示す現在のｂｉｎ（周波数成分）値に対応してＲＡＭ１０３に記憶されている現在フレームパワー配列変数ｄｏＤａｔａ［ｂｉｎ］の値から、現在のフレームの１つ前のフレームにおける現在のｂｉｎ値に対応してＲＡＭ１０３に記憶されている前フレームパワー配列変数ｄｏＤａｔａＢｕｆ［ｂｉｎ］の値を減算する。そして、ＣＰＵ１０１は、その減算結果である差分値をＲＡＭ１０３に記憶される変数である差分値変数ｄｉｖ１に格納する（ステップＳ８０７）。この差分値は、前フレームと現在フレームとの間のパワーの変化分を示している。
【0052】
次に、ＣＰＵ１０１は、ステップＳ８０７で算出した差分値変数ｄｉｖ１の値が、０（ゼロ）よりも大きいか否かを判定する（ステップＳ８０８）。この判定処理では、現在フレームにおける現在のｂｉｎ（周波数成分）値に対応するパワーが、プラス変動になって（増加して）いるか、マイナス変動になって（減少して）いるか（変動無しを含む）が判定される。
【0053】
現在フレームにおける現在のｂｉｎ（周波数成分）値のパワー変動がプラス変動になっておりステップＳ８０８の判定がＹＥＳの場合には、ＣＰＵ１０１は、次の処理を実行する。ＣＰＵ１０１は、現在のｂｉｎ（周波数成分）値のパワーレベル値を示すＲＡＭ１０３に記憶される配列変数であるレベル配列変数Ｌｖ［ｂｉｎ］に、ステップＳ８０７で算出された差分値変数ｄｉｖ１の値をレベル増加分として累算する（ステップＳ８１７）。なお、図８のステップＳ８１７において演算記号「＋＝」は、その左辺の変数にその右辺の値を累算する演算を示している。
【0054】
ステップＳ８１７の処理の後、ＣＰＵ１０１は、ステップＳ８１８でｂｉｎ変数の値を＋１インクリメントし、その後にステップＳ８０６からステップＳ８０７の処理に移行して、次のｂｉｎ（周波数成分）値に対する処理を繰り返し実行する。
【0055】
やがて、現在フレームにおける現在のｂｉｎ（周波数成分）値のパワー変動がマイナス変動（又は変動無し）に転じて差分値変数ｄｉｖ１の値が０以下となり、ステップＳ８０８の判定がＮＯとなる。これは、前述した図３においてフレームｆｐが到来したことを意味する。
【0056】
この時点において、レベル配列変数Ｌｖ［ｂｉｎ］には、現在のｂｉｎ（周波数成分）値に対応する図３に示されるビートレベルｂｅａｔが得られている。ＣＰＵ１０１はまず、このレベル配列変数Ｌｖ［ｂｉｎ］の値を、全周波数帯域である第１周波数帯域の、フレームカウンタ変数ｆｒｍの値が示す現在フレームにおけるビートレベルが格納されるＲＡＭ１０３に記憶される配列変数である第１のビートレベル変動波形配列変数ＢＬ１［ｆｒｍ］に累算する（ステップＳ８０９）。第１周波数帯域は、例えば、サンプリング周波数が４４．１ｋＨｚ（キロヘルツ）の場合、周波数０～２２０５０Ｈｚの範囲を有する。
【0057】
次に、ＣＰＵ１０１は、ｂｉｎ変数の値が示す現在のｂｉｎ（周波数成分）値が、ＢＤ（バスドラム）帯域に属するか否かを判定する（ステップＳ８１０）。ＢＤ帯域は、例えば、周波数２０～１００Ｈｚの範囲を有し、例えば図２の帯域Ａに対応する。
【0058】
ステップＳ８１０の判定がＹＥＳならば、ＣＰＵ１０１は、レベル配列変数Ｌｖ［ｂｉｎ］の値を、第２周波数帯域の一帯域であるＢＤ帯域の、フレームカウンタ変数ｆｒｍの値が示す現在フレームにおけるビートレベルが格納されるＲＡＭ１０３に記憶される配列変数である第２のビートレベル変動波形配列変数ＢＬ２［ｆｒｍ］に累算する（ステップＳ８１１）。その後、ＣＰＵ１０１は、ステップＳ８１６の処理に移行する。
【0059】
ステップＳ８１０の判定がＮＯならば、ＣＰＵ１０１は、ｂｉｎ変数の値が示す現在のｂｉｎ（周波数成分）値が、ＳＤ（スネアドラム）帯域に属するか否かを判定する（ステップＳ８１２）。ＳＤ帯域は、例えば、周波数１２５～２５０Ｈｚの範囲を有し、例えば図２の帯域Ｂに対応する。
【0060】
ステップＳ８１２の判定がＹＥＳならば、ＣＰＵ１０１は、レベル配列変数Ｌｖ［ｂｉｎ］の値を、第２周波数帯域の一帯域であるＳＤ帯域の、フレームカウンタ変数ｆｒｍの値が示す現在フレームにおけるビートレベルが格納されるＲＡＭ１０３に記憶される配列変数である第３のビートレベル変動波形配列変数ＢＬ３［ｆｒｍ］に累算する（ステップＳ８１３）。その後、ＣＰＵ１０１は、ステップＳ８１６の処理に移行する。
【0061】
ステップＳ８１２の判定がＮＯならば、ＣＰＵ１０１は、ｂｉｎ変数の値が示す現在のｂｉｎ（周波数成分）値が、コード帯域に属するか否かを判定する（ステップＳ８１４）。コード帯域は、例えば、周波数３００～６００Ｈｚの範囲を有し、例えば図２の帯域Ｃに対応する。
【0062】
ステップＳ８１４の判定がＹＥＳならば、ＣＰＵ１０１は、レベル配列変数Ｌｖ［ｂｉｎ］の値を、第２周波数帯域の一帯域であるコード帯域の、フレームカウンタ変数ｆｒｍの値が示す現在フレームにおけるビートレベルが格納されるＲＡＭ１０３に記憶される配列変数である第４のビートレベル変動波形配列変数ＢＬ４［ｆｒｍ］に累算する（ステップＳ８１５）。その後、ＣＰＵ１０１は、ステップＳ８１６の処理に移行する。
【0063】
ステップＳ８１１、ステップＳ８１３、又はステップＳ８１５の処理の後、或いはステップＳ８１４の判定がＮＯの場合に、ＣＰＵ１０１は、現在のｂｉｎ（周波数成分）値に対応するレベル配列変数Ｌｖ［ｂｉｎ］の値を０（ゼロ）にクリアする。その後、ＣＰＵ１０１は、ｂｉｎ（周波数成分）値の値を＋１インクリメントし（ステップＳ８１８），ステップＳ８０６の処理からステップＳ８０７の処理に移行して、次のｂｉｎ（周波数成分）値に対応する処理を繰り返し実行する。
【0064】
以上のステップＳ８０６からステップＳ８１８の処理の繰返しにより、全てのｂｉｎ（周波数成分）値に対する処理が終了すると、ステップＳ８０６の判定がＹＥＳとなる。この結果、ＣＰＵ１０１は、ＲＡＭ１０３に記憶されている現在フレームに対応する全てのｂｉｎ（周波数成分）値に対応する現在フレームパワー配列変数ｄｏＤａｔａ［ｎ］を、前フレームパワー配列変数ｄｏＤａｔａＢｕｆ［ｎ］（ｎはすべてのｂｉｎ値を表す）にコピーし、ＲＡＭ１０３に記憶させる（ステップ８１９）。
【0065】
その後、ＣＰＵ１０１は、フレームカウンタ変数ｆｒｍの値を＋１インクリメントし（ステップＳ８２０）、そして、ステップＳ８０２からステップＳ８０３に移行して、次のフレームに対する処理を繰り返し実行する。
【0066】
やがて、楽曲データの末尾までの全てのフレームに対する処理が完了すると、ステップＳ８０２の判定がＹＥＳとなって、図８のメイン処理を終了する。
【0067】
以上のメイン処理により、全帯域に対応する第１のビートレベル変動波形配列変数ＢＬ１［ｆｒｍ］に得られる第１のビートレベル変動波形ＢＬ１［ｆｒｍ］、ＢＤ帯域に対応する第２のビートレベル変動波形配列変数ＢＬ２［ｆｒｍ］に得られる第２のビートレベル変動波形ＢＬ２［ｆｒｍ］、ＳＤ帯域に対応する第３のビートレベル変動波形配列変数ＢＬ３［ｆｒｍ］に得られる第３のビートレベル変動波形ＢＬ３［ｆｒｍ］、及びコード帯域に対応する第４のビートレベル変動波形配列変数ＢＬ４［ｆｒｍ］に得られる第４のビートレベル変動波形ＢＬ４［ｆｒｍ］として、図３（ｂ）の３０２に例示したのと同様の形状の、一楽曲データ分のビートレベル変動波形が得られる。
【0068】
図９及び図１０は、拍解析処理の例を示すフローチャートである。拍解析処理は、メイン処理の場合と同様に、図１のＣＰＵ１０１がＲＯＭ１０２に記憶されている拍解析処理プログラムをＲＡＭ１０３に読み出して実行する処理である。
この拍解析処理では、まず図３で前述したようにして、図８のメイン処理で得られた全帯域に対応する第１のビートレベル変動波形ＢＬ１［ｆｒｍ］、ＢＤ帯域に対応する第２のビートレベル変動波形ＢＬ２［ｆｒｍ］、ＳＤ帯域に対応する第３のビートレベル変動波形ＢＬ３［ｆｒｍ］、及びコード帯域に対応する第４のビートレベル変動波形ＢＬ４［ｆｒｍ］から、重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ（Ｗｅｉｇｈｔｅｄ Ａｒｉｔｈｍｅｔｉｃ Ｍｅａｎ Ｂｅａｔ Ｌｅｖｅｌ）が算出される。
次に、拍解析処理では、図４～図６で前述したようにして、重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬに対して、自己相関が計算される。
更に、拍解析処理では、図７で前述したようにして、自己相関の計算により得られる相関値から相関値のヒストグラム（図７の７００に対応）が算出される。
そして、拍解析処理では、前述したようにして、相関値のヒストグラムのピークに基づいて、楽曲データの拍子を推定する処理が実行される。
この拍解析処理では、まず図３で前述したようにして、図８のメイン処理で得られた全帯域に対応する第１のビートレベル変動波形ＢＬ１［ｆｒｍ］、ＢＤ帯域に対応する第２のビートレベル変動波形ＢＬ２［ｆｒｍ］、ＳＤ帯域に対応する第３のビートレベル変動波形ＢＬ３［ｆｒｍ］、及びコード帯域に対応する第４のビートレベル変動波形ＢＬ４［ｆｒｍ］から、重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ（Ｗｅｉｇｈｔｅｄ Ａｒｉｔｈｍｅｔｉｃ Ｍｅａｎ Ｂｅａｔ Ｌｅｖｅｌ）が算出される。
次に、拍解析処理では、図４～図６で前述したようにして、重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬに対して、自己相関が計算される。
更に、拍解析処理では、図７で前述したようにして、自己相関の計算により得られる相関値から相関値のヒストグラム（図７の７００に対応）が算出される。
そして、拍解析処理では、前述したようにして、相関値のヒストグラムのピークに基づいて、楽曲データの拍子を推定する処理が実行される。
【0069】
上述の重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬの算出処理と、重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬに対する自己相関の計算処理は、図９のステップＳ９０１からステップＳ９０８で実行され、上述の相関値のヒストグラムの算出処理は図９のステップＳ９０９とステップＳ９１０で実行される。更に、上述の相関値のヒストグラムのピークに基づいて楽曲データの拍子を推定する処理は、図１０のステップＳ９１２からステップＳ９２３で実行される。
【0070】
ここで、ＲＡＭ１０３に記憶される比較元先頭位置変数ｄｏＯｒｇとして、図４の比較元データ４０２の先頭位置が、秒を単位とするデータとして格納される。同様に、ＲＡＭ１０３に記憶される比較先先頭位置変数ｄｏＤｓｔとして、図４の比較先データ４０３の先頭位置が、秒を単位とするデータとして格納される。以下の説明では、それぞれの変数に格納される値をその変数名と同じ記号で示す。
比較元先頭位置ｄｏＯｒｇ［秒］は、ＲＯＭ１０２（可変を許容する場合にはＲＡＭ１０３）に記録される比較元時間ステップ幅変数ｄｏＯｒｇＳｔｅｐが示す比較元時間ステップ幅ｄｏＯｒｇＳｔｅｐ［秒］ずつ増加させられながら、楽曲データの先頭から終端までの間で指定される。比較元時間ステップ幅ｄｏＯｒｇＳｔｅｐの値は、例えば２秒である。
また、比較先先頭位置ｄｏＤｓｔは、楽曲データとして指定され得るテンポ範囲が例えば６０から１８０の範囲であるように、値がセットされる。４拍子の楽曲データでは、テンポ６０の場合に小節の長さが４秒、テンポ１８０の場合に小節の長さが１．３３秒となる。つまり、比較先先頭位置ｄｏＤｓｔとしては、１．３３［秒］から４．００［秒］の間の値が、必要な分解能でずらしながら指定される。今、このずらし幅を、ＲＯＭ１０２（可変を許容する場合にはＲＡＭ１０３）に記録される比較先時間ステップ幅変数ｄｏＤｓｔＳｔｅｐが示す比較先時間ステップ幅ｄｏＤｓｔＳｔｅｐ［秒］とする。このとき、比較先先頭位置ｄｏＤｓｔは、下記（１）式で示される演算処理により計算される。
ｄｏＤｓｔ＝ｄｏＯｒｇ＋１．３３＋ｋ×ｄｏＤｓｔＳｔｅｐ ・・・（１）
ここで、ｋは、０以上の整数である。
比較元先頭位置ｄｏＯｒｇ［秒］は、ＲＯＭ１０２（可変を許容する場合にはＲＡＭ１０３）に記録される比較元時間ステップ幅変数ｄｏＯｒｇＳｔｅｐが示す比較元時間ステップ幅ｄｏＯｒｇＳｔｅｐ［秒］ずつ増加させられながら、楽曲データの先頭から終端までの間で指定される。比較元時間ステップ幅ｄｏＯｒｇＳｔｅｐの値は、例えば２秒である。
また、比較先先頭位置ｄｏＤｓｔは、楽曲データとして指定され得るテンポ範囲が例えば６０から１８０の範囲であるように、値がセットされる。４拍子の楽曲データでは、テンポ６０の場合に小節の長さが４秒、テンポ１８０の場合に小節の長さが１．３３秒となる。つまり、比較先先頭位置ｄｏＤｓｔとしては、１．３３［秒］から４．００［秒］の間の値が、必要な分解能でずらしながら指定される。今、このずらし幅を、ＲＯＭ１０２（可変を許容する場合にはＲＡＭ１０３）に記録される比較先時間ステップ幅変数ｄｏＤｓｔＳｔｅｐが示す比較先時間ステップ幅ｄｏＤｓｔＳｔｅｐ［秒］とする。このとき、比較先先頭位置ｄｏＤｓｔは、下記（１）式で示される演算処理により計算される。
ｄｏＤｓｔ＝ｄｏＯｒｇ＋１．３３＋ｋ×ｄｏＤｓｔＳｔｅｐ ・・・（１）
ここで、ｋは、０以上の整数である。
【0071】
図９のフローチャートにおいて、ＣＰＵ１０１はまず、初期化処理を実行し、比較元データの進行を制御するＲＡＭ１０３に記憶される比較元カウンタ変数ｎに０（ゼロ）を設定し、比較元先頭位置変数ｄｏＯｒｇに楽曲データの先頭を表す０．０［秒］を設定する（ステップＳ９０１）。
【0072】
次に、ＣＰＵ１０１は、ステップＳ９１１で比較元カウンタ変数ｎの値を＋１ずつインクリメントすると共に、比較元先頭位置変数ｄｏＯｒｇにＲＯＭ１０２（可変を許容する場合はＲＡＭ１０３）から読み出した比較元時間ステップ幅変数ｄｏＯｒｇＳｔｅｐの値を順次累算しながら、ステップＳ９０２で比較元データの指定が楽曲データの末尾まで達したと判定するまで、ステップＳ９０２からステップＳ９１１までの一連の比較元毎処理を繰り返し実行する。ステップＳ９１１での累算処理の繰返しによって、「ｄｏＯｒｇ＝ｄｏＯｒｇＳｔｅｐ×ｎ」という比較元先頭位置変数ｄｏＯｒｇを進行させる処理が実行されることになり、例えば図４において、楽曲の先頭から末尾までの各経過時刻ｔｎ１、ｔｎ２等が順次指定されてゆくことになる。
【0073】
上述の比較元毎処理の繰返しにおいて、ＣＰＵ１０１はまず、比較先の位置を指定するためのＲＡＭ１０３に記憶される比較先カウンタ変数ｋ（前述の（１）式を参照）に初期値０（ゼロ）を設定する（ステップＳ９０３）。
【0074】
次に、ＣＰＵ１０１は、現在指定されている比較元先頭位置変数ｄｏＯｒｇ（ステップＳ９０１とステップＳ９１１を参照）と、ステップＳ９０３で初期設定された比較先カウンタ変数ｋの値＝０を用いて、前述した（１）式により、比較先先頭位置ｄｏＤｓｔの初期値を算出する（ステップＳ９０４）。
【0075】
そして、ＣＰＵ１０１は、ステップＳ９０８で、比較先カウンタ変数ｋの値を＋１ずつインクリメントすると共に、ステップＳ９０４で初期設定された比較先先頭位置変数ｄｏＤｓｔにＲＯＭ１０２（可変を許容する場合はＲＡＭ１０３）から読み出した比較先時間ステップ幅変数ｄｏＤｓｔＳｔｅｐの値を順次累算しながら、ステップＳ９０５で比較先データの指定が終了したと判定するまで、ステップＳ９０５からステップＳ９０８までの一連の比較先毎処理を繰り返し実行する。
【0076】
上記比較先毎処理の繰返しにおいて、ＣＰＵ１０１はまず、自己相関計算処理を実行する（ステップＳ９０６）。
【0077】
図１１は、図９のステップＳ９０６の自己相関計算処理の詳細例を示すフローチャートである。図１１において、ＣＰＵ１０１は、ステップＳ１１０１で比較元及び比較先の設定時間内の進行を制御するＲＡＭ１０３に記憶される設定時間内カウンタ変数ｉに初期値０（ゼロ）を格納する。その後、ＣＰＵ１０１は、ステップＳ１１０７で設定時間内カウンタ変数ｉの値を＋１ずつインクリメントしながら、ステップＳ１１０２で設定時間内カウンタ変数ｉの値が設定時間Ｔ（図４参照）に対応するＲＯＭ１０２（可変を許容する場合はＲＡＭ１０３）に記憶されている設定時間サンプル数Ｎｕｍに到達したと判定するまで、ステップＳ１１０２からステップＳ１１０７の一連の設定時間内毎処理を繰り返し実行する。
【0078】
上述の設定時間内毎処理の繰返しにおいて、まずＣＰＵ１０１は、下記（２）式及び（３）式で示される各演算処理に基づいて、比較元データの設定時間内及び比較先データの設定時間内の各ｉ番目のデータ位置ｐ０及びｐ１［秒］を計算する（ステップＳ１１０３）。
ｐ０＝ｄｏＯｒｇ＋ｉ×４×（ｄｏＤｓｔ－ｄｏＯｒｇ）／Ｎｕｍ ・・・（２）
ｐ１＝ｄｏＤｓｔ＋ｉ×４×（ｄｏＤｓｔ－ｄｏＯｒｇ）／Ｎｕｍ ・・・（３）
上記（２）式及び（３）式において、設定時間内カウンタ変数ｉの値が図４の設定時間Ｔに相当する設定時間サンプル数Ｎｕｍの値になったときに、上記（２）式及び（３）式で計算される各データ位置は、共に「４×（ｄｏＤｓｔ－ｄｏＯｒｇ）」となる。即ち、設定時間Ｔは、比較先データの比較元データに対するずらし時間の４倍の時間に設定される。
上述の（２）式及び（３）式によって示される演算処理によれば、設定時間Ｔは、固定時間ではなく、比較先データと比較元データのずらし時間の４倍の値｛４×（ｄｏＤｓｔ－ｄｏＯｒｇ）｝になって、ずらし時間に応じた時間範囲になる。このようにして、自己相関演算のためにずらし時間に応じた適切な設定時間Ｔが設定されることになる。
ｐ０＝ｄｏＯｒｇ＋ｉ×４×（ｄｏＤｓｔ－ｄｏＯｒｇ）／Ｎｕｍ ・・・（２）
ｐ１＝ｄｏＤｓｔ＋ｉ×４×（ｄｏＤｓｔ－ｄｏＯｒｇ）／Ｎｕｍ ・・・（３）
上記（２）式及び（３）式において、設定時間内カウンタ変数ｉの値が図４の設定時間Ｔに相当する設定時間サンプル数Ｎｕｍの値になったときに、上記（２）式及び（３）式で計算される各データ位置は、共に「４×（ｄｏＤｓｔ－ｄｏＯｒｇ）」となる。即ち、設定時間Ｔは、比較先データの比較元データに対するずらし時間の４倍の時間に設定される。
上述の（２）式及び（３）式によって示される演算処理によれば、設定時間Ｔは、固定時間ではなく、比較先データと比較元データのずらし時間の４倍の値｛４×（ｄｏＤｓｔ－ｄｏＯｒｇ）｝になって、ずらし時間に応じた時間範囲になる。このようにして、自己相関演算のためにずらし時間に応じた適切な設定時間Ｔが設定されることになる。
【0079】
次に、ＣＰＵ１０１は、上記（２）式及び（３）式の演算で計算された比較元データの設定時間内及び比較先データの設定時間内の各ｉ番目のデータ位置ｐ０及びｐ１［秒］に基づいて、下記（４）式及び（５）式で示される各演算処理を実行する。この結果、ＣＰＵ１０１は、比較元データの設定時間内及び比較先データの設定時間内の各ｉ番目のサンプルデータへのインデックスである、比較元サンプルｉインデックスｉｄｘＯｒｇ＿ｉ及び比較先サンプルｉインデックスｉｄｘＤｓｔ＿ｉを算出する（ステップＳ１１０４）。
ｉｄｘＯｒｇ＿ｉ＝ｐ０×サンプリング周波数［Ｈｚ］ ・・・（４）
ｉｄｘＤｓｔ＿ｉ＝ｐ１×サンプリング周波数［Ｈｚ］ ・・・（５）
ｉｄｘＯｒｇ＿ｉ＝ｐ０×サンプリング周波数［Ｈｚ］ ・・・（４）
ｉｄｘＤｓｔ＿ｉ＝ｐ１×サンプリング周波数［Ｈｚ］ ・・・（５）
【0080】
続いて、ＣＰＵ１０１は、（４）式及び（５）式の各演算で算出された比較元サンプルｉインデックスｉｄｘＯｒｇ＿ｉ及び比較先サンプルｉインデックスｉｄｘＤｓｔ＿ｉが含まれる各フレーム番号である、比較元フレームインデックスｉｄｘＯｒｇ＿ｆ及び比較先フレームインデックスｉｄｘＤｓｔ＿ｆを、下記（６）式及び（７）式で示される演算処理により算出する（ステップＳ１００５）。ただし、ｆｓｉｚｅはフレームサイズ（単位は「サンプル」）である。
ｉｄｘＯｒｇ＿ｆ＝ｉｄｘＯｒｇ＿ｉ／ｆｓｉｚｅ ・・・（６）
ｉｄｘＤｓｔ＿ｆ＝ｉｄｘＤｓｔ＿ｉ／ｆｓｉｚｅ ・・・（７）
ｉｄｘＯｒｇ＿ｆ＝ｉｄｘＯｒｇ＿ｉ／ｆｓｉｚｅ ・・・（６）
ｉｄｘＤｓｔ＿ｆ＝ｉｄｘＤｓｔ＿ｉ／ｆｓｉｚｅ ・・・（７）
【0081】
そして、ＣＰＵ１０１は、上述の（６）式の演算で算出された比較元フレームインデックスｉｄｘＯｒｇ＿ｆをキーとして、下記（８）式で示される演算処理を実行する。
ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］＝Ａ×第１のビートレベル変動波形ＢＬ１［ｉｄｘＯｒｇ＿ｆ］＋Ｂ×第２のビートレベル変動波形ＢＬ２［ｉｄｘＯｒｇ＿ｆ］＋Ｃ×第３のビートレベル変動波形ＢＬ３［ｉｄｘＯｒｇ＿ｆ］＋Ｄ×第４のビートレベル変動波形ＢＬ４［ｉｄｘＯｒｇ＿ｆ］
・・・（８）
この演算処理により、ＣＰＵ１０１は、重み係数Ａ、Ｂ、Ｃ、及びＤを用いて、第１のビートレベル変動波形ＢＬ１［ｉｄｘＯｒｇ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＯｒｇ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＯｒｇ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＯｒｇ＿ｆ］に対する比較元重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］を計算する（ステップＳ１１０６）。ここで、重み係数Ａ、Ｂ、Ｃ、及びＤは、例えばＲＯＭ１０２に記憶されており、又は可変を許容する場合はＲＡＭ１０３に記憶されている。また、第１のビートレベル変動波形ＢＬ１［ｉｄｘＯｒｇ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＯｒｇ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＯｒｇ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＯｒｇ＿ｆ］は夫々、図８のフローチャートのステップＳ８０９、ステップＳ８１１、ステップＳ８１３、及びステップＳ８１５で算出されＲＡＭ１０３に記憶されている。
ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］＝Ａ×第１のビートレベル変動波形ＢＬ１［ｉｄｘＯｒｇ＿ｆ］＋Ｂ×第２のビートレベル変動波形ＢＬ２［ｉｄｘＯｒｇ＿ｆ］＋Ｃ×第３のビートレベル変動波形ＢＬ３［ｉｄｘＯｒｇ＿ｆ］＋Ｄ×第４のビートレベル変動波形ＢＬ４［ｉｄｘＯｒｇ＿ｆ］
・・・（８）
この演算処理により、ＣＰＵ１０１は、重み係数Ａ、Ｂ、Ｃ、及びＤを用いて、第１のビートレベル変動波形ＢＬ１［ｉｄｘＯｒｇ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＯｒｇ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＯｒｇ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＯｒｇ＿ｆ］に対する比較元重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］を計算する（ステップＳ１１０６）。ここで、重み係数Ａ、Ｂ、Ｃ、及びＤは、例えばＲＯＭ１０２に記憶されており、又は可変を許容する場合はＲＡＭ１０３に記憶されている。また、第１のビートレベル変動波形ＢＬ１［ｉｄｘＯｒｇ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＯｒｇ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＯｒｇ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＯｒｇ＿ｆ］は夫々、図８のフローチャートのステップＳ８０９、ステップＳ８１１、ステップＳ８１３、及びステップＳ８１５で算出されＲＡＭ１０３に記憶されている。
【0082】
同様に、ＣＰＵ１０１は、上述の（７）式の演算で算出された比較先フレームインデックスｉｄｘＤｓｔ＿ｆをキーとして、下記（９）式で示される演算処理を実行する。
ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］＝Ａ×第１のビートレベル変動波形ＢＬ１［ｉｄｘＤｓｔ＿ｆ］＋Ｂ×第２のビートレベル変動波形ＢＬ２［ｉｄｘＤｓｔ＿ｆ］＋Ｃ×第３のビートレベル変動波形ＢＬ３［ｉｄｘＤｓｔ＿ｆ］＋Ｄ×第４のビートレベル変動波形ＢＬ４［ｉｄｘＤｓｔ＿ｆ］
・・・（９）
この演算処理により、ＣＰＵ１０１は、上述の重み係数Ａ、Ｂ、Ｃ、及びＤを用いて、第１のビートレベル変動波形ＢＬ１［ｉｄｘＤｓｔ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＤｓｔ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＤｓｔ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＤｓｔ＿ｆ］に対する比較先重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］を計算する（同じくステップＳ１１０６）。ここで、第１のビートレベル変動波形ＢＬ１［ｉｄｘＤｓｔ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＤｓｔ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＤｓｔ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＤｓｔ＿ｆ］は夫々、図８のフローチャートのステップＳ８０９、ステップＳ８１１、ステップＳ８１３、及びステップＳ８１５で算出されＲＡＭ１０３に記憶されている。
ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］＝Ａ×第１のビートレベル変動波形ＢＬ１［ｉｄｘＤｓｔ＿ｆ］＋Ｂ×第２のビートレベル変動波形ＢＬ２［ｉｄｘＤｓｔ＿ｆ］＋Ｃ×第３のビートレベル変動波形ＢＬ３［ｉｄｘＤｓｔ＿ｆ］＋Ｄ×第４のビートレベル変動波形ＢＬ４［ｉｄｘＤｓｔ＿ｆ］
・・・（９）
この演算処理により、ＣＰＵ１０１は、上述の重み係数Ａ、Ｂ、Ｃ、及びＤを用いて、第１のビートレベル変動波形ＢＬ１［ｉｄｘＤｓｔ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＤｓｔ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＤｓｔ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＤｓｔ＿ｆ］に対する比較先重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］を計算する（同じくステップＳ１１０６）。ここで、第１のビートレベル変動波形ＢＬ１［ｉｄｘＤｓｔ＿ｆ］、第２のビートレベル変動波形ＢＬ２［ｉｄｘＤｓｔ＿ｆ］、第３のビートレベル変動波形ＢＬ３［ｉｄｘＤｓｔ＿ｆ］、及び第４のビートレベル変動波形ＢＬ４［ｉｄｘＤｓｔ＿ｆ］は夫々、図８のフローチャートのステップＳ８０９、ステップＳ８１１、ステップＳ８１３、及びステップＳ８１５で算出されＲＡＭ１０３に記憶されている。
【0083】
その後、ＣＰＵ１０１は、設定時間内カウンタ変数ｉを＋１インクリメントし、ステップＳ１１０２を介してステップＳ１１０３の処理に移行し、設定時間Ｔ内の次の位置ｉに対応する演算処理を繰り返す。
【0084】
やがて、設定時間内カウンタ変数ｉの値が設定時間Ｔに対応する末尾位置である設定時間サンプル数Ｎｕｍに到達し、ステップＳ１１０２でｉ≧Ｎｕｍとなってその判定がＹＥＳになると、ＣＰＵ１０１は、次の処理を実行する。ＣＰＵ１０１は、上述のように算出された設定時間Ｔに対応する比較元重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］及び比較先重み付き平均ビートレベル変動波形ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］（０≦ｉ＜Ｎｕｍ）に基づいて、既知の例えば下記（１０）式で示される自己相関演算により、相関値である相関係数ｃｏｒｒを算出する（ステップＳ１１０８）。
ｃｏｒｒ＝Ｃｏｖ（ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］，ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］）／（σ（ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］）×σ（ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］））
但し、（０≦ｉ＜Ｎｕｍ）
・・・（１０）
ここで、Ｃｏｖ（Ｘ，Ｙ）は、値ＸとＹの共分散を算出する関数演算である。また、σ（Ｘ）は、値Ｘの標準偏差を算出する関数演算である。
ｃｏｒｒ＝Ｃｏｖ（ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］，ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］）／（σ（ＷＡＭ＿ＢＬ＿Ｏｒｇ［ｉ］）×σ（ＷＡＭ＿ＢＬ＿Ｄｓｔ［ｉ］））
但し、（０≦ｉ＜Ｎｕｍ）
・・・（１０）
ここで、Ｃｏｖ（Ｘ，Ｙ）は、値ＸとＹの共分散を算出する関数演算である。また、σ（Ｘ）は、値Ｘの標準偏差を算出する関数演算である。
【0085】
その後、ＣＰＵ１０１は、比較先カウンタ変数ｋで示される現在のずれ時間「ｋ×ｄｏＤｓｔＳｔｅｐ」を有する比較元データと比較先データに対する、図１１のフローチャートで示される図９のステップＳ９０６の自己相関計算処理を終了する。
【0086】
図９のフローチャートの説明に戻り、ＣＰＵ１０１は、（１０）式で示される演算処理により計算された相関係数ｃｏｒｒを、比較先カウンタ変数ｋに対応する相関係数配列ｃｏｒｒ［ｋ］に格納し、ＲＡＭ１０３に記憶させる（ステップＳ９０７）。
【0087】
その後、ＣＰＵ１０１は、比較先カウンタ変数ｋの値を＋１インクリメントさせると共に、比較先先頭位置変数ｄｏＤｓｔの値に比較先時間ステップ幅変数ｄｏＤｓｔＳｔｅｐの値を累算する（ステップＳ９０８）。そして、ＣＰＵ１０１は、ステップＳ９０５からＳ９０６の処理に移行し、次の比較先データに対する自己相関計算処理を繰り返す。
【0088】
やがて、全ての比較先カウンタ変数ｋの値に対応する繰返し処理が終了しステップＳ９０５の判定がＹＥＳになると、ＣＰＵ１０１は、次の処理を実行する。ＣＰＵ１０１は、図７で前述したようにして、上述のステップＳ９０２～Ｓ９０８の繰返しにより得られる、比較元カウンタ変数ｎの値によって示される現在の経過時刻における比較元データに対して計算された自己相関波形において、例えば上位５個のピーク位置とそれらの相関値を算出し、夫々ＲＡＭ１０３に記憶される配列変数ＣｏｒｒＰｏｓＦｉｖｅ［ｊ］とＣｏｒｒＦｉｖｅ［ｊ］（０≦ｊ≦４）に格納する（ステップＳ９０９）。
【0089】
続いて、ＣＰＵ１０１は、図７で前述したようにして、ＲＡＭ１０３に記憶される配列変数である相関値のヒストグラムＨｉｓｔの、ステップＳ９０９で算出した上位５個の各ピーク位置ＣｏｒｒＰｏｓＦｉｖｅ［ｊ］に対応する各ビン位置に、ステップＳ９０９で算出した各相関係数ＣｏｒｒＦｉｖｅ［ｊ］を、夫々「Ｈｉｓｔ［ＣｏｒｒＰｏｓＦｉｖｅ［ｊ］］＋＝ＣｏｒｒＦｉｖｅ［ｊ］（０≦ｊ≦４）」として累算する（ステップＳ９１０）。
【0090】
その後、ＣＰＵ１０１は、比較元カウンタ変数ｎの値を＋１インクリメントさせると共に、比較元先頭位置変数ｄｏＯｒｇの値に比較元時間ステップ幅変数ｄｏＯｒｇＳｔｅｐの値を累算する（ステップＳ９１１）。そして、ＣＰＵ１０１は、ステップＳ９０２からＳ９０３の処理に移行し、次の比較元データに対する自己相関計算処理を繰り返す。
【0091】
やがて、全ての比較元データに対する処理が終了しステップＳ９０２の判定がＹＥＳになると、ＣＰＵ１０１は、図１０のステップＳ９１２以降の処理を実行する。
【0092】
まずＣＰＵ１０１は、ヒストグラムＨｉｓｔ［ｋ］から、ヒストグラム値が例えば上位７個までの大きさを有する相関値のヒストグラムのピーク位置を取得する（ステップＳ９１２）。ここで、ヒストグラムＨｉｓｔ［ｋ］は、前述の図９のステップＳ９０９及びステップＳ９１０の処理により得られている。また、ｋは、図４で説明したずらし時間を１．３３秒から４秒まで指定する比較先カウンタ変数値（前述の（１）式参照）である。
【0093】
次に、ＣＰＵ１０１は、ピーク比較元カウンタ変数ｎに対して、ステップＳ９１２で取得した例えば７個のピークのうちの最初のピーク番号＝０をセットする（ステップＳ９１３）。その後、ＣＰＵ１０１は、ステップＳ９２３で＋１ずつインクリメントしながら、ステップＳ９１４で全ての指定が終了したと判定するまで、ピーク比較元カウンタ変数ｎの値を順次指定する。そして、ＣＰＵ１０１は、このピーク比較元カウンタ変数ｎが示すピーク位置（＝ずらし時間長）を、ＲＡＭ１０３に記憶されるピーク比較元長さ変数ｌｅｎ１に格納する（ステップＳ９１５）。
【0094】
続いて、ＣＰＵ１０１は、ピーク比較元カウンタ変数ｎが指定されてピーク比較元が一つ指定される毎に、ピーク比較先カウンタ変数ｋに対して、ステップＳ９１２で取得した例えば７個のピークのうちの最初のピーク番号＝０をセットする（ステップＳ９１６）。その後、ＣＰＵ１０１は、ステップＳ９２２で＋１ずつインクリメントしながら、ステップＳ９１７で全ての指定が終了したと判定するまで、ピーク比較先カウンタ変数ｋの値を順次指定する。
そして、ＣＰＵ１０１はまず、このピーク比較先カウンタ変数ｋが示すピーク位置（＝ずらし時間長）を、ＲＡＭ１０３に記憶されるピーク比較先長さ変数ｌｅｎ２に格納する（ステップＳ９１８）。
そして、ＣＰＵ１０１はまず、このピーク比較先カウンタ変数ｋが示すピーク位置（＝ずらし時間長）を、ＲＡＭ１０３に記憶されるピーク比較先長さ変数ｌｅｎ２に格納する（ステップＳ９１８）。
【0095】
その後、ＣＰＵ１０１は、まずステップＳ９１９、Ｓ９２０、及びＳ９２１で、４拍子、３拍子、及び５拍子を順次仮定する。夫々の仮定のもとで、ＣＰＵ１０１は、ピーク比較元長さ変数ｌｅｎ１の値が小節時間長であると仮定する。そして、ＣＰＵ１０１は、ステップＳ９１８で順次指定されるピーク比較先長さ変数ｌｅｎ２の値を用いて計算される比率ｌｅｎ２／ｌｅｎ１が、前述した４拍子、３拍子、及び５拍子の何れかの倍数関係にあたるか否かを順次判定する。
【0096】
図１２は、図１０のステップＳ９１９の比較元の検証処理である４拍子検証処理の詳細例を示すフローチャートである。ＣＰＵ１０１は、４拍子の倍率を指定するＲＡＭ１０３に記憶される変数ｊを３に初期設定する（ステップＳ１２０１）。その後、ＣＰＵ１０１は、変数ｊの値を＋１ずつインクリメントしながら（ステップＳ１２０６）、ステップＳ１２０２で例えばピーク７個分の比較処理が終了したと判定するまで（図１０のステップＳ９１２を参照）、下記（１１）式で示される演算（ステップＳ１２０３）と、下記（１２）式で示される判定処理（ステップＳ１２０４）を実行する。
ＬＲａｔ［ｎ］＝１００×ｌｅｎ２／ｌｅｎ１－１００×ｊ／４ ・・・（１１）
－１≦ＬＲａｔ［ｎ］≦１ ・・・（１２）
ＬＲａｔ［ｎ］＝１００×ｌｅｎ２／ｌｅｎ１－１００×ｊ／４ ・・・（１１）
－１≦ＬＲａｔ［ｎ］≦１ ・・・（１２）
【0097】
（１１）式で示される演算処理により、ピーク比較元長さ変数ｌｅｎ１の値に対するピーク比較先長さ変数ｌｅｎ２の値の比率と、ステップＳ１２０３の繰返し処理により順次指定される４拍子の各倍率ｊ／４との差分値が算出される。そして、ステップＳ１２０４の繰返し処理により順次実行される（１２）式の判定処理が成立すると、比率ｌｅｎ２／ｌｅｎ１が現在の変数ｊの値に対応する４拍子の倍率ｊ／４に一致又はほぼ一致すると判定される。
【0098】
ステップＳ１２０４の判定がＹＥＳになったら、ＣＰＵ１０１は、ＲＡＭ１０３に記憶される変数ＴｅｍｐｏＯＫ［ｎ］の値をインクリメントする（ステップＳ１２０５）。なお、この変数値は、図１０のステップＳ９１６又はステップＳ９２２でピーク比較先カウンタ変数ｋの値が変更される毎に０（ゼロ）にリセットされる。
【0099】
ステップＳ１２０２の判定がＹＥＳとなって以上の繰返し処理が終了すると、ＣＰＵ１０１は、変数ＴｅｍｐｏＯＫ［ｎ］の値が２以上になっているか否か、即ち、ｌｅｎ２／ｌｅｎ１の比率が４拍子の倍率ｊ／４に一致又はほぼ一致した回数が２回以上あったか否かを判定する（ステップＳ１２０７）。
【0100】
ステップＳ１２０７の判定がＹＥＳならば、ＣＰＵ１０１は、現在仮定されている小節時間長が正しいと判定して、その４拍子の小節時間長を確定させ、同時にテンポも確定させる（ステップＳ１２０８）。
【0101】
ステップＳ１２０７の判定がＮＯならば、ＣＰＵ１０１は、ステップＳ１２０８の処理はスキップして、小節時間長及びテンポは確定させない。
【0102】
その後、ＣＰＵ１０１は、図１２のフローチャートで示される図１０のステップＳ９１９の比較元の検証処理を終了する。ここでＣＰＵ１０１は、上記ステップＳ１２０８を実行することにより小節時間長及びテンポが確定していれば、図９及び図１０のフローチャートで示される拍解析処理を終了して、確定した小節長及びテンポ値を、図１の表示部１０５に表示する。
【0103】
一方、ＣＰＵ１０１は、ステップＳ９１９の検証処理で小節時間長及びテンポが確定しなければ、次に、３拍子の検証処理を実行する（ステップＳ９２０）。図１３は、３拍子の検証処理の詳細例を示すフローチャートである。
【0104】
ＣＰＵ１０１は、３拍子の倍率を指定するＲＡＭ１０３に記憶される変数ｊを３に初期設定する（ステップＳ１３０１）。その後、ＣＰＵ１０１は、変数ｊの値を＋１ずつインクリメントしながら（ステップＳ１３０６）、ステップＳ１３０２で例えばピーク７個分の比較処理が終了したと判定するまで（図１０のステップＳ９１２を参照）、下記（１３）式で示される演算（ステップＳ１３０３）と、前述の（１２）式で示される判定処理（ステップＳ１３０４）を実行する。
ＬＲａｔ［ｎ］＝１００×ｌｅｎ２／ｌｅｎ１－１００×ｊ／３ ・・・（１３）
ＬＲａｔ［ｎ］＝１００×ｌｅｎ２／ｌｅｎ１－１００×ｊ／３ ・・・（１３）
【0105】
（１３）式で示される演算処理により、ピーク比較元長さ変数ｌｅｎ１の値に対するピーク比較先長さ変数ｌｅｎ２の値の比率と、ステップＳ１３０３の繰返し処理により順次指定される３拍子の各倍率ｊ／３との差分値が算出される。そして、ステップＳ１３０４の繰返し処理により順次実行される（１４）式の判定処理が成立すると、比率ｌｅｎ２／ｌｅｎ１が現在の変数ｊの値に対応する３拍子の倍率ｊ／３に一致又はほぼ一致すると判定される。
【0106】
ステップＳ１３０４の判定がＹＥＳになったら、ＣＰＵ１０１は、ＲＡＭ１０３に記憶される変数ＴｅｍｐｏＯＫ［ｎ］の値をインクリメントする（ステップＳ１３０５）。なお、この変数値は、図１０のステップＳ９１６又はステップＳ９２２でピーク比較先カウンタ変数ｋの値が変更される毎に０（ゼロ）にリセットされる。
【0107】
ステップＳ１３０２の判定がＹＥＳとなって以上の繰返し処理が終了すると、ＣＰＵ１０１は、変数ＴｅｍｐｏＯＫ［ｎ］の値が２以上になっているか否か、即ち、ｌｅｎ２／ｌｅｎ１の比率が３拍子の倍率ｊ／３に一致又はほぼ一致した回数が２回以上あったか否かを判定する（ステップＳ１３０７）。
【0108】
ステップＳ１３０７の判定がＹＥＳならば、ＣＰＵ１０１は、現在仮定されている小節時間長が正しいと判定して、その３拍子の小節時間長を確定させ、同時にテンポも確定させる（ステップＳ１３０８）。
【0109】
ステップＳ１３０７の判定がＮＯならば、ＣＰＵ１０１は、ステップＳ１３０８の処理はスキップして、小節時間長及びテンポは確定させない。
【0110】
その後、ＣＰＵ１０１は、図１３のフローチャートで示される図１０のステップＳ９２０の３拍子の検証処理を終了する。ここでＣＰＵ１０１は、上記ステップＳ１３０８を実行することにより小節時間長及びテンポが確定していれば、図９及び図１０のフローチャートで示される拍解析処理を終了して、確定した小節長及びテンポ値を、図１の表示部１０５に表示する。
【0111】
一方、ＣＰＵ１０１は、ステップＳ９２０の検証処理で小節時間長及びテンポが確定しなければ、引き続いて、５拍子の検証処理を実行する（ステップＳ９２１）。図１４は、５拍子の検証処理の詳細例を示すフローチャートである。
【0112】
ＣＰＵ１０１は、５拍子の倍率を指定するＲＡＭ１０３に記憶される変数ｊを３に初期設定する（ステップＳ１４０１）。その後、ＣＰＵ１０１は、変数ｊの値を＋１ずつインクリメントしながら（ステップＳ１４０６）、ステップＳ１４０２で例えばピーク７個分の比較処理が終了したと判定するまで（図１０のステップＳ９１２を参照）、下記（１４）式で示される演算（ステップＳ１４０３）と、前述の（１２）式で示される判定処理（ステップＳ１４０４）を実行する。
ＬＲａｔ［ｎ］＝１００×ｌｅｎ２／ｌｅｎ１－１００×ｊ／５ ・・・（１４）
ＬＲａｔ［ｎ］＝１００×ｌｅｎ２／ｌｅｎ１－１００×ｊ／５ ・・・（１４）
【0113】
（１３）式で示される演算処理により、ピーク比較元長さ変数ｌｅｎ１の値に対するピーク比較先長さ変数ｌｅｎ２の値の比率と、ステップＳ１４０３の繰返し処理により順次指定される５拍子の各倍率ｊ／５との差分値が算出される。そして、ステップＳ１４０４の繰返し処理により順次実行される（１４）式の判定処理が成立すると、比率ｌｅｎ２／ｌｅｎ１が現在の変数ｊの値に対応する５拍子の倍率ｊ／５に一致又はほぼ一致すると判定される。
【0114】
ステップＳ１４０４の判定がＹＥＳになったら、ＣＰＵ１０１は、ＲＡＭ１０３に記憶される変数ＴｅｍｐｏＯＫ［ｎ］の値をインクリメントする（ステップＳ１４０５）。なお、この変数値は、図１０のステップＳ９１６又はステップＳ９２２でピーク比較先カウンタ変数ｋの値が変更される毎に０（ゼロ）にリセットされる。
【0115】
ステップＳ１４０２の判定がＹＥＳとなって以上の繰返し処理が終了すると、ＣＰＵ１０１は、変数ＴｅｍｐｏＯＫ［ｎ］の値が２以上になっているか否か、即ち、ｌｅｎ２／ｌｅｎ１の比率が５拍子の倍率ｊ／５に一致又はほぼ一致した回数が２回以上あったか否かを判定する（ステップＳ１４０７）。
【0116】
ステップＳ１４０７の判定がＹＥＳならば、ＣＰＵ１０１は、現在仮定されている小節時間長が正しいと判定して、その５拍子の小節時間長を確定させ、同時にテンポも確定させる（ステップＳ１４０８）。
【0117】
ステップＳ１４０７の判定がＮＯならば、ＣＰＵ１０１は、ステップＳ１４０８の処理はスキップして、小節時間長及びテンポは確定させない。
【0118】
その後、ＣＰＵ１０１は、図１４のフローチャートで示される図１０のステップＳ９２１の５拍子の検証処理を終了する。ここでＣＰＵ１０１は、上記ステップＳ１４０８を実行することにより小節時間長及びテンポが確定していれば、図９及び図１０のフローチャートで示される拍解析処理を終了して、確定した小節長及びテンポ値を、図１の表示部１０５に表示する。
【0119】
一方、ＣＰＵ１０１は、ステップＳ９２１の検証処理で小節時間長及びテンポが確定しなければ、ステップＳ９２２でピーク比較先カウンタ変数ｋの値を＋１インクリメントする。その後ＣＰＵ１０１は、ステップＳ９１７からステップＳ９１８の処理に移行して、ステップＳ９１２で取得した例えば７個のピークのうちの次のピーク番号に対する上述の処理を繰り返し実行する。
【0120】
やがて、ステップＳ９１２で取得した例えば７個のピークの全てに対して上述の処理が終了してステップＳ９１４の判定がＹＥＳになると、ＣＰＵ１０１は、図１の表示部１０５に、小節長及びテンポを確定できなかった旨をエラー表示して、図９及び図１０のフローチャートで例示される拍解析処理を終了する（ステップＳ９２４）。
【0121】
図１５は、小節線位置特定処理の例を示すフローチャートである。図９から図１４のフローチャートを用いて説明した拍解析処理を実行すると、楽曲全体に対するその楽曲のテンポと拍子と小節時間（小節長）が抽出される。今、上記拍解析処理によって決定された小節時間をｍｅａｓＬｅｎ（単位は「秒」）とし、拍子を例えば４拍子とする。また、楽曲中の小節の順番（何番目の小節か）をｍｅａｓＮｕｍとし、ｍｅａｓＮｕｍ番目の小節の楽曲の先頭からの暫定的な先頭位置（経過時刻、単位は「秒」）をｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］とする。この結果、ｍｅａｓＬｅｎとｍｅａｓＮｕｍを使って、ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］は、下記（１５）式で示される演算により決定される。
ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＝ｍｅａｓＬｅｎ×ｍｅａｓＮｕｍ
・・・（１５）
ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＝ｍｅａｓＬｅｎ×ｍｅａｓＮｕｍ
・・・（１５）
【0122】
上記（１５）式で決定される小節の暫定的な先頭位置ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］は、あくまでも暫定値である。今、小節の正しい先頭位置を小節線位置と呼ぶことにすれば、小節線位置は、時間経過に伴うテンポの変動に伴って各ビートのピーク位置が変化することによって、小節の暫定的な先頭位置ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］からずれる。この小節線ずれ量をｂｅｓｔＰｈａｓｅとすれば、小節線位置は、下記（１６）式で示される演算により決定される。
小節線位置＝ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＋ｂｅｓｔＰｈａｓｅ
・・・（１６）
図１５のフローチャートで示される小節線位置特定処理は、上記（１６）式の小節線ずれ量ｂｅｓｔＰｈａｓｅを算出して、正しい小節線位置を特定するための処理である。
小節線位置＝ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＋ｂｅｓｔＰｈａｓｅ
・・・（１６）
図１５のフローチャートで示される小節線位置特定処理は、上記（１６）式の小節線ずれ量ｂｅｓｔＰｈａｓｅを算出して、正しい小節線位置を特定するための処理である。
【0123】
小節線位置を特定するために、図１５に示されるフローチャートにおいて、ＣＰＵ１０１はまず、小節線ずれ量仮定値ｍｅａｓＰｈａｓｅ（単位は「秒」）を、ステップＳ１５０１で０．０（ゼロ）に初期設定する。その後、ＣＰＵ１０１は、ステップＳ１５０９で例えば０．００５秒（５ミリ秒）ずつ順次インクリメントしながら、ステップＳ１５０２で所定の最大値に達して終了したと判定するまで、ステップＳ１５０２からＳ１５０９までの一連の処理を繰り返し実行する。
【0124】
上記繰返し処理において、ＣＰＵ１０１は更に、ステップＳ１５００３において小節番号ｍｅａｓＮｕｍと誤差合計値ｄｏＶａｌを夫々初期値０（ゼロ）に設定する。その後、ＣＰＵ１０１は、ステップＳ１５０６で小節番号ｍｅａｓＮｕｍを順次インクリメントしながら、ステップＳ１５０４で楽曲の最後の小節に達したと判定するまで、以下に説明するステップＳ１５０５及びＳ１５０６の各処理を実行する。
【0125】
今、現在の小節番号ｍｅａｓＮｕｍと小節線ずれ量仮定値ｍｅａｓＰｈａｓｅに対応する現在評価している小節線位置（現在小節線位置）（単位は「秒」）は、前述の（１６）式と同様にして、下記（１７）式で示される演算により決定される。
現在小節線位置＝ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＋ｍｅａｓＰｈａｓｅ
・・・（１７）
現在小節線位置＝ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＋ｍｅａｓＰｈａｓｅ
・・・（１７）
【0126】
また、現在の小節番号ｍｅａｓＮｕｍに対応する小節を１６分音符単位の位置に分割したときに、先頭の１６分音符単位位置は、上記（１７）式で算出される現在小節線位置に等しい。この現在小節線位置をデジタルのサンプル数で示したものをｉｄｘ［０］とすれば、このｉｄｘ［０］は、上記（１７）式を用いて、下記（１８）式で示される演算により算出される。
ｉｄｘ［０］＝ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＋ｍｅａｓＰｈａｓｅ
×サンプリング周波数［Ｈｚ］ ・・・（１８）
ｉｄｘ［０］＝ｍｅａｓＴｉｍｅ［ｍｅａｓＮｕｍ］＋ｍｅａｓＰｈａｓｅ
×サンプリング周波数［Ｈｚ］ ・・・（１８）
【0127】
また、小節内の１６分音符のサンプル長をｉｄｘ１６とすれば、このｉｄｘ１６は、前述した小節時間ｍｅａｓＬｅｎを用いて、下記（１９）式で示される演算により決定される。
ｉｄｘ１６＝ｍｅａｓＬｅｎ／１６×サンプリング周波数［Ｈｚ］
・・・（１９）
ｉｄｘ１６＝ｍｅａｓＬｅｎ／１６×サンプリング周波数［Ｈｚ］
・・・（１９）
【0128】
上記（１８）式及び（１９）式より、現在の小節番号ｍｅａｓＮｕｍに対応する小節を１６分音符単位の位置に分割したｉｄｘ［０］以外の各サンプル位置ｉｄｘ［ｉ］（１≦ｉ≦１５）は、下記（２０）式で示される演算により決定される。
ｉｄｘ［ｉ］＝ｉｄｘ［ｉ－１］＋ｉｄｘ１６×ｉ ・・・（２０）
ｉｄｘ［ｉ］＝ｉｄｘ［ｉ－１］＋ｉｄｘ１６×ｉ ・・・（２０）
【0129】
更に、現在の小節番号ｍｅａｓＮｕｍに対応する小節を１６分音符単位の位置に分割した各サンプル位置ｉｄｘ［ｉ］（０≦ｉ≦１５）をフレーム番号に変換したフレーム位置ｉｄｘ＿ｆ［ｉ］（０≦ｉ≦１５）は、下記（２１）式で示される演算により決定される。ただし、ｆｓｉｚｅはフレームサイズ（単位は「サンプル」）である。
ｉｄｘ＿ｆ［ｉ］＝ｉｄｘ［ｉ］／ｆｓｉｚｅ ・・・（２１）
ｉｄｘ＿ｆ［ｉ］＝ｉｄｘ［ｉ］／ｆｓｉｚｅ ・・・（２１）
【0130】
以上の（１８）式から（２１）式で示される演算処理により算出されるｉｄｘ＿ｆ［ｉ］（０≦ｉ≦１５）を用いて、下記（２２）式から（２５）式の各演算処理が実行される。これらの演算処理では、現在の小節番号ｍｅａｓＮｕｍに対応する小節を１６分音符単位の位置に分割した各フレーム位置で抽出されるビートレベルＢＬ１、ＢＬ２、ＢＬ３、及びＢＬ４に基づいて、小節線ずれ量仮定値ｍｅａｓＰｈａｓｅに対して算出されるその小節内の１６分音符単位のビートレベル列が算出される。ここで、ＢＬ１、ＢＬ２、ＢＬ３、及びＢＬ４は夫々、前述した図８のフローチャートで示される処理により、全帯域、ＢＤ帯域、ＳＤ帯域、及びコード帯域の各周波数帯域で抽出されるビートレベルである。
｛ＢＬ１［ｉｄｘ＿ｆ［０］］，ＢＬ１［ｉｄｘ＿ｆ［１］］，
ＢＬ１［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ１［ｉｄｘ＿ｆ［１５］］｝
・・・（２２）
｛ＢＬ２［ｉｄｘ＿ｆ［０］］，ＢＬ２［ｉｄｘ＿ｆ［１］］，
ＢＬ２［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ２［ｉｄｘ＿ｆ［１５］］｝
・・・（２３）
｛ＢＬ３［ｉｄｘ＿ｆ［０］］，ＢＬ３［ｉｄｘ＿ｆ［１］］，
ＢＬ３［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ３［ｉｄｘ＿ｆ［１５］］｝
・・・（２４）
｛ＢＬ４［ｉｄｘ＿ｆ［０］］，ＢＬ４［ｉｄｘ＿ｆ［１］］，
ＢＬ４［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ４［ｉｄｘ＿ｆ［１５］］｝
・・・（２５）
｛ＢＬ１［ｉｄｘ＿ｆ［０］］，ＢＬ１［ｉｄｘ＿ｆ［１］］，
ＢＬ１［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ１［ｉｄｘ＿ｆ［１５］］｝
・・・（２２）
｛ＢＬ２［ｉｄｘ＿ｆ［０］］，ＢＬ２［ｉｄｘ＿ｆ［１］］，
ＢＬ２［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ２［ｉｄｘ＿ｆ［１５］］｝
・・・（２３）
｛ＢＬ３［ｉｄｘ＿ｆ［０］］，ＢＬ３［ｉｄｘ＿ｆ［１］］，
ＢＬ３［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ３［ｉｄｘ＿ｆ［１５］］｝
・・・（２４）
｛ＢＬ４［ｉｄｘ＿ｆ［０］］，ＢＬ４［ｉｄｘ＿ｆ［１］］，
ＢＬ４［ｉｄｘ＿ｆ［２］］，・・・，ＢＬ４［ｉｄｘ＿ｆ［１５］］｝
・・・（２５）
【0131】
一方、周波数帯域（全帯域、ＢＤ帯域、ＳＤ帯域、及びコード帯域）の夫々に対応して、夫々が例えば下記（２６）式に例示されるような、１６分音符で１６音符分即ち１小節分のビートを表したビートパターンが用意される。
｛１，０．１，０．３，０．１，
０．５，０．１，０．３，０．１，
０．７，０．１，０．３，０．１，
０．３，０．１，０．３，０．１｝ ・・・（２６）
上記（２６）式において、夫々の行に対応する１６分音符で４音符分の４つずつの数字は、例えば４拍子の場合の各拍を構成する１６分音符×４音符のビートの強さを示している。これらのビートの強さは、最大値が１になるように正規化されている。
｛１，０．１，０．３，０．１，
０．５，０．１，０．３，０．１，
０．７，０．１，０．３，０．１，
０．３，０．１，０．３，０．１｝ ・・・（２６）
上記（２６）式において、夫々の行に対応する１６分音符で４音符分の４つずつの数字は、例えば４拍子の場合の各拍を構成する１６分音符×４音符のビートの強さを示している。これらのビートの強さは、最大値が１になるように正規化されている。
【0132】
例えば、（２６）式の第１行目の４つの数字は、小節内の第１拍目に対応している。その第１拍目において、最初の数字「１」は、第１拍目の１６分音符分解能の１番目（先頭）において、最大ビートが刻まれることを示している。それに続く３つの数字「０．１」「０．３」「０．１」は、第１拍目の１６分音符分解能の２番目、３番目、４番目において、非常に小さな振幅のビートが刻まれることを示している。
【0133】
また例えば、（２６）式の第３行目の４つの数字は、小節内の第３拍目に対応している。その第３拍目において、最初の数字「０．７」は、第３拍目の１６分音符分解能の１番目（先頭）において、大きな振幅のビートが刻まれることを示している。この値は、第１拍目に次いで大きな振幅のビートとなっている。即ち、この例では、第１拍目と第３拍目が、いわゆる強拍に対応していることがわかる。
【0134】
一方例えば、小節内の第２拍目、第４拍目に対応する（２６）式の第２行目及び第４行目の１６分音符分解能の１番目（先頭）の値は夫々、「０．５」及び「０．３」という比較的小さな振幅のビートとなっており、この例では、第２拍目と第４拍目が、いわゆる弱拍に対応していることがわかる。
【0135】
本実施形態では、上記（２６）式に例示されるようなビートパターンが、４つの周波数帯域（全帯域、ＢＤ帯域、ＳＤ帯域、及びコード帯域）の夫々に対応して、全部で４パターン用意される。
【0136】
その後、ステップＳ１５０５において、現在の小節番号ｍｅａｓＮｕｍに対応する現在の小節で、各周波数帯域（全帯域、ＢＤ帯域、ＳＤ帯域、及びコード帯域）に対応して前述の（２２）式から（２５）式で示されるように算出された１６分音符単位の各ビートレベル列の夫々について、（２６）式の形式で用意された４つのビートパターンの夫々に対して、最小二乗誤差が計算される。具体的には、周波数帯域毎に、その周波数帯域に対して（２２）式から（２５）式の何れかにより算出されたビートレベル列とその周波数帯域に対して用意されたビートパターンの各１６分音符単位のデータの差を二乗した値を１６個分加算したものが、その周波数帯域に対する最小二乗誤差として計算される。
【0137】
そして、ステップＳ１５０５において、周波数帯域毎に上述のようにして計算された最小二乗誤差が、４つの周波数帯域分累算されて、その累算結果が変数ｄｏＶに格納される。
【0138】
その後、ステップＳ１５０６において、ステップＳ１５０５にて現在の小節番号ｍｅａｓＮｕｍが示す小節に対して計算された最小二乗誤差ｄｏＶが、楽曲全体の最小二乗誤差累算値を表す変数ｄｏＶａｌに累算される。
【0139】
以上のステップＳ１５０５及びＳ１５０６の各処理が、楽曲全体の小節にわたって実行されることにより、ステップＳ１５０４で楽曲の末端の小節が処理されたと判定された後、ステップＳ１５０８において、上記処理により算出された現在の小節線ずれ量仮定値ｍｅａｓＰｈａｓｅに対応する楽曲全体の最小二乗誤差累算値ｄｏＶａｌが、現在までに得られている誤差最小値ｍｉｎ（初期値はステップＳ１５０１で大きな値に設定されている）よりも小さいか否かが判定される。
【0140】
ステップＳ１５０８の判定がＹＥＳならば、ステップＳ１５０４からＳ１５０７の一連の処理により今回算出された現在の小節線ずれ量仮定値ｍｅａｓＰｈａｓｅに対応する楽曲全体の最小二乗誤差累算値ｄｏＶａｌが、新たな誤差最小値ｍｉｎに設定される。それと共に、今回の小節線ずれ量仮定値ｍｅａｓＰｈａｓｅが新たな小節線ずれ量最適値ｂｅｓｔＰｈａｓｅに設定される。
【0141】
以上の制御処理が、小節線ずれ量仮定値ｍｅａｓＰｈａｓｅが順次更新されながら繰返し実行されることにより、ステップＳ１５０２の判定がＹＥＳとなった時点で、最終的な小節線ずれ量最適値ｂｅｓｔＰｈａｓｅが算出される。そして、この小節線ずれ量最適値ｂｅｓｔＰｈａｓｅを用いて、前述した（１６）式により、各小節番号ｍｅａｓＮｕｍに対応する小節線位置が決定される。
【0142】
図１５に関する上記説明では、前述した（２２）式から（２５）式で示される演算により算出された小節内の１６分音符単位のビートレベル列が使用されたが、小節を１６分音符単位で分割した各小範囲毎にビートレベル値を累算して得られる１６分音符単位のビートレベル列を使用することで、より安定した小節線位置特定処理を実行することができる。
【0143】
以上説明した実施形態において、図１１のステップＳ１１０６の重み付き平均処理で、重み係数Ａ、Ｂ、Ｃ、及びＤの組合せを変えながら、相関が最も高い結果を採用するように制御してもよい。
【0144】
以上の実施形態に関して、更に以下の付記を開示する。
（付記１）
楽曲データの或るタイミングにおける第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、重み付きビートレベルを取得し、
取得された前記重み付きビートレベルの自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
拍子決定方法。
（付記２）
前記重み付きビートレベルに基づいて、前記楽曲データの或るタイミングから設定時間長の比較元と、前記或るタイミングから離れたタイミングから前記設定時間長の比較先を取得し、
取得された前記比較元と前記比較先に基づいて、前記自己相関を計算する、
付記１に記載の拍子決定方法。
（付記３）
前記第１ビートレベルは、前記第１パワーレベルの変動がマイナスからプラスに転じた第１タイミングのパワーレベルと、前記変動が次にプラスからマイナスに転じた第２タイミングのパワーレベルとのレベル差を含み、
前記レベル差を前記第２タイミングにおける前記第１ビートレベルとして取得する、
付記１又は２に記載の拍子決定方法。
（付記４）
前記重み付きビートレベルに基づいて、前記自己相関値のヒストグラムを生成し、
生成された前記ヒストグラムに基づいて、前記拍子を決定する、
付記１乃至３の何れかに記載の拍子決定方法。
（付記５）
前記ヒストグラムから抽出した、値が上位所定数個までの相関値のヒストグラムのピークの夫々が、３拍子、４拍子又は５拍子のいずれの拍子の倍数関係にあるかを判定し、
判定結果に基づいて、前記拍子を決定する、
付記４に記載の拍子決定方法。
（付記６）
前記第１周波数帯域は、全帯域であり、
前記第２周波数帯域は、バスドラム帯域、スネアドラム帯域又はコード帯域の何れかを含む、
付記１乃至５の何れかに記載の拍子決定方法。
（付記７）
楽曲データの或るタイミングにおける第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、重み付きビートレベルを取得し、
取得された前記重み付きビートレベルの自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
処理を実行する拍子決定装置。
（付記８）
コンピュータが、
楽曲データの或るタイミングにおける第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、重み付きビートレベルを取得し、
取得された前記重み付きビートレベルの自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
処理を実行するためのプログラム。
（付記１）
楽曲データの或るタイミングにおける第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、重み付きビートレベルを取得し、
取得された前記重み付きビートレベルの自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
拍子決定方法。
（付記２）
前記重み付きビートレベルに基づいて、前記楽曲データの或るタイミングから設定時間長の比較元と、前記或るタイミングから離れたタイミングから前記設定時間長の比較先を取得し、
取得された前記比較元と前記比較先に基づいて、前記自己相関を計算する、
付記１に記載の拍子決定方法。
（付記３）
前記第１ビートレベルは、前記第１パワーレベルの変動がマイナスからプラスに転じた第１タイミングのパワーレベルと、前記変動が次にプラスからマイナスに転じた第２タイミングのパワーレベルとのレベル差を含み、
前記レベル差を前記第２タイミングにおける前記第１ビートレベルとして取得する、
付記１又は２に記載の拍子決定方法。
（付記４）
前記重み付きビートレベルに基づいて、前記自己相関値のヒストグラムを生成し、
生成された前記ヒストグラムに基づいて、前記拍子を決定する、
付記１乃至３の何れかに記載の拍子決定方法。
（付記５）
前記ヒストグラムから抽出した、値が上位所定数個までの相関値のヒストグラムのピークの夫々が、３拍子、４拍子又は５拍子のいずれの拍子の倍数関係にあるかを判定し、
判定結果に基づいて、前記拍子を決定する、
付記４に記載の拍子決定方法。
（付記６）
前記第１周波数帯域は、全帯域であり、
前記第２周波数帯域は、バスドラム帯域、スネアドラム帯域又はコード帯域の何れかを含む、
付記１乃至５の何れかに記載の拍子決定方法。
（付記７）
楽曲データの或るタイミングにおける第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、重み付きビートレベルを取得し、
取得された前記重み付きビートレベルの自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
処理を実行する拍子決定装置。
（付記８）
コンピュータが、
楽曲データの或るタイミングにおける第１周波数帯域の第１パワーレベルに応じた第１ビートレベルと、前記第１周波数帯域に含まれる第２周波数帯域の第２パワーレベルに応じた第２ビートレベルに基づいて、重み付きビートレベルを取得し、
取得された前記重み付きビートレベルの自己相関値が高い複数のタイミングを取得し、
取得された前記複数のタイミングに基づいて、拍子を決定する、
処理を実行するためのプログラム。
【符号の説明】
【0145】
１０１ ＣＰＵ
１０２ ＲＯＭ
１０３ ＲＡＭ
１０４ 入力部
１０５ 表示部
１０６ 出力部
２０１、２０２、２０３ パワー変動
３０１ パワー変動波形
３０２ ビートレベル変動波形
４０１ 重み付き平均ビートレベル変動波形
４０２ 比較元データ
４０３ 比較先データ
５０１、６０１ 相関値のピーク
７００ 相関値のヒストグラム
７０１、７０２ 相関値のヒストグラムのピーク
１０２ ＲＯＭ
１０３ ＲＡＭ
１０４ 入力部
１０５ 表示部
１０６ 出力部
２０１、２０２、２０３ パワー変動
３０１ パワー変動波形
３０２ ビートレベル変動波形
４０１ 重み付き平均ビートレベル変動波形
４０２ 比較元データ
４０３ 比較先データ
５０１、６０１ 相関値のピーク
７００ 相関値のヒストグラム
７０１、７０２ 相関値のヒストグラムのピーク
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】