酸化剤ガスと燃料ガスとを電気化学反応させて電気エネルギを出力する燃料電池1に適用され、前記燃料電池1の内部を流れる前記酸化剤ガスのガス濃度を測定するガス濃度測定装置であって、前記燃料電池1の出力信号に所定の測定周波数fmおよび前記測定周波数fmよりも高い周波数に設定された高周波数fhを合成した交流信号を印加する交流信号印加手段3と、前記燃料電池1の出力電流を検出する電流検出手段51と、前記燃料電池1の出力電圧を検出する電圧検出手段52と、前記交流信号印加手段3にて前記交流信号を印加した際の前記出力電流および前記出力電圧に基づいて、前記測定周波数fmに対応する内部インピーダンスZおよび前記高周波数に対応する高周波インピーダンスZhを算出するインピーダンス算出手段53と、前記内部インピーダンスZから前記高周波インピーダンスZhを減算することで、前記内部インピーダンスZから前記酸化剤ガスのガス濃度の変化に応じて変化する特定インピーダンスZ ́を算出する特定インピーダンス算出手段54aと、前記特定インピーダンスZ ́の位相角θを算出する位相角算出手段54bと、前記特定インピーダンスZ ́の位相角θに基づいて前記酸化剤ガスのガス濃度を算出するガス濃度算出手段54cと、を備えることを特徴とするガス濃度測定装置。前記測定周波数fmは、前記酸化剤ガスのガス濃度を予め基準ガス濃度とした条件において、前記交流信号印加手段3にて印加する周波数を高周波から低周波まで変化させた際に、前記内部インピーダンスZの実数部および虚数部がなす半円状の軌跡の頂部Pに達したときの周波数を設定することを特徴とする請求項1に記載のガス濃度測定装置。前記燃料電池1の負荷を検出する負荷検出手段53aを備え、前記交流信号印加手段3は、前記燃料電池1の出力信号に、前記測定周波数fm、前記高周波数fh、および前記測定周波数fmと異なる複数の補正用周波数fm1、fm2を合成した交流信号を印加し、前記インピーダンス算出手段53は、前記負荷検出手段53aにて検出された負荷が所定の第1負荷基準値よりも大きい場合に、前記複数の補正用周波数fm1、fm2のうち前記測定周波数fmよりも低い補正用周波数におけるインピーダンスを前記内部インピーダンスZとして算出し、前記負荷検出手段53aにて検出された負荷が前記第1負荷基準値よりも低い値に設定された第2負荷基準値よりも小さい場合に、前記複数の補正用周波数fm1、fm2のうち前記測定周波数fmよりも高い補正用周波数におけるインピーダンスを前記内部インピーダンスZとして算出することを特徴とする請求項1または2に記載のガス濃度測定装置。前記電流検出手段51にて検出した出力電流から直流成分を抽出すると共に、予め前記直流成分と前記出力電圧との関係を関連付けた電流電圧特性に基づいて、前記内部インピーダンスZに含まれる反応抵抗R2を算出する反応抵抗算出手段を備え、前記特定インピーダンス算出手段54aは、前記反応抵抗算出手段により算出された前記反応抵抗R2の影響を除去すると共に、前記内部インピーダンスZから前記高周波インピーダンスZhを減算することで前記特定インピーダンスZ ́を算出することを特徴とする請求項1または2に記載のガス濃度測定装置。前記反応抵抗算出手段は、前記直流成分が予め定められた閾値よりも大きい場合、前記直流成分の増加に伴って前記出力電圧が直線的に減少する前記電流電圧特性に基づいて前記反応抵抗R2を算出し、前記直流成分が前記閾値以下の場合、前記直流成分の増加に伴って前記出力電圧が指数関数的に減少する前記電流電圧特性に基づいて前記反応抵抗R2を算出することを特徴とする請求項4に記載のガス濃度測定装置。前記燃料電池1の温度Tを検出する温度検出手段46を備え、前記交流信号印加手段3は、前記燃料電池1の出力信号に、前記測定周波数fm、前記高周波数fh、および前記測定周波数fmよりも周波数が低い補正用周波数fm3を合成した交流信号を印加し、前記インピーダンス算出手段53は、前記温度検出手段46にて検出された温度Tが所定の基準温度Toよりも低い場合に、前記補正用周波数fm3におけるインピーダンスを前記内部インピーダンスZとして算出することを特徴とする請求項1または2に記載のガス濃度測定装置。前記特定インピーダンスZ ́の位相θおよび前記酸化剤ガスのガス濃度を予め関連付けた制御特性を記憶する記憶手段54dを備え、前記ガス濃度算出手段54cは、前記記憶手段54dに記憶された前記制御特性を参照して、前記位相角算出手段54bにて算出した前記特定インピーダンスZ ́の位相θに対応するガス濃度を算出することを特徴とする請求項1ないし6のいずれか1つに記載のガス濃度測定装置。しかしながら、特許文献1では、燃料電池に印加する交流電流の周波数帯域を特定の帯域に限定しているものの、等価回路の回路定数を推定する場合には、当該特定の帯域において多数の周波数毎に内部インピーダンスを測定する必要があるので、ガス濃度の測定に時間例えば、数秒程度がかかってしまう。この場合、例えば、数ミリ秒~数十ミリ秒といった短時間で酸化剤ガスのガス濃度が変動する車両用の燃料電池システム等では、ガス濃度の測定を適切に行うことができないといった問題がある。本発明は上記点に鑑みて、燃料電池の内部を流れる酸化剤ガスのガス濃度を測定するガス濃度測定装置において、ガス濃度の測定時間の短縮化を図ることを目的とする。上記目的を達成するために、本発明では、燃料電池の内部における酸化剤ガスのガス濃度と燃料電池の出力信号に所定の周波数を印加した際の燃料電池1の内部インピーダンスに含まれる特定のインピーダンス成分の位相角との間の相関関係を利用することで、長時間を要する燃料電池1の等価回路の回路定数を算出することなく、酸化剤ガスのガス濃度を算出する構成としている。すなわち、請求項1に記載の発明では、酸化剤ガスと燃料ガスとを電気化学反応させて電気エネルギを出力する燃料電池1に適用され、燃料電池1の内部を流れる酸化剤ガスのガス濃度を測定するガス濃度測定装置であって、燃料電池1の出力信号に所定の測定周波数fmおよび測定周波数fmよりも高い周波数に設定された高周波数fhを合成した交流信号を印加する交流信号印加手段3と、燃料電池1の出力電流を検出する電流検出手段51と、燃料電池1の出力電圧を検出する電圧検出手段52と、交流信号印加手段3にて交流信号を印加した際の出力電流および出力電圧に基づいて、測定周波数fmに対応する内部インピーダンスZおよび高周波数に対応する高周波インピーダンスZhを算出するインピーダンス算出手段53と、内部インピーダンスZから高周波インピーダンスZhを減算することで、内部インピーダンスZから酸化剤ガスのガス濃度の変化に応じて変化する特定インピーダンスZ ́を算出する特定インピーダンス算出手段54aと、特定インピーダンスZ ́の位相角θを算出する位相角算出手段54bと、特定インピーダンスZ ́の位相角θに基づいて酸化剤ガスのガス濃度を算出するガス濃度算出手段54cと、を備えることを特徴とする。このように、酸化剤ガスのガス濃度および酸化剤ガスのガス濃度の変化に応じて変化する特定インピーダンスZ ́の位相角θとの相関関係を利用して、酸化剤ガスのガス濃度の変化を測定する構成とすれば、燃料電池1の等価回路の回路定数を算出することなく、酸化剤ガスのガス濃度を短時間で算出可能となる。この結果、酸化剤ガスのガス濃度が短時間で変動する燃料電池システムに適用する場合であっても、酸化剤ガスのガス濃度を適切に測定することが可能となる。具体的には、請求項2に記載の発明のように、請求項1に記載のガス濃度測定装置において、測定周波数fmとして、酸化剤ガスのガス濃度を予め基準ガス濃度とした条件において、交流信号印加手段3にて印加する周波数を高周波から低周波まで変化させた際に、内部インピーダンスZの実数部および虚数部がなす半円状の軌跡の頂部Pに達したときの周波数を設定することができる。ここで、燃料電池1の内部インピーダンスZ特定インピーダンスZ ́の位相角θは、酸化剤ガスのガス濃度の変化以外の要因によっても変化してしまうことがある。例えば、燃料電池1の負荷が大きい場合には、内部インピーダンスZが増大して特定インピーダンスZ ́の位相角θが大きくなり、燃料電池1の負荷が小さい場合には、内部インピーダンスZが減少して特定インピーダンスZ ́の位相角θが小さくなるといった傾向がある。すなわち、燃料電池1の負荷変動によって、特定インピーダンスZ ́の位相角θが変化してしまうことがある。そこで、請求項3に記載の発明では、燃料電池1の負荷を検出する負荷検出手段53aを備え、交流信号印加手段3は、燃料電池1の出力信号に、測定周波数fm、高周波数fh、および測定周波数fmと異なる複数の補正用周波数fm1、fm2を合成した交流信号を印加し、インピーダンス算出手段53は、負荷検出手段53aにて検出された負荷が所定の第1負荷基準値よりも大きい場合に、複数の補正用周波数fm1、fm2のうち測定周波数fmよりも低い補正用周波数におけるインピーダンスを内部インピーダンスZとして算出し、負荷検出手段53aにて検出された負荷が第1負荷基準値よりも低い値に設定された第2負荷基準値よりも小さい場合に、複数の補正用周波数fm1、fm2のうち測定周波数fmよりも高い補正用周波数におけるインピーダンスを内部インピーダンスZとして算出することを特徴とする。このように、燃料電池1の負荷変動に応じて、内部インピーダンスZを算出する際の周波数を変更することで、燃料電池1の負荷変動による特定インピーダンスZ ́の位相角θの変化を補正することができる。この結果、酸化剤ガスのガス濃度をより適切に測定することが可能となる。上述の如く、燃料電池1の負荷変動によって特定インピーダンスZ ́の位相角θが変動することがあるが、本発明者らの検討によると、この特定インピーダンスZ ́の位相角θの変動は、燃料電池1の負荷変動に伴う内部インピーダンスZに含まれる反応抵抗R2の変化が要因であることがわかっている。そこで、請求項4に記載の発明では、請求項1または2に記載のガス濃度測定装置において、電流検出手段51にて検出した出力電流から直流成分を抽出すると共に、予め直流成分と出力電圧との関係を関連付けた電流電圧特性に基づいて、内部インピーダンスZに含まれる反応抵抗R2を算出する反応抵抗算出手段を備え、特定インピーダンス算出手段54aは、反応抵抗算出手段により算出された反応抵抗R2の影響を除去すると共に、内部インピーダンスZから高周波インピーダンスZhを減算することで、特定インピーダンスZ ́を算出することを特徴とする。これによると、内部インピーダンスZから燃料電池1の負荷変動によって変動する反応抵抗R2を除去するので、燃料電池1の負荷変動によって特定インピーダンスZ ́の位相角θが変動してしまうことを抑制することができる。この結果、酸化剤ガスのガス濃度をより適切に測定することが可能となる。具体的には、請求項5に記載の発明のように、請求項4に記載のガス濃度測定装置において、反応抵抗算出手段にて、直流成分が予め定められた閾値よりも大きい場合、直流成分の増加に伴って出力電圧が直線的に減少する電流電圧特性に基づいて反応抵抗R2を算出し、直流成分が閾値以下の場合、直流成分の増加に伴って出力電圧が指数関数的に減少する電流電圧特性に基づいて反応抵抗R2を算出する構成とすることができる。また、燃料電池1の温度Tが低い場合には、内部インピーダンスZが増大して特定インピーダンスZ ́の位相角θが大きくなるといった傾向がある。このため、請求項6に記載の発明では、請求項1または2に記載のガス濃度測定装置において、請求項1または2に記載のガス濃度測定装置において、燃料電池1の温度Tを検出する温度検出手段46を備え、交流信号印加手段3は、燃料電池1の出力信号に、測定周波数fm、高周波数fh、および測定周波数fmよりも周波数が低い補正用周波数fm3を合成した交流信号を印加し、インピーダンス算出手段53は、温度検出手段46にて検出された温度Tが所定の基準温度Toよりも低い場合に、補正用周波数fm3におけるインピーダンスを内部インピーダンスZとして算出することを特徴とする。このように、燃料電池1の温度が低い場合には、内部インピーダンスZを算出する際の周波数を低くすることで、燃料電池1の温度が低いことによる特定インピーダンスZ ́の位相角θの変化を補正することができる。この結果、酸化剤ガスのガス濃度をより適切に測定することが可能となる。具体的には、請求項7に記載の発明のように、請求項1ないし6のいずれか1つに記載のガス濃度測定装置において、特定インピーダンスZ ́の位相θおよび酸化剤ガスのガス濃度を予め関連付けた制御特性を記憶する記憶手段54dを備え、ガス濃度算出手段54cでは、記憶手段54dに記憶された制御特性を参照して、位相算出手段54bにて算出した特定インピーダンスZ ́の位相θに対応するガス濃度を算出するようにしてもよい。なお、この欄および特許請求の範囲で記載した各手段の括弧内の符号は、後述する実施形態に記載の具体的手段との対応関係を示すものである。以下、本発明の実施形態について図に基づいて説明する。なお、以下の各実施形態相互において、互いに同一もしくは均等である部分には、図中、同一符号を付してある。第1実施形態以下、本発明の第1実施形態について図1~図8に基づいて説明する。図1は、本実施形態のガス濃度測定装置50を適用した燃料電池システムの全体構成図であり、図2は、本実施形態の電池セルの模式図である。この燃料電池システムは、電気自動車の一種である、いわゆる燃料電池車両に適用されており、車両走行用電動モータ等の電気負荷に電力を供給するものである。まず、燃料電池システムは、図1に示すように、水素と酸素との電気化学反応を利用して電力を発生する燃料電池1を備えている。燃料電池1は、図示しない車両走行用電動モータや2次電池といった各種電気負荷2に供給される電気エネルギを出力するもので、本実施形態では、固体高分子電解質型燃料電池を採用している。より具体的には、燃料電池1は、基本単位となる電池セル10以下、単にセル10と記載する。が複数個、電気的に直列に接続されて構成されたものである。換言すれば、燃料電池1は、複数のセル10が積層配置されて構成されている。図2に示すように、各セル10は、固体高分子からなる電解質膜101aの両側面に一対の電極101b、101cが配置された膜電極接合体MEA:MembraneElectrodeAssembly101と、この膜電極接合体101を狭持する一対のセパレータ102、103で構成されている。一対のセパレータ102、103は、カーボン材や導電性金属よりなる板状プレートからなり、アノード電極101bと対向する面に水素が流れる水素流路図示略が形成され、カソード電極101cと対向する面に空気が流れる空気流路図示略が形成されている。各セル10では、以下に示すように、水素と酸素とを電気化学反応させて、電気エネルギを出力する。負極側:アノード電極H2→2H2e-正極側:カソード電極2H1/2O22e-→H2O図1に戻り、燃料電池1と電気負荷2との間には、双方向に電力を伝達可能なDC-DCコンバータ3を介して電気的に接続されている。このDC-DCコンバータ3は、燃料電池1から電気負荷2、あるいは電気負荷2から燃料電池1への電力の流れを制御するものである。また、各セル10のうち、空気酸化剤ガスのガス濃度以下、単に酸素濃度と記載する。の測定対象となるセル10には、ガス濃度測定装置50が接続されている。このガス濃度測定装置50については後述する。燃料電池1のカソード電極101c側には、酸化剤ガスである空気酸素を燃料電池1に供給するための空気供給配管20、並びに、燃料電池1にて電気化学反応を終えた余剰空気および空気極で生成された生成水を燃料電池1から外気へ排出するための空気排出配管21が接続されている。空気供給配管20の最上流部には、大気中から吸入した空気を燃料電池1に圧送するための空気ポンプ22が設けられ、空気排出配管21には、燃料電池1内の空気の圧力を調整するための空気調圧弁23が設けられている。なお、本実施形態では、空気ポンプ22および空気調圧弁23によって、所定の流量および圧力の空気を燃料電池1に供給する酸化剤ガス側のガス供給手段が構成される。燃料電池1のアノード電極101b側には、燃料ガスである水素を燃料電池1に供給するための水素供給配管30、アノード側に溜まった生成水を微量な水素とともに燃料電池1から外気へ排出するための水素排出配管31が接続されている。水素供給配管30の最上流部には、高圧水素が充填された高圧水素タンク32が設けられ、水素供給配管30における高圧水素タンク32と燃料電池1との間には、燃料電池1に供給される水素の圧力を調整する水素調圧弁33が設けられている。なお、本実施形態では、この水素調圧弁33によって、所定の圧力の水素を燃料電池1に供給する燃料ガス側のガス供給手段が構成される。水素排出配管31には、生成水を微量な水素とともに外気へ排出するために所定の時間間隔で開閉する電磁弁34が設けられている。なお、上述の電気化学反応では、アノード電極101b側において生成水は発生しないものの、アノード電極101b側には、カソード電極101c側から各セル10の電解質膜101aを透過した生成水が溜まるおそれがある。このため、本実施形態では、水素排出配管31および電磁弁34を設けている。ところで、燃料電池1は発電効率確保のために運転中一定温度例えば80°C程度に維持する必要がある。このため、燃料電池1には、燃料電池1を冷却するための冷却水回路40が接続されている。この冷却水回路40には、燃料電池1に冷却水熱媒体を循環させるウォータポンプ41、電動ファン42を備えたラジエータ放熱器43が設けられている。さらに、冷却水回路40には、冷却水を、ラジエータ43を迂回するように流すバイパス流路44が設けられている。冷却水回路40とバイパス流路44との合流点には、バイパス流路44に流れる冷却水流量を調整するための流路切替弁45が設けられている。この流路切替弁45の弁開度が調整されることによって、冷却水回路40の冷却能力が調整される。また、冷却水回路40の燃料電池1の出口側付近には、燃料電池1から流出した冷却水の温度を検出する温度検出手段としての温度センサ46が設けられている。この温度センサ46により冷却水温度を検出することで、燃料電池1の温度を間接的に検出することができる。なお、この温度センサ46の検出信号は、後述する制御装置60に入力される。燃料電池システムには、各種制御を行う発電制御手段としての制御装置ECU60が設けられている。この制御装置60は、入力信号に基づいて、燃料電池システムを構成する各種電気式アクチュエータの作動を制御するもので、CPU、ROM、RAM等からなる周知のマイクロコンピュータとその周辺回路にて構成されている。具体的には、制御装置60の入力側には、ガス濃度測定装置50、温度センサ46の検出信号等以外に、車室内に設けられた車両起動スイッチ60aの操作信号等が入力される。なお、車両起動スイッチ60aは、空気ポンプ22、空気調圧弁23、水素調圧弁33等の作動開始信号を出力する開始信号出力手段の機能を兼ねる。一方、出力側には、上述の空気ポンプ22、空気調圧弁23、水素調圧弁33、電磁弁34、ウォータポンプ41、流路切替弁45等の各種電気式アクチュエータ等が接続されている。次に、本実施形態のガス濃度測定装置50について図3に基づいて説明する。図3は、本実施形態のガス濃度測定装置50の模式図である。本実施形態のガス濃度測定装置50は、測定対象となるセル10測定対象セルにおける空気出口部図2参照付近の酸素濃度を測定するもので、図3に示すように、電流測定装置51、電圧センサ52、上述のDC-DCコンバータ3、信号処理回路53、演算装置54を備えて構成されている。電流測定装置51は、セル面内の局所電流を検出する局所電流測定装置であって、酸素濃度を測定するセル10における空気出口部付近から出力される出力電流を検出する電流検出手段を構成している。本実施形態の電流測定装置51は、測定対象となるセル10に隣接して配置された板状部材を備え、当該セル10を流れる局所電流を検出可能に構成されている。電流測定装置51としては、例えば、板状部材におけるセル10における測定対象部位に対応する局所部位に溝部を形成し、当該溝部に電流センサを配置して構成することができる。なお、電流センサは、シャント抵抗、ホールIC等を利用した周知のセンサを用いることができる。電圧センサ52は、セル10から出力される出力電圧空気出口部付近の一対の電極101b、101c間の電位差を検出する電圧検出手段である。電流測定装置51および電圧センサ52それぞれは、信号処理回路53に接続されており、電流測定装置51および電圧センサ52の各検出信号が信号処理回路53に出力されるように構成されている。DC-DCコンバータ3は、上述のように、電気負荷2および燃料電池1の双方向に電力を伝達可能なものであり、本実施形態では、セル10の出力信号に任意の周波数で、正弦波等の交流信号交流電流を印加する交流信号印加手段を構成している。なお、本実施形態のDC-DCコンバータ3は、異なる複数の周波数を合成した交流信号を燃料電池1の出力信号に印加可能に構成されている。信号処理回路53は、DC-DCコンバータ3にて任意の周波数の交流信号交流電流が印加された際の電流測定装置51および電圧センサ52から出力される各検出信号に基づいて、セル10の内部インピーダンスZを算出するインピーダンス算出手段を構成している。信号処理回路53は、演算装置54に接続されており、信号処理回路53にて算出した内部インピーダンスZの算出結果が演算装置54に出力されるように構成されている。さらに、本実施形態の信号処理回路53は、異なる複数の周波数の交流信号交流電流が印加された場合に、高速フーリエ変換処理等によって、交流信号に含まれる直流成分直流電流および各周波数に対応する交流成分を抽出可能に構成されている。ここで、交流信号に含まれる直流成分直流電流は、燃料電池1の負荷変動に相関関係があるので、信号処理回路53にて交流信号に含まれる直流成分を抽出することで、間接的に燃料電池1の負荷変動を検出することができる。なお、信号処理回路53における交流信号に含まれる直流成分を抽出する構成ハードウェア、ソフトウェアを負荷検出手段53aとする。演算装置54は、各種入力信号に基づいて各種演算処理を実行するもので、CPU、およびROM、RAMといった記憶部記憶手段54d等からなる周知のマイクロコンピュータとその周辺回路により構成されている。例えば、本実施形態の演算装置54は、内部インピーダンスZのうち、セル10内の酸素濃度の変化に応じて変化する特定インピーダンスZ ́を算出する算出処理、当該特定インピーダンスZ ́の位相角θを算出する算出処理、および特定インピーダンスZ ́の位相角θからセル10内の酸素濃度を算出する算出処理を行う。なお、演算装置54にて算出する特定インピーダンスZ ́は、内部インピーダンスZから電解質膜101aの膜抵抗の影響を除去したインピーダンスである。本実施形態では、演算装置54における特定インピーダンスZ ́を算出する構成ハードウェアおよびソフトウェアを特定インピーダンス算出手段54a、特定インピーダンスZ ́の位相角を算出する構成ハードウェアおよびソフトウェアを位相角算出手段54b、セル10内の酸素濃度を算出する構成ハードウェアおよびソフトウェアをガス濃度算出手段54cとする。ここで、セル10の交流インピーダンスZに含まれる特定インピーダンスZ ́の位相角θに基づいてセル10内の酸素濃度を算出する算出方法について、図4~図7に基づいて説明する。図4は、燃料電池1のカソード側電極101cの等価回路を示す回路図である。図4の等価回路におけるR1は、電解質膜101aの膜抵抗、R2は反応抵抗、R3はカソード側過電圧、Cdは電気二重層コンデンサ成分に相当している。図5は、図4に示す等価回路において、高周波から低周波までの交流信号交流電流を印加した場合の燃料電池の内部インピーダンスZを複素平面上に示した特性図Cole-Coleプロット図である。図5に示すように、セル10に印加する交流信号の周波数が無限に大きい場合f=∞の内部インピーダンスZ∞は、電解質膜101aの膜抵抗R1となる。一方、セル10に印加する交流信号の周波数が非常に小さい場合f≒0の内部インピーダンスZoは、膜抵抗R1、反応抵抗R2、カソード側過電圧R3の合算R1R2R3となる。なお、高周波から低周波まで周波数を変化させたときの内部インピーダンスZは、図5に示すように、反応抵抗R2、カソード側過電圧R3の合算値を直径とした半円状の軌跡となる。ところで、燃料電池1は、セル10内の酸素濃度が低下すると、反応抵抗R2およびカソード側過電圧R3が増大する。このため、セル10内の酸素濃度が低下すると、Cole-Coleプロット図における半円状の軌跡の直径R2R3が大きくなる。なお、内部インピーダンスZの実数部から電解質膜101aの膜抵抗R1を減算したインピーダンスが、内部インピーダンスZにおける酸素濃度の変化に応じて変化する特定インピーダンスZ ́となる。図6は、セル10に供給する空気のストイキ比Stを変化させた場合のCole-Coleプロット図を示している。ここで、空気のストイキ比Stとは、燃料電池1に所定値の電流を流す際に必要となる空気の理論供給量と実際の空気の供給量との比を意味する。ストイキ比Stが小さい場合には、セル10内の酸素濃度が低く、ストイキ比Stが大きい場合には、セル10内の酸素濃度が高くなるといった関係がある。なお、図6では、空気のストイキ比Stが1.0の場合の内部インピーダンスZを黒丸、空気のストイキ比Stが1.2の場合の内部インピーダンスZを黒四角、空気のストイキ比Stが1.4の場合の内部インピーダンスZを黒三角、空気のストイキ比Stが1.8の場合の内部インピーダンスZをバツ印で示している。図6に示すように、空気のストイキ比Stが小さい場合には、セル10内の酸素濃度が低下し、空気のストイキ比Stが大きい場合に比べて、反応抵抗R2およびカソード側過電圧R3が増大する。このため、セル10内の酸素濃度の低下に伴って、Cole-Coleプロット図における半円状の軌跡の直径R2R3が相似的に大きくなる。そして、所定の測定周波数fmの交流信号を印加した場合の内部インピーダンスZは、空気のストイキ比Stが小さい場合に、ストイキ比Stが大きい場合に比べて、内部インピーダンスZの虚数部リアクタンス成分が増大する。例えば、ストイキ比Stが1.4の場合に、Cole-Coleプロット図における半円状の軌跡の頂部虚数部のピーク値Pに達する周波数を測定周波数fm20Hzとすると、当該測定周波数fmにおける内部インピーダンスZは、図6の白丸および各白丸を結ぶ一点鎖線で示すように、ストイキ比Stが小さくなるほど内部インピーダンスZの虚数部が増大する。なお、このような関係は、ストイキ比Stが1.1の場合に、Cole-Coleプロット図における半円状の軌跡の頂部虚数部のピークPに達する周波数を測定周波数fm5Hzとしても同様の関係となる図6の白三角および各白三角を結ぶ一点鎖線参照。ここで、特定インピーダンスZ ́の位相角θは、図6に示すように、Cole-Coleプロット図において、セル10に印加する交流信号の周波数が無限に大きい場合f=∞の内部インピーダンスZ∞と、測定周波数の交流信号を印加した際の内部インピーダンスZmとを結ぶ仮想線Lと実軸とのなす角度鋭角に相当している。この特定インピーダンスZ ́の位相角θは、図7の特定インピーダンスZ ́の位相角θとストイキ比Stとの関係を示す特性図に示すように、ストイキ比Stが小さくなるほど大きくなる。つまり、特定インピーダンスZ ́の位相角θとセル10内の酸素濃度との間には、セル10内の酸素濃度の低下に伴って、特定インピーダンスZ ́の位相角θが増大する位相が遅れるといった相関関係がある。このような特定インピーダンスZ ́の位相角θとセル10内の酸素濃度の相関関係を利用すれば、特定インピーダンスZ ́の位相角θを算出することで、セル10内の酸素濃度を算出することが可能となる。上記構成に係る本実施形態のガス濃度測定装置50でのセル10内の酸素濃度を測定する際の作動を、図8のフローチャートにより説明する。図8に示す制御フローは、車両起動スイッチ60aが投入ONされて、燃料電池1の発電状態となるとスタートする。なお、図8では、信号処理回路53および演算装置54にて行う処理を1つのフローチャートで示している。車両が起動すると、フラグ、タイマ等の初期化処理がなされるS10。そして、DC-DCコンバータ3から燃料電池1に所定の交流信号を印加するS20。DC-DCコンバータ3では、予め設定された測定周波数fm、および当該測定周波数fmよりも高い周波数に設定された高周波数fhを合成した交流信号を印加する。具体的には、測定周波数fmとしては、セル10内の酸素濃度を予め設定された基準ガス濃度例えば、ストイキ比=1.4としたときのセル内の酸素濃度とした条件で、印加する周波数を高周波から低周波まで変化させた際に、Cole-Coleプロット図における半円状の軌跡の頂部虚数部のピーク値Pに達したときの周波数を設定することができる。また、高周波数fhとしては、印加する周波数を低周波から高周波まで変化させた際に、内部インピーダンスZの虚数部が0付近に収束したときの周波数を設定することができる。次に、信号処理回路53および演算装置54に出力される各種信号を読み込む。具体的には、電流測定装置51および電圧センサ52にて検出したセル10の出力電流および出力電圧を読み込むS30。次に、S30にて検出した出力電流および出力電圧における測定周波数fmに対応する交流成分を抽出し、抽出した交流成分に基づいて、燃料電池1の内部インピーダンスZを算出するS40。また、S30にて検出した出力電流および出力電圧における高周波数fhに対応する交流成分を抽出し、抽出した交流成分に基づいて、燃料電池1の高周波インピーダンスZhを算出するS50。なお、本実施形態の高周波インピーダンスZhは、セル10に印加する交流信号の周波数が無限に大きい場合f=∞の内部インピーダンスZ∞に相当する。次に、S40で算出した燃料電池1の内部インピーダンスZからS50で算出した高周波インピーダンスZhを減算することで、セル10内の酸素濃度の変化に応じて変化する特定インピーダンスZ ́を算出するS60。具体的には、内部インピーダンスZの実数部から高周波インピーダンスZhの実数部を減算し、内部インピーダンスZから電解質膜101aの膜抵抗R1の影響を除去することで、特定インピーダンスZ ́を算出する。次に、S60にて算出した特定インピーダンスZ ́の位相角θを算出するS70。具体的には、この位相角θは、特定インピーダンスZ ́の実数部をZre ́、虚数部をZim ́としたときに次に示す数式F1を用いて算出することができる。θ=arctanZim ́/Zre ́F1次に、S70にて算出した特定インピーダンスZ ́の位相角θに基づいて、セル10の空気出口部の酸素濃度を算出するS80。セル10の空気出口部における酸素濃度の算出処理では、特定インピーダンスZ ́の位相角θとセル10の空気出口部における酸素濃度との相関関係を予め関連付けた制御特性制御マップを記憶部54dに記憶し、当該制御特性を参照することで、算出することができる。以上説明した本実施形態では、所定の測定周波数fmにおける内部インピーダンスZのうち、セル10内の酸素濃度の変化に応じて変化する特定インピーダンスZ ́の位相角θを利用して、セル10の空気出口部における酸素濃度の変化を測定する構成としている。このため、従来の如く、長時間を要する燃料電池1の等価回路における回路定数の算出処理を行うことなく、セル10の空気出口部における酸素濃度を短時間で算出可能となる。この結果、本実施形態のガス濃度測定装置をセル10の空気出口部における酸素濃度が短時間で変動する車両用の燃料電池システムに適用したとしても、セル10の空気出口部における酸素濃度を適切に測定することができる。さらに、本実施形態では、測定周波数fmおよび高周波数fhを合成した交流信号を印加する構成であるため、測定周波数fmおよび高周波数fhの交流信号を別個に印加する構成に比べて、セル10の空気出口部における酸素濃度を短時間で算出可能となる。なお、酸素濃度の測定では、酸素濃度8程度までの低い濃度側にて測定精度が出るので、この測定精度を出すにはセル10の全体の空気出口部のトータル面積の10程度の領域で酸素濃度を測定することが望ましい。第2実施形態次に、本発明の第2実施形態について図9に基づいて説明する。ここで、図9は、燃料電池1の負荷が変動した場合のCole-Coleプロット図である。なお、図9では、燃料電池1の負荷が変動した場合に、測定周波数fmに対応する内部インピーダンスZの変化を黒丸で示し、後述する第1の補正用周波数fm1に対応する内部インピーダンスZの変化を黒四角で示し、後述する第2の補正用周波数fm2に対応する内部インピーダンスZの変化を黒三角で示している。燃料電池システムでは、電気負荷2からの要求電力量に応じて、燃料電池1の負荷発電量が変動するが、この燃料電池1の負荷が変動によって、特定インピーダンスZ ́の位相角θが変動してしまうことがある。このため、燃料電池1の負荷変動が、ガス濃度測定装置50の酸素濃度の測定精度に影響する可能性がある。例えば、測定周波数fmの交流信号を燃料電池1に印加したときに、図9に示すように、燃料電池1の負荷が通常よりも高くなると、燃料電池1の内部インピーダンスZが増大して、特定インピーダンスZ ́の位相角がθ2からθ3へと大きくなる。一方、燃料電池1の負荷が通常よりも低くなると、燃料電池1の内部インピーダンスZが減少して、特定インピーダンスZ ́の位相角がθ2からθ1へと小さくなる。そこで、本実施形態では、燃料電池1の負荷変動に応じて、燃料電池1の内部インピーダンスZを測定する際の周波数を変更することで、燃料電池1の負荷変動による酸素濃度の測定精度への影響を抑制する構成としている。具体的には、本実施形態では、まず、測定周波数fmおよび高周波数fhに加えて、測定周波数fmと異なる複数の補正用周波数fm1、fm2を合成した交流信号を燃料電池1に印加する。この補正用周波数fm1、fm2としては、例えば、測定周波数例えば20Hzよりも高い周波数例えば30Hzを第1の補正用周波数fm1とし、測定周波数例えば20Hzよりも低い周波数例えば5Hzを第2の補正用周波数fm2とすることができるfm2<fm<fm1。次に、燃料電池1の負荷変動を信号処理回路53における負荷検出手段53aにて検出する。なお、燃料電池1の負荷変動は、電流測定装置51にて検出する出力電流に含まれる直流成分に限らず、例えば、燃料電池1への要求電力量にも相関関係があるので、当該要求電力量から間接的に燃料電池1の負荷変動を検出してもよい。そして、燃料電池1の負荷が予め設定された第1負荷基準値よりも大きい場合には、燃料電池1の負荷の増大によって特定インピーダンスZ ́の位相角θが大きくなっていると判断できるので、測定周波数fmよりも低い第2の補正用周波数fm2に対応するインピーダンスZを燃料電池1の内部インピーダンスZとして算出する。この場合、第2の補正用周波数fm2における内部インピーダンスZ図9の負荷が大きい場合の軌跡における黒三角参照から算出した特定インピーダンスZ ́の位相角θは、測定周波数fmにおける内部インピーダンスZ図9の負荷が大きい場合の軌跡における黒丸参照から算出した特定インピーダンスZ ́の位相角θよりも小さくなる。このため、測定周波数fmよりも低い第2の補正用周波数fm2で燃料電池1の内部インピーダンスZを算出することで、燃料電池1の負荷変動によって増加した特定インピーダンスZ ́の位相角θを小さくする補正を行うことができる。一方、燃料電池1の負荷が予め第1負荷基準値よりも低く値に設定された第2負荷基準値よりも小さい場合には、燃料電池1の負荷の減少によって特定インピーダンスZ ́の位相角θが小さくなっていると判断できるので、測定周波数fmよりも高い第1の補正用周波数fm1に対応するインピーダンスZを燃料電池1の内部インピーダンスZとして算出する。この場合、第1の補正用周波数fm1における内部インピーダンスZ図9の負荷が小さい場合の軌跡における黒四角参照から算出した特定インピーダンスZ ́の位相角θは、測定周波数fmにおける内部インピーダンスZ図9の負荷が小さい場合の軌跡における黒丸参照から算出した特定インピーダンスZ ́の位相角θよりも大きくなる。このため、測定周波数fmよりも高い第1の補正用周波数fm1で燃料電池1の内部インピーダンスZを算出することで、燃料電池1の負荷変動によって減少した特定インピーダンスZ ́の位相角θを大きくする補正を行うことができる。このように本実施形態では、燃料電池1の負荷変動に応じて、内部インピーダンスZを算出する際の周波数を変更する構成としているので、燃料電池1の負荷変動の影響で増減する特定インピーダンスZ ́の位相角θを補正することができる。この結果、酸素濃度をより適切に測定することが可能となる。なお、本実施形態では、測定周波数fm、高周波数fh、および2つの補正用周波数fm1、fm2を合成した交流信号を燃料電池1に印加しているが、測定周波数fm、高周波数fh、および3つ以上の補正用周波数を合成した交流信号を燃料電池1に印加するようにしてもよい。この場合、燃料電池1の負荷レベルを複数の段階に分けて、負荷が高いレベルになるほど、低い周波数の補正用周波数にて内部インピーダンスZを算出し、負荷が低いレベルになるほど、高い周波数の補正用周波数にて内部インピーダンスZを算出するようにすればよい。第3実施形態次に、本発明の第3実施形態について図10に基づいて説明する。ここで、図10は、燃料電池1の温度が低下した場合の特定インピーダンスZ ́の位相θの変化を示す特性図である。なお、図10では、燃料電池1の温度が低温状態から適温状態へと変化した場合に、測定周波数fmに対応する内部インピーダンスZの変化を黒丸で示し、後述する第3の補正用周波数fm3に対応する内部インピーダンスZの変化を黒四角で示している。燃料電池システムでは、燃料電池1の発電開始の初期段階では、通常、燃料電池1の温度が低い低温状態となるが、この燃料電池1の温度が低温状態では、高温状態に比べて、特定インピーダンスZ ́が増大してしまう。このため、燃料電池1の低温状態において、ガス濃度測定装置50にて酸素濃度を測定すると、燃料電池1の温度の影響によって、測定精度が低下する可能性がある。例えば、測定周波数fmの交流信号を燃料電池1に印加したときに、図10に示すように、燃料電池1の温度が適温よりも低い場合には、燃料電池1の内部インピーダンスZが増大して、特定インピーダンスZ ́の位相角がθ3からθ4へと大きくなる。そこで、本実施形態では、燃料電池1の温度に応じて、燃料電池1の内部インピーダンスZを測定する際の周波数を変更することで、燃料電池1の温度変化による酸素濃度の測定精度への影響を抑制する構成としている。具体的には、本実施形態では、測定周波数fmおよび高周波数fhに加えて、測定周波数fmと異なる補正用周波数fm3を合成した交流信号を燃料電池1に印加する。この補正用周波数fm3としては、例えば、測定周波数例えば20Hzよりも低い周波数例えば10Hzを第3の補正用周波数fm3とすることができるfm3<fm。次に、燃料電池1の温度を温度センサ46にて検出し、検出した燃料電池1の温度が予め設定された基準温度Toより低い場合には、特定インピーダンスZ ́の位相角θが大きくなっていると判断できるので、測定周波数fmよりも低い第3の補正用周波数fm3に対応するインピーダンスZを燃料電池1の内部インピーダンスZとして算出する。この場合、第3の補正用周波数fm3における内部インピーダンスZ図10の温度が低い場合の軌跡における黒四角参照から算出した特定インピーダンスZ ́の位相角θは、測定周波数fmにおける内部インピーダンスZ図10の温度が低い場合の軌跡における黒丸参照から算出した特定インピーダンスZ ́の位相角θよりも小さくなる。つまり、測定周波数fmよりも低い第3の補正用周波数fm3で燃料電池1の内部インピーダンスZを算出することで、燃料電池1の温度によって増加した特定インピーダンスZ ́の位相角θを小さくする補正を行うことができる。このように本実施形態では、燃料電池1の温度に応じて、内部インピーダンスZを算出する際の周波数を変更する構成としているので、燃料電池1の温度によって増減する特定インピーダンスZ ́の位相角θを補正することができる。この結果、酸素濃度をより適切に測定することが可能となる。なお、本実施形態では、測定周波数fm、高周波数fh、および1つの補正用周波数fm3を合成した交流信号を燃料電池1に印加しているが、測定周波数fm、高周波数fh、および2つ以上の補正用周波数を合成した交流信号を燃料電池1に印加するようにしてもよい。この場合、燃料電池1の低温域を複数の段階に分けて、温度が小さい低温域になるほど、低い周波数の補正用周波数にて内部インピーダンスZを算出するようにすればよい。第4実施形態次に、本発明の第4実施形態について図11~図13に基づいて説明する。上述の第2実施形態では、燃料電池1の負荷変動によって特定インピーダンスZ ́の位相角θが変動することを説明したが、本発明者らが検証したところ、この特定インピーダンスZ ́の位相角θの変動は、燃料電池1の負荷変動に伴う内部インピーダンスZに含まれる反応抵抗R2の変化が要因であることがわかっている。本実施形態では、信号処理回路53および演算装置54にて、内部インピーダンスZに含まれる反応抵抗R2を算出し、内部インピーダンスZから反応抵抗R2の成分を除去する補正を行うことで、燃料電池1の負荷変動による酸素濃度の測定精度への影響を抑制する構成としている。信号処理回路53および演算装置54にて行う内部インピーダンスZに含まれる反応抵抗R2を算出手法について説明する。まず、DC-DCコンバータ3から燃料電池1に測定周波数fmおよび高周波数fhを合成した交流信号を印加した後、電流測定装置51にてセル10の出力電流を検出する。その後、信号処理回路53にて電流測定装置51で検出したセル10の出力電流から直流成分を抽出する。そして、信号処理回路53にて抽出した直流成分、および予め当該直流成分と出力電圧との関係を関連付けた電流電圧特性に基づいて、反応抵抗R2を算出する。ここで、図11は、燃料電池1における出力電流電流密度と出力電圧との関係を示す電流電圧特性図IV特性図である。図11に示すように、燃料電池1における出力電流電流密度と出力電圧との関係は、電流密度が基準電流密度以下となる領域以降、ターフェル領域と称する。において、電流密度の増加に伴って出力電圧が指数関数的に減少する傾向がある。また、電流密度が基準電流密度より大きくなる領域以降、線形領域と称する。において、電流密度の増加に伴って線形的に減少する傾向がある。このような傾向を示す電流電圧特性は、ターフェル領域を周知のバトラーボルマーの式を基礎とする数式第1の数式により定義することができ、線形領域を予め出力電流電流密度と出力電圧との比例関係を関連付けた数式第2の数式により定義することができる。なお、本実施形態では、上述の第1の数式および第2の数式により電流電圧特性が定義される。このように、燃料電池1の電流電圧特性は、基準電流密度を境に異なる傾向を示す。このため、反応抵抗R2を算出する際には、基準電流密度を閾値として設定して、信号処理回路53にて抽出した直流成分が当該閾値以下であるか否かを判定する。この結果、直流成分が閾値以下と判定された場合には、直流成分の増加に伴って出力電圧が指数関数的に減少する電流電圧特性に相当する第1の数式バトラーボルマーを基礎とする数式に基づいて、セル10の交換電流密度ioを算出し、この交換電流密度ioから反応抵抗R2を算出する。なお、交換電流密度ioは、反応抵抗R2と逆数の関係にあるため、交換電流密度ioを算出することで反応抵抗R2を算出することができる。一方、信号処理回路53にて抽出した直流成分が当該閾値より大きい場合には、直流成分の増加に伴って出力電圧が直線的に減少する電流電圧特性に相当する第2の数式線形式に基づいて、セル10の交換電流密度ioを算出し、この交換電流密度ioから反応抵抗R2を算出する。なお、信号処理回路53および演算装置54における内部インピーダンスZに含まれる反応抵抗R2を算出するためのソフトウェアおよびハードウェアが反応抵抗算出手段を構成する。次に、信号処理回路53および演算装置54にて行う内部インピーダンスZから反応抵抗R2の影響を除去する補正について図12、図13に基づいて説明する。図12は、内部インピーダンスZから反応抵抗R2の影響を除去する際の出力電流の変化を説明する説明図であり、図12のaが反応抵抗R2の影響を除去する前の出力電流を示し、反応抵抗R2の影響を除去した後の出力電流を示している。まず、信号処理回路53にて電流測定装置51で検出したセル10の出力電流から測定周波数fmに対応する交流成分を抽出する。そして、信号処理回路53にて抽出した交流成分、および内部インピーダンスZに含まれる反応抵抗R2により、反応抵抗R2の影響を交流電流に置き換えた関数Htを生成する。その後、信号処理回路53にて電流測定装置51で検出したセル10の出力電流のデータI=HtFtに反応抵抗R2の影響を示す関数Htを畳み込み積分することで、図12に示す反応抵抗R2の影響を除去した出力電流I=Ftを得る。このようにして得られた出力電流I=Ftおよび燃料電池1の出力電圧を信号処理回路53にて高速フーリエ変換処理FFTを行うことで、測定周波数における内部インピーダンスZから反応抵抗R2の影響を除去したインピーダンスを得る。そして、内部インピーダンスZから反応抵抗R2の影響を除去したインピーダンスから高周波インピーダンスZhを減算することで、セル10内の酸素の濃度変化に応じて変化する特定インピーダンスZ ́を算出する。ここで、図13は、燃料電池1の負荷が変動した場合のCole-Coleプロット図であり、図13のaが内部インピーダンスZから反応抵抗R2の影響を除去する前のCole-Coleプロット図を示し、図13のbが内部インピーダンスZから反応抵抗R2の影響を除去した後のCole-Coleプロット図を示している。図13aに示すように、内部インピーダンスZから反応抵抗R2の影響を除去する前は、燃料電池1の負荷に応じて、特定インピーダンスZ ́の位相角θが変動する。これに対して、図13bに示すように、内部インピーダンスZから反応抵抗R2の影響を除去した後は、燃料電池1の負荷によらず、特定インピーダンスZ ́の位相角θが略一定となる。以上説明した本実施形態によれば、内部インピーダンスZから燃料電池1の負荷変動によって変動する反応抵抗R2を除去するので、燃料電池1の負荷変動によって特定インピーダンスZ ́の位相角θが変動してしまうことを抑制することができる。この結果、酸化剤ガスのガス濃度をより適切に測定することが可能となる。ところで、燃料電池1は、発電を繰り返すことで電極101b、101c内の触媒が経時劣化することで、徐々に内部インピーダンスZが増大することがあるが、当該触媒の経時劣化は、反応抵抗R2に相関関係を有することがわかっている。そこで、本実施形態では、反応抵抗R2を算出する度に当該反応抵抗R2の抵抗値を記憶し、記憶した反応抵抗R2の変化出力電圧の変化によって燃料電池1の電極101b、101cにおける触媒の劣化度合いを判定する。これにより、当該触媒の劣化度合いの判定に応じて、内部インピーダンスZに当該触媒の劣化による酸素濃度の測定精度への影響を抑制することが可能となる。他の実施形態以上、本発明の実施例について説明したが、本発明はこれに限定されるものではなく、各請求項に記載した範囲を逸脱しない限り、各請求項の記載文言に限定されず、当業者がそれらから容易に置き換えられる範囲にも及び、かつ、当業者が通常有する知識に基づく改良を適宜付加することができる。例えば、以下のように種々変形可能である。1上述の各実施形態では、特定インピーダンスZ ́の位相角θから酸素濃度を算出する際に、制御特性制御マップを参照しているが、これに限らず、特定インピーダンスZ ́の位相角θとセル10の空気出口部における酸素濃度との相関関係を規定する数式を導出し、導出した数式を用いて酸素濃度を算出するようにしてもよい。2上述の各実施形態では、DC-DCコンバータ3にて燃料電池1の出力信号に交流信号を印加する交流信号印加手段を構成しているが、これに限らず、DC-DCコンバータ3とは別体の発振回路等にて交流信号印加手段を構成してもよい。3上述の各実施形態では、測定周波数fmとして、予め基準ガス濃度とした条件で、印加する周波数を高周波から低周波まで変化させた際に、Cole-Coleプロット図における半円状の軌跡の頂部に達したときの周波数を設定する構成としているが、これに限定されない。測定周波数fmとしては、印加した際に特定インピーダンスZ ́の位相角θが酸素濃度の変化によって変動する周波数であれば任意の周波数に設定してもよい。4上述の各実施形態では、燃料電池システムの制御装置60とガス濃度測定装置50の演算装置54を別体で構成したが、制御装置60に演算装置54としての機能を追加することで、制御装置60と演算装置54とを一体化してもよい。5上述の第2実施形態で説明した燃料電池1の負荷変動に応じて特定インピーダンスZ ́の位相角θの補正と、第3実施形態で説明した燃料電池1の温度に応じて特定インピーダンスZ ́の位相角θの補正とを組み合わせてもよい。例えば、燃料電池1の負荷変動が大きい場合に、第2実施形態で説明した特定インピーダンスZ ́の位相角θの補正を行い、燃料電池1の負荷変動が小さい場合に、第3実施形態で説明した特定インピーダンスZ ́の位相角θの補正を行うといったように、優先順位をつけることで、組み合わせることができる。6上述の第2、第3実施形態では、燃料電池1の負荷変動および燃料電池1の温度に応じて特定インピーダンスZ ́の位相角θを補正する構成としている。これらとは別に、燃料電池1は、発電を繰り返すことで電極101b、101c内の触媒が経時的に劣化して、徐々に特定インピーダンスZ ́が増大してしまうことがある。このため、燃料電池1の電極101b、101c内の触媒が経時的に劣化の影響でガス濃度測定装置50の測定精度が低くなる可能性がある。そこで、セル10内の酸素濃度、燃料電池1の負荷、温度等を予め設定した基準値とした条件で、燃料電池1の電極101b、101c内の触媒の経時劣化による特定インピーダンスZ ́の位相角θと、触媒が経時劣化する前の特定インピーダンスZ ́の位相角θとの変化量Δθを算出し、当該変化量Δθに基づいて、燃料電池1の電極101b、101c内の触媒の経時劣化による酸素濃度の測定精度への影響を抑制する構成としてもよい。7上記各実施形態では、本発明の燃料電池システムを電気自動車に適用した例を説明したが、これに限定されず、船舶及びポータブル発電器等の移動体に適用してもよい。