JPS63278410A - 周波数特性設定装置 - Google Patents
周波数特性設定装置Info
- Publication number
- JPS63278410A JPS63278410A JP62113050A JP11305087A JPS63278410A JP S63278410 A JPS63278410 A JP S63278410A JP 62113050 A JP62113050 A JP 62113050A JP 11305087 A JP11305087 A JP 11305087A JP S63278410 A JPS63278410 A JP S63278410A
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- frequency characteristic
- input
- setting device
- sampling
- 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
- 238000005070 sampling Methods 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
産業上の利用分野
本発明は、任意の振幅又は位相又は群遅延周波数特性な
どの信号の周波数特性を実現するトランスバーサル・フ
ィルタ（以下、ＦＩＲフィルタと呼ぶ）を用いたディジ
タルイコライザの周波数特性を入力し、そのフィルタ係
数を求める周波数特性設定装置に関するものである。
どの信号の周波数特性を実現するトランスバーサル・フ
ィルタ（以下、ＦＩＲフィルタと呼ぶ）を用いたディジ
タルイコライザの周波数特性を入力し、そのフィルタ係
数を求める周波数特性設定装置に関するものである。
従来の技術
近年、音響機器のディジタル化にともない、ＦＩＲフィ
ルタを用いたイコライザの開発が求められている。
ルタを用いたイコライザの開発が求められている。
第６図は従来の周波数特性設定装置のブロック図である
。
。
６は振幅周波数特性を入力する入力回路、７は入力回路
６に入力された振幅周波数特性からＦＩＲフィルタのフ
ィルタ係数を求める演算回路、８は実際に入力された振
幅周波数特性を実現するＦＩＲフィルタである。
６に入力された振幅周波数特性からＦＩＲフィルタのフ
ィルタ係数を求める演算回路、８は実際に入力された振
幅周波数特性を実現するＦＩＲフィルタである。
入力回路６で入力される振幅周波数特性は処理される信
号の標本化周波数をＦ、ＦＩＲフィルタのフィルタ係数
の数をＮとすると、周波数分解能がＦ／Ｎで周波数０か
らＦ／２まで入力される。そして、演算回路７において
入力された振幅周波数特性から折り返しく周波数Ｆ ／
２〜Ｆの帯域の振幅周波数特性値）をつくり、この値
を実数環とし、虚数項を零として、逆フーリエ変換を行
うことにより入力された振幅周波数特性を実現するＦＩ
Ｒフィルタ８のフィルタ係数が求められる。そして、こ
のフィルタ係数がＦＩＲフィルタ８に設定されることに
より実際に入力された振幅周波数特性が実現されること
となる。
号の標本化周波数をＦ、ＦＩＲフィルタのフィルタ係数
の数をＮとすると、周波数分解能がＦ／Ｎで周波数０か
らＦ／２まで入力される。そして、演算回路７において
入力された振幅周波数特性から折り返しく周波数Ｆ ／
２〜Ｆの帯域の振幅周波数特性値）をつくり、この値
を実数環とし、虚数項を零として、逆フーリエ変換を行
うことにより入力された振幅周波数特性を実現するＦＩ
Ｒフィルタ８のフィルタ係数が求められる。そして、こ
のフィルタ係数がＦＩＲフィルタ８に設定されることに
より実際に入力された振幅周波数特性が実現されること
となる。
発明が解決しようとする問題点
しかしながら、第６図で示した従来例では処理を行う周
波数帯域と周波数分解能（ＦＩＲフィルタのフィルタ係
数の数）が決まると演算回路７での演算量（特に逆フー
リエ変換の計算時間）が決まりこれを短縮できないとい
う欠点があった。
波数帯域と周波数分解能（ＦＩＲフィルタのフィルタ係
数の数）が決まると演算回路７での演算量（特に逆フー
リエ変換の計算時間）が決まりこれを短縮できないとい
う欠点があった。
本発明は上記問題点を鑑みなされたもので、フィルタ係
数を求める演算量を少なくし、演算時間を大幅に短縮で
きるようにした周波数特性設定装置を提供するものであ
る。
数を求める演算量を少なくし、演算時間を大幅に短縮で
きるようにした周波数特性設定装置を提供するものであ
る。
問題点を解決するための手段
上記問題点を解決するため本発明は、振幅又は位相又は
群遅延などの周波数特性を入力する入力手段と、入力手
段により入力された周波数特性を複数の異なる標本化周
波数に対応した周波数帯域に分割する分割手段と、分割
手段によ多分割されたそたぞれの周波数帯域の周波数特
性を実現するトランスバーサル・フィルタのフィルタ係
数を求める演算手段によって構成される。
群遅延などの周波数特性を入力する入力手段と、入力手
段により入力された周波数特性を複数の異なる標本化周
波数に対応した周波数帯域に分割する分割手段と、分割
手段によ多分割されたそたぞれの周波数帯域の周波数特
性を実現するトランスバーサル・フィルタのフィルタ係
数を求める演算手段によって構成される。
作 用
本発明は、入力手段により振幅又は位相又は群遅延など
の信号の周波数特性が入力され、分割手段によって異な
る複数の標本化周波数に対応する周波数帯域に入力され
た周波数特性を分割し、演算手段によ多分割されたそれ
ぞれの周波数特性を実現するＦＩＲフィルタのフィルタ
係数を求めることにより、フィルタ係数を求める演算時
間を大幅に短縮できることとなる。
の信号の周波数特性が入力され、分割手段によって異な
る複数の標本化周波数に対応する周波数帯域に入力され
た周波数特性を分割し、演算手段によ多分割されたそれ
ぞれの周波数特性を実現するＦＩＲフィルタのフィルタ
係数を求めることにより、フィルタ係数を求める演算時
間を大幅に短縮できることとなる。
実施例
以下、本発明の実施例について図面を参照しながら説明
する。第１図は、本発明の第１の実施例における周波数
特性設定装置のブロック図を示すものである。
する。第１図は、本発明の第１の実施例における周波数
特性設定装置のブロック図を示すものである。
第１図において、１は任意の振幅周波数特性を２つの標
本化周波数に対する周波数帯域に分けて入力する入力回
路、２は入力された振幅周波数特性を２つの異なる標本
化周波数に対応する周波数帯域に分割する分割回路、３
はそれぞれの標本化周波数の周波数帯域に対応するＦＩ
Ｒフィルタのフィルタ係数を求める演算回路、４と６は
標本化周波数が異なりフィルタ係数の数が同じであるＦ
ＩＲフィルタである。
本化周波数に対する周波数帯域に分けて入力する入力回
路、２は入力された振幅周波数特性を２つの異なる標本
化周波数に対応する周波数帯域に分割する分割回路、３
はそれぞれの標本化周波数の周波数帯域に対応するＦＩ
Ｒフィルタのフィルタ係数を求める演算回路、４と６は
標本化周波数が異なりフィルタ係数の数が同じであるＦ
ＩＲフィルタである。
以下、本実施例の動作について説明する。
入力回路１において、任意の振幅周波数特性を２つの標
本化周波数に対する周波数帯域に分け、その周波数帯域
が重なる部分は周波数分解能が小さい方で、周波数帯域
が重ならない周波数帯域では周波数分解能が大きい方で
振幅周波数特性が入力される。
本化周波数に対する周波数帯域に分け、その周波数帯域
が重なる部分は周波数分解能が小さい方で、周波数帯域
が重ならない周波数帯域では周波数分解能が大きい方で
振幅周波数特性が入力される。
第２図ａに入力回路１で入力される振幅周波数特性を示
す。第２図ａにおいて丸（０）印は、設定された振幅値
を示し、ｆｌ、ｆ２は２つの異なる標本化周波数であり
、ｆ２はｆｌの２倍であるとする。
す。第２図ａにおいて丸（０）印は、設定された振幅値
を示し、ｆｌ、ｆ２は２つの異なる標本化周波数であり
、ｆ２はｆｌの２倍であるとする。
入力回路１で入力できる振幅周波数特性の周波数分解能
はこのｆ、／２で部分される。この時、入力回路１によ
って振幅周波数特性を入力できるｆ１／２までの周波数
帯域の周波数分解能はその周波数帯域に対応するＦＩＲ
フィルタ４のフィルタ係数の数をＮとするとｆ１／Ｎと
なり、ｆ１／２からｆ２／２の周波数帯域における入力
できる振幅周波数特性の周波数分解能はｆ２／Ｎとなる
。但し、実際に標本化周波数ｆ１で信号の処理を行うＦ
ＩＲフィルタ４の低域通過フィルタのカットオフ周波数
付近の周波数が（３であった場合は（周波数０からｆ２
／Ｎの分解能で周波数ポイントをとった時に、低域通過
フィルタのカットオフ周波数と一致するかそれに近い周
波数をｆ３とする）、ｆ３からｆ１／２までの設定でき
る振幅周波数特性の周波数分解能は第３図のようにｆ２
／Ｎとする。なお本実施例はこのような入力方法をとる
ものとする。
はこのｆ、／２で部分される。この時、入力回路１によ
って振幅周波数特性を入力できるｆ１／２までの周波数
帯域の周波数分解能はその周波数帯域に対応するＦＩＲ
フィルタ４のフィルタ係数の数をＮとするとｆ１／Ｎと
なり、ｆ１／２からｆ２／２の周波数帯域における入力
できる振幅周波数特性の周波数分解能はｆ２／Ｎとなる
。但し、実際に標本化周波数ｆ１で信号の処理を行うＦ
ＩＲフィルタ４の低域通過フィルタのカットオフ周波数
付近の周波数が（３であった場合は（周波数０からｆ２
／Ｎの分解能で周波数ポイントをとった時に、低域通過
フィルタのカットオフ周波数と一致するかそれに近い周
波数をｆ３とする）、ｆ３からｆ１／２までの設定でき
る振幅周波数特性の周波数分解能は第３図のようにｆ２
／Ｎとする。なお本実施例はこのような入力方法をとる
ものとする。
以上のようにして入力された振幅周波数特性は次の分割
回路２に入力され、２つの標本化周波数及びＦＩＲフィ
ルタ４，６のフィルタ係数の数に対応した周波数特性第
２図す、ｃに分割される（×印、Δ印がその振幅値）。
回路２に入力され、２つの標本化周波数及びＦＩＲフィ
ルタ４，６のフィルタ係数の数に対応した周波数特性第
２図す、ｃに分割される（×印、Δ印がその振幅値）。
以下、この分割回路２における振幅周波数特性の分割方
法について述べる。
法について述べる。
まず、標本化周波数１１ に対応する帯域０〜ｆ１／２
において、０〜ｆ３までの振幅周波数特性は周波数分解
能ｆ１／Ｎで入力回路１において入力されるのでそのま
まこの値を用いる。標本化周波数ｆ、に対応する周波数
帯域で問題となるのはｆ３〜ｆ、／２までの帯域である
が、この帯域は周波数分解能がｆ２／Ｎで振幅周波数特
性が入力されていることからこの振幅周波数特性の値か
らｆ１／Ｎの周波数分解能の新しい振幅周波数ポイント
を求める必要がある。この新しい振幅周波数ポイントの
求めかたにはいろいろな方法が考えられるが、ここでは
線形補間して求める方法を示す。つまり入力回路１にお
いて入力されたｆ３〜１１／２の振幅周波数の値のうち
１３での値をＡ２次の値をＢとするとＡ、８間にある本
来の標本化周波数ｆ１ での値Ｃは、標本化周波数ｆ２
がｆｌの２倍であるので Ｃ＝（Ｂ Ａ ） ｘ ｆ １／ ｆ ２ ＋ Ａと
して求めることができる。ｆ１／２までの他の値は同様
にして求めることができる。なお、ここでは線形補間を
おこなって標本化周波数ｆ１に対応した周波数帯域の振
幅周波数特性を求めたが、先に述べたように他の方法、
例えば多次関数による方法などいろいろな方法が考えら
れる。
において、０〜ｆ３までの振幅周波数特性は周波数分解
能ｆ１／Ｎで入力回路１において入力されるのでそのま
まこの値を用いる。標本化周波数ｆ、に対応する周波数
帯域で問題となるのはｆ３〜ｆ、／２までの帯域である
が、この帯域は周波数分解能がｆ２／Ｎで振幅周波数特
性が入力されていることからこの振幅周波数特性の値か
らｆ１／Ｎの周波数分解能の新しい振幅周波数ポイント
を求める必要がある。この新しい振幅周波数ポイントの
求めかたにはいろいろな方法が考えられるが、ここでは
線形補間して求める方法を示す。つまり入力回路１にお
いて入力されたｆ３〜１１／２の振幅周波数の値のうち
１３での値をＡ２次の値をＢとするとＡ、８間にある本
来の標本化周波数ｆ１ での値Ｃは、標本化周波数ｆ２
がｆｌの２倍であるので Ｃ＝（Ｂ Ａ ） ｘ ｆ １／ ｆ ２ ＋ Ａと
して求めることができる。ｆ１／２までの他の値は同様
にして求めることができる。なお、ここでは線形補間を
おこなって標本化周波数ｆ１に対応した周波数帯域の振
幅周波数特性を求めたが、先に述べたように他の方法、
例えば多次関数による方法などいろいろな方法が考えら
れる。
また、周波数ｆ４（＞ｆ３）からｆ１／２においては、
ｆ１／２の振幅周波数特性値が零となるように補正を行
った値を求め、これを使用する方法も考えられる。この
方法を示したのが第４図である。このｆ４の値は先に述
べた低域通過フィルタの信号抑制レベルが十分な値（例
えば通過帯域のレベルから１０ｄＢ落ちた周波数）にな
る周波数とする。
ｆ１／２の振幅周波数特性値が零となるように補正を行
った値を求め、これを使用する方法も考えられる。この
方法を示したのが第４図である。このｆ４の値は先に述
べた低域通過フィルタの信号抑制レベルが十分な値（例
えば通過帯域のレベルから１０ｄＢ落ちた周波数）にな
る周波数とする。
周波数１４における入力回路１により入力された振幅値
をＤとすると新しい標本化周波数ｆ１ での値Ｅは Ｅ＝−Ｄ／（ｆ１／２−ｆ、）ｘｆ１／Ｎ＋Ｅとして求
めることができ、多の値も同様にして求めることができ
る。このようにすることにより、先に述べた低域通過フ
ィルタの阻止イｉｉ域の信号を更に抑制することが可能
となる。また、ｆ１／２において振幅周波数特性を零と
することは、演算回路３においてフィルタ係数をリニア
７エーズの条件で求める際に必要な条件でもある。
をＤとすると新しい標本化周波数ｆ１ での値Ｅは Ｅ＝−Ｄ／（ｆ１／２−ｆ、）ｘｆ１／Ｎ＋Ｅとして求
めることができ、多の値も同様にして求めることができ
る。このようにすることにより、先に述べた低域通過フ
ィルタの阻止イｉｉ域の信号を更に抑制することが可能
となる。また、ｆ１／２において振幅周波数特性を零と
することは、演算回路３においてフィルタ係数をリニア
７エーズの条件で求める際に必要な条件でもある。
また同様に標本化周波数１２の帯域も、処理する信号に
影響しない周波数１６から同様な補正を行なうと同様の
効果が得られる。
影響しない周波数１６から同様な補正を行なうと同様の
効果が得られる。
次に、標本化周波数１２に対応する周波数帯域０−ｆ２
／２における振幅周波数特性について述べる。
／２における振幅周波数特性について述べる。
周波数ｆ３からｆ２／２までの値は入力回路１において
周波数分解能ｆ２／Ｎで入力されるためこの値をそのま
ま使用できる。周波数０からｆ３までの値が周波数分解
能がちがうため入力振幅周波数特性をそのまま使用でき
ないため演算を行ってこの帯域の値を求める必要がある
。この演算の方法も色々と考えられるが、ここでは入力
された振幅周波数特性を間引く方法を示す。つまり、こ
こではｆ がｆ の２倍であることから周波数０〜ｆ３
まｒｍ（ヒ周波数ｆ２、フィルタ係数Ｎでの振幅周波数
特性は、周波数分解能ｆ、／Ｎで入力された振幅周波数
特性を２個に１個間引くことにより求めることができる
わけである（第２図Ｃを参照）このように間引きを行う
ことにより入力回路１で周波数分解能１１／Ｎで入力さ
れた振幅周波数特性から標本化周波数１２に対応する周
波数帯域の振幅周波数特性を求めることができる。
周波数分解能ｆ２／Ｎで入力されるためこの値をそのま
ま使用できる。周波数０からｆ３までの値が周波数分解
能がちがうため入力振幅周波数特性をそのまま使用でき
ないため演算を行ってこの帯域の値を求める必要がある
。この演算の方法も色々と考えられるが、ここでは入力
された振幅周波数特性を間引く方法を示す。つまり、こ
こではｆ がｆ の２倍であることから周波数０〜ｆ３
まｒｍ（ヒ周波数ｆ２、フィルタ係数Ｎでの振幅周波数
特性は、周波数分解能ｆ、／Ｎで入力された振幅周波数
特性を２個に１個間引くことにより求めることができる
わけである（第２図Ｃを参照）このように間引きを行う
ことにより入力回路１で周波数分解能１１／Ｎで入力さ
れた振幅周波数特性から標本化周波数１２に対応する周
波数帯域の振幅周波数特性を求めることができる。
なお、ここでは間引きを行って標本化周波数ｆ２に対応
した周波数帯域の振幅周波数特性を求めた力ζ先に述べ
たように他の方法、例えば０〜ｆ６（＜ｆ３）の振幅周
波数特性を零にする、または、ＯｄＢ （＝１．０ ）
にする方法などが考えられる。
した周波数帯域の振幅周波数特性を求めた力ζ先に述べ
たように他の方法、例えば０〜ｆ６（＜ｆ３）の振幅周
波数特性を零にする、または、ＯｄＢ （＝１．０ ）
にする方法などが考えられる。
次に演算回路３において分割回路２で分割された２つの
周波数特性を使用してＦＩＲフィルタ４゜６のフィルタ
係数が求められる。ここでの演算は、それぞれの周波数
帯域の折り返しくｆ １／２〜ｆ１゜ｆ２／２〜ｆ２の
周波数帯域）の振幅周波数特性を求め、リニアフェーズ
の条件から振幅周波数特性を実数項とし、虚数項を零と
した逆フーリエ変換を行いフィルタ係数を求める方法な
どがとられる。
周波数特性を使用してＦＩＲフィルタ４゜６のフィルタ
係数が求められる。ここでの演算は、それぞれの周波数
帯域の折り返しくｆ １／２〜ｆ１゜ｆ２／２〜ｆ２の
周波数帯域）の振幅周波数特性を求め、リニアフェーズ
の条件から振幅周波数特性を実数項とし、虚数項を零と
した逆フーリエ変換を行いフィルタ係数を求める方法な
どがとられる。
また、虚数項を零とはせずに、ある位相周波数特性ある
いは群遅延特性を示す虚数項を与え、逆フーリエ変換を
行う方法なども考えられる。またしベルト変換の関係か
ら最小位相推移を行なう、ようなフィルタ係数演算法な
どもとられる。このようＫすることにより、それぞれの
標本化周波数に対応する周波数帯域に対して、入力され
た振幅周波数特性（または、振幅周波数特性と位相周波
数特性又は群遅延特性）を実現するフィルタ係数が求め
られる。このようにして求められたフィルタ係数は、そ
れぞれの周波数帯域に対応したＦＩＲフィルタ４，５に
入力され、実際に入力回路１に入力された振幅周波数特
性が実現されることとなる。
いは群遅延特性を示す虚数項を与え、逆フーリエ変換を
行う方法なども考えられる。またしベルト変換の関係か
ら最小位相推移を行なう、ようなフィルタ係数演算法な
どもとられる。このようＫすることにより、それぞれの
標本化周波数に対応する周波数帯域に対して、入力され
た振幅周波数特性（または、振幅周波数特性と位相周波
数特性又は群遅延特性）を実現するフィルタ係数が求め
られる。このようにして求められたフィルタ係数は、そ
れぞれの周波数帯域に対応したＦＩＲフィルタ４，５に
入力され、実際に入力回路１に入力された振幅周波数特
性が実現されることとなる。
なお、実施例では入力回路１において振幅周波数特性を
第１種（周波数０からｆ１／Ｎ毎に、１３からはｆ２７
α毎に周波数ポイントをとる）の標本化を行った場合を
示したが、第２種（周波数ｆ ／ＮＸ１／２からｆ／Ｎ
毎に、ｆ３（この場合、周波数ｆ２／Ｎ×１／２からｆ
２／Ｎ毎に周波数ポイントをとっていったときに低域通
過フィルタのカットオフ周波数に一致するかそれに近い
値をｆ３とする）からはｆ２β毎に周波数ポイントをと
る）の標本化を行った場合も、周波数ポイントが周波数
帯域の違いでずれが生じる場合があるが、直線補間又は
多次関数近似などの方法をとることにより同様の考え方
で周波数分割を行うことができ、それぞれの周波数帯域
のフィルタ係数を求めることができる。
第１種（周波数０からｆ１／Ｎ毎に、１３からはｆ２７
α毎に周波数ポイントをとる）の標本化を行った場合を
示したが、第２種（周波数ｆ ／ＮＸ１／２からｆ／Ｎ
毎に、ｆ３（この場合、周波数ｆ２／Ｎ×１／２からｆ
２／Ｎ毎に周波数ポイントをとっていったときに低域通
過フィルタのカットオフ周波数に一致するかそれに近い
値をｆ３とする）からはｆ２β毎に周波数ポイントをと
る）の標本化を行った場合も、周波数ポイントが周波数
帯域の違いでずれが生じる場合があるが、直線補間又は
多次関数近似などの方法をとることにより同様の考え方
で周波数分割を行うことができ、それぞれの周波数帯域
のフィルタ係数を求めることができる。
また、実施例では振幅周波数特性を入力する場。
合を示したが入力を位相周波数特性又は群遅延特性など
の他の周波数特性とした場合も同様の考え方で行うこと
ができる。
の他の周波数特性とした場合も同様の考え方で行うこと
ができる。
第５図は、本発明の第２の実施例における入力回路の振
幅周波数特性の入力方法及び分割回路の周波数分割の方
法を示す図である。第２の実施例の構成は第１図と同様
であシ、入力方法１分割方法のみが異なる。
幅周波数特性の入力方法及び分割回路の周波数分割の方
法を示す図である。第２の実施例の構成は第１図と同様
であシ、入力方法１分割方法のみが異なる。
入力回路１で入力できる振幅周波数特性（第５図ａ）は
、全帯域にわたって処理できる標本化周波数ｆ１． ｆ
２とは異なった（この図の場合は、ｆ１／２Ｎの）周波
数分解能で設定できるようになっている。次の分割回路
２では、実際の信号を処理できる周波数ポイントを得る
ために、入力された周波数ポイントから、１１の標本化
周波数に対応する周波数帯域では、入力された周波数ポ
イントを２個のうち１個を間引く演算を（第６図ｂ）、
ｆ２の標本化周波数に対応する周波数帯域では、入力さ
れた周波数ポイントを４個のうち１個を間引く演算を行
ない（第６図Ｃ）、入力された振幅周波数特性を分割す
る。
、全帯域にわたって処理できる標本化周波数ｆ１． ｆ
２とは異なった（この図の場合は、ｆ１／２Ｎの）周波
数分解能で設定できるようになっている。次の分割回路
２では、実際の信号を処理できる周波数ポイントを得る
ために、入力された周波数ポイントから、１１の標本化
周波数に対応する周波数帯域では、入力された周波数ポ
イントを２個のうち１個を間引く演算を（第６図ｂ）、
ｆ２の標本化周波数に対応する周波数帯域では、入力さ
れた周波数ポイントを４個のうち１個を間引く演算を行
ない（第６図Ｃ）、入力された振幅周波数特性を分割す
る。
以下は第１の実施例と同様の処理を行なう。
本実施例のように、入力回路１においては、処理できる
周波数分解能とは異なった分解能で周波数特性を設定で
きるようにしても、分割回路２において適当な演算を行
なうことにより、実際に処理できる周波数特性を得るこ
とができ、これによｐＦＩＲフィルタ４，６のフィルタ
係数も求める ・ことが可能となる。
周波数分解能とは異なった分解能で周波数特性を設定で
きるようにしても、分割回路２において適当な演算を行
なうことにより、実際に処理できる周波数特性を得るこ
とができ、これによｐＦＩＲフィルタ４，６のフィルタ
係数も求める ・ことが可能となる。
なお、本実施例では、分割回路２では第１種の標本化を
行なった場合を示したが、第２種の標本化を行なう場合
も考えられる。
行なった場合を示したが、第２種の標本化を行なう場合
も考えられる。
また、第１．第２の実施例において、第１種の標本化を
行なった場合は、通常、演算回路３において求めたフィ
ルタ係数に窓をかける必要がある。
行なった場合は、通常、演算回路３において求めたフィ
ルタ係数に窓をかける必要がある。
したがって、第１図の演算回路３において、何らかの特
性をもつ窓関数をかける演算を行なう機能をもたせる構
成も考えられる。ただし、第２種の標本化を行なった場
合は特に窓をかける必要はガいことを明記しておく。
性をもつ窓関数をかける演算を行なう機能をもたせる構
成も考えられる。ただし、第２種の標本化を行なった場
合は特に窓をかける必要はガいことを明記しておく。
発明の効果
以上述べたように、本発明は入力周波数特性を複数標本
化周波数の周波数帯域に分割する構成をとることにより
、それぞれの周波数帯域における処理するデータ数を減
少させることができ、これにより、フィルタ係数を求め
る演算時間を大幅に短縮できることどなる。
化周波数の周波数帯域に分割する構成をとることにより
、それぞれの周波数帯域における処理するデータ数を減
少させることができ、これにより、フィルタ係数を求め
る演算時間を大幅に短縮できることどなる。
第１図は本発明の第１の実施例における周波数特性設定
装置のブロック図、第２図は本発明の第１の実施例にお
ける入力周波数特性及び分割周波数特性を示す特性図、
第３図は本発明の第１の実施例における他の入力周波数
特性を示す特性図、第４図は本発明の第１の実施例にお
ける分割周波数特性を補正する方法を示す特性図、第６
図は本発明の第２の実施例における入力周波数特性及び
分割周波数特性を示す特性図、第６図は従来例を示すブ
ロック図である。 １・・・・・・入力回路、２・・・・・・分割回路、３
・・・・・・演算口ＭＬ ４１５・・・・・・ＦＩＲフ
ィルタ。 代理人の氏名 弁理士 中 尾 敏 男 ほか１名第２
図 第３図 眉城杖
装置のブロック図、第２図は本発明の第１の実施例にお
ける入力周波数特性及び分割周波数特性を示す特性図、
第３図は本発明の第１の実施例における他の入力周波数
特性を示す特性図、第４図は本発明の第１の実施例にお
ける分割周波数特性を補正する方法を示す特性図、第６
図は本発明の第２の実施例における入力周波数特性及び
分割周波数特性を示す特性図、第６図は従来例を示すブ
ロック図である。 １・・・・・・入力回路、２・・・・・・分割回路、３
・・・・・・演算口ＭＬ ４１５・・・・・・ＦＩＲフ
ィルタ。 代理人の氏名 弁理士 中 尾 敏 男 ほか１名第２
図 第３図 眉城杖
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 |
| DE3853372T DE3853372T2 (de) | 1987-03-23 | 1988-01-15 | Filterkoeffizientenberechnung für ein digitales 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 |
| CA000556581A CA1280808C (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 |
| 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 true JPS63278410A (ja) | 1988-11-16 |
| JPH0770956B2 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 (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 | デジタル・グラフイツク・イコライザ |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1280808C (en) | 1987-03-23 | 1991-02-26 | Seiichi Ishikawa | Calculation of filter factors for digital filter |
-
1987
- 1987-05-08 JP JP11305087A patent/JPH0770956B2/ja not_active Expired - Fee Related
Patent Citations (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 | デジタル・グラフイツク・イコライザ |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0770956B2 (ja) | 1995-07-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100462447B1 (ko) | M-비트디지털신호의'n'개의최하위비트를대칭적으로단축시키는장치 | |
| KR100374395B1 (ko) | 디지탈화된신호용등화장치 | |
| US20160037252A1 (en) | Audio filters utilizing sine functions | |
| JPS63278410A (ja) | 周波数特性設定装置 | |
| US5633689A (en) | Apparatus for separating a digital composite video signal into components | |
| CN100391102C (zh) | 音质调整装置以及滤波器装置 | |
| JP2001119270A (ja) | 帯域幅切替型デジタル・フィルター | |
| US10536136B2 (en) | Modified digital filtering with sample zoning | |
| Noceti Filho et al. | Filter design with simultaneous magnitude and phase specifications | |
| JPH05145377A (ja) | 周波数帯域分割フイルタ演算方式 | |
| JPH073935B2 (ja) | 音質調整装置 | |
| JP2979712B2 (ja) | フィルタ装置 | |
| JPH0770951B2 (ja) | オクタ−ブ多重フイルタ | |
| Korohoda et al. | A Practical Multirate FIR-Based Equalizer-A Case Study | |
| JPS6326112A (ja) | トランスバ−サルフイルタ−係数演算装置 | |
| JP2001326991A (ja) | 音声処理装置 | |
| JPH0310469A (ja) | ゴースト除去装置 | |
| JPH0973295A (ja) | 能動制御装置 | |
| JP2629705B2 (ja) | 音質調整装置 | |
| JPS62109421A (ja) | デイジタル・チヤンネル・デバイダ | |
| JPH0738565B2 (ja) | 高速デイジタルフイルタ | |
| CN117784109A (zh) | 一种基于fpga的图像声纳可重构信号预处理方法 | |
| JPS63126311A (ja) | 音質調整装置 | |
| JPH0621983A (ja) | Ｆｓｋ復調回路 | |
| JPH08336157A (ja) | コムフィルター |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |