JP5092748B2 - 雑音抑圧の方法及び装置並びにコンピュータプログラム - Google Patents

Description

技術分野
［０００１］
本発明は、所望の音声信号に重畳されている雑音を抑圧するための雑音抑圧の方法及び装置並びに雑音抑圧に用いるコンピュータプログラムに関する。
背景技術
［０００２］
ノイズサプレッサ（雑音抑圧システム）は、所望の音声信号に重畳されている雑音（ノイズ）を抑圧するシステムであり、一般的に、周波数領域に変換した入力信号を用いて雑音成分のパワースペクトルを推定し、この推定パワースペクトルを入力信号から差し引くことにより、所望の音声信号に混在する雑音を抑圧するように動作する。雑音成分のパワースペクトルを継続的に推定することにより、非定常な雑音の抑圧にも適用することができる。ノイズサプレッサとしては、例えば、北米携帯電話で標準として採用されている非特許文献１（１９９６年１月、テクニカル・リクワイアメント、ＴＩＡ／ＥＩＡ／ＩＳ−１２７−１（Ｔｅｃｈｎｉｃａｌ Ｒｅｑｕｉｒｅｍｅｎｔｓ（ＴＲ４５）．ＥＮＨＡＮＣＥＤ ＶＡＲＩＡＢＬＥ ＲＡＴＥ ＣＯＤＥＣ，ＳＰＥＥＣＨ ＳＥＲＶＩＣＥ ＯＰＴＩＯＮ ３ ＦＯＲＷＩＤＥＢＡＮＤ ＳＰＲＥＡＤ ＳＰＥＣＴＲＵＭ ＤＩＧＩＴＡＬ ＳＹＳＴＥＭＳ，ＴＩＡ／ＥＩＡ／ＩＳ−１２７−１，ＳＥＰ，１９９６））及び特許文献１（特開２００２−２０４１７５号公報）に記載されている方式がある。
［０００３］
通常、音波を収集するマイクロフォンの出力信号をアナログ−ディジタル（ＡＤ）変換したディジタル信号が、入力信号としてノイズサプレッサに供給される。主として、マイクロフォンにおける集音やＡＤ変換の際に付加される低周波成分を抑圧する目的で、一般的に、高域通過フィルタがＡＤ変換とノイズサプレッサの間に配置される。このような構成の例は、例えば特許文献２（米国特許５，６５９，６２２号）に開示されている。
［０００４］
図１に、特許文献１のノイズサプレッサに特許文献２の高域通過フィルタを組み合せた構造を示す。
［０００５］
入力端子１１には、劣化音声信号（所望音声信号と雑音の混在する信号）が、サンプル値系列として供給される。劣化音声信号サンプルは、高域通過フィルタ１７に供給され、低域成分を抑圧された後、フレーム分割部１に供給される。低域成分の抑圧は、入力される劣化音声の線形性を保ち、十分な信号処理性能を発揮するためには、実用上不可欠である。フレーム分割部1は、劣化音声信号サンプルを特定の数を単位としたフレームに分割し、窓掛け処理部2へ伝達する。窓掛け処理部2は、フレームに分割された劣化音声サンプルと窓関数を乗算し、その結果をフーリエ変換部3へ伝達する。
【０００６】
フーリエ変換部3は、窓掛けされた劣化音声サンプルにフーリエ変換を施して複数の周波数成分に分割し、振幅値を多重化して、推定雑音計算部52、雑音抑圧係数生成部82、及び多重乗算部16へ供給する。位相は、逆フーリエ変換部9に伝達する。推定雑音計算部52は、供給された複数の周波数成分それぞれに対して雑音を推定し、雑音抑圧係数生成部82へ伝達する。雑音推定の方式の一例としては、過去の信号対雑音比で劣化音声を重み付けて雑音成分とする方式があり、その詳細は特許文献１に記載されている。
【０００７】
雑音抑圧係数生成部82では、推定した雑音を劣化音声に乗算することにより、雑音が抑圧された強調音声を求めるための雑音抑圧係数を複数の周波数成分それぞれに対して生成する。雑音抑圧係数生成の一例としては、強調音声の平均二乗パワーを最小化する最小平均二乗短時間スペクトル振幅法が広く用いられており、その詳細は特許文献１に記載されている。
【０００８】
周波数別に生成した雑音抑圧係数は多重乗算部16に供給される。多重乗算部16は、フーリエ変換部3から供給された劣化音声と雑音抑圧係数生成部82から供給された雑音抑圧係数を、周波数毎に乗算し、その積を強調音声の振幅として逆フーリエ変換部9に伝達する。逆フーリエ変換部9は、多重乗算部16から供給された強調音声振幅とフーリエ変換部3から供給された劣化音声の位相を合わせて逆フーリエ変換を行い、強調音声信号サンプルとしてフレーム合成部10に供給する。フレーム合成部10では、隣接フレームの強調音声サンプルを用いて当該フレームの出力音声サンプルを合成して出力端子12に供給する。
【発明の開示】
【０００９】
高域通過フィルタ１７は、直流近傍の周波数成分を抑圧するものであり、通常、100Hzから120Hzの周波数以上の成分は抑圧させずにそのまま通過させる。高域通過フィルタ１７の構成は、有限インパルス応答（ＦＩＲ）型または無限インパルス応答（ＩＩＲ）型のフィルタとすることができるが、鋭い通過帯域端特性が必要であるために、通常は後者を用いる。ＩＩＲ型フィルタは、その伝達関数が有理関数で表され、分母係数の感度が極めて高いことが知られている。従って、高域通過フィルタ１７を有限語長演算で実現する際には、十分な精度を達成するために、倍精度演算を多用しなければならず、演算量が多くなるという問題があった。一方、演算量低減のために高域通過フィルタ１７を除去すると、入力信号の線形性を保つことが困難となり、高品質な雑音抑圧が不可能になる。
［００１０］
本発明の目的は、少ない演算量で低域成分を抑圧し、高品質な雑音抑圧を達成することのできる雑音抑圧の方法及び装置を提供することである。
［００１１］
本発明に係る雑音抑圧方法は、入力信号を周波数領域信号に変換し、該周波数領域信号の振幅を補正して振幅補正信号を求め、該振幅補正信号を用いて推定雑音を求め、該推定雑音と前記振幅補正信号を用いて抑圧係数を定め、該抑圧係数で前記振幅補正信号を重みづけしている。
［００１２］
一方、本発明に係る雑音抑圧装置は、入力信号を周波数領域信号に変換する変換部と、該周波数領域信号の振幅を補正して振幅補正信号を求める振幅補正部と、該振幅補正信号を用いて推定雑音を求める雑音推定部と、該推定雑音と前記振幅補正信号を用いて抑圧係数を定める抑圧係数生成部と、該抑圧係数で前記振幅補正信号を重みづけする乗算部とを備えている。
［００１３］
更に、本発明に係る雑音抑圧の信号処理を行なうコンピュータプログラムは、前記入力信号を周波数領域信号に変換する処理と、該周波数領域信号の振幅を補正して振幅補正信号を求める処理と、該振幅補正信号を用いて推定雑音を求める処理と、該推定雑音と前記振幅補正信号を用いて抑圧係数を定める処理と、該抑圧係数で前記振幅補正信号を重みづけする処理とを有している。
［００１４］
特に、本発明に係る雑音抑圧の方法及び装置は、低域成分の抑圧をフーリエ変換後の信号に対して実行することを特徴とする。より具体的には、フーリエ変換出力の振幅に対して低域成分を抑圧するための振幅補正部と、フーリエ変換出力の位相に対して低域成分の振幅変形に対応した位相補正を行う位相補正部とを備えていることを特徴とする。
【００１５】
本発明によれば、周波数領域に変換された信号の振幅に定数を乗算し、位相に定数を加算するので、単精度演算による実現が可能となり、少ない演算量で高品質な雑音抑圧を達成することができる。
【図面の簡単な説明】
【００１６】
【図１】従来の雑音抑圧装置の構成例を示すブロック図である。
【図２】本発明の第1の実施の形態を示すブロック図である。
【図３】本発明の第1の実施の形態に含まれる振幅補正部の構成を示すブロック図である。
【図４】図３に含まれる音声存在確率計算部の構成を示すブロック図である。
【図５】本発明の第２の実施の形態を示すブロック図である。
【図６】本発明の第３の実施の形態を示すブロック図である。
【図７】本発明の第３の実施の形態に含まれる多重乗算部の構成を示すブロック図である。
【図８】本発明の第３の実施の形態に含まれる重みつき劣化音声計算部の構成を示すブロック図である。
【図９】図８に含まれる周波数別SNR計算部の構成を示すブロック図である。
【図１０】図８に含まれる多重非線形処理部の構成を示すブロック図である。
【図１１】非線形処理部における非線形関数の一例を示す図である。
【図１２】本発明の第３の実施の形態に含まれる推定雑音計算部の構成を示すブロック図である。
【図１３】図１２に含まれる周波数別推定雑音計算部の構成を示すブロック図である。
【図１４】図１３に含まれる更新判定部の構成を示すブロック図である。
【図１５】本発明の第３の実施の形態に含まれる推定先天的SNR計算部の構成を示すブロック図である。
【図１６】図１５に含まれる多重値域限定処理部の構成を示すブロック図である。
【図１７】図１５に含まれる多重重みつき加算部の構成を示すブロック図である。
【図１８】図１７に含まれる重みつき加算部の構成を示すブロック図である。
【図１９】本発明の第３の実施の形態に含まれる雑音抑圧係数生成部の構成を示すブロック図である。
【図２０】本発明の第３の実施の形態に含まれる抑圧係数補正部の構成を示すブロック図である。
【図２１】図２０に含まれる周波数別抑圧係数補正部の構成を示すブロック図である。
【符号の説明】
【００１７】
1 フレーム分割部
2,20 窓がけ処理部
3 フーリエ変換部
4,5049 カウンタ
5,52 推定雑音計算部
6,1402 周波数別SNR計算部
7 推定先天的ＳＮＲ計算部
8,82 雑音抑圧係数生成部
9 逆フーリエ変換部
10 フレーム合成部
11 入力端子
12 出力端子
13,16,704,705,1404 多重乗算部
14 重みつき劣化音声計算部
15 抑圧係数補正部
17 高域通過フィルタ
18 振幅補正部
19 位相補正部
21 音声非存在確率記憶部
22 オフセット除去部
501,502,1302,1303,1422,1423,1495,1502,1503,1801,1901,7013,7072,7074 分離部
503,1304,1424,1475,1504,1803,1903,7014,7075 多重化部
504₀〜504_K-1 周波数別推定雑音計算部
520 更新判定部
701 多重値域限定処理部
702 後天的SNR記憶部
703 抑圧係数記憶部
706 重み記憶部
707 多重重みつき加算部
708,5046,7092,7094 加算器
811 MMSE STSA ゲイン関数値計算部
812 一般化尤度比計算部
814 抑圧係数計算部
921 瞬時推定SNR
921₀〜921_K-1 周波数別瞬時推定SNR
922 過去の推定SNR
922₀〜922_K-1 過去の周波数別推定SNR
923 重み
924 推定先天的SNR
924₀〜924_K-1 周波数別推定先天的SNR
1301₀〜1301_K-1,1597,7091,7093 乗算器
1401,5042 推定雑音記憶部
1405 多重非線形処理部
1421₀〜1421_K-1,5048 除算部
1485₀〜1485_K-1 非線形処理部
1501₀〜1501_K-1 周波数別抑圧係数補正部
1591,7012₀〜7012_K-1 最大値選択部
1592 抑圧係数下限値記憶部
1593,5204,5206 閾値記憶部
1594,5203,5205 比較部
1595,5044 スイッチ
1596 修正値記憶部
1802₀〜1802_K-1 重み付け処理部
1902₀〜1902_K-1 位相回転部
5041 レジスタ長記憶部
5045 シフトレジスタ
5047 最小値選択部
5201 論理和計算部
5207 閾値計算部
7011 定数記憶部
7071₀〜7071_K-1 重みつき加算部
7095 定数乗算器
【発明を実施するための最良の形態】
【００１８】
図２は、本発明の第１の形態を示すブロック図である。図２の構成と従来例である図１の構成とは、高域通過フィルタ17、振幅補正部18、位相補正部19、窓がけ処理部20を除いて同一である。以下、これらの相違点を中心に詳細な動作を説明する。
【００１９】
図２では、図１の高域通過フィルタ17が削除され、その代わりに振幅補正部18と位相補正部19と窓がけ処理部20とが設けられている。振幅補正部18と位相補正部19は、高域通過フィルタの周波数応答を周波数領域に変換した信号に対して適用するために設けられている。高域通過フィルタ17の伝達関数にz=exp（j・2πf）を適用して得られるｆの関数の絶対値(振幅周波数応答)を振幅補正部18で入力信号に適用し、位相（位相周波数応答）を位相補正部19で入力信号に適用する。
【００２０】
これらの操作で、高域通過フィルタ17を入力信号に適用したときと同等の効果を得られる。すなわち、高域通過フィルタ17の伝達関数を時間領域で入力信号と畳み込む代わりに、フーリエ変換部３で周波数領域信号に変換された後に周波数応答を乗算することになる。
【００２１】
振幅補正部18の出力は推定雑音計算部52、雑音抑圧係数生成部82、及び多重乗算部16に供給される。位相補正部19の出力は逆フーリエ変換部9に伝達される。
【００２２】
これ以降の動作は、図１を用いて説明した通りである。窓がけ処理部20は、特許文献３（特開2003-131689号公報）に開示されているように、フレーム境界における断続音を抑圧するために設けられている。
【００２３】
図３に、振幅補正部18の構成例を示す。フーリエ変換部3から供給された多重化劣化音声振幅スペクトルは、分離部1801に伝達される。分離部1801は、多重化された劣化音声振幅スペクトルを各周波数成分に分解して、重み付け処理部1802₀〜1802_K-1に伝達する。重み付け処理部1802₀〜1802_K-1はそれぞれ、各周波数成分に分解された劣化音声振幅スペクトルを対応する振幅周波数応答で重み付けし、多重化部1803に伝達する。多重化部1803は、重み付け処理部1802₀〜1802_K-1から伝達された信号を多重化して補正劣化音声振幅スペクトルとして出力する。
【００２４】
図４に、位相補正部19の構成例を示す。フーリエ変換部3から供給された多重化劣化音声位相スペクトルは、分離部1901に伝達される。分離部1901は、多重化された劣化音声位相スペクトルを各周波数成分に分解してそれぞれ位相回転部1902₀〜1902_K-1に伝達する。位相回転部1902₀〜1902_K-1はそれぞれ、各周波数成分に分解された劣化音声位相スペクトルを対応する位相周波数応答に応じて回転させて多重化部1903に伝達する。多重化部1903は、位相回転部1902₀〜1902_K-1から伝達された信号を多重化して、補正劣化音声位相スペクトルとして出力する。位相補正部19の存在は、振幅補正部18ほど重要ではなく、省略することもできる。これは、位相補正部19の有無が出力信号の位相にしか影響を与えず、また、位相情報は音声の内容理解において、振幅情報よりもはるかに重要性が低いことが知られているからである。
【００２５】
図５は、本発明の第２の実施の形態を示すブロック図である。図５の構成と、第１の実施の形態である図２の構成との違いは、オフセット除去部22である。オフセット除去部22は、窓がけ処理された劣化音声に対してオフセットを除去して出力する。オフセット除去の最も簡単な方式は、フレーム毎に劣化音声の平均値を求めてオフセットとし、これを当該フレーム内の全サンプルから差し引くことである。また、フレーム毎の平均値を複数フレームに渡って平均化し、その平均値をオフセットとして差し引いてもよい。オフセット除去によって、次に続くフーリエ変換部3における変換精度が向上し、出力における強調音声の音質を改善することができる。
【００２６】
図６は、本発明の第３の実施の形態を示すブロック図である。入力端子11には、劣化音声信号(所望音声信号と雑音の混在する信号)が、サンプル値系列として供給される。劣化音声信号サンプルは、フレーム分割部1に供給されてK/2サンプル毎のフレームに分割される。ここで、Kは偶数とする。フレームに分割された劣化音声信号サンプルは、窓がけ処理部2に供給され、窓関数w(t)との乗算が行なわれる。第nフレームの入力信号yn(t) (t=0, 1, ..., K/2-1) に対するw(t)で窓がけされた信号yn(t)バーは、次式で与えられる。
【数１】

また、連続する２フレームの一部を重ね合わせ(オーバラップ)して窓がけすることも広く行なわれている。オーバラップ長としてフレーム長の50%を仮定すれば、t=0, 1, ..., K/2-1 に対して、

【数２】

で得られるyn(t)バー (t=0, 1, ..., K-1)が、窓がけ処理部2の出力となる。実数信号に対しては、左右対称窓関数が用いられる。また、窓関数は、抑圧係数を1に設定したときの入力信号と出力信号が計算誤差を除いて一致するように設計される。これは、w(t)+w(t+K/2)=1となることを意味する。
【００２７】
以後、連続する２フレームの50%をオーバラップして窓がけする場合を例として説明を続ける。w(t)としては、例えば次式に示すハニング窓を用いることができる。
【数３】

このほかにも、ハミング窓、ケイザー窓、ブラックマン窓など、様々な窓関数が知られている。窓がけされた出力yn(t)バーは、オフセット除去部22に供給されて、オフセットを除去される。オフセット除去の詳細に関しては、図５を用いて説明した通りである。
【００２８】
オフセット除去後の信号はフーリエ変換部3に供給され、劣化音声スペクトルYn(k)に変換される。劣化音声スペクトルYn(k)は位相と振幅に分離され、劣化音声位相スペクトル arg Yn(k)は、位相補正部19を経て、逆フーリエ変換部9に、劣化音声振幅スペクトル|Yn(k)|は、振幅補正部18を経て、多重乗算部13と多重乗算部16に供給される。位相補正部19と振幅補正部18の動作については、図２を用いて説明した通りである。
【００２９】
多重乗算部13は、振幅補正された劣化音声振幅スペクトルを用いて劣化音声パワースペクトルを計算し、推定雑音計算部5、周波数別SNR(信号対雑音比)計算部6、及び重みつき劣化音声計算部14に伝達する。重みつき劣化音声計算部14は、多重乗算部13から供給された劣化音声パワースペクトルを用いて重みつき劣化音声パワースペクトルを計算し、推定雑音計算部5に伝達する。
【００３０】
推定雑音計算部5は、劣化音声パワースペクトル、重みつき劣化音声パワースペクトル、及びカウンタ4から供給されるカウント値を用いて雑音のパワースペクトルを推定し、推定雑音パワースペクトルとして周波数別SNR計算部6に伝達する。周波数別SNR計算部6は、入力された劣化音声パワースペクトルと推定雑音パワースペクトルを用いて周波数別にSNRを計算し、後天的SNRとして推定先天的SNR計算部7と雑音抑圧係数生成部8に供給する。
【００３１】
推定先天的SNR計算部7は、入力された後天的SNR、及び抑圧係数補正部15から供給された補正抑圧係数を用いて先天的SNRを推定し、推定先天的SNRとして、雑音抑圧係数生成部8に伝達する。雑音抑圧係数生成部8は、入力として供給された後天的SNR、推定先天的SNR及び音声非存在確率記憶部21から供給される音声非存在確率を用いて雑音抑圧係数を生成し、抑圧係数として抑圧係数補正部15に伝達する。抑圧係数補正部15は、入力された推定先天的SNRと抑圧係数を用いて抑圧係数を補正し、補正抑圧係数Gn(k)バーとして多重乗算部16に供給する。多重乗算部16は、フーリエ変換部3から振幅補正部18を経て供給された補正劣化音声振幅スペクトルを、抑圧係数補正部15から供給された補正抑圧係数Gn(k)バーで重み付けすることによって強調音声振幅スペクトル|Xn(k)|バーを求め、逆フーリエ変換部9に伝達する。
【００３２】
|Xn(k)|バーは、次式で与えられる。
【数４】

ここで、Hn(k)は、振幅補正部18における補正利得であり、図１の高域通過フィルタの振幅周波数応答として得られる。
【００３３】
逆フーリエ変換部9は、多重乗算部16から供給された強調音声振幅スペクトル|Xn(k)|バーとフーリエ変換部3から位相補正部19を経て供給された補正劣化音声位相スペクトル arg Yn(k) + arg Hn(k)を乗算して、強調音声Xn(k)バーを求める。すなわち、
【数５】

を実行する。ここで、arg Hn(k)は、位相補正部19における補正位相であり、図１の高域通過フィルタの位相周波数応答として得られる。
【００３４】
逆フーリエ変換部9は、得られた強調音声Xn(k)バーに逆フーリエ変換を施し、１フレームがKサンプルから構成される時間領域サンプル値系列xn(t)バー(t=0, 1, ..., K-1)として窓がけ処理部20に供給する。窓がけ処理部20では、逆フーリエ変換部9から供給された時間領域サンプル値系列xn(t)バーと窓関数w(t)との乗算が行なわれる。第nフレームの入力信号xn(t)(t=0, 1, ..., K/2-1)に対するw(t)で窓がけされた信号xn(t)バーは、次式で与えられる。
【数６】

また、連続する2フレームの一部を重ね合わせ(オーバラップ)して窓がけすることも広く行なわれている。オーバラップ長としてフレーム長の50%を仮定すれば、t=0, 1, ..., K/2-1に対して、
【数７】

で得られるyn(t)バー(t=0, 1, ..., K-1)が、窓がけ処理部20の出力となり、フレーム合成部10に伝達される。
【００３５】
フレーム合成部10は、xn(t)バーの隣接する２フレームからK/2サンプルずつを取り出して重ね合わせ、
【数８】

によって、強調音声xn(t)ハットを得る。得られた強調音声xn(t)ハット (t=0, 1, ..., K-1)が、フレーム合成部10の出力として、出力端子12に伝達される。
【００３６】
図７は、図６に示した多重乗算部13の構成を示すブロック図である。多重乗算部13は、乗算器1301₀〜1301_K-1、分離部1302及び1303、多重化部1304を有する。多重化された状態で図６の振幅補正部18から供給された補正劣化音声振幅スペクトルは、分離部1302及び1303において周波数別のKサンプルに分離され、それぞれ乗算器1301₀〜1301_K-1に供給される。乗算器1301₀〜1301_K-1は、それぞれ入力された信号を２乗し、多重化部1304に伝達する。多重化部1304は、入力された信号を多重化し、劣化音声パワースペクトルとして出力する。
【００３７】
図８は重みつき劣化音声計算部14の構成を示すブロック図である。重みつき劣化音声計算部14は、推定雑音記憶部1401、周波数別SNR計算部1402、多重非線形処理部1405、及び多重乗算部1404を有する。推定雑音記憶部1401は、図６の推定雑音計算部5から供給される推定雑音パワースペクトルを記憶し、1フレーム前に記憶された推定雑音パワースペクトルを周波数別SNR計算部1402へ出力する。
【００３８】
周波数別SNR計算部1402は、推定雑音記憶部1401から供給される推定雑音パワースペクトルと図６の多重乗算部13から供給される劣化音声パワースペクトルを用いてSNRを各周波数毎に求め、多重非線形処理部1405に出力する。多重非線形処理部1405は、周波数別SNR計算部1402から供給されるSNRを用いて重み係数ベクトルを計算し、重み係数ベクトルを多重乗算部 1404に出力する。
【００３９】
多重乗算部1404は、図６の多重乗算部13から供給される劣化音声パワースペクトルと、多重非線形処理部1405から供給される重み係数ベクトルの積を周波数毎に計算し、重みつき劣化音声パワースペクトルを図６の推定雑音記憶部5に出力する。多重乗算部1404の構成は、既に図７を用いて説明した多重乗算部13に等しいので、詳細な説明は省略する。
【００４０】
図９は、図８に含まれる周波数別SNR計算部1402の構成を示すブロック図である。周波数別SNR計算部1402は、除算部1421₀〜1421_K-1、分離部1422及び1423、多重化部1424を有する。図６の多重乗算部13から供給される劣化音声パワースペクトルは、分離部1422に伝達される。図８の推定雑音記憶部1401から供給される推定雑音パワースペクトルは、分離部1423に伝達される。劣化音声パワースペクトルは分離部1422において、推定雑音パワースペクトルは分離部1423において、それぞれ周波数成分に対応したKサンプルに分離され、それぞれ除算部1421₀〜1421_K-1に供給される。
【００４１】
除算部1421₀〜1421_K-1では、次式に従って、供給された劣化音声パワースペクトルを推定雑音パワースペクトルで除算して周波数別SNR γn(k)ハットを求め、多重化部1424に伝達する。
【数９】

ここで、λn-1(k)は1フレーム前に記憶された推定雑音パワースペクトルである。多重化部1424は、伝達されたK個の周波数別SNRを多重化して、図８の多重非線形処理部1405へ伝達する。
【００４２】
次に、図１０を参照しながら、図８の多重非線形処理部1405の構成と動作について詳しく説明する。図１０は、重みつき劣化音声計算部14に含まれる多重非線形処理部1405の構成を示すブロック図である。多重非線形処理部1405は、分離部1495、非線形処理部1485₀〜1485_K-1、及び多重化部1475を有する。分離部1495は、図８の周波数別SNR計算部1402から供給されるSNRを周波数別のSNRに分離し、非線形処理部1485₀〜1485_K-1に出力する。非線形処理部1485₀〜1485_K-1は、それぞれ入力値に応じた実数値を出力する非線形関数を有する。
【００４３】
図１１に、非線形関数の例を示す。f1を入力値としたとき、図１１に示される非線形関数の出力値f2は、
【数１０】

で与えられる。但し、aとbは任意の実数である。
【００４４】
図１０に戻って、非線形処理部1485₀〜1485_K-1は、分離部1495から供給される周波数別SNRを非線形関数によって処理して重み係数を求め、多重化部1475に出力する。すなわち、非線形処理部1485₀〜1485_K-1はSNRに応じた１から０までの重み係数を出力する。SNRが小さい時は１を、大きい時は０を出力する。多重化部1475は、非線形処理部1485₀〜1485_K-1から出力された重み係数を多重化し、重み係数ベクトルとして多重乗算部1404に出力する。
【００４５】
図８の多重乗算部1404で劣化音声パワースペクトルと乗算される重み係数は、SNRに応じた値になっており、SNRが大きい程、すなわち劣化音声に含まれる音声成分が大きい程、重み係数の値は小さくなる。推定雑音の更新には一般に劣化音声パワースペクトルが用いられるが、推定雑音の更新に用いる劣化音声パワースペクトルに対して、SNRに応じた重みづけを行うことで、劣化音声パワースペクトルに含まれる音声成分の影響を小さくすることができ、より精度の高い雑音推定を行うことができる。なお、重み係数の計算に非線形関数を用いた例を示したが、非線形関数以外にも線形関数や高次多項式など、他の形で表されるSNRの関数を用いる事も可能である。
【００４６】
図１２は、図６に示した推定雑音計算部5の構成を示すブロック図である。雑音推定計算部5は、分離部501、502、多重化部503、及び周波数別推定雑音計算部504₀〜504_K-1を有する。
【００４７】
図１２において、分離部501は、図６の重みつき劣化音声計算部14から供給される重みつき劣化音声パワースペクトルを周波数別の重みつき劣化音声パワースペクトルに分離し、周波数別推定雑音計算部504₀〜504_K-1にそれぞれ供給する。分離部502は、図６の多重乗算部13から供給される劣化音声パワースペクトルを周波数別の劣化音声パワースペクトルに分離し、周波数別推定雑音計算部504₀〜504_K-1にそれぞれ出力する。
【００４８】
周波数別推定雑音計算部504₀〜504_K-1は、分離部501から供給される周波数別重みつき劣化音声パワースペクトル、分離部502から供給される周波数別劣化音声パワースペクトル、及び図６のカウンタ4から供給されるカウント値から周波数別推定雑音パワースペクトルを計算し、多重化部503へ出力する。多重化部503は、周波数別推定雑音計算部504₀〜504_K-1から供給される周波数別推定雑音パワースペクトルを多重化し、推定雑音パワースペクトルを図６の周波数別SNR計算部6と重みつき劣化音声計算部14へ出力する。周波数別推定雑音計算部504₀〜504_K-1の構成と動作の詳細な説明は、図１３を参照しながら行う。
【００４９】
図１３は、図１２に示した周波数別推定雑音計算部504₀〜504_K-1の構成を示すブロック図である。周波数別推定雑音計算部504は、更新判定部520、レジスタ長記憶部5041、推定雑音記憶部5042、スイッチ5044、シフトレジスタ5045、加算器5046、最小値選択部5047、除算部5048、カウンタ5049を有する。
【００５０】
スイッチ5044には、図１２の分離部501から、周波数別重みつき劣化音声パワースペクトルが供給されている。スイッチ5044が回路を閉じたときに、周波数別重みつき劣化音声パワースペクトルは、シフトレジスタ5045に伝達される。シフトレジスタ5045は、更新判定部520から供給される制御信号に応じて、内部レジスタの記憶値を隣接レジスタにシフトする。シフトレジスタ長は、後述するレジスタ長記憶部5041に記憶されている値に等しい。シフトレジスタ5045の全レジスタ出力は、加算器5046に供給される。加算器5046は、供給された全レジスタ出力を加算して、加算結果を除算部5048に伝達する。
【００５１】
一方、更新判定部520には、カウント値、周波数別劣化音声パワースペクトル及び周波数別推定雑音パワースペクトルが供給されている。更新判定部520は、カウント値が予め設定された値に到達するまでは常に“１”を、到達した後は入力された劣化音声信号が雑音であると判定されたときに“１”を、それ以外のときに“０”を出力する。更新判定部520の出力は、カウンタ5049、スイッチ5044、及びシフトレジスタ5045に伝達される。
【００５２】
スイッチ5044は、更新判定部520から供給された信号が“１”のときに回路を閉じ、“０”のときに開く。カウンタ5049は、更新判定部520から供給された信号が“１”のときにカウント値を増加し、“０”のときには変更しない。シフトレジスタ5045は、更新判定部520から供給された信号が“１”のときにスイッチ5044から供給される信号サンプルを1サンプル取り込むと同時に、内部レジスタの記憶値を隣接レジスタにシフトする。最小値選択部5047には、カウンタ5049の出力とレジスタ長記憶部5041の出力が供給されている。
【００５３】
最小値選択部5047は、供給されたカウント値とレジスタ長のうち、小さい方を選択して、除算部5048に伝達する。除算部5048は、加算器5046から供給された周波数別劣化音声パワースペクトルの加算値をカウント値又はレジスタ長の小さい方の値で除算し、商を周波数別推定雑音パワースペクトルλn(k)として出力する。Bn(k)(n=0,1, ..., N-1)をシフトレジスタ5045に保存されている劣化音声パワースペクトルのサンプル値とすると、λn(k)は、
【数１１】

で与えられる。
【００５４】
ただし、Ｎはカウント値とレジスタ長のうち、小さい方の値である。カウント値はゼロから始まって単調に増加するので、最初はカウント値で除算が行なわれ、後にはレジスタ長で除算が行なわれる。レジスタ長で除算が行なわれることは、シフトレジスタに格納された値の平均値を求めることになる。最初は、シフトレジスタ5045に十分多くの値が記憶されていないために、実際に値が記憶されているレジスタの数で除算する。実際に値が記憶されているレジスタの数は、カウント値がレジスタ長より小さいときはカウント値に等しく、カウント値がレジスタ長より大きくなるとレジスタ長と等しくなる。
【００５５】
図１４は、図１３に示した更新判定部520の構成を示すブロック図である。更新判定部520は、論理和計算部5201、比較部5203、5205、閾値記憶部5204、5206、閾値計算部5207を有する。
【００５６】
図６のカウンタ4から供給されるカウント値は、比較部5203に伝達される。閾値記憶部5204の出力である閾値も、比較部5203に伝達される。比較部5203は、供給されたカウント値と閾値を比較し、カウント値が閾値より小さいときに“１”を、カウント値が閾値より大きいときに“０”を、論理和計算部5201に伝達する。一方、閾値計算5207は、図１３の推定雑音記憶部5042から供給される周波数別推定雑音パワースペクトルに応じた値を計算し、閾値として閾値記憶部5206に出力する。最も簡単な閾値の計算方法は、周波数別推定雑音パワースペクトルの定数倍である。その他に、高次多項式や非線形関数を用いて閾値を計算することも可能である。
【００５７】
閾値記憶部5206は、閾値計算部5207から出力された閾値を記憶し、1フレーム前に記憶された閾値を比較部5205へ出力する。比較部5205は、閾値記憶部5206から供給される閾値と図１２の分離部502から供給される周波数別劣化音声パワースペクトルを比較し、周波数別劣化音声パワースペクトルが閾値よりも小さければ“１”を、大きければ“０”を論理和計算部5201に出力する。すなわち、推定雑音パワースペクトルの大きさをもとに、劣化音声信号が雑音であるか否かを判別している。論理和計算部5201は、比較部5203の出力値と比較部5205の出力値との論理和を計算し、計算結果を図１３のスイッチ5044、シフトレジスタ5045及びカウンタ5049に出力する。
【００５８】
このように、初期状態や無音区間だけでなく、有音区間でも劣化音声パワーが小さい場合には、更新判定部520は“１”を出力する。すなわち、推定雑音の更新が行われる。閾値の計算は各周波数毎に行われるため、周波数毎に推定雑音の更新を行うことができる。
【００５９】
図１５は、図６に示した推定先天的SNR計算部7の構成を示すブロック図である。推定先天的SNR計算部7は、多重値域限定処理部701、後天的SNR記憶部702、抑圧係数記憶部703、多重乗算部704、705、重み記憶部706、多重重みつき加算部707、加算器708を有する。
【００６０】
図６の周波数別SNR計算部6から供給される後天的SNR γn(k)(k=0,1, ..., K-1)は、後天的SNR記憶部702と加算器708に伝達される。後天的SNR記憶部702は、第nフレームにおける後天的SNR γn(k)を記憶すると共に、第n-1フレームにおける後天的SNR γn-1(k)を多重乗算部705に伝達する。図６の抑圧係数補正部15から供給される補正抑圧係数Gn(k)バー(k=0, 1, ..., K-1)は、抑圧係数記憶部703に伝達される。抑圧係数記憶部703は、第nフレームにおける補正抑圧係数Gn(k)バーを記憶すると共に、第n-1フレームにおける補正抑圧係数Gn-1(k)バーを多重乗算部704に伝達する。
【００６１】
多重乗算部704は、供給されたGn(k)バーを2乗してG2n-1(k)バーを求め、多重乗算部705に伝達する。多重乗算部705は、G2n-1(k)バーとγn-1(k)をk=0, 1,..., K-1に対して乗算してG2n-1(k)バーγn-1(k)を求め、結果を多重重みつき加算部707に過去の推定SNR 922として伝達する。多重乗算部704及び705の構成は、既に図７を用いて説明した多重乗算部13に等しいので、詳細な説明は省略する。
【００６２】
加算器708の他方の端子には-1が供給されており、加算結果γn(k)-1が多重値域限定処理部701に伝達される。多重値域限定処理部701は、加算器708から供給された加算結果γn(k)-1に値域限定演算子P[・]による演算を施し、結果であるP[γn(k)-1]を多重重みつき加算部707に瞬時推定SNR 921として伝達する。ただし、P[x]は次式で定められる。
【数１２】

多重重みつき加算部707には、また、重み記憶部706から重み923が供給されている。多重重みつき加算部707は、これらの供給された瞬時推定SNR 921、過去の推定SNR 922、重み923を用いて推定先天的SNR 924を求める。重み923をαとし、ξn(k)ハットを推定先天的SNRとすると、ξn(k)ハットは、次式によって計算される。
【数１３】

ここで、G2-1(k)γ-1(k)バー=1とする。
【００６３】
図１６は、図１５に示した多重値域限定処理部701の構成を示すブロック図である。多重値域限定処理部701は、定数記憶部7011、最大値選択部7012₀〜7012_K-1、分離部7013、多重化部7014を有する。分離部7013には、図１５の加算器708から、γn(k)-1が供給される。分離部7013は、供給されたγn(k)-1をK個の周波数別成分に分離し、最大値選択部7012₀〜7012_K-1に供給する。最大値選択部7012₀〜7012_K-1の他方の入力には、定数記憶部7011からゼロが供給されている。最大値選択部7012₀〜7012_K-1は、γn(k)-1をゼロと比較し、大きい方の値を多重化部7014へ伝達する。この最大値選択演算は、上述の式１２を実行することに相当する。多重化部7014は、これらの値を多重化して出力する。
【００６４】
図１７は、図１５に示した多重重みつき加算部707の構成を示すブロック図である。多重重みつき加算部707は、重みつき加算部7071₀〜7071_K-1、分離部7072、7074、多重化部7075を有する。分離部7072には、図１５の多重値域限定処理部701から、P[γn(k)-1]が瞬時推定SNR 921として供給される。分離部7072は、P[γn(k)-1]をK個の周波数別成分に分離し、周波数別瞬時推定SNR 921₀〜921_K-1として、重みつき加算部7071₀〜7071_K-1に伝達する。分離部7074には、図１５の多重乗算部705から、G2n-1(k)バーγn-1(k)が過去の推定SNR 922として供給される。
【００６５】
分離部7074は、G2n-1(k)バーγn-1(k)をK個の周波数別成分に分離し、過去の周波数別推定SNR 922₀〜922_K-1として、重みつき加算部7071₀〜7071_K-1に伝達する。一方、重みつき加算部7071₀〜7071_K-1には、重み923も供給される。重みつき加算部7071₀〜7071_K-1は、上述の式１３によって表される重みつき加算を実行し、周波数別推定先天的SNR 924₀〜924_K-1を多重化部7075に伝達する。多重化部7075は、周波数別推定先天的SNR 924₀〜924_K-1を多重化し、推定先天的SNR 924として出力する。重みつき加算部7071₀〜7071_K-1の動作と構成については、次に図１８を参照しながら説明する。
【００６６】
図１８は、図１７に示した重みつき加算部7071の構成を示すブロック図である。重みつき加算部7071は、乗算器7091、7093、定数乗算器7095、加算器7092、7094を有する。図１６の分離部7072から周波数別瞬時推定SNR 921が、図１７の分離部7074から過去の周波数別SNR 922が、図１５の重み記憶部706から重み923が、それぞれ入力として供給される。値αを有する重み923は、定数乗算器7095と乗算器7093に伝達される。定数乗算器7095は入力信号を-1倍して得られた-αを、加算器7094に伝達する。
【００６７】
加算器7094のもう一方の入力としては１が供給されており、加算器7094の出力は両者の和である１-αとなる。１-αは乗算器7091に供給されて、もう一方の入力である周波数別瞬時推定SNR P[γn(k)-1]と乗算され、積である(１-α)P[γn(k)-1]が加算器7092に伝達される。一方、乗算器7093では、重み923として供給されたαと過去の推定SNR 922が乗算され、積であるαG2n-1(k)バーγn-1(k)が加算器7092に伝達される。加算器7092は、(１-α)P[γn(k)-1]とαG2n-1(k)バーγn-1(k)の和を、周波数別推定先天的SNR 904として、出力する。
【００６８】
図１９は、図６に示した雑音抑圧係数生成部8の構成を示すブロック図である。雑音抑圧係数生成部8は、MMSE STSA ゲイン関数値計算部811、一般化尤度比計算部812、及び抑圧係数計算部814を有する。以下、非特許文献２（1984年12月、アイ・イー・イー・イー・トランザクションズ・オン・アクースティクス・スピーチ・アンド・シグナル・プロセシング、第32巻、第6号(IEEETRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING,VOL.32, NO.6,PP.1109-1121, DEC, 1984)、1109〜1121ページ）に記載されている計算式をもとに、抑圧係数の計算方法を説明する。
【００６９】
フレーム番号をn、周波数番号をkとし、γn(k)を図６の周波数別SNR計算部6から供給される周波数別後天的SNR、ξn(k)ハットを図６の推定先天的SNR計算部7から供給される周波数別推定先天的SNR、qを図６の音声非存在確率記憶部21から供給される音声非存在確率とする。また、
ηn(k)= ξn(k)ハット/(1-q)、
vn(k) = (ηn(k)γn(k))/(1+ηn(k))
とする。MMSE STSA ゲイン関数値計算部811は、図６の周波数別SNR計算部6から供給される後天的SNR γn(k)、図６の推定先天的SNR計算部7から供給される推定先天的SNR ξn(k)ハット、及び図６の音声非存在確率記憶部21から供給される音声非存在確率qをもとに、各周波数毎にMMSE STSAゲイン関数値を計算し、抑圧係数計算部814に出力する。
【００７０】
各周波数毎のMMSE STSAゲイン関数値Gn(k)は、
【数１４】

で与えられる。ここで、I0(z)は０次変形ベッセル関数、I1(z)は１次変形ベッセル関数である。変形ベッセル関数については、非特許文献３（1985年、数学辞典、岩波書店、374.Gページ）に記載されている。
【００７１】
一般化尤度比計算部812は、図６の周波数別SNR計算部6から供給される後天的SNR γn(k)、図６の推定先天的SNR計算部7から供給される推定先天的SNR ξn(k)ハット、及び図６の音声非存在確率記憶部21から供給される音声非存在確率qをもとに、周波数毎に一般化尤度比を計算し、抑圧係数計算部814に出力する。
【００７２】
周波数毎の一般化尤度比Λn(k)は、
【数１５】

で与えられる。
【００７３】
抑圧係数計算部814は、MMSE STSA ゲイン関数値計算部811から供給される MMSE STSAゲイン関数値Gn(k)と一般化尤度比計算部812から供給される一般化尤度比Λn(k)から周波数毎に抑圧係数を計算し、図６の抑圧係数補正部15へ出力する。周波数毎の抑圧係数Gn(k)バーは、
【数１６】

で与えられる。周波数別にSNRを計算する代わりに、複数の周波数から構成される帯域に共通なSNRを求めて、これを用いることも可能である。
【００７４】
図２０は、図６に示した抑圧係数補正部15の構成を示すブロック図である。抑圧係数補正部15は、周波数別抑圧係数補正部1501₀〜1501_K-1、分離部1502、1503及び多重化部1504を有する。
【００７５】
分離部1502は、図６の推定先天的SNR計算部7から供給される推定先天的SNRを周波数別成分に分離し、それぞれ周波数別抑圧係数補正部1501₀〜1501_K-1に出力する。分離部1503は、図６の抑圧係数生成部8から供給される抑圧係数を周波数別成分に分離し、それぞれ周波数別抑圧係数補正部1501₀〜1501_K-1に出力する。
【００７６】
周波数別抑圧係数補正部1501₀〜1501_K-1は、分離部1502から供給される周波数別推定先天的SNRと、分離部1503から供給される周波数別抑圧係数から、周波数別補正抑圧係数を計算し、多重化部1504へ出力する。多重化部1504は、周波数別抑圧係数補正部1501₀〜1501_K-1から供給される周波数別補正抑圧係数を多重化し、補正抑圧係数として図６の多重乗算部16と推定先天的SNR計算部7へ出力する。
【００７７】
次に、図２１を参照しながら、周波数別抑圧係数補正部1501₀〜1501_K-1の構成と動作について詳細に説明する。
【００７８】
図２１は、抑圧係数補正部15に含まれる周波数別抑圧係数補正部1501₀〜1501_K-1の構成を示すブロック図である。周波数別抑圧係数補正部1501は、最大値選択部1591、抑圧係数下限値記憶部1592、閾値記憶部1593、比較部1594、スイッチ1595、修正値記憶部1596及び乗算器1597を有する。
【００７９】
比較部1594は、閾値記憶部1593から供給される閾値と、図２０の分離部1502から供給される周波数別推定先天的SNRを比較し、周波数別推定先天的SNRが閾値よりも大きければ“０”を、小さければ“1”をスイッチ1595に供給する。スイッチ1595は、図２０の分離部1503から供給される周波数別抑圧係数を、比較部1594の出力値が“１”のときに乗算器1597に出力し、“０”のときに最大値選択部1591に出力する。すなわち、周波数別推定先天的SNRが閾値よりも小さいときに、抑圧係数の補正が行われる。乗算器1597は、スイッチ1595の出力値と修正値記憶部1596の出力値との積を計算し、最大値選択部1591に出力する。
【００８０】
一方、抑圧係数下限値記憶部1592は、記憶している抑圧係数の下限値を、最大値選択部1591に供給する。最大値選択部1591は、図２０の分離部1503から供給される周波数別抑圧係数、又は乗算器1597で計算された積と、抑圧係数下限値記憶部1592から供給される抑圧係数下限値とを比較し、大きい方の値を図２０の多重化部1504に出力する。すなわち、抑圧係数は抑圧係数下限値記憶部1592が記憶する下限値よりも必ず大きい値になる。
【００８１】
これまで説明した全ての実施の形態では、雑音抑圧の方式として、最小平均２乗誤差短時間スペクトル振幅法を仮定してきたが、その他の方法にも適用することができる。このような方法の例として、非特許文献４（1979年12月、プロシーディングス・オブ・ザ・アイ・イー・イー・イー、第67巻、第12号(PROCEEDINGSOF THE IEEE, VOL.67, NO.12, PP.1586-1604, DEC, 1979)、1586〜1604ページ）に開示されているウィーナーフィルタ法や、非特許文献５（1979年4月、アイ・イー・イー・イー・トランザクションズ・オン・アコースティクス・スピーチ・アンド・シグナル・プロセシング、第27巻、第2号(IEEETRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING,VOL.27, NO.2,PP.113-120, APR, 1979)、113〜120ページ）に開示されているスペクトル減算法などがあるが、これらの詳細な構成例については説明を省略する。
【００８２】
また、上述した各実施形態の雑音抑圧装置は、プログラムなどを蓄積する記憶装置、入力用のキーやスイッチが配置された操作部、ＬＣＤなどの表示装置、操作部からの入力を受け付けて各部の動作を制御する制御装置から構成されるコンピュータ装置によって構成することができる。前述した各実施形態の雑音抑圧装置における動作は、制御装置が記憶装置に格納されたプログラムを実行することで実現される。プログラムは、予め記憶部に格納されていてもよく、また、CD-ROMなどの記録媒体に書き込まれた状態でユーザに提供されてもよい。また、ネットワークを通じて、プログラムを提供することも可能である。

Claims (14)

1. 入力信号に含まれている雑音を抑圧する方法であって、
   入力信号のオフセットを除去してオフセット除去信号を求め、
   該オフセット除去信号を周波数領域信号に変換し、
   該周波数領域信号の振幅を前記オフセット除去の処理と合わせて所望の高域通過特性となるように補正して振幅補正信号を求め、
   該振幅補正信号を用いて推定雑音を求め、
   該推定雑音と前記振幅補正信号を用いて抑圧係数を求め、
   該抑圧係数で前記振幅補正信号を重みづけする、
   ことを特徴とする雑音抑圧の方法。
2. 前記所望の高域通過特性は、直流近傍成分を抑圧し、音声は通過させる、
   ことを特徴とする請求項１に記載の雑音抑圧の方法。
3. 前記周波数領域信号の位相を補正して位相補正信号を求め、
   前記抑圧係数で前記振幅補正信号を重み付けした結果と前記位相補正信号を時間領域信号に変換する、
   ことを特徴とする請求項１または２に記載の雑音抑圧の方法。
4. 前記周波数領域信号の振幅の補正は、前記所望の高域通過特性を得ることが可能な伝達関数に所定の複素指数関数を適用して得られる周波数の関数の絶対値を、前記周波数領域信号の振幅に適用することで、前記オフセット除去の処理と合わせて所望の高域通過特性となるように行なわれる、
   ことを特徴とする請求項２または３に記載の雑音抑圧の方法。
5. 入力信号に含まれている雑音を抑圧する装置であって、
   入力信号のオフセットを除去してオフセット除去信号を求めるオフセット除去部と、
   該オフセット除去信号を周波数領域信号に変換する変換部と、
   該周波数領域信号の振幅を前記オフセット除去の処理と合わせて所望の高域通過特性となるように補正して振幅補正信号を求める振幅補正部と、
   該振幅補正信号を用いて推定雑音を求める雑音推定部と、
   該推定雑音と前記振幅補正信号を用いて抑圧係数を定める抑圧係数生成部と、
   該抑圧係数で前記振幅補正信号を重みづけする乗算部と、
   を有することを特徴とする雑音抑圧の装置。
6. 前記振幅補正部は、前記周波数領域信号の振幅を、前記オフセット除去の処理と合わせて、直流近傍成分を抑圧し、音声は通過させるように補正する、
   ことを特徴とする請求項５に記載の雑音抑圧の装置。
7. 前記周波数領域信号の位相を補正して位相補正信号を求める位相補正部と、
   前記抑圧係数で前記振幅補正信号を重みづけした結果と前記位相補正信号を時間領域信号に変換する逆変換部と、
   を有することを特徴とする請求項５または６に記載の雑音抑圧の装置。
8. 前記振幅補正部は、前記所望の高域通過特性を得ることが可能な伝達関数に所定の複素指数関数を適用して得られる周波数の関数の絶対値を、前記周波数領域信号の振幅に適用することで、前記オフセット除去の処理と合わせて所望の高域通過特性となるように前記周波数領域信号の振幅を補正する、
   ことを特徴とする請求項６または７に記載の雑音抑圧の装置。
9. 入力信号の特定の周波数成分を抑圧するフィルタ処理方法であって、
   入力信号に時間領域で第１のフィルタ処理を行って時間領域フィルタ処理信号を求め、
   該時間領域フィルタ処理信号を複数サンプルから構成されるフレーム毎に周波数領域信号に変換し、
   該周波数領域信号に周波数領域で第２のフィルタ処理を行って周波数領域フィルタ処理信号を求め、
   前記第１のフィルタ処理が少なくとも直流成分を抑圧し、
   前記第１のフィルタ処理と前記第２のフィルタ処理を合わせた特性が、直流近傍成分を抑圧し、音声は通過させる、
   ことを特徴とするフィルタ処理方法。
10. 入力信号の特定の周波数成分を抑圧するフィルタであって、
    入力信号に時間領域で第１のフィルタ処理を行って時間領域フィルタ処理信号を求める第１のフィルタと、
    該時間領域フィルタ処理信号を複数サンプルから構成されるフレーム毎に周波数領域信号に変換する変換部と、
    該周波数領域信号に周波数領域で第２のフィルタ処理を行って周波数領域フィルタ処理信号を求める第２のフィルタと、を少なくとも具備し、
    前記第一のフィルタが少なくとも直流成分を抑圧し、
    前記第１のフィルタ処理と前記第２のフィルタ処理を合わせた特性が、直流近傍成分を抑圧し、音声は通過させる、
    ことを特徴とするフィルタ。
11. 入力信号に含まれている雑音を抑圧する信号処理を行うコンピュータプログラムであって、
    前記入力信号のオフセットを除去してオフセット除去信号を求める処理と、
    該オフセット除去信号を周波数領域信号に変換する処理と、
    該周波数領域信号の振幅を前記オフセット除去の処理と合わせて所望の高域通過特性となるように補正して振幅補正信号を求める処理と、
    該振幅補正信号を用いて推定雑音を求める処理と、
    該推定雑音と前記振幅補正信号を用いて抑圧係数を求める処理と、
    該抑圧係数で前記振幅補正信号を重みづけする処理と、
    をコンピュータに実行させるコンピュータプログラム。
12. 前記振幅補正信号を求める処理は、前記周波数領域信号の振幅を、前記オフセット除去の処理と合わせて、直流近傍成分を抑圧し、音声は通過させるように補正することである、
    請求項１１に記載のコンピュータプログラム。
13. 前記周波数領域信号の位相を補正して位相補正信号を求める処理と、
    前記抑圧係数で前記振幅補正信号を重みづけした結果と前記位相補正信号を時間領域信号に変換する処理と、
    をさらにコンピュータに実行させる、請求項１１または１２に記載のコンピュータプログラム。
14. 前記周波数領域信号の振幅の補正は、前記所望の高域通過特性を得ることが可能な伝達関数に所定の複素指数関数を適用して得られる周波数の関数の絶対値を、前記周波数領域信号の振幅に適用することで、前記オフセット除去の処理と合わせて所望の高域通過特性となるように行なわれる、
    請求項１２または１３に記載のプログラム。

Images