(57)【特許請求の範囲】
【請求項１】
時間領域の入力信号を周波数領域の信号に変換する時間−周波数領域変換手段と、
前記周波数領域の信号を複数のサブベクトルに分割し、分割された各サブベクトルのスペクトル係数を量子化するベクトル量子化手段と、
前記量子化によって得られ、前記複数のサブベクトルのエネルギーの大きさに応じて大きな番号が付与されるコードブック番号を表し、前記コードブック番号が大きいほど使用するビット数が多くなるコードブック指示値について、全ての前記サブベクトルのコードブック指示値を符号化する符号化手段と、
全ての前記コードブック指示値のうち、最多使用ビット数のコードブック指示値となるサブベクトルの位置を特定する特定手段と、
前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定手段と、
前記最多使用ビット数のコードブック指示値が前記符号化手段によって符号化されたビット数を実際値とし、推定された前記最多使用ビット数のコードブック指示値のビット数を推定値として、前記実際値と前記推定値との差分を算出する差分算出手段と、
特定された前記サブベクトルの位置情報と算出された前記差分情報を符号化するパラメータ符号化手段と、
を具備する音声符号化装置。
【請求項２】
前記最多使用ビット数のコードブック指示値の使用ビット数と、所定の閾値との比較結果に応じて、特定された前記サブベクトルの位置情報を符号化するか否かを判定する判定手段を具備する請求項１に記載の音声符号化装置。
【請求項３】
前記差分算出手段は、前記最多使用ビット数のコードブック指示値の使用ビット数が前記閾値よりも大きい場合、前記最多使用ビット数のコードブック指示値が前記符号化手段によって符号化されたビット数を実際値とし、前記実際値と前記推定値との差分を算出する請求項２に記載の音声符号化装置。
【請求項４】
前記推定手段は、前記最多使用ビット数のコードブック指示値の使用ビット数が前記閾値以下の場合、前記入力信号の送信単位において利用可能な総ビット数と、予め定められたサブベクトルのコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記予め定められたサブベクトルのコードブック指示値が使用するビット数を推定し、
前記差分算出手段は、前記予め定められたサブベクトルのコードブック指示値が前記符号化手段によって符号化されたビット数を実際値とし、推定された前記予め定められたサブベクトルのコードブック指示値のビット数を推定値として、前記実際値と前記推定値との差分を算出する請求項２に記載の音声符号化装置。
【請求項５】
ステレオまたはマルチチャネル信号の一つのチャネルまたは複数のチャネルの信号を前記入力信号とする請求項１に記載の音声符号化装置。
【請求項６】
複数フレーム単位または複数サブフレーム単位のスペクトル係数列を入力信号とする請求項１に記載の音声符号化装置。
【請求項７】
請求項１に記載の音声符号化装置より符号化された前記位置情報及び前記差分情報を取得して、復号するパラメータ復号手段と、
前記音声符号化装置より、前記最多使用ビット数のコードブック指示値以外の符号化されたコードブック指示値を取得して、復号するコードブック指示値復号手段と、
前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定手段と、
推定された前記最多使用ビット数のコードブック指示値が使用するビット数と復号された前記差分情報とを加算して前記最多使用ビット数のコードブック指示値を算出する加算手段と、
前記復号された位置情報、前記最多使用ビット数のコードブック指示値以外の復号されたコードブック指示値、及び、算出された前記最多使用ビット数のコードブック指示値を用いて、全てのコードブック指示値を生成するコードブック指示値生成手段と、
生成された前記全てのコードブック指示値を用いて、各サブベクトルのスペクトル係数を逆量子化するベクトル逆量子化手段と、
逆量子化された前記スペクトル係数を時間領域に変換する周波数−時間変換手段と、
を具備する音声復号装置。
【請求項８】
算出された前記最多使用ビット数のコードブック指示値または予め固定された位置のサブベクトルのコードブック指示値の使用ビット数と、所定の閾値との比較結果に応じて、算出された前記最多使用ビット数のコードブック指示値の前記サブベクトルの位置情報を用いて、前記全てのコードブック指示値を生成するか否かを判定する判定手段を具備する請求項７に記載の音声復号装置。
【請求項９】
前記コードブック指示値生成手段は、算出された前記最多使用ビット数のコードブック指示値または前記予め固定された位置のサブベクトルのコードブック指示値の使用ビット数が前記閾値より大きい場合、算出された前記最多使用ビット数のコードブック指示値の前記サブベクトルの位置情報を用いて、前記全てのコードブック指示値を生成する請求項８に記載の音声復号装置。
【請求項１０】
前記コードブック指示値生成手段は、算出された前記最多使用ビット数のコードブック指示値または前記予め固定された位置のサブベクトルのコードブック指示値の使用ビット数が前記閾値以下の場合、前記予め固定された位置のサブベクトルの位置情報を用いて、前記全てのコードブック指示値を生成するか否かを判定する判定手段を具備する請求項８に記載の音声復号装置。
【請求項１１】
復号されたスペクトルを所定数のサブバンドに分割し、分割されたサブバンドをゲイン補正係数によってスケーリングする請求項８に記載の音声復号装置。
【請求項１２】
時間領域の入力信号を周波数領域の信号に変換する時間−周波数領域変換工程と、
前記周波数領域の信号を複数のサブベクトルに分割し、分割された各サブベクトルのスペクトル係数を量子化するベクトル量子化工程と、
前記量子化によって得られ、前記複数のサブベクトルのエネルギーの大きさに応じて大きな番号が付与されるコードブック番号を表し、前記コードブック番号が大きいほど使用するビット数が多くなるコードブック指示値について、全ての前記サブベクトルのコードブック指示値を符号化する符号化工程と、
全ての前記コードブック指示値のうち、最多使用ビット数のコードブック指示値となるサブベクトルの位置を特定する特定工程と、
前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定工程と、
前記最多使用ビット数のコードブック指示値が前記符号化工程において符号化されたビット数を実際値とし、推定された前記最多使用ビット数のコードブック指示値のビット数を推定値として、前記実際値と前記推定値との差分を算出する差分算出工程と、
特定された前記サブベクトルの位置情報と算出された前記差分情報を符号化するパラメータ符号化工程と、
を具備する音声符号化方法。
【請求項１３】
請求項１２に記載の音声符号化方法により符号化された前記位置情報及び前記差分情報を復号するパラメータ復号工程と、
前記音声符号化方法により符号化された、前記最多使用ビット数のコードブック指示値以外のコードブック指示値を復号するコードブック指示値復号工程と、
前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定工程と、
推定された前記最多使用ビット数のコードブック指示値が使用するビット数と復号された前記差分情報とを加算して前記最多使用ビット数のコードブック指示値を算出する加算工程と、
前記復号された位置情報、前記最多使用ビット数のコードブック指示値以外の復号されたコードブック指示値、及び、算出された前記最多使用ビット数のコードブック指示値を用いて、全てのコードブック指示値を生成するコードブック指示値生成工程と、
生成された前記全てのコードブック指示値を用いて、各サブベクトルのスペクトル係数を逆量子化するベクトル逆量子化工程と、
逆量子化された前記スペクトル係数を時間領域に変換する周波数−時間変換工程と、
を具備する音声復号方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ベクトル量子化を利用した音声符号化装置、音声復号装置、音声符号化方法及び音声復号方法に関する。
【背景技術】
【０００２】
オーディオ及び音声の符号化においては、変換符号化と線形予測符号化という二つの主要な符号化方式がある。以下、変換符号化と線形予測符号化について説明する。【０００３】
（１）変換符号化について
変換符号化は、離散フーリエ変換（ＤＦＴ:Discrete Fourier Transform）または修正離散コサイン変換（ＭＤＣＴ:Modified Discrete Cosine Transform）を使用するなどして、時間領域の信号からスペクトル領域の信号に変換し、個々のスペクトル係数を量子化し、符号化する。【０００４】
量子化または符号化の処理では、個々のスペクトル係数の知覚的重要度を決定するために、通常、心理音響モデルが適用され、個々のスペクトル係数は、それらの知覚的重要度に応じて量子化または符号化される。変換符号化は、音楽または一般のオーディオ信号に対して有効である。変換コーデックとしては、例えば、ＭＰ３（MPEG Audio Layer 3）、ＡＡＣ（Advanced Audio Coding）（非特許文献１参照）及びＤｏｌｂｙ ＡＣ３（Audio Code number 3）がある。【０００５】
変換コーデックの簡略な構成を図１に示す。図１（ａ）に示した符号化器１０では、時間−周波数変換部１１において、ＤＦＴまたはＭＤＣＴなどの時間−周波数変換方式によって、時間領域の信号Ｓ（ｎ）が周波数領域の信号Ｓ（ｆ）に変換され、周波数領域の信号Ｓ（ｆ）が心理音響モデル分析部１２及び量子化部１３に出力される。【０００６】
周波数領域の信号Ｓ（ｆ）は、心理音響モデル分析部１２において、心理音響モデル分析が行われ、マスキング曲線が得られる。【０００７】
また、周波数領域の信号Ｓ（ｆ）は、量子化部１３において、量子化ノイズが不可聴となるように、マスキング曲線に従って量子化が行われる。【０００８】
量子化された個々の量子化パラメータは、多重化部１４において多重化され、ビットストリーム情報として復号器側へ送信される。【０００９】
図１（ｂ）に示した復号器２０では、符号化器側から送信された全てのビットストリーム情報が、多重分離部２１において、多重分離される。多重分離された量子化パラメータは、逆量子化部２２において、逆量子化され、周波数領域の信号Ｓ〜（ｆ）に復号される。なお、チルダ（波線符号）「〜」は、図面上では、記号「Ｓ」の上部に付されているが、本明細書では、表記の制約上、記号「Ｓ」の右肩に付している。ここで、チルダは、復号の結果得られた信号であることを表している。【００１０】
復号された周波数領域の信号Ｓ〜（ｆ）は、周波数−時間変換部２３において、逆離散フーリエ変換（ＩＤＦＴ:Inverse Discrete Fourier Transform）または逆修正離散コサイン変換（ＩＭＤＣＴ:Inverse Modified Discrete Cosine Transform）などの周波数−時間変換方式によって、時間領域の信号Ｓ〜（ｎ）に変換される。【００１１】
（２）線形予測符号化について
線形予測符号化は、時間領域における音声信号の予測可能な性質（つまり、時間領域において音声信号が予測可能であるという性質）を利用し、入力された音声信号に対して線形予測を適用することによって、残差信号（または励振信号）を得る。特に、有声範囲の音声信号に対して、この線形予測モデルは、音声を非常に効率的に表現することができる。線形予測の後、残差信号が、主に、ＴＣＸ（Transform Coded eXcitation）またはＣＥＬＰ（Code Excited Linear Prediction）という二つの異なる方式のいずれかによって符号化される。以下、ＴＣＸとＣＥＬＰについて説明する。【００１２】
（２−１）ＴＣＸについて
ＴＣＸ（非特許文献２参照）では、残差信号は、周波数領域において効率的に符号化される。ＴＣＸコーデックとしては、例えば、３ＧＰＰ ＡＭＲ‐ＷＢ＋（Extended Adaptive Multi-Rate Wideband）やＭＰＥＧ ＵＳＡＣ（Unified Speech and Audio Coding）がある。【００１３】
ＴＣＸコーデックの簡略な構成を図２に示す。図２（ａ）に示した符号化器３０では、ＬＰＣ分析部３１において、時間領域における信号の予測可能な性質を利用するために、入力信号Ｓ（ｎ）に対してＬＰＣ分析が行われる。【００１４】
ＬＰＣ分析から生じた個々のＬＰＣパラメータは、量子化部３２において量子化され、量子化インデックスが逆量子化部３３及び多重化部３７に出力される。【００１５】
量子化インデックスは、逆量子化部３３において逆量子化され、ＬＰＣパラメータが復元される。【００１６】
また、入力信号Ｓ（ｎ）は、ＬＰＣ逆フィルタ部３４において、復元されたＬＰＣパラメータを用いたＬＰＣ逆フィルタリングがかけられることによって、時間領域の残差信号Ｓｒ（ｎ）が得られる。【００１７】
時間領域の残差信号Ｓｒ（ｎ）は、時間−周波数変換部３５において、ＤＦＴまたはＭＤＣＴなどの時間−周波数変換方式によって、周波数領域の残差信号Ｓｒ（ｆ）に変換される。【００１８】
周波数領域の残差信号Ｓｒ（ｆ）は、量子化部３６において量子化され、個々の量子化パラメータが多重化部３７に出力される。【００１９】
量子化部３２から出力された量子化インデックスと、量子化部３６から出力された個々の量子化パラメータは、多重化部３７において、多重化され、ビットストリーム情報として復号器側へ送信される。【００２０】
図２（ｂ）に示した復号器４０では、符号化器側から送信された全てのビットストリーム情報が、多重分離部４１において、量子化インデックスと量子化パラメータとに多重分離される。多重分離された量子化インデックスは逆量子化部４４に、多重分離された量子化パラメータは逆量子化部４２に出力される。【００２１】
多重分離された量子化パラメータは、逆量子化部４２において逆量子化され、周波数領域の残差信号Ｓ〜ｒ（ｆ）に復号され、復号された周波数領域の残差信号Ｓ〜ｒ（ｆ）は、周波数−時間変換部４３において、ＩＤＦＴまたはＩＭＤＣＴなどの周波数−時間変換方式によって、時間領域の残差信号Ｓ〜ｒ（ｎ）に変換される。【００２２】
また、多重分離された量子化インデックスは、逆量子化部４４において逆量子化され、ＬＰＣパラメータが得られる。【００２３】
時間領域の残差信号Ｓ〜ｒ（ｎ）は、ＬＰＣ合成フィルタ部４５において、ＬＰＣパラメータを用いた処理が行われ、時間領域の信号Ｓ〜（ｎ）が得られる。【００２４】
（２−２）ＣＥＬＰについて
ＣＥＬＰでは、残差信号は、所定のコードブックを使用して量子化される。また、音声品質をさらに向上させるために、元の信号とＬＰＣ合成後の信号との差分信号を周波数領域に変換して、符号化することがよく行われる。ＣＥＬＰコーデックとしては、例えば、ＩＴＵ‐Ｔ Ｇ．７２９．１（非特許文献３参照）やＩＴＵ‐Ｔ Ｇ．７１８（非特許文献４参照）がある。【００２５】
ＣＥＬＰと変換符号化の階層符号化（またはエンベデッド符号化）の簡略な構成を図３に示す。図３（ａ）に示した符号化器５０では、ＣＥＬＰ符号化部５１において、時間領域における信号の予測可能な性質を利用するために、入力信号Ｓ（ｎ）に対してＣＥＬＰ符号化が行われ、ＣＥＬＰパラメータがＣＥＬＰローカル復号部５２及び多重化部５５に出力される。【００２６】
ＣＥＬＰパラメータは、ＣＥＬＰローカル復号部５２において復号され、合成信号Ｓｓｙｎ（ｎ）が得られる。入力信号Ｓ（ｎ）から合成信号Ｓｓｙｎ（ｎ）を減算することにより、予測誤差信号Ｓｅ（ｎ）が得られる。【００２７】
時間領域の予測誤差信号Ｓｅ（ｎ）は、時間−周波数変換部５３において、ＤＦＴまたはＭＤＣＴなどの時間−周波数変換方式によって、周波数領域の予測誤差信号Ｓｅ（ｆ）に変換される。【００２８】
周波数領域の予測誤差信号Ｓｅ（ｆ）は、量子化部５４において量子化され、個々の量子化パラメータが多重化部５５に出力される。【００２９】
ＣＥＬＰ符号化部５１から出力されたＣＥＬＰパラメータと、量子化部５４から出力された個々の量子化パラメータは、多重化部５５において多重化され、ビットストリーム情報として復号器側へ送信される。【００３０】
図３（ｂ）に示した復号器６０では、符号化器側から送信された全てのビットストリーム情報が、多重分離部６１において、ＣＥＬＰパラメータと個々の量子化パラメータとに多重分離される。多重分離されたＣＥＬＰパラメータはＣＥＬＰ復号部６４に、多重分離された量子化パラメータは逆量子化部６２に出力される。【００３１】
多重分離された量子化パラメータは、逆量子化部６２において、逆量子化され、周波数領域の予測誤差信号Ｓ〜ｅ（ｆ）に復号され、復号された周波数領域の予測誤差信号Ｓ〜ｅ（ｆ）は、周波数−時間変換部６３において、ＩＤＦＴまたはＩＭＤＣＴなどの周波数−時間変換方式によって、時間領域の予測誤差信号Ｓ〜ｅ（ｎ）に変換される。【００３２】
また、多重分離されたＣＥＬＰパラメータは、ＣＥＬＰ復号部６４において復号され、合成信号Ｓｓｙｎ（ｎ）が得られる。【００３３】
予測誤差信号Ｓ〜ｅ（ｎ）と合成信号Ｓｓｙｎ（ｎ）とが加算されることにより、時間領域の信号Ｓ〜（ｎ）が得られる。【００３４】
（３）スプリットマルチレート格子ベクトル量子化について
このような、変換符号化及び線形予測符号化における符号化は、通常、何らかの量子化法を利用している。この量子化法の一つに、スプリットマルチレート格子ベクトル量子化（以下、適宜、「スプリットマルチレート格子ＶＱ」という）（または代数的ベクトル量子化）という方法がある（非特許文献５参照）。【００３５】
ＡＭＲ‐ＷＢ＋（非特許文献６参照）では、ＴＣＸ領分におけるＬＰＣの残差を量子化するために、スプリットマルチレート格子ＶＱが使用される。新たに標準化された音声コーデックであるＩＴＵ‐Ｔ Ｇ．７１８においても、ＭＤＣＴ領分におけるＬＰＣの残差を第３の残差符号化層として量子化するために、スプリットマルチレート格子ＶＱが使用される。【００３６】
スプリットマルチレート格子ＶＱは、格子量子化器に基づいたベクトル量子化法である。具体的に、ＡＭＲ‐ＷＢ＋で使用されるスプリットマルチレート格子ＶＱの場合には、ＲＥ８格子と呼ばれるＧｏｓｓｅｔ格子のサブセットにより構成されるベクトルコードブックを使用して、スペクトルが８個のスペクトル係数のブロックを単位として量子化される（非特許文献５参照）。【００３７】
任意の格子の全ての点は、その格子のいわゆる２乗生成マトリクスＧから、ｃ＝ｓ・Ｇ（ｓは個々の整数値を含む線ベクトル、ｃは生成される格子点である）として生成可能である。【００３８】
あるレート（比率）でのベクトルコードブックを作成するためには、定められた半径のある範囲（８次元）内の格子点のみが採取される。よって、マルチレートコードブックは、それぞれ異なる半径の範囲内の格子点の各サブセットを採取することによって作成される。【００３９】
ＴＣＸコーデックにおいてスプリットマルチレート格子ＶＱを利用した簡略な構成を図４に示す。図４（ａ）に示した符号化器７０では、時間領域における信号の予測可能な性質を利用するために、ＬＰＣ分析部７１において、入力信号Ｓ（ｎ）に対してＬＰＣ分析が行われる。【００４０】
ＬＰＣ分析から生じた個々のＬＰＣパラメータは、量子化部７２において量子化され、量子化インデックスが逆量子化部７３及び多重化部７７に出力される。【００４１】
量子化インデックスは、逆量子化部７３において逆量子化され、ＬＰＣパラメータが復元される。【００４２】
また、入力信号Ｓ（ｎ）は、ＬＰＣ逆フィルタ部７４において、復元されたＬＰＣパラメータを用いたＬＰＣ逆フィルタリングがかけられることによって、残差信号Ｓｒ（ｎ）が得られる。【００４３】
時間領域の残差信号Ｓｒ（ｎ）は、時間−周波数変換部７５において、ＤＦＴまたはＭＤＣＴなどの時間−周波数変換方式によって、周波数領域の残差信号Ｓｒ（ｆ）に変換される。【００４４】
周波数領域の残差信号Ｓｒ（ｆ）は、スプリットマルチレート格子ＶＱ部７６において、スプリットマルチレート格子ＶＱが適用され、個々の量子化パラメータが多重化部７７に出力される。【００４５】
量子化部７２から出力された量子化インデックスと、スプリットマルチレート格子ＶＱ部７６から出力された個々の量子化パラメータは、多重化部７７において多重化され、ビットストリーム情報として復号器側へ送信される。【００４６】
図４（ｂ）に示した復号器８０では、符号化器側から送信された全てのビットストリーム情報が、多重分離部８１において、量子化インデックスと量子化パラメータとに多重分離される。【００４７】
多重分離された量子化パラメータは、スプリットマルチレート格子逆ＶＱ部８２において、スプリットマルチレート格子逆ＶＱが適用され、周波数領域の残差信号Ｓ〜ｒ（ｆ）に復号され、復号された周波数領域の残差信号Ｓ〜ｒ（ｆ）は、周波数−時間変換部８３において、ＩＤＦＴまたはＩＭＤＣＴなどの周波数−時間変換方式によって、時間領域の残差信号Ｓ〜ｒ（ｎ）に変換される。【００４８】
また、多重分離された量子化インデックスは、逆量子化部８４において逆量子化され、ＬＰＣパラメータが得られる。【００４９】
時間領域の残差信号Ｓ〜ｒ（ｎ）に対して、ＬＰＣ合成フィルタ部８５において、ＬＰＣパラメータを用いた処理が行われ、時間領域の信号Ｓ〜（ｎ）が得られる。【００５０】
図５は、スプリットマルチレート格子ＶＱの処理を示すブロック図である。図５において、入力スペクトルＳ（ｆ）は、ブロック分割部９１において、ある数の８次元のブロック（またはベクトル）に分割され、分割された８次元のブロックがマルチレート格子ＶＱ部９２に出力される。【００５１】
分割された８次元の各ブロックは、マルチレート格子ＶＱ部９２において、マルチレート格子ＶＱによって量子化される。この量子化では、まず、スペクトル全体の使用可能なビット数とエネルギーレベルにより、グローバルゲインが算出される。次に、ブロックごとに、元のスペクトルとグローバルゲインとの間の比率を求め、その比率はそれぞれ異なるコードブックによって量子化される。【００５２】
これにより得られる、スプリットマルチレート格子ＶＱの個々の量子化パラメータは、グローバルゲインの量子化インデックス、各ブロックについてのコードブック指示値及び各ブロックについてのコードベクトルインデックスである。【００５３】
図６は、ＡＭＲ‐ＷＢ＋で採用されたスプリットマルチレート格子ＶＱのコードブックリストの概要を示す（非特許文献６参照）。図６では、コードブックＱ０、Ｑ２、Ｑ３またはＱ４が、基本コードブックである。ある格子点がこれらの基本コードブックに含まれていない場合には、基本コードブックのＱ３またはＱ４部分のみを使用して、Ｖｏｒｏｎｏｉ拡張（非特許文献７参照）が適用される。例として、この表中で、Ｑ５はＱ３のＶｏｒｏｎｏｉ拡張であり、Ｑ６はＱ４のＶｏｒｏｎｏｉ拡張である。【００５４】
各コードブックは、ある数のコードベクトルからなり、コードブック中のコードベクトルインデックスは、あるビット数で表現される。このビット数は、下に示す式（１）によって得られる。
【数１】
0006027538.tif 000002
【００５５】
式（１）において、Ｎｂｉｔｓはコードベクトルインデックスを表現するために使用されるビット数、Ｎｃｖはコードブック中のコードベクトル数を表す。【００５６】
なお、コードブックＱ０には、一つのベクトル、すなわち、零ベクトルしかなく、零ベクトルはベクトルの量子化値が０であることを意味する。従って、コードベクトルインデックスのために必要なビットはない。【００５７】
スプリットマルチレート格子ＶＱによって生成された量子化パラメータの３つのセット、すなわち、グローバルゲインのインデックス、コードブックの指示値及びコードベクトルのインデックスからビットストリームを形成するには、二つの方法が考えられる。第１のビットストリーム形成方法を図７に例示し、第２のビットストリーム形成方法を図８に例示する。ここでは、入力スペクトルを６ブロック（ｖ０〜ｖ５）に分割する場合を示す。【００５８】
第１のビットストリーム形成方法では、まず、グローバルゲインＧがスカラー量子化器（図７のＱ）によって量子化される。分割された各ブロックについて、Ｓ（ｆ）／Ｇがマルチレート格子ベクトル量子化器（図７のＶＱ）によって量子化される。図７に示すように、量子化されたグローバルゲインのインデックスは、ビットストリームの先頭の第１領域に配置される。次に、コードブック指示値（Ｃｂ１〜Ｃｂ５）がブロック番号の若い順に先頭側から第２領域に配置され、第２領域に後続して、コードベクトルのインデックスがブロック番号の若い順に先頭側から第３領域に配置される。【００５９】
第２のビットストリーム形成方法では、まず、グローバルゲインＧがスカラー量子化器（図８のＱ）によって量子化される。分割された各ブロックについて、Ｓ（ｆ）／Ｇがマルチレート格子ベクトル量子化器（図８のＶＱ）によって量子化される。図８に示すように、量子化されたグローバルゲインのインデックスは、ビットストリームの先頭の第１領域に配置される。次に、各ベクトルのコードブック指示値とコードベクトルインデックスとを組みとして、第１領域に後続して、ベクトル毎に第２〜第７領域に配置される。
【先行技術文献】
【非特許文献】
【００６０】
【非特許文献１】Karl Heinz Brandenburg, “MP3 and AAC Explained”, AES 17th International Conference, Florence, Italy, September 1999
【非特許文献２】Lefebvre, et al., “High quality coding of wideband audio signals using transform coded excitation (TCX)”, IEEE International Conference on Acoustics, Speech, and Signal Processing, vol. 1, pp. I/193-I/196, Apr. 1994
【非特許文献３】ITU-T Recommendation G.729.1 (2007) “G.729-based embedded variable bit-rate coder: An 8-32kbit/s scalable wideband coder bitstream interoperable with G.729”
【非特許文献４】T. Vaillancourt et al, “ITU-T EV-VBR: A Robust 8-32 kbit/s Scalable Coder for Error Prone Telecommunication Channels”, in Proc. Eusipco, Lausanne, Switzerland, August 2008
【非特許文献５】M. Xie and J.-P. Adoul, “Embedded algebraic vector quantization (EAVQ) with application to wideband audio coding”, IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), Atlanta, GA, U.S.A, 1996, vol. 1, pp. 240-243
【非特許文献６】3GPP TS 26.290 “Extended AMR Wideband Speech Codec (AMR-WB+)”
【非特許文献７】S. Ragot, B. Bessette and R. Lefebvre, “Low-complexity Multi-Rate Lattice Vector Quantization with Application to Wideband TCX Speech Coding at 32kbit/s”, Proc. IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), Montreal, QC, Canada, May, 2004, vol. 1, pp. 501-504
【発明の概要】
【発明が解決しようとする課題】
【００６１】
入力スペクトルは、通常、全てのブロック（以下、適宜、「サブベクトル」という）が均等にエネルギーを有するものではなく、あるサブベクトルにエネルギーが集中する傾向にある。一例として、図９に示したスペクトルの例では、８個のサブベクトルのうち、サブベクトルｖ３が最大のエネルギーを有する。このため、スプリットマルチレート格子ＶＱの処理によって、他のサブベクトルと比べてサブベクトルｖ３のコードブックは、最大のコードブック番号（Ｑｎの整数ｎはコードブック番号と呼ばれる）を有することになる。【００６２】
図６のコードブックリストに示したように、コードブック番号が大きければ大きいほど、コードブック指示値が使用するビット数（つまり、コードブック指示値を表現するために使用されるビット数）はより多くなる。つまり、図９の例では、サブベクトルｖ３はコードブック指示値に最も多くのビット数（以下、適宜、「最多ビット数」という）を使用する。より大きなコードブック番号を有するコードブックについては、コードブック指示値が使用するビット数を減らすことが望ましい。【００６３】
非特許文献７では、サブベクトル数Ｎｓｖ＝８と、使用可能なビット数Ｂｉｔｓａｖａｉｌａｂｌｅ＝１３２ビットの条件下で、スプリットマルチレート格子ＶＱがＴＣＸ音声コーデック用に使用されている。このとき、最大のコードブック番号として１１が測定されたことが言及されている。【００６４】
図９に示したスペクトルの例において、サブベクトルｖ３がＱ１１のコードブックを使用し、サブベクトルｖ４がＱ３のコードブックを使用し、他のすべてのサブベクトルがＱ２のコードブックを使用すると仮定する。この場合の使用ビット数を図１０にまとめる。【００６５】
図１０に示すとおり、サブベクトルｖ３のコードブック指示値は１１ビットを使用し、サブベクトルｖ４のコードブック指示値は３ビットを使用する。その他のサブベクトルのコードブック指示値は２ビットを使用する。すなわち、サブベクトルｖ３のコードブック指示値は、サブベクトルｖ１（ｖ２、ｖ５、ｖ６、ｖ７、またはｖ８）のコードブック指示値よりも５倍も多いビット数を使用することが分かる。【００６６】
非特許文献７に開示の技術では、コードブック指示値とコードベクトルインデックスは２進数に直接変換され、ビットストリームが形成される。従って、全てのベクトルで使用される総ビット数は、次式（２）によって計算される。
【数２】
0006027538.tif 000003
【００６７】
式（２）において、Ｂｉｔｓｔｏｔａｌは総使用ビット数、Ｂｉｔｓｇａｉｎ＿ｑはグローバルゲインの量子化に必要なビット数、Ｂｉｔｓｃｂ＿ｉｎｄｉｃａｔｉｏｎはベクトル当たりのコードブック指示値の使用ビット数、Ｂｉｔｓｃｖ＿ｉｎｄｅｘはベクトル当たりのコードベクトルインデックスの使用ビット数、Ｎはスペクトル全体における総ベクトル数である。【００６８】
図９に示したサブベクトルｖ３のように、より大きなコードブック番号を有するコードブックが選択された場合、このコードブック指示値はより多くのビット数を使用するため、ビットレートが増大するという問題がある。【００６９】
本発明の目的は、スプリットマルチレート格子ベクトル量子化において、より大きなコードブック番号を有するコードブックが選択された場合であっても、ビットレートを低減することができる音声符号化装置、音声復号装置、音声符号化方法及び音声復号方法を提供することである。
【課題を解決するための手段】
【００７０】
本発明の音声符号化装置は、時間領域の入力信号を周波数領域の信号に変換する時間−周波数領域変換手段と、前記周波数領域の信号を複数のサブベクトルに分割し、分割された各サブベクトルのスペクトル係数を量子化するベクトル量子化手段と、前記量子化によって得られた、全ての前記サブベクトルのコードブック指示値を符号化する符号化手段と、全ての前記コードブック指示値のうち、最多使用ビット数のコードブック指示値となるサブベクトルの位置を特定する特定手段と、前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定手段と、前記最多使用ビット数のコードブック指示値が前記符号化手段によって符号化されたビット数を実際値とし、推定された前記最多使用ビット数のコードブック指示値のビット数を推定値として、前記実際値と前記推定値との差分を算出する差分算出手段と、特定された前記サブベクトルの位置情報と算出された前記差分情報を符号化するパラメータ符号化手段と、を具備する構成を採る。【００７１】
本発明の音声復号装置は、上記音声符号化装置より、符号化された前記位置情報及び前記差分情報を取得して、復号するパラメータ復号手段と、前記音声符号化装置より、前記最多使用ビット数のコードブック指示値以外の符号化されたコードブック指示値を取得して、復号するコードブック指示値復号手段と、前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定手段と、推定された前記最多使用ビット数のコードブック指示値が使用するビット数と復号された前記差分情報とを加算して前記最多使用ビット数のコードブック指示値を算出する加算手段と、前記復号された位置情報、前記最多使用ビット数のコードブック指示値以外の復号されたコードブック指示値、及び、算出された前記最多使用ビット数のコードブック指示値を用いて、全てのコードブック指示値を生成するコードブック指示値生成手段と、生成された前記全てのコードブック指示値を用いて、各サブベクトルのスペクトル係数を逆量子化するベクトル逆量子化手段と、逆量子化された前記スペクトル係数を時間領域に変換する周波数−時間変換手段と、を具備する構成を採る。【００７２】
本発明の音声符号化方法は、時間領域の入力信号を周波数領域の信号に変換する時間−周波数領域変換工程と、前記周波数領域の信号を複数のサブベクトルに分割し、分割された各サブベクトルのスペクトル係数を量子化するベクトル量子化工程と、前記量子化によって得られた、全ての前記サブベクトルのコードブック指示値を符号化する符号化工程と、全ての前記コードブック指示値のうち、最多使用ビット数のコードブック指示値となるサブベクトルの位置を特定する特定工程と、前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定工程と、前記最多使用ビット数のコードブック指示値が前記符号化工程において符号化されたビット数を実際値とし、推定された前記最多使用ビット数のコードブック指示値のビット数を推定値として、前記実際値と前記推定値との差分を算出する差分算出工程と、特定された前記サブベクトルの位置情報と算出された前記差分情報を符号化するパラメータ符号化工程と、を具備するようにした。【００７３】
本発明の音声復号方法は、上記音声符号化方法により符号化された前記位置情報及び前記差分情報を復号するパラメータ復号工程と、上記音声符号化方法により符号化された、前記最多使用ビット数のコードブック指示値以外のコードブック指示値を復号するコードブック指示値復号工程と、前記入力信号の送信単位において利用可能な総ビット数と、前記最多使用ビット数のコードブック指示値以外のコードブック指示値の使用ビット数とに基づいて、前記最多使用ビット数のコードブック指示値が使用するビット数を推定する推定工程と、推定された前記最多使用ビット数のコードブック指示値が使用するビット数と復号された前記差分情報とを加算して前記最多使用ビット数のコードブック指示値を算出する加算工程と、前記復号された位置情報、前記最多使用ビット数のコードブック指示値以外の復号されたコードブック指示値、及び、算出された前記最多使用ビット数のコードブック指示値を用いて、全てのコードブック指示値を生成するコードブック指示値生成工程と、生成された前記全てのコードブック指示値を用いて、各サブベクトルのスペクトル係数を逆量子化するベクトル逆量子化工程と、逆量子化された前記スペクトル係数を時間領域に変換する周波数−時間変換工程と、を具備するようにした。
【発明の効果】
【００７４】
本発明によれば、スプリットマルチレート格子ベクトル量子化において、より大きなコードブック番号を有するコードブックが選択された場合であっても、ビットレートを低減することができる。
【図面の簡単な説明】
【００７５】
【図１】変換コーデックの簡略な構成を示す図
【図２】ＴＣＸコーデックの簡略な構成を示す図
【図３】ＣＥＬＰと変換符号化の階層符号化の簡略な構成を示す図
【図４】ＴＣＸコーデックにおいてスプリットマルチレート格子ＶＱを利用した簡略な構成を示す図
【図５】スプリットマルチレート格子ＶＱの処理を示すブロック図
【図６】スプリットマルチレート格子ＶＱのコードブックリストの概要を示す図
【図７】第１のビットストリーム形成方法を示す図
【図８】第２のビットストリーム形成方法を示す図
【図９】入力信号のスペクトルを示す図
【図１０】図９に示したスペクトルに適用するコードブック及び使用ビット数を示す図
【図１１】本発明の実施の形態１に係る符号化装置の構成を示すブロック図
【図１２】図１１に示したコードブック指示値変換部の内部構成を示すブロック図
【図１３】サブベクトルの位置情報と位置情報符号化コードとの対応関係を示す図
【図１４】差分情報と差分情報符号化コードとの対応関係を示す図
【図１５】本発明の実施の形態１に係る復号装置の構成を示すブロック図
【図１６】図１５に示したコードブック指示値逆変換部の内部構成を示すブロック図
【図１７】図９に示したスペクトルに適用するコードブック及びコードブック指示値を示す図
【図１８】本発明の実施の形態１に係るスプリットマルチレート格子ＶＱをＣＥＬＰと変換符号化の階層的符号化に適用したコーデックの構成を示すブロック図
【図１９】本発明の実施の形態１に係るスプリットマルチレート格子ＶＱをＴＣＸコーデックに適用したコーデックの構成を示すブロック図
【図２０】本発明の実施の形態２に係るコードブック指示値変換部の内部構成を示すブロック図
【図２１】図２０に示した第１変換処理部の内部構成を示すブロック図
【図２２】図２０に示した第２変換処理部の内部構成を示すブロック図
【図２３】本発明の実施の形態２に係るコードブック指示値逆変換部の内部構成を示すブロック図
【図２４】図２３に示した第１逆変換処理部の内部構成を示すブロック図
【図２５】図２３に示した第２逆変換処理部の内部構成を示すブロック図
【発明を実施するための形態】
【００７６】
以下、本発明の実施の形態について、図面を参照して詳細に説明する。
【００７７】
（実施の形態１）
図１１は、本発明の実施の形態１に係る符号化装置１００の構成を示すブロック図である。図１１において、時間−周波数変換部１０１は、ＤＦＴまたはＭＤＣＴなどの時間−周波数変換方式によって、時間領域の入力信号Ｓ（ｎ）を周波数領域の入力信号Ｓ（ｆ）に変換し、心理音響モデル分析部１０２及びスプリットマルチレート格子ＶＱ部１０３に周波数領域の入力信号Ｓ（ｆ）を出力する。
【００７８】
心理音響モデル分析部１０２は、時間−周波数変換部１０１から出力された周波数領域の入力信号Ｓ（ｆ）に対して心理音響モデル分析を行い、マスキング曲線を取得する。取得されたマスキング曲線はスプリットマルチレート格子ＶＱ部１０３に出力される。
【００７９】
スプリットマルチレート格子ＶＱ部１０３は、量子化ノイズが不可聴となるように、心理音響モデル分析部１０２から出力されたマスキング曲線に従って、周波数領域の入力信号Ｓ（ｆ）にスプリットマルチレート格子ＶＱを適用する。得られたグローバルゲイン及びコードベクトルインデックスは多重化部１０５に出力される。また、スプリットマルチレート格子ＶＱ部１０３は、図６のコードブックリストを参照し、全てのサブベクトルの各コードブック指示値を算定し、算定した各コードブック指示値をコードブック指示値変換部１０４に出力する。
【００８０】
コードブック指示値変換部１０４は、スプリットマルチレート格子ＶＱ部１０３から出力された全てのサブベクトルの各コードブック指示値を用いて以下のステップ１〜ステップ３の処理を行う。
【００８１】
（ステップ１）スプリットマルチレート格子ＶＱ部１０３から出力された各コードブック指示値のうち、コードブック指示値が最多ビット数を使用するサブベクトルの位置を特定して、位置情報を符号化し、コードブック指示値が最多ビット数を使用するサブベクトルを除くすべてのサブベクトルのコードブック指示値を符号化。
（ステップ２）コードブック指示値が最多ビット数を使用するサブベクトルのコードブック指示値を、それを除くすべてのサブベクトルのコードブック指示値を用いて推定。
（ステップ３）最多ビット数を使用するコードブック指示値の実際の使用ビット数と（ステップ２）にて推定された使用ビット数との差分を符号化。
【００８２】
これによって得られた、位置情報符号化コード、コードブック指示値符号化コード及び差分情報符号化コードが多重化部１０５に出力される。なお、コードブック指示値変換部１０４の詳細については後述する。
【００８３】
多重化部１０５は、スプリットマルチレート格子ＶＱ部１０３から出力されたグローバルゲイン及びコードベクトルインデックスと、コードブック指示値変換部１０４から出力された位置情報符号化コード、コードブック指示値符号化コード及び差分情報符号化コードとを多重化し、多重化したビットストリーム情報を復号装置へ送信する。
【００８４】
次に、上述したコードブック指示値変換部１０４の内部構成について図１２を用いて説明する。スプリットマルチレート格子ＶＱ部１０３から出力されたＮ個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ＝１〜Ｎ）は、コードブック指示値変換部１０４のサブベクトル特定部１２１及びコードブック指示値分離部１２３に入力される。
【００８５】
サブベクトル特定部１２１は、入力されたＮ個のコードブック指示値ｃｂｖｉを参照し、コードブック指示値が最多ビット数を使用するサブベクトルを特定し、特定したサブベクトルの位置情報Ｐｃｂｍａｘを位置情報符号化部１２２及びコードブック指示値分離部１２３に出力する（上記ステップ１に相当）。
【００８６】
位置情報符号化部１２２は、サブベクトル特定部１２１から出力された位置情報Ｐｃｂｍａｘを、図１３の表を用いて位置情報符号化コードを生成し、多重化部１０５に出力する。ここで、図１３は、サブベクトルの位置情報と位置情報符号化コードとを対応付けた表である（上記ステップ１に相当）。
【００８７】
コードブック指示値分離部１２３は、入力されたＮ個のコードブック指示値ｃｂｖｉと、サブベクトル特定部１２１から出力された位置情報Ｐｃｂｍａｘとを参照し、位置情報Ｐｃｂｍａｘに位置するサブベクトルのコードブック指示値ｃｂｍａｘ（またはｃｂｖｉ（ｉ＝Ｐｃｂｍａｘ））を差分算出部１２５に出力する。そして、それ以外のＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｃｂｍａｘ）をコードブック指示値推定部１２４に出力すると共に、コードブック指示値符号化コードとして多重化部１０５に出力する（上記ステップ１に相当）。
【００８８】
コードブック指示値推定部１２４は、最多ビット数を使用する（以下、適宜、「最多使用ビット数」という）コードブック指示値の使用ビット数を、それ以外のＮ−１個のコードブック指示値を用いて推定し、使用ビット数推定値ｃｂ’ｍａｘを生成し、差分算出部１２５に出力する。
【００８９】
最多使用ビット数のコードブック指示値の使用ビット数の推定は、次式（３）に従って実行される。
【数３】
0006027538.tif 000004
【００９０】
式（３）において、右辺のΣは、ｉがＰｃｂｍａｘ以外の全てのＢｉｔｓｃｂｖｉの総和を算出するものとし、ｃｂ’ｍａｘは最多使用ビット数のコードブック指示値の使用ビット数推定値、Ｂｉｔｓａｖａｉｌａｂｌｅは入力信号の送信単位において利用可能な総ビット数、Ｂｉｔｓｃｂｖｉはサブベクトルｖｉのコードブック指示値の使用ビット数をそれぞれ表す。ただし、当該フレームで利用可能な総ビット数は、予め規定の値として設定されていてもよいし、もしくは、規定の値からグローバルゲインのような直前までの符号化処理で使用されたビット数を減じた値として、符号化装置及び復号装置で参照可能な数値であってもよい。
【００９１】
このように、コードブック指示値推定部１２４は、式（３）に示したとおり、総ビット数からＮ−１個のサブベクトルのコードブック指示値で使用されるビット数を減じることにより、最多使用ビット数のコードブック指示値の使用ビット数を推定する（上記ステップ２に相当）。
【００９２】
なお、式（３）において、５で除しているのは、図６に示されるように、コードブック指示値の使用ビット数とコードベクトルインデックスの使用ビット数は１：４の関係にあるためである。よって、減算後のビット数（この時点では、コードブック指示値とコードベクトルインデックスの使用ビット数の両者を含む）を５で除することにより、コードブック指示値の使用ビット数の推定値を算出している。
【００９３】
差分算出部１２５は、コードブック指示値分離部１２３から出力された最多使用ビット数のコードブック指示値の使用ビット数（実際値ｃｂｍａｘ）に対して、コードブック指示値推定部１２４から出力された使用ビット数推定値（ｃｂ’ｍａｘ）を減じて（ｃｂｍａｘ−ｃｂ’ｍａｘ）、差分情報ｃｂｄｉｆｆを算出し、差分符号化部１２６に出力する。
【００９４】
差分符号化部１２６は、図１４に示す表に基づいて、差分情報ｃｂｄｉｆｆから差分情報符号化コードを生成し、多重化部１０５に出力する（上記ステップ３に相当）。ここで、図１４は、差分情報ｃｂｄｉｆｆと差分情報符号化コードとを対応付けた表である。図１４から分かるように、差分情報ｃｂｄｉｆｆは０以下の負の整数である。これは、最多使用ビット数のコードブック指示値の使用ビット数推定値は、利用可能な全てのビット数が量子化に使用されることを前提として算出されるためであり、量子化において利用可能なビット数より多いビット数が使用されることはない。すなわち、推定値は最大可能値であり、推定値が実際値より小さくなることはあり得ない。
【００９５】
図１５は、本発明の実施の形態１に係る復号装置２００の構成を示すブロック図である。図１５において、符号化装置１００から送信された全てのビットストリーム情報が多重分離部２０１に入力される。
【００９６】
多重分離部２０１は、入力されたビットストリーム情報からグローバルゲイン、コードベクトルインデックス、位置情報符号化コード、コードブック指示値符号化コード及び差分情報符号化コードに多重分離する。多重分離されたグローバルゲイン及びコードベクトルインデックスはスプリットマルチレート格子逆ＶＱ部２０３に出力され、多重分離された位置情報符号化コード、コードブック指示値符号化コード及び差分情報符号化コードはコードブック指示値逆変換部２０２に出力される。
【００９７】
コードブック指示値逆変換部２０２は、多重分離部２０１から出力された位置情報符号化コード、コードブック指示値符号化コード及び差分情報符号化コードを用いて、以下のステップ４〜ステップ８の処理を行う。
【００９８】
（ステップ４）最多使用ビット数のコードブック指示値となるサブベクトルの位置を復号。
（ステップ５）その他全てのサブベクトルのコードブック指示値を復号。
（ステップ６）最多使用ビット数のコードブック指示値の使用ビット数を、それを除くすべてのサブベクトルのコードブック指示値の使用ビット数を用いて推定。
（ステップ７）差分情報符号化コードを復号。
（ステップ８）（ステップ６）で推定された使用ビット数推定値と（ステップ７）で復号された差分情報とを加算することによって、最多使用ビット数のコードブック指示値を算出。
【００９９】
復号されたコードブック指示値はスプリットマルチレート格子逆ＶＱ部２０３に出力される。なお、コードブック指示値逆変換部２０２の詳細については後述する。
【０１００】
スプリットマルチレート格子逆ＶＱ部２０３は、多重分離部２０１から出力されたグローバルゲイン及びコードベクトルインデックスと、コードブック指示値逆変換部２０２から出力されたコードブック指示値に対して、スプリットマルチレート格子逆ＶＱを適用し、周波数領域の信号Ｓ〜（ｆ）に復号する。復号された周波数領域の信号Ｓ〜（ｆ）は周波数−時間変換部２０４に出力される。
【０１０１】
周波数−時間変換部２０４は、ＩＤＦＴまたはＩＭＤＣＴなどの周波数−時間変換方式によって、スプリットマルチレート格子逆ＶＱ部２０３から出力された周波数領域の信号Ｓ〜（ｆ）を時間領域の信号Ｓ〜（ｎ）に変換する。
【０１０２】
次に、上述したコードブック指示値逆変換部２０２の内部構成について図１６を用いて説明する。多重分離部２０１から出力された位置情報符号化コードは、位置情報復号部２２１に入力され、コードブック指示値符号化コードはコードブック指示値推定部２２２及びコードブック指示値生成部２２５に入力され、差分情報符号化コードは差分復号部２２３に入力される。
【０１０３】
位置情報復号部２２１は、多重分離部２０１から出力された位置情報符号化コードに対応する位置情報Ｐｃｂｍａｘを、図１３の表を参照して生成し、コードブック指示値生成部２２５に出力する（上記ステップ４に相当）。
【０１０４】
多重分離部２０１から出力されたコードブック指示値符号化コードは、位置情報Ｐｃｂｍａｘのサブベクトルを除いたＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｃｂｍａｘ）を表しており、これらはコードブック指示値生成部２２５及びコードブック指示値推定部２２２に与えられる（上記ステップ５に相当）。
【０１０５】
コードブック指示値推定部２２２は、多重分離部２０１から出力されたＮ−１個のコードブック指示値ｃｂｖｉ（ｉ≠Ｐｃｂｍａｘ）を用いて、位置情報Ｐｃｂｍａｘに位置するサブベクトルのコードブック指示値の使用ビット数推定値ｃｂ’ｍａｘを算出し、加算部２２４に出力する（上記ステップ６に相当）。コードブック指示値推定部２２２は、図１２の符号化装置側と同様の処理を行うので、ここではその詳細な説明は省略する。
【０１０６】
差分復号部２２３は、多重分離部２０１から出力された差分情報符号化コードに対応する差分情報ｃｂｄｉｆｆを、図１４の表を参照して生成し、加算部２２４に出力する（上記ステップ７に相当）。
【０１０７】
加算部２２４は、コードブック指示値推定部２２２から出力された使用ビット数推定値ｃｂ’ｍａｘに差分復号部２２３から出力された差分情報ｃｂｄｉｆｆを加算し、位置情報Ｐｃｂｍａｘに位置するサブベクトルのコードブック指示値ｃｂｍａｘ（またはｃｂｖｉ（ｉ＝Ｐｃｂｍａｘ））を生成し、コードブック指示値生成部２２５に出力する（上記ステップ８に相当）。
【０１０８】
コードブック指示値生成部２２５は、位置情報復号部２２１から出力された位置情報Ｐｃｂｍａｘ、多重分離部２０１から出力されたＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｃｂｍａｘ）、及び、加算部２２４から出力された位置情報Ｐｃｂｍａｘに位置するサブベクトルのコードブック指示値ｃｂｍａｘを用いて、コードブック指示値ｃｂｍａｘが位置情報Ｐｃｂｍａｘに配置されるように、全てのコードブック指示値ｃｂｖｉ（ｉ＝１〜Ｎ）を生成し、スプリットマルチレート格子逆ＶＱ部２０３に出力する。
【０１０９】
次に、図１２に示したコードブック指示値変換部１０４の動作について、図９に示したスペクトルを符号化する場合を例にとって、具体的に説明する。
【０１１０】
図９のスペクトルを符号化装置１００の時間−周波数変換部１０１において変換された周波数領域の入力信号Ｓ（ｆ）とした場合、サブベクトル特定部１２１には、図１７に示すコードブック指示値ｃｂｖｉ（ｉ＝１〜８）が入力されると仮定する。このとき、サブベクトル特定部１２１では、コードブック指示値が最多ビット数を使用するサブベクトルｖ３が特定され、サブベクトルｖ３が位置情報Ｐｃｂｍａｘとして位置情報符号化部１２２及びコードブック指示値分離部１２３に出力される。
【０１１１】
位置情報符号化部１２２では、サブベクトル特定部１２１から出力された位置情報Ｐｃｂｍａｘが、図１３の表を用いて符号化される。位置情報Ｐｃｂｍａｘはサブベクトルｖ３を示すので、図１３の表より、位置情報符号化コード「０１０」が生成される。
【０１１２】
コードブック指示値分離部１２３では、入力された図１７に示すコードブック指示値ｃｂｖｉと、サブベクトル特定部１２１によって特性されたサブベクトルｖ３を示す位置情報Ｐｃｂｍａｘとに基づいて、サブベクトルｖ３のコードブック指示値、すなわち、図１７によると、「１１１１１１１１１１０」（つまり、サブベクトルｖ３のコードブック指示値の使用ビット数は１１ビット）が差分算出部１２５に出力される。また、サブベクトルｖ３以外のサブベクトルｖ１，ｖ２，ｖ４〜ｖ８のコードブック指示値、すなわち、図１７によると、「１０」、「１０」、「１１０」、「１０」、「１０」、「１０」、「１０」がコードブック指示値符号化コードとして多重化部１０５へ出力される。
【０１１３】
コードブック指示値推定部１２４では、最多使用ビット数のコードブック指示値の使用ビット数、ここでは、サブベクトルｖ３のコードブック指示値の使用ビット数が推定される。具体的には、入力信号の送信単位において利用可能な総ビット数を１３２ビットとした場合、サブベクトルｖ３以外のサブベクトルそれぞれのサブベクトル当たりの使用ビット数を図１０から求めると、式（３）に従って、（１３２−１０−１０−１５−１０−１０−１０−１０）／５-≒１１が使用ビット数推定値として求められる。
【０１１４】
差分算出部１２５では、コードブック指示値分離部１２３から出力されたサブベクトルｖ３のコードブック指示値の使用ビット数（ここでは、１１ビット）から、コードブック指示値推定部１２４から出力された使用ビット数推定値（ここでは、１１ビット）が減算され、差分情報ｃｂｄｉｆｆ（ここでは、１１−１１＝０）が算出される。
【０１１５】
差分符号化部１２６では、差分算出部１２５によって算出された差分情報ｃｂｄｉｆｆが０であることから、図１４の表より、差分情報符号化コード「０」が生成される。
【０１１６】
このようにして生成された位置情報符号化コード「０１０」、コードブック指示値符号化コード「１０」、「１０」、「１１０」、「１０」、「１０」、「１０」、「１０」、差分情報符号化コード「０」が多重化部１０５において多重され、復号装置２００に出力される。
【０１１７】
上述したように、図９に示した例では、サブベクトルｖ３に適用される最大のコードブック番号が１１であり、このコードブック指示値に使用されるビット数は１１ビットである。
【０１１８】
位置情報Ｐｃｂｍａｘの符号化では、固定のビット数（Ｂｉｔｓｐｏｓｉｔｉｏｎ＿ｃｂｍａｘ＝ｌｏｇ２（Ｎｓｖ））が使用される。ここでは、サブベクトルの数はＮｓｖ＝８個であるので、位置情報の符号化では、３ビットが使用される。また、差分情報の使用ビット数は、上述したとおり１ビットである。
【０１１９】
すなわち、図９で示した例では、１１ビットで通知されるべきサブベクトルｖ３のコードブック指示値が、３ビットの位置情報と１ビットの差分情報によって通知することができ、７ビットの削減効果が得られる。このようにビットを削減していても情報が失われることが無いため、元のコードブック指示値を完全に復元することができる。
【０１２０】
このように、実施の形態１によれば、所定数のサブベクトルに分割した入力信号のスペクトルのうち、最多ビット数を使用するコードブック指示値を、実際の使用ビット数及び使用ビット数推定値に基づく差分情報と、当該コードブック指示値を使用するサブベクトルの位置情報とに変換することにより、最多ビット数を使用するコードブック指示値の使用ビット数を削減することができ、ビットレートを低減することができる。
【０１２１】
なお、本実施の形態に係るスプリットマルチレート格子ＶＱをＣＥＬＰと変換符号化の階層的符号化に適用してもよく、この場合のコーデックの構成を図１８に示す。図１８（ａ）は、図３（ａ）の符号化器において、量子化部５４をスプリットマルチレート格子ＶＱ部１０３及びコードブック指示値変換部１０４に変更し、図１８（ｂ）は、図３（ｂ）の復号器において、逆量子化部６２をコードブック指示値逆変換部２０２及びスプリットマルチレート格子逆ＶＱ部２０３に変更した点が異なる。
【０１２２】
また、本実施の形態に係るスプリットマルチレート格子ＶＱをＴＣＸコーデックに適用してもよく、この場合のコーデックの構成を図１９に示す。図１９（ａ）は、図４（ａ）の符号化器において、スプリットマルチレート格子ＶＱ部７６をスプリットマルチレート格子ＶＱ部１０３及びコードブック指示値変換部１０４に変更し、図１９（ｂ）は、図４（ｂ）の復号器において、スプリットマルチレート格子逆ＶＱ部８２をコードブック指示値逆変換部２０２及びスプリットマルチレート格子逆ＶＱ部２０３に変更した点が異なる。
【０１２３】
（実施の形態２）
実施の形態１では、入力信号のスペクトルを所定数のサブベクトルに分割し、特定のサブベクトルにエネルギーが集中する場合を想定して説明したが、特定のサブベクトルにエネルギーが集中しない場合も考えられ、このような場合に実施の形態１では、使用ビット数が増大してしまうことがある。
【０１２４】
そこで、本発明の実施の形態２では、特定のサブベクトルにエネルギーが集中しない場合においても、使用ビット数を削減する場合について説明する。
【０１２５】
本発明の実施の形態２に係る符号化装置及び復号装置の構成は、実施の形態１の図１１及び図１５に示した構成のうち、コードブック指示値変換部１０４及びコードブック指示値逆変換部２０２の機能が異なるのみなので、図１１及び図１５を援用し、重複する説明は省略し、異なる機能について説明する。
【０１２６】
図２０は、本発明の実施の形態２に係るコードブック指示値変換部３００の内部構成を示すブロック図である。図２０において、図１２と異なる点は、判定部３０１、第１変換処理部３０２及び第２変換処理部３０３である。
【０１２７】
判定部３０１は、コードブック指示値分離部１２３から出力された最多使用ビットのサブベクトルのコードブック指示値の使用ビット数と、所定の閾値とを比較する。比較の結果、コードブック指示値の使用ビット数が閾値よりも大きい場合、判定部３０１は、スイッチを第１変換処理部３０２に切り替え、第１変換処理部３０２と接続する。一方、コードブック指示値の使用ビット数が閾値以下の場合、判定部３０１は、スイッチを第２変換処理部３０３に切り替え、第２変換処理部３０３と接続する。
【０１２８】
第１変換処理部３０２は、判定部３０１と接続された場合、位置情報Ｐｃｂｍａｘ、最多使用ビット数のサブベクトルのコードブック指示値ｃｂｍａｘ（またはｃｂｖｉ（ｉ＝Ｐｃｂｍａｘ））及び最多使用ビット数のサブベクトル以外のＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｃｂｍａｘ）が入力され、位置情報符号化コード、コードブック指示値符号化コード及び差分情報符号化コードを生成し、出力する。第１変換処理部３０２の内部構成を図２１に示す。図２１において、図１２と共通する符番は図１２の構成と同一なので、重複する説明は省略する。
【０１２９】
再度、図２０を参照するに、第２変換処理部３０３は、判定部３０１と接続された場合、Ｎ個のコードブック指示値ｃｂｖｉ（ｉ＝１〜Ｎ）が入力され、コードブック指示値符号化コード、差分情報符号化コードを生成し、出力する。第２変換処理部３０３の内部構成を図２２に示す。図２２において、図１２と共通する符番は図１２の構成と同一なので、重複する説明は省略する。
【０１３０】
図２２において、固定コードブック指示値分離部３２１は、入力されたＮ個のコードブック指示値ｃｂｖｉ（ｉ＝１〜Ｎ）に対し、予め定められた位置Ｐｆｉｘのサブベクトルのコードブック指示値ｃｂｆｉｘ（またはｃｂｖｉ（ｉ＝Ｐｆｉｘ））とそれ以外のＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｆｉｘ）とに分離し、前者を差分算出部３２３に出力する。そして、後者をコードブック指示値推定部３２２に出力すると共に、コードブック指示値符号化コードとして多重化部１０５に出力する。
【０１３１】
コードブック指示値推定部３２２は、固定コードブック指示値分離部３２１から出力されたＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｆｉｘ）を用いて、位置Ｐｆｉｘのサブベクトルのコードブック指示値ｃｂｆｉｘ（またはｃｂｖｉ（ｉ＝Ｐｆｉｘ））の使用ビット推定値ｃｂ’ｆｉｘを次式（４）に従って算出し、差分算出部３２３に出力する。
【数４】
0006027538.tif 000005
【０１３２】
式（４）において、ｃｂ’ｆｉｘは、位置Ｐｆｉｘのサブベクトルのコードブック指示値の使用ビット推定値、Ｂｉｔｓａｖａｉｌａｂｌｅは利用可能な総ビット数、Ｂｉｔｓｃｂｖｉはサブベクトルｖｉのコードブック指示値の使用ビット数を表す。
【０１３３】
差分算出部３２３は、固定コードブック指示値分離部３２１から出力された位置Ｐｆｉｘのサブベクトルのコードブック指示値ｃｂｆｉｘ（またはｃｂｖｉ（ｉ＝Ｐｆｉｘ））に対して、コードブック指示値推定部３２２から出力された使用ビット数推定値ｃｂ’ｆｉｘを減じて（ｃｂｆｉｘ−ｃｂ’ｆｉｘ）、差分情報ｃｂｄｉｆｆを算出し、差分符号化部１２６に出力する。
【０１３４】
図２３は、本発明の実施の形態２に係るコードブック指示値逆変換部４００の内部構成を示すブロック図である。図２３において、図１６と異なる点は、判定部４０１、第１逆変換処理部４０２及び第２逆変換処理部４０３である。
【０１３５】
判定部４０１は、加算部２２４から出力されたコードブック指示値が入力される。このコードブック指示値は、最多使用ビット数のサブベクトル、もしくは予め定められた固定の位置のサブベクトル、のいずれかのサブベクトルのコードブック指示値である。判定部４０１は、このコードブック指示値の使用ビット数と、所定の閾値（符号化装置と同じ閾値）とを比較する。比較の結果、コードブック指示値の使用ビット数が閾値よりも大きい場合、判定部４０１は、コードブック指示値が最多使用ビット数のサブベクトルのコードブック指示値であると判定し、スイッチを第１逆変換処理部４０２に切り替え、第１逆変換処理部４０２と接続する。一方、コードブック指示値の使用ビット数が閾値以下の場合、判定部４０１は、コードブック指示値が予め定められた固定の位置のサブベクトルのコードブック指示値であると判断し、スイッチを第２逆変換処理部４０３に切り替え、第２逆変換処理部４０３と接続する。
【０１３６】
第１逆変換処理部４０２は、判定部４０１と接続された場合、位置情報Ｐｃｂｍａｘ、最多使用ビット数のサブベクトルのコードブック指示値ｃｂｍａｘ（またはｃｂｖｉ（ｉ＝Ｐｃｂｍａｘ））及び最多使用ビット数のサブベクトル以外のＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｃｂｍａｘ）が入力され、これら情報を用いてＮ個のコードブック指示値ｃｂｖｉ（ｉ＝１〜Ｎ）を生成し、出力する。第１逆変換処理部４０２の内部構成を図２４に示す。図２４において、図１６と共通する符番は図１６の構成と同一なので、重複する説明は省略する。
【０１３７】
再度、図２３を参照するに、第２逆変換処理部４０３は、判定部４０１と接続された場合、予め定められた固定の位置のサブベクトルのコードブック指示値ｃｂｆｉｘ（またはｃｂｖｉ（ｉ＝Ｐｆｉｘ））及びこのサブベクトル以外のＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｆｉｘ）が入力され、これら情報を用いてＮ個のコードブック指示値ｃｂｖｉ（ｉ＝１〜Ｎ）を生成し、出力する。第２逆変換処理部４０３の内部構成を図２５に示す。
【０１３８】
図２５において、固定サブベクトル指定部４２１は、固定されたサブベクトルの位置情報Ｐｆｉｘをコードブック指示値生成部４２２に出力する。
【０１３９】
コードブック指示値生成部４２２は、予め定められた固定の位置のサブベクトルのコードブック指示値ｃｂｆｉｘ（またはｃｂｖｉ（ｉ＝Ｐｆｉｘ））、このサブベクトル以外のＮ−１個のサブベクトルのコードブック指示値ｃｂｖｉ（ｉ≠Ｐｆｉｘ）及び固定の位置情報Ｐｆｉｘを用いて、コードブック指示値ｃｂｆｉｘが位置情報Ｐｆｉｘに配置されるように、Ｎ個のコードブック指示値ｃｂｖｉ（ｉ＝１〜Ｎ）を生成し、出力する。
【０１４０】
このように、実施の形態２によれば、最多使用ビット数のコードブック指示値の使用ビット数が所定の閾値以下の場合には、予め固定された位置のサブベクトルのコードブック指示値を実際のビット数及び使用ビット数推定値の差分情報に変換し、符号化装置と復号装置との間で固定された位置情報を共有することにより、位置情報を省略することができる。これにより、特定のサブベクトルにエネルギーが集中しない場合においても、使用ビット数を削減することができる。
【０１４１】
なお、本実施の形態では、固定された位置のサブベクトルとして、入力スペクトルの特性に応じて決定されてもよい。例えば、最初（または最後）のサブベクトルのコードブックがその他のサブベクトルよりも統計的に大きければ、最初（または最後）のサブベクトルを選択してもよい。要は、符号化装置と復号装置との間で固定された位置情報を共有していればよい。
【０１４２】
（実施の形態３）
先行技術においては、コードブック指示値は、コードブックの使用率に応じて設計されておらず、むしろ単純に、図６に示したコードブックリストが広く使用されている。
【０１４３】
しかしながら、様々なビットレート、様々なサブベクトル数などの状況によって、コードブックの使用についての統計は異なってくる。非特許文献７では、ＲＥ８コードブック各々の使用についてのいくつかの統計がテーブル２（ｃ）にまとめられている（コードブック番号ごとの統計（使用率％））。
【０１４４】
この統計情報によれば、図６のコードブックリストの設計は、非特許文献７における利用法に対しては効率的ではないことが分かる。Ｑ０の指示値は最少ビット数（１ビット）を使用するが、その使用率は３％と非常に低い。また、Ｑ２の使用率は最高（２９％）であるが、その使用ビット数は最少ではない。
【０１４５】
従って、一定の条件（同一のビットレート、同一の量子化すべきサブベクトル数）の下では、各コードブックの使用確率に応じて、コードブック指示値にビット数を割り当てることが考えられる。具体的には、使用確率の高いコードブックにはより少ないビット数を割り当て、使用確率の低いコードブックにはより多くのビット数を割り当てるという、ハフマン（Huffman）テーブル設計手法を使用してコードブック指示値を設計することが望ましい。
【０１４６】
（実施の形態４）
本発明の実施の形態４では、実施の形態１によって削減されたビットを利用して、量子化されたベクトルのゲインの精度を向上させる場合について説明する。
【０１４７】
本実施の形態では、実施の形態１で説明したコードブック指示値変換方法によって削減されたビットは、スペクトルをより小さな帯域に分割し、各帯域に「ゲイン補正係数」を付与することによって、グローバルゲインに対してより細かな分解能を与えるために利用可能である。ゲイン補正係数を送信するために、削減されたビットを利用することによって、量子化性能の向上を図り、音質を向上させることができる。
【０１４８】
また、本実施の形態では、実施の形態１で説明したコードブック指示値変換方法を、ステレオまたはマルチチャネル信号の符号化に応用可能である。例えば、コードブック指示値変換方法を副信号の符号化に適用し、削減されたビットを主信号の符号化に利用する。これは、主信号は副信号よりも知覚的により重要であるから、主観的な質の向上を図ることができる。
【０１４９】
また、本実施の形態では、実施の形態１で説明したコードブック指示値変換方法を、複数フレーム単位（または複数サブフレーム単位）でスペクトル係数列を符号化するコーデックに応用可能である。これによれば、次の符号化段階でのスペクトル係数列または他の何らかのパラメータ列を符号化するために、削減されたビットを蓄積して利用することができる。
【０１５０】
さらに、本実施の形態では、実施の形態１で説明したコードブック指示値変換方法によって削減されたビットは、フレーム損失状況において音質を維持できるように、ＦＥＣ（フレーム消失隠蔽）に利用可能である。
【０１５１】
なお、上記各実施の形態では、スプリットマルチレート格子ＶＱを使用するものとして説明したが、本発明はこれに限定されず、その他のスペクトル係数コーディング手法を用いてもよい。当業者は、本発明の精神から逸脱しない範囲で本発明を修正し、適応させることができると考えられる。
【０１５２】
また、上記各実施の形態では、最多使用ビット数の一つのサブベクトルに対してコードブック指示値を推定する方法について説明したが、本発明はこれに限らず、複数のサブベクトルのコードブック指示値を推定する構成にも本発明を適用できる。
【０１５３】
例えば、コードブック指示値の使用ビット数の上位Ｎ個のサブベクトルに対して本発明を適用する場合、該当するＮ個のサブベクトルの位置情報をそれぞれ復号装置に伝送すると共に、利用可能な総ビット数から当該Ｎ個のサブベクトル以外のサブベクトルの使用ビット数を減じた後、これをＮで除した値を当該Ｎ個のサブベクトルのコードブック指示値の使用ビット数推定値として用いる。この使用ビット数推定値は次式（５）によって算出される。
【数５】
0006027538.tif 000006
【０１５４】
式（５）において、ｃｂ’ｍａｘはコードブック指示値の使用ビット数が上位Ｎ個のサブベクトルに対する使用ビット数推定値、Ｂｉｔｓａｖａｉｌａｂｌｅは利用可能な総ビット数、Ｂｉｔｓｃｂｖｉはサブベクトルｖｉのコードブック指示値に使用されるビット数、ＳＮはコードブック指示値の使用ビット数が上位Ｎ個のサブベクトルを表す集合である。
【０１５５】
符号化装置は、上位Ｎ個のサブベクトルのコードブック指示値の実際値から式（５）によって算出された使用ビット数推定値ｃｂ’ｍａｘを減じ、その差分情報を符号化する。
【０１５６】
また、サブベクトルの位置情報を復号装置へ送信する代わりに、予め定められた位置の複数のサブベクトルに対して式（５）に従ってコードブック指示値を推定してもよい。この場合、サブベクトルの位置情報を復号装置に伝送する必要はない。
【０１５７】
上記各実施の形態では、本発明をハードウェアで構成する場合を例にとって説明したが、本発明はハードウェアとの連携においてソフトウェアによって実現することも可能である。
【０１５８】
また、上記各実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されてもよいし、一部又は全てを含むように１チップ化されてもよい。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
【０１５９】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路又は汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Field Programmable Gate Array）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
【０１６０】
さらには、半導体技術の進歩又は派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。バイオ技術の適用等が可能性としてありえる。
【０１６１】
２０１１年１０月２８日出願の特願２０１１−２３７８２０の日本出願に含まれる明細書、図面および要約書の開示内容は、すべて本願に援用される。
【産業上の利用可能性】
【０１６２】
本発明にかかる音声符号化装置、音声復号装置、音声符号化方法及び音声復号方法は、移動通信システムの無線通信端末装置及び無線通信基地局装置、遠隔会議端末装置、ビデオ会議端末装置及びＶｏＩＰ（Voice over Internet Protocol）端末装置等に適用可能である。
【符号の説明】
【０１６３】
１０１  時間−周波数変換部
１０２  心理音響モデル分析部
１０３  スプリットマルチレート格子ＶＱ部
１０４、３００  コードブック指示値変換部
１０５  多重化部
１２１  サブベクトル特定部
１２２  位置情報符号化部
１２３  コードブック指示値分離部
１２４、２２２、３２２  コードブック指示値推定部
１２５、３２３  差分算出部
１２６  差分符号化部
２０１  多重分離部
２０２、４００  コードブック指示値逆変換部
２０３  スプリットマルチレート格子逆ＶＱ部
２０４  周波数−時間変換部
２２１  位置情報復号部
２２３  差分復号部
２２４  加算部
２２５、４２２  コードブック指示値生成部
３０１、４０１  判定部
３０２  第１変換処理部
３０３  第２変換処理部
３２１  固定コードブック指示値分離部
４０２  第１逆変換処理部
４０３  第２逆変換処理部
４２１  固定サブベクトル指定部
【図１】
0006027538.tif 000007
【図２】
0006027538.tif 000008
【図３】
0006027538.tif 000009
【図４】
0006027538.tif 000010
【図５】
0006027538.tif 000011
【図６】
0006027538.tif 000012
【図７】
0006027538.tif 000013
【図８】
0006027538.tif 000014
【図９】
0006027538.tif 000015
【図１０】
0006027538.tif 000016
【図１１】
0006027538.tif 000017
【図１２】
0006027538.tif 000018
【図１３】
0006027538.tif 000019
【図１４】
0006027538.tif 000020
【図１５】
0006027538.tif 000021
【図１６】
0006027538.tif 000022
【図１７】
0006027538.tif 000023
【図１８】
0006027538.tif 000024
【図１９】
0006027538.tif 000025
【図２０】
0006027538.tif 000026
【図２１】
0006027538.tif 000027
【図２２】
0006027538.tif 000028
【図２３】
0006027538.tif 000029
【図２４】
0006027538.tif 000030
【図２５】
0006027538.tif 000031