特許 6412993 符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6412993

(24)【登録日】2018年10月5日

(45)【発行日】2018年10月24日

(54)【発明の名称】符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、記録媒体

(51)【国際特許分類】

   G10L 19/035 20130101AFI20181015BHJP

   G10L 19/06 20130101ALI20181015BHJP

【ＦＩ】

   G10L19/035 Z

   G10L19/06 Z

【請求項の数】12

【全頁数】28

(21)【出願番号】特願2017-174630(P2017-174630)

(22)【出願日】2017年9月12日

(62)【分割の表示】特願2016-515878(P2016-515878)の分割

【原出願日】2015年2月20日

(65)【公開番号】特開2018-13795(P2018-13795A)

(43)【公開日】2018年1月25日

【審査請求日】2017年9月12日

(31)【優先権主張番号】特願2014-94881(P2014-94881)

(32)【優先日】2014年5月1日

(33)【優先権主張国】JP

(73)【特許権者】

【識別番号】000004226

【氏名又は名称】日本電信電話株式会社

(74)【代理人】

【識別番号】100121706

【弁理士】

【氏名又は名称】中尾 直樹

(74)【代理人】

【識別番号】100128705

【弁理士】

【氏名又は名称】中村 幸雄

(74)【代理人】

【識別番号】100147773

【弁理士】

【氏名又は名称】義村 宗洋

(72)【発明者】

【氏名】守谷 健弘

(72)【発明者】

【氏名】鎌本 優

(72)【発明者】

【氏名】原田 登

【審査官】 間宮 嘉誉

(56)【参考文献】

【文献】 特開平１０−１８７１９７（ＪＰ，Ａ）

【文献】 国際公開第２０１２／１０２１４９（ＷＯ，Ａ１）

【文献】 特開昭５８−１６８０９４（ＪＰ，Ａ）

【文献】 特開昭５９−９４７９５（ＪＰ，Ａ）

【文献】 国際公開第２０１２／１３７６１７（ＷＯ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ１０Ｌ １９／００−２５／９３

Ｓｃｉｅｎｃｅ Ｄｉｒｅｃｔ

ＩＥＥＥ Ｘｐｌｏｒｅ

Ｓｃｏｐｕｓ

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

【請求項1】

所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算部と、
前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算部と、
可変長符号化パラメータを用いて、振幅値に依存する可変長符号化により、前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化部、
を備え、
前記可変長符号化部は、
前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて符号化し、
前記周期性の程度を示す指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて符号化する
符号化装置。

【請求項2】

所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
前記周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化部、
を備え
前記可変長符号化部は、
前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化し、
前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化する
符号化装置。

【請求項3】

線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
入力された指標符号を復号して周期性の程度を示す指標を得る指標復号部と、
前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算部と、
前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算部と、
可変長符号化パラメータを用いて復号する可変長復号部、
を備え、
前記可変長復号部は、
前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて復号し、
前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算部が計算した可変長符号化パラメータを用いて復号する
復号装置。

【請求項4】

線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部と、
可変長符号を復号して周波数領域の系列を得る可変長復号部、
を備え、
前記可変長復号部は、
入力された周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得て、
前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得る
復号装置。

【請求項5】

所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算ステップと、
前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算ステップと、
可変長符号化パラメータを用いて、振幅値に依存する可変長符号化により、前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化ステップ、
を有し、
前記可変長符号化ステップは、
前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて符号化し、
前記周期性の程度を示す指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて符号化する
符号化方法。

【請求項6】

所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、前記入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
前記入力音響信号に由来する周波数領域の系列を符号化する可変長符号化ステップ、
を有し、
前記可変長符号化ステップは、
前記入力音響信号の周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化し、
前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、前記入力音響信号の振幅が大きいことを前提に、前記入力音響信号に由来する周波数領域の系列を符号化する
符号化方法。

【請求項7】

線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
入力された指標符号を復号して周期性の程度を示す指標を得る指標復号ステップと、
前記周期性統合包絡系列から振幅値に依存する可変長符号化パラメータを計算する可変長符号化パラメータ計算ステップと、
前記スペクトル包絡系列から振幅値に依存する可変長符号化パラメータを計算する第２可変長符号化パラメータ計算ステップと、
可変長符号化パラメータを用いて復号する可変長復号ステップ、
を有し、
前記可変長復号ステップは、
前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて復号し、
前記指標が前記の周期性の程度が大きいことを示す範囲ではない場合は、前記第２可変長符号化パラメータ計算ステップで計算した可変長符号化パラメータを用いて復号する
復号方法。

【請求項8】

線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成ステップと、
可変長符号を復号して周波数領域の系列を得る可変長復号ステップ、
を有し、
前記可変長復号ステップは、
入力された周期性の程度を示す指標があらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、前記周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得て、
前記指標があらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、前記スペクトル包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、前記可変長符号を復号して周波数領域の系列を得る
復号方法。

【請求項9】

請求項１または２記載の符号化装置としてコンピュータを機能させるための符号化プログラム。

【請求項10】

請求項３または４記載の復号装置としてコンピュータを機能させるための復号プログラム。

【請求項11】

請求項１または２記載の符号化装置としてコンピュータを機能させるための符号化プログラムを記録したコンピュータが読み取り可能な記録媒体。

【請求項12】

請求項３または４記載の復号装置としてコンピュータを機能させるための復号プログラムを記録したコンピュータが読み取り可能な記録媒体。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、音響信号のスペクトル包絡を利用して音響信号を符号化、復号化する符号化装置、復号装置、符号化方法、復号方法、符号化プログラム、復号プログラム、および記録媒体に関する。

【背景技術】

【0002】

低ビット（例えば10kbit/s〜20kbit/s程度）の音声信号や音響信号の符号化方法として、ＤＦＴ（離散フーリエ変換）やＭＤＣＴ（変形離散コサイン変換）などの直交変換係数に対する適応符号化が知られている。例えば非特許文献１で用いられているＴＣＸ（transform coded excitation：変換符号化励振）符号化方法では、入力された音信号の周波数領域表現である係数列X[1]，…，X[N]から振幅スペクトル包絡の影響を取り除いた系列（正規化係数列X_N[1]，…，X_N[N]）を求め、これを可変長符号化する。ただし、[]内のＮは正整数である。

【0003】

振幅スペクトル包絡は、以下の手順で算出される。
（ｓｔｅｐ１）所定の時間区間であるフレーム単位で、入力された時間領域の音響ディジタル信号（以下、入力音響信号）に対する線形予測分析を行って線形予測係数α₁，…，α_Pを求める。ただし、Ｐは予測次数を示す正整数である。例えば、全極型モデルであるＰ次自己回帰過程により、時刻ｔでの入力音響信号x(t)は、Ｐ時点まで遡った過去の自分自身の値x(t-1)，…，x(t-P)と予測残差e(t)と線形予測係数α₁，…，α_pによって式（１）で表される。
x(t)=α₁x(t-1)+…+α_p x(t-P)+e(t) (1)

【0004】

（ｓｔｅｐ２）線形予測係数α₁，…，α_Pを量子化し、量子化済線形予測係数^α₁，…，^α_Pを求める。量子化済線形予測係数^α₁，…，^α_Pを用いてＮ点の入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、式（２）で求めることができる。ただし、ｎは1≦n≦Nの整数、exp（・）はネイピア数を底とする指数関数、ｊは虚数単位、σは予測残差信号の振幅である。

【数1】

【0005】

なお、本明細書では、右肩に角括弧なしで表記されている記号はべき乗演算を表す。つまり、σ²はσの２乗を表す。また、文中で使用する記号「~」「＾」等は、本来直後の文字の真上に記載されるべきものであるが、テキスト記法の制限により、当該文字の直前に記載する。数式中においてはこれらの記号は本来の位置、すなわち文字の真上に記述している。

【先行技術文献】

【非特許文献】

【0006】

【非特許文献1】Anthony Vetro, “MPEG Unified Speech and Audio Coding”, Industry and Standards, IEEE MultiMedia, April-June, 2013.

【発明の概要】

【発明が解決しようとする課題】

【0007】

音響信号の符号化では、復号側でもスペクトル包絡の情報を得るために、スペクトル包絡に対応する符号を復号側へ送る必要がある。非特許文献１のように線形予測係数によりスペクトル包絡を求める場合には、復号側へ送る「スペクトル包絡に対応する符号」は「線形予測係数に対応する符号」であり、符号量が少なくて済むという利点がある。一方、線形予測係数により求めたスペクトル包絡の情報は、入力音響信号のピッチ周期に起因するピークの付近での近似精度が悪くなることがあり、このことが正規化係数列を可変長符号化する際の符号化効率の低下につながることがある。

【0008】

このような問題に鑑み、本発明では、音響信号のピッチ周期に起因するピークの付近での近似精度がよい周波数領域の包絡系列を少ない符号量で再現できる符号化装置、復号装置を提供する。

【課題を解決するための手段】

【0009】

本発明の符号化装置は、周期性統合包絡生成部と可変長符号化部を備える。周期性統合包絡生成部は、所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する。可変長符号化部は、周期性統合包絡系列の値が大きい周波数ほど、入力音響信号の振幅が大きいことを前提に、入力音響信号に由来する周波数領域の系列を符号化する。本発明の復号装置は、周期性統合包絡生成部と可変長復号部を備える。周期性統合包絡生成部は、線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する。可変長復号部は、周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、可変長符号を復号して周波数領域の系列を得る。

【発明の効果】

【0010】

本発明の符号化装置と復号装置によれば、入力音響信号のピッチ周期に起因するピーク付近での近似精度がよい周波数領域の包絡系列を少ない符号量で再現できる符号化、復号を行うので、音響信号の符号化効率がよい。

【図面の簡単な説明】

【0011】

【図1】実施例１の周期性統合包絡系列生成装置の機能構成例を示す図。

【図2】実施例１の周期性統合包絡系列生成装置の処理フローを示す図。

【図3】周期性包絡系列P[1]，…，P[N]の例を示す図。

【図4A】同じ音響信号に対して生成された系列の違いを説明するための例を示す図であって、数列X[1]，…，X[N]を補間した曲線の形状を示す図。

【図4B】同じ音響信号に対して生成された系列の違いを説明するための例を示す図であって、周期性包絡系列P[1]，…，P[N]を補間した曲線の形状を示す図。

【図4C】同じ音響信号に対して生成された系列の違いを説明するための例を示す図であって、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を補間した曲線の形状を示す図。

【図4D】同じ音響信号に対して生成された系列の違いを説明するための例を示す図であって、周期性統合包絡系列W_M[1]，…，W_M[N]を補間した曲線の形状を示す図。

【図5】実施例２の符号化装置の機能構成例を示す図。

【図6】実施例２の符号化装置の処理フローを示す図。

【図7】実施例２の復号装置の機能構成例を示す図。

【図8】実施例２の復号装置の処理フローを示す図。

【図9】実施例３の符号化装置の機能構成例を示す図。

【図10】実施例３の符号化装置の処理フローを示す図。

【図11】実施例３の復号装置の機能構成例を示す図。

【図12】実施例３の復号装置の処理フローを示す図。

【発明を実施するための形態】

【0012】

以下、本発明の実施の形態について、詳細に説明する。なお、同じ機能を有する構成部には同じ番号を付し、重複説明を省略する。

【実施例1】

【0013】

図１に本発明の周期性統合包絡系列生成装置の機能構成例を、図２に本発明の周期性統合包絡系列生成装置の処理フローを示す。周期性統合包絡系列生成装置１００は、スペクトル包絡系列計算部１２０、周波数領域変換部１１０、周期性分析部１３０、周期性包絡系列生成部１４０、周期性統合包絡生成部１５０を備え、入力された時間領域の音響ディジタル信号を入力音響信号x(t)とし、係数列の周波数成分に基づいて振幅スペクトル包絡系列を変形した周期性統合包絡系列を生成する。

【0014】

＜スペクトル包絡系列計算部１２０＞
スペクトル包絡系列計算部１２０は、入力音響信号x(t)の時間領域の線形予測に基づき、入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を計算する（Ｓ１２０）。ただし、Ｎは正整数である。スペクトル包絡系列計算部１２０は、従来技術と同じであり、以下の手順で計算すればよい。

【0015】

（ｓｔｅｐ１）所定の時間区間であるフレーム単位で、入力音響信号に対する線形予測分析を行って線形予測係数α₁，…，α_Pを求める。ただし、Ｐは予測次数を示す正整数である。例えば、全極型モデルであるＰ次自己回帰過程により、時刻ｔでの入力音響信号x(t)は、Ｐ時点まで遡った過去の自分自身の値x(t-1)，…，x(t-P)と予測残差e(t)と線形予測係数α₁，…，α_pによって式（１）で表される。

【0016】

（ｓｔｅｐ２）線形予測係数α₁，…，α_Pを用いてＮ点の入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、線形予測係数α₁，…，α_Pに対応する量子化済線形予測係数^α₁，…，^α_Pを用いて式（２）で求めることができる。または、振幅スペクトル包絡系列の各値W[n]は、線形予測係数α₁，…，α_Pを用いて、式（２）の^α_pをα_pに置き換えた式で求めることができる。

【0017】

＜周波数領域変換部１１０＞
周波数領域変換部１１０は、所定の時間区間であるフレーム単位で、入力された時間領域の入力音響信号を周波数領域のＮ点の係数列X[1]，…，X[N]に変換して出力する（Ｓ１１０）。周波数領域への変換は、ＭＤＣＴ（変形離散コサイン変換）やＤＦＴ（離散フーリエ変換）などの方法で行えばよい。

【0018】

＜周期性分析部１３０＞
周期性分析部１３０は、係数列X[1]，…，X[N]を入力とし、当該係数列X[1]，…，X[N]の周期Ｔを求め、周期Ｔを出力する（Ｓ１３０）。

【0019】

周期Ｔは、入力音響信号に由来する周波数領域の係数列、例えば、係数列X[1]，…，X[N]、の周期性を有する成分の間隔（係数列が周期的に大きな値となる間隔）に対応する情報である。以下では周期Ｔを間隔Ｔと表現する場合もあるが、表現上の違いだけであり、同じものである。Ｔは正値であり、整数であってもよいし、小数（例えば、5.0、5.25、5.5、5.75）であってもよい。

【0020】

また、周期性分析部１３０は、必要に応じて、係数列X[1]，…，X[N]を入力とし、周期性の程度を示す指標Ｓも求めて出力してもよい。この場合、例えば、係数列X[1]，…，X[N]の周期性を有する成分の部分のエネルギーとそれ以外の部分のエネルギーとの比など基づいて周期性の程度を示す指標Ｓを求める。この場合は、指標Ｓは周波数領域のサンプル列の周期性の程度を示す指標となる。なお、周期性を有する成分の大きさが大きいほど、すなわち、周期Ｔの整数倍のサンプルやその近傍にあるサンプルの振幅（サンプル値の絶対値）が大きいほど、周波数領域のサンプル列の「周期性の程度」は大きい。

【0021】

なお、周期性分析部１３０は、時間領域の入力音響信号から時間領域の周期を求め、求めた時間領域の周期を周波数領域の周期に変換することで、周期Ｔを求めてもよい。また、時間領域の周期を周波数領域の周期に変換したものの定数倍やその近傍の値を周期Ｔとして求めてもよい。同様に、周期性分析部１３０は、時間領域の入力音響信号から、例えば、時間領域の周期分だけ時間がずれた信号列間の相関の大きさ等に基づいて、周期性の程度を示す指標Ｓを求めてもよい。

【0022】

要は、時間領域の入力音響信号やそれに由来する周波数領域係数列から周期Ｔや指標Ｓを求める方法は、従来より様々な方法が存在するので、その何れの方法を選択して利用してもよい。

【0023】

＜周期性包絡系列生成部１４０＞
周期性包絡系列生成部１４０は、間隔Ｔを入力とし、周期性包絡系列P[1]，…，P[N]を出力する（Ｓ１４０）。周期性包絡系列P[1]，…，P[N]は、ピッチ周期に起因する周期でピークを持つ周波数領域の離散系列、すなわち調波モデルに対応する離散系列である。図３に周期性包絡系列P[1]，…，P[N]の例を示す。周期性包絡系列P[1]，…，P[N]は、図３に示された波形のように、間隔Ｔの整数倍の近傍の整数値であるインデックスと、その前後所定数のインデックスに対応する周期性包絡の値のみ正の値を持ち、それ以外は０であるような系列である。間隔Ｔの整数倍の近傍の整数値であるインデックスが周期的に最大値（ピーク）をとり、その前後所定数のインデックスに対応するP[n]の値は、そのインデックスｎがピークに対応するインデックスから離れるにつれて単調減少する関係にある。図３の横軸の1,2,…,は離散化サンプル点のインデックス（以下、「周波数インデックス」）を表す。

【0024】

例えば、ｎを周波数インデックスを表す変数とし、τを極大値（ピーク）に対応する周波数インデックスとして、ピークの形状は以下の関数Q(n)で表せる。ただし、間隔Ｔの小数点以下の桁数がＬ桁であり、間隔Ｔ’をＴ’＝Ｔ×２^Ｌとする。

【数2】

ｈはピークの高さを表し、間隔Ｔが大きいほどピークの高さが高くなる。また、ＰＤはピーク部分の幅を表し、間隔Ｔが大きいほど幅が広くなる。

【0025】

Ｕを１からピークの数までを示す正整数（例えば、図３の場合は１〜１０）とし、ｖを１以上の整数（例えば、１から３程度）とし、floor(・)を小数点以下を切り捨てて整数値を返す関数とすると、周期性包絡系列P[n]は、例えば、

【数3】

のように計算すればよい。ただし、（Ｕ×Ｔ’）／２^Ｌ−ｖ≦ｎ≦（Ｕ×Ｔ’）／２^Ｌ＋ｖである。例えば、L=2の場合、T=20.00であればT’=80、T=20.25であればT’=81、T=20.50であればT’=82、T=20.75であればT’=83である。なお、周期性包絡系列P[n]は、小数点第一位を四捨五入して整数値を返す関数Round(・)を用いて、

【数4】

のように求めてもよい。

【0026】

＜周期性統合包絡生成部１５０＞
周期性統合包絡生成部１５０は、少なくとも、周期性包絡系列P[1]，…，P[N]、振幅スペクトル包絡系列W[1]，…，W[N]を入力とし、周期性統合包絡系列W_M[1]，…，W_M[N]を求める（Ｓ１５０）。具体的には、周期性統合包絡W_M[n]を次式のように求める。

【数5】

なお、δは、周期性統合包絡W_M[n]と係数X[n]の絶対値系列の形状が近くなるように決定される値または予め定めた値である。

【0027】

周期性統合包絡生成部１５０において周期性統合包絡W_M[n]と係数X[n]の絶対値系列の形状が近くなるようにδを決定する場合には、周期性統合包絡生成部１５０は、係数列X[1], …, X[N]も入力とし、決定されたδとそのときの周期性統合包絡系列W_M[1]，…，W_M[N]を出力すればよい。例えば、δは、いくつかのδの候補、例えば、０．４と０．８の２つをδの候補、の中から以下の式により定義されるＥが最小となるδに決めればよい。言い換えると、周期性統合包絡W_M[n]と係数X[n]の絶対値系列の形状が近くなるδに決めればよい。

【数6】

【0028】

δは、周期性統合包絡W_M[n]において周期性包絡P[n]をどの程度考慮するかを決める値である。言い換えれば、δは周期性統合包絡W_M[n]における振幅スペクトル包絡W[n]と周期性包絡P[n]の混合比率を決める値といえる。また、式（９）のＧは係数列X[1]，…，X[N]の各係数X[n]の絶対値の系列と周期性統合包絡系列の逆数の系列との内積である。式（８）の~W_M[n]は、周期性統合包絡の各値W_M[n]をＧで正規化した正規化周期性統合包絡である。式（７）において、係数列X[1]，…，X[N]と正規化周期性統合包絡系列~W_M[1]，…，~W_M[N]の内積の４乗を計算しているのは、特に絶対値の大きい係数X[n]を強調して内積をとった値（距離）を小さくすることを意図している。つまり、係数列X[1]，…，X[N]の中で特に絶対値の大きい係数X[n]と周期性統合包絡W_M[n]が近くなるようにδを決定することを意味している。

【0029】

また、周期性統合包絡生成部１５０において周期性の程度に応じてδの候補数を決定する場合には、周期性統合包絡生成部１５０は、周期性の程度を示す指標Ｓも入力とし、指標Ｓが、周期性が高いことに対応するフレームであることを示している場合には多くの候補数のδの候補の中から式（７）で定義されるＥが最小となるδを選び、指標Ｓが、周期性が低いことに対応するフレームであることを示している場合にはδを予め定めた値としてもよい。すなわち、周期性統合包絡生成部１５０において周期性の程度に応じてδの候補数を決定する場合には、周期性が高いほどδの候補の数を多くすればよい。

【0030】

＜実施例１の発明の効果＞
図４Ａ〜４Ｄに同じ音響信号に対して生成された系列の違いを説明するための例を示す。図４Ａに係数列X[1]，…，X[N]を補間した曲線の形状を、図４Ｂに周期性包絡系列P[1]，…，P[N]を補間した曲線の形状を、図４Ｃに平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を補間した曲線の形状を、図４Ｄに周期性統合包絡系列W_M[1]，…，W_M[N]を補間した曲線の形状を示す。図４Ａ〜４Ｄに示すとおり、周期性統合包絡系列W_M[1]，…，W_M[N]は、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]に比べて、係数列X[1]，…，X[N]に現れる周期的なピークを含んだ形状となっている。また、周期性統合包絡系列W_M[1]，…，W_M[N]は、スペクトル包絡を表す情報である線形予測係数または量子化済線形予測係数の他に、間隔Ｔ、または、間隔Ｔと値δの情報があれば生成できる。したがって、入力音響信号のスペクトル包絡を表す情報に少ない情報量を追加するだけで、入力音響信号のピッチ周期に起因する振幅のピークを、線形予測係数により求まるスペクトル包絡より高精度に表現することができる。すなわち、線形予測係数または量子化済線形予測係数と、間隔Ｔ、または、間隔Ｔと値δと、の少ない情報量で入力音響信号の振幅を高精度に推定することができることになる。なお、平滑化振幅スペクトル包絡~W[n]は次式で表現される包絡であり、γは振幅スペクトル係数を鈍らせる（平滑化する）ための１以下の正の定数である。

【数7】

【0031】

また、本発明の周期性統合包絡系列生成装置を符号化装置と復号装置で用いる場合には、符号化装置に含まれる周期性統合包絡系列生成装置以外の処理部で得られた量子化済線形予測係数^α_ｐを特定する符号（線形予測係数符号C_L）と周期Ｔや時間領域の周期を特定する符号（周期符号C_T）が復号装置に入力されるので、本発明の周期性統合包絡系列生成装置からはδの情報を示す符号を出力すれば、復号側の周期性統合包絡系列生成装置でも符号化側の周期性統合包絡系列生成装置で生成した周期性統合包絡系列と同じ周期性統合包絡系列を生成できる。したがって、符号化装置から復号装置に符号を送る際に増加する符号量は少ない。

【0032】

＜実施例１の発明のポイント＞
実施例１の周期性統合包絡系列生成装置１００では、周期性統合包絡生成部１５０が係数列X[1]，…，X[N]の周期性成分に基づいて、振幅スペクトル包絡系列W[1]，…，W[N]を変形し、周期性統合包絡系列W_M[1]，…，W_M[N]としている点が最も重要なポイントである。特に、係数列X[1]，…，X[N]の周期性の程度が大きいほど、すなわち、周期性を有する成分の大きさが大きいほど、振幅スペクトル包絡系列W[1]，…，W[N]のうち間隔Ｔ（周期）の整数倍およびそれらの近傍のサンプルの値を大きく変更すれば、上記の効果を得やすい。「近傍のサンプル」とは、間隔Ｔの整数倍の近傍の整数値であるインデックスで示されるサンプルである。また、「近傍」とは、例えば、式（３）〜（５）などのあらかじめ定めた方法で決まる範囲とすればよい。

【0033】

また、係数列X[1]，…，X[N]の周期性を有する成分の間隔Ｔが広いほど、式（４）と式（５）に示された周期性包絡系列P[1]，…，P[N]は、大きい値を持ち、広い幅で、すなわち、間隔Ｔ（周期）の整数倍およびそれらの近傍の多くのサンプルで、０以外の値を持つ。つまり、周期性統合包絡生成部１５０は、係数列の周期性を有する成分の間隔Ｔが広いほど、振幅スペクトル包絡系列のうち間隔Ｔ（周期）の整数倍およびそれらの近傍のサンプルの値を大きく変更する。また、周期性統合包絡生成部１５０は、係数列の周期性を有する成分の間隔Ｔが広いほど、振幅スペクトル包絡系列を広い幅で、すなわち、間隔Ｔ（周期）の整数倍およびそれらの近傍の多くのサンプルで、サンプル値を変更する。「近傍の多くのサンプルで」とは、「近傍」に該当する範囲（あらかじめ定めた方法で決まる範囲）に存在するサンプルを多くすることを意味している。つまり、周期性統合包絡生成部１５０は、このように振幅スペクトル包絡系列を変形すれば、上記の効果を得やすい。

【0034】

なお、周期性統合包絡系列が持つ「入力音響信号のピッチ周期に起因する振幅のピークをより高精度に表現することができる。」という特徴を効果的に利用する例としては、符号化装置と復号装置があり、この例を実施例２，３に示している。ただし、周期性統合包絡系列の特徴の利用例は、符号化装置と復号装置以外にも、雑音除去装置やポストフィルタなどがあり得る。したがって、実施例１では周期性統合包絡系列生成装置を説明している。

【0035】

［変形例１］（正規化係数列で周期性分析する例）
変形例１の周期性統合包絡系列生成装置も図１に示す。また、変形例１の周期性統合包絡系列生成装置の処理フローも図２に示す。周期性統合包絡系列生成装置１０１は、周波数領域系列正規化部１１１も備える点と、スペクトル包絡系列計算部１２１、周期性分析部１３１が周期性統合包絡系列生成装置１００と異なり、その他の構成は同じである。以下では相違点についてのみ説明する。

【0036】

＜スペクトル包絡系列計算部１２１＞
スペクトル包絡系列計算部１２１は、振幅スペクトル包絡系列W[1]，…，W[N]だけではなく、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]も求める。

【0037】

具体的には、スペクトル包絡系列計算部１２１は、スペクトル包絡系列計算部１２０で示した（ｓｔｅｐ１），（ｓｔｅｐ２）に加えて、以下の手順の処理を行う。

【0038】

（ｓｔｅｐ３）量子化済線形予測係数^α_ｐのそれぞれにγ^ｐを乗算し、量子化済平滑化線形予測係数^α₁γ，^α₂γ²，…，^α_Pγ^Pを求める。γは平滑化するための１以下の正の定数である。そして、式（１０）によって、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を求める（Ｓ１２１）。もちろん、スペクトル包絡系列計算部１２０と同様に、量子化済線形予測係数^α_ｐに代えて線形予測係数α_ｐを用いてもよい。

【0039】

＜周波数領域系列正規化部１１１＞
周波数領域系列正規化部１１１は、係数列X[1]，…，X[N]の各係数を平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]の各係数で除算して正規化係数列X_N[1]，…，X_N[N]を得る。すなわち、ｎ＝１，…，Ｎに対して
X_N[n]=X[n]/~W[n] （１１）
の計算を行い、正規化係数列X_N[1]，…，X_N[N]を求める（Ｓ１１１）。

【0040】

＜周期性分析部１３１＞
周期性分析部１３１は、正規化係数列X_N[1]，…，X_N[N]を入力とし、当該正規化係数列X_N[1]，…，X_N[N]の周期Ｔを求め、周期Ｔを出力する（Ｓ１３１）。すなわち、本変形例では、入力音響信号に由来する周波数領域の係数列である正規化係数列X_N[1]，…，X_N[N]の周期性を有する成分の間隔を周期Ｔとして求める。また、周期性分析部１３１は、必要に応じて、係数列X[1]，…，X[N]を入力とし、周期性の程度を示す指標Ｓも求めて出力してもよい。

【0041】

その他の処理は周期性統合包絡系列生成装置１００と同じである。したがって、実施例１と同様の効果が得られる。なお、周期性統合包絡系列生成装置１０１の場合は、周期性統合包絡生成部１５０は、振幅スペクトル包絡系列W[1]，…，W[N]の代わりに平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を用いてもよい。この場合は、式（６）の代わりに次式の計算となる。

【数8】

【0042】

［変形例２］（外部から情報が入力される例）
本発明の周期性統合包絡系列生成装置を符号化装置や復号装置が内部に備えている場合には、符号化装置や復号装置に含まれる周期性統合包絡系列生成装置以外の処理部で、係数列X[1]，…，X[N]、正規化係数列X_N[1]，…，X_N[N]、量子化済線形予測係数^α_ｐ、量子化済平滑化線形予測係数^α_ｐγ^ｐ、振幅スペクトル包絡W[1]，…，W[N]、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]、周期Ｔ、指標Ｓなどが求められていることがある。このような場合は、周期性統合包絡系列生成装置に、周波数領域変換部、周波数領域正規化部、スペクトル包絡系列計算部、周期性分析部の少なくとも何れかを備えない構成としてもよい。この場合には、符号化装置内の周期性統合包絡系列生成装置以外の処理部から、量子化済線形予測係数^α_ｐを特定する符号（線形予測係数符号C_L）、周期Ｔや時間領域の周期を特定する符号（周期符号C_T）、指標Ｓを特定する符号、などが出力され、復号装置に入力される。したがって、この場合には、符号化装置内の周期性統合包絡系列生成装置からは、量子化済線形予測係数^α_ｐを特定する符号（線形予測係数符号C_L）、周期Ｔや時間領域の周期を特定する符号（周期符号C_T）、指標Ｓを特定する符号、などを出力する必要がない。

【0043】

また、本発明の周期性統合包絡系列生成装置を符号化装置や復号装置で用いる場合には、符号化装置と復号装置とで同一の周期性統合包絡系列を得られるようにする必要がある。したがって、符号化装置が出力し復号装置に入力される符号から特定可能な情報を用いて周期性統合包絡系列を得る必要がある。たとえば、符号化装置で用いる周期性統合包絡系列生成装置のスペクトル包絡系列計算部では、線形予測係数符号C_Lに対応する量子化済線形予測係数を用いて振幅スペクトル包絡系列を求め、復号装置で用いる周期性統合包絡系列生成装置のスペクトル包絡系列計算部では、符号化装置から出力されて復号装置に入力される線形予測係数符号C_Lに対応する復号線形予測係数を用いて振幅スペクトル包絡系列を求める必要がある。

【0044】

なお、符号化装置や復号装置で周期性統合包絡系列を用いる場合には、上述のように周期性統合包絡系列生成装置を内部に備えるのではなく、周期性統合包絡系列生成装置内の必要な処理部を符号化装置と復号装置に備えるようにすればよい。このような符号化装置や復号装置は実施例２で説明する。

【実施例2】

【0045】

≪符号化装置≫
図５に実施例２の符号化装置の機能構成例を、図６に実施例２の符号化装置の処理フローを示す。符号化装置２００は、スペクトル包絡系列計算部２２１、周波数領域変換部１１０、周波数領域系列正規化部１１１、周期性分析部２３０、周期性包絡系列生成部１４０、周期性統合包絡生成部２５０、可変長符号化パラメータ計算部２６０、可変長符号化部２７０を備える。符号化装置２００は、入力された時間領域の音響ディジタル信号を入力音響信号x(t)とし、少なくとも量子化済線形予測係数^α₁，…，^α_Pを示す符号Ｃ_Ｌ、正規化係数列X_N[1]，…，X_N[N]の周期を表す間隔Ｔの符号Ｃ_Ｔ、正規化係数列X_N[1]，…，X_N[N]を可変長符号化した可変長符号Ｃ_Ｘを出力する。周波数領域系列正規化部１１１は実施例１変形例１と同じである。周波数領域変換部１１０と周期性包絡系列生成部１４０は実施例１と同じである。以下では異なる構成部について説明する。

【0046】

＜スペクトル包絡系列計算部２２１＞
スペクトル包絡系列計算部２２１は、入力音響信号x(t)の時間領域の線形予測に基づき、入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を計算し、計算の過程で得た量子化済線形予測係数^α₁，…，^α_Pを示す符号Ｃ_Ｌも求める（Ｓ２２１）。ただし、Ｎは正整数である。スペクトル包絡系列計算部２２１は、以下の手順で処理すればよい。

【0047】

（ｓｔｅｐ１）所定の時間区間であるフレーム単位で、入力音響信号に対する線形予測分析を行って線形予測係数α₁，…，α_Pを求める。ただし、Ｐは予測次数を示す正整数である。例えば、全極型モデルであるＰ次自己回帰過程により、時刻ｔでの入力音響信号x(t)は、Ｐ時点まで遡った過去の自分自身の値x(t-1)，…，x(t-P)と予測残差e(t)と線形予測係数α₁，…，α_pによって式（１）で表される。

【0048】

（ｓｔｅｐ２）線形予測係数α₁，…，α_Pを符号化して符号Ｃ_Ｌを得て出力するとともに、符号Ｃ_Ｌに対応する量子化済線形予測係数^α₁，…，^α_Pを求める。また、量子化済線形予測係数^α₁，…，^α_Pを用いてＮ点の入力音響信号の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、式（２）で求めることができる。なお、線形予測係数α₁，…，α_Pを符号化して符号Ｃ_Ｌを得る方法は、線形予測係数をLSPパラメータに変換して、LSPパラメータを符号化して符号Ｃ_Ｌを得るなど、線形予測係数に変換可能な係数の何れを符号化して符号Ｃ_Ｌを得る何れの方法を用いてもよい。

【0049】

（ｓｔｅｐ３）量子化済線形予測係数^α_ｐのそれぞれにγ^ｐを乗算し、量子化済平滑化線形予測係数^α₁γ，^α₂γ²，…，^α_Pγ^Pを求める。γはあらかじめ定めた平滑化するための１以下の正の定数である。そして、式（１０）によって、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を求める。

【0050】

＜周期性分析部２３０＞
周期性分析部２３０は、正規化係数列X_N[1]，…，X_N[N]を入力とし、当該正規化係数列X_N[1]，…，X_N[N]の間隔Ｔ（周期的に大きな値となる間隔）を求め、間隔Ｔと間隔Ｔを示す符号Ｃ_Ｔを出力する（Ｓ２３０）。また、周期性分析部２３０は、必要に応じて、周期性の程度を示す指標Ｓ（すなわち、周波数領域のサンプル列の周期性の程度を示す指標）、も求めて出力する。また、周期性分析部２３０は、必要に応じて、指標Ｓを示す符号Ｃ_Ｓも得て出力する。なお、指標Ｓと間隔Ｔ自体は実施例１変形例１の周期性分析部１３１と同じである。

【0051】

＜周期性統合包絡生成部２５０＞
周期性統合包絡生成部２５０は、少なくとも、周期性包絡系列P[1]，…，P[N]、振幅スペクトル包絡系列W[1]，…，W[N]を入力とし、周期性統合包絡系列W_M[1]，…，W_M[N]を求めて周期性統合包絡W_M[n]を出力する。また、周期性統合包絡生成部２５０は、値δとして、予め定めた１つの値ではなく、予め定めた複数の候補値のうちの何れかを選択する場合には、係数列X[1], …, X[N]も入力とし、予め定めた複数の候補値のうち周期性統合包絡W_M[n]と係数X[n]の絶対値系列の形状が近くなる候補値を値δとして求め、値δを示す符号Ｃ_δも出力する（Ｓ２５０）。

【0052】

周期性統合包絡W_M[n]と値δは実施例１と同じであり、周期性統合包絡W_M[n]は式（６），…，（９）のように求めればよい。周期性統合包絡生成部２５０において周期性の程度に応じてδの候補数を決定する場合には、周期性統合包絡生成部２５０は、周期性の程度を示す指標Ｓも入力とし、指標Ｓが周期性が高いことに対応するフレームの場合には多くの候補数のδの候補の中から式（７）で定義されるＥが最小となるδを選び、指標Ｓが周期性が低いことに対応するフレームである場合にはδを１つの予め定めた値としてもよい。なお、δを予め定めた値にする場合は、値δを示す符号Ｃ_δを出力する必要はない。

【0053】

＜可変長符号化パラメータ計算部２６０＞
可変長符号化パラメータ計算部２６０は、周期性統合包絡系列W_M[1]，…，W_M[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と正規化係数列X_N[1]，…，X_N[N]を入力とし、可変長符号化パラメータｒ_ｎを求める（Ｓ２６０）。可変長符号化パラメータ計算部２６０は、周期性統合包絡系列W_M[1]，…，W_M[N]から求めた振幅値に依存して可変長符号化パラメータｒ_ｎを計算することを特徴としている。

【0054】

可変長符号化パラメータは、符号化対象の信号、すなわち、正規化係数列X_N[1]，…，X_N[N]の各係数の振幅の取り得る範囲を特定するパラメータである。例えば、ライス符号化の場合にはライスパラメータが可変長符号化パラメータに相当し、算術符号化の場合は符号化対象の信号の振幅の取り得る範囲が可変長符号化パラメータに相当する。

【0055】

１サンプルごとに可変長符号化を行う場合には、正規化係数列の各係数X_N[n]について可変長符号化パラメータが計算される。複数のサンプルからなるサンプル群ごとに（例えば２サンプルずつ）まとめて可変長符号化を行う場合には、サンプル群ごとに可変長符号化パラメータが計算される。つまり、可変長符号化パラメータ計算部２６０は、正規化係数列の一部である正規化部分係数列ごとに、可変長符号化パラメータｒ_ｎを計算する。ここで、正規化部分係数列は複数個あり、複数個の正規化部分係数列には正規化係数列の係数が重複されずに含まれるものとする。以下に、１サンプルごとにライス符号化を行う場合を例に、可変長符号化パラメータの計算方法を説明する。

【0056】

（ｓｔｅｐ１）正規化係数列X_N[1],…,X_N[N]の各係数の振幅の平均の対数を、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）として次式のように算出する。

【数9】

ｓｂはフレームごとに１度だけ符号化されて、基準となるライスパラメータ（基準となる可変長符号化パラメータ）に対応する符号Ｃ_ｓｂとして復号装置４００に伝送される。あるいは復号装置４００に伝送される別の情報から正規化係数列X_N[1],…,X_N[N]の振幅の平均値を推定できる場合は、符号化装置２００と復号装置４００で共通に振幅の平均値の推定値からｓｂを近似的に決定する方法を決めておいてもよい。例えば、包絡の傾きを表すパラメータ、区分帯域ごとの平均包絡の大きさを表すパラメータを別途使う符号化の場合には、復号装置４００に伝送される別の情報から振幅の平均値を推定できる。この場合は、ｓｂを符号化し、基準となるライスパラメータに対応する符号Ｃ_ｓｂを復号装置４００へ出力しなくてもよい。

【0057】

（ｓｔｅｐ２）下記式により閾値θを算出する。

【数10】

θは、周期性統合包絡系列の各値W_M[n]を平滑化振幅スペクトル包絡系列の各値~W[n]で除算した値の振幅の平均の対数である。

【0058】

（ｓｔｅｐ３） |W_M[n]/~W[n]|がθより大きいほど、正規化係数X_N[n]をライス符号化するためのライスパラメータｒ_ｎをｓｂよりも大きな値として決定する。|W_M[n]/~W[n]|がθより小さいほど、正規化係数X_N[n]をライス符号化するためのライスパラメータｒ_ｎをｓｂよりも小さな値として決定する。

【0059】

（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各X_N[n]についてのライスパラメータｒ_ｎを求める。

【0060】

＜可変長符号化部２７０＞
可変長符号化部２７０は、可変長符号化パラメータ計算部２６０で求めた可変長符号化パラメータｒ_ｎを用いて正規化係数列X_N[1]，…，X_N[N]を可変長符号化し、可変長符号Ｃ_Ｘを出力する（Ｓ２７０）。例えば、可変長符号化部２７０は、可変長符号化パラメータ計算部２６０で求めたライスパラメータｒ_ｎを用いて正規化係数列X_N[1]，…，X_N[N]をライス符号化し、得られた符号を可変長符号Ｃ_Ｘとして出力する。可変長符号化パラメータ計算部２６０で求めたライスパラメータｒ_ｎは、周期性統合包絡系列の振幅値に依存する可変長符号化パラメータであり、周期性統合包絡系列の値が大きい周波数ほど大きな値となっている。ライス符号化は、振幅値に依存する可変長符号化の公知技術のうちの１つであり、ライスパラメータｒ_ｎを用いて振幅値に依存する可変長符号化を行うものである。また、周期性統合包絡生成部２５０で生成した周期性統合包絡系列は、入力音響信号のスペクトル包絡を高精度に表現するものである。すなわち、可変長符号化部２７０は、周期性統合包絡系列の値が大きい周波数ほど、前記入力音響信号の周波数領域の係数列であるX[1]，…，X[N]の振幅が大きいとことを前提に、正規化係数列X_N[1]，…，X_N[N]を可変長符号化していることになり、言い換えれば、可変長符号化パラメータを用いて、振幅値に依存する可変長符号化により、正規化係数列X_N[1]，…，X_N[N]を符号化していることになる。ここでいう振幅値とは、符号化対象の係数列の平均振幅値、係数列に含まれる各係数の振幅の推定値、係数列の振幅の包絡の推定値などである。

【0061】

符号化装置２００は、このような処理によって得られた量子化済線形予測係数^α₁，…，^α_Pを示す符号Ｃ_Ｌ、間隔Ｔを示す符号Ｃ_Ｔ、正規化係数列X_N[1]，…，X_N[N]を可変長符号化した可変長符号Ｃ_Ｘを出力する。また、必要に応じて値δを示す符号Ｃ_δと基準となる可変長符号化パラメータｓｂを示す符号Ｃ_ｓｂも出力する。符号化装置２００から出力された符号は、復号装置４００に入力される。

【0062】

［符号化装置の変形例１］（外部から情報が入力される例）
なお、符号化装置としては、周期性包絡系列生成部１４０と周期性統合包絡生成部２５０と可変長符号化パラメータ計算部２６０と可変長符号化部２７０だけを備え、符号化装置の外部で生成された平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と、正規化係数列X_N[1]，…，X_N[N]、間隔Ｔと、必要に応じて振幅スペクトル包絡系列W[1]，…， W[N]と、必要に応じて指標Ｓとを入力とし、可変長符号Ｃ_Ｘを出力してもよい。

【0063】

［符号化装置の変形例２]（係数列X[n]から間隔Ｔを求める例）
上述の周期性分析部２３０では正規化係数列X_N[1]，…，X_N[N]を入力として間隔Ｔを求めているが、周期性分析部２３０では周波数領域変換部１１０が出力した係数列X[1]，…，X[N]を入力として間隔Ｔを求めてもよい。この場合は、実施例１の周期性分析部１３０と同じ方法で間隔Ｔを求める。

【0064】

≪復号装置≫
図７に実施例２の復号装置の機能構成例を、図８に実施例２の復号装置の処理フローを示す。復号装置４００は、スペクトル包絡系列計算部４２１、周期性包絡系列生成部４４０、周期性統合包絡生成部４５０、可変長符号化パラメータ計算部４６０、可変長復号部４７０、周波数領域系列逆正規化部４１１、周波数領域逆変換部４１０を備える。復号装置４００は、量子化済線形予測係数^α₁，…，^α_Pを示す符号Ｃ_Ｌ、間隔Ｔを示す符号Ｃ_Ｔ、正規化係数列X_N[1]，…，X_N[N]を可変長符号化した可変長符号Ｃ_Ｘを受け取り、音響信号を出力する。なお、必要に応じて値δを示す符号Ｃ_δと基準となる可変長符号化パラメータｓｂを示す符号Ｃ_ｓｂと指標Ｓを示す符号Ｃ_Ｓも受け取る。以下に、各構成部の詳細を示す。

【0065】

＜スペクトル包絡系列計算部４２１＞
スペクトル包絡系列計算部４２１は、符号Ｃ_Ｌを入力とし、振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を計算する（Ｓ４２１）。より具体的には、以下の手順で処理すればよい。

【0066】

（ｓｔｅｐ１）符号Ｃ_Ｌを復号し、復号線形予測係数^α₁，…，^α_Pを得る。

【0067】

（ｓｔｅｐ２）復号線形予測係数^α₁，…，^α_Pを用いてＮ点の振幅スペクトル包絡系列W[1]，…，W[N]を求める。例えば、振幅スペクトル包絡系列の各値W[n]は、式（２）で求めることができる。

【0068】

（ｓｔｅｐ３）復号線形予測係数^α_ｐのそれぞれにγ^ｐを乗算し、復号平滑化線形予測係数^α₁γ，^α₂γ²，…，^α_Pγ^Pを求める。γはあらかじめ定めた平滑化するための１以下の正の定数である。そして、式（１０）によって、平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を求める。

【0069】

＜周期性包絡系列生成部４４０＞
周期性包絡系列生成部４４０は、間隔Ｔを示す符号Ｃ_Ｔを入力とし、符号Ｃ_Ｔを復号し、間隔Ｔを得る。そして、符号化装置２００の周期性包絡系列生成部１４０と同じ方法で周期性包絡系列P[1]，…，P[N]を求め、出力する（Ｓ４４０）。

【0070】

＜周期性統合包絡生成部４５０＞
周期性統合包絡生成部４５０には、周期性包絡系列P[1]，…，P[N]、振幅スペクトル包絡系列W[1]，…，W[N]、符号Ｃ_δ、符号Ｃ_Ｓが入力される。ただし、符号Ｃ_δ、符号Ｃ_Ｓは入力されない場合もある。周期性統合包絡生成部４５０は、符号Ｃ_δを復号し、値δを取得する。ただし、符号Ｃ_δが入力されない場合は、符号Ｃ_δの復号は行わず、周期性統合包絡生成部４５０に予め記憶された値δを取得する。なお、周期性統合包絡生成部４５０は、符号Ｃ_Ｓが入力された場合には、符号Ｃ_Ｓを復号して指標Ｓを取得し、取得した指標Ｓが、周期性が高いことに対応するフレームの場合には符号Ｃ_δを復号して値δを取得し、取得した指標Ｓが、周期性が低いことに対応するフレームである場合には符号Ｃ_δの復号は行わず、周期性統合包絡生成部４５０に予め記憶された値δを取得する。そして、周期性統合包絡生成部４５０は、式（６）によって、周期性統合包絡系列W_M[1]，…，W_M[N]を求める。（Ｓ４５０）

【0071】

＜可変長符号化パラメータ計算部４６０＞
可変長符号化パラメータ計算部４６０は、周期性統合包絡系列W_M[1]，…，W_M[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と符号Ｃ_ｓｂを入力とし、可変長符号化パラメータｒ_ｎを得る（Ｓ４６０）。ただし、復号装置４００に伝送される別の情報から振幅の平均値を推定できる場合は、別の情報から推定した振幅の平均値の推定値からｓｂを近似的に決定する方法を決めておいてもよい。この場合は、符号Ｃ_ｓｂは入力されない。以下に、１サンプルごとにライス復号を行う場合を例に、可変長符号化パラメータの計算方法を説明する。

【0072】

（ｓｔｅｐ１）符号Ｃ_ｓｂを復号して、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）を得る。なお、符号化装置２００と復号装置４００で共通に振幅の平均値の推定値からｓｂを近似的に決定する方法を決めている場合は、その方法で求める。

【0073】

（ｓｔｅｐ２）閾値θを式（１４）で算出する。

【0074】

（ｓｔｅｐ３） |W_M[n]/~W[n]|がθより大きいほど、ライスパラメータｒ_ｎをｓｂよりも大きな値として、符号化装置２００の可変長符号化パラメータ計算部２６０と同じ方法で決定する。|W_M[n]/~W[n]|がθより小さいほど、ライスパラメータｒ_ｎをｓｂよりも小さな値として、符号化装置２００の可変長符号化パラメータ計算部２６０と同じ方法で決定する。

【0075】

（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各X_N[n]についてのライスパラメータｒ_ｎを求める。

【0076】

＜可変長復号部４７０＞
可変長復号部４７０は、可変長符号化パラメータ計算部４６０で求めた可変長符号化パラメータｒ_ｎを用いて可変長符号Ｃ_Ｘを復号して復号正規化係数列^X_N[1]，…，^X_N[N]を得る（Ｓ４７０）。例えば、可変長復号部４７０は、可変長符号化パラメータ計算部４６０で求めたライスパラメータｒ_ｎを用いて可変長符号Ｃ_Ｘを復号して復号正規化係数列^X_N[1]，…，^X_N[N]を得る。可変長復号部４７０の復号方法は、可変長符号化部２７０の符号化方法に対応するものである。

【0077】

＜周波数領域系列逆正規化部４１１＞
周波数領域系列逆正規化部４１１は、復号正規化係数列^X_N[1]，…，^X_N[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]を入力とし、
^X[n]=^X_N[n]・~W[n] （１５）
のように、復号係数列^X[1]，…，^X[N]を求めて出力する（Ｓ４１１）。

【0078】

＜周波数領域逆変換部４１０＞
周波数領域逆変換部４１０は、復号係数列^X[1]，…，^X[N]を入力とし、復号係数列^X[1]，…，^X[N]を所定の時間区間であるフレーム単位の音響信号（時間領域）に変換する（Ｓ４１０）。

【0079】

［復号装置の変形例１］（外部から情報が入力される例）
なお、復号装置としては、周期性包絡系列生成部４４０と周期性統合包絡生成部４５０と可変長符号化パラメータ計算部４６０と可変長復号部４７０だけを備え、復号装置に必要に応じて入力される符号Ｃ_δと符号Ｃ_ｓｂに加えて、復号装置の外部で得られた平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]、振幅スペクトル包絡系列W[1]，…，W[N]、間隔Ｔ、必要に応じて指標Ｓも入力とし、正規化係数列X_N[1]，…，X_N[N]を出力し、外部で平滑化振幅スペクトル包絡系列を乗算して時間領域の音響信号に変換してもよい。

【0080】

＜実施例２の発明の効果＞
可変長符号化は、符号化対象の入力値の振幅の取りうる範囲に合わせて適応的に符号を決定することで符号化効率を向上させる符号化方法である。実施例２では周波数領域の係数列である正規化係数列X_N[1]，…，X_N[N]を符号化対象としているが、符号化対象の係数列に含まれる各係数の振幅の情報をより正確に用いて求めた可変長符号化パラメータを用いて可変長符号化をすれば符号化装置が行う可変長符号化自体の符号化効率は高くなる。しかし、復号装置が可変長符号化パラメータを求めるために、符号化装置から復号装置に対して符号化対象の係数列に含まれる各係数の振幅の情報をより正確に送る必要があり、その分だけ符号化装置から復号装置に送る符号量が増大してしまう。

【0081】

符号量の増大を抑えるためには、少ない符号量の符号から符号化対象の係数列に含まれる各係数の振幅の推定値を得る方法が必要である。実施例２の周期性統合包絡系列W_M[1]，…，W_M[N]は係数列X[1]，…，X[N]を高精度に近似するので、|W_M[1]/~W[1]|，…，|W_M[N]/~W[N]|は可変長符号化対象の係数であるX_N[1], X_N[2],…, X_N[N]の振幅包絡を高精度に近似できる。つまり、|W_M[1]/~W[1]|，…，|W_M[N]/~W[N]|は、符号化対象の各係数の振幅と正の相関を持つ系列となっている。

【0082】

また、|W_M[1]/~W[1]|, |W_M[2]/~W[2]|,…,|W_M[N]/~W[N]|を復号装置側で復元するために必要な情報は、
・量子化済線形予測係数^α₁，…，^α_Pの情報（符号Ｃ_Ｌ）
・間隔Ｔを示す情報（符号Ｃ_Ｔ）
・値δを示す情報（符号Ｃ_δ）
である。すなわち、実施例２の符号化装置と復号装置によれば、符号化装置に入力された入力音響信号のピッチ周期に起因する振幅のピークを含む包絡を、符号Ｃ_Ｌ、符号Ｃ_Ｔ、符号Ｃ_δのみの少ない情報量で、復号装置で再現することが可能となる。

【0083】

なお、実施例２の符号化装置と復号装置は、線形予測やピッチ予測を伴う符号化及び復号を行う符号化装置及び復号装置と併用して用いられることが多い。この場合は、符号Ｃ_Ｌと符号Ｃ_Ｔは、符号化装置２００外にある線形予測やピッチ予測を伴う符号化を行う符号化装置から、復号装置４００外にある線形予測やピッチ予測を伴う復号を行う復号装置に送られている符号である。したがって、符号化装置側に入力された入力音響信号のピッチ周期に起因する振幅のピークを含む包絡を復号装置側で復元するために符号化装置２００から復号装置４００に送る必要があるのは符号Ｃ_δである。符号Ｃ_δの符号量は小さく（それぞれ、せいぜい３ビット程度であり、１ビットでも効果が得られる）、符号化対象の正規化係数列に含まれる部分系列ごとの可変長符号化パラメータに対応する符号の総符号量よりも少ない。

【0084】

よって、実施例２の符号化装置、復号装置によれば、少ない符号量の増加で、符号化効率を向上させることができる。

【0085】

＜実施例２の発明のポイント＞
上述の効果を得るというポイントで実施例２の符号化装置、復号装置を考えると、符号化装置２００が、
・所定時間区間の入力音響信号から求めた線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、入力音響信号から求めた周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部２５０
・周期性統合包絡系列の値が大きい周波数ほど、入力音響信号の振幅が大きいことを前提に、入力音響信号に由来する周波数領域の系列を符号化する可変長符号化部２７０
を有し、復号装置４００が、
・線形予測係数符号に対応する周波数領域の系列であるスペクトル包絡系列と、周期符号に対応する周波数領域の周期と、に基づく周波数領域の系列である周期性統合包絡系列を生成する周期性統合包絡生成部４５０
・周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に、可変長符号を復号して周波数領域の系列を得る可変長復号部４７０、
を有することを特徴とすればよい。なお、「周期性統合包絡系列の値が大きい周波数ほど、入力音響信号の振幅が大きいことを前提に」と「周期性統合包絡系列の値が大きい周波数ほど、音響信号の振幅が大きいことを前提に」とは、周期性統合包絡系列が、入力音響信号または音響信号の振幅の大きい周波数において大きい値になることを特徴としていることを示している。また、「入力音響信号に由来する」とは、入力音響信号から求められることや入力音響信号に対応していることを意味している。例えば、係数列X[1]，…，X[N]や正規化係数列X_N[1]，…，X_N[N]は、入力音響信号に由来する周波数領域の系列である。

【実施例3】

【0086】

≪符号化装置≫
図９に実施例３の符号化装置の機能構成例を、図１０に実施例３の符号化装置の処理フローを示す。符号化装置３００は、スペクトル包絡系列計算部２２１、周波数領域変換部１１０、周波数領域系列正規化部１１１、周期性分析部３３０、周期性包絡系列生成部１４０、周期性統合包絡生成部２５０、可変長符号化パラメータ計算部２６０、第２可変長符号化パラメータ計算部３８０、可変長符号化部３７０を備える。符号化装置３００は、入力された時間領域の音響ディジタル信号を入力音響信号x(t)とし、少なくとも量子化済線形予測係数^α₁，…，^α_Pを示す符号Ｃ_Ｌ、正規化係数列X_N[1]，…，X_N[N]の周期を表す間隔Ｔの符号Ｃ_Ｔ、係数列X[1]，…，X[N]または正規化係数列X_N[1]，…，X_N[N]の周期性の程度を示す所定の指標Ｓと指標Ｓを示す符号Ｃ_Ｓ、正規化係数列X_N[1]，…，X_N[N]を可変長符号化した可変長符号Ｃ_Ｘを出力する。周波数領域系列正規化部１１１は実施例１変形例１と同じである。周波数領域変換部１１０と周期性包絡系列生成部１４０は実施例１と同じである。振幅スペクトル包絡系列計算部２２１、周期性統合包絡生成部２５０、可変長符号化パラメータ計算部２６０は、実施例２と同じである。以下では異なる構成部について説明する。

【0087】

＜周期性分析部３３０＞
周期性分析部３３０は、正規化係数列X_N[1]，…，X_N[N]を入力とし、
当該正規化係数列X_N[1]，…，X_N[N]の周期性の程度を示す指標Ｓと間隔Ｔ（周期的に大きな値となる間隔）とを求め、指標Ｓと指標Ｓを示す符号Ｃ_Ｓと間隔Ｔと間隔Ｔを示す符号Ｃ_Ｔを出力する（Ｓ３３０）。なお、指標Ｓと間隔Ｔ自体は実施例１変形例１の周期性分析部１３１と同じである。

【0088】

そして、符号化装置３００では、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部２６０が可変長符号化パラメータｒ_ｎを計算し、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部３８０が可変長符号化パラメータｒ_ｎを計算する（Ｓ３９０）。「あらかじめ定めた周期性の程度が大きいことを示す範囲」は、例えば、指標Ｓが所定の閾値以上のときとすればよい。

【0089】

＜第２可変長符号化パラメータ計算部３８０＞
第２可変長符号化パラメータ計算部３８０は、振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と正規化係数列X_N[1]，…，X_N[N]を入力とし、可変長符号化パラメータｒ_ｎを求める（Ｓ３８０）。可変長符号化パラメータ計算部２６０は、周期性統合包絡系列W_M[1]，…，W_M[N]から求めた振幅値に依存して可変長符号化パラメータｒ_ｎを計算することを特徴としているのに対して、第２可変長符号化パラメータ計算部３８０は、振幅スペクトル包絡系列から求めた振幅値に依存して可変長符号化パラメータを計算することを特徴としている。以下に、１サンプルごとにライス符号化を行う場合を例に、可変長符号化パラメータの計算方法を説明する。

【0090】

（ｓｔｅｐ１）正規化係数列X_N[1],…,X_N[N]の各係数の振幅の平均の対数を、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）として式（１３）のように算出する。この処理は、可変長符号化パラメータ計算部２６０と同じである。

【0091】

（ｓｔｅｐ２）下記式により閾値θを算出する。

【数11】

θは、振幅スペクトル包絡系列の各値W[n]を平滑化振幅スペクトル包絡系列の各値~W[n]で除算した値の振幅の平均の対数である。

【0092】

（ｓｔｅｐ３） |W[n]/~W[n]|がθより大きいほど、正規化係数X_N[n]をライス符号化するためのライスパラメータｒ_ｎをｓｂよりも大きな値として決定する。|W[n]/~W[n]|がθより小さいほど、正規化係数X_N[n]をライス符号化するためのライスパラメータｒ_ｎをｓｂよりも小さな値として決定する。

【0093】

（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各X_N[n]についてのライスパラメータｒ_ｎを求める。

【0094】

＜可変長符号化部３７０＞
可変長符号化部３７０は、可変長符号化パラメータｒ_ｎを用いて正規化係数列X_N[1],…,X_N[N]を可変長符号化し、可変長符号Ｃ_Ｘを出力する（Ｓ３７０）。ただし、可変長符号化パラメータｒ_ｎは、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部２６０が計算した可変長符号化パラメータｒ_ｎであり、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部３８０が計算した可変長符号化パラメータｒ_ｎである。

【0095】

符号化装置３００は、このような処理によって得られた量子化済線形予測係数^α₁，…，^α_Pを示す符号Ｃ_Ｌ、周期性の程度を示す指標Ｓを示す符号Ｃ_Ｓ、間隔Ｔを示す符号Ｃ_Ｔ、正規化係数列X_N[1],…,X_N[N]を可変長符号化した可変長符号Ｃ_Ｘを出力し、復号側に送信する。また、必要に応じて値δを示す符号Ｃ_δと基準となる可変長符号化パラメータｓｂを示す符号Ｃ_ｓｂも出力し、復号側に送信する。

【0096】

［符号化装置の変形例１］（外部から情報が入力される例）
なお、符号化装置としては、周期性包絡系列生成部１４０と周期性統合包絡生成部２５０と可変長符号化パラメータ計算部２６０と第２可変長符号化パラメータ計算部３８０と可変長符号化部３７０だけを備え、符号化装置の外部で生成された平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と正規化係数列X_N[1],…,X_N[N]、間隔Ｔと、必要に応じて振幅スペクトル包絡系列W[1]，…， W[N]と、必要に応じて指標Ｓとを入力とし、可変長符号Ｃ_Ｘを出力してもよい。

【0097】

［符号化装置の変形例２]（係数列X[n]から間隔Ｔを求める例）
上述の周期性分析部３３０では正規化係数列X_N[1]，…，X_N[N]を入力として間隔Ｔを求めているが、周期性分析部３３０では周波数領域変換部１１０が出力した係数列X [1]，…，X [N]を入力として間隔Ｔを求めてもよい。この場合は、実施例１の周期性分析部１３０と同じ方法で間隔Ｔを求める。

【0098】

≪復号装置≫
図１１に実施例３の復号装置の機能構成例を、図１２に実施例３の復号装置の処理フローを示す。復号装置５００は、スペクトル包絡系列計算部４２１、指標復号部５３０、周期性包絡系列生成部４４０、周期性統合包絡生成部４５０、可変長符号化パラメータ計算部４６０、第２可変長符号化パラメータ計算部５８０、可変長復号部５７０、周波数領域系列逆正規化部４１１、周波数領域逆変換部４１０を備える。復号装置５００は、量子化済線形予測係数^α₁，…，^α_Pを示す符号Ｃ_Ｌ、指標Ｓを示す符号Ｃ_Ｓ、間隔Ｔを示す符号Ｃ_Ｔ、正規化係数列X_N[1]，…，X_N[N]を可変長符号化した可変長符号Ｃ_Ｘを受け取り、音響信号を出力する。なお、必要に応じて値δを示す符号Ｃ_δと基準となる可変長符号化パラメータｓｂを示す符号Ｃ_ｓｂも受け取る。スペクトル包絡系列計算部４２１、周期性包絡系列生成部４４０、周期性統合包絡生成部４５０、可変長符号化パラメータ計算部４６０、周波数領域系列逆正規化部４１１、周波数領域逆変換部４１０は実施例２と同じである。以下では異なる構成部について説明する。

【0099】

＜指標復号部５３０＞
指標復号部５３０は、符号Ｃ_Ｓを復号し、指標Ｓを得る。復号装置５００では、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部４６０が可変長符号化パラメータｒ_ｎを計算し、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部５８０が可変長符号化パラメータｒ_ｎを計算する（Ｓ５９０）。なお、「あらかじめ定めた周期性の程度が大きいことを示す範囲」は、符号化装置３００と同じ範囲である。

【0100】

＜第２可変長符号化パラメータ計算部５８０＞
第２可変長符号化パラメータ計算部５８０は、振幅スペクトル包絡系列W[1]，…，W[N]と平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]と符号Ｃ_ｓｂを入力とし、可変長符号化パラメータｒ_ｎを求める（Ｓ５８０）。ただし、復号装置５００に伝送される別の情報から振幅の平均値を推定できる場合は、別の情報から推定した振幅の平均値の推定値からｓｂを近似的に決定する方法を決めておいてもよい。この場合は、符号Ｃ_ｓｂは入力されない。以下に、１サンプルごとにライス復号を行う場合を例に、可変長符号化パラメータの計算方法を説明する。

【0101】

（ｓｔｅｐ１）符号Ｃ_ｓｂを復号して、基準となるライスパラメータｓｂ（基準となる可変長符号化パラメータ）を得る。なお、符号化装置３００と復号装置５００で共通に振幅の推定値からｓｂを近似的に決定する方法を決めている場合は、その方法で求める。

【0102】

（ｓｔｅｐ２）閾値θを式（１６）で算出する。

【0103】

（ｓｔｅｐ３） |W[n]/~W[n]|がθより大きいほど、ライスパラメータｒ_ｎをｓｂよりも大きな値として、符号化装置３００の第２可変長符号化パラメータ計算部３８０と同じ方法で決定する。|W[n]/~W[n]|がθより小さいほど、ライスパラメータｒ_ｎをｓｂよりも小さな値として、符号化装置３００の第２可変長符号化パラメータ計算部３８０と同じ方法で決定する。

【0104】

（ｓｔｅｐ４）ｓｔｅｐ３の処理を全てのn=1,2,…,Nについて繰り返して、各X_N[n]についてのライスパラメータｒ_ｎを求める。

【0105】

＜可変長復号部５７０＞
可変長復号部５７０は、可変長符号化パラメータｒ_ｎを用いて可変長符号Ｃ_Ｘを復号して復号正規化係数列^X_N[1]，…，^X_N[N]を求める（Ｓ５７０）。ただし、可変長符号化パラメータｒ_ｎは、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲の場合は、可変長符号化パラメータ計算部４６０が計算した可変長符号化パラメータｒ_ｎであり、指標Ｓがあらかじめ定めた周期性の程度が大きいことを示す範囲ではない場合は、第２可変長符号化パラメータ計算部５８０が計算した可変長符号化パラメータｒ_ｎである。

【0106】

［復号装置の変形例１］（外部から情報が入力される例）
なお、復号装置としては、周期性包絡系列生成部４４０と周期性統合包絡生成部４５０と可変長符号化パラメータ計算部４６０と第２可変長符号化パラメータ計算部５８０と可変長復号部５７０だけを備え、復号装置に必要に応じて入力される符号Ｃ_δと符号Ｃ_ｓｂに加えて、復号装置の外部で得られた平滑化振幅スペクトル包絡系列~W[1]，…，~W[N]、振幅スペクトル包絡系列W[1]，…，W[N]、間隔Ｔ、指標Ｓも入力とし、正規化係数列X_N[1]，…，X_N[N]を出力し、外部で平滑化振幅スペクトル包絡系列を乗算して時間領域の音響信号に変換してもよい。

【0107】

＜実施例３の発明の効果＞
入力音響信号の周期性の程度が小さい場合には、入力音響信号のピッチ周期に起因する振幅のピークは小さい。そこで、実施例３の符号化装置、復号装置は、符号化の対象となる音響信号の周期性の程度が大きい場合には周期性統合包絡系列を用いて可変長符号化パラメータを求め、符号化の対象となる音響信号の周期性の程度が大きくない場合には振幅スペクトル包絡系列を用いて可変長符号化パラメータを求めるため、より適した可変長符号化パラメータを用いて可変長符号化でき、符号化精度を上げることができるという効果がある。

【0108】

上述の実施例１〜３では、振幅スペクトル包絡系列、平滑化振幅スペクトル包絡系列、周期性統合包絡系列等について振幅の系列を用いる例を説明したが、振幅の系列に代えてパワーの系列、すなわち、W[n]、~W[n]、W_M[n]としてパワースペクトル包絡系列、平滑化パワースペクトル包絡系列、パワーの系列である周期性統合包絡系列を用いてもよい。

【0109】

［プログラム、記録媒体］
上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。

【0110】

また、上述の構成をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。

【0111】

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。

【0112】

また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。

【0113】

このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。

【0114】

また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。

【符号の説明】

【0115】

１００、１０１ 周期性統合包絡系列生成装置
１１０ 周波数領域変換部 １１１ 周波数領域系列正規化部
１２０、１２１、２２１、４２１ スペクトル包絡系列計算部
１３０，１３１、２３０、３３０ 周期性分析部
１４０、４４０ 周期性包絡系列生成部
１５０、２５０、４５０ 周期性統合包絡生成部
２００、３００ 符号化装置
２６０、３６０、４６０ 可変長符号化パラメータ計算部
２７０、３７０ 可変長符号化部
３８０、５８０ 第２可変長符号化パラメータ計算部
４００、５００ 復号装置
４１０ 周波数領域逆変換部 ４１１ 周波数領域系列逆正規化部
４７０、５７０ 可変長復号部 ５３０ 指標復号部

【図1】

【図2】

【図3】

【図4A】

【図4B】

【図4C】

【図4D】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】