JP6432384B2 - エコー抑圧装置、エコー抑圧プログラム及びエコー抑圧方法 - Google Patents

Description

本発明は、エコー抑圧装置、エコー抑圧プログラム及びエコー抑圧方法に関し、例えば、テレビ会議システムや電話会議システム等において用いられるエコー抑圧装置、エコー抑圧プログラム及びエコー抑圧方法に適用し得るものである。

例えば、テレビ会議システムや電話会議システム等の拡声通話システムでは、スピーカから放音された音（ここで、「音」は音響や音声等を含む。）がマイクに回り込んで送話側に戻る音響エコー信号が発生する。音響エコー信号は、通話の著しい妨げとなるため、音響エコー抑圧方法に関して、これまでも多くの研究、開発が行なわれている。

音響エコー信号を抑圧する１つの手法として、エコーサプレッサー（エコー抑圧装置）を使用する手法がある。エコーサプレッサーとは、遠端信号と近端入力信号とからエコーパス特性、推定エコー信号、エコーサプレスゲインを求めて、近端入力信号とエコーサプレスゲインを乗算することで音響エコー信号を抑圧する手法である。このエコーサプレッサーを使用して音響エコー信号を抑圧することができる。

エコーサプレッサーでエコーパス特性を算出し、更新するときに、重要になるのが遅延量推定である。例えば、遅延量が既知で一定の場合は、遠端信号を遅延量だけ遅延させることで、音響エコー信号を抑圧することができる。

しかし、例えば、ソフトフォン（Ｓｏｆｔｐｈｏｎｅ）やスマートフォン（Ｓｍａｒｔｐｈｏｎｅ）などの音入出力部に、処理待ちデータを一時保持し出力するバッファを備える装置では遅延量が大きくなる。また、ソフトフォンやスマートフォンはエコー抑圧処理以外にも様々な処理（例えば、アプリケーション、プログラム、制御信号処理等）が動作しており、負荷が大きくなるとバッファに音データを記録する処理やバッファから音データを出力する処理が停止するため、遅延量が大きくなり、かつ遅延量の変動が大きくなる。

このような、ソフトフォンやスマートフォンのエコーサプレッサーは遅延量を推定する必要があり、遅延量を推定する手法が特許文献１によって提案されている。

特許文献１に開示される遅延量推定手法は、スピーカからパルス音や報知音を出力し、近端側の空間を伝達してマイクで受音され、スピーカに出力したパルス音や報知音と受音した応答パルス音や報知音応答を解析することで遅延量を推定している。

特開２０１３−２２５７４７号公報

しかしながら、特許文献１の遅延量推定では、パルス音や報知音と似たような音が近端側から入力されると遅延量が正しく推定できない。さらに、遅延量を推定するのにパルス音や報知音を出力すると、非常に耳障りであり、通話の著しい妨げとなる。

パルス音や報知音を出力しない場合は、シングルトーク時に、遠端側から入力されスピーカから出力した過去の遠端信号と、マイクで受音される近端入力信号を例えば相関係数を計算して、相関の高いところがマイクに音響エコー信号として入力されているとし、相関が高い時刻と現時刻との差を遅延量することで遅延量を推定できる。しかし、比較する過去の遠端信号で、同じ音や似た音（例えば同じ母音など）が入力されていると、相関係数がともに高くなるため、遅延量を誤推定し、音響エコー信号を抑圧できない。

そのため、パルス音や報知音などの音を出力しないで、遠端側から同じ音や似た音が入力し、スピーカから出力され近端側の空間を伝達し、音響エコー信号として入力されても、音響エコー信号の遅延量を推定できるエコー抑圧装置、エコー抑圧プログラム及びエコー抑圧方法が望まれている。

本発明は、上記課題に鑑みてなされたものであり、遅延量が大きい機器や遅延量の変動が大きい機器でもパルス音や報知音などの音を出力せずに、遠端側から同じ音や似た音が入力されても、相関が高いフレームから遅延量の候補値を計算し、その遅延量の候補値から過去の信号をと過去の近端入力信号を比較することで、遅延量の候補値から過去の信号を比較して音響エコー信号の遅延量を決定する遅延量推定を提供しようとするものである。

第１の本発明のエコー抑圧装置は、近端入力信号と遠端信号に基づく音響エコー信号を抑圧するエコー抑圧装置において、（１）過去の複数フレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとに基づいてフレーム誤差を計算するフレーム誤差計算手段と、（２）フレーム誤差計算手段で計算したフレーム毎の誤差に基づいて、１又は複数のフレーム遅延量の候補値を計算するフレーム遅延量候補値計算手段と、（３）フレーム遅延量候補値計算手段で判定された各フレーム遅延量の候補値に基づいて、所定フレーム離れたフレームの推定エコー信号と所定フレーム離れたフレームの近端入力信号とを用いて、フレーム遅延量を推定して出力するフレーム遅延量推定手段とを有することを特徴とする。

第２の本発明のエコー抑圧プログラムは、近端入力信号から遠端出力信号に基づく音響エコー信号を抑圧するエコー抑圧プログラムにおいて、コンピュータを、（１）過去の複数フレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとに基づいてフレーム誤差を計算するフレーム誤差計算手段と、（２）フレーム誤差計算手段で計算したフレーム毎の誤差に基づいて、１又は複数のフレーム遅延量の候補値を計算するフレーム遅延量候補値計算手段と、（３）フレーム遅延量候補値計算手段で判定された各フレーム遅延量の候補値に基づいて、所定フレーム離れたフレームの推定エコー信号と所定フレーム離れたフレームの近端入力信号とを用いて、フレーム遅延量を推定して出力するフレーム遅延量推定手段として機能させることを特徴とする。

第３の本発明のエコー抑圧方法は、近端入力信号から遠端出力信号に基づく音響エコー信号を抑圧するエコー抑圧方法において、フレーム誤差計算手段、フレーム遅延量候補値計算手段、及びフレーム遅延量推定手段を有し、（１）フレーム誤差計算手段は、過去の複数フレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとに基づいてフレーム誤差を計算し、（２）フレーム遅延量候補値計算手段は、フレーム誤差計算手段で計算したフレーム毎の誤差に基づいて、１又は複数のフレーム遅延量の候補値を計算し、（３）フレーム遅延量推定手段は、フレーム遅延量候補値計算手段で判定された各フレーム遅延量の候補値に基づいて、所定フレーム離れたフレームの推定エコー信号と所定フレーム離れたフレームの近端入力信号とを用いて、フレーム遅延量を推定して出力することを特徴とする。

本発明によれば、遠端信号に同じ音や似た音が入力されても、その同じ音や似た音の過去の信号を比較し、遅延量を推定し、推定した遅延量を推定エコー信号、推定エコーパス特性、サプレスゲインを算出することで、音響エコー信号を抑圧することができる。

第１の実施形態のエコー抑圧装置の構成を示すブロック図である。 第１の実施形態のフレーム遅延量推定部の処理を示す処理フロー図である。 第２の実施形態のフレーム遅延量推定部の処理を示す処理フロー図である。 第１及び第２の実施形態のフレーム遅延量推定処理の説明図である。 第３の実施形態のフレーム遅延量推定処理の説明図である。 従来のエコー抑圧装置を説明するためのブロック図である。

（Ａ）第１の実施形態
以下では、本発明に係るエコー抑圧装置、エコー抑圧プログラム及びエコー抑圧方法の第１の実施形態を、図面を参照しながら詳細に説明する。

第１の実施形態は、例えば、テレビ会議システムや電話会議システム等の拡声通話システムの音声送受信装置のエコー抑圧装置、エコー抑圧プログラム及びエコー抑圧方法に本発明を適用する場合を例示する。

（Ａ−１）第１の実施形態の構成
図１は、第１の実施形態に係るエコー抑圧装置１００の構成を示すブロック図である。

第１の実施形態に係るエコー抑圧装置１００は、上述した本発明の遅延量推定装置をエコー抑圧装置に実装したものであり、遠端信号が音響エコー信号として近端側のマイクに回り込んで入力されたときに遅延量を推定し、近端入力信号に対してエコー抑圧処理を行い、音響エコー信号を適切に抑圧するものである。

第１の実施形態のエコー抑圧装置１００は、例えば専用ボードとして構築されるようにしても良いし、ＤＳＰ（デジタルシグナルプロセッサ）へのエコー抑圧プログラムの書き込みによって実現されたものであっても良く、ＣＰＵと、ＣＰＵが実行するソフトウェア（エコー抑圧プログラム）によって実現されたものであっても良いが、機能的には、図１で表すことができる。

図１において、第１の実施形態に係るエコー抑圧装置１００は、遠端信号入力端子１０１、出力バッファ１０２、ＤＡ変換器１０３、スピーカ１０４、マイク１０５、ＡＤ変換器１０６、入力バッファ１０７、遠端信号周波数領域変換部１０８、遠端信号振幅スペクトル計算部１０９、エコーパス特性更新部１１０、推定エコー信号計算部１１１、推定エコー信号保持部１１２、近端入力信号周波数領域変換部１１３、近端入力信号振幅スペクトル計算部１１４、フレーム遅延量推定部１１５、推定エコー信号遅延部１１６、エコーサプレスゲイン計算部１１７、エコーサプレス部１１８、近端出力信号時間領域変換部１１９、近端信号出力端子１２０、遠端信号保持部１２１、シングルトーク判定部１２３を有する。

遠端信号入力端子１０１は、例えば、インターネットプロトコル（ＩＰ）網等のネットワークや、携帯電話等の無線ネットワークの電波等に接続されており、接続されている回線を介して遠端側（相手側）の遠端信号が入力される。

遠端信号入力端子１０１に入力された遠端信号は、出力バッファ１０２において一時保持され、ＤＡ変換器１０３の処理準備ができ次第、ＤＡ変換器１０３に遠端信号を出力し、ＤＡ変換器１０３において、デジタル音信号からアナログ音信号に変換され、スピーカ１０４を通して近端側に出力される。

出力バッファ１０２は、エコー抑圧装置１００の処理待ち用のバッファである。

一方、近端側の話者が発した音声等の音信号や、環境音、音響エコー信号（例えば、スピーカ１０４から出力されたアナログ音信号が近端側の空間を伝達して回り込んだ信号）等が重畳したアナログ音信号は、マイク１０５において受音され、ＡＤ変換器１０５おいてデジタル音信号に変換され、入力バッファ１０７おいて一時保持され、エコー抑圧装置１００の準備ができ次第、デジタル音信号を近端入力信号としてエコー抑圧装置１００に入力される。

入力バッファ１０７は、エコー抑圧装置１００の処理待ち用のバッファである。

遠端信号周波数領域変換部１０８は、例えば、高速フーリエ変換（ＦＦＴ）等により、遠端信号を周波数領域の信号に変換し、遠端信号の周波数スペクトルを出力する。

遠端信号振幅スペクトル計算部１０９は、遠端信号の周波数スペクトルに基づいて、遠端信号の振幅スペクトルを算出し、算出した遠端信号の振幅スペクトルを出力する。

エコーパス特性更新部１１０は、前フレームに推定したエコーパス特性を保持しており、保持している推定エコーパス特性に出力する。

さらに、エコーパス特性更新部１１０は、シングルトーク判定部１２３でシングルトークと判定されると、推定したフレーム遅延量だけ遅延させた遠端信号の振幅スペクトルと近端入力信号の振幅スペクトルに基づいて、現フレームのエコーパス特性を算出し、算出した現フレームのエコーパス特性とエコーパス特性更新部１１０に保持している推定エコーパス特性に基づき、推定エコーパス特性を更新し、更新したエコーパス特性をエコーパス特性更新部１１０に保持する。

推定エコー信号計算部１１１は、遠端信号の振幅スペクトルと推定エコーパス特性とを乗じて推定エコー信号の振幅スペクトルを算出し、算出した推定エコー信号の振幅スペクトルを出力する。

推定エコー信号保持部１１２は、推定エコー信号の振幅スペクトルを所定時間保持し、フレーム遅延量推定部１１５でフレーム遅延量を推定するために過去の推定エコー信号を複数フレーム出力する。さらに、推定エコー信号保持部１１２は、後述するフレーム遅延量推定部１１５で求まったフレーム遅延量だけ遅延させた推定エコー信号を出力する。

近端信号周波数領域変換部１１３は、例えば、高速フーリエ変換等により、近端入力信号を周波数領域の信号に変換し、近端入力信号の周波数スペクトルを出力する。

近端入力信号振幅スペクトル計算部１１４は、近端入力信号の周波数スペクトルに基づいて、近端入力信号の振幅スペクトルを算出し、算出した近端入力信号の振幅スペクトルを出力する。

フレーム遅延量推定部１１５は、過去の複数フレームの推定エコー信号の振幅スペクトルと、過去の複数フレームの近端入力信号の振幅スペクトルでフレーム遅延量を推定し、推定した遅延量を出力する。

推定エコー信号遅延部１１６は、フレーム遅延量を用いて推定エコー信号保持部１１２から当該フレーム遅延量だけ遅延させたフレームの推定エコー信号の振幅スペクトルを読み出し、遅延推定エコー信号の振幅スペクトルとして出力する。

エコーサプレスゲイン計算部１１７は、近端入力信号の振幅スペクトルと遅延推定エコー信号の振幅スペクトルとを用いて、近端入力信号に重畳されている音響エコー信号を抑圧するエコーサプレスゲインを算出し、算出したエコーサプレスゲインを出力する。

エコーサプレス部１１８は、エコーサプレスゲインと近端入力信号の周波数スペクトルを乗じることにより、近端入力信号に重畳されている音響エコー信号が抑圧した周波数スペクトルを求め、近端出力信号の周波数スペクトルとして出力する。

近端出力信号時間領域変換部１１９は、近端出力信号の周波数スペクトルを、例えば、逆高速フーリエ変換（ＩｎｖｅｒｓｅＦＦＴ）等により、時間領域のデジタル音信号に変換し、近端出力信号として出力する。

遠端信号保持部１２１は、遠端信号の振幅スペクトルを推定エコー信号保持部１１２と同じ所定時間保持する。そして、遠端信号保持部１２１は、推定したフレーム遅延量だけ遅延させた遠端信号の振幅スペクトルを出力する。

近端出力信号振幅スペクトル計算部１２２は、近端出力信号の周波数スペクトルに基づいて、近端出力信号の振幅スペクトルを算出し、算出した近端出力信号の振幅スペクトルを出力する。

シングルトーク判定部１２３は、近端入力信号の振幅スペクトルと近端出力信号の振幅スペクトル等を用いてシングルトークかシングルトーク以外かを判定を行い、シングルトーク判定結果を出力する。

（Ａ−２）第１の実施形態の動作
次に、第１の実施形態に係るエコーサプレス装置１００におけるエコー抑圧処理及び遅延量推定処理の動作を詳細に説明する。

（Ａ−２−１）エコー抑圧処理
まず、エコー抑圧装置１００の動作開始後、例えば、インターネットプロトコル（ＩＰ）網等のネットワークや、携帯電話等の無線ネットワークの電波等に接続されている回線を介して、遠端側の遠端信号が遠端信号入力端子１０１に入力される。

遠端信号入力端子１０１に入力された遠端信号は、出力バッファ１０２において一時保持され、ＤＡ変換器１０３の処理準備ができ次第、ＤＡ変換器１０３に遠端信号を出力される。ＤＡ変換器１０３において、遠端信号はデジタル音信号からアナログ音信号に変換され、スピーカ１０４を通して近端側に出力される。

一方、近端側の話者が発した音声等の音信号や、環境音、音響エコー信号（例えば、スピーカ１０４から出力されたアナログ音信号が近端側の空間を伝達して回り込んだ信号）等が重畳したアナログ音信号は、マイク１０５において受音され、ＡＤ変換器１０５おいてデジタル音信号に変換される。デジタル音信号は、入力バッファ１０７おいて一時保持され、エコー抑圧装置１００の準備ができ次第、デジタル音信号が近端入力信号としてエコー抑圧装置１００に入力される。

遠端信号周波数領域変換部１０８では、例えば、高速フーリエ変換（ＦＦＴ）等により、遠端信号を周波数領域の信号に変換し、変換された遠端信号の周波数スペクトルＲＯＵＴ（ｉ，ω）を遠端信号振幅スペクトル計算部１０９に出力する。

遠端信号振幅スペクトル計算部１０９では、周波数スペクトルＲＯＵＴ（ｉ，ω）を用いて、（１）式に従い、遠端信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜が求められる。

ここで、ｉはフレーム、ωは周波数ビン、ＲＯＵＴ＿ｒｅａｌ（ｉ，ω）とＲＯＵＴ＿ｉｍａｇｅ（ｉ，ω）は、フレームｉにおける周波数ビンωの遠端信号の周波数スペクトルの実数部と虚数部を示しており、遠端信号の周波数スペクトルＲ０ＵＴ（ｉ，ω）は、（２）式で表すことができる。

（２）式のｊは虚数を表している。そして、遠端信号振幅スペクトル計算部１０９により求められた遠端信号の周波数スペクトル｜ＲＯＵＴ（ｉ，ω）｜は、推定エコー信号計算部１１１、遠端信号保持部１２１に出力する。

エコーパス特性更新部１１０は、保持している前フレームに推定したエコーパス特性｜Ｈ（ｉ−１，ω）｜を推定エコー信号計算部１１１に出力する。

推定エコー信号計算部１１１では、保持しているエコーパス特性｜Ｈ（ｉ−１，ω）｜と、遠端信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜とを用いて、（３）式により、推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜が求められる。

（３）式は遠端信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜に、エコーパス保持部１１０に保持しているエコーパス特性｜Ｈ（ｉ−１，ω）｜の対応する周波数ビンを乗じて、当該周波数ビンの推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜を求める。そして、推定エコー信号計算部１０９により求められた推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜を推定エコー信号保持部１１２に出力する。

推定エコー信号保持部１１２では、推定エコー信号計算部１１１から出力された推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜が書き込まれる。

書込みが完了すればフレーム遅延量推定部１１５に推定エコー信号保持部１１２が保持している過去の推定エコー信号の振幅スペクトルを複数フレーム出力する。

一方、近端入力信号周波数領域変換部１１３では、入力バッファ１０７から出力されたデジタル音信号を近端入力信号として、例えば、高速フーリエ変換（ＦＦＴ）等により、近端入力信号を周波数領域の信号に変換する。変換された近端入力信号の周波数スペクトルＳＩＮ（ｉ，ω）は、近端入力信号振幅スペクトル計算部１１４及びエコーサプレス部１１８に出力される。

近端入力信号振幅スペクトル計算部１１４は、近端入力信号の周波数スペクトルＳＩＮ（ｉ，ω）を用いて、（５）式に従い、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜を求める。

ここで、ＳＩＮ＿ｒｅａ１（ｉ，ω）とＳＩＮ＿ｉｍａｇｅ（ｉ，ω）は、フレームｉにおける周波数ビンωの近端入力信号の周波数スペクトルの実数部と虚数部を示しており、近端入力信号の周波数スペクトルＳＩＮ（ｉ，ω）は、式（６）で表すことができる。

（６）式のｊは虚数を表している。そして、近端入力信号振幅スペクトル計算部１１４により求められた近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜は、フレーム遅延量推定部１１５、エコーサプレスゲイン計算部１１７、シングルトーク判定部１２３、及びエコーパス特性更新部１１０に出力される。

フレーム遅延量推定部１１５は、過去の複数フレームの推定エコー信号の振幅スペクトルと、過去の複数フレームの近端入力信号の振幅スペクトルでフレーム遅延量ｄｅｌａｙ（ｉ）を推定し、推定したフレーム遅延量ｄｅｌａｙ（ｉ）を推定エコー信号遅延部１１６、及び遠端信号保持部１２１に出力する。

推定エコー信号補正部１１６では、フレーム遅延量ｄｅｌａｙ（ｉ）を用いて、推定エコー信号保持部１１２に保持されている過去の推定エコー信号の振幅スペクトルを当該フレーム遅延量だけ遅延させたフレームを読出し、遅延した推定エコー信号とする。この遅延した推定エコー信号は、（７−Ａ）式で表すことができる。

遅延した推定エコー信号の振幅スペクトル｜ＥＣＨＯ＿Ｄｅ１ａｙ（ｉ、ω）｜をエコーサプレスゲイン計算部１１７に出力する。

エコーサプレスゲイン計算部１１７では、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜と、遅延した推定エコー信号の振幅スペクトル｜ＥＣＨＯ＿Ｄｅ１ａｙ（ｉ、ω）｜とを取得して、式（７−Ｂ）を用いて、エコーサプレスゲインＧ（ｉ，ω）を求める。

（７−Ｂ）式は、周波数ビン毎に近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜から、遅延した推定エコー信号の振幅スペクトル｜ＥＣＨＯ＿Ｄｅ１ａｙ（ｉ，ω）｜を差し引いた振幅スペクトルを、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜で除することで、エコーサプレスゲインＧ（ｉ，ω）を求めている。エコーサプレスゲイン計算部１１７により求められたエコーサプレスゲインＧ（ｉ，ω）は、エコーサプレス部１１８に出力する。

エコーサプレス部１１８では、近端入力信号のスペクトルＳＩＮ（ｉ，ω）とエコーサプレスゲインＧ（ｉ，ω）とを用いて、（８）式、（９）式に従い、近端入力信号のスペクトルＳＩＮ（ｉ，ω）に重畳されている音響エコー信号を抑圧する。

ここで、ＳＯＵＴ＿ｒｅａｌ（ｉ，ω）とＳＯＵＴ＿ｉｍａｇｅ（ｉ，ω）は、フレームｉにおける周波数ビンωの近端出力信号の周波数スペクトルの実数部と虚数部を示しており、近端出力信号の周波数スペクトルＳＯＵＴ（ｉ，ω）は、（１０）式で表すことができる。

（１０）式のjは虚数を表している。（８）式と（９）式では、周波数スペクトルの実数部、虚数部にエコーサプレスゲインＧ（ｉ，ω）を周波数ビン毎に乗じて、音響エコー信号を抑圧した近端出力信号の周波数スペクトルを求める。そして、エコーサプレス部１１８により求められた音響エコー信号が抑圧された近端出力信号の周波数スペクトルＳＯＵＴ（ｉ，ω）を近端出力信号時間領域変換部１１９、及び近端出力信号振幅スペクトル計算部１２２に出力する。

近端出力信号時間領域変換部１１９では、近端出力信号のスペクトルＳＯＵＴ（ｉ，ω）を、例えば逆高速フーリエ変換（ＩｎｖｅｒｓｅＦＦＴ）等により、時間領域のデジタル音信号に変換し、近端出力信号を近端信号出力端子１２０に出力する。

近端信号出力端子１２０は、例えば、インターネットプロトコル（ＩＰ）網等のネットワークや、携帯電話等の無線ネットワークの電波等に接続されており、近端出力信号を接続されている回線を介して通話相手である遠端側に出力する。

遠端信号保持部１２１は、遠端信号振幅スペクトル計算部１０９から出力された遠端信号の振幅スペクトル｜ＲＩＮ（ｉ，ω）｜を、推定エコー信号保持部１１２と同じ書込みフレームに書込む。

そして、遠端信号保持部１２１は、（１２）式に従って、フレーム遅延量ｄｅｌａｙ（ｉ）だけ遅延させた遠端信号の振幅スペクトル｜ＲＯＵＴ＿Ｄｅ１ａｙ（ｉ，ω）｜をエコーパス特性更新部１１０に出力する。

近端出力信号振幅スペクトル計算部１２２では、近端出力信号の周波数スペクトルＳＯＵＴ（ｉ，ω）を用いて、（１３）式に従い、近端出力信号の振幅スペクトル｜ＳＯＵＴ（ｉ，ω）｜が求められる。

そして、近端出力信号振幅スペクトル計算部１２２により求められた近端入力信号の振幅スペクトル｜ＳＯＵＴ（ｉ，ω）｜は、シングルトーク判定部１２３に出力される。

シングルトーク判定部１２３では、近端入力信号がシングルトークか又はシングルトーク以外かを、近端入力信号の振幅スペクトルと近端出力信号の振幅スペクトルとを用いて判定する。シングルトークか又はシングルトーク以外かを判定する手法は、例えば、（１４）式に従い、判定する手法がある。

