JP5769670B2 - エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置とプログラム - Google Patents

Description

この発明は、音響再生系を有する通信会議システム等に用いられるエコー消去装置及び、それに適用されるエコー抑圧ゲイン推定方法と、そのプログラムに関する。

短時間スペクトル振幅（STSA：Sort-Time Spectral Amplitude）推定に基づくエコー抑圧処理は、エコーと近端話者音声の無相関を仮定してエコーを抑圧するゲイン係数を推定し、振幅周波数領域でエコーを抑圧する方法であり、例えば非特許文献１に記載されている。これに対し、エコー抑圧後の音声品質向上のため、非特許文献１に記載される上記ゲイン係数の推定方法（以下、「エコー抑圧ゲイン推定方法」と称す）を改善した手法として、特許文献１が提案されている。ここでは、特許文献１に開示されたエコー抑圧ゲイン推定方法を用いたエコー消去装置１０の機能構成例を、図７に示してその動作を簡単に説明する。

エコー消去装置１０は、再生信号周波数分析部８１と、収音信号周波数分析部８２と、エコーパワー推定部８３と、類似度係数計算部１２と、エコー抑圧ゲイン計算部１４と、乗算部８５と、周波数合成部８６とを備える。再生信号ｘ（ｋ）は、例えばサンプリング周波数１６ｋＨｚで離散値とされた信号であり、スピーカ１によって音響信号に変換される。なお、再生信号を離散値化するＡＤ変換器と、その離散値を連続値に変換するＤＡ変換器については省略している。

再生信号周波数分析部８１は、再生信号の離散値を２５６点集めて１フレームとし、１/２オーバーラップ加算による周波数分析を行い、フレーム単位で８ｋＨｚまでの周波数範囲を１２８個の再生信号スペクトルＸ_ωに変換する。ｋはフレーム番号であり、ωはこの例の場合６４Ｈｚの間隔で得られる周波数スペクトルの番号（１〜１２８）である。

マイクロホン２で収音される収音信号ｙ（ｋ）は、再生信号ｘ（ｋ）がスピーカ１から再生されたことによるエコーｄ（ｋ）が、近端話者信号ｓ（ｋ）に重畳した信号である。その収音信号ｙ（ｋ）は、再生信号ｘ（ｋ）と同じように収音信号周波数分析部８２によって収音信号スペクトルＹ_ωに変換される。Ｙ_ωはＹ_ω＝Ｄ_ω＋Ｓ_ω、Ｄ_ωはエコースペクトル、Ｓ_ωは近端話者信号スペクトル、である。

エコーパワー推定部８３は、再生信号スペクトルＸ_ωと収音信号スペクトルＹ_ωとを入力として式（１）に示すエコーパワー推定値｜Ｄ_ω＾｜^２を計算する。＾は推定値であることを表わすが、その表記は式及び図中に示すものが正しい。

ここで、Ｈ_ωはスピーカ１からマイクロホン２に回り込むエコー経路のスペクトル、｜Ｈ_ω＾｜^２は音響結合量の推定値であり、｜Ｈ_ω＾｜^２￣は１フレーム過去の音響結合量の推定値である。ｍｉｎ｛ ， ｝は最小値を選択する関数である。したがって、｜Ｈ_ω＾｜^２は隣接するフレーム間で小さい方の音響結合量がその推定値になる。

エコー抑圧ゲイン計算部１４は、収音信号スペクトルＹ_ωと、エコーパワー推定値｜Ｄ_ω＾｜^２と、類似度係数｜ｒ_ω｜を入力として式（２）に示すゲイン係数Ｇ_ωを出力する。ゲイン係数Ｇ_ωは０〜１の実数値をとり、収音信号スペクトルＹ_ω中にエコー成分が多い場合には小さな値、エコー成分以外の成分が多い場合には大きな値をとる。

乗算部８５は、収音信号スペクトルＹ_ωにゲイン係数Ｇ_ωを乗算する。エコー成分が多い場合のゲイン係数Ｇ_ωは小さな値となるので、乗算部８５の出力信号はエコー成分が抑圧された近端話者信号スペクトル推定値Ｓ_ω＾となる。この各周波数成分ωに対応する近端話者信号スペクトル推定値Ｓ_ω＾は、周波数合成部８６において時間領域の出力信号ｅ（ｋ）に再合成される。

類似度係数計算部１２は、再生信号スペクトルＸ_ωと収音信号スペクトルＹ_ωとを入力として、内積＜Ｘ_ω，Ｙ_ω＞を例えば式（３）で、再生信号スペクトルのノルム値‖Ｘ_ω‖と、収音信号スペクトルのノルム値‖Ｙ_ω‖を、例えば式（４）と式（５）でそれぞれ計算する。

ここで、＊は複素共役、￣は１フレーム過去であることを表わす。εは０＜ε≦１を満たす忘却係数であり、指数関数的な減衰の時定数を決定する。例えばε＝０.０１６とする。εが１に近づくほど現在の再生信号スペクトルＸ_ωと収音信号スペクトルＹ_ωに依存した（重み付けされた）それぞれの値になる。なお、周波数軸方向にも着目した式（６）〜（８）を用いてもよい。

ここでＭ₁〜Ｍ₂は所定の周波数範囲を表わす。式（６）〜（８）は周波数軸方向に相関を取った後に時間軸方向の相関を取る形である。相関を取る順番は逆でもよい。また、式（３）に替えて位相成分を考慮しない絶対値から求める式（９）を用いてもよい。

類似度係数計算部１２は、式（３）〜（５）でそれぞれの値を計算した後に、式（１０）で類似度係数｜ｒ_ω｜を計算する。

ここで類似度係数｜ｒ_ω｜を算出する式（１０）の導出根拠について説明する。
例えばウィナーフィルタ法（Wiener Filter以下「ＷＦ法」と称す）では、式（１１）の評価量εを最小とするゲイン係数Ｇ_ωを推定してエコーを抑圧する。

式（１１）をＧ_ωで微分すると次式で表せる。

式（１２）を０と置いた式からゲイン係数Ｇ_ωを求めると式（１３）となる。

式（１３）を変形するとゲイン係数Ｇ_ωは式（１４）で表せる。

ここでｒ_ωはエコースペクトルＤ_ωと近端話者信号スペクトルＳ_ωの複素コヒーレンスを示す。Ｄ_ω＝Ｈ_ω ^＊Ｘ_ωを用いてｒ_ωから未知のベクトルＤ_ωを消去すると式（１５）で表せる。

ここで＊は複素共役を表す。式（１５）から分かるようにＤ_ωを消去しても未知の変数であるＨ_ωが存在するため、ｒ_ωを求めることは困難である。しかし、その絶対値ならば式（１６）に示すようにＨ_ωが消去可能なことに着目する。

そしてゲイン係数を次式に示すように与える。

この式（１７）から明らかなようにゲイン係数Ｇ_ωは、類似度係数｜ｒ_ω｜が１に近づくほど小さな値になりエコーが抑圧される。と同時にこのようにして求めたゲイン係数は、近端話者信号ｓ（ｋ）の欠損を軽減するように動作する。つまり、＜Ｄ_ω，Ｓ_ω＞＝０となるほどの長時間を要することなく、短時間のスペクトル振幅を用いてエコーを抑圧するので、内積値が０でない分が誤差となって発生するミュージカルノイズの発生を抑圧することができる。

特許第４７８７８５１号

阪内澄宇、羽田陽一、片岡章俊著「ＳＴＳＡ推定に基づくエコー抑圧処理のゲイン強調化方式」信学論（Ａ）,１vol.J88-A,no.6,Jun.2005,p695-703

類似度係数｜ｒ_ω｜は、収音信号中のエコー成分の割合を表す値で０〜１の間の値をとる係数である。したがって、近端話者信号ｓ（ｋ）の欠損を軽減するように動作するが、類似度係数｜ｒ_ω｜は１より大きな値をとることができないため、エコーをより大きく抑圧するようには動作しない。その結果、エコーが十分に抑圧できない場合があり、抑圧し切れなかったエコー成分が音質劣化として聴感上知覚される場合がある。エコーが十分に抑圧できない場合とは、近端話者信号ｓ（ｋ）の相関が残ってエコーが真値より小さく見積もられてしまう場合である。また、近端話者信号ｓ（ｋ）の相関が残る現象は、近端話者信号ｓ（ｋ）が大きい場合に顕著になる。

この発明は、このような課題に鑑みてなされたものであり、エコー成分が残ることに起因する音質劣化を防止するエコー抑圧ゲイン推定方法と、それを用いたエコー消去装置とそのプログラムを提供することを目的とする。

この発明のエコー抑圧ゲイン推定方法は、再生信号周波数分析段階と、収音信号周波数分析段階と、エコーパワー推定段階と、エコーパワー補正係数計算段階と、エコー抑圧ゲイン計算段階と、を有する。再生信号周波数分析段階は、再生信号を周波数領域の再生信号スペクトルＸ_ωに変換する。収音信号周波数分析段階は、収音信号を周波数領域の収音信号スペクトルＹ_ωに変換する。エコーパワー推定段階は、再生信号スペクトルＸ_ωと収音信号スペクトルＹ_ωとを入力としてエコーパワー推定値｜Ｄ_ω＾｜^２を計算する。エコーパワー補正係数計算段階は、収音信号スペクトルＹ_ωとエコーパワー推定値｜Ｄ_ω＾｜^２とを入力として、エコーパワー推定値の平方根｜Ｄ_ω＾｜と収音信号スペクトルの大きさ｜Ｙ_ω｜との内積値を２分のα乗した値を、エコーパワー推定値の平方根のα乗ノルム値‖Ｄ_ω＾‖^αで除した値をエコーパワー補正係数推定値γ_ω＾として計算する。エコー抑圧ゲイン計算段階は、収音信号スペクトルＹ_ωと、エコーパワー補正係数推定値γ_ω＾と、エコーパワー推定値｜Ｄ_ω＾｜^２とを入力として、エコーを抑圧するゲイン係数を上記エコーパワー補正係数推定値γ_ω＾が大きな値をとるほど小さな値になるように計算する。

この発明のエコー抑圧ゲイン推定方法によれば、エコーパワー補正係数推定値が０〜無限大の範囲の値をとり、エコーパワー補正係数推定値が大きな値をとるほどゲイン係数が小さくなりエコーが抑圧され、逆にエコーパワー補正係数推定値が小さな値をとるほどゲイン係数が大きくなり近端話者信号の欠損を軽減するように動作する。したがって、従来法に比べてエコーを十分に抑圧できるように動作し、音声歪みに対してより頑健な音質劣化の少ないエコー抑圧を行うことが可能になる。

この発明のエコー抑圧ゲイン推定方法を用いたエコー消去装置１００の機能構成例を示す図。 この発明のエコー抑圧ゲイン推定方法の動作フローを示す図。 エコーパワー補正係数計算部１１０の機能構成例を示す図。 エコー抑圧ゲイン計算部１２０の機能構成例を示す図。 エコー抑圧ゲイン計算部１２０の動作フローを含めたエコー抑圧ゲイン推定方法の動作フローを示す図。 シミュレーション結果を示す図。 従来のエコー消去装置１０の機能構成を示す図。

以下、この発明の実施の形態を図面を参照して説明する。複数の図面中同一のものには同じ参照符号を付し、説明は繰り返さない。

〔発明の基本的な考え〕
実施例の説明の前に、この発明のエコー抑圧方法の基本的な考えを説明する。この発明は、エコースペクトルＤ_ωと近端話者信号スペクトルＳ_ωとの相関を踏まえた評価量εの最小解を新たに得ることで、エコー抑圧ゲイン推定の高精度化を図る考えである。

