インピーダンスの周波数特性に基づいて燃料電池の特性を評価する燃料電池評価装置において、電流値および測定周波数を変化させながら燃料電池のターフェル領域におけるインピーダンスを取得するインピーダンス取得手段と、前記インピーダンス取得手段により取得された前記インピーダンスから反応抵抗を抽出する抽出手段と、前記抽出手段により抽出された前記反応抵抗と、前記インピーダンス取得手段により当該反応抵抗を含む当該インピーダンスを取得した際の電流値との積を算出する算出手段と、前記算出手段により算出された前記積を、前記測定周波数に対応する当該積の周波数特性として提示する提示手段と、を備えることを特徴とする燃料電池評価装置。前記抽出手段は、前記反応抵抗に起因しない直流抵抗の成分を排除することで前記反応抵抗を抽出することを特徴とする請求項1に記載の燃料電池評価装置。前記抽出手段は、分布定数域に起因する前記直流抵抗の成分を排除することを特徴とする請求項2に記載の燃料電池評価装置。前記抽出手段は、モデルフィッティングを用いて前記直流抵抗の成分を排除することを特徴とする請求項2または3に記載の燃料電池評価装置。前記インピーダンス取得手段、前記算出手段および前記提示手段は、複数の前記ターフェル領域にわたる電流値において測定されて得られる前記インピーダンスを取り扱うことを特徴とする請求項1~4のいずれか1項に記載の燃料電池評価装置。インピーダンスの周波数特性に基づいて燃料電池の特性を評価する燃料電池評価方法において、電流値および測定周波数を変化させながら燃料電池のターフェル領域におけるインピーダンスを取得するインピーダンス取得ステップと、前記インピーダンス取得ステップにより取得された前記インピーダンスから反応抵抗を抽出する抽出ステップと、前記抽出ステップにより抽出された前記反応抵抗と、前記インピーダンス取得ステップにより当該反応抵抗を含む当該インピーダンスを取得した際の電流値との積を算出する算出ステップと、前記算出ステップにより算出された前記積を、前記測定周波数に対応する当該積の周波数特性として提示する提示ステップと、を備えることを特徴とする燃料電池評価方法。電気化学特性に基づいて燃料電池の特性を評価する燃料電池評価装置において、前記燃料電池のターフェルプロットを取得するターフェルプロット取得手段と、前記ターフェルプロット取得手段により取得されたターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出する拡散領域抽出手段と、を備えることを特徴とする燃料電池評価装置。電流値および測定周波数を変化させながら前記拡散領域抽出手段により抽出された前記拡散領域におけるインピーダンスを取得するインピーダンス取得手段と、測定電流値を変化させつつ前記インピーダンス取得手段により取得されるインピーダンスと前記測定電流値との積の差分ΔZiと、前記測定電流値の差分Δiとの比ΔZi/Δiを拡散評価パラメータCとして算出する拡散評価パラメータ算出手段と、を備えることを特徴とする請求項7に記載の燃料電池評価装置。前記拡散領域抽出手段により抽出された前記拡散領域の下限電流値と前記インピーダンス取得手段により取得されたインピーダンスの測定電流値との差分Δ0iと、前記拡散評価パラメータCとを用いて、ターフェル領域におけるインピーダンスを前記下限電流値におけるインピーダンスとして推定する推定手段を備えることを特徴とする請求項8に記載の燃料電池評価装置。前記拡散評価パラメータ算出手段は、前記拡散評価パラメータCを理論式変形することで得られる周波数情報を持った限界拡散電流の評価パラメータを算出することを特徴とする8または9に記載の燃料電池評価装置。電気化学特性に基づいて燃料電池の特性を評価する燃料電池評価方法において、前記燃料電池のターフェルプロットを取得するターフェルプロット取得ステップと、前記ターフェルプロット取得ステップにより取得されたターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出する拡散領域抽出ステップと、をコンピュータが実行することを特徴とする燃料電池評価方法。電流値および測定周波数を変化させながら前記拡散領域抽出手段により抽出された前記拡散領域におけるインピーダンスを取得するインピーダンス取得ステップと、測定電流値を変化させつつ前記インピーダンス取得ステップにより取得されるインピーダンスと前記測定電流値との積の差分ΔZiと、前記測定電流値の差分Δiとの比ΔZi/Δiを拡散評価パラメータCとして算出する拡散評価パラメータ算出ステップと、をコンピュータが実行することを特徴とする請求項11に記載の燃料電池評価方法。前記拡散領域抽出ステップにより抽出された前記拡散領域の下限電流値と前記インピーダンス取得ステップにより取得されたインピーダンスの測定電流値との差分Δ0iと、前記拡散評価パラメータCとを用いて、ターフェル領域におけるインピーダンスを前記下限電流値におけるインピーダンスとして推定する推定ステップをコンピュータが実行することを特徴とする請求項12に記載の燃料電池評価方法。前記拡散評価パラメータ算出ステップでは、前記拡散評価パラメータCを理論式変形することで得られる周波数情報を持った限界拡散電流の評価パラメータを算出することを特徴とする12または13に記載の燃料電池評価方法。発電状態における燃料電池の特性を評価する燃料電池評価装置において、燃料電池の特性計測を行ないながら、電流値を変化させて燃料電池のインピーダンスを取得するインピーダンス取得手段と、前記インピーダンス取得手段により取得された前記インピーダンスと、当該インピーダンスを取得した電流密度との積である規格化インピーダンスを算出する算出手段と、を備えることを特徴とする燃料電池評価装置。前記算出手段は、前記規格化インピーダンスの値が一定となるそれぞれ異なる電流密度の2つの領域を確認することにより、前記燃料電池が安定発電状態であると判断することを特徴とする請求項15に記載の燃料電池評価装置。前記算出手段により算出された前記規格化インピーダンスを、前記燃料電池の発電状態の指標として提示する提示手段を備えることを特徴とする請求項15または16に記載の燃料電池評価装置。発電状態における燃料電池の特性を評価する燃料電池評価方法において、燃料電池の特性計測を行ないながら、電流値を変化させて燃料電池のインピーダンスを取得するインピーダンス取得ステップと、取得された前記インピーダンスと、当該インピーダンスを取得した電流密度との積である規格化インピーダンスを算出する算出ステップと、を有することを特徴とする燃料電池評価方法。前記規格化インピーダンスの値が一定となるそれぞれ異なる電流密度の2つの領域を確認することにより、前記燃料電池が安定発電状態であると判断するステップを有することを特徴とする請求項18に記載の燃料電池評価方法。一方、電極反応を周波数特性から総合的に評価するのに有効な手法としてインピーダンスがある。しかし、電極の反応抵抗Rctは電流と反比例の関係バトラーボルマーの式よりがあるため、測定動作点の電流値の変化に対してインピーダンスの大きさが敏感に変化する。そのため、領域b1、b2の反応過程のそれぞれについて区別してインピーダンスを評価することが必要となる。また、実際に測定されたインピーダンスには直流抵抗接触抵抗膜抵抗などや多孔質電極の特性を表す分布定数域のように反応抵抗に依存しないインピーダンスも含まれている場合がある。ターフェル勾配は、前述のとおり反応抵抗による特性であるため、ターフェル勾配とインピーダンスとの関係を獲得しようとする場合に、反応抵抗のみを抽出してからでなければ、その関係を精度よく求めることができないという事情がある。本発明の一つの目的は、インピーダンスを有用な特性値として用いて電極性能評価を行うことができる燃料電池評価装置および燃料電池評価方法を提供することにある。また、一般的に酸素還元反応は60mVdecと120mVdecの2つのターフェル勾配が存在すると考えられている。しかし、上記の001Acm2以下の領域はb1領域となっている可能性が高く、この技術によって求めたターフェル勾配b1を採用すると、実使用電流密度における活性化過電圧を過小評価してしまう可能性がある。そのため、適切に活性化過電圧を評価するには、2つ目のターフェル勾配を取得する必要がある。しかし、001Acm2以上の負荷領域では濃度過電圧の影響も無視できない可能性があり、どこから濃度過電圧の影響が効いているかを判断するのは難しく、その領域で精度よくターフェル勾配を求めるのは困難である。また、拡散抵抗は動作点の電流値で変化するため、従来から拡散抵抗を評価することは困難とされてきている。本発明の他の一つの目的は、拡散性能の評価を可能とするとともにターフェル領域における反応過程の評価にも寄与しうる燃料電池評価装置等を提供することにある。また、従来の燃料電池の特性評価、特に、触媒評価においては、より確からしい評価を行うために、再現性のある電流-電圧特性ターフェル曲線を取得する必要があった。しかし、ターフェル曲線は、前処理条件や測定の環境条件によって変化することが知られており、燃料電池を安定した発電状態でターフェル曲線を測定するためには前処理条件や測定の環境条件をそろえなければならず、そのため多くの手間、工数がかかっていた。本発明の他の一つの目的は、前処理条件や測定の環境条件を合わせることなく、燃料電池が安定した発電状態にあることを把握することができる燃料電池評価装置および燃料電池評価方法を提供することにある。本発明の燃料電池評価装置は、インピーダンスの周波数特性に基づいて燃料電池の特性を評価する燃料電池評価装置において、電流値および測定周波数を変化させながら燃料電池のターフェル領域におけるインピーダンスを取得するインピーダンス取得手段と、前記インピーダンス取得手段により取得された前記インピーダンスから反応抵抗を抽出する抽出手段と、前記抽出手段により抽出された前記反応抵抗と、前記インピーダンス取得手段により当該反応抵抗を含む当該インピーダンスを取得した際の電流値との積を算出する算出手段と、前記算出手段により算出された前記積を、前記測定周波数に対応する当該積の周波数特性として提示する提示手段と、を備えることを特徴とする。この燃料電池評価装置によれば、抽出された反応抵抗と電流値との積を算出するので、インピーダンスを有用な特性値として用いて電極性能評価を行うことができる。前記抽出手段は、前記反応抵抗に起因しない直流抵抗の成分を排除することで前記反応抵抗を抽出してもよい。前記抽出手段は、分布定数域に起因する前記直流抵抗の成分を排除してもよい。前記抽出手段は、モデルフィッティングを用いて前記直流抵抗の成分を排除してもよい。前記インピーダンス取得手段、前記算出手段および前記提示手段は、複数の前記ターフェル領域にわたる電流値において測定されて得られる前記インピーダンスを取り扱ってもよい。本発明の燃料電池評価方法は、インピーダンスの周波数特性に基づいて燃料電池の特性を評価する燃料電池評価方法において、電流値および測定周波数を変化させながら燃料電池のターフェル領域におけるインピーダンスを取得するインピーダンス取得ステップと、前記インピーダンス取得ステップにより取得された前記インピーダンスから反応抵抗を抽出する抽出ステップと、前記抽出ステップにより抽出された前記反応抵抗と、前記インピーダンス取得ステップにより当該反応抵抗を含む当該インピーダンスを取得した際の電流値との積を算出する算出ステップと、前記算出ステップにより算出された前記積を、前記測定周波数に対応する当該積の周波数特性として提示する提示ステップと、を備えることを特徴とする。この燃料電池評価方法によれば、抽出された反応抵抗と電流値との積を算出するので、インピーダンスを有用な特性値として用いて電極性能評価を行うことができる。本発明の燃料電池評価装置は、電気化学特性に基づいて燃料電池の特性を評価する燃料電池評価装置において、前記燃料電池のターフェルプロットを取得するターフェルプロット取得手段と、前記ターフェルプロット取得手段により取得されたターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出する拡散領域抽出手段と、を備えることを特徴とする。この燃料電池評価装置によれば、ターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出するので、拡散性能の評価とともに、ターフェル領域における反応過程の評価が可能となる。電流値および測定周波数を変化させながら前記拡散領域抽出手段により抽出された前記拡散領域におけるインピーダンスを取得するインピーダンス取得手段と、測定電流値を変化させつつ前記インピーダンス取得手段により取得されるインピーダンスと前記測定電流値との積の差分ΔZiと、前記測定電流値の差分Δiとの比ΔZi/Δiを拡散評価パラメータCとして算出する拡散評価パラメータ算出手段と、を備えてもよい。前記拡散領域抽出手段により抽出された前記拡散領域の下限電流値と前記インピーダンス取得手段により取得されたインピーダンスの測定電流値との差分Δ0iと、前記拡散評価パラメータCとを用いて、ターフェル領域におけるインピーダンスを前記下限電流値におけるインピーダンスとして推定する推定手段を備えてもよい。前記拡散評価パラメータ算出手段は、前記拡散評価パラメータCを理論式変形することで得られる周波数情報を持った限界拡散電流の評価パラメータを算出してもよい。本発明の燃料電池評価方法は、電気化学特性に基づいて燃料電池の特性を評価する燃料電池評価方法において、前記燃料電池のターフェルプロットを取得するターフェルプロット取得ステップと、前記ターフェルプロット取得ステップにより取得されたターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出する拡散領域抽出ステップと、をコンピュータが実行することを特徴とする。この燃料電池評価方法によれば、ターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出するので、拡散性能の評価とともに、ターフェル領域における反応過程の評価が可能となる。電流値および測定周波数を変化させながら前記拡散領域抽出手段により抽出された前記拡散領域におけるインピーダンスを取得するインピーダンス取得ステップと、測定電流値を変化させつつ前記インピーダンス取得ステップにより取得されるインピーダンスと前記測定電流値との積の差分ΔZiと、前記測定電流値の差分Δiとの比ΔZi/Δiを拡散評価パラメータCとして算出する拡散評価パラメータ算出ステップと、をコンピュータが実行してもよい。前記拡散領域抽出ステップにより抽出された前記拡散領域の下限電流値と前記インピーダンス取得ステップにより取得されたインピーダンスの測定電流値との差分Δ0iと、前記拡散評価パラメータCとを用いて、ターフェル領域におけるインピーダンスを前記下限電流値におけるインピーダンスとして推定する推定ステップをコンピュータが実行してもよい。前記拡散評価パラメータ算出ステップでは、前記拡散評価パラメータCを理論式変形することで得られる周波数情報を持った限界拡散電流の評価パラメータを算出してもよい。本発明の燃料電池評価装置は、発電状態における燃料電池の特性を評価する燃料電池評価装置において、燃料電池の特性計測を行ないながら、電流値を変化させて燃料電池のインピーダンスを取得するインピーダンス取得手段と、前記インピーダンス取得手段により取得された前記インピーダンスと、当該インピーダンスを取得した電流密度との積である規格化インピーダンスを算出する算出手段と、を備えることを特徴とする燃料電池評価装置。この燃料電池評価装置によれば、算出された規格化インピーダンスを燃料電池の発電状態の指標として使用できるため、前処理条件や測定の環境条件を合わせることなく、燃料電池の発電状態を把握することができる。また、規格化インピーダンスの情報を有用な解析データとして利用できる。前記算出手段は、前記規格化インピーダンスの値が一定となるそれぞれ異なる電流密度の2つの領域を確認することにより、前記燃料電池が安定発電状態であると判断してもよい。前記算出手段により算出された前記規格化インピーダンスを、前記燃料電池の発電状態の指標として提示する提示手段を備えてもよい。本発明の燃料電池評価方法は、発電状態における燃料電池の特性を評価する燃料電池評価方法において、燃料電池の特性計測を行ないながら、電流値を変化させて燃料電池のインピーダンスを取得するインピーダンス取得ステップと、取得された前記インピーダンスと、当該インピーダンスを取得した電流密度との積である規格化インピーダンスを算出する算出ステップと、を有することを特徴とする。この燃料電池評価方法によれば、算出された規格化インピーダンスを燃料電池の発電状態の指標として使用できるため、前処理条件や測定の環境条件を合わせることなく、燃料電池の発電状態を把握することができる。また、規格化インピーダンスの情報を有用な解析データとして利用できる。前記規格化インピーダンスの値が一定となるそれぞれ異なる電流密度の2つの領域を確認することにより、前記燃料電池が安定発電状態であると判断するステップを有してもよい。本発明の燃料電池評価装置によれば、抽出された反応抵抗と電流値との積を算出するので、インピーダンスを有用な特性値として用いて電極性能評価を行うことができる。本発明の燃料電池評価方法によれば、抽出された反応抵抗と電流値との積を算出するので、インピーダンスを有用な特性値として用いて電極性能評価を行うことができる。本発明の燃料電池評価装置によれば、ターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出するので、拡散性能の評価とともに、ターフェル領域における反応過程の評価が可能となる。本発明の燃料電池評価方法によれば、ターフェルプロットに基づき、ターフェル勾配が電流に比例する比例領域として拡散領域を抽出するので、拡散性能の評価とともに、ターフェル領域における反応過程の評価が可能となる。本発明の燃料電池評価装置によれば、算出された規格化インピーダンスを燃料電池の発電状態の指標として使用できるため、前処理条件や測定の環境条件を合わせることなく、燃料電池の発電状態を把握することができる。また、規格化インピーダンスの情報を有用な解析データとして利用できる。本発明の燃料電池評価方法によれば、前処理条件や測定の環境条件を合わせることなく、燃料電池の発電状態を把握することができる。また、規格化インピーダンスの情報を有用な解析データとして利用できる。以下、本発明の燃料電池評価装置の実施形態実施例1~3について説明する。図1は、実施例1の燃料電池評価装置としての燃料電池セルの測定システムを示すブロック図である。図1に示す燃料電池セル100は単セルであり、Pt/C触媒が表面に修飾された電解質膜101を、ガス拡散層をそれぞれ有するアノード極102およびカソード103で挟み込んだ構造をとる。またアノード102の側には、アノード102から分割された分割電極104が設けられている。アノード102および分割電極104は互いに絶縁された状態にある。測定時には、アノード102に水素、カソード103に空気もしくは酸素を供給し発電状態とする。電気化学測定は3電極系で行い、カソード特性を評価する際は作用極をカソード103、対極をアノード102とする。また、アノード極102およびカソード極103とは絶縁されたアノード側の別の分割電極104を参照極とする。図1に示すように、アノード102、カソード103および分割電極104は、それぞれポテンショスタット1に接続され、ポテンショスタット1による計測結果はレコーダ2に転送される。図1に示すように、ポテンショスタット1およびレコーダ2には、制御演算装置3が接続される。制御演算装置3はポテンショスタット1の制御、レコーダ2からの上記計測結果の取得、および、取得された上記計測結果に基づく演算を実行する。燃料電池セル100のアノード極102には水素が、カソード極103には空気もしくは酸素が供給される。これらの供給ガスは図1に示すガス供給装置4を介して与えられる。ガス供給装置4は制御演算装置3に接続され、燃料電池セル100に供給されるガスのガス濃度、流量、温度等は、ガス供給装置4を介して制御演算装置3により制御される。図1に示すように、制御演算装置3は、ポテンショスタット1を制御する負荷制御部31と、一定のガス供給量のもとで燃料電池セル100の電流値を変化させることによりターフェルプロットを取得するターフェルプロット取得部32と、各電流値の電流に交流を重畳させることで取得される測定データに基づいて燃料電池セル100のインピーダンスを取得するインピーダンス取得部33と、後述するターフェル領域の抽出を行うためのターフェル領域算出部34と、ガス供給装置3を制御するガス供給制御部35と、インピーダンス取得部33により取得されたインピーダンスに対して後述する規格化処理を施す規格化インピーダンス算出部36と、規格化インピーダンス算出部36における算出結果を提示する提示部37と、を備える。ポテンショスタット1を用いた負荷の制御は、電圧規制電流規制のどちらでもよい。電圧規制または電流規制は、負荷制御部31で指定の条件に従って制御され、レコーダ2で測定データが記録される。ターフェルプロットはターフェルプロット取得32によって得られ、ガス濃度の制御はガス供給制御部35およびガス供給装置4により導入ガスに窒素を混合することで行われる。後述するターフェルプロットの直流抵抗補正およびシフトによる濃度補正はターフェル領域算出部34により行われる。燃料電池セル100のインピーダンスはインピーダンス取得部33によって得られる。直流抵抗および分布定数域の補正を含む一連の規格化処理後述するステップS102~S105は規格化インピーダンス算出部36において行われ、周波数特性値として提示手段37によりモニタなどに提示される。また、レコーダ2に、生データのほか、ターフェル領域算出部34および規格化インピーダンス算出部36における算出結果を保存可能としてもよい。図2は、反応抵抗の抽出およびインピーダンスの規格化処理を示すフローチャートである。このフローチャートに示された手順に従って燃料電池セル100の反応抵抗を抽出し、規格化インピーダンスを求めることができる。まず、ステップS101において、ターフェル勾配が異なるb1、b2領域のそれぞれについてインピーダンスを測定する。ここでは、b1、b2領域を予め正しく抽出することが必要となるが、その手順については後述する。図3は、b1、b2領域における測定に基づきインピーダンス取得部33により取得されたインピーダンスを示す図である。図3において、Zi1は電流i1におけるインピーダンスを、Zi2は電流i2におけるインピーダンスを、Zi3は電流i3におけるインピーダンスを、Zi4は電流i4におけるインピーダンスを、それぞれ示す。また、電流i1および電流i2は領域b1に、電流i3および電流i4は領域b2に、それぞれ対応し、i1<i2<i3<i4である。ここでインピーダンスの実軸と交わる点をゼロクロスと定義する。この図3に示す例では、高周波側ゼロクロス値が0とはなっておらず、直流抵抗成分RDCが存在していることが分かる。この直流抵抗成分は反応抵抗とは直接関係しないので、ステップS102において、高周波でゼロクロスした値を直流抵抗値とし、その分を差し引いて原点補正を実施する。図4は、ステップS102において直流補正を施した後のインピーダンスを示す図である。図4において、高周波ゼロクロス値が原点となっていることが確認できる。また、図5は、図4に示すインピーダンス直流抵抗補正済のZf特性、θf特性を示す図である。次に、ステップS103では、分布定数域の判別を行う。図4および図5から判るように、高周波側ではインピーダンスに直線領域が現れていることが分かる。一般的に多孔質電極ではこのような45°の直線が高周波側で現れることが知られている分布定数域。分布定数域に関して成立する等式eq1、eq2、eq3を以下に示す。これらの等式から、分布定数域は反応抵抗Rctには依存していないことが分かる。したがって、分布定数域も除かなければ、後述する規格化によって精度良く反応を分離することができない。図5には、ステップS103において、分布定数域と判断する領域の例が示されている。図5の例では、触媒中のプロトン移動抵抗Rionや電気2重層容量Cdlが反応過程で異なっている可能性も許容し、ゼロクロス周波数からθ-f特性が45°±5°の領域グレー帯の部分かつ同じ反応過程で動作点が異なるインピーダンスZの比較において、低周波ゼロクロス値の差の5以内に該当する周波数までを分布定数域と判断している。ただし、触媒中のプロトン移動抵抗は異なる反応過程でも変化しないという思想や、直線は必ずしも45°とは限らないなどの思想のもと、分布定数域の具体的な判断基準は適宜変更が可能である。次に、ステップS104において、分布定数域を除いたモデルフィッティング処理を行う。図6は、ステップS104におけるモデルフィッティング処理の方法を示す図である。ここでは、フィッティングソフトを用いて、分布定数域以外の周波数帯、すなわち反応抵抗の影響が表れている低周波の周波数帯でのみモデルフィッティングを行う。図6では、Zi3についてのみ例示しており、Zi3についての分布定数域の上限周波数をfixとしている。ここで、分布特性伝送線モデルに起因している直流抵抗成分を総括してRporeと定義する。等価回路モデルとして、Rporeと反応インピーダンスモデルとを直列に接続した回路図6中に記載を使用する。反応インピーダンスモデルとして、実際の反応過程を表現していると考えられるモデルを用いることが望ましく、モデルは適宜選択できる。ここでは、簡易的にRCの2次モデルをZfとして用いてフィッティングを行った例を示す。このフィッティングによりRctとRporeの値が求まり、反応抵抗を反応円弧として簡易的に抽出することができる。図7は、分布定数域分を考慮し原点補正したインピーダンスを示す図である。図7では、図4に示すインピーダンスに対し、上記手順にて求めたRpore分を差し引き、原点を取り直している。ここで、ステップS102およびステップS104の処理は、いずれも直流抵抗分を取り除く処理であるため、ステップS104においてまとめてRporeRDCとして算出してもよい。次に、ステップS105では、直流抵抗分を取り除き反応抵抗のみを抽出したインピーダンスに対して規格化処理を施す。ここでは、原点補正したインピーダンス図7にインピーダンス測定動作点の電流値を乗じることでインピーダンスを規格化する。ここで、インピーダンスとインピーダンスを測定したときの測定電流値との関係を説明する。特に拡散の影響が現れない領域濃度非依存領域では電流Iと電圧Eの間にターフェルの式等式eq4の関係がある。等式eq4中、bはターフェル勾配である。一方、インピーダンス測定では入力電流の交流変調に対する出力電圧の応答に線形関係を仮定して等式eq5にてインピーダンスを算出する。等式eq4式より、濃度非依存領域においては等式eq6の関係がある。したがって、濃度非依存領域では、インピーダンスZに測定電流値Iを乗じた積IZは一定の値となる。この関係は例えば文献多孔質電極内の電流分布とインピーダンスの関係の解析光島重徳、神谷信行、太田健一郎ElectrochemistryP810814、2006に言及されている。この文献中では、周波数を変化させたときのIZ円弧のつぶれ具合を、イオン抵抗/電子抵抗の比を考察することに利用することが記載されている。図8は、抽出された反応抵抗に対して規格化処理して得たインピーダンス特性周波数特性を示す図である。このようなインピーダンス特性は、提示部37により表示される。図8に示すように、規格化インピーダンスが精度良く反応過程領域b1、b2ごとに分離できていることが確認できる。すなわち、Zi1とZi2に対応するプロット、およびZi3とZi4に対応するプロットが、それぞれ、ほぼ同一曲線上に重なり合うことが判る。したがって、インピーダンス特性に基づいて領域b1、b2の反応過程にあるのかが判定できるとともに、異なる電流値におけるインピーダンスを推定可能となる。また、この規格化インピーダンスの形大きさなどから反応過程それぞれの特性について論じることが可能となる。また、燃料電池の性能を多面的に示すパラメータとして広く利用できる。図9は、伝送線モデルを示す図である。このモデルでは、多数の直流抵抗とインピーダンスとを組み合わせることで、プロトン移動抵抗Rionおよび反応インピーダンスが表現される。上記実施形態では、ある周波数で分布定数域と集中定数域を区別し、集中定数域についてのみモデルフィッティングを施すことで簡易的に反応抵抗を抽出している。しかし、実際には分布定数域と集中定数域との境界領域に互いに影響し合っている遷移周波数帯が存在するため、明確に分布定数域と集中定数域を切り分けることはできない。これに対し、図9に示すような一般的な伝送線モデルを用いれば分布定数域遷移領域集中定数域のすべてを表現できるため、全周波数においてフィッティングが適用でき、分布定数域の判別プロセス図2のステップS103を踏むことなく、より正確に反応抵抗を抽出することが可能である。ただし、図9に示すモデルはパラメータ素子が多いため、的確な初期値を与えなければフィッティングの際、発散しやすいことに気をつけなければならない。このように、伝送線モデルを用いてフィッティングを行うことで、より正確に反応抵抗を抽出することができる。その結果、より精度よく反応抵抗を規格化することが可能となる。上記のように、ステップS101において、インピーダンスを測定するに際しては、予めb1、b2領域を正しく抽出することが必要となる。以下、b1、b2領域を抽出する手順について説明する。b1、b2領域の抽出に際し、ガス濃度の制御は、ガス供給制御部35からの指令に従ってガス供給装置4によって行われる。ガス濃度の制御については、圧力変化の影響を避けるため、カソード極103に導入されるガスの総流量を固定し、酸素と窒素の流量割合を変化させることで調節される。このとき、供給酸素不足による拡散性能低下の影響を避けるため、測定時の酸素濃度例えば、後述する酸素濃度C1、C2において利用率が十分低くなるように、あらかじめカソード103に導入されるガス総流量が決められている。すなわち、濃度を低く制御した際でも酸素の絶対量としては十分供給されている状態が確保される。補正勾配導出などの演算処理は制御演算装置3において行われる。図10は、制御演算装置3の動作手順を示すフローチャートである。図10のステップS1~ステップS4は、ターフェルプロットの取得および膜抵抗値の測定、取得のための処理を示している。図10のステップS1では、カソード極103に導入されるガスの酸素濃度を酸素濃度C1に制御し、燃料電池セル100に直流電流負荷I0DCを与え、そのときのセル電圧値ターフェルプロットを取得、保存する。次に、ステップS2では、直流電流負荷はそのままで、交流負荷をさらに重畳させることで、当該直流動作点におけるインピーダンスの測定を実施し、測定結果を保存する。図11は、交流の周波数を変化させながら燃料電池のインピーダンスを測定した場合の周波数特性を例示する図である。図11に示すように、一般的に燃料電池セルのインピーダンスは高周波側で抵抗値が低下しゼロクロスするが、実験的にこのゼロクロスする周波数はおよそ1kHzであることが判明している。このため、ステップS2におけるインピーダンス測定は、例えば1kHzの交流のみで実施するものとし、このときの実数値をその電流負荷における膜抵抗RrInと定義する。これにより、膜抵抗の測定時間を短縮できる。ただし、交流の周波数を変化させながらインピーダンス測定を行い、膜抵抗を正確に算出してもよい。次に、ステップS3では、現在の直流電流負荷IkDCにおけるセル電圧が06VvsRHEに到達したか否か判断し、判断が肯定されれば処理を終了し、判断が否定されればステップS1へ戻る。ステップS4では、直流電流負荷IkDCを所定値ΔIだけ増加させ、ステップS1へ戻る。このとき、Ik1DC-IkDC=ΔI>0とする。ここで、ΔIは電流ステップ直流電流の増加幅を意味し、1mA以下のできるだけ小さいステップを適用することが望ましい。このように、ステップS1~ステップS4の処理では、それぞれの電流負荷において、ターフェルプロットの取得と、膜抵抗値の測定、保存とをセル電圧が06VvsRHEとなるまで複数回繰り返す。ここでカットオフ電圧を06Vとしているのは、2つ目のターフェル勾配勾配b2がおよそ08V以下のカットオフ電圧で観測されるためである。ただし、このカットオフ電圧は状況に応じて燃料電池セル100が劣化しない範囲で下げることが可能である。次に、ステップS1~ステップS4の処理を酸素濃度C2C1C2の条件下において、同様に実施する。図10のステップS11~ステップS13は、ターフェル勾配の算出手順を示している。図10のステップS11では、1式に従い、ステップS2で得られた膜抵抗値からそれぞれの電流負荷における抵抗過電圧ηrInを算出する。ηrIn=InDC×RrIn1式さらに、ステップS11では、酸素濃度C1、C2のそれぞれについて、ステップS2で得られた抵抗過電圧を、ステップS1で得られたターフェルプロットから差し引くことで、補正後のターフェルプロットを算出する。次に、ステップS12では、ステップS11で得られた補正後のターフェルプロットを用いて、移動係数αCを求める。この処理の一例について以下に述べる。図12は、互いに異なる酸素濃度C1、C2でのターフェルプロットを用いて、移動係数αCを求める手法を示す図である。一般的に拡散の影響がないときバルク濃度=電極表面濃度のとき、2式および3式に示すように、ガス濃度に応じて平衡電位や交換電流が変化することが知られている。しかし、ガス濃度によって反応のメカニズム酸素還元反応が変わるわけではないので、酸素濃度C1、C2によってターフェル勾配は変化しない。2式、3式より、ガス濃度に起因するターフェルプロットのシフト変化の傾きは4式のように計算される。そこで、酸素濃度C2の場合の補正後ターフェルプロットに対し、4式を用いて傾き補正を施し、酸素濃度C1の場合の補正後ターフェルプロットに揃える図12。ここでαCは移動係数であり、一般的に計算上は05を用いる場合が多い。しかし、反応や実験条件により異なる場合が多いので、3式を用いて最も両者のプロットが揃うようなαCを最適化して決定する必要がある。この場合、酸素濃度C1、C2以外の濃度酸素濃度C3C4...で得られた補正後ターフェルプロットも使用して、αCのフィッティング精度を高めるようにしてもよい。次に、ステップS13では、ターフェル勾配b1、b2を算出する。この処理の手順については種々の変形が考えられるが、その一例について以下に述べる。図13は濃度非依存領域の決定手法の原理を示す図である。ここで仮に、拡散の影響がなく、活性化のみを観測しているとすれば、酸素濃度C1C2でのプロットは移動係数αCによる移動のみで一致し、また、プロットの全域において良い直線性を示すはずである。しかし、現実には拡散の影響は必ず出てしまい直線性は崩れる。そこで意図的に導入ガス濃度を変化させることにより、ガス拡散係数を変化させつつステップS1~ステップS4を実行してそれぞれのガス拡散係数での補正後ターフェルプロットを取得し、ガス拡散係数の変化に依らずプロットが一致する領域では拡散の影響を受けていないと判断することができる。この手法により決定された拡散の影響がない領域において直線近似を実施し、活性化過電圧に起因する2つのターフェル勾配b1、b2を適切に求めることができる。例えば、図13の例では、特定の電流値を閾値として、電流Iの小さな領域ではプロットはよく一致しているのに対し、電流Iの大きな領域ではプロットが乖離する。したがって、前者の領域においては図13のプロットがターフェル勾配b1、b2を直接的に示していると判断できる。この場合には、上記のようにこの領域において直線近似を実施すればよい。しかし、劣化による拡散性能が低下した場合や、ガス濃度が低すぎるなどの要因により拡散の影響が早く現れ、2つ目のターフェル勾配b2が観測できない場合、もしくは一致したターフェル勾配b2の領域が小さく線形近似の相関が低い場合には、直線近似によりターフェル勾配b2を求めることには無理がある。しかし、この場合には拡散の影響を考慮したターフェル補正式を用いることができる。図14は、比較的電流の小さな領域から、拡散の影響がターフェルプロットに現れている例を示す図である。ここで、拡散の影響が早く現れているということは、図14にも示されるように、限界拡散電流ILも比較的小さい電流負荷で観測が可能である。拡散の影響を考慮したターフェル補正式は5式で表され、右辺第1項は活性化過電圧ηctを、右辺第2項は濃度過電圧ηcを、それぞれ表している。そこで、5式を6式のように変形し、ターフェルプロットのx軸をLogIn/IL-Inにとり直したものが図15に示すプロットである。このように、図14のプロットを図15のプロットへと変換することにより、拡散の影響を考慮した上でターフェル勾配をグラフ上で表現することができる。この補正により、直線近似の相関を高めることが可能となり、精度よくターフェル勾配b2が求められる。このように、ステップS10~ステップS13の手順では、ガス濃度に起因するターフェルプロットのシフトを行い、また、理論式から傾き補正を施してプロットを揃え、プロットが一致した部分から濃度に影響されない領域を抽出している。このため、その領域において直線近似を行うことで適切なターフェル勾配を求めることができるとともに、各領域b1、b2の範囲を正確に抽出することができる。以上のように、本発明の燃料電池評価装置によれば、反応抵抗Rctに起因しないインピーダンスを除去し、反応インピーダンスのみを抽出するプロセスを踏んだ上で規格化処理を施すことで、精度良く反応過程ごとに規格化インピーダンスを分離することが可能となる。また、電流値に対して敏感に変化する反応インピーダンスを上記規格化処理によって領域b1、b2に対して、それぞれ精度良く1つにまとめて表現することができるので、測定動作点の変化の影響を受けずに1つの特性値として詳細な評価を行うことができる。さらに、燃料電池の評価のみならず、溶液系で多孔質電極を用いた場合など、酸素還元反応以外の電極反応においても、ターフェル式に従う反応抵抗を反応過程ごとに分離して周波数情報を持つ特性値反応インピーダンスとして表現し、評価することができる。上記実施例では、提示部37により示される特性として、コールコールプロットを用いる例を示したが、特性を示す方法は任意である。また、コールコールプロット中に、対応する測定周波数を示すようにしてもよい。また、コールコールプロットが描かれる平面に直行する周波数軸を設け、3次元座標上でのプロットとすることにより、測定周波数、実数成分、虚数成分の3つの値を読み取れるようにしてもよい。本発明の適用範囲は上記実施例に限定されることはない。本発明は、インピーダンスの周波数特性に基づいて燃料電池の特性を評価する燃料電池評価装置等に対し、広く適用することができる。図16は、実施例2の燃料電池評価装置としての燃料電池セルの測定システムを示すブロック図である。なお、図16において、図1と同一要素には同一符号を付している。図16に示す燃料電池セル100は単セルであり、Pt/C触媒が表面に修飾された電解質膜101を、ガス拡散層をそれぞれ有するアノード極102およびカソード103で挟み込んだ構造をとる。またアノード102の側には、アノード102から分割された分割電極104が設けられている。アノード102および分割電極104は互いに絶縁された状態にある。測定時には、アノード102に水素、カソード103に空気もしくは酸素を供給し発電状態とする。電気化学測定は3電極系で行い、カソード特性を評価する際は作用極をカソード103、対極をアノード102とする。また、アノード極102およびカソード極103とは絶縁されたアノード側の別の分割電極104を参照極とする。図16に示すように、アノード102、カソード103および分割電極104は、それぞれポテンショスタット1に接続され、ポテンショスタット1による計測結果はレコーダ2に転送される。図16に示すように、ポテンショスタット1およびレコーダ2には、制御演算装置3Aが接続される。制御演算装置30はポテンショスタット1の制御、レコーダ2からの上記計測結果の取得、および、取得された上記計測結果に基づく演算を実行する。燃料電池セル100のアノード極102には水素が、カソード極103には空気もしくは酸素が供給される。これらの供給ガスはガス供給装置不図示を介して与えられる。ガス供給装置は制御演算装置30に接続され、燃料電池セル100に供給されるガスのガス濃度、流量、温度等は、ガス供給装置を介して制御演算装置3Aにより制御される。図16に示すように、制御演算装置30は、ポテンショスタット1を制御する負荷制御部131と、一定のガス供給量のもとで燃料電池セル100の電流値を変化させることによりターフェルプロットを取得するターフェルプロット取得部132と、各電流値の電流に交流を重畳させることで取得される測定データに基づいて燃料電池セル100のインピーダンスを取得するインピーダンス取得部133と、後述する拡散領域の抽出を行うための拡散領域算出部134と、インピーダンス取得部133により取得されたインピーダンスに基づいて後述する拡散評価パラメータを算出する拡散評価パラメータ算出部135と、後述する規格化インピーダンスを算出するための規格化インピーダンス算出部136と、各部における算出結果を提示するための提示部137と、を備える。ポテンショスタット1を用いた負荷の制御は、電圧規制電流規制のどちらでもよい。電圧規制または電流規制は、負荷制御部131で指定の条件に従って制御され、レコーダ2で測定データが記録される。ターフェルプロットはターフェルプロット取得132によって得られ、拡散領域の算出は拡散領域算出部134により行われる。インピーダンスはインピーダンス取得部133によって得られる。拡散評価パラメータの算出は拡散評価パラメータ算出部135において行われ、拡散項の補正を施した規格化インピーダンスの算出は規格化インピーダンス算出部136において行われる。拡散領域算出部134、拡散評価パラメータ算出部135および規格化インピーダンス算出部136における算出結果は、提示部137によりモニタなどに提示される。また、レコーダ2に、生データのほか、拡散領域算出部134、評価パラメータ算出部135および規格化インピーダンス算出部136における算出結果等を保存してもよい。図17は、測定システムの動作を示すフローチャート、図18~図25は、提示部37による表示例を示す図である。図17のステップS201では、電圧掃引もしくは電流掃引によりIV特性を測定し、電流の対数と電圧との関係を表現するターフェルプロットを取得する。この処理は、ターフェルプロット取得部32の機能に相当する。図18は、拡散領域まで測定したターフェルプロットを示す図である。図18に示すように、酸素還元反応には一般的に2つのターフェル勾配が存在することが知られている。拡散領域では直線領域は存在せず、電流が大きくなるにつれ、勾配も変化する。ここで、反応物の拡散を考慮した電極反応式について述べる。分極が十分大きく、逆反応の電流が無視できる場合、電流iおよび過電圧ηは、等式eq11およびeq12のように表すことができる。さらに、等式eq12を次のように変形し、等式eq13を得る。上記のように、等式eq13の左辺はインピーダンスに電流値iを乗じたものとなっており、ターフェル勾配に相当する値である。インピーダンスに電流値を乗じる操作を本明細書では規格化処理と呼ぶ。限界拡散電流iLが測定電流iよりも十分大きい時、等式eq13のように近似することができる。この式から拡散の影響がでている領域では反応のターフェル勾配に拡散起因の電流に比例する項が加わっていることが分かる。ここで、インピーダンスとインピーダンスを測定したときの測定電流値との関係を説明する。特に拡散の影響が現れない領域濃度非依存領域では電流Iと電圧Eの間にターフェルの11式の関係がある。11式中、bはターフェル勾配である。一方、インピーダンス測定では入力電流の交流変調に対する出力電圧の応答に線形関係を仮定して12式にてインピーダンスを算出する。11式より、濃度非依存領域においては13式の関係がある。したがって、濃度非依存領域では、インピーダンスZに測定電流値Iを乗じた積IZは一定の値となる。この関係は例えば文献多孔質電極内の電流分布とインピーダンスの関係の解析光島重徳、神谷信行、太田健一郎ElectrochemistryP810814、2006に言及されている。この文献中では、周波数を変化させたときのIZ円弧のつぶれ具合を、イオン抵抗/電子抵抗の比を考察することに利用することが記載されている。次に、ステップS202では、拡散領域の抽出を行う。この処理は、拡散領域算出部34の機能に相当する。ここでは、まず、図18に示すターフェルプロットを各測定点で微分した値ターフェルプロットの傾きを算出する。図19aは、図18に示すターフェルプロットを各測定点で微分し、その結果の絶対値をターフェル勾配と定義して特性として表したグラフを示している。ターフェルプロット上の直線領域は、ここではプラトー領域として表され、b1、b2の領域を示す2つのプラトーが確認できる。また、等式eq13に示されるように、拡散の影響がでている領域b2を越える領域では電流に比例する項がターフェル勾配に加算され、比例領域を形成する。ただし、等式eq13は限界電流iLが測定電流iよりも十分大きい時に成立する近似式であるので、反応物質の濃度が高く、利用率が小さいほどこの比例領域は広くなる。図19bは、図19aのグラフと同一の特性を電流の対数-ターフェル勾配特性として表現したものを示す。図19では測定電流値がリニア軸表示のため、電流が小さい領域で現れるb1領域が、グラフ上で見えない場合がある。この場合、図20のように測定電流値を対数軸表記することでb1領域を判断し易くなり、例えば的確にインピーダンスの測定動作点を定めることが可能となる。図20は、b2の領域から拡散の影響がでている場合の電流-ターフェル特性を示す図である。図20に示す比例領域を探すことで拡散項を求めることができる。ここで、拡散の影響が出始めている電流をiC0と定義する。図20において、比例領域における拡散項に相当する部分を網掛けで示している。後述の図21においても同様である。図20に示すように、ターフェル勾配は電流iC0を境界として、b2の領域に相当するターフェル勾配に、電流に比例する拡散項が付加された値をとる。図21は、b1の領域から拡散の影響がでている場合の電流-ターフェル特性を示す図である。この場合にも、網掛けで示される比例領域を探すことで拡散項を求めることができる。図21に示すように、ターフェル勾配は電流ic0を境界として、b1の領域に相当するターフェル勾配およびb2の領域に相当するターフェル勾配に、それぞれ電流に比例する拡散項が付加された値をとる。拡散項により生ずる直線の傾きは、すべての領域で一定である。図22は、拡散補正を施した電流-ターフェル勾配特性を示す図である。拡散補正により、図20および図21の拡散項網掛けで示す部分が除去され、b1およびb2の領域を示すプラトーのみが現れる。なお、拡散補正は後述するステップS205における処理に相当する。次に、ステップS203では、インピーダンス測定および測定されたインピーダンスに対する規格化処理を行う。この処理は、規格化インピーダンス算出部136の機能に相当する。図23は、i1~i7の各動作点で測定した規格化インピーダンスZi1~Zi7を示す図である。図19~図22には、個々の規格化インピーダンスZi1~Zi7を測定した動作点が矢印↓により表されている。図19および図20では、規格化インピーダンスZi1~Zi2がb1領域に、規格化インピーダンスZi3~Zi4がb2領域に、規格化インピーダンスZi5~Zi7が領域b2および拡散領域に、それぞれ対応している。また、図21では、規格化インピーダンスZi1がb1の領域に、規格化インピーダンスZi2がb1領域および拡散領域に、規格化インピーダンスZi3~Zi7がb2領域および拡散領域に、それぞれ対応している。なお、図23は、拡散がb2領域の途中から現れる場合図19~図20に相当について示している。図23に示すように、b1の領域に属するインピーダンスおよびb2の領域に属するインピーダンスは、規格化によってそれぞれ共通の曲線群にまとまり、拡散の影響が加わると動作点電流が大きくなるにつれ規格化インピーダンスも大きくなっていることがわかる。すなわち、拡散の影響により、b2の領域の曲線群から離れた位置に曲線が移動する。次に、ステップS204では、拡散性能評価パラメータCを算出する。この処理は、拡散評価パラメータ算出部135の機能に相当する。ここで、拡散性能評価パラメータCの算出方法について述べる。等式eq14に規格化インピーダンスにおける評価パラメータの算出式を示す。この式は、拡散領域にある2点の規格化インピーダンスiZの差を、そのそれぞれの動作点電流値の差で除している。したがって図20あるいは図21に示す比例領域における直線の傾きを求めていることに相当する。すなわち電流に依存しない評価パラメータ値として抽出することができる。図24は、算出された拡散性能評価パラメータを例示する図である。このように拡散のみに起因する周波数情報を持った評価パラメータを求めることで、拡散性能についてより広い視点から議論することが可能となる。また、拡散性能評価パラメータを用いることで、実際に測定していない動作点におけるインピーダンスを推定することもできる。このような処理は推定手段の機能に相当する。次に、ステップS205では、拡散性能評価パラメータを用いて拡散項の補正を行う。この処理は、規格化インピーダンス算出部136の機能に相当する。等式eq15は、拡散性能評価パラメータを用いた拡散項の補正方法を示す補正式である。この式は評価パラメータCに測定動作点の電流値iXと、拡散影響が始まる電流値iC0の差を乗じたものを、元の規格化インピーダンスiZCXから差し引く補正を行うものである。すなわち拡散領域の任意の動作点の拡散項を求めて差し引くことで、元の規格化インピーダンスを補正している。図25は、拡散の補正を施した規格化インピーダンスを示す図である。図23と比較すれば明らかなように、上記手法によって拡散項の影響が排除され、拡散領域での測定結果を示す曲線群がb2の領域の規格化インピーダンスを示す位置に移動することがわかる。これにより、拡散性能が悪くb2領域が狭い場合やb2領域が見えない場合においても、的確にb2領域の規格化インピーダンスを算出し評価を行うことができる。また、等式eq14のCと等式eq13の右辺第2項の-iの係数とは等価直線領域の傾きに相当であるためeq16と表され、eq16を変形して等式eq17を得ることで、限界拡散電流iLについても周波数情報を持った特性値iLとして推測することができる。以上のように、本実施例の燃料電池評価装置によれば、ターフェルプロットの微分値ターフェル勾配を、電流-ターフェル勾配特性として表現し、プラトーおよび電流比例領域を探すことで、b1領域、b2領域、および拡散領域を簡便に分離することができるほか、拡散影響開始電流についても求めることができる。また、電流に依存しない周波数情報を持った拡散性能評価パラメータを算出することで、従来評価が困難であった拡散性能の広い視点からの議論が可能となった。また、限界拡散電流を、周波数情報をもった特性値として推測できることから、最適運転制御時の新たな情報として活用できる可能性がある。さらに、拡散性能評価パラメータを用いて規格化インピーダンスの拡散補正ができることから、拡散性能が悪くb2領域が小さい、もしくはb2領域が見えない電極反応においても、反応過程領域ごとに規格化インピーダンスを分離することができる。さらにまた、ターフェル式に従う反応と拡散の影響とが混合する電極反応において、拡散性能評価パラメータの算出および規格化インピーダンスの拡散補正ができるため、燃料電池の評価のみならず溶液系などの酸素還元反応以外の電極反応においても広く応用できる。本発明の適用範囲は上記実施例に限定されることはない。本発明は、電気化学特性に基づいて燃料電池の特性を評価する燃料電池評価装置等に対し、広く適用することができる。図26は、燃料電池セルアセンブリセルの評価装置を示す図であり、図26aは、燃料電池セルの特性計測時の燃料電池の断面および計測時の接続方法を示す図、図26bは計測評価装置の構成を機能的に示すブロック図である。図26aに示すように、燃料電池セル200は、Pt/C触媒が表面に修飾された電解質膜201と、電解質膜201と電気的に絶縁されたカソード電極202と、カソード電極202と同様に電解質膜201と電気的に絶縁され、電解質膜201を挟んでカソード電極202と対向するアノード電極203とを備える。カソード電極202は、電解質膜201の側から外側に向けて順次、触媒層、拡散層、およびセパレータ不図示を積層して構成される。同様に、アノード電極203は、電解質膜1の側から外側に向けて順次、触媒層、拡散層、およびセパレータ不図示を積層して構成される。カソード電極202およびアノード電極203を構成する上記セパレータには、それぞれ燃料ガスを供給するガス流路が形成されている。また、燃料電池セル200の負荷回路には、カソード電極202とアノード電極203の間の電圧燃料電池電圧Vfcを計測する電圧計測モジュール501と、電子負荷装置40に流れる電流負荷電流を計測する電流計測モジュール502とが接続される。電圧計測モジュール501および電流計測モジュール502による計測結果は、計測評価装置50に与えられる。図26aおよび図26bに示すように、計測評価装置50は、カソード電極202とアノード電極203の間に接続された電子負荷装置40の電流値を制御する負荷制御部51と、電流値を変化させながら燃料電池セル200のインピーダンスを取得するインピーダンス取得手段52と、インピーダンス取得手段52により得られたインピーダンスz’ohmと、インピーダンスを取得した際の電流密度iA/cm2の積である規格化インピーダンスZiA/cm2×z’ohmを算出する算出手段53と、算出手段53により算出された規格化インピーダンスZの値を、インピーダンスを取得した際の電流密度の対数logiに対するグラフで提示するモニタ等の提示手段54と、インピーダンスz’を取得するに際して燃料電池セル200に供給される燃料ガスを制御するガス供給制御部55と、により構成される。次に、規格化インピーダンス計測の手順について説明する。計測評価装置50の負荷制御部51は、電子負荷装置40に対して、インピーダンス計測を行う周波数、直流電流値、重畳交流電流振幅を設定し、電池への電流負荷を制御する。交流成分が重畳された電流負荷に対して、燃料電池電圧Vfcの波形を電圧モジュール501で計測し、インピーダンス取得手段52において、交流電流負荷に対する電圧燃料電池電圧Vfcの交流成分のゲインおよび位相フェーズに基づき、インピーダンスz’ohmを求める。そして、算出手段53は、求めたインピーダンスz’ohmに対し、インピーダンスを取得した電流密度iA/cm2を積算し、規格化インピーダンスZiA/cm2×z’ohmを算出計測する。図27aは、燃料電池セル10の電圧-電流特性を例示する図、図27bは、燃料電池セル200の規格化インピーダンスZ-電流特性を例示する図である。ここで、電流を示す横軸は電流密度の対数logiとされている。計測評価装置50は、負荷制御部51により負荷電流値を徐々に増加させながら、個々の電流値におけるインピーダンスを計測して、規格化インピーダンスZを算出する。このとき、燃料電池セル200は、時間の経過と共に安定発電状態となる。燃料電池セル200が安定発電状態になると、拡散の影響が現れない領域濃度非依存領域である第1ターフェル領域61および第2ターフェル領域62と、拡散の影響が現れる拡散領域濃度依存領域63と、が確認できる。図27bに示すように、規格化インピーダンスZ-電流特性を示すグラフでは、第1ターフェル領域61および第2ターフェル領域62が、それぞれ一定の値の規格化インピーダンスZ1およびZ2をとる領域として認識でき、視覚的にも容易に確認することができる。図27bに示す例では、第1ターフェル領域61では規格化インピーダンスZ1が確認でき、第2ターフェル領域62では規格化インピーダンスZ2が確認できる。つまり、規格化インピーダンスZ-電流特性において、規格化インピーダンスZ1と規格化インピーダンスZ2が確認できれば、燃料電池セル200が安定発電状態となったと判断することができ、再現性のある燃料電池の特性計測ここでは、図27aに示す電圧-電流特性ができたことを知ることができる。図27cは、第1ターフェル領域61および第2ターフェル領域62における規格化インピーダンスZのコールコールプロットを示す図である。図27cに示すコールコールプロット71は、第1ターフェル領域61におけるコールコールプロットであり、図27cに示すコールコールプロット72は、第1ターフェル領域62におけるコールコールプロットである。コールコールプロット71における実数軸の切片座標およびコールコールプロット72における実数軸の切片座標が、それぞれ規格化インピーダンスZ1およびZ2に相当する。また、インピーダンスの規格化により、コールコールプロットに現れる反応円弧は第1ターフェル領域61および第2ターフェル領域でそれぞれ一定形状をとる。なお、第1ターフェル領域61と第2ターフェル領域62との間の遷移領域では、図27cにおけるコールコールプロット74として示すように、反応円弧の形状も両者の中間的な形状となる。このようなコールコールプロットは、インピーダンス取得手段52を介して得られた広い範囲の周波数における計測値に基づいて算出手段53において算出される。しかし、本実施例においては、インピーダンス計測の手順を簡素化することもでき、規格化インピーダンスZ1およびZ2をとる領域が判定できる範囲で、測定周波数を限定してもよい。例えば、既定の単数あるいは複数の周波数でのみ、インピーダンス測定を行ってもよい。図28は、計測評価装置50の動作を示すフローチャートである。計測評価装置50は、負荷電流値を増加させて燃料電池セル200の電流-電圧特性を計測しながら、インピーダンスを計測し規格化インピーダンスZを算出する。まず、計測条件の設定を行うステップS301。例えば、負荷制御部51が電子負荷装置40の負荷電流値の範囲を設定し、インピーダンス取得手段52が膜抵抗値を計測する。これにより、測定されるインピーダンスから膜抵抗値を除外することができる。次に、負荷制御部51は、電子負荷装置4の負荷電流値をステップS301で設定した負荷電流値の初期値ILに設定するステップS302。続いて、インピーダンス取得手段52は、燃料電池セル200のインピーダンスを計測し、算出手段53は規格化インピーダンスZを算出するステップS303。算出された規格化インピーダンスZは、規格化インピーダンスZ-電流特性としてプロットされる。次に、負荷電流値がステップS301で定めた負荷電流値の範囲の最大値imaxに到達したか否かを判断するステップS304。ステップS304の判断が肯定されればステップS306へ進み、判断が否定されればステップS305へ進む。ステップS305では、負荷制御部51は、電子負荷装置40の負荷電流値を一定値だけ増加させ、ステップS303へ戻る。一方、ステップS306では、算出手段53は、ステップS303で算出された規格化インピーダンスZに基づき、規格化インピーダンスZ1が存在するか否かを判定する。ステップS306の判定の結果、規格化インピーダンスZ1が存在しない場合には、ステップS302へ戻る。ステップS304の判定の結果、規格化インピーダンスZ1が存在する場合には、規格化インピーダンスZ1を記憶部に保存する。なお、すでに規格化インピーダンスZ1が保存されている場合には、これを新たに得られた規格化インピーダンスZ1に更新する。次に、算出手段53は、ステップS303で算出された規格化インピーダンスZに基づき、規格化インピーダンスZ2が存在するか否かを判定するステップS308。ステップS308の判定の結果、規格化インピーダンスZ2が存在しない場合には、ステップS302へ戻る。ステップS308の判定の結果、規格化インピーダンスZ2が存在する場合には、規格化インピーダンスZ2を記憶部に保存し、処理を終了する。このように、ステップS302~S309では、規格化インピーダンスZの値が一定となるそれぞれ異なる電流密度の2つの領域、すなわち、図27bに示す第1ターフェル領域61および第2ターフェル領域62が検出されるまで、インピーダンス計測および規格化Zの算出が繰り返される。そして、第1ターフェル領域61および第2ターフェル領域62が検出された状態においてインピーダンス取得手段52において取得されたインピーダンスを、安定発電状態でのインピーダンスであると判断することができる。インピーダンス取得手段52において取得されたインピーダンスは、適宜、計測評価装置50に保存され、燃料電池セル200の評価に用いることができる。なお、上記の手順で安定発電状態が確認された後、あらためてインピーダンス測定、その他の測定を実行することで、安定発電状態における燃料電池セル200について種々の評価を行うことができる。図29は、第1ターフェル領域61および第2ターフェル領域62が検出される以前の、発電状態が安定しない段階で得られる電流-電圧特性および規格化インピーダンスを例示する図である。図29aに示す例では、燃料電池セル200は未だ安定発電状態でないため、電流-電圧特性および規格化インピーダンスとしては拡散領域のものしか現れず、規格化インピーダンスZ1と規格化インピーダンスZ2は両者とも確認できない。図29bはさらに時間が経過した状態を示しており、第1ターフェル領域61の規格化インピーダンスZ1が確認できる。しかし、この段階でも安定した発電状態ではないため、規格化インピーダンスZ2を確認できない。このように、最初に規格化インピーダンスZ1が確認され、さらに安定発電状態に近づくと、規格化インピーダンスZ2が確認できる。燃料電池セル200が安定発電状態となると、規格化インピーダンスZ1と規格化インピーダンスZ2が確認できる図27b。このように、燃料電池セル200の特性計測をしながら、燃料電池セル200のインピーダンスを計測し、インピーダンスと電流密度の積である規格化インピーダンスZを算出して燃料電池セル200の発電状態の指標として使用することにより、燃料電池セル200の発電状態を確認、把握することができる。すなわち、このような規格化インピーダンスZのプロットに基づき、算出手段53における演算によって燃料電池セル200が安定発電状態か否かの判別を自動的に実行することができる。また、電流-電圧特性や規格化インピーダンスを示すグラフ等を提示手段54により表示することで、燃料電池セル200の特性や測定中の状況等を視覚的に確認することもできる。例えば、規格化インピーダンスZを、当該インピーダンスを取得した電流密度の対数に対応付けてプロットしたグラフを表示することで、規格化インピーダンスZの値が一定となるそれぞれ異なる電流密度の2つの領域規格化インピーダンスZ1、Z2を視認でき、容易に燃料電池セル200の安定発電状態を確認、把握することができる。一般に、燃料電池セル200の特性評価のためには、再現性のある電圧-電流特性ターフェル曲線を求めなければならない。しかし、再現性のある電圧-電流特性ターフェル曲線は、測定の環境条件又は前処理条件等によって変化することが知られており、再現性を得るためには測定条件を揃える必要がある。しかし、本実施例の燃料電池評価装置によれば、規格化インピーダンスZのプロットを用いることにより、測定の環境条件又は前処理条件等を合わせる手間なく、容易に燃料電池セル200の発電状態を確認できる。また、規格化インピーダンスZ1、Z2は安定発電状態におけるパラメータとして広く利用できる。上記実施例では、規格化インピーダンスZ1、Z2と燃料電池セル200の安定発電状態との関係について言及しているが、燃料電池の性能解析のため、規格化インピーダンスZの情報を幅広く適用することが可能である。上記実施例では、燃料電池の特性計測として、電流-電圧特性を計測したが、安定発電状態において計測が必要な特性評価であれば、広く適用できる。以上説明したように、本実施例の燃料電池評価装置によれば、インピーダンス取得手段52は、燃料電池の特性計測を行ないながら、電流値を変化させて燃料電池セル200のインピーダンスを取得し、算出手段53は、インピーダンス取得手段52により取得されたインピーダンスと、インピーダンスを取得した電流密度との積である規格化インピーダンスを算出し、提示手段54は、算出手段53が算出した規格化インピーダンスを燃料電池セル200の安定発電状態信号として提示する。これにより、測定の環境条件又は前処理条件等を合わせる手間なく、容易に燃料電池セル200の発電状態の指標として提示することができる。また、規格化インピーダンスの情報を有用な解析データとして利用できる。本発明の適用範囲は上記実施例に限定されることはない。本発明は、燃料電池セルの発電状態における特性を評価する燃料電池評価装置等に対し、広く適用することができる。