WO2026014142A1

WO2026014142A1 - 電池状態算出装置、プログラム、電池状態算出装置の制御方法

Info

Publication number: WO2026014142A1
Application number: PCT/JP2025/021276
Authority: WO
Inventors: 昌明北川; 和生松川; 裕基堀
Original assignee: Denso Corp
Current assignee: Denso Corp
Priority date: 2024-07-08
Filing date: 2025-06-12
Publication date: 2026-01-15
Anticipated expiration: 2027-01-08
Also published as: JP2026009576A

Abstract

電池状態算出装置（２０，１２０，２２０）は、交流駆動部（４２，５１，５６，６０，１３０，１３１）と、蓄電池（１０，１０００）の端子間電圧を検出する電圧検出回路（３１～３３）と、電流を検出する電流検出回路（４３，５２～５４）と、演算部（７０，９３）とを備える。演算部は、交流駆動部によって交流電流が流されている場合における電圧検出回路及び電流検出回路の検出値を取得し、取得した検出値に基づいて、蓄電池の複素インピーダンスを算出し、算出した蓄電池の複素インピーダンスと比較する基準電池の複素インピーダンスと、算出した蓄電池の複素インピーダンスとに基づいて、基準電池の複素インピーダンスに対する、算出した蓄電池の複素インピーダンスの相対誤差を補正する補正値（ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤ）を算出し、算出した補正値に基づいて、算出した蓄電池の複素インピーダンスを補正する。

Description

電池状態算出装置、プログラム、電池状態算出装置の制御方法

関連出願の相互参照

　本出願は、２０２４年７月８日に出願された日本出願番号２０２４－１０９５５３号に基づくもので、ここにその記載内容を援用する。

　本開示は、電池状態算出装置、プログラム、及び電池状態算出装置の制御方法に関する。

　従来、例えば特許文献１に記載されているように、蓄電池の複素インピーダンスを算出する電池状態算出装置が知られている。

国際公開第２０２３／１６２７５１号

　電池状態算出装置が適用されるシステムとして、蓄電池の正極端子と蓄電池の負極端子とを電気的に接続する電気経路を備えるものがある。

　上記装置は、電気経路及び蓄電池を含む閉回路に交流電流を流す交流駆動部と、電圧検出回路と、上記閉回路に流れる電流を検出する電流検出回路と、演算部とを備えている。電圧検出回路は、蓄電池に並列接続され、蓄電池の端子間電圧を検出する。

　演算部は、交流駆動部によって上記閉回路に交流電流が流されている場合における電圧検出回路の検出値及び電流検出回路の検出値を取得し、取得した電圧検出回路の検出値及び電流検出回路の検出値に基づいて、蓄電池の複素インピーダンスを算出する。

　蓄電池の状態（例えば劣化状態）を判定するために、算出した蓄電池の複素インピーダンスと、基準電池の複素インピーダンスとを比較することがある。このため、基準電池についても、上述した方法により複素インピーダンスが算出される。

　上記閉回路に交流電流が流れることに伴い、磁束が発生する。発生した磁束に起因して、電圧検出回路及び電圧検出回路の検出値に誤差が含まれることがある。この場合、蓄電池の複素インピーダンスに含まれる誤差と、基準電池の複素インピーダンスに含まれる誤差との差が大きくなるおそれがある。

　本開示は、蓄電池の複素インピーダンスに含まれる誤差と、基準電池の複素インピーダンスに含まれる誤差との差を低減できる電池状態算出装置、プログラム、及び電池状態算出装置の制御方法を提供することを主たる目的とする。

　本開示は、蓄電池の正極端子と前記蓄電池の負極端子とを電気的に接続する電気経路を備えるシステムに適用される電池状態算出装置において、
　前記電気経路及び前記蓄電池を含む閉回路に交流電流を流す交流駆動部と、
　前記蓄電池に並列接続され、前記蓄電池の端子間電圧を検出する電圧検出回路と、
　前記閉回路に流れる電流を検出する電流検出回路と、
　演算部と、
を備える。

　本開示において、前記演算部は、
　前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得し、取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記蓄電池の複素インピーダンスを算出し、
　算出した前記蓄電池の複素インピーダンスと比較する基準電池の複素インピーダンスと、算出した前記蓄電池の複素インピーダンスとに基づいて、前記基準電池の複素インピーダンスに対する、算出した前記蓄電池の複素インピーダンスの相対誤差を補正する補正値を算出し、
　算出した前記補正値に基づいて、算出した前記蓄電池の複素インピーダンスを補正する。

　これにより、基準電池の複素インピーダンスに対する蓄電池の複素インピーダンスの相対誤差を低減でき、蓄電池の複素インピーダンスに含まれる誤差と、基準電池の複素インピーダンスに含まれる誤差との差を低減することができる。その結果、例えば、基準電池及び蓄電池の複素インピーダンスを比較することにより得られる電池情報（例えば、劣化状態）の精度を高めることができる。

　本開示についての上記目的およびその他の目的、特徴や利点は、添付の図面を参照しながら下記の詳細な記述により、より明確になる。その図面は、
図１は、第１実施形態に係る電池状態算出装置の全体構成図であり、図２は、電圧検出閉回路及び電流検出閉回路を示す図であり、図３は、同一の蓄電池のインピーダンスを第１，第２電池状態算出装置それぞれで測定する場合を示す図であり、図４は、インピーダンス補正処理の手順を示すフローチャートであり、図５は、基準信号の角周波数と補正値との関係を示す図であり、図６は、補正前のコールコールプロットを示す図であり、図７は、第１電池状態測定装置における振幅，位相影響の補正前後のコールコールプロットを示す図であり、図８は、第２電池状態測定装置における振幅，位相影響の補正前後のコールコールプロットを示す図であり、図９は、補正後のコールコールプロットを示す図であり、図１０は、ファラデーの法則の概念図であり、図１１は、Biot-Savartの法則の参考図であり、図１２は、第１実施形態の変形例に係る電池状態算出装置の全体構成図であり、図１３は、第２実施形態に係る電池状態算出装置の全体構成図であり、図１４は、電池ケース等の図示を省略した電池状態算出装置の全体構成図であり、図１５は、第３実施形態に係る電池状態算出装置の全体構成図であり、図１６は、電池ケース等の図示を省略した電池状態算出装置の全体構成図であり、図１７は、インピーダンス補正処理の手順を示すフローチャートであり、図１８は、第３実施形態の変形例に係る電池状態算出装置の全体構成図であり、図１９は、図１８に示すモジュール内の構成を示す図であり、図２０は、第３実施形態の変形例に係る電池状態算出装置の全体構成図であり、図２１は、第３実施形態の変形例に係る電池状態算出装置の全体構成図であり、図２２は、第３実施形態の変形例に係る組電池の平面図であり、図２３は、相互インダクタンスの計算結果の一例を示す図であり、図２４は、第３実施形態の変形例に係る組電池の平面図であり、図２５は、相互インダクタンスの計算結果の一例を示す図であり、図２６は、相互インダクタンスの計算結果の一例を示す図であり、図２７は、第４実施形態に係るインピーダンス補正処理の手順を示すフローチャートであり、図２８は、その他の実施形態に係る補正方法の概要を示す図であり、図２９は、その他の実施形態に係る補正方法の概要を示す図であり、図３０は、その他の実施形態に係る電圧，電流検出部の構成を示す図であり、図３１は、その他の実施形態に係る電気経路の配置態様を示す図であり、図３２は、その他の実施形態に係る電気回路の配置態様を示す図である。

　図面を参照しながら、複数の実施形態を説明する。複数の実施形態において、機能的に及び／又は構造的に対応する部分及び／又は関連付けられる部分には同一の参照符号、又は百以上の位が異なる参照符号が付される場合がある。対応する部分及び／又は関連付けられる部分については、他の実施形態の説明を参照することができる。

　＜第１実施形態＞
　以下、本開示に係る電池状態算出装置を具体化した第１実施形態について、図面を参照しつつ説明する。

　電池状態算出装置２０は、蓄電池１０の状態を算出する装置であり、具体的には、蓄電池１０の内部インピーダンスを計測する装置である。計測された内部インピーダンスに基づいて、蓄電池１０の蓄電状態（具体的にはＳＯＣ）及び劣化状態（ＳＯＨ）等が算出される。本実施形態において、蓄電池１０は、充放電可能な２次電池であり、例えばリチウムイオン蓄電池又はニッケル水素蓄電池である。蓄電池１０は、１つの単位電池、又は複数の単位電池の直列接続体である。単位電池は、１つの電池セル、又は複数の電池セルの直列接続体である。

　電池状態算出装置２０は、電圧処理部３０を備えている。電圧処理部３０は、集積回路（ＩＣ）で構成され、正極側電圧経路３１と、負極側電圧経路３２と、電圧検出部３３とを備えている。正極側電圧経路３１の第１端である正極電圧端子３１ａは、蓄電池１０の正極端子１０ａに接続され、正極側電圧経路３１の第２端は、電圧検出部３３に接続されている。負極側電圧経路３２の第１端である負極電圧端子３２ａは、蓄電池１０の負極端子１０ｂに接続され、負極側電圧経路３２の第２端は、電圧検出部３３に接続されている。電圧検出部３３は、正極側電圧経路３１及び負極側電圧経路３２の間の電圧差を蓄電池１０の端子間電圧として検出する。なお、本実施形態において、電圧処理部３０、正極側電圧経路３１及び負極側電圧経路３２が「電圧検出回路」に相当する。

　電池状態算出装置２０は、蓄電池１０に流れる電流を生成及び検出するための構成として、制限抵抗体４１、電流制御スイッチ４２、シャント抵抗体４３、電流処理部５０、正極側電流経路５３、負極側電流経路５４及び基準信号生成部６０を備えている。

　制限抵抗体４１の第１端には、正極接続部４０ａの第１端が接続されている。正極接続部４０ａの第２端には、蓄電池１０の正極端子１０ａが接続されている。制限抵抗体４１の第２端には、電流制御スイッチ４２を介してシャント抵抗体４３の第１端が接続されている。本実施形態において、電流制御スイッチ４２は、ＮチャネルＭＯＳＦＥＴである。シャント抵抗体４３の第２端には、負極接続部４０ｂの第１端が接続されている。負極接続部４０ｂの第２端には、蓄電池１０の負極端子１０ｂが接続されている。本実施形態において、正極接続部４０ａの第２端は、正極端子１０ａが着脱可能に接続される構成になっており、負極接続部４０ｂの第２端は、負極端子１０ｂが着脱可能に接続される構成になっている。

　なお、正極接続部４０ａの第２端に正極側電圧経路３１が接続され、負極接続部４０ｂの第２端に負極側電圧経路３２が接続されていてもよい。この場合、正極接続部４０ａの第２端が正極端子１０ａに接続されることに伴い、正極側電圧経路３１が正極端子１０ａに電気的に接続される。また、負極接続部４０ｂの第２端が負極端子１０ｂに接続されることに伴い、負極側電圧経路３２が負極端子１０ｂに電気的に接続される。ちなみに、正極接続部４０ａと負極接続部４０ｂとの間における制限抵抗体４１、電流制御スイッチ４２及びシャント抵抗体４３の並び順は、図１に示す並び順に限らない。

　電流処理部５０は、集積回路（ＩＣ）で構成され、電流変調部５１と、電流検出部５２とを備えている。

　電池状態算出装置２０は、正極側電流経路５３、負極側電流経路５４及び基準信号生成部６０を備えている。正極側電流経路５３は、シャント抵抗体４３の第１端と電流検出部５２とを接続する。負極側電流経路５４は、シャント抵抗体４３の第２端と電流検出部５２とを接続する。電流検出部５２は、正極側電流経路５３及び負極側電流経路５４の間の電圧差を蓄電池１０に流れる電流として検出する。電流検出部５２により検出された電流は、電流変調部５１に入力される。

　電流変調部５１は、測定対象である蓄電池１０に所定の交流信号を入力する回路である。基準信号生成部６０は、正弦波の基準信号を生成し、生成した基準信号を電流変調部５１に入力する。電流変調部５１は、電流検出部５２により検出された電流を、基準信号の周波数で変動する正弦波の目標電流に制御すべく、電流制御スイッチ４２のゲート電圧を制御する。交流電流は、例えば、蓄電池１０の充電電流又は放電電流である。なお、目標電流は、正弦波信号に変えて、矩形波信号であってもよい。

　ちなみに、本実施形態において、電流制御スイッチ４２、電流変調部５１及び基準信号生成部６０が「交流駆動部」に相当する。また、シャント抵抗体４３、電流検出部５２、正極側電流経路５３及び負極側電流経路５４が「電流検出回路」に相当する。

　図１では、蓄電池１０の複素インピーダンスの等価回路モデルが示されている。この等価回路モデルでは、蓄電池１０の内部複素インピーダンスは、オーミック抵抗１１と反応抵抗１２との直列接続体により構成されている。オーミック抵抗１１は、蓄電池１０を構成する電極及び電解液における通電抵抗である。また、反応抵抗１２は、電極における電極界面反応による抵抗を表すものであり、抵抗成分１２ａと容量成分１２ｂとの並列接続体として表される。

　電池状態算出装置２０は、温度センサ８０及び圧力センサ８１を備えている。温度センサ８０は、蓄電池１０の温度を検出する。圧力センサ８１は、蓄電池１０の圧力（例えば内圧）を検出する。

　電池状態算出装置２０は、制御装置７０を備えている。制御装置７０は、電子制御装置（ＥＣＵ：Electronic Control Unit）であり、ハードウェアとしてのプロセッサ７１と、記憶部７２とを備えている。制御装置７０において、プロセッサ７１と記憶部７２とは通信バス７３を介して互いに接続されている。

　記憶部７２は、ハードウェアとしてのメモリ及びストレージを含む。メモリは、制御装置７０の処理に使用されるデータを記憶するための記憶装置である。メモリは、例えば、プロセッサ７１が処理を行う場合に一時的に使用するための作業領域をプロセッサ７１に提供する。メモリは、例えば、ＲＯＭ又はＲＡＭを含む。ストレージは、プロセッサ７１が読み込んで実行するための各種プログラム及びデータを記憶する記憶装置であり、非遷移的実体的記録媒体（non-transitory tangible storage medium）である。ストレージは、例えば、ＨＤＤ又はフラッシュメモリを含む。ストレージには、後述する図４等の処理のためのプログラム情報等が記憶されている。

　なお、例えば、記憶部７２には、非遷移的実体的記録媒体に記憶されたプログラム情報がインストールされる。記録媒体は、例えば、ＵＳＢメモリ、ＣＤ－ＲＯＭ又はＤＶＤである。また、例えば、記憶部７２には、例えばＯＴＡ（Over The Air）等、通信ネットワークを介して送信されたプログラム情報がインストールされる。

　制御装置７０には、電圧検出部３３、電流検出部５２、温度センサ８０及び圧力センサ８１の検出値が入力される。制御装置７０は、基準信号生成部６０に、基準信号の生成指令を入力する。また、「演算部」として機能する制御装置７０は、電圧検出部３３及び電流検出部５２の検出値に基づいて、蓄電池１０の複素インピーダンス（Ｚｍ）の情報を算出する。つまり、制御装置７０は、電圧検出部３３及び電流検出部５２の検出値に基づいて、複素インピーダンスの実部（ＲｅＺｍ）、及び複素インピーダンスの虚部（ＩｍＺｍ）を算出する。制御装置７０は、算出結果に基づいて、例えば、複素インピーダンス平面プロット（コールコールプロット）を作成し、電極及び電解質などの特性を把握する。例えば、制御装置７０は、蓄電状態（ＳＯＣ）や劣化状態（ＳＯＨ）を把握する。

　図２には、電圧検出閉回路３４及び電流検出閉回路５５を示す。電圧検出閉回路３４は、正極側電圧経路３１、蓄電池１０の正極端子１０ａから負極端子１０ｂに至る経路、負極側電圧経路３２及び電圧検出部３３を含む閉回路である。電流検出閉回路５５は、シャント抵抗体４３、正極側電流経路５３、電流検出部５２及び負極側電流経路５４を含む閉回路である。

　電流変調部５１の制御により蓄電池１０に交流電流が流れると、磁場が発生する。発生した磁場が電圧検出閉回路３４を鎖交すると、電圧検出閉回路３４に誘導起電力が発生する。この場合、正極側電圧経路３１及び負極側電圧経路３２を介して電圧検出部３３が検出する電圧に検出誤差が生じることとなる。

　また、交流電流が流れることに伴い発生した磁場が電流検出閉回路５５を鎖交すると、電流検出閉回路５５に誘導起電力が発生する。この場合、正極側電流経路５３及び負極側電流経路５４を介して電流検出部５２が検出する電流に検出誤差が生じることとなる。電圧検出部３３及び電流検出部５２の検出値に検出誤差が含まれると、蓄電池１０の複素インピーダンスの算出精度が低下する懸念がある。

　検出誤差の影響は、電池状態算出装置２０を構成する各要素の相対位置に依存する。同じ蓄電池であったとしても、相対位置が異なる２つの電池状態算出装置それぞれにおいて算出された複素インピーダンスは異なることがある。例えば、図３に示すように、第１の組織（例えば、研究機関）が所有する第１電池状態算出装置２０Ａと、第１の組織とは異なる第２の組織（例えば、メーカー企業）が所有する第２電池状態算出装置２０Ｂとで、構成要素の相対位置が異なっているとする。この場合、第１電池状態算出装置２０Ａ及び第２電池状態算出装置２０Ｂの測定対象となる蓄電池が同じ蓄電池１０であったとしても、第１電池状態算出装置２０Ａにより算出された複素インピーダンスと、第２電池状態算出装置２０Ｂにより算出された複素インピーダンスとが異なってしまう。この場合、各組織において算出された複素インピーダンスを適正に比較することができない。そこで、本実施形態では、適正に比較できる技術が採用されている。

　まず、比較するための原理について説明する。

　蓄電池１０に実際に流れる電流変動Ｉｓを下式（ｅｑ１）で表す。下式（ｅｑ１）において、Ｉ０は、蓄電池１０に実際に流れる電流の振幅を示し、Δαは、基準信号生成部６０の基準信号に対する電流Ｉｓの位相差を示し、ｊは虚数単位を示す。

　蓄電池１０に電流Ｉｓが流れる場合に発生する蓄電池１０の実際の端子間電圧変動を下式（ｅｑ２）で表す。下式（ｅｑ２）において、Ｚｂは、蓄電池１０の真の複素インピーダンスを示す。

　電圧検出部３３によって検出された蓄電池１０の端子間電圧変動Ｖｍを下式（ｅｑ３）で表す。下式（ｅｑ３）において、ω０は、基準信号の角周波数を示し、Ｍｖは、メイン閉回路と、電流検出閉回路５５との相互インダクタンスを示す。メイン閉回路は、蓄電池１０の正極端子１０ａから負極端子１０ｂに至る経路、正極接続部４０ａ、制限抵抗体４１、電流制御スイッチ４２、シャント抵抗体４３及び負極接続部４０ｂを含む閉回路である。Δβは、基準信号に対する端子間電圧変動Ｖｍの位相差を示す。

　電流検出部５２によって検出された蓄電池１０の電流変動Ｉｍを下式（ｅｑ４）で表す。下式（ｅｑ４）において、Ｍｉは、上記メイン閉回路と電流検出閉回路５５との相互インダクタンスを示す。Ｒｓは、シャント抵抗体４３の抵抗値を示し、Δγは、基準信号に対する電流変動Ｉｍの位相差を示す。

　上式（ｅｑ４）をＩｍについて解くと、下式（ｅｑ５）が導かれる。

　検出された電圧，電流変動Ｖｍ，Ｉｍに基づいて算出される蓄電池の複素インピーダンスＺｍは、下式（ｅｑ６）により表される。

　上式（ｅｑ６）は、相互インダクタンスの影響と、位相遅延の影響とを含む補正式となる。上式（ｅｑ６）をＺｂについて解くと、下式（ｅｑ７）が導かれる。

　検出された電流，電圧変動Ｉｍ，Ｖｍに基づく複素インピーダンスＺｍに、上式（ｅｑ８）に示す振幅影響係数Ｋａｍｐ、上式（ｅｑ９）に示す位相影響係数Ｋｐｈ及び上式（ｅｑ１０）に示す切片Ｉｎｔｓｐｔを作用させた値が真の複素インピーダンスＺｂとなる。

　なお、相互インダクタンスＭｉは、例えば、磁場解析によるシミュレーションにより算出、θｐ若しくはＫａｍｐを２つ以上のＲｓに対して算出するとともに連立方程式を作成して算出、又はフィッティングにより算出することができる。

　ホワイトノイズ又は不要な信号に起因した誤差の場合、Ｓ／Ｎ比を高くするために測定電流を大きくすれば、誤差を低減できる。これに対し、振幅影響係数Ｋａｍｐ、位相影響係数Ｋｐｈ及び切片Ｉｎｔｓｐｔは、いずれも測定電流に依存しない。つまり、測定電流を大きくすれば低減できる誤差ではない。このため、後述するインピーダンス補正処理が必要となる。

　以下、真の複素インピーダンスがＺｔである同一の蓄電池１０の状態を第１電池状態算出装置２０Ａ及び第２電池状態算出装置２０Ｂそれぞれで測定する場合について説明する。これまでに登場した数式において、第１電池状態算出装置２０Ａにおける検出値，算出値に添え字Ａを付し、第２電池状態算出装置２０Ｂにおける検出値，算出値に添え字Ｂを付す。

　上式（ｅｑ７）を参照すると、第１電池状態算出装置２０Ａにおいて下式（ｅｑ１１）が成立し、第２電池状態算出装置２０Ｂにおいて下式（ｅｑ１２）が成立する。

　上式（ｅｑ１１），（ｅｑ１２）から下式（ｅｑ１３）が導かれる。

　上式（ｅｑ１３）は、振幅影響係数Ｋａｍｐ及び位相影響係数Ｋｐｈを用いると、下式（ｅｑ１４）になる。

　上式（ｅｑ１４）の左辺について説明する。ＺｍＢは、第２電池状態算出装置２０Ｂにおける電圧検出部３３及び電流検出部５２の検出値に基づいて算出される複素インピーダンスである。ＫａｍｐＢは、上式（ｅｑ８）に示したように、基準信号の角周波数ω０、メイン閉回路と電流検出閉回路５５との相互インダクタンスＭｉ、及びシャント抵抗体４３の抵抗値Ｒｓから定まる振幅影響係数である。ＫｐｈＢは、上式（ｅｑ９），（ｅｑ７）に示したように、基準信号の角周波数ω０、メイン閉回路と電流検出閉回路５５との相互インダクタンスＭｉ、シャント抵抗体４３の抵抗値Ｒｓ、及び各位相差Δβ，Δγから定まる位相影響係数である。つまり、振幅影響係数ＫａｍｐＢ及び位相影響係数ＫｐｈＢは、第２電池状態算出装置２０Ｂによって定まる係数である。

　上式（ｅｑ１４）の右辺について説明する。ＺｍＡは、第１電池状態算出装置２０Ａにおける電圧検出部３３及び電流検出部５２の検出値に基づいて算出される複素インピーダンスである。ＫａｍｐＡは、上式（ｅｑ８）に示したように、基準信号の角周波数ω０、メイン閉回路と電流検出閉回路５５との相互インダクタンスＭｉ、及びシャント抵抗体４３の抵抗値Ｒｓから定まる振幅影響係数である。ＫｐｈＡは、上式（ｅｑ９），（ｅｑ７）に示したように、基準信号の角周波数ω０、メイン閉回路と電流検出閉回路５５との相互インダクタンスＭｉ、シャント抵抗体４３の抵抗値Ｒｓ、及び各位相差Δβ，Δγから定まる位相影響係数である。つまり、振幅影響係数ＫａｍｐＡ及び位相影響係数ＫｐｈＡは、第１電池状態算出装置２０Ａによって定まる係数である。

　「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」と「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」との差が、第１，第２電池状態算出装置２０Ａ，２０Ｂの装置構成（例えば、上記相対位置）の相違に起因した複素インピーダンスの誤差となる。「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」と「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」とのうち、いずれか一方が他方に近づけられる補正が実施されることにより、２つの装置構成の相違に起因した複素インピーダンスの相対誤差を低減できる。

　図４は、制御装置７０により実行されるインピーダンス補正処理のフローチャートである。

　ステップＳ１０～Ｓ１３は、第１電池状態算出装置２０Ａの制御装置７０において実行される処理である。

　まず、第１電池状態算出装置２０Ａの正極接続部４０ａ及び負極接続部４０ｂに蓄電池１０の正極端子１０ａ及び負極端子１０ｂが接続される。なお、第１電池状態算出装置２０Ａは、蓄電池１０が収容される恒温槽を備えていてもよい。

　正極接続部４０ａ及び負極接続部４０ｂに蓄電池１０が接続されると、第１電池状態算出装置２０Ａにおいて、蓄電池１０、正極接続部４０ａ、負極接続部４０ｂ、制限抵抗体４１、電流制御スイッチ４２、シャント抵抗体４３、負極接続部４０ｂ、電圧検出閉回路３４及び電流検出閉回路５５の相対位置が固定された状態になる。

　ステップＳ１０において、温度条件が成立しているか否かを判定する。温度条件は、温度センサ８０の検出温度Ｔｂｒが目標温度Ｔｔｇｔであるとの条件、又は検出温度Ｔｂｒが目標温度Ｔｔｇｔと同等であるとの条件である。検出温度Ｔｂｒが目標温度Ｔｔｇｔと同等であるとは、例えば、「Ｔｔｇｔ－５℃≦Ｔｂｒ≦Ｔｔｇｔ＋５℃」、「Ｔｔｇｔ－３℃≦Ｔｂｒ≦Ｔｔｇｔ＋３℃」、「Ｔｔｇｔ－１℃≦Ｔｂｒ≦Ｔｔｇｔ＋１℃」、又は「Ｔｔｇｔ－０．５℃≦Ｔｂｒ≦Ｔｔｇｔ＋０．５℃」のことである。

　温度条件が成立していると判定した場合には、ステップＳ１１に進み、電圧条件が成立しているか否かを判定する。電圧条件は、蓄電池１０に電流が流れていない状態において、電圧検出部３３の検出電圧Ｖｂｒが目標電圧Ｖｔｇｔ（例えば、蓄電池１０の定格電圧）であるとの条件、又は検出電圧Ｖｂｒが目標電圧Ｖｔｇｔと同等であるとの条件である。検出電圧Ｖｂｒが目標電圧Ｖｔｇｔと同等であるとは、例えば、「０．９９７×Ｖｔｇｔ≦Ｖｂｒ≦１．００３×Ｖｔｇｔ」、「０．９９９×Ｖｔｇｔ≦Ｖｂｒ≦１．００１×Ｖｔｇｔ」、又は「０．９９９７×Ｖｔｇｔ≦Ｖｂｒ≦１．０００３×Ｖｔｇｔ」のことである。

　なお、電圧条件において、検出電圧Ｖｂｒに代えて、蓄電池１０のＳＯＣが用いられてもよい。また、温度，電圧条件に加え、蓄電池１０の圧力条件が追加されてもよい。圧力条件は、圧力センサ８１の検出圧力が目標圧力との条件、又は検出圧力が目標圧力と同等であるとの条件である。

　ステップＳ１１において電圧条件が成立していると判定した場合には、ステップＳ１２に進み、第１電池状態算出装置２０Ａが備える電流変調部５１によって電流制御スイッチ４２が制御されるように、基準信号生成部６０に基準信号の生成指令を入力する。ここでは、基準信号の角周波数ω０を徐々に変化（漸増又は漸減）させる周波数スイープを行う。これにより、基準信号の周波数ω０で変動する交流電流が蓄電池１０に流れる。交流電流が流れている状態で、電流検出部５２の検出電流ＩｍＡを取得し、電圧検出部３３の検出電圧ＶｍＡを取得する。詳しくは、角周波数に対応する検出電流ＩｍＡ及び検出電圧ＶｍＡを取得する。

　ステップＳ１３では、取得した検出電流ＩｍＡ及び検出電圧ＶｍＡに基づいて、蓄電池１０の複素インピーダンスＺｍＡ（＝ＶｍＡ／ＩｍＡ）を算出する。その後、算出した複素インピーダンスＺｍＡに、第１電池状態算出装置２０Ａにおける振幅影響係数ＫａｍｐＡ及び位相影響係数ＫｐｈＡを乗算することにより、「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」を算出する。

　振幅影響係数ＫａｍｐＡ及び位相影響係数ＫｐｈＡは、例えば、実験又は計算により予め算出された値であり、記憶部７２に記憶されている。制御装置７０は、「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」の算出に用いる振幅影響係数ＫａｍｐＡ及び位相影響係数ＫｐｈＡを、記憶部７２から読み出して用いる。

　ステップＳ１４～Ｓ１９は、第２電池状態算出装置２０Ｂの制御装置７０において実行される処理である。

　第１電池状態算出装置２０Ａの測定対象になった蓄電池１０が、第１の組織から第２の組織へと運ばれる。運ばれた蓄電池１０の正極端子１０ａ及び負極端子１０ｂに、第２電池状態算出装置２０Ｂの正極接続部４０ａ及び負極接続部４０ｂが接続される。なお、第２電池状態算出装置２０Ｂは、蓄電池１０が収容される恒温槽を備えていてもよい。

　正極接続部４０ａ及び負極接続部４０ｂに蓄電池１０が接続されると、第２電池状態算出装置２０Ｂにおいて、蓄電池１０、正極接続部４０ａ、負極接続部４０ｂ、制限抵抗体４１、電流制御スイッチ４２、シャント抵抗体４３、負極接続部４０ｂ、電圧検出閉回路３４及び電流検出閉回路５５の相対位置が固定された状態になる。

　ステップＳ１４において、温度条件が成立しているか否かを判定する。温度条件は、ステップＳ１０の条件と同じ条件である。

　温度条件が成立していると判定した場合には、ステップＳ１５に進み、電圧条件が成立しているか否かを判定する。電圧条件は、ステップＳ１１の条件と同じ条件である。

　電圧条件が成立していると判定した場合には、ステップＳ１６に進み、第２電池状態算出装置２０Ｂが備える電流変調部５１によって電流制御スイッチ４２が制御されるように、基準信号生成部６０に基準信号の生成指令を入力する。ここでは、基準信号の角周波数ω０を徐々に変化（漸増又は漸減）させる周波数スイープを行う。これにより、基準信号の周波数ω０で変動する交流電流が蓄電池１０に流れる。交流電流が流れている状態で、電流検出部５２の検出電流ＩｍＢを取得し、電圧検出部３３の検出電圧ＶｍＢを取得する。詳しくは、角周波数に対応する検出電流ＩｍＢ及び検出電圧ＶｍＢを取得する。

　ステップＳ１７では、検出電流ＩｍＢ及び検出電圧ＶｍＢに基づいて、蓄電池１０の複素インピーダンスＺｍＢ（＝ＶｍＢ／ＩｍＢ）を算出する。その後、算出した複素インピーダンスＺｍＢに、第２電池状態算出装置２０Ｂにおける振幅影響係数ＫａｍｐＢ及び位相影響係数ＫｐｈＢを乗算することにより、「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」を算出する。

　振幅影響係数ＫａｍｐＢ及び位相影響係数ＫｐｈＢは、例えば、実験又は計算により予め算出された値であり、記憶部７２に記憶されている。制御装置７０は、「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」の算出に用いる振幅影響係数ＫａｍｐＢ及び位相影響係数ＫｐｈＢを、記憶部７２から読み出す。

　ステップＳ１８では、ステップＳ１７で算出した「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」から、ステップＳ１３で算出された「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」を差し引くことにより、補正値ΔＺ（＝ｊ×ω０×（ＭｖＡ－ＭｖＢ））を算出する。補正値ΔＺは、複素インピーダンスの虚部を補正するための値となる。補正値ΔＺは、角周波数ω０が高くなるほど、大きい値になる。算出した補正値ΔＺを、周波数ω０と紐付けて記憶部７２に記憶させる。

　なお、記憶部７２に記憶される情報は、周波数ω０及び補正値ΔＺが紐付けられたテーブル情報又は数式情報であってもよい。数式情報である場合、使用される数式は、例えば、多項式、指数関数、三角関数及びそれらの関数の組み合わせであってもよい。

　ちなみに、上式（ｅｑ１４）から下式（ｅｑＡ）が導かれる。

　図５は、角周波数ω０と補正値ΔＺとの関係を示す。インピーダンス計測における相互インダクタンスの誤差は、角周波数ω０に比例し、相互インダクタンスの差が比例定数になる。２点以上の周波数で測定をし、線形式を解く又はフィッティング等によって相互インダクタンスの差を算出することができる。角周波数ω０と補正値ΔＺとの関係式は、理論的には原点を通るため、原点を通る式を用いてフィッティングするのが好ましい。また、高周波側の測定値を用いた方がΔＺが大きくなるため、ΔＺの精度を向上できる。理想的には、測定周波数範囲において外挿にならないようにすることが望ましい。フィッティングを採用する場合、切片を固定せずにフィッティングする方法、又は切片を０にして（原点を通るようにして）フィッティングする方法がある。なお、測定する周波数は、複数に限らず１点であってもよい。この場合、原点を通るものとして比例定数（傾き）を算出することができる。

　ステップＳ１９では、算出した補正値ΔＺに基づいて、ステップＳ１７で算出した「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」を補正する。詳しくは、「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」から補正値ΔＺを差し引く。

　ステップＳ１８において記憶させた補正値ΔＺは、第２の組織が所有する第２電池状態算出装置２０Ｂにより複素インピーダンスの算出が実施される場合に用いられる。

　なお、第２電池状態算出装置２０Ｂにより算出された「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」を補正することに代えて、第１電池状態算出装置２０Ａにより算出された「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」を補正することもできる。詳しくは、ステップＳ１９において、第１電池状態算出装置２０Ａの制御装置７０は、ステップＳ１３で算出した「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」に、ステップＳ１８で算出された補正値ΔＺを加えることにより、「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」を補正する。

　また、第１，第２電池状態算出装置２０Ａ，２０Ｂにより算出された「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」及び「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」の双方を補正することもできる。詳しくは、ステップＳ１８において、ΔＺ／２を補正値として算出する。ステップＳ１９において、第１電池状態算出装置２０Ａの制御装置７０は、ステップＳ１３で算出した「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」に、算出された補正値ΔＺ／２を加えることにより、「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」を補正する。また、第２電池状態算出装置２０Ｂの制御装置７０は、ステップＳ１７で算出した「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」から補正値ΔＺ／２を差し引くことにより、「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」を補正する。

　ちなみに、ステップＳ１０～Ｓ１３における第１電池状態算出装置２０Ａの上記相対位置と、ステップＳ１４～Ｓ１７における第２電池状態算出装置２０Ｂの上記相対位置とは、同じであってもよいし、異なっていてもよい。

　図６～９を用いて、第１電池状態算出装置２０Ａにより算出された「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」の補正態様について説明する。

　図６は、ステップＳ１３，Ｓ１７において算出された複素インピーダンスＺｍＡ，ＺｍＢを示す。複素インピーダンスＺｍＡ，ＺｍＢは、基準信号の複数の角周波数ω０それぞれに紐付けられた複素インピーダンスの集合である。

　図７は、ステップＳ１３における「ＺｍＡ」及び「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」を示す。図８は、ステップＳ１７における「ＺｍＢ」及び「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」を示す。

　図９は、「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」と、補正された「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」とを示す。

　以上説明した本実施形態によれば、第１電池状態算出装置２０Ａにおいて算出された複素インピーダンスに対する、第２電池状態算出装置２０Ｂにおいて算出された複素インピーダンスの相対誤差を低減できる。これにより、各装置２０Ａ，２０Ｂにおいて算出された複素インピーダンスに含まれる誤差の差を低減することができる。

　上記メイン閉回路に交流電流が流れることに伴い発生した磁束が、電圧検出部３３及び電流検出部５２の検出値に及ぼす影響度合いは、電池状態算出装置２０の構成要素の相対位置が変わると変化する。逆に言えば、メイン閉回路に交流電流が流れることに伴い発生した磁束が、電圧検出部３３及び電流検出部５２の検出値に及ぼす影響の変化度合いは、上記相対位置が固定されていれば小さい。以下、この点について、図１０，１１を用いて説明する。図１０は、ファラデーの法則の概念図である。図１１は、Biot-Savartの法則の参考図である。

　下式（ｅｑＢ）にファラデーの法則を示す。下式（ｅｑＢ）において、「Ｅ（ｘ，ｔ）」は、電場ベクトルを示し、「Ｌ」は線積分の経路を示す。「Ｂ（ｘ，ｔ）」は磁束密度ベクトルを示す。「Ｓ」は左辺の線積分の経路によって囲まれる部分によって閉じる領域を示す。「ｎ」は「Ｓ」上の点の法線ベクトルを示す。「ｘ」は、電流素片からの位置を示すベクトルであり、「ｔ」は時間を示す。つまり、電場ベクトル「Ｅ（ｘ，ｔ）」と磁束密度ベクトル「Ｂ（ｘ，ｔ）」は、場所と時間に依存する値である。「Ｖｉ（ｔ）」は、誘導起電力を示す。

　任意の閉回路Ｃ０を電流Ｉが流れている場合、電流素片Ｉ×ｄｓが任意の点Ｐにつくる磁束密度Ｂは、任意の点Ｏを始点としたベクトルｓ及びｘを用いて、Biot-Savartの法則より、下式（ｅｑＣ）で表される。ここで、「ｔ」は時間を示し、「μ０」は真空の透磁率を示し、「μｒ（ｘ）」は比透磁率を示す。

　上式（ｅｑＣ）を上式（ｅｑＢ）に代入すると、下式（ｅｑＤ）が導かれる。

　上式（ｅｑＤ）の右辺の時間全微分ｄ／ｄｔは、各パラメータに対して偏微分を用いて表すことができる。例えば、時間ｔ及び位置ｘに依存するとすると、下式（ｅｑＥ）が導かれる。

　上式（ｅｑＥ）の右辺の第２項は、例えば、電圧検出閉回路３４（図２参照）の位置が時間によって変化する場合に発生する誘導起電力を表す。ｄｚ／ｄｔが周波数ｆを有する周期関数の場合、周波数ｆで変動する誘導起電力が発生する。同様に、偏微分に用いるパラメータを例えば電流が流れる閉回路（具体的には例えば、上記メイン閉回路）の位置ｓとすると、上式（ｅｑＥ）の右辺の第２項は、ｓが時間によって変化する場合の誘導起電力を表すこともでき、周期性を持つ場合は周期に応じた周波数での誘導起電力が発生する。他のパラメータに関しても同様である。

　ここでは、上式（ｅｑＥ）の右辺第１項に着目し、電流Ｉ（ｔ）のみが時間に依存するとすると、下式（ｅｑＦ）が導かれる。

　上式（ｅｑＦ）の下線部分は、電圧変動測定経路Ｌ（具体的には、電圧検出閉回路３４）と、電流測定経路Ｃ０（具体的には、電流検出閉回路５５）との間の相互インダクタンスＭである。上式（ｅｑＦ）から、電圧検出閉回路３４及び電流検出閉回路の形状が変わらない場合であって、透磁率μｒ（ｘ）が一定の場合に、相互インダクタンスＭは一定値となり変化しないことがわかる。補正値を算出するにあたり、電流検出閉回路５５と電圧検出閉回路３４とを固定することが肝要である。

　この点に鑑み、上記相対位置が固定された状態で、メイン閉回路に交流電流が流される。制御装置７０は、交流電流が流されている場合における電圧検出部３３の検出値及び電流検出部５２の検出値を取得し、取得した検出値に基づいて、複素インピーダンスを算出する。

　第１電池状態算出装置２０Ａで算出された複素インピーダンスに対する、第２電池状態算出装置で算出された複素インピーダンスの相対誤差の変化度合いは、上記相対位置が固定されていれば小さい。このため、本実施形態によれば、上記相対誤差を低減するための補正値の算出精度を高めることができる。

　＜第１実施形態の変形例＞
　・相互インダクタンスＭｉに対してシャント抵抗体４３の抵抗値Ｒｓが大きい場合（上式（ｅｑ８）参照）、インピーダンス補正処理で用いられる振幅影響係数ＫａｍｐＡ，ＫａｍｐＢを１にしてもよい。

　・インピーダンス測定用の電流のエネルギ源は、測定対象となる蓄電池１０に限らず、例えば、図１２に示すように商用電源１７であってもよい。電池状態算出装置２１は、商用電源１７をエネルギ源として、シャント抵抗体４３及び蓄電池１０を含むメイン閉回路に交流電流を流す交流電源５６（「交流駆動部」に相当）と、制御部７４とを備えている。なお、交流電源５６は、例えば、定電流制御又は定電圧制御のどちらで制御されてもよい。制御部７４は、図４と同様に、インピーダンス補正処理を行う。

　＜第２実施形態＞
　以下、第２実施形態について、第１実施形態との相違点を中心に図面を参照しつつ説明する。本実施形態では、図１３，１４を用いて、第１実施形態で説明した電池状態算出装置の具体例について説明する。図１３は、蓄電池１０と、蓄電池１０に取り付けられた電池状態算出装置１２０とを備えるシステムを示す。図１４は、図１３に示す構成のうち、電池ケース及び回路基板等の図示を省略した図である。

　電池状態算出装置１２０は、回路基板９０を備えている。回路基板９０には、制限抵抗体４１、電流制御スイッチ４２、シャント抵抗体４３、正極側電圧経路３１、負極側電圧経路３２、正極側電流経路５３及び負極側電流経路５４が設けられている。回路基板９０において、正極側電圧経路３１、負極側電圧経路３２、正極側電流経路５３及び負極側電流経路５４は、配線パターン及びビア等により構成されている。

　電池状態算出装置１２０は、処理部９１を備えている。処理部９１は、集積回路で構成され、回路基板９０に設けられている。処理部９１は、第１実施形態で説明した電圧検出部３３と、電流処理部５０と、基準信号生成部６０と、プロセッサ及び記憶部を含む制御装置７０との機能を有している。

　蓄電池１０は、捲回体１３、正極側集電体１４ａ、負極側集電体１４ｂ、正極側導電部材１５ａ、負極側導電部材１５ｂ及び電池ケース１６を備えている。捲回体１３は、扁平な直方体形状をなしている。捲回体１３の第１方向（幅方向）の第１端には、正極側集電体１４ａが設けられ、捲回体１３の幅方向の第２端には、負極側集電体１４ｂが設けられている。各集電体１４ａ，１４ｂは、捲回体１３の幅方向及び第２方向（厚さ方向）に直交する第３方向（高さ方向）に延びている。

　正極側集電体１４ａの高さ方向の一端には、正極側導電部材１５ａが設けられ、負極側集電体１４ｂの高さ方向の一端には、負極側導電部材１５ｂが設けられている。正極側導電部材１５ａには、正極端子１０ａが設けられ、負極側導電部材１５ｂには、負極端子１０ｂが設けられている。

　従来は、電池ケース１６の外側の配線に着目した複素インピーダンスの測定系が構築されていた。ここで、蓄電池１０の構成部品のうち核となる部品が捲回体１３である。捲回体１３から電流を引き出す蓄電池の内部構造体も回路の一部として考慮する必要がある。このため、メイン閉回路と電流検出閉回路５５との相互インダクタンスＭｖ、及びメイン閉回路と電流検出閉回路５５との相互インダクタンスＭｉを定量化する場合において、捲回体１３等の蓄電池１０の内部構造体に対する回路基板９０の相対位置も固定されていることが望ましい。そこで、本実施形態の電池状態算出装置１２０を構成する回路基板９０は、蓄電池１０に強固に固定されている。

　図１４に示す例において、蓄電池１０の内部構造を考慮した場合における電圧変動測定経路Ｌ及び電流測定経路Ｃ０について説明する。電圧変動測定経路Ｌ（具体的には、電圧検出閉回路３４）は、正極側電圧経路３１、正極側導電部材１５ａ、正極側集電体１４ａ、捲回体１３、負極側集電体１４ｂ、負極側導電部材１５ｂ、負極側電圧経路３２及び処理部９１を含む閉回路である。電流測定経路Ｃ０（具体的には、電流検出閉回路５５）は、シャント抵抗体４３、正極側電流経路５３、処理部９１及び負極側電流経路５４を含む閉回路である。

　相互インダクタンスＭｖは、上記メイン閉回路に流れる電流が、電圧変動測定経路Ｌにより形成される閉曲面上に発生させる磁場由来の相互インダクタンスである。相互インダクタンスＭｉは、メイン閉回路に流れる電流が、電流測定経路Ｃ０により形成される閉曲面上に発生させる磁場由来の相互インダクタンスである。図１４に示す例において、メイン閉回路、電流検出閉回路５５及び電圧検出閉回路３４それぞれの経路が一部重複しているが、磁場を発生させる電流と、その因果関係としての電圧変動の測定という点で差異があり、経路の重複部分も含めて相互インダクタンスとみなす。

　以上説明した本実施形態によれば、第１実施形態と同様の効果を奏することができる。

　＜第３実施形態＞
　以下、第３実施形態について、第１実施形態との相違点を中心に図面を参照しつつ説明する。本実施形態では、図１５に示すように、電池状態算出装置２２０が車両に備えられている。

　車両に搭載されるシステムは、電池状態算出装置２２０に加え、組電池１０００、インバータ１３０、回転電機１４０及びリレーＳＭＲを備えている。インバータ１３０は、上アームスイッチＳＨ及び下アームスイッチＳＬの直列接続体を複数相分（３相分）備えている。上アームスイッチＳＨ及び下アームスイッチＳＬの接続点には、回転電機１４０の電機子巻線が接続されている。回転電機１４０の回転子は、車両の駆動輪１５０と動力伝達可能になっている。インバータ１３０と組電池１０００とは、リレーＳＭＲを介して接続されている。

　組電池１０００は、複数の蓄電池１００の直列接続体を備えている。本実施形態では、便宜上、組電池１０００は、４つの蓄電池１００である第１，第２，第３，第４蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄを備えていることとする。各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄは、正極端子１０１ａ及び負極端子１０１ｂを備えている。各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄは、同じ仕様の蓄電池であり、具体的には、定格電圧及び満充電容量が同じである。

　システムは、インバータ制御装置１３１を備えている。インバータ制御装置１３１は、プロセッサ及び記憶部を有している。インバータ制御装置１３１は、車両を走行させるためにインバータ１３０の上，下アームスイッチＳＨ，ＳＬのスイッチング制御を行う。

　電池状態算出装置２２０において、制限抵抗体４１の第１端は、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄのうち最も高電位側の第１蓄電池１００Ａの正極端子１０１ａに接続されている。シャント抵抗体４３の第２端は、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄのうち最も低電位側の第４蓄電池１００Ｄの負極端子１０１ｂに接続されている。なお、本実施形態において、制限抵抗体４１、電流制御スイッチ４２、シャント抵抗体４３及び組電池１０００を含む閉回路がメイン閉回路に相当する。

　電池状態算出装置２２０は、監視ＩＣ９２を備えている。監視ＩＣ９２は、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄに対応して個別に設けられた電圧検出部３３、電流変調部５１、電流検出部５２、基準信号生成部６０及び制御部９３を備えている。電圧検出部３３と正極端子１０１ａとは正極側電圧経路３１により接続され、電圧検出部３３と負極端子１０１ｂとは負極側電圧経路３２により接続されている。

　制御部９３には、電圧検出部３３、電流検出部５２、温度センサ８０及び圧力センサ８１の検出値が入力される。温度センサ８０は、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの温度を検出する。圧力センサ８１は、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの圧力を検出する。

　制御部９３は、基準信号生成部６０に、基準信号の生成指令を入力する。また、「演算部」として機能する制御部９３は、電圧検出部３３及び電流検出部５２の検出値に基づいて、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスの実部及び虚部を算出する。制御部９３は、算出結果に基づいて、例えば、複素インピーダンス平面プロット（コールコールプロット）を作成し、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの蓄電状態（ＳＯＣ）や劣化状態（ＳＯＨ）を把握する。

　図１６に、電池状態算出装置２２０の一例を示す。各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄのうち、電位的に隣り合う蓄電池の正極端子１０１ａ及び負極端子１０１ｂは、中間導電部材１０３（例えばバスバー）により接続されている。第１蓄電池１００Ａの正極端子１０１ａと制限抵抗体４１とは、第１端部導電部材１０２（例えばバスバー）により電気的に接続されている。第４蓄電池１００Ｄの負極端子１０１ｂとシャント抵抗体４３とは、第２端部導電部材１０４（例えばバスバー）により電気的に接続されている。本実施形態において、電池状態算出装置２２０、インバータ１３０、回転電機１４０及び組電池１０００を含むシステムの構成部品の相対位置が固定されるようにシステムが構成されている。

　続いて、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスの相対誤差の差を低減するインピーダンス補正処理について説明する。本実施形態では、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄのうち、いずれか１つである基準電池の複素インピーダンスに対する、残りの対象電池の複素インピーダンスの相対誤差の差を低減させる。本実施形態において、基準電池は第１蓄電池１００Ａであり、対象電池は第２，第３，第４蓄電池１００Ｂ，１００Ｃ，１００Ｄである。

　本実施形態では、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの真の複素インピーダンスＺｔが同一であると仮定したインピーダンス補正処理が実行される。例えば、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄが新品の場合、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスが同一となる。以下、これまでに登場した数式において、第１蓄電池１００Ａに関する値に添え字Ａを付し、第２蓄電池１００Ｂに関する値に添え字Ｂを付し、第３蓄電池１００Ｃに関する値に添え字Ｃを付し、第４蓄電池１００Ｄに関する値に添え字Ｄを付す。

　上式（ｅｑ７）を参照すると、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄについて下式（ｅｑ１５）～（ｅｑ１８）が成立する。

　本実施形態では、共通の監視ＩＣ９２に電圧検出部３３、電流検出部５２及び基準信号生成部６０が備えられている。このため、「ΔβＡ，ΔβＢ，ΔβＣ，ΔβＤ＝０」及び「ΔγＡ，ΔγＢ，ΔγＣ，ΔγＤ＝０」であるとする。このため、下式（ｅｑ１９）～（ｅｑ２１）が導かれる。

　なお、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの電流は、共通の電流検出部５２で検出される。このため、上式（ｅｑ１９）～（ｅｑ２１）において、θｐＡ，θｐＢ，θｐＣ，θｐＤは同じ値になり、ＭｉＡ／ＲｓＡ，ＭｉＢ／ＲｓＢ，ＭｉＣ／ＲｓＣ，ＭｉＤ／ＲｓＤは同じ値になる。

　上式（ｅｑ１９）～（ｅｑ２１）において、振幅影響係数ＫａｍｐＡ，ＫａｍｐＢ，ＫａｍｐＣ，ＫａｍｐＤ及び位相影響係数ＫｐｈＡ，ＫｐｈＢ，ＫｐｈＣ，ＫｐｈＤは、電池状態算出装置２２０によって定まる係数である。

　上式（ｅｑ１９）において、ＺｍＡ，ＺｍＢは、第１，第２蓄電池１００Ａ，１００Ｂに対応する電圧検出部３３の検出値、及び電流検出部５２の検出値に基づいて算出される複素インピーダンスである。上式（ｅｑ２０）において、ＺｍＣは、第３蓄電池１００Ｃに対応する電圧検出部３３の検出値、及び電流検出部５２の検出値に基づいて算出される複素インピーダンスである。上式（ｅｑ２１）において、ＺｍＤは、第４蓄電池１００Ｄに対応する電圧検出部３３の検出値、及び電流検出部５２の検出値に基づいて算出される複素インピーダンスである。

　上式（ｅｑ１９）を参照すると、「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」と「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」との差が、第１，第２蓄電池１００Ａ，１００Ｂの複素インピーダンスの誤差となる。上式（ｅｑ２０）を参照すると、「ＺｍＣ×ＫａｍｐＣ×ＫｐｈＣ」と「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」との差が、第１，第３蓄電池１００Ａ，１００Ｃの複素インピーダンスの誤差となる。上式（ｅｑ２１）を参照すると、「ＺｍＤ×ＫａｍｐＤ×ＫｐｈＤ」と「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」との差が、第１，第４蓄電池１００Ａ，１００Ｄの複素インピーダンスの誤差となる。

　このため、「ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」、「ＺｍＣ×ＫａｍｐＣ×ＫｐｈＣ」及び「ＺｍＤ×ＫａｍｐＤ×ＫｐｈＤ」から、「ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」の影響を取り除くことにより、４つの蓄電池の複素インピーダンスの相対誤差を低減できる。

　図１７は、制御部９３により実行されるインピーダンス補正処理のフローチャートである。

　ステップＳ２０では、車両が停車中であるかを判定する。停車中であると判定したことを条件として、ステップＳ２１以降の処理の実行を許可する。これにより、車両の走行に伴い発生する振動が、電池状態算出装置２２０の構成部品の相対位置に及ぼす影響を抑制できる。その結果、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスの相対誤差を低減する補正値の算出精度を向上できる。

　ステップＳ２０において停車中であると判定した場合には、ステップＳ２１に進み、温度条件が成立しているか否かを判定する。温度条件は、温度センサ８０により検出された各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの温度が同一又は同等であるとの条件である。各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの温度が同等であるとは、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの温度のうち、最高温度と最低温度との差ΔＴの絶対値が、例えば、５℃以下、３℃以下、１℃以下又は０．５℃以下のことである。

　温度条件が成立していると判定した場合には、ステップＳ２２に進み、電圧条件が成立しているか否かを判定する。電圧条件は、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄに電流が流れていない状態において、電圧検出部３３により検出された各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの端子間電圧が同一又は同等であるとの条件である。各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの端子間電圧が同等であるとは、例えば、「｜ΔＶ｜／Ｖｓｔｄ×１００≦０．３％」、「｜ΔＶ｜／Ｖｓｔｄ×１００≦０．１％」、又は「｜ΔＶ｜／Ｖｓｔｄ×１００≦０．０３％」のことである。ΔＶは、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの端子間電圧のうち、最高電圧と最低電圧との差ΔＶの絶対値である。Ｖｓｔｄは、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの定格電圧である。

　なお、電圧条件において、検出電圧に代えて、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００ＤのＳＯＣが用いられてもよい。また、温度，電圧条件に加え、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの圧力条件が追加されてもよい。圧力条件は、圧力センサ８１により検出された各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの圧力が同一又は同等であるとの条件である。

　ステップＳ２２において電圧条件が成立していると判定した場合には、ステップＳ２３に進み、電流変調部５１によって電流制御スイッチ４２が制御されるように、基準信号生成部６０に基準信号の生成指令を入力する。ここでは、基準信号の角周波数ω０を徐々に変化（漸増又は漸減）させる周波数スイープを行う。これにより、基準信号の周波数ω０で変動する交流電流が各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄに流れる。交流電流が流れている状態で、電流検出部５２の検出電流Ｉｍを取得し、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄに対応する電圧検出部３３の検出電圧ＶｍＡ，ＶｍＢ，ＶｍＣ，ＶｍＤを取得する。詳しくは、角周波数に対応する検出電流Ｉｍ及び検出電圧ＶｍＡ，ＶｍＢ，ＶｍＣ，ＶｍＤを取得する。

　ステップＳ２４では、取得した検出電流Ｉｍ及び検出電圧ＶｍＡ，ＶｍＢ，ＶｍＣ，ＶｍＤに基づいて、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスＺｍＡ（＝ＶｍＡ／Ｉｍ），ＺｍＢ（＝ＶｍＢ／Ｉｍ），ＺｍＣ（＝ＶｍＣ／Ｉｍ），ＺｍＤ（＝ＶｍＤ／Ｉｍ）を算出する。

　その後、算出した第１蓄電池１００Ａの複素インピーダンスＺｍＡに振幅影響係数ＫａｍｐＡ及び位相影響係数ＫｐｈＡを乗算することにより、「ＺＡ＝ＺｍＡ×ＫａｍｐＡ×ＫｐｈＡ」を算出する。また、算出した第２蓄電池１００Ｂの複素インピーダンスＺｍＢに振幅影響係数ＫａｍｐＢ及び位相影響係数ＫｐｈＢを乗算することにより、「ＺＢ＝ＺｍＢ×ＫａｍｐＢ×ＫｐｈＢ」を算出する。また、算出した第３蓄電池１００Ｃの複素インピーダンスＺｍＣに振幅影響係数ＫａｍｐＣ及び位相影響係数ＫｐｈＣを乗算することにより、「ＺＣ＝ＺｍＣ×ＫａｍｐＣ×ＫｐｈＣ」を算出する。また、算出した第４蓄電池１００Ｄの複素インピーダンスＺｍＤに振幅影響係数ＫａｍｐＤ及び位相影響係数ＫｐｈＤを乗算することにより、「ＺＤ＝ＺｍＤ×ＫａｍｐＤ×ＫｐｈＤ」を算出する。

　振幅影響係数ＫａｍｐＡ，ＫａｍｐＢ，ＫａｍｐＣ，ＫａｍｐＤ及び位相影響係数ＫｐｈＡ，ＫｐｈＢ，ＫｐｈＣ，ＫｐｈＤは、例えば、実験又は計算により予め算出された値であり、制御部９３の記憶部に記憶されている。制御部９３は、ＺＡ～ＺＤの算出に用いる振幅影響係数及び位相影響係数記憶部から読み出す。

　ステップＳ２５では、下式（ｅｑ２２）に示すように、ステップＳ２４で算出したＺＢからＺＡを差し引くことにより、補正値ΔＺｋＢを算出する。

　ステップＳ２５では、下式（ｅｑ２３）に示すように、ステップＳ２４で算出したＺＣからＺＡを差し引くことにより、補正値ΔＺｋＣを算出する。

　ステップＳ２５では、下式（ｅｑ２４）に示すように、ステップＳ２４で算出したＺＤからＺＡを差し引くことにより、補正値ΔＺｋＤを算出する。

　補正値ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤは、複素インピーダンスの虚部を補正するための値となる。補正値ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤは、角周波数ω０が高くなるほど、大きい値になる。

　ステップＳ２６では、算出した補正値ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤに基づいて、ステップＳ２４で算出したＺＢ，ＺＣ，ＺＤを補正する。詳しくは、下式（ｅｑ２５）に示すように、ＺＢから補正値ΔＺｋＢを差し引くことにより、第２蓄電池１００Ｂの補正インピーダンスＺｃＢを算出する。

　また、下式（ｅｑ２６）に示すように、ＺＣから補正値ΔＺｋＣを差し引くことにより、第３蓄電池１００Ｃの補正インピーダンスＺｃＣを算出する。

　また、下式（ｅｑ２７）に示すように、ＺＤから補正値ΔＺｋＤを差し引くことにより、第４蓄電池１００Ｄの補正インピーダンスＺｃＤを算出する。

　ステップＳ２７では、算出した補正値ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤを、角周波数ω０と紐付けて記憶部に記憶させる。なお、記憶部に記憶される情報は、周波数ω０及び補正値ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤが紐付けられたテーブル情報又は数式情報であってもよい。数式情報である場合、使用される数式は、例えば、多項式、指数関数、三角関数及びそれらの関数の組み合わせであってもよい。

　以上説明した本実施形態によれば、第１蓄電池１００Ａの複素インピーダンスに対する第２，第３，第４蓄電池１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスの相対誤差を低減でき、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスに含まれる誤差の差を低減できる。

　＜第３実施形態の変形例＞
　・制御部９３は、図１７のステップＳ２１～Ｓ２６に示したインピーダンス補正処理を、記憶部に記憶されていた補正値の情報が消失したと判定した場合に実行してもよい。

　・制御部９３は、ステップＳ２１～Ｓ２６に示したインピーダンス補正処理を、車両の販売店に車両が保管されている場合に実行してもよい。

　・制御部９３は、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの使用期間が所定期間以下であると判定したことを条件に、インピーダンス補正処理を行ってもよい。各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの使用期間が過度に長くなければ、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄが中古状態になっていたとしても、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの真の複素インピーダンスＺｔのばらつきは小さいと考えられる。

　・図１７のステップＳ２１～Ｓ２６に示したインピーダンス補正処理は、組電池１０００が車両に搭載される前段階において実行されてもよい。例えば、工場における組電池１０００の製造工程において実行されてもよい。

　・図１８，１９に示すように、複数の蓄電池１００Ａ～１００Ｄの直列接続体及び監視ＩＣ９２が一体化されたモジュール３００がシステムに備えられていてもよい。この場合、電圧検出部３３、電流処理部５０及び基準信号生成部６０が所定距離離間して配置されるため、位相差Δα，Δβ，Δγの位相差が０ではなくなることがある。

　・図２０に示すように、シャント抵抗体１４３が第４蓄電池１００ＤとリレーＳＭＲとの間に設けられていてもよい。この場合、メイン閉回路は、破線矢印にて例示するように、各蓄電池１００Ａ～１００Ｄ、リレーＳＭＲ、上アームスイッチＳＨ、回転電機１４０の電機子巻線、下アームスイッチＳＬ、リレーＳＭＲ及びシャント抵抗体４３を含む閉回路となる。また、インピーダンス補正処理において、各１００Ａ～１００Ｄに交流電流を流す制御は、インバータ制御装置１３１によるインバータ１３０のスイッチング制御となる。なお、この場合、インバータ１３０及びインバータ制御装置１３１が「交流駆動部」に相当する。

　・図２１に示すように、電圧検出用のＩＣ（監視ＩＣ９２）とは別に、電流検出用ＩＣ９４が電池状態算出装置２２０に備えられていてもよい。

　車両に搭載されるシステムは、切替スイッチＳＭＣと、ＤＣＤＣコンバータを含む充電器１３２とを備えている。充電器１３２は、切替スイッチＳＭＣを介して蓄電池１００Ａ～１００Ｄの直列接続体に接続されている。なお、本実施形態では、便宜上、充電器１３２及び切替スイッチＳＭＣがインバータ制御装置１３１により制御されるものとする。

　インバータ制御装置１３１は、回転電機１４０のロータを回転させて車両を走行させるべく、インバータ１３０の上，下アームスイッチＳＨ，ＳＬのスイッチング制御を行う。この場合、インバータ制御装置１３１は、リレーＳＭＲをオンし、切替スイッチＳＭＣをオフする第１状態とする。一方、インバータ制御装置１３１は、外部電源１３３から充電器１３２を介して蓄電池１００Ａ～１００Ｄの直列接続体を充電する外部充電を行うべく、充電器１３２を制御する。この場合、インバータ制御装置１３１は、切替スイッチＳＭＣをオンし、リレーＳＭＲをオフする第２状態とする。

　制御部９３は、インバータ制御装置１３１と協働して、第１状態及び第２状態それぞれにおいて、図１７に示したインピーダンス補正処理を行う。これにより、第１状態及び第２状態それぞれに対応した補正値を算出できる。このようにリレーＳＭＲ及び切替スイッチＳＭＣを制御することにより、第１状態では外部電源内のインピーダンスの影響を受けることなく補正値を正確に算出することができる。また、第２状態ではインバータ１３０内の平滑コンデンサ等によるインピーダンスの影響を受けることなく補正値を算出できる。

　・直列接続された複数の蓄電池のうち、最高電位の蓄電池の正極端子と最低電位の蓄電池の負極端子とを接続する電気経路４５の配置態様の一例について説明する。電気経路４５は、制限抵抗体４１、電流制御スイッチ４２及びシャント抵抗体４３を含む。

　図２２には、６つの蓄電池１００Ａ～１００Ｆの例が示されている。この場合における各蓄電池１００Ａ～１００Ｆの相互インダクタンスＭｖの計算結果を図２３に示す。図２４には、図２２とは異なる電気経路４５の配置態様を示す。この場合における各蓄電池１００Ａ～１００Ｆの相互インダクタンスＭｖの計算結果を図２５に示す。各蓄電池１００Ａ～１００Ｆ毎に相互インダクタンスＭｖが異なるため、各蓄電池１００Ａ～１００Ｆのうち対象電池に対応した補正値を算出することが望ましい。

　なお、複数の蓄電池及び電池状態算出装置の相対位置によっては、図２６に示す計算結果になることもある。図２６では、直列接続された複数の蓄電池のうち、最高，最低電池の蓄電池１００Ａ，１００Ｅの相互インダクタンスＭｖが同じ値又は同等の値となり、最高，最低電池の蓄電池１００Ａ，１００Ｅに隣接する蓄電池１００Ｂ，１００Ｄの相互インダクタンスＭｖが同じ値又は同等の値になっている。この場合、複数の蓄電池１００Ａ～１００Ｅを構成する４つの対象電池のうち、一部の対象電池について算出された補正値を、他の対象電池の複素インピーダンスの補正に用いることができる。このため、電池状態算出装置における演算負荷が軽減される。

　＜第４実施形態＞
　以下、第４実施形態について、第３実施形態との相違点を中心に図面を参照しつつ説明する。本実施形態では、研究機関の実験室等、車両外部の電池状態算出装置（以下、外部算出装置）において算出された基準電池の複素インピーダンスに対する、車載蓄電池の複素インピーダンスの相対誤差を低減させる。外部算出装置は、例えば第１実施形態の第１電池状態算出装置２０Ａであり、「第１電池状態算出装置」に相当する。

　車両に備えられた組電池１０００と同一仕様の組電池の複素インピーダンスが、外部算出装置において算出される。この場合において、組電池を構成する４つの蓄電池のうち１つの基準電池（例えば、第１蓄電池）の複素インピーダンスＺｍＳが、外部算出装置において電流検出部５２及び電圧検出部３３の検出値に基づいて算出されたとする。基準電池の真のインピーダンスＺｓは下式（ｅｑ２８）で表される。

　一方、車両に備えられた電池状態算出装置２２０（「第２電池状態算出装置」に相当）において算出された第１，第２，第３，第４蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスと真のインピーダンスとの関係は、下式（ｅｑ２９）～（ｅｑ３２）で表される。

　車両に備えられた第１，第２，第３，第４蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの真のインピーダンスＺＡ，ＺＢ，ＺＣ，ＺＤと、外部算出装置の測定対象となる基準電池（第１蓄電池）の真のインピーダンスＺｓとが同一であると仮定する。この場合、上式（ｅｑ２８）～（ｅｑ３２）を踏まえ、第１，第２，第３，第４蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの補正値ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤが下式（ｅｑ３３）～（ｅｑ３６）のように定まる。

　補正値ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤは、例えば、電池状態算出装置２２の設計時において、外部算出装置における複素インピーダンスの算出結果等に基づいて予め算出される。算出された補正値ＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤは、制御部９３の記憶部に記憶される。

　図２７は、制御部９３により実行されるインピーダンス補正処理のフローチャートである。

　ステップＳ３０では、ステップＳ２０と同様に、車両が停車中であるかを判定する。

　ステップＳ３０において停車中であると判定した場合には、ステップＳ３１に進み、ステップＳ２１と同様に、温度条件が成立しているか否かを判定する。

　温度条件が成立していると判定した場合には、ステップＳ３２に進み、ステップＳ２２と同様に、電圧条件が成立しているか否かを判定する。

　ステップＳ３２において電圧条件が成立していると判定した場合には、ステップＳ３３に進み、ステップＳ２３と同様に、電流変調部５１によって電流制御スイッチ４２が制御されるように、基準信号生成部６０に基準信号の生成指令を入力する。交流電流が各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄに流れている状態で、電流検出部５２の検出電流Ｉｍを取得し、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄに対応する電圧検出部３３の検出電圧ＶｍＡ，ＶｍＢ，ＶｍＣ，ＶｍＤを取得する。詳しくは、角周波数に対応する検出電流Ｉｍ及び検出電圧ＶｍＡ，ＶｍＢ，ＶｍＣ，ＶｍＤを取得する。

　ステップＳ３４では、ステップＳ２４と同様に、取得した検出電流Ｉｍ及び検出電圧ＶｍＡ，ＶｍＢ，ＶｍＣ，ＶｍＤに基づいて、各蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの複素インピーダンスＺｍＡ（＝ＶｍＡ／Ｉｍ），ＺｍＢ（＝ＶｍＢ／Ｉｍ），ＺｍＣ（＝ＶｍＣ／Ｉｍ），ＺｍＤ（＝ＶｍＤ／Ｉｍ）を算出する。

　振幅影響係数ＫａｍｐＡ，ＫａｍｐＢ，ＫａｍｐＣ，ＫａｍｐＤ及び位相影響係数ＫｐｈＡ，ＫｐｈＢ，ＫｐｈＣ，ＫｐｈＤは、例えば、実験又は計算により予め算出された値であり、制御部９３の記憶部に記憶されている。制御部９３は、ＺＡ～ＺＤの算出に用いる振幅影響係数及び位相影響係数を記憶部から読み出す。

　ステップＳ３５では、上式（ｅｑ３３）～（ｅｑ３６）の補正値ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤを記憶部から読み出す。補正値ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤは、複素インピーダンスの虚部を補正するための値となる。補正値ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤは、角周波数ω０が高くなるほど、大きい値になる。

　ステップＳ３６では、読みだした補正値ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤに基づいて、ステップＳ３４で算出したＺＡ，ＺＢ，ＺＣ，ＺＤを補正する。詳しくは、ＺＡ，ＺＢ，ＺＣ，ＺＤから補正値ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤを差し引くことにより、第１，第２，第３，第４蓄電池１００Ａ，１００Ｂ，１００Ｃ，１００Ｄの補正インピーダンスＺｃＡ，ＺｃＢ，ＺｃＣ，ＺｃＤを算出する。

　ステップＳ３７では、算出した補正インピーダンスＺｃＡ，ＺｃＢ，ＺｃＣ，ＺｃＤを、角周波数ω０と紐付けて記憶部に記憶させる。

　以上説明した本実施形態によれば、車両外部において算出された基準電池の複素インピーダンスを基準として、車載蓄電池の複素インピーダンスの相対誤差のばらつきを低減できる。

　＜その他の実施形態＞
　なお、上記各実施形態は、以下のように変更して実施してもよい。

　・第３，第４実施形態において、電圧検出閉回路３４及び電流検出閉回路５５に、インバータ１３０に流れる電流等に起因する磁場が鎖交する場合に誘導起電力が発生し、ノイズとなり得る。インバータ１３０及び各閉回路３４，５５を含むシステムの相対位置を固定する固定機構がシステムに備えられている。ここで、基準信号の角周波数ω０は、システムの固有振動数から所定周波数（例えば数Ｈｚ）ずれた周波数に設定されればよい。

　固有振動数が複数存在する場合、複数の固有振動数のうち１つの固有振動数に対してずらすのではなく、複数の固有振動数に対してずらした値を基準信号の角周波数ω０として選択してもよい。固有振動数に対する角周波数ω０のずらし量は、例えば、インピーダンス計測回路が有するバンドパスフィルタを通過しない程度の値にしてもよい。例えば、バンドパスフィルタ幅が１０Ｈｚの場合、上記ずらし量は、１０Ｈｚ以上とするのが好適である。蓄電池のインピーダンス値は周波数に対して滑らかに変化するのが一般的である。本来用いたい角周波数ω０が固有振動数と一致することを回避するために、測定周波数（角周波数ω０）をずらしたとしても、蓄電池の制御に与える影響が少ないため、測定周波数をずらす方式が成立する。測定周波数をずらすのが困難な場合は、システムの固有振動数をずらしてもよい。例えば、固定機構が抑える箇所の間隔を変えたり、部品の形状や重さを変えたりすることにより、システムの固有振動数をずらすことができる。補正値は、振動が発生していない場合や、固有振動を避けた角周波数ω０を選択する場合に算出するのが好適である。

　電池容量の増加に伴い、蓄電池の複素インピーダンスの虚部の値が小さくなる。電池容量が大きい蓄電池（例えば、１００Ａｈ級の蓄電池）の場合、複素インピーダンスの虚部が例えば１０μΩオーダーとなる。一方、固有振動に伴う虚部の変化量は、例えば約２．５μΩである。この場合、上記変化量を無視できない。そこで、振動数をずらす対策が有効である。

　・例えば第３，第４実施形態において、電池状態算出装置は、各電圧検出回路の正極側電圧経路３１及び負極側電圧経路３２の長さＬと、正極側電圧経路３１及び負極側電圧経路３２の断面積Ａとに基づいて、補正対象となる各蓄電池に対する個別の補正値を算出してもよい。

　上式（ｅｑ１５）を例にすると、長さＬ及び断面積Ａの影響を考慮すると、下式（ｅｑ３７）に示すように、係数Ｇｔが追加される。

　正極側電圧経路３１及び負極側電圧経路３２の長さＬが長いほど、正極側電圧経路３１及び負極側電圧経路３２の抵抗値が大きくなる。その結果、正極側電圧経路３１及び負極側電圧経路３２で発生する電圧降下量も大きくなり、検出電圧と実際の電圧との差が大きくなる。この点を踏まえ、図２３に示すように、電池状態算出装置は、長さＬが長いほど、係数Ｇｔを小さくすればよい。

　正極側電圧経路３１及び負極側電圧経路３２の断面積が小さいほど、正極側電圧経路３１及び負極側電圧経路３２の抵抗値は大きくなる。この点を踏まえ、図２８に示すように、電池状態算出装置は、断面積Ａが小さいほど、係数Ｇｔを小さくすればよい。

　基準信号の角周波数ω０の影響を考慮すると、上式（ｅｑ３７）に示すように、係数Ｓｔが追加される。複素インピーダンスの算出時において、各周波数ω０の周波数スイープが行われる。角周波数ω０の高低に応じて、表皮効果の影響が異なる。高周波ほど抵抗値が大きくなる。この点を踏まえ、電池状態算出装置は、図２９に示すように、角周波数ω０が高いほど、係数Ｓｔを小さくすればよい。なお、係数Ｇｔ，Ｓｔは、例えば、予め計算により算出できる。

　表皮効果による配線の抵抗増加は、配線となる金属体の表面積が大きくなれば低減することができる。そのため、ハーネスとして多芯線を用いることにより、抵抗増加を低減できる。しかし、近年、バッテリーモジュールなどに配線用途として用いられるＦＰＣ基板では、各配線が単心線と同等の構成となっている。そのため、特にモジュールやパックといった配線距離が長くなる対象に対しては、表皮効果の補正を行うことが好適である。

　また、ＦＰＣ配線は単線相当のパターン配線となるが、表面積を大きくするために、例えば同電位の配線を２本以上用いたり、網目構造の配線パターンを用いたりする方法を用いてもよい。

　また、以下の（１），（２）の構成を採用してもよい。

　（１）測定対象となる組電池を構成する電池セルの場所と、電池セルの電流，電圧を測定する測定ＩＣとの距離に応じて、遠い場所ほど同電位の配線数を多くし、近い場所は単線又は遠い場所より少ない本数にするなどの構成にする。

　（２）測定ＩＣからの距離が違い場所は単線で構成し、遠い場所は複数線で構成する。近い場所では、多数の電池セルの線を配置する必要があるので場所が少ないが、遠い場所では、まだ配線していない電池セルの数が少なくなるので複数線とする場所を確保しやすい。このため、ＦＰＣ基板の面積増加を抑制しつつ、遠い場所の電池セルに対して表皮効果の影響を低減できるという効果がある。

　・電圧検出部３３及び電流検出部５２は、図３０に示すように、取得した電圧差が入力されるフィルタ回路４００と、フィルタ回路４００の出力値が入力されるＡＤコンバータ４０１とを備えている。ＡＤコンバータ４０１の出力値が、電圧，電流検出値となる。

　ここで、フィルタ回路４００の時定数に応じて、フィルタ回路４００のゲイン周波数特性及び位相周波数特性が変化する。ここで、電圧検出部３３が備えるフィルタ回路４００におけるゲイン周波数特性の変化量をＧｖ（ω０）とし、電流検出部５２が備えるフィルタ回路４００におけるゲイン周波数特性の変化量をＧｉ（ω０）とする。この場合、上式（ｅｑ１５）を例にすると、変化量の影響を考慮して下式（ｅｑ３８）が導かれる。

　この点を踏まえ、例えば第３，第４実施形態において、電池状態算出装置は、フィルタ回路４００の周波数特性に基づいて、補正対象となる各蓄電池に対する個別の補正値を算出してもよい。なお、係数Ｇｉ，Ｇｖは、例えば、予め計算により算出できる。

　・第３，第４実施形態において、車両に備えられた制御部に代えて、車両外部のサーバ（例えばクラウドサーバ）において、インピーダンス補正処理が実行されてもよい。この場合、サーバは、多数の車両に対応する補正値を収集することができ、収集した補正値をサーバの記憶部に記憶させる。サーバは、例えば、収集した補正値の統計値（例えば平均値）を補正値として、各車両の制御部に送信してもよい。

　・図３１に示すように、蓄電池の正極端子及び負極端子を接続するとともに、電流制御スイッチ４２及びシャント抵抗体４３を含む電気経路において、シャント抵抗体４３の電流流通方向とは逆方向に電流が流れる第１折り返し部４６ａと、第１折り返し部４６ａの電流流通方向とは逆方向に電流が流れる第２折り返し部４６ｂとが設けられてもよい。

　第１折り返し部４６ａに流れる電流によって発生する磁場は、シャント抵抗体４３に流れる電流によって電流検出閉回路５５を鎖交する磁場を低減させる役割を果たす。第２折り返し部４６ｂは、第１折り返し部４６ａによって変更された電流経路を元の経路に極力戻すためのものである。第２折り返し部４６ｂは、電流検出閉回路５５に対して第１折り返し部４６ａよりも遠い位置に配置される。このため、第２折り返し部４６ｂに流れる電流によって電流検出閉回路５５を鎖交する磁場は小さく、鎖交磁場が電流検出誤差に及ぼす影響は小さい。

　・第２実施形態において、図３２に示すように、回路基板９０において、シャント抵抗体４３が備えられる第１板面の裏面である第２板面に、図３１に示した第１折り返し線４６に相当する第１折り返し線４７が備えられてもよい。シャント抵抗体４３と第１折り返し線４７との間に電子部品が介在しないため、シャント抵抗体４３と第１折り返し線４７との電位差を小さくでき、絶縁破壊が発生しにくくなる。なお、第１折り返し部４７は、蓄電池１０の負極端子１０ｂ側の電気経路ではなく、蓄電池１０の正極端子１０ａ側の電気経路に設けられてもよい。

　・車両に備えられる蓄電池は、種々のものを採用することができる。例えば、蓄電池として、長板状のブレードセルが採用されてもよい。

　また、蓄電池としては、複数の単位電池（例えば電池セル）を電池ブロックにまとめた上で、各電池ブロックを直列接続する構成であってもよい。

　また、蓄電池としては、電池ブロックを作成することなく、複数の電池セルの直列接続体が車両のシャーシの収容部に収容されるいわゆるＣＴＰ（Cell to Pack）の構成であってもよい。

　また、ＣＴＰの構成に代えて、車両のシャーシに電池セルを収容する収容部が構成され、収容部に複数の電池セルが収容されるいわゆるＣＴＣ（Cell to Chassis）の構成が用いられてもよい。

　・蓄電池としては、２次電池に限らず、例えば燃料電池であってもよい。

　・電池状態算出装置が搭載される移動体としては、車両に限らず、例えば、航空機又は船舶であってもよい。また、電池状態算出装置の搭載先は、移動体に限らず、定置式の装置であってもよい。

　・本開示は、車両等の移動体に搭載される交換式蓄電池に対しても適用できる。この場合、交換式蓄電池を移動体に搭載した状態でインピーダンス計測を行う場合と、交換式蓄電池を移動体から降ろした状態でインピーダンス計測を行う場合とで、測定系が異なる場合がある。例えば、移動体から降ろした場合に移動体には搭載していない電源を信号源として用いて蓄電池のインピーダンス計測を行う場合などがある。しかし、測定用の信号源、電流検出部及び電圧検出部を備える構成自体は、移動体に搭載した状態と、移動体から降ろした状態とで同じになるため、本開示の補正方法が使用できる。この補正方法を用いることにより、移動体に搭載した場合と、移動体から降ろした場合とのどちらでも比較可能となる。

　・本開示に記載の制御部及びその手法は、コンピュータプログラムにより具体化された一つ乃至は複数の機能を実行するようにプログラムされたプロセッサ及びメモリを構成することによって提供された専用コンピュータにより、実現されてもよい。あるいは、本開示に記載の制御部及びその手法は、一つ以上の専用ハードウェア論理回路によってプロセッサを構成することによって提供された専用コンピュータにより、実現されてもよい。もしくは、本開示に記載の制御部及びその手法は、一つ乃至は複数の機能を実行するようにプログラムされたプロセッサ及びメモリと一つ以上のハードウェア論理回路によって構成されたプロセッサとの組み合わせにより構成された一つ以上の専用コンピュータにより、実現されてもよい。また、コンピュータプログラムは、コンピュータにより実行されるインストラクションとして、コンピュータ読み取り可能な非遷移有形記録媒体に記憶されていてもよい。

　本開示は、実施例に準拠して記述されたが、本開示は当該実施例や構造に限定されるものではないと理解される。本開示は、様々な変形例や均等範囲内の変形をも包含する。加えて、様々な組み合わせや形態、さらには、それらに一要素のみ、それ以上、あるいはそれ以下、を含む他の組み合わせや形態をも、本開示の範疇や思想範囲に入るものである。

Claims

　蓄電池（１０，１０００）の正極端子と前記蓄電池の負極端子とを電気的に接続する電気経路（４１～４３，ＳＭＲ，１３０，１４０）を備えるシステムに適用される電池状態算出装置（２０，１２０，２２０）において、
　前記電気経路及び前記蓄電池を含む閉回路に交流電流を流す交流駆動部（４２，５１，５６，６０，１３０，１３１）と、
　前記蓄電池に並列接続され、前記蓄電池の端子間電圧を検出する電圧検出回路（３１～３３）と、
　前記閉回路に流れる電流を検出する電流検出回路（４３，５２～５４）と、
　演算部（７０，９３）と、
を備え、
　前記演算部は、
　前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得し、取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記蓄電池の複素インピーダンスを算出し、
　算出した前記蓄電池の複素インピーダンスと比較する基準電池の複素インピーダンスと、算出した前記蓄電池の複素インピーダンスとに基づいて、前記基準電池の複素インピーダンスに対する、算出した前記蓄電池の複素インピーダンスの相対誤差を補正する補正値（ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤ）を算出し、
　算出した前記補正値に基づいて、算出した前記蓄電池の複素インピーダンスを補正する、電池状態算出装置。
　前記蓄電池は、複数の単位電池（１００Ａ～１００Ｄ）の直列接続体であり、
　前記電圧検出回路は、前記各単位電池に対応して個別に設けられ、
　前記各単位電池のうち、いずれか１つを前記基準電池（１００Ａ）とし、残りを対象電池（１００Ｂ～１００Ｄ）とする場合において、
　前記演算部（９３）は、
　前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記各電圧検出回路の検出値及び前記電流検出回路の検出値を取得し、取得した前記各電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記各単位電池の複素インピーダンスを算出し、
　算出した前記各単位電池の複素インピーダンスに基づいて、算出した前記基準電池の複素インピーダンスに対する、算出した前記対象電池の複素インピーダンスの相対誤差を補正する前記補正値（ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤ）を算出し、
　算出した前記補正値に基づいて、算出した前記対象電池の複素インピーダンスを補正する、請求項１に記載の電池状態算出装置。
　前記各単位電池、前記電気経路、前記各電圧検出回路及び前記電流検出回路の相対位置が固定されるように前記電池状態算出装置が構成されており、
　前記演算部は、前記相対位置が固定された状態で前記閉回路に交流電流が流されている場合における前記各電圧検出回路の検出値及び前記電流検出回路の検出値を取得する、請求項２に記載の電池状態算出装置。
　前記演算部は、前記基準電池の端子間電圧を検出する前記電圧検出回路及び前記閉回路の相互インダクタンスと、前記対象電池の端子間電圧を検出する前記電圧検出回路及び前記閉回路の相互インダクタンスとの差に応じた前記補正値を算出する、請求項３に記載の電池状態算出装置。
　前記演算部は、前記交流駆動部によって流される交流電流の周波数（ω０）が高いほど、前記補正値を大きく算出する、請求項４に記載の電池状態算出装置。
　前記蓄電池は、３つ以上の前記単位電池の直列接続体であり、
　前記演算部は、前記各対象電池に対する個別の前記補正値を算出する、請求項２～５のいずれか１項に記載の電池状態算出装置。
　前記システムは、車両に搭載され、
　前記演算部は、前記車両が停車中であると判定したことを条件に、前記閉回路に交流電流が流されている場合における前記各電圧検出回路の検出値及び前記電流検出回路の検出値を取得する、請求項２～５のいずれか１項に記載の電池状態算出装置。
　前記各電圧検出回路は、
　前記単位電池の正極端子に電気的に接続される正極側電圧経路（３１）と、
　前記単位電池の負極端子に電気的に接続される負極側電圧経路（３２）と、
を有し、前記正極側電圧経路及び前記負極側電圧経路の電圧差を前記単位電池の端子間電圧として検出し、
　前記演算部は、前記各電圧検出回路の前記正極側電圧経路及び前記負極側電圧経路の長さに基づいて、前記各対象電池に対する個別の前記補正値を算出する、請求項６に記載の電池状態算出装置。
　前記各電圧検出回路は、
　前記単位電池の正極端子に電気的に接続される正極側電圧経路（３１）と、
　前記単位電池の負極端子に電気的に接続される負極側電圧経路（３２）と、
　前記正極側電圧経路及び前記負極側電圧経路の電圧差が入力されるフィルタ回路（４００）と、
を有し、前記フィルタ回路の出力値に基づいて前記単位電池の端子間電圧を検出し、
　前記電流検出回路は、
　前記電気経路に設けられたシャント抵抗体（４３）と、
　前記シャント抵抗体の両端に接続された電流経路（５３，５４）と、
　前記両端に接続された電流経路の電圧差が入力されるフィルタ回路（４００）と、
を有し、前記フィルタ回路の出力値に基づいて前記閉回路に流れる電流を検出し、
　前記演算部は、前記電圧検出回路及び前記電流検出回路それぞれの前記フィルタ回路の周波数特性に基づいて、前記各対象電池に対する個別の前記補正値を算出する、請求項６に記載の電池状態算出装置。
　前記電池状態算出装置は、第１電池状態算出装置（２０Ａ）と、前記第１電池状態算出装置とは異なる場所に設けられる第２電池状態算出装置（２２０）とであり、
　前記第１電池状態算出装置は、
　前記蓄電池の正極端子が接続される正極接続部（４０ａ）と、
　前記蓄電池の負極端子が接続される負極接続部（４０ｂ）と、
　前記正極接続部と前記負極接続部とを接続する前記電気経路と、
を備え、
　前記第１電池状態算出装置において、前記演算部は、前記正極接続部及び前記負極接続部に、前記蓄電池と同一仕様の前記基準電池が接続されるとともに、前記基準電池、前記正極接続部、前記負極接続部、前記電気経路、前記電圧検出回路及び前記電流検出回路の相対位置が固定された状態において、前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得し、取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記基準電池の複素インピーダンスを算出し、
　前記第２電池状態算出装置において、前記演算部は、前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得し、取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記蓄電池の複素インピーダンスを算出し、
　算出した前記蓄電池の複素インピーダンスと、前記第１電池状態算出装置において算出した前記基準電池の複素インピーダンスとに基づいて、算出した前記基準電池の複素インピーダンスに対する、算出した前記蓄電池の複素インピーダンスの相対誤差を補正する前記補正値を算出する、請求項１に記載の電池状態算出装置。
　前記電池状態算出装置は、第１電池状態算出装置（２０Ａ）と、前記第１電池状態算出装置とは異なる場所に設けられる第２電池状態算出装置（２０Ｂ）とであり、
　前記第１電池状態算出装置及び前記第２電池状態算出装置は、
　前記蓄電池（１０）の正極端子（１０ａ）が接続される正極接続部（４０ａ）と、
　前記蓄電池の負極端子（１０ｂ）が接続される負極接続部（４０ｂ）と、
　前記正極接続部と前記負極接続部とを接続する前記電気経路（４１～４３）と、
を備え、
　前記第１電池状態算出装置において、前記演算部は、前記正極接続部及び前記負極接続部に前記蓄電池が接続されるとともに、前記蓄電池、前記正極接続部、前記負極接続部、前記電気経路、前記電圧検出回路及び前記電流検出回路の相対位置が固定された状態において、前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得し、取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記蓄電池の複素インピーダンスを算出し、
　前記第２電池状態算出装置において、前記演算部は、前記正極接続部及び前記負極接続部に、前記第１電池状態算出装置において複素インピーダンスの算出対象になった前記蓄電池が接続されるとともに、前記蓄電池、前記正極接続部、前記負極接続部、前記電気経路、前記電圧検出回路及び前記電流検出回路の相対位置が固定された状態において、前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得し、取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記蓄電池の複素インピーダンスを算出し、
　前記第１電池状態算出装置の前記正極接続部及び前記負極接続部に接続された状態の前記蓄電池を前記基準電池（１０）とし、前記第２電池状態算出装置の前記正極接続部及び前記負極接続部に接続された状態の前記蓄電池を対象電池（１０）とする場合において、前記第１電池状態算出装置及び前記第２電池状態算出装置の少なくとも一方の前記演算部は、前記第１電池状態算出装置において算出した前記基準電池の複素インピーダンスと、前記第２電池状態算出装置において算出した前記対象電池の複素インピーダンスとに基づいて、算出した前記基準電池の複素インピーダンスに対する、算出した前記対象電池の複素インピーダンスの相対誤差を補正する前記補正値を算出する、請求項１に記載の電池状態算出装置。
　蓄電池（１０，１０００）の正極端子と前記蓄電池の負極端子とを電気的に接続する電気経路（４１～４３，ＳＭＲ，１３０，１４０）と、
　前記電気経路及び前記蓄電池を含む閉回路に交流電流を流す交流駆動部（４２，５１，５６，６０，１３０，１３１）と、
　前記蓄電池に並列接続され、前記蓄電池の端子間電圧を検出する電圧検出回路（３１～３３）と、
　前記閉回路に流れる電流を検出する電流検出回路（４３，５２～５４）と、
を備えるシステムに適用されるプログラムにおいて、
　プロセッサ（７１，９３）に、
　前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得する処理と、
　取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記蓄電池の複素インピーダンスを算出する処理と、
　算出した前記蓄電池の複素インピーダンスと比較する基準電池の複素インピーダンスと、算出した前記蓄電池の複素インピーダンスとに基づいて、前記基準電池の複素インピーダンスに対する、算出した前記蓄電池の複素インピーダンスの相対誤差を補正する補正値（ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤ）を算出する処理と、
　算出した前記補正値に基づいて、算出した前記蓄電池の複素インピーダンスを補正する処理と、
を実行させる、プログラム。
　蓄電池（１０，１０００）の正極端子と前記蓄電池の負極端子とを電気的に接続する電気経路（４１～４３，ＳＭＲ，１３０，１４０）と、
　前記電気経路及び前記蓄電池を含む閉回路に交流電流を流す交流駆動部（４２，５１，５６，６０，１３０，１３１）と、
　前記蓄電池に並列接続され、前記蓄電池の端子間電圧を検出する電圧検出回路（３１～３３）と、
　前記閉回路に流れる電流を検出する電流検出回路（４３，５２～５４）と、
を備える電池制御装置（２０，１２０，２２０）の制御方法において、
　前記交流駆動部によって前記閉回路に交流電流が流されている場合における前記電圧検出回路の検出値及び前記電流検出回路の検出値を取得するステップと、
　取得した前記電圧検出回路の検出値及び前記電流検出回路の検出値に基づいて、前記蓄電池の複素インピーダンスを算出するステップと、
　算出した前記蓄電池の複素インピーダンスと比較する基準電池の複素インピーダンスと、算出した前記蓄電池の複素インピーダンスとに基づいて、前記基準電池の複素インピーダンスに対する、算出した前記蓄電池の複素インピーダンスの相対誤差を補正する補正値（ΔＺｋＡ，ΔＺｋＢ，ΔＺｋＣ，ΔＺｋＤ）を算出するステップと、
　算出した前記補正値に基づいて、算出した前記蓄電池の複素インピーダンスを補正するステップと、
を備える、電池状態算出装置の制御方法。