JPH11145839A

JPH11145839A - ディジタルフィルタの量子化係数の決定方法及びその決定プログラムが記録されたコンピュータ上で読み出し可能な記録媒体

Info

Publication number: JPH11145839A
Application number: JP30533597A
Authority: JP
Inventors: Osamu Taniguchi; 修谷口
Original assignee: NEC Corp
Current assignee: NEC Corp
Priority date: 1997-11-07
Filing date: 1997-11-07
Publication date: 1999-05-28
Anticipated expiration: 2017-11-07
Also published as: JP3235546B2

Abstract

(57)【要約】【課題】所定の係数語長でリップル幅の小さいディジ
タルフィルタの量子化係数の決定方法を提供する。【解決手段】ディジタルフィルタの量子化係数を所定
の係数語長にし、且つ、所定の周波数特性のディジタル
フィルタを得るための、ディジタルフィルタの量子化係
数の決定方法において、最適な量子化係数を得るため
に、得られたディジタルフィルタの周波数特性の、予め
定めたリップル幅を越えた量を評価すると共に、前記量
子化係数を変化させながら前記評価を繰り返し実行する
ことで、最適な量子化係数を決定することを特徴とする
ディジタルフィルタの量子化係数の決定方法。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、ディジタルフィル
タの量子化係数の決定方法及びその決定プログラムが記
録されたコンピュータ上で読み出し可能な記録媒体に係
り、特に、目的仕様のリップル幅に対し、ペナルティ関
数と誤差関数とを組み合わせた評価関数を計算すること
により、量子化誤差によるディジタルフィルタ特性の劣
化を減少させ、所定の係数語長で且つ、所定のリップル
内に収めるようにしたディジタルフィルタの量子化係数
の決定方法に関する。

【０００２】

【従来の技術】非巡回形ディジタルフィルタは、動作が
常に安定しており、完全な直線位相特性が可能であり、
Ａ／Ｄコンバータ等に広く使用されている。一方で、急
峻なフィルタ特性を実現するためには、次数が高くなる
だけではなく、係数の量子化誤差を抑えるため係数語長
を大きくとる必要があり、ハードウェア規模が膨大にな
る欠点がある。このため、係数の量子化は、例えば、レ
メツ法で求めた係数を、希望のビット精度で四捨五入す
ることによって行われている。しかし、単に四捨五入し
たのみでは、理想特性からのずれが大きく目的のリップ
ル特性を得ることが困難となる。

【０００３】例えば、８３タップの非巡回形ディジタル
フィルタの係数をレメツ法で算出し、その周波数特性を
示したのが図２である。図に見られるように係数に量子
化の制限が無ければ例えば、仕様が±０．０４ｄＢのリ
ップル幅をとることも問題ではない。しかし、例えば１
０ビットの量子化を行うと、図４のように±０．０４ｄ
Ｂのリップル幅の仕様を満足しない。

【０００４】このようにディジタルフィルタの係数の量
子化誤差による特性の劣化を防ぐため、例えば、特開昭
５８−４７３２２号公報に記載された技術が用いられ
る。この技術は、量子化された係数を用いたフィルタ
に、整数係数を用いたフィルタを従属接続し、前段のフ
ィルタの下位ビットを変化させて最小の誤差評価関数を
求めるものである。詳細に述べると、図６に見られるよ
うに第１のディジタルフィルタ１１及び第２のディジタ
ルフィルタ１２の従属接続によって構成されている。

【０００５】第１のディジタルフィルタ１１の伝達関数
をＦ（ｚ）、第２のディジタルフィルタ１２の伝達関数
をＷ（ｚ）とする。第１のディジタルフィルタ１１は従
来と同様に有限語長に量子化されている。しかし単に四
捨五入ではなく、後述する誤差評価関数を最小化するよ
うにして求められている。第２のフィルタ１２は阻止帯
域に有限個の零点を持つ非巡回形ディジタルフィルタで
あり、その係数は整数である。例えば、図６のように構
成されていて、その伝達係数Ｗ（ｚ）はＷ（ｚ）＝（１＋２ｚ^-1＋ｚ^-2）（１＋ｚ^-1＋ｚ^-2）で正確に与えられ誤差はない。この第１、第２のディジ
タルフィルタを縦列接続した図６に示すディジタルフィ
ルタの伝達係数Ｈ（ｚ）は、Ｈ（ｚ）＝Ｆ（ｚ）・Ｗ（ｚ）で与えられる。従って、第１のフィルタの係数を量子化
することによって生ずる伝達関数Ｆ（ｚ）の誤差をΔＦ
（ｚ）とすると、伝達関数Ｈ（ｚ）の誤差ΔＨ（ｚ）
は、 ΔＨ（ｚ）＝ΔＦ（ｚ）・Ｗ（ｚ）で与えられる。そこで、下記式によって与えられる誤差
評価関数

【０００６】

【数１】

【０００７】即ち、誤差関数εが最小となるように第１
のフィルタ１１の各係数を設定し、第１、第２のディジ
タルフィルタ１１、１２を縦列接続することで所定の特
性のディジタルフィルタを実現している。しかし、上記
した従来のものは、誤差評価関数の結果を小さくするの
みであり、目的仕様のリップル幅に対しては何ら考慮し
ていないから、リップル幅の小さい特性のディジタルフ
ィルタが得られないという欠点があった。

【０００８】

【発明が解決しようとする課題】本発明は、上記した従
来技術の欠点を改良し、別のフィルタ回路を追加するこ
となくリップル幅を小さく、しかも、係数語長が短い、
新規なディジタルフィルタの量子化係数の決定方法、及
びその決定プログラムが記録されたコンピュータ上で読
み出し可能な記録媒体を提供することを目的とする。

【０００９】

【課題を解決するための手段】本発明は、上記した目的
を達成するため、基本的には、以下に記載されたような
技術構成を採用するものである。即ち、本発明に係わる
ディジタルフィルタの量子化係数の決定方法の第１の態
様としては、ディジタルフィルタの量子化係数を所定の
係数語長にし、且つ、所定の周波数特性のディジタルフ
ィルタを得るための、ディジタルフィルタの量子化係数
の決定方法において、最適な量子化係数を得るために、
得られたディジタルフィルタの周波数特性の、予め定め
たリップル幅を越えた量を評価すると共に、前記量子化
係数を変化させながら前記評価を繰り返し実行すること
で、最適な量子化係数を決定するものであり、第２の態
様としては、前記リップル幅を越えた量と、理想周波数
特性と前記得られたディジタルフィルタの周波数特性と
の差とを合計した評価値に基づき評価値を最小にするこ
とで、最適な量子化係数を決定するものであり、第３の
態様としては、前記評価値を得る際、前記リップル幅を
越えた量には重みづけを行い演算するものであり、第４
の態様としては、ディジタルフィルタの量子化係数の決
定方法において、（１）略理想的な量子化係数を求め、
所定の係数語長にするために四捨五入して量子化係数を
得る第１ステップと、（２）理想周波数特性と前記第１
ステップで得られた量子化係数によるフィルタの周波数
特性との第１の誤差を所定周波数にわたって求める第２
ステップと、（３）前記第１ステップで得られた量子化
係数によるフィルタの周波数特性のうち予め定めたリッ
プル幅を越えた部分である第２の誤差を所定の周波数に
わたって求める第３ステップと、（４）前記第１の誤差
と第２の誤差とを合計して量子化係数の最適値を評価す
るための評価値を求める第４ステップと、（５）前記第
４ステップで求められた評価値を、既に求められている
評価値とを比較し、その結果、前記第４ステップで求め
られた評価値が、既に求められている評価値よりよい場
合、その量子化係数を保存する第５ステップとからな
り、量子化係数の所定のビットを変化させながら前記第
１ステップから第５ステップまでを繰り返し実行し、最
適な量子化係数を決定するものであり、第５の態様とし
ては、前記第４ステップにおいて、評価値を求める際、
第２の誤差に重みづけをした後第１の誤差と合計して評
価値を求めるものである。

