特開 2021-6922 情報符号化のコンセプト

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

(12)【公報種別】公開特許公報(A)

(11)【公開番号】特開2021-6922(P2021-6922A)

(43)【公開日】2021年1月21日

(54)【発明の名称】情報符号化のコンセプト

(51)【国際特許分類】

   G10L 19/07 20130101AFI20201218BHJP

【ＦＩ】

   G10L19/07

【審査請求】有

【請求項の数】21

【出願形態】ＯＬ

【全頁数】29

(21)【出願番号】特願2020-164496(P2020-164496)

(22)【出願日】2020年9月30日

(62)【分割の表示】特願2018-192262(P2018-192262)の分割

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

(31)【優先権主張番号】14158396.3

(32)【優先日】2014年3月7日

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

(31)【優先権主張番号】14178789.5

(32)【優先日】2014年7月28日

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

【公序良俗違反の表示】

（特許庁注：以下のものは登録商標）

１．ＭＡＴＬＡＢ

(71)【出願人】

【識別番号】500341779

【氏名又は名称】フラウンホーファー−ゲゼルシャフト・ツール・フェルデルング・デル・アンゲヴァンテン・フォルシュング・アインゲトラーゲネル・フェライン

(74)【代理人】

【識別番号】100205981

【弁理士】

【氏名又は名称】野口 大輔

(72)【発明者】

【氏名】トム・ヴェックストローム

(72)【発明者】

【氏名】クリスティアン・フィッシャー ペデルセン

(72)【発明者】

【氏名】ヨハネス・フィッシャー

(72)【発明者】

【氏名】マティアス・ヒュッテンベルガー

(72)【発明者】

【氏名】アルフォンソ・ピノ

(57)【要約】      （修正有）

【課題】線形予測係数から線スペクトル周波数への変換及びその符号化の改良された手法を提供する。
【解決手段】予測多項式Ａ（ｚ）の線形予測係数を求めるために情報信号を分析するアナライザ（２）と、線形予測係数をスペクトル周波数表現の周波数値へ変換するコンバータ（３）であって、多項式対Ｐ（ｚ）とＱ（ｚ）を分析することにより周波数値を決定するように構成され、Ｐ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）からの厳密な虚数スペクトルを確立しかつＰ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）から導出される厳密な虚数スペクトルのゼロを識別することによって、周波数値を求めるように構成されているコンバータ（３）と、周波数値から量子化された周波数値を求める量子化器（４）と、量子化された周波数値を含むビットストリームＢＳを生成するビットストリーム生成器と、を備える。
【選択図】図１

【特許請求の範囲】

【請求項1】

情報信号（ＩＳ）を符号化するための情報エンコーダであって、前記情報エンコーダ（１）は、
予測多項式Ａ（ｚ）の線形予測係数を得るために前記情報信号（ＩＳ）を分析するアナライザ（２）と、
前記予測多項式Ａ（ｚ）の前記線形予測係数を前記予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値ｆ_１．．．ｆ_ｎへ変換するコンバータ（３）であって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ^−ｍ−ｌＡ（ｚ^−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ^−ｍ−ｌＡ（ｚ^−１）、
ｍは前記予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することにより前記周波数値ｆ_１．．．ｆ_ｎを決定するように構成され、Ｐ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）からの厳密な虚数スペクトル（ＩＥＳ）を確立することによって、及びＰ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）のゼロを識別することによって、前記周波数値（ｆ_１．．．ｆ_ｎ）を得るように構成されたコンバータ（３）と、
前記周波数値（ｆ_１．．．ｆ_ｎ）から量子化された周波数値（ｆ_ｑ１．．．ｆ_ｑｎ）を取得する量子化器（４）と、
前記量子化された周波数値（ｆ_ｑ１．．．ｆ_ｑｎ）を含むビットストリームを生成するビットストリーム生成器（５）と、を備えている情報エンコーダ。

【請求項2】

前記コンバータ（３）は、前記予測多項式Ａ（ｚ）から前記多項式Ｐ（ｚ）及びＱ（ｚ）を決定する決定デバイス（６）を備えている請求項１に記載の情報エンコーダ。

【請求項3】

前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）の前記ゼロを識別するためのゼロ識別子（９）を備えている請求項１又は請求項２に記載の情報エンコーダ。

【請求項4】

前記ゼロ識別子（９）は、
ａ）前記実数スペクトル（ＲＥＳ）からヌル(null)周波数で開始すること、
ｂ）前記実数スペクトル（ＲＥＳ）において符号変化が発見されるまで周波数を増大すること、
ｃ）前記虚数スペクトル（ＩＥＳ）においてさらなる符号変化が発見されるまで周波数を増大すること、及び
ｄ）全てのゼロが発見されるまでステップｂ）とステップｃ）を反復すること、によって前記ゼロを識別するように構成されている請求項３に記載の情報エンコーダ。

【請求項5】

前記ゼロ識別子は、前記ゼロを補間(interpolation)によって識別するように構成されている請求項３又は請求項４に記載の情報エンコーダ。

【請求項6】

前記コンバータ（３）は、長くされた多項式対Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）を生成するように、前記多項式Ｐ（ｚ）及びＱ（ｚ）に対して値「０」を有する１又は複数の係数を加えるゼロパディング・デバイス（１０）を備えている請求項１から５のいずれか一項に記載の情報エンコーダ。

【請求項7】

前記コンバータ（３）は、前記線形予測係数を前記予測多項式Ａ（ｚ）の前記スペクトル周波数表現（ＲＥＳ、ＩＥＳ）の周波数値（ｆ_１．．．ｆ_ｎ）に変換する間に、前記長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）の前記値「０」を有することが知られている係数を用いる演算の少なくとも一部が省略されるように構成されている請求項５又は請求項６に記載の情報エンコーダ。

【請求項8】

前記コンバータ（３）は、前記長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）から合成多項式Ｃ_ｅ（Ｐ_ｅ（ｚ），Ｑ_ｅ（ｚ））を確立するように構成された合成多項式フォーマ（１３）を備えている請求項５から７のいずれか一項に記載の情報エンコーダ。

【請求項9】

前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）からの前記厳密な虚数スペクトル（ＩＥＳ）が、前記合成多項式Ｃ_ｅ（Ｐ_ｅ（ｚ），Ｑ_ｅ（ｚ））を変換することによって、単一のフーリエ変換により、確立されるように構成されている請求項８に記載の情報エンコーダ。

【請求項10】

前記コンバータ（３）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１又は複数の多項式を周波数領域へフーリエ変換するフーリエ変換デバイス（８）と、Ｐ（ｚ）から導出される前記スペクトル（ＲＥＳ）の位相をそれが厳密に実数であるように調整するとともに、Ｑ（ｚ）から導出される前記スペクトル（ＩＥＳ）の位相をそれが厳密に虚数であるように調整する調整デバイス（７、１２）と、を備えている請求項１から９のいずれか一項に記載の情報エンコーダ。

【請求項11】

前記調整デバイス（７、１２）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される前記１若しくは複数の多項式の係数を循環シフトする係数シフタ（７）として構成されている請求項１０に記載の情報エンコーダ。

【請求項12】

前記係数シフタ（７）は、係数列の原初の中間点が前記数列の第１のポジションへシフトされるような方法で係数を循環シフトするように構成されている請求項１１に記載の情報エンコーダ。

【請求項13】

前記調整デバイス（７、１２）は、前記フーリエ変換デバイス（８）の出力の位相をシフトする位相シフタ（１２）として構成されている請求項１０に記載の情報エンコーダ。

【請求項14】

前記位相シフタ（１２）は、ｋ番目の周波数ビンにｅｘｐ（ｉ２πｋｈ／Ｎ）を乗算することによって、前記フーリエ変換デバイス（８）の前記出力の前記位相をシフトするように構成されており、Ｎは、サンプルの長さであり、ｈ＝（ｍ＋ｌ）／２である請求項１３に記載の情報エンコーダ。

【請求項15】

前記コンバータ（３）は、前記多項式対Ｐ（ｚ）とＱ（ｚ）又は前記多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１又は複数の多項式を、Ｐ（ｚ）から導出される前記スペクトル（ＲＥＳ）が厳密に実数であるように、かつＱ（ｚ）から導出される前記スペクトル（ＩＥＳ）が厳密に虚数であるように、ハーフサンプルを用いて周波数領域へフーリエ変換するフーリエ変換デバイス（１４）を備えている請求項１から９のいずれか一項に記載の情報エンコーダ。

【請求項16】

前記コンバータ（３）は、前記多項式Ｐ（ｚ）とＱ（ｚ）から合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を確立するように構成された合成多項式フォーマ（１３）を備えている請求項１から１５のいずれか一項に記載の情報エンコーダ。

【請求項17】

前記コンバータ（３）は、Ｐ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）とＱ（ｚ）からの前記厳密な虚数スペクトル（ＩＥＳ）が、前記合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を変換することによって単一のフーリエ変換により確立されるように構成されている請求項１６に記載の情報エンコーダ。

【請求項18】

前記コンバータ（３）は、前記多項式Ｐ（ｚ）及びＱ（ｚ）の前記スペクトル（ＲＥＳ、ＩＥＳ）の数値範囲を、前記多項式Ｐ（ｚ）及びＱ（ｚ）又は前記多項式Ｐ（ｚ）及びＱ（ｚ）から導出される１つ若しくは複数の多項式にフィルタ多項式Ｂ（ｚ）を乗算することによって制限するための制限デバイス（１１）を備え、前記フィルタ多項式Ｂ（ｚ）は対称性でありかつ単位円上に根を持たない請求項１から１７のいずれか一項に記載の情報エンコーダ。

【請求項19】

前記コンバータ（３）は、前記長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）又は前記長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）から導出される１つ若しくは複数の多項式の前記スペクトル（ＲＥＳ、ＩＥＳ）の数値範囲を、前記長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）にフィルタ多項式Ｂ（ｚ）を乗算することによって制限するための制限デバイス（１１）を備え、前記フィルタ多項式Ｂ（ｚ）は対称性でありかつ単位円上に根を持たない請求項６から１８のいずれか一項に記載の情報エンコーダ。

【請求項20】

