"複数のサンプルをまとめてベクトル量子化してベクトル量子化インデックスと前記複数のサンプルそれぞれの量子化値を求めるベクトル量子化ステップと、複数組の、前記サンプルそれぞれの位置に対応する予め定めた係数による組の中から、前記サンプルの値と、前記サンプルの量子化値と前記サンプルのそれぞれの位置に対応する係数とを乗算した値との間の誤差の全てのサンプル位置についての和を最小にする係数の組を表すインデックス情報を出力する係数組選択ステップと、を有する符号化方法。請求項1の符号化方法であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数に対応する値を第1軸にとり前記係数の値を第2軸にとる平面上の直線上にある係数により構成される。請求項1の符号化方法であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数の値を第1軸にとり前記係数の値を第2軸にとる平面上の直線又は特定の曲線に偏って分布する係数により構成される。請求項1から3の何れかの符号化方法であって、前記係数組選択ステップが出力するインデックス情報のビット数は、ベクトル量子化インデックスに対応する符号用に割り当てられたビット数から前記ベクトル量子化インデックスに対応する符号に実際に用いられたビット数を減算して得られる数以下である。入力されたベクトル量子化インデックスに対応する複数の値を求めて複数のサンプル位置に対応する復号値とするベクトル復号ステップと、入力されたインデックス情報が表す、前記複数のサンプル位置に対応する予め定めた係数による組を用い、前記各サンプル位置に対応する、前記復号値と前記係数とを乗算した値を出力する係数乗算ステップと、を有する復号方法。請求項5の復号方法であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数に対応する値を第1軸にとり前記係数の値を第2軸にとる平面上の直線上にある係数により構成される。請求項5の復号方法であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数の値を第1軸にとり前記係数の値を第2軸にとる平面上の直線又は特定の曲線に偏って分布する係数により構成される。複数のサンプルをまとめてベクトル量子化してベクトル量子化インデックスと前記複数のサンプルそれぞれの量子化値を求めるベクトル量子化部と、複数組の、前記サンプルそれぞれの位置に対応する予め定めた係数による組の中から、前記サンプルの値と、前記サンプルの量子化値と前記サンプルのそれぞれの位置に対応する係数とを乗算した値との間の誤差の全てのサンプル位置についての和を最小にする係数の組を表すインデックス情報を出力する係数組選択部と、を有する符号化装置。請求項8の符号化装置であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数に対応する値を第1軸にとり前記係数の値を第2軸にとる平面上の直線上にある係数により構成される。請求項8の符号化装置であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数の値を第1軸にとり前記係数の値を第2軸にとる平面上の直線又は特定の曲線に偏って分布する係数により構成される。請求項8から10の何れかの符号化装置であって、前記係数組選択ステップが出力するインデックス情報のビット数は、ベクトル量子化インデックスに対応する符号用に割り当てられたビット数から前記ベクトル量子化インデックスに対応する符号に実際に用いられたビット数を減算して得られる数以下である。入力されたベクトル量子化インデックスに対応する複数の値を求めて複数のサンプル位置に対応する復号値とするベクトル復号部と、入力されたインデックス情報が表す、前記複数のサンプル位置に対応する予め定めた係数による組を用い、前記各サンプル位置に対応する、前記復号値と前記係数とを乗算した値を出力する係数乗算部と、を有する復号装置。請求項12の復号装置であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数に対応する値を第1軸にとり前記係数の値を第2軸にとる平面上の直線上にある係数により構成される。請求項12の復号装置であって、前記係数による組は、前記係数が対応するサンプルの位置に対応する時間又は周波数の値を第1軸にとり前記係数の値を第2軸にとる平面上の直線又は特定の曲線に偏って分布する係数により構成される。請求項1の符号化方法の各ステップをコンピュータに実行させるためのプログラム。請求項5の復号方法の各ステップをコンピュータに実行させるためのプログラム。請求項1の符号化方法の各ステップをコンピュータに実行させるためのプログラムを格納したコンピュータ読み取り可能な記録媒体。請求項5の復号方法の各ステップをコンピュータに実行させるためのプログラムを格納したコンピュータ読み取り可能な記録媒体。しかしながら、AVQ法の振幅の量子化精度はSVQ法と比べて荒い。たとえ復号信号から周波数成分が欠損する頻度が低くても振幅の量子化精度が荒ければ、ミュージカルノイズ等が発生する。このような問題はAVQ法に限られるものではなく、量子化精度に基づいてミュージカルノイズ等が発生する場合に共通する。また、このような問題は入力信号が周波数領域の信号である場合だけではなく、入力信号が時間領域の信号である場合にも生じ得る。本発明では、量子化精度に基づいて発生するミュージカルノイズ等を低減する技術が提供される。符号化では、複数組の、サンプルそれぞれの位置に対応する予め定めた係数による組の中から、サンプルの値と、サンプルの量子化値とサンプルのそれぞれの位置に対応する係数とを乗算した値との間の誤差の全てのサンプル位置についての和を最小にする係数の組を表すインデックス情報を出力する。復号では、入力されたベクトル量子化インデックスに対応する複数の値を求めて複数のサンプル位置に対応する復号値とし、入力されたインデックス情報が表す、複数のサンプル位置に対応する予め定めた係数による組を用い、各サンプル位置に対応する、復号値と当該係数とを乗算した値を出力する。符号化においては、複数のサンプルそれぞれの量子化値に乗じられる複数個の係数による組を表すインデックス情報を出力するため、復号時の量子化誤差を低減させ、ミュージカルノイズ等を低減させることができる。復号においては、インデックス情報が表す複数個の係数を複数の復号値に乗じるため、量子化誤差を低減でき、ミュージカルノイズ等を低減することができる。以下、この発明の一実施形態を詳細に説明する。構成図1に例示するように、実施形態の符号化装置11は、正規化値計算部112、正規化値量子化部113、ベクトル量子化部115及び勾配計算部116係数組選択部に相当を例えば含む。図1に例示するように、実施形態の復号装置12は、正規化値復号部121、ベクトル復号部122及び勾配修正部124を例えば含む。必要に応じて、符号化装置11は、周波数領域変換部111を例えば含んでいてもよい。復号装置12は、時間領域変換部125、平滑化部126を例えば含んでいてもよい。符号化処理符号化装置11は図2に例示する符号化方法の各ステップを実行する。入力信号Xkは正規化値計算部112、ベクトル量子化部115及び勾配計算部116に入力される。この例の入力信号Xkは、音響信号などの時系列信号である時間領域信号xnを周波数領域に変換して得られる周波数領域信号である。周波数領域の入力信号Xkが符号化装置11に直接入力されてもよいし、周波数領域変換部111が入力された時間領域の入力信号xnを周波数領域に変換して周波数領域の入力信号Xkが生成されてもよい。周波数領域変換部111が周波数領域の入力信号Xkを生成する場合、周波数領域変換部111は、入力された時間領域の入力信号xnを、例えばMDCTModifiedDiscreteCosineTransformにより周波数領域の入力信号Xkに変換して出力する。nは時間領域での信号の番号離散時間番号であり、kは周波数領域での信号サンプルの番号離散周波数番号である。nの値が大きいほど遅い時間に対応する。kの値が大きいほど高い周波数に対応する。1フレームがL個のサンプルで構成されているとして、時間領域信号xnはフレームごとに周波数領域に変換され、L個の周波数成分を構成する周波数領域の入力信号Xkk=0,1,...,L-1が生成される。Lは所定の正の数、例えば64や80などの1よりも大きな整数である。なお、MDCTを用いる場合には、入力された時系列信号であるL個のサンプルからなるフレームごとに周波数領域へ変換され、この変換対象となるフレームが1/2フレームずつ、すなわち、L/2個のサンプルずつシフトされる。正規化値計算部112は、フレームごとに、入力された入力信号XkのL個のサンプルのうちの所定の個数C0のサンプルを代表する値である正規化値τX0-を計算するステップE1。τX0-は、τX0の上付きバーを意味する。ここで、τは、1フレーム中のL個のサンプルのうちの所定の個数C0のサンプルにより構成される各サブバンドに対して、一意に割り当てられた0以上の整数とする。C0は、L、又は、1及びL以外のLの公約数である。なお、C0をLとすることは、L個のサンプルごとに正規化値を求めることを意味する。C0を1及びL以外のLの公約数にすることは、L個のサンプルをサブバンドに分割して、各サブバンドを構成するC0個のサンプルごとに正規化値を求めることを意味する。例えば、L=64であり、8個の周波数成分でサブバンドを構成するとした場合には、8個のサブバンドが構成され、各サブバンドの正規化値が計算される。また、C0がLである場合τ=0であり、正規化値τX0-はL個のサンプルを代表する値である。すなわち、C0がLである場合、フレームごとに1つの正規化値τX0-が計算される。一方、C0が1及びL以外のLの公約数である場合、τは1フレーム中の各サブバンドに対応する整数τ=0,...,1であり、正規化値τX0-はτに対応するサブフレームに属するC0個のサンプルを代表する値である。すなわち、C0が1及びL以外のLの公約数である場合、フレームごとに個の正規化値τX0-τ=0,...,1が計算される。また、C0の値にかかわらずk=τC0,,τ1C0-1である。正規化値計算部112で計算されたτX0-は、正規化値量子化部113に送られる。[正規化値τX0-の具体例]正規化値τX0-は、C0個のサンプルを代表する値である。言い換えると、正規化値τX0-は、C0個のサンプルに対応する値である。正規化値τX0-の例は、以下のようなC0個のサンプルのパワー平均値に対する平方根である。正規化値τX0-の他の例は、以下のようなC0個のサンプルのパワー合計値に対する平方根をC0で除した値である。正規化値τX0-の他の例は、以下のようなC0個のサンプルの平均振幅値である。正規化値τX0-は上記の例に限定されない[正規化値τX0-の具体例]の説明終わり。正規化値量子化部113は、正規化値τX0-を量子化した量子化正規化値τX-及びその量子化正規化値τX-に対応する正規化値量子化インデックスを求めるステップE2。τX-は、τXの上付きバーを意味する。量子化正規化値τX-はベクトル量子化部115に送られ、正規化値量子化インデックスに対応する符号ビットストリームが復号装置12に送られる。ベクトル量子化部115は、フレームごとに、入力された入力信号XkのL個のサンプルのうちの複数のサンプルXkをまとめてベクトル量子化してベクトル量子化インデックスを生成する。ベクトル量子化インデックスは、量子化代表ベクトルを表わすインデックスである。この例のベクトル量子化部115は、量子化正規化値τX-を用いて複数のXkを正規化し、正規化された複数のサンプルXk’を求める。例えば、ベクトル量子化部115は、XkをτX-で除算する、もしくはXkにτX-の逆数を乗ずることでXk’を求める。ベクトル量子化部115は、例えば、図示していないベクトルコードブック記憶部に記憶された複数の量子化代表ベクトルの中から、複数のXk’を成分とするベクトルに最も近い量子化代表ベクトルを選択し、選択された量子化代表ベクトルを表わすベクトル量子化インデックスを出力することにより、ベクトル量子化を行う。ベクトル量子化部115は、例えばC0個のXk’をまとめてベクトル量子化を行う。ベクトル量子化部115は、例えばAVQ法例えば非特許文献2等参照等のベクトル量子化方法を用いてベクトル量子化を行うが、これ以外のベクトル量子化方法が採用されてもよい。また、C0がサブバンドのサンプル数である場合、人間の聴覚特性を考慮した優先順位が低いサブバンドについてベクトル量子化がなされなくてもよい。優先順位は、例えば、人間が知覚しやすい周波数に対応するサブバンドほど高い。例えば、量子化正規化値τX-が大きいサブバンドほど優先順位を高くする。また、ベクトル量子化によって得られる符号のビット数は、入力信号に応じて様々である。入力信号によっては、ベクトル量子化によって得られる符号ベクトル量子化インデックスなどのビット数が、ベクトル量子化用に割り当てられたビット数未満となり、ベクトル量子化用に割り当てられたビットの一部が余る場合がある。なお、ベクトル量子化用に割り当てられたビットとは、符号化装置11から復号装置12に送られる符号のうちで、ベクトル量子化によって得られる符号ベクトル量子化インデックスに対応する符号のために割り当てられたビットを意味する。ベクトル量子化用に割り当てられたビット数は、ベクトル量子化用に割り当てられたビットのビット数を表す。ベクトル量子化用に割り当てられたビット数は、フレームごとに定められたものであってもよいし、サブバンドごとに定められたものであってもよい。また、ベクトル量子化用に割り当てられたビット数は、入力信号に応じて変動するものであってもよいし、入力信号にかかわらず固定されたものであってもよい。ベクトル量子化部115は、ベクトル量子化用に割り当てられたビットのうち、実際のベクトル量子化で使用されなかったビットのビット数を未使用ビット数Uとして計算する。本形態の例では、1フレームごとL個のサンプルごとに未使用ビット数Uを計算する。例えば、ベクトル量子化部115は、処理対象のフレームでのベクトル量子化用に割り当てられたビット数から、実際にそのフレームに属するL個のサンプルのベクトル量子化によって得られたベクトル量子化インデックスのビット数の総数を減じた値を未使用ビット数Uとする。Uは0以上の整数である。さらに、ベクトル量子化部115は、ベクトル量子化インデックスをローカルデコードした値である複数の量子化値X^kを出力する。例えば、ベクトル量子化部115は、ベクトル量子化インデックスによって表される量子化代表ベクトルの各成分Xk’を、量子化正規化値τX-を用いて逆正規化した値をXkの量子化値X^kとして出力する。例えば、ベクトル量子化部115は、Xk’とτX-との乗算値を量子化値X^kとして出力する。この例の量子化値X^kは復号装置12で得られる復号値X^kと等しい。ベクトル量子化がなされていないサブバンドの量子化値X^kは0となる。なお、X^は、Xの上付きハットを意味する。ベクトル量子化部115は、ベクトル量子化インデックス、未使用ビット数U及び量子化値X^kを勾配計算部116に送るステップE3。勾配計算部116は、MMAX組の、C0個の勾配係数傾き補正ゲインにより構成される組を、例えば図示していない記憶部に保持する。MMAXは2以上の整数である。例えば、勾配計算部116は、C0個の勾配係数傾き補正ゲインγmkk=0C0-1を要素とする勾配係数ベクトルγm=[γm0γmC0-1]複数個の勾配係数による組をmm=0MMAX-1行目の行ベクトルとする式1に示す勾配行列γを保持する。勾配計算部116は、フレームごとに、入力された入力信号XkのL個のサンプルのうちのC0個のサンプルXkk=0C0-1に対応する値を要素とする第1ベクトルと、これらC0個のサンプルXkの量子化値X^kk=0C0-1のそれぞれを勾配係数ベクトルγmの各要素γmkで修正したC0個の修正値に対応する値を要素とする第2ベクトルとの誤差を最小化する例えば最小にする勾配係数ベクトルの行番号m’を求め、行番号m’を表すインデックス情報idxをベクトル量子化用に割り当てられたビットのうちで使用されなかったビットの領域未使用ビット領域と呼ぶに書き込む。言い換えると、勾配計算部116は、複数の勾配係数γmkを要素とする勾配係数ベクトルγmを行ベクトルとする勾配行列γの中から、複数のサンプルXkに対応する値を要素とする第1ベクトルと、複数の量子化値X^kを勾配係数ベクトルγmの各要素で修正した複数の修正値に対応する値を要素とする第2ベクトルとの間の誤差を最小化する勾配係数ベクトルの行番号m’を表すインデックス情報idxを出力し、ベクトル量子化インデックスに対応する符号ビットストリームの未使用ビット領域に配置する。すなわち、勾配計算部116は、例えば、複数組の、サンプルXkそれぞれの位置に対応する予め定めた係数γmkによる組の中から、サンプルXkの値と、サンプルの量子化値X^kとサンプルのそれぞれの位置に対応する係数γmkとを乗算した値との間の誤差の全てのサンプル位置についての和を最小にする係数γmの組を表すインデックス情報idxを出力する。なお、本形態のサンプルXkそれぞれの位置は、各離散周波数番号kに対応する周波数軸上の位置であるステップE4。これにより、符号化装置11は、未使用ビット領域を有効利用して振幅の量子化誤差を修正するための情報を復号装置12に送ることができ、量子化精度に基づいて発生するミュージカルノイズ等を低減させることができる。勾配係数ベクトルγmの要素であるC0個の勾配係数γm0,γmC0-1は互いに相関を持つ。言い換えると、各勾配係数ベクトルγmは互いに相関を持つ複数の勾配係数γm0,γmC0-1を要素とするベクトルである。X0,XC0-1は、kを第1軸にとり、Xkを第2軸にとるkXk平面上の直線や曲線に偏って分布する場合が多い。このようなX0,XC0-1の特性を考慮した勾配係数γm0,γmC0-1を要素とする勾配係数ベクトルγmを用いることで、量子化誤差を高い精度で修正できる。例えば、同一の行番号mに対応する各勾配係数γm0,γmC0-1が、k勾配係数γmkが乗じられる量子化値X^kに対応する周波数に対応する値、言い換えると、勾配係数γmkに対応する周波数に対応する値を第1軸にとりγmk勾配係数の値を第2軸にとるkγmk平面上の直線又は特定の曲線に偏って分布するものとする。言い換えると、例えば列番号に対応する値kを第1軸にとり当該列番号での勾配係数γmkを第2軸にとるkγmk平面上の直線又は特定の曲線に偏って分布する複数の勾配係数γm0,γmC0-1を要素とするベクトルを勾配係数ベクトルγmとする。より具体的には、例えば同一の行番号mに対応する各勾配係数γm0,γmC0-1が、kγmk平面上の直線上又は特定の曲線上に位置するものとする。言い換えると、kγmk平面上の直線又は特定の曲線上の各勾配係数γm0,γmC0-1を要素とするベクトルを勾配係数ベクトルγmとした。kγmk平面上の直線上又は特定の曲線は、例えば行番号mごとに異なる。以下に勾配行列γを例示する。この例は、C0=8、MMAX=3の場合のものである。この例では、各勾配係数γm0,γm7が行番号mm=012ごとに定められた直線上に配置される。第1ベクトルの例は、C0個のXkを要素とするベクトル、C0個のXkの大きさ|Xk|を要素とするベクトル、C0個のXk又は|Xk|に定数や変数を乗じた値を要素とするベクトルなどである。第2ベクトルの例は、C0個の修正値を要素とするベクトル、C0個の修正値の大きさを要素とするベクトル、C0個の修正値又は大きさに定数や変数を乗じた値を要素とするベクトルなどである。修正値の例は、X^kとγmkとの積、X^kの大きさ|X^k|とγmkとの積、X^kとγmkとの積の大きさ、X^kとγmkとの積の大きさを表す値、その他X^kとγmkとの積に対応する値などである。第1ベクトルと第2ベクトルとの間の誤差の例は、第1ベクトルと第2ベクトルとの間の距離である。距離の定義に限定はなく、マンハッタン距離、ユークリッド距離、それらの変形等どのようなものでもよい。第1ベクトルと第2ベクトルとの誤差を最小化する勾配係数ベクトルの例は、第1ベクトルと第2ベクトルとの誤差を最小にする勾配係数ベクトルや、或る探索条件や探索範囲内で第1ベクトルと第2ベクトルとの誤差を最小にする勾配係数ベクトルなどである。未使用ビット領域は、例えば、定められた未使用ビット領域の基準位置例えば先頭アドレスと入力された未使用ビット数Uとによって特定できる。ただし、未使用ビット領域に書き込むことが可能なインデックス情報idxのビット数の上限は未使用ビット数Uである。よって、必ずしもあらゆる行番号に対応するインデックス情報idxを未使用ビット領域に書き込めるとは限らない。そこで勾配計算部116は、例えば、未使用ビット領域に書き込むことが可能なインデックス情報idxで識別可能な行番号の範囲を探索範囲とし、行番号m’を選択する。言い換えると、勾配計算部116は、未使用ビット領域に書き込むことが可能なインデックス情報idxが表す行番号m’のみを選択する。すなわち、勾配計算部116は、ベクトル量子化インデックスに対応する符号用に割り当てられたビット数うち実際にベクトル量子化インデックスに対応する符号に用いられなかったビット数で表現可能なインデックス情報idxで識別可能な行番号m’のみを選択する。例えば、勾配計算部116は、未使用ビット領域に書き込むことが可能なインデックス情報idxで識別可能なmMAX種類の行番号m=0,mMAX-1から、以下のような行番号m’を特定し、行番号m’に対応するインデックス情報idxを未使用ビット領域に書き込む。なお、∥∥はのノルムであり、argminm∥∥は∥∥を最小にするmをm’とすることを意味し、argminmはargminの下付きのmを意味し、χ=[X0,,XC0-1]、χ^=[X^0,,X^C0-1]であり、Amは以下のような行番号mに対応する勾配係数ベクトルγm=[γm0,γmC0-1]を対角成分に持つ対角行列を意味する。上述のようなインデックス情報idxのビット数は、ベクトル量子化インデックスに対応する符号用に割り当てられたビット数からベクトル量子化インデックスに対応する符号に実際に用いられたビット数を減算して得られる数以下である。以上より、未使用ビット領域のみを用いてインデックス情報idxを伝送できる。[ステップE4の具体例]この例の勾配計算部116は、図3に記載された各処理を行い、選択した勾配係数ベクトルの行番号を表すインデックス情報idxを、未使用ビット領域に書き込む。C0がLである場合には、フレームごとに図3のステップE4の処理が実行され、C0が1及びL以外のLの公約数である場合には、1フレーム中の各サブバンドについて図3のステップE4の処理が繰り返し実行される。勾配計算部116は、入力された未使用ビット数Uと0とを比較しステップE40、U>0でなければ、以下に示すように、勾配計算部116は、入力された複数の量子化値X^0,,X^C0-1を更新することなく、ステップE4の処理を終えるステップE412。U>0であれば、勾配計算部116は、m=0,idx=0とすることでm及びidxの値を初期化しステップE41、ステップE42に進む。ステップE42では、勾配計算部116は、未使用ビット数Uを用い、未使用ビット領域に書き込むことが可能なインデックス情報idxで識別可能な行番号の範囲を探索範囲として定め、この探索範囲行番号の範囲を特定するための探索範囲特定値mMAXを特定する。言い換えると、勾配計算部116は、未使用ビット領域に書き込むことが可能なインデックス情報idxで識別可能な行番号の種類を特定するmMAXを求めるステップE42。通常、未使用ビット数Uのビットでは2U種類の行番号を識別可能である。そのため、2U種類の行番号の範囲を探索範囲としてもよい。しかしながら、この具体例では勾配係数ベクトルγm’を用いた補正を行わないことを示す値をインデックス情報idxが取り得る値の1つに割り当て、残りの2U-1個の値で識別可能な2U-1種類の行番号を探索範囲とする。探索範囲行番号の範囲とmMAXとの関係は予め定められる必要がある。この具体例のmMAXは1以上2U-1以下なおかつMMAX以下の整数であり、探索範囲が行番号0,mMAX-1である。勾配計算部116は、例えば以下の式によってmMAXを求める。mMAX=max[min{2U-1,MMAX},1]勾配計算部116は、以下の式に表される計算を行うステップE43。勾配計算部116は、mとmMAXとを比較ステップE44し、m<mMAXであれば、勾配計算部116は、以下のようにeを計算しステップE45、eMINとeとを比較しステップE46、eMIN>eでなければ勾配計算部116は、mを1だけインクリメントしステップE48、ステップE44に進む。eMIN>eであれば勾配計算部116は、idx=m1及びeMIN=eとすることでidx及びeMINを更新しステップE47、mを1だけインクリメントしてステップE48、ステップE44に進む。一方、ステップE44でm<mMAXでないと判定された場合、勾配計算部116は、未使用ビット領域にidxを書き込みステップE49。この例の勾配計算部116は、復号装置12が、mMAXに基づいて、必要なidxが未使用ビット領域のどこに配置されているのかを特定できるようにしておく。例えば、未使用ビット領域でのidxの格納位置がmMAXによって定まる。次に勾配計算部116は、idx>0であるか否かidx=0であるか否かを判定するステップE410。idx>0であればidx=0でなければ、勾配計算部116は、ローカルデコード値である複数の量子化値X^bC0,...,X^b1C0-1を、以下のように更新しステップE411、ステップE4の処理を終える。idx>0でなければidx=0であれば、以下に示すように、勾配計算部116は、ローカルデコード値である複数の量子化値X^bC0,...,X^b1C0-1を更新することなくステップE412、ステップE4の処理を終える。ただし、C0がLである場合b=0である。C0が1及びL以外のLの公約数である場合、bは1フレーム中の各サブバンドに対応する整数0,...,1の何れかであり、例えば最も低い周波数のサブバンドに対応する整数0である[ステップE4の具体例]の説明終わり。ベクトル量子化インデックスと未使用ビット領域に書き込まれたインデックス情報idxとを含む修正後ベクトル量子化インデックスに対応する符号ビットストリームは復号装置12に送られる。復号処理復号装置12は図4に例示する復号方法の各ステップを実行する。正規化値復号部121は、復号装置12に入力された正規化値量子化インデックスに対応する復号正規化値τX-を求めるステップD1。復号正規化値τX-は、ベクトル復号部122に送られる。図示していないコードブック記憶部に複数の正規化値量子化インデックスのそれぞれに対応する正規化値が記憶されているとする。正規化値復号部121は、入力された正規化量子化インデックスをキーとしてそのコードブック記憶部を参照して、その正規化量子化インデックスに対応する正規化値を取得して、復号正規化値τX-とする。ベクトル復号部122は、復号装置12に入力された修正後ベクトル量子化インデックスが含むベクトル量子化インデックスに対応する複数の値を求めて複数の復号値X^kとする。また、ベクトル復号部122は、ベクトル量子化インデックスを用いて未使用ビット数Uを計算するステップD2。本形態では、図示していないベクトルコードブック記憶部に複数のベクトル量子化インデックスのそれぞれに対応する量子化代表ベクトルが記憶されているとする。ベクトル復号部122は、入力されたベクトル量子化インデックスに対応する量子化代表ベクトルをキーとしてそのベクトルコードブック記憶部を参照して、そのベクトル量子化インデックスに対応する量子化代表ベクトルを取得する。ベクトル復号部122は、量子化代表ベクトルの成分Xk’を量子化正規化値τX-で逆正規化した復号値X^kを出力する。例えば、ベクトル復号部122は、Xk’とτX-との乗算値を復号値X^kとして出力する。また、ベクトル復号部122は、ベクトル量子化用に割り当てられたビットのうち、実際のベクトル量子化で使用されなかったビットのビット数を未使用ビット数Uとして計算する。本形態の例では、1フレームごとL個のサンプルごとに未使用ビット数Uを計算する。例えば、ベクトル復号部122は、処理対象のフレームでのベクトル量子化用に割り当てられたビット数から、そのフレームに対応するベクトル量子化インデックスのビット数の総数を減じた値を未使用ビット数Uとする。復号値E^k及び未使用ビット数Uは、勾配修正部124に送られる。勾配修正部124は、符号化装置11で使用したのと同じ勾配行列γ式1参照を、例えば図示していない記憶部に保持する。勾配修正部124は、復号装置12に入力された修正後ベクトル量子化インデックスが含む未使用ビット領域からインデックス情報idxを読み出し、idxが表す行番号m’の勾配係数ベクトルγm’の各要素γm’kk=0C0-1を用い、C0個の復号値X^kk=0C0-1を修正する。言い換えると、勾配修正部124は、複数の勾配係数γmkを要素とする勾配係数ベクトルγmを行ベクトルとする勾配行列γ中のインデックス情報idxが表す行番号m’の勾配係数ベクトルγm’の各要素γm’kを用い、複数の復号値X^kを修正するステップD3。例えば、勾配修正部124は、インデックス情報idxが表す行番号m’の勾配係数ベクトルγm’の要素γm’kと復号値X^kとの乗算値を復号値X^kの修正値X^UDkとする。勾配修正部124は修正値X^UDkを出力する。すなわち、勾配修正部124は、入力されたインデックス情報idxが表す、複数のサンプル位置に対応する予め定めた係数γm’kによる組を用い、各サンプル位置に対応する復号値X^kと係数γm’kとを乗算した値を出力する。[ステップD3の具体例]この例の勾配修正部124は、図5に記載された各処理を行い、復号値X^kを修正する。勾配修正部124は、入力された未使用ビット数Uと0とを比較しステップD30、U>0でなければ、以下に示すように、勾配計算部116は、入力された複数の量子化値X^0,,X^C0-1を更新することなくステップD36、ステップD3の処理を終える。U>0であれば、勾配修正部124は、前述のステップE42と同じ方法によってmMAXを設定するステップE32。勾配修正部124は、例えば以下の式によってmMAXを求める。mMAX=max[min{2U-1,MMAX},1]勾配修正部124は、mMAXに基づき、修正後ベクトル量子化インデックスの未使用ビット領域からインデックス情報idxを読み込むステップD33。例えば、勾配修正部124は、mMAXに基づいてインデックス情報idxの格納位置を特定し、インデックス情報idxを読み込む。勾配修正部124は、idx>0であるか否かidx=0であるか否かを判定するステップD34。idx>0であればidx=0でなければ、勾配修正部124は、複数の復号値X^bC0,...,X^b1C0-1を、以下のように更新しステップD35、ステップD3の処理を終える。idx>0でなければidx=0であれば、以下に示すように、勾配修正部124は、複数の復号値X^bC0,...,X^b1C0-1を更新することなくステップD36、ステップD3の処理を終える。以上で[ステップD3の具体例]の説明を終える。時間領域での復号信号が必要な場合、勾配修正部124から出力された修正値X^UDkは時間領域変換部125に入力され、時間領域変換部125はX^UDkを例えば逆フーリエ変換により時間領域信号znに変換して出力する。本形態の特徴このように本形態では、復号装置12が符号化装置11で選択された勾配係数ベクトルを用いて複数の復号値X^kを修正するため、量子化誤差に基づくミュージカルノイズ等を軽減できる。また本形態では、互いに相関を持つ勾配係数γm0,γmC0-1を要素とするベクトルを勾配係数ベクトルγmとした。例えば、kγmk平面上の直線又は特定の曲線に偏って分布する複数の勾配係数γm0,γmC0-1を要素とするベクトルを勾配係数ベクトルγmとした。音声信号や音響信号などの入力信号のエンベロープは直線や曲線をなす場合が多く、このような入力信号の特性を反映させた勾配係数ベクトルγm用いることにより、インデックス情報idxの情報量を抑制しつつ量子化誤差を高精度に修正できる。図6の例の場合、k=0,...,63での入力信号の大きさ|Xk|はkが大きいほど小さい。そのため、kγmk平面上の傾きが負となる直線に偏って分布する勾配係数γm0,γm63を要素とする勾配係数ベクトルγmを用いて|X0|,|X63|を修正することで、量子化値の大きさ|X^0|,|X^63|との誤差を小さくできる。このようにサブバンドごとに入力信号の特性に適した勾配係数ベクトルγmを用いることで量子化誤差を効率的に削減することができる。また、符号化装置11で選択された勾配係数ベクトルγm’を特定するためのインデックス情報idxは未使用ビット領域を有効利用して伝送されるため、インデックス情報idxを伝送するための付加領域が不要である。変形例:なお、本発明は上述の実施の形態に限定されるものではない。例えば、復号装置12が平滑化部126を含む場合、平滑化部126は、ステップD3図4で得られた修正値X^UDkを入力とし、修正値X^UDkよりも過去の修正値X^UDk’が0でない場合、過去の修正値X^UDk’と修正値X^UDkとを重み付き加算した値を平滑化値X^POSTkとして出力する。X^UDk’が0の場合、平滑化部126は、修正値の重み付き加算、すなわち修正値の平滑化を行わず、X^UDkをX^POSTkとして出力する図4/ステップD4。過去の修正値X^UDk’の例は、修正値X^UDkに対応するフレームの1フレーム過去のステップD3で得られた修正値や、修正値X^UDkに対応するフレームの1フレーム過去のステップD4で得られた平滑化値などである。X^POSTkは次式のように表わされる。αとβは調整係数であり、求める性能及び仕様に応じて適宜決定される。例えば、α=0.85、β=0.15である。α、βは、求められる性能及び仕様に応じて適宜変更されてもよい。φはの正負符号を表す。これにより、X^UDkの振幅特性の時間軸方向の不連続性によって生じるミュージカルノイズ等を軽減することができる。時間領域での復号信号が必要な場合、平滑化部126から出力されたX^POSTkは時間領域変換部125に入力され、時間領域変換部125はX^POSTkを例えば逆フーリエ変換により時間領域信号znに変換して出力する。また、入力信号Xkは、周波数領域信号である必要はなく、時間領域信号等の任意の信号であってもよい。すなわち、この発明は、周波数領域信号以外の任意の信号に対する符号化、復号に用いることができる。この場合、同一の行番号mに対応する各勾配係数γm0,γmC0-1が、例えばk勾配係数γmkが乗じられる量子化値X^kに対応する時間に対応する値、言い換えると、勾配係数γmkに対応する時間に対応する値を第1軸にとりγmk勾配係数の値を第2軸にとるkγmk平面上の直線又は特定の曲線に偏って分布するものとする。より具体的には、例えば同一の行番号mに対応する各勾配係数γm0,γmC0-1が、kγmk平面上の直線上又は特定の曲線上に位置するものとする。なお、この変形例ではkが離散時間に対応する離散時間番号であり、サンプルXkそれぞれの位置は、離散時間番号kに対応する時間軸上の位置である。kが離散時間番号である場合、kの値が大きいほど遅い時間に対応する。また、フレームごとに入力信号Xkに対する正規化値FGAINが定められ、ベクトル量子化部115が入力信号の各サンプルの値Xkの代わりにXkを正規化値FGAINで正規化した値を用い、量子化正規化値τX-の代わりにτX-を正規化値FGAINで正規化した値を用い、ステップE3の処理が実行されてもよい。例えば、XkがXk/FGAINに置換され、τX-がτX-/FGAINに置換されて、ステップE3の処理が実行されてもよい。また、この場合、正規化値計算部112が存在せず、正規化値量子化部113に量子化正規化値τX-の代わりにXkを正規化値FGAINで正規化した値が入力されてもよい。この場合、ベクトル量子化部115が、量子化正規化値τX-の代わりにXkを正規化値FGAINで正規化した値の量子化値を用いてステップE3の処理を行ってもよい。正規化値量子化インデックスが正規化値FGAINで正規化した値の量子化値に対応するものであってもよい。また、上記実施形態では、符号化装置11の勾配計算部116が、idx>0であるか否かを判定し、idx>0であれば複数の量子化値X^bC0,...,X^b1C0-1を更新し、idx>0でなければ更新しないこととした図3/ステップE410~E412。また、復号装置12の勾配修正部124が、idx>0であるか否かを判定し、idx>0であれば複数の量子化値X^bC0,...,X^b1C0-1を更新し、idx>0でなければ更新しないこととした図5/ステップD34~D36。これらの変形として、式1に示した勾配行列γに1の要素のみからなる行番号m=-1の行ベクトル勾配係数ベクトルγ-1=[γ-10γ-1C0-1]=[11]を付加し、idx>0であるか否かにかかわらず勾配計算部116や勾配修正部124が、以下のような演算を行ってもよい。また、例示した行番号mやインデックス情報idxの具体的な数値は本発明を限定するものではなく、例示したmやidxの番号が繰り上がっても良いし、繰り下がってもよいし、一部の番号が使用されなくてもよい。また、上記の実施形態ではインデックス情報idxを未使用ビット数Uの未使用ビット領域に格納する例を示したが、インデックス情報idxが未使用ビット領域に格納されない形態であってもよい。また、上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。ハードウェア、プログラム及び記録媒体:上述した符号化装置11や復号装置12は、例えば、CPUやRAMなどからなる公知又は専用のコンピュータと上述の処理内容が記述された特別なプログラムとから構成される。この場合、特別なプログラムはCPUに読み込まれ、CPUが特別なプログラムを実行することによって各機能が実現される。また、特別なプログラムは、単一のプログラム列によって構成されてもよいし、他のプログラムやライブラリを読み出して目的の機能を達成するものであってもよい。このようなプログラムは、コンピュータで読み取り可能記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。コンピュータで読み取り可能な記録媒体の例は非一時的なnontransitory記録媒体である。プログラムの流通は、例えば、そのプログラムを記録したDVD、CD-ROM等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。また、プログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、プログラムを流通させてもよい。このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、符号化装置11や復号装置12の処理部の少なくとも一部が特別な集積回路であってもよい。"