(57)【特許請求の範囲】
【請求項１】
音響波形を入力する入力手段と、
前記入力された音響波形のサンプルデータについて、所定の時間範囲を対象に自己相関を求める第１の自己相関手段と、
前記第１の自己相関手段の計算結果に基づき、前記入力された音響波形の属する音域を判定する音域判定手段と、
前記音域判定手段で判定した音域に対応する時間範囲を対象に、前記音響波形について自己相関を求める第２の自己相関手段と、
前記第２の自己相関手段の計算結果に基づき音名を決定する音名決定手段と
を具えることを特徴とする音名検出器。
【請求項２】
入力される音響波形の音名を検出する処理をコンピュータに実行させるためのソフトウェアプログラムであって、
音響波形を入力するステップと、
前記入力された音響波形のサンプルデータについて、所定の時間範囲を対象に自己相関を求める第１の自己相関ステップと、
前記第１の自己相関ステップの計算結果に基づき、前記入力された音響波形の属する音域を判定する音域判定ステップと、
前記音域判定ステップで判定した音域に対応する時間範囲を対象に、前記音響波形について自己相関を求める第２の自己相関ステップと、
前記第２の自己相関手段の計算結果に基づき音名を決定する音名決定ステップと
を含むことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、外部から入力される音響信号のピッチを測定し、その音響波形の音名を検出する音名検出器及びプログラムである。
【背景技術】
【０００２】
周知の通り、調律器はピアノ等の自然楽器の調律を行う際に利用される装置である。調律器として従来から知られるものには、楽器から発音された音響信号の音高（ピッチ）を検出し、当該音響信号の音名及びオクターブと、該音響信号の音高が調律の基準となる音高（調律対象周波数）からどれだけずれているかを表示する機能を有するものがある。この種の調律器において音高を測定する方法には、次のような方法が従来から知られる。すなわち、入力された音響信号（入力波形）をそのゼロクロス点毎に「１」又は「０」に反転してなる矩形パルス波に整形することで２値化した参照波形を作成し、該作成した参照波形を前記各ゼロクロス点の時間間隔ずつ遅延させた遅延波形を作成し、参照波形と遅延波形の相関を求めて、相関が高い遅延量を入力波形の周期とすることで、入力波形の音高を測定する方法が知られる（例えば下記特許文献１参照）。
また、上記参照波形と遅延波形の相関を求めるピッチ測定方法の改良として、所定の閾値に基づき入力波形のレベルの上限ピーク値及び下限ピーク値を判定し、参照波形と遅延波形との間でピーク値の存在点の相関を求める処理を行い、参照波形と遅延波形の相関と、ピーク値存在点の相関との双方の相関が高い遅延波形の遅延量を入力波形の周期とする方法があった（例えば、下記特許文献２参照）。
【特許文献１】特公平３−４２４１２号公報
【特許文献２】特開平９−２５７５５８号公報
【０００３】
しかし、上記特許文献１又は特許文献２のいずれに開示された方法においても、入力波形をそのゼロクロス点毎に「１」又は「０」に反転してなる矩形パルス波に整形することで２値化した参照波形を作成しているので、ゼロクロス点付近で生じたノイズに大きく影響されてしまい、ゼロクロスの周期決定が難しく、ピッチ測定が困難であるという不都合があった。また、入力波形に強い倍音成分が含まれる場合（基音が弱い場合）には、その倍音の周期で高い相関がとれてしまう等、ピッチ測定の精度が低下してしまう。また、低音域では、参照波形と遅延波形の遅延量が大きくなるため、上記の方法ではピッチ測定が困難になり、ピッチ測定精度が不十分だった。
【０００４】
上記特許文献２においては、ピッチ測定の別の方法として、入力波形全体の自己相関を求め、これにより入力波形のピッチを測定する方法が存在することが記述されている。しかしながら、入力波形全体の自己相関を求める方法では、同文献２において指摘されている通り、非常に膨大な計算量を必要とする。従って、この方法によって、或る楽器の全音域にわたり高い精度の自己相関関数の計算を行うことは実用的ではなかった。このことから、この種のピッチ測定方法を適用した製品（調律器）においては、例えば低音域のピッチ測定ができない等の性能上の制約があった。
【０００５】
【発明の開示】
【発明が解決しようとする課題】
【０００６】
この発明は、基音のレベルが弱い場合や低音域の楽音等であっても高精度なピッチ測定を行えるようにすると共に、ピッチ測定のための計算量を少なくした音名検出器を提供することを目的とする。
【課題を解決するための手段】
【０００７】
この発明は、音響波形を入力する入力手段と、前記入力された音響波形のサンプルデータについて所定の時間範囲を対象に自己相関を求める第１の自己相関手段と、前記第１の自己相関手段の計算結果に基づき、前記入力された音響波形の属する音域を判定する音域判定手段と、前記音域判定手段で判定した音域に対応する時間範囲を対象に、前記音響波形について自己相関を求める第２の自己相関手段と、前記第２の自己相関手段の計算結果に基づき音名を決定する音名決定手段とを具えることを特徴とする音名検出器である。
【０００８】
また、入力される音響波形の音名を検出する処理をコンピュータに実行させるためのソフトウェアプログラムであって、音響波形を入力するステップと、前記入力された音響波形のサンプルデータについて、所定の時間範囲を対象に自己相関を求める第１の自己相関ステップと、前記第１の自己相関ステップの計算結果に基づき、前記入力された音響波形の属する音域を判定する音域判定ステップと、前記音域判定ステップで判定した音域に対応する時間範囲を対象に、前記音響波形について自己相関を求める第２の自己相関ステップと、前記第２の自己相関手段の計算結果に基づき音名を決定する音名決定ステップとを含むことを特徴とするプログラムとして構成してもよい。
【発明の効果】
【０００９】
この発明によれば、入力された音響波形のサンプルデータについて、所定の時間範囲を対象に自己相関を求める第１の自己相関手段が、所定の時間範囲のみを対象にした自己相関の計算を行い、該第１の自己相関手段の計算結果に基づき、音域判定手段は音響波形の属する音域を判定する。前記所定の時間範囲は、音域判定手段において音域の判定が可能な計算結果を得られる範囲であればよく、これは極短い時間範囲に設定できる。従って、第１の自己相関手段における計算量は極めて少量であってよい。第２の自己相関手段では、音域判定手段で判定した音域に対応する時間範囲を対象に音響波形について自己相関を求め、音名決定手段はこの第２の自己相関手段の計算結果に基づき音名を決定する。第２の自己相関手段においては、音域判定手段により既にピッチを測定すべき音域が絞り込まれているので、限られた音域のみを対象にして自己相関の計算を行うことができる。従って、第２の自己相関手段の自己相関処理もまた、限られた処理対象範囲のみについて計算すればよいことになるので、その計算量を削減することができる。自己相関によるピッチ測定は、ノイズに対して比較的安定したピッチ測定を行うことができ、また、基音のレベルが弱い場合や低音域の楽音等においても比較的高精度なピッチ測定が可能である。従って、この発明によれば、自己相関処理の計算量を大幅に少なくし、安定性並びに精度の高いピッチ測定を効率的に行うことができるという優れた効果を奏する。
【発明を実施するための最良の形態】
【００１０】
以下添付図面を参照して、この発明の一実施例ついて説明する。以下の実施例では、この発明に係る音名検出器を、アコースティックピアノを調律するための調律器の形態で構成及び実施する例について説明する。
【００１１】
図１は、この発明の一実施例に係る調律器の使用時の概要を説明するための概念図である。図１において符号１はこの実施例に係る調律器であり、符号２は調律すべき楽器の一例としてアコースティックな発音機構を具えた鍵盤楽器（アップライトピアノ）を示している。調律器１は、例えばＰＤＡ（Ｐｅｒｓｏｎａｌ  Ｄｉｇｉｔａｌ  Ａｓｓｉｓｔａｎｔｓ：情報携帯端末）で構成され、この実施例に係る調律機能の実現するためのソフトウェアプログラムを実行可能である。また、調律器１は、周知のＰＤＡに一般的であるのと同様に、情報の表示及び画面上からの操作入力が可能な表示器（例えば液晶スクリーン等）３を具える。調律器１には、ピアノ２において発音されたピアノ演奏音を取得するためのマイクロフォン４が接続されている。該マイクロフォン４によってピックアップした音が調律器１に対する入力波形として供給される。詳しくは後述する通り、この実施例に係る調律器１は、マイクロフォン４を介して供給される入力波形のピッチ（音高）を検出し、該検出したピッチに基づく音名を表示器３に表示する機能を有する。また、調律器１は、マイクロフォン４を介して供給される入力波形の周波数と比較基準用の周波数（調律対象周波数）の位相関係つまり周波数もしくは周期関係を示す画像（表示パターン）を表示器３に表示する機能を有する。なお、図１においては、図示及び説明の便宜上、調律器１に対してマイク４を外部接続するように描かれているが、マイク４は調律器１に内蔵される構成であっても差し支えない。
以下の説明において、最初に入力波形の周波数と比較基準用の周波数（調律対象周波数）の位相関係つまり周波数もしくは周期関係を示す画像（表示パターン）を表示する表示制御について述べて、次に、ピッチ検出処理についての説明を行うものとする。
【００１２】
図２は調律器１の電気的ハードウェア構成の概略を示すブロック図である。調律器１は、ＣＰＵ１０、ＲＯＭ１１、ＲＡＭ１２を含むマイクロコンピュータと、マイクロフォン４から供給される入力波形を取り込むためのインターフェース１３、表示制御部１４及び操作検出部１５とから構成され、各装置間がバス１６を介して接続される。インターフェース１３は、アンプ及びＡＤ変換器を含み、マイク４を介して入力されアナログオーディオ波形信号を適宜増幅及びディジタル波形信号に変換して信号処理系に供給する。表示制御部１４はＣＰＵ１０からの指示に基づき表示器３の表示を制御する。ＰＤＡにおける表示器３の表示は、ソフトウェアによって概ね１５〜２０Ｈｚの周期で更新されており、この周期は当該調律器１に入力されうるピアノ音の最低音の周波数よりも遅いものである。また、操作検出部１５を介して表示器３上での操作入力やその他の操作子の入力に応じた指示がＣＰＵ１に与えられる。
ＣＰＵ１０は調律器１の全体的な動作を制御すると共に、この実施例に係る調律機能、すなわち、音名検出機能及び入力波形の周波数と調律対象周波数のずれに応じた表示パターンを表示する機能を実現するためのソフトウェアプログラムを実行する。前記各機能を実現するためのソフトウェアプログラムは例えばＲＯＭ１１等のメモリに記憶されていてよい。ＲＡＭ１２は、ＣＰＵ１０が信号処理を実行する際に、各種パラメータや各種データを記憶するワークエリア、あるいは、マイク４を介して取り込んだ波形データを格納するメモリ領域として使用される。
【００１３】
図３（ａ），（ｂ）は、表示器３に読み出される表示画面の一例を示す図である。表示画面３０には、調律対象周波数と入力波形の周波数とのずれ（両者の位相関係）を示す画像「表示パターン３２」を表示するための波形表示部３１が設けられている。図３（ａ）は調律対象周波数と入力波形（ピアノ音）の周波数が一致している場合の表示パターン３２を示し、同図（ｂ）は調律対象周波数と入力波形の周波数がずれている場合の表示パターン３２を示している。なお、この表示パターン３２の表示制御の原理や表示動作については後述する。
【００１４】
また、表示画面３０には、数値や音名（Ａ〜Ｆ）等の情報を入力するためのボタン画像群からなる操作入力部３３が具わる。ユーザは、画面上に表示されたボタン画像にペン等の操作ツールを当てて、操作入力することができる。ユーザは、操作入力部３３から調律対象となる音高の指定や、各種パラメータの設定等を行う。表示画面３０の設定値表示部３４には、調律対象となる音高等の各種パラメータの設定値が表示される。
設定値表示部３４に表示されるパラメータには、図に示すとおり、音名（ｎｏｔｅ）、オクターブ（ｏｃｔ．）、キー番号（ｋｅｙＮｏ．）、セント値（ｃｅｎｔ）、及び、調律対象周波数（ｆｒｅｑ）がある。ユーザは、操作入力部３３のボタン画像等を用いて音名／数値入力により音名（ｎｏｔｅ）及びオクターブ（ｏｃｔ．）の値を指定できる。キー番号（ｋｅｙＮｏ．）は、周知の通り、ピアノ鍵盤に具わる８８鍵の各鍵に対して、最低音の鍵から最高音の鍵の順に１〜８８の番号を順次割り当てられた、各鍵に固有の番号であって、音名とオクターブにより音高が特定されれば、当該音高に対応するキー番号が一義に定まる。すなわち、ユーザは、音名及びオクターブの指定により、調律対象とする調律対象音高を指定することができる。調律対象周波数は、基本的には、該指定された調律対象音高に対応する周波数に相当するものである。音高に対応する正規の周波数は、各音高毎に予め設定された所定値であり、例えばテーブル参照等により設定されてよい。以下、この実施例において、該調律対象周波数の値を示すパラメータを「Ｈｚ値」と称する。また、セント値は、前記Ｈｚ値をユーザが任意に修正するためのパラメータである。従って、設定値表示部３４の調律対象周波数（ｆｒｅｑ）の欄には、指定された調律対象音高の周波数をセント値で修正したＨｚ値が表示されることになる。例えば、同図に示すとおり、指定された音高が第５オクターブのＡ（ｋｅｙＮｏ．＝４９）で、セント値が０であれば、調整対象周波数は第５オクターブのＡに相当する周波数（この例では４４０Ｈｚ）になる。なお、セント値は、周知の通り、音高差を対数表現で表した値であって、１００セントが１２音階の半音に相当する。
【００１５】
図４及び図５は、表示パターン３２の表示制御の原理を説明するための概念図である。図４（ａ），（ｂ）において、符号４０ａ，４０ｂで示す波形は調律器１に入力された入力波形（ピアノ音）である。この実施例に係る表示パターン３２の表示制御においては、入力された入力波形から、調律対象周波数Ｈｚの周期毎に、後述する窓幅パラメータ「Ｗ値」に基づく時間範囲（位相範囲）にわたる波形区間の波形を抽出する処理を行い、該抽出された各波形区間毎に、該区間内の波形の周期性を明示する表示情報を作成する。前記表示情報は、当該波形をそのゼロクロス点毎に「１」又は「０」に反転してなる２値化情報（後述の図７を参照）であり、該区間毎の波形の周期性を２段階の濃度で画像化しうるデータである。そして、表示画面３０（表示器３）の表示更新周期の１周期内で作成される前記各波形区間毎の各表示情報を重ね合わせて、該重ね合わせた多値の情報に応じた表示形態の表示パターン３２を表示器３に表示させる。以上がこの実施例に係る表示パターンの表示制御の原理の大略である。
【００１６】
図４（ａ）において、（１）は調律対象周波数Ｈｚに一致する周波数の入力波形４０ａが入力された場合の波形図を示しており、（２）は該（１）の入力波形４０ａから調律対象周波数Ｈｚの周期ごとに抽出される波形区間毎の波形４１ａ〜４１ｅであって、同図において該各波形区間毎の波形４１ａ〜４１ｅは、各波形のゼロクロス点毎に「１」又は「０」に反転した２値の表示情報によって各波形の周期性を２段階の濃度で画像化（周期パターン）されたものにより示されている。入力波形４０ａの周波数と調律対象周波数Ｈｚが一致しているので、各波形区間の波形４１ａ〜４１ｅの抽出は、入力波形４０ａの周期に同期したタイミングで行われることになる。従って、各波形区間の波形４１ａ〜４１ｅとして、各々の位相関係が一致つまり同じ周期特性を持つものが抽出される。これら各波形区間の波形４１ａ〜４１ｅの周期パターンの画像を重ね合わせると、図５（ａ）に示すような表示パターン３２ａを得ることができる。各波形４１ａ〜４１ｅの周期性が一致しているので、これらを重ね合わせた表示パターン３２ａは、各波形４１ａ〜４１ｅと同じ周期性を持つ２段階の濃淡画像となることが図５（ａ）から見て取れる。
一方、図４（ｂ）において、（３）は調律対象周波数Ｈｚと異なる周波数の入力波形４０ｂが入力された場合（Ｈｚ値と入力波形の周波数がずれている場合）の波形図を示し、（４）は該（３）の入力波形４０ｂから調律対象周波数Ｈｚの周期ごとに抽出される波形区間毎の波形４１ｆ〜４１ｊを２段階の濃度で画像化した周期パターンを示している。この場合、入力波形４０ｂと調律対象周波数Ｈｚがずれているので、各波形区間の波形４１ｆ〜４１ｊとして、それぞれ位相がずれた周期パターンが抽出されることになる。これら各波形区間の波形４１ｆ〜４１ｊの周期パターンの画像を重ね合わせると、図５（ｂ）に示すような表示パターン３２ｂを得ることができる。各波形４１ｆ〜４１ｊの位相関係がずれているので、これらを重ね合わせた表示パターン３２ｂは、該各波形４１ｆ〜４１ｊの位相関係のずれに応じた多段階（少なくとも３段階以上）の濃淡を持つ多値の濃淡画像となる。
【００１７】
図５（ａ），（ｂ）に示す表示パターン３２ａ，３２ｂは、表示更新周期の１周期内で抽出された各波形区間毎の波形４１ａ〜４１ｅ，４１ｆ〜４１ｊを重ね合わせたものである（図４（ａ），（ｂ）参照）。すなわち、この実施例では、表示更新周期の１周期内で可能な数（図の例では５つ）だけ波形区間毎の波形の抽出を行い、該表示更新周期の１周期毎に、抽出した各波形区間毎の波形の表示情報を重ね合わせて、表示パターン３２の表示態様を更新する。図５（ａ）のように、入力波形４０ａの周波数と調律対象周波数Ｈｚが一致している場合には、各表示更新周期毎に更新される表示パターン３２ａは、同じ周期性を持つ２段階の濃淡画像となるので、見かけ上、表示パターン３２ａの表示態様で停止しているように見える。これに対して、図５（ｂ）のように、入力波形４０ｂの周波数と調律対象周波数Ｈｚがずれている場合には、各表示更新周期毎に更新される表示パターン３２ｂの周期性が一定でないため、見かけ上、表示パターン３２ｂが波形表示部３１（図３参照）上を流れるように見え、表示態様が一定状態で停止することがない。また、表示器の表示更新周期が入力波形４０ｂの周波数の周期と調律対象周波数Ｈｚの周期と公倍数の関係になる場合に、表示パターン３２ｂが停止しているように見えてしまったとしても、表示パターン３２ｂは多段階（少なくとも３段階以上）の濃淡を持つ多値の濃淡画像によって入力波形４０ｂの周波数と調律対象周波数Ｈｚのずれを明示することができるようになる。なお、表示パターン３２の表示制御の動作の詳細については後述する。
【００１８】
入力波形から各波形区間を抽出する際の時間範囲（位相範囲）すなわち入力波形に窓をかける幅（窓幅）は、窓幅パラメータ「Ｗ値」に基づき設定される。Ｗ値は、調律対象周波数Ｈｚの波形を何周期分抽出するかによって、各波形区間毎の波形を抽出する時間範囲（位相範囲）を規定するパラメータである。例えば、図４乃至図５に示す例では、波形区間の時間範囲として、調律対象周波数Ｈｚの波形を２．５波形分抽出するようＷ値が設定されている。Ｈｚ値の周期と入力波形４０ａの周期が一致していれば、波形４１ａ〜４１ｅとしてちょうど２．５周期分の波形が抽出され、これらを重ね合わせた表示パターン３２ａとして２．５周期分の周期パターンの画像が表示される（図４乃至図５の（ａ）参照）。また、Ｈｚ値の周期と入力波形４０ｂの周期がずれている場合には、一致していなければ、図４乃至図５の（ｂ）に示す通り、抽出される各波形４１ｆ〜４１ｊは２．５周期分の波形にならない。
Ｗ値の設定を変更することで、波形表示部３１における表示パターン３２（図３参照）の表示サイズを変更可能である。各波形区間毎の波形を抽出する時間範囲を小さくすれば、波形表示部３１における表示パターン３２の表示サイズは相対的に大きくなることになる。表示パターン３２の表示サイズを拡大することで、ユーザはＨｚ値の周期と入力波形４０ａの周期の微小なずれまで視認できるようになり、より精密な調律が行えるようになるという利点がある。
【００１９】
次に、この実施例に係る調律器１において調律作業を実行するための作業手順等について、図６、図７及び図９のフローチャートを参照して説明する。
調律器１の電源投入に応じて、調律器１では図６のフローチャートに示す処理が実行される。図６のステップＳ１では各種パラメータの初期設定を行っている。ここで初期設定されるパラメータは、基準ピッチ、調律対象周波数（Ｈｚ値）、セント値及びＷ値等である。各パラメータの初期設定値は、例えば、基準ピッチ＝４４０Ｈｚ、調律対象周波数Ｈｚ＝４４０Ｈｚ、セント値＝０ｃｅｎｔ、Ｗ値＝Ｈｚ値の２．５周期分、にそれぞれ設定されるものとする。ステップＳ２において、作業者は、当該調律器１において音階を設定するための基準ピッチ（スタンダードピッチ）を入力する。すなわち、基準ピッチとなるピアノ中央のＡ音（キー番号４９のＡ）の周波数を、４４０Ｈｚ、４４２Ｈｚ或るいは４３９Ｈｚ等、幾つかの候補うちからを任意に選択することができる。続いて、調律器１では、表示画面３０（図３参照）上の波形表示部３１において表示パターン３２の表示が開始されると共に、図７に示す表示制御の動作（タスク）が開始される（ステップＳ３及びステップＳ４）。なお、図６では、図示及び説明の便宜上、ステップＳ３の「表示開始」及びステップＳ４の「タスク開始」を別々のステップとして描いているが、双方とも、実質的には、後述図７及び図９の処理の開始に照応している。
ステップＳ５においては、調律カーブの選択を受け付けている。調律カーブは、ピアノの８８鍵の各鍵に対応する音高の周波数を定めたデータテーブルである。調律カーブとして、ピアノの種類（グランドピアノ／アップライトピアノ）や、大きさ等に応じた複数種類のデータテーブルがＲＯＭ２乃至ＲＡＭ３等適宜のメモリ内に記憶されおり、作業者は調律するピアノの種類に応じた調律カーブを選択できてよい。該調律カーブに記述された各音高の周波数は、ピアノの特性に鑑みて、高音側のピッチが平均率による周波数の理論値よりも高めに設定されている。上記ステップＳ２において基準ピッチとして設定された周波数に基づき、平均率による各音高の周波数の理論値を計算することはできる。しかし、実際のピアノの調律においては、その理論値をそのまま各音高の周波数として適用するのは不適当な場合がある。その場合には、調律カーブを使用することで、ピアノの特性や、使用するピアノの種類或いは大きさ等に適った各音高毎の周波数を得ることができる。
【００２０】
上記ステップＳ１〜Ｓ５は、調律器１の使用に際しての初期設定に相当する処理である。すなわち、これら各ステップＳ１〜Ｓ５の実行順序は図示の順に限らない。ユーザは、以下に述べるステップＳ６〜Ｓ１１により各種パラメータの設定を行う。
【００２１】
ユーザの操作入力（前記図３等を参照）により、調律対象の音高が指定された場合には（ステップＳ６のｙｅｓ）、ステップＳ７において、前記ステップＳ５で選択した調律カーブ又は平均律による理論値に基づき、指定された音高の周波数を求めて、この値を調律対象周波数のパラメータ「Ｈｚ値」に設定する。また、ユーザの操作入力により、セント値が入力された場合には（ステップＳ８のｙｅｓ）、ステップＳ９において、該入力されたセント値に応じて調律対象周波数Ｈｚを修正する。また、ステップＳ１０では、Ｗ値の変更による表示パターンの表示サイズの選択を受け付けている。この実施例によれば、ユーザは表示パターンの表示サイズとして通常のサイズと拡大サイズのいずれかを選択できる。すなわち、Ｗ値が変更された場合（ステップＳ１０のｙｅｓ）、ステップＳ１１では前記変更に応じてＷ値を設定する。Ｗ値が変更されると、入力波形から波形区間を抽出する時間範囲が変更されるので、表示パターンの表示サイズを任意に変更できる。表示サイズ変更の使い方としては、最初に、通常サイズの表示を見ながら大まかに調律を行った後、表示サイズを拡大して微妙な周波数のずれを見る場合等に使用できる。
以降、調律器１の電源がオンされている間は、ステップＳ６〜Ｓ１１を繰り返すことで、ユーザによるＨｚ値変更及び表示サイズ変更を受け付けることができるようになっている。
【００２２】
図７は、上記図６のステップＳ３及びＳ４において起動開始する表示パターン表示制御の動作の手順を示すフローチャートである。
図７に示す処理は、当該調律器１の表示器３（表示画面３０）の表示更新周期に応じた起動タイミング毎に起動するタイマ処理であり、この処理の起動機会毎に、波形表示部３１の表示パターン３２（図３参照）が更新される。この実施例では一例として調律器１をＰＤＡで構成する例を示している。前述の通りＰＤＡにおける表示器３の表示は、ソフトウェアによって概ね１５〜２０Ｈｚの周期で更新される。従って、当該処理もまた１５〜２０Ｈｚ周期程度で起動する。
ステップＳ２０において、マイクロフォン４を介して入力された入力波形（ピアノ音）を所定のサンプリング周期でサンプリングすることでピアノ音を調律器１に取り込む。サンプリング周波数は例えば４４．１ｋＨｚとする。ステップＳ２１では、各サンプリングタイミング毎にサンプリングされた入力波形（ディジタル波形信号）をＲＡＭ１２上のメモリ領域に書き込む。ステップＳ２２では、前記ＲＡＭ１２上のメモリ領域に所定サンプル数以上の波形信号のサンプルデータが読み込まれたかどうかを判断する。波形信号のサンプルデータの読み込みが該所定サンプル数に満たない場合（ステップＳ２２のｎｏ）は、波形信号のサンプルデータが該所定サンプル数に達するまで、波形信号の取り込み（上記ステップＳ２０及びＳ２１）を繰り返す。一方、ＲＡＭ１２上のメモリ領域に読み込んだ波形信号のサンプルデータが所定サンプル数に達したら（ステップＳ２２のｙｅｓ）、次のステップＳ２３に処理を進める。なお、前記所定サンプル数は、概ね、１０２４〜２０４８サンプル程度とする。
【００２３】
ステップＳ２３において、前記図５のステップＳ７又はＳ９において設定された調律対象周波数のパラメータＨｚ値に基づきバンドパスフィルタのフィルタ係数（パスする帯域幅及び中心周波数）を設定し、ステップＳ２４において、ＲＡＭ２４のメモリ領域に読み込んだ波形に対して、前記設定したフィルタ係数によりバンドパスフィルタ処理を施す。このフィルタ処理により、入力波形に含まれる倍音成分等を除去し、調律対象周波数の成分を抽出することができる。
【００２４】
ステップＳ２５において、前記図５のステップＳ１１で設定されたＷ値に基づき、波形区間毎の波形として何周期分の波形を抽出するか（窓幅）を設定する。Ｗ値が２．５周期分に設定されていれば、窓幅は調律対象周波数Ｈｚの周期の２．５周期分に相当する時間範囲になる。
【００２５】
ステップＳ２６において、入力波形から、調律対象周波数Ｈｚ値の周期毎に、前記ステップＳ２５で設定した窓幅の波形区間の波形を抽出する（窓をかける）。これにより、Ｗ値に応じた時間範囲の波形区間毎の波形が、表示更新周期の１周期内でＨｚ値の周期が到来する回数に応じた複数区間分抽出される。
ステップＳ２７では、抽出された各波形区間毎に、該区間内の波形を、そのゼロクロス点毎に値が"０箔１"で反転する２値の情報に整形することで、その周期性を２値の濃淡によって明示する表示情報（「２値化情報」）を作成する。図８（ａ），（ｂ）において、抽出された各波形区間毎の波形の２値化情報（ゼロクロス点毎に値が"０箔１"で反転する２値の情報）の構成例を示すと共に、各２値化情報の右側に、該各２値化情報を２段階の濃淡からなる周期パターン化したものを示す。図８（ａ），（ｂ）に示す通り、２値化情報は、その各アドレスが当該波形区間における位相（時間）に対応付けられたデータで構成されている。すなわち、同図の例では、２値化情報は、Ｈｚ値の周期の２．５サイクル分の時間範囲を２５段階のアドレスに分割し、各アドレスの値"０"又は"１"により、当該波形区間内の波形の周期性を２段階の濃度（「０」または「１」）で表現するデータである。同図（ａ）は入力波形とＨｚ値の周期が一致している場合、（ｂ）は入力波形とＨｚ値の周期がずれている場合を示しており、図示及び説明の便宜上、各波形区間毎の波形に対して、図４，５と同様に、符号４１ａ〜４１ｅ（図８（ａ））、符号４１ｆ〜４１ｊ（図８（ｂ））を付与している。
【００２６】
ステップＳ２８では、表示更新周期の１周期内で作成した前記各波形区間毎の２値化情報の相加平均を求めることで、該各２値化情報を重ね合わせた多値の情報（多値化情報）を作成する。そして、ステップＳ２９では、前記ステップＳ２８において作成された多値化情報に基づく表示形態で、波形表示部３１に表示パターン３２を表示させる。表示パターン３２は、多値化情報に基づく濃淡画像を輝度乃至色相の違いによって示す画像である（図３（ａ），（ｂ）の表示画面３０を参照）。
図８（ａ）に示す通り、入力波形の周波数と調律対象周波数Ｈｚが一致している場合には、２値化情報４１ａ〜４１ｅの相加平均である多値化情報４２ａは、各波形４１ａ〜４１ｅと同じ周期性を持つ２値（図においては"０"又は"５"）のデータとなる。従って、この多値化情報４２ａに基づき表示される表示パターン３２ａは２段階の濃淡画像となる。一方、図８（ｂ）に示す通り、調律対象周波数Ｈｚ値と入力波形の周波数が不一致の場合には、各波形区間毎の２値化情報４１ｆ〜４１ｊの相加平均として求められる多値化情報４２ｂは多段階（３段階以上；図においては"１"，"２"，"３"又は"４"）の値を持つ情報となる。従って、この多値化情報４２ｂに基づき表示される表示パターン３２ｂは、多段階（３段階以上）の濃淡変化を持つ濃淡画像となる。
【００２７】
上記ステップＳ２０〜Ｓ２９の処理を表示画面３０の表示更新周期に応じて実行することで、表示パターン３２の表示形態が表示更新周期ごとに更新される。これにより、調律対象周波数Ｈｚと入力波形の周波数が一致している場合（図８（ａ）参照）には、表示更新周期毎に、周期性の一致した２段階の濃淡を持つ表示パターン３２ａが繰り返し表示されるため、見かけ上、２段階の濃淡からなる表示形態の画像が一定位置で停止しているように見えるようになる。これに対して、調律対象周波数Ｈｚと入力波形の周波数が不一致の場合（図８（ｂ）参照）、表示パターン３２は大略以下の２通りの見え方になる。すなわち、各表示更新周期毎に更新される表示パターン３２ｂの周期性が一定しない（各各表示更新周期毎に更新される表示パターンの位相がずれる）ことから、見かけ上、表示パターン３２ｂが波形表示部３１（図３参照）上を流れているように見え、その表示態様が一定状態で停止しない。或いは、表示更新周期が入力波形の周波数の周期と調律対象周波数Ｈｚの周期と公倍数の関係になる場合には、入力波形の周波数と調律対象周波数Ｈｚが不一致であっても、表示パターン３２ｂが一定位置に停止しているように見えてしまうが、この実施例によれば、表示パターン３２ｂのように、濃淡画像の濃淡が多段階（３段階以上）で表示されることにより、入力波形の周波数と調律対象周波数Ｈｚがずれていることを明示できる。
【００２８】
次に、調律器１における音名検出（ピッチ測定）機能について図９及び図１０のフローチャートを参照して説明する。
図９に示す音名検出処理は、所定の起動周期毎に繰り返し実行される処理であり、該起動周期は比較的遅いサイクル（１秒間に数回程度）であってよい。図９のステップＳ３０において、マイクロフォン４を介して入力された入力波形（ピアノ音）を所定のサンプリング周期でサンプリングすることでピアノ音を調律器１に取り込む。入力波形を取り込む際のサンプリング周波数は例えば４４．１ｋＨｚとする。ステップＳ３１では、各サンプリングタイミング毎にサンプリングされた入力波形のレベルが所定の音量レベル（音量閾値）以上であるかどうか判断する。入力波形のレベルが前記音量閾値以下であれば（ステップＳ３１のｎｏ）、以下の処理を行わずにリターンする。一方、音量閾値以上の波形が入力されていれば（ステップＳ３１のｙｅｓ）、ステップＳ３２において入力波形（ディジタル信号）をＲＡＭ１２上のメモリ領域に書き込む。そして、ステップＳ３３において、該メモリ領域に書き込んだ入力波形に対して図１０に示す自己相関処理を実行する。
【００２９】
周知の通り、１つの波形信号を自己相関処理することで、該波形信号自身の持つ周期性を調べることができる。すなわち、或る入力波形ｘ（ｋ）について自己相関Ｒ（ｍ）を求め、自己相関処理の結果、高い相関を得ることができた遅延量ｍを該入力波形ｘ（ｋ）の周期として推定することができる。ここで、自己相関処理によるピッチ測定をピアノ演奏音に特化する場合には、測定対象とするピッチの範囲は、ピアノの８８鍵の各々に対応するピッチ範囲に限定されて差し支えない。従って、ピアノ演奏音のピッチ測定を行う場合には、最大でもピアノの８８鍵の各々の周波数の周期（周波数の逆数）に対応する各時刻について自己相関関数Ｒ（ｍ）を計算すればよい。すなわち、従来から知られる方法では、ピアノの８８鍵の各々の周波数の周期に対応する８８点の時刻を遅れ時間を変数ｍに設定して、この８８点の各遅れ時間について自己相関の計算を行う必要があった。
これに対して、この実施例によれば、先ず第１の自己相関手段による自己相関の計算の結果に基づき入力波形の音域を大雑把に判定し、第２の自己相関手段において前記判定された音域のみについてピッチ測定のための自己相関を求めることで、自己相関処理の計算量を大幅に減少せしめ、効率的にピッチ測定（音名検出）が行えるようになる。この実施例に係る自己相関処理の具体的な動作の一例については図１０等を参照して以下に説明する。
【００３０】
図１０は前記ステップＳ３３における自己相関処理の手順の一例を示すフローチャートである。ステップＳ４０、ステップＳ４１及びＳ４２において、前記図９のステップＳ３２においてメモリ領域に取り込んだ入力波形ｘ（ｋ）について、所定の時間範囲を対象に自己相関関数Ｒ（ｍ）を求める処理を行う（第１の自己相関手段）。
ステップＳ４０では、図１１に示すデータテーブル５０を参照して、自己相関の計算に使用する変数ｍ（遅れ時間）を値の小さいものから順に設定する。データテーブル５０は、計算に使用する遅れ時間ｍとして所定の時間範囲内の複数ポイントの時刻のデータから構成されるもので、ＲＯＭ１１又はＲＡＭ１２等適宜のメモリに記憶されていてよい。このデータテーブル５０に記憶されている各遅れ時間が「第１の自己相関手段」における自己相関処理の対象となる時間範囲に対応する。この実施例では一例としてｍ１＝６，ｍ２＝１２，ｍ３＝２５，ｍ４＝５０の４点の遅れ時間のデータがデータテーブル５０に記憶されている。なお、この実施例において遅れ時間ｍはミリ秒（ｍｓｅｃ）単位で表現されたデータとする。
ステップＳ４１では、ステップＳ４０で設定された遅れ時間ｍについて下記に示す計算式１により入力波形ｘ（ｋ）について自己相関関数Ｒ（ｍ）を求める計算を行う。
【数１】
0004470823.tif 000002
この計算式１は、Ｍ個のサンプリングデータからなる或る入力波形ｘ（ｋ）と、それを遅れ時間ｍだけずらした入力波形ｘ（ｋ−ｍ）に対して、それぞれかけ合せて累積し、平均化するという処理を表している。この自己相関処理の計算結果は適宜のバッファメモリに記録する。計算式１による計算は後述する通り自己相関関数Ｒ（ｍ）の特性の大まかな傾向が判れば良いため、当該計算に使用するサンプルデータ数Ｍは５００〜１０００サンプル程度であってよい。
ステップＳ４２では、データテーブル５０のすべての遅れ時間ｍ１〜ｍ４について自己相関を計算を行ったかどうか調べ、未だすべての変数ｍについて計算が終わっていなければ（ステップＳ４２のｎｏ）、ステップＳ４０に戻り、テーブル５０を参照して次の遅れ時間ｍを設定すると共に、該設定された変数ｍについて計算式１による計算を行うという処理を繰り返す。このように第１の自己相関手段では、上記ステップＳ４０〜Ｓ４２の処理を各遅れ時間ｍ（この例ではｍ１〜ｍ４の４点）毎に行うことで、所定の時間範囲（ｍ１〜ｍ４の範囲）を処理対象範囲として入力波形についての自己相関Ｒ（ｍ）を求める。
【００３１】
図１２（ａ）はステップＳ４０〜Ｓ４２における自己相関処理の計算結果を模式的に示す図である。同図（ａ）において縦軸に自己相関関数Ｒ（ｍ）の値をとり、横軸に遅れ時間ｍをとる。また、入力波形が低い音であった場合の自己相関関数Ｒ（ｍ）の特性の一例を実線で示し、高い音であった場合の自己相関関数Ｒ（ｍ）の特性の一例を破線で示している。図１２（ａ）において、処理対象の時間範囲（ｍ１〜ｍ４）における自己相関関数Ｒ（ｍ）を見てみると、低音（実線）の場合には当該時間範囲においてはＲ（ｍ）が正の値のままであり、また、高音（破線）の場合には時間範囲の最後（ｍ４＝５０ｍｓｅｃ）の時点ではＲ（ｍ）が負の値になっている。このことから明らかなように、ステップＳ４０〜Ｓ４２において計算式１により処理対象の時間範囲（ｍ１〜ｍ４）の４点の遅れ時間ｍについて自己相関を求めることで、自己相関関数Ｒ（ｍ）の特性の大まかな傾向から遅れ時間ｍの値が小さいうちに入力波形の属する音域を判断することが可能である。
ステップＳ４３では、上記ステップＳ４０〜Ｓ４２による自己相関処理（第１の自己相関手段）の計算結果にづき、最後の変数ｍ４についての計算結果Ｒ（ｍ）の値まで負であれば（ステップＳ４３のｙｅｓ）、当該入力波形の音域を高音域と判断し、ステップＳ４４の「高音域用計算式２」による自己相関処理を行う。また、同計算結果Ｒ（ｍ）の値が正であれば（ステップＳ４３のｎｏ）、当該入力波形の音域を低音域と判断し、ステップＳ４５の「低音域用計算式２」による自己相関処理を行う。ステップＳ４４又はＳ４５における自己相関処理が第２の自己相関手段に相当する。
【００３２】
ステップＳ４４では、下記に示す高音域用計算式２により入力波形ｘ（ｉ）について自己相関関数Ｒ（ｍ）を求める計算を行う。
【数２】
0004470823.tif 000003
高音域用計算式２では、波形の周期が短いものと断定できるので、計算に使用する入力波形ｘ（ｉ）のサンプルデータ数Ｍ1を小さい値に設定してよい。この実施例では一例としてサンプルデータ数Ｍ1＝５１２とする。変数ｍには高音域に対応する各鍵に割り当てられた音高周波数の逆数（周期）に対応する各値が設定される。この実施例において、例えば鍵盤中央に位置するキー番号４５の鍵から最高音のキー番号８８の鍵の４４鍵に対応する音域を「高音域」とする。すなわち、この実施例では、第４オクターブの音名Ｆ〜第８オクターブの音名Ｃ（最高音）の音域が「高音域」となる。従って、変数ｍとしては、前記高音域の各鍵に割り当てられた音高（第４オクターブの音名Ｆ〜第８オクターブの音名Ｃ）の周波数の逆数、すなわち、変数ｍ＝１０．５４，１１．１６，１１．８３・・・１１２．５，１１９．２，１２６．３の４４個の各値が設定されうる。言い換えれば、高音域は変数ｍの値が比較的小さい範囲である。なお、ここで挙げた変数ｍの数値例は、ピアノの基準ピッチ（第４オクターブの音名Ａ）の周波数を４４０Ｈｚとし、当該調律器１のサンプリングレートを４４．１ｋＨｚとした場合のものであって、変数ｍはｍ＝（１／周波数）＊４４．１ｋなる演算によって算出できる（但し「ｋ」は１０００を表す）。
上記計算式２により、Ｍ1個（５１２個）のサンプリングデータからなる入力波形ｘ（ｉ）を上記４４個の各遅れ時間ｍについて自己相関処理を行うことで、高音域に属する各音高周波数（第４オクターブの音名Ｆ〜第８オクターブの音名Ｃ）の各周期に対応する各遅れ時間ｍについて自己相関関数Ｒ（ｍ）を求める。
【００３３】
ステップＳ４５では、下記に示す低音域用計算式３により入力波形ｘ（ｊ）について自己相関関数Ｒ（ｍ）を求める計算を行う。
【数３】
0004470823.tif 000004
この低音域用計算式３では、「Ｍ2−ｍ／２」とすることで計算に使用する入力波形ｘ（ｊ）のサンプルデータ数Ｍ2を小さい値に設定できるよう工夫している。この実施例では一例としてサンプルデータ数Ｍ2＝２０４８とする。変数ｍには低音域に対応する各鍵に割り当てられた音高周波数の逆数（周期）に対応する各値が設定される。この実施例において、例えばキー番号１の最低音の鍵から鍵盤中央に位置するキー番号４４の鍵までの４４鍵に対応する音域を「低音域」とする。すなわち、この実施例では、第０オクターブの音名Ａ（最低音）〜第４オクターブの音名Ｅの音域が「低音域」となる。従って、変数ｍとしては、前記低音域の各鍵に割り当てられた音高（第０オクターブの音名Ａ〜第４オクターブの音名Ｅ）の周波数の逆数、すなわち、変数ｍ＝１３３．８，１４１．７，１５０．２・・・１４２８，１５１３，１６０３の４４個の各値が設定されうる。言い換えれば、低音域は変数ｍの値が比較的大きい範囲である。なお、ここで挙げた変数ｍの数値例の算出条件は上記と同様である。
上記計算式３により、Ｍ2個（２０４８個）のサンプリングデータからなる入力波形ｘ（ｊ）を上記４４個の各遅れ時間ｍについて自己相関処理を行うこととで、低音域に属する各音高周波数（第０オクターブの音名Ａ〜第４オクターブの音名Ｅ）の各周期に対応する各遅れ時間ｍについて自己相関関数Ｒ（ｍ）を求める。
【００３４】
ステップＳ４６では、上記ステップＳ４４又はステップＳ４５（第２の自己相関手段）による計算結果として出力された自己相関関数Ｒ（ｍ）のうちで、Ｒ（ｍ）が極大値となった遅れ時間ｍを求める。図１２（ｂ）はステップＳ４４における高音域用計算式２の計算結果、同図（ｃ）はステップＳ４５における低音域用計算式３の計算結果を、それぞれ模式的に示す図であって、両図において縦軸は自己相関関数Ｒ（ｍ）、横軸は遅れ時間ｍをそれぞれ示している。Ｒ（ｍ）が極大値になるポイントは、最も相関が高い遅れ時間ｍである。従って、Ｒ（ｍ）が極大値となった遅れ時間ｍは、当該入力波形の周期に対応する時間であると判断することができる。予め上記ステップＳ４０〜ステップＳ４３の処理により入力波形の音域を大まかに判別しているので、後段の計算（ステップＳ４４又はステップＳ４５）においては、高音域又は低音域のいずれかに対応する処理範囲の各変数ｍについて自己相関計算をするだけで、Ｒ（ｍ）が極大値になるポイントを抽出できる。
【００３５】
ステップＳ４７では、前記ステップＳ４６で求めたＲ（ｍ）が極大値となった遅れ時間ｍに基づき当該入力波形の音高（ピッチ）を決定する。或る波形信号の自己相関を求め、相関が高い遅延量を該波形信号の周期に推定できるのは周知の通りである。前述の通り、遅れ時間ｍはピアノの各鍵に割り当てれた音高の周波数の逆数に対応する値に設定されている。従って、遅れ時間ｍの値を周波数に換算すれば、当該遅れ時間ｍの周期に相当するピアノ鍵に割り当てれた音高の周波数を得ることができる。
【００３６】
ステップＳ４８では、前記ステップＳ４７で決定された音高（ピッチ）に対応する音名を表示する制御を行う。音名の表示は例えば表示画面３０の設定値表示部３４（図３参照）における音名（ｎｏｔｅ）欄になされるものとする。音名のみならずオクターブ（ｏｃｔ．）やキー番号（ｋｅｙＮｏ．）を併せて表示してもよい。
【００３７】
以上説明した通り、この実施例によれば、先ずステップＳ４０〜Ｓ４２の処理（第１の自己相関手段）において、所定の時間範囲の数点（この実施例では４点）の変数ｍについて自己相関処理を行い、ステップＳ４３において前記計算の結果に基づき入力波形の音域（高音域又は低音域のいずれか）を判定することで、極少ない計算量で入力波形の大雑把に音域別に振り分ける。その後、ステップＳ４４又はＳ４５の処理（第２の自己相関手段）により、前記音域のいずれかについてピッチ測定のための自己相関処理を行う。ステップＳ４４又はＳ４５のいずれの計算においても、ピッチ測定の対象となる音域を限定しているので、計算に使用するデータ数Ｍ1乃至Ｍ2や変数ｍの数を少なくすることができる。よって、第２の自己相関手段の計算量も削減することができる。従って、この実施例によれば、自己相関処理の計算量を大幅に少なくすることで、比較的ノイズに強く、基音が弱い場合等にも正確な測定が可能な精度の高いピッチ測定を効率的に行うことができるようになる。
【００３８】
なお、上記自己相関処理における遅れ時間ｍはサンプルデータの参照位置をいくつずらすかに対応している。このことは該自己相関処理が所定サンプリングレート（上記の例では４４．１ｋＨｚ）でサンプリングされたディジタル信号を扱っていることから理解される。従って、上記図１０のステップＳ４２，ステップＳ４４又はステップＳ４５の計算において、遅れ時間ｍが小数点以下を含む値の場合には、入力波形を遅れ時間ｍずらした波形（計算式１ではｘ（ｋ−ｍ），計算式２ではｘ（ｉ−ｍ），計算式３ではｘ（ｊ−ｍ））は適宜の補間演算（例えば１次補間等）によるリサンプリングを行うものとする。
【００３９】
なお、上記図１０のステップＳ４５の低音域用計算式３においては、更に、計算時に入力波形ｘ（ｊ）の変数ｊを適宜間引いて計算しても差し支えなく、これによりより一層計算量を少なくすることができる。
【００４０】
また、上記図１０のステップＳ４１の自己相関の計算式１では、変数ｍは図１１のデータテーブル５０に示すｍ１〜ｍ４の４点について計算するものとしたが、変数ｍ１〜ｍ４の各数値は一例であって図示の例に限定されない。また、変数ｍの点数も４点に限らず更に細かく時間を設定して、計算点数を増やしても良い。また、自己相関の計算式１で使用する遅れ時間の時間範囲は、ステップＳ４３において大まかな音域判別が可能な計算結果を得ることができる時間範囲であればよく、これも上記の例に限らない。すなわち、自己相関の計算式１の要点は、所定の短い時間範囲を処理対象に絞った少ない計算量で入力波形について自己相関を求め、その結果から当該入力波形の音域を大雑把に判定できるものであればよい。また、上記図１０のステップＳ４１では、データテーブル５０を参照して変数ｍを設定するものとしたが、これに限らず、適宜の演算により変数ｍを設定してもよい。
【００４１】
また、上記実施例において、自己相関の計算式１〜３に使用する変数ｍの具体例として挙げた各数値は、調律の基準ピッチ（第４オクターブのＡ）を４４０Ｈｚ、調律器のサンプリングレートを４４．１ｋＨｚと仮定した場合の例である。従って、上記の条件がことなれば各数値はそれに応じて変動するものである。
【００４２】
なお、上記図１０のステップＳ４４及びＳ４５において、キー番号１〜４４を低音域に割り当て、キー番号４５〜８８を高音域に割り当てる例について説明したが、低音域及び高音域の音高割り当ての形態はこの限りではない。また、例えば、低音域にキー番号１〜５０を割り当て、高音域にキー番号３９〜８８を割り当てる等、低音域と高音域とで一部の音高が重複するよう割り当てるようにしても良い。このように構成することで、計算量が若干増えるものの誤検出の防止効果が高まる。なお、低音域及び高音域に対する音高の割り当ては、各音域の自己相関処理における変数ｍの数値範囲に対応する。
【００４３】
また、上記の実施例では、図１０のステップＳ４３においては、自己相関の計算式１の結果に基づき入力波形を高音域と低音域の何れか一方に判別する例を示したが、これに限らず、ステップＳ４３において、自己相関の計算式１の結果に基づき入力波形を高音域、低音域及び中音域の３領域のいずれかに判別するよう構成してもよい。この場合の音域判別構成としては、自己相関の計算式１の処理対象範囲中の最後の変数ｍ４におけるＲ（ｍ）の値が正か負かという判定要素に加えて、例えば、所定の極小値痾とＲ（ｍ）が最初に負になる時点との比較を判定要素とする。そして、両判定要素の結果の組み合わせに応じて、入力波形を高音域、中音域及び低音域の３領域で判別してもよい。この場合は、ステップＳ４３から中音域用の計算式に分岐する処理が更に加わる。中音域用の自己相関処理では、中音域に該当する音高の周波数の逆数を変数ｍとする自己相関の計算が行われる。
【００４４】
なお、上記実施例においては、この発明に係る調律器をピアノの調律に利用する例について述べたが、該調律器はピアノ以外の楽器の調律にも適用可能である。また、この発明に係る調律器をＰＤＡ（携帯情報端末）によって構成する例について説明したが、これに限らず、上述した調律機能を実現するソフトウェアプログラムを実行可能な装置でさえあれば、パーソナルコンピュータ、その他適宜の装置により構成及び実施して差し支えない。また、上述した調律機能を実現する信号処理回路を備えた調律器などで構成してもよい。また、表示器は、液晶画面に限らず、複数のＬＥＤを所定の時間範囲（位相範囲）に対応付けて列状に並べたもの等で構成してもよい。また、この発明は、コンピュータにおいて、上述した調律機能を実現するソフトウェアプログラムとして構成及び実施されてもよい。
【図面の簡単な説明】
【００４５】
【図１】この発明の一実施例に係る調律器を使用する際の全体像を示す外観図。
【図２】同実施例に係る調律器の電気的ハードウェア構成を示すブロック図。
【図３】同実施例に係る調律器の表示画面の一例を示す図であって、（ａ）は調律対象周波数と入力波形の周波数が一致している状態、（ｂ）は調律対象周波数と入力波形の周波数が不一致な状態。
【図４】同実施例に係る調律器における表示パターンの表示制御原理の概要を説明するための図であって、（ａ）は調律対象周波数と入力波形の周波数が一致している状態、（ｂ）は調律対象周波数と入力波形の周波数が不一致な状態。
【図５】図４に示す抽波形区間毎に抽出した入力波形を重ねたものを表す図。
【図６】同実施例に係る調律器における使用手順の概要を示すフローチャート。
【図７】同実施例に係る調律器における表示制御の手順を示すフローチャート。
【図８】同実施例に係る調律器において入力波形に基づき作成される２値化情報及び該２値化情報に基づき多値化情報を説明するための概念図であって、（ａ）は調律対象周波数と入力波形の周波数が一致している状態、（ｂ）は調律対象周波数と入力波形の周波数が不一致な状態。
【図９】同実施例に係る調律器における音名検出処理の手順を示すフローチャート。
【図１０】前記図９の音名検出処理における自己相関処理の手順を示すフローチャート。
【図１１】同実施例に係る自己相関の計算に用いる変数ｍのデータテーブル。
【図１２】（ａ）は同実施例に係る自己相関の計算式１の結果、（ｂ）は同実施例に係る自己相関の計算式２の結果、（２）は同実施例に係る自己相関の計算式３の結果をそれぞれ模擬的に示す図。
【符号の説明】
【００４６】
１  調律器、２  ピアノ、３  表示器、４  マイクロフォン、１０  ＣＰＵ、１１  ＲＯＭ、１２  ＲＡＭ、１３  インターフェース、１４  表示制御部、１５  操作検出部、３０  表示画面、３１  波形表示部、３２  表示パターン
【図１】
0004470823.tif 000005
【図２】
0004470823.tif 000006
【図３】
0004470823.tif 000007
【図４】
0004470823.tif 000008
【図５】
0004470823.tif 000009
【図６】
0004470823.tif 000010
【図７】
0004470823.tif 000011
【図８】
0004470823.tif 000012
【図９】
0004470823.tif 000013
【図１０】
0004470823.tif 000014
【図１１】
0004470823.tif 000015
【図１２】
0004470823.tif 000016