(57)【要約】  音声信号に適した符号化と音楽信号に適した符号化とを階層構造にして組み合わせた符号化方式において、符号化の品質劣化を抑えつつ、符号化装置における処理量を削減することができる符号化装置。この装置において、本選択候補限定部（１０９）は、入力信号のスペクトルと残差スペクトルとを用いて予備選択された所定の数の抑圧係数をCELP成分抑圧部（１０４）に対して指示し、変換符号化部（１１０）は指示された抑圧係数をCELP成分抑圧部（１０４）で用いて生成された抑圧スペクトルをCELP残差信号スペクトル算出部（１０５）に入力して算出された残差スペクトルを用いて第２符号化を行い、歪評価部（１１２）は第２符号化により得られた第２符号を復号して生成された第２復号信号のスペクトルと抑圧スペクトルと入力信号のスペクトルとを用いて、指示された抑圧係数の中から一つの抑圧係数を決定する。
2012032759.tif 000002
【特許請求の範囲】
【請求項１】
入力信号に対する第１の符号化により得られた第１符号を復号して生成された第１復号信号のスペクトルを出力する第１符号化部と、
前記第１復号信号のスペクトルの振幅を、複数の抑圧係数の中から指示された抑圧係数を用いて抑圧して抑圧スペクトルを生成する抑圧部と、
前記入力信号のスペクトルと前記抑圧スペクトルとを用いて残差スペクトルを算出する残差スペクトル算出部と、
前記入力信号のスペクトルと前記残差スペクトルとを用いて、所定の数の抑圧係数を予備選択し、前記予備選択された抑圧係数を前記抑圧部に対して指示する予備選択部と、
前記指示された抑圧係数を前記抑圧部で用いて生成された抑圧スペクトルを前記残差スペクトル算出部に入力して算出された残差スペクトルを用いて第２の符号化を行い、前記第２の符号化により得られた第２符号を復号して生成された第２復号信号のスペクトルと、前記抑圧スペクトルと、前記入力信号のスペクトルと、を用いて、前記指示された抑圧係数の中から一つの抑圧係数を決定する第２符号化部と、
を具備する符号化装置。
【請求項２】
前記第２符号化部は、
前記残差スペクトルに対して立てられたパルスを前記第２の符号化で符号化し、前記第２の符号化による符号化歪が最小となる前記抑圧係数を探索し、
前記予備選択部は、
前記残差スペクトルを用いて前記パルスの位置を推定する推定手段と、
前記残差スペクトルのうち、推定された前記パルスの位置における振幅を減衰させて推定残差スペクトルを生成する減衰手段と、
前記推定残差スペクトルと前記入力信号のスペクトルとを用いて、前記符号化歪の推定エネルギである推定歪エネルギを算出する算出手段と、
前記推定歪エネルギに基づいて、前記複数の抑圧係数のうち、前記所定の数の抑圧係数を予備選択する候補限定手段と、
を具備する請求項１記載の符号化装置。
【請求項３】
前記複数の抑圧係数は、抑圧の程度に関して昇順又は降順でインデックスが付与され、
前記候補限定手段は、
最大のインデックス及び最小のインデックスに対応する前記抑圧係数のうち、前記推定歪エネルギが大きい方の抑圧係数を、前記所定の数の抑圧係数から除外する、
請求項２記載の符号化装置。
【請求項４】
前記複数の抑圧係数は、抑圧の程度に関して昇順又は降順でインデックスが付与され、
前記候補限定手段は、
前記複数の抑圧係数のうちの前記推定歪エネルギが最も小さい抑圧係数と、前記推定歪エネルギが最も小さい抑圧係数に付与されたインデックスの前後のインデックスに対応する２つの抑圧係数と、を前記所定の数の抑圧係数として予備選択する、
請求項２記載の符号化装置。
【請求項５】
前記複数の抑圧係数は、抑圧の程度に関して昇順又は降順でインデックスが付与され、
前記候補限定手段は、
前記複数の抑圧係数のうちの前記推定歪エネルギが最も小さい第１の抑圧係数と、前記第１の抑圧係数に付与されたインデックスの前後のインデックスに対応する２つの抑圧係数のうち前記推定歪エネルギが小さい方の第２の抑圧係数と、を前記所定の数の抑圧係数として予備選択する、
請求項２記載の符号化装置。
【請求項６】
前記推定手段は、
前記残差スペクトルの振幅の統計量に基づいて算出される閾値と、前記残差スペクトルの振幅とを比較して、前記パルスの位置を推定する、
請求項２記載の符号化装置。
【請求項７】
前記統計量は、少なくとも前記振幅の標準偏差を含む、
請求項６記載の符号化装置。
【請求項８】
前記減衰手段は、
推定された前記パルスの位置におけるスペクトルの振幅に、０以上１未満の値を有する係数を乗算して前記振幅を減衰させる、
請求項２記載の符号化装置。
【請求項９】
前記推定手段は、
前記残差スペクトルの特徴に応じて、推定する前記パルスの本数を設定し、設定された本数の前記パルスの位置を推定する、
請求項２記載の符号化装置。
【請求項１０】
前記特徴は、前記残差スペクトルの各帯域における振幅のばらつきであり、
前記推定手段は、
前記パルスの本数を、前記ばらつきが大きくなるほど少なくするように設定する、
請求項９記載の符号化装置。
【請求項１１】
前記特徴は、前記残差スペクトルのトーン性であり、
前記推定手段は、
前記パルスの本数を、前記トーン性が高くなるほど少なくするように設定する、
請求項９記載の符号化装置。
【請求項１２】
前記特徴は、前記残差スペクトルの雑音性であり、
前記推定手段は、
前記パルスの本数を、前記雑音性が高くなるほど多くするように設定する、
請求項９記載の符号化装置。
【請求項１３】
前記減衰手段は、
前記残差スペクトルの特徴に応じて、推定された前記パルスの位置におけるスペクトルの振幅を減衰させる、
請求項２記載の符号化装置。
【請求項１４】
前記特徴は、前記残差スペクトルの各帯域における振幅のばらつきであり、
前記減衰手段は、
前記スペクトルの減衰度合を、前記ばらつきが大きくなるほど大きくするように設定する、
請求項１３記載の符号化装置。
【請求項１５】
前記特徴は、前記残差スペクトルのトーン性であり、
前記減衰手段は、
前記スペクトルの減衰度合を、前記トーン性が高くなるほど大きくするように設定する、
請求項１３記載の符号化装置。
【請求項１６】
前記特徴は、前記残差スペクトルの雑音性であり、
前記減衰手段は、
前記スペクトルの減衰度合を、前記雑音性が高くなるほど小さくするように設定する、
請求項１３記載の符号化装置。
【請求項１７】
入力信号に対する第１の符号化により得られた第１符号を復号して生成された第１復号信号のスペクトルを出力する第１符号化ステップと、
前記第１復号信号のスペクトルの振幅を、複数の抑圧係数の中から指示された抑圧係数を用いて抑圧して抑圧スペクトルを生成する抑圧ステップと、
前記入力信号のスペクトルと前記抑圧スペクトルとを用いて残差スペクトルを算出する残差スペクトル算出ステップと、
前記入力信号のスペクトルと前記残差スペクトルとを用いて、前記抑圧ステップで用いる所定の数の抑圧係数を予備選択し、前記予備選択された抑圧係数を前記指示された抑圧係数に設定する予備選択ステップと、
前記指示された抑圧係数を前記抑圧ステップで用いて生成された抑圧スペクトルを用いて前記残差スペクトル算出ステップで算出された残差スペクトルを用いて第２の符号化を行い、前記第２の符号化により得られた第２符号を復号して生成された第２復号信号のスペクトルと、前記抑圧スペクトルと、前記入力信号のスペクトルと、を用いて、前記指示された抑圧係数の中から一つの抑圧係数を決定する第２符号化ステップと、
を有する符号化方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、符号化装置及び符号化方法に関する。
【背景技術】
【０００２】
音声及び音楽等を、低ビットレートかつ高音質で圧縮できる符号化方式として、音声信号に適したCELP（Code Excited Linear Prediction）符号化方式と、音楽信号に適した変換符号化方式とを階層構造にして組み合わせた符号化方式が提案されている（例えば、非特許文献１参照）。なお、以下においては、音声信号と音楽信号とを総称して音響信号と称することもある。【０００３】
この符号化方式では、符号化装置は、まず、CELP符号化方式で入力信号を符号化してCELP符号化データを生成する。次いで、符号化装置は、入力信号とCELP復号信号（CELP符号化データの復号結果）との残差信号（以下、CELP残差信号と呼ぶ）を周波数領域に変換して得られる残差スペクトルを変換符号化することにより、高音質化を図っている。変換符号化方式としては、残差スペクトルのエネルギが大きい周波数にパルスを立てて、そのパルスの情報を符号化する方式が提案されている（非特許文献１参照）。【０００４】
しかしながら、CELP符号化方式は、音声信号の符号化には適しているが、音楽信号に対しては符号化モデルが異なるので音質が悪くなる。そのため、上記符号化方式で音楽信号を符号化した場合、CELP残差信号の成分が大きくなるので、変換符号化によりCELP残差信号（残差スペクトル）を符号化しても音質が向上しにくいという課題がある。【０００５】
この課題を解決するために、CELP復号信号の周波数成分（以下、CELP成分と呼ぶ）の振幅を抑圧した結果を用いて算出される残差スペクトルを変換符号化することで高音質化を図る符号化方式（CELP成分抑圧方法）が提案されている（例えば、特許文献１及び非特許文献１（section 6.11.6.2）参照）。【０００６】
非特許文献１に開示されたCELP成分抑圧方法では、入力信号のサンプリング周波数が16kHzの場合、0.8kHz〜5.5kHzの中帯域のみでCELP成分の振幅の抑圧（以下、CELP抑圧と呼ぶ）が行われる。ただし、非特許文献１では、符号化装置は、CELP残差信号に対して変換符号化を直接行うのではなく、その前に別の変換符号化方式（例えば、非特許文献１（Section 6.11.6.1）参照）によってCELP成分の残差信号を小さくしている。このため、符号化装置は、中帯域であっても上記別の変換符号化方式によって符号化された周波数成分に対してはCELP抑圧を行わない。また、中帯域内のCELP抑圧を行わない周波数以外の他の周波数では、CELP抑圧の程度（強さ）を示すCELP抑圧係数は一様である。CELP抑圧係数は、CELP抑圧の強度別にコードブック（以下、CELP抑圧係数コードブックと呼ぶ）に格納されている。CELP抑圧係数コードブックには、CELP成分を全く抑圧しないことを意味する係数（＝1.0）も格納されている。【０００７】
符号化装置は、変換符号化を行う前に、CELP成分（CELP復号信号）と、CELP抑圧係数コードブックに格納されているCELP抑圧係数とを乗じることでCELP抑圧を行ってから、入力信号とCELP復号信号（CELP抑圧後のCELP復号信号）との残差スペクトルを求め、残差スペクトルを変換符号化する。この変換符号化は、全てのCELP抑圧係数に対して行われる。そして、符号化装置は、変換符号化データの復号信号とCELP成分が抑圧されたCELP復号信号とを加算した信号と、入力信号との残差信号を算出し、残差信号のエネルギ（以下、符号化歪と呼ぶ）が最小となるCELP抑圧係数を決定して、探索したCELP抑圧係数（符号化歪が最小となるCELP抑圧係数）を符号化する。これにより、符号化装置では、帯域全体として符号化歪を最小にした変換符号化を行うことができる。以下では、CELP抑圧係数毎に変換符号化を行い、符号化歪（残差信号のエネルギ）が最小となるCELP抑圧係数を決定する一連の処理を「本選択」と呼ぶこととする。【０００８】
一方、復号装置は、符号化装置から送信されるCELP抑圧係数を用いて、CELP復号信号のCELP成分を抑圧し、CELP成分が抑圧されたCELP復号信号に変換符号化の復号信号を加算する。これにより、復号装置では、CELP符号化と変換符号化とを階層構造にして組み合わせた符号化を行う際のCELP符号化による音質の劣化を抑えた復号信号を得ることができる。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】米国特許出願公開第２００９／０１１２６０７号明細書
【非特許文献】
【００１０】
【非特許文献１】Recommendation ITU-T G.718,2008年6月
【発明の概要】
【発明が解決しようとする課題】
【００１１】
しかしながら、上述したCELP成分抑圧方法により、CELP抑圧係数コードブックに格納されているCELP抑圧係数毎に変換符号化を行うことで、符号化歪の評価（以下、歪評価と呼ぶことがある）を行う場合には、CELP抑圧係数の全ての候補、つまり、CELP抑圧係数コードブックに格納されている全てのCELP抑圧係数に対して変換符号化を行う必要があるため、符号化装置における処理量が非常に大きくなってしまうという課題がある。【００１２】
本発明の目的は、CELP抑圧係数毎に生成される、変換符号化処理に対する入力信号（以下、ターゲット信号と呼ぶ）の中から一部を選択（以下、「予備選択」と呼ぶ）して、本選択において変換符号化を行う対象を限定することで、符号化の品質劣化を抑えつつ、符号化装置における処理量を削減することができる符号化装置及び符号化方法を提供することである。
【課題を解決するための手段】
【００１３】
本発明の一態様に係る符号化装置は、入力信号に対する第１の符号化により得られた第１符号を復号して生成された第１復号信号のスペクトルを出力する第１符号化部と、前記第１復号信号のスペクトルの振幅を、複数の抑圧係数の中から指示された抑圧係数を用いて抑圧して抑圧スペクトルを生成する抑圧部と、前記入力信号のスペクトルと前記抑圧スペクトルとを用いて残差スペクトルを算出する残差スペクトル算出部と、前記入力信号のスペクトルと前記残差スペクトルとを用いて、所定の数の抑圧係数を予備選択し、前記予備選択された抑圧係数を前記抑圧部に対して指示する予備選択部と、前記指示された抑圧係数を前記抑圧部で用いて生成された抑圧スペクトルを前記残差スペクトル算出部に入力して算出された残差スペクトルを用いて第２の符号化を行い、前記第２の符号化により得られた第２符号を復号して生成された第２復号信号のスペクトルと、前記抑圧スペクトルと、前記入力信号のスペクトルと、を用いて、前記指示された抑圧係数の中から一つの抑圧係数を決定する第２符号化部と、を具備する。【００１４】
本発明の一態様に係る符号化方法は、入力信号に対する第１の符号化により得られた第１符号を復号して生成された第１復号信号のスペクトルを出力する第１符号化ステップと、前記第１復号信号のスペクトルの振幅を、複数の抑圧係数の中から指示された抑圧係数を用いて抑圧して抑圧スペクトルを生成する抑圧ステップと、前記入力信号のスペクトルと前記抑圧スペクトルとを用いて残差スペクトルを算出する残差スペクトル算出ステップと、前記入力信号のスペクトルと前記残差スペクトルとを用いて、前記抑圧ステップで用いる所定の数の抑圧係数を予備選択し、前記予備選択された抑圧係数を前記指示された抑圧係数に設定する予備選択ステップと、前記指示された抑圧係数を前記抑圧ステップで用いて生成された抑圧スペクトルを用いて前記残差スペクトル算出ステップで算出された残差スペクトルを用いて第２の符号化を行い、前記第２の符号化により得られた第２符号を復号して生成された第２復号信号のスペクトルと、前記抑圧スペクトルと、前記入力信号のスペクトルと、を用いて、前記指示された抑圧係数の中から一つの抑圧係数を決定する第２符号化ステップと、を有する。
【発明の効果】
【００１５】
本発明によれば、音声信号に適した符号化と音楽信号に適した符号化とを階層構造にして組み合わせた符号化方式において、全てのCELP抑圧係数候補に対して変換符号化を逐次行う方法と比較して、符号化の品質劣化を抑えつつ、符号化装置における処理量を削減することができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の実施の形態１に係る符号化装置の構成を示すブロック図
【図２】本発明の実施の形態１に係る復号装置の構成を示すブロック図
【図３】本発明の実施の形態２に係る符号化装置の構成を示すブロック図
【発明を実施するための形態】
【００１７】
以下、本発明の各実施の形態について、図面を参照して詳細に説明する。なお、本発明に係る符号化装置及び復号装置として、音響符号化装置及び音響復号装置を例にとって説明する。なお、上述のように、音声信号と音楽信号とを総称して音響信号と称することとする。すなわち、音響信号は、実質的に音声信号のみ、実質的に音楽信号のみ、音声信号及び音楽信号が混在した信号、のいずれの信号をも表すものとする。
【００１８】
また、本発明に係る符号化装置及び復号装置は、少なくとも２つの符号化を行う階層を有する。以下の説明においては、音声信号に適した符号化としてCELP符号化を、音楽信号に適した符号化として変換符号化を、それぞれ代表して用いるものとし、符号化装置及び復号装置は、CELP符号化と変換符号化とを階層構造にして組み合わせた符号化方式を用いる。
【００１９】
（実施の形態１）
図１は、本発明の実施の形態１に係る符号化装置１００の主要な構成を示すブロック図である。符号化装置１００は、音声及び音楽等の入力信号を、CELP符号化と変換符号化とを階層構造にして組み合わせた符号化方式を用いて符号化して、符号化データを出力する。図１に示すように、符号化装置１００は、MDCT（Modified Discrete Cosine Transform：修正離散コサイン変換）部１０１、CELP符号化部１０２、MDCT部１０３、CELP成分抑圧部１０４、CELP残差信号スペクトル算出部１０５、パルス位置推定部１０６、推定パルス減衰部１０７、推定歪評価部１０８、本選択候補限定部１０９、変換符号化部１１０、加算部１１１、歪評価部１１２及び多重化部１１３を具備する。各部は以下の動作を行う。
【００２０】
図１に示す符号化装置１００において、MDCT部１０１は、入力信号に対してMDCT処理を行って入力信号スペクトルを生成する。そして、MDCT部１０１は、生成した入力信号スペクトルをCELP残差信号スペクトル算出部１０５、歪評価部１１２及び推定歪評価部１０８に出力する。
【００２１】
CELP符号化部１０２は、入力信号をCELP符号化方式により符号化してCELP符号化データを生成する。また、CELP符号化部１０２は、生成したCELP符号化データを復号（ローカルデコード）してCELP復号信号を生成する。そして、CELP符号化部１０２は、CELP符号化データを多重化部１１３に出力し、CELP復号信号をMDCT部１０３に出力する。
【００２２】
MDCT部１０３は、CELP符号化部１０２から入力されるCELP復号信号に対してMDCT処理を行ってCELP復号信号スペクトルを生成する。そして、MDCT部１０３は、生成したCELP復号信号スペクトルをCELP成分抑圧部１０４に出力する。
【００２３】
このように、例えば、CELP符号化部１０２及びMDCT部１０３は、入力信号に対する第１の符号化により得られた第１符号を復号して生成された第１復号信号のスペクトルを出力する第１符号化部として動作する。
【００２４】
CELP成分抑圧部１０４は、CELP抑圧の程度（強さ）を示すCELP抑圧係数が格納されたCELP抑圧係数コードブックを具備する。例えば、CELP抑圧係数コードブックには、抑圧しないことを意味する1.0から、CELP成分の振幅を半分にする0.5までの４種類のCELP抑圧係数が格納されている。つまり、CELP抑圧係数は、CELP抑圧の程度（強さ）が大きいほど値がより小さくなる。また、ここでのCELP抑圧係数コードブックでは、CELP抑圧係数がCELP抑圧の程度（強さ）の昇順又は降順で格納されているものとする。また、各CELP抑圧係数には、CELP抑圧の程度（強さ）に関して昇順又は降順でインデックス（CELP抑圧係数インデックス）が付与されているものとする。
【００２５】
まず、CELP成分抑圧部１０４は、推定歪評価部１０８、本選択候補限定部１０９又は歪評価部１１２から入力されるCELP抑圧係数インデックスに従って、CELP抑圧係数コードブックの中からCELP抑圧係数を選択する。そして、CELP成分抑圧部１０４は、選択したCELP抑圧係数を、MDCT部１０３から入力されるCELP復号信号スペクトルの周波数成分毎に乗じて、CELP成分抑圧スペクトルを算出する。そして、CELP成分抑圧部１０４は、CELP成分抑圧スペクトルをCELP残差信号スペクトル算出部１０５及び加算部１１１に出力する。
【００２６】
CELP残差信号スペクトル算出部１０５は、MDCT部１０１から入力される入力信号スペクトルと、CELP成分抑圧部１０４から入力されるCELP成分抑圧スペクトルとの差分であるCELP残差信号スペクトルを算出する。具体的には、CELP残差信号スペクトル算出部１０５は、入力信号スペクトルからCELP成分抑圧スペクトルを減じることで、CELP残差信号スペクトルを得る。そして、CELP残差信号スペクトル算出部１０５は、CELP残差信号スペクトルを変換符号化部１１０、パルス位置推定部１０６及び推定パルス減衰部１０７に出力する。
【００２７】
パルス位置推定部１０６は、CELP残差信号スペクトル算出部１０５から入力されるCELP残差信号スペクトル（変換符号化対象の信号。以下、ターゲット信号と呼ぶことがある。）を用いて、変換符号化部１１０で符号化されるパルス位置（例えば、CELP残差信号スペクトルの振幅が大きい周波数）を推定する。そして、パルス位置推定部１０６は、推定したパルス位置（推定パルス位置）を推定パルス減衰部１０７に出力する。
【００２８】
推定パルス減衰部１０７は、CELP残差信号スペクトル算出部１０５から入力されるCELP残差信号スペクトルのうち、パルス位置推定部１０６から入力される推定パルス位置における振幅を減衰させる。そして、推定パルス減衰部１０７は、減衰後のスペクトルを変換符号化推定残差スペクトルとして推定歪評価部１０８に出力する。
【００２９】
推定歪評価部１０８は、MDCT部１０１から入力される入力信号スペクトル、及び、推定パルス減衰部１０７から入力される変換符号化推定残差スペクトルを用いて、変換符号化による符号化歪（歪エネルギ）の推定値である推定歪エネルギを算出する。そして、推定歪評価部１０８は、推定歪エネルギを本選択候補限定部１０９に出力する。
【００３０】
また、推定歪評価部１０８は、後述する予備選択探索において評価対象のCELP抑圧係数に対応する変換符号化推定残差スペクトルを得るために、評価対象のCELP抑圧係数インデックスをCELP成分抑圧部１０４に出力する。例えば、推定歪評価部１０８は、CELP抑圧係数インデックスｊ＝１に対する推定歪エネルギを算出する際には、CELP抑圧係数インデックスｊ＝１をCELP成分抑圧部１０４に出力する。そして、推定歪評価部１０８は、CELP成分抑圧部１０４、CELP残差信号スペクトル算出部１０５、パルス位置推定部１０６、推定パルス減衰部１０７で順次処理された結果である変換符号化推定残差スペクトル（CELP抑圧係数インデックスｊ＝１に対応）に対する推定歪エネルギを算出する。
【００３１】
本選択候補限定部１０９は、推定歪評価部１０８から入力される推定歪エネルギの分布に基づいて、CELP抑圧コードブックに格納されているCELP抑圧係数のうち、後述する本選択探索で探索するCELP抑圧係数（変換符号化に用いるCELP抑圧係数）の候補を限定する。そして、本選択候補限定部１０９は、限定されたCELP抑圧係数の候補を示すCELP抑圧係数インデックスをCELP成分抑圧部１０４に出力する。なお、以下において、ここで限定されたCELP抑圧係数の候補をまとめてCELP抑圧係数群、また、限定されたCELP抑圧係数の候補に対応するCELP抑圧係数インデックスをまとめてCELP抑圧係数インデックス群、と呼ぶことがある。
【００３２】
このように、例えば、パルス位置推定部１０６、推定パルス減衰部１０７、推定歪評価部１０８及び本選択候補限定部１０９は、入力信号スペクトルとCELP残差信号スペクトルとを用いて、所定の数のCELP抑圧係数を予備選択し、予備選択されたCELP抑圧係数をCELP成分抑圧部１０４に対して指示する予備選択部として動作する。
【００３３】
なお、図１に示す符号化装置１００において、CELP成分抑圧部１０４、CELP残差信号スペクトル算出部１０５、パルス位置推定部１０６、推定パルス減衰部１０７、推定歪評価部１０８及び本選択候補限定部１０９は、閉ループを構成する。この閉ループを構成する各構成部は、CELP成分抑圧部１０４が具備するCELP抑圧コードブックに格納されているCELP抑圧係数のうち、推定歪評価部１０８から指示されるCELP抑圧係数インデックスに対応するCELP抑圧係数を用いて、後述する本選択探索において探索対象となる候補（CELP抑圧係数インデックス）を探索する。以下、この探索処理を、「予備選択探索」と呼ぶ。
【００３４】
変換符号化部１１０は、CELP残差信号スペクトル算出部１０５から入力されるCELP残差信号スペクトル（ターゲット信号）を変換符号化により符号化して、変換符号化データを生成する。また、変換符号化部１１０は、生成した変換符号化データを復号（ローカルデコード）して、変換符号化復号信号スペクトルを生成する。このとき、変換符号化部１１０は、CELP残差信号スペクトルと変換符号化復号信号スペクトルとの歪が小さくなるように符号化を行う。例えば、変換符号化部１１０は、CELP残差信号スペクトルの振幅（エネルギ）が大きい周波数に、パルスを立てることで上記歪を小さくするように符号化を行う。そして、変換符号化部１１０は、符号化により得られた変換符号化データを歪評価部１１２に出力し、変換符号化復号信号スペクトルを加算部１１１に出力する。
【００３５】
加算部１１１は、CELP成分抑圧部１０４から入力されるCELP成分抑圧スペクトルと、変換符号化部１１０から入力される変換符号化復号信号スペクトルとを加算して復号信号スペクトルを算出し、復号信号スペクトルを歪評価部１１２に出力する。
【００３６】
歪評価部１１２は、CELP成分抑圧部１０４が備えるCELP抑圧係数コードブックに格納されたCLEP抑圧係数のうち、一部のインデックス（本選択候補限定部１０９で限定されたCELP抑圧係数インデックス）を走査して、MDCT部１０１から入力される入力信号スペクトルと加算部１１１から入力される復号信号スペクトルとの歪（すなわち、変換符号化による符号化歪）が最小となるCELP抑圧係数インデックスを探索する。つまり、歪評価部１１２は、上記一部のインデックスに対応するCELP抑圧係数を用いてCELP抑圧を行うようにCELP成分抑圧部１０４を制御する（CELP抑圧係数インデックスを出力する）。そして、歪評価部１１２は、算出した歪が最小となるCELP抑圧係数インデックスを、CELP抑圧係数最適インデックスとして多重化部１１３に出力し、変換符号化部１１０から入力される変換符号化データのうちCELP抑圧係数最適インデックスに対応する変換符号化データ（歪最小時の変換符号化データ）を多重化部１１３に出力する。
【００３７】
このように、例えば、変換符号化部１１０、加算部１１１及び歪評価部１１２は、上述した予備選択部から指示されたCELP抑圧係数をCELP成分抑圧部１０４で用いて生成されたCELP抑圧スペクトルをCELP残差信号スペクトル算出部１０５に入力して算出されたCELP残差信号スペクトルを用いて変換符号化（第２の符号化）を行い、変換符号化により得られた変換符号化データ（第２符号）を復号して生成された変換符号化復号信号スペクトル（第２復号信号のスペクトル）と、CELP抑圧スペクトルと、入力信号スペクトルと、を用いて、指示されたCELP抑圧係数の中から一つのCELP抑圧係数を決定する第２符号化部として動作する。
【００３８】
なお、図１に示す符号化装置１００において、CELP成分抑圧部１０４、CELP残差信号スペクトル算出部１０５、変換符号化部１１０、加算部１１１及び歪評価部１１２は、閉ループを構成する。この閉ループを構成する各構成部は、CELP成分抑圧部１０４が具備するCELP抑圧コードブックに格納されている複数のCELP抑圧係数のうち、本選択候補限定部１０９から指示されるCELP抑圧係数インデックスに対応するCELP抑圧係数を用いて復号信号スペクトルを生成し、入力信号スペクトルと復号信号スペクトルとの歪（変換符号化による符号化歪）が最小となる候補（CELP抑圧係数インデックス）を探索する。以下、この探索処理を、「本選択探索」と呼ぶ。
【００３９】
多重化部１１３は、CELP符号化部１０２から入力されるCELP符号化データ、歪評価部１１２から入力される変換符号化データ（歪最小時の変換符号化データ）及びCELP抑圧係数最適インデックスを多重化して、多重化結果を符号化データとして復号装置へ送信する。
【００４０】
次に、復号装置２００について説明する。復号装置２００は、符号化装置１００から送信される符号化データを復号して、復号信号を出力する。
【００４１】
図２は、復号装置２００の主要な構成を示すブロック図である。復号装置２００は、分離部２０１、変換符号化復号部２０２、CELP復号部２０３、MDCT部２０４、CELP成分抑圧部２０５、加算部２０６、IMDCT（Inverse Modified Discrete Cosine Transform：逆修正離散コサイン変換）部２０７を具備する。各部は以下の動作を行う。
【００４２】
図２に示す復号装置２００において、分離部２０１は、CELP符号化データと、変換符号化データと、CELP抑圧係数最適インデックスとを含む符号化データを、符号化装置１００（図１）から伝送路（図示せず）を介して受信する。分離部２０１は、符号化データを、CELP符号化データと、変換符号化データと、CELP抑圧係数最適インデックスとに分離する。そして、分離部２０１は、CELP符号化データをCELP復号部２０３に出力し、変換符号化データを変換符号化復号部２０２に出力し、CELP抑圧係数最適インデックスをCELP成分抑圧部２０５に出力する。
【００４３】
変換符号化復号部２０２は、分離部２０１から入力される変換符号化データを復号して、変換符号化復号信号スペクトルを生成し、変換符号化復号信号スペクトルを加算部２０６に出力する。
【００４４】
CELP復号部２０３は、分離部２０１から入力されるCELP符号化データを復号して、CELP復号信号をMDCT部２０４に出力する。
【００４５】
MDCT部２０４は、CELP復号部２０３から入力されるCELP復号信号に対して、MDCT処理を行ってCELP復号信号スペクトルを生成する。そして、MDCT部２０４は、生成したCELP復号信号スペクトルをCELP成分抑圧部２０５に出力する。
【００４６】
CELP成分抑圧部２０５は、CELP成分抑圧部１０４が具備するCELP抑圧係数コードブックと同様のCELP抑圧係数コードブックを具備する。CELP成分抑圧部２０５が具備するCELP抑圧係数コードブックは、基本的にはCELP成分抑圧部１０４が具備するCELP抑圧係数コードブックと全く同じCELP抑圧係数コードブックであればよいが、他の何らかの調整等も含めて抑圧する場合には、必ずしも同じでなくてもよい。CELP成分抑圧部２０５は、分離部２０１から入力されるCELP抑圧係数最適インデックスに対応するCELP抑圧係数を、MDCT部２０４から入力されるCELP復号信号スペクトルの周波数成分毎に乗ずることにより、CELP復号信号スペクトル（CELP成分）が抑圧されたCELP成分抑圧スペクトルを算出する。そして、CELP成分抑圧部２０５は、算出したCELP成分抑圧スペクトルを加算部２０６に出力する。
【００４７】
加算部２０６は、符号化装置１００の加算部１１１と同様にして、CELP成分抑圧部２０５から入力されるCELP成分抑圧スペクトルと、変換符号化復号部２０２から入力される変換符号化復号信号スペクトルとを加算して、復号信号スペクトルを算出する。そして、加算部２０６は、算出した復号信号スペクトルをIMDCT部２０７に出力する。
【００４８】
IMDCT部２０７は、加算部２０６から入力される復号信号スペクトルに対して、IMDCT処理を行って復号信号を出力する。
【００４９】
次に、符号化装置１００（図１）における予備選択探索処理の詳細について説明する。
【００５０】
まず、パルス位置推定部１０６における、推定パルス位置の推定方法の一例について説明する。
【００５１】
一般に、変換符号化では、入力信号（ここでは、CELP残差信号スペクトル）の振幅が大きい周波数にパルスを立てるように符号化が行われる。このとき、立てられるパルスの本数、及び、パルスの振幅と入力信号との誤差は、設定されたビットレート又は信号の周波数特性により異なる。そのため、変換符号化における符号化歪は実際に符号化を行わないと正確に求めることができない。ただし、変換符号化において符号化されるパルス位置は、統計的手法を用いることにより推定することが可能である。
【００５２】
ここで、CELP残差信号スペクトルが正規分布であると仮定する。また、変換符号化では振幅がより大きい周波数でパルスが立ち、パルスの情報が符号化されるとする。例えば、符号化装置１００は、CELP残差信号スペクトルのうち、振幅が大きい上位１０％の周波数でパルスが符号化されると仮定して、変換符号化部１１０で符号化されるパルス位置を判定するための閾値（振幅の閾値）を算出する。
【００５３】
具体的には、まず、CELP残差信号スペクトルの絶対値平均Iavg[j]が、次式（１）に従って算出される。
【数１】
2012032759.tif 000003
【００５４】
ここで、Iavg[j]はCELP抑圧係数インデックスｊにおけるCELP残差信号スペクトルの絶対値平均を表し、ｉは周波数サンプルの番号を表し、CrはCELP残差信号スペクトルの振幅を表す。また、CELP抑圧係数インデックスの総数をＭ個とし、周波数サンプルの総数をＮ個とする。
【００５５】
次いで、CELP抑圧係数インデックスｊにおけるCELP残差信号スペクトルの標準偏差σ[j]が、次式（２）に従って算出される。
【数２】
2012032759.tif 000004
【００５６】
そして、式（１）により算出された絶対値平均Iavg[j]及び式（２）により算出された標準偏差σ[j]を用いて閾値Ithrは、例えば、次式（３）に従って算出される。
【数３】
2012032759.tif 000005
【００５７】
ここで、βは閾値Ithrの値を制御する定数である。例えば、CELP残差信号スペクトルのうち、振幅が大きい上位１０％の周波数が選択されるように閾値を設定する際には、βの値を約１．６に設定する。また、例えば、CELP残差信号スペクトルのうち、振幅が大きい上位５％の周波数が選択されるように閾値を設定する際には、βの値を約２．０に設定する。なお、βの設定値は正規分布表に従って求めることができる。
【００５８】
パルス位置推定部１０６は、式（３）に示す閾値Ithrを用いることで、変換符号化部１１０で符号化されるパルス位置（推定パルス位置）を推定する。具体的には、パルス位置推定部１０６は、次式（４）に従って、CELP抑圧係数インデックスｊにおいて、変換符号化部１１０で符号化されるパルス位置を推定する。
【数４】
2012032759.tif 000006
【００５９】
ここで、Iep[j][i]は、CELP抑圧係数インデックスｊの各周波数サンプルｉ（１≦ｉ≦Ｎ）においてパルスが立てられるか否かの推定結果を示す。すなわち、式（４）に示すように、CELP抑圧係数インデックスｊにおいて、パルスが立てられると推定された周波数サンプルｉではIep[j][i]=1.0となり、それ以外の周波数サンプルではIep[j][i]=0.0となる。すなわち、パルス位置推定部１０６は、Iep[j][i]=1.0となる周波数サンプルを、推定パルス位置とする。
【００６０】
このように、パルス位置推定部１０６は、CELP残差信号スペクトル(ターゲット信号)の分布特性に基づき、変換符号化部１１０での符号化の結果として求められるパルスの位置を低演算量で効率的に推定している。具体的には、パルス位置推定部１０６は、CELP残差信号スペクトル（ターゲット信号）の振幅又は絶対値の統計量に基づいて算出される閾値（Ithr）と、CELP残差信号スペクトルの振幅とを比較して、変換符号化部１１０で符号化されるパルス（推定パルス位置）を推定する。これにより、パルス位置推定部１０６では、振幅の閾値判定を行うのみでよく、変換符号化部１１０で符号化されると推定されるパルス位置を、変換符号化部１１０での処理量よりも少ない処理量で特定することが可能となる。また、パルス位置推定部１０６で用いられる上記統計量として、標準偏差σを少なくとも含むようにすればよい。このようにターゲット信号の振幅又は絶対値のばらつきの度合いを定量的に表す標準偏差を用いて閾値を算出することにより、少ない演算量でパルス位置の推定精度の高い閾値を算出することが可能となる。
【００６１】
次いで、推定パルス減衰部１０７は、パルス位置推定部１０６で推定された推定パルス位置（Iep[j][i]=1.0に対応する帯域）の振幅を減衰させて、変換符号化推定残差スペクトルを生成する。
【００６２】
例えば、ここでは、簡単のため、推定パルス減衰部１０７でのスペクトル減衰の結果、推定パルス位置（Iep[j][i]=1.0に対応する帯域）では、CELP残差信号スペクトルの振幅に対して或る一定の比率の誤差が残り、他のパルス位置（Iep[j][i]=0.0に対応する帯域）では、CELP残差信号スペクトルが誤差としてそのまま残るものとする。具体的には、推定パルス減衰部１０７は、次式（５）に従って、変換符号化推定残差スペクトルCraを算出する。
【数５】
2012032759.tif 000007
【００６３】
ここで、αは推定パルス位置においてCELP残差信号スペクトルの振幅をどの程度誤差として残すかを示す（つまり、減衰度合を示す）、０以上１未満の定数（以後、推定残差係数と呼ぶ）を表す。例えば、推定パルス位置における誤差を零と見なす場合には、α＝０．０に設定され、推定パルス位置において１０％の誤差を見込む場合には、α＝０．１に設定される。すなわち、推定パルス減衰部１０７は、CELP残差信号スペクトルの振幅に、推定残差係数（０以上１未満の値）を乗算することで、変換符号化推定残差スペクトル（つまり、復号信号スペクトルの推定値）を算出する。このように、０以上１未満の定数をCELP残差信号スペクトルに乗じて変換符号化による誤差を推定することは、変換符号化により所定のSNR(Signal Noise Ratio)が得られるように誤差を算出していることになる。このときのSNRは次式（６）で表される。
【数６】
2012032759.tif 000008
【００６４】
次いで、推定歪評価部１０８は、次式（７）に従って、入力信号スペクトル及び変換符号化推定残差スペクトルを用いて、変換符号化による符号化歪（歪エネルギ）の推定値である推定歪エネルギＥｅを算出する（以下、推定歪評価と呼ぶことがある）。
【数７】
2012032759.tif 000009
【００６５】
ここで、Ｓは入力信号スペクトルを表す。また、θはCELP抑圧係数毎に設定される一定値を表し、CELP抑圧係数間の推定歪エネルギの調整機能を有する。例えば、CELP抑圧係数（インデックスｊ）が零のときはθ[j]=1.0に設定され、CELP抑圧係数（インデックスｊ）が大きいほど、θ[j]=0.0に近づくように調整される。
【００６６】
このように、推定歪評価部１０８は、推定パルス位置におけるスペクトルの振幅を０以上１未満の比率に減衰させた変換符号化推定残差スペクトルに対する推定歪エネルギを算出する。これにより、推定歪評価部１０８では、変換符号化部１１０で符号化されると推定されたパルス位置での推定歪エネルギを、変換符号化部１１０での処理量よりも少ない処理量で推定することが可能となる。
【００６７】
なお、予備選択探索において、全てのCELP抑圧係数で推定歪評価を行う場合には、推定歪評価部１０８は、CELP抑圧係数インデックスを全て走査するように動作する。すなわち、推定歪評価部１０８は、CELP抑圧係数インデックスを全てCELP成分抑圧部１０４に出力する。一方、予備選択探索において、推定歪評価を行うCELP抑圧係数の候補を限定することも可能である。
【００６８】
例えば、CELP抑圧係数インデックスの総数がＭ＝４の場合に３候補のみを予備選択探索する場合を説明する。この時、最も強く抑圧する係数と最も弱く抑圧する係数とのうちいずれかを本選択探索から除外することで候補を絞る。まず、CELP抑圧係数インデックスｊ＝１及びｊ＝４に対する推定歪エネルギ（つまり、Ee[1]及びEe[4]）を算出する。次いで、推定歪評価部１０８は、Ee[1]がEe[4]よりも小さい場合には、CELP抑圧係数インデックスｊ＝２に対する推定歪エネルギ（つまり、Ee[2]）を算出し、Ee[4]がEe[1]よりも小さい場合には、CELP抑圧係数インデックスｊ＝３に対する推定歪エネルギ（つまり、Ee[3]）を算出する。すなわち、ｊ＝１、４及び（２又は３のいずれか一方）の３種類のCELP抑圧係数に限定して推定歪評価が行われ、予備選択探索が完了する。よって、推定歪評価部１０８は、３つのCELP抑圧係数に対してのみ推定歪評価を行えばよく、ｊ＝１〜４の４つのCELP抑圧係数を全て評価する場合と比べて、予備選択探索に要する処理量を約３／４に抑えることができる。
【００６９】
次いで、本選択候補限定部１０９は、推定歪エネルギの分布に基づいて、本選択探索の探索対象であるCELP抑圧係数（変換符号化に用いるCELP抑圧係数）の候補を限定する。つまり、本選択候補限定部１０９は、推定歪エネルギに基づいて、CELP抑圧係数コードブックに格納されている複数のCELP抑圧係数のうち、所定の数のCELP抑圧係数を予備選択する。以下、本選択候補限定部１０９での本選択探索の限定方法１及び２について説明する。なお、以下では、一例として、Ｍ＝４（ｊ＝１〜４）の場合について説明する。
【００７０】
＜方法１＞
方法１では、CELP抑圧係数の最も大きい係数と最も小さい係数とについて予備選択探索を行い、推定歪エネルギが大きい方は本選択探索で選択される可能性が小さいと判断し、そのCELP抑圧係数を本選択探索から除外することで、本選択探索の処理量を減らす。
【００７１】
上記を実現する方法を以下に説明する。まず、本選択候補限定部１０９には、CELP抑圧係数インデックスｊ＝１及びｊ＝４に対する推定歪エネルギ（つまり、Ee[1]及びEe[4]）が入力される。
【００７２】
（１）本選択候補限定部１０９は、Ee[1]とEe[4]とを比較する。
【００７３】
（２）Ee[1]がEe[4]よりも小さい場合、本選択候補限定部１０９は、本選択探索をｊ＝１，２，３の３種類のCELP抑圧係数に限定する。一方、Ee[4]がEe[1]よりも小さい場合、本選択候補限定部１０９は、本選択探索をｊ＝２，３，４の３種類のCELP抑圧係数に限定する。
【００７４】
本選択探索では、このようにして限定された３つのCELP抑圧係数（CELP抑圧係数インデックス）を用いる。
【００７５】
つまり、本選択候補限定部１０９は、CELP成分抑圧部１０４に格納されている複数のCELP抑圧係数のうち、最大値を用いた場合の推定歪エネルギと、最小値を用いた場合の推定歪エネルギとを比較（上記例では、最小のインデックスｊ＝１及び最大のインデックスｊ＝４を比較）して、推定歪エネルギが大きい方のCELP抑圧係数を本選択探索の対象（本選択探索のCELP抑圧係数群）から除外する。つまり、予備選択探索を行うことで、本選択探索における探索対象候補が１つ削減される。
【００７６】
このとき、符号化装置１００において、予備選択探索での演算回数（推定歪評価の回数）は２回（上記例ではｊ＝１，４の２回）となり、本選択探索での演算回数が３回（ｊ＝１，２，３又はｊ＝２，３，４）となる。このとき、本選択探索での変換符号化の１回の処理量（削減分）の方が、予備選択探索での２回の演算における処理量よりも大きい場合には、符号化装置１００全体での処理量は削減される。
【００７７】
このようにして、方法１では、必要最小限のCELP抑圧係数（ここでは、最大値と最小値との２個のCELP抑圧係数）についてのみで予備選択探索が行われる。また、方法１では、推定歪エネルギの大きいCELP抑圧係数が本選択探索の対象から除外される。これにより、本選択探索において全てのCELP抑圧係数を探索する場合と比較して、符号化の品質劣化を抑えつつ、符号化装置１００における処理量を削減することができる。
【００７８】
＜方法２＞
方法２では、全てのCELP抑圧係数で予備選択探索を行い、推定歪エネルギから本選択探索でも選択される可能性の高いCELP抑圧係数を限定することで、本選択探索の処理量を減らす。この時最も推定歪エネルギが小さい候補は必ず本選択探索の候補として残すようにする。そして、残された候補に付与されたCELP抑圧係数インデックスに隣接するインデックス（片方または両方）のCELP抑圧係数も本選択探索の候補として残すようにする。これは、CELP抑圧係数インデックスが抑圧の程度に関して昇順または降順に配置されている場合に、これらCELP抑圧係数候補が本選択探索時に歪エネルギが最も小さい候補として選択される可能性が、推定歪エネルギが最小の候補およびそれに隣接する候補以外のCELP抑圧係数候補よりも高いからである。
【００７９】
上記を実現する方法として、本選択探索で２種類のCELP抑圧係数を探索対象とする場合について説明する。
【００８０】
本選択候補限定部１０９には、全てのCELP抑圧係数（ｊ＝１〜４）に対する推定歪エネルギ（つまり、Ee[1]〜Ee[4]）が入力される。
【００８１】
（１）本選択候補限定部１０９は、推定歪エネルギEe[1]〜Ee[4]のうち、最小の推定歪エネルギを探索し、最小の推定歪エネルギに対応するCELP抑圧係数インデックスを保存する。
【００８２】
（２）本選択候補限定部１０９は、保存したCELP抑圧係数インデックス（つまり、最小の推定歪エネルギに対応するCELP抑圧係数インデックス）の前後（両端）のCELP抑圧係数インデックスに対応する推定歪エネルギを比較し、推定歪エネルギが小さい方のCELP抑圧係数インデックスを保存する。
【００８３】
（３）本選択候補限定部１０９は、（１）の処理で保存したCELP抑圧係数インデックス（つまり、最小の推定歪エネルギに対応するCELP抑圧係数インデックス）、及び、（２）の処理で保存したCELP抑圧係数インデックスの２種類のCELP抑圧係数を、本選択探索のCELP抑圧係数群として限定する。
【００８４】
本選択探索では、このようにして限定された２つのCELP抑圧係数（CELP抑圧係数インデックス）を用いる。
【００８５】
つまり、本選択候補限定部１０９は、CELP成分抑圧部１０４に格納されている複数のCELP抑圧係数のうち、推定歪エネルギが最小のCELP抑圧係数（第１のCELP抑圧係数）、及び、推定歪エネルギが最小のCELP抑圧係数の前後のCELP抑圧係数インデックスに対応するCELP抑圧係数のうち推定歪エネルギが小さいCELP抑圧係数（第２のCELP抑圧係数）を、本選択探索の対象として特定する。すなわち、本選択候補限定部１０９は、複数のCELP抑圧係数のうちの推定歪エネルギが最も小さいCELP抑圧係数（第１のCELP抑圧係数）と、推定歪エネルギが最も小さいCELP抑圧係数に付与されたCELP抑圧係数インデックスの前後のCELP抑圧係数インデックスに対応する２つのCELP抑圧係数のうち推定歪エネルギが小さい方のCELP抑圧係数（第２のCELP抑圧係数）と、を所定の数のCELP抑圧係数として予備選択する。
【００８６】
このとき、符号化装置１００において、予備選択探索での演算回数（推定歪評価の回数）は４回（ｊ＝１〜４）となり、本選択探索での演算回数が２回となる。このとき、本選択探索での変換符号化の２回の処理量（削減分）の方が、予備選択探索での４回の演算における処理量よりも大きい場合には、符号化装置１００全体での処理量は削減される。すなわち、方法１と同様、本選択探索での変換符号化の１回の処理量の方が、予備選択探索での２回の演算における処理量よりも大きい場合には、符号化装置１００全体での処理量は削減される。
【００８７】
このようにして、方法２では、全てのCELP抑圧係数を対象として予備選択探索が行われるものの、方法１と比較して、本選択探索対象であるCELP抑圧係数群をより狭く限定する。これにより、本選択探索における処理量を方法１よりも削減することができる。
【００８８】
また、方法２では、推定歪エネルギが最小のCELP抑圧係数、及び、当該CELP抑圧係数の両端のCELP抑圧係数インデックスに対応するCELP抑圧係数のうち推定歪エネルギがより小さいCELP抑圧係数が、本選択探索の対象となる。つまり、予備選択探索では、本選択探索において最適なCELP抑圧係数（歪エネルギが最小のCELP抑圧係数）として決定される可能性が高いCELP抑圧係数が探索される。よって、方法２では、本選択探索において全てのCELP抑圧係数を探索する場合と比較して、符号化の品質劣化を抑えつつ、符号化装置１００における処理量を削減することができる。
【００８９】
なお、方法２において、本選択候補限定部１０９は、CELP成分抑圧部１０４に格納されている複数のCELP抑圧係数のうち、推定歪エネルギが最小のCELP抑圧係数（例えば、CELP抑圧係数インデックスｊ）、及び、推定歪エネルギが最小のCELP抑圧係数の前後のCELP抑圧係数インデックスに対応するCELP抑圧係数群（例えば、CELP抑圧係数インデックス［ｊ−１］及び［ｊ＋１］）を、本選択探索の対象として特定してもよい。すなわち、本選択候補限定部１０９は、複数のCELP抑圧係数のうちの推定歪エネルギが最も小さいCELP抑圧係数と、推定歪エネルギが最も小さいCELP抑圧係数に付与されたインデックスの前後のインデックスに対応する２つのCELP抑圧係数と、を所定の数のCELP抑圧係数として予備選択してもよい。
【００９０】
以上、本選択候補限定部１０９での本選択探索の対象となるCELP抑圧係数群の限定方法１及び２について説明した。このように、方法１では、方法２と比較して、本選択探索の対象を広くすることで、本選択探索の対象を限定することによる本選択探索の性能劣化をより小さくすることができる。一方、方法２では、方法１と比較して本選択探索での処理量をより削減することができる。
【００９１】
このように、符号化装置１００では、予備選択探索において、推定歪評価部１０８が、予備選択探索で探索対象とするCELP抑圧係数インデックスをCELP成分抑圧部１０４に出力する。これにより、推定歪評価部１０８には、CELP抑圧係数インデックス毎に変換符号化推定残差スペクトルが入力され、推定歪評価部１０８は、CELP抑圧係数インデックスにそれぞれ対応する推定歪エネルギを算出する。そして、本選択候補限定部１０９は、推定歪エネルギに基づいて、実際に変換符号化を用いて歪評価を行う本選択探索で探索対象とするCELP抑圧係数インデックスを限定する。すなわち、符号化装置１００では、予備選択探索において、本選択探索での変換符号化の歪エネルギがより小さくなると見込まれる（推定される）CELP抑圧係数を特定する。
【００９２】
次いで、符号化装置１００では、本選択探索において、本選択候補限定部１０９から指示されるCELP抑圧係数インデックス群のみを用いて、変換符号化部１１０で変換符号化が行われ、歪評価部１１２で歪エネルギが最小となるCELP抑圧係数の探索が行われる。そして、歪エネルギが最小となるCELP抑圧係数に対応するCELP抑圧係数インデックスが多重化部１１３に出力され、当該CELP抑圧係数インデックスは、符号化装置１００の符号化データの一部として復号装置２００へ送信される。
【００９３】
つまり、本実施の形態では、符号化装置１００は、変換符号化で符号化されるパルス位置を統計的に推定し、推定したパルス位置で推定される推定歪エネルギを算出し、推定歪エネルギのより小さいCELP抑圧係数を、本選択探索の対象となるCELP抑圧係数群として限定する（予備選択探索）。そして、符号化装置１００は、予備選択探索にて候補が限定されたCELP抑圧係数毎に変換符号化を行い、残差信号のエネルギ（歪エネルギ）が最小となるCELP抑圧係数を決定する（本選択探索）。
【００９４】
こうすることで、符号化装置１００は、予備選択探索において、歪エネルギが小さいと見込まれるCELP抑圧係数のみを本選択探索の対象とすることで、変換符号化を行う回数を削減する。ここで、予備選択探索では、前述したように、パルス位置推定部１０６でのパルス位置の推定、推定パルス減衰部１０７での変換符号化推定残差スペクトルの算出、及び、推定歪評価部１０８での歪エネルギの算出を、それぞれ変換符号化部１１０での処理よりも少ない処理量で行うことが可能となる。よって、予備選択探索において本選択探索の対象となるCELP抑圧係数群を予め限定することにより、全てのCELP抑圧係数に対して変換符号化を逐次行う場合と比較して、符号化装置１００での処理量を削減することができる。
【００９５】
また、予備選択探索では、本選択探索の対象として、推定歪エネルギが小さいと見込まれるCELP抑圧係数、すなわち、本選択探索において歪エネルギが最小として評価される可能性が高いCELP抑圧係数のみに候補を限定する。これにより、本選択探索の対象となるCELP抑圧係数群を限定することによる、符号化の品質劣化を抑えることができる。
【００９６】
よって、本実施の形態によれば、音声信号に適した符号化と音楽信号に適した符号化とを階層構造にして組み合わせた符号化方式において、全てのCELP抑圧係数候補に対して変換符号化を逐次行う方法と比較して、符号化の品質劣化を抑えつつ、符号化装置における処理量を削減することができる。
【００９７】
なお、本実施の形態において、予備選択探索時に算出された値のうち、本選択探索時にも使用される値（例えば、CELP残差信号スペクトル等）については、本選択探索時に再度算出せずに、予備選択探索時に算出された値を利用してもよい。これにより、符号化装置では、本選択探索時の処理量を更に削減することができる。
【００９８】
（実施の形態２）
図３は、本発明の実施の形態２に係る符号化装置３００の主要な構成を示すブロック図である。なお、図３において、実施の形態１（図１）と同一の構成要素には同一の符号を付しその説明を省略する。図３に示す符号化装置３００では、図１に示す符号化装置１００に対してターゲット信号特徴抽出部３０１が追加される点が異なる。また、パルス位置推定部３０２及び推定パルス減衰部３０３には、ターゲット信号特徴抽出部３０１から出力される特徴情報が入力信号として追加される点が実施の形態１と異なる。
【００９９】
図３に示す符号化装置３００において、ターゲット信号特徴抽出部３０１は、CELP残差信号スペクトル算出部１０５から入力されるCELP残差信号スペクトル（ターゲット信号）を用いて、当該ターゲット信号の特徴を抽出する。
【０１００】
ここで、一例として、変換符号化としてFPC（Factorial Pulse Coding）を用いる場合について説明する。FPCでは、符号化対象（ここではCELP残差信号スペクトル）のスペクトルの振幅のばらつきが小さいときには符号化できるパルス本数がより多くなり、符号化対象のスペクトルの振幅のばらつきが大きいときには符号化できるパルス本数がより少なくなる、という特徴がある。例えば、或る帯域にエネルギが集中するターゲット信号では、FPCで符号化されるパルス本数は少なくなり、全帯域にエネルギが分散しているターゲット信号では、FPCで符号化されるパルス本数は多くなる。
【０１０１】
すなわち、符号化装置３００では、ターゲット信号（CELP残差信号スペクトル）の上記特徴を抽出して、抽出した特徴に基づいてFPCで符号化されるパルス本数を予測することができる。つまり、予備選択探索において、ターゲット信号のパルス位置を正確に推定することが可能となる。
【０１０２】
本実施の形態では、ターゲット信号特徴抽出部３０１は、ターゲット信号の振幅の平均値と、振幅の最大値との比を、ターゲット信号の特徴として抽出する。具体的には、ターゲット信号特徴抽出部３０１は、式（１）に従って、ターゲット信号の振幅の平均値Iavgを算出する。また、ターゲット信号特徴抽出部３０１は、ターゲット信号の絶対値振幅の最大値をtmaxとする。ここで、tmax/Iavgの値が大きいほど、或る特定の帯域にエネルギが集中している可能性が高い。つまり、tmax/Iavgの値が大きいほど、スペクトルのばらつきが大きい可能性が高い。
【０１０３】
よって、ターゲット信号特徴抽出部３０１は、tmax/Iavgの値が大きいほど、予備選択探索において推定するターゲット信号のパルス本数を少なくすべきであると判定する。一方、ターゲット信号特徴抽出部３０１は、tmax/Iavgの値が小さいほど、帯域全体にエネルギが分散している可能性が高いので、予備選択探索において推定するターゲット信号のパルス本数を多くすべきであると判定する。そこで、ターゲット信号特徴抽出部３０１は、tmax/Iavgの値に応じて、次式（８）に従って、ターゲット信号の特徴に基づいて予測されるターゲット信号のパルス本数に関する情報を特徴情報Ｋとして生成する。
【数８】
2012032759.tif 000010
【０１０４】
ここで、κhは予備選択探索（パルス位置推定部３０２）において推定されるパルスの本数を減少させるか否かを判定するために予め設定された閾値であり、κlは予備選択探索において推定されるパルスの本数を増加させるか否かを判定するために予め設定された閾値である。
【０１０５】
パルス位置推定部３０２は、CELP残差信号スペクトル算出部１０５から入力されるCELP残差信号スペクトル（ターゲット信号）、及び、ターゲット信号特徴抽出部３０１から入力される特徴情報Ｋを用いて、変換符号化部１１０で符号化されるパルス位置（推定パルス位置）を推定する。具体的には、パルス位置推定部３０２は、実施の形態１（パルス位置推定部１０６）で用いた式（３）の代わりに、次式（９）に示す閾値Ithr[j]を用いる。
【数９】
2012032759.tif 000011
【０１０６】
すなわち、式（９）では、特徴情報Ｋ（０．９，１．０，１．１）の値に応じてβの値がフレーム毎に適応的に補正され、パルス位置推定部３０２で選択されるパルス本数が適応的に制御される。換言すると、パルス位置推定部３０２は、式（９）に示すように、実施の形態１（式（３））を、ターゲット信号特徴抽出部３０１から入力される特徴情報Ｋを用いて補正する。
【０１０７】
これにより、パルス位置推定部３０２では、ターゲット信号において或る特定の帯域にエネルギが集中している可能性が高い場合（式（８）においてtmax/Iavg＞κhの場合）、特徴情報Ｋ＝１．１であるので、「β」が「β＊１．１」となり閾値Ithr[j]はより大きくなるように制御される。よって、パルス位置推定部３０２では、閾値Ithr[j]を超えるパルス本数がより少なくなる。
【０１０８】
一方、パルス位置推定部３０２では、ターゲット信号の全帯域にエネルギが分散している可能性が高い場合（式（８）においてtmax/Iavg＜κlの場合）、特徴情報Ｋ＝０．９であるので、「β」が「β＊０．９」となり閾値Ithr[j]はより小さくなるように制御される。よって、パルス位置推定部３０２では、閾値Ithr[j]を超えるパルス本数がより多くなる。
【０１０９】
すなわち、パルス位置推定部３０２は、式（８）においてtmax/Iavg＞κhの場合（スペクトルのばらつきが大きい場合）には、推定するパルス本数を少なく設定し、式（８）においてtmax/Iavg＜κlの場合（スペクトルのばらつきが小さい場合）には、推定するパルス本数を多く設定する。つまり、パルス位置推定部３０２は、CELP残差信号スペクトルの特徴に応じて、推定するパルスの本数を設定し、設定された本数のパルスの位置を推定する。例えば、パルス位置推定部３０２は、パルスの本数を、CELP残差信号スペクトルの各帯域における振幅のばらつきが大きくなるほど少なくするように設定する。
【０１１０】
推定パルス減衰部３０３は、ターゲット信号特徴抽出部３０１から入力される特徴情報を用いて、CELP残差信号スペクトル算出部１０５から入力されるCELP残差信号スペクトルのうち、パルス位置推定部３０２から入力される推定パルス位置のスペクトルを減衰させる。
【０１１１】
具体的には、推定パルス減衰部３０３は、実施の形態１（推定パルス減衰部１０７）で用いた式（５）の代わりに、次式（１０）に従って、変換符号化推定残差スペクトルＣｒａを算出する。
【数１０】
2012032759.tif 000012
【０１１２】
すなわち、式（１０）では、特徴情報Ｋ（０．９，１．０，１．１）の値に応じて推定残差計数αの値がフレーム毎に適応的に補正され、推定パルス減衰部３０３での減衰度合（推定誤差量）が適応的に制御される。換言すると、推定パルス減衰部３０３は、式（１０）に示すように、実施の形態１（式（５））を、ターゲット信号特徴抽出部３０１から入力される特徴情報Ｋを用いて補正する。
【０１１３】
これにより、推定パルス減衰部３０３では、ターゲット信号において或る特定の帯域にエネルギが集中している可能性が高い場合（式（８）においてtmax/Iavg＞κhの場合）、特徴情報Ｋ＝１．１であるので、「α」が「α／１．１」となり、推定パルス位置における誤差がより小さくなるように制御される。一方、推定パルス減衰部３０３では、ターゲット信号において全帯域にエネルギが分散している可能性が高い場合（式（８）においてtmax/Iavg＞κhの場合）、特徴情報Ｋ＝０．９であるので、「α」が「α／０．９」となり、推定パルス位置における誤差がより大きくなるように制御される。
【０１１４】
すなわち、推定パルス減衰部３０３は、式（８）においてtmax/Iavg＞κhの場合（スペクトルの振幅のばらつきが大きい場合）には、スペクトルの減衰度合を大きくし、式（８）においてtmax/Iavg＜κlの場合（スペクトルの振幅のばらつきが小さい場合）には、スペクトルの減衰度合を小さくする。すなわち、推定パルス減衰部３０３は、CELP残差信号スペクトルの減衰度合を、CELP残差信号スペクトルの各帯域における振幅のばらつきが大きくなるほど大きくするように設定する。
【０１１５】
換言すれば、スペクトルの振幅のばらつきに応じて、変換符号化の誤差の推定値により算出されるSNRが適応的に変化することになる。そのときのSNRは次式（１１）で表される。
【数１１】
2012032759.tif 000013
【０１１６】
このように、符号化装置３００は、ターゲット信号（CELP残差信号スペクトル）の特徴（ここでは、スペクトルの振幅のばらつき（tmax/Iavg））に応じて、変換符号化部１１０で符号化されるパルス本数及びパルスの誤差（推定パルス減衰部３０３での減衰度合）を適応的に制御する。これにより、符号化装置３００では、変換符号化部１１０で符号化されると推定されるパルス位置での歪エネルギを、実施の形態１よりも精度良く推定することができる。また、実施の形態１と同様、符号化装置３００では、推定パルス位置の推定、推定パルス減衰部１０７での変換符号化推定残差スペクトルの算出、及び、推定歪評価部１０８での歪エネルギの算出を、それぞれ変換符号化部１１０での処理よりも少ない処理量で行うことができる。
【０１１７】
よって、本実施の形態によれば、音声信号に適した符号化と音楽信号に適した符号化とを階層構造にして組み合わせた符号化方式において、実施の形態１と比較して、符号化の品質劣化を更に抑えつつ、全てのCELP抑圧係数候補に対して変換符号化を逐次行う方法と比較して、符号化装置における処理量を削減することができる。
【０１１８】
なお、本実施の形態では、ターゲット信号の特徴として、スペクトルの振幅のばらつきを用いる場合について説明したが、本発明は、ターゲット信号の特徴としてスペクトルの振幅のばらつきを用いる場合に限定されない。例えば、ターゲット信号の特徴として、ターゲット信号のトーン性を用いてもよい。ここでいうトーン性とは、スペクトルのピークの大きさ、若しくはダイナミックレンジの大きさを示す指標である。例えば、ターゲット信号又はその絶対値の算術平均に対する幾何平均の比を測定し、この比が０に近いときはトーン性が高いと判定することができる。具体的には、図３に示す符号化装置３００において、ターゲット信号特徴抽出部３０１は、ターゲット信号のトーン性を測定する。そして、パルス位置推定部３０２は、パルスの本数を、トーン性が高くなるほど少なくするように設定する。例えば、パルス位置推定部３０２は、ターゲット信号のトーン性が高い場合には閾値を大きく設定して、推定パルス本数が少なくなるように制御し、ターゲット信号のトーン性が低い場合には閾値を小さくして、推定パルス本数が多くなるように制御すればよい。また、推定パルス減衰部３０３は、CELP残差信号スペクトルの減衰度合を、トーン性が高くなるほど大きくするように設定する。つまり、推定パルス減衰部３０３は、ターゲット信号のトーン性が高い場合には推定残差係数を小さくして（減衰度合を大きくして）、残差信号（誤差）が小さくなるように制御し、ターゲット信号のトーン性が低い場合には推定残差係数を大きくして（減衰度合を小さくして）、残差信号（誤差）が大きくなるように制御すればよい。このように、ターゲット信号の特徴としてトーン性を用いる場合でも、本実施の形態と同様の効果を得ることができる。
【０１１９】
また、例えば、ターゲット信号の特徴として、ターゲット信号の雑音性を用いてもよい。ここでいう雑音性とはターゲット信号のエネルギの偏りの少なさを示す指標である。例えば、ターゲット信号をいくつかの帯域で区切って帯域毎のエネルギを測定し、帯域毎のエネルギの分散が小さいときは雑音性が高いと判定することができる。具体的には、図３に示す符号化装置３００において、ターゲット信号特徴抽出部３０１は、ターゲット信号の雑音性を測定する。そして、パルス位置推定部３０２は、パルスの本数を、雑音性が高くなるほど多くするように設定する。例えば、パルス位置推定部３０２は、ターゲット信号の雑音性が高い場合には閾値を小さく設定して、推定パルス本数が多くなるように制御し、ターゲット信号の雑音性が低い場合には閾値を大きくして、推定パルス本数が少なくなるように制御すればよい。また、推定パルス減衰部３０３は、CELP残差信号スペクトルの減衰度合を、雑音性が高くなるほど小さくするように設定する。つまり、推定パルス減衰部３０３は、ターゲット信号の雑音性が高い場合には推定残差係数を大きくして（減衰度合を小さくして）、残差信号（誤差）が大きくなるように制御し、ターゲット信号の雑音性が低い場合には推定残差係数を小さくして（減衰度合を大きくして）、残差信号（誤差）が小さくなるように制御すればよい。このように、ターゲット信号の特徴として雑音性を用いる場合でも、本実施の形態と同様の効果を得ることができる。
【０１２０】
以上、本発明の各実施の形態について説明した。
【０１２１】
なお、上記各実施の形態では、パルス位置推定部において、変換符号化部への入力信号（CELP残差信号スペクトル）が正規分布であると仮定し、振幅の大きい上位周波数を選択するための閾値（Ithr）を設定する場合について説明した。しかし、パルス位置推定部は、変換符号化部への入力信号（CELP残差信号スペクトル）が正規分布以外の他の分布を仮定できる場合には、当該分布モデルに応じて閾値（Ithr）を設定してもよい。
【０１２２】
また、上記各実施の形態では、パルス位置推定部において、変換符号化部で符号化されるパルス数の上限値を上回るパルス本数を推定する場合があり得る。これに対して、パルス位置推定部は、当該上限値を用いて、推定されるパルス数を制御してもよい。このとき、パルス位置推定部は、振幅がより小さいパルスを除外したり、より高域側のパルスを除外したりしてもよい。又は、パルス位置推定部は、上述した振幅及び周波数帯域の条件に加え、信号の特徴から算出できる他の条件を組み合わせて、除外するパルスを決定してもよい。
【０１２３】
また、上記各実施の形態では、CELP抑圧係数コードブックに格納されるCELP抑圧係数が、CELP抑圧の程度の昇順又は降順で格納されている場合について説明した。しかし、抑圧係数の候補を限定する方法として、格納される順序によらない方法を用いる場合には、必ずしも昇順又は降順としなくても良い。
【０１２４】
また、上記各実施の形態では、音声信号に適した符号化の一例としてCELP符号化を用いて説明したが、本発明はADPCM（Adaptive Differential Pulse Code Modulation)、APC（Adaptive Prediction Coding)、ATC（Adaptive Transform Coding）、TCX（Transform Coded Excitation）等を用いても実現可能であり、同様の効果が得られる。
【０１２５】
また、上記各実施の形態では、音楽信号に適した符号化の一例として変換符号化を用いて説明したが、音声信号に適した符号化方式の復号信号と入力信号との残差信号を周波数領域で効率良く符号化できる方式であれば良い。このような方式として、FPC（Factorial Pulse Coding）及びAVQ（Algebraic Vector Quantization）などがあり、同様の効果を得ることができる。
【０１２６】
また、以上の説明では、符号化装置１００、３００から出力された符号化データを復号装置２００で受信するとしたが、これに限るものではない。すなわち、復号装置２００は、符号化装置１００、３００の構成において生成された符号化データでなくても、復号に必要な符号化データを有する符号化データを生成可能な符号化装置により出力された符号化データであれば、復号可能である。
【０１２７】
また、上記各実施の形態では、本発明をハードウェアで構成する場合を例にとって説明したが、本発明はハードウェアとの連係においてソフトウェアでも実現することも可能である。
【０１２８】
また、上記各実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されてもよいし、一部または全てを含むように１チップ化されてもよい。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
【０１２９】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Field Programmable Gate Array）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル／プロセッサを利用してもよい。
【０１３０】
さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。バイオ技術の適用等が可能性としてありえる。
【０１３１】
２０１０年９月１０日出願の特願２０１０−２０３６５７の日本出願に含まれる明細書、図面および要約書の開示内容は、すべて本願に援用される。
【産業上の利用可能性】
【０１３２】
本発明は、符号化の品質劣化を抑えつつ装置全体としての演算量を削減することができ、例えば、パケット通信システム、移動通信システムなどに適用できる。
【符号の説明】
【０１３３】
１００，３００  符号化装置
２００  復号装置
１０１，１０３，２０４  MDCT部
１０２  CELP符号化部
１０４，２０５  CELP成分抑圧部
１０５  CELP残差信号スペクトル算出部
１０６，３０２  パルス位置推定部
１０７，３０３  推定パルス減衰部
１０８  推定歪評価部
１０９  本選択候補限定部
１１０  変換符号化部
１１１，２０６  加算部
１１２  歪評価部
１１３  多重化部
２０１  分離部
２０２  変換符号化復号部
２０３  CELP復号部
２０７  IMDCT部
３０１  ターゲット信号特徴抽出部
【図１】
2012032759.tif 000014
【図２】
2012032759.tif 000015
【図３】
2012032759.tif 000016
【国際調査報告】
2012032759.tif 000017
2012032759.tif 000018
2012032759.tif 000019
2012032759.tif 000020
2012032759.tif 000021
2012032759.tif 000022