(57)【要約】      （修正有）【課題】情報信号（ＩＳ）を符号化するための情報エンコーダを提供する。【解決手段】情報エンコーダ（１）は、予測多項式Ａ（ｚ）の線形予測係数を求めるためにＩＳを分析するアナライザ（２）と、予測多項式Ａ（ｚ）から求まる線スペクトル対の多項式対Ｐ（ｚ）とＱ（ｚ）を分析することにより周波数値ｆ１．．．ｆｎを決定し、Ｐ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）からの厳密な虚数スペクトル（ＩＥＳ）を確立しかつＰ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）から導出される厳密な虚数スペクトル（ＩＥＳ）のゼロを識別して、周波数値（ｆ１．．．ｆｎ）を求めるコンバータ（３）と、周波数値（ｆ１．．．ｆｎ）から量子化された周波数値（ｆｑ１．．．ｆｑｎ）を求める量子化器（４）と、量子化された周波数値を含むビットストリームを生成するビットストリーム生成器を備える。【選択図】図１
2019049729.tif 000002
【特許請求の範囲】
【請求項１】
情報信号（ＩＳ）を符号化するための情報エンコーダであって、前記情報エンコーダ（１）は、
予測多項式Ａ（ｚ）の線形予測係数を得るために前記情報信号（ＩＳ）を分析するアナライザ（２）と、
前記予測多項式Ａ（ｚ）の前記線形予測係数を前記予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値ｆ１．．．ｆｎへ変換するコンバータ（３）であって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）、
ｍは前記予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することにより前記周波数値ｆ１．．．ｆｎを決定するように構成され、Ｐ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）からの厳密な虚数スペクトル（ＩＥＳ）を確立することによって、及びＰ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）のゼロを識別することによって、前記周波数値（ｆ１．．．ｆｎ）を得るように構成されたコンバータ（３）と、
前記周波数値（ｆ１．．．ｆｎ）から量子化された周波数値（ｆｑ１．．．ｆｑｎ）を取得する量子化器（４）と、
前記量子化された周波数値（ｆｑ１．．．ｆｑｎ）を含むビットストリームを生成するビットストリーム生成器（５）と、を備えている情報エンコーダ。
【請求項２】
前記コンバータ（３）は、前記予測多項式Ａ（ｚ）から前記多項式Ｐ（ｚ）及びＱ（ｚ）を決定する決定デバイス（６）を備えている請求項１に記載の情報エンコーダ。
【請求項３】
前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）の前記ゼロを識別するためのゼロ識別子（９）を備えている請求項１又は請求項２に記載の情報エンコーダ。
【請求項４】
前記ゼロ識別子（９）は、
ａ）前記実数スペクトル（ＲＥＳ）からヌル(null)周波数で開始すること、
ｂ）前記実数スペクトル（ＲＥＳ）において符号変化が発見されるまで周波数を増大すること、
ｃ）前記虚数スペクトル（ＩＥＳ）においてさらなる符号変化が発見されるまで周波数を増大すること、及び
ｄ）全てのゼロが発見されるまでステップｂ）とステップｃ）を反復すること、によって前記ゼロを識別するように構成されている請求項３に記載の情報エンコーダ。
【請求項５】
前記ゼロ識別子は、前記ゼロを補間(interpolation)によって識別するように構成されている請求項３又は請求項４に記載の情報エンコーダ。
【請求項６】
前記コンバータ（３）は、長くされた多項式対Ｐｅ（ｚ）及びＱｅ（ｚ）を生成するように、前記多項式Ｐ（ｚ）及びＱ（ｚ）に対して値「０」を有する１又は複数の係数を加えるゼロパディング・デバイス（１０）を備えている請求項１から５のいずれか一項に記載の情報エンコーダ。
【請求項７】
前記コンバータ（３）は、前記線形予測係数を前記予測多項式Ａ（ｚ）の前記スペクトル周波数表現（ＲＥＳ、ＩＥＳ）の周波数値（ｆ１．．．ｆｎ）に変換する間に、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）の前記値「０」を有することが知られている係数を用いる演算の少なくとも一部が省略されるように構成されている請求項５又は請求項６に記載の情報エンコーダ。
【請求項８】
前記コンバータ（３）は、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）から合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を確立するように構成された合成多項式フォーマ（１３）を備えている請求項５から７のいずれか一項に記載の情報エンコーダ。
【請求項９】
前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）からの前記厳密な虚数スペクトル（ＩＥＳ）が、前記合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を変換することによって、単一のフーリエ変換により、確立されるように構成されている請求項８に記載の情報エンコーダ。
【請求項１０】
前記コンバータ（３）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１又は複数の多項式を周波数領域へフーリエ変換するフーリエ変換デバイス（８）と、Ｐ（ｚ）から導出される前記スペクトル（ＲＥＳ）の位相をそれが厳密に実数であるように調整するとともに、Ｑ（ｚ）から導出される前記スペクトル（ＩＥＳ）の位相をそれが厳密に虚数であるように調整する調整デバイス（７、１２）と、を備えている請求項１から９のいずれか一項に記載の情報エンコーダ。
【請求項１１】
前記調整デバイス（７、１２）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される前記１若しくは複数の多項式の係数を循環シフトする係数シフタ（７）として構成されている請求項１０に記載の情報エンコーダ。
【請求項１２】
前記係数シフタ（７）は、係数列の原初の中間点が前記数列の第１のポジションへシフトされるような方法で係数を循環シフトするように構成されている請求項１１に記載の情報エンコーダ。
【請求項１３】
前記調整デバイス（７、１２）は、前記フーリエ変換デバイス（８）の出力の位相をシフトする位相シフタ（１２）として構成されている請求項１０に記載の情報エンコーダ。
【請求項１４】
前記位相シフタ（１２）は、ｋ番目の周波数ビンにｅｘｐ（ｉ２πｋｈ／Ｎ）を乗算することによって、前記フーリエ変換デバイス（８）の前記出力の前記位相をシフトするように構成されており、Ｎは、サンプルの長さであり、ｈ＝（ｍ＋ｌ）／２である請求項１３に記載の情報エンコーダ。
【請求項１５】
前記コンバータ（３）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１又は複数の多項式を、Ｐ（ｚ）から導出される前記スペクトル（ＲＥＳ）が厳密に実数であるように、かつＱ（ｚ）から導出される前記スペクトル（ＩＥＳ）が厳密に虚数であるように、ハーフサンプルを用いて周波数領域へフーリエ変換するフーリエ変換デバイス（１４）を備えている請求項１から９のいずれか一項に記載の情報エンコーダ。
【請求項１６】
前記コンバータ（３）は、前記多項式Ｐ（ｚ）とＱ（ｚ）から合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を確立するように構成された合成多項式フォーマ（１３）を備えている請求項１から１５のいずれか一項に記載の情報エンコーダ。
【請求項１７】
前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）からの前記厳密な虚数スペクトル（ＩＥＳ）が、前記合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を変換することによって単一のフーリエ変換により確立されるように構成されている請求項１６に記載の情報エンコーダ。
【請求項１８】
前記コンバータ（３）は、前記多項式Ｐ（ｚ）及びＱ（ｚ）の前記スペクトル（ＲＥＳ、ＩＥＳ）の数値範囲を、前記多項式Ｐ（ｚ）及びＱ（ｚ）又は前記多項式Ｐ（ｚ）及びＱ（ｚ）から導出される１つ若しくは複数の多項式にフィルタ多項式Ｂ（ｚ）を乗算することによって制限するための制限デバイス（１１）を備え、前記フィルタ多項式Ｂ（ｚ）は対称性でありかつ単位円上に根を持たない請求項１から１７のいずれか一項に記載の情報エンコーダ。
【請求項１９】
前記コンバータ（３）は、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）又は前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）から導出される１つ若しくは複数の多項式の前記スペクトル（ＲＥＳ、ＩＥＳ）の数値範囲を、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）にフィルタ多項式Ｂ（ｚ）を乗算することによって制限するための制限デバイス（１１）を備え、前記フィルタ多項式Ｂ（ｚ）は対称性でありかつ単位円上に根を持たない請求項６から１８のいずれか一項に記載の情報エンコーダ。
【請求項２０】
情報信号（ＩＳ）を符号化するための情報エンコーダ（１）を作動するための方法であって、
予測多項式Ａ（ｚ）の線形予測係数を求めるために前記情報信号（ＩＳ）を分析するステップと、
前記予測多項式Ａ（ｚ）の前記線形予測係数を前記予測多項式Ａ（ｚ）のスペクトル周波数表現（ＲＥＳ、ＩＥＳ）の周波数値（ｆ１．．．ｆｎ）へ変換するステップであって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）
ｍは前記予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することによって、前記周波数値（ｆ１．．．ｆｎ）を決定し、Ｐ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）からの厳密な虚数スペクトル（ＩＥＳ）を確立しかつＰ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）のゼロを識別することによって、前記周波数値（ｆ１．．．ｆｎ）を得るステップと、
前記周波数値（ｆ１．．．ｆｎ）から量子化された周波数値（ｆｑ１．．．ｆｑｎ）を求めるステップと、
前記量子化された周波数値（ｆｑ１．．．ｆｑｎ）を含むビットストリーム（ＢＳ）を生成するステップと、を含む方法。
【請求項２１】
プロセッサ上で実行されることにより請求項２０に記載の方法を実行するコンピュータプログラム。
【発明の詳細な説明】
【背景技術】
【０００１】
音声符号化において最も頻繁に使用されるパラダイムは、代数符号励振線形予測（ＡＣＥＬＰ）であり、これはＡＭＲ−ファミリ、Ｇ．７１８やＭＰＥＧ  ＵＳＡＣ［１−３］等のような規格に使用されている。これは、ソースモデルを用いる音声モデリングを基礎としている。そのようなソースモデルは、スペクトル包絡線をモデリングするための線形予測子（ＬＰ）、基本周波数をモデリングするための長時間予測子（ＬＴＰ）及び残余（residual）のための代数コードブックにより構成されている。
【０００２】
線形予測モデルの係数は、量子化に対する感受性が極めて高く、よって通常これらは、量子化に先立って、まず線スペクトル周波数（ＬＳＦ）又はイミッタンススペクトル周波数（ＩＳＦ）に変換される。ＬＳＦ／ＩＳＦ領域は、量子化誤差に対して堅牢（robust）であり、これらの領域では、予測子の安定性を容易に保全することができ、よって、量子化のための適切な領域が提供される［４］。
【０００３】
以下、周波数値と称するＬＳＦ／ＩＳＦは、次のように、次数ｍの線形予測多項式Ａ（ｚ）から求めることができる。線スペクトル対の多項式は、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）        （１）
と定義される。ここで、線スペクトル対表現ではｌ＝１であり、かつイミッタンススペクトル対表現ではｌ＝０であるが、原則的には、任意のｌ≧０が有効である。したがって、以後、単にｌ≧０であるものと仮定する。
【０００４】
元の予測子は、Ａ（ｚ）＝１／２［Ｐ（ｚ）＋Ｑ（ｚ）］を用いていつでも復元できることに留意されたい。したがって、多項式Ｐ（ｚ）及びＱ（ｚ）は、Ａ（ｚ）の全情報を包含する。
【０００５】
ＬＳＰ／ＩＳＰ多項式の中心的性質は、Ａ（ｚ）がその全ての根を単位円内に有する場合に限り、Ｐ（ｚ）及びＱ（ｚ）の根が単位円上でインタレースされる、というものである。Ｐ（ｚ）及びＱ（ｚ）の根が単位円上に存在するため、それらを角度のみで表現することができる。これらの角度は、周波数に対応し、またＰ（ｚ）及びＱ（ｚ）のスペクトルは、根に対応する周波数においてその対数振幅スペクトルにおける垂線を有することから、これらの根を周波数値と称する。
【０００６】
周波数値は、結果的に、予測子Ａ（ｚ）の全情報を符号化することになる。さらに、周波数値のうちの１つにおける僅かな誤差がもたらす、対応する周波数の近くにスペクトルが位置決めされる復元された予測子のスペクトル誤差が小さくなるように、周波数値は、量子化誤差に対して堅牢（robust）であることが分かっている。これらの好ましい性質に起因して、ＬＳＦ又はＩＳＦ領域における量子化は、全てのメインストリーム音声コーデックにおいて使用されている［１−３］。
【０００７】
しかしながら、周波数値を用いる場合の課題の１つは、多項式Ｐ（ｚ）及びＱ（ｚ）の係数からその位置を効率的に求めることにある。結局のところ、多項式の根を求めることは、古典的かつ困難な課題である。このタスクに関してこれまでに提案されている方法は、下記の手法を含んでいる：
・  初期の手法のうちの１つは、ゼロが単位円上に存在し、よって、これが振幅スペクトルにおいてゼロとして現出する、ということを用いている［５］。したがって、Ｐ（ｚ）とＱ（ｚ）の係数の離散型フーリエ変換を利用することにより、振幅スペクトルの谷を探索することができる。各谷は、根の位置を示し、スペクトルが充分にアップサンプリングされていれば、全ての根を求めることができる。しかしながら、この方法では、谷の位置からの正確なポジションの決定が困難であるため、近似ポジションしか得られない。
・  最も高頻度で使用される手法は、チェビシェフ多項式を基礎とするものであり、［６］に提示されている。これは、多項式Ｐ（ｚ）及びＱ（ｚ）はそれぞれが対称かつ逆対称であり、それによって多くの冗長な情報を含む、という認識に依存する。ｚ＝±１における自明な零点を除去することにより、かつ置換ｘ＝ｚ＋ｚ−１（チェビシェフ変換として知られる）を用いれば、多項式を代替表現ＦＰ（ｘ）及びＦＱ（ｘ）へ変換することができる。これらの多項式の次数は、Ｐ（ｚ）とＱ（ｚ）の半分であり、よってこれらは、−２から＋２までの範囲にのみ実根を有する。多項式ＦＰ（ｘ）及びＦＱ（ｘ）は、ｘが実数のとき実数値であることに留意されたい。さらに、これらの根は、単根であることから、ＦＰ（ｘ）とＦＱ（ｘ）は、その各根において零交差を有する。
【０００８】
ＡＭＲ−ＷＢ等の音声コーデックにおいて、この手法は、実軸上の固定格子上で多項式ＦＰ（ｘ）及びＦＱ（ｘ）を評価して全ての零交差を求めるように適用される。根位置は、零交差の周囲での線形補間によってさらに精緻化される。この手法の優位点は、冗長係数の省略により複雑さが低減されることにある。
【０００９】
上述の方法は、既存のコーデックにおいて充分に機能するが、幾つかの問題点もある。
【先行技術文献】
【非特許文献】
【００１０】
【非特許文献１】[1] B.  Bessette, R. Salami, R. Lefebvre, M. Jelinek, J. Rotola-Pukkila, J. Vainio, H. Mikkola, and K. Jaervinen, “The adaptive multirate wideband speech codec (AMR-WB)”, Speech and Audio Processing, IEEE Transac- tions on, vol. 10, no. 8, pp. 620-636, 2002.
【非特許文献２】[2] ITU-T G.718, “Frame error robust narrow-band and wideband embedded variable bit-rate coding of speech and audio from 8-32 kbit/s”, 2008.
【非特許文献３】[3] M. Neuendorf, P. Gournay, M. Multrus, J. Lecomte, B. Bessette, R. Geiger, S. Bayer, G. Fuchs, J. Hilpert, N. Rettelbach, R. Salami, G. Schuller, R. Lefebvre, and B. Grill, “Unified speech and audio coding scheme for high quality at low bitrates”, in Acoustics, Speech and Signal Processing. ICASSP 2009. IEEE Int Conf, 2009, pp. 1-4.
【非特許文献４】[4] T. Baeckstroem and C. Magi, “Properties of line spectrum pair polynomials - a review”, Signal Processing, vol. 86, no. 11, pp. 3286-3298, November 2006.
【非特許文献５】[5] G. Kang and L. Fransen, “Application of line-spectrum pairs to low-bit- rate speech encoders”, in Acoustics, Speech, and Signal Processing, IEEE International Conference on ICASSP’85., vol. 10. IEEE, 1985, pp. 244-247.
【非特許文献６】[6] P. Kabal and R. P. Ramachandran, “The computation of line spectral frequencies using Chebyshev polynomials”, Acoustics, Speech and Signal Processing, IEEE Transactions on, vol. 34, no. 6, pp. 1419-1426, 1986.
【非特許文献７】[7] 3GPP TS 26.190 V7.0.0, “Adaptive multi-rate (AMR-WB) speech codec”, 2007.
【非特許文献８】[8] T. Baeckstroem, C. Magi, and P. Alku, “Minimum separation of line spec- tral frequencies”, IEEE Signal Process. Lett., vol. 14, no. 2, pp. 145-147, February 2007.
【非特許文献９】[9] T. Baeckstroem, “Vandermonde factorization of Toeplitz matrices and applications in filtering and warping,” IEEE Trans. Signal Process., vol. 61, no. 24, pp. 6257-6263, 2013.
【非特許文献１０】[10] V. F. Pisarenko, “The retrieval of harmonics from a covariance function”, Geophysical Journal of the Royal Astronomical Society, vol. 33, no. 3, pp. 347-366, 1973.
【非特許文献１１】[11] E. Durand, Solutions Numeriques des Equations Algebriques. Paris: Masson, 1960.
【非特許文献１２】[12] I. Kerner, “Ein Gesamtschrittverfahren zur Berechnung der Nullstellen von Polynomen”, Numerische Mathematik, vol. 8, no. 3, pp. 290-294, May 1966.
【非特許文献１３】[13] O. Aberth, “Iteration methods for finding all zeros of a polynomial simultaneously”, Mathematics of Computation, vol. 27, no. 122, pp. 339-344, April 1973.
【非特許文献１４】[14] L. Ehrlich, “A modified newton method for polynomials”, Communications of the ACM, vol. 10, no. 2, pp. 107-108, February 1967.
【非特許文献１５】[15] D. Starer and A. Nehorai, “Polynomial factorization algorithms for adaptive root estimation”, in Int. Conf. on Acoustics, Speech, and Signal Processing, vol. 2.  Glasgow, UK: IEEE, May 1989, pp. 1158-1161.
【非特許文献１６】[16] --, “Adaptive polynomial factorization by coefficient matching”, IEEE Transactions on Signal Processing, vol. 39, no. 2, pp. 527-530, February 1991.
【非特許文献１７】[17] G. H. Golub and C. F. van Loan, Matrix Computations, 3rd ed. John Hopkins University Press, 1996.
【非特許文献１８】[18] T. Saramaeki, “Finite impulse response filter design”, Handbook for Digital Signal Processing, pp. 155-277, 1993.
【発明の概要】
【発明が解決しようとする課題】
【００１１】
解決すべき課題は、情報を符号化するための改良された概念を提供することである。
【課題を解決するための手段】
【００１２】
第１の態様において、この課題は、情報信号を符号化するための情報エンコーダによって解決される。情報エンコーダは、
予測多項式Ａ（ｚ）の線形予測係数を得るために情報信号を分析するアナライザと、
予測多項式Ａ（ｚ）の線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値へ変換するコンバータであって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）、
ｍは予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することにより前記周波数値を決定するように構成され、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）からの厳密な虚数スペクトルを確立することによって、及びＰ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）から導出される厳密な虚数スペクトルのゼロを識別することによって、前記周波数値を得るように構成されたコンバータと、
前記周波数値から量子化された周波数値を取得する量子化器と、
前記量子化された周波数値を含むビットストリームを生成するビットストリーム生成器と、を備える。
【００１３】
本発明による情報エンコーダは零交差探索を用いるのに対し、従来技術による根を求めるためのスペクトル手法は、振幅スペクトル中の谷の検知に依存する。しかしながら、谷を探索する場合の確度は、零交差を探索する場合より劣る。例えば、数列［４，２，１，２，３］を考察されたい。明らかに、最小値は第３の要素であり、よって、ゼロは、第２の要素と第４の要素との間のどこかに存在することになる。言い替えれば、ゼロが第３の要素の右側に存在するか、左側に存在するかを決定することはできない。しかしながら、数列［４，２，１，−２，−３］について考察する場合は、直ちに、零交差が第３の要素と第４の要素との間に存在することが分かり、よって、許容誤差が半減する。振幅スペクトル手法の場合、零交差探索の場合と同じ確度を達成するには、分析点の数を２倍にする必要があることになる。
【００１４】
振幅│Ｐ（ｚ）│と│Ｑ（ｚ）│の評価を比較すれば、零交差手法は、確度において著しい優位点を有する。例えば、数列３，２，−１，−２を考察されたい。零交差手法では、ゼロは、２と−１との間に存在することが明らかである。しかしながら、対応する振幅列３，２，１，２を考察する場合、ゼロは、第２の要素と最後の要素との間のどこかに存在する、としか結論することができない。言い替えれば、零交差手法の場合の確度は、振幅ベースの手法の２倍である。
【００１５】
さらに、本発明による情報エンコーダは、ｍ＝１２８等の長い予測子を用いてもよい。これに対して、チェビシェフ変換は、Ａ（ｚ）の長さが比較的短い、例えばｍ≦２０である場合にしか十分に機能しない。長い予測子の場合、チェビシェフ変換は、数値的に不安定であり、よって、アルゴリズムの実用的実装は不可能である。
【００１６】
したがって、提案する情報エンコーダの主要な特性は、零交差が探索されることに起因して、そして時間領域から周波数領域への変換が行われ、それによって極めて低い計算複雑性によりゼロを発見し得ることに起因して、チェビシェフベースの方法のように高い確度、又はより優れた確度を達成し得ることである。
【００１７】
結果的に、本発明による情報エンコーダは、ゼロ（根）をより正確に決定するだけでなく、低い計算複雑性で決定する。
【００１８】
本発明による情報エンコーダは、数列の線スペクトルを決定する必要があるあらゆる信号処理アプリケーションにおいて使用可能である。本明細書では、例示的に、情報エンコーダを音声符号化の側面で論じる。本発明は、音声信号又は汎用オーディオ信号等の入力信号から線スペクトルを決定するための方法を必要とし、かつ入力信号がデジタルフィルタ又は他の数列として表現される、スペクトル振幅包絡線、知覚周波数マスキングしきい値、時間振幅包絡線、知覚時間マスキングしきい値、又は他の包絡線形状、若しくは線スペクトルを用いて包絡線情報を表現する、符号化、分析又は処理のための自己相関信号等の包絡線形状と等価の他の表現、をモデリングするための線形予測子を使用する、音声、オーディオ及び／又はビデオ符号化デバイス又はアプリケーションにおいて適用可能である。
【００１９】
情報信号は、例えば、オーディオ信号又はビデオ信号であってもよい。周波数値は、線スペクトル周波数であっても、イミッタンススペクトル周波数であってもよい。ビットストリーム内で伝送される量子化された周波数値は、オーディオ信号又はビデオ信号を再現するために、デコーダがビットストリームを復号できるようにする。
【００２０】
本発明の好適な一実施形態によれば、コンバータは、予測多項式Ａ（ｚ）から多項式Ｐ（ｚ）とＱ（ｚ）を決定するための決定デバイスを備えている。
【００２１】
本発明の好適な実施形態によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）から導出される厳密な虚数スペクトルのゼロを識別するためのゼロ識別子を備えている。
【００２２】
本発明の好適な一実施形態によれば、ゼロ識別子は、
ａ）実数スペクトルをヌル周波数において開始することと、
ｂ）実数スペクトルにおいて符号変化が発見されるまで周波数を増大することと、
ｃ）虚数スペクトルにおいてさらなる符号変化が発見されるまで周波数を増大することと、
ｄ）全てのゼロが発見されるまでステップｂ）及びステップｃ）を反復すること、によってゼロを識別するように構成されている。
【００２３】
Ｑ（ｚ）ひいてはスペクトルの虚数部は、ヌル周波数において常にゼロを有することに留意されたい。根は重なり合うことから、Ｐ（ｚ）ひいてはスペクトルの実数部は、ヌル周波数において常に非ゼロになる。したがって、ヌル周波数における実数部から開始して、最初の零交差ひいては最初の周波数値を示す最初の符号変化が発見されるまで、周波数を増大させることができる。
【００２４】
根はインタレースされるので、Ｑ（ｚ）のスペクトルは次の符号変化を有することとなる。したがって、Ｑ（ｚ）のスペクトルの符号変化が発見されるまで、周波数を増大させることができる。次に、このプロセスは、全ての周波数値が発見されるまで、Ｐ（ｚ）及びＱ（ｚ）のスペクトル間で交互に反復されてもよい。したがって、スペクトルにおける零交差の位置決めに使用される手法は、チェビシェフ領域において適用される手法に類似する［６、７］。
【００２５】
Ｐ（ｚ）及びＱ（ｚ）のゼロはインタレースされることから、１パスで全てのゼロを発見するように実数部及び複素数部上のゼロの探索を交互に行い、よって、複雑性を完全探索の半分に減らすことができる。
【００２６】
本発明の好適な一実施形態によれば、ゼロ識別子は、ゼロを補間によって識別するように構成される。
【００２７】
零交差手法に加えて、ゼロのポジションを、例えば［７］のような従来方法において行われるようなより高い確度で推定できるように、補間を容易に適用することができる。
【００２８】
本発明の好適な一実施形態によれば、コンバータは、長くされた多項式対Ｐｅ（ｚ）及びＱｅ（ｚ）を生成するための多項式Ｐ（ｚ）及びＱ（ｚ）へ値「０」を有する１つ又は複数の係数を加えるゼロパディング・デバイスを備えている。確度は、評価されるスペクトルの長さを拡大することによってさらに向上させることができる。システムに関する情報を基礎として、事例によっては、実質上、周波数値間の最小距離を決定し、ひいては、全ての周波数値をそれにより発見することができるスペクトルの最小長さを決定することが可能である［８］。
【００２９】
本発明の好適な一実施形態によれば、コンバータは、線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値へ変換する間に、長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）の値「０」を有することが知られている係数を用いる演算の少なくとも一部が省略されるように構成される。
【００３０】
しかしながら、スペクトルの長さが増せば、計算複雑性も増大する。複雑性の最大の誘因は、Ａ（ｚ）の係数の、時間領域から周波数領域への高速フーリエ変換などの変換である。しかしながら、これは、係数ベクトルが所望される長さにゼロパディングされていることから、極めて疎である。複雑性の低減には、この事実を容易に使用することができる。これは、どの係数がゼロであるかが正確に分かり、よって、高速フーリエ変換の各反復において単純にゼロを包含する演算を省くことができるという意味において、どちらかといえば単純な問題である。このように疎な高速フーリエ変換の適用は、単純であり、よって当業者たるプログラマは誰でもこれを実装することができる。このような実装の複雑性は、Ｏ（Ｎｌｏｇ２（１＋ｍ＋ｌ））であり、ここで、Ｎは、スペクトルの長さであり、ｍ及びｌは、先に定義した通りである。
【００３１】
本発明の好適な一実施形態によれば、コンバータは、長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）から合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を確立するように構成された合成多項式フォーマを備えている。
【００３２】
本発明の好適な一実施形態によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）からの厳密な虚数スペクトルが、合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を変換することによって単一のフーリエ変換により、確立されるように構成される。
【００３３】
本発明の好適な一実施形態によれば、コンバータは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式を周波数領域へフーリエ変換するフーリエ変換デバイスと、Ｐ（ｚ）から導出されるスペクトルの位相を、それが厳密に実数であるように調整するとともに、Ｑ（ｚ）から導出されるスペクトルの位相を、それが厳密に虚数であるように調整する調整デバイスと、を備えている。フーリエ変換デバイスは、高速フーリエ変換を基礎とするものであっても、離散型フーリエ変換を基礎とするものであってもよい。
【００３４】
本発明の好適な一実施形態によれば、調整デバイスは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式の係数を循環シフトする係数シフタとして構成されている。
【００３５】
本発明の好適な一実施形態によれば、係数シフタは、係数列の原初の中間点がその数列の第１のポジションへシフトされるような方法で、係数を循環シフトするように構成されている。
【００３６】
理論的には、対称数列のフーリエ変換が実数値をとり、かつ逆対称の数列が純虚数フーリエスペクトルを有することは周知である。本事例における入力数列は、長さｍ＋ｌの多項式Ｐ（ｚ）又はＱ（ｚ）の係数であるが、これより遙かに長い長さＮ＞＞（ｍ＋ｌ）の離散型フーリエ変換を有する方が好まれると思われる。より長いフーリエスペクトルを生成する従来的手法は、入力信号のゼロパディングである。しかしながら、数列のゼロパディングは、対称性が保たれるように慎重に実装されなければならない。
【００３７】
まず、係数、
［ｐ０，ｐ１，ｐ２，ｐ１，ｐ０］
を有する多項式Ｐ（ｚ）について考察する。
【００３８】
通常、ＦＦＴアルゴリズムを適用する方法は、対称点が第１の要素であることを要求し、よって、例えば、ＭＡＴＬＡＢに適用する場合、
ｆｆｔ（［ｐ２，ｐ１，ｐ０，ｐ０，ｐ１］）
のように書いて実数値の出力を求めることができる。具体的には、循環シフトが適用されてもよく、よって、中間点要素に対応する対称点、即ち係数ｐ２は、第１のポジションになるように左へシフトされる。次に、ｐ２の左側にあった係数は、数列の終わりに付加される。
【００３９】
ゼロパディングされた数列、
［ｐ０，ｐ１，ｐ２，ｐ１，ｐ０，０，０．．．０］
の場合も、同じプロセスを適用することができる。よって、数列、
［ｐ２，ｐ１，ｐ０，０，０．．．０，ｐ０，ｐ１］
が実数値の離散型フーリエ変換を有する。ここで、入力数列におけるゼロの数は、スペクトルの所望される長さをＮとすれば、Ｎ−ｍ−ｌである。
【００４０】
同様に、多項式Ｑ（ｚ）に対応する係数、
［ｑ０，ｑ１，０，−ｑ１，−ｑ０］
について考察されたい。前述の中間点が第１のポジションにくるように循環シフトを適用すれば、
［０，−ｑ１，−ｑ０，ｑ０，ｑ１］
が得られ、これは、純虚数の離散型フーリエ変換を有する。次には、この数列にゼロパディング変換を行なうことができる。
［０，−ｑ１，−ｑ０，０，０．．．０，ｑ０，ｑ１］
【００４１】
上記は、数列の長さが奇数である事例にのみ当てはまり、よって、ｍ＋ｌは偶数であることに留意されたい。ｍ＋ｌが奇数である事例に関しては、２つのオプションがある。即ち、周波数領域において循環シフトを実装することができ、あるいは、ハーフサンプルでＤＦＴを適用することができる（下記参照）。
【００４２】
本発明の好適な一実施形態によれば、調整デバイスは、フーリエ変換デバイスの出力の位相をシフトするための位相シフタとして構成される。
【００４３】
本発明の好適な一実施形態によれば、位相シフタは、フーリエ変換デバイスの出力の位相を、ｋ番目の周波数ビンにｅｘｐ（ｉ２πｋｈ／Ｎ）を乗算することによってシフトするように構成される。ここで、Ｎは、サンプルの長さであり、かつｈ＝（ｍ＋ｌ）／２である。
【００４４】
時間領域における循環シフトが周波数領域における位相回転と同一であることはよく知られている。具体的には、時間領域におけるｈ＝（ｍ＋ｌ）／２ステップのシフトは、ｋ番目の周波数ビンとｅｘｐ（−ｉ２πｋｈ／Ｎ）との乗算に一致する。ここで、Ｎはスペクトルの長さである。したがって、循環シフトの代わりに、周波数領域における乗算を適用すれば、全く同じ結果を得ることができる。ただし、この手法では、複雑性が僅かに増す。ｈ＝（ｍ＋ｌ）／２は、ｍ＋ｌが偶数であるときにのみ整数であることに留意されたい。ｍ＋ｌが奇数のとき、循環シフトは、有理数のステップ分の遅延を必要とすることになり、これを直に実装することは困難である。代わりに、上述の相回転によって、周波数領域における対応するシフトを適用することができる。
【００４５】
本発明の好適な実施形態によれば、コンバータは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式を、Ｐ（ｚ）から導出されるスペクトルが厳密に実数であるように、かつＱ（ｚ）から導出されるスペクトルが厳密に虚数であるように、ハーフサンプルを用いて周波数領域へフーリエ変換するフーリエ変換デバイスを備えている。
【００４６】
ある代替例は、ハーフサンプルを有するＤＦＴを実装するものである。具体的には、従来のＤＦＴを、
2019049729.tif 000003
で定義できるのに対して、ハーフサンプルＤＦＴは、次式のように定義することができる。
2019049729.tif 000004
【００４７】
この公式化に対しては、ＦＦＴのような高速実装を容易に考案することができる。
【００４８】
この公式化の利点は、いまや対称点が通常のｎ＝１ではなく、ｎ＝１／２に存在することにある。よって、このハーフサンプルＤＦＴを用いれば、数列、
［２，１，０，０，１，２］
で、実数値のフーリエスペクトルが得られる。
【００４９】
よって、ｍ＋ｌが奇数である場合、係数ｐ０，ｐ１，ｐ２，ｐ２，ｐ１，ｐ０を有する多項式Ｐ（ｚ）に関して、ハーフサンプルＤＦＴ及びゼロパディングを用いれば、入力数列が、
［ｐ２，ｐ１，ｐ０，０，０．．．０，ｐ０，ｐ１，ｐ２］
のとき、実数値のスペクトルを求めることができる。
【００５０】
同様に、多項式Ｑ（ｚ）の場合、数列、
［−ｑ２，−ｑ１，−ｑ０，０，０．．．０，ｑ０，ｑ１，ｑ２］
にハーフサンプルＤＦＴを適用して、純虚数スペクトルを求めることができる。
【００５１】
これらの方法により、ｍとｌの如何なる組合せに対しても、多項式Ｐ（ｚ）の実数値スペクトルと任意のＱ（ｚ）の純虚数スペクトルを求めることができる。実際には、Ｐ（ｚ）とＱ（ｚ）のスペクトルが各々純実数と純虚数であることから、これらを単一の複素スペクトルに格納することができ、よって、これは、Ｐ（ｚ）＋Ｑ（ｚ）＝２Ａ（ｚ）のスペクトルに一致する。係数２によるスケーリングでは、根の位置に変化がなく、よってこれを無視することができる。したがって、単一のＦＦＴを用いてＡ（ｚ）のスペクトルのみを評価すれば、Ｐ（ｚ）とＱ（ｚ）のスペクトルを求めることができる。先に説明したように、Ａ（ｚ）の係数に循環シフトを適用するだけでよい。
【００５２】
例えば、ｍ＝４でありｌ＝０である場合、Ａ（ｚ）の係数は、
［ａ０，ａ１，ａ２，ａ３，ａ４］
であり、これは、ゼロパディングによって、次のように任意の長さＮにすることができる。
［ａ０，ａ１，ａ２，ａ３，ａ４，０，０．．．０］
【００５３】
次に、（ｍ＋ｌ）／２＝２ステップの循環シフトを適用すれば、次のようになる。
［ａ２，ａ３，ａ４，０，０．．．０，ａ０，ａ１］
【００５４】
この数列のＤＦＴを採用すれば、スペクトルの実数部及び複素数部にＰ（ｚ）とＱ（ｚ）のスペクトルが存在する。
【００５５】
本発明の好適な一実施形態によれば、コンバータは、多項式Ｐ（ｚ）とＱ（ｚ）から合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を確立するように構成された合成多項式フォーマを備えている。
【００５６】
本発明の好適な一実施形態によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）からの厳密な虚数スペクトルが、例えば高速フーリエ変換（ＦＦＴ）により合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を変換することによって、単一のフーリエ変換により確立されるように構成されている。
【００５７】
多項式Ｐ（ｚ）とＱ（ｚ）は、各々、ｚ−（ｍ＋ｌ）／２における対称軸と対称でありかつ逆対称である。よって、単位円ｚ＝ｅｘｐ（ｉθ）上で評価されるｚ−（ｍ＋ｌ）／２Ｐ（ｚ）とｚ−（ｍ＋ｌ）／２Ｑ（ｚ）のスペクトルは、各々、実数値及び複素数値ということになる。ゼロは、単位円上に存在することから、零交差を探索すれば、ゼロを発見することができる。さらに、単位円上の評価は、単純に高速フーリエ変換によって実装することができる。
【００５８】
ｚ−（ｍ＋ｌ）／２Ｐ（ｚ）とｚ−（ｍ＋ｌ）／２Ｑ（ｚ）に対応するスペクトルは、各々実数及び複素数であることから、これらを単一の高速フーリエ変換によって実装することができる。具体的には、和ｚ−（ｍ＋ｌ）／２（Ｐ（ｚ）＋Ｑ（ｚ））を求めれば、スペクトルの実数部と複素数部は、各々、ｚ−（ｍ＋ｌ）／２Ｐ（ｚ）とｚ−（ｍ＋ｌ）／２Ｑ（ｚ）に対応する。さらに、
ｚ−（ｍ＋ｌ）／２（Ｐ（ｚ）＋Ｑ（ｚ））＝２ｚ−（ｍ＋ｌ）／２Ａ（ｚ）        （４）
であることから、Ｐ（ｚ）とＱ（ｚ）を明示的に決定することなく、２ｚ−（ｍ＋ｌ）／２Ａ（ｚ）のＦＦＴを直に取り込んでｚ−（ｍ＋ｌ）／２Ｐ（ｚ）とｚ−（ｍ＋ｌ）／２Ｑ（ｚ）に対応するスペクトルを求めることができる。知りたいものはゼロの位置だけであるため、スカラ２による乗算を省略し、代わりにＦＦＴによるｚ−（ｍ＋ｌ）／２Ａ（ｚ）を評価することができる。Ａ（ｚ）は、ｍ＋１個の非ゼロ係数しか持たないことから、ＦＦＴ枝刈りを用いて複雑性を低減できることを観察されたい［１１］。全ての根が発見されることを保証するためには、２つのゼロ毎にその間の少なくとも１つの周波数上でスペクトルが評価されるに足る長い長さＮのＦＦＴを用いなければならない。
【００５９】
本発明の好適な一実施形態によれば、コンバータは、多項式Ｐ（ｚ）とＱ（ｚ）のスペクトルの数値範囲を、多項式Ｐ（ｚ）とＱ（ｚ）又は多項式Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイスを備えている。ここで、フィルタ多項式Ｂ（ｚ）は、対称性であって、単位円上に根を持たない。
【００６０】
音声コーデックは、限定リソースを有するモバイルデバイス上で実装されることが多く、よって、数値演算は、固定小数点表示によって実装されなければならない。したがって、実装されるアルゴリズムは、その範囲が限定される数値表現を用いて動作することが極めて重要である。しかしながら、一般的な音声スペクトル包絡線の場合、フーリエスペクトルの数値範囲があまりに大きく、よって、零交差位置の保持を保証するためには、ＦＦＴの３２ビット実装が必要である。
【００６１】
一方で、１６ビットのＦＦＴは、より低い複雑性で実装可能である場合が多く、よって、この１６ビット範囲内に適合するようにスペクトル値の範囲を限定することが有益であると考えられる。方程式｜Ｐ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜と｜Ｑ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜から、Ｂ（ｚ）Ａ（ｚ）の数値範囲を限定することによってＢ（ｚ）Ｐ（ｚ）とＢ（ｚ）Ｑ（ｚ）の数値範囲も限定されることが分かる。Ｂ（ｚ）が単位円上にゼロを持たなければ、Ｂ（ｚ）Ｐ（ｚ）及びＢ（ｚ）Ｑ（ｚ）は、単位円上にＰ（ｚ）及びＱ（ｚ）と同じ零交差を有する。さらに、Ｂ（ｚ）は、ｚ−（ｍ＋ｌ＋ｎ）／２Ｐ（ｚ）Ｂ（ｚ）及びｚ−（ｍ＋ｌ＋ｎ）／２Ｑ（ｚ）Ｂ（ｚ）がそれぞれ対称性と逆対称を保ち、かつそのスペクトルが純実数及び純虚数であるように、対称性でなければならない。したがって、ｚ（ｎ＋ｌ）／２Ａ（ｚ）のスペクトルを評価する代わりに、ｚ（ｎ＋ｌ＋ｎ）／２Ａ（ｚ）Ｂ（ｚ）を評価することができる。ここで、Ｂ（ｚ）は、単位円上に根を持たない次数ｎの対称多項式である。言い替えれば、先に述べたものと同じ手法を適用することができるが、まずは、Ａ（ｚ）をフィルタＢ（ｚ）で乗算し、次に修正した位相シフトｚ−（ｍ＋ｌ＋ｎ）／２を適用する。
【００６２】
残りのタスクは、Ａ（ｚ）Ｂ（ｚ）の数値範囲が、Ｂ（ｚ）は必ず対称性でありかつ単位円上に根を持たないという制約によって限定されるように、フィルタＢ（ｚ）を設計することである。この要件を満たす最も単純なフィルタは、次数２の線形位相フィルタ、
Ｂ１（ｚ）＝β０＋β１ｚ−１＋β２ｚ−２          （５）
である。ここで、βｋ∈Ｒはパラメータであり、かつ│β２│＞２│β１│である。βｋを調整すれば、スペクトル傾斜を修正することができ、積Ａ（ｚ）Ｂ１（ｚ）の数値範囲を縮小することができる。計算上極めて効率的な手法は、βを、０周波数及びナイキストにおける振幅が等しくなるように、即ち│Ａ（１）Ｂ１（１）│＝│Ａ（−１）Ｂ１（−１）│であるように選択することであり、例えば、次式を選ぶことができる。
β０＝Ａ（１）−Ａ（−１）、
β１＝２（Ａ（１）＋Ａ（−１））      （６）
【００６３】
この手法は、略平坦なスペクトルを提供する。
【００６４】
Ａ（ｚ）が高域通過特性を有するのに対して、Ｂ１（ｚ）は低域通過であり（図５も参照されたい）、積Ａ（ｚ）Ｂ１（ｚ）は、予期されるように、０周波数及びナイキスト周波数において同じ振幅を有し、かつこれは、多かれ少なかれ平坦である。Ｂ１（ｚ）は１自由度しか持たないので、積が完全に平坦になると予期できないことは明らかである。さらに、Ｂ１（ｚ）Ａ（ｚ）の最高ピークと最低の谷との比は、Ａ（ｚ）のそれより遙かに小さいものであり得ることを観察されたい。これは、Ｂ１（ｚ）Ａ（ｚ）の数値範囲がＡ（ｚ）のそれより遙かに小さいという望ましい効果が達成されていることを意味する。
【００６５】
第２の、僅かに複雑性の高い方法は、Ａ（０．５ｚ）のインパルス応答の自己相関ｒｋを計算するものである。ここでは、０．５を乗算することによってＡ（ｚ）のゼロが原点方向へ移動し、それによってスペクトル振幅が約半分に縮小される。自己相関ｒｋにレビンソン−ダービンを適用すれば、最小位相である次数ｎのフィルタＨ（ｚ）が得られる。次に、Ｂ２（ｚ）＝ｚ−ｎＨ（ｚ）Ｈ（ｚ−１）を定義して、略一定である│Ｂ２（ｚ）Ａ（ｚ）│を求めることができる。│Ｂ２（ｚ）Ａ（ｚ）│の範囲が│Ｂ１（ｚ）Ａ（ｚ）│のそれより小さいことに留意され得る。Ｂ（ｚ）を設計するためのさらなる手法は、ＦＩＲ設計に関する古典的文献［１８］において容易に見出すことができる。
【００６６】
本発明の好適な一実施形態によれば、コンバータは、長くされた多項式Ｐｅ（ｚ）とＱｅ（ｚ）又は長くされた多項式Ｐｅ（ｚ）とＱｅ（ｚ）から導出される１つ又は複数の多項式のスペクトルの数値範囲を、長くされた多項式Ｐｅ（ｚ）とＱｅ（ｚ）にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイスを備えている。ここで、フィルタ多項式Ｂ（ｚ）は、対称性であって、単位円上に根を持たない。Ｂ（ｚ）は、先に説明した通りに求めることができる。
【００６７】
さらなる態様では、この課題は、情報信号を符号化するための情報エンコーダを作動するための方法によって解決される。本方法は、
予測多項式Ａ（ｚ）の線形予測係数を求めるべく情報信号を分析するステップと、
予測多項式Ａ（ｚ）の線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値ｆ１．．．ｆｎへ変換するステップであって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）
ｍは予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
のように定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することによって前記周波数値ｆ１．．．ｆｎを決定し、、Ｐ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）からの厳密な虚数スペクトルを確立しかつＰ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）から導出される厳密な虚数スペクトルのゼロを識別することによって、前記周波数値ｆ１．．．ｆｎを得るステップと、
前記周波数値ｆ１．．．ｆｎから量子化された周波数値ｆｑ１．．．ｆｑｎを得るステップと、
前記量子化された周波数値ｆｑ１．．．ｆｑｎを含むビットストリームを生成するステップと、を含む。
【００６８】
さらに、本プログラムは、プロセッサ上で実行されることにより本発明による方法を実行するコンピュータプログラムによっても注目される。
【図面の簡単な説明】
【００６９】
【図１】図１は、本発明による情報エンコーダの一実施形態を示す略図である。
【図２】図２は、Ａ（ｚ）、Ｐ（ｚ）及びＱ（ｚ）の例示的な関係を示す。
【図３】図３は、本発明による情報エンコーダのコンバータの第１の実施形態を示す略図である。
【図４】図４は、本発明による情報エンコーダのコンバータの第２の実施形態を示す略図である。
【図５】図５は、予測子Ａ（ｚ）、対応する平坦化フィルタＢ１（ｚ）及びＢ２（ｚ）及び積Ａ（ｚ）Ｂ１（ｚ）及びＡ（ｚ）Ｂ２（ｚ）の例示的な振幅スペクトルを示す。
【図６】図６は、本発明による情報エンコーダのコンバータの第３の実施形態を示す略図である。
【図７】図７は、本発明による情報エンコーダのコンバータの第４の実施形態を示す略図である。
【図８】図８は、本発明による情報エンコーダのコンバータの第５の実施形態を示す略図である。
【発明を実施するための形態】
【００７０】
続いて、添付の図面を参照し、本発明の好適な実施形態について論じる。
【００７１】
図１は、本発明による情報エンコーダ１の一実施形態を示す略図である。
【００７２】
情報信号ＩＳを符号化するための情報エンコーダ１は、
予測多項式Ａ（ｚ）の線形予測係数を求めるべく情報信号ＩＳを分析するアナライザ２と、
予測多項式Ａ（ｚ）の線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現ＲＥＳ、ＩＥＳの周波数値ｆ１．．．ｆｎへ変換するコンバータ３であって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）、
ｍは予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
のように定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することにより、前記周波数値ｆ１．．．ｆｎを決定するように構成され、前記コンバータ３は、Ｐ（ｚ）から導出される厳密な実数スペクトルＲＥＳ及びＱ（ｚ）からの厳密な虚数スペクトルＩＥＳを確立しかつＰ（ｚ）から導出される厳密な実数スペクトルＲＥＳ及びＱ（ｚ）から導出される厳密な虚数スペクトルＩＥＳのゼロを識別することによって、前記周波数値ｆ１．．．ｆｎを得るように構成されたコンバータ３と、
前記周波数値ｆ１．．．ｆｎから量子化された周波数値ｆｑ１．．．ｆｑｎを得る量子化器４と、
前記量子化された周波数値ｆｑ１．．．ｆｑｎを含むビットストリームＢＳを生成するビットストリーム生成器５と、を備えている。
【００７３】
本発明による情報エンコーダ１は、零交差探索を用いるのに対して、従来技術による根を求めるためのスペクトル手法は、振幅スペクトルにおける谷の発見に依存する。しかしながら、谷を探索する場合の確度は、零交差を探索する場合より劣る。例えば、数列［４，２，１，２，３］を考察されたい。明らかに、最小値は、第３の要素であり、ゼロは、第２の要素と第４の要素との間のどこかに存在することになる。言い替えれば、ゼロが第３の要素の右側に存在するか、左側に存在するかを決定することはできない。しかしながら、数列［４，２，１，−２，−３］について考察すれば、直ちに、零交差が第３の要素と第４の要素との間に存在することが分かり、許容誤差が半減する。振幅スペクトル手法の場合、零交差探索の場合と同じ確度を達成するには、分析点の数を２倍にする必要があることになる。
【００７４】
振幅│Ｐ（ｚ）│と│Ｑ（ｚ）│の評価を比較すると、零交差手法は、確度において著しい優位点を有する。例えば、数列３，２，−１，−２を考察されたい。零交差手法では、ゼロは、２と−１との間に存在することが明らかである。しかしながら、対応する振幅列３，２，１，２を考察する場合、ゼロは、第２の要素と最後の要素との間のどこかに存在する、としか結論することができない。言い替えれば、零交差手法の場合の確度は、振幅ベースの手法の２倍である。
【００７５】
さらに、本発明による情報エンコーダは、ｍ＝１２８等の長い予測子を用いてもよい。これに対して、チェビシェフ変換は、Ａ（ｚ）の長さが比較的短い、例えばｍ≦２０である場合にしか十分に機能しない。長い予測子の場合、チェビシェフ変換は、数値的に不安定であり、アルゴリズムの実用的実装は不可能である。
【００７６】
したがって、提案する情報エンコーダ１の主要な特性は、零交差が探索されることに起因するとともに、時間領域から周波数領域への変換が行われ、極めて低い計算複雑性によりゼロを発見し得ることに起因して、チェビシェフベースの方法と同じく高い確度、又はより優れた確度を達成し得ることにある。
【００７７】
結果的に、本発明による情報エンコーダ１は、ゼロ（根）をより正確に決定するだけでなく、低い計算複雑性で決定する。
【００７８】
本発明による情報エンコーダ１は、数列の線スペクトルを決定する必要があるあらゆる信号処理アプリケーションにおいて使用可能である。本明細書では、例示的に、情報エンコーダ１を音声符号化の側面で論じる。本発明は、音声信号又は汎用オーディオ信号等の入力信号から線スペクトルを決定するための方法を必要とし、かつ入力信号がデジタルフィルタ又は他の数列として表現される、スペクトル振幅包絡線、知覚周波数マスキングしきい値、時間振幅包絡線、知覚時間マスキングしきい値、又は他の包絡線形状、又は線スペクトルを用いて包絡線情報を表現する、符号化、分析又は処理のための自己相関信号等の包絡線形状と等価の他の表現、をモデリングするための線形予測子を使用する、音声、オーディオ及び／又はビデオ符号化デバイス又はアプリケーションにおいて適用可能である。
【００７９】
情報信号ＩＳは、例えば、オーディオ信号又はビデオ信号であってもよい。
【００８０】
図２は、Ａ（ｚ）、Ｐ（ｚ）及びＱ（ｚ）の例示的な関係を示す。垂直な点線は、周波数値ｆ１．．．ｆ６を描いている。零交差が常に見えるように、振幅は、デシベルスケールではなく直線軸上に表されていることに留意されたい。線スペクトル周波数は、Ｐ（ｚ）とＱ（ｚ）の零交差において発生することが分かる。さらに、Ｐ（ｚ）及びＱ（ｚ）の振幅は、どこにおいても２│Ａ（ｚ）│以下、即ち、｜Ｐ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜及び｜Ｑ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜である。
【００８１】
図３は、本発明による情報エンコーダのコンバータの第１の実施形態を示す略図である。
【００８２】
本発明の好適な一実施形態によれば、コンバータ３は、予測多項式Ａ（ｚ）から多項式Ｐ（ｚ）とＱ（ｚ）を決定する決定デバイス６を備えている。
【００８３】
本発明の好適な一実施形態によれば、コンバータは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式を周波数領域へフーリエ変換するフーリエ変換デバイス８と、Ｐ（ｚ）から導出されるスペクトルＲＥＳの位相をそれが厳密に実数であるように調整するとともに、Ｑ（ｚ）から導出されるスペクトルＩＥＳの位相をそれが厳密に虚数であるように調整する調整デバイス７と、を備えている。フーリエ変換デバイス８は高速フーリエ変換を基礎とするものであってもよいし、離散型フーリエ変換を基礎とするものであってもよい。
【００８４】
本発明の好適な一実施形態によれば、調整デバイス７は、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式の係数を循環シフトする係数シフタ７として構成されている。
【００８５】
本発明の好適な一実施形態によれば、係数シフタ７は、係数列の原初の中間点がその数列の第１のポジションへシフトされるような方法で、係数を循環シフトするように構成されている。
【００８６】
理論的に、対称数列のフーリエ変換が実数値をとり、かつ逆対称の数列が純虚数フーリエスペクトルを有することが知られている。本事例における入力数列は、長さｍ＋ｌの多項式Ｐ（ｚ）又はＱ（ｚ）の係数であるが、これより遙かに長い長さＮ＞＞（ｍ＋ｌ）の離散型フーリエ変換を有する方が好まれると考えられる。より長いフーリエスペクトルを生成する従来的手法は、入力信号のゼロパディングである。しかし、数列のゼロパディングは、対称性が保たれるように慎重に実装されなければならない。
【００８７】
まず、係数、
［ｐ０，ｐ１，ｐ２，ｐ１，ｐ０］
を有する多項式Ｐ（ｚ）について考察する。
【００８８】
高速フーリエ変換アルゴリズムが通常適用される方法では、対称点が第１の要素であるということが要求され、それによって、例えばＭＡＴＬＡＢに適用する場合、実数値の出力を求めるために
ｆｆｔ（［ｐ２，ｐ１，ｐ０，ｐ０，ｐ１］）
のように書くことができる。具体的には、循環シフトが適用されてもよく、中間点要素に対応する対称点、すなわち係数ｐ２は、第１のポジションになるように左へシフトされる。次に、ｐ２の左側にあった係数は、数列の終わりに付加される。
【００８９】
ゼロパディングされた数列、
［ｐ０，ｐ１，ｐ２，ｐ１，ｐ０，０，０．．．０］
の場合も、同じプロセスを適用することができる。よって、数列、
［ｐ２，ｐ１，ｐ０，０，０．．．０，ｐ０，ｐ１］
が実数値の離散型フーリエ変換を有する。ここで、入力数列におけるゼロの数は、スペクトルの所望される長さをＮとすれば、Ｎ−ｍ−ｌである。
【００９０】
同様に、多項式Ｑ（ｚ）に対応する係数、
［ｑ０，ｑ１，０，−ｑ１，−ｑ０］
について考察されたい。前述の中間点が第１のポジションにくるように循環シフトを適用すれば、
［０，−ｑ１，−ｑ０，ｑ０，ｑ１］
が得られ、これは、純虚数の離散型フーリエ変換を有する。次には、この数列にゼロパディング変換を行なうことができる。
［０，−ｑ１，−ｑ０，０，０．．．０，ｑ０，ｑ１］
【００９１】
上記は、数列の長さが奇数である事例にのみ当てはまり、よって、ｍ＋ｌは偶数であることに留意されたい。ｍ＋ｌが奇数である事例に関しては、２つのオプションがある。すなわち、周波数領域において循環シフトを実装することができ、あるいは、ハーフサンプルでＤＦＴを適用することができる。
【００９２】
本発明の好適な実施形態によれば、コンバータ３は、Ｐ（ｚ）から導出される厳密な実数スペクトルＲＥＳとＱ（ｚ）から導出される厳密な虚数スペクトルＩＥＳのゼロを識別するためのゼロ識別子９を備えている。
【００９３】
本発明の好適な一実施形態によれば、ゼロ識別子９は、
ａ）実数スペクトルＲＥＳをヌル周波数において開始することと、
ｂ）実数スペクトルＲＥＳにおいて符号変化が発見されるまで周波数を増大することと、
ｃ）虚数スペクトルＩＥＳにおいてさらなる符号変化が発見されるまで周波数を増大することと、
ｄ）全てのゼロが発見されるまでステップｂ）及びステップｃ）を反復すること、によってゼロを識別するように構成されている。
【００９４】
Ｑ（ｚ）、ひいてはスペクトルの虚数部ＩＥＳは、ヌル周波数において常にゼロを有することに留意されたい。根は重なり合うことから、Ｐ（ｚ）、ひいてはスペクトルの実数部ＲＥＳは、ヌル周波数において常に非ゼロになる。したがって、ヌル周波数における実数部から開始して、最初の零交差、ひいては最初の周波数値ｆ１を示す最初の符号変化が発見されるまで、周波数を増大させることができる。
【００９５】
根はインタレースされることから、Ｑ（ｚ）のスペクトルＩＥＳは、次の符号変化を有する。したがって、Ｑ（ｚ）のスペクトルＩＥＳの符号変化が発見されるまで、周波数を増大させることができる。次に、このプロセスは、全ての周波数値ｆ１．．．ｆｎが発見されるまで、Ｐ（ｚ）とＱ（ｚ）のスペクトル間で交互に反復されてもよい。したがって、スペクトルＲＥＳとＩＥＳにおける零交差の位置決めに使用される手法は、チェビシェフ領域において適用される手法に類似する［６、７］。
【００９６】
Ｐ（ｚ）とＱ（ｚ）のゼロはインタレースされることから、１パスで全てのゼロを発見するように実数部ＲＥＳと複素数部ＩＥＳ上のゼロの探索を交互に行い、複雑性を完全探索の半分に減らすことができる。
【００９７】
本発明の好適な一実施形態によれば、ゼロ識別子９は補間によってゼロを識別するように構成されている。
【００９８】
零交差手法に加えて、ゼロのポジションを、例えば［７］のような従来方法において行われるようなより高い確度で推定できるように、補間を容易に適用することができる。
【００９９】
図４は、本発明による情報エンコーダ１のコンバータ３の第２の実施形態を示す略図である。
【０１００】
本発明の好適な一実施形態によれば、コンバータ３は、長くされた多項式対Ｐｅ（ｚ）とＱｅ（ｚ）を生成するために多項式Ｐ（ｚ）とＱ（ｚ）に対して値「０」を有する１つ又は複数の係数を加えるゼロパディング・デバイス１０を備えている。確度は、評価されるスペクトルＲＥＳ、ＩＥＳの長さを拡大することによってさらに向上させることができる。システムに関する情報を基礎として、事例によっては、実質上、周波数値ｆ１．．．ｆｎ間の最小距離を決定し、ひいては、全ての周波数値ｆ１．．．ｆｎをそれで発見することができるスペクトルＲＥＳ、ＩＥＳの最小長さを決定することが可能である［８］。
【０１０１】
本発明の好適な一実施形態によれば、コンバータ３は、線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現ＲＥＳ、ＩＥＳの周波数値ｆ１．．．ｆｎへ変換する間に、長くされた多項式Ｐｅ（ｚ）とＱｅ（ｚ）の値「０」を有することが知られている係数を用いる演算の少なくとも一部が省略されるように構成される。
【０１０２】
しかしながら、スペクトルの長さが増せば、計算複雑性も増大する。複雑性の最大の誘因は、Ａ（ｚ）の係数の、高速フーリエ変換等の時間領域から周波数領域への変換である。しかしながら、これは、係数ベクトルが所望される長さにゼロパディングされていることから、極めて疎である。複雑性を低減するためにこの事実を容易に使用することができる。これは、どの係数がゼロであるかが正確に分かり、高速フーリエ変換の各反復において単純にゼロを包含する演算を省くことができるという意味において、どちらかといえば単純な問題である。このように疎な高速フーリエ変換の適用は単純であり、当業者であるプログラマは誰でもこれを実装することができる。このような実装の複雑性は、Ｏ（Ｎｌｏｇ２（１＋ｍ＋ｌ））である。ここで、Ｎはスペクトルの長さ、ｍ及びｌは先に定義した通りである。
【０１０３】
本発明の好適な一実施形態によれば、コンバータは、長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）又は長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）から導出される１つ又は複数の多項式のスペクトルの数値範囲を、長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイス１１を備えている。ここで、フィルタ多項式Ｂ（ｚ）は対称性であって単位円上に根を持たない。Ｂ（ｚ）は、先に説明した通りに求めることができる。
【０１０４】
図５は、予測子Ａ（ｚ）、対応する平坦化フィルタＢ１（ｚ）とＢ２（ｚ）及び積Ａ（ｚ）Ｂ１（ｚ）とＡ（ｚ）Ｂ２（ｚ）の例示的な振幅スペクトルを示す。水平の点線は、０周波数とナイキスト周波数におけるＡ（ｚ）Ｂ１（ｚ）のレベルを示す。
【０１０５】
本発明の好適な一実施形態（不図示）によれば、コンバータ３は、多項式Ｐ（ｚ）及びＱ（ｚ）のスペクトルＲＥＳ、ＩＥＳの数値範囲を、多項式Ｐ（ｚ）及びＱ（ｚ）又は多項式Ｐ（ｚ）及びＱ（ｚ）から導出される１つ又は複数の多項式にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイス１１を備えている。ここで、フィルタ多項式Ｂ（ｚ）は対称性であって、単位円上に根を持たない。
【０１０６】
音声コーデックは、限定リソースを有するモバイルデバイス上で実装されることが多く、数値演算は固定小数点表示によって実装されなければならない。したがって、実装されるアルゴリズムは、範囲が限定される数値表現を用いて動作することが極めて重要である。しかしながら、一般的な音声スペクトル包絡線の場合、フーリエスペクトルの数値範囲があまりに大きく、零交差位置の保持を保証するためには、ＦＦＴの３２ビット実装が必要である。
【０１０７】
一方で、１６ビットのＦＦＴはより低い複雑性で実装可能である場合が多く、この１６ビット範囲内に適合するようにスペクトル値の範囲を限定することが有益であると考えられる。方程式｜Ｐ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜及び｜Ｑ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜からは、Ｂ（ｚ）Ａ（ｚ）の数値範囲を限定することにより、Ｂ（ｚ）Ｐ（ｚ）及びＢ（ｚ）Ｑ（ｚ）の数値範囲も限定されることが分かる。Ｂ（ｚ）が単位円上にゼロを持たなければ、Ｂ（ｚ）Ｐ（ｚ）及びＢ（ｚ）Ｑ（ｚ）は、単位円上にＰ（ｚ）及びＱ（ｚ）と同じ零交差を有する。さらに、Ｂ（ｚ）は、ｚ−（ｍ＋ｌ＋ｎ）／２Ｐ（ｚ）Ｂ（ｚ）及びｚ−（ｍ＋ｌ＋ｎ）／２Ｑ（ｚ）Ｂ（ｚ）が各々対称性及び逆対称を保ちかつそのスペクトルが純実数及び純虚数であるように、対称でなければならない。したがって、ｚ（ｎ＋ｌ）／２Ａ（ｚ）のスペクトルを評価する代わりに、ｚ（ｎ＋ｌ＋ｎ）／２Ａ（ｚ）Ｂ（ｚ）を評価することができる。ここで、Ｂ（ｚ）は、単位円上に根を持たない次数ｎの対称多項式である。言い替えれば、先に述べたものと同じ手法を適用することができるが、まずは、Ａ（ｚ）をフィルタＢ（ｚ）で乗算し、次に修正した位相シフトｚ−（ｍ＋ｌ＋ｎ）／２を適用する。
【０１０８】
残りのタスクは、Ａ（ｚ）Ｂ（ｚ）の数値範囲が、Ｂ（ｚ）は必ず対称性でありかつ単位円上に根を持たないという制約によって限定されるように、フィルタＢ（ｚ）を設計することである。この要件を満たす最も単純なフィルタは、次数２の線形位相フィルタ、Ｂ１（ｚ）＝β０＋β１ｚ−１＋β２ｚ−２である。ここで、βｋ∈Ｒはパラメータであり、かつ│β２│＞２│β１│である。βｋを調整すればスペクトル傾斜を修正することができ、積Ａ（ｚ）Ｂ１（ｚ）の数値範囲を縮小することができる。計算上極めて効率的な手法は、βを０周波数及びナイキストにおける振幅が等しくなるように、すなわち│Ａ（１）Ｂ１（１）│＝│Ａ（−１）Ｂ１（−１）│であるように選択することであり、例えば、β０＝Ａ（１）−Ａ（−１）及びβ１＝２（Ａ（１）＋Ａ（−１））を選ぶことができる。
【０１０９】
この手法は、略平坦なスペクトルを提供する。
【０１１０】
Ａ（ｚ）が高域通過特性を有するのに対して、Ｂ１（ｚ）は低域通過であり、積Ａ（ｚ）Ｂ１（ｚ）は、予期されるように、０周波数及びナイキスト周波数において同じ振幅を有し、かつこれは多かれ少なかれ平坦であるということが図５から観察される。Ｂ１（ｚ）は１自由度しか持たないので、積が完全に平坦になると予期できないことは明らかである。さらに、Ｂ１（ｚ）Ａ（ｚ）の最高ピークと最低の谷との比は、Ａ（ｚ）のそれより遙かに小さいものであり得ることを観察されたい。これは、Ｂ１（ｚ）Ａ（ｚ）の数値範囲がＡ（ｚ）のそれより遙かに小さいという望ましい効果が達成されていることを意味する。
【０１１１】
僅かに複雑性の高い、第２の方法は、Ａ（０．５ｚ）のインパルス応答の自己相関ｒｋを計算するものである。ここでは、０．５を乗算することによってＡ（ｚ）のゼロが原点方向へ移動し、それによってスペクトル振幅が約半分に縮小される。自己相関ｒｋにレビンソン−ダービンを適用すれば、最小位相である次数ｎのフィルタＨ（ｚ）が得られる。次には、Ｂ２（ｚ）＝ｚ−ｎＨ（ｚ）Ｈ（ｚ−１）を定義して、略一定である│Ｂ２（ｚ）Ａ（ｚ）│を求めることができる。│Ｂ２（ｚ）Ａ（ｚ）│の範囲は、│Ｂ１（ｚ）Ａ（ｚ）│のそれより小さいことが留意されるであろう。Ｂ（ｚ）を設計するためのさらなる手法は、ＦＩＲ設計に関する古典的文献［１８］において容易に見出すことができる。
【０１１２】
図６は、本発明による情報エンコーダ１のコンバータ３の第３の実施形態を示す略図である。
【０１１３】
本発明の好適な一実施形態によれば、調整デバイス１２は、フーリエ変換デバイス８の出力の位相をシフトするための位相シフタ１２として構成されている。
【０１１４】
本発明の好適な一実施形態によれば、位相シフタ１２は、フーリエ変換デバイス８の出力の位相を、ｋ番目の周波数ビンにｅｘｐ（ｉ２πｋｈ／Ｎ）を乗算することによってシフトするように構成されている。ここで、Ｎはサンプルの長さであり、かつｈ＝（ｍ＋ｌ）／２である。
【０１１５】
時間領域における循環シフトが周波数領域における位相回転と同一であることはよく知られている。具体的には、時間領域におけるｈ＝（ｍ＋ｌ）／２ステップのシフトは、ｋ番目の周波数ビンとｅｘｐ（−ｉ２πｋｈ／Ｎ）との乗算に一致する。ただし、Ｎはスペクトルの長さである。したがって、循環シフトの代わりに周波数領域における乗算を適用すれば、全く同じ結果を得ることができる。ただし、この手法では、複雑性が僅かに増す。ｈ＝（ｍ＋ｌ）／２は、ｍ＋ｌが偶数であるときにのみ整数であることに留意されたい。ｍ＋ｌが奇数のとき、循環シフトは有理数のステップ分の遅延を必要とすることになり、これを直に実装ことは困難である。代わりに、上述の相回転によって、周波数領域における対応するシフトを適用することができる。
【０１１６】
図７は、本発明による情報エンコーダ１のコンバータ３の第４の実施形態を示す略図である。
【０１１７】
本発明の好適な一実施形態によれば、コンバータ３は、多項式Ｐ（ｚ）及びＱ（ｚ）から合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を確立するように構成された合成多項式フォーマ１３を備えている。
【０１１８】
本発明の好適な一実施形態によれば、コンバータ３は、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）からの厳密な虚数スペクトルが、単一のフーリエ変換により、例えば高速フーリエ変換（ＦＦＴ）により合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を変換することによって確立されるように構成されている。
【０１１９】
多項式Ｐ（ｚ）及びＱ（ｚ）は、各々、ｚ−（ｍ＋ｌ）／２における対称軸と対称であり逆対称である。よって、単位円ｚ＝ｅｘｐ（ｉθ）上で評価されるｚ−（ｍ＋ｌ）／２Ｐ（ｚ）及びｚ−（ｍ＋ｌ）／２Ｑ（ｚ）のスペクトルは、各々、実数値及び複素数値ということになる。ゼロは単位円上に存在することから、零交差を探索すればゼロを発見することができる。さらに、単位円上の評価は、単純に高速フーリエ変換によって実行することができる。
【０１２０】
ｚ−（ｍ＋ｌ）／２Ｐ（ｚ）及びｚ−（ｍ＋ｌ）／２Ｑ（ｚ）に対応するスペクトルは、各々実数と複素数であることから、これらを単一の高速フーリエ変換によって実行することができる。具体的には、和ｚ−（ｍ＋ｌ）／２（Ｐ（ｚ）＋Ｑ（ｚ））を求めれば、スペクトルの実数部と複素数部は、各々、ｚ−（ｍ＋ｌ）／２Ｐ（ｚ）及びｚ−（ｍ＋ｌ）／２Ｑ（ｚ）に対応する。さらに、ｚ−（ｍ＋ｌ）／２（Ｐ（ｚ）＋Ｑ（ｚ））＝２ｚ−（ｍ＋ｌ）／２Ａ（ｚ）であることから、Ｐ（ｚ）及びＱ（ｚ）を明示的に決定することなく、２ｚ−（ｍ＋ｌ）／２Ａ（ｚ）のＦＦＴを直に取り込んでｚ−（ｍ＋ｌ）／２Ｐ（ｚ）及びｚ−（ｍ＋ｌ）／２Ｑ（ｚ）に対応するスペクトルを求めることができる。知りたいものはゼロの位置のみであることから、スカラ２による乗算を省略して、代わりにＦＦＴによるｚ−（ｍ＋ｌ）／２Ａ（ｚ）を評価することができる。Ａ（ｚ）は、ｍ＋１個の非ゼロ係数しか持たないことから、ＦＦＴ枝刈り（FFT pruning）を用いて複雑性を低減できることを観察されたい［１１］。全ての根が発見されることを保証するためには、２つのゼロ毎にその間の少なくとも１つの周波数上でスペクトルが評価されるに足る長い長さＮのＦＦＴを用いなければならない。
【０１２１】
本発明の好適な一実施形態（不図示）によれば、コンバータ３は、長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）から合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を確立するように構成された合成多項式フォーマを備えている。
【０１２２】
本発明の好適な一実施形態（不図示）によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）からの厳密な虚数スペクトルが、合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を変換することによって単一のフーリエ変換により確立されるように構成されている。
【０１２３】
図８は、本発明による情報エンコーダ１のコンバータ３の第５の実施形態を示す略図である。
【０１２４】
本発明の好適な実施形態によれば、コンバータ３は、多項式対Ｐ（ｚ）及びＱ（ｚ）又は多項式対Ｐ（ｚ）及びＱ（ｚ）から導出される１つ又は複数の多項式を、Ｐ（ｚ）から導出されるスペクトルが厳密に実数であり、かつＱ（ｚ）から導出されるスペクトルが厳密に虚数であるように、ハーフサンプルを用いて周波数領域へフーリエ変換するフーリエ変換デバイス１４を備えている。
【０１２５】
ある代替例は、ハーフサンプルを有するＤＦＴを設けたものである。具体的には、従来のＤＦＴが、
2019049729.tif 000005
と定義されるのに対して、ハーフサンプルＤＦＴは、次式のように定義することができる。
2019049729.tif 000006
【０１２６】
この公式化に対しては、ＦＦＴのような高速実装を容易に考案することができる。
【０１２７】
この公式化の利点は、対称点が通常のｎ＝１ではなく、ｎ＝１／２に存在することにある。よって、このハーフサンプルＤＦＴを用いれば、数列、
［２，１，０，０，１，２］
で、実数値のフーリエスペクトルＲＥＳが得られる。
【０１２８】
よって、ｍ＋ｌが奇数である場合、係数ｐ０，ｐ１，ｐ２，ｐ２，ｐ１，ｐ０を有する多項式Ｐ（ｚ）に関して、ハーフサンプルＤＦＴ及びゼロパディングを用いれば、次のような入力数列の場合に実数値スペクトルＲＥＳを求めることができる。
［ｐ２，ｐ１，ｐ０，０，０．．．０，ｐ０，ｐ１，ｐ２］
【０１２９】
同様に、多項式Ｑ（ｚ）の場合、数列、
［−ｑ２，−ｑ１，−ｑ０，０，０．．．０，ｑ０，ｑ１，ｑ２］
にハーフサンプルＤＦＴを適用して、純虚数スペクトルＩＥＳを求めることができる。
【０１３０】
これらの方法により、ｍ及びｌの如何なる組合せに対しても、多項式Ｐ（ｚ）の実数値スペクトル及び任意のＱ（ｚ）の純虚数スペクトルを求めることができる。実際には、Ｐ（ｚ）及びＱ（ｚ）のスペクトルが各々純実数と純虚数であることから、これらを単一の複素スペクトルに記憶することができ、これは、Ｐ（ｚ）＋Ｑ（ｚ）＝２Ａ（ｚ）のスペクトルに一致する。係数２によるスケーリングでは、根の位置に変化がなく、よってこれを無視することができる。したがって、単一のＦＦＴを用いてＡ（ｚ）のスペクトルのみを評価すれば、Ｐ（ｚ）及びＱ（ｚ）のスペクトルを求めることができる。先に説明したように、Ａ（ｚ）の係数に循環シフトを適用するだけでよい。
【０１３１】
例えば、ｍ＝４であり、ｌ＝０である場合、Ａ（ｚ）の係数は、
［ａ０，ａ１，ａ２，ａ３，ａ４］
であり、これは、ゼロパディングによって、次のように任意の長さＮにすることができる。
［ａ０，ａ１，ａ２，ａ３，ａ４，０，０．．．０］
【０１３２】
次に、（ｍ＋ｌ）／２＝２ステップの循環シフトを適用すれば、次のようになる。
［ａ２，ａ３，ａ４，０，０．．．０，ａ０，ａ１］
【０１３３】
この数列のＤＦＴを採用すれば、スペクトルの実数部ＲＥＳと複素数部ＩＥＳにＰ（ｚ）とＱ（ｚ）のスペクトルが存在する。
【０１３４】
ｍ＋ｌが偶数の事例における全体的アルゴリズムは、次のように述べることができる。ａｋで示されるＡ（ｚ）の係数は、長さＮのバッファ内に存在するものとする。
【０１３５】
１．左へ（ｍ＋ｌ）／２ステップのａｋに循環シフトを適用する。
【０１３６】
２．数列ａｋの高速フーリエ変換を計算し、これをＡｋで示す。
【０１３７】
３．全ての周波数値が発見されるまで、ｋ＝０を起点とし、次の（ａ）、（ｂ）を交互に行う。
（ａ）符号（実数（Ａｋ））＝符号（実数（Ａｋ＋１））である間に、ｋ：＝ｋ＋１を増大する。零交差が発見されれば、ｋを周波数値のリストに記憶する。
（ｂ）符号（虚数（Ａｋ））＝符号（虚数（Ａｋ＋１））である間に、ｋ：＝ｋ＋１を増大する。零交差が発見されれば、ｋを周波数値のリストに記憶する。
【０１３８】
４．各周波数値について、ＡｋとＡｋ＋１との間を補間し、正確なポジションを決定する。
【０１３９】
ここで、関数符号（ｘ）、実数（ｘ）及び虚数（ｘ）は、各々、ｘの符号、ｘの実数部及びｘの虚数部を指す。
【０１４０】
ｍ＋ｌが奇数の事例では、循環シフトが僅かに左へ（ｍ＋ｌ−１）／２ステップ低減され、通常の高速フーリエ変換がハーフサンプル高速フーリエ変換によって置換される。
【０１４１】
あるいは、循環シフトと第１のフーリエ変換との組合せを、常に、高速フーリエ変換及び周波数領域における位相シフトと置換することもできる。
【０１４２】
より正確な根位置を求めるためには、上述の提案方法を用いて最初の推定を行い、次に、根軌跡を精緻化（refine）する第２のステップを適用することが可能である。精緻化（refinement）のためには、デュラン−ケルナー法、アバース−エールリッヒ法、ラゲールのガウス−ニュートン法又はその他等［１１−１７］のあらゆる古典的な多項式根探索方法を適用することができる。
【０１４３】
ある公式化において、提示の本方法は、次のようなステップより成る。
【０１４４】
（ａ）長さＮまでゼロパディングされた長さｍ＋ｌ＋１の数列、ただしｍ＋ｌは偶数、に対して、バッファ長さがＮとなって出力スペクトルの所望される長さに一致するように、左へ（ｍ＋ｌ）／２ステップの循環シフトを適用する、又は、
長さＮまでゼロパディングされた長さｍ＋ｌ＋１の数列、ただしｍ＋ｌは奇数、に対して、バッファ長さがＮとなって出力スペクトルの所望される長さに一致するように、左へ（ｍ＋ｌ−１）／２ステップの循環シフトを適用する。
【０１４５】
（ｂ）ｍ＋ｌが偶数であれば、数列に通常のＤＦＴを適用する。ｍ＋ｌが奇数であれば、方程式３又は等価の表現が記述しているように、数列にハーフサンプルＤＦＴを適用する。
【０１４６】
（ｃ）入力信号が対称又は逆対称であれば、周波数領域表現の零交差を探索し、その位置をリストに記憶する。
【０１４７】
入力信号が合成列Ｂ（ｚ）＝Ｐ（ｚ）＋Ｑ（ｚ）であれば、周波数領域表現の実数部と虚数部の双方において零交差を探索し、その位置をリストに記憶する。入力信号が、合成列Ｂ（ｚ）＝Ｐ（ｚ）＋Ｑ（ｚ）であり、かつＰ（ｚ）とＱ（ｚ）の根が交互するか、又はそれと同様の構造を有する場合には、周波数領域表現の実数部と虚数部を交互して零交差を探索し、その位置をリストに記憶する。
【０１４８】
別の公式化では、本提示方法は、次のようなステップより成る。
【０１４９】
（ａ）形式が先行点と同じである入力信号に対して、入力数列にＤＦＴを適用する。
【０１５０】
（ｂ）周波数領域値に相回転を適用する。これは、入力数列の左へ（ｍ＋ｌ）／２ステップによる循環シフトに等しい。
【０１５１】
（ｃ）先行点で行われたような零交差探索を適用する。
【０１５２】
以下は、記述した実施形態によるエンコーダ１及び方法に関連する言及である。
【０１５３】
幾つかの態様は、装置のコンテキストにおいて記述されているが、これらの態様が、対応する方法の説明をも表現し、ブロック又はデバイスが方法ステップ又は方法ステップの特徴に対応していることは明らかである。同様に、方法ステップのコンテキストにおいて記述されている態様は、対応する装置の対応するブロック又はアイテム又は特徴の説明をも表現している。
【０１５４】
所定の実装要件によっては、本発明の実施形態は、ハードウェアにおいてもソフトウェアにおいても実装することができる。その実装は、個々の方法が実行されるようにプログラム可能コンピュータシステムと共働する（又は、共働することができる）電子読取り可能制御信号を格納しているデジタル記憶媒体、例えばフロッピーディスク、ＤＶＤ、ＣＤ、ＲＯＭ、ＰＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ又はフラッシュメモリを用いて実行することができる。
【０１５５】
本発明による幾つかの実施形態は、本明細書に記述している方法のうちの１つが実行されるようにプログラム可能コンピュータシステムと共働することができる電子読取り可能制御信号を有するデータキャリアを備える。
【０１５６】
概して、本発明の実施形態は、プログラムコードを有するコンピュータ・プログラム・プロダクトとしての実装が可能であり、プログラムコードは、コンピュータ上でコンピュータ・プログラム・プロダクトが実行されると本発明方法のうちの１つを実行するように作動する。プログラムコードは、例えば、機械読取り可能キャリア上に記憶されてもよい。
【０１５７】
他の実施形態は、機械読取り可能キャリア、又は非一時的記憶媒体上に記憶される、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムを含む。
【０１５８】
したがって、言い替えれば、本発明方法の一実施形態は、コンピュータ上でコンピュータプログラムが実行されると本明細書に記述している方法のうちの１つを実行するためのプログラムコードを有するコンピュータプログラムである。
【０１５９】
したがって、本発明方法のさらなる実施形態は、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムを記憶して備えるデータキャリア（又はデジタル記憶媒体、若しくはコンピュータ読取り可能媒体）である。
【０１６０】
したがって、本発明方法のさらなる実施形態は、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムを表すデータストリーム又は信号シーケンスである。データストリーム又は信号シーケンスは、例えば、データ通信接続を介して、例えばインターネットを介して伝送されるように構成されてもよい。
【０１６１】
さらなる実施形態は、本明細書に記述している方法のうちの１つを実行するように構成され、又は実行するように調整される処理手段、例えばコンピュータ、又はプログラマブル論理デバイス、を含む。
【０１６２】
さらなる実施形態は、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムをインストールしているコンピュータを含む。
【０１６３】
実施形態によっては、プログラム可能な論理デバイス（例えば、フィールド・プログラマブル・ゲート・アレイ）は、本明細書に記述している方法の機能のうちの一部又は全てを実行するために使用されてもよい。実施形態によっては、フィールド・プログラマブル・ゲート・アレイは、本明細書に記述している方法のうちの１つを実行するために、マイクロプロセッサと共働してもよい。概して、本方法は、効果的には、あらゆるハードウェア装置によって実行される。
【０１６４】
本発明を、幾つかの実施形態に関連して記述してきたが、本発明の範囲に含まれる変更、置換及び均等物が存在する。本発明の方法及び組成を実装する方法には、多くの代替方法が存在することも留意されるべきである。したがって、添付した以下のクレームの意図は、このような変更、置換及び均等物を本発明の精神及び範囲に包含されるものとして解釈することにある。
【符号の説明】
【０１６５】
１  情報エンコーダ
２  アナライザ
３  コンバータ
４  量子化器
５  ビットストリーム生成器
６  決定デバイス
７  係数シフタ
８  フーリエ変換デバイス
９  ゼロ識別子
１０  ゼロパディング・デバイス
１１  制限デバイス
１２  位相シフタ
１３  合成多項式フォーマ
１４  ハーフサンプル・フーリエ変換デバイス
ＩＳ  情報信号
ＲＥＳ  実数スペクトル
ＩＥＳ  虚数スペクトル
ｆ１．．．ｆｎ  周波数値
ｆｑ１．．．ｆｑｎ  量子化された周波数値
ＢＳ  ビットストリーム
【図１】
2019049729.tif 000007
【図２】
2019049729.tif 000008
【図３】
2019049729.tif 000009
【図４】
2019049729.tif 000010
【図５】
2019049729.tif 000011
【図６】
2019049729.tif 000012
【図７】
2019049729.tif 000013
【図８】
2019049729.tif 000014
【手続補正書】
【提出日】平成30年11月12日(2018.11.12)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
情報信号（ＩＳ）を符号化するための情報エンコーダであって、前記情報エンコーダ（１）は、
予測多項式Ａ（ｚ）の線形予測係数を得るために前記情報信号（ＩＳ）を分析するアナライザ（２）と、
前記予測多項式Ａ（ｚ）の前記線形予測係数を前記予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値ｆ１．．．ｆｎへ変換するコンバータ（３）であって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）、
ｍは前記予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することにより前記周波数値ｆ１．．．ｆｎを決定するように構成され、Ｐ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）からの厳密な虚数スペクトル（ＩＥＳ）を確立することによって、及びＰ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）のゼロを識別することによって、前記周波数値（ｆ１．．．ｆｎ）を得るように構成されたコンバータ（３）と、
前記周波数値（ｆ１．．．ｆｎ）から量子化された周波数値（ｆｑ１．．．ｆｑｎ）を取得する量子化器（４）と、
前記量子化された周波数値（ｆｑ１．．．ｆｑｎ）を含むビットストリームを生成するビットストリーム生成器（５）と、を備えている情報エンコーダ。
【請求項２】
前記コンバータ（３）は、前記予測多項式Ａ（ｚ）から前記多項式Ｐ（ｚ）及びＱ（ｚ）を決定する決定デバイス（６）を備えている請求項１に記載の情報エンコーダ。
【請求項３】
前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）の前記ゼロを識別するためのゼロ識別子（９）を備えている請求項１又は請求項２に記載の情報エンコーダ。
【請求項４】
前記ゼロ識別子（９）は、
ａ）前記実数スペクトル（ＲＥＳ）からヌル(null)周波数で開始すること、
ｂ）前記実数スペクトル（ＲＥＳ）において符号変化が発見されるまで周波数を増大すること、
ｃ）前記虚数スペクトル（ＩＥＳ）においてさらなる符号変化が発見されるまで周波数を増大すること、及び
ｄ）全てのゼロが発見されるまでステップｂ）とステップｃ）を反復すること、によって前記ゼロを識別するように構成されている請求項３に記載の情報エンコーダ。
【請求項５】
前記ゼロ識別子は、前記ゼロを補間(interpolation)によって識別するように構成されている請求項３又は請求項４に記載の情報エンコーダ。
【請求項６】
前記コンバータ（３）は、長くされた多項式対Ｐｅ（ｚ）及びＱｅ（ｚ）を生成するように、前記多項式Ｐ（ｚ）及びＱ（ｚ）に対して値「０」を有する１又は複数の係数を加えるゼロパディング・デバイス（１０）を備えている請求項１から５のいずれか一項に記載の情報エンコーダ。
【請求項７】
前記コンバータ（３）は、前記線形予測係数を前記予測多項式Ａ（ｚ）の前記スペクトル周波数表現（ＲＥＳ、ＩＥＳ）の周波数値（ｆ１．．．ｆｎ）に変換する間に、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）の前記値「０」を有することが知られている係数を用いる演算の少なくとも一部が省略されるように構成されている請求項５又は請求項６に記載の情報エンコーダ。
【請求項８】
前記コンバータ（３）は、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）から合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を確立するように構成された合成多項式フォーマ（１３）を備えている請求項５から７のいずれか一項に記載の情報エンコーダ。
【請求項９】
前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）からの前記厳密な虚数スペクトル（ＩＥＳ）が、前記合成多項式Ｃｅ（Ｐｅ（ｚ），Ｑｅ（ｚ））を変換することによって、単一のフーリエ変換により、確立されるように構成されている請求項８に記載の情報エンコーダ。
【請求項１０】
前記コンバータ（３）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１又は複数の多項式を周波数領域へフーリエ変換するフーリエ変換デバイス（８）と、Ｐ（ｚ）から導出される前記スペクトル（ＲＥＳ）の位相をそれが厳密に実数であるように調整するとともに、Ｑ（ｚ）から導出される前記スペクトル（ＩＥＳ）の位相をそれが厳密に虚数であるように調整する調整デバイス（７、１２）と、を備えている請求項１から９のいずれか一項に記載の情報エンコーダ。
【請求項１１】
前記調整デバイス（７、１２）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される前記１若しくは複数の多項式の係数を循環シフトする係数シフタ（７）として構成されている請求項１０に記載の情報エンコーダ。
【請求項１２】
前記係数シフタ（７）は、係数列の原初の中間点が前記数列の第１のポジションへシフトされるような方法で係数を循環シフトするように構成されている請求項１１に記載の情報エンコーダ。
【請求項１３】
前記調整デバイス（７、１２）は、前記フーリエ変換デバイス（８）の出力の位相をシフトする位相シフタ（１２）として構成されている請求項１０に記載の情報エンコーダ。
【請求項１４】
前記位相シフタ（１２）は、ｋ番目の周波数ビンにｅｘｐ（ｉ２πｋｈ／Ｎ）を乗算することによって、前記フーリエ変換デバイス（８）の前記出力の前記位相をシフトするように構成されており、Ｎは、サンプルの長さであり、ｈ＝（ｍ＋ｌ）／２である請求項１３に記載の情報エンコーダ。
【請求項１５】
前記コンバータ（３）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１又は複数の多項式を、Ｐ（ｚ）から導出される前記スペクトル（ＲＥＳ）が厳密に実数であるように、かつＱ（ｚ）から導出される前記スペクトル（ＩＥＳ）が厳密に虚数であるように、ハーフサンプルを用いて周波数領域へフーリエ変換するフーリエ変換デバイス（１４）を備えている請求項１から９のいずれか一項に記載の情報エンコーダ。
【請求項１６】
前記コンバータ（３）は、前記多項式Ｐ（ｚ）とＱ（ｚ）から合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を確立するように構成された合成多項式フォーマ（１３）を備えている請求項１から１５のいずれか一項に記載の情報エンコーダ。
【請求項１７】
前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）からの前記厳密な虚数スペクトル（ＩＥＳ）が、前記合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を変換することによって単一のフーリエ変換により確立されるように構成されている請求項１６に記載の情報エンコーダ。
【請求項１８】
前記コンバータ（３）は、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）又は前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）から導出される１つ若しくは複数の多項式の前記スペクトル（ＲＥＳ、ＩＥＳ）の数値範囲を、前記長くされた多項式Ｐｅ（ｚ）及びＱｅ（ｚ）にフィルタ多項式Ｂ（ｚ）を乗算することによって制限するための制限デバイス（１１）を備え、前記フィルタ多項式Ｂ（ｚ）は対称性でありかつ単位円上に根を持たない請求項６から１７のいずれか一項に記載の情報エンコーダ。
【請求項１９】
情報信号（ＩＳ）を符号化するための情報エンコーダ（１）を作動するための方法であって、
予測多項式Ａ（ｚ）の線形予測係数を求めるために前記情報信号（ＩＳ）を分析するステップと、
前記予測多項式Ａ（ｚ）の前記線形予測係数を前記予測多項式Ａ（ｚ）のスペクトル周波数表現（ＲＥＳ、ＩＥＳ）の周波数値（ｆ１．．．ｆｎ）へ変換するステップであって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ−ｍ−ｌＡ（ｚ−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ−ｍ−ｌＡ（ｚ−１）
ｍは前記予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することによって、前記周波数値（ｆ１．．．ｆｎ）を決定し、Ｐ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）からの厳密な虚数スペクトル（ＩＥＳ）を確立しかつＰ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）のゼロを識別することによって、前記周波数値（ｆ１．．．ｆｎ）を得るステップと、
前記周波数値（ｆ１．．．ｆｎ）から量子化された周波数値（ｆｑ１．．．ｆｑｎ）を求めるステップと、
前記量子化された周波数値（ｆｑ１．．．ｆｑｎ）を含むビットストリーム（ＢＳ）を生成するステップと、を含む方法。
【請求項２０】
プロセッサ上で実行されることにより請求項１９に記載の方法を実行するコンピュータプログラム。