JP5325134B2 - 反響消去方法、反響消去装置、そのプログラムおよび記録媒体 - Google Patents
反響消去方法、反響消去装置、そのプログラムおよび記録媒体 Download PDFInfo
- Publication number
- JP5325134B2 JP5325134B2 JP2010023554A JP2010023554A JP5325134B2 JP 5325134 B2 JP5325134 B2 JP 5325134B2 JP 2010023554 A JP2010023554 A JP 2010023554A JP 2010023554 A JP2010023554 A JP 2010023554A JP 5325134 B2 JP5325134 B2 JP 5325134B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- frequency component
- echo
- signal
- coefficient
- 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 - Fee Related
Links
Landscapes
- Circuit For Audible Band Transducer (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Telephone Function (AREA)
Description
次に、非特許文献2のエコーキャンセラの動作原理について説明する。非特許文献2のエコーキャンセラでは、線形の適応フィルタ
w(i)=0.54−0.46cos(2πn/L)
を用いて、
〔参考文献1〕E.R.Ferrara,"Fast Implimentation of LMS Adaptive Filters", IEEE Transactions on Acoustics, Speech, and Signal Processing, 1980, Vol.28, No.4, p.474-475
ブロック抽出部270は、周波数逆変換部260で得られた時間領域の疑似エコーの要素のうち後半のL個
差信号生成部250は周波数成分ごとに合計K個設けられ、ブロック抽出部270で抽出したL個の要素を、収音信号
〔参考文献2〕牧野昭二、「サブバンド信号処理 −実時間動作化の奥の手−」、日本音響学会誌、2000年、56巻12号、p.845-851
乗算係数計算部420は、帯域成分xk(r)ごとにMk個の乗算係数βk,m(m=1,2,・・・, Mk)を、各帯域成分の大きさに応じて生成する(S22)。具体的には以下のように生成する。
差信号生成部450は帯域成分ごとに合計K個設けられ、それぞれ収音信号の帯域成分yk(r)から疑似エコーy'k(r)を差し引いて差信号ek(r)を生成する(S25)。
ここで、μは更新量を制御するステップサイズ(0<μ<2)でり、δは微小な正の定数である。
本発明が有効である理由の説明のため、以下に実験結果を示す。
まず、非線形歪がインパルス応答に及ぼす影響を調べるために、振幅の異なる複数のM系列を用いてインパルス応答の測定を行った。図11にM系列を用いて測定した5つのインパルス応答を示す。振幅の小さいM系列で測定したものから順にインパルス応答1〜5と表記する。各インパルス応答は同じ概形であるが正確な値は異なっていることがわかる。
Claims (7)
- 時間領域の受話信号
(nは時刻を表す整数、Lはフレーム長)をK個の周波数成分X(k) (k=1,2,・・・,K、Kは2以上の整数)に変換する第1周波数変換ステップと、
上記周波数成分X(k)ごとにMk個(Mkはkの値に応じた正の整数)の乗算係数βk,m(m=1,2,・・・, Mk)を、各周波数成分の大きさに応じて生成する乗算係数生成ステップと、
上記周波数成分X(k)ごとにそれぞれ用意されたMk個の適応フィルタにおいて、当該適応フィルタのフィルタ係数Hk,mと上記乗算係数βk,mと上記周波数成分X(k)とを用いて、適応フィルタごとの疑似エコーY'(k,1), Y'(k,2),・・・, Y'(k, Mk)を求め、これらを合算して周波数成分X(k)に対するトータルでの疑似エコーY'(k)を生成する疑似エコー生成ステップと、
時間領域の収音信号
をK個の周波数成分Y(k)に変換する第2周波数変換ステップと、
上記周波数成分Y(k)から上記疑似エコーY'(k)を差し引いて差信号E(k)を生成する差信号生成ステップと、
上記差信号E(1), E(2),・・・,E(L/2+1),E*(L/2),E*(L/2−1),・・・,E*(2)(*は複素共役を表す)を合成し時間領域の送話信号
に変換する周波数逆変換ステップと、
上記適応フィルタのフィルタ係数Hk,mを、上記差信号E(k)と上記周波数成分X(k)と上記乗算係数βk,mとを用いて更新するフィルタ係数更新ステップと、
を実行する反響消去方法。 - 時間領域の受話信号
(nは時刻を表す整数、Lはフレーム長)をK個の周波数成分X(k) (k=1,2,・・・,K、Kは2以上の整数)に変換する第1周波数変換ステップと、
上記周波数成分X(k)ごとにMk個(Mkはkの値に応じた正の整数)の乗算係数βk,m(m=1,2,・・・, Mk)を、各周波数成分の大きさに応じて生成する乗算係数生成ステップと、
上記周波数成分X(k)ごとにそれぞれ用意されたMk個の適応フィルタにおいて、当該適応フィルタのフィルタ係数Hk,mと上記乗算係数βk,mと上記周波数成分X(k)とを用いて、適応フィルタごとの疑似エコーY'(k,1), Y'(k,2),・・・, Y'(k, Mk)を求め、これらを合算して周波数成分X(k)に対するトータルでの疑似エコーY'(k)を生成する疑似エコー生成ステップと、
K個の上記疑似エコーY'(1), Y'(2),・・・, Y'(K)を合成し時間領域の疑似エコー
に変換する周波数逆変換ステップと、
上記周波数逆変換ステップで得られた時間領域の疑似エコーの要素のうち後半のL個
を抽出するブロック抽出ステップと、
上記ブロック抽出ステップで抽出したL個の要素を、収音信号
の対応する各要素からそれぞれ差し引くことにより、送話信号
を生成する送話信号生成ステップと、
上記送話信号の先頭にL個のゼロを付加し、ゼロ付加送話信号
を生成するゼロ付加ステップと、
上記ゼロ付加送話信号を、K個の周波数成分E'(k)に変換する第2周波数変換ステップと、
上記周波数成分X(k)と上記ゼロ付加送話信号の周波数成分E'(k)とから更新係数S(k)を計算する更新係数計算ステップと、
上記適応フィルタのフィルタ係数Hk,mを、上記更新係数S(k)と上記周波数成分X(k)と上記乗算係数βk,mとを用いて更新するフィルタ係数更新ステップと、
を実行する反響消去方法。 - 請求項1又は2に記載の反響消去方法であって、
上記乗算係数を、複数の周波数成分X(k)の大きさをまとめて評価することにより生成することを特徴とする反響消去方法。 - 時間領域の受話信号
(nは時刻を表す整数、Lはフレーム長)をK個の周波数成分X(k) (k=1,2,・・・,K、Kは2以上の整数)に変換する第1周波数変換部と、
上記周波数成分X(k)ごとにM k 個(M k はkの値に応じた正の整数)の乗算係数β k,m (m=1,2,・・・, M k )を、各周波数成分の大きさに応じて生成する乗算係数計算部と、
上記周波数成分X(k)ごとにそれぞれ用意されたM k 個の適応フィルタにおいて、当該適応フィルタのフィルタ係数H k,m と上記乗算係数β k,m と上記周波数成分X(k)とを用いて、適応フィルタごとの疑似エコーY'(k,1), Y'(k,2),・・・, Y'(k, M k )を求め、これらを合算して周波数成分X(k)に対するトータルでの疑似エコーY'(k)を生成する複合適応フィルタ部と、
時間領域の収音信号
をK個の周波数成分Y(k)に変換する第2周波数変換部と、
上記周波数成分Y(k)から上記疑似エコーY'(k)を差し引いて差信号E(k)を生成する差信号生成部と、
上記差信号E(1), E(2),・・・,E(L/2+1),E * (L/2),E * (L/2−1),・・・,E * (2)(*は複素共役を表す)を合成し時間領域の送話信号
に変換する周波数逆変換部と、
を備え、
上記複合適応フィルタ部は、上記適応フィルタのフィルタ係数H k,m を、上記差信号E(k)と上記周波数成分X(k)と上記乗算係数β k,m とを用いて更新する
反響消去装置。 - 時間領域の受話信号
(nは時刻を表す整数、Lはフレーム長)をK個の周波数成分X(k) (k=1,2,・・・,K、Kは2以上の整数)に変換する第1周波数変換部と、
上記周波数成分X(k)ごとにM k 個(M k はkの値に応じた正の整数)の乗算係数β k,m (m=1,2,・・・, M k )を、各周波数成分の大きさに応じて生成する乗算係数計算部と、
上記周波数成分X(k)ごとにそれぞれ用意されたM k 個の適応フィルタにおいて、当該適応フィルタのフィルタ係数H k,m と上記乗算係数β k,m と上記周波数成分X(k)とを用いて、適応フィルタごとの疑似エコーY'(k,1), Y'(k,2),・・・, Y'(k, M k )を求め、これらを合算して周波数成分X(k)に対するトータルでの疑似エコーY'(k)を生成する複合適応フィルタ部と、
K個の上記疑似エコーY'(1), Y'(2),・・・, Y'(K)を合成し時間領域の疑似エコー
に変換する周波数逆変換部と、
上記周波数逆変換部で得られた時間領域の疑似エコーの要素のうち後半のL個
を抽出するブロック抽出部と、
上記ブロック抽出部で抽出したL個の要素を、収音信号
の対応する各要素からそれぞれ差し引くことにより、送話信号
を生成する送話信号生成部と、
上記送話信号の先頭にL個のゼロを付加し、ゼロ付加送話信号
を生成するゼロ付加部と、
上記ゼロ付加送話信号を、K個の周波数成分E'(k)に変換する第2周波数変換部と、
上記周波数成分X(k)と上記ゼロ付加送話信号の周波数成分E'(k)とから更新係数S(k)を計算する更新係数計算部と、
を備え、
上記複合適応フィルタ部は、上記適応フィルタのフィルタ係数H k,m を、上記更新係数S(k)と上記周波数成分X(k)と上記乗算係数β k,m とを用いて更新する
反響消去装置。 - 請求項1乃至3のいずれかに記載の反響消去方法をコンピュータに実行させるためのプログラム。
- 請求項1乃至3のいずれかに記載の反響消去方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010023554A JP5325134B2 (ja) | 2010-02-04 | 2010-02-04 | 反響消去方法、反響消去装置、そのプログラムおよび記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010023554A JP5325134B2 (ja) | 2010-02-04 | 2010-02-04 | 反響消去方法、反響消去装置、そのプログラムおよび記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011166239A JP2011166239A (ja) | 2011-08-25 |
JP5325134B2 true JP5325134B2 (ja) | 2013-10-23 |
Family
ID=44596454
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010023554A Expired - Fee Related JP5325134B2 (ja) | 2010-02-04 | 2010-02-04 | 反響消去方法、反響消去装置、そのプログラムおよび記録媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5325134B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3667662B1 (en) * | 2018-12-12 | 2022-08-10 | Panasonic Intellectual Property Corporation of America | Acoustic echo cancellation device, acoustic echo cancellation method and acoustic echo cancellation program |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03243020A (ja) * | 1990-02-21 | 1991-10-30 | Fujitsu Ltd | サブバンド型エコーキャンセラ |
JP2000353989A (ja) * | 1999-06-11 | 2000-12-19 | Matsushita Electric Ind Co Ltd | エコーキャンセラ |
US9318119B2 (en) * | 2005-09-02 | 2016-04-19 | Nec Corporation | Noise suppression using integrated frequency-domain signals |
US8275120B2 (en) * | 2006-05-30 | 2012-09-25 | Microsoft Corp. | Adaptive acoustic echo cancellation |
JP4413205B2 (ja) * | 2006-06-21 | 2010-02-10 | 日本電信電話株式会社 | エコー抑圧方法、装置、エコー抑圧プログラム、記録媒体 |
JP5016500B2 (ja) * | 2007-01-17 | 2012-09-05 | ティーオーエー株式会社 | 未知系同定システム |
-
2010
- 2010-02-04 JP JP2010023554A patent/JP5325134B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2011166239A (ja) | 2011-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5075042B2 (ja) | エコー消去装置、エコー消去方法、そのプログラム、記録媒体 | |
JP4787851B2 (ja) | エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置と、装置プログラムと記録媒体 | |
JP2019078864A (ja) | 楽音強調装置、畳み込みオートエンコーダ学習装置、楽音強調方法、プログラム | |
JP5161157B2 (ja) | 周波数領域エコー除去装置、周波数領域エコー除去方法、プログラム | |
JP5769671B2 (ja) | エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置とプログラム | |
JP5889224B2 (ja) | エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置とプログラム | |
JP5325134B2 (ja) | 反響消去方法、反響消去装置、そのプログラムおよび記録媒体 | |
JP5769670B2 (ja) | エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置とプログラム | |
JP5769672B2 (ja) | エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置とプログラム | |
JP5524316B2 (ja) | パラメータ推定装置、エコー消去装置、パラメータ推定方法、及びプログラム | |
JP5583181B2 (ja) | 縦続接続型伝達系パラメータ推定方法、縦続接続型伝達系パラメータ推定装置、プログラム | |
JP5438629B2 (ja) | ステレオ反響消去方法、ステレオ反響消去装置、ステレオ反響消去プログラム | |
JP5044594B2 (ja) | 多チャネルエコー消去装置とその方法、そのプログラム | |
JP5562451B1 (ja) | エコー抑圧ゲイン推定方法とそれを用いたエコー消去装置とプログラム | |
JP5826712B2 (ja) | マルチチャネルエコー消去装置、マルチチャネルエコー消去方法、およびプログラム | |
JP2004274412A (ja) | 反響消去装置、方法、及び反響消去プログラム | |
JP6827908B2 (ja) | 音源強調装置、音源強調学習装置、音源強調方法、プログラム | |
JP5925149B2 (ja) | 音響結合量推定装置、エコー消去装置、その方法及びプログラム | |
JP5526053B2 (ja) | 反響消去方法とその装置とプログラム | |
JP7235117B2 (ja) | エコー消去装置、エコー消去方法、プログラム | |
JP5086969B2 (ja) | エコーキャンセル装置、その方法、そのプログラム、およびその記録媒体 | |
JP5097148B2 (ja) | 音響結合量算出装置とその方法と、プログラム | |
JP5596724B2 (ja) | 縦続接続型伝達系パラメータ推定方法、縦続接続型伝達系パラメータ推定装置、プログラム | |
WO2021100094A1 (ja) | 音源信号推定装置、音源信号推定方法、プログラム | |
JP5889233B2 (ja) | 音響結合量推定装置、エコー消去装置、その方法及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20110624 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120116 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130423 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130613 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20130709 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130719 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5325134 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
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 |
|
LAPS | Cancellation because of no payment of annual fees |