（１４）式の Ｆｓはサンプリング周波数、ＴＨ１は閾値である。つまり、（１４）式の条件が真のときはシングルトークと判定し、偽のときはシングルトーク以外として判定する。閾値ＴＨ１は、（１４）式の場合、シングルトーク時は（１４）式の左辺が小さい値になるので、小さい固定値（例えばＴＨ１＝０．３等）やフレームで変化する変数などにしても良い。なお、シングルトークか又はシングルトーク以外かの判定の手段は、種々の方法を広く適用することができ、例えば、推定したフレーム遅延量だけ遅延させた遠端信号の振幅スペクトルと、近端信号の振幅スペクトルとの相関を求め、その相関値が高いときはシングルトークとする方法で判定しても良い。シングルトーク判定部１２３は、シングルトーク判定結果をエコーパス特性更新部１１０に出力する。

エコーパス特性更新部１１０では、シングルトーク判定部１２３でシングルトークと判定されたフレームで（１５）式に従い、現フレームの推定エコーパス特性｜Ｈ１（ｉ，ω）｜を遠端信号保持部１２１からの遅延させた遠端信号の振幅スペクトル｜ＲＩＮ＿Ｄｅ１ａｙ（ｉ，ω）｜と、近端出力信号振幅スペクトル計算部１１４からの近端出力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜を用いて求める。

現フレームのエコーパス特性｜Ｈ１（ｉ，ω）｜が求まれば、エコーパス特性更新部１１０に保持されている１フレーム前の推定エコーパス特性｜Ｈ（ｉ−１，ω）｜を読み出し、｜Ｈ（ｉ−１，ω）｜と、｜Ｈ１（ｉ，ω）｜とを用いてエコーパス特性を更新する。具体的には、エコーパス特性更新部１１０は、（１６）式に従って、エコーパス特性｜Ｈ（ｉ，ω）｜を更新する。

（１６）式において、ａは時定数フィルタの係数であり、ａは、０以上１より小さい値であって、エコーパス特性の更新を遅くしたい場合、ａは１に近い値が望ましく（例えばａ＝０．９９等の値）、更新を早くしたい場合、ａは０に近い値が望ましい（例えばａ＝０．０１等の値）。更新した推定エコーパス特性はエコーパス特性更新部１１０に保持される。

（Ａ−２−２）フレーム遅延量推定処理
次に、第１の実施形態に係るフレーム遅延量推定部１１５におけるフレーム遅延量推定処理の動作を、図面を参照しながら詳細に説明する。

図２は、第１の実施形態に係るエコー抑圧装置１００のフレーム遅延量推定部１１５の処理フローを示す処理フロー図である。

図２において、本発明のエコー抑圧装置１００のフレーム遅延量推定部１１５は、フレーム誤差計算２０１、フレーム遅延量候補計算２０２、フレーム遅延量推定２０３を行う。

フレーム誤差計算２０１では、近端入力信号の振幅スペクトルと過去の複数の推定エコー信号の振幅スペクトルとの誤差をフレーム毎に計算し、フレーム毎の誤差を出力する。

フレーム遅延量候補値計算２０２では、フレーム毎の誤差から、フレーム遅延量の候補値を計算して出力する。

フレーム遅延量推定２０３では、入力されたフレーム遅延量の候補値から推定エコー信号遅延部と遠端信号保持部１２１で使用するフレーム遅延量を推定し、推定したフレーム遅延量を出力する。

（Ａ−２−３）フレーム遅延量推定の詳細な処理
フレーム誤差計算２０１において、フレーム遅延量推定部１１５は、推定エコー信号保持部１１２から出力される過去の推定エコー信号の振幅スペクトルの複数フレームＥＣＨＯ＿Ｂｕｆｆｅｒ（ｎ−ｋ，ω）と近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜とのフレーム誤差Ｅｒｒｏｒ（ｋ）がフレーム毎に（１７）式に従い求められる。

Ｆｓはサンプリング周波数、Ｎは探索フレーム数である。そして、計算されたフレーム誤差Ｅｒｒｏｒ（ｋ）をフレーム遅延量候補値計算２０２に出力される。

フレーム遅延量候補値計算２０２では、フレーム誤差計算２０１で求められたフレーム誤差Ｅｒｒｏｒ（ｋ）に基づきフレーム遅延量の候補値が計算される。フレーム遅延量の候補値の計算手法は、例えば、（１８）式と（１９）式に従って、フレーム誤差Ｅｒｒｏｒ（ｋ）が最小になるフレーム遅延量と、２番目に最小になるフレーム遅延量の候補値とする。

（１８）式と（１９）式のａｒｇｍｉｎ（ｘ（ｋ））という関数は、ｘ（ｋ）が最小になる引数を出力する関数であり、（１８）式はフレーム誤差Ｅｒｒｏｒ（ｋ）が最小になるフレーム遅延量を出力し、ｄｅｌａｙｌ（ｉ）に代入するという式であり、（１９）式はフレーム誤差Ｅｒｒｏｒ（ｋ）が2番目に小さくになるフレーム遅延量を出力し、ｄｅｌａｙ２（ｉ）に代入するという式である。なお、フレーム遅延量の候補値計算の手法は、種々の方法を広く適用することができ、例えば、現フレームまでのフレーム遅延量の候補値の平均値を計算し、フレーム誤差Ｅｒｒｏｒ（ｋ）が小さく平均値に近い値をフレーム遅延量の候補値とする方法でも良い。フレーム遅延量候補値計算２０２では、計算したフレーム遅延量の候補値ｄｅｌａｙ１（ｉ）、ｄｅｌａｙ２（ｉ）をフレーム遅延量推定２０３に出力する。

フレーム遅延量推定２０３では、フレーム遅延量の候補値から推定エコー信号の算出やエコーパス特性の算出に使用するフレーム遅延量を推定し出力する。

図４は、第１の実施形態に係るフレーム遅延量を推定する方法を説明する説明図である。フレーム遅延量の推定方法は、図４に示すようにフレーム遅延量の候補値であるフレームからＬフレーム過去のフレームの推定エコー信号の振幅スペクトルと、Ｌフレーム過去のフレームの近端入力信号の振幅スペクトルとのフレーム誤差を、（２０）式と（２１）式に従い求める。

（２０）式と（２１）式は、フレーム遅延量の候補値からＬフレーム過去のフレームのフレーム誤差を計算するという式である。図４において、フレーム誤差が最小となるフレーム遅延量の候補値（図４では遅延量候補値１と表記）と、２番目にフレーム誤差が小さくなるフレーム遅延量の候補値（図４では遅延量候補値２と表記）としている。

この場合、フレーム遅延量候補値１のフレームからＬフレーム過去のフレーム（実線部分）の推定エコー信号の振幅スペクトルと、現フレームからＬフレーム過去のフレームの近端入力信号の振幅スペクトルとのフレーム誤差Ｅｒｒｏｒ１を計算する。同様に、遅延量候補値２のフレームからＬフレーム過去フレーム（実線部分）の推定エコー信号の振幅スペクトルと、現フレームからＬフレーム過去のフレームの近端入力信号の振幅スペクトルとの誤差を計算する。

ここで、Ｌは、例えば、約１音素分の長さが望ましい。つまり、フレーム遅延量の候補値のフレームからＬフレーム離れたフレーム（換言すると、非連続のフレーム）である。これは、フレーム遅延量の候補値から所定長（例えば１音素分の長さ程度）離れたフレームの振幅スペクトルの誤差を見ることでフレーム遅延量を精度良く推定することができるからである。

遅延量候補値１のフレームからＬフレーム過去のフレームの推定エコー信号の振幅スペクトルと現フレームからＬフレーム過去のフレームの近端入力信号の振幅スペクトルとのフレーム誤差Ｅｒｒｏｒ１と、遅延量候補値２のフレームからＬフレーム過去のフレームの推定エコー信号の振幅スペクトルと現フレームからＬフレーム過去のフレームの近端入力信号の振幅スペクトルとのフレーム誤差Ｅｒｒｏｒ２とを用いて、（２２）式でフレーム遅延量を決定する。

（２２）式は、Ｅｒｒｏｒ１と、Ｅｒｒｏｒ２とを比較し、Ｅｒｒｏｒ１＜Ｅｒｏｒｒ２の場合、遅延量候補値１をフレーム遅延量として決定し、そうでない場合、遅延量候補値２をフレーム遅延量として決定する。これは、Ｌフレームだけ過去のフレームの近端入力信号の振幅スペクトルと推定エコー信号の振幅スペクトルとの誤差を求めることで、より高い確度でフレーム遅延量の候補値を決定し、精度の高いフレーム遅延量を決定している。

そして、フレーム遅延量推定部２０３は、決定したフレーム遅延量を推定エコー信号遅延部１１６と遠端信号保持部１２１に出力する。

（Ａ−３）第１の実施形態の効果
以上のように、第１の実施形態によれば、パルス音や報知音などの音を出力しなくてもフレーム遅延量を推定でき、音響エコー信号を抑圧することができ、エコー抑圧性能を維持する。

（Ｂ）第２の実施形態
次に、本発明に係るエコー抑圧装置、エコー抑圧プログラム及びエコー抑圧方法の第２の実施形態を、図面を参照しながら詳細に説明する。

（Ｂ−１）第２の実施形態の構成
第２の実施形態では、フレーム遅延量推定部１１５の処理が、第１の実施形態に係るエコー抑圧装置１００と異なるのみであり、それ以外の構成要素は第１の実施形態に係る図１のエコー抑圧装置１００の構成要素と同一又は対応するものである。

（Ｂ−２）第２の実施形態の動作
（Ｂ−２−１）フレーム遅延量推定処理
第２の実施形態のフレーム遅延量推定部の処理フローを、図面を参照しながら詳細に説明する。

図３は、第２の実施形態に係るエコー抑圧装置２００のフレーム遅延量推定部１１５の処理フローを示す処理フロー図である。

図３において、第２の実施形態に係るエコー抑圧装置２００のフレーム遅延量推定部１１５は、フレーム誤差計算２０１、フレーム遅延量候補計算２０２、フレーム遅延量推定判定３０１、フレーム遅延量推定２０３、フレーム遅延量決定３０２を行う。

フレーム遅延量推定判定３０１は、入力されたフレーム遅延量の候補値からフレーム遅延量を推定するかしないかを判定する。

フレーム遅延量決定３０２は、フレーム遅延量推定判定３０１による結果に基づいて、推定しないと判定すると近端入力信号の振幅スペクトルと当該フレーム遅延量の候補値のフレーム誤差に基づいてフレーム遅延量の候補値から推定エコー信号の算出やエコーパス特性の算出に使用するフレーム遅延量を決定し、決定したフレーム遅延量を出力する。

（Ｂ−２−２）フレーム遅延量推定処理の詳細な動作
フレーム遅延量推定判定３０１では、入力されたフレーム遅延量の候補値からフレーム遅延量を推定するかしないかを判定する。フレーム遅延量を推定するかしないかを判定する手法は、例えば、（２１）式に従い、判定する。

（２３）式のＴＨ２は閾値である。（２３）式の条件が真のとき、フレーム遅延量の候補値が近い値であることを意味しているので、フレーム遅延量の推定は行わないと判定する。一方、（２３）式の条件が偽のとき、フレーム遅延量の候補値がまったく異なる値であることを意味しているので、どちらを用いるかフレーム遅延量を推定する必要があると判定する。閾値ＴＨ２は固定閾値（例えばＴＨ２＝５やＴＨ２＝Ｌなど）としても良く、変動閾値（例えば、現フレームまでの推定したフレーム遅延量の推定値の平均値）としても良い。なお、フレーム遅延量を推定するかしないかを判定の手段は、種々の方法を広く適用することができ、例えば、現フレームまで推定したフレーム遅延量の平均値を求め、その平均値にフレーム遅延量の候補値が近いかどうかを判定し、両方とも近い値になったときはフレーム遅延量を推定しないと判定し、どちらか一方、または両方ともまったく異なる値になったときはフレーム遅延量を推定する方法で判定しでも良い。そして、フレーム遅延量推定判定３０１は判定結果に基づいてフレーム遅延量推定２０３かフレーム遅延量決定３０２を行う。

フレーム遅延量決定３０２は、フレーム遅延量推定判定３０１でフレーム遅延量を推定しないと判定されると、フレーム遅延量の候補値から推定エコー信号の算出やエコーパス特性の算出に使用するフレーム遅延量を決定し出力する。フレーム遅延量を決定する手法は、例えば、フレーム誤差計算部２０１で算出した当該フレーム遅延量の候補値のフレーム誤差に基づいて（２４）式に従い決定する手法がある。

（２２）式は、当該フレーム遅延量の候補値の誤差が小さいフレーム遅延量を推定エコー信号遅延部と遠端信号保持部１２１で使用するフレーム遅延量とするという式である。フレーム遅延量推定判定３０１でフレーム遅延量を推定しないと判定されるということは、当該フレーム遅延量の候補値の誤差が小さい値になっているということなので、フレーム誤差Ｅｒｒｏｒ（ｋ）が小さいほうをフレーム遅延量としている。なお、フレーム遅延量を決定する手法は、種々の方法を広く適用することができ、例えば，当該フレーム遅延量の候補値の平均値を計算しフレーム遅延量とする方法で判定しでも良い。そして、フレーム遅延量決定３０２は、決定したフレーム遅延量を推定エコー信号遅延部１１６と遠端信号保持部１２１に出力する。

フレーム遅延量推定２０３は、フレーム遅延量推定判定３０１でフレーム遅延量を推定すると判定されると、第１の実施形態と同様にフレーム遅延量の候補値から推定エコー信号の算出やエコーパス特性の算出に使用するフレーム遅延量を推定し出力する。

（Ｂ−３）第２の実施形態の効果
以上のように、第２の実施形態によれば、パルス音や報知音などの音を出力しなくてもフレーム遅延量を推定でき、フレーム遅延量の候補値が大きく異なっているときにのみ、フレーム遅延量を推定し、音響エコー信号を抑圧することができ、エコー抑圧性能を維持する。

（Ｃ）第３の実施形態
次に、本発明に係るエコー抑圧装置、エコー抑圧プログラム、及びエコー抑圧方法の第３の実施形態を、図面を参照しながら詳細に説明する。

（Ｃ−１）第３の実施形態の構成
第３の実施形態では、フレーム遅延量推定部１１５のフレーム遅延量推定２０３が、第１及び第２の実施形態に係るエコー抑圧装置１００及び２００と異なるのみであり、それ以外の構成要素は、第１及び第２の実施形態に係る図１のエコー抑圧装置１００及び２００の構成要素と同一又は対応するものである。

（Ｃ−２）第３の実施形態の動作
（Ｃ−２−１）フレーム遅延量推定処理
第３の実施形態のフレーム遅延量推定部１１５の処理フローを、図面を参照しながら詳細に説明する。

図５は、第３の実施形態に係るエコー抑圧装置３００のフレーム遅延量推定部１１５のフレーム処理量推定を示す図である。

図５において、本発明のエコー抑圧装置３００のフレーム遅延量推定部１１５のフレーム遅延量推定は、フレーム遅延量の候補値から過去の複数フレームの推定エコー信号と過去の複数フレームの近端入力信号を用いて入力されたフレーム遅延量の候補値からフレーム遅延量を推定し、推定したフレーム遅延量を出力する。

（Ｃ−２−２）フレーム遅延量推定処理の詳細な動作
フレーム遅延量推定２０３は、フレーム遅延量の候補値から推定エコー信号の算出やエコーパス特性の算出に使用するフレーム遅延量を推定し出力する。推定方法は、図５に示すようにフレーム遅延量の候補値であるフレームからＬｍフレーム過去のフレームの推定エコー信号と過去の近端入力信号の振幅スペクトルとのフレーム誤差を式（２５）と（２６）に従い求める。

（２５）式と（２６）式はフレーム遅延量の候補値からＬｍフレーム過去のフレームのフレーム誤差を計算するという式である。Ｌｍは、例えば、約ｍ音素分の長さ（すなわち、Ｌｍ＝ｍ×Ｌ）が望ましい。フレーム遅延量の候補値であるフレームからＬｍフレーム過去のフレーム誤差Ｅｒｒｏｒ１（ｍ）とＥｒｒｏｒ２（ｍ）を用いて、例えば（２７）式でフレーム遅延量を決定する。

（２７）式はフレーム遅延量の候補値のうちフレーム誤差Ｅｒｒｏｒ１（ｍ）とＥｒｒｏｒ２（ｍ）の合計値が小さいほうをフレーム遅延量として使用するという式である。なお、フレーム遅延量の推定の手法は、種々の方法を広く適用することができ、例えば、フレーム誤差Ｅｒｒｏｒ１（ｍ）とＥｒｒｏｒ２（ｍ）を１≦ｍ≦Ｍで比較し、小さい回数が多かったほうをフレーム遅延量とするとする方法で判定しでも良い。そして、フレーム遅延量推定部２０３は、決定したフレーム遅延量を推定エコー信号遅延部１１６と遠端信号保持部１２１に出力する。

（Ｃ−３）第３の実施形態の効果
以上のように、第３の実施形態によれば、パルス音や報知音などの音を出力しなくてもフレーム遅延量を推定でき、複数のフレームでフレーム遅延量を推定することで、音響エコー信号を抑圧することができ、エコー抑圧性能を維持する。

（Ｄ）他の実施形態
上述した各実施形態においても、種々の変形実施形態を説明したが、本発明は以下の変形実施形態についても適用することができる。

（Ｄ−１）上述した各実施形態で説明したエコー抑圧装置は、例えば、テレビ会議システムや電話会議システム等に用いられる音声通信装置を含む装置に搭載されるようにしても良い。また、携帯電話機やスマートフォン等の携帯端末に本発明のエコー抑圧装置は搭載されるようにしても良い。

（Ｄ−２）上述した第１の実施形態では、フレーム遅延量推定部１１５が、候補値のフレームからＬフレーム過去のフレームの、推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとの誤差を求めて、フレーム遅延量を決定する場合を例示した。しかし、近端入力信号を一時的に保持することができるのであれば、候補値のフレームからＬフレーム未来のフレームの、推定エコー信号の振幅スペクトルと、現フレームからＬフレーム未来のフレームの近端入力信号の振幅スペクトルとの誤差に基づいて、フレーム遅延量を決定するようにしてもよい。

（Ｄ−３）上述した第１〜第３の実施形態では、候補値のフレームからＬフレーム過去のフレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとの誤差を、所定の演算式を用いて求めて、候補値を決定してフレーム遅延量を求める場合を例示した。しかし、候補値のフレームからＬフレーム過去のフレームの推定エコー信号の振幅スペクトルと、現フレームの近端入力信号の振幅スペクトルとの相関を求めて遅延量の候補値を決定するようにしても良い。

（Ｄ−４）上述した第１〜第３の実施形態では、候補値のフレームからＬフレーム過去のフレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとの誤差を求めるようにしたが、近端入力信号を一時的に保持することができるのであれば、候補値のフレームからＬフレーム未来のフレームの推定エコー信号の振幅スペクトルと、現フレームからＬフレーム未来のフレームの近端入力信号の振幅スペクトルとの誤差を求めるようにしても良い。

１００…エコー抑圧装置、１０１…遠端信号入力端子、１０２…出力バッファ１０３…ＤＡ変換器、１０４…スピーカ、１０５…マイク、１０６…ＡＤ変換器、１０７…入力バッファ、１０８…遠端信号周波数領域変換算部、１０９…遠端信号振幅スペクトル計算部、１１０…エコーパス特性更新部、１１１…推定エコー信号計算部、１１２…推定エコー信号保持部、１１３…近端入力信号周波数領域変換部、１１４…近端入力信号振幅スペクトル計算部、１１５フレーム遅延量推定部、１１６推定エコー信号遅延部、１１７…エコーサプレスゲイン計算部、１１８…エコーサプレス部、１１９…近端出力信号時間領域変換部、１２０…近端信号出力端子、１２１…遠端信号保持部、１２２…近端出力信号振幅スペクトル計算部、１２３…シングルトーク判定部、２０１…フレーム誤差計算、２０２…フレーム遅延量候補値計算、２０３…フレーム遅延量推定、３０１…フレーム遅延量推定判定、３０２…フレーム遅延量決定部、４００…エコーキャンセラ、４０１…スピーカ、４０２…マイク、４０３…適応フィルタ部、４０４…エコー抑圧処理部、４０５…遅延量推定部、４０６…遅延処理部。

Claims (5)

1. 近端入力信号と遠端信号に基づく音響エコー信号を抑圧するエコー抑圧装置において、
   過去の複数フレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとに基づいてフレーム誤差を計算するフレーム誤差計算手段と、
   上記フレーム誤差計算手段で計算したフレーム毎の誤差に基づいて、１又は複数のフレーム遅延量の候補値を計算するフレーム遅延量候補値計算手段と、
   上記フレーム遅延量候補値計算手段で判定された上記各フレーム遅延量の候補値に基づいて、上記遅延量候補値から所定フレーム離れたフレームの推定エコー信号と現フレームから所定フレーム離れたフレームの近端入力信号とを用いて、フレーム遅延量を推定して出力するフレーム遅延量推定手段と
   を有することを特徴とするエコー抑圧装置。
2. 上記フレーム遅延量推定手段が、上記各フレーム遅延量の候補値を用いて、フレーム遅延量を推定するか推定しないかを判定し、フレーム遅延量を推定しないと判定されたとき上記フレーム遅延量の候補値のうち上記フレーム誤差計算手段で計算したフ上記フレーム誤差が小さいほうをフレーム遅延量として決定して出力し、フレーム遅延量を推定する判定されたとき上記フレーム遅延量候補値計算手段でフレーム遅延量の候補値を用いてフレーム遅延量を推定して出力することを特徴とする請求項１に記載のエコー抑圧装置。
3. 上記フレーム遅延量推定手段が、フレーム遅延量を推定するとき上記フレーム遅延量候補値計算手段によるフレーム遅延量の候補値に基づいて、過去又は未来の複数フレームの推定エコー信号と、過去又は未来の複数フレームの近端入力信号とを用いてフレーム遅延量を推定して出力することを特徴とする請求項１又は２に記載のエコー抑圧装置。
4. 近端入力信号から遠端出力信号に基づく音響エコー信号を抑圧するエコー抑圧プログラムにおいて、
   コンピュータを、
   過去の複数フレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとに基づいてフレーム誤差を計算するフレーム誤差計算手段と、
   上記フレーム誤差計算手段で計算したフレーム毎の誤差に基づいて、１又は複数のフレーム遅延量の候補値を計算するフレーム遅延量候補値計算手段と、
   上記フレーム遅延量候補値計算手段で判定された上記各フレーム遅延量の候補値に基づいて、上記遅延量候補値から所定フレーム離れたフレームの推定エコー信号と現フレームから所定フレーム離れたフレームの近端入力信号とを用いて、フレーム遅延量を推定して出力するフレーム遅延量推定手段と
   して機能させることを特徴とするエコー抑圧プログラム。
5. 近端入力信号から遠端出力信号に基づく音響エコー信号を抑圧するエコー抑圧方法において、
   フレーム誤差計算手段、フレーム遅延量候補値計算手段、及びフレーム遅延量推定手段を有し、
   上記フレーム誤差計算手段は、過去の複数フレームの推定エコー信号の振幅スペクトルと、近端入力信号の振幅スペクトルとに基づいてフレーム誤差を計算し、
   上記フレーム遅延量候補値計算手段は、上記フレーム誤差計算手段で計算したフレーム毎の誤差に基づいて、１又は複数のフレーム遅延量の候補値を計算し、
   上記フレーム遅延量推定手段は、上記フレーム遅延量候補値計算手段で判定された上記各フレーム遅延量の候補値に基づいて、上記遅延量候補値から所定フレーム離れたフレームの推定エコー信号と現フレームから所定フレーム離れたフレームの近端入力信号とを用いて、フレーム遅延量を推定して出力する
   ことを特徴とするエコー抑圧方法。
   

Images