JP4156428B2 - Echo canceling method, echo canceling device, echo canceling program - Google Patents

Echo canceling method, echo canceling device, echo canceling program Download PDF

Info

Publication number
JP4156428B2
JP4156428B2 JP2003108875A JP2003108875A JP4156428B2 JP 4156428 B2 JP4156428 B2 JP 4156428B2 JP 2003108875 A JP2003108875 A JP 2003108875A JP 2003108875 A JP2003108875 A JP 2003108875A JP 4156428 B2 JP4156428 B2 JP 4156428B2
Authority
JP
Japan
Prior art keywords
filter
error
signal
simulated
echo
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
Application number
JP2003108875A
Other languages
Japanese (ja)
Other versions
JP2004320204A (en
Inventor
末廣 島内
陽一 羽田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2003108875A priority Critical patent/JP4156428B2/en
Publication of JP2004320204A publication Critical patent/JP2004320204A/en
Application granted granted Critical
Publication of JP4156428B2 publication Critical patent/JP4156428B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Telephone Function (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、スピーカのような音響再生手段からマイクロホンのような音響収音手段へと回り込む反響を消去するための反響消去装置に関する。
【0002】
【従来の技術】
音響再生手段から音響収音手段へ回り込む反響を消去する反響消去装置は、図6のように接続される。図6に示す10は従来の反響消去装置を示す。従来の反響消去装置10内では、音響再生手段1と音響収音手段2間の反響路のインパルス応答hを推定し、推定したインパルス応答h´と再生信号入力端子3に入力された再生信号xの畳込み演算h´*xを生成し、実際の反響信号yから減算することで、反響消去信号出力端子4に反響消去信号eを得る。
しかし、推定したインパルス応答と再生信号畳込み演算には、多くの演算量を必要とし、実装上の問題となっている。
【0003】
近年、この問題を解決するため、再生信号や反響信号を一旦周波数領域変換し、反響路のインパルス応答の周波数領域変換に対応したパラメータを推定し、畳込みの代わりに乗算処理を用いたり(非特許文献1)、あるいはより小さい畳込み演算に分割したりする(特許文献1)などして、演算量を削減する方法が提案されている。周波数領域変換の例としては、(高速)離散フーリエ変換、(高速)離散コサイン変換、(高速)ハートレー変換などがある。
図6で音響再生手段1としてスピーカを挙げているが、音響再生手段1としては、再生前段の増幅器やバッファも含む。また、同様に音響収音手段2は、マイクロホンの後段の増幅器やバッファも含む。
【0004】
また、(非特許文献1)などに記載されている構成を模式的に図7に示す。図7の構成は以下の手段を含んでいる。すなわち、再生信号入力端子3から入力されて音響再生手段1へと出力される再生信号を入力し一定時間蓄積し再生信号列を得る再生信号入力手段101と、音響再生手段1と同一空間に存在する音響収音手段2から収音信号を入力し一定時間蓄積し収音信号列を得る収音信号入力手段102と、再生信号列を周波数領域変換し再生信号変換列を得る再生信号変換手段103と、再生信号変換列を周波数領域毎に入力しフィルタ処理することにより、音響再生手段1から音響収音手段2へと回り込む反響信号の周波数領域変換列を模擬する模擬反響信号変換列を生成するフィルタ処理手段104と、模擬反響信号変換列と収音信号列とを入力し模擬反響信号変換列の模擬誤差を出力する模擬誤差出力手段105と、再生信号変換列と模擬誤差を入力し、フィルタ処理手段104のフィルタ係数の周波数領域毎の特性誤差を計算するフィルタ誤差計算手段106と、フィルタ誤差計算手段106が計算した特性誤差に周波数領域毎に第一調整係数を乗じる第一更新量調整手段107と、第一調整係数を乗じられた特性誤差を加えることによりフィルタ処理手段104のフィルタ係数を更新するフィルタ更新手段108とを含む。
【0005】
【特許文献1】
特開平9−116472号(図6)
【非特許文献1】
Simon Haykin 著 「適応フィルタ理論」 科学技術出版、2001年1月10日、pp.500−541
【0006】
【発明が解決しようとする課題】
図7の構成において、フィルタ誤差計算手段106が、反響消去のための真の特性とフィルタ処理手段104で用いるフィルタ係数との特性誤差を正しく計算し出力できれば、フィルタ誤差計算手段106から出力される特性誤差に対して、第一更新量調整手段107が乗じる第一調整係数は1でよい。すなわち最も速く特性誤差を修正できる。
然し乍ら、実際にフィルタ誤差計算手段106が特性誤差を正しく計算し出力することは困難である。その原因の一つは、音響収音手段2が収音した収音信号に含まれる反響信号以外の周囲騒音、回路雑音などの外乱成分の影響である。この外乱成分が、フィルタ誤差計算手段106が出力する特性誤差の計算に悪影響を及ぼす。従って、この影響を避けるため、第一更新量調整手段107を導入し、第一調整係数を小さく設定したのが従来の技術である。
【0007】
従来は第一係数を小さく設定するため、フィルタ処理手段104で用いるフィルタ係数が反響消去のための真の特性に近づく速度は低くなり、完全な反響消去を達成するまでに長い時間を要していた。
本発明の目的は、上記のように、音響収音手段2が収音した収音信号に含まれる外乱成分の影響を最小限とし、第一更新調整手段107が極力大きな第一調整係数をとることができ、フィルタ処理手段104で用いるフィルタ係数が反響消去のための真の特性に近づく速度を高めることにある。
【0008】
【課題を解決するための手段】
この発明の請求項1では音響再生手段へと出力される再生信号を入力し一定時間蓄積し再生信号列を得る再生信号入力処理と、音響再生手段と同一空間に存在する音響収音手段から収音信号を入力し一定時間蓄積し収音信号列を得る収音信号入力処理と、再生信号列を周波数領域変換し再生信号変換列を得る再生信号変換処理と、再生信号変換列を周波数領域毎に入力しフィルタ処理することにより音響再生手段から音響収音手段へと回り込む反響信号の周波数領域変換列を模擬する模擬反響信号変換列を生成するフィルタ処理と、模擬反響信号変換列と収音信号列とを入力し模擬反響信号変換列の模擬誤差を出力する模擬誤差出力処理と、再生信号変換列と模擬誤差を入力し前記フィルタ処理のフィルタ係数の周波数領域毎の特性誤差を計算するフィルタ誤差計算処理と、フィルタ誤差計算処理が計算した特性誤差に周波数領域毎に第一調整係数を乗じる第一更新量調整処理と、第一調整係数を乗じられた特性誤差を加えることによりフィルタ処理のフィルタ係数を更新するフィルタ更新処理とを実行する反響消去方法において、模擬誤差出力処理で出力される模擬誤差の大きさに対する収音信号に含まれる反響信号以外の外乱信号成分の大きさの事前に測定した推定値の比率を外乱比率として計算し、外乱比率が0から所定値の範囲を逸脱した場合は範囲内になるように数値の打ち切りを行うことにより、外乱比率を周波数領域毎に逐次推定する外乱比率推定処理と、第一更新量調整処理の前段又は後段において、フィルタ係数の特性誤差又は第一調整係数を乗じたフィルタ係数の特性誤差に、0から所定値の間で、外乱比率の大きさに応じ、外乱比率が大きい程小さい値をとる第二調整係数を周波数領域毎に乗じる第二更新量調整処理とを実行する反響消去方法を提案する。
【0009】
この発明の請求項2では音響再生手段へと出力される再生信号を入力し一定時間蓄積し再生信号列を得る再生信号入力手段と、音響再生手段と同一空間に存在する音響収音手段から収音信号を入力し一定時間蓄積し収音信号列を得る収音信号入力手段と、再生信号列を周波数領域変換し再生信号変換列を得る再生信号変換手段と、再生信号変換列を周波数領域毎に入力しフィルタ処理することにより音響再生手段から音響収音手段へと回り込む反響信号の周波数領域変換列を模擬する模擬反響信号変換列を生成するフィルタ処理手段と、模擬反響信号変換列と収音信号列とを入力し模擬反響信号変換列の模擬誤差を出力する模擬誤差出力手段と、再生信号変換列と模擬誤差を入力しフィルタ処理手段のフィルタ係数の周波数領域毎の特性誤差を計算するフィルタ誤差計算手段と、フィルタ誤差計算手段が計算した特性誤差に周波数領域毎に第一調整係数を乗じる第一更新量調整手段と、第一調整係数を乗じられた特性誤差を加えることによりフィルタ処理手段のフィルタ係数を更新するフィルタ更新手段とを有する反響消去装置において、模擬誤差出力手段が出力する模擬誤差の大きさに対する収音信号に含まれる反響信号以外の外乱信号成分の大きさの事前に測定した推定値の比率を外乱比率として計算し、外乱比率が0から所定値の範囲を逸脱した場合は範囲内になるように数値の打ち切りを行うことにより、外乱比率を周波数領域毎に逐次推定する外乱比率推定手段と、第一更新量調整手段の前段又は後段において、フィルタ係数の特性誤差又は第一調整係数を乗じたフィルタ係数の特性誤差に、0から所定値の間で外乱比率の大きさに応じ、外乱比率が大きい程小さい値をとる第二調整係数を周波数領域毎に乗じる第二更新量調整手段とを有する反響消去装置を提案する。
【0010】
この発明の請求項ではコンピュータが解読可能な符号列によって記述され、コンピュータに請求項1又は3記載の反響消去方法を実行させる反響消去プログラムを提案する。
作用
この発明の構成によれば音響収音手段2が収音した収音信号に含まれる外乱成分が、フィルタ誤差計算手段が出力するフィルタ係数の特性誤差に及ぼした比率に応じて、計算されたフィルタ係数の特性誤差の不確かな部分を小さくすることができる。このため、第一更新量調整手段では極力大きな第一調整係数をとることができ、フィルタ処理手段で用いるフィルタ係数が反響消去のための真の近づく速度を高めることができる。
【0011】
外乱比率推定手段における、模擬誤差出力手段が出力する模擬誤差の大きさに対する音響収音手段からの収音信号に含まれる反響信号以外の外乱信号成分の大きさの比率を表わす外乱比率の推定は、次のように実現できる。ここで、N分割した周波数領域ω(n=1…N)における模擬誤差出力手段が出力する模擬誤差の大きさを|E(ωn)|とし、周波数領域ωnにおける外乱の大きさの推定値を|N(ωn)|とする。外乱推定値|N(ωn)|は、無音声区間等に事前測定した値などを用いるものとする。このとき外乱比率DR(ωn)の計算は、
DR(ωn)=|N(ωn)|/|E(ωn)| …(1)
として計算できる。ここで、外乱比率DR(ωn)の値は、0から所定値Tで例えばT=1の値をとるものとし、実使用環境で、この範囲を逸脱した場合は、0から所定値の範囲になるよう数値の打ち切りを行なう。DR(ωn)の値は、|E(ωn)|に対する|N(ωn)|の割合が大きくなるほど、大きな値を与えるようにすれば、(1)式以外に基づいてもよい。
【0012】
また、第二更新量調整手段の第二調整係数β(ωn)は0から所定値例えば1までの値をとり、上記の外乱比率DR(ωn)を用いて、
β(ωn)=(T−DR(ωnp1/p …(2)
のように決定できる。ここで、pは正の実数とする。また、(2)式の他、図4のAに示すように、DR(ωn)の値に対して滑らかに対応させたり、図4のBにしめすように、不連続に対応するようなものであっても、DR(ωn)の値が大きくなるほど、β(ωn)の値が小さくなるようなものであれば、本発明の効果を得られる実施例に含まれる。DR(ωn)が(1)式以外の形で求められた場合でも、ここに示したDR(ωn)とβ(ωn)の対応の中から適切なものを選択可能である。
【0013】
【発明の実施の形態】
図1にこの発明による反響消去装置の一実施例を示す。図1に示す反響消去装置を反響消去装置方法と共に説明する。
図1に示す100はこの発明による反響消去装置を示す。この発明による反響消去装置100は図6に示した従来の反響消去装置10に以下の手段が付加される。すなわち、模擬誤差出力手段105が出力する模擬誤差の大きさに対する音響収音手段2からの収音信号に含まれる反響信号以外の外乱信号成分の大きさの比率を表わす外乱比率DR(ωn)を周波数領域毎に逐次推定する外乱比率推定手段111と、実際にフィルタ誤差計算手段106が出力するフィルタ係数の特性誤差に、外乱比率推定手段111が推定した外乱比率の大きさに応じ、外乱比率が0パーセントの時は所定値、例えば1、外乱比率が100パーセントの時は0、その間の外乱比率の時は0から所定値の間の値をとる第二調整係数β(ωn)を周波数領域毎に乗じる第二更新量調整手段110とを含む。
【0014】
この構成により、音響収音手段2が収音した収音信号に含まれる外乱成分が、フィルタ誤差計算手段106が出力するフィルタ係数の特性誤差に及ぼした比率に応じて、計算されたフィルタ係数の特性誤差の不確かな部分の値を小さくすることができる。このため、第一更新量調整手段107では極力大きな第一調整係数をとることができ、フィルタ処理手段104で用いるフィルタ係数が反響消去のための真の特性に近づく速度を高めることができる。
【0015】
発明の変形実施例
図1に示した実施例では第二更新量調整手段110を第一更新量調整手段107の前段に設けた例を説明したが、図2に示すように第二更新量調整手段110を第一更新量調整手段107の後段に配置しても図1に示した実施例と同様の作用効果を得ることができる。
【0016】
また、図3に示すように、フィルタ誤差計算手段106とフィルタ更新手段108の間に第三更新量調整手段112を配置し、この第三更新量調整手段112で第一調整係数と第二調整係数との例えば積を求め、この積の値を第三調整係数とし、この第三調整係数を周波数領域毎のフィルタ係数の特性誤差に乗じる構成としても図1に示した実施例と同様の作用効果を得ることができる。
以上説明したこの発明による反響消去装置および反響消去方法はコンピュータが解読可能な符号列によって記述された音響消去プログラムをコンピュータに実行させて実現される。図5にこの発明による音響消去プログラムをコンピュータ20にインストールした状況を示す。
【0017】
コンピュータ20はよく知られているように、中央演算処理装置(CPU)21と、コンピュータの立上げおよび立下げ等を制御する基本プログラム等を記録した読み出専用メモリ(ROM)22と、データを一時記録する他に、この発明の反響消去方法を実行するためのプログラムを格納する読み出し、読み出し・書き込み可能なメモリ(RAM)23と、入力ポート24、出力ポート25等により構成することができる。
入力ポート24には再生信号入力端子3と、音響収音手段2が接続され、再生信号Xと、収音信号とが入力される。また、出力ポート25には音響再生手段1と反響消去信号出力端子4とが接続され、再生信号に対応する音響と、反響消去信号e=y−y´が出力される。
【0018】
読み出し、書き込み可能なメモリ23にはデータ格納領域23Aの他に、再生信号入力手段101として動作する再生信号入力処理プログラム23Bと、収音信号入力手段102として動作する収音信号入力処理プログラム23Cと、再生信号変換手段103として動作する再生信号変換処理プログラム23Dと、フィルタ処理手段104として動作するフィルタ処理プログラム23Eと、模擬誤差出力手段105として動作する模擬誤差出力処理プログラム23Fと、フィルタ誤差計算手段106として動作するフィルタ誤差計算処理プログラム23Gと、第一更新量調整手段107として動作する第一更新量調整処理プログラム23Hと、フィルタ更新手段108として動作するフィルタ更新処理プログラム23Iと、第二更新量調整手段110として動作する第二更新量調整処理プログラム23Jと、外乱比率推定手段111として動作する外乱比率推定処理プログラム23Kとがインストールされ、これらの各プログラムが中央演算処理装置21によって解読されて反響消去動作を実行する。読み出し、書き込み可能なメモリ23にインストールされた各プログラム23B〜23Kは予めコンピュータが読み出し可能な磁気ディスク或はCD−ROM等に記録され、これらの記録媒体からインストールされるか又は通信回路を通じてインストールされる。
【0019】
【発明の効果】
本発明による反響消去装置は、フィルタ処理手段で用いるフィルタ係数が反響消去のための真の特性に近づく速度の低下を抑えながら、音響収音手段が収音した収音信号に含まれる外乱成分の影響を低減でき、反響消去性能を高めることができる。
【図面の簡単な説明】
【図1】この発明による反響消去装置の一実施例を説明するためのブロック図。
【図2】この発明による反響消去装置の他の実施例を説明するためのブロック図。
【図3】この発明による反響消去装置の更に他の実施例を説明するためのブロック図。
【図4】この発明の要部の動作を説明するための特性曲線図。
【図5】この発明による反響消去装置をコンピュータで実現した場合の構成を説明するための構成概念図。
【図6】反響消去装置の概要を説明するためのブロック図。
【図7】従来の反響消去装置を説明するためのブロック図。
【符号の説明】
1 音響再生手段 105 模擬誤差出力手段
2 音響収音手段 106 フィルタ誤差計算手段
3 再生信号入力端子 107 第一更新量調整手段
4 反響消去信号出力端子 108 フィルタ更新手段
100 反響消去装置 110 第二更新量調整手段
101 再生信号入力手段 111 外乱比率推定手段
102 収音信号入力手段
103 再生信号変換手段
104 フィルタ処理手段
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an echo canceling apparatus for canceling echo that circulates from an acoustic reproduction means such as a speaker to an acoustic sound collection means such as a microphone.
[0002]
[Prior art]
The echo canceling device that cancels the echo that circulates from the sound reproducing means to the sound collecting means is connected as shown in FIG. Reference numeral 10 shown in FIG. 6 denotes a conventional echo canceling apparatus. In the conventional echo canceling apparatus 10, the impulse response h of the echo path between the sound reproducing means 1 and the sound collecting means 2 is estimated, and the estimated impulse response h ′ and the reproduced signal x input to the reproduced signal input terminal 3 are estimated. generating a convolution h'* x, it is subtracted from the actual echo signal y, to obtain a echo cancellation signal e to the echo cancellation signal output terminal 4.
However, the estimated impulse response and the reproduction signal convolution calculation require a large amount of calculation, which is a problem in implementation.
[0003]
In recent years, in order to solve this problem, a reproduction signal and a reverberation signal are once subjected to frequency domain transformation, a parameter corresponding to the frequency domain transformation of the impulse response of the echo path is estimated, and multiplication processing is used instead of convolution (non-convolution). Patent Document 1) or a method of reducing the amount of calculation has been proposed by dividing into smaller convolution operations (Patent Document 1). Examples of frequency domain transforms include (fast) discrete Fourier transform, (fast) discrete cosine transform, and (fast) Hartley transform.
In FIG. 6, a speaker is used as the sound reproducing unit 1, but the sound reproducing unit 1 includes an amplifier and a buffer in the previous stage of reproduction. Similarly, the sound pickup means 2 includes an amplifier and a buffer after the microphone.
[0004]
Moreover, the structure described in (nonpatent literature 1) etc. is typically shown in FIG. The configuration of FIG. 7 includes the following means. That is, the reproduction signal input means 101 that inputs the reproduction signal that is input from the reproduction signal input terminal 3 and is output to the sound reproduction means 1 and accumulates for a certain period of time to obtain a reproduction signal string, and the sound reproduction means 1 exist in the same space. A sound collection signal input means 102 for inputting a sound collection signal from the sound sound collection means 2 and storing it for a certain period of time to obtain a sound collection signal string; and a reproduction signal conversion means 103 for frequency domain conversion of the reproduction signal string to obtain a reproduction signal conversion string. Then, a reproduction signal conversion sequence is input for each frequency domain and filtered to generate a simulated reverberation signal conversion sequence that simulates the frequency domain conversion sequence of the reverberation signal that circulates from the sound reproduction means 1 to the sound collection means 2. Filter processing means 104, simulated error output means 105 for inputting a simulated echo signal conversion sequence and a collected sound signal sequence and outputting a simulated error of the simulated echo signal conversion sequence, reproduction signal conversion sequence and simulated error The filter error calculation means 106 that calculates the characteristic error of each filter domain of the filter coefficient of the filter processing means 104 is input, and the characteristic error calculated by the filter error calculation means 106 is multiplied by the first adjustment coefficient for each frequency domain. An update amount adjustment unit 107 and a filter update unit 108 that updates the filter coefficient of the filter processing unit 104 by adding a characteristic error multiplied by the first adjustment coefficient are included.
[0005]
[Patent Document 1]
JP-A-9-116472 (FIG. 6)
[Non-Patent Document 1]
Simon Haykin "Adaptive Filter Theory" Science and Technology Publishing, January 10, 2001, pp. 500-541
[0006]
[Problems to be solved by the invention]
In the configuration of FIG. 7, if the filter error calculation unit 106 can correctly calculate and output the characteristic error between the true characteristic for echo cancellation and the filter coefficient used in the filter processing unit 104, the filter error calculation unit 106 outputs the characteristic error. The first adjustment coefficient multiplied by the first update amount adjustment means 107 may be 1 for the characteristic error. That is, the characteristic error can be corrected most quickly.
However, it is difficult for the filter error calculation means 106 to actually calculate and output the characteristic error correctly. One of the causes is the influence of disturbance components such as ambient noise and circuit noise other than the echo signal included in the collected sound signal collected by the sound collecting means 2. This disturbance component adversely affects the calculation of the characteristic error output from the filter error calculation means 106. Therefore, in order to avoid this influence, the conventional technique is to introduce the first update amount adjusting means 107 and set the first adjustment coefficient small.
[0007]
Conventionally, since the first coefficient is set small, the speed at which the filter coefficient used in the filter processing unit 104 approaches the true characteristic for echo cancellation is low, and it takes a long time to achieve complete echo cancellation. It was.
As described above, the object of the present invention is to minimize the influence of disturbance components included in the collected sound signal picked up by the sound pickup means 2 and the first update adjustment means 107 takes the largest first adjustment coefficient. It is possible to increase the speed at which the filter coefficient used in the filter processing unit 104 approaches the true characteristic for echo cancellation.
[0008]
[Means for Solving the Problems]
According to the first aspect of the present invention, a reproduction signal input process for obtaining a reproduction signal sequence by inputting a reproduction signal to be output to the sound reproduction means and obtaining a reproduction signal string, and collecting from a sound collection means existing in the same space as the sound reproduction means. Sound signal input processing that inputs sound signals and accumulates them for a certain period of time to obtain a sound collection signal sequence, reproduction signal conversion processing that converts a reproduction signal sequence to a frequency domain and obtains a reproduction signal conversion sequence, and a reproduction signal conversion sequence for each frequency domain Filter processing for generating a simulated reverberation signal conversion sequence for simulating the frequency domain conversion sequence of the reverberation signal that circulates from the sound reproduction means to the sound collection means by being input to the filter, and the simulated reverberation signal conversion sequence and the sound collection signal And a simulated error output process for outputting a simulated error of the simulated echo signal conversion sequence, and a reproduction signal conversion sequence and the simulated error are input, and a characteristic error for each frequency domain of the filter coefficient of the filter processing is measured. A filter error calculation process, a first update amount adjustment process for multiplying the characteristic error calculated by the filter error calculation process by a first adjustment coefficient for each frequency domain, and a characteristic error multiplied by the first adjustment coefficient. In the echo cancellation method for executing the filter update process for updating the filter coefficient of the process, the magnitude of the disturbance signal component other than the echo signal included in the collected sound signal with respect to the magnitude of the simulated error output in the simulated error output process The ratio of the estimated value measured in advance is calculated as the disturbance ratio, and when the disturbance ratio deviates from the predetermined value range from 0, the disturbance ratio is set for each frequency domain by censoring the numerical value so that it falls within the range. Filter factor multiplied by the characteristic error of the filter coefficient or the first adjustment coefficient before or after the disturbance ratio estimation process for successive estimation and the first update amount adjustment process A second update amount adjustment process is performed by multiplying the characteristic error of 0 by a predetermined value between 0 and a predetermined value in accordance with the magnitude of the disturbance ratio, and by multiplying a second adjustment coefficient that takes a smaller value as the disturbance ratio increases for each frequency region. We propose an echo cancellation method.
[0009]
According to the second aspect of the present invention, the reproduction signal input means for inputting the reproduction signal output to the sound reproduction means and accumulating it for a predetermined time to obtain a reproduction signal string, and the sound collection means existing in the same space as the sound reproduction means are collected. A sound collection signal input means for inputting a sound signal and accumulating for a certain time to obtain a sound collection signal string, a reproduction signal conversion means for obtaining a reproduction signal conversion string by converting the reproduction signal string in a frequency domain, and a reproduction signal conversion string for each frequency domain Filter processing means for generating a simulated reverberation signal conversion sequence for simulating the frequency domain conversion sequence of the reverberation signal that circulates from the sound reproduction means to the sound collection means by being input to the filter, and a simulated reverberation signal conversion sequence and sound collection A simulated error output means for inputting a signal string and outputting a simulated error of the simulated echo signal conversion string, and a characteristic error for each frequency domain of the filter coefficient of the filter processing means by inputting the reproduction signal conversion string and the simulated error A filter error calculating means for calculating, a first update amount adjusting means for multiplying the characteristic error calculated by the filter error calculating means by a first adjustment coefficient for each frequency region, and a characteristic error multiplied by the first adjustment coefficient. In the echo canceller having the filter update means for updating the filter coefficient of the filter processing means, the magnitude of the disturbance signal component other than the echo signal included in the collected sound signal with respect to the magnitude of the simulated error output by the simulated error output means. The ratio of the estimated value measured in advance is calculated as the disturbance ratio, and when the disturbance ratio deviates from the predetermined value range from 0, the disturbance ratio is set for each frequency domain by censoring the numerical value so that it falls within the range. The filter multiplied by the characteristic error of the filter coefficient or the first adjustment coefficient in the previous stage or subsequent stage of the disturbance ratio estimation means for sequentially estimating and the first update amount adjustment means And a second update amount adjusting means for multiplying the characteristic error of the coefficient by a second adjustment coefficient that takes a smaller value as the disturbance ratio increases in accordance with the magnitude of the disturbance ratio between 0 and a predetermined value. An erasing device is proposed.
[0010]
According to a fifth aspect of the present invention, there is proposed an echo canceling program which is described by a code string readable by a computer and causes the computer to execute the echo canceling method according to the first or third aspect .
Disturbance component contained in the action collected sound signal sound pickup unit 2 is picked up according to the construction of this invention, in accordance with the ratio of the filter error calculation means had the characteristic error of the filter coefficients to be output, is calculated The uncertain part of the characteristic error of the filter coefficient can be reduced. For this reason, the first update amount adjusting means can take as large a first adjustment coefficient as possible, and the true speed at which the filter coefficient used in the filter processing means approaches for echo cancellation can be increased.
[0011]
In the disturbance ratio estimation means, the estimation of the disturbance ratio representing the ratio of the magnitude of the disturbance signal component other than the echo signal included in the sound pickup signal from the sound pickup means to the magnitude of the simulation error output by the simulation error output means is It can be realized as follows. Here, the magnitude of the simulated error output by the simulated error output means in the frequency domain ω (n = 1... N) divided into N is assumed to be | E (ω n ) |, and the magnitude of the disturbance in the frequency domain ω n is estimated. Let the value be | N (ω n ) |. As the disturbance estimated value | N (ω n ) |, a value measured in advance in a silent section or the like is used. At this time, the disturbance ratio DR (ω n ) is calculated as follows:
DR (ω n ) = | N (ω n ) | / | E (ω n ) | (1)
Can be calculated as Here, the value of the disturbance ratio DR (ω n ) is a value from 0 to a predetermined value T, for example, T = 1. If the value deviates from this range in an actual use environment, the range is from 0 to a predetermined value. Truncate the numerical value so that The value of DR (ω n ) may be based on other than equation (1) as long as the ratio of | N (ω n ) | to | E (ω n ) |
[0012]
Further, the second adjustment coefficient β (ω n ) of the second update amount adjusting means takes a value from 0 to a predetermined value, for example 1, and uses the above disturbance ratio DR (ω n ),
β (ω n ) = (T−DR (ω n ) p ) 1 / p (2)
Can be determined as follows. Here, p is a positive real number. Further, in addition to the expression (2), as shown in A of FIG. 4, the value of DR (ω n ) is made to correspond smoothly, or as shown in FIG. Even if the value of DR (ω n ) increases, the value of β (ω n ) decreases as the value of DR (ω n ) increases. Even when DR (ω n ) is obtained in a form other than equation (1), an appropriate one can be selected from the correspondence between DR (ω n ) and β (ω n ) shown here.
[0013]
DETAILED DESCRIPTION OF THE INVENTION
FIG. 1 shows an embodiment of an echo canceling apparatus according to the present invention. The echo canceling apparatus shown in FIG. 1 will be described together with the echo canceling apparatus method.
Reference numeral 100 shown in FIG. 1 represents an echo canceling apparatus according to the present invention. In the echo canceling apparatus 100 according to the present invention, the following means are added to the conventional echo canceling apparatus 10 shown in FIG. That is, a disturbance ratio DR (ω n ) representing the ratio of the magnitude of the disturbance signal component other than the echo signal included in the sound pickup signal from the sound pickup means 2 with respect to the magnitude of the simulation error output from the simulation error output means 105. The disturbance ratio estimating means 111 for successively estimating the disturbance ratio according to the magnitude of the disturbance ratio estimated by the disturbance ratio estimating means 111 to the characteristic error of the filter coefficient actually output by the filter error calculating means 106 The frequency of the second adjustment coefficient β (ω n ), which is a predetermined value when 1 is 0%, for example, 1 when the disturbance ratio is 100%, and 0 when the disturbance ratio is between 0 and the predetermined value. Second update amount adjustment means 110 for multiplying each area.
[0014]
With this configuration, the disturbance coefficient included in the collected sound signal picked up by the sound pickup means 2 has the calculated filter coefficient according to the ratio of the characteristic error of the filter coefficient output from the filter error calculation means 106. The value of the uncertain part of the characteristic error can be reduced. For this reason, the first update amount adjusting means 107 can take as large a first adjustment coefficient as possible, and the speed at which the filter coefficient used in the filter processing means 104 approaches the true characteristic for echo cancellation can be increased.
[0015]
Modified embodiment of the invention In the embodiment shown in Fig. 1, the second update amount adjusting means 110 is provided in the preceding stage of the first update amount adjusting means 107, but as shown in Fig. 2, Even if the two update amount adjusting means 110 is arranged at the subsequent stage of the first update amount adjusting means 107, the same effect as that of the embodiment shown in FIG. 1 can be obtained.
[0016]
Further, as shown in FIG. 3, a third update amount adjustment unit 112 is arranged between the filter error calculation unit 106 and the filter update unit 108, and the first adjustment coefficient and the second adjustment are adjusted by the third update amount adjustment unit 112. For example, a product with a coefficient is obtained, the value of this product is set as a third adjustment coefficient, and the characteristic of the filter coefficient for each frequency domain is multiplied by the same function as that of the embodiment shown in FIG. An effect can be obtained.
The echo canceling apparatus and the echo canceling method according to the present invention described above are realized by causing a computer to execute a sound canceling program described by a computer-readable code string. FIG. 5 shows a state in which the sound erasing program according to the present invention is installed in the computer 20.
[0017]
As is well known, the computer 20 has a central processing unit (CPU) 21, a read-only memory (ROM) 22 in which basic programs for controlling the startup and shutdown of the computer, etc. are recorded, In addition to the temporary recording, the read / write memory (RAM) 23 storing the program for executing the echo canceling method of the present invention, the input port 24, the output port 25, and the like can be used.
The reproduction signal input terminal 3 and the sound pickup means 2 are connected to the input port 24, and the reproduction signal X and the sound collection signal are input thereto. Further, the sound reproducing means 1 and the echo canceling signal output terminal 4 are connected to the output port 25, and the sound corresponding to the reproducing signal and the echo canceling signal e = yy ′ are output.
[0018]
In addition to the data storage area 23A, the readable / writable memory 23 includes a reproduction signal input processing program 23B that operates as the reproduction signal input means 101, and a sound collection signal input processing program 23C that operates as the sound collection signal input means 102. A reproduction signal conversion processing program 23D operating as the reproduction signal conversion means 103, a filter processing program 23E operating as the filter processing means 104, a simulation error output processing program 23F operating as the simulation error output means 105, and a filter error calculation means A filter error calculation processing program 23G that operates as 106, a first update amount adjustment processing program 23H that operates as first update amount adjustment means 107, a filter update processing program 23I that operates as filter update means 108, and a second update amount Adjusting hand The second update amount adjustment processing program 23J that operates as 110 and the disturbance ratio estimation processing program 23K that operates as the disturbance ratio estimation means 111 are installed, and each of these programs is decoded by the central processing unit 21 to perform echo cancellation operation Execute. The programs 23B to 23K installed in the readable and writable memory 23 are recorded in advance on a magnetic disk or CD-ROM that can be read by a computer, and installed from these recording media or installed through a communication circuit. The
[0019]
【The invention's effect】
The echo canceling apparatus according to the present invention suppresses a decrease in the speed at which the filter coefficient used in the filter processing means approaches the true characteristic for echo canceling, while reducing the disturbance component included in the collected sound signal collected by the acoustic sound collecting means. The influence can be reduced and the echo canceling performance can be enhanced.
[Brief description of the drawings]
FIG. 1 is a block diagram for explaining an embodiment of an echo canceling apparatus according to the present invention;
FIG. 2 is a block diagram for explaining another embodiment of the echo canceling apparatus according to the present invention;
FIG. 3 is a block diagram for explaining still another embodiment of the echo canceling apparatus according to the present invention.
FIG. 4 is a characteristic curve diagram for explaining the operation of the main part of the present invention.
FIG. 5 is a conceptual diagram for explaining the configuration when the echo canceling apparatus according to the present invention is realized by a computer.
FIG. 6 is a block diagram for explaining the outline of the echo canceling apparatus.
FIG. 7 is a block diagram for explaining a conventional echo canceling apparatus.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 Sound reproduction | regeneration means 105 Simulated error output means 2 Sound pickup means 106 Filter error calculation means 3 Playback signal input terminal 107 1st update amount adjustment means 4 Echo cancellation signal output terminal 108 Filter update means 100 Echo cancellation apparatus 110 2nd update amount Adjustment means 101 Reproduction signal input means 111 Disturbance ratio estimation means 102 Sound pickup signal input means 103 Reproduction signal conversion means 104 Filter processing means

Claims (5)

音響再生手段へと出力される再生信号を入力し一定時間蓄積し再生信号列を得る再生信号入力処理と、前記音響再生手段と同一空間に存在する音響収音手段から収音信号を入力し一定時間蓄積し収音信号列を得る収音信号入力処理と、前記再生信号列を周波数領域変換し再生信号変換列を得る再生信号変換処理と、前記再生信号変換列を周波数領域毎に入力しフィルタ処理することにより前記音響再生手段から前記音響収音手段へと回り込む反響信号の周波数領域変換列を模擬する模擬反響信号変換列を生成するフィルタ処理と、前記模擬反響信号変換列と前記収音信号列とを入力し模擬反響信号変換列の模擬誤差を出力する模擬誤差出力処理と、前記再生信号変換列と前記模擬誤差を入力し前記フィルタ処理のフィルタ係数の周波数領域毎の特性誤差を計算するフィルタ誤差計算処理と、前記フィルタ誤差計算処理が計算した特性誤差に周波数領域毎に第一調整係数を乗じる第一更新量調整処理と、前記第一調整係数を乗じられた特性誤差を加えることにより前記フィルタ処理のフィルタ係数を更新するフィルタ更新処理とを実行する反響消去方法において、
前記模擬誤差出力処理で出力される前記模擬誤差の大きさに対する前記収音信号に含まれる前記反響信号以外の外乱信号成分の大きさの事前に測定した推定値の比率を外乱比率として計算し、前記外乱比率が0から所定値の範囲を逸脱した場合は前記範囲内になるように数値の打ち切りを行うことにより、外乱比率を周波数領域毎に逐次推定する外乱比率推定処理と、前記第一更新量調整処理の前段又は後段において、前記フィルタ係数の特性誤差又は前記第一調整係数を乗じた前記フィルタ係数の特性誤差に、0から所定値の間で、前記外乱比率の大きさに応じ、外乱比率が大きい程小さい値をとる第二調整係数を周波数領域毎に乗じる第二更新量調整処理とを実行することを特徴とする反響消去方法。
A reproduction signal input process for inputting a reproduction signal output to the sound reproduction means and accumulating a reproduction signal sequence for a certain period of time, and a sound collection signal inputted from the sound collection means existing in the same space as the sound reproduction means A sound collection signal input process for accumulating time to obtain a sound collection signal string, a reproduction signal conversion process for obtaining a reproduction signal conversion string by converting the reproduction signal string in a frequency domain, and a filter for inputting the reproduction signal conversion string for each frequency domain Filter processing for generating a simulated echo signal conversion sequence that simulates a frequency domain conversion sequence of an echo signal that circulates from the acoustic reproduction means to the acoustic sound collection means by processing, the simulated echo signal conversion sequence, and the collected sound signal A simulated error output process for inputting a sequence and outputting a simulated error of the simulated echo signal conversion sequence; and for each frequency domain of the filter coefficient of the filter processing by inputting the reproduction signal conversion sequence and the simulated error A filter error calculation process for calculating a characteristic error; a first update amount adjustment process for multiplying the characteristic error calculated by the filter error calculation process by a first adjustment coefficient for each frequency domain; and a characteristic multiplied by the first adjustment coefficient. In an echo canceling method for executing a filter update process for updating a filter coefficient of the filter process by adding an error,
Calculating a ratio of a pre-measured estimated value of the magnitude of a disturbance signal component other than the echo signal included in the collected sound signal to the magnitude of the simulated error output in the simulated error output process, as a disturbance ratio, When the disturbance ratio deviates from a predetermined value range from 0, a disturbance ratio estimation process for sequentially estimating the disturbance ratio for each frequency domain by performing numerical truncation so that the disturbance ratio falls within the range, and the first update Before or after the amount adjustment processing, the filter coefficient characteristic error or the filter coefficient characteristic error multiplied by the first adjustment coefficient is a disturbance between 0 and a predetermined value depending on the magnitude of the disturbance ratio. And a second update amount adjustment process in which a second adjustment coefficient that takes a smaller value as the ratio increases is multiplied for each frequency region.
音響再生手段へと出力される再生信号を入力し一定時間蓄積し再生信号列を得る再生信号入力手段と、前記音響再生手段と同一空間に存在する音響収音手段から収音信号を入力し一定時間蓄積し収音信号列を得る収音信号入力手段と、前記再生信号列を周波数領域変換し再生信号変換列を得る再生信号変換手段と、前記再生信号変換列を周波数領域毎に入力しフィルタ処理することにより前記音響再生手段から前記音響収音手段へと回り込む反響信号の周波数領域変換列を模擬する模擬反響信号変換列を生成するフィルタ処理手段と、前記模擬反響信号変換列と前記収音信号列とを入力し模擬反響信号変換列の模擬誤差を出力する模擬誤差出力手段と、前記再生信号変換列と前記模擬誤差を入力し前記フィルタ処理手段のフィルタ係数の周波数領域毎の特性誤差を計算するフィルタ誤差計算手段と、前記フィルタ誤差計算手段が計算した特性誤差に周波数領域毎に第一調整係数を乗じる第一更新調整手段と、前記第一調整係数を乗じられた特性誤差を加えることにより前記フィルタ処理手段のフィルタ係数を更新するフィルタ更新手段とを有する反響消去装置において、
前記模擬誤差出力手段が出力する前記模擬誤差の大きさに対する前記収音信号に含まれる前記反響信号以外の外乱信号成分の大きさの事前に測定した推定値の比率を外乱比率として計算し、前記外乱比率が0から所定値の範囲を逸脱した場合は前記範囲内になるように数値の打ち切りを行うことにより、外乱比率を周波数領域毎に逐次推定する外乱比率推定手段と、前記第一更新量調整手段の前段又は後段において、前記フィルタ係数の特性誤差又は前記第一調整係数を乗じた前記フィルタ係数の特性誤差に、0から所定値の間で前記外乱比率の大きさに応じ、外乱比率が大きい程小さい値をとる第二調整係数を周波数領域毎に乗じる第二更新量調整手段とを有することを特徴とする反響消去装置。
A reproduction signal input means for inputting a reproduction signal output to the sound reproduction means and accumulating for a certain period of time to obtain a reproduction signal sequence; and a sound collection signal inputted from the sound collection means existing in the same space as the sound reproduction means and constant Sound collection signal input means for accumulating time and obtaining a sound collection signal string, reproduction signal conversion means for obtaining a reproduction signal conversion string by frequency domain conversion of the reproduction signal string, and inputting and reproducing the reproduction signal conversion string for each frequency domain Filter processing means for generating a simulated reverberation signal conversion sequence for simulating a frequency domain conversion sequence of an echo signal that circulates from the sound reproduction means to the sound collection means by processing, the simulated reverberation signal conversion sequence, and the sound collection A simulated error output means for inputting a signal string and outputting a simulated error of the simulated echo signal conversion string; a frequency of a filter coefficient of the filter processing means for inputting the reproduction signal conversion string and the simulated error; A filter error calculating means for calculating a characteristic error for each region; a first update adjusting means for multiplying the characteristic error calculated by the filter error calculating means for each frequency region; and the first adjustment coefficient. In the echo canceller having the filter update means for updating the filter coefficient of the filter processing means by adding the characteristic error,
Calculating a ratio of a pre-measured estimated value of the magnitude of a disturbance signal component other than the echo signal included in the collected sound signal to the magnitude of the simulated error output by the simulated error output means, as the disturbance ratio, When the disturbance ratio deviates from a predetermined value range from 0, a disturbance ratio estimation unit that sequentially estimates the disturbance ratio for each frequency region by performing numerical truncation so that the disturbance ratio falls within the range, and the first update amount Before or after the adjustment means, a disturbance ratio is determined according to the magnitude of the disturbance ratio between 0 and a predetermined value to the filter coefficient characteristic error or the filter coefficient characteristic error multiplied by the first adjustment coefficient. A reverberation canceling device comprising: a second update amount adjusting unit that multiplies a second adjustment coefficient that takes a smaller value as the value increases for each frequency region.
請求項1記載の反響消去方法であって、The echo cancellation method according to claim 1,
各前記周波数領域をωEach said frequency domain is ω n (n=1,2、・・・、N)とし、前記第二調整係数をβ(ω(N = 1, 2,..., N) and the second adjustment coefficient is β (ω n )とし、前記所定値をTとし、前記外乱比率をDR(ω), The predetermined value is T, and the disturbance ratio is DR (ω n )としたとき、前記第二調整係数β(ω), The second adjustment coefficient β (ω n )を)
β(ω      β (ω n )=(T−DR(ω) = (T-DR (ω n ) p ) 1/p1 / p (ただし、pは正の実数)      (Where p is a positive real number)
により求めることを特徴とする反響消去方法。The echo canceling method characterized by calculating | requiring by.
請求項2記載の反響消去装置であって、The echo canceling device according to claim 2,
各前記周波数領域をωEach said frequency domain is ω n (n=1,2、・・・、N)とし、前記第二調整係数をβ(ω(N = 1, 2,..., N) and the second adjustment coefficient is β (ω n )とし、前記所定値をTとし、前記外乱比率をDR(ω), The predetermined value is T, and the disturbance ratio is DR (ω n )としたとき、前記第二調整係数β(ω), The second adjustment coefficient β (ω n )を)
β(ω      β (ω n )=(T−DR(ω) = (T-DR (ω n ) p ) 1/p1 / p (ただし、pは正の実数)      (Where p is a positive real number)
により求めることを特徴とする反響消去装置。An echo canceling device characterized by the above.
コンピュータが解読可能な符号列によって記述され、コンピュータに請求項1又は3記載の反響消去方法を実行させる反響消去プログラム。An echo canceling program, which is described by a computer-readable code string, and causes the computer to execute the echo canceling method according to claim 1 or 3.
JP2003108875A 2003-04-14 2003-04-14 Echo canceling method, echo canceling device, echo canceling program Expired - Fee Related JP4156428B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003108875A JP4156428B2 (en) 2003-04-14 2003-04-14 Echo canceling method, echo canceling device, echo canceling program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003108875A JP4156428B2 (en) 2003-04-14 2003-04-14 Echo canceling method, echo canceling device, echo canceling program

Publications (2)

Publication Number Publication Date
JP2004320204A JP2004320204A (en) 2004-11-11
JP4156428B2 true JP4156428B2 (en) 2008-09-24

Family

ID=33470212

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003108875A Expired - Fee Related JP4156428B2 (en) 2003-04-14 2003-04-14 Echo canceling method, echo canceling device, echo canceling program

Country Status (1)

Country Link
JP (1) JP4156428B2 (en)

Also Published As

Publication number Publication date
JP2004320204A (en) 2004-11-11

Similar Documents

Publication Publication Date Title
JP5452655B2 (en) Multi-sensor voice quality improvement using voice state model
JP5634959B2 (en) Noise / dereverberation apparatus, method and program thereof
KR20070055599A (en) Method of processing a noisy sound signal and device for implementing said method
JP2008197284A (en) Filter coefficient calculation device, filter coefficient calculation method, control program, computer-readable recording medium, and audio signal processing apparatus
CN113436643A (en) Method, device, equipment and storage medium for training and applying speech enhancement model
US11238882B2 (en) Dry sound and ambient sound separation
CN112201273B (en) Noise power spectral density calculation method, system, equipment and medium
JP5151483B2 (en) Coefficient measuring device, effect applying device, and musical sound generating device
JP4787851B2 (en) Echo suppression gain estimation method, echo canceller using the same, device program, and recording medium
JP5161157B2 (en) Frequency domain echo removal apparatus, frequency domain echo removal method, program
JP3920795B2 (en) Echo canceling apparatus, method, and echo canceling program
JP4156428B2 (en) Echo canceling method, echo canceling device, echo canceling program
CN111035323A (en) Method, equipment and system for actively reducing noise of sweeping robot
JP2010044150A (en) Reverberation removing device and reverberation removing method, and program and recording medium thereof
JP5562451B1 (en) Echo suppression gain estimation method, echo canceller and program using the same
JP5583181B2 (en) Cascade connection type transmission system parameter estimation method, cascade connection type transmission system parameter estimation device, program
JP5325134B2 (en) Echo canceling method, echo canceling apparatus, program thereof, and recording medium
JP3917116B2 (en) Echo canceling apparatus, method, echo canceling program, and recording medium recording the program
JP4242320B2 (en) Voice recognition method, apparatus and program thereof, and recording medium thereof
JP4094522B2 (en) Echo canceling apparatus, method, echo canceling program, and recording medium recording the program
JP4094523B2 (en) Echo canceling apparatus, method, echo canceling program, and recording medium recording the program
JP6226465B2 (en) Audio signal processing apparatus, recording / reproducing apparatus, and program
JP4729072B2 (en) Echo erasing device, echo erasing method, echo erasing program, and recording medium
JP5086969B2 (en) Echo canceling apparatus, method thereof, program thereof, and recording medium thereof
CN116129928A (en) Near-end speech intelligibility enhancement method and system for broadcast communication scene

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050721

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20050721

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20070420

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080325

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080521

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: 20080701

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080709

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110718

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120718

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130718

Year of fee payment: 5

LAPS Cancellation because of no payment of annual fees