(57)【要約】【課題】  適応デジタルフィルタを用いた二次電池の内部状態を正確に推定する。【解決手段】  電池ＥＣＵは、ＦＦＴ等を用いて、二次電池の充放電電流が計算精度の低い波形であるか否かを判断するステップ（Ｓ１００）と、二次電池の充放電電流が計算精度の低い波形であると（Ｓ１００にてＹＥＳ）、モータＥＣＵに駆動指令を出力して、二次電池の充放電電流を微小電流分だけ変化させるステップ（Ｓ２００）とを含む、プログラムを実行する。【選択図】        図７
2006105823.tif 000002
【特許請求の範囲】
【請求項１】
電池モデルを用いた適応デジタルフィルタ演算で開路電圧を推定して、その開路電圧から二次電池の内部状態を推定する装置であって、
前記二次電池の電流値を検知するための検知手段と、
前記検知された電流値が入力され、パラメータを逐次推定する、電池モデルを用いた適応デジタルフィルタと、
前記適応デジタルフィルタにより推定されたパラメータに基づいて、開路電圧推定値を算出するための算出手段と、
予め定められた関係と前記開路電圧推定値とに基づいて、前記電池内部状態を推定するための推定手段と、
前記電流値が予め定められた条件を満足すると、電流値に変動を与えるための変動付与手段とを含む、推定装置。
【請求項２】
前記条件は、前記電流値の波形が正弦波であるという、請求項１に記載の推定装置。
【請求項３】
前記条件は、前記電流値の波形が単調増加であるという、請求項１に記載の推定装置。
【請求項４】
前記推定手段は、前記電池の残存容量を推定するための手段を含む、請求項１〜３のいずれかに記載の推定装置。
【請求項５】
前記推定手段は、前記電池に入出力可能な電力値を推定するための手段を含む、請求項１〜３のいずれかに記載の推定装置。
【請求項６】
電池モデルを用いた適応デジタルフィルタ演算で開路電圧を推定して、その開路電圧から二次電池の内部状態を推定する方法であって、
前記二次電池の電流値を検知するステップと、
パラメータを逐次推定する、電池モデルを用いた適応デジタルフィルタに、前記検知された電流値を入力するステップと、
前記適応デジタルフィルタにより推定されたパラメータに基づいて、開路電圧推定値を算出するステップと、
予め定められた関係と前記開路電圧推定値とに基づいて、前記電池内部状態を推定する推定ステップと、
前記電流値が予め定められた条件を満足すると、電流値に変動を与えるステップとを含む、推定方法。
【請求項７】
前記条件は、前記電流値の波形が正弦波であるという、請求項６に記載の推定方法。
【請求項８】
前記条件は、前記電流値の波形が単調増加であるという、請求項６に記載の推定方法。
【請求項９】
前記推定ステップは、前記電池の残存容量を推定するステップを含む、請求項６〜８のいずれかに記載の推定方法。
【請求項１０】
前記推定ステップは、前記電池に入出力可能な電力値を推定するステップを含む、請求項６〜８のいずれかに記載の推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、二次電池の残存容量（ＳＯＣ（States Of Charge））を含む二次電池の内部状態を算出する技術に関し、特に、適応デジタルフィルタを用いて電池の内部状態を推定する技術に関する。
【背景技術】
【０００２】
エンジン（たとえば、ガソリンエンジンやディーゼルエンジンなどの公知の機関を用いることが考えられる。）と電気モータとを組合せたハイブリッドシステムと呼ばれるパワートレインを搭載した車両が開発され、実用化されている。このような車両においては、運転者のアクセル操作量に関係なく、エンジンによる運転と電気モータとによる運転とが自動的に切換えられて、最も効率が良くなるように制御される。たとえば、エンジンが、定常状態で運転されて二次電池（以下、バッテリや電池と記載するが全て同義である）を充電する発電機を回すために運転される場合、あるいは二次電池のＳＯＣなどに応じて走行中に間欠的に運転される場合などは、運転者によるアクセルの操作量とは無関係にエンジンの運転および停止を繰返す。つまりエンジンと電気モータとをそれぞれ単独、または協同して動作させることにより、燃料消費向上や排気ガスを大幅に抑制することが可能になる。
【０００３】
このように、ハイブリッド車両のエンジンは、走行中においても間欠駆動が行なわれることになり、頻繁に停止制御が行なわれるようになる。エンジンの停止制御を行なう場合、エンジンの回転エネルギ（運動エネルギ）をジェネレータに与えて、ジェネレータにより電気エネルギに変換して、車両を停止させることがある。このときに発電された電力は二次電池を充電するために用いられる。この二次電池のＳＯＣが高いときに充電を許可すると過充電となり二次電池の温度上昇を招き、二次電池の寿命を縮めることになりかねない。そこで、二次電池のＳＯＣをできる限り正確に推定することが求められる。また、制動時に回生制動を行なう時にも、回収した電力を二次電池に充電できるか否かを判断するためにもＳＯＣをできる限り正確に推定することが求められる。このようなことは、走行用の二次電池を有する、ハイブリッド車両、電気自動車、燃料電池車のいずれにも共通する。
【０００４】
非特許文献１は、二次電池の内部状態量（ＳＯＣを含む）を電池の端子電圧、電流温度などを用いて推定する際における、電池の過渡特性の影響による推定精度の悪化を改善するために適応デジタルフィルタを適用した、電池内部状態量の推定手法を開示する。
【０００５】
この電池内部状態量の推定手法は、リチウムイオン電池などの二次電池の電流と端子電圧との関係を、線形パラメータ可変（ＬＰＶ（Linear Parameter Varying））モデルで記述することにより、「電池内部状態量の推定手法」を、「適応デジタルフィルタによるモデルパラメータ逐次推定問題」に帰着させている。特に、物理的な意味を持つ連続時間形式電池モデルのパラメータを、適応デジタルフィルタで直接同定する「ハイブリッド同定方式」を適用し、また、常に変化するモデルパラメータを精度良く逐次推定する必要性から、パラメータ同定方法として特に「両限トレースゲイン型最小二乗法」を適用している。このようにして、動作条件（ＳＯＣ、温度、劣化度など）によって常に変動する電池モデルパラメータ（内部抵抗、時定数など）と、電池内部状態量の一例である、開路電圧（開放電圧（ＯＣＶ：Open Circuit Voltage）と同じ）、ＳＯＣ、入出力可能パワーを精度良く逐次推定できる推定システムを実現している。
【非特許文献１】「１３８  適応デジタルフィルタ理論を用いた電池内部状態の推定手法」、自動車技術会、学術講演会前刷集、Ｎｏ．３０−０３、２００３
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、上述した文献に開示された電池内部状態量の推定手法においては、計測された電圧から適応デジタルフィルタを用いて開路電圧を算出し、開路電圧からＳＯＣを算出しているが、電池の電流波形に依存して推定精度が変動するという問題点を有する。たとえば、その理由は明確ではないが、電池の電流波形が正弦波や単調増加電流である場合には、電池の内部状態量の推定精度が悪化する。
【０００７】
本発明は、上述の課題を解決するためになされたものであって、その目的は、このような推定精度の悪化を回避することができる、適応デジタルフィルタを用いた二次電池の内部状態を推定する装置およびその方法を提供することである。
【課題を解決するための手段】
【０００８】
第１の発明に係る推定装置は、電池モデルを用いた適応デジタルフィルタ演算で開路電圧を推定して、その開路電圧から二次電池の内部状態を推定する。この推定装置は、二次電池の電流値を検知するための検知手段と、検知された電流値が入力され、パラメータを逐次推定する、電池モデルを用いた適応デジタルフィルタと、適応デジタルフィルタにより推定されたパラメータに基づいて、開路電圧推定値を算出するための算出手段と、予め定められた関係と開路電圧推定値とに基づいて、電池内部状態を推定するための推定手段と、電流値が予め定められた条件を満足すると、電流値に変動を与えるための変動付与手段とを含む。
【０００９】
また、第６の発明に係る推定方法は、電池モデルを用いた適応デジタルフィルタ演算で開路電圧を推定して、その開路電圧から二次電池の内部状態を推定する。この推定方法は、二次電池の電流値を検知するステップと、パラメータを逐次推定する、電池モデルを用いた適応デジタルフィルタに、検知された電流値を入力するステップと、適応デジタルフィルタにより推定されたパラメータに基づいて、開路電圧推定値を算出するステップと、予め定められた関係と開路電圧推定値とに基づいて、電池内部状態を推定する推定ステップと、電流値が予め定められた条件を満足すると、電流値に変動を与えるステップとを含む。
【００１０】
これらの発明によると、二次電池の電流値が予め定められた、推定精度を悪化させる波形等であるという条件を満足すると、電流値に変動を与える。この与えられた変動によりたとえば、二次電池の開路電圧推定値の精度が向上して、この開路電圧を用いて算出されるＳＯＣや入出力可能パワー等の内部状態を正確に算出することができる。その結果、このような推定精度の悪化を回避することができる、適応デジタルフィルタを用いた二次電池の内部状態を推定する装置およびその方法を提供することができる。
【００１１】
第２の発明に係る推定装置においては、第１の発明の構成に加えて、条件は、電流値の波形が正弦波であり、第７の発明に係る推定方法においては、第６の発明の構成に加えて、条件は、電流値の波形が正弦波であるというものである。
【００１２】
これらの発明によると、二次電池の電流波形が正弦波であるときには、開路電圧等の二次電池の内部状態の推定値の精度が悪化するので、電流値に変動を与える。この与えられた変動により、たとえば、二次電池の開路電圧推定値の精度が向上する。
【００１３】
第３の発明に係る推定装置においては、第１の発明の構成に加えて、条件は、電流値の波形が単調増加であり、第８の発明に係る推定方法においては、第６の発明の構成に加えて、条件は、電流値の波形が単調増加であるというものである。
【００１４】
これらの発明によると、二次電池の電流波形が単調増加しているときには、開路電圧等の二次電池の内部状態の推定値の精度が悪化するので、電流値に変動を与える。この与えられた変動により、たとえば、二次電池の開路電圧推定値の精度が向上する。
【００１５】
第４の発明に係る推定装置においては、第１〜３のいずれかの発明の構成に加えて、推定手段は、電池の残存容量を推定するための手段を含み、第９の発明に係る推定方法においては、第６〜８のいずれかの発明の構成に加えて、推定ステップは、電池の残存容量を推定するステップを含む。
【００１６】
これらの発明によると、二次電池の電流波形が推定精度を悪化させる場合には、開路電圧等の二次電池の内部状態の推定値の精度が悪化するので、電流値に変動を与え、開路電圧推定値を算出する。この与えられた変動により、二次電池の開路電圧推定値の精度が向上するので、これを用いて算出する電池の残存容量（ＳＯＣ）の推定精度も向上する。
【００１７】
第５の発明に係る推定装置においては、第１〜３のいずれかの発明の構成に加えて、推定手段は、電池に入出力可能な電力値を推定するための手段を含み、第１０の発明に係る推定方法においては、第６〜８のいずれかの発明の構成に加えて、推定ステップは、電池に入出力可能な電力値を推定するステップを含む。
【００１８】
これらの発明によると、二次電池の電流波形が推定精度を悪化させる場合には、開路電圧等の二次電池の内部状態の推定値の精度が悪化するので、電流値に変動を与え、開路電圧推定値を算出する。この与えられた変動により、二次電池の開路電圧推定値の精度が向上するので、これを用いて算出する電池に入出力可能な電力値の推定精度も向上する。
【発明を実施するための最良の形態】
【００１９】
以下、図面を参照しつつ、本発明の実施の形態について説明する。以下の説明では、同一の部品には同一の符号を付してある。それらの名称および機能も同じである。したがってそれらについての詳細な説明は繰返さない。
【００２０】
まず、図１〜５を参照して本発明の前提となる適応デジタルフィルタについて説明する。
【００２１】
図１は、未知プラントが以下の式（１）で表わされるモデルである。
【００２２】
【数１】
2006105823.tif 000003
【００２３】
式（１）において、ｙは計測可能なプラント出力、ωは計測可能なプラント入力ｕと、入力ｕと出力ｙをフィルタに通した出力からなるベクトル、θは計測不可能なパラメータベクトルである。これらの信号は離散時間または連続時間信号である。
【００２４】
プラント出力ｙと、モデルを用いたプラント出力の推定値ｙ'との誤差が０に近づくように、式（１）で示すモデルのパラメータ推定値ベクトルθ'を、最小二乗法などを用いて逐次的に調整する方法が「適応デジタルフィルタ（逐次型モデルパラメータ同定）である。計測可能なプラント入力ｕ、出力ｙ、プラントモデルの構造情報（次数など）だけで、プラントの特性変化に容易に対応できる。
【００２５】
二次電池の電流と端子電圧との関係を、電池構造をもとに図２に示す簡易モデルで表現する。充放電電流Ｉ（ｔ）、開路電圧ＶO（ｔ）と端子電圧Ｖ（ｔ）との関係は式（２）または式（３−１）で表わされる。
【００２６】
【数２】
2006105823.tif 000004
【００２７】
ここで、ｓは微分オペレータ、Ｋは内部抵抗、Ｔ1およびＴ2は時定数であり、式（３−２）のように与えられる。
【００２８】
【数３】
2006105823.tif 000005
【００２９】
さらに、式中の開路電圧は、充放電電流の積算値に相関が深く、実験結果から式（４）で近似する。ただし、ｈは開路電圧用定数である。
【００３０】
【数４】
2006105823.tif 000006
【００３１】
式（４）を式（３）に代入して式（５）が得られる。
【００３２】
【数５】
2006105823.tif 000007
【００３３】
以上のように、電流を入力、端子電圧を出力として電池をモデル化している。なお、電池内部で起こる現象を詳細に記述したモデルではなく、特定の動作条件（ＳＯＣ，温度，劣化度など）での入出力関係を線形モデルで近似したものである。したがって各モデルパラメータ（内部抵抗Ｋ，時定数Ｔ1，Ｔ2，開路電圧用定数ｈ）は固定値ではなく、動作条件の変化に応じて常に変化することが前提である。つまり、線形パラメータ可変（ＬＰＶ）システムとして定義することで問題を簡略化している。上記モデルは離散化すれば、当然回帰型モデルとなる。
【００３４】
図３に、開路電圧とＳＯＣとの定常的な相関特性を示す。一般に、この特性は電池固有の特性であって、温度や劣化度には影響されないとされている。
【００３５】
連続時間形式で記述された式（５）の電池モデルを、離散化せず、適応したデジタルフィルタ（逐次型モデルパラメータ同定）にダイレクトに帰着させるハイブリッド同定方式を適用する。これは、１）連続時間系モデルのパラメータが物理的な意味をもち、推定精度が把握しやすいこと、２）離散時間同定法の方が一般的にパラメータの収束性がよいこと、３）さらに同定のための信号として入出力信号のフィルタ出力を用いるための観測ノイズの影響を緩和できるからである。
【００３６】
以下、適応デジタルフィルタが適用可能なモデル標準形への定式化を示す。まず、式（５）を式（６）のように書き換える。
【００３７】
【数６】
2006105823.tif 000008
【００３８】
ここで、ローパスフィルタ１／Ｇlp（ｓ）は、式（６）の両辺が厳密にプロパーとなるように式（７）のように選ぶ（τは時定数）。
【００３９】
【数７】
2006105823.tif 000009
【００４０】
電流および電圧のフィルタ出力を式（８）で定義する。
【００４１】
【数８】
2006105823.tif 000010
【００４２】
このように定義すると、式（６）は式（９）のように書ける。
【００４３】
【数９】
2006105823.tif 000011
【００４４】
式（９）は、適応デジタルフィルタのモデル標準形であり適応デジタルフィルタが適用可能である。なお、θは未知パラメータであり、それ以外は計測不可能な信号をもとに式（８）のフィルタ処理で得られる値である。
【００４５】
パラメータ推定値をθ'、パラメータ推定値θ'を用いて推定される出力推定値をｙ'とすると、式（１０）に示す誤差方程式が得られる。
【００４６】
【数１０】
2006105823.tif 000012
【００４７】
式（１０）の誤差方程式に対して、逐次型のパラメータ同定アルゴリズムを適用する。特に、変動する未知パラメータθを精度良く推定するために、最小二乗法を発展させた両限トレースゲイン方式を用いた、離散時間形式で記述されたアルゴリズムを式（１１）に示す。ただし、λ3，α1は、０＜λ3＜１，０＜α1＜∞の範囲で設定する適応推定の重みであり、γU，γLは、適応推定の上下限を設定する値である。
【００４８】
【数１１】
2006105823.tif 000013
【００４９】
適応デジタルフィルタにより逐次推定されたパラメータθ'（ｔ）から開路電圧推定値Ｖ'0（ｔ）を算出する。ただし、適応デジタルフィルタで用いた式（４）から算出すると、誤差が積分的に溜まるので、式（４）ではなく、積分器を含まない式（３）の電池モデルを直接用いて開路電圧推定値Ｖ'0（ｔ）を算出する。前述の処理と同様に、電池モデルを表わす式（３）の分母多項式を払い、両辺に式（７）のローパスフィルタを乗じて整理すると式（１２）が得られる。
【００５０】
【数１２】
2006105823.tif 000014
【００５１】
式（１２）の左辺は、開路電圧Ｖ0（ｔ）にローパスフィルタ処理を施した値であるが、開路電圧Ｖ0（ｔ）は、変化が比較的緩やかであるため、これを開路電圧推定値Ｖ'0（ｔ）とみなす。さらにＫ等のパラメータを適応デジタルフィルタで算出される連続時間系モデルの推定パラメータθ'（ｔ）すなわち、Ｋ'（ｔ），Ｔ1'（ｔ），Ｔ2'（ｔ）に置き換えて、式（８）に示すフィルタ処理後の電流Ｉ1（ｔ），Ｉ2（ｔ）および端子電圧Ｖ1（ｔ），Ｖ2（ｔ）を式（１２）に代入して、式（１３）のように開路電圧推定値Ｖ'0（ｔ）を算出する。なお、同定機構を図４に示す。
【００５２】
【数１３】
2006105823.tif 000015
【００５３】
適応デジタルフィルタにより逐次推定された電池モデルパラメータθ'（ｔ）と、既に算出された内部状態量である開路電圧推定値Ｖ'0（ｔ）から、充電率（ＳＯＣ）推定値Ｓ'OC（ｔ），入出力可能パワー推定値Ｐ'IN（ｔ），Ｐ'OUT（ｔ）を算出する。充電率推定値Ｓ'OC（ｔ）は、予め計測した開路電圧ＶO（ｔ）とＳＯＣとの相関特性マップ（図３）を用いて、逐次推定された開路電圧推定値Ｖ'0（ｔ）から換算する。入出力パワーは、逐次推定された開路電圧推定値Ｖ'0（ｔ）と内部抵抗推定値Ｋ'（ｔ）から電流Ｉ（ｔ）と端子電圧Ｖ（ｔ）の定常特性である式（１４）を用いて算出する。
【００５４】
【数１４】
2006105823.tif 000016
【００５５】
入力可能パワー推定値Ｐ'IN（ｔ）は、電池の上限端子電圧ＶMAXに至る限界値と定義すれば式（１５）で算出できる。
【００５６】
【数１５】
2006105823.tif 000017
【００５７】
出力可能パワー推定値Ｐ'OUT（ｔ）は、電池の下限端子電圧ＶMINに至る限界値と定義すれば式（１６）で算出できる。
【００５８】
【数１６】
2006105823.tif 000018
【００５９】
内部状態量の算出方法を図示すると図５のようになる。本実施の形態においては、サンプリング周期を５０ｍｓ、演算周期を２００ｍｓとした。このような構成を有する本実施の形態に係る二次電池内部状態推定装置が適用される車両について図６を用いて説明する。
【００６０】
この車両は、エンジン（ＥＮＧ）とモータとを車両の走行源として搭載する。エンジン（ＥＮＧ）から変速機（ＡＭＴ）を介して駆動輪（たとえば前輪）を駆動するとともに、バッテリからインバータ（ＩＮＶ２）を介してモータに電力が供給され、モータから発生するトルクによりエンジン（ＥＮＧ）をアシストする。また、エンジン（ＥＮＧ）により３相発電機が駆動されて発電された電力がインバータ（ＩＮＶ１）または整流器を介して直流変換されて、その変換された電力によりバッテリが充電される。
【００６１】
なお、このような構成を有する車両に限定して、本発明の実施の形態にかかる二次電池内部状態推定装置が適用されるものではない。たとえば、このようなハイブリッド車両ではなく、エンジンを有さない車両（電気自動車）であっても、エネルギ源としてのエンジンの代わりに燃料電池を有する車両であってもよい。
【００６２】
このような構成を有する本実施の形態に係る二次電池内部状態推定装置における、精度向上を実現するプログラムの制御構造についてフローチャートを用いて説明する。なお、二次電池内部状態推定装置は、たとえば電池ＥＣＵ（Electronic Control Unit）により実現される。
【００６３】
図７を参照して、ステップＳ１００にて、電池ＥＣＵは、計算精度の悪い電流波形であるか否かを判断する。バッテリの内部状態の推定開始後に全く電流が流れないときを含めて、特定の正弦波の波形を有する電流、特にオフセットのない正弦波形状の負荷電流の場合に誤差が大きくなる傾向を有する。これを以下のようにして検知する。たとえば、
１）バッテリ電流に対してＦＦＴ（Fast Fourier Transform）処理を行ない、周波数分散の急峻度により、
２）ＦＦＴで最大振幅の基本波振幅を抽出し、その振幅２乗の時間平均と信号全体の２乗値時間平均と比率（歪率）により、
３）特定の周波数の正弦波電流が大きい場合（ＦＦＴ処理）は、時定数の長い１次遅れフィルターを通して、オフセット成分を抽出して、抽出されたオフセットの大きさにより、
計算精度の悪い電流波形であるか否かを判断する。
【００６４】
なお、これらの判断手法は一例である。計算精度の悪い電流波形であると（ステップＳ１００にてＹＥＳ）、処理はステップＳ２００へ移される。もしそうでないと（ステップＳ１００にてＮＯ）、特に何も処理されないで、この処理は終了する。すなわち、計算精度の悪い電流波形でないので、推定結果に問題を生じないからである。
【００６５】
ステップＳ２００にて、電池ＥＣＵは、微小電流変化処理を実行する。この処理においては、
１）バッテリに対する負荷であるモータに対して、ユーザの違和感が発生しない範囲でトルク指令を変動させたり（ゼロトルク電流指令を電池ＥＣＵがモータＥＣＵに送信して、モータ側でトルク変動レスになるように処理しても良い）
２）エアコンディショナ作動時であればファンのオンオフ制御を切り換えたり、デフォッガ作動時であればそのオンオフ制御を切り換えるようにしたり
３）オルタネータの発電電力をオンオフ制御を切り換えたり（またはＨｉ／Ｌｏ切り換え）
することによりバッテリの電流が微小変化される。
【００６６】
以上のような構造およびフローチャートを有する二次電池内部状態推定装置の動作について説明する。
【００６７】
バッテリの電流を監視して、たとえばＦＦＴ処理の結果、計算精度が悪い波形の電流波形であると（ステップＳ１００にてＹＥＳ）、微小電流変化処理が実行される（ステップＳ２００）。微小電流変化処理が実行されると、バッテリの電流波形が変化して、計算結果が向上する。
【００６８】
図８に本実施の形態に係る二次電池内部状態推定装置を用いた場合（すなわち、微小電流変化処理を行なった場合）、図９に従来の推定装置を用いた場合（すなわち、微小電流変化処理を行なわない場合）の、それぞれの時間変化を示す。なお、図８および図９においては、１００秒後に突然バッテリの温度変化などでモデル（係数、パラメータ）が変化した場合の同定性能を評価するための図である。１００秒よりも前が「Model１」であって、それより後が「Model２」である。
【００６９】
図９に示すように、従来の推定装置においては、３段目左欄の適応推定誤差の部分に表わされるように、出力Ｙの推定誤差（ｙ−ｙ＿ｈ）が０であっても、右欄に示すように、推定結果と真値とは一致しない（一致しているのは２段目の「Model１」のみである）。なお、太線が推定演算した開路電圧ＯＣＶ（式（１３）のＶ'O（ｔ））を示し、細線がe(t)＝y-y_hを一次遅れ処理（時定数：１０秒）した演算結果である「適応推定誤差」を示す。
【００７０】
ところが、図８に示すように、微小電流変化処理を行なった場合には、３段目左欄の適応推定誤差の部分に表わされるように、出力Ｙの推定誤差（ｙ−ｙ＿ｈ）が０になると、右欄に示すように、推定結果と真値とは１〜３段目の全てにおいて「Model１」および「Model２」で一致している。すなわち、バッテリ電流が正弦波（０．１Ｈｚ）の場合に、方形波を追加することで微小電流変化処理を行なうと、推定結果と真値とが一致するのである。
【００７１】
なお、開路電圧ＯＣＶの推定が向上することにより、ＳＯＣも入出力可能パワーも精度良く推定することができるようになる（式（１５）および式（１６）参照）。また、二次電池の内部状態量の推定値の精度が低下するのは、電流波形が正弦波である場合に限定されない。この他にも単調増加波形や、オフセットを有する波形等の場合に、推定精度が悪化する。
【００７２】
以上のようにして、本実施の形態に係る二次電池内部状態推定装置によると、二次電池の電流値が予め定められた、推定精度を悪化させる波形等である、たとえば正弦波であると、電流値に変動を与える。この与えられた変動により、二次電池の開路電圧推定値の精度が向上して、この開路電圧を用いて算出されるＳＯＣや入出力可能パワー等の内部状態を正確に算出することができる。
【００７３】
今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【図面の簡単な説明】
【００７４】
【図１】適応デジタルフィルタの概略構成図である。
【図２】電池ＬＰＶモデルを示す図である。
【図３】開路電圧とＳＯＣとの関係を示す図である。
【図４】パラメータ同定機構の概略構成図である。
【図５】内部状態量の算出機構の概略構成図である。
【図６】二次電池内部状態推定装置が適用される車両のモデルである。
【図７】本実施の形態に係る二次電池内部状態推定装置において実行されるプログラムの制御構造を示すフローチャートである。
【図８】本実施の形態に係る二次電池内部状態推定装置の効果を説明するための図（その１）である。
【図９】本実施の形態に係る二次電池内部状態推定装置の効果を説明するための図（その２）である。
【図１】
2006105823.tif 000019
【図２】
2006105823.tif 000020
【図３】
2006105823.tif 000021
【図４】
2006105823.tif 000022
【図５】
2006105823.tif 000023
【図６】
2006105823.tif 000024
【図７】
2006105823.tif 000025
【図８】
2006105823.tif 000026
【図９】
2006105823.tif 000027