(57)【要約】【課題】  予測符号化方式により符号化された音響信号を復号化して再生する際に、量子化雑音を有効に除去して再生音の主観的品質を向上させる。【解決手段】  量子化雑音除去のためのフィルタは、タップ長をパラメータとするフィルタ関数の族でその特性が規定される複数のディジタルフィルタであってフィルタ係数が正負で対称であるようなディジタルフィルタから選定され、ここで、フィルタ関数の族は、その伝達関数において、その利得が、零周波数における１からナイキスト周波数における０まで単調に減少する曲線となるようにフィルタ係数が決定されている。また、フィルタ関数の族からの特定のフィルタ関数の選定は、原信号と量子化雑音について、タップ長をパラメータとする関数の族で規定される複数のフィルタのそれぞれを適用した場合に除去されるエネルギーを算出し、各タップ長ごとに求めたそれらのエネルギー差に基づいて行われる。【選択図】  図１
2015087616.tif 000002
【特許請求の範囲】
【請求項１】
予測符号化に基づく符号化方式により符号化された音響信号を復号化して再生する音響再生装置における復号化後の信号に含まれる量子化雑音を除去するためのフィルタの設計方法であって、
前記フィルタは、タップ長をパラメータとするフィルタ関数の族でその特性が規定される複数のディジタルフィルタであってフィルタ係数が正負で対称であるようなディジタルフィルタから選定され、前記フィルタ関数の族は、その伝達関数において、その利得が、零周波数における１からナイキスト周波数における０まで単調に減少する曲線となるように前記フィルタ係数が決定されていることを特徴とするフィルタの設計方法。
【請求項２】
前記フィルタ関数の族ｈＬ（ｎ）の前記伝達関数ＨＬ（ｆ）を、下の第一式で規定し、
前記零周波数における１からナイキスト周波数における０まで単調に減少する曲線を規定するために、下の第二及び第三式を前記第一式に対する条件式として規定することを特徴とする請求項１に記載のフィルタの設計方法。
【数６】
2015087616.tif 000013
ここで、ｆは、Ｆ／Ｆｓ∈［０，１］、すなわち、周波数Ｆ［Hz］をサンプリング周波数Ｆｓ［Hz］で除した正規化周波数であり、Ｌは、タップ長であり、正の奇数である。
ＨＬ（１／２）＝０
【数７】
2015087616.tif 000014
【請求項３】
前記フィルタ関数の族からの特定のフィルタ関数の選定は、個々の音響について、その原信号と、その原信号から得られた再生信号との差である前記量子化雑音を求め、前記原信号と前記量子化雑音の両者について、前記タップ長をパラメータとする関数の族で規定される複数のディジタルフィルタのそれぞれを適用した場合に除去されるエネルギーを算出し、個々のタップ長ごとに、前記原信号に係る除去エネルギーと前記量子化雑音に係る除去エネルギーとの差を求め、その差が最も小さいタップ長に係る前記フィルタ関数をその特性として有するディジタルフィルタを最適なフィルタとして選定することを特徴とする請求項２に記載のフィルタの設計方法。
【請求項４】
前記除去されるエネルギーは、前記原信号又は前記量子化雑音に対して、注目するタップ長に係るフィルタの関数を周波数領域上で反転して得たハイパスフィルタを畳み込んだ二乗和として得ることを特徴とする請求項３に記載のフィルタの設計方法。
【請求項５】
前記タップ長Ｌ＝３，５，７，９について前記算出処理を行うことを特徴とする請求項３又は４に記載のフィルタの設計方法。
【請求項６】
前記ディジタルフィルタは、ＦＩＲディジタルフィルタであることを特徴とする請求項１乃至５のいずれかに記載のフィルタの設計方法。
【請求項７】
予測符号化に基づく符号化方式により符号化された音響信号を復号化して再生信号を生成する復号化部と、
前記復号化部から出力された前記再生信号に含まれる量子化雑音を除去するためのフィルタと、
を備える音響再生装置であって、
前記フィルタは、タップ長をパラメータとするフィルタ関数の族でその特性が規定される複数のディジタルフィルタであってフィルタ係数が正負で対称であるようなディジタルフィルタから選定され、前記フィルタ関数の族は、その伝達関数において、その利得が、零周波数における１からナイキスト周波数における０まで単調に減少する曲線となるように前記フィルタ係数が決定されていることを特徴とする音響再生装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、フィルタの設計方法及びそのフィルタを備えた音響再生装置に関し、特に、音響信号の符号化方式として適応的差分パルス符号変調方式を採用する場合のフィルタの設計方法及びそのフィルタを備えた音響再生装置に関する。
【背景技術】
【０００２】
時間軸方向に強い相関を有するディジタル信号を符号化する方式としては、予測値との差分を符号化することで符号量を削減する予測符号化方式がある。音響信号もかかる性質を有する信号であり、その符号化には、予測符号化方式の一つである適応的差分パルス符号変調方式（Adaptive Differential Pulse Code Modulation）（以下、「ＡＤＰＣＭ符号化方式」と称す）が広く利用されている。このＡＤＰＣＭ符号化方式は、予測符号化において、予測係数と量子化幅を適応的に変更する方式である。
【０００３】
図５（ａ）は、ＡＤＰＣＭ符号化方式に基づき、原信号ｓ（ｎ）からＡＤＰＣＭ出力信号ｘ（ｎ）を生成するＡＤＰＣＭ符号化装置の典型的な構成を示すブロック図である。同図に示すＡＤＰＣＭ符号化装置は、減算器１１、量子化器１２、逆量子化器１３、加算器１４、及び予測器１５を備えている。
【０００４】
図５（ａ）において、音響信号のサンプルｎについて、減算器１１により、原信号ｓ（ｎ）から予測値ｐ（ｎ）を差し引くことにより、差分信号ｄ（ｎ）が求められる。なお、予測値ｐ（ｎ）については後述する。量子化器１２は、入力する差分信号ｄ（ｎ）に対して適応量子化関数Ｑを施すことにより、符号化装置の出力としてのＡＤＰＣＭ出力信号ｘ（ｎ）を出力する。逆量子化器１３は、入力されるＡＤＰＣＭ出力信号ｘ（ｎ）に対して適応逆量子化関数Ｑ−１を施すことにより、量子化差分信号Ｑ−１（ｘ（ｎ））を出力する。逆量子化器１３から出力された量子化差分信号Ｑ−１（ｘ（ｎ））は、加算器１４及び予測器１５に入力される。加算器１４は、量子化差分信号Ｑ−１（ｘ（ｎ））と予測値ｐ（ｎ）とを加算して再生信号ｙ（ｎ）として出力する。予測器１５は、次の原信号ｓ（ｎ）を予測するもので、ＡＤＰＣＭ符号化方式においては、入力する再生信号ｙ（ｎ）について、少なくとも直近の過去２サンプルの再生信号ｙ（ｎ−１）及びｙ（ｎ−２）を記憶しており、それらと量子化差分信号Ｑ−１（ｘ（ｎ））に基づき、式（１）により得られた結果を前述の予測値ｐ（ｎ）として出力する。  
ｐ（ｎ）＝α１×ｙ（ｎ−１）＋α２×ｙ（ｎ−２）      ・・・（１）
ここで、α１及びα２は予測係数である。
【０００５】
図５（ａ）に示されたＡＤＰＣＭ符号化装置から出力されたＡＤＰＣＭ出力信号ｘ（ｎ）を復号化するためのＡＤＰＣＭ復号化装置については、同図（ｂ）に示している。但し、同図（ａ）を見れば分かる通り、同図（ｂ）に示したＡＤＰＣＭ復号化装置の構成は、同図（ａ）に示したＡＤＰＣＭ符号化装置の一部として含まれる構成と同一である。従って、その説明を省略する。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００８−２９８８８１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上述のＡＤＰＣＭ符号化方式を一例とする予測符号化方式は、高速な復号化処理が可能で、ハードウェア実装時の回路規模も小さくなるので、処理コスト等の制約から変換符号化が使用できない場面で特に有効な方式である。しかしながら、同程度の圧縮比で比較した場合、予測符号化方式は、変換符号化方式に比べて主観的品質が劣るという課題がある。その主要因は、量子化器の処理により本質的に生ずる量子化雑音による。
【０００８】
ここで、図５（ａ）における原信号ｓ（ｎ）と再生信号ｙ（ｎ）の差である量子化雑音ｎｑ（ｎ）は、式（２）に示すように、結果的に差分信号（予測誤差）ｄ（ｎ）に対する量子化誤差そのものに帰着する。  
ｎｑ（ｎ）＝ｓ（ｎ）−ｙ（ｎ）
＝ｄ（ｎ）＋ｐ（ｎ）−（Ｑ−１（ｘ（ｎ））＋ｐ（ｎ））
＝ｄ（ｎ）−Ｑ−１（Ｑ（ｄ（ｎ）））          ・・・（２）
【０００９】
そこで、量子化器１２における量子化係数をｑとした場合、量子化雑音ｎｑ（ｎ）は、剰余ｄ（ｎ）mod ｑで表現できるので、量子化係数ｑに対して差分信号（予測誤差）ｄ（ｎ）の振幅幅が大きく、ある程度の複雑性を有していれば、量子化雑音ｎｑ（ｎ）は、周波数領域全域に一様なエネルギーを有する信号とみなせる。一方、符号化の対象となる楽音の原信号ｓ（ｎ）は、一般にエネルギーが低周波域に集中し、高周波になるに従って減衰する傾向がある。
【００１０】
そこで、原信号ｓ（ｎ）のエネルギーと量子化雑音ｎｑ（ｎ）のエネルギーとを、想定する周波数全域に渡って比較したグラフを図６に示す。同図で明らかなように、楽音によっては、高周波域で量子化雑音ｎｑ（ｎ）のエネルギーが原信号ｓ（ｎ）のエネルギーを上回る場合があり、これが上述の主観的品質の劣化に直接影響していると考えられる。
【００１１】
なお、特許文献１には、ＡＤＰＣＭ符号化方式で符号化された音声信号を復号化して再生する音声再生装置において、再生時にノイズを抑制して再生する技術が開示されているが、その目的は、特定の処理、具体的には１６ｋｂｐｓＡＤＰＣＭ符号化におけるノイズフロアや無線エラー時のノイズを低減することであり（段落［００１５］）、量子化雑音に焦点を当てるものではない。また、その構成も、ノイズ判定の結果であるノイズの有無に応じて、フィルタの遮断周波数を変更（例えば、3000Hzから2000Hzに下げる）する、という単純なものである（段落［００２２］〜［００２３］、［００５２］〜［００５３］）。
【００１２】
本発明は上述のような事情から為されたものであり、本発明の目的は、予測符号化方式により符号化された音響信号を復号化して再生する際に、その再生信号に含まれる量子化雑音を有効に除去して再生音の主観的品質を向上させることができるフィルタの設計方法及びそのフィルタを備えた音響再生装置を提供することにある。
【課題を解決するための手段】
【００１３】
上記目的を達成するため、本発明のフィルタの設計方法は、予測符号化に基づく符号化方式により符号化された音響信号を復号化して再生する音響再生装置における復号化後の信号に含まれる量子化雑音を除去するためのフィルタの設計方法であって、前記フィルタは、タップ長をパラメータとするフィルタ関数の族でその特性が規定される複数のディジタルフィルタであってフィルタ係数が正負で対称であるようなディジタルフィルタから選定され、前記フィルタ関数の族は、その伝達関数において、その利得が、零周波数における１からナイキスト周波数における０まで単調に減少する曲線となるように前記フィルタ係数が決定されていることを要旨とする。なお、前記ディジタルフィルタとしては、ＦＩＲディジタルフィルタが好適である。
【００１４】
このとき、前記フィルタ関数の族ｈＬ（ｎ）の前記伝達関数ＨＬ（ｆ）を、下の第一式で規定し、前記零周波数における１からナイキスト周波数における０まで単調に減少する曲線を規定するために、下の第二及び第三式を前記第一式に対する条件式として規定することが好適である。
【数１】
2015087616.tif 000003
ここで、ｆは、Ｆ／Ｆｓ∈［０，１］、すなわち、周波数Ｆ［Hz］をサンプリング周波数Ｆｓ［Hz］で除した正規化周波数であり、Ｌは、タップ長であり、正の奇数である。
ＨＬ（１／２）＝０
【数２】
2015087616.tif 000004
【００１５】
更に、前記フィルタ関数の族からの特定のフィルタ関数の選定は、個々の音響について、その原信号と、その原信号から得られた再生信号との差である前記量子化雑音を求め、前記原信号と前記量子化雑音の両者について、前記タップ長をパラメータとする関数の族で規定される複数のディジタルフィルタのそれぞれを適用した場合に除去されるエネルギーを算出し、個々のタップ長ごとに、前記原信号に係る除去エネルギーと前記量子化雑音に係る除去エネルギーとの差を求め、その差が最も小さいタップ長に係る前記フィルタ関数をその特性として有するディジタルフィルタを最適なフィルタとして選定することを要旨とする。
【００１６】
このとき、前記除去されるエネルギーは、前記原信号又は前記量子化雑音に対して、注目するタップ長に係るフィルタの関数を周波数領域上で反転して得たハイパスフィルタを畳み込んだ二乗和として得られる。
【００１７】
また、前記タップ長Ｌ＝３，５，７，９について前記算出処理を行うことが好適である。
【００１８】
また、上記目的を達成するため、本発明の音響再生装置は、予測符号化に基づく符号化方式により符号化された音響信号を復号化して再生信号を生成する復号化部と、前記復号化部から出力された前記再生信号に含まれる量子化雑音を除去するためのフィルタと、を備える音響再生装置であって、前記フィルタは、タップ長をパラメータとするフィルタ関数の族でその特性が規定される複数のディジタルフィルタであってフィルタ係数が正負で対称であるようなディジタルフィルタから選定され、前記フィルタ関数の族は、その伝達関数において、その利得が、零周波数における１からナイキスト周波数における０まで単調に減少する曲線となるように前記フィルタ係数が決定されていることを要旨とする。
【発明の効果】
【００１９】
本発明によれば、予測符号化に基づく符号化方式により符号化された音響信号を復号化して再生する際に、その再生信号に含まれる量子化雑音を有効に除去して再生音の主観的品質を向上させることができる。
【００２０】
特に、上述の具体的なフィルタ関数の群を規定する式は、原信号と量子化雑音の周波数領域におけるエネルギー分布の特徴を鑑みて定めているので、確実に有効である。
【００２１】
また、上述のフィルタ関数の族からの特定のフィルタ関数の選定においては、個々の音響について実際に量子化雑音を求めて、その量子化雑音に基づいてフィルタの選定を行っているので、最適なフィルタが期待できる。
【図面の簡単な説明】
【００２２】
【図１】本発明のフィルタの設計方法の一実施形態の手順を示すフローチャートである。
【図２】最適なフィルタ関数を選定する処理の機能構成図である。
【図３】表１に示したフィルタ係数ｈＬ（ｎ）により定まる関数に基づくフィルタの周波数特性を、各タップ長Ｌについて示す図である。
【図４】本発明のフィルタの設計方法により選定されたフィルタを後段に含む音響再生装置の構成ブロック図である。
【図５】ＡＤＰＣＭ符号化装置及びＡＤＰＣＭ復号化装置の構成を示すブロック図である。
【図６】原信号ｓ（ｎ）のエネルギーと量子化雑音ｎｑ（ｎ）のエネルギーとを、想定する周波数全域に渡って比較したグラフを示す図である。
【発明を実施するための形態】
【００２３】
以下、図面を参照して、本発明の実施の形態について詳細に説明する。  
想定する周波数領域における原信号ｓ（ｎ）のエネルギーと量子化雑音ｎｑ（ｎ）のエネルギーとの関係は、図６に示したようになっているので、予測符号化の一例であるＡＤＰＣＭ符号化方式により符号化された音響信号を復号して再生する音響再生装置における復号化の後処理用に設けられるフィルタを適切に設計することにより、量子化雑音ｎｑ（ｎ）を有効に除去し得ることに着目した。
【００２４】
まず、前提として、通過域の品質を維持するために、かかるフィルタは、通過域の利得と位相を変えない特性を有していることが望ましい。また、図６に示すように、音響信号の周波数スペクトルは緩やかに変化するので、どの周波数以上を除去すべきか、その閾値が曖昧であることから、位相が保たれる特性（線形位相特性）を考慮する必要がある。従って、フィルタとしては、柔軟な特性が設計できるよう、線形位相特性のディジタルフィルタ（特に、有限インパルス応答（Finite Impulse Response：ＦＩＲ）フィルタ）により実現することが好適である。
【００２５】
図１は、本発明のフィルタの設計方法の一実施形態の手順を示すフローチャートである。  
そこで、かかるフィルタの設計方法としては、同図（ａ）を参照して、まず、フィルタ関数族を規定する（ステップＳ１）。具体的には、タップ長ＬをパラメータとするディジタルフィルタｈＬ（ｎ）を考える。ここで、制約として、ｈＬ（ｎ）が零位相特性を有するよう、フィルタ係数を正負で対称となるようにする。すなわち、ｈＬ（ｎ）＝ｈＬ（−ｎ）とする。すると、当該フィルタの伝達関数ＨＬ（ｆ）は、式（３）のように表現される。なお、ここでＬは正の奇数であり、ｆは、Ｆ／Ｆｓ∈［０，１］、すなわち、周波数Ｆ［Hz］をサンプリング周波数Ｆｓ［Hz］で除した正規化周波数である。  
【数３】
2015087616.tif 000005
【００２６】
式（３）に対して更に以下の２つの制約条件を導入する。  
条件［１］：式（４）で示すように、ナイキスト周波数ｆ＝１／２における利得を零にする。  
ＨＬ（１／２）＝０          ・・・（４）
条件［２］：低周波域の特性を平坦にするため、式（５）に示すように、零周波数ｆ＝０に関するｍ階微分までの条件を導入する。  
【数４】
2015087616.tif 000006
【００２７】
条件［１］における式（４）と条件［２］における式（５）とを式（３）に代入したものを連立させて解くことで、一意にフィルタ係数ｈＬ（ｎ）を求めることができる。このようにステップＳ１においては、各タップ長Ｌについて求められたフィルタ係数ｈＬ（ｎ）により定まるフィルタ関数の各タップ長についての集合体をフィルタ関数族として規定する。例として、タップ長Ｌ＝｛３，５，７，９｝について求めたフィルタ係数ｈＬ（ｎ）の一例を表１に示す。  
【表１】
2015087616.tif 000007
また、表１に示したフィルタ係数ｈＬ（ｎ）により定まる関数に基づくフィルタの周波数特性を、各タップ長Ｌについて、図３に示す。
【００２８】
図１（ａ）に戻り、次に、ステップＳ１で規定されたフィルタ関数から、個々の音響信号について、量子化雑音ｎｑ（ｎ）を除去するのに最適なフィルタ関数（特定のタップ数）を選定する（ステップＳ２）。このステップＳ２の処理の詳細を示すフローチャートを図１（ｂ）に示し、その処理の機能構成図を図２に示す。概要としては、量子化雑音ｎｑ（ｎ）のエネルギーが原信号ｓ（ｎ）のエネルギーを上回る帯域に対して、より除去機能が働くフィルタが品質向上に有効であるので、原信号ｓ（ｎ）と量子化雑音ｎｑ（ｎ）の周波数分布を比較することにより、タップ長Ｌを決定して特定のフィルタ関数を選定する。
【００２９】
詳細には、図１（ｂ）及び図２を参照して、まず、フィルタ関数を選定しようとする音響を選択してくる（ステップＳ２１）。次に、選択した音響の原信号ｓ（ｎ）から、図５に示すような使用する符号化装置により再生信号ｙ（ｎ）を求め、更に、図２に示すように、減算器２により原信号ｓ（ｎ）から再生信号ｙ（ｎ）を減算して量子化雑音ｎｑ（ｎ）を求める（ステップＳ２２）。
【００３０】
次に、想定するタップ長Ｌのうちから１つのタップ長Ｌを選択することにより、ステップＳ１において規定したフィルタ関数族の内から１つの関数を選択して１つのフィルタｈＬ（ｎ）を選択する（ステップＳ２３）。そして、選択したフィルタｈＬ（ｎ）により除去される成分を得るために、その通過・遮断特性を逆にしたハイパスフィルタｈＬｒ（ｎ）を導入する（ステップＳ２４）。なお、表記の便宜上、フィルタｈに対応するハイパスフィルタをｈｒと以後表記する。そこで、フィルタｈＬ（ｎ）に対応するハイパスフィルタｈＬｒ（ｎ）は、それぞれの伝達関数の関係が式（６）で表されることから、求められる。  
ＨＬｒ（ｆ）＝１−ＨＬ（ｆ）          ・・・（６）
図３には、ハイパスフィルタの例として、タップ長Ｌ＝９のフィルタに対応するハイパスフィルタが示されている。
【００３１】
選択したタップ長ＬにおけるハイパスフィルタｈＬｒ（ｎ）が得られると、次に、原信号ｓ（ｎ）及び量子化雑音ｎｑ（ｎ）について、選択したフィルタｈＬ（ｎ）により除去される除去されるエネルギーＥｓ及びエネルギーＥｎを求める（ステップＳ２５）。各エネルギーＥｓ及びＥｎは、ステップＳ２４において求めたハイパスフィルタｈＬｒ（ｎ）を利用して、それぞれ式（７）、式（８）により算出できる。すなわち、ハイパスフィルタｈＬｒ（ｎ）を畳み込んだ［ｓ＊ｈＬｒ］（ｎ）及び［ｎｑ＊ｈＬｒ］（ｎ）の二乗和である  。
【数５】
2015087616.tif 000008
【００３２】
次に、ステップＳ２５で求めたそれぞれのエネルギーから、それらの差の絶対値｜Ｅｎ−Ｅｓ｜を求める（ステップＳ２６）。
【００３３】
そして、ステップＳ２３〜ステップＳ２６までの処理を、想定するタップ長Ｌ、言い換えれば想定するフィルタ関数すべてについて行う（ステップＳ２７）。つまり、ステップＳ２７において、まだ検討すべきフィルタ関数が残っていればステップＳ２３に戻って次のタップ長Ｌを選択し、一方、すべて終了したならば、ステップＳ２８に移行する。なお、想定するタップ長Ｌ、すなわちいくつまでのタップ長について検討するかの範囲については、最適なタップ長Ｌが得られる期待と、得られた最適なタップ長Ｌによるフィルタのための計算量又は回路規模とのトレードオフの観点から決めればよい。つまり、大きなタップ長Ｌに処理を延ばせば延ばすほど最適なタップ長Ｌが求まる可能性は高いが、逆にそのタップ長Ｌを求める計算量又は回路規模は大きくなり、またそれを実現するフィルタリング処理の計算量又は回路規模は大きくなってしまう。
【００３４】
ところで、エネルギーＥｓよりエネルギーＥｎが小さい場合、選択したフィルタｈＬ（ｎ）により除去される帯域では原信号ｓ（ｎ）よりも量子化雑音ｎｑ（ｎ）のエネルギーが大きいということであり、より低域までフィルタリングを行ってもよいということが分かる。逆に、エネルギーＥｓよりエネルギーＥｎが大きい場合、原信号ｓ（ｎ）よりも量子化雑音ｎｑ（ｎ）のエネルギーが小さいということであり、過剰にフィルタリングをしているといえる。かかる観点から、ステップＳ２６で求めた、各タップ長Ｌ、言い換えれば各フィルタ、についての各エネルギー差絶対値｜Ｅｎ−Ｅｓ｜が、小さいほど最適にフィルタリングを行っていることとなる。従って、ステップＳ２８では、各エネルギー差絶対値｜Ｅｎ−Ｅｓ｜を比較する。そして、それが最も小さいタップ長で定まるフィルタ関数を選定する（ステップＳ２９）。
【００３５】
次に、上述したフィルタの設計方法により選定されたフィルタを後処理フィルタとして含む音響再生装置の一実施形態について説明する。図４は、かかる音響再生装置の構成ブロック図である。同図に示した音響再生装置は、逆量子化器１３、加算器１４、予測器１５、フィルタ１、及びスピーカ５０を少なくとも備えている。なお、増幅器等は省略している。そこで、逆量子化器１３は、入力されるＡＤＰＣＭ出力信号ｘ（ｎ）に対して逆量子化を施して量子化差分信号Ｑ−１（ｘ（ｎ））を出力する。予測器１５は、入力される量子化差分信号Ｑ−１（ｘ（ｎ））及び再生信号ｙ（ｎ）に基づいて、予測値ｐ（ｎ）を生成する。加算器１４は、逆量子化器１３から出力された量子化差分信号Ｑ−１（ｘ（ｎ））と、予測器１５から出力された予測値ｐ（ｎ）を加算して再生信号ｙ（ｎ）として出力する。フィルタ１は、前述のように、特定の音源に対して最適に設計されたフィルタであり、入力される再生信号ｙ（ｎ）に対してフィルタリング処理を行う。そして、最終的に音響としてスピーカ５０から出力される。
【実施例１】
【００３６】
具体的ないくつかの音響について評価を実施した。なお、この実施例では、音質の評価としてＰＥＡＱ客観評価法（ＩＴＵ−Ｒ勧告BS.1387-1）を採用した。ＰＥＡＱは、基準信号と評価対象信号を入力とし、基準信号からの音質劣化度をＯＤＧ（Objective Difference Grade）で表すというものである。このＯＤＧと評価語の対応を表２に示す。  
【表２】
2015087616.tif 000009
【００３７】
今回、具体的な音響信号として、表３に示すような性質の異なる６つの信号（音響番号Ｔ１〜Ｔ６）を例として採用した。ここで、これらの音響信号はすべて、サンプリング周波数44.1[kHz]の16[bit]で直線量子化された２[channel]のディジタル信号である。また、図１の手順における想定するタップ長Ｌとしては、｛３，５，７，９｝の範囲とした。  
【表３】
2015087616.tif 000010
【００３８】
各音響信号Ｔ１〜Ｔ６について、フィルタなしの場合の再生信号ｙ（ｎ）についてのＯＤＧと、Ｌ＝３，５，７，９の４つのフィルタを採用した場合の再生信号ｙ（ｎ）についてのＯＤＧを表４に示す。なお、音質比較に使用した基準信号はすべて符号化前の原信号である。表４から分かることは、フィルタを採用した場合には、フィルタなしの場合と比較してＯＤＧが概して向上している。故に、まず、本発明のフィルタ設計方法におけるフィルタ関数族の規定のしかたは有意に効果があると評価できる。詳細には、表４にあるように、音響番号Ｔ１についてＬ＝９が最適であり、また、Ｔ２についてはＬ＝３、Ｔ３についてはＬ＝７、Ｔ４についてはＬ＝５、Ｔ５についてはＬ＝７、Ｔ６についてはＬ＝９が最適と評価された。特に、音響番号Ｔ５については、Ｌ＝７のとき、ＯＤＧが2.13上昇している。  
【表４】
2015087616.tif 000011
【００３９】
一方、本発明のフィルタ設計方法による、タップ長Ｌ、すなわち特定のフィルタ関数、の選定の結果を表５に示す。これによれば、音響番号Ｔ３及びＴ６については、本発明により選定されたタップ長Ｌと最適なＯＤＧを示すタップ長Ｌ（表４）が一致した。また、他の音響番号についても近接したタップ長Ｌが選定されており、表４の最適なＯＤＧと、本発明により選定されたＬについて表４が示すＯＤＧとの差は、いずれも僅かである。また、本発明により選定されたタップ長Ｌについて表４が示すＯＤＧと、フィルタなしの場合のＯＤＧとから、ＯＤＧは平均0.76程度向上していることが分かる。  
【表５】
2015087616.tif 000012
【００４０】
以上のように、上述の本発明の実施形態によれば、音響再生装置の復号化処理の後処理として入れるフィルタの特性を表す関数として、式（３）〜（５）により規定される関数族を採用しているので、例えば単純な遮断周波数を有するフィルタと比較して、個々の音響に応じて量子化雑音をきめ細やかに排除して主観的品質を向上させることができる。また、短いタップ長のフィルタで実装できるので、ＡＤＰＣＭ符号化方式を代表とする符号化方式の特徴である処理速度や回路規模の利点を生かしつつ、再生音の主観的品質を向上させることができる。
【００４１】
特に、式（３）〜（５）は、原信号と量子化雑音の周波数領域におけるエネルギー分布の特徴を鑑みて定めているので、極めて有効である。
【００４２】
また、フィルタ関数の族からの特定のフィルタ関数の選定においては、個々の音響について実際に量子化雑音を求めて、その量子化雑音に基づいてフィルタの選定を行っているので、最適なフィルタが期待できる。
【産業上の利用可能性】
【００４３】
本発明のフィルタの設計方法及びそのフィルタを備えた音響再生装置は、音響信号を、例えばＡＤＰＣＭ符号化方式をはじめとする予測符号化方式で符号化する場合に採用される。
【符号の説明】
【００４４】
１  フィルタ
２  加算器
３  ハイパスフィルタ
１１  減算器
１２  量子化器
１３  逆量子化器
１４  加算器
１５  予測器
５０  スピーカ
【図１】
2015087616.tif 000015
【図２】
2015087616.tif 000016
【図３】
2015087616.tif 000017
【図４】
2015087616.tif 000018
【図５】
2015087616.tif 000019
【図６】
2015087616.tif 000020