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

Abstract

本発明は、変換素子を通じて入力した信号を処理するにあたり、変換素子の性能の違いや個体差による出力のバラツキを補正し、精度の高い信号処理を行なう。変換素子を通じて入力信号を入力する入力手段と、基準変換素子を通じて入力した基準信号の最小値を記憶する記憶手段と、前記入力信号の最小値と前記基準信号の最小値とを比較する比較手段と、前記比較手段による比較の結果に応じて、前記入力信号に補正を加える補正手段と、を含む。

Description

本発明は、信号を処理して所望の出力を得るための信号処理技術に関する。

変換素子を用いて入力した信号を処理して所望の出力を得るための信号処理技術が知られている。例えば、劣化信号（所望の信号に雑音が重畳された信号）から、雑音を抑圧し、強調信号（所望の信号を強調した信号）を出力する雑音抑圧技術（ｎｏｉｓｅ ｓｕｐｐｒｅｓｓｉｎｇ ｔｅｃｈｎｏｌｏｇｙ）が存在する。所望の音声信号に重畳されている雑音（ノイズ）を抑圧するノイズサプレッサは、携帯電話など様々な音声端末において利用されている。
この種の技術の一例として、特許文献１は、入力信号に１より小さな抑圧係数を乗算することによって、ノイズを抑圧する方法を開示する。特許文献２は、推定された雑音を劣化信号から直接減算することによって、雑音を抑圧する方法を開示する。また、特許文献３は、雑音が所望信号に対して十分に小さいという条件が満たされない場合にも、十分な雑音抑圧効果と強調信号における小さな歪とを実現できる雑音抑圧システムを開示する。特許文献３は、所望信号に混入する雑音の特性が事前にある程度わかる場合を想定している。特許文献３に記載の技術は、事前に記録しておいた雑音情報（雑音の特性に関する情報）を、劣化信号から減算することで、雑音を抑圧する。

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

しかしながら、上述の特許文献１乃至３に開示された構成では、変換素子の性能の違いや個体差によって、出力にバラツキが出てしまい、精度の高い信号処理を行なうことができなかった。
以上を踏まえ、本発明は、上述の課題を解決する信号処理技術を提供することを目的とする。

上記目的を達成するため、本発明に係る装置は、変換素子を通じて入力信号を入力する入力手段と、基準変換素子を通じて入力した基準信号の最小値を記憶する記憶手段と、前記入力信号の最小値と前記基準信号の最小値とを比較する比較手段と、前記比較手段による比較の結果に応じて、前記入力信号に補正を加える補正手段と、を含む。
上記目的を達成するため、本発明に係る方法は、変換素子を通じて入力信号を入力し、基準変換素子を通じて入力した基準信号の最小値と、入力信号の最小値とを比較し、前記比較の結果に応じて、前記入力信号に補正を加える。
上記目的を達成するため、本発明に係るプログラム記録媒体に格納されるプログラムは、変換素子を通じて入力信号を入力するステップと、基準変換素子を通じて入力した基準信号の最小値と、入力信号の最小値とを比較するステップと、前記比較の結果に応じて、前記入力信号に補正を加えるステップと、をコンピュータに実行させる。

本発明によれば、変換素子の性能の違いや個体差による出力のバラツキを補正し、精度の高い信号処理を行なう信号処理技術を提供することができる。

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

以下に、図面を参照して、本発明の実施の形態について例示的に詳しく説明する。ただし、以下の実施の形態に記載されている構成要素はあくまで例示であり、本発明の技術範囲をそれらのみに限定する趣旨のものではない。なお、以下の説明で、「変換素子」とは、いわゆるトランスデューサーである。具体的には、「変換素子」は、測定、情報転送を含む様々な目的のために、ある種類のエネルギーを別のものに変える電気的、電子的な素子または電気機械である。「変換素子」は、たとえばセンサやマイクロフォン（以下、マイク）のようにある測定量を電気信号に変える素子、機器などを含む。
（第１実施形態）
本発明の第１実施形態としての信号処理装置１００について、図１を用いて説明する。
信号処理装置１００は、入力部１０１と、基準最小値記憶部１０２と、比較部１０３と、補正部１０４とを含む。入力部１０１は、変換素子１１１を通じて入力信号１２０を、比較部１０３及び補正部１０４に入力する。基準最小値記憶部１０２は、基準変換素子を通じて入力した基準信号の最小値（基準最小値）を記憶する。更に、比較部１０３は、入力信号１２０の最小値と基準最小値とを比較する。補正部１０４は、比較部１０３による比較の結果に応じて、入力信号１２０に補正を加える。
以上の構成により、本実施形態に係る信号処理装置１００は、変換素子の性能の違いや個体差による出力のバラツキを補正し、精度の高い信号処理を行なうことができる。
（第２実施形態）
本発明に係る信号処理方法を実現する第２実施形態として、雑音抑圧装置２００について説明する。図２は、雑音抑圧装置２００の全体構成を示すブロック図である。雑音抑圧装置２００は、例えばデジタルカメラ、ノートパソコン、携帯電話などといった装置の一部としても機能するが、本発明はこれに限定されるものではない。雑音抑圧装置２００は、入力信号からのノイズ除去を要求されるあらゆる信号処理装置に適用可能である。
＜全体構成＞
図２に示すように、雑音抑圧装置２００は、入力部２０１と、最小値記憶部２０２と、利得算出部２０３と、補正部２０４と、出力部２０５とを含む。これらのうち、入力部２０１は、変換素子としてのマイク２１１と、マイク２１１の出力に変換処理を施す変換部２１２とを含む。入力部２０１は、音声信号を周波数成分に分解し、比較手段としての利得算出部２０３と補正部２０４とに供給する。
マイク２１１は、劣化信号（所望信号と雑音の混在する信号）が、サンプル値系列として供給される。マイク２１１に劣化信号が供給されると、変換部２１２は、供給された劣化信号にフーリエ変換などの変換を施して、複数の周波数成分に分割する。変換部２１２は、複数の周波数成分のうち振幅スペクトル２２０を利得算出部２０３及び利得制御部２４１に供給する。変換部２１２は、複数の周波数成分のうち位相スペクトル２３０を、逆変換部２５２に伝達する。利得制御部２４１は、変換部２１２から振幅スペクトルを受ける。利得制御部２４１は、振幅スペクトルに利得を乗算して、その結果を雑音抑圧部２４２へ供給する。なお、ここでは、変換部２１２は、利得制御部２４１を介して雑音抑圧部２４２に振幅スペクトル２２０を供給しているが、本発明はこれに限定されるものではない。変換部２１２は、利得制御部２４１を介して振幅スペクトル２２０の二乗に相当するパワースペクトルを雑音抑圧部２４２に供給しても良い。
最小値記憶部２０２は、半導体メモリなどの記憶素子を含む。最小値記憶部２０２は、雑音に関する基準最小値を記憶している。基準最小値は、基準変換素子の一例としての基準となるマイクで本装置が抑圧しようとする雑音のみを静かな部屋で録音することにより決定されても良い。例えば、本実施形態に係る雑音抑圧装置２００がデジタルカメラに実装された場合を考える。この場合、基準最小値は、雑音抑圧装置２００が実装されたデジタルカメラの電源が入れられた状態で発生しているノイズを基準マイクがひろった値でも良い。雑音抑圧装置２００は、入力部２０１から周波数成分ごとの音声信号が入力される。そのため、本実施形態においては、基準最小値も周波数成分ごとに用意するものとする。しかしながら本発明の実施形態はこれに限定されるものではない。
利得算出部２０３は、変換部２１２から出力された音声信号の各周波数成分の最小値を抽出する最小値抽出部２３１を含む。また、利得算出部２０３は、抽出した最小値を、最小値記憶部２０２から読出した基準最小値と比較する最小値比較部２３２を含む。利得算出部２０３は、抽出した最小値と基準最小値との比を用いて（例えば、抽出した最小値が基準最小値と一致するように）入力信号に適用すべき周波数成分ごとの利得制御値（補正係数）を算出する。最小値抽出部２３１は、変換部２１２から供給された劣化信号振幅（またはパワースペクトル）を１サンプルごと、または数サンプルごとに解析して、最小値を導き出す。最小値抽出部２３１は、解析する度に最小値を更新し、過去の全入力中の最小値を抽出する。つまり抽出が長時間になればなるほど、最小値は小さくなる。具体的には、最小値抽出部２３１は、例えば１サンプル目の最小値と２サンプル目の最小値を比較し、更に３サンプル目の最小値を比較して更新する。このため長くサンプリングすればするほどどんどん最小値は小さくなる。
また、最小値抽出部２３１は、最小値を一定時間毎にリセットしてもよい。リセットの間隔が長くなるほど、最小値は劣化信号中の最小の成分を表すようになる。劣化信号が所望信号と雑音から構成されて、雑音の方が所望信号よりも信号レベルが低いときは、劣化信号の最小値は雑音の最小値になる。最小値記憶部２０２は、基準最小値として、雑音だけを静かな環境で録音して得られた最小値を記憶している。従って、利得算出部２０３は、同じ雑音の最小値を比較することになり、利得制御の基本データを得ることができる。
利得制御部２４１は、利得算出部２０３で算出した利得に基づいて、利得の制御を行なう。利得制御のタイミングは、１サンプルごとでもよいし、所定サンプル数ごとでもよい。なお、雑音抑圧装置２００は、全周波数に同じ利得を用いて調整しても良い。つまり変換部２１２でフーリエ変換を行なう前に、変換部２１２が最小値によって利得調整しても良い。
雑音情報記憶部２０７は、半導体メモリなどの記憶素子を含む。雑音情報記憶部２０７は、雑音情報（雑音の特性に関する情報）を記憶している。雑音情報は、例えば、雑音のスペクトルの形でも良い。雑音情報は、スペクトルの形の他、位相の周波数特性、特定の周波数における強弱や時間変化などの特徴量でも良い。雑音情報は、その他、統計量（最大、最小、分散、メジアン）などでも良い。スペクトルが１０２４の周波数成分で表わされている場合、雑音情報記憶部２０７は、１０２４の振幅（またはパワー）データを記憶する。雑音情報記憶部２０７は、１０２４の振幅（またはパワー）データに代えて、複数の周波数成分を統合して得られたサブバンドのデータを記憶していても良い。サブバンドを用いると、雑音抑圧装置２００は、必要なメモリ量と演算量を削減することができる。そして最小値記憶部２０２は、それぞれのスペクトルについての最小値を記憶している。雑音情報記憶部２０７に記録された雑音情報は、雑音情報調整部２４３に供給される。雑音情報調整部２４３は、倍率係数を乗算することにより雑音情報を補正し、補正雑音情報として雑音抑圧部２４２に供給する。
雑音抑圧部２４２は、利得制御部２４１から供給された劣化信号振幅スペクトルと雑音情報調整部２４３から供給された補正雑音情報２６０とを用いて、各周波数で雑音を抑圧する。雑音抑圧部２４２は、雑音抑圧結果としての強調信号振幅スペクトル２４０を逆変換部２５２に伝達する。
雑音抑圧部２４２は、同時に、強調信号振幅スペクトル２４０を雑音情報調整部２４３に伝達する。雑音情報調整部２４３は、雑音抑圧結果としての強調信号振幅スペクトル２４０に基づいて、雑音情報を補正する。
逆変換部２５２は、雑音抑圧部２４２から供給された強調信号振幅スペクトル２４０と変換部２１２から供給された劣化信号の位相スペクトル２３０とを合わせて逆変換を行い、強調信号サンプルとして、出力端子２５１に供給する。
＜変換部２１２の構成＞
図３は、変換部２１２の内部構成を示すブロック図である。図３に示すように、変換部２１２はフレーム分割部３０１、窓がけ処理部（ｗｉｎｄｏｗｉｎｇ ｕｎｉｔ）３０２、及びフーリエ変換部３０３を含む。劣化信号サンプルは、フレーム分割部３０１に供給され、Ｋ／２サンプル毎のフレームに分割される。ここで、Ｋは偶数とする。フレームに分割された劣化信号サンプルは、窓がけ処理部３０２に供給され、窓関数（ｗｉｎｄｏｗ ｆｕｎｃｔｉｏｎ）であるｗ（ｔ）との乗算が行なわれる。第ｎフレームの入力信号ｙｎ（ｔ）（ｔ＝０，１，．．．，Ｋ／２−１）に対するｗ（ｔ）で窓がけ（ｗｉｎｄｏｗｉｎｇ）された信号は、次式（１）で与えられる。

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

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

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

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

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

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

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

