交流インピーダンス測定データを、抵抗RとCPEConstantPhaseElementとが並列接続された回路ブロックを1つ以上有する等価回路モデルにフィッティングし、前記等価回路モデルの回路定数を求めるフィッティング部と、基準となる電池の交流インピーダンス測定データを、前記等価回路モデルにフィッティングして得られたCPE指数P値を保管するP値保管部と、電池の劣化度と前記等価回路モデルの回路定数との相関を記録したデータベースを参照し、判定対象電池の交流インピーダンス測定データを、前記等価回路モデルに前記P値を固定値としてフィッティングを行なって得られた回路定数に基づいて、前記判定対象電池の劣化判定を行なう劣化判定部と、を備えたことを特徴とする電池劣化判定装置。前記劣化判定部は、前記得られた回路定数によるフィッティング誤差があらかじめ設定した基準値以上の場合は、前記判定対象電池が使用不適当である旨の判定を行なうことを特徴とする請求項1に記載の電池劣化判定装置。前記フィッティング部によるフィッティングに先立ち、回路定数の初期値を設定するフィッティング初期値設定部をさらに備え、前記フィッティング初期値設定部は、交流インピーダンス測定データをナイキスト線図に表わした場合に、Zim=0となるゼロクロスポイントと、dZim/dZreの周波数特性線図の変曲点およびゼロを用いた楕円近似を利用して回路定数の初期値を設定することを特徴とする請求項1または2に記載の電池劣化判定装置。基準となる電池の交流インピーダンス測定データを、抵抗RとCPEConstantPhaseElementとが並列接続された回路ブロックを1つ以上有する等価回路モデルにフィッティングして得られたCPE指数P値を保管するP値保管ステップと、電池の劣化度と前記等価回路モデルの回路定数との相関を記録したデータベースを参照し、判定対象電池の交流インピーダンス測定データを前記等価回路モデルに前記P値を固定値としてフィッティングを行なって得られた回路定数に基づいて、前記判定対象電池の劣化判定を行なう劣化判定ステップと、を含むことを特徴とする電池劣化判定方法。既知の劣化度の電池の交流インピーダンス測定データを前記等価回路モデルに前記P値を固定値としてフィッティングを行なって得られた回路定数と前記既知の劣化度との相関を複数の劣化度について記録することで、前記データベースを作成するステップをさらに含むことを特徴とする請求項4に記載の電池劣化判定方法。一般に、RC並列回路ブロックを3段以上接続した等価回路を用いることで、つぶれた円弧状の複素インピーダンス特性に対してフィッティング誤差を小さくできる。しかし、その分フィッティングすべき未知の回路定数が増加することになる。また、RとCPEとの並列回路ブロックを1つ以上用いることで、つぶれた円弧状の交流インピーダンス特性に対してフィッティング誤差を小さくできることが知られている。ここで、CPEは、[数1]に示すようなインピーダンスZCPEで表わされる非線形要素であり、回路定数としてCPE指数のPとCPE定数のTとが用いられる。例えば、図11に示すようなRとCPEとの並列回路ブロックを2段接続した等価回路モデルを用いることで、図12に示すように、つぶれた円弧状の複素インピーダンス特性に対してフィッティング誤差を小さくできる。なお、図12において、破線は実際に測定された複素インピーダンス特性を示し、実線はRとCPEとの並列回路ブロックを2段接続した等価回路モデルにフィッティングして得られた回路定数による複素インピーダンス特性を示している。しかし、RとCPEとの並列回路ブロックは1段についてフィッティングすべき未知の回路定数はRとPとTの3つになり、2段モデルであっても、RC2段モデルと比較してフィッティングすべき未知の回路定数が増加する。フィッティングすべき未知の回路定数が多くなると、フィッティング処理が煩雑化し、得られる解が複数存在して劣化判定が不能となったり困難となる場合がある。このため、等価回路モデルの未知の回路定数を増やすことなくフィッティング誤差を小さくする技術が望まれている。そこで、本発明は、電池の劣化判定に用いる等価回路モデルの未知の回路定数を増やすことなくフィッティング誤差を小さくする技術を提供することを目的とする。上記課題を解決するため、本発明の第1の態様である電池劣化判定装置は、交流インピーダンス測定データを、抵抗RとCPEConstantPhaseElementとが並列接続された回路ブロックを1つ以上有する等価回路モデルにフィッティングし、前記等価回路モデルの回路定数を求めるフィッティング部と、基準となる電池の交流インピーダンス測定データを、前記等価回路モデルにフィッティングして得られたCPE指数P値を保管するP値保管部と、電池の劣化度と前記等価回路モデルの回路定数との相関を記録したデータベースを参照し、判定対象電池の交流インピーダンス測定データを、前記等価回路モデルに前記P値を固定値としてフィッティングを行なって得られた回路定数に基づいて、前記判定対象電池の劣化判定を行なう劣化判定部と、を備えたことを特徴とする。ここで、前記劣化判定部は、前記得られた回路定数によるフィッティング誤差があらかじめ設定した基準値以上の場合は、前記判定対象電池が使用不適当である旨の判定を行なうことができる。また、前記フィッティング部によるフィッティングに先立ち、回路定数の初期値を設定するフィッティング初期値設定部をさらに備え、前記フィッティング初期値設定部は、交流インピーダンス測定データをナイキスト線図に表わした場合に、Zim=0となるゼロクロスポイントと、dZim/dZreの周波数特性線図の変曲点およびゼロを用いた楕円近似を利用して回路定数の初期値を設定するようにしてもよい。上記課題を解決するため本発明の第2の態様である電池劣化判定方法は、基準となる電池の交流インピーダンス測定データを、抵抗RとCPEConstantPhaseElementとが並列接続された回路ブロックを1つ以上有する等価回路モデルにフィッティングして得られたCPE指数P値を保管するP値保管ステップと、電池の劣化度と前記等価回路モデルの回路定数との相関を記録したデータベースを参照し、判定対象電池の交流インピーダンス測定データを前記等価回路モデルに前記P値を固定値としてフィッティングを行なって得られた回路定数に基づいて、前記判定対象電池の劣化判定を行なう劣化判定ステップと、を含むことを特徴とする。ここで、既知の劣化度の電池の交流インピーダンス測定データを前記等価回路モデルに前記P値を固定値としてフィッティングを行なって得られた回路定数と前記既知の劣化度との相関を複数の劣化度について記録することで、前記データベースを作成するステップをさらに含めるようにしてもよい。本発明によれば、電池の劣化判定に用いる等価回路モデルの未知の回路定数を増やすことなくフィッティング誤差を小さくする技術が提供される。本発明の実施の形態について図面を参照して説明する。図1は、本実施形態に係る電池劣化判定システムの構成を示すブロック図である。本図に示すように電池劣化判定システム10は、電池劣化判定装置100と、回路定数-劣化度相関データベース200と、交流インピーダンス測定装置300を備えて構成される。電池劣化判定システム10は、これらを1台の装置として構成してもよいし、それぞれ独立した装置として構成してもよい。劣化判定対象の電池は、リチウムイオン電池等の二次電池が好適であるが、本発明は、燃料電池その他の電池の劣化判定に適用することができる。さらには、電池のSOCStateOfCharge判定に適用することも可能である。本実施形態においては、抵抗RとCPEConstantPhaseElementとが並列接続された回路ブロックを1つ以上有する等価回路モデルをフィッティングに用いるものとする。これにより、つぶれた円弧状の複素インピーダンス特性に対してフィッティング誤差を小さくする。具体的には、図11に示したように、抵抗R1と、抵抗R2とCPE2とが並列になった回路ブロックと、抵抗R3とCPE3とが並列になった回路ブロックと、CPE4と、抵抗R5とコイルL5とが並列になった回路ブロックとが接続された等価回路モデルを用いるものとする。ただし、本発明は、この等価回路モデルに限られず、抵抗RとCPEとが並列になった回路ブロックを1つ以上含んだ等価回路モデルであればよい。交流インピーダンス測定装置300は、測定対象の電池400の交流インピーダンスを測定し、測定結果を記録する装置であり、交流インピーダンス測定部310、交流インピーダンス測定条件設定部320、測定データ格納部330を備えている。交流インピーダンス測定部は、交流インピーダンス測定条件設定部320がユーザから受け付けた測定条件にしたがって、電池400に交流電圧あるいは交流電流を印加し、電流あるいは電圧を測定することで交流インピーダンスを算出する。測定は複数の周波数について行ない、記憶領域である測定データ格納部330に、例えば、ナイキスト線図データとして格納する。ただし、交流インピーダンスの測定方法は、上述の例に限られず種々の手法を用いることができる。例えば、複数の周波数信号が重畳された電圧波形あるいは電流波形を印加して、電流波形あるいは電圧波形を測定し、電圧波形、電流波形をそれぞれ離散フーリエ変換DFTして、周波数成分ごとの比を求めるようにしてもよい。また、交流インピーダンスを測定せずに別途取得したナイキスト線図データを測定データ格納部330に格納するようにしてもよい。回路定数-劣化度相関データベース200は、等価回路パラメータの回路定数と電池の劣化度容量劣化、出力劣化等との相関関係を記録するデータベースである。劣化判定対象の電池の交流インピーダンスを測定し、得られた複素インピーダンス特性を等価回路モデルにフィッティングすることにより回路定数を取得し、回路定数-劣化度相関データベース200を参照することで、電池の劣化判定を行なうことができる。電池劣化判定装置100は、コンピュータプログラムにしたがって動作を行なうPC等の情報処理装置を用いて構成することができ、制御部110、入出力部120、フィッティング初期値設定部130、フィッティング部140、P値保管部150、劣化判定部160を備えている。制御部110は、電池劣化判定装置100の各機能部を後述するような手順で動作させることにより、電池劣化判定処理を制御する。入出力部120は、回路定数-劣化度相関データベース200に対するデータ入出力、交流インピーダンス測定装置300の測定データ格納部330からのデータ読込を制御するとともに、ユーザから操作設定を受け付けたり、図示しない表示装置等に電池劣化判定結果等を出力する。フィッティング初期値設定部130は、測定データ格納部330から取得したナイキスト線図データから、所定のアルゴリズムに基づいて等価回路モデルの回路定数のフィッティング初期値を設定する。フィッティング部140は、フィッティング初期値設定部130が設定した初期値を用いて、等価回路モデルの回路定数のフィッティングを行なう。フィッティングアルゴリズムは既存の種々のアルゴリズムを採用することができる。P値保管部150は、判定対象の電池の劣化判定処理に先立ち、あらかじめ所定の条件で取得したP値を保管する記憶領域である。劣化判定部160は、回路定数-劣化度相関データベース200を参照して、劣化度測定対象の電池の複素インピーダンス特性を等価回路モデルにフィッティングして得られた回路定数から電池の劣化度を判定する。劣化度は、容量劣化量、出力劣化量等種々の態様で表現することができる。次に、本実施形態の電池劣化判定装置100の電池劣化判定動作について説明する。図2は、電池劣化判定動作の概要を示すフローチャートである。本実施形態では、まず、電池劣化判定用等価回路モデルのP値を設定するS10。すなわち、実際の電池劣化判定に先立ち、等価回路モデルのP値を設定しておき、設定されたP値を固定値として等価回路モデルのフィッティングを行なう。これにより、CPEの未知の回路定数はCPE定数Tのみとなるため、RとCPEとの並列回路ブロックを含んだ等価回路モデルであっても、未知の回路定数の数をCR並列回路ブロックと同数にすることができる。したがって、未知の回路定数を増やすことなくフィッティング誤差を小さくすることができる。ここで、あらかじめ設定されたP値を用いて、劣化判定対象電池のフィッティングを行なうようにしているのは、P値は、同種の電池であれば、電池としての特性が大きく損なわれない程度の劣化度、SOC、温度等にかかわらず、一定の値とみなせることが、実験的に確かめられたためである。次に、劣化度が既知の電池の回路定数を、設定されたP値を固定値としたフィッティングにより求めて、等価回路モデルの回路定数と劣化程度との相関データを作成するS20。相関データは、異なる劣化度について作成し、回路定数-劣化度相関データベース200に記録する。以降は、設定されたP値と、回路定数-劣化度相関データベース200に記録された相関データに基づいて、判定対象電池の劣化判定を繰り返すことができるS30。これらの処理の詳細な手順について説明する。まず、電池劣化判定用等価回路モデルのP値を設定する処理S10の詳細な手順について図3のフローチャートを参照して説明する。本処理では、P値を設定するための基準となる電池を用意し、解が複数になったり発散せずにP値を特定できる複素インピーダンスが得られる電池状態にコンディションニングするS101。解が複数になったり発散せずにP値を特定できる複素インピーダンス特性は、例えば、低周波領域方向で線が広がるような形状であり、一般に、所定の劣化状態や、温度が低い場合、SOCStateOfChargeが低い場合等に広がるような特性になる。上述のように、これらの条件が変化した場合であっても、P値は一定とみなせるため、フィッティングを容易に行なえる電池状態を任意に設定することができる。逆に、円弧が近接した複素インピーダンス特性は、フィッティングが発散する場合があるため、このような複素インピーダンス特性となる電池状態は避けるようにする。そして、交流インピーダンス測定装置300を用いて基準電池の交流インピーダンスを測定するS102。得られた測定データは、測定データ格納部330にナイキスト線図として格納する。次に、得られた基準電池のナイキスト線図について、フィッティング初期値設定部130を用いて、フィッティング初期値を設定するS103。フィッティング初期値は、例えば、以下のようなアルゴリズムで設定することができる。なお、基準電池の交流インピーダンスの測定の結果、図4に示すようなナイキスト線図が得られたものとする。R1初期値については、ナイキスト線図のゼロクロスポイントにおけるZreであるZre0を求め、Zre0をR1初期値とする。P2初期値、P3初期値については、図5に示すようなdZimf/DZrefの周波数特性を求め、周波数の高い方からdZimf/DZrefのゼロ、あるいは変曲点を探し、第1の周波数をf1、第2の周波数をf2とする。ただし、先に変曲点が現れた場合は、変曲点の周波数を使用し、ゼロは1回無視する。そして、[数2]に示すように、dZimf1とDZref1-Zre0との比から、P2初期値を求める。次いで、Zreが、Zref2-Zimf2×3/4、Zref2-Zimf2/4、Zref2Zimf2/4、Zref2Zimf2×3/4となる点のZimを求める。データ点がない場合は直線補間すればよい。求められた4つのデータ点から、長軸がZre方向と平行となる楕円を既存の方法で近似する。そして、求められた楕円の長軸長さをb、短軸長さをaとし、[数3]からP3初期値を求める。また、R2初期値は、R2=Zref1-Zre0×2で求め、T2初期値は、T2=1/[2πf1^P2×R2]で求め、R3初期値は、R3=bで求め、T3初期値は、T3=1/[2πf2^P3×R3]で求める。もちろん、他のアルゴリズムで各回路定数のフィッティング初期値を設定するようにしてもよい。図3の説明に戻って、フィッティング用初期値が設定されると、フィッティング部140が、設定された初期値を用いて等価回路モデルのフィッティングを行なうS104。低周波領域方向で広がるような形状のナイキスト線図が得られるように基準電池をコンディショニングしているため、既存のアルゴリズムを用いてフィッティングすることにより、各回路定数を容易に求めることができる。そして、フィッティングにより得られた回路定数のうち、P値P2、P3を、P値保管部150に格納しS105、電池劣化判定用等価回路モデルのP値を設定する処理S10を終了する。次に、等価回路モデルの回路定数と劣化程度の相関データを作成する処理S20の詳細な手順について図6のフローチャートを参照して説明する。まず、容量劣化や出力劣化等の劣化度が既知の電池を採択するS201。そして、採択された電池を劣化判定用のインピーダンス測定条件にコンディションニングしS202、交流インピーダンス測定装置300を用いて交流インピーダンスを測定するS203。劣化判定用のインピーダンス測定条件は、所定の温度、所定のSOC等とすることができる。得られた測定データは、測定データ格納部330にナイキスト線図として格納する。次に、得られた既知劣化度の電池のナイキスト線図について、フィッティング初期値設定部130を用いて、フィッティング初期値を設定するS204。フィッティング初期値は、上述の処理S103と同様のアルゴリズムで設定することができる。ただし、P値については固定値とするため初期値を算出する必要がない。そして、フィッティング部140が、P値保管部150に保管されているP値を固定の回路定数とし、設定された初期値を用いて等価回路モデルのフィッティングを行なうS205。P値を固定の回路定数としているため、未知の回路定数がRC2段モデルと同数である。したがって、既存のアルゴリズムを用いてフィッティングすることにより、各回路定数を容易に求めることができる。求められた回路定数は、電池の既知の劣化度と関連付けて回路定数-劣化度相関データベース200に相関データとして登録するS206。以上の処理を、異なる劣化度の電池に対して繰り返し行ないS207:No、登録された相関データの個数が充足するとS207:Yes、等価回路モデルの回路定数と劣化程度の相関データを作成する処理S20を終了する。図7は、回路定数R3と出力劣化量との相関例を示す図である。本図に示すように、等価回路モデルのRは、充放電サイクルを繰り返すことによって大きくなることが知られている。次に、判定対象電池の劣化判定処理S30の詳細な手順について図8のフローチャートを参照して説明する。まず、劣化判定対象電池を劣化判定用のインピーダンス測定条件にコンディションニングするS301。具体的には、温度、SOC等を、等価回路モデルの回路定数と劣化程度の相関データを作成する処理S20時と同等にコンディショニングする。そして、交流インピーダンス測定装置300を用いて判定対象電池の交流インピーダンスを測定するS302。得られた測定データは、測定データ格納部330にナイキスト線図として格納する。次に、得られた判定対象電池のナイキスト線図について、フィッティング初期値設定部130を用いて、フィッティング初期値を設定するS303。フィッティング初期値は、上述の処理S103と同様のアルゴリズムで設定することができる。ただし、P値については固定値とするため初期値を算出する必要がない。そして、フィッティング部140が、P値保管部150に保管されているP値を固定の回路定数とし、設定された初期値を用いて等価回路モデルのフィッティングを行なうS304。P値を固定の回路定数としているため、未知の回路定数がRC2段モデルと同数である。したがって、既存のアルゴリズムを用いてフィッティングすることにより、各回路定数を容易に求めることができる。そして、劣化判定部160が、回路定数-劣化度相関データベース200を参照して、求められた回路定数から判定対象電池の劣化度を判定するS305。判定結果は、図示しない表示装置や印刷装置等に出力しS306、判定対象電池の劣化判定処理S30を終了する。なお、本実施形態ではフィッティング誤差を小さくすることができるが、電池の劣化が進んで使用不能状態となると、P値を含んだ特性が変動し、フィッティング誤差が大きくなる。このため、劣化判定部160は、フィッティングの結果得られた回路定数を用いて算出された複素インピーダンス特性と、測定の結果得られた複素インピーダンス特性との2乗誤差を求め、2乗誤差があらかじめ設定した基準値以上となった場合には、その電池が使用不適当である旨の判定結果を、入出力部120を介して出力するようにしてもよい。以上説明したように、本実施形態の劣化度判定は、抵抗RとCPEとが並列になった回路ブロックを含んだ等価回路モデルを用いてフィッティングを行なっているため、フィッティング誤差を小さくすることができる。したがって、判定結果の妥当性を高めることができる。また、CPEの回路定数のうち、P値を固定値としているため、フィッティングすべき未知の回路定数をRCモデルと同等となり、容易にフィッティングを行なうことができる。