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
Application number
JP1321265A
Other languages
English (en)
Other versions
JP2960086B2 (ja
Inventor
Tim A Williams
ティム・エー・ウィリアムズ
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.)
Motorola Solutions Inc
Original Assignee
Motorola Inc
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 Motorola Inc filed Critical Motorola Inc
Publication of JPH02214322A publication Critical patent/JPH02214322A/ja
Application granted granted Critical
Publication of JP2960086B2 publication Critical patent/JP2960086B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00Adaptive networks
    • H03H21/0012Digital 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

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 [産業上の利用分野] 本発明は、−数的にはデータプロセッサに関し、より特
定的には適応フィルタ型データプロセッサに関する。
[従来の技術] 一般に、伝達関数はリニアな、時間不変システム、ある
いは時間変化、適応型システムを表わす。
リニアな時間不変システムに関しては、所定のインパル
ス関数が励起およびインパルス応答結果のために用いら
れる。何らかの池の入力波形が与えられた場合、インパ
ルス応答を入力波形とコンボルブすることにより該イン
パルス応答を出力を決定するために利用することができ
る。これに対して、適応型システムにおいては、インパ
ルス関数は各データサンプル点においであるいはデータ
サンプルの各組において変化している。従って、適応型
システムは時間変化しておりあるいはノンリニアである
。デジタルフィルタがフィルタの所望の関数に応じて適
応プロセスにより特定の伝達関数に向けられている。デ
ジタルフィルタの1つの応用は所定の伝達関数により1
組のデータをろ波する場合における最小2乗平均誤差を
得ることである。この関数は最小2乗平均(L M S
 : Ieast−mean−square )アルゴ
リズムとして知られている。
この場合のエラーは何らかの所望の目標とフィルタか実
際に出力しているものとの間の誤差を表わしかつフィル
タの誤調整の測定値となる。LMSアルゴリズムは、(
1)部分的に先のエラーに基づくフィルタの係数の更新
である適応計算、(2)所定の係数値およびデータ値の
積の累積であるコンボリューション計算、(3)そして
エラー計算を含む、エラー計算はコンボリューション計
算の出力値を所定の所望値から減算することである。
商業的に入手可能な適応フィルタ型プロセッサは富士退
社のMB87XXAFPが例にあげられ、これはLMS
適応FIRフィルタアルゴリズムをサポートすることが
可能である。知られたLMSアルゴリズムシステムにお
いては、コンボリューション計算が最初に行なわれそこ
からエラーが決定される0次に前記エラーが利用される
、更新計算がなされる。非常に高度な多重タップのフィ
ルタに対しては、このことはすべてのタップに対するコ
ンボリューション計算が最初に計算されなければならず
、そして次にすべてのタップに対するエラーおよび更新
計算がなされなければならない。
[発明が解決しようとする課題〕 この知られたLMSアルゴリズムの設備に対する不都合
は多重タップフィルタを実施するために非常に多くのプ
ロセッサのクロックサイクルが必要となり、係数、デー
タおよびエラー値を格納するために多重メモリ構造を使
用する必要があり、そしてパイプライン式多重かつ分離
されたフィルタにおける効率の悪さを含む、典型的には
、パイプライン様式で多重フィルタ機能を実施する従来
のシステムはカスケード構造で動作する個々の集積回路
製品を使用した。
従って、本発明の目的は、最小平均2乗(LMS)アル
ゴリズムを実施するための改良されたパイプライン式プ
ロセッサを提供することにある。
本発明の他の目的は、適応型フィルタを実施するための
改良されたデジタルフィルタ構造を提供することにある
本発明のさらに他の目的は、最小平均2乗(LMS)ア
ルゴリズムを実施するための改良された方法を提供する
ことにある。
し課題を解決するための手段および作用]本発明におけ
るこれらおよび他の目的を実施するにあたり、1つの形
式として、Nを整数とした場合に、Nタップのフィルタ
における最小2乗平均(LMS)アルゴリズムを実施す
るための方法およびパイプライン式プロセッサが提供さ
れる。
一連のデータ値および係数値のコンボリューションがN
のタップに対して計算される。実質的に同時に、フィル
タの各タップに対する係数値を更新することによりプロ
セッサは一連の係数値の適応計算を連成する。パイプラ
イン式プロセッサの更新部は所定の時間期間の開所定数
の更新計算を行なう、コンボリューション部が更新部に
結合され同じ時間期間の開所定数のコンボリューション
積を連成する。記憶部が更新部およびコンボリューショ
ン部に結合され最小平均2乗アルゴリズムを実施するた
めに更新部およびコンボリューション部の双方によって
要求されるデータ値および係数値を選択的に提供する。
これらおよび他の目的、特徴および利点は添付の図面と
組合わせて以下の詳細な説明からより明瞭に理解できる
であろう。
[実施例] 一般的に、LMSアルゴリズムは2つの相互依存する方
程式によって構成される。2つの方程式の内の1つは典
型的には更新(update)方程式と称されかつ他方
の方程式はコンボリューション方程式として言及される
。更新方程式はデジタルフィルタの適応部分であり伝統
的なニュートンーラルフソン反復アルゴリズムと同様の
様式でフィルタ係数を修正する。コンボリューション方
程式は伝統的なFIRフィルタ方程式である。双方の方
程式は各データサンプル期間°にフィルタの各タップに
おいて連成される。
第1図にはLMSアルゴリズムによって要求されるコン
ボリューションを実施するための時間不変フィルタ構造
10が示されている。このコンボリューションは次の方
程式で表わされる。
N=1 Y(k)=ΣX  (j、k)*W(J、  k)J=
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
の入力に接続されている。
遅延回路14の出力はタッグ2および時間kにおけるデ
ータ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は所
定の係数値およびタップデータ値の積を表わす個々の積
を受ける。加算器回路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)は固定された、変化
しない係数でありこれがインパルス応答が時間的に不変
であるという理由である。
第2図には知られたフィルタ構造20が示されており、
この構造は方程式1で表わ、されるように、適応計算お
よびコンボリューション演算の双方が行なわれる適応型
または時間変化型フィルタの単一フィルタタップを表わ
す、言替えれば、第1図の係数値W(0)、W(1)お
よびW(2)が各時間期間に形成され時間的に変化する
。このL MSアルゴリズムに対する適応計算は次のよ
うに表わされる。
W(J、 k)=w(、J、 k−1)十βme(k−
1)11X(j、に−1)   ・・−(2)上式にお
いて、j=0.・・・、N−1であり、かつに=o、1
,2.・・・、■。
遅延回路21はデータ入力X(1,k)を受けるための
入力を有する。遅延回路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の累積を与える。
動作においては、データ人力X(1,k)がタップ位置
1から遅延回路21によって受信される。
遅延回路21はデータ人力X(2,k)を提供しこれは
乗算器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のような構造が各タップに対し必要
とされるから、マルチタップのフィルタに対し要求され
る回路の量は非常に大きい。
第3図には第2図の構造20の変形であるフィルタ構造
30が示されている。構造30においては、データ人力
X(1,k)が遅延回路31の入力に接続されている。
遅延回路31の出力はデータ人力X(2,k)を提供し
これは遅延回路33の入力および乗算器37の第1の入
力の双方に結合されている。遅延回路33の出力は遅延
されたデータ人力X(2,に−1>を提供しかつ乗算器
34の第1の入力に接続されている1乗算器34の第2
の入力はエラー環βe(k−1)を受ける。
乗算器34の出力は加算器回路35の第1の入力に接続
されている。加算器回路35の出力は係数値W(2,k
)を提供しかつ遅延回路36の入力および乗算器37の
第2の入力の双方に接続されている。遅延回路36の出
力は係数値W(2,に1)を提供しかつ加算器回路35
の第2の入力に接続されている。乗算器37の出力は加
算器回路38の第1の入力に接続されている。加算器回
路38の第2の入力は論理ゼロの値に接続されている。
加算器回路38の出力は方程式(1)のコンボリューシ
ョンアルゴリズムを実施するうえで信号Y(k)の第1
の累積を提供する。
動作においては、フィルタ構造30は第2図の遅延回路
24の位置を再配置することによりフィルタ構造20と
かなり構造的に異なっているが機能的には等価である。
第2図の遅延回路24の機能を加算器23の前に位置づ
けることによって、同じ機能を連成することができる。
従って、第3図においては、第2図の遅延回路24によ
って実施された遅延機能は遅延回路33および36そし
てエラー環βe(k−1)を乗算器34に結合すること
により連成される。フィルタ構造20のように、フィル
タ構造30は方程式1の適応規則および方程式2のコン
ボリューションアルゴリズムの双方を行なうことによっ
てLMSアルゴリズムを実施する。
第4図にはフィルタ構造30の延長でありがっLMSア
ルゴリズムを実施するための3タツプフイルタであるフ
ィルタ構造40が示されている。
本発明は、任意の数のフィルタタップを実施するために
も利用できかつ第4図における3タツ1の実施例は限定
ではなく説明のためだけにすぎない。
データ人力X(0,k)が遅延回路42の入力および乗
算器44の第1の入力に結合されている。
遅延回路42の出力は遅延された入力データX(1,k
)を提供しかつ乗算器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の入力に接続されている。
加算器回F#158の出力は更新された係数値W(1゜
k)を提供しかつ乗算器52の第2の入力および遅延回
路60の入力の双方に接続されている。遅延回路60の
出力は加算器回路58の第2の入力に接続されている。
遅延回路54の出力はまた乗算器64の第1の入力およ
び遅延回路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の入力に接続されている。
加算器回路76の出力はLMSアルゴリズムのコンボリ
ューション出力信号Y(k)を提供する。
動作においては、時間変化LMSフィルタ槽遣40にお
ける各タップ計算は更新計算(方程式2)およびコンボ
リューション計算(方程式1)を必要とする。これらの
2つの計算を連成するために各々の新しいデータサング
ルに対しフィルタタッグごとに2つの乗算/加算操作が
必要とされる。
フィルタ構造40は乗算および加算の操作のみならず更
新およびコンポリ−シュン計算のパイプライン化をも可
能とする。各加算器および各乗算器はそのそれぞれの計
算を形成するため完全な半サイクルを持つ事を許容され
るであろう、タッフ′データ人力Xおよび係数値Wが示
すように、各タップ位置におけるすべての操作はサイク
ルにのような単一時間サイクルの間に実質的に同時に起
っている0例えば、タップ位置3においては、乗算器6
4および66が時間サイクルにの最初の半分の間に機能
している9時間サイクルにの第2の半分の間に、加算器
68および72が機能している。
同様に、タップ位22の乗算器52および56とタップ
位置1の乗算器44および46もまた時間サイクルにの
最初の半分の間に機能している。タップ位置2の加算器
58および74とタップ位置1の加算器48および76
もまた時間サイクルにの第2の半分の間に機能している
時間サイクルにの最初の半分の間、乗算器64は方程式
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に形成される。
時間サイクルにの第2の半分の間、加算器72は先に形
成されたW(2,k)およびX(2,k)の積を強制さ
れた2進ゼロと加算している。引き続きかつ順次に加算
器74および76は引き続くタッグに対するトータルの
累積を完了し方程式1のコンボリューションアルゴリズ
ムの実行を完了する0時間サイクルにの第2の半分の間
に同時に加算器68は遅延回路70によって提供された
先の係数値W(2,k)を乗算器66により提供された
積と加算し次の時間サイクル(k+1>に対する更新さ
れた係数値W(2,に+1)を提供している。各タップ
位置に対する計算の各々は実質的に同時に起っているが
、計算は順次パイプライン化されかつ右側の最も高いタ
ッグ位置から左側の最も低いタップ位置へと計算できる
ことに注目すべきである。
図示された形式では、第4図に示された3タツプの各タ
ップのフィルタ構造が第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つの乗算器の同時的な再符号化能力は速度、サイズ
およびパイプライン化の利点を提供する。
第5図には所定の整数Nのタッグに対する第4図のフィ
ルタ構造40を実施するパイプライン式グロセッサ79
が示されている。特に、プロセッサ79はLMSアルゴ
リズムを実行するよう機能する。プロセッサ79は双方
向データバス80により通信する。一般に、プロセッサ
79は時間変化する係数値を与えるための更新部81、
コンボリューション部82、そして制御部83を有する
ものと考えることができる。ラッチ85の入力はデータ
バス80に接続されている。ラッチ85の出力はマルチ
プレクサ86の第1の入力に接続されている。制御信号
C4はマルチプレクサ86の制御入力に接続されている
。マルチプレクサ86の第2の入力は2進ゼロの値に接
続されている。
マルチプレクサ86の出力は乗算器87の第1の入力に
#枕されている0乗算器87の出力はラッチ88の第1
の入力に接続されている。ラッチ89の入力はデータバ
ス80に接続されている。ラッチ89の出力はマルチプ
レクサ90の第1の入力に接続されている。rcl、と
名付けられた制御信号はマルチプレクサ90の制御入力
に結合されている。マルチプレクサ90の出力はレコー
ダ(再符号化器)回路92および93の各々の入力に接
続されておりかつラッチ94の入力に接続されている。
レコーダ回路92の第1の出力は乗算器87の第2の入
力に接続されている。レコーダ回路92の第2の出力は
ラッチ88の第2の入力に接続されている。ラッチ88
の第1の出力は「C」と名付けられた信号を提供しかつ
加算器96の第1の入力に接続され、そしてラッチ88
の第2の出力は加算器96の制御入力に接続されている
。加算器96の出力はラッチ97の入力に接続されてい
る。ラッチ97の出力はラッチ98の入力に接続されて
いる。ラッチ94および98の各々の出力はランダムア
クセスメモリ(RAM)99の入力に接続されている。
制御信号「C2」および「C3」はRAM99のそれぞ
れの制御入力に結合されている。RAM99の出力はラ
ッチ101およびラッチ102の双方の入力に接続され
ている。ラッチ102の出力はマルチプレクサ90の第
2の入力に接続されており、かつラッチ101の出力は
ラッチ104の入力に接続されている。ラッチ104の
出力は加算器96の第2の入力に接続されている。
コンボリューション部82においては、乗算器106の
第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に接続されている。
動作においては、プロセッサ79は更新およびコンボリ
ューション操作を連続して連成するためにパイプライン
様式で機能する。プロセッサ7つはまたLMSアルゴリ
ズムを実施するために必要とされる乗算および加算操作
をパイプライン処理する。このようにして計算を組織化
することにより各加算器および乗算器がその結果を形成
するために完全な半サイクルを持つことを許容する。ま
た、単一のRAM構造がタップ重み付けおよびフィルタ
データの双方を保持するために利用されている。係数お
よびデータに対する独立のRAMの使用を除去すること
によりかなりの量の回路領域を節約する。
第6図を参照すると、プロセッサ79の動作かより容易
に理解できる。第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から読まれる。
係数値W(2,に−1)が正しい値であるという理由は
第3のタップ位置において第4図のフィルタ構造40か
ら確認することができる。
第1のクロックサイクルの第2のクロック位相において
は、乗算器87は出力信号C(2,k>を提供するよう
機能しこの出力信号は引き続き加算器96によって使用
される。デー°少入力X(3゜k)がRAM99からマ
ルチプレクサ90を介してレコーダ92に結合される。
エラー項βe(k−1)がデータバスからラッチ85お
よびマルチプレクサ86を介して提供される。!&初は
、マルチプレクサ86は制御信号C4によって制御され
エラー更新がないことを表わすゼロ値入力を提供する。
LMSアルゴリズムのY(k)の引き続く計算のために
、図示しない回路か乗算器87による使用のためエラー
項βe(k−1)を提供する。レコーダ92は乗算器8
7による使用のための所定の再符号化アルゴリズムを実
施するよう11能する。
レコーダ92と共に付加的なレコーダを用いてより高い
精度のデータに備えることができることを理解すべきで
ある。第1のクロックサイクルにおいては、乗算器10
6は出力を提供しないが、これは加算器96がいまだW
オペランドに対する値を提供していないからである。第
1のクロックサイクルにおいてRAM99のアドレスA
がアドレスされ、かつクロックサイクルφlの間にRA
M99が読み取られることに注目すべきである。また第
1のクロックサイクルの第2のクロック位相の間ラッチ
114はゼロにリセットされる。
第2のクロックサイクルの第1のクロック位相において
は、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から読ま
れる。
第2のクロックサイクルの第2のクロック位相において
は、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に応答して選
択される。
第3のクロックサイクルの第1のクロック位相において
は、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に結合される。
同時に、コンボリューション部82の加算器110が乗
算器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)を計算する。
第4のクロックサイクルの第1のクロック位相において
は、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に出力される。
図示しない回路が次に出力されたY(k)の値を利用し
かつこの値を所定の正しい値から減算して後続の(k+
1)時間フレームにおいて使用されるべき新しいエラー
環βe (k)を計算する。
第4のクロックサイクルの第2のクロック位相において
は、RAM99が第1のクロック位相で計算された係数
値W(0,k)を格納するために書き込まれる。レコー
ダ93にはデータ入力×(0,k)が与えられ、かつ乗
算器106には係数IW(0,k>が与えられる。更新
部81の乗算器87は第4のクロックサイクルの間は機
能しないが、それは最低のタップ位置、ゼロ、に対する
Cオペランドは先に計算されたからである。第4のクロ
ックサイクルの間、RAM99のアドレスA位置は読み
出しおよび書き込みのため制御信号C2に応答して再び
選択され得る。
第4のクロックサイクルが終了した後、プロセッサ79
の動作は繰返される。即ち、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が実施できることを意味する。
本発明の他の見地はLMSアルゴリズムを実施する多重
フィルタのパイプライン式動作を含む。
典型的には、他の者がNタップフィルタを用いあるプロ
セッサまたは集積回路でNのコンボリューション操作を
完全に連成しかつ次にコンボリューションおよび更新の
双方の操作を行なうな、めの他のプロセッサまたは単一
ハードウェア回路の時分割多重によりNの更新操作を行
なうことによりLMSアルゴリズムを実施している。こ
の場合、Nの更新操作に関連するオペランドの各々は更
新操作のそれぞれのオペランドに選択的に結合されなけ
ればならなかった。これに対し、本発明においては、更
新構造は本質的にコンボリューション構造において与え
られる。言い換えれば、本発明におけるNタップフィル
タはコンボリューションおよび更新の双方をN回計算す
ることにより連成して実施される。これはプロセッサの
実行時間をかなり短縮しかつNタップフィルタを実施す
るために使用される回路の量を最小化する。多重フィル
タがこのようにして独立のフィルタ間の遷移の間クロッ
クサイクルの損失なしに実施できる6本発明のプロセッ
サ構造は制御信号C4を介しマルチプレクサ86を制御
可能とすることにより係数更新操作の外部制御を提供す
る。いずれの新しいデータサンプル対しても、外部制御
回路はフィルタの係数が更新されているか否かを決定す
ることかできる。伝統的な外・部制御回路もまたマルチ
プレクサ86において適応ベータ(β)を選択すること
により更新レートを決定することができる。外部回路は
またもしそれが望ましい場合にはタップデータを知られ
た状態にすることができる。
本発明が好ましい実施例に関連して説明されたが、当業
者には本発明は多くの方法で変形することができかつ上
に説明されたもの以外の多くの実施例を含むことができ
ることは明らかである。従って、添付の請求の範囲は本
発明の真の精神および範囲内にある発明のすべての変形
をも含むことを意図している。
【図面の簡単な説明】
第1図は、LMSアルゴリズムに対するコンボリューシ
ョン方程式を実施するための知られたフィルタ構造を示
すブロック図、 第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. 【特許請求の範囲】 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に記載の方法。
JP1321265A 1988-12-12 1989-12-11 最小2乗平均アルゴリズムを実施するためのパイプライン式プロセッサ Expired - Fee Related JP2960086B2 (ja)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Cited By (1)

* Cited by examiner, † Cited by third party
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