JPH06318846A - 適応フィルタ装置における係数更新方法 - Google Patents
適応フィルタ装置における係数更新方法Info
- Publication number
- JPH06318846A JPH06318846A JP34817793A JP34817793A JPH06318846A JP H06318846 A JPH06318846 A JP H06318846A JP 34817793 A JP34817793 A JP 34817793A JP 34817793 A JP34817793 A JP 34817793A JP H06318846 A JPH06318846 A JP H06318846A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- input signal
- input
- coefficient
- reference signal
- 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.)
- Granted
Links
Landscapes
- Complex Calculations (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
(57)【要約】
【目的】 係数更新の演算量が少ないのにも拘らず収束
性が良く、且つ適応誤差が小さい適応フィルタを提供す
る。 【構成】 適応フィルタの係数更新の演算のために、正
規化ＬＭＳ（学習同定法）の変形である符号法を更に変
形した方法を採用する。新しい方法では、サンプリング
時点の異なる複数の入力信号の絶対値の平均値を求め
る。入力信号が平均値よりも小さい時には符号法に従う
符号関数をゼロとする。入力信号が平均値以上の時には
従来の符号法と同様に処理する。
性が良く、且つ適応誤差が小さい適応フィルタを提供す
る。 【構成】 適応フィルタの係数更新の演算のために、正
規化ＬＭＳ（学習同定法）の変形である符号法を更に変
形した方法を採用する。新しい方法では、サンプリング
時点の異なる複数の入力信号の絶対値の平均値を求め
る。入力信号が平均値よりも小さい時には符号法に従う
符号関数をゼロとする。入力信号が平均値以上の時には
従来の符号法と同様に処理する。
Description
【０００１】
【産業上の利用分野】本発明は、ノイズキャンセリング
やシステムの同定等に使用されるディジタルフィルタの
１種である適応フィルタ装置における係数更新方法に関
する。
やシステムの同定等に使用されるディジタルフィルタの
１種である適応フィルタ装置における係数更新方法に関
する。
【０００２】
【従来の技術及び発明が解決しようとする課題】トラン
スバーサル型適応フィルタはサンプリング時点の異なる
複数の入力信号に対して係数をそれぞれ乗算する手段を
有する。各係数は所定のアルゴリズムに従って更新され
る。係数更新のための代表的なアルゴリズムとしてＬＭ
Ｓ（Least Mean Square）アルゴリズムがある。この
典型的なＬＭＳアルゴリズムの変形として正規化ＬＭＳ
（ＮＬＭＳ）又は学習同定法と呼ばれているアルコリズ
ムが例えば昭和４３年９月の「計測と制御」第７巻第９
号又は「ＩＥＥＥ トランザクション・オートマチック
・コントロール」Vol.ＡＣ−１２ １９６７年６月の野
田、南雲の論文で知られている。この正規化ＬＭＳアル
ゴリズムは収束速度の向上をはかることができるという
特長を有する反面、係数を求めるための計算量が多いと
いう欠点を有する。また、正規化ＬＭＳの変形として、
上記論文に“符号法”による正規化ＬＭＳが記載されて
いる。この符号法においては、サンプリング時点の異な
る複数の入力信号が０より大きい値の時には＋１の符
号、０より小さい値の時には−１の符号、０の時には０
の符号に置き換える処理を含むので、係数更新の計算量
が大幅に少なくなる。しかし、収束速度が遅くなり且つ
収束誤差も大きくなる。
スバーサル型適応フィルタはサンプリング時点の異なる
複数の入力信号に対して係数をそれぞれ乗算する手段を
有する。各係数は所定のアルゴリズムに従って更新され
る。係数更新のための代表的なアルゴリズムとしてＬＭ
Ｓ（Least Mean Square）アルゴリズムがある。この
典型的なＬＭＳアルゴリズムの変形として正規化ＬＭＳ
（ＮＬＭＳ）又は学習同定法と呼ばれているアルコリズ
ムが例えば昭和４３年９月の「計測と制御」第７巻第９
号又は「ＩＥＥＥ トランザクション・オートマチック
・コントロール」Vol.ＡＣ−１２ １９６７年６月の野
田、南雲の論文で知られている。この正規化ＬＭＳアル
ゴリズムは収束速度の向上をはかることができるという
特長を有する反面、係数を求めるための計算量が多いと
いう欠点を有する。また、正規化ＬＭＳの変形として、
上記論文に“符号法”による正規化ＬＭＳが記載されて
いる。この符号法においては、サンプリング時点の異な
る複数の入力信号が０より大きい値の時には＋１の符
号、０より小さい値の時には−１の符号、０の時には０
の符号に置き換える処理を含むので、係数更新の計算量
が大幅に少なくなる。しかし、収束速度が遅くなり且つ
収束誤差も大きくなる。
【０００３】そこで、本発明は、係数更新のための計算
量を符号法と同様に少なくし、収束速度を正規化ＬＭＳ
と同様に速くし、収束誤差を正規化ＬＭＳと符号法の中
間程度にすることができる適応フィルタ装置における係
数更新方法を提供することにある。
量を符号法と同様に少なくし、収束速度を正規化ＬＭＳ
と同様に速くし、収束誤差を正規化ＬＭＳと符号法の中
間程度にすることができる適応フィルタ装置における係
数更新方法を提供することにある。
【０００４】
【課題を解決するための手段】上記目的を達成するため
の本発明は、所定のサンプリング周期で基準信号ｄ(t)
を供給するための基準信号供給手段と、前記基準信号に
関係を有する入力信号ｘ(t) を所定のサンプリング周期
で供給するための入力信号供給手段と、前記基準信号供
給手段と前記入力信号供給手段に接続され、１サンプリ
ング周期を単位遅延時間として前記入力信号を順次に遅
延することによってサンプリング時点の異なる０番目か
らＬ−１番目までのＬ個（但し、Ｌは２よりも大きい正
の整数）の入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ
−Ｌ＋１）を同時に得、前記０番目からＬ番目の入力信
号に０番目からＬ−１番目の係数ｗ0 (t) 、ｗ1(t) 、
…ｗL-1 (t) を夫々乗算し、これによって得られたＬ個
の乗算出力を加算した出力ｙ(t) を得、前記出力ｙ(t)
と前記基準信号ｄ(t) との誤差信号ｅ(t)を求め、前記
入力信号ｘ(t) と前記誤差信号ｅ(t) とに基づいて前記
０番目からＬ−１番目の係数の更新値ｗ0 （ｔ＋１）、
ｗ1 （ｔ＋１）…ｗL-1 （ｔ＋１）を決定するためのデ
ィジタルフィルタ信号処理回路とを有する適応フィルタ
装置における係数更新方法であって、前記更新係数ｗ0
（ｔ＋１）、ｗ1 （ｔ＋１) 、…ｗL-1 （ｔ＋１）を、 ｗi （ｔ＋１）＝ｗi (t) ＋β［ｅ(t) ／｛Σj=0 L-1
x ′（ｔ−ｊ）ｘ（ｔ−ｊ）］ｘ′（ｔ−ｉ） ［ここで、各記号は次のものを示す。ｉ及びｊは０から
Ｌ−１までの正の整数；ｔはサンプリング時刻を示す記
号；ｗi （ｔ＋１）はＬ個の更新係数ｗ0 （ｔ＋１）、
ｗ1 （ｔ＋１）、…ｗL-1（ｔ＋１）；ｗi (t) はＬ個
の更新前の係数ｗ0 (t) 、ｗ1 (t) 、…ｗL-1 (t) ；β
は任意の定数；ｅ(t) はサンプリング時刻ｔにおける誤
差信号；ｘ（ｔ−ｊ）はサンプリング時点の異なるＬ個
の入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋
１）；ｘ′（ｔ−ｉ）及びｘ′（ｔ−ｊ）はサンプリン
グ時点の異なるＬ個の入力信号ｘ（ｔ）、ｘ（ｔ−
１）、…ｘ（ｔ−Ｌ＋１）に対応した変形入力信号ｘ′
（ｔ）、ｘ′（ｔ−１）、…ｘ′（ｔ−Ｌ＋１）、｛但
し、このｘ′（ｔ−ｉ）及びｘ′（ｔ−ｊ）は、 ｜ｘ（ｔ−ｉ）｜≧（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （ｃは１又はこれ以外の定数）が成立する時にｓｇｎ
［ｘ（ｔ−ｉ）］（ここで、ｓｇｎ［ｘ（ｔ−ｉ）］は
ｘ（ｔ−ｉ）＞０の時に＋１、ｘ（ｔ−ｉ）＜０の時は
−１、ｘ（ｔ−ｉ）＝０の時は０を示す。）であり、 ｜ｘ（ｔ−ｉ）｜＜（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にゼロである。］に従って決定することを
特徴とする適応フィルタ装置における係数更新方法に係
わるものである。また、請求項２の発明は、所定のサン
プリング周期で基準信号ｄ(t) を供給するための基準信
号供給手段と、前記基準信号に関係を有する入力信号ｘ
(t) を所定のサンプリング周期で供給するための入力信
号供給手段と、前記基準信号供給手段と前記入力信号供
給手段に接続され、１サンプリング周期を単位遅延時間
として前記入力信号を順次に遅延することによってサン
プリング時点の異なる０番目からＬ−１番目までのＬ個
（但し、Ｌは２よりも大きい正の整数）の入力信号ｘ
（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋１）を同時に
得、前記０番目からＬ番目の入力信号に０番目からＬ−
１番目の係数ｗ0 (t) 、ｗ1(t) 、…ｗL-1 (t) を夫々
乗算し、これによって得られたＬ個の乗算出力を加算し
た出力ｙ(t) を得、前記出力ｙ(t) と前記基準信号ｄ
(t) との誤差信号ｅ(t)を求め、前記入力信号ｘ(t) と
前記誤差信号ｅ(t) とに基づいて前記０番目からＬ−１
番目の係数の更新値ｗ0 （ｔ＋１）、ｗ1 （ｔ＋１）…
ｗL-1 （ｔ＋１）を決定するためのディジタルフィルタ
信号処理回路とを有する適応フィルタ装置における係数
更新方法であって、前記更新係数ｗ0 （ｔ＋１）、ｗ1
（ｔ＋１) 、…ｗL-1 （ｔ＋１）を、 ｗi （ｔ＋１）＝ｗi （ｔ）＋βｅ（ｔ）x ′（ｔ−
ｉ）ｘ ［ここで、各記号は次のものを示す。ｉは０からＬ−１
までの正の整数；ｔはサンプリング時刻を示す記号；ｗ
i （ｔ＋１）はＬ個の更新係数ｗ0 （ｔ＋１）、ｗ1
（ｔ＋１）、…ｗL-1（ｔ＋１）；ｗi (t) はＬ個の更
新前の係数ｗ0 (t) 、ｗ1 (t) 、…ｗL-1 (t) ；βは任
意の定数；ｅ(t) はサンプリング時刻ｔにおける誤差信
号；ｘ′（ｔ−ｉ）はサンプリング時点の異なるＬ個の
入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋１）
に対応した変形入力信号ｘ′（ｔ）、ｘ′（ｔ−１）、
…ｘ′（ｔ−Ｌ＋１）、｛但し、このｘ′（ｔ−ｉ）
は、 ｜ｘ′（ｔ−ｉ）｜≧（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （ｃは定数、ｊは０からＬ−１までの正の整数）が成立
する時にｓｇｎ［ｘ（ｔ−ｉ）］（ここで、ｓｇｎ［ｘ
（ｔ−ｉ）］はｘ（ｔ−ｉ）＞０の時に＋１、ｘ（ｔ−
ｉ）＜０の時は−１、ｘ（ｔ−ｉ）＝０の時は０を示
す。）であり、 ｜ｘ′（ｔ−ｉ）｜＜（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にゼロである。］に従って決定することを
特徴とする適応フィルタ装置における係数更新方法に係
わるものである。なお、本願において、Σj=0 L-1 は
ｘ′（ｔ−ｊ）ｘ（ｔ−ｊ）のｊを０からＬ−１まで変
化させたＬ個の値の総和を求めることを示し、またΣ
j=0 L-1 ｜ｘ（ｔ−ｊ）｜は｜ｘ（ｔ−ｊ）｜のｊを０
からＬ−１まで変化させたＬ個の値の総和を求めること
を示す。
の本発明は、所定のサンプリング周期で基準信号ｄ(t)
を供給するための基準信号供給手段と、前記基準信号に
関係を有する入力信号ｘ(t) を所定のサンプリング周期
で供給するための入力信号供給手段と、前記基準信号供
給手段と前記入力信号供給手段に接続され、１サンプリ
ング周期を単位遅延時間として前記入力信号を順次に遅
延することによってサンプリング時点の異なる０番目か
らＬ−１番目までのＬ個（但し、Ｌは２よりも大きい正
の整数）の入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ
−Ｌ＋１）を同時に得、前記０番目からＬ番目の入力信
号に０番目からＬ−１番目の係数ｗ0 (t) 、ｗ1(t) 、
…ｗL-1 (t) を夫々乗算し、これによって得られたＬ個
の乗算出力を加算した出力ｙ(t) を得、前記出力ｙ(t)
と前記基準信号ｄ(t) との誤差信号ｅ(t)を求め、前記
入力信号ｘ(t) と前記誤差信号ｅ(t) とに基づいて前記
０番目からＬ−１番目の係数の更新値ｗ0 （ｔ＋１）、
ｗ1 （ｔ＋１）…ｗL-1 （ｔ＋１）を決定するためのデ
ィジタルフィルタ信号処理回路とを有する適応フィルタ
装置における係数更新方法であって、前記更新係数ｗ0
（ｔ＋１）、ｗ1 （ｔ＋１) 、…ｗL-1 （ｔ＋１）を、 ｗi （ｔ＋１）＝ｗi (t) ＋β［ｅ(t) ／｛Σj=0 L-1
x ′（ｔ−ｊ）ｘ（ｔ−ｊ）］ｘ′（ｔ−ｉ） ［ここで、各記号は次のものを示す。ｉ及びｊは０から
Ｌ−１までの正の整数；ｔはサンプリング時刻を示す記
号；ｗi （ｔ＋１）はＬ個の更新係数ｗ0 （ｔ＋１）、
ｗ1 （ｔ＋１）、…ｗL-1（ｔ＋１）；ｗi (t) はＬ個
の更新前の係数ｗ0 (t) 、ｗ1 (t) 、…ｗL-1 (t) ；β
は任意の定数；ｅ(t) はサンプリング時刻ｔにおける誤
差信号；ｘ（ｔ−ｊ）はサンプリング時点の異なるＬ個
の入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋
１）；ｘ′（ｔ−ｉ）及びｘ′（ｔ−ｊ）はサンプリン
グ時点の異なるＬ個の入力信号ｘ（ｔ）、ｘ（ｔ−
１）、…ｘ（ｔ−Ｌ＋１）に対応した変形入力信号ｘ′
（ｔ）、ｘ′（ｔ−１）、…ｘ′（ｔ−Ｌ＋１）、｛但
し、このｘ′（ｔ−ｉ）及びｘ′（ｔ−ｊ）は、 ｜ｘ（ｔ−ｉ）｜≧（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （ｃは１又はこれ以外の定数）が成立する時にｓｇｎ
［ｘ（ｔ−ｉ）］（ここで、ｓｇｎ［ｘ（ｔ−ｉ）］は
ｘ（ｔ−ｉ）＞０の時に＋１、ｘ（ｔ−ｉ）＜０の時は
−１、ｘ（ｔ−ｉ）＝０の時は０を示す。）であり、 ｜ｘ（ｔ−ｉ）｜＜（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にゼロである。］に従って決定することを
特徴とする適応フィルタ装置における係数更新方法に係
わるものである。また、請求項２の発明は、所定のサン
プリング周期で基準信号ｄ(t) を供給するための基準信
号供給手段と、前記基準信号に関係を有する入力信号ｘ
(t) を所定のサンプリング周期で供給するための入力信
号供給手段と、前記基準信号供給手段と前記入力信号供
給手段に接続され、１サンプリング周期を単位遅延時間
として前記入力信号を順次に遅延することによってサン
プリング時点の異なる０番目からＬ−１番目までのＬ個
（但し、Ｌは２よりも大きい正の整数）の入力信号ｘ
（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋１）を同時に
得、前記０番目からＬ番目の入力信号に０番目からＬ−
１番目の係数ｗ0 (t) 、ｗ1(t) 、…ｗL-1 (t) を夫々
乗算し、これによって得られたＬ個の乗算出力を加算し
た出力ｙ(t) を得、前記出力ｙ(t) と前記基準信号ｄ
(t) との誤差信号ｅ(t)を求め、前記入力信号ｘ(t) と
前記誤差信号ｅ(t) とに基づいて前記０番目からＬ−１
番目の係数の更新値ｗ0 （ｔ＋１）、ｗ1 （ｔ＋１）…
ｗL-1 （ｔ＋１）を決定するためのディジタルフィルタ
信号処理回路とを有する適応フィルタ装置における係数
更新方法であって、前記更新係数ｗ0 （ｔ＋１）、ｗ1
（ｔ＋１) 、…ｗL-1 （ｔ＋１）を、 ｗi （ｔ＋１）＝ｗi （ｔ）＋βｅ（ｔ）x ′（ｔ−
ｉ）ｘ ［ここで、各記号は次のものを示す。ｉは０からＬ−１
までの正の整数；ｔはサンプリング時刻を示す記号；ｗ
i （ｔ＋１）はＬ個の更新係数ｗ0 （ｔ＋１）、ｗ1
（ｔ＋１）、…ｗL-1（ｔ＋１）；ｗi (t) はＬ個の更
新前の係数ｗ0 (t) 、ｗ1 (t) 、…ｗL-1 (t) ；βは任
意の定数；ｅ(t) はサンプリング時刻ｔにおける誤差信
号；ｘ′（ｔ−ｉ）はサンプリング時点の異なるＬ個の
入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋１）
に対応した変形入力信号ｘ′（ｔ）、ｘ′（ｔ−１）、
…ｘ′（ｔ−Ｌ＋１）、｛但し、このｘ′（ｔ−ｉ）
は、 ｜ｘ′（ｔ−ｉ）｜≧（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （ｃは定数、ｊは０からＬ−１までの正の整数）が成立
する時にｓｇｎ［ｘ（ｔ−ｉ）］（ここで、ｓｇｎ［ｘ
（ｔ−ｉ）］はｘ（ｔ−ｉ）＞０の時に＋１、ｘ（ｔ−
ｉ）＜０の時は−１、ｘ（ｔ−ｉ）＝０の時は０を示
す。）であり、 ｜ｘ′（ｔ−ｉ）｜＜（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にゼロである。］に従って決定することを
特徴とする適応フィルタ装置における係数更新方法に係
わるものである。なお、本願において、Σj=0 L-1 は
ｘ′（ｔ−ｊ）ｘ（ｔ−ｊ）のｊを０からＬ−１まで変
化させたＬ個の値の総和を求めることを示し、またΣ
j=0 L-1 ｜ｘ（ｔ−ｊ）｜は｜ｘ（ｔ−ｊ）｜のｊを０
からＬ−１まで変化させたＬ個の値の総和を求めること
を示す。
【０００５】
【発明の作用及び効果】本発明における基本的な演算方
法は符号法と同一であり、符号法と同様にｓｇｎ［ｘ
（ｔ−ｉ）］の符号関数を使用する。しかし、すべての
入力信号に対してｓｇｎ［ｘ（ｔ−ｉ）］を使用しな
い。上記の（１／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−ｊ）｜はサ
ンプリング時点の異なるＬ個の入力信号ｘ（ｔ）、ｘ
（ｔ−１）、…ｘ（ｔ−Ｌ＋１）の絶対値の平均値を示
す。この平均値をしきい値としてこれよりも小さい入力
信号の時にはｘ'(ｔ−i)の値を符号関数Ｓｇ＋［ｘ（ｔ
−ｉ）］を使用して決定しないでゼロとする。このよう
に処理すると、収束速度が従来の符号法よりも速くな
る。また、収束誤差（適応誤差）は従来の正規化ＬＭＳ
よりは大きいが、従来の符号法よりも小さくなる。な
お、計算量は従来の符号法とほぼ同じである。従って、
本発明に従う方法は、収束誤差が従来の正規化ＬＭＳと
従来の符号法の中間程度であってもかまわないが、少な
い計算量と速い収束速度を要求する適応フィルタ装置に
好適であり、例えば低価格であるにも拘らず性能の良い
ノイズキャンセラー等を提供することができる。また、
請求項２の発明は、請求項１の発明と同様な技術思想を
ＬＭＳ（LeastMean Square ：最小２乗平均アルゴリズ
ム）に応用したものであり、請求項１の発明よりも更に
演算を簡略化することができる。
法は符号法と同一であり、符号法と同様にｓｇｎ［ｘ
（ｔ−ｉ）］の符号関数を使用する。しかし、すべての
入力信号に対してｓｇｎ［ｘ（ｔ−ｉ）］を使用しな
い。上記の（１／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−ｊ）｜はサ
ンプリング時点の異なるＬ個の入力信号ｘ（ｔ）、ｘ
（ｔ−１）、…ｘ（ｔ−Ｌ＋１）の絶対値の平均値を示
す。この平均値をしきい値としてこれよりも小さい入力
信号の時にはｘ'(ｔ−i)の値を符号関数Ｓｇ＋［ｘ（ｔ
−ｉ）］を使用して決定しないでゼロとする。このよう
に処理すると、収束速度が従来の符号法よりも速くな
る。また、収束誤差（適応誤差）は従来の正規化ＬＭＳ
よりは大きいが、従来の符号法よりも小さくなる。な
お、計算量は従来の符号法とほぼ同じである。従って、
本発明に従う方法は、収束誤差が従来の正規化ＬＭＳと
従来の符号法の中間程度であってもかまわないが、少な
い計算量と速い収束速度を要求する適応フィルタ装置に
好適であり、例えば低価格であるにも拘らず性能の良い
ノイズキャンセラー等を提供することができる。また、
請求項２の発明は、請求項１の発明と同様な技術思想を
ＬＭＳ（LeastMean Square ：最小２乗平均アルゴリズ
ム）に応用したものであり、請求項１の発明よりも更に
演算を簡略化することができる。
【０００６】
【第１の実施例】次に、図１〜図６を参照して本発明の
第１の実施例に係わる適応フィルタ装置を説明する。図
１は未知のシステムのモデリングを行うための対応フィ
ルタ装置であって、入力信号ｘ(t) を所定のサンプリン
グ周期で供給するための入力端子１に未知の動的システ
ム２が接続されている。未知の動的システム２は等価的
に２次の遅延素子３と、ノイズ源４と、遅延素子３の出
力とノイズ源４のノイズを加算する加算器５とを有す
る。この未知の動的システム２は所定サンプリング周期
で望みの応答即ち基準信号ｄ(t) を出力する。
第１の実施例に係わる適応フィルタ装置を説明する。図
１は未知のシステムのモデリングを行うための対応フィ
ルタ装置であって、入力信号ｘ(t) を所定のサンプリン
グ周期で供給するための入力端子１に未知の動的システ
ム２が接続されている。未知の動的システム２は等価的
に２次の遅延素子３と、ノイズ源４と、遅延素子３の出
力とノイズ源４のノイズを加算する加算器５とを有す
る。この未知の動的システム２は所定サンプリング周期
で望みの応答即ち基準信号ｄ(t) を出力する。
【０００７】入力端子１はトランスバーサル型適応フィ
ルタ６にも接続されている。トランスバーサル型適応フ
ィルタ６は、１サンプリング周期を単位遅延時間として
入力信号ｘ(t) を順次に遅延させるための遅延手段とし
てのＬ−１個（３個）の遅延素子７、８、９と、入力端
子１に接続されたタップライン１０及び遅延素子７、
８、９の出力端子に接続されたタップライン１１、１
２、１３と、各タップライン１０〜１３に接続されたＬ
個（４個）の係数乗算器１４、１５、１６、１７と、各
乗算器１４、１５、１６、１７の出力を順次に加算する
加算器１８、１９、２０と、係数更新手段２１とから成
る。
ルタ６にも接続されている。トランスバーサル型適応フ
ィルタ６は、１サンプリング周期を単位遅延時間として
入力信号ｘ(t) を順次に遅延させるための遅延手段とし
てのＬ−１個（３個）の遅延素子７、８、９と、入力端
子１に接続されたタップライン１０及び遅延素子７、
８、９の出力端子に接続されたタップライン１１、１
２、１３と、各タップライン１０〜１３に接続されたＬ
個（４個）の係数乗算器１４、１５、１６、１７と、各
乗算器１４、１５、１６、１７の出力を順次に加算する
加算器１８、１９、２０と、係数更新手段２１とから成
る。
【０００８】最初のタップライン１０にはサンプリング
時刻ｔにおける入力信号ｘ(t) がそのまま送られる。遅
延素子７、８、９の出力端子には１つ前のサンプリング
時刻（ｔ−１）の入力信号ｘ（ｔ−１）、２つ前のサン
プリング時刻（ｔ−２）の入力信号ｘ（ｔ−２）、３つ
前のサンプリング時刻（ｔ−３）の入力信号ｘ（ｔ−
３）が得られる。タップラインが全部でＬ個の場合に
は、Ｌ番目のラインにｘ（ｔ−Ｌ＋１）の信号が得られ
る。
時刻ｔにおける入力信号ｘ(t) がそのまま送られる。遅
延素子７、８、９の出力端子には１つ前のサンプリング
時刻（ｔ−１）の入力信号ｘ（ｔ−１）、２つ前のサン
プリング時刻（ｔ−２）の入力信号ｘ（ｔ−２）、３つ
前のサンプリング時刻（ｔ−３）の入力信号ｘ（ｔ−
３）が得られる。タップラインが全部でＬ個の場合に
は、Ｌ番目のラインにｘ（ｔ−Ｌ＋１）の信号が得られ
る。
【０００９】乗算器１４、１５、１６、１７は係数更新
手段２１に接続されており、ここから与えられた係数ｗ
0(t)、ｗ1(t)、ｗ2(t)、ｗ3(t)をサンプリング時刻の異
なる４つの入力信号ｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−
２）、ｘ（ｔ−３）にそれぞれ乗算し、ｗ0(t)ｘ(t) 、
ｗ1(t)ｘ（ｔ−１）、ｗ2(t)ｘ（ｔ−２）、ｗ3(t)ｘ
（ｔ−３）を送出する。
手段２１に接続されており、ここから与えられた係数ｗ
0(t)、ｗ1(t)、ｗ2(t)、ｗ3(t)をサンプリング時刻の異
なる４つの入力信号ｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−
２）、ｘ（ｔ−３）にそれぞれ乗算し、ｗ0(t)ｘ(t) 、
ｗ1(t)ｘ（ｔ−１）、ｗ2(t)ｘ（ｔ−２）、ｗ3(t)ｘ
（ｔ−３）を送出する。
【００１０】３つの加算器１８、１９、２０は４つの乗
算器１４、１５、１６、１７の出力の総和を求め、フィ
ルタ出力として ｙ(t) ＝ｗ0(t)ｘ（ｔ）＋ｗ1(t)ｘ（ｔ−１）＋ｗ2(t)ｘ（ｔ−２）＋ｗ3(t) ｘ（ｔ−３） 式（１） を出力する。
算器１４、１５、１６、１７の出力の総和を求め、フィ
ルタ出力として ｙ(t) ＝ｗ0(t)ｘ（ｔ）＋ｗ1(t)ｘ（ｔ−１）＋ｗ2(t)ｘ（ｔ−２）＋ｗ3(t) ｘ（ｔ−３） 式（１） を出力する。
【００１１】未知の動的システム２と加算器２０に接続
された誤差演算回路２２はｄ(t) −ｙ(t) の演算を実行
して誤差（エラー）信号ｅ(t) を出力する。
された誤差演算回路２２はｄ(t) −ｙ(t) の演算を実行
して誤差（エラー）信号ｅ(t) を出力する。
【００１２】係数更新手段２１は、係数決定に使用する
情報を得るために誤差演算回路２２と入力端子１と各遅
延素子７、８、９の出力端子に接続されている。
情報を得るために誤差演算回路２２と入力端子１と各遅
延素子７、８、９の出力端子に接続されている。
【００１３】図１ではトランスバーサル型適応フィルタ
６の理解を容易にするために、これを個別回路によって
機能的（等価的）に示したが、実際にはフィルタ６及び
誤差演算回路２２はディジタル信号処理回路又はマイク
ロプロセッサ（マイクロコンピュータ）で構成されてい
る。
６の理解を容易にするために、これを個別回路によって
機能的（等価的）に示したが、実際にはフィルタ６及び
誤差演算回路２２はディジタル信号処理回路又はマイク
ロプロセッサ（マイクロコンピュータ）で構成されてい
る。
【００１４】係数更新手段は、タップライン１０〜１３
及び乗算器１４〜１７を一般的Ｌ個とした場合におい
て、次式を演算するように構成されている。 ｗi （ｔ＋１）＝ｗi (t) ＋β［ｅ(t) ／｛Σj=0 L-1 x ′（ｔ−ｊ）ｘ（ ｔ−ｊ）］ｘ′（ｔ−ｊ）］ｘ′（ｔ−ｉ） 式（２） この式（２）の各記号は特許請求の範囲及び課題を解決
する手段の欄において定義した通りである。なお、実施
例ではｉ及びｊは０、１、２、３の値をとる。
及び乗算器１４〜１７を一般的Ｌ個とした場合におい
て、次式を演算するように構成されている。 ｗi （ｔ＋１）＝ｗi (t) ＋β［ｅ(t) ／｛Σj=0 L-1 x ′（ｔ−ｊ）ｘ（ ｔ−ｊ）］ｘ′（ｔ−ｊ）］ｘ′（ｔ−ｉ） 式（２） この式（２）の各記号は特許請求の範囲及び課題を解決
する手段の欄において定義した通りである。なお、実施
例ではｉ及びｊは０、１、２、３の値をとる。
【００１５】式（２）を適応フィルタの分野で通常使用
されているベクトル形式で示すと次式になる。 Ｗ（ｔ＋１）＝Ｗ(t) ＋β［ｅ(t) ／｛Ｘ′T (t) Ｘ(t) ｝］Ｘ′(t) 式（３） ここで、ｔはサンプリング時刻を示す記号、Ｗ(t) はサ
ンプリング時刻ｔにおけるＬ個の係数を示す係数ベクト
ル、Ｗ（ｔ＋１）はサンプリング時刻ｔ＋１におけるＬ
個の係数を示す係数ベクトル、βは任意の定数、ｅ(t)
はサンプリング時刻ｔにおける前記誤差信号、Ｘ(t) は
サンプリング時点の異なるＬ個の入力信号を示すベクト
ル、Ｘ′(t) はサンプリング時点の異なるＬ個の入力信
号｛ｘ（ｔ）、ｘ（ｔ−１）…ｘ（ｔ−Ｌ＋１）｝に対
応したＬ個の変形入力信号｛ｘ′（ｔ）、ｘ′（ｔ−
１）…ｘ′（ｔ−Ｌ＋１）｝を示すベクトル｛但し、前
記Ｌ個の変形入力信号から選ばれた任意の変形入力信号
をｘ′（ｔ−ｉ）（但し、ｉは０からＬ−１までの正の
整数である。）で表わした場合において、 ｜ｘ（ｔ−ｉ）｜≧（１／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （但し、ここでｊは０からＬ−１までの数値である。）
が成立する時には、 ｘ′（ｔ−ｉ）＝ｓｇｎ［ｘ（ｔ−ｉ）］ （ここで、ｓｇｎ［ｘ（ｔ−ｉ）］は、ｘ（ｔ−ｉ）＞
０の時に＋１、ｘ（ｔ−ｉ）＜０の時に−１、ｘ（ｔ−
ｉ）＝０の時に０の符号を示す。）であり、 ｜ｘ（ｔ−ｉ）｜＜（１／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にはｘ′（ｔ−ｉ）＝０である。｝Ｘ′T
(t) はベクトルＸ′(t) の転置ベクトルである。
されているベクトル形式で示すと次式になる。 Ｗ（ｔ＋１）＝Ｗ(t) ＋β［ｅ(t) ／｛Ｘ′T (t) Ｘ(t) ｝］Ｘ′(t) 式（３） ここで、ｔはサンプリング時刻を示す記号、Ｗ(t) はサ
ンプリング時刻ｔにおけるＬ個の係数を示す係数ベクト
ル、Ｗ（ｔ＋１）はサンプリング時刻ｔ＋１におけるＬ
個の係数を示す係数ベクトル、βは任意の定数、ｅ(t)
はサンプリング時刻ｔにおける前記誤差信号、Ｘ(t) は
サンプリング時点の異なるＬ個の入力信号を示すベクト
ル、Ｘ′(t) はサンプリング時点の異なるＬ個の入力信
号｛ｘ（ｔ）、ｘ（ｔ−１）…ｘ（ｔ−Ｌ＋１）｝に対
応したＬ個の変形入力信号｛ｘ′（ｔ）、ｘ′（ｔ−
１）…ｘ′（ｔ−Ｌ＋１）｝を示すベクトル｛但し、前
記Ｌ個の変形入力信号から選ばれた任意の変形入力信号
をｘ′（ｔ−ｉ）（但し、ｉは０からＬ−１までの正の
整数である。）で表わした場合において、 ｜ｘ（ｔ−ｉ）｜≧（１／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （但し、ここでｊは０からＬ−１までの数値である。）
が成立する時には、 ｘ′（ｔ−ｉ）＝ｓｇｎ［ｘ（ｔ−ｉ）］ （ここで、ｓｇｎ［ｘ（ｔ−ｉ）］は、ｘ（ｔ−ｉ）＞
０の時に＋１、ｘ（ｔ−ｉ）＜０の時に−１、ｘ（ｔ−
ｉ）＝０の時に０の符号を示す。）であり、 ｜ｘ（ｔ−ｉ）｜＜（１／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にはｘ′（ｔ−ｉ）＝０である。｝Ｘ′T
(t) はベクトルＸ′(t) の転置ベクトルである。
【００１６】図１の４個の係数 ｗ0 (t) 、ｗ1 (t) 、
ｗ2 (t) 、ｗ3 (t) の更新係数を式（２）に従って示す
と次の通りである。 ｗ0 （ｔ＋１）＝ｗ0 (t) ＋β［ｅ(t) ／Ａ］ｘ0 (t) ｗ1 （ｔ＋１）＝ｗ1 (t) ＋β［ｅ(t) ／Ａ］ｘ1 (t−１) ｗ2 （ｔ＋１）＝ｗ2 (t) ＋β［ｅ(t) ／Ａ］ｘ2 (t−２) ｗ3 （ｔ＋１）＝ｗ3 (t) ＋β［ｅ(t) ／Ａ］ｘ3 (t−３) 式（４） 但し、式（４）におけるＡは｛Σj=0 L-1 ｘ′（ｔ−
ｊ）ｘ（ｔ−ｊ）｝を示す。
ｗ2 (t) 、ｗ3 (t) の更新係数を式（２）に従って示す
と次の通りである。 ｗ0 （ｔ＋１）＝ｗ0 (t) ＋β［ｅ(t) ／Ａ］ｘ0 (t) ｗ1 （ｔ＋１）＝ｗ1 (t) ＋β［ｅ(t) ／Ａ］ｘ1 (t−１) ｗ2 （ｔ＋１）＝ｗ2 (t) ＋β［ｅ(t) ／Ａ］ｘ2 (t−２) ｗ3 （ｔ＋１）＝ｗ3 (t) ＋β［ｅ(t) ／Ａ］ｘ3 (t−３) 式（４） 但し、式（４）におけるＡは｛Σj=0 L-1 ｘ′（ｔ−
ｊ）ｘ（ｔ−ｊ）｝を示す。
【００１７】図２は式（４）の演算を実行するための図
１の係数更新手段２１を機能的（等価的）に示すブロッ
ク図である。この係数更新手段２１は、図１の入力端子
１及び遅延素子７、８、９の出力端子に接続される４つ
の入力端子３１、３２、３３、３４と誤差演算回路２２
に接続されるエラー信号入力端子３５とを有する。４つ
の入力端子３１、３２、３３、３４に接続された平均値
算出回路３６は、入力端子３１〜３４の信号ｘ（ｔ）、
ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ（ｔ−３）を入力とし
てこれ等の絶対値の平均値を算出する。即ち、（１／
Ｌ）Σj=0 L-1 ｜ｘ（ｔ−ｊ）｜＝ｘa (t) 演算を実行
する。なお、１／Ｌの代わりにｃ／Ｌを使用し、ｃを任
意の定数とすることもできる。
１の係数更新手段２１を機能的（等価的）に示すブロッ
ク図である。この係数更新手段２１は、図１の入力端子
１及び遅延素子７、８、９の出力端子に接続される４つ
の入力端子３１、３２、３３、３４と誤差演算回路２２
に接続されるエラー信号入力端子３５とを有する。４つ
の入力端子３１、３２、３３、３４に接続された平均値
算出回路３６は、入力端子３１〜３４の信号ｘ（ｔ）、
ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ（ｔ−３）を入力とし
てこれ等の絶対値の平均値を算出する。即ち、（１／
Ｌ）Σj=0 L-1 ｜ｘ（ｔ−ｊ）｜＝ｘa (t) 演算を実行
する。なお、１／Ｌの代わりにｃ／Ｌを使用し、ｃを任
意の定数とすることもできる。
【００１８】入力端子３１、３２、３３、３４と平均値
算出回路３６とに接続された４つの判定回路３７、３
８、３９、４０は入力信号ｘ（ｔ）、ｘ（ｔ−１）、ｘ
（ｔ−２）、ｘ（ｔ−３）が平均値ｘa （ｔ）以上か、
又は小さいかを判定し、更に、入力信号が平均値以上の
場合には入力信号が正の値か、又は負の値か又はゼロか
を判定し、＋１又は−１又は０の符号を発生する。即
ち、判定回路３７〜４０は、ｘ0 ′(t) 、ｘ1 ′(t) 、
ｘ2 ′(t) 、ｘ3 '(t) を決定する回路であり、 ｘ（ｔ）≧ｘa （ｔ） ｘ（ｔ−１）≧ｘa （ｔ） ｘ（ｔ−２）≧ｘa （ｔ） ｘ（ｔ−３）≧ｘa （ｔ） であり、且つｘ（ｔ−ｉ）が正の値の時には＋１を発生
し、ｘ（ｔ−ｉ）が負の値の時には−１を発生し、ｘ
（ｔ−ｉ）がゼロの時にはゼロを発生する。また、ｘ
（ｔ−ｉ）＜ｘa (t) の時にはゼロを発生する。
算出回路３６とに接続された４つの判定回路３７、３
８、３９、４０は入力信号ｘ（ｔ）、ｘ（ｔ−１）、ｘ
（ｔ−２）、ｘ（ｔ−３）が平均値ｘa （ｔ）以上か、
又は小さいかを判定し、更に、入力信号が平均値以上の
場合には入力信号が正の値か、又は負の値か又はゼロか
を判定し、＋１又は−１又は０の符号を発生する。即
ち、判定回路３７〜４０は、ｘ0 ′(t) 、ｘ1 ′(t) 、
ｘ2 ′(t) 、ｘ3 '(t) を決定する回路であり、 ｘ（ｔ）≧ｘa （ｔ） ｘ（ｔ−１）≧ｘa （ｔ） ｘ（ｔ−２）≧ｘa （ｔ） ｘ（ｔ−３）≧ｘa （ｔ） であり、且つｘ（ｔ−ｉ）が正の値の時には＋１を発生
し、ｘ（ｔ−ｉ）が負の値の時には−１を発生し、ｘ
（ｔ−ｉ）がゼロの時にはゼロを発生する。また、ｘ
（ｔ−ｉ）＜ｘa (t) の時にはゼロを発生する。
【００１９】入力端子３１、３２、３３、３４と判定回
路３７、３８、３９、４０にそれぞれ接続された４つの
乗算器４１、４２、４３、４４は式（２）のｘ′（ｔ−
ｉ）ｘ（ｔ−ｉ）の演算を実行するものであり、ｘ′
(t) ｘ(t) 、ｘ′（ｔ−１）ｘ（ｔ−１）、ｘ′（ｔ−
２）ｘ（ｔ−２）、ｘ′（ｔ−３）ｘ（ｔ−３）を出力
する。なお、ｘ′（ｔ）、ｘ′（ｔ−１）、ｘ′（ｔ−
２）、ｘ′（ｔ−３）は＋１又は−１又は０であるの
で、この乗算は簡単に行うことができる。
路３７、３８、３９、４０にそれぞれ接続された４つの
乗算器４１、４２、４３、４４は式（２）のｘ′（ｔ−
ｉ）ｘ（ｔ−ｉ）の演算を実行するものであり、ｘ′
(t) ｘ(t) 、ｘ′（ｔ−１）ｘ（ｔ−１）、ｘ′（ｔ−
２）ｘ（ｔ−２）、ｘ′（ｔ−３）ｘ（ｔ−３）を出力
する。なお、ｘ′（ｔ）、ｘ′（ｔ−１）、ｘ′（ｔ−
２）、ｘ′（ｔ−３）は＋１又は−１又は０であるの
で、この乗算は簡単に行うことができる。
【００２０】乗算器４１〜４４に接続された３つの加算
器４５、４６、４７は式（４）のＡを求めるものであ
り、Ａとしてｘ′(t) ｘ(t) ＋ｘ′（ｔ−１）（ｔ−
１）＋ｘ′（ｔ−２）ｘ（ｔ−２）＋ｘ′（ｔ−３）ｘ
（ｔ−３）を除算器４８に送る。
器４５、４６、４７は式（４）のＡを求めるものであ
り、Ａとしてｘ′(t) ｘ(t) ＋ｘ′（ｔ−１）（ｔ−
１）＋ｘ′（ｔ−２）ｘ（ｔ−２）＋ｘ′（ｔ−３）ｘ
（ｔ−３）を除算器４８に送る。
【００２１】エラー信号入力端子３５に接続されたβ乗
算器４９はエラー信号ｅ(t) に定数βを乗算して除算器
４８に送る。βは例えば０．５である。
算器４９はエラー信号ｅ(t) に定数βを乗算して除算器
４８に送る。βは例えば０．５である。
【００２２】除算器４８は式（２）のβ［ｅ(t) ／Ａ］
を求めて乗算器５０、５１、５２、５３に送る。
を求めて乗算器５０、５１、５２、５３に送る。
【００２３】乗算器５０、５１、５２、５３は判定回路
３７、３８、３９、４０に接続されていると共に除算器
４８に接続され、式（２）の右辺の第２項のｘ′（ｔ−
ｉ）の乗算を行って、 ｛βｅ(t) ／Ａ｝ｘ′（ｔ） ｛βｅ(t) ／Ａ｝ｘ′（ｔ−１） ｛βｅ(t) ／Ａ｝ｘ′（ｔ−２） ｛βｅ(t) ／Ａ｝ｘ′（ｔ−３） を出力する。
３７、３８、３９、４０に接続されていると共に除算器
４８に接続され、式（２）の右辺の第２項のｘ′（ｔ−
ｉ）の乗算を行って、 ｛βｅ(t) ／Ａ｝ｘ′（ｔ） ｛βｅ(t) ／Ａ｝ｘ′（ｔ−１） ｛βｅ(t) ／Ａ｝ｘ′（ｔ−２） ｛βｅ(t) ／Ａ｝ｘ′（ｔ−３） を出力する。
【００２４】４個のメモリ５４、５５、５６、５７は現
在時点ｔの係数を記憶するものである。即ち、式（２）
のｗi (t) に対応するｗ0 (t) 、ｗ1 (t) 、ｗ2 (t) 、
ｗ3(t) を記憶するものである。
在時点ｔの係数を記憶するものである。即ち、式（２）
のｗi (t) に対応するｗ0 (t) 、ｗ1 (t) 、ｗ2 (t) 、
ｗ3(t) を記憶するものである。
【００２５】４つの加算器５８、５９、６０、６１は乗
算器５０、５１、５２、５３の出力とメモリ５４、５
５、５６、５７の出力とに接続され、これ等の出力の加
算値を出力端子６２、６３、６４、６５に送出するもの
である。４つの出力端子６２、６３、６４、６５には式
（２）の更新係数ｗi （ｔ＋１）に対応するｗ0 （ｔ＋
１）、ｗ1 （ｔ＋１）、ｗ2 （ｔ＋２）、ｗ3 （ｔ＋
３）が得られ、これ等が図１の乗算器１４、１５、１
６、１７に次のサンプリング時点（ｔ＋１）の係数とし
て送られる。
算器５０、５１、５２、５３の出力とメモリ５４、５
５、５６、５７の出力とに接続され、これ等の出力の加
算値を出力端子６２、６３、６４、６５に送出するもの
である。４つの出力端子６２、６３、６４、６５には式
（２）の更新係数ｗi （ｔ＋１）に対応するｗ0 （ｔ＋
１）、ｗ1 （ｔ＋１）、ｗ2 （ｔ＋２）、ｗ3 （ｔ＋
３）が得られ、これ等が図１の乗算器１４、１５、１
６、１７に次のサンプリング時点（ｔ＋１）の係数とし
て送られる。
【００２６】図３は係数更新手段２１における係数更新
処理の所定のプログラムに従って実行する手順を示す。
まず、ブロック７０でスタートすると、ブロック７１で
入力信号ｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ
（ｔ−３）を得る。次に、ブロック７２で各入力信号の
絶対値の平均値ｘa （ｔ）を求める。。
処理の所定のプログラムに従って実行する手順を示す。
まず、ブロック７０でスタートすると、ブロック７１で
入力信号ｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ
（ｔ−３）を得る。次に、ブロック７２で各入力信号の
絶対値の平均値ｘa （ｔ）を求める。。
【００２７】次に、ブロック７３に示すように入力信号
の絶対値｜ｘ(t) ｜が平均値ｘa (ｔ) 以上か否かを判
定する。この値が平均値ｘa (t) よりも大きい時にはブ
ロック７４に示すようにｘ(t) がゼロよりも大きいか小
さいかを判断する。ｘ(t) ＞０の時にはブロック７５に
示すようにｘ′(t) として＋１を出力し、ｘ(t) ＜０の
時にはブロック７６に示すようにｘ′(t) として−１を
出力する。一方、ブロック７３において｜ｘ(t) ｜が平
均値ｘa （t)未満であることが判定された時にはブロッ
ク７７に示すようにｘ′(t) としてゼロを出力する。な
お、ｘ′(t) の値はメモリで一時保持する。
の絶対値｜ｘ(t) ｜が平均値ｘa (ｔ) 以上か否かを判
定する。この値が平均値ｘa (t) よりも大きい時にはブ
ロック７４に示すようにｘ(t) がゼロよりも大きいか小
さいかを判断する。ｘ(t) ＞０の時にはブロック７５に
示すようにｘ′(t) として＋１を出力し、ｘ(t) ＜０の
時にはブロック７６に示すようにｘ′(t) として−１を
出力する。一方、ブロック７３において｜ｘ(t) ｜が平
均値ｘa （t)未満であることが判定された時にはブロッ
ク７７に示すようにｘ′(t) としてゼロを出力する。な
お、ｘ′(t) の値はメモリで一時保持する。
【００２８】ブロック７８、７９、８０に示すように、
ｘ′（ｔ−１）、ｘ′（ｔ−２）、ｘ′（ｔ−３）の値
をブロック７３〜７７と同一の方法で順次に求める。
ｘ′（ｔ−１）、ｘ′（ｔ−２）、ｘ′（ｔ−３）の値
をブロック７３〜７７と同一の方法で順次に求める。
【００２９】次に、図４のブロック８１に示すように
ｘ′(t) ＝１か否かを判定する。ｘ′(t) ＝１の時には
ｘ′(t) ｘ(t) の値としてブロック８２に示すようにｘ
(t) を得る。一方、ｘ′(t) ＝１でない場合にはブロッ
ク８３においてｘ′(t) ＝−１か否かを判定し、ＹＥＳ
の場合にはブロック８４に示すようにｘ′(t) ｘ(t) と
して−ｘ(t) を得る。また、ブロック８３の出力がＮＯ
の場合はブロック８５に示すようにｘ′(t) ｘ(t) とし
てゼロを得る。
ｘ′(t) ＝１か否かを判定する。ｘ′(t) ＝１の時には
ｘ′(t) ｘ(t) の値としてブロック８２に示すようにｘ
(t) を得る。一方、ｘ′(t) ＝１でない場合にはブロッ
ク８３においてｘ′(t) ＝−１か否かを判定し、ＹＥＳ
の場合にはブロック８４に示すようにｘ′(t) ｘ(t) と
して−ｘ(t) を得る。また、ブロック８３の出力がＮＯ
の場合はブロック８５に示すようにｘ′(t) ｘ(t) とし
てゼロを得る。
【００３０】次に、ブロック８６、８７、８８に示すよ
うに、ｘ′（ｔ−１）ｘ（ｔ−１）の値、ｘ′（ｔ−
２）ｘ（ｔ−２）の値、ｘ′（ｔ−３）ｘ（ｔ−３）の
値を、ブロック８１〜８５と同様な方法で求める。
うに、ｘ′（ｔ−１）ｘ（ｔ−１）の値、ｘ′（ｔ−
２）ｘ（ｔ−２）の値、ｘ′（ｔ−３）ｘ（ｔ−３）の
値を、ブロック８１〜８５と同様な方法で求める。
【００３１】次に、ブロック８１〜８５で求めたｘ′
(t) ｘ(t) の値とブロック８６、８７、８８で求めた値
とをブロック８９に示すように加算して値Ａを求める。
(t) ｘ(t) の値とブロック８６、８７、８８で求めた値
とをブロック８９に示すように加算して値Ａを求める。
【００３２】次に、ブロック９０に示すようにβｅ(t)
／Ａを求める。次に、ブロック９１に示すようにβｅ
(t) ／Ａに対してｘ′（ｔ）、ｘ′（ｔ−１）、ｘ′
（ｔ−２）、ｘ′（ｔ−３）の値（＋１、又は−１又は
０）を乗算してβｅ(t) ／Ａ｛ｘ′(t) ｝、βｅ(t) ／
Ａ｛ｘ′（ｔ−１）｝、βｅ(t) ／Ａ｛ｘ′（ｔ−
２）｝、βｅ(t) ／Ａ｛ｘ′（ｔ−３）｝求める。
／Ａを求める。次に、ブロック９１に示すようにβｅ
(t) ／Ａに対してｘ′（ｔ）、ｘ′（ｔ−１）、ｘ′
（ｔ−２）、ｘ′（ｔ−３）の値（＋１、又は−１又は
０）を乗算してβｅ(t) ／Ａ｛ｘ′(t) ｝、βｅ(t) ／
Ａ｛ｘ′（ｔ−１）｝、βｅ(t) ／Ａ｛ｘ′（ｔ−
２）｝、βｅ(t) ／Ａ｛ｘ′（ｔ−３）｝求める。
【００３３】次に、ブロック９２に示すようにブロック
９１で求めた値にｗ0 （ｔ）、ｗ0（ｔ−１）、ｗ0
（ｔ−２）、ｗ0 （ｔ−３）を加算して更新係数ｗ0
（ｔ＋１）、ｗ1 （ｔ＋１）、ｗ2 （ｔ＋１）、ｗ3
（ｔ＋１）を求め、ブロック９３に示すようにメモリに
格納し、図１の乗算器１４〜１７におけるｔ＋１時点の
新しい係数として使用する。その後、図３のブロック７
１に戻り、次のサンプリング周期の動作を繰返す。
９１で求めた値にｗ0 （ｔ）、ｗ0（ｔ−１）、ｗ0
（ｔ−２）、ｗ0 （ｔ−３）を加算して更新係数ｗ0
（ｔ＋１）、ｗ1 （ｔ＋１）、ｗ2 （ｔ＋１）、ｗ3
（ｔ＋１）を求め、ブロック９３に示すようにメモリに
格納し、図１の乗算器１４〜１７におけるｔ＋１時点の
新しい係数として使用する。その後、図３のブロック７
１に戻り、次のサンプリング周期の動作を繰返す。
【００３４】上述から明らかなように、本実施例ではサ
ンプリング時点の異なる入力信号の絶対値の平均値ｘa
(t) を求めるが、この平均値ｘa (t)の計算は比較的簡
単に達成される。即ち、平均値を計算する場合には４つ
の入力端子３１、３２、３３、３４の合計を求めること
が必要になるが、１サンプリング時点毎に４つの入力信
号ｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ（ｔ−
３）の内の最も古い１つの入力信号ｘ（ｔ−３）が次の
新しい入力信号に入れ替わるのみであるから、この平均
値計算は簡単に達成される。また、ｘ′（ｔ）、ｘ′
（ｔ−１）、ｘ′（ｔ−２）、ｘ′（ｔ−３）は＋１又
は−１又は０であるので、式（２）の右辺の演算を容易
に達成することができる。即ち、本発明に従う方法の係
数更新の演算量は従来の符号法におけるそれとほぼ同等
である。
ンプリング時点の異なる入力信号の絶対値の平均値ｘa
(t) を求めるが、この平均値ｘa (t)の計算は比較的簡
単に達成される。即ち、平均値を計算する場合には４つ
の入力端子３１、３２、３３、３４の合計を求めること
が必要になるが、１サンプリング時点毎に４つの入力信
号ｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ（ｔ−
３）の内の最も古い１つの入力信号ｘ（ｔ−３）が次の
新しい入力信号に入れ替わるのみであるから、この平均
値計算は簡単に達成される。また、ｘ′（ｔ）、ｘ′
（ｔ−１）、ｘ′（ｔ−２）、ｘ′（ｔ−３）は＋１又
は−１又は０であるので、式（２）の右辺の演算を容易
に達成することができる。即ち、本発明に従う方法の係
数更新の演算量は従来の符号法におけるそれとほぼ同等
である。
【００３５】次に、本発明の方法と従来の正規化ＬＭＳ
及び符号法との関係を更に詳しく説明する。式（３）に
対応させて従来の正規化ＬＭＳの係数更新のための演算
式を示すと次の通りになる。 Ｗ（ｔ＋１）＝Ｗ(t) ＋β［ｅ(t) ／｛‖Ｘ(t) ‖2 ｝
Ｘ(t) この演算式では‖Ｘ(t) ‖2 の演算があるので、演算量
が本発明の方法よりも多い。また、従来の符号法の係数
更新のための演算式を示すと次の通りになる。 Ｗ（ｔ＋１）＝Ｗ(t) ＋β［ｅ(t) ／｛Ｘ′T (t) Ｘ
(t) ｝Ｘ′(t) ］ 但し、ここで、Ｘ′(t) ＝ｓｇｎｘ（ｔ−ｉ）であり、
ｘ（ｔ−ｉ）＞０の時に＋１、ｘ（ｔ−ｉ）＜０の時に
−１、ｘ（ｔ−ｉ）＝０の時に０である。符号法と本発
明の方法とを比べると、本発明では平均値ｘa （t)を求
め、この平均値ｘa (t) よりも入力信号が小さい時には
Ｘ′(t) をゼロとするステップを新たに導入した点で両
者は異なる。従って、演算量においては前述したように
本発明と符号法とで大差ない。
及び符号法との関係を更に詳しく説明する。式（３）に
対応させて従来の正規化ＬＭＳの係数更新のための演算
式を示すと次の通りになる。 Ｗ（ｔ＋１）＝Ｗ(t) ＋β［ｅ(t) ／｛‖Ｘ(t) ‖2 ｝
Ｘ(t) この演算式では‖Ｘ(t) ‖2 の演算があるので、演算量
が本発明の方法よりも多い。また、従来の符号法の係数
更新のための演算式を示すと次の通りになる。 Ｗ（ｔ＋１）＝Ｗ(t) ＋β［ｅ(t) ／｛Ｘ′T (t) Ｘ
(t) ｝Ｘ′(t) ］ 但し、ここで、Ｘ′(t) ＝ｓｇｎｘ（ｔ−ｉ）であり、
ｘ（ｔ−ｉ）＞０の時に＋１、ｘ（ｔ−ｉ）＜０の時に
−１、ｘ（ｔ−ｉ）＝０の時に０である。符号法と本発
明の方法とを比べると、本発明では平均値ｘa （t)を求
め、この平均値ｘa (t) よりも入力信号が小さい時には
Ｘ′(t) をゼロとするステップを新たに導入した点で両
者は異なる。従って、演算量においては前述したように
本発明と符号法とで大差ない。
【００３６】本発明に従う適応フィルタの効果を確認す
るために、図１の装置の入力端子１に白色雑音を加え、
ノイズ源４からも白色雑音を加えて収束特性を調べたと
ころ、図５に示す結果が得られた。図５において、Ｃ1
が本発明に従う方法の特性、Ｃ2 が従来の正規化ＬＭＳ
の特性、Ｃ3 が従来の符号法の特性を示す。この結果か
ら明らかなように、本発明によれば収束性即ち推定エラ
ー（収束誤差）が正規化ＬＭＳとほぼ同等となり、従来
の符号法よりは僅かに改善される。なお、図５において
横軸はサンプル番号を示し、縦軸はＥ＝｛‖Ｈ−ｗ(t)
‖2 ｝／‖Ｈ‖2 で示される推定エラーを示す。但し、
上式におけるＨは未知パラメータを表わし、Ｗはその推
定値を表わし、推定値の初期値はゼロである。また、図
１の入力端子１に白色ノイズを１／（１−０．９５
ｚ-1）で示されるＩＩＲ型ディジタルフィルタを通した
得た有色ノイズを加えて図５と同様に収束性を調べたと
ころ、図６に示す結果が得られた。図６から明らかなよ
うに、本発明の方法は入力信号が有色雑音の場合に従来
の符号法よりも収束性が優れている。
るために、図１の装置の入力端子１に白色雑音を加え、
ノイズ源４からも白色雑音を加えて収束特性を調べたと
ころ、図５に示す結果が得られた。図５において、Ｃ1
が本発明に従う方法の特性、Ｃ2 が従来の正規化ＬＭＳ
の特性、Ｃ3 が従来の符号法の特性を示す。この結果か
ら明らかなように、本発明によれば収束性即ち推定エラ
ー（収束誤差）が正規化ＬＭＳとほぼ同等となり、従来
の符号法よりは僅かに改善される。なお、図５において
横軸はサンプル番号を示し、縦軸はＥ＝｛‖Ｈ−ｗ(t)
‖2 ｝／‖Ｈ‖2 で示される推定エラーを示す。但し、
上式におけるＨは未知パラメータを表わし、Ｗはその推
定値を表わし、推定値の初期値はゼロである。また、図
１の入力端子１に白色ノイズを１／（１−０．９５
ｚ-1）で示されるＩＩＲ型ディジタルフィルタを通した
得た有色ノイズを加えて図５と同様に収束性を調べたと
ころ、図６に示す結果が得られた。図６から明らかなよ
うに、本発明の方法は入力信号が有色雑音の場合に従来
の符号法よりも収束性が優れている。
【００３７】本発明に従う方法、従来の正規化ＬＭＳ、
従来の符号法のそれぞれについて、３２回の測定におけ
る収束後に適応誤差の２乗平均値を求めたところ、白色
ノイズの場合は、０．０１３、 ０．０１０、 ０．０
２２であり、有色ノイズの場合は０．０１２、 ０．０
１１、 ０．０１８であった。この結果から明らかなよ
うに適応誤差においては、本発明の方法は従来の正規化
ＬＭＳとほぼ同等であり、従来の符号法よりは優れてい
る。
従来の符号法のそれぞれについて、３２回の測定におけ
る収束後に適応誤差の２乗平均値を求めたところ、白色
ノイズの場合は、０．０１３、 ０．０１０、 ０．０
２２であり、有色ノイズの場合は０．０１２、 ０．０
１１、 ０．０１８であった。この結果から明らかなよ
うに適応誤差においては、本発明の方法は従来の正規化
ＬＭＳとほぼ同等であり、従来の符号法よりは優れてい
る。
【００３８】
【第２の実施例】次に、図７を参照して、本発明の第２
の実施例に係わる係数更新手段及び方法を説明する。第
２の実施例の適応フィルタ装置の基本的構成は図１と同
一であり、係数更新手段２１が図７に示すように変形さ
れている。図７の係数更新手段において、図２と共通す
る部分には同一の符号を付してその説明を省略する。図
２と図７との比較から明らかなように、図７の回路は、
図２の回路から乗算器４１、４２、４３、４４、加算器
４５、４６、４７、除算器４８を省いたものに相当す
る。
の実施例に係わる係数更新手段及び方法を説明する。第
２の実施例の適応フィルタ装置の基本的構成は図１と同
一であり、係数更新手段２１が図７に示すように変形さ
れている。図７の係数更新手段において、図２と共通す
る部分には同一の符号を付してその説明を省略する。図
２と図７との比較から明らかなように、図７の回路は、
図２の回路から乗算器４１、４２、４３、４４、加算器
４５、４６、４７、除算器４８を省いたものに相当す
る。
【００３９】図７の係数更新手段は、 Ｗi (t＋1 ）＝Ｗｉ（ｔ）＋βｅ（ｔ）ｘ′（ｔ−ｉ） を求めるように構成されている。この式における各記号
は特許請求の範囲で定義されている。上記を求めるため
に係数演算手段にはｘ（ｔ−ｉ）として端子３１〜３４
にｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ（ｔ−
３）が入力し、端子３５にｅ（ｔ）が入力している。
は特許請求の範囲で定義されている。上記を求めるため
に係数演算手段にはｘ（ｔ−ｉ）として端子３１〜３４
にｘ（ｔ）、ｘ（ｔ−１）、ｘ（ｔ−２）、ｘ（ｔ−
３）が入力し、端子３５にｅ（ｔ）が入力している。
【００４０】第２の実施例の係数更新は、ＬＭＳ（最小
２乗平均アルゴリズム）に基づくものであり、第１の実
施例よりも演算が更に簡略化されている。従って、第１
の実施例と同様に低価格であるにも係わらず性能の良い
ノイズキャンセラ−等を提供することができる。
２乗平均アルゴリズム）に基づくものであり、第１の実
施例よりも演算が更に簡略化されている。従って、第１
の実施例と同様に低価格であるにも係わらず性能の良い
ノイズキャンセラ−等を提供することができる。
【００４１】
【変形例】本発明は上述の実施例に限定されるものでな
く、例えば次の変形が可能なものである。 （１） 適応フィルタを図１に示す個別回路の組み合せ
で構成することができる。 （２） 未知の動的システム２を入力信号１に接続しな
い構成にすることもできる。
く、例えば次の変形が可能なものである。 （１） 適応フィルタを図１に示す個別回路の組み合せ
で構成することができる。 （２） 未知の動的システム２を入力信号１に接続しな
い構成にすることもできる。
【図１】本発明の第１の実施例に従う適応フィルタ装置
を機能的に示すブロック図である。
を機能的に示すブロック図である。
【図２】図１の係数更新手段を機能的に示すブロック図
である。
である。
【図３】係数更新動作を示す流れ図である。
【図４】図３に続く係数更新動作を示す流れ図である。
【図５】白色雑音を入力させた時の収束性を示す図であ
る。
る。
【図６】有色雑音を入力させた時の収束性を示す図であ
る。
る。
【図７】第２の実施例の係数更新手段を機能的に示すブ
ロック図である。
ロック図である。
１ 入力端子 ７〜９ 遅延素子 ２２ 誤差演算回路
Claims (2)
- 【請求項１】 所定のサンプリング周期で基準信号ｄ
(t) を供給するための基準信号供給手段と、 前記基準信号に関係を有する入力信号ｘ(t) を所定のサ
ンプリング周期で供給するための入力信号供給手段と、 前記基準信号供給手段と前記入力信号供給手段に接続さ
れ、１サンプリング周期を単位遅延時間として前記入力
信号を順次に遅延することによってサンプリング時点の
異なる０番目からＬ−１番目までのＬ個（但し、Ｌは２
よりも大きい正の整数）の入力信号ｘ（ｔ）、ｘ（ｔ−
１）、…ｘ（ｔ−Ｌ＋１）を同時に得、前記０番目から
Ｌ番目の入力信号に０番目からＬ−１番目の係数ｗ0
(t) 、ｗ1(t) 、…ｗL-1 (t) を夫々乗算し、これによ
って得られたＬ個の乗算出力を加算した出力ｙ(t) を
得、前記出力ｙ(t) と前記基準信号ｄ(t) との誤差信号
ｅ(t)を求め、前記入力信号ｘ(t) と前記誤差信号ｅ(t)
とに基づいて前記０番目からＬ−１番目の係数の更新
値ｗ0 （ｔ＋１）、ｗ1 （ｔ＋１）…ｗL-1 （ｔ＋１）
を決定するためのディジタルフィルタ信号処理回路とを
有する適応フィルタ装置における係数更新方法であっ
て、 前記更新係数ｗ0 （ｔ＋１）、ｗ1 （ｔ＋１) 、…ｗL-
1 （ｔ＋１）を、 ｗi （ｔ＋１）＝ｗi （ｔ）＋β［ｅ（ｔ）／｛Σj=0
L-1 x ′（ｔ−ｊ）ｘ（ｔ−ｊ）］ｘ′（ｔ−ｉ） ［ここで、各記号は次のものを示す。ｉ及びｊは０から
Ｌ−１までの正の整数；ｔはサンプリング時刻を示す記
号；ｗi （ｔ＋１）はＬ個の更新係数ｗ0 （ｔ＋１）、
ｗ1 （ｔ＋１）、…ｗL-1（ｔ＋１）；ｗi (t) はＬ個
の更新前の係数ｗ0 (t) 、ｗ1 (t) 、…ｗL-1 (t) ；β
は任意の定数；ｅ(t) はサンプリング時刻ｔにおける誤
差信号；ｘ（ｔ−ｊ）はサンプリング時点の異なるＬ個
の入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋
１）；ｘ′（ｔ−ｉ）及びｘ′（ｔ−ｊ）はサンプリン
グ時点の異なるＬ個の入力信号ｘ（ｔ）、ｘ（ｔ−
１）、…ｘ（ｔ−Ｌ＋１）に対応した変形入力信号ｘ′
（ｔ）、ｘ′（ｔ−１）、…ｘ′（ｔ−Ｌ＋１）、｛但
し、このｘ′（ｔ−ｉ）及びｘ′（ｔ−ｊ）は、 ｜ｘ（ｔ−ｉ）｜≧（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （ｃは定数）が成立する時にｓｇｎ［ｘ（ｔ−ｉ）］
（ここで、ｓｇｎ［ｘ（ｔ−ｉ）］はｘ（ｔ−ｉ）＞０
の時に＋１、ｘ（ｔ−ｉ）＜０の時は−１、ｘ（ｔ−
ｉ）＝０の時は０を示す。）であり、 ｜ｘ（ｔ−ｉ）｜＜（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にゼロである。］に従って決定することを
特徴とする適応フィルタ装置における係数更新方法。 - 【請求項２】 所定のサンプリング周期で基準信号ｄ
(t) を供給するための基準信号供給手段と、 前記基準信号に関係を有する入力信号ｘ(t) を所定のサ
ンプリング周期で供給するための入力信号供給手段と、 前記基準信号供給手段と前記入力信号供給手段に接続さ
れ、１サンプリング周期を単位遅延時間として前記入力
信号を順次に遅延することによってサンプリング時点の
異なる０番目からＬ−１番目までのＬ個（但し、Ｌは２
よりも大きい正の整数）の入力信号ｘ（ｔ）、ｘ（ｔ−
１）、…ｘ（ｔ−Ｌ＋１）を同時に得、前記０番目から
Ｌ番目の入力信号に０番目からＬ−１番目の係数ｗ0
(t) 、ｗ1(t) 、…ｗL-1 (t) を夫々乗算し、これによ
って得られたＬ個の乗算出力を加算した出力ｙ(t) を
得、前記出力ｙ(t) と前記基準信号ｄ(t) との誤差信号
ｅ(t)を求め、前記入力信号ｘ(t) と前記誤差信号ｅ(t)
とに基づいて前記０番目からＬ−１番目の係数の更新
値ｗ0 （ｔ＋１）、ｗ1 （ｔ＋１）…ｗL-1 （ｔ＋１）
を決定するためのディジタルフィルタ信号処理回路とを
有する適応フィルタ装置における係数更新方法であっ
て、 前記更新係数ｗ0 （ｔ＋１）、ｗ1 （ｔ＋１) 、…ｗL-
1 （ｔ＋１）を、 ｗi （ｔ＋１）＝ｗi （ｔ）＋βｅ（ｔ）x ′（ｔ−
ｉ）ｘ ［ここで、各記号は次のものを示す。ｉは０からＬ−１
までの正の整数；ｔはサンプリング時刻を示す記号；ｗ
i （ｔ＋１）はＬ個の更新係数ｗ0 （ｔ＋１）、ｗ1
（ｔ＋１）、…ｗL-1（ｔ＋１）；ｗi (t) はＬ個の更
新前の係数ｗ0 (t) 、ｗ1 (t) 、…ｗL-1 (t) ；βは任
意の定数；ｅ(t) はサンプリング時刻ｔにおける誤差信
号；ｘ′（ｔ−ｉ）はサンプリング時点の異なるＬ個の
入力信号ｘ（ｔ）、ｘ（ｔ−１）、…ｘ（ｔ−Ｌ＋１）
に対応した変形入力信号ｘ′（ｔ）、ｘ′（ｔ−１）、
…ｘ′（ｔ−Ｌ＋１）、｛但し、このｘ′（ｔ−ｉ）
は、 ｜ｘ′（ｔ−ｉ）｜≧（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ （ｃは定数、ｊは０からＬ−１までの正の整数）が成立
する時にｓｇｎ［ｘ（ｔ−ｉ）］（ここで、ｓｇｎ［ｘ
（ｔ−ｉ）］はｘ（ｔ−ｉ）＞０の時に＋１、ｘ（ｔ−
ｉ）＜０の時は−１、ｘ（ｔ−ｉ）＝０の時は０を示
す。）であり、 ｜ｘ′（ｔ−ｉ）｜＜（ｃ／Ｌ）Σj=0 L-1 ｜ｘ（ｔ−
ｊ）｜ が成立する時にゼロである。］に従って決定することを
特徴とする適応フィルタ装置における係数更新方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05348177A JP3092647B2 (ja) | 1992-12-30 | 1993-12-24 | 適応フィルタ装置 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4-360631 | 1992-12-30 | ||
| JP36063192 | 1992-12-30 | ||
| JP05348177A JP3092647B2 (ja) | 1992-12-30 | 1993-12-24 | 適応フィルタ装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06318846A true JPH06318846A (ja) | 1994-11-15 |
| JP3092647B2 JP3092647B2 (ja) | 2000-09-25 |
Family
ID=26578691
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP05348177A Expired - Fee Related JP3092647B2 (ja) | 1992-12-30 | 1993-12-24 | 適応フィルタ装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3092647B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100937379B1 (ko) * | 2001-04-26 | 2010-01-20 | 톰슨 라이센싱 | 고선명 텔레비전 신호와 같은 수신된 동기 데이터 신호를 처리하기 위한 적응 등화기 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005229312A (ja) | 2004-02-12 | 2005-08-25 | Daihen Corp | 適応型ディジタルフィルタ |
-
1993
- 1993-12-24 JP JP05348177A patent/JP3092647B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100937379B1 (ko) * | 2001-04-26 | 2010-01-20 | 톰슨 라이센싱 | 고선명 텔레비전 신호와 같은 수신된 동기 데이터 신호를 처리하기 위한 적응 등화기 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3092647B2 (ja) | 2000-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3008763B2 (ja) | 適応フィルタによるシステム同定の方法および装置 | |
| JP2885269B2 (ja) | 適応制御ろ波器 | |
| JP3068996B2 (ja) | 適応フィルタ装置及びそのステップゲインの発生方法 | |
| JP2947093B2 (ja) | 適応フィルタによるシステム同定の方法および装置 | |
| JPH04259111A (ja) | 適応等化器 | |
| EP0878920A2 (en) | Echo canceler employing segment-based reinitialization | |
| US6223194B1 (en) | Adaptive filter, step size control method thereof, and record medium therefor | |
| JPH06318846A (ja) | 適応フィルタ装置における係数更新方法 | |
| US5260896A (en) | Adaptive filter and a method of preventing divergent behavior of the adaptive filter | |
| JP3707443B2 (ja) | 適応忘却係数制御適応フィルタ、および忘却係数適応制御方法 | |
| JP2541044B2 (ja) | 適応フィルタ装置 | |
| Gannot et al. | On the application of the unscented Kalman filter to speech processing | |
| JP3180739B2 (ja) | 適応フィルタによる未知システム同定の方法及び装置 | |
| JPH03266516A (ja) | アダプティブ・フィルタ適応化方法及び装置 | |
| Aboulnasr et al. | Efficient implementation of echo cancellers for ISDN applications | |
| JP3236226B2 (ja) | エコーキャンセラ | |
| JP3708040B2 (ja) | エコーキャンセラ | |
| JP3147864B2 (ja) | 適応ステップサイズ制御適応フィルタ、及び適応ステップサイズ制御方法 | |
| JPH1168518A (ja) | 係数更新回路 | |
| JP3121998B2 (ja) | 音響反響除去装置 | |
| JPS61273008A (ja) | 適応形トランスバ−サルフイルタ | |
| JP3217619B2 (ja) | 音響反響除去装置 | |
| JPH05252000A (ja) | 信号処理装置 | |
| JP2541040B2 (ja) | 適応フィルタ装置における係数更新方法 | |
| JP2001217685A (ja) | 適応フィルタ、その制御方法及びプログラムを記憶した記憶媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |