JP2871120B2 - 自動採譜装置 - Google Patents
自動採譜装置Info
- Publication number
- JP2871120B2 JP2871120B2 JP1143291A JP1143291A JP2871120B2 JP 2871120 B2 JP2871120 B2 JP 2871120B2 JP 1143291 A JP1143291 A JP 1143291A JP 1143291 A JP1143291 A JP 1143291A JP 2871120 B2 JP2871120 B2 JP 2871120B2
- Authority
- JP
- Japan
- Prior art keywords
- peak
- power spectrum
- fundamental frequency
- frequency
- peaks
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Description
【０００１】
【産業上の利用分野】本発明は、音楽を楽譜もしくは楽
譜に相当する符号に変換する自動採譜装置に係り、さら
に詳細には、実際に鳴っている音である基本周波数の抽
出のための構成に関する。
譜に相当する符号に変換する自動採譜装置に係り、さら
に詳細には、実際に鳴っている音である基本周波数の抽
出のための構成に関する。
【０００２】
【従来の技術】従来、複数の楽器によって演奏された音
楽の採譜は、音楽的知識を有する採譜者によって行われ
てきた。また、採譜を人が行うのでなく装置により自動
的に行うものでは、単音からなる音楽を採譜したり、あ
るいは鍵盤の押された情報から採譜を行う装置が提案さ
れている。
楽の採譜は、音楽的知識を有する採譜者によって行われ
てきた。また、採譜を人が行うのでなく装置により自動
的に行うものでは、単音からなる音楽を採譜したり、あ
るいは鍵盤の押された情報から採譜を行う装置が提案さ
れている。
【０００３】
【発明が解決しようとする課題】しかしながら、従来の
採譜装置では、楽器の種類や楽器数に制約があり、採譜
の範囲が限定され、一般的な音楽の演奏をそのまま採譜
できないという問題があった。また、複数の楽器によっ
て演奏された音楽をＡ／Ｄ変換した信号をＦＦＴ等によ
り周波数分析しただけでは、各楽器の倍音が多数観測さ
れ、どの音が基本周波数であるのか容易に判断できない
という問題があった。本発明は、上述した問題点を解決
するもので、音楽をＡ／Ｄ変換した信号を周波数分析
し、その結果から倍音を除去し、基本周波数のみを取り
出すことによって、楽器の種類や楽器数に制約を受けな
い自動採譜装置を提供することを目的とする。
採譜装置では、楽器の種類や楽器数に制約があり、採譜
の範囲が限定され、一般的な音楽の演奏をそのまま採譜
できないという問題があった。また、複数の楽器によっ
て演奏された音楽をＡ／Ｄ変換した信号をＦＦＴ等によ
り周波数分析しただけでは、各楽器の倍音が多数観測さ
れ、どの音が基本周波数であるのか容易に判断できない
という問題があった。本発明は、上述した問題点を解決
するもので、音楽をＡ／Ｄ変換した信号を周波数分析
し、その結果から倍音を除去し、基本周波数のみを取り
出すことによって、楽器の種類や楽器数に制約を受けな
い自動採譜装置を提供することを目的とする。
【０００４】
【課題を解決するための手段】上記目的を達成するため
に請求項１の発明は、音楽信号を楽譜もしくは楽譜に相
当する符号に変換する採譜装置であって、音楽信号を取
り込みＡ／Ｄ変換する信号取り込み部と、前記信号取り
込み部で取り込んだデータに対して、一定時間毎に周波
数解析を行うことにより、一定時間内での周波数方向の
パワー・スペクトルを計算する周波数解析処理部と、パ
ワー・スペクトルより基本周波数候補を抽出する基本周
波数候補抽出部とを備え、前記基本周波数候補抽出部
は、パワー・スペクトルのピークからパワーのしきい値
を計算するパワーのしきい値計算手段と、前記パワー・
スペクトルのピークのうちの一つが演奏音の基本周波数
かどうかを、それより高い周波数の他のピークのなかで
元のピークの倍音となっているピークと、元のピーク
と、両者の間にあるピークの周波数とパワーに基づいて
判定する基本周波数らしさ判定手段と、前記パワー・ス
ペクトルのピークのうちの一つがそれより低い周波数の
他のピークの中のどれかの倍音となっているかを検索す
る基本周波数検索手段と、前記基本周波数検索手段で検
索されたピークの倍音系列が、パワー・スペクトルのピ
ークの中にどの程度含まれているかどうかを算出する倍
音系列含有度計算手段と、前記パワー・スペクトルのピ
ークのうちの一つが他のピークの倍音であると判断され
た時に、倍音とされた方のピークをパワー・スペクトル
から除去する倍音除去手段とからなるものである。
に請求項１の発明は、音楽信号を楽譜もしくは楽譜に相
当する符号に変換する採譜装置であって、音楽信号を取
り込みＡ／Ｄ変換する信号取り込み部と、前記信号取り
込み部で取り込んだデータに対して、一定時間毎に周波
数解析を行うことにより、一定時間内での周波数方向の
パワー・スペクトルを計算する周波数解析処理部と、パ
ワー・スペクトルより基本周波数候補を抽出する基本周
波数候補抽出部とを備え、前記基本周波数候補抽出部
は、パワー・スペクトルのピークからパワーのしきい値
を計算するパワーのしきい値計算手段と、前記パワー・
スペクトルのピークのうちの一つが演奏音の基本周波数
かどうかを、それより高い周波数の他のピークのなかで
元のピークの倍音となっているピークと、元のピーク
と、両者の間にあるピークの周波数とパワーに基づいて
判定する基本周波数らしさ判定手段と、前記パワー・ス
ペクトルのピークのうちの一つがそれより低い周波数の
他のピークの中のどれかの倍音となっているかを検索す
る基本周波数検索手段と、前記基本周波数検索手段で検
索されたピークの倍音系列が、パワー・スペクトルのピ
ークの中にどの程度含まれているかどうかを算出する倍
音系列含有度計算手段と、前記パワー・スペクトルのピ
ークのうちの一つが他のピークの倍音であると判断され
た時に、倍音とされた方のピークをパワー・スペクトル
から除去する倍音除去手段とからなるものである。
【０００５】
【作用】請求項１の構成によれば、信号取り込み部は、
音楽をＡ／Ｄ変換し、計算機内で扱えるデジタル信号に
する。周波数解析処理部は、信号取り込み部で取り込ん
だデータに対して、ある一定時間毎に周波数解析を行う
ことにより、前記一定時間内での周波数方向のパワー・
スペクトルを計算する。基本周波数候補抽出部はこのパ
ワー・スペクトルより基本周波数候補を抽出する。そし
て、この基本周波数候補抽出部を構成する各手段は以下
のような作用をする。パワーのしきい値計算手段は、パ
ワー・スペクトルのピークからパワーのしきい値を計算
する。基本周波数らしさ判定手段は、パワー・スペクト
ルのピークのうちの一つが基本周波数かどうかを、それ
より高い周波数の他のピークのなかで元のピークの倍音
となっているピークと、元のピークと、両者の間にある
ピークの周波数とパワーとから判定する。基本周波数検
索手段は、前記パワー・スペクトルのピークのうちの一
つが、それより低い周波数の他のピークの中のどれかの
倍音となっているかを検索する。倍音系列含有度計算手
段は、前記基本周波数検索手段で検索されたピークの倍
音系列が、パワー・スペクトルのピークの中にどの程度
含まれているかどうかを算出する。倍音除去手段は、前
記パワー・スペクトルのピークのうちの一つが他のピー
クの倍音であると判断された場合に、倍音とされた方の
ピークをパワー・スペクトルから除去する。
音楽をＡ／Ｄ変換し、計算機内で扱えるデジタル信号に
する。周波数解析処理部は、信号取り込み部で取り込ん
だデータに対して、ある一定時間毎に周波数解析を行う
ことにより、前記一定時間内での周波数方向のパワー・
スペクトルを計算する。基本周波数候補抽出部はこのパ
ワー・スペクトルより基本周波数候補を抽出する。そし
て、この基本周波数候補抽出部を構成する各手段は以下
のような作用をする。パワーのしきい値計算手段は、パ
ワー・スペクトルのピークからパワーのしきい値を計算
する。基本周波数らしさ判定手段は、パワー・スペクト
ルのピークのうちの一つが基本周波数かどうかを、それ
より高い周波数の他のピークのなかで元のピークの倍音
となっているピークと、元のピークと、両者の間にある
ピークの周波数とパワーとから判定する。基本周波数検
索手段は、前記パワー・スペクトルのピークのうちの一
つが、それより低い周波数の他のピークの中のどれかの
倍音となっているかを検索する。倍音系列含有度計算手
段は、前記基本周波数検索手段で検索されたピークの倍
音系列が、パワー・スペクトルのピークの中にどの程度
含まれているかどうかを算出する。倍音除去手段は、前
記パワー・スペクトルのピークのうちの一つが他のピー
クの倍音であると判断された場合に、倍音とされた方の
ピークをパワー・スペクトルから除去する。
【０００６】
【実施例】以下、本発明を具体化した一実施例を図面を
参照して説明する。図１は本発明による自動採譜装置の
ブロック図である。本装置は、音楽信号が入力されるオ
ーディオ・アンプ１と、ローパス・フィルター２と、Ａ
／Ｄ変換装置３と、Ｉ／Ｏポート４と、基本周波数候補
を抽出処理するＣＰＵ５と、ＲＡＭ６と、ＲＯＭ７と、
抽出した結果を表示するデイスプレイ８とから構成され
ている。
参照して説明する。図１は本発明による自動採譜装置の
ブロック図である。本装置は、音楽信号が入力されるオ
ーディオ・アンプ１と、ローパス・フィルター２と、Ａ
／Ｄ変換装置３と、Ｉ／Ｏポート４と、基本周波数候補
を抽出処理するＣＰＵ５と、ＲＡＭ６と、ＲＯＭ７と、
抽出した結果を表示するデイスプレイ８とから構成され
ている。
【０００７】次に、上記図１に図２、図３を加えて、本
自動採譜装置により行われる基本周波数抽出処理のため
の機能構成を説明する。本装置は、機能構成要素とし
て、信号取り込み部２１と高速フーリエ変換（ＦＦＴ）
処理部２２と基本周波数候補抽出部２３とからなる。信
号取り込み部２１は、オーディオ・アンプ１、ローパス
・フィルタ２、Ａ／Ｄ変換装置３、Ｉ／Ｏポート４、Ｃ
ＰＵ５等により構成される。入力された音楽信号はオー
ディオ・アンプ１により増幅され、この増幅された信号
は、ローパス・フィルタ２に入力され、例えば、５．５
ｋＨｚ以下の周波数成分のみが通過し、標本化時の折返
し歪みを抑えている。この出力信号は、Ａ／Ｄ変換装置
３により１２ｋＨｚ、１６ｂｉｔで標本化される。標本
化されたデータは、Ｉ／Ｏポート４を介し、ＣＰＵ５に
取り込まれ、ＲＡＭ６に記憶される。
自動採譜装置により行われる基本周波数抽出処理のため
の機能構成を説明する。本装置は、機能構成要素とし
て、信号取り込み部２１と高速フーリエ変換（ＦＦＴ）
処理部２２と基本周波数候補抽出部２３とからなる。信
号取り込み部２１は、オーディオ・アンプ１、ローパス
・フィルタ２、Ａ／Ｄ変換装置３、Ｉ／Ｏポート４、Ｃ
ＰＵ５等により構成される。入力された音楽信号はオー
ディオ・アンプ１により増幅され、この増幅された信号
は、ローパス・フィルタ２に入力され、例えば、５．５
ｋＨｚ以下の周波数成分のみが通過し、標本化時の折返
し歪みを抑えている。この出力信号は、Ａ／Ｄ変換装置
３により１２ｋＨｚ、１６ｂｉｔで標本化される。標本
化されたデータは、Ｉ／Ｏポート４を介し、ＣＰＵ５に
取り込まれ、ＲＡＭ６に記憶される。
【０００８】ＦＦＴ処理部２２は、ＣＰＵ５、ＲＡＭ６
等でなり、ＣＰＵ５はＲＡＭ６より標本化されたデータ
を読み出し、２５ｍｓｅｃ毎を１フレームとし、１フレ
ーム毎に８５．３ｍｓｅｃハミング窓を掛けた後、ＦＦ
Ｔ分析により対数パワー・スペクトルが算出される。次
に、ＣＰＵ５は、算出された対数パワー・スペクトルか
ら放射線内挿処理によりピーク周波数を求める。図５は
以上のようにして求めたピーク・スペクトルを、時間軸
を横軸に、ピーク周波数を鍵盤番号に変換して縦軸にと
り、強度を濃淡で示したものである。同図は倍音除去前
のデータである。基本周波数候補抽出部２３は、ＣＰＵ
５でなり、図３のようにパワーのしきい値計算手段３１
と、基本周波数らしさ判定手段３２と、基本周波数検索
手段３３と、倍音系列含有度計算手段３４と、倍音除去
手段３５とから構成されている。
等でなり、ＣＰＵ５はＲＡＭ６より標本化されたデータ
を読み出し、２５ｍｓｅｃ毎を１フレームとし、１フレ
ーム毎に８５．３ｍｓｅｃハミング窓を掛けた後、ＦＦ
Ｔ分析により対数パワー・スペクトルが算出される。次
に、ＣＰＵ５は、算出された対数パワー・スペクトルか
ら放射線内挿処理によりピーク周波数を求める。図５は
以上のようにして求めたピーク・スペクトルを、時間軸
を横軸に、ピーク周波数を鍵盤番号に変換して縦軸にと
り、強度を濃淡で示したものである。同図は倍音除去前
のデータである。基本周波数候補抽出部２３は、ＣＰＵ
５でなり、図３のようにパワーのしきい値計算手段３１
と、基本周波数らしさ判定手段３２と、基本周波数検索
手段３３と、倍音系列含有度計算手段３４と、倍音除去
手段３５とから構成されている。
【０００９】基本周波数候補抽出部２３では、ＦＦＴ処
理部２２の各フレーム毎に処理を進める。いま、あるフ
レームに対して、ＦＦＴ処理部２２で求めたスペクトル
のピークがＮ個あるとし、以下、図４のフローチャート
に従って説明する。図４は、基本周波数候補抽出部２３
での処理を、ある１フレーム分について示したものであ
る。Ｎ個のピークについては、各々周波数Ｆとパワーの
値Ｐが周波数によってソートされてＲＡＭ６に記憶され
ている。このＮ個のピークのパワーの値から、パワーに
よるしきい値Ｓｐを計算する（ステップ＃１）。Ｎ個の
ピークについては、周波数の高い方から低い方へ、順に
処理を進める。Ｎ個のピークの、周波数の高い方からｎ
番目のピークの周波数をＦ（ｎ）、パワーをＰ（ｎ）と
する。Ｆ（ｎ）が基本周波数であるかどうかを、Ｆ
（ｎ）よりも高い周波数Ｆ（ｍ）およびパワーＰ（ｍ）
（１≦ｍ＜ｎ）から判定する（＃２，＃３）。この判定
は以下のようにして行う。まず、Ｆ（ｍ）の中にＦ
（ｎ）の第Ｘ次倍音（Ｘ≧２）となっているＦ（ｘ）が
あるかどうか探す。Ｆ（ｘ）が見つかった場合、Ｆ
（ｎ）とＦ（ｘ）の関係は、Ｆ（ｎ）が基本周波数で、
Ｆ（ｘ）が第Ｘ次倍音となっている場合（これをケース
１とする）と、Ｆ（ｎ）、Ｆ（ｘ）がより低い、両者の
公約数となっている周波数の、それぞれれＬ次（Ｌ≧
２）、Ｌ×Ｘ次倍音となっている場合（これをケース２
とする）とがあり得る。
理部２２の各フレーム毎に処理を進める。いま、あるフ
レームに対して、ＦＦＴ処理部２２で求めたスペクトル
のピークがＮ個あるとし、以下、図４のフローチャート
に従って説明する。図４は、基本周波数候補抽出部２３
での処理を、ある１フレーム分について示したものであ
る。Ｎ個のピークについては、各々周波数Ｆとパワーの
値Ｐが周波数によってソートされてＲＡＭ６に記憶され
ている。このＮ個のピークのパワーの値から、パワーに
よるしきい値Ｓｐを計算する（ステップ＃１）。Ｎ個の
ピークについては、周波数の高い方から低い方へ、順に
処理を進める。Ｎ個のピークの、周波数の高い方からｎ
番目のピークの周波数をＦ（ｎ）、パワーをＰ（ｎ）と
する。Ｆ（ｎ）が基本周波数であるかどうかを、Ｆ
（ｎ）よりも高い周波数Ｆ（ｍ）およびパワーＰ（ｍ）
（１≦ｍ＜ｎ）から判定する（＃２，＃３）。この判定
は以下のようにして行う。まず、Ｆ（ｍ）の中にＦ
（ｎ）の第Ｘ次倍音（Ｘ≧２）となっているＦ（ｘ）が
あるかどうか探す。Ｆ（ｘ）が見つかった場合、Ｆ
（ｎ）とＦ（ｘ）の関係は、Ｆ（ｎ）が基本周波数で、
Ｆ（ｘ）が第Ｘ次倍音となっている場合（これをケース
１とする）と、Ｆ（ｎ）、Ｆ（ｘ）がより低い、両者の
公約数となっている周波数の、それぞれれＬ次（Ｌ≧
２）、Ｌ×Ｘ次倍音となっている場合（これをケース２
とする）とがあり得る。
【００１０】このどちらの場合であるかを判断するため
に、Ｆ（ｘ）とＦ（ｎ）の間に存在する周波数Ｆ
（ｙ）、およびパワーＰ（ｙ）（ｘ＜ｙ＜ｎ）につい
て、以下のような処理をする。Ｐ（ｙ）＞（Ｐ（ｘ）＋
Ｐ（ｎ））／２であれば、Ｆ（ｙ）とＦ（ｘ）、Ｆ
（ｎ）の周波数差を計算する。いま、Ｆｘｙ＝Ｆ（ｘ）
−Ｆ（ｙ）、Ｆｙｎ＝Ｆ（ｙ）−Ｆ（ｎ）とする。Ｆｘ
ｙとＦｙｎの小さい方をＦ０とする。Ｆ０がＦ（ｎ）、
Ｆ（ｘ）を両方とも倍音とする場合は、前述のケース２
である可能性が高いので、Ｆ（ｎ）は基本周波数でない
と判断する。また、Ｆ（ｘ）が見つからなかった場合も
Ｆ（ｎ）は基本周波数でないと判断する。これ以外の場
合、すなわちＦ（ｍ）の中に存在するＦ（ｎ）の第Ｘ次
倍音（Ｘ≧２）となっている全てのＦ（ｘ）について、
Ｆ（ｘ）とＦ（ｎ）の間に存在するＦ（ｙ）から計算さ
れた全てのＦ０の中に、Ｆ（ｎ）、Ｆ（ｘ）を両方とも
倍音とするものがひとつもない場合、Ｆ（ｎ）は基本周
波数であると判断される。
に、Ｆ（ｘ）とＦ（ｎ）の間に存在する周波数Ｆ
（ｙ）、およびパワーＰ（ｙ）（ｘ＜ｙ＜ｎ）につい
て、以下のような処理をする。Ｐ（ｙ）＞（Ｐ（ｘ）＋
Ｐ（ｎ））／２であれば、Ｆ（ｙ）とＦ（ｘ）、Ｆ
（ｎ）の周波数差を計算する。いま、Ｆｘｙ＝Ｆ（ｘ）
−Ｆ（ｙ）、Ｆｙｎ＝Ｆ（ｙ）−Ｆ（ｎ）とする。Ｆｘ
ｙとＦｙｎの小さい方をＦ０とする。Ｆ０がＦ（ｎ）、
Ｆ（ｘ）を両方とも倍音とする場合は、前述のケース２
である可能性が高いので、Ｆ（ｎ）は基本周波数でない
と判断する。また、Ｆ（ｘ）が見つからなかった場合も
Ｆ（ｎ）は基本周波数でないと判断する。これ以外の場
合、すなわちＦ（ｍ）の中に存在するＦ（ｎ）の第Ｘ次
倍音（Ｘ≧２）となっている全てのＦ（ｘ）について、
Ｆ（ｘ）とＦ（ｎ）の間に存在するＦ（ｙ）から計算さ
れた全てのＦ０の中に、Ｆ（ｎ）、Ｆ（ｘ）を両方とも
倍音とするものがひとつもない場合、Ｆ（ｎ）は基本周
波数であると判断される。
【００１１】以上により、Ｐ（ｎ）＞ＳｐかつＦ（ｎ）
が基本周波数であると判断された場合は（＃４，＃５で
ＹＥＳ）、そのピークは基本周波数であるとして次のピ
ークに処理を移す（＃６，＃７）。そうでない場合は、
Ｆ（ｎ）よりも低い周波数Ｆ（ｐ）（ｎ＜ｐ≦Ｎ）の中
に、Ｆ（ｎ）を第Ｚ次倍音（Ｚ≧２）とするＦ（ｚ）が
あるか探す（＃８，＃９）。Ｆ（ｚ）が見つかった場合
には、Ｆ（ｚ）のＺ次までの倍音系列をＦ（ｑ）（ｚ≦
ｑ≦ｎ）の中で探す。Ｆ（ｑ）の中にＺ個含まれるはず
のＦ（ｚ）の倍音系列のうち、実際にＦ（ｑ）の中で見
つかった個数をＷ個とし、倍音系列含有度Ｗ／Ｚを計算
する（＃１０）。Ｗ／Ｚがあるしきい値以上の場合には
（＃１１でＹＥＳ）、Ｆ（ｎ）はＦ（ｚ）の倍音である
として除去する（＃１２）。条件を満たすＦ（ｚ）が見
つからなかった場合には、Ｆ（ｎ）を除去しない。な
お、Ｆ（ｎ）を除去した場合でも以降の処理に支障をき
たすことのないように、Ｆ（ｎ）にピークが存在したと
いう情報は残しておく。以上の処理をＮ個のピーク全て
について行ったら、次のフレームに処理を移す。ＦＦＴ
処理部２２でパワー・スペクトルを算出した全てのフレ
ームについて以上の処理を終了した時、基本周波数候補
抽出部２３の処理が終了する。
が基本周波数であると判断された場合は（＃４，＃５で
ＹＥＳ）、そのピークは基本周波数であるとして次のピ
ークに処理を移す（＃６，＃７）。そうでない場合は、
Ｆ（ｎ）よりも低い周波数Ｆ（ｐ）（ｎ＜ｐ≦Ｎ）の中
に、Ｆ（ｎ）を第Ｚ次倍音（Ｚ≧２）とするＦ（ｚ）が
あるか探す（＃８，＃９）。Ｆ（ｚ）が見つかった場合
には、Ｆ（ｚ）のＺ次までの倍音系列をＦ（ｑ）（ｚ≦
ｑ≦ｎ）の中で探す。Ｆ（ｑ）の中にＺ個含まれるはず
のＦ（ｚ）の倍音系列のうち、実際にＦ（ｑ）の中で見
つかった個数をＷ個とし、倍音系列含有度Ｗ／Ｚを計算
する（＃１０）。Ｗ／Ｚがあるしきい値以上の場合には
（＃１１でＹＥＳ）、Ｆ（ｎ）はＦ（ｚ）の倍音である
として除去する（＃１２）。条件を満たすＦ（ｚ）が見
つからなかった場合には、Ｆ（ｎ）を除去しない。な
お、Ｆ（ｎ）を除去した場合でも以降の処理に支障をき
たすことのないように、Ｆ（ｎ）にピークが存在したと
いう情報は残しておく。以上の処理をＮ個のピーク全て
について行ったら、次のフレームに処理を移す。ＦＦＴ
処理部２２でパワー・スペクトルを算出した全てのフレ
ームについて以上の処理を終了した時、基本周波数候補
抽出部２３の処理が終了する。
【００１２】図６は、図５のデータに基本周波数候補抽
出部２３による倍音除去の処理を施した後のデータであ
る。図６では、図５に比べて、倍音除去により不要なデ
ータが除去されていることが分かる。本発明は上記の実
施例に限定するものではなく、その趣旨を逸脱しない範
囲において種々の変更を加えることができる。例えば、
本実施例においては、ＦＦＴ処理部２２のフレーム毎に
基本周波数候補抽出部２３の処理を行ったが、予め近接
するフレーム間でパワー・スペクトルの安定部が続く区
間を一つの分析区間として設定し、その分析区間毎に以
降の処理を行うことも可能である。
出部２３による倍音除去の処理を施した後のデータであ
る。図６では、図５に比べて、倍音除去により不要なデ
ータが除去されていることが分かる。本発明は上記の実
施例に限定するものではなく、その趣旨を逸脱しない範
囲において種々の変更を加えることができる。例えば、
本実施例においては、ＦＦＴ処理部２２のフレーム毎に
基本周波数候補抽出部２３の処理を行ったが、予め近接
するフレーム間でパワー・スペクトルの安定部が続く区
間を一つの分析区間として設定し、その分析区間毎に以
降の処理を行うことも可能である。
【００１３】
【発明の効果】以上のように本発明によれば、Ａ／Ｄ変
換して取り込んだ音楽信号を周波数解析し、その結果か
ら倍音を除去することによって、実際に鳴っている音の
基本周波数のみを取り出すようにしているので、楽器の
種類や楽器数に制約を受けずに基本周波数を抽出するこ
とができ、採譜することが容易となる。
換して取り込んだ音楽信号を周波数解析し、その結果か
ら倍音を除去することによって、実際に鳴っている音の
基本周波数のみを取り出すようにしているので、楽器の
種類や楽器数に制約を受けずに基本周波数を抽出するこ
とができ、採譜することが容易となる。
【図１】 本発明の一実施例による自動採譜装置のブロ
ック図である。
ック図である。
【図２】 自動採譜装置の基本周波数抽出までの機能構
成図である。
成図である。
【図３】 基本周波数候補抽出部２３の構成図である。
【図４】 基本周波数候補抽出部２３における処理を示
すフローチャートである。
すフローチャートである。
【図５】 ＦＦＴ処理部２２により算出された、倍音除
去処理前のパワー・スペクトルの説明図である。
去処理前のパワー・スペクトルの説明図である。
【図６】 図５に基本周波数候補抽出部２による倍音除
去処理を施した後のパワー・スペクトルの説明図であ
る。
去処理を施した後のパワー・スペクトルの説明図であ
る。
３ Ａ／Ｄ変換装置 ５ ＣＰＵ ６ ＲＡＭ ２１ 信号取り込み部 ２２ ＦＦＴ処理部 ２３ 基本周波数候補抽出部 ３１ パワーのしきい値計算手段 ３２ 基本周波数らしさ判定手段 ３３ 基本周波数検索手段 ３４ 倍音系列含有度計算手段 ３５ 倍音除去手段
Claims (1)
- 【請求項１】 音楽信号を楽譜もしくは楽譜に相当する
符号に変換する採譜装置において、 音楽信号を取り込みＡ／Ｄ変換する信号取り込み部と、 前記信号取り込み部で取り込んだデータに対して、一定
時間毎に周波数解析を行うことにより、一定時間内での
周波数方向のパワー・スペクトルを計算する周波数解析
処理部と、 パワー・スペクトルより基本周波数候補を抽出する基本
周波数候補抽出部とを備え、 前記基本周波数候補抽出部は、パワー・スペクトルのピ
ークからパワーのしきい値を計算するパワーのしきい値
計算手段と、前記パワー・スペクトルのピークのうちの
一つが演奏音の基本周波数かどうかを、それより高い周
波数の他のピークのなかで元のピークの倍音となってい
るピークと、元のピークと、両者の間にあるピークの周
波数とパワーに基づいて判定する基本周波数らしさ判定
手段と、前記パワー・スペクトルのピークのうちの一つ
がそれより低い周波数の他のピークの中のどれかの倍音
となっているかを検索する基本周波数検索手段と、前記
基本周波数検索手段で検索されたピークの倍音系列が、
パワー・スペクトルのピークの中にどの程度含まれてい
るかどうかを算出する倍音系列含有度計算手段と、前記
パワー・スペクトルのピークのうちの一つが他のピーク
の倍音であると判断された時に、倍音とされた方のピー
クをパワー・スペクトルから除去する倍音除去手段とか
らなる ことを特徴とする自動採譜装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1143291A JP2871120B2 (ja) | 1991-01-07 | 1991-01-07 | 自動採譜装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1143291A JP2871120B2 (ja) | 1991-01-07 | 1991-01-07 | 自動採譜装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04261591A JPH04261591A (ja) | 1992-09-17 |
| JP2871120B2 true JP2871120B2 (ja) | 1999-03-17 |
Family
ID=11777922
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1143291A Expired - Fee Related JP2871120B2 (ja) | 1991-01-07 | 1991-01-07 | 自動採譜装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2871120B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4697919B2 (ja) * | 2001-06-29 | 2011-06-08 | 大日本印刷株式会社 | 音響信号の符号化方法 |
| JP4672474B2 (ja) * | 2005-07-22 | 2011-04-20 | 株式会社河合楽器製作所 | 自動採譜装置及びプログラム |
| JP4953068B2 (ja) * | 2007-02-26 | 2012-06-13 | 独立行政法人産業技術総合研究所 | 和音判別装置、和音判別方法およびプログラム |
| JP6307814B2 (ja) * | 2013-08-26 | 2018-04-11 | カシオ計算機株式会社 | 基音可視化装置、基音可視化方法およびプログラム |
| JP6904198B2 (ja) * | 2017-09-25 | 2021-07-14 | 富士通株式会社 | 音声処理プログラム、音声処理方法および音声処理装置 |
| CN111542874B (zh) | 2017-11-07 | 2023-09-01 | 雅马哈株式会社 | 数据生成装置及记录介质 |
-
1991
- 1991-01-07 JP JP1143291A patent/JP2871120B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH04261591A (ja) | 1992-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109147796B (zh) | 语音识别方法、装置、计算机设备及计算机可读存储介质 | |
| JP3114975B2 (ja) | 音素推定を用いた音声認識回路 | |
| US8543387B2 (en) | Estimating pitch by modeling audio as a weighted mixture of tone models for harmonic structures | |
| CN113012720B (zh) | 谱减法降噪下多语音特征融合的抑郁症检测方法 | |
| RU2427909C2 (ru) | Способ формирования отпечатка для звукового сигнала | |
| CN111583957B (zh) | 基于五音阶乐律声谱图和级联神经网络的戏曲分类方法 | |
| CN110248300A (zh) | 一种基于自主学习的啸叫抑制方法及扩声系统 | |
| JP2871120B2 (ja) | 自動採譜装置 | |
| CN110379438B (zh) | 一种语音信号基频检测与提取方法及系统 | |
| JP2779325B2 (ja) | ボコーダーにおける前処理の相関関係式を用いたピッチ検索時間短縮方法 | |
| CN110415730B (zh) | 一种音乐分析数据集构建方法及基于其的音高、时值提取方法 | |
| JP5924968B2 (ja) | 楽譜位置推定装置、及び楽譜位置推定方法 | |
| KR100766170B1 (ko) | 다중 레벨 양자화를 이용한 음악 요약 장치 및 방법 | |
| CN113066512B (zh) | 佛教音乐识别方法、装置、设备及存储介质 | |
| CN114996489A (zh) | 新闻数据的违规检测方法、装置、设备及存储介质 | |
| CN114420099A (zh) | 一种多声源融合场景的人声检测方法及装置 | |
| JP2806048B2 (ja) | 自動採譜装置 | |
| JP3232112B2 (ja) | 小節境界時刻抽出装置 | |
| JP2004102023A (ja) | 特定音響信号検出方法、信号検出装置、信号検出プログラム及び記録媒体 | |
| JP3102089B2 (ja) | 自動採譜装置 | |
| CN113744715A (zh) | 声码器语音合成方法、装置、计算机设备及存储介质 | |
| JP3919359B2 (ja) | 楽音信号のアタック位置検出装置 | |
| JP3012994B2 (ja) | 音韻識別方法 | |
| JP3423233B2 (ja) | オーディオ信号処理方法および装置 | |
| JPH04323696A (ja) | 自動採譜装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19981208 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080108 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090108 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090108 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100108 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |