(57)【要約】【課題】エコー成分が残ることに起因する音質劣化を防止すると共に、ミュージカルノイズ発生も抑圧する。【解決手段】この発明によるエコー抑圧ゲイン推定方法は、再生信号スペクトルの大きさ｜Ｘω｜と収音信号スペクトルの大きさ｜Ｙω｜との内積値＜｜Ｘω｜，｜Ｙω｜＞を、再生信号スペクトルの二乗ノルム値‖Ｘω‖２に音響結合量の推定値の平方根｜Ｈω＾｜を乗じた値で除した値をエコーパワー補正係数推定値として計算するエコーパワー補正係数計算段階と、収音信号スペクトルとエコーパワー補正係数推定値とエコーパワー推定値とを入力として、エコーを抑圧するゲイン係数を上記エコーパワー補正係数推定値が大きな値をとるほど小さな値になるように計算するエコー抑圧ゲイン計算段階と、を有し、０〜無限大の範囲の値をとるエコーパワー補正係数推定値が、近端話者信号の欠損を軽減するように動作する。【選択図】図１
2014017697.tif 000002
【特許請求の範囲】
【請求項１】
再生信号を周波数領域の再生信号スペクトルＸωに変換する再生信号周波数分析段階と、
収音信号を周波数領域の収音信号スペクトルＹωに変換する収音信号周波数分析段階と、
上記再生信号スペクトルＸωと上記収音信号スペクトルＹωとを入力として、エコーパワー推定値｜Ｄω＾｜２と音響結合量の推定値｜Ｈω＾｜２を計算するエコーパワー推定段階と、
上記再生信号スペクトルＸωと上記収音信号スペクトルＹωと上記音響結合量の推定値｜Ｈω＾｜２とを入力として、強調係数αとしてエコーパワー補正係数推定値γω＾を次式で計算する
【数２５】
2014017697.tif 000027
エコーパワー補正係数計算段階と、
上記収音信号スペクトルＹωと上記エコーパワー補正係数推定値γω＾と上記エコーパワー推定値｜Ｄω＾｜２とを入力として、エコーを抑圧するゲイン係数を上記エコーパワー補正係数推定値γω＾が大きな値をとるほど小さな値になるように計算するエコー抑圧ゲイン計算段階と、
を有するエコー抑圧ゲイン推定方法。
【請求項２】
請求項１に記載のエコー抑圧ゲイン推定方法において、
上記エコー抑圧ゲイン計算段階は、上記エコーパワー推定値の平方根｜Ｄω＾｜を上記αでべき乗した値に、上記エコーパワー補正係数推定値γω＾を乗算したエコーパワー調整値を生成するエコーパワー調整過程を、含むことを特徴とするエコー抑圧ゲイン推定方法。
【請求項３】
請求項２に記載のエコー抑圧ゲイン推定方法において、
上記エコー抑圧ゲイン計算段階は、
上記エコーパワー調整過程と、
上記収音信号スペクトルの大きさ｜Ｙω｜を上記αでべき乗した値から上記エコーパワー調整値を減算する減算過程と、
上記減算過程の出力信号を上記収音信号スペクトルの大きさ｜Ｙω｜を上記αでべき乗した値で除算する除算過程と、
から成ることを特徴とするエコー抑圧ゲイン推定方法。
【請求項４】
再生信号を周波数領域の再生信号スペクトルＸωに変換する再生信号周波数分析部と、
収音信号を周波数領域の収音信号スペクトルＹωに変換する収音信号周波数分析部と、
上記再生信号スペクトルＸωと上記収音信号スペクトルＹωとを入力として、エコーパワー推定値｜Ｄω＾｜２と音響結合量の推定値｜Ｈω＾｜２を計算するエコーパワー推定部と、
上記再生信号スペクトルＸωと上記収音信号スペクトルＹωと上記音響結合量の推定値｜Ｈω＾｜２とを入力として、強調係数αとしてエコーパワー補正係数推定値γω＾を次式で計算する
【数２６】
2014017697.tif 000028
エコーパワー補正係数計算部と、
上記収音信号スペクトルＹωと上記エコーパワー補正係数推定値γω＾と上記エコーパワー推定値｜Ｄω＾｜２とを入力として、エコーを抑圧するゲイン係数を上記エコーパワー補正係数推定値γω＾が大きな値をとるほど小さな値になるように計算するエコー抑圧ゲイン計算部と、
を具備するエコー抑圧ゲイン推定装置。
【請求項５】
請求項１乃至３の何れかに記載したエコー抑圧ゲイン推定方法を、コンピュータで処理するためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、音響再生系を有する通信会議システム等に用いられるエコー消去装置及び、それに適用されるエコー抑圧ゲイン推定方法と、そのプログラムに関する。
【背景技術】
【０００２】
短時間スペクトル振幅（STSA：Sort-Time Spectral Amplitude）推定に基づくエコー抑圧処理は、エコーと近端話者音声の無相関を仮定してエコーを抑圧するゲイン係数を推定し、振幅周波数領域でエコーを抑圧する方法であり、例えば非特許文献１に記載されている。これに対し、エコー抑圧後の音声品質向上のため、非特許文献１に記載される上記ゲイン係数の推定方法（以下、「エコー抑圧ゲイン推定方法」と称す）を改善した手法として、特許文献１が提案されている。ここでは、特許文献１に開示されたエコー抑圧ゲイン推定方法を用いたエコー消去装置１０の機能構成例を、図６に示してその動作を簡単に説明する。
【０００３】
エコー消去装置１０は、再生信号周波数分析部８１と、収音信号周波数分析部８２と、エコーパワー推定部８３と、類似度係数計算部１２と、エコー抑圧ゲイン計算部１４と、乗算部８５と、周波数合成部８６とを備える。再生信号ｘ（ｋ）は、例えばサンプリング周波数１６ｋＨｚで離散値とされた信号であり、スピーカ１によって音響信号に変換される。なお、再生信号を離散値化するＡＤ変換器と、その離散値を連続値に変換するＤＡ変換器については省略している。
【０００４】
再生信号周波数分析部８１は、再生信号の離散値を２５６点集めて１フレームとし、１/２オーバーラップ加算による周波数分析を行い、フレーム単位で８ｋＨｚまでの周波数範囲を１２８個の再生信号スペクトルＸωに変換する。ｋはフレーム番号であり、ωはこの例の場合６４Ｈｚの間隔で得られる周波数スペクトルの番号（１〜１２８）である。
【０００５】
マイクロホン２で収音される収音信号ｙ（ｋ）は、再生信号ｘ（ｋ）がスピーカ１から再生されたことによるエコーｄ（ｋ）が、近端話者信号ｓ（ｋ）に重畳した信号である。その収音信号ｙ（ｋ）は、再生信号ｘ（ｋ）と同じように収音信号周波数分析部８２によって収音信号スペクトルＹωに変換される。ＹωはＹω＝Ｄω＋Ｓω、Ｄωはエコースペクトル、Ｓωは近端話者信号スペクトル、である。
【０００６】
エコーパワー推定部８３は、再生信号スペクトルＸωと収音信号スペクトルＹωとを入力として式（１）に示すエコーパワー推定値｜Ｄω＾｜２を計算する。＾は推定値であることを表わすが、その表記は式及び図中に示すものが正しい。
【数１】
2014017697.tif 000003
【０００７】
ここで、Ｈωはスピーカ１からマイクロホン２に回り込むエコー経路のスペクトル、｜Ｈω＾｜２は音響結合量の推定値であり、｜Ｈω＾｜２￣は１フレーム過去の音響結合量の推定値である。ｍｉｎ｛  ，  ｝は最小値を選択する関数である。したがって、｜Ｈω＾｜２は隣接するフレーム間で小さい方の音響結合量がその推定値になる。
【０００８】
エコー抑圧ゲイン計算部１４は、収音信号スペクトルＹωと、エコーパワー推定値｜Ｄω＾｜２と、類似度係数｜ｒω｜を入力として式（２）に示すゲイン係数Ｇωを出力する。ゲイン係数Ｇωは０〜１の実数値をとり、収音信号スペクトルＹω中にエコー成分が多い場合には小さな値、エコー成分以外の成分が多い場合には大きな値をとる。
【数２】
2014017697.tif 000004
【０００９】
乗算部８５は、収音信号スペクトルＹωにゲイン係数Ｇωを乗算する。エコー成分が多い場合のゲイン係数Ｇωは小さな値となるので、乗算部８５の出力信号はエコー成分が抑圧された近端話者信号スペクトル推定値Ｓω＾となる。この各周波数成分ωに対応する近端話者信号スペクトル推定値Ｓω＾は、周波数合成部８６において時間領域の出力信号ｅ（ｋ）に再合成される。
【００１０】
類似度係数計算部１２は、再生信号スペクトルＸωと収音信号スペクトルＹωとを入力として、内積＜Ｘω，Ｙω＞を例えば式（３）で、再生信号スペクトルのノルム値‖Ｘω‖と、収音信号スペクトルのノルム値‖Ｙω‖を、例えば式（４）と式（５）でそれぞれ計算する。
【数３】
2014017697.tif 000005
【００１１】
ここで、＊は複素共役、￣は１フレーム過去であることを表わす。εは０＜ε≦１を満たす忘却係数であり、指数関数的な減衰の時定数を決定する。例えばε＝０.０１６とする。εが１に近づくほど現在の再生信号スペクトルＸωと収音信号スペクトルＹωに依存した（重み付けされた）それぞれの値になる。なお、周波数軸方向にも着目した式（６）〜（８）を用いてもよい。
【数４】
2014017697.tif 000006
【００１２】
ここでＭ1〜Ｍ2は所定の周波数範囲を表わす。式（６）〜（８）は周波数軸方向に相関を取った後に時間軸方向の相関を取る形である。相関を取る順番は逆でもよい。また、式（３）に替えて位相成分を考慮しない絶対値から求める式（９）を用いてもよい。
【数５】
2014017697.tif 000007
【００１３】
類似度係数計算部１２は、式（３）〜（５）でそれぞれの値を計算した後に、式（１０）で類似度係数｜ｒω｜を計算する。
【数６】
2014017697.tif 000008
【００１４】
ここで類似度係数｜ｒω｜を算出する式（１０）の導出根拠について説明する。
例えばウィナーフィルタ法（Wiener Filter以下「ＷＦ法」と称す）では、式（１１）の評価量εを最小とするゲイン係数Ｇωを推定してエコーを抑圧する。
【数７】
2014017697.tif 000009
【００１５】
式（１１）をＧωで微分すると次式で表せる。
【数８】
2014017697.tif 000010
【００１６】
式（１２）を０と置いた式からゲイン係数Ｇωを求めると式（１３）となる。
【数９】
2014017697.tif 000011
【００１７】
式（１３）を変形するとゲイン係数Ｇωは式（１４）で表せる。
【数１０】
2014017697.tif 000012
【００１８】
ここでｒωはエコースペクトルＤωと近端話者信号スペクトルＳωの複素コヒーレンスを示す。Ｄω＝Ｈω＊Ｘωを用いてｒωから未知のベクトルＤωを消去すると式（１５）で表せる。
【数１１】
2014017697.tif 000013
【００１９】
ここで＊は複素共役を表す。式（１５）から分かるようにＤωを消去しても未知の変数であるＨωが存在するため、ｒωを求めることは困難である。しかし、その絶対値ならば式（１６）に示すようにＨωが消去可能なことに着目する。
【数１２】
2014017697.tif 000014
【００２０】
そしてゲイン係数を次式に示すように与える。
【数１３】
2014017697.tif 000015
【００２１】
この式（１７）から明らかなようにゲイン係数Ｇωは、類似度係数｜ｒω｜が１に近づくほど小さな値になりエコーが抑圧される。と同時にこのようにして求めたゲイン係数は、近端話者信号ｓ（ｋ）の欠損を軽減するように動作する。つまり、＜Ｄω，Ｓω＞＝０となるほどの長時間を要することなく、短時間のスペクトル振幅を用いてエコーを抑圧するので、内積値が０でない分が誤差となって発生するミュージカルノイズの発生を抑圧することができる。
【先行技術文献】
【特許文献】
【００２２】
【特許文献１】特許第４７８７８５１号
【非特許文献】
【００２３】
【非特許文献１】阪内澄宇、羽田陽一、片岡章俊著「ＳＴＳＡ推定に基づくエコー抑圧処理のゲイン強調化方式」信学論（Ａ）,１vol.J88-A,no.6,Jun.2005,p695-703
【発明の概要】
【発明が解決しようとする課題】
【００２４】
類似度係数｜ｒω｜は、収音信号中のエコー成分の割合を表す値で０〜１の間の値をとる係数である。したがって、近端話者信号ｓ（ｋ）の欠損を軽減するように動作するが、類似度係数｜ｒω｜は１より大きな値をとることができないため、エコーをより大きく抑圧するようには動作しない。その結果、エコーが十分に抑圧できない場合があり、抑圧し切れなかったエコー成分が音質劣化として聴感上知覚される場合がある。エコーが十分に抑圧できない場合とは、近端話者信号ｓ（ｋ）の相関が残ってエコーが真値より小さく見積もられてしまう場合である。また、近端話者信号ｓ（ｋ）の相関が残る現象は、近端話者信号ｓ（ｋ）が大きい場合に顕著になる。
【００２５】
この発明は、このような課題に鑑みてなされたものであり、エコー成分が残ることに起因する音質劣化を防止するエコー抑圧ゲイン推定方法と、それを用いたエコー消去装置とそのプログラムを提供することを目的とする。
【課題を解決するための手段】
【００２６】
この発明のエコー抑圧ゲイン推定方法は、再生信号周波数分析段階と、収音信号周波数分析段階と、エコーパワー推定段階と、エコーパワー補正係数計算段階と、エコー抑圧ゲイン計算段階と、を有する。再生信号周波数分析段階は、再生信号を周波数領域の再生信号スペクトルＸωに変換する。収音信号周波数分析段階は、収音信号を周波数領域の収音信号スペクトルＹωに変換する。エコーパワー推定段階は、再生信号スペクトルＸωと収音信号スペクトルＹωとを入力としてエコーパワー推定値｜Ｄω＾｜２と音響結合量の推定値｜Ｈω＾｜２を計算する。エコーパワー補正係数計算段階は、再生信号スペクトルＸωと収音信号スペクトルＹωと音響結合量の推定値｜Ｈω＾｜２とを入力として、強調係数αとしてエコーパワー補正係数推定値γω＾を次式で計算する。
【００２７】
【数１４】
2014017697.tif 000016
【００２８】
エコー抑圧ゲイン計算段階は、収音信号スペクトルＹωと、エコーパワー補正係数推定値γω＾と、エコーパワー推定値｜Ｄω＾｜２とを入力として、エコーを抑圧するゲイン係数をエコーパワー補正係数推定値γω＾が大きな値をとるほど小さな値になるように計算する。
【発明の効果】
【００２９】
この発明のエコー抑圧ゲイン推定方法によれば、エコーパワー補正係数推定値が０〜無限大の範囲の値をとり、エコーパワー補正係数推定値が大きな値をとるほどゲイン係数が小さくなりエコーが抑圧され、逆にエコーパワー補正係数推定値が小さな値をとるほどゲイン係数が大きくなり近端話者信号の欠損を軽減するように動作する。したがって、従来法に比べてエコーを十分に抑圧できるように動作し、音声歪みに対してより頑健な音質劣化の少ないエコー抑圧を行うことが可能になる。
【図面の簡単な説明】
【００３０】
【図１】この発明のエコー抑圧ゲイン推定方法を用いたエコー消去装置１００の機能構成例を示す図。
【図２】この発明のエコー抑圧ゲイン推定方法の動作フローを示す図。
【図３】エコーパワー補正係数計算部１１０の機能構成例を示す図。
【図４】エコー抑圧ゲイン計算部１２０の機能構成例を示す図。
【図５】エコー抑圧ゲイン計算部１２０の動作フローを含めたエコー抑圧ゲイン推定方法の動作フローを示す図。
【図６】従来のエコー消去装置１０の機能構成を示す図。
【発明を実施するための形態】
【００３１】
以下、この発明の実施の形態を図面を参照して説明する。複数の図面中同一のものには同じ参照符号を付し、説明は繰り返さない。
【００３２】
〔発明の基本的な考え〕
実施例の説明の前に、この発明のエコー抑圧方法の基本的な考えを説明する。この発明は、エコースペクトルＤωと近端話者信号スペクトルＳωとの相関を踏まえた評価量εの最小解を新たに得ることで、エコー抑圧ゲイン推定の高精度化を図る考えである。
【００３３】
上記した式（１３）を変形するとゲイン係数Ｇωｐは次式で表せる。なお、Ｇωｐのｐは、上記した従来技術のゲイン係数Ｇωと区別するための記号である。
【数１５】
2014017697.tif 000017
【００３４】
ここで、γωはエコーパワー補正係数であり、エコースペクトルＤωと収音信号スペクトルＹωの内積値を、エコースペクトルの二乗ノルム値‖Ｄω‖２で除した値を示す。式（１８）から分かるように未知の変数であるＤωが存在するため、このままではゲインを求めることは困難である。
【００３５】
そこで、この発明では、上記した式（１）の｜Ｄω＾｜２＝｜Ｈω＾｜２｜Ｘω｜２の関係を用いてエコーパワー補正係数推定値γω＾を推定する。エコーパワー補正係数推定値γω＾は式（１９）に示すように与える。
【数１６】
2014017697.tif 000018
【００３６】
エコーパワー補正係数推定値γω＾は、再生信号スペクトルの大きさ｜Ｘω｜と収音信号スペクトルの大きさ｜Ｙω｜との内積値＜｜Ｘω｜，｜Ｙω｜＞を、再生信号スペクトルの二乗ノルム値‖Ｘω‖２に音響結合量の推定値の平方根｜Ｈω＾｜を乗じた値で除した値である。
【００３７】
ここで内積＜｜Ｘω｜，｜Ｙω｜＞を例えば式（２０）で、再生信号スペクトルの二乗ノルム値‖Ｘω‖２を例えば式（２１）で、それぞれ計算する。
【数１７】
2014017697.tif 000019
【００３８】
なお、式（２０）、式（２１）に置き換えて、周波数軸方向にも着目した式（２２）、式（２２）を用いてもよい。
【数１８】
2014017697.tif 000020
【００３９】
そして、ゲイン係数を次式に示すように与える。
【数１９】
2014017697.tif 000021
【００４０】
式（１９）を一般化すると次式で表せる。
【数２０】
2014017697.tif 000022
【００４１】
つまり式（１９）は、α＝２とした場合であり、α＝２である式（２４）は、ウィナーフィルタに近いフィルタであることを表している。αは強調係数、その値は０より大きな整数である。なお、式（２４）を一般化すると次式で表せる。
【数２１】
2014017697.tif 000023
【００４２】
式（２６）のエコーパワー補正係数推定値γω＾は、上記式（２５）に示した一般化したγω＾である。エコーパワー補正係数推定値γω＾は０〜無限大の範囲で値をとる。よって、エコーパワー補正係数推定値γω＾が大きければゲイン係数Ｇωｐが小さくなりエコーが抑圧される。また、エコーパワー補正係数推定値γω＾が小さければゲイン係数Ｇωｐが大きくなり近端話者信号の欠損を軽減するように動作する。その結果、エコーを十分に抑圧できるように動作し、音声歪みに対してより頑健な音質劣化の少ないエコー抑圧を行うことが可能になる。
【実施例１】
【００４３】
図１に、この発明のエコー消去装置１００の機能構成例を示す。その動作フローを図２に示す。エコー消去装置１００は、再生信号周波数分析部８１と、収音信号周波数分析部８２と、エコーパワー推定部８３と、エコーパワー補正係数計算部１１０と、エコー抑圧ゲイン計算部１２０と、乗算部８５と、周波数合成部８６と、を具備する。エコー消去装置１００は、例えばＲＯＭ、ＲＡＭ、ＣＰＵ等で構成されるコンピュータに所定のプログラムが読み込まれて、ＣＰＵがそのプログラムを実行することで実現されるものである。
【００４４】
再生信号周波数分析部８１は、再生信号ｘ（ｋ）を周波数領域の再生信号スペクトルＸωに変換する（ステップＳ８１）。収音信号周波数分析部８２は、収音信号ｙ（ｋ）を周波数領域の収音信号スペクトルＹωに変換する（ステップＳ８２）。
【００４５】
エコーパワー推定部８３は、再生信号スペクトルＸωと収音信号スペクトルＹωとを入力としてエコーパワー推定値｜Ｄω＾｜２と音響結合量の推定値｜Ｈω＾｜２を、上記式（１）で計算する（ステップＳ８３）。エコーパワー補正係数計算部１１０は、再生信号スペクトルＸωと収音信号スペクトルＹωと音響結合量の推定値｜Ｈω＾｜２とを入力として、エコーパワー補正係数推定値γω＾を計算する（ステップＳ１１０）。
【００４６】
エコー抑圧ゲイン計算部１２０は、収音信号スペクトルＹωと、エコーパワー補正係数推定値γω＾と、エコーパワー推定値｜Ｄω＾｜２とを入力として、エコーを抑圧するゲイン係数Ｇωｐをエコーパワー補正係数推定値γω＾が大きな値をとるほど小さな値になるように計算する（ステップＳ１２０）。なお、乗算部８５と周波数合成部８６の処理ステップは省略している。
【００４７】
エコー消去装置１００は、従来のエコー消去装置１０（図６）の類似度係数計算部１２がエコーパワー補正係数計算部１１０に、エコー抑圧ゲイン計算部１４がエコー抑圧ゲイン計算部１２０に置き換わった以外は、従来技術と同じ構成である。この新しい構成のエコーパワー補正係数計算部１１０とエコー抑圧ゲイン計算部１２０について、更に詳しく説明する。
【００４８】
〔エコーパワー補正係数計算部〕
図３に、エコーパワー補正係数計算部１１０のより具体的な機能構成例を示す。エコーパワー補正係数計算部１１０は、内積演算手段１１０ａと、ノルム演算手段１１０ｂと、除算手段１１０ｃと、レジスタ１１０ｄと、レジスタ１１０ｅと、を備える。
【００４９】
内積演算手段１１０ａは、再生信号スペクトルＸωと収音信号スペクトルＹωとを入力として、式（２７）でその内積を演算する。
【数２２】
2014017697.tif 000024
【００５０】
ここでαは強調係数であり、その値はレジスタ１１０ｅに記憶されている。また、εは上記した忘却係数であり、例えばε＝０．０１６といった値が、予めレジスタ１１０ｄに記憶されている。εが１に近づくほど現在の再生信号スペクトルの大きさ｜Ｘω｜と収音信号スペクトルの大きさ｜Ｙω｜に依存した（重み付けされた）それぞれの値になる。
【００５１】
ノルム演算手段１１０ｂは、式（２８）で再生信号スペクトルのα乗ノルム値‖Ｘω‖αを演算する。
【数２３】
2014017697.tif 000025
【００５２】
除算手段１１０ｃは、内積値と再生信号スペクトルのα乗ノルム値‖Ｘω‖αと音響結合量の推定値｜Ｈω＾｜２とを入力として、上記した式（２５）の演算をしてエコーパワー補正係数推定値γω＾を出力する。
【００５３】
〔エコー抑圧ゲイン計算部〕
図４に、エコー抑圧ゲイン計算部１２０のより具体的な機能構成例を示す。エコー抑圧ゲイン計算部１２０は、エコーパワー調整手段１２０ａと、減算手段１２０ｂと、除算手段１２０ｃと、エコーの消し残りを無くす定数Ｃを記録したレジスタ１２０ｄと、強調係数αを記録したレジスタ１２０ｅとを備える。
【００５４】
エコーパワー調整手段１２０ａは、エコーパワー推定値｜Ｄω＾｜２とエコーパワー補正係数推定値γω＾を入力として、エコーパワー推定値の平方根｜Ｄω＾｜を強調係数αでべき乗した後にお互いを乗算し、さらにその値に定数Ｃを乗算してエコーパワー調整値を計算する（図５のステップＳ１２０ａ）。このエコーパワー調整過程を含むことで、エコー成分を十分に抑圧することが可能になる。
【００５５】
減算手段１２０ｂは、エコーパワー調整値と収音信号スペクトルＹωとを入力とし、収音信号スペクトルＹωの大きさを強調係数αでべき乗した値からエコーパワー調整値を減算する（ステップＳ１２０ｂ）。
【００５６】
除算手段１２０ｃは、減算手段１２０ｂの出力信号を、収音信号スペクトルＹωの大きさを強調係数αでべき乗した値で除算してゲイン係数Ｇωｐを出力する（ステップＳ１２０ｃ）。つまり、エコー抑圧ゲイン計算部１２０は、式（２５）を計算してゲイン係数Ｇωｐを出力する。
【数２４】
2014017697.tif 000026
【００５７】
式（２９）のエコーパワー補正係数推定値γω＾は、上記式（２５）に示した一般化したγω＾である。ここで強調係数αは正の数であるとして例えば２とする。定数Ｃは正の数であるとして例えば１.０とする。強調係数αと、定数Ｃは、ゲイン係数Ｇωｐを適当な値に調整するものである。強調係数α＝２とするとウィナーフィルタに近い式になる。
【００５８】
以上のようにして求めたゲイン係数Ｇωｐは、近端話者信号ｓ（ｋ）の欠損を軽減するようにエコー消去装置１００を動作させる。したがって、従来法に比べて音声歪みに対して頑健なミュージカルノイズ発生の少ないエコー抑圧を行うことが可能になる。
【００５９】
なお、この発明の方法及び装置は上述の実施形態に限定されるものではなく、この発明の趣旨を逸脱しない範囲で適宜変更が可能である。例えば、音響結合量の推定値｜Ｈω＾｜は、エコーパワー推定部８３で計算したものを、エコーパワー補正係数計算部１１０で用いる例で説明を行ったが、再生信号スペクトルＸωと収音信号スペクトルＹωとが入力されるエコーパワー補正係数計算部１１０で音響結合量の推定値｜Ｈω＾｜を求めるようにしても良い。
【００６０】
また、上記方法及び装置において説明した処理は、記載の順に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されるとしてもよい。
【００６１】
また、上記装置における処理手段をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、各装置における処理手段がコンピュータ上で実現される。
【００６２】
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。具体的には、例えば、磁気記録装置として、ハードディスク装置、フレキシブルディスク、磁気テープ等を、光ディスクとして、ＤＶＤ（Digital Versatile Disc）、ＤＶＤ−ＲＡＭ（Random Access Memory）、ＣＤ−ＲＯＭ（Compact Disc Read Only Memory）、ＣＤ−
Ｒ（Recordable）/ＲＷ（ReWritable）等を、光磁気記録媒体として、ＭＯ（Magneto Optical disc）等を、半導体メモリとしてＥＥＰ−ＲＯＭ（Electronically Erasable and Programmable-Read Only Memory）等を用いることができる。
【００６３】
また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記録装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。
【００６４】
また、各手段は、コンピュータ上で所定のプログラムを実行させることにより構成することにしてもよいし、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
【図１】
2014017697.tif 000029
【図２】
2014017697.tif 000030
【図３】
2014017697.tif 000031
【図４】
2014017697.tif 000032
【図５】
2014017697.tif 000033
【図６】
2014017697.tif 000034