JPH0770956B2 - 周波数特性設定装置 - Google Patents
周波数特性設定装置Info
- Publication number
- JPH0770956B2 JPH0770956B2 JP11305087A JP11305087A JPH0770956B2 JP H0770956 B2 JPH0770956 B2 JP H0770956B2 JP 11305087 A JP11305087 A JP 11305087A JP 11305087 A JP11305087 A JP 11305087A JP H0770956 B2 JPH0770956 B2 JP H0770956B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- frequency characteristic
- sampling
- input
- setting device
- 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
Landscapes
- Filters That Use Time-Delay Elements (AREA)
- Tone Control, Compression And Expansion, Limiting Amplitude (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Description
【発明の詳細な説明】 産業上の利用分野 本発明は、任意の振幅又は位相又は群遅延周波数特性な
どの信号の周波数特性を実現するトランスバーサル・フ
ィルタ（以下、FIRフィルタと呼ぶ）を用いたディジタ
ルイコライザの周波数特性を入力し、そのフィルタ係数
を求める周波数特性設定装置に関するものである。
どの信号の周波数特性を実現するトランスバーサル・フ
ィルタ（以下、FIRフィルタと呼ぶ）を用いたディジタ
ルイコライザの周波数特性を入力し、そのフィルタ係数
を求める周波数特性設定装置に関するものである。
従来の技術 近年、音響機器のディジタル化にともない、FIRフィル
タを用いたイコライザの開発が求められている。
タを用いたイコライザの開発が求められている。
第６図は従来の周波数特性設定装置のブロック図であ
る。
る。
６は振幅周波数特性を入力する入力回路、７は入力回路
６に入力された振幅周波数特性からFIRフィルタのフィ
ルタ係数を求める演算回路、８は実際に入力された振幅
周波数特性を実現するFIRフィルタである。
６に入力された振幅周波数特性からFIRフィルタのフィ
ルタ係数を求める演算回路、８は実際に入力された振幅
周波数特性を実現するFIRフィルタである。
入力回路６で入力される振幅周波数特性は処理される信
号の標本化周波数をＦ、FIRフィルタのフィルタ係数の
数をＮとすると、周波数分解能がF/Nで周波数０からF/2
まで入力される。そして、演算回路７において入力され
た振幅周波数特性から折り返し（周波数F/2〜Ｆの帯域
の振幅周波数特性値）をつくり、この値を実数項とし、
虚数項を零として、逆フーリエ変換を行うことにより入
力された振幅周波数特性を実現するFIRフィルタ８のフ
ィルタ係数が求められる。そして、このフィルタ係数が
FIRフィルタ８に設定されることにより実際に入力され
た振幅周波数特性が実現されることとなる。
号の標本化周波数をＦ、FIRフィルタのフィルタ係数の
数をＮとすると、周波数分解能がF/Nで周波数０からF/2
まで入力される。そして、演算回路７において入力され
た振幅周波数特性から折り返し（周波数F/2〜Ｆの帯域
の振幅周波数特性値）をつくり、この値を実数項とし、
虚数項を零として、逆フーリエ変換を行うことにより入
力された振幅周波数特性を実現するFIRフィルタ８のフ
ィルタ係数が求められる。そして、このフィルタ係数が
FIRフィルタ８に設定されることにより実際に入力され
た振幅周波数特性が実現されることとなる。
発明が解決しようとする問題点 しかしながら、第６図で示した従来例では処理を行う周
波数帯域と周波数分解能（FIRフィルタのフィルタ係数
の数）が決まると演算回路７での演算量（特に逆フーリ
エ変換の計算時間）が決まりこれを短縮できないという
欠点があった。
波数帯域と周波数分解能（FIRフィルタのフィルタ係数
の数）が決まると演算回路７での演算量（特に逆フーリ
エ変換の計算時間）が決まりこれを短縮できないという
欠点があった。
本発明は上記問題点を鑑みなされたもので、フィルタ係
数を求める演算量を少なくし、演算時間を大幅に短縮で
きるようにした周波数特性設定装置を提供するものであ
る。
数を求める演算量を少なくし、演算時間を大幅に短縮で
きるようにした周波数特性設定装置を提供するものであ
る。
問題点を解決するための手段 上記問題点を解決するため本発明は、振幅又は位相又は
群遅延などの周波数特性を入力する入力手段と、入力手
段により入力された周波数特性を複数の異なる標本化周
波数に対応した周波数帯域に分割する分割手段と、分割
手段により分割されたそれぞれの周波数帯域の周波数特
性を実現するトランスバーサル・フィルタのフィルタ係
数を求める演算手段によって構成される。
群遅延などの周波数特性を入力する入力手段と、入力手
段により入力された周波数特性を複数の異なる標本化周
波数に対応した周波数帯域に分割する分割手段と、分割
手段により分割されたそれぞれの周波数帯域の周波数特
性を実現するトランスバーサル・フィルタのフィルタ係
数を求める演算手段によって構成される。
作用 本発明は、入力手段により振幅又は位相又は群遅延など
の信号の周波数特性が入力され、分割手段によって異な
る複数の標本化周波数に対応する周波数帯域に入力され
た周波数特性を分割し、演算手段により分割されたそれ
ぞれの周波数特性を実現するFIRフィルタのフィルタ係
数を求めることにより、フィルタ係数を求める演算時間
を大幅に短縮できることとなる。
の信号の周波数特性が入力され、分割手段によって異な
る複数の標本化周波数に対応する周波数帯域に入力され
た周波数特性を分割し、演算手段により分割されたそれ
ぞれの周波数特性を実現するFIRフィルタのフィルタ係
数を求めることにより、フィルタ係数を求める演算時間
を大幅に短縮できることとなる。
実 施 例 以下、本発明の実施例について図面を参照しながら説明
する。第１図は、本発明の第１の実施例における周波数
特性設定装置のブロック図を示すものである。
する。第１図は、本発明の第１の実施例における周波数
特性設定装置のブロック図を示すものである。
第１図において、１は任意の振幅周波数特性を２つの標
本化周波数に対する周波数帯域に分けて入力する入力回
路、２は入力された振幅周波数特性を２つの異なる標本
化周波数に対応する周波数帯域に分割する分割回路、３
はそれぞれの標本化周波数の周波数帯域に対応するFIR
フィルタのフィルタ係数を求める演算回路、４と５は標
本化周波数が異なりフィルタ係数の数が同じであるFIR
フィルタである。
本化周波数に対する周波数帯域に分けて入力する入力回
路、２は入力された振幅周波数特性を２つの異なる標本
化周波数に対応する周波数帯域に分割する分割回路、３
はそれぞれの標本化周波数の周波数帯域に対応するFIR
フィルタのフィルタ係数を求める演算回路、４と５は標
本化周波数が異なりフィルタ係数の数が同じであるFIR
フィルタである。
以下、本実施例の動作について説明する。
入力回路１において、任意の振幅周波数特性を２つの標
本化周波数に対する周波数帯域に分け、その周波数帯域
が重なる部分は周波数分解能が小さい方で、周波数帯域
が重ならない周波数帯域では周波数分解能が大きい方で
振幅周波数特性が入力される。
本化周波数に対する周波数帯域に分け、その周波数帯域
が重なる部分は周波数分解能が小さい方で、周波数帯域
が重ならない周波数帯域では周波数分解能が大きい方で
振幅周波数特性が入力される。
第２図ａに入力回路１で入力される振幅周波数特性を示
す。第２図ａにおいて丸（○）印は、設定された振幅値
を示し、f1,f2は２つの異なる標本化周波数であり、f2
はf1の２倍であるとする。入力回路１で入力できる振幅
周波数特性の周波数分解能はこのf1/2で二分される。こ
の時、入力回路１によって振幅周波数特性を入力できる
f1/2までの周波数帯域の周波数分解能はその周波数帯域
に対応するFIRフィルタ４のフィルタ係数の数をＮとす
るとf1/Nとなり、f1/2からf2/2の周波数帯域における入
力できる振幅周波数特性の周波数分解能はf2/Nとなる。
但し、実際に標本化周波数f1で信号の処理を行うFIRフ
ィルタ４の低域通過フィルタのカットオフ周波数付近の
周波数がf3であった場合は（周波数０からf2/Nの分解能
で周波数ポイントをとった時に、低域通過フィルタのカ
ットオフ周波数と一致するかそれに近い周波数をf3とす
る）、f3からf1/2までの設定できる振幅周波数特性の周
波数分解能は第３図のようにf2/Nとする。なお本実施例
はこのような入力方法をとるものとする。
す。第２図ａにおいて丸（○）印は、設定された振幅値
を示し、f1,f2は２つの異なる標本化周波数であり、f2
はf1の２倍であるとする。入力回路１で入力できる振幅
周波数特性の周波数分解能はこのf1/2で二分される。こ
の時、入力回路１によって振幅周波数特性を入力できる
f1/2までの周波数帯域の周波数分解能はその周波数帯域
に対応するFIRフィルタ４のフィルタ係数の数をＮとす
るとf1/Nとなり、f1/2からf2/2の周波数帯域における入
力できる振幅周波数特性の周波数分解能はf2/Nとなる。
但し、実際に標本化周波数f1で信号の処理を行うFIRフ
ィルタ４の低域通過フィルタのカットオフ周波数付近の
周波数がf3であった場合は（周波数０からf2/Nの分解能
で周波数ポイントをとった時に、低域通過フィルタのカ
ットオフ周波数と一致するかそれに近い周波数をf3とす
る）、f3からf1/2までの設定できる振幅周波数特性の周
波数分解能は第３図のようにf2/Nとする。なお本実施例
はこのような入力方法をとるものとする。
以上のようにして入力された振幅周波数特性は次の分割
回路２に入力され、２つの標本化周波数及びFIRフィル
タ4,5のフィルタ係数の数に対応した周波数特性第２図
b,cに分割される（×印，△印がその振幅値）。以下、
この分解回路２における振幅周波数特性の分割方法につ
いて述べる。
回路２に入力され、２つの標本化周波数及びFIRフィル
タ4,5のフィルタ係数の数に対応した周波数特性第２図
b,cに分割される（×印，△印がその振幅値）。以下、
この分解回路２における振幅周波数特性の分割方法につ
いて述べる。
まず、標本化周波数f1に対応する帯域０〜f1/2におい
て、０〜f3までの振幅周波数特性は周波数分解能f1/Nで
入力回路１において入力されるのでそのままこの値を用
いる。標本化周波数f1に対応する周波数帯域で問題とな
るのはf3〜f1/2までの帯域であるが、この帯域は周波数
分解能がf2/Nで振幅周波数特性が入力されていることか
らこの振幅周波数特性の値からf1/Nの周波数分解能の新
しい振幅周波数ポイントを求める必要がある。この新し
い振幅周波数ポイントの求めかたにはいろいろな方法が
考えられるが、ここでは線形補間して求める方法を示
す。つまり入力回路１において入力されたf3〜f1/2の振
幅周波数の値のうちf3での値をA,次の値をＢとするとA,
B間にある本来の標本化周波数f1での値Ｃは、標本化周
波数f2がf1の２倍であるので Ｃ＝（Ｂ−Ａ）×f1/f2＋Ａ として求めることができる。f1/2までの他の値は同様に
して求めることができる。なお、ここでは線形補間をお
こなって標本化周波数f1に対応した周波数帯域の振幅周
波数特性を求めたが、先に述べたように他の方法、例え
ば多次関数による方法などいろいろな方法が考えられ
る。
て、０〜f3までの振幅周波数特性は周波数分解能f1/Nで
入力回路１において入力されるのでそのままこの値を用
いる。標本化周波数f1に対応する周波数帯域で問題とな
るのはf3〜f1/2までの帯域であるが、この帯域は周波数
分解能がf2/Nで振幅周波数特性が入力されていることか
らこの振幅周波数特性の値からf1/Nの周波数分解能の新
しい振幅周波数ポイントを求める必要がある。この新し
い振幅周波数ポイントの求めかたにはいろいろな方法が
考えられるが、ここでは線形補間して求める方法を示
す。つまり入力回路１において入力されたf3〜f1/2の振
幅周波数の値のうちf3での値をA,次の値をＢとするとA,
B間にある本来の標本化周波数f1での値Ｃは、標本化周
波数f2がf1の２倍であるので Ｃ＝（Ｂ−Ａ）×f1/f2＋Ａ として求めることができる。f1/2までの他の値は同様に
して求めることができる。なお、ここでは線形補間をお
こなって標本化周波数f1に対応した周波数帯域の振幅周
波数特性を求めたが、先に述べたように他の方法、例え
ば多次関数による方法などいろいろな方法が考えられ
る。
また、周波数f4（＞f3）からf1/2においては、f1/2の振
幅周波数特性値が零となるように補正を行った値を求
め、これを使用する方法も考えられる。この方法を示し
たのが第４図である。このf4の値は先に述べた低域通過
フィルタの信号抑制レベルが十分な値（例えば通過帯域
のレベルから10dB落ちた周波数）になる周波数とする。
周波数f4における入力回路１により入力された振幅値を
Ｄとすると新しい標本化周波数f1での値Ｅは Ｅ＝−D/（f1/2−f4）×f1/N＋Ｄ として求めることができ、多の値も同様にして求めるこ
とができる。このようにすることにより、先に述べた低
域通過フィルタの阻止帯域の信号を更に抑制することが
可能となる。また、f1/2において振幅周波数特性を零と
することは、演算回路３においてフィルタ係数をリニア
フェーズの条件で求める際に必要な条件でもある。
幅周波数特性値が零となるように補正を行った値を求
め、これを使用する方法も考えられる。この方法を示し
たのが第４図である。このf4の値は先に述べた低域通過
フィルタの信号抑制レベルが十分な値（例えば通過帯域
のレベルから10dB落ちた周波数）になる周波数とする。
周波数f4における入力回路１により入力された振幅値を
Ｄとすると新しい標本化周波数f1での値Ｅは Ｅ＝−D/（f1/2−f4）×f1/N＋Ｄ として求めることができ、多の値も同様にして求めるこ
とができる。このようにすることにより、先に述べた低
域通過フィルタの阻止帯域の信号を更に抑制することが
可能となる。また、f1/2において振幅周波数特性を零と
することは、演算回路３においてフィルタ係数をリニア
フェーズの条件で求める際に必要な条件でもある。
また同様に標本化周波数f2の帯域も、処理する信号に影
響しない周波数f5から同様な補正を行なうと同様の効果
が得られる。
響しない周波数f5から同様な補正を行なうと同様の効果
が得られる。
次に、標本化周波数f2に対応する周波数帯域０〜f2/2に
おける振幅周波数特性について述べる。
おける振幅周波数特性について述べる。
周波数f3からf2/2までの値は入力回路１において周波数
分解能f2/Nで入力されるためこの値をそのまま使用でき
る。周波数０からf3までの値が周波数分解能がちがうた
め入力振幅周波数特性をそのまま使用できないため演算
を行ってこの帯域の値を求める必要がある。この演算の
方法も色々と考えられるが、ここでは入力された振幅周
波数特性を間引く方法を示す。つまり、ここではf2がf1
の２倍であることから周波数０〜f3までの標本化周波数
f2、フィルタ係数Ｎでの振幅周波数特性は、周波数分解
能f1/Nで入力された振幅周波数特性を２個に１個間引く
ことにより求めることができるわけである（第２図Ｃを
参照）このように間引きを行うことにより入力回路１で
周波数分解能f1/Nで入力された振幅周波数特性から標本
化周波数f2に対応する周波数帯域の振幅周波数特性を求
めることができる。
分解能f2/Nで入力されるためこの値をそのまま使用でき
る。周波数０からf3までの値が周波数分解能がちがうた
め入力振幅周波数特性をそのまま使用できないため演算
を行ってこの帯域の値を求める必要がある。この演算の
方法も色々と考えられるが、ここでは入力された振幅周
波数特性を間引く方法を示す。つまり、ここではf2がf1
の２倍であることから周波数０〜f3までの標本化周波数
f2、フィルタ係数Ｎでの振幅周波数特性は、周波数分解
能f1/Nで入力された振幅周波数特性を２個に１個間引く
ことにより求めることができるわけである（第２図Ｃを
参照）このように間引きを行うことにより入力回路１で
周波数分解能f1/Nで入力された振幅周波数特性から標本
化周波数f2に対応する周波数帯域の振幅周波数特性を求
めることができる。
なお、ここでは間引きを行って標本化周波数f2に対応し
た周波数帯域の振幅周波数特性を求めたが、先に述べた
ように他の方法、例えば０〜f6（＜f3）の振幅周波数特
性を零にする、または、0db（＝1.0）にする方法などが
考えられる。
た周波数帯域の振幅周波数特性を求めたが、先に述べた
ように他の方法、例えば０〜f6（＜f3）の振幅周波数特
性を零にする、または、0db（＝1.0）にする方法などが
考えられる。
次に演算回路３において分割回路２で分割された２つの
周波数特性を使用してFIRフィルタ4,5のフィルタ係数が
求められる。ここでの演算は、それぞれの周波数帯域の
折り返し（f1/2〜f1,f2/2〜f2の周波数帯域）の振幅周
波数特性を求め、リニアフェーズの条件から振幅周波数
特性を実数項とし、虚数項を零とした逆フーリエ変換を
行いフィルタ係数を求める方法などがとられる。また、
虚数項を零とはせずに、ある位相周波数特性あるいは群
遅延特性を示す虚数項を与え、逆フーリエ変換を行う方
法なども考えられる。またヒルベルト変換の関係から最
小位相推移を行なうようなフィルタ係数演算法などもと
られる。このようにすることにより、それぞれの標本化
周波数に対応する周波数帯域に対して、入力された振幅
周波数特性（または、振幅周波数特性と位相周波数特性
又は群遅延特性）を実現するフィルタ係数が求められ
る。このようにして求められたフィルタ係数は、それぞ
れの周波数帯域に対応したFIRフィルタ4,5に入力され、
実際に入力回路１に入力された振幅周波数特性が実現さ
れることとなる。
周波数特性を使用してFIRフィルタ4,5のフィルタ係数が
求められる。ここでの演算は、それぞれの周波数帯域の
折り返し（f1/2〜f1,f2/2〜f2の周波数帯域）の振幅周
波数特性を求め、リニアフェーズの条件から振幅周波数
特性を実数項とし、虚数項を零とした逆フーリエ変換を
行いフィルタ係数を求める方法などがとられる。また、
虚数項を零とはせずに、ある位相周波数特性あるいは群
遅延特性を示す虚数項を与え、逆フーリエ変換を行う方
法なども考えられる。またヒルベルト変換の関係から最
小位相推移を行なうようなフィルタ係数演算法などもと
られる。このようにすることにより、それぞれの標本化
周波数に対応する周波数帯域に対して、入力された振幅
周波数特性（または、振幅周波数特性と位相周波数特性
又は群遅延特性）を実現するフィルタ係数が求められ
る。このようにして求められたフィルタ係数は、それぞ
れの周波数帯域に対応したFIRフィルタ4,5に入力され、
実際に入力回路１に入力された振幅周波数特性が実現さ
れることとなる。
なお、実施例では入力回路１において振幅周波数特性を
第１種（周波数０からf1/N毎に、f3からはf2/N毎に周波
数ポイントをとる）の標本化を行った場合を示したが、
第２種（周波数f1/N×1/2からf1/N毎に、f3（この場
合、周波数f2/N×1/2からf2/N毎に周波数ポイントをと
っていったときに低域通過フィルタのカットオフ周波数
に一致するかそれに近い値をf3とする）からはf2/N毎に
周波数ポイントをとる）の標本化を行った場合も、周波
数ポイントが周波数帯域の違いでずれが生じる場合があ
るが、直線補間又は多次関数近似などの方法をとること
により同様の考え方で周波数分割を行うことができ、そ
れぞれの周波数帯域のフィルタ係数を求めることができ
る。
第１種（周波数０からf1/N毎に、f3からはf2/N毎に周波
数ポイントをとる）の標本化を行った場合を示したが、
第２種（周波数f1/N×1/2からf1/N毎に、f3（この場
合、周波数f2/N×1/2からf2/N毎に周波数ポイントをと
っていったときに低域通過フィルタのカットオフ周波数
に一致するかそれに近い値をf3とする）からはf2/N毎に
周波数ポイントをとる）の標本化を行った場合も、周波
数ポイントが周波数帯域の違いでずれが生じる場合があ
るが、直線補間又は多次関数近似などの方法をとること
により同様の考え方で周波数分割を行うことができ、そ
れぞれの周波数帯域のフィルタ係数を求めることができ
る。
また、実施例では振幅周波数特性を入力する場合を示し
たが入力が位相周波数特性又は群遅延特性などの他の周
波数特性とした場合も同様の考え方で行うことができ
る。
たが入力が位相周波数特性又は群遅延特性などの他の周
波数特性とした場合も同様の考え方で行うことができ
る。
第５図は、本発明の第２の実施例における入力回路の振
幅周波数特性の入力方法及び分割回路の周波数分割の方
法を示す図である。第２の実施例の構成は第１図と同様
であり、入力方法，分割方法のみが異なる。
幅周波数特性の入力方法及び分割回路の周波数分割の方
法を示す図である。第２の実施例の構成は第１図と同様
であり、入力方法，分割方法のみが異なる。
入力回路１で入力できる振幅周波数特性（第５図ａ）
は、全帯域にわたって処理できる標本化周波数f1,f2と
は異なった（この図の場合は、f1/2Nの）周波数分解能
で設定できるようになっている。次の分割回路２では、
実際の信号を処理できる周波数ポイントを得るために、
入力された周波数ポイントから、f1の標本化周波数に対
応する周波数帯域では、入力された周波数ポイントを２
個のうち１個を間引く演算を（第５図ｂ）、f2の標本化
周波数に対応する周波数帯域では、入力された周波数ポ
イントを４個のうち１個を間引く演算を行ない（第５図
ｃ）、入力された振幅周波数特性を分割する。
は、全帯域にわたって処理できる標本化周波数f1,f2と
は異なった（この図の場合は、f1/2Nの）周波数分解能
で設定できるようになっている。次の分割回路２では、
実際の信号を処理できる周波数ポイントを得るために、
入力された周波数ポイントから、f1の標本化周波数に対
応する周波数帯域では、入力された周波数ポイントを２
個のうち１個を間引く演算を（第５図ｂ）、f2の標本化
周波数に対応する周波数帯域では、入力された周波数ポ
イントを４個のうち１個を間引く演算を行ない（第５図
ｃ）、入力された振幅周波数特性を分割する。
以下は第１の実施例と同様の処理を行なう。
本実施例のように、入力回路１においては、処理できる
周波数分解能とは異なった分解能で周波数特性を設定で
きるようにしても、分割回路２において適当な演算を行
なうことにより、実際に処理できる周波数特性を得るこ
とができ、これによりFIRフィルタ4,5のフィルタ係数も
求めることが可能となる。
周波数分解能とは異なった分解能で周波数特性を設定で
きるようにしても、分割回路２において適当な演算を行
なうことにより、実際に処理できる周波数特性を得るこ
とができ、これによりFIRフィルタ4,5のフィルタ係数も
求めることが可能となる。
なお、本実施例では、分割回路２では第１種の標本化を
行なった場合を示したが、第２種の標本化を行なう場合
も考えられる。
行なった場合を示したが、第２種の標本化を行なう場合
も考えられる。
また、第1,第２の実施例において、第１種の標本化を行
ない、第１種の標本化周波数によるIFFTを行なった場合
は、通常、演算回路３において求めたフィルタ係数に窓
をかける必要がある。したがって、第１図の演算回路３
において、何らかの特性をもつ窓関数をかける演算を行
なう機能をもたせる構成も考えられる。ただし、第２種
の標本化を行ない、第１種の標本化周波数によるIFFTを
行なった場合は、特に窓をかける必要はないことを明記
しておく。
ない、第１種の標本化周波数によるIFFTを行なった場合
は、通常、演算回路３において求めたフィルタ係数に窓
をかける必要がある。したがって、第１図の演算回路３
において、何らかの特性をもつ窓関数をかける演算を行
なう機能をもたせる構成も考えられる。ただし、第２種
の標本化を行ない、第１種の標本化周波数によるIFFTを
行なった場合は、特に窓をかける必要はないことを明記
しておく。
発明の効果 以上述べたように、本発明の周波数特性設定装置は、入
力された任意の振幅又は位相又は群遅延などの周波数特
性を、複数の異なる標本化周波数に対応した周波数帯域
に分割する構成をとることにより、それぞれの周波数帯
域における処理するデータ数を減少させることができ、
これにより、フィルタ係数を求める演算時間を大幅に短
縮できることとなる。
力された任意の振幅又は位相又は群遅延などの周波数特
性を、複数の異なる標本化周波数に対応した周波数帯域
に分割する構成をとることにより、それぞれの周波数帯
域における処理するデータ数を減少させることができ、
これにより、フィルタ係数を求める演算時間を大幅に短
縮できることとなる。
第１図は本発明の第１の実施例における周波数特性設定
装置のブロック図、第２図は本発明の第１の実施例にお
ける入力周波数特性及び分割周波数特性を示す特性図、
第３図は本発明の第１の実施例における他の入力周波数
特性を示す特性図、第４図は本発明の第１の実施例にお
ける分割周波数特性を補正する方法を示す特性図、第５
図は本発明の第２の実施例における入力周波数特性及び
分割周波数特性を示す特性図、第６図は従来例を示すブ
ロック図である。 １……入力回路、２……分割回路、３……演算回路、4,
5……FIRフィルタ。
装置のブロック図、第２図は本発明の第１の実施例にお
ける入力周波数特性及び分割周波数特性を示す特性図、
第３図は本発明の第１の実施例における他の入力周波数
特性を示す特性図、第４図は本発明の第１の実施例にお
ける分割周波数特性を補正する方法を示す特性図、第５
図は本発明の第２の実施例における入力周波数特性及び
分割周波数特性を示す特性図、第６図は従来例を示すブ
ロック図である。 １……入力回路、２……分割回路、３……演算回路、4,
5……FIRフィルタ。
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭59−52912（ＪＰ，Ａ) 特開 昭62−82708（ＪＰ，Ａ) 特開 昭62−95008（ＪＰ，Ａ) 特開 昭63−246915（ＪＰ，Ａ) 特開 昭63−246914（ＪＰ，Ａ) 特開 昭63−244924（ＪＰ，Ａ) 特開 昭63−234617（ＪＰ，Ａ) 特公 平６−91416（ＪＰ，Ｂ２) 米国特許4896285（ＵＳ，Ａ) 欧州特許出願公開284175（ＥＰ，Ａ) 欧州特許出願公開544647（ＥＰ，Ａ)
Claims (11)
- 【請求項１】任意の振幅又は位相又は群遅延などの周波
数特性を入力する入力手段と、前記入力手段により入力
された周波数特性を複数の異なる標本化周波数に対応し
た周波数帯域に分割する分割手段と、前記分割手段によ
り分割されたそれぞれの周波数帯域の周波数特性を実現
するトランスバーサル・フィルタのフィルタ係数を求め
る演算手段を具備した周波数特性設定装置。 - 【請求項２】入力手段は、それぞれのトランスバーサル
・フィルタのフィルタ係数の数に対応した周波数分解能
で周波数特性を入力することを特徴とする特許請求の範
囲第１項記載の周波数特性設定装置。 - 【請求項３】入力手段は、複数の異なる標本化周波数に
対応した周波数帯域が重なる部分においては、それぞれ
の周波数帯域に対するトランスバーサル・フィルタのフ
ィルタ係数の数に対応した周波数分解能が小さい方で入
力し、入力された周波数特性を分割する分割手段は、周
波数帯域が重なる周波数分解能が大きい周波数特性を、
前記入力手段により入力された重なった周波数分解能が
小さい周波数特性から演算を行なうことにより求め、周
波数特性を分割することを特徴とする特許請求の範囲第
１項記載の周波数特性設定装置。 - 【請求項４】入力手段は、複数の異なる標本化周波数に
対応した周波数帯域が重なる部分の境界付近において
は、それぞれの周波数帯域に対するトランスバーサル・
フィルタのフィルタ係数の数に対応した周波数分解能が
大きい方で入力し、入力された周波数特性を分割する分
割手段は、周波数帯域が重なる周波数分解能が小さい周
波数特性のうち他の、周波数帯域と重なる部分の境界付
近の周波数特性を前記入力手段により入力された境界付
近の重なった周波数分解能が大きい周波数特性から演算
を行うことにより求め、周波数特性を分割することを特
徴とする特許請求の範囲第１項記載の周波数特性設定装
置。 - 【請求項５】入力手段は、トランスバーサル・フィルタ
のフィルタ係数の数に無関係な周波数分解能で入力し、
入力された周波数特性をそれぞれの前記トランスバーサ
ル・フィルタのフィルタ係数の数に対応した周波数分解
能での値になるように演算して設定することを特徴とす
る特許請求の範囲第１項記載の周波数特性設定装置。 - 【請求項６】入力手段により入力された周波数特性を分
割する分割手段は、実際に信号を処理するそれぞれのト
ランスバーサル・フィルタに入力される信号を周波数分
割する帯域通過フィルタの高域のカットオフ周波数以上
の周波数から前記トランスバーサル・フィルタに対応す
る周波数帯域のナイキスト周波数（標本化周波数の半分
の周波数）に向かって周波数特性が零となるように補正
を行い、分割することを特徴とする特許請求の範囲第１
項記載の周波数特性設定装置。 - 【請求項７】入力手段により入力された周波数特性を分
割する分割手段は、異なる複数の標本化周波数に対応す
る周波数帯域のナイキスト周波数（標本化周波数の半分
の周波数）において周波数特性が零となるように補正を
行い、分割することを特徴とする特許請求の範囲第１項
記載の周波数特性設定装置。 - 【請求項８】振幅又は位相又は群遅延などの周波数特性
を入力する入力手段は、第１種（周波数の標本化ポイン
トを周波数０から標本化周波数とトランスバーサル・フ
ィルタのフィルタ係数できまる周波数分解能毎にとって
いく方法）の標本化で行うことを特徴とする特許請求の
範囲第１項記載の周波数特性設定装置。 - 【請求項９】振幅又は位相又は群遅延などの周波数特性
を入力する入力手段は、第２種（周波数の標本化ポイン
トを周波数０からではなく、標本化周波数とトランスバ
ーサル・フィルタのフィルタ係数できまる周波数分解能
の半分の周波数からその周波数分解能毎にとっていく方
法）の標本化で行うことを特徴とする特許請求の範囲第
１項記載の周波数特性設定装置。 - 【請求項１０】振幅又は位相又は群遅延などの周波数特
性を入力する入力手段は、第１種，第２種の標本化に無
関係に入力し、分割手段は入力手段により入力された周
波数特性から、第１種あるいは第２種の標本化ポイント
を演算して求め、分割することを特徴とする特許請求の
範囲第１項記載の周波数特性設定装置。 - 【請求項１１】入力手段又は、分割手段により設定され
た周波数特性が第１種の標本化で行なわれていた場合
は、演算手段において、求められたフィルタ係数に対し
て何らかの特性を持つ窓関数をかけて、この窓をかけた
ものをフィルタ係数とすることを特徴とする特許請求の
範囲第１項記載の周波数特性設定装置。
Priority Applications (8)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11305087A JPH0770956B2 (ja) | 1987-05-08 | 1987-05-08 | 周波数特性設定装置 |
| EP88300324A EP0284175B1 (en) | 1987-03-23 | 1988-01-15 | Calculation of filter factors for digital filter |
| DE3856040T DE3856040T2 (de) | 1987-03-23 | 1988-01-15 | Filterkoeffizientenberechnung für ein digitales Filter |
| EP93102770A EP0544647B1 (en) | 1987-03-23 | 1988-01-15 | Calculation of filter factors for digital filter |
| US07/144,477 US4896285A (en) | 1987-03-23 | 1988-01-15 | Calculation of filter factors for digital filter |
| DE3853372T DE3853372T2 (de) | 1987-03-23 | 1988-01-15 | Filterkoeffizientenberechnung für ein digitales Filter. |
| CA000556581A CA1280808C (en) | 1987-03-23 | 1988-01-15 | Calculation of filter factors for digital filter |
| KR1019880003042A KR910007021B1 (ko) | 1986-07-18 | 1988-03-22 | 필터계수 연산장치 및 필터계수의 연산설정방법 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11305087A JPH0770956B2 (ja) | 1987-05-08 | 1987-05-08 | 周波数特性設定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS63278410A JPS63278410A (ja) | 1988-11-16 |
| JPH0770956B2 true JPH0770956B2 (ja) | 1995-07-31 |
Family
ID=14602225
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11305087A Expired - Fee Related JPH0770956B2 (ja) | 1986-07-18 | 1987-05-08 | 周波数特性設定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0770956B2 (ja) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4896285A (en) | 1987-03-23 | 1990-01-23 | Matsushita Electric Industrial Co., Ltd. | Calculation of filter factors for digital filter |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5952912A (ja) * | 1982-09-20 | 1984-03-27 | Nec Corp | デイジタル周波数特性可変装置 |
| JPS6282708A (ja) * | 1985-10-07 | 1987-04-16 | Victor Co Of Japan Ltd | デジタル・グラフイツク・イコライザ |
| JPS6295008A (ja) * | 1985-10-22 | 1987-05-01 | Victor Co Of Japan Ltd | デジタル・グラフイツク・イコライザ |
-
1987
- 1987-05-08 JP JP11305087A patent/JPH0770956B2/ja not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4896285A (en) | 1987-03-23 | 1990-01-23 | Matsushita Electric Industrial Co., Ltd. | Calculation of filter factors for digital filter |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS63278410A (ja) | 1988-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4649507A (en) | Segmented transversal filter | |
| JPH0759186A (ja) | 音響信号の線形歪補償方法及びその装置 | |
| US4689759A (en) | Process and installation for the analysis and retrieval of a sampling and interpolation signal | |
| EP0450720B1 (en) | Ghost cancellation circuit | |
| JPH09205345A (ja) | オーディオ信号のフィルタリングのための方法及び装置 | |
| JPH0770956B2 (ja) | 周波数特性設定装置 | |
| JP3109389B2 (ja) | 適応フィルタシステム | |
| JP2979712B2 (ja) | フィルタ装置 | |
| JP2961194B2 (ja) | ゴースト除去回路 | |
| JPH0771359B2 (ja) | マルチウエイスピ−カ装置用ネツトワ−ク | |
| JP3414426B2 (ja) | 騒音制御装置 | |
| JPH1051269A (ja) | ローパスフィルタ | |
| JPH04357800A (ja) | 群遅延特性補正装置 | |
| US7262717B2 (en) | Sample rate conversion combined with filter | |
| JP3055560B2 (ja) | デジタルフィルタのフィルタ係数変更方法 | |
| JP4214391B2 (ja) | デジタル・フィルタの設計方法 | |
| JP2629705B2 (ja) | 音質調整装置 | |
| JP2625015B2 (ja) | ゴースト除去装置 | |
| JP3297326B2 (ja) | ディジタルフイルター装置 | |
| JP2932761B2 (ja) | ディジタル信号等化器 | |
| JPH08139601A (ja) | 信号処理回路 | |
| JPH0770957B2 (ja) | Ｆｉｒフィルタ | |
| JPH073935B2 (ja) | 音質調整装置 | |
| JPH0728198B2 (ja) | オクタ−ブ多重フイルタ | |
| JPH0511448B2 (ja) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |