JPH02293800A - ピツチ関連遅延値を導出する方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 Ａ．産業上の利用分野 本発明は音声信号を効率的にコーディングする方法に関
する。

Ｂ．従来の技術及びその課題 音声信号の効率的なコーディングという場合、これは音
声信号の高品質のデジタルエンコーディングを得るとい
うことを意味するだけでなく、コスト及びコーダの複雑
さの最適化をも意味している。

既に知られたコーダのうちあるものは、もとの音声信号
を処理し、そこから音声代表残差信号を導出し、遅延装
置を調整するのに使用される検出されたピッチ関連デー
タで調節されるＬＴＰ（Ｌｏｎｇ−Ｔｅｒｍ　Ｐｒｅｄ
ｉｃｔｉｏｎ）手段を用いて残差予測信号を計算し、現
在の残差と予測された残差を組合せて残差誤差信号を生
成し、最後にこの残差誤差信号を低ビットレートでコー
ド化する。、上述のタイプのコーディング手法の効率に
対する大きな改善は欧州特許出願第ＥＰ８７４３０００
６．４号（特開昭６３−２２３７９９号）において示さ
れるように，ピッチ又はそのピッチの調波（以下、単に
ピッチ、ピッチ代表データ又はピッチ関連データともい
う）を次のようにして検出することによって達成された
。すなわち、この方法は、第１のステップでゼロクロス
及びピークの検出による粗いピッチの測定を行い、第２
のステップで検出されたピッチのピークについて遂行さ
れる自己相関に基づいてピッチを徴調するものである。

上記のピッチ検出方法は特に有益なものであるけれども
、コーダの全体的な計算負荷に比してかなり大きい計算
負荷をもたらす。したがって本発明の目的は従来に比べ
て計算負荷の少ないピンチ検出の技術を提供することに
ある。

Ｃ．課題を解決するための手段 この目的を達成するため、音声から導出されたデジタル
信号ｒ　（ｎ）を低ビットレー１・の信号に変換する長
区間予測に基づく音声コーダで使用される長区間予測フ
ィルタを調整するためのピッチ関連遅延値Ｍを導出する
本発明の方法は下記の（ａ）ないし（ｈ）のステップを
有することを特徴としている。

（ａ）上記信号ｒ　（ｎ）を各々がＮ個のサンプルから
成るセグメントに分割するステップ。

（ｂ）各セグメントをｊ個（ｊは所定の整数）のサ一３ ブセグメントに分割するステップ。

（ｃ）現信号の最初のサブセグメントと復元された信号
のセグメントとの相互相関をとって下記の相ただしｋ’
　＝Ｎ／ｊ （ｄ）ピーク位置Ｒ（Ｍｌ）についてＲ　（ｎ）の値を
選出し、上記フィルタの遅延値をＭ１に設定し、信号を
サブセグメント１個分についてシフトするステップ。

（ｅ）Ｍｌの高調波及び低調波のところに位置する所定
数のサンプルについてのサンプル指標ｎを計算するステ
ップ。

（ｆ）ステップ（ｅ）で定義されたｎについて相互相関
関数の値Ｒ　（ｎ）を計算するステップ。

（ｇ）新しい遅延値Ｍ２を導出するためピーク位置につ
いてＲ　（ｎ）の値を選出するステップ。

（ｈ）上記（ｅ）から（ｇ）までのステップをＭ１の代
わりにＭ２について実行し、以下同様にＭｊになるまで
（ｅ）から（ｇ）までのステップを繰返すステップ。

なお、ステップ（ｅ）において、ｐ及びｋを所定の整数
としてｎ　＝　ｐ　Ｍ　１　＋　ｋであるようなＭｌ／
ｐ，・・・・・・、Ｍ１／３、Ｍ１／２、Ｍ１、２Ｍ１
、３Ｍ１、・・・・・・、ｐＭ１となる所定数のサンプ
ルについてｎを計算する。

以下、本発明の作用を実施例とともに説明する。

Ｄ．実施例 はじめに本実施例を概説する。本実施例は長区間予測に
基づく音声コーダにおいて遅延データとして用いられる
ピッチ関連データを最小の計算負荷で高速に導出する方
法を提供するものである。

