WO2011055830A1 - 信号処理方法、情報処理装置、及び信号処理プログラム - Google Patents

Abstract

　劣化信号中の衝撃音を十分に抑圧すること。　劣化信号中の衝撃音を抑圧する。そのために、劣化信号中において衝撃音を検出する。そして、検出された衝撃音の位相情報を、位相情報の変化量が小さくなるように劣化信号中の衝撃音以外の信号の位相情報で処理することを特徴とする。

Description

信号処理方法、情報処理装置、及び信号処理プログラム

　本発明は、劣化信号中の雑音を抑圧して所望の信号を強調するための信号処理技術に関する。

　劣化信号（所望の信号に雑音が重畳された信号）から、雑音を抑圧し、強調信号（所望の信号を強調した信号）を出力する雑音抑圧技術（ｎｏｉｓｅ　ｓｕｐｐｒｅｓｓｉｎｇ　ｔｅｃｈｎｏｌｏｇｙ）が知られている。例えば、ノイズサプレッサは、所望の音声信号に重畳されている雑音（ノイズ）を抑圧するシステムであり、携帯電話など様々な音声端末において利用されている。
　この種の技術に関し、特許文献１には、入力信号に１より小さな抑圧係数を乗算することによって、ノイズを抑圧する方法が開示されている。また、特許文献２には、推定された雑音を劣化信号から直接減算することによって、雑音を抑圧する方法が開示されている。しかし、特許文献１及び２に記載の技術は、雑音の推定に平均操作が含まれており、衝撃音のような突発的な雑音を十分に抑圧することができない。
　これに対し、非特許文献１は、劣化信号パワースペクトルの平坦度と過去からの増分に基づいて衝撃音を検出する雑音抑圧システムを開示している。非特許文献１に記載のシステムは、非音声区間で衝撃音が検出されなかったときに、背景雑音を推定する。非特許文献１に記載のシステムは、非音声区間で衝撃音が検出されたときには、劣化信号を背景雑音の推定値で置換することで衝撃音を抑圧し、劣化信号と背景雑音との差を用いて衝撃音推定値を更新する。非特許文献１に記載のシステムは、音声区間で衝撃音が検出されたときには、衝撃音推定値を劣化信号から減算することで衝撃音を抑圧する。

特許第４２８２２２７号 特開平８−２２１０９２号

Ａ．Ｓｕｇｉｙａｍａ，Ｓｉｎｇｌｅ−ｃｈａｎｎｅｌ　ｉｍｐａｃｔ−ｎｏｉｓｅ　ｓｕｐｐｒｅｓｓｉｏｎ　ｗｉｔｈ　ｎｏ　ａｕｘｉｌｉａｒｙ　ｉｎｆｏｒｍａｔｉｏｎ　ｆｏｒ　ｉｔｓ　ｄｅｔｅｃｔｉｏｎ，″Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ＷＡＳＰＡＡ２００７，ｐｐ．１２７−−１３０，Ｏｃｔ．２００７．

　しかしながら、上述の非特許文献１に開示された構成では、位相に関して衝撃音抑圧処理を適用しておらず、位相の不連続性がそのまま存在した。その結果、ユーザにとって、衝撃音の抑圧が十分に感じられない場合があった。
　以上を踏まえ、本発明は、上述の課題を解決する信号処理技術を提供することを目的とする。

　上記目的を達成するため、本発明に係る信号処理方法は、劣化信号中の衝撃音を抑圧するために、前記劣化信号中において衝撃音を検出し、検出された衝撃音の位相情報を、位相情報の変化量が小さくなるように前記劣化信号中の前記衝撃音以外の信号の位相情報で処理する。
　上記目的を達成するため、本発明に係る情報処理装置は、劣化信号中の衝撃音を抑圧する情報処理装置であって、前記劣化信号中において衝撃音を検出する検出手段と、検出された衝撃音の位相情報を、前記劣化信号中の前記衝撃音以外の信号の位相情報で処理する位相処理手段と、を備える。
　上記目的を達成するため、本発明に係るプログラム記録媒体に格納される信号処理プログラムは、劣化信号中の衝撃音を抑圧する信号処理プログラムであって、コンピュータに、前記劣化信号中において衝撃音を検出する工程と、検出された衝撃音の位相情報を、前記劣化信号中の前記衝撃音以外の信号の位相情報を用いて処理する工程と、を実行させる。

　本発明によれば、劣化信号中の位相情報に関して衝撃音抑圧処理を適用することにより、位相に起因する信号の不連続性を低減し、衝撃音を十分に抑圧することができる。

本発明の第１実施形態としての雑音抑圧装置１００の概略構成を示すブロック図。 本発明の第１実施形態としての雑音抑圧装置１００に含まれる変換部２の構成を示すブロック図。 本発明の第１実施形態としての雑音抑圧装置１００に含まれる逆変換部４の構成を示すブロック図。 本発明の第１実施形態としての雑音抑圧装置１００に含まれる衝撃音抑圧部１１の構成を示すブロック図。 本発明の第２実施形態としての雑音抑圧装置１００に含まれる衝撃音検出部１０の構成を示すブロック図。 本発明の第３実施形態としての雑音抑圧装置３００の概略構成を示すブロック図。 本発明の第３実施形態としての雑音抑圧装置３００に含まれる衝撃音抑圧部１２の構成を示すブロック図。 本発明の第４実施形態としての雑音抑圧装置４００の概略構成を示すブロック図。 本発明の第５実施形態としての雑音抑圧装置に含まれる衝撃音抑圧部５２の構成を示すブロック図。 本発明の第６実施形態としての雑音抑圧装置の概略構成を示すブロック図。 本発明の第６実施形態としての雑音抑圧装置に含まれる衝撃音抑圧部６２の構成を示すブロック図。 本発明の第７実施形態としての雑音抑圧装置に含まれる衝撃音抑圧部７２の構成を示すブロック図。 本発明の第８実施形態としての雑音抑圧装置８００の概略構成を示すブロック図。 本発明の第９実施形態としての雑音抑圧装置９００の概略構成を示すブロック図。 本発明の第１０実施形態としての雑音抑圧装置１０００の概略構成を示すブロック図。 本発明の他の実施形態としての信号処理プログラムを実行するコンピュータ１１００の概略構成図。 雑音抑圧装置１００の別の構成例を示すブロック図。

　以下に、図面を参照して、本発明の実施の形態について例示的に詳しく説明する。ただし、以下の実施の形態に記載されている構成要素はあくまで例示であり、本発明の技術範囲をそれらのみに限定する趣旨のものではない。
　なお、本明細書中の「雑音」は、処理対象となる情報以外の不要な情報一般を示し、音に限定されるものではない。また、本明細書中の「衝撃音」は、雑音の一種であり、短時間に急激な変化を示す情報一般を示し、音に限定されるものではない。
　（第１実施形態）
　＜全体構成＞
　本発明に係る信号処理方法を実現する第１実施形態として、雑音抑圧装置について説明する。図１は、雑音抑圧装置１００の全体構成を示すブロック図である。図１６は、雑音抑圧装置１００の別の構成例を示すブロック図である。雑音抑圧装置１００は、例えばデジタルカメラ、ノートパソコン、携帯電話などといった装置の一部としても機能するが、本発明はこれに限定されるものではなく、入力信号からのノイズ除去を要求されるあらゆる情報処理装置に適用可能である。
　入力端子１には、劣化信号（所望信号と雑音の混在する信号）が、サンプル値系列として供給される。入力端子１に供給された劣化信号は、変換部２においてフーリエ変換などの変換を施されて複数の周波数成分に分割される。複数の周波数成分のうち、振幅は振幅スペクトルとして多重化され、衝撃音検出部１０と逆変換部４に伝達される。一方、位相は、位相スペクトルとして衝撃音抑圧部１１へ供給される。
　衝撃音検出部１０は、劣化信号スペクトルの周波数特性と時間特性に基づいて、衝撃音の存在を検出する。衝撃音検出部１０は、検出に当って、周波数特性と時間特性のいずれかを用いてもよいし、双方を用いてもよい。また、双方を用いる際に、衝撃音検出部１０は、それぞれの特性評価結果の重み付き和またはより複雑な関数で表現される統合結果を用いることもできる。衝撃音抑圧部１１は、変換部２から供給された劣化信号に対し、衝撃音検出部１０から供給された衝撃音検出情報に基づいて、各周波数で衝撃音を抑圧し、衝撃音抑圧結果を強調信号位相スペクトルとして逆変換部４に伝達する。
　逆変換部４は、衝撃音抑圧部１１から供給された強調信号位相スペクトルと変換部２から供給された劣化信号振幅スペクトルとを合わせて逆変換を行い、強調信号サンプルとして、出力端子５に供給する。
　＜変換部２の構成＞
　図２は、変換部２の構成を示すブロック図である。図２に示すように、変換部２はフレーム分割部２１、窓がけ処理部（ｗｉｎｄｏｗｉｎｇ　ｕｎｉｔ）２２、及びフーリエ変換部２３を含む。劣化信号サンプルは、フレーム分割部２１に供給され、Ｋ／２サンプル毎のフレームに分割される。ここで、Ｋは偶数とする。フレームに分割された劣化信号サンプルは、窓がけ処理部２２に供給され、窓関数（ｗｉｎｄｏｗ　ｆｕｎｃｔｉｏｎ）であるｗ（ｔ）との乗算が行なわれる。第ｎフレームの入力信号ｙｎ（ｔ）（ｔ＝０，１，．．．，Ｋ／２−１）に対するｗ（ｔ）で窓がけ（ｗｉｎｄｏｗｉｎｇ）された信号は、次式（１）で与えられる。

　また、窓がけ処理部２２は、連続する２フレームの一部を重ね合わせ（オーバラップ）して窓がけしてもよい。オーバラップ長としてフレーム長の５０％を仮定すれば、ｔ＝０，１，．．．，Ｋ／２−１に対して、以下の式（２）で得られる左辺が、窓がけ処理部２２の出力となる。

　窓がけ処理部２２は、実数信号に対しては、左右対称窓関数を用いてもよい。また、窓関数は、ＭＭＳＥ　ＳＴＳＡ法における抑圧係数を１に設定したとき、又はＳＳ法においてゼロを減算したときの入力信号と出力信号が計算誤差を除いて一致するように設計される。これは、ｗ（ｔ）＋ｗ（ｔ＋Ｋ／２）＝１となることを意味する。
　以後、連続する２フレームの５０％をオーバラップして窓がけする場合を例として説明を続ける。窓かけ処理部２２は、ｗ（ｔ）として、例えば、次式（３）に示すハニング窓を用いてもよい。

　このほかにも、ハミング窓、ケイザー窓、ブラックマン窓など、様々な窓関数が知られている。窓がけされた出力はフーリエ変換部２３に供給され、劣化信号スペクトルＹｎ（ｋ）に変換される。劣化信号スペクトルＹｎ（ｋ）は位相と振幅に分離され、劣化信号位相スペクトルａｒｇ　Ｙｎ（ｋ）は衝撃音抑圧部１１に、劣化信号振幅スペクトル｜Ｙｎ（ｋ）｜は衝撃音検出部１０と逆変換部４に供給される。既に説明したように、変換部２は、振幅スペクトルの代わりにパワースペクトルを利用することもできる。
　＜逆変換部４の構成＞
　図３は、逆変換部４の構成を示すブロック図である。図３に示すように、逆変換部４は逆フーリエ変換部４３、窓がけ処理部４２、及び、フレーム合成部４１を含む。逆フーリエ変換部４３は、衝撃音抑圧部１１から供給された強調信号位相スペクトルと変換部２から供給された劣化信号振幅スペクトルとを組み合わせて、強調信号（以下の式（４）の左辺）を求める。

　逆フーリエ変換部４３は、得られた強調信号に逆フーリエ変換を施す。逆フーリエ変換された強調信号は、１フレームがＫサンプルを含む時間領域サンプル値系列ｘｎ（ｔ）（ｔ＝０，１，．．．，Ｋ−１）として、窓がけ処理部４２に供給され、窓関数ｗ（ｔ）との乗算が行なわれる。第ｎフレームの入力信号ｘｎ（ｔ）（ｔ＝０，１，．．．，Ｋ／２−１）に対してｗ（ｔ）で窓がけされた信号は、次式（５）の左辺で与えられる。

　また、窓がけ処理部４２は、連続する２フレームの一部を重ね合わせ（オーバラップ）して窓がけしてもよい。フレーム長の５０％をオーバラップ長として仮定すれば、ｔ＝０，１，．．．，Ｋ／２−１に対して、以下の式の左辺が、窓がけ処理部４２の出力となり、フレーム合成部４１に伝達される。

　フレーム合成部４１は、窓がけ処理部４２からの隣接する２フレームの出力を、Ｋ／２サンプルずつ取り出して重ね合わせ、以下の式（７）によって、ｔ＝０，１，．．．，Ｋ−１における出力信号（式（７）の左辺）を得る。得られた出力信号は、フレーム合成部４１から出力端子５に伝達される。

　なお、図２と図３において変換部２と逆変換部４における変換をフーリエ変換として説明した。変換部２、逆変換部４は、フーリエ変換に代えて、コサイン変換、修正コサイン変換、アダマール変換、ハール変換、ウェーブレット変換など、他の変換を用いることもできる。例えば、コサイン変換や修正コサイン変換は、変換結果として振幅だけしか得られない。このため、図１における変換部２から逆変換部４に至る経路は不要になる。また、雑音情報記憶部６に記録する雑音情報も、振幅（又はパワー）だけとなり、記憶容量の削減、雑音抑圧処理における演算量の削減に貢献する。変換部２、逆変換部４がハール変換を用いた場合には、乗算が不要となり、ＬＳＩ化したときの面積を小さくすることができる。変換部２、逆変換部４がウェーブレット変換を用いた場合には、周波数によって時間解像度を異なったものに変更できるために、雑音抑圧効果の向上が期待できる。
　また、変換部２が周波数成分を複数統合してから、衝撃音抑圧部１１が実際の抑圧を行うこともできる。その際、変換部２は、聴覚特性の弁別能力が高い低周波領域から、能力が低い高周波領域に向かって、よりたくさんの周波数成分を統合することにより、高い音質を達成することができる。このように、複数の周波数成分を統合してから衝撃音抑圧を実行すると、雑音抑圧を適用する周波数成分の数が少なくなり、雑音抑圧装置１００は、全体の演算量を削減することができる。
　＜衝撃音抑圧部１１の構成＞
　図４は、衝撃音抑圧部１１の内部構成を示すブロック図である。図４に示すように、衝撃音抑圧部１１は、遅延部１１１と合成部１１２とを含む。遅延部１１１は、入力である劣化信号位相スペクトルを遅延させる。遅延量は一つである必要はなく、遅延部１１１は、複数の遅延量で入力を遅延して、複数の遅延信号を生成することもできる。合成部１１２は、劣化信号位相スペクトルと遅延部１１１から供給される遅延した劣化信号位相スペクトルを用いて、強調信号位相スペクトルを合成する。
　合成部１１２は、位相処理を、衝撃音検出部１０から衝撃音の検出が伝達されたときだけ行う。合成部１１２は、位相処理として、過去（衝撃音発生前）の値を用いて、以下の数式（８）に示す処理を位相に適用することができる。

　ここで、Ｎは周波数の範囲を、Ｍはフレーム番号の範囲を表わし、ｃｌｐは周波数ｐ、ｌフレーム前の劣化信号位相スペクトルに対する係数を表わす。
　すなわち、合成部１１２は、周波数範囲ｋ−Ｎ／２からｋ＋Ｎ／２、フレーム番号ｎからｎ−Ｍ＋１の範囲にある劣化信号位相スペクトルの線形結合で、強調信号の位相を計算する。最も簡単な例が、各周波数における１フレーム前の位相との平均である。或いは、合成部１１２は、１フレーム前の位相と同じ位相を適用（置換）しても良い。これにより、現在の位相そのものよりも過去の位相との差が小さくなり、衝撃音として知覚されにくくなる。この考え方を拡張すると、衝撃音抑圧部１１は、信号全体を遅延させ、衝撃音に続く未来の信号成分の位相を、過去の信号成分の位相と同様に用いて、位相の変化を抑圧することによって、衝撃音抑圧効果を向上させることもできる。この位相処理による衝撃音抑圧効果は非常に大きく、パワー制御や振幅制御を行なわなくても、位相処理だけで、衝撃音抑圧効果を得ることができる。
　さらに、衝撃音抑圧部１１は、位相に過去の値とは無関係な成分を加算することもできる。このような成分の例としては、ランダム位相がある。さらに、衝撃音抑圧部１１は、ランダム位相が４５度以下であるというように、ランダム位相の範囲に制約を加えることもできる。衝撃音抑圧部１１は、位相に過去の値と無関係な成分を付加することで、衝撃音を効果的に抑圧できる。
　以上説明したように、本実施形態にあっては、雑音抑圧装置１００は、劣化信号中の衝撃音を抑圧する際に、劣化信号中において衝撃音を検出し、検出された衝撃音の位相成分を、劣化信号中の衝撃音以外の信号の位相成分を用いて処理した。これにより、雑音抑圧装置１００は、衝撃音を一層効果的に抑圧することができる。
　（第２実施形態）
　次に、図５を用いて本発明の第２実施形態について説明する。本実施形態は衝撃音の検出方法に特徴を有する衝撃音抑圧装置及び方法に関するものである。従来の衝撃音検出方法では十分な検出精度を得られていなかったが、本実施形態のように衝撃音の検出を行なうことで、非常に精度良く劣化信号中の衝撃音を検出できるようになる。
　本実施形態における衝撃音検出部１０は、劣化信号振幅スペクトルの周波数特性と時間特性に基づいて、衝撃音の存在を検出する。衝撃音検出部１０は、検出に当って、周波数特性と時間特性のいずれかを用いてもよいし、双方を用いてもよい。また、衝撃音検出部１０は、双方を用いる際に、それぞれの特性評価結果の重み付き和またはより複雑な関数で表現される統合結果を用いることもできる。
　＜衝撃音検出部の構成＞
　図５は、衝撃音検出部１０の構成を示すブロック図である。図５を参照すると、衝撃音検出部１０は、スペクトル周波数特性評価部１０１、スペクトル時間特性評価部１０２、及び統合部１０３を含む。
　スペクトル周波数特性評価部１０１は、スペクトルの周波数方向変化に関する特性を評価して、統合部１０３に供給する。スペクトル周波数特性評価部１０１は、スペクトルの周波数方向変化に関する特性として、スペクトルの周波数方向平坦度を評価する。スペクトル周波数特性評価部１０１は、スペクトル平坦度として、隣接周波数点におけるスペクトルの差分絶対値の総和を用いることができる。周波数ｋ、フレームｎにおける劣化信号振幅スペクトル｜Ｙｎ（ｋ）｜を用いると、スペクトル周波数特性評価部１０１は、フレームｎにおけるスペクトル平坦度Ｆｍｆ（ｎ）を、次式（９）で求めることができる。

　また、スペクトル周波数特性評価部１０１は、スペクトル平坦度として、平均スペクトルとの差分絶対値の総和を用いることもできる。フレームｎにおける平均劣化信号振幅スペクトル｜Ｙｎ｜バーを用いると、スペクトル周波数特性評価部１０１は、フレームｎにおけるスペクトル平坦度Ｆｍｆ（ｎ）を、次式（１０）で求めることができる。

　スペクトル周波数特性評価部１０１は、ｋによって平坦度計算の対象とする周波数を制限することもできる。特に、衝撃音スペクトルは高域で強く、通常信号のスペクトルは低域で強いので、スペクトル周波数特性評価部１０１は、ｋの範囲を高周波領域に限定することで、検出精度を高くすることができる。また、スペクトル周波数特性評価部１０１は、平坦度を複数のサブバンドで個別に求め、これらの線形または非線形結合によって、全体的な平坦度を求めてもよい。スペクトル周波数特性評価部１０１は、サブバンド処理を、衝撃音と摩擦音の識別に利用することもできる。衝撃音と摩擦音はいずれも広帯域にわたって平坦なスペクトル特性を有するが、一般的に摩擦音の方が狭帯域であり、低域パワーが弱い。このような特性の違いを識別するために、サブバンド処理と複数のサブバンド平坦度の結合は有効である。
　スペクトル周波数特性評価部１０１は、このようにして求めた平坦度を閾値と比較して、平坦度のスコアを求める。平坦度のスコアは、どの程度平坦であるかを表す指標で、例えば、１と０の間に正規化した値として表現することができる。平坦度の上限閾値をσＨ、下限閾値をσＬ、これらに対応した平坦度をＦＨ、ＦＬとすると、スペクトル周波数特性評価部１０１は、平坦度のスコアＳｆ（ｎ）を、式（１１）で定めることができる。

　式（１１）においては、上限閾値と下限閾値の間を線形補間しているが、これは任意の関数または多項式などによる補間を適用することができる。
　閾値は、例えば、予め定められた値の他、平坦度の過去の平均や中央値など、あるいはそれらに基づいて計算した値などである。また、スペクトル周波数特性評価部１０１は、複数の閾値を用意しておき、劣化信号スペクトルを分析した結果に基づいて、選択的に用いることもできる。このような分析結果の例は、劣化信号振幅スペクトル、パワースペクトル、これらの統計量（平均値、中央値、最大値、最小値、分散）などである。
　一方、スペクトル時間特性評価部１０２は、スペクトルの時間方向変化に関する特性を評価して、統合部１０３に供給する。スペクトル時間特性評価部１０２は、スペクトルの時間方向変化として、振幅又はパワースペクトルの増分を用いることができる。スペクトル時間特性評価部１０２は、時間方向変化の評価を、各周波数点で行う。スペクトル時間特性評価部１０２は、これらの線形または非線形結合によって、全体的な変化を求めてもよい。また、スペクトル時間特性評価部１０２は、時間方向変化をサブバンドで求めることもできる。例えば、スペクトル時間特性評価部１０２は、次式（１２）によって、一つのサブバンドにおける時間方向の変化Ｆｍｔ（ｎ）を求めることができる。

　ｋの値によって、サブバンド下限が決定される。また、スペクトル時間特性評価部１０２は、Ｎ−１の代わりに特定の周波数番号を用いれば、サブバンド上限を指定できる。また、スペクトル時間特性評価部１０２は、時間方向変化を複数のサブバンドで個別に求め、これらの線形または非線形結合によって、全体的な変化を求めてもよい。衝撃音スペクトルは高域で強く、通常信号のスペクトルは低域で強いので、スペクトル時間特性評価部１０２は、高周波領域に限定して変化を評価することで、検出精度を高くすることができる。
　また、スペクトル時間特性評価部１０２は、振幅又はパワースペクトルの周波数方向の統計量（平均値、中央値、最大値、最小値、分散）、あるいはそれらのうちいくつかの組合せを用いることもできる。例えば、スペクトル時間特性評価部１０２は、最小値を用いた場合、次式（１３）によって時間変化を求めることができる。

　スペクトル時間特性評価部１０２がこのような最小値の時間変化を用いることで、衝撃音検出部１０は、非常に的確に衝撃音を検出することができる。これは、劣化信号の周波数方向の統計量は、通常、幅広い値を取り得るが、衝撃音に関しては、周波数方向の最小値が大きい傾向があるからである。
　特に、これらの統計量のいずれかが小さな分散を有している場合には、スペクトル時間特性評価部１０２は、分散の小さい統計量を利用することで、検出の精度を高くすることができる。
　なお、数式（９）~（１２）に対しては、劣化信号振幅スペクトル｜Ｙｎ（ｋ）｜に代えてパワースペクトル｜Ｙｎ（ｋ）｜２を用いることもできる。
　スペクトル時間特性評価部１０２は、このようにして求めた時間変化を閾値と比較して、時間変化のスコアを求める。時間変化のスコアは、どの程度の時間変化が存在するかを表す指標で、例えば、１と０の間に正規化した値として表現することができる。時間変化の上限閾値、下限閾値、これらに対応した時間変化量を用いて、スペクトル時間特性評価部１０２は、時間変化のスコアＳｔ（ｎ）を、式（１１）と同様にして定めることができる。スペクトル時間特性評価部１０２は、平坦度のスコアと同じように、線形補間の代わりに任意の関数または多項式などによる補間を適用することができる。
　スペクトル時間特性評価部１０２は、閾値として予め定められた値の他、時間変化の過去の平均や中央値など、あるいはそれらに基づいて計算した値などを用いてもよい。また、スペクトル時間特性評価部１０２は、複数の閾値を用意しておき、劣化信号振幅スペクトルを分析した結果に基づいて、選択的に用いることもできる。このような分析結果の例は、劣化信号振幅スペクトル、パワースペクトル、これらの統計量（平均値、中央値、最大値、最小値、分散）などである。
　統合部１０３は、スペクトル周波数特性評価部１０１から供給されたスペクトルの周波数方向変化に関する特性とスペクトル時間特性評価部１０２から供給されたスペクトルの時間方向変化に関する特性を統合して、衝撃音データを生成し、これを出力する。衝撃音データは、例えば、０と１の間の値に正規化された衝撃音らしさである。例えば、衝撃音データが１であることは、１００％の確信をもって衝撃音と判定することであり、０．８であることは、２０％の不確定性をもって衝撃音と判定することである。
　前記特性の統合の最も簡単な方法は、前記平坦度のスコアと時間変化のスコアの論理積である。統合部１０３は、２つのスコアの双方が１に等しいときに衝撃音データを１とする。また、統合部１０３は、論理積の代わりに論理和を用いることもできる。統合部１０３は、いずれか一方のスコアが１であれば、衝撃音データを１とする。
　統合部１０３は、これらのスコアを統合した統合スコアを用いて衝撃音データを計算することもできる。例えば、統合部１０３は、両者の和を衝撃音データとすると、論理積や論理和より不確実な場合も、衝撃音を１以上に設定することができる。スコアの統合は、両者の単純な和だけでなく、線形、非線形関数を含む様々な統合でもよい。この統合に用いる関数によって、周波数特性または時間特性のどちらをどの程度重視するのかを調整することができる。
　このようにして求めた衝撃音データが１以上であれば、雑音抑圧装置は、衝撃音が確実に存在すると判定して、衝撃音を完全に抑圧する。衝撃音データが１以下であれば、雑音抑圧装置は、その値に応じて衝撃音抑圧の程度を弱くする。
　以上のように、本実施形態では、雑音抑圧装置は、劣化信号中の衝撃音を抑圧する際に、劣化信号から振幅又はパワー成分を抽出し、その振幅又はパワー成分の時間方向変化の統計量を用いて、衝撃音を検出する。これにより、雑音抑圧装置は、より正確に衝撃音を検出することが可能となる。
　なお、本実施形態では、第１実施形態の一部としての衝撃音検出部１０について説明したが、本実施形態の衝撃音検出方法は、第１実施形態に記載の衝撃音抑圧方法に限定されるものではなく、どのように衝撃音を抑圧するかについては問わない。すなわち、雑音抑圧装置は、本実施形態の方法で検出した衝撃音に対し、第１実施形態で説明したように、位相処理によって衝撃音を抑圧してもよいし、振幅やパワーを制御することによって衝撃音を抑圧してもよい。
　（第３実施形態）
　ここで、本発明の第３実施形態としての雑音抑圧装置について説明する。図６は、本実施形態に係る雑音抑圧装置３００を示す図である。雑音抑圧装置３００は、第１衝撃音抑圧部１１と第２衝撃音抑圧部１２とを備える。衝撃音検出部１０と第２衝撃音抑圧部１２には、変換部２から劣化信号振幅スペクトル｜Ｙｎ（ｋ）｜が供給される。逆変換部４には、第１衝撃音抑圧部から強調信号位相スペクトルａｒｇ　Ｘｎ（ｋ）バーが、第２衝撃音抑圧部１２から強調信号振幅スペクトル｜Ｘｎ（ｋ）バー｜が供給される。第１衝撃音抑圧部１１の構成は、第１実施形態で説明した衝撃音抑圧部と同じであるため、ここでは詳しい説明を省略する。
　図７は、第２衝撃音抑圧部１２の内部構成を示すブロック図である。図７に示すように、第２衝撃音抑圧部１２は、遅延部１２１と合成部１２２とを含む。遅延部１２１は、入力である劣化信号振幅スペクトルを遅延させる。遅延量は一つである必要はない。遅延部１２１は、複数の遅延量で入力を遅延して、複数の遅延信号を生成することもできる。合成部１２２は、入力した劣化信号振幅スペクトルと遅延部１２１から供給される遅延した劣化信号振幅スペクトルとを合成して、強調信号振幅スペクトルを生成する。合成部１２２は、遅延信号との合成処理を、衝撃音検出部１０において衝撃音が検出されたときだけ行う。
　合成部１２２は、合成処理として、過去の値を用いて、第１実施形態の数式（８）に示す処理を適用したように、次式（１４）に示す処理を適用することができる。

　数式（１４）においてｃｌｐは、周波数ｐ、ｌフレーム遅延した劣化信号振幅スペクトルに対する係数を表わす。すなわち、合成部１２２は、周波数範囲がｋ−Ｎ／２からｋ＋Ｎ／２、フレーム番号ｎからｎ−Ｍ＋１の範囲にある劣化信号振幅スペクトの線形結合で、強調信号振幅スペクトルを計算する。最も簡単な例が、各周波数における１フレーム前のサンプルとの平均である。平均を用いることで、現在のサンプル単独の場合よりも過去のサンプルとの差が小さくなり、衝撃音として知覚されにくくなる。
　合成の別の例として、合成部１２２は、過去のサンプルから求めた値（例えば平均値、最大値）を上限として、現在のサンプルに制限をかけてもよい。この合成法も、現在のサンプル単独の場合よりも過去のサンプルとの差を小さくし、衝撃音として知覚されにくくする。また、第２衝撃音抑圧部１２は、位相に対する処理と同様に、信号全体を遅延させ、衝撃音に続く未来の信号成分の振幅スペクトルを、過去の信号成分の振幅スペクトルと同様に用いて、振幅スペクトルの変化を抑圧することによって、衝撃音抑圧効果を向上させることもできる。なお、これらの合成において、劣化信号振幅スペクトルに代えて劣化信号パワースペクトルを利用できることは、すでに他の説明で述べた通りである。
　以上のように本実施形態は、検出された衝撃音の振幅又はパワー成分が小さくなるように、劣化信号中の衝撃音以外の信号を用いて処理する。このように位相と振幅又はパワーの両面から衝撃音を処理することで、雑音抑圧装置３００は、より一層効果的に衝撃音を抑圧することが可能となる。
　（第４実施形態）
　次に、本発明の第４実施形態としての雑音抑圧装置４００について図８を用いて説明する。本実施形態の雑音抑圧装置４００は、第１実施形態の雑音抑圧装置１００に加えて、雑音存在情報を入力する入力端子９を有する。図８における衝撃音抑圧部３１は、入力端子９から供給された雑音存在情報を用いて、各周波数で、第１実施形態で説明した位相処理を行なうことにより衝撃音を抑圧する。そして、衝撃音抑圧部３１は、衝撃音抑圧結果を強調信号スペクトルとして逆変換部４に伝達する。強調信号位相スペクトルは、雑音存在情報が雑音の存在を示すときには劣化信号位相スペクトルに対して第１実施形態で説明した位相処理を施し衝撃音を抑圧したものとなる。一方、雑音存在情報が雑音の不存在を示すときには劣化信号位相スペクトルそのものとなる。
　これにより、より効率の良い衝撃音抑圧を行なうことができる。
　（第５実施形態）
　次に本発明の第５実施形態としての雑音抑圧装置について説明する。本実施形態は、図６を用いて説明した第３実施形態に係る雑音抑圧装置をベースにしたものであり、第３実施形態とは、第２衝撃音抑圧部１２の内部構成が異なるものである。その他の構成及び動作は、第３実施形態と同様であるため、ここでは詳しい説明を省略する。
　本実施形態に係る第２衝撃音抑圧部５２の内部構成について、図９に示す。図９は、衝撃音抑圧部５２の構成を示すブロック図である。図９に示すように、第２衝撃音抑圧部５２は、遅延部１２１、合成部１２４に加えて、重要度評価部１２３を備えている。遅延部１２１については、第３実施形態において図７を用いて説明した通りの構成であるためここでは説明を省略する。
　重要度評価部１２３は、重要度に応じた処理を遂行するための情報（重要度情報）を生成し、合成部１２４に供給する。合成部１２４は、強調信号スペクトル合成処理に加えて、重要度評価部１２３から供給された重要度情報に従って、重要度に応じた処理を行う。
　重要度評価部１２３の生成する重要度情報の第１の例は、劣化信号振幅スペクトルのピークである。重要度評価部１２３は、スペクトルのピークを、各周波数点におけるスペクトルを隣接する周波数点におけるスペクトルと比較して、十分に大きいかどうかを評価して検出することができる。最も簡単な例は、各周波数点のスペクトルをその両隣（低域側及び高域側）のスペクトルと比較して、その差が閾値より大きいときにピークと判定することである。差に対する閾値は、両側のスペクトルに対して等しい必要はない。日本工業規格　ＪＩＳ×４３３２−３「音響映像オブジェクトの符号化　—第３部　音響—」、２００２年３月には、高域側の差分閾値を低域側差分閾値よりも小さくすることが、聴覚特性に合致すると記載されている。同様にして、重要度評価部１２３は、低域側及び高域側の複数の周波数点に対して差分を求め、これらの情報を総合してピークを検出することもできる。すなわち、すぐ隣の周波数点に対しては差分が大きいが、それよりも離れた隣接周波数点同志においては差分が小さい周波数点を検出すれば、それがピークとなる。重要度評価部１２３は、このようにして検出したピークの位置（周波数）と大きさ（重要度）を、合成部１２４へ供給する。
　重要度評価部１２３の生成する重要度情報の第２の例は、劣化信号振幅スペクトルの大きさである。重要度評価部１２３は、スペクトルがピークを形成しなくても、その値が大きいときには、その周波数を大振幅として検出する。例えば、大きな値のスペクトルが周波数方向に連続すると、ピークとしては検出されない。しかし、このような部分は、聴覚にとって重要である。そこで、重要度評価部１２３は、検出した大振幅の位置（周波数）と大きさ（重要度）を、合成部１２４へ供給する。
　重要度評価部１２３の生成する重要度情報の第３の例は、劣化信号振幅スペクトルの雑音らしさである。前述のピーク検出を行い、検出されたピークのうち、特に低域に存在するピークは雑音の可能性が低い。また、スペクトル値が小さくピークでない位置では雑音らしさが高い。すなわち、ピークは雑音らしさが低く、非ピークでスペクトル値が小さいときは雑音らしさが高い。重要度評価部１２３は、これらのピークの位置（周波数）と大きさ（重要度）を、合成部１２４へ供給する。
　重要度評価部１２３の生成する重要度情報は、既に説明したピーク、大振幅、及び雑音らしさを適切に組み合わせてもよい。例えば、大振幅のスペクトルに対してピーク検出の閾値を低くして、振幅が大きい帯域では小さなピークも検出されるように制御することなどがその例である。重要度評価部１２３は、指標を組み合わせて用いることで、より正確な重要度情報を得ることができる。また、これまでの他の説明のように、重要度評価部１２３は、処理を特定の周波数帯域に限定する、サブバンド処理などを適用することが可能である。
　具体的に合成部１２４は、重要度評価部１２３から供給された周波数点以外では、図７で説明した合成部１２２と同じ強調信号スペクトル合成処理を行う。重要度評価部１２３から供給された周波数点には重要な信号成分が存在し、これらは強調信号の音質に重要な役割を果たす。そこで、衝撃音抑圧部５２は、これらの周波数点においてはその重要度に応じた抑圧を適用する。つまり、衝撃音抑圧部５２は、重要度が高い場合には弱い抑圧を、重要度が低い場合には強い抑圧を適用する。
　以上、本実施形態によれば、雑音の振幅又はパワースペクトルに対して、重要度を加味した抑圧を行なうことができ、より高品質の出力を得ることができる。
　（第６実施形態）
　次に本発明の第６実施形態としての雑音抑圧装置について説明する。本実施形態は、図６を用いて説明した第３実施形態に係る雑音抑圧装置をベースにしたものであり、第３実施形態とは、第２衝撃音抑圧部１２の内部構成が異なるものである。その他の構成及び動作は、第３実施形態と同様であるため、ここでは詳しい説明を省略する。
　図１０Ａは、本実施形態に係る雑音抑圧装置の全体構成図である。図６の構成によく似ているが、第２衝撃音抑圧部６２が、入力端子９から雑音存在情報が供給される点で異なる。その他の構成及び動作については、図６と同様であるためここでは詳細な説明を省略する。
　図１０Ｂは、第２衝撃音抑圧部６２の内部構成を示すブロック図である。図１０Ｂに示すように、衝撃音抑圧部６２は、遅延部１２１、合成部１３４及び背景音推定部１２５を含んでいる。遅延部１２１については、図７で説明したものと同様であるのでここでは説明を省略する。背景音推定部１２５は、劣化信号振幅スペクトルを変換部２から受け、雑音存在情報を入力端子９から受けて背景音レベルを推定し、背景音レベル推定値として合成部１３４に供給する。背景音推定部１２５は、背景音レベル推定値を、入力として劣化信号振幅スペクトルが供給されるときには背景音振幅スペクトルの推定値、劣化信号パワースペクトルが供給されるときには背景音パワースペクトルの推定値として求める。背景音推定部１２５は、背景音の推定を、雑音存在情報によって雑音が存在するとされるときだけ行い、背景音の推定値を更新する。合成部１３４は、合成部１２２と同じ強調信号スペクトル合成処理に加えて、背景音推定部１２５から供給された背景音推定値によって異なる処理を行う。
　合成部１３４は、入力端子９から供給される雑音存在情報が雑音の存在を表すときに、背景音推定部１２５から供給された背景音推定値を下限値とした抑圧を行う。すなわち、合成部１３４は、合成した結果が背景音推定値より小さいときには、背景音推定値と等しくなるまで抑圧を弱くし、それを強調信号スペクトルとして出力する。合成部１３４は、合成結果が背景音推定値に等しいかそれよりも大きいときは、合成結果をそのまま強調信号スペクトルとして出力する。入力端子９から供給される雑音存在情報が雑音の非存在を表すときには、合成部１３４は背景音推定値を下限値とした処理を行わず、合成結果をそのまま強調信号スペクトルとして出力する。
　以上のように、衝撃音抑圧部６２は、背景音推定値を下限値とした抑圧を行うことにより、過剰抑圧を回避し、自然な聴感を与える強調信号を得ることができる。
　（第７実施形態）
　次に、本発明の第７実施形態としての雑音抑圧装置について図１１を用いて説明する。本実施形態は、図６を用いて説明した第３実施形態に係る雑音抑圧装置をベースにしたものであり、第３実施形態とは、第２衝撃音抑圧部７２の内部構成が異なるものである。その他の構成及び動作は、第３実施形態と同様であるため、ここでは詳しい説明を省略する。
　図１１は、第２衝撃音抑圧部７２の内部構成を示すブロック図である。図１１に示すように、第２衝撃音抑圧部７２は、遅延部１２１、合成部１２２及び白色化処理部１２７を含んでいる。遅延部１２１と合成部１２２との関係については、図５から図７で説明した通りであるためここでは説明を省略する。白色化処理部１２７は、合成部１２２から強調信号スペクトルを供給され、これを白色化して、白色化強調信号スペクトルとして出力する。
　白色化処理部１２７は、強調信号振幅スペクトルの平均値を求め、この平均値からの分散を基準値以下にする。具体的には、白色化処理部１２７は、平均値＋εを超える振幅スペクトル値を平均値＋εで置き換える。また、白色化処理部１２７は、平均値−εより小さい振幅スペクトル値を平均値−εで置き換える。白色化処理部１２７は、それ以外の強調信号振幅スペクトルについては、そのままとする。また、白色化処理部１２７は、平均値±εで置き換える代わりに、平均値±εの範囲の乱数で置き換えても良い。例えば、白色化処理部１２７は、平均値＋εを超える振幅スペクトル値は平均値＋εと平均値の間の乱数で置き換える。また、白色化処理部１２７は、平均値−εより小さい振幅スペクトル値は平均値−εと平均値の間の乱数で置き換える。白色化処理によって振幅スペクトルの値が均一化し、雑音が知覚されにくくなる。
　さらに、図１１の構成に加えて、図９を用いて説明した重要度評価部１２３を設けても良い。その場合、白色化処理部１２７は、白色化処理に重要度評価部１２３の出力を利用することもできる。重要度評価部１２３が雑音らしさを求め、雑音らしさが大きいときだけ、白色化処理部１２７が白色化処理を適用する。このようにすることで、所望信号成分が少ないときに強調信号が白色信号に近くなり、雑音として知覚されにくくなる。
　白色化処理部１２７は、これら白色化処理において、複数のサブバンドで個別の処理を行うこともできる。また、白色化処理部１２７は、特定のサブバンドで白色化を回避することもできる。この場合、白色化処理部１２７は、サブバンド毎に異なる平均値を用いるので、自然な聴感を与える強調信号を得ることができる。
　（第８実施形態）
　図１２は、本発明の第８実施形態に係る雑音抑圧装置の構成を示すブロック図である。本実施形態を、第１実施形態と比べた場合、雑音抑圧部３が追加されている他は、同様の構成である。したがって、同じ構成については同じ符号を付してその説明を省略する。
　雑音抑圧部３は、変換部２から供給された劣化信号振幅スペクトルと入力した雑音情報（外部から供給される雑音に関する情報）とを用いて、各周波数で雑音を抑圧し、雑音抑圧結果としての強調信号振幅スペクトルを逆変換部４に伝達する。
　以上の構成により、雑音抑圧装置は、衝撃音以外のノイズをも的確に抑圧することが可能となる。
　（第９実施形態）
　図１３は、本発明の第９実施形態に係る雑音抑圧装置の構成を示すブロック図である。本実施形態を、第８実施形態と比べた場合、衝撃音検出部９０が、雑音抑圧部３によって雑音を抑圧された結果を用いて衝撃音の検出を行なう点で異なる。その他は、同様の構成であるため、同じ構成については同じ符号を付してその説明を省略する。
　雑音抑圧部３からの出力は衝撃音検出部９０に入力される。衝撃音検出部９０の構成は、第１実施形態で説明した衝撃音検出部１０と同じであるため、ここでは詳しい説明を省略する。
　以上の構成により、衝撃音検出部９０は、雑音抑圧部３によって雑音を抑圧された結果を用いて、より正確に衝撃音を検出することが可能となる。
　（第１０実施形態）
　図１４は、本発明の第１０実施形態に係る雑音抑圧装置の構成を示すブロック図である。本実施形態を、第８実施形態と比べた場合、衝撃音検出部９１が、雑音情報を用いて衝撃音の検出を行なう点で異なる。衝撃音検出部９１は、供給された雑音情報（例えば、雑音の存在を示す情報（雑音存在情報）、スペクトル形状等に関する情報を含む雑音情報）を利用して、衝撃音の検出を行う。その他は、同様の構成であるため、同じ構成については同じ符号を付してその説明を省略する。
　衝撃音検出部９１は、変換部２から供給された劣化信号振幅スペクトルと入力した雑音情報とを用いて、雑音情報が雑音の存在を表すときに衝撃音を検出する。
　以上の構成により、雑音抑圧装置は、正確に衝撃音を検出し、これを抑圧することが可能となる。
　（他の実施形態）
　以上説明してきた第１乃至第１０実施形態では、それぞれ別々の特徴を持つ雑音抑圧装置について説明したが、それらの特徴を如何様に組み合わせた雑音抑圧装置も、本発明の範疇に含まれる。
　また、本発明は、複数の機器から構成されるシステムに適用しても良いし、単体の装置に適用しても良い。さらに、本発明は、実施形態の機能を実現するソフトウェアの信号処理プログラムが、システム或いは装置に直接或いは遠隔から供給される場合にも適用可能である。したがって、本発明の機能をコンピュータで実現するために、コンピュータにインストールされるプログラム、或いはそのプログラムを格納した媒体、そのプログラムをダウンロードさせるＷＷＷサーバも、本発明の範疇に含まれる。
　図１５は、第１実施形態を信号処理プログラムにより構成する場合に、その信号処理プログラムを実行するコンピュータ１１００の構成図である。コンピュータ１１００は、入力部１１０１と、ＣＰＵ１１０２と、出力部１１０３と、メモリ１１０４と、通信制御部１１０６とを含む。
　ＣＰＵ１１０２は、信号処理プログラムを読み込むことにより、コンピュータ１１００の動作を制御する。すなわち、信号処理プログラムを実行したＣＰＵ１１０２は、劣化信号中から衝撃音を検出する（Ｓ８０１）。次に、ＣＰＵ１１０２は、劣化信号中において検出した衝撃音の位相情報を、衝撃音以外の信号の位相情報を用いて処理する（Ｓ８０２）。
　これにより、第１実施形態と同様の効果を得ることができる。
　以上、実施形態を参照して本願発明を説明したが、本願発明は以上の実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で同業者が理解し得る様々な変更をすることができる。
　この出願は、２００９年１１月９日に出願された日本出願特願２００９−２５６５９６を基礎とする優先権を主張し、その開示の全てをここに取り込む。

Claims (20)

1. 　劣化信号中の衝撃音を抑圧するために、
   　前記劣化信号中において衝撃音を検出し、
   　検出された衝撃音の位相情報を、該位相情報の変化量が小さくなるように前記劣化信号中の前記衝撃音以外の信号の位相情報で処理する信号処理方法。
2. 　前記衝撃音の位相情報を、前記劣化信号中の衝撃音発生前の信号の位相情報で処理する請求項１に記載の信号処理方法。
3. 　前記衝撃音の位相情報を、前記劣化信号中の衝撃音発生前の信号の位相情報で置換する請求項２に記載の信号処理方法。
4. 　前記衝撃音の位相情報を、前記劣化信号中の衝撃音発生前の信号の位相情報と、前記衝撃音の位相情報との平均値で置換する請求項２に記載の信号処理方法。
5. 　前記劣化信号を遅延させ、前記衝撃音の位相情報を、前記劣化信号中の衝撃音発生前の信号の位相情報で処理する請求項２乃至４の何れか１項に記載の信号処理方法。
6. 　前記劣化信号を、前記位相情報と振幅又はパワー情報とに分解し、
   　前記振幅又はパワー情報を用いて、前記劣化信号中の衝撃音を検出する請求項１乃至５の何れか１項に記載の信号処理方法。
7. 　検出された衝撃音の前記振幅又はパワー情報が小さくなるように、前記劣化信号中の前記衝撃音以外の信号の前記振幅又はパワー情報を用いて処理する請求項６に記載の信号処理方法。
8. 　検出された衝撃音の前記振幅又はパワー情報を、前記劣化信号中の前記衝撃音が発生する前の信号の振幅又はパワー情報と結合させる請求項７に記載の信号処理方法。
9. 　検出された衝撃音の前記振幅又はパワー情報を、前記劣化信号中の前記衝撃音が発生する前の信号の振幅又はパワー情報を用いて平均化する請求項７に記載の信号処理方法。
10. 　前記劣化信号中の前記衝撃音が発生する前の信号の振幅又はパワー情報を用いて、検出された衝撃音の前記振幅又はパワー情報に制限をかける請求項６に記載の信号処理方法。
11. 　前記劣化信号を遅延させ、前記衝撃音の振幅又はパワー情報を、前記劣化信号中の衝撃音発生前後の信号の振幅又はパワー情報で処理する請求項８乃至１０の何れか１項に記載の信号処理方法。
12. 　雑音存在情報を入力し、該雑音存在情報が雑音の存在を示すときに前記衝撃音の抑圧を行なう請求項１乃至１１の何れか１項に記載の信号処理方法。
13. 　前記劣化信号中の重要度を評価し、
    　前記劣化信号中の重要度が高い部分は弱く、それ以外は強く前記衝撃音を抑圧する請求項１乃至１２の何れか１項に記載の信号処理方法。
14. 　前記劣化信号中の背景音を推定し、
    　前記劣化信号中の背景音推定値を下限値として、前記衝撃音を抑圧する請求項１乃至１１の何れか１項に記載の信号処理方法。
15. 　前記振幅又はパワー情報の平均値を求め、該平均値からの分散を基準値以下にする請求項６に記載の信号処理方法。
16. 　前記振幅又はパワー情報中の雑音を、雑音情報を用いて抑圧し、その結果を用いて、前記衝撃音を検出する請求項６に記載の信号処理方法。
17. 　前記振幅又はパワー情報中の雑音を、雑音情報を用いて抑圧し、
    　前記雑音情報を用いて、前記衝撃音を検出する請求項６に記載の信号処理方法。
18. 　入力信号を周波数領域信号に変換し、
    　該周波数領域信号を用いて衝撃音を検出し、
    　前記衝撃音が検出されたときに、
    　振幅と位相の変化分が小さくなるように衝撃音を抑圧することを特徴とする雑音抑圧の方法。
19. 　劣化信号中の衝撃音を抑圧する情報処理装置であって、
    　前記劣化信号中において衝撃音を検出する検出手段と、
    　検出された衝撃音の位相情報を、該位相情報の変化量が小さくなるように前記劣化信号中の前記衝撃音以外の信号の位相情報で処理する位相処理手段と、
    　を備える情報処理装置。
20. 　劣化信号中の衝撃音を抑圧する信号処理プログラムであって、
    　コンピュータに、
    　前記劣化信号中において衝撃音を検出する工程と、
    　検出された衝撃音の位相情報を、該位相情報の変化量が小さくなるように前記劣化信号中の前記衝撃音以外の信号の位相情報を用いて処理する工程と、
    　を実行させる信号処理プログラムを格納するプログラム記録媒体。

Images