Claims (12)

1. 変換素子を通じて入力信号を入力する入力手段と、
   基準変換素子を通じて入力した基準信号の最小値を記憶する記憶手段と、
   前記入力信号の最小値と前記基準信号の最小値とを比較する比較手段と、
   前記比較手段による比較の結果に応じて、前記入力信号に補正を加える補正手段と、
   を含む信号処理装置。
2. 前記比較手段は、前記入力信号の最小値と前記基準信号の最小値との比を求め、
   前記補正手段は、前記比較手段が求めた比に応じて、前記入力信号を利得制御する、
   請求項１に記載の信号処理装置。
3. 前記補正手段は、前記入力信号の最小値と前記基準信号の最小値とが一致するように補正係数を定め、前記補正係数を用いて、前記変換素子の出力を補正する、
   請求項１または２に記載の信号処理装置。
4. 前記補正手段は、
   雑音情報を用いて劣化信号中の雑音を抑圧する雑音抑圧手段と、
   前記比較手段による比較の結果に応じて、前記雑音情報を調整して前記雑音抑圧手段に供給する雑音情報調整手段と、
   を含む請求項１、２または３に記載の信号処理装置。
5. 前記雑音情報調整手段は、更に、前記劣化信号中の雑音を抑圧した結果に基づいて前記雑音情報を調整する、
   請求項４に記載の信号処理装置。
6. 前記補正手段は、前記雑音情報調整手段に供給すべき前記雑音情報を記憶する雑音情報記憶手段を更に含む、
   請求項４または５に記載の信号処理装置。
7. 前記補正手段は、雑音源から前記雑音情報を入力して、前記雑音抑圧に用いる、
   請求項４または５に記載の信号処理装置。
8. 前記補正手段は、前記入力信号中に雑音が存在するか否かを示す情報を入力し、前記入力信号中に雑音が存在している場合に、前記入力信号の補正を行なう、
   請求項４乃至７の何れか１項に記載の信号処理装置。
9. 前記補正手段は、前記入力信号中に所望信号がどの程度存在しているかを判定し、その判定結果に基づいて、前記雑音情報を調整する、
   請求項４乃至８の何れか１項に記載の信号処理装置。
10. 前記変換素子はマイクロフォンである、
    請求項１乃至９の何れか１項に記載の信号処理装置。
11. 変換素子を通じて入力信号を入力し、
    基準変換素子を通じて入力した基準信号の最小値と、入力信号の最小値とを比較し、
    前記比較の結果に応じて、前記入力信号に補正を加える、
    信号処理方法。
12. 変換素子を通じて入力信号を入力するステップと、
    基準変換素子を通じて入力した基準信号の最小値と、入力信号の最小値とを比較するステップと、
    前記比較の結果に応じて、前記入力信号に補正を加えるステップと、
    をコンピュータに実行させる信号処理プログラムを格納したプログラム記録媒体。

Images