上記した式（１３）を変形するとゲイン係数Ｇ_ω ^ｐは次式で表せる。なお、Ｇ_ω ^ｐのｐは、上記した従来技術のゲイン係数Ｇ_ωと区別するための記号である。

ここで、γ_ωはエコーパワー補正係数であり、エコースペクトルＤ_ωと収音信号スペクトルＹ_ωの内積値を、エコースペクトルの二乗ノルム値‖Ｄ_ω‖^２で除した値を示す。式（１８）から分かるように未知の変数であるＤ_ωが存在するため、このままではゲインを求めることは困難である。

そこで、この発明では、エコースペクトルＤ_ωと収音信号スペクトルＹ_ωの内積値を、エコーパワー推定値の平方根｜Ｄ_ω＾｜と収音信号スペクトルの大きさ｜Ｙ_ω｜の内積値で代用すると共に、エコースペクトルの二乗ノルム値‖Ｄ_ω‖^２をエコーパワー推定値の平方根の二乗ノルム値‖Ｄ_ω＾‖^２で代用する。つまり、エコーパワー補正係数推定値γ_ω＾を式（１９）で計算する。

ここで内積＜｜Ｄ_ω＾｜，｜Ｙ_ω｜＞を例えば式（２０）で、エコーパワー推定値の平方根の二乗ノルム値‖Ｄ_ω＾‖^２を例えば式（２１）で、それぞれ計算する。

なお、式（２０）、式（２１）に置き換えて、周波数軸方向にも着目した式（２２）、式（２３）を用いてもよい。

そして、ゲイン係数を次式に示すように与える。

式（１９）は一般化すると次式で表せる。

エコーパワー補正係数推定値γ_ω＾は、エコーパワー推定値の平方根｜Ｄ_ω＾｜と収音信号スペクトルの大きさ｜Ｙ_ω｜との内積値を２分のα乗した値を、エコーパワー推定値の平方根のα乗ノルム値で除した値である。

つまり式（１９）は、α＝２とした場合であり、α＝２である式（２４）は、ウィナーフィルタに近いフィルタであることを表している。αは強調係数、その値は０より大きな整数である。なお、式（２４）を一般化すると次式で表せる。

式（２６）のエコーパワー補正係数推定値γ_ω＾は、上記式（２５）に示した一般化したγ_ω＾である。このエコーパワー補正係数推定値γ_ω＾は０〜無限大の範囲で値をとる。よって、エコーパワー補正係数推定値γ_ω＾が大きければゲイン係数Ｇ_ω ^ｐが小さくなりエコーが抑圧される。また、エコーパワー補正係数推定値γ_ω＾が小さければゲイン係数Ｇ_ω ^ｐが大きくなり近端話者信号の欠損を軽減するように動作する。その結果、エコーを十分に抑圧できるように動作し、音声歪みに対してより頑健な音質劣化の少ないエコー抑圧を行うことが可能になる。

図１に、この発明のエコー消去装置１００の機能構成例を示す。その動作フローを図２に示す。エコー消去装置１００は、再生信号周波数分析部８１と、収音信号周波数分析部８２と、エコーパワー推定部８３と、エコーパワー補正係数計算部１１０と、エコー抑圧ゲイン計算部１２０と、乗算部８５と、周波数合成部８６と、を具備する。エコー消去装置１００は、例えばＲＯＭ、ＲＡＭ、ＣＰＵ等で構成されるコンピュータに所定のプログラムが読み込まれて、ＣＰＵがそのプログラムを実行することで実現されるものである。

再生信号周波数分析部８１は、再生信号ｘ（ｋ）を周波数領域の再生信号スペクトルＸ_ωに変換する（ステップＳ８１）。収音信号周波数分析部８２は、収音信号ｙ（ｋ）を周波数領域の収音信号スペクトルＹ_ωに変換する（ステップＳ８２）。

エコーパワー推定部８３は、再生信号スペクトルＸ_ωと収音信号スペクトルＹ_ωとを入力としてエコーパワー推定値｜Ｄ_ω＾｜^２を計算する（ステップＳ８３）。エコーパワー補正係数計算部１１０は、エコーパワー推定値｜Ｄ_ω＾｜^２と収音信号スペクトルＹ_ωとを入力として、エコーパワー推定値の平方根｜Ｄ_ω＾｜と収音信号スペクトルの大きさ｜Ｙ_ω｜との内積値を２分のα乗した値を、エコーパワー推定値の平方根のα乗ノルム値‖Ｄ_ω＾‖^αで除した値をエコーパワー補正係数推定値γ_ω＾として計算する（ステップＳ１１０）。

エコー抑圧ゲイン計算部１２０は、収音信号スペクトルＹ_ωと、エコーパワー補正係数推定値γ_ω＾と、エコーパワー推定値｜Ｄ_ω＾｜^２とを入力として、エコーを抑圧するゲイン係数Ｇ_ω ^ｐをエコーパワー補正係数推定値γ_ω＾が大きな値をとるほど小さな値になるように計算する（ステップＳ１２０）。なお、乗算部８５と周波数合成部８６の処理ステップは省略している。

エコー消去装置１００は、従来のエコー消去装置１０（図７）の類似度係数計算部１２がエコーパワー補正係数計算部１１０に、エコー抑圧ゲイン計算部１４がエコー抑圧ゲイン計算部１２０に置き換わった以外は、従来技術と同じ構成である。この新しい構成のエコーパワー補正係数計算部１１０とエコー抑圧ゲイン計算部１２０について、更に詳しく説明する。

〔エコーパワー補正係数計算部〕
図３に、エコーパワー補正係数計算部１１０のより具体的な機能構成例を示す。エコーパワー補正係数計算部１１０は、内積演算手段１１０ａと、ノルム演算手段１１０ｂと、除算手段１１０ｃと、レジスタ１１０ｄと、レジスタ１１０ｅと、を備える。

内積演算手段１１０ａは、エコーパワー推定部８３が出力するエコーパワー推定値｜Ｄ_ω＾｜^２と収音信号周波数分析部８２が出力する収音信号スペクトルＹ_ωとを入力として、式（２７）でその内積を演算する。

ここでαは強調係数であり、その値はレジスタ１１０ｅに記憶されている。また、εは上記した忘却係数であり、例えばε＝０．０１６といった値が、予めレジスタ１１０ｄに記憶されている。εが１に近づくほど現在のエコーパワー推定値の平方根｜Ｄ_ω＾｜と収音信号スペクトルの大きさ｜Ｙ_ω｜に依存した（重み付けされた）それぞれの値になる。

ノルム演算手段１１０ｂは、式（２８）でエコーパワー推定値の平方根のα乗ノルム値‖Ｄ_ω＾‖^αを演算する。

除算手段１１０ｃは、その内積値を、エコーパワー推定値の平方根のα乗ノルム値‖Ｄ_ω＾‖^αで除してエコーパワー補正係数推定値γ_ω＾を演算する。エコーパワー補正係数推定値γ_ω＾は、式（２５）で表せる。

〔エコー抑圧ゲイン計算部〕
図４に、エコー抑圧ゲイン計算部１２０のより具体的な機能構成例を示す。エコー抑圧ゲイン計算部１２０は、エコーパワー調整手段１２０ａと、減算手段１２０ｂと、除算手段１２０ｃと、エコーの消し残りを無くす定数Ｃを記録したレジスタ１２０ｄと、強調係数αを記録したレジスタ１２０ｅとを備える。

エコーパワー調整手段１２０ａは、エコーパワー推定値｜Ｄ_ω＾｜^２とエコーパワー補正係数推定値γ_ω＾を入力として、エコーパワー推定値の平方根｜Ｄ_ω＾｜を強調係数αでべき乗した後にお互いを乗算し、さらにその値に定数Ｃを乗算してエコーパワー調整値を計算する（図５のステップＳ１２０ａ）。このエコーパワー調整過程を含むことで、エコー成分を十分に抑圧することが可能になる。

減算手段１２０ｂは、エコーパワー調整値と収音信号スペクトルＹ_ωとを入力とし、収音信号スペクトルの大きさ｜Ｙ_ω｜を強調係数αでべき乗した値からエコーパワー調整値を減算する（ステップＳ１２０ｂ）。

除算手段１２０ｃは、減算手段１２０ｂの出力信号を、収音信号スペクトルの大きさ｜Ｙ_ω｜を強調係数αでべき乗した値で除算してゲイン係数Ｇ_ω ^ｐを出力する（ステップＳ１２０ｃ）。つまり、エコー抑圧ゲイン計算部１２０は、式（２９）を計算してゲイン係数Ｇ_ω ^ｐを出力する。

式（２９）のエコーパワー補正係数推定値γ_ω＾は、上記式（２５）に示した一般化したγ_ω＾である。ここで強調係数αは正の数であるとして例えば２とする。定数Ｃは正の数であるとして例えば１.０とする。強調係数αと、定数Ｃは、ゲイン係数Ｇ_ω ^ｐを適当な値に調整するものである。強調係数α＝２とするとウィナーフィルタに近い式になる。

以上のようにして求めたゲイン係数Ｇ_ω ^ｐは、近端話者信号ｓ（ｋ）の欠損を軽減するようにエコー消去装置１００を動作させる。したがって、従来法に比べて音声歪みに対して頑健なミュージカルノイズ発生の少ないエコー抑圧を行うことが可能になる。
〔シミュレーション結果〕
この発明のエコー抑圧ゲイン推定方法を、短時間スペクトル振幅エコー抑圧処理に適用して従来法と性能比較を行った。従来法は、非特許文献１に記載されたＷＦ法とした。性能の比較は、音声特徴量の差（ベクトル距離）であるケプストラム距離を用いて行った。シミュレーション条件は、サンプリング周波数１６ｋＨｚ、周波数分析長２５６点、１/２オーバーラップ加算による周波数分析合成、エコー経路ｄ（ｋ）の残響時間は３００ｍｓの部屋とした。

シミュレーション結果を図６に示す。図６の横軸は時間[秒]、縦軸はケプストラム距離である。試験で用いた収音信号のうち、０秒から４秒までの区間はエコーのみの（近端話者信号が存在しない）区間、４秒から８秒までの区間は近端話者信号のみの（エコーが存在しない）区間、８秒から１２秒までの区間はエコーと近端話者信号が混在する区間である。このうち、図６には、８秒から１２秒までのエコーと近端話者信号が混在する区間を示した。破線がＷＦ法、実線が本発明である。図６より、この発明のエコー抑圧ゲイン推定方法の方が、ケプストラム距離の値が小さいことが分かる。つまり、全体を通して近端話者信号を再現できていることが分かる。

以上のようにこの発明のエコー抑圧ゲイン推定方法を用いることでエコー成分が残ることに起因する音質劣化を防止すると共に、高いエコー抑圧性能を達成しながらミュージカルノイズ発生も軽減することができる。

なお、この発明の方法及び装置は上述の実施形態に限定されるものではなく、この発明の趣旨を逸脱しない範囲で適宜変更が可能である。また、上記方法及び装置において説明した処理は、記載の順に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されるとしてもよい。

また、上記装置における処理手段をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、各装置における処理手段がコンピュータ上で実現される。

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。具体的には、例えば、磁気記録装置として、ハードディスク装置、フレキシブルディスク、磁気テープ等を、光ディスクとして、ＤＶＤ（Digital Versatile Disc）、ＤＶＤ−ＲＡＭ（Random Access Memory）、ＣＤ−ＲＯＭ（Compact Disc Read Only Memory）、ＣＤ−
Ｒ（Recordable）/ＲＷ（ReWritable）等を、光磁気記録媒体として、ＭＯ（Magneto Optical disc）等を、半導体メモリとしてＥＥＰ−ＲＯＭ（Electronically Erasable and Programmable-Read Only Memory）等を用いることができる。

また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記録装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。

また、各手段は、コンピュータ上で所定のプログラムを実行させることにより構成することにしてもよいし、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。

Claims (5)

1. 再生信号を周波数領域の再生信号スペクトルＸ_ωに変換する再生信号周波数分析段階と、
   収音信号を周波数領域の収音信号スペクトルＹ_ωに変換する収音信号周波数分析段階と、
   上記再生信号スペクトルＸ_ωと上記収音信号スペクトルＹ_ωとを入力としてエコーパワー推定値｜Ｄ_ω＾｜^２を計算するエコーパワー推定段階と、
   上記収音信号スペクトルＹ_ωと上記エコーパワー推定値｜Ｄ_ω＾｜^２とを入力として、上記エコーパワー推定値の平方根｜Ｄ_ω＾｜と上記収音信号スペクトルの大きさ｜Ｙ_ω｜との内積値を２分のα乗した値を、上記エコーパワー推定値の平方根のα乗ノルム値‖Ｄ_ω＾‖^αで除した値をエコーパワー補正係数推定値γ_ω＾として計算するエコーパワー補正係数計算段階と、
   上記収音信号スペクトルＹ_ωと上記エコーパワー補正係数推定値γ_ω＾と上記エコーパワー推定値｜Ｄ_ω＾｜^２とを入力として、エコーを抑圧するゲイン係数を上記エコーパワー補正係数推定値γ_ω＾が大きな値をとるほど小さな値になるように計算するエコー抑圧ゲイン計算段階と、
   を有するエコー抑圧ゲイン推定方法。
2. 請求項１に記載のエコー抑圧ゲイン推定方法において、
   上記エコー抑圧ゲイン計算段階は、上記エコーパワー推定値の平方根｜Ｄ_ω＾｜を上記αでべき乗した値に、上記エコーパワー補正係数推定値γ_ω＾を乗算したエコーパワー調整値を生成するエコーパワー調整過程を、含むことを特徴とするエコー抑圧ゲイン推定方法。
3. 請求項２に記載のエコー抑圧ゲイン推定方法において、
   上記エコー抑圧ゲイン計算段階は、
   上記エコーパワー調整過程と、
   上記収音信号スペクトルの大きさ｜Ｙ_ω｜を上記αでべき乗した値から上記エコーパワー調整値を減算する減算過程と、
   上記減算過程の出力信号を上記収音信号スペクトルの大きさ｜Ｙ_ω｜を上記αでべき乗した値で除算する除算過程と、
   から成ることを特徴とするエコー抑圧ゲイン推定方法。
4. 再生信号を周波数領域の再生信号スペクトルＸ_ωに変換する再生信号周波数分析部と、
   収音信号を周波数領域の収音信号スペクトルＹ_ωに変換する収音信号周波数分析部と、
   上記再生信号スペクトルＸ_ωと上記収音信号スペクトルＹ_ωとを入力としてエコーパワー推定値｜Ｄ_ω＾｜^２を計算するエコーパワー推定部と、
   上記収音信号スペクトルＹ_ωと上記エコーパワー推定値｜Ｄ_ω＾｜^２とを入力として、上記エコーパワー推定値の平方根｜Ｄ_ω＾｜と上記収音信号スペクトルの大きさ｜Ｙ_ω｜との内積値を２分のα乗した値を、上記エコーパワー推定値の平方根のα乗ノルム値‖Ｄ_ω＾‖^αで除した値をエコーパワー補正係数推定値γ_ω＾として計算するエコーパワー補正係数計算部と、
   上記収音信号スペクトルＹ_ωと上記エコーパワー補正係数推定値γ_ω＾と上記エコーパワー推定値｜Ｄ_ω＾｜^２とを入力として、エコーを抑圧するゲイン係数を上記エコーパワー補正係数推定値γ_ω＾が大きな値をとるほど小さな値になるように計算するエコー抑圧ゲイン計算部と、
   を具備するエコー抑圧ゲイン推定装置。
5. 請求項１乃至３の何れかに記載したエコー抑圧ゲイン推定方法を、コンピュータで処理するためのプログラム。

Images