(57)【特許請求の範囲】
【請求項１】
二次電池の温度を検出する電池温度検出ステップと、
前記二次電池の充電率を求める充電率取得ステップと、
前記電池温度検出ステップにて検出した前記二次電池の温度と、前記充電率取得ステップにて取得した前記二次電池の充電率とに基づいて、前記二次電池に印加する電力の大きさを制御する制御ステップと、
前記制御ステップにて制御された同じ大きさで正負に変化する電力を前記二次電池に印加して、前記二次電池の応答電力の大きさが異なる状態の各々について、前記二次電池の電流値および電圧値を取得するデータ取得ステップと、
前記データ取得ステップにて取得した、前記二次電池の電流値および電圧値に基づいて、前記二次電池の回路モデルのモデル定数を算出するモデル定数算出ステップと、
を具備することを特徴とするモデル定数取得方法。
【請求項２】
前記データ取得ステップでは、前記二次電池の応答電力の大きさが異なる複数周波数の正弦波交流電流を印加し、各周波数について、前記二次電池の電流値および電圧値を取得することを特徴とする請求項１に記載のモデル定数取得方法。
【請求項３】
前記データ取得ステップでは、電流の向きが変化する矩形波の電力を印加し、前記二次電池の応答電力の大きさが異なるタイミングの各々で、前記二次電池の電流値および電圧値を取得することを特徴とする請求項１に記載のモデル定数取得方法。
【請求項４】
前記モデル定数算出にて、様々な劣化の程度の二次電池について前記モデル定数を算出し、得られたモデル定数に基づいて、二次電池の劣化の程度に応じて前記モデル定数の補間係数を示す補間式を算出することを特徴とする請求項１から３のいずれか一項に記載のモデル定数取得方法。
【請求項５】
二次電池の温度を検出する電池温度検出部と、
前記二次電池の充電率を求める充電率取得部と、
前記電池温度検出部が検出した前記二次電池の温度と、前記充電率取得部が取得した前記二次電池の充電率とに基づいて、前記電力印加部が前記二次電池に印加する電力の大きさを制御する制御部と、
前記制御部によって制御された同じ大きさで正負に変化する電力を前記二次電池に印加する電力印加部と、
前記二次電池の応答電力の大きさが異なる状態の各々について、前記二次電池の電流値および電圧値を取得するデータ取得部と、
前記データ取得部が取得した、前記二次電池の電流値および電圧値に基づいて、前記二次電池の回路モデルのモデル定数を算出するモデル定数算出部と、
を具備することを特徴とするモデル定数取得装置。
【請求項６】
前記電力印加部は、前記二次電池の応答電力の大きさが異なる複数周波数の正弦波交流電流を印加し、
前記データ取得部は、各周波数について、前記二次電池の電流値および電圧値を取得することを特徴とする請求項５に記載のモデル定数取得装置。
【請求項７】
前記電力印加部は、電流の向きが変化する矩形波の電力を印加し、
前記データ取得部は、前記二次電池の応答電力の大きさが異なるタイミングの各々で、前記二次電池の電流値および電圧値を取得する
ことを特徴とする請求項５に記載のモデル定数取得装置。
【請求項８】
前記モデル定数算出部は、様々な劣化の程度の二次電池について前記モデル定数を算出し、得られたモデル定数に基づいて、二次電池の劣化の程度に応じて前記モデル定数の補間係数を示す補間式を算出することを特徴とする請求項５から７のいずれか一項に記載のモデル定数取得装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、モデル定数取得方法およびモデル定数取得装置に関する。
【背景技術】
【０００２】
電池の充電率推定は、電池を運用する上で極めて重要であり、高い推定精度が求められる。この充電率推定法の１つに電流積算法があるが、測定誤差や容量誤差による誤差が蓄積されるため、簡易等価回路モデルを適用して充電率（ＳＯＣ；State Of Charge）推定を行う手法である回路モデル法が主流となっている。回路モデル法では、事前に充電率（ＳＯＣ）または開放電圧（ＯＣＶ； Open Circuit Voltage）と温度とをパラメータに、回路定数を同定算出する必要があり、代表的な方法に交流インピーダンス方がある。
【０００３】
また、特許文献１に記載の電池モデル同定方法では、Ｍ系列入力電流作成部２により周波数成分の異なるＭ系列信号を電流入力として電池に入力し、その際の電池の端子電圧を電圧センサで測定し、パラメータ推定部が、測定結果に基づいてシステム同定を行い、電池の周波数特性を算出し、算出した周波数特性に基づいて電池モデルのパラメータを同定する。
これにより、入力信号を実際の大きさとしつつ、電池の状態を変化させないように電池モデルの同定を行うことができる、とされている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−１５７４９２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
交流インピーダンス法では、モデル定数の精度が得られない、電池状態が変化する、同定に時間を要するなどの欠点がある。
また、特許文献１に記載の電池モデル同定方法では、Ｍ系列信号を生成する必要があり、装置構成が複雑になってしまう。
【０００６】
本発明は、このような事情を考慮してなされたものであり、その目的は、より簡単な構成の装置を用いて、電池状態を変化させずに、二次電池の回路モデルのモデル定数を取得することができるモデル定数取得方法およびモデル定数取得装置を提供することにある。
【課題を解決するための手段】
【０００７】
この発明は上述した課題を解決するためになされたもので、本発明の一態様によるモデル定数取得方法は、二次電池の温度を検出する電池温度検出ステップと、前記二次電池の充電率を求める充電率取得ステップと、前記電池温度検出ステップにて検出した前記二次電池の温度と、前記充電率取得ステップにて取得した前記二次電池の充電率とに基づいて、前記二次電池に印加する電力の大きさを制御する制御ステップと、前記制御ステップにて制御された同じ大きさで正負に変化する電力を前記二次電池に印加して、前記二次電池の応答電力の大きさが異なる状態の各々について、前記二次電池の電流値および電圧値を取得するデータ取得ステップと、前記データ取得ステップにて取得した、前記二次電池の電流値および電圧値に基づいて、前記二次電池の回路モデルのモデル定数を算出するモデル定数算出ステップと、を具備することを特徴とする。
【０００９】
また、本発明の一態様によるモデル定数取得方法は、上述のモデル定数取得方法であって、前記データ取得ステップでは、前記二次電池の応答電力の大きさが異なる複数周波数の正弦波交流電流を印加し、各周波数について、前記二次電池の電流値および電圧値を取得することを特徴とする。
【００１０】
また、本発明の一態様によるモデル定数取得方法は、上述のモデル定数取得方法であって、前記データ取得ステップでは、電流の向きが変化する矩形波の電力を印加し、前記二次電池の応答電力の大きさが異なるタイミングの各々で、前記二次電池の電流値および電圧値を取得することを特徴とする。
【００１１】
また、本発明の一態様によるモデル定数取得方法は、上述のモデル定数取得方法であって、前記モデル定数算出にて、様々な劣化の程度の二次電池について前記モデル定数を算出し、得られたモデル定数に基づいて、二次電池の劣化の程度に応じて前記モデル定数の補間係数を示す補間式を算出することを特徴とする。
【００１２】
また、本発明の一態様によるモデル定数取得装置は、二次電池の温度を検出する電池温度検出部と、前記二次電池の充電率を求める充電率取得部と、前記電池温度検出部が検出した前記二次電池の温度と、前記充電率取得部が取得した前記二次電池の充電率とに基づいて、前記電力印加部が前記二次電池に印加する電力の大きさを制御する制御部と、前記制御部によって制御された同じ大きさで正負に変化する電力を前記二次電池に印加する電力印加部と、前記二次電池の応答電力の大きさが異なる状態の各々について、前記二次電池の電流値および電圧値を取得するデータ取得部と、前記データ取得部が取得した、前記二次電池の電流値および電圧値に基づいて、前記二次電池の回路モデルのモデル定数を算出するモデル定数算出部と、を具備することを特徴とする。
【００１４】
また、本発明の一態様によるモデル定数取得装置は、上述のモデル定数取得装置であって、前記電力印加部は、前記二次電池の応答電力の大きさが異なる複数周波数の正弦波交流電流を印加し、前記データ取得部は、各周波数について、前記二次電池の電流値および電圧値を取得することを特徴とする。
【００１５】
また、本発明の一態様によるモデル定数取得装置は、上述のモデル定数取得装置であって、前記電力印加部は、電流の向きが変化する矩形波の電力を印加し、前記データ取得部は、前記二次電池の応答電力の大きさが異なるタイミングの各々で、前記二次電池の電流値および電圧値を取得することを特徴とする。
【００１６】
また、本発明の一態様によるモデル定数取得装置は、上述のモデル定数取得装置であって、前記モデル定数算出部は、様々な劣化の程度の二次電池について前記モデル定数を算出し、得られたモデル定数に基づいて、二次電池の劣化の程度に応じて前記モデル定数の補間係数を示す補間式を算出することを特徴とする。
【発明の効果】
【００１７】
本発明によれば、より簡単な構成の装置を用いて、電池状態を変化させずに、二次電池の回路モデルのモデル定数を取得することができる。
【図面の簡単な説明】
【００１８】
【図１】本発明の一実施形態におけるモデル定数取得装置の概略構成を示す構成図である。
【図２】同実施形態における、二次電池を模擬する電池モデルの例を示す説明図である。
【図３】同実施形態において、電力印加部が二次電池に印加する高周波の正弦波交流、および、二次電池の応答電圧の例を示す説明図である。
【図４】同実施形態において、電力印加部が二次電池に印加する低周波の正弦波交流、および、二次電池の応答電圧の例を示す説明図である。
【図５】同実施形態の電池モデルにおける、インピーダンスおよびキャパシタンスを示す変数名を表記した説明図である。
【図６】同実施形態において、電力印加部が高周波電力を印加した状態での、回路モデルの電流値および電圧値を示す変数名を表記した説明図である。
【図７】同実施形態において、電力印加部が低周波電力を印加した状態での、回路モデルの電流値および電圧値を示す変数名を表記した説明図である。
【図８】同実施形態における波高値テーブル記憶部が記憶する波高値テーブルの例を示す説明図である。
【図９】同実施形態において、モデル定数取得装置が正弦波交流電力を印加してモデル定数を算出する処理手順を示すフローチャートである。
【図１０】同実施形態において、電力印加部が印加するステップ電力の電流値と、二次電池の応答電圧値との例を示す説明図である。
【図１１】同実施形態において、電力印加部がステップ電力を印加する場合に、波高値テーブル記憶部が記憶する波高値テーブルの例を示す説明図である。
【図１２】同実施形態において、モデル定数取得装置がステップ電力を印加してモデル定数を算出する処理手順を示すフローチャートである。
【図１３】同実施形態の実験結果を示す説明図である。
【図１４】同実施形態における、二次電池の劣化度合いの各段階での抵抗値ＤＣＲに対する補正係数の例を示す説明図である。
【図１５】同実施形態において、モデル定数算出部が求めた補正係数の平均値と、容量劣化率との関係を示す説明図である。
【図１６】同実施形態において、モデル定数取得装置が、モデル定数の補間式を算出する処理手順を示すフローチャートである。
【図１７】同実施形態において、モデル定数算出部が１つの補間式のみを算出する例を示す説明図である。
【発明を実施するための形態】
【００１９】
以下、図面を参照して、本発明の実施の形態について説明する。図１は、本発明の一実施形態におけるモデル定数取得装置の概略構成を示す構成図である。同図において、モデル定数取得装置１は、電流センサ１０１と、温度センサ１０３と、電流検出部１１１と、電圧検出部１１２と、電池温度検出部１１３と、データ取得部１２０と、モデル定数算出部１３０と、制御部１４０と、電力印加部１５０と、波高値テーブル記憶部１９０と、恒温器２００とを具備する。また、モデル定数取得装置１には二次電池９００が設置されている。
【００２０】
二次電池９００は、例えばリチウムイオン電池または鉛電池などの二次電池である。二次電池９００は、モデル定数取得対象であり、二次電池９００のモデル定数をより正確に取得することで、二次電池９００の充電率（State Of Charge；ＳＯＣ）をより正確に求めることができる。このモデル定数について、図２を参照して説明する。
【００２１】
図２は、二次電池９００を模擬する電池モデルの例を示す説明図である。同図に示す電池モデルは、抵抗９０３およびキャパシタ（コンデンサ）９０４が並列に接続された回路（以下、「並列回路部分」と称する）と、抵抗９０２と、直流電源９０１とが直列に接続されて構成されている。
直流電源９０１は、抵抗値を有さない仮想的な電池であり、二次電池９００の開放電圧を示す。
【００２２】
抵抗９０２は、二次電池９００の内部抵抗のうち、直流電流と交流電流とのいずれにも作用する成分を模擬する。以下では、抵抗９０２の抵抗値をＤＣＲ［Ω（オーム）］と表記する。
抵抗９０３およびキャパシタ９０４が並列に接続された回路は、二次電池９００の内部抵抗のうち、直流電流と交流電流とで作用の度合いが異なる成分を示す。かかる抵抗成分は、例えば二次電池９００の濃度分極現象に起因して生じる。この並列接続回路は、直流電流や比較的周波数の低い電流に対しては比較的大きい抵抗値を示し、比較的周波数の高い電流に対しては比較的小さい抵抗値を示す。以下では、抵抗９０３の抵抗値をＡＣＲ［Ω］と表記し、キャパシタ９０４のキャパシタンスをＡＣＣ［Ｆ（ファラッド）］と表記する。
【００２３】
電力印加部１５０は、制御部１４０の制御に従って、二次電池９００に電力を印加する。特に、電力印加部１５０は、同じ大きさで正負に変化する電力を二次電池９００に印加する。より具体的には、電力印加部１５０は、二次電池９００の応答電力が異なる複数周波数の正弦波交流電流を印加する。
【００２４】
なお、以下では、電力印加部１５０が定電流にて二次電池９００に電力を印加して、電圧検出部１１２が二次電池９００の両端電圧値を検出することで二次電池９００の応答電力を検出する場合について説明するが、本発明の適用範囲はこれに限らない。例えば、二次電池９００が定電圧にて二次電池９００に電力を印加し、電流検出部１１１が充放電電流値（二次電池９００を流れる電流の値）を検出することで二次電池９００の応答電力を検出するようにしてもよい。
【００２５】
電流センサ１０１は、例えば電流検出トランス（Current Trans；ＣＴ）を有し、充放電電流値に応じた電流を電流検出部１１１に出力する。
温度センサ１０３は、二次電池９００に接して設置され、二次電池９００の温度（以下、「電池温度」と称する）に応じた電流を電池温度検出部１１３に出力する。
【００２６】
電流検出部１１１は、電流センサ１０１からの電流に基づいて充放電電流値を検出（測定）し、得られた電流値をデータ取得部１２０に出力する。
電圧検出部１１２は、二次電池９００の両端電圧値を検出（測定）し、得られた電圧値をデータ取得部１２０に出力する。
電池温度検出部１１３は、温度センサ１０３からの電流に基づいて二次電池９００の温度を検出し、得られた温度（を示すデータ）をデータ取得部１２０に出力する。
【００２７】
データ取得部１２０は、二次電池９００の応答電力が異なる状態の各々について、二次電池９００の電流値および電圧値を取得する。より具体的には、データ取得部１２０は、電力印加部１５０が印加する各周波数について、二次電池９００の電流値および電圧値を取得する。
【００２８】
また、データ取得部１２０は、本発明における充電率取得部の一例であり、二次電池９００の充電率を求める。例えば、データ取得部１２０は、充放電していない状態で、一定時間経過後の電池電圧（ＯＣＶ）より、二次電池９００の充電率を算出する。なお、データ取得部１２０が求める充電率は、制御部１４０が電力印加部１５０に印加させる電流値を決定するためのものである。
【００２９】
モデル定数算出部１３０は、データ取得部１２０が取得した、二次電池９００の電流値（充放電電流値）および電圧値（両端電圧値）に基づいて、二次電池９００の回路モデルのモデル定数を算出する。
制御部１４０は、電池温度検出部１１３が検出した二次電池９００の温度と、データ取得部１２０が取得した二次電池９００の充電率とに基づいて、電力印加部１５０が二次電池９００に印加する電力の大きさを制御する。
【００３０】
波高値テーブル記憶部１９０は、電力印加部１５０が二次電池９００に印加する電流の波高値を示す波高値テーブルを予め記憶している。制御部１４０は、この波高値テーブルを参照して、電力印加部１５０が二次電池９００に印加する電力の大きさを決定する。
恒温器２００は、恒温器２００自らの内部の温度を設定温度に保つ。特に、恒温器２００は、二次電池９００を格納して、二次電池９００の温度を一定に保つ。
【００３１】
次に図３および図４を参照して、電力印加部１５０が二次電池９００に印加する交流電流について説明する。
図３は、電力印加部１５０が二次電池９００に印加する比較的高周波（以下では、単に「高周波」と称する）の正弦波交流、および、二次電池９００の応答電圧の例を示す説明図である。ここで、電圧検出部１１２が検出する二次電池９００の両端電圧値が、二次電池９００の応答電圧値に該当する。
同図のグラフの横軸は、基準時刻からの経過時間にて時刻を示し、縦軸は、電力印加部１５０が印加する正弦波交流の電流値および二次電池９００からの応答の電圧値を示す。
【００３２】
また、同図の線Ｌ１１は、電力印加部１５０が印加する正弦波交流の電流値を示している。図３の例では、電力印加部１５０は、波高値が１００Ａの正弦波交流を印加している。
また、線Ｌ１２は、電圧検出部１１２が検出する二次電池９００の両端電圧値（応答電圧値）を示している。
【００３３】
図４は、電力印加部１５０が二次電池９００に印加する比較的低周波（以下では、単に「低周波」と称する）の正弦波交流、および、二次電池９００の応答電圧の例を示す説明図である。図３の場合と同様、同図のグラフの横軸は、基準時刻からの経過時間にて時刻を示し、縦軸は、電力印加部１５０が印加する正弦波交流の電流値および二次電池９００からの応答の電圧値を示す。
【００３４】
また、同図の線Ｌ２１は、電力印加部１５０が印加する正弦波交流の電流値を示している。図３場合と同様、図４の例では、電力印加部１５０は、波高値が１００Ａの正弦波交流を印加している。
また、線Ｌ１２は、電圧検出部１１２が検出する二次電池９００の両端電圧値（応答電圧値）を示している。
図３と図４とを比較すると、周波数の高い図３の場合のほうが、応答電圧値が小さくなっている。二次電池９００のキャパシタ成分（図２のキャパシタ９０４）が電流を流すことで、二次電池９００全体でのインピーダンスが小さくなったためである。
【００３５】
この電力印加部１５０が印加する正弦波交流の周波数について図５を参照して説明する。
図５は、図２で説明した電池モデルにおける、インピーダンスおよびキャパシタンスを示す変数名を表記した説明図である。
図２で説明したように、抵抗９０２の抵抗値をＤＣＲ［Ω］とし、抵抗９０３の抵抗値をＡＣＲ［Ω］とし、キャパシタ９０４のキャパシタンスをＡＣＣ［Ｆ］とする。これらと、ラプラス演算子ｓと、並列回路部分の時定数τとを用いて、等価回路モデルの電池インピーダンスＺ［Ω］は式（１）のように示される。
【００３６】
【数１】
0005717671.tif 000002
【００３７】
この時定数τよりも応答が充分短い高周波電力を二次電池９００に印加した場合、ｓ→∞と考えることができ、電池インピーダンスＺは式（２）のように示される。
【００３８】
【数２】
0005717671.tif 000003
【００３９】
一方、時定数τよりも応答が充分長い低周波電力を二次電池９００に印加した場合、ｓ→０と考えることができ、電池インピーダンスＺは式（３）のように示される。
【００４０】
【数３】
0005717671.tif 000004
【００４１】
ここで、時定数τは、二次電池（二次電池９００）の特性に応じて定まる。例えば、時定数τが１０秒程度である場合、時定数τの１０倍の周期を有する０．０１ヘルツ（Ｈｚ）の電力を比較的低周波の電力として用い、時定数τの１０分の１の周期を有する１ヘルツの電力を比較的高周波の電力として用いることが考えられる。
低周波、高周波それぞれの周波数を決定して制御部１４０に記憶させておき、制御部１４０は、当該周波数にて電力印加部１５０の出力する電力を制御する。
【００４２】
あるいは、二次電池の特性が不明の場合、様々な周波数の電力を二次電池９００に印加し、低周波側、高周波側の各々で、周波数の印加に対して応答電圧値の変化が検出されない周波数を取得することで、低周波、高周波それぞれの周波数を決定するようにしてもよい。
【００４３】
次に、図６および図７を参照して、モデル定数算出部１３０が行うモデル定数（抵抗値ＤＣＲおよびＡＣＲと、キャパシタンスＡＣＣと）の算出について説明する。
図６は、電力印加部１５０が高周波電力を印加した状態における、回路モデルの電流値および電圧値を示す変数名を表記した説明図である。同図に示すように、直流電源９０１の電圧値をＯＣＶ１１とし、抵抗９０２の電圧値をＶＤとし、抵抗９０３とキャパシタンス９０４とで構成する並列回路の電圧値をＶＣＲとする。また、二次電池９００の両端電圧値をＥ１１とし、充放電電流値をＩ１１とする。
これらの変数を用いて、電圧値ＶＤは、式（４）のように示される。
【００４４】
【数４】
0005717671.tif 000005
【００４５】
また、電圧値ＶＣＲは、式（５）のように示され、上記のように、高周波においてはｓ→∞と考えることができるので、電圧値ＶＣＲ→０と考えることができる。
【００４６】
【数５】
0005717671.tif 000006
【００４７】
そうすると、電圧値ＶＤは、式（６）のように示される。
【００４８】
【数６】
0005717671.tif 000007
【００４９】
そこで、モデル定数算出部１３０は、式（７）および式（４）に基づいて、抵抗値ＶＤＲを算出する。
【００５０】
【数７】
0005717671.tif 000008
【００５１】
但し、ＭＡＸ（ＶＤ）は、電圧値ＶＤの最大値（波高値）を示し、ＭＡＸ（Ｉ１１）は、電流値Ｉ１１の最大値（波高値）を示す
ここで、モデル定数算出部１３０は、電源電圧の電圧値ＯＣＶ１１を、以下の手順で求める。
（１）高周波電圧Ｅ１１に対応するＳＯＣ（充電率）を、ＯＣＶ−ＳＯＣテーブルを参照して求め、初期ＳＯＣとする。
（２）初期ＳＣ＋（（Ｉ１１の積分値／３６００）×１００×計算周期）／電池定格容量にて、ＳＯＣを算出する。
（３）ＳＯＣ−ＯＣＶテーブルを参照して、ＯＣＶ１１を求める。
【００５２】
抵抗値ＡＣＲおよびキャパシタンスＡＣＣについては、モデル定数算出部１３０は、低周波電力のデータから求める。
図７は、電力印加部１５０が低周波電力を印加した状態における、回路モデルの電流値および電圧値を示す変数名を表記した説明図である。同図に示すように、直流電源９０１の電圧値をＯＣＶ１２とし、抵抗９０２の電圧値をＶＤとし、抵抗９０３とキャパシタンス９０４とで構成する並列回路の電圧値をＶＣＲとする。また、充放電電流値をＩ１２とする。
これらの変数を用いて、電圧値ＶＣＲは、式（８）のように示される。
【００５３】
【数８】
0005717671.tif 000009
【００５４】
ここで、ＶＣＲの移動平均値（例えば前後１秒）をＶＣＲ２とすると、ＶＣＲ２は、式（９）のように示される。
【００５５】
【数９】
0005717671.tif 000010
【００５６】
そこで、モデル定数算出部１３０は、移動平均値ＶＣＲ２と電流値Ｉ１２とを比較して、以下の手順で抵抗値ＡＣＲおよびキャパシタンスＡＣＣを求める。
（１）移動平均ＶＣＲ２と電流値Ｉ１２との位相Ｔ１（ＡＣＲ×ＡＣＣ）を算出する。
（２）移動平均ＶＣＲ２と｛１／（１＋Ｔ１×ｓ）｝×Ｉ１２との振幅比を算出し、抵抗値ＡＣＲとする。
（３）位相Ｔと（ＡＣＲ×ＡＣＣ）／ＡＣＲより、ＡＣＣを算出する。
【００５７】
なお、モデル定数算出部１３０は、電圧値ＯＣＶ１１を求めた手順と同様にしてＯＣＶ１２を求める。
また、低周波、高周波の各々について、１サイクル目では電圧変動が大きく、正確な測定を行えないことが考えられる。そこで、低周波、高周波の各々について、モデル定数算出部１３０は、２サイクル目以降のデータにてモデル定数を算出する。
【００５８】
次に、図８を参照して、電力印加部１５０が印加する電力（本実施形態では電流）の大きさについて説明する。
図８は、波高値テーブル記憶部１９０が記憶する波高値テーブルの例を示す説明図である。同図に示す波高値テーブルは、充電率（ＳＯＣ）と、電池温度の各値について、低周波（同図の例では１００ｓ、すなわち周期１００秒）電力として印加する電流値と、高周波（同図の例では１ｓ、すなわち周期１秒）電力として印加する電流値とを格納している。
【００５９】
ここで、二次電池は、一般に、充電率と電池温度とによって印加可能な電力の大きさが異なる。一方、測定時のノイズの影響を低減させるためには、二次電池９００に、なるべく大きい電流を印加することが好ましい。そこで、波高値テーブル記憶部１９０が、充電率および電池温度毎に、二次電池９００に印加可能な電流値を波高値テーブルとして記憶しておき、制御部１４０は、電池温度検出部１１３が検出した電池温度と、データ取得部１２０が求めた充電率とに基づいて、波高値テーブルから電流値を読み出し、当該電流値の電力を電力印加部１５０に印加させる。
【００６０】
次に、図９を参照して、モデル定数取得装置１の動作について説明する。
図９は、モデル定数取得装置１が正弦波交流電力を印加してモデル定数を算出する処理手順を示すフローチャートである。モデル定数取得装置１は、二次電池９００を設置された状態で、モデル定数算出を指示するユーザ操作を検出すると、同図の処理を開始する。
【００６１】
同図の処理において、まず、電池温度検出部１１３が、温度センサ１０３からの電流に基づいて電池温度を検出し、データ取得部１２０を介して制御部１４０に出力する（ステップＳ１０１）。
また、データ取得部１２０は、二次電池９００の充電率を算出して制御部１４０に出力する（ステップＳ１０２）。例えば、データ取得部１２０は、充放電していない状態で、一定時間経過後の電池電圧（ＯＣＶ）より、二次電池９００の充電率を算出する。
【００６２】
次に、制御部１４０は、電池温度検出部１１３が検出した電池温度と、データ取得部１２０が算出した充電率とに基づいて、波高値テーブル記憶部１９０から電流値を読み出すことで、電力印加部１５０に印加させる電流値（波高値）を決定する（ステップＳ１０３）。
【００６３】
そして、制御部１４０は、電力印加部１５０を制御して高周波の正弦波交流電力を二次電池９００に印加させ、データ取得部１２０は、電流値や電圧値のデータを取得する（ステップＳ１０４）。より具体的には、制御部１４０は、ステップＳ１０３で決定した電流値、および、予め設定されている高周波の周波数にて正弦波交流電力を二次電池９００に印加するよう電力印加部１５０を制御し、電力印加部１５０は、制御部１４０の制御に従って二次電池９００に電力を印加する。そして、データ取得部１２０は、電流検出部１１１が検出する充放電電流値および電圧検出部１１２が検出する二次電池９００の両端電圧値を取得する。
【００６４】
データ取得部１２０がデータを取得すると、制御部１４０は、電力印加部１５０を制御して低周波の正弦波交流電力を二次電池９００に印加させ、データ取得部１２０は、電流値や電圧値のデータを取得する（ステップＳ１０５）。より具体的には、制御部１４０は、ステップＳ１０３で決定した電流値、および、予め設定されている低周波の周波数にて正弦波交流電力を二次電池９００に印加するよう電力印加部１５０を制御し、電力印加部１５０は、制御部１４０の制御に従って二次電池９００に電力を印加する。そして、データ取得部１２０は、電流検出部１１１が検出する充放電電流値および電圧検出部１１２が検出する二次電池９００の両端電圧値を取得する。
【００６５】
データ取得部１２０がデータを取得すると、モデル定数算出部１３０は、データ取得部１２０が取得したデータに基づいて、上述したようにモデル定数を算出する（ステップＳ１０６）。
その後、同図の処理を終了する。モデル定数算出部１３０が算出したモデル定数は、例えば、モデル定数取得装置１のユーザが、二次電池９００のＢＭＵ（Battery Management Unit）のＲＯＭ（Read Only Memory）に搭載する（記憶させる）。
【００６６】
以上のように、電力印加部１５０は、正弦波交流を印加し、モデル定数算出部１３０は、当該正弦波交流を印加した状態における測定値に基づいてモデル定数を算出する。
従って、電力印加部１５０として、正弦波交流を印加可能な電源回路を用いればよく、装置構成を簡単にすることができる。
【００６７】
また、電力印加部１５０が、同じ大きさで正負に変化する電力の一種である正弦波交流を二次電池９００に印加するので、充電電流と放電電流とが同じになり、二次電池９００の状態（特に、充電率や開放電圧値）が変化しない。従って、より正確に測定を行うことができ、より高精度なモデル定数を得ることができる。
例えば、当該モデル定数を、二次電池９００のＢＭＵのＲＯＭに搭載しておくことで、当該ＢＭＵは、より正確に二次電池９００の充電率を算出することができる。
【００６８】
また、制御部１４０が、電池温度および充電率に基づいて、電力印加部１５０が印加する電流の大きさを制御するので、二次電池９００の許容範囲内で比較的大電力を印加することができる。これにより、ノイズの影響を低減させ、より高精度なモデル定数を得ることができる。
その際、恒温器２００が電池温度を一定に保つことで、制御部１４０は、より正確に電力印加部１５０が印加する電流の大きさを制御できる。
【００６９】
なお、電力印加部１５０が印加する電力は正弦波交流に限らない。例えば、電力印加部１５０が、同じ波高値および時間幅を有する正負それぞれのステップ波の電力（以下、「ステップ電力」と称する）を交互に印加するようにしてもよい。
【００７０】
図１０は、電力印加部１５０が印加するステップ電力の電流値と、二次電池９００の応答電圧値との例を示す説明図である。同図（Ａ）および（Ｂ）のグラフの横軸は、基準時刻（電力印加部１５０がステップ電力を印加したタイミング）からの経過時間にて時刻を示す。また、同図（Ａ）の縦軸は電流値を示し、同図（Ｂ）の縦軸は電圧値を示す。
【００７１】
また、同図（Ａ）は、電力印加部１５０が印加するステップ電力の電流値を示している。
電力印加部１５０は、時刻０（基準時刻）以前の時間帯において、電流値０のステップ電力を二次電池９００に印加している（すなわち、二次電池９００に電力を印加していない）。以下では、この基準時刻以前の時間帯において電力印加部１５０が印加する電流値０のステップ電力、および、基準時刻以前の時間帯を、いずれもステップ番号Ｎ１にて表記する。
【００７２】
その後、電力印加部１５０は、時刻０から時刻ｔ１までの時間帯において、電流値Ｉ２１のステップ電力を二次電池９００に印加している。以下では、この時刻０から時刻ｔ１までの時間帯において電力印加部１５０が印加する電流値Ｉ２１のステップ電力、および、時刻０から時刻ｔ１までの時間帯を、いずれもステップ番号Ｎ２にて表記する。ステップ番号Ｎ２の電力は、電流値Ａが正の電力、すなわち充電電力となっている。
【００７３】
さらにその後、電力印加部１５０は、時刻ｔ１以後の時間帯において、電流値０のステップ電力を二次電池９００に印加している（すなわち、二次電池９００に電力を印加していない）。以下では、この時刻ｔ１以後の時間帯において電力印加部１５０が印加する電流値０のステップ電力および、時刻ｔ１以後の時間帯を、いずれもステップ番号Ｎ３にて表記する。
【００７４】
なお、電力印加部１５０は、ステップ番号Ｎ３の時間帯にて二次電池９００の応答電圧値が安定した後、ステップ番号Ｎ２で印加したのと同じ大きさの電力かつ同じ時間幅で、電流が逆向きの電力を二次電池９００に印加する（すなわち、二次電池９００を放電させる）。そして、モデル定数算出部１３０は、放電側についても、以下に説明する充電側の場合と同様にモデル定数を算出し、充電側と放電側とでモデル定数の平均を取る。
【００７５】
一方、図１０（Ｂ）は、二次電池９００の応答電圧値を示している。
同図に示すように、電力印加部１５０が電力を印加すると、印加を開始したタイミングで応答電圧値が上昇し、その後も応答電圧値が上昇した後安定状態に至る。二次電池９００のキャパシタ成分（図２のキャパシタ９０４で模擬される）に電位が蓄えられるためである。
なお、二次電池９００の両端電圧値にノイズ（スパイク成分）が含まれている場合、例えば、電圧検出部１１２がローパスフィルタ（Low Pass Filter）を用いてノイズを除去する。
【００７６】
二次電池９００の応答電圧値について、より詳細に説明すると、ステップ番号Ｎ１では、二次電池９００が前回放電を行ってから時間が充分経過しており、二次電池９００の応答電圧値（両端電圧値）は電圧値ＯＣＶ２０にて安定している。
その後、電力印加部１５０がステップ番号Ｎ２の電力の印加を開始した時刻０において、二次電池９００の応答電圧値は電圧値Ｅ２１まで上昇し、その後さらに電圧値Ｅ２２まで上昇して安定する。
【００７７】
そして、電力印加部１５０がステップ番号Ｎ２の電力の印加を終了した時刻ｔ１において、二次電池９００の応答電圧は、時刻０における電圧上昇分と同じ分だけ降下して電圧値Ｅ２３となり（すなわち、Ｅ２１−ＯＣＶ２０＝Ｅ２２−Ｅ２３）、その後も応答電圧値が降下した後、電圧値Ｅ２４にて安定している。ここで、電圧値Ｅ２４は、二次電池９００のキャパシタ成分（図２のキャパシタ９０４で模擬される）に残存する電位の分だけ電圧値ＯＣＶ２０よりも大きくなっている。
これらの電圧値および電流値を用いて、抵抗値ＤＣＲは、式（１０）のように示される
【００７８】
【数１０】
0005717671.tif 000011
【００７９】
そこで、モデル定数算出部１３０は、電圧検出部１１２が検出した電圧値ＯＣＶ２０およびＥ２１と、ステップ番号Ｎ１における電流値Ｉ２１とを式（１０）に代入して、抵抗値ＤＣＲを算出する。
【００８０】
また、図１０に示す時定数τは、式（１）における並列回路部分の時定数τと同一であり、時間ｔ１から時刻ｔ２までの時間（τ＝ｔ２−ｔ１）を示す。ここで、時刻ｔ２は、ステップ番号Ｎ３の時間帯における応答電圧の降下分Ｖｐ１（Ｖｐ１＝Ｅ２３−Ｅ２４）のうち、１／ｅに相当する電圧値Ｖｐ２の電位が二次電池９００のキャパシタ成分に蓄電されている状態となった時刻である。すなわち、時刻ｔ２は、二次電池９００の両端電圧値が、Ｅ２４＋Ｖｐ１×（１／ｅ）に等しくなった時刻である。ここで、「ｅ」はネイピア数（すなわち、自然対数の底）を示す。
データ取得部１２０は、電圧検出部１１２が検出する二次電池９００の両端電圧に基づいて、時刻ｔ２を取得（測定）する。
【００８１】
また、時定数τを用いて、抵抗９０３（図２）の抵抗値ＡＣＲとキャパシタ９０４のキャパシタンスＡＣＣとの関係は、式（１１）のように示される。
【００８２】
【数１１】
0005717671.tif 000012
【００８３】
また、ステップ番号Ｎ３の時間帯における応答電圧の降下分Ｖｐ１は、式（１２）のように示される。
【００８４】
【数１２】
0005717671.tif 000013
【００８５】
これら式（１１）および（１２）に基づいて、抵抗値ＡＣＲは、式（１３）のように示される。
【００８６】
【数１３】
0005717671.tif 000014
【００８７】
また、式（１１）より、キャパシタンスＡＣＣは、式（１４）のように示される。
【００８８】
【数１４】
0005717671.tif 000015
【００８９】
そこで、モデル定数算出部１３０は、データ取得部１２０が取得する電圧値Ｅ２２およびＥ２３に基づいて応答電圧の降下分Ｖｐ１を算出し、このＶｐ１と、ステップ番号Ｎ２における電流値Ｉ１２と、ステップ番号Ｎ２からＮ３への切り替わりの時刻ｔ１と、データ取得部１２０が検出する時刻ｔ２とを、式（１３）に代入して抵抗値ＡＣＲを算出する。さらに、モデル定数算出部１３０は、時刻ｔ１およびｔ２と、算出した抵抗値ＡＣＲとを式（１４）に代入してキャパシタンスＡＣＣを算出する。
【００９０】
また、電力印加部１５０が交流電力を印加する場合と同様、電力印加部１５０は波高値テーブル記憶部１９０は、波高値テーブルを記憶し、制御部１４０は、当該波高値テーブルを参照して電力印加部１５０が印加する電圧の電流値を制御する。
図１１は、電力印加部１５０がステップ電力を印加する場合に、波高値テーブル記憶部１９０が記憶する波高値テーブルの例を示す説明図である。同図に示す波高値テーブルは、充電率（ＳＯＣ）と、電池温度の各値について、ステップ電力として印加する電流値を格納している。図８で説明したのと同様、電力印加部１５０がステップ電力を印加する場合も、波高値テーブル記憶部１９０が、充電率および電池温度毎に、二次電池９００に印加可能な電流値を波高値テーブルとして記憶しておく。そして、制御部１４０は、電池温度検出部１１３が検出した電池温度と、データ取得部１２０が求めた充電率とに基づいて、波高値テーブルから電流値を読み出し、当該電流値の電力を電力印加部１５０に印加させることで、二次電池９００に対して可及的に大きい電力を印加して、ノイズの影響を低減させる。
【００９１】
次に、図１２を参照して、モデル定数取得装置１の動作について説明する。
図１２は、モデル定数取得装置１がステップ電力を印加してモデル定数を算出する処理手順を示すフローチャートである。モデル定数取得装置１は、二次電池９００を設置された状態で、モデル定数算出を指示するユーザ操作を検出すると、同図の処理を開始する。
【００９２】
同図の処理において、ステップＳ２０１〜Ｓ２０３は、図１のステップＳ１０１〜Ｓ１０３と同様である。なお、データ取得部１２０は、制御部１４０が決定した電流値を電流値Ｉ２１として取得する。あるいは、データ取得部１２０が、電流検出部１１１の検出した電流値を電流値Ｉ２１として取得するようにしてもよい。
【００９３】
次に、電力印加部１５０は、制御部１４０の制御に従ってステップ電力（図１０のステップ番号Ｎ２のステップ電力に相当する）の印加を開始する（ステップＳ２０４）。
そして、二次電池９００の両端電圧値が安定（静定）した後、電力印加部１５０は、制御部１４０の制御に従って、ステップＳ２０４で開始したステップ電力の印加を終了し、データ取得部１２０は、印加終了時刻ｔ１を取得する（ステップＳ２０５）。
【００９４】
その後、データ取得部１２０は、二次電池９００の両端電圧値を時刻と対応付けて記録し、二次電池９００の両端電圧が安定（静定）した後、電圧値の降下分Ｖｐ１と、時定数τに対応する時刻ｔ２（図１０参照）とを算出する（ステップＳ２０６）。あるいは、モデル定数取得装置１が電圧値の降下分Ｖｐ１を予め測定しておくなど取得しておき、時刻ｔ２をリアルタイムで検出するようにしてもよい。
【００９５】
次に、モデル定数算出部１３０は、データ取得部１２０が取得したデータに基づいて、上述したようにモデル定数を算出する（ステップＳ２０７）。
次に、モデル定数取得装置１は、逆向きの充放電電流についても同様にモデル定数を算出する（ステップＳ２０８）。
【００９６】
そして、モデル定数取得装置１（モデル定数算出部１３０）は、ステップＳ２０７で算出したモデル定数と、ステップＳ２０８で算出したモデル定数との平均値を算出して、最終的なモデル定数とする（ステップＳ２０９）。
その後、同図の処理を終了する。
【００９７】
以上のように、電力印加部１５０は、電流の向きが変化する矩形波の電力を印加し、データ取得部１２０は、二次電池９００の応答電力の大きさが異なるタイミングの各々で、二次電池９００の電流値および電圧値を取得する。そして、モデル定数算出部１３０は、データ取得部１２０が取得した電流値および電圧値に基づいて、モデル定数を算出する。
従って、電力印加部１５０として、ステップ電力を印加可能な電源回路を用いればよく、装置構成を簡単にすることができる。
【００９８】
また、電力印加部１５０が、同じ大きさで正負に変化するステップ電力を二次電池９００に印加するので、充電電流と放電電流とが同じになり、二次電池９００の状態（特に、充電率や開放電圧値）が変化しない。従って、より正確に測定を行うことができ、より高精度なモデル定数を得ることができる。
例えば、当該モデル定数を、二次電池９００のＢＭＵのＲＯＭに搭載しておくことで、当該ＢＭＵは、より正確に二次電池９００の充電率を算出することができる。
また、モデル定数算出部１３０が、充電電流、放電電流の各々についてモデル定数を算出して平均を取ることで、より高精度なモデル定数を得ることができる。
【００９９】
また、制御部１４０が、電池温度および充電率に基づいて、電力印加部１５０が印加する電流の大きさを制御するので、二次電池９００の許容範囲内で比較的大電力を印加することができる。これにより、ノイズの影響を低減させ、より高精度なモデル定数を得ることができる。
その際、恒温器２００が電池温度を一定に保つことで、制御部１４０は、より正確に電力印加部１５０が印加する電流の大きさを制御できる。
【０１００】
なお、本願発明者は、正弦波交流電力の印加によるパラメータ定数の算出と、ステップ電力の印加によるパラメータ定数の算出について実験を行った。
図１３は、本実施形態の実験結果を示す説明図である。
同図（Ａ）、（Ｂ）、（Ｃ）、（Ｄ）は、正弦波交流電力の印加における、それぞれ、抵抗値ＤＣＲ、抵抗値ＡＣＲ、キャパシタンスＡＣＣ、時定数τを示している。また、同図（Ｅ）、（Ｆ）、（Ｇ）、（Ｈ）は、ステップ電力の印加における、それぞれ、抵抗値ＤＣＲ、抵抗値ＡＣＲ、キャパシタンスＡＣＣ、時定数τを示している。
【０１０１】
同図において、抵抗値ＤＣＲ、抵抗値ＡＣＲ、キャパシタンスＡＣＣおよび時定数τのいずれについても、正弦波交流電力の印加とステップ電力の印加とで同様の傾向が得られているように、正弦波交流電力の印加とステップ電力の印加とのそれぞれで良好な結果が得られ、本発明が有効であることが確認された。
【０１０２】
なお、モデル定数算出部１３０が、二次電池９００の劣化に応じたモデル定数の補正係数を算出するようにしてもよい。
例えば、モデル定数取得装置１のユーザが、モデル定数取得対象の二次電池９００と同一規格（同レシピ）かつ製造ロッド同一の二次電池を複数用意し、寿命サイクル試験を行って段階的な劣化度合いの二次電池を生成しておく。そして、モデル定数取得装置１は、劣化度合いの各段階の二次電池に対して、上述した正弦波電力ないしステップ電力を印加してのモデル定数算出を行い、得られたモデル定数に基づいて補正係数を算出する。
例えば、モデル定数算出部１３０は、二次電池の劣化によるモデル定数の変化において支配的である抵抗値ＤＣＲについて補正係数を算出する。
【０１０３】
図１４は、二次電池の劣化度合いの各段階における抵抗値ＤＣＲに対する補正係数の例を示す説明図である。同図（Ａ）、（Ｂ）、（Ｃ）は、それぞれ、容量劣化率６．６％（パーセント）、１６．４％、２８．１％の各々の二次電池における補正係数を示している。ここで、容量劣化率とは、二次電池の初期状態（劣化していない状態）における容量を１００％として、減少した容量（劣化した二次電池の容量から、初期状態における容量を減算した差）の割合を％で示した値である。
【０１０４】
同図（Ａ）、（Ｂ）、（Ｃ）の各々のグラフの横軸は、いずれも温度を示し、縦軸は、いずれも抵抗値ＤＣＲに対する補正係数（抵抗値ＤＣＲに乗算すべき係数）の値を示している。モデル定数算出部１３０は、劣化した状態におけるＤＣＲを、初期状態における抵抗値ＤＣＲで除算して、当該補正係数を算出する。
また、同図（Ａ）、（Ｂ）、（Ｃ）の各々には、充電率２０％から８０％までの１０％刻みの各値について、電池温度１５℃（度）、２５℃、４０℃の各々における補正係数と、当該補正係数の温度毎の平均値が示されている。線Ｌ３１、Ｌ３２、Ｌ３３は、いずれも温度毎の平均値を結んだ線である。
【０１０５】
図１４に示すように、容量劣化率いのいずれにおいても、また、電池温度のいずれにおいても、補正係数のばらつきは比較的小さい。そこで、モデル定数算出部１３０は、容量劣化率および電池温度の各々について補正係数の平均を算出し、当該平均値にて補正係数を近似する。
【０１０６】
さらに、モデル定数算出部１３０は、補正係数の直線近似を行う。
図１５は、モデル定数算出部１３０が求めた補正係数の平均値と、容量劣化率との関係を示す説明図である。同図のグラフの横軸は容量劣化率を示し、縦軸は抵抗値ＤＣＲに対する補正係数の値を示している。また、線Ｌ４１、Ｌ４２、Ｌ４３は、それぞれ、電池温度１５℃、２５℃、４０℃における補正係数の平均値を結んだ線である。
【０１０７】
同図に示すように、線Ｌ４１、Ｌ４２、Ｌ４３は、いずれもおおよそ直線的であり、モデル定数算出部１３０は、これら線Ｌ４１、Ｌ４２、Ｌ４３の各々について（すなわち、電池温度毎に）、例えば最小二乗法を用いて直線近似を行う。同図の線Ｌ５１、Ｌ５２、Ｌ５３は、それぞれ、線Ｌ４１、Ｌ４２、Ｌ４３を直線近似して得られた線の例を示している。
【０１０８】
以上のようにして、モデル定数算出部１３０は、電池温度１５℃の場合の補間式として、式Ｌ５１を示す式（１５）を算出する。
【０１０９】
【数１５】
0005717671.tif 000016
【０１１０】
この式（１５）の変数ｘに二次電池の容量劣化率を代入することで、変数ｙの値として、抵抗値ＤＣＲに対する補正係数を得られる。
また、モデル定数算出部１３０は、電池温度２５℃の場合の補間式として、式Ｌ５２を示す式（１６）を算出する。
【０１１１】
【数１６】
0005717671.tif 000017
【０１１２】
また、モデル定数算出部１３０は、電池温度４０℃の場合の補間式として、式Ｌ５３を示す式（１７）を算出する。
【０１１３】
【数１７】
0005717671.tif 000018
【０１１４】
以上のように、モデル定数算出部１３０は、二次電池の劣化によるモデル定数の変化において支配的である抵抗値ＤＣＲについて補正係数を算出する。従って、二次電池の容量が劣化している場合にも、より正確な補正係数を得ることができる。
例えば、当該モデル定数を、二次電池９００のＢＭＵのＲＯＭに搭載しておき、ＢＭＵが当該補正係数に基づいてモデル定数を補正することで、当該ＢＭＵは、より正確に二次電池９００の充電率を算出することができる。
【０１１５】
また、モデル定数算出部１３０は、式（１５）〜（１７）のように直線を表す式にて補間式を取得する。従って、例えば二次電池９００のＢＭＵのＲＯＭ当該補間式を搭載する場合などに、傾きと切片を記憶する容量があればよく、記憶容量が少なくて済む。
また、補間式から補間係数を算出する際も、簡単な式で算出できるので、例えば二次電池９００のＢＭＵなど補間係数を算出する装置は、軽い計算量で済み、また、より迅速に補間係数を得ることができる。
【０１１６】
なお、電池温度に合致する補間式が無い場合、例えば二次電池９００のＢＭＵなど補間係数を算出する装置は、補間式から得られる補正係数を、電池温度に関して補間して用いる。
例えば、ＢＭＵが電池温度１５℃、２５℃、４０℃の各々について補間式を記憶しており、測定した電池温度が３５℃の場合、ＢＭＵは、まず、３５℃に近い温度２５℃および４０℃について補間係数を算出する。そして、２５℃の補間係数と、４０℃の補間係数とを、式（１８）に示すように、それぞれ、温度差に反比例した比で重み付け平均して補間係数を算出する。
【０１１７】
【数１８】
0005717671.tif 000019
【０１１８】
次に、図１６を参照して、モデル定数取得装置１の動作について説明する。
図１６は、モデル定数取得装置１が、モデル定数の補間式を算出する処理手順を示すフローチャートである。モデル定数取得装置１は、補間式の算出を指示するユーザ操作を検出すると、同図の処理を開始する。
【０１１９】
同図の処理において、まず、モデル定数取得装置１は、容量劣化率が異なる二次電池の設置を順に受けて、それぞれについて、また、設定された電池温度および充電率の各々について、上述した正弦波交流電力の印加、あるいは、ステップ電力の印加にて、補正係数を算出する（ステップＳ３０１）。
【０１２０】
次に、モデル定数取得装置１（モデル定数算出部１３０）は、容量劣化率毎、かつ、電池温度毎に、補正係数の平均値（充電率間における平均値）を算出する（ステップＳ３０２）。
さらに、モデル定数取得装置１（モデル定数算出部１３０）は、電池温度毎に、容量劣化率と補正係数の関係を直線近似し、近似した直線を示す式を算出して補正係数の補間式とする（ステップＳ３０３）。
その後、同図の処理を終了する。
【０１２１】
なお、空調の効いた室内など、環境温度の変化が小さい用途の二次電池に対しては、モデル定数算出部１３０が当該環境温度に応じた１つの補間式のみを算出するようにしてもよい。
図１７は、モデル定数算出部１３０が１つの補間式のみを算出する例を示す説明図である。同図のグラフの横軸は容量劣化率を示し、縦軸は抵抗値ＤＣＲに対する補正係数の値を示している。また、容量劣化率−７％、−１６％、−２８％の各々について、環境温度の範囲で３点の補正係数が得られているのに対し、モデル定数算出部１３０は、１本の線Ｌ６１を示す式（１９）を算出している。
【０１２２】
【数１９】
0005717671.tif 000020
【０１２３】
このように、モデル定数算出部１３０が、１つの補間式を算出することで、例えば二次電池９００のＢＭＵなど補間係数を算出する装置が当該補間式を記憶する際に必要な記憶容量が少なくて済む。
また、補間式から補間係数を算出する際も、簡単な式で算出できるので、例えば二次電池９００のＢＭＵなど補間係数を算出する装置は、軽い計算量で済み、また、より迅速に補間係数を得ることができる。
【０１２４】
なお、データ取得部１２０やモデル定数算出部１３０や制御部１４０の全部または一部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各部の処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【０１２５】
以上、本発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計変更等も含まれる。
【符号の説明】
【０１２６】
１  モデル定数取得装置
１０１  電流センサ
１０３  温度センサ
１１１  電流検出部
１１２  電圧検出部
１１３  電池温度検出部
１２０  データ取得部
１３０  モデル定数算出部
１４０  制御部
１５０  電力印加部
１９０  波高値テーブル記憶部
２００  恒温器
９００  二次電池
【図１】
0005717671.tif 000021
【図２】
0005717671.tif 000022
【図３】
0005717671.tif 000023
【図４】
0005717671.tif 000024
【図５】
0005717671.tif 000025
【図６】
0005717671.tif 000026
【図７】
0005717671.tif 000027
【図８】
0005717671.tif 000028
【図９】
0005717671.tif 000029
【図１０】
0005717671.tif 000030
【図１１】
0005717671.tif 000031
【図１２】
0005717671.tif 000032
【図１３】
0005717671.tif 000033
【図１４】
0005717671.tif 000034
【図１５】
0005717671.tif 000035
【図１６】
0005717671.tif 000036
【図１７】
0005717671.tif 000037