(57)【要約】  線形予測符号化装置は、時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る線形予測分析部２２１と、符号帳記憶部２２２に記憶された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補と、線形予測分析部２２１が得た線形予測係数に変換可能な係数と、のηの値を適合させる適合部２２Ａと、ηの値が適合された線形予測係数に変換可能な係数の複数個の候補と線形予測係数に変換可能な係数とを用いて、線形予測分析部２２１が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得る符号化部２２４と、を備えている。
2016167215.tif 000002
【特許請求の範囲】
【請求項１】
パラメータηを正の数として、時系列信号に対応するパラメータηを、その時系列信号に対応する周波数領域サンプル列の絶対値のη乗をパワースペクトルと見做すことにより推定されたスペクトル包絡で上記周波数領域サンプル列を除算した系列である白色化スペクトル系列のヒストグラムを近似する一般化ガウス分布の形状パラメータとし、η1はパラメータηの所定の値であるとして、
時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る線形予測分析部と、
Ｎ種類（Ｎは１以上の整数）のパラメータηのそれぞれに対応するＮ個の符号帳が記憶され、各符号帳にはそれぞれのパラメータηに対応する線形予測係数に変換可能な係数の候補が複数個格納された符号帳記憶部と、
上記符号帳記憶部に記憶された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補と、上記線形予測分析部が得た線形予測係数に変換可能な係数と、のηの値を適合させる適合部と、
上記ηの値が適合された線形予測係数に変換可能な係数の複数個の候補と線形予測係数に変換可能な係数とを用いて、上記線形予測分析部が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得る符号化部と、
を含む線形予測符号化装置。
【請求項２】
請求項１の線形予測符号化装置であって、
上記適合部は、上記符号帳記憶部に記憶された線形予測係数に変換可能な係数の候補に対して、η1に応じた第一線形変換を行い、第一線形変換後の線形予測係数に変換可能な係数の複数個の候補を得る線形変換部を含み、
上記符号化部は、上記線形予測分析部が得た線形予測係数に変換可能な係数と、上記適合部が得た上記第一線形変換後の線形予測係数に変換可能な係数の複数個の候補と、を用いて、上記線形予測分析部が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得る、
線形予測符号化装置。
【請求項３】
請求項１の線形予測符号化装置であって、
上記適合部は、上記線形予測分析部が得た線形予測係数に変換可能な係数に対して、η1に応じた第二線形変換を行い、第二線形変換後の線形予測係数に変換可能な係数を得る線形変換部を含み、
上記符号化部は、上記適合部が得た上記第二線形変換後の線形予測係数に変換可能な係数と、上記符号帳に格納された線形予測係数に変換可能な係数の複数個の候補と、を用いて、上記線形予測分析部が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得る、
線形予測符号化装置。
【請求項４】
請求項１の線形予測符号化装置であって、
η2,η3はパラメータηの所定の値であるとして、
上記符号帳記憶部には、η2に対応する符号帳が記憶されており、
上記適合部は、上記符号帳記憶部に記憶された線形予測係数に変換可能な係数の複数個の候補に対して、η3に応じた第一線形変換を行い、第一線形変換後の線形予測係数に変換可能な係数の複数個の候補を得、上記線形予測分析部が得た線形予測係数に変換可能な係数に対して、η3に応じた第二線形変換を行い、第二線形変換後の線形予測係数に変換可能な係数を得る、線形変換部であり、
上記符号化部は、上記適合部が得た上記第二線形変換後の線形予測係数に変換可能な係数と、上記適合部が得た上記第一線形変換後の線形予測係数に変換可能な係数の複数個の候補と、を用いて、上記線形予測分析部が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得る、
線形予測符号化装置。
【請求項５】
請求項１の線形予測符号化装置であって、
η2はパラメータηの所定の値であるとして、
上記符号帳記憶部には、複数の符号帳が記憶されており、
上記適合部は、上記符号帳記憶部に記憶された複数の符号帳の中から上記η2に応じて符号帳を選択する符号帳選択部と、上記線形予測分析部で得られた線形予測係数に変換可能な係数に対する、η2に応じた第二線形変換を行う線形変換部とであり、
上記符号化部は、上記第二線形変換後の線形予測係数に変換可能な係数について、上記選択された符号帳を用いて符号化して線形予測係数符号を得る、
線形予測符号化装置。
【請求項６】
請求項１の線形予測符号化装置であって、
η2はパラメータηの所定の値であるとして、
上記符号帳記憶部には、複数の符号帳が記憶されており、
上記適合部は、上記符号帳記憶部に記憶された複数の符号帳の中から上記η2に応じて符号帳を選択する符号帳選択部と、上記選択された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補に対する、η1に応じた第一線形変換を行う線形変換部とであり、
上記符号化部は、上記線形予測分析部で得られた線形予測係数に変換可能な係数について、上記第一線形変換後の線形予測係数に変換可能な係数の候補を用いて符号化して線形予測係数符号を得る、
線形予測符号化装置。
【請求項７】
請求項１の線形予測符号化装置であって、
η2,η3はパラメータηの所定の値であるとして、
上記符号帳記憶部には、複数の符号帳が記憶されており、
上記適合部は、上記符号帳記憶部に記憶された複数の符号帳の中から上記η3に応じて符号帳を選択する符号帳選択部と、上記選択された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補に対する、η2に応じた第一線形変換を行うと共に、上記線形予測分析部で得られた線形予測係数に変換可能な係数に対する、η2に応じた第二線形変換を行う線形変換部とであり、
上記符号化部は、上記第二線形変換後の線形予測係数に変換可能な係数について、上記第一線形変換後の線形予測係数に変換可能な係数の候補を用いて符号化して線形予測係数符号を得る、
線形予測符号化装置。
【請求項８】
請求項２の線形予測符号化装置において、
上記線形変換部は、上記η1が小さいほど上記第一線形変換後の線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列が平坦になるように上記第一線形変換を行う、
線形予測符号化装置。
【請求項９】
請求項２から８の何れかの線形予測符号化装置において、
pを線形予測係数に変換可能な係数の次数とし、上記線形予測係数に変換可能な係数又は上記線形予測係数に変換可能な係数の候補を^ω[k][k=1,2,…,p]とし、上記第一線形変換後及び上記第二線形変換後の線形予測係数に変換可能な係数又は上記線形予測係数に変換可能な係数の候補を~ω[k][k=1,2,…,p]とし、x1,x2,…xp,y1,y2,…yp-1,z2,z3,…zpを所定の非負の数とし、y1,y2,…yp-1,z2,z3,…zpの少なくとも１つは所定の正の数であるとし、Kをx1,x2,…xp,y1,y2,…yp-1,z2,z3,…zp以外の要素が０である行列として、
上記線形変換部は、下記式により上記第一線形変換と上記第二線形変換との少なくとも一方を行う、
【数１７】
2016167215.tif 000019
線形予測符号化装置。
【請求項１０】
請求項２の線形予測符号化装置において、
上記線形変換部は、上記η1が小さいほど上記第一線形変換後の線形予測係数に変換可能な係数の候補の次数が小さくなるように上記第一線形変換を行う、
線形予測符号化装置。
【請求項１１】
請求項１，２，３，４の何れかの線形予測符号化装置であって、
上記符号帳記憶部には、複数の符号帳が記憶されており、
上記適合部は、上記符号帳記憶部に記憶された複数の符号帳の中から上記η1に応じて符号帳を選択する符号帳選択部を含み、
上記符号化部は、上記線形予測分析部が得た線形予測係数に変換可能な係数と、上記適合部が得た線形予測係数に変換可能な係数の複数個の候補と、を用いて、上記線形予測分析部が得た上記線形予測係数に変換可能な係数に対応する線形予測係数符号を得る、
線形予測符号化装置
【請求項１２】
請求項１１の線形予測符号化装置において、
上記符号帳記憶部には、線形予測係数に変換可能な係数の候補数が異なる複数の符号帳が記憶されており、
上記符号帳選択部は、上記η1が大きいほど、上記符号帳記憶部に記憶された複数の符号帳の中から、線形予測係数に変換可能な係数の候補数が多い符号帳を選択する、
線形予測符号化装置。
【請求項１３】
請求項１１又は１２の線形予測符号化装置において、
上記符号帳記憶部には、符号帳に記憶された線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列の平坦度合いが異なる複数の符号帳が記憶されており、
上記符号帳選択部は、上記η1が小さいほど、上記符号帳記憶部に記憶された複数の符号帳の中から、符号帳に記憶された線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列がより平坦である符号帳を選択する、
線形予測符号化装置。
【請求項１４】
請求項１１又は１２の線形予測符号化装置において、
上記符号帳記憶部には、線形予測係数に変換可能な係数の候補間の間隔が異なる複数の符号帳が記憶されており、
上記符号帳選択部は、上記ηが小さいほど、上記符号帳記憶部に記憶された複数の符号帳の中から、線形予測係数に変換可能な係数の候補間の間隔が狭い符号帳を選択する、
線形予測符号化装置。
【請求項１５】
パラメータηを正の数として、時系列信号に対応するパラメータηを、その時系列信号に対応する周波数領域サンプル列の絶対値のη乗をパワースペクトルと見做すことにより推定されたスペクトル包絡で上記周波数領域サンプル列を除算した系列である白色化スペクトル系列のヒストグラムを近似する一般化ガウス分布の形状パラメータとし、η1はパラメータηの所定の値であるとして、
時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る線形予測分析部と、
符号帳が記憶された符号帳記憶部と、
入力されたη1に基づいて、上記符号帳記憶部に記憶された符号帳と上記線形予測係数に変換可能な係数との少なくとも一方を適合させる適合部と、
上記符号帳又は上記適合された符号帳を用いて、上記線形予測係数に変換可能な係数又は上記適合された線形予測係数に変換可能な係数を符号化する符号化部と、
を含む線形予測符号化装置。
【請求項１６】
符号帳が記憶された符号帳記憶部と、
η1を正の数として、入力されたη1に基づいて、上記符号帳記憶部に記憶された符号帳と、上記符号帳に格納された複数個の線形予測係数に変換可能な係数の候補のうち、入力された線形予測係数符号に対応する線形予測係数に変換可能な係数の候補との少なくとも一方を適合させる適合部を含み、
上記線形予測係数に変換可能な係数は、上記線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列を得るために用いられる、
線形予測復号装置。
【請求項１７】
請求項１６の線形予測復号装置であって、
上記符号帳に格納された複数個の線形予測係数に変換可能な係数の候補のうち、入力された線形予測係数符号に対応する線形予測係数に変換可能な係数の候補を線形予測係数に変換可能な係数として得る復号部を更に含み、
上記適合部は、上記復号部で得られた線形予測係数に変換可能な係数に対して、所定の正の数であるη1に応じた線形変換をして線形予測係数に変換可能な係数を得る線形変換部である、
線形予測復号装置。
【請求項１８】
請求項１６の線形予測復号装置であって、
上記符号帳には、複数の符号帳が記憶されており、
η2を正の数として、上記適合部は、上記符号帳記憶部に記憶された複数の符号帳の中から上記η２に応じて符号帳を選択する符号帳選択部と、復号部で得られた線形予測係数に変換可能な係数に対して、所定の正の数であるη1に応じた線形変換をして線形予測係数に変換可能な係数を得る線形変換部とであり、
上記選択された符号帳に格納された複数個の線形予測係数に変換可能な係数の候補のうち、入力された線形予測係数符号に対応する線形予測係数に変換可能な係数の候補を線形予測係数に変換可能な係数として得る上記復号部を更に含む、
線形予測復号装置。
【請求項１９】
請求項１７の線形予測復号装置において、
上記線形変換部は、上記η1が小さいほど上記線形変換部で得られた線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列が平坦になるように上記線形変換を行う、
線形予測復号装置。
【請求項２０】
請求項１７から１９の何れかの線形予測復号装置において、
pを線形予測係数に変換可能な係数の次数とし、上記復号部で得られた線形予測係数に変換可能な係数を^ω[k][k=1,2,…,p]とし、上記線形変換後の線形予測係数に変換可能な係数を~ω[k][k=1,2,…,p]とし、x1,x2,…xp,y1,y2,…yp-1,z2,z3,…zpを所定の非負の数とし、y1,y2,…yp-1,z2,z3,…zpの少なくとも１つは所定の正の数であるとし、Kをx1,x2,…xp,y1,y2,…yp-1,z2,z3,…zp以外の要素が０である行列として、
上記線形変換部は、下記式により線形変換を行う、
【数１８】
2016167215.tif 000020
線形予測復号装置。
【請求項２１】
請求項１７の線形予測復号装置において、
上記線形変換部は、上記η1が小さいほど上記線形変換後の線形予測係数に変換可能な係数の次数が小さくなるように上記線形変換を行う、
線形予測復号装置。
【請求項２２】
請求項１６の線形予測復号装置であって、
上記符号帳には、複数の符号帳が記憶されており、
上記適合部は、上記符号帳記憶部に記憶された複数の符号帳の中から上記η1に応じて符号帳を選択する符号帳選択部であり、上記選択された符号帳を用いて、入力された線形予測係数符号を復号して線形予測係数に変換可能な係数を得る復号部を更に含む
線形予測復号装置。
【請求項２３】
請求項２２の線形予測復号装置において、
上記符号帳記憶部には、線形予測係数に変換可能な係数の候補数が異なる複数の符号帳が記憶されており、
上記符号帳選択部は、上記η1が大きいほど、上記符号帳記憶部に記憶された複数の符号帳の中から、線形予測係数に変換可能な係数の候補数が多い符号帳を選択する、
線形予測復号装置。
【請求項２４】
請求項２２又は２３の線形予測復号装置において、
上記符号帳記憶部には、符号帳に記憶された線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列の平坦度合いが異なる複数の符号帳が記憶されており、
上記符号帳選択部は、上記η1が小さいほど、上記符号帳記憶部に記憶された複数の符号帳の中から、符号帳に記憶された線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列が平坦である符号帳を選択する、
線形予測復号装置。
【請求項２５】
請求項２２又は２３の線形予測復号装置において、
上記符号帳記憶部には、線形予測係数に変換可能な係数の候補間の間隔が異なる複数の符号帳が記憶されており、
上記符号帳選択部は、上記η1が小さいほど、上記符号帳記憶部に記憶された複数の符号帳の中から、線形予測係数に変換可能な係数の候補間の間隔が狭い符号帳を選択する、
線形予測復号装置。
【請求項２６】
パラメータηを正の数として、時系列信号に対応するパラメータηを、その時系列信号に対応する周波数領域サンプル列の絶対値のη乗をパワースペクトルと見做すことにより推定されたスペクトル包絡で上記周波数領域サンプル列を除算した系列である白色化スペクトル系列のヒストグラムを近似する一般化ガウス分布の形状パラメータとし、η1はパラメータηの所定の値であるとして、
線形予測分析部が、時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る線形予測分析ステップと、
適合部が、Ｎ種類（Ｎは１以上の整数）のパラメータηのそれぞれに対応するＮ個の符号帳が記憶され、各符号帳にはそれぞれのパラメータηに対応する線形予測係数に変換可能な係数の候補が複数個格納された符号帳記憶部に記憶された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補と、上記線形予測分析ステップが得た線形予測係数に変換可能な係数と、のηの値を適合させる適合ステップと、
符号化部が、上記ηの値が適合された線形予測係数に変換可能な係数の複数個の候補と線形予測係数に変換可能な係数とを用いて、上記線形予測分析部が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得る符号化ステップと、
を含む線形予測符号化方法。
【請求項２７】
パラメータηを正の数として、時系列信号に対応するパラメータηを、その時系列信号に対応する周波数領域サンプル列の絶対値のη乗をパワースペクトルと見做すことにより推定されたスペクトル包絡で上記周波数領域サンプル列を除算した系列である白色化スペクトル系列のヒストグラムを近似する一般化ガウス分布の形状パラメータとし、η1はパラメータηの所定の値であるとして、
時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る線形予測分析ステップ、
入力されたη1に基づいて、符号帳記憶部に記憶された符号帳と上記線形予測係数に変換可能な係数との少なくとも一方を適合させる適合ステップと、
上記符号帳又は上記適合された符号帳を用いて、上記線形予測係数に変換可能な係数又は上記適合された線形予測係数に変換可能な係数を符号化する符号化ステップと、
を含む線形予測符号化方法。
【請求項２８】
η1を正の数として、入力されたη1に基づいて、符号帳記憶部に記憶された符号帳と、上記符号帳に格納された複数個の線形予測係数に変換可能な係数の候補のうち、入力された線形予測係数符号に対応する線形予測係数に変換可能な係数の候補との少なくとも一方を適合させる適合ステップを含み、
上記線形予測係数に変換可能な係数は、上記線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列を得るために用いられる、
線形予測復号方法。
【請求項２９】
請求項１から１５の何れかの線形予測符号化装置又は請求項１６から２５の何れかの線形予測復号装置の各部としてコンピュータを機能させるためのプログラム。
【請求項３０】
請求項１から１５の何れかの線形予測符号化装置又は請求項１６から２５の何れかの線形予測復号装置の各部としてコンピュータを機能させるためのプログラムが記録されたコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、線形予測係数に変換可能な係数を符号化又は復号する技術に関する。
【背景技術】
【０００２】
線形予測係数に変換可能な係数の１つであるＬＳＰパラメータの量子化技術として、ベクトル量子化等の手法が知られている（例えば、非特許文献１参照）。【０００３】
ところで、公知とはなっていないが、発明者によりパラメータηが提案されている。このパラメータηは、例えば3GPP EVS（Enhanced Voice Services）規格で使われているような線形予測包絡を利用する周波数領域の係数の量子化値を算術符号化する符号化方式において、算術符号の符号化対象の属する確率分布を定める形状パラメータである。パラメータηは、符号化対象の分布と関連性を有しており、パラメータηを適宜定めると効率の良い符号化及び復号を行うことが可能である。【０００４】
また、パラメータηは、時系列信号の特徴を表す指標と成り得る。このため、パラメータηを適宜用いると、ＬＳＰパラメータ等の線形予測係数に変換可能な係数を効率の良く符号化及び復号を行うことが可能である。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】守谷健弘，「高圧縮音声符号化の必須技術：線スペクトル対（LSP）」，NTT技術ジャーナル，２０１４年９月，Ｐ．５８−６０
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、パラメータηを用いた線形予測係数に変換可能な係数の符号化及び復号技術は知られていなかった。【０００７】
本発明は、パラメータηを用いて線形予測係数に変換可能な係数の符号化又は復号を行う線形予測符号化装置、線形予測復号装置、これらの方法、プログラム及び記録媒体を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の一態様による線形予測符号化装置によれば、パラメータηを正の数として、時系列信号に対応するパラメータηを、その時系列信号に対応する周波数領域サンプル列の絶対値のη乗をパワースペクトルと見做すことにより推定されたスペクトル包絡で周波数領域サンプル列を除算した系列である白色化スペクトル系列のヒストグラムを近似する一般化ガウス分布の形状パラメータとし、η1はパラメータηの所定の値であるとして、時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る線形予測分析部と、Ｎ種類（Ｎは１以上の整数）のパラメータηのそれぞれに対応するＮ個の符号帳が記憶され、各符号帳にはそれぞれのパラメータηに対応する線形予測係数に変換可能な係数の候補が複数個格納された符号帳記憶部と、符号帳記憶部に記憶された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補と、線形予測分析部が得た線形予測係数に変換可能な係数と、のηの値を適合させる適合部と、ηの値が適合された線形予測係数に変換可能な係数の複数個の候補と線形予測係数に変換可能な係数とを用いて、線形予測分析部が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得る符号化部と、を備えている。【０００９】
本発明の一態様による線形予測符号化装置によれば、パラメータηを正の数として、時系列信号に対応するパラメータηを、その時系列信号に対応する周波数領域サンプル列の絶対値のη乗をパワースペクトルと見做すことにより推定されたスペクトル包絡で周波数領域サンプル列を除算した系列である白色化スペクトル系列のヒストグラムを近似する一般化ガウス分布の形状パラメータとし、η1はパラメータηの所定の値であるとして、時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る線形予測分析部と、符号帳が記憶された符号帳記憶部と、入力されたη1に基づいて、符号帳記憶部に記憶された符号帳と線形予測係数に変換可能な係数との少なくとも一方を適合させる適合部と、符号帳又は適合された符号帳を用いて、線形予測係数に変換可能な係数又は適合された線形予測係数に変換可能な係数を符号化する符号化部と、を備えている。【００１０】
本発明の一態様による線形予測復号装置によれば、符号帳が記憶された符号帳記憶部と、η1を正の数として、入力されたη1に基づいて、符号帳記憶部に記憶された符号帳と、符号帳に格納された複数個の線形予測係数に変換可能な係数の候補のうち、入力された線形予測係数符号に対応する線形予測係数に変換可能な係数の候補との少なくとも一方を適合させる適合部を備えており、線形予測係数に変換可能な係数は、線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列を得るために用いられる。
【発明の効果】
【００１１】
パラメータηを用いて線形予測係数に変換可能な係数の符号化又は復号を行うことができる。
【図面の簡単な説明】
【００１２】
【図１】線形予測符号化装置の例を説明するためのブロック図。
【図２】線形予測符号化装置の例を説明するためのブロック図。
【図３】線形予測符号化装置の例を説明するためのブロック図。
【図４】線形予測符号化方法の例を説明するためのフローチャート。
【図５】ＬＳＰパラメータとηとの関係の例を説明するための図。
【図６】線形予測復号装置の例を説明するためのブロック図。
【図７】線形予測復号方法の例を説明するためのフローチャート。
【図８】符号化装置の例を説明するためのブロック図。
【図９】符号化方法の例を説明するためのフローチャート。
【図１０】符号化部の例を説明するためのブロック図。
【図１１】符号化部の例を説明するためのブロック図。
【図１２】符号化部の処理の例を説明するためのフローチャート。
【図１３】復号装置の例を説明するためのブロック図。
【図１４】復号方法の例を説明するためのフローチャート。
【図１５】復号部の処理の例を説明するためのフローチャート。
【図１６】符号化装置の例を説明するためのブロック図。
【図１７】符号化方法の例を説明するためのフローチャート。
【図１８】パラメータ決定装置の例を説明するためのブロック図。
【図１９】パラメータ決定方法の例を説明するためのフローチャート。
【図２０】一般化ガウス分布を説明するための図。
【図２１】線形予測符号化装置の例を説明するためのブロック図。
【図２２】線形予測符号化方法の例を説明するためのフローチャート。
【図２３】線形予測復号装置の例を説明するためのブロック図。
【図２４】線形予測復号方法の例を説明するためのフローチャート。
【図２５】線形予測符号化装置の例を説明するためのブロック図。
【図２６】線形予測符号化装置の例を説明するためのブロック図。
【図２７】線形予測符号化装置の例を説明するためのブロック図。
【図２８】線形予測復号装置の例を説明するためのブロック図。
【発明を実施するための形態】
【００１３】
［線形予測符号化装置、線形予測復号装置及びこれらの方法］
以下、線形予測符号化装置、線形予測復号装置及びこれらの方法を用いた符号化装置、復号装置及びこれらの方法の例について説明する。
【００１４】
［線形予測符号化装置、線形予測復号装置及びこれらの方法の第一実施形態］
（符号化）
第一実施形態の線形予測符号化装置及び方法の一例について説明する。
【００１５】
第一実施形態の線形予測符号化装置は、図１、図２又は図３に示すように、線形予測分析部２２１、符号帳記憶部２２２、符号化部２２４及び線形変換部２２５を例えば備えている。図１、図２又は図３の例では線形予測符号化装置の外部に周波数領域変換部２２０が設けられているが、線形予測符号化装置が周波数領域変換部２２０を更に備えていてもよい。線形予測符号化装置の各部が、図４に例示する各処理を行うことにより線形予測符号化方法が実現される。
【００１６】
＜周波数領域変換部２２０＞
周波数領域変換部２２０には、時系列信号である時間領域の音信号が入力される。
【００１７】
周波数領域変換部４１は、所定の時間長のフレーム単位で、入力された時間領域の音信号を周波数領域のN点のMDCT係数列X(0),X(1),…,X(N-1)に変換する。Nは正の整数である。
【００１８】
得られたMDCT係数列X(0),X(1),…,X(N-1)は、線形予測分析部２２１に出力される。
【００１９】
特に断りがない限り、以降の処理はフレーム単位で行われるものとする。
【００２０】
このようにして、周波数領域変換部２２０は、時系列信号に対応する、例えばMDCT係数列である周波数領域サンプル列を求める。
【００２１】
＜線形予測分析部２２１＞
線形予測分析部２２１には、例えばMDCT係数列X(0),X(1),…,X(N-1)である周波数領域サンプル列及びその周波数領域サンプル列に対応するパラメータη1が入力される。
【００２２】
パラメータη1は、正の数である。パラメータη1は、例えば、後述するパラメータ決定部２７，２７’により決定される。パラメータη1は、例えば3GPP EVS（Enhanced Voice Services）規格で使われているような線形予測包絡を利用する周波数領域の係数の量子化値を算術符号化する符号化方式において、算術符号の符号化対象の属する確率分布を定めるパラメータηである。パラメータηは、時系列信号の特徴を表す指標と成り得るものである。後に出てくるパラメータη2,η3も、パラメータηである。η1,η2,η3は、パラメータηの所定の値とも言える。
【００２３】
なお、パラメータη1についての情報は、線形予測復号装置に送信されるとする。例えば、パラメータη1を表すパラメータ符号が線形予測復号装置に送信される。
【００２４】
線形予測分析部２２１は、MDCT係数列X(0),X(1),…,X(N-1)及びη1を用いて、以下の式（A7)により定義される~R(0),~R(1),…,~R(N-1)を用いて線形予測分析を行って線形予測係数係数に変換可能な係数を生成する（ステップＤＥ１）。
【００２５】
【数１】
2016167215.tif 000003
【００２６】
生成された線形予測係数係数に変換可能な係数は、符号化部２２４に出力される。
【００２７】
具体的には、線形予測分析部２２は、まずMDCT係数列X(0),X(1),…,X(N-1)の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換に相当する演算、すなわち式(A7)の演算を行うことにより、MDCT係数列X(0),X(1),…,X(N-1)の絶対値のη1乗に対応する時間領域の信号列である擬似相関関数信号列~R(0),~R(1),…,~R(N-1)を求める。そして、線形予測分析部２２は、求まった擬似相関関数信号列~R(0),~R(1),…,~R(N-1)を用いて線形予測分析を行って、線形予測係数係数に変換可能な係数を生成する。
【００２８】
このようにして、線形予測分析部２２１は、η1を正の数として、時系列信号に対応する周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る。
【００２９】
線形予測係数に変換可能な係数とは、例えばLSP、PARCOR係数、ISP等である。線形予測係数に変換可能な係数は、線形予測係数自体であってもよい。
【００３０】
pを所定の正の数とし、線形予測係数に可能な係数の次数をp次とする。
【００３１】
＜符号帳記憶部２２２＞
符号帳記憶部２２２には、パラメータη2に対応する線形予測係数に変換可能な係数の候補が複数個格納された符号帳が記憶されている。
【００３２】
以下、線形予測係数に変換可能な係数の候補と、その線形予測係数に変換可能な係数の候補に対応する符号とのペアを、候補符号ペアと呼ぶことにする。符号帳には、複数個の候補符号ペアが記憶されている。言い換えると、Nを所定の２以上の数とすると、符号帳には、N個の候補ペアが記憶されている。線形予測係数に変換可能な係数の候補に対応する符号のそれぞれには、所定の数のビットが割り当てられている。各符号は、割り当てられた所定の数のビットで表現される。
【００３３】
線形予測係数に変換可能な係数の次数がpであるため、線形予測係数に変換可能な係数の各候補はp個の値から構成される。
【００３４】
パラメータη2に対応する線形予測係数に変換可能な係数の候補とは、パラメータηの値がη2である周波数領域サンプル列に対応する線形予測係数に変換可能な係数を符号化するために最適化された線形予測係数に変換可能な係数の候補である。
【００３５】
＜線形変換部２２５＞
線形変換部２２５には、線形予測分析部２２１が得た線形予測係数に変換可能な係数と、その線形予測係数に変換可能な係数に対応するパラメータη1とが入力される。パラメータη1は、例えば、後述するパラメータ決定部２７，２７’により決定される。
【００３６】
線形変換部２２５は、第一線形変換部２２５１及び第二線形変換部２２５２の少なくとも一方を備えている。
【００３７】
以下、（１）図１に示すように線形変換部２２５が第一線形変換部２２５１を備えている場合を第１の場合とし、（２）図２に示すように線形変換部２２５が第二線形変換部２２５２を備えている場合を第２の場合とし、（３）図３に示すように線形変換部２２５が第一線形変換部２２５１及び第二線形変換部２２５２を備えている場合を第３の場合として、各場合について説明する。
【００３８】
（１）第１の場合
この場合、線形変換部２２５の第一線形変換部２２５１は、符号帳記憶部２２２に記憶された線形予測係数に変換可能な係数の候補に対し、少なくとも入力されたパラメータη1に応じた第一線形変換を行う（ステップＤＥ２）。
【００３９】
例えば、第一線形変換部２２５１は、入力されたパラメータη1と符号帳記憶部２２２に格納された線形予測係数に変換可能な係数の候補に対応するパラメータη2とに応じた第一線形変換により、符号帳記憶部２２２から読み込んだパラメータη2に対応する線形予測係数に変換可能な係数の候補を、パラメータη1に対応する線形予測係数に変換可能な係数の候補に変換する。
【００４０】
パラメータη1に対応する線形予測係数に変換可能な係数の候補とは、パラメータηの値がη1である周波数領域サンプル列に対応する線形予測係数に変換可能な係数を符号化するために最適化された線形予測係数に変換可能な係数の候補である。
【００４１】
第一線形変換後の線形予測係数に変換可能な係数の候補は、符号化部２２４に出力される。
【００４２】
なお、パラメータη1の値とパラメータη2の値とが同一である場合には、第一線形変換部２２５１は、第一線形変換をしなくてもよい。
【００４３】
また、例えば、線形変換部２２５の第一線形変換部２２５１は、入力されたパラメータη1に応じて、入力されたパラメータη1が小さいほど、第一線形変換後の線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列が平坦になるように、符号帳記憶部２２２から読み込んだ線形予測係数に変換可能な係数の候補に対して第一線形変換を行い、変換後の線形予測係数に変換可能な係数の候補を出力する。
【００４４】
一般にパラメータηが小さいほど、非平滑化スペクトル包絡系列は平坦になる傾向があり、線形予測係数に変換可能な係数はより同じような値を取る傾向がある。例えば線形予測係数に変換可能な係数がLSPである場合には、パラメータηが小さいほど、LSPである線形予測係数に変換可能な係数は０からπまでを均等分割した値により近づく傾向がある。
【００４５】
図５に、パラメータηが各値を取るときのLSPパラメータの値の例を示す。図５の横軸はパラメータηであり、縦軸はLSPパラメータである。図５をみると、パラメータηが小さいほどLSPパラメータは０からπまでを均等分割した値に近づく傾向があることがわかる。
【００４６】
この傾向を用いて、パラメータη１が小さいほど、非平滑化スペクトル包絡系列がより平坦な場合に対応するように線形予測係数に変換可能な係数の候補を変換したものを用いて符号化及び復号を行うことにより量子化性能を向上させることができる。
【００４７】
（２）第２の場合
この場合、線形変換部２２５の第二線形変換部２２５２は、線形予測分析部２２１で得られた線形予測係数に変換可能な係数に対し、少なくとも入力されたパラメータη1に応じた第二線形変換を行う（ステップＤＥ２）。
【００４８】
例えば、第二線形変換部２２５２は、線形予測分析部２２１で得られたパラメータη1に対応する線形予測係数に変換可能な係数を、符号帳記憶部２２２に格納された線形予測係数に変換可能な係数の候補に対応するようにするために、パラメータη2に対応する線形予測係数に変換可能な係数に、第二線形変換する。
【００４９】
第二線形変換後の線形予測係数に変換可能な係数は、符号化部２２４に出力される。
【００５０】
なお、パラメータη1の値とパラメータη2の値とが同一である場合には、第二線形変換部２２５２は、第二線形変換をしなくてもよい。
【００５１】
または、例えば、線形変換部２２５の第二線形変換部２２５２は、入力されたパラメータη1に応じて、入力されたパラメータη1が小さいほど、第二線形変換後の線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列が平坦になるように、入力された線形予測係数に変換可能な係数に対して第二線形変換を行い、変換後の線形予測係数に変換可能な係数を出力する。
【００５２】
（３）第３の場合
この場合、線形変換部２２５の第一線形変換部２２５１は、符号帳記憶部２２２に記憶された線形予測係数に変換可能な係数の候補に対し、少なくともパラメータη3に応じた第一線形変換を行う。パラメータη3は、正の値であり、パラメータη2とは異なる値を予め定めておくか、線形予測係数符号化装置の外部から入力されるものである。
【００５３】
例えば、第一線形変換部２２５１は、パラメータη3と符号帳記憶部２２２に格納された線形予測係数に変換可能な係数の候補に対応するパラメータη2とに応じた第一線形変換により、符号帳記憶部２２２から読み込んだパラメータη2に対応する線形予測係数に変換可能な係数の候補を、パラメータη3に対応する線形予測係数に変換可能な係数の候補に変換する。
【００５４】
パラメータη3に対応する線形予測係数に変換可能な係数の候補とは、パラメータηの値がη3である周波数領域サンプル列に対応する線形予測係数に変換可能な係数を符号化するために最適化された線形予測係数に変換可能な係数の候補である。
【００５５】
第一線形変換後の線形予測係数に変換可能な係数の候補は、符号化部２２４に出力される。
【００５６】
なお、パラメータη2の値とパラメータη3の値とが同一である場合には、第一線形変換部２２５１は、第一線形変換をしなくてもよい。
【００５７】
また、例えば、線形変換部２２５の第一線形変換部２２５１は、パラメータη3が小さいほど、第一線形変換後の線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡が平坦になるように、符号帳記憶部２２２から読み込んだ線形予測係数に変換可能な係数の候補に対して第一線形変換を行い、変換後の線形予測係数に変換可能な係数の候補を出力する。
【００５８】
また、この第３の場合、線形変換部２２５の第二線形変換部２２５２は、線形予測分析部２２１で得られた線形予測係数に変換可能な係数に対し、少なくともパラメータη1に応じた第二線形変換を行う。
【００５９】
例えば、第二線形変換部２２５２は、線形予測分析部２２１で得られたパラメータη1に対応する線形予測係数に変換可能な係数を、パラメータη3に対応する線形予測係数に変換可能な係数に、第二線形変換する。
【００６０】
第二線形変換後の線形予測係数に変換可能な係数の候補は、符号化部２２４に出力される。
【００６１】
なお、パラメータη1の値とパラメータη3の値とが同一である場合には、第二線形変換部２２５２は、第二線形変換をしなくてもよい。
【００６２】
または、例えば、線形変換部２２５の第二線形変換部２２５２は、入力されたパラメータη1に応じて、入力されたパラメータη1が小さいほど、第二線形変換後の線形予測係数に変換可能な係数に対応する振幅スペクトル包絡が平坦になるように、入力された線形予測係数に変換可能な係数に対して第二線形変換を行い、変換後の線形予測係数に変換可能な係数を出力する。
【００６３】
このようにして、（３）第３の場合には、線形変換部２２５は、符号帳記憶部２２２に記憶された線形予測係数に変換可能な係数の候補に対する、η3に応じた第一線形変換と、線形予測分析部２２１で得られた線形予測係数に変換可能な係数に対する、η3に応じた第二線形変換との少なくとも一方を行う（ステップＤＥ２）。
【００６４】
＜符号化部２２４＞
符号化部２２４の処理は、線形変換部２２５の構成に応じて異なる。このため、線形変換部２２５が（１）第１の場合、（２）第２の場合及び（３）第３の場合のそれぞれ場合の符号化部２２４の処理について以下に説明する。
【００６５】
（１）第１の場合
線形変換部２２が（１）第１の場合には、符号化部２２４には、線形予測分析部２２１が得た線形予測係数に変換可能な係数と、線形変換部２２５の第一線形変換部２２５１が得た第一線形変換後の線形予測係数に変換可能な係数の候補とが入力される。
【００６６】
符号化部２２４は、線形予測係数に変換可能な係数について、第一線形変換後の線形予測係数に変換可能な係数の候補を用いて符号化して線形予測係数符号を得る（ステップＤＥ３）。
【００６７】
具体的には、符号化部２２４は、複数個の、第一線形変換後の線形予測係数に変換可能な係数の候補の中で、線形予測係数に変換可能な係数に最も近いものを選択し、その選択された候補に対応する符号を線形予測係数符号とする。
【００６８】
得られた線形予測係数符号は、復号装置に出力される。
【００６９】
（２）第２の場合
線形変換部２２が（２）第２の場合には、符号化部２２４には、線形予測分析部２２１の第二線形変換部２２５２が得た線形予測係数に変換可能な係数と、符号帳記憶部２２２に記憶された線形予測係数に変換可能な係数の候補とが入力される。
【００７０】
符号化部２２４は、第二線形変換後の線形予測係数に変換可能な係数について、線形予測係数に変換可能な係数の候補を用いて符号化して線形予測係数符号を得る（ステップＤＥ３）。
【００７１】
具体的には、符号化部２２４は、複数個の、線形予測係数に変換可能な係数の候補の中で、第二線形変換後の線形予測係数に変換可能な係数に最も近いものを選択し、その選択された候補に対応する符号を線形予測係数符号とする。
【００７２】
得られた線形予測係数符号は、復号装置に出力される。
【００７３】
（３）第３の場合
線形変換部２２が（３）第３の場合には、符号化部２２４には、線形予測分析部２２１の第二線形変換部２２５２が得た線形予測係数に変換可能な係数と、線形予測分析部２２１の第一線形変換部２２５１が得た線形予測係数に変換可能な係数の候補とが入力される。
【００７４】
符号化部２２４は、第二線形変換後の線形予測係数に変換可能な係数について、第一線形変換後の線形予測係数に変換可能な係数の候補を用いて符号化して線形予測係数符号を得る（ステップＤＥ３）。
【００７５】
具体的には、符号化部２２４は、複数個の、第一線形変換後の線形予測係数に変換可能な係数の候補の中で、第二線形変換後の線形予測係数に変換可能な係数に最も近いものを選択し、その選択された候補に対応する符号を線形予測係数符号とする。
【００７６】
得られた線形予測係数符号は、復号装置に出力される。
【００７７】
このように、線形予測係数に変換可能な係数を線形予測係数に変換可能な係数の候補を用いて符号化する際に、線形予測係数に変換可能な係数に対応するパラメータηと線形予測係数に変換可能な係数の候補に対応するパラメータηとが同じ値または近い値となるように、線形予測係数に変換可能な係数と線形予測係数に変換可能な係数の候補の少なくとも何れかに対して線形変換を行ったものを符号化に用いることにより、符号化歪を小さくすることができる及び／又は線形予測係数符号の符号量を小さくすることができる。
【００７８】
（復号）
第一実施形態の線形予測復号装置及び方法の一例について説明する。
【００７９】
第一実施形態の線形予測復号装置は、図６に示すように、符号帳記憶部３１１、復号部３１３及び線形変換部３１４を例えば備えている。線形予測復号装置の各部が、図７に例示する各処理を行うことにより線形予測復号方法が実現される。
【００８０】
＜符号帳記憶部３１１＞
符号帳記憶部３１１には、符号帳記憶部２２２に記憶されている符号帳と同じ符号帳が記憶されている。すなわち、符号帳記憶部３１１には、パラメータη2に対応する線形予測係数に変換可能な係数の候補が複数個格納された符号帳が記憶されている。
【００８１】
＜復号部３１３＞
復号部３１３には、線形予測符号化装置が出力した線形予測係数符号が入力される。
【００８２】
復号部３１３は、符号帳記憶部３１１に記憶された複数個の線形予測係数に変換可能な係数の候補のうち、入力された線形予測係数符号に対応する線形予測係数に変換可能な係数の候補を線形予測係数に変換可能な係数として得る（ステップＤＤ１）。
【００８３】
得られた線形予測係数に変換可能な係数は、線形変換部３１４に出力される。
【００８４】
得られた線形予測係数に変換可能な係数は、符号帳記憶部３１１に記憶されたパラメータη2に対応する複数個の線形予測係数に変換可能な係数の候補の何れか１つである。このため、復号部３１３で得られた線形予測係数に変換可能な係数は、パラメータη2に対応する線形予測係数に変換可能な係数となる。
【００８５】
＜線形変換部３１４＞
線形変換部３１４には、復号部３１３で得られたパラメータη2に対応する線形予測係数に変換可能な係数と、パラメータη1とが入力される。このパラメータη1は、例えば線形予測符号化装置から受信したパラメータ符号を復号することにより得られるものである。
【００８６】
線形変換部３１４は、パラメータη2に対応する線形予測係数に変換可能な係数に対して、少なくともパラメータη1に応じた線形変換をして線形変換後の線形予測係数に変換可能な係数を得る。
【００８７】
例えば、線形変換部３１４は、入力されたパラメータη1と線形予測係数に変換可能な係数に対応するパラメータη2とに応じた線形変換により、パラメータη2に対応する線形予測係数に変換可能な係数を、パラメータη1に対応する線形予測係数に変換可能な係数に変換する。
【００８８】
得られた線形変換後の線形予測係数に変換可能な係数は、線形予測復号装置又は方法による復号結果として出力される。
【００８９】
なお、パラメータη1の値とパラメータη2の値とが同一である場合には、線形変換部３１４は、線形変換をしなくてもよい。
【００９０】
また、線形変換部３１４は、パラメータη2に対応する線形予測係数に変換可能な係数を線形変換してパラメータη1に対応する線形予測係数に変換可能な係数を得る際に、パラメータη1ともパラメータη2とも異なるパラメータη4を用いて、線形変換を複数回行う構成としてもよい。
【００９１】
例えば、線形変換を２回行う場合について説明する。この場合、線形変換部３１４は、パラメータη2に対応する線形予測係数に変換可能な係数を線形変換してパラメータη4に対応する線形予測係数に変換可能な係数を得る。また、線形変換部３１４は、得られたパラメータη4に対応する線形予測係数に変換可能な係数を線形変換してパラメータη1に対応する線形予測係数に変換可能な係数を得る。ここで、パラメータη4を線形予測係数符号化装置が用いたパラメータη3と同一の値とすれば、２つの線形変換に、線形予測係数符号化装置の線形変換部２２５の第３の場合におけるパラメータη2に対応する線形予測係数に変換可能な係数の候補からパラメータη3に対応する線形予測係数に変換可能な係数の候補を得る線形変換と、線形予測係数符号化装置の線形変換部２２５の第３の場合におけるパラメータη1に対応する線形予測係数に変換可能な係数をパラメータη3に対応する線形予測係数に変換可能な係数を得る線形変換と、同一の線形変換を用いることができる。
【００９２】
なお、線形変換部３１４は、パラメータη2からパラメータη3への線形変換と、パラメータη3からパラメータη1への線形変換とを合成した１つの線形変換を、パラメータη2に対応する線形予測係数に変換可能な係数に対してすることにより、パラメータη1に対応する線形予測係数に変換可能な係数を得てもよい。
得られたパラメータη1に対応する線形予測係数に変換可能な係数は、線形予測復号装置又は方法による復号結果として出力される。
【００９３】
また、例えば、線形変換部３１４は、線形予測符号化装置の線形変換部２２５と同様に、入力されたη1が小さいほど、線形変換後の線形予測係数に変換可能な係数に対応する振幅スペクトル包絡が平坦になるように、復号部３１３で得られた線形予測係数に変換可能な係数を線形変換して線形変換後の線形予測係数に変換可能な係数を得てもよい。
【００９４】
これは、一般にパラメータηが小さいほど、非平滑化スペクトル包絡系列は平坦になるという傾向に基づくものである。
【００９５】
線形変換部３１４で得られた線形変換後の線形予測係数に変換可能な係数は、線形変換部３１４で得られた線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列を得るために用いられる。
【００９６】
［線形変換］
以下、第一線形変換及び第二線形変換等の線形変換の例について説明する。
【００９７】
線形変換前の線形予測係数に変換可能な係数又は線形予測係数に変換可能な係数の候補を^ω[k][k=1,2,…,p]とし、線形変換後の線形予測係数に変換可能な係数又は上記線形予測係数に変換可能な係数の候補を~ω[k][k=1,2,…,p]とする。また、線形変換前の線形予測係数に変換可能な係数はLSPであるとする。このとき、第一線形変換部２２５１、第二線形変換部２２５２、逆線形変換部２２６及び線形変換部３１４は、例えば以下の式に示される線形変換を行う。
【００９８】
【数２】
2016167215.tif 000004
【００９９】
ここで、x1,x2,…xp,y1,y2,…yp-1,z2,z3,…zpを所定の非負の数とし、y1,y2,…yp-1,z2,z3,…zpの少なくとも１つは所定の正の数であるとし、Kをx1,x2,…xp,y1,y2,…yp-1,z2,z3,…zp以外の要素が０である行列とする。
【０１００】
x1,x2,…xp,y1,y2,…yp-1,z2,z3,…zpの具体的な値は、線形変換前の線形予測係数に変換可能な係数又は線形予測係数に変換可能な係数の候補に対応するパラメータη（以下、線形変換前パラメータηAとする）の値と、線形変換後の線形予測係数に変換可能な係数又は線形予測係数に変換可能な係数の候補に対応するパラメータη（以下、線形変換後パラメータηBとする）の値とに基づいて適宜定まるものである。
【０１０１】
異なる複数の、線形変換前パラメータηAと線形変換後パラメータηBとの組に対応するx1,x2,…xp,y1,y2,…yp-1,z2,z3,…zpの具体的な値を図示していない記憶部に予め記憶しておく。第一線形変換部２２５１、第二線形変換部２２５２、逆線形変換部２２６及び線形変換部３１４は、線形変換をするときに、その線形変換における線形変換前パラメータηAと線形変換後パラメータηBとの組に対応するx1,x2,…xp,y1,y2,…yp-1,z2,z3,…zpの具体的な値を読み込み、読み込んだこれらの値を用いて上記式による線形変換を行えばよい。
【０１０２】
ところで、パラメータη1が大きい場合には、線形予測係数に変換可能な係数を使って計算したスペクトル包絡の変動は大きい傾向がある。このため、次数が大きい線形予測係数に変換可能な係数の候補を用いて符号化及び復号をすることが望ましい。
【０１０３】
逆に、パラメータη1が小さい場合には、線形予測係数に変換可能な係数を使って計算したスペクトル包絡の変動は小さい傾向がある。このため、次数が小さい線形予測係数に変換可能な係数の候補を用いて符号化及び復号をしても量子化歪は小さいため符号化及び復号の精度はそれほど悪くならない。
【０１０４】
このため、線形変換部２２５の第一線形変換部２２５１は、パラメータη1が小さいほど第一線形変換後の線形予測係数に変換可能な係数の候補の次数が小さくなるように第一線形変換を行ってもよい。
【０１０５】
同様に、線形変換部３１４は、パラメータη1が小さいほど線形変換後の線形予測係数に変換可能な係数の次数が小さくなるように線形変換を行ってもよい。
【０１０６】
このように、線形変換前の線形変換前の線形予測係数に変換可能な係数又は線形予測係数に変換可能な係数の候補の次数と、線形変換後の線形予測係数に変換可能な係数又は線形予測係数に変換可能な係数の候補の次数とが異なるように線形変換が行われてもよい。
【０１０７】
なお、第一線形変換部２２５１は、線形変換前の次数と線形変換後の次数とが同じである線形変換を行った後に線形変換後の線形予測係数に変換可能な係数の候補の次数を減らしてもよい。また、第一線形変換部２２５１は、線形変換後の線形予測係数に変換可能な係数の候補の次数を減らした後に線形変換前の次数と線形変換後の次数とが同じである線形変換を行ってもよい。
【０１０８】
同様に、線形変換部３１４は、線形変換前の次数と線形変換後の次数とが同じである線形変換を行った後に線形変換後の線形予測係数に変換可能な係数の次数を減らしてもよい。また、線形変換部３１４は、線形変換後の線形予測係数に変換可能な係数の次数を減らした後に線形変換前の次数と線形変換後の次数とが同じである線形変換を行ってもよい。
【０１０９】
また、第一線形変換部２２５１は、パラメータη1が小さい場合には、線形変換後の線形予測係数に変換可能な係数の複数の候補を統合することにより、パラメータη1が小さいほど線形変換後の線形予測係数に変換可能な係数の複数の候補数を減らしてもよい。
【０１１０】
［線形予測符号化装置、線形予測復号装置及びこれらの方法の第二実施形態］
（符号化）
第二実施形態の線形予測符号化装置及び方法の一例について説明する。
【０１１１】
第二実施形態の線形予測符号化装置は、図２１に示すように、線形予測分析部２２１、符号帳記憶部２２２、符号帳選択部２２３及び符号化部２２４を例えば備えている。図２１の例では線形予測符号化装置の外部に周波数領域変換部２２０が設けられているが、線形予測符号化装置が周波数領域変換部２２０を更に備えていてもよい。線形予測符号化装置の各部が、図２２に例示する各処理を行うことにより線形予測符号化方法が実現される。
【０１１２】
第二実施形態では、「パラメータη1」のことを「パラメータη」と表記する。
【０１１３】
＜周波数領域変換部２２０＞
周波数領域変換部２２０には、時系列信号である時間領域の音信号が入力される。
【０１１４】
周波数領域変換部４１は、所定の時間長のフレーム単位で、入力された時間領域の音信号を周波数領域のN点のMDCT係数列X(0),X(1),…,X(N-1)に変換する。Nは正の整数である。
【０１１５】
得られたMDCT係数列X(0),X(1),…,X(N-1)は、線形予測分析部２２１に出力される。
【０１１６】
特に断りがない限り、以降の処理はフレーム単位で行われるものとする。
【０１１７】
このようにして、周波数領域変換部２２０は、時系列信号に対応する、例えばMDCT係数列である周波数領域サンプル列を求める。
【０１１８】
＜線形予測分析部２２１＞
線形予測分析部２２１には、例えばMDCT係数列X(0),X(1),…,X(N-1)である周波数領域サンプル列及びその周波数領域サンプル列に対応するパラメータηが入力される。
【０１１９】
パラメータηは、正の数である。パラメータηは、例えば、後述するパラメータ決定部２７，２７’により決定される。パラメータηは、例えば3GPP EVS（Enhanced Voice Services）規格で使われているような線形予測包絡を利用する周波数領域の係数の量子化値を算術符号化する符号化方式において、算術符号の符号化対象の属する確率分布を定める形状パラメータである。パラメータηは、時系列信号の特徴を表す指標と成り得るものである。
【０１２０】
線形予測分析部２２１は、線形予測分析部２２は、MDCT係数列X(0),X(1),…,X(N-1)及びηを用いて、以下の式（A7)により定義される~R(0),~R(1),…,~R(N-1)を用いて線形予測分析行って線形予測係数に変換可能な係数を生成する（ステップＤＥ１）。
【０１２１】
【数３】
2016167215.tif 000005
【０１２２】
生成された線形予測係数に変換可能な係数は、符号化部２２４に出力される。
【０１２３】
具体的には、線形予測分析部２２は、まずMDCT係数列X(0),X(1),…,X(N-1)の絶対値のη乗をパワースペクトルと見做した逆フーリエ変換に相当する演算、すなわち式(A7)の演算を行うことにより、MDCT係数列X(0),X(1),…,X(N-1)の絶対値のη乗に対応する時間領域の信号列である擬似相関関数信号列~R(0),~R(1),…,~R(N-1)を求める。そして、線形予測分析部２２は、求まった擬似相関関数信号列~R(0),~R(1),…,~R(N-1)を用いて線形予測分析を行って、線形予測係数に変換可能な係数を生成する。
【０１２４】
このようにして、線形予測分析部２２１は、ηを正の数として、時系列信号に対応する周波数領域サンプル列の絶対値のη乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得る。
【０１２５】
線形予測係数に変換可能な係数とは、例えばLSP,PARCOR係数、ISP等である。線形予測係数に変換可能な係数は、線形予測係数自体であってもよい。
【０１２６】
pを所定の正の数とし、線形予測係数に可能な係数の次数をp次とする。
【０１２７】
＜符号帳記憶部２２２＞
符号帳記憶部２２２には、複数の符号帳が記憶されている。
【０１２８】
以下、線形予測係数に変換可能な係数の候補と、その線形予測係数に変換可能な係数の候補に対応する符号とのペアを、候補符号ペアと呼ぶことにする。各符号帳には、複数の候補符号ペアが記憶されている。言い換えると、Iを所定の２以上の数として、Niをiに応じて定まる所定の２以上の数とすると、符号帳i(i=1,2,…,I)のそれぞれには、Ni個の候補ペアが記憶されている。線形予測係数に変換可能な係数の候補に対応する符号のそれぞれには、所定の数のビットが割り当てられている。各符号は、割り当てられた所定の数のビットで表現される。
【０１２９】
線形予測係数に変換可能な係数の次数がpであるため、、線形予測係数に変換可能な係数の各候補はp個の値から構成される。
【０１３０】
符号帳記憶部２２２に記憶されている複数の符号帳は、符号帳選択部２２３の符号帳の選択方法によって異なる。このため、符号帳記憶部２２２に記憶されている複数の符号帳の例は、後述する符号帳選択部２２３の例と合わせて説明する。
【０１３１】
＜符号帳選択部２２３＞
符号帳選択部２２３には、パラメータηが入力される。
【０１３２】
符号帳選択部２２３は、符号帳記憶部２２２に記憶された複数の符号帳の中から入力されたηに応じて符号帳を選択する（ステップＤＥ２）。選択された符号帳についての情報は、符号化部２２４に出力される。
【０１３３】
以下、符号帳記憶部２２２に記憶された複数の符号帳の例及び符号帳選択部２２３による符号帳の選択基準の例について説明する。
【０１３４】
（１）第一の方法
第一の方法では、符号帳記憶部２２２には、線形予測係数に変換可能な係数の候補数が異なる複数の符号帳が記憶されている。また、符号帳選択部２２３は、パラメータηが大きいほど、符号帳記憶部２２２に記憶された複数の符号帳の中から、線形予測係数に変換可能な係数の候補数が多い符号帳を選択する。
【０１３５】
パラメータηが大きい場合には、線形予測係数に変換可能な係数の取り得る範囲は広い傾向があるため、線形予測係数に変換可能な係数を表現するために必要な線形予測係数に変換可能な係数の候補数は多くなる。このため、パラメータηが大きい場合には、線形予測係数に変換可能な係数の候補数が多い符号帳を用いて符号化及び復号をすることが望ましい。
【０１３６】
逆に、パラメータηが小さい場合には、線形予測係数に変換可能な係数の取り得る範囲は狭い傾向があるため、少ない個数の線形予測係数に変換可能な係数の候補で線形予測係数に変換可能な係数を表現することができる。このため、パラメータが小さい場合には、線形予測係数に変換可能な係数の候補数が少ない符号帳を用いて符号化及び復号をしても量子化歪は小さいため符号化及び復号の精度はそれほど悪くならない。
【０１３７】
このため、第一の方法では、符号帳選択部２２３は、パラメータηが大きいほど、符号帳記憶部２２２に記憶された複数の符号帳の中から、線形予測係数に変換可能な係数の候補数が多い符号帳を選択する。
【０１３８】
パラメータηの大きさについての判断は、言い換えれば適切な符号帳の選択は、閾値に基づいて行うことができる。例えば、第一符号帳の線形予測係数に変換可能な係数の候補数の方が、第二符号帳の線形予測係数に変換可能な係数の候補数よりもよりも少ないとする。この場合、パラメータηの閾値を１つ予め定めておき、入力されたパラメータηが閾値よりも小さい場合はパラメータηが小さいと判断し第一符号帳を選択する。入力されたパラメータηが閾値以上である場合はパラメータηが大きいと判断し第二符号帳を選択する。符号帳の数が３以上である場合には、符号帳の数から１を減算した値の個数の閾値を用いてこれと同様に符号帳を選択すればよい。
【０１３９】
なお、符号帳が多層構造を有しており、パラメータηに応じてどの層まで用いるのかを決定してもよい。例えば、p=16であり、16次の線形予測係数に変換可能な係数を２層の符号帳で符号化する例について説明する。この符号帳の第一層には10ビット、第二層には5ビットの量子化ビット数が割り当てられているとする。これにより、第一層には210=1024個の、線形予測係数に変換可能な係数の候補である16次元ベクトルとその候補に対応する符号とのペアが格納され、第二層には25=32個の、線形予測係数に変換可能な係数の候補である16次元ベクトルとその候補に対応する符号とのペアが格納されているとする。
【０１４０】
この場合、パラメータηが大きい場合には、第一層及び第二層を用いることにし、パラメータηが小さい場合には第一層のみを用いることにする。パラメータηが大きいか小さいかの判断は、上記と同様に閾値に基づいて行うことができる。
【０１４１】
パラメータηが大きい場合には、まず第一層の線形予測係数に変換可能な係数の候補の中で、入力された線形予測係数に変換可能な係数に最も近いもの及び対応する符号を選択する。次に選択された線形予測係数に変換可能な係数の候補の値を入力された線形予測係数に変換可能な係数から減算し、第二層の線形予測係数に変換可能な係数の候補の中で、その減算値と最も近いもの及び対応する符号を選択する。この場合、第一層及び第二層で選択された２個の符号が線形予測係数符号となる。すなわち、線形予測係数符号は15ビットで表現される。また、第一層及び第二層で選択された線形予測係数に変換可能な係数の候補の和が、入力された線形予測係数に変換可能な係数の量子化結果となる。
【０１４２】
パラメータηが小さい場合には、第一層の線形予測係数に変換可能な係数の候補の中で、入力された線形予測係数に変換可能な係数に最も近いもの及び対応する符号を選択する。この場合、第一層で選択された符号が線形予測係数符号となる。すなわち、線形予測係数符号は10ビットで表現される。また、第一層で選択された線形予測係数に変換可能な係数の候補が、入力された線形予測係数に変換可能な係数の量子化結果となる。
【０１４３】
第一層から構成される符号帳と、第一層及び第二層から構成される符号帳とを異なる符号帳と考えると、この例も（１）第一の方法の一例と言える。
【０１４４】
この多層構造を有する符号帳の例にように、１つの符号帳の中の候補符号ペアの数が可変である場合には、言い換えれば１つの符号帳の中の候補符号ペアの探索範囲が可変である場合には、パラメータηが小さいほど、候補符号ペアの探索範囲を狭くしてもよい。探索範囲が異なる候補符号ペアの集合を異なる符号帳と考えれば、この例も（１）第一の方法の一例と言える。
【０１４５】
（２）第二の方法
第二の方法では、符号帳記憶部２２２には、符号帳に記憶された線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η乗した系列である非平滑化スペクトル包絡系列の平坦度合いが異なる複数の符号帳が記憶されている。また、符号帳選択部２２３は、ηが小さいほど、符号帳記憶部２２２に記憶された複数の符号帳の中から、符号帳に記憶された線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η乗した系列である非平滑化スペクトル包絡系列がより平坦である符号帳を選択する。
【０１４６】
一般にパラメータηが小さいほど、非平滑化スペクトル包絡系列は平坦になる傾向があり、線形予測係数に変換可能な係数はより同じような値を取る傾向がある。例えば線形予測係数に変換可能な係数がLSPである場合には、パラメータηが小さいほど、LSPパラメータである線形予測係数に変換可能な係数は０からπまでを均等分割した値により近づく傾向がある。
【０１４７】
図５に、パラメータηが各値を取るときのLSPパラメータの値の例を示す。図５の横軸はパラメータηであり、縦軸はLSPパラメータである。図５をみると、パラメータηが小さいほどLSPパラメータは０からπまでを均等分割した値に近づく傾向があることがわかる。
【０１４８】
線形予測係数に変換可能な係数がISPパラメータの場合にも、同様の傾向がある。すなわち、線形予測係数に変換可能な係数がISPパラメータの場合、パラメータηが小さいほど、ISPパラメータである線形予測係数に変換可能な係数は０からπまでを均等分割した値により近づく傾向がある。
線形予測係数に変換可能な係数がPARCOR係数の場合には、パラメータηが小さいほど、PARCOR係数である線形予測係数に変換可能な係数は全体的に値が小さくなる傾向がある。
【０１４９】
第二の方法は、これらの傾向を用いて、パラメータηが小さいほど、非平滑化スペクトル包絡系列がより平坦な場合に対応する線形予測係数に変換可能な係数の候補を用いて符号化及び復号を行うことにより量子化性能を向上させようとするものである。
【０１５０】
線形予測係数に変換可能な係数がLSP又はPARCOR係数であるとして、符号帳i(i=1,2,…,I)の線形予測係数に変換可能な係数の候補を^ωn[1],^ωn[2],…,^ωn[p](n=1,2,…,Ni)と表記する。また、非平滑化スペクトル包絡が最も平坦な場合に対応する線形予測係数に変換可能な係数をωF[1],ωF[2],…,ωF[p]と表記する。
【０１５１】
この場合、第二の方法は、例えば、符号帳記憶部２２２には、以下のSi1の値が異なる複数の符号帳i(i=1,2,…,I)が記憶されているとし、符号帳選択部２２３が、ηが小さいほど、以下のSi1の値が小さい符号帳iを選択することにより実現される。
【０１５２】
Si1=(1/pNi)Σn=1NiΣk=1p|^ωn[k]-ωF[k]|
第二の方法においても、適切な符号帳の選択を閾値に基づいて行ってもよい。例えば、第一符号帳の線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η乗した系列である非平滑化スペクトル包絡系列の方が、第二符号帳の線形予測係数に変換可能な係数の候補に対応する振幅スペクトル包絡の系列を１／η乗した系列である非平滑化スペクトル包絡系列よりも平坦であるとする。この場合、パラメータηの閾値を１つ予め定めておき、入力されたパラメータηが閾値よりも小さい場合はパラメータηが小さいと判断し第一符号帳を選択する。入力されたパラメータηが閾値以上である場合はパラメータηが大きいと判断し第二符号帳を選択する。符号帳の数が３以上である場合には、符号帳の数から１を減算した値の個数の閾値を用いてこれと同様に符号帳を選択すればよい。
【０１５３】
（３）第三の方法
第三の方法では、符号帳記憶部２２２には、線形予測係数に変換可能な係数の候補間の間隔が異なる複数の符号帳が記憶されている。また、符号帳選択部２２３は、ηが小さいほど、符号帳記憶部２２２に記憶された複数の符号帳の中から、線形予測係数に変換可能な係数の候補間の間隔が狭い符号帳を選択する。
【０１５４】
線形予測係数に変換可能な係数の候補間の間隔とは、その符号帳に含まれる線形予測係数に変換可能な係数の候補間の間隔の広さを表す指標であればどのようなものであってもよい。例えば、線形予測係数に変換可能な係数の候補間の間隔は、その符号帳に含まれる、ある線形予測係数に変換可能な係数の候補と、別のある線形予測係数に変換可能な係数の候補との距離の平均値であってもよいし、その距離の最大値、最小値又は中央値であってもよい。
【０１５５】
第一の方法で述べたように、パラメータηが大きい場合には、線形予測係数に変換可能な係数の変動は大きい傾向がある。このため、線形予測係数に変換可能な係数の候補間の間隔が広い符号帳を用いて符号化及び復号をすることが望ましい。
【０１５６】
逆に、パラメータηが小さい場合には、線形予測係数に変換可能な係数の変動は小さい傾向がある。このため、線形予測係数に変換可能な係数の候補間の間隔が狭い符号帳を用いて符号化及び復号をしても量子化歪は小さいため符号化及び復号の精度はそれほど悪くならない。
【０１５７】
第三の方法は、この傾向を利用したものである。
【０１５８】
符号帳i(i=1,2,…,I)の線形予測係数に変換可能な係数の候補を^ωn[1],^ωn[2],…,^ωn[p](n=1,2,…,Ni)と表記する。
【０１５９】
この場合、第三の方法は、例えば、符号帳記憶部２２２には、以下のSi2の値が異なる複数の符号帳i(i=1,2,…,I)が記憶されているとし、符号帳選択部２２３が、ηが小さいほど、以下のSi2の値が小さい符号帳iを選択することにより実現される。
【０１６０】
Si2=(1/Ni)Σn=1Ni-1(Σk=1p(^ωn[k]-^ωn+1[k])|2)1/2
この例のように、また、線形予測係数に変換可能な係数の候補間の間隔は、その符号帳に含まれる、隣接する２個の線形予測係数に変換可能な係数の候補の距離の平均値であってもよい。
【０１６１】
第三の方法においても、適切な符号帳の選択を閾値に基づいて行ってもよい。例えば、第一符号帳の線形予測係数に変換可能な係数の候補間の間隔の方が、第二符号帳の線形予測係数に変換可能な係数の候補間の間隔よりも狭いとする。この場合、パラメータηの閾値を１つ予め定めておき、入力されたパラメータηが閾値よりも小さい場合はパラメータηが小さいと判断し第一符号帳を選択する。入力されたパラメータηが閾値以上である場合はパラメータηが大きいと判断し第二符号帳を選択する。符号帳の数が３以上である場合には、符号帳の数から１を減算した値の個数の閾値を用いてこれと同様に符号帳を選択すればよい。
【０１６２】
＜符号化部２２４＞
符号化部２２４には、線形予測分析部２２１が得た線形予測係数に変換可能な係数及び符号帳選択部２２３が得た選択された符号帳についての情報が入力される。
【０１６３】
符号化部２２４は、選択された符号帳を用いて、線形予測係数に変換可能な係数を符号化して線形予測係数符号を得る（ステップＤＥ３）。得られた線形予測係数符号は、復号装置に出力される。
【０１６４】
（復号）
第二実施形態の線形予測復号装置及び方法の一例について説明する。
【０１６５】
第二実施形態の線形予測復号装置は、図２３に示すように、符号帳記憶部３１１、符号帳選択部３１２及び復号部３１３を例えば備えている。線形予測復号装置の各部が、図２４に例示する各処理を行うことにより線形予測復号方法が実現される。
【０１６６】
第二実施形態では、「パラメータη1」のことを「パラメータη」と表記する。
【０１６７】
＜符号帳記憶部３１１＞
符号帳記憶部３１１には、複数の符号帳が記憶されている。
【０１６８】
以下、線形予測係数に変換可能な係数の候補と、その線形予測係数に変換可能な係数の候補に対応する符号とのペアを、候補符号ペアと呼ぶことにする。各符号帳には、複数の候補符号ペアが記憶されている。言い換えると、Iを所定の２以上の数として、Niをiに応じて定まる所定の２以上の数とすると、符号帳i(i=1,2,…,I)には、Ni個の候補ペアが記憶されている。線形予測係数に変換可能な係数の候補に対応する符号のそれぞれには、所定の数のビットが割り当てられている。各符号は、割り当てられた所定の数のビットで表現される。
【０１６９】
pを所定の正の数とし、線形予測係数に変換可能な係数の次元がpであるとすると、各線形予測係数に変換可能な係数の候補はp個の値から構成される。
【０１７０】
符号帳記憶部３１１に記憶されている複数の符号帳は、符号帳選択部３１２の符号帳の選択方法によって異なる。このため、符号帳記憶部３１１に記憶されている複数の符号帳の例は、後述する符号帳選択部３１２の例と合わせて説明する。
【０１７１】
なお、符号帳記憶部３１１には、符号帳記憶部２２２に記憶されている複数の符号帳と同じ符号帳が記憶されている。
【０１７２】
＜符号帳選択部３１２＞
符号帳選択部３１２には、パラメータηが入力される。パラメータηは、パラメータ符号を復号することにより得られる。パラメータηは、符号化装置及び復号装置で予め定められた同一の数であってもよい。
【０１７３】
符号帳選択部３１２は、符号帳記憶部３１１に記憶された複数の符号帳の中から入力されたηに応じて符号帳を選択する（ステップＤＤ１）。選択された符号帳についての情報は、復号部３１３に出力される。
【０１７４】
符号帳記憶部３１１には、符号帳記憶部２２２に記憶された複数の符号帳と同じ符号帳が記憶されているとする。また、符号帳選択部３１２には、符号化装置の符号帳選択部２２３による符号帳の選択基準と同じ選択基準が予め定められているとする。これにより、符号側で選択される符号帳と同じ内容の符号帳が復号側でも選択されることになる。
【０１７５】
符号帳の選択基準については、符号化側で説明したため、ここでは重複説明を省略する。
【０１７６】
＜復号部３１３＞
復号部３１３には、符号化装置が出力した線形予測係数符号及び符号帳選択部３１２が得た選択された符号帳についての情報が入力される。また、復号部３１３は、選択された符号帳についての情報により特定される符号帳を符号帳記憶部３１１により読み込む。
【０１７７】
復号部３１３は、選択された符号帳を用いて、線形予測係数符号を復号して線形予測係数に変換可能な係数を得る（ステップＤＤ２）。
【０１７８】
線形予測係数に変換可能な係数は、線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列を１／η乗した系列である非平滑化スペクトル包絡系列を得るために用いられる。
【０１７９】
［線形予測符号化装置、線形予測復号装置及びこれらの方法の変形例］
図１から図３、図２１及び図２５から図２７に一点鎖線で示すように、適合部２２Ａが符号帳選択部２２３及び線形変換部２２５の少なくとも一方から構成されているとすると、適合部２２Ａは、入力されたη1に基づいて、符号帳記憶部２２２に記憶された符号帳と、線形予測分析部２２１により生成された線形予測係数に変換可能な係数との少なくとも一方を適合させていると言える。言い換えれば、適合部２２Ａは、符号帳記憶部２２に記憶された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補と、線形予測分析部２２１が得た線形予測係数に変換可能な係数と、のηの値を適合させていると言える。適合部２２Ａは、例えば、適合前の「符号帳記憶部２２２に記憶されている符号帳、つまり線形予測係数に変換可能な係数の複数個の候補に対応するパラメータηの値と、線形予測分析部２２１により生成された線形予測係数に変換可能な係数に対応するパラメータηの値との差」に比べて、適合後の２つのパラメータηの値の差が小さくなるように、少なくとも一方の線形予測係数に変換可能な係数を変形しているとも言える。なお、適合部２２Ａは、適合後には２つのパラメータηの値がほぼ同じ値になるように適合を行っているとも言える。。第一実施形態で説明した線形変換部２２５の第一線形変換部２２５１の処理及び第二実施形態で説明した符号帳選択部２２３の処理は、符号帳記憶部２２２に記憶された符号帳の適合の一例である。第二実施形態で説明した線形変換部２２５の第二線形変換部２２５２の処理は、線形予測分析部２２１により生成された線形予測係数に変換可能な係数の適合の一例である。
【０１８０】
この場合、符号化部２２４は、適合部２２Ａにより適合された少なくとも一方の符号帳及び線形予測係数に変換可能な係数を用いて、符号化を行っていると言える。言い換えれば、符号化部２２４は、符号帳選択部２２３で選択された符号帳又は適合部２２Ａにより適合された符号帳を用いて、線形予測分析部２２１により線形予測係数に変換可能な係数又は適合部２２Ａにより適合された線形予測係数に変換可能な係数を符号化していると言える。さらに、言い換えれば、符号化部２２４は、ηの値が適合された線形予測係数に変換可能な係数の複数個の候補と線形予測係数に変換可能な係数とを用いて、線形予測分析部２２１が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得ていると言える。
【０１８１】
第一実施形態の（１）第１の場合の適合部２２Ａは、符号帳記憶部２２２に記憶された線形予測係数に変換可能な係数の候補に対して、η1に応じた第一線形変換を行い、第一線形変換後の線形予測係数に変換可能な係数の複数個の候補を得る線形変換部２２５を備えていると言える。この場合、符号化部２２４は、線形予測分析部２２１が得た線形予測係数に変換可能な係数と、適合部２２Ａが得た第一線形変換後の線形予測係数に変換可能な係数の複数個の候補と、を用いて、線形予測分析部２２１が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得ていると言える。
【０１８２】
第一実施形態の（２）第２の場合の適合部２２Ａは、線形予測分析部２２１が得た線形予測係数に変換可能な係数に対して、η1に応じた第二線形変換を行い、第二線形変換後の線形予測係数に変換可能な係数を得る線形変換部２２５を備えていると言える。この場合、符号化部２２４は、適合部２２Ａが得た第二線形変換後の線形予測係数に変換可能な係数と、符号帳に格納された線形予測係数に変換可能な係数の複数個の候補と、を用いて、線形予測分析部２２１が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得ていると言える。
【０１８３】
第一実施形態の（３）第３の場合の適合部２２Ａは、符号帳記憶部２２２には、η2に対応する符号帳が記憶されているとして、符号帳記憶部２２２に記憶された線形予測係数に変換可能な係数の複数個の候補に対して、η3に応じた第一線形変換を行い、第一線形変換後の線形予測係数に変換可能な係数の複数個の候補を得、線形予測分析部２２１が得た線形予測係数に変換可能な係数に対して、η3に応じた第二線形変換を行い、第二線形変換後の線形予測係数に変換可能な係数を得ていると言える。この場合、符号化部２２４は、適合部２２Ａが得た第二線形変換後の線形予測係数に変換可能な係数と、適合部２２Ａが得た第一線形変換後の線形予測係数に変換可能な係数の複数個の候補と、を用いて、線形予測分析部が得た線形予測係数に変換可能な係数に対応する線形予測係数符号を得ていると言える。
【０１８４】
適合部２２Ａは、例えば図２５に示す符号帳選択部２２３及び第二線形変換部２２５２により、符号帳の適合を行ってもよい。例えば、パラメータη2は所定のパラメータηであるとして、符号帳選択部２２３は、符号帳記憶部２２２に記憶された複数の符号帳の中からパラメータη2に応じて符号帳を選択する。そして、第二線形変換部２２５２は、線形予測分析部２２１で得られた線形予測係数に変換可能な係数に対する、η2に応じた第二線形変換を行う。この場合、符号化部２２４は、第二線形変換後の線形予測係数に変換可能な係数について、選択された符号帳を用いて符号化して線形予測係数符号を得る。
【０１８５】
適合部２２Ａは、例えば図２６に示す符号帳選択部２２３及び第一線形変換部２２５１により、符号帳の適合を行ってもよい。例えば、パラメータη2は所定のパラメータηであるとして、符号帳選択部２２３は、符号帳記憶部２２２に記憶された複数の符号帳の中からパラメータη2に応じて符号帳を選択する。そして、第一線形変換部２２５１は、選択された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補に対する、η1に応じた第一線形変換を行う。この場合、符号化部２２４は、線形予測分析部２２１で得られた線形予測係数に変換可能な係数について、第一線形変換後の線形予測係数に変換可能な係数の候補を用いて符号化して線形予測係数符号を得る。
【０１８６】
適合部２２Ａは、例えば図２７に示す符号帳選択部２２３、第一線形変換部２２５１及び第二変換部２２５２により、符号帳の適合を行ってもよい。例えば、パラメータη2,η3は所定のパラメータηであるとして、符号帳選択部２２３は、符号帳記憶部２２２に記憶された複数の符号帳の中からパラメータη3に応じて符号帳を選択する。そして、第一線形変換部２２５１は、選択された符号帳に格納された線形予測係数に変換可能な係数の複数個の候補に対する、η2に応じた第一線形変換を行う。そして、第二線形変換部２２５２は、線形予測分析部２２１で得られた線形予測係数に変換可能な係数に対する、η2に応じた第二線形変換を行う。この場合、符号化部２２４は、第二線形変換後の線形予測係数に変換可能な係数について、第一線形変換後の線形予測係数に変換可能な係数の候補を用いて符号化して線形予測係数符号を得る。
【０１８７】
図６、図２３及び図２８に一点鎖線で示すように、適合部３１Ａが符号帳選択部３１２及び線形変換部３１４の少なくとも一方と、復号部３１３とから構成されているとすると、適合部３１Ａは、η1を正の数として、入力されたη1に基づいて、符号帳記憶部３１１に記憶された符号帳と、符号帳に格納された複数個の線形予測係数に変換可能な係数の候補のうち、入力された線形予測係数符号に対応する線形予測係数に変換可能な係数の候補との少なくとも一方を適合させていると言える。
【０１８８】
適合部３１Ａは、例えば図２８に示す符号帳選択部３１２及び線形変換部３１４の両方において適合の処理を行ってもよい。例えば、η2を正の数として、符号帳選択部３１２は、符号帳記憶部３１１に記憶された複数の符号帳の中からパラメータη２に応じて符号帳を選択する。そして、線形変換部３１４は、復号部３１３で得られた線形予測係数に変換可能な係数に対して、所定の正の数であるη1に応じた線形変換をして線形予測係数に変換可能な係数を得る。
【０１８９】
［符号化装置、復号装置及びこれらの方法］
以下、線形予測符号化装置、線形予測復号装置及びこれらの方法を用いた符号化装置、復号装置及びこれらの方法の例について説明する。
【０１９０】
［符号化装置、復号装置及びこれらの方法の第一実施形態］
（符号化）
第一実施形態の符号化装置の構成例を図８に示す。第一実施形態の符号化装置は、図８に示すように、周波数領域変換部２１と、線形予測分析部２２と、非平滑化振幅スペクトル包絡系列生成部２３と、平滑化振幅スペクトル包絡系列生成部２４と、包絡正規化部２５と、符号化部２６と、パラメータ決定部２７とを例えば備えている。この符号化装置により実現される第一実施形態の符号化方法の各処理の例を図９に示す。
【０１９１】
以下、図８の各部について説明する。
【０１９２】
＜パラメータ決定部２７＞
第一実施形態では、所定の時間区間ごとに複数のパラメータηの何れかがパラメータ決定部２７により選択可能とされている。
【０１９３】
パラメータ決定部２７には、複数のパラメータηがパラメータηの候補として記憶されているとする。パラメータ決定部２７は、複数のパラメータの中の１つのパラメータηを順次読み出し、線形予測分析部２２、非平滑化振幅スペクトル包絡系列生成部２３及び復号化部２６に出力する（ステップＡ０）。
【０１９４】
周波数領域変換部２１、線形予測分析部２２、非平滑化振幅スペクトル包絡系列生成部２３、平滑化振幅スペクトル包絡系列生成部２４、包絡正規化部２５及び符号化部２６は、パラメータ決定部２７が順次読み出した各パラメータηに基づいて、例えば以下に説明するステップＡ１からステップＡ６の処理を行い同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対して符号を生成する。一般に、パラメータηを所与として、同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対して２個以上の符号が得られる場合がある。この場合、同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対する符号は、これらの得られた２個以上の符号をまとめたものである。この例では、符号は、線形予測係数符号と、利得符号と、整数信号符号とを合わせたものである。これにより、同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対する各パラメータηごとの符号が得られる。
【０１９５】
ステップＡ６の処理の後に、パラメータ決定部２７は、同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対して各パラメータηごとに得られた符号の中から１つの符号を選択し、選択された符号に対応するパラメータηを決定する（ステップＡ７）。この決定されたパラメータηが、その同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対するパラメータηとなる。そして、パラメータ決定部２７は、選択された符号及び決定されたパラメータηを表す符号を復号装置に出力する。パラメータ決定部２７によるステップＡ７の処理の詳細については後述する。
【０１９６】
以下では、パラメータ決定部２７により１つのパラメータη1が読み出されており、この読み出された１つのパラメータη1について処理が行われるとする。
【０１９７】
＜周波数領域変換部２１＞
周波数領域変換部２１には、時間領域の時系列信号である音信号が入力される。音信号の例は、音声ディジタル信号又は音響ディジタル信号である。
【０１９８】
周波数領域変換部２１は、所定の時間長のフレーム単位で、入力された時間領域の音信号を周波数領域のN点のMDCT係数列X(0),X(1),…,X(N-1)に変換する（ステップＡ１）。Nは正の整数である。
【０１９９】
得られたMDCT係数列X(0),X(1),…,X(N-1)は、線形予測分析部２２と包絡正規化部２５に出力される。
【０２００】
特に断りがない限り、以降の処理はフレーム単位で行われるものとする。
【０２０１】
このようにして、周波数領域変換部２１は、音信号に対応する、例えばMDCT係数列である周波数領域サンプル列を求める。
【０２０２】
＜線形予測分析部２２＞
線形予測分析部２２には、周波数領域変換部２１が得たMDCT係数列X(0),X(1),…,X(N-1)が入力される。
【０２０３】
線形予測分析部２２は、［線形予測符号化装置、線形予測復号装置及びこれらの方法］で説明した図１から図３、図２１の何れかの線形予測符号化装置である。［符号化装置、復号装置及びこれらの方法］及び図８では、［線形予測符号化装置、線形予測復号装置及びこれらの方法］で説明した図１から図３、図２１の何れかの線形予測符号化装置のことを「線形予測分析部２２」と表記する。なお、線形予測分析部２２は、図２５から図２７の何れかの線形予測符号化装置であってもよい。
【０２０４】
線形予測分析部２２は、［線形予測符号化装置、線形予測復号装置及びこれらの方法］で説明した処理と同様の処理により、例えばMDCT係数列である周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得て、得られた線形予測係数に変換可能な係数を符号化して線形予測係数符号を得る。
【０２０５】
得られた線形予測係数符号は、パラメータ決定部２７及び復号装置に出力される。
【０２０６】
また、線形予測符号化装置の線形変換部２２５が（１）第１の場合には、符号化部２２４で得られた線形予測係数符号に対応する、パラメータη1に対応する線形予測係数に変換可能な係数が量子化線形予測係数^β1,^β2,…,^βpとして、非平滑化スペクトル包絡系列生成部２３と平滑化振幅スペクトル包絡系列生成部２４に出力される。
【０２０７】
線形予測符号化装置の線形変換部２２５が（２）第２の場合には、符号化部２２４で得られた線形予測係数符号に対応する、パラメータη2に対応する線形予測係数に変換可能な係数が、図２に破線で示す逆線形変換部２２６に入力される。逆線形変換部２２６は、線形予測係数符号に対応する、パラメータη2に対応する線形予測係数に変換可能な係数に対して第二線形変換部２２５２が行った第二線形変換の逆の線形変換を行い、パラメータη1に対応する線形予測係数に変換可能な係数とする。このパラメータη1に対応する線形予測係数に変換可能な係数が、量子化線形予測係数^β1,^β2,…,^βpとして、非平滑化スペクトル包絡系列生成部２３と平滑化振幅スペクトル包絡系列生成部２４に出力される。なお、パラメータη1の値とパラメータη2の値とが同一である場合には、逆線形変換部２２６は、線形変換をしなくてもよい。
【０２０８】
線形予測符号化装置の線形変換部２２５が（３）第３の場合には、符号化部２２４で得られた線形予測係数符号に対応する、パラメータη3に対応する線形予測係数に変換可能な係数が、図３に破線で示す逆線形変換部２２６に入力される。逆線形変換部２２６は、線形予測係数符号に対応する、パラメータη3に対応する線形予測係数に変換可能な係数に対して第二線形変換部２２５２が行った第二線形変換の逆の線形変換を行い、パラメータη1に対応する線形予測係数に変換可能な係数とする。このパラメータη1に対応する線形予測係数に変換可能な係数が、量子化線形予測係数^β1,^β2,…,^βpとして、非平滑化スペクトル包絡系列生成部２３と平滑化振幅スペクトル包絡系列生成部２４に出力される。なお、パラメータη1の値とパラメータη3の値とが同一である場合には、逆線形変換部２２６は、線形変換をしなくてもよい。
【０２０９】
なお、線形予測分析処理の過程で予測残差のエネルギーσ2が算出される。この場合、算出された予測残差のエネルギーσ2は、符号化部２６の分散パラメータ決定部２６８に出力される。
【０２１０】
＜非平滑化振幅スペクトル包絡系列生成部２３＞
非平滑化振幅スペクトル包絡系列生成部２３には、線形予測分析部２２が生成した量子化線形予測係数^β1,^β2,…,^βpが入力される。
【０２１１】
非平滑化振幅スペクトル包絡系列生成部２３は、量子化線形予測係数^β1,^β2,…,^βpに対応する振幅スペクトル包絡の系列である非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)を生成する（ステップＡ３）。
【０２１２】
生成された非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)は、符号化部２６に出力される。
【０２１３】
非平滑化振幅スペクトル包絡系列生成部２３は、量子化線形予測係数^β1,^β2,…,^βpを用いて、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)として、式(A2)により定義される非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)を生成する。
【０２１４】
【数４】
2016167215.tif 000006
【０２１５】
このようにして、非平滑化振幅スペクトル包絡系列生成部２３は、線形予測分析部２２により生成された線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列を得ることによりスペクトル包絡の推定を行う。ここで、ｃを任意の数として、複数の値から構成される系列をｃ乗した系列とは、複数の値のそれぞれをｃ乗した値から構成される系列のことである。例えば、振幅スペクトル包絡の系列を１／η1乗した系列とは、振幅スペクトル包絡の各係数を１／η1乗した値から構成される系列のことである。
【０２１６】
非平滑化振幅スペクトル包絡系列生成部２３による１／η1乗の処理は、線形予測分析部２２で行われた周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した処理に起因するものである。すなわち、非平滑化振幅スペクトル包絡系列生成部２３による１／η1乗の処理は、線形予測分析部２２で行われた周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した処理によりη1乗された値を元の値に戻すために行われる。
【０２１７】
＜平滑化振幅スペクトル包絡系列生成部２４＞
平滑化振幅スペクトル包絡系列生成部２４には、線形予測分析部２２が生成した量子化線形予測係数^β1,^β2,…,^βpが入力される。
【０２１８】
平滑化振幅スペクトル包絡系列生成部２４は、量子化線形予測係数^β1,^β2,…,^βpに対応する振幅スペクトル包絡の系列の振幅の凸凹を鈍らせた系列である平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)を生成する（ステップＡ４）。
【０２１９】
生成された平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)は、包絡正規化部２５及び符号化部２６に出力される。
【０２２０】
平滑化振幅スペクトル包絡系列生成部２４は、量子化線形予測係数^β1,^β2,…,^βpと補正係数γを用いて、平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)として、式(A3)により定義される平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)を生成する。
【０２２１】
【数５】
2016167215.tif 000007
【０２２２】
ここで、補正係数γは予め定められた１未満の定数であり非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)の振幅の凹凸を鈍らせる係数、言い換えれば非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)を平滑化する係数である。
【０２２３】
＜包絡正規化部２５＞
包絡正規化部２５には、周波数領域変換部２１が得たMDCT係数列X(0),X(1),…,X(N-1)及び平滑化振幅スペクトル包絡生成部２４が生成した平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)が入力される。
【０２２４】
包絡正規化部２５は、MDCT係数列X(0),X(1),…,X(N-1)の各係数を、対応する平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)の各値で正規化することにより、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)を生成する（ステップＡ５）。
【０２２５】
生成された正規化MDCT係数列は、符号化部２６に出力される。
【０２２６】
包絡正規化部２５は、例えば、k=0,1,…,N-1として、MDCT係数列X(0),X(1),…,X(N-1)の各係数X(k)を平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)で除算することにより、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)の各係数XN(k)を生成する。すなわち、k=0,1,…,N-1として、XN(k)=X(k)/^Hγ(k)である。
【０２２７】
＜符号化部２６＞
符号化部２６には、包絡正規化部２５が生成した正規化MDCT係数列XN(0),XN(1),…,XN(N-1)、非平滑化振幅スペクトル包絡生成部２３が生成した非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)、平滑化振幅スペクトル包絡生成部２４が生成した平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)及び線形予測分析部２２が算出した予測残差のエネルギーσ2が入力される。
【０２２８】
符号化部２６は、図１２に示すステップＡ６１からステップＡ６５の処理を例えば行うことにより符号化を行う（ステップＡ６）。
【０２２９】
符号化部２６は、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)に対応するグローバルゲインgを求め（ステップＡ６１）、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)の各係数をグローバルゲインgで割り算した結果を量子化した整数値による系列である量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)を求め（ステップＡ６２）、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の各係数に対応する分散パラメータφ(0),φ(1),…,φ(N-1)をグローバルゲインgと非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)と平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)と平均残差のエネルギーσ2とから式(A1)により求め（ステップＡ６３）、分散パラメータφ(0),φ(1),…,φ(N-1)を用いて量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)を算術符号化して整数信号符号を得（ステップＡ６４）、グローバルゲインgに対応する利得符号を得る（ステップＡ６５）。
【０２３０】
【数６】
2016167215.tif 000008
【０２３１】
ここで、上記の式(A1)における正規化振幅スペクトル包絡系列^HN(0),^HN(1),…,^HNは、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)の各値を、対応する平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)の各値で除算したもの、すなわち、以下の式(A8)により求まるものである。
【０２３２】
【数７】
2016167215.tif 000009
【０２３３】
生成された整数信号符号と利得符号は正規化MDCT係数列に対応する符号として、パラメータ決定部２７に出力される。
【０２３４】
符号化部２６は、ステップＡ６１からステップＡ６５により、整数信号符号のビット数が、予め配分されたビット数である配分ビット数B以下、かつ、なるべく大きな値となるようなグローバルゲインgを決定し、決定されたグローバルゲインgに対応する利得符号と、この決定されたグローバルゲインgに対応する整数信号符号とを生成する機能を実現している。
【０２３５】
符号化部２６が行うステップＡ６１からステップＡ６５のうち、の特徴的な処理が含まれるのはステップＡ６３であり、グローバルゲインgと量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)のそれぞれを符号化することにより正規化MDCT係数列に対応する符号を得る符号化処理自体には、非特許文献１に記載された技術を含む様々な公知技術が存在する。以下では符号化部２６が行う符号化処理の具体例を２つ説明する。
【０２３６】
［符号化部２６が行う符号化処理の具体例１］
符号化部２６が行う符号化処理の具体例１として、ループ処理を含まない例について説明する。
【０２３７】
具体例１の符号化部２６の構成例を図１０に示す。具体例１の符号化部２６は、図１０に示すように、利得取得部２６１と、量子化部２６２と、分散パラメータ決定部２６８と、算術符号化部２６９と、利得符号化部２６５とを例えば備えている。以下、図１０の各部について説明する。
【０２３８】
＜利得取得部２６１＞
利得取得部２６１には、包絡正規化部２５が生成した正規化MDCT係数列XN(0),XN(1),…,XN(N-1)が入力される。
【０２３９】
利得取得部２６１は、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)から、整数信号符号のビット数が、予め配分されたビット数である配分ビット数B以下、かつ、なるべく大きな値となるようなグローバルゲインgを決定して出力する（ステップＳ２６１）。利得取得部２６１は、例えば、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)のエネルギーの合計の平方根と配分ビット数Bと負の相関のある定数との乗算値をグローバルゲインgとして得て出力する。または、利得取得部２６１は、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)のエネルギーの合計と、配分ビット数Bと、グローバルゲインgと、の関係を予めテーブル化しておき、そのテーブルを参照することによりグローバルゲインgを得て出力してもよい。
【０２４０】
このようにして、利得取得部２６１は、例えば正規化MDCT係数列である正規化周波数領域サンプル列の全サンプルを除算するための利得を得る。
【０２４１】
得られたグローバルゲインgは、量子化部２６２及び分散パラメータ決定部２６８に出力される。
【０２４２】
＜量子化部２６２＞
量子化部２６２には、包絡正規化部２５が生成した正規化MDCT係数列XN(0),XN(1),…,XN(N-1)及び利得取得部２６１が得たグローバルゲインgが入力される。
【０２４３】
量子化部２６２は、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)の各係数をグローバルゲインgで割り算した結果の整数部分による系列である量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)を得て出力する（ステップＳ２６２）。
【０２４４】
このようにして、量子化部２６２は、例えば正規化MDCT係数列である正規化周波数領域サンプル列の各サンプルを、利得で除算するとともに量子化して量子化正規化済係数系列を求める。
【０２４５】
得られた量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)は、算術符号化部２６９に出力される。
【０２４６】
＜分散パラメータ決定部２６８＞
分散パラメータ決定部２６８には、パラメータ決定部２７が読み出したパラメータη1、利得取得部２６１が得たグローバルゲインg、非平滑化振幅スペクトル包絡生成部２３が生成した非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)、平滑化振幅スペクトル包絡生成部２４が生成した平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)及び線形予測分析部２２が得た予測残差のエネルギーσ2が入力される。
【０２４７】
分散パラメータ決定部２６８は、グローバルゲインgと、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)と、平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)と、予測残差のエネルギーσ2とから、上記の式(A1),式(A8)により分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータを得て出力する（ステップＳ２６８）。
【０２４８】
得られた分散パラメータ系列φ(0),φ(1),…,φ(N-1)は、算術符号化部２６９に出力される。
【０２４９】
＜算術符号化部２６９＞
算術符号化部２６９には、パラメータ決定部２７が読み出したパラメータη1、量子化部２６２が得た量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)及び分散パラメータ決定部２６８が得た分散パラメータ系列φ(0),φ(1),…,φ(N-1)が入力される。
【０２５０】
算術符号化部２６９は、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の各係数に対応する分散パラメータとして分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータを用いて、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)を算術符号化して整数信号符号を得て出力する（ステップＳ２６９）。
【０２５１】
算術符号化部２６９は、算術符号化の際に、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の各係数が一般化ガウス分布fGG(X|φ(k),η1)に従うときに最適になるような算術符号を構成し、この構成に基づく算術符号により符号化を行う。この結果、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の各係数へのビット割り当ての期待値が分散パラメータ系列φ(0),φ(1),…,φ(N-1)で決定されることになる。
【０２５２】
得られた整数信号符号は、パラメータ決定部２７に出力される。
【０２５３】
量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の中の複数の係数に跨って算術符号化が行われてもよい。この場合、分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータは、式(A1),式(A8)からわかるように、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)に基づいているため、算術符号化部２６９は、推定されたスペクトル包絡（非平滑化振幅スペクトル包絡）を基に実質的にビット割り当てが変わる符号化を行っていると言える。
【０２５４】
＜利得符号化部２６５＞
利得符号化部２６５には、利得取得部２６１が得たグローバルゲインｇが入力される。
【０２５５】
利得符号化部２６５は、グローバルゲインｇを符号化して利得符号を得て出力する（ステップＳ２６５）。
【０２５６】
生成された整数信号符号と利得符号は正規化MDCT係数列に対応する符号として、パラメータ決定部２７に出力される。
【０２５７】
本具体例１のステップＳ２６１，Ｓ２６２，Ｓ２６８，Ｓ２６９，Ｓ２６５がそれぞれ上記のステップＡ６１，Ａ６２，Ａ６３，Ａ６４，Ａ６５に対応する。
【０２５８】
［符号化部２６が行う符号化処理の具体例２］
符号化部２６が行う符号化処理の具体例２として、ループ処理を含む例について説明する。
【０２５９】
具体例２の符号化部２６の構成例を図１１に示す。具体例２の符号化部２６は、図１１に示すように、利得取得部２６１と、量子化部２６２と、分散パラメータ決定部２６８と、算術符号化部２６９と、利得符号化部２６５と、判定部２６６と、利得更新部２６７とを例えば備えている。以下、図１１の各部について説明する。
【０２６０】
＜利得取得部２６１＞
利得部２６１には、包絡正規化部２５が生成した正規化MDCT係数列XN(0),XN(1),…,XN(N-1)が入力される。
【０２６１】
利得取得部２６１は、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)から、整数信号符号のビット数が、予め配分されたビット数である配分ビット数B以下、かつ、なるべく大きな値となるようなグローバルゲインgを決定して出力する（ステップＳ２６１）。利得取得部２６１は、例えば、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)のエネルギーの合計の平方根と配分ビット数Bと負の相関のある定数との乗算値をグローバルゲインgとして得て出力する。
【０２６２】
得られたグローバルゲインgは、量子化部２６２及び分散パラメータ決定部２６８に出力される。
【０２６３】
利得取得部２６１が得たグローバルゲインgは、量子化部２６２及び分散パラメータ決定部２６８で用いられるグローバルゲインの初期値となる。
【０２６４】
＜量子化部２６２＞
量子化部２６２には、包絡正規化部２５が生成した正規化MDCT係数列XN(0),XN(1),…,XN(N-1)及び利得取得部２６１又は利得更新部２６７が得たグローバルゲインgが入力される。
【０２６５】
量子化部２６２は、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)の各係数をグローバルゲインgで割り算した結果の整数部分による系列である量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)を得て出力する（ステップＳ２６２）。
【０２６６】
ここで、量子化部２６２が初回に実行される際に用いられるグローバルゲインgは、利得取得部２６１が得たグローバルゲインg、すなわちグローバルゲインの初期値である。また、量子化部２６２が２回目以降に実行される際に用いられるグローバルゲインgは、利得更新部２６７が得たグローバルゲインg、すなわちグローバルゲインの更新値である。
【０２６７】
得られた量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)は、算術符号化部２６９に出力される。
【０２６８】
＜分散パラメータ決定部２６８＞
分散パラメータ決定部２６８には、パラメータ決定部２７が読み出したパラメータη1、利得取得部２６１又は利得更新部２６７が得たグローバルゲインg、非平滑化振幅スペクトル包絡生成部２３が生成した非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)、平滑化振幅スペクトル包絡生成部２４が生成した平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)及び線形予測分析部２２が得た予測残差のエネルギーσ2が入力される。
【０２６９】
分散パラメータ決定部２６８は、グローバルゲインgと、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)と、平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)と、予測残差のエネルギーσ2とから、上記の式(A1),式(A8)により分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータを得て出力する（ステップＳ２６８）。
【０２７０】
ここで、分散パラメータ決定部２６８が初回に実行される際に用いられるグローバルゲインgは、利得取得部２６１が得たグローバルゲインg、すなわちグローバルゲインの初期値である。また、分散パラメータ決定部２６８が２回目以降に実行される際に用いられるグローバルゲインgは、利得更新部２６７が得たグローバルゲインg、すなわちグローバルゲインの更新値である。
【０２７１】
得られた分散パラメータ系列φ(0),φ(1),…,φ(N-1)は、算術符号化部２６９に出力される。
【０２７２】
＜算術符号化部２６９＞
算術符号化部２６９には、パラメータ決定部２７が読み出したパラメータη1、量子化部２６２が得た量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)及び分散パラメータ決定部２６８が得た分散パラメータ系列φ(0),φ(1),…,φ(N-1)が入力される。
【０２７３】
算術符号化部２６９は、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の各係数に対応する分散パラメータとして分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータを用いて、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)を算術符号化して、整数信号符号と整数信号符号のビット数である消費ビット数Ｃとを得て出力する（ステップＳ２６９）。
【０２７４】
算術符号化部２６９は、算術符号化の際に、量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の各係数が一般化ガウス分布fGG(X|φ(k),η1)に従うときに最適になるようなビット割り当てを算術符号により行い、行われたビット割り当てに基づく算術符号により符号化を行う。
【０２７５】
得られた整数信号符号及び消費ビット数Ｃは、判定部２６６に出力される。
【０２７６】
量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の中の複数の係数に跨って算術符号化が行われてもよい。この場合、分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータは、式(A1),式(A8)からわかるように、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)に基づいているため、算術符号化部２６９は、推定されたスペクトル包絡（非平滑化振幅スペクトル包絡）を基に実質的にビット割り当てが変わる符号化を行っていると言える。
【０２７７】
＜判定部２６６＞
判定部２６６には、算術符号化部２６９が得た整数信号符号が入力される。
【０２７８】
判定部２６６は、利得の更新回数が予め定めた回数の場合には、整数信号符号を出力するとともに、利得符号化部２６５に対し利得更新部２６７が得たグローバルゲインｇを符号化する指示信号を出力し、利得の更新回数が予め定めた回数未満である場合には、利得更新部２６７に対し、算術符号化部２６４が計測した消費ビット数Ｃを出力する（ステップＳ２６６）。
【０２７９】
＜利得更新部２６７＞
利得更新部２６７には、算術符号化部２６４が計測した消費ビット数Ｃが入力される。
【０２８０】
利得更新部２６７は、消費ビット数Ｃが配分ビット数Ｂより多い場合にはグローバルゲインｇの値を大きな値に更新して出力し、消費ビット数Ｃが配分ビット数Ｂより少ない場合にはグローバルゲインｇの値を小さな値に更新し、更新後のグローバルゲインｇの値を出力する（ステップＳ２６７）。
【０２８１】
利得更新部２６７が得た更新後のグローバルゲインgは、量子化部２６２及び利得符号化部２６５に出力される。
【０２８２】
＜利得符号化部２６５＞
利得符号化部２６５には、判定部２６６からの出力指示及び利得更新部２６７が得たグローバルゲインｇが入力される。
【０２８３】
利得符号化部２６５は、指示信号に従って、グローバルゲインｇを符号化して利得符号を得て出力する（ステップ２６５）。
【０２８４】
判定部２６６が出力した整数信号符号と、利得符号化部２６５が出力した利得符号は、正規化MDCT係数列に対応する符号として、パラメータ決定部２７に出力される。
【０２８５】
すなわち、本具体例２においては、最後に行われたステップＳ２６７が上記のステップＡ６１に対応し、ステップＳ２６２，Ｓ２６３，Ｓ２６４，Ｓ２６５がそれぞれ上記のステップＡ６２，Ａ６３，Ａ６４，Ａ６５に対応する。
【０２８６】
なお、符号化部２６が行う符号化処理の具体例２については、国際公開公報WO2014/054556などに更に詳細に説明されている。
【０２８７】
［符号化部２６の変形例］
符号化部２６は、例えば以下の処理を行うことにより、推定されたスペクトル包絡（非平滑化振幅スペクトル包絡）を基にビット割り当てを変える符号化を行ってもよい。
【０２８８】
符号化部２６は、まず、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)に対応するグローバルゲインgを求め、正規化MDCT係数列XN(0),XN(1),…,XN(N-1)の各係数をグローバルゲインgで割り算した結果を量子化した整数値による系列である量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)を求める。
【０２８９】
この量子化正規化済係数系列XQ(0),XQ(1),…,XQ(N-1)の各係数に対応する量子化ビットは、XQ(k)の分布がある範囲内で一様であると仮定して、その範囲を包絡の推定値から決めることができる。複数のサンプルごとの包絡の推定値を符号化することもできるが、符号化部２６は、例えば以下の式(A9)のように線形予測に基づく正規化振幅スペクトル包絡系列の値^HN(k)を使用してXQ(k)の範囲を決めることができる。
【０２９０】
【数８】
2016167215.tif 000010
【０２９１】
あるkにおけるXQ(k)を量子化するときに、XQ(k)の二乗誤差を最小とするために
【０２９２】
【数９】
2016167215.tif 000011
【０２９３】
の制約のもとに、割り当てるビット数b(k)
【０２９４】
【数１０】
2016167215.tif 000012
【０２９５】
を設定することができる。Bは予め定められた正の整数である。この際にb(k)が整数となるように四捨五入するとか、０より小さくなる場合にはb(k)=0とするなどして、b(k)の再調整の処理を符号化部２６は行ってもよい。
【０２９６】
また、符号化部２６は、サンプルごとの割り当てでなく、複数のサンプルをまとめて配分ビット数を決めて、量子化にもサンプルごとのスカラ量子化でなく、複数のサンプルをまとめたベクトルごとの量子化をすることも可能である。
【０２９７】
サンプルkのXQ(k)の量子化ビット数b(k)が上記で与えられ、サンプルごとに符号化するとすると、XQ(k)は-2b(k)-1から2b(k)-1までの2b(k)種類の整数を取り得る。符号化部２６は、b(k)ビットで各サンプルを符号化して整数信号符号を得る。
【０２９８】
生成された整数信号符号は、復号装置に出力される。例えば、生成されたXQ(k)に対応するb(k)ビットの整数信号符号は、k=0から順次復号装置に出力される。
【０２９９】
もし、XQ(k)が上記の-2b(k)-1から2b(k)-1までの範囲をこえる場合には最大値、または最小値に置き換える。
【０３００】
ｇが小さすぎるとこの置き換えで量子化歪が発生し、ｇが大きすぎると量子化誤差は大きくなり、XQ(k)のとりうる範囲がb(k)に比べて小さすぎて、情報の有効利用ができないことになる。このため、ｇの最適化を行ってもよい。
【０３０１】
符号化部２６は、グローバルゲインgを符号化して利得符号を得て出力する。
【０３０２】
この符号化部２６の変形例のように、符号化部２６は算術符号化以外の符号化を行ってもよい。
【０３０３】
＜パラメータ決定部２７＞
ステップＡ１からステップＡ６の処理により、同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対して各パラメータη1ごとに生成された符号（この例では、線形予測係数符号、利得符号及び整数信号符号）は、パラメータ決定部２７に入力される。
【０３０４】
パラメータ決定部２７は、同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対して各パラメータη1ごとに得られた符号の中から１つの符号を選択し、選択された符号に対応するパラメータηを決定する（ステップＡ７）。この決定されたパラメータηが、その同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対するパラメータηとなる。そして、パラメータ決定部２７は、選択された符号及び決定されたパラメータηを表すパラメータ符号を復号装置に出力する。符号の選択は、符号の符号量及び符号に対応する符号化歪の少なくとも一方に基づいて行われる。例えば、符号量が最も小さい符号又は符号化歪が最も小さい符号が選択される。
【０３０５】
ここで、符号化歪みとは、入力信号から得られる周波数領域サンプル列と、生成された符号をローカルデコードすることにより得られる周波数領域サンプル列との誤差のことである。符号化装置は、符号化歪みを計算するための符号化歪計算部を備えていてもよい。この符号化歪計算部は、以下に述べる復号装置と同様の処理を行う復号部を備え、この復号部が生成された符号をローカルデコードする。その後、符号化歪計算部は、入力信号から得られる周波数領域サンプル列と、ローカルデコードすることにより得られた周波数領域サンプル列との誤差を計算し、符号化歪とする。
【０３０６】
（復号）
符号化装置に対応する復号装置の構成例を図１３に示す。第一実施形態の復号装置は、図１３に示すように、線形予測係数復号部３１と、非平滑化振幅スペクトル包絡系列生成部３２と、平滑化振幅スペクトル包絡系列生成部３３と、復号部３４と、包絡逆正規化部３５と、時間領域変換部３６と、パラメータ復号部３７とを例えば備えている。この復号装置により実現される第一実施形態の復号方法の各処理の例を図１４に示す。
【０３０７】
復号装置には、符号化装置が出力した、パラメータ符号、正規化MDCT係数列に対応する符号及び線形予測係数符号が少なくとも入力される。
【０３０８】
以下、図１３の各部について説明する。
【０３０９】
＜パラメータ復号部３７＞
パラメータ復号部３７には、符号化装置が出力したパラメータ符号が入力される。
【０３１０】
パラメータ復号部３７は、パラメータ符号を復号することにより復号パラメータηを求める。求まった復号パラメータηは、線形予測係数復号部３１、非平滑化振幅スペクトル包絡系列生成部３２、平滑化振幅スペクトル包絡系列生成部３３及び復号部３４に出力される。パラメータ復号部３７には、複数の復号パラメータηが候補として記憶されいる。パラメータ復号部３７は、パラメータ符号に対応する復号パラメータηの候補を復号パラメータηとして求める。パラメータ復号部３７に記憶されている複数の復号パラメータηは、符号化装置のパラメータ決定部２７に記憶された複数のパラメータηと同じである。
【０３１１】
＜線形予測係数復号部３１＞
線形予測係数復号部３１には、符号化装置が出力した線形予測係数符号及びパラメータ復号部３７により得られた復号パラメータηが入力される。
【０３１２】
線形予測係数復号部３１は、［線形予測符号化装置、線形予測復号装置及びこれらの方法］で説明した図６、図２１を用いて上記説明した線形予測復号装置である。［符号化装置、復号装置及びこれらの方法］及び図１３では、［線形予測符号化装置、線形予測復号装置及びこれらの方法］で説明した図６、図２１の線形予測符号化装置のことを「線形予測係数復号部３１」と表記する。なお、線形予測係数復号部３１は、図２８の線形予測復号装置であってもよい。
【０３１３】
線形予測係数復号部３１は、復号パラメータηをパラメータη1とする［線形予測符号化装置、線形予測復号装置及びこれらの方法］で説明した処理と同様の処理により、入力された線形予測係数符号を復号することにより、復号された線形予測係数に変換可能な係数である復号線形予測係数^β1,^β2,…, ^βpを得る（ステップＢ１）。
【０３１４】
得られた復号線形予測係数^β1,^β2,…, ^βpは、非平滑化振幅スペクトル包絡系列生成部３２及び非平滑化振幅スペクトル包絡系列生成部３３に出力される。
【０３１５】
＜非平滑化振幅スペクトル包絡系列生成部３２＞
非平滑化振幅スペクトル包絡系列生成部３２には、パラメータ復号部３７が求めた復号パラメータη及び線形予測係数復号部３１が得た復号線形予測係数^β1,^β2,…,^βpが入力される。
【０３１６】
非平滑化振幅スペクトル包絡系列生成部３２は、復号線形予測係数^β1,^β2,…,^βpに対応する振幅スペクトル包絡の系列である非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)を上記の式(A2)により生成する（ステップＢ２）。
【０３１７】
生成された非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)は、復号部３４に出力される。
【０３１８】
このようにして、非平滑化振幅スペクトル包絡系列生成部３２は、線形予測係数復号部３１により生成された線形予測係数に変換可能な係数に対応するに対応する振幅スペクトル包絡の系列を１／η乗した系列である非平滑化スペクトル包絡系列を得る。
【０３１９】
＜平滑化振幅スペクトル包絡系列生成部３３＞
平滑化振幅スペクトル包絡系列生成部３３には、パラメータ復号部３７が求めた復号パラメータη及び線形予測係数復号部３１が得た復号線形予測係数^β1,^β2,…,^βpが入力される。
【０３２０】
平滑化振幅スペクトル包絡系列生成部３３は、復号線形予測係数^β1,^β2,…,^βpに対応する振幅スペクトル包絡の系列の振幅の凹凸を鈍らせた系列である平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)を上記の式A(3)により生成する（ステップＢ３）。
【０３２１】
生成された平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)は、復号部３４及び包絡逆正規化部３５に出力される。
【０３２２】
＜復号部３４＞
復号部３４には、パラメータ復号部３７が求めた復号パラメータη、符号化装置が出力した正規化MDCT係数列に対応する符号、非平滑化振幅スペクトル包絡生成部３２が生成した非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)及び平滑化振幅スペクトル包絡生成部３３が生成した平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)が入力される。
【０３２３】
復号部３４は、分散パラメータ決定部３４２を備えている。
【０３２４】
復号部３４は、図１５に示すステップＢ４１からステップＢ４４の処理を例えば行うことにより復号を行う（ステップＢ４）。すなわち、復号部３４は、フレームごとに、入力された正規化MDCT係数列に対応する符号に含まれる利得符号を復号してグローバルゲインｇを得る（ステップＢ４１）。復号部３４の分散パラメータ決定部３４２は、グローバルゲインgと非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)と平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)とパラメータηとから上記の式(A1)により分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータを求める（ステップＢ４２）。復号部３４は、正規化MDCT係数列に対応する符号に含まれる整数信号符号を分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータに対応する算術復号の構成に従い、算術復号して復号正規化済係数系列^XQ(0),^XQ(1),…,^XQ(N-1)を得（ステップＢ４３）、復号正規化済係数系列^XQ(0),^XQ(1),…,^XQ(N-1)の各係数にグローバルゲインｇを乗算して復号正規化MDCT係数列^XN(0),^XN(1),…,^XN(N-1)を生成する（ステップＢ４４）。このように、復号部３４は、非平滑化スペクトル包絡系列に基づいて実質的に変わるビット割り当てに従って、入力された整数信号符号の復号を行ってもよい。
【０３２５】
なお、［符号化部２６の変形例］に記載された処理により符号化が行われた場合には、復号部３４は例えば以下の処理を行う。復号部３４は、フレームごとに、入力された正規化MDCT係数列に対応する符号に含まれる利得符号を復号してグローバルゲインｇを得る。復号部３４の分散パラメータ決定部３４２は、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)と平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)とから上記の式(A9)により分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータを求める。復号部３４は、分散パラメータ系列φ(0),φ(1),…,φ(N-1)の各分散パラメータφ(k)に基づいて式(A10)によりb(k)を求めることができ、XQ(k)の値をそのビット数b(k)で順次復号して、復号正規化済係数系列^XQ(0),^XQ(1),…,^XQ(N-1)を得て、復号正規化済係数系列^XQ(0),^XQ(1),…,^XQ(N-1)の各係数にグローバルゲインｇを乗算して復号正規化MDCT係数列^XN(0),^XN(1),…,^XN(N-1)を生成する。このように、復号部３４は、非平滑化スペクトル包絡系列に基づいて変わるビット割り当てに従って、入力された整数信号符号の復号を行ってもよい。
【０３２６】
生成された復号正規化MDCT係数列^XN(0),^XN(1),…,^XN(N-1)は、包絡逆正規化部３５に出力される。
【０３２７】
＜包絡逆正規化部３５＞
包絡逆正規化部３５には、平滑化振幅スペクトル包絡生成部３３が生成した平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)及び復号部３４が生成した復号正規化MDCT係数列^XN(0),^XN(1),…,^XN(N-1)が入力される。
【０３２８】
包絡逆正規化部３５は、平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)を用いて、復号正規化MDCT係数列^XN(0),^XN(1),…,^XN(N-1)を逆正規化することにより、復号MDCT係数列^X(0),^X(1),…,^X(N-1)を生成する（ステップＢ５）。
【０３２９】
生成された復号MDCT係数列^X(0),^X(1),…,^X(N-1)は、時間領域変換部３６に出力される。
【０３３０】
例えば、包絡逆正規化部３５は、k=0,1,…,N-1として、復号正規化MDCT係数列^XN(0),^XN(1),…,^XN(N-1)の各係数^XN(k)に、平滑化振幅スペクトル包絡系列^Hγ(0),^Hγ(1),…,^Hγ(N-1)の各包絡値^Hγ(k)を乗じることにより復号MDCT係数列^X(0),^X(1),…,^X(N-1)を生成する。すなわち、k=0,1,…,N-1として、^X(k)=^XN(k)×^Hγ(k)である。
【０３３１】
＜時間領域変換部３６＞
時間領域変換部３６には、包絡逆正規化部３５が生成した復号MDCT係数列^X(0),^X(1),…,^X(N-1)が入力される。
【０３３２】
時間領域変換部３６は、フレームごとに、包絡逆正規化部３５が得た復号MDCT係数列^X(0),^X(1),…,^X(N-1)を時間領域に変換してフレーム単位の音信号（復号音信号）を得る（ステップＢ６）。
【０３３３】
このようにして、復号装置は、周波数領域での復号により時系列信号を得る。
【０３３４】
［符号化装置、復号装置及びこれらの方法の第二実施形態］
第一実施形態の符号化装置及び方法は、複数のパラメータηのそれぞれについて符号化を行い符号を生成し、パラメータηごとに生成された符号の中から最適な符号を選択し、選択された符号及び選択された符号に対応するパラメータ符号を出力するものであった。
【０３３５】
これに対して、第二実施形態の符号化装置及び方法は、まずパラメータ決定部２７がパラメータηを決定し、決定されたパラメータηに基づいて符号化を行い符号を生成し出力するものである。第二実施形態では、所定の時間区間ごとにパラメータηがパラメータ決定部２７により可変とされている。ここで、所定の時間区間ごとにパラメータηが可変とは、所定の時間区間が変わればパラメータηも変わり得ることを意味し、同一の時間区間ではパラメータηの値は変わらないとする。
【０３３６】
以下、第一実施形態と異なる部分を中心に説明する。第一実施形態と同様の部分については重複説明を省略する。
【０３３７】
（符号化）
第二実施形態の符号化装置の構成例を図１６に示す。符号化装置は、図１６に示すように、周波数領域変換部２１と、線形予測分析部２２と、非平滑化振幅スペクトル包絡系列生成部２３と、平滑化振幅スペクトル包絡系列生成部２４と、包絡正規化部２５と、符号化部２６と、パラメータ決定部２７’とを例えば備えている。この符号化装置により実現される符号化方法の各処理の例を図１７に示す。
【０３３８】
以下、図１６の各部について説明する。
【０３３９】
＜パラメータ決定部２７’＞
パラメータ決定部２７’には、時系列信号である時間領域の音信号が入力される。音信号の例は、音声ディジタル信号又は音響ディジタル信号である。
【０３４０】
パラメータ決定部２７’は、入力された時系列信号に基づいて、後述する処理により、パラメータηを決定する（ステップＡ７’）。以下、パラメータ決定部２７’により決定されたパラメータηをパラメータη1とする。
【０３４１】
パラメータ決定部２７’により決定されたη1は、線形予測分析部２２、非平滑化振幅スペクトル包絡推定部２３、及び平滑化振幅スペクトル包絡推定部２４及び符号化部２６に出力される。
【０３４２】
また、パラメータ決定部２７’は、決定されたη1を符号化することによりパラメータ符号を生成する。生成されたパラメータ符号は、復号装置に送信される。
【０３４３】
パラメータ決定部２７’の詳細については後述する。
【０３４４】
周波数領域変換部２１、線形予測分析部２２、非平滑化振幅スペクトル包絡系列生成部２３、平滑化振幅スペクトル包絡系列生成部２４、包絡正規化部２５及び符号化部２６は、パラメータ決定部２７が決定したパラメータη1に基づいて、第一実施形態と同様の処理により符号を生成する（ステップＡ１からステップＡ６）。この例では、符号は、線形予測係数符号と、利得符号と、整数信号符号とを合わせたものである。生成された符号は、復号装置に送信される。
【０３４５】
パラメータ決定部２７’の構成例を図１８に示す。パラメータ決定部２７’は、図１８に示すように、周波数領域変換部４１と、スペクトル包絡推定部４２と、白色化スペクトル系列生成部４３と、パラメータ取得部４４とを例えば備えている。スペクトル包絡推定部４２は、線形予測分析部４２１及び非平滑化振幅スペクトル包絡系列生成部４２２を例えば備えている。例えばこのパラメータ決定部２７’により実現されるパラメータ決定方法の各処理の例を図１９に示す。
【０３４６】
以下、図１８の各部について説明する。
【０３４７】
＜周波数領域変換部４１＞
周波数領域変換部４１には、時系列信号である時間領域の音信号が入力される。音信号の例は、音声ディジタル信号又は音響ディジタル信号である。
【０３４８】
周波数領域変換部４１は、所定の時間長のフレーム単位で、入力された時間領域の音信号を周波数領域のN点のMDCT係数列X(0),X(1),…,X(N-1)に変換する。Nは正の整数である。
【０３４９】
得られたMDCT係数列X(0),X(1),…,X(N-1)は、スペクトル包絡推定部４２及び白色化スペクトル系列生成部４３に出力される。
【０３５０】
特に断りがない限り、以降の処理はフレーム単位で行われるものとする。
【０３５１】
このようにして、周波数領域変換部４１は、音信号に対応する、例えばMDCT係数列である周波数領域サンプル列を求める（ステップＣ４１）。
【０３５２】
＜スペクトル包絡推定部４２＞
スペクトル包絡推定部４２には、周波数領域変換部２１が得たMDCT係数列X(0),X(1),…,X(N-1)が入力される。
【０３５３】
スペクトル包絡推定部４２は、所定の方法で定められるパラメータη0に基づいて、時系列信号に対応する周波数領域サンプル列の絶対値のη0乗をパワースペクトルとして用いたスペクトル包絡の推定を行う（ステップＣ４２）。
【０３５４】
推定されたスペクトル包絡は、白色化スペクトル系列生成部４３に出力される。
【０３５５】
スペクトル包絡推定部４２は、例えば以下に説明する線形予測分析部４２１及び非平滑化振幅スペクトル包絡系列生成部４２２の処理により、非平滑化振幅スペクトル包絡系列を生成することによりスペクトル包絡の推定を行う。
【０３５６】
パラメータη0は所定の方法で定められるとする。例えば、η0を０より大きい所定の数とする。例えば、η0＝１とする。また、現在パラメータηを求めようとしているフレームよりも前のフレームで求まったηを用いてもよい。現在パラメータηを求めようとしているフレーム（以下、現フレームとする。）よりも前のフレームとは、例えば現フレームのよりも前のフレームであって現フレームの近傍のフレームである。現フレームの近傍のフレームは、例えば現フレームの直前のフレームである。
【０３５７】
＜線形予測分析部４２１＞
線形予測分析部４２１には、周波数領域変換部４１が得たMDCT係数列X(0),X(1),…,X(N-1)が入力される。
【０３５８】
線形予測分析部４２１は、MDCT係数列X(0),X(1),…,X(N-1)を用いて、以下の式（C1)により定義される~R(0),~R(1),…,~R(N-1)を用いて線形予測分析を行った線形予測係数β1,β2,…,βpを生成し、生成された線形予測係数β1,β2,…,βpを符号化して線形予測係数符号と線形予測係数符号に対応する量子化された線形予測係数である量子化線形予測係数^β1,^β2,…,^βpとを生成する。
【０３５９】
【数１１】
2016167215.tif 000013
【０３６０】
生成された量子化線形予測係数^β1,^β2,…,^βpは、非平滑化スペクトル包絡系列生成部４２２に出力される。
【０３６１】
具体的には、線形予測分析部４２１は、まずMDCT係数列X(0),X(1),…,X(N-1)の絶対値のη0乗をパワースペクトルと見做した逆フーリエ変換に相当する演算、すなわち式(C1)の演算を行うことにより、MDCT係数列X(0),X(1),…,X(N-1)の絶対値のη0乗に対応する時間領域の信号列である擬似相関関数信号列~R(0),~R(1),…,~R(N-1)を求める。そして、線形予測分析部４２１は、求まった擬似相関関数信号列~R(0),~R(1),…,~R(N-1)を用いて線形予測分析を行って、線形予測係数β1,β2,…,βpを生成する。そして、線形予測分析部４２１は、生成された線形予測係数β1,β2,…,βpを符号化することにより、線形予測係数符号と、線形予測係数符号に対応する量子化線形予測係数^β1,^β2,…,^βpとを得る。
【０３６２】
線形予測係数β1,β2,…,βpは、MDCT係数列X(0),X(1),…,X(N-1)の絶対値のη0乗をパワースペクトルと見做したときの時間領域の信号に対応する線形予測係数である。
【０３６３】
線形予測分析部４２１による線形予測係数符号の生成は、例えば従来的な符号化技術によって行われる。従来的な符号化技術とは、例えば、線形予測係数そのものに対応する符号を線形予測係数符号とする符号化技術、線形予測係数をLSPパラメータに変換してLSPパラメータに対応する符号を線形予測係数符号とする符号化技術、線形予測係数をPARCOR係数に変換してPARCOR係数に対応する符号を線形予測係数符号とする符号化技術などである。
【０３６４】
このようにして、線形予測分析部４２１は、例えばMDCT係数列である周波数領域サンプル列の絶対値のη0乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を生成する（ステップＣ４２１）。
【０３６５】
なお、線形予測分析部４２１は、［線形予測符号化装置、線形予測復号装置及びこれらの方法］の欄で説明した方法により、線形予測係数符号を得て、得られた線形予測係数符号に対応する線形予測係数に変換可能な係数を量子化線形予測係数^β1,^β2,…,^βpとしてもよい。
【０３６６】
＜非平滑化振幅スペクトル包絡系列生成部４２２＞
非平滑化振幅スペクトル包絡系列生成部４２２には、線形予測分析部４２１が生成した量子化線形予測係数^β1,^β2,…,^βpが入力される。
【０３６７】
非平滑化振幅スペクトル包絡系列生成部４２２は、量子化線形予測係数^β1,^β2,…,^βpに対応する振幅スペクトル包絡の系列である非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)を生成する。
【０３６８】
生成された非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)は、白色化スペクトル系列生成部４３に出力される。
【０３６９】
非平滑化振幅スペクトル包絡系列生成部４２２は、量子化線形予測係数^β1,^β2,…,^βpを用いて、非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)として、式(C2)により定義される非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)を生成する。
【０３７０】
【数１２】
2016167215.tif 000014
【０３７１】
このようにして、非平滑化振幅スペクトル包絡系列生成部４２２は、疑似相関関数信号列に対応する振幅スペクトル包絡の系列を１／η0乗した系列である非平滑化スペクトル包絡系列を線形予測分析部４２１により生成された線形予測係数に変換可能な係数に基づいて得ることによりスペクトル包絡の推定を行う（ステップＣ４２２）。
【０３７２】
＜白色化スペクトル系列生成部４３＞
白色化スペクトル系列生成部４３には、周波数領域変換部４１が得たMDCT係数列X(0),X(1),…,X(N-1)及び非平滑化振幅スペクトル包絡生成部４２２が生成した非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)が入力される。
【０３７３】
白色化スペクトル系列生成部４３は、MDCT係数列X(0),X(1),…,X(N-1)の各係数を、対応する非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)の各値で除算することにより、白色化スペクトル系列XW(0),XW(1),…,XW(N-1)を生成する。
【０３７４】
生成された白色化スペクトル系列XW(0),XW(1),…,XW(N-1)は、パラメータ取得部４４に出力される。
【０３７５】
白色化スペクトル系列生成部４３は、例えば、k=0,1,…,N-1として、MDCT係数列X(0),X(1),…,X(N-1)の各係数X(k)を非平滑化振幅スペクトル包絡系列^H(0),^H(1),…,^H(N-1)の各値^H(k)で除算することにより、白色化スペクトル系列XW(0),XW(1),…,XW(N-1)の各値XW(k)を生成する。すなわち、k=0,1,…,N-1として、XW(k)=X(k)/^H(k)である。
【０３７６】
このようにして、白色化スペクトル系列生成部４３は、例えば非平滑化振幅スペクトル包絡系列であるスペクトル包絡で例えばMDCT係数列である周波数領域サンプル列を除算した系列である白色化スペクトル系列を得る（ステップＣ４３）。
【０３７７】
＜パラメータ取得部４４＞
パラメータ取得部４４には、白色化スペクトル系列生成部４３が生成した白色化スペクトル系列XW(0),XW(1),…,XW(N-1)が入力される。
【０３７８】
パラメータ取得部４４は、パラメータηを形状パラメータとする一般化ガウス分布が白色化スペクトル系列XW(0),XW(1),…,XW(N-1)のヒストグラムを近似するパラメータηを求める（ステップＣ４４）。言い換えれば、パラメータ取得部４４は、パラメータηを形状パラメータとする一般化ガウス分布が白色化スペクトル系列XW(0),XW(1),…,XW(N-1)のヒストグラムの分布に近くなるようなパラメータηを決定する。
【０３７９】
パラメータηを形状パラメータとする一般化ガウス分布は、例えば以下のように定義される。Γは、ガンマ関数である。
【０３８０】
【数１３】
2016167215.tif 000015
【０３８１】
一般化ガウス分布は、形状パラメータであるηを変えることにより、図２０のようにη＝１の時はラプラス分布、η＝２の時はガウス分布、といったように様々な分布を表現することができるものである。ηは、０より大きい所定の数である。ηは、０より大きい２以外の所定の数であってもよい。具体的には、ηは、２未満の所定の正の数であってよい。φは分散に対応するパラメータである。
【０３８２】
ここで、パラメータ取得部４４が求めるηは、例えば以下の式(C3)により定義される。F-1は、関数Fの逆関数である。この式は、いわゆるモーメント法により導出されるものである。
【０３８３】
【数１４】
2016167215.tif 000016
【０３８４】
逆関数F-1が定式化されている場合には、パラメータ取得部４４は、定式化された逆関数F-1にm1/((m2)1/2)の値を入力したときの出力値を計算することによりパラメータηを求めることができる。
【０３８５】
逆関数F-1が定式化されていない場合には、パラメータ取得部４４は、式(C3)で定義されるηの値を計算するために、例えば以下に説明する第一方法又は第二方法によりパラメータηを求めてもよい。
【０３８６】
パラメータηを求めるための第一方法について説明する。第一の方法では、パラメータ取得部４４は、白色化スペクトル系列に基づいてm1/((m2)1/2)を計算し、予め用意しておいた異なる複数の、ηと対応するＦ（η）のペアを参照して、計算されたm1/((m2)1/2)に最も近いＦ（η）に対応するηを取得する。
【０３８７】
予め用意しておいた異なる複数の、ηと対応するＦ（η）のペアは、パラメータ取得部４４の記憶部４４１に予め記憶しておく。パラメータ取得部４４は、記憶部４４１参照して、計算されたm1/((m2)1/2)に最も近いＦ（η）を見つけ、見つかったＦ（η）に対応するηを記憶部４４１から読み込み出力する。
【０３８８】
計算されたm1/((m2)1/2)に最も近いＦ（η）とは、計算されたm1/((m2)1/2)との差の絶対値が最も小さくなるＦ（η）のことである。
【０３８９】
パラメータηを求めるための第二方法について説明する。第二の方法では、逆関数F-1の近似曲線関数を例えば以下の式(C3’)で表される~F-1として、パラメータ取得部４４は、白色化スペクトル系列に基づいてm1/((m2)1/2)を計算し、近似曲線関数~F-1に計算されたm1/((m2)1/2)を入力したときの出力値を計算することによりηを求める。この近似曲線関数~F-1は使用する定義域において出力が正値となる単調増加関数であればよい。
【０３９０】
【数１５】
2016167215.tif 000017
【０３９１】
なお、パラメータ取得部４４が求めるηは、式(C3)ではなく、式(C3'')のように予め定めた正の整数q1及びq2を用いて（ただしq1<q2）式(C3)を一般化した式により定義されてもよい。
【０３９２】
【数１６】
2016167215.tif 000018
【０３９３】
なお、ηが式(C3'')により定義される場合も、ηが式(C3)により定義されている場合と同様の方法により、ηを求めることができる。すなわち、パラメータ取得部４４が、白色化スペクトル系列に基づいてそのq1次モーメントであるmq1とそのq2次モーメントであるmq2とに基づく値mq1/((mq2)q1/q2)を計算した後、例えば上記の第一及び第二の方法と同様、予め用意しておいた異なる複数の、ηと対応するＦ’（η）のペアを参照して、計算されたmq1/((mq2)q1/q2)に最も近いＦ’（η）に対応するηを取得するか、逆関数F’-1の近似曲線関数を~F’-1として、近似曲線関数~F-1に計算されたmq1/((mq2)q1/q2)を入力したときの出力値を計算してηを求めることができる。
【０３９４】
このようにηは次元が異なる２つの異なるモーメントmq1,mq2に基づく値であるとも言える。例えば、次元が異なる２つの異なるモーメントmq1,mq2のうち、次元が低い方のモーメントの値又はこれに基づく値（以下、前者とする。）と次元が高い方のモーメントの値又はこれに基づく値（以下、後者とする）との比の値、この比の値に基づく値、又は、前者を後者で割って得られる値に基づき、ηを求めてもよい。モーメントに基づく値とは、例えば、そのモーメントをmとしQを所定の実数としてmQのことである。また、これらの値を近似曲線関数~F-1に入力してηを求めてもよい。この近似曲線関数~F’-1は上記同様、使用する定義域において出力が正値となる単調増加関数であればよい。
【０３９５】
パラメータ決定部２７’は、ループ処理によりパラメータηを求めてもよい。すなわち、パラメータ決定部２７’は、パラメータ取得部４４で求まるパラメータηを所定の方法で定められるパラメータη0とする、スペクトル包絡推定部４２、白色化スペクトル系列生成部４３及びパラメータ取得部４４の処理を更に１回以上行ってもよい。
【０３９６】
この場合、例えば、図１８で破線で示すように、パラメータ取得部４４で求まったパラメータηは、スペクトル包絡推定部４２に出力される。スペクトル包絡推定部４２は、パラメータ取得部４４で求まったηをパラメータη0として用いて、上記説明した処理と同様の処理を行いスペクトル包絡の推定を行う。白色化スペクトル系列生成部４３は、新たに推定されたスペクトル包絡に基づいて、上記説明した処理と同様の処理を行い白色化スペクトル系列を生成する。パラメータ取得部４４は、新たに生成された白色化スペクトル系列に基づいて、上記説明した処理と同様の処理を行いパラメータηを求める。
【０３９７】
例えば、スペクトル包絡推定部４２、白色化スペクトル系列生成部４３及びパラメータ取得部４４の処理は、所定の回数であるτ回だけ更に行われてもよい。τは所定の正の整数であり、例えばτ＝１又はτ＝２である。
【０３９８】
また、スペクトル包絡推定部４２は、今回求まったパラメータηと前回求まったパラメータηとの差の絶対値が所定の閾値以下となるまで、スペクトル包絡推定部４２、白色化スペクトル系列生成部４３及びパラメータ取得部４４の処理を繰り返してもよい。
【０３９９】
（復号）
第二実施形態の復号装置及び方法は、第一実施形態と同様であるため重複説明を省略する。
【０４００】
［符号化装置、復号装置及びこれらの方法の変形例］
線形予測分析部２２及び非平滑化振幅スペクトル包絡系列生成部２３を１つのスペクトル包絡推定部２Ａとして捉えると、このスペクトル包絡推定部２Ａは、時系列信号に対応する例えばMDCT係数列である周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做したスペクトル包絡（非平滑化振幅スペクトル包絡系列）の推定を行っていると言える。ここで、「パワースペクトルと見做した」とは、パワースペクトルを通常用いるところに、η1乗のスペクトルを用いることを意味する。
【０４０１】
この場合、スペクトル包絡推定部２Ａの線形予測分析部２２は、例えばMDCT係数列である周波数領域サンプル列の絶対値のη1乗をパワースペクトルと見做した逆フーリエ変換を行うことにより得られる疑似相関関数信号列を用いて線形予測分析を行い線形予測係数に変換可能な係数を得ていると言える。また、スペクトル包絡推定部２Ａの非平滑化振幅スペクトル包絡系列生成部２３は、線形予測分析部２２により得られた線形予測係数に変換可能な係数に対応する振幅スペクトル包絡の系列を１／η1乗した系列である非平滑化スペクトル包絡系列を得ることによりスペクトル包絡の推定を行っていると言える。
【０４０２】
また、平滑化振幅スペクトル包絡系列生成部２４、包絡正規化部２５及び符号化部２６を１つの符号化部２Ｂとして捉えると、この符号化部２Ｂは、スペクトル包絡推定部２Ａにより推定されたスペクトル包絡（非平滑化振幅スペクトル包絡系列）を基にビット割り当てを変える又は実質的にビット割り当てが変わる符号化を時系列信号に対応する例えばMDCT係数列である周波数領域サンプル列の各係数に対して行っていると言える。
【０４０３】
復号部３４及び包絡逆正規化部３５を１つの復号部３Ａとして捉えると、この復号部３Ａは、非平滑化スペクトル包絡系列に基づいて変わるビット割り当て又は実質的に変わるビット割り当てに従って、入力された整数信号符号の復号を行うことにより時系列信号に対応する周波数領域サンプル列を得ていると言える。
【０４０４】
符号化部２Ｂは、スペクトル包絡（非平滑化振幅スペクトル包絡系列）を基にビット割り当てを変える又は実質的にビット割り当てが変わる符号化を行うのであれば、上記説明した算術符号化以外の符号化処理を行ってもよい。この場合、復号部３Ａは、符号化部２Ｂが行った符号化処理に対応する復号処理を行う。
【０４０５】
例えば、符号化部２Ｂは、スペクトル包絡（非平滑化振幅スペクトル包絡系列）に基づいて決定されたRiceパラメータを用いて周波数領域サンプル列をGolomb-Rice符号化してもよい。この場合、復号部３Ａは、スペクトル包絡（非平滑化振幅スペクトル包絡系列）に基づいて決定されたRiceパラメータを用いてGolomb-Rice復号してもよい。
【０４０６】
第一実施形態において、符号化装置は、パラメータηを決定する際に符号化処理を最後まで行わなくてもよい。言い換えれば、パラメータ決定部２７は、推定符号量に基づいてパラメータηを決定してもよい。この場合、符号化部２Ｂは、複数のパラメータηのそれぞれを用いて同一の所定の時間区間の時系列信号に対応する周波数領域サンプル列に対する上記と同様の符号化処理により得られる符号の推定符号量を得る。パラメータ決定部２７は、得られた推定符号量に基づいて複数のパラメータηの何れか１つを選択する。例えば、推定符号量が最も小さいパラメータηを選択する。符号化部２Ｂは、選択されたパラメータηを用いて上記と同様の符号化処理を行うことにより符号を得て出力する。
【０４０７】
上記説明した処理は、記載の順にしたがって時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。
【０４０８】
［プログラム及び記録媒体］
また、各装置又は各方法における各部をコンピュータによって実現してもよい。その場合、各装置又は各方法の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、各装置又は各方法における各部がコンピュータ上で実現される。
【０４０９】
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。
【０４１０】
また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させてもよい。
【０４１１】
このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶部に格納する。そして、処理の実行時、このコンピュータは、自己の記憶部に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実施形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよい。さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、プログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。
【０４１２】
また、コンピュータ上で所定のプログラムを実行させることにより、各装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
【図１】
2016167215.tif 000021
【図２】
2016167215.tif 000022
【図３】
2016167215.tif 000023
【図４】
2016167215.tif 000024
【図５】
2016167215.tif 000025
【図６】
2016167215.tif 000026
【図７】
2016167215.tif 000027
【図８】
2016167215.tif 000028
【図９】
2016167215.tif 000029
【図１０】
2016167215.tif 000030
【図１１】
2016167215.tif 000031
【図１２】
2016167215.tif 000032
【図１３】
2016167215.tif 000033
【図１４】
2016167215.tif 000034
【図１５】
2016167215.tif 000035
【図１６】
2016167215.tif 000036
【図１７】
2016167215.tif 000037
【図１８】
2016167215.tif 000038
【図１９】
2016167215.tif 000039
【図２０】
2016167215.tif 000040
【図２１】
2016167215.tif 000041
【図２２】
2016167215.tif 000042
【図２３】
2016167215.tif 000043
【図２４】
2016167215.tif 000044
【図２５】
2016167215.tif 000045
【図２６】
2016167215.tif 000046
【図２７】
2016167215.tif 000047
【図２８】
2016167215.tif 000048
【国際調査報告】
2016167215.tif 000049
2016167215.tif 000050