JP4848049B2 - 符号化方法、復号方法、それらの装置、プログラム及び記録媒体 - Google Patents

符号化方法、復号方法、それらの装置、プログラム及び記録媒体 Download PDF

Info

Publication number
JP4848049B2
JP4848049B2 JP2010542108A JP2010542108A JP4848049B2 JP 4848049 B2 JP4848049 B2 JP 4848049B2 JP 2010542108 A JP2010542108 A JP 2010542108A JP 2010542108 A JP2010542108 A JP 2010542108A JP 4848049 B2 JP4848049 B2 JP 4848049B2
Authority
JP
Japan
Prior art keywords
prediction
prediction order
order
coefficient
code
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2010542108A
Other languages
English (en)
Other versions
JPWO2010067800A1 (ja
Inventor
健弘 守谷
登 原田
優 鎌本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2010542108A priority Critical patent/JP4848049B2/ja
Application granted granted Critical
Publication of JP4848049B2 publication Critical patent/JP4848049B2/ja
Publication of JPWO2010067800A1 publication Critical patent/JPWO2010067800A1/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
    • H03M7/40—Conversion to or from variable length codes, e.g. Shannon-Fano code, Huffman code, Morse code
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction

Description

本発明は、時系列信号を線形予測分析して符号化する技術に関し、特に、予測次数を適応的に選択して線形予測符号化を行う符号化方法、その復号方法、それらの装置、プログラム及び記録媒体に関する。
音響信号や映像情報などの時系列信号を通信路によって伝送したり、情報記録媒体に記録する場合、時系列信号を圧縮符号に変換してから伝送したり記録する方法が、伝送効率や記録効率の点で有効である。また、近年のブロードバンドの普及や記憶装置の容量増加に伴い、圧縮率の高さを優先する非可逆圧縮符号化（ロッシー符号化：Lossy Coding）方式よりも原信号の完全再生を条件とした可逆圧縮符号化（ロスレス符号化：Lossless Coding）方式が重視されつつある（例えば、非特許文献１参照）。そのような中、線形予測分析技術を用いて音響信号を可逆圧縮符号化する線形予測符号化方式がMPEG（Moving Picture Expert Group）の国際標準規格「MPEG-4 ALS」として承認されている（例えば、非特許文献２参照）。
図１は、線形予測符号化を行う従来の符号化装置１１００の機能構成を説明するためのブロック図である。図２は、図１の符号化装置１１００で生成された符号を復号する復号装置１２００の機能構成を説明するためのブロック図である。また、図３は、可逆圧縮符号化方式における予測次数と符号量（「符号長」と呼ぶ場合もある）の関係を説明するためのグラフである。なお、図３の横軸は予測次数を示し、縦軸は符号量を示す。
＜符号化方法＞
符号化装置１１００のフレームバッファ１１１１には、標本化・量子化されたPCM (pulse code modulation)形式の時系列信号x(n)（nは離散時間を示すインデックス）が入力される。フレームバッファ１１１１は、予め定められた時間区間（以下、「フレーム」と呼ぶ）分の時系列信号x(n)（n=1,...,N）（Ｎは正の整数）をバッファし、符号化装置１１００は、フレーム毎に時系列信号x(n)を符号化する。
まず、１フレーム分の時系列信号x(n)（n=1,...,N）が線形予測分析部１１１２に送られ、線形予測分析部１１１２は、線形予測分析によって１次からPopt次までのPARCOR係数k(m)（m=1,...,Popt）を算出する。
［線形予測分析と最適予測次数］
線形予測分析では、ある時点nの時系列信号x(n)と、その時点nよりも過去のP個（Pは正の整数である。Pを「予測次数」と呼ぶ。）の時点n-1,n-2,...,n-Pの時系列信号x(n-1),x(n-2),...,x(n-P)をそれぞれ係数α(m)(m=1,...,P)（「線形予測係数」と呼ぶ）で重み付けしたものと、予測残差e(n)（「予測誤差」と呼ぶ場合もある）と間に線形1次結合が成り立つと仮定する。この仮定に基づいた線形予測モデルは式(1)のようになる。線形予測分析では、入力された時系列信号x(n)（n=1,...,N）に対し、式(1)を満たす線形予測係数α(m)（m=1,...,P）又はそれに変換可能なPARCOR係数k(m)（m=1,...,P）などの係数を算出する（予測次数Pについての線形予測分析）。
e(n)=x(n)+α(1)・x(n-1)+α(2)・x(n-2)+...+α(P)・x(n-P) …(1)
また、ある時点nの時系列信号y(n)を、その時点nよりも過去のP個の時点n-1,n-2,...,n-Pの時系列信号x(n-1),x(n-2),...,x(n-P)を用いて推定する式(2)の線形FIR(Finite Impulse Response)フィルタを「線形予測フィルタ」と呼ぶ。
y(n)=-{α(1)・x(n-1)+α(2)・x(n-2)+...+α(P)・x(n-P)} …(2)
また、Poptは最適な予測次数Pを示す正の整数であり、これを「最適予測次数」と呼ぶ。非特許文献２の方式では、ＭＤＬ原理(Minimum Description Length Principle)に基づいて、最適予測次数Poptを決定する。ＭＤＬ原理では、
（符号語長）＝（モデルの記述長）＋（そのモデルによるデータの記述長）
を最小にするモデルを最良とする。すなわち、非特許文献２の方式では、
（ロスレスで復号するために必要な符号量）＝（PARCOR係数のために必要な符号量）＋（予測残差のために必要な符号量） …(3)
を最小にする予測次数Pを最適予測次数Poptとする。
図３の直線３Ａで模式的に示すようにPARCOR係数のために必要な符号量は予測次数に比例して増加する。また、一般に予測次数が高くなれば予測残差のエネルギーは小さくなり、予測残差をエントロピー符号化した場合の符号量は、曲線３Ｂで模式的に示すように対数的に小さくなる。そのため、直線３Ａと曲線３Ｂの和である曲線３Ｃで模式的に示すように、ロスレスで復号するために必要な符号量は、予測次数が高くなるほど小さくなるのではなく、或る予測次数で最小となる。線形予測分析部１１１２は、予め定められた最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数を探索範囲とし、ロスレス復号するために必要な符号量が最小となる予測次数を最適予測次数Poptとして選択する（［線形予測分析と最適予測次数］の説明終わり）。
算出された１次からPopt次までのPARCOR係数k(m)（m=1,...,Popt）は、量子化部１１１３に送られ、量子化されて１次からPopt次までの量子化PARCOR係数i(m)（m=1,...,Popt）が生成される。なお、「量子化PARCOR係数」は、PARCOR係数の量子化値そのものであってもよいし、当該PARCOR係数の量子化値に付されたインデックスであってもよい。１次からPopt次までの量子化PARCOR係数i(m)（m=1,...,Popt）は、係数符号化部１１１４に送られ、そこでエントロピー符号化されて係数符号Ckが生成される。また、１次からPopt次までの量子化PARCOR係数i(m)（m=1,...,Popt）は、線形予測係数変換部１１１５にも送られる。線形予測係数変換部１１１５は、これらのPARCOR係数i(m)（m=1,...,Popt）と線形予測分析部１１１２で選択された最適予測次数Poptとを用い、最適予測次数Poptの線形予測フィルタの各線形予測係数α(m)（m=1,...,Popt）を算出する。
線形予測部１１１６は、１フレーム分の時系列信号x(n)（n=1,...,N）と各線形予測係数α(m)（m=1,...,Popt）と線形予測分析部１１１２で選択された最適予測次数Poptとを用い、最適予測次数Poptの線形予測フィルタによって線形予測値y(n)（n=1,...,N）を生成し、減算部１１１７は、時系列信号x(n)から線形予測値y(n)を減算した予測残差e(n)を算出する（予測フィルタ処理）。算出された予測残差e(n)（n=1,...,N）は、残差符号化部１１２０に送られ、そこでエントロピー符号化されて残差符号Ceが生成される。
線形予測分析部１１１２で選択された最適予測次数Poptと、係数符号化部１１１４で生成された係数符号Ckと、残差符号化部１１２０で生成された残差符号Ceは、合成部１１３０に送られ、そこで合成されて符号Cgが生成される。
＜復号方法＞
復号装置１２００に入力された符号Cgは、分離部１２１０で、最適予測次数Poptと係数符号Ckと残差符号Ceとに分離される。係数符号Ckと残差符号Ceとは、それぞれ係数復号部１２３１と残差復号部１２２０で復号され、１次からPopt次までの量子化PARCOR係数i(m)（m=1,...,Popt）と、予測残差e(n)（n=1,...,N）とが生成される。
最適予測次数Poptと量子化PARCOR係数i(m)（m=1,...,Popt）は、線形予測係数変換部１２３２に送られ、線形予測係数変換部１２３２は、これらを用いて最適予測次数Poptの線形予測フィルタの各線形予測係数α(m)（m=1,...,Popt）を算出する。線形予測部１２３３は、１次からPopt次までの各線形予測係数α(m)（m=1,...,Popt）をフィルタ係数とする、最適予測次数Poptの線形予測フィルタによって線形予測値y(n)（n=1,...,N）を生成し、加算部１２３４は、線形予測値y(n)と予測残差e(n)とをそれぞれ加算して時系列信号のロスレス復元値x(n)（n=1,...,N）を生成する（逆予測フィルタ処理）。
MatHans, "Lossless Compression of Digital Audio", IEEE SIGNAL PROCESSING MAGAZINE, July 2001, pp.21-32. ISO/IEC 14496-3 AMENDMENT 2: Audio Lossless Cording (ALS), new audio profiles and BSAC extensions.
従来方式では、復号装置が最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数を探索範囲とし、最適予測次数Poptを選択していた。
この場合、最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数が最適予測次数Poptの候補となるため、探索に必要な演算量が増加する。
また、符号化装置で選択された最適予測次数Poptは復号装置に送られるが、最適予測次数Poptとなり得る整数の候補が多くなるにつれて、それらを指定するための符号の平均符号長も長くなり、符号圧縮率が低下する。
本発明はこのような点に鑑みてなされたものであり、最適予測次数の探索に必要な演算量を削減し、符号圧縮率を向上させることが可能な技術を提供することを目的とする。
本発明では、最小予測次数を予め定められた０以上の整数とし、最大予測次数を当該最小予測次数よりも大きな予め定められた整数とした場合における、当該最小予測次数以上かつ当該最大予測次数以下のすべての整数を要素とする整数集合から一部の要素を間引いた部分集合の各要素から、入力された時系列信号を線形予測符号化するための最適予測次数を選択し、当該最適予測次数を指定するための符号を出力する。
これにより、最小予測次数以上かつ最大予測次数以下のすべての整数を最適予測次数の候補とする場合に比べ、最適予測次数の探索に必要な演算量を削減できる。また、最適予測次数となり得る整数の候補を間引くことにより、最適予測次数を指定するための符号の平均符号長も短くなり、符号圧縮率が向上する。
本発明では、最適予測次数の探索に必要な演算量を削減し、符号圧縮率を向上させることができる。
線形予測符号化を行う従来の符号化装置の機能構成を説明するためのブロック図。 図１の符号化装置で生成された符号を復号する復号装置の機能構成を説明するためのブロック図。 可逆圧縮符号化方式における予測次数と符号量の関係を説明するためのグラフ。 図４（Ａ）は、各予測次数で線形予測符号化を行った場合の圧縮率を説明するためのグラフであり、図４（Ｂ）は、最小予測次数Pmin=1及び最大予測次数Pmax=12とし、整数集合Sの全要素を探索範囲として適応的に最適予測次数Poptを選択した場合における、最適予測次数Poptの頻度分布を説明するためのグラフである。 第１実施形態の符号化装置の機能構成を説明するためのブロック図。 図５の線形予測分析部の機能構成を説明するためのブロック図。 第１実施形態の復号装置の機能構成を説明するためのブロック図。 第１実施形態の符号化方法を説明するためのフローチャート。 第１実施形態の復号方法を説明するためのフローチャート。 図１０（Ａ）（Ｂ）は可変長符号化のための符号表の例示である。 第２実施形態の線形予測分析部の機能構成を説明するためのブロック図。 第２実施形態の符号化方法の一部を説明するためのフローチャート。
以下では本発明の原理を説明した後、各実施形態の説明を行っていく。
〔原理〕
図４（Ａ）は、各予測次数で線形予測符号化を行った場合の圧縮率を説明するためのグラフである。ここで、図４（Ａ）の横軸は予測次数を示し、縦軸は各予測次数で線形予測符号化を行った場合の平均圧縮率を示す。なお、圧縮率とは、（線形予測符号のビット長／符号化対象のビット長）×１００を意味する。
図４（Ａ）に示すように、任意の２つの予測次数が互いに近ければ、それらにそれぞれ対応する平均圧縮率も互いに近い。特に、平均圧縮率を最小にする予測次数の近傍では、予測次数の変化に伴う平均圧縮率の変動が小さい。すなわち、最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数を探索範囲として選択された最適予測次数を用いなくても、圧縮性能は大きく劣化しない。本形態では、最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数を探索範囲として最適予測次数Poptを選択するのではなく、最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数Pmin, Pmin+1,..., Pmax-1, Pmaxを要素とする整数集合Sから一部の要素を間引いた部分集合Ssub⊂Sの各要素から、入力された時系列信号を線形予測符号化するための最適予測次数を選択し、当該最適予測次数を指定するための符号を生成する。これにより、最適予測次数Poptの探索に必要な演算量を削減できるとともに、最適予測次数Poptを指定するための符号の符号長も短くなり、符号圧縮率が向上する。なお、最小予測次数Pminは０以上の整数であり、最大予測次数Pmaxは最小予測次数Pmin以上の整数であり、何れも予め定められた値である。例えば、最小予測次数Pminや最大予測次数Pmaxは、処理単位となる時間区間（フレーム又はフレームを複数の時間区分に区分したサブフレーム）の長さごとに予め定められた値である。また、最適予測次数を指定するための符号の例は、各最適予測次数に付された最適予測次数ごとに異なるインデックスや、各最適予測次数を特定する情報を可変長符号化した符号である。また、可変長符号化される各最適予測次数を特定する情報は、各最適予測次数に付された最適予測次数ごとに異なるインデックスであってもよいし、各最適予測次数そのものであってもよい。また、可変長符号化の例は、アルファ符号化、デルタ符号化、ハフマン符号化（Huffman Coding）、ゴロム・ライス符号化（Golomb-Rice Coding）、ゴロム符号化（Golomb Coding）、及びその他のエントロピー符号化である。本形態では最適予測次数となり得る範囲が限定されるため、いずれの場合も各最適予測次数Poptを区別して表現するために必要な符号の平均符号量を削減できる。
＜間引き方の例＞
以下に、最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数Pmin, Pmin+1,..., Pmax-1, Pmaxを要素とする整数集合S=(Pmin, Pmin+1,..., Pmax-1, Pmax)から一部の要素を間引いた部分集合Ssubの設定例を示す。なお、floor{・}は・以下の最大の整数を示す。
(例1)整数集合Sの全体にわたって等間隔で間引きを行う。これは、本来の最適予測次数が整数集合S全体に一様に分布する場合に有効な設定方法である。
・Ssub=(Pmax，Pmax-2, Pmax-4,..., Pmin)
・Ssub=(Pmax, floor{(Pmax-Pmin)/2})
(例2)整数集合Sのうち所定値以上の整数を等間隔で間引き、当該所定値未満の整数は部分集合Ssubの要素としない。これは、本来の最適予測次数がPmaxに近い領域に多く分布する場合に有効な設定方法である。
・Ssub=(Pmax，Pmax-1, Pmax-2)
・Ssub=(Pmax，Pmax-2, Pmax-4,..., floor{(Pmax-Pmin)/2})
(例3)整数集合Sのうち所定値以上の整数を等間隔で間引き、当該所定値未満の整数はPmin近傍を除いて部分集合Ssubの要素としない。これは、本来の最適予測次数がPmaxに近い領域とPminに多く分布する場合に有効な設定方法である。
・Ssub=(Pmax，Pmax-1, Pmax-2, Pmin)
・Ssub=(Pmax，Pmax-2, Pmax-4,..., floor{Pmax/2}, Pmin)
・Ssub=(Pmax，Pmax-1, Pmax-2, Pmin+1,Pmin)
・Ssub=(Pmax-1, Pmax-2, Pmin+1)
(例4)Pmaxに近い領域ほど予め定められた整数区間中の要素数が多くなるように、整数集合Sの要素から不等間隔で要素を間引いて部分集合Ssubを設定する。これは、本来の最適予測次数がPmaxに近い領域に多く分布する場合に有効な設定方法である。
・Ssub=(Pmax，Pmax-1, Pmax-3, Pmax-6, Pmax-9,..., Pmin)
・Ssub=(Pmax，Pmax-1, Pmax-3, Pmax-6, Pmax-9)
・Ssub=(Pmax，Pmax-1, Pmax-3, Pmax-6, Pmax-9, Pmin)
＜最適予測次数Poptの選択方法の例＞
本形態の最適予測次数Poptの選択は、例えば、(I) 入力された時系列信号から線形予測係数若しくは当該線形予測係数に変換可能な係数を算出するステップと、(II) 予測次数を指定するための符号の符号量と、当該予測次数での線形予測係数若しくは当該線形予測係数に変換可能な係数の符号量又は推定符号量と、当該線形予測係数を用いて時系列信号の線形予測を行った場合の予測残差の符号量又は推定符号量と、の合計値SUMを指標とし、各予測次数から最適予測次数Poptを選択するステップと、を含む処理によって行なわれる。このように、予測次数を指定するための符号の符号量をも考慮して最適予測次数Poptを選択することで、特にフレーム中のサンプル数Ｎがさほど大きくない場合に、圧縮率をより向上させることが可能な最適予測次数Poptを選択することができる。なお、ステップ(II)の処理の例は以下の通りである。
(例11)部分集合Ssubの全要素を探索範囲とし、合計値SUMが最小となる予測次数を最適予測次数Poptとして選択する。
(例12)部分集合Ssubの要素を予め定められた順序（例えば昇降順）で探索し、合計値SUMが最初に予め定められた閾値以下となった予測次数を最適予測次数Poptとして選択する。
(例13)最初の処理単位となる時間区間（フレーム又はフレームを複数の時間区分に区分したサブフレーム）での最適予測次数Poptの選択は上記(11)又は(12)によって行う。それ以後の処理単位となる時間区間では、１つ前の時間区間で選択された最適予測次数よりも大きな最小の部分集合Ssubの要素（存在しない場合はPmax）を予測次数として選択し、新たに選択した予測次数に対応する新たな合計値SUMを求める。そして、新たな合計値SUMが１つ前の時間区間で求めた合計値SUMよりも小さい場合には、新たに選択した予測次数を新たな最適予測次数Poptとし、新たな合計値SUMが１つ前の時間区間で求めた合計値SUM以上の場合には、１つ前の時間区間で選択された最適予測次数を新たな最適予測次数Poptとする。
(例14)なお、予測次数Pでの線形予測係数若しくは線形予測係数に変換可能な係数の符号量又は推定符号量と、当該線形予測係数を用いて時系列信号の線形予測を行った場合の予測残差の符号量又は推定符号量との合計値を、当該予測次数Pに対応するSUMとして上記の選択を行ってもよい。
＜複数種類設定された予備集合Spreから部分集合Ssubを選択＞
また、整数集合Sから一部の要素を間引いた予備集合Spre⊂Sを複数種類設定しておき、それらの予備集合Spreの１つを部分集合Ssubとして選択してもよい。なお、予備集合Spreの具体例は、上述の部分集合Ssubの例と同様である。
また、予備集合Spreから部分集合Ssubを選択する基準としては、入力された時系列信号が予め定められた閾値より大きいか否かに応じて、予備集合Spreから部分集合Ssubを選択するものを例示できる。その具体例は以下の通りである。
(例21)入力された時系列信号が予め定められた閾値を超える場合、時系列信号が無音区間を含む可能性が低いと判断して(例2)のようにPminを要素に含まない予備集合Spreの一つを部分集合Ssubと選択し、入力された時系列信号が予め定められた閾値以下である場合、時系列信号が無音区間を含む可能性が高いと判断して(例3)のようにPminを要素に含む予備集合Spreの一つを部分集合Ssubと選択する。
(例22)入力された時系列信号が予め定められた閾値を超える場合、(例4)のようにPmaxに近い領域ほど予め定められた整数区間中の要素数が多くなるように、整数集合Sの要素から不当間隔で要素を間引いて生成された予備集合Spreを部分集合Ssubとして選択し、入力された時系列信号が予め定められた閾値以下である場合、(例3)のようにPminを要素に含む予備集合Spreの一つを部分集合Ssubと選択する。
＜最適予測次数Poptを特定する情報の可変長符号化＞
図４（Ｂ）は、最小予測次数Pmin=1及び最大予測次数Pmax=12とし、整数集合Sの全要素を探索範囲として適応的に最適予測次数Poptを選択した場合における、最適予測次数Poptの頻度分布を説明するためのグラフである。ここで、図４（Ｂ）の横軸は選択された最適予測次数Poptを示し、縦軸はその頻度を示す。なお、図４（Ｂ）の頻度は、それらの合計が１になるように正規化されている。
図４（Ｂ）の例では、最小予測次数Pmin=1及び最大予測次数Pmax=12とが最適予測次数Poptに選択される頻度が高く、それについで最大予測次数Pmax=12に近い予測次数が最適予測次数Poptに選択される頻度が高い。音響信号などの時系列信号の場合、適応的に選択される最適予測次数Poptが整数集合S全体にわたって一様に分布することは稀であり、選択される最適予測次数Poptの頻度分布に偏りが生じることが多い。特に、最大予測次数Pmaxを小さめに設定した場合には、選択される最適予測次数Poptが最大予測次数Pmax近傍に分布する頻度が高く、さらに時系列信号に無音区間が多く含まれている場合には、最小予測次数Pminの頻度も高くなる。
そのため、第１符号化対象が第２符号化対象よりも予め定められた特定の値に近い場合に、当該第１符号化対象に当該第２符号化対象の符号よりも短い符号長の符号が割り当てられる頻度が、当該第１符号化対象に当該第２符号化対象の符号よりも長い符号長の符号が割り当てられる頻度よりも高い可変長符号化方法を用い、選択された最適予測次数Poptを特定する情報を可変長符号化し、それを、最適予測次数Poptを指定するための符号とする。これにより、最適予測次数Poptを指定するための符号の符号量を削減できる。ここで、予め定められた特定の値は、例えば、最大予測次数Pmaxを特定する情報である。なお、可変長符号化される各最適予測次数Poptを特定する情報の例は、最適予測次数Poptに付されたインデックスや、最適予測次数Poptそのものである。また、このような可変長符号化の例外処理として、最適予測次数Poptが最小予測次数Pminである場合を特別扱いし、最適予測次数Poptが最小予測次数Pminである場合に短い符号長の符号（例えば、符号長が最短の符号）を割り当ててもよい。
また、複数種類設定された予備集合Spreから部分集合Ssubを選択する場合、部分集合Ssubの要素から選択される最適予測次数Poptを特定する情報を可変長符号化するための符号化規則（例えば、符号表）を予備集合Spreごとに独立に設定しておき、それらの符号化規則を、部分集合Ssubとして選択された予備集合Spreごとに使い分けることが望ましい。これにより、すべての予備集合Spreに共通の符号化規則を用いる場合に比べて符号長を短くできる。
〔第１実施形態〕
次に、本発明の第１実施形態を説明する。本形態は、最適予測次数を指定するための符号が、最適予測次数を特定する情報を可変長符号化した符号である形態である。
図５は、第１実施形態の符号化装置１００の機能構成を説明するためのブロック図であり、図６は、図５の線形予測分析部１１２の機能構成を説明するためのブロック図である。また、図７は、第１実施形態の復号装置２００の機能構成を説明するためのブロック図である。なお、これらの図において図１や図２と同じ構成については、図１や図２と同じ符号を用いる。
図５に示すように、本形態の符号化装置１００は、フレームバッファ１１１１、線形予測分析部１１２、量子化部１１１３、係数符号化部１１１４、線形予測係数変換部１１１５、線形予測部１１１６、減算部１１１７、最適予測次数符号化部１１８、残差符号化部１１２０、及び合成部１１３０を有する。また、図６に示すように、線形予測分析部１１２は、係数算出部１１２ａ、予測次数指定符号量算出部１１２ｂ、係数符号量算出部１１２ｃ、予測残差符号量算出部１１２ｄ、最適予測次数探索部１１２ｅ、及び選択部１１２ｆを有する。
図７に示すように、本形態の復号装置２００は、最適予測次数復元部２３５、分離部１２１０、係数復号部１２３１、線形予測係数変換部１２３２、残差復号部１２２０、線形予測部１２３３、及び加算部１２３４を有する。
なお、本形態の符号化装置１００や復号装置２００は、例えば、ＣＰＵ（central processing unit）、ＲＡＭ（random-access memory）、ＲＯＭ（read-only memory）等を備えた公知のコンピュータ又は専用のコンピュータに所定のプログラムが読み込まれ、ＣＰＵがそれを実行することによって構成された特別な装置である。すなわち、符号化装置１００のフレームバッファ１１１１は、例えば、ＲＡＭ、キャッシュメモリ、レジスタ等のメモリであり、線形予測分析部１１２、量子化部１１１３、係数符号化部１１１４、線形予測係数変換部１１１５、線形予測部１１１６、減算部１１１７、最適予測次数符号化部１１８、残差符号化部１１２０、及び合成部１１３０は、例えば、ＣＰＵが所定のプログラムを実行することで構築される処理部である。また、復号装置２００の最適予測次数復元部２３５、分離部１２１０、係数復号部１２３１、線形予測係数変換部１２３２、残差復号部１２２０、線形予測部１２３３、及び加算部１２３４は、例えば、ＣＰＵが所定のプログラムを実行することで構築される処理部である。また、これらの処理部の少なくとも一部が集積回路等の電子回路によって構成されてもよい。さらに、必要に応じ、符号化装置１００や復号装置２００に、各処理部の処理によって出力されたデータを格納し、各処理部の別の処理時にデータが読み出される一時メモリを設けてもよい。また、このような各処理部の実現方法は、以下の各実施形態やその変形例でも同様である。
＜前処理＞
前処理として、線形予測分析部１１２の係数算出部１１２ａは、上述したように間引きされた部分集合Ssubの各要素を利用できるように設定されている。
また、符号化装置１００の最適予測次数符号化部１１８は、最適予測次数Poptを特定する情報を前述のように可変長符号化するための符号表が設定されている。また、復号装置２００の最適予測次数復元部２３５は、最適予測次数符号化部１１８で可変長符号化された符号を復号できるように当該符号表が設定されている。図１０（Ａ）は、この可変長符号化のための符号表の例示である。この例の符号表では、部分集合Ssubの各要素である間引かれた予測次数Pとそれに対応する符号と符号長とが対応付けられている。
＜符号化方法＞
図８は、第１実施形態の符号化方法を説明するためのフローチャートである。以下、図８を用いて本形態の符号化方法を説明する。なお、以下では１フレーム分の処理のみを説明するが、実際は同様な処理が各フレームについて実行される。
符号化装置１００（図５）のフレームバッファ１１１１には、標本化・量子化されたPCM形式の時系列信号x(n)が入力される。これらの時系列信号x(n)は、線形量子化（「一様量子化」と呼ぶ場合もある）されたものであってもよいし、圧伸量子化（例えば、ITU-T Recommendation G.711, “Pulse Code Modulation (PCM) of Voice Frequencies”参照）のような非線形量子化（「非一様量子化」と呼ぶ場合もある）されたものであってもよい。また、時系列信号x(n)がPCM形式の信号ではなく、量子化されていない信号であってもよい。フレームバッファ１１１１は、１フレーム分の時系列信号x(n)（n=1,...,N）をバッファし、符号化装置１００は、フレーム毎に時系列信号x(n)（n=1,...,N）を符号化する。
１フレーム分の時系列信号x(n)（n=1,...,N）は、線形予測分析部１１２（図６）に送られる。線形予測分析部１１２の係数算出部１１２ａは、Pmin以上かつPmax以下のすべての整数を要素とする整数集合Sから一部の要素を間引いた部分集合Ssub⊂Sの各要素を各予測次数Pとする。係数算出部１１２ａは、入力された時系列信号x(n)（n=1,...,N）の線形予測分析を、当該各予測次数Pについてそれぞれ行う。これにより、係数算出部１１２ａは、各予測次数Pにそれぞれ対応する、1次からP次までの各PARCOR係数k(m)(m=1,...,P)を算出して出力する（ステップＳ１１）。なお、係数算出部１１２ａは、時系列信号x(n)（n=1,...,N）をそのまま線形予測分析する構成であってもよいし、非線形量子化されて入力された時系列信号x(n)（n=1,...,N）を線形量子化や他の非線形量子化にマッピングしてから線形予測分析を行う構成であってもよい。すなわち、「時系列信号（n=1,...,N）の線形予測分析」とは、時系列信号x(n)（n=1,...,N）をそのまま線形予測分析すること、又は、非線形量子化されて入力された時系列信号x(n)（n=1,...,N）を線形量子化や他の非線形量子化にマッピングしてから線形予測分析することである。また、PARCOR係数k(m)(m=1,...,P)の算出は、レビンソン・ダービン（Levinson-Durbin）法やバーグ（Burg）法などの逐次的方法によって行われてもよいし、自己相関法や共分散法のように予測次数ごとに連立方程式（予測残差を最小にする線形予測係数を解とする連立方程式）を解くことによって行われてもよい（「守谷健弘著、“音声符号化”、社団法人電子情報通信学会、１９９８」「日野幹雄著、“スペクトル解析”、朝倉書店、１９７９」など参照）。
また、係数算出部１１２ａで用いられた各予測次数Pは予測次数指定符号量算出部１１２ｂに送られ、予測次数指定符号量算出部１１２ｂは各予測次数Pを指定するための符号の符号量apを算出する（ステップＳ１２）。本形態の場合、予測次数指定符号量算出部１１２ｂは、例えば、最適予測次数符号化部１１８で用いるのと同じ符号表（例えば、図１０（Ａ））を用いて各予測次数Pを可変長符号化した符号の符号長を抽出し、それを符号量apとする。また、係数算出部１１２ａで各予測次数Pについてそれぞれ算出された1次からP次までの各PARCOR係数k(m)(m=1,...,P)は、係数符号量算出部１１２ｃに送られる。係数符号量算出部１１２ｃは、予測次数Pごとに、PARCOR係数k(m)(m=1,...,P)の符号量akpを算出する（ステップＳ１２）。言い換えると、係数符号量算出部１１２ｃは、予測次数Pごとに、当該予測次数Pに対応する1次からP次までのPARCOR係数k(m)(m=1,...,P)の符号量akpを算出する。具体的には、係数符号量算出部１１２ｃは、例えば、係数符号化部１１１４と同じ符号化方法によってPARCOR係数k(m)(m=1,...,P)を符号化し、予測次数Pごとに、当該予測次数Pに対応する1次からP次までのPARCOR係数k(m)(m=1,...,P)の符号量akpを算出する。また、PARCOR係数k(m)(m=1,...,P)の振幅値やエネルギーを用いて、予測次数Pごとに、当該予測次数Pに対応する1次からP次までのPARCOR係数k(m)(m=1,...,P)の推定符号量akpを求めてもよい。また、符号量や推定符号量akpは、予測次数Pごとに各PARCOR係数k(1),...,k(P)をまとめて符号化した場合のものであってもよいし、各PARCOR係数k(1),...,k(P)ごとの符号量又は推定符号量を予測次数Pごとに合算したものでもよい。言い換えると、各予測次数Pにそれぞれ対応する符号量や推定符号量akpは、当該予測次数Pに対応する１次からP次までの各PARCOR係数k(1),...,k(P)をまとめて１つの符号化対象として符号化して得られる符号の符号量やその推定量であってもよいし、当該予測次数Pに対応する１次からP次までのPARCOR係数k(1),...,k(P)をそれぞれ別個の符号化対象として符号化して得られる各符号の符号量又は推定符号量を合算したものであってもよい。
また、係数算出部１１２ａで算出された各予測次数Pでの各PARCOR係数k(m)（m=1,...,P）は予測残差符号量算出部１１２ｄにも送られ、予測残差符号量算出部１１２ｄは、予測次数Pごとに１フレーム分の予測残差e(n)（n=1,...,N）の符号量akpを算出する（ステップＳ１３）。具体的には、予測残差符号量算出部１１２ｄは、例えば、１フレーム分の時系列信号x(n)（n=1,...,N）と各PARCOR係数k(m)（m=1,...,P）とを用い、実際に、予測次数Pに対応する１フレーム分の予測残差e(n)（n=1,...,N）を求め、残差符号化部１１２０と同じ符号化方法によってそれらの符号化を行って、当該予測次数Pに対応する１フレーム分の予測残差e(n)（n=1,...,N）の符号量akpを算出する。なお、１フレーム分の予測残差e(n)（n=1,...,N）の符号化は、各予測残差e(n)をそれぞれ別個に符号化対象とするものであってもよいし、１フレーム分の予測残差e(n)（n=1,...,N）をまとめて１つの符号化対象とするものであってもよい。
また、予測残差符号量算出部１１２ｄが、予測残差e(n)の振幅値やエネルギーを用いて推定符号量akpを算出してもよい。例えば、時系列入力信号x(n)の１フレームあたりのエネルギーE(0)は
Figure 0004848049
となる。１次のPARCOR係数k(1)を用いて、１次の線形予測を行ったときの予測残差のエネルギーE(1)は、
E(1)=E(0){1-k(1)2} …(5)
となる。２次の線形予測を行ったときの予測残差のエネルギーE(2)は、
E(2)=E(1){1-k(2)2} …(6)
となる。これを予測次数Pまで繰り返すと、P次の線形予測を行ったときの予測残差のエネルギーE(P)は、
Figure 0004848049
となる。予測残差符号量算出部１１２ｄは、この残差エネルギーを用いて予測残差e(n)（n=1,...,N）の推定符号量akpを算出することができる。例えば、残差符号化部１２０がアルファ符号で符号化を行うのであれば、正の予測残差e(n)の符号量は2・e(n)+1であるため、１フレーム分の予測残差e(n)（n=1,...,N）の推定符号量akpは
Figure 0004848049
で算出できる。また、残差符号化部１２０がゴロム・ライス符号（Golomb-Rice Coding）で符号化を行う場合であっても、ゴロム・ライス符号の前段ビットはアルファ符号であり、後段ビットはライスパラメータに応じた有限長の符号であるため、式(8)で１フレーム分の予測残差e(n)の推定符号量akpを算出することは妥当である。
ステップＳ１２からＳ１４で算出された予測次数Pを指定するための符号の符号量apと、PARCOR係数k(m)（m=1,...,P）の符号量又は推定符号量akpと、予測残差e(n)（n=1,...,N）の符号量又は推定符号量aepとは、最適予測次数探索部１１２ｅに入力される。最適予測次数探索部１１２ｅは、これらの合計値SUMを指標として、各予測次数Pから最適予測次数Poptを選択する（ステップＳ１５）。この選択方法の具体例は前述した通り、例えば、合計値SUMを最小にする予測次数Pを最適予測次数Poptとして選択する。
ステップＳ１５で選択された最適予測次数Poptは最適予測次数符号化部１１８に入力され、最適予測次数符号化部１１８は、それを可変長符号化した符号CPopt（「最適予測次数Poptを指定するための符号」に相当）を生成して出力する（ステップＳ１６）。前述した通り、この可変長符号化方法は、第１符号化対象が第２符号化対象よりも予め定められた特定の値に近い場合に、当該第１符号化対象に当該第２符号化対象の符号よりも短い符号長の符号が割り当てられる頻度が、当該第１符号化対象に当該第２符号化対象の符号よりも長い符号長の符号が割り当てられる頻度よりも高い可変長符号化方法である。また、予め定められた特定の値の一例は最大予測次数Pmaxである。このような可変長符号化方法の具体例としては、ハフマン符号化（Huffman Coding）などを例示できる。
また、ステップＳ１５で選択された最適予測次数Poptは選択部１１２ｆに入力される。選択部１１２ｆは、係数算出部１１２ａで算出された各予測次数Pに対応するPARCOR係数k(m)(m=1,...,P)から、最適予測次数Poptに対応するPARCOR係数k(m)(m=1,...,Popt)を選択し、出力する（ステップＳ１７）。
ステップＳ１７で選択された最適予測次数Poptに対応するPARCOR係数k(m)(m=1,...,Popt)は、量子化部１１１３に送られ、量子化されて１次からPopt次までの量子化PARCOR係数i(m)(m=1,...,Popt)が生成される（ステップＳ１８）。
次に、係数符号化部１１１４が、ステップＳ１８で生成された１次からPopt次までの量子化PARCOR係数i(m)(m=1,...,Popt)に対応する係数符号Ckを生成して出力する（ステップＳ１９）。この符号化は、ゴロム・ライス符号化方法やハフマン符号化方法などのエントロピー符号化方法によって行われる。
また、ステップＳ１５で選択された最適予測次数Poptや、ステップＳ１７で選択された最適予測次数PoptでのPARCOR係数k(m)(m=1,...,Popt)は線形予測係数変換部１１１５にも送られる。線形予測係数変換部１１１５は、これらを用いて予測次数Poptの線形予測フィルタの各線形予測係数α(m)（m=1,...,Popt）を算出する（ステップＳ２０）。この算出は、例えば、量子化PARCOR係数i(m)を逆量子化した後、レビンソン（Levinson）・アルゴリズムを適用して行う。なお、量子化PARCOR係数i(m)を逆量子化する処理は、量子化PARCOR係数i(m)に対応するPARCOR係数k(m)の範囲のうちの予め定めた何れかの値k'(m)を求める処理である。例えば、η1≦k(m)＜η2のPARCOR係数k(m)に対応する量子化ＰＡＲＣＯＲ係数がi(m)である場合、量子化PARCOR係数i(m)を逆量子化したPARCOR係数k'(m)の例はη1とη2の平均値である。
線形予測部１１１６は、入力された１フレーム分の時系列信号x(n)（n=1,...,N）と各線形予測係数α(m)(m=1,...,Popt)と最適予測次数Poptとを用い、各線形予測係数α(m)(m=1,...,Popt)をフィルタ係数とする、最適予測次数Poptでの線形予測フィルタy(n)（n=1,...,N）を生成して出力する。減算部１１１７は、入力された線形予測フィルタy(n)（n=1,...,N）と、時系列信号x(n)とを用い、予測残差e(n)を算出する（ステップＳ２１）。なお、線形予測部１１１６及び減算部１１１７は、例えば、時系列信号x(n)（n=1,...,N）をそのまま用いて予測フィルタ処理を行う。また、時系列信号x(n)（n=1,...,N）が非線形量子化された信号である場合には、線形予測部１１１６及び減算部１１１７は、時系列信号x(n)（n=1,...,N）を線形量子化領域や他の非線形量子領域にマッピングしてから予測フィルタ処理を行ってもよい。また、この場合に、線形予測部１１１６及び減算部１１１７が、各線形予測係数α(m)と時系列信号x(n)とを線形予測フィルタ（式(2)）に適用して算出された線形予測値y(n)を非線形量子領域にマッピングしてからから予測フィルタ処理を行ってもよい。
算出された予測残差e(n)（n=1,...,N）は残差符号化部１１２０に送られ、そこでエントロピー符号化され、予測残差e(n)（n=1,...,N）に対応する残差符号Ceが生成される（ステップＳ２２）。最適予測次数符号化部１１８で生成された符号CPoptと、係数符号化部１１１４で生成された係数符号Ckと、残差符号化部１１２０で生成された残差符号Ceは、合成部１１３０に送られ、そこで合成されて符号Cgが生成される（ステップＳ２３）。そして、符号化装置１００は、生成した符号Cgを出力する。
＜復号方法＞
図９は、第１実施形態の復号方法を説明するためのフローチャートである。以下、図９を用いて本形態の復号方法を説明する。なお、以下では１フレーム分の処理のみを説明するが、実際は同様な処理が各フレームについて実行される。
復号装置２００の分離部１２１０は、入力された符号Cgを分離し、最適予測次数Poptを指定するための符号CPoptと、１次からPopt次までの量子化PARCOR係数i(m)(m=1,...,Popt)に対応する係数符号Ckと、予測残差e(n)(n=1,...,N)に対応する残差符号Ceとを生成して出力する（ステップＳ３１）。最適予測次数復元部２３５は、ステップＳ１６の可変長符号化方式に対応する復号方式によって符号CPoptを復号し、最適予測次数Poptを生成して出力する（ステップＳ３２）。また、残差復号部１２２０は、分離部１０２１から出力された残差符号Ce（予測次数が最適予測次数である線形予測分析を時系列信号に対して行った結果から得られた予測残差の符号）を復号し、予測残差e(n)(n=1,...,N)を生成して出力する（ステップＳ３３）。また、係数復号部１２３１は、分離部１２１０から出力された係数符号Ck（予測次数が最適予測次数である線形予測分析を時系列信号に対して行って得られた線形予測係数若しくは当該線形予測係数に変換可能な係数の符号）を復号し、量子化PARCOR係数i(m)(m=1,...,Popt)を生成して出力する（ステップＳ３４）。最適予測次数復元部２３５で生成された最適予測次数Poptと、係数復号部１２３１で生成された量子化PARCOR係数i(m)(m=1,...,Popt)とは、線形予測係数変換部１２３２に入力される。線形予測係数変換部１２３２は、量子化PARCOR係数i(m)(m=1,...,Popt)を逆量子化したPARCOR係数k'(m)(m=1,...,Popt)を用い、最適予測次数Poptでの各線形予測係数α(m)（m=1,...,Popt）を算出して出力する（ステップＳ３５）。この処理はレビンソン・アルゴリズムを適用して行う。
線形予測部１２３３と加算部１２３４は、最適予測次数復元部２３５で生成された最適予測次数Poptと、線形予測係数変換部１２３２で生成された各線形予測係数α(m)（m=1,...,Popt）と、残差復号部１２２０から出力された予測残差e(n)（n=1,...,N）とを用い、逆予測フィルタ処理によって時系列信号のロスレス復元値x(n)（n=1,...,N）を生成して出力する（ステップＳ３６）。言い換えると、線形予測部１２３３と加算部１２３４は、各線形予測係数α(m)（m=1,...,Popt）をフィルタ係数とする線形予測フィルタと、予測残差e(n)（n=1,...,N）とを用い、逆予測フィルタ処理によって時系列信号のロスレス復元値x(n)（n=1,...,N）を生成して出力する。より具体的には、例えば、線形予測部１２３３が、１次からPopt次までの各線形予測係数α(m)（m=1,...,Popt）をフィルタ係数とする、最適予測次数Poptの線形予測フィルタによって線形予測値y(n)（n=1,...,N）を生成する。そして、加算部１２３４が、線形予測値y(n)と予測残差e(n)とをそれぞれ加算して時系列信号のロスレス復元値x(n)（n=1,...,N）を生成する。
〔第１実施形態の変形例〕
第１実施形態の変形例では、最小予測次数Pmin以上かつ最大予測次数Pmax以下のすべての整数を要素とする整数集合Sから一部の要素を間引いた予備集合Spre⊂Sが複数種類設定されており、複数種類設定された予備集合Spreの１つを部分集合Ssubとして選択する。
この場合、符号化装置１００の線形予測分析部１１２は、さらに、複数種類設定された予備集合Spreから部分集合Ssubを選択する探索次数選択部１１２ｇを有する（図６）。また、前処理において、線形予測分析部１１２の係数算出部１１２ａには、複数種類の予備集合Spreが設定され、探索次数選択部１１２ｇの指示に従って予備集合Spreの１つを部分集合Ssubとして選択可能にされる。
この場合の符号化方法では、ステップＳ１１の前に、探索次数選択部１１２ｇが複数種類設定された予備集合Spreの１つを部分集合Ssubとして選択する。この選択はフレームやサブフレーム単位で行なわれてもよいし、複数のフレーム単位で行なわれてもよい。また、探索次数選択部１１２ｇが行う部分集合Ssubの選択は、外部から探索次数選択部１１２ｇに入力された指示信号（例えば、フレーム長を特定する信号）に応じて行なわれてもよいし、前述のように入力された時系列信号x(n)に応じて行なわれてもよい。選択された部分集合Ssubを特定する識別子subは、係数算出部１１２ａに送られ、係数算出部１１２ａはそれに従って部分集合Ssubを選択し、ステップＳ１１以降の処理が実行される。
ここで、最適予測次数符号化部１１８が、すべての予備集合Spreに共通の符号化規則で部分集合Ssubの要素から選択された最適予測次数Poptを特定する情報を可変長符号化するのであれば、符号化方法の他の処理、及び復号方法の処理は第１実施形態と同じである。一方、部分集合Ssubを可変長符号化するための符号化規則を予備集合Spreごとに独立に設定しておき、最適予測次数符号化部１１８が、それらの符号化規則を、部分集合Ssubとして選択された予備集合Spreごとに使い分ける場合には以下の相違点がある。
まず、ステップＳ１６の処理において、探索次数選択部１１２ｇが選択した部分集合Ssubを特定する識別子subを最適予測次数符号化部１１８に送り、最適予測次数符号化部１１８は、識別子subに応じて最適予測次数Poptを可変長符号化するための符号化規則（例えば、符号表）を選択し、選択した符号化規則に従って最適予測次数Poptを特定する情報を可変長符号化して出力する。例えば、最適予測次数符号化部１１８が図１０（Ａ）（Ｂ）の２つの符号表を利用できるように設定しておく。そして、例えば、識別子subが部分集合Ssub=(15,12,7,2,0)を示すものであったならば、最適予測次数符号化部１１８は、図１０（Ａ）の符号表に従って最適予測次数Poptを特定する情報を可変長符号化する。一方、識別子subが部分集合Ssub=(15,13,10,7,4)を示すものであったならば、最適予測次数符号化部１１８は、図１０（Ｂ）の符号表に従って最適予測次数Poptを特定する情報を可変長符号化する。
また、このような場合には、復号装置２００でも、最適予測次数Poptを特定する情報がどのような符号化規則に従って符号化されたかを知る必要がある。そのため、ステップＳ２３の処理において、合成部１１３０は、最適予測次数符号化部１１８で生成された符号CPoptと、係数符号化部１１１４で生成された係数符号Ckと、残差符号化部１１２０で生成された残差符号Ceと、識別子subとを合成した符号Cgを生成して出力する。そして、復号装置２００の分離部１２１０は、ステップＳ３１において、入力された符号Cgを分離し、最適予測次数Poptを指定するための符号CPoptと、１次からPopt次までの量子化PARCOR係数i(m)（m=1,...,Popt）に対応する係数符号Ckと、予測残差e(n)（n=1,...,N）に対応する残差符号Ceと、識別子subとを生成して出力する。最適予測次数復元部２３５は、ステップＳ３２において、識別子subを用いてステップＳ１６の可変長符号化方式に対応する復号方式を特定し、それによって符号CPoptを復号し、最適予測次数Poptを生成して出力する。
〔第２実施形態〕
本形態は、第１実施形態の変形例である。第１実施形態では、符号化の際、予測次数Pを指定するための符号の符号量apと、PARCOR係数k(m)（m=1,...,P）の符号量又は推定符号量akpと、予測残差e(n)（n=1,...,N）の符号量又は推定符号量aepとの合計値SUMを指標として、各予測次数Pから最適予測次数Poptを選択していた（ステップＳ１５）。しかし、第２実施形態では、符号化の際、PARCOR係数k(m)（m=1,...,P）の符号量又は推定符号量akpと、予測残差e(n)（n=1,...,N）の符号量又は推定符号量aepとの合計値SUMを指標として、各予測次数Pから最適予測次数Poptを選択する。
図１１は、第２実施形態の線形予測分析部３１２の機能構成を説明するためのブロック図である。また、図１２は、第２実施形態の符号化方法の一部を説明するためのフローチャートである。
図１１に示すように、本形態の線形予測分析部３１２は、予測次数指定符号量算出部１１２ｂが存在せず、最適予測次数探索部１１２ｅが最適予測次数探索部３１２ｅに置換された点で、第１実施形態の線形予測分析部１１２と相違する。そして、図１２に示すように、本形態の符号化方法は、ステップＳ１２の処理が実行されず、ステップＳ１５の代わりに、最適予測次数探索部３１２ｅが、PARCOR係数k(m)（m=1,...,P）の符号量又は推定符号量akpと、予測残差e(n)（n=1,...,N）の符号量又は推定符号量akpとの合計値SUM指標とし、各予測次数Pから最適予測次数Popを選択する。
なお、復号装置の構成や復号方法は第１実施形態と同じである。また、第１実施形態の変形例を第２実施形態に適用することも可能である。
〔その他の変形例〕
なお、本発明は上述の実施の形態に限定されるものではない。例えば、第１実施形態では、線形予測分析部の係数算出部が、各予測次数Pでの1次からP次までのPARCOR係数k(m)（m=1,...,P）を生成し、各予測次数Pを指定するための符号の符号量と、各予測次数PでのPARCOR係数k(m)（m=1,...,P）の符号量又は推定符号量と、残差e(n)（n=1,...,N）の符号量又は推定符号量との合計値SUMを指標とし、各予測次数Pから最適予測次数Poptを選択していた。また、例えば、第２実施形態では、線形予測分析部の係数算出部が、各予測次数Pでの1次からP次までのPARCOR係数k(m)（m=1,...,P）を生成し、各予測次数PでのPARCOR係数k(m)（m=1,...,P）の符号量又は推定符号量と、残差e(n)（n=1,...,N）の符号量又は推定符号量との合計値SUMを指標とし、各予測次数Pから最適予測次数Poptを選択していた。しかし、線形予測分析部の係数算出部が、各予測次数Pでの線形予測係数α(m)（m=1,...,P）を生成し、最適予測次数探索部が、各予測次数PでのPARCOR係数k(m)（m=1,...,P）の符号量又は推定符号量の代わりに、各予測次数Pでの線形予測係数α(m)（m=1,...,P）の符号量又は推定符号量を用いて各予測次数Pから最適予測次数Poptを選択してもよい。
また、線形予測分析部の係数算出部で、部分集合Ssubの各要素である各予測次数Pでの線形予測係数若しくは当該線形予測係数に変換可能な係数をそれぞれ算出し、各予測次数Pに対し、それぞれ予測次数Pの符号と、係数の符号と、予測残差e(n)（n=1,...,N）の符号とをすべて生成し、最適予測次数探索部が実際の合計符号量を比較して各予測次数Pから最適予測次数Poptを選択してもよい。
また、上記の実施形態では、最適予測次数符号化部１１８が、最適予測次数Poptを特定する情報を可変長符号化することとした。しかし、最適予測次数Poptを特定する情報を可変長符号化することなく、最適予測次数Poptを特定する情報がそのまま符号Cgに含まれもよい。
また、上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。
また、上述の構成をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。
また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。
このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。
また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
本発明の産業上の利用分野としては、例えば、音響信号の可逆圧縮符号化・復号化技術を例示できる。また、本発明は、音響信号以外にも、映像信号、生体信号、地震波信号、センサーアレイ信号などの可逆圧縮符号化・復号化技術にも応用可能である。
１００，１１００ 符号化装置
２００，１２００ 復号装置

Claims (13)

  1. (A) 入力された時系列信号から、線形予測係数若しくは当該線形予測係数に変換可能な係数を算出するステップと、
    (B) 最小予測次数を予め定められた０以上の整数とし、最大予測次数を当該最小予測次数よりも大きな予め定められた整数とした場合における、当該最小予測次数以上かつ当該最大予測次数以下のすべての整数を要素とする整数集合から当該最大予測次数に近い領域ほど要素数が多くなるように不等間隔で一部の要素を間引いた部分集合の要素である各予測次数について、当該予測次数までの前記線形予測係数若しくは当該線形予測係数に変換可能な係数の符号量又は推定符号量と、当該予測次数までの前記線形予測係数を用いて前記時系列信号の線形予測を行った場合の予測残差の符号量又は推定符号量と、の合計値に相当する値を算出するステップと、
    (C) 前記合計値に相当する値が最小となる予測次数を最適予測次数として求め、当該最適予測次数、を指定するための符号を出力するステップと、
    を有する符号化方法。
  2. (A) 入力された時系列信号から、線形予測係数若しくは当該線形予測係数に変換可能な係数を算出するステップと、
    (B) 最小予測次数を予め定められた０以上の整数とし、最大予測次数を当該最小予測次数よりも大きな予め定められた整数とした場合における、当該最小予測次数以上かつ当該最大予測次数以下のすべての整数を要素とする整数集合から当該最大予測次数に近い領域ほど要素数が多くなるように不等間隔で一部の要素を間引いた部分集合の要素である各予測次数について、当該予測次数を指定するための符号の符号量と、当該予測次数までの前記線形予測係数若しくは当該線形予測係数に変換可能な係数の符号量又は推定符号量と、当該予測次数までの前記線形予測係数を用いて前記時系列信号の線形予測を行った場合の予測残差の符号量又は推定符号量と、の合計値に相当する値を算出するステップと、
    (C) 前記合計値に相当する値が最小となる予測次数を最適予測次数として求め、当該最適予測次数を指定するための符号を出力するステップと、
    を有する符号化方法。
  3. 請求項１又は２の符号化方法であって、
    前記最小予測次数以上かつ前記最大予測次数以下のすべての整数を要素とする整数集合から当該最大予測次数に近い領域ほど要素数が多くなるように不等間隔で一部の要素を間引いた予備集合が複数種類設定されており、
    前記複数種類設定された予備集合の１つを前記部分集合として選択するステップをさらに有する、
    ことを特徴とする符号化方法。
  4. 請求項１又は２の符号化方法であって、
    第１符号化対象が第２符号化対象よりも予め定められた特定の値に近い場合に、当該第１符号化対象に当該第２符号化対象の符号よりも短い符号長の符号が割り当てられる頻度が、当該第１符号化対象に当該第２符号化対象の符号よりも長い符号長の符号が割り当てられる頻度よりも高い可変長符号化方法を用い、前記最適予測次数を特定する情報を可変長符号化するステップをさらに有し、
    前記最適予測次数を指定するための符号は、前記最適予測次数を特定する情報を可変長符号化した符号である、
    ことを特徴とする符号化方法。
  5. 請求項４の符号化方法であって、
    前記予め定められた特定の値は、前記最大予測次数を特定する情報である、
    ことを特徴とする符号化方法。
  6. (A) 最小予測次数を予め定められた０以上の整数とし、最大予測次数を当該最小予測次数よりも大きな予め定められた整数とした場合における、当該最小予測次数以上かつ当該最大予測次数以下のすべての整数を要素とする整数集合から当該最大予測次数に近い領域ほど要素数が多くなるように不等間隔で一部の要素を間引いた部分集合の要素である最適予測次数を指定するための符号が入力され、当該最適予測次数を得るステップと、
    (B) 残差符号を復号して予測残差を得るステップと、
    (C)係数符号を復号して、前記最適予測次数までの線形予測係数若しくは当該線形予測係数に変換可能な係数を得るステップと、
    (D) 前記線形予測係数をフィルタ係数とする線形予測フィルタと前記予測残差とを用いて時系列信号を合成して出力するステップと、
    を有する復号方法。
  7. 入力された時系列信号から、線形予測係数若しくは当該線形予測係数に変換可能な係数を算出する係数算出部と、
    最小予測次数を予め定められた０以上の整数とし、最大予測次数を当該最小予測次数よりも大きな予め定められた整数とした場合における、当該最小予測次数以上かつ当該最大予測次数以下のすべての整数を要素とする整数集合から当該最大予測次数に近い領域ほど要素数が多くなるように不等間隔で一部の要素を間引いた部分集合の要素である各予測次数について、当該予測次数までの前記線形予測係数若しくは当該線形予測係数に変換可能な係数の符号量又は推定符号量と、当該予測次数までの前記線形予測係数を用いて前記時系列信号の線形予測を行った場合の予測残差の符号量又は推定符号量と、の合計値に相当する値を算出し、前記合計値に相当する値が最小となる予測次数を最適予測次数として求め、当該最適予測次数を指定するための符号を出力する最適予測次数探索部と、
    を有する符号化装置。
  8. 入力された時系列信号から、線形予測係数若しくは当該線形予測係数に変換可能な係数を算出する係数算出部と、
    最小予測次数を予め定められた０以上の整数とし、最大予測次数を当該最小予測次数よりも大きな予め定められた整数とした場合における、当該最小予測次数以上かつ当該最大予測次数以下のすべての整数を要素とする整数集合から当該最大予測次数に近い領域ほど要素数が多くなるように不等間隔で一部の要素を間引いた部分集合の要素である各予測次数について、当該予測次数を指定するための符号の符号量と、当該予測次数までの前記線形予測係数若しくは当該線形予測係数に変換可能な係数の符号量又は推定符号量と、当該予測次数までの前記線形予測係数を用いて前記時系列信号の線形予測を行った場合の予測残差の符号量又は推定符号量と、の合計値に相当する値を算出し、前記合計値に相当する値が最小となる予測次数を最適予測次数として求め、当該最適予測次数を指定するための符号を出力する最適予測次数探索部と、
    を有する符号化装置。
  9. 最小予測次数を予め定められた０以上の整数とし、最大予測次数を当該最小予測次数よりも大きな予め定められた整数とした場合における、当該最小予測次数以上かつ当該最大予測次数以下のすべての整数を要素とする整数集合から当該最大予測次数に近い領域ほど要素数が多くなるように不等間隔で一部の要素を間引いた部分集合の要素である最適予測次数を指定するための符号が入力され、当該最適予測次数を得る最適予測次数復元部と、
    残差符号を復号して予測残差を得る予測残差復号部と、
    係数符号を復号して、前記最適予測次数までの線形予測係数若しくは当該線形予測係数に変換可能な係数を得る係数復号部と、
    前記線形予測係数をフィルタ係数とする線形予測フィルタと前記予測残差とを用いて時系列信号を合成して出力する復元処理部と、
    を有する復号装置。
  10. 請求項１又は２の符号化方法の処理をコンピュータに実行させるためのプログラム。
  11. 請求項６の復号方法の処理をコンピュータに実行させるためのプログラム。
  12. 請求項１又は２の符号化方法の処理をコンピュータに実行させるためのプログラムを格納したコンピュータ読み取り可能な記録媒体。
  13. 請求項６の復号方法の処理をコンピュータに実行させるためのプログラムを格納したコンピュータ読み取り可能な記録媒体。
JP2010542108A 2008-12-09 2009-12-08 符号化方法、復号方法、それらの装置、プログラム及び記録媒体 Active JP4848049B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010542108A JP4848049B2 (ja) 2008-12-09 2009-12-08 符号化方法、復号方法、それらの装置、プログラム及び記録媒体

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2008313530 2008-12-09
JP2008313530 2008-12-09
PCT/JP2009/070551 WO2010067800A1 (ja) 2008-12-09 2009-12-08 符号化方法、復号方法、それらの装置、プログラム及び記録媒体
JP2010542108A JP4848049B2 (ja) 2008-12-09 2009-12-08 符号化方法、復号方法、それらの装置、プログラム及び記録媒体

Publications (2)

Publication Number Publication Date
JP4848049B2 true JP4848049B2 (ja) 2011-12-28
JPWO2010067800A1 JPWO2010067800A1 (ja) 2012-05-17

Family

ID=42242789

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010542108A Active JP4848049B2 (ja) 2008-12-09 2009-12-08 符号化方法、復号方法、それらの装置、プログラム及び記録媒体

Country Status (2)

Country Link
JP (1) JP4848049B2 (ja)
WO (1) WO2010067800A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110444215A (zh) * 2014-05-01 2019-11-12 日本电信电话株式会社 编码装置、编码方法、程序、记录介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019198383A1 (ja) * 2018-04-13 2019-10-17 日本電信電話株式会社 符号化装置、復号装置、符号化方法、復号方法、プログラム、および記録媒体

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0511799A (ja) * 1991-07-08 1993-01-22 Fujitsu Ltd 音声符号化方式
JP2007279393A (ja) * 2006-04-06 2007-10-25 Nippon Telegr & Teleph Corp <Ntt> Ｐａｒｃｏｒ係数算出装置、ｐａｒｃｏｒ係数算出方法、そのプログラムおよびその記録媒体
JP2007286200A (ja) * 2006-04-13 2007-11-01 Nippon Telegr & Teleph Corp <Ntt> 適応ブロック長符号化装置、その方法、プログラム及び記録媒体
JP2007286146A (ja) * 2006-04-13 2007-11-01 Nippon Telegr & Teleph Corp <Ntt> 適応ブロック長符号化装置、その方法、プログラム及び記録媒体
JP2007318691A (ja) * 2006-05-29 2007-12-06 Nippon Telegr & Teleph Corp <Ntt> 線形予測モデル次数決定装置、線形予測モデル次数決定方法、そのプログラムおよび記録媒体
JP2008026460A (ja) * 2006-07-19 2008-02-07 Nippon Telegr & Teleph Corp <Ntt> 符号化装置、それらの方法、プログラム及び記録媒体

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4838774B2 (ja) * 2007-07-18 2011-12-14 日本電信電話株式会社 マルチチャネル線形予測符号化の予測係数決定方法及び装置、プログラム、記録媒体
JP4838773B2 (ja) * 2007-07-18 2011-12-14 日本電信電話株式会社 線形予測符号化の予測次数決定方法及びそれを使った予測係数決定方法及び装置、プログラム、その記録媒体
JP4691082B2 (ja) * 2007-09-11 2011-06-01 日本電信電話株式会社 線形予測モデル次数決定装置、線形予測モデル次数決定方法、そのプログラムおよび記録媒体

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0511799A (ja) * 1991-07-08 1993-01-22 Fujitsu Ltd 音声符号化方式
JP2007279393A (ja) * 2006-04-06 2007-10-25 Nippon Telegr & Teleph Corp <Ntt> Ｐａｒｃｏｒ係数算出装置、ｐａｒｃｏｒ係数算出方法、そのプログラムおよびその記録媒体
JP2007286200A (ja) * 2006-04-13 2007-11-01 Nippon Telegr & Teleph Corp <Ntt> 適応ブロック長符号化装置、その方法、プログラム及び記録媒体
JP2007286146A (ja) * 2006-04-13 2007-11-01 Nippon Telegr & Teleph Corp <Ntt> 適応ブロック長符号化装置、その方法、プログラム及び記録媒体
JP2007318691A (ja) * 2006-05-29 2007-12-06 Nippon Telegr & Teleph Corp <Ntt> 線形予測モデル次数決定装置、線形予測モデル次数決定方法、そのプログラムおよび記録媒体
JP2008026460A (ja) * 2006-07-19 2008-02-07 Nippon Telegr & Teleph Corp <Ntt> 符号化装置、それらの方法、プログラム及び記録媒体

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110444215A (zh) * 2014-05-01 2019-11-12 日本电信电话株式会社 编码装置、编码方法、程序、记录介质
CN110444216A (zh) * 2014-05-01 2019-11-12 日本电信电话株式会社 解码装置、解码方法、程序、记录介质
CN110444217A (zh) * 2014-05-01 2019-11-12 日本电信电话株式会社 解码装置、解码方法、程序、记录介质
CN110444215B (zh) * 2014-05-01 2022-10-21 日本电信电话株式会社 编码装置、编码方法、记录介质
CN110444216B (zh) * 2014-05-01 2022-10-21 日本电信电话株式会社 解码装置、解码方法、记录介质

Also Published As

Publication number Publication date
JPWO2010067800A1 (ja) 2012-05-17
WO2010067800A1 (ja) 2010-06-17

Similar Documents

Publication Publication Date Title
JP5337235B2 (ja) 符号化方法、復号方法、符号化装置、復号装置、プログラム及び記録媒体
CN103325377B (zh) 音频编码方法
JP5688861B2 (ja) レベル・モードとラン・レングス／レベル・モードの間での符号化を適応させるエントロピー符号化
JPH11143499A (ja) 切替え型予測量子化の改良された方法
JP4866484B2 (ja) パラメータ選択方法、パラメータ選択装置、プログラム及び記録媒体
JP5303074B2 (ja) 符号化方法、復号方法、それらの装置、プログラム及び記録媒体
JP2007510170A (ja) オーディオ符号化における高速コードブック選択方法
JP4848049B2 (ja) 符号化方法、復号方法、それらの装置、プログラム及び記録媒体
JP4918103B2 (ja) 符号化方法、復号方法、それらの装置、プログラム及び記録媒体
JP4834179B2 (ja) 符号化方法、その装置、プログラム及び記録媒体
JP4191503B2 (ja) 音声楽音信号符号化方法、復号化方法、符号化装置、復号化装置、符号化プログラム、および復号化プログラム
JP4091506B2 (ja) ２段音声画像符号化方法、その装置及びプログラム及びこのプログラムを記録した記録媒体
US8576927B2 (en) Encoding method, encoding device, decoding method, decoding device, program, and recording medium
JP5281485B2 (ja) 双方向予測符号化装置、双方向予測復号装置、それらの方法、それらのプログラム及びその記録媒体
JP4918108B2 (ja) 符号化方法、符号化装置、復号方法、復号装置、プログラム及び記録媒体
JP4787889B2 (ja) 符号化方法、符号化装置、復号方法、復号装置、プログラム及び記録媒体
JP2011223615A (ja) 符号化方法、符号化装置、復号方法、復号装置及びプログラム
JP2001175295A (ja) 音声符号化装置
JP2001188592A (ja) 音声符号化方法
JP2001195094A (ja) 音声符号化方法及び音声復号化方法

Legal Events

Date Code Title Description
TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20111004

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20111014

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20141021

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4848049

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350