(57)【要約】スピーチ符号化アルゴリズムのコードブックベクトルを決定することによって、スピーチ信号を符号化する装置が提供される。この装置は、自己相関行列Ｒを決定する行列決定部（１１０）と、自己相関行列Ｒに依存してコードブックベクトルを決定するコードブックベクトル決定部（１２０）とを含む。行列決定部（１１０）は、ベクトルｒのベクトル係数を決定することによって、自己相関行列Ｒを決定するよう構成され、自己相関行列Ｒは複数の行と複数の列とを含み、ベクトルｒは自己相関行列Ｒの列の１つ又は行の１つを示し、ここでR(i , j) = r(|i− j|)であり、R(i , j)は自己相関行列Ｒの係数を示し、ｉは自己相関行列Ｒの複数の行の１つを示す第１インデックスであり、ｊは自己相関行列Ｒの複数の列の１つを示す第２インデックスである。【選択図】  図１
2015532456.tif 000002
【特許請求の範囲】
【請求項１】
スピーチ符号化アルゴリズムのコードブックベクトルを決定することによって、スピーチ信号を符号化する装置であって、
自己相関行列Ｒを決定する行列決定部（１１０）と、
前記自己相関行列Ｒに依存して前記コードブックベクトルを決定するコードブックベクトル決定部（１２０）とを含み、
前記行列決定部（１１０）は、ベクトルｒのベクトル係数を決定することによって、前記自己相関行列Ｒを決定するよう構成され、前記自己相関行列Ｒは複数の行と複数の列とを含み、前記ベクトルｒは前記自己相関行列Ｒの列の１つ又は行の１つを示し、
R(i , j) = r(|i− j|)
であり、ここでR(i,j)は前記自己相関行列Ｒの係数を示し、ｉは前記自己相関行列Ｒの複数の行の１つを示す第１インデックスであり、ｊは前記自己相関行列Ｒの複数の列の１つを示す第２インデックスである、装置。
【請求項２】
請求項１に記載の装置において、
前記行列決定部（１１０）は次式を適用することによって前記ベクトルｒのベクトル係数を決定するよう構成され、
2015532456.tif 000060
ここで、ｈ（ｋ）は線形予測モデルの知覚的に重み付けられたインパルス応答を示し、ｋは整数であるインデックスであり、ｌは整数であるインデックスである、装置。
【請求項３】
請求項１又は２に記載の装置において、
前記行列決定部（１１０）は知覚的に重み付けられた線形予測子に依存して、前記自己相関行列Ｒを決定するよう構成されている、装置。
【請求項４】
請求項１乃至３のいずれかに記載の装置において、
前記コードブックベクトル決定部（１２０）は、次式を適用することにより前記コードブックベクトルを決定するよう構成され、
2015532456.tif 000061
ここで、Ｒは自己相関行列であり、
2015532456.tif 000062
は前記スピーチ符号化アルゴリズムのコードブックベクトルの１つであり、
2015532456.tif 000063
は正規化された相関である、装置。
【請求項５】
請求項４に記載の装置において、
前記コードブックベクトル決定部（１２０）は、前記正規化された相関
2015532456.tif 000064
を最大化するスピーチ符号化アルゴリズムのコードブックベクトル
2015532456.tif 000065
を決定するよう構成されている、装置。
【請求項６】
請求項１乃至５のいずれかに記載の装置において、
前記コードブックベクトル決定部（１２０）は、行列分解を行うことによって前記自己相関行列Ｒを分解するよう構成されている、装置。
【請求項７】
請求項６に記載の装置において、
前記コードブックベクトル決定部（１２０）は、前記コードブックベクトルを決定するための対角行列Ｄを決定するために、前記行列分解を行うよう構成されている、装置。
【請求項８】
請求項７に記載の装置において、
前記コードブックベクトル決定部（１２０）は、次式を用いることによって前記コードブックベクトルを決定するよう構成され、
2015532456.tif 000066
ここで、Ｄは前記対角行列であり、ｆは第１ベクトルであり、
2015532456.tif 000067
は第２ベクトルである、装置。
【請求項９】
請求項７又は８に記載の装置において、
前記コードブックベクトル決定部（１２０）は、前記自己相関行列Ｒに対してファンデルモンド因数分解を行って前記自己相関行列Ｒを分解するよう構成され、かつ前記行列分解を行って前記コードブックベクトルを決定するための前記対角行列Ｄを決定するよう構成されている、装置。
【請求項１０】
請求項７乃至９のいずれかに記載の装置において、
前記コードブックベクトル決定部（１２０）は、前記コードブックベクトルを決定するために次式を用いるよう構成され、
2015532456.tif 000068
ここで、Ｃは畳み込み行列を示し、Ｖはフーリエ変換を示し、ｘはスピーチ信号を示す、装置  。
【請求項１１】
請求項７乃至１０のいずれかに記載の装置において、
前記コードブックベクトル決定部（１２０）は、前記自己相関行列Ｒに対して特異値分解を行って前記自己相関行列Ｒを分解するよう構成され、かつ前記行列分解を行って前記コードブックベクトルを決定するための前記対角行列Ｄを決定するよう構成されている、装置。
【請求項１２】
請求項７乃至１０のいずれかに記載の装置において、
前記コードブックベクトル決定部（１２０）は、前記自己相関行列Ｒに対してコレスキー分解を行って前記自己相関行列Ｒを分解するよう構成され、かつ前記行列分解を行って前記コードブックベクトルを決定するための前記対角行列Ｄを決定するよう構成されている、装置。
【請求項１３】
請求項１乃至１２のいずれかに記載の装置において、
前記コードブックベクトル決定部（１２０）は、前記スピーチ信号のゼロインパルス応答に依存して前記コードブックベクトルを決定するよう構成されている、装置。
【請求項１４】
請求項１乃至１３のいずれかに記載の装置において、
前記装置は代数符号励振線形予測スピーチ符号化を用いて前記スピーチ信号を符号化する符号器であり、
前記コードブックベクトル決定部（１２０）は、代数コードブックのコードブックベクトルとして、前記自己相関行列Ｒに基づいて前記コードブックベクトルを決定するよう構成されている、装置。
【請求項１５】
スピーチ符号化アルゴリズムのコードブックベクトルを決定することによって、スピーチ信号を符号化する方法であって、
自己相関行列Ｒを決定するステップと、
前記自己相関行列Ｒに依存して前記コードブックベクトルを決定するステップとを含み、
前記自己相関行列Ｒを決定するステップは、ベクトルｒのベクトル係数を決定するステップを含み、
前記自己相関行列Ｒは複数の行と複数の列とを含み、前記ベクトルｒは前記自己相関行列Ｒの列の１つ又は行の１つを示し、
R(i , j) = r(|i− j|)
であり、ここでR(i , j)は前記自己相関行列Ｒの係数を示し、ｉは前記自己相関行列Ｒの複数の行の１つを示す第１インデックスであり、ｊは前記自己相関行列Ｒの複数の列の１つを示す第２インデックスである、方法。
【請求項１６】
請求項１に記載の装置によって符号化された符号化済みスピーチ信号を復号化して、復号化済みスピーチ信号を得る復号器（２２０）。
【請求項１７】
請求項１５に記載の方法によって符号化された符号化済みスピーチ信号を復号化して、復号化済みスピーチ信号を得る方法。
【請求項１８】
入力スピーチ信号を符号化して、符号化済みスピーチ信号を得る請求項１乃至１４のいずれかに記載の装置（２１０）と、
前記符号化済みスピーチ信号を復号化して、復号化済みスピーチ信号を得る請求項１６に記載の復号器（２２０）と、
を含むシステム。
【請求項１９】
請求項１５に記載の方法に従って入力スピーチ信号を符号化し、符号化済みスピーチ信号を得るステップと、
請求項１７に記載の方法に従って前記符号化済みスピーチ信号を復号化し、復号化済みスピーチ信号を得るステップと、を含む方法。
【請求項２０】
コンピュータ又は信号プロセッサ上で実行されたとき、請求項１５、１７又は１９に記載の方法を実行するためのコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、オーディオ信号符号化に関し、特に、自己相関(autocorrelation)ドメインにおけるＡＣＥＬＰを用いたスピーチ信号の符号化装置に関する。
【背景技術】
【０００２】
ＣＥＬＰ(Code Excited Linear Prediction:符号励振線形予測)によるスピーチ符号化において、スピーチ信号のスペクトル包絡（又は等価的に短時間の時間構造）は線形予測（ＬＰ）モデルによって表され、予測残余は長時間予測子（ＬＴＰ，適応型コードブックとしても知られている）によってモデル化され、残余信号はコードブック（固定コードブックとしても知られている）によって表現される。後者、つまり固定コードブックは、一般に代数コードブックとして適用され、ここではそのコードブックが代数公式又はアルゴリズムによって表現されるため、そのアルゴリズムを除いて全体のコードブックを記憶しておく必要がない。また同時に、速やかな探索アルゴリズムを可能にする。残余に対して代数コードブックを適用するＣＥＬＰコーデックは、代数符号励振線形予測（ＡＣＥＬＰ）コーデック（非特許文献１、２、３、４を参照）として知られている。
【０００３】
スピーチ符号化において、代数残余コードブックを用いることは、例えば非特許文献１７、１３、１８のように、主流のコーデックの中で選択された一手法と言える。ＡＣＥＬＰは、線形予測（ＬＰ）フィルタによるスペクトル包絡のモデリングと、長時間予測子（ＬＴＰ）による音声の基本周波数のモデリングと、代数コードブックによる予測残余のモデリングとに基づいている。ＬＴＰと代数コードブックパラメータは、ある知覚ドメインにおいて最小二乗アルゴリズムによって最適化され、その知覚ドメインはあるフィルタによって規定される。
【０００４】
ＡＣＥＬＰ型アルゴリズムの演算的に最も複雑な部分つまりボトルネックは、残余コードブックの最適化である。現在知られている唯一の最適なアルゴリズムは、各サブフレームについてサイズＮpのスペースの全数探索(exhaustive search)であろう。そこでは、全てのポイントで
2015532456.tif 000003
の複雑性の評価が求められる。代表的な値はＰ＝８パルスを持つサブフレーム長さＮ＝６４（すなわち、５ｍｓ）であるから、このことは毎秒当り１０20回以上の演算を意味している。明らかに、このことは実行可能な選択肢ではない。ハードウエアの条件によって設定される複雑性の制限内に止めるために、コードブックの最適化手法は、非最適な反復的アルゴリズムで作動する必要がある。最適化プロセスへの多数のそのようなアルゴリズム及び改善が、例えば非特許文献１７、１９、２０、２１、２２等によって過去に提案されてきた。
【０００５】
明らかに、ＡＣＥＬＰ最適化は、線形予測モデルの出力としてスピーチ信号ｘ（ｎ）を記述することに基づいており、その場合、推定されたスピーチ信号は、
2015532456.tif 000004
となり、ここで、ａ（ｋ）はＬＰ係数であり、
2015532456.tif 000005
は残余信号である。ベクトル形式では、この式は以下のように表すことができる。
2015532456.tif 000006
ここで、行列Ｈは対角要素ｈ（０）と下対角要素ｈ（１），．．．，ｈ（３９）とを持つ下三角のテプリッツ畳み込み行列(lower triangular Toeplitz convolution matrix)として定義され、ベクトルｈ（ｋ）はＬＰモデルのインパルス応答である。この表記法では、知覚モデル（通常は重み付きＬＰモデルと対応している）は省略されている点に注意すべきであるが、知覚モデルはインパルス応答ｈ（ｋ）内に含まれていると考えられる。この省略により、一般的な結果に影響を及ぼすことなく、表記を簡素化できる。この知覚モデルの包含は非特許文献１内で適用されている。
【０００６】
モデルの適合性は二乗誤差によって測定される。すなわち、
2015532456.tif 000007
この二乗誤差は最適なモデルパラメータを発見するために使用される。ここで、ＬＴＰとパルスコードブックは共にベクトルｅをモデル化するために使用されると考えられる。実際の適用は、対応する文献（非特許文献１〜４参照）の中で見ることができる。
【０００７】
実際のところ、上述の適合性の尺度は以下のように簡略化され得る。行列Ｂ＝ＨTＨがｈ（ｎ）の相関を含み、ｃkがｋ番目の固定コードブックベクトルであり、
2015532456.tif 000008
と設定する。ここで、ｇはゲインファクタである。ｇが最適に選択されると仮定すると、コードブックは次式の探索基準を最大化することによって、探索される。
2015532456.tif 000009
ここで、d = HTxは目標ベクトルとインパルス応答ｈ（ｎ）との間の相関を含むベクトルであり、上付き文字Ｔは転置を示す。ベクトルｄと行列Ｂはコードブック探索の前に計算される。この公式は、ＬＴＰとパルスコードブックの両方の最適化において共通に使用される。
【０００８】
上記公式の使用を最適化するにあたり、多くの研究がなされてきた。例えば、
１）探索アルゴリズムによって実際にアクセスされる行列Ｂの要素だけが計算される。又は、
２）パルス探索の試行錯誤アルゴリズムは、以前の予備選別（例えば非特許文献１、５参照）に基づいて、高い成功確率を持つコードブックベクトルだけを試すように、削減される。
【０００９】
ＡＣＥＬＰアルゴリズムの実際の詳細は、ゼロインパルス応答（ＺＩＲ）の概念に関連している。この概念は、合成された残余に対比してオリジナル領域の合成信号を考慮するときに現れる。残余はフレーム又はサブフレームサイズに対応したブロック内に符号化される。しかし、式１のＬＰモデルを用いてオリジナル領域信号を合成する際に、固定長の残余は、ＬＰフィルタのインパルス応答に対応した無限長の「テール」を持つであろう。すなわち、残余コードブックベクトルは有限長であるが、現在のフレーム又はサブフレームをはるかに超えて、合成信号に影響を及ぼすであろう。１つのフレームの将来に対する影響は、コードブックベクトルをゼロで拡張し、この拡張信号について式１の合成出力を計算することによって、計算され得る。合成された信号のこの拡張は、ゼロインパルス応答として知られている。次に、現在のフレームを符号化する際に以前のフレームの影響を考慮するために、以前のフレームのＺＩＲが現在のフレームの目標値から減算される。よって、現在のフレームを符号化する際、その信号の、以前のフレームによって既にモデル化されていなかった部分だけが考慮される。
【００１０】
実際上、ＺＩＲは次のようにして考慮される。（サブ）フレームＮ−１の符号化が完了したとき、量子化された残余は次の（サブ）フレームＮの長さまでゼロで拡張されている。拡張された量子化済み残余はＬＰによってフィルタリングされ、量子化済み信号のＺＩＲが取得される。量子化済み信号のＺＩＲは、次にオリジナル（量子化されていない）信号から減算され、この修正済み信号が（サブ）フレームＮを符号化する際に目標信号を形成する。このようにして、（サブ）フレームＮ−１内に形成された全ての量子化誤差は、（サブ）フレームＮを量子化する際に考慮されるであろう。このやり方は、出力信号の知覚的品質をかなり改善できる。
【００１１】
しかしながら、オーディオ符号化にとってさらに改善された概念が提供されるならば、一層望ましい。
【先行技術文献】
【非特許文献】
【００１２】
【非特許文献１】[1] Salami, R. and Laflamme, C. and Bessette, B. and Adoul, J.P., “ITU-T G. 729 Annex A: reduced complexity 8 kb/s CS-ACELP codec for digital simultaneous voice and data“, Communications Magazine, IEEE, vol 35, no 9, pp 56-63, 1997.
【非特許文献２】[2] 3GPP TS 26.190 V7.0.0 , “Adaptive Multi-Rate (AMR-WB) speech codec“, 2007.
【非特許文献３】[3] 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.
【非特許文献４】[4] Schroeder, M. and Atal, B., “Code-excited linear prediction (CELP): High-quality speech at very low bit rates“, Acoustics, Speech, and Signal Processing, IEEE Int Conf, pp 937-940, 1985.
【非特許文献５】[5] Byun, K.J. and Jung, H.B. and Hahn, M. and Kim, K.S., “A fast ACELP codebook search method“, Signal Processing, 2002 6th International Conference on, vol 1, pp 422-425, 2002.
【非特許文献６】[6] G. H. Golub and C. F. van Loan, “Matrix Computations“, 3rd Edition, John Hopkins University Press, 1996.
【非特許文献７】[7] Boley, D.L. and Luk, F.T. and Vandevoorde, D., “Vandermonde factorization of a Hankel matrix“, Scientific computing, pp 27-39, 1997.
【非特許文献８】[8] Baeckstroem, T. and Magi, C., “Properties of line spectrum pair polynomials - A review”, Signal processing, vol. 86, no. 11, pp. 3286-3298, 2006.
【非特許文献９】[9] A. Haermae, M. Karjalainen, L. Savioja, V. Vaelimaeki, U. Laine, and J. Huopaniemi, “Frequencywarped signal processing for audio applications,” J. Audio Eng. Soc, vol. 48, no. 11, pp. 1011-1031, 2000.
【非特許文献１０】[10] T. Laakso, V. Vaelimaeki, M. Karjalainen, and U. Laine, “Splitting the unit delay [FIR/all pass filters design],” IEEE Signal Process. Mag.,vol. 13, no. 1, pp. 30-60, 1996.
【非特許文献１１】[11] J. Smith III and J. Abel, “Bark and ERB bilinear transforms,” IEEE Trans. Speech Audio Process., vol. 7, no. 6, pp. 697-708, 1999.
【非特許文献１２】[12] R. Schappelle, “The inverse of the confluent Vandermonde matrix,” IEEE Trans. Autom. Control, vol. 17, no. 5, pp. 724-725, 1972.
【非特許文献１３】[13] B. Bessette, R. Salami, R. Lefebvre, M. Jelinek, J. Rotola-Pukkila, J. Vainio, H. Mikkola, and K. Jarvinen, “The adaptive multirate wideband speech codec (AMR-WB),” Speech and Audio Processing, IEEE Transactions on, vol. 10, no. 8, pp. 620-636, 2002.
【非特許文献１４】[14] M. Bosi and R. E. Goldberg, Introduction to Digital Audio Coding and Standards. Dordrecht, The Netherlands: Kluwer Academic Publishers, 2003.
【非特許文献１５】[15] B. Edler, S. Disch, S. Bayer, G. Fuchs, and R. Geiger, “A time-warped MDCT approach to speech transform coding,” in Proc 126th AES Convention, Munich, Germany, May 2009.
【非特許文献１６】[16] J. Makhoul, “Linear prediction: A tutorial review,” Proc. IEEE, vol. 63, no. 4, pp. 561-580, April 1975.
【非特許文献１７】[17] J.-P. Adoul, P. Mabilleau, M. Delprat, and S. Morissette, “Fast CELP coding based on algebraic codes,” in Acoustics, Speech, and Signal Processing, IEEE Int Conf (ICASSP’87), April 1987, pp. 1957-1960.
【非特許文献１８】[18] ISO/IEC 23003-3:2012, “MPEG-D (MPEG audio technologies), Part 3: Unified speech and audio coding,” 2012.
【非特許文献１９】[19] F.-K. Chen and J.-F. Yang, “Maximum-take-precedence ACELP: a low complexity search method,” in Acoustics, Speech, and Signal Processing, 2001. Proceedings.(ICASSP’01). 2001 IEEE International Conference on, vol. 2. IEEE, 2001, pp. 693-696.
【非特許文献２０】[20] R. P. Kumar, “High computational performance in code exited linear prediction speech model using faster codebook search techniques,” in Proceedings of the International Conference on Computing: Theory and Applications. IEEE Computer Society, 2007, pp. 458-462.
【非特許文献２１】[21] N. K. Ha, “A fast search method of algebraic codebook by reordering search sequence,” in Acoustics, Speech, and Signal Processing, 1999. Proceedings., 1999 IEEE International Conference on, vol. 1. IEEE, 1999, pp. 21-24.
【非特許文献２２】[22] M. A. Ramirez and M. Gerken, “Efficient algebraic multipulse search,” in Telecommunications Symposium, 1998. ITS’98 Proceedings. SBT/IEEE International. IEEE, 1998, pp. 231-236.
【非特許文献２３】[23] ITU-T Recommendation G.191, “Software tool library 2009 user’s manual,” 2009.
【非特許文献２４】[24] ITU-T Recommendation P.863, “Perceptual objective listening quality assessment,” 2011.
【非特許文献２５】[25] T. Thiede, W. Treurniet, R. Bitto, C. Schmidmer, T. Sporer, J. Beerends, C. Colomes, M. Keyhl, G. Stoll, K. Brandeburg et al., “PEAQ - the ITU standard for objective measurement of perceived audio quality,” Journal of the Audio Engineering Society, vol. 48, 2012.
【非特許文献２６】[26] ITU-R Recommendation BS.1534-1, “Method for the subjective assessment of intermediate quality level of coding systems,” 2003.
【発明の概要】
【発明が解決しようとする課題】
【００１３】
本発明の目的は、オーディオオブジェクト符号化についてそのような改善された概念を提供することである。
【課題を解決するための手段】
【００１４】
本発明の目的は、請求項１に係る装置と、請求項１５に係る方法と、請求項１６に係る復号器と、請求項１７に係る復号化方法と、請求項１８に係るシステムと、請求項１９に係る方法と、請求項２０に係るコンピュータプログラムとによって達成される。
【００１５】
スピーチ符号化アルゴリズムのコードブックベクトルを決定することによって、スピーチ信号を符号化する装置が提供される。この装置は、自己相関行列Ｒを決定する行列決定部と、自己相関行列Ｒに依存してコードブックベクトルを決定するコードブックベクトル決定部とを含む。行列決定部は、ベクトルｒのベクトル係数を決定することによって、自己相関行列Ｒを決定するよう構成され、自己相関行列Ｒは複数の行と複数の列とを含み、ベクトルｒは自己相関行列Ｒの列の１つ又は行の１つを示し、R(i , j) = r(|i− j|)であり、R(i , j)は自己相関行列Ｒの係数を示し、ｉは自己相関行列Ｒの複数の行の１つを示す第１インデックスであり、ｊは自己相関行列Ｒの複数の列の１つを示す第２インデックスである。
【００１６】
この装置は、スピーチ信号を符号化するためにコードブックベクトルを使用するよう構成されている。例えば、この装置は符号化済みスピーチ信号を生成することができ、この符号化済みスピーチ信号は、複数の線形予測係数と、有声音の基本周波数の指示（例えばピッチパラメータ）と、コードブックベクトルの指示、例えばコードブックベクトルのインデックスとを含む。
【００１７】
さらに、上述の実施形態に従う装置によって符号化された符号化済みスピーチ信号を復号化して復号化済みスピーチ信号を得る復号器が提供される。
【００１８】
さらに、あるシステムが提供される。このシステムは、上述の実施形態に従い入力スピーチ信号を符号化して符号化済みスピーチ信号を得る装置を含む。さらに、このシステムは、上述の実施形態に従い符号化済みスピーチ信号を復号化して復号化済みスピーチ信号を得る復号器を含む。
【００１９】
スピーチ符号化アルゴリズムＡＣＥＬＰの目的関数についての改善された概念が提供され、そこでは、現在のフレームのパラメータを最適化するに際し、以前のフレームのインパルス応答の現在のフレームに対する影響だけでなく、現在のフレームのインパルス応答の次のフレームに対する効果をも考慮に入れる。幾つかの実施形態が、相関行列を変更することによって、このような改善を達成でき、その相関行列は自己相関行列への従来のＡＣＥＬＰ最適化の中核にあり、エルミート・テプリッツ構造（Hermitian Toeplitz structure）を有する。この構造を用いることにより、演算上の複雑さとメモリ要求の観点からみてＡＣＥＬＰ最適化をさらに効率化することができる。同時に、適用される知覚モデルはさらに一貫性がよくなり、フレーム間依存性が回避され得るので、パケットロスの影響下での性能を改善できる。
【００２０】
ＡＣＥＬＰ方式を用いたスピーチ符号化は、知覚ドメインでの最小二乗アルゴリズムに基づいており、その知覚ドメインはフィルタによって特定される。本発明の実施形態によれば、最小二乗問題の従来の定義の演算上の複雑さは、ゼロインパルス応答の次のフレームへの影響を考慮することで軽減され得る。提案の修正法では、目的関数の中で現れる相関行列にテプリッツ構造を導入し、その目的関数とは構造を簡素化し演算量を削減するものである。提案の概念は、知覚的品質を低下させずに１７％まで演算上の複雑さを低減できる。
【００２１】
本発明の実施形態は、目的関数のわずかな修正によって、残余コードブックの最適化における演算量をさらに削減できるという知見に基づいている。この演算量の削減は、知覚的品質における低下を招かずに達成される。代替法として、ＡＣＥＬＰ残余最適化は反復的な探索アルゴリズムに基づいていることから、提案の修正法を用いて、演算量を増加させずに反復回数を増加させることが可能となり、この方法によって改善された知覚的品質を得ることも可能になる。
【００２２】
既存、及び修正された目的関数は、知覚をモデル化し、知覚歪みを最小化するよう努めるものである。しかし、従来の手法に対する最適解は修正された目的関数に関して必ずしも最適ではなく、逆もまた然りである。このことだけで、１つの手法が他の手法より優れていることを意味するものではないが、分析的な議論によると、修正された目的関数の方が一層一貫性を有することが示されている。具体的には、従来の目的関数とは対照的に、提案の概念は、１つのサブフレーム内の全てのサンプルを、一貫しかつ良好に定義された知覚的及び信号モデルを用いて、同等に取り扱う。
【００２３】
いくつかの実施形態において、提案の修正法は、残余コードブックの最適化だけを変更するように適用可能である。したがって、この方法はビットストリーム構造を変更せず、既存のＡＣＥＬＰコーデックに対する後方互換性をもって適用可能である。
【００２４】
さらに、スピーチ符号化アルゴリズムのコードブックベクトルを決定することにより、スピーチ信号を符号化する方法が提案される。この方法は、
−自己相関行列Ｒを決定するステップと、
−自己相関行列Ｒに依存してコードブックベクトルを決定するステップと、を含む。
【００２５】
自己相関行列Ｒを決定するステップは、ベクトルｒのベクトル係数を決定するステップを含む。自己相関行列Ｒは複数の行と複数の列とを含む。ベクトルｒは自己相関行列Ｒの前記列の１つ又は前記行の１つを示し、ここで
R(i , j) = r(|i− j|)
である。
【００２６】
Ｒ（ｉ，ｊ）は自己相関行列Ｒの係数を示し、ｉは自己相関行列Ｒの複数の行の１つを示す第１インデックスであり、ｊは自己相関行列Ｒの複数の列の１つを示す第２インデックスである。
【００２７】
さらに、上述の実施形態に従ってスピーチ信号を符号化する方法に従って符号化された符号化済みスピーチ信号を復号化し、復号化済みスピーチ信号を得る方法が提案される。
【００２８】
さらに、次の方法が提案される。この方法は、
−スピーチ信号を符号化する上述の方法に従って入力スピーチ信号を符号化し、符号化済みスピーチ信号を得るステップと、
−スピーチ信号を復号化する上述の方法に従って符号化済みスピーチ信号を復号化して復号化済みスピーチ信号を得るステップと、を含む。
【００２９】
さらに、コンピュータ又は信号プロセッサ上で実行されたとき、上述の方法を実行するためのコンピュータプログラムが提供される。
【００３０】
好適な実施形態は従属項において提供される。
【００３１】
以下に、本発明の好ましい実施形態について図面を参照しながら説明する。
【図面の簡単な説明】
【００３２】
【図１】一実施形態に従い、スピーチ符号化アルゴリズムのコードブックベクトルを決定することにより、スピーチ信号を符号化する装置を示す。
【図２】一実施形態に従う復号器を示す。
【図３】一実施形態に従ってスピーチ信号を符号化する装置と復号器とを含むシステムを示す。
【発明を実施するための形態】
【００３３】
図１は、一実施形態に従いスピーチ符号化アルゴリズムのコードブックベクトルを決定することにより、スピーチ信号を符号化する装置を示す。
【００３４】
この装置は、自己相関行列Ｒを決定するための行列決定部（１１０）と、自己相関行列Ｒに依存してコードブックベクトルを決定するためのコードブックベクトル決定部（１２０）とを備えている。
【００３５】
行列決定部（１１０）は、ベクトルｒのベクトル係数を決定することにより、自己相関行列Ｒを決定するよう構成されている。
【００３６】
自己相関行列Ｒは複数の行と複数の縦列とを含み、ベクトルｒは自己相関行列Ｒの列の１つと行の１つとを示し、R(i , j) = r(|i− j|)である。
【００３７】
Ｒ（ｉ，ｊ）は自己相関行列Ｒの係数を示し、ｉは自己相関行列Ｒの複数の行の１つを示す第１インデックスであり、ｊは自己相関行列Ｒの複数の列の１つを示す第２インデックスである。
【００３８】
この装置は、スピーチ信号を符号化するためにコードブックベクトルを使用するよう構成されている。例えば、この装置は符号化済みスピーチ信号を生成することができ、この符号化済みスピーチ信号は、複数の線形予測係数と、有声音の基本周波数の指示（例えばピッチパラメータ）と、コードブックベクトルの指示とを含む。
【００３９】
例えば、スピーチ信号を符号化する特定の実施形態によれば、本装置は、スピーチ信号に依存して複数の線形予測係数（ａ（ｋ））を決定するよう構成されてもよい。さらに、本装置は、複数の線形予測係数（ａ（ｋ））に依存して残余信号を決定するよう構成される。さらに加えて、行列決定部１１０は、残余信号に依存して自己相関行列Ｒを決定するよう構成されてもよい。
【００４０】
以下では、本願発明の幾つかのさらなる実施形態について説明する。
【００４１】
式３と式４とに戻って、式３は次のように知覚モデルの適合性を示す二乗誤差を定義しており、
2015532456.tif 000010
式４は最大化されるべき探索基準を示している。
2015532456.tif 000011
【００４２】
ＡＣＥＬＰアルゴリズムは、式４に中心を合わせたものであり、一方、式４は式３に基づいている。
【００４３】
幾つかの実施形態は、これら式の分析により、量子化された残余値ｅ（ｋ）が、インデックスｋに依存して誤差エネルギーε2に対して非常に異なる効果を持つことが明らかである、という知見に基づいている。例えば、インデックスｋ＝１及びｋ＝Ｎの場合を考慮した場合、もし残余コードブックの唯一の非ゼロ値がｋ＝１のときに現れたならば、次に誤差エネルギーε2は次式のようになり、
2015532456.tif 000012
一方、ｋ＝Ｎの場合には、誤差エネルギーε2は次式のようになる。
2015532456.tif 000013
【００４４】
換言すれば、ｅ（１）は範囲１〜Ｎのインパルス応答ｈ（ｋ）で重み付けられ、一方ｅ（Ｎ）はｈ（１）だけで重み付けられる。スペクトル重みの点からみれば、このことは、極端な場合にはｅ（Ｎ）が線形的に重み付けられるように、各ｅ（ｋ）が異なるスペクトル重み関数で重み付けられることを意味している。知覚的モデリングの観点からみれば、１つのフレームの中の全てのサンプルに対して同じ知覚的重みを適用することは、有意義であろう。よって、式３はＺＩＲを次のフレームの中へと考慮するように拡張されるべきである。ここで、特に、先行技術との違いは、前のフレームからのＺＩＲと次のフレーム内へのＺＩＲとの両方が考慮されるという点に注目すべきである。
【００４５】
ｅ（ｋ）をオリジナルの量子化されていない残余とし、
2015532456.tif 000014
を量子化された残余と仮定する。さらに、両方の残余を１〜Ｎの範囲内で非ゼロとし、その他の範囲ではゼロと仮定する。そこで、次式となる。
2015532456.tif 000015
【００４６】
同様に、同じ関係式を行列形式で次のように表すことができる。
2015532456.tif 000016
ここで、
2015532456.tif 000017
はインパルス応答ｈ（ｋ）に対応した無限次元の畳み込み行列である。これを式３に代入すると、
2015532456.tif 000018
となり、ここで
2015532456.tif 000019
はｈ（ｎ）の自己相関に対応した有限サイズのエルミート・テプリッツ行列である。式４についても同様の導出により、次の目的関数が得られる。
2015532456.tif 000020
【００４７】
この目的関数は、式４と非常に似ている。主たる相違は、相関行列Ｂに代えてエルミート・テプリッツ行列Ｒを分母に用いた点である。
【００４８】
上述の通り、この新規な公式化は、１つのフレーム内の残余ｅの全てのサンプルが同じ知覚的重み付けを受けることになる、という利点を有する。しかしながら、重要なことは、この公式化は演算上の複雑さと共にメモリ要求についても、かなりの利点をもたらすことである。Ｒはエルミート・テプリッツ行列であるから、第１列ｒ（０）．．．ｒ（Ｎ−１）はこの行列を完全に定義している。換言すれば、完全なＮ×Ｎの行列を記憶する代わりに、Ｎ×１のベクトルｒ（ｋ）を記憶するだけで十分であるから、メモリ割り当てにおいてかなりの節約になるということである。さらに、全てのＮ×Ｎ要素を決定する必要がないだけでなく、第１のＮ×１列を決定するだけでよいので、演算上の複雑さも同様に低減される。また、要素（ｉ，ｊ）がR(i , j) = r(|i− j|)によって発見できるので、その行列内での指標付けも簡易になる。
【００４９】
式１０における目的関数が式４と非常に類似しているので、一般のＡＣＥＬＰの構造は維持され得る。特に、以下の如何なる演算もいずれかの目的関数を用いて、アルゴリズムのわずかな修正だけで実行できる。
１．ＬＴＰ遅れ（適応型コードブック）の最適化
２．残余のモデル化のためのパルスコードブック（固定されたコードブック）の最適化
３．個別的又は連帯的な、ＬＴＰとパルスのゲインの最適化
４．性能が式３の二乗誤差によって測定される他のパラメータの最適化
【００５０】
従来のＡＣＥＬＰアプリケーションにおいて修正されるべき唯一の部分は、行列Ｒで置き換えられる相関行列Ｂと、次のフレームへＺIRを包含させなければならない目標値との取り扱いである。
【００５１】
幾つかの実施例では、相関行列ＢがＡＣＥＬＰアルゴリズム内のどの箇所に表れるかに拘わらず、相関行列Ｂを自己相関行列Ｒによって置き換えることにより、本発明の概念を利用する。行列Ｂの全ての出現が無くなれば、その値を計算することが回避され得る。
【００５２】
例えば、自己相関行列Ｒは、自己相関行列Ｒの第１列r(0), .., r(N-1)の係数を決定することによって、決定される。
【００５３】
行列Ｒは式９においてＲ＝ＨTＨによって定義される。従って、その要素Ｒij＝ｒ（ｉ−ｊ）は次式によって計算され得る。
2015532456.tif 000021
【００５４】
すなわち、系列r(k)はh(k)の自己相関である。
【００５５】
しかしながら、ｒ（ｋ）は、さらに効率的な手段によって取得され得る場合も多い。特に、ＡＭＲやＧ．７１８などのスピーチ符号化標準では、系列ｈ（ｋ）は、プリエンファシスを含めるべく用いられる知覚的重み関数Ｗ（ｚ）によってフィルタリングされた線形予測フィルタＡ（ｚ）のインパルス応答である。換言すれば、ｈ（ｋ）は線形予測モデルの知覚的に重み付けられたインパルス応答を示している。
【００５６】
フィルタＡ（ｚ）は、通常、スピーチ信号ｒX（ｋ）の自己相関から推定され、ここで、ｒX（ｋ）は既知である。Ｈ（ｚ）＝Ａ-1（ｕ）Ｗ（ｚ）であるから、自己相関系列ｒ（ｋ）は次式によりｗ（ｋ）の自己相関を計算することによって決定され得る。
2015532456.tif 000022
そのため、ｈ（ｋ）の自己相関は次式のようになる。
2015532456.tif 000023
【００５７】
幾つかの実施例では、これら式は全体システムの設計に依存して修正されてもよい。
【００５８】
次に、コードブックのコードブックベクトルは、例えば自己相関行列Ｒに基づいて決定されてもよい。特に、幾つかの実施例に従えば、式１０はコードブックのコードブックベクトルを決定するために使用されてもよい。
【００５９】
この文脈において、式１０は、
2015532456.tif 000024
の形式で目的関数を定義しており、ここでは行列Ｒが対称的なテプリッツ構造を有するという点を除けば、この形式はスピーチ符号化標準ＡＭＲやＧ．７１８におけるのと同様な形式である。目的関数は、基本的に、目標ベクトルｄとコードブックベクトル
2015532456.tif 000025
との間の正規化された相関であり、最も良好な可能性のあるコードブックベクトルは、正規化された相関
2015532456.tif 000026
について最も高い値を与える、つまり正規化された相関
2015532456.tif 000027
を最大化するものである。
【００６０】
よって、コードブックベクトルは上述の標準における幾つかの手法と同様の手法を用いて最適化され得る。特に、以下に説明するように、例えば残余に対して最良の代数コードブック（即ち固定のコードブック）ベクトル
2015532456.tif 000028
を発見する非常に簡素なアルゴリズムが適用され得る。しかしながら、効率的な探索アルゴリズムの設計にかなりの努力が投下されてきたこと（ＡＭＲ及びＧ．７１８を参照）、及びこの探索アルゴリズムが唯一の例示し得る適用例であることに注目すべきである。
【００６１】
１．初期コードブックベクトル
2015532456.tif 000029
を定義し、パルス数をｐ＝０に設定する。
２．初期コードブックベクトル品質尺度をｆ0＝０に設定する。
３．一時的コードブック品質尺度を
2015532456.tif 000030
に設定する。
４．コードブックベクトルの各位置ｋについて、
（i）ｐを１ずつ増加する
（ii）もし位置ｋが負のパルスを既に含んでいる場合には、ステップ（vii）を続ける
（iii）一時的コードブックベクトル
2015532456.tif 000031
を生成して、位置ｋにおいて正のパルスを追加する
（iv）その一時的コードブックベクトルの品質をｆ（εp+）によって評価する
（v）その一時的コードブックベクトルがそれ以前の如何なるコードブックよりも良好である場合
2015532456.tif 000032
には、そのコードブックベクトルを保存し、
2015532456.tif 000033
に設定し、次の反復へと続く
（vi）もし位置ｋが正のパルスを既に含んでいる場合には、次の反復へと続く
（vii）一時的コードブックベクトル
2015532456.tif 000034
を生成し、位置ｋにおいて負のパルスを追加する
（viii）その一時的コードブックベクトルの品質をｆ（εp-）によって評価する
（ix）その一時的コードブックベクトルがそれ以前の如何なるコードブックよりも良好である場合
2015532456.tif 000035
には、このコードブックベクトルを保存し、
2015532456.tif 000036
に設定し、次の反復へと続く
５．コードブックベクトル
2015532456.tif 000037
を、保存されたコードブックベクトルの中で最後（つまり最良）であると定義する。
６．もしパルスｐの数が所望のパルスの数に到達すれば、次に出力ベクトルを
2015532456.tif 000038
と定義し、停止する。もしそうでなければ、ステップ４を継続する。
【００６２】
既に指摘したように、従来のＡＣＥＬＰアプリケーションと比べて、幾つかの実施例では、ＺＩＲを次のフレームに含めるように目標値が修正される。
【００６３】
式１は、ＡＣＥＬＰ型コーデックにおいて使用される線形予測モデルを表している。現在のフレーム（及び全ての将来のフレーム）の残余がゼロに設定された場合に、ゼロインパルス応答（ＺＩＲ、またゼロ入力応答としても知られている）は線形予測モデルの出力に当てはまる。このＺＩＲは、位置Ｎから前方（将来）へゼロである残余を次式のように定義することによって、容易に計算され得る。
2015532456.tif 000039
これにより、ＺＩＲは次のように定義され得る。
2015532456.tif 000040
【００６４】
このＺＩＲを入力信号から差し引くことにより、現在のフレームより前方へのみからの残余に依存する信号が得られる。
【００６５】
等価的に、ＺＩＲは過去の入力信号を次式のようにフィルタリングすることにより決定し得る。
2015532456.tif 000041
【００６６】
ＺＩＲが取り除かれている入力信号は、目標値としてよく知られており、位置Ｋで始まるフレームについてｄ（ｎ）＝ｘ（ｎ）−ＺＩＲK（ｎ）として定義され得る。この目標値は、原則としてＡＭＲ及びＧ．７１８標準における目標値と完全に等しい。信号を量子化する際、量子化された信号
2015532456.tif 000042
はフレームの持続時間Ｋ≦ｎ＜Ｋ＋Ｎについてｄ（ｎ）と比較される。
【００６７】
逆に、現在のフレームの残余は後続のフレームに対して影響を及ぼすので、信号を量子化する際には、その影響を考慮することが有用である。つまり、現在のフレームを超えてｎ＞Ｋ＋Ｎにおいても偏差
2015532456.tif 000043
を評価したいと望む可能性がある。しかしながら、それを実行するために、後続のフレームの残余をゼロに設定することによってのみ、現在のフレームの残余の影響を考慮したいと望む可能性もある。そのため、次のフレーム内への
2015532456.tif 000044
のＺＩＲが比較されてもよい。換言すれば、修正された目標値が次式で得られる。
2015532456.tif 000045
【００６８】
等価的に、Ａ（ｚ）のインパルス応答ｈ（ｎ）を使用することで、次式となる。
2015532456.tif 000046
【００６９】
この式は、簡便な行列形式ｄ'＝Ｈｅで記載することができ、ここでＨ及びｅは式２と同様に定義される。修正された目標値は、式２の正にｘであると見なすことができる。
【００７０】
行列Ｒの計算において、理論上、インパルス応答ｈ（ｋ）は無限数列であり、現実のシステムでは実現性のないものであることに留意されたい。
【００７１】
しかしながら、次のいずれかが可能である。
１）インパルス応答を有限長まで切り詰め又は窓掛けし、切り詰められたインパルス応答の自己相関を決定すること。
２）関連するＬＰ及び知覚フィルタのフーリエスペクトルを使用してインパルス応答のパワースペクトルを計算し、逆フーリエ変換によって自己相関を得ること。
【００７２】
ここで、ＬＴＰを用いた拡張方法について説明する。
【００７３】
長時間予測子（ＬＴＰ：long-time predictor）は、また事実上、線形予測子である。
【００７４】
ある実施形態によれば、行列決定部１１０は知覚的に重み付けされた線形予測子、例えば長時間予測子に基づいて、自己相関行列Ｒを決定するよう構成されてもよい。
【００７５】
ＬＰ及びＬＴＰは１つの結合予測子へと畳み込まれることができ、その結合予測子はスペクトル包絡形状とハーモニック構造との両方を含んでいる。そのような予測子のインパルス応答は非常に長くなるであろうし、そのため従来技術では取り扱うことが非常に難しいものである。しかしながら、線形予測子の自己相関が既知であれば、結合予測子の自己相関は、ＬＴＰを用いて前方及び後方に単に自己相関をフィルタリングすること、又は周波数領域での同様なプロセスによって、計算することができる。
【００７６】
ＬＴＰを用いる従来方法では、ＬＴＰ遅れがフレーム長さより短い場合に、問題があることに注意すべきである。なぜなら、ＬＴＰがそのフレーム内でフィードバック・ループを引き起こす恐れがあるからである。前記目的関数内にＬＴＰが含まれる長所は、ＬＴＰの遅れがフレーム長さよりも短い場合に、このフィードバックが最適化において明確に考慮される点である。
【００７７】
以下では、非相関領域における高速最適化のための拡張について説明する。
【００７８】
ＡＣＥＬＰシステムの設計における中心的な課題は、演算上の複雑性を削減することである。ＡＣＥＬＰシステムは、ＬＰによるフィルタリングが残余サンプル同士間の複雑な相関を引き起こすために複雑であり、それら残余サンプルとは行列Ｂ又は現在の文脈においては行列Ｒによって記述されるものである。ｅ（ｎ）のサンプル同士が相関しているので、ｅ（ｎ）を所望の精度で正確に量子化することは不可能である。しかし、式３又は式１０の目的関数に関してそれぞれ最適な量子化を発見するために、試行錯誤の手法を用いて異なる量子化の多くの組合せが試されなければならない。
【００７９】
行列Ｒの導入によって、これら相関に対する新たな視点が得られる。すなわち、Ｒはエルミート・テプリッツ構造を有するので、ハンケル行列の特異値分解、コレスキー分解（Cholesky decomposition）、ファンデルモンド分解（Vandermonde decomposition）などの複数の効率的な行列分解が適用され得る（非特許文献６及び７参照）（ハンケル行列はテプリッツ行列を逆転したものであり、そのため同じ分解がテプリッツ及びハンケル行列に適用され得る）。Ｒ＝ＥＤＥHをＲの分解とし、ＤがＲと同じサイズ、同じランクの対角行列とする。その場合、式９は次に以下のように修正され得る。
2015532456.tif 000047
ここで、
2015532456.tif 000048
である。Ｄは対角であるから、ｆ（ｋ）の各サンプルの誤差は他のサンプルｆ（ｉ）に対して独立している。式１０において、コードブックベクトルは最適ゲインによってスケールされ、それにより新たな目的関数は次のようになると考えられる。
2015532456.tif 000049
【００８０】
ここでも、（１つのラインの量子化を変更することは、全てのラインについて最適なゲインを変化させるので）サンプル同士は相関しているが、しかし、式１０と比較すれば、相関の影響はここでは限定されている。しかしながら、たとえ相関が考慮されても、この目的関数の最適化は式３又は式１０の最適化に比べてずっと簡単である。
【００８１】
この分解手法を使用して次のことが可能になる。
１．所望の精度で任意の通常のスカラー又はベクトル量子化技術を適用すること。
２．任意の従来のＡＣＥＬＰパルス探索アルゴリズムを用いて、目的関数として式１２を使用すること。
【００８２】
両方の手法が、式１２について最適に近い量子化を与える。従来の量子化技術が一般にいかなる強引な方法（可能なレートループを除いて）を必要とせず、かつ行列ＤがＢとＲのいずれよりも簡素であるから、両方の量子化方法は従来のＡＣＥＬＰパルス探索アルゴリズムに比べてより複雑性が低い。そのため、この手法における演算上の複雑性の主な原因は、行列分解の演算となる。
【００８３】
幾つかの実施例では、コードブックのコードブックベクトルを決定するために式１２を用いる。
【００８４】
例えば、Ｒ＝ＥHＤＥの形式のＲについて、幾つかの行列因数分解が存在する。例えば、
（ａ）例えばGNU Scientific Library (http://www.gnu.org/software/gsl/manual/html_node/Real-Symmetric-Matrices.html)を用いることで、固有値分解が計算され得る。行列Ｒは実数でかつ対称形（テプリッツと同様に）であり、それにより関数“gsl_eigen_symm()”が行列ＥとＤを決定するために使用され得る。同じ固有値分解の他の構成は非特許文献６において容易に取得できる。
（ｂ）非特許文献８に記載されたアルゴリズムを用いて、テプリッツ行列のファンデルモンド因数分解（非特許文献７）が使用され得る。このアルゴリズムは、行列ＥとＤとをもたらすものであり、Ｅはファンデルモンド行列であり、そのファンデルモンド行列は不均一な周波数分布を持つ離散フーリエ変換と等価である。
【００８５】
このような因数分解を使用することで、残余ベクトルｅはｆ＝ＥHｅ又はｆ'＝Ｄ1/2ＥHｅによって変換領域へ変換されることができる。この領域においてはいかなる一般的な量子化も適用可能であり、例えば
１．ベクトルｆ’はＡＣＥＬＰの一般的な構成とまさに同様に代数コードブックによって量子化され得る。しかし、ｆ’の成分同士は非相関であるから、ＡＣＥＬＰにおけるような複雑な探索関数は不要であり、以下のような単純なアルゴリズムが適用され得る。
（ａ）初期ゲインをｇ＝１に設定する。
（ｂ）f’を
2015532456.tif 000050
によって量子化する。
（ｃ）f’におけるパルスの数が所定の値ｐより大きく、
2015532456.tif 000051
の場合には、ゲインｇを増加させ、ステップｂに戻る。
（ｄ）その他の場合、すなわち
2015532456.tif 000052
におけるパルスの数が所定の値ｐより小さく、
2015532456.tif 000053
の場合には、ゲインｇを減少させ、ステップｂに戻る。
（ｅ）その他の場合、すなわち
2015532456.tif 000054
におけるパルスの数が所定の値ｐと等しく、
2015532456.tif 000055
の場合には、処理を停止させ得る。
２．標準ＡＭＲ−ＷＢ＋又はＭＰＥＧ  ＵＳＡＣのＴＣＸにおいてスペクトルラインの量子化に使用されているのと同様に、算術符号器が使用され得る。
【００８６】
ｆ’の要素が対角（式１２から分かるように）であり、それら要素が式１２の目的関数において同じ重みを持つので、それら要素は個別にかつ同じ量子化ステップサイズで量子化されることができる、という点に注意すべきである。その量子化は、式１２における目的関数の最適な（最大の）値を自動的に発見するであろうし、その量子化精度をもつことも可能である。換言すると、上述の量子化アルゴリズムは、共に式１２に関して最適な量子化をもたらすことになるであろう。
【００８７】
最適化のこの利点は、ｆ’の要素が別々に処理され得るという事実につながっている。もしあるコードブック手法が使用され、そこではコードブックベクトルｃkが非自明（１つ以上の非ゼロ要素を持つ）であるならば、これらコードブックベクトルはもはや独立した要素を持たず、行列因数分解の利点が失われるであろう。
【００８８】
テプリッツ行列のファンデルモンド因数分解が選択されることができ、そのファンデルモンド行列がフーリエ変換行列ではあるが、偏在する周波数を持つ点に注目されたい。換言すると、ファンデルモンド行列は周波数ワープされたフーリエ変換に対応する。つまりこの場合には、ベクトルｆはワープされた周波数スケール（非特許文献８の“root-exchange property”を参照）上での残余信号の周波数領域表現に対応する。
【００８９】
重要な点は、この結果は周知ではないということである。実際上、この結果は、もし信号ｘが畳み込み行列Ｃを用いてフィルタリングされた場合に、次式となることを示し、
2015532456.tif 000056
ここで、Ｖは（例えばワープされた）フーリエ変換（これは単位円上の要素を持つファンデルモンド行列である）であり、Ｄは対角行列である。すなわち、フィルタリングされた信号のエネルギーを測定することが必要ならば、周波数ワープされた信号のエネルギーが等価的に測定され得ることになる。逆に、ワープされたフーリエ領域においてなされるべき如何なる評価も、フィルタリングされた時間領域でも等価的に実行され得る。時間と周波数の二重性により、時間領域窓掛けと時間ワーピングとの間の等価性もまた存在する。しかし、実際上の問題は、上述の関係を満足する畳み込み行列Ｃを発見することが数値的に敏感な問題であり、それにより、畳み込み行列Ｃに代えて近似解
2015532456.tif 000057
を発見する方が容易であるという点である。
【００９０】
関係式||Ｃｘ||2＝||ＤＶｘ||2は、コードブックのコードブックベクトルを決定するために用いられ得る。
【００９１】
そのため、ここでＨによって、式２における畳み込み行列と同様の行列がＣに代えて示されるであろうということにまず注目すべきである。次に、もし量子化ノイズ
2015532456.tif 000058
を最小化したいと望む場合には、そのエネルギーは次式で求められる。
2015532456.tif 000059
【００９２】
次に、フレーム独立性についての拡張を説明する。
【００９３】
符号化されたスピーチ信号が、ラジオ電波のような不完全な伝送ラインを用いて不変に伝送された場合、データのパケットは時折消失するであろう。もしパケットＮ−１を完全に復号するためにはパケットＮが必要となるように、フレーム同士が互いに依存性を有するならば、パケットＮ−１の損失はパケットＮ−１とＮの両方の合成を損なうであろう。一方、もしフレーム同士が独立しているならば、パケットＮ−１の損失はパケットＮ−１の合成だけを損なうであろう。そのため、フレーム間依存性の影響を受けない装置および方法は重要である。
【００９４】
従来のＡＣＥＬＰシステムでは、フレーム間依存性の主たる原因はＬＴＰであり、ある程度までＬＰにも関係している。特に、両者は無限インパルス応答（ＩＩＲ）フィルタであるから、損傷されたフレームは損傷サンプルの「無限の」テールをもたらすであろう。実際のところ、そのようなテールは幾つかのフレーム長を持ち得るので、知覚的に煩わしいものである。
【００９５】
本発明の枠組みを使用することで、フレーム間依存性が生成される経路が、現在のフレームから次のフレームへとＺＩＲによって定量化されうる。このフレーム間依存性を回避するため、従来のＡＣＥＬＰに対する３つの修正法を実行する必要がある。
１．前のフレームから現在の（サブ）フレームまでのＺＩＲを計算する際、量子化された残余ではなく、ゼロで拡張されたオリジナル（量子化されていない）残余から、ＺＩＲが計算されるべきである。このようにして、前の（サブ）フレームからの量子化誤差は現在の（サブ）フレームに対して伝播しない。
２．現在のフレームを量子化する際、オリジナル信号と量子化された信号との間の、次のフレームへのＺＩＲにおける誤差が考慮されなければならない。このことは、上述のように相関行列Ｂを自己相関行列Ｒによって置き換えることによってなし得る。これにより、次のフレームへのＺＩＲにおける誤差が現在のフレーム内の誤差と共に最小化されることを確実にする。
３．誤差伝播がＬＰとＬＴＰとに起因しているので、両成分はＺＩＲに含まれていなければならない。これは、ＺＩＲがＬＰだけのために計算される従来の手法とは異なる点である。
【００９６】
もし、現在のフレームを量子化する際に前のフレームの量子化誤差が考慮されなければ、出力の知覚的品質における効率が低下する。したがって、誤差伝播の恐れがないときには、前の誤差を考慮するよう選択することが可能である。例えば、従来のＡＣＥＬＰシステムは、全ての２０ｍｓフレームが４個又は５個のサブフレームにサブ分割されるフレーミングを適用する。ＬＴＰと残余は、各サブフレームについて個別に量子化され、符号化されるが、全体のフレームはデータの１つのブロックとして伝送される。したがって、個別のサブフレームは消失されることはないが、完全なフレームが消失することはあり得る。さらに、フレーム境界においてのみフレーム独立性のＺＩＲを使用する必要があるが、残りのサブフレーム間ではＺＩＲがフレーム間依存性を持って使用され得る、ということになる。
【００９７】
幾つかの実施形態では、現在のフレームのインパルス応答の効果を次のフレーム内、つまり現在のフレームの目的関数内に含めることによって、従来のＡＣＥＬＰアルゴリズムを修正する。最適化問題の目的関数において、この修正は相関行列を、エルミート・テプリッツ構造を有する自己相関行列に置き換えることに相当する。この修正は次のような利点がある。
１．自己相関行列の追加されたエルミート・テプリッツ構造に起因して、演算の複雑性とメモリ要求とが低減される。
２．同じ知覚モデルが全てのサンプルに適用されるであろう。よって、知覚モデルの設計と調整がより簡易となり、その適用がより効率的にかつ一貫性を有することになる。
３．前のフレームからの量子化されていないインパルス応答と、次のフレームへの量子化されたインパルス応答とだけを考慮することにより、現在のフレームの量子化におけるフレーム間相関を完全に回避できる。このことは、パケットロスが予想されるようなシステムのロバスト性を改善できる。
【００９８】
図２は、上記実施例に従う装置によって符号化された符号化済みスピーチ信号を復号化して、復号化済みスピーチ信号を得るための復号器２２０を示す。この復号器２２０は、符号化済みスピーチ信号を受信するよう構成され、この符号化済みスピーチ信号は、上記実施例の１つに従いスピーチ信号を符号化する装置によって決定された、コードブックベクトルの指示、例えば決定されたコードブックベクトルのインテックスを含む。さらに、復号器２２０は符号化済みスピーチ信号を復号化して、前記コードブックベクトルに依存して復号化済みスピーチ信号を得るよう構成されている。
【００９９】
図３は、一実施例に従うシステムを示す。このシステムは、入力スピーチ信号を符号化して符号化済みスピーチ信号を得る上記実施例の１つに従う装置２１０を含む。符号化済みスピーチ信号は、スピーチ信号を符号化する装置２１０によって決定された決定済みコードブックベクトルの指示、例えばコードブックベクトルのインデックスを含む。さらに、このシステムは、符号化済みスピーチ信号を復号化して、復号化済みスピーチ信号を得る上記実施例に従う復号器２２０を含む。復号器２２０は符号化済みスピーチ信号を受信するよう構成されている。さらに、復号器２２０は符号化済み信号を復号化して、前記決定済みコードブックベクトルに依存して復号化済み信号を得るよう構成されている。
【０１００】
これまで装置の文脈で幾つかの態様を説明してきたが、これらの態様は対応する方法の説明でもあることは明らかであり、そのブロック又は装置が方法ステップ又は方法ステップの特徴に対応することは明らかである。同様に、方法ステップを説明する文脈で示した態様もまた、対応する装置の対応するブロック、項目又は特徴を表している。
【０１０１】
本発明の分解された信号は、デジタル記憶媒体に格納されてもよく、またインターネットのような無線伝送媒体や有線伝送媒体のような伝送媒体によって伝送されてもよい。
【０１０２】
構成要件にも依るが、本発明の実施形態は、ハードウエア又はソフトウエアにおいて構成可能である。この構成は、その中に格納される電子的に読み取り可能な制御信号を有し、本発明の各方法が実行されるようにプログラム可能なコンピュータシステムと協働する（又は協働可能な）、デジタル記憶媒体、例えばフレキシブルディスク，ＤＶＤ，ＣＤ，ＲＯＭ，ＰＲＯＭ，ＥＰＲＯＭ，ＥＥＰＲＯＭ又はフラッシュメモリなどを使用して実行することができる。
【０１０３】
本発明に従う幾つかの実施形態は、上述した方法の１つを実行するようプログラム可能なコンピュータシステムと協働可能で、電子的に読み取り可能な制御信号を有する非一時的データキャリアを含んでも良い。
【０１０４】
一般的に、本発明の実施例は、プログラムコードを有するコンピュータプログラム製品として構成することができ、このプログラムコードは当該コンピュータプログラム製品がコンピュータ上で作動するときに、本発明の方法の一つを実行するよう作動する。そのプログラムコードは例えば機械読み取り可能なキャリアに記憶されても良い。
【０１０５】
本発明の他の実施形態は、上述した方法の１つを実行するための、機械読み取り可能なキャリアに格納されたコンピュータプログラムを含む。
【０１０６】
換言すれば、本発明の方法の一実施形態は、そのコンピュータプログラムがコンピュータ上で作動するときに、上述した方法の１つを実行するためのプログラムコードを有する、コンピュータプログラムである。
【０１０７】
本発明の他の実施形態は、上述した方法の１つを実行するために記録されたコンピュータプログラムを含む、データキャリア（又はデジタル記憶媒体又はコンピュータ読み取り可能な媒体）である。
【０１０８】
本発明の他の実施形態は、上述した方法の１つを実行するためのコンピュータプログラムを表現するデータストリーム又は信号列である。そのデータストリーム又は信号列は、例えばインターネットを介するデータ通信接続を介して伝送されるように構成されても良い。
【０１０９】
他の実施形態は、上述した方法の１つを実行するように構成又は適用された、例えばコンピュータ又はプログラム可能な論理デバイスのような処理手段を含む。
【０１１０】
他の実施形態は、上述した方法の１つを実行するためのコンピュータプログラムがインストールされたコンピュータを含む。
【０１１１】
幾つかの実施形態において、（例えば書換え可能ゲートアレイのような）プログラム可能な論理デバイスは、上述した方法の幾つか又は全ての機能を実行するために使用されても良い。幾つかの実施形態において、書換え可能ゲートアレイは、上述した方法の１つを実行するためにマイクロプロセッサと協働しても良い。一般的に、そのような方法は、好適には任意のハードウエア装置によって実行される。
【０１１２】
上述した実施形態は、本発明の原理を単に例示的に示したにすぎない。本明細書に記載した構成及び詳細について修正及び変更が可能であることは、当業者にとって明らかである。従って、本発明は、本明細書に実施形態の説明及び解説の目的で提示した具体的詳細によって限定されるものではなく、添付した特許請求の範囲によってのみ限定されるべきである。
【図１】
2015532456.tif 000069
【図２】
2015532456.tif 000070
【図３】
2015532456.tif 000071
【国際調査報告】
2015532456.tif 000072
2015532456.tif 000073