(57)【要約】【課題】  時間的な定常性が高い定常ノイズ成分を、精度良く検出する。【解決手段】  時間周波数変換部１２１が、入力信号ＩＮＤに対して、フレーム期間ごとに時間周波数変換を施して、変換結果を、振幅スペクトル情報ＦＱＤとして、パワースペクトル算出部１２２へ送る。また、時間周波数変換部１２１は、最新フレーム期間における入力信号ＩＮＤの各周波数成分の位相情報ＰＨＤを算出する。そして、パワースペクトル算出部１２２が、振幅スペクトル情報ＦＱＤに基づいて、パワースペクトル分布ＰＳＤを算出する。引き続き、推定部１２４が、周波数毎に、記憶部１２３に記憶された最近の所定期間で得られたパワースペクトル分布において最小値であったパワースペクトルを抽出し、抽出されたパワースペクトルに基づいて、定常ノイズ成分のパワースペクトル分布ＮＳＤを推定する。【選択図】  図１
2012177828.tif 000002
【特許請求の範囲】
【請求項１】
信号成分よりも周波数毎のパワーの時間的な定常性が高い定常ノイズ成分を検出するノイズ検出装置であって、
入力信号に対して、フレーム期間ごとに時間周波数変換を施す時間周波数変換部と；
前記時間周波数変換部による変換結果に基づいて、前記入力信号の前記フレーム期間ごとのパワースペクトル分布を算出するパワースペクトル算出部と；
前記信号成分のいずれの周波数成分も、前記信号成分に連続的には含まれることがない蓋然性が高い期間長を有する最近の期間における、前記パワースペクトル算出部により算出されたパワースペクトル分布を記憶する記憶部と；
周波数毎に、前記記憶部に記憶されたパワースペクトル分布において最小値であったパワースペクトルを抽出し、前記抽出されたパワースペクトルに基づいて、前記定常ノイズ成分のパワースペクトル分布を推定する推定部と；
を備えることを特徴とするノイズ検出装置。
【請求項２】
前記推定部は、
前記抽出されたパワースペクトルが形成する分布が、周波数の変化に応じて連続な変化をしている場合には、前記抽出されたパワースペクトルが形成する分布を、前記定常ノイズ成分のパワースペクトル分布として推定し、
前記抽出されたパワースペクトルが形成する分布が、周波数の変化に応じて不連続な変化をしている場合には、前記抽出されたパワースペクトルが形成する分布を、周波数の変化に応じて滑らかな変化とする連続化処理を施したものを、前記定常ノイズ成分のパワースペクトル分布として推定する、
ことを特徴とする請求項１に記載のノイズ検出装置。
【請求項３】
請求項１又は２に記載のノイズ検出装置と；
前記ノイズ検出装置により検出された定常ノイズ成分のパワースペクトル分布を利用して、前記ノイズ検出装置への入力信号から前記定常ノイズ成分を除去するノイズ低減部と；
を備えることを特徴とするノイズ低減装置。
【請求項４】
前記記憶部は、最新フレーム期間における前記入力信号の時間周波数変換に対応して得られる前記最新フレーム期間における各周波数成分の位相情報を更に記憶し、
前記ノイズ低減部は、前記最新フレーム期間における前記入力信号のパワースペクトル分布から前記検出された定常ノイズ成分のパワースペクトル分布を差し引き、前記差し引いたパワースペクトル分布と、前記位相情報とに基づいて、前記最新フレーム期間における前記入力信号から前記検出された定常ノイズ成分を除去するノイズ低減処理を施した信号を生成する、ことを特徴とする請求項３に記載のノイズ低減装置。
【請求項５】
前記定常ノイズ成分のスペクトル分布は、所定周波数付近にピークを有するスペクトル分布であることが予測され、
前記ノイズ低減部は、前記検出された定常ノイズ成分のパワースペクトル分布が、前記所定周波数付近に有意なピークを有していない場合には、前記ノイズ低減処理を実行しない、
ことを特徴とする請求項３又は４に記載のノイズ低減装置。
【請求項６】
信号成分よりも周波数毎のパワーの時間的な定常性が高い定常ノイズ成分を検出するノイズ検出装置において使用されるノイズ検出方法であって、
入力信号に対して、フレーム期間ごとに時間周波数変換を施す時間周波数変換工程と；
前記時間周波数変換工程における変換結果に基づいて、前記フレーム期間ごとの前記入力信号のパワースペクトル分布を算出するパワースペクトル算出工程と；
前記信号成分のいずれの周波数成分も、前記信号成分に連続的には含まれることがない蓋然性が高い期間長を有する最近の期間における、前記パワースペクトル算出工程において算出されたパワースペクトル分布を記憶部に記憶させる記憶工程と；
周波数毎に、前記記憶部に記憶されたパワースペクトル分布において最小値であったパワースペクトルを抽出し、前記抽出されたパワースペクトルに基づいて、前記定常ノイズ成分のパワースペクトル分布を推定する推定工程と；
を備えることを特徴とするノイズ検出方法。
【請求項７】
請求項６に記載のノイズ検出方法を演算手段に実行させる、ことを特徴とするノイズ検出プログラム。
【請求項８】
請求項７に記載のノイズ検出プログラムが、演算手段により読み取り可能に記録されている、ことを特徴とする記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ノイズ検出装置、ノイズ低減装置、ノイズ検出方法、ノイズ検出プログラム、及び、当該ノイズ検出プログラムが記録された記録媒体に関する。
【背景技術】
【０００２】
従来から、ラジオ受信機や電話機を代表とする受信信号から音声を再生する信号処理装置が広く普及している。こうした装置のうち、自動車等の車両に搭載されるものについては、送電周波数（日本においては、５０Ｈｚ又は６０Ｈｚ）の妨害波を定常的に受信することがあり、ピーク周波数及びピークパワーが音声信号と比べて定常的なノイズ（いわゆるハムノイズ）が信号に混入することがある。また、建物内の固定的な位置に配置されている装置についても、増幅器やＤＡ（Digital to Analogue）変換器、ＡＤ（Analogue to Digital）変換器の回路の不具合により、ピーク周波数又はピークパワーが音声信号と比べて定常的なノイズの影響等を受けることがある。このため、こうした定常的なノイズを検出し、その検出結果に基づいてノイズを低減させて、音声品質を維持するための技術が提案されている。
【０００３】
かかる提案技術の一つとして、周波数領域におけるピーク周波数の定常性及びパワーの定常性を有するノイズを検出し、その検出結果に基づいてノイズ低減を図る技術がある（特許文献１参照：以下、「従来例」と呼ぶ）。この従来例の技術では、入力信号を時間周波数変換してスペクトルを算出し、当該スペクトルから各周波数成分のパワーを算出する。引き続き、当該パワーのピーク周波数が定常的な周波数、又は、当該パワーの大きさが定常的な周波数を抽出し、当該抽出結果に基づいて、ノイズの有無を判定する。そして、当該判定結果に基づいて、定常的なノイズを除去するようになっている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−１５４０９２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上述した従来例の技術では、ノイズのピーク周波数の定常性又はパワーの定常性を有するノイズを検出する。しかしながら、スピーチ等の音声信号における各周波数成分のパワーが大きく変化する場合には、当該音声信号成分が含まれる信号においては、定常的なノイズのピーク周波数の定常性やパワーの定常性が失われてしまうことがある。この結果、定常的なノイズが有する周波数を抽出できず、再生音声品質を維持することができなくなってしまう場合があった。
【０００６】
このため、音声信号成分の周波数帯域における定常性を有するノイズ成分を検出し、当該ノイズ成分を低減することのできる技術が待望されている。かかる要請に応えることが、本発明が解決すべき課題の一つとして挙げられる。
【０００７】
本発明は、上記の事情を鑑みてなされたものであり、信号成分の周波数帯域における周波数毎のパワーの時間的な定常性が高い定常ノイズ成分を検出することができる新たなノイズ検出装置及びノイズ検出方法を提供することを目的とする。また、本発明は、当該定常ノイズ成分を低減し、ＳＮ比の高い信号を得ることができる新たなノイズ低減装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】
請求項１に記載の発明は、信号成分よりも周波数毎のパワーの時間的な定常性が高い定常ノイズ成分を検出するノイズ検出装置であって、入力信号に対して、フレーム期間ごとに時間周波数変換を施す時間周波数変換部と；前記時間周波数変換部による変換結果に基づいて、前記入力信号の前記フレーム期間ごとのパワースペクトル分布を算出するパワースペクトル算出部と；前記信号成分のいずれの周波数成分も、前記信号成分に連続的には含まれることがない蓋然性が高い期間長を有する最近の期間における、前記パワースペクトル算出部により算出されたパワースペクトル分布を記憶する記憶部と；周波数毎に、前記記憶部に記憶されたパワースペクトル分布において最小値であったパワースペクトルを抽出し、前記抽出されたパワースペクトルに基づいて、前記定常ノイズ成分のパワースペクトル分布を推定する推定部と；を備えることを特徴とするノイズ検出装置である。
【０００９】
請求項３に記載の発明は、請求項１又は２に記載のノイズ検出装置と；前記ノイズ検出装置により検出された定常ノイズ成分のパワースペクトル分布を利用して、前記ノイズ検出装置への入力信号から前記定常ノイズ成分を除去するノイズ低減部と；を備えることを特徴とするノイズ低減装置である。
【００１０】
請求項６に記載の発明は、信号成分よりも周波数毎のパワーの時間的な定常性が高い定常ノイズ成分を検出するノイズ検出装置において使用されるノイズ検出方法であって、入力信号に対して、フレーム期間ごとに時間周波数変換を施す時間周波数変換工程と；前記時間周波数変換工程における変換結果に基づいて、前記フレーム期間ごとの前記入力信号のパワースペクトル分布を算出するパワースペクトル算出工程と；前記信号成分のいずれの周波数成分も、前記信号成分に連続的には含まれることがない蓋然性が高い期間長を有する最近の期間における、前記パワースペクトル算出工程において算出されたパワースペクトル分布を記憶部に記憶させる記憶工程と；周波数毎に、前記記憶部に記憶されたパワースペクトル分布において最小値であったパワースペクトルを抽出し、前記抽出されたパワースペクトルに基づいて、前記定常ノイズ成分のパワースペクトル分布を推定する推定工程と；を備えることを特徴とするノイズ検出方法である。
【００１１】
請求項７に記載の発明は、請求項６に記載のノイズ検出方法を演算手段に実行させる、ことを特徴とするノイズ検出プログラムである。
【００１２】
請求項８に記載の発明は、請求項７に記載のノイズ検出プログラムが、演算手段により読み取り可能に記録されている、ことを特徴とする記録媒体である。
【図面の簡単な説明】
【００１３】
【図１】本発明の一実施形態に係るノイズ低減装置の構成を概略的に示すブロック図である。
【図２】図１の装置によるノイズ検出処理及びノイズ低減処理を説明するためのフローチャートである。
【図３】記憶部内に記憶される入力信号のパワースペクトル分布の例を示す図である。
【図４】定常ノイズ成分のパワースペクトル分布の例を示す図である。
【図５】推定された定常ノイズ成分のパワースペクトル分布に不連続がある場合に、連続化処理を施したパワースペクトル分布の例を説明するための図である。
【図６】ノイズ低減部によるノイズ低減処理を説明するためのパワースペクトル分布の例を説明するための図である。
【発明を実施するための形態】
【００１４】
以下、本発明の一実施形態を、図１〜図６を参照して説明する。なお、以下の説明及び図面においては、同一又は同等の要素には同一の符号を付し、重複する説明を省略する。
【００１５】
［構成］
図１には、一実施形態に係るノイズ低減装置１００の概略的な構成がブロック図にて示されている。このノイズ低減装置１００は、外部からの受信信号ＲＣＤを、入力端子１９１を介して受信し、ノイズ検出処理及びノイズ検出結果に基づくノイズ低減処理を施して、出力端子１９２からノイズ低減信号ＮＲＤを出力する装置である。
【００１６】
ここで、ノイズ低減装置１００を備える装置がラジオ受信機である場合には、外部からの受信信号ＲＣＤは、ラジオ放送波である。
【００１７】
図１に示されるように、ノイズ低減装置１００は、生成部１１０と接続されている。そして、生成部１１０は、ノイズ低減装置１００に入力信号ＩＮＤを供給するようになっている。
【００１８】
ここで、ノイズ低減装置１００を備える装置がラジオ受信機である場合には、生成部１１０には、例えば、アンテナ部、選局部（チューナ部）及び検波部が含まれる。そして、検波部による検波結果である検波信号が、入力信号ＩＮＤとして、ノイズ低減装置１００に供給される。
【００１９】
ノイズ低減装置１００は、ノイズ検出装置１２０と、ノイズ低減部１３０とを備えている。
【００２０】
上記のノイズ検出装置１２０は、生成部１１０から送られた入力信号ＩＮＤを受ける。そして、ノイズ検出装置１２０は、入力信号ＩＮＤから、信号成分よりも周波数毎のパワーの時間的な定常性が高い定常ノイズ成分を検出する。かかる機能を有するノイズ検出装置１２０は、時間周波数変換部１２１と、パワースペクトル算出部１２２と、記憶部１２３と、推定部１２４とを備えている。
【００２１】
上記の時間周波数変換部１２１は、生成部１１０から送られた入力信号ＩＮＤを受ける。そして、時間周波数変換部１２１は、入力信号ＩＮＤに対して、所定時間長を有するフレーム期間ごとに時間周波数変換を施す。フレーム期間ごとの入力信号ＩＮＤに対する時間周波変換結果は、各周波数成分の振幅スペクトル情報ＦＱＤとして、パワースペクトル算出部１２２へ送られる。また、時間周波数変換部１２１は、入力信号ＩＮＤの時間周波数変換に対応して得られるフレーム期間における各周波数成分の位相情報ＰＨＤを、記憶部１２３へ送る。
【００２２】
上記のパワースペクトル算出部１２２は、時間周波数変換部１２１から送られた各周波数成分の振幅スペクトル情報ＦＱＤを受ける。そして、パワースペクトル算出部１２２は、当該振幅スペクトル情報ＦＱＤに基づいて、入力信号ＩＮＤのパワースペクトル分布（以下、「入力信号スペクトル分布」とも記す）ＰＳＤを算出する。こうして算出された入力信号スペクトル分布ＰＳＤは、記憶部１２３へ送られるとともに、ノイズ低減部１３０へ送られる。
【００２３】
上記の記憶部１２３は、メモリ素子を備えて構成され、パワースペクトル算出部１２２から送られた入力信号ＩＮＤのパワースペクトル分布ＰＳＤが、最近の所定期間分、記憶される。かかるパワースペクトル分布ＰＳＤの記憶領域は、いわゆるリングバッファ構成とされており、当該最近の所定期間分に含まれるＭ個のフレーム期間ごとのパワースペクトル分布ＰＳＤが記憶されるようになっている。また、記憶部１２３には、最新フレーム期間における各周波数成分の位相情報ＰＨＤが記憶される。この記憶部１２３には、時間周波数変換部１２１、パワースペクトル算出部１２２及び推定部１２４がアクセスできるようになっている。
【００２４】
ここで、「所定期間」は、スピーチ等の音声は一定時間内に有音と無音を繰り返すという特徴に着目して、信号成分のいずれの周波数成分も、信号成分に連続的には含まれることがない蓋然性が高い期間長を有している。また、記憶部１２３は、最近の所定期間分の入力信号ＩＮＤのパワースペクトル分布として、パワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMを記憶しているものとする。
【００２５】
上記の推定部１２４は、記憶部１２３にアクセスし、記憶部１２３に記憶されているパワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMを、スペクトル分布データＳＤＤとして読み取る。そして、推定部１２４は、スペクトル分布データＳＤＤに基づいて、周波数毎に、パワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMにおいて最小値であったパワースペクトルを抽出する。引き続き、推定部１２４は、抽出された周波数毎のパワースペクトルに基づいて、定常ノイズ成分のパワースペクトル分布（以下、「ノイズスペクトル分布」とも記す）ＮＳＤを推定する。こうして推定されたノイズスペクトル分布ＮＳＤは、ノイズ低減部１３０へ送られる。
【００２６】
上記のノイズ低減部１３０は、パワースペクトル算出部１２２から送られた入力信号ＩＮＤのパワースペクトル分布ＰＳＤを受けるとともに、推定部１２４から送られた定常ノイズ成分のパワースペクトル分布ＮＳＤを受ける。また、ノイズ低減部１３０は、記憶部１２３にアクセスし、記憶部１２３に記憶されている最新フレーム期間における各周波数成分の位相情報ＰＨＤを、最新位相情報ＰＨＩとして読み取る。そして、ノイズ低減部１３０は、入力信号スペクトル分布ＰＳＤ、ノイズスペクトル分布ＮＳＤ、及び、最新位相情報ＰＨＩに基づいて、最新のフレーム期間における入力信号ＩＮＤから定常ノイズ成分を除去するノイズ低減処理を施す。ノイズ低減部１３０による処理結果は、ノイズ低減信号ＮＲＤとして、出力端子１９２を介して、外部（ノイズ低減装置１００を備える装置がラジオ受信機である場合には、例えば、復調部）へ出力される。ノイズ低減部１３０によるノイズ低減処理については、後述する。
【００２７】
［動作］
以上のようにして構成されたノイズ低減装置１００の動作について、ノイズ検出処理、及び、当該ノイズ検出処理結果に基づくノイズ低減処理に主に着目して説明する。
【００２８】
受信信号ＲＣＤが、入力端子１９１を介して生成部１１０で受信されると、生成部１１０は、ノイズ低減装置１００で処理可能な形態の入力信号ＩＮＤを生成する。そして、生成部１１０は、生成された入力信号ＩＮＤを、ノイズ検出装置１２０の時間周波数変換部１２１へ送る（図１参照）。
【００２９】
なお、記憶部１２３には、最近の所定期間における入力信号スペクトル分布が記憶されているものとする。
【００３０】
このような状況のもとで、本実施形態のノイズ検出及びノイズ低減処理は、図２に示されるように、まず、ステップＳ１１において、時間周波数変換部１２１が、高速フーリエ変換等の手法を用いて、入力信号ＩＮＤに対して、時間周波数変換を施す。そして、時間周波数変換部１２１は、変換結果として得られた最新フレーム期間における各周波数成分の振幅スペクトル情報ＦＱＤをパワースペクトル算出部１２２へ送る。
【００３１】
また、時間周波数変換部１２１は、入力信号ＩＮＤの時間周波数変換に対応して得られる最新フレーム期間における各周波数成分の位相情報ＰＨＤを算出する。そして、時間周波数変換部１２１は、位相情報ＰＨＤを記憶部１２３へ送る。この結果、記憶部１２３に記憶される位相情報ＰＨＩが、当該新たに算出された最新フレーム期間における各周波数成分の位相情報に更新される。
【００３２】
次に、ステップＳ１２において、パワースペクトル算出部１２２が、最新フレーム期間における各周波数成分の振幅スペクトル情報ＦＱＤに基づいて、入力信号スペクトル分布ＰＳＤを算出する。そして、パワースペクトル算出部１２２は、新たに算出された入力信号スペクトル分布ＰＳＤを記憶部１２３へ送る。この結果、記憶部１２３に記憶されるパワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMが、当該新たに算出された入力信号スペクトル分布ＰＳＤを含む最近の所定期間におけるＭ個のスペクトル分布に更新される。
【００３３】
なお、記憶部１２３内に記憶されたパワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMの例が、図３に示されている。
【００３４】
また、パワースペクトル算出部１２２は、新たに算出された入力信号スペクトル分布ＰＳＤを、ノイズ低減部１３０へ送る。
【００３５】
次いで、ステップＳ１３において、推定部１２４が、記憶部１２３にアクセスして、記憶部１２３に記憶されているパワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMを、スペクトル分布データＳＤＤとして読み取る。そして、推定部１２４は、周波数毎に、パワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMにおいて最小値であったパワースペクトルを抽出し、抽出された周波数毎のパワースペクトルに基づいて、定常ノイズ成分のパワースペクトル分布ＮＳＤを推定する。
【００３６】
図４には、図３に示されるパワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMに基づいて推定された定常ノイズ成分のパワースペクトル分布ＮＳＤの例が示されている。引き続き、推定部１２４は、ノイズスペクトル分布ＮＳＤをノイズ低減部１３０へ送る。この後、処理はステップＳ１４へ進む。
【００３７】
なお、推定部１２４は、推定された定常ノイズ成分のパワースペクトル分布が、図５（Ａ）に示されるように、周波数の変化に応じて不連続に変化している場合には、周波数の変化に応じて滑らかな変化とする連続化処理を施す。そして、推定部１２４は、図５（Ｂ）に示されるように、当該連続化処理を施したものを、定常ノイズ成分のパワースペクトル分布ＮＳＤとして推定する。ここで、記憶部１２３に記憶されているパワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMが得られた最近の所定期間は、上述したように、信号成分のいずれの周波数成分も、信号成分に連続的には含まれることがない蓋然性が高い期間長を有している。このため、推定部１２４による推定結果は、定常ノイズ成分のパワースペクトル分布の精度良い推定結果となっている。
【００３８】
ステップＳ１４では、ノイズ低減部１３０が、まず、最新フレーム期間における入力信号ＩＮＤのパワースペクトル分布ＰＳＤから定常ノイズ成分のパワースペクトル分布ＮＳＤを差し引く。図６（Ａ）には、最新フレーム期間における入力信号ＩＮＤのパワースペクトル分布の例が示され、図６（Ｂ）には、定常ノイズ成分のパワースペクトル分布ＮＳＤの例が示されている。そして、図６（Ｃ）には、入力信号ＩＮＤのパワースペクトル分布ＰＳＤから定常ノイズ成分のパワースペクトル分布ＮＳＤを差し引いたパワースペクトル分布の例が示されている。
【００３９】
そして、ノイズ低減部１３０は、記憶部１２４にアクセスして、最新フレーム期間における入力信号ＩＮＤの最新位相情報ＰＨＩを読み出し、当該差し引いたパワースペクトル分布と読み出した入力信号ＩＮＤの最新位相情報ＰＨＩとに基づいて、逆高速フーリエ変換等の手法を用いて、周波数時間変換を行う。この結果、最新フレーム期間における入力信号ＩＮＤから定常ノイズ成分を除去したノイズ低減信号ＮＳＤが生成されて、出力端子１９２を介して、外部へ出力される。
【００４０】
この後、処理はステップＳ１１へ戻る。以後、ステップＳ１１〜Ｓ１４の処理が繰り返され、入力信号ＩＮＤに対するノイズ検出処理及びノイズ低減処理が行われる。
【００４１】
以上説明したように、本実施形態では、受信信号ＲＣＤを受けた生成部１１０が、ノイズ低減装置１００で処理可能な形態の入力信号ＩＮＤを生成し、ノイズ検出装置１２０の時間周波数変換部１２１へ送る。そして、時間周波数変換部１２１が、フレーム期間における入力信号ＩＮＤに対して、時間周波数変換を施して得られた、最新フレーム期間における各周波数成分の振幅スペクトル情報ＦＱＤをパワースペクトル算出部１２２へ送る。また、時間周波数変換部１２１は、入力信号ＩＮＤの時間周波数変換に対応して得られる最新フレーム期間における各周波数成分の位相情報ＰＨＤを算出し、記憶部１２３に記憶させる。そして、パワースペクトル算出部１２２が、各周波数成分の振幅スペクトル情報ＦＱＤに基づいて、入力信号ＩＮＤのパワースペクトル分布ＰＳＤを算出し、新たに算出された入力信号ＩＮＤのパワースペクトル分布ＰＳＤを含む最近の所定期間における各フレーム期間で得られたＭ個のパワースペクトル分布ＰＳＤを記憶部１２３に記憶させる。
【００４２】
引き続き、推定部１２４が、記憶部１２３に記憶されているパワースペクトル分布ＰＳＤ1，ＰＳＤ2，…，ＰＳＤMを読み取り、周波数毎に、最小値であったパワースペクトルを抽出し、抽出された周波数毎のパワースペクトルに基づいて、定常ノイズ成分のパワースペクトル分布ＮＳＤを推定する。そして、ノイズ低減部１３０が、まず、最新フレーム期間における入力信号ＩＮＤのパワースペクトル分布ＰＳＤから定常ノイズ成分のパワースペクトル分布ＮＳＤを差し引く。引き続き、ノイズ低減部１３０は、当該差し引いたパワースペクトル分布と、最新フレーム期間における入力信号ＩＮＤの最新位相情報ＰＨＩとに基づいて、周波数時間変換を行い、最新フレーム期間における入力信号ＩＮＤから定常ノイズ成分を除去したノイズ低減信号ＮＳＤを生成する。
【００４３】
したがって、本実施形態によれば、周波数毎のパワーの時間的な定常性が高い定常ノイズ成分を検出、低減することができ、出力信号のＳＮ比を向上させることができる。
【００４４】
［実施形態の変形］
本発明は、上記の実施形態に限定されるものではなく、様々な変形が可能である。
【００４５】
例えば、上記の実施形態では、最新フレーム期間における入力信号のパワースペクトル分布から定常ノイズ成分のパワースペクトル分布を差し引いてノイズ低減信号を生成するというスペクトルサブトラクション法を用いてノイズ低減処理を行った。これに対して、ノイズ低減手法として、例えば、定常ノイズ成分のパワースペクトルに基づいて時間領域の定常ノイズ信号を生成し、当該定常ノイズ信号の位相を反転させた信号を、入力信号に加算して、ノイズ低減信号を生成するようにしてもよい。
【００４６】
また、上記の実施形態では、常に、入力信号のパワースペクトル分布から定常ノイズ成分のパワースペクトル分布を差し引いてノイズ低減信号を生成した。これに対して、定常ノイズ成分のスペクトル分布が、所定周波数付近にピークを有するスペクトル分布であることが予測されるときに、定常ノイズ成分のパワースペクトル分布が、当該所定周波数付近に有意なピークを有していない場合には、ノイズ低減処理を実行しないようにしてもよい。
【００４７】
なお、定常ノイズ成分のスペクトル分布が、所定周波数付近にピークを有するスペクトル分布であることが予測されるときに、定常ノイズ成分のパワースペクトル分布が、当該所定周波数付近に有意なピークを有している場合には、当該所定周波数付近のみのノイズ成分を除去するノイズ低減処理を行うようにしてもよいし、又、全周波数範囲においてノイズ成分を除去するノイズ低減処理を行うようにしてもよい。
【００４８】
また、上記の実施形態では、ラジオ受信機に本発明を適用したが、固定電話機や携帯電話機等の音声通信装置に本発明を適用してもよい。
【００４９】
なお、上記の実施形態におけるノイズ低減装置１００の一部又は全部を中央処理装置（ＣＰＵ：Central Processing Unit）、読出専用メモリ（ＲＯＭ：Read Only Memory）、ランダムアクセスメモリ（ＲＡＭ：Random Access Memory）等を備えた演算手段としてのコンピュータとして構成し、予め用意されたプログラムを当該コンピュータで実行することにより、上記の実施形態における処理の一部又は全部を実行するようにしてもよい。また、このプログラムは、ＣＤ−ＲＯＭ、ＤＶＤ等の可搬型記録媒体に記録された形態で取得されるようにしてもよいし、インターネットなどのネットワークを介した配信の形態で取得されるようにしてもよい。
【符号の説明】
【００５０】
１００  …  ノイズ低減装置
１２０  …  ノイズ検出装置
１２１  …  時間周波数変換部
１２２  …  パワースペクトル算出部
１２３  …  記憶部
１２４  …  推定部
１３０  …  ノイズ低減部
【図１】
2012177828.tif 000003
【図２】
2012177828.tif 000004
【図３】
2012177828.tif 000005
【図４】
2012177828.tif 000006
【図５】
2012177828.tif 000007
【図６】
2012177828.tif 000008