JPH0740656B2 - 無限長インパルス応答フイルタ - Google Patents
無限長インパルス応答フイルタInfo
- Publication number
- JPH0740656B2 JPH0740656B2 JP61227028A JP22702886A JPH0740656B2 JP H0740656 B2 JPH0740656 B2 JP H0740656B2 JP 61227028 A JP61227028 A JP 61227028A JP 22702886 A JP22702886 A JP 22702886A JP H0740656 B2 JPH0740656 B2 JP H0740656B2
- Authority
- JP
- Japan
- Prior art keywords
- filter
- words
- equation
- signal
- word
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
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
- H03H17/04—Recursive filters
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Complex Calculations (AREA)
- Networks Using Active Elements (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
- Input Circuits Of Receivers And Coupling Of Receivers And Audio Equipment (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Electrophonic Musical Instruments (AREA)
- Analogue/Digital Conversion (AREA)
Description
【発明の詳細な説明】 〔産業上の利用分野〕 この発明は、無限長のインパルス応答、(以下、IIRと
略称する)フィルタに関する。
略称する)フィルタに関する。
当業者には周知の如く、IIRフィルタはディジタルフィ
ルタの一種であり、時間的に間隔Tだけ隔置された一連
のディジタルワードから成る(例えばアナログ信号のサ
ンプル値から成ってもよい)入力信号の濾波を行うのに
使用できる。フィルタはTに等しい複数の各期間(クロ
ック周期)の間働いて入力信号の相対的に遅延したワー
ドの組を処理して濾波され、出力信号のワードを発生す
る。IIRフィルタは巡回型である点でその他の各種のデ
ィジタルフィルタと区別される。すなわち、各出力ワー
ド(またはサンプル)は過去すなわち前に発生されたワ
ード又はワード群並びに現在及び過去の入力ワードの関
数である。従って、IIRフィルタは少くとも1ワードを
他の信号と加算する加算手段に帰還する帰還ループを有
する。帰還ループでは帰還されたワードは重み係数を乗
算される。乗算動作は時間を要する。そこ連続した加算
動作を行う。フィルタが動作中ならば、加算及び乗算に
要する全時間はTを越えてはならない。すなわちデータ
レート(ワードまたはサンプル周波数)fsが非常に高す
ぎてT(=1/fs)が上述の乗算及び加算の全時間より大
きい入力信号をフィルタは処理できない。これは最も速
く(且つ最も高価な)乗算器及び加算手段が十分に速く
ない場合でも高速の入力信号を処理することを意図した
フィルタの設計者にとって厳しい制約を課することにな
る。
ルタの一種であり、時間的に間隔Tだけ隔置された一連
のディジタルワードから成る(例えばアナログ信号のサ
ンプル値から成ってもよい)入力信号の濾波を行うのに
使用できる。フィルタはTに等しい複数の各期間(クロ
ック周期)の間働いて入力信号の相対的に遅延したワー
ドの組を処理して濾波され、出力信号のワードを発生す
る。IIRフィルタは巡回型である点でその他の各種のデ
ィジタルフィルタと区別される。すなわち、各出力ワー
ド(またはサンプル)は過去すなわち前に発生されたワ
ード又はワード群並びに現在及び過去の入力ワードの関
数である。従って、IIRフィルタは少くとも1ワードを
他の信号と加算する加算手段に帰還する帰還ループを有
する。帰還ループでは帰還されたワードは重み係数を乗
算される。乗算動作は時間を要する。そこ連続した加算
動作を行う。フィルタが動作中ならば、加算及び乗算に
要する全時間はTを越えてはならない。すなわちデータ
レート(ワードまたはサンプル周波数)fsが非常に高す
ぎてT(=1/fs)が上述の乗算及び加算の全時間より大
きい入力信号をフィルタは処理できない。これは最も速
く(且つ最も高価な)乗算器及び加算手段が十分に速く
ない場合でも高速の入力信号を処理することを意図した
フィルタの設計者にとって厳しい制約を課することにな
る。
従って、この発明の目的は改善されたIIRフィルタを提
供するにある。
供するにある。
この発明の他の目的はその速度が余り早すぎて従来のII
Rフィルタでは処理出来ない入力信号を処理することが
できるIIRフィルタを提供するにある。
Rフィルタでは処理出来ない入力信号を処理することが
できるIIRフィルタを提供するにある。
この発明の更に他の目的は従来のフィルタでは入力信号
の速度すなわちデータレートが所定値に制限されていた
けれども、所定値よりかなり大きな値の速度すなわちデ
ータレートを有する入力信号を処理できる乗算器及び加
算手段を用いたIIRフィルタを提供するにある。
の速度すなわちデータレートが所定値に制限されていた
けれども、所定値よりかなり大きな値の速度すなわちデ
ータレートを有する入力信号を処理できる乗算器及び加
算手段を用いたIIRフィルタを提供するにある。
またこの発明の他の目的は高速入力信号を処理できしか
も比較的低速(従って比較的廉価)な部品を使用できる
IIRフィルタを提供するにある。
も比較的低速(従って比較的廉価)な部品を使用できる
IIRフィルタを提供するにある。
この発明は時間的にTだけ離れた一連のディジタルワー
ドから成る入力信号を濾波するIIRフィルタを提供す
る。このフィルタはp・T(但しpは2の倍数である整
数)に等しい複数の各期間の間夫々動作し、入力信号の
隣接ワードの各組(各組は相互に1ワードだけオフセッ
トされている)を処理してp・Tに等しい各期間の間p
個の出力ワードを発生するp(例えば2)個のIIRフィ
ルタ部を有する。またフィルタはp・Tに等しい各期間
中に発生されたp個の出力ワードをマルチプレクスして
時間的に間隔Tだけ離れた一連の出力ワードから成る濾
波出力信号を形成するマルチプレクスを有する。p個の
処理ユニットがp・Tに等しい期間中にp個の出力ワー
ドを発生するように働く、すなわち各処理ユニットがp
で分割された入力信号のデータレートに等しいデータレ
ートで動作できることによって、各出力ワードを形成す
るのに利用できる時間は、pに等しい関数だけ(上述し
た従来のフィルタに比較して)増大する。従って、同じ
動作速度の処理回転(乗算器及び加算器)を使用すれ
ば、このIIRフィルタの最大動作速度は上述した従来の
タイプのIIRフィルタに対してpの関数だけ増大する。
ドから成る入力信号を濾波するIIRフィルタを提供す
る。このフィルタはp・T(但しpは2の倍数である整
数)に等しい複数の各期間の間夫々動作し、入力信号の
隣接ワードの各組(各組は相互に1ワードだけオフセッ
トされている)を処理してp・Tに等しい各期間の間p
個の出力ワードを発生するp(例えば2)個のIIRフィ
ルタ部を有する。またフィルタはp・Tに等しい各期間
中に発生されたp個の出力ワードをマルチプレクスして
時間的に間隔Tだけ離れた一連の出力ワードから成る濾
波出力信号を形成するマルチプレクスを有する。p個の
処理ユニットがp・Tに等しい期間中にp個の出力ワー
ドを発生するように働く、すなわち各処理ユニットがp
で分割された入力信号のデータレートに等しいデータレ
ートで動作できることによって、各出力ワードを形成す
るのに利用できる時間は、pに等しい関数だけ(上述し
た従来のフィルタに比較して)増大する。従って、同じ
動作速度の処理回転(乗算器及び加算器)を使用すれ
ば、このIIRフィルタの最大動作速度は上述した従来の
タイプのIIRフィルタに対してpの関数だけ増大する。
この発明の上述した目的、その他の目的、特徴及び利点
は以下に添付図面を参照して詳述する実施例の説明から
明らかになろう。
は以下に添付図面を参照して詳述する実施例の説明から
明らかになろう。
ディジタルフィルタの設計は当業者には周知の如く(エ
ア・アール・ラビナ(L.R.Rabiner)及びビー・ゴール
ド(B.GOLD)著、1975年プレンチス−ホール(Prentice
−Hall)社発行、“ディジタル信号処理の理論と応用”
を参照)、1極(1次)IIRフィルタのZ伝達特性H
(z)(これから周波数応答を導出できる)は次式によ
って表わすことができる。
ア・アール・ラビナ(L.R.Rabiner)及びビー・ゴール
ド(B.GOLD)著、1975年プレンチス−ホール(Prentice
−Hall)社発行、“ディジタル信号処理の理論と応用”
を参照)、1極(1次)IIRフィルタのZ伝達特性H
(z)(これから周波数応答を導出できる)は次式によ
って表わすことができる。
ただし、 X(z)=ディジタル入力シーケンスすなわち信号x
(n)のz変換、 Y(z)=入力シーケンスすなわち信号x(n)に応答
してフィルタにより発生された濾波出力シーケンスすな
わち信号のz変数、 z=ejωT,すなわちcos(ωT)+jsin(ωT)に
等しい複素指数(但し、ωは入力シーケンスすなわち信
号x(n)に等価なサンプルされたシヌソイド(サンプ
ル間隔T)の角周波数)、 a0,a1,a1=定数(フィルタで行われる乗算動作の重み係
数)、及び Z-1=1サンプリング間隔Tの遅延に対応した遅延演算
子(その間隔は入力信号を作るディジタルワードの時間
間隔に等しい。) 当業者には周知の如く、1次IIRフィルタにより行われ
る動作は、次の1次差分方程式(これは伝達特性を描
く)により表わすことができる。
(n)のz変換、 Y(z)=入力シーケンスすなわち信号x(n)に応答
してフィルタにより発生された濾波出力シーケンスすな
わち信号のz変数、 z=ejωT,すなわちcos(ωT)+jsin(ωT)に
等しい複素指数(但し、ωは入力シーケンスすなわち信
号x(n)に等価なサンプルされたシヌソイド(サンプ
ル間隔T)の角周波数)、 a0,a1,a1=定数(フィルタで行われる乗算動作の重み係
数)、及び Z-1=1サンプリング間隔Tの遅延に対応した遅延演算
子(その間隔は入力信号を作るディジタルワードの時間
間隔に等しい。) 当業者には周知の如く、1次IIRフィルタにより行われ
る動作は、次の1次差分方程式(これは伝達特性を描
く)により表わすことができる。
y(n)=a0x(n)+a1x(n−1) −b1y(n−1)……(2) 式(1)及び(2)は第4図に示すような構成によりハ
ードウエアで実現できる。第4図において、そのz変換
がX(z)である入力信号x(n)は入力端子(10)に
供給される。入力信号は間隔Tを有する一連のディジタ
ルワード(例えば、8ビットワードであるが、必ずしも
これに限らない)である。各ワードは周波数fs(=1/
T)でサンプルされているアナログ信号のディジタル化
されたサンプル値であってもよい。各ワードのビットは
並列に到来するものと考えることができる。また各ワー
ドは第4図の回路間で並列に伝送されるので、第4図
(及びその他の図)において種々のラインは適当なビッ
ト容量を有するバスすなわちハイウエイであることがわ
かる。
ードウエアで実現できる。第4図において、そのz変換
がX(z)である入力信号x(n)は入力端子(10)に
供給される。入力信号は間隔Tを有する一連のディジタ
ルワード(例えば、8ビットワードであるが、必ずしも
これに限らない)である。各ワードは周波数fs(=1/
T)でサンプルされているアナログ信号のディジタル化
されたサンプル値であってもよい。各ワードのビットは
並列に到来するものと考えることができる。また各ワー
ドは第4図の回路間で並列に伝送されるので、第4図
(及びその他の図)において種々のラインは適当なビッ
ト容量を有するバスすなわちハイウエイであることがわ
かる。
入力信号は乗算器(12)へ供給され、ここで各連続する
ワードに重み係数a0が乗算される。また入力信号はZ-1
遅延要素(14)に供給され、ここでTに等しい遅延すな
わち1ワードまたはサンプル周期の遅延を受ける。遅延
要素(14)はワード当りのビットの数に等しい多数のラ
ッチで構成してよく、各ラッチはクロックパルス発生回
路(図示せず)から周波数fs(1/T)でクロックするす
なわち同期信号を供給され、クロックパルス発生回路は
また第4図のその他の回路にも接続され、その動作を円
滑化されるべくそれ等の動作を同期化している。(当業
者には周知の如く、第4図の回路は更に複数のラッチ
(図示せず)を有してもよく、各ラッチはクロック信号
を供給されて適切な動作を達成する。簡単化のため、こ
れ等のラッチは第4図から削減している。
ワードに重み係数a0が乗算される。また入力信号はZ-1
遅延要素(14)に供給され、ここでTに等しい遅延すな
わち1ワードまたはサンプル周期の遅延を受ける。遅延
要素(14)はワード当りのビットの数に等しい多数のラ
ッチで構成してよく、各ラッチはクロックパルス発生回
路(図示せず)から周波数fs(1/T)でクロックするす
なわち同期信号を供給され、クロックパルス発生回路は
また第4図のその他の回路にも接続され、その動作を円
滑化されるべくそれ等の動作を同期化している。(当業
者には周知の如く、第4図の回路は更に複数のラッチ
(図示せず)を有してもよく、各ラッチはクロック信号
を供給されて適切な動作を達成する。簡単化のため、こ
れ等のラッチは第4図から削減している。
遅延要素(14)より発生された出力信号は乗算器(16)
に供給され、ここでその各連続するワードに重み係数a1
が乗算される。
に供給され、ここでその各連続するワードに重み係数a1
が乗算される。
夫々重み係数a0及びa1を乗算された入力信号の一対の隣
接ワード(すなわち現在のワード及び前のワード)に対
応する乗算器(12)及び(16)からの各出力信号は、各
クロック周期の間、加算器または加算手段(18)の各入
力側に供給される。適当に共通接続された複数の個々の
加算ユニットで構成される加算器(18)は出力信号すな
わちシーケンスy(n)を発生し、そのz変換は出力端
子(20)でY(z)である。
接ワード(すなわち現在のワード及び前のワード)に対
応する乗算器(12)及び(16)からの各出力信号は、各
クロック周期の間、加算器または加算手段(18)の各入
力側に供給される。適当に共通接続された複数の個々の
加算ユニットで構成される加算器(18)は出力信号すな
わちシーケンスy(n)を発生し、そのz変換は出力端
子(20)でY(z)である。
上述の如く、IIRフィルタは、各出力ワード(またはサ
ンプル)が過去すなわち前に発生されたワードまたはワ
ート群並びに現在及び過去の入力ワードまたはサンプル
の関数である点で巡回型である。このために、第4図の
フィルタは加算器(18)の出力側と加算器(18)の別な
入力側の間に接続された帰還ループを有し、加算器(18
は帰還ループを介して帰還された出力ワードと乗算器
(12)及び(16)からの各ワードを加算するように働
く。帰還ループはZ-1遅延要素(22)(遅延時間T)と
乗算器(24)を有し、乗算器(24)は遅延要素(22)か
らの出力信号を重み係数−b1を乗算する。(乗算器(1
2)及び(16)からの入力が正で、乗算器(24)からの
入力が負であるように考慮して、加算器が加算動作を行
うように構成されていれば、(−b1よりむしろ)値+b1
を代わりに乗算器(24)に供給することができる。) 第4図の回路が上述の式(1)に従って動作すること
は、第4図の回路の簡単な数学的分析によって容易に理
解できる。各クロック周期(T)の間、第4図のフィル
タは現在及び前の入力ワードと前の出力ワードを処理し
て新しい出力ワードを発生する。斯る各処理中、遅延要
素(22)により遅延され、乗算器(24)で重み係数−b1
を乗算された前の出力は、加算器(18)で乗算器(12)
及び(16)から到来するワードに加算される。乗算及び
加算動作は共に時間を要する。例えば乗算器(24)での
乗算動作に100nsを要し、加算器(18)での加算動作に4
0nsを要するものとする。この場合、フィルタが動作す
るには、クロック周期Tは140ns以下であってはならな
い。従って、フィルタはデータレート(入力ワードまた
はサンプル周波数)fsが1/140.10-9=7.14MHzより大き
な入力信号を処理できない。乗算器及び/又は加算器の
速度を増大することのみによって高速の信号を処理でき
るが、これは困難で及び/または非常に高価で及び/ま
たは不可能である。
ンプル)が過去すなわち前に発生されたワードまたはワ
ート群並びに現在及び過去の入力ワードまたはサンプル
の関数である点で巡回型である。このために、第4図の
フィルタは加算器(18)の出力側と加算器(18)の別な
入力側の間に接続された帰還ループを有し、加算器(18
は帰還ループを介して帰還された出力ワードと乗算器
(12)及び(16)からの各ワードを加算するように働
く。帰還ループはZ-1遅延要素(22)(遅延時間T)と
乗算器(24)を有し、乗算器(24)は遅延要素(22)か
らの出力信号を重み係数−b1を乗算する。(乗算器(1
2)及び(16)からの入力が正で、乗算器(24)からの
入力が負であるように考慮して、加算器が加算動作を行
うように構成されていれば、(−b1よりむしろ)値+b1
を代わりに乗算器(24)に供給することができる。) 第4図の回路が上述の式(1)に従って動作すること
は、第4図の回路の簡単な数学的分析によって容易に理
解できる。各クロック周期(T)の間、第4図のフィル
タは現在及び前の入力ワードと前の出力ワードを処理し
て新しい出力ワードを発生する。斯る各処理中、遅延要
素(22)により遅延され、乗算器(24)で重み係数−b1
を乗算された前の出力は、加算器(18)で乗算器(12)
及び(16)から到来するワードに加算される。乗算及び
加算動作は共に時間を要する。例えば乗算器(24)での
乗算動作に100nsを要し、加算器(18)での加算動作に4
0nsを要するものとする。この場合、フィルタが動作す
るには、クロック周期Tは140ns以下であってはならな
い。従って、フィルタはデータレート(入力ワードまた
はサンプル周波数)fsが1/140.10-9=7.14MHzより大き
な入力信号を処理できない。乗算器及び/又は加算器の
速度を増大することのみによって高速の信号を処理でき
るが、これは困難で及び/または非常に高価で及び/ま
たは不可能である。
最大動作速度に対して上述の制約をかなり緩和した第4
図のIIRフィルタの変形例を第1図に示す。しかし、第
1図の回路の説明する前に、その導出に至る数学的動作
を説明する。
図のIIRフィルタの変形例を第1図に示す。しかし、第
1図の回路の説明する前に、その導出に至る数学的動作
を説明する。
上記式(1)の分子と分母の両方が同じ量だけ乗算され
れば、伝達応答は変化しない。従って、分子及び分母が
共に分母の複数共役(極位置の共役)すなわち(1−b1
Z-1)を乗算されれば、式(1)は次の如くなる。
れば、伝達応答は変化しない。従って、分子及び分母が
共に分母の複数共役(極位置の共役)すなわち(1−b1
Z-1)を乗算されれば、式(1)は次の如くなる。
ここで、Z-2は2つのサンプリング間隔(2T)の遅延に
対応した遅延演算子であり、その他の符号は上述で定義
したとおりである。
対応した遅延演算子であり、その他の符号は上述で定義
したとおりである。
式(5)の伝送特性に対応した差分子方程式(2次)は
次のとおりである。
次のとおりである。
y(n)=a0x(n)+(a1−a0b1)x(n−1) −b1a1x(n−2)+b1 2y(n−2)……(6) 式(4)の分母のZ-1の関数である2つの項は相殺さ
れ、式(5)の分母には1とZ-2の関数の項とが残って
いることがわかる。式(1)に戻り、式(1)の分母は
数学的に第4図の帰還ループを表わしていることが第4
図の回路と式(1)の比較からわかる。従って全般に第
4図の回路と同じであるが、式(5)に従って動作する
IIRフィルタを実現することが可能であれば、分母はZ-2
の関数の項(しかしZ-1の関数の項でない)を含むの
で、各出力ワードを形成するのに利用できる時間を倍加
することは可能である。従って、同じ速度で動作する乗
算器及び加算器にとって、フィルタの最大動作速度は倍
加する。
れ、式(5)の分母には1とZ-2の関数の項とが残って
いることがわかる。式(1)に戻り、式(1)の分母は
数学的に第4図の帰還ループを表わしていることが第4
図の回路と式(1)の比較からわかる。従って全般に第
4図の回路と同じであるが、式(5)に従って動作する
IIRフィルタを実現することが可能であれば、分母はZ-2
の関数の項(しかしZ-1の関数の項でない)を含むの
で、各出力ワードを形成するのに利用できる時間を倍加
することは可能である。従って、同じ速度で動作する乗
算器及び加算器にとって、フィルタの最大動作速度は倍
加する。
このようなIIRフィルタを実現することが可能な一方法
を第1図に示す。第1図に示すIIRフィルタは入力端子
(30)に入力信号x(n)が供給される。上述の如く、
入力信号x(n)は時間的に間隔(サンプリング間隔)
Tだけ隔置された一連のディジタルワードから成る。第
1図の点Aに現われる入力信号を第2図(A)に示す。
入力信号は一連のワードX-2,X-1,X0,X1,X2等から成る。
を第1図に示す。第1図に示すIIRフィルタは入力端子
(30)に入力信号x(n)が供給される。上述の如く、
入力信号x(n)は時間的に間隔(サンプリング間隔)
Tだけ隔置された一連のディジタルワードから成る。第
1図の点Aに現われる入力信号を第2図(A)に示す。
入力信号は一連のワードX-2,X-1,X0,X1,X2等から成る。
入力信号はカスケードまたはタンデム接続された3つの
Z-1遅延要素(32),(34),(36)から成る遅延手段
を通され、Tに等しい時間の増加分だけ順次遅延され
る。これにより一連の4つの相対的に遅延した信号が発
生され各信号は前の信号に対してZ-1だけ連続して遅延
される。すなわち、第1図の点Aにおける入力信号は遅
延を受けず(第2図(A))、第1図の点Bにおける入
力信号はZ-1だけ遅延され(第2図(B))、第1図の
点Cにおける入力信号はZ-2だけ遅延され(第2図
(C))、第1図の点Dにおける入力信号はZ-3だけ遅
延される(第2図(D))。
Z-1遅延要素(32),(34),(36)から成る遅延手段
を通され、Tに等しい時間の増加分だけ順次遅延され
る。これにより一連の4つの相対的に遅延した信号が発
生され各信号は前の信号に対してZ-1だけ連続して遅延
される。すなわち、第1図の点Aにおける入力信号は遅
延を受けず(第2図(A))、第1図の点Bにおける入
力信号はZ-1だけ遅延され(第2図(B))、第1図の
点Cにおける入力信号はZ-2だけ遅延され(第2図
(C))、第1図の点Dにおける入力信号はZ-3だけ遅
延される(第2図(D))。
2つのIIR処理ユニット(38),(40)は接続されて点
A〜Dで得られる一連の4つの相対的に遅延した信号の
各グループを受ける。特に、ユニット(38)は点A〜C
の信号のグループを受け、ユニット(40)は点B〜Dの
信号のグループを受け、従って2つのグループは相互に
対して1ワードだけオフセットされている。
A〜Dで得られる一連の4つの相対的に遅延した信号の
各グループを受ける。特に、ユニット(38)は点A〜C
の信号のグループを受け、ユニット(40)は点B〜Dの
信号のグループを受け、従って2つのグループは相互に
対して1ワードだけオフセットされている。
2つのIIR処理ユニット(38),(40)は図に示すよう
に、実質的に同一の構成のものが好ましい。各ユニット
は3つのラッチ(42)を有し、各ラッチは入力信号の各
ワードのビット数に適切なビット容量を有し、クロック
信号CK/2(第2図)の制御のもとに、一連の4つの相対
的に遅延した信号の適当なグループ(A〜C又はB〜
D)からのワードを各ユニット(38),(40)内にラッ
チする。各ユニット(38),(40)内では、ラッチ(4
2)の各出力側が各乗算器(46)を介して加算器(44)
の各入力側に接続され、各乗算器(46)は重み係数a0,
(a1−b1a0)及び−b1a1を有し、これ等の重み係数は式
(5)の分子のものと対応している。
に、実質的に同一の構成のものが好ましい。各ユニット
は3つのラッチ(42)を有し、各ラッチは入力信号の各
ワードのビット数に適切なビット容量を有し、クロック
信号CK/2(第2図)の制御のもとに、一連の4つの相対
的に遅延した信号の適当なグループ(A〜C又はB〜
D)からのワードを各ユニット(38),(40)内にラッ
チする。各ユニット(38),(40)内では、ラッチ(4
2)の各出力側が各乗算器(46)を介して加算器(44)
の各入力側に接続され、各乗算器(46)は重み係数a0,
(a1−b1a0)及び−b1a1を有し、これ等の重み係数は式
(5)の分子のものと対応している。
第4図と同様に、各加算器(44)の出力側は遅延要素
(48)及び乗算器(50)から成る帰還ループを介して加
算器(44)の別な入力側に接続されている。しかし、こ
の場合、遅延要素(48)はZ-2遅延要素であり、従って2
Tの遅延が生ずることに注意されたい。また、乗算器(5
0)の重み係数は+b1 2であり、これは差分子方程式すな
わち式(6)の最後の項に現われていることに注意され
たい。
(48)及び乗算器(50)から成る帰還ループを介して加
算器(44)の別な入力側に接続されている。しかし、こ
の場合、遅延要素(48)はZ-2遅延要素であり、従って2
Tの遅延が生ずることに注意されたい。また、乗算器(5
0)の重み係数は+b1 2であり、これは差分子方程式すな
わち式(6)の最後の項に現われていることに注意され
たい。
加算器(44)の出力側は(クロック信号CK/2が供給され
る)各ラッチ(52)を介してマルチプレクサ(54)の各
入力側に接続され、マルチプレクサ(54)は出力端子
(56)を有し、ここにフィルタの濾波出力信号y(n)
が発生される。
る)各ラッチ(52)を介してマルチプレクサ(54)の各
入力側に接続され、マルチプレクサ(54)は出力端子
(56)を有し、ここにフィルタの濾波出力信号y(n)
が発生される。
第1図のIIRフィルタの動作を次に説明する。ラッチ(4
2),(52)及び処理ユニット(38),(40)の他の回
路に供給されるクロック信号CK/2は入力信号の周波数す
なわちデータレートの1/2の周波数である。すなわち、
クロック信号CK/2の周期は2Tに等しい。従って、遅延要
素(32),(34),(36)がZ-1の遅延を生じるため
に、それらはクロック信号CK/2の周波数の2倍の周波数
(fsに等しい)のクロック信号CKを供給される。遅延要
素(32),(34),(36)はラッチ(42),(52)と同
じ構成とすることができ、これ等は処理ユニット(3
8),(40)を通る種々の信号のパイプラインとして使
用され、両者の唯一の機能的差異はそれ等が異なった周
波数のクロック信号(夫々CK及びCK/2)を供給されるこ
とである。
2),(52)及び処理ユニット(38),(40)の他の回
路に供給されるクロック信号CK/2は入力信号の周波数す
なわちデータレートの1/2の周波数である。すなわち、
クロック信号CK/2の周期は2Tに等しい。従って、遅延要
素(32),(34),(36)がZ-1の遅延を生じるため
に、それらはクロック信号CK/2の周波数の2倍の周波数
(fsに等しい)のクロック信号CKを供給される。遅延要
素(32),(34),(36)はラッチ(42),(52)と同
じ構成とすることができ、これ等は処理ユニット(3
8),(40)を通る種々の信号のパイプラインとして使
用され、両者の唯一の機能的差異はそれ等が異なった周
波数のクロック信号(夫々CK及びCK/2)を供給されるこ
とである。
次に述べるように、ラッチ(42)は、遅延要素(32),
(34),(36)と関連して、入力信号を効果的にデマル
チプレクス(demuliplex)するように働き、クロック信
号CK/2の各クロック周期(2T)の間処理ユニット(3
8),(40)は、1ワードだけオフセットされており、
そして第4図の従来のIIRフィルタではTに等しい連続
するクロック周期の間処理(2倍の速度で)されなけれ
ばならなかった入力信号のワードの各組を同時にすなわ
ち並列に処理できる。その分(帰還ループ)にZ-1の項
のない式(5)を満たすことにより達成される処理ユニ
ット(38),(40)のデータレートを半減するこの効果
により、フィルタの最大動作速度の制約を半減すること
ができる。
(34),(36)と関連して、入力信号を効果的にデマル
チプレクス(demuliplex)するように働き、クロック信
号CK/2の各クロック周期(2T)の間処理ユニット(3
8),(40)は、1ワードだけオフセットされており、
そして第4図の従来のIIRフィルタではTに等しい連続
するクロック周期の間処理(2倍の速度で)されなけれ
ばならなかった入力信号のワードの各組を同時にすなわ
ち並列に処理できる。その分(帰還ループ)にZ-1の項
のない式(5)を満たすことにより達成される処理ユニ
ット(38),(40)のデータレートを半減するこの効果
により、フィルタの最大動作速度の制約を半減すること
ができる。
一連のクロックパルスを有するクロック信号CK/2を第2
図に示す。上述の如く、それは2Tに等しい周期を有し、
従って、入力信号のデータレートfsの1/2の周波数を有
する。第2図に矢印で示すように、クロック信号はその
クロックパルスの前縁でラッチ(42)(及びこれに接続
されたその他の回路)を動作状態とする。時間t0で生ず
るこのような縁の丁度1つ前で、ワードX0,X-1及びX-2
がユニット(38)のラッチ(42)に記憶され、ワードX
-1,X-2及びX-3がユニット(40)のラッチ(42)に記憶
されるものとする。時間t0で、これ等種々のワードはラ
ッチ(42)から開放され、ユニット(38),(40)で処
理される。従って、このときに、ワードX0が第1図の点
E(第2図(E)参照)に現われ、ワードX-1が第1図
の点G(第2図(G)参照)に現われる。
図に示す。上述の如く、それは2Tに等しい周期を有し、
従って、入力信号のデータレートfsの1/2の周波数を有
する。第2図に矢印で示すように、クロック信号はその
クロックパルスの前縁でラッチ(42)(及びこれに接続
されたその他の回路)を動作状態とする。時間t0で生ず
るこのような縁の丁度1つ前で、ワードX0,X-1及びX-2
がユニット(38)のラッチ(42)に記憶され、ワードX
-1,X-2及びX-3がユニット(40)のラッチ(42)に記憶
されるものとする。時間t0で、これ等種々のワードはラ
ッチ(42)から開放され、ユニット(38),(40)で処
理される。従って、このときに、ワードX0が第1図の点
E(第2図(E)参照)に現われ、ワードX-1が第1図
の点G(第2図(G)参照)に現われる。
第1図の説明から理解できるように、連続するクロック
周期(時点t0とt1の間)の間処理ユニット(38)は出力
ワードY0(第2図(F))を発生し、この出力ワードY0
は入力ワードX0,X-1及びX-2の組に基づいており、これ
は式(5)に適合する。同様に、同じクロック周期中、
処理ユニット(40)は出力ワードY-1(第2図(H))
を発生し、この出力ワードY-1は入力ワードX-1,X-2及び
X-3の組に基づいており、これはまた式(4)に適合す
る。出力ワードY-1及びY0は2Tに等しいクロック周期中
同時に(並列に)発生されるので、処理ユニットの回路
は2つのワードがTに等しい接続したクロック周期で連
続的に発生される第4図の従来の回路の場合よりその動
作を行うために2倍の時間を持っている。上述の如く、
帰還ループはZ-1遅延要素を組み込む必要はないが、こ
の実施例では単一のZ-2遅延要素(48)を備えることに
より、速度の軽減を達成している。(Z-2遅延要素(4
8)はラッチ(42),(52)と同じラッチで構成し、同
じクロック信号C/2を供給するようにしてもよい)。
周期(時点t0とt1の間)の間処理ユニット(38)は出力
ワードY0(第2図(F))を発生し、この出力ワードY0
は入力ワードX0,X-1及びX-2の組に基づいており、これ
は式(5)に適合する。同様に、同じクロック周期中、
処理ユニット(40)は出力ワードY-1(第2図(H))
を発生し、この出力ワードY-1は入力ワードX-1,X-2及び
X-3の組に基づいており、これはまた式(4)に適合す
る。出力ワードY-1及びY0は2Tに等しいクロック周期中
同時に(並列に)発生されるので、処理ユニットの回路
は2つのワードがTに等しい接続したクロック周期で連
続的に発生される第4図の従来の回路の場合よりその動
作を行うために2倍の時間を持っている。上述の如く、
帰還ループはZ-1遅延要素を組み込む必要はないが、こ
の実施例では単一のZ-2遅延要素(48)を備えることに
より、速度の軽減を達成している。(Z-2遅延要素(4
8)はラッチ(42),(52)と同じラッチで構成し、同
じクロック信号C/2を供給するようにしてもよい)。
処理ユニット(38),(40)により夫々発生された出力
ワードY0,Y-1はラッチ(52)を介してマルチプレクサ
(54)に供給され、ここでそれ等は適当な順位で前デー
タレートfsで共にマルチプレクスされて出力ワードシー
ケンスy(n)を発生する(第2図参照)。
ワードY0,Y-1はラッチ(52)を介してマルチプレクサ
(54)に供給され、ここでそれ等は適当な順位で前デー
タレートfsで共にマルチプレクスされて出力ワードシー
ケンスy(n)を発生する(第2図参照)。
第1図のフィルタの入力信号が実質的にデマルチプレッ
クスされるにも拘わらず、これは入力信号のどのワード
もスキッピング(Skipping)を含まないことに注意され
たい。これは受け入れることはできない。第1図のフィ
ルタは第4図のフィルタと同じ方法で入力信号を効果的
に濾波する。すなわちそれは同じ伝達応答を有するが、
ワードの組を直列に(全レートで)処理するよりむしろ
ワードのペアを並列に(減衰したレートで)処理するよ
うにする。
クスされるにも拘わらず、これは入力信号のどのワード
もスキッピング(Skipping)を含まないことに注意され
たい。これは受け入れることはできない。第1図のフィ
ルタは第4図のフィルタと同じ方法で入力信号を効果的
に濾波する。すなわちそれは同じ伝達応答を有するが、
ワードの組を直列に(全レートで)処理するよりむしろ
ワードのペアを並列に(減衰したレートで)処理するよ
うにする。
要するに、第1図のIIRフィルムは各出力ワードを発生
するのに利用できる時間を第4図の従来のIIRフィルム
より2倍もっていると云う利点を呈し、その結果関連す
る回路構成の同じ動作時間中に、それが処理できる最大
信号周波数すなわちデータレートは倍加する。
するのに利用できる時間を第4図の従来のIIRフィルム
より2倍もっていると云う利点を呈し、その結果関連す
る回路構成の同じ動作時間中に、それが処理できる最大
信号周波数すなわちデータレートは倍加する。
概して、フィルムの動作を決定する伝達特性を表わす式
の分母X(z)がZ-pの関数を作ると云う手法を用いる
ことにより、第1図のフィルタは入力データレートより
遅いP倍(この場合P=2)働く。同じ手法をPより大
きな値に適用できる(多分2の整数累乗に対応するとき
のみ)。例えば伝達特性方程式の分母X(z)はZ-4の
関数に減少されるかもしれないが、この場合データレー
トは4(P=4)の関数だけ減少され、4つの処理ユニ
ットを用いて4組の入力ワードを並列にすなわち、同時
に処理する。同様に、分母はZ-8の関数に減少されるか
もしれないが、この場合、データレートは8(P=8)
の関数に減少され、8つの処理ユニットを用いて8組の
入力ワードを並列に処理する。概して伝達特性の分母X
(z)は次の形式の式に減少され、 X(z)=1+(−1)P/2・b1 P・Z-P それによって、データレートはPの関数(Pは2の整数
倍数)だけ減少でき、P個の処理ユニットを用いてP組
の入力ワードを並列に処理する。しかしながら、明らか
に2より大きなPの値に対して、伝達特性方程式の分子
Y(z)は同じフィルタ応答を得るのにもっと複雑とな
る。
の分母X(z)がZ-pの関数を作ると云う手法を用いる
ことにより、第1図のフィルタは入力データレートより
遅いP倍(この場合P=2)働く。同じ手法をPより大
きな値に適用できる(多分2の整数累乗に対応するとき
のみ)。例えば伝達特性方程式の分母X(z)はZ-4の
関数に減少されるかもしれないが、この場合データレー
トは4(P=4)の関数だけ減少され、4つの処理ユニ
ットを用いて4組の入力ワードを並列にすなわち、同時
に処理する。同様に、分母はZ-8の関数に減少されるか
もしれないが、この場合、データレートは8(P=8)
の関数に減少され、8つの処理ユニットを用いて8組の
入力ワードを並列に処理する。概して伝達特性の分母X
(z)は次の形式の式に減少され、 X(z)=1+(−1)P/2・b1 P・Z-P それによって、データレートはPの関数(Pは2の整数
倍数)だけ減少でき、P個の処理ユニットを用いてP組
の入力ワードを並列に処理する。しかしながら、明らか
に2より大きなPの値に対して、伝達特性方程式の分子
Y(z)は同じフィルタ応答を得るのにもっと複雑とな
る。
この明細書の上述した部分はこの発明を1極IIRフィル
タに適用した場合を扱ったが、この発明はもっと複雑な
特性のIIRフィルタにも適用できる。次の一例は、いか
にこの発明を2局IIRフィルタに適用できるかと云うこ
とである。
タに適用した場合を扱ったが、この発明はもっと複雑な
特性のIIRフィルタにも適用できる。次の一例は、いか
にこの発明を2局IIRフィルタに適用できるかと云うこ
とである。
2極(2次)IIRフィルタの伝達特性は次式により表わ
すことができる。
すことができる。
ここでa2及びb2は定数(重み係数)を表わし、その他の
符号は上述に定義したとおりである。
符号は上述に定義したとおりである。
式(7)の伝達特性に対応する差分方程式(2次)は以
下の如くである。
下の如くである。
y(n)=a0x(n)+a1x(n−1)+a2x(n−2) −b1y(n−1)−b2y(n−2)……(8) 式(7)及び(8)がハードウエア信号で満たすことが
できる構成を一例を第5図に示す。この構成は大部分第
4図のものと同じであり、従ってここでは第4図のもの
と異なる点のみ説明する。第5図では、第4図の回路と
同じ部分は同じ参照符号で表わしている。
できる構成を一例を第5図に示す。この構成は大部分第
4図のものと同じであり、従ってここでは第4図のもの
と異なる点のみ説明する。第5図では、第4図の回路と
同じ部分は同じ参照符号で表わしている。
第5図の構成は次の(i),(ii)以外は第4図のもの
と同じである。
と同じである。
(i)更に図に示すように、Z-1遅延要素(60)及び乗
算器(62)(重み係数a2)が接続され、式(7)の分子
の第3項を発生する。及び (ii)帰還ループは更に図に示すように接続されたZ-1
遅延要素(64)及び乗算器(64)(重み係数−b2)を有
し、式(7)の分母の第3項を発生する。
算器(62)(重み係数a2)が接続され、式(7)の分子
の第3項を発生する。及び (ii)帰還ループは更に図に示すように接続されたZ-1
遅延要素(64)及び乗算器(64)(重み係数−b2)を有
し、式(7)の分母の第3項を発生する。
低いデータレートで動作するようにデマルチプレックス
できる第5図の2極フィルタの見解を得るためには、式
(7)の分子及び分母の各々と乗算するときに、Z-1の
関数の項を含まない分子を持つ伝達特性方程式(或いは
回路はこれを基礎のなすことができる)を生じる式を見
い出す必要がある。特にZ-2及びZ-4のみから成る分母の
式は適当である。このような式の1つは(1−b1Z-1+b
2Z-2)である。
できる第5図の2極フィルタの見解を得るためには、式
(7)の分子及び分母の各々と乗算するときに、Z-1の
関数の項を含まない分子を持つ伝達特性方程式(或いは
回路はこれを基礎のなすことができる)を生じる式を見
い出す必要がある。特にZ-2及びZ-4のみから成る分母の
式は適当である。このような式の1つは(1−b1Z-1+b
2Z-2)である。
式(7)の分子と分母に(1−b1Z-1+b2Z-2)を乗算す
ることにより次式が得られる。
ることにより次式が得られる。
Z-1及びZ-3の関数である式(9)の分母の項は相殺さ
れ、1とZ-4のみが残っていることがわかる。
れ、1とZ-4のみが残っていることがわかる。
上述の相殺を行い、定数を簡単化して式(9)を変形す
ることにより次式が得られる。
ることにより次式が得られる。
ただし w0=a0; w1=a1−a0b1; w2=a0b2−a1b1+a2; w3=a1b2−a2b1; w4=b2a2; v1=2b2−b1 2;及び v2=b2 2・ 式(9)の伝達特性に対応する差分子方程式(4次)は
次式の如くである。
次式の如くである。
y(n)=w0x(n)+w1x(n−1) +w2x(n−2)+wZx(n−3) +w4x(n−4)−v1y(n−2) −v2y(n−4)……(11) 上述からわかるように、式(10)は帰還ループがZ-2及
びZ-4(実際にはカスケート接続された2つのZ-2遅延要
素で実現)の関数のみを含むような方法でハードウエア
語で満たすことができ、第1図の1極IIRフィルタの完
成が1/2のデータレートで動作する第4図のデマルチプ
レックス化の見解に基づいていると云う同じ方法で、1/
2のデータレートで動作する第5図のデマルチプレック
ス化の見解に基づいて2極IIRフィルタを完成できる。
びZ-4(実際にはカスケート接続された2つのZ-2遅延要
素で実現)の関数のみを含むような方法でハードウエア
語で満たすことができ、第1図の1極IIRフィルタの完
成が1/2のデータレートで動作する第4図のデマルチプ
レックス化の見解に基づいていると云う同じ方法で、1/
2のデータレートで動作する第5図のデマルチプレック
ス化の見解に基づいて2極IIRフィルタを完成できる。
式(10)及び(11)を満たす構成の一例を第3図に示
す。第3図の構成は第1図のものと大部分同じであり、
従ってここでは第1図のものと異なる点のみを説明す
る。第3図において、第1図のフィルタの回路と同じ部
分は同一の参照符号で表わしている。
す。第3図の構成は第1図のものと大部分同じであり、
従ってここでは第1図のものと異なる点のみを説明す
る。第3図において、第1図のフィルタの回路と同じ部
分は同一の参照符号で表わしている。
第3図の構成は、IIR処理ユニット(38),(40)にお
いて、次の(i)〜(v)以外は第1図のものと同じで
ある。
いて、次の(i)〜(v)以外は第1図のものと同じで
ある。
(i)式(10)の分子はZ-3及びZ-4の関数である余分の
項をもっていることで式(7)の分子と異なると云うこ
とを考慮して、遅延手段は遅延要素(32),(34),
(36)とカスケードまたはタンデム接続されたZ-1遅延
要素(70),(72)を更に有する。
項をもっていることで式(7)の分子と異なると云うこ
とを考慮して、遅延手段は遅延要素(32),(34),
(36)とカスケードまたはタンデム接続されたZ-1遅延
要素(70),(72)を更に有する。
(ii)上側の3つの乗算器(46)に対する重み係数は、
式(5)の重み係数a0,a1及びa2よりむしろ式(10)の
重み係数w0,w1及びw2である。
式(5)の重み係数a0,a1及びa2よりむしろ式(10)の
重み係数w0,w1及びw2である。
(iii)更に2つのラッチ(24)及び乗算器(46)(重
み係数w3及びw4)が図に示すように接続される。
み係数w3及びw4)が図に示すように接続される。
(iv)帰還ループの乗算器(48)の重み係数は式(10)
の重み係数−V1である、及び (v)帰還ループはZ-2遅延要素(48)とカスケードま
たはタンデム接続されたZ-2遅延要素(74)及び図に示
すように接続された乗算器(76)(重み係数−v2)を更
に有し、式(10)の分母の第3項を発生する。
の重み係数−V1である、及び (v)帰還ループはZ-2遅延要素(48)とカスケードま
たはタンデム接続されたZ-2遅延要素(74)及び図に示
すように接続された乗算器(76)(重み係数−v2)を更
に有し、式(10)の分母の第3項を発生する。
第3図のフィルタは第1図のフィルタと非常に似た方法
で働き、第5図の従来の2極フィルタの応答を繰り返え
すが、斯る方法では2Tに等しい各クロック周期の間、ペ
アの組ワードがTに等しい各連続したクロック周期で連
続的に(全データレートで)処理される代りに、並列に
(半分のデータレートで)処理される。従って、第3図
のフィルタは第1図のフィルタが第4図にまさる利点を
呈したように第5図のフィルタにまさる利点を呈する。
すなわち、第3図のフィルタは出力信号ワードを発生す
るのに第5図のフィルタより2倍の時間を有、従って、
関連する回路構成の同じ動作時間中に処理できる最大入
力信号周波数すなわちデータレートは倍加する。
で働き、第5図の従来の2極フィルタの応答を繰り返え
すが、斯る方法では2Tに等しい各クロック周期の間、ペ
アの組ワードがTに等しい各連続したクロック周期で連
続的に(全データレートで)処理される代りに、並列に
(半分のデータレートで)処理される。従って、第3図
のフィルタは第1図のフィルタが第4図にまさる利点を
呈したように第5図のフィルタにまさる利点を呈する。
すなわち、第3図のフィルタは出力信号ワードを発生す
るのに第5図のフィルタより2倍の時間を有、従って、
関連する回路構成の同じ動作時間中に処理できる最大入
力信号周波数すなわちデータレートは倍加する。
第3図のフィルタの設計の基礎をなす原理は、全(入
力)データレートよりもっと遅いP倍(Pは2より大き
くかつ多分2の整数累乗である)で動作するフィルタを
作るのに拡大できる。例えばフィルタの動作を決定する
伝達特性方程式の分母X(z)はZ-4及びZ-8の関数に減
少するかもしれないが、この場合、データレートは4
(8,等)の関数だけ減少し、4(8,等)つの処理ユニッ
トを用いて4(8,等)つの入力ワードを並列に同時に処
理する。概して、分母X(z)は、Zに対する−P及び
−2Pの累乗の項を含むように、次数2PのZの関数を作る
ことができる。しかしながら、明らかに、2より大きな
Pの値のため、同じフィルタ応答を得るには式の分子は
非常に複雑になる。
力)データレートよりもっと遅いP倍(Pは2より大き
くかつ多分2の整数累乗である)で動作するフィルタを
作るのに拡大できる。例えばフィルタの動作を決定する
伝達特性方程式の分母X(z)はZ-4及びZ-8の関数に減
少するかもしれないが、この場合、データレートは4
(8,等)の関数だけ減少し、4(8,等)つの処理ユニッ
トを用いて4(8,等)つの入力ワードを並列に同時に処
理する。概して、分母X(z)は、Zに対する−P及び
−2Pの累乗の項を含むように、次数2PのZの関数を作る
ことができる。しかしながら、明らかに、2より大きな
Pの値のため、同じフィルタ応答を得るには式の分子は
非常に複雑になる。
第1図及び第3図の各実施例において、IIR処理ユニッ
ト(38),(40)の各々は、遅延手段(32),(34),
(36)または(32),(34),(36),(70),(72)
の各部分と関連して、式(5)または式(10)を満足す
る各デマルチプレクス化IIRフィルタまたはフィルタ部
を構成するが、2つのユニットで処理される隣接ワード
のグループで1ワードだけ相互に対してオフセットされ
ていることがわかる。例えば第1図の場合には、処理ユ
ニット(38)と遅延要素(32),(34)は一方のIIRフ
ィルタ部を形成し、処理ユニット(40)と遅延要素(3
2),(34),(36)は他方のIIRフィルタ部を形成す
る。第1図及び第3図のフィルタでは、遅延手段は各II
Rフィルタ部を分担してハードウエア的に節約してい
る。何故なら第1図では3つの遅延要素(32),(3
4),(36)を5つの遅延要素(ユニット(38)に2
つ、ユニット(40)に3つ)で置換しなければならず、
第3図では5つの遅延要素(32),(34),(36),
(70),(72)を9つの遅延要素(ユニット(38)で4
つ、ユニット(40)で5つ)で置換しなければならない
点で、各IIRフィルタ部の個々の遅延手段を組込むこと
は遅延要素の数を増大することになるからである。それ
でもなお、IIR処理ユニット(38),(40)及び分担し
た遅延手段は各個別の遅延手段を含む夫々のIIRフィル
タ部て置換できることがわかる。このような構成は、回
路基板の如き夫々実質的に同一の回路アセンブリィとし
て形成されるかもしれない2以上の実質的に同一の部分
のようにフィルタをハードウエア語で満たすことが可能
な場合には好ましいものと考えられる。
ト(38),(40)の各々は、遅延手段(32),(34),
(36)または(32),(34),(36),(70),(72)
の各部分と関連して、式(5)または式(10)を満足す
る各デマルチプレクス化IIRフィルタまたはフィルタ部
を構成するが、2つのユニットで処理される隣接ワード
のグループで1ワードだけ相互に対してオフセットされ
ていることがわかる。例えば第1図の場合には、処理ユ
ニット(38)と遅延要素(32),(34)は一方のIIRフ
ィルタ部を形成し、処理ユニット(40)と遅延要素(3
2),(34),(36)は他方のIIRフィルタ部を形成す
る。第1図及び第3図のフィルタでは、遅延手段は各II
Rフィルタ部を分担してハードウエア的に節約してい
る。何故なら第1図では3つの遅延要素(32),(3
4),(36)を5つの遅延要素(ユニット(38)に2
つ、ユニット(40)に3つ)で置換しなければならず、
第3図では5つの遅延要素(32),(34),(36),
(70),(72)を9つの遅延要素(ユニット(38)で4
つ、ユニット(40)で5つ)で置換しなければならない
点で、各IIRフィルタ部の個々の遅延手段を組込むこと
は遅延要素の数を増大することになるからである。それ
でもなお、IIR処理ユニット(38),(40)及び分担し
た遅延手段は各個別の遅延手段を含む夫々のIIRフィル
タ部て置換できることがわかる。このような構成は、回
路基板の如き夫々実質的に同一の回路アセンブリィとし
て形成されるかもしれない2以上の実質的に同一の部分
のようにフィルタをハードウエア語で満たすことが可能
な場合には好ましいものと考えられる。
この発明は第1図及び第3図と関連して上述した方法よ
り他の方法で実施できることがわかる。例えば、第1図
及び第3図に示した回路は置き換え(流れ図反転)手法
(エー・ヴィ・オッペンヘイム(A.V.Oppenheim)及び
アール・ダブリュ・シャフア(R.W.Schafer)著、1975
年プレンチス−ホール社発行、“ディジタル信号処理”
第130〜160頁参照)を用いて他の構成に置き換えること
ができる。置き換えにより得られた他の構成は第1図及
び第3図に示す構成と実質的に異なるけれども、それ等
は等価的に働き、第1図及び第3図の構成の同じ伝達特
性を有する。
り他の方法で実施できることがわかる。例えば、第1図
及び第3図に示した回路は置き換え(流れ図反転)手法
(エー・ヴィ・オッペンヘイム(A.V.Oppenheim)及び
アール・ダブリュ・シャフア(R.W.Schafer)著、1975
年プレンチス−ホール社発行、“ディジタル信号処理”
第130〜160頁参照)を用いて他の構成に置き換えること
ができる。置き換えにより得られた他の構成は第1図及
び第3図に示す構成と実質的に異なるけれども、それ等
は等価的に働き、第1図及び第3図の構成の同じ伝達特
性を有する。
更に、かつて伝達特性方程式用の適当な形の分母X
(z)を導出する数学的方法は完成したが、式を満たす
方法は上述したものと異なったルートを取ることができ
る。例えば第3図の2極デマルチプレクス化フィルタを
考える。式(10)を満足する第3図の構成は、式(7)
を満足する第5図の非デマルチプレクス化2極フィルタ
の構成に基づいている。しかし、当業者には周知の如
く、式(7)は第5図に示すものと異なった形で満たす
ことができる。従って、第3図のものと等価的に働くフ
ィルタを構成することは可能であり、それは式(10)を
満足するがむしろ異なった構成であり、第5図に示すも
のにより何らほかの式(7)を満足する構成に基づいて
いるからである。
(z)を導出する数学的方法は完成したが、式を満たす
方法は上述したものと異なったルートを取ることができ
る。例えば第3図の2極デマルチプレクス化フィルタを
考える。式(10)を満足する第3図の構成は、式(7)
を満足する第5図の非デマルチプレクス化2極フィルタ
の構成に基づいている。しかし、当業者には周知の如
く、式(7)は第5図に示すものと異なった形で満たす
ことができる。従って、第3図のものと等価的に働くフ
ィルタを構成することは可能であり、それは式(10)を
満足するがむしろ異なった構成であり、第5図に示すも
のにより何らほかの式(7)を満足する構成に基づいて
いるからである。
従って、概して、この発明で実施するフィルタは次の手
法によって設計できる。基本的または非デマルチプレク
ス化した形が選択され、基本的または非デマルチプレク
ス化した形の伝達特性方程式の重み係数に対する適当な
値を計算することによりフィルタの所望の応答が選択さ
れる。次に式(例えばX(z)の極位置の共役)は基本
的伝達特性方程式の分子Y(z)及び分母X(z)の両
方と乗算されるとき、基底方程式の変換を次式の分母X
(z)を持つ変形方程式に導くことがわかる。
法によって設計できる。基本的または非デマルチプレク
ス化した形が選択され、基本的または非デマルチプレク
ス化した形の伝達特性方程式の重み係数に対する適当な
値を計算することによりフィルタの所望の応答が選択さ
れる。次に式(例えばX(z)の極位置の共役)は基本
的伝達特性方程式の分子Y(z)及び分母X(z)の両
方と乗算されるとき、基底方程式の変換を次式の分母X
(z)を持つ変形方程式に導くことがわかる。
ここで、Pは2に等しくすなわち2の整数累乗であり、
フィルタがデマルチプレクスされる方法の数すなわち出
力信号のP個のワードを同時に処理するのに用いられる
IIR処理部の数を表わし、Lは基本的すなわち元の伝達
特性の極の数に等しい。変形方程式は次数m.P(ここで
mはフィルタの基本的または非デマルチプレクス化した
形の次数である)を有するデマルチプレクス化フィルタ
またはフィルタ部に対応する。次に変形方程式は任意の
所望の方法で実行され、IIRフィルタ部を形成するか、
またはその遅延手段をその(P−1)個の相対物と分担
するIIR処理ユニットを形成する。出来れば、P個の各
処理ユニットまたはフィルタ部に対して同じ実行が使用
され、それ等は全て実質的に同一の構成となる。しか
し、原理的には変形方程式は1以上の方法で実行され、
異なった構成の処理ユニットまたはフィルタ部を形成す
ることができる。次に乗算器はP個のフィルタ部(すな
わち共通の遅延手段を分担するP個の個別のフィルタ部
又はP個の処理ユニット)と組み合わされて、完全なデ
マルチプレクス化フィルタを形成する。
フィルタがデマルチプレクスされる方法の数すなわち出
力信号のP個のワードを同時に処理するのに用いられる
IIR処理部の数を表わし、Lは基本的すなわち元の伝達
特性の極の数に等しい。変形方程式は次数m.P(ここで
mはフィルタの基本的または非デマルチプレクス化した
形の次数である)を有するデマルチプレクス化フィルタ
またはフィルタ部に対応する。次に変形方程式は任意の
所望の方法で実行され、IIRフィルタ部を形成するか、
またはその遅延手段をその(P−1)個の相対物と分担
するIIR処理ユニットを形成する。出来れば、P個の各
処理ユニットまたはフィルタ部に対して同じ実行が使用
され、それ等は全て実質的に同一の構成となる。しか
し、原理的には変形方程式は1以上の方法で実行され、
異なった構成の処理ユニットまたはフィルタ部を形成す
ることができる。次に乗算器はP個のフィルタ部(すな
わち共通の遅延手段を分担するP個の個別のフィルタ部
又はP個の処理ユニット)と組み合わされて、完全なデ
マルチプレクス化フィルタを形成する。
この発明の実施例を添付図面を参照して詳細に説明した
けれども、この発明はこれ等の実施例に限定されること
なく、種々の変更、変形がこの発明の要旨を逸脱するこ
となくなし得ることが当業者には理解されるであろう。
けれども、この発明はこれ等の実施例に限定されること
なく、種々の変更、変形がこの発明の要旨を逸脱するこ
となくなし得ることが当業者には理解されるであろう。
上述の如くこの発明によれば、従来のIIRフィルタでは
処理できない高速入力信号を処理でき比較的低速用の部
品を用いて低廉化がはかれる。
処理できない高速入力信号を処理でき比較的低速用の部
品を用いて低廉化がはかれる。
第1図はこの発明で実施される1極IIRフィルタを示す
図、第2図は第1図のフィルタの種々の箇所に現われる
信号を示す図、第3図はこの発明で実施される2極IIR
フィルタを示す図、第4図は1極(1次)IIRフィルタ
の従来例を示す図、第5図は2極(2次)IIRフィルタ
の従来例を示す図である。 (32),(34),(36)は遅延要素、(38),(40)は
IIR処理ユニット、(54)はマルチプレクサである。
図、第2図は第1図のフィルタの種々の箇所に現われる
信号を示す図、第3図はこの発明で実施される2極IIR
フィルタを示す図、第4図は1極(1次)IIRフィルタ
の従来例を示す図、第5図は2極(2次)IIRフィルタ
の従来例を示す図である。 (32),(34),(36)は遅延要素、(38),(40)は
IIR処理ユニット、(54)はマルチプレクサである。
フロントページの続き (56)参考文献 特開 昭62−84611(JP,A) 特開 昭57−60773(JP,A) 特公 平3−10262(JP,B2) 特公 平1−20805(JP,B2) 米国特許4811263(US,A) 米国特許4328426(US,A) 米国特許477612(US,A) 英国特許2181008(GB,A) 欧州特許公開218396(EP,A) 欧州特許公開137464(EP,A) 欧州特許公開45596(EP,A)
Claims (1)
- 【請求項1】時間間隔Tだけ離れた一連のディジタルワ
ードを含む入力信号を濾波するための無限長インパルス
応答フィルタであって、 p・T(但しpは整数で且つ2の倍数)に等しい複数の
期間中の各々の期間に動作して入力信号の隣接ワードの
夫々のセットを処理し、それらのセットをお互いに1ワ
ードだけオフセットし、従って、前記各期間にp個の出
力ワードを作る、p個の巡回型無限長インパルス応答フ
ィルタ部分と、 前記p・Tに等しい各期間中に前記p個の出力ワードを
まとめてマルチプレクスし、時間的に前記期間Tだけ離
れた一連の出力ワードを含む濾波された出力信号を構成
する手段と、 を備えた無限長インパルス応答フィルタ。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| GB8523655 | 1985-09-25 | ||
| GB8523655A GB2181008B (en) | 1985-09-25 | 1985-09-25 | Infinite impulse response filters |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6272218A JPS6272218A (ja) | 1987-04-02 |
| JPH0740656B2 true JPH0740656B2 (ja) | 1995-05-01 |
Family
ID=10585697
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP61227028A Expired - Fee Related JPH0740656B2 (ja) | 1985-09-25 | 1986-09-25 | 無限長インパルス応答フイルタ |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US4811263A (ja) |
| EP (1) | EP0216580B1 (ja) |
| JP (1) | JPH0740656B2 (ja) |
| AT (1) | ATE65351T1 (ja) |
| CA (1) | CA1265588A (ja) |
| DE (1) | DE3680270D1 (ja) |
| GB (1) | GB2181008B (ja) |
Families Citing this family (31)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2181008B (en) | 1985-09-25 | 1989-09-20 | Sony Corp | Infinite impulse response filters |
| JPH01248816A (ja) * | 1988-03-30 | 1989-10-04 | Toshiba Corp | デジタルフィルタ |
| JP2870756B2 (ja) * | 1988-04-20 | 1999-03-17 | 株式会社リコー | 空間フィルタ画像処理装置 |
| JPH0828649B2 (ja) * | 1989-02-16 | 1996-03-21 | 日本電気株式会社 | ディジタルフィルタ |
| US4980849A (en) * | 1989-02-28 | 1990-12-25 | Raytheon Company | Method and apparatus for autoregressive model simulation |
| US5089981A (en) * | 1989-04-24 | 1992-02-18 | Audio Precision, Inc. | Hybrid form digital filter |
| US5049885A (en) * | 1989-10-13 | 1991-09-17 | Cincinnati Microwave, Inc. | Long range police radar warning receiver |
| US5134406A (en) * | 1989-10-13 | 1992-07-28 | Cincinnati Microwave, Inc. | Long range police radar warning receiver with multiple array averaging |
| US5079553A (en) * | 1989-10-13 | 1992-01-07 | Cincinnati Microwave, Inc. | Long range police radar warning receiver |
| US4954828A (en) * | 1989-10-13 | 1990-09-04 | Cincinnati Microwave, Inc. | Long range police radar warning receiver |
| JPH0537299A (ja) * | 1991-04-12 | 1993-02-12 | Sony Corp | 巡回形デイジタルフイルタ |
| US6108680A (en) * | 1991-10-30 | 2000-08-22 | Texas Instruments Incorporated | System and method for filtering using an interleaved/retimed architecture |
| JP3297880B2 (ja) * | 1992-03-18 | 2002-07-02 | テキサス インスツルメンツ インコーポレイテツド | Iirディジタル・フィルタ |
| US5432723A (en) * | 1993-12-02 | 1995-07-11 | Motorola, Inc. | Parallel infinite impulse response (IIR) filter with low quantization effects and method therefor |
| US5745396A (en) * | 1995-04-28 | 1998-04-28 | Lucent Technologies Inc. | Pipelined adaptive IIR filter |
| AUPO170196A0 (en) * | 1996-08-16 | 1996-09-12 | University Of Alberta | A finite-dimensional filter |
| US6058404A (en) * | 1997-04-11 | 2000-05-02 | Texas Instruments Incorporated | Apparatus and method for a class of IIR/FIR filters |
| US6446193B1 (en) | 1997-09-08 | 2002-09-03 | Agere Systems Guardian Corp. | Method and apparatus for single cycle processing of data associated with separate accumulators in a dual multiply-accumulate architecture |
| CA2269562C (en) * | 1998-04-29 | 2008-07-29 | Lockheed Martin Corporation | Drive-line damping control for an electric vehicle |
| US6529926B1 (en) * | 2000-09-20 | 2003-03-04 | Santel Networks, Inc. | Analog discrete-time FIR filter |
| US6553398B2 (en) | 2000-09-20 | 2003-04-22 | Santel Networks, Inc. | Analog fir filter with parallel interleaved architecture |
| US7042958B2 (en) * | 2003-06-04 | 2006-05-09 | Tropian, Inc. | Digital time alignment in a polar modulator |
| US8725785B1 (en) | 2004-08-09 | 2014-05-13 | L-3 Communications Corp. | Parallel infinite impulse response filter |
| ES2297964B1 (es) * | 2005-03-03 | 2009-08-25 | Angel Carpeño Velayos | Filtro incremental inteligente. |
| US7809927B2 (en) * | 2007-09-11 | 2010-10-05 | Texas Instruments Incorporated | Computation parallelization in software reconfigurable all digital phase lock loop |
| EP2559254B1 (en) * | 2010-04-13 | 2019-01-02 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Interpolation filter in a hybrid video decoder and encoder |
| US9053045B1 (en) | 2011-09-16 | 2015-06-09 | Altera Corporation | Computing floating-point polynomials in an integrated circuit device |
| US8949298B1 (en) | 2011-09-16 | 2015-02-03 | Altera Corporation | Computing floating-point polynomials in an integrated circuit device |
| US8543634B1 (en) * | 2012-03-30 | 2013-09-24 | Altera Corporation | Specialized processing block for programmable integrated circuit device |
| US9207909B1 (en) | 2012-11-26 | 2015-12-08 | Altera Corporation | Polynomial calculations optimized for programmable integrated circuit device structures |
| US9348795B1 (en) | 2013-07-03 | 2016-05-24 | Altera Corporation | Programmable device using fixed and configurable logic to implement floating-point rounding |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US477612A (en) | 1892-06-21 | Electric switch | ||
| US4328426A (en) | 1980-08-04 | 1982-05-04 | Xerox Corporation | Filter for image pixels |
| GB2181008A (en) | 1985-09-25 | 1987-04-08 | Sony Corp | Infinite impulse response filters |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2133118A5 (ja) * | 1971-04-08 | 1972-11-24 | Trt Telecom Radio Electr | |
| US4317092A (en) * | 1980-06-30 | 1982-02-23 | Hewlett-Packard Company | Recursive low pass digital filter |
| JPS5834615A (ja) * | 1981-08-24 | 1983-03-01 | Victor Co Of Japan Ltd | Iirデイジタルフイルタ |
| FR2546008B1 (fr) * | 1983-05-11 | 1985-07-12 | Labo Electronique Physique | Circuit d'egalisation adaptative et de demodulation conjointes |
| US4777612A (en) * | 1983-10-05 | 1988-10-11 | Nec Corporation | Digital signal processing apparatus having a digital filter |
-
1985
- 1985-09-25 GB GB8523655A patent/GB2181008B/en not_active Expired
-
1986
- 1986-09-11 AT AT86307025T patent/ATE65351T1/de not_active IP Right Cessation
- 1986-09-11 DE DE8686307025T patent/DE3680270D1/de not_active Expired - Lifetime
- 1986-09-11 EP EP86307025A patent/EP0216580B1/en not_active Expired - Lifetime
- 1986-09-23 CA CA000518802A patent/CA1265588A/en not_active Expired - Lifetime
- 1986-09-24 US US06/911,134 patent/US4811263A/en not_active Expired - Lifetime
- 1986-09-25 JP JP61227028A patent/JPH0740656B2/ja not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US477612A (en) | 1892-06-21 | Electric switch | ||
| US4328426A (en) | 1980-08-04 | 1982-05-04 | Xerox Corporation | Filter for image pixels |
| GB2181008A (en) | 1985-09-25 | 1987-04-08 | Sony Corp | Infinite impulse response filters |
| US4811263A (en) | 1985-09-25 | 1989-03-07 | Sony Corporation | Infinite impulse response filters |
Also Published As
| Publication number | Publication date |
|---|---|
| GB8523655D0 (en) | 1985-10-30 |
| CA1265588A (en) | 1990-02-06 |
| EP0216580A2 (en) | 1987-04-01 |
| GB2181008B (en) | 1989-09-20 |
| US4811263A (en) | 1989-03-07 |
| EP0216580B1 (en) | 1991-07-17 |
| DE3680270D1 (de) | 1991-08-22 |
| EP0216580A3 (en) | 1989-03-01 |
| ATE65351T1 (de) | 1991-08-15 |
| JPS6272218A (ja) | 1987-04-02 |
| GB2181008A (en) | 1987-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0740656B2 (ja) | 無限長インパルス応答フイルタ | |
| US9432043B2 (en) | Sample rate converter, an analog to digital converter including a sample rate converter and a method of converting a data stream from one data rate to another data rate | |
| US5367476A (en) | Finite impulse response digital filter | |
| JP3584027B2 (ja) | デジタルフィルタ | |
| JPH0126204B2 (ja) | ||
| JPH0834407B2 (ja) | 入力加重形トランスバーサルフィルタ | |
| JPH0421218B2 (ja) | ||
| EP0426294B1 (en) | Decimation-interpolation circuits | |
| JPH0998069A (ja) | Fir型ディジタルフィルタ | |
| JPH0590897A (ja) | オーバーサンプリングフイルタ回路 | |
| JP2865681B2 (ja) | フィルタの重み付け係数の調整方法 | |
| JPH10509011A (ja) | 改良されたディジタルフィルタ | |
| JPH1131945A (ja) | 非対称周波数応答特性を有する有限インパルス応答フィルタ | |
| WO2005002051A1 (ja) | デジタルフィルタ | |
| JP3097599B2 (ja) | ディジタルフィルタ | |
| KR100270165B1 (ko) | 디지털 인터폴레이션 기능이 있는 병렬 fir 필터 시스템 | |
| JPH0423609A (ja) | Firディジタルフィルタ | |
| JP3258938B2 (ja) | デシメーションフィルタ | |
| JPH0120805B2 (ja) | ||
| Noll | High throughput digital filters | |
| JPH04349710A (ja) | デジタル信号の出力ワード率を増加する補間回路 | |
| JPH02149011A (ja) | サンプリング周波数変換装置 | |
| RU1786638C (ru) | Цифровой нерекурсивный фильтр | |
| JPH05299973A (ja) | デシメーションフィルタ | |
| JP3041932B2 (ja) | サンプルレート変換回路 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |