JPH0844399A - 音響信号変換符号化方法および復号化方法 - Google Patents

音響信号変換符号化方法および復号化方法

Info

Publication number
JPH0844399A
JPH0844399A JP7052389A JP5238995A JPH0844399A JP H0844399 A JPH0844399 A JP H0844399A JP 7052389 A JP7052389 A JP 7052389A JP 5238995 A JP5238995 A JP 5238995A JP H0844399 A JPH0844399 A JP H0844399A
Authority
JP
Japan
Prior art keywords
coefficient
residual
outline
spectrum
acoustic signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP7052389A
Other languages
English (en)
Other versions
JP3087814B2 (ja
Inventor
Naoki Iwagami
直樹 岩上
Takehiro Moriya
健弘 守谷
Satoshi Miki
聡 三樹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP07052389A priority Critical patent/JP3087814B2/ja
Publication of JPH0844399A publication Critical patent/JPH0844399A/ja
Application granted granted Critical
Publication of JP3087814B2 publication Critical patent/JP3087814B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Transmission Systems Not Characterized By The Medium Used For Transmission (AREA)

Abstract

(57)【要約】 【目的】 音響信号の周波数特性を能率的に平坦化し、
少ない情報量で符号化する。 【構成】 入力音響信号を変形離散コサイン変換してそ
のスペクトラム特性を求める。他方、入力音響信号から
線形予測分析部１７で予測係数を求め、その予測係数を
スペクトラム概形計算部２１でフーリエ変換して音響信
号のスペクトラム特性の概形を求め、正規化部２２でそ
のスペクトラム特性をその概形で正規化して残差係数を
得る。残差概形正規化部２６はその残差係数を、残差概
形計算部２３で予測した残差概形により正規化して微細
構造係数を得て、量子化部２５はその微細構造係数をベ
クトル量子化する。逆正規化部３１は量子化微細構造を
逆正規化して残差係数を再生し、残差概形計算部２３は
その再生残差係数から残次フレームの残差係数の概形を
予測する。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】この発明はオーディオ信号、特に
音楽信号や音声信号などの音響信号を、周波数領域係数
に変換し、その周波数領域係数をできるだけ少ない情報
量でディジタル符号化する符号化方法、及びその符号化
音響信号を復号化する復号化方法に関する。
【０００２】
【従来の技術】現在、オーディオ信号を高能率に符号化
する方法として、原音をフレームと呼ばれる５〜５０ｍ
ｓ程度の一定間隔の区間に分割し、その１フレームの信
号に時間−周波数変換（フーリエ変換）を行って得た周
波数領域係数（周波数軸上のそれぞれの点におけるサン
プル値）を、その周波数特性の包絡形状（スペクトラム
概形）と、周波数領域係数をスペクトラム概形で平坦化
して得られる残差係数との２つの情報に分離し、それぞ
れを符号化することが提案されている。このような符号
化法として、適応スペクトラム聴感制御エントロピー符
号化法（ASPEC ：Adaptive Spectral Perceptual Entro
py Coding ）や、重み付きベクトル量子化による変換符
号化法（TCWVQ：Transform Coding with Weighted Vect
or Quantization)、エムペグ−オーディオ・レイヤ３方
式（MPEG-Audio Layer III）などがある。それぞれの技
術については、K.Brandenburg, J.Herre, J.D.Johnston
etal：“ASPEC：Adaptive spectral entropy coding o
f high quality music signals”, Proc. AES'91及び
T.Moriya, H.Suda：“An 8 Kbit/s transform coderfor
noisy channels,”Proc.ICASSP'89, pp196-199 及びIS
O/IEC標準IS-11172-3に述べられている。
【０００３】これらの符号化法では、高能率な符号化を
実現するためには、残差係数はできるだけ周波数特性が
平坦であることが望ましい。このため、上述のASPECやM
PEG-Audio Layer III では、周波数領域係数をいくつか
の小帯域に分割し、小帯域内の信号を、帯域の強さを表
すスケーリングファクタと呼ばれる値で割ることによっ
て正規化する事によって周波数特性の平坦化をはかる。
即ち図１に示すように入力端子１１から入力されたディ
ジタル化された音響入力信号は時間−周波数変換部（変
形離散コサイン変換：ＭＤＣＴ）２により周波数領域係
数に変換され、この周波数領域係数は帯域分割部３によ
り複数の小帯域に分割され、これら小帯域信号はそれぞ
れ代表値計算・量子化部４1 〜４n でその平均値又は最
大値などの帯域の強さを表す代表値（スケーリングファ
クタ）が計算され、かつそのスケーリングファクタは量
子化されて、全体として周波数領域係数の概形が得られ
る。前記分割された各小帯域信号は正規化部５1 〜５n
でそれぞれ対応する帯域の前記量子化されたスケーリン
グファクタで正規化され帯域残差係数とされ、これら正
規化によって得られた帯域残差係数は量子化部６で帯域
合成されスペクトラム残差を得る。その結果、時間−周
波数変換部２により得られた前記周波数領域係数はその
周波数特性の概形が取り除かれ、平坦化された残差係数
となり、その残差係数は量子化される。この残差係数の
量子化を示すインデックスＩR と、前記各代表値を量子
化したインデックスとがそれぞれ復号器へ送出される。
【０００４】この方法よりも高能率な周波数領域係数の
平坦化方法として、線形予測分析を用いる方法がある。
周知のように、線形予測係数は入力信号の周波数特性を
平坦化するように動作する線形予測フィルタ（逆フィル
タと呼ばれている）のインパルス応答を表している。こ
の方法では、図２に示すように端子１１に与えられたデ
ィジタル音響信号を線形予測分析・予測係数量子化部７
で線形予測し、得られた線形予測係数を線形予測分析フ
ィルタ、いわゆる逆フィルタ８にフィルタ係数として設
定し、この線形予測分析フィルタ８を端子１１からの入
力信号で駆動することによって周波数特性の平坦化され
た残差信号を得る。この残差信号を時間−周波数変換部
（離散コサイン変換：ＤＣＴ）２で周波数領域の信号、
即ち残差係数に変換し、残差量子化部６で量子化し、そ
の量子化を表すインデックスIRと線形予測係数を量子化
したインデックスIPとを復号器へ送出する。この方法
は、TCWVQ で用いられている。
【０００５】
【発明が解決しようとする課題】上述の何れの方法にお
いても、周波数特性の全体的な概形を正規化するにとど
まり、オーディオ信号に含まれるピッチ成分のような微
視的な周波数特性の凹凸を能率良く除去することはでき
ない。このことが障害となり、ピッチ成分の強い音楽信
号や音声信号などを符号化する際の性能の向上が困難で
あった。
【０００６】線形予測分析については斉藤、中田“音声
情報処理の基礎”第６章（オーム社）に、ＤＣＴについ
てはK.R.Rao, P.Yip著安田、藤原訳“画像符号化技術−
ＤＣＴとその国際標準”第２章（オーム社）に、および
ＭＤＣＴについてはISO/IEC標準IS-11172-3に述べられ
ている。この発明の目的は、変換符号化において、入力
音響信号の周波数特性をその概形で正規化して得られる
残差係数にピッチ成分が含まれる場合でも少ない情報量
で能率良く符号化可能な符号化方法及びその符号化音響
信号を復号する方法を提供することにある。
【０００７】
【課題を解決するための手段】この発明の符号化方法に
よれば、入力音響信号を周波数領域に変換して符号化す
る方法において、入力音響信号についてフレーム単位
で、その周波数特性の概形が平坦化された残差係数を得
る第１の段階と、その第１の段階で得られた現フレーム
の残差係数の概形を現在、または過去のフレームの残差
係数から予測し、予測残差概形を生成する第２の段階
と、第１の段階で得られた現フレームの上記残差係数
を、第２の段階で得られた予測残差概形で正規化して微
細構造係数を得る第３段階と、その第３段階で得られた
微細構造係数を量子化し、その量子化微細構造を表すイ
ンデックスを上記音響信号の符号化出力の少なくとも１
部として出力する第４の段階とを有している。
【０００８】第１の段階で残差係数を得るには、入力音
響信号を周波数領域係数に変換した後にその周波数特性
の概形を取り除くか、入力音響信号を時間領域でその周
波数特性の概形を平均化した後に、周波数領域係数に変
換してもよい。第２の段階で予測残差概形を生成するに
は、量子化微細構造係数を逆正規化して再生残差係数を
生成し、その再生残差係数からそのスペクトラム概形を
求め、そのスペクトラム概形から次フレームの残差係数
に対する予測残差概形を合成する。
【０００９】上記符号化方法の第２段階において、予測
残差概形が現フレームにおける残差係数のスペクトラム
概形に最も近くなるよう、そのスペクトラム概形を量子
化し、その量子化を表すインデックスを符号化出力の他
の一部として出力してもよい。その場合、現フレームの
残差係数のスペクトラム概形と、過去のフレームの量子
化スペクトラム概形とを予め決めた予測係数を使って線
形合成し、その線形合成値が上記スペクトラム概形と最
も近くなるように上記量子化スペクトラム概形を決定
し、その時の上記線形合成値を上記予測残差概形として
もよい。あるいは、現フレームの量子化スペクトラム概
形と、過去のフレームの上記予測残差概形とを線形合成
し、その線形合成値が上記スペクトラム概形と最も近く
なるように上記量子化スペクトラム概形を決定し、その
時の上記線形合成値を上記予測残差概形として得てもよ
い。
【００１０】上記符号化方法において、重ね合わせ直交
変換を使って入力音響信号を周波数領域係数に変換して
もよい。その場合は、入力音響信号を線形予測分析して
得た線形予測係数のスペクトラム振幅を上記周波数領域
係数の概形として求め、それを使って周波数領域係数を
正規化するのが好ましい。この発明の復号化方法によれ
ば、周波数領域に変換されて符号化された音響信号を復
号する方法において、入力された第１の量子化インデッ
クスから再生した微細構造係数を、過去のフレームの情
報から合成した残差概形で逆正規化して現フレームの再
生残差係数を得る第１の段階と、その第１の段階で得ら
れた残差係数から元の音響信号の周波数特性の概形が与
えられた音響信号を再生する第２の段階とを有する。
【００１１】第１の段階は、上記再生残差係数から次フ
レームに対する残差係数の概形を合成する第３段階を含
む。第３段階は更に、上記再生残差係数のスペクトラム
概形を計算する第４段階と、現フレームに対し、予め決
めた１又は連続する複数の過去のフレームの上記スペク
トラム概形にそれぞれ予測係数を乗残して線形合成によ
り現フレームの残差係数の概形を得る第５段階を含む。
【００１２】第２の段階において周波数特性の概形が与
えられた音響信号を再生するには、周波数領域の再生残
差係数に与える場合と、再生残差係数を時間領域に変換
して得た残差信号に与える場合とがある。上記復号方法
において、上記残差概形は、符号側から送られたインデ
ックスを逆量子化して得た現フレームと過去のフレーム
の量子化スペクトラム概形を線形合成して得てもよい。
【００１３】あるいは、上記残差概形は、過去のフレー
ムにおける残差概形と、符号側から送られたインデック
スを逆量子化して得た現フレームの量子化スペクトラム
概形を線形合成して得てもよい。
【００１４】
【作用】一般に、周波数領域係数をスペクトラム概形で
正規化して得られる残差係数は、ピッチ成分を含んでお
り、全体のパワーに対してエネルギーの大きいスパイク
となって現れる。ピッチ成分は長時間持続することから
スパイクは複数フレームにわたって同じ位置に現れるの
で、残差係数のパワーはフレーム間で相関が高い。この
発明では、前フレームの残差係数の振幅またはその概形
と現フレームの相関を用いて残差係数の冗長度を取り除
くため、つまり前記スパイクを取り除き、残差係数より
も平坦化された微細構造係数を量子化するため高い量子
化能率が得られる。また、入力信号が複数のピッチをも
っている場合にも周波数領域上ではピッチは分離されて
いるので対応できる。
【００１５】
【実施例】図３にこの発明による符号化方法を適用した
符号器１０と復号化方法を適用した復号器５０とを示
し、図中のＡ，Ｂ，…，Ｅで示す出力の波形例を図５Ａ
〜５Ｅに示す。この発明においても、まず、入力音響信
号の符号化に必要なビット数を低減するため、入力信号
に対し、まず平坦化されたスペクトラム概形を有する残
差係数を求めるが、その方法としては例えば次の２つが
ある。
【００１６】(a) 入力信号を周波数領域係数に変換する
と共に、入力信号のスペクトラム概形を求め、前記周波
数領域係数をそのスペクトラム概形で正規化して残差係
数を得る。 (b) 線形予測係数で制御される逆フィルタにより時間領
域で入力信号を処理して残差信号を求め、その残差信号
を周波数領域係数に変換し、残差係数を得る。
【００１７】上記方法(a) において、入力信号のスペク
トラム概形を得る方法として次の３つの方法が考えられ
る。 (c) 上述の事実を応用して、入力信号の線形予測係数を
フーリエ変換して求める。 (d) 図１で説明したように、入力信号を変換して得た周
波数領域係数を複数の帯域に分割し、それぞれの帯域の
スケーリングファクタをスペクトラム概形として得る。
【００１８】(e) 入力信号を変換して得た周波数領域係
数の絶対値を逆変換して得た時間領域信号の線形予測係
数を求め、その線形予測係数をフーリエ変換して求め
る。方法(c) と(e) は次の事実に基づいている。前述の
ように、線形予測係数は、入力信号の周波数特性を平坦
化するように動作する逆フィルタのインパルス応答（周
波数特性）を表しており、従って線形予測係数のスペク
トラム概形は、入力信号のスペクトラム概形に対応して
いる。詳しく言えば、線形予測係数をフーリエ変換して
得られるスペクトラム振幅は、入力信号のスペクトラム
概形の逆数となっている。
【００１９】この発明は上記方法(a) と(c),(d),(e) の
どの組み合わせを使ってもよいし、方法(b) のみを使っ
てもよい。 図３の実施例は上記方法(a) と(c) の組み
合わせを使った場合である。符号器１０には入力端子１
１からディジタル信号とされた音響信号が入力され、フ
レーム分割部１４でそのＮ入力サンプルごとに過去２×
Ｎサンプルの入力系列を抽出し、この２×Ｎ個のサンプ
ルの系列を重ね合わせ直交変換（ＬＯＴ：Lapped Ortho
gonal Transform） 処理用フレームとする。このＬＯＴ
処理用フレームは時間窓掛部１５で時間窓が掛けられ
る。ＬＯＴについては例えばH.S.Malvar,“Signal Proc
essing with Lapped Transform, ”ArtechHouseに述べ
られている。例えば０から始まってｎ番目の窓関数の値
W(n)は式(１)で表わされるものを用いるのが一般的であ
り、この実施例でもこれを用いる。
【００２０】 W(n)＝ sin｛π(n+0.5)／(2N)｝ (1) この時間窓が掛けられた信号は一方で重ね合わせ直交変
換の一種であるＮ次のＭＤＣＴ(Modified Discrete Cos
ine Transform:変形離散コサイン変換）部１６で変形離
散コサイン変換されて周波数領域係数（周波数軸上のそ
れぞれの点におけるサンプル値）に変換され、例えば図
４Ａに示すようなスペクトラム振幅が得られる。時間窓
掛け部１５の出力は他方で線形予測分析部１７で線形予
測分析され、Ｐ次の予測係数α0,…,αPが求められる。
この予測係数α0,…,αPを量子化部１８で例えばＬＳＰ
パラメータ、あるいはｋパラメータに変換してから量子
化してスペクトラム概形を示すインデックスIPを得る。
【００２１】前記予測係数α0,…,αPのスペクトラム概
形がスペクトラム概形計算部２１で求められる。図４Ｂ
は得られたスペクトラム概形の例を示す。予測係数のス
ペクトラム概形を得るには、例えば図５Ａに示すよう
に、Ｐ＋１個の量子化予測係数（αパラメータ）の後
に、(4×N-P-1)個の０をつなげて得た長さ４×Ｎのサン
プル系列を離散フーリエ変換（例えば高速フーリエ変換
ＦＦＴ）し、更にその２×Ｎ次のパワースペクトラムを
計算し、このスペクトラムの各奇数次をそれぞれ取り出
し、それらについてそれぞれ平方根をとる。このように
して得られたＮ点のスペクトラム振幅は、スペクトラム
概形の逆数を表している。あるいは図５Ｂに示すように
Ｐ＋１個の量子化予測係数（αパラメータ）の後に(2×
N-P-1)個の０をつなげた長さ２×Ｎのサンプル系列をＦ
ＦＴ分析し、その結果についてＮ次のパワースペクトラ
ムを計算する。０番目から始まってｉ番目のスペクトラ
ム概形の逆数は、ｉ＝Ｎ−１以外ではｉ＋１番目とｉ番
目の各パワースペクトラムの平方根を平均して、つまり
補間して得る。Ｎ−１番目のスペクトラム概形の逆数
は、Ｎ−１番目のパワースペクトラムの平方根をとって
得る。正規化部２２はこのようにして得られたスペクト
ラム概形で、ＭＤＣＴ部１６からのスペクトラム振幅を
各対応サンプルごとに割って正規化し、図４Ｃに示すよ
うな現フレームＦの残差係数R(F)を得る。ただし、前述
のように量子化予測係数αをフーリエ変換して直接得ら
れるものはスペクトラム概形の逆数であり、従って実際
には正規化部２２はＭＤＣＴ部１６の出力とスペクトラ
ム概形計算部２１の出力（スペクトラム概形の逆数）を
単に乗算すればよい。しかしながら、以降の説明におい
ても便宜上、スペクトラム概形計算部２１はスペクトラ
ム概形を出力するものとする。
【００２２】従来においては上述とは異なる方法で得た
残差係数を量子化してそのインデックスを送出するが、
一般に音響信号（特に音声信号、音楽信号）の残差係数
は図４Ｃに示したようにピッチ成分などの比較的大きな
変動を含んでいる。そこでこの発明では、過去又は現在
のフレームの残差係数から予測した現フレームの残差係
数R(F)の概形ER(F) で、現フレームの残差係数R(F)を正
規化して微細構造係数を得、その微細構造係数を量子化
する。この実施例では正規化により得られた微細構造係
数をそのレベルが大きい成分程、重要視するように重み
付き量子化する場合であり、スペクトラム概形計算部２
１で得られたスペクトラム概形と残差概形計算部２３で
得られた残差係数概形ER(F) とを重み計算部２４で各対
応サンプルごとに掛け合わせて重み付け係数w1,…,w
N（ベクトルＷで表す）を得て量子化部２５へ供給す
る。この重み付け係数に聴感制御を施すこともできる。
この実施例では、重み付け係数に０.6前後の定数をべき
乗する。この他の聴感制御方法として、エムペグ−オー
ディオ方式で用いられている聴覚モデルによって求めた
各サンプルごとに必要なＳＮＲ（信号対雑音比）を非対
数化して重み付け係数と掛け合わせる方法などもある。
この方法では、入力信号を分析して得られた周波数特性
から、各周波数サンプルごとに聴感的にノイズが検知で
きる最小のＳＮＲを、聴覚モデルによってマスキング量
を推定することによって計算する。このＳＮＲが各サン
プルごとに必要なＳＮＲである。エムペグ−オーディオ
における聴覚モデルの技術についてはＩＳＯ／ＩＥＣ標
準IS-11172-3に述べられている。
【００２３】正規化部２２からの現フレームＦの残差係
数R(F)を残差係数概形計算部２３からの残差係数概形で
割って正規化して微細構造係数を得ることが残差概形正
規化部２６で行われる。その現フレームの微細構造係数
はパワー正規化部２７でその振幅の平均値またはパワー
の平均値の平方根である正規化ゲインg(F)で割算すこと
により正規化され、正規化微細構造係数X(F)＝(x1,…,x
N)として量子化部２５へ供給される。このパワー正規化
の際の正規化ゲインg(F)は逆正規化部３１に与えられる
と共に量子化され、その量子化ゲインを示すインデック
スIGが出力される。
【００２４】量子化部２５では正規化微細構造係数X(F)
に対し重み付け係数Ｗにより重み付けを行った後、ベク
トル量子化する。この例ではインターリーブ型重み付き
ベクトル量子化する場合で、まず、Ｎ個のサンプルから
成る正規化微細構造係数xj(j=1,…,N)と重み付け係数wj
(j=1,…,N)の系列を、それぞれ同じように次のようにＮ
／Ｍ個ずつのＭ個の小系列にインターリーブして再配列
する。即ちそれぞれの係数のｋ番目の小系列のｉ番目の
サンプル値xk i ，wk i と、もとの系列のｊサンプル目の
値xj，wjとの関係を次式(２)に示す。
【００２５】 xk i＝xiM+k， wk i＝wiM+k， (２) 即ち、ｊ＝iM＋ｋの関係がある。ただし、ｋ＝0,1,…,M
-1であり，ｉ＝0,1,…,(N/M)-1である。Ｎ＝１６、Ｍ＝
４の場合の式(２)のインターリーブ法で正規化微細構造
係数xj(j=1,…,N)を再配列する場合の再配列前と後の関
係を図６に示す。重み付け係数wjの場合も同様である。
このようにして得られた微細構造係数と重み付け係数の
Ｍ個の小系列対の各々について重み付きベクトル量子化
する。インターリーブ後のｋ番目の小系列微細構造係数
のサンプル値をxk i ，重み付け係数の値をwk i ，符号帳
のインデックスｍのベクトルC(m)のｉ番目の要素の値を
ci(m) としたときのベクトル量子化の際の重み付き距離
尺度dk(m) は次式 dk(m)＝Σ〔wk i｛xk i−ci(m)｝〕2 (３) で規定され、Σはｉ＝０から(N/M)-1 までの加算オペレ
ータである。この距離尺度dk(m) が最小になる符号ベク
トルC(mk) を探索する事をｋ＝1,…,Mに付いて行い、そ
れぞれの符号ベクトルのインデックスm1,…,mMから量子
化インデックスImを得る。
【００２６】図７は上述のインターリーブ型重み付きベ
クトル量子化を行う量子化部２５の構成を示す。この図
を参照してｋ番目の小系列xk i に付いての量子化手順を
説明する。入力微細構造係数xjと重み付け係数wj（j=1,
…,N）はインターリーブ部２５Ａで式(２)のように再配
列され、ｋ番目の小系列xk i ，wk i がそれぞれ減算部２
５Ｂ及び２乗器２５Ｅに与えられる。符号帳２５Ｃから
選択したベクトルC(m)の要素系列ci(m) と微細構造係数
小系列xk i との差が減算部２５Ｂで求められ、その差が
２乗器２５Ｄで２乗される。一方、重み付け係数小系列
wk i は２乗器２５Ｅで２乗され、その出力と２乗器２５
Ｄの出力との内積が内積計算部２５Ｆで求められ、その
内積値dk i が最小となるベクトルC(mk) を符号帳２５Ｃ
から探索することが最適符号帳探索部２５Ｇで行われ、
そのdk i が最小となるベクトルC(mk) を示すインデック
スmkが出力される。
【００２７】このようにして、量子化部２５で量子化に
より得られたＭ個のベクトルC(m1),C(m2),…,C(mM)を構
成する要素系列である量子化小系列C(m)は、逆正規化部
３１で式(２)に従って元の配列に並べかえられて量子化
正規化微細構造係数とされ、これに対しパワー正規化部
２７で得た正規化ゲインg(F)で逆正規化し、更に残差概
形計算部２３からの残差係数概形が乗算されて量子化残
差係数Rq(F) が再生される。その量子化残差係数の概形
が残差概形計算部２３により求められる。
【００２８】残差概形計算部２３の具体例を図８を参照
して説明する。この例では、残差概形正規化部２６へ入
力された現フレームＦの残差信号R(F)に対し、残差概形
計算部２３が前フレームＦ−１の残差係数R(F-1)を使っ
て決定した予測係数β1(F-1)〜β4(F-1)を使って合成し
た残差概形ER(F) で正規化する場合である。残差概形計
算部２３の線形合成部３７は、この例では４つの縦続接
続された１フレーム遅延段３５1〜３５4と、それぞれの
遅延段の出力Ｅ1〜Ｅ4に予測係数β1〜β4を乗算する乗
算器３６1〜３６4と、全ての乗算結果の対応するサンプ
ルをそれぞれ加算し、加算結果を合成残差概形係数ER"
(F)（Ｎサンプル）として出力する加算器３４とから構
成されている。現フレームＦにおいて、遅延段３５1 〜
３５4 はそれらの出力E1(F)〜E4(F)として過去のフレー
ム(F-1)〜(F-4)に測定した残差スペクトラム概形E(F-1)
〜E(F-4)を出力しており、予測係数β1〜β4としては、
前フレーム(F-1) で決定された値β1(F-1)〜β4(F-1)が
設定されている。従って、現フレームの加算器３４の出
力ER"(F)は次式 ER"(F)＝β1(F-1)E(F-1)＋β2(F-1)E(F-2)＋…＋β4(F-1)E(F-4) で表される。
【００２９】図８の例では、加算器３４の出力ＥR は必
要に応じて定数加算部３８にて同一定数が各サンプルに
加算されて予測残差概形係数ER' を得る。定数加算部３
８で定数を加算するのは加算器３４の出力として得られ
る予測残差概形係数ＥR の予測が大きく間違っていた場
合に、その影響を小とするためである。この加算定数
は、例えば加算器３４の出力の１フレームの平均パワー
に0.05を乗算した値であり、加算器３４の出力である予
測残差概形係数ＥR の振幅の平均が１０２４のとき、前
記加算定数は５０程度とされる。定数加算部３８の出力
ER' は必要に応じて正規化部３９で１フレーム（Ｎ点）
のパワーの平均が１になるように正規化して最終的な現
フレームＦの予測残差概形（以降単に残差概形とも呼
ぶ）ER(F) を得る。
【００３０】この様にして得られる残差概形ER(F) は、
例えば図４Ｄに示すように、正規化部２２からの図４Ｃ
に示す残差係数R(F)中の強いピッチ成分と対応する位置
に単極性のインパルスを有している。一般に、オーディ
オ信号では隣接フレーム間でピッチ成分の周波数位置が
ほとんど変わらないので、残差概形正規化部２６におい
てこの残差概形ER(F) により入力残差係数R(F)を割算す
ることにより、ピッチ成分レベルが抑圧され、図４Ｅに
示すようなランダム成分が主体の微細構造係数が得られ
る。この正規化により得られた微細構造係数が前述のよ
うにパワー正規化部２７、量子化部２５で順次処理さ
れ、それぞれ正規化ゲインg(F)と量子化小系列ベクトル
C(m)が逆正規化部３１に与えられる。逆正規化部３１で
は、再生部３１Ａでは、量子化小系列ベクトルC(m)を再
配列して量子化正規化微細構造係数Xq(F)を求め、それ
に正規化ゲインg(F)を乗算して量子化微細構造係数を再
生し、その再生出力に乗算器３１Ｂで現フレームＦの残
差概形ER(F) を乗算することにより量子化残差係数R
q(F) を再生する。現フレームＦにおいてこの再生され
た量子化残差係数（再生残差係数）Rq(F) が残差概形計
算部２３のスペクトラム振幅計算部３２に与えられる。
【００３１】スペクトラム振幅計算部３２は、逆正規化
部３１からの再生された量子化残差係数Rq(F) のＮ個の
サンプルのスペクトラムの振幅を計算する。その計算さ
れたＮ個のスペクトラム振幅に対し、窓関数畳み込み部
３３で周波数窓関数が畳み込まれて現フレームＦの再生
残差係数Rq(F) の振幅の概形、即ち残差係数概形E(F)を
得て線形合成部３７に与える。スペクトラム振幅計算部
３２では、例えば、再生残差係数Rq(F) の各サンプルの
絶対値をスペクトラム振幅とする方法や、再生残差係数
Rq(F) の各サンプルの２乗値と１つ前のフレーム(F-1)
の残差係数Rq(F-1) の対応するサンプルの２乗値との和
の平方根をスペクトラム振幅とする方法などがある。ま
た、このスペクトラム振幅を対数で表してもよい。窓関
数の畳み込み方法は、幅が３〜９サンプル程度で、窓関
数の形状として、三角窓、ハミング窓、ハニング窓、指
数関数窓などを用いることもできるし、窓形状を適応的
に可変にしても良い。窓関数として指数関数窓を使う場
合の具体例としては、ｇを１以上の予め決めた整数とす
ると、窓関数を例えば次式
【００３２】
【数１】 で規定し、ａの値として例えばａ＝０.5が用いられる。
上式の場合の窓の幅は2g+1である。窓関数の畳み込みに
より、周波数軸上の各点におけるサンプル値は、その正
方向に隣接するｇ個のサンプル値と負方向に隣接するｇ
個のサンプル値の影響を受けた値に変換される。これに
より、残差概形計算部２３での残差概形の予測効果が過
敏にならないようにでき、従って復号音に異常音が生じ
るのを抑えることができる。窓の幅を１２サンプル以上
とすると、残差係数概形中のピッチ成分にもとづく変動
が不明となり、好ましくない。
【００３３】窓関数が畳み込まれて得られたスペクトラ
ム概形E(F)は、現フレームのスペクトラム概形E0(F) と
して線形合成部３７に与えられると共に、予測係数計算
部４０にも与えられる。予測係数計算部４０は線形合成
部３７への入力E0(F) と遅延段３５1 〜３５4 の出力E1
=E(F-1) 〜E4=E(F-4) が与えられ、後述のようにスペク
トラム概形E0(F) に対する加算器３４の出力ER" の自乗
誤差が最小となるように予測係数β1(F)〜β4(F)を適応
的に決定する。その後、遅延段３５1〜３５4は与えられ
ているスペクトラム概形E0〜E3をそれぞれ取り込んで、
更新されたスペクトラム概形E1〜E4として出力し、１フ
レームに付いての処理サイクルを終了する。その結果得
られる加算器３４の出力（合成残差概形）ER" に基づい
て、前述と同様に次フレーム(F+1) の残差係数R(F+1)に
対する予測残差概形係数ER(F+1)が生成される。
【００３４】予測係数β1〜β4は以下のようにして求め
ることができる。図８では予測次数は４次となっている
が、ここでは一般化のためにＱ次とする。ｑを１≦ｑ≦
Ｑを満足する任意の整数とし、ｑ段目における予測係数
の値をβq とする。乗算器３６1〜３６Q（Ｑ＝４）に対
する予測係数（乗算係数）をβ1〜βQとし、また、ｑ段
目出力の係数系列をベクトルＥq で表し、遅延段３５1
〜３５Qの出力をＥ1 ，Ｅ2 ，…，ＥQ とし、さらに、
窓関数畳み込み部３３の出力であるスペクトラム概形の
係数系列（現フレームの残差係数概形係数）E(F)をベク
トルＥ0 で表す。この時、次式(４)で求められる相互相
関係数ｒを用いて、連立一次方程式(５)をβ1〜βQにつ
いて解くことによって加算器３４の出力ＥR のＥ0 に対
する自乗誤差（予測誤差）が最小になる予測係数β1 〜
βQ が求められる。
【００３５】 ｒi,j ＝Ｅi ・Ｅj (４)
【００３６】
【数２】 線形合成部３７において参照する過去のフレームのスペ
クトラム概形は４フレーム前までに限ることなく、１フ
レーム前だけでも、それ以上前でもよく、従って遅延段
の数Ｑは１以上任意の数でよい。
【００３７】この様に、図８の残差概形計算部２３を使
った符号化方法によれば、正規化部２２からの残差係数
R(F)は、その過去の残差係数から推定された残差概形ER
(F)で正規化されるため、その正規化された微細構造係
数は残差係数R(F)よりも概形の変化が少ない平坦なもの
になる。従って、それだけ少ないビット数で量子化する
ことができる。また残差係数のスペクトラム振幅系列に
対し窓関数畳み込み部３３で周波数窓関数を畳み込んで
求めたスペクトラム概形E(F)から予測した残差概形E
R(F) で残差係数R(F)を正規化するので、残差係数R(F)
中の例えばピッチ成分とそれぞれ対応した位置に現れる
強いパルスに対し、残差概形の推定が周波数軸方向に１
サンプル程度ずれても大きな予測誤差を生ずることなく
動作する。窓関数畳み込みを用いない場合は、推定誤り
が生じると大きな予測誤差が生じる原因となる。
【００３８】図３において、符号器１０から線形予測係
数の量子化値を示すインデックスIPと、微細構造係数の
パワー正規化利得の量子化値を示すインデックスIGと、
微細構造係数の量子化値を示すインデックスImとが出力
される。復号器５０においてはインデックスIP，IG，Im
が入力され、再生部５１で入力されたインデックスImか
ら正規化微細構造係数が再生され、正規化ゲイン再生部
５２で入力されたインデックスIGから正規化ゲインが再
生され、パワー逆正規化部５３で前記再生された正規化
微細構造係数が前記再生された正規化ゲインにより逆正
規化されて微細構造係数とされる。その微細構造係数
は、残差逆正規化部５４で、残差概形計算部５５からの
残差概形ＥR が掛算されて逆正規化され、残差係数R(F)
が再生される。この再生された残差係数R(F)の概形が残
差概形計算部５５で符号器１０の残差概形計算部２３と
同じ方法で計算される。
【００３９】一方、再生・スペクトラム概形計算部５６
で入力されたインデックスIPから線形予測係数α0〜αP
が再生され、更にその線形予測係数から、符号器１０の
スペクトラム概形計算部２１と同じ方法でスペクトラム
概形が計算され、その計算されたスペクトラム概形が逆
正規化部５７で残差概形逆正規化部５４からの再生残差
係数R(F)に対して掛算されて逆正規化され、周波数領域
係数が再生される。その周波数領域係数はフレーム毎に
逆ＭＤＣＴ部５８でＮ次の逆変形離散コサイン変換がさ
れて２Ｎサンプルの時間領域信号（逆ＬＯＴ処理用フレ
ームと呼ぶ）に変換される。この時間領域信号は窓掛部
５９でフレームごとに例えば式(１)で表わされる形状の
時間窓がかけられる。その窓かけ出力はフレーム重ね合
わせ部６１で長さ２×Ｎサンプルの逆ＬＯＴ処理用現フ
レームの前半Ｎサンプルと前フレームの後半Ｎサンプル
とが互いに加算され、得られたＮサンプルを現フレーム
の再生音響信号として出力端子９１に出力される。
【００４０】上述においてＰ，Ｎ，Ｍの各値は、Ｐ＝６
０前後、Ｎ＝５１２，Ｍ＝６４程度を目安に自由に選べ
るが、Ｐ＋１＜Ｎ×４でなくてはならない。また、上記
実施例では、図６で説明した残差概形係数のインターリ
ーブベクトル量子化の際の係数系列の分割数Ｍの値は、
Ｎ／Ｍの値が整数であるものとして説明したが、Ｍの値
は、必ずしもＮ／Ｍが整数となるように設定する必要は
ない。整数とならない場合には、分割したそれぞれの小
系列の一部を１サンプルずつ長くして、不足サンプル数
を補えばよい。
【００４１】残差概形計算部２３，５５としては図９に
図８と対応する部分に同一符号を付けて示すように、窓
関数畳み込み部３３の出力を平均計算部４１で例えば１
０フレームにわたる平均を対応するサンプル毎に計算
し、又は１フレーム内平均をフレーム毎に求め、つまり
直流成分を検出し、その直流成分を窓関数畳み込みの出
力から減算器４２で差引き、得られた変動分のみを遅延
段３５1 へ供給し、加算器３４の出力に加算器４３で同
じ平均計算部４１の出力を加えてもよい。加算器３４の
出力ER" が減算器４２の出力Ｅ0 になるべく近ずくよう
に予測係数β1 〜βQ を決定する。このような予測係数
β1〜βQは、前述と同様に式(４)、(５)により決定する
ことができる。しかし図９の構成によれば、変動分につ
いてのみ予測するため予測の能率がより向上する。
【００４２】更に図１０に示すように図９中の減算器４
２の出力のフレーム内の各サンプルの２乗の平均値の平
方根、つまり標準偏差を振幅検出部４４で計算し、この
標準偏差で減算器４２の出力を割算器４５において割る
ことにより正規化して変動を平坦化したスペクトラム概
形Ｅ0 を遅延段３５1 へ供給し、加算器３４の出力ER"
が割算器４５の出力になるべく一致するように式(４)、
(５)により予測係数β 1 〜βQ を決定し、加算器３４の
出力に振幅検出部４４の出力である標準偏差を乗算器４
６で乗算して逆正規化し、その逆正規化出力を加算器４
３へ供給して残差概形係数ER(F) を得る。ただし、図１
０のように処理すると、図８に示した構成における予測
係数β1〜βQを求める方程式(５)は、次式(６)
【００４３】
【数３】 のように近似することが可能である。つまり、線形合成
部３７に与えられるスペクトラム概形のパワーが正規化
されているので、式(５)の左辺第１項中の対角要素
r1,1，r2,2，…は互いに等しくなり、またri,j＝rj,iと
なる。式(６)中の行列はテプリッツ型なので、この方程
式はLevinson-Durbin アルゴリズムで高速に解くことが
できる。また、図８及び９では、Ｑ×Ｑ個の相関係数を
計算しなくてはならなかったが、図１１ではＱ個の相関
係数を計算するだけでよいので、予測係数β1 〜βQ を
求めるのに少ない演算量ですむ。また、相関係数r
0,jは、式(４)のように求めても良いが、下式(７)のよ
うに互いにｊフレーム離れた係数ベクトルＥi とＥi+j
の内積をｉ＝０からnMAXに渡って加算して求めると、よ
り安定性が増す。
【００４４】 r0,j ＝ (1/S)ΣＥi・Ｅi+j (７) ここで、Σはｉ＝０からnMAXまでの総和演算子であり、
Ｓは平均化のための定数であり、Ｓ≧Ｑである。nMAXの
値は、Ｓ−１でも良いし、(S-j-1) でも良い。なお、Le
vinson-Durbin アルゴリズムの詳細については、「音声
情報処理の基礎」（斉藤、中田共著、オーム社）に記述
してある。
【００４５】図１０において、振幅検出部４４で標準偏
差を求める代りに、各サンプルの絶対値の平均値を用い
てもよい。図８及び９における予測係数β1 〜βQ の計
算について、相関係数ri,jを下式(８)のように計算する
こともできる。 ri,j ＝ (1/S)ΣＥn+i・Ｅn+i+j (８) ここで、Σはｎ＝０からｎMAX までの総和演算子であ
り、Ｓは平均化のための定数であり、Ｓ≧Ｑである。ｎ
MAX の値は、S-1 でも良いし、S-j-1 でも良い。このよ
うに相関係数を計算すると、ＳがＱよりも十分大きい場
合には、ri,j＝r0,jと近似でき、予測係数を求める方程
式(５)は、式(６)と同一に近似でき、Levinson-Durbin
アルゴリズムで高速に解くことができる。
【００４６】残差概形計算部２３，５５における残差概
形の予測係数β1〜βQを全帯域一括で決定していたが、
残差概形計算部２３，５５の入力をまず小帯域に分割
し、これら各小帯域毎に独立に予測係数を設定してもよ
い。この小帯域の分割方法としては、周波数領域におい
て線形尺度で等分割、対数尺度で等分割、バーク尺度で
等分割などを用いることができる。
【００４７】残差概形計算部２３，５５において、予測
係数β1 〜βQ の予測誤差の影響が小さくなるように、
窓関数畳み込み部３３における窓幅、窓中心、場合によ
っては窓形状を変更しても良い。また、窓関数の畳込み
と、予測係数β1 〜βQ による線形合成とを一括して行
ってもよい。その例を図１１に示す。この例は予測次数
Ｑを４、窓幅Ｔを３とした場合で、遅延段３５1 〜３５
4 の出力は、それぞれ周波数軸に沿って１サンプル正側
にシフトするシフタ７p1〜７p4と、１サンプル負側にシ
フトするシフタ７n1〜７n4とにそれぞれ供給され、正側
シフタ７p1〜７p4の各出力はそれぞれ乗算器８p1〜８p4
を通じて加算器３４へ供給され、負側シフタ７n1〜７n4
の各出力はそれぞれ乗算器８n1〜８n4を通じて加算器３
４へ供給される。いま乗算器３６1，８n1，８p1，３
６2，８n2，８p2，…，８p4の各乗算係数をそれぞれ
β1,β2,β3,β4,β5,β6,…,βu（この例ではｕ＝１
２）とし、その各入力スペクトラム概形ベクトルを
Ｅ1，Ｅ2，Ｅ3，Ｅ4，…，Ｅu とし、スペクトラム振幅
計算部３２の出力をＥ0 とする時、加算器３４の出力Ｅ
R のＥ0に対する２乗誤差が最小となる予測係数β1〜β
uは、予測係数計算部４０において次の速度一次方程式
（10）を解くことにより求めることができる。
【００４８】ri,j ＝Ｅi・Ｅj
【００４９】
【数４】 この様にして決定された予測係数β1〜βuに対し得られ
る加算器３４の出力Ｅ R は、図８と同様に必要に応じて
定数加算し、更に正規化して現フレームＦの残差概形ER
(F) とされ、残差概形正規化部２６で残差係数R(F)に対
する概形正規化のために使用される。このように窓関数
も適応化することは、図９及び１０の例にも適用でき
る。
【００５０】図３、８〜１１の各例では、残差概形正規
化部２６に入力された現フレームＦの残差係数R(F)に対
し、残差概形計算部２３において１つ前のフレームＦ−
１の残差係数R(F-1)を使って決定した予測係数β1(F-1)
〜βQ(F-1)（又はβu） を使って得られた予測残差係数
概形ER(F) を使って残差係数R(F)の正規化を行い、それ
によって微細構造係数を求める場合を説明した。しか
し、残差概形正規化部２６に入力された現フレームＦの
残差係数R(F)に対し、残差概形計算部２３において現フ
レームで予測係数β1(F)〜βQ(F)（図１１の場合はβu
であるが、以降全てβQで代表する） を決定し、次式 ER"(F)＝β1(F)E1(F)＋β2(F)E2(F)＋…＋βQ(F)EQ(F) により合成残差概形ER"(F)を求め、それから得られた予
測残差概形ER(F) を使って現フレームの残差係数R(F)を
正規化するように構成してもよい。その場合は、現フレ
ームＦの残差係数R(F)を直接使って予測係数β1〜βQを
決定するため、図３中に破線で示すように正規化部２２
の出力である現フレームの残差係数R(F)を直接残差概形
計算部２３に与える。この方法は図８〜１１のいずれの
残差概形計算部２３にも適用できるが、代表して図８に
適用した場合の構成を図１２に示す。
【００５１】図１２において、図８と対応する部分には
同じ番号を付けてある。図８と異なる点は、残差概形計
算部２３にもう一組のスペクトラム振幅計算部３２' と
窓関数畳み込み部３３' を設け、現フレームの残差係数
R(F)を直接スペクトラム振幅計算部３２' に与えてその
スペクトラム振幅概形を求め、それに窓関数畳み込み部
３３' で窓関数を畳み込んでスペクトラム概形Et 0(F)を
求め、それを予測係数計算部４０に与えていることであ
る。従って、再生残差係数Rq(F) から求めた現フレーム
のスペクトラム概形E0(F) は線形合成部３７の第１遅延
段３５1 のみに与えられる。
【００５２】まず、正規化部２２（図３参照）から残差
概形正規化部２６へ与えられた現フレームＦでの入力残
差係数R(F)は、残差概形計算部２３のスペクトラム振幅
計算部３２' 及び窓関数畳み込み部３３' でスペクトラ
ム振幅計算部３２及び窓関数計算部３３と同様の処理を
受け、残差係数R(F)のスペクトラム概形Et 0(F)が求めら
れ、それが予測係数計算部４０に与えられる。予測係数
計算部４０は、与えられた係数ベクトルEt 0 に対し加算
器３４の出力ER" の２乗誤差が最小となるように、図８
の場合と同様に式(４)、(５)により予測係数β1〜β4を
計算する。その係数β1〜β4が乗算器３６1〜３６4に与
えられた時の加算器３４の出力が現フレームＦの合成残
差概形ER"(F) として得られる。
【００５３】合成残差概形ER"(F)は必要に応じて定数加
算部３８、正規化部３９で図８の場合と同様に処理され
て現フレームの残差概形ER(F) として残差概形正規化部
２６に与え、現フレームＦの入力残差係数R(F)を正規化
し、微細構造係数を得る。微細構造係数は図３で説明し
たようにパワー正規化部２７でパワーが正規化され、更
に量子化部２５で重み付きベクトル量子化され、パワー
正規化部２７での正規化ゲインの量子化インデックスIG
と量子化部２５での量子化インデックスImが復号器５０
に供給される。一方、量子化部２５から出力されたイン
ターリーブ型重み付きベクトルC(m)を逆正規化部３１で
再配列し、正規化ゲインg(F)で逆正規化して得られた再
生残差係数Rq(F) は、残差概形計算部２３のスペクトラ
ム振幅計算部３２に与えられてＮ個のサンプル点のスペ
クトラム振幅をそれぞれ計算し、得られた残差係数振幅
に対し窓関数畳み込み部３３で窓関数を畳み込んで残差
係数の概形E0(F) を得る。このスペクトラム概形E0(F)
を線形合成部３７に対する現フレームＦの入力係数ベク
トルＥ0 として与え、遅延段３５1〜３５4はそれぞれ与
えられたスペクトラム概形E0〜E3を取り込んで、更新さ
れたスペクトラム概形E1〜E4として出力し、１フレーム
に付いての処理サイクルを終了する。
【００５４】図１２の実施例では、現フレームの残差係
数R(F)からβ1〜β4を決定し、これを用いて現フレーム
の残差概形ER(F) を予測合成するが、図３の復号器５０
においては、現フレームの再生残差係数Rq(F) は残差概
形逆正規化部５４においてパワー逆正規化部５３からの
現フレームの微細構造係数と残差概形計算部５５からの
現フレームの残差概形係数とを使って求めるべきもので
あり、従って残差概形計算部５５は現フレームの予測係
数β1〜β4を決定するための現フレームの残差係数R(F)
を与えられていない。従って、式(４)、(５)により予測
係数β1〜β4を決定することはできない。そこで、符号
器１０において図１２の残差概形計算部２３を使った場
合、対応する復号器５０の残差概形計算部５５には、符
号器１０側の予測係数計算部４０で決定した現フレーム
の予測係数β1〜β4を量子化してその量子化インデック
スIBを供給し、復号器５０では与えられたインデックス
IBから再生した係数β1〜β4を使って現フレームの残差
概形を計算する。
【００５５】即ち、図１３に復号器５０の残差概形計算
部５５を示すように、符号器１０の予測係数計算部４０
から与えられた現フレームの予測係数β1〜β4の量子化
インデックスIBは残差概形計算部５５の逆量子化部６０
で逆量子化されて予測係数β 1〜β4が再生され、線形合
成部６２の乗算器６６1〜６６4に設定される。これらの
予測係数β1〜β4は遅延段６５1〜６５4の出力とそれぞ
れ乗算され、それら乗算結果は加算器６７で加算され残
差概形ＥR が合成される。この合成残差概形Ｅ R は符号
器側と同様に定数加算部６８、正規化部６９で処理さ
れ、現フレームの残差概形ER(F) として残差概形逆正規
化部５４に与えられる。残差概形逆正規化部５４はパワ
ー逆正規化部５３からの現フレームの微細構造係数と上
記残差概形ER(F) を乗算して現フレームの残差係数R(F)
を求め、スペクトラム振幅計算部６３に与えると共に、
逆正規化部５７（図３）に与える。与えられた再生残差
係数Rq(F) に対しスペクトラム振幅計算部６３及び窓関
数畳み込み部６４は図１３の対応する部分と同様の処理
を行って残差係数のスペクトラム概形を得て、線形合成
部６２に与える。従って、図１２の残差概形計算部２３
に対応する復号器の残差概形計算部５５は予測係数計算
部を有していない。図１２の予測係数計算部４０におけ
る予測係数の量子化方法としては、例えば予測係数をＬ
ＳＰパラメータに変換した後、例えばフレーム間差分ベ
クトル量子化等の量子化を施すＬＳＰ量子化法を用いる
ことができる。
【００５６】図８〜１０及び１２の残差概形計算部２３
において、それぞれ加算器３４の現フレームの出力であ
る合成残差概形ＥR に対する１〜４フレーム前の残差係
数のスペクトラム概形Ｅ1 〜Ｅ4 の寄与の程度に応じ
て、それぞれ乗算器３６1 〜３６4 の乗算係数β1 〜β
4 を予め固定的に決めてもよく、例えば古いフレーム
程、重み（乗算係数）が小とされる。あるいは同一重
み、この例では１／４として４つのフレームのサンプル
係数の平均を用いてもよい。この様に係数β1 〜β4を
固定した場合は、式(４)、(５)の計算を行う予測係数計
算部４０は不要である。この場合、復号器５０側の対応
する残差概形計算部５５でも同じ係数β1 〜β 4 を使え
ばよいので係数β1 〜β4 を復号器５０に転送する必要
はない。図１１の場合も同様に係数β1 〜β12を固定し
てもよい。
【００５７】図８〜１０及び１２の残差概形計算部２３
の各例において、その構成を最も簡単にするには、例え
ば図８において、加算器３４，遅延段３５2〜３５4，乗
算器３６2〜３６4を省略して、乗算器３６1 の出力を定
数加算部３８へ直接供給し、前フレーム(F-1) のスペク
トラム概形E1＝E(F-1)からのみ、残差概形係数ER(F)を
推定してもよい。この変形は図１０にも同様に適用で
き、その場合、乗算器３６1 、８p1、８n1の出力のみを
加算器３４に与えればよい。
【００５８】図３及び８〜１２の各例では、残差概形計
算部２３は、合成残差概形係数ER"が、入力された再生
残差係数Rq(F) 又は残差係数R(F)から求めたスペクトラ
ム概形E(F)とできるだけ一致するように線形予測により
予測係数β（β1,β2,…）を決めて予測残差概形係数ER
(F) を求めていたが、この様な線形予測を行わないで、
残差概形を決定する実施例を図１４、１５及び１６を参
照して説明する。
【００５９】図１４は図３と対応させて符号器１０と復
号器５０の全体の構成を示し、残差概形計算部２３に対
する接続は図３における破線の接続に対応する。従っ
て、図１２の場合と同様な逆正規化部３１は設けられて
いない。図３及び１２と異なる点は、図１４の実施例に
おける残差概形計算部２３は、線形合成すべき残差概形
ＥR が、与えられた残差係数R(F)のスペクトラム概形と
できるだけ一致するようにスペクトラム概形を量子化
し、その時の線形合成出力ＥR を残差概形ER(F) として
使うと共に、その時の量子化インデックスIQを復号器５
０に供給する。復号器５０は、与えられたスペクトラム
概形量子化インデックスIQから残差概形計算部５５で逆
量子化処理によりスペクトラム概形E(F)を再生し、残差
概形逆正規化部５４に与える。その他の各部の処理は図
３の場合と同様であり、説明を省略する。
【００６０】図１４の実施例における符号器１０と復号
器５０の残差概形計算部２３及び５５の構成例を図１５
に示す。残差概形計算部２３は残差係数R(F)が与えら
れ、そのＮサンプル点のスペクトラム振幅をそれぞれ計
算するスペクトラム振幅計算部３２と、得られたＮ点の
スペクトラム振幅に窓関数を畳み込み演算し、スペクト
ラム概形E(F)を求める窓関数畳み込み部３３と、得られ
たスペクトラム概形E(F)を量子化する量子化部３０と、
量子化スペクトラム概形が現フレームの量子化スペクト
ラム概形係数Ｅq0として入力され、過去のフレームの量
子化スペクトラム概形係数と線形合成を行う線形合成部
３７とから構成されている。線形合成部３７は、図１２
に示すものとほぼ同様に、遅延段３５1〜３５4と、乗算
器３６1 〜３６4 と、加算器３４とから構成されてい
る。図１５の場合、１〜４フレーム前の量子化スペクト
ラム概形係数Ｅq1〜Ｅq4をそれぞれ予測係数β1 〜β4
で乗算した結果だけでなく、乗算器３６0 で現フレーム
の入力量子化スペクトラム概形係数Ｅq0に予測係数β0
を乗算した結果も加算器３４に与えて合成し、予測残差
概形ER(F) として出力する。また、予測係数β0 〜β4
は予め決めた固定値である。量子化部３０は残差概形ER
(F) の、入力スペクトラム概形E(F)に対する２乗誤差が
最小となるようにスペクトラム概形E(F)を量子化し、得
られた量子化スペクトラム概形係数Eq0 を線形合成部３
７に与えると共に、その量子化インデックスIQを復号器
５０の残差概形計算部５５に供給する。
【００６１】残差概形計算部５５の再生部６０は、与え
られた量子化インデックスIQから現フレームの量子化ス
ペクトラム概形係数を再生し、符号器１０側と同様に遅
延段６５1 〜６５4 と乗算器６６0 〜６６4 と加算器６
７とから構成された線形合成部６２は再生部６０からの
現フレームの量子化スペクトラム概形係数と、遅延段６
５1 〜６５4 からの過去のフレームの量子化スペクトラ
ム概形係数とを線形合成する。加算器６７から合成され
た残差概形ER(F) が出力され、残差概形逆正規化部５４
に与えられる。なお、乗算器６６0 〜６６4 には符号器
１０側と同じ係数β0〜β4が与えられている。また、符
号器１０の量子化部３０における量子化法として、スカ
ラ量子化を行ってもよいし、ベクトル量子化を行っても
よい。後者の場合、図７で説明したようなインターリー
ブされた係数系列のベクトル量子化を使ってもよい。
【００６２】図１６は図１５の変形実施例を示し、図１
５と対応する部分には同じ参照番号を付けてある。この
実施例では、量子化部３０においてスペクトラム概形E
(F)に対する予測残差概形（加算器３４の出力）ER(F)
の２乗誤差が最小となるように量子化を行う点は同じで
あるが、線形合成部３７の構成が異なっている。即ち、
縦続接続された遅延段３５1〜３５4の入力として予測残
差概形ER(F) が与えられ、従って、遅延段３５1 〜３５
4 からそれぞれ１〜４フレーム前の予測残差概形ER(F-
1)〜ER(F-4)が出力されている。また、量子化部３０か
ら出力される量子化スペクトラム概形Eq(F) が直接加算
器３４に入力されている。従って、線形合成部３７は過
去のフレーム(F-1)〜(F-4)における予測残差概形ER(F-
1)〜ER(F-4)と現フレームＦにおける量子化概形係数と
を線形合成し、現フレームの予測残差概形ER(F) を出力
する。復号器５０側の線形合成部６２も同様の構成とさ
れ、過去のフレームにおける合成残差概形と、再生され
た現フレームにおける量子化概形係数とを線形合成する
ことにより現フレームの残差概形を得る。
【００６３】図８〜１２、１５及び１６の各残差概形計
算部２３において、帯域処理部を設け、窓関数畳み込み
部３３からのスペクトラム概形を複数の帯域に区分し、
変動成分の少ない高次帯域に付いてそのスペクトラム概
形を振幅が一定の平坦な概形に近似してもよい。図１７
は、例えば図８における畳み込み部３３と遅延部３５と
の間に挿入するそのような帯域処理部４７の例を示す。
この例では窓関数畳み込み部３３の出力E(F)は帯域処理
部４７に入力され、分割部４７Ａで最低次（最低周波
数）から約２／３のサンプル点を中心に例えば５０次程
度の幅の狭い中間帯域の成分EB(F) と、これより高次の
高次帯域成分EH(F) と低次の低次帯域成分EL(F) に分割
し、高次帯域成分EH(F) を平均部４７Ｂへ供給して、そ
の高次帯域のスペクトラム振幅の平均値を求めて、高次
帯域成分EH(F) を全てその平均値に置き換えて出力し、
低次帯域成分EL(F) はそのまゝ出力する。また中間帯域
の成分EB(F) については、漸近部４７Ｃにより中間帯域
の高低両端におけるスペクトラム振幅がそれぞれ平均部
４７Ｂで得た平均値及び低次帯域EL(F)の最高次のスペ
クトラム振幅に一致し、その間で直線的に変化するよう
に各スペクトラム点の振幅が修正されて出力される。つ
まり、高周波成分は変動が少ないので高次帯域内のスペ
クトラム振幅を一定値、ここでは平均値、に近似してい
る。
【００６４】図８〜１２の各残差概形計算部２３におい
て、入力音響信号の典型的な複数の状態に応じた好まし
い予測係数β1 〜βQ （又はβu） の組をそれぞれイン
デックスに対応した係数ベクトルとして符号帳に用意し
ておき、もっとも良く残差概形を予測する係数ベクトル
を符号帳から選択して使用し、その係数ベクトルを示す
インデックスを復号器５０の残差概形計算部５５に転送
してもよい。
【００６５】図８〜１１の各実施例のように、過去の残
差係数概形から現在の残差概形を予測する線形予測モデ
ルにおいては、その系の安定性を判断するためパラメー
タｋが用いられるが、この発明においても、予測係数を
ｋパラメータに変換して、その絶対値が１.0に近い、或
いは１.0を越えている場合、強制的に予め決められた係
数に設定したり、または別の残差概形構成方法に、例え
ば図８から図９へ切り替えたり、あるいは予め決められ
た残差概形（例えば凹凸のない平坦な信号）に切り替え
る、などの処理を施して系の安定性が増すようにするこ
ともできる。
【００６６】図３及び１４の実施例の符号器１０では、
線形予測分析部１７で線形予測分析の際、窓掛け部１５
からの入力音響信号の自己相関係数を用いて予測係数を
求めた。しかし、図１８に示すようにＭＤＣＴ部１６で
求めた周波数領域係数の各サンプル（スペクトラム）の
絶対値を絶対値部８１で求め、その絶対値出力を逆フー
リエ変換部８２で逆フーリエ変換して自己相関係数と
し、その自己相関係数を線形予測分析部１７で線形予測
分析してもよい。この場合は、その分析に先立って相関
を求める必要はない。
【００６７】更に、図３及び１４の各実施例において
は、符号器側では入力信号の線形予測係数α0〜αPを量
子化し、その量子化予測係数をフーリエ変換して入力信
号のスペクトラム概形（周波数特性の概形）を求め、そ
の概形により入力信号の周波数特性を正規化して残差係
数を得た。また、その量子化予測係数のインデックスIP
を復号器に転送し、復号器でそのインデックスIPから線
形予測係数α0 〜αP を得、更にその線形予測係数から
周波数特性の概形を求めた。しかし、以下のようにして
もよい。即ち、図１９に図３と対応する部分に同一符号
を付けて示すように、ＭＤＣＴ部１６からの周波数領域
係数をスケーリングファクタ計算量子化部１９に分岐入
力して、その周波数領域係数をいくつかの小帯域に分割
し、その各小帯域ごとのサンプルの絶対値の平均値、ま
たは最大値をスケーリングファクタとして求め、このス
ケーリングファクタを量子化し、そのインデックスISを
復号器５０へ送出する。また前記求めたスケーリングフ
ァクタでＭＤＣＴ部１６からの周波数領域係数を正規化
部２２において、各対応小帯域ごとに割算して残差係数
R(F)を得て残差概形正規化部２６へ供給する。更に、前
記求めたスケーリングファクタと、残差概形計算部２３
からの残差係数の概形の対応する小帯域のサンプルとを
重み計算部２４で乗算して重み付け係数Ｗ(w1,…,wN)を
得て量子化部２５へ供給する。復号器５０では入力され
たインデックスISからスケーリングファクタ再生部７１
でスケーリングファクタを再生し、そのスケーリングフ
ァクタを逆正規化部５７で残差概形逆正規化部５４から
の再生残差係数に乗算して周波数領域係数を再生して逆
ＭＤＣＴ部５８へ供給する。
【００６８】上述では入力音響信号を周波数領域係数に
変換した後、残差係数を得たが、入力音響信号から時間
領域でスペクトラム概形を取り除いた残差信号を求め、
その残差信号を周波数領域の残差係数に変換してもよ
い。即ち、図２０に図３と対応する部分に同一符号を付
けて示すように、入力端子１１からの音響入力信号は線
形予測分析部１７で線形予測分析され、得られた線形予
測係数α0〜αPは量子化部１８で量子化され、その量子
化された線形予測係数が逆フィルタ２８にフィルタ係数
として設定され、この逆フィルタ２８に端子１１からの
入力音響信号が通されて、平坦化された周波数特性を有
する時間領域の残差信号を得る。この残差信号はコサイ
ン変換部２９で離散コサイン変換されて周波数領域の残
差係数R(F)に変換されて残差概形正規化部２６へ供給さ
れる。一方、量子化部１８からの量子化線形予測係数が
スペクトラム概形計算部２１へ供給され、入力信号の周
波数特性の概形が求められ、重み計算部２４に与えられ
る。符号器１０におけるその他の部分の処理は図３と同
様である。
【００６９】復号器５０では残差概形逆正規化部５４か
らの再生残差係数Rq(F) は逆コサイン変換部７２で逆離
散コサイン変換されて時間領域の残差信号とされて合成
フィルタ７３へ供給される。一方、入力されたインデッ
クスIPは再生部７４で線形予測係数α0〜αPに再生さ
れ、その線形予測係数が合成フィルタ７３のフィルタ係
数として設定される。合成フィルタ７３は逆コサイン変
換部７２からの残差信号が与えられ、音響信号を合成し
て出力端子９１へ出力する。図２０に示した例では時間
周波数変換はＭＤＣＴよりもＤＣＴの方が適する。
【００７０】図３、１４、１９及び２０の実施例におけ
る量子化部２５として、図２１に示すように構成し、図
２２に示す処理手順に従って量子化を行ってもよい。即
ち、まずスカラー量子化部２５Ａにおいて量子化幅制御
部２５Ｄから与えられた予め決められた最大量子化幅で
微細構造係数X(F)をスカラー量子化する(S1)。次に誤差
計算部２５Ｂにおいて入力された微細構造係数X(F)に対
する量子化された微細構造係数Xq(F) の誤差を計算する
(S2)。誤差としては、例えば重み付け係数Ｗを使った重
み付け２乗誤差を使う。次に量子化ループ制御部２５Ｃ
においてこの量子化誤差が聴感的に許容される所定値よ
り小さいかを判断する(S3)。この時、量子化誤差が所定
値より小さければ、その量子化微細構造係数Xq(F) とそ
れを表すインデックスImを出力すると共に、使用された
量子化幅を表すインデックスIDを量子化幅制御部２５Ｄ
から出力し、量子化処理を終了する。量子化ループ制御
部２５Ｃは、ステップＳ３で量子化誤差が所定値より大
きいと判断した場合は、量子化微細構造係数Xq(F) に使
用されるビット数が使用可能最大ビット数を越えていな
いか判断する(S4)。使用されるビット数が使用可能最大
ビット数を越えていなければ、処理ループを継続すべき
と判断し、量子化幅制御部２５Ｄに対し前回より小さい
所定の量子化幅をスカラー量子化部２５Ａに与えさせ(S
5)、スカラー量子化部２５Ａで再量子化する。以下同様
の処理手順を繰り帰る。ステップＳ４で使用ビット数が
使用可能最大ビット数以上の場合は、前回のループによ
る量子化微細構造係数Xq(F) とそのインデックスImと共
に量子化幅インデックスIDを出力して量子化処理を終了
する。
【００７１】この量子化部２５に対応する復号器５０の
再生部５１（図３、１４、１９、２０参照）には、量子
化インデックスImと量子化幅インデックスIDを供給し、
再生部５１はこれらのインデックスに基づいて微細構造
係数を再生する。
【００７２】
【発明の効果】以上述べたように、この発明によれば、
ピッチ成分を含む信号が入力された場合にみられる周波
数領域の残差係数におけるフレーム間の高い相関を利用
して、その残差係数の概形を正規化して、より平坦化し
た微細構造係数として量子化するため、高い量子化能率
が得られる。また、複数ピッチが存在する場合にも周波
数領域ではピッチが分離されるので問題を生じない。更
に、残差係数の概形を適応的に決定するので、ピッチ成
分の変化傾向に追従することもできる。
【００７３】また、ＭＤＳＴのような重ね合わせ直交変
換を使って入力音響信号を周波数領域係数に変換し、そ
れを、音響信号の線形予測係数から求めたスペクトラム
振幅（即ち入力音響信号の周波数特性の概形）で周波数
領域において正規化する実施例では、フレーム間ノイズ
を生じさせないで周波数領域係数の高能率な平坦化を実
現することが可能である。
【００７４】Ｐ＝６０，Ｎ＝５１２，Ｍ＝６４，Ｑ＝２
とし、線形予測係数α0〜αPと正規化ゲインの量子化の
ための情報量を十分大きくとり、微細構造係数を２ビッ
ト／サンプルの情報量でベクトル量子化する条件で、図
８に示した残差概形計算部２３を用いて様々な音楽ソー
スを符号化復号化した場合、残差概形計算部２３，５５
を取り除いて符号化復号化した場合よりも平均約５dB，
最大約１０dBのセグメンタルＳＮＲの改善がみられる。
また、聴感的にもより自然性の高い音が得られる。
【図面の簡単な説明】
【図１】入力信号の周波数特性をスケーリングファクタ
によって平坦化する従来の符号器を示すブロック図。
【図２】入力信号の周波数特性を線形予測分析フィルタ
で平坦化する従来の符号器を示すブロック図。
【図３】この発明の符号化方法が適用された符号器、及
びこの発明の復号化方法が適用された復号器の各実施例
を示すブロック図。
【図４】Ａは図３におけるＭＤＣＴ部で求めた周波数領
域波形の例、Ｂはスペクトラム概形計算部により求めた
スペクトラム概形の例、Ｃは正規化部２２により求めた
残差係数の例、Ｄは残差概形計算部２３で求めた残差概
形の例、Ｅは残差概形正規化部２６で求めた微細構造係
数の例をそれぞれ示す図。
【図５】Ａは予測係数から周波数特性の概形を得る方法
を示す図、Ｂは予測係数から周波数特性の概形を得る他
の方法を示す図。
【図６】ベクトル量子化の際の信号系列と分割小系列と
の関係例を示す図。
【図７】量子化部２５の構成例を示すブロック図。
【図８】図３中の残差概形計算部２３，５５の具体例を
示すブロック図。
【図９】図８中の残差概形計算部２３，５５の変形例を
示すブロック図。
【図１０】図９中の残差概形計算部２３、５５の変形例
を示すブロック図。
【図１１】図３中の残差概形計算部２３，５５で窓関数
及び予測係数の両者を適応的に制御するようにした例を
示すブロック図。
【図１２】図３中の残差概形計算部２３の更に他の例を
示すブロック図。
【図１３】図１２に対応する復号器側の残差概形計算部
５５の例を示すブロック図。
【図１４】この発明の符号器と復号器の他の実施例を示
すブロック図。
【図１５】図１４における残差概形計算部２３、５５の
具体例を示すブロック図。
【図１６】図１４における残差概形計算部２３、５５の
具体例を示すブロック図。
【図１７】残差概形計算部２３においてスペクトラム概
形の高次帯域部を一定値に近似した処理をする帯域処理
部の構成を示すブロック図。
【図１８】図３の符号器における部分変形例を示すブロ
ック図。
【図１９】この発明による符号化方法及び復号化方法を
適用した符号器及び復号器の他の例を示すブロック図。
【図２０】時間領域で残差信号を得る構成とした符号器
とそれに対応する復号器の例を示すブロック図。
【図２１】図３、１４、１９及び２０の実施例における
量子化部２５の他の構成例を示すブロック図。
【図２２】図２１の量子化部による量子化処理手順を示
すフロー図。

Claims (45)

    【特許請求の範囲】
  1. 【請求項１】 入力音響信号を周波数領域に変換して符
    号化する音響信号変換符号化方法において、 入力音響信号についてフレーム単位に、その周波数特性
    の概形が平坦化された周波数領域の残差係数を得る第１
    の段階と、 上記第１の段階で得られた現フレームの上記残差係数の
    概形を、現在または過去のフレームの上記残差係数から
    予測し、予測残差概形を生成する第２の段階と、 上記第１の段階で得られた現フレームの上記残差係数
    を、上記第２の段階で得られた上記予測残差概形で正規
    化して微細構造係数を得る第３の段階と、 上記第３の段階で得られた上記微細構造係数を量子化
    し、その量子化微細構造を表すインデックスを上記音響
    信号の符号化出力の少なくとも一部として出力する第４
    の段階、とを含む。
  2. 【請求項２】 請求項１記載の符号化方法において、上
    記第２の段階は、現フレームの上記予測残差概形によ
    り、上記量子化微細構造係数を逆正規化して再生残差係
    数を生成する第５の段階と、上記再生残差係数からその
    スペクトラム概形を求める第６段階と、上記スペクトラ
    ム概形から次フレームの上記残差係数に対する上記予測
    残差概形を合成する第７の段階と、を含む。
  3. 【請求項３】 請求項２記載の符号化方法において、上
    記第７の段階は、現フレームに対し予め決めた１つ又は
    連続する複数の過去のフレームの上記再生残差係数のス
    ペクトラム概形から線形合成により上記予測残差概形を
    合成する処理を含む。
  4. 【請求項４】 請求項３記載の符号化方法において、上
    記第２段階は、上記過去のフレームの上記再生残差係数
    のスペクトラム概形から合成する上記予測残差概形が、
    現フレームの残差係数の概形を目標としてそれに近づく
    よう上記過去のフレームの上記スペクトラム概形に対す
    る上記線形合成を制御する第８段階を含む。
  5. 【請求項５】 請求項４記載の符号化方法において、上
    記線形合成の最適制御は、現フレームの上記再生残差係
    数のスペクトラム概形を目標に決定され、決定された最
    適制御は次フレームにおける上記線形合成に適用され
    る。
  6. 【請求項６】 請求項４記載の符号化方法において、上
    記線形合成の最適制御は、現フレームの上記残差係数の
    スペクトラム概形を目標に決定され、決定された最適制
    御は現フレームにおける上記予測残差概形の線形合成に
    適用される。
  7. 【請求項７】 請求項５又は６記載の符号化方法におい
    て、上記第７段階の上記線形合成は、上記過去のフレー
    ムの上記再生残差係数のスペクトラム概形に対してそれ
    ぞれ予測係数を乗算し、それによって得られた乗算結果
    を加算して上記予測残差概形を得る処理であり、上記第
    ８段階は、上記加算結果が上記目標に近づくよう上記予
    測係数を決定する処理を含む。
  8. 【請求項８】 請求項７記載の符号化方法において、上
    記第８段階は上記予測係数決定のための上記目標が現フ
    レームの上記残差係数のスペクトラム概形の場合、上記
    予測係数を量子化したインデックスを上記符号化出力の
    他の一部の符号として出力する第９段階を含む。
  9. 【請求項９】 請求項７又は８記載の符号化方法におい
    て、上記第７段階の上記線形合成は各上記過去のフレー
    ムのサンプル群に対し、周波数軸上において少くとも１
    サンプル正方向及び負方向にそれぞれずれた第１シフト
    サンプル群と第２シフトサンプル群を生成し、上記第１
    及び第２シフトサンプル群にそれぞれ予測係数を乗算
    し、これら乗算結果を上記過去のフレームに対する予測
    係数の乗算結果と共に全て加算することによって上記予
    測残差概形を得る処理を含む。
  10. 【請求項１０】 請求項３、５乃至９のいずれかに記載
    の符号化方法において、上記第６の段階は、上記再生残
    差係数から求めた上記スペクトラム概形の現在及び過去
    の複数フレーム間平均値を各サンプルごとに、または現
    フレーム内サンプルの平均値を求める第１０の段階と、
    上記平均値を現フレームの上記スペクトラム概形から減
    算し、その減算結果を上記スペクトラム概形として上記
    第７段階に与える第１１段階とを含み、上記第７段階は
    上記線形合成の結果に上記平均値を加算し、その加算結
    果から上記予測残差概形を得る第１２の段階を含む。
  11. 【請求項１１】 請求項１０記載の符号化方法におい
    て、上記第６の段階は、上記第１１の段階の減算結果の
    フレーム内平均振幅を計算する第１３の段階と、上記第
    １３の段階で得た上記減算結果の平均振幅で上記第１１
    段階の上記減算結果を割算し、その割算結果を上記スペ
    クトラム概形として上記第７段階に与える第１４段階を
    含み、上記第７段階は上記線形合成の結果に上記第１３
    段階の上記減算結果の平均振幅を乗算し、その乗算結果
    を上記線形合成の結果として上記第１２段階に与える第
    １４段階を含む。
  12. 【請求項１２】 請求項３、５乃至１１のいずれかに記
    載の符号化方法において、上記第６の段階は上記再生残
    差係数の上記スペクトラム概形に窓関数を畳み込む処理
    を含み、上記第７段階はその畳み込み結果を上記スペク
    トラム概形として使い線形合成を行う処理を含む。
  13. 【請求項１３】 請求項３、５乃至１２のいずれかに記
    載の符号化方法において、上記第７の段階は上記線形合
    成の結果に予め決めた定数を加算して上記予測残差概形
    を得る処理を含む。
  14. 【請求項１４】 請求項４乃至９のいずれかに記載の符
    号化方法において、上記第８の段階における上記線形合
    成の制御は、上記目標の周波数領域係数及び上記再生残
    差係数のスペクトラム概形をそれぞれ複数の小帯域に分
    割し、この小帯域ごとに行う処理を含む。
  15. 【請求項１５】 請求項１記載の符号化方法において、
    上記第２の段階は、上記予測残差概形が現フレームにお
    ける上記残差係数のスペクトラム概形に最も近くなるよ
    うに上記スペクトラム概形を量子化し、その量子化を表
    すインデックスを上記符号化出力の他の一部として出力
    する処理を含む。
  16. 【請求項１６】 請求項１５の符号化方法において、上
    記第２段階は現フレームの上記量子化スペクトラム概形
    と、過去のフレームの量子化スペクトラム概形とを予め
    決めた予測係数を使って線形合成し、その線形合成値が
    上記スペクトラム概形と最も近くなるように上記量子化
    スペクトラム概形を決定し、その時の上記線形合成値を
    上記予測残差概形として得る処理を含む。
  17. 【請求項１７】 請求項１５の符号化方法において、上
    記第２段階は現フレームの量子化スペクトラム概形と、
    過去のフレームの上記予測残差概形とを線形合成し、そ
    の線形合成値が上記スペクトラム概形と最も近くなるよ
    うに上記量子化スペクトラム概形を決定し、その時の上
    記線形合成値を上記予測残差概形として得る処理を含
    む。
  18. 【請求項１８】 請求項１乃至１７のいずれかに記載の
    符号化方法において、上記第１段階は、上記入力音響信
    号を周波数領域係数に変換し、上記入力音響信号をフレ
    ーム毎に線形予測分析して線形予測係数を求め、上記線
    形予測係数を周波数領域係数に変換して上記音響信号の
    スペクトラム概形を求め、上記音響信号の周波数領域係
    数を上記スペクトラム概形で正規化して上記残差係数を
    得る処理を含む。
  19. 【請求項１９】 請求項１乃至１７のいずれかに記載の
    符号化方法において、上記第１段階は、上記入力音響信
    号を周波数領域係数に変換し、上記周波数領域係数のス
    ペクトラム振幅を時間領域信号に逆変換し、上記時間領
    域信号を線形予測分析して線形予測係数を求め、上記線
    形予測係数を周波数領域係数に変換して上記音響信号の
    スペクトラム概形を求め、上記音響信号の周波数領域係
    数を上記スペクトラム概形で正規化して上記残差係数を
    得る処理を含む。
  20. 【請求項２０】 請求項１８又は１９記載の符号化方法
    において、上記線形予測係数を周波数領域係数に変換す
    る処理は、上記線形予測係数を量子化して量子化線形予
    測係数を求め、上記量子化線形予測係数を上記線形予測
    係数として上記周波数領域係数に変換すると共に、上記
    量子化線形予測係数を表すインデックスを上記符号化出
    力の他の一部として出力する処理を含む。
  21. 【請求項２１】 請求項１乃至１７のいずれかに記載の
    符号化方法において、上記第１段階は、上記入力音響信
    号を周波数領域係数に変換し、上記周波数領域係数を複
    数の小帯域に分割してそれぞれのスケーリングファクタ
    を求め、上記入力音響信号の周波数領域係数を上記スケ
    ーリングファクタで正規化して上記残差係数を得る処理
    を含む。
  22. 【請求項２２】 請求項１乃至１７のいずれかに記載の
    符号化方法において、上記第１段階は、上記入力音響信
    号を線形予測分析して線形予測係数を求め、上記線形予
    測係数で制御された逆フィルタに上記音響信号を通して
    残差信号を得て、上記残差信号を周波数領域係数に変換
    して上記残差係数を得る処理を含む。
  23. 【請求項２３】 請求項２２記載の符号化方法におい
    て、上記残差信号を得る処理は、上記線形予測係数を量
    子化して得られた量子化線形予測係数を上記線形予測係
    数として上記逆フィルタに与えてそれを制御すると共
    に、上記量子化線形予測係数を表すインデックスを上記
    符号化出力の他の一部として出力する処理を含む。
  24. 【請求項２４】 請求項１乃至２３のいずれかに記載の
    符号化方法において、上記入力音響信号を周波数領域係
    数に変換する処理は、上記入力音響信号をフレーム単位
    に重ね合わせ直交変換する処理を含む。
  25. 【請求項２５】 フレーム毎に予め決めた複数のサンプ
    ルの周波数領域係数に変換されて符号化された音響信号
    を復号する復号化方法において、 入力された第１量子化インデックスから再生した微細構
    造係数を、過去のフレームの情報から予測した残差係数
    の概形で逆正規化して現在フレームにおける再生残差係
    数を得る第１の段階と、 上記第１段階で得られた上記再生残差係数から音響信号
    の周波数特性の概形が与えられた音響信号を再生する第
    ２の段階と、を含む。
  26. 【請求項２６】 請求項２５記載の復号方法において、
    上記第１段階は上記再生残差係数から次フレームに対す
    る上記残差係数の概形を合成する第３段階を含む。
  27. 【請求項２７】 請求項２６記載の復号方法において、
    上記第３の段階は、上記再生残差係数のスペクトラム概
    形を計算する第４の段階と、現フレームに対し、予め決
    めた１又は連続する複数の過去のフレームの上記スペク
    トラム概形にそれぞれ予測係数を乗算して線形合成によ
    り現フレームの上記残差係数の概形を得る第５の段階が
    含まれる。
  28. 【請求項２８】 請求項２７記載の復号方法において、
    上記第５の段階は、上記線形合成により得られる上記残
    差概形が現在フレームにおける上記再生残差係数の概形
    にもっとも近くなるように上記線形合成を適応的に制御
    する第６の段階が含まれる。
  29. 【請求項２９】 請求項２８記載の復号方法において、
    上記第６の段階における上記線形合成の制御は、上記残
    差係数のスペクトラム概形を複数の小帯域に分割し、こ
    の小帯域ごとに行う処理を含む。
  30. 【請求項３０】 請求項２７、２８又は２９記載の復号
    方法において、上記第４の段階は、上記再生残差係数か
    ら求めた上記スペクトラム概形の現在及び過去の複数フ
    レーム間平均値を各サンプルごとに、またはフレーム内
    サンプルの平均値を求める第７の段階と、上記平均値を
    現フレームの上記スペクトラム概形から減算し、その減
    算結果を上記スペクトラム概形として上記第５段階に与
    える第８段階を含み、上記第５段階は上記線形合成の結
    果に上記平均値を加算して上記予測残差概形を得る第９
    の段階を含む。
  31. 【請求項３１】 請求項３０記載の復号方法において、
    上記第３の段階は、上記第８の段階の減算結果のフレー
    ム内平均振幅を計算する第１０の段階と、上記平均振幅
    で上記第８の段階の減算結果を割算し、その割算結果を
    上記スペクトラム概形として上記第５段階に与える第１
    １段階とを含み、上記第５段階は上記線形合成の結果に
    上記減算結果の平均振幅を乗算し、その乗算結果を上記
    線形合成の結果として上記第９段階に与える第１２段階
    を含む。
  32. 【請求項３２】 請求項２７乃至３１のいずれかに記載
    の復号方法において、上記第４の段階は上記再生残差係
    数のスペクトラム概形に窓関数を畳み込む処理を含み、
    上記第５段階はその畳み込み結果を上記スペクトラム概
    形として使って上記線形合成を行う処理を含む。
  33. 【請求項３３】 請求項２７乃至３１のいずれかに記載
    の復号方法において、上記第５段階は上記線形合成は各
    上記過去のフレームのサンプル群に対し、周波数軸上に
    おいて少くとも１サンプル正方向及び負方向にずれた第
    １シフトサンプル群と第２シフトサンプル群を生成し、
    上記第１及び第２シフトサンプル群にそれぞれ予測係数
    を乗算し、これら乗算結果を上記過去のフレームに対す
    る予測係数の乗算結果と共に全て加算することにより上
    記予測残差概形を得る処理を含む。
  34. 【請求項３４】 請求項２７乃至３３のいずれかに記載
    の復号方法において、上記第５の段階は上記線形合成の
    出力に予め決めた定数を加算して上記予測残差概形を得
    る処理を含む。
  35. 【請求項３５】 請求項２６記載の復号方法において、
    上記第３の段階は、上記再生残差係数のスペクトラム概
    形を計算する第４の段階と、現フレームに対し、予め決
    めた１又は連続する複数の過去のフレームの上記スペク
    トラム概形のそれぞれに、入力された第３の量子化イン
    デックスにより指定された上記予測係数を乗算し、これ
    ら乗算結果を加算することにより現フレームの上記再生
    残差係数の概形を得る第５の段階を含む。
  36. 【請求項３６】 請求項２５又は３５記載の復号方法に
    おいて、上記第１段階における上記再生残差概形は、符
    号側から送られたインデックスをそれぞれ逆量子化して
    得た現フレームと過去のフレームの量子化スペクトラム
    概形を線形合成して得る。
  37. 【請求項３７】 請求項２５又は３５記載の復号方法に
    おいて、上記第１段階における上記再生残差概形は、過
    去のフレームにおける合成残差概形と、符号側から送ら
    れたインデックスを逆量子化して得た現フレームの量子
    化スペクトラム概形を線形合成して得る。
  38. 【請求項３８】 請求項２５乃至３７の何れかに記載の
    復号方法において、上記第２段階は、入力された第２量
    子化インデックスを逆量子化して上記音響信号の周波数
    特性の概形情報を再生する処理と、上記周波数特性の概
    形情報に基づいて上記周波数特性の概形が与えられた上
    記音響信号を再生する処理とを含む。
  39. 【請求項３９】 請求項３８記載の復号方法において、
    上記第２段階は、入力された上記第２インデックスから
    上記周波数特性の概形情報として上記音響信号の線形予
    測係数を再生し、上記再生された線形予測係数から上記
    音響信号の周波数特性の概形を求め、上記第１段階で得
    た上記再生残差係数を上記音響信号の周波数特性の概形
    で逆正規化して上記周波数領域係数を得て、上記周波数
    領域係数を時間領域信号に変換して上記音響信号を得る
    処理を含む。
  40. 【請求項４０】 請求項３９記載の復号方法において、
    上記周波数特性の概形を求める処理は、上記線形予測係
    数をフーリエ変換し、その結果得られたスペクトラム振
    幅を上記周波数特性の概形として得る処理を含む。
  41. 【請求項４１】 請求項３８記載の復号方法において、
    上記第２段階は、上記第１段階で得られた上記再生残差
    係数を時間領域の残差信号に変換する処理と、入力され
    た第２量子化インデックスから上記周波数特性の概形情
    報として上記音響信号の線形予測係数を再生する処理
    と、上記線形予測係数をフィルタ係数として上記残差信
    号を逆フィルタ処理して上記音響信号を再生する処理を
    含む。
  42. 【請求項４２】 請求項３８記載の復号方法において、
    上記第２段階は、上記第１段階で得られた上記再生残差
    係数を複数の小帯域に分割し、入力された量子化スケー
    ルファクタインデックスからそれぞれの小帯域に対応す
    るスケールファクタを上記周波数特性の概形情報として
    再生し、それぞれの小帯域の上記再生残差係数を対応す
    る上記スケールファクタで逆正規化して上記周波数特性
    の概形が付与された周波数領域係数を求め、上記周波数
    領域係数を時間領域信号に変換して上記音響信号を再生
    する処理を含む。
  43. 【請求項４３】 請求項３９又は４０記載の復号方法に
    おいて、上記周波数領域係数から上記時間領域信号への
    変換は逆重ね合わせ直交変換である。
  44. 【請求項４４】 請求項３８記載の復号方法において、
    上記第２の段階は上記周波数特性の概形情報に基づいて
    上記再生残差計数に上記周波数特性の概形を与えて周波
    数領域係数を生成する処理と、上記周波数領域係数を時
    間領域信号に変換することにより再生された上記音響信
    号を得る処理を含む。
  45. 【請求項４５】 請求項４４記載の復号方法において、
    上記周波数領域係数から上記時間領域信号への変換は逆
    重ね合わせ直交変換である。
JP07052389A 1994-03-17 1995-03-13 音響信号変換符号化装置および復号化装置 Expired - Lifetime JP3087814B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP07052389A JP3087814B2 (ja) 1994-03-17 1995-03-13 音響信号変換符号化装置および復号化装置

Applications Claiming Priority (7)

Application Number Priority Date Filing Date Title
JP4723594 1994-03-17
JP4844394 1994-03-18
JP11119294 1994-05-25
JP6-48443 1994-05-25
JP6-111192 1994-05-25
JP6-47235 1994-05-25
JP07052389A JP3087814B2 (ja) 1994-03-17 1995-03-13 音響信号変換符号化装置および復号化装置

Publications (2)

Publication Number Publication Date
JPH0844399A true JPH0844399A (ja) 1996-02-16
JP3087814B2 JP3087814B2 (ja) 2000-09-11

Family

ID=27462012

Family Applications (1)

Application Number Title Priority Date Filing Date
JP07052389A Expired - Lifetime JP3087814B2 (ja) 1994-03-17 1995-03-13 音響信号変換符号化装置および復号化装置

Country Status (1)

Country Link
JP (1) JP3087814B2 (ja)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11327600A (ja) * 1997-10-03 1999-11-26 Matsushita Electric Ind Co Ltd オーディオ信号圧縮方法、オーディオ信号圧縮装置、音声信号圧縮方法、音声信号圧縮装置，音声認識方法および音声認識装置
WO2000057399A1 (fr) * 1999-03-19 2000-09-28 Sony Corporation Procede d'inclusion d'information additionnelle et dispositif a cet effet, et procede de decodage d'information additionnelle et dispositif a cet effet
US6345246B1 (en) 1997-02-05 2002-02-05 Nippon Telegraph And Telephone Corporation Apparatus and method for efficiently coding plural channels of an acoustic signal at low bit rates
US6477490B2 (en) 1997-10-03 2002-11-05 Matsushita Electric Industrial Co., Ltd. Audio signal compression method, audio signal compression apparatus, speech signal compression method, speech signal compression apparatus, speech recognition method, and speech recognition apparatus
US6549147B1 (en) 1999-05-21 2003-04-15 Nippon Telegraph And Telephone Corporation Methods, apparatuses and recorded medium for reversible encoding and decoding
JP2007327980A (ja) * 2005-05-11 2007-12-20 Toa Corp デジタルフィルタ、周期性騒音低減装置および騒音低減装置
WO2012005210A1 (ja) * 2010-07-05 2012-01-12 日本電信電話株式会社 符号化方法、復号方法、装置、プログラムおよび記録媒体
WO2012005209A1 (ja) * 2010-07-05 2012-01-12 日本電信電話株式会社 符号化方法、復号方法、装置、プログラム及び記録媒体
WO2012005211A1 (ja) * 2010-07-05 2012-01-12 日本電信電話株式会社 符号化方法、復号方法、符号化装置、復号装置、プログラム、及び記録媒体
JP2013008040A (ja) * 2008-01-24 2013-01-10 Nippon Telegr & Teleph Corp <Ntt> 符号化方法、復号化方法、それらの装置、及びそれらのプログラムと記録媒体
JP5419876B2 (ja) * 2008-08-08 2014-02-19 パナソニック株式会社 スペクトル平滑化装置、符号化装置、復号装置、通信端末装置、基地局装置及びスペクトル平滑化方法
JP2016514857A (ja) * 2013-04-05 2016-05-23 ドルビー・インターナショナル・アーベー オーディオ・エンコーダおよびデコーダ

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6345246B1 (en) 1997-02-05 2002-02-05 Nippon Telegraph And Telephone Corporation Apparatus and method for efficiently coding plural channels of an acoustic signal at low bit rates
US6477490B2 (en) 1997-10-03 2002-11-05 Matsushita Electric Industrial Co., Ltd. Audio signal compression method, audio signal compression apparatus, speech signal compression method, speech signal compression apparatus, speech recognition method, and speech recognition apparatus
JPH11327600A (ja) * 1997-10-03 1999-11-26 Matsushita Electric Ind Co Ltd オーディオ信号圧縮方法、オーディオ信号圧縮装置、音声信号圧縮方法、音声信号圧縮装置，音声認識方法および音声認識装置
WO2000057399A1 (fr) * 1999-03-19 2000-09-28 Sony Corporation Procede d'inclusion d'information additionnelle et dispositif a cet effet, et procede de decodage d'information additionnelle et dispositif a cet effet
US7299189B1 (en) 1999-03-19 2007-11-20 Sony Corporation Additional information embedding method and it's device, and additional information decoding method and its decoding device
US6549147B1 (en) 1999-05-21 2003-04-15 Nippon Telegraph And Telephone Corporation Methods, apparatuses and recorded medium for reversible encoding and decoding
JP2007327980A (ja) * 2005-05-11 2007-12-20 Toa Corp デジタルフィルタ、周期性騒音低減装置および騒音低減装置
JP4690243B2 (ja) * 2005-05-11 2011-06-01 ティーオーエー株式会社 デジタルフィルタ、周期性騒音低減装置および騒音低減装置
JP2013008040A (ja) * 2008-01-24 2013-01-10 Nippon Telegr & Teleph Corp <Ntt> 符号化方法、復号化方法、それらの装置、及びそれらのプログラムと記録媒体
JP5419876B2 (ja) * 2008-08-08 2014-02-19 パナソニック株式会社 スペクトル平滑化装置、符号化装置、復号装置、通信端末装置、基地局装置及びスペクトル平滑化方法
WO2012005211A1 (ja) * 2010-07-05 2012-01-12 日本電信電話株式会社 符号化方法、復号方法、符号化装置、復号装置、プログラム、及び記録媒体
WO2012005209A1 (ja) * 2010-07-05 2012-01-12 日本電信電話株式会社 符号化方法、復号方法、装置、プログラム及び記録媒体
JP5325340B2 (ja) * 2010-07-05 2013-10-23 日本電信電話株式会社 符号化方法、復号方法、符号化装置、復号装置、プログラム、及び記録媒体
JP5331249B2 (ja) * 2010-07-05 2013-10-30 日本電信電話株式会社 符号化方法、復号方法、装置、プログラムおよび記録媒体
JP5331248B2 (ja) * 2010-07-05 2013-10-30 日本電信電話株式会社 符号化方法、復号方法、装置、プログラム及び記録媒体
WO2012005210A1 (ja) * 2010-07-05 2012-01-12 日本電信電話株式会社 符号化方法、復号方法、装置、プログラムおよび記録媒体
US9319645B2 (en) 2010-07-05 2016-04-19 Nippon Telegraph And Telephone Corporation Encoding method, decoding method, encoding device, decoding device, and recording medium for a plurality of samples
JP2016514857A (ja) * 2013-04-05 2016-05-23 ドルビー・インターナショナル・アーベー オーディオ・エンコーダおよびデコーダ
US10043528B2 (en) 2013-04-05 2018-08-07 Dolby International Ab Audio encoder and decoder
US10515647B2 (en) 2013-04-05 2019-12-24 Dolby International Ab Audio processing for voice encoding and decoding
US11621009B2 (en) 2013-04-05 2023-04-04 Dolby International Ab Audio processing for voice encoding and decoding using spectral shaper model

Also Published As

Publication number Publication date
JP3087814B2 (ja) 2000-09-11

Similar Documents

Publication Publication Date Title
US5684920A (en) Acoustic signal transform coding method and decoding method having a high efficiency envelope flattening method therein
US5371853A (en) Method and system for CELP speech coding and codebook for use therewith
US6871106B1 (en) Audio signal coding apparatus, audio signal decoding apparatus, and audio signal coding and decoding apparatus
US5517595A (en) Decomposition in noise and periodic signal waveforms in waveform interpolation
US5950155A (en) Apparatus and method for speech encoding based on short-term prediction valves
EP0770989B1 (en) Speech encoding method and apparatus
EP0858067B1 (en) Multichannel acoustic signal coding and decoding methods and coding and decoding devices using the same
EP0772186A2 (en) Speech encoding method and apparatus
KR19980024885A (ko) 벡터양자화 방법, 음성부호화 방법 및 장치
EP0966793A1 (en) Audio coding method and apparatus
KR19980024519A (ko) 벡터양자화 방법, 음성부호화 방법 및 장치
US6593872B2 (en) Signal processing apparatus and method, signal coding apparatus and method, and signal decoding apparatus and method
JPH10124092A (ja) 音声符号化方法及び装置、並びに可聴信号符号化方法及び装置
JP3541680B2 (ja) 音声音楽信号の符号化装置および復号装置
JP3087814B2 (ja) 音響信号変換符号化装置および復号化装置
JP3186007B2 (ja) 変換符号化方法、復号化方法
JP2003050600A (ja) 線スペクトル平方根を発生し符号化するための方法と装置
JP4359949B2 (ja) 信号符号化装置及び方法、並びに信号復号装置及び方法
JP4281131B2 (ja) 信号符号化装置及び方法、並びに信号復号装置及び方法
JP3237178B2 (ja) 符号化方法及び復号化方法
US6098037A (en) Formant weighted vector quantization of LPC excitation harmonic spectral amplitudes
JP4618823B2 (ja) 信号符号化装置及び方法
JP3099876B2 (ja) 多チャネル音声信号符号化方法及びその復号方法及びそれを使った符号化装置及び復号化装置
JP3163206B2 (ja) 音響信号符号化装置
Ooi et al. A computationally efficient wavelet transform CELP coder

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070714

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20080714

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20080714

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20090714

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20090714

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20100714

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20100714

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20110714

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20120714

Year of fee payment: 12

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

Free format text: PAYMENT UNTIL: 20130714

Year of fee payment: 13

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

EXPY Cancellation because of completion of term