これは、処理すべき信号を各々Ｎ個のサンプルから成る
連続的なセグメントに分割し、各６セグメントをｊ個の
サブセグメントに分割し、現在の最初のサブセグメント
のサンプルと既にデコードされたセグメントとの相互相
関をとって相互相関関数を導出するとともに第１の遅延
Ｍ１として用いられる相互相関のピーク位置の指標を導
出し、音声コーダのＬＴＰループに対してＭ１を設定し
、上記第１の遅延Ｍ］の高調波及び低調波のところのサ
ンプル指標を計算し、上記指標されたサンプルについて
新しい相互相関関数を計算して新しい遅延データＭ２を
導出し、これを最後のサブセグメン１−まで続行し、さ
らに次の信号セグメントについてこのプロセスを繰返す
ことによって達成される。以下、図面を参照しながら本
実施例を詳細に説明する。

第１図は本発明を適用することのできる音声コーダを示
す図である。もとの音声信号ｓ　（ｎ）はＡ／Ｄコンバ
ータ（図示せず）においてナイキスト周波数でサンプリ
ングされ１サンプルにっき１２ビッ１〜でＰＣＭコード
化される。そのようなコーダ（ＲＰＥ／ＬＴＰ）は中程
度のビットレ−１・で市外品質の音声コード化圧縮を実
現できるけれども、圧縮すべき信号が連続的成分を呈す
る場合には可聴周波のノイズトーンが生成されることが
ある。この場合、これはＡ／Ｄコンバータの使用に起因
するものである。ＲＰＥ／ＬＴＰコーダ又はデコーダに
おいては、高周波成分が生成されることが必要であり、
これはベースバンドフォールディング（ｂａｓｅｂａｎ
ｄ　ｆｏｌｄｉｎｇ）によって達成される。その結果、
音声信号が高レベルのオフセッ１・を含んでいる場合は
、ベースバンド信号もまたこのオフセットを含み、余分
に復元された信号がミラー周波数で純音を呈することに
なる。オフセットの検出はＣ　Ｅ　Ｐ　Ｔ　（Ｅｕｒｏ
ｐｅａｎ　Ｃｏｍｍｉｓｓｉｏｎ　ｆｏｒＰｏｓｔ　ａ
ｎｄ　Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ）のＧＳＭ
Ｏ６．１０で定義されるようなノッチハイパスフィルタ
を使用する装置９で行われる。

要約するに、ＤＣ成分を除去するために作られたこのフ
ィルタは固定係数の再帰的デジタルフィルタから成る（
この係数はＣＥＰＴによって無線電話について定義され
ているものである）。

オフセット検出のためのより簡単な別のアルゴリズムは
以下のようにＬＴＰループにおいて（すなわち、装置２
２の出力に関して）実施するものである。

デコードされた信号のＤＣ成分は残差誤差信号ｅＪ（ｎ
）から除去されてオフセットのない新しい信号ｅ’　（
ｎ）が以下の計算を実行することにより得られる。

ｘｏｆ（１）＝ｘ’Ｌ（］．）−Σ（ｘ’Ｌ（ｉ））／
Ｃ　　１＝１、Ｃｊ＝１ ここでｘ　’　Ｌ　（　１　）はＲＩ’Ｅで選択した遅
延Ｌについてデコードされたパルスの振幅を表わし、Ｃ
はこれらのパルスの数を表わす。

次に、信号Ｘ。ｆ（ｎ）がゼロの値のサンプルをインタ
ーリーブすることによってオーバーサンプリングされ、
オフセットのない全帯域信号ｅ’（ｎ）が生成される。

受信側では、同様なオペレーションがデコードされたベ
ースバンド信号について行われる。

再び第１図の説明に戻る。装置９によって供給される事
前処理された信号が短区間予測フィルタ１０に送られる
。

短区間予測フィルタはラティスデジタルフィル゛タから
成り、このフィルタのタップ係数はＬＰＧ分析によって
装置１１において信号から動的に導出される。この目的
のため、事前処理される信号は各セグメントがオーバー
ラップしない１６０個のサンプルに分けられる（各々は
２０ミリ秒の信号である）。ＬＰＣ分析はＳｈｕｒ再帰
アルゴリズムを用いて８つの反射係数を計算することに
よって各セグメン１へについて行われる。Ｓｈｕｒアル
ゴリズムの詳細が必要であれば、前掲のＧＳＭＯ６．１
０の仕様を参照されたい。

反射係数はＬ　Ａ　Ｒ　（ｌｏｇ　ａｒｅａ　ｒａｔｉ
ｏ）係数に変換される。これは区分的かつ線形的に３２
ビッ１・（６、５、５、４、３、３、３、３）で量子化
され、ｓ　（ｎ）の再合成の間に使用するためにコード
化される。

短区間分析フィルタの８つの係数は次のように処理され
る。まず、量子化されコード化されたＬＡＲ係数がデコ
ードされる。次に、直前のＴ−　Ａ　Ｒ係数のセツ１〜
が５ミリ秒の遷移期間内で線形的に補間され、スプリア
スな過渡状態が防止される。

最後に、補間されたＬＡＲがラティスフィルタの反射係
数に再変換される。このフィルタは音声から導出された
（すなわち、残差）信号ｒ（ｎ）の１６０個のサンプル
を生成する。信号ｒ　（ｎ）はピッチ関連周波数で一定
の冗長性を有する比較的平坦な周波数スペクトラムを示
す。

装置１２はこの残差信号を処理してそこからピッチ、調
波を表わすデータ、換言すればピッチ関連情報Ｍ及びゲ
インパラメータｂを導出する。Ｍ及びｂは次の式で表わ
されるような２領域における演算を実行する長区間予測
フィルタ１４を調整するのに使用される。

一阿 Ｒ”（ｚ）＝ｂ．ｚ　　Ｒ’　（ｚ）　　　　　（１）
ここで、Ｒ’（ｚ）及びＲ　ＩＩ　（　ｚ）はそれぞれ
時間領域信号ｒ’　（ｎ）及びｒｕ（ｎ）の２領域変換
を表わす。

したがって式（１）の演算を実行する装置は動的にＭに
調整すべき長さ（ピッチ又は調波に関連する遅延データ
）を有する遅延線と、ゲイン装置とを具備しなければな
らない。これについては後でさらに詳しく説明する。

ｂ及びＭを効率的に測定することはコーダにとっては重
要な事項である。というのは、Ｍで調整された長区間予
測フィルタからの予測残差信号出力ｒ″′（ｎ）は長区
間相関解除された予測誤差信号ｅ　（ｎ）を導出するた
めに残差信号から減ずる必要があり、．（ｎ）はＲ　Ｐ
　Ｅ　（Ｒｅｇｕｌａｒ　ＰｕｌｓｅＥｘｃｉｔａｔｉ
ｏｎ　）手法を用いてパルスｘ　（ｎ）のシーケンスに
コード化されるからである。換言すれば、ＲＰＥ装置１
６を使用して、たとえば連続的なＰＣＭコード化された
ｅ　（ｎ）のサンプルのサブセグメントをＡＰＣＭ量子
化器２０で次に量子化されるより少ない数（たとえば１
５未満）の最も有意なパルスに変換する。以上の考察に
より、フィルタ１４の正確な調整したがってｂ及びＭの
良好な算定の重要性が理解されよう。

簡単にいうと、ＲＰＥ手法を用いる場合、４０個のｅ　
（ｎ）のサンプルから成る各サブグループはインターリ
ーブされたシーケンスに分けられる。

たとえば、２つの１３サンプルと１つの１４サンプルが
インターリーブされる。ＲＰＥ１６はこの３つのインタ
ーリーブされたシーケンスの中から１つを選択する。こ
れは、もとのシーケンスと比較された場合に最小２乗誤
差を与えるものである。

２ビット（Ｌ）で選択されたシーケンスを識別すること
はデータシーケンスｘ　（ｎ）の適切な同期Ｌ に役立つ。

ＲＰＥコーディングの動作についてさらに詳細が必要で
あれば、Ｐ．　Ｋｒｏｏｎらによる”Ｒｅｇｕ’ｌａｒ
Ｐｕｌｓｅ　Ｅｘｃｉｔａｔｉｏｎ，　ａ　Ｎｏｖｅｌ
　Ａｐｐｒｏａｃｈ　ｔｏＥｆｆｅｃｔｉｖｅ　ａｎｄ
　Ｅｆｆｉｃｉｅｎｔ　Ｍｕｌｔｉｐｌｅ　Ｃｏｄｉｎ
ｇ　ａＳｐｅｅｃｈ”　（ＩＥＥＨ　Ｔｒａｎｓａｃｔ
ｉｏｎｓ　ａｎｄ　ＡｃｏｕｓｔｉｃｓＳｐｅｅｃｈ　
ａｎｄ　Ｓｉｇｎａｌ　Ｐｒｏｃｅｓｓｉｎｇ　ＡＳＳ
Ｐ３４　　第５号１９８６年１０月）を参照されたい。

ＲＰＥによる長区間予測を使うと，特に次の場合に全体
的ビットレート対品質パラメータを最適化することがで
きる。すなわち、長区間予測フイルタ１４にｒ　（ｎ）
にできるだけ近いパルス列ｒ′（ｎ）を供給した場合、
換言すれば，装置１６及び量子化器２０によって供給さ
れるコード化ノイズ及び量子化ノイズが補償された場合
である。このため、デコーディングオペレーションは装
置２２で行われる。装置２２の出力ｅ’　（ｎ）に予測
一１２ 残差ｒ″′（ｎ）が加えられて復元された残差ｒ′（ｎ
）が生成される。さらに、ＲＰＥコーダのまわりの閉ル
ープは最小の範囲をピッチ関連データ検出ウィンドウに
設定することによって実時間で動作することができる。

第１図の長区間予測フィルタ１４のインプリメンテーシ
ミンが第２図に示されている。復元された残差信号は１
２０個分のサンプルの長さ（Ｍについての最大値は１２
０）の遅延線（又はシフ１〜レジスタ）に供給される。

この遅延線の出方はＬＴＰ係数計算手段１２に与えられ
、後の処理のため，係数ｂ及びＭが導出される。遅延線
のタップは前に計算されたＭの値に調整される。ゲイン
係数ｂは、ｒ（ｎ）からｒＩＩ（ｎ）を減じてｅ　（ｎ
）を生成する前に上記タップで利用可能なデータに適用
される。

長区間予測残差信号をこうして残差信号から減じて誤差
信号ｅ（ｎ）を導出する。誤差信号ｅ（ｎ）は量子化器
２０による量子化の前にＲＰＥ１６によってコード化さ
れる。

このコーダのアーキテクチャの重要な利点はＭが装置１
２で正確に測定されている限りそれはＳ（ｎ）のピッチ
又はピッチの調波を表わす遅延となるはずであるという
事実から導かれる。

このため、遅延Ｍは５ミリ秒（サンプル４０個分）ごと
に計算される。信号ｒ　（ｎ）は１６０個分のサンプル
の長さの連続的なセグメントに分けられる。なお、各セ
グメントはｊ（たとえばｊ＝４）個のサブセグメントに
さらに分けられる。

ｒ　（ｎ）のサンプルのうちの最初のサブセグメント及
び前に復元された励起（　ｅｘｃｉｔａｔｉｏｎ　）セ
グメントｙ　（ｎ）は以下のようにして相互相関が計算
される。

計算されたＲ　（ｎ）の値はピークの場所に関して分類
（選出）され以下の式によって第１の最適遅延値Ｍ１が
導出される。

Ｒ　（Ｍｌ）　＝Ｍａｘ　（Ｒ　（ｎ））　　ｎ　＝４
０，１２０　　（３）対応するゲインの値ｂ１は以下の
式から導出される。

ＬＴＰフィルタはｂ１及びＭ１で調整され、その信号が
１サブセグメント（４０個分のサンプル）にわたっーで
シフトされる。

次のサブセグメン１・の場合、ピッチ関連遅延値は次の
ようにして算定される。

ｐが所定の整数値（たとえばｐ＝３）であるようなＭ１
、２Ｍ１、３Ｍ１、・・・・・・、ｐＭｌ、Ｍ１／２、
Ｍ１／３、・・・・・・、Ｍｌ／ｐを導出するため第１
のＭ１の倍数及び約数が計算される。次にｋを所定の整
数値（たとえばｋ＝５）として以下のようなサンプルの
指標ｎが定義される。

ｎ　＝　（Ｍｌ−ｋ）、　（Ｍｌ−ｋ−１）、　・・・
・・　　、　（Ｍ１）、　（Ｍｌ＋ｋ−１）、　（Ｍｌ
十ｋ）ｎ　＝　（２Ｍ１−ｋ）、（２Ｍ１−ｋ−１）、
・・・・・・、（２Ｍ１）、（２Ｍ１４ｋ−１）、（２
Ｍ１＋ｋ）ｎ　＝　（ｐＭ１−ｋ）、（ｐＭ１−ｋ−１
）、・曲・、（ｐＭ１）、・・曲、（ｐＭ１＋ｋ−１）
、（ｐＭ１＋ｋ） ｎ　＝　（（Ｍｌ／２）−ｋ）、（（Ｍｌ／２）−ｋ−
１）、・・・・・・、（Ｍｌ／２）、・・曲（（Ｍｌ／
２）＋ｋ−１）、　（（Ｍｌ／２）＋ｋ）ｎ　＝　（（
Ｍｌ／３）−ｋ）、（（Ｍｌ／３）−ｋ−１）、・・・
・・・、（Ｍｌ／３）、・・・・・・（（Ｍｌ／３）＋
ｋ−１）、（（Ｍｌ／３）＋ｋ）ｎ＝（（阿１／ｐ）−
ｋ）、　（（Ｍｌ／ｐ）−ｋ−１）、　・・・・・・、
　（Ｍｌ／ｐ）、　・・・・・・（（Ｍｌ／ｐ）＋ｋ−
１）、（（Ｍｌ／ｐ）＋ｋ）ただし、３　９　＜　ｎ　
＜　１　２　１換言すれば、上記の計算されたｎの値は
Ｍ１の倍数及び約数に選択すべきピッチ関連値に関して
位置するサンプルについてのサンプル指標である。

相互相関関数（２）は上記のように定義され指標された
サンプルについて計算され、そのようにして計算された
Ｒ　（ｎ）の値は再びピークの場所に関して分類（選出
）される。これにより、第２のサブセグメントについて
の新しい最適の遅延Ｍ２が導出される。

同じアルゴリズムがＭ１の代わりにＭ２に対して繰返さ
れ、次の遅延Ｍ３が計算される。これは直前の現サブセ
グメントＭｊに到達するまで続けられる。そして全体的
なプロセスを次のサンプルセグメントに対して繰返すこ
とができる。

各々のＭの値について、式（４）に基づき対応するゲイ
ンｂが計算される。これらのＬＴＰパラメータは２ビッ
ト及び７ビットでそれぞれコード化することができる。

第３Ａ図、第３Ｂ図及び第４図は以上の高速ピッチ検出
プロセスのアルゴリズムを示す図である。

このアルゴリズムはマイクロプロセッサで走行するプロ
グラムに容易に変換することができる。この例はｊ＝４
個のサブセグメントに分けられる１６０個分のサンプル
の長さを有するセグメン１・の処理に関するものである
。音声コーディング分析の場合、ｓ　（ｎ）のフローは
１６０個分のサンプルの長さのセグメントに分けられ、
まず、オフセット検出処理及び１６０個の１１　８．　
ＩＩサンプルの生成が行われる。次にサンプル“Ｓｏ′
はＬＡＲデータに量子化される８つのＰＡＲＣＯＲ係数
ｋｉを生成するＬＰＧ分析に付される。

ＰＡＲＣＯＲ係数ｋｉを用いてＬＰＧ短区間フィルタを
調整し、１６０個のサンプルＩｔ　ｓｏＩ＋を処理して
残差信号ｒ　（ｎ）を導出する。上記ｒ（ｎ）のサンプ
ルセグメントは４０個分のサンプルの長さのサブセグメ
ントに分けられる。これらのサブセグメントの各々は前
に導出された１２０個分のサンプルの長さのｙセグメン
トを用いるＬＴＰ係数の計算のために処理される。ＬＴ
Ｐ係数の計算によってサブセグメントの合成のために量
子化されたｂ及びＭが与えられる。これらのｂ及びＭの
データはＬＴＰフィルタの調整のために用いられる量子
化（Ｑ）の前に一度量子化解除（Ｑ）される（又は直接
的に選択される）。次に、ｒ（ｎ）から上記ＬＴＰフィ
ルタの出力を減じると、ｅ　（ｎ）が得られる。

４０個の連続的なｅ（ｎ）のサンプルはより少ないＸＬ
個のサンプルのセット及びセット参照しにコード化され
る。Ｘ及びＬはそれぞれ量子化Ｌ （Ｑ）される。量子化解除（Ｑ）されオーバーサンプリ
ングされたサンプルのサブセグメンｌ−（ｅ’（ｎ））
はＬＴＰ係数の計算から始まるオペレーションを繰返す
ことによって完全なセグメン１−までのＬＴＰの合成及
び遅延線の更新のために使用される。

相関的な音声合成（すなわち、デコーディング）は以下
のオペレーションを有する。

・１６０個のｅ′のサンプルを生成するための量子化解
除されたパラメータＸＬ及びＬを用いるＲＰＥデコーデ
ィング ・量子化解除されたＬＴＰパラメータを用いるＬＴＰ合
成及び遅延線の更新、１６０個の復元された残差サンプ
ルｒ″の導出 ・合成された残差信号サンプルについてのＬＰＧ合成及
び合成された音声信号Ｓ′の生成ＬＴＰ係数の計算ステ
ップをさらに詳しく説明する（第４図）。Ｍ１のＨ１算
のためにバッファされた最初の入力サンプルは現在のｙ
の１２０個のサン−プル（Ｏないし１１９）とｒの４０
個のサンプル（０ないし３９）である。次にＲ　（ｎ）
の値が式（３）に従って分類され、Ｍ１が導出される。

このＭｌｋ用いて式（４）に従ってｂ１が計算され、こ
れに応じてＬＴＰフィルタがセットされ、信号が１サブ
セグメント（すなわち、サンプル４０個分）だけシフト
される。次に以下の式に従ってサンプル指標をセットす
ることによってＭ２が計算される。

ｎ＝ｐＭｊ−１＋ｋ　　　　　　　．　　　　　　　　
　（５）ｐ＝　（１／３、１／２、１、２、３）．ｋ＝
−５、−４，・・一・、＋５及び３　９　＜　ｎ　＜　
１　２　１換言すれば、ピッチ関連データＭの高調波及
び低調波に関して位置するサンプルについてサンプル指
標ｎをセッ１へすることである。次に、を計算し、Ｍ２
及びｂ２の導出のためＲ（ｎ）の分類（選出）にもどる
。

最終的に、式（５）で始まるプロセスはＭ３、ｂ３及び
Ｍ４、ｂ４を導出するために繰返される。

以上本発明のプロセスを．ＲＰＥ手法を用いることによ
って低レー１〜が達成されるような特定のコーダの例に
関して説明したが、本発明はＭ　Ｐ　Ｅ（ｍｕｌｔｉｐ
ｌｅ　ｅｘｃｉｔａｔｉｏｎ）やＣ　Ｅ　Ｌ　Ｐ　（ｃ
ｏｄｅｅｘｃｉｔｅｄ　ｌｉｎｅａｒ　ｐｒｅｄｉｃｔ
ｉｖｅ　ｃｏｄｉｎｇ）のような他の低レートコード化
手法にも適用できる。また、ｒ　（ｎ）はフルバンド残
差又はベースバンド残差のいずれでもよい。

Ｅ．発明の効果 以上説明したように本発明よれば、ピーク検出に関する
計算負荷を減ずることができる。

【図面の簡単な説明】

第１図及び第２図は本発明を適用することのできる音声
コーダを示す図、第３Ａ図、第３Ｂ図及び第４図は本発
明に基づくプロセスを表わすアルゴリズムを説明する図
である。 出願人　　インターナショナル・ビジネス・マシーンズ
・コーポレーション 代理人　　弁理士　　頓　　宮　　孝 （外１名）

Claims (1)

1. 【特許請求の範囲】 音声から導出されたデジタル信号ｒ（ｎ）を低ビットレ
   ートの信号に変換する長区間予測に基づく音声コーダで
   使用される長区間予測フィルタを調整するためのピッチ
   関連遅延値Ｍを導出する下記の（ａ）ないし（ｈ）のス
   テップを有する方法。 （ａ）上記信号ｒ（ｎ）を各々がＮ個のサンプルから成
   るセグメントに分割するステップ。 （ｂ）各セグメントをｊ個（ｊは所定の整数）のサブセ
   グメントに分割するステップ。 （ｃ）現信号の最初のサブセグメントと復元された信号
   のセグメントと相互相関をとって下記の相互相関関数Ｒ
   （ｎ）を導出するステップ。 Ｒ（ｎ）＝Σ＾ｋ’＾−＾１＿ｉ＿＝＿０ｒ（ｉ）ｙ（
   ｉ−ｎ）ｎ＝ｋ’ないしＮただしｋ’＝Ｎ／ｊ （ｄ）ピーク位置Ｒ（Ｍ１）についてＲ（ｎ）の値を選
   出し、上記フィルタの遅延値をＭ１に設定し、信号をサ
   ブセグメント１個分についてシフトするステップ。 （ｅ）Ｍ１の高調波及び低調波のところに位置する所定
   数のサンプルについてのサンプル指標ｎを計算するステ
   ップ。 （ｆ）ステップ（ｅ）で定義されたｎについて相互相関
   関数の値Ｒ（ｎ）を計算するステップ。 （ｇ）新しい遅延値Ｍ２を導出するためピーク位置につ
   いてＲ（ｎ）の値を選出するステップ。 （ｈ）上記（ｅ）から（ｇ）までのステップをＭ１の代
   わりにＭ２について実行し、以下同様にＭｊになるまで
   （ｅ）から（ｇ）までのステップを繰返すステップ。