【００１０】又、本発明に係わるディジタルフィルタの
量子化係数の決定プログラムがコンピュータ上で読み出
し可能な記録媒体は、第１の態様としては、シーケンサ
にコントロールされてＲＪＭに格納された量子化係数と
ＲＡＭに格納された畳み込みの演算結果とを取り出して
所定の乗算を行い、この結果をアキュムレータに格納す
ると共に、畳み込みの結果を前記ＲＡＭに格納するハー
ドウェアを含み、ディジタルフィルタの量子化係数を所
定の係数語長にし、且つ、所定の周波数特性のディジタ
ルフィルタを得るための、ディジタルフィルタの量子化
係数の決定方法において、最適な量子化係数を得るため
に、得られたディジタルフィルタの周波数特性の、予め
定めたリップル幅を越えた量を評価すると共に前記量子
化係数を変化させながら前記評価を繰り返し実行するこ
とで、最適な量子化係数を決定するものであり、第２の
態様としては、前記リップル幅を越えた量と、理想周波
数特性と前記得られたディジタルフィルタの周波数特性
との差とを合計した評価値に基づき評価値を最小にする
ことで、最適な量子化係数を決定するものであり、第３
の態様としては、前記評価値を得る際、前記リップル幅
を越えた量には重みづけを行い演算するものであり、第
４の態様としては、シーケンサにコントロールされてＲ
ＯＭに格納された量子化係数とＲＡＭに格納された畳み
込みの演算結果とを取り出して所定の乗算を乗算器で行
い、この結果をアキュムレータに格納すると共に、畳み
込みの演算結果を前記ＲＡＭに格納するディジタルフィ
ルタの前記量子化係数の決定方法において、（１）略理
想的な量子化係数を求め、所定の係数語長にするために
四捨五入して量子化係数を得る第１ステップと、（２）
理想周波数特性を前記第１ステップで得られた量子化係
数によるフィルタの周波数特性との第１の誤差を所定周
波数にわたって求める第２ステップと、（３）前記第１
ステップで得られた量子化係数によるフィルタの周波数
特性のうち予め定めたリップル幅を越えた部分である第
２の誤差を所定の周波数にわたって求める第３ステップ
と、（４）前記第１の誤差と第２の誤差とを合計して量
子化係数の最適値を評価するための評価値を求める第４
ステップと、（５）前記第４ステップで求められた評価
値と、既に求められている評価値とを比較し、その結
果、前記第４ステップで求められた評価値が、既に求め
られている評価値よりよい場合、その量子化係数を保存
する第５ステップとからなり、量子化係数の所定のビッ
トを変化させながら前記第１ステップから第５ステップ
まで繰り返し実行し、最適な量子化係数を決定するもの
であり、第５の態様としては、前記第４ステップ１にお
いて、評価値を求める際、第２の誤差に重みづけをした
後、第１の誤差と合計して評価値を求めるものである。

【００１１】

【発明の実施の形態】本発明に係るディジタルフィルタ
の量子化係数の決定方法は、（１）略理想的な量子化係
数を求め、所定の係数語長にするために四捨五入して量
子化係数を得る第１ステップと、（２）理想周波数特性
と前記第１ステップで得られた量子化係数によるフィル
タの周波数特性との第１の誤差を所定周波数にわたって
求める第２ステップと、（３）前記第１ステップで得ら
れた量子化係数によるフィルタの周波数特性のうち予め
定めたリップル幅を越えた部分である第２の誤差を所定
の周波数にわたって求める第３ステップと、（４）前記
第１の誤差と第２の誤差とを合計して量子化係数の最適
値を評価するための評価値を求める第４ステップと、
（５）前記第４ステップで求められた評価値と、既に求
められている評価値とを比較し、その結果、前記第４ス
テップで求められた評価値が、既に求められている評価
値よりよい場合、その量子化係数を保存する第５ステッ
プとからなり、量子化係数の所定のビットを変化させな
がら前記第１ステップから第５ステップまでを繰り返し
実行し、最適な量子化係数を決定するものであるから、
別のフィルタ回路を追加することなく、リップル幅の小
さいディジタルフィルタが得られる。

【００１２】しかも、係数語長は所定のビット数で短い
から、回路規模は大きくならない。

【００１３】

【実施例】以下に、本発明に係るディジタルフィルタの
量子化係数の決定方法の具体例を図面を参照しながら詳
細に説明する。図１乃至図５は本発明の具体例を示す図
であり、図には、ディジタルフィルタの量子化係数を所
定の係数語長にし、且つ、所定の周波数特性のディジタ
ルフィルタを得るための、ディジタルフィルタの量子化
係数の決定方法において、最適な量子化係数を得るため
に、得られたディジタルフィルタの周波数特性の、予め
定めたリップル幅を越えた量を評価する（Ｓ１０３，１
０４，１０５）と共に、前記量子化係数を変化させなが
ら前記評価を繰り返し実行することで、最適な量子化係
数を決定することを特徴とするディジタルフィルタの量
子化係数の決定方法が示され、又、前記リップル幅を越
えた量と、理想周波数特性と前記得られたディジタルフ
ィルタの周波数特性との差とを合計した評価値に基づき
評価値を最小にする（ステップＳ１０４、Ｓ１０５）こ
とで、最適な量子化係数を決定するディジタルフィルタ
の量子化係数の決定方法が示され、更に、前記評価値を
得る際、前記リップル幅を越えた量には重みづけを行い
演算するディジタルフィルタの量子化係数の決定方法が
示されている。

【００１４】図１を参照すると、本実施例は、非巡回形
ディジタルフィルタの係数をレメツ法で求め、その係数
を量子化して四捨五入する計算部Ｓ１０１と、誤差関数
Ｅ（ｆ，｛ｈ｝）を算出する計算部Ｓ１０２と、ペナル
ティ関数Ｐ（ｆ，｛ｈ｝）を算出する計算部Ｓ１０３
と、ペナルティ関数に重みＷを乗じたものと誤差関数と
を加算したものを評価関数として計算する計算部Ｓ１０
４と、この評価関数の値が既に得られている結果と比較
し小さければ、その状態と評価関数の値とを保持する保
持部Ｓ１０５と、すべての係数パターンが終了したかを
判別する判定部１０６と、各係数の下位ビットを１ビッ
トだけ修正させる修正部Ｓ１０７と、結果を出力する出
力部Ｓ１０８から構成される。

【００１５】誤差関数Ｅ（ｆ，｛ｈ｝）を算出する計算
部Ｓ１０２の計算手法について以下に詳細に述べる。図
３を参照すると、レメツ法で求めた係数を四捨五入して
得られた量子化係数での非巡回ディジタルフィルタの周
波数特性をＲ（ｆ，｛ｈ｝）とし、ｆは周波数、ｈは量
子化された係数、｛ｈ｝は量子化された係数の一つの組
み合わせを示す。またＨ（ｆ）はリプル幅の無い理想の
周波数特性を示すとする。誤差関数（第１の誤差）を計
算するということは、

【００１６】

【数２】

【００１７】を計算するということである。ここでｆｅ
ｎｄはリプル帯域の周波数である。誤差関数は理想関数
Ｈ（ｆ）から量子化された関数Ｒ（ｆ，｛ｈ｝）の誤差
の全体の累積の量をあらわしたものである。次にペナル
ティ関数Ｐ（ｆ，｛ｈ｝）（第２の誤差）を算出する計
算部Ｓ１０３について説明する。図４を参照すると、ペ
ナルティ関数は目的仕様の値とするリプル幅δからのず
れの量をリプル帯域内にわたって積分された量を意味す
る。即ち、図３の斜線部に対応する。ペナルティ関数は
以下のように示される。

【００１８】

【数３】

【００１９】ただし、関数Ｇ（ｆ，｛ｈ｝）は

【００２０】

【数４】

【００２１】である。次に、評価関数Ｑ（｛ｈ｝）を算
出する計算部Ｓ１０４について説明する。評価関数はペ
ナルティ関数（２）式に重みＷを乗じたものと誤差関数
（１）式を加算したものである。即ち、以下の（３）式
になる。

【００２２】

【数５】

【００２３】また、保持部Ｓ１０５は既に得られた過去
の評価関数の結果を保持し、現在の結果と比較し評価関
数が小さければ保持する部分である。また、判定部Ｓ１
０６は予定した係数の組み合わせの全パターンが終了し
たかを判別する部分であり、本発明では組み合わせの大
きさを考慮して下位１ビットのみの変化の組み合わせと
した。また、修正部Ｓ１０７は量子化された係数のなか
で一つの係数に着目し、下位ビットのなかの一つのビッ
トを修正させる部分である。例えば、着目した係数が１
０１１０１１１０１の場合下位ビットの１つを変化させ
１０１１０１１１００と修正することである。

【００２４】次に、本発明を更に、具体的に説明する
と、まず、計算部Ｓ１０１において、非巡回形ディジタ
ルフィルタの係数をレメツ法で求め、その係数を所定の
係数語長で量子化して四捨五入を行う。この場合、量子
化される前の係数では図２に見られるように、リプル幅
は十分目的の仕様を満たしているが、係数が量子化され
ると図４に見られるように、特性が悪化し、目的の仕様
であるリプル幅±δを超えてしまう。次に、この状態で
の誤差関数を求める（Ｓ１０２）。この場合、（１）式
を使用して計算するが、実際は周波数特性の周波数の幅
δｆを適当にとり図３で見られるように、量子化された
係数での非巡回ディジタルフィルタの周波数特性Ｒ
（ｆ，｛ｈ｝）とリプル幅の無い理想の周波数特性を持
つＨ（ｆ）との差をリプル帯域内で総和をとることにな
る。式にすると、

【００２５】

【数６】

【００２６】となる。ペナルティ関数を算出する計算部
（Ｓ１０３）についても、同様に計算を行い

【００２７】

【数７】

【００２８】となる。次に評価関数Ｑ（｛ｈ｝）を算出
する（Ｓ１０４）。次に、保持部（Ｓ１０５）で評価関
数の大きさの評価及び保持を行う。次に、評価関数の、
最小値を求めるためには修正部Ｓ１０７の量子化された
係数を修正する項ですべてのビットパターンについて係
数を修正させ、すべての係数のビットパターンに対しそ
れぞれ評価関数を算出して最小値を求めることが望まし
い。しかし、すべてのビットパターンについて行うと、
１０ビットの係数で構成され、８３タップのハーフバン
ド非巡回ディジタルフィルタの場合、１０２４²¹＝１．
６５×１０⁶³であり非現実的である。

【００２９】そこで、レメツ法を四捨五入することで量
子化した係数は、最適係数の組み合わせに対してそれほ
ど変化しなく、フィルタの全体の特性もそれほど変わら
ないと考え下位１ビットのみを変化させ、すべてのパタ
ーンについてシミュレーションを行うと、２²¹＝２．１
０×１０⁶パターンであり実行にそれほど問題はない。

【００３０】図５は本発明を適用するためのハードウェ
アを示し、ハードウェア１はシーケンサ２にコントロー
ルされて、ＲＯＭ３に格納された量子化係数と、ＲＡＭ
４に格納された畳み込みの演算結果とを取り出して所定
の乗算を乗算器５で行い、この結果をアキュムレータ６
に格納すると共に、畳み込みの演算結果をＲＡＭ４に格
納してディジタルフィルタを得るものである。

