JP2003058149A

JP2003058149A - 周波数解析方法および音響信号の符号化方法

Info

Publication number: JP2003058149A
Application number: JP2001243179A
Authority: JP
Inventors: Toshio Motegi; 敏雄茂出木
Original assignee: Dai Nippon Printing Co Ltd
Current assignee: Dai Nippon Printing Co Ltd
Priority date: 2001-08-10
Filing date: 2001-08-10
Publication date: 2003-02-28
Anticipated expiration: 2021-08-10
Also published as: JP4662407B2

Abstract

(57)【要約】【課題】標準周波数間の境界付近に信号成分が存在す
る場合であっても、正確にその信号成分を抽出すること
ができ、さらには、精度の高い符号化を行うことが可能
な周波数解析方法および音響信号の符号化方法を提供す
る。【解決手段】音響信号に対して単位区間を設定して区
間信号の抽出を行い（Ｓ１）、区間信号に対して短時間
フーリエ変換を用いて、微細な間隔で設定された各周波
数に対応する調和信号との相関演算を行う（Ｓ２）。相
関演算の結果、ピークとなる微細周波数を選出し（Ｓ
３）、選出されたピーク周波数に対応する調和信号に限
って区間信号との相関を、一般化調和解析の手法により
求め、求めた相関値、および対応するピーク周波数に最
も近い標準周波数の組を周波数成分として抽出する（Ｓ
４）。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、放送メディア（ラジ
オ、テレビ）、通信メディア（ＣＳ映像・音声配信、イ
ンターネット音楽配信、通信カラオケ）、パッケージメ
ディア（ＣＤ、ＭＤ、カセット、ビデオ、ＬＤ、ＣＤ−
ＲＯＭ、ゲームカセット、携帯音楽プレーヤ向け固体メ
モリ媒体）などで提供する各種オーディオコンテンツの
制作、並びに、音楽演奏録音信号から楽譜出版、通信カ
ラオケ配信用ＭＩＤＩデータ、演奏ガイド機能付き電子
楽器向け自動演奏データ、携帯電話・ＰＨＳ・ポケベル
などの着信メロディデータを自動的に作成する自動採譜
技術に関する。

【０００２】

【従来の技術】音響信号に代表される時系列信号には、
その構成要素として複数の周期信号が含まれている。こ
のため、与えられた時系列信号にどのような周期信号が
含まれているかを解析する手法は、古くから知られてい
る。例えば、フーリエ解析は、与えられた時系列信号に
含まれる周波数成分を解析するための方法として広く利
用されている。

【０００３】このような時系列信号の解析方法を利用す
れば、音響信号を符号化することも可能である。コンピ
ュータの普及により、原音となるアナログ音響信号を所
定のサンプリング周波数でサンプリングし、各サンプリ
ング時の信号強度を量子化してデジタルデータとして取
り込むことが容易にできるようになってきており、こう
して取り込んだデジタルデータに対してフーリエ解析な
どの手法を適用し、原音信号に含まれていた周波数成分
を抽出すれば、各周波数成分を示す符号によって原音信
号の符号化が可能になる。

【０００４】また、電子楽器による楽器音を符号化しよ
うという発想から生まれたＭＩＤＩ（Musical Instrume
nt Digital Interface）規格も、パーソナルコンピュー
タの普及とともに盛んに利用されるようになってきてい
る。このＭＩＤＩ規格による符号データ（以下、ＭＩＤ
Ｉデータという）は、基本的には、楽器のどの鍵盤キー
を、どの程度の強さで弾いたか、という楽器演奏の操作
を記述したデータであり、このＭＩＤＩデータ自身に
は、実際の音の波形は含まれていない。そのため、実際
の音を再生する場合には、楽器音の波形を記憶したＭＩ
ＤＩ音源が別途必要になるが、その符号化効率の高さが
注目を集めており、ＭＩＤＩ規格による符号化および復
号化の技術は、現在、パーソナルコンピュータを用いて
楽器演奏、楽器練習、作曲などを行うソフトウェアに広
く採り入れられている。

【０００５】そこで、音響信号に代表される時系列信号
に対して、所定の手法で解析を行うことにより、その構
成要素となる周期信号を抽出し、抽出した周期信号をＭ
ＩＤＩデータを用いて符号化しようとする提案がなされ
ている。例えば、特開平１０−２４７０９９号公報、特
開平１１−７３１９９号公報、特開平１１−７３２００
号公報、特開平１１−９５７５３号公報、特開２０００
−９９００９号公報、特開２０００−９９０９２号公
報、特開２０００−９９０９３号公報、特開２０００−
２６１３２２号公報、特開２００１−５４５０号公報、
特開２００１−１４８６３３号公報には、任意の時系列
信号について、構成要素となる周波数を解析し、その解
析結果からＭＩＤＩデータを作成することができる種々
の方法が提案されている。

【０００６】上記公報に記載された発明では、音響信号
との相関を求める際に、音階で特定される標準周波数よ
り細かい間隔で微細周波数を設定し、この微細周波数に
対応した周期関数を用いて行っている。音響信号として
再生可能とするためには、標準周波数に対応して符号化
する必要があるので、具体的には、相関値が高い微細周
波数を求め、その相関値をその微細周波数に最も近い標
準周波数の値として符号化している。この際、相関値が
高い微細周波数が隣接する標準周波数間の境界付近に存
在する場合には、他の信号成分による影響を受けるもの
と判断し、符号化対象として抽出しない処理を行ってい
る。

【０００７】

【発明が解決しようとする課題】しかしながら、上記の
ように隣接する標準周波数間の境界付近に相関値の大き
い微細周波数が存在した場合に、その微細周波数を抽出
しないようにすると、実際に境界部分に信号成分が存在
していたとしても削除してしまうことになり、逆にこの
微細周波数を抽出するようにすると、他の信号成分によ
る影響であって実際には存在していなくても抽出してし
まうことになるという問題が生じている。

【０００８】上記のような点に鑑み、本発明は、標準周
波数間の境界付近に信号成分が存在する場合であって
も、正確にその信号成分を抽出することができ、さらに
は、精度の高い符号化を行うことが可能な周波数解析方
法および音響信号の符号化方法を提供することを課題と
する。

【０００９】

【課題を解決するための手段】上記課題を解決するた
め、本発明では、時系列信号から所定の標準周波数に対
応した信号成分を分離するための周波数解析方法とし
て、前記標準周波数の間隔よりも微細な間隔で複数の微
細周波数を設定し、各微細周波数に対応する複数の調和
信号集合を準備する調和信号準備段階と、前記各標準周
波数に対応する相関値を格納するための主相関配列と、
各微細周波数に対応する相関値を格納するための微細相
関配列を準備する配列準備段階と、前記時系列信号に所
定の長さの単位区間を設定し、区間信号を抽出するため
の区間信号抽出段階と、前記全ての調和信号と前記区間
信号との相関を演算して、各微細周波数に対応する相関
値を算出し、前記微細相関配列に算出された相関値を設
定する相関演算段階と、前記微細相関配列中の相関配列
より周囲の微細周波数における相関値に比べ局所的にピ
ークとなる微細周波数を複数抽出するピーク周波数選出
段階と、前記各ピークとなる微細周波数に対応する調和
信号と前記区間信号との相関を再度演算し、前記微細周
波数の最近傍の標準周波数に対応する主相関配列に設定
し、前記再度算出された相関値と調和信号との積で生成
される含有信号を前記区間信号から減じ、差分信号を新
たに区間信号とすることにより区間信号を更新する処理
を抽出された全てのピークに対して順次行って主相関配
列の値を決定するスペクトル算出段階を有することを特
徴とする。

【００１０】本発明によれば、標準周波数より微細な間
隔で微細周波数を設定し、各微細周波数の区間信号に対
する相関を求め、相関値が局所的にピークとなっている
微細周波数（以下、ピーク周波数という）を複数選出
し、選出された複数のピーク周波数と区間信号との相関
を一般化調和解析の手法により求め、算出した相関値を
対応するピーク周波数に最も近い標準周波数の相関値と
して、抽出すべき周波数成分を決定するようにしたの
で、隣接する標準周波数間の境界付近に位置する微細周
波数がピークとなる場合であっても、正確に周波数成分
を抽出することが可能となる。

【００１１】

【発明の実施の形態】以下、本発明の実施形態について
図面を参照して詳細に説明する。

【００１２】（1.基本原理）はじめに、本発明に係る周
波数解析方法および音響信号符号化方法の基本原理を述
べておく。この基本原理は、前掲の各公報に開示されて
いるので、ここではその概要のみを簡単に述べることに
する。

【００１３】図１（ａ）に示すように、時系列信号とし
てアナログ音響信号が与えられたものとする。図１の例
では、横軸に時間ｔ、縦軸に振幅（強度）をとって、こ
の音響信号を示している。ここでは、まずこのアナログ
音響信号を、デジタルの音響データとして取り込む処理
を行う。これは、従来の一般的なＰＣＭの手法を用い、
所定のサンプリング周波数でこのアナログ音響信号をサ
ンプリングし、振幅を所定の量子化ビット数を用いてデ
ジタルデータに変換する処理を行えば良い。ここでは、
説明の便宜上、ＰＣＭの手法でデジタル化した音響デー
タの波形も図１（ａ）のアナログ音響信号と同一の波形
で示すことにする。

【００１４】続いて、この解析対象となる音響信号の時
間軸上に、複数の単位区間を設定する。図１（ａ）に示
す例では、時間軸ｔ上に等間隔に６つの時刻ｔ１〜ｔ６
が定義され、これら各時刻を始点および終点とする５つ
の単位区間ｄ１〜ｄ５が設定されている。図１の例で
は、全て同一の区間長をもった単位区間が時間軸上で重
複せずに設定されているが、隣接する単位区間が時間軸
上で部分的に重なり合うような区間設定を行ってもかま
わない。

【００１５】こうして単位区間が設定されたら、各単位
区間ごとの音響信号（以下、区間信号と呼ぶことにす
る）について、それぞれ代表周波数を選出する。各区間
信号には、通常、様々な周波数成分が含まれているが、
例えば、その中で成分の強度割合の大きな周波数成分を
代表周波数として選出すれば良い。ここで、代表周波数
とはいわゆる基本周波数が一般的であるが、音声のフォ
ルマント周波数などの倍音周波数や、ノイズ音源のピー
ク周波数も代表周波数として扱うことがある。代表周波
数は１つだけ選出しても良いが、音響信号によっては複
数の代表周波数を選出した方が、より精度の高い符号化
が可能になる。図１（ｂ）には、個々の単位区間ごとに
それぞれ３つの代表周波数を選出し、１つの代表周波数
を１つの代表符号（図では便宜上、音符として示してあ
る）として符号化した例が示されている。ここでは、代
表符号（音符）を収容するために３つのトラックＴ１，
Ｔ２，Ｔ３が設けられているが、これは個々の単位区間
ごとに選出された３つずつの代表符号を、それぞれ異な
るトラックに収容するためである。

【００１６】例えば、単位区間ｄ１について選出された
代表符号ｎ（ｄ１，１），ｎ（ｄ１，２），ｎ（ｄ１，
３）は、それぞれトラックＴ１，Ｔ２，Ｔ３に収容され
ている。ここで、各符号ｎ（ｄ１，１），ｎ（ｄ１，
２），ｎ（ｄ１，３）は、ＭＩＤＩ符号におけるノート
ナンバーを示す符号である。ＭＩＤＩ符号におけるノー
トナンバーは、０〜１２７までの１２８通りの値をと
り、それぞれピアノの鍵盤の１つのキーを示すことにな
る。具体的には、例えば、代表周波数として４４０Ｈｚ
が選出された場合、この周波数はノートナンバーｎ＝６
９（ピアノの鍵盤中央の「ラ音（Ａ３音）」に対応）に
相当するので、代表符号としては、ｎ＝６９が選出され
ることになる。もっとも、図１（ｂ）は、上述の方法に
よって得られる代表符号を音符の形式で示した概念図で
あり、実際には、各音符にはそれぞれ強度に関するデー
タも付加されている。例えば、トラックＴ１には、ノー
トナンバーｎ（ｄ１，１），ｎ（ｄ２，１）・・・なる
音高を示すデータとともに、ｅ（ｄ１，１），ｅ（ｄ
２，１）・・・なる強度を示すデータが収容されること
になる。この強度を示すデータは、各代表周波数の成分
が、元の区間信号にどの程度の度合いで含まれていたか
によって決定される。具体的には、各代表周波数をもっ
た周期関数の区間信号に対する相関値に基づいて強度を
示すデータが決定されることになる。また、図１（ｂ）
に示す概念図では、音符の横方向の位置によって、個々
の単位区間の時間軸上での位置が示されているが、実際
には、この時間軸上での位置を正確に数値として示すデ
ータが各音符に付加されていることになる。

【００１７】音響信号を符号化する形式としては、必ず
しもＭＩＤＩ形式を採用する必要はないが、この種の符
号化形式としてはＭＩＤＩ形式が最も普及しているた
め、実用上はＭＩＤＩ形式の符号データを用いるのが好
ましい。ＭＩＤＩ形式では、「ノートオン」データもし
くは「ノートオフ」データが、「デルタタイム」データ
を介在させながら存在する。「ノートオン」データは、
特定のノートナンバーＮとベロシティーＶを指定して特
定の音の演奏開始を指示するデータであり、「ノートオ
フ」データは、特定のノートナンバーＮとベロシティー
Ｖを指定して特定の音の演奏終了を指示するデータであ
る。また、「デルタタイム」データは、所定の時間間隔
を示すデータである。ベロシティーＶは、例えば、ピア
ノの鍵盤などを押し下げる速度（ノートオン時のベロシ
ティー）および鍵盤から指を離す速度（ノートオフ時の
ベロシティー）を示すパラメータであり、特定の音の演
奏開始操作もしくは演奏終了操作の強さを示すことにな
る。

【００１８】前述の方法では、第ｉ番目の単位区間ｄｉ
について、代表符号としてＪ個のノートナンバーｎ（ｄ
ｉ，１），ｎ（ｄｉ，２），・・・，ｎ（ｄｉ，Ｊ）が
得られ、このそれぞれについて強度ｅ（ｄｉ，１），ｅ
（ｄｉ，２），・・・，ｅ（ｄｉ，Ｊ）が得られる。そ
こで、次のような手法により、ＭＩＤＩ形式の符号デー
タを作成することができる。まず、「ノートオン」デー
タもしくは「ノートオフ」データの中で記述するノート
ナンバーＮとしては、得られたノートナンバーｎ（ｄ
ｉ，１），ｎ（ｄｉ，２），・・・，ｎ（ｄｉ，Ｊ）を
そのまま用いれば良い。一方、「ノートオン」データも
しくは「ノートオフ」データの中で記述するベロシティ
ーＶとしては、得られた強度ｅ（ｄｉ，１），ｅ（ｄ
ｉ，２），・・・，ｅ（ｄｉ，Ｊ）を所定の方法で規格
化した値を用いれば良い。また、「デルタタイム」デー
タは、各単位区間の長さに応じて設定すれば良い。

【００１９】（2.1.周期関数との相関を求める具体的な
方法）上述した基本原理の基づく方法では、区間信号に
対して、１つまたは複数の代表周波数が選出され、この
代表周波数をもった周期信号によって、当該区間信号が
表現されることになる。ここで、選出される代表周波数
は、文字どおり、当該単位区間内の信号成分を代表する
周波数である。この代表周波数を選出する具体的な方法
には、後述するように、短時間フーリエ変換を利用する
方法と、一般化調和解析の手法を利用する方法とがあ
る。いずれの方法も、基本的な考え方は同じであり、あ
らかじめ周波数の異なる複数の周期関数を用意してお
き、これら複数の周期関数の中から、当該単位区間内の
区間信号に対する相関が高い周期関数を見つけ出し、こ
の相関の高い周期関数の周波数を代表周波数として選出
する、という手法を採ることになる。すなわち、代表周
波数を選出する際には、あらかじめ用意された複数の周
期関数と、単位区間内の区間信号との相関を求める演算
を行うことになる。そこで、ここでは、周期関数との相
関を求める具体的な方法を述べておく。

【００２０】複数の周期関数として、図２に示すような
三角関数が用意されているものとする。これらの三角関
数は、同一周波数をもった正弦関数と余弦関数との対か
ら構成されており、１２８通りの標準周波数ｆ（０）〜
ｆ（１２７）のそれぞれについて、正弦関数および余弦
関数の対が定義されていることになる。ここでは、同一
の周波数をもった正弦関数および余弦関数からなる一対
の関数を、当該周波数についての周期関数として定義す
ることにする。すなわち、ある特定の周波数についての
周期関数は、一対の正弦関数および余弦関数によって構
成されることになる。このように、一対の正弦関数と余
弦関数とにより周期関数を定義するのは、信号に対する
周期関数の相関値を求める際に、相関値が位相の影響を
受ける事を考慮するためである。なお、図２に示す各三
角関数内の変数Ｆおよびｋは、区間信号Ｘについてのサ
ンプリング周波数Ｆおよびサンプル番号ｋに相当する変
数である。例えば、周波数ｆ（０）についての正弦波
は、ｓｉｎ（２πｆ（０）ｋ／Ｆ）で示され、任意のサ
ンプル番号ｋを与えると、区間信号を構成する第ｋ番目
のサンプルと同一時間位置における周期関数の振幅値が
得られる。

【００２１】ここでは、１２８通りの標準周波数ｆ
（０）〜ｆ（１２７）を図３に示すような式で定義した
例を示すことにする。すなわち、第ｎ番目（０≦ｎ≦１
２７）の標準周波数ｆ（ｎ）は、以下に示す〔数式１〕
で定義されることになる。

【００２２】〔数式１〕ｆ（ｎ）＝４４０×２^γ ⁽ⁿ⁾ γ（ｎ）＝（ｎ−６９）／１２

【００２３】このような式によって標準周波数を定義し
ておくと、最終的にＭＩＤＩデータを用いた符号化を行
う際に便利である。なぜなら、このような定義によって
設定される１２８通りの標準周波数ｆ（０）〜ｆ（１２
７）は、等比級数をなす周波数値をとることになり、Ｍ
ＩＤＩデータで利用されるノートナンバーに対応した周
波数になるからである。したがって、図２に示す１２８
通りの標準周波数ｆ（０）〜ｆ（１２７）は、対数尺度
で示した周波数軸上に等間隔（ＭＩＤＩにおける半音単
位）に設定した周波数ということになる。

【００２４】続いて、任意の区間の区間信号に対する各
周期関数の相関の求め方について、具体的な説明を行
う。例えば、図４に示すように、ある単位区間ｄについ
て区間信号Ｘが与えられていたとする。ここでは、区間
長Ｌをもった単位区間ｄについて、サンプリング周波数
Ｆでサンプリングが行なわれており、全部でｗ個のサン
プル値が得られているものとし、サンプル番号を図示の
ように、０，１，２，３，・・・，ｋ，・・・，ｗ−
２，ｗ−１とする（白丸で示す第ｗ番目のサンプルは、
右に隣接する次の単位区間の先頭に含まれるサンプルと
する）。この場合、任意のサンプル番号ｋについては、
Ｘ（ｋ）なる振幅値がデジタルデータとして与えられて
いることになる。短時間フーリエ変換においては、Ｘ
（ｋ）に対して各サンプルごとに中央の重みが１に近
く、両端の重みが０に近くなるような窓関数Ｗ（ｋ）を
乗ずることが通常である。すなわち、Ｘ（ｋ）×Ｗ
（ｋ）をＸ（ｋ）と扱って以下のような相関計算を行う
もので、窓関数の形状としては余弦波形状のハミング窓
が一般に用いられている。ここで、ｗは以下の記述にお
いても定数のような記載をしているが、一般にはｎの値
に応じて変化させ、区間長Ｌを超えない範囲で最大とな
るＦ／ｆ（ｎ）の整数倍の値に設定することが望まし
い。

