(57)【要約】【課題】二次電池の特性、特に内部パラメータを高精度に推定し得る装置を提供する。【解決手段】推定装置１０は、二次電池の内部パラメータを用いた二次電池の電気化学反応モデルに基づくシミュレーションモデル、及び内部パラメータのとり得る値の範囲を示す事前分布を記憶する。内部パラメータを設定し、シミュレーションモデルを用いて電流データに対応する電圧データを算出して事前分布及び実測電圧データを用いて評価することで、実測電圧データに適合する内部パラメータを推定する。内部パラメータは、物理量に基づくパラメータと、関数に基づくパラメータを含み、物理量に基づくパラメータの事前分布は確率分布関数とし、関数に基づくパラメータの事前分布は実測値に対する関数のフィッティング度合いとする。【選択図】図１
2015135286.tif 000002
【特許請求の範囲】
【請求項１】
二次電池の内部パラメータを用いた二次電池の電気化学反応モデルに基づくシミュレーションモデルを記憶するモデル記憶手段と、
前記内部パラメータのとり得る値の範囲を示す事前分布を記憶する事前分布記憶手段と、
前記内部パラメータを設定し、前記シミュレーションモデルを用いて電流データに対応する電圧データを算出して前記事前分布及び実測電圧データを用いて評価することで、前記実測電圧データに適合する前記内部パラメータを推定する処理手段と、
を備えることを特徴とする二次電池の特性推定装置。
【請求項２】
前記内部パラメータは、物理量に基づくパラメータと、関数に基づくパラメータを含み、
前記物理量に基づくパラメータの事前分布は確率分布関数であり、関数に基づくパラメータの事前分布は実測値に対する関数のフィッティング度合いであることを特徴とする請求項１記載の二次電池の特性推定装置。
【請求項３】
前記処理手段は、前記内部パラメータのベクトルｘを設定したときに前記シミュレーションモデルを用いて算出された電圧データをｇ（ｘ）、実測電圧データをｙｔ、内部パラメータの前記事前分布をｐ（ｘ）としたときに、
スコア＝ｐ（ｙ｜ｘ）ｐ（ｘ）
但し、σを実測電圧データの標準偏差として、
【数１】
2015135286.tif 000019
を算出することで評価し、前記スコアが最大となる前記ベクトルｘを推定することを特徴とする請求項１，２のいずれかに記載の二次電池の特性推定装置。
【請求項４】
前記内部パラメータは、開回路電位、イオン導電率、交換電流密度、イオン拡散定数、初期イオン濃度、最大イオン濃度の少なくともいずれかであることを特徴とする請求項１〜３のいずれかに記載の二次電池の特性推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は二次電池の特性推定装置に関し、特に二次電池の内部パラメータを推定する装置に関する。
【背景技術】
【０００２】
従来から、コンピュータを用いた電池シミュレーションにより、二次電池の内部パラメータを推定する技術が提案されている。
【０００３】
下記の非特許文献１には、電気化学反応モデルに基づく電池シミュレーションにおいて、モデルの一部のパラメータ（初期ＳＯＣ、反応速度に関わる係数、拡散係数）を、充放電データを用いて出力電圧の推定値と実測値の差が最小になるように非線形最小２乗法で推定することが記載されている。
【０００４】
また、非特許文献２には、電気化学反応モデルに基づく電池シミュレーションにおいて、電気化学反応モデルの全パラメータを、充放電データを用いて出力電圧と実測値の差が最小になるように遺伝的アルゴリズム（Genetic Algorithm）により推定することが記載されている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】S. Santhanagopalan, Q. Guo, and R.E. White “Parameter Estimation and Model Discrimination for a Lithium-Ion Cell,” Journal of ElectroChem., 154(3), A198-A206(2007)
【非特許文献２】J.C. Forman, S.J. Moura, J.L. Stein, and H.K. Fathy, “Genetic identification and fisher identifiability analysis of the Doyle-Fuller-Newman model from experimental cycling of a LiFePO4 Cell,” Journal of Power Sources, 210, 263-275, (2012)
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、非特許文献１では、電池の電気化学反応モデルの一部のパラメータを推定できるものの、全てのパラメータを非線形最小２乗法で同様に推定すると、パラメータの範囲を規定する制約がないためパラメータが想定外の値に収束してしまう恐れがあり、推定精度が担保されない問題がある。
【０００７】
また、非特許文献２では、電池の電気化学反応モデルの全てのパラメータを推定できるものの、ＯＣＰ（開回路電位）−ＳＯＣ（充電状態）曲線のパラメータは関数のパラメータではなく、曲線を構成する点を直接推定しており、各点に対してとり得る範囲を設定できるが点の数だけ推定すべきパラメータ数が多くなり計算リソースを多く要してしまう。さらに、各点に対する拘束条件が考慮されておらず、適切な値に収束する保証もない。
【０００８】
本発明の目的は、二次電池の特性、特に内部パラメータを高精度に推定し得る装置を提供することにある。
【課題を解決するための手段】
【０００９】
本発明は、二次電池の内部パラメータを用いた二次電池の電気化学反応モデルに基づくシミュレーションモデルを記憶するモデル記憶手段と、前記内部パラメータのとり得る値の範囲を示す事前分布を記憶する事前分布記憶手段と、前記内部パラメータを設定し、前記シミュレーションモデルを用いて電流データに対応する電圧データを算出して前記事前分布及び実測電圧データを用いて評価することで、前記実測電圧データに適合する前記内部パラメータを推定する処理手段とを備えることを特徴とする二次電池の特性推定装置である。
【００１０】
本発明の１つの実施形態では、前記内部パラメータは、物理量に基づくパラメータと、関数に基づくパラメータを含み、前記物理量に基づくパラメータの事前分布は確率分布関数であり、関数に基づくパラメータの事前分布は実測値に対する関数のフィッティング度合いであることを特徴とする。
【００１１】
本発明では、二次電池の内部パラメータが想定外の値をとって振動し、収束しない問題を解消するために、内部パラメータに対する拘束条件を事前分布として与える。このとき、内部パラメータには物理量に基づくパラメータと関数に基づくパラメータが存在し得るが、物理量に基づくパラメータについてはガウス分布等の確率分布関数で事前分布を規定し、関数に基づくパラメータについては関数のフィッティング度合いで事前分布を規定することで、内部パラメータを適切な値に収束させることができる。
【００１２】
なお、内部パラメータは複数存在し得るところ、これを内部パラメータベクトルｘとし、シミュレーションモデルを用いて算出された電圧データをｇ（ｘ）、実測電圧データをｙｔとすると、
ｆ（ｘ）＝Σ（ｙｔ−ｇ（ｘ））２
が最小となるように内部パラメータベクトルｘを推定すればよいが、これは、ｙが生じた後のｘの事後確率ｐ（ｘ｜ｙ）が最大となるようにｘを推定することに等しい。ベイズの定理に従えば、ｐ（ｘ｜ｙ）はｐ（ｙ｜ｘ）ｐ（ｘ）に比例するから、ｐ（ｙ｜ｘ）ｐ（ｘ）が最大となるように推定すればよく、ｐ（ｘ）として内部パラメータベクトルｘの事前分布が与えられる。
【発明の効果】
【００１３】
本発明によれば、事前分布を用いることにより二次電池の内部パラメータを適切な値に収束させ、高精度に推定することができる。
【図面の簡単な説明】
【００１４】
【図１】実施形態の概念構成図である。
【図２】実施形態の概念構成図である。
【図３】実施形態の処理説明図である。
【図４】実施形態の推定装置の構成ブロック図である。
【図５】実施形態の処理フローチャートである。
【図６】物理量に基づくパラメータの事前分布を示す説明図である。
【図７】関数に基づくパラメータの事前分布を示す説明図である。
【発明を実施するための形態】
【００１５】
以下、図面に基づき本発明の実施形態について説明する。
【００１６】
図１に、本実施形態における二次電池の特性を推定する概略構成図を示す。二次電池の特性としての内部パラメータを推定するための推定装置１０はコンピュータで構成され、この推定装置１０には、電池シミュレーションモデルと内部パラメータ事前分布がプログラムとして与えられる。電池シミュレーションモデル及び内部パラメータ事前分布は、推定装置１０のプログラムメモリに記憶される。
【００１７】
他方、推定装置１０には、電池充放電データが入力される。推定装置１０は、与えられた電池シミュレーションモデルと内部パラメータ事前分布を用いて、電池充放電データから二次電池の内部パラメータを推定し、その推定結果を出力する。
【００１８】
二次電池の充放電データから二次電池の内部パラメータを同時に推定する場合、各パラメータの値が収束しないおそれがあるが、本実施形態では、電池シミュレーションモデルを与える際に、内部パラメータの事前分布も与えるため、この内部パラメータ事前分布が制約条件ないし拘束条件として機能し、各パラメータの値が振動する等して収束しない事態を解消する。
【００１９】
二次電池の内部パラメータとしては、開回路電位（ＯＣＰ）Ｕｃ、イオン導電率κ、交換電流密度ｉｏ、イオン拡散定数Ｄｓ、初期イオン濃度ｃｓｏ、最大イオン濃度ｃｓｍａｘがあるが、本実施形態では内部パラメータを物理量に基づくパラメータと関数に基づくパラメータに分類し、それぞれに対して事前分布を与えることで適切な制約を課している。
【００２０】
図２に、本実施形態における内部パラメータ事前分布を具体的に示す。内部パラメータは、物理量に基づくパラメータと、関数に基づくパラメータに分類され、物理量に基づくパラメータのとり得る値の確率分布、及び計測結果に基づいて近似された関数のフィッティング度合いを事前分布として与える。物理量に基づくパラメータは、交換電流密度ｉｏ等であり、関数、すなわち計測結果に基づいて近似された関数に基づくパラメータは、開回路電位Ｕｃ等である。物理量に基づくパラメータのみならず、関数に基づくパラメータについても、計測結果に対する関数のフィッティング度合いを確率分布化して事前分布を与えている点に留意されたい。非特許文献２では、関数曲線を構成する点を直接推定しているため計算リソースを多く要するとともに、各点に対する拘束条件が課されていないため適切な値に収束する保証がないが、本実施形態では関数全体に拘束条件を課すことができるので計算リソースが削減されるとともに、適切な値に収束し得る。
【００２１】
図３に、推定装置１０における推定アルゴリズムを模式的に示す。推定装置１０は、二次電池の内部パラメータとして内部パラメータベクトルｘを設定し、電池シミュレーションモデルに与える。また、電池シミュレーションモデルの入力として、電流・温度パターンデータを与える。電池シミュレーションモデルは、内部パラメータベクトルｘを用いて、電流・温度パターンデータに対応する電圧データを推定して出力する。電池シミュレーションモデルは、公知の電気化学反応モデルに基づくシミュレーションモデルを用いることができ、例えば非特許文献２と同様にＤｏｙｌｅ−Ｆｕｌｌｅｒ−Ｎｅｗｍａｎ（ＤＦＮ）モデルを用いることができる。推定装置１０は、電池シミュレーションモデルにより推定された電圧データと、実電池で得られた実測の電圧データ（実測値）とを比較し、推定された電圧データを評価する。内部パラメータベクトルｘで電池シミュレーションモデルを用いて推定された電圧データをｇ（ｘ）、実測値をｙｔとすると、２乗誤差
ｆ（ｘ）＝Σ（ｙｔ−ｇ（ｘ））２
が最小になるように推定する。
【００２２】
パラメータベクトルｘの各要素は独立ではなく相互に関連し、例えば推定中にあるパラメータが大きくなると別のパラメータが小さくなり、その逆も生じるためパラメータが収束しない問題が生じ得る。この問題の多くはパラメータが想定外の値をとることにより生じる。本実施形態では、各パラメータに対する拘束条件を事前分布として導入するために、出力電圧パターンｙが得られた場合に、
ｐ（ｘ｜ｙ）∝ｐ（ｙ｜ｘ）ｐ（ｘ）
が最大になるように推定する。ここで、ｐ（ｘ｜ｙ）は、事象ｙが生じた後での、事象ｘの確率（事後確率あるいは条件付確率）であり、上式はベイズの定理として知られている。ｐ（ｘ）はパラメータベクトルｘの事前分布である。パラメータベクトルｘの事前分布を用いて拘束しつつ、事後確率ｐ（ｘ｜ｙ）が最大となるように内部パラメータベクトルｘを算出して推定結果を出力する。
【００２３】
以下、本実施形態における推定装置１０について、より詳しく説明する。
【００２４】
図４に、推定装置１０の構成ブロック図を示す。推定装置１０は、コンピュータから構成され、ＣＰＵ１２、プログラムメモリ１４、ワーキングメモリ１６、入出力インタフェースＩ／Ｆ１８、及びデータ記憶部２０を備える。なお、これら以外にも、キーボードやマウス等の入力装置、ディスプレイ等の出力装置を備えてもよい。
【００２５】
ＣＰＵ１２は、電池シミュレーションモデルを用いて電圧データを算出し、実測値と比較することで内部パラメータを評価する。ＣＰＵ１２は、単一である必要はなく、複数個のＣＰＵが協働して分散処理しても良く、あるいは専用回路で処理を実行しても良い。
【００２６】
プログラムメモリ１４には、電池シミュレーションモデルを記述するプログラム、及び内部パラメータの事前分布が記憶される。内部パラメータ事前分布は、物理量に基づくパラメータについてはそのとり得る値の確率分布であり、関数に基づくパラメータについてはその関数のフィッティング度合いである。シミュレーションモデルを記述するプログラムをプログラムメモリ１４に記憶し、内部パラメータの事前分布をデータ記憶部２０に記憶してもよい。
【００２７】
物理量に基づくパラメータは、二次電池をリチウムイオン二次電池とした場合、交換電流密度ｉｏ、イオン拡散定数Ｄｓ、初期イオン濃度ｃｓｏ、最大イオン濃度ｃｓｍａｘであり、事前分布はそれぞれ次のように設定してプログラムメモリ１４に記憶する。
【００２８】
＜交換電流密度ｉｏ＞
交換電流密度ｉｏの事前分布は、ガウス分布を仮定して、
【数１】
2015135286.tif 000003
と定義する。ここで、ｉｏｉｎｉは初期値、σは標準偏差であり、例えば１０％の誤差があるとして、σ＝ｉｏｉｎｉ×０．１とする。
【００２９】
＜イオン拡散定数Ｄｓ＞
イオン拡散係数Ｄｓの事前分布は、ガウス分布を仮定して、
【数２】
2015135286.tif 000004
と定義する。ここで、Ｄｓｉｎｉは初期値であり、標準偏差σは例えば１０％の誤差があるとして、σ＝Ｄｓｉｎｉ×０．１とする。
【００３０】
＜初期イオン濃度ｃｓｏ＞
初期イオン濃度ｃｓｏの事前分布は、ガウス分布を仮定して、
【数３】
2015135286.tif 000005
と定義する。ここで、ｃｓｏｉｎｉは初期値であり、標準偏差σは例えば１０％の誤差があるとして、σ＝ｃｓｏｉｎｉ×０．１とする。
【００３１】
＜最大イオン濃度ｃｓｍａｘ＞
最大イオン濃度ｃｓｍａｘの事前分布は、ガウス分布を仮定して、
【数４】
2015135286.tif 000006
と定義する。ここで、ｃｓｍａｘｉｎｉは初期値であり、標準偏差σは例えば１０％の誤差があるとして、σ＝ｃｓｍａｘｉｎｉ×０．１とする。
【００３２】
また、関数に基づくパラメータは、開回路電位ＯＣＰ、イオン導電率κであり、事前分布はそれぞれ次のように設定する。
【００３３】
＜開回路電位ＯＣＰＵｃ＞
正極の開回路電位曲線Ｕｃは、
【数５】
2015135286.tif 000007
の曲線でフィッティングする。ここで、係数ａｃ＝（ａ１・・・ａ１４）の初期値ａｃｏは、
【数６】
2015135286.tif 000008
と設定する。ａｃ自体は物理量ではなく、物理的意味を有しないため、事前分布は関数Ｕｃのフィッティング度合い
【数７】
2015135286.tif 000009
と定義する。ここで、Θ＝（０．１，０．２，０．３，・・・，０．９）とし、各局所ＳＯＣθｋでの標準偏差σｋは、例えば開回路電位曲線に１％の誤差があると想定して、σｋ＝Ｕｃ（θｋ｜ａｃｏ）×０．０１とする。
【００３４】
負極の開回路電位曲線Ｕａは、
【数８】
2015135286.tif 000010
の曲線でフィッティングする。ここで、係数ａａ＝（ａ１・・・ａ１４）の初期値ａａｏは、
【数９】
2015135286.tif 000011
と設定する。ａａ自体は物理量ではなく、物理的意味を有しないため、事前分布は関数Ｕａのフィッティング度合い
【数１０】
2015135286.tif 000012
と定義する。ここで、Θ＝（０．１，０．２，０．３，・・・，０．９）とし、各局所ＳＯＣθｋでの標準偏差σｋは、例えば開回路電位曲線に１％の誤差があると想定して、σｋ＝Ｕｃ（θｋ｜ａａｏ）×０．０１とする。
【００３５】
＜イオン導電率κ＞
イオン導電率κは、
【数１１】
2015135286.tif 000013
の曲線でフィッティングする。ここで、係数ｂ＝（ｂ１ｂ２・・・ｂ５）の初期値ｂｏは、
【数１２】
2015135286.tif 000014
と設定する。ｂ自体は物理量ではなく、物理的意味を有しないため、事前分布は関数κとのフィッティング度合い
【数１３】
2015135286.tif 000015
と定義する。ここで、ｃ＝（０．０００５，０．０００８，０．００１０，０．００１２，０．００１５，０．００２０）とし、各イオン濃度ｃκでの標準偏差σκは、例えばイオン導電率に１％の誤差があるとしてσκ＝κ（ｃ｜ｂｏ）×０．０１とする。
【００３６】
入出力Ｉ／Ｆ１８は、充放電データを入力してデータ記憶部２０に記憶する。充放電データは、図３に示すように電流・温度パターンデータであるが、本実施形態では、温度は２５度で固定し、電流パターンを供給する。具体的には、例えば０．１Ｃ，０．５Ｃ，１．０Ｃ，２．０Ｃの４通りの定電流でＳＯＣ１００％から０％まで放電し、逆に、０％から１００％まで充電したときの電流値の時系列データを供給する。
【００３７】
図５に、推定装置１０における推定処理のフローチャートを示す。遺伝的アルゴリズム（Genetic Algorithm）を用いた処理である。
【００３８】
まず、内部パラメータの初期値を設定する（S１０１）。初期値は、例えば次のように設定する。
＜交換電流密度ｉｏ＞
正極ｉｏｃｉｎｉ＝０．０８ｍＡ／ｃｍ２
負極ｉｏａｉｎｉ＝０．１１ｍＡ／ｃｍ２
＜イオン拡散係数Ｄｓ＞
正極Ｄｓｃｉｎｉ＝１．０×１０−９ｃｍ２／ｓ
負極Ｄｓａｉｎｉ＝３．９×１０−１０ｃｍ２／ｓ
＜初期イオン濃度ｃｓｏ＞
正極ｃｓｏｃｉｎｉ＝３．９０（ｍｏｌ／ｄｍ３）
負極ｃｓｏａｉｎｉ＝１４．８７（ｍｏｌ／ｄｍ３）
＜最大イオン濃度ｃｓｍａｘ＞
正極ｃｓｍａｘｃｉｎｉ＝２２．８６（ｍｏｌ／ｄｍ３）
負極ｃｓｍａｘａｉｎｉ＝２６．３９（ｍｏｌ／ｄｍ３）
＜開回路電位Ｕｃ＞
正極の係数ａｃの初期値ａｃｏ及び負極の係数ａａの初期値ａａｏは上記した通りである。
＜イオン導電率κ＞
係数ｂの初期値ｂｏは上記した通りである。
【００３９】
設定された内部パラメータの初期値は、入出力Ｉ／Ｆ１８を介してプログラムメモリ１４あるいはデータ記憶部２０に記憶される。
【００４０】
次に、複数仮説を作成する（Ｓ１０２）。具体的には、ＣＰＵ１２は、プログラムメモリ１４あるいはデータ記憶部２０に記憶されている各内部パラメータの初期値を読み出し、これらの初期値にノイズを加えて複数、例えば１０００個の仮説を作成する。交換電流密度ｉｏの場合、
ｉｏｃ＝ｉｏｃ＋０．０１×ｎｒａｎｄ×ｉｏｃ
等である。ここで、ｎｒａｎｄは平均０、分散１の正規分布に従う乱数である。
【００４１】
次に、作成した複数の仮説それぞれをパラメータとして、充放電データと同様の電流値、すなわち０．１Ｃ、０．５Ｃ、１Ｃ、２Ｃに対する電圧値を電池シミュレーションモデルを用いて算出し、以下の式でスコアｐ（ｙ｜ｘ）を算出する。
【数１４】
2015135286.tif 000016
ここで、ｇ（ｘ）は電池シミュレーションモデルを用いて算出された時系列の電圧値であり、ｙｔはデータ記憶部２０に記憶されている実測値の電圧値である。また、σは実測値の標準偏差である。
【００４２】
次に、ＣＰＵ１２は、各内部パラメータの事前分布をプログラムメモリ１４あるいはデータ記憶部２０から読み出し、各仮説のパラメータの事前分布ｐ（ｘ）を、それぞれのパラメータの事前分布を用いて以下のように算出する。
【数１５】
2015135286.tif 000017
【００４３】
この式は、各パラメータの事前分布が互いに独立と仮定したものである。
【００４４】
そして、これらの算出結果より、仮説毎に事後確率ｐ（ｘ｜ｙ）に比例するスコアを算出することで仮説を評価する（Ｓ１０３）。スコアは、ベイズの定理に従い、
ｓｃｏｒｅｉ＝ｐ（ｙ｜ｘ）ｐ（ｘ）
により算出する。このスコアｓｃｏｒｅｉが最大となる内部パラメータベクトルｘが実測値に最も適合する内部パラメータベクトルと推定される。
【００４５】
次に、ＣＰＵ１２は、遺伝的アルゴリズムにおける交叉、複製、突然変異の確率をそれぞれ０．７，０．２，０．１として動作を選択する（Ｓ１０４）。また、ＣＰＵ１２は、各仮説を、
【数１６】
2015135286.tif 000018
の確率で、選択するか否かを決定する（ルーレット選択）。そして、ＣＰＵ１２は、選択した交叉、複製、突然変異のいずれかの動作を、上記の確率で選択した仮説に対して実行する。仮設の数が所定数、例えば１０００個になるまで動作の選択及び仮説の選択を繰り返す。
【００４６】
最後に、ＣＰＵ１２は、世代数（繰り返し回数）が所定の上限値、例えば１００００に達したか、あるいは前世代と現世代の仮説の平均ベクトルの２乗誤差が０．００００００１以下になったか否かを判定し（Ｓ１０５）、条件が満たされた場合に繰り返し処理を終了する。前世代と現世代の仮説の平均ベクトルの２乗誤差が０．００００００１以下とは、要するに前回と今回の推定結果に変化がないことを意味する。ＣＰＵ１２は、最終的に得られた内部パラメータベクトルｘを最終的な推定結果として入出力Ｉ／Ｆ１８を介して外部機器、例えばディスプレイに出力する。
【００４７】
図６に、交換電流密度ｉｏの事前分布の一例を示す。事前分布は、図示のようにガウス分布で与えられる。また、図７に、開回路電位Ｕｃの事前分布の一例を示す。事前分布は、関数のフィッティング度合いで与えられる。
【００４８】
このように、本実施形態では、二次電池の内部パラメータを物理量に基づくパラメータと関数に基づくパラメータに分類し、前者に関してはガウス分布として事前分布を規定し、後者に関しては関数のフィッティング度合いとして事前分布を規定することで制約条件を課すことで、それぞれのパラメータが想定外の値に収束してしまう事態を確実に防止できる。
【００４９】
以上、本発明の実施形態について説明したが、本発明はこれに限定されず、種々の変形が可能である。
【００５０】
例えば、本実施形態では、物理量に基づくパラメータの事前分布としてガウス分布を設定したが、ガウス分布に限らず他の確率分布関数を用いても良い。
【００５１】
また、本実施形態では、内部パラメータの推定に遺伝的アルゴリズムを用いているが、遺伝的アルゴリズムに限らず他のアルゴリズムを用いても良い。
【００５２】
また、本実施形態では、二次電池の内部パラメータとして、開回路電位、イオン導電率、交換電流密度、イオン拡散定数、初期イオン濃度、最大イオン濃度を示したが、必ずしもこれら全てを推定する必要はなく、いずれか一部、例えばイオン拡散定数とイオン導電率と開回路電位を推定する構成としてもよい。
【００５３】
さらに、本実施形態では、二次電池としてリチウムイオン二次電池を例示したが、これに限定されるものではなく任意の二次電池に適用し得る。リチウムイオン二次電池は、具体的には負極にＬｉｘＣ６、正極にＬｉｙＭｎ２Ｏ４、電解液にエチレンカーボネートＥＣとジメチルカーボネートＤＭＣ、混合体積比ＥＣ：ＤＭＣ＝１：２としたものであるが、これに限定されない。初期値に関しては、計測して求めることができるが、例えば交換電流密度に関しては交流インピーダンス法で求めた内部抵抗値から算出することができ、イオン拡散定数についてはＮＭＲ法等を用いることができる。
【符号の説明】
【００５４】
１０  推定装置、１２  ＣＰＵ、１４  プログラムメモリ、１６  ワーキングメモリ、１８  入出力Ｉ／Ｆ、２０  データ記憶部。
【図１】
2015135286.tif 000020
【図２】
2015135286.tif 000021
【図３】
2015135286.tif 000022
【図４】
2015135286.tif 000023
【図５】
2015135286.tif 000024
【図６】
2015135286.tif 000025
【図７】
2015135286.tif 000026