雑音のある音声信号を処理する方法および該方法を実行する装置
本発明は、雑音のある音声信号の処理の方法と、該方法を実行する装置に関する。本発明による方法は、雑音のある音声信号の各フレーム（ｘ（ｋ，ｎ））に対して実行される次のステップ、すなわち、周波数領域変換を適用するステップと、フレームに対して雑音の電力スペクトル密度（ＰＳＤ）を推定するステップと、雑音の推定されたＰＳＤと、フレームに対応する求める信号のＰＳＤの推定値から第１の雑音低減フィルタを計算するステップと、フレームに対する第１の雑音修正された推定値を得るために、計算された雑音低減フィルタを利用してフレーム（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））を濾波するステップと、前述のフレームの第１の雑音修正された推定値から、雑音のある音声信号のフレームに対応する求める信号と本質的に同じ位置にある高調波を有する、第２の信号のフレームを得るステップと、を有する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、雑音のある音声信号の処理に関する。本発明は特に、そのような信号に存在する雑音の低減に関する。
【背景技術】
【０００２】
音声信号の中の雑音、すなわち妨害信号、を低減する技術が公知である。それらは、音声信号の品質と明瞭性を改善するために、音声信号が出現する音響的環境を考慮することを目的としている。これらの技術は、考慮している音声信号から、この雑音のある信号に処理を行うことによって、有用な情報を抽出することにある。このような技術は、音声信号が何人かの発言者の間で送信される、電話、電話会議、ビデオ会議のような用途における、たとえば会話通信に適用される。このような技術は、さらに、雑音のある環境での音声収集の用途に、あるいは、音声信号が雑音に満ちた環境で発音された場合に性能が大幅に変化してしまう音声認識の用途に適用される。
【０００３】
これらの技術は通常、雑音低減フィルタの伝達関数を推定し、次に、スペクトル領域での乗算にもとづいて、濾波・プロセスを実行することにある。これらの技術は、「短期間スペクトル減衰による雑音低減」と呼ばれる方法に含まれる。
【０００４】
これらの技術によれば、考慮している音声信号ｘ（ｎ）は、有用な信号成分ｓ（ｎ）と雑音成分ｂ（ｎ）を含み、ここでｎは離散時間における時間インデックスを表している。しかし、連続した時間での信号の表現も採用することも可能なことが気づくであろう。信号ｘ（ｎ）は、一定の長さと、インデックスｋの連続した複数のフレームｘ（ｎ、ｋ）として構成される。これらのフレームのそれぞれは、まず重みづけウインドウと乗算され、雑音低減フィルタの計算に必要なスペクトル量の、その後の推定を改善することを可能にする。このようにウインドウ化された各フレームは次に、たとえば離散的あるいは高速フーリエ変換を用いて、スペクトル領域において分析される、この動作は、短期間フーリエ変換（ＳＴＦＴ）と呼ばれる。
【０００５】
注目する信号の、このようにして得られた周波数表現Ｘ（ｋ，ｆ）（ここで、ｆは周波数インデックス）は同一の時間において、雑音低減フィルタの伝達関数Ｈ（ｋ、ｆ）を推定することと、この伝達関数と雑音のある信号の短期間スペクトルの簡単な乗算によって、このフィルタをスペクトル領域で用いることを可能にする。したがって、濾波の結果は次式
【０００６】
【数１】
【０００７】
のように、記述することができる。
【０００８】
次に、得られた信号の時間領域への戻りが、逆スペクトル変換によって実行される。対応する時間信号は、ブロック・オーバーラップ・アッド技術（ＯＬＡは「オーバーラップ・アッド」を表わす）によって、あるいは、ブロック・セーブ技術（ＯＬＳは「オーバーラップ・セーブ」を表わす）によって、最終的に合成される。信号を時間領域においで再構築するこの動作は、逆短期間フーリエ変換（ＩＳＴＦＴ）と呼ばれる。
【０００９】
短期間スペクトル減衰方法の詳細な説明が、下記の参考文献に見られる。J.S. Lim、A.V. Oppenheim、「雑音のある通話のエンハンスメントと帯域幅圧縮」、Proceedings of the IEEE、Vol. 67、pp. 1586-1604、1979、および、 R.E. Crochiere、L.R. Rabiner、「マルチレート・ディジタル信号処理」、プレンティス・ホール、１９８３。
【００１０】
周波数成分ｆを有する、インデックスｋの時間セグメントにわたる観察信号Ｘ（ｋ、ｆ）に適用される、短期間スペクトル減衰Ｈ（ｋ、ｆ）は通常、ローカルな信号対雑音比ＳＮＲ（ｋ、ｆ）の推定値にもとづいて求められる。すべての抑圧規則に共通する特性は、次式
【００１１】
【数２】
【００１２】
によって与えられる、それらの漸近的な挙動にある。
【００１３】
多くの技術において、次の仮定がなされている。すなわち、雑音と有用な信号は統計的に無相関であり、有用な信号は断続的（無音期間の存在）であり、人間の耳は信号の位相（通常、処理によって修正されない）に気づかない。
【００１４】
一般に使用される抑圧規則のうちで、電力スペクトル減算、振幅スペクトル減算、およびウイーナー・フィルタの直接実施を、例として挙げることができる。これらの規則に対して、有用な音声信号(speech signal)の周波数成分ｆの短期間の推定値は、それぞれ次のように記述することができる。
【００１５】
電力スペクトル減算に対して（ J.S. Lim および A.V. Oppenheimによる前述の論文を参照のこと）
【００１６】
【数３】
【００１７】
振幅スペクトル減算に対して（ S.F. Boll「スペクトル減算を用いた、会話における音響雑音の抑圧」、IEEE Trans. on Audio, Speech and Signal Processing、Vol. 27、No. 2、pp. 113-120、April 1979、参照）。
【００１８】
【数４】
【００１９】
ウイーナー・濾波に対して（ J.S. Lim および A.V. Oppenheimによる前述の論文を参照のこと）
【００２０】
【数５】
【００２１】
これらの式において、
【００２２】
【数６】
【００２３】
と
【００２４】
【数７】
【００２５】
は、それぞれ、インデックスｋの時間ウインドウにわたる、観察信号Ｘ（ｋ、ｆ）の周波数成分ｆ内に存在する有用な信号の電力スペクトル密度と、雑音の電力スペクトル密度を表す。
【００２６】
雑音のある信号に適用されたスペクトル減衰の挙動を、上記の式にもとづいて、所与の周波数成分ｆに対して測定したローカルな信号対雑音比の関数として、検討することが可能である。上述の短期間の３つの抑圧規則について、これらの曲線が図１にプロットされている。ローカルな信号対雑音比が大きい場合には（図１の右側の部分）、各規則は実質的に同じ減衰を与えることがわかる。ガウス・モデルに対して最大の尤度の意味で最適な、最適な電力減算規則（ 0. Cappe、「 Ephraim および Malah 雑音サプレッサーを使用するミュージカル・ノイズ現象の除去」、IEEE Trans. on Speech and Audio Processing、Vol. 2、No. 2、pp. 345-349、April 1994、参照）は、依然として、雑音の電力が処理の出力において最上位の規則である。３つの抑圧規則に対して、カットオフ値に関するローカルな信号対雑音比の僅かな変化が、完全な減衰（Ｈ（ｋ、ｆ）≒０）から、スペクトルの無視できる修正（Ｈ（ｋ、ｆ）≒１）への切り替えで十分であることを指摘できる。
【００２７】
この後者の性質は、「音楽的雑音(musical noise)」と呼ばれる現象の原因の１つである。具体的にいうと、同時に決定論的、かつランダムな成分を有する周囲騒音は、ボーカル非アクティビティ（vocal non-activity）の期間のみを特徴とすることができる。ランダム成分の存在のために、ボーカルアクティビティ（vocal activity ）の期間の雑音の周波数成分ｆの実際の寄与と、ボーカル非アクティビティの各瞬間の間のいくつかのフレームにわたって行われたその平均の推定値との間には非常に大きい違いがある。この違いのために、ローカルな信号対雑音比の推定値は、カットオフ・レベルの付近で変動することがあり、したがって、現れては消え、かつ平均の寿命が、考慮している分析ウインドウの大きさを統計的に超えないスペクトル成分を処理の出力に生ずることがある。通過帯域の全体にわたるこの挙動の一般化は、可聴で不快な残留雑音を招く。
【００２８】
この残留雑音の影響を減らすいくつかの研究が試みられてきた。推薦された解決策は、いくつか方法を取っている。すなわち、短期間の推定値の平均化（ S.F. Boll による前述の論文を参照）、雑音電力スペクトルの過大推定（ M. Berouti 他、「音響雑音によって損なわれた通話のエンハンスメント」、Int. Conf. on Speech, Signal Processing, pp. 208-211, 1979、および、P. Lockwood、 J. Boudy、「非線形スペクトル減算器、隠れマルコフ・モデルおよび射影演算を使用する自動車内の頑強な音声認識のための実験」Proc. of EUSIPCO'91, pp. 79-82, 1991、参照）、あるいは、雑音スペクトル密度の最小値の追跡（ R. Martin, 「最小の統計にもとづくスペクトル減算」 Signal Processing VII: Theories and Applications, EUSIPCO'94, pp. 1182-1185, September 1994、参照）。
【００２９】
音楽的雑音を抑制する比較的に有効な解決策は、「有向決定」（directed-decision）と呼ばれる有用な信号の電力スペクトル密度の推定値から成る。（ Y. Ephraim, and D. Malah、「最小平均２乗誤差短期間スペクトル振幅推定値を用いた通話向上」、IEEE Trans. on Audio, Speech and Signal Processing、Vol. 32、No. 6、pp. 1109-1121、1984、および O.Cappe による前述の論文を参照）。この推定値は、有用な信号の瞬間的電力スペクトル密度と長期間の電力スペクトル密度の間で妥協を行い、その結果、音楽的雑音を効果的に除去するのが可能になる。この推定値に固有の遅れを補償することによって、この解決策を改善することも公知である（フランス特許出願公開明細書第２８２０２２７号、および C. Plapous、C. Marro、L. Mauuary、P. Scalart、「２ステップ雑音低減技術」、ICASSP, May 2004、参照）。
【００３０】
いくつかの研究も、通話と付加的な雑音信号の統計モデルにもとづく新しい抑圧規則を確立することに関している。これらの研究は、従来の方法に対して付加的な自由度を有するので、「ソフト・ディシジョン」アルゴリズムと呼ばれる新しいアルゴリズムの導入を可能にした。（ R.J. Mac Aulay、M.L. Malpass、「ソフト・ディシジョン雑音抑圧フィルタを用いた通話向上」、IEEE trans. on Audio, Speech and Signal Processing、Vol. 28、No. 2、pp. 137-145、April 1980、Y. Ephraim、 D. Malah、「最適な非線形スペクトル振幅推定値を使用する通話向上」、Int. Conf. on Speech, Signal Processing、pp. 1118-1121、1983、および、Y. Ephraim、D. Malah の上述の論文、「最小平均２乗誤差短期間スペクトル振幅推定値を使用する通話向上」、参照）。
【００３１】
上述したように、短期間スペクトル減衰の計算は、スペクトル成分のそれぞれに対する信号対雑音比の推定に依存している。例として、上記した式は、それぞれ次の量を含んでいる。
【００３２】
【数８】
【００３３】
このように、特に、歪みに関して、また、雑音レベルの効果的低減に関しての、雑音低減技術の性能は、信号対雑音比のこの推定値の妥当性に左右される。
【００３４】
この欠点は、公知の通話雑音除去システムの主要な限界となっている。具体的にいうと、現在の雑音除去システムは、信号対雑音比が低すぎることを特徴とする高調波の雑音除去は不可能である。実際には、雑音除去アルゴリズムは、各周波数に通話成分が存在するか、あるいは存在しないかを検出するために、ＳＮＲを使用している。推定ＳＮＲがあまりにも好ましくないものであれば、アルゴリズムは信号成分がないとみなし、それを抑圧する。このように、公知の雑音除去システムによって高調波が損なわれることがあるが、このような高調波が存在しなければならないことは、先験的に公知である。ここで、大部分の言語において、口に出された音（voiced sounds）（高調波）は、発声された音（sounds uttered）の非常に大きい部分を表していることに、注目すべきである。
【発明の開示】
【発明が解決しようとする課題】
【００３５】
本発明の目的は、公知の雑音除去システムの限界を克服することである。
【００３６】
本発明の他の目的は、雑音除去方法の性能を改善することである。
【００３７】
本発明の他の目的は、信号を過度に歪ませない音声信号処理を提案することである。特に、実行される信号の処理は、この信号に含まれる高調波のすべてあるいは一部の保存を可能にする。
【００３８】
本発明の他の目的は、音声信号処理の完了に際して、音楽的雑音の出現を制限することである。
【００３９】
本発明の他の目的は、有用な信号の高調波コム（harmonic comb）の良好な推定値を得ることである。
【課題を解決するための手段】
【００４０】
本発明は、連続したフレームとして構成された雑音のある音声信号の処理方法を提案する。該処理は、前記フレームの少なくとも１つに関して、
雑音のある音声信号の前記フレームに、周波数領域への変換を適用するステップと、
前記フレームに対する雑音の電力スペクトル密度を推定するステップと、
雑音の推定された電力スペクトル密度と、フレームに対応する有用な信号の電力スペクトル密度の推定値とにもとづいて、第１の雑音低減フィルタを計算するステップと、
前記フレームの第１の、雑音が除去された推定値を得るために、計算された前記第１の雑音低減フィルタを利用して、雑音のある音声信号の前記フレームを濾波するステップと、
前記雑音のある音声信号の前記フレームに対応する前記有用な信号と実質的に同じ位置にある高調波を有する第２の信号のフレームを、雑音のある音声信号の前記フレームの第１の、雑音が除去された推定値にもとづいて得るステップと、
を有する。
【００４１】
したがって、雑音のある音声信号のこのような処理は、信号に適用された第１の雑音低減フィルタの出力において高調波性の再生を得ることを可能にする。第２の信号のこのようにして得られたフレームは、第１の濾波の間に出現した可能性のある、雑音が除去された信号のフレームの歪みを除去するように作られる。
【００４２】
この処理は、後での利用の主題である、有用な信号の高調波コムの良好な推定値を得ることを可能にする。
【００４３】
有利なことに、雑音除去の領域において、本方法はさらに次のステップ、すなわち、
前記雑音の推定された電力スペクトル密度と、前記フレームの第１の、雑音が除去された推定値の電力と得られた第２の信号の前記フレームの電力の組み合わせとにもとづいて、第２の雑音低減フィルタを計算するステップと、
前記フレームの第２の、雑音が除去された推定値を得るために、計算された第２の雑音低減フィルタを利用して、前記雑音のある音声信号の前記フレームを濾波するステップと、
前記フレームの第２の雑音が除去された推定値を合成するステップと
を有する。
【００４４】
この実施態様においては、第２の雑音低減フィルタは、それが、高調波再生から生ずる信号によって駆動されるので、高調波を保存するように計算される。第２の濾波の完了時に得られた前記フレームの第２の雑音が除去された推定値は、したがって、入力信号の高調波が破壊されるかあるいは少なくとも損なわれる従来の雑音除去システムで得られた推定値よりも有効である。
【００４５】
第１の雑音低減フィルタと、必要ならば、第２の雑音低減フィルタの計算は、短期間スペクトル減衰技術を実行する第１のパスを含んでもよい。実例として、たとえば、次の技術、電力スペクトル減算、振幅スペクトル減算、および開ループ・ウイーナー・フィルタ等を引用することができる。有利なことに、それは、たとえば、電力スペクトル減算、振幅スペクトル減算、および開ループ・ウイーナー・フィルタ等の短期間スペクトル減衰技術を実行する第２のパスをさらに有してもよく、また、第１のパスの間に実行される計算を考慮して、前記フレームに対応する有用な信号の電力スペクトル密度の推定値を使用してもよい。
【００４６】
有利なことに、雑音のある音声信号の前記フレームに対応する有用な信号と実質的に同じ位置にある高調波を有する第２の信号のフレームを得ることは、前記第１の推定値が時間領域内にある場合には、雑音のある音声信号の前記フレームの第１の、雑音が除去された推定値への非線形の関数の適用を含み、前記第１の推定値が周波数領域内にある場合には、雑音のある音声信号の前記フレームの第１の、雑音が除去された推定値と非線形の関数の周波数領域への変換の第１の、雑音が除去された推定値との間の円形の畳み込みの適用を含む。
【００４７】
前記非線形の関数は、たとえば、単一波の修正関数、絶対値、雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値と閾値との間の最大値、雑音のある音声信号の前記フレームの前記第１の雑音が除去された推定値と閾値との間の最小値、の中の１つであってもよい。
【００４８】
本発明は、上述の方法を実行するように設計された手段を有する、雑音のある音声信号を処理する装置をさらに提案する。
【００４９】
本発明は、プログラムが計算手段によってロードされ実行されるときに、上述の方法を実行するように作られた命令を有する、情報媒体上のコンピュータ・プログラムをさらに提案する。
【発明を実施するための最良の形態】
【００５０】
本発明による他の特徴と利点は、添付図面を参照した非限定的な代表的実施形態の下記の説明において明らかとなろう。
【００５１】
図２は、本発明による一実施形態による装置１を示す。雑音のある音声信号の現在のフレーム２が、処理を受けるために装置１に入力される。このフレームは、前述の表記法によって、Ｘ（ｋ，ｎ）と示される。
【００５２】
装置１は、役割が、雑音のある信号ｘ（ｋ、ｎ）の現在のフレームを周波数領域Ｘ（ｋ、ｆ）に切り替えることである分析ユニット３を有する。まず、フレームは、信号ｘw（
ｋ，ｎ）
【００５３】
【数９】
【００５４】
を送出する重みづけウインドウｗ（ｎ）と乗算される。
【００５５】
次に、離散的フーリエ変換（ＤＦＴ）、すなわち
【００５６】
【数１０】
【００５７】
を使用して、周波数領域への切り替えが実行される。
【００５８】
当該ＤＦＴは、高速フーリエ変換（すなわちＦＦＴ）によって実行するのが有利である。しかし、ウェーブレット変換のような、周波数領域への他の変換も可能である。同じことが、対応する逆動作、すなわち、後に説明する、時間領域に戻すための逆離散的フーリエ変換（ＩＤＦＴ）にも当てはまる。
【００５９】
さらに、ボーカルアクティビティを検出する機能４（ＤＶＡ）を、現在のフレーム２で実行するのが有利である。ＤＶＡは、雑音の電力スペクトル密度（ＰＳＤ）の推定値をいつ更新するべきかを知ることを可能にする。したがって、ＤＶＡによって検出された各「雑音のみの」フレームｋbに対して、雑音電力スペクトル密度
【００６０】
【数１１】
【００６１】
が下記の再帰式
【００６２】
【数１２】
【００６３】
にしたがって、機能ブロック５によって推定される。
【００６４】
パラメータα（ｋb）は平滑化係数を制御する。該パラメータ長時間にわたって変化することがある。
【００６５】
現在のフレーム２が雑音フレームとして検出されなければ、雑音の電力スペクトル密度の推定値は固定される、すなわち、以前の雑音フレームに対して得られたＰＳＤの最後の値が保持される。
【００６６】
【数１３】
【００６７】
の推定値は、この指数関数的平滑化推定値（exponential-smoothing estimator）に限定されず、任意の他の電力スペクトル密度推定値を使用できることを指摘しておく。
【００６８】
時間領域で動作するか、周波数領域で動作するかにかかわらず、任意の種類のＤＶＡを使用することができる。ボーカルアクティビティのこのような検出を不要にすることも可能である。
【００６９】
次に、第１の雑音低減フィルタの推定が、図２の機能ブロック７によって実行される。この第１の雑音低減フィルタの伝達関数は、２パス技術（two-pass technique）によって周波数領域で有利なことに推定される。（フランス特許出願公開明細書第２８２０２２７号、および、上述の C. Plapous、C. Marro、L. Mauuary、P. Scalart、「２ステップの雑音除去技術」、ICASSP, May 2004、を参照）。第１のパスにおいて、次の伝達関数を有するフィルタが計算される。
【００７０】
【数１４】
【００７１】
フィルタの伝達関数に対するこの式は、有用な信号
【００７２】
【数１５】
【００７３】
のＰＳＤと、雑音
【００７４】
【数１６】
【００７５】
のＰＳＤに依存する数式である。この意味で、
【００７６】
【数１７】
【００７７】
が、たとえば、電力スペクトル減算規則（式（１））、振幅スペクトル減算規則（式（２））、開ループ・ウイーナー・フィルタ規則（式（３））等を実行するように、関数ｆstep1を短期間スペクトル減衰技術にしたがって選択してもよい。周波数ベースの任意の他の雑音抑圧規則も、
【００７８】
【数１８】
【００７９】
を推定するために、実行してもよい。
【００８０】
雑音電力スペクトル密度
【００８１】
【数１９】
【００８２】
の計算を、以上に詳細に説明した。スペクトル量
【００８３】
【数２０】
【００８４】
に関しては、それは、ボーカルアクティビティ期間の間、信号と雑音の混合のために、直接得ることはできない。それを計算するためには、次式
【００８５】
【数２１】
【００８６】
にしたがって有向決定推定値（directed-decision estimator）が使用される。（上述の、Y. Ephraim、 D. Malah、「最小平均２乗誤差短期スペクトル振幅推定値を使用する通話エ向上」を参照）。ここで、β（ｋ）は、時間と共に変化することがある重心パラメータであり、
【００８７】
【数２２】
【００８８】
は、インデックスｋ−１の以前のフレームに関して推定された有用な信号のスペクトルである。推定誤差の場合には負になる可能性がある量
【００８９】
【数２３】
【００９０】
の閾値処理（thresholding）を保証する関数Ｐは、次式によって与えられる。
【００９１】
【数２４】
【００９２】
【数２５】
【００９３】
の推定は、この有向決定推定値に限定されないことを指摘しておく。具体的にいうと、指数関数的平滑化推定値あるいは任意の他の電力スペクトル密度推定値を使用してもよい。
【００９４】
伝達関数
【００９５】
【数２６】
【００９６】
はその後、有用な信号
【００９７】
【数２７】
【００９８】
のＰＳＤの推定を改良するために再使用してもよい。次に、量
【００９９】
【数２８】
【０１００】
を、次式によって得る。
【０１０１】
【数２９】
【０１０２】
次に、第２のパスは、
【０１０３】
【数３０】
【０１０４】
にもとづいて、第１の雑音低減フィルタの伝達関数の推定値
【０１０５】
【数３１】
【０１０６】
、すなわち、
【０１０７】
【数３２】
【０１０８】
を計算することである。
【０１０９】
独特の特徴が、有用な信号
【０１１０】
【数３３】
【０１１１】
のＰＳＤの「高速な」更新にある、この２パス計算は、第１の雑音低減フィルタ
【０１１２】
【数３４】
【０１１３】
に２つの利点をもたらす。一方では、有用な信号の非定常性を高速に、特にその時間的包絡線（たとえば、無音／通話遷移の間の通話信号の開始あるいは終了（attacks or extinctions））の高速な変化の間に、追跡することが可能である。他方では、雑音低減フィルタはより良好に推定され、これは、本方法の性能の向上によって明らかにされる。
【０１１４】
【数３５】
【０１１５】
および
【０１１６】
【数３６】
【０１１７】
は、たとえば、電力スペクトル減算規則（式（１））、振幅スペクトル減算規則（式（２））、あるいは、開ループ・ウイーナー・フィルタ規則（式（３））のような、短期間スペクトル減衰技術を実行してもよい。周波数ベースの何らかの他の雑音抑圧規則も、
【０１１８】
【数３７】
【０１１９】
および
【０１２０】
【数３８】
【０１２１】
を推定するために、実行してもよい。勿論、第２のバスを実行しないで、第１のパスのみに制限することも可能である。
【０１２２】
次に、第１の計算されたフィルタ
【０１２３】
【数３９】
【０１２４】
にしたがって、装置１の機能ブロック６によって濾波が実行される。これは周波数領域で実行してもよく、これは２つのスペクトルを乗算することに対応する。そのままでは、これは時間領域における円形の畳み込みの動作と等価である。したがって、フレームのテンポでのクリックによって（by clicks in tempo with the frames）聞くときに明らかにされる時間的エイリアシング(aliasing)による歪みを避けるために、いくつかの予防措置をとることが必要である。したがって、線形の畳み込み制約条件を満たすために、各入力フレームにある数の空サンプルを加えること（いわゆる「ゼロ・パディング」技術）と、雑音低減フィルタのインパルス応答の時間的サポートを制限すること（これは時間領域あるいは周波数領域で実行される）が、共に必要である。
【０１２５】
雑音低減フィルタのインパルス応答の時間的サポートを制限するために、時間領域に制約条件を導入することが可能であることに留意すべきであり、これは下記のことを必要とする。
ｉ）フィルタ
【０１２６】
【数４０】
【０１２７】
の伝達関数の知識にもとづいて、インパルス応答ｈ1（ｋ，ｎ）を得ることを可能にする第１の「逆」スペクトルの変換
ｉｉ）切り取られた（truncated）時間フィルタ
【０１２８】
【数４１】
【０１２９】
を得るために、このインパルス応答の点の数の制限。
ｉｉｉ）強制(constrained)インパルス応答
【０１３０】
【数４２】
【０１３１】
にもとづいて、フィルタ
【０１３２】
【数４３】
【０１３３】
の修正された伝達関数を得ることを可能にする第２の「直接」スペクトル変換。
【０１３４】
雑音低減フィルタ
【０１３５】
【数４４】
【０１３６】
の（周波数）伝達関数が利用できるので、（時間的）インパルス応答
【０１３７】
【数４５】
【０１３８】
は、逆離散的フーリエ変換（ＩＤＦＴ）によって得られる。すなわち、
【０１３９】
【数４６】
【０１４０】
このインパルス応答は次に、最上位のＬfilt1係数を選択し、ウインドウｗfilt1によって重みづけすることによって、時間的に制限される。すなわち
【０１４１】
【数４７】
【０１４２】
雑音低減フィルタの時間的サポートのこの制限は、二重の利点をもたらす。一方では、時間エイリアシング（線形畳み込みに従う）の問題の回避を可能にする。他方では、平滑化を保証し、あまりにも攻撃的なフィルタの作用の回避を可能にする。
【０１４３】
雑音低減フィルタ
【０１４４】
【数４８】
【０１４５】
の伝達関数は、したがって、離散的フーリエ変換（ＤＦＴ）
【０１４６】
【数４９】
【０１４７】
によって得られる。
【０１４８】
最後に、雑音のある信号の第１の、雑音が除去されたバージョンは、周波数濾波によって、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算すること、すなわち、
【０１４９】
【数５０】
【０１５０】
によって得られる
この周波数濾波ステップは、時間的濾波によって、すなわち、ｘw（ｋ，ｎ）を
【０１５１】
【数５１】
【０１５２】
によって直接濾波し、次にその結果の離散的フーリエ変換（ＤＦＴ）を行うことによって、等価な方法で実行してもよい。
【０１５３】
有用な信号
【０１５４】
【数５２】
【０１５５】
のこの第１の推定値は、「背景技術」で説明したように、雑音低減フィルタ
【０１５６】
【数５３】
【０１５７】
によってある高調波が抑圧されているので、可聴の歪みを含んでいる。勿論、歪みのレベルはＳＮＲに直接関係している。周波数領域の雑音が強ければ強いほど、この領域の高調波は劣化しやすい。
【０１５８】
この欠点を緩和するために、図２に示す装置１の機能ブロック８が、破壊され、すなわち非常に大きく減衰させられた高調波を有する疑似信号
【０１５９】
【数５４】
【０１６０】
を生成する。
【０１６１】
信号
【０１６２】
【数５５】
【０１６３】
は、時間領域であるいは周波数領域で得てもよい。時間領域では、信号に非線形性、たとえば、単一波の修正関数、絶対値、閾値に対して最小値／最大値等の中の関数、を加えることによって高調波性（harmonicity）を回復することが可能である。図３のステップ１４に示すように、ＮＬが、使用される非線形の関数を表すとすれば、
【０１６４】
【数５６】
【０１６５】
が得られる。ここで、
【０１６６】
【数５７】
【０１６７】
は、図３のステップ１３にしたがって、逆離散的フーリエ変換（ＩＤＦＴ）によって得られる。すなわち、
【０１６８】
【数５８】
【０１６９】
したがって、図３のステップ１５に示すように、離散的フーリエ変換（ＤＦＴ）によって、
【０１７０】
【数５９】
【０１７１】
を得る。すなわち、
【０１７２】
【数６０】
【０１７３】
高調波は、有用な信号の位置と同じ位置に再現されることがわかるであろう。これは、高調波を復活させるために時間的非線形性が使用されるということによって、暗示的に保証される。
【０１７４】
さらに、非線形性が時間領域で用いられる場合には、xw(ｋ,ｎ)と
【０１７５】
【数６１】
【０１７６】
の畳み込みによって直接得られた信号
【０１７７】
【数６２】
【０１７８】
に非線形性を適用することによって、周波数領域への切り替えを回避することが可能であることがわかるであろう。
【０１７９】
図４は非線形性の効果と利点を示している。第１のスペクトル（曲線１６）は、雑音を含まない有用な信号のフレームに対応し、基準の役割をする。第２のスペクトル（曲線１７）は、信号
【０１８０】
【数６２】
【０１８１】
に対応する。この第２のスペクトルにおいて、いくつかの高調波は完全に破壊され、他のものは劣化したことが明らかである。最後のスペクトル（曲線１８）は、上述の信号
【０１８２】
【数６３】
【０１８３】
に対応する。信号
【０１８４】
【数６４】
【０１８５】
に適用された非線形性が、信号
【０１８６】
【数６５】
【０１８７】
内の失われた高調波の再現を効果的に可能にしたことが明らかである。この信号
【０１８８】
【数６６】
【０１８９】
のスペクトルは、したがって、従来の雑音除去技術が破壊する高調波を維持できるであろう、特に、新しい雑音低減フィルタを計算するのに非常に有用な高調波性情報を含んでいる。
【０１９０】
時間的非線形性は、信号
【０１９１】
【数６７】
【０１９２】
のスペクトルと、有用な信号と同じ高調波構造（高調波の位置のレベルにおいて）を含んでいる信号のスペクトルの間の円形の畳み込みを実行することにある等価な周波数を有している。このようなスペクトルは、たとえば、
【０１９３】
【数６８】
【０１９４】
の関数に対応する時間信号ｄ（ｋ、ｎ）によって得られる。したがって、信号
【０１９５】
【数６９】
【０１９６】
を周波数領域で、次のように、得ることも可能である。
【０１９７】
【数７０】
【０１９８】
ここで、記号
【０１９９】
【数７１】
【０２００】
は、円形の畳み込み演算子に対応する。選択された時間的非線形性が絶対値関数であれば、たとえば、信号ｄ（ｋ、ｎ）は次式
【０２０１】
【数７２】
【０２０２】
のように表現することが可能であり、関数ｓｇｎはそれが適用される値の符号を指している。
【０２０３】
したがって、信号
【０２０４】
【数７３】
【０２０５】
は、それが時間的非線形性によってあるいは周波数の円形の畳み込みによって得られる場合には、
【０２０６】
【数７４】
【０２０７】
の高調波が破壊され、あるいは劣化した高調波を有している。それは、雑音が除去された信号の高調波コム（harmonic comb）の本質的に良好な推定値を与え、高調波は有用な信号と同じ位置に維持されている。高調波コムのこの推定値は特に有用であり、特に、しかし限定されるものではなく、雑音除去の目的のために利用することができる。
【０２０８】
有利なことに、この信号は、以下に説明するように、従来のアルゴリズムによって通常破壊されてしまう高調波を維持できる雑音低減フィルタの計算を可能にすることができる。
【０２０９】
本発明は、信号
【０２１０】
【数７５】
【０２１１】
を得るために、以上に説明した例に限定されるものではないことがわかるであろう。実際に本発明は、雑音低減フィルタ
【０２１２】
【数７６】
【０２１３】
によって抑圧された高調波の再現を目的とする任意の他の方法に拡大される。
【０２１４】
有利なことに、第２のフィルタが図２の機能ブロック９によって、この後推定される。この第２の雑音低減フィルタ
【０２１５】
【数７７】
【０２１６】
の伝達関数は、図５に示すように、次のように、周波数領域で推定され、ｇは関数を表している。
【０２１７】
【数７８】
【０２１８】
フィルタの伝達関数に対するこの式は、雑音
【０２１９】
【数７９】
【０２２０】
のＰＳＤと高調波修復ステップによって推定された有用な信号のＰＳＤに依存する数式である。すなわち。
【０２２１】
【数８０】
【０２２２】
パラメータρ（ｋ、ｆ）は、信号
【０２２３】
【数８１】
【０２２４】
を生成するように選択された非線形性ＮＬの関数として、ＰＳＤ
【０２２５】
【数８２】
【０２２６】
の計算において、信号
【０２２７】
【数８３】
【０２２８】
の再投入のレベルを制御するために使用される。このパラメータは一定でも、あるいは、周波数および／または時間に依存してもよいことがわかる。
【０２２９】
関数ｇは、
【０２３０】
【数８４】
【０２３１】
が、たとえば、電力スペクトル減算規則（式（１））、振幅スペクトル減算規則（式（２））、あるいは、開ループ・ウイーナー・フィルタ規則（式（３））を実行するように、選択してよい。周波数ベースの任意の他の雑音抑圧規則を、
【０２３２】
【数８５】
【０２３３】
を推定するために、勿論実行してもよい。
【０２３４】
この雑音除去伝達関数
【０２３５】
【数８６】
【０２３６】
は、従来の雑音除去技術によって通常は破壊されてしまう信号の高調波を維持する利点を持っている。したがって、これは、雑音が除去された信号の歪みを制限する。このフィルタは、フィルタ
【０２３７】
【数８７】
【０２３８】
の良好な特性、すなわち、非定常性の良好な追跡、ならびに、残留音楽的雑音が非常に少ないこと、を保持していることが、さらに指摘されるべきである。
【０２３９】
この後に、第２の推定雑音低減フィルタにもとづいて、第２の濾波が、装置１の機能ブロック１０によって実行される。したがって、第１の伝達関数については、雑音低減フィルタ
【０２４０】
【数８８】
【０２４１】
の（周波数）伝達関数と等価な（時間的）インパルス応答
【０２４２】
【数８９】
【０２４３】
が、逆離散的フーリエ変換（ＩＤＦＴ）によって得られる。すなわち、
【０２４４】
【数９０】
【０２４５】
このインパルス応答はこの後、最上位のＬfilt2係数を選択し、ウインドウｗfilt2によって重みづけすることによって、時間的に有利に制限される。すなわち、
【０２４６】
【数９１】
【０２４７】
雑音低減フィルタ
【０２４８】
【数９２】
【０２４９】
の伝達関数はこの後、離散的フーリエ変換（ＤＦＴ）によって得られる。すなわち、
【０２５０】
【数９３】
【０２５１】
最後に、雑音が除去された周波数フレームが、周波数濾波によって、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算することによって得られる。すなわち、
【０２５２】
【数９４】
【０２５３】
この周波数濾波ステップは、時間的濾波によって、すなわち、ｘw（ｋ,ｎ）を
【０２５４】
【数９５】
【０２５５】
によって直接濾波することによって、等価に実行することもできる。この場合、時間領域で信号
【０２５６】
【数９６】
【０２５７】
を直接得る。
【０２５８】
したがって、濾波を周波数領域で実行する場合、時間領域に戻るのが望ましい。これは、逆スペクトル変換（ＩＤＦＴ）を適用することによって、装置１の機能ブロック１１によって実行される。すなわち、
【０２５９】
【数９７】
【０２６０】
雑音が除去された出力信号はこの後、ブロック・オーバーラップおよびアッド(block overlap and add)技術（ＯＬＡは「オーバーラップ・アッド」を意味する）によって、あるいは、ブロック・セーブ(block save)技術（ＯＬＳは「オーバーラップ・セーブ」を意味する）によって、有利なことに合成される。信号を時間領域において再構成するこの動作は、逆短期フーリエ変換（ＩＳＴＦＴ）と呼ばれる。
【０２６１】
これらのステップが完了すると、フレーム１２が装置１の出力に得られる。この出力フレーム１２は、歪みが抑えられている、入力フレーム２の、雑音が除去された推定値である。
【０２６２】
以降の説明では、本発明による具体化の有利な選択を説明する。これらの選択は代表的な実施形態であり、限定する性質を有すると見なされるべきではない。
【０２６３】
装置１の機能ブロック３および１１によって実行される分析と合成の具体化例に対しては、ＯＬＡ技術が選択されるのが有利である。さらに、分析されるフレーム同士の間にオーバーラップが存在する。たとえば、５０％のオーバーラップ率を選択してもよい。長さＬの、雑音のある信号ｘ（ｋ,ｎ）｛ｎ＝０，．．．．．．，Ｌ-１｝の現在のフレーム２はまず、信号ｘw（ｋ、ｎ）を送出する重みづけウインドウｗ（n）と乗算される。すなわち、
【０２６４】
【数９８】
【０２６５】
この例では、使用されたウインドウはサイズＬのハニング(Hanning)ウインドウである。すなわち、
【０２６６】
【数９９】
【０２６７】
ここで、Ｌ＝２５６であり、これは、サンプリング周波数Ｆｅ＝８kHzヘルツにおいて、長さ３２ｍｓのフレームに対応する。
【０２６８】
周波数領域への切り替えは、離散的フーリエ変換（ＤＦＴ）の数値的に最適化されたバージョンである高速フーリエ変換（ＦＦＴ）を使用して有利に実行される。より良好な周波数分解能を得るために、また、線形の畳み込み制約条件を満たすために、長さＬFFT＝５１２の高速フーリエ変換が使用される。したがって、ウインドウ化されたブロックｘw
（ｋ,ｎ）を２５６個の空サンプル（「ゼロパディング」）によって長くすることが、あらかじめ必要である。すなわち、
【０２６９】
【数１００】
【０２７０】
入力フレームのＦＦＴ Ｘ（ｋ,ｆ）が次に、次式
【０２７１】
【数１０１】
【０２７２】
によって得られる。
【０２７３】
ボーカルアクティビティの検出は、機能ブロック４によって、任意に実行することができる。
【０２７４】
機能ブロック５によって実行される、雑音電力スペクトル密度
【０２７５】
【数１０２】
【０２７６】
の推定に関する限り、次の方法で有利に進むことが可能である。ＤＶＡによって検出される「雑音のみ」の各フレームｋbに対して、雑音電力スペクトル密度
【０２７７】
【数１０３】
【０２７８】
が、次の再帰式
【０２７９】
【数１０４】
【０２８０】
によって推定される。
【０２８１】
現在のフレームが雑音フレームとして検出されなければ、雑音の電力スペクトル密度の推定は固定される。平滑化量（smoothing quantity）αは一定に選択され、α＝０．８８２５に等しい。これは、１２８ｍｓの時定数（５０％オーバーラップと、サンプリング周波数Ｆｅ＝８ｋHzに関係する）に対応し、信頼できる推定と雑音の統計の時間的変化の追跡との間の妥協を保証するに十分であると判断される。
【０２８２】
機能ブロック７によって推定された第１の雑音低減フィルタの伝達関数
【０２８３】
【数１０５】
【０２８４】
は、
【０２８５】
【数１０６】
【０２８６】
のように、有利なことに開ループ・ウイーナー・フィルタの伝達関数である。
【０２８７】
ここで、η（ｋ,ｆ）は、有用な信号
【０２８８】
【数１０７】
【０２８９】
のＰＳＤと雑音
【０２９０】
【数１０８】
【０２９１】
のＰＳＤとの推定値の比率によって理論的に定められた事前の信号対雑音比（ＳＮＲ）を表す。すなわち、
【０２９２】
【数１０９】
【０２９３】
有用な信号
【０２９４】
【数１１０】
【０２９５】
のＰＳＤは、有向決定推定値（directed-decision estimator）、すなわち、
【０２９６】
【数１１１】
【０２９７】
を使用して有利に得られる。ここで、重心パラメータβは一定に選択され、０．９８に等しい。上述したように、関数Ｐは、（１−β）によって重みづけされた量の閾値化（thresholding）を保証する。すなわち、
【０２９８】
【数１１２】
【０２９９】
伝達関数
【０３００】
【数１１３】
【０３０１】
はその後、有用な信号
【０３０２】
【数１１４】
【０３０３】
のＰＳＤの推定値を改良するために、有利に再使用される。したがって、
【０３０４】
【数１１５】
【０３０５】
によって量
【０３０６】
【数１１６】
【０３０７】
を得、これは、
【０３０８】
【数１１７】
【０３０９】
によって与えられる事前の信号対雑音比η2（ｋ,ｆ）の第２の推定を行うために使用される。
【０３１０】
次に、第２のパスが、η2（ｋ,ｆ）にもとづいて第１の雑音低減フィルタの伝達関数
【０３１１】
【数１１８】
【０３１２】
の推定値の計算に有利に存在する。すなわち、
【０３１３】
【数１１９】
【０３１４】
機能ブロック６によって実行される濾波は次に、以下に示すように実行される。雑音低減フィルタ
【０３１５】
【数１２０】
【０３１６】
の（周波数）伝達関数が利用可能であるので、対応する時間応答
【０３１７】
【数１２１】
【０３１８】
は、逆離散的フーリエ変換（ＩＤＦＴ）の数値的に最適化されたバージョンである逆高速フーリエ変換（ＩＦＦＴ）を使用することによって、有利に得られる。すなわち、
【０３１９】
【数１２２】
【０３２０】
信号
【０３２１】
【数１２３】
【０３２２】
が実数であるので、対応する時間フィルタがまず原因である(causal)とされる。意図された用途にとって重要なサンプルに対応するこのフィルタＬfilt1＝２５６係数がその後選択される。その後に、このようにして得られたインパルス応答は、長さＬfilt1のハニングウインドウによって重み付けされる。最後に、重み付けされたインパルス応答は、エイリアシング無しの線形の畳み込み制約条件を満たすインパルス応答
【０３２３】
【数１２４】
【０３２４】
を与えるために、ＬFFT−Ｌfilt1個の零で補充される。
【０３２５】
雑音低減フィルタの伝達関数
【０３２６】
【数１２５】
【０３２７】
はその後、離散的フーリエ変換（ＤＦＴ）の数値的に最適化されたバージョンである高速フーリエ変換（ＦＦＴ）を有利に使用することで得られる。すなわち、
【０３２８】
【数１２６】
【０３２９】
最後に、雑音のある信号の第１の、雑音が除去されたバージョンは、周波数濾波、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算することによって得られる。すなわち、
【０３３０】
【数１２７】
【０３３１】
機能ブロック８によって実行される、失われた高調波を回復するステップにおいて、信号
【０３３２】
【数１２８】
【０３３３】
が疑似信号
【０３３４】
【数１２９】
【０３３５】
を生成するためにその後に使用される。まず、時間的信号
【０３３６】
【数１３０】
【０３３７】
が、逆高速フーリエ変換（ＦＦＴＩ）を有利に使用することによって得られる。すなわち、
【０３３８】
【数１３１】
【０３３９】
その後、単一波の非線形修正関数が、信号に高調波性を回復させる目的で、この信号に適用される。すなわち、
【０３４０】
【数１３２】
【０３４１】
勿論、単一波の修正以外の（絶対値、閾値に対して最大／最小、等）任意の他の非線形性を、上に示したように、使用することができる。最後に、信号
【０３４２】
【数１３３】
【０３４３】
が、高速フーリエ変換（ＦＦＴ）を有利に使用することによって得られる。すなわち、
【０３４４】
【数１３４】
【０３４５】
この信号は、第２の雑音低減フィルタの伝達関数を有利に計算するために、その後、使用される。
【０３４６】
この代表的な実施形態において、このような第２のフィルタの推定は、装置１の機能ブロック９によって、次のように実行される。第２の雑音低減フィルタの伝達関数
【０３４７】
【数１３５】
【０３４８】
は、
【０３４９】
【数１３６】
【０３５０】
であるように、開ループ・ウイーナー・フィルタの伝達関数である。ここで、ηharmo（ｋ,ｆ）は、有用な信号
【０３５１】
【数１３７】
【０３５２】
のＰＳＤの推定値と雑音
【０３５３】
【数１３８】
【０３５４】
のＰＳＤの推定値との比率によって理論的に定められた事前の信号対雑音比（ＳＮＲ）を表す。すなわち、
【０３５５】
【数１３９】
【０３５６】
ここで、有用な信号のＰＳＤ、
【０３５７】
【数１４０】
【０３５８】
は
【０３５９】
【数１４１】
【０３６０】
のように得られる。
【０３６１】
信号
【０３６２】
【数１４２】
【０３６３】
の再投入のレベルを制御するパラメータρ（ｋ,ｆ）は、一定として選択され、意図された用途に対して０．５に等しい。
【０３６４】
雑音低減フィルタ
【０３６５】
【数１４３】
【０３６６】
の（周波数）伝達関数が利用可能であるので、対応する時間応答
【０３６７】
【数１４４】
【０３６８】
は、逆離散的フーリエ変換（ＩＤＦＴ）の数値的に最適化されたバージョンである逆高速フーリエ変換（ＦＦＴＩ）を有利に使用することによって得られる。すなわち、
【０３６９】
【数１４５】
【０３７０】
信号
【０３７１】
【数１４６】
【０３７２】
が実数であるので、対応する時間フィルタはまず原因であるとされる。この用途に対して重要なサンプルに対応するこのフィルタのＬfilt2＝２５６個の係数は、その後に選択される。このようにして得られたインパルス応答はその後長さＬfilt2のハニングウインドウによって重み付けされる。最後に、この重み付けされたインパルス応答は、エイリアシング無し線形の畳み込み制約条件を満たすインパルス応答
【０３７３】
【数１４７】
【０３７４】
を与えるために、ＬFFT＝Ｌfilt2個の零で補充される。
【０３７５】
雑音低減フィルタの伝達関数
【０３７６】
【数１４８】
【０３７７】
はその後に、離散的フーリエ変換（ＤＦＴ）の数値的に最適化されたバージョンである高速フーリエ変換（ＦＦＴ）を有利に使用することによって得られる。すなわち、
【０３７８】
【数１４９】
【０３７９】
雑音低減フィルタ
【０３８０】
【数１５０】
【０３８１】
の係数の値は、雑音除去が強力過ぎることを回避するために、任意に限界値を設定する（thresheld）ことが可能である。
【０３８２】
最後に、雑音が除去された周波数フレームは、機能ブロック１０によって実行される周波数濾波によって、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算することによって得られる。すなわち、
【０３８３】
【数１５１】
【０３８４】
このようにして得られた信号の時間領域への戻りは、逆高速フーリエ変換（ＦＦＴＩ）を有利に使用することによって、逆スペクトル変換によって、機能ブロック１１のレベルにおいて、その後に実行される。すなわち、
【０３８５】
【数１５２】
【０３８６】
雑音が除去された出力信号はその後、ブロック・オーバーラップおよびアッド技術ＯＬＡによって合成される。
【０３８７】
さらに、図２に示す装置１は、１つの音声ピックアップ装置、１つの通信装置、あるいは１つの音声認識装置のような、意図される用途の機能として、特定の装置の中に設置してもよいことに留意されるであろう。
【０３８８】
さらに、本装置によって実行される処理は、たとえば通信端末の場合のようにローカルに、あるいは、たとえば音声認識サーバの場合にネットワーク内に集中化して実行してもよいことが留意される。
【０３８９】
さらに、実際には、本発明による雑音のある音声信号を処理する装置１は、本発明による信号を処理する方法を実施するのに適した（電子的）ハードウェア手段および／またはソフトウェア手段から成る。
【０３９０】
好ましい具体化例によれば、本発明による雑音のある音声信号を処理する方法のステップは、本発明によるこのような装置の１つで使用されるコンピュータ・プログラムの命令によって決定される。
【０３９１】
本発明による方法は、前述のプログラムが、装置に組み込まれた計算手段にロードされたときに実行され、次にその動作はプログラムの実行によって制御される。
【０３９２】
ここで、表現「コンピュータ・プログラム」は、（ソフトウェア）集合を形成する１つまたは複数のコンピュータ・プログラムを意味すると理解される。ソフトウェア集合の目的は、適切な計算手段によって実行される時に、本発明を実施することである。
【０３９３】
したがって、本発明の主題も、このようなコンピュータ・プログラムであり、特に、情報媒体に格納された１つのソフトウェアの形式のコンピュータ・プログラムである。このような情報媒体は、本発明によるプログラムを格納することが可能な任意の実体あるいは装置からなってよい。
【０３９４】
たとえば、当該媒体は、ＲＯＭ、たとえばＣＤＲＯＭ、あるいはマイクロ電子回路ＲＯＭ、あるいは磁気記録手段、たとえばハードディスク、のようなハードウェア記憶手段を有してもよい。変形例として、情報媒体は、プログラムが組み込まれ、当該方法の実行に適した、あるいは当該方法の実行に使用されることに適した集積回路であってもよい。
【０３９５】
さらに、情報媒体は、無線によってあるいは他の手段によって、電気ケーブルあるいは光ケーブルによって伝送可能な電気信号あるいは光信号のような伝送可能な無形の媒体であってもよい。本発明によるプログラムは、特に、インターネット形式のネットワークからダウンロードすることができる。
【０３９６】
設計の視点から、本発明によるコンピュータ・プログラムは、任意のプログラミング言語を使用することが可能であり、ソース・コード、オブジェクトコード、あるいは、ソース・コードとオブジェクトコードの間の中間コード（たとえば部分的にコンパイルされた形式）の形式、あるいは、本発明による方法を実施するのに望ましい任意の他の形式とすることが可能である。
【図面の簡単な説明】
【０３９７】
【図１】すでに説明済みであるが、３つの公知の短期間の抑圧規則の挙動を比較する図である。
【図２】本発明による実施形態を実行することが可能な、さまざまな機能ブロックを示す図である。
【図３】雑音が除去された信号内の失われた高調波の回復の例を示す図である。
【図４】本発明の一実施形態による信号内の高調波の回復の効果を示す図である。
【図５】本発明の一実施形態による第２の雑音低減フィルタの推定の例を示す図である。
【技術分野】
【０００１】
本発明は、雑音のある音声信号の処理に関する。本発明は特に、そのような信号に存在する雑音の低減に関する。
【背景技術】
【０００２】
音声信号の中の雑音、すなわち妨害信号、を低減する技術が公知である。それらは、音声信号の品質と明瞭性を改善するために、音声信号が出現する音響的環境を考慮することを目的としている。これらの技術は、考慮している音声信号から、この雑音のある信号に処理を行うことによって、有用な情報を抽出することにある。このような技術は、音声信号が何人かの発言者の間で送信される、電話、電話会議、ビデオ会議のような用途における、たとえば会話通信に適用される。このような技術は、さらに、雑音のある環境での音声収集の用途に、あるいは、音声信号が雑音に満ちた環境で発音された場合に性能が大幅に変化してしまう音声認識の用途に適用される。
【０００３】
これらの技術は通常、雑音低減フィルタの伝達関数を推定し、次に、スペクトル領域での乗算にもとづいて、濾波・プロセスを実行することにある。これらの技術は、「短期間スペクトル減衰による雑音低減」と呼ばれる方法に含まれる。
【０００４】
これらの技術によれば、考慮している音声信号ｘ（ｎ）は、有用な信号成分ｓ（ｎ）と雑音成分ｂ（ｎ）を含み、ここでｎは離散時間における時間インデックスを表している。しかし、連続した時間での信号の表現も採用することも可能なことが気づくであろう。信号ｘ（ｎ）は、一定の長さと、インデックスｋの連続した複数のフレームｘ（ｎ、ｋ）として構成される。これらのフレームのそれぞれは、まず重みづけウインドウと乗算され、雑音低減フィルタの計算に必要なスペクトル量の、その後の推定を改善することを可能にする。このようにウインドウ化された各フレームは次に、たとえば離散的あるいは高速フーリエ変換を用いて、スペクトル領域において分析される、この動作は、短期間フーリエ変換（ＳＴＦＴ）と呼ばれる。
【０００５】
注目する信号の、このようにして得られた周波数表現Ｘ（ｋ，ｆ）（ここで、ｆは周波数インデックス）は同一の時間において、雑音低減フィルタの伝達関数Ｈ（ｋ、ｆ）を推定することと、この伝達関数と雑音のある信号の短期間スペクトルの簡単な乗算によって、このフィルタをスペクトル領域で用いることを可能にする。したがって、濾波の結果は次式
【０００６】
【数１】
【０００７】
のように、記述することができる。
【０００８】
次に、得られた信号の時間領域への戻りが、逆スペクトル変換によって実行される。対応する時間信号は、ブロック・オーバーラップ・アッド技術（ＯＬＡは「オーバーラップ・アッド」を表わす）によって、あるいは、ブロック・セーブ技術（ＯＬＳは「オーバーラップ・セーブ」を表わす）によって、最終的に合成される。信号を時間領域においで再構築するこの動作は、逆短期間フーリエ変換（ＩＳＴＦＴ）と呼ばれる。
【０００９】
短期間スペクトル減衰方法の詳細な説明が、下記の参考文献に見られる。J.S. Lim、A.V. Oppenheim、「雑音のある通話のエンハンスメントと帯域幅圧縮」、Proceedings of the IEEE、Vol. 67、pp. 1586-1604、1979、および、 R.E. Crochiere、L.R. Rabiner、「マルチレート・ディジタル信号処理」、プレンティス・ホール、１９８３。
【００１０】
周波数成分ｆを有する、インデックスｋの時間セグメントにわたる観察信号Ｘ（ｋ、ｆ）に適用される、短期間スペクトル減衰Ｈ（ｋ、ｆ）は通常、ローカルな信号対雑音比ＳＮＲ（ｋ、ｆ）の推定値にもとづいて求められる。すべての抑圧規則に共通する特性は、次式
【００１１】
【数２】
【００１２】
によって与えられる、それらの漸近的な挙動にある。
【００１３】
多くの技術において、次の仮定がなされている。すなわち、雑音と有用な信号は統計的に無相関であり、有用な信号は断続的（無音期間の存在）であり、人間の耳は信号の位相（通常、処理によって修正されない）に気づかない。
【００１４】
一般に使用される抑圧規則のうちで、電力スペクトル減算、振幅スペクトル減算、およびウイーナー・フィルタの直接実施を、例として挙げることができる。これらの規則に対して、有用な音声信号(speech signal)の周波数成分ｆの短期間の推定値は、それぞれ次のように記述することができる。
【００１５】
電力スペクトル減算に対して（ J.S. Lim および A.V. Oppenheimによる前述の論文を参照のこと）
【００１６】
【数３】
【００１７】
振幅スペクトル減算に対して（ S.F. Boll「スペクトル減算を用いた、会話における音響雑音の抑圧」、IEEE Trans. on Audio, Speech and Signal Processing、Vol. 27、No. 2、pp. 113-120、April 1979、参照）。
【００１８】
【数４】
【００１９】
ウイーナー・濾波に対して（ J.S. Lim および A.V. Oppenheimによる前述の論文を参照のこと）
【００２０】
【数５】
【００２１】
これらの式において、
【００２２】
【数６】
【００２３】
と
【００２４】
【数７】
【００２５】
は、それぞれ、インデックスｋの時間ウインドウにわたる、観察信号Ｘ（ｋ、ｆ）の周波数成分ｆ内に存在する有用な信号の電力スペクトル密度と、雑音の電力スペクトル密度を表す。
【００２６】
雑音のある信号に適用されたスペクトル減衰の挙動を、上記の式にもとづいて、所与の周波数成分ｆに対して測定したローカルな信号対雑音比の関数として、検討することが可能である。上述の短期間の３つの抑圧規則について、これらの曲線が図１にプロットされている。ローカルな信号対雑音比が大きい場合には（図１の右側の部分）、各規則は実質的に同じ減衰を与えることがわかる。ガウス・モデルに対して最大の尤度の意味で最適な、最適な電力減算規則（ 0. Cappe、「 Ephraim および Malah 雑音サプレッサーを使用するミュージカル・ノイズ現象の除去」、IEEE Trans. on Speech and Audio Processing、Vol. 2、No. 2、pp. 345-349、April 1994、参照）は、依然として、雑音の電力が処理の出力において最上位の規則である。３つの抑圧規則に対して、カットオフ値に関するローカルな信号対雑音比の僅かな変化が、完全な減衰（Ｈ（ｋ、ｆ）≒０）から、スペクトルの無視できる修正（Ｈ（ｋ、ｆ）≒１）への切り替えで十分であることを指摘できる。
【００２７】
この後者の性質は、「音楽的雑音(musical noise)」と呼ばれる現象の原因の１つである。具体的にいうと、同時に決定論的、かつランダムな成分を有する周囲騒音は、ボーカル非アクティビティ（vocal non-activity）の期間のみを特徴とすることができる。ランダム成分の存在のために、ボーカルアクティビティ（vocal activity ）の期間の雑音の周波数成分ｆの実際の寄与と、ボーカル非アクティビティの各瞬間の間のいくつかのフレームにわたって行われたその平均の推定値との間には非常に大きい違いがある。この違いのために、ローカルな信号対雑音比の推定値は、カットオフ・レベルの付近で変動することがあり、したがって、現れては消え、かつ平均の寿命が、考慮している分析ウインドウの大きさを統計的に超えないスペクトル成分を処理の出力に生ずることがある。通過帯域の全体にわたるこの挙動の一般化は、可聴で不快な残留雑音を招く。
【００２８】
この残留雑音の影響を減らすいくつかの研究が試みられてきた。推薦された解決策は、いくつか方法を取っている。すなわち、短期間の推定値の平均化（ S.F. Boll による前述の論文を参照）、雑音電力スペクトルの過大推定（ M. Berouti 他、「音響雑音によって損なわれた通話のエンハンスメント」、Int. Conf. on Speech, Signal Processing, pp. 208-211, 1979、および、P. Lockwood、 J. Boudy、「非線形スペクトル減算器、隠れマルコフ・モデルおよび射影演算を使用する自動車内の頑強な音声認識のための実験」Proc. of EUSIPCO'91, pp. 79-82, 1991、参照）、あるいは、雑音スペクトル密度の最小値の追跡（ R. Martin, 「最小の統計にもとづくスペクトル減算」 Signal Processing VII: Theories and Applications, EUSIPCO'94, pp. 1182-1185, September 1994、参照）。
【００２９】
音楽的雑音を抑制する比較的に有効な解決策は、「有向決定」（directed-decision）と呼ばれる有用な信号の電力スペクトル密度の推定値から成る。（ Y. Ephraim, and D. Malah、「最小平均２乗誤差短期間スペクトル振幅推定値を用いた通話向上」、IEEE Trans. on Audio, Speech and Signal Processing、Vol. 32、No. 6、pp. 1109-1121、1984、および O.Cappe による前述の論文を参照）。この推定値は、有用な信号の瞬間的電力スペクトル密度と長期間の電力スペクトル密度の間で妥協を行い、その結果、音楽的雑音を効果的に除去するのが可能になる。この推定値に固有の遅れを補償することによって、この解決策を改善することも公知である（フランス特許出願公開明細書第２８２０２２７号、および C. Plapous、C. Marro、L. Mauuary、P. Scalart、「２ステップ雑音低減技術」、ICASSP, May 2004、参照）。
【００３０】
いくつかの研究も、通話と付加的な雑音信号の統計モデルにもとづく新しい抑圧規則を確立することに関している。これらの研究は、従来の方法に対して付加的な自由度を有するので、「ソフト・ディシジョン」アルゴリズムと呼ばれる新しいアルゴリズムの導入を可能にした。（ R.J. Mac Aulay、M.L. Malpass、「ソフト・ディシジョン雑音抑圧フィルタを用いた通話向上」、IEEE trans. on Audio, Speech and Signal Processing、Vol. 28、No. 2、pp. 137-145、April 1980、Y. Ephraim、 D. Malah、「最適な非線形スペクトル振幅推定値を使用する通話向上」、Int. Conf. on Speech, Signal Processing、pp. 1118-1121、1983、および、Y. Ephraim、D. Malah の上述の論文、「最小平均２乗誤差短期間スペクトル振幅推定値を使用する通話向上」、参照）。
【００３１】
上述したように、短期間スペクトル減衰の計算は、スペクトル成分のそれぞれに対する信号対雑音比の推定に依存している。例として、上記した式は、それぞれ次の量を含んでいる。
【００３２】
【数８】
【００３３】
このように、特に、歪みに関して、また、雑音レベルの効果的低減に関しての、雑音低減技術の性能は、信号対雑音比のこの推定値の妥当性に左右される。
【００３４】
この欠点は、公知の通話雑音除去システムの主要な限界となっている。具体的にいうと、現在の雑音除去システムは、信号対雑音比が低すぎることを特徴とする高調波の雑音除去は不可能である。実際には、雑音除去アルゴリズムは、各周波数に通話成分が存在するか、あるいは存在しないかを検出するために、ＳＮＲを使用している。推定ＳＮＲがあまりにも好ましくないものであれば、アルゴリズムは信号成分がないとみなし、それを抑圧する。このように、公知の雑音除去システムによって高調波が損なわれることがあるが、このような高調波が存在しなければならないことは、先験的に公知である。ここで、大部分の言語において、口に出された音（voiced sounds）（高調波）は、発声された音（sounds uttered）の非常に大きい部分を表していることに、注目すべきである。
【発明の開示】
【発明が解決しようとする課題】
【００３５】
本発明の目的は、公知の雑音除去システムの限界を克服することである。
【００３６】
本発明の他の目的は、雑音除去方法の性能を改善することである。
【００３７】
本発明の他の目的は、信号を過度に歪ませない音声信号処理を提案することである。特に、実行される信号の処理は、この信号に含まれる高調波のすべてあるいは一部の保存を可能にする。
【００３８】
本発明の他の目的は、音声信号処理の完了に際して、音楽的雑音の出現を制限することである。
【００３９】
本発明の他の目的は、有用な信号の高調波コム（harmonic comb）の良好な推定値を得ることである。
【課題を解決するための手段】
【００４０】
本発明は、連続したフレームとして構成された雑音のある音声信号の処理方法を提案する。該処理は、前記フレームの少なくとも１つに関して、
雑音のある音声信号の前記フレームに、周波数領域への変換を適用するステップと、
前記フレームに対する雑音の電力スペクトル密度を推定するステップと、
雑音の推定された電力スペクトル密度と、フレームに対応する有用な信号の電力スペクトル密度の推定値とにもとづいて、第１の雑音低減フィルタを計算するステップと、
前記フレームの第１の、雑音が除去された推定値を得るために、計算された前記第１の雑音低減フィルタを利用して、雑音のある音声信号の前記フレームを濾波するステップと、
前記雑音のある音声信号の前記フレームに対応する前記有用な信号と実質的に同じ位置にある高調波を有する第２の信号のフレームを、雑音のある音声信号の前記フレームの第１の、雑音が除去された推定値にもとづいて得るステップと、
を有する。
【００４１】
したがって、雑音のある音声信号のこのような処理は、信号に適用された第１の雑音低減フィルタの出力において高調波性の再生を得ることを可能にする。第２の信号のこのようにして得られたフレームは、第１の濾波の間に出現した可能性のある、雑音が除去された信号のフレームの歪みを除去するように作られる。
【００４２】
この処理は、後での利用の主題である、有用な信号の高調波コムの良好な推定値を得ることを可能にする。
【００４３】
有利なことに、雑音除去の領域において、本方法はさらに次のステップ、すなわち、
前記雑音の推定された電力スペクトル密度と、前記フレームの第１の、雑音が除去された推定値の電力と得られた第２の信号の前記フレームの電力の組み合わせとにもとづいて、第２の雑音低減フィルタを計算するステップと、
前記フレームの第２の、雑音が除去された推定値を得るために、計算された第２の雑音低減フィルタを利用して、前記雑音のある音声信号の前記フレームを濾波するステップと、
前記フレームの第２の雑音が除去された推定値を合成するステップと
を有する。
【００４４】
この実施態様においては、第２の雑音低減フィルタは、それが、高調波再生から生ずる信号によって駆動されるので、高調波を保存するように計算される。第２の濾波の完了時に得られた前記フレームの第２の雑音が除去された推定値は、したがって、入力信号の高調波が破壊されるかあるいは少なくとも損なわれる従来の雑音除去システムで得られた推定値よりも有効である。
【００４５】
第１の雑音低減フィルタと、必要ならば、第２の雑音低減フィルタの計算は、短期間スペクトル減衰技術を実行する第１のパスを含んでもよい。実例として、たとえば、次の技術、電力スペクトル減算、振幅スペクトル減算、および開ループ・ウイーナー・フィルタ等を引用することができる。有利なことに、それは、たとえば、電力スペクトル減算、振幅スペクトル減算、および開ループ・ウイーナー・フィルタ等の短期間スペクトル減衰技術を実行する第２のパスをさらに有してもよく、また、第１のパスの間に実行される計算を考慮して、前記フレームに対応する有用な信号の電力スペクトル密度の推定値を使用してもよい。
【００４６】
有利なことに、雑音のある音声信号の前記フレームに対応する有用な信号と実質的に同じ位置にある高調波を有する第２の信号のフレームを得ることは、前記第１の推定値が時間領域内にある場合には、雑音のある音声信号の前記フレームの第１の、雑音が除去された推定値への非線形の関数の適用を含み、前記第１の推定値が周波数領域内にある場合には、雑音のある音声信号の前記フレームの第１の、雑音が除去された推定値と非線形の関数の周波数領域への変換の第１の、雑音が除去された推定値との間の円形の畳み込みの適用を含む。
【００４７】
前記非線形の関数は、たとえば、単一波の修正関数、絶対値、雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値と閾値との間の最大値、雑音のある音声信号の前記フレームの前記第１の雑音が除去された推定値と閾値との間の最小値、の中の１つであってもよい。
【００４８】
本発明は、上述の方法を実行するように設計された手段を有する、雑音のある音声信号を処理する装置をさらに提案する。
【００４９】
本発明は、プログラムが計算手段によってロードされ実行されるときに、上述の方法を実行するように作られた命令を有する、情報媒体上のコンピュータ・プログラムをさらに提案する。
【発明を実施するための最良の形態】
【００５０】
本発明による他の特徴と利点は、添付図面を参照した非限定的な代表的実施形態の下記の説明において明らかとなろう。
【００５１】
図２は、本発明による一実施形態による装置１を示す。雑音のある音声信号の現在のフレーム２が、処理を受けるために装置１に入力される。このフレームは、前述の表記法によって、Ｘ（ｋ，ｎ）と示される。
【００５２】
装置１は、役割が、雑音のある信号ｘ（ｋ、ｎ）の現在のフレームを周波数領域Ｘ（ｋ、ｆ）に切り替えることである分析ユニット３を有する。まず、フレームは、信号ｘw（
ｋ，ｎ）
【００５３】
【数９】
【００５４】
を送出する重みづけウインドウｗ（ｎ）と乗算される。
【００５５】
次に、離散的フーリエ変換（ＤＦＴ）、すなわち
【００５６】
【数１０】
【００５７】
を使用して、周波数領域への切り替えが実行される。
【００５８】
当該ＤＦＴは、高速フーリエ変換（すなわちＦＦＴ）によって実行するのが有利である。しかし、ウェーブレット変換のような、周波数領域への他の変換も可能である。同じことが、対応する逆動作、すなわち、後に説明する、時間領域に戻すための逆離散的フーリエ変換（ＩＤＦＴ）にも当てはまる。
【００５９】
さらに、ボーカルアクティビティを検出する機能４（ＤＶＡ）を、現在のフレーム２で実行するのが有利である。ＤＶＡは、雑音の電力スペクトル密度（ＰＳＤ）の推定値をいつ更新するべきかを知ることを可能にする。したがって、ＤＶＡによって検出された各「雑音のみの」フレームｋbに対して、雑音電力スペクトル密度
【００６０】
【数１１】
【００６１】
が下記の再帰式
【００６２】
【数１２】
【００６３】
にしたがって、機能ブロック５によって推定される。
【００６４】
パラメータα（ｋb）は平滑化係数を制御する。該パラメータ長時間にわたって変化することがある。
【００６５】
現在のフレーム２が雑音フレームとして検出されなければ、雑音の電力スペクトル密度の推定値は固定される、すなわち、以前の雑音フレームに対して得られたＰＳＤの最後の値が保持される。
【００６６】
【数１３】
【００６７】
の推定値は、この指数関数的平滑化推定値（exponential-smoothing estimator）に限定されず、任意の他の電力スペクトル密度推定値を使用できることを指摘しておく。
【００６８】
時間領域で動作するか、周波数領域で動作するかにかかわらず、任意の種類のＤＶＡを使用することができる。ボーカルアクティビティのこのような検出を不要にすることも可能である。
【００６９】
次に、第１の雑音低減フィルタの推定が、図２の機能ブロック７によって実行される。この第１の雑音低減フィルタの伝達関数は、２パス技術（two-pass technique）によって周波数領域で有利なことに推定される。（フランス特許出願公開明細書第２８２０２２７号、および、上述の C. Plapous、C. Marro、L. Mauuary、P. Scalart、「２ステップの雑音除去技術」、ICASSP, May 2004、を参照）。第１のパスにおいて、次の伝達関数を有するフィルタが計算される。
【００７０】
【数１４】
【００７１】
フィルタの伝達関数に対するこの式は、有用な信号
【００７２】
【数１５】
【００７３】
のＰＳＤと、雑音
【００７４】
【数１６】
【００７５】
のＰＳＤに依存する数式である。この意味で、
【００７６】
【数１７】
【００７７】
が、たとえば、電力スペクトル減算規則（式（１））、振幅スペクトル減算規則（式（２））、開ループ・ウイーナー・フィルタ規則（式（３））等を実行するように、関数ｆstep1を短期間スペクトル減衰技術にしたがって選択してもよい。周波数ベースの任意の他の雑音抑圧規則も、
【００７８】
【数１８】
【００７９】
を推定するために、実行してもよい。
【００８０】
雑音電力スペクトル密度
【００８１】
【数１９】
【００８２】
の計算を、以上に詳細に説明した。スペクトル量
【００８３】
【数２０】
【００８４】
に関しては、それは、ボーカルアクティビティ期間の間、信号と雑音の混合のために、直接得ることはできない。それを計算するためには、次式
【００８５】
【数２１】
【００８６】
にしたがって有向決定推定値（directed-decision estimator）が使用される。（上述の、Y. Ephraim、 D. Malah、「最小平均２乗誤差短期スペクトル振幅推定値を使用する通話エ向上」を参照）。ここで、β（ｋ）は、時間と共に変化することがある重心パラメータであり、
【００８７】
【数２２】
【００８８】
は、インデックスｋ−１の以前のフレームに関して推定された有用な信号のスペクトルである。推定誤差の場合には負になる可能性がある量
【００８９】
【数２３】
【００９０】
の閾値処理（thresholding）を保証する関数Ｐは、次式によって与えられる。
【００９１】
【数２４】
【００９２】
【数２５】
【００９３】
の推定は、この有向決定推定値に限定されないことを指摘しておく。具体的にいうと、指数関数的平滑化推定値あるいは任意の他の電力スペクトル密度推定値を使用してもよい。
【００９４】
伝達関数
【００９５】
【数２６】
【００９６】
はその後、有用な信号
【００９７】
【数２７】
【００９８】
のＰＳＤの推定を改良するために再使用してもよい。次に、量
【００９９】
【数２８】
【０１００】
を、次式によって得る。
【０１０１】
【数２９】
【０１０２】
次に、第２のパスは、
【０１０３】
【数３０】
【０１０４】
にもとづいて、第１の雑音低減フィルタの伝達関数の推定値
【０１０５】
【数３１】
【０１０６】
、すなわち、
【０１０７】
【数３２】
【０１０８】
を計算することである。
【０１０９】
独特の特徴が、有用な信号
【０１１０】
【数３３】
【０１１１】
のＰＳＤの「高速な」更新にある、この２パス計算は、第１の雑音低減フィルタ
【０１１２】
【数３４】
【０１１３】
に２つの利点をもたらす。一方では、有用な信号の非定常性を高速に、特にその時間的包絡線（たとえば、無音／通話遷移の間の通話信号の開始あるいは終了（attacks or extinctions））の高速な変化の間に、追跡することが可能である。他方では、雑音低減フィルタはより良好に推定され、これは、本方法の性能の向上によって明らかにされる。
【０１１４】
【数３５】
【０１１５】
および
【０１１６】
【数３６】
【０１１７】
は、たとえば、電力スペクトル減算規則（式（１））、振幅スペクトル減算規則（式（２））、あるいは、開ループ・ウイーナー・フィルタ規則（式（３））のような、短期間スペクトル減衰技術を実行してもよい。周波数ベースの何らかの他の雑音抑圧規則も、
【０１１８】
【数３７】
【０１１９】
および
【０１２０】
【数３８】
【０１２１】
を推定するために、実行してもよい。勿論、第２のバスを実行しないで、第１のパスのみに制限することも可能である。
【０１２２】
次に、第１の計算されたフィルタ
【０１２３】
【数３９】
【０１２４】
にしたがって、装置１の機能ブロック６によって濾波が実行される。これは周波数領域で実行してもよく、これは２つのスペクトルを乗算することに対応する。そのままでは、これは時間領域における円形の畳み込みの動作と等価である。したがって、フレームのテンポでのクリックによって（by clicks in tempo with the frames）聞くときに明らかにされる時間的エイリアシング(aliasing)による歪みを避けるために、いくつかの予防措置をとることが必要である。したがって、線形の畳み込み制約条件を満たすために、各入力フレームにある数の空サンプルを加えること（いわゆる「ゼロ・パディング」技術）と、雑音低減フィルタのインパルス応答の時間的サポートを制限すること（これは時間領域あるいは周波数領域で実行される）が、共に必要である。
【０１２５】
雑音低減フィルタのインパルス応答の時間的サポートを制限するために、時間領域に制約条件を導入することが可能であることに留意すべきであり、これは下記のことを必要とする。
ｉ）フィルタ
【０１２６】
【数４０】
【０１２７】
の伝達関数の知識にもとづいて、インパルス応答ｈ1（ｋ，ｎ）を得ることを可能にする第１の「逆」スペクトルの変換
ｉｉ）切り取られた（truncated）時間フィルタ
【０１２８】
【数４１】
【０１２９】
を得るために、このインパルス応答の点の数の制限。
ｉｉｉ）強制(constrained)インパルス応答
【０１３０】
【数４２】
【０１３１】
にもとづいて、フィルタ
【０１３２】
【数４３】
【０１３３】
の修正された伝達関数を得ることを可能にする第２の「直接」スペクトル変換。
【０１３４】
雑音低減フィルタ
【０１３５】
【数４４】
【０１３６】
の（周波数）伝達関数が利用できるので、（時間的）インパルス応答
【０１３７】
【数４５】
【０１３８】
は、逆離散的フーリエ変換（ＩＤＦＴ）によって得られる。すなわち、
【０１３９】
【数４６】
【０１４０】
このインパルス応答は次に、最上位のＬfilt1係数を選択し、ウインドウｗfilt1によって重みづけすることによって、時間的に制限される。すなわち
【０１４１】
【数４７】
【０１４２】
雑音低減フィルタの時間的サポートのこの制限は、二重の利点をもたらす。一方では、時間エイリアシング（線形畳み込みに従う）の問題の回避を可能にする。他方では、平滑化を保証し、あまりにも攻撃的なフィルタの作用の回避を可能にする。
【０１４３】
雑音低減フィルタ
【０１４４】
【数４８】
【０１４５】
の伝達関数は、したがって、離散的フーリエ変換（ＤＦＴ）
【０１４６】
【数４９】
【０１４７】
によって得られる。
【０１４８】
最後に、雑音のある信号の第１の、雑音が除去されたバージョンは、周波数濾波によって、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算すること、すなわち、
【０１４９】
【数５０】
【０１５０】
によって得られる
この周波数濾波ステップは、時間的濾波によって、すなわち、ｘw（ｋ，ｎ）を
【０１５１】
【数５１】
【０１５２】
によって直接濾波し、次にその結果の離散的フーリエ変換（ＤＦＴ）を行うことによって、等価な方法で実行してもよい。
【０１５３】
有用な信号
【０１５４】
【数５２】
【０１５５】
のこの第１の推定値は、「背景技術」で説明したように、雑音低減フィルタ
【０１５６】
【数５３】
【０１５７】
によってある高調波が抑圧されているので、可聴の歪みを含んでいる。勿論、歪みのレベルはＳＮＲに直接関係している。周波数領域の雑音が強ければ強いほど、この領域の高調波は劣化しやすい。
【０１５８】
この欠点を緩和するために、図２に示す装置１の機能ブロック８が、破壊され、すなわち非常に大きく減衰させられた高調波を有する疑似信号
【０１５９】
【数５４】
【０１６０】
を生成する。
【０１６１】
信号
【０１６２】
【数５５】
【０１６３】
は、時間領域であるいは周波数領域で得てもよい。時間領域では、信号に非線形性、たとえば、単一波の修正関数、絶対値、閾値に対して最小値／最大値等の中の関数、を加えることによって高調波性（harmonicity）を回復することが可能である。図３のステップ１４に示すように、ＮＬが、使用される非線形の関数を表すとすれば、
【０１６４】
【数５６】
【０１６５】
が得られる。ここで、
【０１６６】
【数５７】
【０１６７】
は、図３のステップ１３にしたがって、逆離散的フーリエ変換（ＩＤＦＴ）によって得られる。すなわち、
【０１６８】
【数５８】
【０１６９】
したがって、図３のステップ１５に示すように、離散的フーリエ変換（ＤＦＴ）によって、
【０１７０】
【数５９】
【０１７１】
を得る。すなわち、
【０１７２】
【数６０】
【０１７３】
高調波は、有用な信号の位置と同じ位置に再現されることがわかるであろう。これは、高調波を復活させるために時間的非線形性が使用されるということによって、暗示的に保証される。
【０１７４】
さらに、非線形性が時間領域で用いられる場合には、xw(ｋ,ｎ)と
【０１７５】
【数６１】
【０１７６】
の畳み込みによって直接得られた信号
【０１７７】
【数６２】
【０１７８】
に非線形性を適用することによって、周波数領域への切り替えを回避することが可能であることがわかるであろう。
【０１７９】
図４は非線形性の効果と利点を示している。第１のスペクトル（曲線１６）は、雑音を含まない有用な信号のフレームに対応し、基準の役割をする。第２のスペクトル（曲線１７）は、信号
【０１８０】
【数６２】
【０１８１】
に対応する。この第２のスペクトルにおいて、いくつかの高調波は完全に破壊され、他のものは劣化したことが明らかである。最後のスペクトル（曲線１８）は、上述の信号
【０１８２】
【数６３】
【０１８３】
に対応する。信号
【０１８４】
【数６４】
【０１８５】
に適用された非線形性が、信号
【０１８６】
【数６５】
【０１８７】
内の失われた高調波の再現を効果的に可能にしたことが明らかである。この信号
【０１８８】
【数６６】
【０１８９】
のスペクトルは、したがって、従来の雑音除去技術が破壊する高調波を維持できるであろう、特に、新しい雑音低減フィルタを計算するのに非常に有用な高調波性情報を含んでいる。
【０１９０】
時間的非線形性は、信号
【０１９１】
【数６７】
【０１９２】
のスペクトルと、有用な信号と同じ高調波構造（高調波の位置のレベルにおいて）を含んでいる信号のスペクトルの間の円形の畳み込みを実行することにある等価な周波数を有している。このようなスペクトルは、たとえば、
【０１９３】
【数６８】
【０１９４】
の関数に対応する時間信号ｄ（ｋ、ｎ）によって得られる。したがって、信号
【０１９５】
【数６９】
【０１９６】
を周波数領域で、次のように、得ることも可能である。
【０１９７】
【数７０】
【０１９８】
ここで、記号
【０１９９】
【数７１】
【０２００】
は、円形の畳み込み演算子に対応する。選択された時間的非線形性が絶対値関数であれば、たとえば、信号ｄ（ｋ、ｎ）は次式
【０２０１】
【数７２】
【０２０２】
のように表現することが可能であり、関数ｓｇｎはそれが適用される値の符号を指している。
【０２０３】
したがって、信号
【０２０４】
【数７３】
【０２０５】
は、それが時間的非線形性によってあるいは周波数の円形の畳み込みによって得られる場合には、
【０２０６】
【数７４】
【０２０７】
の高調波が破壊され、あるいは劣化した高調波を有している。それは、雑音が除去された信号の高調波コム（harmonic comb）の本質的に良好な推定値を与え、高調波は有用な信号と同じ位置に維持されている。高調波コムのこの推定値は特に有用であり、特に、しかし限定されるものではなく、雑音除去の目的のために利用することができる。
【０２０８】
有利なことに、この信号は、以下に説明するように、従来のアルゴリズムによって通常破壊されてしまう高調波を維持できる雑音低減フィルタの計算を可能にすることができる。
【０２０９】
本発明は、信号
【０２１０】
【数７５】
【０２１１】
を得るために、以上に説明した例に限定されるものではないことがわかるであろう。実際に本発明は、雑音低減フィルタ
【０２１２】
【数７６】
【０２１３】
によって抑圧された高調波の再現を目的とする任意の他の方法に拡大される。
【０２１４】
有利なことに、第２のフィルタが図２の機能ブロック９によって、この後推定される。この第２の雑音低減フィルタ
【０２１５】
【数７７】
【０２１６】
の伝達関数は、図５に示すように、次のように、周波数領域で推定され、ｇは関数を表している。
【０２１７】
【数７８】
【０２１８】
フィルタの伝達関数に対するこの式は、雑音
【０２１９】
【数７９】
【０２２０】
のＰＳＤと高調波修復ステップによって推定された有用な信号のＰＳＤに依存する数式である。すなわち。
【０２２１】
【数８０】
【０２２２】
パラメータρ（ｋ、ｆ）は、信号
【０２２３】
【数８１】
【０２２４】
を生成するように選択された非線形性ＮＬの関数として、ＰＳＤ
【０２２５】
【数８２】
【０２２６】
の計算において、信号
【０２２７】
【数８３】
【０２２８】
の再投入のレベルを制御するために使用される。このパラメータは一定でも、あるいは、周波数および／または時間に依存してもよいことがわかる。
【０２２９】
関数ｇは、
【０２３０】
【数８４】
【０２３１】
が、たとえば、電力スペクトル減算規則（式（１））、振幅スペクトル減算規則（式（２））、あるいは、開ループ・ウイーナー・フィルタ規則（式（３））を実行するように、選択してよい。周波数ベースの任意の他の雑音抑圧規則を、
【０２３２】
【数８５】
【０２３３】
を推定するために、勿論実行してもよい。
【０２３４】
この雑音除去伝達関数
【０２３５】
【数８６】
【０２３６】
は、従来の雑音除去技術によって通常は破壊されてしまう信号の高調波を維持する利点を持っている。したがって、これは、雑音が除去された信号の歪みを制限する。このフィルタは、フィルタ
【０２３７】
【数８７】
【０２３８】
の良好な特性、すなわち、非定常性の良好な追跡、ならびに、残留音楽的雑音が非常に少ないこと、を保持していることが、さらに指摘されるべきである。
【０２３９】
この後に、第２の推定雑音低減フィルタにもとづいて、第２の濾波が、装置１の機能ブロック１０によって実行される。したがって、第１の伝達関数については、雑音低減フィルタ
【０２４０】
【数８８】
【０２４１】
の（周波数）伝達関数と等価な（時間的）インパルス応答
【０２４２】
【数８９】
【０２４３】
が、逆離散的フーリエ変換（ＩＤＦＴ）によって得られる。すなわち、
【０２４４】
【数９０】
【０２４５】
このインパルス応答はこの後、最上位のＬfilt2係数を選択し、ウインドウｗfilt2によって重みづけすることによって、時間的に有利に制限される。すなわち、
【０２４６】
【数９１】
【０２４７】
雑音低減フィルタ
【０２４８】
【数９２】
【０２４９】
の伝達関数はこの後、離散的フーリエ変換（ＤＦＴ）によって得られる。すなわち、
【０２５０】
【数９３】
【０２５１】
最後に、雑音が除去された周波数フレームが、周波数濾波によって、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算することによって得られる。すなわち、
【０２５２】
【数９４】
【０２５３】
この周波数濾波ステップは、時間的濾波によって、すなわち、ｘw（ｋ,ｎ）を
【０２５４】
【数９５】
【０２５５】
によって直接濾波することによって、等価に実行することもできる。この場合、時間領域で信号
【０２５６】
【数９６】
【０２５７】
を直接得る。
【０２５８】
したがって、濾波を周波数領域で実行する場合、時間領域に戻るのが望ましい。これは、逆スペクトル変換（ＩＤＦＴ）を適用することによって、装置１の機能ブロック１１によって実行される。すなわち、
【０２５９】
【数９７】
【０２６０】
雑音が除去された出力信号はこの後、ブロック・オーバーラップおよびアッド(block overlap and add)技術（ＯＬＡは「オーバーラップ・アッド」を意味する）によって、あるいは、ブロック・セーブ(block save)技術（ＯＬＳは「オーバーラップ・セーブ」を意味する）によって、有利なことに合成される。信号を時間領域において再構成するこの動作は、逆短期フーリエ変換（ＩＳＴＦＴ）と呼ばれる。
【０２６１】
これらのステップが完了すると、フレーム１２が装置１の出力に得られる。この出力フレーム１２は、歪みが抑えられている、入力フレーム２の、雑音が除去された推定値である。
【０２６２】
以降の説明では、本発明による具体化の有利な選択を説明する。これらの選択は代表的な実施形態であり、限定する性質を有すると見なされるべきではない。
【０２６３】
装置１の機能ブロック３および１１によって実行される分析と合成の具体化例に対しては、ＯＬＡ技術が選択されるのが有利である。さらに、分析されるフレーム同士の間にオーバーラップが存在する。たとえば、５０％のオーバーラップ率を選択してもよい。長さＬの、雑音のある信号ｘ（ｋ,ｎ）｛ｎ＝０，．．．．．．，Ｌ-１｝の現在のフレーム２はまず、信号ｘw（ｋ、ｎ）を送出する重みづけウインドウｗ（n）と乗算される。すなわち、
【０２６４】
【数９８】
【０２６５】
この例では、使用されたウインドウはサイズＬのハニング(Hanning)ウインドウである。すなわち、
【０２６６】
【数９９】
【０２６７】
ここで、Ｌ＝２５６であり、これは、サンプリング周波数Ｆｅ＝８kHzヘルツにおいて、長さ３２ｍｓのフレームに対応する。
【０２６８】
周波数領域への切り替えは、離散的フーリエ変換（ＤＦＴ）の数値的に最適化されたバージョンである高速フーリエ変換（ＦＦＴ）を使用して有利に実行される。より良好な周波数分解能を得るために、また、線形の畳み込み制約条件を満たすために、長さＬFFT＝５１２の高速フーリエ変換が使用される。したがって、ウインドウ化されたブロックｘw
（ｋ,ｎ）を２５６個の空サンプル（「ゼロパディング」）によって長くすることが、あらかじめ必要である。すなわち、
【０２６９】
【数１００】
【０２７０】
入力フレームのＦＦＴ Ｘ（ｋ,ｆ）が次に、次式
【０２７１】
【数１０１】
【０２７２】
によって得られる。
【０２７３】
ボーカルアクティビティの検出は、機能ブロック４によって、任意に実行することができる。
【０２７４】
機能ブロック５によって実行される、雑音電力スペクトル密度
【０２７５】
【数１０２】
【０２７６】
の推定に関する限り、次の方法で有利に進むことが可能である。ＤＶＡによって検出される「雑音のみ」の各フレームｋbに対して、雑音電力スペクトル密度
【０２７７】
【数１０３】
【０２７８】
が、次の再帰式
【０２７９】
【数１０４】
【０２８０】
によって推定される。
【０２８１】
現在のフレームが雑音フレームとして検出されなければ、雑音の電力スペクトル密度の推定は固定される。平滑化量（smoothing quantity）αは一定に選択され、α＝０．８８２５に等しい。これは、１２８ｍｓの時定数（５０％オーバーラップと、サンプリング周波数Ｆｅ＝８ｋHzに関係する）に対応し、信頼できる推定と雑音の統計の時間的変化の追跡との間の妥協を保証するに十分であると判断される。
【０２８２】
機能ブロック７によって推定された第１の雑音低減フィルタの伝達関数
【０２８３】
【数１０５】
【０２８４】
は、
【０２８５】
【数１０６】
【０２８６】
のように、有利なことに開ループ・ウイーナー・フィルタの伝達関数である。
【０２８７】
ここで、η（ｋ,ｆ）は、有用な信号
【０２８８】
【数１０７】
【０２８９】
のＰＳＤと雑音
【０２９０】
【数１０８】
【０２９１】
のＰＳＤとの推定値の比率によって理論的に定められた事前の信号対雑音比（ＳＮＲ）を表す。すなわち、
【０２９２】
【数１０９】
【０２９３】
有用な信号
【０２９４】
【数１１０】
【０２９５】
のＰＳＤは、有向決定推定値（directed-decision estimator）、すなわち、
【０２９６】
【数１１１】
【０２９７】
を使用して有利に得られる。ここで、重心パラメータβは一定に選択され、０．９８に等しい。上述したように、関数Ｐは、（１−β）によって重みづけされた量の閾値化（thresholding）を保証する。すなわち、
【０２９８】
【数１１２】
【０２９９】
伝達関数
【０３００】
【数１１３】
【０３０１】
はその後、有用な信号
【０３０２】
【数１１４】
【０３０３】
のＰＳＤの推定値を改良するために、有利に再使用される。したがって、
【０３０４】
【数１１５】
【０３０５】
によって量
【０３０６】
【数１１６】
【０３０７】
を得、これは、
【０３０８】
【数１１７】
【０３０９】
によって与えられる事前の信号対雑音比η2（ｋ,ｆ）の第２の推定を行うために使用される。
【０３１０】
次に、第２のパスが、η2（ｋ,ｆ）にもとづいて第１の雑音低減フィルタの伝達関数
【０３１１】
【数１１８】
【０３１２】
の推定値の計算に有利に存在する。すなわち、
【０３１３】
【数１１９】
【０３１４】
機能ブロック６によって実行される濾波は次に、以下に示すように実行される。雑音低減フィルタ
【０３１５】
【数１２０】
【０３１６】
の（周波数）伝達関数が利用可能であるので、対応する時間応答
【０３１７】
【数１２１】
【０３１８】
は、逆離散的フーリエ変換（ＩＤＦＴ）の数値的に最適化されたバージョンである逆高速フーリエ変換（ＩＦＦＴ）を使用することによって、有利に得られる。すなわち、
【０３１９】
【数１２２】
【０３２０】
信号
【０３２１】
【数１２３】
【０３２２】
が実数であるので、対応する時間フィルタがまず原因である(causal)とされる。意図された用途にとって重要なサンプルに対応するこのフィルタＬfilt1＝２５６係数がその後選択される。その後に、このようにして得られたインパルス応答は、長さＬfilt1のハニングウインドウによって重み付けされる。最後に、重み付けされたインパルス応答は、エイリアシング無しの線形の畳み込み制約条件を満たすインパルス応答
【０３２３】
【数１２４】
【０３２４】
を与えるために、ＬFFT−Ｌfilt1個の零で補充される。
【０３２５】
雑音低減フィルタの伝達関数
【０３２６】
【数１２５】
【０３２７】
はその後、離散的フーリエ変換（ＤＦＴ）の数値的に最適化されたバージョンである高速フーリエ変換（ＦＦＴ）を有利に使用することで得られる。すなわち、
【０３２８】
【数１２６】
【０３２９】
最後に、雑音のある信号の第１の、雑音が除去されたバージョンは、周波数濾波、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算することによって得られる。すなわち、
【０３３０】
【数１２７】
【０３３１】
機能ブロック８によって実行される、失われた高調波を回復するステップにおいて、信号
【０３３２】
【数１２８】
【０３３３】
が疑似信号
【０３３４】
【数１２９】
【０３３５】
を生成するためにその後に使用される。まず、時間的信号
【０３３６】
【数１３０】
【０３３７】
が、逆高速フーリエ変換（ＦＦＴＩ）を有利に使用することによって得られる。すなわち、
【０３３８】
【数１３１】
【０３３９】
その後、単一波の非線形修正関数が、信号に高調波性を回復させる目的で、この信号に適用される。すなわち、
【０３４０】
【数１３２】
【０３４１】
勿論、単一波の修正以外の（絶対値、閾値に対して最大／最小、等）任意の他の非線形性を、上に示したように、使用することができる。最後に、信号
【０３４２】
【数１３３】
【０３４３】
が、高速フーリエ変換（ＦＦＴ）を有利に使用することによって得られる。すなわち、
【０３４４】
【数１３４】
【０３４５】
この信号は、第２の雑音低減フィルタの伝達関数を有利に計算するために、その後、使用される。
【０３４６】
この代表的な実施形態において、このような第２のフィルタの推定は、装置１の機能ブロック９によって、次のように実行される。第２の雑音低減フィルタの伝達関数
【０３４７】
【数１３５】
【０３４８】
は、
【０３４９】
【数１３６】
【０３５０】
であるように、開ループ・ウイーナー・フィルタの伝達関数である。ここで、ηharmo（ｋ,ｆ）は、有用な信号
【０３５１】
【数１３７】
【０３５２】
のＰＳＤの推定値と雑音
【０３５３】
【数１３８】
【０３５４】
のＰＳＤの推定値との比率によって理論的に定められた事前の信号対雑音比（ＳＮＲ）を表す。すなわち、
【０３５５】
【数１３９】
【０３５６】
ここで、有用な信号のＰＳＤ、
【０３５７】
【数１４０】
【０３５８】
は
【０３５９】
【数１４１】
【０３６０】
のように得られる。
【０３６１】
信号
【０３６２】
【数１４２】
【０３６３】
の再投入のレベルを制御するパラメータρ（ｋ,ｆ）は、一定として選択され、意図された用途に対して０．５に等しい。
【０３６４】
雑音低減フィルタ
【０３６５】
【数１４３】
【０３６６】
の（周波数）伝達関数が利用可能であるので、対応する時間応答
【０３６７】
【数１４４】
【０３６８】
は、逆離散的フーリエ変換（ＩＤＦＴ）の数値的に最適化されたバージョンである逆高速フーリエ変換（ＦＦＴＩ）を有利に使用することによって得られる。すなわち、
【０３６９】
【数１４５】
【０３７０】
信号
【０３７１】
【数１４６】
【０３７２】
が実数であるので、対応する時間フィルタはまず原因であるとされる。この用途に対して重要なサンプルに対応するこのフィルタのＬfilt2＝２５６個の係数は、その後に選択される。このようにして得られたインパルス応答はその後長さＬfilt2のハニングウインドウによって重み付けされる。最後に、この重み付けされたインパルス応答は、エイリアシング無し線形の畳み込み制約条件を満たすインパルス応答
【０３７３】
【数１４７】
【０３７４】
を与えるために、ＬFFT＝Ｌfilt2個の零で補充される。
【０３７５】
雑音低減フィルタの伝達関数
【０３７６】
【数１４８】
【０３７７】
はその後に、離散的フーリエ変換（ＤＦＴ）の数値的に最適化されたバージョンである高速フーリエ変換（ＦＦＴ）を有利に使用することによって得られる。すなわち、
【０３７８】
【数１４９】
【０３７９】
雑音低減フィルタ
【０３８０】
【数１５０】
【０３８１】
の係数の値は、雑音除去が強力過ぎることを回避するために、任意に限界値を設定する（thresheld）ことが可能である。
【０３８２】
最後に、雑音が除去された周波数フレームは、機能ブロック１０によって実行される周波数濾波によって、すなわち、雑音のある信号のスペクトルと雑音低減フィルタの伝達関数を乗算することによって得られる。すなわち、
【０３８３】
【数１５１】
【０３８４】
このようにして得られた信号の時間領域への戻りは、逆高速フーリエ変換（ＦＦＴＩ）を有利に使用することによって、逆スペクトル変換によって、機能ブロック１１のレベルにおいて、その後に実行される。すなわち、
【０３８５】
【数１５２】
【０３８６】
雑音が除去された出力信号はその後、ブロック・オーバーラップおよびアッド技術ＯＬＡによって合成される。
【０３８７】
さらに、図２に示す装置１は、１つの音声ピックアップ装置、１つの通信装置、あるいは１つの音声認識装置のような、意図される用途の機能として、特定の装置の中に設置してもよいことに留意されるであろう。
【０３８８】
さらに、本装置によって実行される処理は、たとえば通信端末の場合のようにローカルに、あるいは、たとえば音声認識サーバの場合にネットワーク内に集中化して実行してもよいことが留意される。
【０３８９】
さらに、実際には、本発明による雑音のある音声信号を処理する装置１は、本発明による信号を処理する方法を実施するのに適した（電子的）ハードウェア手段および／またはソフトウェア手段から成る。
【０３９０】
好ましい具体化例によれば、本発明による雑音のある音声信号を処理する方法のステップは、本発明によるこのような装置の１つで使用されるコンピュータ・プログラムの命令によって決定される。
【０３９１】
本発明による方法は、前述のプログラムが、装置に組み込まれた計算手段にロードされたときに実行され、次にその動作はプログラムの実行によって制御される。
【０３９２】
ここで、表現「コンピュータ・プログラム」は、（ソフトウェア）集合を形成する１つまたは複数のコンピュータ・プログラムを意味すると理解される。ソフトウェア集合の目的は、適切な計算手段によって実行される時に、本発明を実施することである。
【０３９３】
したがって、本発明の主題も、このようなコンピュータ・プログラムであり、特に、情報媒体に格納された１つのソフトウェアの形式のコンピュータ・プログラムである。このような情報媒体は、本発明によるプログラムを格納することが可能な任意の実体あるいは装置からなってよい。
【０３９４】
たとえば、当該媒体は、ＲＯＭ、たとえばＣＤＲＯＭ、あるいはマイクロ電子回路ＲＯＭ、あるいは磁気記録手段、たとえばハードディスク、のようなハードウェア記憶手段を有してもよい。変形例として、情報媒体は、プログラムが組み込まれ、当該方法の実行に適した、あるいは当該方法の実行に使用されることに適した集積回路であってもよい。
【０３９５】
さらに、情報媒体は、無線によってあるいは他の手段によって、電気ケーブルあるいは光ケーブルによって伝送可能な電気信号あるいは光信号のような伝送可能な無形の媒体であってもよい。本発明によるプログラムは、特に、インターネット形式のネットワークからダウンロードすることができる。
【０３９６】
設計の視点から、本発明によるコンピュータ・プログラムは、任意のプログラミング言語を使用することが可能であり、ソース・コード、オブジェクトコード、あるいは、ソース・コードとオブジェクトコードの間の中間コード（たとえば部分的にコンパイルされた形式）の形式、あるいは、本発明による方法を実施するのに望ましい任意の他の形式とすることが可能である。
【図面の簡単な説明】
【０３９７】
【図１】すでに説明済みであるが、３つの公知の短期間の抑圧規則の挙動を比較する図である。
【図２】本発明による実施形態を実行することが可能な、さまざまな機能ブロックを示す図である。
【図３】雑音が除去された信号内の失われた高調波の回復の例を示す図である。
【図４】本発明の一実施形態による信号内の高調波の回復の効果を示す図である。
【図５】本発明の一実施形態による第２の雑音低減フィルタの推定の例を示す図である。
【特許請求の範囲】
【請求項１】
連続した複数のフレームとして構成された、雑音のある音声信号（ｘ（ｎ））を処理する方法であって、前記フレーム（ｘ（ｋ,ｎ））の少なくとも１つに関して、
前記雑音のある音声信号の前記フレームに、周波数領域への変換を適用するステップと、
前記フレームに対する前記雑音の電力スペクトル密度
【数１】
を推定するステップと、
前記雑音の前記推定された電力スペクトル密度と、前記フレームに対応する有用な信号の前記電力スペクトル密度の推定値とにもとづいて、第１の雑音低減フィルタ
【数２】
を計算するステップと、
前記フレームの第１の、雑音が除去された推定値
【数３】
を得るために、計算された前記第１の雑音低減フィルタを利用して、前記雑音のある音声信号（ｘ(ｋ,ｎ);Ｘ(ｋ、ｆ)）の前記フレームを濾波するステップと、
前記雑音のある音声信号の前記フレームに対応する前記有用な信号と実質的に同じ位置にある高調波を有する、第２の信号
【数４】
のフレームを、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値にもとづいて得るステップと、
を有する、連続したフレームとして構成された雑音のある音声信号を処理する方法。
【請求項２】
前記雑音の前記推定された電力スペクトル密度と、前記フレームの第１の、雑音が除去された推定値
【数５】
の電力と得られた前記第２の信号
【数６】
の前記フレームの電力の組み合わせとにもとづいて、第２の雑音低減フィルタ
【数７】
を計算するステップと、
前記フレームの第２の、雑音が除去された
【数８】
推定値を得るために、計算された前記第２の雑音低減フィルタを利用して、前記雑音のある音声信号（ｘ(ｋ,ｎ);Ｘ(ｋ、ｆ)）の前記フレームを濾波するステップと、
前記フレームの前記第２の雑音が除去された推定値を合成するステップと
をさらに有する、請求項１に記載の方法。
【請求項３】
前記第２の雑音低減フィルタを利用して実行された前記濾波の結果が前記周波数領域にある場合、前記フレームの前記第２の、雑音が除去された推定値を合成する前に、前記濾波の前記結果に前記時間領域への変換を適用する、請求項２に記載の方法。
【請求項４】
前記雑音のある音声信号（ｘ（ｋ,ｎ））の前記フレームに、それに対して前記周波数領域への変換を適用する前に、重みづけウインドウ（ｗ（ｎ））を乗算する、請求項１から請求項３までのいずれか１項に記載の方法。
【請求項５】
前記重みづけウインドウは、サイズが前記雑音のある音声信号の前記フレームの長さＬであるハニングウインドウである、請求項４に記載の方法。
【請求項６】
前記重みづけウインドウは、ｗ（ｎ）＝０．５−０．５．ｃｏｓ（πｎ／Ｌ）の形式であり、ここで、ｎは０からＬ−１の範囲の整数であり、Ｌは２５６に等しい、請求項５に記載の方法。
【請求項７】
前記雑音のある音声信号（ｘ（ｋ，ｎ））の前記フレームに適用される、前記周波数領域への前記変換は、長さ５１２の高速フーリエ変換である、請求項１から請求項６までのいずれか１項に記載の方法。
【請求項８】
ボーカルアクティビティを前記雑音のある音声信号（ｘ（ｋ，ｎ））の前記フレーム内で検索し、前記フレームに対する前記雑音の前記電力スペクトル密度
【数９】
の推定値は、前記フレーム内でボーカルアクティビティが検出されない場合には、前記フレームに含まれる前記雑音の前記電力スペクトル密度の推定値に対応し、前記フレーム内でボーカルアクティビティが検出された場合には、前記フレームに先行する、前記雑音のある音声信号の少なくとも１つのフレームに含まれる、前記雑音の前記電力スペクトル密度の推定値に対応する、請求項１から請求項７までのいずれか１項に記載の方法。
【請求項９】
前記雑音のある音声信号のインデックスｋbの、ボーカルアクティビティが検出されない少なくとも１つのフレームに含まれる前記雑音の前記電力スペクトル密度を、
【数１０】
の形式の指数関数による平滑化によって求め、ここで、Ｘ（ｋb，ｆ）は前記雑音のある音声信号のインデックスｋbのフレームの前記周波数領域への変換を表し、αは平滑化量である、請求項８に記載の方法。
【請求項１０】
αは実質的に０．８８２５に等しい、請求項９に記載の方法。
【請求項１１】
前記第１の雑音低減フィルタの前記計算は、短期スペクトル減衰の技術を実行する第１のパスを有する、請求項１から請求項１０までのいずれか１項に記載の方法。
【請求項１２】
前記第１の雑音低減フィルタの前記計算は、
【数１１】
の形式の伝達関数を有する開ループ・ウイーナー・フィルタを実行する第１のパスを含み、ここで、ｋは前記雑音のある音声信号の前記フレームのインデックスであり、η（ｋ，ｆ）は前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値と前記雑音の前記推定電力スペクトル密度の比率を表す、請求項１１に記載の方法。
【請求項１３】
前記フレームに対応する有用な信号の前記電力スペクトル密度の前記推定値は、形式
【数１２】
の有向決定推定値を有し、ここで、βは重心パラメータ、あるいは指数関数平滑化推定値である、請求項１２に記載の方法。
【請求項１４】
前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値は、実質的に０．９８に等しい重心パラメータβを使用する有向決定推定値を含む、請求項１３に記載の方法。
【請求項１５】
前記第１の雑音低減フィルタの前記計算は短期スペクトル減衰の技術を実行する第２のパスをさらに含み、前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値は前記第１のパスの間に実行される前記計算を考慮している、請求項１１から１４のいずれか１項に記載の方法。
【請求項１６】
前記第１のパスの間に実行される前記計算を考慮した、前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値は
【数１３】
の形式であって、ここで、
【数１４】
は前記第１のパスの間に計算された伝達関数を表し、Ｘ（ｋ，ｆ）は前記雑音のある音声信号のインデックスｋの前記フレームの前記周波数領域への変換を表し、また、前記第１の雑音低減フィルタの前記計算は形式
【数１５】
の伝達関数を有する開ループ・ウイーナー・フィルタを実行する第２のパスを有し、ここで、η2（ｋ，ｆ）は前記第１のパスの間に計算された前記計算を考慮した前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値と前記雑音の前記推定電力スペクトル密度の比率を表す、請求項１５に記載の方法。
【請求項１７】
前記第１の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波を、前記フレーム（Ｘ（ｋ，ｆ））の周波数領域への変換で実行する、請求項１から請求項１６までのいずれか１項に記載の方法。
【請求項１８】
前記第１の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波は、前記第１の雑音低減フィルタの係数の選択と、前記第１の雑音低減フィルタの前記時間的インパルス応答の重みづけを含む、請求項１から請求項１７までのいずれか１項に記載の方法。
【請求項１９】
前記第１の雑音低減フィルタの２５６個の係数を、ハニングウインドウによって選択し、重み付けする、請求項１８に記載の方法。
【請求項２０】
前記雑音のある音声信号の前記フレームに適用される前記周波数領域への変換は、長さが５１２の高速フーリエ変換であり、前記第１の雑音低減フィルタ（ｘw（ｋ，ｎ））の前記時間的インパルス応答を、２５６個の零で補充する、請求項１９に記載の方法。
【請求項２１】
前記雑音のある音声信号の前記フレームに対応する前記有用な信号と実質的に同じ位置にある高調波を有する第２の信号
【数１６】
のフレームを得ることは、前記第１の推定値が前記時間領域内にある場合には、前記雑音のある音声信号
【数１７】
の前記フレームの前記第１の、雑音が除去された推定値への非線形の関数の適用を含み、前記第１の推定値が前記周波数領域内にある場合には、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値
【数１８】
と、非線形の関数（ｄ（ｋ，ｎ））の前記周波数領域への変換の前記第１の、雑音が除去された推定値との間の円形の畳み込みの適用を含む、請求項１から請求項２０までのいずれか１項に記載の方法。
【請求項２２】
前記非線形の関数（ｄ（ｋ，ｎ））は、単一波の修正関数、絶対値、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値と閾値との間の最大値、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値と閾値との間の最小値、の中の１つである、請求項２１記載の方法。
【請求項２３】
前記第２の雑音低減フィルタ
【数１９】
の前記計算は、短期スペクトル減衰の技術を実行する、請求項２から請求項２２までのいずれか１項に記載の方法。
【請求項２４】
前記第２の雑音低減フィルタの前記計算は、
【数２０】
および
【数２１】
である
【数２２】
の伝達関数を有する開ループ・ウイーナー・フィルタを実行し、ここで、
【数２３】
はインデックスｋの前記フレームの前記第１の、雑音が除去された推定値を表し、
【数２４】
は前記第２の信号の前記フレームを表し、ρ（ｋ，ｆ）は前記第２の信号に対する再投入パラメータである、
請求項２３に記載の方法。
【請求項２５】
ρ（ｋ，ｆ）は、周波数および／または時間の関数として変化する、請求項２４に記載の方法。
【請求項２６】
ρ（ｋ，ｆ）は実質的に一定であり、０．５に等しい、請求項２４に記載の方法。
【請求項２７】
前記第２の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波は、前記フレーム（Ｘ（ｋ，ｆ））の周波数領域への変換で実行される、請求項１から請求項２６までのいずれか１項に記載の方法。
【請求項２８】
前記第２の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波は、前記第２の雑音低減フィルタの時間的インパルス応答の選択と重みづけを含む、請求項２から請求項２７までのいずれか１項に記載の方法。
【請求項２９】
前記第２の雑音低減フィルタの２５６個の係数をハニングウインドウによって選択し、重み付けする、請求項２８記載の方法。
【請求項３０】
前記雑音のある音声信号の前記フレームに適用される前記周波数領域への前記変換は、長さが５１２の高速フーリエ変換であり、前記第２の雑音低減フィルタの前記時間的インパルス応答を２５６個の零で補充する、請求項２９記載の方法。
【請求項３１】
前記フレームの第２の、雑音が除去された推定値の前記合成は、ＯＬＡブロック・オーバーラップおよびアッド、あるいは、ＯＬＳブロック・セーブを使用する、請求項２から請求項３０までのいずれか１項に記載の方法。
【請求項３２】
請求項１から請求項３１までのいずれか１項に記載の方法を実行するように設計された手段を有する、雑音のある音声信号を処理する装置（１）。
【請求項３３】
情報媒体上のコンピュータ・プログラムにおいて、前記プログラムが計算手段によってロードされ、実行される場合に、請求項１から請求項３１までのいずれか１項に記載の方法を実行するように作られた命令を有することを特徴とする、情報媒体上のコンピュータ・プログラム。
【請求項１】
連続した複数のフレームとして構成された、雑音のある音声信号（ｘ（ｎ））を処理する方法であって、前記フレーム（ｘ（ｋ,ｎ））の少なくとも１つに関して、
前記雑音のある音声信号の前記フレームに、周波数領域への変換を適用するステップと、
前記フレームに対する前記雑音の電力スペクトル密度
【数１】
を推定するステップと、
前記雑音の前記推定された電力スペクトル密度と、前記フレームに対応する有用な信号の前記電力スペクトル密度の推定値とにもとづいて、第１の雑音低減フィルタ
【数２】
を計算するステップと、
前記フレームの第１の、雑音が除去された推定値
【数３】
を得るために、計算された前記第１の雑音低減フィルタを利用して、前記雑音のある音声信号（ｘ(ｋ,ｎ);Ｘ(ｋ、ｆ)）の前記フレームを濾波するステップと、
前記雑音のある音声信号の前記フレームに対応する前記有用な信号と実質的に同じ位置にある高調波を有する、第２の信号
【数４】
のフレームを、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値にもとづいて得るステップと、
を有する、連続したフレームとして構成された雑音のある音声信号を処理する方法。
【請求項２】
前記雑音の前記推定された電力スペクトル密度と、前記フレームの第１の、雑音が除去された推定値
【数５】
の電力と得られた前記第２の信号
【数６】
の前記フレームの電力の組み合わせとにもとづいて、第２の雑音低減フィルタ
【数７】
を計算するステップと、
前記フレームの第２の、雑音が除去された
【数８】
推定値を得るために、計算された前記第２の雑音低減フィルタを利用して、前記雑音のある音声信号（ｘ(ｋ,ｎ);Ｘ(ｋ、ｆ)）の前記フレームを濾波するステップと、
前記フレームの前記第２の雑音が除去された推定値を合成するステップと
をさらに有する、請求項１に記載の方法。
【請求項３】
前記第２の雑音低減フィルタを利用して実行された前記濾波の結果が前記周波数領域にある場合、前記フレームの前記第２の、雑音が除去された推定値を合成する前に、前記濾波の前記結果に前記時間領域への変換を適用する、請求項２に記載の方法。
【請求項４】
前記雑音のある音声信号（ｘ（ｋ,ｎ））の前記フレームに、それに対して前記周波数領域への変換を適用する前に、重みづけウインドウ（ｗ（ｎ））を乗算する、請求項１から請求項３までのいずれか１項に記載の方法。
【請求項５】
前記重みづけウインドウは、サイズが前記雑音のある音声信号の前記フレームの長さＬであるハニングウインドウである、請求項４に記載の方法。
【請求項６】
前記重みづけウインドウは、ｗ（ｎ）＝０．５−０．５．ｃｏｓ（πｎ／Ｌ）の形式であり、ここで、ｎは０からＬ−１の範囲の整数であり、Ｌは２５６に等しい、請求項５に記載の方法。
【請求項７】
前記雑音のある音声信号（ｘ（ｋ，ｎ））の前記フレームに適用される、前記周波数領域への前記変換は、長さ５１２の高速フーリエ変換である、請求項１から請求項６までのいずれか１項に記載の方法。
【請求項８】
ボーカルアクティビティを前記雑音のある音声信号（ｘ（ｋ，ｎ））の前記フレーム内で検索し、前記フレームに対する前記雑音の前記電力スペクトル密度
【数９】
の推定値は、前記フレーム内でボーカルアクティビティが検出されない場合には、前記フレームに含まれる前記雑音の前記電力スペクトル密度の推定値に対応し、前記フレーム内でボーカルアクティビティが検出された場合には、前記フレームに先行する、前記雑音のある音声信号の少なくとも１つのフレームに含まれる、前記雑音の前記電力スペクトル密度の推定値に対応する、請求項１から請求項７までのいずれか１項に記載の方法。
【請求項９】
前記雑音のある音声信号のインデックスｋbの、ボーカルアクティビティが検出されない少なくとも１つのフレームに含まれる前記雑音の前記電力スペクトル密度を、
【数１０】
の形式の指数関数による平滑化によって求め、ここで、Ｘ（ｋb，ｆ）は前記雑音のある音声信号のインデックスｋbのフレームの前記周波数領域への変換を表し、αは平滑化量である、請求項８に記載の方法。
【請求項１０】
αは実質的に０．８８２５に等しい、請求項９に記載の方法。
【請求項１１】
前記第１の雑音低減フィルタの前記計算は、短期スペクトル減衰の技術を実行する第１のパスを有する、請求項１から請求項１０までのいずれか１項に記載の方法。
【請求項１２】
前記第１の雑音低減フィルタの前記計算は、
【数１１】
の形式の伝達関数を有する開ループ・ウイーナー・フィルタを実行する第１のパスを含み、ここで、ｋは前記雑音のある音声信号の前記フレームのインデックスであり、η（ｋ，ｆ）は前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値と前記雑音の前記推定電力スペクトル密度の比率を表す、請求項１１に記載の方法。
【請求項１３】
前記フレームに対応する有用な信号の前記電力スペクトル密度の前記推定値は、形式
【数１２】
の有向決定推定値を有し、ここで、βは重心パラメータ、あるいは指数関数平滑化推定値である、請求項１２に記載の方法。
【請求項１４】
前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値は、実質的に０．９８に等しい重心パラメータβを使用する有向決定推定値を含む、請求項１３に記載の方法。
【請求項１５】
前記第１の雑音低減フィルタの前記計算は短期スペクトル減衰の技術を実行する第２のパスをさらに含み、前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値は前記第１のパスの間に実行される前記計算を考慮している、請求項１１から１４のいずれか１項に記載の方法。
【請求項１６】
前記第１のパスの間に実行される前記計算を考慮した、前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値は
【数１３】
の形式であって、ここで、
【数１４】
は前記第１のパスの間に計算された伝達関数を表し、Ｘ（ｋ，ｆ）は前記雑音のある音声信号のインデックスｋの前記フレームの前記周波数領域への変換を表し、また、前記第１の雑音低減フィルタの前記計算は形式
【数１５】
の伝達関数を有する開ループ・ウイーナー・フィルタを実行する第２のパスを有し、ここで、η2（ｋ，ｆ）は前記第１のパスの間に計算された前記計算を考慮した前記フレームに対応する前記有用な信号の前記電力スペクトル密度の前記推定値と前記雑音の前記推定電力スペクトル密度の比率を表す、請求項１５に記載の方法。
【請求項１７】
前記第１の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波を、前記フレーム（Ｘ（ｋ，ｆ））の周波数領域への変換で実行する、請求項１から請求項１６までのいずれか１項に記載の方法。
【請求項１８】
前記第１の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波は、前記第１の雑音低減フィルタの係数の選択と、前記第１の雑音低減フィルタの前記時間的インパルス応答の重みづけを含む、請求項１から請求項１７までのいずれか１項に記載の方法。
【請求項１９】
前記第１の雑音低減フィルタの２５６個の係数を、ハニングウインドウによって選択し、重み付けする、請求項１８に記載の方法。
【請求項２０】
前記雑音のある音声信号の前記フレームに適用される前記周波数領域への変換は、長さが５１２の高速フーリエ変換であり、前記第１の雑音低減フィルタ（ｘw（ｋ，ｎ））の前記時間的インパルス応答を、２５６個の零で補充する、請求項１９に記載の方法。
【請求項２１】
前記雑音のある音声信号の前記フレームに対応する前記有用な信号と実質的に同じ位置にある高調波を有する第２の信号
【数１６】
のフレームを得ることは、前記第１の推定値が前記時間領域内にある場合には、前記雑音のある音声信号
【数１７】
の前記フレームの前記第１の、雑音が除去された推定値への非線形の関数の適用を含み、前記第１の推定値が前記周波数領域内にある場合には、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値
【数１８】
と、非線形の関数（ｄ（ｋ，ｎ））の前記周波数領域への変換の前記第１の、雑音が除去された推定値との間の円形の畳み込みの適用を含む、請求項１から請求項２０までのいずれか１項に記載の方法。
【請求項２２】
前記非線形の関数（ｄ（ｋ，ｎ））は、単一波の修正関数、絶対値、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値と閾値との間の最大値、前記雑音のある音声信号の前記フレームの前記第１の、雑音が除去された推定値と閾値との間の最小値、の中の１つである、請求項２１記載の方法。
【請求項２３】
前記第２の雑音低減フィルタ
【数１９】
の前記計算は、短期スペクトル減衰の技術を実行する、請求項２から請求項２２までのいずれか１項に記載の方法。
【請求項２４】
前記第２の雑音低減フィルタの前記計算は、
【数２０】
および
【数２１】
である
【数２２】
の伝達関数を有する開ループ・ウイーナー・フィルタを実行し、ここで、
【数２３】
はインデックスｋの前記フレームの前記第１の、雑音が除去された推定値を表し、
【数２４】
は前記第２の信号の前記フレームを表し、ρ（ｋ，ｆ）は前記第２の信号に対する再投入パラメータである、
請求項２３に記載の方法。
【請求項２５】
ρ（ｋ，ｆ）は、周波数および／または時間の関数として変化する、請求項２４に記載の方法。
【請求項２６】
ρ（ｋ，ｆ）は実質的に一定であり、０．５に等しい、請求項２４に記載の方法。
【請求項２７】
前記第２の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波は、前記フレーム（Ｘ（ｋ，ｆ））の周波数領域への変換で実行される、請求項１から請求項２６までのいずれか１項に記載の方法。
【請求項２８】
前記第２の雑音低減フィルタを利用した、前記雑音のある音声信号（ｘ（ｋ，ｎ）；Ｘ（ｋ，ｆ））の前記フレームの前記濾波は、前記第２の雑音低減フィルタの時間的インパルス応答の選択と重みづけを含む、請求項２から請求項２７までのいずれか１項に記載の方法。
【請求項２９】
前記第２の雑音低減フィルタの２５６個の係数をハニングウインドウによって選択し、重み付けする、請求項２８記載の方法。
【請求項３０】
前記雑音のある音声信号の前記フレームに適用される前記周波数領域への前記変換は、長さが５１２の高速フーリエ変換であり、前記第２の雑音低減フィルタの前記時間的インパルス応答を２５６個の零で補充する、請求項２９記載の方法。
【請求項３１】
前記フレームの第２の、雑音が除去された推定値の前記合成は、ＯＬＡブロック・オーバーラップおよびアッド、あるいは、ＯＬＳブロック・セーブを使用する、請求項２から請求項３０までのいずれか１項に記載の方法。
【請求項３２】
請求項１から請求項３１までのいずれか１項に記載の方法を実行するように設計された手段を有する、雑音のある音声信号を処理する装置（１）。
【請求項３３】
情報媒体上のコンピュータ・プログラムにおいて、前記プログラムが計算手段によってロードされ、実行される場合に、請求項１から請求項３１までのいずれか１項に記載の方法を実行するように作られた命令を有することを特徴とする、情報媒体上のコンピュータ・プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公表番号】特表２００８−５１３８１９（Ｐ２００８−５１３８１９Ａ）
【公表日】平成２０年５月１日（２００８．５．１）
【国際特許分類】
【出願番号】特願２００７−５３１７９７（Ｐ２００７−５３１７９７）
【出願日】平成１７年９月１４日（２００５．９．１４）
【国際出願番号】ＰＣＴ／ＦＲ２００５／００２２８４
【国際公開番号】ＷＯ２００６／０３２７６０
【国際公開日】平成１８年３月３０日（２００６．３．３０）
【出願人】（５９１０３４１５４）フランス テレコム (290)
【公表日】平成２０年５月１日（２００８．５．１）
【国際特許分類】
【出願日】平成１７年９月１４日（２００５．９．１４）
【国際出願番号】ＰＣＴ／ＦＲ２００５／００２２８４
【国際公開番号】ＷＯ２００６／０３２７６０
【国際公開日】平成１８年３月３０日（２００６．３．３０）
【出願人】（５９１０３４１５４）フランス テレコム (290)
[ Back to top ]