(57)【要約】【課題】本発明は阻止域が可変なフィルタ装置及びディジタルフィルタの設計方法に関し、簡単な構成で、瞬時に阻止域を変更し、効率よくノイズを低減し、高速、高精度にフィルタリングを行なうことができるフィルタ装置及びディジタルフィルタの設計方法を提供することを目的とする。【解決手段】本発明は、係数に応じた複数の阻止帯域を有するディジタルフィルタの複数の阻止域のうち少なくとも一つの阻止域が所望の周波数と一致するようにディジタルフィルタの係数を最小自乗近似問題又は半正定値計画問題に定式化して求めることを特徴とする。【選択図】図２
2007142563.tif 000002
【特許請求の範囲】
【請求項１】
係数に応じた複数の阻止帯域を有するディジタルフィルタと、
前記ディジタルフィルタの複数の阻止域のうち少なくとも一つの阻止域が所望の周波数と一致するように前記ディジタルフィルタの係数を求める係数設定手段を有することを特徴とするフィルタ装置。
【請求項２】
前記係数設定手段は、最小自乗近似問題に定式化して係数を求めることを特徴とする請求項１記載のフィルタ装置。
【請求項３】
前記係数設定手段は、ωを正規化角周波数、φ1,φ2,・・・,φKを大きな阻止域減衰量を必要とする各帯域の中心正規化角周波数、τ1,τ2,・・・,τKを大きな阻止域減衰量を必要とする各帯域幅、δ1,δ2,・・・,δKを大きな阻止域減衰量を必要とする各帯域の中心正規化角周波数φ1,φ2,・・・,φKのスペクトルパラメータとし、
所望の周波数特性を
【数１】
2007142563.tif 000032
伝達関数を
【数２】
2007142563.tif 000033
重み関数を
【数３】
2007142563.tif 000034
としたときの評価関数が
【数４】
2007142563.tif 000035
で表し、
評価関数が最小となるように係数
【数５】
2007142563.tif 000036
を求めることを特徴とする請求項２記載のフィルタ装置。
【請求項４】
前記係数設定手段は、半正定値計画問題に定式化して係数を求めることを特徴とする請求項１記載のフィルタ装置。
【請求項５】
係数に応じた複数の阻止帯域を有するディジタルフィルタの設計方法であって、
前記ディジタルフィルタの複数の阻止域のうち少なくとも一つの阻止域が所望の周波数と一致するように前記ディジタルフィルタの係数を求めることを特徴とするディジタルフィルタの設計方法。
【請求項６】
最小自乗近似問題に定式化して前記係数を求めることを特徴とする請求項５記載のディジタルフィルタの設計方法。
【請求項７】
ωを正規化角周波数、φ1,φ2,・・・,φKを大きな阻止域減衰量を必要とする各帯域の中心正規化角周波数、τ1,τ2,・・・,τKを大きな阻止域減衰量を必要とする各帯域幅、δ1,δ2,・・・,δKを大きな阻止域減衰量を必要とする各帯域の中心正規化角周波数φ1,φ2,・・・,φKのスペクトルパラメータとし、
所望の周波数特性を
【数６】
2007142563.tif 000037
伝達関数を
【数７】
2007142563.tif 000038
重み関数を
【数８】
2007142563.tif 000039
としたときの評価関数が
【数９】
2007142563.tif 000040
で表し、
評価関数が最小となるように係数
【数１０】
2007142563.tif 000041
を求めることを特徴とする請求項６記載のディジタルフィルタの設計方法。
【請求項８】
前記係数設定手段は、半正定値計画問題に定式化して係数を求めることを特徴とする請求項５記載のディジタルフィルタの設計方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はフィルタ装置及びディジタルフィルタの設計方法に係り、特に、阻止域が可変なフィルタ装置及びディジタルフィルタの設計方法に関する。
【背景技術】
【０００２】
ディジタル信号処理の分野において、処理目的に応じて瞬時に周波数特性を変化させたいことがある。これを実現するために可変フィルタが提案されている。可変フィルタは、信号処理を行う過程でディジタルフィルタの周波数特性を瞬時に変化可能なフィルタである。近年、可変ディジタルフィルタは、通信、計測、音響、画像処理などの分野で注目されている。可変ディジタルフィルタとしては、安定性に優れたＦＩＲディジタルフィルタが使われている。
【０００３】
しかしながら、一般にＦＩＲディジタルフィルタは、大きな阻止域減衰量を得るためには高い次数が必要となり、遅延時間が増大するため、高速高精度計測を行う際にしばしば問題となる。そこで、大きなノイズを発生する主要な部分に対して、部分的に大きな阻止域減衰量を持つ直線位相ＦＩＲディジタルフィルタを実現することで、阻止域全体が大きな阻止域減衰量をもつフィルタより、フィルタ次数の低減を行い遅延時間の短縮を行う方法が考えられる。この方法において、大きなノイズが発生する部分が変化した場合、それに併せてフィルタの特性を変化させることが重要となる。
【０００４】
そこで、本出願人は、スペクトルパラメータを用いて阻止域の一部が可変可能なＦＩＲディジタルフィルタの設計法を提案した（特許文献１、非特許文献１）、従来と同精度の計測を従来よりも高速に行うことが出来ることを示した。  
【特許文献１】特開２００４−１５０８８３号公報
【非特許文献１】新保友成、杉田泰則、相川直幸、森地徹、若狭由喜夫、木村隆、”阻止域の一部が可変なＦＩＲフィルタの設計法”電子情報通信学会論文誌（Ａ），Vol.J87-A，No.12，pp.1511-1517，(2004，12)
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、大きなノイズが発生する原因が複数あり、かつ、それらが変化する場合、従来のフィルタを適用することが出来なかった。すなわち、阻止域全体、あるいは、大きなノイズが発生する可能性のある帯域全体の減衰量を大きくする必要があり、フィルタの次数の低減効果はあまり望めない。その結果、効率よくノイズを低減し、高速、高精度にフィルタリングを行なうことができないなどの課題があった。
【０００６】
本発明は上記の点に鑑みてなされたもので、簡単な構成で、瞬時に阻止域を変更し、効率よくノイズを低減し、高速、高精度にフィルタリングを行なうことができるフィルタ装置及びディジタルフィルタの設計方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、係数に応じた複数の阻止帯域を有するディジタルフィルタと、ディジタルフィルタの複数の阻止域のうち少なくとも一つの阻止域が所望の周波数と一致するようにディジタルフィルタの係数を求める係数設定手段を有することを特徴とする。
【０００８】
また、係数設定手段は、最小自乗近似問題又は半正定値計画問題に定式化して係数を求めることを特徴とする。
【発明の効果】
【０００９】
本発明によれば、係数に応じた複数の阻止帯域を有するディジタルフィルタの複数の阻止域のうち少なくとも一つの阻止域が所望の周波数と一致するようにディジタルフィルタの係数を求めることにより、簡単な構成で、瞬時に阻止域を変更し、効率よくノイズを低減し、高速、高精度にフィルタリングを行なうことができる。また、このとき、最小自乗近似問題又は半正定値計画問題に定式化して係数を求めることにより、フィルタの次数を挙げることなく、効率よく係数の演算を行なえる。
【発明を実施するための最良の形態】
【００１０】
〔第１実施例〕
図１は本発明の第１実施例のブロック構成図を示す。
【００１１】
本実施例のフィルタ装置１００は、アンプ１１１、アナログフィルタ１１２、アナログ−ディジタル変換器１１３、ディジタル信号処理部１１４、係数演算部１１５から構成されている。アンプ１１１は、入力端子Ｔinに供給された入力信号を増幅する。アンプ１１１で増幅された入力信号は、アナログ−ディジタル変換器１１３に供給される。アナログ−ディジタル変換器１１３は、アンプ１１１からのアナログ信号をディジタル信号に変換する。
【００１２】
アナログ−ディジタル変換器１１３で変換されたディジタル信号は、ディジタル信号処理部１１４に供給される。ディジタル信号処理部１１４は、ＦＩＲフィルタであり、複数の阻止域を有する周波数特性を有する。ディジタル信号処理部１１４の周波数特性は、係数演算部１１５から供給される係数によって決定されている。ディジタル信号処理部１１４でフィルタリングされたディジタルデータは、出力端子Ｔoutから出力される。また、係数演算部１１５は、コントロール端子Ｔcntから供給されるコントロール信号に応じて係数を決定する。コントロール端子Ｔcntに供給されるコントロール信号は、ノイズなどの周波数に応じて決定される。係数演算部１１５は、ディジタル信号処理部１１４の複数の阻止域のうちいずれかの阻止域がノイズなどの周波数と一致するように係数を演算し、ディジタル信号処理部１１４に供給する。
【００１３】
図２は本発明の第１実施例の処理フローチャートを示す。
【００１４】
本実施例のフィルタ装置１００は、ステップＳ１−１で係数演算部１１５によって演算された係数をディジタル信号処理部１１４に設定する。フィルタ装置１００は、ステップＳ１−２で設定された係数に応じた周波数特性により信号処理を実行する。フィルタ装置１００はステップＳ１−３でコントロール端子Ｔcntからコントロール信号を読み取って、ステップＳ１−４で最小自乗法により係数を演算する。  
なお、本実施例では、ディジタル信号処理部１１４に設定されるＦＩＲディジタルフィルタは、複数の大きな阻止域減衰量をもつ帯域が、スペクトルパラメータを用いることにより変化可能に設計されている。また、このとき、本実施例では、伝達関数がスペクトルパラメータ数の次元を有する多項式で表されており、伝達関数を周波数領域設計問題における重み付き最小自乗近似問題に定式化し、多項式係数を求める。さらに、得られた多項式係数より係数を求める。
【００１５】
フィルタ装置１００は、ステップＳ１−５でフィルタリング処理が終了するまで、ステップＳ１−１〜Ｓ１−４の動作を繰り返す。
〔複数の阻止域可変フィルタ〕
次にステップＳ１−４における係数の演算方法について説明する。
【００１６】
図３は複数の阻止域が可変な理想振幅特性を有する周波数特性図を示す。
【００１７】
図３に示す複数の阻止帯域が可変な理想振幅周波数特性の関数Ｄは、
【００１８】
【数１】
2007142563.tif 000003
で表せる。式（１）において、ωpは通過域端正規化角周波数、ωsは阻止域端正規化角周波数、φ1,φ2,・・・,φKは大きな阻止域減衰量を必要とする各帯域の中心正規化角周波数、τ1,τ2,・・・,τKは大きな阻止域減衰量を必要とする各帯域幅、δ1,δ2,・・・,δKは大きな阻止域減衰量を必要とする各帯域の中心正規化角周波数φ1,φ2,・・・,φKのスペクトルパラメータを表す。このとき、スペクトルパラメータδ1,δ2,・・・,δKは、
【００１９】
【数２】
2007142563.tif 000004
とする。
【００２０】
また、重み関数Ｗを
【００２１】
【数３】
2007142563.tif 000005
とする。ここで、ｇ1,ｇ2,・・・,ｇKは、大きな阻止域減衰量を必要とする各帯域の重みで正の実数値である。
【００２２】
また、上述の所望特性と重みを有する特性を近似する２Ｎ次の複数の阻止域可変直線位相ＦＩＲディジタルフィルタの伝達関数の周波数応答をスペクトルパラメータδ1,δ2,・・・,δKを用いて
【００２３】
【数４】
2007142563.tif 000006
と表す。ここで、
【００２４】
【数５】
2007142563.tif 000007
はスペクトルパラメータが変化するとその値が変化するので、参考文献「新保  友成、杉田  泰則、相川  直幸、木村  隆、若狭  由喜夫、森地  徹”阻止域の一部が可変なＦＩＲフィルタの設計法”、信学論、vol.j85-a,no.12,dec.2002」に示される記述を拡張して以下のように各スペクトルパラメータに対して、それぞれＬｉ 次の多項式を用いて係数を近似することを考えと、
【００２５】
【数６】
2007142563.tif 000008
で表すことができる。
【００２６】
なお、以下、議論を簡単にするためにＫ＝２の場合について説明を行なう。
【００２７】
式（５）を式（４）に代入した伝達関数の振幅特性
【００２８】
【数７】
2007142563.tif 000009
が、式（１）に示す理想振幅特性の関数と式（３）に示す重み関数との間で、最小自乗の意味で最適となる評価関数を
【００２９】
【数８】
2007142563.tif 000010
と与える。ここで、重みの正の実数値ｇ(i,l1,l2)は、評価関数Ｊが最小となるように求めることである。評価関数Ｊにおいて各スペクトルパラメータδ1,δ2,・・・,δKの変化量を
【００３０】
【数９】
2007142563.tif 000011
個に等分し、０からπまでの周波数領域を
【００３１】
【数１０】
2007142563.tif 000012
個に等分すると、離散した評価関数
【００３２】
【数１１】
2007142563.tif 000013
が得られる。式（７）で表される離散した関数を最小化するためには、係数ｇ(p,q,r)に関数の偏微分を求め、零とすればよい。したがって、
【００３３】
【数１２】
2007142563.tif 000014
で表せる。ただし、
【００３４】
【数１３】
2007142563.tif 000015
である。
【００３５】
式（８）−（１１）を整理すると、
【００３６】
【数１４】
2007142563.tif 000016
となる。そこで、式（１２）においてｐ＝０,１,・・・,Ｎ、ｑ＝０,１,・・・,Ｌ1、ｒ＝０,１,・・・,Ｌ2を代入し、スペクトルパラメータに関する多次元配列を１次元配列に変換すると、以下の行列が得られる。
【００３７】
【数１５】
2007142563.tif 000017
行列Ａのi,j番目の要素a(i,j)は、
【００３８】
【数１６】
2007142563.tif 000018
であり、ｊは以下のように(N+1)(L1+1)(L2+1)まで連続的に変化し、それに対応する
【００３９】
【数１７】
2007142563.tif 000019
は以下のように変化する。
【００４０】
【数１８】
2007142563.tif 000020
以上により、最適な係数は、
【００４１】
【数１９】
2007142563.tif 000021
となる。
〔設計例〕
次に設計例として、大きな減衰量をもつ阻止域が２つ可変可能な仕様を考える．
［設計例の仕様］
フィルタ次数：60
各多項式次数：L1=L2=6
通過域：0-0.1π
阻止域：0.2π-π
重み：１
また、大きな減衰量をもつ阻止域について
中心角周波数：φ1=0.35π,φ2=0.75π
帯域幅：τ1=τ2=0.1
角周波数変化量：δ1=δ2=-0.05〜0.05
各帯域の重み：10
とし、通過域を20等分、阻止域全体を500等分に離散化し、また、大きな減衰量の角周波数変化量についてもそれぞれ0.05刻みに離散化した仕様とした。
【００４２】
図４はスペクトルパラメータδ1,δ2に対する小さな減衰量の最大誤差、大きな減衰量の最大誤差を示す図、図５は図４に示す各スペクトルパラメータδ1,δ2における周波数特性図を示す。
【００４３】
図５に示すようにスペクトルパラメータδ1,δ2が変化したときに、振幅特性がそれぞれ変化させることができることが分かる。よって、スペクトルパラメータδ1,δ2を変化させることによって、所望の振幅特性に近似させることができる。
【００４４】
なお、このとき、図５（Ｊ）に示すようにスペクトルパラメータδ1=0，δ2=0.01の離散化していない点についても小さな減衰量のところで誤差が0.0034(-49dB)、大きな減衰量のところで誤差が0.000013(-97dB)とよい振幅特性が得られていることが分かる。この結果、連続的に複数の阻止域を変化することが可能であることが分かる。
【００４５】
〔第２実施例〕
なお、第１実施例では、最小自乗法により複数の阻止域を可変したが、半正定値計画（semi definite programming,ＳＤＰ）法によって複数の阻止域を可変するようにしてもよい。
【００４６】
次にフィルタ係数の半正定値計画法による近似解法について説明する。
【００４７】
Ｈを周波数応答、ａをフィルタ係数、ｇを第１のスペクトルパラメータ、ｑを第２のスペクトルパラメータ、Ｎをフィルタ次数、ωを角周波数とすると、
周波数応答Ｈ(ω,ｇ，ｑ)は、
【００４８】
【数２０】
2007142563.tif 000022
で表される。
【００４９】
ここで、行列は、
【００５０】
【数２１】
2007142563.tif 000023
で表される。なお、ここで、求めたいフィルタ係数は、行列
【００５１】
【数２２】
2007142563.tif 000024
で表されている。
【００５２】
ここで、Wを重み関数、Dを理想特性の関数、eを誤差関数とすると、
誤差関数ｅは、
【００５３】
【数２３】
2007142563.tif 000025
で表せる。よって、誤差関数ｅの最大点を最小にする条件、
【００５４】
【数２４】
2007142563.tif 000026
を満足するように、係数を決定する。
【００５５】
ここで、最大許容誤差をλとすると、各評価点において
【００５６】
【数２５】
2007142563.tif 000027
を満足すればよい。ここで、
【００５７】
【数２６】
2007142563.tif 000028
とする。
【００５８】
すべての評価点について考えると、評価関数Ｆは、
【００５９】
【数２７】
2007142563.tif 000029
で表される。また、このときの条件は、
【００６０】
【数２８】
2007142563.tif 000030
で表せる。
【００６１】
なお、このとき、
【００６２】
【数２９】
2007142563.tif 000031
である。
【００６３】
〔設計例〕
図６は本発明の第２実施例の図４に示す各スペクトルパラメータδ1,δ2における周波数特性図を示す。
【００６４】
図６に示すように半正定値計画問題を解くことにより複数の阻止域を有するＦＩＲディジタルフィルタを設計した場合においても、スペクトルパラメータδ1,δ2が変化したときに、振幅特性がそれぞれ変化させることができることが分かる。よって、スペクトルパラメータδ1,δ2を変化させることによって、所望の振幅特性に近似させることができる。
【００６５】
上記の条件から半正定値計画によってパラメータを決定することが可能となる。
【００６６】
半正定値計画を用いることにより、簡単な処理で所望の周波数応答を得ることが可能となり、短時間で複数の阻止域を有するＦＩＲフィルタの振幅特性を可変することができる。
【図面の簡単な説明】
【００６７】
【図１】本発明の第１実施例のブロック構成図である。
【図２】本発明の第１実施例の処理フローチャートである。
【図３】複数の阻止域が可変な理想振幅特性を有する周波数特性図である。
【図４】スペクトルパラメータδ1,δ2に対する小さな減衰量の最大誤差、大きな減衰量の最大誤差を示すである。
【図５】図４に示す各スペクトルパラメータδ1,δ2における周波数特性図である。
【図６】本発明の第２実施例の図４に示す各スペクトルパラメータδ1,δ2における周波数特性図である。
【符号の説明】
【００６８】
１００  フィルタ装置
１１１  アンプ、１１２  アナログフィルタ、１１３  アナログ−ディジタル変換器
１１４  ディジタル信号処理部、１１５  係数演算部
【図１】
2007142563.tif 000042
【図２】
2007142563.tif 000043
【図３】
2007142563.tif 000044
【図４】
2007142563.tif 000045
【図５】
2007142563.tif 000046
【図６】
2007142563.tif 000047