JP3420705B2 - エコー抑圧方法及び装置並びにエコー抑圧プログラムが記憶されたコンピュータに読取り可能な記憶媒体 - Google Patents
エコー抑圧方法及び装置並びにエコー抑圧プログラムが記憶されたコンピュータに読取り可能な記憶媒体Info
- Publication number
- JP3420705B2 JP3420705B2 JP25441198A JP25441198A JP3420705B2 JP 3420705 B2 JP3420705 B2 JP 3420705B2 JP 25441198 A JP25441198 A JP 25441198A JP 25441198 A JP25441198 A JP 25441198A JP 3420705 B2 JP3420705 B2 JP 3420705B2
- Authority
- JP
- Japan
- Prior art keywords
- echo
- signal
- short
- transmission signal
- time spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Circuit For Audible Band Transducer (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Telephonic Communication Services (AREA)
Description
【０００１】
【発明の属する技術分野】この発明は、エコー抑圧方法
及び装置並びにエコー抑圧プログラムが記憶されたコン
ピュータに読取り可能な記憶媒体に係わり、例えば２線
４線変換系及び拡声通話系などにおいて、ハウリングの
原因及び聴覚上の障害となるエコー信号を抑圧し、送話
者音声信号を強調する技術に関する。
及び装置並びにエコー抑圧プログラムが記憶されたコン
ピュータに読取り可能な記憶媒体に係わり、例えば２線
４線変換系及び拡声通話系などにおいて、ハウリングの
原因及び聴覚上の障害となるエコー信号を抑圧し、送話
者音声信号を強調する技術に関する。
【０００２】
【従来の技術】まず、このようなハウリングの原因及び
聴覚上の障害となるエコー信号について、図８に示す拡
声通話系を参照して説明する。この図において、符号
１，３は送話用マイクロホン、２，４は受話スピーカ、
５，７は送話信号増幅器、６，８は受話信号増幅器、９
は伝送路、１０は送話者、１１は受話者である。
聴覚上の障害となるエコー信号について、図８に示す拡
声通話系を参照して説明する。この図において、符号
１，３は送話用マイクロホン、２，４は受話スピーカ、
５，７は送話信号増幅器、６，８は受話信号増幅器、９
は伝送路、１０は送話者、１１は受話者である。
【０００３】送話者１０が発声した送話音声は、送話用
マイクロホン１、送話信号増幅器５、伝送路９、受話信
号増幅器８及び受話スピーカ４を経て受話者１１に伝送
される。このような拡声通話系は、従来の電話通話系の
ように送受話器を手に持つ必要がないため、作業をしな
がらの通話が可能であり、また自然な対面通話が実現で
きるという長所を持ち、通信会議やテレビ電話、拡声電
話機などに広く利用が進められている。
マイクロホン１、送話信号増幅器５、伝送路９、受話信
号増幅器８及び受話スピーカ４を経て受話者１１に伝送
される。このような拡声通話系は、従来の電話通話系の
ように送受話器を手に持つ必要がないため、作業をしな
がらの通話が可能であり、また自然な対面通話が実現で
きるという長所を持ち、通信会議やテレビ電話、拡声電
話機などに広く利用が進められている。
【０００４】しかしながら、拡声通話系の欠点として、
エコーの存在が間題となっている。すなわち、スピーカ
４から受話側に伝わった音声が、マイクロホン３によっ
て受音され、送話信号増幅器７、伝送路９、受話信号増
幅器６及びスピーカ２を経て送話者１０側に伝送され
る。送話者１０とって、この現象は、自分の発声した音
声がスピーカ２から再生されるというエコー現象であ
り、音響エコーなどと呼ばれている。このエコー現象
は、拡声通話系において通話の障害や不快感などの悪影
響を生じる。
エコーの存在が間題となっている。すなわち、スピーカ
４から受話側に伝わった音声が、マイクロホン３によっ
て受音され、送話信号増幅器７、伝送路９、受話信号増
幅器６及びスピーカ２を経て送話者１０側に伝送され
る。送話者１０とって、この現象は、自分の発声した音
声がスピーカ２から再生されるというエコー現象であ
り、音響エコーなどと呼ばれている。このエコー現象
は、拡声通話系において通話の障害や不快感などの悪影
響を生じる。
【０００５】さらに、スピーカ２から発声された音声
は、マイクロホン１によって受音され、音声信号の閉ル
ープが形成される。そして、この閉ループのループゲイ
ンが１より大きい場合には、ハウリング現象が発生し
て、通話不能状態に至る。
は、マイクロホン１によって受音され、音声信号の閉ル
ープが形成される。そして、この閉ループのループゲイ
ンが１より大きい場合には、ハウリング現象が発生し
て、通話不能状態に至る。
【０００６】このような拡声通話系の問題点を克服する
ために、エコー消去装置（エコーキャンセラ）が利用さ
れている。図９は、従来のフルバンド方式のエコーキャ
ンセラの構成例を示すブロック図である。なお、この図
において、上記図８と同一の構成要素については同一符
号を付している。
ために、エコー消去装置（エコーキャンセラ）が利用さ
れている。図９は、従来のフルバンド方式のエコーキャ
ンセラの構成例を示すブロック図である。なお、この図
において、上記図８と同一の構成要素については同一符
号を付している。
【０００７】この図において、符号２１はエコーキャン
セラ、２２は疑似エコー経路畳み込み演算回路、２３は
エコー経路推定回路、２４は減算器、２５は非線形エコ
ー抑圧処理部である。また、ｘ（ｋ）は受話信号、ｈ
（ｋ）は受話スピーカ４と送話用マイクロホン３の問の
エコー経路伝達特性（インパルス応答）、ｙ（ｋ）はエ
コー信号、ｙe（ｋ）は疑似エコー信号，ｈe（ｋ）はエ
コー経路インパルス応答の推定値、ｂ（ｋ）は誤差信
号、ｓ（ｋ）は近端話者の送話信号、ｚ（ｋ）はマイク
ロホン出力信号を表している。
セラ、２２は疑似エコー経路畳み込み演算回路、２３は
エコー経路推定回路、２４は減算器、２５は非線形エコ
ー抑圧処理部である。また、ｘ（ｋ）は受話信号、ｈ
（ｋ）は受話スピーカ４と送話用マイクロホン３の問の
エコー経路伝達特性（インパルス応答）、ｙ（ｋ）はエ
コー信号、ｙe（ｋ）は疑似エコー信号，ｈe（ｋ）はエ
コー経路インパルス応答の推定値、ｂ（ｋ）は誤差信
号、ｓ（ｋ）は近端話者の送話信号、ｚ（ｋ）はマイク
ロホン出力信号を表している。
【０００８】エコーキャンセラ２１では、エコー経路推
定回路２３においてエコー経路のインパルス応答を推定
し、その推定値ｈe（ｋ）を疑似エコー経路畳み込み演
算回路２２に転送する。そして、疑似エコー経路畳み込
み演算回路２２において、上記推定値ｈe（ｋ）と受話
信号ｚ（ｋ）との畳み込み演算を実行して、疑似エコー
信号ｙe（ｋ）を合成し、減算器２４においてマイクロ
ホン３の出力信号から疑似エコー信号ｙe（ｋ）を減算
する。この場合、エコー経路インパルス応答の推定が良
好に行なわれていれば、エコー信号ｙ（ｋ）と疑似エコ
ー信号ｙe（ｋ）とはほぼ等しいものとなり、上記減算
の結果、マイクロホン３の出力に含まれるエコー信号ｙ
（ｋ）は除去される。
定回路２３においてエコー経路のインパルス応答を推定
し、その推定値ｈe（ｋ）を疑似エコー経路畳み込み演
算回路２２に転送する。そして、疑似エコー経路畳み込
み演算回路２２において、上記推定値ｈe（ｋ）と受話
信号ｚ（ｋ）との畳み込み演算を実行して、疑似エコー
信号ｙe（ｋ）を合成し、減算器２４においてマイクロ
ホン３の出力信号から疑似エコー信号ｙe（ｋ）を減算
する。この場合、エコー経路インパルス応答の推定が良
好に行なわれていれば、エコー信号ｙ（ｋ）と疑似エコ
ー信号ｙe（ｋ）とはほぼ等しいものとなり、上記減算
の結果、マイクロホン３の出力に含まれるエコー信号ｙ
（ｋ）は除去される。
【０００９】ここで、疑似エコー経路畳み込み演算回路
２２は、エコー経路インパルス応答ｈ（ｋ）の経時変動
に追従する必要がある。そのため、エコー経路推定回路
２３では、適応アルゴリズムを用いてエコー経路インパ
ルス応答ｈ（ｋ）の推定を行なう。この推定動作は、受
話状態、すなわちｓ（ｋ）≒０であり、ｚ（ｋ）≒ｙ
（ｋ）とみなせる時に実行される。
２２は、エコー経路インパルス応答ｈ（ｋ）の経時変動
に追従する必要がある。そのため、エコー経路推定回路
２３では、適応アルゴリズムを用いてエコー経路インパ
ルス応答ｈ（ｋ）の推定を行なう。この推定動作は、受
話状態、すなわちｓ（ｋ）≒０であり、ｚ（ｋ）≒ｙ
（ｋ）とみなせる時に実行される。
【００１０】適応アルゴリズムとは、疑似エコー経路に
おけるインパルス応答の推定値ｈe（ｋ）を求めるアル
ゴリズムであって、ＬＭＳ法、学習同定法あるいはＥＳ
法などが知られている。ここで、疑似エコー経路の値が
真のエコー経路の値に近く、疑似エコー信号ｙe（ｋ）
がエコー信号ｙ（ｋ）にほぼ等しくなった状態を「収束
した」と呼ぶ。また、疑似エコー経路畳み込み演算回路
２２とエコー経路推定回路２３を合わせて、ここでは適
応フィルタと呼ぶことにする。
おけるインパルス応答の推定値ｈe（ｋ）を求めるアル
ゴリズムであって、ＬＭＳ法、学習同定法あるいはＥＳ
法などが知られている。ここで、疑似エコー経路の値が
真のエコー経路の値に近く、疑似エコー信号ｙe（ｋ）
がエコー信号ｙ（ｋ）にほぼ等しくなった状態を「収束
した」と呼ぶ。また、疑似エコー経路畳み込み演算回路
２２とエコー経路推定回路２３を合わせて、ここでは適
応フィルタと呼ぶことにする。
【００１１】
【発明が解決しようとする課題】ところで、上記従来の
エコーキャンセラ（すなわち適応フィルタ）は、受話状
態すなわち近端話者送話信号が存在しないときにのみ、
エコー経路インパルス応答の推定を行なう。遠端話者と
近端話者が同時に会話するダブルトークの際、つまりエ
コー信号と送話信号が一緒にマイクロホン３に入力され
ている状態では、推定動作を行なうことができない。こ
のため、ダブルトーク中に人が動くなどのエコー経路の
変動があると、真のエコー信号を推定できないためエコ
ーを十分に消去することができず、結果として通話品質
に劣化させてしまうという問題点がある。
エコーキャンセラ（すなわち適応フィルタ）は、受話状
態すなわち近端話者送話信号が存在しないときにのみ、
エコー経路インパルス応答の推定を行なう。遠端話者と
近端話者が同時に会話するダブルトークの際、つまりエ
コー信号と送話信号が一緒にマイクロホン３に入力され
ている状態では、推定動作を行なうことができない。こ
のため、ダブルトーク中に人が動くなどのエコー経路の
変動があると、真のエコー信号を推定できないためエコ
ーを十分に消去することができず、結果として通話品質
に劣化させてしまうという問題点がある。
【００１２】また、適応フィルタは、エコー経路インパ
ルス応答の振幅及び位相の両方を推定するため、収束し
た場合にはエコーの十分な抑圧が可能となる反面、より
多くの情報を推定しなければならないために、先に述べ
たようなエコー経路が変動した場合に、瞬間的に消し切
れないエコー（残留エコー）が増大し、耳障りの原因と
なる。この残留エコーを抑圧するために、一般にエコー
キャンセラは先の適応フィルクに加えて、エコーサプレ
ッサやセンタクリッパのような非線形エコー抑圧処理部
２５を併せ持つ。
ルス応答の振幅及び位相の両方を推定するため、収束し
た場合にはエコーの十分な抑圧が可能となる反面、より
多くの情報を推定しなければならないために、先に述べ
たようなエコー経路が変動した場合に、瞬間的に消し切
れないエコー（残留エコー）が増大し、耳障りの原因と
なる。この残留エコーを抑圧するために、一般にエコー
キャンセラは先の適応フィルクに加えて、エコーサプレ
ッサやセンタクリッパのような非線形エコー抑圧処理部
２５を併せ持つ。
【００１３】この非線形エコー抑圧処理部２５による処
理は、エコー経路の変動などには影響を受けないため、
ロバストにエコーを抑圧できるという利点がある。しか
し、エコーと同時に送話音声が存在する場合には、それ
らを区別なく抑圧してしまうため、送話音声に歪みや音
の途切れを引き起こすといった問題点がある。
理は、エコー経路の変動などには影響を受けないため、
ロバストにエコーを抑圧できるという利点がある。しか
し、エコーと同時に送話音声が存在する場合には、それ
らを区別なく抑圧してしまうため、送話音声に歪みや音
の途切れを引き起こすといった問題点がある。
【００１４】本発明は、このような問題点に鑑みてなさ
れたもので、以下の点を目的とするものである。 （１）ダブルトーク中のエコー経路の変動に対して、送
話音声を歪ませたりあるいは音切れさせることなく、聴
感上十分にエコーを抑圧して良好な通話品質を確保す
る。 （２）送話音声を歪ませたりあるいは音切れさせること
なく、適応フィルタのみでは完全に消去できない残留エ
コーを聴感上十分に抑圧して良好な通話品質を確保す
る。
れたもので、以下の点を目的とするものである。 （１）ダブルトーク中のエコー経路の変動に対して、送
話音声を歪ませたりあるいは音切れさせることなく、聴
感上十分にエコーを抑圧して良好な通話品質を確保す
る。 （２）送話音声を歪ませたりあるいは音切れさせること
なく、適応フィルタのみでは完全に消去できない残留エ
コーを聴感上十分に抑圧して良好な通話品質を確保す
る。
【００１５】
【課題を解決するための手段】上記目的を達成するため
に、本発明では、エコー抑圧方法に係わる第１の手段と
して、受話信号ｘ（ｋ）が時間変動するエコー経路伝達
関数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ
（ｋ）として送話信号ｓ（ｋ）に重畳された送信信号ｙ
（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエコー
信号ｂ（ｋ）を抑圧する方法において、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルにそれ
ぞれ変換し、送話信号ｓ（ｋ）が含まれない状態におけ
る送信信号ｙ（ｋ）の短時間スペクトル及び前記受話信
号ｘ（ｋ）の短時間スペクトルを用いてエコー経路伝達
関数ｈ（ｋ）の振幅周波数関数を算出し、該振幅周波数
関数並びに送信信号ｙ（ｋ）の短時間スペクトル及び前
記受話信号ｘ（ｋ）の短時間スペクトルを用いて、以降
の時刻におけるエコー信号ｂ（ｋ）の短時間スペクトル
の振幅特性値及び送話信号ｓ（ｋ）の短時間スペクトル
の振幅特性値を推定し、該エコー信号ｂ（ｋ）の短時間
スペクトルの振幅特性値及び送話信号ｓ（ｋ）の短時間
スペクトルの振幅特性値及び送信信号ｙ（ｋ）の短時間
スペクトルの振幅特性値に基づく演算式を用いてエコー
抑圧ゲインＧ（ω）を算出し、該エコー抑圧ゲインＧ
（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗算す
ることによりエコー信号ｂ（ｋ）を抑圧するという手段
を採用する。
に、本発明では、エコー抑圧方法に係わる第１の手段と
して、受話信号ｘ（ｋ）が時間変動するエコー経路伝達
関数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ
（ｋ）として送話信号ｓ（ｋ）に重畳された送信信号ｙ
（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエコー
信号ｂ（ｋ）を抑圧する方法において、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルにそれ
ぞれ変換し、送話信号ｓ（ｋ）が含まれない状態におけ
る送信信号ｙ（ｋ）の短時間スペクトル及び前記受話信
号ｘ（ｋ）の短時間スペクトルを用いてエコー経路伝達
関数ｈ（ｋ）の振幅周波数関数を算出し、該振幅周波数
関数並びに送信信号ｙ（ｋ）の短時間スペクトル及び前
記受話信号ｘ（ｋ）の短時間スペクトルを用いて、以降
の時刻におけるエコー信号ｂ（ｋ）の短時間スペクトル
の振幅特性値及び送話信号ｓ（ｋ）の短時間スペクトル
の振幅特性値を推定し、該エコー信号ｂ（ｋ）の短時間
スペクトルの振幅特性値及び送話信号ｓ（ｋ）の短時間
スペクトルの振幅特性値及び送信信号ｙ（ｋ）の短時間
スペクトルの振幅特性値に基づく演算式を用いてエコー
抑圧ゲインＧ（ω）を算出し、該エコー抑圧ゲインＧ
（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗算す
ることによりエコー信号ｂ（ｋ）を抑圧するという手段
を採用する。
【００１６】エコー抑圧方法に係わる第２の手段とし
て、受話信号ｘ（ｋ）が時間変動するエコー経路伝達関
数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ（ｋ）
として送話信号ｓ（ｋ）に重畳された送信信号ｙ（ｋ）
に対し、該送信信号ｙ（ｋ）に重畳されたエコー信号ｂ
（ｋ）を抑圧する方法において、受話信号ｘ（ｋ）及び
送信信号ｙ（ｋ）を短時間スペクトルにそれぞれ変換
し、ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変
動を検出し、エコー低減信号ｅ（ｋ）と受話信号ｘ
（ｋ）とに基づいてエコー経路伝達関数ｈ（ｋ）の推定
値である疑似エコー経路伝達関数ｈe（ｋ）を算出し、
該疑似エコー経路伝達関数ｈe（ｋ）を用いて送信信号
ｙ（ｋ）に含まれるエコー信号ｂ（ｋ）を抑圧して送話
信号ｓ（ｋ）に残留エコー信号ｂd（ｋ）が含まれた前
記エコー低減信号ｅ（ｋ）を生成し、ダブルトーク中の
エコー経路伝達関数ｈ（ｋ）の変動が検出された場合
に、該変動検出前の疑似エコー経路伝達関数ｈe（ｋ）
から疑似振幅周波数関数｜Ｈe（ω）｜2を算出し、ダブ
ルトーク中のエコー経路伝達関数ｈ（ｋ）の変動が検出
されない場合に、送話信号ｓ（ｋ）が含まれない状態に
おける送信信号ｙ（ｋ）の短時間スペクトル及び前記受
話信号ｘ（ｋ）の短時間スペクトルを用いて残留振幅周
波数関数｜Ｈd（ω）｜2を算出し、択一的に算出された
前記疑似振幅周波数関数｜Ｈe（ω）｜2あるいは残留振
幅周波数関数｜Ｈd（ω）｜2に基づいてエコー信号ｂ
（ｋ）あるいは残留エコー信号ｂd（ｋ）の短時間スペ
クトルの振幅特性値を推定し、該エコー信号ｂ（ｋ）あ
るいは残留エコー信号ｂd（ｋ）の短時間スペクトルの
振幅特性値及び送話信号ｓ（ｋ）の短時間スペクトルの
振幅特性値及び送信信号ｙ（ｋ）の短時間スペクトルの
振幅特性値に基づく演算式を用いてエコー抑圧ゲインＧ
（ω）を算出し、ダブルトーク中のエコー経路伝達関数
ｈ（ｋ）の変動が検出されない場合には、前記エコー抑
圧ゲインＧ（ω）をエコー低減信号ｅ（ｋ）の短時間ス
ペクトルに乗算することにより残留エコー信号ｂd
（ｋ）を抑圧し、ダブルトーク中のエコー経路伝達関数
ｈ（ｋ）の変動が検出された場合には、エコー抑圧ゲイ
ンＧ（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗
算することによりエコー信号ｂ（ｋ）を抑圧するという
手段を採用する。
て、受話信号ｘ（ｋ）が時間変動するエコー経路伝達関
数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ（ｋ）
として送話信号ｓ（ｋ）に重畳された送信信号ｙ（ｋ）
に対し、該送信信号ｙ（ｋ）に重畳されたエコー信号ｂ
（ｋ）を抑圧する方法において、受話信号ｘ（ｋ）及び
送信信号ｙ（ｋ）を短時間スペクトルにそれぞれ変換
し、ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変
動を検出し、エコー低減信号ｅ（ｋ）と受話信号ｘ
（ｋ）とに基づいてエコー経路伝達関数ｈ（ｋ）の推定
値である疑似エコー経路伝達関数ｈe（ｋ）を算出し、
該疑似エコー経路伝達関数ｈe（ｋ）を用いて送信信号
ｙ（ｋ）に含まれるエコー信号ｂ（ｋ）を抑圧して送話
信号ｓ（ｋ）に残留エコー信号ｂd（ｋ）が含まれた前
記エコー低減信号ｅ（ｋ）を生成し、ダブルトーク中の
エコー経路伝達関数ｈ（ｋ）の変動が検出された場合
に、該変動検出前の疑似エコー経路伝達関数ｈe（ｋ）
から疑似振幅周波数関数｜Ｈe（ω）｜2を算出し、ダブ
ルトーク中のエコー経路伝達関数ｈ（ｋ）の変動が検出
されない場合に、送話信号ｓ（ｋ）が含まれない状態に
おける送信信号ｙ（ｋ）の短時間スペクトル及び前記受
話信号ｘ（ｋ）の短時間スペクトルを用いて残留振幅周
波数関数｜Ｈd（ω）｜2を算出し、択一的に算出された
前記疑似振幅周波数関数｜Ｈe（ω）｜2あるいは残留振
幅周波数関数｜Ｈd（ω）｜2に基づいてエコー信号ｂ
（ｋ）あるいは残留エコー信号ｂd（ｋ）の短時間スペ
クトルの振幅特性値を推定し、該エコー信号ｂ（ｋ）あ
るいは残留エコー信号ｂd（ｋ）の短時間スペクトルの
振幅特性値及び送話信号ｓ（ｋ）の短時間スペクトルの
振幅特性値及び送信信号ｙ（ｋ）の短時間スペクトルの
振幅特性値に基づく演算式を用いてエコー抑圧ゲインＧ
（ω）を算出し、ダブルトーク中のエコー経路伝達関数
ｈ（ｋ）の変動が検出されない場合には、前記エコー抑
圧ゲインＧ（ω）をエコー低減信号ｅ（ｋ）の短時間ス
ペクトルに乗算することにより残留エコー信号ｂd
（ｋ）を抑圧し、ダブルトーク中のエコー経路伝達関数
ｈ（ｋ）の変動が検出された場合には、エコー抑圧ゲイ
ンＧ（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗
算することによりエコー信号ｂ（ｋ）を抑圧するという
手段を採用する。
【００１７】エコー抑圧方法に係わる第３の手段とし
て、上記第２の手段において、送信信号ｙ（ｋ）とエコ
ー低減信号ｅ（ｋ）とからエコー信号ｂ（ｋ）の低減量
を求め、受話信号ｘ（ｋ）が存在し、かつ前記低減量が
単調増大している場合に送話信号ｓ（ｋ）が含まれてい
ない状態であると判断するという手段を採用する。
て、上記第２の手段において、送信信号ｙ（ｋ）とエコ
ー低減信号ｅ（ｋ）とからエコー信号ｂ（ｋ）の低減量
を求め、受話信号ｘ（ｋ）が存在し、かつ前記低減量が
単調増大している場合に送話信号ｓ（ｋ）が含まれてい
ない状態であると判断するという手段を採用する。
【００１８】エコー抑圧方法に係わる第４の手段とし
て、上記第１〜第３いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値と１単位時
間前の該エコー信号ｂ（ｋ）の短時間スペクトルの振幅
特性値のうち、何れか大きい方を当該エコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値の推定値とする
という手段を採用する。
て、上記第１〜第３いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値と１単位時
間前の該エコー信号ｂ（ｋ）の短時間スペクトルの振幅
特性値のうち、何れか大きい方を当該エコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値の推定値とする
という手段を採用する。
【００１９】エコー抑圧方法に係わる第５の手段とし
て、上記第１〜第４いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値に聴感補正
を施した値をエコー信号ｂ（ｋ）の短時間スペクトルの
振幅特性値の推定値とするという手段を採用する。
て、上記第１〜第４いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値に聴感補正
を施した値をエコー信号ｂ（ｋ）の短時間スペクトルの
振幅特性値の推定値とするという手段を採用する。
【００２０】エコー抑圧方法に係わる第６の手段とし
て、上記第５の手段において、受話音声によるエコー信
号のマスキングレベルと人間の可聴レベルと伝送遅延量
に基づいて聴感補正ゲインを設定し、該聴感補正ゲイン
をエコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
に乗算して聴感補正するという手段を採用する。
て、上記第５の手段において、受話音声によるエコー信
号のマスキングレベルと人間の可聴レベルと伝送遅延量
に基づいて聴感補正ゲインを設定し、該聴感補正ゲイン
をエコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
に乗算して聴感補正するという手段を採用する。
【００２１】エコー抑圧方法に係わる第７の手段とし
て、上記第１〜第６いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値と１単位時
間前の送話信号ｓ（ｋ）の短時間スペクトルの振幅特性
値と送信信号ｙ（ｋ）の短時間スペクトル並びにスムー
ジングファクタβと信号の正の成分だけを取り出す半波
整流関数Ｐからなる下式（１）に基づいて送話信号ｓ
（ｋ）の短時間スペクトルの振幅特性値の推定値を算出
するという手段を採用する。
て、上記第１〜第６いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値と１単位時
間前の送話信号ｓ（ｋ）の短時間スペクトルの振幅特性
値と送信信号ｙ（ｋ）の短時間スペクトル並びにスムー
ジングファクタβと信号の正の成分だけを取り出す半波
整流関数Ｐからなる下式（１）に基づいて送話信号ｓ
（ｋ）の短時間スペクトルの振幅特性値の推定値を算出
するという手段を採用する。
【００２２】
【数１９】
【００２３】エコー抑圧方法に係わる第８の手段とし
て、上記第７の手段において、スージングファクタβを
０．９〜０．９９９の範囲内の値とするという手段を採
用する。
て、上記第７の手段において、スージングファクタβを
０．９〜０．９９９の範囲内の値とするという手段を採
用する。
【００２４】エコー抑圧方法に係わる第９の手段とし
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（２）の値ＧPE（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用するという手段を
採用する。
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（２）の値ＧPE（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用するという手段を
採用する。
【００２５】
【数２０】
【００２６】エコー抑圧方法に係わる第１０の手段とし
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（３）の値ＧW（ω）を、前
記エコー抑圧ゲインＧ（ω）に適用するという手段を採
用する。
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（３）の値ＧW（ω）を、前
記エコー抑圧ゲインＧ（ω）に適用するという手段を採
用する。
【００２７】
【数２１】
【００２８】エコー抑圧方法に係わる第１１の手段とし
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（４）の値ＧML（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用するという手段を
採用する。
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（４）の値ＧML（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用するという手段を
採用する。
【００２９】
【数２２】
【００３０】エコー抑圧方法に係わる第１２の手段とし
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2と送信信号パワー｜Ｙt（ω）｜2からなる式
（５）に基づいて中間変数υ（ω）を算出し、該中間変
数υ（ω）とエコー信号パワー｜Ｅte（ω）｜2と送信
信号パワー｜Ｙt（ω）｜2並びに中間変数υ（ω）に対
する０次及び１次の修正ベッセル関数Ｉ0，Ｉ1からなる
演算式（６）の値ＧEM（ω）を、前記エコー抑圧ゲイン
Ｇ（ω）に適用するという手段を採用する。
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2と送信信号パワー｜Ｙt（ω）｜2からなる式
（５）に基づいて中間変数υ（ω）を算出し、該中間変
数υ（ω）とエコー信号パワー｜Ｅte（ω）｜2と送信
信号パワー｜Ｙt（ω）｜2並びに中間変数υ（ω）に対
する０次及び１次の修正ベッセル関数Ｉ0，Ｉ1からなる
演算式（６）の値ＧEM（ω）を、前記エコー抑圧ゲイン
Ｇ（ω）に適用するという手段を採用する。
【００３１】
【数２３】
【００３２】
【数２４】
【００３３】一方、本発明では、エコー抑圧装置に係わ
る第１の手段として、受話信号ｘ（ｋ）が時間変動する
エコー経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、エ
コー信号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳された
送信信号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳さ
れたエコー信号ｂ（ｋ）を抑圧するエコー抑圧装置にお
いて、受話信号ｘ（ｋ）を短時間スペクトルに変換する
受話信号変換手段と、送信信号ｙ（ｋ）を短時間スペク
トルに変換する送信信号変換手段と、受話信号ｘ（ｋ）
と送信信号ｙ（ｋ）とに基づいて送話信号ｓ（ｋ）が含
まれていない状態を検出する片側発話状態検出部と、送
話信号ｓ（ｋ）が含まれない状態における送信信号ｙ
（ｋ）の短時間スペクトル及び前記受話信号ｘ（ｋ）の
短時間スペクトルを用いてエコー経路伝達関数ｈ（ｋ）
の振幅周波数関数を算出するエコー経路伝達関数算出手
段と、該振幅周波数関数並びに送信信号ｙ（ｋ）の短時
間スペクトル及び前記受話信号ｘ（ｋ）の短時間スペク
トルを用いて、以降の時刻におけるエコー信号ｂ（ｋ）
の短時間スペクトルの振幅特性値を推定するエコー信号
振幅周波数特性推定手段と、該振幅周波数関数並びに送
信信号ｙ（ｋ）の短時間スペクトル及び前記受話信号ｘ
（ｋ）の短時間スペクトルを用いて、以降の時刻におけ
る送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値を
推定する送話信号振幅周波数特性推定手段と、該エコー
信号ｂ（ｋ）の短時間スペクトルの振幅特性値及び送話
信号ｓ（ｋ）の短時間スペクトルの振幅特性値及び送信
信号ｙ（ｋ）の短時間スペクトルの振幅特性値に基づく
演算式を用いてエコー抑圧ゲインＧ（ω）を算出するエ
コー抑圧ゲイン決定部と、該エコー抑圧ゲインＧ（ω）
を送信信号ｙ（ｋ）の短時間スペクトルに乗算するエコ
ー抑圧処理部と、該エコー抑圧処理部の出力を逆変換し
て伝送路に出力する送話信号逆変換手段とを具備する手
段を採用する。
る第１の手段として、受話信号ｘ（ｋ）が時間変動する
エコー経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、エ
コー信号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳された
送信信号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳さ
れたエコー信号ｂ（ｋ）を抑圧するエコー抑圧装置にお
いて、受話信号ｘ（ｋ）を短時間スペクトルに変換する
受話信号変換手段と、送信信号ｙ（ｋ）を短時間スペク
トルに変換する送信信号変換手段と、受話信号ｘ（ｋ）
と送信信号ｙ（ｋ）とに基づいて送話信号ｓ（ｋ）が含
まれていない状態を検出する片側発話状態検出部と、送
話信号ｓ（ｋ）が含まれない状態における送信信号ｙ
（ｋ）の短時間スペクトル及び前記受話信号ｘ（ｋ）の
短時間スペクトルを用いてエコー経路伝達関数ｈ（ｋ）
の振幅周波数関数を算出するエコー経路伝達関数算出手
段と、該振幅周波数関数並びに送信信号ｙ（ｋ）の短時
間スペクトル及び前記受話信号ｘ（ｋ）の短時間スペク
トルを用いて、以降の時刻におけるエコー信号ｂ（ｋ）
の短時間スペクトルの振幅特性値を推定するエコー信号
振幅周波数特性推定手段と、該振幅周波数関数並びに送
信信号ｙ（ｋ）の短時間スペクトル及び前記受話信号ｘ
（ｋ）の短時間スペクトルを用いて、以降の時刻におけ
る送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値を
推定する送話信号振幅周波数特性推定手段と、該エコー
信号ｂ（ｋ）の短時間スペクトルの振幅特性値及び送話
信号ｓ（ｋ）の短時間スペクトルの振幅特性値及び送信
信号ｙ（ｋ）の短時間スペクトルの振幅特性値に基づく
演算式を用いてエコー抑圧ゲインＧ（ω）を算出するエ
コー抑圧ゲイン決定部と、該エコー抑圧ゲインＧ（ω）
を送信信号ｙ（ｋ）の短時間スペクトルに乗算するエコ
ー抑圧処理部と、該エコー抑圧処理部の出力を逆変換し
て伝送路に出力する送話信号逆変換手段とを具備する手
段を採用する。
【００３４】エコー抑圧装置に係わる第２の手段とし
て、受話信号ｘ（ｋ）が時間変動するエコー経路伝達関
数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ（ｋ）
として送話信号ｓ（ｋ）に重畳された送信信号ｙ（ｋ）
に対し、該送信信号ｙ（ｋ）に重畳されたエコー信号ｂ
（ｋ）を抑圧するエコー抑圧装置であって、受話信号ｘ
（ｋ）を短時間スペクトルに変換する受話信号変換手段
と、送信信号ｙ（ｋ）を短時間スペクトルに変換する送
信信号変換手段と、受話信号ｘ（ｋ）と送信信号ｙ
（ｋ）とに基づいて送話信号ｓ（ｋ）が含まれていない
状態を検出する片側発話状態検出部と、ダブルトーク中
のエコー経路伝達関数ｈ（ｋ）の変動を検出するエコー
経路変動検出部と、適応フィルタを用いてエコー経路伝
達関数ｈ（ｋ）の推定値である疑似エコー経路伝達関数
ｈe（ｋ）を生成し、該疑似エコー経路伝達関数ｈe
（ｋ）を用いて送信信号ｙ（ｋ）に含まれるエコー信号
ｂ（ｋ）を抑圧して送話信号ｓ（ｋ）に残留エコー信号
ｂd（ｋ）が含まれたエコー低減信号ｅ（ｋ）を出力す
ると共に、前記エコー経路変動検出部によってダブルト
ーク中のエコー経路伝達関数ｈ（ｋ）の変動が検出され
たときに、該変動検出前の疑似エコー経路伝達関数ｈe
（ｋ）を出力すると共に送信信号ｙ（ｋ）をそのまま出
力するエコーキャンセラと、 前記エコー経路変動検出
部によってダブルトーク中のエコー経路伝達関数ｈ
（ｋ）の変動が検出されない場合に、送話信号ｓ（ｋ）
が含まれない場合の前記エコー低減信号ｅ（ｋ）の短時
間スペクトルと受話信号ｘ（ｋ）の短時間スペクトルと
に基づいて残留振幅周波数関数｜Ｈd（ω）｜2を算出す
るエコー経路伝達関数算出手段と、前記エコー経路変動
検出部によってダブルトーク中のエコー経路伝達関数ｈ
（ｋ）の変動が検出された場合に、エコーキャンセラか
ら出力された疑似エコー経路伝達関数ｈe（ｋ）から疑
似振幅周波数関数｜Ｈe（ω）｜2を算出する第２のエコ
ー経路伝達関数算出手段と、エコー経路伝達関数算出手
段あるいは第２のエコー経路伝達関数算出手段から択一
的に入力される前記疑似振幅周波数関数｜Ｈe（ω）｜2
あるいは残留振幅周波数関数｜Ｈd（ω）｜2に基づいて
エコー信号ｂ（ｋ）あるいは残留エコー信号ｂd（ｋ）
の短時間スペクトルの振幅周波数特性を推定するエコー
信号振幅周波数特性推定手段と、該振幅周波数特性並び
にエコー低減信号ｅ（ｋ）あるいは送信信号ｙ（ｋ）の
短時間スペクトル及び前記受話信号ｘ（ｋ）の短時間ス
ペクトルを用いて、以降の時刻における残留エコー信号
ｂd（ｋ）あるいはエコー信号ｂ（ｋ）の短時間スペク
トルの振幅周波数特性を推定するエコー信号振幅周波数
特性推定手段と、該振幅周波数特性並びにエコー低減信
号ｅ（ｋ）あるいは送信信号ｙ（ｋ）の短時間スペクト
ル及び前記受話信号ｘ（ｋ）の短時間スペクトルを用い
て、以降の時刻における送話信号ｓ（ｋ）の短時間スペ
クトルの振幅周波数特性を推定する送話信号振幅周波数
特性推定手段と、前記残留エコー信号ｂd（ｋ）あるい
はエコー信号ｂ（ｋ）の短時間スペクトルの振幅周波数
特性及び送話信号ｓ（ｋ）の短時間スペクトルの振幅周
波数特性及びエコー低減信号ｅ（ｋ）あるいは送信信号
ｙ（ｋ）の短時間スペクトルの振幅周波数特性に基づく
演算式を用いてエコー抑圧ゲインＧ（ω）を算出するエ
コー抑圧ゲイン決定部と、該エコー抑圧ゲインＧ（ω）
をエコー低減信号ｅ（ｋ）あるいは送信信号ｙ（ｋ）の
短時間スペクトルに乗算するエコー抑圧処理部と、該エ
コー抑圧処理部の出力を逆変換して伝送路に出力する送
話信号逆変換手段とを具備する手段を採用する。
て、受話信号ｘ（ｋ）が時間変動するエコー経路伝達関
数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ（ｋ）
として送話信号ｓ（ｋ）に重畳された送信信号ｙ（ｋ）
に対し、該送信信号ｙ（ｋ）に重畳されたエコー信号ｂ
（ｋ）を抑圧するエコー抑圧装置であって、受話信号ｘ
（ｋ）を短時間スペクトルに変換する受話信号変換手段
と、送信信号ｙ（ｋ）を短時間スペクトルに変換する送
信信号変換手段と、受話信号ｘ（ｋ）と送信信号ｙ
（ｋ）とに基づいて送話信号ｓ（ｋ）が含まれていない
状態を検出する片側発話状態検出部と、ダブルトーク中
のエコー経路伝達関数ｈ（ｋ）の変動を検出するエコー
経路変動検出部と、適応フィルタを用いてエコー経路伝
達関数ｈ（ｋ）の推定値である疑似エコー経路伝達関数
ｈe（ｋ）を生成し、該疑似エコー経路伝達関数ｈe
（ｋ）を用いて送信信号ｙ（ｋ）に含まれるエコー信号
ｂ（ｋ）を抑圧して送話信号ｓ（ｋ）に残留エコー信号
ｂd（ｋ）が含まれたエコー低減信号ｅ（ｋ）を出力す
ると共に、前記エコー経路変動検出部によってダブルト
ーク中のエコー経路伝達関数ｈ（ｋ）の変動が検出され
たときに、該変動検出前の疑似エコー経路伝達関数ｈe
（ｋ）を出力すると共に送信信号ｙ（ｋ）をそのまま出
力するエコーキャンセラと、 前記エコー経路変動検出
部によってダブルトーク中のエコー経路伝達関数ｈ
（ｋ）の変動が検出されない場合に、送話信号ｓ（ｋ）
が含まれない場合の前記エコー低減信号ｅ（ｋ）の短時
間スペクトルと受話信号ｘ（ｋ）の短時間スペクトルと
に基づいて残留振幅周波数関数｜Ｈd（ω）｜2を算出す
るエコー経路伝達関数算出手段と、前記エコー経路変動
検出部によってダブルトーク中のエコー経路伝達関数ｈ
（ｋ）の変動が検出された場合に、エコーキャンセラか
ら出力された疑似エコー経路伝達関数ｈe（ｋ）から疑
似振幅周波数関数｜Ｈe（ω）｜2を算出する第２のエコ
ー経路伝達関数算出手段と、エコー経路伝達関数算出手
段あるいは第２のエコー経路伝達関数算出手段から択一
的に入力される前記疑似振幅周波数関数｜Ｈe（ω）｜2
あるいは残留振幅周波数関数｜Ｈd（ω）｜2に基づいて
エコー信号ｂ（ｋ）あるいは残留エコー信号ｂd（ｋ）
の短時間スペクトルの振幅周波数特性を推定するエコー
信号振幅周波数特性推定手段と、該振幅周波数特性並び
にエコー低減信号ｅ（ｋ）あるいは送信信号ｙ（ｋ）の
短時間スペクトル及び前記受話信号ｘ（ｋ）の短時間ス
ペクトルを用いて、以降の時刻における残留エコー信号
ｂd（ｋ）あるいはエコー信号ｂ（ｋ）の短時間スペク
トルの振幅周波数特性を推定するエコー信号振幅周波数
特性推定手段と、該振幅周波数特性並びにエコー低減信
号ｅ（ｋ）あるいは送信信号ｙ（ｋ）の短時間スペクト
ル及び前記受話信号ｘ（ｋ）の短時間スペクトルを用い
て、以降の時刻における送話信号ｓ（ｋ）の短時間スペ
クトルの振幅周波数特性を推定する送話信号振幅周波数
特性推定手段と、前記残留エコー信号ｂd（ｋ）あるい
はエコー信号ｂ（ｋ）の短時間スペクトルの振幅周波数
特性及び送話信号ｓ（ｋ）の短時間スペクトルの振幅周
波数特性及びエコー低減信号ｅ（ｋ）あるいは送信信号
ｙ（ｋ）の短時間スペクトルの振幅周波数特性に基づく
演算式を用いてエコー抑圧ゲインＧ（ω）を算出するエ
コー抑圧ゲイン決定部と、該エコー抑圧ゲインＧ（ω）
をエコー低減信号ｅ（ｋ）あるいは送信信号ｙ（ｋ）の
短時間スペクトルに乗算するエコー抑圧処理部と、該エ
コー抑圧処理部の出力を逆変換して伝送路に出力する送
話信号逆変換手段とを具備する手段を採用する。
【００３５】エコー抑圧装置に係わる第３の手段とし
て、上記第２の手段において、送信信号ｙ（ｋ）とエコ
ー低減信号ｅ（ｋ）とからエコー信号ｂ（ｋ）の低減量
を求め、受話信号ｘ（ｋ）が存在し、かつ前記低減量が
単調増大している場合に送話信号ｓ（ｋ）が含まれてい
ない状態であると判断するように片側発話状態検出部を
構成するという手段を採用する。
て、上記第２の手段において、送信信号ｙ（ｋ）とエコ
ー低減信号ｅ（ｋ）とからエコー信号ｂ（ｋ）の低減量
を求め、受話信号ｘ（ｋ）が存在し、かつ前記低減量が
単調増大している場合に送話信号ｓ（ｋ）が含まれてい
ない状態であると判断するように片側発話状態検出部を
構成するという手段を採用する。
【００３６】エコー抑圧装置に係わる第４の手段とし
て、上記第１〜第３いずれかの手段において、振幅周波
数関数並びに送信信号ｙ（ｋ）の短時間スペクトル及び
前記受話信号ｘ（ｋ）の短時間スペクトルを用いて、以
降の時刻におけるエコー信号ｂ（ｋ）の短時間スペクト
ルの振幅特性値を推定するエコー信号パワー計算部と、
この振幅特性値を１単位時間遅延させて出力する第１の
遅延演算子と、エコー信号パワー計算部と第１の遅延演
算子の出力のうち、大きい方を当該エコー信号ｂ（ｋ）
の短時間スペクトルの振幅特性値の推定値とするエコー
信号パワー決定部とからエコー信号振幅周波数特性推定
手段を構成するという手段を採用する。
て、上記第１〜第３いずれかの手段において、振幅周波
数関数並びに送信信号ｙ（ｋ）の短時間スペクトル及び
前記受話信号ｘ（ｋ）の短時間スペクトルを用いて、以
降の時刻におけるエコー信号ｂ（ｋ）の短時間スペクト
ルの振幅特性値を推定するエコー信号パワー計算部と、
この振幅特性値を１単位時間遅延させて出力する第１の
遅延演算子と、エコー信号パワー計算部と第１の遅延演
算子の出力のうち、大きい方を当該エコー信号ｂ（ｋ）
の短時間スペクトルの振幅特性値の推定値とするエコー
信号パワー決定部とからエコー信号振幅周波数特性推定
手段を構成するという手段を採用する。
【００３７】エコー抑圧装置に係わる第５の手段とし
て、上記第１〜第４いずれかの手段において、エコー信
号振幅周波数特性推定手段に、エコー信号ｂ（ｋ）の短
時間スペクトルの振幅特性値に聴感補正を施すエコー信
号聴感補正手段を備えるという手段を採用する。
て、上記第１〜第４いずれかの手段において、エコー信
号振幅周波数特性推定手段に、エコー信号ｂ（ｋ）の短
時間スペクトルの振幅特性値に聴感補正を施すエコー信
号聴感補正手段を備えるという手段を採用する。
【００３８】エコー抑圧装置に係わる第６の手段とし
て、上記第５の手段において、受話音声によるエコー信
号のマスキングレベルと人間の可聴レベルと伝送遅延量
に基づいて設定された聴感補正ゲインを記憶する聴感補
正ゲイン記憶部と、前記聴感補正ゲインをエコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値に乗算する聴感
補正部とからエコー信号聴感補正手段を構成するという
手段を採用する。
て、上記第５の手段において、受話音声によるエコー信
号のマスキングレベルと人間の可聴レベルと伝送遅延量
に基づいて設定された聴感補正ゲインを記憶する聴感補
正ゲイン記憶部と、前記聴感補正ゲインをエコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値に乗算する聴感
補正部とからエコー信号聴感補正手段を構成するという
手段を採用する。
【００３９】エコー抑圧装置に係わる第７の手段とし
て、上記第１〜第６いずれかの手段において、スムージ
ングファクタβを記憶するスムージングファクタ記憶部
と、送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値
を１単位時間遅延させる第２の遅延演算子と、該第２の
遅延演算子の出力とエコー信号ｂ（ｋ）の短時間スペク
トルの振幅特性値と送信信号ｙ（ｋ）の短時間スペクト
ルとスムージングファクタβと信号の正の成分だけを取
り出す半波整流関数Ｐからなる下式（１）に基づいて送
話信号ｓ（ｋ）の短時間スペクトルの振幅特性値の推定
値を算出する送信信号パワー計算部とから送話信号振幅
周波数特性推定手段を構成するという手段を採用する。
て、上記第１〜第６いずれかの手段において、スムージ
ングファクタβを記憶するスムージングファクタ記憶部
と、送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値
を１単位時間遅延させる第２の遅延演算子と、該第２の
遅延演算子の出力とエコー信号ｂ（ｋ）の短時間スペク
トルの振幅特性値と送信信号ｙ（ｋ）の短時間スペクト
ルとスムージングファクタβと信号の正の成分だけを取
り出す半波整流関数Ｐからなる下式（１）に基づいて送
話信号ｓ（ｋ）の短時間スペクトルの振幅特性値の推定
値を算出する送信信号パワー計算部とから送話信号振幅
周波数特性推定手段を構成するという手段を採用する。
【００４０】
【数２５】
【００４１】エコー抑圧装置に係わる第８の手段とし
て、上記第７の手段において、スージングファクタβを
０．９〜０．９９９の範囲内の値とするという手段を採
用する。
て、上記第７の手段において、スージングファクタβを
０．９〜０．９９９の範囲内の値とするという手段を採
用する。
【００４２】エコー抑圧装置に係わる第９の手段とし
て、上記第１〜第８いずれかの手段において、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（２）に基づいて値ＧPE（ω）を算出するようにエコ
ー抑圧ゲイン決定部を構成する手段を採用する。
て、上記第１〜第８いずれかの手段において、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（２）に基づいて値ＧPE（ω）を算出するようにエコ
ー抑圧ゲイン決定部を構成する手段を採用する。
【００４３】
【数２６】
【００４４】エコー抑圧装置に係わる第１０の手段とし
て、上記第１〜第８いずれかの手段において、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（３）に基づいて値ＧW（ω）を算出するようにエコ
ー抑圧ゲイン決定部を構成すれるという手段を採用す
る。
て、上記第１〜第８いずれかの手段において、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（３）に基づいて値ＧW（ω）を算出するようにエコ
ー抑圧ゲイン決定部を構成すれるという手段を採用す
る。
【００４５】
【数２７】
【００４６】エコー抑圧装置に係わる第１１の手段とし
て、上記第１〜第８いずれかの手段において、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（４）に基づいて値ＧML（ω）を算出するようにエコ
ー抑圧ゲイン決定部を構成するという手段を採用する。
て、上記第１〜第８いずれかの手段において、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（４）に基づいて値ＧML（ω）を算出するようにエコ
ー抑圧ゲイン決定部を構成するという手段を採用する。
【００４７】
【数２８】
【００４８】エコー抑圧装置に係わる第１２の手段とし
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2と送信信号パワー｜Ｙt（ω）｜2からなる式
（５）に基づいて中間変数υ（ω）を算出し、該中間変
数υ（ω）とエコー信号パワー｜Ｅte（ω）｜2と送信
信号パワー｜Ｙt（ω）｜2並びに中間変数υ（ω）に対
する０次及び１次の修正ベッセル関数Ｉ0，Ｉ1からなる
演算式（６）に基づいて値ＧEM（ω）をエコー抑圧ゲイ
ンＧ（ω）として算出するようにエコー抑圧ゲイン決定
部を構成するという手段を採用する。
て、上記第１〜第８いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2と送信信号パワー｜Ｙt（ω）｜2からなる式
（５）に基づいて中間変数υ（ω）を算出し、該中間変
数υ（ω）とエコー信号パワー｜Ｅte（ω）｜2と送信
信号パワー｜Ｙt（ω）｜2並びに中間変数υ（ω）に対
する０次及び１次の修正ベッセル関数Ｉ0，Ｉ1からなる
演算式（６）に基づいて値ＧEM（ω）をエコー抑圧ゲイ
ンＧ（ω）として算出するようにエコー抑圧ゲイン決定
部を構成するという手段を採用する。
【００４９】
【数２９】
【００５０】
【数３０】
【００５１】さらに、本発明では、エコー抑圧プログラ
ムが記憶されたコンピュータに読取り可能な記憶媒体係
わる第１の手段として、受話信号ｘ（ｋ）が時間変動す
るエコー経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、
エコー信号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳され
た送信信号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳
されたエコー信号ｂ（ｋ）を抑圧するエコー抑圧プログ
ラムが記憶されたコンピュータに読取り可能な記憶媒体
において、受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）を短
時間スペクトルにそれぞれ変換する処理と、送話信号ｓ
（ｋ）が含まれない状態における送信信号ｙ（ｋ）の短
時間スペクトル及び前記受話信号ｘ（ｋ）の短時間スペ
クトルを用いてエコー経路伝達関数ｈ（ｋ）の振幅周波
数関数を算出する処理と、該振幅周波数関数並びに送信
信号ｙ（ｋ）の短時間スペクトル及び前記受話信号ｘ
（ｋ）の短時間スペクトルを用いて、以降の時刻におけ
るエコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
及び送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値
を推定する処理と、該エコー信号ｂ（ｋ）の短時間スペ
クトルの振幅特性値及び送話信号ｓ（ｋ）の短時間スペ
クトルの振幅特性値及び送信信号ｙ（ｋ）の短時間スペ
クトルの振幅特性値に基づく演算式を用いてエコー抑圧
ゲインＧ（ω）を算出する処理と、該エコー抑圧ゲイン
Ｇ（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗算
することによりエコー信号ｂ（ｋ）を抑圧するする処理
とからなるエコー抑圧プログラムを記憶するという手段
を採用する。
ムが記憶されたコンピュータに読取り可能な記憶媒体係
わる第１の手段として、受話信号ｘ（ｋ）が時間変動す
るエコー経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、
エコー信号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳され
た送信信号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳
されたエコー信号ｂ（ｋ）を抑圧するエコー抑圧プログ
ラムが記憶されたコンピュータに読取り可能な記憶媒体
において、受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）を短
時間スペクトルにそれぞれ変換する処理と、送話信号ｓ
（ｋ）が含まれない状態における送信信号ｙ（ｋ）の短
時間スペクトル及び前記受話信号ｘ（ｋ）の短時間スペ
クトルを用いてエコー経路伝達関数ｈ（ｋ）の振幅周波
数関数を算出する処理と、該振幅周波数関数並びに送信
信号ｙ（ｋ）の短時間スペクトル及び前記受話信号ｘ
（ｋ）の短時間スペクトルを用いて、以降の時刻におけ
るエコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
及び送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値
を推定する処理と、該エコー信号ｂ（ｋ）の短時間スペ
クトルの振幅特性値及び送話信号ｓ（ｋ）の短時間スペ
クトルの振幅特性値及び送信信号ｙ（ｋ）の短時間スペ
クトルの振幅特性値に基づく演算式を用いてエコー抑圧
ゲインＧ（ω）を算出する処理と、該エコー抑圧ゲイン
Ｇ（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗算
することによりエコー信号ｂ（ｋ）を抑圧するする処理
とからなるエコー抑圧プログラムを記憶するという手段
を採用する。
【００５２】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第２の手段とし
て、上記第１の手段において、エコー信号ｂ（ｋ）の短
時間スペクトルの振幅特性値と１単位時間前の該エコー
信号ｂ（ｋ）の短時間スペクトルの振幅特性値のうち、
何れか大きい方を当該エコー信号ｂ（ｋ）の短時間スペ
クトルの振幅特性値の推定値とする処理を含むエコー抑
圧プログラムを記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第２の手段とし
て、上記第１の手段において、エコー信号ｂ（ｋ）の短
時間スペクトルの振幅特性値と１単位時間前の該エコー
信号ｂ（ｋ）の短時間スペクトルの振幅特性値のうち、
何れか大きい方を当該エコー信号ｂ（ｋ）の短時間スペ
クトルの振幅特性値の推定値とする処理を含むエコー抑
圧プログラムを記憶するという手段を採用する。
【００５３】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第３の手段とし
て、上記第１または第２の手段において、エコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値に聴感補正を施
した値をエコー信号ｂ（ｋ）の短時間スペクトルの振幅
特性値の推定値とする処理を含むエコー抑圧プログラム
を記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第３の手段とし
て、上記第１または第２の手段において、エコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値に聴感補正を施
した値をエコー信号ｂ（ｋ）の短時間スペクトルの振幅
特性値の推定値とする処理を含むエコー抑圧プログラム
を記憶するという手段を採用する。
【００５４】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第４の手段とし
て、上記第３の手段において、受話音声によるエコー信
号のマスキングレベルと人間の可聴レベルと伝送遅延量
に基づいて聴感補正ゲインを設定し、該聴感補正ゲイン
をエコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
に乗算して聴感補正する処理を含むエコー抑圧プログラ
ムを記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第４の手段とし
て、上記第３の手段において、受話音声によるエコー信
号のマスキングレベルと人間の可聴レベルと伝送遅延量
に基づいて聴感補正ゲインを設定し、該聴感補正ゲイン
をエコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
に乗算して聴感補正する処理を含むエコー抑圧プログラ
ムを記憶するという手段を採用する。
【００５５】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第５の手段とし
て、上記第１〜第４いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値と１単位時
間前の送話信号ｓ（ｋ）の短時間スペクトルの振幅特性
値と送信信号ｙ（ｋ）の短時間スペクトル並びにスムー
ジングファクタβと信号の正の成分だけを取り出す半波
整流関数Ｐからなる下式（１）に基づいて送話信号ｓ
（ｋ）の短時間スペクトルの振幅特性値の推定値を算出
する処理を含むエコー抑圧プログラムを記憶するという
手段を採用する。
ュータに読取り可能な記憶媒体係わる第５の手段とし
て、上記第１〜第４いずれかの手段において、エコー信
号ｂ（ｋ）の短時間スペクトルの振幅特性値と１単位時
間前の送話信号ｓ（ｋ）の短時間スペクトルの振幅特性
値と送信信号ｙ（ｋ）の短時間スペクトル並びにスムー
ジングファクタβと信号の正の成分だけを取り出す半波
整流関数Ｐからなる下式（１）に基づいて送話信号ｓ
（ｋ）の短時間スペクトルの振幅特性値の推定値を算出
する処理を含むエコー抑圧プログラムを記憶するという
手段を採用する。
【００５６】
【数３１】
【００５７】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第６の手段とし
て、上記第５の手段において、スムージングファクタβ
を０．９〜０．９９９の範囲内の値とする処理を含むエ
コー抑圧プログラムを記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第６の手段とし
て、上記第５の手段において、スムージングファクタβ
を０．９〜０．９９９の範囲内の値とする処理を含むエ
コー抑圧プログラムを記憶するという手段を採用する。
【００５８】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第７の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（２）の値ＧPE（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用する処理を含むエ
コー抑圧プログラムを記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第７の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（２）の値ＧPE（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用する処理を含むエ
コー抑圧プログラムを記憶するという手段を採用する。
【００５９】
【数３２】
【００６０】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第８の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（３）の値ＧW（ω）を、前
記エコー抑圧ゲインＧ（ω）に適用する処理を含むエコ
ー抑圧プログラムを記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第８の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（３）の値ＧW（ω）を、前
記エコー抑圧ゲインＧ（ω）に適用する処理を含むエコ
ー抑圧プログラムを記憶するという手段を採用する。
【００６１】
【数３３】
【００６２】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第９の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（４）の値ＧML（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用する処理を含むエ
コー抑圧プログラムを記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第９の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2とからなる演算式（４）の値ＧML（ω）を、
前記エコー抑圧ゲインＧ（ω）に適用する処理を含むエ
コー抑圧プログラムを記憶するという手段を採用する。
【００６３】
【数３４】
【００６４】エコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体係わる第１０の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2と送信信号パワー｜Ｙt（ω）｜2からなる式
（５）に基づいて中間変数υ（ω）を算出し、該中間変
数υ（ω）とエコー信号パワー｜Ｅte（ω）｜2と送信
信号パワー｜Ｙt（ω）｜2並びに中間変数υ（ω）に対
する０次及び１次の修正ベッセル関数Ｉ0，Ｉ1からなる
演算式（６）の値ＧEM（ω）を、前記エコー抑圧ゲイン
Ｇ（ω）に適用する処理を含むエコー抑圧プログラムを
記憶するという手段を採用する。
ュータに読取り可能な記憶媒体係わる第１０の手段とし
て、上記第１〜第６いずれかの手段において、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2と送信信号パワー｜Ｙt（ω）｜2からなる式
（５）に基づいて中間変数υ（ω）を算出し、該中間変
数υ（ω）とエコー信号パワー｜Ｅte（ω）｜2と送信
信号パワー｜Ｙt（ω）｜2並びに中間変数υ（ω）に対
する０次及び１次の修正ベッセル関数Ｉ0，Ｉ1からなる
演算式（６）の値ＧEM（ω）を、前記エコー抑圧ゲイン
Ｇ（ω）に適用する処理を含むエコー抑圧プログラムを
記憶するという手段を採用する。
【００６５】
【数３５】
【００６６】
【数３６】
【００６７】
【発明の実施の形態】以下、図１〜図６を参照して、本
発明に係わるエコー抑圧方法及び装置並びにエコー抑圧
プログラムが記憶されたコンピュータに読取り可能な記
憶媒体の一実施形態について説明する。なお、以下の説
明では、上記図８及び図９において既に説明した構成要
素については、同一符号を付してその説明を省略する。
発明に係わるエコー抑圧方法及び装置並びにエコー抑圧
プログラムが記憶されたコンピュータに読取り可能な記
憶媒体の一実施形態について説明する。なお、以下の説
明では、上記図８及び図９において既に説明した構成要
素については、同一符号を付してその説明を省略する。
【００６８】〔第１実施形態〕まず、図１〜図３を参照
して、本発明の第１実施形態について説明する。図１
は、本実施形態の機能構成を示すブロック図である。こ
の図において、符号３１は高速フーリエ変換部（受話信
号変換手段）、３２は高速フーリエ変換部（送信信号変
換手段）、３３は高速逆フーリエ変換部（送話信号逆変
換手段）、３４は片側発話状態検出部、３５はエコー経
路結合量計算部、３６はエコー信号パワー計算部、３
７，４３は遅延演算子、３８はエコー信号パワー決定
部、３９は聴感補正ゲイン記憶部、４０は聴感補正部、
４１はスムージングファクタ記憶部、４２は送話信号パ
ワー計算部、４４はエコー抑圧ゲイン決定部、４５はエ
コー抑圧処理部、４６はエコー経路結合量計算用受話信
号選択部、４７はエコー経路結合量計算用送信信号選択
部である。なお、この図では、受話スピーカ４と送話用
マイクロホン３との問のエコー経路のエコー経路伝達特
性（インパルス応答）をｈ（ｋ）としている。
して、本発明の第１実施形態について説明する。図１
は、本実施形態の機能構成を示すブロック図である。こ
の図において、符号３１は高速フーリエ変換部（受話信
号変換手段）、３２は高速フーリエ変換部（送信信号変
換手段）、３３は高速逆フーリエ変換部（送話信号逆変
換手段）、３４は片側発話状態検出部、３５はエコー経
路結合量計算部、３６はエコー信号パワー計算部、３
７，４３は遅延演算子、３８はエコー信号パワー決定
部、３９は聴感補正ゲイン記憶部、４０は聴感補正部、
４１はスムージングファクタ記憶部、４２は送話信号パ
ワー計算部、４４はエコー抑圧ゲイン決定部、４５はエ
コー抑圧処理部、４６はエコー経路結合量計算用受話信
号選択部、４７はエコー経路結合量計算用送信信号選択
部である。なお、この図では、受話スピーカ４と送話用
マイクロホン３との問のエコー経路のエコー経路伝達特
性（インパルス応答）をｈ（ｋ）としている。
【００６９】ここで、上記エコー経路結合量計算部３５
とエコー経路結合量計算用受話信号選択部４６とエコー
経路結合量計算用送信信号選択部４７とは、エコー経路
伝達関数算出手段を構成するものである。エコー信号パ
ワー計算部３６と遅延演算子３７とエコー信号パワー決
定部３８と聴感補正ゲイン記憶部３９と聴感補正部４０
とは、エコー信号振幅周波数特性推定手段を構成するも
のである。聴感補正ゲイン記憶部３９と聴感補正部４０
とは、エコー信号聴感補正手段を構成するものである。
また、スムージングファクタ記憶部４１と送話信号パワ
ー計算部４２と遅延演算子４３とは、送話信号振幅周波
数特性推定手段を構成するものである。
とエコー経路結合量計算用受話信号選択部４６とエコー
経路結合量計算用送信信号選択部４７とは、エコー経路
伝達関数算出手段を構成するものである。エコー信号パ
ワー計算部３６と遅延演算子３７とエコー信号パワー決
定部３８と聴感補正ゲイン記憶部３９と聴感補正部４０
とは、エコー信号振幅周波数特性推定手段を構成するも
のである。聴感補正ゲイン記憶部３９と聴感補正部４０
とは、エコー信号聴感補正手段を構成するものである。
また、スムージングファクタ記憶部４１と送話信号パワ
ー計算部４２と遅延演算子４３とは、送話信号振幅周波
数特性推定手段を構成するものである。
【００７０】高速フーリエ変換部３１は、受話信号ｘ
（ｋ）を高速フーリエ変換し、短時間スペクトルＸ
t（ω）（ｔ：時間を示す添字）としてエコー信号パワ
ー計算部３６及びエコー経路結合量計算用受話信号選択
部４６に出力する。
（ｋ）を高速フーリエ変換し、短時間スペクトルＸ
t（ω）（ｔ：時間を示す添字）としてエコー信号パワ
ー計算部３６及びエコー経路結合量計算用受話信号選択
部４６に出力する。
【００７１】高速フーリエ変換部３２は、送信信号ｙ
（ｋ）を高速フーリエ変換し、短時間スペクトルＹ
t（ω）として送話信号パワー計算部４２、エコー抑圧
ゲイン決定部４４及びエコー抑圧処理部４５及びエコー
経路結合量計算用送信信号選択部４７に出力する。な
お、上記送信信号ｙ（ｋ）は、送話信号ｓ（ｋ）にエコ
ー信号ｂ（ｋ）が重畳されたもの、すなわちｙ（ｋ）＝
ｓ（ｋ）＋ｂ（ｋ）として与えられるものである。
（ｋ）を高速フーリエ変換し、短時間スペクトルＹ
t（ω）として送話信号パワー計算部４２、エコー抑圧
ゲイン決定部４４及びエコー抑圧処理部４５及びエコー
経路結合量計算用送信信号選択部４７に出力する。な
お、上記送信信号ｙ（ｋ）は、送話信号ｓ（ｋ）にエコ
ー信号ｂ（ｋ）が重畳されたもの、すなわちｙ（ｋ）＝
ｓ（ｋ）＋ｂ（ｋ）として与えられるものである。
【００７２】片側発話状態検出部３４は、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を用いて、送信信号ｙ
（ｋ）に送話信号ｓ（ｋ）が含まれない状態、つまり送
信信号ｙ（ｋ）がエコー信号ｂ（ｋ）のみから構成され
る状態（時刻Ｔ）を検出し、その検出結果をエコー経路
結合量計算用受話信号選択部４６及びエコー経路結合量
計算用送信信号選択部４７に出力する。
（ｋ）及び送信信号ｙ（ｋ）を用いて、送信信号ｙ
（ｋ）に送話信号ｓ（ｋ）が含まれない状態、つまり送
信信号ｙ（ｋ）がエコー信号ｂ（ｋ）のみから構成され
る状態（時刻Ｔ）を検出し、その検出結果をエコー経路
結合量計算用受話信号選択部４６及びエコー経路結合量
計算用送信信号選択部４７に出力する。
【００７３】エコー経路結合量計算用受話信号選択部４
６（以後、単に受話信号選択部４６という）は、片側発
話状態検出部３４の検出結果に基づいて高速フーリエ変
換部３１から時々刻々と入力される受話信号ｘ（ｋ）の
短時間スペクトルＸｔ（ω）のうち、上記時刻Ｔにおけ
る短時間スペクトルＸｔ＝Ｔ（ω）を選択してエコー経
路結合量計算部３５に出力する。
６（以後、単に受話信号選択部４６という）は、片側発
話状態検出部３４の検出結果に基づいて高速フーリエ変
換部３１から時々刻々と入力される受話信号ｘ（ｋ）の
短時間スペクトルＸｔ（ω）のうち、上記時刻Ｔにおけ
る短時間スペクトルＸｔ＝Ｔ（ω）を選択してエコー経
路結合量計算部３５に出力する。
【００７４】エコー経路結合量計算用送信信号選択部４
７（以後、単に送信信号選択部４７という）は、片側発
話状態検出部３４の検出結果に基づいて高速フーリエ変
換部３２から時々刻々と入力される送信信号ｙ（ｋ）の
短時間スペクトルＹｔ（ω）のうち、上記時刻Ｔにおけ
る短時間スペクトルＹｔ＝Ｔ（ω）を選択してエコー経
路結合量計算部３５に出力する。
７（以後、単に送信信号選択部４７という）は、片側発
話状態検出部３４の検出結果に基づいて高速フーリエ変
換部３２から時々刻々と入力される送信信号ｙ（ｋ）の
短時間スペクトルＹｔ（ω）のうち、上記時刻Ｔにおけ
る短時間スペクトルＹｔ＝Ｔ（ω）を選択してエコー経
路結合量計算部３５に出力する。
【００７５】エコー経路結合量計算部３５は、時刻Ｔに
おける受話信号ｘ（ｋ）の短時間スペクトルＸ
t=T（ω）と同じく時刻Ｔにおける送信信号ｙ（ｋ）の
短時間スペクトルＹt=T（ω）とに基づいて、短時間ス
ペクトルＸt=T（ω），Ｙt=T（ω）の各振幅特性値を受
話信号パワー｜Ｘt（ω）｜2及び送信信号パワー｜Ｙ
t=T（ω）｜2（＝｜Ｅt=T（ω）｜2）として算出し、該
受話信号パワー｜Ｘt（ω）｜2と送信信号パワー｜Ｙ
t=T（ω）｜2とからエコー経路伝達特性ｈ（ｋ）の振幅
周波数関数を振幅特性値の２乗値つまりエコー経路結合
量｜Ｈe（ω）｜2として計算し、エコー信号パワー計算
部３６に出力する。
おける受話信号ｘ（ｋ）の短時間スペクトルＸ
t=T（ω）と同じく時刻Ｔにおける送信信号ｙ（ｋ）の
短時間スペクトルＹt=T（ω）とに基づいて、短時間ス
ペクトルＸt=T（ω），Ｙt=T（ω）の各振幅特性値を受
話信号パワー｜Ｘt（ω）｜2及び送信信号パワー｜Ｙ
t=T（ω）｜2（＝｜Ｅt=T（ω）｜2）として算出し、該
受話信号パワー｜Ｘt（ω）｜2と送信信号パワー｜Ｙ
t=T（ω）｜2とからエコー経路伝達特性ｈ（ｋ）の振幅
周波数関数を振幅特性値の２乗値つまりエコー経路結合
量｜Ｈe（ω）｜2として計算し、エコー信号パワー計算
部３６に出力する。
【００７６】エコー信号パワー計算部３６は、現時刻の
受話信号パワー｜Ｘt（ω）｜2に上記エコー経路結合量
｜Ｈe（ω）｜2を乗じて、エコー信号ｂ（ｋ）の短時間
スペクトルＥt（ω）の振幅特性値を予測エコー信号パ
ワー｜Ｅte（ω）｜2として算出し、遅延演算子３７及
びエコー信号パワー決定部３８に出力する。
受話信号パワー｜Ｘt（ω）｜2に上記エコー経路結合量
｜Ｈe（ω）｜2を乗じて、エコー信号ｂ（ｋ）の短時間
スペクトルＥt（ω）の振幅特性値を予測エコー信号パ
ワー｜Ｅte（ω）｜2として算出し、遅延演算子３７及
びエコー信号パワー決定部３８に出力する。
【００７７】遅延演算子３７は、予測エコー信号パワー
｜Ｅte（ω）｜2を単位時間遅延させるものであり、予
測エコー信号パワー｜Ｅte（ω）｜2の入力に対して１
単位時間前（ｔ−１）の予測エコー信号パワー｜Ｅt-1e
（ω）｜2をエコー信号パワー決定部３８に出力する。
｜Ｅte（ω）｜2を単位時間遅延させるものであり、予
測エコー信号パワー｜Ｅte（ω）｜2の入力に対して１
単位時間前（ｔ−１）の予測エコー信号パワー｜Ｅt-1e
（ω）｜2をエコー信号パワー決定部３８に出力する。
【００７８】エコー信号パワー決定部３８は、上記予測
エコー信号パワー｜Ｅte（ω）｜2と１単位時間前の予
測エコー信号パワー｜Ｅt-1e（ω）｜2とを比較し、よ
り大きい方の値を予測エコー信号パワー｜Ｅ’te（ω）
｜2に決定し、聴感補正部４０に出力する。
エコー信号パワー｜Ｅte（ω）｜2と１単位時間前の予
測エコー信号パワー｜Ｅt-1e（ω）｜2とを比較し、よ
り大きい方の値を予測エコー信号パワー｜Ｅ’te（ω）
｜2に決定し、聴感補正部４０に出力する。
【００７９】一方、聴感補正ゲイン記憶部３９は、聴感
補正ゲインＦ（ω）を記憶するものである。この聴感補
正ゲインＦ（ω）は、エコーの物理的（客観的）な周波
数特性を人間の主観的な（聴感上聞こえる）周波数特性
に補正するゲインである。上記予測エコー信号パワー｜
Ｅ’te（ω）｜2は客観量であり、実際に送信されて相
手側で拡声され、人によって検知されるエコーの主観量
とは異なる。
補正ゲインＦ（ω）を記憶するものである。この聴感補
正ゲインＦ（ω）は、エコーの物理的（客観的）な周波
数特性を人間の主観的な（聴感上聞こえる）周波数特性
に補正するゲインである。上記予測エコー信号パワー｜
Ｅ’te（ω）｜2は客観量であり、実際に送信されて相
手側で拡声され、人によって検知されるエコーの主観量
とは異なる。
【００８０】これは、相手側で話している声やこちらか
ら送られる声にマスキングされるといった人間の聴覚特
性上の影響を受けるためである。この聴覚特性は通話状
態に応じて変化するため、聴感補正ゲインＦ（ω）を一
義的に決定することはできないが、その代表値を決定す
ることは可能である。例えば、特願平９−５９５４９号
の所要エコー抑圧量の周波数特性の実験値から聴感補正
ゲインＦ（ω）を計算することができる。
ら送られる声にマスキングされるといった人間の聴覚特
性上の影響を受けるためである。この聴覚特性は通話状
態に応じて変化するため、聴感補正ゲインＦ（ω）を一
義的に決定することはできないが、その代表値を決定す
ることは可能である。例えば、特願平９−５９５４９号
の所要エコー抑圧量の周波数特性の実験値から聴感補正
ゲインＦ（ω）を計算することができる。
【００８１】聴感補正部４０は、上記予測エコー信号パ
ワー｜Ｅ’te（ω）｜2に聴感補正ゲインＦ（ω）を乗
算し、得られた値Ｆ（ω）｜Ｅ’te（ω）｜2を送話信
号パワー計算部４２及びエコー抑圧ゲイン決定部４４に
出力する。
ワー｜Ｅ’te（ω）｜2に聴感補正ゲインＦ（ω）を乗
算し、得られた値Ｆ（ω）｜Ｅ’te（ω）｜2を送話信
号パワー計算部４２及びエコー抑圧ゲイン決定部４４に
出力する。
【００８２】スムージングファクタ記憶部４１は、スム
ージングファクタβを記憶するものである。例えば、こ
のスムージングファクタ記憶部４１には、スムージング
ファクタβとして、０．９〜０．９９９の範囲内の値、
好ましくは「０．９９」が記憶される。このスムージン
グファクタβは、以下に説明する式（１）に示すよう
に、予測送話信号パワー｜Ｓte（ω）｜2を算出する際
に１単位時間前の予測送話信号パワー｜Ｓt-1e（ω）｜
2の影響度を設定する値である。スムージングファクタ
βが大きい程、過去の推定値すなわち予測送話信号パワ
ー｜Ｓt-1e（ω）｜2が支配的となり、予測送話信号パ
ワー｜Ｓte（ω）｜2は時間的により平均化されること
になる。
ージングファクタβを記憶するものである。例えば、こ
のスムージングファクタ記憶部４１には、スムージング
ファクタβとして、０．９〜０．９９９の範囲内の値、
好ましくは「０．９９」が記憶される。このスムージン
グファクタβは、以下に説明する式（１）に示すよう
に、予測送話信号パワー｜Ｓte（ω）｜2を算出する際
に１単位時間前の予測送話信号パワー｜Ｓt-1e（ω）｜
2の影響度を設定する値である。スムージングファクタ
βが大きい程、過去の推定値すなわち予測送話信号パワ
ー｜Ｓt-1e（ω）｜2が支配的となり、予測送話信号パ
ワー｜Ｓte（ω）｜2は時間的により平均化されること
になる。
【００８３】送話信号パワー計算部４２は、このスムー
ジングファクタβと上記聴感補正部４０から入力された
値Ｆ（ω）｜Ｅ’te（ω）｜2と遅延演算子４３から入
力された予測送話信号パワー｜Ｓt-1e（ω）｜2とを用
いて、送話信号ｓ（ｋ）の短時間スペクトルＳt（ω）
の振幅特性値を予測送話信号パワー｜Ｓte（ω）｜2と
して算出し、該予測送話信号パワー｜Ｓte（ω）｜2を
遅延演算子４３及びエコー抑圧ゲイン決定部４４に出力
する。
ジングファクタβと上記聴感補正部４０から入力された
値Ｆ（ω）｜Ｅ’te（ω）｜2と遅延演算子４３から入
力された予測送話信号パワー｜Ｓt-1e（ω）｜2とを用
いて、送話信号ｓ（ｋ）の短時間スペクトルＳt（ω）
の振幅特性値を予測送話信号パワー｜Ｓte（ω）｜2と
して算出し、該予測送話信号パワー｜Ｓte（ω）｜2を
遅延演算子４３及びエコー抑圧ゲイン決定部４４に出力
する。
【００８４】遅延演算子４３は、予測送話信号パワー｜
Ｓte（ω）｜2を単位時間遅延させるものであり、予測
送話信号パワー｜Ｓte（ω）｜2が入力されると、１単
位時間前（ｔ−１）の予測送話信号パワー｜Ｓt-1e
（ω）｜2を出力する。
Ｓte（ω）｜2を単位時間遅延させるものであり、予測
送話信号パワー｜Ｓte（ω）｜2が入力されると、１単
位時間前（ｔ−１）の予測送話信号パワー｜Ｓt-1e
（ω）｜2を出力する。
【００８５】エコー抑圧ゲイン決定部４４は、予測送話
信号パワー｜Ｓte（ω）｜2と上記聴感補正部４０から
入力された値Ｆ（ω）｜Ｅ’te（ω）｜2と高速フーリ
エ変換部３２から入力された短時間スペクトルＹ
t（ω）とに基づいてエコー抑圧ゲインＧ（ω）を算出
し、エコー抑圧処理部４５に出力する。
信号パワー｜Ｓte（ω）｜2と上記聴感補正部４０から
入力された値Ｆ（ω）｜Ｅ’te（ω）｜2と高速フーリ
エ変換部３２から入力された短時間スペクトルＹ
t（ω）とに基づいてエコー抑圧ゲインＧ（ω）を算出
し、エコー抑圧処理部４５に出力する。
【００８６】エコー抑圧処理部４５は、短時間スペクト
ルＹt（ω）にエコー抑圧ゲインＧ（ω）を乗算し、エ
コーを抑圧した処理信号Ｓte（ω）を高速逆フーリエ変
換部３３に出力する。
ルＹt（ω）にエコー抑圧ゲインＧ（ω）を乗算し、エ
コーを抑圧した処理信号Ｓte（ω）を高速逆フーリエ変
換部３３に出力する。
【００８７】また、高速逆フーリエ変換部３３は、処理
信号Ｙte（ω）を逆フーリエ変換することにより、エコ
ー信号ｂ（ｋ）が抑圧され送話信号ｓ（ｋ）が強調され
た時間信号ｓe（ｋ）を伝送路９に出力する。
信号Ｙte（ω）を逆フーリエ変換することにより、エコ
ー信号ｂ（ｋ）が抑圧され送話信号ｓ（ｋ）が強調され
た時間信号ｓe（ｋ）を伝送路９に出力する。
【００８８】本実施形態では、上記機能構成のエコー抑
圧装置をエコー抑圧プログラムをコンピュータ装置に組
み込むことによって実現する。この場合、コンピュータ
装置は、受話信号ｘ（ｋ）と送信信号ｙ（ｋ）とを入力
信号とし、該受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）に
エコー抑圧プログラムに基づいて上記図１に示す流れで
処理することにより、エコー信号ｂ（ｋ）が抑圧された
時間信号ｓe（ｋ）を生成する。
圧装置をエコー抑圧プログラムをコンピュータ装置に組
み込むことによって実現する。この場合、コンピュータ
装置は、受話信号ｘ（ｋ）と送信信号ｙ（ｋ）とを入力
信号とし、該受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）に
エコー抑圧プログラムに基づいて上記図１に示す流れで
処理することにより、エコー信号ｂ（ｋ）が抑圧された
時間信号ｓe（ｋ）を生成する。
【００８９】次に、上記エコー抑圧装置の動作について
詳しく説明する。まず、片側発話状態検出部３４は、受
話信号ｘ（ｋ）と送信信号ｙ（ｋ）とを常時監視し、受
話信号だけの状態つまり送話信号ｓ（ｋ）が送信されて
いない状態を検出すると、該検出の時刻Ｔを受話信号選
択部４６に出力する。この結果、受話信号選択部４６
は、受話信号ｘ（ｋ）の短時間スペクトルＸt（ω）の
うち、当該時刻Ｔにおける短時間スペクトルＸ
t=T（ω）を選択しエコー経路結合量計算部３５に出力
する。
詳しく説明する。まず、片側発話状態検出部３４は、受
話信号ｘ（ｋ）と送信信号ｙ（ｋ）とを常時監視し、受
話信号だけの状態つまり送話信号ｓ（ｋ）が送信されて
いない状態を検出すると、該検出の時刻Ｔを受話信号選
択部４６に出力する。この結果、受話信号選択部４６
は、受話信号ｘ（ｋ）の短時間スペクトルＸt（ω）の
うち、当該時刻Ｔにおける短時間スペクトルＸ
t=T（ω）を選択しエコー経路結合量計算部３５に出力
する。
【００９０】エコー経路結合量計算部３５は、このよう
にして選択的に入力された時刻Ｔの受話信号ｘ（ｋ）の
短時間スペクトルＸt=T（ω）について受話信号パワー
｜Ｘt（ω）｜2を算出すると共に、送信信号ｙ（ｋ）の
送信信号パワー｜Ｙt=T（ω）｜2、すなわちエコー信号
パワー｜Ｅt=T（ω）｜2を算出し、該エコー信号パワー
｜Ｅt=T（ω）｜2と受話信号パワー｜Ｘt=T（ω）｜2と
に基づいてエコー経路結合量｜Ｈe（ω）｜2を算出す
る。
にして選択的に入力された時刻Ｔの受話信号ｘ（ｋ）の
短時間スペクトルＸt=T（ω）について受話信号パワー
｜Ｘt（ω）｜2を算出すると共に、送信信号ｙ（ｋ）の
送信信号パワー｜Ｙt=T（ω）｜2、すなわちエコー信号
パワー｜Ｅt=T（ω）｜2を算出し、該エコー信号パワー
｜Ｅt=T（ω）｜2と受話信号パワー｜Ｘt=T（ω）｜2と
に基づいてエコー経路結合量｜Ｈe（ω）｜2を算出す
る。
【００９１】エコー信号ｂ（ｋ）は、受話信号ｘ（ｋ）
がエコー経路伝達特性（インパルス応答）ｈ（ｋ）とし
て与えられるエコー経路（受話スピーカ４と送話用マイ
クロホン３の問の伝達特）を伝達した信号、つまりｂ
（ｋ）＝ｈ（ｋ）・ｘ（ｋ）として与えられる。したが
って、エコー経路伝達特性ｈ（ｋ）のパワーであるエコ
ー経路結合量｜Ｈe（ω）｜2は、エコー信号パワー｜Ｅ
t=T（ω）｜2を受話信号パワー｜Ｘt=T（ω）｜2によっ
て除算することによって算出される。
がエコー経路伝達特性（インパルス応答）ｈ（ｋ）とし
て与えられるエコー経路（受話スピーカ４と送話用マイ
クロホン３の問の伝達特）を伝達した信号、つまりｂ
（ｋ）＝ｈ（ｋ）・ｘ（ｋ）として与えられる。したが
って、エコー経路伝達特性ｈ（ｋ）のパワーであるエコ
ー経路結合量｜Ｈe（ω）｜2は、エコー信号パワー｜Ｅ
t=T（ω）｜2を受話信号パワー｜Ｘt=T（ω）｜2によっ
て除算することによって算出される。
【００９２】これ以降の処理では、送話信号ｓ（ｋ）が
存在しない状態における（すなわち時刻Ｔにおける）受
話信号パワー｜Ｘt=T（ω）｜2とエコー信号パワー｜Ｅ
t=T（ω）｜2とから算出されたエコー経路結合量｜Ｈe
（ω）｜2に基づいて各種演算が行われる。
存在しない状態における（すなわち時刻Ｔにおける）受
話信号パワー｜Ｘt=T（ω）｜2とエコー信号パワー｜Ｅ
t=T（ω）｜2とから算出されたエコー経路結合量｜Ｈe
（ω）｜2に基づいて各種演算が行われる。
【００９３】まず、エコー信号パワー計算部３６は、現
在時刻ｔにおける受話信号ｘ（ｋ）の短時間スペクトル
Ｘｔ（ω）から受話信号パワー｜Ｘｔ（ω）｜２を算出
し、該受話信号パワー｜Ｘｔ（ω）｜２にエコー経路結
合量｜Ｈe（ω）｜２を乗じることによりエコー信号パ
ワー｜Ｅｔe（ω）｜２を算出する。この予測エコー信
号パワー｜Ｅｔe（ω）｜２は、時々刻々と変化するエ
コー経路伝達特性ｈ（ｋ）に対して、送話信号ｓ（ｋ）
が存在しない上記時刻Ｔについて求めたエコー経路結合
量｜Ｈe（ω）｜２を用いて現在時刻ｔにおけるエコー
信号のパワーを推定するものである。
在時刻ｔにおける受話信号ｘ（ｋ）の短時間スペクトル
Ｘｔ（ω）から受話信号パワー｜Ｘｔ（ω）｜２を算出
し、該受話信号パワー｜Ｘｔ（ω）｜２にエコー経路結
合量｜Ｈe（ω）｜２を乗じることによりエコー信号パ
ワー｜Ｅｔe（ω）｜２を算出する。この予測エコー信
号パワー｜Ｅｔe（ω）｜２は、時々刻々と変化するエ
コー経路伝達特性ｈ（ｋ）に対して、送話信号ｓ（ｋ）
が存在しない上記時刻Ｔについて求めたエコー経路結合
量｜Ｈe（ω）｜２を用いて現在時刻ｔにおけるエコー
信号のパワーを推定するものである。
【００９４】エコー信号パワー決定部３８は、このよう
にして推定されたエコー信号パワー｜Ｅte（ω）｜2と
該エコー信号パワー｜Ｅte（ω）｜2に対して１単位時
間前のエコー信号パワー｜Ｅt-1e（ω）｜2とを比較
し、より大きい方の値を予測エコー信号パワー｜Ｅ’te
（ω）｜2に決定する。
にして推定されたエコー信号パワー｜Ｅte（ω）｜2と
該エコー信号パワー｜Ｅte（ω）｜2に対して１単位時
間前のエコー信号パワー｜Ｅt-1e（ω）｜2とを比較
し、より大きい方の値を予測エコー信号パワー｜Ｅ’te
（ω）｜2に決定する。
【００９５】この予測エコー信号パワー｜Ｅ’te（ω）
｜2は、聴感補正部４０において聴聴感補正ゲインＦ
（ω）が乗算され聴感補正値Ｆ（ω）｜Ｅ’te（ω）｜
2に変換される。そして、送話信号パワー計算部４２で
は、上記聴感補正値Ｆ（ω）｜Ｅ’te（ω）｜2とスム
ージングファクタβと１単位時間前の予測送話信号パワ
ー｜Ｓt-1e（ω）｜2とを用いて、予測送話信号パワー
｜Ｓte（ω）｜2を算出する。この場合、送話信号パワ
ー計算部４２は、下式（１）に基づいて予測送話信号パ
ワー｜Ｓte（ω）｜2を算出する。ここで、Ｐは、信号
の正の成分だけを取り出す半波整流関数である。
｜2は、聴感補正部４０において聴聴感補正ゲインＦ
（ω）が乗算され聴感補正値Ｆ（ω）｜Ｅ’te（ω）｜
2に変換される。そして、送話信号パワー計算部４２で
は、上記聴感補正値Ｆ（ω）｜Ｅ’te（ω）｜2とスム
ージングファクタβと１単位時間前の予測送話信号パワ
ー｜Ｓt-1e（ω）｜2とを用いて、予測送話信号パワー
｜Ｓte（ω）｜2を算出する。この場合、送話信号パワ
ー計算部４２は、下式（１）に基づいて予測送話信号パ
ワー｜Ｓte（ω）｜2を算出する。ここで、Ｐは、信号
の正の成分だけを取り出す半波整流関数である。
【００９６】
【数３７】
【００９７】このようにして予測送話信号パワー｜Ｓte
（ω）｜2が算出されると、エコー抑圧ゲイン決定部４
４は、該予測送話信号パワー｜Ｓte（ω）｜2と上記聴
感補正値Ｆ（ω）｜Ｅ’te（ω）｜2及び送信信号ｙ
（ｋ）の短時間スペクトルＹt（ω）に基づいてエコー
抑圧ゲインＧ（ω）を算出する。例えば、エコー抑圧ゲ
イン決定部４４は、エコー抑圧ゲインＧ（ω）を下式
（２）に基づいてＧPE（ω）を算出し、このＧPE（ω）
をエコー抑圧ゲインＧ（ω）とする。
（ω）｜2が算出されると、エコー抑圧ゲイン決定部４
４は、該予測送話信号パワー｜Ｓte（ω）｜2と上記聴
感補正値Ｆ（ω）｜Ｅ’te（ω）｜2及び送信信号ｙ
（ｋ）の短時間スペクトルＹt（ω）に基づいてエコー
抑圧ゲインＧ（ω）を算出する。例えば、エコー抑圧ゲ
イン決定部４４は、エコー抑圧ゲインＧ（ω）を下式
（２）に基づいてＧPE（ω）を算出し、このＧPE（ω）
をエコー抑圧ゲインＧ（ω）とする。
【００９８】
【数３８】
【００９９】また、エコー抑圧ゲインＧ（ω）として
は、式（３）に基づいて得られたＧW（ω）を適用する
ことが可能であると共に、式（４）に基づいて得られた
ＧML（ω）をエコー抑圧ゲインＧ（ω）に適用すること
も考えられる。
は、式（３）に基づいて得られたＧW（ω）を適用する
ことが可能であると共に、式（４）に基づいて得られた
ＧML（ω）をエコー抑圧ゲインＧ（ω）に適用すること
も考えられる。
【０１００】
【数３９】
【０１０１】
【数４０】
【０１０２】さらに、エコー抑圧ゲインＧ（ω）として
は、予測送話信号パワー｜Ｓte（ω）｜2と予測エコー
信号パワー｜Ｅ’te（ω）｜2を用いて式（５）によっ
て中間変数υ（ω）を算出し、０次及び１次の修正ベッ
セル関数Ｉ0，Ｉ1を用いて式（６）によって算出された
ＧEM（ω）を適用することが考えられる。
は、予測送話信号パワー｜Ｓte（ω）｜2と予測エコー
信号パワー｜Ｅ’te（ω）｜2を用いて式（５）によっ
て中間変数υ（ω）を算出し、０次及び１次の修正ベッ
セル関数Ｉ0，Ｉ1を用いて式（６）によって算出された
ＧEM（ω）を適用することが考えられる。
【０１０３】
【数４１】
【０１０４】
【数４２】
【０１０５】このようにして、送話信号ｓ（ｋ）が存在
しない時刻Ｔにおけるエコー経路結合量｜Ｈe（ω）｜2
を用いて得られた予測送話信号パワー｜Ｓte（ω）｜2
と聴感補正値Ｆ（ω）｜Ｅ’te（ω）｜2、及び現在時
刻ｔにおける送信信号ｙ（ｋ）の短時間スペクトルＹt
（ω）に基づいて、エコー抑圧ゲインＧ（ω）が算出さ
れる。
しない時刻Ｔにおけるエコー経路結合量｜Ｈe（ω）｜2
を用いて得られた予測送話信号パワー｜Ｓte（ω）｜2
と聴感補正値Ｆ（ω）｜Ｅ’te（ω）｜2、及び現在時
刻ｔにおける送信信号ｙ（ｋ）の短時間スペクトルＹt
（ω）に基づいて、エコー抑圧ゲインＧ（ω）が算出さ
れる。
【０１０６】エコー抑圧処理部４５では、上記エコー抑
圧ゲインＧ（ω）を現在時刻ｔの送信信号ｙ（ｋ）の短
時間スペクトルＹt（ω）に乗算することにより、当該
短時間スペクトルＹt（ω）に含まれるエコー信号成分
が抑圧されると共に送話信号成分が強調された処理信号
Ｓte（ω）を生成する。そして、この処理信号Ｓte
（ω）は、高速逆フーリエ変換部３３によって時間信号
ｓe（ｋ）に変換されて伝送路９に送出される。
圧ゲインＧ（ω）を現在時刻ｔの送信信号ｙ（ｋ）の短
時間スペクトルＹt（ω）に乗算することにより、当該
短時間スペクトルＹt（ω）に含まれるエコー信号成分
が抑圧されると共に送話信号成分が強調された処理信号
Ｓte（ω）を生成する。そして、この処理信号Ｓte
（ω）は、高速逆フーリエ変換部３３によって時間信号
ｓe（ｋ）に変換されて伝送路９に送出される。
【０１０７】このように、本実施形態では、エコー抑圧
ゲインＧ（ω）を決定するエコー経路結合量｜Ｈe
（ω）｜2を、送話信号ｓ（ｋ）が存在しない時刻Ｔに
おける送信信号ｙ（ｋ）つまりエコー信号ｂ（ｋ）に基
づいて算出する。すなわち、従来の適応フィルクは、エ
コー経路のインパルス応答の振幅と位相の両方を推定す
るのに対して、本実施形態では、エコー経路結合量｜Ｈ
e（ω）｜2つまり振幅の二乗値のみを用いる。このた
め、エコー経路の変動つまりエコー経路伝達特性ｈ
（ｋ）の変動に対して、適応フィルタによる処理に比ベ
て位相値の変動を無視する分、ロバストにエコーを抑圧
することができる。
ゲインＧ（ω）を決定するエコー経路結合量｜Ｈe
（ω）｜2を、送話信号ｓ（ｋ）が存在しない時刻Ｔに
おける送信信号ｙ（ｋ）つまりエコー信号ｂ（ｋ）に基
づいて算出する。すなわち、従来の適応フィルクは、エ
コー経路のインパルス応答の振幅と位相の両方を推定す
るのに対して、本実施形態では、エコー経路結合量｜Ｈ
e（ω）｜2つまり振幅の二乗値のみを用いる。このた
め、エコー経路の変動つまりエコー経路伝達特性ｈ
（ｋ）の変動に対して、適応フィルタによる処理に比ベ
て位相値の変動を無視する分、ロバストにエコーを抑圧
することができる。
【０１０８】また、上記中間変数υ（ω）を用いて得ら
れたＧEM（ω）をエコー抑圧ゲインＧ（ω）に適用する
ことにより、送信信号ｙ（ｋ）に対するエコー信号ｂ
（ｋ）の割合に基づいてエコー抑圧ゲインＧ（ω）が設
定されることになる。これによって、送信信号ｙ（ｋ）
に含まれたエコー信号ｂ（ｋ）のみを抑圧することが可
能である。したがって、従来の非線形エコー抑圧処理で
問題となっている送話音声の歪みや音の途切れを最小限
に抑えてエコーのみの抑圧を行うことが可能である。
れたＧEM（ω）をエコー抑圧ゲインＧ（ω）に適用する
ことにより、送信信号ｙ（ｋ）に対するエコー信号ｂ
（ｋ）の割合に基づいてエコー抑圧ゲインＧ（ω）が設
定されることになる。これによって、送信信号ｙ（ｋ）
に含まれたエコー信号ｂ（ｋ）のみを抑圧することが可
能である。したがって、従来の非線形エコー抑圧処理で
問題となっている送話音声の歪みや音の途切れを最小限
に抑えてエコーのみの抑圧を行うことが可能である。
【０１０９】さらに、エコー信号ｂ（ｋ）のパワーをエ
コー経路結合量｜Ｈe（ω）｜２と受話信号パワー｜Ｘ
ｔ（ω）｜２から決定しているが、これだけではエコー
の残響成分の影響を加味していない。そこで、１単位時
間前（ｔ−１）に推定されたエコー信号パワー｜Ｅ
ｔ−１e（ω）｜２と現在時刻ｔに推定されたエコー信
号パワー｜Ｅｔe（ω）｜２とを比較し、これらのうち
より大きい値をエコー信号パワーに決定している。した
がって、１単位時間前（ｔ−１）に発生したエコーの残
響成分をも考慮に入れてエコー信号パワーを決定できる
ため、エコーをより十分に抑圧することが可能である。
コー経路結合量｜Ｈe（ω）｜２と受話信号パワー｜Ｘ
ｔ（ω）｜２から決定しているが、これだけではエコー
の残響成分の影響を加味していない。そこで、１単位時
間前（ｔ−１）に推定されたエコー信号パワー｜Ｅ
ｔ−１e（ω）｜２と現在時刻ｔに推定されたエコー信
号パワー｜Ｅｔe（ω）｜２とを比較し、これらのうち
より大きい値をエコー信号パワーに決定している。した
がって、１単位時間前（ｔ−１）に発生したエコーの残
響成分をも考慮に入れてエコー信号パワーを決定できる
ため、エコーをより十分に抑圧することが可能である。
【０１１０】発話者（受話信号側の室内）に返ってくる
エコーは、その発話者の発する音声にマスキングされる
ため、その周波数特性は聴感上平坦にはならない。そこ
で、上記エコーに対して特願平９−５９５４９号で提案
されている所要エコー抑圧量に基づく聴感補正ゲインを
周波数軸上でエコー信号パワーに乗じることにより、聴
感上の重みを付けている。このため、聴感上聞こえ易い
周波数帯域のエコー成分はより大きく抑圧することが可
能となる。反対に聴感上あまり気にならないエコー成分
を含む周波数帯域を不必要に抑圧しないことにより、送
話音声の歪みを最小限に抑えることが可能である。
エコーは、その発話者の発する音声にマスキングされる
ため、その周波数特性は聴感上平坦にはならない。そこ
で、上記エコーに対して特願平９−５９５４９号で提案
されている所要エコー抑圧量に基づく聴感補正ゲインを
周波数軸上でエコー信号パワーに乗じることにより、聴
感上の重みを付けている。このため、聴感上聞こえ易い
周波数帯域のエコー成分はより大きく抑圧することが可
能となる。反対に聴感上あまり気にならないエコー成分
を含む周波数帯域を不必要に抑圧しないことにより、送
話音声の歪みを最小限に抑えることが可能である。
【０１１１】したがって、上記４つの理由からダブルト
ーク中のエコー経路の変動に対しても送話音声を歪ませ
たり、音切れさせることなく、聴感上十分にエコーを抑
圧することが可能である。
ーク中のエコー経路の変動に対しても送話音声を歪ませ
たり、音切れさせることなく、聴感上十分にエコーを抑
圧することが可能である。
【０１１２】図２は、本実施形態の効果を確認するため
の主観評価実験システムの構成を示すものである。この
主観評価実験では、通話環境条件は拡声通話系における
ダブルトークを模擬した設定とし、評定者１２は、その
ときの相手側音声の品質について評価する。具体的に
は、評価側５１では話者１０が女声の短文を拡声し、ス
ピーカ２からはそれとは異なる相手側５２の話者１１が
発した男声の短文を拡声した。
の主観評価実験システムの構成を示すものである。この
主観評価実験では、通話環境条件は拡声通話系における
ダブルトークを模擬した設定とし、評定者１２は、その
ときの相手側音声の品質について評価する。具体的に
は、評価側５１では話者１０が女声の短文を拡声し、ス
ピーカ２からはそれとは異なる相手側５２の話者１１が
発した男声の短文を拡声した。
【０１１３】この際、男声には、相手側室内５２でマイ
クロホン３に回り込んだ評価側女声のエコーが重畳され
ており、エコー抑圧処理３０は。このエコーを抑圧す
る。なお、音声の周波数帯域は１００Ｈｚ〜７ｋＨｚと
し、受話音量などはＩＴＵ勧告Ｐ．３４に従った。ま
た、相手側室内５２の残響時間は、約１５０ｍｓであ
る。
クロホン３に回り込んだ評価側女声のエコーが重畳され
ており、エコー抑圧処理３０は。このエコーを抑圧す
る。なお、音声の周波数帯域は１００Ｈｚ〜７ｋＨｚと
し、受話音量などはＩＴＵ勧告Ｐ．３４に従った。ま
た、相手側室内５２の残響時間は、約１５０ｍｓであ
る。
【０１１４】図３は、上記主観評価実験システムにおい
て、従来及び本実施形態によって処理を行なた評価音声
に対する実験結果のオピニオン平均値（ＭＯＳ）を示す
ものである。ここで、本実施形態の処理を行う際のＦＦ
Ｔ点数は５１２、通信条件としての伝送遅延は０ｍｓ
（処理遅延のみ）、音響結合量は−３ｄＢに設定した。
この図から明らかなように、本実施形態は、従来方式に
比ベてエコー抑圧の処理品質が主観評価の上においても
改善されていることが分かる。
て、従来及び本実施形態によって処理を行なた評価音声
に対する実験結果のオピニオン平均値（ＭＯＳ）を示す
ものである。ここで、本実施形態の処理を行う際のＦＦ
Ｔ点数は５１２、通信条件としての伝送遅延は０ｍｓ
（処理遅延のみ）、音響結合量は−３ｄＢに設定した。
この図から明らかなように、本実施形態は、従来方式に
比ベてエコー抑圧の処理品質が主観評価の上においても
改善されていることが分かる。
【０１１５】なお、本実施形態は、エコー抑圧プログラ
ムに基づくソフトウエア的な手段によってエコー抑圧装
置を構成したが、本発明はこれに限定されるものではな
い。例えば、図１に示す機能構成を回路素子によってハ
ードウエア的に構成することも可能である。
ムに基づくソフトウエア的な手段によってエコー抑圧装
置を構成したが、本発明はこれに限定されるものではな
い。例えば、図１に示す機能構成を回路素子によってハ
ードウエア的に構成することも可能である。
【０１１６】また、上記実施形態では、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルに変換
する手段として高速フーリエ変換（ＦＦＴ）を採用して
いる。しかし、本発明はこれに限定されるものではな
く、離散コサイン変換、ウェーブレット変換あるいは周
波数分割回路を用いて各週端数帯域に分割することによ
り、受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）を短時間ス
ペクトルに変換することが考えられる。
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルに変換
する手段として高速フーリエ変換（ＦＦＴ）を採用して
いる。しかし、本発明はこれに限定されるものではな
く、離散コサイン変換、ウェーブレット変換あるいは周
波数分割回路を用いて各週端数帯域に分割することによ
り、受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）を短時間ス
ペクトルに変換することが考えられる。
【０１１７】〔第２実施形態〕次に、図４〜図７を参照
して、本発明の第２実施形態について説明する。この第
２実施形態は、上記第１実施形態に図９に示したエコー
キャンセラ２１とほぼ同等のエコーキャンセラを組み合
わせることにより、残留エコーをさらに抑圧するもので
ある。なお、以下の説明では、図１及び図９に示した構
成要素と同等の構成要素については説明を省略する。
して、本発明の第２実施形態について説明する。この第
２実施形態は、上記第１実施形態に図９に示したエコー
キャンセラ２１とほぼ同等のエコーキャンセラを組み合
わせることにより、残留エコーをさらに抑圧するもので
ある。なお、以下の説明では、図１及び図９に示した構
成要素と同等の構成要素については説明を省略する。
【０１１８】図４は、本実施形態におけるエコー抑圧装
置の機能構成を示すブロック図である。この図におい
て、符号６０はエコーキャンセラ、６１はエコー抑圧
部、６２はエコー経路変動検出部、６３はエコー経路結
合量計算部（第２のエコー経路伝達関数算出手段）、６
４は片側発話状態検出部である。
置の機能構成を示すブロック図である。この図におい
て、符号６０はエコーキャンセラ、６１はエコー抑圧
部、６２はエコー経路変動検出部、６３はエコー経路結
合量計算部（第２のエコー経路伝達関数算出手段）、６
４は片側発話状態検出部である。
【０１１９】エコーキャンセラ６０は、図５に示すよう
に、図９に示した従来のエコーキャンセラ２１とほぼ同
等の構成を備えるものであるが、エコー経路推定回路２
３’の構成が従来のエコー経路推定回路２３とは異なっ
ている。すなわち、本実施形態のエコー経路推定回路２
３’は、一種の適応フィルタであり、エコー経路変動検
出部６２から出力される制御信号に基づいて疑似エコー
経路ｈe（ｋ）を推定すると共に、該疑似エコー経路ｈe
（ｋ）を疑似エコー経路畳み込み演算回路２２に出力す
るのみではなく、エコー経路結合量計算部６３にも出力
するように構成されている。
に、図９に示した従来のエコーキャンセラ２１とほぼ同
等の構成を備えるものであるが、エコー経路推定回路２
３’の構成が従来のエコー経路推定回路２３とは異なっ
ている。すなわち、本実施形態のエコー経路推定回路２
３’は、一種の適応フィルタであり、エコー経路変動検
出部６２から出力される制御信号に基づいて疑似エコー
経路ｈe（ｋ）を推定すると共に、該疑似エコー経路ｈe
（ｋ）を疑似エコー経路畳み込み演算回路２２に出力す
るのみではなく、エコー経路結合量計算部６３にも出力
するように構成されている。
【０１２０】このようなエコーキャンセラ６０は、エコ
ー信号ｂ（ｋ）と送話信号ｓ（ｋ）との和として与えら
れる送信信号ｙ（ｋ）（＝ｂ（ｋ）＋ｓ（ｋ））から疑
似エコー信号ｂe（ｋ）を減算することにより、エコー
信号ｂ（ｋ）を減衰させたエコー低減信号ｅ（ｋ）（＝
ｙ（ｋ）−ｂe（ｋ）＝ｂd（ｋ）＋ｓ（ｋ））を上記エ
コー抑圧部６１とエコー経路変動検出部６２と片側発話
状態検出部６４とに出力する。ここで、信号ｂd（ｋ）
は、エコー信号ｂ（ｋ）から疑似エコー信号ｂe（ｋ）
を差し引いて得られる残留エコー信号である。
ー信号ｂ（ｋ）と送話信号ｓ（ｋ）との和として与えら
れる送信信号ｙ（ｋ）（＝ｂ（ｋ）＋ｓ（ｋ））から疑
似エコー信号ｂe（ｋ）を減算することにより、エコー
信号ｂ（ｋ）を減衰させたエコー低減信号ｅ（ｋ）（＝
ｙ（ｋ）−ｂe（ｋ）＝ｂd（ｋ）＋ｓ（ｋ））を上記エ
コー抑圧部６１とエコー経路変動検出部６２と片側発話
状態検出部６４とに出力する。ここで、信号ｂd（ｋ）
は、エコー信号ｂ（ｋ）から疑似エコー信号ｂe（ｋ）
を差し引いて得られる残留エコー信号である。
【０１２１】エコー抑圧部６１は、第１実施形態におけ
る高速フーリエ変換部３１，３２、高速逆フーリエ変換
部３３、エコー信号パワー計算部３６、遅延演算子３
７，４３、エコー信号パワー決定部３８、聴感補正ゲイ
ン記憶部３９、聴感補正部４０、スムージングファクタ
記憶部４１、送話信号パワー計算部４２、エコー抑圧ゲ
イン決定部４４、エコー抑圧処理部４５、エコー経路結
合量計算用受話信号選択部４６、エコー経路結合量計算
用送信信号選択部４７をまとめたものであると共に、上
記エコー経路結合量計算部３５とは異なる新たなエコー
経路結合量計算部３５’を具備するものである。
る高速フーリエ変換部３１，３２、高速逆フーリエ変換
部３３、エコー信号パワー計算部３６、遅延演算子３
７，４３、エコー信号パワー決定部３８、聴感補正ゲイ
ン記憶部３９、聴感補正部４０、スムージングファクタ
記憶部４１、送話信号パワー計算部４２、エコー抑圧ゲ
イン決定部４４、エコー抑圧処理部４５、エコー経路結
合量計算用受話信号選択部４６、エコー経路結合量計算
用送信信号選択部４７をまとめたものであると共に、上
記エコー経路結合量計算部３５とは異なる新たなエコー
経路結合量計算部３５’を具備するものである。
【０１２２】本実施形態のエコー経路結合量計算部３
５’は、図６に示すように、エコー経路変動検出部６２
から入力される制御信号に基づいて、ダブルトーク状態
でエコー経路ｈ（ｋ）が変動する状態以外（つまり受話
状態、送話状態あるいはダブルトーク状態でエコー経路
が変動しない状態）における残留エコー経路結合量｜Ｈ
d（ω）｜2（残留振幅周波数関数）を受話信号ｘ（ｋ）
の短時間スペクトルＸt= T（ω）とエコーキャンセラ６
０の出力信号ｅ（ｋ）の短時間スペクトルＥt=T（ω）
（つまり残留エコー信号ｂd（ｋ）の短時間スペクトル
Ｂd（ω））とを用いて算出するものである。
５’は、図６に示すように、エコー経路変動検出部６２
から入力される制御信号に基づいて、ダブルトーク状態
でエコー経路ｈ（ｋ）が変動する状態以外（つまり受話
状態、送話状態あるいはダブルトーク状態でエコー経路
が変動しない状態）における残留エコー経路結合量｜Ｈ
d（ω）｜2（残留振幅周波数関数）を受話信号ｘ（ｋ）
の短時間スペクトルＸt= T（ω）とエコーキャンセラ６
０の出力信号ｅ（ｋ）の短時間スペクトルＥt=T（ω）
（つまり残留エコー信号ｂd（ｋ）の短時間スペクトル
Ｂd（ω））とを用いて算出するものである。
【０１２３】また、エコー経路結合量計算部３５’は、
上記各短時間スペクトルＸｔ＝Ｔ（ω），Ｅ
ｔ＝Ｔ（ω）が入力されている期間においては、その振
幅特性｜Ｘｔ＝Ｔ（ω）｜，｜Ｅｔ＝Ｔ（ω）｜から残
留エコー経路結合量｜Ｈd（ω）｜２を順次算出してエ
コー信号パワー計算部３６に出力する。また、これ以外
の期間において、エコー経路結合量計算部３５’は、エ
コー経路変動検出部６２から制御信号が入力されない限
り、最も最近計算した残留エコー経路結合量｜Ｈd
（ω）｜２をエコー信号パワー計算部３６に出力する。
上記各短時間スペクトルＸｔ＝Ｔ（ω），Ｅ
ｔ＝Ｔ（ω）が入力されている期間においては、その振
幅特性｜Ｘｔ＝Ｔ（ω）｜，｜Ｅｔ＝Ｔ（ω）｜から残
留エコー経路結合量｜Ｈd（ω）｜２を順次算出してエ
コー信号パワー計算部３６に出力する。また、これ以外
の期間において、エコー経路結合量計算部３５’は、エ
コー経路変動検出部６２から制御信号が入力されない限
り、最も最近計算した残留エコー経路結合量｜Ｈd
（ω）｜２をエコー信号パワー計算部３６に出力する。
【０１２４】また、エコー経路結合量計算部３５’は、
エコー経路変動検出部６２からダブルトーク状態かつエ
コー経路が変動する状態である旨の制御信号が入力され
ると、全ての処理を中止してエコー信号パワー計算部３
６への残留エコー経路結合量｜Ｈd（ω）｜2の出力を停
止する。したがって、残留エコー経路結合量｜Ｈd
（ω）｜2は、ダブルトーク状態かつエコー経路ｈ
（ｋ）が変動する状態以外、つまり受話状態、送話状態
あるいはダブルトーク状態でエコー経路ｈ（ｋ）が変動
しない場合のみに、エコー信号パワー計算部３６に供給
されることになる。
エコー経路変動検出部６２からダブルトーク状態かつエ
コー経路が変動する状態である旨の制御信号が入力され
ると、全ての処理を中止してエコー信号パワー計算部３
６への残留エコー経路結合量｜Ｈd（ω）｜2の出力を停
止する。したがって、残留エコー経路結合量｜Ｈd
（ω）｜2は、ダブルトーク状態かつエコー経路ｈ
（ｋ）が変動する状態以外、つまり受話状態、送話状態
あるいはダブルトーク状態でエコー経路ｈ（ｋ）が変動
しない場合のみに、エコー信号パワー計算部３６に供給
されることになる。
【０１２５】このようなエコー経路結合量計算部３５’
を具備するエコー抑圧部６１は、受話信号ｘ（ｋ）及び
エコーキャンセラ６０から出力されたエコー低減信号ｅ
（ｋ）に基づいて時間信号ｓe（ｋ）を生成して伝送路
に出力する。すなわち、エコー抑圧部６１は、エコー低
減信号ｅ（ｋ）を第１実施形態の送信信号ｙ（ｋ）と見
なし、またエコー低減信号ｅ（ｋ）に重畳されている残
留エコー信号ｂd（ｋ）を第１実施形態のエコー信号ｂ
（ｋ）と見なし、さらにエコー経路ｈ（ｋ）にエコーキ
ャンセラ６０を含めた経路ｇ（ｋ）を第１実施形態のエ
コー経路ｈ（ｋ）と見なして、第１実施形態と同様に残
留エコー信号ｂd（ｋ）を抑圧した時間信号ｓe（ｋ）を
生成する。
を具備するエコー抑圧部６１は、受話信号ｘ（ｋ）及び
エコーキャンセラ６０から出力されたエコー低減信号ｅ
（ｋ）に基づいて時間信号ｓe（ｋ）を生成して伝送路
に出力する。すなわち、エコー抑圧部６１は、エコー低
減信号ｅ（ｋ）を第１実施形態の送信信号ｙ（ｋ）と見
なし、またエコー低減信号ｅ（ｋ）に重畳されている残
留エコー信号ｂd（ｋ）を第１実施形態のエコー信号ｂ
（ｋ）と見なし、さらにエコー経路ｈ（ｋ）にエコーキ
ャンセラ６０を含めた経路ｇ（ｋ）を第１実施形態のエ
コー経路ｈ（ｋ）と見なして、第１実施形態と同様に残
留エコー信号ｂd（ｋ）を抑圧した時間信号ｓe（ｋ）を
生成する。
【０１２６】片側発話状態検出部６４は、送信信号ｙ
（ｋ）とエコー低減信号ｅ（ｋ）との比であるエコー低
減量を算出する。そして、このエコー低減量に基づいて
受話信号ｘ（ｋ）が存在する状態においてエコー低減量
が時間的に単調増大していると判断した場合は、送信信
号ｙ（ｋ）中に送話信号ｓ（ｋ）が含まれていない状態
（片側発話状態）であると判定し、これ以外の状態の場
合には送話信号ｓ（ｋ）が含まれる状態（ダブルトーク
状態）であると判定して、この判定結果を制御信号とし
て上記エコー抑圧部６１内のエコー経路結合量計算用受
話信号選択部４６とエコー経路結合量計算用送信信号選
択部４７とに出力する。
（ｋ）とエコー低減信号ｅ（ｋ）との比であるエコー低
減量を算出する。そして、このエコー低減量に基づいて
受話信号ｘ（ｋ）が存在する状態においてエコー低減量
が時間的に単調増大していると判断した場合は、送信信
号ｙ（ｋ）中に送話信号ｓ（ｋ）が含まれていない状態
（片側発話状態）であると判定し、これ以外の状態の場
合には送話信号ｓ（ｋ）が含まれる状態（ダブルトーク
状態）であると判定して、この判定結果を制御信号とし
て上記エコー抑圧部６１内のエコー経路結合量計算用受
話信号選択部４６とエコー経路結合量計算用送信信号選
択部４７とに出力する。
【０１２７】本実施形態では、このような片側発話状態
検出部６４を備えているので、上述した第１実施形態に
比較して、より精度良く送信信号ｙ（ｋ）に送話信号ｓ
（ｋ）が含まれていない状態を検出することができる。
検出部６４を備えているので、上述した第１実施形態に
比較して、より精度良く送信信号ｙ（ｋ）に送話信号ｓ
（ｋ）が含まれていない状態を検出することができる。
【０１２８】エコー経路変動検出部６２は、受話信号ｘ
（ｋ）、送信信号ｙ（ｋ）及びエコー減衰信号ｅ（ｋ）
に基づいて、受話信号ｘ（ｋ）と送話信号ｓ（ｋ）が同
時に存在するダブルトーク状態状態におけるエコー経路
ｈ（ｋ）の変動を検出し、該検出結果を制御信号として
エコー経路結合量計算部３５’とエコーキャンセラ６０
内のエコー経路推定回路２３’とに出力する。
（ｋ）、送信信号ｙ（ｋ）及びエコー減衰信号ｅ（ｋ）
に基づいて、受話信号ｘ（ｋ）と送話信号ｓ（ｋ）が同
時に存在するダブルトーク状態状態におけるエコー経路
ｈ（ｋ）の変動を検出し、該検出結果を制御信号として
エコー経路結合量計算部３５’とエコーキャンセラ６０
内のエコー経路推定回路２３’とに出力する。
【０１２９】エコー経路推定回路２３’は、上記エコー
経路変動検出部６２からダブルトーク状態においてエコ
ー経路ｈ（ｋ）が変動した旨の制御信号が入力される
と、自らが推定した疑似エコー経路ｈe（ｋ）の全ての
値を「０」とする。この結果、疑似エコー経路畳み込み
演算回路２２は疑似エコー信号ｂd（ｋ）の生成を停止
するので、減算器２４は、送話用マイクロホン３から入
力された送信信号ｙ（ｋ）をそのままエコー抑圧部６１
等に出力する。
経路変動検出部６２からダブルトーク状態においてエコ
ー経路ｈ（ｋ）が変動した旨の制御信号が入力される
と、自らが推定した疑似エコー経路ｈe（ｋ）の全ての
値を「０」とする。この結果、疑似エコー経路畳み込み
演算回路２２は疑似エコー信号ｂd（ｋ）の生成を停止
するので、減算器２４は、送話用マイクロホン３から入
力された送信信号ｙ（ｋ）をそのままエコー抑圧部６１
等に出力する。
【０１３０】また、エコー経路推定回路２３’は、エコ
ー経路変動検出部６２からダブルトーク状態においてエ
コー経路ｈ（ｋ）が変動した旨の制御信号が入力される
と、当該検出結果が入力される直前に推定した疑似エコ
ー経路ｈe（ｋ）をエコー経路結合量計算部６３に出力
する。
ー経路変動検出部６２からダブルトーク状態においてエ
コー経路ｈ（ｋ）が変動した旨の制御信号が入力される
と、当該検出結果が入力される直前に推定した疑似エコ
ー経路ｈe（ｋ）をエコー経路結合量計算部６３に出力
する。
【０１３１】エコー経路結合量計算部６３は、上記疑似
エコー経路ｈe（ｋ）を高速フーリエ変換して得られた
疑似エコー経路周波数特性Ｈe（ω）の振幅特性の２乗
値である疑似エコー経路結合量｜Ｈe（ω）｜２（疑似
振幅周波数関数）を上記エコー抑圧部６１内のエコー信
号パワー計算部３６に出力する。すなわち、エコー経路
結合量計算部６３は、ダブルトーク状態においてエコー
経路ｈ（ｋ）が変動した場合、当該変動の直前に推定さ
れた疑似エコー経路ｈe（ｋ）から疑似エコー経路結合
量｜Ｈe（ω）｜２を算出してエコー信号パワー計算部
３６に出力する。
エコー経路ｈe（ｋ）を高速フーリエ変換して得られた
疑似エコー経路周波数特性Ｈe（ω）の振幅特性の２乗
値である疑似エコー経路結合量｜Ｈe（ω）｜２（疑似
振幅周波数関数）を上記エコー抑圧部６１内のエコー信
号パワー計算部３６に出力する。すなわち、エコー経路
結合量計算部６３は、ダブルトーク状態においてエコー
経路ｈ（ｋ）が変動した場合、当該変動の直前に推定さ
れた疑似エコー経路ｈe（ｋ）から疑似エコー経路結合
量｜Ｈe（ω）｜２を算出してエコー信号パワー計算部
３６に出力する。
【０１３２】この結果、エコー信号パワー計算部３６に
は、受話状態、送話状態あるいはダブルトーク状態でエ
コー経路ｈ（ｋ）が変動しないときに残留エコー経路結
合量｜Ｈd（ω）｜２が入力され、ダブルトーク状態に
おいてエコー経路ｈ（ｋ）が変動するときには、当該変
動の直前にエコー経路推定回路２３’によって推定され
た疑似エコー経路ｈe（ｋ）に基づく疑似エコー経路結
合量｜Ｈe（ω）｜２が択一的に入力される。エコー信
号パワー計算部３６は、通話状況に応じて択一的に入力
される残留エコー経路結合量｜Ｈd（ω）｜２あるいは
疑似エコー経路結合量｜Ｈe（ω）｜２に基づいて予測
エコー信号パワー｜Ｅｔe（ω）｜２を算出する。
は、受話状態、送話状態あるいはダブルトーク状態でエ
コー経路ｈ（ｋ）が変動しないときに残留エコー経路結
合量｜Ｈd（ω）｜２が入力され、ダブルトーク状態に
おいてエコー経路ｈ（ｋ）が変動するときには、当該変
動の直前にエコー経路推定回路２３’によって推定され
た疑似エコー経路ｈe（ｋ）に基づく疑似エコー経路結
合量｜Ｈe（ω）｜２が択一的に入力される。エコー信
号パワー計算部３６は、通話状況に応じて択一的に入力
される残留エコー経路結合量｜Ｈd（ω）｜２あるいは
疑似エコー経路結合量｜Ｈe（ω）｜２に基づいて予測
エコー信号パワー｜Ｅｔe（ω）｜２を算出する。
【０１３３】また、ダブルトーク状態でエコー経路ｈ
（ｋ）が変動するときには、疑似エコー経路畳み込み演
算回路２２の動作は停止されて、エコー抑圧部６１には
送話用マイクロホン３からの送信信号ｙ（ｋ）がそのま
ま入力される。すなわち、エコー抑圧部６１は、エコー
キャンセラ６０をパスした状態で動作することになる。
（ｋ）が変動するときには、疑似エコー経路畳み込み演
算回路２２の動作は停止されて、エコー抑圧部６１には
送話用マイクロホン３からの送信信号ｙ（ｋ）がそのま
ま入力される。すなわち、エコー抑圧部６１は、エコー
キャンセラ６０をパスした状態で動作することになる。
【０１３４】図７は、本実施形態の効果を示すシミュレ
ーション結果である。このシミュレーションでは、エコ
ー経路推定回路２３’（適応フィルタ）にはＮＬＭＳア
ルゴリズムを用い、そのタップ数は５１２とした。ま
た、エコー抑圧部６１の各高速フーリエ変換におけるＦ
ＦＴ点数も５１２とした。
ーション結果である。このシミュレーションでは、エコ
ー経路推定回路２３’（適応フィルタ）にはＮＬＭＳア
ルゴリズムを用い、そのタップ数は５１２とした。ま
た、エコー抑圧部６１の各高速フーリエ変換におけるＦ
ＦＴ点数も５１２とした。
【０１３５】この図７において、（ａ）は受話信号ｘ
（ｋ）、（ｂ）は送話信号ｓ（ｋ）、（ｃ）はエコー低
減信号ｅ（ｋ）（＝ｂd（ｋ）＋ｓ（ｋ））から送話信
号ｓ（ｋ）を減じた信号（すなわち残留エコー信号ｂd
（ｋ）、（ｄ）はエコー抑圧部６１の出力信号ｓe
（ｋ）から送話信号ｓ（ｋ）を減じた信号を示してい
る。また、これらの（ａ）〜（ｄ）は、各々に時間の経
過とともに、受話状態（Reseive）→送話状態（Send）
→ダブルトーク状態（Double-talk）→エコー経路の変
動（Path Change）→受話状態（Reseive）に移行する場
合を示している。
（ｋ）、（ｂ）は送話信号ｓ（ｋ）、（ｃ）はエコー低
減信号ｅ（ｋ）（＝ｂd（ｋ）＋ｓ（ｋ））から送話信
号ｓ（ｋ）を減じた信号（すなわち残留エコー信号ｂd
（ｋ）、（ｄ）はエコー抑圧部６１の出力信号ｓe
（ｋ）から送話信号ｓ（ｋ）を減じた信号を示してい
る。また、これらの（ａ）〜（ｄ）は、各々に時間の経
過とともに、受話状態（Reseive）→送話状態（Send）
→ダブルトーク状態（Double-talk）→エコー経路の変
動（Path Change）→受話状態（Reseive）に移行する場
合を示している。
【０１３６】この図に示されるように、本実施形態のエ
コー抑圧装置によれば、各通話状態においてエコー信号
を有効に抑圧していることが解る。特に、（ｃ）と
（ｄ）とを比較すると、エコー抑圧部６１は、ダブルト
ーク状態においてエコー経路の変動がない状態でエコー
低減信号ｅ（ｋ）に含まれる残留エコー信号ｂd（ｋ）
を十分に低減させると共に、ダブルトーク状態における
エコー経路の変動に対しても残留エコー信号ｂd（ｋ）
の増大を抑えていることが解る。
コー抑圧装置によれば、各通話状態においてエコー信号
を有効に抑圧していることが解る。特に、（ｃ）と
（ｄ）とを比較すると、エコー抑圧部６１は、ダブルト
ーク状態においてエコー経路の変動がない状態でエコー
低減信号ｅ（ｋ）に含まれる残留エコー信号ｂd（ｋ）
を十分に低減させると共に、ダブルトーク状態における
エコー経路の変動に対しても残留エコー信号ｂd（ｋ）
の増大を抑えていることが解る。
【０１３７】
【発明の効果】以上説明したように、本発明に係わるエ
コー抑圧方法及び装置並びにエコー抑圧プログラムが記
憶されたコンピュータに読取り可能な記憶媒体によれ
ば、以下のような効果を奏する。 （１）受話信号ｘ（ｋ）が時間変動するエコー経路伝達
関数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ
（ｋ）として送話信号ｓ（ｋ）に重畳された送信信号ｙ
（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエコー
信号ｂ（ｋ）を抑圧する方法において、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルにそれ
ぞれ変換し、送話信号ｓ（ｋ）が含まれない状態におけ
る送信信号ｙ（ｋ）の短時間スペクトル及び前記受話信
号ｘ（ｋ）の短時間スペクトルを用いてエコー経路伝達
関数ｈ（ｋ）の振幅周波数関数を算出し、該振幅周波数
関数並びに送信信号ｙ（ｋ）の短時間スペクトル及び前
記受話信号ｘ（ｋ）の短時間スペクトルを用いて、以降
の時刻におけるエコー信号ｂ（ｋ）の短時間スペクトル
の振幅特性値及び送話信号ｓ（ｋ）の短時間スペクトル
の振幅特性値を推定し、該エコー信号ｂ（ｋ）の短時間
スペクトルの振幅特性値及び送話信号ｓ（ｋ）の短時間
スペクトルの振幅特性値及び送信信号ｙ（ｋ）の短時間
スペクトルの振幅特性値に基づく演算式を用いてエコー
抑圧ゲインＧ（ω）を算出し、該エコー抑圧ゲインＧ
（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗算す
ることによりエコー信号ｂ（ｋ）を抑圧するので、ダブ
ルトーク中のエコー経路の変動に対して、送話音声を歪
ませたり、音切れさせることなく、聴感上十分にエコー
を抑圧して通話品質を良好に保つことが可能である。 （２）受話信号ｘ（ｋ）が時間変動するエコー経路伝達
関数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ
（ｋ）として送話信号ｓ（ｋ）に重畳された送信信号ｙ
（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエコー
信号ｂ（ｋ）を抑圧する方法において、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルにそれ
ぞれ変換し、ダブルトーク中のエコー経路伝達関数ｈ
（ｋ）の変動を検出し、エコー低減信号ｅ（ｋ）と受話
信号ｘ（ｋ）とに基づいてエコー経路伝達関数ｈ（ｋ）
の推定値である疑似エコー経路伝達関数ｈe（ｋ）を算
出し、該疑似エコー経路伝達関数ｈe（ｋ）を用いて送
信信号ｙ（ｋ）に含まれるエコー信号ｂ（ｋ）を抑圧し
て送話信号ｓ（ｋ）に残留エコー信号ｂd（ｋ）が含ま
れた前記エコー低減信号ｅ（ｋ）を生成し、ダブルトー
ク中のエコー経路伝達関数ｈ（ｋ）の変動が検出された
場合に、該変動検出前の疑似エコー経路伝達関数ｈe
（ｋ）から疑似振幅周波数関数｜Ｈe（ω）｜2を算出
し、ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変
動が検出されない場合に、送話信号ｓ（ｋ）が含まれな
い状態における送信信号ｙ（ｋ）の短時間スペクトル及
び前記受話信号ｘ（ｋ）の短時間スペクトルを用いて残
留振幅周波数関数｜Ｈd（ω）｜2を算出し、択一的に算
出された前記疑似振幅周波数関数｜Ｈe（ω）｜2あるい
は残留振幅周波数関数｜Ｈd（ω）｜2に基づいてエコー
信号ｂ（ｋ）あるいは残留エコー信号ｂd（ｋ）の短時
間スペクトルの振幅特性値を推定し、該エコー信号ｂ
（ｋ）あるいは残留エコー信号ｂd（ｋ）の短時間スペ
クトルの振幅特性値及び送話信号ｓ（ｋ）の短時間スペ
クトルの振幅特性値及び送信信号ｙ（ｋ）の短時間スペ
クトルの振幅特性値に基づく演算式を用いてエコー抑圧
ゲインＧ（ω）を算出し、ダブルトーク中のエコー経路
伝達関数ｈ（ｋ）の変動が検出されない場合には、前記
エコー抑圧ゲインＧ（ω）をエコー低減信号ｅ（ｋ）の
短時間スペクトルに乗算することにより残留エコー信号
ｂd（ｋ）を抑圧し、ダブルトーク中のエコー経路伝達
関数ｈ（ｋ）の変動が検出された場合には、エコー抑圧
ゲインＧ（ω）を送信信号ｙ（ｋ）の短時間スペクトル
に乗算することによりエコー信号ｂ（ｋ）を抑圧するの
で、ダブルトーク中のエコー経路の変動に対して、送話
音声を歪ませたり音切れさせることなく、エコー信号を
さらに抑圧することができる。したがって、通話品質を
さらに高品質にすることが可能である。
コー抑圧方法及び装置並びにエコー抑圧プログラムが記
憶されたコンピュータに読取り可能な記憶媒体によれ
ば、以下のような効果を奏する。 （１）受話信号ｘ（ｋ）が時間変動するエコー経路伝達
関数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ
（ｋ）として送話信号ｓ（ｋ）に重畳された送信信号ｙ
（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエコー
信号ｂ（ｋ）を抑圧する方法において、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルにそれ
ぞれ変換し、送話信号ｓ（ｋ）が含まれない状態におけ
る送信信号ｙ（ｋ）の短時間スペクトル及び前記受話信
号ｘ（ｋ）の短時間スペクトルを用いてエコー経路伝達
関数ｈ（ｋ）の振幅周波数関数を算出し、該振幅周波数
関数並びに送信信号ｙ（ｋ）の短時間スペクトル及び前
記受話信号ｘ（ｋ）の短時間スペクトルを用いて、以降
の時刻におけるエコー信号ｂ（ｋ）の短時間スペクトル
の振幅特性値及び送話信号ｓ（ｋ）の短時間スペクトル
の振幅特性値を推定し、該エコー信号ｂ（ｋ）の短時間
スペクトルの振幅特性値及び送話信号ｓ（ｋ）の短時間
スペクトルの振幅特性値及び送信信号ｙ（ｋ）の短時間
スペクトルの振幅特性値に基づく演算式を用いてエコー
抑圧ゲインＧ（ω）を算出し、該エコー抑圧ゲインＧ
（ω）を送信信号ｙ（ｋ）の短時間スペクトルに乗算す
ることによりエコー信号ｂ（ｋ）を抑圧するので、ダブ
ルトーク中のエコー経路の変動に対して、送話音声を歪
ませたり、音切れさせることなく、聴感上十分にエコー
を抑圧して通話品質を良好に保つことが可能である。 （２）受話信号ｘ（ｋ）が時間変動するエコー経路伝達
関数ｈ（ｋ）のエコー経路を伝搬し、エコー信号ｂ
（ｋ）として送話信号ｓ（ｋ）に重畳された送信信号ｙ
（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエコー
信号ｂ（ｋ）を抑圧する方法において、受話信号ｘ
（ｋ）及び送信信号ｙ（ｋ）を短時間スペクトルにそれ
ぞれ変換し、ダブルトーク中のエコー経路伝達関数ｈ
（ｋ）の変動を検出し、エコー低減信号ｅ（ｋ）と受話
信号ｘ（ｋ）とに基づいてエコー経路伝達関数ｈ（ｋ）
の推定値である疑似エコー経路伝達関数ｈe（ｋ）を算
出し、該疑似エコー経路伝達関数ｈe（ｋ）を用いて送
信信号ｙ（ｋ）に含まれるエコー信号ｂ（ｋ）を抑圧し
て送話信号ｓ（ｋ）に残留エコー信号ｂd（ｋ）が含ま
れた前記エコー低減信号ｅ（ｋ）を生成し、ダブルトー
ク中のエコー経路伝達関数ｈ（ｋ）の変動が検出された
場合に、該変動検出前の疑似エコー経路伝達関数ｈe
（ｋ）から疑似振幅周波数関数｜Ｈe（ω）｜2を算出
し、ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変
動が検出されない場合に、送話信号ｓ（ｋ）が含まれな
い状態における送信信号ｙ（ｋ）の短時間スペクトル及
び前記受話信号ｘ（ｋ）の短時間スペクトルを用いて残
留振幅周波数関数｜Ｈd（ω）｜2を算出し、択一的に算
出された前記疑似振幅周波数関数｜Ｈe（ω）｜2あるい
は残留振幅周波数関数｜Ｈd（ω）｜2に基づいてエコー
信号ｂ（ｋ）あるいは残留エコー信号ｂd（ｋ）の短時
間スペクトルの振幅特性値を推定し、該エコー信号ｂ
（ｋ）あるいは残留エコー信号ｂd（ｋ）の短時間スペ
クトルの振幅特性値及び送話信号ｓ（ｋ）の短時間スペ
クトルの振幅特性値及び送信信号ｙ（ｋ）の短時間スペ
クトルの振幅特性値に基づく演算式を用いてエコー抑圧
ゲインＧ（ω）を算出し、ダブルトーク中のエコー経路
伝達関数ｈ（ｋ）の変動が検出されない場合には、前記
エコー抑圧ゲインＧ（ω）をエコー低減信号ｅ（ｋ）の
短時間スペクトルに乗算することにより残留エコー信号
ｂd（ｋ）を抑圧し、ダブルトーク中のエコー経路伝達
関数ｈ（ｋ）の変動が検出された場合には、エコー抑圧
ゲインＧ（ω）を送信信号ｙ（ｋ）の短時間スペクトル
に乗算することによりエコー信号ｂ（ｋ）を抑圧するの
で、ダブルトーク中のエコー経路の変動に対して、送話
音声を歪ませたり音切れさせることなく、エコー信号を
さらに抑圧することができる。したがって、通話品質を
さらに高品質にすることが可能である。
【図１】 本発明の第１実施形態に係わるエコー抑圧装
置の機能構成を示すブロック図である。
置の機能構成を示すブロック図である。
【図２】 本発明の第１実施形態における主観評価実験
システムの構成を示す説明図である。
システムの構成を示す説明図である。
【図３】 本発明の第１実施形態に主観評価実験の実験
結果を示す特性図である。
結果を示す特性図である。
【図４】 本発明の第２実施形態に係わるエコー抑圧装
置の機能構成を示すブロック図である。
置の機能構成を示すブロック図である。
【図５】 本発明の第２実施形態に係わるエコーキャン
セラの機能構成を示すブロック図である。
セラの機能構成を示すブロック図である。
【図６】 本発明の第２実施形態に係わる残留エコー経
路結合量と疑似エコー経路結合量の生成を示すブロック
図である。
路結合量と疑似エコー経路結合量の生成を示すブロック
図である。
【図７】 本発明の第２実施形態の効果を示す計算機シ
ミュレーションの実験結果である。
ミュレーションの実験結果である。
【図８】 拡声通話系の構成を示す説明図である。
【図９】 従来のフルバンド方式のエコーキャンセラの
構成例を示すブロック図である。
構成例を示すブロック図である。
１，３……送話用マイクロホン
２，４……受話スピーカ
５，７……送話信号増幅器
６，８……受話信号増幅器
９……伝送路
１０……送話者
１１……受話者
２１……エコーキャンセラ
２２……疑似エコー経路畳み込み演算回路
２３……エコー経路推定回路
２４……減算器
２５……非線形エコー抑圧処理部
３１……高速フーリエ変換部（受話信号変換手段）
３２……高速フーリエ変換部（送信信号変換手段）
３３……高速逆フーリエ変換部（送話信号逆変換手段）
３４……片側発話状態検出部
３５……エコー経路結合量計算部（エコー経路伝達関数
算出手段） ３６……エコー信号パワー計算部（エコー信号振幅周波
数特性推定手段） ３７……遅延演算子（エコー信号振幅周波数特性推定手
段） ３８……エコー信号パワー決定部（エコー信号振幅周波
数特性推定手段） ３９……聴感補正ゲイン記憶部（エコー信号振幅周波数
特性推定手段，エコー信号聴感補正手段） ４０……聴感補正部（エコー信号振幅周波数特性推定手
段，エコー信号聴感補正手段） ４１……スムージングファクタ記憶部（送話信号振幅周
波数特性推定手段） ４２……送話信号パワー計算部（送話信号振幅周波数特
性推定手段） ４３……遅延演算子（送話信号振幅周波数特性推定手
段） ４４……エコー抑圧ゲイン決定部（送話信号振幅周波数
特性推定手段） ４５……エコー抑圧処理部 ４６……エコー経路結合量計算用受話信号選択部（エコ
ー経路伝達関数算出手段） ４７……エコー経路結合量計算用送信信号選択部（エコ
ー経路伝達関数算出手段） ２３’……エコー経路推定回路（適応フィルタ） ３５’……エコー経路結合量計算部（エコー経路伝達関
数算出手段） ６０……エコーキャンセラ ６１……エコー抑圧部 ６２……エコー経路変動検出部 ６３……エコー経路結合量計算部（第２のエコー経路伝
達関数算出手段） ６４……片側発話状態検出部 ｂ（ｋ）……エコー信号 ｂd（ｋ）……残留エコー信号 ｅ（ｋ）……エコー低減信号 ｜Ｅt=Te（ω）｜2……時刻Ｔにおけるエコー信号パワ
ー ｜Ｅte（ω）｜2， ｜Ｅ’te（ω）｜2……予測エコー
信号パワー Ｆ（ω）……聴感補正ゲイン Ｇ（ω）……エコー抑圧ゲイン ｜Ｈe（ω）｜2……エコー経路結合量，疑似エコー経路
結合量 ｜Ｈd（ω）｜2……残留エコー経路結合量 ｓ（ｋ）……送話信号 ｓe（ｋ）……時間信号 Ｓte（ω）……処理信号 ｜Ｓte（ω）｜2……予測送話信号パワー ｘ（ｋ）……受話信号 Ｘt（ω）……時刻ｔにおける受話信号ｘ（ｋ）の短時
間スペクトル ｜Ｘt（ω）｜2……受話信号パワー ｙ（ｋ）……送信信号 Ｙt（ω）……時刻ｔにおける送信信号ｙ（ｋ）の短時
間スペクトル ｜Ｙt（ω）｜2……送信信号のパワー β……スムージングファクタ
算出手段） ３６……エコー信号パワー計算部（エコー信号振幅周波
数特性推定手段） ３７……遅延演算子（エコー信号振幅周波数特性推定手
段） ３８……エコー信号パワー決定部（エコー信号振幅周波
数特性推定手段） ３９……聴感補正ゲイン記憶部（エコー信号振幅周波数
特性推定手段，エコー信号聴感補正手段） ４０……聴感補正部（エコー信号振幅周波数特性推定手
段，エコー信号聴感補正手段） ４１……スムージングファクタ記憶部（送話信号振幅周
波数特性推定手段） ４２……送話信号パワー計算部（送話信号振幅周波数特
性推定手段） ４３……遅延演算子（送話信号振幅周波数特性推定手
段） ４４……エコー抑圧ゲイン決定部（送話信号振幅周波数
特性推定手段） ４５……エコー抑圧処理部 ４６……エコー経路結合量計算用受話信号選択部（エコ
ー経路伝達関数算出手段） ４７……エコー経路結合量計算用送信信号選択部（エコ
ー経路伝達関数算出手段） ２３’……エコー経路推定回路（適応フィルタ） ３５’……エコー経路結合量計算部（エコー経路伝達関
数算出手段） ６０……エコーキャンセラ ６１……エコー抑圧部 ６２……エコー経路変動検出部 ６３……エコー経路結合量計算部（第２のエコー経路伝
達関数算出手段） ６４……片側発話状態検出部 ｂ（ｋ）……エコー信号 ｂd（ｋ）……残留エコー信号 ｅ（ｋ）……エコー低減信号 ｜Ｅt=Te（ω）｜2……時刻Ｔにおけるエコー信号パワ
ー ｜Ｅte（ω）｜2， ｜Ｅ’te（ω）｜2……予測エコー
信号パワー Ｆ（ω）……聴感補正ゲイン Ｇ（ω）……エコー抑圧ゲイン ｜Ｈe（ω）｜2……エコー経路結合量，疑似エコー経路
結合量 ｜Ｈd（ω）｜2……残留エコー経路結合量 ｓ（ｋ）……送話信号 ｓe（ｋ）……時間信号 Ｓte（ω）……処理信号 ｜Ｓte（ω）｜2……予測送話信号パワー ｘ（ｋ）……受話信号 Ｘt（ω）……時刻ｔにおける受話信号ｘ（ｋ）の短時
間スペクトル ｜Ｘt（ω）｜2……受話信号パワー ｙ（ｋ）……送信信号 Ｙt（ω）……時刻ｔにおける送信信号ｙ（ｋ）の短時
間スペクトル ｜Ｙt（ω）｜2……送信信号のパワー β……スムージングファクタ
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 阪内 澄宇，羽田 陽一，“短時間ス
ペクトラル振幅推定に基づいた非線形エ
コー抑圧処理の検討”，日本音響学会
平成10年度春季研究発表会 講演論文集
Ｉ，1998年 ３月17日，ｐｐ．551−
552，（３−５−11)
阪内 澄宇，佐々木 潤子，羽田 陽
一，“ＳＴＳＡ推定に基づいた非線形エ
コー抑圧と適応フィルタの組合せ”，日
本音響学会 平成10年度春秋季研究発表
会 講演論文集 Ｉ，1998年 ９月24
日，ｐｐ．607−608，（１−Ｐ−６)
(58)調査した分野(Int.Cl.7，ＤＢ名)
H04B 3/00
H04B 7/00
H04M 3/00
H04R 3/00
Claims (34)
- 【請求項１】 受話信号ｘ（ｋ）が時間変動するエコー
経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、エコー信
号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳された送信信
号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエ
コー信号ｂ（ｋ）を抑圧する方法であって、 受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）を短時間スペク
トルにそれぞれ変換し、 送話信号ｓ（ｋ）が含まれない状態における送信信号ｙ
（ｋ）の短時間スペクトル及び前記受話信号ｘ（ｋ）の
短時間スペクトルを用いてエコー経路伝達関数ｈ（ｋ）
の振幅周波数関数を算出し、 該振幅周波数関数並びに送信信号ｙ（ｋ）の短時間スペ
クトル及び前記受話信号ｘ（ｋ）の短時間スペクトルを
用いて、以降の時刻におけるエコー信号ｂ（ｋ）の短時
間スペクトルの振幅特性値及び送話信号ｓ（ｋ）の短時
間スペクトルの振幅特性値を推定し、 該エコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
及び送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値
及び送信信号ｙ（ｋ）の短時間スペクトルの振幅特性値
に基づく演算式を用いてエコー抑圧ゲインＧ（ω）を算
出し、 該エコー抑圧ゲインＧ（ω）を送信信号ｙ（ｋ）の短時
間スペクトルに乗算することによりエコー信号ｂ（ｋ）
を抑圧することを特徴とするエコー抑圧方法。 - 【請求項２】 受話信号ｘ（ｋ）が時間変動するエコー
経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、エコー信
号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳された送信信
号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳されたエ
コー信号ｂ（ｋ）を抑圧する方法であって、 受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）を短時間スペク
トルにそれぞれ変換し、 ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変動を
検出し、 エコー低減信号ｅ（ｋ）と受話信号ｘ（ｋ）とに基づい
てエコー経路伝達関数ｈ（ｋ）の推定値である疑似エコ
ー経路伝達関数ｈe（ｋ）を算出し、 該疑似エコー経路伝達関数ｈe（ｋ）を用いて送信信号
ｙ（ｋ）に含まれるエコー信号ｂ（ｋ）を抑圧して送話
信号ｓ（ｋ）に残留エコー信号ｂd（ｋ）が含まれた前
記エコー低減信号ｅ（ｋ）を生成し、 ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変動が
検出された場合に、該変動検出前の疑似エコー経路伝達
関数ｈe（ｋ）から疑似振幅周波数関数｜Ｈe（ω）｜2
を算出し、 ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変動が
検出されない場合に、送話信号ｓ（ｋ）が含まれない状
態における送信信号ｙ（ｋ）の短時間スペクトル及び前
記受話信号ｘ（ｋ）の短時間スペクトルを用いて残留振
幅周波数関数｜Ｈd（ω）｜2を算出し、 択一的に算出された前記疑似振幅周波数関数｜Ｈe
（ω）｜2あるいは残留振幅周波数関数｜Ｈd（ω）｜2
に基づいてエコー信号ｂ（ｋ）あるいは残留エコー信号
ｂd（ｋ）の短時間スペクトルの振幅特性値を推定し、 該エコー信号ｂ（ｋ）あるいは残留エコー信号ｂd
（ｋ）の短時間スペクトルの振幅特性値及び送話信号ｓ
（ｋ）の短時間スペクトルの振幅特性値及び送信信号ｙ
（ｋ）の短時間スペクトルの振幅特性値に基づく演算式
を用いてエコー抑圧ゲインＧ（ω）を算出し、 ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変動が
検出されない場合には、前記エコー抑圧ゲインＧ（ω）
をエコー低減信号ｅ（ｋ）の短時間スペクトルに乗算す
ることにより残留エコー信号ｂd（ｋ）を抑圧し、ダブ
ルトーク中のエコー経路伝達関数ｈ（ｋ）の変動が検出
された場合には、エコー抑圧ゲインＧ（ω）を送信信号
ｙ（ｋ）の短時間スペクトルに乗算することによりエコ
ー信号ｂ（ｋ）を抑圧することを特徴とするエコー抑圧
方法。 - 【請求項３】 送信信号ｙ（ｋ）とエコー低減信号ｅ
（ｋ）とからエコー信号ｂ（ｋ）の低減量を求め、受話
信号ｘ（ｋ）が存在し、かつ前記低減量が単調増大して
いる場合に送話信号ｓ（ｋ）が含まれていない状態であ
ると判断することを特徴とする請求項２記載のエコー抑
圧方法。 - 【請求項４】 エコー信号ｂ（ｋ）の短時間スペクトル
の振幅特性値と１単位時間前の該エコー信号ｂ（ｋ）の
短時間スペクトルの振幅特性値のうち、何れか大きい方
を当該エコー信号ｂ（ｋ）の短時間スペクトルの振幅特
性値の推定値とすることを特徴とする請求項１〜３いず
れかに記載のエコー抑圧方法。 - 【請求項５】 エコー信号ｂ（ｋ）の短時間スペクトル
の振幅特性値に聴感補正を施した値をエコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値の推定値とする
ことを特徴とする請求項１〜４いずれか記載のエコー抑
圧方法。 - 【請求項６】 受話音声によるエコー信号のマスキング
レベルと人間の可聴レベルと伝送遅延量に基づいて聴感
補正ゲインを設定し、該聴感補正ゲインをエコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値に乗算して聴感
補正することを特徴とする請求項５記載のエコー抑圧方
法。 - 【請求項７】 エコー信号ｂ（ｋ）の短時間スペクトル
の振幅特性値と１単位時間前の送話信号ｓ（ｋ）の短時
間スペクトルの振幅特性値と送信信号ｙ（ｋ）の短時間
スペクトル並びにスムージングファクタβと信号の正の
成分だけを取り出す半波整流関数Ｐからなる下式（１）
に基づいて送話信号ｓ（ｋ）の短時間スペクトルの振幅
特性値の推定値を算出することを特徴とする請求項１〜
６いずれかに記載のエコー抑圧方法。 【数１】 - 【請求項８】 スムージングファクタβを０．９〜０．
９９９の範囲内の値とすることを特徴とする請求項７記
載のエコー抑圧方法。 - 【請求項９】 送話信号パワー｜Ｓte（ω）｜2とエコ
ー信号パワー｜Ｅte（ω）｜2とからなる演算式（２）
の値ＧPE（ω）を、前記エコー抑圧ゲインＧ（ω）に適
用することを特徴とする請求項１〜８いずれかに記載の
エコー抑圧方法。 【数２】 - 【請求項１０】 送話信号パワー｜Ｓte（ω）｜2とエ
コー信号パワー｜Ｅte（ω）｜2とからなる演算式
（３）の値ＧW（ω）を、前記エコー抑圧ゲインＧ
（ω）に適用することを特徴とする請求項１〜８いずれ
かに記載のエコー抑圧方法。 【数３】 - 【請求項１１】 送話信号パワー｜Ｓte（ω）｜2とエコ
ー信号パワー｜Ｅte（ω）｜2とからなる演算式（４）
の値ＧML（ω）を、前記エコー抑圧ゲインＧ（ω）に適
用することを特徴とする請求項１〜８いずれかに記載の
エコー抑圧方法。 【数４】 - 【請求項１２】 送話信号パワー｜Ｓte（ω）｜2とエ
コー信号パワー｜Ｅte（ω）｜2と送信信号パワー｜Ｙt
（ω）｜2からなる式（５）に基づいて中間変数υ
（ω）を算出し、該中間変数υ（ω）とエコー信号パワ
ー｜Ｅte（ω）｜2と送信信号パワー｜Ｙt（ω）｜2並
びに中間変数υ（ω）に対する０次及び１次の修正ベッ
セル関数Ｉ0，Ｉ1からなる演算式（６）の値ＧEM（ω）
を、前記エコー抑圧ゲインＧ（ω）に適用することを特
徴とする請求項１〜８いずれかに記載のエコー抑圧方
法。 【数５】 【数６】 - 【請求項１３】 受話信号ｘ（ｋ）が時間変動するエコ
ー経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、エコー
信号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳された送信
信号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳された
エコー信号ｂ（ｋ）を抑圧するエコー抑圧装置であっ
て、 受話信号ｘ（ｋ）を短時間スペクトルに変換する受話信
号変換手段と、 送信信号ｙ（ｋ）を短時間スペクトルに変換する送信信
号変換手段と、 受話信号ｘ（ｋ）と送信信号ｙ（ｋ）とに基づいて送話
信号ｓ（ｋ）が含まれていない状態を検出する片側発話
状態検出部と、 送話信号ｓ（ｋ）が含まれない状態における送信信号ｙ
（ｋ）の短時間スペクトル及び前記受話信号ｘ（ｋ）の
短時間スペクトルを用いてエコー経路伝達関数ｈ（ｋ）
の振幅周波数関数を算出するエコー経路伝達関数算出手
段と、 該振幅周波数関数並びに送信信号ｙ（ｋ）の短時間スペ
クトル及び前記受話信号ｘ（ｋ）の短時間スペクトルを
用いて、以降の時刻におけるエコー信号ｂ（ｋ）の短時
間スペクトルの振幅特性値を推定するエコー信号振幅周
波数特性推定手段と、 該振幅周波数関数並びに送信信号ｙ（ｋ）の短時間スペ
クトル及び前記受話信号ｘ（ｋ）の短時間スペクトルを
用いて、以降の時刻における送話信号ｓ（ｋ）の短時間
スペクトルの振幅特性値を推定する送話信号振幅周波数
特性推定手段と、 該エコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
及び送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値
及び送信信号ｙ（ｋ）の短時間スペクトルの振幅特性値
に基づく演算式を用いてエコー抑圧ゲインＧ（ω）を算
出するエコー抑圧ゲイン決定部と、 該エコー抑圧ゲインＧ（ω）を送信信号ｙ（ｋ）の短時
間スペクトルに乗算するエコー抑圧処理部と、 該エコー抑圧処理部の出力を逆変換して伝送路に出力す
る送話信号逆変換手段と、 を具備することを特徴とするエコー抑圧装置。 - 【請求項１４】 受話信号ｘ（ｋ）が時間変動するエコ
ー経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、エコー
信号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳された送信
信号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳された
エコー信号ｂ（ｋ）を抑圧するエコー抑圧装置であっ
て、 受話信号ｘ（ｋ）を短時間スペクトルに変換する受話信
号変換手段と、 送信信号ｙ（ｋ）を短時間スペクトルに変換する送信信
号変換手段と、 受話信号ｘ（ｋ）と送信信号ｙ（ｋ）とに基づいて送話
信号ｓ（ｋ）が含まれていない状態を検出する片側発話
状態検出部と、 ダブルトーク中のエコー経路伝達関数ｈ（ｋ）の変動を
検出するエコー経路変動検出部と、 適応フィルタを用いてエコー経路伝達関数ｈ（ｋ）の推
定値である疑似エコー経路伝達関数ｈe（ｋ）を生成
し、該疑似エコー経路伝達関数ｈe（ｋ）を用いて送信
信号ｙ（ｋ）に含まれるエコー信号ｂ（ｋ）を抑圧して
送話信号ｓ（ｋ）に残留エコー信号ｂd（ｋ）が含まれ
たエコー低減信号ｅ（ｋ）を出力すると共に、前記エコ
ー経路変動検出部によってダブルトーク中のエコー経路
伝達関数ｈ（ｋ）の変動が検出されたときに、該変動検
出前の疑似エコー経路伝達関数ｈe（ｋ）を出力すると
共に送信信号ｙ（ｋ）をそのまま出力するエコーキャン
セラと、 前記エコー経路変動検出部によってダブルト
ーク中のエコー経路伝達関数ｈ（ｋ）の変動が検出され
ない場合に、送話信号ｓ（ｋ）が含まれない場合の前記
エコー低減信号ｅ（ｋ）の短時間スペクトルと受話信号
ｘ（ｋ）の短時間スペクトルとに基づいて残留振幅周波
数関数｜Ｈd（ω）｜2を算出するエコー経路伝達関数算
出手段と、 前記エコー経路変動検出部によってダブルトーク中のエ
コー経路伝達関数ｈ（ｋ）の変動が検出された場合に、
エコーキャンセラから出力された疑似エコー経路伝達関
数ｈe（ｋ）から疑似振幅周波数関数｜Ｈe（ω）｜2を
算出する第２のエコー経路伝達関数算出手段と、 エコー経路伝達関数算出手段あるいは第２のエコー経路
伝達関数算出手段から択一的に入力される前記疑似振幅
周波数関数｜Ｈe（ω）｜2あるいは残留振幅周波数関数
｜Ｈd（ω）｜2に基づいてエコー信号ｂ（ｋ）あるいは
残留エコー信号ｂd（ｋ）の短時間スペクトルの振幅周
波数特性を推定するエコー信号振幅周波数特性推定手段
と、 該振幅周波数特性並びにエコー低減信号ｅ（ｋ）あるい
は送信信号ｙ（ｋ）の短時間スペクトル及び前記受話信
号ｘ（ｋ）の短時間スペクトルを用いて、以降の時刻に
おける残留エコー信号ｂd（ｋ）あるいはエコー信号ｂ
（ｋ）の短時間スペクトルの振幅周波数特性を推定する
エコー信号振幅周波数特性推定手段と、 該振幅周波数特性並びにエコー低減信号ｅ（ｋ）あるい
は送信信号ｙ（ｋ）の短時間スペクトル及び前記受話信
号ｘ（ｋ）の短時間スペクトルを用いて、以降の時刻に
おける送話信号ｓ（ｋ）の短時間スペクトルの振幅周波
数特性を推定する送話信号振幅周波数特性推定手段と、 前記残留エコー信号ｂd（ｋ）あるいはエコー信号ｂ
（ｋ）の短時間スペクトルの振幅周波数特性及び送話信
号ｓ（ｋ）の短時間スペクトルの振幅周波数特性及びエ
コー低減信号ｅ（ｋ）あるいは送信信号ｙ（ｋ）の短時
間スペクトルの振幅周波数特性に基づく演算式を用いて
エコー抑圧ゲインＧ（ω）を算出するエコー抑圧ゲイン
決定部と、 該エコー抑圧ゲインＧ（ω）をエコー低減信号ｅ（ｋ）
あるいは送信信号ｙ（ｋ）の短時間スペクトルに乗算す
るエコー抑圧処理部と、 該エコー抑圧処理部の出力を逆変換して伝送路に出力す
る送話信号逆変換手段と、 を具備することを特徴とするエコー抑圧装置。 - 【請求項１５】 片側発話状態検出部は、送信信号ｙ
（ｋ）とエコー低減信号ｅ（ｋ）とからエコー信号ｂ
（ｋ）の低減量を求め、受話信号ｘ（ｋ）が存在し、か
つ前記低減量が単調増大している場合に送話信号ｓ
（ｋ）が含まれていない状態であると判断することを特
徴とする請求項１４記載のエコー抑圧装置。 - 【請求項１６】 エコー信号振幅周波数特性推定手段
は、 振幅周波数関数並びに送信信号ｙ（ｋ）の短時間スペク
トル及び前記受話信号ｘ（ｋ）の短時間スペクトルを用
いて、以降の時刻におけるエコー信号ｂ（ｋ）の短時間
スペクトルの振幅特性値を推定するエコー信号パワー計
算部と、 この振幅特性値を１単位時間遅延させて出力する第１の
遅延演算子と、 エコー信号パワー計算部と第１の遅延演算子の出力のう
ち、大きい方を当該エコー信号ｂ（ｋ）の短時間スペク
トルの振幅特性値の推定値とするエコー信号パワー決定
部とからなることを特徴とする請求項１３〜１５いずれ
かに記載のエコー抑圧装置。 - 【請求項１７】 エコー信号振幅周波数特性推定手段
は、エコー信号ｂ（ｋ）の短時間スペクトルの振幅特性
値に聴感補正を施すエコー信号聴感補正手段を備えるこ
とを特徴とする請求項１３〜１６いずれかに記載のエコ
ー抑圧装置。 - 【請求項１８】 エコー信号聴感補正手段は、 受話音声によるエコー信号のマスキングレベルと人間の
可聴レベルと伝送遅延量に基づいて設定された聴感補正
ゲインを記憶する聴感補正ゲイン記憶部と、 前記聴感補正ゲインをエコー信号ｂ（ｋ）の短時間スペ
クトルの振幅特性値に乗算する聴感補正部とからなるこ
とを特徴とする請求項１７記載のエコー抑圧装置。 - 【請求項１９】 送話信号振幅周波数特性推定手段は、 スムージングファクタβを記憶するスムージングファク
タ記憶部と、 送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値を１
単位時間遅延させる第２の遅延演算子と、 該第２の遅延演算子の出力とエコー信号ｂ（ｋ）の短時
間スペクトルの振幅特性値と送信信号ｙ（ｋ）の短時間
スペクトルとスムージングファクタβと信号の正の成分
だけを取り出す半波整流関数Ｐからなる下式（１）に基
づいて送話信号ｓ（ｋ）の短時間スペクトルの振幅特性
値の推定値を算出する送信信号パワー計算部とからなる
ことを特徴とする請求項１３〜１８いずれかに記載のエ
コー抑圧装置。 【数７】 - 【請求項２０】 スムージングファクタβを０．９〜
０．９９９の範囲内の値とすることを特徴とする請求項
１９記載のエコー抑圧装置。 - 【請求項２１】 エコー抑圧ゲイン決定部は、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（２）に基づいて値ＧPE（ω）を算出することを特徴
とする請求項１３〜１９いずれかに記載のエコー抑圧装
置。 【数８】 - 【請求項２２】 エコー抑圧ゲイン決定部は、エコー抑
圧ゲインＧ（ω）として、送話信号パワー｜Ｓte（ω）
｜2とエコー信号パワー｜Ｅte（ω）｜2とからなる演算
式（３）に基づいて値ＧW（ω）を算出することを特徴
とする請求項１３〜１９いずれかに記載のエコー抑圧装
置。 【数９】 - 【請求項２３】 エコー抑圧ゲイン決定部は、エコー
抑圧ゲインＧ（ω）として、送話信号パワー｜Ｓte
（ω）｜2とエコー信号パワー｜Ｅte（ω）｜2とからな
る演算式（４）に基づいて値ＧML（ω）を算出すること
を特徴とする請求項１３〜１９いずれかに記載のエコー
抑圧装置。 【数１０】 - 【請求項２４】 エコー抑圧ゲイン決定部は、送話信号
パワー｜Ｓte（ω）｜2とエコー信号パワー｜Ｅte
（ω）｜2と送信信号パワー｜Ｙt（ω）｜2からなる式
（５）に基づいて中間変数υ（ω）を算出し、該中間変
数υ（ω）とエコー信号パワー｜Ｅte（ω）｜2と送信
信号パワー｜Ｙt（ω）｜2並びに中間変数υ（ω）に対
する０次及び１次の修正ベッセル関数Ｉ0，Ｉ1からなる
演算式（６）に基づいて値ＧEM（ω）をエコー抑圧ゲイ
ンＧ（ω）として算出することを特徴とする請求項１３
〜１９いずれかに記載のエコー抑圧装置。 【数１１】 【数１２】 - 【請求項２５】 受話信号ｘ（ｋ）が時間変動するエコ
ー経路伝達関数ｈ（ｋ）のエコー経路を伝搬し、エコー
信号ｂ（ｋ）として送話信号ｓ（ｋ）に重畳された送信
信号ｙ（ｋ）に対し、該送信信号ｙ（ｋ）に重畳された
エコー信号ｂ（ｋ）を抑圧するエコー抑圧プログラムが
記憶されたコンピュータに読取り可能な記憶媒体であっ
て、 受話信号ｘ（ｋ）及び送信信号ｙ（ｋ）を短時間スペク
トルにそれぞれ変換する処理と、 送話信号ｓ（ｋ）が含まれない状態における送信信号ｙ
（ｋ）の短時間スペクトル及び前記受話信号ｘ（ｋ）の
短時間スペクトルを用いてエコー経路伝達関数ｈ（ｋ）
の振幅周波数関数を算出する処理と、 該振幅周波数関数並びに送信信号ｙ（ｋ）の短時間スペ
クトル及び前記受話信号ｘ（ｋ）の短時間スペクトルを
用いて、以降の時刻におけるエコー信号ｂ（ｋ）の短時
間スペクトルの振幅特性値及び送話信号ｓ（ｋ）の短時
間スペクトルの振幅特性値を推定する処理と、 該エコー信号ｂ（ｋ）の短時間スペクトルの振幅特性値
及び送話信号ｓ（ｋ）の短時間スペクトルの振幅特性値
及び送信信号ｙ（ｋ）の短時間スペクトルの振幅特性値
に基づく演算式を用いてエコー抑圧ゲインＧ（ω）を算
出する処理と、 該エコー抑圧ゲインＧ（ω）を送信信号ｙ（ｋ）の短時
間スペクトルに乗算することによりエコー信号ｂ（ｋ）
を抑圧するする処理と、 からなることを特徴とするエコー抑圧プログラムが記憶
されたコンピュータに読取り可能な記憶媒体。 - 【請求項２６】 エコー信号ｂ（ｋ）の短時間スペクト
ルの振幅特性値と１単位時間前の該エコー信号ｂ（ｋ）
の短時間スペクトルの振幅特性値のうち、何れか大きい
方を当該エコー信号ｂ（ｋ）の短時間スペクトルの振幅
特性値の推定値とする処理を含むことを特徴とする請求
項２５記載のエコー抑圧プログラムが記憶されたコンピ
ュータに読取り可能な記憶媒体。 - 【請求項２７】 エコー信号ｂ（ｋ）の短時間スペクト
ルの振幅特性値に聴感補正を施した値をエコー信号ｂ
（ｋ）の短時間スペクトルの振幅特性値の推定値とする
処理を含むことを特徴とする請求項２５または２６記載
のエコー抑圧プログラムが記憶されたコンピュータに読
取り可能な記憶媒体。 - 【請求項２８】 受話音声によるエコー信号のマスキン
グレベルと人間の可聴レベルと伝送遅延量に基づいて聴
感補正ゲインを設定し、該聴感補正ゲインをエコー信号
ｂ（ｋ）の短時間スペクトルの振幅特性値に乗算して聴
感補正する処理を含むことを特徴とする請求項２７記載
のエコー抑圧プログラムが記憶されたコンピュータに読
取り可能な記憶媒体。 - 【請求項２９】 エコー信号ｂ（ｋ）の短時間スペクト
ルの振幅特性値と１単位時間前の送話信号ｓ（ｋ）の短
時間スペクトルの振幅特性値と送信信号ｙ（ｋ）の短時
間スペクトル並びにスムージングファクタβと信号の正
の成分だけを取り出す半波整流関数Ｐからなる下式
（１）に基づいて送話信号ｓ（ｋ）の短時間スペクトル
の振幅特性値の推定値を算出する処理を含むことを特徴
とする請求項２５〜２８いずれかに記載のエコー抑圧プ
ログラムが記憶されたコンピュータに読取り可能な記憶
媒体。 【数１３】 - 【請求項３０】 スムージングファクタβを０．９〜
０．９９９の範囲内の値とする処理を含むことを特徴と
する請求項２９記載のエコー抑圧プログラムが記憶され
たコンピュータに読取り可能な記憶媒体。 - 【請求項３１】 送話信号パワー｜Ｓte（ω）｜2とエ
コー信号パワー｜Ｅte（ω）｜2とからなる演算式
（２）の値ＧPE（ω）を、前記エコー抑圧ゲインＧ
（ω）に適用する処理を含むことを特徴とする請求項２
５〜３０いずれかに記載のエコー抑圧プログラムが記憶
されたコンピュータに読取り可能な記憶媒体。 【数１４】 - 【請求項３２】 送話信号パワー｜Ｓte（ω）｜2とエ
コー信号パワー｜Ｅte（ω）｜2とからなる演算式
（３）の値ＧW（ω）を、前記エコー抑圧ゲインＧ
（ω）に適用する処理を含むことを特徴とする請求項２
５〜３０いずれかに記載のエコー抑圧プログラムが記憶
されたコンピュータに読取り可能な記憶媒体。 【数１５】 - 【請求項３３】 送話信号パワー｜Ｓte（ω）｜2とエ
コー信号パワー｜Ｅte（ω）｜2とからなる演算式
（４）の値ＧML（ω）を、前記エコー抑圧ゲインＧ
（ω）に適用する処理を含むことを特徴とする請求項２
５〜３０いずれかに記載のエコー抑圧プログラムが記憶
されたコンピュータに読取り可能な記憶媒体。 【数１６】 - 【請求項３４】 送話信号パワー｜Ｓte（ω）｜2とエ
コー信号パワー｜Ｅte（ω）｜2と送信信号パワー｜Ｙt
（ω）｜2からなる式（５）に基づいて中間変数υ
（ω）を算出し、該中間変数υ（ω）とエコー信号パワ
ー｜Ｅte（ω）｜2と送信信号パワー｜Ｙt（ω）｜2並
びに中間変数υ（ω）に対する０次及び１次の修正ベッ
セル関数Ｉ0，Ｉ1からなる演算式（６）の値ＧEM（ω）
を、前記エコー抑圧ゲインＧ（ω）に適用する処理を含
むことを特徴とする請求項２５〜３０いずれかに記載の
エコー抑圧プログラムが記憶されたコンピュータに読取
り可能な記憶媒体。 【数１７】 【数１８】
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25441198A JP3420705B2 (ja) | 1998-03-16 | 1998-09-08 | エコー抑圧方法及び装置並びにエコー抑圧プログラムが記憶されたコンピュータに読取り可能な記憶媒体 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10-65923 | 1998-03-16 | ||
| JP6592398 | 1998-03-16 | ||
| JP25441198A JP3420705B2 (ja) | 1998-03-16 | 1998-09-08 | エコー抑圧方法及び装置並びにエコー抑圧プログラムが記憶されたコンピュータに読取り可能な記憶媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11331046A JPH11331046A (ja) | 1999-11-30 |
| JP3420705B2 true JP3420705B2 (ja) | 2003-06-30 |
Family
ID=26407092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25441198A Expired - Lifetime JP3420705B2 (ja) | 1998-03-16 | 1998-09-08 | エコー抑圧方法及び装置並びにエコー抑圧プログラムが記憶されたコンピュータに読取り可能な記憶媒体 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3420705B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009148049A1 (ja) | 2008-06-02 | 2009-12-10 | ヤマハ株式会社 | 音響エコーキャンセラ及び音響エコーキャンセル方法 |
| US8335311B2 (en) | 2005-07-28 | 2012-12-18 | Kabushiki Kaisha Toshiba | Communication apparatus capable of echo cancellation |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1684442A4 (en) * | 2003-11-11 | 2009-04-15 | Mitsubishi Electric Corp | ECHO SUPPRESSOR |
| JP4448424B2 (ja) * | 2004-10-25 | 2010-04-07 | 日本電信電話株式会社 | 音声スイッチ方法、この方法を実施する装置、プログラムおよびその記録媒体 |
| JP4448423B2 (ja) * | 2004-10-25 | 2010-04-07 | 日本電信電話株式会社 | エコー抑圧方法、この方法を実施する装置、プログラムおよびその記録媒体 |
| EP1853087B1 (en) | 2005-02-21 | 2014-07-30 | Fujitsu Ltd. | Echo canceller |
| JP4818014B2 (ja) * | 2005-07-28 | 2011-11-16 | 株式会社東芝 | 信号処理装置 |
| JP4916394B2 (ja) | 2007-07-03 | 2012-04-11 | 富士通株式会社 | エコー抑圧装置、エコー抑圧方法及びコンピュータプログラム |
| JP4978352B2 (ja) * | 2007-07-11 | 2012-07-18 | ヤマハ株式会社 | エコーキャンセラ |
| JP5071480B2 (ja) | 2007-08-24 | 2012-11-14 | 富士通株式会社 | エコー抑制装置、エコー抑制システム、エコー抑制方法及びコンピュータプログラム |
| WO2009047858A1 (ja) | 2007-10-12 | 2009-04-16 | Fujitsu Limited | エコー抑圧システム、エコー抑圧方法、エコー抑圧プログラム、エコー抑圧装置、音出力装置、オーディオシステム、ナビゲーションシステム及び移動体 |
| JP2009130377A (ja) * | 2007-11-19 | 2009-06-11 | Oki Semiconductor Co Ltd | 受信信号調整装置、及びエコーキャンセラ |
| CN102165707B (zh) * | 2008-09-24 | 2013-12-04 | 三菱电机株式会社 | 回波消除装置 |
| JP5087024B2 (ja) * | 2009-02-10 | 2012-11-28 | 日本電信電話株式会社 | エコー消去装置とその方法と、プログラム |
| JP5466581B2 (ja) * | 2010-06-04 | 2014-04-09 | 日本電信電話株式会社 | 反響消去方法、反響消去装置及び反響消去プログラム |
| JP5562451B1 (ja) * | 2013-01-31 | 2014-07-30 | 日本電信電話株式会社 | エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置とプログラム |
| JP6111795B2 (ja) * | 2013-03-28 | 2017-04-12 | 富士通株式会社 | 信号処理装置、及び信号処理方法 |
| JP6363429B2 (ja) * | 2014-08-25 | 2018-07-25 | 日本電信電話株式会社 | データ構造、データ生成装置、データ生成方法、およびプログラム |
| CN105280195B (zh) * | 2015-11-04 | 2018-12-28 | 腾讯科技（深圳）有限公司 | 语音信号的处理方法及装置 |
| JP6878776B2 (ja) | 2016-05-30 | 2021-06-02 | 富士通株式会社 | 雑音抑圧装置、雑音抑圧方法及び雑音抑圧用コンピュータプログラム |
| CN114859180B (zh) * | 2022-05-30 | 2024-06-25 | 海南电网有限责任公司乐东供电局 | 一种基于连续波的电缆故障双端定位方法 |
-
1998
- 1998-09-08 JP JP25441198A patent/JP3420705B2/ja not_active Expired - Lifetime
Non-Patent Citations (2)
| Title |
|---|
| 阪内 澄宇，佐々木 潤子，羽田 陽一，"ＳＴＳＡ推定に基づいた非線形エコー抑圧と適応フィルタの組合せ"，日本音響学会 平成１０年度春秋季研究発表会 講演論文集 Ｉ，１９９８年 ９月２４日，ｐｐ．６０７−６０８，（１−Ｐ−６） |
| 阪内 澄宇，羽田 陽一，"短時間スペクトラル振幅推定に基づいた非線形エコー抑圧処理の検討"，日本音響学会 平成１０年度春季研究発表会 講演論文集Ｉ，１９９８年 ３月１７日，ｐｐ．５５１−５５２，（３−５−１１） |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8335311B2 (en) | 2005-07-28 | 2012-12-18 | Kabushiki Kaisha Toshiba | Communication apparatus capable of echo cancellation |
| WO2009148049A1 (ja) | 2008-06-02 | 2009-12-10 | ヤマハ株式会社 | 音響エコーキャンセラ及び音響エコーキャンセル方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH11331046A (ja) | 1999-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3420705B2 (ja) | エコー抑圧方法及び装置並びにエコー抑圧プログラムが記憶されたコンピュータに読取り可能な記憶媒体 | |
| US8311234B2 (en) | Echo canceller and communication audio processing apparatus | |
| CN110838300B (zh) | 回声消除的处理方法及处理系统 | |
| EP1698159B1 (en) | System and method for enhanced stereo audio | |
| EP3791565B1 (en) | Method and apparatus utilizing residual echo estimate information to derive secondary echo reduction parameters | |
| US8160239B2 (en) | Echo canceller and speech processing apparatus | |
| US20130301840A1 (en) | Methods for processing audio signals and circuit arrangements therefor | |
| JP4568439B2 (ja) | エコー抑圧装置 | |
| JP3507020B2 (ja) | 反響抑圧方法、反響抑圧装置及び反響抑圧プログラム記憶媒体 | |
| JP3607625B2 (ja) | 多チャネル反響抑圧方法、その装置、そのプログラム及びその記録媒体 | |
| US7403608B2 (en) | Echo processing devices for single-channel or multichannel communication systems | |
| US6834108B1 (en) | Method for improving acoustic noise attenuation in hand-free devices | |
| JP3753996B2 (ja) | 反響抑圧装置、反響抑圧方法及びプログラム | |
| JP3407392B2 (ja) | ステレオエコーキャンセラ | |
| US8369511B2 (en) | Robust method of echo suppressor | |
| JP2001044896A (ja) | 通話装置および通話方法 | |
| JP2002223182A (ja) | 反響消去方法、その装置、そのプログラム及びその記録媒体 | |
| JPH09116613A (ja) | 反響消去装置 | |
| JP3220979B2 (ja) | 音声スイッチ | |
| JP2007504687A (ja) | ネットワークにおける通信リンクの送話品質の測定 | |
| US7023986B2 (en) | Echo canceller in a communication system at a terminal | |
| JP4877083B2 (ja) | 残留エコー抑圧制御装置、方法及びプログラム | |
| US7711107B1 (en) | Perceptual masking of residual echo | |
| US20020012429A1 (en) | Interference-signal-dependent adaptive echo suppression | |
| JP3355594B2 (ja) | エコーキャンセラ装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090418 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090418 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100418 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100418 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110418 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120418 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130418 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140418 Year of fee payment: 11 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |