(57)【特許請求の範囲】
【請求項１】
第１の音声信号と、音声出力部により再生された該第１の音声信号を音声入力部により集音した音を表す第２の音声信号とから、該音声出力部から該音声入力部への音の経路のインパルス応答を求め、該インパルス応答に応じて残響特性を算出する残響特性算出部と、
音源からの音を前記音声入力部により集音した音を表す第３の音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて前記音声入力部から前記音源までの距離を推定する距離推定部と、
前記距離の推定値が長くなるほど、前記残響特性を大きくするよう補正する残響特性補正部と、
前記補正された残響特性に応じて前記第３の音声信号に含まれる残響成分の周波数スペクトルを推定し、周波数ごとに、前記第３の音声信号の周波数スペクトルと該残響成分の周波数スペクトルの差を算出し、当該差が小さい周波数ほど前記第３の音声信号を減衰させる残響抑制部と、
を有する残響抑制装置。
【請求項２】
前記距離推定部は、前記第３の音声信号の周波数特性として、前記第３の音声信号が持つ周波数帯域全体のパワースペクトルに対する所定の周波数以下の周波数帯域のパワースペクトルの比を算出し、当該比が大きいほど前記距離の推定値を短くする、請求項１に記載の残響抑制装置。
【請求項３】
前記距離推定部は、前記音声入力部を有する機器について接話設定されている第１の期間における前記第３の音声信号の音量の平均値を算出し、前記機器が接話設定されていない第２の期間における前記第３の音声信号の音量を前記平均値から減じた値が大きいほど、当該第２の期間における前記距離の推定値を長くする、請求項１に記載の残響抑制装置。
【請求項４】
前記残響特性補正部は、第１の周波数における前記残響特性を除する係数を、前記第１の周波数よりも高い第２の周波数における前記残響特性を除する係数よりも大きくする、請求項１〜３の何れか一項に記載の残響抑制装置。
【請求項５】
第１の音声信号と、音声出力部により再生された該第１の音声信号を音声入力部により集音した音を表す第２の音声信号とから、該音声出力部から該音声入力部への音の経路のインパルス応答を求め、該インパルス応答に応じて残響特性を算出し、
音源からの音を前記音声入力部により集音した音を表す第３の音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて前記音声入力部から前記音源までの距離を推定し、
前記距離の推定値が長くなるほど、前記残響特性を大きくするよう補正し、
前記補正された残響特性に応じて前記第３の音声信号に含まれる残響成分の周波数スペクトルを推定し、
周波数ごとに、前記第３の音声信号の周波数スペクトルと該残響成分の周波数スペクトルの差を算出し、当該差が小さい周波数ほど前記第３の音声信号を減衰させる、
ことを含む残響抑制方法。
【請求項６】
第１の音声信号と、音声出力部により再生された該第１の音声信号を音声入力部により集音した音を表す第２の音声信号とから、該音声出力部から該音声入力部への音の経路のインパルス応答を求め、該インパルス応答に応じて残響特性を算出し、
音源からの音を前記音声入力部により集音した音を表す第３の音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて前記音声入力部から前記音源までの距離を推定し、
前記距離の推定値が長くなるほど、前記残響特性を大きくするよう補正し、
前記補正された残響特性に応じて前記第３の音声信号に含まれる残響成分の周波数スペクトルを推定し、
周波数ごとに、前記第３の音声信号の周波数スペクトルと該残響成分の周波数スペクトルの差を算出し、当該差が小さい周波数ほど前記第３の音声信号を減衰させる、
ことをコンピュータに実行させる残響抑制用コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば、マイクロホンを用いて集音した音声信号に含まれる残響成分を抑制する残響抑制装置、残響抑制方法及び残響抑制用コンピュータプログラムに関する。
【背景技術】
【０００２】
音を反射しやすい物体に囲まれた空間に配置されたマイクロホンでそのマイクロホンの周囲の音を集音する場合、音源から直接マイクロホンに到達する音と、音源から発した音がその物体で１回以上反射されてから間接的にマイクロホンに到達する音がある。この間接的にマイクロホンに到達する音は残響音と呼ばれる。残響音がマイクロホンに達するタイミングは、音源からの音が直接マイクロホンに達するタイミングよりも遅いため、残響音が存在すると、マイクロホンで集音された音と残響音とが重畳されて、その音は聞き取り難くなる。特に近年、防水機能を備えた携帯電話機が利用可能となっている。このような携帯電話機が浴室といった残響が生じ易い空間で使用されると、その携帯電話機により集音された音は残響音のために非常に聞き取り難くなるおそれがある。そこで、残響音を抑制する技術が研究されている（例えば、特許文献１を参照）。
【０００３】
例えば、特許文献１に開示された残響除去方法は、スピーカとマイクロホンとの音響結合によって形成される帰還経路のインパルス応答を適応フィルタにより適応的に同定してマイクロホンで集音した残響音声信号から帰還経路のエコー成分を推定する。そしてこの残響除去方法は、推定されたエコー成分を帰還経路の出力信号より減算し、その減算結果に含まれたエコー成分推定値の推定誤差が最小となるように適応フィルタのフィルタ係数を更新する。さらにこの残響除去方法は、エコー成分推定値の推定誤差が最小となったときのフィルタ係数を残響空間のインパルス応答に代用し、そのフィルタ係数から求めた残響空間の伝達関数とマイクロホンで集音した残響音声信号との演算から元の音声信号を求める。
【０００４】
特許文献１に開示された技術は、上記のように、スピーカから発してマイクロホンにより集音される音の残響特性によって、音源からマイクロホンに達する音の残響特性を近似している。しかし、携帯電話機のハンズフリー機能を用いてユーザが通話する場合のように、音源からマイクロホンまでの距離がスピーカからマイクロホンまでの距離に対して非常に大きくなったり、逆に小さくなることがある。このような場合、スピーカから発してマイクロホンにより集音される音の残響特性は、音源であるユーザの口からマイクロホンに達する音の残響特性と大きく異なり、その結果として、音源から発した音の残響成分が適切に除去されないことがあった。
【０００５】
一方、マイクロホンアレーを用いて音源との間の距離を推定する技術も提案されている（例えば、特許文献２を参照）。特許文献２に開示された音源距離測定装置は、複数のマイクロホンで受音された受音信号をそれぞれ周波数領域の信号に変換し、その周波数領域の信号をベクトル化して空間相関行列を算出する。そしてこの音源距離測定装置は、空間相関行列とマイクロホンの配置情報とから直接音のパワーと残響音のパワーを求め、距離−直間比データベースを参照して、直接音パワーを残響音パワーで除した直間比と対応する音源距離推定値を求める。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００６−６７１２７号公報
【特許文献２】特開２０１１−５３０６２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
特許文献２に開示された技術では、マイクロホンから音源までの距離を推定するために、複数のマイクロホンが用いられる。しかし、携帯電話機のような携帯可能な小型の機器では、筺体内のスペースが限られており、その距離を推定するために用いられる複数のマイクロホンを筺体内に配置することが困難なことがある。
【０００８】
そこで本明細書は、一つのマイクロホンを用いて、マイクロホンから音源までの距離が変動してもそのマイクロホンで集音した音に含まれる残響成分を抑制可能な残響抑制装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
一つの実施形態によれば、残響抑制装置が提供される。この残響抑制装置は、第１の音声信号と、音声出力部により再生された第１の音声信号を音声入力部により集音した音を表す第２の音声信号とから、音声出力部から音声入力部への音の経路のインパルス応答を求め、そのインパルス応答に応じて残響特性を算出する残響特性算出部と、音源からの音を音声入力部により集音した音を表す第３の音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて音声入力部から音源までの距離を推定する距離推定部と、距離の推定値が長くなるほど、残響特性を大きくするよう補正する残響特性補正部と、補正された残響特性に応じて第３の音声信号に含まれる残響成分の周波数スペクトルを推定し、第３の音声信号の周波数スペクトルと残響成分の周波数スペクトルの差が小さいほど第３の音声信号を減衰させる残響抑制部とを有する。
【００１０】
本発明の目的及び利点は、請求項において特に指摘されたエレメント及び組み合わせにより実現され、かつ達成される。
上記の一般的な記述及び下記の詳細な記述の何れも、例示的かつ説明的なものであり、請求項のように、本発明を限定するものではないことを理解されたい。
【発明の効果】
【００１１】
本明細書に開示された残響抑制装置は、一つのマイクロホンを用いて、マイクロホンから音源までの距離が変動してもそのマイクロホンで集音した音に含まれる残響成分を抑制できる。
【図面の簡単な説明】
【００１２】
【図１】第１の実施形態による残響抑制装置が実装された携帯電話機の概略構成図である。
【図２】第１の実施形態による残響抑制装置の概略構成図である。
【図３】音の経路のインパルス応答の一例を示す図である。
【図４】入力音量及び低域パワー比率と距離との関係を示す参照テーブルの一例である。
【図５】（ａ）及び（ｂ）は、それぞれ、音源からマイクロホンに達する音の経路の一例を示す図である。
【図６】（ａ）及び（ｂ）は、それぞれ、マイクロホンから音源までの推定距離と補正ゲインとの関係の一例を示す図である。
【図７】入力音のパワースペクトルと残響スペクトルの差分に対する抑制ゲインの関係の一例を示す図である。
【図８】第１の実施形態による残響抑制処理の動作フローチャートである。
【図９】接話時の平均音量とハンズフリー時の入力音量の差と、マイクロホンから音源までの距離との関係の一例を示す図である。
【図１０】第２の実施形態による残響抑制装置の距離推定部により実行される距離推定処理の動作フローチャートである。
【図１１】各実施形態またはその変形例による残響抑制装置の各部の機能を実現するコンピュータプログラムが動作することにより、残響抑制装置として動作するコンピュータの構成図である。
【発明を実施するための形態】
【００１３】
以下、図を参照しつつ、様々な実施形態による残響抑制装置について説明する。
この残響抑制装置は、スピーカから出力される再生音をマイクロホンで集音した音を用いて、音源から発してマイクロホンで集音される入力音の残響特性を近似的に求める。さらにこの残響抑制装置は、音源からマイクロホンまでの距離に応じて変化する入力音の特徴量を求めて、その特徴量により音源からマイクロホンまでの距離を推定し、その推定距離に応じて残響特性を補正する。そしてこの残響抑制装置は、補正された残響特性に従って入力音に含まれる残響成分を抑制する。
【００１４】
図１は、第１の実施形態による残響抑制装置が実装された携帯電話機の概略構成図である。図１に示されるように、携帯電話機１は、制御部２と、通信部３と、マイクロホン４と、アナログ／デジタル変換器５と、残響抑制装置６と、デジタル／アナログ変換器７と、スピーカ８とを有する。
このうち、制御部２、通信部３及び残響抑制装置６は、それぞれ別個の回路として形成される。あるいはこれらの各部は、その各部に対応する回路が集積された一つの集積回路として携帯電話機１に実装されてもよい。さらに、これらの各部は、携帯電話機１が有するプロセッサ上で実行されるコンピュータプログラムにより実現される、機能モジュールであってもよい。
【００１５】
制御部２は、少なくとも一つのプロセッサ、不揮発性のメモリ及び揮発性のメモリ及びその周辺回路を有する。制御部２は、携帯電話機１が有するキーパッドなどの操作部（図示せず）を介した操作により通話が開始されると、携帯電話機１と基地局装置（図示せず）との間における、無線接続、切断などの呼制御処理を携帯電話機１が準拠する通信規格に従って実行する。そして制御部２は、その呼制御処理の結果に応じて、通信部３に対して音声通話の開始または終了を指示する。さらに、制御部２は、通信部３を介して基地局装置から受信したダウンリンク信号に含まれる符号化された音声信号を取り出し、その音声信号を復号する。そして制御部２は、復号した音声信号を再生音声信号として、デジタル／アナログ変換器７及び残響抑制装置６へ出力する。
【００１６】
また制御部２は、残響抑制装置６がマイクロホン４を介して入力された入力音声信号に含まれる残響成分を抑制することで生成した出力音声信号を符号化し、その符号化された出力音声信号を含むアップリンク信号を生成する。そして制御部２は、そのアップリンク信号を通信部３へ渡す。なお、音声信号に対する符号化方式としては、例えば、Third Generation Partnership Project(3GPP)により標準化されたAdaptive Multi-Rate-NarrowBand(AMR-NB)方式、またはAdaptive Multi-Rate-WideBand(AMR-WB)方式などが用いられる。
【００１７】
通信部３は、基地局装置との間で無線通信する。そして通信部３は、基地局装置から無線信号を受信して、その無線信号をベースバンド周波数を持つダウンリンク信号に変換する。そして通信部３は、ダウンリンク信号に対して分離、復調、誤り訂正復号などの受信処理を行った後、そのダウンリンク信号を制御部２へ渡す。また通信部３は、制御部２から受け取ったアップリンク信号に対して誤り訂正符号化、変調及び多重化などの送信処理を行った後、そのアップリンク信号を無線周波数を持つ搬送波に重畳して基地局装置へ送信する。
【００１８】
マイクロホン４は、音声入力部の一例であり、携帯電話機１の周囲の音を集音し、その音の強度に応じたアナログ音声信号を生成する。マイクロホン４にて集音される音には、例えば、音源から直接マイクロホン４に達する直接音だけでなく、音源からの音声が携帯電話機１の周囲の壁などによって反射されてからマイクロホン４に間接的に達する残響音が含まれることがある。なお、音源は、例えば、携帯電話機１を用いて通話するユーザの口である。そしてマイクロホン４は、そのアナログ音声信号をアナログ／デジタル変換器５へ出力する。
【００１９】
アナログ／デジタル変換器５は、マイクロホン４から受け取ったアナログ音声信号を所定のサンプリングピッチでサンプリングすることによりデジタル化された入力音声信号を生成する。また、アナログ／デジタル変換器５は、増幅器を有し、アナログ音声信号を増幅した後にデジタル化してもよい。
アナログ／デジタル変換器５は、入力音声信号を残響抑制装置６へ出力する。
【００２０】
残響抑制装置６は、入力音声信号に含まれる残響成分を抑制する。そして残響抑制装置６は、残響成分が抑制された出力音声信号を制御部２へ出力する。なお、残響抑制装置６の詳細については後述する。
【００２１】
デジタル／アナログ変換器７は、制御部２から受け取った再生音声信号をデジタル−アナログ変換することでアナログ化する。なお、デジタル／アナログ変換器７は、増幅器を有し、その増幅器により、アナログ化された再生音声信号を増幅してもよい。そしてデジタル／アナログ変換器７は、アナログ化された再生音声信号をスピーカ８へ出力する。
スピーカ８は、音声出力部の一例であり、デジタル／アナログ変換器７から受け取った再生音声信号を再生する。
【００２２】
以下、残響抑制装置６の詳細について説明する。
図２は、第１の実施形態による残響抑制装置６の概略構成図である。残響抑制装置６は、残響特性算出部１１と、距離推定部１２と、残響特性補正部１３と、残響抑制部１４とを有する。
残響抑制装置６が有するこれらの各部は、それぞれ、別個の回路として残響抑制装置６に実装されてもよく、あるいはそれらの各部の機能を実現する一つの集積回路であってもよい。
【００２３】
残響特性算出部１１は、再生音声信号と、その再生音声信号がスピーカ８により再生され、マイクロホン４により集音されることにより得られた入力音声信号とに基づいて、スピーカ８からマイクロホン４への経路を含む音の経路のインパルス応答を求める。そして残響特性算出部１１は、そのインパルス応答に応じて、携帯電話機１の外部に位置する音源からの音の残響特性を近似的に求める。
本実施形態では、残響特性算出部１１は、有限インパルス応答型の適応フィルタを用いて、音の経路のインパルス応答を算出する。インパルス応答を表す適応フィルタの係数を、最小二乗法を用いて更新する場合、インパルス応答は次式に従って算出される。
【数１】
0005834948.tif 000002
ここで、s(t)は、再生音声信号を表し、x(t)は、入力音声信号を表す。そしてw(t)は音の経路のインパルス応答を表す。またe(t)は、誤差信号を表す。さらに係数μは、インパルス応答w(t)の更新速度を決めるための収束係数であり、例えば、0.01〜0.1に設定される。演算子'*'は、畳み込み演算を表す。
【００２４】
図３は、音の経路のインパルス応答の一例を示す図である。図３において、横軸は音がマイクロホン４に入力されてからの経過時間を表し、縦軸は音のパワーを表す。そしてグラフ３００がインパルス応答w(t)を表す。グラフ３００に示されるように、時刻0〜時刻t1までは音のパワーは急激に減少する。これは、音源から直接マイクロホン４に達する直接音の経路についてのインパルス応答に相当する。これに対し、時刻t1以降、時刻t2までの期間では、音源から発し、壁などで反射されてからマイクロホン４に達する残響音が含まれるため、パワーは緩やかに減少する。そして時刻t2において残響音がほとんど消失する。
【００２５】
なお、時刻t1は、例えば、以下のように、予め実験的に決定される。残響がほとんど生じない環境下に携帯電話機１を置いて（１）式により求めたインパルス応答w(t)に対して、残響が生じ易い環境下に環境下に携帯電話機１を置いて（１）式により求めたインパルス応答w(t)が大きくなる最初の時刻がt1として決定される。また時刻t2も、残響が生じ易い環境下に環境下に携帯電話機１を置いて（１）式により求めたインパルス応答w(t)がほぼ無視できる大きさとなる最初の時刻として実験的に定められる。例えば、時刻t1=50msec、時刻t2=400msecである。
【００２６】
上記のように、インパルス応答の成分のうち、時刻t1までの成分については直接音が支配的であり、時刻t1以降の成分については残響音が支配的である。そこで残響特性算出部１１は、インパルス応答w(t)のうちの時刻0〜t1までの成分w1(t)を直接音の成分、時刻t1〜t2までの成分w2(t)を残響音の成分とする。そして残響特性算出部１１は、直接音の成分w1(t)及び残響音の成分w2(t)をそれぞれ時間周波数変換することにより、直接音スペクトルW1(f)及び残響音スペクトルW2(f)を求める。なお、時間周波数変換として、残響特性算出部１１は、例えば、高速フーリエ変換または修正離散コサイン変換を用いることができる。
【００２７】
残響特性算出部１１は、直接音スペクトルW1(f)の実数成分の自乗と虚数成分の自乗とを加算することで、直接音のパワー|W1(f)|2を求める。同様に、残響特性算出部１１は、残響音スペクトルW2(f)の実数成分の自乗と虚数成分の自乗とを加算することで、残響音のパワー|W2(f)|2を求める。そして残響特性算出部１１は、次式に従って直接音スペクトルの振幅の絶対値と残響音スペクトルの振幅の絶対値の比を求めることで残響特性H(f)を算出する。
【数２】
0005834948.tif 000003
【００２８】
残響特性算出部１１は、インパルス応答w(t)を更新する度に、残響特性H(f)も更新する。そして残響特性算出部１１は、残響特性H(f)を残響特性補正部１３へ出力する。
【００２９】
距離推定部１２は、再生音声信号の音量が無視できるほど小さいときの入力音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて、マイクロホン４から音源までの距離を推定する。
【００３０】
一般に、音源がマイクロホンに近づくほど、マイクロホンで集音されるその音源からの音のパワーが大きくなる。そこで距離推定部１２は、次式に従って、再生音声信号及び入力音声信号の音量を、それぞれ、所定長を持つフレームごとに算出する。
【数３】
0005834948.tif 000004
ここでLは一つのフレームに含まれるサンプリング点の数である。サンプリング点の数Lは、例えば、16msec、32msecあるいは64msecを持つフレームに対して、128、256あるいは1024に設定される。またnは、0以上の整数であり、入力音声信号が入力開始されてからのフレーム数を表す。x(t)、s(t)は、それぞれ、入力音声信号、再生音声信号を表す。そしてPx(n)、Ps(n)は、それぞれ、フレームnにおける入力音声信号の音量、再生音声信号の音量を表す。
【００３１】
また、音源がマイクロホンに近づくほど、マイクロホンで集音されるその音源からの音のスペクトルのパワーに占める、低域成分のパワーの比が大きくなることが、マイクロホンの近接効果として知られている。そのため、距離推定部１２は、このパワーの比をマイクロホン４から音源までの距離を推定するために利用できる。
【００３２】
そこで距離推定部１２は、上記の比を求めるために、フレームごとに、入力音声信号を周波数領域の入力音スペクトルに変換する。なお、入力音スペクトルを算出するために、距離推定部１２は、例えば、高速フーリエ変換あるいは修正離散コサイン変換といった時間周波数変換を用いることができる。
【００３３】
距離推定部１２は、さらに、入力音スペクトルに基づいて、次式に従って入力音の全周波数帯域のパワースペクトルを算出する。
【数４】
0005834948.tif 000005
ここでX(n,f)は、n番目のフレームについての入力音スペクトルであり、PX(n,f)は、n番目のフレームの入力音のパワースペクトルである。またfは周波数を表す。
【００３４】
距離推定部１２は、全周波数帯域についての入力音のパワースペクトルの平均値に対する、低周波数帯域についての入力音のパワースペクトルの平均値の比である低域パワー比率rp(n)を、例えば、次式に従って算出する。
【数５】
0005834948.tif 000006
ここで、Mは、周波数帯域の総数であり、mlowは、低周波数帯域の上限の周波数帯域を表す。例えば、入力音スペクトルが128個の周波数帯域の成分で表される場合(M=128)、mlowは16に設定される。すなわち、入力音スペクトルが0〜4kHzの周波数成分を有する場合、低周波数帯域は、0〜500Hzの周波数帯域に相当する。mlowは、例えば、マイクロホン４の近接効果によって入力音声信号のうちのパワーが増幅する周波数成分に対応する周波数の上限値に相当する値に設定される。
【００３５】
距離推定部１２は、フレームごとに、再生音声信号の音量Ps(n)が所定の音量閾値Thp未満か否か判定する。そして距離推定部１２は、音量Ps(n)が音量閾値Thp以上である場合、マイクロホン４から音源までの距離を推定しない。これにより、距離推定部１２は、スピーカ８により再生された再生音声信号がマイクロホン４によって集音されたことによってマイクロホン４から音源までの距離の推定値の誤差が大きくなることを防止できる。
【００３６】
一方、距離推定部１２は、音量Ps(n)が音量閾値Thp未満である場合、入力音声信号の音量Px(n)及び低域パワー比率rp(n)と、マイクロホン４から音源までの距離との関係を表す参照テーブルを参照する。これにより、距離推定部１２は、そのフレームの入力音声信号の音量Px(n)及び低域パワー比率rp(n)に対応するマイクロホン４から音源までの距離を推定できる。なお、音量閾値Thpは、スピーカ８から出力された再生音声信号がマイクロホン４を介して入力されても、音源までの距離の推定に影響しない再生音声信号の音量の上限値、例えば、-48dBovに設定される。
【００３７】
図４は、入力音量及び低域パワー比率と距離との関係を示す参照テーブルの一例である。参照テーブル４００において、一番上の行４０１の各欄は入力音声信号の音量Px(n)の値の範囲を表し、二番目の行４０２の各欄は低域パワー比率rp(n)の値の範囲を表す。そして三番目の行４０３の各欄は、その欄と同じ列に示された入力音声信号の音量Px(n)の値の範囲及び低域パワー比率rp(n)の値の範囲に対応する、マイクロホン４から音源までの推定距離を表す。なお、参照テーブル４００は、例えば、距離推定部１２が有する、不揮発性の半導体メモリに予め記憶される。
例えば、入力音声信号の音量Px(n)が-25dBovであり、かつ、低域パワー比率rp(n)が1dBであれば、参照テーブル４００を参照することにより、マイクロホン４から音源までの距離は10cmと推定される。また、入力音声信号の音量Px(n)が-31dBovであり、かつ、低域パワー比率rp(n)が-5dBであれば、参照テーブル４００を参照することにより、マイクロホン４から音源までの距離は40cmと推定される。
【００３８】
なお、距離推定部１２は、マイクロホン４の近接効果による低域パワー比率の変化が大きい場合、低域パワー比率のみに基づいてマイクロホン４から音源までの距離を推定してもよい。この場合も、距離推定部１２は、低域パワー比率とマイクロホン４から音源までの距離との関係を表す参照テーブルを参照することにより、低域パワー比率に応じてマイクロホン４から音源までの距離を推定すればよい。
【００３９】
距離推定部１２は、マイクロホン４から音源までの推定距離を求める度に、その推定距離を残響特性補正部１３に通知する。
【００４０】
残響特性補正部１３は、マイクロホン４から音源までの推定距離に応じて残響特性を補正する。
【００４１】
図５（ａ）及び図５（ｂ）は、それぞれ、音源からマイクロホンに達する音の経路の一例を示す図である。図５（ａ）では、音源５０１から発し、マイクロホン４に入力される直接音の経路５０２の距離は0.1mである。一方、音の経路のインパルス応答に含まれる直接音の成分と残響音の成分との区切りの時刻t1が上記のように50msecである場合、空気中の音速は約340m/secであるので、音源５０１から発した音は時刻t1までに約17m進む。そのため、音源５０１から発し、マイクロホン４の周囲の壁面５０３で１回以上反射してからマイクロホン４に入力される残響音の経路５０４の長さは約17.1mとなる。
【００４２】
一方、図５（ｂ）では、音源５１１から発し、マイクロホン４に入力される直接音の経路５１２の距離は1mである。また時刻t1が上記のように50msecである場合、音源５０１から発し、マイクロホン４の周囲の壁面５１３で１回以上反射してからマイクロホン４に入力される残響音の経路５１４の長さは約18mとなる。
【００４３】
ここで、点音源からの距離r0の地点における、その点音源からの音に対する、点音源からの距離r1(ただし、r1>r0)の地点における、その点音源からの音の減衰量Aは、次式で表される。
【数６】
0005834948.tif 000007
【００４４】
そこで、図５（ａ）に示された音源の位置から図５（ｂ）に示された音源の位置へ音源が移動した場合、マイクロホン４に達する直接音は、（６）式により、-20dB減衰する。一方、図５（ａ）に示された音源の位置から図５（ｂ）に示された音源の位置へ音源が移動しても、マイクロホン４に達する残響音は、（６）式により、-0.45dBしか減少しない。したがって、マイクロホン４から音源までの距離が長いほど、残響成分の影響が大きくなる。
したがって、残響特性補正部１３は、マイクロホン４から音源までの推定距離が長くなるほど、残響特性が大きくなるように残響特性を補正する。
【００４５】
残響特性補正部１３は、 マイクロホン４から音源までの推定距離が長くなるほど、値が小さくなり、かつ正の値を持つ補正ゲインαを求める。例えば、マイクロホン４から音源までの推定距離と補正ゲインαとの関係を表す関係式の係数が、残響特性補正部１３が有する不揮発性の半導体メモリに予め記憶される。そして残響特性補正部１３は、その関係式にしたがって、その推定距離に対応する補正ゲインαの値を求める。あるいは、マイクロホン４から音源までの推定距離と補正ゲインαとの関係を表す参照テーブルが、残響特性補正部１３が有する不揮発性の半導体メモリに予め記憶される。そして残響特性補正部１３は、その参照テーブルを参照することにより、推定距離に対応する補正ゲインαの値を求めてもよい。
【００４６】
図６（ａ）及び図６（ｂ）は、それぞれ、マイクロホンから音源までの推定距離と補正ゲインαとの関係の一例を示す図である。図６（ａ）において、横軸はマイクロホン４から音源までの推定距離を表し、縦軸は補正ゲインαの値を表す。そしてグラフ６０１は、それぞれ、推定距離に対する補正ゲインαの値を表す。グラフ６０１に示されるように、推定距離が長くなるにつれて単調減少するように、補正ゲインαの値は決定される。
【００４７】
図６（ｂ）は、周波数帯域ごとに補正ゲインαを決定する場合の推定距離と補正ゲインαの関係を示す。図６（ｂ）においても、横軸はマイクロホン４から音源までの推定距離を表し、縦軸は補正ゲインαの値を表す。そしてグラフ６０２〜６０４は、それぞれ、周波数50Hz、100Hz、及び1kHzにおける、推定距離に対する補正ゲインαの値を表す。マイクロホン４から音源までの距離が近いと、直接音に占める低域成分がより増幅されるので、残響音の影響も低下する。そこでグラフ６０２〜６０４に示されるように、残響特性補正部１３は、周波数が低いほど、補正ゲインαの値を大きくする。
【００４８】
残響特性補正部１３は、次式に従って、推定された残響特性を補正する。
【数７】
0005834948.tif 000008
ここで、H(f)は、残響特性算出部１１により推定された残響特性であり、H'(f)は、補正後の残響特性である。またαは補正ゲインである。（７）式より、マイクロホン４から音源までの推定距離が長くなり、補正ゲインαが小さくなるほど、補正後の残響特性が大きくなることが分かる。
【００４９】
なお、残響特性補正部１３は、一旦特定した補正ゲインαを、残響特性補正部１３が有する揮発性または不揮発性の半導体メモリに記憶しておいてもよい。そして残響特性補正部１３は、次に距離推定部１２から推定距離が通知され、補正ゲインαが再度求められるまで、記憶された補正ゲインαを用いて残響特性を補正してもよい。
残響特性補正部１３は、補正された残響特性H'(f)を残響抑制部１４へ出力する。
【００５０】
残響抑制部１４は、残響特性補正部１３から受け取った、補正された残響特性H'(f)に基づいて、入力音声信号に含まれる残響成分を抑制する。
残響抑制部１４は、フレームごとに、入力音声信号を周波数領域の入力音スペクトルに変換する。なお、入力音スペクトルを算出するために、残響抑制部１４は、距離推定部１２と同様に、例えば、高速フーリエ変換あるいは修正離散コサイン変換を用いることができる。
【００５１】
次に、残響抑制部１４は、次式のように、入力音スペクトルに補正された残響特性を乗じることにより、残響音のスペクトルを推定する。
【数８】
0005834948.tif 000009
ここで、X(n,f)は、n番目のフレームについての入力音スペクトルであり、H'(f)は、補正された残響特性である。そしてR(n,f)は、n番目のフレームについての推定された残響音のスペクトルである。
【００５２】
残響抑制部１４は、入力音スペクトルから推定された残響音のスペクトルを減算した差分値(δ(n,f)=X(n,f)-R(n,f))に基づいて、入力音スペクトルの減衰量を表す抑制ゲインG(f)を決定する。例えば、残響抑制部１４は、その差分値δ(n,f)が小さいほど、すなわち、入力音スペクトルに占める残響成分が多いほど、抑制ゲインを大きくする。
【００５３】
図７は、入力音スペクトルと推定された残響のスペクトル間の差分値δ(n,f)に対する抑制ゲインG(f)の関係の一例を示す図である。図７において、横軸は差分値δ(n,f)を表し、縦軸は抑制ゲインG(f)を表す。グラフ７００は、差分値δ(n,f)と抑制ゲインG(f)の関係を表す。
グラフ７００に示されるように、例えば、差分値δ(n,f)が6dB以下であれば、抑制ゲインG(f)は10dBに設定される。そして差分値δ(n,f)が6dBから12dBの間であれば、その差分値δ(n,f)が増加するにつれて抑制ゲインG(f)は線形に減少し、その差分値δ(n,f)が12dB以上となると、抑制ゲインG(f)は0dB、すなわち、入力音スペクトルを全く減衰させない値に設定される。
【００５４】
例えば、差分値δ(n,f)と抑制ゲインG(f)との関係を表す関係式の係数が、残響抑制部１４が有する不揮発性の半導体メモリに予め記憶される。そして残響抑制部１４は、その関係式にしたがって、差分値δ(n,f)に対応する抑制ゲインG(f)を決定する。あるいは、差分値δ(n,f)と抑制ゲインG(f)との関係を表す参照テーブルが、残響抑制部１４が有する不揮発性の半導体メモリに予め記憶される。そして残響抑制部１４は、その参照テーブルを参照することにより、差分値δ(n,f)に対応する抑制ゲインG(f)を決定してもよい。
【００５５】
残響抑制部１４は、例えば、次式に従って、入力音スペクトルのうちの残響成分を抑制したスペクトルである出力音スペクトルを算出する。
【数９】
0005834948.tif 000010
ここで、Y(n,f)は出力音スペクトルである。（９）式から明らかなように、抑制ゲインG(f)が大きいほど、出力音スペクトルY(n,f)も小さくなる。
【００５６】
残響抑制部１４は、出力音スペクトルY(n,f)を、入力音スペクトルを算出する際に用いた時間周波数変換の逆変換を用いて時間領域の信号に変換することにより、出力音声信号を得る。
そして残響抑制部１４は、出力音声信号を携帯電話機１の制御部２へ出力する。
【００５７】
図８は、残響抑制装置６により実行される残響抑制処理の動作フローチャートである。
残響特性算出部１１は、再生音声信号と、スピーカ８で再生された再生音声信号をマイクロホン４により集音して得られた入力音声信号とから、適応フィルタを用いてスピーカ８からマイクロホン４の音の経路のインパルス応答w(t)を求める（ステップＳ１０１）。残響特性算出部１１は、インパルス応答w(t)のうち、直接音に相当する成分と残響音に相当する成分とをそれぞれ時間周波数変換して得られた直接音スペクトルW1(f)と残響音スペクトルW2(f)の絶対値の比|W2(f)/W1(f)|を残響特性H(f)として推定する（ステップＳ１０２）。残響特性算出部１１は、残響特性H(f)を残響特性補正部１３へ出力する。
【００５８】
距離推定部１２は、再生音声信号が無視できるほど小さいときの、携帯電話機１の外部の音源からの音をマイクロホン４で集音した入力音声信号の音量と周波数特性とを求める（ステップＳ１０３）。そして距離推定部１２は、その音量と周波数特性に基づいてマイクロホン４から音源までの距離を推定する（ステップＳ１０４）。距離推定部１２は、その推定距離を残響特性補正部１３へ通知する。
【００５９】
残響特性補正部１３は、推定距離が長くなるほど、残響特性H(f)が大きくなるように、残響特性H(f)を補正する（ステップＳ１０５）。そして残響特性補正部１３は、補正された残響特性H'(f)を残響抑制部１４へ出力する。
【００６０】
残響抑制部１４は、入力音声信号を時間周波数変換することで入力音スペクトルを求め、その入力音スペクトルに補正された残響特性H'(f)を乗じることで残響音のスペクトルを推定する（ステップＳ１０６）。そして残響抑制部１４は、入力音スペクトルと推定された残響音スペクトルの差が小さいほど、抑制ゲインG(f)を大きな値に設定する（ステップＳ１０７）。残響抑制部１４は、抑制ゲインG(f)が大きいほど、入力音スペクトルを減衰させることで、出力音スペクトルを算出する（ステップＳ１０８）。そして残響抑制部１４は、出力音スペクトルを周波数時間変換することで、出力音声信号を得る（ステップＳ１０９）。そして残響抑制装置６は、その出力音声信号を出力し、残響抑制処理を終了する。
【００６１】
以上に説明してきたように、この残響抑制装置は、マイクロホンから音源までの距離を推定し、その推定距離に応じて、再生音声信号に基づいて近似的に算出された残響特性を補正する。そのため、この残響抑制装置は、マイクロホンから音源までの距離が変化しても、残響特性を適切に求めることができるので、入力音声信号中の残響成分を抑制できる。またこの残響抑制装置は、マイクロホンから音源までの距離を推定するのに入力音声信号の音量及び周波数特性を利用するので、１本のマイクロホンだけでその距離を推定できる。
【００６２】
次に、第２の実施形態による残響抑制装置について説明する。第２の実施形態による残響抑制装置は、その残響抑制装置が実装された携帯電話機の制御部から、接話状態かハンズフリー状態かを表す情報を取得する。そしてこの残響抑制装置は、接話状態における入力音声信号の音量を基準として、ハンズフリー状態における入力音声信号の音量との差に基づいてマイクロホンから音源までの距離を推定する。
【００６３】
第２の実施形態による残響抑制装置は、第１の実施形態による残響抑制装置と比較して、距離推定部１２の処理が異なる。そこで以下では、距離推定部１２について説明する。第２の実施形態による残響抑制装置のその他の構成要素については、第１の実施形態による残響抑制装置の対応する構成要素についての説明を参照されたい。また以下では、第２の実施形態による残響抑制装置も、図１に示された携帯電話機１に実装されたものとして説明する。
【００６４】
携帯電話機１の制御部２は、キーパッドまたはタッチパネルといった携帯電話機１の操作部（図示せず）を介した操作により、通話が開始されると、その開始時点での通話距離設定を表す通話距離設定情報を残響抑制装置６へ通知する。なお、通話距離設定として、例えば、接話設定とハンズフリー設定とが含まれる。接話設定は、スピーカ８から出力される音の音量が、ユーザがスピーカ８に耳を近づけなければ聞こえないほど小さく設定されることで、ユーザにマイクロホン４に近接した状態で通話させる設定である。一方、ハンズフリー設定は、スピーカ８から出力される音の音量が大きく、ユーザが携帯電話機１から離れていても通話できる設定である。そのため、ハンズフリー設定では、ユーザとマイクロホン４間の距離は、接話設定時のその距離よりも大きくなる可能性がある。
また制御部２は、操作部を介して接話設定とハンズフリー設定とが切り替えられると、切替後の通話距離設定を表す通話距離設定情報を残響抑制装置６へ通知する。
【００６５】
距離推定部１２は、第１の実施形態と同様に、入力音声信号の音量を、（３）式に従って、所定長を持つフレームごとに算出する。
また距離推定部１２は、制御部２から通知された通話距離設定情報により、接話設定がなされていると判断している期間中、マイクロホン４から音源までの推定距離として、予め設定された値、例えば、5cmを残響特性補正部１３へ出力する。この値は、例えば、複数のユーザのそれぞれについて、接話時における、マイクロホン４からそのユーザの口までの距離を測定し、その測定値の平均値を求めることにより決定され、距離推定部１２が有する不揮発性の半導体メモリに予め記憶される。
【００６６】
さらに、距離推定部１２は、接話設定がなされていると判断している期間中、次式に従って、接話時の入力音声信号の平均音量を算出する。接話時の入力音声信号の平均音量は、ハンズフリー時における、マイクロホン４から音源までの距離を推定するための基準値として利用される。
【数１０】
0005834948.tif 000011
ここで、Pca(n)、Pca(n-1)は、それぞれ、n番目、(n-1)番目のフレームにおける、接話時の入力音声信号の平均音量である。またPx(n)は、n番目のフレームにおける、接話時の入力音声信号の音量である。そしてcは、平滑化係数であり、例えば、0.01〜0.1に設定される。
距離推定部１２は、接話時の入力音声信号の平均音量が更新される度に、その平均音量を距離推定部１２が有する半導体メモリに記憶する。
【００６７】
距離推定部１２は、制御部２から通知された通話距離設定情報により、ハンズフリー設定がなされていると判断している期間中、メモリに記憶しておいた接話時の入力音声信号の平均音量からフレームごとの入力音声信号の音量を引いた差分音量Pdを求める。そして距離推定部１２は、差分音量Pdとマイクロホン４から音源までの距離の関係を表す参照テーブルを参照することにより、その差分音量Pdに対応する、マイクロホン４から音源までの推定距離を決定する。なお、参照テーブルは、例えば、距離推定部１２が有する、不揮発性の半導体メモリに予め記憶される。
距離推定部１２は、その推定距離を残響特性補正部１３へ出力する。
【００６８】
図９は、差分音量Pdとマイクロホン４から音源までの距離の関係の一例を表す図である。図９において、横軸は差分音量Pdを表し、縦軸はマイクロホン４から音源までの推定距離を表す。そしてグラフ９００は、差分音量Pdと推定距離との関係を表すグラフである。グラフ９００に示されるように、差分音量Pdが大きくなるほど、すなわち、接話時の入力音声信号の平均音量に比べてハンズフリー時の入力音声信号の音量が小さいほど、推定距離が長くなる。
【００６９】
図１０は、距離推定部１２により実行される距離推定処理の動作フローチャートである。この距離推定処理は、図８に示された残響抑制処理の動作フローチャートにおける、ステップＳ１０３及びＳ１０４の代わりに実行される。
距離推定部１２は、現在のフレームについて入力音声信号の音量を求める（ステップＳ２０１）。そして距離推定部１２は、携帯電話機１の制御部２から受け取った通話距離設定情報に基づいて、現在の通話距離設定が接話設定か否か判定する（ステップＳ２０２）。現在の通話距離設定が接話設定である場合（ステップＳ２０２−Ｙｅｓ）、距離推定部１２は、接話設定に対して予め定められた距離を残響特性補正部１３へ出力する（ステップＳ２０３）。また距離推定部１２は、現在のフレームについての入力音声信号の音量と、直前までのフレームにおける、接話時の入力音声信号の平均音量とに基づいて、その平均音量を更新し、その平均音量を記憶する（ステップＳ２０４）。
【００７０】
一方、現在の通話距離設定がハンズフリー設定である場合（ステップＳ２０２−Ｎｏ）、距離推定部１２は、接話時の入力音声信号の平均音量から現在のフレームについての入力音声信号の音量を減じた差分音量を求める（ステップＳ２０５）。距離推定部１２は、差分音量と、マイクロホン４から音源までの距離との関係を表す参照テーブルを参照することにより、ステップＳ２０５で求められた差分音量に対応するマイクロホン４から音源までの推定距離を決定し、その推定距離を残響特性補正部１３へ出力する（ステップＳ２０６）。
【００７１】
ステップＳ２０４またはＳ２０６の後、距離推定部１２は、距離推定処理を終了する。
【００７２】
この第２の実施形態によれば、接話時の入力音声信号の平均音量とハンズフリー時の入力音声信号の音量との差に基づいてマイクロホンから音源までの距離を推定するので、ユーザの発声音量の個人差によらず、マイクロホンから音源までの距離を適切に推定できる。
【００７３】
なお、変形例によれば、残響抑制装置は、距離推定部及び残響抑制部とは別個に、入力音声信号をフレーム単位で時間周波数変換することにより入力音スペクトルを算出する時間周波数変換部（図示せず）を有してもよい。この場合、時間周波数変換部により生成されたフレームごとの入力音スペクトルは、距離推定部及び残響抑制部に入力される。これにより、入力音スペクトルを算出するための演算量または回路規模が削減可能となる。
【００７４】
また他の変形例によれば、残響特性算出部は、（３）式に従って再生音声信号の音量を算出し、その音量が音量閾値Thp以上である場合に限り、残響特性を算出してもよい。この場合、フレームごとに、図８に示した動作フローチャートにおける、ステップＳ１０１、Ｓ１０２とステップＳ１０３、Ｓ１０４とが再生音声信号の音量に応じて選択的に実行される。
【００７５】
さらに、上記の各実施形態による残響抑制装置の各部が有する各機能をコンピュータに実現させるコンピュータプログラムは、磁気記録媒体あるいは光記録媒体といった、コンピュータによって読み取り可能な媒体に記録された形で提供されてもよい。
【００７６】
図１１は、上記の実施形態またはその変形例による残響抑制装置の各部の機能を実現するコンピュータプログラムが動作することにより、残響抑制装置として動作するコンピュータの構成図である。
コンピュータ１００は、ユーザインターフェース部１０１と、通信インターフェース部１０２と、記憶部１０３と、記憶媒体アクセス装置１０４と、プロセッサ１０５とを有する。プロセッサ１０５は、ユーザインターフェース部１０１、通信インターフェース部１０２、記憶部１０３及び記憶媒体アクセス装置１０４と、例えば、バスを介して接続される。
【００７７】
ユーザインターフェース部１０１は、例えば、キーボードとマウスなどの入力装置と、液晶ディスプレイといった表示装置とを有する。または、ユーザインターフェース部１０１は、タッチパネルディスプレイといった、入力装置と表示装置とが一体化された装置を有してもよい。そしてユーザインターフェース部１０１は、例えば、ユーザの操作に応じて、残響抑制処理を開始させる操作信号をプロセッサ１０５へ出力する。
【００７８】
通信インターフェース部１０２は、コンピュータ１００を、マイクロホン及びスピーカと接続するためのオーディオインターフェース及びその制御回路を有してもよい。
さらに、通信インターフェース部１０２は、イーサネット（登録商標）などの通信規格に従った通信ネットワークに接続するための通信インターフェース及びその制御回路を有してもよい。
この場合には、通信インターフェース部１０２は、通信ネットワークに接続された他の機器から、再生音声信号を取得し、プロセッサ１０５へ渡す。また通信インターフェース部１０２は、プロセッサ１０５から受け取った、残響抑制された出力音声信号を通信ネットワークを介して他の機器へ出力してもよい。
【００７９】
記憶部１０３は、例えば、読み書き可能な半導体メモリと読み出し専用の半導体メモリとを有する。そして記憶部１０３は、プロセッサ１０５上で実行される、残響抑制処理を実行するためのコンピュータプログラム、及び残響抑制処理で利用される、各参照テーブルといったデータを記憶する。
【００８０】
記憶媒体アクセス装置１０４は、例えば、磁気ディスク、半導体メモリカード及び光記憶媒体といった記憶媒体１０６にアクセスする装置である。記憶媒体アクセス装置１０４は、例えば、記憶媒体１０６に記憶されたプロセッサ１０５上で実行される、残響抑制処理用のコンピュータプログラムを読み込み、プロセッサ１０５に渡す。
【００８１】
プロセッサ１０５は、上記の各実施形態の何れかまたは変形例による残響抑制処理用コンピュータプログラムを実行することにより、マイクロホンを介して入力された音声信号に含まれる残響成分を抑制する。そしてプロセッサ１０５は、残響抑制された入力音声信号を通信インターフェース部１０２を介して他の機器へ出力する。
【００８２】
ここに挙げられた全ての例及び特定の用語は、読者が、本発明及び当該技術の促進に対する本発明者により寄与された概念を理解することを助ける、教示的な目的において意図されたものであり、本発明の優位性及び劣等性を示すことに関する、本明細書の如何なる例の構成、そのような特定の挙げられた例及び条件に限定しないように解釈されるべきものである。本発明の実施形態は詳細に説明されているが、本発明の精神及び範囲から外れることなく、様々な変更、置換及び修正をこれに加えることが可能であることを理解されたい。
【００８３】
以上説明した実施形態及びその変形例に関し、更に以下の付記を開示する。
（付記１）
第１の音声信号と、音声出力部により再生された該第１の音声信号を音声入力部により集音した音を表す第２の音声信号とから、該音声出力部から該音声入力部への音の経路のインパルス応答を求め、該インパルス応答に応じて残響特性を算出する残響特性算出部と、
音源からの音を前記音声入力部により集音した音を表す第３の音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて前記音声入力部から前記音源までの距離を推定する距離推定部と、
前記距離の推定値が長くなるほど、前記残響特性を大きくするよう補正する残響特性補正部と、
前記補正された残響特性に応じて前記第３の音声信号に含まれる残響成分の周波数スペクトルを推定し、前記第３の音声信号の周波数スペクトルと該残響成分の周波数スペクトルの差が小さいほど前記第３の音声信号を減衰させる残響抑制部と、
を有する残響抑制装置。
（付記２）
前記残響特性算出部は、前記インパルス応答に含まれる、前記音声出力部から前記音声入力部へ直接達する直接音に対応する期間の成分のスペクトルに対する前記音声出力部から前記音声入力部へ間接的に達する残響音に対応する期間の成分のスペクトルの比を前記残響特性とする、付記１に記載の残響抑制装置。
（付記３）
前記距離推定部は、前記第３の音声信号の周波数特性として、前記第３の音声信号が持つ周波数帯域全体のパワースペクトルに対する所定の周波数以下の周波数帯域のパワースペクトルの比を算出し、当該比が大きいほど前記距離の推定値を短くする、付記１または２に記載の残響抑制装置。
（付記４）
前記所定の周波数は、前記音声入力部に対する近接効果により、前記第３の音声信号のうちのパワーが増幅する周波数成分に対応する周波数の上限値である、付記３に記載の残響抑制装置。
（付記５）
前記距離推定部は、前記音声入力部を有する機器について接話設定されている第１の期間における前記第３の音声信号の音量の平均値を算出し、前記機器が接話設定されていない第２の期間における前記第３の音声信号の音量を前記平均値から減じた値が大きいほど、当該第２の期間における前記距離の推定値を長くする、付記１または２に記載の残響抑制装置。
（付記６）
前記残響特性補正部は、第１の周波数における前記残響特性を除する係数を、前記第１の周波数よりも高い第２の周波数における前記残響特性を除する係数よりも大きくする、付記１〜５の何れか一項に記載の残響抑制装置。
（付記７）
第１の音声信号と、音声出力部により再生された該第１の音声信号を音声入力部により集音した音を表す第２の音声信号とから、該音声出力部から該音声入力部への音の経路のインパルス応答を求め、該インパルス応答に応じて残響特性を算出し、
音源からの音を前記音声入力部により集音した音を表す第３の音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて前記音声入力部から前記音源までの距離を推定し、
前記距離の推定値が長くなるほど、前記残響特性を大きくするよう補正し、
前記補正された残響特性に応じて前記第３の音声信号に含まれる残響成分の周波数スペクトルを推定し、
前記第３の音声信号の周波数スペクトルと該残響成分の周波数スペクトルの差が小さいほど前記第３の音声信号を減衰させる、
ことを含む残響抑制方法。
（付記８）
第１の音声信号と、音声出力部により再生された該第１の音声信号を音声入力部により集音した音を表す第２の音声信号とから、該音声出力部から該音声入力部への音の経路のインパルス応答を求め、該インパルス応答に応じて残響特性を算出し、
音源からの音を前記音声入力部により集音した音を表す第３の音声信号の音量及び周波数特性のうちの少なくとも一方に基づいて前記音声入力部から前記音源までの距離を推定し、
前記距離の推定値が長くなるほど、前記残響特性を大きくするよう補正し、
前記補正された残響特性に応じて前記第３の音声信号に含まれる残響成分の周波数スペクトルを推定し、
前記第３の音声信号の周波数スペクトルと該残響成分の周波数スペクトルの差が小さいほど前記第３の音声信号を減衰させる、
ことをコンピュータに実行させる残響抑制用コンピュータプログラム。
【符号の説明】
【００８４】
１    携帯電話機
２    制御部
３    通信部
４    マイクロホン
５    アナログ／デジタル変換器
６    残響抑制装置
７    デジタル／アナログ変換器
８    スピーカ
１１    残響特性算出部
１２    距離推定部
１３    残響特性補正部
１４    残響抑制部
１００    コンピュータ
１０１    ユーザインターフェース部
１０２    通信インターフェース部
１０３    記憶部
１０４    記憶媒体アクセス装置
１０５    プロセッサ
１０６    記憶媒体
【図１】
0005834948.tif 000012
【図２】
0005834948.tif 000013
【図３】
0005834948.tif 000014
【図４】
0005834948.tif 000015
【図５】
0005834948.tif 000016
【図６】
0005834948.tif 000017
【図７】
0005834948.tif 000018
【図８】
0005834948.tif 000019
【図９】
0005834948.tif 000020
【図１０】
0005834948.tif 000021
【図１１】
0005834948.tif 000022