【００２５】このような区間信号Ｘに対して、第ｎ番目
の標準周波数ｆ（ｎ）をもった正弦関数Ｒｎとの相関値
を求める原理を示す。両者の相関値Ａ（ｎ）は、図５の
第１の演算式によって定義することができる。ここで、
Ｘ（ｋ）は、図４に示すように、区間信号Ｘにおけるサ
ンプル番号ｋの振幅値であり、ｓｉｎ（２πｆ（ｎ）ｋ
／Ｆ）は、時間軸上での同位置における正弦関数Ｒｎの
振幅値である。この第１の演算式は、単位区間ｄ内の全
サンプル番号ｋ＝０〜ｗ−１の次元について、それぞれ
区間信号Ｘの振幅値と正弦関数Ｒｎの振幅ベクトルの内
積を求める式ということができる。

【００２６】同様に、図５の第２の演算式は、区間信号
Ｘと、第ｎ番目の標準周波数ｆ（ｎ）をもった余弦関数
との相関値を求める式であり、両者の相関値はＢ（ｎ）
で与えられる。なお、相関値Ａ（ｎ）を求めるための第
１の演算式も、相関値Ｂ（ｎ）を求めるための第２の演
算式も、最終的に２／ｗが乗ぜられているが、これは相
関値を規格化するためのものでり、前述のとおりｗはｎ
に依存して変化させるのが一般的であるため、この係数
もｎに依存する変数である。

【００２７】区間信号Ｘと標準周波数ｆ（ｎ）をもった
標準周期関数との相関実効値は、図５の第３の演算式に
示すように、正弦関数との相関値Ａ（ｎ）と余弦関数と
の相関値Ｂ（ｎ）との二乗和平方根値Ｅ（ｎ）によって
示すことができる。この相関実効値の大きな標準周期関
数の周波数を代表周波数として選出すれば、この代表周
波数を用いて区間信号Ｘを符号化することができる。

【００２８】すなわち、この相関値Ｅ（ｎ）が所定の基
準以上の大きさとなる１つまたは複数の標準周波数を代
表周波数として選出すれば良い。なお、ここで「相関値
Ｅ（ｎ）が所定の基準以上の大きさとなる」という選出
条件は、例えば、何らかの閾値を設定しておき、相関値
Ｅ（ｎ）がこの閾値を超えるような標準周波数ｆ（ｎ）
をすべて代表周波数として選出する、という絶対的な選
出条件を設定しても良いが、例えば、相関値Ｅ（ｎ）の
大きさの順にＱ番目までを選出する、というような相対
的な選出条件を設定しても良い。

【００２９】（2.2.一般化調和解析の手法）ここでは、
本発明に係る音響信号の符号化を行う際に有用な一般化
調和解析の手法について説明する。既に説明したよう
に、音響信号を符号化する場合、個々の単位区間内の区
間信号について、相関値の高いいくつかの代表周波数を
選出することになる。一般化調和解析は、より高い精度
で代表周波数の選出を可能にする手法であり、その基本
原理は次の通りである。

【００３０】図６（ａ）に示すような単位区間ｄについ
て、信号Ｓ（ｊ）なるものが存在するとする。ここで、
ｊは後述するように、繰り返し処理のためのパラメータ
である（ｊ＝１〜Ｊ）。まず、この信号Ｓ（ｊ）に対し
て、図２に示すような１２８通りの周期関数すべてにつ
いての相関値を求める。そして、最大の相関値が得られ
た１つの周期関数の周波数を代表周波数として選出し、
当該代表周波数をもった周期関数を要素関数として抽出
する。続いて、図６（ｂ）に示すような含有信号Ｇ
（ｊ）を定義する。この含有信号Ｇ（ｊ）は、抽出され
た要素関数に、その振幅として、当該要素関数の信号Ｓ
（ｊ）に対する相関値を乗じることにより得られる信号
である。例えば、周期関数として図２に示すように、一
対の正弦関数と余弦関数とを用い、周波数ｆ（ｎ）が代
表周波数として選出された場合、振幅Ａ（ｎ）をもった
正弦関数Ａ（ｎ）ｓｉｎ（２πｆ（ｎ）ｋ／Ｆ）と、振
幅Ｂ（ｎ）をもった余弦関数Ｂ（ｎ）ｃｏｓ（２πｆ
（ｎ）ｋ／Ｆ）との和からなる信号が含有信号Ｇ（ｊ）
ということになる（図６（ｂ）では、図示の便宜上、一
方の関数しか示していない）。ここで、Ａ（ｎ），Ｂ
（ｎ）は、図５の式で得られる規格化された相関値であ
るから、結局、含有信号Ｇ（ｊ）は、信号Ｓ（ｊ）内に
含まれている周波数ｆ（ｎ）をもった信号成分というこ
とができる。

【００３１】こうして、含有信号Ｇ（ｊ）が求まった
ら、信号Ｓ（ｊ）から含有信号Ｇ（ｊ）を減じることに
より、差分信号Ｓ（ｊ＋１）を求める。図６（ｃ）は、
このようにして求まった差分信号Ｓ（ｊ＋１）を示して
いる。この差分信号Ｓ（ｊ＋１）は、もとの信号Ｓ
（ｊ）の中から、周波数ｆ（ｎ）をもった信号成分を取
り去った残りの信号成分からなる信号ということができ
る。そこで、パラメータｊを１だけ増加させることによ
り、この差分信号Ｓ（ｊ＋１）を新たな信号Ｓ（ｊ）と
して取り扱い、同様の処理を、パラメータｊをｊ＝１〜
Ｊまで１ずつ増やしながらＪ回繰り返し実行すれば、Ｊ
個の代表周波数を選出することができる。

【００３２】このような相関計算の結果として出力され
るＪ個の含有信号Ｇ（１）〜Ｇ（Ｊ）は、もとの区間信
号Ｘの構成要素となる信号であり、もとの区間信号Ｘを
符号化する場合には、これらＪ個の含有信号の周波数を
示す情報および振幅（強度）を示す情報を符号データと
して用いるようにすれば良い。尚、Ｊは代表周波数の個
数であると説明してきたが、標準周波数ｆ（ｎ）の個数
と同一すなわちＪ＝１２８であってもよく、周波数スペ
クトルを求める目的においてはそのように行うのが通例
である。

【００３３】以上のような処理により、各単位区間につ
いて、各周波数に対する強度値の集合である周波数群が
得られることになる。このようにして所定数の周波数群
が選出されたら、この周波数群の各周波数に対応する
「音の高さを示す情報」、選出された各周波数の信号強
度に対応する「音の強さを示す情報」、当該単位区間の
始点に対応する「音の発音開始時刻を示す情報」、当該
単位区間に後続する単位区間の始点に対応する「音の発
音終了時刻を示す情報」、の４つの情報を含む所定数の
符号データを作成すれば、当該単位区間内の区間信号Ｘ
を所定数の符号データにより符号化することができる。
符号データとして、ＭＩＤＩデータを作成するのであれ
ば、「音の高さを示す情報」としてノートナンバーを用
い、「音の強さを示す情報」としてベロシティーを用
い、「音の発音開始時刻を示す情報」としてノートオン
時刻を用い、「音の発音終了時刻を示す情報」としてノ
ートオフ時刻を用いるようにすれば良い。

【００３４】（3.1.本発明に係る周波数解析方法および
音響信号の符号化方法）以下、本発明に係る周波数解析
方法および音響信号の符号化方法について説明してい
く。上述のように、従来の手法では、微細周波数に対す
る相関値を求め、各標準周波数の範囲において、最大の
相関値を当該標準周波数の相関値としている。例えば、
原音響信号のスペクトルが図７（ａ）に示すような状態
となるものを考えてみる。なお、図７において、横軸は
周波数、縦軸は強度値もしくは相関値である。図７
（ａ）においては、ノートナンバーｎの範囲（標準周波
数ｆ（ｎ）の範囲）に含まれる音と、ノートナンバーｎ
＋１の範囲に含まれる音が発せられていることになる。
これを短時間フーリエ変換の手法により、各微細周波数
についての相関を求めると、図７（ｂ）に示すようにな
る。図７（ｂ）においては、各標準周波数間を１３に分
割した間隔で微細周波数が設定されている。従来の符号
化では、ここで各ノートナンバーの範囲内において、最
大の相関値となるものを抽出する。この際、各ノートナ
ンバーの範囲において、境界部分に存在するものは、他
の信号成分の影響であるとみなして抽出しないようにし
ている。そのため、図７（ｃ）に示すようにノートナン
バーｎからは最大の相関値が抽出されるが、ノートナン
バーｎ＋１からは抽出されないことになる。

【００３５】これに対して本発明では、各ノートナンバ
ーの範囲において最大の相関値を抽出するのではなく、
ピークとなる相関値を抽出するようにする。ピークとな
る相関値とは、自身の両側の微細周波数の相関値が自身
よりも小さくなっているものをいう。ピークとなる相関
値を抽出した場合には、それ自身がピークであるために
他の信号成分の影響を受けているということは考えにく
い。そのため、ノートナンバーの範囲の境界部分に存在
していても抽出するものとする。したがって、図７
（ｂ）に示したような相関値の集合から図７（ｄ）に示
すような信号成分が抽出されることになる。

【００３６】また、２つ目の例として図８に示すような
ものを考えてみる。図８（ａ）に示すような２つの音が
発せられたとすると、短時間フーリエ変換の手法によ
り、図８（ｂ）に示すような各微細周波数に対応する相
関値が得られる。従来の手法では、図８（ｃ）に示すよ
うに、ノートナンバーｎからは抽出されるが、ノートナ
ンバーｎ＋１からは抽出されないということになるが、
本発明によれば、図８（ｄ）に示すようにそれぞれのノ
ートナンバーに対応した周波数成分が抽出されることに
なる。

【００３７】続いて、本発明の周波数解析方法につい
て、音響信号を符号化する場合を例にとって、図９に示
すフローチャートを用いて説明する。まず、時系列の音
響信号に対して単位区間を設定して区間信号の抽出を行
う（ステップＳ１）。これは、上記（基本原理）の項に
おいて、図１を用いて説明したように、所定の長さの単
位区間が設定され、抽出が行われる。

【００３８】次に、抽出された区間信号とあらかじめ準
備された調和信号との相関を算出する（ステップＳ
２）。この調和信号とは、上記（周期関数との相関を求
める具体的な方法）の項で説明した周期関数と同じもの
である。ただし、図２の例では、ノートナンバーに対応
した１２８個の標準周波数を設定したが、本発明では、
各標準周波数間に等比級数の間隔で１２個の微細周波数
を設定し、全部で１２８×１３個の調和信号を用意す
る。例えば、隣接する標準周波数ｆ（ｎ）とｆ（ｎ＋
１）間では、ｆ（ｎ＋１／１３）〜ｆ（ｎ＋１２／１
３）の１２個の微細周波数が設定されることになる。こ
のステップＳ２においては、区間信号と１２８×１３個
の調和信号との相関を短時間フーリエ変換の手法により
求める。この結果、１２８×１３個の各微細周波数に対
応する相関値が得られる。この相関値は、あらかじめ準
備された微細相関配列に格納される。

【００３９】続いて、算出された相関値に基づいてピー
クとなる微細周波数（以下、ピーク周波数という）を選
出する（ステップＳ３）。具体的には、算出された１２
８×１３個の相関値を基に、隣接する両側の微細周波数
の相関値よりも、自身の相関値が大きくなっている微細
周波数を選出する。

【００４０】次に、選出された各ピーク周波数に対応す
る調和信号と、区間信号との相関を一般化調和解析の手
法により求める（ステップＳ４）。具体的には、上記
（一般化調和解析の手法）の項で説明したのと同様に、
まず、選出されたピーク周波数のうち、相関値が最大の
ものに対応する調和信号と区間信号Ｓ（１）（＝Ｘ）と
の相関を求める。続いて、このときの調和信号に、その
振幅として算出された相関値を乗じることにより得られ
る信号を含有信号Ｇ（１）とし、この含有信号Ｇ（１）
を区間信号Ｓ（１）から減じることにより差分信号Ｓ
（２）を得る。次に、ステップＳ３において選出された
ピーク周波数のうち、区間信号Ｘとの相関値が２番目に
大きいピーク周波数の調和信号と、差分信号Ｓ（２）と
の相関を求め、含有信号Ｇ（２）を区間信号から減じる
ことにより差分信号Ｓ（３）を得る。このようにして、
一般化調和解析の手法により、選出された全てのピーク
周波数の相関値を求めていく。このピーク周波数は、微
細周波数に対応しているため、符号化の際には、ノート
ナンバーに対応した標準周波数に戻す必要がある。その
ため、算出された相関値を、ピーク周波数に最も近い標
準周波数に対応する相関値として、主相関配列に設定す
る。これにより、標準周波数とその相関値の組が複数得
られることになり、この標準周波数をノートナンバーと
し、相関値をベロシティとすることにより、ＭＩＤＩ規
格の符号化データが得られる。

【００４１】以上の処理を全単位区間により行うことに
より、音響信号全体の符号化が可能となる。なお、単位
区間については、図１に示したように、前の単位区間の
終了時刻と後の単位区間の開始時刻が同一となるように
設定されるのは、まれであり、前の単位区間と後の単位
区間が重なるように設定されることの方が多い。

