"電力系統の交流信号の周波数を検出する周波数検出装置であって、前記交流信号に含まれる不平衡成分と高調波成分を除去する複素係数フィルタからなるフィルタ手段と、前記フィルタ手段から出力される交流信号を用いて前記電力系統の交流信号の基本波成分の周波数を検出する周波数検出手段と、を備えたことを特徴とする周波数検出装置。前記フィルタ手段は、中心周波数が前記電力系統の交流信号の周波数に設定された帯域通過型の複素係数フィルタである、請求項1に記載の周波数検出装置。前記周波数検出手段は、中心周波数が前記電力系統の交流信号の周波数に設定され、通過帯域で位相差が直線的に変化する位相特性を有する通過帯域型の複数係数フィルタからなる第2のフィルタ手段と、前記第2のフィルタ手段に入力される入力信号と前記第2のフィルタ手段から出力される出力信号との位相差を算出する位相差算出手段と、前記第2のフィルタ手段の前記位相特性を用いて前記位相差算出手段で算出された位相差から前記交流信号の基本波成分の周波数を算出する周波数算出手段と、を含む、請求項1又は2に記載の周波数検出装置。前記周波数算出手段は、前記位相特性に基づき前記通過帯域における位相差から周波数を求めるための所定の関係式が設定されており、前記位相差算出手段で算出された算出値に対して前記関係式の演算処理をすることにより前記交流信号の基本波成分の周波数を算出する、請求項3に記載の周波数検出装置。前記関係式は、周波数fが位相差ψに対してf=pψqp:係数、q:定数で表わされる関係式である、請求項4に記載の周波数検出装置。前記周波数算出手段は、前記位相特性に基づき前記通過帯域における位相差と周波数との関係を示す所定のテーブルが設定されており、そのテーブルから前記位相差算出手段で算出された算出値に対応する周波数を読み出すことにより前記交流信号の基本波成分の周波数を求める、請求項1又は2に記載の周波数検出装置。前記周波数検出手段は、位相差が中心周波数では零で、当該中心周波数より大きい周波数領域では負になり、小さい周波数領域では正になる位相特性を有し、かつ、前記中心周波数が変更可能な通過帯域型の複数係数フィルタからなる第2のフィルタ手段と、前記第2のフィルタ手段に入力される入力信号と前記第2のフィルタ手段から出力される出力信号との位相差を算出する位相差算出手段と、前記位相差算出手段で算出される位相差が正の場合は、その位相差が零になるまで前記第2のフィルタ手段の中心周波数を所定の変化量で低下させ、前記位相差算出手段で算出される位相差が負の場合は、その位相差が零になるまで前記第2のフィルタ手段の中心周波数を前記変化量で上昇させ、変化後の中心周波数を周波数の検出値として出力する中心周波数制御手段と、を含む、請求項1又は2に記載の周波数検出装置。前記フィルタ手段は、阻止周波数が前記電力系統の交流信号の周波数の不平衡成分と所定の次数の高調波成分とに設定され、当該阻止周波数が変更可能な通過阻止型の複数係数フィルタであり、前記中心周波数制御手段から出力される中心周波数に基づいて、当該中心周波数に対して不平衡成分と所定次数の高調波成分となるように前記フィルタ手段の阻止周波数を制御する阻止周波数制御手段を更に備える、請求項7に記載の周波数検出装置。前記フィルタ手段の前記阻止周波数制御手段には、更に中心周波数が前記電力系統の系統周波数に設定された帯域通過型の複素係数フィルタが組み合わされている、請求項8に記載の周波数検出装置。前記周波数検出手段は、前記フィルタ手段から出力される交流信号のレベルがゼロレベルを交差するゼロクロスタイミングを検出する検出手段と、前記検出手段により検出される前記ゼロクロスタイミングを用いて前記ゼロクロスタイミング間における所定のクロックパルスのパルス数をカウントし、そのカウント値を用いて前記電力系統の交流信号の基本波成分の周波数を算出する周波数算出手段と、を含む、請求項1又は2に記載の周波数検出装置。電力系統の交流信号には系統周波数fsの基本波正相成分以下、単に基本波成分という。の他に基本波逆相成分以下、不平衡成分という。や低次の高調波成分例えば、5次,7次,11次の高調波成分。以下、単に高調波成分という。が含まれることが多い。従来のゼロクロス点間カウント法による周波数検出方法では、これらの成分が検出電圧に含まれていると、検出電圧の波形が正確に基本波成分の波形にならず波形歪が生じ、基本波成分の波形のゼロクロス点を正確に検出できないので、十分な系統周波数fsの検出精度を得ることができないという問題がある。乗算式PLL法を用いた位相検出装置内の角周波数を利用して周波数を検出する方法でも同様の問題がある。特に、乗算式PLL法を用いた位相検出装置は、装置内で位相を生成し、その位相と入力される交流信号の位相との位相差を算出し、その位相差がゼロとなるように生成位相を制御することによって入力される交流信号の位相を検出するから、位相差が生じた場合は、位相検出装置内で生成される位相が検出対象の位相からずれることになる。従って、位相検出装置内で生成される位相から周波数を算出する方法では、電力系統の交流信号の周波数が変動した場合の検出精度や応答速度の点で問題が生じる。本発明は、上記した事情のもとで考え出されたものであって、不平衡成分や高調波成分の影響を受けることなく高い精度で交流信号の周波数を検出することができる周波数検出装置を提供することを目的とする。請求項1に記載の周波数検出装置は、電力系統の交流信号の周波数を検出する周波数検出装置であって、前記交流信号に含まれる不平衡成分と高調波成分を除去する複素係数フィルタからなるフィルタ手段と、前記フィルタ手段から出力される交流信号を用いて前記電力系統の交流信号の基本波成分の周波数を検出する周波数検出手段と、を備えたことを特徴とする。請求項1に記載の周波数検出装置において、前記フィルタ手段は、中心周波数が前記電力系統の交流信号の周波数に設定された帯域通過型の複素係数フィルタであるとよい請求項2。また、前記周波数検出手段は、中心周波数が前記電力系統の交流信号の周波数に設定され、通過帯域で位相差が直線的に変化する位相特性を有する通過帯域型の複数係数フィルタからなる第2のフィルタ手段と、前記第2のフィルタ手段に入力される入力信号と前記第2のフィルタ手段から出力される出力信号との位相差を算出する位相差算出手段と、前記第2のフィルタ手段の前記位相特性を用いて前記位相差算出手段で算出された位相差から前記交流信号の基本波成分の周波数を算出する周波数算出手段と、を含む構成にするとよい請求項3。また、前記周波数算出手段は、前記位相特性に基づき前記通過帯域における位相差から周波数を求めるための所定の関係式が設定されており、前記位相差算出手段で算出された算出値に対して前記関係式の演算処理をすることにより前記交流信号の基本波成分の周波数を算出するとよい請求項4。例えば、前記関係式は、周波数fが位相差ψに対してf=pψqp:係数、q:定数で表わされる関係式である請求項5。また、前記周波数算出手段は、前記位相特性に基づき前記通過帯域における位相差と周波数との関係を示す所定のテーブルが設定されており、そのテーブルから前記位相差算出手段で算出された算出値に対応する周波数を読み出すことにより前記交流信号の基本波成分の周波数を求めるとよい請求項6。また、請求項1又は2に記載の周波数検出装置において、前記周波数検出手段は、位相差が中心周波数では零で、当該中心周波数より大きい周波数領域では負になり、小さい周波数領域では正になる位相特性を有し、かつ、前記中心周波数が変更可能な通過帯域型の複数係数フィルタからなる第2のフィルタ手段と、前記第2のフィルタ手段に入力される入力信号と前記第2のフィルタ手段から出力される出力信号との位相差を算出する位相差算出手段と、前記位相差算出手段で算出される位相差が正の場合は、その位相差が零になるまで前記第2のフィルタ手段の中心周波数を所定の変化量で低下させ、前記位相差算出手段で算出される位相差が負の場合は、その位相差が零になるまで前記第2のフィルタ手段の中心周波数を前記変化量で上昇させ、変化後の中心周波数を周波数の検出値として出力する中心周波数制御手段と、を含む構成にするとよい請求項7。請求項7に記載の周波数検出装置において、前記フィルタ手段は、阻止周波数が前記電力系統の交流信号の周波数の不平衡成分と所定の次数の高調波成分とに設定され、当該阻止周波数が変更可能な通過阻止型の複数係数フィルタであり、前記中心周波数制御手段から出力される中心周波数に基づいて、当該中心周波数に対して不平衡成分と所定次数の高調波成分となるように前記フィルタ手段の阻止周波数を制御する阻止周波数制御手段を更に備えるとよい請求項8。この場合、前記フィルタ手段の前記阻止周波数制御手段には、更に中心周波数が前記電力系統の交流信号の周波数に設定された帯域通過型の複素係数フィルタが組み合わされているとよい請求項9。また、請求項1又は2に記載の周波数検出装置において、前記周波数検出手段は、前記フィルタ手段から出力される交流信号のレベルがゼロレベルを交差するゼロクロスタイミングを検出する検出手段と、前記検出手段により検出される前記ゼロクロスタイミングを用いて前記ゼロクロスタイミング間における所定のクロックパルスのパルス数をカウントし、そのカウント値を用いて前記電力系統の交流信号の基本波成分の周波数を算出する周波数算出手段と、を含む構成にするとよい請求項10。請求項1又は2に記載の周波数検出装置によれば、周波数検出手段の前段に交流信号に含まれる不平衡成分と高調波成分を除去する複数係数フィルタからなるフィルタ手段を設けているので、入力される電力系統の交流信号に含まれている不平衡成分と高調波成分を複数係数フィルタによって好適に除去することができる。従って、高速かつ高精度で電力系統の周波数を検出することができる。また、請求項3乃至6に記載の周波数検出装置によれば、電力系統の交流信号は、所定のサンプリング周期で、通過帯域で位相差が直線的に変化する位相特性位相差ψと周波数fとがf=pψqの関係を有する特性。p,qは定数を有する複素係数バンドパスフィルタに入力され、その入力信号と出力信号とから位相差ψが算出される。そして、その位相差ψに対してf=pψqの関係式を演算して周波数fが算出され、その算出値が周波数検出値として出力される。複素係数バンドパスフィルタに電力系統の交流信号の基本波成分のみを入力し、その入力信号と出力信号とから位相差ψを算出し、その位相差ψを用いてリニアな位相特性から周波数fを算出するので、電力系統の交流信号の入力に対して高い精度で周波数を連続的に検出することができる。また、請求項7に記載の周波数検出装置によれば、複素係数バンドパスフィルタの入力信号と出力信号とから算出された位相差ψがψ>0であれば、複素係数バンドパスフィルタの中心周波数を低下させ、ψ<0であれば、複素係数バンドパスフィルタの中心周波数を上昇させて、位相差ψがゼロとなるように複素係数バンドパスフィルタの中心周波数が制御される。複素係数バンドパスフィルタの中心周波数がψ=0となる周波数に制御されると、複素係数バンドパスフィルタの中心周波数は入力信号の周波数と一致するから、その制御値が周波数検出値として出力される。複素係数バンドパスフィルタに電力系統の交流信号を入力し、その入力信号と出力信号とから位相差ψを算出し、その位相差ψに基づき複素係数バンドパスフィルタの中心周波数を電力系統の交流信号の基本波成分の周波数に一致させるように制御し、その制御値を周波数検出値として出力するので、電力系統の交流信号の入力に対して高い精度で周波数を連続的に検出することができる。また、請求項8又は9に記載の周波数検出装置によれば、第2のフィルタ手段を複素係数ノッチフィルタとし、位相差ψに基づく複素係数バンドパスフィルタの中心周波数の制御に連動して複素係数ノッチフィルタの阻止周波数を制御するので、不平衡成分や高調波成分を好適に除去することができる。また、請求項10に記載の周波数検出装置によれば、複素係数フィルタにより基本波成分のみが抽出された交流信号を用いてゼロクロス点間カウント法により当該交流信号の周波数を検出するので、従来のゼロクロス点間カウント法による周波数検出装置よりも高速かつ高精度で周波数を検出することができる。以下、本発明の実施の形態を、添付図面を参照して具体的に説明する。図1は、本発明に係る周波数検出装置の第1実施形態のブロック構成を示す図である。図1に示す周波数検出装置1Aは、電力系統の三相交流電圧U,V,Wの各相の相電圧vu,vv,vwを検出した電圧に含まれる不平衡成分や高調波成分など周波数検出処理で外乱となる成分を除去する外乱除去部2Aと、帯域通過型の複素係数フィルタ複素係数バンドパスフィルタの位相特性を用いて外乱除去部2Aから出力される交流信号電力系統の三相交流の基本波成分の周波数を算出する周波数検出部3Aで構成されている。外乱除去部2Aは、電力系統に設けられた電圧検出器から入力される三相交流電圧vu,vv,vw所定のサンプリング周期で入力される瞬時値を互いに直交する二相交流電圧vα,vβに変換する三相/二相変換部11と、三相/二相変換部11から出力される二相交流電圧vα,vβに含まれる不平衡成分周波数-fsの成分と所定次数の高調波成分主として5次高調波成分-5fs、7次高調波成分7fs、11次高調波成分-11fsなどの高調波成分を除去する複素係数フィルタを用いた第1複素係数フィルタ部12と、第1複素係数フィルタ部12から出力される不平衡成分及び高調波成分を含まない二相交流電圧vr,vjを正規化する正規化部13と、を含む。なお、第1複素係数フィルタ部12のゲインを調整することにより正規化部13を省略することができる。三相/二相変換部11は、電圧検出器から入力される三相交流電圧vu,vv,vwを二相交流電圧vα,vβに変換する。対称三相交流電圧vu,vv,vwは、図4aに示すように、u軸を位相θの基準軸とし、u軸に対して±2π/3の角度で開いた方向にv軸とw軸を反時計回りに配置したuvw座標系において、角速度ωで反時計回りに回転する電圧ベクトルVのu軸成分、v軸成分及びw軸成分と考えることができる。また、二相交流電圧vα,vβは、同図aに示すように、uvw座標系にu軸に沿うα軸とそのα軸に直交するβ軸とによるαβ座標系を設定した場合の電圧ベクトルVのα軸成分とβ軸成分と考えることができる。従って、三相二相変換処理は、三相交流電圧vu,vv,vwを二相交流電圧vα,vβに変換する処理である。三相交流電圧vu,vv,vwを二相交流電圧vα,vβに変換する変換式は、である。従って、三相/二相変換部11は、1式,2式の演算を行うことにより、電圧検出器から入力される三相交流電圧vu,vv,vwを二相交流電圧vα,vβに変換する。図2は、三相二相変換部11の演算回路を示すブロック図である。同図に示すように、三相二相変換部11は、5個の乗算器11a~11eと2個の加算器11f,11gで構成される。乗算器11a,11b,11dは、それぞれ1式の各項を演算する演算器であり、乗算器11c,11eは、それぞれ2式の各項を演算する演算器である。また、加算器11fは1式の各項を加算する演算器であり、加算器11gは2式の各項を加算する演算器である。電圧検出器で検出される電力系統の三相交流電圧vu,vv,vwは、一般に、基本波成分以外に不平衡成分や3次、5次、7次、11次などの奇数次の高調波成分図3の周波数成分参照が含まれる非対称三相交流電圧である。従って、三相二相変換部11からはこれらの成分についても三相二相変換した成分を含む二相交流電圧vα’,vβ’が出力される。三相交流電圧vu,vv,vwの基本波成分vsu,vsv,vswは、基本波の電圧ベクトルVsのuvw座標系におけるu、v、wの各軸方向の成分として定義されるものである。一方、三相交流電圧vu,vv,vwの不平衡成分vsu’,vsv’,vsw’は、図4bに示すように、uvw座標系に対してu軸、v軸及びw軸の配列順が逆になっているuwv座標系において、各不平衡成分vsu’,vsv’,vsw’が基本波の電圧ベクトルVsのu、v、wの各軸方向の成分として定義されるものである。従って、基本波成分vsu,vsv,vswを、vsu=Ascos...3Avsv=Ascos...3Bvsw=Ascos...3B但し、As;基本波成分の振幅とすると、不平衡成分vsu’,vsv’,vsw’は、vsu’=As’cos...4Avsv’=As’cos...4Bvsw’=As’cos...4C但し、As’;不平衡成分の振幅で表わされる。3A式~3C式を1式と2式に代入すると、三相/二相変換部11から出力される基本波成分の二相交流電圧vsα,vsβは、vsα=√Ascos...5vsβ=√Assin...6となる。また、4A式~4C式を1式と2式に代入すると、三相/二相変換部11から出力される不平衡成分の二相交流電圧vsα’,vsβ’は、vsα’=√As’cos...7vsβ’=-√As’sin...8となる。cos=cos、-sin=sinであるから、これらを7式、8式に代入すると、vsα’=√As’cos...7’vsβ’=√As’sin...8’となる。7’式及び8’式と5式及び6式を比較すると、基本波成分の角周波数はωsであるのに対し、不平衡成分の角周波数は-ωsである点が相違する。角周波数ωsを正の周波数とすると、不平衡成分の角周波数-ωsは負の周波数となるから、不平衡成分vsu’,vsv’,vsw’を三相二相変換して得られる二相交流電圧vsα’,vsβ’は負の周波数を有する電圧ということができる。二相交流電圧vsα,vsβを成分とする電圧ベクトルVsをVs=Asexpで表わすと、二相交流電圧vsα’,vsβ’を成分とする電圧ベクトルVs’はVs’=As’expで表わされる。図3において、基本波成分を正の周波数領域の周波数fsの位置に表示し、不平衡成分の周波数を-fsとして不平衡成分を負の周波数領域の周波数-fsの位置に表示しているのは上記の周波数の関係を示している。なお、図3には、周波数検出に影響のある5次、7次、11次の高調波成分のみを描いている。3の整数倍の高調波成分は線間電圧には表れず、相電圧でもΔ結線のトランスで除去され、11次よりも大きい奇数次の高調波成分はレベルが小さく、無視し得るから、図3には記載していない。不平衡成分が負の周波数になるのは、不平衡成分の相順uvwが時計回りの順が基本波成分の相順uvwが反時計回りの順に対して逆になるからである。従って、基本波成分の周波数fsをn倍n:2以上の整数したn次高調波成分を三相二相変換した二相交流電圧vnα,vnβ添え字のnは次数。以下、同じの角周波数ωnが基本波成分の二相交流電圧vsα,vsβと同じの符号になる場合は、そのn次高調波成分の周波数fnは正の周波数となり、n次高調波成分vnu,vnv,vnwの相順は基本波成分vsu,vsv,vswの相順と同一になる。逆に、二相交流電圧vnα,vnβの角周波数ωnが基本波成分の二相交流電圧vsα,vsβと逆の符号になる場合は、そのn次高調波成分の周波数fnは負の周波数となり、n次高調波成分vnu,vnv,vnwの相順は基本波成分vsu,vsv,vswの相順と逆になる。5次、7次、11次の高調波成分は、U,V,Wの各相の電圧ベクトルをVnu,Vnv,Vnwと表記し、電圧ベクトルVnをVn=AnexpAn:n次高調波成分の振幅とすると、Vnu=Ancos...9Vnv=Ancos...10Vnw=Ancos...11但し、n=5,7,11で表わされる。9式~11式を1式、2式に代入して5次,7次,11次の高調波成分の二相交流電圧V5α,V5β,V7α,V7β,V11α,V11βを算出すると、v5α=√A5cos...12v5β=√A5sin...13v7α=√A7cos...14v7β=√A7sin...15v11α=√A11cos...16v11β=√A11sin...17となる。12式~17式より、5次,11次の高調波成分は負の周波数を有し、7次の高調波成分は正の周波数を有するから、図3では、5次高調波成分と11次高調波成分は、負の周波数領域の周波数-5fsと-11fsの位置にそれぞれ表示され、7次高調波成分は正の周波数領域の周波数7fsの位置に表示されている。従って、三相/二相変換部11から複素係数フィルタ部12には、5式、6式、7’式、8’式、12式~17式で表わされる基本波成分、不平衡成分及び5次、7次、11次の高調波成分の二相交流電圧vsα,vsβ,vsα’,vsβ’,vnα,vnβn=5,7,11を含む二相交流電圧vα’,vβ’が出力される。三相交流電圧vu,vv,vwに、不平衡成分や5次,7次,11次の高調波成分が含まれている場合、加算器11fからは5式,7’式,12式,14式,16式の演算結果の合成値が出力され、加算器11gからは6式,8’式,13式,15式,17式の演算結果の合成値が出力される。第1複素係数フィルタ部12は、z変換表現による伝達関数Hが下記の18式で表される1次のIIRフィルタからなる複素係数バンドパスフィルタBPFで構成される。18式において、複素係数a1におけるfd[Hz]は、通過帯域の中心周波数f0をサンプリング周波数で正規化した正規化周波数である。また、Ωd[rad/s]は、正規化角周波数である。例えば、サンプリング周波数をfsrとし、中心周波数f0を系統周波数fsに設定すると、fdはfs/fsr、Ωdは2πfd=2πfs/fsrとなる。なお、正規化した角周波数Ωdは、-π<Ωd<πである。また、rは、通過帯域の帯域幅を決めるパラメータ0<r<1である。図5は、上記18式の演算処理を行う処理回路を示すブロック図である。同図に示すように、第1複素係数フィルタ部12は、18式の分母の演算処理がフィードバック回路で構成され、そのフィードバック回路の出力に分子の係数b0を乗算する回路によって構成される。図5に示すブロック図において、ukk:離散時間を表すインデックス番号は入力データ、xkは第1複素係数フィルタ部12の状態データ、ykは第1複素係数フィルタ部12の出力データである。入力データuk、状態データxk及び出力データykの間には、xk=rexpjΩdxk1uk...19yk=1-rxk...20が成立する。第1複素係数フィルタ部12は複素係数バンドパスフィルタで構成されるので、入力データukが複素データか実データ複素データの虚数部が0のデータかに関わらず、状態データxk及び出力データukが複素信号のデータとなる。従って、入力データuk、状態データxk及び出力データykをそれぞれuk=urkjujk、xk=xrkjxjk、yk=yrkjyjkの複素データ、複素係数a1をa1=rexpjΩd=arjaj=rcosjrsinとして19式と20式に代入し、実数部と虚数部の関係式に分けると、xrk=rcosxrk1-rsinxjk1urk...21xjk=rcosxjk1rsinxrk1ujk...22yrk=1-rxrk...23yjk=1-rxjk...24となる。図6は、21式~24式に基づき第1複素係数フィルタ部12の複素演算処理を行う回路構成を示す図である。同図において、係数arと係数ajはそれぞれ複素係数a1=rexpの実数部と虚数部であり、ar=rcos、aj=rsinである。同図に示すように、第1複素係数フィルタ部12は、6個の乗算器12a~12fと、2個の加算器12g,12hと、2個の遅延回路12i,12jで構成される。遅延回路12iは、状態データの実数部xrk1を生成する回路であり、遅延回路12jは、状態データの虚数部xjk1を生成する回路である。乗算器12a,12bはそれぞれ21式の第1項と第2項負の符号を含むを演算する演算器であり、加算器12gは21式の第1項と第2項と第3項を加算する演算器である。従って、加算器12gから21式で示す状態データの実数部xrkが出力される。一方、乗算器12c,12dはそれぞれ22式の第1項と第2項を演算する演算器であり、加算器12hは22式の第1項と第2項と第3項を加算する演算器である。従って、加算器12hから22式で示す状態データの虚数部xjkが出力される。また、乗算器12e,12fはそれぞれ23式と24式を演算する演算器である。本実施形態では、三相二相変換部11を設け、三相交流電圧vu,vv,vwを互いに直交する二相交流電圧vα,vβに変換しているが、二相交流電圧vα,vβは、それぞれ複素データurjujの実数部と虚数部に対応させることができるので、二相交流電圧vαのサンプリングデータを入力データの実数部urkとして加算器12gに入力し、二相交流電圧vβのサンプリングデータを入力データの虚数部ujkとして加算器12hに入力している。二相交流電圧vαのサンプリングデータが第1複素係数フィルタ部12に入力される毎に、遅延回路12i、乗算器12a,12b,12e及び加算器12gで21式及び23式の演算処理が繰り返され、これにより、乗算器12eからは5式で示される基本波成分を三相二相変換した二相交流電圧vsαのみの出力データyrkが出力される。また、二相交流電圧vβのサンプリングデータが第1複素係数フィルタ部12に入力される毎に、遅延回路12j、乗算器12c,12d,12f及び加算器12hで22式及び24式の演算処理が繰り返され、これにより、乗算器12fからは9式で示される基本波成を三相二相変換した二相交流電圧vsβのみの出力データyjkが出力される。バンドパスフィルタを実係数の2次IIRフィルタで構成した場合、その2次IIRフィルタの伝達関数Hz=expjωは、H=rcosz1z1r2z2で表わされる。この伝達関数Hの振幅特性Mを求めると、r2=0を満たすωで極が表れるから、2次IIRフィルタはその極の周波数を通過させる特性を有する。r≒1とすると、cos≒1より、2次IIRフィルタを通過させる正規化周波数fdはfd=±Ωd/2πとなるから、正規化角周波数Ωdを基本波成分の角周波数に設定した実係数の2次IIRフィルタでは、不平衡成分fd=-Ωd/2πの成分も通過させることになる。18式に示す伝達関数Hの振幅特性M求めると、M=1-r/√{1-2rcosr2}となり、r2=0を満たすωだけに極が表れるから、正規化角周波数Ωdを基本波成分の角周波数ω=Ωdに設定した複素係数の1次IIRフィルタでは、基本波成分fd=Ωd/2πの成分だけを通過させ、不平衡成分や高調波成分を通過させることはない。図7は、不平衡成分と5次、7次、11次の高調波成分を5含む三相交流電圧vu,vv,vwの波形を示す図である。図8は、図7に示す三相交流電圧vu,vv,vwを三相二相変換部11及び中心周波数f0を系統周波数fs60[Hz]に設定した第1複素係数フィルタ部12によってフィルタリング処理をした後、二相交流電圧を三相交流電圧に逆変換して得た三相交流電圧vu’,vv’,vw’をシミュレートした波形を示す図である。中心周波数f0を系統周波数fsに設定した第1複素係数フィルタ部12は、図3に示す周波数特性を有するから、図8に示されるように、三相交流電圧vu,vv,vwに含まれる不平衡成分と5次、7次、11次の高調波成分を第1複素係数フィルタ部12によりを好適に除去できることが分かる。なお、図1では、第1複素係数フィルタ部12から出力データyrk,yjkによって出力される二相交流電圧を、二相交流電圧vα,vβと区別するため、それぞれvr,vjと表記している。第1複素係数フィルタ部12から出力される、5式,6式で表わされる二相交流電圧vr,vjは、U相の電圧vuの位相角φを0とした場合であるが、電圧系統9のU相の電圧vuの位相がずれ、位相角φ≠0の場合は、第1複素係数フィルタ部12から出力される二相交流電圧vr,vjは、vr=Ascos、vj=Assinとなる。正規化部13は、第1複素係数フィルタ部12から出力される二相交流電圧vr,vjのレベルを1に正規化する演算処理を行う。第1複素係数フィルタ部12から出力される二相交流電圧vr,vjは振幅が同一の正弦波と余弦波で、√vr2vj2を演算することにより振幅が求められるから、正規化部13では、出力データyrk,yjkに対してそれぞれyrk/√yrk2yjk2とyjk/√yrk2yjk2の演算処理を行って二相交流電圧vr,vjの正規化処理が行われる。従って、正規化部13からは、vr’=cosとvj’=sinで表わされる信号のデータが出力される。周波数検出部3Aは、通過帯域における位相特性がリニアな複素係数バンドパスフィルタからなる第2複素係数フィルタ部14と、この第2複素係数フィルタ部14に入力される電圧ベクトルVin=expを用いて両電圧ベクトルの位相差ψを算出する位相差算出部15と、位相差算出部15で算出される位相差ψを所定の演算式を用いて電圧ベクトルVinの周波数f=の共役な電圧ベクトルVinは、Vin=expであり、この電圧ベクトルVinと出力電圧ベクトルVout=expjとの乗算を行うと、VinVout=exp=expより位相差ψのベクトルが得られる。電圧ベクトルVin=exp=Rin-jXin、電圧ベクトルVout=expj=RoutjXoutとすると、VinVout=RinRoutXinXoutjRinXout-XinRout=cosjsinであるから、sin=RinXout-XinRoutより、位相差算出部15は、sin1RinXout-XinRoutの演算処理をすることにより位相差ψを算出する。なお、cos=RinRoutXinXoutより、cos1RinRoutXinXoutの演算処理をすることにより位相差ψを算出してもよく、tan=RinXout-XinRout/RinRoutXinXoutより、tan1RinXout-XinRoutRinRoutXinXoutの演算処理をすることにより位相差ψを算出してもよい。図10は、位相差算出部15の演算処理を行う回路構成を示す図である。位相差算出部15には第2複素係数フィルタ部14に入力される電圧ベクトルVinの実数部Rinと虚数部Xin、第2複素係数フィルタ部14から出力される電圧ベクトルVoutの実数部Routと虚数部Xoutの2つのデータが入力される。電圧ベクトルVinの実数部Rinは、正規化部13から出力される電圧vr’=cosのデータであり、虚数部Xinは、正規化部13から出力される電圧vj’=sinのデータである。一方、電圧ベクトルVoutの実数部Routは第2複素係数フィルタ部14から出力される電圧yr=cosのデータであり、虚数部Xoutは第2複素係数フィルタ部14から出力される電圧yj=sinωstψのデータである。従って、位相差ψの演算式は、vr’yj-vj’yr=Ψとすると、ψ=sin1とこの第2複素係数フィルタ部14から出力される電圧ベクトルVout=expを用いて両電圧ベクトルの位相差ψを算出する位相差算出部15と、位相差算出部15で算出される位相差ψに基づいて中心周波数f0を変更する中心周波数制御部17とを含む。位相差算出部15は、第1実施形態の位相差算出部15と同一の構成である。従って、ここでは、第2係数フィルタ部14’と中心周波数制御部17の構成及び演算処理について説明する。第2係数フィルタ部14’に設けられる複素係数バンドパスフィルタは、第1実施形態の第2係数フィルタ部14に設けられる複素係数バンドパスフィルタに対して複素係数a1が中心周波数制御部17からフィードバックされる正規化角周波数Ωd’=2πfd’/fsrを用いて変更可能になされている点が異なる。すなわち、第1実施形態では、図6に示す複素係数バンドパスフィルタの演算回路において、乗算器12a,12dに入力される係数ar複素係数a1の実数部と乗算器12b,12cに入力される係数aj複素係数a1の虚数部が固定であったが、第2実施形態では、乗算器12a,12dに入力される係数arと乗算器12b,12cに入力される係数ajが中心周波数制御部17によって変更設定される正規化角周波数Ωd’によって制御される。図15は、第2複素係数フィルタ部14’における複素演算処理の回路を示す図であるが、同図は、図6に示す第2複素係数フィルタ部14の複素演算処理の回路に対して、複素係数a1の実数部の係数arを演算する係数実数部演算回路12lと虚数部の係数ajを演算する係数虚数部演算回路12kが追加されている点が異なる。従って、係数実数部演算回路12l及び係数虚数部演算回路12k以外の部分の演算処理は、図6に示した回路と同じであるので、ここでは詳細な説明は省略する。係数実数部演算回路12lでは、中心周波数制御部17からフィードバックされる正規化角周波数Ωd’を用いて、ar=rcos...26の演算式により複素係数a1の実数部の係数arが算出される。係数実数部演算回路12lで算出された係数arは乗算器12aと乗算器12dに入力される。また、係数虚数部演算回路12kでは、中心周波数制御部17からフィードバックされる正規化角周波数Ωd’を用いて、aj=rsin...27の演算式により複素係数a1の虚数部の係数ajが算出される。係数虚数部演算回路12kで算出された係数ajは乗算器12bと乗算器12cに入力される。中心周波数制御部17は、位相差算出部15から入力される位相差ψを0と比較し、ψ>0、ψ<0及びψ=0のいずれかの比較結果に応じて第2複素係数フィルタ部14’の正規化角周波数Ωd’を変化させる。Ωd’=2πfd’/fsrで、fd’は第2複素係数フィルタ部14’の中心周波数f0に相当するから、正規化角周波数Ωd’を変化させることにより第2複素係数フィルタ部14’の中心周波数f0が変化することになる。従って、中心周波数制御部17は、正規化角周波数Ωd’を制御することにより実質的に第2複素係数フィルタ部14’の中心周波数f0を制御している。ψ>0であれば、図9の位相特性より第2係数フィルタ部14’に入力される電圧ベクトルVinの周波数fは現在の第2係数フィルタ部14’の中心周波数f0より低い周波数領域にあり、ψ<0であれば、電圧ベクトルVinの周波数fは中心周波数f0より高い周波数領域にあり、ψ=0であれば、電圧ベクトルVinの周波数fは中心周波数f0と同一であることを示すから、ψ>0、ψ<0、ψ=0の各比較結果は、それぞれ低い周波数領域、高い周波数領域、中心周波数と同一の判別結果に対応する。中心周波数制御部17は、位相差算出部15から出力される位相差ψを用いて所定の正規化角周波数Ωd0に対する変化量ΔΩdを算出し、その変化量ΔΩdを正規化角周波数Ωd0に加算又は減算して正規化角周波数Ωdの変化値Ωd’を算出する。所定の正規化角周波数Ωd0は、例えば、fd=60[Hz]系統周波数に設定したときの値である。具体的には、中心周波数制御部17は、位相差算出部15から出力される位相差ψ[rad]に対して、ψ×180/πの演算を行って角度の単位の位相差ψ’[deg]に変換した後、-K×ψ’/fsrK:0~1の所定のゲイン。例えば、K=01、fsr:サンプリング周波数。の演算を行って変化量ΔΩdを算出し、その変化量ΔΩdをΩd0=2πf0/fsrに加算する処理を行う。位相差ψには、25式の演算処理でψ<0又は0≦ψに応じて正負の符号が付されるので、変化量ΔΩdのΩd0への加算処理で、ψ<0の場合には変化量ΔΩdをΩd0から減算し、0≦ψの場合には変化量ΔΩdをΩd0に加算する処理が行われることになる。なお、上記の正規化角周波数Ωdを変化させるための演算処理は一例であり、他の演算方法により位相差ψの符号に応じて正規化角周波数Ωdを増加又は減少させてもよい。周波数検出装置1Bでは、中心周波数制御部17で変更設定された正規化角周波数Ωd’が第2複素係数フィルタ部14’にフィードバックされるとともに、正規化角周波数Ωd’に対してΩd’×fsr/2πの演算処理を行って周波数の変化値fd’が算出され、その変化値fd’が周波数検出値f0’として出力される。中心周波数制御部17は、位相差算出部15から入力される位相差ψが0でなければ、位相差ψが0となる方向に第2複素係数フィルタ部14’の中心周波数を変化させ、位相差ψが0になると、その時の中心周波数を保持するように制御するから、第2複素係数フィルタ部14’の中心周波数は、当該第2複素係数フィルタ部14’の入力信号の周波数fに一致するように制御される。従って、中心周波数制御部17から出力される周波数検出値f0’は、周波数検出装置1Bに入力される三相交流電圧vu,vv,vwの周波数fが安定していれば、その周波数fを示し、周波数fが変動すれば、その変動に追従するように変化し、入力される三相交流電圧vu,vv,vwの周波数fを示すものとなる。図16は、本発明に係る周波数検出装置の第3実施形態のブロック構成を示す図である。第3実施形態に係る周波数検出装置1Cは、第2実施形態に係る周波数検出装置1Bにおいて、第1複素係数フィルタ部12を帯域阻止型の複素係数フィルタ複素係数ノッチフィルタに変更したものである。なお、複素係数ノッチフィルタを用いた場合、入力信号の基本波成分の周波数がずれると、それに伴い不平衡成分と高調波成分の周波数もずれるため、複素係数ノッチフィルタの阻止周波数をずらす必要がある。このため、第3実施形態では、中心周波数制御部17から出力される正規化角周波数Ωd’を外乱除去部2B内の第1複素係数フィルタ部12’にもフィードバックし、複素係数ノッチフィルタの阻止周波数を決定する複素係数を変更するようにしている。図3に示すように、電力系統の交流信号には基本波成分fs以外に不平衡成分-fsや5次、7次、11次の高調波成分-5fs,7fs,-11fsが含まれることが分かっているので、図16に示す周波数検出装置1Cの第1複素係数フィルタ部12’は、それらの周波数毎にz変換表現による伝達関数Hzが下記の28式で表される複素係数ノッチフィルタを設け、それらを多段に接続することによって図18に示す周波数特性を有する複素係数ノッチフィルタを構成している。なお、Ωdは、Ωd=2πnfs/fsrであり、n=-1に設定すると、不平衡成分-fsを阻止周波数とする複素係数ノッチフィルタとなる。また、n=-5、n=7、n=-11に設定すると、それぞれ5次、7次、11次の高調波成分を阻止周波数とする複素係数ノッチフィルタとなる。従って、第1複素係数フィルタ部12’は、図17に示すように、-fs、-5fs、7fs、-11fsを阻止周波数とする複素係数ノッチフィルタ121,122,123,124を縦続接続した構成となっている。また、上記28式の演算処理を行う処理回路のブロック図は、図19に示す構成となり、複素係数ノッチフィルタNFを用いた複素係数フィルタ部12の複素演算処理を行う回路は、図20に示す構成となる。図19は、図5に示すブロック図に対して、入力データukから状態データykを減算し、その減算値を出力データekとして出力する回路を追加したものである。また、図20は、図15に示すブロック図に対して、係数実数部演算回路12lと係数虚数部演算回路12kを不平衡成分又は高調波成分の周波数を中心周波数とする係数arを演算する係数実数部演算回路12l’と係数ajを演算する係数虚数部演算回路12k’に変更している。そして、実数部の乗算器12eの後段に加算器12nを追加し、当該加算器12nで入力データの実数部urkから状態データykの実数部yrkを減算して出力データの実数部erkを出力する構成としている。また、虚数部の乗算器12fの後段に加算器12oを追加し、当該加算器12oで入力データの虚数部ujkから状態データykの虚数部yjkを減算して出力データの虚数部ejkを出力する構成としている。係数実数部演算回路12l’では、中心周波数制御部17からフィードバックされる正規化角周波数Ωd’を用いて、ar=rcos...29の演算式により複素係数a1の実数部の係数arが算出される。また、係数虚数部演算回路12k’では、中心周波数制御部17からフィードバックされる正規化角周波数Ωd’を用いて、aj=rsin...30の演算式により複素係数a1の虚数部の係数ajが算出される。なお、nは、阻止周波数の次数を示し、図17に示された複素係数ノッチフィルタ121ではn=-1、複素係数ノッチフィルタ122,123,124ではそれぞれn=-5,7,-11となる。係数実数部演算回路12l’で算出された係数arは乗算器12aと乗算器12dに入力され、係数虚数部演算回路12k’で算出された係数ajは乗算器12bと乗算器12cに入力される。図20に示す回路は、図15に示す回路に対して、上述した係数実数部演算回路12l’及び係数虚数部演算回路12k’の演算内容の他は、加算器12n,12oでの減算処理が追加された点が異なるだけであるから、図20に示す回路の演算処理の詳細説明は省略する。第3実施形態に係る周波数検出装置1Cでは、複素係数ノッチフィルタ121~124によって不平衡成分-fsの成分や高調波成分-5fs,7fs,-11fsの成分のレベルがそれぞれ抑制されるので、これらの成分が周波数検出部3Bに入力されることを好適に阻止することができる。従って、第3実施形態に係る周波数検出装置1Cでは安定かつ高い精度で電力系統の電圧信号の周波数fを検出することができる。図21は、第3実施形態に係る周波数検出装置1Cで検出される周波数fの応答特性をシミュレーションした結果である。また、図22は、シミュレーション開始から0.3秒後に周波数検出装置1Cから出力される周波数fの変動状態を拡大した図である。シミュレーションの条件は、図12,図13に示すシミュレーション結果と同じである。図21に示すように、第3実施形態の場合もシミュレーション開始から0.2秒後に電力系統の周波数fsを瞬時的に60.0[Hz]から60.3[Hz]に上昇させると、周波数検出装置1Cから出力される周波数fはパルス状に急変するが、周波数急変時時刻0.2秒から凡そ0.05秒後にはほぼ60.3[Hz]に収束することが分かる。また、第3実施形態では、周波数急変時から0.1秒経過した時時刻0.3秒には、図22に示すように、周波数検出装置1Cから出力される周波数fのリップルはなく、第1実施形態よりも高速かつ高精度に電力系統の電圧信号の周波数fsを検出することができることが分かる。図23は、本発明に係る周波数検出装置の第4実施形態のブロック構成を示す図である。第4実施形態に係る周波数検出装置1Dは、周波数検出部3Cをゼロクロス点間カウント法により系統周波数fsを検出する構成としたものである。周波数検出部3Cには電圧検出器で検出された三相交流電圧vu,vv,vwから複素係数バンドパスフィルタによって不平衡成分と高調波成分を除去した基本波成分の三相交流電圧vsu,vsv,vswを入力するため、周波数検出装置1Dの外乱除去部2Cは、第1,第2実施形態に係る周波数検出装置1A,1Bの外乱除去部2Aに対し正規化部13に代えて二相交流電圧vr,vjを三相交流電圧vsu,vsv,vswに変換する二相三相変換部18が設けられている。第1複素係数フィルタ部12から出力される二相交流電圧vr,vjを三相交流電圧vu,vv,vwに変換する変換式は、である。従って、二相/三相変換部18は、31式~33式の演算を行うことにより、第1複素係数フィルタ部12から入力される二相交流電圧vr,vjを三相交流電圧vsu,vsv,vswに変換する。図24は、二相三相変換部18の演算回路を示すブロック図である。同図に示すように、二相三相変換部18は、4個の乗算器18a~18dと2個の加算器18e,18fで構成される。乗算器18aは31式を演算する演算器であり、乗算器18bは32式,33式の第1項を演算する演算器であり、乗算器18c,18dはそれぞれ32式,33式の第2項を演算する演算器である。また、加算器18eは32式の各項を加算する演算器であり、加算器18fは33式の各項を加算する演算器である。周波数検出部3Cは、二相/三相変換部18から入力される三相交流電圧vsu,vsv,vswのレベルをそれぞれ基準レベル0vを比較するレベル比較部19と、レベル比較部19の比較結果に基づいて各三相交流電圧vsu,vsv,vswが基準レベル0vを交差するタイミングゼロクロスタイミングを検出するゼロクロス検出部20と、系統周波数fsよりも高周波のクロックパルスを発生するパルス発生部21と、ゼロクロス検出部20から出力されるゼロクロスタイミングによってパルス発生部21から出力されるクロックパルスのカウントを制御することより各三相交流電圧vsu,vsv,vswの半周期若しくは1周期の時間を計測し、その時間を用いて系統周波数fsを算出する周波数算出部22と、を含む。レベル比較部19は、所定のサンプリング周期で入力される三相交流電圧vsuの電圧レベルを基準レベルゼロレベルと比較する。ゼロクロス検出部20は、レベル比較部19の比較結果に基づいて、電圧レベルがゼロレベルを超えると、ハイレベルになり、電圧レベルがゼロレベル以下になると、ローレベルになるパルスPuを周波数算出部22に出力する。レベル比較部19は、三相交流電圧vsv,vswについても同様の処理を行い、それぞれパルスPvとパルスPwを周波数算出部22に出力する。図25に示すように、パルスPuの立上がりタイミングは、三相交流電圧vsuが負レベルから正レベルにゼロレベルをクロスする点tuに対応し、パルスPuの立下がりタイミングは、三相交流電圧vsuが正レベルから負レベルにゼロレベルをクロスする点tuに対応する。同様に、パルスPvの立上がりタイミングと立下がりタイミングは、三相交流電圧vsvが負レベルから正レベルにゼロレベルをクロスする点tvと三相交流電圧vsvが正レベルから負レベルにゼロレベルをクロスする点tvとに対応し、パルスPwの立上がりタイミングと立下がりタイミングは、三相交流電圧vswが負レベルから正レベルにゼロレベルをクロスする点twと三相交流電圧vswが正レベルから負レベルにゼロレベルをクロスする点twとに対応する。周波数算出部22にはパルス発生部21から系統周波数fsの数百倍~千数百倍の周波数を有するクロックパルスCKが入力されている。図8に示したように、二相/三相変換部18から入力される三相交流電圧vsu,vsv,vswは電力系統の三相交流電圧の基本波成分対称三相交流であるので、周波数算出部22で算出される三相交流電圧vsu,vsv,vswの周波数fsu,fsv,fswは同一になる。従って、周波数算出部22にはレベル比較部19から三相交流電圧vsu,vsv,vswのゼロクロスタイミングの検出信号が電力系統の交流信号の周期の1/6の周期で入力される。従って、周波数算出部22は、図25に示すように、ゼロクロス検出部20から入力されるパルスPuの立上がりタイミングtuでクロックパルスCKのカウントを開始し、その立上がりタイミングtuに続くパルスPwの立下がりタイミングtwでクロックパルスCKのカウントを停止し、そのカウント値Nを一時保存する。周波数算出部22は、カウント値NにクロックパルスCKの周期τ[sec]を乗じてゼロクロス点tu-tw間の時間T=τN[sec]を算出し、そのゼロクロス点間の時間Tの6倍の時間の逆数を演算して電力系統の交流信号の周波数fs=1/6T=1/6τNを算出し、その算出結果を周波数検出値として出力する。また、周波数算出部22は、ゼロクロス検出部20から入力されるパルスPwの立下がりタイミングtwでクロックパルスCKのカウントを開始し、その立下がりタイミングtwに続くパルスPvの立上がりタイミングtvでクロックパルスCKのカウントを停止し、そのカウント値NについてもクロックパルスCKの周期τ[sec]を乗じてゼロクロス点tw-tv間の時間T[sec]を算出する。以下、周波数算出部22は、パルスPv,Pu,Pwの立上がりタイミングと立下りタイミングが順番に検出される毎にゼロクロス点tv-tu間,tu-tw間,tw-tv間,tv-tu間,tu-tw間,...の各時間T[sec]を算出し、各ゼロクロス点間の時間Tの6倍の時間の逆数を演算して電力系統の交流信号の周波数fsを算出し、その算出結果を周波数検出値として出力する。すなわち、周波数検出部22は、レベル比較部19からゼロクロスタイミングの検出信号が入力される毎に電力系統の交流信号の周波数fsを算出する。なお、ゼロクロス検出部20から入力されるパルスPuの立上がりタイミングtuと立下がりタイミングtuとの間のクロックパルスCKのパルス数3Nをカウントし、そのカウント値3NにクロックパルスCKの周期τ[sec]を乗じてパルスPuのハイレベル期間Tu[sec]を算出し、そのハイレベル期間Tuの2倍の期間の逆数を演算して電力系統の交流信号の周波数fs=1/2Tu=1/6Nτを算出するようにしてもよい。ゼロクロス検出部20から入力されるパルスPv,Pwについても同様である。また、例えば、パルスPuの立上がりタイミングtuから次の立上がりタイミングtuまでクロックパルスCKをカウントし、そのカウント値Nu=6NにクロックパルスCKの周期τ[sec]を乗じて電力系統の交流信号の周期2T[sec]を算出し、その算出値の逆数1/6Nτを演算して電力系統の交流信号の周波数fsを算出するようにしてもよい。パルスPv,Pwについても同様である。ゼロクロス点間カウント法は、三相交流電圧vu,vv,vwのレベルを所定の周期でサンプリングするため、サンプリング周期の間にゼロクロス点が入る場合は線形補間などによってそのゼロクロス点を検出する必要がある。また、三相交流電圧vu,vv,vwに不平衡成分や高調波成分が含まれる場合は電圧波形の歪みによりゼロクロス点の周期が正確に系統周波数fsに対応する周期にはならないので、周波数fsを数周期分にわたって算出し、それらの算出結果の平均値を取るなどして不平衡成分や高調波成分の悪影響を抑制する必要がある。第4実施形態に係るに係る周波数検出装置1Dでは、周波数算出部3Cに入力される三相交流電圧が基本波成分の三相交流電圧vsu,vsv,vswであるので、ゼロクロス点の周期は正確に系統周波数fsに対応する周期の6倍となっている。従って、第4実施形態に係るに係る周波数検出装置1Dによれば、ゼロクロス点とサンプリング点との誤差を無視し得る程度にサンプリング周波数を高くすることにより、従来のゼロクロス点間カウント法よりも高速かつ高精度で周波数fsを検出することができる。以上より、本発明に係る周波数検出装置1A,1B,1Cによれば、電力系統の電圧信号の周波数fsが急峻に変動した場合でも0.05秒以内にその変動に追従し、高速かつ高精度で電力系統の三相交流電圧の周波数fsを連続的に検出することができるという優れた効果を奏する。また、本発明に係る周波数検出装置1Dは、三相交流電圧のゼロクロス点間系統周波数fsに対応する周期の1/6の時間を要するので、従来の従来のゼロクロス点間カウント法による周波数検出装置に比べて高速かつ高精度で三相交流電圧の周波数fsを検出することかできるという優れた効果を奏する。なお、第1複素係数フィルタ部12には複素係数バンドパスフィルタ又は複素係数ノッチフィルタのいずれを用いてもよいが、好ましくは複素係数バンドパスフィルタよりも複素係数ノッチフィルタを用いたほうが高速かつ高精度の位相検出特性を得ることができる。また、複素係数ノッチフィルタと複素係数バンドパスフィルタとを組み合わせれば、両者の特性の相乗効果を期待することができ、より高速かつ高精度の位相検出特性を得ることができる。また、周知のように、複素係数ノッチフィルタ及び複素係数バンドパスフィルタを多段構成とすれば、急峻なフィルタ特性とすることができるとともに、不平衡成分や高調波成分の除去特性と応答性を容易に調整できるので、実装する場合は適当な段数の多段構成にするとよい。例えば、系統連系インバータを連系させる電力系統が2次の高調波成分を多く含む系統の場合は、2次高調波を除去する複素係数ノッチフィルタと複素係数バンドパスフィルタを組み合わせればよく、不平衡成分や高調波成分をあまり多く含まない電力系統であれば、応答速度の速いフィルタ構成にすればよい。第1実施形態に係る周波数検出装置1Aでは、第1複素係数フィルタ部12を複素係数バンドパスフィルタで構成する場合について説明したが、第3実施形態のように第1複素係数フィルタ部12を複素係数ノッチフィルタで構成してもよい。但し、この場合は、第3実施形態で説明したように、複素係数ノッチフィルタの阻止周波数を周波数算出部16で算出される周波数fに対して-f,-5f,7f,-11fの周波数に制御する必要がある。従って、この場合は、位相差算出部15の後段に正規角周波数の変化値Ωd’を算出する演算処理部を設け、その演算処理部で算出した正規角周波数の変化値Ωd’を複素係数ノッチフィルタで構成した第1複素係数フィルタ部12’にフィードバックさせるようにするとよい。また、第1実施形態の第1複素係数フィルタ部12を複素係数ノッチフィルタで構成する場合にはその複素係数ノッチフィルタに複素係数バンドパスフィルタを組み合わせた構成にして不平衡成分や高調波成分の除去効果を高めるようにしてもよい。この場合は、複素係数バンドパスフィルタの中心周波数を固定にしてもよいが、複素係数バンドパスフィルタの中心周波数を可変にし、複素係数ノッチフィルタの周波数を制御する機構例えば、上記の正規角周波数の変化値Ωd’を算出する演算処理部を利用して複素係数バンドパスフィルタの中心周波数を制御するようにしてもよい。また、第1実施形態の第1複素係数フィルタ部12を中心周波数が可変の複素係数バンドパスフィルタだけにしてもよい。この場合は、系統の周波数の変化に応じて複素係数バンドパスフィルタの中心周波数を変化させる構成を設ける必要がある。また、第2実施形態では、第1複素係数フィルタ部12の複素係数バンドパスフィルタの中心周波数を固定としたが、第2実施形態でも複素係数バンドパスフィルタの中心周波数を可変にし、電力系統の周波数の変化に応じて複素係数バンドパスフィルタの中心周波数を変化させる構成としてもよく、その複素係数バンドパスフィルタに第3実施形態の複素係数ノッチフィルタを組み合わせた構成としてもよい。これらの場合は、第2複素係数フィルタ部14’の中心周波数を制御する構成を利用して複素係数バンドパスフィルタの中心周波数や複素係数ノッチフィルタの阻止周波数を制御するようにすればよい。もちろん、中心周波数固定の複素係数バンドパスフィルタに複素係数ノッチフィルタを組み合わせた構成としてもよい。なお、第3実施形態で、第1複素係数フィルタ部12’を中心周波数固定の複素係数バンドパスフィルタで構成したり、中心周波数可変の複素係数バンドパスフィルタで構成したり、これらの複素係数バンドパスフィルタを複素係数ノッチフィルタに組み合わせた構成したりすることは、実質的に上述した第2実施形態の第1複素係数フィルタ部12の変形例と等価な構成となる。また、上記の説明では入力信号として三相交流電圧vu,vv,vwが入力される場合について説明したが、本発明に係る周波数検出装置1A~1Dは、単相交流電圧vが入力される場合にも適用することができる。この場合は、図1、図14、図16、図23に示すブロック図において、三相/二相変換部11、二相/三相変換部18を除いた構成にすればよい。単相の場合は、単相交流電圧vが1つしかないので、その単相交流電圧vのサンプリングデータが入力データの実数部urkとして第1複素係数フィルタ部12に入力され、入力データの虚数部ujkには0が入力される。なお、図1、図14、図16、図23に示す周波数検出装置1A~1Dにおいて、三相/二相変換部11と二相/三相変換部18を除去し、U,V,Wのいずれかの相の交流電圧vのサンプリングデータを入力データの実数部urkとして第1複素係数フィルタ部12に入力し、入力データの虚数部ujkには0を入力するようにしてもよい。複素係数フィルタを用いた第1複素係数フィルタ部12では、単相交流電圧が入力された場合でも三相交流電圧の場合と同様に互いに直交する二相交流電圧vr,vj正弦波と余弦波の信号が出力されるので、複素係数フィルタ部12,正規化部13及び周波数算出部1A,1B,1C,1Dは、図1、図14、図16、図23に示す三相用の周波数検出装置1A~1Dと同様の構成で実現することができる。従って、本発明は、簡単な構成で種々の系統連系インバータの周波数検出装置に広く適用することができる。"