(57)【特許請求の範囲】
【請求項１】
入力信号を複数の帯域に分割し、帯域信号を出力する周波数分割手段と、
前記帯域信号の振幅成分を求める振幅算出手段と、
前記入力信号に含まれる雑音の振幅成分を推定して推定雑音振幅成分を帯域毎に求める雑音推定手段と、
帯域毎に異なる重み係数を発生する重み係数発生手段と、
前記複数の帯域信号それぞれに対応する前記重み係数を乗算して前記複数の帯域信号の振幅成分を時間的に平滑化した平滑化振幅成分を求める振幅平滑化手段と、
帯域毎に前記平滑化振幅成分と前記推定雑音振幅成分から抑圧係数を求める抑圧量算出手段と、
帯域毎に前記帯域信号を前記抑圧係数に基づいて抑圧する雑音抑圧手段と、
前記雑音抑圧手段が出力する複数の帯域の雑音抑圧後の帯域信号を合成して出力する周波数合成手段
を有することを特徴とする雑音抑圧装置。
【請求項２】
入力信号を複数の帯域に分割し、帯域信号を出力する周波数分割手段と、
前記帯域信号の振幅成分を求める振幅算出手段と、
前記入力信号に含まれる雑音の振幅成分を推定して推定雑音振幅成分を帯域毎に求める雑音推定手段と、
帯域毎に異なる重み係数を時間的に変化させ、出力する重み係数発生手段と、
前記複数の帯域信号それぞれに対応する前記重み係数を乗算して前記複数の帯域信号の振幅成分を時間的に平滑化した平滑化振幅成分を求める振幅平滑化手段と、
帯域毎に前記平滑化振幅成分と前記推定雑音振幅成分から抑圧係数を求める抑圧量算出手段と、
帯域毎に前記帯域信号を前記抑圧係数に基づいて抑圧する雑音抑圧手段と、
前記雑音抑圧手段が出力する複数の帯域の雑音抑圧後の帯域信号を合成して出力する周波数合成手段
を有することを特徴とする雑音抑圧装置。
【請求項３】
請求項１または２記載の雑音抑圧装置において、
前記重み係数発生手段は、予め設定された重み係数を出力することを特徴とする雑音抑圧装置。
【請求項４】
請求項１または２記載の雑音抑圧装置において、
前記重み係数発生手段は、前記入力信号の振幅成分に基づいて帯域毎に重み係数を算出することを特徴とする雑音抑圧装置。
【請求項５】
請求項１または２記載の雑音抑圧装置において、
前記重み係数発生手段は、前記平滑化振幅成分に基づいて帯域毎に重み係数を算出することを特徴とする雑音抑圧装置。
【請求項６】
請求項１または２記載の雑音抑圧装置において、
前記重み係数発生手段は、前記入力信号の振幅成分と前記推定雑音振幅成分の比に基づいて帯域毎に重み係数を算出することを特徴とする雑音抑圧装置。
【請求項７】
請求項１または２記載の雑音抑圧装置において、
前記重み係数発生手段は、前記平滑化振幅成分と前記推定雑音振幅成分の比に基づいて帯域毎に重み係数を算出することを特徴とする雑音抑圧装置。
【請求項８】
請求項１乃至７のいずれか記載の雑音抑圧装置において、
前記重み係数発生手段は、時間的な総和が１となる重み係数を発生することを特徴とする雑音抑圧装置。
【請求項９】
請求項１乃至８のいずれか記載の雑音抑圧装置において、
前記周波数分割手段は、高速フーリエ変換器であり、
前記周波数合成手段は、高速逆フーリエ変換器であることを特徴とする雑音抑圧装置。
【請求項１０】
請求項１乃至８のいずれか記載の雑音抑圧装置において、
前記周波数分割手段は、複数のバンドパスフィルタで構成され、
前記周波数合成手段は、加算回路で構成されたことを特徴とする雑音抑圧装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、雑音抑圧装置に関し、雑音が重畳した音声信号から雑音成分を低減させる雑音抑圧装置に関する。
【背景技術】
【０００２】
携帯電話システムやＩＰ（Ｉｎｔｅｒｎｅｔ  Ｐｒｏｔｏｃｏｌ）電話システム等において、マイクには話者の音声に加えて環境騒音が入力される。その結果、音声信号が劣化し、音声の明瞭感が損なわれる。そこで、劣化した音声信号から雑音成分を低減させ、通話品質を高める技術が従来から開発されている（例えば、非特許文献１及び特許文献１参照）。【０００３】
図１は、従来の雑音抑圧装置の一例のブロック図を示す。同図中、時間周波数変換部１０は、単位時間（フレーム）毎に、現フレームｎの入力信号ｘｎ（ｋ）を時間領域ｋから周波数領域ｆに変換し、入力信号の周波数領域信号Ｘｎ（ｆ）を求める。振幅算出部１１は周波数領域信号Ｘｎ（ｆ）から入力信号の振幅成分｜Ｘｎ（ｆ）｜（以下、「入力振幅成分」という）を求める。雑音推定部１２は話者の音声が無い場合の入力振幅成分｜Ｘｎ（ｆ）｜から推定雑音の振幅成分μｎ（ｆ）（以下、「推定雑音振幅成分」という）を求める。【０００４】
抑圧係数算出部１３は｜Ｘｎ（ｆ）｜とμｎ（ｆ）から（１）式にしたがって抑圧係数Ｇｎ（ｆ）を求める。【０００５】
【数１】
0004423300.tif 000002
雑音抑圧部１４はＸｎ（ｆ）とＧｎ（ｆ）から（２）式にしたがって雑音抑圧後の振幅成分Ｓ＊ｎ（ｆ）を求める。【０００６】
【数２】
0004423300.tif 000003
時間周波数変換部１５はＳ＊ｎ（ｆ）を周波数領域から時間領域に変換し、雑音抑圧後の信号ｓ＊ｎ（ｋ）を求める。【０００７】
（非特許文献１）Ｓ．Ｆ．Ｂｏｌｌ，”Ｓｕｐｒｅｓｓｉｏｎ  ｏｆ  Ａｃｏｕｓｔｉｃ  Ｎｏｉｓｅ  ｉｎ  Ｓｐｅｅｃｈ  Ｕｓｉｎｇ  Ｓｐｅｃｔｒａｌ  Ｓｕｂｔｒａｃｔｉｏｎ”，  ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎ  ｏｎ  Ａｃｏｕｓｔｉｃｓ，Ｓｐｅｅｃｈ，ａｎｄ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ，  ＡＳＳＰ−３３，  ｖｏｌ．２７，ｐｐ．１１３−１２０，１９７９
（特許文献１）特開２００４−２０６７９
【発明の開示】
【発明が解決しようとする課題】
【０００８】
図１では、推定雑音振幅成分μｎ（ｆ）は、例えば、過去の話者の音声が含まれないフレームにおける入力信号の振幅成分を平均することで求める。このように背景雑音の平均的な（長期的な）傾向は過去の入力振幅成分に基づいて推定する。【０００９】
図２は、従来の抑圧係数算出方法の一例の原理図を示す。同図中、抑圧係数算出部１６では、現フレームｎの振幅成分｜Ｘｎ（ｆ）｜と推定雑音振幅成分μｎ（ｆ）に基づいて抑圧係数Ｇｎ（ｆ）を算出し、この抑圧係数を入力振幅成分に乗算することで、入力信号に含まれる雑音成分を抑圧する。【００１０】
しかしながら、現フレームに重畳している（短期的な）雑音の振幅成分を正確に求めることは困難である。すなわち、現フレームに重畳している雑音の振幅成分と推定雑音振幅成分との間には推定誤差（以下、雑音推定誤差）が生じる。このため、図３に示すように、実線で示す雑音の振幅成分と、破線で示す推定雑音振幅成分との差である雑音推定誤差が大きくなる。【００１１】
この結果、雑音抑圧装置において上記の雑音推定誤差は過剰抑圧や抑圧不足を引き起こす。さらに、雑音推定誤差がフレーム毎に大きく変動するので過剰抑圧や抑圧不足も変動し、雑音抑圧性能に時間的なムラが発生してしまう。この雑音抑圧性能の時間的なムラがミュージカルノイズ（ｍｕｓｉｃａｌ  ｎｏｉｓｅ）として知られる異音を生じさせる。【００１２】
図４は、従来の抑圧係数算出方法の他の一例の原理図を示す。このものは、雑音抑圧装置における過剰抑圧や抑圧不足に伴って発生する異音を抑止することを目的とした平均化雑音抑圧技術である。同図中、振幅平滑化部１７では現フレームｎの振幅成分｜Ｘｎ（ｆ）｜の平滑化を行い、抑圧係数算出部１８は、平滑化後の入力信号の振幅成分Ｐｎ（ｆ）（以下、「平滑化振幅成分」という）と、推定雑音振幅成分μｎ（ｆ）に基づいて抑圧係数Ｇｎ（ｆ）を求める。【００１３】
振幅成分の平滑化方法としては、次の２つの方法が用いられる。
（第１の平滑化方法）
現フレームと過去数フレーム分の入力振幅成分の平均値を平滑化振幅成分Ｐｎ（ｆ）とする。この方法は単純な平均化であり、平滑化振幅成分は（３）式により求めることができる。【００１４】
【数３】
0004423300.tif 000004
（第２の平滑化方法）
現フレームの振幅成分｜Ｘｎ（ｆ）｜と、直前フレームの平滑化振幅成分Ｐｎ−１（ｆ）との荷重平均値を平滑化振幅成分Ｐｎ（ｆ）とする。これは指数平滑化と呼ばれ、平滑化振幅成分は（４）式により求めることができる。【００１５】
【数４】
0004423300.tif 000005
図４の抑圧係数算出方法では、抑圧係数を算出する前に入力振幅成分を平均化または指数平滑化することにより、話者の音声が入力されないときは、図５に示すように、実線で示す雑音の振幅成分と、破線で示す推定雑音振幅成分との差である雑音推定誤差を小さくできる。この結果、図２の抑圧係数算出で問題であった、雑音入力時の過剰抑圧や抑圧不足を抑えることができ、ミュージカルノイズを抑止することができる。【００１６】
しかしながら、話者の音声が入力されるときには、図６に示すように平滑化振幅成分が鈍ってしまい、破線で示す音声信号の振幅成分と、実線で示す平滑化振幅成分との誤差（以下、「音声推定誤差」という）が大きくなる。【００１７】
この結果、音声推定誤差の大きな平滑化振幅成分と推定雑音振幅に基づいて抑圧係数を求め、入力振幅成分に抑圧係数を乗算するので、入力信号に含まれる音声成分を誤って抑圧してしまい、音質の劣化を引き起こすという問題があった。この現象は特に音声の話頭（音声の始まりの区間）において顕著である。【００１８】
本発明は、上記の点に鑑みなされたものであり、ミュージカルノイズの発生を抑えつつ、音声への影響を最小限にし、安定した雑音抑圧性能を実現する雑音抑圧装置を提供することを総括的な目的とする。
【課題を解決するための手段】
【００１９】
この目的を達成するため、本発明は、入力信号の振幅成分を帯域毎に求める振幅算出手段と、前記入力信号の雑音の振幅成分を推定して推定雑音振幅成分を帯域毎に求める雑音推定手段と、帯域毎に異なる重み係数を発生する重み係数発生手段と、前記複数の帯域信号それぞれに対応する前記帯域毎に異なる重み係数を乗算して前記複数の帯域信号の振幅成分を時間的に平滑化した平滑化振幅成分を求める振幅平滑化手段と、帯域毎に前記平滑化振幅成分と前記推定雑音振幅成分から抑圧係数を求める抑圧量算出手段と、帯域毎に前記帯域信号を前記抑圧係数に基づいて抑圧する雑音抑圧手段と、前記雑音抑圧手段が出力する複数の帯域の雑音抑圧後の帯域信号を合成して出力する周波数合成手段を有し構成する。
【発明の効果】
【００２０】
このような雑音抑圧装置によれば、ミュージカルノイズの発生を抑えつつ、音声への影響を最小限にし、安定した雑音抑圧性能を実現できる。
【図面の簡単な説明】
【００２１】
【図１】従来の雑音抑圧装置の一例のブロック図である。
【図２】従来の抑圧係数算出方法の一例の原理図である。
【図３】従来の雑音推定誤差を説明するための図である。
【図４】従来の抑圧係数算出方法の他の一例の原理図
【図５】従来の雑音推定誤差を説明するための図である。
【図６】従来の音声推定誤差を説明するための図である。
【図７】本発明における抑圧係数算出の原理図である。
【図８】本発明における抑圧係数算出の原理図である。
【図９】ＦＩＲフィルタを使用する場合の振幅平滑化部の構成図である。
【図１０】ＩＩＲフィルタを使用する場合の振幅平滑化部の構成図である。
【図１１】本発明における重み係数の一例を示す図である。
【図１２】平滑化振幅成分と推定雑音振幅成分から抑圧係数を求める関係式を示す図である。
【図１３】本発明の雑音推定誤差を説明するための図である。
【図１４】本発明の音声推定誤差を説明するための図である。
【図１５】雑音が重畳した音声の入力信号の波形図である。
【図１６】従来の雑音抑圧装置の出力音声信号の波形図である。
【図１７】本発明の雑音抑圧装置の出力音声信号の波形図である。
【図１８】本発明の雑音抑圧装置の第１実施形態のブロック図である。
【図１９】本発明の雑音抑圧装置の第２実施形態のブロック図である。
【図２０】本発明の雑音抑圧装置の第３実施形態のブロック図である。
【図２１】非線形関数ｆｕｎｃを示す図である。
【図２２】本発明の雑音抑圧装置の第４実施形態のブロック図である。
【図２３】信号対雑音比と重み係数の関係を示す図である。
【図２４】本発明の雑音抑圧装置の第５実施形態のブロック図である。
【図２５】本発明装置を適用した携帯電話の一実施形態のブロック図である。
【図２６】本発明装置を適用した携帯電話の他の実施形態のブロック図である。
【符号の説明】
【００２２】
２１  振幅平滑化部
２２  抑圧係数算出部
２３  重み係数算出部
３０  ＦＦＴ部
３１，４１  振幅算出部
３２，４２  雑音推定部
３３  振幅平滑化部
３４  振幅保持部
３５  重み係数保持部
３６，４６  抑圧係数算出部
３７，４７  雑音抑圧部
４０  チャンネル分割部
４３  振幅平滑化部
４４  振幅保持部
４５  重み係数算出部
４８  チャンネル合成部
【発明を実施するための最良の形態】
【００２３】
以下、図面に基づいて本発明の実施形態について説明する。【００２４】
本発明における抑圧係数算出の原理図を図７および図８に示す。本発明では図４と同様に、抑圧係数を算出する前に入力振幅成分の平滑化を行う。【００２５】
図７において、振幅平滑化部２１では現フレームｎの振幅成分｜Ｘｎ（ｆ）｜と重み係数ｗｍ（ｆ）を用いて平滑化振幅成分Ｐｎ（ｆ）を得る。抑圧係数算出部２２は、平滑化振幅成分Ｐｎ（ｆ）と、推定雑音振幅成分μｎ（ｆ）に基づいて抑圧係数Ｇｎ（ｆ）を求める。【００２６】
図８において、重み係数算出部２３は、入力振幅成分から特徴量（信号対雑音比や入力信号の振幅等）を算出し、特徴量に基づいて重み係数ｗｍ（ｆ）を適応的に制御する。振幅平滑化部２１では現フレームｎの振幅成分｜Ｘｎ（ｆ）｜と重み係数算出部２３からの重み係数ｗｍ（ｆ）を用いて平滑化振幅成分Ｐｎ（ｆ）を得る。抑圧係数算出部２２は、平滑化振幅成分Ｐｎ（ｆ）と、推定雑音振幅成分μｎ（ｆ）に基づいて抑圧係数Ｇｎ（ｆ）を求める。【００２７】
平滑化の方法としてはＦＩＲフィルタを使用する方法とＩＩＲフィルタを使用する方法があり、本発明ではいずれの平滑化方法を選択してもよい。
（ＦＩＲフィルタを使用する場合）
ＦＩＲフィルタを使用する場合の振幅平滑化部２１の構成を図９に示す。同図中、振幅保持部２５では、過去Ｎフレーム分の入力振幅成分（平滑化前の振幅成分）を保持する。さらに平滑化部２６において、過去Ｎフレーム分の平滑化前の振幅成分と、現在の振幅成分から（５）式にしたがって平滑化後の振幅成分を求める。【００２８】
【数５】
0004423300.tif 000006
（ＩＩＲフィルタを使用する場合）
ＩＩＲフィルタを使用する場合の振幅平滑化部の構成を図１０に示す。同図中、振幅保持部２７では、過去Ｎフレーム分の平滑化後の振幅成分を保持する。さらに平滑化部２８において、過去Ｎフレーム分の平滑化後の振幅成分と、現在の振幅成分から（６）式にしたがって平滑化後の振幅成分を求める。【００２９】
【数６】
0004423300.tif 000007
上記の（５），（６）式において、ｍはフィルタを構成する遅延素子数であり、ｗ０（ｆ）〜ｗｍ（ｆ）はフィルタを構成するｍ＋１個の乗算器それぞれの重み係数であり、この値を調整することで、入力信号を平滑化する際の平滑化の強さを制御することができる。【００３０】
従来は（３），（４）式から明らかなように全周波数帯域において同じで重み係数を使用していたが、本発明では（５），（６）式のように、重み係数ｗｍ（ｆ）は、周波数の関数として表わされ、帯域毎に異なる値を使用することを特徴としている。【００３１】
図１１に本発明における重み係数ｗ０（ｆ）の一例を示す。図１１では、入力信号の性質が低周波数帯域は変動しにくく高周波数帯域は変動しやすいことを想定しており、現フレームの振幅成分｜Ｘｎ（ｆ）｜にかかる重み係数ｗ０（ｆ）を実線に示すように低域で大きく高域で小さくなる値とすることで、高周波数帯域の変動に追随し、かつ、低周波数帯域には平滑化をより強くかけている。なお、帯域毎に重み係数の時間的な総和は１とされており、ｗ１（ｆ）＝１−ｗ０（ｆ）とした場合、ｗ１（ｆ）は一点鎖線に示すようになる。【００３２】
また、従来の（４）式は重み係数としての平滑化係数αが定数であったが、本発明では、重み係数ｗｍ（ｆ）を変数として、図８に示す重み係数算出部２３で入力振幅成分から信号対雑音比や入力信号の振幅などの特徴量を算出し、特徴量に基づいて重み係数を適応的に制御する。【００３３】
平滑化振幅成分Ｐｎ（ｆ）と推定雑音振幅成分μｎ（ｆ）から抑圧係数Ｇｎ（ｆ）を求める際の関係式としては、任意のものを選択可能である。例えば（１）式を用いても良く、さらには、図１２に示すような関係式を適用しても良い。図１２では、Ｐｎ（ｆ）／μｎ（ｆ）が小さくなるほどＧｎ（ｆ）を小さくしている。【００３４】
本発明の雑音抑圧装置では、抑圧係数を算出する前に入力振幅成分を平滑化するため、話者の音声が入力されないときは、図１３に示すように、実線で示す雑音の振幅成分と、破線で示す推定雑音振幅成分との差である雑音推定誤差を小さくできる。【００３５】
さらに、話者の音声が入力されるときにも、図１４に示すように、破線で示す音声信号の振幅成分と、実線で示す平滑化振幅成分との差である音声推定誤差を小さくできる。この結果、音声への影響を最小限にしつつ、ミュージカルノイズの発生を抑え安定した雑音抑圧性能を実現できる。【００３６】
ここで、図１５に示すように、雑音が重畳した音声の入力信号が供給された場合、図４の抑圧係数算出方法を用いた従来の雑音抑圧装置の出力音声信号は図１６に示す波形となり、本発明の雑音抑圧装置の出力音声信号は図１７に示す波形となる。【００３７】
図１６の波形と図１７の波形を比べてみると、話頭の区間τにおいて、図１７の波形の劣化が小さいことが分かる。それぞれの出力音を比較するために、雑音入力時の抑圧性能を音声が無い区間で測定し、音声入力時の音質劣化を話頭の区間で測定した結果を以下に示す。【００３８】
雑音入力時の抑圧性能（非音声区間で測定）は、従来の雑音抑圧装置が約１４ｄＢで、本発明の雑音抑圧装置が約１４ｄＢである。音声入力時の音質劣化（音声の話頭区間で測定）は、従来の雑音抑圧装置が約４ｄＢであるのに対し、本発明の雑音抑圧装置は約１ｄＢとなり、約３ｄＢ改善される。これにより、本発明は音声入力時に音声成分の抑圧を小さくして音質劣化を軽減することが可能となる。【００３９】
図１８は、本発明の雑音抑圧装置の第１実施形態のブロック図を示す。この実施形態はチャネル分割・合成にＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）／ＩＦＦＴ（Ｉｎｖｅｒｓｅ  ＦＦＴ）を使用し、ＦＩＲフィルタによる平滑化方法を採用し、抑圧係数の算出には（１）式を採用している。【００４０】
同図中、ＦＦＴ部３０は、単位時間（フレーム）毎に、現フレームｎの入力信号ｘｎ（ｋ）を時間領域ｋから周波数領域ｆに変換し、入力信号の周波数領域信号Ｘｎ（ｆ）を求める。なお、添字ｎはフレーム番号を表している。【００４１】
振幅算出部３１は周波数領域信号Ｘｎ（ｆ）から入力振幅成分｜Ｘｎ（ｆ）｜を求める。雑音推定部３２は音声区間検出を行い、話者音声の非検出時に入力振幅成分｜Ｘｎ（ｆ）｜から（７）式にしたがって推定雑音振幅成分μｎ（ｆ）を求める。【００４２】
【数７】
0004423300.tif 000008
振幅平滑化部３３は、入力振幅成分｜Ｘｎ（ｆ）｜と、振幅保持部３４で保持している直前フレームの入力振幅成分｜Ｘｎ−１（ｆ）｜と、重み係数保持部３５で保持している重み係数ｗｍ（ｆ）から（８）式にしたがって平均化振幅成分Ｐｎ（ｆ）を求める。ただし、ｆＳは音声をデジタル化する際のサンプリング周波数である。また、重み係数ｗｍ（ｆ）は図１１に示すものとする。【００４３】
【数８】
0004423300.tif 000009
抑圧係数算出部３６は、平均化振幅成分Ｐｎ（ｆ）と推定雑音振幅成分μｎ（ｆ）から（９）式にしたがって抑圧係数Ｇｎ（ｆ）を求める。
【００４４】
【数９】
0004423300.tif 000010
雑音抑圧部３７は、Ｘｎ（ｆ）とＧｎ（ｆ）から（１０）式にしたがって雑音抑圧後の振幅成分Ｓ＊ｎ（ｆ）を求める。【００４５】
【数１０】
0004423300.tif 000011
ＩＦＦＴ部３７は、振幅成分Ｓ＊ｎ（ｆ）を周波数領域から時間領域に変換し、雑音抑圧後の信号ｓ＊ｎ（ｋ）を求める。【００４６】
図１９は、本発明の雑音抑圧装置の第２実施形態のブロック図を示す。この実施形態はチャネル分割・合成にバンドパスフィルタを使用し、ＦＩＲフィルタによる平滑化方法を採用し、抑圧係数の算出には（１）式を採用している。【００４７】
同図中、チャンネル分割部４０は、帯域フィルタ（ＢＰＦ）を用いて入力信号ｘｎ（ｋ）を（１１）式にしたがって帯域信号ｘＢＰＦ（ｉ，ｋ）に分割する。なお、添字ｉはチャンネル番号を表す。【００４８】
【数１１】
0004423300.tif 000012
振幅算出部４１は、帯域信号ｘＢＰＦ（ｉ，ｋ）から各フレームにおいて（１２）式にしたがって帯域別入力振幅Ｐｏｗ（ｉ，ｎ）を算出する。なお、添字ｎはフレーム番号を表す。【００４９】
【数１２】
0004423300.tif 000013
雑音推定部４２は音声区間検出を行い話者音声の非検出時に帯域別入力振幅成分Ｐｏｗ（ｉ，ｎ）から（１３）式にしたがって推定雑音の振幅成分μ（ｉ，ｎ）を求める。【００５０】
【数１３】
0004423300.tif 000014
重み係数算出部４５は、帯域別入力振幅成分Ｐｏｗ（ｉ，ｎ）を所定の閾値ＴＨＲ１と比較して重み係数ｗ（ｉ，ｍ）を算出する。ただし、ｍ＝０，１，２とする。
Ｐｏｗ（ｉ，ｎ）≧ＴＨＲ１のとき、
ｗ（ｉ，０）＝０．７
ｗ（ｉ，１）＝０．２
ｗ（ｉ，２）＝０．１
Ｐｏｗ（ｉ，ｎ）＜ＴＨＲ１のとき、
ｗ（ｉ，０）＝０．４
ｗ（ｉ，１）＝０．３
ｗ（ｉ，２）＝０．３
つまり、チャンネル毎に重み係数の時間的な総和は１とされている。【００５１】
振幅平滑化部４３は、振幅保持部４４で保持している帯域別入力振幅成分Ｐｏｗ（ｉ，ｎ−１），Ｐｏｗ（ｉ，ｎ−２）と、振幅算出部４１からの帯域別入力振幅成分Ｐｏｗ（ｉ，ｎ）と、重み係数ｗ（ｉ，ｍ）から（１４）式にしたがって平滑化入力振幅成分ＰｏｗＡＶ（ｉ，ｎ）を算出する。【００５２】
【数１４】
0004423300.tif 000015
抑圧係数算出部４６は、平滑化入力振幅成分ＰｏｗＡＶ（ｉ，ｎ）と推定雑音振幅成分μ（ｉ，ｎ）から（１５）式により抑圧係数Ｇ（ｉ，ｎ）を算出する。【００５３】
【数１５】
0004423300.tif 000016
雑音抑圧部４７において、帯域信号ｘＢＰＦ（ｉ，ｋ）と抑圧係数Ｇ（ｉ，ｎ）から（１６）式にしたがって雑音抑圧後の帯域信号ｓ＊ＢＰＦ（ｉ，ｋ）を求める。【００５４】
【数１６】
0004423300.tif 000017
チャンネル合成部４８は加算回路で構成され、帯域信号ｓ＊ＢＰＦ（ｉ，ｋ）を（１７）式にしたがって加算合成して出力音声信号ｓ＊（ｋ）を求める。【００５５】
【数１７】
0004423300.tif 000018
図２０は、本発明の雑音抑圧装置の第３実施形態のブロック図を示す。この実施形態はチャネル分割・合成にＦＦＴ／ＩＦＦＴを使用し、ＩＩＲフィルタによる平滑化方法を採用し、抑圧係数の算出には非線形関数を採用している。【００５６】
同図中、ＦＦＴ部３０は、単位時間（フレーム）毎に、現フレームｎの入力信号ｘｎ（ｋ）を時間領域ｋから周波数領域ｆに変換し、入力信号の周波数領域信号Ｘｎ（ｆ）を求める。なお、添字ｎはフレーム番号を表している。【００５７】
振幅算出部３１は周波数領域信号Ｘｎ（ｆ）から入力振幅成分｜Ｘｎ（ｆ）｜を求める。雑音推定部３２は音声区間検出を行い、話者音声の非検出時に入力振幅成分｜Ｘｎ（ｆ）｜から（７）式にしたがって推定雑音振幅成分μｎ（ｆ）を求める。【００５８】
振幅平滑化部５１は、入力振幅成分｜Ｘｎ（ｆ）｜と、振幅保持部５２で保持している過去２フレームの平均化振幅成分Ｐｎ−１（ｆ），Ｐｎ−２（ｆ）と、重み係数保持部５３で保持している重み係数ｗｍ（ｆ）から（１８）式にしたがって平均化振幅成分Ｐｎ（ｆ）を求める。【００５９】
【数１８】
0004423300.tif 000019
重み係数算出部５３は、平均化振幅成分Ｐｎ（ｆ）を所定の閾値ＴＨＲ２と比較して重み係数ｗｍ（ｆ）を算出する。ただし、ｍ＝０，１，２とする。
Ｐｎ（ｆ）≧ＴＨＲ２のとき、
ｗ０（ｆ）＝１．０
ｗ１（ｆ）＝０．０
ｗ２（ｆ）＝０．０
Ｐｎ（ｆ）＜ＴＨＲ２のとき、
ｗ０（ｆ）＝０．６
ｗ１（ｆ）＝０．２
ｗ２（ｆ）＝０．２
つまり、帯域毎に重み係数の時間的な総和は１とされている。
【００６０】
抑圧係数算出部５４は、平均化振幅成分Ｐｎ（ｆ）と推定雑音振幅成分μｎ（ｆ）から（１９）式に示す非線形関数ｆｕｎｃを用いて抑圧係数Ｇｎ（ｆ）を求める。なお、非線形関数ｆｕｎｃを図２１に示す。【００６１】
【数１９】
0004423300.tif 000020
雑音抑圧部３７は、Ｘｎ（ｆ）とＧｎ（ｆ）から（１０）式にしたがって雑音抑圧後の振幅成分Ｓ＊ｎ（ｆ）を求める。ＩＦＦＴ部３７は、振幅成分Ｓ＊ｎ（ｆ）を周波数領域から時間領域に変換し、雑音抑圧後の信号ｓ＊ｎ（ｋ）を求める。【００６２】
このように、重み係数を平滑化後の振幅成分に基づいて制御することにより、非定常な雑音に強固で安定した制御が可能となる。【００６３】
図２２は、本発明の雑音抑圧装置の第４実施形態のブロック図を示す。この実施形態はチャネル分割・合成にＦＦＴ／ＩＦＦＴを使用し、ＦＩＲフィルタによる平滑化方法を採用し、抑圧係数の算出には非線形関数を採用している。【００６４】
同図中、ＦＦＴ部３０は、単位時間（フレーム）毎に、現フレームｎの入力信号ｘｎ（ｋ）を時間領域ｋから周波数領域ｆに変換し、入力信号の周波数領域信号Ｘｎ（ｆ）を求める。なお、添字ｎはフレーム番号を表している。【００６５】
振幅算出部３１は周波数領域信号Ｘｎ（ｆ）から入力振幅成分｜Ｘｎ（ｆ）｜を求める。雑音推定部３２は音声区間検出を行い、話者音声の非検出時に入力振幅成分｜Ｘｎ（ｆ）｜から（７）式にしたがって推定雑音振幅成分μｎ（ｆ）を求める。【００６６】
信号対雑音比計算部５６は、現フレームの入力振幅成分｜Ｘｎ（ｆ）｜と、推定雑音振幅成分μｎ（ｆ）から、（２０）式にしたがって帯域毎に信号対雑音比ＳＮＲｎ（ｆ）を求める。【００６７】
【数２０】
0004423300.tif 000021
重み係数算出部５７は、信号対雑音比ＳＮＲｎ（ｆ）から重み係数ｗ０（ｆ）を求める。なお、ＳＮＲｎ（ｆ）とｗ０（ｆ）の関係を図２３に示す。また、ｗ０（ｆ）からｗ１（ｆ）を（２１）式にしたがって算出する。つまり、帯域毎に重み係数の時間的な総和は１とされている。【００６８】
【数２１】
0004423300.tif 000022
振幅平滑化部５８は、現フレームの入力振幅成分｜Ｘｎ（ｆ）｜と、振幅保持部３４で保持している直前フレームの入力振幅成分｜Ｘｎ−１（ｆ）｜と、重み係数算出部５７からの重み係数ｗｍ（ｆ）即ちｗ０（ｆ），ｗ１（ｆ）から（２２）式にしたがって平均化振幅成分Ｐｎ（ｆ）を求める。【００６９】
【数２２】
0004423300.tif 000023
抑圧係数算出部３６は、平均化振幅成分Ｐｎ（ｆ）と推定雑音振幅成分μｎ（ｆ）から（９）式にしたがって抑圧係数Ｇｎ（ｆ）を求める。雑音抑圧部３７は、Ｘｎ（ｆ）とＧｎ（ｆ）から（１０）式にしたがって雑音抑圧後の振幅成分Ｓ＊ｎ（ｆ）を求める。ＩＦＦＴ部３７は、振幅成分Ｓ＊ｎ（ｆ）を周波数領域から時間領域に変換し、雑音抑圧後の信号ｓ＊ｎ（ｋ）を求める。【００７０】
このように、重み係数を信号対雑音比に基づいて制御することにより、マイクの音量によらず安定した制御が可能となる。【００７１】
図２４は、本発明の雑音抑圧装置の第５実施形態のブロック図を示す。この実施形態はチャネル分割・合成にＦＦＴ／ＩＦＦＴを使用し、ＩＩＲフィルタによる平滑化方法を採用し、抑圧係数の算出には非線形関数を採用している。【００７２】
同図中、ＦＦＴ部３０は、単位時間（フレーム）毎に、現フレームｎの入力信号ｘｎ（ｋ）を時間領域ｋから周波数領域ｆに変換し、入力信号の周波数領域信号Ｘｎ（ｆ）を求める。なお、添字ｎはフレーム番号を表している。【００７３】
振幅算出部３１は周波数領域信号Ｘｎ（ｆ）から入力振幅成分｜Ｘｎ（ｆ）｜を求める。雑音推定部３２は音声区間検出を行い、話者音声の非検出時に入力振幅成分｜Ｘｎ（ｆ）｜から（７）式にしたがって推定雑音振幅成分μｎ（ｆ）を求める。【００７４】
振幅平滑化部５１は、入力振幅成分｜Ｘｎ（ｆ）｜と、振幅保持部５２で保持している過去２フレームの平均化振幅成分Ｐｎ−１（ｆ），Ｐｎ−２（ｆ）と、重み係数算出部６１からの重み係数ｗｍ（ｆ）から（１８）式にしたがって平均化振幅成分Ｐｎ（ｆ）を求める。
【００７５】
信号対雑音比計算部６０において、平滑化振幅成分Ｐｎ（ｆ）と、推定雑音振幅成分μｎ（ｆ）から、（２３）式にしたがって帯域毎に信号対雑音比ＳＮＲｎ（ｆ）を算出する。【００７６】
【数２３】
0004423300.tif 000024
重み係数算出部６１は、信号対雑音比ＳＮＲｎ（ｆ）から重み係数ｗ０（ｆ）を求める。なお、ＳＮＲｎ（ｆ）とｗ０（ｆ）の関係を図２３に示す。また、ｗ０（ｆ）からｗ１（ｆ）を（２１）式にしたがって算出する。【００７７】
抑圧係数算出部５４は、平均化振幅成分Ｐｎ（ｆ）と推定雑音振幅成分μｎ（ｆ）から（１９）式に示す非線形関数ｆｕｎｃを用いて抑圧係数Ｇｎ（ｆ）を求める。雑音抑圧部３７は、Ｘｎ（ｆ）とＧｎ（ｆ）から（１０）式にしたがって雑音抑圧後の振幅成分Ｓ＊ｎ（ｆ）を求める。ＩＦＦＴ部３７は、振幅成分Ｓ＊ｎ（ｆ）を周波数領域から時間領域に変換し、雑音抑圧後の信号ｓ＊ｎ（ｋ）を求める。【００７８】
このように、重み係数を平滑化後の信号対雑音比に基づいて制御することにより、非定常な雑音に強固で安定した制御が可能となり、マイクの音量によらず安定した制御が可能る。【００７９】
図２５は、本発明装置を適用した携帯電話の一実施形態のブロック図を示す。同図中、マイクロホン７１の出力音声信号は、本発明の雑音抑圧装置７０にて雑音抑圧されたのちエンコーダ７２で符号化され、送信部７３から公衆網７４に送出される。【００８０】
図２６は、本発明装置を適用した携帯電話の他の実施形態のブロック図を示す。同図中、公衆網７４から送信された信号は受信部７５で受信され、デコーダ７６にて復号され、本発明の雑音抑圧装置７０にて雑音抑圧される。こののち、スピーカ７７に供給され発音される。【００８１】
なお、図２５と図２６を複合化して送話系と受話系の双方に本発明の雑音抑圧装置７０を設けても良い。【００８２】
なお、振幅算出部３１，４１が請求項記載の振幅算出手段に対応し、雑音推定部３２，４２が雑音推定手段に対応し、重み係数保持部３５，重み係数算出部４５，信号対雑音比計算部５６，６０が重み係数発生手段に対応し、振幅平滑化部３３，４３が振幅平滑化手段に対応し、抑圧係数算出部３６，４６が抑圧量算出手段に対応し、雑音抑圧部３７，４７が雑音抑圧手段に対応し、ＦＦＴ部３０，チャンネル分割部４０が周波数分割手段に対応し、ＩＦＦＴ部３８，チャンネル合成部４８が周波数合成手段に対応する。
【図１】
0004423300.tif 000025
【図２】
0004423300.tif 000026
【図３】
0004423300.tif 000027
【図４】
0004423300.tif 000028
【図５】
0004423300.tif 000029
【図６】
0004423300.tif 000030
【図７】
0004423300.tif 000031
【図８】
0004423300.tif 000032
【図９】
0004423300.tif 000033
【図１０】
0004423300.tif 000034
【図１１】
0004423300.tif 000035
【図１２】
0004423300.tif 000036
【図１３】
0004423300.tif 000037
【図１４】
0004423300.tif 000038
【図１５】
0004423300.tif 000039
【図１６】
0004423300.tif 000040
【図１７】
0004423300.tif 000041
【図１８】
0004423300.tif 000042
【図１９】
0004423300.tif 000043
【図２０】
0004423300.tif 000044
【図２１】
0004423300.tif 000045
【図２２】
0004423300.tif 000046
【図２３】
0004423300.tif 000047
【図２４】
0004423300.tif 000048
【図２５】
0004423300.tif 000049
【図２６】
0004423300.tif 000050