JP3233184B2

JP3233184B2 - 音声符号化方法

Info

Publication number: JP3233184B2
Application number: JP05300695A
Authority: JP
Inventors: 一則間野
Original assignee: Nippon Telegraph and Telephone Corp
Current assignee: Nippon Telegraph and Telephone Corp
Priority date: 1995-03-13
Filing date: 1995-03-13
Publication date: 2001-11-26
Anticipated expiration: 2016-11-26
Also published as: JPH08248995A

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、音声信号を４ｋｂｉｔ
／ｓ程度の情報量で符号化する高能率ディジタル音声符
号化、特に、符号励振線形予測符号化の励振信号符号化
等において用いて好適な音声符号化方法に関するもので
ある。

【０００２】

【従来の技術】ディジタル移動通信や、音声蓄積サービ
スでは、電波や記憶媒体の効率的利用を図るために、種
々の高能率音声符号化法が用いられている。その代表的
な基本技術として、符号励振線形予測符号化（Code Exc
ited Linear Prediction，CELP）がある。この技術の詳
細は、M.R.Schroeder and B.S.Atal：“Code-Excited L
inear Prediction (CELP)：High-quality Speech at Ve
ry Low Bit Rates”，Proc. IEEE ICASSP-85，25.1.1，
pp.937-940，（1985年）に述べられている。

【０００３】ＣＥＬＰ符号化おいては、５ｍｓから４０
ｍｓ程度を１フレームとし、あるいは、フレームをさら
に分割してサブフレームとし、フレーム、あるいはサブ
フレームの信号を１ベクトルとして取り扱う。そして、
過去の励振信号からなるピッチ適応符号帳中の一つの適
応符号ベクトルと、あらかじめ蓄積しておいた固定的な
雑音又はパルス列からなる雑音符号帳の雑音符号ベクト
ルとの重み付き和を励振信号とする。この励振信号を線
形予測合成フィルタに通した合成波形と入力音声との聴
覚重みつき波形歪みを最小とするように、適応符号、雑
音符号、利得符号を決定する。

【０００４】ＣＥＬＰ符号化を基本として、さらに音声
のピッチ周期性を雑音符号帳に反映させた方法としてピ
ッチ同期雑音励振源符号励振線形予測符号化（Pitch Sy
nchronous Innovation CELP, PSI-CELP）がある。この
技術の詳細は、三樹、守谷、間野、大室：“ピッチ同期
雑音励振源をもつＣＥＬＰ符号化(PSI-CELP)”、電子情
報通信学会論文誌Ａ，Vol.J77-A, No.3, pp.285-292（1
994年３月）に記述されている。前述したＣＥＬＰの場
合には、雑音符号帳から得られる雑音符号ベクトルが常
に固定されたものであったが、ＰＳＩ−ＣＥＬＰでは、
適応符号帳から得られるピッチ周期を用い、適応符号帳
と同様にして、雑音符号帳に格納されたベクトルを先頭
から周期化することによって、雑音符号ベクトルを得て
おり、この点が両者で異なっている。

【０００５】図３に従来のＰＳＩ−ＣＥＬＰ符号化法の
基本的なブロック図を示す。まず、入力端子１から音声
を入力する。線形予測分析部２において音声の線形予測
分析を行い、量子化された予測係数Ａが合成フィルタ３
の係数となる。４は第１の符号帳であり、直前の過去の
合成フィルタ３への入力として使用された励振音源（Ａ
ＣＢ）が蓄えられたバッファである適応符号帳４ａと、
固定の雑音符号帳４ｂ（以下、固定符号帳と称する。）
からなる。第１の符号帳４へ入力される符号Ｌは、適応
符号帳４ａ内の適応符号ベクトルが選択された場合に
は、音声波形のピッチ周期に対応する。一方、固定符号
帳４ｂ内の固定符号ベクトルが選択された場合、符号Ｌ
は固定符号帳４ｂのベクトルの番号を示す。ただし、こ
れらの第１の符号帳４から出力される適応符号ベクトル
と固定符号ベクトルは、切り換え部１３で排他的に選択
された後、乗算部７へ供給される。

【０００６】なお、第１の符号帳４からは、合成波形歪
みが最小となように各符号ベクトルが選択されるが、適
応符号帳４ａから選択された符号ベクトルは、その後、
図４（１）に示すように符号Ｌに対応するピッチ周期Ｔ
（Ｌ）で過去の励振信号を図３では図示を省略した周期
化部４ａ−１によって繰り返し、周期化することによっ
て上述した適応符号ベクトルとされ、他方、固定符号帳
からは、図４（２）に示すように、格納されているベク
トルそのものが選択される。なお、図４に示すＮはサン
プリング周期を単位とするサブフレームの長さ、Ｆは任
意のベクトルの番号を示す。

【０００７】一方、図３に示す５は、第２の符号帳であ
り、通常、雑音符号帳と呼ばれる。ＰＳＩ−ＣＥＬＰで
は、雑音符号帳５についても上記と同様に合成波形歪み
最小となるベクトルが選択されるが、その際、第１の符
号帳４のベクトルとして適応符号ベクトルが選択された
場合には、図４（３）に示すように、入力符号Ｃに対応
する雑音符号帳５のＣ番目に格納されている雑音符号ベ
クトルが選択され、さらに周期化部６において符号Ｌに
対応するピッチ周期Ｔ（Ｌ）で周期化され、この周期化
されたベクトルが励振信号となる。他方、符号Ｌが固定
符号帳４ｂを指定した場合には、周期化は行われずに、
ベクトルそのものとして選択され、切り換え部１４へ供
給される。

【０００８】第１の符号帳４と第２の符号帳（雑音符号
帳）５から選択された第１、第２符号帳ベクトルには、
それぞれ、乗算部７、８において、波形歪みが最小とな
るように最適化された利得符号Ｇ₀、Ｇ₁が乗算される。
そして、この乗算部７、８の出力は、加算部１１で合成
された後、合成フィルタ部３へ入力される。

【０００９】また、９は、入力波形と合成波形との聴覚
重み付きの歪みを求める歪み計算部であり、１０は、最
適な符号検索を行い、送出符号を決定する符号帳検索部
である。なお、ＰＳＩ−ＣＥＬＰに対して従来のＣＥＬ
Ｐ符号化とは、図３で第１の符号帳４を適応符号帳４ａ
のみとし、第２の符号帳５から選択された符号ベクトル
が周期化されない場合に相当する。

【００１０】次に、ＰＳＩ−ＣＥＬＰの第２の符号帳
（雑音符号帳）５の検索について述べる。入力音声信号
から、減算部１２において、合成フィルタ３の過去のフ
レームからの零入力応答と第１の符号帳４から選択され
た符号ベクトルに最適ゲインＧ₀を乗じて合成フィルタ
３を通した合成波形を引いたものが目標ベクトルとな
る。そして、符号帳検索部１０では、この目標ベクトル
と、第２の符号帳５のベクトルに利得Ｇ₁を乗じて合成
フィルタ３を通した信号との歪が小さくなるように、符
号ベクトルが選択される。なお、一般に、歪の計算で
は、聴感重み付けフィルタによる重みが目標ベクトルと
合成音声の両方に掛けられるが、重み付けを行うか否か
は、以下の説明では特に問題とはならないので、以下、
重み付けについて、省略して説明を行う。

【００１１】ここで、サブフレーム長をＮ、目標ベクト
ルをＸ＝（ｘ₀，ｘ_l，……，ｘ_N-1）^tとし、第２の雑音
符号帳５の符号ベクトルをＥ＝（ｅ₀，ｅ₁，…，
ｅ_N-1）^t、利得をｇ、合成フィルタ３のインパルス応答
｛ｈ₀，ｈ₁，…，ｈ_N-1｝を用いたインパルス応答行列
Ｈを

【００１２】

【数２】

【００１３】とすると、第２の雑音符号Ｃの決定時の歪
尺度Ｄは次式となる。Ｄ＝｜Ｘ−ｇＨＥ｜² ＝｜Ｘ｜²−２ｇ（Ｘ^tＨＥ）十ｇ²｜ＨＥ｜² （２）ただし、利得ｇは第２の雑音符号帳ベクトルＥを定めた
後で決定するので、この段階では任意にとる。すなわ
ち、歪尺度Ｄを最小にする利得ｇは、∂Ｄ／∂ｇ＝０よ
り、ｇ＝（Ｘ^tＨＥ）／｜ＨＥ｜²である。また、目標ベ
クトルＸは各サブフレーム内で共通である。従って、歪
尺度Ｄを最小にする雑音符号ベクトルＥの検索には、Ｄ’＝（Ｘ^tＨＥ）²／（Ｅ^tＨ^tＨＥ）（３）を最大とするベクトルＥを第２の符号帳５から選択すれ
ば良い。

【００１４】従来のＣＥＬＰ符号化において、第２の
（雑音）符号帳から選択される雑音符号ベクトルＥを、
少数のパルスからなるもの、例えば、零要素９０％程
度、非零要素１０％程度のスパース性を有するものとす
ると、（３）式を計算する場合には、雑音符号ベクトル
Ｅの非零の要素に限って、行列とベクトルの演算を行え
ば良いので、符号帳検索時の演算量を非常に小さくする
ことができ、かつ、符号帳を格納するためのメモリ量を
小さくすることができる。

【００１５】例えば、雑音符号ベクトルＥとして非零要
素が２つのパルス性ベクトルＥ＝（０，……，０，
ｅ_i，０，……，０，ｅ_j，０，……，０）^tとする。こ
の場合、ｄ^t＝Ｘ^tＨ＝（ｄ₀，ｄ₁，…，ｄ_N-1）をあら
かじめ計算しておけば、（３）式の分子中のＸ^tＨＥ
は、Ｘ^tＨＥ＝ｄ^tＥ＝ｄ_iｅ_i＋ｄ_jｅ_j （４）によって計算できるので、Ｎ次元のベクトルの内積が２
項の積和で計算できることになる。

【００１６】また、（３）式の分母については、Ｈ^tＨ
＝Ｆ＝｛ｆ_ij｝（この場合ｆ_ij＝ｆ_j _i）をあらかじめ計
算しておくと、Ｅ^tＨ^tＨＥ＝Ｅ^tＦＥ＝ｅ² _iｆ_ii＋２ｅ_iｅ_jｆ_ij＋ｅ² _jｆ_jj （５）によって計算できるので、これによっても大きく演算量
を低減することができる。

【００１７】

【発明が解決しようとする課題】しかし、ＰＳＩ−ＣＥ
ＬＰ符号化に関して、従来のＣＥＬＰ符号化と同様の少
数パルスの格納ベクトルからなる雑音符号ベクトルによ
る符号化方法を実現するには、次のような問題点があ
る。ＰＳＩ−ＣＥＬＰ符号化においては、雑音符号ベク
トルが、適応符号のピッチ周期に応じて周期化される。
このピッチ周期化によって、１サブフレームあたりのパ
ルス数がピッチ周期に依存して増減したり、又は非整数
ピッチ周期での周期化では第１ピッチでのパルスが次の
ピッチ周期で非整数値サンプルの補間によって非パルス
的となったりして、スパース性がなくなるので、インパ
ルス応答行列演算による演算量低減の利点が失われてし
まうという問題があった。

【００１８】この発明はこのような背景の下になされた
ものであり、その目的は、ＰＳＩ−ＣＥＬＰ符号化に関
して、少数パルスからなる雑音ベクトル符号帳構成で、
低演算量での符号帳検索を可能とする低ビットレートの
音声符号化方法を提供することにある。

【００１９】

【課題を解決するための手段】上記の課題を解決するた
めに、本発明は、第２の符号帳の各符号ベクトルが最低
１個の零であるベクトル要素を含んでなり、第１の符号
帳を構成する適応符号帳から第１の符号ベクトルが選択
された場合、ピッチ周期に対応して第２の符号帳の格納
ベクトルをピッチ周期化して第２の符号ベクトルとし、
第２の符号ベクトルの整数ピッチ周期長に相当する部分
ベクトルと、線形予測合成フィルタのインパルス応答に
対してピッチ周期ごとにその成分を加え合わせたピッチ
周期性を持つ応答行列を含む評価値を作り、該評価値が
最大となるように、第２の符号帳から第２の符号ベクト
ルを探し出すことを特徴とする。

【００２０】

【作用】第２の（雑音）符号帳検索での符号帳を構成す
る格納ベクトルとして零要素を多く含むベクトルをＣ＝
（ｃ₀，ｃ₁，…，ｃ_N-1）とする。ここで、Ｃの要素の
大部分は零である。また、第１の符号帳検索でピッチ周
期性のあるベクトルが選択された場合に、従来のＰＳＩ
−ＣＥＬＰでは、格納ベクトルＣをピッチ周期化して、
周期化ベクトルＥを計算する。このとき、ピッチ周期化
処理は、次の行列Ｐ（Ｎ×Ｎ）（周期化行列Ｐ）で表現
できる。

【００２１】

【数３】

【００２２】ここで、ｉｎｔ［Ｔ（Ｌ）］は、ピッチ周
期Ｔ（Ｌ）を越えない最大整数値である。また、ベクト
ル（…，ｆ_-2，ｆ_-1，ｆ₀，ｆ₁，ｆ₂，…）は、２ピッ
チ周期めの補間関数の係数、ベクトル（…，ｆ'_-2，ｆ'
_-1，ｆ'₀，ｆ'₁，ｆ'₂，…）は、３ピッチ周期めの補間
関数の係数である。

【００２３】このときのピッチ周期化された雑音符号ベ
クトルＥは、Ｅ＝ＰＣ＝Ｐ_subＣ_sub （７）で表現される。ここで、Ｐ_subは行列Ｐの左側（Ｎ×ｉ
ｎｔ［Ｔ（Ｌ）］）行列であり、Ｃ_subは、Ｃ_sub＝（ｃ
₀，ｃ₁，…，ｃ_int[T(L)]）であり、ベクトルＣの整数
ピッチ周期次元の部分ベクトルである。ただし、周期性
の無い場合には、Ｃ_sub＝Ｃとする。

【００２４】従来のＰＳＩ−ＣＥＬＰでは、ベクトルＣ
_subとしてスパースであっても（６）式の行列Ｐによ
り、雑音符号ベクトルＥはスパースではない全次元に非
零要素をもつ可能性のあるベクトルとなる。このような
雑音符号ベクトルＥを（３）式に代入して計算を行う方
式では雑音符号帳の格納に用いるメモリ量は節約できる
が、Ｅは必ずしもスパースとは限らないので符号帳検索
の演算量は節約できない。

【００２５】そこで、本発明では、雑音符号ベクトルＥ
の整数ピッチ周期長に相当する部分のベクトルと、合成
フィルタのインパルス応答にピッチ周期ごとにその成分
を加え合わせたピッチ周期性をもつ応答行列との積によ
って得られる音声信号に基づいて、歪最小化となる第２
の符号帳検索を行なうようにする。すなわち、この場
合、第２の雑音符号Ｃを決定する際の歪尺度を評価値Ｄ
₂とすると、（３）式と（７）式より、

【００２６】Ｄ₂＝（Ｘ^tＨＥ）²／（Ｅ^tＨ^tＨＥ）＝（Ｘ^tＨ(ＰＣ_sub)）²／（(ＰＣ_sub)^tＨ^tＨ(ＰＣ_sub)）＝（(Ｘ^tＨＰ)Ｃ_sub）²／（Ｃ_sub ^t((ＨＰ)^t(ＨＰ))Ｃ_sub）（８）

【００２７】となるので、従来のインパルス応答行列Ｈ
に対して、行列Ｐを右から乗じて周期性をもたせたピッ
チ周期性のある応答行列ＨＰを最初に計算し、それを用
いて、スパースベクトルＣ_subとの分子及び分母のベク
トル演算を行うとスパース性を利用できるので演算量を
低減できる。

【００２８】これを言い換えると、ＰＣ_subは、非零成
分の多いベクトルなので、Ｈ・（ＰＣ_sub）の演算には
スパース性を利用できないが、符号帳の各パターンであ
るＣ_s _ubをスパースにしておいて、（ＨＰ）を予め計算
してから、（ＨＰ）Ｃ_subの計算をするようにすれば、
（４）、（５）式の効果が得られるということである。

【００２９】一方、第１の符号帳検索でピッチ周期性の
ないベクトルが選択された場合には、周期化を行わない
ので、スパース化が保たれた雑音符号ベクトルと、合成
フィルタのインパルス応答系列からなる行列との積によ
って得られる音声信号に基づいて、歪最小化となる演算
量を低減した第２の符号帳検索を行う。こうして、本発
明の目的であるスパース符号帳による符号帳検索の演算
量と雑音符号帳格納のメモリ量を低減したＰＳＩ−ＣＥ
ＬＰ符号化を実現することが可能となる。

【００３０】

【実施例】図１に、本発明の実施例としてスパース雑音
符号帳をもつ音声符号化のブロック図を示す。図３と同
じ番号のモジュールは、図３と同じものである。まず、
入力端子１から音声ベクトルＸ_inputを入力する。線形
予測分析部２において音声の線形予測分析を行い、量子
化された予測係数Ａあるいは、予測係数Ａから求めたイ
ンパルス応答行列Ｈが合成フィルタ３ａの係数となる。
４は第１の符号帳であり、適応符号帳４ａと固定雑音符
号帳４ｂからなる。符号Ｌによって、第１の符号帳４か
ら適応符号ベクトルが選択される場合には、符号Ｌが音
声波形のピッチ周期に対応する。一方、固定符号ベクト
ルが選択された場合には、符号Ｌは固定符号帳のベクト
ルの番号を示す。この第１の符号帳選択は図３に示す従
来のＣＥＬＰで行われるものと同じである。

【００３１】入力音声ベクトルＸ_inputから、過去のフ
レームからの零入力応答ベクトルＸ_z _irを引き、さら
に、第１の符号帳４から選択されたベクトルに最適利得
Ｇ₀を乗じ、合成フィルタ３ａを通した合成波形Ｘ_adpを
引いたものが目標ベクトルＸである。

【００３２】他方、この図に示す５０は、第２の符号帳
であり、この場合、これをスパース雑音符号帳とする。
ただし、本実施例では、第２の符号帳５０をスパースな
雑音符号帳として構成するときに、各ベクトルを、非零
成分の位置情報と振幅情報とによって表現し、バッファ
メモリに格納する。例えば、ベクトルの次元を８０とす
ると、仮に、全次元とも非零成分をとりうるとした場合
に必要な１ベクトル当たりのメモリ量は、１成分当たり
１ワードとすると８０ワードである。これに対して、本
実施例では、第２の符号帳５０に必要なメモリ量は、例
えば非零成分の個数が８、その位置情報と振幅情報を表
現するのにそれぞれ１ワードとすれば、１６ワードのメ
モリ量で済む。

【００３３】また、スパース雑音符号帳５０の検索時に
は、スイッチ１００を接点１０２の方に接続する。第１
の符号帳４の符号Ｌによって周期化行列部２００で周期
化行列Ｐを構成し、線形予測分析部２から供給されるイ
ンパルス応答行列Ｈとともに、ピッチ周期化合成フィル
タ計算部２０１において行列ＨＰを生成する。行列ＨＰ
の計算は、実際には、（Ｎ×Ｎ）行列どうしの全行ベク
トルと全列ベクトルの内積のすべてを計算する必要はな
く、行列Ｐの各要素と行列Ｈのインパルス応答系列ｈ_n
を用いてｈｐ_nを計算し、それを行列Ｈのｈ_nの代わりに
配置することによって行われる。例えば、（１）、
（６）式より２ピッチ周期めのピッチ周期成分を含む応
答行列ｈｐ_nは、次のようにして得られる。