【００４２】（3.2.抽出周波数成分決定の２段階処理）
上記ステップＳ３・Ｓ４においては、１回だけピーク周
波数を選出し、選出された各ピーク周波数について一般
化調和解析の手法により相関値を算出したが、同様の処
理を再度行うことにより解析精度および符号化精度を高
めることができる。このような場合について、以下説明
する。上記ステップＳ４のように区間信号Ｘから含有信
号を減じる処理を繰り返し行っていくことにより、異な
った区間信号Ｘ’が得られることになる。この区間信号
Ｘ’に対してステップＳ２〜ステップＳ４の処理を繰り
返す。

【００４３】区間信号Ｘ’と相関をとる微細周波数とし
ては、既に主相関配列に相関値が格納された標準周波数
に該当するものは除かれる。例えば、標準周波数ｆ
（ｎ）に対応する相関値が主相関配列に設定されている
場合には、微細周波数ｆ（ｎ−６／１３）〜ｆ（ｎ＋６
／１３）は、区間信号Ｘ’との相関算出の対象から外さ
れる。残った微細周波数に対応する調和信号と区間信号
Ｘ’との相関を短時間フーリエ変換の手法により求めら
れ、算出された相関値は微細相関配列に格納される。

【００４４】続いて、ステップＳ３において、算出され
た相関値に基づいて再びピーク周波数を選出した後、ス
テップＳ４において、選出された各ピーク周波数と区間
信号Ｘ’との相関を求めた後、算出された相関値を主相
関配列に設定する。

【００４５】ステップＳ２〜ステップＳ４の処理は、何
回でも繰り返すことが可能であるが、音響信号の符号化
処理を行う場合には、２回程度で十分である。２回行う
と、ＭＩＤＩ規格で同時発音可能な６４音程度は抽出で
きるためである。

【００４６】（3.3.同一標準周波数内における抽出周波
数成分の決定）上記ステップＳ４における標準周波数の
決定および相関値の算出は、ピーク周波数の近傍にある
標準周波数を、算出された相関値に対応するものとして
決定した。これは、１つの標準周波数の範囲に含まれる
ピーク周波数が１つだけである場合には問題ないが、１
つの標準周波数の範囲に含まれるピーク周波数が複数あ
る場合には、そのどちらを選ぶかが問題となる。そのた
め、後者の場合には、１つの標準周波数に含まれるピー
ク周波数のうち、短時間フーリエ変換による相関値が最
大のものだけをピーク周波数として選出し、一般化調和
解析の手法により区間信号Ｘとの相関を求めるようにす
る。

【００４７】（3.4.区間更新幅の決定処理）以上のよう
な周波数成分の抽出処理および抽出した周波数成分の符
号化は各単位区間ごとに行われるが、上述のように単位
区間の設定は固定ではなく、状況に応じて変更すること
が可能である。このとき、１つの単位区間の長さをどの
ように設定するかは、解析結果に大きな影響を与える重
要な問題である。すなわち、区間長を短く設定すればす
るほど、時間分解能が向上するため、音声信号（ボーカ
ル信号）などの周波数変化を解析するには好都合である
が、逆に、周波数分解能は低下するため、解析精度は低
下するという問題が生じることになる。そこで従来は、
１つの単位区間の中に、解析に最も都合が良いと思われ
るサンプル数（例えば、１０２４サンプル）が含まれる
ように、標準的な区間長をあらかじめ定めておき、この
固定の区間長に基づいて時間軸を機械的に区切ってゆく
という方法を採るのが一般的であった。

【００４８】しかしながら、固定長の単位区間を機械的
に設定するという方法では、解析対象となる時系列信号
がどのような信号であるかに関わらず、常に一義的な区
間設定が行われることになるため、単位区間の境界部分
において不連続性が生じるおそれがある。すなわち、隣
接する単位区間についての周波数解析の結果に大きな差
が生じてしまうことがある。一般に、得られた周波数解
析の結果を、スペクトルのグラフとして観察するために
利用するような場合には、単位区間の境界において解析
結果に差が生じていても、大きな問題にはならない。と
ころが、音響信号に対する解析結果を利用して、原音響
信号に対応する符号データを作成するような場合、境界
部分における不連続性は大きな問題となる。このような
方法で符号化されたデータに基づいて、原音響信号を再
生しようとすると、単位区間の境界において、音程の急
激な変動が生じたり、音とびが生じたりすることにな
る。

【００４９】このような境界における不連続性を避ける
ために、時間軸上において単位区間を部分的に重複させ
て設定する方法も提案されているが、この場合も重複さ
せる量が固定であるため、本来重複させる必要のない箇
所まで重複して処理することになり、無駄な演算処理が
増加するという問題が生じることになる。本発明では、
時系列信号を解析する際に、単位区間の境界において、
音の連続性を保つと共に、単位区間の重複を最小限に抑
えるようにした単位区間の設定についても工夫を行って
いるので、以下に説明する。

【００５０】特に本発明においては、上記ステップＳ２
における短時間フーリエ変換により得られた相関値を利
用することにより単位区間の設定を行うようにしたこと
を特徴としている。具体的には、ある単位区間において
周波数成分の抽出を行う際に、ステップＳ２において最
初に得られる微細相関配列の値を、区間決定のために用
いる直前相関配列にコピーしておく。そして、ステップ
Ｓ２〜ステップＳ４の処理を所定の回数だけ行って、１
つの単位区間において抽出すべき周波数成分を決定した
後、次の単位区間の位置を設定する。本発明において
は、単位区間は固定長であるため、単位区間の開始位置
が定まれば、単位区間が決定する。従って本発明におい
ては、前の単位区間の開始位置から後の単位区間の開始
位置までの送り幅を定めることにより、後の単位区間が
決定されることになる。

【００５１】まず、既に確定している前の単位区間と、
これから位置を決定する単位区間の類似度を求める（ス
テップＳ６）。具体的には、ピーク周波数選出のため
に、ステップＳ２で行った短時間フーリエ変換の結果得
られた微細相関配列の値を利用して以下の〔数式２〕に
より算出する。

【００５２】〔数式２〕ＳＭ＝100−100×［Σk｛Ｓ_dn-1(k)−Ｓ_dn(k)｝²］^1/2
／｛ΣkＳ_dn-1(k)²｝^1/ ²

【００５３】上記〔数式２〕において、Ｓ_dn(k)は、決
定しようとしている単位区間ｄnの区間信号に対する微
細周波数ｋの相関値であり、Ｓ_dn-1(k)は、直前の単位
区間ｄn-1の区間信号に対する微細周波数ｋの相関値で
ある。ただし、ｋは周波数そのものを示すものではな
く、ノートナンバーｎに対応する標準周波数を１３分割
した値で表現され、ｋ＝−６〜１２７×１３＋６の
値をとる。類似度ＳＭの値が大きい程両区間信号は類似
していることになる。

【００５４】次に、この類似度ＳＭの値と閾値を比較す
る。類似度ＳＭが閾値２より小さい場合は（ステップＳ
７）、設定されている送り幅を半分にする（ステップＳ
８）。送り幅を半分にした位置に仮単位区間の設定をし
直し（ステップＳ１）、新たな仮単位区間の区間信号と
調和信号との相関を求める（ステップＳ２）。ステップ
Ｓ１、Ｓ２およびステップＳ６〜Ｓ８の処理は、類似度
ＳＭが閾値２以上となるまで繰り返される。すなわち、
閾値２以上となるまで、後の単位区間を前の単位区間と
重複させる範囲を大きくしていくことになる。

【００５５】類似度ＳＭが閾値２以上となったら、類似
度ＳＭと閾値１との比較を行う。類似度が閾値１未満で
ある場合は、仮単位区間を単位区間として確定し、ピー
ク位置の選出処理へ進む（ステップＳ３）。類似度が閾
値１以上である場合は、送り幅を２倍に設定する（ステ
ップＳ１０）と共に、仮単位区間を単位区間として確定
し、ピーク位置の選出処理へ進む（ステップＳ３）。こ
こで、２倍に変更された送り幅は、次の単位区間の送り
幅として利用されることになる。単位区間の確定後は、
ステップＳ３、Ｓ４の処理は同様に行われ、抽出周波数
成分が決定される。なお、当然のことながら、上記２つ
の閾値は、閾値１＞閾値２の関係を満たすものとなって
いる。

【００５６】単位区間についての抽出周波数が決定した
ら、この時点で設定されている送り幅だけ移動した位置
を開始位置として（ステップＳ５）、仮単位区間を設定
する（ステップＳ１）。このとき、通常はあらかじめ設
定されている標準送り幅が利用され、ステップＳ１０に
より送り幅が２倍に設定されている場合には、標準送り
幅の２倍を送り幅とする。

【００５７】ここで、設定される単位区間を、図１０に
示す具体的な例を用いて説明する。図１０（ａ）は、時
系列の音響信号を示す図である。図１０（ｂ）に示すよ
うに音響信号の先頭位置から固定長の単位区間ｄ１が設
定され、ステップＳ２からステップＳ４の処理が行われ
ることにより、抽出周波数が決定される。次の単位区間
ｄ２の設定にあたっては、ステップＳ５において、あら
かじめ設定されている標準送り幅だけ後に送られて仮設
定される。本実施形態では、標準送り幅を固定長ｄの１
／２としている。よって、単位区間ｄ１の先頭からｄ／
２だけずれた位置が単位区間ｄ２の先頭となるように単
位区間ｄ２が仮設定される。このとき仮設定された単位
区間ｄ２の状態を図１０（ｃ）に示す。

【００５８】そして、ステップＳ１において、仮設定さ
れた単位区間ｄ２の位置の音響信号を区間信号として抽
出し、この区間信号との相関をステップＳ２の手法によ
り求める。続いてステップＳ６において上記〔数式２〕
に示した類似度を算出する。算出した類似度が閾値２未
満である場合、すなわち、両区間信号があまり似ていな
い場合は、ステップＳ８に進んで送り幅が１／２にされ
る。ここでは、ｄ／２に設定されているので、ｄ／４に
設定されるということになる。この場合、ステップＳ１
に戻って、単位区間ｄ１の先頭からｄ／４だけずれた位
置が単位区間ｄ２の先頭となるように単位区間ｄ２が仮
設定し直される。このとき仮設定された単位区間ｄ２の
状態を図１０（ｄ）に示す。単位区間ｄ２が仮設定し直
されると、再びステップＳ２、ステップＳ６の処理が繰
り返される。さらに、ステップＳ７において、類似度が
閾値２未満である場合は、ステップＳ８において送り幅
がｄ／８に設定され、ステップＳ１に戻って単位区間ｄ
１の先頭からｄ／８だけずれた位置が単位区間ｄ２の先
頭となるように単位区間ｄ２が仮設定し直される。この
とき仮設定された単位区間ｄ２の状態を図１０（ｅ）に
示す。

【００５９】すなわち、図１０（ｃ）〜図１０（ｅ）に
示すように、先行する単位区間ｄ１と後続の単位区間ｄ
２における区間信号の類似度が閾値２以上となるまで、
単位区間ｄ２は単位区間ｄ１との重複区間を増やしてい
くことになる。ただし、このままだと、単位区間ｄ１と
単位区間ｄ２が果てしなく近付いてしまう可能性がある
ので、送り幅の最小距離をあらかじめ設定しておき、送
り幅が最小距離以下となった場合は、その時点で仮設定
されている単位区間を正式な単位区間として確定し、ス
テップＳ３の処理に進むようにする。

【００６０】一方、単位区間ｄ１の区間信号と、仮設定
された単位区間ｄ２の区間信号の類似度が閾値２以上で
閾値１未満となった場合には、仮設定されている単位区
間ｄ２を確定して、ステップＳ３以降の処理を行う。ま
た、類似度が閾値１以上となった場合も、仮設定されて
いる単位区間ｄ２が確定されて、ステップＳ３以降の処
理が行われることになるが、この場合、送り幅を２倍に
設定する処理も平行して行う。そして、単位区間ｄ２に
ついての抽出周波数成分が決定された後、ステップＳ５
において、単位区間の先頭から２倍に設定された送り幅
だけ移動した位置を単位区間ｄ３の先頭位置とする。例
えば、単位区間ｄ２が図１０（ｃ）に示すような位置で
確定された場合、単位区間ｄ３は、図１０（ｆ）に示す
位置に仮設定されることになる。

【００６１】以上、本発明の好適な実施形態について説
明したが、上記周波数解析方法および符号化方法は、コ
ンピュータ等で実行されることは当然である。具体的に
は、図９のフローチャートに示したようなステップを上
記手順で実行するためのプログラムをコンピュータに搭
載しておく。そして、音響信号等の時系列信号をＰＣＭ
方式等でデジタル化した後、コンピュータに取り込み、
ステップＳ１〜ステップＳ１０の処理を行った後、抽出
した周波数成分もしくはＭＩＤＩ形式等の符号データを
デジタルデータとしてコンピュータより出力する。出力
された符号データは、例えば、ＭＩＤＩデータの場合、
ＭＩＤＩシーケンサ、ＭＩＤＩ音源を用いて音声として
再生される。

【００６２】

【発明の効果】以上、説明したように本発明によれば、
時系列信号から所定の標準周波数に対応した信号成分を
分離するための周波数解析方法として、前記時系列信号
に所定の長さの単位区間を設定して区間信号を抽出し、
前記全ての調和信号と前記区間信号との相関を演算して
各微細周波数に対応する相関値を算出し、前記算出され
た相関値に対応する微細周波数のうち、周囲の微細周波
数における相関値に比べ局所的にピークとなる微細周波
数を複数抽出するピーク周波数選出段階と、前記各ピー
クとなる微細周波数に対応する調和信号と前記区間信号
との相関を再度演算し、前記微細周波数の最近傍の標準
周波数に対応する相関値として設定し、前記再度算出さ
れた相関値と調和信号との積で生成される含有信号を前
記区間信号から減じ、差分信号を新たに区間信号とする
ことにより区間信号を更新する処理を抽出された全ての
ピークに対して順次行って標準周波数に対応する相関値
を決定するようにしたので、隣接する標準周波数間の境
界付近に位置する微細周波数がピークとなる場合であっ
ても、正確に周波数成分を抽出することが可能となると
いう効果を奏する。さらに、抽出された周波数成分をＭ
ＩＤＩ等の形式に符号化することにより、原音響信号を
忠実に再現することが可能な符号データの作成が可能と
なる。

【図面の簡単な説明】

【図１】本発明に係る周波数解析方法および音響信号符
号化方法の基本原理を示す図である。

【図２】本発明で利用される周期関数の一例を示す図で
ある。

【図３】図２に示す各周期関数の周波数とＭＩＤＩノー
トナンバーｎとの関係式を示す図である。

【図４】解析対象となる信号と周期信号との相関計算の
手法を示す図である。

【図５】図４に示す相関計算を行うための計算式を示す
図である。

【図６】一般化調和解析の基本的な手法を示す図であ
る。

【図７】周波数成分の選出に関して従来との比較を示す
第１の例を示す図である。

【図８】周波数成分の選出に関して従来との比較を示す
第２の例を示す図である。

【図９】本発明に係る周波数解析方法および音響信号符
号化方法のフローチャートである。

【図１０】単位区間の設定手順を示す図である。

【符号の説明】

ｄ１〜ｄ５・・・単位区間ｄ・・・単位区間長Ｘ・・・初期区間信号

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.⁷ 識別記号ＦＩテーマコート゛(参考）Ｇ１０Ｌ 9/18 Ｍ

Claims

【特許請求の範囲】

【請求項１】時系列信号から所定の標準周波数に対応し
た信号成分を分離するための周波数解析方法であって、前記標準周波数の間隔よりも微細な間隔で複数の微細周
波数を設定し、各微細周波数に対応する複数の調和信号
集合を準備する調和信号準備段階と、前記各標準周波数に対応する相関値を格納するための主
相関配列と、各微細周波数に対応する相関値を格納する
ための微細相関配列を準備する配列準備段階と、前記時系列信号に所定の長さの単位区間を設定し、区間
信号を抽出するための区間信号抽出段階と、前記全ての調和信号と前記区間信号との相関を演算し
て、各微細周波数に対応する相関値を算出し、前記微細
相関配列に算出された相関値を設定する相関演算段階
と、前記微細相関配列中の相関配列より周囲の微細周波数に
おける相関値に比べ局所的にピークとなる微細周波数を
複数抽出するピーク周波数選出段階と、前記各ピークとなる微細周波数に対応する調和信号と前
記区間信号との相関を再度演算し、前記微細周波数の最
近傍の標準周波数に対応する主相関配列に設定し、前記
再度算出された相関値と調和信号との積で生成される含
有信号を前記区間信号から減じ、差分信号を新たに区間
信号とすることにより区間信号を更新する処理を抽出さ
れた全てのピークに対して順次行って主相関配列の値を
決定するスペクトル算出段階と、を有することを特徴とする周波数解析方法。
【請求項２】前記相関演算段階、前記ピーク周波数選出
段階、および前記スペクトル算出段階を複数回実行する
ようにしていることを特徴とする請求項１に記載の周波
数解析方法。
【請求項３】さらに、類似度算出段階を有し、前記配列準備段階が、直前の単位区間において得られた
微細相関配列を保持するための区間決定用相関配列を準
備し、前記類似度算出段階が、直前の単位区間における微細相
関配列の値を複写して得られる区間決定用相関配列の値
と、仮設定された現単位区間における微細相関配列の値
の類似度を算出し、当該類似度に基づいて単位区間の送
り幅を決定するものであり、前記区間信号抽出段階は、前記決定された送り幅に基づ
く位置に設定される単位区間から区間信号を抽出するも
のであるものであることを特徴とする請求項１に記載の
周波数解析方法。
【請求項４】前記類似度算出段階により算出される類似
度が所定値より低い場合には送り幅を小さくなるように
設定するものであることを特徴とする請求項３に記載の
周波数解析方法。
【請求項５】前記スペクトル算出段階は、１つの標準周
波数が最近傍となる複数の微細周波数がピークとして選
出された場合に、そのうち最大の相関値を有する微細周
波数についてのみ、区間信号との相関を求めるようにす
るものであることを特徴とする請求項１に記載の周波数
解析方法。
【請求項６】前記時系列信号は音響信号であって、前記スペクトル算出段階によって選出された標準周波数
に対応する音の高さ情報と、相関値に対応する音の強さ
情報と、各単位区間の開始点に対応する音の発音開始時
刻と、各単位区間の終了点に対応する音の発音終了時刻
からなる４つの情報に所定の変換を施して符号データを
生成する符号化段階をさらに有することを特徴とする請
求項１から請求項５のいずれかに記載の音響信号の符号
化方法。
【請求項７】コンピュータに、前記標準周波数の間隔よ
りも微細な間隔で複数の微細周波数を設定して各微細周
波数に対応する複数の調和信号集合を準備する調和信号
準備段階、前記各標準周波数に対応する相関値を格納す
るための主相関配列と各微細周波数に対応する相関値を
格納するための微細相関配列を準備する配列準備段階、
前記時系列信号に所定の長さの単位区間を設定して区間
信号を抽出するための区間信号抽出段階、前記全ての調
和信号と前記区間信号との相関を演算して、各微細周波
数に対応する相関値を算出し、前記微細相関配列に算出
された相関値を設定する相関演算段階、前記微細相関配
列中の相関配列より周囲の微細周波数における相関値に
比べ局所的にピークとなる微細周波数を複数抽出するピ
ーク周波数選出段階、前記各ピークとなる微細周波数に
対応する調和信号と前記区間信号との相関を再度演算
し、前記微細周波数の最近傍の標準周波数に対応する主
相関配列に設定し、前記再度算出された相関値と調和信
号との積で生成される含有信号を前記区間信号から減
じ、差分信号を新たに区間信号とすることにより区間信
号を更新する処理を抽出された全てのピークに対して順
次行って主相関配列の値を決定するスペクトル算出段
階、を実行させるためのプログラム。