(57)【要約】【課題】遠端出力信号の周波数成分が小さいときのエコーパス特性の誤学習を防ぎ、エコー抑圧性能を維持する。【解決手段】本発明は、遠端出力信号の振幅スペクトルと閾値とを周波数ビン毎に比較して、周波数ビン毎に周波数成分の有無を判定する周波数ビン成分有無判定手段と、周波数ビン成分有無判定手段の判定結果により周波数成分があると判断された周波数ビンの遠端出力信号の振幅スペクトルと、近端入力信号の対応する周波数ビンの振幅スペクトルとに基づいて、当該周波数ビンのエコーパス特性を推定する周波数ビンエコーパス特性推定手段と、周波数ビン毎のエコーパス特性と遠端入力信号の振幅スペクトルとに基づいて推定エコー信号を算出し、近端入力信号の振幅スペクトルから推定エコー信号を抑圧するエコー抑圧装置である。【選択図】      図１
2015070290.tif 000002
【特許請求の範囲】
【請求項１】
入力された遠端出力信号を周波数領域に変換して、上記遠端出力信号の振幅スペクトルを求める遠端出力信号振幅スペクトル算出手段と、
入力された近端入力信号を周波数領域に変換して、上記近端入力信号の振幅スペクトルを求める近端入力信号振幅スペクトル算出手段と、
上記遠端出力信号の振幅スペクトルと閾値とを周波数ビン毎に比較して、周波数ビン毎に周波数成分の有無を判定する周波数ビン成分有無判定手段と、
上記周波数ビン成分有無判定手段の判定結果により周波数成分があると判断された周波数ビンの上記遠端出力信号の振幅スペクトルと、上記近端入力信号の振幅スペクトルの対応する周波数ビンとに基づいて、当該周波数ビンのエコーパス特性を推定する周波数ビンエコーパス特性推定手段と、
上記周波数ビンエコーパス特性推定手段により推定された周波数ビン毎のエコーパス特性と上記遠端出力信号の振幅スペクトルとに基づいて推定エコー信号を算出する推定エコー信号算出手段と、
上記近端入力信号の振幅スペクトルから上記推定エコー信号を抑圧するエコー抑圧手段と
を備えることを特徴とするエコー抑圧装置。
【請求項２】
上記エコーパス特性を周波数ビン毎に保持するエコーパス特性保持手段を備え、
上記周波数ビンエコーパス特性推定手段が、上記エコーパス特性保持手段に保持されるエコーパス特性を上記周波数ビン毎に更新するものである
ことを特徴とする請求項１に記載のエコー抑圧装置。
【請求項３】
上記周波数ビンエコーパス特性推定手段が、上記遠端出力信号の振幅スペクトルの周波数成分があると判断された周波数ビンについて上記エコーパス特性を更新し、上記遠端出力信号の振幅スペクトルの周波数成分がないと判断された周波数ビンについて上記エコーパス特性を更新しないことを特徴とする請求項２に記載のエコー抑圧装置。
【請求項４】
上記周波数ビンエコーパス特性推定手段が、上記遠端出力信号の振幅スペクトルの周波数成分がある周波数ビンと、その次の上記遠端出力信号の振幅スペクトルの周波数成分がある周波数ビンとの間であって、上記遠端出力信号の振幅スペクトルの成分がないと判断された周波数ビンの上記エコーパス特性を補間し、更新することを特徴とする請求項２又は３に記載のエコー抑圧装置。
【請求項５】
コンピュータを、
入力された遠端出力信号を周波数領域に変換して、上記遠端出力信号の振幅スペクトルを求める遠端出力信号振幅スペクトル算出手段と、
入力された近端入力信号を周波数領域に変換して、上記近端入力信号の振幅スペクトルを求める近端入力信号振幅スペクトル算出手段と、
上記遠端出力信号の振幅スペクトルと閾値とを周波数ビン毎とに比較して、周波数ビン毎に周波数成分の有無を判定する周波数ビン成分有無判定手段と、
上記周波数ビン成分有無判定手段の判定結果により周波数成分があると判断された上記遠端出力信号の振幅スペクトルと、上記近端入力信号の対応する周波数ビンの振幅スペクトルに基づいて、当該周波数ビンのエコーパス特性を推定する周波数ビンエコーパス特性推定手段と、
上記周波数ビンエコーパス特性推定手段により推定された周波数ビン毎のエコーパス特性と上記遠端出力信号の振幅スペクトルとに基づいて推定エコー信号を算出する推定エコー信号算出手段と、
上記近端入力信号の振幅スペクトルから上記推定エコー信号を抑圧するエコー抑圧手段と
して機能させることを特徴とするエコー抑圧プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エコー抑圧装置及びエコー抑圧プログラムに関し、例えば、テレビ会議システムや電話会議システム等において用いられるエコー抑圧装置及びエコー抑圧プログラムに適用し得るものである。
【背景技術】
【０００２】
例えば、テレビ会議システムや電話会議システム等の拡声通話システムでは、スピーカから放音された音（ここで、音は音響や音声等を含む。）がマイクに回り込んで送話側に戻る音響エコー信号が発生する。音響エコー信号は、通話の著しい妨げとなるため、音響エコー抑圧方法に関して、これまでも多くの研究、開発が行なわれている。
【０００３】
音響エコー抑圧方法の１つとして、エコーサプレッサー（エコー抑圧装置）を使用する手法がある。エコーサプレッサーとは、遠端入力信号と近端入力信号とからエコーパス特性、推定エコー信号、エコーサプレスゲインを求めて、近端入力信号とエコーサプレスゲインを乗算することで音響エコー信号を抑圧する手法である。このエコーサプレッサーを使用して音響エコー信号を抑圧するエコー抑圧装置が非特許文献１によって提案されている。
【０００４】
非特許文献１に開示されるエコー抑圧装置は、過去のフレームの遠端出力信号と近端入力信号とに基づいてエコーパス特性を求める。そして、求めたエコーパス特性と遠端出力信号とを乗算した信号を推定エコー信号とし、エコー抑圧装置は、近端入力信号と推定エコー信号とに基づいてエコーサプレスゲインを求め、エコーサプレスゲインを近端入力信号と乗算してエコーを抑圧する。
【０００５】
非特許文献１に記載のエコー抑圧装置について、図４を参照しながら説明する。図４は、非特許文献１に記載の従来のエコー抑圧装置を説明するためのブロック図である。
【０００６】
非特許文献１に記載のエコー抑圧装置４００は、遠端入力信号入力部４０１、デジタル−アナログ（ＤＡ）変換器４０２、スピーカ４０３、マイク４０４、アナログ−デジタル（ＡＤ）変換器４０５、第１周波数分析部４０６、第２周波数分析部４０７、音響結合量計算部４０８、推定エコー信号計算部４０９、ゲイン計算部４１０、積算部４１１、周波数−時間変換部４１２、近端出力信号出力部４１３を有する。
【０００７】
遠端入力部４０１は、通話先から、例えば、ネットワークや電波を介してエコー抑圧装置４００に入力される遠端側（相手側）のデジタル音信号が遠端入力信号として入力され、遠端出力信号としてデジタル−アナログ（ＤＡ）変換器４０２と第１周波数分析部４０６に出力する。ＤＡ変換器４０２は、遠端出力信号として入力されたデジタル音信号をアナログ音信号に変換してスピーカ４０３に出力する。スピーカ４０３は、入力されたアナログ音信号を近端側（自分側）に出力する。
【０００８】
マイク４０４は、近端側の話者が発した音声等の音信号や環境音や音響エコー信号（スピーカ４０３から出力された音信号が近端側の空間を伝達しマイクに入力される音）等の音信号が重畳して入力され、これらをアナログ音信号に変換し、ＡＤ変換器４０５に出力する。ＡＤ変換器４０５は、入力されたアナログ音信号をデジタル音信号に変換し、近端入力信号として第２周波数分析部４０７に出力する。
【０００９】
第１周波数分析部４０６は、遠端入力端子４０１に遠端出力信号として入力されたデジタル音信号を入力として、そのデジタル音信号を周波数スペクトルに変換し遠端入力信号のスペクトルとして出力する。第２周波数分析部４０７は、ＡＤ変換器４０５から近端入力信号として出力されたデジタル音信号を入力として、そのデジタル音信号を周波数スペクトルに変換し近端入力信号のスペクトルとして出力する。音響結合量計算部４０８は、遠端出力信号のスペクトルと近端入力信号のスペクトルとを入力として、エコーパス特性の推定値を出力する。推定エコー信号計算部４０９は、遠端出力信号のスペクトルとエコーパス特性の推定値とに基づいて、推定エコー信号を計算する。ゲイン計算部４１０は、推定エコー信号と近端入力信号のスペクトルを入力として、サプレスゲイン係数を計算する。
【００１０】
積算部４１１は、近端入力信号のスペクトルにサプレスゲイン係数を乗じて近端入力信号からエコーを抑圧したスペクトルを出力する。周波数−時間変換部４１２は、近端入力信号から音響エコー信号を抑圧したスペクトルを時間領域の信号に変換し、近端出力信号（例えば、ネットワークや電波を介して出力される音信号）として近端出力信号出力部４１３に出力される。
【先行技術文献】
【非特許文献】
【００１１】
【非特許文献１】Ｃ．Ｆａｌ１ｅｒ  ａｎｄ  Ｃ．Ｔｕｒｎｅｒｙ，“ＥＳＴＩＭＡＴＩＯＮ  ＴＨＥ  ＤＥＬＡＹ  ＡＮＤ  ＣＯＬＯＲＡＴＩＯＮ  ＥＦＦＥＣＴ  ＯＦ  ＴＨＥ  ＡＣＯＵＳＴＩＣ  ＥＣＨＯ  ＰＡＴＨ  ＦＯＲ  ＬＯＷ  ＣＯＭＰＬＥＸＩＴＹ  ＥＣＨＯ  ＳＵＰＰＲＥＳＳＩＯＮ”，Ｐｒｏｃ．Ｉｎｔｌ．Ｗｏｒｋｓ．ｏｎ  Ａｃｏｕｓｔ．Ｅｃｈｏ  ａｎｄ  Ｎｏｉｓｅ  Ｃｏｎｔｒｏｌ（ＩＷＡＥＮＣ）２００５，ｐｐ．５３−５６，Ｏｃｔ．２００５
【発明の概要】
【発明が解決しようとする課題】
【００１２】
しかしながら、非特許文献１に記載のエコー抑圧装置では、遠端側から入力されるデジタル音信号には、音声信号の微細構造の谷になる周波数ビンや、音信号のスペクトル崩落の形状などで遠端入力信号の周波数成分が小さい周波数ビンが存在する。
【００１３】
そのため、遠端入力信号の小さい周波数ビンの周波数成分が、背景ノイズの同じ周波数ビンの周波数成分に埋もれてしまい、その周波数ビンのエコーパス特性を求めても、実際のエコーパス特性とは全く異なる値が算出されてエコーパス特性を更新することになる。このことから、推定されたエコーパス特性が実際のものと異なるエコーパス特性になり、音響エコー信号を適切に抑圧することができないという問題が生じ得る。
【００１４】
そこで、本発明は、上記課題に鑑みてなされたものであり、遠端入力信号の周波数成分の小さい周波数ビンについては、近端入力信号と遠端入力信号とを用いてエコーパス特性を算出せずに、過去のエコーパス特性を用いてエコーパス特性を更新することで、エコーパス特性の誤学習を防ぎ、エコー抑圧性能を維持することができるエコー抑圧装置及びエコー抑圧プログラムを提供しようとするものである。
【課題を解決するための手段】
【００１５】
かかる課題を解決するために、第１の本発明は、（１）入力された遠端出力信号を周波数スペクトルに変換して、遠端出力信号の振幅スペクトルを求める遠端出力信号振幅スペクトル算出手段と、（２）入力された近端入力信号を周波数スペクトルに変換して、近端入力信号の振幅スペクトルを求める近端入力信号振幅スペクトル算出手段と、（３）遠端出力信号の振幅スペクトルと閾値とを周波数ビン毎に比較して、周波数ビン毎に周波数成分の有無を判定する周波数ビン成分有無判定手段と、（４）周波数ビン成分有無判定手段の判定結果により遠端出力信号の周波数成分があると判断された周波数ビンの振幅スペクトルと、近端入力信号の対応する周波数ビンの振幅スペクトルとに基づいて、当該周波数ビンのエコーパス特性を推定する周波数ビンエコーパス特性推定手段と、（５）周波数ビンエコーパス特性推定手段により推定された周波数ビン毎のエコーパス特性と遠端出力信号の振幅スペクトルとに基づいて推定エコー信号を算出する推定エコー信号算出手段と、（６）近端入力信号の振幅スペクトルからサプレスゲインを求め、音響エコー信号を抑圧するエコー抑圧手段とを備えることを特徴とするエコー抑圧装置である。
【００１６】
第２の本発明は、コンピュータを、（１）入力された遠端出力信号を周波数スペクトルに変換して、遠端出力信号の振幅スペクトルを求める遠端出力信号振幅スペクトル算出手段と、（２）入力された近端入力信号を周波数領域に変換して、近端入力信号の振幅スペクトルを求める近端入力信号振幅スペクトル算出手段と、（３）遠端出力信号の振幅スペクトルと閾値とを周波数ビン毎に比較して、周波数ビン毎に周波数成分の有無を判定する周波数ビン成分有無判定手段と、（４）周波数ビン成分有無判定手段の判定結果により遠端出力信号の周波数成分があると判断された周波数ビンの振幅スペクトルと、近端入力信号の対応する周波数ビンの振幅スペクトルとに基づいて、当該周波数ビンのエコーパス特性を推定する周波数ビンエコーパス特性推定手段と、（５）周波数ビンエコーパス特性推定手段により推定された周波数ビン毎のエコーパス特性と遠端出力信号の振幅スペクトルとに基づいて推定エコー信号を算出する推定エコー信号算出手段と、（６）近端入力信号の振幅スペクトルからサプレスゲインを求め、音響エコー信号を抑圧するエコー抑圧手段として機能させることを特徴とするエコー抑圧プログラムである。
【発明の効果】
【００１７】
本発明によれば、遠端入力信号の周波数成分の小さい周波数ビンについて、過去のエコーパス特性からエコーパス特性を算出し更新することで、エコーパス特性の誤学習を防止し、エコー抑圧性能を維持することができる。
【図面の簡単な説明】
【００１８】
【図１】第１の実施形態のエコー抑圧装置の構成を示すブロック図である。
【図２】第１の実施形態のエコー抑圧装置のエコーパス特性更新部の詳細な構成を示すブロック図である。
【図３】第２の実施形態のエコー抑圧装置のエコーパス特性更新部の詳細な構成を示すブロック図である。
【図４】従来のエコー抑圧装置を説明するためのブロック図である。
【発明を実施するための形態】
【００１９】
（Ａ）第１の実施形態
以下では、本発明のエコー抑圧装置及びエコー抑圧プログラムの第１の実施形態を、図面を参照しながら詳細に説明する。
【００２０】
第１の実施形態は、例えば、テレビ会議システムや電話会議システム等の拡声通話システムの音声送受信装置のエコー抑圧装置及びプログラムに本発明を適用する場合を例示する。
【００２１】
（Ａ−１）第１の実施形態の構成
図１は、第１の実施形態に係るエコー抑圧装置１００の構成を示すブロック図である。
【００２２】
第１の実施形態に係るエコー抑圧装置１００は、上述した非特許文献１に記載されるエコー抑圧装置を実装したものであり、遠端出力信号が音響エコー信号として近端側のマイクに回り込んで入力されたときでも、近端入力信号に対してエコーサプレス処理（エコー抑圧処理）を行い、音響エコー信号を適切に抑圧するものである。
【００２３】
第１の実施形態のエコー抑圧装置１００は、例えば専用ボードとして構築されるようにしても良いし、ＤＳＰ（デジタルシグナルプロセッサ）へのエコー抑圧プログラムの書き込みによって実現されたものであっても良く、ＣＰＵと、ＣＰＵが実行するソフトウェア（エコー抑圧プログラム）によって実現されたものであっても良いが、機能的には、図１で表すことができる。
【００２４】
図１において、第１の実施形態に係るエコー抑圧装置１００は、遠端入力信号入力部１０１、デジタル−アナログ（ＤＡ）変換器１０２、スピーカ１０３、マイクロホン（以下、マイクと呼ぶ。）１０４、アナログ−デジタル（ＡＤ）変換器１０５、遠端信号周波数領域変換部１０６、遠端出力信号振幅スペクトル計算部１０７、エコーパス特性保持部１０８、推定エコー信号計算部１０９、近端入力信号周波数領域変換部１１０、近端入力信号振幅スペクトル計算部１１１、エコーサプレスゲイン計算部１１２、エコーサプレス部１１３、近端出力信号振幅スペクトル計算部１１４、エコーパス特性更新部１１５、近端出力信号時間領域変換部１１６、近端出力信号出力部１１７を有する。
【００２５】
遠端入力信号入力部１０１は、遠端側からデジタル音信号が遠端入力信号として入力され、その入力された遠端入力信号を遠端出力信号としてＤＡ変換器１０２及び遠端出力信号周波数領域変換部１０６に与えるものである。遠端入力信号は、エコー抑圧装置１００の近端に対して遠端に位置する装置から有線回線又は無線回線を通じて入力されたデジタル音信号である。
【００２６】
ＤＡ変換器１０２は、入力された遠端出力信号に基づく音をスピーカ１０３から出力するために、遠端出力信号（デジタル信号）をアナログ信号に変換してスピーカ１０３に出力するものである。
【００２７】
スピーカ１０３は、ＤＡ変換器１０２から出力されたアナログ音信号を音信号として出力するものである。
【００２８】
マイク１０４は、近端話者が発した音声や環境音や音響エコー信号（スピーカ１０３から出力された音信号が近端側の空間を伝達してマイク１０４に入力される音）等の音信号が重畳して入力され、アナログ音信号（電気信号）に変換され、その変換したアナログ音信号をＡＤ変換器１０５に出力するものである。
【００２９】
ＡＤ変換器１０５は、マイク１０４からのアナログ音信号をデジタル音信号に変換して、その変換後のデジタル音信号を近端入力信号として近端入力信号周波数領域変換部１１０に出力するものである。
【００３０】
遠端出力信号周波数領域変換部１０６は、遠端入力信号入力部１０１からの遠端出力信号を、例えば、高速フーリエ変換等により、周波数スペクトルに変換し遠端出力信号の周波数スペクトルを求め、変換した遠端出力信号の周波数スペクトルを遠端出力力信号振幅スペクトル計算部１０７に出力する。
【００３１】
遠端出力信号振幅スペクトル計算部１０７は、遠端出力信号周波数領域変換部１０６から出力された遠端出力信号の周波数スペクトルに基づいて、遠端出力信号の振幅スペクトルを算出するものである。遠端出力信号振幅スペクトル計算部１０７は、算出した遠端出力信号の振幅スペクトルをエコーパス特性更新部１１５及び推定エコー信号計算部１０９に出力する。
【００３２】
エコーパス特性保持部１０８は、推定したエコーパス特性を周波数ビン毎に保持するものである。
【００３３】
推定エコー信号計算部１０９は、遠端出力信号振幅スペクトル計算部１０７から出力された遠端出力信号の振幅スペクトルと、エコーパス特性保持部１０８に保持されている推定したエコーパス特性とを乗じて推定エコー信号を算出するものである。推定エコー信号計算部１０９は、算出した推定エコー信号をエコーサプレスゲイン計算部１１２に出力する。
【００３４】
近端入力信号周波数領域変換部１１０は、ＡＤ変換器１０５から出力される近端入力信号を、例えば、高速フーリエ変換等により、周波数スペクトルに変換し近端入力信号の周波数スペクトルを求め、変換した近端入力信号の周波数スペクトルを近端入力信号振幅スペクトル計算部１１１及びエコーサプレス部１１３に出力する。
【００３５】
近端入力信号振幅スペクトル計算部１１１は、近端入力信号周波数領域変換部１１０から出力された近端入力信号の周波数スペクトルに基づいて、近端入力信号の振幅スペクトルを算出するものである。近端入力信号振幅スペクトル計算部１１１は、算出した近端入力信号の振幅スペクトルをエコーパス特性更新部１１５及びエコーサプレスゲイン計算部１１２に出力する。
【００３６】
エコーサプレスゲイン計算部１１２は、近端入力信号からエコーを抑圧するゲインを算出するために、推定エコー信号計算部１０９により算出された推定エコー信号の振幅スペクトルと、近端入力信号振幅スペクトル計算部１１１により算出された近端入力信号の振幅スペクトルとを用いてエコーサプレスゲインを算出するものである。エコーサプレスゲイン計算部１１２は、算出したエコーサプレスゲインをエコーサプレス部１１３に出力する。
【００３７】
エコーサプレス部１１３は、エコーサプレスゲイン計算部１１２により算出されたエコーサプレスゲインと、近端入力信号周波数領域変換部１１０から出力された近端入力信号の周波数スペクトルを乗じることにより、エコーが抑圧された近端入力信号の周波数スペクトルを求めるものである。エコーサプレス部１１３は、エコーを抑圧した近端入力信号の周波数スペクトルを近端出力信号として、近端出力信号振幅スペクトル計算部１１４及び近端出力信号時間領域変換部１１６に出力する。
【００３８】
近端出力信号振幅スペクトル計算部１１４は、エコーサプレス部１１３から出力されたエコー抑圧後の近端出力信号の周波数スペクトルに基づいて、近端出力信号の振幅スペクトルを算出するものである。近端出力信号振幅スペクトル計算部１１４は、算出した近端出力信号の振幅スペクトルをエコーパス特性更新部１１５に出力する。
【００３９】
エコーパス特性更新部１１５は、遠端出力信号振幅スペクトル計算部１０７からの遠端出入力信号の振幅スペクトルと、近端入力信号振幅スペクトル計算部１１１からの近端入力信号の振幅スペクトルと、近端出力信号振幅スペクトル計算部１１４からの近端出力信号の振幅スペクトルとに基づいてシングルトーク（近端入力信号に音響エコー信号が存在し、近端話者の音声が存在していない状態）かダブルトーク（近端入力信号に音響エコー信号が存在し、近端話者の音声が存在している状態）を判定する。
【００４０】
シングルトークと判定されたときに、そのシングルトークと判定された処理フレームの遠端出力信号の振幅スペクトルと閾値とを周波数ビン毎に比較する。そして、周波数ビン毎の遠端出力信号の振幅スペクトルが閾値より大きい場合、エコーパス特性更新部１１５は、その周波数ビンの遠端出力信号の振幅スペクトルと近端入力信号の振幅スペクトルとに基づいて、当該周波数ビンの現時刻のエコーパス特性を算出する。そして、エコーパス特性更新部１１５は、エコーパス特性保持部１０８に保持されている推定エコーパス特性と算出した現時刻のエコーパス特性とを用いて推定エコーパス特性を更新し、エコーパス特性保持部１０８に保持させる。なお、エコーパス特性更新部１１５は、エコーパス特性推定手段に含まれるものである。
【００４１】
つまり、エコーパス特性更新部１１５は、シングルトークのときに、エコーパス特性を周波数ビン毎に更新するものであって、閾値よりも大きい遠端出力信号の周波数ビンの振幅スペクトルについては現時刻のエコーパス特性を算出し、その算出した現時刻のエコーパス特性に更新する。閾値以下の遠端出力信号の周波数ビンの振幅スペクトルについてはエコーパス特性の更新を行わず、過去のエコーパス特性を用いてエコーパス特性を求める。すなわち、遠端出力信号の周波数成分が閾値以下となる周波数ビンについては、その周波数ビンのエコーパス特性の更新しないため、エコーパス特性の学習を防止できる。
【００４２】
近端出力信号時間領域変換部１１６は、エコーサプレス部１１３から出力された近端出力信号の周波数スペクトルを、例えば、逆高速フーリエ変換等により、時間領域の近端出力信号に変換するものである。近端出力信号時間領域変換部１１６は変換した近端出力信号を近端出力信号出力部１１７に出力する。
【００４３】
近端出力信号出力部１１７は、近端出力信号時間領域変換部１１６から出力される近端出力信号（デジタル信号）を、通信相手である遠端側に向けて出力するものである。
【００４４】
図２は、第１の実施形態に係るエコーパス特性更新部１１５の詳細な構成を示す構成図である。
【００４５】
図２において、第１の実施形態に係るエコーパス特性更新部１１５は、シングルトーク判定部２０１、周波数ビン信号有無判定部２０２、周波数ビンエコーパス特性計算部２０３、周波数ビンエコーパス特性更新部２０４を有する。
【００４６】
シングルトーク判定部２０１は、シングルトークとダブルトークを判断するものである。ここで、シングルトークとは、遠端側のみ音信号が発生し、その音がスピーカ１０３から出力され近端側の空間を伝達し音響エコー信号として回り込んでマイク１０４に受音されている状態である。またダブルトークは、遠端側と近端側との両方が音声信号を発生している状態である。
【００４７】
シングルトーク判定部２０１は、シングルトークとダブルトークを判断するために、例えば、近端入力信号振幅スペクトル計算部１１１からの近端入力信号の振幅スペクトルと、近端出力信号スペクトル計算部１１４からの近端出力信号の振幅スペクトルとを処理フレーム毎に比較し、近端入力信号の振幅スペクトルに対して近端出力信号の振幅スペクトルが非常に小さい場合はシングルトーク、それ以外の場合はダブルトーク、又は近端側の音声のみ入力されている状態と判定する。そして、シングルトーク判定部２０１は、シングルトークであると判定すると、シングルトークであると判定した処理フレームを周波数ビンエコーパス特性計算部２０３に与える。
【００４８】
なお、シングルトーク判定部２０１によるシングルトークの判定手法は、種々の方法を広く適用することができ、例えば、シングルトーク判定部２０１が、遠端出力信号の振幅スペクトルと、近端入力信号の振幅スペクトルとを入力し、相関が強い場合に、シングルトークと判定するとしても良い。
【００４９】
周波数ビン信号有無判定部２０２は、遠端出力信号振幅スペクトル計算部１０７からの遠端出力信号の振幅スペクトルに基づいて、周波数ビン毎に信周波数成分の有無を判定するものである。つまり、周波数ビン信号有無判定部２０２は、遠端出力信号の振幅スペクトルと閾値とを周波数ビン毎に比較する。そして、遠端出力信号の振幅スペクトルが閾値より大きい場合、当該周波数ビンの周波数成分があると判定する。一方、遠端出力信号の振幅スペクトルが閾値以下の場合、当該周波数ビンの周波数成分がないと判定する。周波数ビン信号有無判定部２０２は、周波数ビン毎に信号有無判定結果を周波数ビンエコーパス特性計算部２０３に与える。
【００５０】
周波数ビンエコーパス特性計算部２０３は、シングルトークの場合に、エコーパス特性を周波数ビン毎にエコーパス特性を算出するものである。周波数ビンエコーパス特性計算部２０３は、シングルトーク判定部２０１からシングルトークであると判定した処理フレームの、周波数ビン信号有無判定部２０２からの周波数ビン信号有無判定で周波数成分があると判定された周波数ビンに対して、当該処理フレームの当該周波数ビンのエコーパス特性を、遠端出力信号振幅スペクトル計算部１０７からの遠端出力信号の振幅スペクトルと、近端出力信号振幅スペクトル計算部１１４からの近端出力信号の振幅スペクトルとを用いて算出するものである。周波数ビンエコーパス特性計算部２０３は、算出した周波数ビンのエコーパス特性を周波数ビンエコーパス特性更新部２０４に与える。
【００５１】
周波数ビンパス特性更新部２０４は、周波数ビンエコーパス特性計算部２０３により算出された周波数ビン毎のエコーパス特性を取得すると、エコーパス特性保持部１０８に保持されるエコーパス特性と、周波数ビンエコーパス特性計算部２０３から取得した周波数ビンのエコーパス特性を用いて更新するものである。このとき、周波数ビン信号有無判定部２０２で周波数成分があると判断された周波数ビンのみ更新し、周波数成分がないと判断された周波数ビンは更新しない。これにより、シングルトークのときに、遠端出力信号の振幅スペクトルが閾値以上と判断された周波数ビン毎のエコーパス特性を更新することができる。
【００５２】
（Ａ−２）第１の実施形態の動作
次に、第１の実施形態に係るエコー抑圧装置１００におけるエコー抑圧処理の動作を詳細に説明する。
【００５３】
まず、エコー抑圧装置１００の動作開始後、例えばインターネットプロトコル（ＩＰ）網等のネットワークや、携帯電話等の無線ネットワークの電波を介して、遠端側からの音信号（デジタル音信号）が、遠端入力信号として遠端入力信号入力部１０１に入力される。
【００５４】
遠端入力信号は、遠端出力信号としてＤＡ変換器１０２及び遠端出力信号周波数領域変換部１０６に与えられる。ＤＡ変換器１０２では、遠端入力信号入力部１０１から入力されたデジタル信号である遠端出力信号がアナログ信号に変換されてスピーカ１０３に出力され、スピーカ１０３から遠端側の音信号が出力される。
【００５５】
また、マイク１０４は近端話者が発した音声や環境音、音響エコー信号（スピーカ１０３から出力された音信号が近端側の空間を伝達してマイク１０４に入力される音）等の音信号が重畳して受音され、その受音された音はマイク１０４によりアナログ音信号（電気信号）に変換されて、アナログ音信号が近端入力信号としてＡＤ変換器１０５に与えられる。
【００５６】
マイク１０４から出力されたアナログ音信号はＡＤ変換器１０５に与えられる。ＡＤ変換器１０５では、アナログ音信号が、標本化、量子化、符号化の処理がなされてデジタル音信号に変換され、近端入力信号として近端入力信号周波数領域変換部１１０に与えられる。
【００５７】
遠端出力信号周波数領域変換部１０６では、例えば、高速フーリエ変換等により、遠端出力信号が周波数スペクトルＲＯＵＴ（ｉ，ω）に変換され、その変換された遠端出力信号の周波数スペクトルＲＯＵＴ（ｉ，ω）は遠端出力信号振幅スペクトル計算部１０７に与えられる。
【００５８】
遠端出力信号振幅スペクトル計算部１０７では、遠端出力信号の周波数スペクトルＲＯＵＴ（ｉ，ω）を用いて、（１）式に従い、遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜が求められる。
【数１】
2015070290.tif 000003
【００５９】
ここで、ｉはフレーム、ωは周波数ビン、ＲＯＵＴ＿ｒｅａｌ（ｉ，ω）とＲＯＵＴ＿ｉｍａｇｅ（ｉ，ω）は、フレームｉにおける周波数ビンωの遠端出力信号の周波数スペクトルの実数部と虚数部を示しており、遠端出力信号の周波数スペクトルＲＯＵＴ（ｉ，ω）は（２）で表すことができる。（２）のｊは虚数を表している。
【数２】
2015070290.tif 000004
【００６０】
遠端出力信号振幅スペクトル計算部１０７により求められた遠端入力信号の周波数スペクトルは、エコーパス特性更新部１１５及び推定エコー信号計算部１０９に与えられる。
【００６１】
推定エコー信号計算部１０９では、エコーパス保持部１０８に保持されているエコーパス特性｜Ｈ（ｉ，ω）｜と、遠端出力信号振幅スペクトル計算部１０７からの遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜とを用いて、（３）式により、推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜が求められる。つまり、推定エコー信号計算部１０９は、遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜に、エコーパス保持部１０８に保持されているエコーパス特性｜Ｈ（ｉ，ω）｜の対応する周波数ビンを乗じて、当該周波数ビンの推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜を求める。
【数３】
2015070290.tif 000005
【００６２】
一方、近端入力信号周波数領域変換部１１０では、例えば、高速フーリエ変換等により、近端入力信号が周波数スペクトルに変換され、その変換された近端入力信号の周波数スペクトルは、近端入力信号振幅スペクトル計算部１１１及びエコーサプレス部１１３に与えられる。
【００６３】
近端入力信号振幅スペクトル計算部１１１は、近端入力信号周波数領域変換部１１０からの近端入力信号の周波数スペクトルＳＩＮ（ｉ，ω）を用いて、（４）式に従い、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜が求められる。
【数４】
2015070290.tif 000006
【００６４】
ここで、ＳＩＮ＿ｒｅａｌ（ｉ，ω）とＳＩＮ＿ｉｍａｇｅ（ｉ，ω）は、フレームｉにおける周波数ビンωの近端入力信号の周波数スペクトルの実数部と虚数部を示しており、近端入力信号の周波数スペクトルＳＩＮ（ｉ，ω）は（５）式で表すことができる。（５）のｊは虚数を表している。
【数５】
2015070290.tif 000007
【００６５】
そして、近端入力信号振幅スペクトル計算部１１１により求められた近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜は、エコーサプレスゲイン計算部１１２及びエコーパス特性更新部１１５に与えられる。
【００６６】
エコーサプレスゲイン計算部１１２では、推定エコー信号計算部１０９からの推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜と、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜とを用いて、（６）式に従って、エコーサプレスゲインＧ（ｉ，ω）が求められる。つまり、エコーサプレスゲイン計算部１１２は、周波数ビン毎に、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜から推定エコー信号の振幅スペクトル｜ＥＣＨＯ（ｉ，ω）｜を差し引いた振幅スペクトルを、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜で割ることで、エコーサプレスゲインＧ（ｉ，ω）を求める。
【数６】
2015070290.tif 000008
【００６７】
エコーサプレス部１１３では、近端入力信号周波数領域変換部１１０からの近端入力信号のスペクトルＳＩＮ（ｉ，ω）の実数部と虚数部と、エコーサプレスゲイン計算部１１２からのエコーサプレスゲインＧ（ｉ，ω）とを用いて、（７）式、（８）式に従い、近端入力信号のスペクトルＳＩＮ（ｉ，ω）に重畳されている音響エコー信号を抑圧する。つまり、エコーサプレス部１１３は、周波数スペクトルの実数部、虚数部にエコーサプレスゲインＧ（ｉ，ω）を周波数ビン毎に乗じて、音響エコー信号を抑圧した近端出力信号の周波数スペクトルを求める。
【数７】
2015070290.tif 000009
【００６８】
ここで、ＳＯＵＴ＿ｒｅａｌ（ｉ，ω）とＳＯＵＴ＿ｉｍａｇｅ（ｉ，ω）は、時刻ｉにおける周波数ビンωの近端出力信号の周波数スペクトルの実数部と虚数部を示しており、近端出力信号の周波数スペクトルＳＯＵＴ（ｉ，ω）は（９）式で表すことができる。（９）のｊは虚数を表している。
【数８】
2015070290.tif 000010
【００６９】
そして、音響エコー信号が抑圧された近端出力信号の周波数スペクトルは、近端出力信号振幅スペクトル計算部１１４に与えられる。近端出力信号振幅スペクトル計算部１１４では、近端出力信号の周波数スペクトルを用いて、（１０）式に従って、近端出力信号の振幅スペクトル｜ＳＯＵＴ(ｉ，ω)｜が求められる。
【数９】
2015070290.tif 000011
【００７０】
エコーパス特性更新部１１５では、シングルトーク判定部２０１が、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜と、近端出力信号の振幅スペクトル｜ＳＯＵＴ（ｉ，ω）｜とを処理フレーム毎に比較して、シングルトークを判定する。例えば、シングルトーク判定部２０１は、近端入力信号振幅スペクトル計算部１１１からの近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜と、近端出力信号スペクトル計算部１１４からの近端出力信号の振幅スペクトル｜ＳＯＵＴ（ｉ，ω）｜とを処理フレーム毎に比較し、近端入力信号の振幅スペクトルに対して近端出力信号の振幅スペクトルが非常に小さい場合はシングルトーク、それ以外の場合はダブルトーク、又は近端側の音声のみ入力されている状態と判定する。
【００７１】
また、周波数ビン信号有無判定部２０２は、シングルトークと判断されたときに、そのシングルトークと判断された処理フレームの遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜と閾値ＴＨ（ｉ，ω）を周波数ビン毎に比較し周波数ビン毎に周波数成分の有無を判定する。遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜が閾値ＴＨ（ｉ，ω）より大きい場合、当該周波数ビンの周波数成分があると判定し、遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜が閾値ＴＨ（ｉ，ω）以下の場合、当該周波数ビンの周波数成分がないと判定し、その判定結果を周波数ビンエコーパス特性計算部２０３に与える。
【００７２】
閾値ＴＨ（ｉ，ω）は、遠端出力信号の背景雑音を周波数成分がありと判別するのを回避するために、遠端出力信号の背景雑音より大きい値が望ましい(例えば遠端出力信号の背景雑音のフレーム平均値の２倍など)。また、閾値ＴＨ（ｉ，ω）は、周波数ビン毎に異なる値を設定するようにしても良い。
【００７３】
周波数ビンエコーパス特性計算部２０３では、シングルトーク判定部２０１でシングルトークと判断された処理フレームの、周波数ビン信号有無判定部２０２で遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜の周波数成分があると判定された周波数ビンに対して、（１１）式に従い、当該処理フレームの当該周波数ビンのエコーパス特性｜Ｈ１（ｉ，ω）｜を遠端出力信号振幅スペクトル計算部１０７からの遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜と、近端出力信号振幅スペクトル計算部１１４からの近端出力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜を用いて求める。
【数１０】
2015070290.tif 000012
【００７４】
周波数ビンエコーパス特性更新部２０４は、エコーパス特性保持部１０８に保持されているエコーパス特性｜Ｈ（ｉ，ω）｜を読み出し、エコーパス特性保持部１０８に保持されるエコーパス特性｜Ｈ（ｉ，ω）｜と、周波数ビンエコーパス特性計算部２０３から取得したエコーパス特性｜Ｈ１（ｉ，ω）｜を用いて更新する。このとき、周波数ビン信号有無判定部２０２で周波数成分があると判断された周波数ビンのみエコーパス特性を更新し、周波数成分がないと判断された周波数ビンはエコーパス特性の更新を行わず、エコーパス特性保持部１０８に保持されている周波数ビンのエコーパス特性｜Ｈ（ｉ，ω）｜を用いる。
【００７５】
具体的には、周波数ビンエコーパス特性構成部２０４は、（１２）式に従って、エコーパス特性｜Ｈ（ｉ，ω）｜を更新する。
【数１１】
2015070290.tif 000013
【００７６】
ａは時定数フィルタの係数であり、ａは１以下の値であって、１に近い値が望ましく、例えばａ＝０．９９等の値とすることができる。
【００７７】
つまり、遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜の周波数成分がある周波数ビンについて、周波数ビンエコーパス特性構成部２０４は、エコーパス特性保持部１０８に保持されている当該周波数ビンのエコーパス特性｜Ｈ（ｉ，ω）｜と、周波数ビンエコーパス特性計算部２０３により求められた当該周波数ビンのエコーパス特性｜Ｈ１（ｉ，ω）｜とを用いて、次のフレームで用いる周波数ビンのエコーパス特性｜Ｈ（ｉ＋１，ω）｜を求める。また、遠端出力信号の振幅スペクトルＲＯＵＴ（ｉ，ω）の周波数成分がない周波数ビンについて、周波数ビンエコーパス特性構成部２０４は、エコーパス特性保持部１０８に保持されている当該周波数ビンのエコーパス特性｜Ｈ（ｉ，ω）｜を用いる。
【００７８】
エコーサプレス部１１３により音響エコー信号が抑圧された近端出力信号の周波数スペクトルＳＯＵＴ（ｉ，ω）を（９）に従い、ＳＯＵＴ＿ｒｅａｌ（ｉ，ω）とＳＯＵＴ＿ｉｍａｇｅ（ｉ，ω）から求めて、近端出力信号の周波数スペクトルＳＯＵＴ（ｉ，ω）を近端出力信号時間領域変換部１１６に与えられる。
【００７９】
近端出力信号時間領域変換部１１６では、近端出力信号のスペクトルＳＯＵＴ（ｉ，ω）が、例えば、逆高速フーリエ変換等により、時間領域の信号に変換され、変換された時間領域の近端出力信号が近端出力信号出力部１１７に与えられる。近端出力信号出力部１１７は、近端信号時間領域変換部１１６からの近端出力信号をデジタル信号でネットワーク網（ＩＰ網）や携帯電話等の無線ネットワークで使用される電波を介して遠端側に出力する。
【００８０】
（Ａ−３）第１の実施形態の効果
以上のように、第１の実施形態によれば、遠端出力信号の周波数ビンの周波数成分が小さいときでもエコーパス特性を更新しないことでエコーパス特性の誤学習を防止し、エコーを抑圧することができる。
【００８１】
（Ｂ）第２の実施形態
次に、本発明のエコー抑圧装置及びエコー抑圧プログラムの第２の実施形態を、図面を参照しながら説明する。
【００８２】
（Ｂ−１）第２の実施形態の構成
第２の実施形態では、エコーパス特性更新部の構成が、第１の実施形態に係るエコー抑圧装置１００と異なるのみであり、それ以外の構成要素は第１の実施形態に係る図１のエコー抑圧装置１００の構成要素と同一又は対応するものである。そこで、以下では、第２の実施形態に係るエコーパス特性更新部３０１の構成を図１及び図３を参照しながら詳細に説明する。
【００８３】
図３は、第２の実施形態に係るエコーパス特性更新部３０１の詳細な構成を示す構成図である。
【００８４】
図３において、第２の実施形態に係るエコーパス特性更新部３０１は、第１の実施形態で説明した、シングルトーク判定部２０１、周波数ビン信号有無判定部２０２、周波数ビンエコーパス特性計算部２０３、周波数ビンエコーパス特性更新部２０４に加えて、エコーパス特性補間部３０２を有する。
【００８５】
なお、図３において、第１の実施形態に係るエコーパス特性更新部１１５の構成要素と同一又は対応するものについては同一の符号を付している。また、第１の実施形態と同一又は対応する構成要素の詳細な説明は重複することになるためここでは省略する。
【００８６】
エコーパス特性補間部３０２は、周波数ビンエコーパス特性計算部２０３により求められた周波数ビン毎のエコーパス特性を補間して、補間後の周波数ビンのエコーパス特性を周波数ビンエコーパス特性更新部２０４に与えるものである。
【００８７】
これにより、周波数ビンエコーパス特性更新部２０４は、エコーパス特性補間部３０２により補正された周波数ビンのエコーパス特性を用いて、エコーパス特性保持部１０８に保持されている周波数ビン毎のエコーパス特性を更新する。
【００８８】
（Ｂ−２）第２の実施形態の動作
次に、第２の実施形態のエコー抑圧装置１００におけるエコー抑圧処理の動作を詳細に説明する。ここでは、第２の実施形態に係るエコー抑圧装置１００は、エコーパス特性更新部３０１における処理動作を詳細に説明する。
【００８９】
エコーパス特性更新部３０１では、シングルトーク判定部２０１が、第１の実施形態と同様に、近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ、ω）｜と近端出力信号近端出力信号の振幅スペクトル｜ＳＯＵＴ（ｉ、ω）｜とを用いて処理フレーム毎にシングルトークか否かの判定を行う。シングルトークと判断されたとき、周波数ビン信号有無判定部２０２が遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ、ω）｜から周波数ビン毎の周波数成分の有無を判定する。そして、周波数ビンエコーパス特性計算部２０３は、第１の実施形態と同様に、シングルトークと判断された処理フレームで、遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ、ω）｜の周波数成分がある周波数ビンについて、遠端出力信号の振幅スペクトル｜ＲＯＵＴ（ｉ，ω）｜と近端入力信号の振幅スペクトル｜ＳＩＮ（ｉ，ω）｜とを用いて、（１１）式に従って、周波数ビン毎のエコーパス特性｜Ｈ１（ｉ、ω１）｜を求める。
【００９０】
周波数ビンエコーパス特性計算部２０３により求められた周波数ビン毎のエコーパス特性｜Ｈ１（ｉ、ω１）｜は、エコーパス特性補間部３０２に与えられる。
【００９１】
エコーパス特性補間部３０２は、周波数ビンエコーパス特性計算部２０３から周波数ビンのエコーパス特性｜Ｈ１（ｉ、ω１）｜を取得する。このとき、周波数ビン信号有無判定部２０２により遠端出力信号の振幅スペクトルで周波数成分があると判断された周波数ビン（ここでは第１の周波数ビンω１と呼ぶ。）と、その周波数ビンの次に周波数成分があると判断された周波数ビン（ここでは第２の周波数ビンω２と呼ぶ。）との間で、周波数成分がないと判断された周波数ビンについて、第１の周波数ビンのエコーパス特性｜Ｈ１（ｉ、ω１）｜と、第２の周波数ビンのエコーパス特性｜Ｈ１（ｉ、ω２）｜とに基づいて、当該信周波数成分がないと判断した周波数ビンのエコーパス特性｜Ｈ１（ｉ、ω）｜を補間する。
【００９２】
ここで、エコーパス特性補間部３０２による補間方法は、種々の方法を適用できるが、例えば、第１の周波数ビンのエコーパス特性｜Ｈ１（ｉ、ω１）｜分と、第２の周波数ビンのエコーパス特性｜Ｈ１（ｉ、ω２）｜とを用いて、（１３）式に従って、線形補間により周波数成分がないと判断された周波数ビンのエコーパス特性｜Ｈ１（ｉ、ω）｜を求めるようにしても良い。
【数１２】
2015070290.tif 000014
【００９３】
そして、周波数ビンエコーパス特性更新部２０４では、第１の実施形態と同様に、エコーパス保持部１０８に保持されているエコーパス特性｜Ｈ（ｉ，ω）｜とエコーパス特性補間部３０２で補間されたエコーパス特性｜Ｈ１（ｉ，ω）｜を用いて、（１４）式に従って、周波数ビン毎のエコーパス特性を更新する。
【数１３】
2015070290.tif 000015
【００９４】
（Ｂ−３）第２の実施形態の効果
以上のように、第２の実施形態によれば、第１の実施形態の効果に加えて、遠端出力信号の振幅スペクトルが小さい周波数ビンでも、算出したエコーパス特性を補間することができ、エコーパス特性をその補間した周波数成分が小さい周波数ビンのエコーパス特性を更新することにより、近端入力信号の周波数成分が小さい周波数ビンのエコーパス特性の更新し、音響エコー信号を抑圧することができる。
【００９５】
（Ｃ）他の実施形態
上述した各実施形態においても、種々の変形実施形態を説明したが、本発明は以下の変形実施形態についても適用することができる。
【００９６】
上述した各実施形態で説明したエコー抑圧装置は、例えば、テレビ会議システムや電話会議システム等に用いられる音声通信装置を含む装置に搭載されるようにしても良い。また、携帯電話機やスマートホン等の携帯端末に本発明のエコー抑圧装置は搭載されるようにしても良い。
【００９７】
上述した各実施形態では、シングルトークの場合に周波数ビン毎のエコーパス特性の更新を行うことを説明した。シングルトーク判定部がシングルトークではないと判定した場合には、周波数ビン毎のエコーパス特性の更新を行わないようにする。
【符号の説明】
【００９８】
１００…エコー抑圧装置、１０１…遠端入力信号入力部、１０２…ＤＡ変換器、１０３…スピーカ、１０４…マイク、１０５…ＡＤ変換器、１０６…遠端出力信号周波数領域変換算部、１０７…遠端出力信号振幅スペクトル計算部、１０８…エコーパス特性保持部、１０９…推定エコー信号計算部、１１０…近端入力信号周波数領域変換部、１１１…近端入力信号振幅スペクトル計算部、１１２…エコーサプレスゲイン計算部、１１３…エコーサプレス部、１１４…近端出力信号振幅スペクトル計算部、１１５、３０１…エコーパス特性更新部、１１６…近端出力信号時間領域変換部、１１７…近端出力信号出力部、
２０１…シングルトーク判定部、２０２…周波数ビン信号有無判定部、２０３…周波数ビンエコーパス特性計算部、２０４…周波数ビンエコーパス特性更新部、３０２…エコーパス特性補間部。
【図１】
2015070290.tif 000016
【図２】
2015070290.tif 000017
【図３】
2015070290.tif 000018
【図４】
2015070290.tif 000019