JP4784016B2

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

Info

Publication number: JP4784016B2
Application number: JP2001244084A
Authority: JP
Inventors: 敏雄茂出木
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: 2011-09-28
Anticipated expiration: 2021-08-10
Also published as: JP2003058191A

Description

【０００１】
【発明の属する技術分野】
本発明は周波数解析方法および音響信号の符号化方法に関し、特に、生体信号の解析、医療診断・カルテ作成支援、医学教育等の医療関連、エンジン音の解析等の故障診断、自然音を音楽素材として取り込んだ音響コンテンツ制作支援への応用に適した技術に関する。
【０００２】
【従来の技術】
音響信号に代表される時系列信号には、その構成要素として複数の周期信号が含まれている。このため、与えられた時系列信号にどのような周期信号が含まれているかを解析する手法は、古くから知られている。例えば、フーリエ解析は、与えられた時系列信号に含まれる周波数成分を解析するための方法として広く利用されている。
【０００３】
このような時系列信号の解析方法を利用すれば、音響信号を符号化することも可能である。コンピュータの普及により、原音となるアナログ音響信号を所定のサンプリング周波数でサンプリングし、各サンプリング時の信号強度を量子化してデジタルデータとして取り込むことが容易にできるようになってきており、こうして取り込んだデジタルデータに対してフーリエ解析などの手法を適用し、原音信号に含まれていた周波数成分を抽出すれば、各周波数成分を示す符号によって原音信号の符号化が可能になる。
【０００４】
一方、電子楽器による楽器音を符号化しようという発想から生まれたＭＩＤＩ（Musical Instrument Digital Interface）規格も、パーソナルコンピュータの普及とともに盛んに利用されるようになってきている。このＭＩＤＩ規格による符号データ（以下、ＭＩＤＩデータという）は、基本的には、楽器のどの鍵盤キーを、どの程度の強さで弾いたか、という楽器演奏の操作を記述したデータであり、このＭＩＤＩデータ自身には、実際の音の波形は含まれていない。そのため、実際の音を再生する場合には、楽器音の波形を記憶したＭＩＤＩ音源が別途必要になるが、その符号化効率の高さが注目を集めており、ＭＩＤＩ規格による符号化および復号化の技術は、現在、パーソナルコンピュータを用いて楽器演奏、楽器練習、作曲などを行うソフトウェアに広く採り入れられている。
【０００５】
そこで、音響信号に代表される時系列信号に対して、所定の手法で解析を行うことにより、その構成要素となる周期信号を抽出し、抽出した周期信号をＭＩＤＩデータを用いて符号化しようとする提案がなされている。例えば、特開平１０−２４７０９９号公報、特開平１１−７３１９９号公報、特開平１１−７３２００号公報、特開平１１−９５７５３号公報、特開２０００−９９００９号公報、特開２０００−９９０９２号公報、特開２０００−９９０９３号公報、特開２０００−２６１３２２号公報、特開２００１−５４５０号公報、特開２００１−１４８６３３号公報には、任意の時系列信号について、構成要素となる周波数を解析し、その解析結果からＭＩＤＩデータを作成することができる種々の方法が提案されている。
【０００６】
【発明が解決しようとする課題】
上記各公報において提案してきたＭＩＤＩ符号化方式により、演奏録音等から得られる音響信号の効率的な符号化が可能になった。上記のように提案してきたＭＩＤＩ符号化方式は、大きく２通りに分けられる。１つは、特開平１０−２４７０９９号公報に代表されるような時系列信号の波形から強度値のピークを抽出すると共に抽出したピークの出現頻度に基づいて周波数を決定する手法である。もう１つは、特開２０００−２６１３２２号公報に代表されるような短時間フーリエ変換や一般化調和解析により時系列信号の解析を行なって、周波数成分を抽出する手法である。
【０００７】
前者は、所定の強度値以上をとるピークに基づいて抽出する周波数成分を決定していくため、時間分解能が高いと共に演算処理にかかる負荷が小さいという利点があるが、所定の区間において１つの周波数成分しか抽出できないという欠点がある。後者は、所定の区間において複数の周波数成分を抽出することができるが、周期関数との多数回の積和演算を必要とするため演算負荷が膨大であり、また、解析区間をオーバーラップさせながら演算するために充分な時間分解能が得られないという欠点がある。
【０００８】
特に、生体信号の一つである心音を解析しようとする場合、心音の波形成分は時間長が短いため、時間分解能が低い手法であるフーリエ変換や一般化調和解析を用いた手法は用いられることは少ない。この場合、信号強度のピークに基づいて周波数成分を抽出する手法が用いられることになるが、上述のように同時に１つの周波数成分しか抽出できないため、特徴的な１つの音しか判断することができず、隠れた疾患から生ずる僅かな音を検出することができないという問題がある。
【０００９】
上記のような点に鑑み、本発明は、所定の区間において複数の周波数成分を抽出することができ、かつ、演算負荷を抑え、時間分解能を高めることが可能であると共に、特に心音の解析を行なう場合に隠れた音成分も抽出することが可能な周波数解析方法および音響信号の符号化方法を提供することを課題とする。
【００１０】
【課題を解決するための手段】
上記課題を解決するため、本発明では、正負両極性の時系列信号から所定の主周波数の含有信号成分を分離するための周波数解析方法として、前記時系列信号に対して信号が０レベルと交差する複数の主ゼロ交差点を求めるための主ゼロ交差点検出段階と、前記主ゼロ交差点の中から、主ゼロ交差点ｔ１、ｔ２が、差ｔ２−ｔ１で定まる周期と主ゼロ交差点２つ分に相当する単位周期の関係が所定の条件を満たす範囲で、差ｔ２−ｔ１が最大となる場合に、前記主ゼロ交差点ｔ１を始点、前記主ゼロ交差点ｔ２を終点として選出することにより、全区間に渡って複数の単位区間を設定する単位区間設定段階と、前記各単位区間の始点となる主ゼロ交差点ｔ１と、終点となる主ゼロ交差点ｔ２の間にｋ個の主ゼロ交差点が存在するとき、（ｔ２−ｔ１）×２／ｋを主周期として算出し、前記単位区間における時系列信号の最大振幅を主振幅として求める主成分算出段階と、前記各単位区間の隣接主ゼロ交差点間における前記時系列信号の直流成分を除去して補正時系列信号を得る時系列信号補正段階と、前記補正時系列信号に対して信号が０レベルと交差する複数の副ゼロ交差点を求めるための副ゼロ交差点算出段階と、前記各単位区間における副ゼロ交差点の配置から副周期を算出し、当該単位区間における補正時系列信号から副振幅を求める副成分算出段階と、を有し、前記単位区間設定段階において設定された全ての単位区間について、前記主成分算出段階以降を繰り返すことにより、前記時系列信号における各単位区間に対応する主周期と主振幅および副周期と副振幅を算出するようにしたことを特徴とする。
【００１１】
本発明によれば、時系列信号に対して主ゼロ交差点を求め、この主ゼロ交差点を基に単位区間を設定し、各単位区間について主周期と主振幅を算出した後、時系列信号の補正を行って補正後の補正時系列信号から副周期と副振幅を算出するようにしたので、１つの区間において、複数の周波数成分を抽出することが可能となる。
【００１２】
【発明の実施の形態】
以下、本発明を図示する実施形態に基づいて説明する。
（１．本発明の基本原理）
はじめに、本発明に係る周波数解析方法および音声信号の符号化方法の基本原理を図１を参照しながら説明する。図１の上段に示すように、時系列の強度信号としてアナログ音響信号が与えられたものとする。図示の例では、横軸に時間軸ｔ、縦軸に信号強度Ａをとってこの音響信号を示している。本発明では、まずこのアナログ音響信号を、デジタルの音響信号として取り込む処理を行う。これは、従来の一般的なＰＣＭの手法を用い、所定のサンプリング周波数でこのアナログ音響信号をサンプリングし、信号強度Ａを所定の量子化ビット数を用いてデジタルデータに変換する処理を行えばよい。ここでは、説明の便宜上、ＰＣＭの手法でデジタル化した音響信号の波形も、図１の上段のアナログ音響信号と同一の波形で示すことにする。
【００１３】
次に、このデジタル音響信号の時間軸ｔ上に複数の単位区間を設定する。図示の例では、６つの単位区間Ｕ１〜Ｕ６が設定されている。第ｉ番目の単位区間Ｕｉは、時間軸ｔ上の始端ｓｉおよび終端ｅｉの座標値によって、その時間軸ｔ上での位置と長さとが示される。たとえば、単位区間Ｕ１は、始端ｓ１〜終端ｅ１までの（ｅ１−ｓ１）なる長さをもつ区間である。
【００１４】
こうして、複数の単位区間が設定されたら、個々の単位区間内の音響信号に基づいて、個々の単位区間を代表する所定の代表周波数および代表強度を定義する。ここでは、第ｉ番目の単位区間Ｕｉについて、代表周波数Ｆｉおよび代表強度Ａｉが定義された状態が示されている。たとえば、第１番目の単位区間Ｕ１については、代表周波数Ｆ１および代表強度Ａ１が定義されている。代表周波数Ｆ１は、始端ｓ１〜終端ｅ１までの区間に含まれている音響信号の周波数成分の代表値であり、代表強度Ａｉは、同じく始端ｓ１〜終端ｅ１までの区間に含まれている音響信号の信号強度の代表値である。単位区間Ｕ１内の音響信号に含まれる周波数成分は、通常、単一ではなく、信号強度も変動するのが一般的である。従来の手法とも共通する基本原理のポイントは、１つの単位区間について、１つもしくは複数の代表周波数と代表強度を定義し、これら代表値を用いて符号化を行う点にある。
【００１５】
すなわち、個々の単位区間について、それぞれ代表周波数および代表強度が定義されたら、時間軸ｔ上での個々の単位区間の始端位置および終端位置を示す情報と、定義された代表周波数および代表強度を示す情報と、により符号データを生成し、個々の単位区間の音響信号を個々の符号データによって表現するのである。単一の周波数をもち、単一の信号強度をもった音響信号が、所定の期間だけ持続する、という事象を符号化する手法として、ＭＩＤＩ規格に基づく符号化を利用することができる。ＭＩＤＩ規格による符号データ（ＭＩＤＩデータ）は、いわば音符によって音を表現したデータということができ、図１では、下段に示す音符によって、最終的に得られる符号データの概念を示している。
【００１６】
結局、各単位区間内の音響信号は、代表周波数Ｆ１に相当する音程情報（ＭＩＤＩ規格におけるノートナンバー）と、代表強度Ａ１に相当する強度情報（ＭＩＤＩ規格におけるベロシティー）と、単位区間の長さ（ｅ１−ｓ１）に相当する長さ情報（ＭＩＤＩ規格におけるデルタタイム）と、をもった符号データに変換されることになる。このようにして得られる符号データの情報量は、もとの音響信号のもつ情報量に比べて、著しく小さくなり、飛躍的な符号化効率が得られることになる。なお、代表周波数をノートナンバーｎの関数ｆ（ｎ）としたときに、代表周波数ｆ（ｎ）とノートナンバーｎ（０≦ｎ≦１２７）との関係は、以下に示す〔数式１〕で定義されることになる。
【００１７】
〔数式１〕
ｆ（ｎ）＝４４０×２^γ ⁽ⁿ⁾
γ（ｎ）＝（ｎ−６９）／１２
【００１８】
また、周波数と周期は互いに逆数の関係にあり、一方が定まれば他方も一義的に定まり、上記〔数式１〕により周波数が定まればノートナンバーも定まる。したがって、本発明においては、周波数・周期・ノートナンバーのいずれか１つの用語を用いている場合においても、他の２つについても同時に定まっているものとして説明していく。
【００１９】
（本発明に係る音響信号の符号化方法の具体的な手順）
上記基本原理をふまえて、以下、本発明に係る周波数解析方法を音響信号の符号化と共に説明する。図２は、本発明のより具体的な手順を示すフローチャートである。まず、符号化対象となる原音響信号をデジタルデータの形式で入力する（ステップＳ１）。
【００２０】
続いて、入力された原音響信号に対して信号から主ゼロ交差点を検出する（ステップＳ２）。主ゼロ交差点とは、音響信号が０レベルになる地点であるゼロ交差点のうち、補正処理が行われていない原音響信号に対して検出されるものを示す。ここで、図１に示した原音響信号の一部を時間軸方向に拡大したものを図３（ａ）に示す。ステップＳ２においては、信号が０レベルとなるゼロ交差点を検出する。このゼロ交差点は、時刻で表現されることになる。図３（ａ）の原音響信号から検出された主ゼロ交差点を図３（ｂ）に示す。
【００２１】
次に、ステップＳ２において検出された主ゼロ交差点にしたがって、単位区間を決定する（ステップＳ３）。単位区間とは、周波数成分の抽出および符号データを作成するための基本となる区間であり、図１の単位区間Ｕに相当する。単位区間は、始点となる主ゼロ交差点ｔ１から終点となる主ゼロ交差点ｔ２までの区間として設定される。始点となる主ゼロ交差点ｔ１としては、単位区間として未だ設定されていない区間において、時間的に先頭の主ゼロ交差点が選出される。終点となる主ゼロ交差点ｔ２は、以下のようにして選出される。
【００２２】
図３（ａ）に示すように原音響信号の波形は、始点となるゼロ交差点から２つ移動した位置のゼロ交差点で１周期となっている。そのため、この整数倍のセロ交差点を移動した位置を終点とするのが好ましい。また、安定した周波数値を得るため、終点をできる始点から離れた位置にとり、その平均周期をこの単位区間の周期とできるようにすることが好ましい。そのために、まず始点となる主ゼロ交差点ｔ１から２つ移動した主ゼロ交差点までの距離を単位周期Ｔ_Bと設定する。終点となる主ゼロ交差点ｔ２としては、始点となる主ゼロ交差点ｔ１から２の整数倍だけ離れた主ゼロ交差点のうち以下の〔数式２〕を満たす範囲で、ｔ２−ｔ１が最大となるものを選出するようにする。
【００２３】
〔数式２〕
｜Ｎ_A−Ｎ_B｜ ≦ １
ただし、Ｎ＝４０×ｌｏｇ₁₀（１／４４０Ｔ）＋６９
Ｔ＝（ｔ２−ｔ１）×２／ｋ
【００２４】
上記〔数式２〕において、Ｎ_Aはｔ２−ｔ１間の平均周期Ｔ_Aにより定まるノートナンバー、Ｎ_Bは単位周期Ｔ_Bにより定まるノートナンバーである。また、ｋは、始点となる主ゼロ交差点から終点となる主ゼロ交差点までに移動した主ゼロ交差点の数であり、２の整数倍となっている。なお、上記〔数式２〕において、第２式は、周期ＴからノートナンバーＮを求めるためのものであり、第３式は、平均周期Ｔ_Aを求めるためのものであり、ｋ＝２とした場合に単位周期Ｔ_Bが得られる。上記〔数式２〕においては、終点となる主ゼロ交差点ｔ２を、始点となる主ゼロ交差点ｔ１から、４個分、６個分、８個分…というように２の整数倍だけ離れた主ゼロ交差点に仮設定して判断していく。結局、ステップＳ２における処理は、単位周期Ｔ_Bで定まる周波数と、平均周期Ｔ_Aで定まる周波数をそれぞれノートナンバーに変換した際に、半音の差（ノートナンバーとしては±１の差）以内に収まるような主ゼロ交差点ｔ２を終点として選出する処理を行うことになる。単位区間の設定は音響信号の全区間に渡って行われる。後続する単位区間においては、先行する単位区間の終点となる主ゼロ交差点の次の主ゼロ交差点を、始点となる主ゼロ交差点として、上記説明と同様に処理を行っていく。
【００２５】
全単位区間の設定が行われたら、各単位区間について、主周波数成分の算出を行う（ステップＳ４）。周波数成分とは、周期（周波数）および振幅を意味するものである。周期としては、単位区間における平均周期Ｔ_Aが設定され、振幅としては、単位区間内において、信号の絶対値が最大となる値Ａが設定される。なお、主周波数成分とは、周波数成分のうち、補正が行われていない原音響信号から抽出されたものをいい、ステップＳ７で算出される副周波数成分と区別している。
【００２６】
上記ステップＳ３およびステップＳ４による処理、すなわち、ゼロ交差点を用いて単位区間を決定し、単位区間に対応する周波数成分を算出する処理については、従来から行われているものと同様である。ただし、これでは、上記発明が解決しようとする課題の段落で説明したように、１つの単位区間について１つの周波数成分しか抽出できない。本発明では、１つの単位区間から複数の周波数成分を抽出するために後述するようなステップＳ５以下の処理を行っている。
【００２７】
ステップＳ４において主周波数成分が算出されたら、次にその単位区間における音響信号から直流成分の除去を行う（ステップＳ５）。ここでの、直流成分の除去は単位区間全体に渡って一律に行うのではなく、音響信号が０以上になっている部分については、当該部分の平均強度を各時刻の信号強度から減じるようにし、音響信号が０以下になっている部分についても、当該部分の平均強度を正負の符号をそのままにして信号強度から減じる（実際には平均強度は増加していることになる）ようにする。これにより、振幅の大きな部分が抑えられた補正音響信号が得られることになる。図３（ｃ）に示した原音響信号を補正した補正音響信号を図３（ｄ）に示す。
【００２８】
続いて、補正音響信号に対して、信号が０レベルと交差する副ゼロ交差点を検出する（ステップＳ６）。これは、上記ステップＳ２において、原音響信号に対して主ゼロ交差点を検出したのと全く同様な処理で行われる。ただし、図３（ｄ）に示したように、補正音響信号では、原音響信号に比べて明らかにゼロ交差点が増えているので、多くの副ゼロ交差点が検出されることになる。図３（ｄ）に示した補正音響信号から検出された副ゼロ交差点を図３（ｅ）に示す。
【００２９】
続いて、補正音響信号および検出された副ゼロ交差点に従って、副周波数成分の算出を行う（ステップＳ７）。副周波数成分のうち、周期としては、単位区間における始点となる副ゼロ交差点ｔ３と終点となる副ゼロ交差点ｔ４で定まる平均周期が与えられる。また、副周波数成分のうち、振幅としては、単位区間内において、補正音響信号の絶対値が最大となる値が設定される。
【００３０】
上記ステップＳ５〜ステップＳ７の処理は、所定の回数だけ繰り返され、繰り返された回数分の副周波数成分が抽出されることになる。
【００３１】
全ての単位区間について、主周波数成分、副周波数成分の抽出が行われたら、抽出された周波数成分を基に符号化を行う（ステップＳ８）。符号化の段階では、特に主周波数成分と副周波数成分は区別せず、１つの単位区間に複数の周波数成分があるものとして扱われる。こうして、各単位区間について、所定数の周波数成分が得られたら、各周波数成分の周期（周波数）に対応する「音の高さを示す情報」、各周波数成分の振幅に対応する「音の強さを示す情報」、当該単位区間の始点（主ゼロ交差点ｔ１）に対応する「音の発音開始時刻を示す情報」、当該単位区間の終点（主ゼロ交差点ｔ２）に対応する「音の発音終了時刻を示す情報」、の４つの情報を含む所定数の符号データを作成すれば、当該単位区間内の区間信号を所定数の符号データにより符号化することができる。符号データとして、ＭＩＤＩデータを作成するのであれば、「音の高さを示す情報」としてノートナンバーを用い、「音の強さを示す情報」としてベロシティーを用い、「音の発音開始時刻を示す情報」としてノートオン時刻を用い、「音の発音終了時刻を示す情報」としてノートオフ時刻を用いるようにすれば良い。
【００３２】
以上、本発明の好適な実施形態について説明したが、上記周波数解析方法および音響信号の符号化方法は、コンピュータ等で実行されることは当然である。具体的には、図１のフローチャートに示したようなステップを上記手順で実行するためのプログラムをコンピュータに搭載しておく。そして、音響信号をＰＣＭ方式等でデジタル化した後、コンピュータに取り込み、ステップＳ１〜ステップＳ８の処理を行った後、周波数解析を行った場合はその周波数成分を抽出し、符号化まで行った場合はＭＩＤＩ形式等の符号データをコンピュータより出力する。符号化を行った場合は、出力された符号データは、例えば、ＭＩＤＩデータの場合、ＭＩＤＩシーケンサ、ＭＩＤＩ音源を用いて音声として再生される。
【００３３】
【発明の効果】
以上、説明したように本発明によれば、時系列信号に対して信号が０レベルと交差する複数の主ゼロ交差点を求めるための主ゼロ交差点検出段階と、前記主ゼロ交差点の中から始点となる主ゼロ交差点と終点となる主ゼロ交差点を選出することにより、全区間に渡って複数の単位区間を設定する単位区間設定段階と、前記各単位区間における主ゼロ交差点の配置から主周期を算出し、前記単位区間における時系列信号から主振幅を求める主成分算出段階と、前記各単位区間における隣接主ゼロ交差点間に対応する前記時系列信号の直流成分を除去して補正時系列信号を得る時系列信号補正段階と、前記補正時系列信号に対して信号が０レベルと交差する複数の副ゼロ交差点を求めるための副ゼロ交差点算出段階と、前記各単位区間における副ゼロ交差点の配置から副周期を算出し、当該単位区間における補正時系列信号から副振幅を求める副成分算出段階とを有し、全単位区間について、前記主成分算出段階以降を繰り返すことにより、前記時系列信号における各単位区間に対応する主周期と主振幅および副周期と副振幅を算出するようにしたので、１つの区間において、複数の周波数成分を抽出することが可能となるという効果を奏する。特に生体信号の１つである心音の解析に用いた場合に隠れた疾患に基づく音成分も抽出することが可能となるという効果も奏する。
【図面の簡単な説明】
【図１】本発明に係る周波数解析方法および音響信号の符号化方法の基本原理を示す図である。
【図２】本発明に係る周波数解析方法および音響信号の符号化方法のフローチャートである。
【図３】音響信号における単位区間の設定および周期・振幅の決定を説明するための図である。
【符号の説明】
Ｔ_A ・・・平均周期
Ｔ_B ・・・単位周期
ｔ１・・・始点となる主ゼロ交差点
ｔ２・・・終点となる主ゼロ交差点
ｔ３・・・始点となる副ゼロ交差点
ｔ４・・・終点となる副ゼロ交差点

Claims

正負両極性の時系列信号から所定の主周波数の含有信号成分を分離するための周波数解析方法であって、
前記時系列信号に対して信号が０レベルと交差する複数の主ゼロ交差点を求めるための主ゼロ交差点検出段階と、
前記主ゼロ交差点の中から、主ゼロ交差点ｔ１、ｔ２が、差ｔ２−ｔ１で定まる周期と主ゼロ交差点２つ分に相当する単位周期の関係が所定の条件を満たす範囲で、差ｔ２−ｔ１が最大となる場合に、前記主ゼロ交差点ｔ１を始点、前記主ゼロ交差点ｔ２を終点として選出することにより、全区間に渡って複数の単位区間を設定する単位区間設定段階と、
前記各単位区間の始点となる主ゼロ交差点ｔ１と、終点となる主ゼロ交差点ｔ２の間にｋ個の主ゼロ交差点が存在するとき、（ｔ２−ｔ１）×２／ｋを主周期として算出し、前記単位区間における時系列信号の最大振幅を主振幅として求める主成分算出段階と、
前記各単位区間の隣接主ゼロ交差点間における前記時系列信号の直流成分を除去して補正時系列信号を得る時系列信号補正段階と、
前記補正時系列信号に対して信号が０レベルと交差する複数の副ゼロ交差点を求めるための副ゼロ交差点算出段階と、
前記各単位区間における副ゼロ交差点の配置から副周期を算出し、当該単位区間における補正時系列信号から副振幅を求める副成分算出段階と、を有し、
前記単位区間設定段階において設定された全ての単位区間について、前記主成分算出段階以降を繰り返すことにより、前記時系列信号における各単位区間に対応する主周期と主振幅および副周期と副振幅を算出するようにしていることを特徴とする周波数解析方法。
前記副成分算出段階が、前記単位区間内の最初の副ゼロ交差点ｔ３と、末端となる副ゼロ交差点ｔ４の間にｋ個の副ゼロ交差点が存在するとき、前記副周期として（ｔ４−ｔ３）×２／ｋを与え、副振幅として前記単位区間における補正時系列信号の最大振幅を与えることを特徴とする請求項１に記載の周波数解析方法。
前記時系列信号補正段階から副成分算出段階までを所定の回数繰り返し、１つの単位区間に対して複数の副周期と副振幅を算出するようにすることを特徴とする請求項１または請求項２に記載の周波数解析方法。
時系列信号として与えられた音響信号を符号化する方法であって、
前記音響信号に対して信号が０レベルと交差する複数の主ゼロ交差点を求めるための主ゼロ交差点検出段階と、
前記主ゼロ交差点の中から、主ゼロ交差点ｔ１、ｔ２が、差ｔ２−ｔ１で定まる周期と主ゼロ交差点２つ分に相当する単位周期の関係が所定の条件を満たす範囲で、差ｔ２−ｔ１が最大となる場合に、前記主ゼロ交差点ｔ１を始点、前記主ゼロ交差点ｔ２を終点として選出して単位区間を設定する単位区間設定段階と、
前記各単位区間の始点となる主ゼロ交差点ｔ１と、終点となる主ゼロ交差点ｔ２の間にｋ個の主ゼロ交差点が存在するとき、（ｔ２−ｔ１）×２／ｋを主周期として算出し、前記単位区間における音響信号の最大振幅を主振幅として求める主成分算出段階と、
前記単位区間の隣接主ゼロ交差点間における前記音響信号の直流成分を除去し、補正音響信号を得るための音響信号補正段階と、
前記補正音響信号に対して信号が０レベルと交差する複数の副ゼロ交差点を求めるための副ゼロ交差点算出段階と、
前記単位区間内の最初の副ゼロ交差点ｔ３と、末端となる副ゼロ交差点ｔ４の間にｋ個の主ゼロ交差点が存在するとき、（ｔ４−ｔ３）×２／ｋを副周期として算出し、当該単位区間における補正音響信号の最大振幅を副振幅として求める副成分算出段階と、を有し、
前記単位区間設定段階において設定された全ての単位区間について、前記主成分算出段階以降を繰り返すことにより、前記音響信号における各単位区間に対応する主周期と主振幅および副周期と副振幅を算出し、
前記算出された主周期または副周期に対応する音の高さ情報と、主振幅または副振幅に対応する音の強さ情報と、各単位区間の始点となる主ゼロ交差点に対応する発音開始時刻と、各単位区間の終点となる主ゼロ交差点に対応する発音終了時刻とからなる情報に変換を施して符号データを生成することを特徴とする音響信号の符号化方法。
前記音響信号補正段階から副成分算出段階までを所定の回数繰り返し、１つの単位区間に対して複数の副周期と副振幅を算出するようにすることを特徴とする請求項４に記載の音響信号の符号化方法。
コンピュータに、
時系列信号に対して信号が０レベルと交差する複数の主ゼロ交差点を求めるための主ゼロ交差点検出段階、
前記主ゼロ交差点の中から、主ゼロ交差点ｔ１、ｔ２が、差ｔ２−ｔ１で定まる周期と主ゼロ交差点２つ分に相当する単位周期の関係が所定の条件を満たす範囲で、差ｔ２−ｔ１が最大となる場合に、前記主ゼロ交差点ｔ１を始点、前記主ゼロ交差点ｔ２を終点として選出することにより、全区間に渡って複数の単位区間を設定する単位区間設定段階、
前記各単位区間の始点となる主ゼロ交差点ｔ１と、終点となる主ゼロ交差点ｔ２の間にｋ個の主ゼロ交差点が存在するとき、（ｔ２−ｔ１）×２／ｋを主周期として算出し、前記単位区間における時系列信号の最大振幅を主振幅として求める主成分算出段階、
前記各単位区間の隣接主ゼロ交差点間における前記時系列信号の直流成分を除去して補正時系列信号を得る時系列信号補正段階、
前記補正時系列信号に対して信号が０レベルと交差する複数の副ゼロ交差点を求めるための副ゼロ交差点算出段階、
前記各単位区間内の最初の副ゼロ交差点ｔ３と、末端となる副ゼロ交差点ｔ４の間にｋ個の主ゼロ交差点が存在するとき、（ｔ４−ｔ３）×２／ｋを副周期として算出し、当該単位区間における補正時系列信号の最大振幅を副振幅として求める副成分算出段階、
を実行させ、
前記単位区間設定段階において設定された全ての単位区間について、前記主成分算出段階以降を繰り返し実行させることにより、前記時系列信号における各単位区間に対応する主周期と主振幅および副周期と副振幅を算出するためのプログラム。
前記時系列信号補正段階から副成分算出段階までを所定の回数繰り返し、１つの単位区間に対して複数の副周期と副振幅を算出するようにすることを特徴とする請求項６に記載のプログラム。