【００３１】そして、上記ＲＯＭ３の量子化係数を、
（１）略理想的な量子化係数を求め、所定の係数語長に
するために四捨五入して量子化係数を得る第１ステップ
Ｓ１０１と、（２）理想周波数特性と前記第１ステップ
で得られた量子化係数によるフィルタの周波数特性との
第１の誤差を所定周波数にわたって求める第２ステップ
Ｓ１０２と、（３）前記第１ステップで得られた量子化
係数によるフィルタの周波数特性のうち予め定めたリッ
プル幅を越えた部分である第２の誤差を所定の周波数に
わたって求める第３ステップＳ１０３と、（４）前記第
１の誤差と第２の誤差とを合計して量子化係数の最適値
を評価するための評価値を求める第４ステップＳ１０４
と、（５）前記第４ステップで求められた評価値と、既
に求められている評価値とを比較し、その結果、前記第
４ステップで求められた評価値が、既に求められている
評価値よりよい場合、その量子化係数を保存する第５ス
テップＳ１０５とを、量子化係数の所定のビットを変化
させながら（Ｓ１０７）、前記第１ステップから第５ス
テップまでを繰り返し実行し、最適な量子化係数を得
る。

【００３２】なお上記ステップを記録したコンピュータ
上で読み出し可能な記録媒体にすることで、短時間で、
最適な量子化係数を得ることが出来るようになる。な
お、上記具体例では非巡回形ディジタルフィルタを例に
説明したが巡回形ディジタルフィルタ等あらゆるディジ
タルフィルタに適用出来る。

【００３３】

【発明の効果】本発明に係るディジタルフィルタの量子
化係数の決定方法は上述のように構成したので、フィル
タ回路を追加することなく、リップル幅の小さいディジ
タルフィルタが得られる。しかも、係数語長は所定のビ
ット数で短いから、回路規模は大きくならない。

【図面の簡単な説明】

【図１】本発明の処理ステップを示す流れ図である。

【図２】量子化される前の係数による特性の一例を示す
グラフである。

【図３】誤差関数を説明するグラフである。

【図４】ペナルティ関数を説明するグラフである。

【図５】ハードウエアの構成を示す図である。

【図６】従来技術を示すブロック図である。

【図７】従来技術を示すブロック図である。

【符号の説明】

２シーケンサ３ＲＯＭ４ＲＡＭ５乗算器６アキュムレータＳ１０１第１の工程Ｓ１０２第２の工程Ｓ１０３第３の工程Ｓ１０４第４の工程Ｓ１０５第５の工程

Claims

【特許請求の範囲】

【請求項１】ディジタルフィルタの量子化係数を所定
の係数語長にし、且つ、所定の周波数特性のディジタル
フィルタを得るための、ディジタルフィルタの量子化係
数の決定方法において、最適な量子化係数を得るため
に、得られたディジタルフィルタの周波数特性の、予め
定めたリップル幅を越えた量を評価すると共に前記量子
化係数を変化させながら前記評価を繰り返し実行するこ
とで、最適な量子化係数を決定することを特徴とするデ
ィジタルフィルタの量子化係数の決定方法。
【請求項２】前記リップル幅を越えた量と、理想周波
数特性と前記得られたディジタルフィルタの周波数特性
との差とを合計した評価値に基づきこの評価値を最小に
することで、最適な量子化係数を決定することを特徴と
する請求項１記載のディジタルフィルタの量子化係数の
決定方法。
【請求項３】前記評価値を得る際、前記リップル幅を
越えた量には重みづけを行い演算することを特徴とする
請求項１又は２記載のディジタルフィルタの量子化係数
の決定方法。
【請求項４】ディジタルフィルタの量子化係数の決定
方法において、（１）略理想的な量子化係数を求め、所
定の係数語長にするために四捨五入して量子化係数を得
る第１ステップと、（２）理想周波数特性と前記第１ス
テップで得られた量子化係数によるフィルタの周波数特
性との第１の誤差を所定周波数にわたって求める第２ス
テップと、（３）前記第１ステップで得られた量子化係
数によるフィルタの周波数特性のうち予め定めたリップ
ル幅を越えた部分である第２の誤差を所定の周波数にわ
たって求める第３ステップと、（４）前記第１の誤差と
第２の誤差とを合計して量子化係数の最適値を評価する
ための評価値を求める第４ステップと、（５）前記第４
ステップで求められた評価値と、既に求められている評
価値とを比較し、その結果、前記第４ステップで求めら
れた評価値が、既に求められている評価値よりよい場合
その量子化係数を保存する第５ステップとからなり、量子化係数の所定のビットを変化させながら前記第１ス
テップから第５ステップまでを繰り返し実行し、最適な
量子化係数を決定することを特徴とするディジタルフィ
ルタの量子化係数の決定方法。
【請求項５】前記第４ステップにおいて、評価値を求
める際、第２の誤差に重みづけをした後第１の誤差と合
計して評価値を求めることを特徴とする請求項４記載の
ディジタルフィルタの量子化係数の決定方法。
【請求項６】シーケンサにコントロールされてＲＯＭ
に格納された量子化係数と、ＲＡＭに格納された畳み込
みの演算結果とを取り出して所定の乗算を行い、この結
果をアキュムレータに格納すると共に、畳み込みの結果
を前記ＲＡＭに格納するハードウェアを含み、ディジタ
ルフィルタの量子化係数を所定の係数語長にし、且つ、
所定の周波数特性のディジタルフィルタを得るための、
ディジタルフィルタの量子化係数の決定方法において、
最適な量子化係数を得るために、得られたディジタルフ
ィルタの周波数特性の、予め定めたリップル幅を越えた
量を評価すると共に前記量子化係数を変化させながら前
記評価を繰り返し実行することで、最適な量子化係数を
決定することを特徴とするディジタルフィルタの量子化
係数の決定プログラムが記録されたコンピュータ上で読
み出し可能な記録媒体。
【請求項７】前記リップル幅を越えた量と、理想周波
数特性と前記得られたディジタルフィルタの周波数特性
との差とを合計した評価値に基づき評価値を最小にする
ことで、最適な量子化係数を決定することを特徴とする
請求項６記載のディジタルフィルタの量子化係数の決定
プログラムが記録されたコンピュータ上で読み出し可能
な記録媒体。
【請求項８】前記評価値を得る際、前記リップル幅を
越えた量には重みづけを行い演算することを特徴とする
請求項６又は７記載のディジタルフィルタの量子化係数
の決定プログラムが記載されたコンピュータ上で読み出
し可能な記録媒体。
【請求項９】シーケンサにコントロールされてＲＯＭ
に格納された量子化係数と、ＲＡＭに格納された畳み込
みの演算結果とを取り出して所定の乗算を乗算器で行
い、この結果をアキュムレータに格納すると共に、畳み
込みの演算結果を前記ＲＡＭに格納するディジタルフィ
ルタの前記量子化係数の決定方法において、（１）略理
想的な量子化係数を求め、所定の係数語長にするために
四捨五入して量子化係数を得る第１ステップと、（２）
理想周波数特性と前記第１ステップで得られた量子化係
数によるフィルタの周波数特性との第１の誤差を所定周
波数にわたって求める第２ステップと、（３）前記第１
ステップで得られた量子化係数によるフィルタの周波数
特性のうち予め定めたリップル幅を越えた部分である第
２の誤差を所定の周波数にわたって求める第３ステップ
と、（４）前記第１の誤差と第２の誤差とを合計して量
子化係数の最適値を評価するための評価値を求める第４
ステップと、（５）前記第４ステップで求められた評価
値と、既に求められている評価値とを比較し、その結
果、前記第４ステップで求められた評価値が、既に求め
られている評価値よりよい場合その量子化係数を保存す
る第５ステップとからなり、量子化係数の所定のビットを変化させながら前記第１ス
テップから第５ステップまでを繰り返し実行し、最適な
量子化係数を決定することを特徴とするディジタルフィ
ルタの量子化係数の決定プログラムが記録されたコンピ
ュータ上で読み出し可能な記録媒体。
【請求項１０】前記第４ステップにおいて、評価値を
求める際、第２の誤差に重みづけをした後、第１の誤差
と合計して評価値を求めることを特徴とする請求項９記
載のディジタルフィルタの量子化係数の決定プログラム
が記録されたコンピュータ上で読み出し可能な記録媒
体。