JPWO2003023960A1 - デジタルフィルタおよびその設計方法 - Google Patents
デジタルフィルタおよびその設計方法 Download PDFInfo
- Publication number
- JPWO2003023960A1 JPWO2003023960A1 JP2003527888A JP2003527888A JPWO2003023960A1 JP WO2003023960 A1 JPWO2003023960 A1 JP WO2003023960A1 JP 2003527888 A JP2003527888 A JP 2003527888A JP 2003527888 A JP2003527888 A JP 2003527888A JP WO2003023960 A1 JPWO2003023960 A1 JP WO2003023960A1
- Authority
- JP
- Japan
- Prior art keywords
- sequence
- coefficient group
- filter coefficient
- filter
- value
- 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
Images
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/06—Non-recursive filters
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Picture Signal Circuits (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Compression Or Coding Systems Of Tv Signals (AREA)
- Television Systems (AREA)
- Transmission Systems Not Characterized By The Medium Used For Transmission (AREA)
Abstract
Description
本発明は、デジタルフィルタおよびその設計方法に関し、特に、複数の遅延器から成るタップ付き遅延線を備え、各タップの信号をそれぞれ数倍した後、加算して出力するFIRフィルタに用いて好適なものである。
背景技術
通信、計測、音声・画像信号処理、医療、地震学などの様々な分野で提供されている種々の電子機器においては、その内部で何らかのデジタル信号処理を行っているのが通常である。デジタル信号処理の最も重要な基本操作に、各種の信号や雑音が混在している入力信号の中から、必要なある周波数帯域の信号のみを取り出すフィルタリング処理がある。このために、デジタル信号処理を行う電子機器では、デジタルフィルタが用いられることが多い。
デジタルフィルタとしては、IIR(Infinite Impulse Response:無限長インパルス応答)フィルタやFIR(Finite Impulse Response:有限長インパルス応答)フィルタが多く用いられる。このうちFIRフィルタは、次のような利点を持つ。第1に、FIRフィルタの伝達関数の極はz平面の原点のみにあるため、回路は常に安定である。第2に、完全に正確な直線位相特性を実現することができる。
フィルタを通過域と阻止域との配置から分類すると、ローパスフィルタ、ハイパスフィルタ、帯域通過フィルタ、帯域消去フィルタの4つに分けられる。IIRフィルタやFIRフィルタで基本となるのはローパスフィルタであり、その他のハイパスフィルタ、帯域通過フィルタ、帯域消去フィルタは、ローパスフィルタから周波数変換等の処理を行うことによって導かれる。
電子機器の用途によっては、互いに対称的な周波数特性を有するローパスフィルタとハイパスフィルタとの対から成る対称型FIRフィルタを用いることがある。そのような対称型FIRフィルタは、まず基本となるローパスフィルタを設計し、これを周波数変換することにより、ローパスフィルタと特性が対称的なハイパスフィルタを設計する。
ところで、FIRフィルタは、有限時間長で表されるインパルス応答がそのままフィルタの係数となっている。したがって、FIRフィルタを設計するということは、希望の周波数特性が得られるようにフィルタ係数群を決定するということである。
従来、FIR型のローパスフィルタからハイパスフィルタ等のフィルタ係数を導く周波数変換では、フィルタのカットオフ周波数を変換する処理が行われていた。具体的には、サンプリング周波数とカットオフ周波数との比率をもとに、窓関数やチェビシェフ近似法などを用いた畳み込み演算等を行うことにより、フィルタの伝達関数を求め、それを更に周波数成分に置き換える処理を行っていた。
しかしながら、窓関数やチェビシェフ近似法などを用いた周波数変換は、その計算が非常に複雑である。そのため、これをソフトウェアで実現すると処理負荷が重くなり、ハードウェアで実現すると回路規模が大きくなるという問題があった。
また、従来の設計法で得られるフィルタの周波数特性は、窓関数や近似式に依存するので、これらをうまく設定しないと、良好な周波数特性を得ることができない。しかし、窓関数や近似式を適当に設定することは一般に困難であり、所望の周波数特性を有するフィルタを設計することは、非常に大変であった。
また、周波数変換によらず、ハイパスフィルタなどのフィルタ係数を直接的に求める方法も考えられる。しかし、この場合は、所望の周波数特性を得るのに必要なフィルタ係数を試行錯誤して求めなければならず、簡単には設計できないという問題があった。
本発明は、このような問題を解決するために成されたものであり、所望の周波数特性を有するFIRデジタルフィルタを簡易的に設計できるようにすることを目的とする。
発明の開示
本発明のデジタルフィルタは、複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、上記フィルタ係数群は、その数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定したものであることを特徴とする。
例えば、この条件を満たすフィルタ係数群の数列は、−1,0,9,16,9,0,−1の比率から成る。
このようなデジタルフィルタを複数段縦続接続しても良い。
本発明の他の態様では、複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、上記フィルタ係数群は、その数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定したものであることを特徴とする。
例えば、この条件を満たすフィルタ係数群の数列は、1,0,−9,16,−9,0,1の比率から成る。
このようなデジタルフィルタを複数段縦続接続しても良い。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1のフィルタ係数群を基本とし、上記第1のフィルタ係数群の各値と合算した結果が基準値となるように各値が設定された数列から成る第2のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2のフィルタ係数群を基本とし、上記第2のフィルタ係数群の各値と合算した結果が基準値となるように各値が設定された数列から成る第1のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とする。
また、本発明によるデジタルフィルタの設計方法は、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明の他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明のその他の態様では、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上述と同様の操作を行うことにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列を変換して、変換後の数列の合計値がゼロで、上記変換後の数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列を変換して、変換後の数列の合計値が非ゼロで、上記変換後の数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタをn段(n≧1)縦続接続し、これによって得られたn次のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記n次のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記n次のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタをn段(n≧1)縦続接続し、これによって得られたn次のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記n次のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記n次のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせるようにしたことを特徴とする。
本発明のその他の態様では、請求の範囲第30項、第33項、第38項、第39項の何れか1項に記載の設計方法により求められたフィルタ係数群を用いるフィルタをm段(m≧1)縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とする。
本発明のその他の態様では、請求の範囲第30項、第33項、第38項、第39項の何れか1項に記載の設計方法により求められたフィルタ係数群を用いるフィルタと、請求の範囲第30項、第33項、第38項、第39項の他の何れか1項に記載の設計方法により求められた上記フィルタ係数群を用いるフィルタとを任意に縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1の基本のフィルタ係数群を用いるように成された第1の基本フィルタをn段(n≧1)縦続接続し、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせ、これによってn次のフィルタ係数群を得るとともに、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2の基本のフィルタ係数群を用いるように成された第2の基本フィルタをm段(m≧1)縦続接続し、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせ、これによってm次のフィルタ係数群を得て、上記n次のフィルタ係数群を用いるフィルタと、上記m次のフィルタ係数群を用いるフィルタとを縦続接続し、さらに、当該縦続接続によって得られたフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した最初のフィルタ係数群、または、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した最初のフィルタ係数群をもとに、上記最初のフィルタ係数群を用いるように成されたフィルタを複数段縦続接続した場合に得られるフィルタ係数群であって、接続段数を1つずつ増やしていったときに変化していく周波数が目的の周波数を飛び超える段数だけ上記基本フィルタを縦続接続した場合に得られる2番目のフィルタ係数群を求め、上記求めた2番目のフィルタ係数群をもとに、当該2番目のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、転写フィルタ係数群を求めて、上記転写フィルタ係数群をもとに、上記転写フィルタ係数群を用いるように成された中間成果物のフィルタを複数段縦続接続した場合に得られるフィルタ係数群であって、接続段数を1つずつ増やしていったときに変化していく周波数が上記目的の周波数を飛び超える段数だけ上記中間成果物のフィルタを縦続接続した場合に得られる3番目のフィルタ係数群を更に求め、以下同様に、上記求めた3番目のフィルタ係数群をもとに、転写フィルタ係数群と、当該転写フィルタ係数群を用いるように成された中間成果物のフィルタを縦続接続した場合のフィルタ係数群とを繰り返し求めて、これによって得られる周波数と上記目的の周波数との誤差を小さくしていき、最終的なフィルタ係数群を得るようにしたことを特徴とする。
発明を実施するための最良の形態
以下、本発明の一実施形態を図面に基づいて説明する。
本実施形態のデジタルフィルタは、複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するトランスバーサル型のFIRフィルタである。FIRフィルタでは、これに与える複数のフィルタ係数の数列が対称型であれば、位相特性は直線となり、かつ、回路が常に安定となる。
本実施形態では、このFIRフィルタに対して与える対称型の基本的なフィルタ係数群として、図2の▲4▼に示すような数列{−1,0,9,16,9,0,−1}を用いる。この▲4▼に示すフィルタ係数群は、その数列が対称型であり、数列の合計値が非ゼロで、数列の1つ飛びの合計値が同符号で互いに等しくなるという性質を持っている(−1+9+9+(−1)=16、0+16+0=16)。
上記図2の▲4▼に示すフィルタ係数は、同図の▲3▼に示すデジタル基本関数を1回移動平均演算することによって得られるものである。このデジタル基本関数は、1クロック毎にデータ値が−1,1,8,8,1,−1と変化する関数である。このデジタル基本関数の数列は、同図の▲1▼に示す数列を2回移動平均演算することによって得られる。
ここで、上記デジタル基本関数の数列の意味について、図19を用いて説明する。
図19は、上記デジタル基本関数を1回移動平均して得られる上述の数値列{−1,0,9,16,9,0,−1}に対して、4倍のオーバーサンプリングとコンボリューション演算とを施した結果を示す図である。
図19(a)において、一番左の列に示される一連の数値列は、元の数値列{−1,0,9,16,9,0,−1}に対して4倍のオーバーサンプリングを行った値である。また、一番左から右に向かって4列分の数値列は、一番左の列に示される数値列を1つずつ下方向にシフトしていったものである。図19(a)の列方向は時間軸を示しており、数値列を下方向にシフトするということは、一番左の列に示される数値列を徐々に遅延させていくことに対応する。
すなわち、左から2列目の数値列は、一番左の列に示される数値列を4倍周波数のクロック4CKの1/4位相分だけずらした数値列であることを示す。また、左から3列目の数値列は、左から2列目に示される数値列を4倍周波数のクロック4CKの1/4位相分だけずらした数値列、左から4列目の数値列は、左から3列目に示される数値列を4倍周波数のクロック4CKの1/4位相分だけ更にずらした数値列であることを示す。
また、左から5列目の数値列は、1〜4列目の各数値列を対応する行どうしで加算して4で割った値である。この左から5列目までの処理により、4相のコンボリューション演算を伴う4倍のオーバーサンプリングがデジタル的に実行されることになる。
上記5列目から右に向かって4列分の数値列は、5列目に示される数値列を1つずつ下方向にシフトしていったものである。また、左から9列目の数値列は、5〜8列目の各数値列を対応する行どうしで加算して4で割った値である。この左から9列目までの処理により、4相のコンボリューション演算を伴う4倍のオーバーサンプリングがデジタル的に2回実行されることになる。
また、左から10列目の数値列は、9列目に示される数値列を1つ下方向にシフトしたものである。また、左から11列目(一番右の列)の数値列は、9列目の数値列と10列目の数値列とを対応する行どうしで加算して2で割った値である。
この図19(a)の一番右の列に示される最終的に得られた数値列をグラフ化したのが、図19(b)である。この図19(b)に示すような波形を有する関数は、横軸に沿った標本位置がt1からt4の間にあるときにのみ“0”以外の有限な値を有し、それ以外の領域では値が全て“0”となる関数、つまり標本位置t1,t4において値が“0”に収束する関数である。このように関数の値が局所的な領域で“0”以外の有限の値を有し、それ以外の領域で“0”となる場合を「有限台」と称する。
また、この図19(b)に示す関数は、中央の標本位置t5においてのみ極大値をとり、t1,t2,t3,t4の4つの標本位置において値が“0”になるという特徴を有する標本化関数であり、滑らかな波形のデータを得るために必要なサンプル点は全て通る。
図3は、図2の▲1▼〜▲4▼に示す数列をFFT(Fast Fourier Transfer:高速フーリエ変換)した結果の周波数−ゲイン特性を示す図である。なお、この図3では、ゲインを“1”で基準化している。この特性図から分かるように、▲4▼の数列をフィルタ係数に用いた場合には、中心周波数においてゲインが0.5となり、かつ、低周波領域でのオーバーシュートや高周波領域でのリンギングも存在しない良好なローパスフィルタ特性が得られる。
このようなローパスフィルタ特性を実現する上述した{−1,0,9,16,9,0,−1}の数値列は、図19(b)に示した有限台の標本化関数の基礎となるものである。従来一般的に用いられていた標本化関数はt=±∞の標本位置で“0”に収束するのに対し、図19(b)に示した標本化関数は、t=t1,t4の有限の標本位置で“0”に収束する。
そのため、上述の数値列をFFT変換した場合、t=t1〜t4の範囲内に相当するデータだけが意味を持つ。t=t1〜t4の範囲外に相当するデータについては、本来これを考慮すべきであるのに無視しているという訳ではなく、理論的に考慮する必要がないため、打ち切り誤差は発生しない。したがって、上述の数値列をフィルタ係数として用いれば、窓関数を用いて係数の打ち切りを行う必要もなく、良好なローパスフィルタ特性を得ることができる。よって、これを基本ローパスフィルタとする。
図1は、上記基本ローパスフィルタの構成例を示す図である。このフィルタでは、縦続接続された6個のD型フリップフロップ11〜16によって入力信号を1クロックCKずつ順次遅延させる。そして、D型フリップフロップ11〜16の入出力タップから取り出した信号に対し、5個の係数器21〜25により上記▲4▼の数列から成るフィルタ係数をそれぞれ乗算し、それらの乗算結果を全て加算器31〜34で加算して出力する。
加算器34の出力段に設けられている1/32倍の乗算器35は、各タップの信号と各フィルタ係数とを積和演算することによって大きくなった振幅を元の振幅に戻すために用いる。つまり、ゲインを“1”で基準化している。“32”という値は、全てのフィルタ係数値を加算した結果である(−1+0+9+16+9+0+(−1)=32)。ここでは、図2の▲4▼に示す整数列を基本のフィルタ係数群として用い、加算出力を1/32倍しているが、フィルタ係数群の値自体を1/32倍しても良い。
図4は、図1に示す基本ローパスフィルタをn段(n≧1)縦続接続した場合に、それによって得られるn次のフィルタ係数群をそれぞれ示す図である。ここでは代表として、1次、2次、4次のフィルタ係数群を示している。なお、ここに示したフィルタ係数群は、係数値自体を1/32倍して10−3で丸めたものである。
上述のように、1次のローパスフィルタ係数群は、7個のフィルタ係数の数列{−1,0,9,16,9,0,−1}/32から成る。これに対し、基本ローパスフィルタを2段縦続接続した場合の2次のローパスフィルタ係数群は13個のフィルタ係数、4段縦続接続した場合の4次のローパスフィルタ係数群は25個のフィルタ係数、…というように、接続する段数に応じてフィルタ係数の数が増えていく。
図4に示した1次、2次、4次の各フィルタ係数群のうち、太枠で示した部分は、それぞれのフィルタ係数群の中央値を示す。図4において、2次および4次のフィルタ係数群におけるR[i]列C[j]行のフィルタ係数値は、以下の式により求められる。
(−(R[i−8]C[j−1]+R[i−2]C[j−1])+(R[i−6]C[j−1]+R[i−4]C[j−1])*9+R[i−5]C[j−1]*16)/32
図5は、n次(n=1,2,4,8,16)ローパスフィルタのフィルタ係数群の数列をFFTした結果の周波数−ゲイン特性を示す図である。図5の例では、44.1KHzのサンプリング周波数で量子化した信号を、同じく44.1KHzの遅延レードを持つ(D型フリップフロップ11〜16に供給するクロックCKの周波数が44.1KHz)ローパスフィルタに通した場合のフィルタ特性を示している。
図5に示すように、基本ローパスフィルタをn段縦続接続して次数を大きくしていくことにより、通過域から阻止域へと移る中間の遮断域(過渡域)の傾斜およびカットオフ周波数を段階的に変化させることができる。すなわち、基本ローパスフィルタを任意の段数だけ縦続接続することにより、遮断域の傾斜およびカットオフ周波数として所望の特性を有するローパスフィルタを簡単に設計することができる。なお、図5から分かるように、遮断域の傾斜およびカットオフ周波数の変化幅は、次数を大きくするに従って徐々に小さくなっていく。
図6は、1次〜16次のフィルタ係数群の最大値(中央値)と、ゲインが平坦な部分より−3db下がるカットオフ周波数のバンド幅との相関関係を示す図である。これから分かるように、両者間には相関があり、係数最大値およびカットオフ周波数の双方とも、次数を大きくするほど小さくなる。カットオフ周波数が段階的に小さくなる様子は、図5にも表されている。カットオフ周波数を微細に調整することも可能であるが、その方法については後述する。
本実施形態では、基本のフィルタ係数群として、{−1,0,9,16,9,0,−1}の比率から成る数列を用いている。これにより、基本ローパスフィルタを何段接続しても、通過域が平坦な特性を持ち、オーバーシュートが全く存在しない良好なローパスフィルタ特性を得ることができる。本実施形態で用いる基本のフィルタ係数群がこのような特性を持つことにより、当該基本ローパスフィルタをもとに、任意の周波数特性を持つ各種のデジタルフィルタを簡易的に設計することが可能となる(これについては後述する)。
また、本実施形態では、D型フリップフロップ11〜16の遅延レート(クロックレート)を変えることによって、ローパスフィルタの周波数帯域をシフトさせることができる。例えば、44.1KHzのサンプリングクロックで音声信号を量子化した場合に、ローパスフィルタの通過域を0〜22KHzとして、量子化信号から23〜44.1KHzの周波数成分を除去するためには、D型フリップフロップ11〜16のクロック周波数を88.2KHz以上に大きくする必要がある(クロック周波数が44.1KHzだと、図5に示したように通過域を0〜22KHzにできない)。
逆に、D型フリップフロップ11〜16のクロックレートをサンプリング周波数の1/nに小さくすることにより、異なる周波数特性を持つローパスフィルタを得ることができる。図7は、クロックレートを1/4にした場合のローパスフィルタの構成例を示す図であり、図1に示した構成要素と同一の機能を有する構成要素には同一の符号を付している。なお、ここで用いるクロックレートが1/4のクロックは、図7中に示したものに限らず、デューティ比が50%のものを用いても良い。
図8は、図7に示したローパスフィルタおよびこれを4段縦続接続した4次のローパスフィルタの周波数−ゲイン特性を示す図である。この図8と図5との比較から分かるように、クロックレートを小さくすると、通過域のバンド幅が減少し、遮断域の傾斜も急になる。これにより、クロックレートを変えるだけで、遮断域の傾斜およびカットオフ周波数として所望の特性を有するローパスフィルタを簡単に設計することができる。また、このような特性を利用して、任意の周波数特性を持つ各種のデジタルフィルタを簡易的に設計することも可能となる(これについても後述する)。
さらに、上述のような基本ローパスフィルタをもとに、ハイパスフィルタを簡単に導くこともできる。ここで求めるハイパスフィルタは、その周波数−ゲイン特性が基本ローパスフィルタと比べて、ある周波数軸に対して左右に対称な特性を持つ対称型のハイパスフィルタと、基本ローパスフィルタとハイパスフィルタとのゲイン出力を合算すると、全周波数帯域でゲインが完全に1(基準値)となる補完型のハイパスフィルタとが考えられる。
図1に示す基本ローパスフィルタに対して、ゲインが1に基準化された単一パルスを入力すると、各係数器21〜25を通して基本のフィルタ係数群{−1,0,9,16,9,0,−1}が出力される。これは、このフィルタ係数群の数列そのものがローパスフィルタの特性を表していると考えられる。したがって、基本ローパスフィルタのフィルタ係数群の数列を変えることにより、ハイパスフィルタを作ることができる。
例えば、基本ローパスフィルタとの対称性を保つために、数列の絶対値はそのままにして中央値以外を符号変換する。基本のフィルタ係数群に対してこのような条件を満たすフィルタ係数群の数列は、{1,0,−9,16,−9,0,1}となる。このようにすると、サンプリング周波数の1/4,3/4の周波数軸に対して左右に対称な特性を持つハイパスフィルタを得ることができる。
図9は、基本ローパスフィルタを4段縦続接続した4次のクロックレート1(基準のクロックCK)のローパスフィルタおよび、基本のフィルタ係数群を符号変換して生成した転写フィルタ係数群を用いるハイパスフィルタを4段縦続接続した4次のクロックレート1の対称型ハイパスフィルタの周波数−ゲイン特性を示す図である。
この特性図から分かるように、4次のローパスフィルタと比べて、サンプリング周波数(44.1KHz)の1/4,3/4の周波数軸に対して左右に対称な特性を持つ4次のハイパスフィルタを得ることができる。なお、4次ローパスフィルタのフィルタ係数群を求めた後、その数列を転写フィルタ係数群に変換した場合でも、1次ローパスフィルタのフィルタ係数群を転写フィルタ係数群に変換した後、当該転写フィルタ係数群から4次ハイパスフィルタのフィルタ係数群を求めた場合でも、その数列は全く同じものとなる。
また、他の例では、基本ローパスフィルタのフィルタ係数群の数列{−1,0,9,16,9,0,−1}と対応する要素どうしで加算した場合に、中央値どうしの加算結果が32(これは数列の各値を合算した基準値であり、各フィルタ係数の値が1/32倍されているときは、基準値は1である)で、それ以外が0となる数列を新しいフィルタ係数群とする。このようにすると、得られるフィルタ係数群はその数列の合計値が0、数列の1つ飛びの合計値が逆符号で等しい値となり、ローパスフィルタのゲインとのトータルが全周波数帯域で完全に1となるという特性を持つハイパスフィルタを得ることができる。
図10は、4次のクロックレート1のローパスフィルタに用いるフィルタ係数群から、4次のクロックレート1のハイパスフィルタに用いるフィルタ係数群を求める転写について説明するための図である。なお、ここでは、フィルタ係数群の値自体を1/32倍し、各フィルタ係数値を10−3で丸めた例について示しており、丸めた結果として値が0となる部分の図示は省略している。
ここでは、まず、4次のローパスフィルタのフィルタ係数群を小数点以下3桁で取り出す(表の第1列)。次に、当該取り出したフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転する。数列の中央値については、基準値(この場合は1)から中央値を減算する(表の第2列)。これによって得られた転写フィルタ係数群(表の第3列)が、求める4次のハイパスフィルタのフィルタ係数群である。
図11は、基本ローパスフィルタを4段縦続接続した4次のクロックレート1のローパスフィルタおよび、基本のフィルタ係数群を符号変換等して生成した転写フィルタ係数群を用いるハイパスフィルタを4段縦続接続した4次のクロックレート1の補完型ハイパスフィルタの周波数−ゲイン特性を示す図である。
この特性図から分かるように、4次ローパスフィルタと比べて、当該4次ローパスフィルタのゲインとのトータルが全周波数帯域で完全に1になるという特性を持つ4次ハイパスフィルタを得ることができる。なお、4次ローパスフィルタのフィルタ係数群を求めた後、その数列を転写フィルタ係数群に変換した場合でも、1次ローパスフィルタのフィルタ係数群を転写フィルタ係数群に変換した後、当該転写フィルタ係数群から4次ハイパスフィルタのフィルタ係数群を求めた場合でも、その数列は全く同じものとなる。
以上のように、本実施形態によれば、FIRローパスフィルタに与える基本的なフィルタ係数群として、所定の条件を満たす数列を与えるようにしたので、ローパスフィルタに対するフィルタ係数群の符号を変える程度の簡単な操作だけで、ハイパスフィルタを設計することができる。
なお、以上では、{−1,0,9,16,9,0,−1}の比率から成る数列をフィルタ係数群として利用するローパスフィルタを基本として、当該フィルタ係数群を転写することにより、ハイパスフィルタのフィルタ係数群を求める例について説明した。これとは逆に、{1,0,−9,16,−9,0,1}の比率から成る数列をフィルタ係数群として利用するハイパスフィルタを基本として、当該フィルタ係数群を転写することにより、ローパスフィルタのフィルタ係数群を求めることも可能である。
次に、以上に説明したフィルタの縦続接続、クロックレートの変換、フィルタ係数群の転写を駆使して、所望の周波数特性を有するFIRデジタルフィルタを設計する例について説明する。ここでは一例として、サンプリング周波数44.1KHz、−3db周波数(カットオフ周波数)が20KHzのローパスフィルタを設計する場合について説明する。
まず、基本ローパスフィルタをもとに、n次のクロックレート1/2(クロックCK/2)のローパスフィルタLPF2(n)を作る。このとき、図12に示すように、次数nを調整して、フィルタの通過域を必要なバンド幅にする。図12の例では、基本ローパスフィルタを8段縦続接続して8次のローパスフィルタを作っている。
これとは別に、基本ハイパスフィルタ(基本ローパスフィルタのフィルタ係数群を転写したもの)をもとに、n次のクロックレート1のハイパスフィルタHPF1(n)を作る。このとき、図13の実線に示すように、次数nを調整して、フィルタの通過域のバンド幅が十分に小さくなるようにする。図13の例では、基本ハイパスフィルタを21段縦続接続して21次のハイパスフィルタを作っている。なお、図13中には、比較のために、図12に示した8次のローパスフィルタの特性も併せて示している(点線)。
次に、上述のように作成した8次のローパスフィルタLPF2(8)と、21次のハイパスフィルタHPF1(21)とを縦続接続することにより、それぞれの通過域のオーバーラップした部分を取り出す。図14は、その場合の周波数−ゲイン特性、周波数−位相特性を示す図である。これから分かるように、通過域のバンド幅が十分に小さく、遮断域の傾斜が十分に大きい周波数特性を持ったハイパスフィルタを得ることができる。また、ほぼ直線的な位相特性も得ることができる。
なお、ここでは8次のローパスフィルタと21次のハイパスフィルタとの合成により通過域のバンド幅を小さくしているが、ハイパスフィルタの次数をより大きくすることにより、通過城のバンド幅を小さくしていくことも可能である。しかし、上述したように、遮断域の傾斜およびカットオフ周波数の変化幅は、次数を大きくするに従って徐々に小さくなっていく。したがって、次数を大きくするだけでは、通過域のバンド幅を小さくするのに限界がある。また、あまり次数を大きくし過ぎると、デジタルフィルタの段数が増えて、回路規模が増大する。これに対して、上述のようにローパスフィルタとハイパスフィルタとを合成する方法によれば、次数をそれほど大きくしなくても、通過域のバンド幅を十分に小さくすることができる。
次に、合成により得たハイパスフィルタの係数を10−3で丸めた後、それを補完型の転写係数に変換する。その結果を示したのが図15である。図15において、左側の数列は、合成フィルタ係数を10−3で丸めた結果を示すものであり、右側の数列は、補完型の転写係数に変換した結果を示すものである。このように、ハイパスフィルタのフィルタ係数群を転写フィルタ係数群に変換することにより、ローパスフィルタのフィルタ係数群を得ることができる。そこで、図15の右側に示す転写フィルタ係数の数列を新たなフィルタ係数群として、求めるローパスフィルタを構成する。さらに、当該ローパスフィルタを複数段縦続接続して次数を調整し、通過域のバンド幅を調整する。
図16および図17は、このようにして作成したローパスフィルタの周波数−ゲイン特性および周波数−位相特性を示す図である。これから分かるように、カットオフ周波数がほぼ20KHz、通過域がほぼ平坦で、遮断域の傾斜が極めて大きい理想型に近いローパスフィルタを得ることができる。また、直線位相特性も実現することができる。
以上のように、本実施形態によれば、{−1,0,9,16,9,0,−1}の比率から成る数列をフィルタ係数群として用いる基本ローパスフィルタを利用して、フィルタの縦続接続、クロックレートの変換、フィルタ係数群の転写を任意に組み合わせて適用することにより、所望の周波数特性を有するデジタルフィルタを極めて簡単に設計することができる。
次に、フィルタの縦続接続およびフィルタ係数群の転写を行うことにより、カットオフ周波数を微細に調整する方法について説明する。
上述したように、ローパスフィルタの次数を大きくしていくと、カットオフ周波数は小さくなる方へ変化していく。逆に、ハイパスフィルタの次数を大きくしていくと、カットオフ周波数は大きくなる方へ変化していく。しかも、これらの変化幅は、次数を大きくするに従って徐々に小さくなっていく。
このような性質を利用して、図18に示すように、まず1次のローパスフィルタの接続段数を1つずつ増やしていく。このとき、縦続接続によって小さくなる方へ変化していくカットオフ周波数が、目的とするカットオフ周波数を飛び超えるまで、接続段数を増やしていく(a)。そして、ローパスフィルタをn段接続して目的のカットオフ周波数を飛び超えたら、当該縦続接続によって得られたフィルタ係数群を、補完型の転写係数に変換する。これにより、n次のハイパスフィルタのフィルタ係数群が得られる。
次に、このようにして得たn次のハイパスフィルタの接続段数を1つずつ増やしていく。ここでは、n次のハイパスフィルタを縦続接続していくので、カットオフ周波数の変化幅は、1次のローパスフィルタをn段接続したときの変化幅に続けて更に小さくなっていく。そして、このとき大きい方へ変化していくカットオフ周波数が、目的とするカットオフ周波数を先ほどとは逆方向から再び飛び超えるまで、接続段数を増やしていく(b)。
そして、n次のハイパスフィルタをm段接続して目的のカットオフ周波数を飛び超えたら、当該縦続接続によって得られたフィルタ係数群を、補完型の転写係数に変換する。これにより、nm次のローパスフィルタのフィルタ係数群が得られる。次に、このようにして得たnm次のローパスフィルタの接続段数を1つずつ増やしていく。そして、小さい方へ変化していくカットオフ周波数が、目的とするカットオフ周波数を再び飛び超えるまで、接続段数を増やしていく(c)。
以下同様にして、当該縦続接続によって得られたフィルタ係数群をもとに、その転写フィルタ係数群の生成と、当該転写フィルタ係数群を用いる中間成果物のフィルタの縦続接続とを繰り返し行う。これにより、得られるカットオフ周波数と目的のカットオフ周波数との誤差を小さくしていき、最終的なフィルタ係数群を得る。
このような周波数の微調整方法を用いれば、任意のカットオフ周波数を特性として持つローパスフィルタあるいはハイパスフィルタを、極めて簡単な操作で設計することができる。
以上に説明した本実施形態のデジタルフィルタおよびフィルタ設計の手法は、ハードウェア構成、DSP、ソフトウェアの何れによっても実現することが可能である。例えばソフトウェアによって実現する場合、本実施形態のデジタルフィルタおよびフィルタ設計装置は、実際にはコンピュータのCPUあるいはMPU、RAM、ROMなどで構成され、RAMやROMに記憶されたプログラムが動作することによって実現できる。
したがって、コンピュータが上記本実施形態の機能を果たすように動作させるプログラムを例えばCD−ROMのような記録媒体に記録し、コンピュータに読み込ませることによって実現できるものである。上記プログラムを記録する記録媒体としては、CD−ROM以外に、フレキシブルディスク、ハードディスク、磁気テープ、光ディスク、光磁気ディスク、DVD、不揮発性メモリカード等を用いることができる。また、上記プログラムをインターネット等のネットワークを介してコンピュータにダウンロードすることによっても実現できる。
また、コンピュータが供給されたプログラムを実行することにより上述の実施形態の機能が実現されるだけでなく、そのプログラムがコンピュータにおいて稼働しているOS(オペレーティングシステム)あるいは他のアプリケーションソフト等と共同して上述の実施形態の機能が実現される場合や、供給されたプログラムの処理の全てあるいは一部がコンピュータの機能拡張ボードや機能拡張ユニットにより行われて上述の実施形態の機能が実現される場合も、かかるプログラムは本発明の実施形態に含まれる。
なお、上記実施形態では、基本ローパスフィルタに対するフィルタ係数群の数列として{−1,0,9,16,9,0,−1}を用い、基本ハイパスフィルタに対するフィルタ係数群の数列として{1,0,−9,16,−9,0,1}を用いたが、上記実施形態で説明した条件を満たす数列であれば、これ以外の数列をフィルタ係数群として用いても良い。
また、上記実施形態では、補完型の転写フィルタ係数群の例として、転写前後の各フィルタ係数値を対応する要素どうしで加算した場合、数列の中央値どうしの加算結果が1となり、それ以外の加算結果が0となる場合について説明したが、係数値の加算結果がトータルで1となるようにしても良い。例えば、中央値どうしの加算結果が0.5で、その両隣の値どうしの加算結果が0.25ずつである場合などが該当する。
このような条件を満たす転写フィルタ係数群を求める方法は、基本とするフィルタ係数群の数列の中央値については0.5から中央値を減算し、中央値の両隣の値については0.25からその両隣の値を減算する。そして、これ以外については絶対値をそのままにして符号反転すればよい。
また、上述した実施形態は単なる一例を示したに過ぎないものであって、本発明はこれに限定されるものではない。すなわち、フィルタの縦続接続、クロックレートの変換、フィルタ係数群の転写を任意に組み合わせて適用することにより、任意の周波数特性を有するFIRデジタルフィルタを簡易的に設計することが可能である。
その他、上記実施形態は、何れも本発明を実施するにあたっての具体化の一例を示したものに過ぎず、これによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその精神、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
以上説明したように本発明によれば、所望の周波数特性を有するFIRデジタルフィルタを極めて簡易的に設計することができる。
産業上の利用可能性
本発明は、所望の周波数特性を有するFIRデジタルフィルタを簡易的に設計できるようにするのに有用である。
【図面の簡単な説明】
図1は、基本ローパスフィルタの構成例を示す図である。
図2は、基本のフィルタ係数群の説明図である。
図3は、図2に示したフィルタ係数群の周波数−ゲイン特性を示す図である。
図4は、図1の基本ローパスフィルタをn段(n≧1)縦続接続した場合に得られるn次のフィルタ係数群の数列を示す図である。
図5は、n次(n=1,2,4,8,16)のフィルタ係数群の数列をFFTした結果の周波数−ゲイン特性を示す図である。
図6は、1次〜16次のフィルタ係数群の最大値とカットオフ周波数のバンド幅との相関関係を示す図である。
図7は、クロックレートを1/4にした場合のローパスフィルタの構成例を示す図である。
図8は、図7に示したローパスフィルタおよびこれを4段縦続接続した4次ローパスフィルタの周波数−ゲイン特性を示す図である。
図9は、4次のクロックレート1のローパスフィルタおよび4次のクロックレート1の対称型ハイパスフィルタの周波数−ゲイン特性を示す図である。
図10は、4次のクロックレート1のローパスフィルタ係数群から、4次のクロックレート1のハイパスフィルタ係数群を求める転写を説明するための図である。
図11は、4次のクロックレート1のローパスフィルタおよび4次のクロックレート1の補完型ハイパスフィルタの周波数−ゲイン特性を示す図である。
図12は、8次のクロックレート1/2のローパスフィルタの周波数−ゲイン特性を示す図である。
図13は、21次のクロックレート1のハイパスフィルタの周波数−ゲイン特性を示す図である。
図14は、図12に示す8次のローパスフィルタと、図13に示す21次のハイパスフィルタとを合成した場合に得られる周波数−ゲイン特性および周波数−位相特性を示す図である。
図15は、8次ローパスフィルタと21次ハイパスフィルタとの合成フィルタ係数から補完型の転写係数への転写を説明するための図である。
図16は、最終成果物であるローパスフィルタの周波数−ゲイン特性を示す図である。
図17は、最終成果物であるローパスフィルタの周波数−位相特性を示す図である。
図18は、カットオフ周波数を微細に調整するための方法を説明するための図である。
図19は、本実施形態で用いるデジタル基本関数の意味を説明するための図である。
は一般に困難であり、所望の周波数特性を有するフィルタを設計することは、非常に大変であった。
また、周波数変換によらず、ハイパスフィルタなどのフィルタ係数を直接的に求める方法も考えられる。しかし、この場合は、所望の周波数特性を得るのに必要なフィルタ係数を試行錯誤して求めなければならず、簡単には設計できないという問題があった。
本発明は、このような問題を解決するために成されたものであり、所望の周波数特性を有するFIRデジタルフィルタを簡易的に設計できるようにすることを目的とする。
発明の開示
本発明のデジタルフィルタは、
複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1のフィルタ係数群を基本とし、上記第1のフィルタ係数群の各値と合算した結果が基準値となるように各値が設定された数列から成る第2のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2のフィルタ係数群を基本とし、上記第2のフィルタ係数群の各値と合算した結果が基準値となるように各値が設定された数列から成る第1のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とする。
また、本発明によるデジタルフィルタの設計方法は、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるよ
うに値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明の他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明のその他の態様では、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上述と同様の操作を行うことにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタをn段(n≧1)縦続接続し、これによって得られたn次のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記n次のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタをn段(n≧1)縦続接続し、これによって得られたn次のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記n次のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする。
本発明のその他の態様では、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせるようにしたことを特徴とする。
本発明のその他の態様では、請求の範囲第31項、第32項、第34項、第35項の何れか1項に記載の設計方法により求められたフィルタ係数群を用いるフィルタをm段(m≧1)縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とする。
本発明のその他の態様では、複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1の基本のフィルタ係数群を用いるように成された第1の基本フィルタをn段(n≧1)縦続接続し、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせ、これによってn次のフィルタ係数群を得るとともに、数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2の基本のフィルタ係数群を用いるように成された第2の基本フィルタをm段(m≧1)縦続接続し、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせ、これによってm次のフィルタ係数群を得て、上記n次のフィルタ係数群を用いるフィルタと、上記m次のフィルタ係数群を用いるフィルタとを縦続接続し、さらに、当該縦続接続によって得られたフィルタ係数群の数列の
値列を下方向にシフトするということは、一番左の列に示される数値列を徐々に遅延させていくことに対応する。
すなわち、左から2列目の数値列は、一番左の列に示される数値列を4倍周波数のクロック4CKの1/4位相分だけずらした数値列であることを示す。また、左から3列目の数値列は、左から2列目に示される数値列を4倍周波数のクロック4CKの1/4位相分だけずらした数値列、左から4列目の数値列は、左から3列目に示される数値列を4倍周波数のクロック4CKの1/4位相分だけ更にずらした数値列であることを示す。
また、左から5列目の数値列は、1〜4列目の各数値列を対応する行どうしで加算した値である。この左から5列目までの処理により、4相のコンボリューション演算を伴う4倍のオーバーサンプリングがデジタル的に実行されることになる。
上記5列目から右に向かって4列分の数値列は、5列目に示される数値列を1つずつ下方向にシフトしていったものである。また、左から9列目の数値列は、5〜8列目の各数値列を対応する行どうしで加算した値である。この左から9列目までの処理により、4相のコンボリューション演算を伴う4倍のオーバーサンプリングがデジタル的に2回実行されることになる。
また、左から10列目の数値列は、9列目に示される数値列を1つ下方向にシフトしたものである。また、左から11列目(一番右の列)の数値列は、9列目の数値列と10列目の数値列とを対応する行どうしで加算した値である。
この図19(a)の一番右の列に示される最終的に得られた数値列をグラフ化したのが、図19(b)である。この図19(b)に示すような波形を有する関数は、横軸に沿った標本位置がt1からt4の間にあ
Claims (46)
- 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
上記フィルタ係数群は、その数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定したものであることを特徴とするデジタルフィルタ。 - 上記フィルタ係数群の数列が−1,0,9,16,9,0,−1の比率から成ることを特徴とする請求の範囲第1項に記載のデジタルフィルタ。
- 請求の範囲第2項に記載のデジタルフィルタを複数段縦続接続したことを特徴とするデジタルフィルタ。
- 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
上記フィルタ係数群は、その数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定したものであることを特徴とするデジタルフィルタ。 - 上記フィルタ係数群の数列が1,0,−9,16,−9,0,1の比率から成ることを特徴とする請求の範囲第4項に記載のデジタルフィルタ。
- 請求の範囲第5項に記載のデジタルフィルタを複数段縦続接続したことを特徴とするデジタルフィルタ。
- 請求の範囲第2項に記載のデジタルフィルタと請求の範囲第5項に記載のデジタルフィルタとを縦続接続したことを特徴とするデジタルフィルタ。
- 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1のフィルタ係数群を基本とし、上記第1のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転することによって得た第2のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1のフィルタ係数群を基本とし、上記第1のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することによって得た第2のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 上記第1のフィルタ係数群は、−1,0,9,16,9,0,−1の比率から成る数列をフィルタ係数群として用いるように成された基本フィルタを1段もしくは複数段縦続接続した場合に得られる数列から成ることを特徴とする請求の範囲第8項に記載のデジタルフィルタ。
- 上記第1のフィルタ係数群は、−1,0,9,16,9,0,−1の比率から成る数列をフィルタ係数群として用いるように成された基本フィルタを1段もしくは複数段縦続接続した場合に得られる数列から成ることを特徴とする請求の範囲第9項に記載のデジタルフィルタ。
- 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1のフィルタ係数群を基本とし、上記第1のフィルタ係数群の各値と合算した結果が基準値となるように各値が設定された数列から成る第2のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2のフィルタ係数群を基本とし、上記第2のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転することによって得た第1のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2のフィルタ係数群を基本とし、上記第2のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することによって得た第1のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 上記第2のフィルタ係数群は、1,0,−9,16,−9,0,1の比率から成る数列をフィルタ係数群として用いるように成された基本フィルタを1段もしくは複数段縦続接続した場合に得られる数列から成ることを特徴とする請求の範囲第13項に記載のデジタルフィルタ。
- 上記第2のフィルタ係数群は、1,0,−9,16,−9,0,1の比率から成る数列をフィルタ係数群として用いるように成された基本フィルタを1段もしくは複数段縦続接続した場合に得られる数列から成ることを特徴とする請求の範囲第14項に記載のデジタルフィルタ。
- 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2のフィルタ係数群を基本とし、上記第2のフィルタ係数群の各値と合算した結果が基準値となるように各値が設定された数列から成る第1のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1のフィルタ係数群を基本とし、上記第1のフィルタ係数群の数列を変換して、変換後の数列の合計値がゼロで、上記変換後の数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 複数の遅延器から成るタップ付き遅延線を備え、各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタであって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2のフィルタ係数群を基本とし、上記第2のフィルタ係数群の数列を変換して、変換後の数列の合計値が非ゼロで、上記変換後の数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として与えるように設計されたことを特徴とするデジタルフィルタ。 - 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 上記基本のフィルタ係数群は、−1,0,9,16,9,0,−1の比率の数列から成ることを特徴とする請求の範囲第20項に記載のデジタルフィルタの設計方法。
- 上記基本のフィルタ係数群は、−1,0,9,16,9,0,−1の比率の数列から成ることを特徴とする請求の範囲第21項に記載のデジタルフィルタの設計方法。
- 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 上記基本のフィルタ係数群は、1,0,−9,16,−9,0,1の比率の数列から成ることを特徴とする請求の範囲第24項に記載のデジタルフィルタの設計方法。
- 上記基本のフィルタ係数群は、1,0,−9,16,−9,0,1の比率の数列から成ることを特徴とする請求の範囲第25項に記載のデジタルフィルタの設計方法。
- 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列を変換して、変換後の数列の合計値がゼロで、上記変換後の数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群をもとに、上記基本のフィルタ係数群の数列を変換して、変換後の数列の合計値が非ゼロで、上記変換後の数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定することにより、上記各タップの信号に対するフィルタ係数群を求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタをn段(n≧1)縦続接続し、これによって得られたn次のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 上記n次のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする請求の範囲第30項に記載のデジタルフィルタの設計方法。
- 上記n次のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする請求の範囲第30項に記載のデジタルフィルタの設計方法。
- 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタをn段(n≧1)縦続接続し、これによって得られたn次のフィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 上記n次のフィルタ係数群の数列の絶対値をそのままにして上記数列の中央値以外を符号反転し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする請求の範囲第33項に記載のデジタルフィルタの設計方法。
- 上記n次のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算し、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とする請求の範囲第33項に記載のデジタルフィルタの設計方法。
- 上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせるようにしたことを特徴とする請求の範囲第30項に記載のデジタルフィルタの設計方法。
- 上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせるようにしたことを特徴とする請求の範囲第33項に記載のデジタルフィルタの設計方法。
- 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタの、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせるようにしたことを特徴とするデジタルフィルタの設計方法。 - 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した基本のフィルタ係数群を用いるように成された基本フィルタの、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせるようにしたことを特徴とするデジタルフィルタの設計方法。 - 請求の範囲第30項に記載の設計方法により求められたフィルタ係数群を用いるフィルタをm段(m≧1)縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とするデジタルフィルタの設計方法。
- 請求の範囲第33項に記載の設計方法により求められたフィルタ係数群を用いるフィルタをm段(m≧1)縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とするデジタルフィルタの設計方法。
- 請求の範囲第38項に記載の設計方法により求められたフィルタ係数群を用いるフィルタをm段(m≧1)縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とするデジタルフィルタの設計方法。
- 請求の範囲第39項に記載の設計方法により求められたフィルタ係数群を用いるフィルタをm段(m≧1)縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とするデジタルフィルタの設計方法。
- 請求の範囲第30項、第33項、第38項、第39項に記載の設計方法により求められたフィルタ係数群を用いるフィルタと、請求の範囲第30項、第33項、第38項、第39項の他の何れか1項に記載の設計方法により求められた上記フィルタ係数群を用いるフィルタとを任意に縦続接続することにより、フィルタの通過周波数帯域を調整するようにしたことを特徴とするデジタルフィルタの設計方法。
- 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した第1の基本のフィルタ係数群を用いるように成された第1の基本フィルタをn段(n≧1)縦続接続し、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせ、これによってn次のフィルタ係数群を得るとともに、
数列が対称型であり、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した第2の基本のフィルタ係数群を用いるように成された第2の基本フィルタをm段(m≧1)縦続接続し、上記遅延器の遅延レートを変えることによってフィルタの通過周波数帯域をシフトさせ、これによってm次のフィルタ係数群を得て、
上記n次のフィルタ係数群を用いるフィルタと、上記m次のフィルタ係数群を用いるフィルタとを縦続接続し、
さらに、当該縦続接続によって得られたフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、これによって得られた転写フィルタ係数群を、上記各タップの信号に対するフィルタ係数群として求めるようにしたことを特徴とするデジタルフィルタの設計方法。 - 複数の遅延器から成るタップ付き遅延線における各タップの信号を、与えられるフィルタ係数群によりそれぞれ数倍した後、加算して出力するデジタルフィルタを設計する方法であって、
数列が対称型であり、上記数列の合計値が非ゼロで、上記数列の1つ飛びの合計値が同符号で互いに等しくなるように値を設定した最初のフィルタ係数群、または、上記数列の合計値がゼロで、上記数列の1つ飛びの合計値が逆符号で互いに等しくなるように値を設定した最初のフィルタ係数群をもとに、上記最初のフィルタ係数群を用いるように成されたフィルタを複数段縦続接続した場合に得られるフィルタ係数群であって、接続段数を1つずつ増やしていったときに変化していく周波数が目的の周波数を飛び超える段数だけ上記基本フィルタを縦続接続した場合に得られる2番目のフィルタ係数群を求め、
上記求めた2番目のフィルタ係数群をもとに、当該2番目のフィルタ係数群の数列の中央値以外については絶対値をそのままにして符号反転するとともに、上記数列の中央値については基準値から上記中央値を減算することにより、転写フィルタ係数群を求めて、
上記転写フィルタ係数群をもとに、上記転写フィルタ係数群を用いるように成された中間成果物のフィルタを複数段縦続接続した場合に得られるフィルタ係数群であって、接続段数を1つずつ増やしていったときに変化していく周波数が上記目的の周波数を飛び超える段数だけ上記中間成果物のフィルタを縦続接続した場合に得られる3番目のフィルタ係数群を更に求め、
以下同様に、上記求めた3番目のフィルタ係数群をもとに、転写フィルタ係数群と、当該転写フィルタ係数群を用いるように成された中間成果物のフィルタを縦続接続した場合のフィルタ係数群とを繰り返し求めて、これによって得られる周波数と上記目的の周波数との誤差を小さくしていき、最終的なフィルタ係数群を得るようにしたことを特徴とするデジタルフィルタの設計方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001273345 | 2001-09-10 | ||
JP2001273345 | 2001-09-10 | ||
PCT/JP2002/009046 WO2003023960A1 (fr) | 2001-09-10 | 2002-09-05 | Filtre numerique et son procede de conception |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2003023960A1 true JPWO2003023960A1 (ja) | 2004-12-24 |
JP4300272B2 JP4300272B2 (ja) | 2009-07-22 |
Family
ID=19098574
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003527888A Expired - Fee Related JP4300272B2 (ja) | 2001-09-10 | 2002-09-05 | デジタルフィルタおよびその設計方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US7330865B2 (ja) |
EP (1) | EP1441440A1 (ja) |
JP (1) | JP4300272B2 (ja) |
KR (1) | KR20040029471A (ja) |
CN (1) | CN100358239C (ja) |
TW (1) | TW569523B (ja) |
WO (1) | WO2003023960A1 (ja) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030200243A1 (en) * | 2002-04-19 | 2003-10-23 | Matsushita Electric Industrial Co., Ltd. | Finite impulse response filter, communication transmission apparatus and communication reception apparatus |
JPWO2004079905A1 (ja) * | 2003-03-03 | 2006-06-08 | 有限会社ニューロソリューション | デジタルフィルタの設計方法および装置、デジタルフィルタ設計用プログラム、デジタルフィルタ |
TW200501565A (en) | 2003-05-15 | 2005-01-01 | Neuro Solution Corp | Digital filter and design method, design apparatus, and digital filter design program thereof |
CN1938947A (zh) * | 2004-02-17 | 2007-03-28 | 神经网路处理有限公司 | 数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器 |
JP2005287820A (ja) * | 2004-03-31 | 2005-10-20 | U-Medica Inc | 生体現象測定記録装置及びノイズ成分除去方法 |
JPWO2006048958A1 (ja) * | 2004-11-05 | 2008-05-22 | 有限会社ニューロソリューション | デジタルフィルタおよびその設計方法、設計装置、デジタルフィルタ設計用プログラム |
JP4982080B2 (ja) * | 2005-12-16 | 2012-07-25 | パナソニック株式会社 | デジタルフィルタ |
KR100871690B1 (ko) * | 2006-07-12 | 2008-12-05 | 삼성전자주식회사 | 디지털 제어가 가능한 필터 시스템 및 필터링 방법 |
IL178744A0 (en) * | 2006-10-19 | 2007-09-20 | Eci Telecom Ltd | Method for estimating bandwidth limiting effects in transmission communication systems |
CN101299600B (zh) * | 2007-04-30 | 2010-08-25 | 大唐移动通信设备有限公司 | 一种基于ofdm系统的数字滤波器的设计方法 |
US8949303B2 (en) * | 2008-06-10 | 2015-02-03 | Japanese Science And Technology Agency | Filter |
DE102009038594A1 (de) | 2009-08-26 | 2011-03-03 | Sew-Eurodrive Gmbh & Co. Kg | Verfahren zur Verarbeitung eines analog bestimmten Messwertsignals, Resolveranordnung zur Durchführung des Verfahrens und Verfahren zur Bestimmung eines Ausgangsstromes eines Umrichters |
CN102388534B (zh) * | 2010-02-20 | 2014-04-30 | 华为技术有限公司 | 滤波设备和用于提供滤波设备的方法 |
CN105591630B (zh) * | 2014-10-20 | 2018-06-29 | 普天信息技术有限公司 | 基于fpga的插值滤波器优化的方法及装置 |
JP6356584B2 (ja) * | 2014-11-28 | 2018-07-11 | 株式会社デンソー | 予測進路推定装置 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4542475A (en) * | 1982-11-23 | 1985-09-17 | Rca Corporation | Sampled data filter system as for a digital TV |
JP2544095B2 (ja) * | 1984-08-17 | 1996-10-16 | ヤマハ株式会社 | 電子楽器 |
JPH05243908A (ja) * | 1991-12-06 | 1993-09-21 | Nec Corp | ろ波器 |
GB9226536D0 (en) * | 1992-12-21 | 1993-02-17 | Unilever Plc | Foodstuffs and other compositions |
WO1997009780A1 (en) * | 1995-09-07 | 1997-03-13 | Philips Electronics N.V. | Improved digital filter |
JP4023842B2 (ja) * | 1995-09-28 | 2007-12-19 | ソニー株式会社 | ディジタルフィルタ及び音響再生装置 |
CA2177664C (en) * | 1996-05-29 | 2003-01-28 | Andrew G. Deczky | Digital receive filter for communications systems |
JPH10149348A (ja) * | 1996-07-01 | 1998-06-02 | Sun Microsyst Inc | フィルタ対称を使用した信号内挿及びデシメーション |
US6662200B2 (en) * | 2001-01-03 | 2003-12-09 | Intel Corporation | Multiplierless pyramid filter |
-
2002
- 2002-09-05 CN CNB028176138A patent/CN100358239C/zh not_active Expired - Fee Related
- 2002-09-05 WO PCT/JP2002/009046 patent/WO2003023960A1/ja active Application Filing
- 2002-09-05 KR KR10-2004-7003450A patent/KR20040029471A/ko not_active Application Discontinuation
- 2002-09-05 JP JP2003527888A patent/JP4300272B2/ja not_active Expired - Fee Related
- 2002-09-05 EP EP02798034A patent/EP1441440A1/en not_active Withdrawn
- 2002-09-09 TW TW091120479A patent/TW569523B/zh not_active IP Right Cessation
-
2004
- 2004-03-01 US US10/708,405 patent/US7330865B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN100358239C (zh) | 2007-12-26 |
WO2003023960A1 (fr) | 2003-03-20 |
KR20040029471A (ko) | 2004-04-06 |
US7330865B2 (en) | 2008-02-12 |
TW569523B (en) | 2004-01-01 |
CN1554147A (zh) | 2004-12-08 |
JP4300272B2 (ja) | 2009-07-22 |
US20040172434A1 (en) | 2004-09-02 |
EP1441440A1 (en) | 2004-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4300272B2 (ja) | デジタルフィルタおよびその設計方法 | |
US7529788B2 (en) | Digital filter design method and device, digital filter design program, and digital filter | |
JPWO2006048958A1 (ja) | デジタルフィルタおよびその設計方法、設計装置、デジタルフィルタ設計用プログラム | |
US20050289206A1 (en) | Digital filter design method and device, digital filter design program, digital filter | |
KR20050075365A (ko) | 디지털 필터의 설계 방법 및 장치, 디지털 필터 설계용프로그램, 디지털 필터 | |
JPWO2005078925A1 (ja) | デジタルフィルタの設計方法および装置、デジタルフィルタ設計用プログラム、デジタルフィルタ | |
US7590675B2 (en) | Digital filter, design method thereof, design device, and digital filter design program | |
KR20050071604A (ko) | 디지털 필터의 설계 방법 및 장치, 디지털 필터 설계용프로그램, 디지털 필터 | |
EP1533898A1 (en) | Digital filter designing method, digital filter designing program, digital filter | |
WO2003047097A1 (fr) | Procede de conception d'un filtre numerique, appareil de projet, programme de projet pour la realisation d'un tel filtre, et filtre numerique y relatif | |
WO2007102611A1 (ja) | 補間関数生成回路 | |
JP4265119B2 (ja) | デジタルフィルタ | |
JP4300273B2 (ja) | 音質調整装置及びこれに用いるフィルタ装置、音質調整方法、フィルタの設計方法 | |
US20050171988A1 (en) | Digital filter design method and device, digital filter design program, and digital filter | |
JP2005020554A (ja) | デジタルフィルタ | |
US20050120067A1 (en) | Digital filter designing method, digital filter designing program, digital filter | |
Mottaghi-Kashtiban et al. | FIR filters involving shifts and only two additions, efficient for short word-length signal processing | |
Kostić et al. | Sensitivity of notch filters realized by parallel connection of allpass subfilters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20040621 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050822 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050822 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081021 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20081118 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090303 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20090319 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090319 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120501 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |