燃料电池的状态检测装置以及方法与流程

本发明涉及一种燃料电池的状态检测装置以及方法。

背景技术：

已知一种测定燃料电池的电压值和阻抗值、基于这些值来检测燃料电池的内部状态的燃料电池的状态检测装置。

例如，日本专利第4640661号提出了以下方案：分别计算与电解质膜电阻对应的第一频域下的第一阻抗以及与电解质膜电阻同催化剂层电阻的合计值对应的低于第一频域的第二频域下的第二阻抗，基于作为第二阻抗与第一阻抗之差的差分阻抗来计算催化剂层的含水量。

另外，日本特开2005-285614号公报中记载了以下方案：获取与燃料电池的复阻抗曲线(cole-cole图)的同实轴的交点处的频率f1、表示氧发生反应时的反应电阻(阴极的反应电阻)的第一区域内的频率f2以及表示与氧的扩散有关的电阻的第二区域内的频率f3对应的复阻抗，根据获取到的复阻抗来求出内部电阻值。

技术实现要素：

然而，在日本专利第4640661号中，无法分别掌握阳极(anode)的状态量和阴极(cathode)的状态量。另外，在日本特开2005-285614号公报中也是，在阻抗曲线中阳极的状态和阴极的状态混在一起，难以个别地掌握阳极的状态量和阴极的状态量。

本发明是着眼于这种问题而完成的，其目的在于提供一种能够个别地检测燃料电池中的阳极的状态量、阴极的状态量等内部状态量的燃料电池的状态检测装置以及方法。

用于解决问题的方案

根据本发明的某个方式，提供一种接受阳极气体和阴极气体的供给来进行发电的燃料电池的状态检测装置。更详细地说，该状态检测装置具备阻抗获取单元，该阻抗获取单元获取高频阻抗和低频阻抗，该高频阻抗是基于从至少包括对阳极的状态量表现出响应性的频带的高频带选择出的频率的阻抗，该低频阻抗是基于从至少包括对阴极的状态量表现出响应性的频带的低频带选择出的频率的阻抗。另外，状态检测装置具备内部状态量估计单元，该内部状态量估计单元将获取到的高频阻抗与低频阻抗进行组合，分别估计作为燃料电池的内部状态的阳极的状态量和阴极的状态量。

附图说明

图1是本发明的实施方式的燃料电池单元的立体图。

图2是图1的燃料电池的ii-ii截面图。

图3是本发明的实施方式的燃料电池系统的概要结构图。

图4a是表示在输入了低频带的交流电压的情况下流过燃料电池的简易等效电路模型的电流的路径的图。

图4b是表示在输入了比图4a的情况高的频带的交流电压的情况下流过燃料电池的简易等效电路模型的电流的路径的图。

图4c是表示在输入了比图4b的情况高的频带的交流电压的情况下流过燃料电池简易等效电路模型的电流的路径的图。

图4d是表示在输入了高频带的交流电压的情况下流过燃料电池的简易等效电路模型的电流的路径的图。

图5是表示一个实施方式所涉及的状态量的估计的流程的流程图。

图6是表示一个实施方式所涉及的状态量的估计的流程的流程图。

图7是分别示出稳定时和非稳定时的燃料电池的i-v特性线的图。

图8是表示一个实施方式所涉及的状态量的估计的流程的流程图。

图9示出了阴极的双电层电容的候选的频率响应。

图10a示出了阳极的双电层电容的候选的频率响应。

图10b示出了阳极112的反应电阻值的候选的频率响应。

图11是表示一个实施方式所涉及的状态量的估计的流程的流程图。

图12示出了稳定时的燃料电池1的i-v特性线。

图13是说明用于进行i-v特性线的斜率δv/δi的计算的一组电流和电压的设定方法的一例的图。

图14是概要性地示出在一个实施方式所涉及的燃料电池系统中阻抗测量所涉及的重要部分的框图。

具体实施方式

下面，参照附图等来说明本发明的实施方式。

燃料电池单元构成为用作为燃料极的阳极和作为氧化剂极的阴极将电解质膜夹在中间。在燃料电池单元中，向阳极供给含有氢的阳极气体，另一方面向阴极供给含有氧的阴极气体，通过使用这些气体来进行发电。阳极和阴极这两个电极处进行的电极反应如下。

阳极：2h2→4h⁺+4e^-

阴极：4h⁺+4e^-+o2→2h2o

图1和图2是用于说明本发明的一个实施方式的燃料电池单元10的结构的图。图1是燃料电池单元10的立体图，图2是图1的燃料电池单元10的ii-ii截面图。

如图1和图2所示，燃料电池单元10具备膜电极组件(mea)11以及以将mea11夹在中间的方式配置的阳极隔板12和阴极隔板13。

mea11由电解质膜111、阳极112以及阴极113构成。mea11在电解质膜111的其中一面侧具有阳极112，在另一面侧具有阴极113。

电解质膜111是由氟系树脂形成的质子传导性的离子交换膜。电解质膜111在湿润状态下表现出良好的电传导性。此外，作为电解质膜111，也可以根据设想的燃料电池的对应，例如使用使磷酸(h3po4)浸渍于规定的基质所得到的材料等其它材料。

阳极112具备催化剂层112a和气体扩散层112b。催化剂层112a是由铂或承载有铂等的炭黑粒子形成的构件，设置成与电解质膜111接触。气体扩散层112b配置于催化剂层112a的外侧。气体扩散层112b是由具有气体扩散性和导电性的碳布形成的构件，设置成与催化剂层112a及阳极隔板12接触。

与阳极112同样地，阴极113也具备催化剂层113a和气体扩散层113b。催化剂层113a配置于电解质膜111与气体扩散层113b之间，气体扩散层113b配置于催化剂层113a与阴极隔板13之间。

阳极隔板12配置于气体扩散层112b的外侧。阳极隔板12具备用于向阳极112供给阳极气体(氢气)的多个阳极气体流路121。阳极气体流路121形成为槽状通路。

阴极隔板13配置于气体扩散层113b的外侧。阴极隔板13具备用于向阴极113供给阴极气体(空气)的多个阴极气体流路131。阴极气体流路131形成为槽状通路。

阳极隔板12和阴极隔板13构成为使在阳极气体流路121中流动的阳极气体的流动方向与在阴极气体流路131中流动的阴极气体的流动方向互为反向。此外，阳极隔板12和阴极隔板13也可以构成为使这些气体的流动方向为向相同方向流动。

在将这种燃料电池单元10用作汽车用电源的情况下，由于所要求的电力大，因此作为将数百块燃料电池单元10层叠而得到的燃料电池堆来使用。然后，构成向燃料电池堆供给阳极气体和阴极气体的燃料电池系统，取出用于驱动车辆的电力。此外，在本实施方式中，以层叠燃料电池单元10所得到的燃料电池堆为单位来进行后述的阻抗测定，但是也可以是以一块燃料电池单元10为单位或以燃料电池堆的一部分(例如数十块单元)为单位来进行阻抗测定。

