JP2013148724A - 雑音抑圧装置、雑音抑圧方法およびプログラム - Google Patents

雑音抑圧装置、雑音抑圧方法およびプログラム Download PDF

Info

Publication number
JP2013148724A
JP2013148724A JP2012009240A JP2012009240A JP2013148724A JP 2013148724 A JP2013148724 A JP 2013148724A JP 2012009240 A JP2012009240 A JP 2012009240A JP 2012009240 A JP2012009240 A JP 2012009240A JP 2013148724 A JP2013148724 A JP 2013148724A
Authority
JP
Japan
Prior art keywords
noise
band
unit
noise suppression
band power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2012009240A
Other languages
English (en)
Inventor
Kenichi Makino
堅一 牧野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sony Corp
Original Assignee
Sony Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sony Corp filed Critical Sony Corp
Priority to JP2012009240A priority Critical patent/JP2013148724A/ja
Priority to US13/719,696 priority patent/US20130191118A1/en
Priority to CN201310009827.4A priority patent/CN103220440A/zh
Publication of JP2013148724A publication Critical patent/JP2013148724A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L21/0216—Noise filtering characterised by the method used for estimating noise
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L21/0216—Noise filtering characterised by the method used for estimating noise
    • G10L21/0232—Processing in the frequency domain

Landscapes

  • Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Quality & Reliability (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Noise Elimination (AREA)

Abstract

【課題】入力信号から雑音信号の推定を行って雑音信号を選択的に低減する際の音質の向上を図る。
【解決手段】雑音帯域パワー推定部により、帯域パワー算出部で得られた各帯域分割信号の帯域パワーと雑音判定部の判定結果から、各帯域の雑音の帯域パワーを推定する。雑音抑圧ゲイン決定部により、各帯域分割信号の帯域パワーと推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する。雑音帯域パワー推定部は、非定常雑音における雑音変化への追従速度を定常雑音における雑音変化への追従速度より上げる。非定常雑音への追従性能を改善でき、定常雑音のみならず、非定常雑音にも効果的な雑音抑圧を実現できる。
【選択図】図４

Description

本技術は、雑音抑圧装置、雑音抑圧方法およびプログラムに関し、特に、入力信号から雑音信号の推定を行って雑音信号を選択的に低減した出力信号を得る雑音抑圧装置等に関する。
近年、ＶｏＩＰ（Voice over Internet Protocol）や携帯電話などの通信機器やＩＣレコーダなど、マイクロホンで収音した人間の音声をＡＤ（Analog to Digital）変換し、デジタル信号として伝送・記録したうえで再生を行う電子機器が広く普及している。これらの電子機器の使用時には周囲環境から発せられる音がマイクに混入し音声を聞き取る妨げとなる。
そこで、従来、携帯電話などにおいては、入力信号から雑音信号の推定を行い、雑音信号を選択的に低減する雑音抑圧技術が導入されてきた。この種の雑音抑圧技術は、例えば、非特許文献１などに開示されている。
Yariv Ephraim and David Malarah, "SpeechEnhancementUsing a Minimum Mean Square Error Short-Time Spectral Amplitude Estimator",IEEE Transactions on Accoustics, Speech, And Signal Processing,Vol. ASSP-32, No. 6, December1994 pp1109-1121.
雑音には、パワー変動を伴わない定常雑音の他に、衣ずれ音や紙のこすれ音などの摩擦雑音や風切り音といった、雑音性のスペクトル形状を持ちながらもパワー変動を伴う非定常雑音が存在する。
本技術の目的は、定常雑音のみならず、非定常雑音にも効果的な雑音抑圧を実現することにある。
本技術の概念は、
入力信号を所定フレーム長のフレームに分割してフレーム化するフレーム化部と、
上記フレーム化部で得られたフレーム化信号を複数の帯域に分割して帯域分割信号を得る帯域分割部と、
上記帯域分割部で得られた各帯域分割信号から帯域パワーを得る帯域パワー算出部と、 上記フレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定部と、
上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音判定部の判定結果から、各帯域の雑音の帯域パワーを推定する雑音帯域パワー推定部と、
上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する雑音抑圧ゲイン決定部と、
上記帯域分割部で得られた各帯域分割信号に、上記雑音抑圧ゲイン決定部で決定された各帯域の雑音抑圧ゲインを適用して雑音が抑圧された帯域分割信号を得る雑音抑圧部と、
上記雑音抑圧部で得られた各帯域分割信号を帯域合成して雑音抑圧されたフレーム化信号を得る帯域合成部と、
上記帯域合成部で得られた各フレームのフレーム化信号をフレーム合成して雑音を抑圧した出力信号を得るフレーム合成部とを備え、
上記雑音帯域パワー推定部は、
上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
雑音抑圧装置にある。
この発明において、フレーム化部により、入力信号は所定時間長のフレームに分割されてフレーム化される。そして、帯域分割部により、フレーム化信号が複数の帯域に分割されて帯域分割信号が得られる。例えば、この帯域分割部では、フレーム化信号に対して高速フーリエ変換が施されて周波数領域の信号とされて、複数帯域に分割される。
帯域パワー算出部により、帯域分割部で得られた各帯域分割信号から帯域パワーが得られる。この場合、例えば、フーリエ変換で得られる複素スペクトルからパワースペクトルが算出され、パワースペクトルの帯域内の最大値あるいは平均値などが代表値、つまり帯域パワーとされる。
雑音判定部により、フレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かが判定される。つまり、この雑音判定部により、各帯域が、定常雑音であるか、非定常雑音であるか、あるいは音声であるかが判定される。例えば、各帯域が順次判定帯域とされ、この判定帯域の帯域分割信号の現在のフレームと過去のフレームとの帯域パワーが比較され、該帯域パワーの変動が閾値内にあるとき、この判定帯域は定常雑音と判定される。この判定は、雑音のパワーはフレーム間で一定であり、逆にパワー変動が大きい信号は雑音ではないとの仮定に基づいている。また、例えば、各帯域が順次判定帯域とされ、フレーム化信号が非定常雑音の特性を有し、判定帯域内に音声由来のピークが存在しないとき、この判定帯域は非定常雑音と判定される。
雑音帯域パワー推定部により、帯域パワー算出部で得られた各帯域分割信号の帯域パワーと雑音判定部の判定結果から、各帯域の雑音帯域パワーが推定される。この場合、非定常雑音における雑音変化への追従速度は、定常雑音における雑音変化への追従速度より上げられる。例えば、雑音帯域パワー推定部は、帯域毎に、帯域パワー算出部で得られた現在フレームの帯域パワーと１フレーム前に推定された雑音の帯域パワーとを加重加算することにより、現在フレームの推定された雑音のパワーを得、非定常雑音における現在フレームの帯域パワーの重みは、定常雑音における上記現在フレームの帯域パワーの重みより大きくされる。
雑音抑圧ゲイン決定部により、帯域パワー算出部で得られた各帯域分割信号の帯域パワーと雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインが決定される。そして、雑音抑圧部により、帯域分割部で得られた各帯域分割信号に、雑音抑圧ゲイン決定部で決定された各帯域の雑音抑圧ゲインが適用されて、雑音が抑圧された帯域分割信号が得られる。そして、帯域合成部により、雑音抑圧部で得られた各帯域分割信号が帯域合成されて雑音抑圧されたフレーム化信号が得られ、フレーム合成部により、帯域合成部で得られた各フレームのフレーム化信号がフレーム合成されて、雑音が抑圧された出力信号が得られる。
このように、本技術においては、雑音帯域パワー推定部で各帯域の雑音帯域パワーを推定する際に、非定常雑音における雑音変化への追従速度は、定常雑音における雑音変化への追従速度より上げられる。非定常雑音は定常雑音と比べて信号の変化の速度が速いが、非定常雑音の場合はノイズの追従速度が速められるので、こと非定常雑音への追従性能が改善される。従って、定常雑音のみならず、非定常雑音にも効果的な雑音抑圧を実現できる。
なお、本技術において、例えば、雑音抑圧ゲイン決定部は、帯域毎に、帯域パワー算出部で得られた各帯域分割信号の帯域パワーと雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーとからＳＮＲを算出するＳＮＲ算出部と、帯域毎に、ＳＮＲ算出部で算出されたＳＮＲを平滑化するＳＮＲ平滑化部とを有する構成とされてもよい。
この場合、雑音抑圧ゲイン決定部では、ＳＮＲ平滑部で平滑化された各帯域のＳＮＲに基づいて、各帯域の雑音抑圧ゲインが決定される。また、この場合、平滑化の係数は、雑音判定部の判定結果および周波数帯域に基づいて変更される。例えば、雑音抑圧ゲイン決定部では、ＳＮＲ平滑化部で平滑化された各帯域のＳＮＲと共に、ＳＮＲ算出部で算出されたＳＮＲに基づいて、各帯域の雑音抑圧ゲインが決定される、ようにしてもよい。
また、例えば、雑音抑圧ゲイン決定部では、帯域毎に、現在のフレームの信号の帯域パワーと推定された雑音の帯域パワーの比が第１のＳＮＲとされ、直前フレームの信号の帯域のパワーと雑音抑圧ゲインが乗算された量と直前フレームの推定された雑音の帯域パワーの比が第２のＳＮＲとされる。そして、雑音抑圧ゲイン決定部では、第１のＳＮＲと第２のＳＮＲとが用いられて雑音抑圧ゲインが決定される。
このように雑音抑圧ゲイン決定部において、例えば、帯域毎に、平滑化ＳＮＲに基づいて雑音抑圧ゲインが決定されるが、その平滑化係数が雑音判定部の判定結果および帯域に基づいて変更される。例えば、各フレーム、各帯域で、非雑音と判定される場合の平滑化係数（α）は小さい値へと変更され、雑音と判定される場合の平滑化係数（α）は大きい値へと変更される。これにより、信号の時間変化が激しい区間での平滑化ＳＮＲの追従性を改善でき、また、信号の時間変化が少ない区間での平滑化ＳＮＲの無用な変化を抑制できる。そのため、各帯域の雑音抑圧ゲインの精度を高めることができ、音質の劣化を少なく抑えることができる。
また、本技術において、例えば、雑音抑圧ゲイン決定部で決定された雑音抑圧ゲインが予め設定された下限値より小さくなるとき、雑音抑圧ゲインの値をこの下限値に修正する雑音抑圧ゲイン修正部をさらに備え、雑音抑圧部は、雑音抑圧ゲイン修正部で修正された雑音抑圧ゲインを用いる、ようにされてもよい。
この場合、下限値は帯域別に設定される。例えば、非雑音の信号が音声である場合に、音声信号が含まれる確率の高い帯域については雑音抑圧ゲインの下限値が高めの値に設定される。そして、雑音抑圧ゲイン決定部で決定された雑音抑圧ゲインが下限値を下回った場合は下限値で置き換えられる。これにより、雑音抑圧ゲイン決定部で決定された雑音抑圧ゲインの誤差があっても聴感上の音質劣化が少なくされる。
また、本技術の他の概念は、
複数チャネルのそれぞれ所定フレーム長のフレームに分割してフレーム化する複数のフレーム化部と、
上記複数のフレーム化部で得られたフレーム化信号をそれぞれ複数の帯域に分割して帯域分割信号を得る複数の帯域分割部と、
上記複数の帯域分割部で得られた各帯域分割信号からそれぞれ帯域パワーを得る複数の帯域パワー算出部と、
上記複数チャネルのフレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定部と、
上記複数の帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音判定部の判定結果からそれぞれ各帯域の雑音の帯域パワーを推定する複数の雑音帯域パワー推定部と、
上記複数の帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記複数の雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーに基づいてそれぞれ各帯域の雑音抑圧ゲインを決定する複数の雑音抑圧ゲイン決定部と、
上記複数の帯域分割部で得られた各帯域分割信号に、上記複数の雑音抑圧ゲイン決定部で決定された各帯域の雑音抑圧ゲインを適用してそれぞれ雑音が抑圧された帯域分割信号を得る複数の雑音抑圧部と、
上記複数の雑音抑圧部で得られた各帯域分割信号を帯域合成してそれぞれ雑音抑圧されたフレーム化信号を得る複数の帯域合成部と、
上記複数の帯域合成部で得られた各フレームのフレーム化信号をフレーム合成してそれぞれ雑音を抑圧した出力信号を得るフレーム合成部とを備え、
上記雑音帯域パワー推定部は、
上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
雑音抑圧装置にある。
本技術においては、各チャネルにおいて、各帯域の雑音抑圧ゲインが決定されて、雑音抑圧処理が行われる。複数チャネルのフレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かが判定される。例えば、各帯域が順次判定帯域とし、判定帯域について各チャネルで定常雑音および非定常雑音のいずれであるか否かが判定され、全てのチャネルで定常雑音と判定するときこの判定帯域は定常雑音であると判定され、全てのチャネルで非定常雑音と判定するときこの判定帯域は非定常雑音であると判定される。各チャネルにおいて、フレーム毎に各帯域の雑音抑圧ゲインを決定する際には、雑音判定部の判定結果が共通に用いられる。
このように本技術においては、複数チャネル（例えば、ステレオ信号の場合の左右チャネル）における雑音の帯域パワーの推定誤差により複数チャネルの雑音抑圧ゲインに意図しない振幅差が発生することが抑制され、複数チャネルの不整合による定位の破壊が回避される。
本技術によれば、定常雑音のみならず、非定常雑音にも効果的な雑音抑圧を実現できる。
本技術の基本的な雑音低減の方針を示す図である。 雑音のみが存在するフレームにおける雑音低減の効果を説明するための図である。 雑音と音声が重複しているフレームにおける雑音低減の効果を説明するための図である。 本技術の第１の実施の形態としての雑音抑圧装置の構成例を示すブロック図である。 有声音検出部のゼロクロス幅計算部における計算動作を説明するための図である。 フレーム化信号が音声（非雑音）である場合における、信号波形（各サンプルの振幅）と、ゼロクロス幅のヒストグラムの一例を示す図である。 フレーム化信号が音声（雑音）である場合における、信号波形（各サンプルの振幅）と、ゼロクロス幅のヒストグラムの一例を示す図である。 有声帯域判定部が実行する判定処理の一例を示すフローチャートである。 非定常雑音判定部が実行するノイズテンプレートＢN(rmin,b)を求める処理の一例を示すフローチャートである。 非定常雑音判定部が実行する非定常雑音フラグＦnsn(u)の出力処理の一例を示すフローチャートである。 雑音・非雑音判定部の判定処理の手順を説明するためのフローチャートである。 α算出部で算出される重み係数α(k,b)の推移例を示す図である。 本技術の第２の実施の形態としての雑音抑圧装置の構成例を示すブロック図である。 雑音抑圧装置を構成する雑音抑圧ゲイン生成部の構成例を示すブロック図である。 雑音・非雑音判定部の判定処理の手順を説明するためのフローチャートである。 雑音抑圧処理をソフトウェアで行うコンピュータ装置の構成例を示す図である。
以下、発明を実施するための形態（以下、「実施の形態」とする）について説明する。なお、説明を以下の順序で行う。
１．第１の実施の形態
２．第２の実施の形態
３．変形例
図１は、本技術の基本的な雑音低減の方針を示している。雑音のみが含まれるフレームについては、全帯域一律に振幅を下げることによって、雑音低減の効果を得る。一方、音声と雑音が混合されたフレームについては、音声に由来するスペクトルの山は保存し、谷を下げる（掘る）ことによって、雑音低減の効果を得る。
また、本技術では、定常雑音の抑圧を行うスペクトラルサブトラクションの枠組に、非定常雑音の帯域パワーを推定する手段を追加する。非定常雑音は定常雑音と比べて信号の変化の速度が速いため、定常雑音と同様の方法で推定値の更新を行うのでは雑音の変化に追従できなり。そこで、当該フレームの雑音が定常雑音であるか非定常雑音であるかを判定し、非定常雑音の場合はノイズの追従速度を速めることで雑音への追従性能を改善する。
非定常雑音の帯域パワーの推定は、定常雑音の場合と同様、帯域毎に各フレームで信号の状態を監視して雑音と非雑音の判定を行い、雑音と判定されたフレームで雑音の推定値を逐次更新することによって行われる。
雑音のみが存在するフレームについては、図２に示すように、雑音から全ての帯域において雑音推定値を差し引くことで、雑音低減の効果を得るようにされる。ただし、非定常雑音の場合には、定常雑音の場合と同様の追従速度では雑音の振幅変化に追従できずに雑音推定誤差が大きくなり、それが原因で出力の残留雑音が大きくなる結果につながる。そのため、雑音推定の追従速度が上げられる。
一方、雑音と音声が重複しているフレームでは、図３に示すように、非定常スペクトル上で雑音と音声を分離することが困難であるため、スペクトルのピークを音声信号由来であると仮定し、スペクトルのピーク以外の部分、すなわち谷の部分を抑圧することによって雑音抑圧効果を狙う。これを実現するために、スペクトルのピークを検出したうえで、ピーク以外の谷の部分で雑音推定値の更新を行うようにする。この場合も、非定常雑音の場合には、雑音推定の追従速度が上げられる。
ここで、スペクトルのピーク検出を行う際、単純にピークを検出するだけでは偽のピークを検出する恐れがある。そのため、ピークの周波数軸上での間隔が一定に揃っているかチェックするなど、より確実に音声由来のピークを捉えることによって、雑音の推定精度を向上させることができる。
＜１．第１の実施の形態＞
［雑音抑圧装置の構成］
図４は、本技術の第１の実施の形態としての雑音抑圧装置１０の構成例を示している。この雑音抑圧装置１０は、信号入力端子１１と、フレーム化部１２と、窓がけ部１３と、高速フーリエ変換部１４と、雑音抑圧ゲイン生成部１５を有している。また、この雑音抑圧装置１０は、フーリエ係数修正部１６と、逆高速フーリエ変換部１７と、窓がけ部１８と、オーバーラップ加算部１９と、信号出力端子２０を有している。
信号入力端子１１は、入力信号ｙ(n)を供給する端子である。この入力信号ｙ(n)は、標本化周波数がｆsのデジタル信号である。フレーム化部１２は、信号入力端子１１に供給される入力信号ｙ(n)を、フレーム毎の処理を行うために、所定フレーム長、例えばフレーム長がＮfサンプルのフレームに分割してフレーム化する。例えば、ｕ番目のフレームの信号のｎ番目のサンプルは、ｙf(u,n)と表記される。フレーム化部１２におけるフレーム化処理においては、隣接フレームをオーバーラップさせてもよい。
窓がけ部１３は、フレーム化信号ｙf(u,n)に分析窓wana(n)による窓がけを行う。窓がけ部１３は、分析窓wana(n)として、例えば、以下の（１）式で定義されるものを用いる。Ｎｗは、窓長である。
Figure 2013148724
高速フーリエ変換部１４は、窓がけ部１３で窓がけされたフレーム化信号ｙf(u,n)に対して、高速フーリエ変換（ＦＦＴ：Fast Fouriertransform）処理を施し、時間領域信号を周波数領域信号に変換する。雑音抑圧ゲイン生成部１５は、フレーム化処理で得られたフレーム化信号ｙf(u,n)と、高速フーリエ変換処理で得られた各フーリエ係数（各周波数スペクトル）とに基づいて、各フーリエ係数に対応した雑音抑圧ゲインを生成する。この各フーリエ係数に対応した雑音抑圧ゲインは、周波数軸上のフィルタを構成する。この雑音抑圧ゲイン生成部１５の詳細については、後述する。
フーリエ係数修正部１６は、高速フーリエ変換処理で得られた各フーリエ係数と、雑音抑圧ゲイン生成部１５で生成された各フーリエ係数に対応した雑音抑圧ゲインの積をとることで係数修正を行う。つまり、フーリエ係数修正部１６は、周波数軸上で、雑音を抑圧するためのフィルタ計算を行う。
逆高速フーリエ変換部１７は、係数修正された各フーリエ係数に対して、逆高速フーリエ変換（ＩＦＦＴ：Inverse Fast Fourier transform）処理を施す。この逆高速フーリエ変換部１７は、上述の高速フーリエ変換部１４とは逆の処理を行い、周波数領域信号を時間領域信号に変換する。
窓がけ部１８は、逆高速フーリエ変換部１７で得られた雑音抑圧されたフレーム化信号に合成窓wsyn(n)による窓がけを行う。窓がけ部１８は、合成窓wsyn(n)として、例えば、以下の（２）式で定義されるものを用いる。
Figure 2013148724
なお、窓がけ部１３における分析窓wana(n)および窓がけ部１８における合成窓wsyn(n)の形状は任意でよい。しかし、一連の分析・合成系において、完全再構成条件を満たすものを用いることが望ましい。
オーバーラップ加算部１９は、窓がけ部１８で窓がけされた各フレームのフレーム化信号のフレーム境界部分について重ね合わせを行って、雑音が抑圧された出力信号を得る。信号出力端子２０は、オーバーラップ加算部１９で得られた出力信号を出力する。
雑音抑圧装置１０の動作を簡単に説明する。信号入力端子１１に入力信号ｙ(n)が供給され、この入力信号ｙ(n)は、フレーム化部１２に供給される。このフレーム化部１２では、フレーム毎に処理を行うために、入力信号ｙ(n)がフレーム化される。つまり、このフレーム化部１２では、入力信号ｙ(n)が、所定フレーム長、例えばフレーム長がＮfサンプルのフレームに分割される。フレーム毎のフレーム化信号ｙf(u,n)は、窓がけ部１３に、順次供給される。
窓がけ部１３では、後述する高速フーリエ変換部１４で安定したフーリエ係数を得るために、フレーム化信号ｙf(u,n)に分析窓wana(n)による窓がけが行われる。このように窓がけされたフレーム化信号ｙf(u,n)は、高速フーリエ変換部１４に供給される。この高速フーリエ変換部１４では、窓がけされたフレーム化信号ｙf(u,n)に対して、高速フーリエ変換処理が施され、時間領域信号から周波数領域信号に変換される。高速フーリエ変換処理で得られた各フーリエ係数（各周波数スペクトル）は、フーリエ係数修正部１６に供給される。
フレーム化部１２で得られたフレーム毎のフレーム化信号ｙf(u,n)は、雑音抑圧ゲイン生成部１５に供給される。また、高速フーリエ変換部１４で得られたフレーム毎の各フーリエ係数は、雑音抑圧ゲイン生成部１５に供給される。雑音抑圧ゲイン生成部１５では、フレーム毎に、フレーム化信号ｙf(u,n)と各フーリエ係数とに基づいて、各フーリエ係数に対応した雑音抑圧ゲインが生成される。この各フーリエ係数に対応した雑音抑圧ゲインは、フーリエ係数修正部１６に供給される。
フーリエ係数修正部１６では、フレーム毎に、高速フーリエ変換部１４で高速フーリエ変換処理されて得られた各フーリエ係数と、雑音抑圧ゲイン生成部１５で生成された各フーリエ係数に対応した雑音抑圧ゲインの積が取られて係数修正が行われる。つまり、このフーリエ係数修正部１６では、周波数軸上で、雑音を抑圧するためのフィルタ計算が行われる。係数修正された各フーリエ係数は、逆高速フーリエ変換部１７に供給される。
逆高速フーリエ変換部１７では、フレーム毎に、係数修正された各フーリエ係数に対して、逆高速フーリエ変換処理が施され、周波数領域信号が時間領域信号に変換される。逆高速フーリエ変換部１７で得られたフレーム化信号は、窓がけ部１８に供給される。この窓がけ部１８では、フレーム毎に、逆高速フーリエ変換部１７で得られた雑音抑圧されたフレーム化信号に合成窓wsyn(n)による窓がけが行われる。
窓がけ部１８で窓がけされた各フレームのフレーム化信号は、オーバーラップ加算部１９に供給される。このオーバーラップ加算部１９では、各フレームのフレーム化信号のフレーム境界部分について重ね合わせが行われて、雑音が抑圧された出力信号が得られる。そして、この出力信号は、信号出力端子２０に出力される。
［雑音抑圧ゲイン生成部］
雑音抑圧ゲイン生成部１５の詳細を説明する。この雑音抑圧ゲイン生成部１５は、基本的には、上述した非特許文献１などに開示されている雑音抑圧技術を用いて、雑音抑圧ゲインを生成する。最初に、この雑音抑圧技術の概要を以下に説明する。
この雑音抑圧技術では、第ｕフレーム、ｂ番目の帯域の入力帯域信号をＹ(u,b)としたとき、以下の（３）式に示されるように、雑音抑圧ゲインＧ(u,b)が用いられて、雑音を抑圧した帯域信号Ｘ(u,b)が得られる。雑音抑圧ゲインＧ(u,b)は、事前ＳＮＲ「ξ(u,b)」および事後ＳＮＲ「γ(u,b)」から計算される。
Ｘ(u,b)＝Ｇ(u,b)Ｙ(u,b) ・・・（３）
事後ＳＮＲ「γ(u,b)」は、入力信号の帯域パワーをＢ(u,b)、雑音の推定帯域パワーをＤ(u,b)とするとき、以下の（４）式により計算される。
γ(u,b)＝Ｂ(u,b)／Ｄ(u,b) ・・・（４）
事前ＳＮＲ「ξ(u, b)」は、重み係数（平滑化係数）αを用いて、以下の（５）式により計算される。ここで、Ｐ[・]は、以下の（６）式のように定義される演算子である。
ξ(u,b)＝αＧ２(u-1,b)γ(u-1,b)＋(1-α)Ｐ[γ(u,b)-1] ・・・（５）
Figure 2013148724
雑音抑圧ゲインＧ(u,b)は、事前ＳＮＲ「ξ(u, b)」および事後ＳＮＲ「γ(u,b)」を用いて、以下の（７）式のように計算される。In(x) は第一種変形ベッセル関数である。
Figure 2013148724
雑音抑圧ゲインは事前ＳＮＲおよび事後ＳＮＲの推定値から計算されるため、推定精度は雑音抑圧の適切さに直接的な影響を及ぼす。中でも、雑音の帯域パワー推定値Ｄ(u,b)は、ＳＮＲ推定値全てに影響を及ぼすため、この推定精度向上は装置全体の性能向上を目指すうえで重要な課題となる。
雑音の帯域パワーに推定誤差が無いと仮定した場合も、上述の事前ＳＮＲの計算方法（（５）式参照）において、非特許文献１では、α＝０．９８と固定値で扱うことを推奨しており、速い信号変化に推定が追従できない。結果として、雑音抑圧ゲインＧ(u,b)の推定誤差が生じ、音声の冒頭がつぶれてしまうなど音質劣化の原因となる。一方、追従速度を速めるためにαに小さな値を用いると、こんどはミュージカルノイズと呼ばれる聴感上耳障りな副作用が生じ、音質が劣化してしまうという問題がある。
雑音抑圧ゲイン生成部１５は、基本的には、上述した非特許文献１などに開示されている雑音抑圧技術を用いる。しかし、雑音の帯域パワーを精度良く推定するとともに、信号の状態に応じて適応的な係数変更を行うことで、最適な雑音抑圧ゲインＧ(u,b)の生成を行い得るものとされる。
雑音抑圧ゲイン生成部１５は、帯域分割部２１と、帯域パワー算出部２２と、有声音検出部２３と、有声帯域判定部３５と、非定常雑音判定部３６と、雑音・非雑音判定部２７と、雑音帯域パワー推定部２８を有している。また、雑音抑圧ゲイン生成部１５は、事後ＳＮＲ算出部２９と、α算出部３０と、事前ＳＮＲ算出部３１と、雑音抑圧ゲイン算出部３２と、雑音抑圧ゲイン修正部３３と、フィルタ構成部３４を有している。
帯域分割部２１は、高速フーリエ変換部１４で高速フーリエ変換処理されて得られた各周波数スペクトル（各フーリエ係数）を、所定数Ｎbの周波数帯域、例えば、２５の周波数帯域に分割する。表１は、帯域分割の一例を示している。バンド番号は、各帯域を識別するために付された番号である。各周波数帯域は、人間の聴覚システムが、高域ほど知覚分解能が劣化するという聴覚心理の研究から得られた知見に基づいている。
Figure 2013148724
帯域パワー算出部２２は、帯域分割部２１で分割された帯域毎に、周波数スペクトルから帯域パワーＢ(u,b)を算出する。ここで、(u,b)は、第ｕフレーム、ｂ番目の帯域を示している。帯域パワー算出部２２は、帯域パワーＢ(u,b)を算出する方法として、各周波数スペクトルからパワースペクトルをそれぞれ算出し、周波数レンジ内で最大値を求め、この最大値を代表値としてＢ(u,b)とする方法を用いる。なお、帯域パワー算出部２２は、帯域パワーＢ(u,b)を算出する方法として、各周波数スペクトルからパワースペクトルをそれぞれ算出し、周波数レンジ内における平均値を求め、この平均値を代表値としてＢ(u,b)とする方法を用いてもよい。
有声音検出部２３は、フレーム化部１２で得られたフレーム化信号ｙf(u,n)に基づいて、フレーム毎に、有声音が含まれているか否かを示す有声音フラグＦv(u)を出力する。この有声音検出部２３は、ゼロクロス幅計算部２４、ヒストグラム計算部２５および有声音フラグ算出部２６を有している。
ゼロクロス幅計算部２４は、フレーム化された連続するサンプル間で符号が、例えば正から負、あるいは負から正というように逆転する箇所、または、逆の符号を有するサンプル間で０という値を有するサンプルが存在する箇所をゼロクロス点として検出する。そして、ゼロクロス幅計算部２４は、図５に示すように、隣接するゼロクロス点の間のサンプル数を計算してＬz(0)，Ｌz(1)，・・・，Ｌz(m)のようにゼロクロス幅として記録する。
ヒストグラム計算部２５は、ゼロクロス幅計算部２４からゼロクロス幅Ｌz(p)を受け取って、フレーム内における分布を調べる。例えば、１０サンプルおきに２０領域で統計を取る場合、ヒストグラム計算部２５は、初期値として、Ｈz(q)＝０（０≦ｑ＜２０）とする。そして、ヒストグラム計算部２５は、以下の（８）式のように、ヒストグラムＨz(q)を得る。
Figure 2013148724
有声音フラグ算出部２６は、ヒストグラム計算部２５で得られた度数Ｈz(q)が最大値となるインデックス（階級）ｑpeakを求める。そして、有声音フラグ算出部２６は、インデックスｑpeakの度数Ｈz(q)を、そのインデックスｑpeakのしきい値Ｔh(q)と比較し、以下の（９）式に示すように、有声音フラグＦv(u)を設定する。ここで、各インデックスは、各ゼロクロス幅レンジを示す。
Figure 2013148724
図６（ａ），（ｂ）は、フレーム化信号ｙf(u,n)が音声（非雑音）である場合における、信号波形（各サンプルの振幅）と、ゼロクロス幅のヒストグラムの一例を示している。音声（非雑音）の場合、同様の波形が繰り返され、所定のゼロクロス幅レンジの度数が大きくなる。そのため、Ｈz(q)＞Ｔh(q)となり、有声音フラグＦv(u)は、Ｆv(u)＝１に設定される。ここで、しきい値Ｔh(q)は、ゼロクロス幅レンジ（インデックス）毎に設定されており、ゼロクロス幅の小さなゼロクロス幅レンジに対応するＴh(q)ほど大きな値とされている。
一方、図７（ａ），（ｂ）は、フレーム化信号ｙf(u,n)が雑音である場合における、信号波形（各サンプルの振幅）と、ゼロクロス幅のヒストグラムの一例を示している。雑音の場合、ゼロクロス幅の小さなゼロクロス幅レンジの度数が大きくなる。そのため、Ｈz(q)≦Ｔh(q)となり、有声音フラグＦv(u)は、Ｆv(u)＝０に設定される。
有声帯域判定部３５は、フレーム毎に、有声音検出部２３で得られた有声音フラグＦv(u)と、高速フーリエ変換部１４で高速フーリエ変換処理されて得られた各周波数スペクトル（各フーリエ係数）を用いて、各帯域の有声帯域フラグＰv(u,b)を設定する。有声帯域判定部３５は、第ｕフレームの入力フーリエ係数Ｙ(u,k)の振幅を調べ、帯域毎に、帯域内に音声に由来するスペクトルのピークが存在するか否かを確認し、以下の（１０）式に示すように、有声帯域フラグＰv(u, b)を設定する。
Figure 2013148724
音声由来のピークが存在するかどうかは、例えば以下の（１）、（２）の条件で判定することができる。
（１）有声音フラグＦv(u) がセットされている
（２）フーリエ係数の振幅の極大点の値が帯域内平均値のＭt倍以上である(Ｍtはしきい値)
有声帯域判定部３５は、図８のフローチャートに示す判定処理を、フレーム毎に、各帯域において、実行する。有声帯域判定部３５は、ステップＳＴ２１において、処理を開始し、その後に、ステップＳＴ２２の処理に移る。このステップＳＴ２２において、有声帯域判定部３５は、有声音フラグＦv(u)が０より大きいか否か、つまり有声音フラグＦv(u)がセットされているか否かを判断する。
Ｆv(u)＞０でなく有声音フラグＦv(u)がセットされていないとき、有声帯域判定部３５は、ステップＳＴ２３の処理に進み、Ｐv(u,b)＝０に設定し、ステップＳＴ２４において、処理を終了する。一方、Ｆv(u)＞０であって有声音フラグＦv(u)がセットされているとき、有声帯域判定部３５は、音声由来のピークが存在するか否かを判断するための処理に移る。
有声帯域判定部３５は、ステップＳＴ２５において、有声帯域判定部３５は、ｋ＝Ｋbstart、Ｂs＝０に初期化する。ここで、「Ｋbstart」は帯域内のフーリエ係数の最初の番号であり、「Ｋbend」は帯域内のフーリエ係数の最後の番号である。次に、有声帯域判定部３５は、ステップＳＴ２６において、Ｂs＝Ｂs＋|Ｙ(u,k)|の演算を行うと共に、ｋの値を１だけ増加させる。そして、有声帯域判定部３５は、ステップＳＴ２７において、ｋがＫbendより小さいかを判断する。ｋがＫbendより小さいとき、有声帯域判定部３５は、ステップＳＴ２６に戻り、上述したと同様の処理を繰り返し、帯域内のフーリエ係数Ｙ(u,k)の絶対値和を求める。ｋがＫbendとなるとき、有声帯域判定部３５は、ステップＳＴ２８の処理に移る。
このステップＳＴ２８において、有声帯域判定部３５は、Ｂm＝Ｂs/(Ｋbend−Ｋbstart＋1)の演算を行って、帯域内平均値Ｂmを求める。次に、有声帯域判定部３５は、ステップＳＴ２９において、ｋ＝Ｋbstart＋1とする。そして、有声帯域判定部３５は、ステップＳＴ３０において、フーリエ係数Ｙ(u,k)が極大点であるか否かを判断する。つまり、有声帯域判定部３５は、|Ｙ(u,k-1)|＜|Ｙ(u,k)|かつ、|Ｙ(u,k+1)|＜|Ｙ(u,k)|という極大点の条件を満たすか否かを判断する。
極大点の条件を満たさないとき、有声帯域判定部３５は、ステップＳＴ３１において、ｋを１だけ増加させる。そして、有声帯域判定部３５は、ステップＳＴ３２において、ｋがＫbend−１より小さいか否かを判断する。ｋがＫbend−１以下のとき、有声帯域判定部３５は、ステップＳＴ３０に戻り、次のフーリエ係数Ｙ(u,k)が極大点であるか否かを判断する。ステップＳＴ３２でｋがＫbend−１より大きくなるとき、つまり帯域内に極大点がなかったとき、有声帯域判定部３５は、ステップＳＴ２３の処理に進み、Ｐv(u, b)＝０に設定し、ステップＳＴ２４において、処理を終了する。
ステップＳＴ３０でｋ番目のフーリエ係数Ｙ(u,k)が極大点の条件を満たすとき、有声帯域判定部３５は、ステップＳＴ３３の処理に移る。このステップＳＴ３３において、有声帯域判定部３５は、極大点の値が帯域内平均値ＢmのＭt倍以上であるか否かを判断する。つまり、有声帯域判定部３５は、Ｂm＊Ｍt＜|Ｙ(u,k)|の条件を満たすか否かを判断する。
この条件を満たさないとき、有声帯域判定部３５は、ステップＳＴ２３の処理に進み、Ｐv(u,b)＝０に設定し、ステップＳＴ２４において、処理を終了する。一方、この条件を満たすとき、有声帯域判定部３５は、ステップＳＴ３４の処理に進み、Ｐv(u,b)＝１に設定し、ステップＳＴ２４において、処理を終了する。
図４に戻って、非定常雑音判定部３６は、有声帯域判定部３５でＰv(u,b)＝０と判定された帯域の信号が非定常雑音の特性を持つか否かを判定する。すなわち、非定常雑音判定部３６は、フレーム毎に、有声帯域判定部３５で得られた有声帯域フラグＰv(u,b)と、帯域パワー算出部２２で算出された帯域パワーＢ(u,b)を用いて、非定常雑音フラグＦnsn(u)を出力する。
非定常雑音判定部３６は、まず、現在フレームの帯域パワーＢ(u,b)に関して、対象とするノイズに対応したノイズテンプレートＢN(r,b)を、(１≦ｒ≦Ｎr)の範囲で探索し、最も近いノイズテンプレートＢN(rmin,b)を求める。図９のフローチャートは、そのノイズテンプレートＢN(rmin,b)を求める処理の一例を示している。
非定常雑音判定部３６は、ステップＳＴ４１において、処理を開始し、その後に、ステップＳＴ４２の処理に移る。このステップＳＴ４２において、非定常雑音判定部３６は、ｒ＝１、ｃmin＝＋∞、ｒmin＝０にセットする。また、非定常雑音判定部３６は、ステップＳＴ４３において、ｂ＝１、ｄ＝０、ｐ＝０、ｐＮ＝０にセットする。
次に、非定常雑音判定部３６は、ステップＳＴ４４において、有声帯域フラグＰv(u,b)が０より大きいか否か、つまり有声帯域フラグＰv(u,b)がセットされているか否かを判断する。Ｐv(u,b)＞０でなく有声帯域フラグＰv(u,b)がセットされていないとき、非定常雑音判定部３６は、ステップＳＴ４５の処理に移る。このステップＳＴ４５において、非定常雑音判定部３６は、ｄ＝ｄ＋Ｂ(u,b)・ＢN(r,b)、ｐ＝ｐ＋Ｂ(u,b)・Ｂ(u,b)、ｐN＝ｐN＋ＢN(r,b)・ＢN(r,b)の演算を行う。
このステップＳＴ４５の処理の後、非定常雑音判定部３６は、ステップＳＴ４６の処理に移る。上述のステップＳＴ４４でＰv(u,b)＞０であって有声帯域フラグＰv(u,b)がセットされているときも、このステップＳＴ４６の処理に移る。このステップＳＴ４６において、非定常雑音判定部３６は、ｂを１だけ増加させる。
次に、非定常雑音判定部３６は、ステップＳＴ４７において、ｂ≦Ｎbであるか否かを判断する。ｂ≦Ｎbであるとき、非定常雑音判定部３６は、ステップＳＴ４４の処理に戻って、上述したと同様の処理を繰り返す。一方、ｂ≦Ｎbでないとき、非定常雑音判定部３６は、ステップＳＴ４８の処理に移る。このステップＳＴ４８において、非定常雑音判定部３６は、c＝ｄ／√（ｐ・ｐＮ）の演算を行う。
次に、非定常雑音判定部３６は、ステップＳＴ４９において、ｃ＜ｃminであるか否かを判断する。ｃ＜ｃminであるとき、非定常雑音判定部３６は、ステップＳＴ５０において、ｃmin＝ｃ、ｒmin＝c,ｒmim＝ｒとする。そして、ステップＳＴ５１において、ｒを１だけ増加させる。ステップＳＴ４９でｃ＜ｃminでないとき、非定常雑音判定部３６は、直ちにステップＳＴ５１に進み、ｒを１だけ増加させる。
次に、非定常雑音判定部３６は、ステップＳ５２において、ｒ≦Ｎrであるか否かを判断する。ｒ≦Ｎrであるとき、非定常雑音判定部３６は、ステップＳＴ４３に戻り、上述したと同様の動作を繰り返す。一方、ｒ≦Ｎrでないとき、非定常雑音判定部３６は、ステップＳＴ５３において、処理を終了する。
上述の図９のフローチャートの処理により、帯域パワーＢ(u,b)に関して、最も近いノイズテンプレートＢN(rmin,b)が求められる。
次に、非定常雑音判定部３６は、当該フレームに非定常雑音が存在するか判定する。現在フレームの前後±Ｓフレームについて、上述で求めたテンプレートＢN(ｒmin,b)と帯域パワーＢ(u+s,b)の相関ｌ(u+s)とゲイン係数gN(u+s)を求める(−Ｓ≦ｓ≦Ｓ)。そして、非定常雑音判定部３６は、以下の（１）、（２）の条件に基づいて判定し、非定常雑音フラグＦnsn(u)を出力する。
（１）相関l(u + s) がlMAX を越えていない
（２）ゲイン係数ｇN(u+s)の分散が閾値ＧNTを越えている
図１０のフローチャートは、その非定常雑音フラグＦnsn(u)の出力する処理の一例を示している。非定常雑音判定部３６は、ステップＳＴ６１において、処理を開始し、その後に、ステップＳＴ６２の処理に移る。このステップＳＴ６２において、非定常雑音判定部３６は、ｓ＝−Ｓにセットする。また、非定常雑音判定部３６は、ステップＳＴ６３において、ｂ＝１、ｄ＝０、ｐ＝０、ｐＮ＝０にセットする。
次に、非定常雑音判定部３６は、ステップＳＴ６４において、有声帯域フラグＰv(u,b)が０より大きいか否か、つまり有声帯域フラグＰv(u,b)がセットされているか否かを判断する。Ｐv(u,b)＞０でなく有声帯域フラグＰv(u,b)がセットされていないとき、非定常雑音判定部３６は、ステップＳＴ６５の処理に移る。このステップＳＴ６５において、非定常雑音判定部３６は、ｄ＝ｄ＋Ｂ(u+s,b)・ＢN(rmin,b)、ｐ＝ｐ＋Ｂ(u+s,b)・Ｂ(u,b)、ｐN＝ｐN＋ＢN(rmin,b)・ＢN(rmin,b)の演算を行う。
このステップＳＴ６５の処理の後、非定常雑音判定部３６は、ステップＳＴ６６の処理に移る。上述のステップＳＴ６４でＰv(u,b)＞０であって有声帯域フラグＰv(u,b)がセットされているときも、このステップＳＴ６６の処理に移る。このステップＳＴ６６において、非定常雑音判定部３６は、ｂを１だけ増加させる。
次に、非定常雑音判定部３６は、ステップＳＴ６７において、ｂ≦Ｎbであるか否かを判断する。ｂ≦Ｎbであるとき、非定常雑音判定部３６は、ステップＳＴ６４の処理に戻って、上述したと同様の処理を繰り返す。一方、ｂ≦Ｎbでないとき、非定常雑音判定部３６は、ステップＳＴ６８の処理に移る。このステップＳＴ６８において、非定常雑音判定部３６は、ｌ＝ｄ／√（ｐ・ｐＮ）、ｇN(u+s)＝√（ｐ・ｐＮ）の演算を行う。
次に、非定常雑音判定部３６は、ステップＳＴ６９において、ｌ＜ｌMAXであるか否かを判断する。ｌ＜ｌMAXであるとき、非定常雑音判定部３６は、ステップＳＴ７０において、ｓを１だけ増加させる。そして、非定常雑音判定部３６は、ステップＳＴ７１において、ｓ≦Ｓであるか否かを判断する。ｓ≦Ｓであるとき、非定常雑音判定部３６は、ステップＳＴ６３に戻り、上述したと同様の動作を繰り返す。一方、ｓ≦Ｓでないとき、ステップＳＴ７２の処理に移る。
このステップＳＴ７２において、非定常雑音判定部３６は、ゲイン係数ｇN(u+s)の分散が閾値ＧNTを越えているか否かを判断する。分散が閾値ＧNTを越えてとき、非定常雑音判定部３６は、ステップＳＴ７３において、Ｆnsn(u)＝１に設定し、その後、ステップＳＴ７４において、処理を終了する。
一方、ステップＳＴ７２で分散が閾値ＧNTを越えていないとき、非定常雑音判定部３６は、ステップＳＴ７５ｓにおいて、Ｆnsn(u)＝０に設定し、その後、ステップＳＴ７４において、処理を終了する。また、上述のステップＳＴ６９でｌ＜ｌMAXでないとき、非定常雑音判定部３６は、ステップＳＴ７５ｓにおいて、Ｆnsn(u)＝０に設定し、その後、ステップＳＴ７４において、処理を終了する。
上述の図１０のフローチャートの処理により、第ｕフレームに非定常雑音が存在するか否かを示す非定常雑音フラグＦnsn(u)の設定が行われる。
図４に戻って、雑音・非雑音判定部２７は、フレーム毎に、各帯域の雑音帯域フラグＦnz(u,b)を設定する。この場合、雑音・非雑音判定部２７は、有声音検出部２３からの有声音フラグＦv(u)、有声帯域判定部３５からの有声帯域フラグＰv(u,b)、非定常雑音判定部３６からの非定常雑音フラグＦnsn(u)および帯域パワー算出部２２からの帯域パワーＢ(u,b)を用いる。雑音・非雑音判定部２７は、図１１のフローチャートに示す判定処理を、フレーム毎に、各帯域において、実行する。
雑音・非雑音判定部２７は、ステップＳＴ１において、判定処理を開始し、システムの初期化を行う。この初期化において、雑音・非雑音判定部２７は、雑音候補フレーム連続カウンタＣn(b)を、Ｃn(b)＝０に初期化しておく。
次に、雑音・非雑音判定部２７は、ステップＳＴ２の処理に移る。このステップＳＴ２において、雑音・非雑音判定部２７は、非定常雑音フラグＦnsn(u)が０より大きいか否か、つまりＦnsn(u)＝１であるか否かを判断する。Ｆnsn(u)＝１でないとき、雑音・非雑音判定部２７は、ステップＳＴ３の処理に移る。
このステップＳＴ３において、雑音・非雑音判定部２７は、雑音・非雑音判定部２７は、有声音フラグＦv(u)が０より大きいか否か、つまり、Ｆv(u)＝１であるか否かを判断する。Ｆv(u)＝１であるとき、つまり現在フレームｕが有声音であるとき、雑音・非雑音判定部２７は、ステップＳＴ４において、雑音候補フレーム連続カウンタＣn(b)をクリアし、Ｃn(b)＝０とする。そして、雑音・非雑音判定部２７は、現在帯域ｂは雑音でないと判定し、ステップＳＴ５において、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝０とし、その後に、ステップＳＴ６において、判定処理を終了する。
ステップＳＴ３でＦv(u)＝０であるとき、つまり現在フレームｕが有声音でないとき、雑音・非雑音判定部２７は、ステップＳＴ７の処理に移る。このステップＳＴ７において、雑音・非雑音判定部２７は、現在フレームｕの帯域パワーＢ(u,b)と、直前フレームｕ−１の帯域パワーＢ(u-1,b)のパワー比を求める。そして、雑音・非雑音判定部２７は、このステップＳＴ７において、パワー比が、低レベル側のしきい値ＴpL(b)と高レベル側のしきい値ＴpH(b)との間に収まっているか否かを判断する。
雑音・非雑音判定部２７は、パワー比がしきい値間に収まっているとき、現在帯域ｂを雑音の候補とし、パワー比がしきい値間に収まっていないとき、現在帯域ｂは雑音でないと判定する。この判定は、雑音信号のパワーは一定であり、逆にパワー変動が激しい信号は雑音でないとの仮定に基づいている。
パワー比がしきい値間に収まっていないとき、つまり、現在帯域ｂは雑音でないと判定するとき、雑音・非雑音判定部２７は、ステップＳＴ４において、雑音候補フレーム連続カウンタＣn(b)をクリアし、Ｃn(b)＝０とする。そして、雑音・非雑音判定部２７は、ステップＳＴ５において、Ｆnz(u,b)＝０とし、その後に、ステップＳＴ６において、判定処理を終了する。
一方、パワー比がしきい値間に収まっているとき、つまり、現在帯域ｂを雑音の候補とするとき、雑音・非雑音判定部２７は、ステップＳＴ８の処理に移る。このステップＳＴ８において、雑音・非雑音判定部２７は、雑音候補フレーム連続カウンタＣn(b)を、１だけカウントアップする。
そして、雑音・非雑音判定部２７は、ステップＳＴ９において、雑音候補フレーム連続カウンタＣn(b)が、しきい値Ｔcを超えたか否かを判定する。Ｃn(b)＞Ｔcでないとき、雑音・非雑音判定部２７は、現在帯域ｂは雑音でないと判定し、ステップＳＴ５において、Ｆnz(k,b)＝０とし、その後に、ステップＳＴ６において、判定処理を終了する。
一方、Ｃn(b)＞Ｔcであるとき、雑音・非雑音判定部２７は、ステップＳＴ１０の処理に移る。このステップＳＴ１０において、雑音・非雑音判定部２７は、現在帯域ｂが雑音（定常雑音）であると判定し、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝１にセットし、その後に、ステップＳＴ６において、判定処理を終了する。
また、ステップＳＴ２でＦnsn(u)＝１であるとき、雑音・非雑音判定部２７は、ステップＳＴ１１の処理に移る。このステップＳＴ１１において、雑音・非雑音判定部２７は、有声帯域フラグＰv(u,b)が０より大きいか否か、つまりＰv(u,b)＝１であるか否かを判断する。
Ｐv(u,b)＝１であるとき、雑音・非雑音判定部２７は、現在帯域ｂは雑音でないと判定し、ステップＳＴ５において、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝０とし、その後に、ステップＳＴ６において、判定処理を終了する。一方、Ｐv(u,b)＝１でないとき、雑音・非雑音判定部２７は、現在帯域ｂは雑音（非定常雑音）であると判定し、ステップＳＴ１２において、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝２とし、その後に、ステップＳＴ６において、判定処理を終了する。
上述の図１１のフローチャートの判定処理では、定常雑音の判定に関しては、有声音検出部２３で得られた有声音フラグＦv(u)を用いてフレーム全体で一つの雑音・非雑音判定を行い、これと帯域毎の判定を組み合わせることで、最終的な判定結果としている。これは、帯域毎の信号の状態を監視して行う判定だけでは不十分な場合があるからである。例えば、帯域パワーの定常性を検出して雑音と判定しようとする場合、特に帯域分割の帯域幅が広いようなケースではトーン性の信号と雑音の区別がつかない。したがって、図１１のフローチャートの判定処理を行うことで、定常雑音の判定に関して帯域毎の雑音判定精度を高めることができる。
図４に戻って、雑音帯域パワー推定部２８は、フレーム毎に、各帯域の雑音帯域パワー推定値Ｄ(u,b)を推定する。雑音帯域パワー推定部２８は、雑音・非雑音判定部２７で設定された雑音帯域フラグＦnz(u,b)に基づいて、雑音の帯域についてのみ、雑音帯域パワー推定値Ｄ(u,b)の更新を行う。つまり、雑音帯域パワー推定部２８は、Ｆnz(u,b)＝１である定常雑音帯域と、Ｆnz(u,b)＝２である非定常雑音帯域とで、雑音帯域パワー推定値Ｄ(u,b)の更新を行う。
雑音帯域パワー推定部２８における雑音帯域パワー推定値Ｄ(k,b)の更新方法の例としては、例えば、以下の（１１）式に示すように、帯域パワーＢ(u,b)を使用し、指数重みμnzを用いて更新する方法が考えられる。この際、雑音帯域パワー推定部２８は、帯域毎に、帯域パワー算出部２２で得られた現在フレームの帯域パワーと１フレーム前に推定された雑音の帯域パワーとを加重加算することにより、現在フレームの推定された雑音のパワーを得るものである。この場合、定常雑音の場合と非定常雑音の場合とでは指数重みμnzの値は異なるものとされる。
Figure 2013148724
定常雑音の場合は、雑音の振幅変動が小さいため、μnz の値が小さくても十分雑音の変化に追従できる。これに対して、非定常雑音の場合は、雑音の振幅変動が大きく、μnz の値が大きいままでは変動に追従できず雑音の推定誤差が大きくなることから、雑音の低減が十分にできなかったり、音声に副作用を及ぼしたりする。そのため、雑音の特性に応じて、指数重みが切り替えられる。すなわち、非定常雑音における現在フレームの帯域パワーの重みは、定常雑音における現在フレームの帯域パワーの重みより大きくされる。
Ｆnz(u,b)＝１で定常雑音の場合、μnz＝μnz1とされる。このμnz1は、例えば、０．９から１．０程度の値で、雑音帯域パワー推定値Ｄ(k,b)が実際の雑音の変化に追従し、なおかつ、聴感上違和感が無い程度に、設定されることが望ましい。また、Ｆnz(u,b)＝２で非定常雑音の場合、μnz＝μnz2とされる。このμnz2は、μnz1よりも小さく、例えば、０．７から０．８程度の値で、比較的小さい値であることが望ましい。また、μnz1およびμnz2は、それぞれ想定している雑音の特性に応じて、雑音の変化に追従し、聴感上違和感が無い値に調節されることが望ましい。
事後ＳＮＲ算出部２９は、以下の（１２）式に基づき、入力信号の帯域パワーＢ(u,b)と、雑音帯域パワーの推定値Ｄ(u,b)を用いて、フレーム毎に、各帯域の事後ＳＮＲ「γ(u,b)」を算出する。なお、この（１２）式は、上述した（４）式と同じものである。この事後ＳＮＲ算出部２９は、ＳＮＲ算出部を構成している。
γ(u,b)＝Ｂ(u,b)／Ｄ(u,b) ・・・（１２）
事前ＳＮＲ算出部３１は、以下の（１３）式に基づき、フレーム毎に、各帯域の事前ＳＮＲ「ξ(u,b)」を算出する。この場合、事前ＳＮＲ算出部３１は、直前フレームおよび現在フレームの事後ＳＮＲ「γ(u-1,b)，γ(u,b)」と、直前フレームの雑音抑圧ゲインＧ′(u-1,b)と、重み係数αを用いる。なお、この（１３）式は、雑音抑圧ゲインＧ(u-1,b)がリミッタ処理による修正後の雑音抑圧ゲインＧ′(u-1,b)に変わっていることを除き、上述した（５）式と同じものである。
ξ(u,b)＝αＧ′２(u-1,b)γ(u-1,b)＋(１-α)Ｐ[γ(u,b)-1] ・・・（１３）
α算出部３０は、上述の（１３）式における重み係数αを、定数ではなく、フレームと周波数帯域で変動する重み係数α(u,b)として、（１４）式に基づいて、算出する。αMAX(b)およびαMIN(b)は、それぞれ、帯域ごとに設定された重み係数α(u,b)の最大値および最小値である。重み係数α(u,b)を（１４）式に基づいて算出する場合、雑音と判定される帯域ｂにあっては、重み係数α(u,b)は最大値αMAX(b)に近づいていき、非雑音と判定される帯域ｂにあっては、最小値αMIN(b)となる。図１２は、重み係数α(u,b)の推移例を示している。
Figure 2013148724
上述の（１３）式におけるαを、上述のα(u,b)を用いた形に書き直すと、以下の（１５）式のようになる。
ξ(u,b)＝α(u-1,b)Ｇ′２(u-1,b)γ(u-1,b)＋(1-α(u,b))Ｐ[γ(u,b)-1]
・・・（１５）
事前ＳＮＲ算出部３１は、上述の（１５）式に基づき、事前ＳＮＲ「ξ(u,b)」の算出を行う。上述した重み係数α(u,b)の算出の仕組みにより、音声など一般的に変化の激しい非雑音に対しては追従が早く、一方、定常性を仮定している雑音に対してはゆっくり追従するよう、事前ＳＮＲ「ξ(u,b)」の計算がされるようになる。この事前ＳＮＲ算出部３１は、ＳＮＲ平滑部を構成している。
雑音抑圧ゲイン算出部３２は、事後ＳＮＲ算出部２９で算出された事後ＳＮＲ「γ(u,b)」と事前ＳＮＲ算出部３１で算出された事前ＳＮＲ「ξ(u,b)」とから、以下の（１６）式に基づいて、フレーム毎に、各帯域の雑音抑圧ゲインＧ(u,b)を算出する。なお、この（１６）式は、上述した（７）式と同じものである。
Figure 2013148724
雑音抑圧ゲイン修正部３３は、雑音抑圧ゲイン算出部３２で算出された雑音抑圧ゲインＧ(u,b)に、帯域毎に予め設定してある雑音抑圧ゲインの下限値ＧMIN(b)に基づいてリミッタをかけ、修正された雑音抑圧ゲインＧ′(u,b)を算出する。以下の（１７）式は、雑音抑圧ゲイン修正部３３におけるリミッタ処理を表している。
Figure 2013148724
この雑音抑圧ゲイン修正部３３は、聴感上の雑音低減量を最大化しつつも、雑音推定における過大見積りに起因して雑音抑圧ゲインが小さくなりすぎないようにするために設けられている。ここで、下限値ＧMIN(b)は、対象音源の性質や聴覚心理に基づいて、帯域別に設定される。例えば、非雑音の信号が音声である場合に、音声信号が含まれる確率の高い帯域については雑音抑圧ゲインの下限値が高めの値に設定される。雑音抑圧ゲインＧ(u,b)が下限値ＧMIN(b)を下回った場合は、この下限値ＧMIN(b)で置き換えられる。これにより、雑音抑圧ゲインＧ(u,b)の誤差があっても聴感上の音質劣化が少なくされる。
フィルタ構成部３４は、雑音抑圧ゲイン修正部３３で修正されたフレーム毎の各帯域の雑音抑圧ゲインＧ′(u,b)から、フレーム毎に、各フーリエ係数に対応した雑音抑圧ゲインを算出し、周波数軸上のフィルタを構成する。算出方法は、帯域分割部２１でフーリエ係数を帯域分割したのと逆写像して得られたものをそのまま用いる単純な方法でもよいが、ゲインが周波数軸上で不連続にならないように上述の方法で得られたものをさらに周波数軸上でスムージングする方法でもよい。
雑音抑圧ゲイン生成部１５の動作を簡単に説明する。高速フーリエ変換部１４でフレーム毎に高速フーリエ変換処理されて得られた各周波数スペクトル（各フーリエ係数）は、帯域分割部２１および有声帯域判定部３５に供給される。帯域分割部２１では、フレーム毎に、各周波数スペクトルが、所定数Ｎb、例えば２５の周波数帯域に分割される（表１参照）。
帯域分割部２１で帯域分割されて得られた各帯域の周波数スペクトルは、フレーム毎に、帯域パワー算出部２２に供給される。この帯域パワー算出部２２では、フレーム毎に、各帯域の帯域パワーＢ(u,b)が算出される。例えば、帯域ｂ内の各周波数スペクトルに対応したパワースペクトルがそれぞれ算出され、その最大値あるいは平均値が帯域パワーＢ(u,b)とされる。この帯域パワーＢ(u,b)は、非定常雑音判定部３６、雑音・非雑音判定部２７、雑音帯域パワー推定部２８および事後ＳＮＲ算出部２９に供給される。
また、フレーム化部１２で得られたフレーム化信号ｙf(u,n)は、有声音検出部２３に供給される。この有声音検出部２３では、フレーム化信号ｙf(u,n)に基づいて、フレーム毎に、有声音が含まれているか否かを示す有声音フラグＦv(u)が得られる。有声音検出部２３では、フレーム全体の雑音・非雑音の判定が行われ、非雑音と判定されるときはＦv(u)＝１とされ、雑音と判定されるときはＦv(u)＝０とされる。ここで、有声音検出部２３における雑音・非雑音の判定は、フレーム化信号ｙf(u,n)に基づいてゼロクロス幅が検出され、このゼロクロス幅のヒストグラムが計算されることで行われる。
また、有声音検出部２３で得られた有声音フラグＦv(u)は、有声帯域判定部３５に供給される。この有声帯域判定部３５では、フレーム毎に、有声音フラグＦv(u)と、高速フーリエ変換部１４で得られた各周波数スペクトル（各フーリエ係数）が用いられ、各帯域の有声帯域フラグＰv(u,b)が設定される。この場合、第ｕフレームの入力フーリエ係数Ｙ(u,k)の振幅が調べられ、帯域毎に、帯域内に音声に由来するスペクトルのピークが存在するか否かが確認されて、有声帯域フラグＰv(u, b)が設定される。
また、有声音検出部２３で得られた有声音フラグＦv(u)および有声帯域判定部３５で得られた有声帯域フラグＰv(u, b)は、非定常雑音判定部３６に供給される。この非定常雑音判定部３６では、有声帯域判定部３５でＰv(u,b)＝０と判定された帯域の信号が非定常雑音の特性を持つか否かが判定される。この場合、まず、現在フレームの帯域パワーＢ(u,b)に関して、対象とするノイズに対応したノイズテンプレートＢN(r,b)が探索されて、最も近いノイズテンプレートＢN(rmin,b)が求められる。
その後、当該フレームに非定常雑音が存在するかが判定される。この場合、現在フレームの前後±Ｓフレームについて、上述で求めたテンプレートＢN(ｒmin,b)と帯域パワーＢ(u+s,b)の相関ｌ(u+s)とゲイン係数gN(u+s)が求められる。そして、相関l(u+s)がlMAXを越えていないか、ゲイン係数ｇN(u+s)の分散が閾値ＧNTを越えているか、の条件に基づいて判定され、非定常雑音フラグＦnsn(u)が出力される。
また、有声音検出部２３で得られたフレーム毎の有声音フラグＦv(u)、有声帯域判定部３５で得られた有声帯域フラグＰv(u, b)および非定常雑音判定部３６で得られた非定常雑音フラグＦnsn(u)は、雑音・非雑音判定部２７に供給される。この雑音・非雑音判定部２７では、フレーム毎に、各フラグと、各帯域の帯域パワーＢ(u,b)が用いられて、各帯域の雑音帯域フラグＦnz(u,b)が設定される（図１１参照）。
この場合、非定常雑音フラグＦnsn(u)が０であり、かつ有声音フラグＦv(u)が１であってフレーム全体として非雑音と判定されているときは、全ての帯域が雑音でないと判定され、全ての帯域でＦnz(u,b)＝０とされる。
また、非定常雑音フラグＦnsn(u)が０であるが、有声音フラグＦv(u)が０であってフレーム全体として雑音と判定されているときは、帯域毎に、帯域パワーの定常性検出により雑音か非雑音かの判定が行われる。帯域パワーに定常性があって、その帯域が雑音候補と判定されるとき、その帯域の雑音候補フレーム連続カウンタＣn(b)がカウントアップされる。そして、そのカウント値がしきい値Ｔｃを超えるとき、その帯域は雑音（定常性）であると判定され、Ｆnz(u,b)＝１とされる。
一方、帯域パワーに定常性がなく、その帯域が非雑音と判定されるとき、Ｆnz(u,b)＝０とされる。また、帯域パワーに定常性があって、その帯域が雑音候補と判定されても、雑音候補フレーム連続カウンタＣn(b)のカウント値がしきい値Ｔｃ以下のとき、その帯域は非雑音であると判定され、Ｆnz(u,b)＝０とされる。
また、非定常雑音フラグＦnsn(u)が１であるが、有声帯域フラグＰv(u,b)が１であるときは、その帯域は雑音でないと判定され、Ｆnz(u,b)＝０とされる。また、非定常雑音フラグＦnsn(u)が１であり、有声帯域フラグＰv(u,b)が０であるときは、その帯域は雑音（非定常雑音）であると判定され、Ｆnz(u,b)＝２とされる。
雑音・非雑音判定部２７でフレーム毎に設定された各帯域の雑音帯域フラグＦnz(u,b)は雑音帯域パワー推定部２８に供給される。また、帯域パワー算出部２２でフレーム毎に算出された各帯域の帯域パワーＢ(u,b)は、雑音帯域パワー推定部２８に供給される。この雑音帯域パワー推定部２８では、フレーム毎に、各帯域の雑音帯域パワー推定値Ｄ(u,b)が推定される。
この雑音帯域パワー推定部２８では、雑音帯域フラグＦnz(u,b)に基づいて、Ｆnz(u,b)＝１，２である帯域、つまり雑音の帯域についてのみ、雑音帯域パワー推定値Ｄ(u,b)の更新が行われる。例えば、帯域パワーＢ(u,b)が使用され、指数重みμnzが用いられて更新される（（１１）式参照）。この場合、定常雑音の場合と非定常雑音の場合とでは指数重みμnzの値は異なるものとされる。
すなわち、Ｆnz(u,b)＝１で定常雑音の場合、μnz＝μnz1とされる。このμnz1は、例えば、０．９から１．０程度の値で、雑音帯域パワー推定値Ｄ(k,b)が実際の雑音の変化に追従し、なおかつ、聴感上違和感が無い程度に、設定される。また、Ｆnz(u,b)＝２で非定常雑音の場合、μnz＝μnz2とされる。このμnz2は、μnz1よりも小さく、例えば、０．７から０．８程度の値で、比較的小さい値に設定される。これにより、非定常雑音における雑音変化への追従速度が定常雑音における雑音変化への追従速度より上げられることから、雑音の低減が十分にできなかったり、音声に副作用を及ぼしたりする、という不都合を回避可能となる。
雑音帯域パワー推定部２８でフレーム毎に推定された各帯域の雑音帯域パワー推定値Ｄ(u,b)は、事後ＳＮＲ算出部２９に供給される。また、帯域パワー算出部２２でフレーム毎に算出された各帯域の帯域パワーＢ(u,b)は、事後ＳＮＲ算出部２９に供給される。事後ＳＮＲ算出部２９では、フレーム毎に、帯域パワーＢ(u,b)と、雑音帯域パワーの推定値Ｄ(u,b)が用いられて、各帯域の事後ＳＮＲ「γ(u,b)」が算出される（（１２）式参照）。
雑音・非雑音判定部２７でフレーム毎に設定された各帯域の雑音帯域フラグＦnz(u,b)は、α算出部３０に供給される。このα算出部３０では、フレーム毎に、各帯域の事前ＳＮＲ「ξ(u,b)」の算出（（１５）式参照）のための重み係数α(u,b)が算出される。重み係数α(u,b)は、雑音と判定される帯域ｂにあっては、最大値αMAX(b)に近づいていくように更新され、非雑音と判定される帯域ｂにあっては、直ちに最小値αMIN(b)とされる（（１４）式、図１２参照）。
事後ＳＮＲ算出部２９でフレーム毎に算出された各帯域の事後ＳＮＲ「γ(u,b)」は、事前ＳＮＲ算出部３１に供給される。また、α算出部３０でフレーム毎に算出された各帯域の重み係数α(u,b)は、事前ＳＮＲ算出部３１に供給される。さらに、雑音抑圧ゲイン修正部３３で修正された直前フレームの各帯域の雑音抑圧ゲインＧ′(u,b)は、事前ＳＮＲ算出部３１に供給される。事前ＳＮＲ算出部３１では、フレーム毎に、各帯域の事前ＳＮＲ「ξ(u, b)」が算出される（（１５）式参照）。この場合、直前フレームおよび現在フレームの事後ＳＮＲ「γ(u?1,b)，γ(u,b)」と、直前フレームの雑音抑圧ゲインＧ′(u-1,b)と、重み係数α(u,b)が用いられる。
上述したように、α算出部３０で算出される各帯域の重み係数α(u,b)は、雑音と判定される帯域ｂにあっては、最大値αMAX(b)に近づいていくように更新され、非雑音と判定される帯域ｂにあっては、直ちに最小値αMIN(b)とされる。そのため、音声など一般的に変化の激しい非雑音に対しては追従が早く、一方、定常性を仮定している雑音に対してはゆっくり追従するよう、事前ＳＮＲ「ξ(u,b)」が計算される。
事後ＳＮＲ算出部２９でフレーム毎に算出された各帯域の事後ＳＮＲ「γ(u,b)」は、雑音抑圧ゲイン算出部３２に供給される。また、事前ＳＮＲ算出部３１でフレーム毎に算出された各帯域の事前ＳＮＲ「ξ(u, b)」は、雑音抑圧ゲイン算出部３２に供給される。雑音抑圧ゲイン算出部３２では、フレーム毎に、事後ＳＮＲ「γ(u,b)」と事前ＳＮＲ「ξ(u,b)」とから、各帯域の雑音抑圧ゲインＧ(u,b)が算出される（（１６）式参照）。
雑音抑圧ゲイン算出部３２でフレーム毎に算出された各帯域の雑音抑圧ゲインＧ(u,b)は、雑音抑圧ゲイン修正部３３に供給される。この雑音抑圧ゲイン修正部３３では、フレーム毎に、各帯域の雑音抑圧ゲインＧ(u,b)に対して、帯域毎に予め設定してある雑音抑圧ゲインの下限値ＧMIN(b)に基づいてリミッタがかけられ、修正された雑音抑圧ゲインＧ′(u,b)が算出される。
雑音抑圧ゲイン修正部３３でフレーム毎に修正された各帯域の雑音抑圧ゲインＧ′(u,b)は、フィルタ構成部３４に供給される。このフィルタ構成部３４では、フレーム毎に、各帯域の雑音抑圧ゲインＧ′(u,b)から、各フーリエ係数に対応した雑音抑圧ゲインが算出される。このようにフィルタ構成部３４でフレーム毎に算出された各フーリエ係数に対応した雑音抑圧ゲインは、雑音抑圧ゲイン生成部１５の出力としてフーリエ係数修正部１６に供給される。
上述したように、図４に示す雑音抑圧装置１０において、雑音抑圧ゲイン生成部１５の雑音・非雑音判定部３６では、帯域毎に、雑音か非雑音の判定がされる他、さらに雑音が定常雑音であるか非定常雑音であるかの判定も行われて、雑音帯域フラグＦnz(u,b)が設定される。そして、雑音帯域パワー推定部２８では、フレーム毎に、各帯域の雑音帯域パワー推定値Ｄ(u,b)が推定されるが、雑音帯域フラグＦnz(u,b)に基づき、雑音の帯域についてのみ、雑音帯域パワー推定値Ｄ(u,b)の更新が行われる。
この場合、非定常雑音の場合の指数重みμnz2は定常雑音の場合の指数重みμnz1より小さく設定される。そのため、非定常雑音における雑音変化への追従速度が定常雑音における雑音変化への追従速度より上げられる。したがって、雑音が非定常雑音であった場合に、その低減が十分にできなかったり、音声に副作用を及ぼしたりする、という不都合を回避可能となる。
また、図４に示す雑音抑圧装置１０において、雑音抑圧ゲイン生成部１５の雑音抑圧ゲイン算出部３２では、事後ＳＮＲ「γ(u,b)」と事前ＳＮＲ「ξ(u,b)」とから、各帯域の雑音抑圧ゲインＧ(u,b)が算出される。また、事前ＳＮＲ算出部３１では、各帯域の事前ＳＮＲ「ξ(u,b)」が算出される。この場合、直前フレームおよび現在フレームの事後ＳＮＲ「γ(u-1,b)，γ(u,b)」と、直前フレームの雑音抑圧ゲインＧ′(u-1,b)と、重み係数α(u,b)が用いられる。
α算出部３０で算出される各帯域の重み係数α(u,b)は、信号の状態に応じて適応的に変化させられる。つまり、重み係数α(u,b)は、雑音と判定される帯域ｂ（Ｆnz(u,b)＝１）にあっては最大値αMAX(b)に近づいていくように更新され、非雑音と判定される帯域ｂ（Ｆnz(u,b)＝０）にあっては直ちに最小値αMIN(b)とされる。そのため、音声など一般的に変化の激しい非雑音に対しては追従が早く、一方、定常性を仮定している雑音に対してはゆっくり追従するよう、事前ＳＮＲ「ξ(u,b)」が計算される。
そのため、雑音抑圧ゲイン生成部１５で算出される各帯域の雑音抑圧ゲインＧ(u,b)の精度（追従性）を高めることができる。したがって、例えば、音声信号の冒頭など信号の変化が激しい箇所で発生する音質劣化を抑制でき、定常雑音の区間など信号の変化がゆるやかな箇所ではミュージカルノイズを抑制でき、音質向上を図ることができる。
また、上述したように、図４に示す雑音抑圧装置１０において、雑音抑圧ゲイン生成部１５の雑音・非雑音判定部２７では、有声音フラグＦv(u)と、各帯域の帯域パワーＢ(u,b)を用いて、各帯域の雑音帯域フラグＦnz(u,b)が設定される。つまり、雑音と非雑音が混合した信号においても非雑音と重複しない帯域の雑音に関しては検出できる。また、雑音帯域パワー推定部２８では、雑音帯域フラグＦnz(u,b)に基づいて、Ｆnz(u,b)＝１，２である帯域、つまり雑音の帯域についてのみ、雑音帯域パワー推定値Ｄ(u,b)の更新が行われる。そのため、雑音帯域パワー推定値Ｄ(u,b)の推定における時間追従性を改善できると共に、推定精度を高めることができる。したがって、結果的に雑音抑圧ゲインの精度を高めることができ、音質向上を図ることができる。
また、上述したように、図４に示す雑音抑圧装置１０において、雑音抑圧ゲイン生成部１５の雑音・非雑音判定部２７では、有声音フラグＦv(u)と、各帯域の帯域パワーＢ(u,b)を用いて、各帯域の雑音帯域フラグＦnz(u,b)が設定される。つまり、雑音・非雑音判定部２７では、有声音フラグＦv(u)によるフレーム全体での雑音・非雑音判定が行われ、これと帯域パワーの定常性検出による帯域毎の判定が組み合わせられることで、最終的な判定結果が得られる。したがって、帯域毎の雑音・非雑音の判定精度を高めることができる。
また、上述したように、図４に示す雑音抑圧装置１０において、雑音抑圧ゲイン生成部１５の雑音抑圧ゲイン修正部３３では、修正された雑音抑圧ゲインＧ′(u,b)が算出される。この場合、各帯域の雑音抑圧ゲインＧ(u,b)に対して、帯域毎に予め設定してある雑音抑圧ゲインの下限値ＧMIN(b)に基づいてリミッタがかけられて、修正が行われる。したがって、聴感上の雑音低減量を最大化しつつも、推定誤差などによる音質劣化を最小限に抑えることができる。
なお、図４に示す雑音抑圧装置１０において、雑音抑圧ゲイン生成部１５の雑音・非雑音判定部２７では、有声音フラグＦv(u)と、各帯域の帯域パワーＢ(u,b)を用いて、各帯域の雑音帯域フラグＦnz(u,b)が設定される。しかし、雑音・非雑音判定部２７では、フレーム毎に、有声音フラグＦv(u)あるいは各帯域の帯域パワーＢ(u,b)のいずれかのみを用いて、各帯域の雑音帯域フラグＦnz(u,b)を設定することも考えられる。
有声音フラグＦv(u)のみを用いて各帯域の雑音帯域フラグＦnz(u,b)を設定する場合、雑音・非雑音判定部２７では、例えば、図１１のフローチャートにおいて、ステップＳＴ７の処理を除いた判定処理が行われる。一方、各帯域の帯域パワーＢ(u,b)のみを用いて各帯域の雑音帯域フラグＦnz(k,b)を設定する場合、雑音・非雑音判定部２７では、例えば、図１１のフローチャートにおいて、ステップＳＴ３の処理を除いた判定処理が行われる。
＜２．第２の実施の形態＞
［雑音抑圧装置］
図１３は、第２の実施の形態としての雑音抑圧装置１０Ｓの構成例を示している。図４に示す雑音抑圧装置１０がモノラル信号の雑音抑圧に適用した場合の構成例であるのに対して、この雑音抑圧装置１０Ｓはステレオ信号の雑音抑圧に適用した場合の構成例である。この図１３において、図４と対応する部分には、同一符号、あるいは同一符号に「Ｌ」、「Ｒ」の文字を添えて示し、適宜、その詳細説明を省略する。ステレオ信号に適用する際には、基本的にはモノラル信号の処理を各チャネルに行えばよい。しかし、ステレオ信号の場合は推定誤差などの要因により処理結果の定位が崩れてしまうなどの弊害が発生する。そのため、ステレオ信号に対応するための対策が施される。
この雑音抑圧装置１０Ｓは、左チャネル（Ｌch）処理系１００Ｌと、右チャネル（Ｒch）処理系１００Ｌと、雑音抑圧ゲイン生成部１５Ｓとから構成されている。左チャネル処理系１００Ｌおよび右チャネル処理系１００Ｌは、それぞれ、図４に示す雑音抑圧装置１０における信号入力端子１１から信号出力端子２０までの処理系と同様に構成されている。
すなわち、左チャネル処理系１００Ｌは、信号入力端子１１Ｌと、フレーム化部１２Ｌと、窓がけ部１３Ｌと、高速フーリエ変換部１４Ｌを有している。また、この左チャネル処理系１００Ｌは、フーリエ係数修正部１６Ｌと、逆高速フーリエ変換部１７Ｌと、窓がけ部１８Ｌと、オーバーラップ加算部１９Ｌと、信号出力端子２０Ｌを有している。
また、右チャネル処理系１００Ｒは、信号入力端子１１Ｒと、フレーム化部１２Ｒと、窓がけ部１３Ｒと、高速フーリエ変換部１４Ｒを有している。また、この右チャネル処理系１００Ｒは、フーリエ係数修正部１６Ｒと、逆高速フーリエ変換部１７Ｒと、窓がけ部１８Ｒと、オーバーラップ加算部１９Ｒと、信号出力端子２０Ｒを有している。
雑音抑圧ゲイン生成部１５Ｓは、フレーム毎に、左チャネル処理系１００Ｌの各フーリエ係数に対応した雑音抑圧ゲインと、右チャネル処理系１００Ｒの各フーリエ係数に対応した雑音抑圧ゲインを生成する。この雑音抑圧ゲイン生成部１５Ｓは、左チャネル処理系１００Ｌおよび右チャネル処理系１００Ｒの各フーリエ係数に対応した雑音抑圧ゲインＧfL(u,f)，ＧfR(u,f)を生成する。この場合、雑音抑圧ゲイン生成部１５Ｓは、フレーム化信号および各フーリエ係数（各周波数スペクトル）に基づいて、それぞれのチャネルの雑音抑圧ゲインＧfL(u,f)，ＧfR(u,f)を生成する。この雑音抑圧ゲイン生成部１５Ｓの詳細については、後述する。
雑音抑圧装置１０Ｓの動作を簡単に説明する。左チャネル処理系１００Ｌにおいて、信号入力端子１１Ｌに左チャネルの入力信号ｙL(n)が供給され、この入力信号ｙL(n)は、フレーム化部１２Ｌに供給される。このフレーム化部１２Ｌでは、フレーム毎に処理を行うために、入力信号ｙL(n)がフレーム化される。つまり、このフレーム化部１２Ｌでは、入力信号ｙL(n)が、所定フレーム長、例えばフレーム長がＮfサンプルのフレームに分割される。フレーム毎のフレーム化信号ｙfL(u,n)は、窓がけ部１３Ｌに、順次供給される。
窓がけ部１３Ｌでは、後述する高速フーリエ変換部１４Ｌで安定したフーリエ係数を得るために、フレーム化信号ｙfL(u,n)に分析窓wana(n)による窓がけが行われる。このように窓がけされたフレーム化信号ｙfL(u,n)は、高速フーリエ変換部１４Ｌに供給される。この高速フーリエ変換部１４Ｌでは、窓がけされたフレーム化信号ｙfL(u,n)に対して、高速フーリエ変換処理が施され、時間領域信号から周波数領域信号に変換される。高速フーリエ変換処理で得られた各フーリエ係数（各周波数スペクトル）ＹfL(u,f)は、フーリエ係数修正部１６Ｌに供給される。なお、(u,f)は、ｕ番目のフレームのｆ番目の周波数を示している。
また、右チャネル処理系１００Ｒにおいて、信号入力端子１１Ｒに右チャネルの入力信号ｙR(n)が供給され、この入力信号ｙR(n)は、フレーム化部１２Ｒに供給される。このフレーム化部１２Ｒでは、フレーム毎に処理を行うために、入力信号ｙR(n)がフレーム化される。つまり、このフレーム化部１２Ｒでは、入力信号ｙR(n)が、所定フレーム長、例えばフレーム長がＮfサンプルのフレームに分割される。フレーム毎のフレーム化信号ｙfR(u,n)は、窓がけ部１３Ｒに、順次供給される。
窓がけ部１３Ｒでは、後述する高速フーリエ変換部１４Ｒで安定したフーリエ係数を得るために、フレーム化信号ｙfR(u,n)に分析窓wana(n)による窓がけが行われる。このように窓がけされたフレーム化信号ｙfR(u,n)は、高速フーリエ変換部１４Ｒに供給される。この高速フーリエ変換部１４Ｒでは、窓がけされたフレーム化信号ｙfR(u,n)に対して、高速フーリエ変換処理が施され、時間領域信号から周波数領域信号に変換される。高速フーリエ変換処理で得られた各フーリエ係数（各周波数スペクトル）ＹfR(u,f)は、フーリエ係数修正部１６Ｒに供給される。なお、(u,f)は、ｕ番目のフレームのｆ番目の周波数を示している。
フレーム化部１２Ｌ，１２Ｒで得られたフレーム毎のフレーム化信号ｙfL(u,n)，ｙfR(u,n)は、雑音抑圧ゲイン生成部１５Ｓに供給される。また、高速フーリエ変換部１４Ｌ，１４Ｒで得られたフレーム毎のフーリエ係数ＹfL(u,n)，ＹfR(u,n)は、雑音抑圧ゲイン生成部１５Ｓに供給される。雑音抑圧ゲイン生成部１５Ｓでは、フレーム毎に、フレーム化信号ｙfL(u,n)，ｙfR(u,n)とフーリエ係数ＹfL(u,n)，ＹfR(u,n)に基づいて、左右チャネルに共通の、各フーリエ係数に対応した雑音抑圧ゲインが生成される。
また、左チャネル処理系１００Ｌにおいて、フーリエ係数修正部１６Ｌでは、フレーム毎に、高速フーリエ変換部１４Ｌで高速フーリエ変換処理されて得られた各フーリエ係数ＹfL(u,n)の修正が行われる。この場合、各フーリエ係数ＹfL(u,n)と、雑音抑圧ゲイン生成部１５Ｓで生成された各フーリエ係数に対応した雑音抑圧ゲインＧfL(u,f)の積が取られて係数修正が行われる。つまり、このフーリエ係数修正部１６Ｌでは、周波数軸上で、雑音を抑圧するためのフィルタ計算が行われる。係数修正された各フーリエ係数は、逆高速フーリエ変換部１７Ｌに供給される。
逆高速フーリエ変換部１７Ｌでは、フレーム毎に、係数修正された各フーリエ係数に対して、逆高速フーリエ変換処理が施され、周波数領域信号が時間領域信号に変換される。逆高速フーリエ変換部１７Ｌで得られたフレーム化信号は、窓がけ部１８Ｌに供給される。この窓がけ部１８Ｌでは、フレーム毎に、逆高速フーリエ変換部１７Ｌで得られたフレーム化信号に合成窓wsyn(n)による窓がけが行われる。
窓がけ部１８Ｌで窓がけされた各フレームのフレーム化信号は、オーバーラップ加算部１９Ｌに供給される。このオーバーラップ加算部１９Ｌでは、各フレームのフレーム化信号のフレーム境界部分について重ね合わせが行われて、雑音が抑圧された出力信号が得られる。そして、この出力信号は、左チャネル処理系１００Ｌの信号出力端子２０Ｌに出力される。
また、右チャネル処理系１００Ｒにおいて、フーリエ係数修正部１６Ｒでは、フレーム毎に、高速フーリエ変換部１４Ｒで高速フーリエ変換処理されて得られた各フーリエ係数ＹfR(u,n)の修正が行われる。この場合、各フーリエ係数ＹfR(u,n)と、雑音抑圧ゲイン生成部１５Ｓで生成された各フーリエ係数に対応した雑音抑圧ゲインＧfR(u,f)の積が取られて係数修正が行われる。つまり、このフーリエ係数修正部１６Ｒでは、周波数軸上で、雑音を抑圧するためのフィルタ計算が行われる。係数修正された各フーリエ係数は、逆高速フーリエ変換部１７Ｒに供給される。
逆高速フーリエ変換部１７Ｒでは、フレーム毎に、係数修正された各フーリエ係数に対して、逆高速フーリエ変換処理が施され、周波数領域信号が時間領域信号に変換される。逆高速フーリエ変換部１７Ｒで得られたフレーム化信号は、窓がけ部１８Ｒに供給される。この窓がけ部１８Ｒでは、フレーム毎に、逆高速フーリエ変換部１７Ｒで得られたフレーム化信号に合成窓wsyn(n)による窓がけが行われる。
窓がけ部１８Ｒで窓がけされた各フレームのフレーム化信号は、オーバーラップ加算部１９Ｒに供給される。このオーバーラップ加算部１９Ｒでは、各フレームのフレーム化信号のフレーム境界部分について重ね合わせが行われて、雑音が抑圧された出力信号が得られる。そして、この出力信号は、右チャネル処理系１００Ｒの信号出力端子２０Ｒに出力される。
［雑音抑圧ゲイン生成部］
雑音抑圧ゲイン生成部１５Ｓの詳細を説明する。図１４は、雑音抑圧ゲイン生成部１５Ｓの構成例を示している。この図１４において、図４と対応する部分には、同一符号、あるいは同一符号に「Ｌ」、「Ｒ」、「Ｓ」の文字を添えて示し、適宜、その詳細説明を省略する。ここで、「Ｌ」は左チャネル側の処理部分を示し、「Ｒ」は右チャネル側の処理部分を示し、「Ｓ」は左右チャネルの共通処理部分を示す。
雑音抑圧ゲイン生成部１５Ｓは、帯域分割部２１Ｌ，２１Ｒと、帯域パワー算出部２２Ｌ，２２Ｒと、有声音検出部２３Ｌ，２３Ｒと、有声帯域判定部３５Ｌ，３５Ｒと、非定常雑音判定部３６Ｌ，３６Ｒを有している。また、雑音抑圧ゲイン生成部１５Ｓは、雑音・非雑音判定部２７Ｓと、雑音帯域パワー推定部２８Ｌ，２８Ｒを有している。また、雑音抑圧ゲイン生成部１５Ｓは、事後ＳＮＲ算出部２９Ｌ，２９Ｒと、α算出部３０Ｓと、事前ＳＮＲ算出部３１Ｌ，３１Ｒと、雑音抑圧ゲイン算出部３２Ｌ，３２Ｒと、雑音抑圧ゲイン修正部３３Ｌ，３３Ｒと、フィルタ構成部３４Ｌ，３４Ｒを有している。
帯域分割部２１Ｌ，２１Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の帯域分割部２１と同様に構成されている。帯域分割部２１Ｌ，２１Ｒは、高速フーリエ変換部１４Ｌ，１４Ｒで得られた各周波数スペクトル（各フーリエ係数）ＹfL(u,f)，ＹfR(u,f)を、例えば、２５の周波数帯域に分割する（表１参照）。帯域パワー算出部２２Ｌ，２２Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の帯域パワー算出部２２と同様に構成されている。帯域パワー算出部２２Ｌ，２２Ｒは、帯域分割部２１Ｌ，２１Ｒで分割された帯域毎に、周波数スペクトルから帯域パワーＢL(u,b)，ＢR(u,b)を算出する。
有声音検出部２３Ｌ，２３Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の有声音検出部２３と同様に構成されている。有声音検出部２３Ｌ，２３Ｒは、フレーム化部１２Ｌ，１２Ｒで得られたフレーム化信号ｙfL(u,n)，ｙfR(u,n)に基づいて、フレーム毎に、有声音が含まれているか否かを示す有声音フラグＦvL(u)，ＦvR(u)を出力する。
有声帯域判定部３５Ｌ，３５Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の有声帯域判定部３５と同様に構成されている。有声帯域判定部３５Ｌ，３５Ｒは、有声音検出部２３Ｌ，２３Ｒで得られた有声音フラグＦvL(u)，ＦvR(u)と、帯域パワー算出部２２Ｌ，２２Ｒで算出された各帯域の帯域パワーＢL(u,b)，ＢR(u,b)に基づいて、フレーム毎かつ帯域毎に、有声帯域であるか否かを示す有声帯域フラグＰvL(u,b)，ＰvR(u,b)を出力する。
非定常雑音判定部３６Ｌ，３６Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の非定常雑音判定部３６と同様に構成されている。非定常雑音判定部３６Ｌ，３６Ｒは、有声帯域判定部３５Ｌ，３５Ｒで得られた有声帯域フラグＰvL(u,b)，ＰvR(u,b)と、帯域パワー算出部２２Ｌ，２２Ｒで算出された各帯域の帯域パワーＢL(u,b)，ＢR(u,b)に基づいて、フレーム毎に、非定常雑音が存在するフレームであるか否かを示す非定常雑音フラグＦnsnL(u)，ＦnsnR(u)を出力する。
雑音・非雑音判定部２７Ｓは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の雑音・非雑音判定部２７とほぼ同様に構成されている。この雑音・非雑音判定部２７Ｓは、ステレオ対応とされており、フレーム毎に、左右チャネル共通の各帯域の雑音帯域フラグＦnz(u,b)を設定する。
雑音・非雑音判定部２７Ｓは、各帯域の雑音帯域フラグＦnz(u,b)を設定する。この場合、雑音・非雑音判定部２７Ｓは、有声音検出部２３Ｌ，２３Ｒで得られた有声音フラグＦvL(k)，ＦvR(k)と、帯域パワー算出部２２Ｌ，２２Ｒで算出された各帯域の帯域パワーＢL(k,b)，ＢR(k,b)を用いる。さらに、雑音・非雑音判定部２７Ｓは、有声帯域判定部３５Ｌ，３５Ｒで得られた有声帯域フラグＰvL(u,b)，ＰvR(u,b)と、非定常雑音判定部３６Ｌ，３６Ｒで得られた非定常雑音フラグＦnsnL(u)，ＦnsnR(u)を用いる。雑音・非雑音判定部２７Ｓは、図１５のフローチャートに示す判定処理を、フレーム毎に、各帯域において、実行する。
雑音・非雑音判定部２７Ｓは、ステップＳＴ１１１において、判定処理を開始し、システムの初期化を行う。この初期化において、雑音・非雑音判定部２７Ｓは、雑音候補フレーム連続カウンタＣn(b)を、Ｃn(b)＝０に初期化しておく。
次に、雑音・非雑音判定部２７Ｓは、ステップＳＴ１１２の処理に移る。このステップＳＴ１１２において、雑音・非雑音判定部２７Ｓは、非定常雑音フラグＦnsnL(u)，ＦnsnR(u)が０より大きいか否か、つまりＦnsnL(u)，ＦnsnR(u)が１であるか否かを判断する。ＦnsnL(u)＝１かつＦnsnR(u)＝１でないとき、つまり現在フレームｕの少なくとも左右チャネルのいずれかに非定常雑音が存在しないとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１６の処理に移る。雑音・非雑音判定部２７Ｓは、ステップＳＴ１１３の処理に移る。
このステップＳＴ１１３において、雑音・非雑音判定部２７Ｓは、有声音フラグＦvL(u)，ＦvR(u)が０より大きいか否か、つまり、ＦvL(u)，ＦvR(u)が１であるか否かを判断する。ＦvL(u)＝１、かつＦvR(u)＝１であるとき、つまり現在フレームｕが左右チャネル共に有声音であるとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１１４において、雑音候補フレーム連続カウンタＣn(b)をクリアし、Ｃn(b)＝０とする。そして、雑音・非雑音判定部２７Ｓは、現在帯域ｂは雑音でないと判定し、ステップＳＴ１１５において、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝０とし、その後に、ステップＳＴ１１６において、判定処理を終了する。
ステップＳＴ１１３でＦvL(u)＝１、かつＦvR(u)＝１でないとき、つまり現在フレームｕの少なくとも左右チャネルのいずれかが有声音でないとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１１７の処理に移る。このステップＳＴ１１７において、雑音・非雑音判定部２７Ｓは、左チャネル側の現在フレームｕの帯域パワーＢL(u,b)と、直前フレームｕ−１の帯域パワーＢL(u-1,b)のパワー比を求める。また、このステップＳＴ１１７において、雑音・非雑音判定部２７Ｓは、右チャネル側の現在フレームｕの帯域パワーＢR(u,b)と、直前フレームｕ−１の帯域パワーＢR(u-1,b)のパワー比を求める。
そして、雑音・非雑音判定部２７Ｓは、このステップＳＴ１１７において、左右チャネルのパワー比の双方が低レベル側のしきい値ＴpL(b)と高レベル側のしきい値ＴpH(b)との間に収まっているか否かを判断する。すなわち、ＴpL(b)＜ＢL(u,b)/ＢL(u-1,b)＜ＴpH(b)、かつＴpL(b)＜ＢR(u,b)/ＢR(u-1,b)＜ＴpH(b)であるか否かを判断する。
雑音・非雑音判定部２７Ｓは、左右チャネルのパワー比の双方がしきい値間に収まっているとき、現在帯域ｂを雑音の候補とし、左右チャネルのパワー比の双方がしきい値間に収まっていないとき、現在帯域ｂは雑音でないと判定する。この判定は、雑音信号のパワーは一定であり、逆にパワー変動が激しい信号は雑音でないとの仮定に基づいている。
左右チャネルのパワー比の双方がしきい値間に収まっていないとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１３において、雑音候補フレーム連続カウンタＣn(b)をクリアし、Ｃn(b)＝０とする。そして、雑音・非雑音判定部２７Ｓは、現在帯域ｂは雑音でないと判定し、ステップＳＴ１４において、Ｆnz(k,b)＝０とし、その後に、ステップＳＴ１５において、判定処理を終了する。
一方、左右チャネルのパワー比の双方がしきい値間に収まっているとき、つまり、現在帯域ｂを雑音の候補とするとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１１８の処理に移る。このステップＳＴ１１８において、雑音・非雑音判定部２７Ｓは、雑音候補フレーム連続カウンタＣn(b)を、１だけカウントアップする。
そして、雑音・非雑音判定部２７Ｓは、ステップＳＴ１１９において、雑音候補フレーム連続カウンタＣn(b)が、しきい値Ｔcを超えたか否かを判定する。Ｃn(b)＞Ｔcでないとき、雑音・非雑音判定部２７Ｓは、現在帯域ｂは雑音でないと判定し、ステップＳＴ１１５において、Ｆnz(u,b)＝０とし、その後に、ステップＳＴ１１６において、判定処理を終了する。
一方、Ｃn(b)＞Ｔcであるとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１２０の処理に移る。このステップＳＴ１２０において、雑音・非雑音判定部２７Ｓは、現在帯域ｂが雑音であると判定し、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝１にセットし、その後に、ステップＳＴ１１６において、判定処理を終了する。
また、ステップＳＴ１１２でＦnsnL(u)＝１かつＦnsnR(u)＝１であるとき、つまり現在フレームｕの左右チャネルの双方に非定常雑音が存在するとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１２１の処理に移る。このステップＳＴ１１３において、雑音・非雑音判定部２７は、有声帯域フラグＰvL(u,b)，ＰvR(u,b)が０より大きいか否か、つまり、有声帯域フラグＰvL(u,b)，ＰvR(u,b)が１であるか否かを判断する。
ＰvL(u,b)＝１かつＰvR(u,b)＝１であるとき、つまり左右チャネルの双方ともに有声帯域であるとき、雑音・非雑音判定部２７Ｓは、ステップＳＴ１１５において、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝０とし、その後に、ステップＳＴ１１６において、判定処理を終了する。一方、ＰvL(u,b)，ＰvR(u,b)のいずれかが０であるとき、雑音・非雑音判定部２７は、現在帯域ｂは雑音（非定常雑音）であると判定し、ステップＳＴ１２２において、雑音帯域フラグＦnz(u,b)を、Ｆnz(u,b)＝２とし、その後に、ステップＳＴ１１６において、判定処理を終了する。
図１４に戻って、雑音帯域パワー推定部２８Ｌ，２８Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の雑音帯域パワー推定部２８と同様に構成されている。雑音帯域パワー推定部２８Ｌ，２８Ｒは、フレーム毎に、各帯域の雑音帯域パワー推定値ＤL(u,b)，ＤR(u,b)を推定する。雑音帯域パワー推定部２８Ｌ，２８Ｒは、Ｆnz(u,b)＝１である帯域、つまり雑音の帯域についてのみ、雑音帯域パワー推定値ＤL(u,b)，ＤR(u,b)の更新を行う（（１１）式参照）。この場合、雑音帯域パワー推定部２８Ｌ，２８Ｒは、雑音・非雑音判定部２７Ｓで設定された左右チャネル共通の雑音帯域フラグＦnz(u,b)に基づいて、処理を行う。
事後ＳＮＲ算出部２９Ｌ，２９Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の事後ＳＮＲ算出部２９と同様に構成されている。事後ＳＮＲ算出部２９Ｌ，２９Ｒは、フレーム毎に、各帯域の事後ＳＮＲ「γL(u,b)，γR(u,b)」を算出する(（１２）式参照)。この場合、事後ＳＮＲ算出部２９Ｌ，２９Ｒは、入力信号の帯域パワーＢL(u,b)，ＢR(u,b)と、雑音帯域パワーの推定値ＤL(u,b)，ＤR(u,b)を用いる。
事前ＳＮＲ算出部３１Ｌ，３１Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の事前ＳＮＲ算出部３１と同様に構成されている。事前ＳＮＲ算出部３１Ｌ，３１Ｒは、フレーム毎に、各帯域の事前ＳＮＲ「ξL(u,b)，ξR(u,b)」を算出する(（１５）式参照)。
ここで、事前ＳＮＲ算出部３１Ｌは、各帯域の事前ＳＮＲ「ξL(u,b)」を算出する。この場合、事前ＳＮＲ算出部３１Ｌは、直前フレームおよび現在フレームの事後ＳＮＲ「γL(u-1,b)，γL(u,b)」と、直前フレームの雑音抑圧ゲインＧ′L(u-1,b)と、左右チャネル共通の重み係数α(u,b)を用いる。また、事前ＳＮＲ算出部３１Ｒは、各帯域の事前ＳＮＲ「ξR(u,b)」を算出する。この場合、事前ＳＮＲ算出部３１Ｒは、直前フレームおよび現在フレームの事後ＳＮＲ「γR(u-1,b)，γR(u,b)」と、直前フレームの雑音抑圧ゲインＧ′R(u-1,b)と、左右チャネル共通の重み係数α(u,b)を用いる。
α算出部３０Ｓは、図４に示す雑音抑圧装置１０におけるα算出部３０と同様に構成され、事前ＳＮＲ算出部３１Ｌ，３１Ｒで用いられる左右チャネル共通の重み係数α(u,b)を算出する。このα算出部３０Ｓは、定数ではなく、フレームと帯域で変動する重み係数α(u,b)として算出する（（１４）式参照）。この重み係数α(u,b)は、雑音と判定される帯域ｂ（Ｆnz(u,b)＝１，２）にあっては、重み係数α(k,b)は最大値αMAX(b)に近づいていき、非雑音と判定される帯域ｂ（Ｆnz(u,b)＝０）にあっては、最小値αMIN(b)となる。
雑音抑圧ゲイン算出部３２Ｌ，３２Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の雑音抑圧ゲイン算出部３２と同様に構成されている。雑音抑圧ゲイン算出部３２Ｌ，３２Ｒは、フレーム毎に、各帯域の雑音抑圧ゲインＧL(u,b)，ＧR(u,b)を算出する((１６）式参照)。この場合、雑音抑圧ゲイン算出部３２Ｌ，３２Ｒは、事後ＳＮＲ「γL(u,b)，γR(u,b)」と、事前ＳＮＲ「ξL(u,b)，ξR(u,b)」とから、各帯域の雑音抑圧ゲインＧL(u,b)，ＧR(u,b)を算出する。
雑音抑圧ゲイン修正部３３Ｌ，３３Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５の雑音抑圧ゲイン修正部３３と同様に構成されている。雑音抑圧ゲイン修正部３３Ｌ，３３Ｒは、フレーム毎に、雑音抑圧ゲイン算出部３２Ｌ，３２Ｒで算出された雑音抑圧ゲインＧL(u,b)，ＧR(u,b)の修正を行う。つまり、雑音抑圧ゲイン修正部３３Ｌ，３３Ｒは、修正された雑音抑圧ゲインＧ′L(u,b)，Ｇ′R(u,b)を算出する（（１７）式参照）。この場合、雑音抑圧ゲイン修正部３３Ｌ，３３Ｒは、雑音抑圧ゲインＧL(u,b)，ＧR(u,b)に、帯域毎に予め設定してある雑音抑圧ゲインの下限値ＧMIN(b)に基づいてリミッタをかける。
フィルタ構成部３４Ｌ，３４Ｒは、図４に示す雑音抑圧装置１０における雑音抑圧ゲイン生成部１５のフィルタ構成部３４と同様に構成されている。フィルタ構成部３４Ｌ，３４Ｒは、フレーム毎に、雑音抑圧ゲイン修正部３３で修正された各帯域の雑音抑圧ゲインＧ′L(u,b)，Ｇ′R(u,b)から、各フーリエ係数に対応した雑音抑圧ゲインＧfL(u,f)，ＧfR(u,f)を算出する。この場合、フィルタ構成部３４Ｌ，３４Ｒは、周波数軸上のフィルタを構成する。
雑音抑圧ゲイン生成部１５Ｓの動作を簡単に説明する。高速フーリエ変換部１４Ｌ，１４Ｒでフレーム毎に高速フーリエ変換処理されて得られた各周波数スペクトル（各フーリエ係数）ＹfL(u,f)，ＹfR(u,f)は、帯域分割部２１Ｌ，２１Ｒに供給される。この帯域分割部２１Ｌ，２１Ｒでは、フレーム毎に、各周波数スペクトルＹfL(u,f)，ＹfR(u,f)が、所定数Ｎb、例えば２５の周波数帯域に分割される（表１参照）。
帯域分割部２１Ｌ，２１Ｒで帯域分割されて得られた各帯域の周波数スペクトルは、フレーム毎に、帯域パワー算出部２２Ｌ，２２Ｒに供給される。この帯域パワー算出部２２Ｌ，２２Ｒでは、フレーム毎に、各帯域の帯域パワーＢL(u,b)，ＢR(u,b)が算出される。例えば、帯域ｂ内の各周波数スペクトルに対応したパワースペクトルがそれぞれ算出され、その最大値あるいは平均値が帯域パワーＢL(u,b)，ＢR(u,b)とされる。
また、フレーム化部１２Ｌ，１２Ｒで得られたフレーム化信号ｙfL(u,n)，ｙfR(u,n)は、有声音検出部２３Ｌ，２３Ｒに供給される。この有声音検出部２３Ｌ，２３Ｒでは、フレーム化信号ｙfL(u,n)，ｙfR(u,n)に基づいて、フレーム毎に、有声音が含まれているか否かを示す有声音フラグＦvL(u)，ＦvR(u)が得られる。有声音検出部２３Ｌ，２３Ｒでは、フレーム全体の雑音・非雑音の判定が行われ、非雑音と判定されるときはＦvL(u)，ＦvR(u)＝１とされ、雑音と判定されるときはＦvL(u)，ＦvR(u)＝０とされる。ここで、有声音検出部２３Ｌ，２３Ｒにおける雑音・非雑音の判定は、フレーム化信号ｙfL(u,n)，ｙfR(u,n)に基づいてゼロクロス幅が検出され、このゼロクロス幅のヒストグラムが計算されることで行われる。
また、有声音検出部２３Ｌ，２３Ｒで得られた有声音フラグＦvL(u)，ＦvR(u)は、有声帯域判定部３５Ｌ，３５Ｒに供給される。この有声帯域判定部３５Ｌ，３５Ｒでは、フレーム毎に、有声音フラグＦvL(u)，ＦvR(u)と、高速フーリエ変換部１４Ｌ，１４Ｒで得られた各周波数スペクトル（各フーリエ係数）が用いられ、各帯域の有声帯域フラグＰvL(u,b)，ＰvR(u,b)が設定される。この場合、第ｕフレームの入力フーリエ係数ＹfL(u,k)，ＹfR(u,k)の振幅が調べられ、帯域毎に、帯域内に音声に由来するスペクトルのピークが存在するか否かが確認されて、有声帯域フラグＰvL(u,b)，ＰvR(u,b)が設定される。
また、有声帯域判定部３５Ｌ，３５Ｒで得られた有声帯域フラグＰvL(u,b)，ＰvR(u,b)は、非定常雑音判定部３６Ｌ，３６Ｒに供給される。この非定常雑音判定部３６Ｌ，３６Ｒでは、フレーム毎に、有声帯域フラグＰvL(u,b)，ＰvR(u,b)高速フーリエ変換部１４Ｌ，１４Ｒで得られた各周波数スペクトル（各フーリエ係数）が用いられ、非定常雑音フラグＦnsnL(u)，ＦnsnR(u)が設定される。
この場合、有声帯域判定部３５Ｌ，３５ＲでＰvL(u,b)，ＰvR(u,b)＝０とされた帯域の信号が非定常雑音の特性を持つか否かが判定される。この場合、まず、現在フレームの帯域パワーＢL(u,b)，ＢR(u,b)に関して、対象とするノイズに対応したノイズテンプレートＢN(r,b)が探索されて、最も近いノイズテンプレートＢNL(rmin,b)，ＢNR(rmin,b)が求められる。
その後、当該フレームに非定常雑音が存在するかが判定される。この場合、現在フレームの前後±Ｓフレームについて、上述で求めたテンプレートＢNL(rmin,b)，ＢNR(rmin,b)と帯域パワーＢ(u+s,b)の相関ｌ(u+s)とゲイン係数gN(u+s)が求められる。そして、相関l(u+s)がlMAXを越えていないか、ゲイン係数ｇN(u+s)の分散が閾値ＧNTを越えているか、の条件に基づいて判定され、非定常雑音フラグＦnsnL(u)，ＦnsnR(u)が得られる。
有声音検出部２３Ｌ，２３Ｒで得られたフレーム毎の有声音フラグＦvL(k)，ＦvR(k)は、雑音・非雑音判定部２７Ｓに供給される。また、有声音検出部２３Ｌ，２３Ｒで得られた有声音フラグＦvL(u)，ＦvR(u)は、雑音・非雑音判定部２７Ｓに供給される。また、有声帯域判定部３５Ｌ，３５Ｒで得られた有声帯域フラグＰvL(u,b)，ＰvR(u,b)は、雑音・非雑音判定部２７Ｓに供給される。さらに、また、帯域パワー算出部２２Ｌ，２２Ｒで算出されたフレーム毎の各帯域の帯域パワーＢL(u,b)，ＢR(u,b)は、雑音・非雑音判定部２７Ｓに供給される。この雑音・非雑音判定部２７Ｓでは、フレーム毎に、各フラグと、各帯域の帯域パワーＢL(u,b)，ＢR(u,b)を用いて、左右チャネル共通の各帯域の雑音帯域フラグＦnz(u,b)が、設定される（図１５参照）。
この場合、ＦvL(u)＝１、かつＦvR(u)＝１であって、左右チャネルの双方でフレーム全体として非雑音と判定されているときは、全ての帯域が雑音でないと判定され、全ての帯域でＦnz(u,b)＝０とされる。
また、ＦvL(u)＝１、かつＦvR(u)＝１ではなく、左右チャネルの双方でフレーム全体として非雑音と判定されていないときは、帯域毎に、帯域パワーの定常性検出により雑音か非雑音かの判定が行われる。左右チャネルの双方で帯域パワーに定常性があって、その帯域が雑音候補と判定されるとき、その帯域の雑音候補フレーム連続カウンタＣn(b)がカウントアップされる。そして、そのカウント値がしきい値Ｔｃを超えるとき、その帯域は雑音であると判定され、Ｆnz(u,b)＝１とされる。
一方、左右チャネルの双方、あるいはいずれかで帯域パワーに定常性がなく、その帯域が非雑音と判定されるとき、Ｆnz(u,b)＝０とされる。また、左右チャネルの双方で帯域パワーに定常性があって、その帯域が雑音候補と判定されても、雑音候補フレーム連続カウンタＣn(b)のカウント値がしきい値Ｔｃ以下のとき、その帯域は非雑音であると判定され、Ｆnz(u,b)＝０とされる。
また、ＦnsnL(u)＝１、かつＦnsnR(u)＝１でなく、ＰvL(u,b)＝１、かつＰvR(u,b)＝１であるときは、その帯域は雑音でないと判定され、Ｆnz(u,b)＝０とされる。また、ＦnsnL(u)＝１、かつＦnsnR(u)＝１でなく、また、ＰvL(u,b)＝１、かつＰvR(u,b)＝１でないときは、その帯域は雑音（非定常雑音）であると判定され、Ｆnz(u,b)＝２とされる。
雑音・非雑音判定部２７Ｓでフレーム毎に設定された左右チャネル共通の各帯域の雑音帯域フラグＦnz(u,b)は、α算出部３０Ｓに供給される。このα算出部３０Ｓでは、フレーム毎に、各帯域の事前ＳＮＲ「ξL(u,b)，ξR(u,b)」の算出のための、左右チャネル共通の重み係数α(u,b)が算出される（（１４）式参照）。この場合、雑音と判定される帯域ｂ（Ｆnz(u,b)＝１，２）にあっては、重み係数α(k,b)は最大値αMAX(b)に近づいていくように更新され、非雑音と判定される帯域ｂ（Ｆnz(u,b)＝０）にあっては、直ちに最小値αMIN(b)とされる。
雑音・非雑音判定部２７Ｓでフレーム毎に設定された左右チャネル共通の各帯域の雑音帯域フラグＦnz(u,b)は雑音帯域パワー推定部２８Ｌ，２８Ｒに供給される。また、帯域パワー算出部２２Ｌ，２２Ｒでフレーム毎に算出された各帯域の帯域パワーＢL(u,b)，ＢR(u,b)は、雑音帯域パワー推定部２８Ｌ，２８Ｒに供給される。この雑音帯域パワー推定部２８Ｌ，２８Ｒでは、フレーム毎に、各帯域の雑音帯域パワー推定値ＤL(u,b)，ＤR(u,b)が推定される。
この雑音帯域パワー推定部２８Ｌ，２８Ｒでは、雑音帯域フラグＦnz(u,b)に基づいて、Ｆnz(u,b)＝１，２である帯域、つまり雑音の帯域についてのみ、雑音帯域パワー推定値Ｄ(u,b)の更新が行われる。例えば、帯域パワーＢ(u,b)が使用され、指数重みμnzが用いられて更新される（（１１）式参照）。この場合、定常雑音の場合と非定常雑音の場合とでは指数重みμnzの値は異なるものとされる。
すなわち、Ｆnz(u,b)＝１で定常雑音の場合、μnz＝μnz1とされる。このμnz1は、例えば、０．９から１．０程度の値で、雑音帯域パワー推定値Ｄ(k,b)が実際の雑音の変化に追従し、なおかつ、聴感上違和感が無い程度に、設定される。また、Ｆnz(u,b)＝２で非定常雑音の場合、μnz＝μnz2とされる。このμnz2は、μnz1よりも小さく、例えば、０．７から０．８程度の値で、比較的小さい値に設定される。これにより、非定常雑音における雑音変化への追従速度が定常雑音における雑音変化への追従速度より上げられることから、雑音の低減が十分にできなかったり、音声に副作用を及ぼしたりする、という不都合を回避可能となる。
雑音帯域パワー推定部２８Ｌ，２８Ｒでフレーム毎に推定された各帯域の雑音帯域パワー推定値ＤL(u,b)，ＤR(u,b)は、事後ＳＮＲ算出部２９Ｌ，２９Ｒに供給される。また、帯域パワー算出部２２Ｌ，２２Ｒでフレーム毎に算出された各帯域の帯域パワーＢL(u,b)，ＢR(u,b)は、事後ＳＮＲ算出部２９Ｌ，２９Ｒに供給される。事後ＳＮＲ算出部２９Ｌ，２９Ｒでは、フレーム毎に、各帯域の事後ＳＮＲ「γL(u,b)，γR(u,b)」が算出される（（１２）式参照）。この場合、帯域パワーＢL(u,b)，ＢR(u,b)と、雑音帯域パワーの推定値ＤL(u,b)，ＤR(u,b)が用いられる。
事後ＳＮＲ算出部２９Ｌ，２９Ｒでフレーム毎に算出された各帯域の事後ＳＮＲ「γL(u,b)，γR(u,b)」は、事前ＳＮＲ算出部３１Ｌ，３１Ｒに供給される。また、α算出部３０Ｓでフレーム毎に算出された左右チャネル共通の各帯域の重み係数α(u,b)は、事前ＳＮＲ算出部３１Ｌ，３１Ｒに供給される。さらに、雑音抑圧ゲイン修正部２３Ｌ，２３Ｒで修正された直前フレームの各帯域の雑音抑圧ゲインＧ′L(u,b)，Ｇ′R(u,b)は、事前ＳＮＲ算出部３１Ｌ，３１Ｒに供給される。
事前ＳＮＲ算出部３１Ｌ，３１Ｒでは、各帯域の事前ＳＮＲ「ξL(u, b)，ξR(u, b)」が算出される（（１５）式参照）。事前ＳＮＲ算出部３１Ｌでは、フレーム毎に、各帯域の事前ＳＮＲ「ξL(u, b)」が算出される。この場合、直前フレームおよび現在フレームの事後ＳＮＲ「γL(u-1,b)，γL(u,b)」と、直前フレームの雑音抑圧ゲインＧ′L(u-1,b)と、重み係数α(u,b)が用いられる。また、事前ＳＮＲ算出部３１Ｒでは、各帯域の事前ＳＮＲ「ξR(u,b)」が算出される。この場合、フレーム毎に、直前フレームおよび現在フレームの事後ＳＮＲ「γR(u-1,b)，γR(u,b)」と、直前フレームの雑音抑圧ゲインＧ′R(u-1,b)と、重み係数α(u,b)が用いられる。
上述したように、左右チャネル共通の各帯域の重み係数α(u,b)は、雑音と判定される帯域ｂにあっては、最大値αMAX(b)に近づいていくように更新され、非雑音と判定される帯域ｂにあっては、直ちに最小値αMIN(b)とされる。そのため、音声など一般的に変化の激しい非雑音に対しては追従が早く、一方、定常性を仮定している雑音に対してはゆっくり追従するよう、事前ＳＮＲ「ξL(u,b)，ξR(u,b)」が計算される。
事後ＳＮＲ算出部２９Ｌ，２９Ｒでフレーム毎に算出された各帯域の事後ＳＮＲ「γL(u,b)，γR(u,b)」は、雑音抑圧ゲイン算出部３２Ｌ，３２Ｒに供給される。また、事前ＳＮＲ算出部３１Ｌ，３１Ｒでフレーム毎に算出された各帯域の事前ＳＮＲ「ξL(u,b)，ξR(u,b)」は、雑音抑圧ゲイン算出部３２Ｌ，３２Ｒに供給される。雑音抑圧ゲイン算出部３２Ｌ，３２Ｒでは、フレーム毎に、事後ＳＮＲ「γL(u,b)，γR(u,b)」と事前ＳＮＲ「ξL(u,b)，ξR(u,b)」とから、各帯域の雑音抑圧ゲインＧL(u,b)，ＧR(u,b)が算出される（（１６）式参照）。
雑音抑圧ゲイン算出部３２Ｌ，３２Ｒでフレーム毎に算出された各帯域の雑音抑圧ゲインＧL(u,b)，ＧR(u,b)は、雑音抑圧ゲイン修正部３３Ｌ，３３Ｒに供給される。この雑音抑圧ゲイン修正部３３Ｌ，３３Ｒでは、フレーム毎に、修正された雑音抑圧ゲインＧ′L(u,b)，Ｇ′R(u,b)が算出される。この場合、各帯域の雑音抑圧ゲインＧL(u,b)，ＧR(u,b)に対して、帯域毎に予め設定してある雑音抑圧ゲインの下限値ＧMIN(b)に基づいてリミッタがかけられる。
雑音抑圧ゲイン修正部３３Ｌ，３３Ｒでフレーム毎に修正された各帯域の雑音抑圧ゲインＧ′L(u,b)，Ｇ′R(u,b)は、フィルタ構成部３４Ｌ，３４Ｒに供給される。このフィルタ構成部３４Ｌ，３４Ｒでは、フレーム毎に、各帯域の雑音抑圧ゲインＧ′L(u,b)，Ｇ′R(u,b)から、各フーリエ係数に対応した雑音抑圧ゲインＧfL(u,f)，ＧfR(u,f)が算出される。このようにフィルタ構成部３４Ｌ，３４Ｒでフレーム毎に算出された各フーリエ係数に対応した雑音抑圧ゲインは、雑音抑圧ゲイン生成部１５Ｓの出力としてフーリエ係数修正部１６Ｌ，１６Ｒに供給される。
上述したように、図１３に示す雑音抑圧装置１０Ｓはステレオ信号に適用した場合の構成例であるが、雑音抑圧ゲイン生成部１５Ｓは、基本的に、図４に示す雑音抑圧装置１０の雑音抑圧ゲイン生成部１５と同様に構成されている。したがって、図１３に示す雑音抑圧装置１０Ｓにおいても、図４に示す雑音抑圧装置１０と同様の効果を得ることができる。
また、図１３に示す雑音抑圧装置１０Ｓにおいて、雑音抑圧ゲイン生成部１５Ｓの雑音・非雑音判定部２７Ｓでは、フレーム毎に、左右チャネル共通の各帯域の雑音帯域フラグＦnz(u,b)が、設定される。この場合、有声音フラグＦvL(u)，ＦvR(u)と、各帯域の帯域パワーＢL(u,b)，ＢR(u,b)を用いられる。そして、雑音帯域パワー推定部２８Ｌ，２８Ｒでは、フレーム毎に、雑音・非雑音判定部２７Ｓで設定された左右チャネル共通の各帯域の雑音帯域フラグＦnz(u,b)が使用されて、各帯域の雑音帯域パワー推定値ＤL(u,b)，ＤR(u,b)が推定される。
このように、左右チャネルの雑音・非雑音の判定を共通にし、雑音帯域パワー推定部２８Ｌ，２８Ｒで共通の判定結果が用いられる。したがって、図１３に示す雑音抑圧装置１０Ｓにおいて、雑音抑圧ゲイン生成部１５Ｓでは、左右チャネルの雑音帯域パワー推定値ＤL(u,b)，ＤR(u,b)の推定誤差により雑音抑圧ゲインＧL(u,b)，ＧR(u,b)に意図しない振幅差が発生することを抑制できる。そして、これにより、左右チャネルの不整合による定位の破壊を回避できる。
なお、図１３に示す雑音抑圧装置１０Ｓは、ステレオ信号の雑音抑圧に適用した場合の構成例である。詳細説明は省略するが、３チャネル以上の多チャンネル信号の雑音抑圧に適用する雑音抑圧装置も、各チャネルの雑音・非雑音の判定を共通にして、同様に構成できることは勿論である。
＜３．変形例＞
なお、上述実施の形態における雑音抑圧装置１０，１０Ｓは、ハードウェアで構成できる他、同様の処理をソフトウェアで行うこともできる。図１６は、ソフトウェアで処理を行うコンピュータ装置５０の構成例を示している。このコンピュータ装置５０は、ＣＰＵ１８１、ＲＯＭ１８２、ＲＡＭ１８３およびデータ入出力部（データＩ／Ｏ）１８４により構成されている。
ＲＯＭ１８２には、ＣＰＵ１８１の処理プログラム、その他の必要なデータが格納されている。ＲＡＭ１８３は、ＣＰＵ１８１のワークエリアとして機能する。ＣＰＵ１８１は、ＲＯＭ１８２に格納されている処理プログラムを必要に応じて読み出し、読み出した処理プログラムをＲＡＭ１８３に転送して展開し、当該展開された処理プログラムを読み出して、雑音抑圧処理を実行する。
このコンピュータ装置５０においては、入力信号（モノラル信号、ステレオ信号）は、データＩ／Ｏ１８４を介して入力され、ＲＡＭ１８３に蓄積される。このＲＡＭ１８３に蓄積された入力信号に対して、ＣＰＵ１８１により、上述実施の形態と同様の雑音抑圧処理が行われる。そして、処理結果としての雑音が抑圧された出力信号は、データＩ／Ｏ１８４を介して外部に出力される。
また、本技術は、以下のような構成を取ることもできる。
（１）入力信号を所定フレーム長のフレームに分割してフレーム化するフレーム化部と、
上記フレーム化部で得られたフレーム化信号を複数の帯域に分割して帯域分割信号を得る帯域分割部と、
上記帯域分割部で得られた各帯域分割信号から帯域パワーを得る帯域パワー算出部と、
上記フレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定部と、
上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音判定部の判定結果から、各帯域の雑音の帯域パワーを推定する雑音帯域パワー推定部と、
上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する雑音抑圧ゲイン決定部と、
上記帯域分割部で得られた各帯域分割信号に、上記雑音抑圧ゲイン決定部で決定された各帯域の雑音抑圧ゲインを適用して雑音が抑圧された帯域分割信号を得る雑音抑圧部と、
上記雑音抑圧部で得られた各帯域分割信号を帯域合成して雑音抑圧されたフレーム化信号を得る帯域合成部と、
上記帯域合成部で得られた各フレームのフレーム化信号をフレーム合成して雑音を抑圧した出力信号を得るフレーム合成部とを備え、
上記雑音帯域パワー推定部は、
上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
雑音抑圧装置。
（２）上記雑音帯域パワー推定部は、帯域毎に、上記帯域パワー算出部で得られた現在フレームの帯域パワーと１フレーム前に推定された雑音の帯域パワーとを加重加算することにより、現在フレームの推定された雑音のパワーを得、
上記非定常雑音における上記現在フレームの帯域パワーの重みは、上記定常雑音における上記現在フレームの帯域パワーの重みより大きくされる
前記（１）に記載の雑音抑圧装置。
（３）上記雑音判定部は、
所定の帯域が雑音であるか否かの判定に、音声に由来するスペクトルのピークが当該帯域内に存在しないことを条件として用いる
前記（１）または（２）に記載の雑音抑圧装置。
（４）上記雑音抑圧ゲイン決定部は、
帯域毎に、上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーとからＳＮＲを算出するＳＮＲ算出部と、
帯域毎に、上記ＳＮＲ算出部で算出されたＳＮＲを平滑化するＳＮＲ平滑化部とを有し、
上記ＳＮＲ平滑部で平滑化された各帯域のＳＮＲに基づいて、各帯域の雑音抑圧ゲインを決定し、
上記ＳＮＲ平滑化部は、
平滑化の係数を、上記雑音判定部の判定結果および周波数帯域に基づいて変更する
前記（１）から（３）いずれかに記載の雑音抑圧装置。
（５）上記雑音抑圧ゲイン決定部は、
上記ＳＮＲ平滑化部で平滑化された各帯域のＳＮＲと共に、上記ＳＮＲ算出部で算出されたＳＮＲに基づいて、各帯域の雑音抑圧ゲインを決定する
前記（４）に記載の雑音抑圧装置。
（６）上記雑音抑圧ゲイン決定部は、
帯域毎に、現在のフレームの信号の帯域パワーと上記推定された雑音の帯域パワーの比を第１のＳＮＲとし、直前フレームの信号の帯域のパワーと雑音抑圧ゲインを乗算した量と直前フレームの推定された雑音の帯域パワーの比を第２のＳＮＲとし、上記第１のＳＮＲと上記第２のＳＮＲとを用いて雑音抑圧ゲインを決定する
前記（４）に記載の雑音抑圧装置。
（７）上記雑音抑圧ゲイン決定部で決定された雑音抑圧ゲインが予め設定された下限値より小さくなるとき、雑音抑圧ゲインの値を該下限値に修正する雑音抑圧ゲイン修正部をさらに備え、
上記雑音抑圧部は、上記雑音抑圧ゲイン修正部で修正された雑音抑圧ゲインを用いる
前記（４）から（６）のいずれかに記載の雑音抑圧装置。
（８）複数チャネルのそれぞれ所定フレーム長のフレームに分割してフレーム化する複数のフレーム化部と、
上記複数のフレーム化部で得られたフレーム化信号をそれぞれ複数の帯域に分割して帯域分割信号を得る複数の帯域分割部と、
上記複数の帯域分割部で得られた各帯域分割信号からそれぞれ帯域パワーを得る複数の帯域パワー算出部と、
上記複数チャネルのフレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定部と、
上記複数の帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音判定部の判定結果からそれぞれ各帯域の雑音の帯域パワーを推定する複数の雑音帯域パワー推定部と、
上記複数の帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記複数の雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーに基づいてそれぞれ各帯域の雑音抑圧ゲインを決定する複数の雑音抑圧ゲイン決定部と、
上記複数の帯域分割部で得られた各帯域分割信号に、上記複数の雑音抑圧ゲイン決定部で決定された各帯域の雑音抑圧ゲインを適用してそれぞれ雑音が抑圧された帯域分割信号を得る複数の雑音抑圧部と、
上記複数の雑音抑圧部で得られた各帯域分割信号を帯域合成してそれぞれ雑音抑圧されたフレーム化信号を得る複数の帯域合成部と、
上記複数の帯域合成部で得られた各フレームのフレーム化信号をフレーム合成してそれぞれ雑音を抑圧した出力信号を得るフレーム合成部とを備え、
上記雑音帯域パワー推定部は、
上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
雑音抑圧装置。
（９）上記雑音判定部は、各帯域を順次判定帯域とし、該判定帯域について各チャネルで定常雑音および非定常雑音のいずれであるか否かを判定し、全てのチャネルで定常雑音と判定するとき該判定帯域は定常雑音であると判定し、全てのチャネルで非定常雑音と判定するとき該判定帯域は非定常雑音であると判定する
前記（８）に記載の雑音抑圧装置。
（１０）入力信号を所定フレーム長のフレームに分割してフレーム化するフレーム化ステップと、
上記フレーム化ステップで得られたフレーム化信号を複数の帯域に分割して帯域分割信号を得る帯域分割ステップと、
上記帯域分割ステップで得られた各帯域分割信号から帯域パワーを得る帯域パワー算出ステップと、
上記フレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定ステップと、
上記帯域パワー算出ステップで得られた各帯域分割信号の帯域パワーと上記雑音判定ステップの判定結果から、各帯域の雑音の帯域パワーを推定する雑音帯域パワー推定ステップと、
上記帯域パワー算出ステップで得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定ステップで推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する雑音抑圧ゲイン決定ステップと、
上記帯域分割ステップで得られた各帯域分割信号に、上記雑音抑圧ゲイン決定ステップで決定された各帯域の雑音抑圧ゲインを適用して雑音が抑圧された帯域分割信号を得る雑音抑圧ステップと、
上記雑音抑圧ステップで得られた各帯域分割信号を帯域合成して雑音抑圧されたフレーム化信号を得る帯域合成ステップと、
上記帯域合成ステップで得られた各フレームのフレーム化信号をフレーム合成して雑音を抑圧した出力信号を得るフレーム合成ステップとを備え、
上記雑音帯域パワー推定ステップでは、
上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
雑音抑圧方法。
（１１）コンピュータを、
入力信号を所定フレーム長のフレームに分割してフレーム化するフレーム化手段と、
上記フレーム化手段で得られたフレーム化信号を複数の帯域に分割して帯域分割信号を得る帯域分割手段と、
上記帯域分割手段で得られた各帯域分割信号から帯域パワーを得る帯域パワー算出手段と、
上記フレーム化信号の特性に基づいて、帯域毎に定常雑音および非定常雑音のいずれであるか否かを判定する雑音判手段と、
上記帯域パワー算出手段で得られた各帯域分割信号の帯域パワーと上記雑音判定手段の判定結果から、各帯域の雑音の帯域パワーを推定する雑音帯域パワー推定手段と、
上記帯域パワー算出手段で得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定手段で推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する雑音抑圧ゲイン決定手段と、
上記帯域分割手段で得られた各帯域分割信号に、上記雑音抑圧ゲイン決定手段で決定された各帯域の雑音抑圧ゲインを適用して雑音が抑圧された帯域分割信号を得る雑音抑圧手段と、
上記雑音抑圧手段で得られた各帯域分割信号を帯域合成して雑音抑圧されたフレーム化信号を得る帯域合成手段と、
上記帯域合成手段で得られた各フレームのフレーム化信号をフレーム合成して雑音を抑圧した出力信号を得るフレーム合成手段として機能させ、
上記雑音帯域パワー推定手段は、
上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
プログラム。
１０，１０Ｓ・・・雑音抑圧装置
１１，１１Ｌ，１１Ｒ・・・信号入力端子
１２，１２Ｌ，１２Ｒ・・・フレーム化部
１３，１３Ｌ，１３Ｒ・・・窓がけ部
１４，１４Ｌ，１４Ｒ・・・高速フーリエ変換部
１５，１５Ｓ・・・雑音抑圧ゲイン生成部
１６，１６Ｌ，１６Ｒ・・・フーリエ係数修正部
１７，１７Ｌ，１７Ｒ・・・逆高速フーリエ変換部
１８，１８Ｌ，１８Ｒ・・・窓がけ部
１９，１９Ｌ，１９Ｒ・・・オーバーラップ加算部
２０，２０Ｌ，２０Ｒ・・・信号出力端子
２１，２１Ｌ，２１Ｒ・・・帯域分割部
２２，２２Ｌ，２２Ｒ・・・帯域パワー算出部
２３，２３Ｌ，２３Ｒ・・・有声音検出部
２４・・・ゼロクロス幅計算部
２５・・・ヒストグラム計算部
２６・・・有声音フラグ算出部
２７，２７Ｓ・・・雑音・非雑音判定部
２８，２８Ｌ，２８Ｒ・・・雑音帯域パワー推定部
２９，２９Ｌ，２９Ｒ・・・事後ＳＮＲ算出部
３０，３０Ｓ・・・α算出部
３１，３１Ｌ，３１Ｒ・・・事前ＳＮＲ算出部
３２，３２Ｌ，３２Ｒ・・・雑音抑圧ゲイン算出部
３３，３３Ｌ，３３Ｒ・・・雑音抑圧ゲイン修正部
３４，３４Ｌ，３４Ｒ・・・フィルタ構成部
３５，３５Ｌ，３５Ｒ・・・有声帯域判定部
３６，３６Ｌ，３６Ｒ・・・非定常雑音判定部
５０・・・コンピュータ装置
１００Ｌ・・・左チャネル処理系
１００Ｒ・・・右チャネル処理系

Claims (11)

  1. 入力信号を所定フレーム長のフレームに分割してフレーム化するフレーム化部と、
    上記フレーム化部で得られたフレーム化信号を複数の帯域に分割して帯域分割信号を得る帯域分割部と、
    上記帯域分割部で得られた各帯域分割信号から帯域パワーを得る帯域パワー算出部と、
    上記フレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定部と、
    上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音判定部の判定結果から、各帯域の雑音の帯域パワーを推定する雑音帯域パワー推定部と、
    上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する雑音抑圧ゲイン決定部と、
    上記帯域分割部で得られた各帯域分割信号に、上記雑音抑圧ゲイン決定部で決定された各帯域の雑音抑圧ゲインを適用して雑音が抑圧された帯域分割信号を得る雑音抑圧部と、
    上記雑音抑圧部で得られた各帯域分割信号を帯域合成して雑音抑圧されたフレーム化信号を得る帯域合成部と、
    上記帯域合成部で得られた各フレームのフレーム化信号をフレーム合成して雑音を抑圧した出力信号を得るフレーム合成部とを備え、
    上記雑音帯域パワー推定部は、
    上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
    雑音抑圧装置。
  2. 上記雑音帯域パワー推定部は、帯域毎に、上記帯域パワー算出部で得られた現在フレームの帯域パワーと１フレーム前に推定された雑音の帯域パワーとを加重加算することにより、現在フレームの推定された雑音のパワーを得、
    上記非定常雑音における上記現在フレームの帯域パワーの重みは、上記定常雑音における上記現在フレームの帯域パワーの重みより大きくされる
    請求項１に記載の雑音抑圧装置。
  3. 上記雑音判定部は、
    所定の帯域が雑音であるか否かの判定に、音声に由来するスペクトルのピークが当該帯域内に存在しないことを条件として用いる
    請求項１に記載の雑音抑圧装置。
  4. 上記雑音抑圧ゲイン決定部は、
    帯域毎に、上記帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーとからＳＮＲを算出するＳＮＲ算出部と、
    帯域毎に、上記ＳＮＲ算出部で算出されたＳＮＲを平滑化するＳＮＲ平滑化部とを有し、
    上記ＳＮＲ平滑部で平滑化された各帯域のＳＮＲに基づいて、各帯域の雑音抑圧ゲインを決定し、
    上記ＳＮＲ平滑化部は、
    平滑化の係数を、上記雑音判定部の判定結果および周波数帯域に基づいて変更する
    請求項１に記載の雑音抑圧装置。
  5. 上記雑音抑圧ゲイン決定部は、
    上記ＳＮＲ平滑化部で平滑化された各帯域のＳＮＲと共に、上記ＳＮＲ算出部で算出されたＳＮＲに基づいて、各帯域の雑音抑圧ゲインを決定する
    請求項４に記載の雑音抑圧装置。
  6. 上記雑音抑圧ゲイン決定部は、
    帯域毎に、現在のフレームの信号の帯域パワーと上記推定された雑音の帯域パワーの比を第１のＳＮＲとし、直前フレームの信号の帯域のパワーと雑音抑圧ゲインを乗算した量と直前フレームの推定された雑音の帯域パワーの比を第２のＳＮＲとし、上記第１のＳＮＲと上記第２のＳＮＲとを用いて雑音抑圧ゲインを決定する
    請求項４に記載の雑音抑圧装置。
  7. 上記雑音抑圧ゲイン決定部で決定された雑音抑圧ゲインが予め設定された下限値より小さくなるとき、雑音抑圧ゲインの値を該下限値に修正する雑音抑圧ゲイン修正部をさらに備え、
    上記雑音抑圧部は、上記雑音抑圧ゲイン修正部で修正された雑音抑圧ゲインを用いる
    請求項４に記載の雑音抑圧装置。
  8. 複数チャネルのそれぞれ所定フレーム長のフレームに分割してフレーム化する複数のフレーム化部と、
    上記複数のフレーム化部で得られたフレーム化信号をそれぞれ複数の帯域に分割して帯域分割信号を得る複数の帯域分割部と、
    上記複数の帯域分割部で得られた各帯域分割信号からそれぞれ帯域パワーを得る複数の帯域パワー算出部と、
    上記複数チャネルのフレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定部と、
    上記複数の帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記雑音判定部の判定結果からそれぞれ各帯域の雑音の帯域パワーを推定する複数の雑音帯域パワー推定部と、
    上記複数の帯域パワー算出部で得られた各帯域分割信号の帯域パワーと上記複数の雑音帯域パワー推定部で推定された各帯域の雑音の帯域パワーに基づいてそれぞれ各帯域の雑音抑圧ゲインを決定する複数の雑音抑圧ゲイン決定部と、
    上記複数の帯域分割部で得られた各帯域分割信号に、上記複数の雑音抑圧ゲイン決定部で決定された各帯域の雑音抑圧ゲインを適用してそれぞれ雑音が抑圧された帯域分割信号を得る複数の雑音抑圧部と、
    上記複数の雑音抑圧部で得られた各帯域分割信号を帯域合成してそれぞれ雑音抑圧されたフレーム化信号を得る複数の帯域合成部と、
    上記複数の帯域合成部で得られた各フレームのフレーム化信号をフレーム合成してそれぞれ雑音を抑圧した出力信号を得るフレーム合成部とを備え、
    上記雑音帯域パワー推定部は、
    上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
    雑音抑圧装置。
  9. 上記雑音判定部は、各帯域を順次判定帯域とし、該判定帯域について各チャネルで定常雑音および非定常雑音のいずれであるか否かを判定し、全てのチャネルで定常雑音と判定するとき該判定帯域は定常雑音であると判定し、全てのチャネルで非定常雑音と判定するとき該判定帯域は非定常雑音であると判定する
    請求項８に記載の雑音抑圧装置。
  10. 入力信号を所定フレーム長のフレームに分割してフレーム化するフレーム化ステップと、
    上記フレーム化ステップで得られたフレーム化信号を複数の帯域に分割して帯域分割信号を得る帯域分割ステップと、
    上記帯域分割ステップで得られた各帯域分割信号から帯域パワーを得る帯域パワー算出ステップと、
    上記フレーム化信号の特性に基づいて、各帯域が定常雑音および非定常雑音のいずれであるか否かを判定する雑音判定ステップと、
    上記帯域パワー算出ステップで得られた各帯域分割信号の帯域パワーと上記雑音判定ステップの判定結果から、各帯域の雑音の帯域パワーを推定する雑音帯域パワー推定ステップと、
    上記帯域パワー算出ステップで得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定ステップで推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する雑音抑圧ゲイン決定ステップと、
    上記帯域分割ステップで得られた各帯域分割信号に、上記雑音抑圧ゲイン決定ステップで決定された各帯域の雑音抑圧ゲインを適用して雑音が抑圧された帯域分割信号を得る雑音抑圧ステップと、
    上記雑音抑圧ステップで得られた各帯域分割信号を帯域合成して雑音抑圧されたフレーム化信号を得る帯域合成ステップと、
    上記帯域合成ステップで得られた各フレームのフレーム化信号をフレーム合成して雑音を抑圧した出力信号を得るフレーム合成ステップとを備え、
    上記雑音帯域パワー推定ステップでは、
    上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
    雑音抑圧方法。
  11. コンピュータを、
    入力信号を所定フレーム長のフレームに分割してフレーム化するフレーム化手段と、
    上記フレーム化手段で得られたフレーム化信号を複数の帯域に分割して帯域分割信号を得る帯域分割手段と、
    上記帯域分割手段で得られた各帯域分割信号から帯域パワーを得る帯域パワー算出手段と、
    上記フレーム化信号の特性に基づいて、帯域毎に定常雑音および非定常雑音のいずれであるか否かを判定する雑音判手段と、
    上記帯域パワー算出手段で得られた各帯域分割信号の帯域パワーと上記雑音判定手段の判定結果から、各帯域の雑音の帯域パワーを推定する雑音帯域パワー推定手段と、
    上記帯域パワー算出手段で得られた各帯域分割信号の帯域パワーと上記雑音帯域パワー推定手段で推定された各帯域の雑音の帯域パワーに基づいて、各帯域の雑音抑圧ゲインを決定する雑音抑圧ゲイン決定手段と、
    上記帯域分割手段で得られた各帯域分割信号に、上記雑音抑圧ゲイン決定手段で決定された各帯域の雑音抑圧ゲインを適用して雑音が抑圧された帯域分割信号を得る雑音抑圧手段と、
    上記雑音抑圧手段で得られた各帯域分割信号を帯域合成して雑音抑圧されたフレーム化信号を得る帯域合成手段と、
    上記帯域合成手段で得られた各フレームのフレーム化信号をフレーム合成して雑音を抑圧した出力信号を得るフレーム合成手段として機能させ、
    上記雑音帯域パワー推定手段は、
    上記非定常雑音における雑音変化への追従速度を上記定常雑音における雑音変化への追従速度より上げる
    プログラム。
JP2012009240A 2012-01-19 2012-01-19 雑音抑圧装置、雑音抑圧方法およびプログラム Pending JP2013148724A (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2012009240A JP2013148724A (ja) 2012-01-19 2012-01-19 雑音抑圧装置、雑音抑圧方法およびプログラム
US13/719,696 US20130191118A1 (en) 2012-01-19 2012-12-19 Noise suppressing device, noise suppressing method, and program
CN201310009827.4A CN103220440A (zh) 2012-01-19 2013-01-11 噪声抑制装置、噪声抑制方法以及程序

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012009240A JP2013148724A (ja) 2012-01-19 2012-01-19 雑音抑圧装置、雑音抑圧方法およびプログラム

Publications (1)

Publication Number Publication Date
JP2013148724A true JP2013148724A (ja) 2013-08-01

Family

ID=48797948

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012009240A Pending JP2013148724A (ja) 2012-01-19 2012-01-19 雑音抑圧装置、雑音抑圧方法およびプログラム

Country Status (3)

Country Link
US (1) US20130191118A1 (ja)
JP (1) JP2013148724A (ja)
CN (1) CN103220440A (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015034898A (ja) * 2013-08-09 2015-02-19 キヤノン株式会社 音声処理装置及び撮像装置
JP2018072593A (ja) * 2016-10-31 2018-05-10 沖電気工業株式会社 雑音推定装置、プログラム及び方法
JP2018207315A (ja) * 2017-06-05 2018-12-27 キヤノン株式会社 音声処理装置及びその制御方法
JP2020134887A (ja) * 2019-02-25 2020-08-31 富士通株式会社 音信号処理プログラム、音信号処理方法及び音信号処理装置
WO2023228615A1 (ja) * 2022-05-25 2023-11-30 パナソニックＩｐマネジメント株式会社 音声特徴量算出方法、音声特徴量算出装置、及び、口腔機能評価装置

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6361156B2 (ja) * 2014-02-10 2018-07-25 沖電気工業株式会社 雑音推定装置、方法及びプログラム
JP6337519B2 (ja) 2014-03-03 2018-06-06 富士通株式会社 音声処理装置、雑音抑圧方法、およびプログラム
US9721580B2 (en) * 2014-03-31 2017-08-01 Google Inc. Situation dependent transient suppression
EP3155618B1 (en) 2014-06-13 2022-05-11 Oticon A/S Multi-band noise reduction system and methodology for digital audio signals
US10181329B2 (en) * 2014-09-05 2019-01-15 Intel IP Corporation Audio processing circuit and method for reducing noise in an audio signal
US9928978B1 (en) 2015-03-30 2018-03-27 Sean Butler Device monitoring prevention in power systems
US10242689B2 (en) 2015-09-17 2019-03-26 Intel IP Corporation Position-robust multiple microphone noise estimation techniques
DE102016104665A1 (de) * 2016-03-14 2017-09-14 Ask Industries Gmbh Verfahren und Vorrichtung zur Aufbereitung eines verlustbehaftet komprimierten Audiosignals
US10789967B2 (en) 2016-05-09 2020-09-29 Harman International Industries, Incorporated Noise detection and noise reduction
DK3252766T3 (da) 2016-05-30 2021-09-06 Oticon As Audiobehandlingsanordning og fremgangsmåde til estimering af signal-til-støj-forholdet for et lydsignal
US11483663B2 (en) 2016-05-30 2022-10-25 Oticon A/S Audio processing device and a method for estimating a signal-to-noise-ratio of a sound signal
US10861478B2 (en) 2016-05-30 2020-12-08 Oticon A/S Audio processing device and a method for estimating a signal-to-noise-ratio of a sound signal
US10433076B2 (en) 2016-05-30 2019-10-01 Oticon A/S Audio processing device and a method for estimating a signal-to-noise-ratio of a sound signal
AU2017286519B2 (en) * 2016-06-13 2020-05-07 Med-El Elektromedizinische Geraete Gmbh Recursive noise power estimation with noise model adaptation
US10418015B2 (en) * 2017-10-02 2019-09-17 GM Global Technology Operations LLC System for spectral shaping of vehicle noise cancellation
CN107819964B (zh) * 2017-11-10 2021-04-06 Oppo广东移动通信有限公司 提高通话质量的方法、装置、终端和计算机可读存储介质
CN108169533B (zh) * 2017-12-20 2020-08-11 郭伟 基于频谱分割变换的反馈式光纤电流互感器
CN109616135B (zh) * 2018-11-14 2021-08-03 腾讯音乐娱乐科技（深圳）有限公司 音频处理方法、装置及存储介质
US11146607B1 (en) * 2019-05-31 2021-10-12 Dialpad, Inc. Smart noise cancellation
CN111142084B (zh) * 2019-12-11 2023-04-07 中国电子科技集团公司第四十一研究所 一种微动太赫兹谱识别检测算法
CN113516992A (zh) * 2020-08-21 2021-10-19 腾讯科技（深圳）有限公司 一种音频处理方法、装置、智能设备及存储介质
CN112863534B (zh) * 2020-12-31 2022-05-10 思必驰科技股份有限公司 噪声音频消除方法、语音识别方法
JP2022156943A (ja) * 2021-03-31 2022-10-14 富士通株式会社 雑音判定プログラム、雑音判定方法及び雑音判定装置

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL84948A0 (en) * 1987-12-25 1988-06-30 D S P Group Israel Ltd Noise reduction system
FI100840B (fi) * 1995-12-12 1998-02-27 Nokia Mobile Phones Ltd Kohinanvaimennin ja menetelmä taustakohinan vaimentamiseksi kohinaises ta puheesta sekä matkaviestin
US6415253B1 (en) * 1998-02-20 2002-07-02 Meta-C Corporation Method and apparatus for enhancing noise-corrupted speech
US6289309B1 (en) * 1998-12-16 2001-09-11 Sarnoff Corporation Noise spectrum tracking for speech enhancement
JP3454206B2 (ja) * 1999-11-10 2003-10-06 三菱電機株式会社 雑音抑圧装置及び雑音抑圧方法
FI116643B (fi) * 1999-11-15 2006-01-13 Nokia Corp Kohinan vaimennus
JP4282227B2 (ja) * 2000-12-28 2009-06-17 日本電気株式会社 ノイズ除去の方法及び装置
EP1376539B8 (en) * 2001-03-28 2010-12-15 Mitsubishi Denki Kabushiki Kaisha Noise suppressor
KR100400226B1 (ko) * 2001-10-15 2003-10-01 삼성전자주식회사 음성 부재 확률 계산 장치 및 방법과 이 장치 및 방법을이용한 잡음 제거 장치 및 방법
US7593851B2 (en) * 2003-03-21 2009-09-22 Intel Corporation Precision piecewise polynomial approximation for Ephraim-Malah filter
US7133825B2 (en) * 2003-11-28 2006-11-07 Skyworks Solutions, Inc. Computationally efficient background noise suppressor for speech coding and speech recognition
JP4520732B2 (ja) * 2003-12-03 2010-08-11 富士通株式会社 雑音低減装置、および低減方法
US7492889B2 (en) * 2004-04-23 2009-02-17 Acoustic Technologies, Inc. Noise suppression based on bark band wiener filtering and modified doblinger noise estimate
US7454332B2 (en) * 2004-06-15 2008-11-18 Microsoft Corporation Gain constrained noise suppression
EP1760696B1 (en) * 2005-09-03 2016-02-03 GN ReSound A/S Method and apparatus for improved estimation of non-stationary noise for speech enhancement
JP4863713B2 (ja) * 2005-12-29 2012-01-25 富士通株式会社 雑音抑制装置、雑音抑制方法、及びコンピュータプログラム
US7885810B1 (en) * 2007-05-10 2011-02-08 Mediatek Inc. Acoustic signal enhancement method and apparatus
CN101802909B (zh) * 2007-09-12 2013-07-10 杜比实验室特许公司 通过噪声水平估计调整进行的语音增强
US20100008520A1 (en) * 2008-07-09 2010-01-14 Yamaha Corporation Noise Suppression Estimation Device and Noise Suppression Device
US8571231B2 (en) * 2009-10-01 2013-10-29 Qualcomm Incorporated Suppressing noise in an audio signal
JP5528538B2 (ja) * 2010-03-09 2014-06-25 三菱電機株式会社 雑音抑圧装置
JP2012058358A (ja) * 2010-09-07 2012-03-22 Sony Corp 雑音抑圧装置、雑音抑圧方法およびプログラム

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015034898A (ja) * 2013-08-09 2015-02-19 キヤノン株式会社 音声処理装置及び撮像装置
JP2018072593A (ja) * 2016-10-31 2018-05-10 沖電気工業株式会社 雑音推定装置、プログラム及び方法
JP2018207315A (ja) * 2017-06-05 2018-12-27 キヤノン株式会社 音声処理装置及びその制御方法
JP2020134887A (ja) * 2019-02-25 2020-08-31 富士通株式会社 音信号処理プログラム、音信号処理方法及び音信号処理装置
JP7156084B2 (ja) 2019-02-25 2022-10-19 富士通株式会社 音信号処理プログラム、音信号処理方法及び音信号処理装置
WO2023228615A1 (ja) * 2022-05-25 2023-11-30 パナソニックＩｐマネジメント株式会社 音声特徴量算出方法、音声特徴量算出装置、及び、口腔機能評価装置

Also Published As

Publication number Publication date
US20130191118A1 (en) 2013-07-25
CN103220440A (zh) 2013-07-24

Similar Documents

Publication Publication Date Title
JP2013148724A (ja) 雑音抑圧装置、雑音抑圧方法およびプログラム
JP2012058358A (ja) 雑音抑圧装置、雑音抑圧方法およびプログラム
JP3454206B2 (ja) 雑音抑圧装置及び雑音抑圧方法
US7454332B2 (en) Gain constrained noise suppression
EP2546831B1 (en) Noise suppression device
JP3591068B2 (ja) 音声信号の雑音低減方法
EP2362389B1 (en) Noise suppressor
US20120057722A1 (en) Noise removing apparatus and noise removing method
CN110739005A (zh) 一种面向瞬态噪声抑制的实时语音增强方法
CN101593522A (zh) 一种全频域数字助听方法和设备
KR20030040399A (ko) 음성 강화 시스템
WO2019205798A1 (zh) 语音增强的方法、装置及设备
WO2006123721A1 (ja) 雑音抑圧方法およびその装置
JP3588030B2 (ja) 音声区間判定装置及び音声区間判定方法
CN103544961B (zh) 语音信号处理方法及装置
JP2008065090A (ja) ノイズサプレス装置
JPWO2010046954A1 (ja) 雑音抑圧装置および音声復号化装置
CN111508514A (zh) 基于补偿相位谱的单通道语音增强算法
KR20110068637A (ko) 잡음 환경의 입력신호로부터 잡음을 제거하는 방법 및 그 장치
KR20100045933A (ko) 잡음 억압 장치 및 잡음 억압 방법
WO2019205796A1 (zh) 减少频域处理量的方法、装置及设备
JP2004341339A (ja) 雑音抑圧装置
KR101052432B1 (ko) 음성 품질 평가를 위한 발음 종속 조음 보상
JP2008116686A (ja) 雑音抑圧装置
JP2005258158A (ja) ノイズ除去装置