【００３４】

【数４】

【００３５】この行列ＨＰと、目標ベクトルＸと、各ス
パース雑音ベクトルＣの候補に関して、符号帳検索部２
０２で歪最小となる雑音符号Ｃを求める。

【００３６】最終的な音声出力を求める場合には、第２
の符号帳５０の出力をスイッチ１００を接点１０１の方
に接続する。図３と同様に、符号Ｌの値に応じて適応符
号帳４ａが選択された場合には、周期化部６で周期化
し、固定符号帳４ｂが選択された場合には、周期化せず
に、乗算部８で利得Ｇ₁を乗じて第２の符号ベクトルを
生成する。そして、第１の符号帳４からの符号ベクトル
との和によって励振信号ベクトルを生成し、合成フィル
タ３ｂを通して音声出力を得る。

【００３７】ここで、第２の符号帳５０の雑音符号Ｃを
検索する場合について、そのフローを図２を参照して説
明する。なお、ここで、検索とは、（８）式の評価値Ｄ
₂の計算を各雑音符号Ｃについて実行していくことであ
る。まず、図１に示す周期化行列部２００、ピッチ周期
化合成フィルタ計算部２０１によって、ピッチ周期化用
の合成フィルタ行列ＨＰを算出する（ステップ４０
０）。そして、符号帳検索部２０２において、ステップ
４０１からの以下の処理を行っていく。ステップ４０１
では、Ｙ^t＝Ｘ^tＨＰによって（８）式に示す評価値Ｄ₂
の分子の共通項を求める。次に、ステップ４０２で、雑
音符号ベクトルＣの全候補（候補数Ｋ）、ベクトルＣ
_sub ⁽ⁱ⁾，（ｉ＝０，…，Ｋ−１）に対して、評価値Ｄ₂
の分子Ｊ_N ⁽ⁱ⁾＝（Ｙ^tＣ_sub ⁽ⁱ⁾）²を計算する。この場
合、ベクトルＣ_s _ub ⁽ⁱ⁾は、スパースなため、非零成分に
ついてのみ計算を行うようにすれば、この計算はフルパ
ルスの場合に比較して高速に計算できることができる。

【００３８】次に、ステップ４０３で、ステップ４０２
で求めたＪ_N ⁽ⁱ⁾について、大きい順にＭ個の候補｛Ｊ_N
⁽ⁱ⁾；ｉ＝ｉ₀，…，ｉ_M-1｝を予備選択する。ステップ
４０４では、予備選択されたＭ個の候補について評価値
Ｄ₂の分母Ｊ_D ⁽ⁱ⁾＝（Ｃ_sub ⁽ⁱ⁾）^t（ＨＰ）^t（ＨＰ）Ｃ
_sub ⁽ⁱ⁾の計算を行う。そして、ステップ４０５におい
て、Ｍ個の評価値Ｄ₂ ⁽ⁱ⁾＝Ｊ_N ⁽ⁱ⁾／Ｊ_D ⁽ⁱ⁾を最大とする
ベクトルＣ_sub ⁽ⁱ⁾を決定する。

【００３９】

【発明の効果】以上のように、本発明では、ＰＳＩ−Ｃ
ＥＬＰ音声符号化において、第２の符号帳の各符号ベク
トルが最低１個の零であるベクトル要素を含んでなり、
第１の符号帳を構成する適応符号帳から第１の符号ベク
トルが選択された場合、第２の符号ベクトルの整数ピッ
チ周期長に相当する部分ベクトル（Ｃ_sub）と、線形予
測合成フィルタのインパルス応答に対してピッチ周期ご
とにその成分を加え合わせたピッチ周期性を持つ応答行
列（Ｈ・Ｐ）を含む評価値（Ｄ₂）を作り、該評価値が
最大となるように、第２の符号帳から第２の符号ベクト
ルを探し出すので、メモリ量と演算量を低減した低ビッ
トレート音声符号化を実現することができる。

【図面の簡単な説明】

【図１】本発明の実施例によるスパース雑音符号帳をも
つ音声符号化のブロック図である。

【図２】本発明の実施例によるスパース雑音符号帳検索
のフロー図である。

【図３】従来のＰＳＩ−ＣＥＬＰ音声符号化のブロック
図である。

【図４】ＰＳＩ−ＣＥＬＰ音声符号化の励振信号生成の
説明図である。

【符号の説明】

２線形予測分析部３ａ、３ｂ合成フィルタ部４第１の符号帳４ａ適応符号帳４ｂ固定雑音符号帳５０第２の符号帳２００周期化行列２０１ピッチ周期化合成フィルタ計算部２０２符号帳検索部

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.⁷，ＤＢ名) G10L 19/00 - 19/14 H04B 14/04 H03M 7/30

Claims

(57)【特許請求の範囲】

【請求項１】過去の駆動音源ベクトルをピッチ周期で
繰り返す適応符号帳と予め蓄積してある固定の雑音符号
帳により構成され、前記適応符号帳と前記雑音符号帳の
何れかからひとつ第１の符号ベクトルが選択される第１
の符号帳と、該第１の符号帳における前記固定の雑音符
号帳とは異なる雑音符号帳である第２の符号帳とから音
源ベクトル符号帳が成り立ち、前記第１の符号帳において前記適応符号帳から符号ベク
トルが選択された場合、前記ピッチ周期に対応して前記
第２の符号帳の格納ベクトルをピッチ周期化して第２の
符号ベクトルとし、又は前記固定の雑音符号帳から符号
ベクトルが選択された場合、前記第２の符号帳の格納ベ
クトルをピッチ周期化せずに第２の符号ベクトルとし、前記第１の符号ベクトルと前記第２の符号ベクトルの重
み付き和によって、線形予測合成フィルタを駆動して、
フレーム単位に音声信号を合成し、入力信号と前記合成された音声信号との歪みを最小とす
るように、前記第１の符号帳と前記第２の符号帳からそ
れぞれ符号ベクトルを探し出す音声符号化方法におい
て、前記第２の符号帳の各符号ベクトルが最低１個の零であ
るベクトル要素を含んでなり、前記第１の符号帳を構成する前記適応符号帳から前記第
１の符号ベクトルが選択された場合、前記第２の符号ベ
クトルの整数ピッチ周期長に相当する部分ベクトル（Ｃ
_sub）と、前記線形予測合成フィルタのインパルス応答
に対して前記ピッチ周期ごとにその成分を加え合わせた
ピッチ周期性を持つ応答行列（Ｈ・Ｐ）を含む評価値
（Ｄ₂）を作り、該評価値が最大となるように、前記第
２の符号帳から前記第２の符号ベクトルを探し出すこと
を特徴とする音声符号化方法。
【請求項２】前記応答行列（Ｈ・Ｐ）は、前記線形予
測合成フィルタのインパルス応答行列Ｈに、下式で示す
周期化行列Ｐを右から乗じて求めた行列であることを特
徴とする請求項１記載の音声符号化方法。【数１】ただし、Ｎは前記フレーム単位の長さ、ｉｎｔ［Ｔ
（Ｌ）］は前記ピッチ周期長Ｔ（Ｌ）を越えない最大整
数値、ベクトル（…，ｆ_-2，ｆ_-1，ｆ₀，ｆ₁，ｆ₂，
…）は２ピッチ周期めの補間関数の係数、ベクトル
（…，ｆ'_-2，ｆ'_-1，ｆ'₀，ｆ'₁，ｆ'₂，…）は３ピッ
チ周期めの補間関数の係数である。