JPH02214322A - 最小2乗平均アルゴリズムを実施するためのパイプライン式プロセッサ - Google Patents
最小2乗平均アルゴリズムを実施するためのパイプライン式プロセッサInfo
- Publication number
- JPH02214322A JPH02214322A JP1321265A JP32126589A JPH02214322A JP H02214322 A JPH02214322 A JP H02214322A JP 1321265 A JP1321265 A JP 1321265A JP 32126589 A JP32126589 A JP 32126589A JP H02214322 A JPH02214322 A JP H02214322A
- Authority
- JP
- Japan
- Prior art keywords
- multiplier
- convolution
- input
- output
- data
- 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
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 230000003044 adaptive effect Effects 0.000 claims abstract description 24
- 238000009825 accumulation Methods 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 7
- 230000035508 accumulation Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 10
- 230000004044 response Effects 0.000 description 7
- 230000003111 delayed effect Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 101100493710 Caenorhabditis elegans bath-40 gene Proteins 0.000 description 1
- 235000010724 Wisteria floribunda Nutrition 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Complex Calculations (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
め要約のデータは記録されません。
Description
定的には適応フィルタ型データプロセッサに関する。
いは時間変化、適応型システムを表わす。
ス関数が励起およびインパルス応答結果のために用いら
れる。何らかの池の入力波形が与えられた場合、インパ
ルス応答を入力波形とコンボルブすることにより該イン
パルス応答を出力を決定するために利用することができ
る。これに対して、適応型システムにおいては、インパ
ルス関数は各データサンプル点においであるいはデータ
サンプルの各組において変化している。従って、適応型
システムは時間変化しておりあるいはノンリニアである
。デジタルフィルタがフィルタの所望の関数に応じて適
応プロセスにより特定の伝達関数に向けられている。デ
ジタルフィルタの1つの応用は所定の伝達関数により1
組のデータをろ波する場合における最小2乗平均誤差を
得ることである。この関数は最小2乗平均(L M S
: Ieast−mean−square )アルゴ
リズムとして知られている。
際に出力しているものとの間の誤差を表わしかつフィル
タの誤調整の測定値となる。LMSアルゴリズムは、(
1)部分的に先のエラーに基づくフィルタの係数の更新
である適応計算、(2)所定の係数値およびデータ値の
積の累積であるコンボリューション計算、(3)そして
エラー計算を含む、エラー計算はコンボリューション計
算の出力値を所定の所望値から減算することである。
社のMB87XXAFPが例にあげられ、これはLMS
適応FIRフィルタアルゴリズムをサポートすることが
可能である。知られたLMSアルゴリズムシステムにお
いては、コンボリューション計算が最初に行なわれそこ
からエラーが決定される0次に前記エラーが利用される
、更新計算がなされる。非常に高度な多重タップのフィ
ルタに対しては、このことはすべてのタップに対するコ
ンボリューション計算が最初に計算されなければならず
、そして次にすべてのタップに対するエラーおよび更新
計算がなされなければならない。
は多重タップフィルタを実施するために非常に多くのプ
ロセッサのクロックサイクルが必要となり、係数、デー
タおよびエラー値を格納するために多重メモリ構造を使
用する必要があり、そしてパイプライン式多重かつ分離
されたフィルタにおける効率の悪さを含む、典型的には
、パイプライン様式で多重フィルタ機能を実施する従来
のシステムはカスケード構造で動作する個々の集積回路
製品を使用した。
ゴリズムを実施するための改良されたパイプライン式プ
ロセッサを提供することにある。
改良されたデジタルフィルタ構造を提供することにある
。
ルゴリズムを実施するための改良された方法を提供する
ことにある。
るこれらおよび他の目的を実施するにあたり、1つの形
式として、Nを整数とした場合に、Nタップのフィルタ
における最小2乗平均(LMS)アルゴリズムを実施す
るための方法およびパイプライン式プロセッサが提供さ
れる。
のタップに対して計算される。実質的に同時に、フィル
タの各タップに対する係数値を更新することによりプロ
セッサは一連の係数値の適応計算を連成する。パイプラ
イン式プロセッサの更新部は所定の時間期間の開所定数
の更新計算を行なう、コンボリューション部が更新部に
結合され同じ時間期間の開所定数のコンボリューション
積を連成する。記憶部が更新部およびコンボリューショ
ン部に結合され最小平均2乗アルゴリズムを実施するた
めに更新部およびコンボリューション部の双方によって
要求されるデータ値および係数値を選択的に提供する。
組合わせて以下の詳細な説明からより明瞭に理解できる
であろう。
程式によって構成される。2つの方程式の内の1つは典
型的には更新(update)方程式と称されかつ他方
の方程式はコンボリューション方程式として言及される
。更新方程式はデジタルフィルタの適応部分であり伝統
的なニュートンーラルフソン反復アルゴリズムと同様の
様式でフィルタ係数を修正する。コンボリューション方
程式は伝統的なFIRフィルタ方程式である。双方の方
程式は各データサンプル期間°にフィルタの各タップに
おいて連成される。
ボリューションを実施するための時間不変フィルタ構造
10が示されている。このコンボリューションは次の方
程式で表わされる。
0 に=o、 1. 2.−00
(1)上述の式において、jはタップ番号値を表わし、
kは時間そしてNはフィルタにおけるタップ数を表わす
、フィルタ構造10の入力端子はタップ0および時間k
におけるデータ値Xを表わすX(0゜k)と称されるデ
ータ入力を受ける。データ入力値は乗算器11の第1の
端子に接続されかつZ領域におけるz−1で表わされる
遅延機能を有する遅延回路12の入力に接続されている
。W(0)で表わされる係数値が乗算器11の第2の入
力に接続されている。遅延回路12の出力はタップ1お
よび時間kにおけるデータ値Xを表わすデータ人力X(
1,k)を与える。遅延図8?112の出力は乗算器1
3の第1の入力および遅延回路14の入力に接続されて
いる。W(1)で表わされる係数値が乗算器13の第2
の入力に接続されている。
ータfiiXを表わすデータ入力X(2,k>を与える
。遅延図f#114の出力は乗算器15の第1の入力に
接続されている0乗算器15の第2の入力はW(2)で
表わされる係数値を受ける。乗算器15の出力は加算器
回路16の第1の入力に接続されている。加算器回路1
6の第2の入力は2進ゼロ値を受ける。加算器回路16
の出力は加算器回路17の第1の入力に接続されている
9乗算器13の出力は加算器回路17の第2の入力に接
続されている。加算器回路17の出力は加算器回路18
の第1の入力に接続されている0乗算器11の出力は加
算器回路18の第2の入力に接続されている。加算器回
路18の出力は方程式(1)で規定されるコンボリュー
ション出力信号Y(k)を与える。
定の係数値およびタップデータ値の積を表わす個々の積
を受ける。加算器回路16.17および18はこれらの
積を加算して方程式lで規定される時間kにおける累積
された出力信号Y(k)を提供するl初は、タップ2に
おけるデータ入力はW(2)で乗算されかつゼロと加算
される0次に同じ時間インターバルkにおいて、最初の
積がW(1)により乗算されたタップ1におけるデータ
入力の積と加算される。加算器回路17の加算結果はW
(0)により乗算されたタップゼロにおける出力データ
から得られる積と加算されY (k)を与える。係数値
W(0)、W(1)およびW(2)は固定された、変化
しない係数でありこれがインパルス応答が時間的に不変
であるという理由である。
この構造は方程式1で表わ、されるように、適応計算お
よびコンボリューション演算の双方が行なわれる適応型
または時間変化型フィルタの単一フィルタタップを表わ
す、言替えれば、第1図の係数値W(0)、W(1)お
よびW(2)が各時間期間に形成され時間的に変化する
。このL MSアルゴリズムに対する適応計算は次のよ
うに表わされる。
1)11X(j、に−1) ・・−(2)上式にお
いて、j=0.・・・、N−1であり、かつに=o、1
,2.・・・、■。
入力を有する。遅延回路21の出力はデ−少入力X(2
,k)を提供しかつ乗算器22の第1の入力および乗算
器25の第1の入力の双方に接続されている0乗算器2
2の第2の入力はエラー信号βe (k)に接続されて
いる。βの項はフィルタの用途に応じて選択される所定
の定数値である0乗算器22の出力は加算器23の第1
の入力に接続されている。加算器23の出力は遅延回路
24の入力に接続されている。遅延回路24の出力は加
算器回路23の第2の入力および乗算器25の第2の入
力の双方に接続されている。遅延回路24は更新された
係数値W(2,k)を提供する0乗算器25の出力は加
算揺回111g26の第1の入力に接続されている。加
算器回路26の第2の入力は2進上口入力を受ける。加
算揺回126の出力は式1のコンボリューションアルゴ
リズムにおける信号Y(k)の第1の累積を与える。
1から遅延回路21によって受信される。
乗算器22および乗算器25の双方に結合される0乗算
器25は同時に係数値W(2,k)を受けかつ方程式1
に従って乗算器25の出力に初期コンボリューション積
を提供する。同時に、係数値W(2,k)が[X (2
,k)*βe (K)]の積と加算され乗算回路25に
よる引き続く使用のために更新された係数値W(2,に
+1)が提供される。このようにして、タップ位置2に
対する係数値が遅延回路21およびエラー環βe(k>
により提供されるタップ値の関数として常に変更される
。フィルタ構造20のような構造が各タップに対し必要
とされるから、マルチタップのフィルタに対し要求され
る回路の量は非常に大きい。
30が示されている。構造30においては、データ人力
X(1,k)が遅延回路31の入力に接続されている。
これは遅延回路33の入力および乗算器37の第1の入
力の双方に結合されている。遅延回路33の出力は遅延
されたデータ人力X(2,に−1>を提供しかつ乗算器
34の第1の入力に接続されている1乗算器34の第2
の入力はエラー環βe(k−1)を受ける。
されている。加算器回路35の出力は係数値W(2,k
)を提供しかつ遅延回路36の入力および乗算器37の
第2の入力の双方に接続されている。遅延回路36の出
力は係数値W(2,に1)を提供しかつ加算器回路35
の第2の入力に接続されている。乗算器37の出力は加
算器回路38の第1の入力に接続されている。加算器回
路38の第2の入力は論理ゼロの値に接続されている。
ョンアルゴリズムを実施するうえで信号Y(k)の第1
の累積を提供する。
24の位置を再配置することによりフィルタ構造20と
かなり構造的に異なっているが機能的には等価である。
けることによって、同じ機能を連成することができる。
って実施された遅延機能は遅延回路33および36そし
てエラー環βe(k−1)を乗算器34に結合すること
により連成される。フィルタ構造20のように、フィル
タ構造30は方程式1の適応規則および方程式2のコン
ボリューションアルゴリズムの双方を行なうことによっ
てLMSアルゴリズムを実施する。
ルゴリズムを実施するための3タツプフイルタであるフ
ィルタ構造40が示されている。
も利用できかつ第4図における3タツ1の実施例は限定
ではなく説明のためだけにすぎない。
算器44の第1の入力に結合されている。
)を提供しかつ乗算器46の第1の入力に接続されてい
る0乗算器46の第2の入力は工ラー項βe(k−1>
に接続されている0乗算器46の出力は加算器回路48
の第1の入力に接続されている。加算器回路48の出力
は更新された係数値W(0,k)を提供しかつ乗算器4
4の第2の入力および遅延回路50の入力の双方に接続
されている。遅延回路50の出力は加算器回路48の第
2の入力に接続されている。遅延回路42の出力はまた
乗算器52の第1の入力および遅延回路54の入力の双
方に#続されている。遅延回路54の出力は遅延された
データ人力X(2,k>を提供しかつ乗算器56の第1
の入力に接続されている0乗算器56の第2の入力はエ
ラー環βe(k−1)に接続されている0乗算器56の
出力は加算器回路58の第1の入力に接続されている。
k)を提供しかつ乗算器52の第2の入力および遅延回
路60の入力の双方に接続されている。遅延回路60の
出力は加算器回路58の第2の入力に接続されている。
び遅延回路62の入力の双方に接続されている。遅延回
路62の出力は遅延されたデータ人力X(3,k)を提
供しかつ乗算器66の第1の入力に接続されている0乗
算器66の第2の入力はエラー環βe(k−1)に接続
されている0乗算器66の出力は加算器回路68の第1
の入力に接続されている。加算器回路68の出力は更新
された係数値W(2,K)を提供しかつ乗算器64の第
2の入力および遅延回路70の入力の双方に接続されて
いる。遅延回路70の出力は加算器回路68の第2の入
力に接続されている1乗算器64の出力は加算器回路7
2の第1の入力に接続されている。加算器回#172の
第2の入力は2進ゼロの値に接続されている。加算器回
#I72の出力は加算器回路74の第1の入力に接続さ
れている0乗算器52の出力は加算器回路74の第2の
入力に接続されている。加算器回路74の出力は加算器
回路76の第1の入力に接続されている0乗算器44の
出力は加算器回路76の第2の入力に接続されている。
ューション出力信号Y(k)を提供する。
ける各タップ計算は更新計算(方程式2)およびコンボ
リューション計算(方程式1)を必要とする。これらの
2つの計算を連成するために各々の新しいデータサング
ルに対しフィルタタッグごとに2つの乗算/加算操作が
必要とされる。
新およびコンポリ−シュン計算のパイプライン化をも可
能とする。各加算器および各乗算器はそのそれぞれの計
算を形成するため完全な半サイクルを持つ事を許容され
るであろう、タッフ′データ人力Xおよび係数値Wが示
すように、各タップ位置におけるすべての操作はサイク
ルにのような単一時間サイクルの間に実質的に同時に起
っている0例えば、タップ位置3においては、乗算器6
4および66が時間サイクルにの最初の半分の間に機能
している9時間サイクルにの第2の半分の間に、加算器
68および72が機能している。
位置1の乗算器44および46もまた時間サイクルにの
最初の半分の間に機能している。タップ位置2の加算器
58および74とタップ位置1の加算器48および76
もまた時間サイクルにの第2の半分の間に機能している
。
1のコンボリューションアルゴリズムの稜部分を実施す
るためW(2,k)およびX(2゜k)を乗算している
。同時に、乗算器66はX(3,k>およびβe(k−
1)の積を形成しておりこれは方程式2の適応アルゴリ
ズムの部分である。データ人力X(3,k>はX(2,
に−1>と等価であることに注意すべきであるが、ここ
で(k−1)は時間サイクルにの前の時間サイクルであ
る。従って、時間サイクルkに対する更新された係数値
Wは先の時間サイクル(k−1>からのデータおよびエ
ラー情報でもって時間サイクルkに形成される。
成されたW(2,k)およびX(2,k)の積を強制さ
れた2進ゼロと加算している。引き続きかつ順次に加算
器74および76は引き続くタッグに対するトータルの
累積を完了し方程式1のコンボリューションアルゴリズ
ムの実行を完了する0時間サイクルにの第2の半分の間
に同時に加算器68は遅延回路70によって提供された
先の係数値W(2,k)を乗算器66により提供された
積と加算し次の時間サイクル(k+1>に対する更新さ
れた係数値W(2,に+1)を提供している。各タップ
位置に対する計算の各々は実質的に同時に起っているが
、計算は順次パイプライン化されかつ右側の最も高いタ
ッグ位置から左側の最も低いタップ位置へと計算できる
ことに注目すべきである。
ップのフィルタ構造が第3図の個々のタツ7梢造と異な
ることに注目すべきである。特に、第3図のフィルタ構
造30の遅延回路33は第4図のフィルタ構造40の各
フィルタのタッグ位置において削除されている0例えば
、フィルタ構造40のフィルタタップ位置0においては
、この遅延機能は遅延回路42によって提供される。フ
ィルタタップ位置1においては、該遅延機能は遅延回路
54により提供され、かつフィルタタップ位置2におい
ては遅延機能は遅延回路62によって提供される。多重
タップフィルタを実施するのに必要とされる遅延回路の
数を最小化することに加え、フィルタ構造40はまた遅
延回路42の出力から双方の乗算器46および52を同
時に駆動または再符号化(recode) Lかつ遅延
回路54の出力から双方の乗算器56および64を駆動
または再符号化するために単一のデータ値を用いること
ができるという利点を提供する。以下に示されるように
、2つの乗算器の同時的な再符号化能力は速度、サイズ
およびパイプライン化の利点を提供する。
ルタ構造40を実施するパイプライン式グロセッサ79
が示されている。特に、プロセッサ79はLMSアルゴ
リズムを実行するよう機能する。プロセッサ79は双方
向データバス80により通信する。一般に、プロセッサ
79は時間変化する係数値を与えるための更新部81、
コンボリューション部82、そして制御部83を有する
ものと考えることができる。ラッチ85の入力はデータ
バス80に接続されている。ラッチ85の出力はマルチ
プレクサ86の第1の入力に接続されている。制御信号
C4はマルチプレクサ86の制御入力に接続されている
。マルチプレクサ86の第2の入力は2進ゼロの値に接
続されている。
#枕されている0乗算器87の出力はラッチ88の第1
の入力に接続されている。ラッチ89の入力はデータバ
ス80に接続されている。ラッチ89の出力はマルチプ
レクサ90の第1の入力に接続されている。rcl、と
名付けられた制御信号はマルチプレクサ90の制御入力
に結合されている。マルチプレクサ90の出力はレコー
ダ(再符号化器)回路92および93の各々の入力に接
続されておりかつラッチ94の入力に接続されている。
力に接続されている。レコーダ回路92の第2の出力は
ラッチ88の第2の入力に接続されている。ラッチ88
の第1の出力は「C」と名付けられた信号を提供しかつ
加算器96の第1の入力に接続され、そしてラッチ88
の第2の出力は加算器96の制御入力に接続されている
。加算器96の出力はラッチ97の入力に接続されてい
る。ラッチ97の出力はラッチ98の入力に接続されて
いる。ラッチ94および98の各々の出力はランダムア
クセスメモリ(RAM)99の入力に接続されている。
れの制御入力に結合されている。RAM99の出力はラ
ッチ101およびラッチ102の双方の入力に接続され
ている。ラッチ102の出力はマルチプレクサ90の第
2の入力に接続されており、かつラッチ101の出力は
ラッチ104の入力に接続されている。ラッチ104の
出力は加算器96の第2の入力に接続されている。
第1の入力がr W 、と名付けられた信号を受信する
ためにラッチ97の出力に接続されている。レコーダ9
3の第1の出力は乗算器106の第2の入力に接続され
ている0乗算器106の出力はラッチ108の第1の入
力に接続されている。レコーダ93の第2の出力はラッ
チ108の第2の入力に接続されている。ラッチ108
の・出力は「T」と名付けられた信号を提供しかつ加算
器110の第1の入力に接続されている。加算器110
の加算または減算動作モードのいずれかを選択するため
のラッチ108の制御出力は加算器110に接続されて
いる。加算器110の出力はラッチ112の入力に接続
されている。ラッチ112の出力はラッチ114および
116の双方の入力に接続されている。ラッチ114の
出力は加算器110の第2の入力に接続されている。リ
セット信号がラッチ114のリセット端子に接続されて
おり、かつ制御信号C5がラッチ116の制御入力に結
合されている。ラッチ116の出力はバス駆動回路11
8の入力に接続されている。バス駆動回路118の出力
はデータバス80に接続されている。
ューション操作を連続して連成するためにパイプライン
様式で機能する。プロセッサ7つはまたLMSアルゴリ
ズムを実施するために必要とされる乗算および加算操作
をパイプライン処理する。このようにして計算を組織化
することにより各加算器および乗算器がその結果を形成
するために完全な半サイクルを持つことを許容する。ま
た、単一のRAM構造がタップ重み付けおよびフィルタ
データの双方を保持するために利用されている。係数お
よびデータに対する独立のRAMの使用を除去すること
によりかなりの量の回路領域を節約する。
に理解できる。第6図には1つのコンボリューション出
力Y(k)を提供するために発生するマシン操作を図示
するタイミングシーケンス図か示されている0図示され
たステップは次に後続の時間変化する出力信号Y (k
)を提供するために繰返される1図示された形態におい
ては、第4図に示されるような3タツプフイルタが示さ
れそこでは出力を提供するために4マシンクロツクサイ
クルが必要とされる。−数的に、本発明はNタップのフ
ィルりを(N+1)クロックサイクルで実施するが、こ
こでNは整数である。各クロックサイクルは2つの位相
、即ちφ1およびφ2を有し、これらはほぼ50%のデ
イ−ティサイクルを有するオーバラッグしないクロック
である。クロック位相φ1の間は、更新部81の加算器
96およびコンボリューション部82の加算器110が
動作している。クロック位相φ2の間は、更新部81の
乗算器87およびコンボリューション部82の乗算器1
06が動作している。第1のクロックサイクルの第1の
クロック位相においては加算器96はいまだ必要とされ
ずかつ動作していない、RAM99はデータバス80か
らラッチ89およびマルチプレクサ90を介して先に受
信された入力データX(3,k)を提供するために読ま
れる。加算器110はラッチ114を通りフィードバッ
ク経路を介して与えられるY(k)の現在の値を時間期
間(k−1)の間に乗算器106によって与えられる最
後の乗算器出力値T(0,に−1)と加算することによ
り先の時間期間(k−1)から更新されたY(k>の値
を提供するよう機能する。
る。同時に、係数値W(2,に−1>が加算器96の第
1の入力においてクロックサイクル2での引き続く使用
のためRAM99から読まれる。
第3のタップ位置において第4図のフィルタ構造40か
ら確認することができる。
は、乗算器87は出力信号C(2,k>を提供するよう
機能しこの出力信号は引き続き加算器96によって使用
される。デー°少入力X(3゜k)がRAM99からマ
ルチプレクサ90を介してレコーダ92に結合される。
よびマルチプレクサ86を介して提供される。!&初は
、マルチプレクサ86は制御信号C4によって制御され
エラー更新がないことを表わすゼロ値入力を提供する。
、図示しない回路か乗算器87による使用のためエラー
項βe(k−1)を提供する。レコーダ92は乗算器8
7による使用のための所定の再符号化アルゴリズムを実
施するよう11能する。
精度のデータに備えることができることを理解すべきで
ある。第1のクロックサイクルにおいては、乗算器10
6は出力を提供しないが、これは加算器96がいまだW
オペランドに対する値を提供していないからである。第
1のクロックサイクルにおいてRAM99のアドレスA
がアドレスされ、かつクロックサイクルφlの間にRA
M99が読み取られることに注目すべきである。また第
1のクロックサイクルの第2のクロック位相の間ラッチ
114はゼロにリセットされる。
は、RAM99は乗算器87による後続の使用のために
データ人力X(2,k)を提供するために読み取られる
。係数値W(2,に−1)およびオペランドC(2,k
)はそれぞれRAM 99および乗算器87から各々加
算器96の第1および第2の入力に接続される。また、
第2のクロックサイクルの第1の位相の間は、加算器1
10は動作していないが、これは乗算器106が加算器
110にTオへランドを提供するために第1のクロック
サイクルにおいて動作していなかったからである。係数
値W(1,に−1)はまた加算器96による後続の使用
のために第1のクロック位相の間にRAM99から読ま
れる。
は、RAM99は第1のクロック位相で計算された係数
値W(2,k)および新しいデータ入力X(3,に+1
)の双方を格納するために書き込まれる。レコーダ92
にはデータ入力X(2,k>が与えられ、即ちレコーダ
92にはデータ入力X(2,k>が与えられかつ乗算器
106には係数値W(2,k)が与えられる0乗算器8
7はC(1,k>を計算し、かつ乗算器106はT(2
,k)を計算する。第2のクロックサイクルの間、R,
AM99のアドレスB位置が制御信号C2に応答して選
択される。
は、RAM99が乗算器87による後の使用のためにX
(1,k)を提供しかつW(Ok−1)を提供するため
に読み取られる。オペランドW(1,に−1>およびC
(1゜k)がそれぞれRAM99および乗算器87から
加算器96の第1および第2の入力にそれぞれ接続され
る。加算器96はW(1,k)を計算しこれはラッチ9
7により格納されかつ後続のクロック位相においてRA
M99に書き込まれると共に乗算器99に結合される。
算器106によって先のクロック位相において提供され
たオペランドT(2,k)をゼロオペランド値と加算し
オペランドY(k)を提供する。第3のクロックサイク
ルの間、RAM99のアドレスC位置が制御信号C2に
応じて選択される。 第3のクロックサイクルの第2の
クロック位相においては、RAM99が第3のクロック
サイクルの第1のクロック位相において計算された係数
オペランドW(1,k)および新しいデータ値X(2,
k>の双方を格納するために書き込まれる。データ人力
X(1,k)がレコーダ92および93の双方に結合さ
れ、かつ係数オペランドW(1,k)が乗算器106に
接続される0乗算器87はC(0,k)を計算し、がっ
乗算器106はT(2,k)を計算する。
は、RAM99は読み取られない、先に読まれた係数オ
ペランドW(0,に−1)および先に計算されたC(0
,k)はそれぞれ加算器96の第1および第2の入力に
接続される。また、第4のクロックサイクルの第1の位
相の間、加算器110は先のY (k)の値を第3のク
ロックサイクルの第2のクロック位相において乗算器1
06により計算されたT(1,k)と加算することによ
り更新されたY(k)の値を計算する。i&後の値Y(
k>がラッチ116に格納されかつ制御信号C5に応答
してデータバス80に出力される。
かつこの値を所定の正しい値から減算して後続の(k+
1)時間フレームにおいて使用されるべき新しいエラー
環βe (k)を計算する。
は、RAM99が第1のクロック位相で計算された係数
値W(0,k)を格納するために書き込まれる。レコー
ダ93にはデータ入力×(0,k)が与えられ、かつ乗
算器106には係数IW(0,k>が与えられる。更新
部81の乗算器87は第4のクロックサイクルの間は機
能しないが、それは最低のタップ位置、ゼロ、に対する
Cオペランドは先に計算されたからである。第4のクロ
ックサイクルの間、RAM99のアドレスA位置は読み
出しおよび書き込みのため制御信号C2に応答して再び
選択され得る。
の動作は繰返される。即ち、Y(k十1)のようなYオ
ペランドに対する後続の値が計算されかつLMSアルゴ
リズムが実施される。また、多重フィルタが単一の時間
フレームにおいてパイプライン様式で実施できこの場合
冬時間フレームは複数のクロックビートを含み、各複数
のものは図示されたクロックの1つに等価である0以上
の説明により第4図の3タツプフイルタ構造がプロセッ
サ79により4クロツクサイクルで実施されることが明
らかであろう、プロセッサ79はそれぞれ各形式のオペ
ランドを格納するために2つのRAMを使用することに
反して、フィルタ係数オペランド、W、およびデータオ
ペランド、X、の双方を格納するために単一のRAMが
利用できるように実施される。また、乗算器56および
乗算器64の双方が同時に同じデータオペランドに結合
される第4図のフィルタ構造はレコーダ92および93
が共にデータバス80またはRAM99から同じデータ
オペランドによって同時に駆動できるようにプロセッサ
79が実施できることを意味する。
フィルタのパイプライン式動作を含む。
セッサまたは集積回路でNのコンボリューション操作を
完全に連成しかつ次にコンボリューションおよび更新の
双方の操作を行なうな、めの他のプロセッサまたは単一
ハードウェア回路の時分割多重によりNの更新操作を行
なうことによりLMSアルゴリズムを実施している。こ
の場合、Nの更新操作に関連するオペランドの各々は更
新操作のそれぞれのオペランドに選択的に結合されなけ
ればならなかった。これに対し、本発明においては、更
新構造は本質的にコンボリューション構造において与え
られる。言い換えれば、本発明におけるNタップフィル
タはコンボリューションおよび更新の双方をN回計算す
ることにより連成して実施される。これはプロセッサの
実行時間をかなり短縮しかつNタップフィルタを実施す
るために使用される回路の量を最小化する。多重フィル
タがこのようにして独立のフィルタ間の遷移の間クロッ
クサイクルの損失なしに実施できる6本発明のプロセッ
サ構造は制御信号C4を介しマルチプレクサ86を制御
可能とすることにより係数更新操作の外部制御を提供す
る。いずれの新しいデータサンプル対しても、外部制御
回路はフィルタの係数が更新されているか否かを決定す
ることかできる。伝統的な外・部制御回路もまたマルチ
プレクサ86において適応ベータ(β)を選択すること
により更新レートを決定することができる。外部回路は
またもしそれが望ましい場合にはタップデータを知られ
た状態にすることができる。
者には本発明は多くの方法で変形することができかつ上
に説明されたもの以外の多くの実施例を含むことができ
ることは明らかである。従って、添付の請求の範囲は本
発明の真の精神および範囲内にある発明のすべての変形
をも含むことを意図している。
ョン方程式を実施するための知られたフィルタ構造を示
すブロック図、 第2図は、LMSアルゴリズムに対する適応法則および
コンボリューションを実施するための他の知られたフィ
ルタ構造を示すブロック図、第3図は、LMSアルゴリ
ズムに対する適応ルールおよびコンボリューションを実
施するための他の知られたフィルタ構造を示すブロック
図、第4図は、本発明に係わるLMSアルゴリズムの適
応ルールを実施するためのフィルタ構造を示すブロック
図、 第5図は、本発明に係わる第4図のフィルタ構造を実施
するためのパイプライン式プロセッサを示すブロック図
、そして 第6図は、本発明のプロセッサのパイプライン式動作を
説明するための図式的説明図である。 40:フィルタ構造、 42、 50. 54,60,62.70:遅延回路、 44.46,52,56,64,66:乗算器、48.
58.68,72,74,76二加算器、79:パイプ
ライン式プロセッサ、 80:データパス、 8に更新部、 82:コンボリューション部、 83:制御部、85.
88,89,94,97,98,101104.108
,112,114 116:ラッチ、 86.90:マルチブレフサ、 87.106:乗算器、 92.93:レコーダ、 96.110:加算器、 99 : RAM、118:
バスドライバ回路。
Claims (1)
- 【特許請求の範囲】 1、Nを整数とし、Nタップに対する一連のデータ値お
よび係数値のコンボリューションを計算しかつほぼ同時
にフィルタの各タップに対する係数値を更新することに
より該一連の係数値の適応計算を行なうことにより、N
タップフィルタにおける最小2乗平均(LMS)アルゴ
リズムを実施するためのパイプライン式プロセッサであ
って、該プロセッサは、 LMSアルゴリズムを実施するのに必要とされる係数値
に対し所定数の更新計算を行なうための更新手段であっ
て、前記更新計算は所定の期間の間に行なわれるもの、 前記更新手段に結合され同じ期間の間に各々所定の係数
値と所定のデータ値を乗算することにより形成される所
定数のコンボリューション積を加算し前記コンボリュー
ションを表わす出力信号を提供するためのコンボリュー
ション手段、そして前記更新手段およびコンボリューシ
ョン手段の双方に結合され前記更新手段およびコンボリ
ューション手段の双方によって必要とされるデータ値お
よび係数値を選択的に提供し最小2乗平均アルゴリズム
を実施するための格納手段、 を具備することを特徴とするパイプライン式プロセッサ
。 2、前記所定の期間は所定の周波数のクロック信号の(
N+1)クロックサイクルを有する請求項1に記載のパ
イプライン式プロセッサ。 3、前記更新手段およびコンボリューション手段は各々
再符号化された乗算器を具備しかつ前記格納手段は選択
的に同時に前記更新手段およびコンボリューション手段
の各々の再符号化された乗算器に同じデータ値を提供す
る請求項1に記載のパイプライン式プロセッサ。 4、データ値および係数値の双方を格納する前記格納手
段は単一のメモリ構造である請求項1に記載のパイプラ
イン式プロセッサ。 5、前記更新手段はさらに、 被乗数オペランドを受けるための第1の入力、再符号化
された乗算器オペランドを受けるための第2の入力、そ
して出力を有する乗算器であって、該乗算器は前記フィ
ルタのために選択された所定のエラーオペランドおよび
再符号化されたデータオペランドを乗算して更新された
係数値を提供するもの、そして 前記乗算器の出力に結合された第1の入力、そして前記
格納手段に結合された第2の入力を有する加算器であっ
て、該加算器は前記フィルタの各タップに対し係数更新
値および前記格納手段により与えられた所定の先の係数
の和を提供し、それにより適応計算を行なうもの、 を具備する請求項1に記載のパイプライン式プロセッサ
。 6、前記コンボリューション手段はさらに、被乗数オペ
ランドを受けるための第1の入力、再符号化された乗算
器オペランドを受けるための第2の入力、そして出力を
有する乗算器であって、前記乗算器は前記フィルタの各
タップに対する所定の更新された係数値および再符号化
されたデータオペランドを乗算しコンボリューション項
を提供するもの、そして 前記乗算器の出力に結合された第1の入力、そしてその
出力に結合された第2の入力を有する加算器であって、
該加算器は前記フィルタのNのタップに対し一連のデー
タ値および係数値のコンボリューションを実施するため
にNの累算を行なうもの、 を具備する請求項1に記載のパイプライン式プロセッサ
。 7、Nを整数とすると、Nのタップに対する一連のデー
タ値および係数値のコンボリューションを計算しかつほ
ぼ同時にフィルタの各タップに対する係数値を更新する
ことにより前記一連の係数値の適応計算を行なうことに
よりNタップのフィルタ機能を連成するパイプライン式
プロセッサにおける最小2乗平均(LMS)アルゴリズ
ムを実施するための方法であって、該方法は、 所定の時間フレームの間に所定数の更新計算を行なうた
めにプロセッサの更新回路部を時間多重する段階、 前記と同じ期間の間に所定数のコンボリューション積を
行ない、最終コンボリューション積はLMSアルゴリズ
ムの完了した計算を表わすコンボリューション回路部を
時間多重する段階、そして単一の格納手段を更新回路部
およびコンボリューション回路部の双方に結合して最小
2乗平均を実施するために更新回路部およびコンボリュ
ーション回路部の双方によって必要とされるデータ値お
よび係数値を提供する段階、 を具備することを特徴とする前記方法。 8、さらに、 前記所定の時間フレームを規定するために所定の周波数
のクロック信号を利用する段階であって、前記最小2乗
平均アルゴリズムは時間フレームごとにNの更新および
コンボリューション積を計算するために(N+1)サイ
クルのクロック信号を必要とするもの、 を具備する請求項7に記載の方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US283,101 | 1988-12-12 | ||
| US07/283,101 US4947363A (en) | 1988-12-12 | 1988-12-12 | Pipelined processor for implementing the least-mean-squares algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH02214322A true JPH02214322A (ja) | 1990-08-27 |
| JP2960086B2 JP2960086B2 (ja) | 1999-10-06 |
Family
ID=23084522
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1321265A Expired - Fee Related JP2960086B2 (ja) | 1988-12-12 | 1989-12-11 | 最小2乗平均アルゴリズムを実施するためのパイプライン式プロセッサ |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4947363A (ja) |
| EP (1) | EP0373468B1 (ja) |
| JP (1) | JP2960086B2 (ja) |
| DE (1) | DE68926154T2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6745218B1 (en) | 1999-03-16 | 2004-06-01 | Matsushita Electric Industrial Co., Ltd. | Adaptive digital filter |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5230007A (en) * | 1991-06-28 | 1993-07-20 | Motorola, Inc. | Method for optimizing an adaptive filter update coefficient |
| US5450339A (en) * | 1991-10-10 | 1995-09-12 | Harris Corp | Noncanonic fully systolic LMS adaptive architecture |
| US5535150A (en) * | 1993-04-20 | 1996-07-09 | Massachusetts Institute Of Technology | Single chip adaptive filter utilizing updatable weighting techniques |
| US5517346A (en) * | 1994-09-16 | 1996-05-14 | Varian Associates | Spectral modification through phase modulation with spatial extent |
| GB9418755D0 (en) * | 1994-09-16 | 1994-11-02 | Ionica L3 Limited | Filter |
| US5907497A (en) * | 1995-12-28 | 1999-05-25 | Lucent Technologies Inc. | Update block for an adaptive equalizer filter configuration |
| US6125438A (en) * | 1997-04-21 | 2000-09-26 | Matsushita Electrical Industrial Co., Ltd. | Data processor |
| US6285412B1 (en) | 1997-07-23 | 2001-09-04 | Harris Corporation | Adaptive pre-equalization apparatus for correcting linear distortion of a non-ideal data transmission system |
| US6519010B2 (en) | 1998-06-26 | 2003-02-11 | Harris Corporation | Broadcast transmission system with sampling and correction arrangement for correcting distortion caused by amplifying and signal conditioning components |
| US6865588B2 (en) * | 2002-01-03 | 2005-03-08 | Intel Corporation | Adaptive filtering with tap leakage using error filtering |
| EP1841065A1 (en) * | 2006-03-29 | 2007-10-03 | Mitel Networks Corporation | Modified least-mean-squares method with reduced computational complexity |
| US8352215B2 (en) * | 2009-10-23 | 2013-01-08 | Sas Institute Inc. | Computer-implemented distributed iteratively reweighted least squares system and method |
| CN102394592B (zh) * | 2011-10-18 | 2013-12-11 | 北京理工大学 | 一种基于Backlash算子的自适应滤波器 |
| US10560313B2 (en) | 2018-06-26 | 2020-02-11 | Sas Institute Inc. | Pipeline system for time-series data forecasting |
| US10685283B2 (en) | 2018-06-26 | 2020-06-16 | Sas Institute Inc. | Demand classification based pipeline system for time-series data forecasting |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4321686A (en) * | 1980-01-24 | 1982-03-23 | Communications Satellite Corporation | Correction processor of self-adaptive filters |
| JPS61102813A (ja) * | 1984-10-25 | 1986-05-21 | Hitachi Denshi Ltd | 適応デイジタルフイルタ |
| US4829463A (en) * | 1985-03-27 | 1989-05-09 | Akai Electric Co. Ltd. | Programmed time-changing coefficient digital filter |
| US4802111A (en) * | 1986-03-10 | 1989-01-31 | Zoran Corporation | Cascadable digital filter processor employing moving coefficients |
| US4726036A (en) * | 1987-03-26 | 1988-02-16 | Unisys Corporation | Digital adaptive filter for a high throughput digital adaptive processor |
-
1988
- 1988-12-12 US US07/283,101 patent/US4947363A/en not_active Expired - Fee Related
-
1989
- 1989-12-04 DE DE68926154T patent/DE68926154T2/de not_active Expired - Fee Related
- 1989-12-04 EP EP89122308A patent/EP0373468B1/en not_active Expired - Lifetime
- 1989-12-11 JP JP1321265A patent/JP2960086B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6745218B1 (en) | 1999-03-16 | 2004-06-01 | Matsushita Electric Industrial Co., Ltd. | Adaptive digital filter |
Also Published As
| Publication number | Publication date |
|---|---|
| EP0373468A2 (en) | 1990-06-20 |
| DE68926154T2 (de) | 1996-10-24 |
| EP0373468A3 (en) | 1991-07-03 |
| JP2960086B2 (ja) | 1999-10-06 |
| EP0373468B1 (en) | 1996-04-03 |
| US4947363A (en) | 1990-08-07 |
| DE68926154D1 (de) | 1996-05-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5287299A (en) | Method and apparatus for implementing a digital filter employing coefficients expressed as sums of 2 to an integer power | |
| EP0042452B1 (en) | Signal processor computing arrangement and method of operating said arrangement | |
| JPH02214322A (ja) | 最小2乗平均アルゴリズムを実施するためのパイプライン式プロセッサ | |
| Hartley | Optimization of canonic signed digit multipliers for filter design | |
| US6665696B2 (en) | Delayed adaptive least-mean-square digital filter | |
| FI97002C (fi) | Suora FIR-suodatin, menetelmä pistetulon laskemiseksi FIR-suodattimessa ja menetelmä suoran FIR-suodattimen suunnittelemiseksi | |
| US6009448A (en) | Pipelined parallel-serial architecture for a modified least mean square adaptive filter | |
| JP2002152014A (ja) | 正規最小平均二乗アルゴリズムに基づいた係数適応用ハードウエアアクセリレータ | |
| JPH039471A (ja) | 移動平均処理装置 | |
| JPS62286307A (ja) | 多重ステージデジタル信号乗算加算装置 | |
| US5367476A (en) | Finite impulse response digital filter | |
| JPS6037513B2 (ja) | デジタル回路 | |
| US4802111A (en) | Cascadable digital filter processor employing moving coefficients | |
| JPH05216627A (ja) | 乗算器および乗算方法 | |
| US4939684A (en) | Simplified processor for digital filter applications | |
| JPH11266140A5 (ja) | ||
| US6788738B1 (en) | Filter accelerator for a digital signal processor | |
| JPS63278411A (ja) | 多段デジタル・フィルタ | |
| JP4388141B2 (ja) | ディジタルフィルタ用共有リソース | |
| US5781462A (en) | Multiplier circuitry with improved storage and transfer of booth control coefficients | |
| US5650952A (en) | Circuit arrangement for forming the sum of products | |
| KR100337716B1 (ko) | 곱의합을형성하는회로 | |
| JP2864597B2 (ja) | ディジタル演算回路 | |
| JPS61213926A (ja) | Dsp演算処理方式 | |
| SU1562904A1 (ru) | Устройство дл умножени на коэффициенты |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070730 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080730 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090730 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090730 Year of fee payment: 10 |
|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: R3D03 |
|
| LAPS | Cancellation because of no payment of annual fees |