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

Description

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

劣化信号（所望の信号と雑音とが混合された信号）から、雑音の一部又は全部を抑圧し、強調信号（所望の信号を強調した信号）を出力する信号処理技術として、雑音抑圧技術（ｎｏｉｓｅ ｓｕｐｐｒｅｓｓｉｎｇ ｔｅｃｈｎｏｌｏｇｙ）が知られている。例えば、ノイズサプレッサは、所望の音声信号に重畳されている雑音（ノイズ）を抑圧するシステムであり、携帯電話など様々な音声端末において利用されている。
この種の技術に関し、特許文献１には、入力信号に１より小さな抑圧係数を乗算することによって、ノイズを抑圧する方法が開示されており、特許文献２には、推定された雑音を劣化信号から直接減算することによって、雑音を抑圧する方法が開示されている。
特許文献１及び２に記載の技術は、既に雑音が混合されて劣化している所望信号から、雑音を推定しなければならない。しかし、劣化信号だけから正確に雑音を推定することには限界があり、特許文献１及び２に記載された方法は、一般的に、雑音が所望信号に対して十分小さい場合のみ有効である。雑音が所望信号に対して十分に小さいという条件が満たされない場合は、雑音推定値の精度が低いため、特許文献１及び２に記載された方法では、十分な雑音抑圧の効果が得られず、さらに強調信号に大きな歪が含まれていた。
一方、特許文献３には、所望信号に混入する雑音の特性が事前にある程度わかる場合に、事前に記録しておいた雑音情報（雑音の特性に関する情報）を、劣化信号から減算することで、雑音を抑圧する技術が開示されている。またここでは、入力信号を分析して得られた入力信号パワーが大きいときは大きな係数を、その入力信号パワーが小さいときは小さな係数を、雑音情報に積算して、その積算結果を劣化信号から減算する方法も開示されている。

特許第４２８２２２７号公報 特開平８−２２１０９２号公報 特開２００６−２７９１８５号公報

しかしながら、未知の特性を有する雑音と既知の特性を有する雑音とを両方抑圧する技術はなかった。例えば、携帯電話を操作しながら音声通話をおこなう場合や、インタビューの録画を行なう場合など、人の話す音声を強調して、他の機械ノイズや背景ノイズなどを抑圧することができなかった。
以上を踏まえ、本発明は、上述の課題を解決する信号処理技術を提供することを目的とする。

上記目的を達成するため、本発明に係る装置は、入力された劣化信号に含まれた雑音を抑圧する雑音抑圧手段を備え、前記雑音抑圧手段は、既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する第１出力手段と、未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力する第２出力手段と、を含み、前記既知雑音情報と前記未知雑音情報とを用いて雑音抑圧処理を行なう。
上記目的を達成するため、本発明に係る方法は、既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力し、未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力し、前記既知雑音情報と前記未知雑音情報とを用いて、入力された劣化信号に含まれた雑音を抑圧する。
上記目的を達成するため、プログラム記録媒体は、既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する第１出力ステップと、未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力する第２出力ステップと、前記既知雑音情報と前記未知雑音情報とを用いて、入力された劣化信号に含まれた雑音を抑圧する雑音抑圧ステップと、をコンピュータに実行させる信号処理プログラムを格納する。

本発明は、未知の特性を有する雑音と既知の特性を有する雑音とを両方抑圧することができる信号処理技術を提供する。

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

以下に、図面を参照して、本発明の実施の形態について例示的に詳しく説明する。ただし、以下の実施の形態に記載されている構成要素はあくまで例示であり、本発明の技術範囲をそれらのみに限定する趣旨のものではない。
（第１実施形態）
本発明に係る情報処理装置の第１実施形態として、劣化信号（所望の信号と雑音とが混合された信号）から、雑音の一部又は全部を抑圧し、強調信号（所望の信号を強調した信号）を出力する装置について説明する。
本実施形態では、携帯電話で音声通話しつつ撮影または機械操作する場合やＩＣレコーダで音声録音中に機械操作を行なう場合や、インタビュー録画中に機械を操作したり機械が動作したりする場合など、を前提としている。そして、本実施形態によれば、装置周辺で発生した雑音（使用者や話者の周囲にある背景ノイズ）の抑圧と、装置自体で発生した雑音（オートフォーカス駆動ノイズ、ズーミングノイズ、シャッター音など既知のノイズ）の抑圧とをバランスよく両立させることができる。
情報処理装置としては、例えばデジタルカメラ、デジタルビデオカメラ、ノートパソコン、携帯電話、ＩＣレコーダ、テレビ会議システムなどといった装置が挙げられるが、本発明はこれに限定されるものではなく、ノイズ除去を要求されるあらゆる装置に適用可能である。
図１は、情報処理装置に含まれる雑音抑圧部３の構成を示すブロック図である。図に示すように、雑音抑圧部３は、既知雑音情報出力部３０１と、未知雑音情報出力部３０３とを含み、既知雑音情報出力部３０１から出力された既知雑音情報と未知雑音情報出力部３０３から出力された未知雑音情報とを用いて雑音抑圧処理を行なう。また、図１６は、情報処理装置Ａの構成を示す他のブロック図である。該情報処理装置Ａは、雑音抑圧部３を備える。雑音抑圧部３は、既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する第１出力部３０１（既知雑音情報出力部３０１）と、未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して未知雑音情報を出力する第２出力部３０３（未知雑音情報出力部３０３）とを備える。以下では図１を用いて説明する。
既知雑音情報出力部３０１は、既知の性質を有する既知雑音を抑圧するため、記憶部３１１に予め記憶された雑音情報を用いて既知雑音情報を出力する。記憶部３１１は、半導体メモリなどの記憶素子を含み、抑圧対象としての既知の雑音の特性に関する情報（雑音情報）を記憶している。抑圧対象として記憶される既知の雑音は、例えば、シャッター音、モータ駆動音、ズーム音、オートフォーカスのフォーカシングノイズ（カチカチという音）等である。一方、未知雑音情報出力部３０３は、未知の性質を有する未知雑音を抑圧するため、劣化信号に含まれる未知雑音を推定する推定部３３１を含む。
以上の構成により、本実施形態の情報処理装置は、未知の特性を有する雑音と既知の特性を有する雑音とを両方抑圧することができる。
（第２実施形態）
図２は、本発明の第２実施形態としての情報処理装置に含まれる雑音抑圧部２３の構成を示すブロック図である。第１実施形態と同様に、雑音抑圧部２３は、既知雑音情報出力部２３１と、未知雑音情報出力部２３３とを含み、既知雑音情報出力部２３１から出力された既知雑音情報と未知雑音情報出力部２３３から出力された未知雑音情報とを用いて雑音抑圧処理を行なう。また、図示していないが、既知雑音情報出力部２３１は、予め雑音情報を記憶する記憶部を含み、未知雑音情報出力部２３３は、劣化信号に含まれる未知雑音を推定する推定部を含む。
ただし、本実施形態の雑音抑圧部２３は、既知雑音情報出力部２３１から出力された既知雑音情報を用いて雑音抑圧処理を行なった信号に対して、未知雑音情報出力部２３３から出力された未知雑音情報を用いて更に雑音抑圧処理を加える。つまり、既知雑音情報出力部２３１から出力された既知雑音情報は、既知雑音抑圧部２３２に供給され、入力された劣化信号に含まれる既知雑音の抑圧に用いられる。更に、既知雑音を抑圧済みの信号は、未知雑音情報出力部２３３に供給され、未知雑音情報出力部２３３は、既に既知雑音の抑圧処理が施された信号を用いて、未知雑音の推定を行なう。そして、推定の結果として出力された未知雑音情報は、未知雑音抑圧部２３４に供給され、既知雑音を抑圧済みの信号に対して推定された未知雑音の抑圧処理が加えられて強調信号として出力される。
以上の構成により、本実施形態の情報処理装置は、未知の特性を有する雑音と既知の特性を有する雑音との両方抑圧を実現する。
（第３実施形態）
図３は、本発明の第２実施形態としての情報処理装置に含まれる雑音抑圧部３３の構成を示すブロック図である。本実施形態に係る雑音抑圧部３３の構成及び動作は、第２実施形態の雑音抑圧部２３とほぼ同様であるが、入力された劣化信号が未知雑音情報出力部３３３に供給される点で異なる。
本実施形態では、未知雑音情報出力部３３３は、入力された劣化信号を用いて、未知雑音の推定を行なう。そして、推定の結果として出力された未知雑音情報は、未知雑音抑圧部２３４に供給され、既知雑音を抑圧済みの信号に対して推定された未知雑音の抑圧処理が加えられて強調信号として出力される。
以上の構成により、既知雑音情報出力部２３１及び既知雑音抑圧部２３２の雑音抑圧精度が低い場合にでも、未知雑音の推定を高い精度で行なうことが可能となる。
（第４実施形態）
図４は、本発明の第４実施形態としての情報処理装置に含まれる雑音抑圧部４３の構成を示すブロック図である。本実施形態に係る雑音抑圧部４３の構成及び動作は、第３実施形態の雑音抑圧部３３とほぼ同様であるが、入力された劣化信号が既知雑音情報出力部４３１にも供給される点で異なる。
本実施形態では、既知雑音情報出力部３３１は、入力された劣化信号を用いて、記憶部に記憶した雑音情報から既知雑音情報を生成する。例えば、既知雑音情報出力部４３１は、入力した劣化信号を分析し、分析結果に応じた混合方法によって予め記憶された雑音情報を混合して、既知雑音情報としての混合雑音情報（擬似雑音情報）を生成し出力する。混合対象となる複数の雑音情報のうち少なくとも１つは予め記憶部に記憶されたものである。
より具体的な例を示すと、記憶部に予め複数の雑音情報が記憶されており、既知雑音情報出力部４３１は、それらを組み合わせて混合する。混合対象は、例えば、既知の雑音の最大雑音情報と平均雑音情報の組合せ、最大雑音情報と平均雑音情報と最小雑音情報の組合せ、ピーク成分雑音情報とそれ以外の組合せ、衝撃成分雑音情報とそれ以外の組合せ等である。
なお、平均雑音情報、最大雑音情報、最小雑音情報、ピーク成分雑音情報、衝撃成分雑音情報の定義は以下の通りである。
平均雑音情報：既知の雑音の全体（複数フレーム）についてフーリエ変換により導きだされる複数のスペクトルについて同じ周波数成分の振幅（またはパワー）を平均したもの。いわゆる時間方向に平均した平均スペクトル。
最大雑音情報：既知の雑音の全体（複数フレーム）についてフーリエ変換により導きだされる複数のスペクトルの周波数成分ごとの振幅（またはパワー）の最大値、いわゆる最大スペクトル。
最小雑音情報：既知の雑音の全体（複数フレーム）についてフーリエ変換により導きだされる複数のスペクトルの周波数成分ごとの振幅（またはパワー）の最小値、いわゆる最小スペクトル。
ピーク成分雑音情報：既知の雑音の全体（複数フレーム）についてフーリエ変換により導きだされるスペクトルにおいて、周波数方向に順に振幅を比較したときに、近傍の値に比べて突出して大きい値を有する周波数成分。
衝撃成分雑音情報：衝撃音の全体についてフーリエ変換により導きだされる複数のスペクトルの平均、いわゆる衝撃音の平均スペクトル。衝撃音自体は、フーリエ変換前の音声信号の時間変化を観測したときに、極めて短時間だけ大きな値を持つものであるが、フーリエ変換後のスペクトルは、所定の周波数帯域に渡って振幅がほぼ一定という特徴を有する。
以上のように、入力した劣化信号の分析結果に応じた混合方法によって予め記憶された雑音情報を混合して、既知雑音情報としての混合雑音情報（擬似雑音情報）を生成すれば、変動の激しい信号特性に対応した雑音抑圧を実現することができる。
（第５実施形態）
図５は、本発明の第５実施形態としての情報処理装置に含まれる雑音抑圧部５３の構成を示すブロック図である。本実施形態に係る雑音抑圧部５３は、既知雑音情報出力部２３１から出力された既知雑音情報と、未知雑音情報出力部３３３から出力された未知雑音とを選択的に用いて雑音抑圧処理を行なう点で、上記第２乃至第４実施形態と異なる。既知雑音情報出力部２３１の構成及び動作は、第２実施形態と同様であり、未知雑音情報出力部３３３の構成及び動作は、第３実施形態と同様であるため同じ符号を用いている。
雑音抑圧部５３において選択部５３５は、例えば、既知雑音情報に関する利得ｇ１（０＜ｇ１＜１）と、未知雑音情報に関する利得ｇ２（０＜ｇ２＜１）とをそれぞれ算出し、何れか小さい方または大きい方を選択して抑圧部５３６に供給する。抑圧部５３６は、選択部５３５から供給された利得を劣化信号に積算することにより雑音抑圧処理を行なう。ここで、雑音抑圧部５３の雑音抑圧処理は、劣化信号の周波数成分ごとに行なわれる。つまり、雑音抑圧部５３は、周波数ごとに、既知雑音と未知雑音の何れがより深刻かをみて、深刻な方に合わせて抑圧する。
なお、選択部５３５に対して劣化信号を入力し、選択部５３５が劣化信号を分析して、既知雑音情報と未知雑音情報との何れを抑圧に用いるかを決定してもよい。この場合、抑圧部５３６は、選択部５３５から供給された既知雑音情報または未知雑音情報を、劣化信号から減算することにより雑音抑圧処理を行なう。
本実施形態によれば、より効率的に雑音抑圧を実現することができる。例えば、選択部５３５が利得の最小値を選択すれば、より深刻な雑音を効果的に抑圧できる。つまり、本実施形態によれば、既知雑音が未知雑音より全体的に大きい状況では、既知雑音が大きい周波数では既知雑音を優先的に抑圧し、既知雑音が小さい周波数では未知雑音を優先的に抑圧することができる。一方、選択部５３５が利得の最大値を選択すれば、雑音抑圧による歪みの発生を効果的に防止でき高い音質を実現できる。さらに雑音抑圧部５３は、歪みが発生しない利得の閾値Ｇｔｈが分かっている場合には、その閾値Ｇｔｈよりも利得ｇ１、ｇ２の両方が大きい場合には小さい方の利得を用いて抑圧を行ない、何れかの利得が閾値Ｇｔｈよりも小さい場合には、その閾値Ｇｔｈを用いて抑圧を行なってもよい。
（第６実施形態）
図６は、本発明の第６実施形態としての情報処理装置に含まれる雑音抑圧部６３の構成を示すブロック図である。本実施形態に係る雑音抑圧部６３は、既知雑音情報出力部４３１が劣化信号の分析結果に応じて既知雑音情報を生成する点で、第５実施形態と異なる。そのため、既知雑音情報出力部４３１には劣化信号が供給されている。劣化信号を分析して、その分析結果に応じて既知雑音情報を生成する過程については、図４を用いて第４実施形態で説明した通りであるので、その説明を省略する。また、その他の構成及び動作は第５実施形態と同様であるため、同じ構成については同じ符号を用いてその説明を省略する。
本実施形態のように、入力した劣化信号の分析結果に応じて既知雑音情報を生成すれば、第５実施形態の効果に加えて、変動の激しい信号特性に対応した雑音抑圧を実現することができる。
（第７実施形態）
図７は、本発明の第７実施形態としての情報処理装置に含まれる雑音抑圧部７３の構成を示すブロック図である。本実施形態に係る雑音抑圧部７３は、既知雑音情報出力部２３１から出力された既知雑音情報と、未知雑音情報出力部３３３から出力された未知雑音とを混合して雑音抑圧処理を行なう点で、上記第５実施形態と異なる。他の構成は第５実施形態と同様であるため同じ符号を用いて説明を省略する。
雑音抑圧部７３において混合部７３５は、例えば、既知雑音情報に関する利得ｇ１（０＜ｇ１＜１）と、未知雑音情報に関する利得ｇ２（０＜ｇ２＜１）とをそれぞれ算出し、それらを混合した値（例えば中間値）を抑圧部５３６に供給する。抑圧部５３６は、混合部７３５から供給された利得を劣化信号に積算することにより雑音抑圧処理を行なう。ここで、雑音抑圧部７３の雑音抑圧処理は、劣化信号の周波数成分ごとに行なわれる。つまり、周波数ごとに、既知雑音と未知雑音の抑圧が行なわれる。
なお、混合部７３５に対して劣化信号を入力し、混合部７３５が劣化信号を分析して、既知雑音情報と未知雑音情報とを分析結果に応じて混合して、抑圧に用いてもよい。つまり、混合部７３５が、既知雑音情報と、未知雑音情報とを、入力された劣化信号の分析結果に応じた重み付けで混合してもよい。
本実施形態によれば、抑圧の利得、あるいは抑圧のために減算するべき成分をより正確に求めることができ、効果的に雑音抑圧を実現することができる。その結果、雑音の抑圧と、それによる歪みの発生防止とをバランスよく実現することが可能となる。
（第８実施形態）
図８は、本発明の第８実施形態としての情報処理装置に含まれる雑音抑圧部８３の構成を示すブロック図である。本実施形態に係る雑音抑圧部８３は、既知雑音情報出力部４３１が劣化信号の分析結果に応じて既知雑音情報を生成する点で、第７実施形態と異なる。そのため、既知雑音情報出力部４３１には、劣化信号が供給されている。劣化信号を分析して、その分析結果に応じて既知雑音情報を生成する過程については、図４を用いて第４実施形態で説明した通りであるので、その説明を省略する。また、その他の構成及び動作は第７実施形態と同様であるため、同じ構成については同じ符号を用いてその説明を省略する。
本実施形態のように、入力した劣化信号の分析結果に応じて既知雑音情報を生成すれば、第７実施形態の効果に加えて、変動の激しい信号特性に対応した雑音抑圧を実現することができる。
（第９実施形態）
本実施形態では、図９を用いて、第１乃至第８実施形態で説明した雑音抑圧部３、２３、３３、４３、５３、６３、７３、８３の周辺構成の一例を説明する。図９及び以下の説明において、代表して雑音抑圧部３について示す。
図９に示すように、雑音抑圧部３の周囲には、入力端子１、変換部２、逆変換部４、出力端子５が設けられている。
入力端子１には、劣化信号が、サンプル値系列として供給される。入力端子１に供給された劣化信号は、変換部２においてフーリエ変換などの変換を施されて複数の周波数成分に分割される。複数の周波数成分の振幅スペクトルは雑音抑圧部３へ供給され、位相スペクトルは、逆変換部４に伝達される。なお、ここでは、雑音抑圧部３に振幅スペクトルが供給されているが、本発明はこれに限定されるものではなく、その二乗に相当するパワースペクトルが雑音抑圧部３に供給されても良い。
雑音抑圧部３は、変換部２から供給された劣化信号振幅スペクトルの各周波数で雑音を抑圧し、雑音抑圧結果としての強調信号振幅スペクトルを逆変換部４に伝達する。逆変換部４は、雑音抑圧部３から供給された強調信号振幅スペクトルと変換部２から供給された劣化信号の位相スペクトルとを合わせて逆変換を行い、強調信号サンプルとして、出力端子５に供給する。
［変換部の構成］
図１０は、変換部２の構成を示すブロック図である。図１０に示すように、変換部２はフレーム分割部１２１、窓がけ処理部（ｗｉｎｄｏｗｉｎｇ ｕｎｉｔ）１２２、及びフーリエ変換部１２３を含む。劣化信号サンプルは、フレーム分割部１２１に供給され、Ｋ／２サンプル毎のフレームに分割される。ここで、Ｋは偶数とする。フレームに分割された劣化信号サンプルは、窓がけ処理部１２２に供給され、窓関数（ｗｉｎｄｏｗ ｆｕｎｃｔｉｏｎ）であるｗ（ｔ）との乗算が行なわれる。第ｎフレームの入力信号ｙｎ（ｔ）（ｔ＝０，１，．．．，Ｋ／２−１）に対するｗ（ｔ）で窓がけ（ｗｉｎｄｏｗｉｎｇ）された信号は、次式（１）で与えられる。

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

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

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

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

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

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

なお、図１０と図１１において変換部２と逆変換部４における変換をフーリエ変換として説明したが、変換部２と逆変換部４は、フーリエ変換に代えて、コサイン変換、修正コサイン変換、アダマール変換、ハール変換、ウェーブレット変換など、他の変換を用いることもできる。例えば、コサイン変換や修正コサイン変換は、変換結果として振幅だけしか得ないため、図９における変換部２から逆変換部４に至る経路は不要になる。また、雑音情報記憶部６に記録する雑音情報は振幅（又はパワー）だけとなるため、記憶容量の削減、雑音抑圧処理における演算量の削減に貢献する。ハール変換は乗算を不要とし、ＬＳＩ化したときの面積を小さくすることができる。ウェーブレット変換は、周波数によって時間解像度を異なったものに変更できるため、雑音抑圧効果の向上が期待できる。
また、変換部２が周波数成分を複数統合してから、雑音抑圧部３が実際の抑圧を行うこともできる。その際、変換部２は、聴覚特性の弁別能力が高い低周波領域から、能力が低い高周波領域に向かって、よりたくさんの周波数成分を統合することにより、高い音質を達成することができる。このように、複数の周波数成分を統合してから雑音抑圧が実行されることで、雑音抑圧を適用する周波数成分の数が少なくなり、全体の演算量が削減される。
［雑音抑圧部の処理］
雑音抑圧部３においては、様々な抑圧を行うことが可能である。抑圧方法には、代表的なものとして、ＳＳ（Ｓｐｅｃｔｒｕｍ Ｓｕｂｔｒａｃｔｉｏｎ：スペクトル減算）法とＭＭＳＥ ＳＴＳＡ（Ｍｉｎｉｍｕｍ Ｍｅａｎ−Ｓｑｕａｒｅ Ｅｒｒｏｒ Ｓｈｏｒｔ−Ｔｉｍｅ Ｓｐｅｃｔｒａｌ Ａｍｐｌｉｔｕｄｅ Ｅｓｔｉｍａｔｏｒ：最小二乗平均誤差短時間振幅スペクトル推定）法とがある。ＳＳ法は、雑音情報を、変換部２から供給された劣化信号振幅スペクトルから減算する。ＭＭＳＥ ＳＴＳＡ法は、雑音情報と変換部２から供給された劣化信号振幅スペクトルとを用いて、複数の周波数成分それぞれに対して抑圧係数を計算し、この抑圧係数を劣化信号振幅スペクトルに乗算する。この抑圧係数は、強調信号の平均二乗パワーを最小化するように決定される。
雑音抑圧部３は、雑音の抑圧に際しては、過剰な抑圧を避けるために、フロアリングを適用することができる。フロアリングとは、最大抑圧量を超える抑圧を避ける方法であり、フロアリングパラメータが最大抑圧量を決定する。ＳＳ法は、補正雑音情報を劣化信号振幅スペクトルから減算した結果が、フロアリングパラメータより小さくならないように制約をかける。具体的には、ＳＳ法は、減算結果がフロアリングパラメータよりも小さいときには、減算結果をフロアリングパラメータで置換する。また、ＭＭＳＥ ＳＴＳＡ法は、補正雑音情報と劣化信号振幅スペクトルから求めた抑圧係数が、フロアリングパラメータよりも小さいときに、抑圧係数をフロアリングパラメータで置換する。フロアリングの詳細に関しては、文献「Ｍ．Ｂｅｒｏｕｔｉ，Ｒ．Ｓｃｈｗａｒｔｚ ａｎｄ Ｊ．Ｍａｋｈｏｕｌ，″Ｅｎｈａｎｃｅｍｅｎｔ ｏｆ ｓｐｅｅｃｈ ｃｏｒｒｕｐｔｅｄ ｂｙ ａｃｏｕｓｔｉｃ ｎｏｉｓｅ，″Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＩＣＡＳＳＰ’７９，ｐｐ．２０８−−２１１，Ａｐｒ．１９７９」に開示されている。雑音抑圧部３は、フロアリングパラメータを導入することによって、過剰な抑圧を生じることがなく、強調信号の歪が大きくなることを防止することができる。
雑音抑圧部３は、雑音情報の周波数成分数を劣化信号スペクトルの周波数成分数よりも小さく設定することもできる。このとき、複数の雑音情報が複数の周波数成分に対して共用される。劣化信号スペクトルと雑音情報の双方に対して、複数の周波数成分を統合する場合と比べて、劣化信号スペクトルの周波数分解能が高いので、雑音抑圧部３は、周波数成分の統合が全くない場合よりも少ない演算量で、高い音質を達成することができる。劣化信号スペクトルの周波数成分数よりも少ない周波数成分数の雑音情報を用いた抑圧の詳細は、特開２００８−２０３８７９号に開示されている。
（第１０実施形態）
本発明の第１０実施形態について、図１２を用いて説明する。第９実施形態と比べた場合、本実施形態に係る雑音抑圧部３に、雑音抑圧結果がフィードバックされている点で異なる。他の構成は第９実施形態と同様であるため、ここでは同じ構成については同じ符号を付してその説明を省略する。
本実施形態において、雑音抑圧部３は、雑音抑圧結果としての強調信号振幅スペクトルに基づいた倍率係数を抑圧すべき雑音の雑音情報（例えば利得）に乗算することにより雑音情報を補正する。
雑音抑圧部３は、帰還（ｆｅｅｄｂａｃｋ）された雑音抑圧結果を用いて倍率係数を更新するときには、所望信号が入力されていないタイミングでの雑音抑圧結果が大きいほど（抑圧されずに残った雑音が大きいほど）補正後の雑音情報が大きくなるように、倍率係数を更新する。所望信号が入力されていないタイミングでの雑音抑圧結果が大きいということは、抑圧が不十分であることを示し、倍率係数を変更することによって補正後の雑音情報を大きくすることが望ましいからである。補正後の雑音情報が大きいときには、ＳＳ法では減算する値が大きくなり、雑音抑圧結果は小さくなる。また、ＭＭＳＥ ＳＴＳＡ法のような乗算型の抑圧では、抑圧係数の計算に用いる信号対雑音比の推定値が小さくなり、小さな抑圧係数が得られる。これは、より強力な雑音抑圧をもたらす。倍率係数を更新するにあたって、複数の方法が考えられる。例として、再計算法及び逐次更新法について説明する。
雑音抑圧結果としては、雑音が完全に抑圧された状態が理想である。このため、雑音抑圧部３は、例えば、劣化信号の振幅又はパワーが小さいときに、雑音が完全に抑圧されるように、倍率係数を再計算又は逐次更新することができる。劣化信号の振幅又はパワーが小さいときには、抑圧しようとする雑音以外の信号のパワーも小さい確率が高いからである。雑音抑圧部３は、劣化信号の振幅又はパワーが小さいことを、劣化信号の振幅又はパワーが閾値よりも小さいことを用いて検出できる。
また、雑音抑圧部３は、劣化信号の振幅又はパワーが小さいことを、劣化信号の振幅又はパワーと記憶部３１１に記録されている雑音情報との差分が、閾値より小さいことを用いても検出できる。すなわち、雑音抑圧部３は、劣化信号の振幅又はパワーが雑音情報と似ているときに、劣化信号における雑音情報の占有率が高い（信号対雑音比が低い）ことを利用する。特に、雑音抑圧部３は、複数の周波数点における情報を複合的に用いることにより、スペクトル概形を比較することが可能となり、検出精度を高くすることができる。
ＳＳ法における倍率係数は、各周波数において、補正雑音情報が、所望信号が入力されていないタイミングでの劣化信号スペクトルに等しくなるように、再計算される。言い換えれば、雑音抑圧部３は、雑音だけが入力された時点で変換部２から供給された劣化信号振幅スペクトル｜Ｙｎ（ｋ）｜と、倍率係数αｎと雑音情報ν（ｋ）との積が一致するように倍率係数αｎを求める。ここでｎはフレーム番号、ｋは、周波数番号である。すなわち、倍率係数αｎ（ｋ）は次式（８）で計算される。
αｎ（ｋ）＝｜Ｙｎ（ｋ）｜／νｎ（ｋ） ・・・（８）
一方、ＳＳ法における倍率係数の逐次更新は、各周波数において、所望信号が入力されていないタイミングでの強調信号振幅スペクトルがゼロに近づくように、倍率係数を少しずつ更新する。雑音抑圧部３は、逐次更新に最小二乗平均（Ｌｅａｓｔ Ｓｑｕａｒｅｓ Ｍｅｔｈｏｄ，ＬＭＳ）アルゴリズムを用いる場合には、ｎ番目フレーム、周波数番号ｋの誤差ｅｎ（ｋ）を用いて、αｎ＋１（ｋ）を次式（９）で計算する。
αｎ＋１（ｋ）＝αｎ（ｋ）＋μｅｎ（ｋ）／νｎ（ｋ）・・・（９）
ただし、μはステップサイズと呼ばれる微小定数である。
雑音抑圧部３は、計算して得られた倍率係数αｎ（ｋ）を直ちに利用するときには、数式（９）の代わりに以下の数式（１０）を用いる。
αｎ（ｋ）＝αｎ−１（ｋ）＋μｅｎ（ｋ）／νｎ（ｋ）・・・（１０）
すなわち、雑音抑圧部３は、現在の誤差を用いて現在の倍率係数αｎ（ｋ）を計算し、直ちに適用する。雑音抑圧部３は、倍率係数を直ちに更新することにより、リアルタイムで高精度の雑音抑圧を実現できる。
雑音抑圧部３は、正規化最小二乗平均（Ｎｏｒｍａｌｉｚｅｄ Ｌｅａｓｔ Ｓｑｕａｒｅｓ Ｍｅｔｈｏｄ，ＮＬＭＳ）アルゴリズムを用いる場合には、上述の誤差ｅｎ（ｋ）を用いて、倍率係数αｎ＋１（ｋ）を次式（１１）で計算する。
αｎ＋１（ｋ）＝αｎ（ｋ）＋μｅｎ（ｋ）νｎ（ｋ）／σｎ（ｋ）^２・・・（１１）
σｎ（ｋ）^２は、雑音情報νｎ（ｋ）の平均パワーであり、ＦＩＲフィルタに基づく平均（スライド窓を用いた移動平均）やＩＩＲフィルタに基づく平均（漏れ積分）などを用いて計算できる。
また、雑音抑圧部３は、摂動法を用いて、以下の式（１２）によって倍率係数αｎ＋１（ｋ）を計算しても良い。
αｎ＋１（ｋ）＝αｎ（ｋ）＋μｅｎ（ｋ）・・・（１２）
また、雑音抑圧部３は、誤差の符号だけ表わす符号関数ｓｇｎ｛ｅｎ（ｋ）｝を用いて、以下の式（１３）によって倍率係数αｎ＋１（ｋ）を計算しても良い。
αｎ＋１（ｋ）＝αｎ（ｋ）＋μ・ｓｇｎ｛ｅｎ（ｋ）｝・・・（１３）
同様に、雑音抑圧部３は、最小二乗アルゴリズム（Ｌｅａｓｔ Ｓｑｕａｒｅｓ，ＬＳ）アルゴリズムやその他の適応アルゴリズムを用いてもよい。また、雑音抑圧部３は、更新した倍率係数を直ちに適用することも可能であり、数式（９）から数式（１０）への変更を参照して、数式（１１）〜数式（１３）を変形して、倍率係数をリアルタイム更新してもよい。
ＭＭＳＥ ＳＴＳＡ法は、倍率係数を逐次更新する。雑音抑圧部３は、各周波数において、数式（８）から数式（１３）を用いて説明した方法と同様の方法で、倍率係数αｎ（ｋ）を更新する。
上述した倍率係数の更新方法としての再計算法と逐次更新法について、再計算法は追従速度が速く、逐次更新法は精度が高いという特徴がある。これらの特徴を活かすために、雑音抑圧部３は、最初は再計算法を用い、後に逐次更新法を用いる、というように更新方法を変更することも可能である。雑音抑圧部３は、更新方法の変更のタイミングを決定するにあたり、倍率係数が最適値に十分近くなったことを条件として更新方法を変更することもできる。また、例えば、雑音抑圧部３は、予め定められた時間が経過したときに更新方法を変更してもよい。またさらに、雑音抑圧部３は、倍率係数の補正量が予め定められた閾値よりも小さくなったときに変更することもできる。
以上、本実施形態によれば、雑音抑圧に用いられる雑音情報を雑音抑圧結果に基づいて更新するので、未知な雑音を含む多種多様な雑音を抑圧することができる。
（第１１実施形態）
本発明の第１１実施形態について、図１３を用いて説明する。第１１実施形態としての情報処理装置に含まれる雑音抑圧部３には、入力端子９から、入力した劣化信号中に既知雑音が存在するか否かを示す情報（既知雑音存在情報）が供給される。図５乃至図８の例では、選択部５３５または混合部７３５に既知雑音存在情報を供給し、選択部５３５または混合部７３５は既知雑音存在情報に応じて選択または混合を行なう。
これにより、特定の雑音が存在しているタイミングで、確実に雑音を抑圧することができる。
（第１２実施形態）
本発明の第１２実施形態について、図１４を用いて説明する。第１２実施形態としての情報処理装置に含まれる雑音抑圧部３には、所望信号存在判定部８からの出力が供給される。所望信号存在判定部８には、変換部２からの劣化信号振幅スペクトルが伝達され、そこで、劣化信号振幅スペクトルを解析し、所望信号が存在するか否か、或いは、どの程度存在するのかを判定する。図５乃至図８の例では、選択部５３５または混合部７３５に所望信号存在判定部８からの出力が供給され、選択部５３５または混合部７３５は、供給された判定結果に応じて選択または混合を行なう。
本実施形態によれば、劣化信号中の雑音の割合に応じて雑音の抑圧を行なうので、結果的により精度の高い雑音抑圧結果を得ることができる。
（他の実施形態）
以上説明してきた第１乃至第１２実施形態は、それぞれ別々の特徴を持つ情報処理装置について説明したが、それらの特徴を如何様に組み合わせた情報処理装置も、本発明の範疇に含まれる。
また、本発明は、複数の機器から構成されるシステムに適用されても良いし、単体の装置に適用されても良い。さらに、本発明は、実施形態の機能を実現するソフトウェアの信号処理プログラムが、システム或いは装置に直接或いは遠隔から供給される場合にも適用可能である。したがって、本発明の機能をコンピュータで実現するために、コンピュータにインストールされるプログラム、或いはそのプログラムを格納した媒体、そのプログラムをダウンロードさせるＷＷＷサーバも、本発明の範疇に含まれる。
図１５は、第１実施形態を信号処理プログラムにより構成した場合に、その信号処理プログラムを実行するコンピュータ１５００の構成図である。コンピュータ１５００は、入力部１５０１と、ＣＰＵ１５０２と、雑音情報記憶部１５０３と、出力部１５０４と、メモリ１５０５とを含む。
ＣＰＵ１５０２は、メモリ１５０５に格納された信号処理プログラムを読み込むことにより、コンピュータ１５００全体の動作を制御する。すなわち、信号処理プログラムを実行したＣＰＵ１５０２は、既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する（Ｓ１５２１）。次に、ＣＰＵ１５０２は、未知の性質を有する未知雑音を抑圧するため、未知雑音を推定して未知雑音情報を出力する（Ｓ１５２２）。更に、ＣＰＵ１５０２は、既知雑音情報と未知雑音情報とを用いて入力部１５０１で入力された劣化信号に含まれた雑音を抑圧し（Ｓ１５２３）、出力部１５０４から強調信号を出力する。これにより、第１実施形態と同様の効果を得ることができる。
［実施形態の他の表現］
上記の実施形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。
（付記１）
入力された劣化信号に含まれた雑音を抑圧する雑音抑圧手段を備えた情報処理装置であって、
前記雑音抑圧手段は、
既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する第１出力手段と、
未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力する第２出力手段と、
を含み、前記既知雑音情報と前記未知雑音情報とを用いて雑音抑圧処理を行なうことを特徴とする情報処理装置。
（付記２）
前記雑音抑圧手段は、
前記第１出力手段から出力された既知雑音情報を用いて雑音抑圧処理を行なった信号に対して、前記第２出力手段から出力された未知雑音情報を用いて更に雑音抑圧処理を加えることを特徴とする付記１に記載の情報処理装置。
（付記３）
前記第２出力手段は、入力された前記劣化信号に基づいて、前記未知雑音を推定することを特徴とする付記１または２に記載の情報処理装置。
（付記４）
前記第１出力手段は、入力された前記劣化信号に基づいて、記憶した前記雑音情報から既知雑音情報を生成することを特徴とする付記１、２または３に記載の情報処理装置。
（付記５）
前記雑音抑圧手段は、
前記第１出力手段から出力された前記既知雑音情報及び前記第２出力手段から出力された前記未知雑音情報の何れか一方を選択する選択手段を更に含み、
前記選択手段によって選択された雑音情報を用いて、雑音抑圧処理を行なうことを特徴とする付記１に記載の情報処理装置。
（付記６）
前記選択手段は、
前記劣化信号に基づいて、前記既知雑音情報及び前記未知雑音情報の何れか一方を選択することを特徴とする付記５に記載の情報処理装置。
（付記７）
前記雑音抑圧手段は、
前記第１出力手段から出力された前記既知雑音情報及び前記第２出力手段から出力された前記未知雑音情報を混合する混合手段を更に含み、
前記混合手段による混合によって生成された混合雑音情報を用いて、雑音抑圧処理を行なうことを特徴とする付記１に記載の情報処理装置。
（付記８）
前記混合手段は、
前記劣化信号に基づいて、前記既知雑音情報及び前記未知雑音情報を混合することを特徴とする付記７に記載の情報処理装置。
（付記９）
入力した音声信号をスペクトルに変換して前記劣化信号を生成する変換手段と、
前記雑音抑圧手段によって雑音を抑圧されたスペクトルを逆変換して出力する音声信号を生成する逆変換手段と、
を更に備えたことを特徴とする付記１乃至８の何れかに記載の情報処理装置。
（付記１０）
前記第１出力手段は、前記雑音抑圧手段による抑圧結果に応じて、既知雑音情報を補正して出力することを特徴とする付記１乃至９の何れかに記載の情報処理装置。
（付記１１）
劣化信号中に抑圧すべきでない所望信号がどの程度存在するか否かを判定する判定手段を更に有し、
前記雑音抑圧手段は、前記判定手段による判定結果に応じて雑音抑圧処理を行なうことを特徴とする付記１乃至１０の何れかに記載の情報処理装置。
（付記１２）
入力された劣化信号に含まれた雑音を抑圧する信号処理方法であって、
既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力し、
未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力し、
前記既知雑音情報と前記未知雑音情報とを用いて前記劣化信号に含まれた雑音を抑圧することを特徴とする信号処理方法。
（付記１３）
入力された劣化信号に含まれた雑音を抑圧する信号処理プログラムであって、
既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する第１出力ステップと、
未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力する第２出力ステップと、
前記既知雑音情報と前記未知雑音情報とを用いて前記劣化信号に含まれた雑音を抑圧する雑音抑圧ステップと、
をコンピュータに実行させることを特徴とする信号処理プログラム。
以上、実施形態を参照して本願発明を説明したが、本願発明は上記実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解しうる様々な変更をすることができる。
この出願は、２０１０年５月２５日に出願された日本出願特願２０１０−１１９４９５を基礎とする優先権を主張し、その開示の全てをここに取り込む。

Claims (13)

1. 入力された劣化信号に含まれた既知の雑音と、該既知の雑音とともに存在する未知の雑音を抑圧するための雑音抑圧手段を備え、
   前記雑音抑圧手段は、
   既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する第１出力手段と、
   未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力する第２出力手段と、
   を含み、前記既知雑音情報と前記未知雑音情報とを用いて雑音抑圧処理を行なう情報処理装置。
2. 前記雑音抑圧手段は、
   前記第１出力手段から出力された既知雑音情報を用いて雑音抑圧処理を行なった信号に対して、前記第２出力手段から出力された未知雑音情報を用いて更に雑音抑圧処理を加える請求項１に記載の情報処理装置。
3. 前記第２出力手段は、入力された前記劣化信号に基づいて、前記未知雑音を推定する請求項１または２に記載の情報処理装置。
4. 前記第１出力手段は、入力された前記劣化信号に基づいて、記憶した前記雑音情報から既知雑音情報を生成する請求項１、２または３に記載の情報処理装置。
5. 前記雑音抑圧手段は、
   前記第１出力手段から出力された前記既知雑音情報及び前記第２出力手段から出力された前記未知雑音情報の何れか一方を選択する選択手段を更に含み、
   前記選択手段によって選択された雑音情報を用いて、雑音抑圧処理を行なう請求項１に記載の情報処理装置。
6. 前記選択手段は、
   前記劣化信号に基づいて、前記既知雑音情報及び前記未知雑音情報の何れか一方を選択する請求項５に記載の情報処理装置。
7. 前記雑音抑圧手段は、
   前記第１出力手段から出力された前記既知雑音情報及び前記第２出力手段から出力された前記未知雑音情報を混合する混合手段を更に含み、
   前記混合手段による混合によって生成された混合雑音情報を用いて、雑音抑圧処理を行なう請求項１に記載の情報処理装置。
8. 前記混合手段は、
   前記劣化信号に基づいて、前記既知雑音情報及び前記未知雑音情報を混合する請求項７に記載の情報処理装置。
9. 入力した劣化信号をスペクトルに変換して前記劣化信号を生成する変換手段と、
   前記雑音抑圧手段によって雑音を抑圧されたスペクトルを逆変換して出力する強調信号を生成する逆変換手段と、
   を更に備える請求項１乃至８の何れか１項に記載の情報処理装置。
10. 前記第１出力手段は、前記雑音抑圧手段による抑圧結果に応じて、既知雑音情報を補正して出力する請求項１乃至９の何れか１項に記載の情報処理装置。
11. 劣化信号中に抑圧すべきでない所望信号がどの程度存在するか否かを判定する判定手段を更に有し、
    前記雑音抑圧手段は、前記判定手段による判定結果に応じて雑音抑圧処理を行なう請求項１乃至１０の何れか１項に記載の情報処理装置。
12. 入力された劣化信号に含まれた既知の雑音と、該既知の雑音とともに存在する未知の雑音を抑圧するための信号処理方法であって、
    既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力し、
    未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力し、
    前記既知雑音情報と前記未知雑音情報とを用いて、入力された劣化信号に含まれた雑音を抑圧する信号処理方法。
13. 入力された劣化信号に含まれた既知の雑音と、該既知の雑音とともに存在する未知の雑音を抑圧するための信号処理プログラムであって、
    既知の性質を有する既知雑音を抑圧するため、予め記憶された雑音情報を用いて既知雑音情報を出力する第１出力ステップと、
    未知の性質を有する未知雑音を抑圧するため、前記未知雑音を推定して、未知雑音情報を出力する第２出力ステップと、
    前記既知雑音情報と前記未知雑音情報とを用いて、入力された劣化信号に含まれた雑音を抑圧する雑音抑圧ステップと、
    をコンピュータに実行させる信号処理プログラム。

Images