JPH01162409A - フィルタ係数演算装置 - Google Patents

Abstract

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

Description

【発明の詳細な説明】 産業上の利用分野 本発明は、任意の振幅周波数特性及び位相周波数特性を
実現するトランスバーサル・フィルタＣ以下、ＦＩＲフ
ィルタと呼ぶ）を用いたイコライザのフィルタ係数を求
めるフィルタ係数演算装置に関するものである。

従来の技術 近年、ＣＤ　、ＤＡＴなどの音響機器のディジタル化に
伴ない、ＦＩＲ型のディジタルフィルタを用いたイコラ
イザの開発が求められている。

しかし従来、振幅周波数特性と位相周波数特性を一つの
ＦＩＲフィルタにより、それぞれ独立に制御することは
できなかった。

第５図に、従来の振幅周波数特性のみを制御することが
できるＦＩＲフィルタのフィルタ係数を５−一。

求めるフィルタ係数演算装置のブロック図を示す。

第５図において、２１は任意の振幅周波数特性ＩＨ（ω
）１を入力する入力部、２２は入力された振幅周波数特
性を伝達関数としてこの伝達関数を逆フーリエ変換する
ことによりフィルタ係数を求める逆フーリエ変換部、２
３は求１ったフィルタ係数をＦＩＲフィルタに設定する
設定回路、２４は実際に与えられた振幅周波数特性を実
現するＦ工Ｒフィルタである。

希望する振幅周波数特性ＩＨ（ω）１は、入力部２１に
より入力される（但し、ω−π−２πの特性はω−Ｑ−
πの入力特性を折り返して制作している）。

次に、逆フーリエ変換部２２において伝達関数Ｈ（ω）
−１［ω）１　　・・印・・川・・・・印・・・・ｆｌ
Ｊ］としてＨ（ω）を逆フーリエ変換することによりフ
ィルタ係数（Ｈ（ω）に対するインパルス応答）請求め
ることができる。

逆フーリエ変換は、次式のように実行する。

（ω＝２πに／Ｎ○≦ｎ≦ｎ−１） ６へ−７ （１２）式で求まったｈ（ｎ）はフィルタ係数として設
定回路２３によって、ＦＩＲフィルタ２４に設定され、
ここで与えられた振幅周波数特性が実現されることとな
る。

位相周波数特性は（１１）式で伝達関数を与えたことに
より直線位相となる。

発明が解決しようとする問題点 しかしながら第５図で示した従来例では、最初に述べた
ように任意の振幅周波数特性は実現できるが、位相周波
数特性はＦＩＲフィルタのタップ係数の数によシー義的
に決まる直線位相となり、位相周波数特性は任意に設定
することはできないという欠点があった。

本発明は、上記問題点に鑑み、任意の振幅周波数特性と
任意の位相周波数特性をそれぞれ独立に設定し、それを
実印できるＦＩＲフィルタのフィルタ係数を求めるフィ
ルタ係数演算装置を提供するものである。

問題点を解決するための手段 上記問題点を解決するため、本発明は、振幅特７ベー／ 性及び位相特性をそれぞれ独立に設定する振幅設定手段
及び位相設定手段と、伝達関数を求める第１の伝達関数
演算手段と、求まった伝達関数を逆フーリエ変換する第
１の逆フーリエ変換手段と、逆フーリエ変換した結果の
虚数項を零としてフーリエ変換するフーリエ変換手段と
、フーリエ変換した結果のスカラー値を求めるスカラー
演算手段と、振幅値に対する補正値を求める補正値演算
手段と、補正値演算手段によって求まった値をある値と
比較する第１の比較手段と、第１の伝達関数演算手段に
設定されている振幅特性と求１つだ補正値を乗算する第
１の乗算手段と、ツーＩＪ　工変換手段により求まった
虚数項から位相特性を求める位相演算手段と、位相設定
手段により設定された位相特性と位相演算手段により求
まった位相特性の差の比較を行う第２の比較手段と、第
２の比較手段の差がある値よシ大きい場合は、その差を
位相設定手段に設定された位相特性に加算する加算手段
と、差がある値より小さい場合は第１の伝達関数演算手
段に設定されている振幅特性と振幅設定手段により設定
された振幅特性を乗算する第２の乗算手段と、乗算され
た振幅特性を振幅項とし、第１の伝達関数演算手段に設
定されている位相特性を位相項としその伝達関数を求め
る第２の伝達関数演算手段と、求まった伝達関数を逆フ
ーリエ変換する第２の逆フーリエ変換手段と、逆フーリ
エ変換により求まった結果の実数項のみを外部のＦＩＲ
フィルタのフィルタ係数として設定する設定手段によっ
て構成される。

作用 本発明は上記した構成によシ、振幅設定手段及び位相設
定手段によシ希望する振幅周波数特性。

位相周波数特性が設定され、この振幅周波数特性及び位
相周波数特性をもつ伝達関数を逆フーリエ変換する逆フ
ーリエ変換手段によりインパルス応答が計算され、求ま
ったインパルス応答ヲフーリ工変換し振幅特性及び位相
特性を演算によシ求め、振幅設定手段及び位相設定手段
に設定された特性と演算によυ求まった特性を比較、補
正する演算を繰り返すことによシ正確な特性（インパル
ス応９、−７ 答）を実現させ、設定手段により求まったインパルス応
答がフィルタ係数として外部のＦＩＲフィルタに設定す
ることにより希望の特性が実現される。

実施例 以下、本発明の実施例について図面を参照しながら説明
する。

第１図は、本発明のフィルタ係数演算装置における一実
施例のブロック図を示すものである。

第１図において、１は希望の振幅周波数特性Ａｉ　を第
２種の標本化周波数で設定する振幅設定回路で、振幅設
定回路１では、ω−π−２πの特性は伝達関数を求める
際にω＝０−πの特性を折り返して制作している。２は
希望の群遅延歪特性τｉを第２種の標本化周波数で設定
し、それを周波数で積分して位相特性θｉを演算する位
相設定回路で、ω−π−２πの特性は伝達関数を求める
際にω＝０−πの特性を折り返して制作している。

３は位相設定回路２で求められた位相特性を位相項とし
、振幅項が初期値１の平坦特性であるもの１０・・−ノ から伝達関数Ｈ１を求める伝達関数演算回路、４は乗算
回路３で求まった伝達関数ａｉ、を第１種の標本化周波
数で逆フーリエ変換する逆フーリエ変換回路、６は逆フ
ーリエ変換回路４で逆フーリエ変換した結果の虚数項を
強制的に零として第１種の標本化周波数でフーリエ変換
するフーリエ変換回路、６はフーリエ変換で求まった実
数項と虚数項のスカラーｃ４を求めるスカラー演算回路
、７はスカラー変換回路６で求まった結果のそれぞれの
値Ｃ１に対してその値と次の値の乗算を行い、それの平
方根を求め更にその値の逆数Ｂ１を求める補正値演算回
路、８は補正値演算回路７で求まったＢｊの大きさをチ
エツクする比較回路、９は補正値演算回路７で求ま−た
Ｂｉを伝達関数演算回路３に設定されている振幅項と乗
算しその値を新しい振幅項として設定する乗算回路、１
ｏはフーリエ変換回路５において求まった虚数項の自然
対数を求め、符号を変え、各々の値に対し次の値との平
均値を求めてそれを位相特性とする位相演算回路、１１
は位相演算回路１ｏにおいて求まつ１１６−ジ た位相特性と、位相設定装置２によって求まった位相特
性を比較チエツクし、その差を求める比較回路、１２は
比較回路１１において求まった位相特性の差を伝達関数
演算回路３に設定されている位相項に加算しその値を新
し力位和項として設定する加算回路、１３は比較回路１
１において位相特性の差がある値よシ小さくなった場合
に振幅設定回路１に設定された振幅特性を伝達関数演算
回路３に設定されている振幅項と乗算する乗算回路、１
４は乗算回路１３で求まった振幅特性を振幅項とし、加
算回路１２によって調整された位相特性を位相項として
伝達関数を求める伝達関数演算回路、１５は求まった伝
達関数を逆フーリエ変換する逆フーリエ変換回路、１６
は逆フーリエ変換回路１５で求まった実数環をフィルタ
係数としてＦＩＲフィルタ１７に設定する設定回路、１
７は実際に信号のフィルタリングを行うＦＩＲフィルタ
である。

（なお、同様の機能を持つものは同じ名称を用いている
。） 以下、本実施例の動作について、図面に従って説明する
。

振幅設定回路１においてω＝Ｏ〜πの周波数範囲で希望
の振幅周波数特性Ａｉが第２種の標本化周波数毎に設定
される。第２図（２Ｌ）に第２種の標本化周波数で設定
された振幅特性の例を、第２図（１））に第１種の標本
化周波数で設定された振幅特性の例を示す。第２図に於
て、黒丸が設定される振幅値のポイントを示している。

第２図において、Ｎは実際のフィルタリングを行うＦＩ
Ｒフィルタ１７のフィルタ係数の数である。

同様に位相設定回路２においてω＝０〜πの周波数範囲
で希望の群遅延歪特性τｉが第２種の標本化周波数毎に
設定され、これを周波数で積分して位相特性θ１が求め
られる。

θ１を求める積分は、例えば次式のように行う。

１３、。

伝達関数演算回路３は振幅項を初期値１とし、位相項は
初期値を位相設定回路２で求められた位相特性を用いて
次式のようにω＝○〜２πの周波数範囲の伝達関数Ｈ１
を算出する。

ｉ＝ｏ　ＮＮ／２−１ （Ｒｅａｌ（）は実数環をｘｍａ、ｇ（）は虚数項を表
わす）伝達関数演算回路３で演算された伝達関数を逆フ
ーリエ変換回路４により第１種の標本化周波数で逆フー
リエ変換を行う（［１２１式参照）。（２）式で求まる
位相項が奇関数になれば逆フーリエ変換した結果は実数
のみとなり、これをそのままＦＩＲフィルタ１７に設定
すれば希望の周波数特性が実現されるが、位相設定回路
２で任意の群遅延歪特性が設定されるので、（２）式に
よって求まる位相特性は奇関数にはならない（ω＝πに
おいて不連続に１４Ａ−） なる）。これによシ逆フーリエ変換した結果は実数環の
みとはならず虚数項がでてくる。この結果の実数環をＦ
ＩＲフィルタ１７に設定すると虚数項の影響で実現され
る周波数特性（振幅及び位相）は設定した周波数特性と
は違ったものとなる。従って以下のような補正を行なう
。

まず、フーリエ変換回路５において逆フーリエ変換によ
って求まった結果の虚数項を強制的に零として第１種の
標本化周波数をもちいてフーリエ変換を行う。次にスカ
ラー演算回路６においてフーリエ変換によって求まった
結果の実数環Ｒ１と虚数項工ｉのスカラー（振幅値）Ｃ
１を求める。

ａｉ＝Ｊ而「「汀ｒ　　　・・・・・・・・・・・・・
・・・・・・・・（３）この値Ｃ１が逆フーリエ変換４
で求まった結果の実数環を、そのままＦＩＲフィルタに
設定したときの振幅特性に対応する（振幅設定回路１で
設定された振幅特性が平坦特性の場合）。この値Ｃ１が
振幅設定回路１で設定された振幅値と等しくなるように
伝達関数演算回路３に設定される振幅値を補正すればい
いわけである。但Ｌ７、先に述べた１５ヘーン が逆フーリエ変換した結果、そのまま実数項をＦＩＲフ
ィルタに設定したときの実現される振幅特性と設定した
振幅特性との差は設定された振幅特性が１の平坦でない
場合も平坦の場合も位相特性が等しい場合は同じである
ことが実験によってわかっているため、この補正値は振
幅特性が１の平坦な特性で求め、最終的に設定された特
性に補正値をかけたものを使用してフィルタ係数を求め
る。逆に、設定された振幅特性を目標特性として補正を
行うと正しい値に収束しない場合がある。

従って、目標特性を１としてｃｉとの差を求めれば良い
わけであるが、振幅設定及び群遅延歪の設定が第２種の
標本化周波数で設定され、逆フーリエ変換及びフーリエ
変換が第１種の標本化周波数で行われていることからそ
のままｃｌを使用すると、対応する周波数位置のずれが
生じ、正しい補正が行われない。では、なぜ第２種と第
１種の標本化周波数を混在させているかを簡単に説明す
る。通常、ある振幅特性及び位相特性から伝達関数を求
め、これを逆フーリエ変換してＦ工Ｒフィルタ１７のフ
ィルタ係数を求める場合、振幅特性上にリップルが生じ
るのを防ぐため、最終的に求まった値（逆フーリエ変換
によって求まった結果）にある適当な特性を持つ究を掛
ける。また、この窓の特性は設定される振幅特性によっ
て最適なものを使う必要がある。しかし、これは第１種
の標本化周波数で特性を設定し、第１種の標本化周波数
で逆フーリエ変換を行う場合であシ、第２種の標本化周
波数で特性を設定し、第１種の標本化周波数で逆フーリ
エ変換を行う場合（その逆も可）は、その周波数のずれ
により逆フーリエ変換を行うだけで、窓が掛かったのと
同様の効果が得られ、実現される特性も良好な（リップ
ルがない）ものとなる。従って、本発明の実施例では第
２種と第１種を混在させて演算を行っている。

そこで、次の補正値演算回路７において次式のような計
算を行い周波数ずれを考慮した補正値Ｂｉの算出を行っ
ている。

Ｂｉ＝１／＄〒　・川・川・・・・川・・（４）（ｉ＝
ｏ〜Ｎ−１、但しｉ　＝　Ｎ　−１のとき１７７、−７ ＣＭ−ｃＮ−１と１〜で計算する） この補正値の意味はこの値を対数をとって見てみると理
解しやすい。これを示したのが（５）式である。

ｂｉ＝ｏ−（Ｏｉ＋Ｏｉ＋１）／２　　・川・・・・印
・（５）つまシ、求まった振幅値ａ４の次の値との平均
値を求め、これをｏｄＢから引いて、これを補正値とす
るということである。この様にｃｉをそのまま使用せず
に次の値との平均値を用いることにより、周波数ずれの
保証が可能となる。

この補正値の求め方を示した図を第３図に示す。

第３図に於て、白丸が目標とする振幅値、黒丸がフーリ
エ変換して求めた振幅値Ｃ１、かけ印（×印）が求める
補正値Ｂｉである。

次に、比較回路８においてこのようにして求まった補正
値をチエツクし、この値がある範囲におさまっていない
場合は、乗算回路９によりこの補正値を現在伝達関数演
算回路３に設定されている振幅特性と乗算し、伝達関数
演算回路３に設定１〜なおす。第１回目の演算では初期
値として伝達量１８、、−７ 数演算回路３には振幅特性として１が設定されているこ
とからＢｉ自身が振、幅特性となるわけである。

第２回目からは、前の補正値Ｂｉに新しく求まった補正
値Ｂｉが足されていく形となる。この演算（伝達関数演
算回路３から、比較回路８及び乗算回路９までの演算）
は、比較回路８での補正値の値がある範囲以下になるま
で繰り返す（演算を繰シ返す回数を決めて行っても良い
だろう）。このようにして補正値がある範囲におさまる
と、次に位相特性の補正を開始する。

振幅特性の補正と同様に行うが、先ず位相演算回路１ｏ
においてフーリエ変換回路５において求まった結果の虚
数項の自然対数をとシ、符号を変え、各々の値に対し次
の値との平均値を求める。

この値が実現される位相特性に対応した値となる。

フーリエ変換により求まった値を、実数項Ｒｉと虚数項
工１とするとこのときの伝達関数Ｈ′ｌばｕ／ｉ＝　Ｒ
１＋コニｔ　　　・・・・・・・・山中・・・・・・・
・・印・（６）となる。また、このときの振幅値をＡ　
Ｉｉ、位相をθ′１とおくと １９ヘ−リ Ｈ’ｉ　＝　Ａ’ｉ　ｅＸｐ　（−ｊθ’ｉ　）　　・
川＝　−−−（７）と表される。

（ア）式の両辺の自然対数をとると ＩｎＨ’ｉ−：［ｎＡ’ｉ−コθ′ｉ　　　・・・・・
・・・・・・・・・・（８）同様に、（６）式の両辺の
自然対数をとるとＩｎＨ’ｉ−島Ｒｉ＋ｊ工酊ｉ　・・
曲・・・・・印・（９）となる。（８）式と（９）式を
比較すると、θ’ｉ＝　−ＩｎＩｉ　　　　・・・・・
印印・・中国・山・（１■であることがわかる。

また、振幅特性の補正の場合と同様に、第２種。

第１種の標本化周波数の関係から生じる周波数ずれを考
慮し、 θ′ｉ−（θｉ＋θｉ＋＋）／２　　・・・・・印・・
・印・（１１）（但し、１＝Ｎ−１のときθＮ＝ｏとし
て計算する） とする（以上すべて１＝ｏ−Ｎ−１）。

従って、位相演算回路１ｏのように演算を行えば、位相
特性が求まること力３分かる。次に、比較回路１１にお
いて位相演算回路１ｏによって求まった位相特性と現在
伝達関数演算回路３に設定されている位相特性を比較チ
エツクし、これらの差がある範囲におさまっていない場
合は、加算回路１２によりこの差を現在伝達関数演算回
路３に設定されている位相特性に加算し、伝達関数演算
回路３に設定しなおす。

そして、以上の演算を比較回路１１での位相特性の差７
３；ある範囲以下になるまで繰り返す（とれも演算回数
だけを決めて行っても良いだろう）。

以上の演算を繰り返し、比較回路１１での差の値がある
範囲以下になったとき、乗算回路１３において求脣った
振幅補正値と振幅設定回路１に設定された振幅特性との
乗算を行ない、これを新たな振幅項として設定し、また
このとき伝達関数演算回路３に設定されている位相特性
を位相項として次の伝達関数演算回路１４に設定する。

伝達関数演算回路１４では設定された特性から伝達関数
を求め、これを逆フーリエ変換回路１５において第１種
の標本化周波数で逆フーリエ変換し、求まった結果の実
数項を設定回路１６によってＦＩＲフィルタ１７にフィ
ルタ係数として設定２１６−ノ する。

振幅及び位相とも補正を繰り返して求まったものである
ので、振幅設定回路１及び位相設定回路２により設定さ
れた特性が正確に設定されることとなる。

第４図に以上の演算及び処理方法を示すフローチャート
を示す。

なお、本実施例では位相設定回路２は群遅延歪を入力し
て、これを積分し位相特性を求めるようにしたが、位相
特性を直接設定するようにすると、積分機能を省くこと
が可能である。また、実施例では第２種の標本化周波数
で振幅及び位相０群遅延歪）特性を設定し、第１種の標
本化周波数で逆フーリエ変換及びフーリエ変換を行うよ
うにしたが、逆に第１種の標本化周波数で振幅及び位相
０群遅延歪）特性を設定し、第２種の標本化周波数で逆
フーリエ変換及びフーリエ変換を行うようにしてもどう
ようの効果が得られる。

１だ、本実施例では先に振幅特性の補正を行い、次に位
相特性の補正を行う構成にしたが、逆に、２２、、。

先に位相特性の補正を行い、次に振幅特性の補正を行う
構成にしても効果は同様である。

発明の効果 本発明は、第２種または第１種の標本化周波数で振幅特
性及び位相特性を独立に設定する振幅設定手段及び位相
設定手段と、振幅項は平坦特性を初期値とし、位相項は
初期値を位相設定手段にょシ設定されたものとして伝達
関数を求める第１の伝達関数演算手段と、求まった伝達
関数を第１種または第２種の標本化周波数で逆フーリエ
変換する第１の逆フーリエ変換手段と、逆フーリエ変換
した結果の虚数項を零として第１種または第２種の標本
化周波数でフーリエ変換するフーリエ変換手段と、フー
リエ変換した結果のスカラー値を求めるスカラー演算手
段と、振幅特性に対する補正値を求める補正値演算手段
と、補正値演算手段によって求まった値をある値と比較
する第１の比較手段と、第１の伝達関数演算手段に設定
されている振幅特性と求まった補正値の乗算を行なう第
１の乗算手段と、フーリエ変換手段により求まった２３
へ−７ 虚数環から位相特性を求める位相演算手段と、位相設定
手段によシ設定された位相特性と位相演算手段によシ求
まった位相特性の差の比較を行う第２の比較手段と、第
２の比較手段の差がある値よシ大きい場合は、その差を
位相設定手段に設定された位相特性に加算する加算手段
と、差がある値よシ小さい場合は第１の伝達関数演算手
段に設定されている振幅特性と振幅設定手段によシ設定
された振幅特性を乗算する第２の乗算手段と、乗算され
た振幅特性を振幅項とし、第１の伝達関数演算手段に設
定されている位相特性を位相項としその伝達関数を求め
る第２の伝達関数演算手段と、求まった伝達関数を第１
種または第２種の標本化周波数で逆フーリエ変換する第
２の逆フーリエ変換手段と、逆フーリエ変換によシ求ま
った結果の実数項のみを外部のＦＩＲフィルタのフィル
タ係数として設定手段を具備することにより、振幅周波
数特性と位相周波数特性を独立に寸だ正確に制御するこ
とができるフィルタ係数を算出するフィルタ係数演算装
置を提供することができる。

【図面の簡単な説明】

第１図は本発明の実施例のフィルタ係数演算装置のブロ
ック図、第２図は第２種及び第１種の標本化周波数での
振幅特性の設定方法を示した特性図、第３図は本発明の
補正値演算手段における振幅特性の補正値を求める方法
を説明した周波数特性図、第４図は本発明の演算方法を
示しだフローチャート、第５図は従来のフィルタ係数演
算装置のブロック図である。 １・・・・・・振幅設定回路、２・・・・・・位相設定
回路、３・・・・・・伝達関数演算回路、４・・・・・
・逆フーリエ変換回路、５・・・・・・フーリエ変換回
路、６・・・・・・スカラー演算回路、７・・・・・・
補正値演算回路、８・・・・・・比較回路、１０・・・
・・・位相演算回路、９・・・・・・乗算回路−１１・
・・・・・比較回路、１２・・・・・・加算回路、１３
・・・・・・乗算回路、１４・・・・・・伝達関数演算
回路、１５・・・・・・逆フーリエ変換回路、１６・・
・・・・設定回路、１７・・・・・・ＦＩＲフィルタ。 代理人の氏名　弁理士　中　尾　敏　男　ほか１名算　
す

Claims (2)

【特許請求の範囲】
1. （１）第２種の標本化周波数（ωＫ＝２π（Ｋ＋１／２
   ）／Ｎ，Ｋ＝Ｏ〜Ｎ−１）又は第１種の標本化周波数（
   ωＫ＝２πＫ／Ｎ，Ｋ＝Ｏ〜Ｎ−１）で振幅特性及び位
   相特性をそれぞれ独立に設定する振幅設定手段及び位相
   設定手段と、位相項は初期値として前記位相設定手段に
   より設定された位相特性を、振幅項は初期値を平坦特性
   として設定し伝達関数を求める第１の伝達関数演算手段
   と、振幅特性及び位相特性を第２種の標本化周波数で設
   定した場合は第１種で、振幅特性及び位相特性を第１種
   の標本化周波数で設定した場合は第２種の標本化周波数
   で、求まった伝達関数を逆フーリエ変換する第１の逆フ
   ーリエ変換手段と、前記逆フーリエ変換した結果の虚数
   項を零として振幅特性及び位相特性を第２種の標本化周
   波数で設定した場合は第１種で、振幅特性及び位相特性
   を第１種の標本化周波数で設定した場合は第２種の標本
   化周波数で、フーリエ変換するフーリエ変換手段と、前
   記フーリエ変換した結果の実数項と虚数項のスカラー値
   を求めるスカラー演算手段と、求まったスカラー値全部
   に対して各々の標本化周波数に対するスカラー値とその
   次の標本化周波数のスカラー値の積を行い、その結果の
   平方根を求め、更にその逆数を求める補正値演算手段と
   、前記補正値演算手段によって求まった値をある値と比
   較する第１の比較手段と、補正値がある値より大きい場
   合は第１の伝達関数演算手段に設定されている振幅特性
   と求まった補正値を乗算する第１の乗算手段と、乗算さ
   れた振幅値を振幅項として第１の伝達関数演算手段に設
   定し再び上記の演算を繰り返し、補正値がある値より小
   さい場合は前記フーリエ変換手段により求まった結果の
   虚数項の自然対数をとり、符号を変え、更に各々の標本
   化周波数に対する位相値に対しその次の標本化周波数の
   値との平均値を求め、それを位相特性とする位相演算手
   段と、前記位相設定手段により設定された位相特性と位
   相演算手段により求まった位相特性の差を求め、その差
   の比較を行う第２の比較手段と、前記第２の比較手段の
   差がある値より大きい場合は、その差を位相設定手段に
   設定された位相特性に加算する加算手段と、前記加算手
   段により求まった位相特性を位相項として第１の伝達関
   数演算手段に設定し再び上記の演算を繰り返し、差があ
   る値より小さい場合は第１の伝達関数演算手段に設定さ
   れている振幅特性と振幅設定手段により設定された振幅
   特性を乗算する第２の乗算手段と、乗算された振幅特性
   を振幅項とし、前記第１の伝達関数演算手段に設定され
   ている位相特性を位相項とし、その伝達関数を求める第
   ２の伝達関数演算手段と、求まった伝達関数を逆フーリ
   エ変換する第２の逆フーリエ変換手段と、逆フーリエ変
   換により求まった結果の実数項のみを外部のＦＩＲフィ
   ルタのフィルタ係数として設定する設定手段とを具備し
   たフィルタ係数演算装置。
2. （２）位相設定手段は、群遅延歪特性を先ず設定し、そ
   れを周波数で積分してその値を位相特性とすることを特
   徴とする特許請求の範囲第１項記載のフィルタ係数演算
   装置。