情報信号（ＩＳ）を符号化するための情報エンコーダ（１）を作動するための方法であって、
予測多項式Ａ（ｚ）の線形予測係数を求めるために前記情報信号（ＩＳ）を分析するステップと、
前記予測多項式Ａ（ｚ）の前記線形予測係数を前記予測多項式Ａ（ｚ）のスペクトル周波数表現（ＲＥＳ、ＩＥＳ）の周波数値（ｆ_１．．．ｆ_ｎ）へ変換するステップであって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ^−ｍ−ｌＡ（ｚ^−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ^−ｍ−ｌＡ（ｚ^−１）
ｍは前記予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することによって、前記周波数値（ｆ_１．．．ｆ_ｎ）を決定し、Ｐ（ｚ）から導出される厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）からの厳密な虚数スペクトル（ＩＥＳ）を確立しかつＰ（ｚ）から導出される前記厳密な実数スペクトル（ＲＥＳ）及びＱ（ｚ）から導出される前記厳密な虚数スペクトル（ＩＥＳ）のゼロを識別することによって、前記周波数値（ｆ_１．．．ｆ_ｎ）を得るステップと、
前記周波数値（ｆ_１．．．ｆ_ｎ）から量子化された周波数値（ｆ_ｑ１．．．ｆ_ｑｎ）を求めるステップと、
前記量子化された周波数値（ｆ_ｑ１．．．ｆ_ｑｎ）を含むビットストリーム（ＢＳ）を生成するステップと、を含む方法。

【請求項21】

プロセッサ上で実行されることにより請求項２０に記載の方法を実行するコンピュータプログラム。

【発明の詳細な説明】

【背景技術】

【0001】

音声符号化において最も頻繁に使用されるパラダイムは、代数符号励振線形予測（ＡＣＥＬＰ）であり、これはＡＭＲ−ファミリ、Ｇ．７１８やＭＰＥＧ ＵＳＡＣ［１−３］等のような規格に使用されている。これは、ソースモデルを用いる音声モデリングを基礎としている。そのようなソースモデルは、スペクトル包絡線をモデリングするための線形予測子（ＬＰ）、基本周波数をモデリングするための長時間予測子（ＬＴＰ）及び残余（residual）のための代数コードブックにより構成されている。

【0002】

線形予測モデルの係数は、量子化に対する感受性が極めて高く、よって通常これらは、量子化に先立って、まず線スペクトル周波数（ＬＳＦ）又はイミッタンススペクトル周波数（ＩＳＦ）に変換される。ＬＳＦ／ＩＳＦ領域は、量子化誤差に対して堅牢（robust）であり、これらの領域では、予測子の安定性を容易に保全することができ、よって、量子化のための適切な領域が提供される［４］。

【0003】

以下、周波数値と称するＬＳＦ／ＩＳＦは、次のように、次数ｍの線形予測多項式Ａ（ｚ）から求めることができる。線スペクトル対の多項式は、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ^−ｍ−ｌＡ（ｚ^−１）
Ｑ（ｚ）＝Ａ（ｚ）−ｚ^−ｍ−ｌＡ（ｚ^−１） （１）
と定義される。ここで、線スペクトル対表現ではｌ＝１であり、かつイミッタンススペクトル対表現ではｌ＝０であるが、原則的には、任意のｌ≧０が有効である。したがって、以後、単にｌ≧０であるものと仮定する。

【0004】

元の予測子は、Ａ（ｚ）＝１／２［Ｐ（ｚ）＋Ｑ（ｚ）］を用いていつでも復元できることに留意されたい。したがって、多項式Ｐ（ｚ）及びＱ（ｚ）は、Ａ（ｚ）の全情報を包含する。

【0005】

ＬＳＰ／ＩＳＰ多項式の中心的性質は、Ａ（ｚ）がその全ての根を単位円内に有する場合に限り、Ｐ（ｚ）及びＱ（ｚ）の根が単位円上でインタレースされる、というものである。Ｐ（ｚ）及びＱ（ｚ）の根が単位円上に存在するため、それらを角度のみで表現することができる。これらの角度は、周波数に対応し、またＰ（ｚ）及びＱ（ｚ）のスペクトルは、根に対応する周波数においてその対数振幅スペクトルにおける垂線を有することから、これらの根を周波数値と称する。

【0006】

周波数値は、結果的に、予測子Ａ（ｚ）の全情報を符号化することになる。さらに、周波数値のうちの１つにおける僅かな誤差がもたらす、対応する周波数の近くにスペクトルが位置決めされる復元された予測子のスペクトル誤差が小さくなるように、周波数値は、量子化誤差に対して堅牢（robust）であることが分かっている。これらの好ましい性質に起因して、ＬＳＦ又はＩＳＦ領域における量子化は、全てのメインストリーム音声コーデックにおいて使用されている［１−３］。

【0007】

しかしながら、周波数値を用いる場合の課題の１つは、多項式Ｐ（ｚ）及びＱ（ｚ）の係数からその位置を効率的に求めることにある。結局のところ、多項式の根を求めることは、古典的かつ困難な課題である。このタスクに関してこれまでに提案されている方法は、下記の手法を含んでいる：
・ 初期の手法のうちの１つは、ゼロが単位円上に存在し、よって、これが振幅スペクトルにおいてゼロとして現出する、ということを用いている［５］。したがって、Ｐ（ｚ）とＱ（ｚ）の係数の離散型フーリエ変換を利用することにより、振幅スペクトルの谷を探索することができる。各谷は、根の位置を示し、スペクトルが充分にアップサンプリングされていれば、全ての根を求めることができる。しかしながら、この方法では、谷の位置からの正確なポジションの決定が困難であるため、近似ポジションしか得られない。
・ 最も高頻度で使用される手法は、チェビシェフ多項式を基礎とするものであり、［６］に提示されている。これは、多項式Ｐ（ｚ）及びＱ（ｚ）はそれぞれが対称かつ逆対称であり、それによって多くの冗長な情報を含む、という認識に依存する。ｚ＝±１における自明な零点を除去することにより、かつ置換ｘ＝ｚ＋ｚ^−１（チェビシェフ変換として知られる）を用いれば、多項式を代替表現ＦＰ（ｘ）及びＦＱ（ｘ）へ変換することができる。これらの多項式の次数は、Ｐ（ｚ）とＱ（ｚ）の半分であり、よってこれらは、−２から＋２までの範囲にのみ実根を有する。多項式ＦＰ（ｘ）及びＦＱ（ｘ）は、ｘが実数のとき実数値であることに留意されたい。さらに、これらの根は、単根であることから、ＦＰ（ｘ）とＦＱ（ｘ）は、その各根において零交差を有する。

【0008】

ＡＭＲ−ＷＢ等の音声コーデックにおいて、この手法は、実軸上の固定格子上で多項式ＦＰ（ｘ）及びＦＱ（ｘ）を評価して全ての零交差を求めるように適用される。根位置は、零交差の周囲での線形補間によってさらに精緻化される。この手法の優位点は、冗長係数の省略により複雑さが低減されることにある。

【0009】

上述の方法は、既存のコーデックにおいて充分に機能するが、幾つかの問題点もある。

【先行技術文献】

【非特許文献】

【0010】

【非特許文献1】[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】[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】[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】[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】[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】[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】[7] 3GPP TS 26.190 V7.0.0, “Adaptive multi-rate (AMR-WB) speech codec”, 2007.

【非特許文献8】[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】[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】[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】[11] E. Durand, Solutions Numeriques des Equations Algebriques. Paris: Masson, 1960.

【非特許文献12】[12] I. Kerner, “Ein Gesamtschrittverfahren zur Berechnung der Nullstellen von Polynomen”, Numerische Mathematik, vol. 8, no. 3, pp. 290-294, May 1966.

【非特許文献13】[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】[14] L. Ehrlich, “A modified newton method for polynomials”, Communications of the ACM, vol. 10, no. 2, pp. 107-108, February 1967.

【非特許文献15】[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】[16] --, “Adaptive polynomial factorization by coefficient matching”, IEEE Transactions on Signal Processing, vol. 39, no. 2, pp. 527-530, February 1991.

【非特許文献17】[17] G. H. Golub and C. F. van Loan, Matrix Computations, 3rd ed. John Hopkins University Press, 1996.

【非特許文献18】[18] T. Saramaeki, “Finite impulse response filter design”, Handbook for Digital Signal Processing, pp. 155-277, 1993.

【発明の概要】

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

【0011】

解決すべき課題は、情報を符号化するための改良された概念を提供することである。

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

【0012】

第１の態様において、この課題は、情報信号を符号化するための情報エンコーダによって解決される。情報エンコーダは、
予測多項式Ａ（ｚ）の線形予測係数を得るために情報信号を分析するアナライザと、
予測多項式Ａ（ｚ）の線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値へ変換するコンバータであって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ^−ｍ−ｌＡ（ｚ^−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ^−ｍ−ｌＡ（ｚ^−１）、
ｍは予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
として定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することにより前記周波数値を決定するように構成され、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）からの厳密な虚数スペクトルを確立することによって、及びＰ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）から導出される厳密な虚数スペクトルのゼロを識別することによって、前記周波数値を得るように構成されたコンバータと、
前記周波数値から量子化された周波数値を取得する量子化器と、
前記量子化された周波数値を含むビットストリームを生成するビットストリーム生成器と、を備える。

【0013】

本発明による情報エンコーダは零交差探索を用いるのに対し、従来技術による根を求めるためのスペクトル手法は、振幅スペクトル中の谷の検知に依存する。しかしながら、谷を探索する場合の確度は、零交差を探索する場合より劣る。例えば、数列［４，２，１，２，３］を考察されたい。明らかに、最小値は第３の要素であり、よって、ゼロは、第２の要素と第４の要素との間のどこかに存在することになる。言い替えれば、ゼロが第３の要素の右側に存在するか、左側に存在するかを決定することはできない。しかしながら、数列［４，２，１，−２，−３］について考察する場合は、直ちに、零交差が第３の要素と第４の要素との間に存在することが分かり、よって、許容誤差が半減する。振幅スペクトル手法の場合、零交差探索の場合と同じ確度を達成するには、分析点の数を２倍にする必要があることになる。

【0014】

振幅│Ｐ（ｚ）│と│Ｑ（ｚ）│の評価を比較すれば、零交差手法は、確度において著しい優位点を有する。例えば、数列３，２，−１，−２を考察されたい。零交差手法では、ゼロは、２と−１との間に存在することが明らかである。しかしながら、対応する振幅列３，２，１，２を考察する場合、ゼロは、第２の要素と最後の要素との間のどこかに存在する、としか結論することができない。言い替えれば、零交差手法の場合の確度は、振幅ベースの手法の２倍である。

【0015】

さらに、本発明による情報エンコーダは、ｍ＝１２８等の長い予測子を用いてもよい。これに対して、チェビシェフ変換は、Ａ（ｚ）の長さが比較的短い、例えばｍ≦２０である場合にしか十分に機能しない。長い予測子の場合、チェビシェフ変換は、数値的に不安定であり、よって、アルゴリズムの実用的実装は不可能である。

【0016】

したがって、提案する情報エンコーダの主要な特性は、零交差が探索されることに起因して、そして時間領域から周波数領域への変換が行われ、それによって極めて低い計算複雑性によりゼロを発見し得ることに起因して、チェビシェフベースの方法のように高い確度、又はより優れた確度を達成し得ることである。

【0017】

結果的に、本発明による情報エンコーダは、ゼロ（根）をより正確に決定するだけでなく、低い計算複雑性で決定する。

【0018】

本発明による情報エンコーダは、数列の線スペクトルを決定する必要があるあらゆる信号処理アプリケーションにおいて使用可能である。本明細書では、例示的に、情報エンコーダを音声符号化の側面で論じる。本発明は、音声信号又は汎用オーディオ信号等の入力信号から線スペクトルを決定するための方法を必要とし、かつ入力信号がデジタルフィルタ又は他の数列として表現される、スペクトル振幅包絡線、知覚周波数マスキングしきい値、時間振幅包絡線、知覚時間マスキングしきい値、又は他の包絡線形状、若しくは線スペクトルを用いて包絡線情報を表現する、符号化、分析又は処理のための自己相関信号等の包絡線形状と等価の他の表現、をモデリングするための線形予測子を使用する、音声、オーディオ及び／又はビデオ符号化デバイス又はアプリケーションにおいて適用可能である。

【0019】

情報信号は、例えば、オーディオ信号又はビデオ信号であってもよい。周波数値は、線スペクトル周波数であっても、イミッタンススペクトル周波数であってもよい。ビットストリーム内で伝送される量子化された周波数値は、オーディオ信号又はビデオ信号を再現するために、デコーダがビットストリームを復号できるようにする。

【0020】

本発明の好適な一実施形態によれば、コンバータは、予測多項式Ａ（ｚ）から多項式Ｐ（ｚ）とＱ（ｚ）を決定するための決定デバイスを備えている。

【0021】

本発明の好適な実施形態によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）から導出される厳密な虚数スペクトルのゼロを識別するためのゼロ識別子を備えている。

【0022】

本発明の好適な一実施形態によれば、ゼロ識別子は、
ａ）実数スペクトルをヌル周波数において開始することと、
ｂ）実数スペクトルにおいて符号変化が発見されるまで周波数を増大することと、
ｃ）虚数スペクトルにおいてさらなる符号変化が発見されるまで周波数を増大することと、
ｄ）全てのゼロが発見されるまでステップｂ）及びステップｃ）を反復すること、によってゼロを識別するように構成されている。

【0023】

Ｑ（ｚ）ひいてはスペクトルの虚数部は、ヌル周波数において常にゼロを有することに留意されたい。根は重なり合うことから、Ｐ（ｚ）ひいてはスペクトルの実数部は、ヌル周波数において常に非ゼロになる。したがって、ヌル周波数における実数部から開始して、最初の零交差ひいては最初の周波数値を示す最初の符号変化が発見されるまで、周波数を増大させることができる。

【0024】

根はインタレースされるので、Ｑ（ｚ）のスペクトルは次の符号変化を有することとなる。したがって、Ｑ（ｚ）のスペクトルの符号変化が発見されるまで、周波数を増大させることができる。次に、このプロセスは、全ての周波数値が発見されるまで、Ｐ（ｚ）及びＱ（ｚ）のスペクトル間で交互に反復されてもよい。したがって、スペクトルにおける零交差の位置決めに使用される手法は、チェビシェフ領域において適用される手法に類似する［６、７］。

【0025】

Ｐ（ｚ）及びＱ（ｚ）のゼロはインタレースされることから、１パスで全てのゼロを発見するように実数部及び複素数部上のゼロの探索を交互に行い、よって、複雑性を完全探索の半分に減らすことができる。

【0026】

本発明の好適な一実施形態によれば、ゼロ識別子は、ゼロを補間によって識別するように構成される。

【0027】

零交差手法に加えて、ゼロのポジションを、例えば［７］のような従来方法において行われるようなより高い確度で推定できるように、補間を容易に適用することができる。

【0028】

本発明の好適な一実施形態によれば、コンバータは、長くされた多項式対Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）を生成するための多項式Ｐ（ｚ）及びＱ（ｚ）へ値「０」を有する１つ又は複数の係数を加えるゼロパディング・デバイスを備えている。確度は、評価されるスペクトルの長さを拡大することによってさらに向上させることができる。システムに関する情報を基礎として、事例によっては、実質上、周波数値間の最小距離を決定し、ひいては、全ての周波数値をそれにより発見することができるスペクトルの最小長さを決定することが可能である［８］。

【0029】

本発明の好適な一実施形態によれば、コンバータは、線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値へ変換する間に、長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）の値「０」を有することが知られている係数を用いる演算の少なくとも一部が省略されるように構成される。

【0030】

しかしながら、スペクトルの長さが増せば、計算複雑性も増大する。複雑性の最大の誘因は、Ａ（ｚ）の係数の、時間領域から周波数領域への高速フーリエ変換などの変換である。しかしながら、これは、係数ベクトルが所望される長さにゼロパディングされていることから、極めて疎である。複雑性の低減には、この事実を容易に使用することができる。これは、どの係数がゼロであるかが正確に分かり、よって、高速フーリエ変換の各反復において単純にゼロを包含する演算を省くことができるという意味において、どちらかといえば単純な問題である。このように疎な高速フーリエ変換の適用は、単純であり、よって当業者たるプログラマは誰でもこれを実装することができる。このような実装の複雑性は、Ｏ（Ｎｌｏｇ_２（１＋ｍ＋ｌ））であり、ここで、Ｎは、スペクトルの長さであり、ｍ及びｌは、先に定義した通りである。

【0031】

本発明の好適な一実施形態によれば、コンバータは、長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）から合成多項式Ｃ_ｅ（Ｐ_ｅ（ｚ），Ｑ_ｅ（ｚ））を確立するように構成された合成多項式フォーマを備えている。

【0032】

本発明の好適な一実施形態によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）からの厳密な虚数スペクトルが、合成多項式Ｃ_ｅ（Ｐ_ｅ（ｚ），Ｑ_ｅ（ｚ））を変換することによって単一のフーリエ変換により、確立されるように構成される。

【0033】

本発明の好適な一実施形態によれば、コンバータは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式を周波数領域へフーリエ変換するフーリエ変換デバイスと、Ｐ（ｚ）から導出されるスペクトルの位相を、それが厳密に実数であるように調整するとともに、Ｑ（ｚ）から導出されるスペクトルの位相を、それが厳密に虚数であるように調整する調整デバイスと、を備えている。フーリエ変換デバイスは、高速フーリエ変換を基礎とするものであっても、離散型フーリエ変換を基礎とするものであってもよい。

【0034】

本発明の好適な一実施形態によれば、調整デバイスは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式の係数を循環シフトする係数シフタとして構成されている。

【0035】

本発明の好適な一実施形態によれば、係数シフタは、係数列の原初の中間点がその数列の第１のポジションへシフトされるような方法で、係数を循環シフトするように構成されている。

【0036】

理論的には、対称数列のフーリエ変換が実数値をとり、かつ逆対称の数列が純虚数フーリエスペクトルを有することは周知である。本事例における入力数列は、長さｍ＋ｌの多項式Ｐ（ｚ）又はＱ（ｚ）の係数であるが、これより遙かに長い長さＮ＞＞（ｍ＋ｌ）の離散型フーリエ変換を有する方が好まれると思われる。より長いフーリエスペクトルを生成する従来的手法は、入力信号のゼロパディングである。しかしながら、数列のゼロパディングは、対称性が保たれるように慎重に実装されなければならない。

【0037】

まず、係数、
［ｐ_０，ｐ_１，ｐ_２，ｐ_１，ｐ_０］
を有する多項式Ｐ（ｚ）について考察する。

【0038】

通常、ＦＦＴアルゴリズムを適用する方法は、対称点が第１の要素であることを要求し、よって、例えば、ＭＡＴＬＡＢに適用する場合、
ｆｆｔ（［ｐ_２，ｐ_１，ｐ_０，ｐ_０，ｐ_１］）
のように書いて実数値の出力を求めることができる。具体的には、循環シフトが適用されてもよく、よって、中間点要素に対応する対称点、即ち係数ｐ_２は、第１のポジションになるように左へシフトされる。次に、ｐ_２の左側にあった係数は、数列の終わりに付加される。

【0039】

ゼロパディングされた数列、
［ｐ_０，ｐ_１，ｐ_２，ｐ_１，ｐ_０，０，０．．．０］
の場合も、同じプロセスを適用することができる。よって、数列、
［ｐ_２，ｐ_１，ｐ_０，０，０．．．０，ｐ_０，ｐ_１］
が実数値の離散型フーリエ変換を有する。ここで、入力数列におけるゼロの数は、スペクトルの所望される長さをＮとすれば、Ｎ−ｍ−ｌである。

【0040】

同様に、多項式Ｑ（ｚ）に対応する係数、
［ｑ_０，ｑ_１，０，−ｑ_１，−ｑ_０］
について考察されたい。前述の中間点が第１のポジションにくるように循環シフトを適用すれば、
［０，−ｑ_１，−ｑ_０，ｑ_０，ｑ_１］
が得られ、これは、純虚数の離散型フーリエ変換を有する。次には、この数列にゼロパディング変換を行なうことができる。
［０，−ｑ_１，−ｑ_０，０，０．．．０，ｑ_０，ｑ_１］

【0041】

上記は、数列の長さが奇数である事例にのみ当てはまり、よって、ｍ＋ｌは偶数であることに留意されたい。ｍ＋ｌが奇数である事例に関しては、２つのオプションがある。即ち、周波数領域において循環シフトを実装することができ、あるいは、ハーフサンプルでＤＦＴを適用することができる（下記参照）。

【0042】

本発明の好適な一実施形態によれば、調整デバイスは、フーリエ変換デバイスの出力の位相をシフトするための位相シフタとして構成される。

【0043】

本発明の好適な一実施形態によれば、位相シフタは、フーリエ変換デバイスの出力の位相を、ｋ番目の周波数ビンにｅｘｐ（ｉ２πｋｈ／Ｎ）を乗算することによってシフトするように構成される。ここで、Ｎは、サンプルの長さであり、かつｈ＝（ｍ＋ｌ）／２である。

【0044】

時間領域における循環シフトが周波数領域における位相回転と同一であることはよく知られている。具体的には、時間領域におけるｈ＝（ｍ＋ｌ）／２ステップのシフトは、ｋ番目の周波数ビンとｅｘｐ（−ｉ２πｋｈ／Ｎ）との乗算に一致する。ここで、Ｎはスペクトルの長さである。したがって、循環シフトの代わりに、周波数領域における乗算を適用すれば、全く同じ結果を得ることができる。ただし、この手法では、複雑性が僅かに増す。ｈ＝（ｍ＋ｌ）／２は、ｍ＋ｌが偶数であるときにのみ整数であることに留意されたい。ｍ＋ｌが奇数のとき、循環シフトは、有理数のステップ分の遅延を必要とすることになり、これを直に実装することは困難である。代わりに、上述の相回転によって、周波数領域における対応するシフトを適用することができる。

【0045】

本発明の好適な実施形態によれば、コンバータは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式を、Ｐ（ｚ）から導出されるスペクトルが厳密に実数であるように、かつＱ（ｚ）から導出されるスペクトルが厳密に虚数であるように、ハーフサンプルを用いて周波数領域へフーリエ変換するフーリエ変換デバイスを備えている。

【0046】

ある代替例は、ハーフサンプルを有するＤＦＴを実装するものである。具体的には、従来のＤＦＴを、

で定義できるのに対して、ハーフサンプルＤＦＴは、次式のように定義することができる。

【0047】

この公式化に対しては、ＦＦＴのような高速実装を容易に考案することができる。

【0048】

この公式化の利点は、いまや対称点が通常のｎ＝１ではなく、ｎ＝１／２に存在することにある。よって、このハーフサンプルＤＦＴを用いれば、数列、
［２，１，０，０，１，２］
で、実数値のフーリエスペクトルが得られる。

【0049】

よって、ｍ＋ｌが奇数である場合、係数ｐ_０，ｐ_１，ｐ_２，ｐ_２，ｐ_１，ｐ_０を有する多項式Ｐ（ｚ）に関して、ハーフサンプルＤＦＴ及びゼロパディングを用いれば、入力数列が、
［ｐ_２，ｐ_１，ｐ_０，０，０．．．０，ｐ_０，ｐ_１，ｐ_２］
のとき、実数値のスペクトルを求めることができる。

【0050】

同様に、多項式Ｑ（ｚ）の場合、数列、
［−ｑ_２，−ｑ_１，−ｑ_０，０，０．．．０，ｑ_０，ｑ_１，ｑ_２］
にハーフサンプルＤＦＴを適用して、純虚数スペクトルを求めることができる。

【0051】

これらの方法により、ｍとｌの如何なる組合せに対しても、多項式Ｐ（ｚ）の実数値スペクトルと任意のＱ（ｚ）の純虚数スペクトルを求めることができる。実際には、Ｐ（ｚ）とＱ（ｚ）のスペクトルが各々純実数と純虚数であることから、これらを単一の複素スペクトルに格納することができ、よって、これは、Ｐ（ｚ）＋Ｑ（ｚ）＝２Ａ（ｚ）のスペクトルに一致する。係数２によるスケーリングでは、根の位置に変化がなく、よってこれを無視することができる。したがって、単一のＦＦＴを用いてＡ（ｚ）のスペクトルのみを評価すれば、Ｐ（ｚ）とＱ（ｚ）のスペクトルを求めることができる。先に説明したように、Ａ（ｚ）の係数に循環シフトを適用するだけでよい。

【0052】

例えば、ｍ＝４でありｌ＝０である場合、Ａ（ｚ）の係数は、
［ａ_０，ａ_１，ａ_２，ａ_３，ａ_４］
であり、これは、ゼロパディングによって、次のように任意の長さＮにすることができる。
［ａ_０，ａ_１，ａ_２，ａ_３，ａ_４，０，０．．．０］

【0053】

次に、（ｍ＋ｌ）／２＝２ステップの循環シフトを適用すれば、次のようになる。
［ａ_２，ａ_３，ａ_４，０，０．．．０，ａ_０，ａ_１］

【0054】

この数列のＤＦＴを採用すれば、スペクトルの実数部及び複素数部にＰ（ｚ）とＱ（ｚ）のスペクトルが存在する。

【0055】

本発明の好適な一実施形態によれば、コンバータは、多項式Ｐ（ｚ）とＱ（ｚ）から合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を確立するように構成された合成多項式フォーマを備えている。

【0056】

本発明の好適な一実施形態によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）からの厳密な虚数スペクトルが、例えば高速フーリエ変換（ＦＦＴ）により合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を変換することによって、単一のフーリエ変換により確立されるように構成されている。

【0057】

多項式Ｐ（ｚ）とＱ（ｚ）は、各々、ｚ^{−（ｍ＋ｌ）／２}における対称軸と対称でありかつ逆対称である。よって、単位円ｚ＝ｅｘｐ（ｉθ）上で評価されるｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）とｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）のスペクトルは、各々、実数値及び複素数値ということになる。ゼロは、単位円上に存在することから、零交差を探索すれば、ゼロを発見することができる。さらに、単位円上の評価は、単純に高速フーリエ変換によって実装することができる。

【0058】

ｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）とｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）に対応するスペクトルは、各々実数及び複素数であることから、これらを単一の高速フーリエ変換によって実装することができる。具体的には、和ｚ^{−（ｍ＋ｌ）／２}（Ｐ（ｚ）＋Ｑ（ｚ））を求めれば、スペクトルの実数部と複素数部は、各々、ｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）とｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）に対応する。さらに、
ｚ^{−（ｍ＋ｌ）／２}（Ｐ（ｚ）＋Ｑ（ｚ））＝２ｚ^{−（ｍ＋ｌ）／２}Ａ（ｚ） （４）
であることから、Ｐ（ｚ）とＱ（ｚ）を明示的に決定することなく、２ｚ^{−（ｍ＋ｌ）／２}Ａ（ｚ）のＦＦＴを直に取り込んでｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）とｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）に対応するスペクトルを求めることができる。知りたいものはゼロの位置だけであるため、スカラ２による乗算を省略し、代わりにＦＦＴによるｚ^{−（ｍ＋ｌ）／２}Ａ（ｚ）を評価することができる。Ａ（ｚ）は、ｍ＋１個の非ゼロ係数しか持たないことから、ＦＦＴ枝刈りを用いて複雑性を低減できることを観察されたい［１１］。全ての根が発見されることを保証するためには、２つのゼロ毎にその間の少なくとも１つの周波数上でスペクトルが評価されるに足る長い長さＮのＦＦＴを用いなければならない。

【0059】

本発明の好適な一実施形態によれば、コンバータは、多項式Ｐ（ｚ）とＱ（ｚ）のスペクトルの数値範囲を、多項式Ｐ（ｚ）とＱ（ｚ）又は多項式Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイスを備えている。ここで、フィルタ多項式Ｂ（ｚ）は、対称性であって、単位円上に根を持たない。

【0060】

音声コーデックは、限定リソースを有するモバイルデバイス上で実装されることが多く、よって、数値演算は、固定小数点表示によって実装されなければならない。したがって、実装されるアルゴリズムは、その範囲が限定される数値表現を用いて動作することが極めて重要である。しかしながら、一般的な音声スペクトル包絡線の場合、フーリエスペクトルの数値範囲があまりに大きく、よって、零交差位置の保持を保証するためには、ＦＦＴの３２ビット実装が必要である。

【0061】

一方で、１６ビットのＦＦＴは、より低い複雑性で実装可能である場合が多く、よって、この１６ビット範囲内に適合するようにスペクトル値の範囲を限定することが有益であると考えられる。方程式｜Ｐ（ｅ^ｉθ）｜≦２｜Ａ（ｅ^ｉθ）｜と｜Ｑ（ｅ^ｉθ）｜≦２｜Ａ（ｅ^ｉθ）｜から、Ｂ（ｚ）Ａ（ｚ）の数値範囲を限定することによってＢ（ｚ）Ｐ（ｚ）とＢ（ｚ）Ｑ（ｚ）の数値範囲も限定されることが分かる。Ｂ（ｚ）が単位円上にゼロを持たなければ、Ｂ（ｚ）Ｐ（ｚ）及びＢ（ｚ）Ｑ（ｚ）は、単位円上にＰ（ｚ）及びＱ（ｚ）と同じ零交差を有する。さらに、Ｂ（ｚ）は、ｚ^{−（ｍ＋ｌ＋ｎ）／２}Ｐ（ｚ）Ｂ（ｚ）及びｚ^{−（ｍ＋ｌ＋ｎ）／２}Ｑ（ｚ）Ｂ（ｚ）がそれぞれ対称性と逆対称を保ち、かつそのスペクトルが純実数及び純虚数であるように、対称性でなければならない。したがって、ｚ^{（ｎ＋ｌ）／２}Ａ（ｚ）のスペクトルを評価する代わりに、ｚ^{（ｎ＋ｌ＋ｎ）／２}Ａ（ｚ）Ｂ（ｚ）を評価することができる。ここで、Ｂ（ｚ）は、単位円上に根を持たない次数ｎの対称多項式である。言い替えれば、先に述べたものと同じ手法を適用することができるが、まずは、Ａ（ｚ）をフィルタＢ（ｚ）で乗算し、次に修正した位相シフトｚ^{−（ｍ＋ｌ＋ｎ）／２}を適用する。

【0062】

残りのタスクは、Ａ（ｚ）Ｂ（ｚ）の数値範囲が、Ｂ（ｚ）は必ず対称性でありかつ単位円上に根を持たないという制約によって限定されるように、フィルタＢ（ｚ）を設計することである。この要件を満たす最も単純なフィルタは、次数２の線形位相フィルタ、
Ｂ_１（ｚ）＝β_０＋β_１ｚ^−１＋β_２ｚ^−２ （５）
である。ここで、β_ｋ∈Ｒはパラメータであり、かつ│β_２│＞２│β_１│である。β_ｋを調整すれば、スペクトル傾斜を修正することができ、積Ａ（ｚ）Ｂ_１（ｚ）の数値範囲を縮小することができる。計算上極めて効率的な手法は、βを、０周波数及びナイキストにおける振幅が等しくなるように、即ち│Ａ（１）Ｂ_１（１）│＝│Ａ（−１）Ｂ_１（−１）│であるように選択することであり、例えば、次式を選ぶことができる。
β_０＝Ａ（１）−Ａ（−１）、
β_１＝２（Ａ（１）＋Ａ（−１）） （６）

【0063】

この手法は、略平坦なスペクトルを提供する。

【0064】

Ａ（ｚ）が高域通過特性を有するのに対して、Ｂ_１（ｚ）は低域通過であり（図５も参照されたい）、積Ａ（ｚ）Ｂ_１（ｚ）は、予期されるように、０周波数及びナイキスト周波数において同じ振幅を有し、かつこれは、多かれ少なかれ平坦である。Ｂ_１（ｚ）は１自由度しか持たないので、積が完全に平坦になると予期できないことは明らかである。さらに、Ｂ_１（ｚ）Ａ（ｚ）の最高ピークと最低の谷との比は、Ａ（ｚ）のそれより遙かに小さいものであり得ることを観察されたい。これは、Ｂ_１（ｚ）Ａ（ｚ）の数値範囲がＡ（ｚ）のそれより遙かに小さいという望ましい効果が達成されていることを意味する。

【0065】

第２の、僅かに複雑性の高い方法は、Ａ（０．５ｚ）のインパルス応答の自己相関ｒ_ｋを計算するものである。ここでは、０．５を乗算することによってＡ（ｚ）のゼロが原点方向へ移動し、それによってスペクトル振幅が約半分に縮小される。自己相関ｒ_ｋにレビンソン−ダービンを適用すれば、最小位相である次数ｎのフィルタＨ（ｚ）が得られる。次に、Ｂ_２（ｚ）＝ｚ^−ｎＨ（ｚ）Ｈ（ｚ^−１）を定義して、略一定である│Ｂ_２（ｚ）Ａ（ｚ）│を求めることができる。│Ｂ２（ｚ）Ａ（ｚ）│の範囲が│Ｂ_１（ｚ）Ａ（ｚ）│のそれより小さいことに留意され得る。Ｂ（ｚ）を設計するためのさらなる手法は、ＦＩＲ設計に関する古典的文献［１８］において容易に見出すことができる。

【0066】

本発明の好適な一実施形態によれば、コンバータは、長くされた多項式Ｐ_ｅ（ｚ）とＱ_ｅ（ｚ）又は長くされた多項式Ｐ_ｅ（ｚ）とＱ_ｅ（ｚ）から導出される１つ又は複数の多項式のスペクトルの数値範囲を、長くされた多項式Ｐ_ｅ（ｚ）とＱ_ｅ（ｚ）にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイスを備えている。ここで、フィルタ多項式Ｂ（ｚ）は、対称性であって、単位円上に根を持たない。Ｂ（ｚ）は、先に説明した通りに求めることができる。

【0067】

さらなる態様では、この課題は、情報信号を符号化するための情報エンコーダを作動するための方法によって解決される。本方法は、
予測多項式Ａ（ｚ）の線形予測係数を求めるべく情報信号を分析するステップと、
予測多項式Ａ（ｚ）の線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現の周波数値ｆ_１．．．ｆ_ｎへ変換するステップであって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ^−ｍ−ｌＡ（ｚ^−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ^−ｍ−ｌＡ（ｚ^−１）
ｍは予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
のように定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することによって前記周波数値ｆ_１．．．ｆ_ｎを決定し、、Ｐ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）からの厳密な虚数スペクトルを確立しかつＰ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）から導出される厳密な虚数スペクトルのゼロを識別することによって、前記周波数値ｆ_１．．．ｆ_ｎを得るステップと、
前記周波数値ｆ_１．．．ｆ_ｎから量子化された周波数値ｆ_ｑ１．．．ｆ_ｑｎを得るステップと、
前記量子化された周波数値ｆ_ｑ１．．．ｆ_ｑｎを含むビットストリームを生成するステップと、を含む。

【0068】

さらに、本プログラムは、プロセッサ上で実行されることにより本発明による方法を実行するコンピュータプログラムによっても注目される。

【図面の簡単な説明】

【0069】

【図1】図１は、本発明による情報エンコーダの一実施形態を示す略図である。

【図2】図２は、Ａ（ｚ）、Ｐ（ｚ）及びＱ（ｚ）の例示的な関係を示す。

【図3】図３は、本発明による情報エンコーダのコンバータの第１の実施形態を示す略図である。

【図4】図４は、本発明による情報エンコーダのコンバータの第２の実施形態を示す略図である。

【図5】図５は、予測子Ａ（ｚ）、対応する平坦化フィルタＢ_１（ｚ）及びＢ_２（ｚ）及び積Ａ（ｚ）Ｂ_１（ｚ）及びＡ（ｚ）Ｂ_２（ｚ）の例示的な振幅スペクトルを示す。

【図6】図６は、本発明による情報エンコーダのコンバータの第３の実施形態を示す略図である。

【図7】図７は、本発明による情報エンコーダのコンバータの第４の実施形態を示す略図である。

【図8】図８は、本発明による情報エンコーダのコンバータの第５の実施形態を示す略図である。

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

【0070】

続いて、添付の図面を参照し、本発明の好適な実施形態について論じる。

【0071】

図１は、本発明による情報エンコーダ１の一実施形態を示す略図である。

【0072】

情報信号ＩＳを符号化するための情報エンコーダ１は、
予測多項式Ａ（ｚ）の線形予測係数を求めるべく情報信号ＩＳを分析するアナライザ２と、
予測多項式Ａ（ｚ）の線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現ＲＥＳ、ＩＥＳの周波数値ｆ_１．．．ｆ_ｎへ変換するコンバータ３であって、
Ｐ（ｚ）＝Ａ（ｚ）＋ｚ^−ｍ−ｌＡ（ｚ^−１）、及び、
Ｑ（ｚ）＝Ａ（ｚ）−ｚ^−ｍ−ｌＡ（ｚ^−１）、
ｍは予測多項式Ａ（ｚ）の次数、ｌはゼロ以上、
のように定義される多項式対Ｐ（ｚ）及びＱ（ｚ）を分析することにより、前記周波数値ｆ_１．．．ｆ_ｎを決定するように構成され、前記コンバータ３は、Ｐ（ｚ）から導出される厳密な実数スペクトルＲＥＳ及びＱ（ｚ）からの厳密な虚数スペクトルＩＥＳを確立しかつＰ（ｚ）から導出される厳密な実数スペクトルＲＥＳ及びＱ（ｚ）から導出される厳密な虚数スペクトルＩＥＳのゼロを識別することによって、前記周波数値ｆ_１．．．ｆ_ｎを得るように構成されたコンバータ３と、
前記周波数値ｆ_１．．．ｆ_ｎから量子化された周波数値ｆ_ｑ１．．．ｆ_ｑｎを得る量子化器４と、
前記量子化された周波数値ｆ_ｑ１．．．ｆ_ｑｎを含むビットストリームＢＳを生成するビットストリーム生成器５と、を備えている。

【0073】

本発明による情報エンコーダ１は、零交差探索を用いるのに対して、従来技術による根を求めるためのスペクトル手法は、振幅スペクトルにおける谷の発見に依存する。しかしながら、谷を探索する場合の確度は、零交差を探索する場合より劣る。例えば、数列［４，２，１，２，３］を考察されたい。明らかに、最小値は、第３の要素であり、ゼロは、第２の要素と第４の要素との間のどこかに存在することになる。言い替えれば、ゼロが第３の要素の右側に存在するか、左側に存在するかを決定することはできない。しかしながら、数列［４，２，１，−２，−３］について考察すれば、直ちに、零交差が第３の要素と第４の要素との間に存在することが分かり、許容誤差が半減する。振幅スペクトル手法の場合、零交差探索の場合と同じ確度を達成するには、分析点の数を２倍にする必要があることになる。

【0074】

振幅│Ｐ（ｚ）│と│Ｑ（ｚ）│の評価を比較すると、零交差手法は、確度において著しい優位点を有する。例えば、数列３，２，−１，−２を考察されたい。零交差手法では、ゼロは、２と−１との間に存在することが明らかである。しかしながら、対応する振幅列３，２，１，２を考察する場合、ゼロは、第２の要素と最後の要素との間のどこかに存在する、としか結論することができない。言い替えれば、零交差手法の場合の確度は、振幅ベースの手法の２倍である。

【0075】

さらに、本発明による情報エンコーダは、ｍ＝１２８等の長い予測子を用いてもよい。これに対して、チェビシェフ変換は、Ａ（ｚ）の長さが比較的短い、例えばｍ≦２０である場合にしか十分に機能しない。長い予測子の場合、チェビシェフ変換は、数値的に不安定であり、アルゴリズムの実用的実装は不可能である。

【0076】

したがって、提案する情報エンコーダ１の主要な特性は、零交差が探索されることに起因するとともに、時間領域から周波数領域への変換が行われ、極めて低い計算複雑性によりゼロを発見し得ることに起因して、チェビシェフベースの方法と同じく高い確度、又はより優れた確度を達成し得ることにある。

【0077】

結果的に、本発明による情報エンコーダ１は、ゼロ（根）をより正確に決定するだけでなく、低い計算複雑性で決定する。

【0078】

本発明による情報エンコーダ１は、数列の線スペクトルを決定する必要があるあらゆる信号処理アプリケーションにおいて使用可能である。本明細書では、例示的に、情報エンコーダ１を音声符号化の側面で論じる。本発明は、音声信号又は汎用オーディオ信号等の入力信号から線スペクトルを決定するための方法を必要とし、かつ入力信号がデジタルフィルタ又は他の数列として表現される、スペクトル振幅包絡線、知覚周波数マスキングしきい値、時間振幅包絡線、知覚時間マスキングしきい値、又は他の包絡線形状、又は線スペクトルを用いて包絡線情報を表現する、符号化、分析又は処理のための自己相関信号等の包絡線形状と等価の他の表現、をモデリングするための線形予測子を使用する、音声、オーディオ及び／又はビデオ符号化デバイス又はアプリケーションにおいて適用可能である。

【0079】

情報信号ＩＳは、例えば、オーディオ信号又はビデオ信号であってもよい。

【0080】

図２は、Ａ（ｚ）、Ｐ（ｚ）及びＱ（ｚ）の例示的な関係を示す。垂直な点線は、周波数値ｆ_１．．．ｆ_６を描いている。零交差が常に見えるように、振幅は、デシベルスケールではなく直線軸上に表されていることに留意されたい。線スペクトル周波数は、Ｐ（ｚ）とＱ（ｚ）の零交差において発生することが分かる。さらに、Ｐ（ｚ）及びＱ（ｚ）の振幅は、どこにおいても２│Ａ（ｚ）│以下、即ち、｜Ｐ（ｅ^ｉθ）｜≦２｜Ａ（ｅ^ｉθ）｜及び｜Ｑ（ｅ^ｉθ）｜≦２｜Ａ（ｅ^ｉθ）｜である。

【0081】

図３は、本発明による情報エンコーダのコンバータの第１の実施形態を示す略図である。

【0082】

本発明の好適な一実施形態によれば、コンバータ３は、予測多項式Ａ（ｚ）から多項式Ｐ（ｚ）とＱ（ｚ）を決定する決定デバイス６を備えている。

【0083】

本発明の好適な一実施形態によれば、コンバータは、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式を周波数領域へフーリエ変換するフーリエ変換デバイス８と、Ｐ（ｚ）から導出されるスペクトルＲＥＳの位相をそれが厳密に実数であるように調整するとともに、Ｑ（ｚ）から導出されるスペクトルＩＥＳの位相をそれが厳密に虚数であるように調整する調整デバイス７と、を備えている。フーリエ変換デバイス８は高速フーリエ変換を基礎とするものであってもよいし、離散型フーリエ変換を基礎とするものであってもよい。

【0084】

本発明の好適な一実施形態によれば、調整デバイス７は、多項式対Ｐ（ｚ）とＱ（ｚ）又は多項式対Ｐ（ｚ）とＱ（ｚ）から導出される１つ又は複数の多項式の係数を循環シフトする係数シフタ７として構成されている。

【0085】

本発明の好適な一実施形態によれば、係数シフタ７は、係数列の原初の中間点がその数列の第１のポジションへシフトされるような方法で、係数を循環シフトするように構成されている。

【0086】

理論的に、対称数列のフーリエ変換が実数値をとり、かつ逆対称の数列が純虚数フーリエスペクトルを有することが知られている。本事例における入力数列は、長さｍ＋ｌの多項式Ｐ（ｚ）又はＱ（ｚ）の係数であるが、これより遙かに長い長さＮ＞＞（ｍ＋ｌ）の離散型フーリエ変換を有する方が好まれると考えられる。より長いフーリエスペクトルを生成する従来的手法は、入力信号のゼロパディングである。しかし、数列のゼロパディングは、対称性が保たれるように慎重に実装されなければならない。

【0087】

まず、係数、
［ｐ_０，ｐ_１，ｐ_２，ｐ_１，ｐ_０］
を有する多項式Ｐ（ｚ）について考察する。

【0088】

高速フーリエ変換アルゴリズムが通常適用される方法では、対称点が第１の要素であるということが要求され、それによって、例えばＭＡＴＬＡＢに適用する場合、実数値の出力を求めるために
ｆｆｔ（［ｐ_２，ｐ_１，ｐ_０，ｐ_０，ｐ_１］）
のように書くことができる。具体的には、循環シフトが適用されてもよく、中間点要素に対応する対称点、すなわち係数ｐ_２は、第１のポジションになるように左へシフトされる。次に、ｐ_２の左側にあった係数は、数列の終わりに付加される。

【0089】

ゼロパディングされた数列、
［ｐ_０，ｐ_１，ｐ_２，ｐ_１，ｐ_０，０，０．．．０］
の場合も、同じプロセスを適用することができる。よって、数列、
［ｐ_２，ｐ_１，ｐ_０，０，０．．．０，ｐ_０，ｐ_１］
が実数値の離散型フーリエ変換を有する。ここで、入力数列におけるゼロの数は、スペクトルの所望される長さをＮとすれば、Ｎ−ｍ−ｌである。

【0090】

同様に、多項式Ｑ（ｚ）に対応する係数、
［ｑ_０，ｑ_１，０，−ｑ_１，−ｑ_０］
について考察されたい。前述の中間点が第１のポジションにくるように循環シフトを適用すれば、
［０，−ｑ_１，−ｑ_０，ｑ_０，ｑ_１］
が得られ、これは、純虚数の離散型フーリエ変換を有する。次には、この数列にゼロパディング変換を行なうことができる。
［０，−ｑ_１，−ｑ_０，０，０．．．０，ｑ_０，ｑ_１］

【0091】

上記は、数列の長さが奇数である事例にのみ当てはまり、よって、ｍ＋ｌは偶数であることに留意されたい。ｍ＋ｌが奇数である事例に関しては、２つのオプションがある。すなわち、周波数領域において循環シフトを実装することができ、あるいは、ハーフサンプルでＤＦＴを適用することができる。

【0092】

本発明の好適な実施形態によれば、コンバータ３は、Ｐ（ｚ）から導出される厳密な実数スペクトルＲＥＳとＱ（ｚ）から導出される厳密な虚数スペクトルＩＥＳのゼロを識別するためのゼロ識別子９を備えている。

【0093】

本発明の好適な一実施形態によれば、ゼロ識別子９は、
ａ）実数スペクトルＲＥＳをヌル周波数において開始することと、
ｂ）実数スペクトルＲＥＳにおいて符号変化が発見されるまで周波数を増大することと、
ｃ）虚数スペクトルＩＥＳにおいてさらなる符号変化が発見されるまで周波数を増大することと、
ｄ）全てのゼロが発見されるまでステップｂ）及びステップｃ）を反復すること、によってゼロを識別するように構成されている。

【0094】

Ｑ（ｚ）、ひいてはスペクトルの虚数部ＩＥＳは、ヌル周波数において常にゼロを有することに留意されたい。根は重なり合うことから、Ｐ（ｚ）、ひいてはスペクトルの実数部ＲＥＳは、ヌル周波数において常に非ゼロになる。したがって、ヌル周波数における実数部から開始して、最初の零交差、ひいては最初の周波数値ｆ_１を示す最初の符号変化が発見されるまで、周波数を増大させることができる。

【0095】

根はインタレースされることから、Ｑ（ｚ）のスペクトルＩＥＳは、次の符号変化を有する。したがって、Ｑ（ｚ）のスペクトルＩＥＳの符号変化が発見されるまで、周波数を増大させることができる。次に、このプロセスは、全ての周波数値ｆ_１．．．ｆ_ｎが発見されるまで、Ｐ（ｚ）とＱ（ｚ）のスペクトル間で交互に反復されてもよい。したがって、スペクトルＲＥＳとＩＥＳにおける零交差の位置決めに使用される手法は、チェビシェフ領域において適用される手法に類似する［６、７］。

【0096】

Ｐ（ｚ）とＱ（ｚ）のゼロはインタレースされることから、１パスで全てのゼロを発見するように実数部ＲＥＳと複素数部ＩＥＳ上のゼロの探索を交互に行い、複雑性を完全探索の半分に減らすことができる。

【0097】

本発明の好適な一実施形態によれば、ゼロ識別子９は補間によってゼロを識別するように構成されている。

【0098】

零交差手法に加えて、ゼロのポジションを、例えば［７］のような従来方法において行われるようなより高い確度で推定できるように、補間を容易に適用することができる。

【0099】

図４は、本発明による情報エンコーダ１のコンバータ３の第２の実施形態を示す略図である。

【0100】

本発明の好適な一実施形態によれば、コンバータ３は、長くされた多項式対Ｐ_ｅ（ｚ）とＱ_ｅ（ｚ）を生成するために多項式Ｐ（ｚ）とＱ（ｚ）に対して値「０」を有する１つ又は複数の係数を加えるゼロパディング・デバイス１０を備えている。確度は、評価されるスペクトルＲＥＳ、ＩＥＳの長さを拡大することによってさらに向上させることができる。システムに関する情報を基礎として、事例によっては、実質上、周波数値ｆ_１．．．ｆ_ｎ間の最小距離を決定し、ひいては、全ての周波数値ｆ_１．．．ｆ_ｎをそれで発見することができるスペクトルＲＥＳ、ＩＥＳの最小長さを決定することが可能である［８］。

【0101】

本発明の好適な一実施形態によれば、コンバータ３は、線形予測係数を予測多項式Ａ（ｚ）のスペクトル周波数表現ＲＥＳ、ＩＥＳの周波数値ｆ_１．．．ｆ_ｎへ変換する間に、長くされた多項式Ｐ_ｅ（ｚ）とＱ_ｅ（ｚ）の値「０」を有することが知られている係数を用いる演算の少なくとも一部が省略されるように構成される。

【0102】

しかしながら、スペクトルの長さが増せば、計算複雑性も増大する。複雑性の最大の誘因は、Ａ（ｚ）の係数の、高速フーリエ変換等の時間領域から周波数領域への変換である。しかしながら、これは、係数ベクトルが所望される長さにゼロパディングされていることから、極めて疎である。複雑性を低減するためにこの事実を容易に使用することができる。これは、どの係数がゼロであるかが正確に分かり、高速フーリエ変換の各反復において単純にゼロを包含する演算を省くことができるという意味において、どちらかといえば単純な問題である。このように疎な高速フーリエ変換の適用は単純であり、当業者であるプログラマは誰でもこれを実装することができる。このような実装の複雑性は、Ｏ（Ｎｌｏｇ_２（１＋ｍ＋ｌ））である。ここで、Ｎはスペクトルの長さ、ｍ及びｌは先に定義した通りである。

【0103】

本発明の好適な一実施形態によれば、コンバータは、長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）又は長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）から導出される１つ又は複数の多項式のスペクトルの数値範囲を、長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイス１１を備えている。ここで、フィルタ多項式Ｂ（ｚ）は対称性であって単位円上に根を持たない。Ｂ（ｚ）は、先に説明した通りに求めることができる。

【0104】

図５は、予測子Ａ（ｚ）、対応する平坦化フィルタＢ_１（ｚ）とＢ_２（ｚ）及び積Ａ（ｚ）Ｂ_１（ｚ）とＡ（ｚ）Ｂ_２（ｚ）の例示的な振幅スペクトルを示す。水平の点線は、０周波数とナイキスト周波数におけるＡ（ｚ）Ｂ_１（ｚ）のレベルを示す。

【0105】

本発明の好適な一実施形態（不図示）によれば、コンバータ３は、多項式Ｐ（ｚ）及びＱ（ｚ）のスペクトルＲＥＳ、ＩＥＳの数値範囲を、多項式Ｐ（ｚ）及びＱ（ｚ）又は多項式Ｐ（ｚ）及びＱ（ｚ）から導出される１つ又は複数の多項式にフィルタ多項式Ｂ（ｚ）を乗算することによって制限する制限デバイス１１を備えている。ここで、フィルタ多項式Ｂ（ｚ）は対称性であって、単位円上に根を持たない。

【0106】

音声コーデックは、限定リソースを有するモバイルデバイス上で実装されることが多く、数値演算は固定小数点表示によって実装されなければならない。したがって、実装されるアルゴリズムは、範囲が限定される数値表現を用いて動作することが極めて重要である。しかしながら、一般的な音声スペクトル包絡線の場合、フーリエスペクトルの数値範囲があまりに大きく、零交差位置の保持を保証するためには、ＦＦＴの３２ビット実装が必要である。

【0107】

一方で、１６ビットのＦＦＴはより低い複雑性で実装可能である場合が多く、この１６ビット範囲内に適合するようにスペクトル値の範囲を限定することが有益であると考えられる。方程式｜Ｐ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜及び｜Ｑ（ｅｉθ）｜≦２｜Ａ（ｅｉθ）｜からは、Ｂ（ｚ）Ａ（ｚ）の数値範囲を限定することにより、Ｂ（ｚ）Ｐ（ｚ）及びＢ（ｚ）Ｑ（ｚ）の数値範囲も限定されることが分かる。Ｂ（ｚ）が単位円上にゼロを持たなければ、Ｂ（ｚ）Ｐ（ｚ）及びＢ（ｚ）Ｑ（ｚ）は、単位円上にＰ（ｚ）及びＱ（ｚ）と同じ零交差を有する。さらに、Ｂ（ｚ）は、ｚ^{−（ｍ＋ｌ＋ｎ）／２}Ｐ（ｚ）Ｂ（ｚ）及びｚ^{−（ｍ＋ｌ＋ｎ）／２}Ｑ（ｚ）Ｂ（ｚ）が各々対称性及び逆対称を保ちかつそのスペクトルが純実数及び純虚数であるように、対称でなければならない。したがって、ｚ^{（ｎ＋ｌ）／２}Ａ（ｚ）のスペクトルを評価する代わりに、ｚ^{（ｎ＋ｌ＋ｎ）／２}Ａ（ｚ）Ｂ（ｚ）を評価することができる。ここで、Ｂ（ｚ）は、単位円上に根を持たない次数ｎの対称多項式である。言い替えれば、先に述べたものと同じ手法を適用することができるが、まずは、Ａ（ｚ）をフィルタＢ（ｚ）で乗算し、次に修正した位相シフトｚ^{−（ｍ＋ｌ＋ｎ）／２}を適用する。

【0108】

残りのタスクは、Ａ（ｚ）Ｂ（ｚ）の数値範囲が、Ｂ（ｚ）は必ず対称性でありかつ単位円上に根を持たないという制約によって限定されるように、フィルタＢ（ｚ）を設計することである。この要件を満たす最も単純なフィルタは、次数２の線形位相フィルタ、Ｂ_１（ｚ）＝β_０＋β_１ｚ^−１＋β_２ｚ^−２である。ここで、β_ｋ∈Ｒはパラメータであり、かつ│β_２│＞２│β_１│である。β_ｋを調整すればスペクトル傾斜を修正することができ、積Ａ（ｚ）Ｂ_１（ｚ）の数値範囲を縮小することができる。計算上極めて効率的な手法は、βを０周波数及びナイキストにおける振幅が等しくなるように、すなわち│Ａ（１）Ｂ_１（１）│＝│Ａ（−１）Ｂ_１（−１）│であるように選択することであり、例えば、β_０＝Ａ（１）−Ａ（−１）及びβ_１＝２（Ａ（１）＋Ａ（−１））を選ぶことができる。

【0109】

この手法は、略平坦なスペクトルを提供する。

【0110】

Ａ（ｚ）が高域通過特性を有するのに対して、Ｂ_１（ｚ）は低域通過であり、積Ａ（ｚ）Ｂ_１（ｚ）は、予期されるように、０周波数及びナイキスト周波数において同じ振幅を有し、かつこれは多かれ少なかれ平坦であるということが図５から観察される。Ｂ_１（ｚ）は１自由度しか持たないので、積が完全に平坦になると予期できないことは明らかである。さらに、Ｂ_１（ｚ）Ａ（ｚ）の最高ピークと最低の谷との比は、Ａ（ｚ）のそれより遙かに小さいものであり得ることを観察されたい。これは、Ｂ_１（ｚ）Ａ（ｚ）の数値範囲がＡ（ｚ）のそれより遙かに小さいという望ましい効果が達成されていることを意味する。

【0111】

僅かに複雑性の高い、第２の方法は、Ａ（０．５ｚ）のインパルス応答の自己相関ｒ_ｋを計算するものである。ここでは、０．５を乗算することによってＡ（ｚ）のゼロが原点方向へ移動し、それによってスペクトル振幅が約半分に縮小される。自己相関ｒ_ｋにレビンソン−ダービンを適用すれば、最小位相である次数ｎのフィルタＨ（ｚ）が得られる。次には、Ｂ_２（ｚ）＝ｚ^−ｎＨ（ｚ）Ｈ（ｚ^−１）を定義して、略一定である│Ｂ_２（ｚ）Ａ（ｚ）│を求めることができる。│Ｂ２（ｚ）Ａ（ｚ）│の範囲は、│Ｂ_１（ｚ）Ａ（ｚ）│のそれより小さいことが留意されるであろう。Ｂ（ｚ）を設計するためのさらなる手法は、ＦＩＲ設計に関する古典的文献［１８］において容易に見出すことができる。

【0112】

図６は、本発明による情報エンコーダ１のコンバータ３の第３の実施形態を示す略図である。

【0113】

本発明の好適な一実施形態によれば、調整デバイス１２は、フーリエ変換デバイス８の出力の位相をシフトするための位相シフタ１２として構成されている。

【0114】

本発明の好適な一実施形態によれば、位相シフタ１２は、フーリエ変換デバイス８の出力の位相を、ｋ番目の周波数ビンにｅｘｐ（ｉ２πｋｈ／Ｎ）を乗算することによってシフトするように構成されている。ここで、Ｎはサンプルの長さであり、かつｈ＝（ｍ＋ｌ）／２である。

【0115】

時間領域における循環シフトが周波数領域における位相回転と同一であることはよく知られている。具体的には、時間領域におけるｈ＝（ｍ＋ｌ）／２ステップのシフトは、ｋ番目の周波数ビンとｅｘｐ（−ｉ２πｋｈ／Ｎ）との乗算に一致する。ただし、Ｎはスペクトルの長さである。したがって、循環シフトの代わりに周波数領域における乗算を適用すれば、全く同じ結果を得ることができる。ただし、この手法では、複雑性が僅かに増す。ｈ＝（ｍ＋ｌ）／２は、ｍ＋ｌが偶数であるときにのみ整数であることに留意されたい。ｍ＋ｌが奇数のとき、循環シフトは有理数のステップ分の遅延を必要とすることになり、これを直に実装ことは困難である。代わりに、上述の相回転によって、周波数領域における対応するシフトを適用することができる。

【0116】

図７は、本発明による情報エンコーダ１のコンバータ３の第４の実施形態を示す略図である。

【0117】

本発明の好適な一実施形態によれば、コンバータ３は、多項式Ｐ（ｚ）及びＱ（ｚ）から合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を確立するように構成された合成多項式フォーマ１３を備えている。

【0118】

本発明の好適な一実施形態によれば、コンバータ３は、Ｐ（ｚ）から導出される厳密な実数スペクトルとＱ（ｚ）からの厳密な虚数スペクトルが、単一のフーリエ変換により、例えば高速フーリエ変換（ＦＦＴ）により合成多項式Ｃ（Ｐ（ｚ），Ｑ（ｚ））を変換することによって確立されるように構成されている。

【0119】

多項式Ｐ（ｚ）及びＱ（ｚ）は、各々、ｚ^{−（ｍ＋ｌ）／２}における対称軸と対称であり逆対称である。よって、単位円ｚ＝ｅｘｐ（ｉθ）上で評価されるｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）及びｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）のスペクトルは、各々、実数値及び複素数値ということになる。ゼロは単位円上に存在することから、零交差を探索すればゼロを発見することができる。さらに、単位円上の評価は、単純に高速フーリエ変換によって実行することができる。

【0120】

ｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）及びｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）に対応するスペクトルは、各々実数と複素数であることから、これらを単一の高速フーリエ変換によって実行することができる。具体的には、和ｚ^{−（ｍ＋ｌ）／２}（Ｐ（ｚ）＋Ｑ（ｚ））を求めれば、スペクトルの実数部と複素数部は、各々、ｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）及びｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）に対応する。さらに、ｚ^{−（ｍ＋ｌ）／２}（Ｐ（ｚ）＋Ｑ（ｚ））＝２ｚ^{−（ｍ＋ｌ）／２}Ａ（ｚ）であることから、Ｐ（ｚ）及びＱ（ｚ）を明示的に決定することなく、２ｚ^{−（ｍ＋ｌ）／２}Ａ（ｚ）のＦＦＴを直に取り込んでｚ^{−（ｍ＋ｌ）／２}Ｐ（ｚ）及びｚ^{−（ｍ＋ｌ）／２}Ｑ（ｚ）に対応するスペクトルを求めることができる。知りたいものはゼロの位置のみであることから、スカラ２による乗算を省略して、代わりにＦＦＴによるｚ^{−（ｍ＋ｌ）／２}Ａ（ｚ）を評価することができる。Ａ（ｚ）は、ｍ＋１個の非ゼロ係数しか持たないことから、ＦＦＴ枝刈り（FFT pruning）を用いて複雑性を低減できることを観察されたい［１１］。全ての根が発見されることを保証するためには、２つのゼロ毎にその間の少なくとも１つの周波数上でスペクトルが評価されるに足る長い長さＮのＦＦＴを用いなければならない。

【0121】

本発明の好適な一実施形態（不図示）によれば、コンバータ３は、長くされた多項式Ｐ_ｅ（ｚ）及びＱ_ｅ（ｚ）から合成多項式Ｃ_ｅ（Ｐ_ｅ（ｚ），Ｑ_ｅ（ｚ））を確立するように構成された合成多項式フォーマを備えている。

【0122】

本発明の好適な一実施形態（不図示）によれば、コンバータは、Ｐ（ｚ）から導出される厳密な実数スペクトル及びＱ（ｚ）からの厳密な虚数スペクトルが、合成多項式Ｃ_ｅ（Ｐ_ｅ（ｚ），Ｑ_ｅ（ｚ））を変換することによって単一のフーリエ変換により確立されるように構成されている。

【0123】

図８は、本発明による情報エンコーダ１のコンバータ３の第５の実施形態を示す略図である。

【0124】

本発明の好適な実施形態によれば、コンバータ３は、多項式対Ｐ（ｚ）及びＱ（ｚ）又は多項式対Ｐ（ｚ）及びＱ（ｚ）から導出される１つ又は複数の多項式を、Ｐ（ｚ）から導出されるスペクトルが厳密に実数であり、かつＱ（ｚ）から導出されるスペクトルが厳密に虚数であるように、ハーフサンプルを用いて周波数領域へフーリエ変換するフーリエ変換デバイス１４を備えている。

【0125】

ある代替例は、ハーフサンプルを有するＤＦＴを設けたものである。具体的には、従来のＤＦＴが、

と定義されるのに対して、ハーフサンプルＤＦＴは、次式のように定義することができる。

【0126】

この公式化に対しては、ＦＦＴのような高速実装を容易に考案することができる。

【0127】

この公式化の利点は、対称点が通常のｎ＝１ではなく、ｎ＝１／２に存在することにある。よって、このハーフサンプルＤＦＴを用いれば、数列、
［２，１，０，０，１，２］
で、実数値のフーリエスペクトルＲＥＳが得られる。

【0128】

よって、ｍ＋ｌが奇数である場合、係数ｐ_０，ｐ_１，ｐ_２，ｐ_２，ｐ_１，ｐ_０を有する多項式Ｐ（ｚ）に関して、ハーフサンプルＤＦＴ及びゼロパディングを用いれば、次のような入力数列の場合に実数値スペクトルＲＥＳを求めることができる。
［ｐ_２，ｐ_１，ｐ_０，０，０．．．０，ｐ_０，ｐ_１，ｐ_２］

【0129】

同様に、多項式Ｑ（ｚ）の場合、数列、
［−ｑ_２，−ｑ_１，−ｑ_０，０，０．．．０，ｑ_０，ｑ_１，ｑ_２］
にハーフサンプルＤＦＴを適用して、純虚数スペクトルＩＥＳを求めることができる。

【0130】

これらの方法により、ｍ及びｌの如何なる組合せに対しても、多項式Ｐ（ｚ）の実数値スペクトル及び任意のＱ（ｚ）の純虚数スペクトルを求めることができる。実際には、Ｐ（ｚ）及びＱ（ｚ）のスペクトルが各々純実数と純虚数であることから、これらを単一の複素スペクトルに記憶することができ、これは、Ｐ（ｚ）＋Ｑ（ｚ）＝２Ａ（ｚ）のスペクトルに一致する。係数２によるスケーリングでは、根の位置に変化がなく、よってこれを無視することができる。したがって、単一のＦＦＴを用いてＡ（ｚ）のスペクトルのみを評価すれば、Ｐ（ｚ）及びＱ（ｚ）のスペクトルを求めることができる。先に説明したように、Ａ（ｚ）の係数に循環シフトを適用するだけでよい。

【0131】

例えば、ｍ＝４であり、ｌ＝０である場合、Ａ（ｚ）の係数は、
［ａ_０，ａ_１，ａ_２，ａ_３，ａ_４］
であり、これは、ゼロパディングによって、次のように任意の長さＮにすることができる。
［ａ_０，ａ_１，ａ_２，ａ_３，ａ_４，０，０．．．０］

【0132】

次に、（ｍ＋ｌ）／２＝２ステップの循環シフトを適用すれば、次のようになる。
［ａ_２，ａ_３，ａ_４，０，０．．．０，ａ_０，ａ_１］

【0133】

この数列のＤＦＴを採用すれば、スペクトルの実数部ＲＥＳと複素数部ＩＥＳにＰ（ｚ）とＱ（ｚ）のスペクトルが存在する。

【0134】

ｍ＋ｌが偶数の事例における全体的アルゴリズムは、次のように述べることができる。ａ_ｋで示されるＡ（ｚ）の係数は、長さＮのバッファ内に存在するものとする。

【0135】

１．左へ（ｍ＋ｌ）／２ステップのａ_ｋに循環シフトを適用する。

【0136】

２．数列ａ_ｋの高速フーリエ変換を計算し、これをＡ_ｋで示す。

【0137】

３．全ての周波数値が発見されるまで、ｋ＝０を起点とし、次の（ａ）、（ｂ）を交互に行う。
（ａ）符号（実数（Ａ_ｋ））＝符号（実数（Ａ_ｋ＋１））である間に、ｋ：＝ｋ＋１を増大する。零交差が発見されれば、ｋを周波数値のリストに記憶する。
（ｂ）符号（虚数（Ａ_ｋ））＝符号（虚数（Ａ_ｋ＋１））である間に、ｋ：＝ｋ＋１を増大する。零交差が発見されれば、ｋを周波数値のリストに記憶する。

【0138】

４．各周波数値について、Ａ_ｋとＡ_ｋ＋１との間を補間し、正確なポジションを決定する。

【0139】

ここで、関数符号（ｘ）、実数（ｘ）及び虚数（ｘ）は、各々、ｘの符号、ｘの実数部及びｘの虚数部を指す。

【0140】

ｍ＋ｌが奇数の事例では、循環シフトが僅かに左へ（ｍ＋ｌ−１）／２ステップ低減され、通常の高速フーリエ変換がハーフサンプル高速フーリエ変換によって置換される。

【0141】

あるいは、循環シフトと第１のフーリエ変換との組合せを、常に、高速フーリエ変換及び周波数領域における位相シフトと置換することもできる。

【0142】

より正確な根位置を求めるためには、上述の提案方法を用いて最初の推定を行い、次に、根軌跡を精緻化（refine）する第２のステップを適用することが可能である。精緻化（refinement）のためには、デュラン−ケルナー法、アバース−エールリッヒ法、ラゲールのガウス−ニュートン法又はその他等［１１−１７］のあらゆる古典的な多項式根探索方法を適用することができる。

【0143】

ある公式化において、提示の本方法は、次のようなステップより成る。

【0144】

（ａ）長さＮまでゼロパディングされた長さｍ＋ｌ＋１の数列、ただしｍ＋ｌは偶数、に対して、バッファ長さがＮとなって出力スペクトルの所望される長さに一致するように、左へ（ｍ＋ｌ）／２ステップの循環シフトを適用する、又は、
長さＮまでゼロパディングされた長さｍ＋ｌ＋１の数列、ただしｍ＋ｌは奇数、に対して、バッファ長さがＮとなって出力スペクトルの所望される長さに一致するように、左へ（ｍ＋ｌ−１）／２ステップの循環シフトを適用する。

【0145】

（ｂ）ｍ＋ｌが偶数であれば、数列に通常のＤＦＴを適用する。ｍ＋ｌが奇数であれば、方程式３又は等価の表現が記述しているように、数列にハーフサンプルＤＦＴを適用する。

【0146】

（ｃ）入力信号が対称又は逆対称であれば、周波数領域表現の零交差を探索し、その位置をリストに記憶する。

【0147】

入力信号が合成列Ｂ（ｚ）＝Ｐ（ｚ）＋Ｑ（ｚ）であれば、周波数領域表現の実数部と虚数部の双方において零交差を探索し、その位置をリストに記憶する。入力信号が、合成列Ｂ（ｚ）＝Ｐ（ｚ）＋Ｑ（ｚ）であり、かつＰ（ｚ）とＱ（ｚ）の根が交互するか、又はそれと同様の構造を有する場合には、周波数領域表現の実数部と虚数部を交互して零交差を探索し、その位置をリストに記憶する。

【0148】

別の公式化では、本提示方法は、次のようなステップより成る。

【0149】

（ａ）形式が先行点と同じである入力信号に対して、入力数列にＤＦＴを適用する。

【0150】

（ｂ）周波数領域値に相回転を適用する。これは、入力数列の左へ（ｍ＋ｌ）／２ステップによる循環シフトに等しい。

【0151】

（ｃ）先行点で行われたような零交差探索を適用する。

【0152】

以下は、記述した実施形態によるエンコーダ１及び方法に関連する言及である。

【0153】

幾つかの態様は、装置のコンテキストにおいて記述されているが、これらの態様が、対応する方法の説明をも表現し、ブロック又はデバイスが方法ステップ又は方法ステップの特徴に対応していることは明らかである。同様に、方法ステップのコンテキストにおいて記述されている態様は、対応する装置の対応するブロック又はアイテム又は特徴の説明をも表現している。

【0154】

所定の実装要件によっては、本発明の実施形態は、ハードウェアにおいてもソフトウェアにおいても実装することができる。その実装は、個々の方法が実行されるようにプログラム可能コンピュータシステムと共働する（又は、共働することができる）電子読取り可能制御信号を格納しているデジタル記憶媒体、例えばフロッピーディスク、ＤＶＤ、ＣＤ、ＲＯＭ、ＰＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ又はフラッシュメモリを用いて実行することができる。

【0155】

本発明による幾つかの実施形態は、本明細書に記述している方法のうちの１つが実行されるようにプログラム可能コンピュータシステムと共働することができる電子読取り可能制御信号を有するデータキャリアを備える。

【0156】

概して、本発明の実施形態は、プログラムコードを有するコンピュータ・プログラム・プロダクトとしての実装が可能であり、プログラムコードは、コンピュータ上でコンピュータ・プログラム・プロダクトが実行されると本発明方法のうちの１つを実行するように作動する。プログラムコードは、例えば、機械読取り可能キャリア上に記憶されてもよい。

【0157】

他の実施形態は、機械読取り可能キャリア、又は非一時的記憶媒体上に記憶される、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムを含む。

【0158】

したがって、言い替えれば、本発明方法の一実施形態は、コンピュータ上でコンピュータプログラムが実行されると本明細書に記述している方法のうちの１つを実行するためのプログラムコードを有するコンピュータプログラムである。

【0159】

したがって、本発明方法のさらなる実施形態は、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムを記憶して備えるデータキャリア（又はデジタル記憶媒体、若しくはコンピュータ読取り可能媒体）である。

【0160】

したがって、本発明方法のさらなる実施形態は、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムを表すデータストリーム又は信号シーケンスである。データストリーム又は信号シーケンスは、例えば、データ通信接続を介して、例えばインターネットを介して伝送されるように構成されてもよい。

【0161】

さらなる実施形態は、本明細書に記述している方法のうちの１つを実行するように構成され、又は実行するように調整される処理手段、例えばコンピュータ、又はプログラマブル論理デバイス、を含む。

【0162】

さらなる実施形態は、本明細書に記述している方法のうちの１つを実行するためのコンピュータプログラムをインストールしているコンピュータを含む。

【0163】

実施形態によっては、プログラム可能な論理デバイス（例えば、フィールド・プログラマブル・ゲート・アレイ）は、本明細書に記述している方法の機能のうちの一部又は全てを実行するために使用されてもよい。実施形態によっては、フィールド・プログラマブル・ゲート・アレイは、本明細書に記述している方法のうちの１つを実行するために、マイクロプロセッサと共働してもよい。概して、本方法は、効果的には、あらゆるハードウェア装置によって実行される。

【0164】

本発明を、幾つかの実施形態に関連して記述してきたが、本発明の範囲に含まれる変更、置換及び均等物が存在する。本発明の方法及び組成を実装する方法には、多くの代替方法が存在することも留意されるべきである。したがって、添付した以下のクレームの意図は、このような変更、置換及び均等物を本発明の精神及び範囲に包含されるものとして解釈することにある。

【符号の説明】

【0165】

１ 情報エンコーダ
２ アナライザ
３ コンバータ
４ 量子化器
５ ビットストリーム生成器
６ 決定デバイス
７ 係数シフタ
８ フーリエ変換デバイス
９ ゼロ識別子
１０ ゼロパディング・デバイス
１１ 制限デバイス
１２ 位相シフタ
１３ 合成多項式フォーマ
１４ ハーフサンプル・フーリエ変換デバイス
ＩＳ 情報信号
ＲＥＳ 実数スペクトル
ＩＥＳ 虚数スペクトル
ｆ_１．．．ｆ_ｎ 周波数値
ｆ_ｑ１．．．ｆ_ｑｎ 量子化された周波数値
ＢＳ ビットストリーム

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】