另外，在燃料电池堆中，通过将一块燃料电池单元10中的阳极112、阴极113以及电解质膜111串联配置多块来构成为作为总和的阳极、阴极以及电解质膜。然而，下面为了便于说明，对该作为总和的阳极、阴极以及电解质膜也标注与单元单体的阳极112、阴极113以及电解质膜111相同的标记。

图3是本发明的一个实施方式的燃料电池系统100的概要图。

燃料电池系统100具备燃料电池1、阴极气体供排装置2、阳极气体供排装置3、电力系统5以及控制器6。

燃料电池1是如上所述那样层叠多块燃料电池单元10(单位电池)而成的层叠电池。燃料电池1接受阳极气体和阴极气体的供给来发出车辆行驶所需的电力。燃料电池1具有阳极侧端子1a和阴极侧端子1b作为取出电力的输出端子。

阴极气体供排装置2向燃料电池1供给阴极气体，并且将从燃料电池1排出的阴极排气排出到外部。阴极气体供排装置2具备阴极气体供给通路21、阴极气体排出通路22、过滤器23、气流传感器24、阴极压缩机25、阴极压力传感器26、水分回收装置(wrd；waterrecoverydevice)27以及阴极压力调节阀28。

阴极气体供给通路21是流通向燃料电池1供给的阴极气体的通路。阴极气体供给通路21的一端连接于过滤器23，另一端连接于燃料电池1的阴极气体入口部。

阴极气体排出通路22是流通从燃料电池1排出的阴极排气的通路。阴极气体排出通路22的一端连接于燃料电池1的阴极气体出口部，另一端形成为开口端。阴极排气是包含阴极气体、通过电极反应而产生的水蒸气等的混合气体。

过滤器23是将取入到阴极气体供给通路21的阴极气体中含有的尘、埃等去除的构件。

阴极压缩机25设置于比过滤器23更靠下游侧的阴极气体供给通路21。阴极压缩机25加压输送阴极气体供给通路21内的阴极气体来供给到燃料电池1。

气流传感器24设置于过滤器23与阴极压缩机25之间的阴极气体供给通路21。气流传感器24检测供给到燃料电池1的阴极气体的流量。

阴极压力传感器26设置于阴极压缩机25与wrd27之间的阴极气体供给通路21。阴极压力传感器26检测供给到燃料电池1的阴极气体的压力。由阴极压力传感器26检测出的阴极气体压力代表包括燃料电池1的阴极气体流路等在内的整个阴极系统的压力。

wrd27横跨阴极气体供给通路21和阴极气体排出通路22地将它们连接。wrd27是如下的装置：回收在阴极气体排出通路22中流动的阴极排气中的水分，使用所回收的该水分来加湿在阴极气体供给通路21中流动的阴极气体。

阴极压力调节阀28设置于比wrd27更靠下游的阴极气体排出通路22。阴极压力调节阀28由控制器6来控制开闭，对供给到燃料电池1的阴极气体的压力进行调整。

接着，说明阳极气体供排装置3。

阳极气体供排装置3向燃料电池1供给阳极气体，并且将从燃料电池1排出的阳极排气排出到阴极气体排出通路22。阳极气体供排装置3具备高压罐31、阳极气体供给通路32、阳极压力调节阀33、阳极压力传感器34、阳极气体排出通路35、缓冲罐36、放气通路37以及放气阀38。

高压罐31是将要向燃料电池1供给的阳极气体保持为高压状态来贮存的容器。

阳极气体供给通路32是将从高压罐31排出的阳极气体供给到燃料电池1的通路。阳极气体供给通路32的一端连接于高压罐31，另一端连接于燃料电池1的阳极气体入口部。

阳极压力调节阀33设置于比高压罐31更靠下游的阳极气体供给通路32。阳极压力调节阀33由控制器6来控制开闭，对供给到燃料电池1的阳极气体的压力进行调整。

阳极压力传感器34设置于比阳极压力调节阀33更靠下游的阳极气体供给通路32。阳极压力传感器34检测供给到燃料电池1的阳极气体的压力。由阳极压力传感器34检测出的阳极气体压力代表包括缓冲罐36、燃料电池1的阳极气体流路等在内的整个阳极系统的压力。

阳极气体排出通路35是流通从燃料电池1排出的阳极排气的通路。阳极气体排出通路35的一端连接于燃料电池1的阳极气体出口部，另一端连接于缓冲罐36。阳极排气中包含电极反应中未被使用的阳极气体、从阴极气体流路131向阳极气体流路121泄漏过来的氮等杂质气体、水分等。

缓冲罐36是暂时蓄积通过阳极气体排出通路35流过来的阳极排气的容器。积存在缓冲罐36中的阳极排气在放气阀38被打开时通过放气通路37排出到阴极气体排出通路22。

放气通路37是用于排出阳极排气的通路。放气通路37的一端连接于阳极气体排出通路35，另一端连接于比阴极压力调节阀28更靠下游的阴极气体排出通路22。

放气阀38设置于放气通路37。放气阀38由控制器6来控制开闭，对从阳极气体排出通路35排出到阴极气体排出通路22的阳极排气的放气流量进行控制。

当执行放气阀38为开阀状态的放气控制时，阳极排气通过放气通路37和阴极气体排出通路22排出到外部。此时，阳极排气在阴极气体排出通路22内与阴极排气混合。通过像这样使阳极排气与阴极排气混合后排出到外部，混合气体中的阳极气体浓度(氢浓度)被设定为排出容许浓度以下的值。

电力系统5具备电流传感器51、电压传感器52、行驶马达53、逆变器54、蓄电池55以及dc/dc转换器56。

电流传感器51检测从燃料电池1取出的输出电流。电压传感器52检测燃料电池1的输出电压、也就是说阳极侧端子1a与阴极侧端子1b之间的端子间电压。电压传感器52既可以构成为检测每块燃料电池单元10的电压，也可以构成为检测每多块燃料电池单元10的电压。

行驶马达53是三相交流同步马达，是用于驱动车轮的驱动源。行驶马达53具有作为电动机的功能和作为发电机的功能，该作为电动机的功能是从燃料电池1和蓄电池55接受电力的供给来进行旋转驱动，该作为发电机的功能是通过被外力驱动旋转来进行发电。

逆变器54由多个igbt等半导体开关构成。逆变器54的半导体开关由控制器6控制开关，由此将直流电力变换为交流电力，或者将交流电力变换为直流电力。在使行驶马达53作为电动机而发挥功能的情况下，逆变器54将燃料电池1的输出电力与蓄电池55的输出电力的合成直流电力变换为三相交流电力来供给到行驶马达53。与此相对，在使行驶马达53作为发电机而发挥功能的情况下，逆变器54将行驶马达53的再生电力(三相交流电力)变换为直流电力来供给到蓄电池55。

蓄电池55构成为被充入燃料电池1的输出电力的剩余部分和行驶马达53的再生电力。充入到蓄电池55的电力根据需要而被供给到阴极压缩机25等辅机类、行驶马达53。

dc/dc转换器56是使燃料电池1的输出电压升降的双向性的电压变换机。通过利用dc/dc转换器56对燃料电池1的输出电压进行控制，来调整燃料电池1的输出电流等。

控制器6由具备中央运算装置(cpu)、只读存储器(rom)、随机存取存储器(ram)以及输入输出接口(i/o接口)的微计算机构成。除了来自电流传感器51、电压传感器52等各种传感器的信号以外，来自检测加速踏板的踏下量的加速行程传感器(未图示)等传感器的信号也被输入到控制器6。

控制器6根据燃料电池系统100的运转状态来控制阳极压力调节阀33、阴极压力调节阀28、阴极压缩机25等，对供给到燃料电池1的阳极气体、阴极气体的压力、流量进行调整。

另外，控制器6基于行驶马达53的要求电力、阴极压缩机25等辅机类的要求电力、蓄电池55的充放电要求等来计算目标输出电力。控制器6基于目标输出电力，参照预先决定的燃料电池1的iv特性(电流电压特性)来计算燃料电池1的目标输出电流。然后，控制器6进行以下控制：通过dc/dc转换器56对燃料电池1的输出电压进行控制以使燃料电池1的输出电流为目标输出电流，供给行驶马达53、辅机类所需的电流。

另外，控制器6对阴极压缩机25等进行控制以使燃料电池1的各电解质膜111的湿润度(含水量)为适于发电的状态。

另外，在后述的第一实施方式～第六实施方式中，控制器6将在燃料电池1的输出电压上叠加规定频率的交流信号所得到的电压值的振幅值除以同样叠加交流信号所得到的电流值的振幅值，来计算规定频率下的燃料电池1的阻抗z。

在如上所述那样说明的燃料电池系统100中，由控制器6、电流传感器51、电压传感器52以及dc/dc转换器56来构成燃料电池1的状态检测装置。

在本实施方式中，在燃料电池1中，设定考虑到作为阳极112的状态量的反应电阻值ra和双电层电容ca、作为阴极113的状态量的反应电阻值rc和双电层电容cc以及作为电解质膜111的状态量的电解质膜电阻值rm的简易的等效电路模型，基于该简易等效电路模型来进行燃料电池1的状态估计。

此外，电解质膜电阻值rm是其值根据电解质膜111的湿润度决定的状态量。通常，具有电解质膜111越干燥则电解质膜电阻值rm越高的趋势。

另外，阳极112的反应电阻值ra与阳极112处的阳极气体的反应相应地增减，例如，当存在阳极气体不足等无法顺畅地进行该反应的因素时，与之相应地反应电阻值ra上升。

并且，阳极112的双电层电容ca是以表示燃料电池1中阳极112所具有的电容的方式进行模型化而得到的。因而，双电层电容ca是基于构成阳极112的材料、大小等各种要素而决定的。

另外，阴极113的反应电阻值rc与阴极113处的阴极气体的反应相应地增减，例如，当存在阴极气体不足等无法顺畅地进行该反应的因素时，与之相应地反应电阻值rc上升。

并且，阴极113的双电层电容cc是以表示阴极113所具有的电容的方式进行模型化而得到的。因而，双电层电容成分的值cc是基于构成阴极113的材料、大小等各种要素而决定的。

在此，本发明人们发现：在燃料电池1的简易等效电路模型中，叠加于燃料电池1的输出电流的交流信号(交流电流)在燃料电池内部流动的路径存在频率依赖特性。下面，说明交流电流在燃料电池内部流动的路径的频率依赖特性。

图4a～图4d是按交流电流的频带来示意性地示出在本实施方式所涉及的燃料电池1的简易等效电路模型中叠加于燃料电池1的输出电流的交流电流在燃料电池内部流动的路径的图。

在图4a中，例如示出了属于0hz附近的低频带(下面，也记载为第一频带)的频率的交流电流的路径。另外，在图4b中，例如示出了属于与几hz左右的第一频带相比稍高的频带(下面，也记载为第二频带)的频率的交流电流的路径。并且，在图4c中，例如示出了属于与几十hz～几千hz的第二频带相比稍高的频带(下面，也记载为第三频带)的频率的交流电流的路径。另外，在图4d中，例如示出了属于几万hz以上的最高频带(下面，也记载为第四频带)的频率的交流电流的路径。此外，在图4a～图4d中，以粗线示出交流电流的路径。

首先，关于图4a所示的属于第一频带的频率的交流电流，由于频率低，因此值的变动平缓，电流的值具有与作为固定值的直流相近的性质。因而，这样与直流相近的性质的交流电流不流入阳极112的双电层电容和阴极113的双电层电容侧部分，或者即使流入阳极112的双电层电容和阴极113的双电层电容侧部分，其大小也小至能够忽视的程度。即，如图所示，交流电流实质上仅流过阳极112的反应电阻、电解质膜电阻以及阴极113的反应电阻的部分。

接着，关于图4b所示的属于第二频带的频率的交流电流，与上述属于第一频带的频率的交流电流相比，值的变动变大，作为交流的性质变得更强。因而，认为如图所示那样，交流电流也开始流过阴极113的双电层电容侧。

另一方面，已知阳极112的反应电阻值ra与阴极113的反应电阻值rc相比取非常小的值，因此电流比较容易流向阳极112的反应电阻侧。因而，认为在第二频带的频率程度的交流电流中，电流依然不流入阳极112的双电层电容侧部分，或者即使流入阳极112的双电层电容侧部分，其大小也小至能够忽视的程度。

并且，关于图4c所示的属于第三频带的频率的交流电流，与上述属于第二频带的频率的交流电流相比，值的变动变得更大，因此作为交流的性质进一步变强。因而，认为阳极112的双电层电容的影响也变得无法忽视，电流也流过阳极112的双电层电容部分。

另一方面，在该第三频带，阴极113处的氧化还原反应变得无法追随上述交流电流的值的变动速度，产生表观上不再发生该氧化还原反应的状态。

因而，阴极113处的阴极气体的反应实质上不发生，因此能够忽视上述氧化还原反应对阴极113的反应电阻的影响。

即，在第三频带，如图4c所示，交流电流不流入阴极113的反应电阻，或者即使流入阴极113的反应电阻、其大小也小至能够忽视的程度，因此认为实质上仅流过双电层电容成分。

此外，在阳极112处，氧化还原反应对交流电流的值的变动的追随性能比较高，该氧化还原反应能够在第三频带时仍追随交流电流的值的变动。因而，认为如图所示那样，属于第三频带的频率的交流电流依然流过阳极112的反应电阻。

然后，关于图4d所示的属于第四频带的频率的交流电流，与上述属于第三频带的频率的交流电流相比，值的变动进一步变大，因此不仅是阴极113、阳极112处的氧化还原反应也变得无法追随该交流电流的值的变动。

因而，除了阴极113以外，阳极112处的反应也实质上不发生，能够忽视阴极113的反应电阻和阳极112的反应电阻双方的影响。

即，在第四频带，交流电流不流入阴极113和阳极112双方的反应电阻，或者即使流入阴极113和阳极112双方的反应电阻，其大小也小至能够忽视的程度。因而，如图所示，属于第四频带的频率的交流电流仅流过阴极113和阳极112各自的双电层电容侧。

根据以上的说明可以了解的是，对于上述的从第一频带选择的频率的交流电流、从第二频带选择的频率的交流电流、从第三频带选择的频率的交流电流以及从第四频带选择的频率的交流电流而言，流过燃料电池的简易等效电路中的各要素的路径不同。

因而，本发明人们想到了：利用像这样与频率相应的交流电流的路径的不同，参照基于简易等效电路而得到的阻抗的式子，

[式1]

(其中，j表示虚数单位。)

根据基于属于各频带的频率的阻抗来个别地估计各种状态量。

例如，从上述第四频带(下面也记载为“电解质膜响应频带”)选择的频率的交流电流流过电解质膜电阻、阳极112的双电层电容以及阴极113的双电层电容的部分，因此基于从该电解质膜响应频带选择的频率的阻抗(下面也记载为“电解质膜响应阻抗”)包含电解质膜电阻值rm的信息。

此外，该电解质膜响应频带是用于所谓的hfr(highfrequencyresistance：高频电阻)的频带。因而，在阻抗的式(1)中当设ω→∞时，能够视作阻抗z与电解质膜电阻值rm大致一致。

另外，从第三频带(下面也记载为“阳极响应频带”)选择的频率的交流电流流过电解质膜电阻、阳极112的反应电阻、阳极112的双电层电容以及阴极113的双电层电容的部分，因此基于从该阳极响应频带选择的频率的阻抗(下面记载为“阳极响应阻抗”)至少包含阳极112的反应电阻值ra和阳极112的双电层电容值ca的信息。

特别是在该情况下，在图4c所示的等效电路中能够忽视阴极113的反应电阻，因此阻抗的式子呈现为下式。

[式2]

并且，从第二频带选择的频率的交流电流流过电解质膜电阻、阳极112的反应电阻、阴极113的反应电阻以及阴极113的双电层电容的部分，因此基于从该第二频带选择的频率的阻抗包含电解质膜电阻值、阳极112的反应电阻值、阴极113的反应电阻值rc以及阴极113的双电层电容值cc的信息来作为状态量。

另外，从作为最低频带的第一频带(下面也记载为“低频带”)选择的频率的交流电流流过电解质膜电阻、阳极112的反应电阻以及阴极113的反应电阻的部分，因此基于从该低频带选择的频率的阻抗(下面记载为低频阻抗)至少包含阴极113的反应电阻值rc的信息。

下面，在各实施方式中，说明使用上述电解质膜响应频带、阳极响应频带以及低频带燃料中的至少两个来进行的各状态量的估计的详情。

此外，一般已知的是，“频率f”与“角频率ω”之间存在ω＝2πf的关系，它们之间只存在乘以无量纲的常数2π的差异，因此在各实施方式中将“频率”与“角频率”同等看待，无论表示哪个都使用“ω”的符号，以简化说明。

(第一实施方式)

下面，说明第一实施方式。

图5是表示本实施方式所涉及的状态量的估计的流程的流程图。

如图所示，首先，在步骤s101中，选择电解质膜响应频带中的一个点的频率ωh，求出基于频率ωh的阻抗z(ωh)。

具体地说，在阻抗测量定时，控制器6控制dc/dc转换器56，使得在从燃料电池1输出的输出电压和输出电流上叠加电解质膜响应频带的频率ωh的交流信号。

并且，控制器6对由电压传感器52测定出的输出电压的值v实施傅立叶变换来得到电压振幅值v(ωh)，对由电流传感器51测定出的输出电流的值i实施傅立叶变换处理来得到电压振幅值i(ωh)，求出它们之比v(ωh)/i(ωh)来作为阻抗z(ωh)。此外，测量阻抗z(ωh)的手法在针对从电解质膜响应频带以外的阳极响应频带、低频带选择出的频率进行的情况下也是同样的，因此以后省略其详细说明。

接着，在步骤s102中，控制器6根据所得到的阻抗z(ωh)来估计电解质膜电阻值rm。具体地说，如上所述，电解质膜响应频带是在所谓的hfr测量中使用的频带，基于从该高频带选择出的频率ωh的阻抗z(ωh)或其实部成分zr(ωh)大致与电解质膜电阻值rm一致。即，将阻抗z(ωh)或其实部成分zr(ωh)的值直接估计为电解质膜电阻值rm。

在步骤s103中，控制器6选择阳极响应频带中的两个点的频率ω1、ω2，求出基于该频率ω1、ω2的阳极响应阻抗z(ω1)、z(ω2)。

在步骤s104中，控制器6根据估计出的电解质膜电阻值rm以及所得到的两个阻抗z(ω1)、z(ω2)，来估计阳极112的反应电阻值ra和阳极112的双电层电容值ca。

具体说明该估计的方式。首先，在选择阳极响应频带中的两个点的频率ω1、ω2的情况下，如上所述那样能够忽视阴极113的反应电阻，因而，作为阻抗的式子，能够使用从基于简易等效电路的阻抗的式(1)去除阴极113的反应电阻值ra后的式(2)。

在此，在式(2)中，代入作为已知的值的两个点的频率ω1、ω2以及基于它们的阻抗z(ω1)、z(ω2)的组合，取其实部zr(ω1)和zr(ωm2)。而且，考虑到估计出的电解质膜电阻值rm已知，能够得到未知数为ra和ca的两个方程式。因而，只要对得到的两个方程式进行求解就能够求出ra和ca。

示出求出未知数ra和ca的方法的一例。首先，当取式(2)的实部来进行变形时，为下式。

[式3]

考虑横轴为ω²、纵轴为1/zr的平面，在该平面中，式(3)示为直线，其斜率mr呈现为下式。

[式4]

在此，两个点的频率ω1、ω2已知，因此当将这两个点的频率ω1、ω2以及与它们对应的阻抗测量值的实部zr(ω1)、zr(ω2)绘制于上述平面时，将它们连接的直线可定，斜率mr的值可定。即，式(4)的未知数是ra和ca。

接着，式(3)所表示的直线的截距a呈现为下式。

[式5]

在此，与斜率mr的值同样地，根据点的频率ω1、ω2以及与它们对应的阻抗测量值的实部zr1、zr2，截距a的值也可定。而且，zr与阻抗测量值的实部zr1及zr2相当，因此式(5)的未知数只有ra。

因而，根据式(5)，能够将阳极112的反应电阻值ra以下式的方式求出。

[式6]

另外，通过将根据式(6)决定的ra代入到式(4)，能够将阳极112的双电层电容值ca以下式的方式求出。

[式7]

此外，求出ra和ca的计算方法不限于上述的计算方法，能够使用各种适当的计算方法。

接着，在步骤s105中，控制器6选择低频带中的一个点的频率ωl，测量基于该频率ωl的阻抗z(ωl)。

在步骤s106中，控制器6使用已估计出的电解质膜电阻值rm、阳极112的反应电阻值ra和阳极112的双电层电容值ca以及测量出的阻抗z(ωl)，来估计阴极113的双电层电容值cc。

具体说明该估计的方式。如上所述，低频带的频率ωl的交流电流流过燃料电池1的简易等效电路中的全部电路要素、即阳极112的反应电阻和双电层电容、电解质膜电阻以及阴极113的反应电阻和双电层电容的部分。因而，基于频率ωl得到的低频阻抗z(ωl)包含阳极112的反应电阻值ra和双电层电容值ca、电解质膜电阻值rm以及阴极113的反应电阻值rc和双电层电容值cc的信息。因而，作为阻抗的式子，需要使用考虑到上述全部电路要素的式(1)。

在式(1)中，代入作为已知的值的频率ωl以及基于频率ωl的阻抗z(ωl)，取其实部zr(ωl)和虚部zi(ωl)。而且，考虑到估计出的电解质膜电阻值rm、阳极112的反应电阻值ra以及阳极112的双电层电容值ca已知，能够得到未知数为rc和cc的两个方程式。因而，只要对这两个方程式进行求解就能够求出未知数rc和cc。

示出像这样求出未知数rc和cc的方法的一例。首先，当取式(1)的实部来进行变形时，为下式。

[式8]

另外，当取式(1)的虚部来进行变形时，为下式。

[式9]

在此，频率ωl、与频率ωl对应的阻抗测量值的实部zr(ωl)和虚部zi(ωl)以及ra和ca已知，因此当将它们代入到式(8)和式(9)来进行变形时，阴极113的双电层电容值cc为下式。

[式10]

其中，在式(10)中ω是ωl，a能够如下述的式(11)那样定义。

[式11]

并且，阴极113的反应电阻值rc以下式的方式被求出。

[式12]

其中，式(12)中的a如上述式(11)那样定义，式(12)中的b如下述式(13)那样定义。

[式13]

如以上那样，通过步骤s101～步骤s106的步骤，作为燃料电池1的状态量，能够估计出电解质膜电阻值rm、阳极112的反应电阻值ra、阳极112的双电层电容值ca、阴极113的反应电阻值rc以及阴极113的双电层电容值cc。

根据上述的本实施方式，能够得到以下的效果。在本实施方式中，由控制器6、电流传感器51、电压传感器52以及dc/dc转换器56来构成状态检测装置。另外，阻抗获取单元和内部状态量估计单元由控制器6构成。

根据本实施方式，接受阳极气体和阴极气体的供给来进行发电的燃料电池1的状态检测装置的阻抗获取单元获取基于从至少包括对阳极112的状态量ra、ca表现出响应性的频带的高频带(阳极响应频带和电解质膜响应频带)选择出的频率ωh、ω1、ω2的高频阻抗z(ωh)、z(ω1)、z(ω2)以及基于从至少包括对阴极的状态量rc、cc表现出响应性的频带的低频带选择出的频率ωl的低频阻抗z(ωl)(步骤s101、步骤s103、步骤s105)。

然后，燃料电池1的状态检测装置的内部状态量估计单元将获取到的高频阻抗z(ωh)、z(ω1)、z(ω2)与低频阻抗z(ωl)进行组合，分别估计作为燃料电池1的内部状态的阳极112的状态量ra、ca和阴极113的状态量rc、cc。

据此，能够基于获取到的高频阻抗z(ωh)、z(ω1)、z(ω2)和低频阻抗z(ωl)这样的从不同频带得到的阻抗信息，利用与频率的大小相应的阳极112的反应和阴极113的反应对于电流变动的追随速度差，来至少分别个别地探测阳极112的状态量ra、ca和阴极113的状态量rc、cc。因而，能够得到高精度的阳极112的状态量ra、ca和阴极113的状态量(rc、cc)的信息，其结果是能够使利用这些状态量进行的燃料电池1的动作控制更适当。

并且，根据本实施方式，内部状态量估计单元基于高频阻抗z(ωh)、z(ω1)、z(ω2)来估计某个内部状态量rm、ra、ca，基于估计出的该内部状态量rm、ra、ca以及低频阻抗z(ωl)来估计其它内部状态量rc、cc。

由此，对于仅利用作为一个频带的低频带的低频阻抗z(ωl)无法确定的内部状态量rc、cc，能够基于根据作为其它频带的高频带的高频阻抗z(ωh)、z(ω1)、z(ω2)估计出的内部状态量rm、ra、ca来使其确定。即，能够更可靠地进行多种内部状态量rm、ra、ca、rc、cc彼此的区分。

此外，反之也可以是，内部状态量估计单元基于低频阻抗z(ωl)来估计某个内部状态量，基于估计出的内部状态量以及高频阻抗z(ωh)、z(ω1)、z(ω2)来估计其它内部状态量。

另外，根据本实施方式，上述高频带(阳极响应频带和电解质膜响应频带)包括阳极响应频带和电解质膜响应频带，该阳极响应频带是对燃料电池1的阳极112的状态量ra、ca表现出响应性的频带，该电解质膜响应频带是频率比阳极响应频带高的频带，对燃料电池1的电解质膜的状态量rm表现出响应性。而且，阻抗获取单元获取基于从阳极响应频带选择出的频率的阳极响应阻抗z(ω1)、z(ω2)以及基于从电解质膜响应频带选择出的频率的电解质膜响应阻抗z(ωh)这两方，来作为高频阻抗z(ωh)、z(ω1)、z(ω2)(步骤s101、步骤s103)。

由此，能够基于电解质膜响应阻抗z(ωh)和阳极响应阻抗z(ω1)、z(ω2)来分别估计燃料电池1的电解质膜111的状态量rm和阳极112的状态量ra、ca。

并且，根据本实施方式，内部状态量估计单元基于电解质膜响应阻抗z(ωh)来估计电解质膜111的状态量rm(步骤s102)，基于估计出的该电解质膜111的状态量rm以及阳极响应频带阻抗z(ω1)、z(ω2)来估计阳极112的状态量ra、ca(步骤s104)。

由此，能够基于估计出的该电解质膜111的状态量rm以及阳极响应频带阻抗z(ω1)、z(ω2)，更可靠地与其它状态量分开地估计阳极112的状态量ra、ca。

特别是，在本实施方式中，阳极112的状态量ra、ca包含该阳极112的反应电阻值ra和双电层电容值ca，阴极113的状态量rc、cc包含该阴极113的反应电阻值rc和双电层电容值cc。而且，内部状态量估计单元基于阳极响应阻抗z(ω1)、z(ω2)来估计阳极112的反应电阻值ra和阳极112的双电层电容值ca(步骤s104)。另外，内部状态量估计单元基于估计出的电解质膜111的状态量rm、阳极112的反应电阻值ra和阳极112的双电层电容值ca以及低频阻抗z(ωl)来估计阴极113的反应电阻值rc(步骤s106)。

据此，对于包含除阴极113的反应电阻值rc以外的全部信息的低频带的低频阻抗z(ωl)，能够应用基于阳极响应阻抗(z(ω1)、z(ω2))估计出的阳极112的反应电阻值ra和双电层电容值ca、基于电解质膜响应阻抗z(ωh)估计出的电解质膜111的状态量rm。

因而，能够根据包含除目标的状态量rc以外的信息的低频带的低频阻抗z(ωl)，来适当地分开估计该目标的状态量rc。

(第二实施方式)

下面，说明第二实施方式。此外，对与已经说明的第一实施方式的要素相同的要素标注同一标记。

图6是表示第二实施方式所涉及的状态量的估计的流程的流程图。图6中的步骤s101～步骤s104与图5中的s101～步骤s104相同，因此省略其详细说明。在第二实施方式中，将预先设定的燃料电池1的i-v(电流电压)特性线图(i-v特性图)中的特性线的直线部分的斜率视作低频阻抗来获取该斜率，以代替在低频带的频率下测量低频阻抗。

如图所示，在经过步骤s101～步骤s104之后、即获取到阳极112的反应电阻值ra和双电层电容值ca的估计值之后，在步骤s205中，将燃料电池1的i-v特性图中的特性线的直线部分的斜率δv/δi视作低频阻抗z(ωl)，获取该斜率δv/δi。

图7中分别示出稳定时和非稳定时的燃料电池1的i-v特性线。此外，该燃料电池1的i-v特性线是预先基于实验等而决定的。特性线cv1表示稳定时的i-v特性，特性线cv2表示非稳定时的i-v特性。在此，稳定时是指不是起步时、停车时等突然加速状态的、稳定行驶时的燃料电池1的输出特性。

特别是，根据图可以了解的是，在稳定时的特性线cv1的稳定区域p，斜率δv/δi的变动小而大致取固定值，呈直线形态。因而，在稳定区域p，不管输出电流i如何都能够将斜率δv/δi视作固定值。

这样，δv/δi的值固定的稳定区域p是使稳定时的特性线cv1的δv/δi的值为规定值以下的横轴(输出电流i)的区间。

在本实施方式中，控制器6将该稳定区域p的δv/δi的值预先存储到未图示的存储器等，在获取低频阻抗z(ωl)的定时从该存储器读出δv/δi的值，视作低频阻抗z(ωl)。这样得到的低频阻抗z(ωl)与现实的值良好匹配。

然后，在步骤s206中，使用作为低频阻抗z(ωl)而获取的δv/δi的值来进行阴极113的反应电阻值rc的估计。

具体地说明。当在上述的式(1)中假定ω为低频率(ω→0)时，认为下式成立。

[式14]

因而，在式(14)中，当将阻抗z置换为δv/δi时，为下式。

[式15]

由此，通过将通过步骤s101～步骤s104的过程而估计出的电解质膜电阻值rm和阳极112的反应电阻值ra代入到式(15)，能够计算出阴极113的反应电阻值rc。

根据以上说明的本实施方式所涉及的燃料电池1的状态检测装置，作为阻抗获取单元的控制器6获取燃料电池1的i-v特性线的斜率δv/δi来作为低频阻抗z(ω1)。即，不直接测量就能够获取低频阻抗z(ω1)。

此外，例如，也可以是，利用获取i-v特性线的斜率δv/δi的值来作为低频阻抗z(ω1)以及通过测量来获取低频阻抗z(ω1)这两方的方法，来获取低频阻抗z(ω1)，将对利用这两方的方法得到的低频阻抗z(ω1)进行相互比较/校正等而获取到的更高精度的低频阻抗z(ω1)用在阴极113的反应电阻值rc的估计中。

并且，在本实施方式中，作为阻抗获取单元的控制器6在燃料电池1的i-v特性线cv1的斜率的值的变动为规定值以下的稳定区域p获取斜率δv/δi来作为低频阻抗z(ω1)。

在像这样斜率δv/δi的变动比较小的稳定区域p，无论输出电流i的测量值如何，都可以将斜率δv/δi的值视作固定，因此无需针对每个输出电压v、输出电流i的测量值来计算斜率δv/δi的值，能够减少运算量。

(第三实施方式)

说明第三实施方式。此外，对与已说明的实施方式的要素相同的要素标注同一标记。

图8是表示本实施方式所涉及的状态量的估计的流程的流程图。如图所示，在本实施方式中，省略了与图5所示的步骤s101及步骤s102相当的使用电解质膜响应频带的频率进行的电解质膜电阻值rm的估计。

特别是，在本实施方式中，在特有的步骤s304中，使用在阳极响应频带的两个点的频率ω1、ω2获取到的阳极响应阻抗z(ω1)和z(ω2)，来估计作为状态量的阳极112的反应电阻值ra、阳极112的双电层电容值ca、阴极113的双电层电容值cc以及电解质膜电阻值rm(步骤s304)。

下面，说明步骤s304中的状态量估计的一个方式。

在本实施方式中，也基于上述的阻抗的式(2)来进行计算。取式(2)的实部来得到式(3)、基于式(3)得到式(4)的步骤与第一实施方式所涉及的估计阳极112的反应电阻值ra和阳极112的双电层电容值ca的情况相同。

然后，如果对式(4)进行变形，则能够得到下式。

[式16]

此外，如上所述，mr是将两个阻抗z(ω1)和z(ω2)连接的直线的斜率，是已知的值。

另一方面，当取式(2)的虚部时为下式。

[式17]

在此，当将式(16)的ra代入到上述式(17)、使两边乘以ω时，为下式。

[式18]

然后，当将上述已知的频率ω1和ω2以及与它们对应的阻抗测量值的虚数成分zi1和zi2分别代入到式(18)来得到两个式子、取这两个式子之差来消除阴极的双电层电容cc时，能够得到与作为未知数的阳极的双电层电容ca有关的4次方程式。

[式19]

对式(19)的4次方程式进行求解，考虑到ca不能取虚数值，作为阳极的双电层电容ca的候选，能够得到两个解。

[式20]

[式21]

此外，式(19)的4次方程式能够使用本领域技术人员所知的各种解法。

其中，t1是如下所述那样定义的常数。

[式22]

以上，说明了本发明的实施方式，但是上述实施方式不过示出了本发明的应用例的一部分，其宗旨并不在于将本发明的保护范围限定于上述实施方式的具体结构。

并且，式中的a2、a1以及a0分别为下式。

[式23]

并且，通过将ca1和ca2分别代入到上述式(16)，能够与该ca1和ca2对应地决定ra1和ra2作为反应电阻的估计值的候选。估计值的候选ra1和ra2如下述那样。

[式24]

[式25]

在此，需要根据上述的阳极112的双电层电容值的候选ca1和ca2以及反应电阻值的候选ra1和ra2来决定适于现实特性的真正的估计值。说明其方法的一例。

在本实施方式中，对于该真正的估计值的决定，并非仅根据ca1、ra1、ca2以及ra2的值来判断，而是对上述式(17)中的阻抗虚部的式子进行变形而得到的阴极113的双电层电容值cc的式子。

[式26]

图9中示出了阴极113的双电层电容值的候选cc1、cc2的频率响应。此外，该图表是基于预先通过实验等计算出的使频率ω1和ω2在阳极响应频带的范围连续变化而得到的双电层电容值的候选cc1、cc2的数据的图表。

此外，在该图表中，将cc1描绘的线表示为虚线，将cc2描绘的线表示为实线。另外，频率ωd是对于阳极112的反应电阻值和双电层电容值的候选的组(ca1、ra1)和(ca2、ra2)而言(ca1、ra1)＝(ca2、ra2)的频率。即，在频率ωd下，ca1、ra1、ca2以及ra2中的上述式(20)、(21)、(24)以及(25)的根号内部为0。

如图所示，在频率ω<ωd的区域，双电层电容值的估计值候选cc2基本取0以下的值，在ωd之前cc2的值相对于频率变化极端敏感，因此在频率ω<ωd的区域，cc1是应该现实采用的真正的估计值。

因而，关于阴极113的双电层电容值和反应电阻值，在频率ω<ωd的区域，也分别采用与上述cc1对应的ca1和ra1。

另一方面，在ω>ωd的区域，仅观察阴极113的双电层电容值的候选(cc1、cc2)的变化是难以判断应该采用cc1和cc2中的哪一个的。因此，通过直接研究阳极112的反应电阻值和双电层电容值的候选的组(ca1、ra1)和(ca2、ra2)来进行该判断。

图10a示出了阳极112的双电层电容的候选ca1、ca2的频率响应。另外，图10b示出了阳极112的反应电阻值的候选ra1、ra2的频率响应。此外，这些图表也是基于预先通过实验等计算出的使频率ω1和ω2在阳极响应频带的范围连续变化而得到的候选的组(ca1、ra1)和(ca2、ra2)的数据的图表。

当参照图10a时，在ω>ωd的区域，阳极112的双电层电容值的候选ca1的频率极端敏感。因而，在ω>ωd的区域，作为阳极112的双电层电容值的真正的估计值，ca2是应该现实采用的值。因而，在频率ω>ωd的区域，应该分别采用ca2以及与其对应的ra1。

此外，参照图10b可以了解的是，在比频率ωd小的ω<ωd的区域，反应电阻值的候选ra2对于频率变化极端敏感，因此判断出反应电阻值的候选ra1是应该现实采用的真正的估计值。因而，在该频率ω<ωd的区域，应该分别采用与ra1对应的ca1和ra1，可知这一点与基于阴极113的双电层电容值的频率响应进行的考察一致。

另外，在频率ω＝ωd时，(ca1、ra1)＝(ca2、ra2)，因此可以将这些候选的组中的任一组采用为真正的候选的组。

基于以上的考察可知，在决定真正的估计值时，应该从候选的组(ca1、ra1)和(ca2、ra2)中决定的对象是根据频率而变化的。具体地说，根据阳极响应频带的两个点的频率ω1、ω2以及频率ωd的大小，从候选的组(ca1、ra1)和(ca2、ra2)决定适当的组。并且，如果将所决定的阳极112的双电层电容值ca和反应电阻值ra的估计值代入到式(3)，则由于频率ω和阻抗测量值的实部zr已知而可以求出电解质膜电阻值rm。

使用这样求出的阳极112的双电层电容值ca、反应电阻值ra以及电解质膜电阻值rm的估计值，与第一实施方式同样地进行以后的步骤s105和步骤s106，将阴极113的反应电阻值rc也估计出来。

根据以上说明的本实施方式所涉及的燃料电池1的状态判定，作为阻抗获取单元和内部状态量估计单元的控制器6仅获取阳极响应阻抗z(ω1)和z(ω2)来作为高频阻抗，基于阳极响应频带阻抗z(ω1)和z(ω2)来估计阳极112的状态量ca和ra。

由此，尽管省略基于电解质膜响应阻抗的测量的电解质膜电阻值rm的估计来减轻对控制器6的负荷，也能够估计阳极112的状态量ca和ra，最终能够将作为阴极113的状态量的反应电阻值rc也估计出来。

(第四实施方式)

说明第四实施方式。此外，对与已说明的实施方式的要素相同的要素标注同一标记。

图11是表示本实施方式所涉及的状态量的估计的流程的流程图。如图所示，在本实施方式中，与第三实施方式同样地，在步骤s103求出阳极响应阻抗z(ω1)、z(ω2)，在步骤s304中求出阳极112的反应电阻值ra和双电层电容值ca、阴极113的双电层电容值cc以及电解质膜电阻值rm的估计值。

之后，与第二实施方式的情况同样地，在步骤s205中基于燃料电池1的i-v特性来获取低频阻抗δv/δi，在步骤s206中根据这样获取到的低频阻抗δv/δi以及电解质膜电阻值rm的估计值来估计阴极113的反应电阻值rc。

因而，根据本实施方式所涉及的燃料电池1的状态判定，不直接测量低频阻抗z(ω1)就能够估计，并且能够省略基于电解质膜响应阻抗的测量进行的电解质膜电阻值rm的估计，因此能够进一步减轻对控制器6的负荷。

(第五实施方式)

说明第五实施方式。此外，对与已说明的实施方式的要素相同的要素标注同一标记。

在本实施方式中，在第二实施方式和第四实施方式所涉及的步骤s205中，取代将图7的稳定时的特性线cv1的稳定区域p的δv/δi的值事先存储到存储器的方式，而使用实际的输出电压v和输出电流i的测量值以计算δv/δi的值。

图12示出了稳定时的燃料电池1的i-v特性线。特别是，在本实施方式中，对于在规定的测量定时由电流传感器51测定出的输出电流i1、i2以及在该定时由电压传感器52测定出的输出电压v1、v2，通过计算-(v1-v2)/(i1-i2)来计算斜率δv/δi。

即，根据输出电压和输出电流的测量值来决定视作低频阻抗的斜率δv/δi。

在本实施方式中，像这样基于电流和电压的两组测量值(i1、v1)、(i2、v2)来计算燃料电池1的i-v特性线的斜率δv/δi。由此，与使用在稳定区域p视作固定值而决定的斜率δv/δi的情况相比，能够得到更高精度地反映出实际的特性的δv/δi的值。作为结果，将该δv/δi的值视作低频阻抗来计算出的阴极113的反应电阻值rc的估计值的精度也提高。

(第六实施方式)

说明第六实施方式。此外，对与已说明的实施方式的要素相同的要素标注同一标记。

在本实施方式中，取代如第五实施方式那样测量输出电流和输出电压的两组测量值(i1、v1)、(i2、v2)以求出i-v特性线的斜率δv/δi，而是使用输出电流和输出电压的一个测量值(i3、v3)以及事先设定的一个点(iset、vset)来进行i-v特性线的斜率δv/δi的计算。

图13是说明用于进行i-v特性线的斜率δv/δi的计算的一组电流和电压的设定方法的一例的图。此外，在该图中为了使附图清楚，以虚线表示稳定时的特性线cv1。如图所示，在本实施方式中，图中涂黑四方形所表示的点与上述(iset、vset)相当。特别是，iset＝0。

因而，根据上述测量值(i3、v3)和事先设定值(iset、vset)，通过计算-(vset-v3)/(iset-i3)来计算斜率δv/δi的值。

如上所述，根据本实施方式，基于电流和电压的一组测量值(i3、v3)和事先设定的一组电流和电压的值(iset、vset)来计算i-v特性线的斜率δv/δi的值。

因而，对于燃料电池1的i-v特性线的斜率δv/δi，在计算该斜率的值时使用的i-v特性线上的两个点中，作为一个点使用事先设定的(iset、vset)来抑制运算量，作为另一个点使用测量值(i3、v3)，由此能够将计算的精度也确保为一定以上。

(第七实施方式)

下面，说明第七实施方式。此外，对与已说明的实施方式的要素相同的要素标注同一标记。

在第一实施方式等中测量燃料电池1的阻抗时，测定叠加有交流信号的输出电流i和输出电压v，在本实施方式中，取代这种结构，进行所谓的激励电流施加法，即，从规定的测定用电流源向燃料电池1提供电流i，基于该供给电流i以及输出的电压v来计算阻抗z＝v/i。

图14是概要性地示出本实施方式的燃料电池系统100中阻抗测量所涉及的重要部分的框图。

如图所示，在本实施方式所涉及的燃料电池系统100中，设置有一边调整交流电流一边向燃料电池1施加该交流电流的施加交流电流调整部200。

施加交流电流调整部200除了与构成为堆的燃料电池1的正极端子(阴极侧端子)1b及负极端子(阳极侧端子)1a连接以外，还与中途端子1c连接。此外，与中途端子1c连接的部分如图所示那样接地。

而且，施加交流电流调整部200具有：正极侧电压测定传感器210，其测定正极端子1b相对于中途端子1c的正极侧交流电位差v1；以及负极侧电压测定传感器212，其测定负极端子1a相对于中途端子1c的负极侧交流电位差v2。

并且，施加交流电流调整部200具有：正极侧交流电源部214，其向包括正极端子1b和中途端子1c的电路施加交流电流i1；负极侧交流电源部216，其向包括负极端子1a和中途端子1c的电路施加交流电流i2；控制器218，其对这些交流电流i1和交流电流i2的振幅、相位进行调整；以及运算部220，其基于正极侧交流电位差v1、v2和交流电流i1、i2来进行燃料电池1的阻抗z的运算。

在本实施方式中，控制器218以使正极侧交流电位差v1与负极侧交流电位差v2相等的方式调节交流电流i1和交流电流i2的振幅和相位。此外，该控制器218也可以由图3所示的控制器6构成。

另外，运算部220包括未图示的ad变换器、微机芯片等硬件以及计算阻抗的程序等软件结构，将正极侧交流电位差v1除以交流电流i1来计算从中途端子1c到正极端子1b的阻抗z1，将负极侧交流电位差v2除以交流电流i2来计算从中途端子1c到负极端子1a的阻抗z2。并且，运算部220通过取阻抗z1与阻抗z2之和来计算燃料电池1的整体阻抗z。

根据上述的本实施方式所涉及的燃料电池的状态估计装置，能够得到以下的效果。

本实施方式所涉及的燃料电池的状态估计装置具有：交流电源部214、216，该交流电源部214、216与燃料电池1连接，向该燃料电池1输出交流电流i1、i2；作为交流调整部的控制器218，其基于作为从燃料电池1的正极侧1b的电位减去中途部分1c的电位而求出的电位差的正极侧交流电位差v1以及作为从燃料电池1的负极侧1a的电位减去中途部分1c的电位而求出的电位差的负极侧交流电位差v2，来调整交流电流i1、i2；以及阻抗运算部220，其基于调整后的交流电流i1、i2以及正极侧交流电位差v1和负极侧交流电位差v2来运算燃料电池1的阻抗z。

控制器218以使燃料电池1的正极侧的正极侧交流电位差v1与负极侧的负极侧交流电位差v2实质上一致的方式，对由正极侧交流电源部214施加的交流电流i1和由负极侧交流电源部216施加的交流电流i2的振幅和相位进行调节。由此，正极侧交流电位差v1的振幅与负极侧交流电位差v2的振幅变得相等，因此正极端子1b与负极端子1a实质上为相等电位。因而，防止用于阻抗测量的交流电流i1、i2流向负载53，因此防止对燃料电池1的发电产生影响。

另外，在燃料电池1处于发电状态的情况下执行上述阻抗测量的情况下，测量用交流电位会叠加在通过该发电产生的电压上，因此正极侧交流电位差v1和负极侧交流电位差v2的值本身变大，但是正极侧交流电位差v1和负极侧交流电位差v2的相位、振幅本身并不改变，因此能够与燃料电池1未处于发电状态的情况同样地执行高精度的阻抗测量。

以上，说明了本发明的实施方式，但是上述实施方式不过示出了本发明的应用例的一部分，其宗旨并不在于将本发明的保护范围限定于上述实施方式的具体结构。例如，各实施方式中的获取阳极响应阻抗、电解质膜响应阻抗以及低频阻抗的步骤(步骤s101、步骤s103以及步骤s105)等不限定于在各实施方式中说明的步骤顺序，能够任意地变更。

例如，也可以在将获取阳极响应阻抗、电解质膜响应阻抗以及低频阻抗的步骤全部进行完之后，进行各状态量的估计。

另外，在燃料电池1中估计多个内部状态量的方式不仅限于在上述各实施方式中说明的方式。

例如，也可以取代第一实施方式、第三实施方式中的步骤s105中的从低频带选择一个频率ωl的方式，而在低频带选择两个频率ωl1、ωl2来求出低频阻抗z(ωl1)和z(ωl2)。由此，不仅最终求出阴极113的反应电阻rc，还能够求出阴极113的双电层电容ca的估计值。

另外，燃料电池1的简易等效电路的方式也不限定于上述各实施方式中使用的方式。例如，也可以设定除了包括在上述各实施方式中说明的各极的反应电阻、双电层电容等电路元件以外、还包括扩散电阻、电子输送电阻以及离聚物电阻等其它要素的等效电路，以作为基于这些其它要素的内部状态量的扩散电阻值、电子输送电阻值以及离聚物电阻值等为估计的对象。