JPH06196969A - 有限インパルス応答フィルタ - Google Patents

有限インパルス応答フィルタ

Info

Publication number
JPH06196969A
JPH06196969A JP5219131A JP21913193A JPH06196969A JP H06196969 A JPH06196969 A JP H06196969A JP 5219131 A JP5219131 A JP 5219131A JP 21913193 A JP21913193 A JP 21913193A JP H06196969 A JPH06196969 A JP H06196969A
Authority
JP
Japan
Prior art keywords
sum
input
filter
output
subtractor
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
JP5219131A
Other languages
English (en)
Other versions
JP3041563B2 (ja
Inventor
Aaru Surabuin Kiisu
キース・アール・スラヴィン
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.)
Tektronix Japan Ltd
Original Assignee
Sony Tektronix 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 Sony Tektronix Corp filed Critical Sony Tektronix Corp
Publication of JPH06196969A publication Critical patent/JPH06196969A/ja
Application granted granted Critical
Publication of JP3041563B2 publication Critical patent/JP3041563B2/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
    • H03H17/06Non-recursive 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)【要約】 【目的】 FIRフィルタのゲート数を減らし、集積回
路チップ上で必要となるシリコンの量を低減するととも
に高速化する。 【構成】 2乗の項を展開する恒等式を利用したFIR
フィルタの畳み込みの式を実現するハードウエアを構成
する。複数のフィルタ段14に入力信号に係数を加算し
て2乗しその総和を求め、減算器24において入力信号
の2乗の和との差をとる。この結果と係数の2乗の和と
の差を出力減算器26において取る。この結果を2で割
れば、求める出力信号z(t)が得られる。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、電気信号のデジタル的
な濾波(フィルタリング)に関し、特に有限インパルス
応答(FIR)フィルタに関する。
【0002】
【従来の技術及び発明が解決しようとする課題】標準的
な有限インパルス応答(FIR)フィルタの畳み込み
は、次の数1を満たす。
【0003】
【数1】 z(t)=SUM{y(t−i)*C(T−1−i),(i,0,T−1)}
【0004】ここでz(t)は時間tにおける濾波され
た出力、CはT段夫々のフィルタ(濾波)係数、y
(t)は時間tにおける入力、SUMはiの範囲を0か
らT−1までとしたときのyとCの積の総和を意味す
る。なお、”*”は積を意味する。このような構成のF
IRフィルタにおいて、ゲート数の複雑さG(M,N)
は被乗数間の交差項の数に比例し、それは乗算のM+N
ビットの結果に対しておよそ次の数2で与えられる。
【0005】
【数2】G(M,N)=11M*N−10N−5M
【0006】ここでNは1以上、Mは2以上である。ゲ
ートの数によって、FIRフィルタに必要な集積回路中
のシリコン領域の広さが決まる。
【0007】そこで本発明の目的は、構成のより簡単な
有限インパルス応答(FIR)フィルタを提供すること
である。本発明の他の目的は、必要なゲート数を減少さ
せ、よって集積回路中のシリコンの量を低減させたFI
Rフィルタを提供することである。
【0008】
【課題を解決するための手段】本発明は、2乗(スクエ
アリング)の技術に基づく有限インパルス応答(FI
R)フィルタを提供する。これは、おおまかに言って2
つの要因(2つのS)に基づき、集積回路の形成に必要
なゲート数を低減している。各フィルタ段14のフィル
タ係数Cが入力信号と加算され、その和が2乗され、各
フィルタ段14からの結果は累算される。入力信号の2
乗の和は、最後のフィルタ段14の出力で減算される。
さらに、全てのフィルタ係数Cの2乗の和で減算され、
その結果を2で割ることにより出力信号zを生成する。
この構成は、従来のバイナリ又は剰余数システム(RN
S)演算のどちらかを用いたハードウエアで実現しても
良い。バイナリ演算は、剰余数が1つのRNSの特殊な
場合である。ここで法は、nをビット数としたときに2
^n(2のn乗)である。
【0009】
【実施例】図1は、本発明による有限インパルス応答
(FIR)フィルタの一実施例のブロック図である。入
力信号y(t)は、2乗回路12及び複数T個のフィル
タ段14に印加される。分離フィルタ係数Cに対応する
信号(係数信号)が各フィルタ段14に印加される。各
フィルタ段14は、入力加算回路141を有し、これに
データ入力線16から入力信号y(t)が入力される。
また、適切な係数Cが第2の入力として印加され、入力
信号に加算される。入力加算回路141の出力は2乗回
路143に入力され、入力信号と係数の和の2乗が生成
される。2乗回路143の出力は、アキュムレータ(累
算器)145に入力される。アキュムレータ145は、
アキュムレータ・レジスタ149とこれに直列な第2加
算回路147を有し、システム・クロックCKでクロッ
クされる。複数あるフィルタ段のアキュムレータ145
は直列に縦続接続されるので、第1フィルタ段14
(0)を除けば、その一方の入力は2乗回路143の出
力で、他方の入力は前段のアキュムレータの出力であ
る。
【0010】最初の2乗回路12からの2乗された入力
信号は、入力アキュムレータ18に入力される。入力ア
キュムレータ18は、加算回路20とレジスタ22を有
し、システム・クロックCKでクロックされる。アキュ
ムレータ18の出力は、このアキュムレータ18の他方
の入力として印加されるので、アキュムレータ18の出
力は入力信号の2乗の和である。ただし、加算回路20
の2つの入力の間にはクロックCKに関して1クロック
の時間的なずれがある。アキュムレータ18の出力は、
第1フィルタ段14(0)のアキュムレータ145の他
の入力としても入力される。そして、アキュムレータ1
8の出力は、最後のフィルタ段14(T−1)のアキュ
ムレータ145からの出力と共に減算回路24に印加さ
れる。ここでの差が出力減算回路26に入力される。全
ての係数Cの2乗の和も出力減算回路26に入力され
る。濾波された(フィルタのかかった)出力信号z
(t)の2倍の出力(2*z(t))が、出力減算回路
26の出力として得られる。これを割算回路(図示せ
ず)により2で割れば、求めるz(t)が得られる。濾
波されたデジタル出力信号を生成するには、出力減算回
路26の出力の最小位ビット(LSB)を捨てれば良
い。
【0011】この構成によれば2乗の際に交差項が重複
しているので、必要なゲート数は従来設計に比べて半分
で済む。Nビットの2乗器(結果は2Nビット)の場
合、ゲート数G(N)は、次の数3で与えれる。
【0012】
【数3】 G(N)=11*N^2/2−31*N/2−5
【0013】これは、MとNが大きく且つほぼ等しい場
合のG(M,N)のほぼ半分に大変近い。なお、”^”
は累乗を意味する。本発明では、数4に示す恒等式を利
用している。
【0014】
【数4】 y(i)*C(j)={(y(i)+C(j))^2 −y(i)^2−C(j)^2}/2
【0015】さらには、ハードウエアの付加は最小限に
してy^2の和を算出する方法を用い、y(i)*C
(j)ではなく(y(i)+C(j))^2の和を算出
する。C(j)^2も、最小限ハードウエアの付加で非
実時間軸形式で算出しても良いし、又は予め算出してお
き各係数の組み合わせをロードした様にして外部から取
り入れても良い。
【0016】従来のTタップ(段)FIR畳み込みの式
は次の数5で表される。
【0017】
【数5】 z(t)=SUM{y(t−i)*C(i),(i,0,T−1)}
【0018】2乗の手法に基づくTタップの設計では、
0からTのnに対してT+1個の蓄積要素S(n,t)
の組がある。ハードウエア上では、これら蓄積要素と係
数は次の関係で結合される。
【0019】
【数6】 z(t)={S(T,t)−S(0,t)−K}/2
【数7】 K=SUM{C(i)^2,(i,0,T−1)}
【0020】図1の一実施例のおいては、S(T,t)
に対応する信号は複数フィルタ段14で生成され、S
(0,t)に対応する信号は入力アキュムレータ18か
ら出力されている。また、Kに対応する信号は出力減算
器26のB入力に印加されている。2つの帰納式Sは次
の数8及び数9で示される。
【数8】 S(n,t)= S(n−1,t−1)+(y(t)+C(T−n))^2 (n>0のとき)
【数9】 S(n,t)= S(n,t−1)+y(t)^2 (n=0のとき)
【0021】ここでy(t)は時間tにおけるフィルタ
への入力である。n=0を数9に代入してr回反復適用
すれば数10が得られる。
【0022】
【数10】 S(0,t)=S(0,t−r) +SUM{y(t−i)^2,(i,0,r−1)}
【0023】数10においてtをt−jに、rをm−j
に置き換えると数11を得る。
【0024】
【数11】 S(0,t−j)=S(0,t−m) +SUM{y(t−j−i)^2,(i,0,m−j−1)} =S(0,t−m) +SUM{y(t−i)^2,(i,j,m−1)}
【0025】数11は、最後のm−jの入力信号(サン
プル)の2乗の和に任意の定数S(0,t−m)を加え
たものである。数8の帰納式をn=Tで開始してT回用
いると数12が得られる。
【0026】
【数12】 S(T,t)=S(0,t−T) +SUM{(y(t−i)+C(i))^2,(i,0,T−1)}
【0027】2乗の項を展開し、j=Tとして数11を
S(0,t−T)に代入すると、
【0028】
【数13】 S(T,t)=S(0,t−m) +SUM{y(t−i)^2,(i,T,m−1)} +SUM{y(t−i)^2,(i,0,T−1)} +2*SUM{y(t−i)*C(i),(i,0,T−1)} +SUM{C(i)^2,(i,0,T−1)}
【0029】数13中のy(t−i)^2に関する和を
まとめると、
【0030】
【数14】 S(T,t)=S(0,t−m) +SUM{y(t−i)^2,(i,0,m−1)} +2*SUM{y(t−i)*C(i),(i,0,T−1)} +SUM{C(i)^2,(i,0,T−1)}
【0031】数7をKに、数14をS(T,t)に、そ
してj=0として数11をS(0,t)に夫々数6対し
て代入すると数15が得られる。
【0032】
【数15】 z(t)=[{S(0,t−m) +SUM{y(t−i)^2,(i,0,m−1)} +2*SUM{y(t−i)*C(i),(i,0,T−1)} +SUM{C(i)^2,(i,0,T−1)}} −{S(0,t−m) +SUM{y(t−i)^2,(i,0,m−1)}} −SUM{C(i)^2,(i,0,T−1)}]/2 =SUM{y(t−i)*C(i),(i,0,T−1)}
【0033】数15は、従来のTタップFIR畳み込み
表現である。この結果はまた数11における蓄積要素S
(0,t−m)の最初の状態は結果に影響せず、よって
指数nの値を大きくしたときも数9における累積の反復
を清算する必要はない。
【0034】数9は、零でない入力に対してS(0,
t)のオーバーフローが結果的に生じることを示唆して
いる。よって、システムが線形オーバーフロー特性、つ
まり、計算におけるいくつかの段での加算によるオーバ
ーフロー又はアンダーフローは、後で同じ大きさの減算
によって正しい値に回復する(数6を参照)。
【0035】オーバーフローは元来線形なので、剰余数
システム(Residue Number System、RNS)は自ずと
この要求を満たす。法(Modulus又はMod)が小さければ
2乗器も大変小さいので、経済的なハードウエアを実現
でき、動作速度も本来的に大変高速である。ダイナミッ
ク・レンジRは、システム中の全ての法の積で与えられ
る。法が1つしかないRNSの特殊な場合を考えれば、
2乗技術に基づくFIRフィルタは、従来の演算構造を
用いて実現しても良い。選択した法は通常R=2^nの
形式を取る。ここでnは整数、つまり、ビット数であ
る。他の形式では効率が悪くなる傾向にある。後の減算
がなければ元の状態に回復しないので、それ自身のはし
ご(ラダー)状の累算をオーバーフローなく処理するに
は充分なビット数を維持する必要がある。もし次の数1
6に示す関係が保証されれば、どちらの手法でもオーバ
ーフローがない。
【0036】
【数16】 R>=SUM{ABS(C(i)),(i,0,T−1)}*Rin*2
【0037】ここでRinは許容されるダイナミック入
力レンジ、つまり、入力データとフィルタ係数の間で可
能な最大の相関関係を示すRである。ABSは、C
(i)の絶対値を返す関数である。
【0038】図2は、RNSシステムにおける2乗型F
IRフィルタの実施例のより具体的なブロック図を示し
ている。デジタル入力信号はRNS分配器(スプリッ
タ)30に入力され、デジタル入力信号が選択された法
(Modulus又はMod)に応じて対応する剰余数形式に分離
される。各剰余数信号は、図1において説明した如くF
IRフィルタ32で夫々独立に処理される。独立処理し
た信号はRNS混合器34に入力され、バイナリ(2
値)形式のデジタル出力信号が生成される。簡素なRN
Sシステムの1例を以下に説明する。RNSシステムの
詳細については、IEEE(アメリカ電子・電気通信学
会)出版の再版シリーズ「剰余数システム演算(Residu
e Number System Arithmetic):現代デジタル信号処理
の応用(Modern Application in Digital Signal Proce
ssing)」(1986)を参照されたい。
【0039】0から14までの範囲(レンジ)をバイナ
リ(2値)で扱うには4ビット必要である。3又は5を
法として用いれば、夫々2ビット又は3ビットのシステ
ムが必要となる。以下の表1は、バイナリ値と剰余値と
の間の変換を示している。
【0040】
【表1】
【0041】法演算を用いれば、2乗回路に入力するビ
ット数は2乗回路から出力されるビット数と等しい。こ
のRNSシステムで示されるバイナリ値の範囲は、2つ
の法の積、つまり、15(=3*5)に等しい。より多
い法を用いれば、処理できる値の範囲は拡大する。
【0042】
【発明の効果】上述のように本発明は、2乗の手法に基
づくFIRフィルタを提供する。これは、高速でゲート
数が少なく、よって実施の際に集積回路チップ上に必要
となるシリコンの量を低減することができる。
【図面の簡単な説明】
【図1】本発明によるFIRフィルタの一実施例のブロ
ック図である。
【図2】本発明によるFIRフィルタのより具体的な実
施例のブロック図である。
【符号の説明】 12 2乗回路 14 フィルタ段 16 データ入力線 18 入力アキュムレータ 20 加算回路 22 レジスタ 24 減算回路 26 出力減算回路 30 分配器 32 モジュール 34 混合器 141 入力加算回路 143 2乗回路 145 アキュムレータ 147 第2加算回路 149 アキュムレータ・レジスタ

Claims (1)

    【特許請求の範囲】
  1. 【請求項1】 入力信号の2乗を出力する2乗手段と、 上記入力信号の2乗を上記2乗手段の出力と加算する加
    算手段と、 上記入力信号と係数信号の和の2乗の和を出力する累算
    手段と、 上記加算手段の出力から上記累算手段の出力を減算する
    第1減算手段と、 第1減算手段から上記係数信号の2乗の和を減算する第
    2減算手段とを具えることを特徴とする有限インパルス
    応答フィルタ。
JP5219131A 1992-08-13 1993-08-11 有限インパルス応答フィルタ Expired - Fee Related JP3041563B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US07/929,039 US5561616A (en) 1992-08-13 1992-08-13 Fir filter based upon squaring
US929039 1992-08-13

Publications (2)

Publication Number Publication Date
JPH06196969A true JPH06196969A (ja) 1994-07-15
JP3041563B2 JP3041563B2 (ja) 2000-05-15

Family

ID=25457226

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5219131A Expired - Fee Related JP3041563B2 (ja) 1992-08-13 1993-08-11 有限インパルス応答フィルタ

Country Status (2)

Country Link
US (1) US5561616A (ja)
JP (1) JP3041563B2 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5999656A (en) * 1997-01-17 1999-12-07 Ricoh Co., Ltd. Overlapped reversible transforms for unified lossless/lossy compression
US6233593B1 (en) 1998-11-18 2001-05-15 Tektronix, Inc. Fir filter based upon squaring
US7476777B2 (en) * 2002-09-17 2009-01-13 Ceres, Inc. Biological containment system

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5444547A (en) * 1977-09-16 1979-04-09 Fuji Xerox Co Ltd Heat fixing roll for electrophotography

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4817025A (en) * 1984-02-03 1989-03-28 Sharp Kabushiki Kaisha Digital filter
JPH0828649B2 (ja) * 1989-02-16 1996-03-21 日本電気株式会社 ディジタルフィルタ
JPH10294646A (ja) * 1990-02-16 1998-11-04 Sony Corp サンプリングレート変換装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5444547A (en) * 1977-09-16 1979-04-09 Fuji Xerox Co Ltd Heat fixing roll for electrophotography

Also Published As

Publication number Publication date
JP3041563B2 (ja) 2000-05-15
US5561616A (en) 1996-10-01

Similar Documents

Publication Publication Date Title
US3665171A (en) Nonrecursive digital filter apparatus employing delayedadd configuration
EP0948133A2 (en) Digital filter with efficient quantization circuitry
JPS6360927B2 (ja)
JPH0656942B2 (ja) デジタル信号用遅延回路
US6304133B1 (en) Moving average filter
US5798954A (en) Digital filter device having a bit shifter unit
JPH06196969A (ja) 有限インパルス応答フィルタ
US8090013B2 (en) Method and system of providing a high speed Tomlinson-Harashima Precoder
EP0791242B1 (en) Improved digital filter
Lehto et al. Synthesis of narrowband linear-phase FIR filters with a piecewise-polynomial impulse response
US6058404A (en) Apparatus and method for a class of IIR/FIR filters
Ohlsson et al. Implementation of low complexity FIR filters using a minimum spanning tree
Shanthi et al. HIGH SPEED AND AREA EFFICIENT FPGA IMPLEMENTATION OF FIR FILTER USING DISTRIBUTED ARITHMETIC.
WO2005002051A1 (ja) デジタルフィルタ
Narasimha et al. Implementation of LOW Area and Power Efficient Architectures of Digital FIR filters
JP3090043B2 (ja) ディジタル補間フィルタ回路
JPH0342729B2 (ja)
EP0820145A2 (en) Interpolation filter
US6757702B1 (en) Adaptive filter
Living et al. High performance distributed arithmetic FPGA decimators for video-frequency applications
JPH036690B2 (ja)
JPS6165362A (ja) デジタル信号処理装置
JPS6097743A (ja) 適応型線形予測装置
JPH08335850A (ja) 簡易デジタルフィルタ
JPS5833754A (ja) デイジタル乗算回路

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees