JP5036733B2 - 推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステム及び方法 - Google Patents

推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステム及び方法 Download PDF

Info

Publication number
JP5036733B2
JP5036733B2 JP2008557199A JP2008557199A JP5036733B2 JP 5036733 B2 JP5036733 B2 JP 5036733B2 JP 2008557199 A JP2008557199 A JP 2008557199A JP 2008557199 A JP2008557199 A JP 2008557199A JP 5036733 B2 JP5036733 B2 JP 5036733B2
Authority
JP
Japan
Prior art keywords
battery
vector
estimated
predicted
parameter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2008557199A
Other languages
English (en)
Other versions
JP2009528662A (ja
Inventor
グレゴリー、エル.プレット
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
LG Chem Ltd
Original Assignee
LG Chem Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by LG Chem Ltd filed Critical LG Chem Ltd
Publication of JP2009528662A publication Critical patent/JP2009528662A/ja
Application granted granted Critical
Publication of JP5036733B2 publication Critical patent/JP5036733B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables

Description

本発明は、デジタルフィルタリング(digital filtering)技術を活用して推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステム及び方法に関する。
再充電可能なバッテリーパックの技術分野において、直接的には測定できないが、現在のバッテリーパック状態を示すことができる量的要素(quantities)が推定できることがある応用例で要求される。このような量的要素の中で、一部は上記パックの充電状態(SOC,State Of Charge)のように急速に変化され得、全体の領域範囲を数分内に旋回することもできる。他の一部はセル容量(cell capacity)のように非常に徐々に変化され得、10年以上の通常の使用により20%程度の小さい範囲内で変化され得る。急速に変化する傾向のある量的要素はシステムの状態を構成し、遅く変化する傾向のある量的要素はシステムの時変(time varying)パラメータを構成する。
バッテリーシステムの分野、特に例えば、ハイブリッド電気輸送手段(Hybrid Electric Vehicle:HEV)、バッテリー電気輸送手段(Battery Electric Vehicles:BEV)、ラップトップコンピュータのバッテリー、移動可能な装備のバッテリーパックなどのように、バッテリーの生命に悪影響を及ぼさずにできれば活発に長い時間作動する必要があるものにおいて、どれくらいのバッテリーエネルギーが現在作動できるかなどを推定するのに、急変するパラメータ(例えば、SOC)に対する情報が用いられることが望ましい。さらに、有用なサービス時間を拡張しながらパックの寿命期間に先行演算(prior calculations)を正確に維持し、パックの健康状態(state of health:SOH)を判断するのに活用するために、緩変(slowly varying)するパラメータ(例えば、全体容量)に対する情報を確認することが望ましい。
しかし、本発明者は、ここで数学的アルゴリズムがバッテリーの内部バッテリーパラメータに対する高度で、且つ正確な推定を提供することができないということを認識している。なぜなら、それらは非線形動作特性(non‐linear operational characteristics)を持つバッテリーに対して充分に最適化されないからである。バッテリーは、一般的に非線形動作特性を持つので、より正確な方法が求められる。
結果的に、本発明者は、推定バッテリー状態と推定バッテリーパラメータとをより正確に定めるシステム及び方法に対する必要性を認識する。特に本発明者は、ここで推定バッテリー状態と推定バッテリーパラメータとを推定する結合された動作(combined behavior)と関連して、二つの連結されたフィルタ(two coupled filters)を活用する必要性を認識する。
本発明の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法が提供される。上記推定バッテリー状態ベクトルは、第1所定時間でバッテリー状態を示す。上記推定バッテリーパラメータベクトルは、上記第1所定時間で上記バッテリーパラメータを示す。上記方法は、上記第1所定時間以前の第2所定時間で上記バッテリー状態、バッテリー入力ノイズ、バッテリー出力変数を測定するセンサーに関わるセンサーノイズ、上記バッテリー状態の不確実性、上記バッテリー入力ノイズの不確実性及び上記センサーノイズの不確実性を示す複数個の推定増加バッテリー状態ベクトルを定める段階を含む。上記方法はまた、上記複数個の推定増加バッテリー状態ベクトルに基づいて上記第1所定時間で上記バッテリー状態と上記バッテリー状態の不確実性を示す複数個の予測バッテリー状態ベクトルを定める段階を含む。上記方法は、上記第1所定時間で上記バッテリーパラメータと上記バッテリーパラメータの不確実性を示す複数個の予測バッテリーパラメータベクトルを定める段階を含む。上記方法はまた、上記複数個の予測バッテリーパラメータベクトルに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第1複数個の予測バッテリー出力ベクトルを定める段階を含む。上記方法はまた、上記複数個の推定増加バッテリー状態ベクトルと上記複数個の予測バッテリー状態ベクトルとに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第2複数個の予測バッテリー出力ベクトルとを定める段階を含む。上記方法は、上記第1所定時間で獲得されたバッテリー出力変数の少なくとも一つの測定値を持つバッテリー出力ベクトルを定める段階を含む。上記方法は、上記複数個の予測バッテリーパラメータベクトル、第1複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルを定める段階を含む。上記方法はまた、上記複数個の予測バッテリー状態ベクトル、上記第2複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリー状態を示す推定バッテリー状態ベクトルを定める段階を含む。
本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステムが提供される。上記推定バッテリー状態ベクトルは、第1所定時間でバッテリー状態を示す。上記推定バッテリーパラメータベクトルは、第1所定時間で上記バッテリーパラメータを示す。上記システムは、上記バッテリーの出力変数を示す第1信号を生成するように構成されるセンサーを含む。上記システムはまた、上記センサーと動作可能に連結されたコンピュータを含む。上記コンピュータは、上記第1所定時間以前の第2所定時間でバッテリー状態、バッテリー入力ノイズ、バッテリー出力変数を測定するセンサーに関わるセンサーノイズ、上記バッテリー状態の不確実性、上記バッテリー入力ノイズの不確実性及び上記センサーノイズの不確実性を示す複数個の推定増加バッテリー状態ベクトルを定めるように構成される。上記コンピュータはまた、上記複数個の推定増加バッテリー状態ベクトルに基づいて上記第1所定時間で上記バッテリー状態と上記バッテリー状態の不確実性を示す複数個の予測バッテリー状態ベクトルを定めるように構成される。上記コンピュータはまた、上記第1所定時間で上記バッテリーパラメータと上記バッテリーパラメータの不確実性を示す複数個の予測バッテリーパラメータベクトルを定めるように構成される。上記コンピュータはまた、上記複数個の予測バッテリーパラメータベクトルに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第1複数個の予測バッテリー出力ベクトルを定めるように構成される。上記コンピュータはまた、上記複数個の推定増加バッテリー状態ベクトルと上記複数個の予測バッテリー状態ベクトルとに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第2複数個の予測バッテリー出力ベクトルを定めるように構成される。上記コンピュータはまた、上記第1信号に基づいてバッテリー出力ベクトルを定めるように構成される。上記コンピュータはまた、上記複数個の予測バッテリーパラメータベクトル、第1複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルを定めるように構成される。上記コンピュータはまた、上記複数個の予測バッテリー状態ベクトル、上記第2複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリー状態を示す推定バッテリー状態ベクトルを定めるように構成される。
上記実施例による他のシステム及び/または方法は以下の図面と詳細な説明に対する検討によって当業者に明らかであるか明らかになるであろう。そういう全ての追加的なシステムと方法は本発明の範囲にあり添付した請求項によって保護される。
図1を参照すれば、バッテリー12の状態を示す推定バッテリー状態ベクトルとバッテリー12のパラメータを示す推定バッテリーパラメータベクトルとを定めるシステム10が示される。上記バッテリー12は、少なくとも一つのバッテリーセル14を含む。もちろん、上記バッテリー12は複数個のバッテリーセルをさらに含むことができる。各バッテリーセルは再充電可能なセルであるか、或いは再充電できないセルであり得る。また、各バッテリーセルは本発明の技術分野の当業者に知られた電気‐化学構成を持つ陽極と陰極を用いて製作できる。
入力変数は、特定時間のバッテリー入力信号の値に定義される。例えば、上記入力変数は、バッテリーに入って行く電流中の一つと上記バッテリー温度とを含むことができる。出力変数は、特定時間のバッテリー出力信号の値に定義される。例えば、出力変数は、バッテリー出力電圧中の一つとバッテリー圧力(battery pressure)とを含むことができる。
上記システム10は、一つ以上の電圧センサー20、負荷回路26及びコンピュータのような演算ユニット28を含み、また一つ以上の温度センサー22及び電流センサー24を含むことができる。
上記電圧センサー20は、上記バッテリー12の一つ以上のバッテリーセルによって算出される電圧を示す第1出力信号を生成する。上記電圧センサー20は、上記コンピュータ28の入出力インターフェース(I/O interface)46と上記バッテリー12との間で電気的に連結される。上記電圧センサー20は、上記第1出力信号を上記コンピュータ28に伝送する。表現の明確性のために、ここでは一つの電圧センサーが示される。しかし、上記システム10の入れ替え可能な実施例においては複数個の電圧センサー(例えば、バッテリーセルごとに一つの電圧センサー)がシステム10に活用できることに注目すべきである。
上記温度センサー22は、上記バッテリー12の一つ以上の温度を示す第2出力信号を生成する。上記温度センサー22は、バッテリー12の近くに位置し上記コンピュータ28の入出力インターフェース46に電気的に連結される。上記温度センサー22は、上記コンピュータ28に上記第2出力信号を伝送する。表現の明確性のために、ここでは一つの温度センサーが示される。しかし、上記システム10の入れ替え可能な実施例においては複数個の温度センサー(例えば、バッテリーセルごとに一つの温度センサー)がシステム10に活用できることに注目すべきである。
上記電流センサー24は、上記バッテリー12のバッテリーセルによって供給されるか抜けた電流を示す第3出力信号を生成する。上記電流センサー24は、上記負荷回路26と上記バッテリー12との間に電気的に連結される。上記電流センサー24は、上記コンピュータ28の入出力インターフェース46とも電気的に連結される。上記電流センサー24は上記コンピュータ28に上記第3出力信号を伝送する。
上記負荷回路26は、電気的に上記電流センサー24と連結され上記バッテリー12から電流が供給されるか抜けるようになる。上記負荷回路26は、上記バッテリー12と電気的に連結できる全ての電気装置を含む。
上記コンピュータ28は、以下にさらに詳しく説明するように、上記バッテリー12の状態を示す推定バッテリー状態ベクトルと上記バッテリー(12)のパラメータを示す推定バッテリーパラメータベクトルを定めるために提供される。上記コンピュータ28は、中央処理装置(CPU)40、ROM44、RAM45のような揮発性メモリ及び入出力インターフェース46を含む。上記CPU40は、上記ROM44、RAM45及び入出力インターフェース46と作動可能に連結される。上記CPU40は、クロック42を含む。ROM44とRAM45を含むコンピュータで読み取り可能な媒体は、PROMs、EPROMs、EEPROMs、フラッシュメモリまたはデータの貯蔵が可能な何れの他の電気的、磁気的、光学的または組み合わせメモリ、CPU40によって使われた実行可能な命令を表現するもののような知られた多くのメモリ装置を使用することで具現される。
理解を助けるために、以下の方法の方程式で活用される表記法を説明する。上記サーカムフレックスシンボル(circumflex symbol)は、推定された大きさ(quantity)を示す(例えば、
Figure 0005036733
は、実際の大きさxに対する推定値を示す)。スーパースクリプトシンボル(superscript symbol)「−」は、アプリオリ(a priori)推定値を示す(すなわち、過去データに基づいた大きさの現在値を予測)。スーパースクリプトシンボル「+」は、アポステリオリ(a posteriori)推定値を示す(例えば、
Figure 0005036733
は、時間kを含んで全ての測定された値に基づいたタイムインデックスkにおける実際の大きさxに対する推定値である)。上記ティルデシンボル(tilde symbol)は、推定大きさ(estimated quantity)の誤差を示す
Figure 0005036733
上記バッテリー12に関わるバッテリーパラメータベクトルを定める方法論に対する詳しい説明を提供するに先立って、一般的な概括説明が提供される。バッテリー状態ベクトルは、上記バッテリー12に関わる少なくとも一つの充電状態(SOC)を含む。上記SOC値は動作するのに用いられ得る上記バッテリーの現在利用可能な容量を意味する0から100パーセントまで値である。
バッテリーセル動作(battery cell behavior)の数学的モデルが上記バッテリー12の状態ベクトルの推定値を演算するために上記方法で用いられる。上記バッテリーセルダイナミックス(battery cell dynamics)に対する数学的モデルは知られており、また以下に記述するように状態方程式(state equation)と出力方程式(output equation)とを含む離散時間状態空間モデル(discrete‐time state‐space model)を用いて表現できる。
上記バッテリー12に関わる上記状態ベクトルを定めるために用いられる上記状態方程式は次のようである。

=f(xk−1,uk−1,wk−1,k−1,k)

上記数式において、xは、タイムインデックス(time index)kにおいて上記バッテリー12に関わる状態ベクトルであり、uは、上記バッテリー12に対する知られた決定要素入力値(known/deterministic input)、wは、上記システムの状態に影響を及ぼす何れの測定されない入力値をモデリングするプロセスノイズ(process noise)または撹乱(disturbance)であり、そして、f(xk−1,uk−1,wk−1,k−1,k)は、状態変異関数(state transition function)である。
上記バッテリー12に関わる出力ベクトル(output vector)は、次の数式を利用して定められる。

=h(x,u,v,k)

上記数式において、h(x,u,v,k)は、測定関数(measurement function)であり、vは、非メモリモード(memory‐less mode)で上記バッテリー12の出力測定に影響を及ぼすが、上記バッテリー12の状態ベクトルには影響を及ぼさないセンサーノイズ(sensor noise)である。
システム状態xは、少なくとも現在出力を予測するために必要なセルの数学的モデルと現在入力とともに最小量の情報を含む。セル14の場合、状態は、例えばSOC、他の時定数に対する分極電圧(polarization voltage)レベル、ヒステリシスレベルを含む。上記システム外部からの入力uは、現在セルの電流iの最小値を含み、選択的にセル温度(もし状態に温度変化自体が設計されないなら)を含むことができる。システムパラメータθは、システム測定された入力と出力とが分かることで直接的に定められないことから、時間につれて単に徐々に変わる値である。このようなものとしては、セル容量、抵抗、分極電圧時定数、分極電圧調和要素、ヒステリシス調和要素、ヒステリシス率定数、効率要素などを含む。モデル出力yは、物理的に測定可能なセル量的要素または最小値で測定された量的要素から直接的に計算できるもの(例えば、負荷セル電圧)に対応する。
動的パラメータ(parameter dynamics)の数学的モデルも活用される。例示的なモデルは次の形態を持つ。
Figure 0005036733
第1方程式は、パラメータθが優先的に定数であることを示すが、rで表示される「ノイズ」プロセスによってモデルされたこの例から時間が経つにつれて徐々に変化できる。
出力dは、何れの推定誤差(estimation error)eが勘案されたg(・,・,・)によってモデルされた最適の動的パラメータ(optimum parameter dynamics)の関数である。上記最適の動的パラメータg(・,・,・)は、上記システム状態x、外部入力u及び時変数パラメータθセットに対する関数である。
図2ないし図8を参照しながら、本発明の望ましい実施例による上記バッテリー12の状態を示す推定バッテリー状態ベクトルと上記バッテリー(12)のパラメータを示す推定バッテリーパラメータベクトルとを定める方法を説明する。上記方法は、上記コントローラ28によって実行されるソフトウェアアルゴリズムを活用して具現できる。上記ソフトウェアアルゴリズムは、上記ROM44や上記RAM45または当業者に知られたコンピュータで読み取り可能な媒体に貯蔵される。
60段階において、上記コンピュータ28は、第1所定時間で獲得されたバッテリー入力変数のうち少なくとも一つの測定された値を持つバッテリー入力ベクトルuを生成する。
62段階において、上記コンピュータ28は、上記第1所定時間で獲得されたバッテリー出力変数のうち少なくとも一つの測定された値を持つバッテリー出力ベクトルyを生成する。
Figure 0005036733
66段階において、上記コンピュータ28は、以下の数式2を活用して上記第2所定時間でバッテリーパラメータ共分散の共分散を示す推定バッテリーパラメータベクトル共分散マトリックス(estimated battery parameter vector covariance matrix)
Figure 0005036733
とパラメータノイズに該当する共分散マトリックスに基づいて上記第1所定時間でバッテリーパラメータ共分散を示す予測バッテリーパラメータベクトル共分散マトリックス(predicted battery parameter vector covariance matrix)
Figure 0005036733
とを定める。
Figure 0005036733
68段階において、上記コンピュータ28は、以下の数式3を活用して第2所定時間で上記バッテリー12の状態を示す推定増加バッテリー状態ベクトル(estimated augmented battery state vector)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
は、上記第2所定時間でバッテリー入力ノイズベクトル(battery input noise vector)の予測値に該当し;
Figure 0005036733
は、上記第2所定時間でセンサーノイズベクトル(sensor noise vector)の予測値に該当する;そして、
Tは、マトリックス/ベクトル転置演算子(matrix/vectortranspose operator)である。
70段階において、上記コンピュータ28は、以下の数式4を活用して上記推定増加バッテリー状態ベクトルに関わる推定増加バッテリー状態ベクトル共分散マトリックス(estimated augmented battery state vector covariance matrix)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
diag()は、入力独立変数(input arguments)からブロック対角線行列(block‐diagonal matrix)を構成する関数である。
Figure 0005036733
74段階において、上記コンピュータ28は、以下の数式6を活用して上記第1所定時間でバッテリー状態をそれぞれ示す複数個の予測バッテリー状態ベクトル(predicted battery state vectors)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
80段階において、上記コンピュータ28は、以下の数式9を活用して上記第2所定時間でバッテリー12のパラメータをそれぞれ示す複数個の予測バッテリーパラメータベクトル
Figure 0005036733
を定める。
Figure 0005036733
82段階において、上記コンピュータ28は、以下の数式10を活用して上記第1所定時間で上記バッテリー12の出力をそれぞれ示す第1複数個の予測バッテリー出力ベクトル(a first plurality of predicted battery output vectors)
Figure 0005036733
を定める。
Figure 0005036733
上記数式において、
Figure 0005036733
Figure 0005036733
Figure 0005036733
86段階において、上記コンピュータ28は、以下の数式12を活用して上記第1所定時間でバッテリー出力をそれぞれ示す第2複数個の予測バッテリー出力ベクトル(a second plurality of predicted battery output vectors)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
88段階において、上記コンピュータ28は、以下の数式13を活用して上記第1所定時間でバッテリー出力値を示す第2予測バッテリー出力ベクトル(a second predicted batteryoutput vector)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
110段階において、上記コンピュータ28は、新しい第1及び第2所定時間を選択し、上記110段階以後、上記方法は60段階に戻る。
図9ないし図16を参照して、本発明の他の望ましい実施例によってスクエア‐ルートシグマ‐ポイントカルマンフィルタ(Square‐Root Sigma‐Point Kalman filters)を活用して上記バッテリー12に関わる推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法を説明する。上記方法は、上記コントローラ28によって実行されるソフトウェアアルゴリズムを活用して具現できる。上記ソフトウェアアルゴリズムは、上記ROM44やRAM45または技術分野の当業者に知られた他のコンピュータで読み取り可能な媒体に貯蔵される。
120段階において、上記コンピュータ28は、第1所定時間で獲得されたバッテリー入力変数のうち少なくとも一つの測定された値を持つバッテリー入力ベクトルu を生成する。
122段階において、上記コンピュータ28は、上記第1所定時間で獲得されたバッテリー出力変数のうち少なくとも一つの測定された値を持つバッテリー出力ベクトルy を生成する。
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
132段階において、上記コンピュータ28は、以下の数式28を活用して上記推定増加バッテリー状態ベクトルに関わる推定増加バッテリー状態ベクトルスクエア‐ルート共分散マトリックス(estimated augmented battery state vector square‐root covariance matrix)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
Figure 0005036733
136段階において、上記コンピュータ28は、以下の数式30を活用して上記第1所定時間でバッテリー12の状態をそれぞれ示す複数個の予測バッテリー状態ベクトル(predicted battery state vectors)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
140段階において、上記コンピュータ28は、以下の数式32を活用して予測バッテリー状態ベクトルスクエア‐ルート共分散マトリックス
Figure 0005036733
を定める。
Figure 0005036733
142段階において、上記コンピュータ28は、以下の数式33を活用して上記第1所定時間でバッテリー12のパラメータをそれぞれ示す複数個の予測バッテリーパラメータベクトル
Figure 0005036733
を定める。
Figure 0005036733
144段階において、上記コンピュータ28は、以下の数式34を活用して上記第1所定時間で上記バッテリー12の出力をそれぞれ示す第1複数個の予測バッテリー出力ベクトル(a first plurality of predicted battery output vectors)
Figure 0005036733
を定める。
Figure 0005036733
上記数式において、
Figure 0005036733
Figure 0005036733
Figure 0005036733
148段階において、上記コンピュータ28は、数式36を活用して上記第1所定時間でバッテリー出力をそれぞれ示す第2複数個の予測バッテリー出力ベクトル(a second plurality of predicted batteryoutput vectors)
Figure 0005036733
を定める。
Figure 0005036733
上記数式において、
Figure 0005036733
150段階において、上記コンピュータ28は、数式37を活用して上記第1所定時間でバッテリー出力値を示す第2予測バッテリー出力ベクトル(a second predicted batteryoutput vector)
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
上記数式において、
Figure 0005036733
160段階において、上記コンピュータ28は、以下の数式42を活用して第2予測クロス‐共分散マトリックス
Figure 0005036733
を定める。
Figure 0005036733
Figure 0005036733
Figure 0005036733
Figure 0005036733
上記数式において、downdate{ }は、第2独立変数(second argument)を使用して第1独立変数に対するマトリックスダウンデートオペレーション(matrix downdate operation)を演算する。
168段階において、上記コンピュータ28は、以下の数式46を活用して上記第1所定時間で上記バッテリー12のパラメータを示す推定バッテリーパラメータベクトル
Figure 0005036733
を定める。
Figure 0005036733
170段階において、上記コンピュータ28は、以下の数式47を活用して上記推定バッテリーパラメータベクトルに関わる推定バッテリーパラメータベクトルスクエア‐ルート共分散マトリックス(estimated battery parameter vector square‐root covariance matrix)
Figure 0005036733
を定める。
Figure 0005036733
172段階において、上記コンピュータ28は、新しい第1及び第2所定時間を選択する。172段階以後、上記方法は120段階に戻る。
推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルの両方を定める上記システム及び方法は、他のシステムと方法とに比べて実質的な利点を提供する。特に、上記システム及び方法は、非線形動作特性を持つバッテリーに対して上記推定バッテリー状態ベクトルと上記推定バッテリーパラメータベクトルとをさらに正確に定めることができる技術的な効果を提供する。
上記方法は、プロッピィーディスク、CD‐ROM、ハードディスクまたはコンピュータで認識可能な貯蔵媒体のような媒体で具体化される命令語を含むコンピュータプログラムコードの形態でも実現することができる。ここで上記コンピュータプログラムコードは、コンピュータ内にローディングされコンピュータによって実行され、上記コンピュータは上記方法の実行が可能な装備になる。本方法は、コンピュータプログラムコードの形態、例えば、貯蔵媒体に貯蔵された形態、コンピュータによってローディングおよび/または実行される形態、変調されたキャリアウエーブの形態、伝送媒体、電気配線またはケーブル、光ファイバ、電磁気放射の経由などの形態でも実現することができる。上記コンピュータプログラムコードはコンピュータ内にローディングされコンピュータにより実行される場合、 上記コンピュータは上記方法の実行が可能な装備になる。汎用マイクロプロセッサーで具現される場合、上記コンピュータプログラムコードセグメントは上記マイクロプロセッサーが特定の論理回路を作るように構成される。
望ましい実施例を参照しながら本発明が説明されているが、当業者により多様な変形例が可能であり、均等な要素が本発明の範囲を逸脱しない領域で本発明の要素に代替できることが理解されるべきである。また、多くの変形例が本発明の本質的領域を逸脱することなく本発明の教示に対する特別な状況と内容に適用されることができる。従って、本発明は、本発明を実施するために開示される実施例に限られず、本発明は添付された請求の範囲内にある全ての実施例を含む。さらに、第1、第2などに対する用語の使用は重要さの手順を意味するよりは、一つの要素を他の要素と区別するために使われる。
本発明の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステムの概路図である。 本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明の他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。 本発明のさらに他の望ましい実施例による推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法の流れ図である。

Claims (28)

  1. 第1所定時間でバッテリー状態を示す推定バッテリー状態ベクトルと上記第1所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルとを定める方法であって、
    上記第1所定時間以前の第2所定時間で上記バッテリー状態、バッテリー入力ノイズ、バッテリー出力変数を測定するセンサーに関わるセンサーノイズ、上記バッテリー状態の不確実性、上記バッテリー入力ノイズの不確実性及び上記センサーノイズの不確実性を示す複数個の推定増加バッテリー状態ベクトルを定める段階;
    上記複数個の推定増加バッテリー状態ベクトルに基づいて上記第1所定時間で上記バッテリー状態と上記バッテリー状態の不確実性を示す複数個の予測バッテリー状態ベクトルを定める段階;
    上記第1所定時間で上記バッテリーパラメータと上記バッテリーパラメータの不確実性を示す複数個の予測バッテリーパラメータベクトルを定める段階;
    上記複数個の予測バッテリーパラメータベクトルに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第1複数個の予測バッテリー出力ベクトルを定める段階;
    上記複数個の推定増加バッテリー状態ベクトルと上記複数個の予測バッテリー状態ベクトルとに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第2複数個の予測バッテリー出力ベクトルを定める段階;
    上記第1所定時間で獲得されたバッテリー出力変数の少なくとも一つの測定値を持つバッテリー出力ベクトルを定める段階;
    上記複数個の予測バッテリーパラメータベクトル、第1複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルを定める段階;及び
    上記複数個の予測バッテリー状態ベクトル、上記第2複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリー状態を示す推定バッテリー状態ベクトルを定める段階を含むことを特徴とする推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  2. 上記複数個の推定増加バッテリー状態ベクトルを定める段階は、
    メモリーから上記第2所定時間で上記バッテリー状態を示す推定バッテリー状態ベクトルを検索(retrieving)する段階;
    上記推定バッテリー状態ベクトル、予測バッテリー入力ノイズベクトル及び予測センサーノイズベクトルに基づいて上記第2所定時間で上記バッテリー状態、上記バッテリー入力ノイズ及び上記センサーノイズを示す第1推定増加バッテリー状態ベクトルを定める段階;
    上記メモリーから上記第2所定時間で上記バッテリー状態の不確実性を示す推定バッテリー状態ベクトル共分散マトリックスを検索する段階;
    上記推定バッテリー状態ベクトル共分散マトリックス、上記バッテリー入力ノイズの不確実性を示す共分散マトリックス、上記センサーノイズの不確実性を示す共分散マトリックスに基づいて上記バッテリー状態の不確実性、上記バッテリー入力ノイズ及び上記センサーノイズを示す第1推定増加バッテリー状態ベクトル共分散マトリックスを定める段階;及び
    上記第1推定増加バッテリー状態ベクトルと上記第1推定増加バッテリー状態ベクトル共分散マトリックスに基づいて上記複数個の推定増加バッテリー状態ベクトルを演算する段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  3. 上記複数個の推定増加バッテリー状態ベクトルを演算する段階は、
    上記複数個の推定増加バッテリー状態ベクトルのうち一つのメンバーを上記第1推定増加バッテリー状態ベクトルと同一にセットする段階;
    上記複数個の推定増加バッテリー状態ベクトルのLメンバーの第1追加セット(a first additional set)を上記第1推定増加バッテリー状態ベクトル共分散マトリックスのマトリックススクエア‐ルートから抽出された各列によってかけられた定数に足された上記第1推定増加バッテリー状態ベクトルと同一にセットする段階;及び
    上記第1複数個の推定増加バッテリー状態ベクトルのLメンバーの第2追加セットを上記第1推定増加バッテリー状態ベクトルから引いた(subtracted)上記第1推定増加バッテリー状態ベクトル共分散マトリックスの上記マトリックススクエア‐ルートから抽出された各列によってかけられた定数と同一にセットする段階をさらに含み、
    上記Lは、上記第1推定増加バッテリー状態ベクトルの長さであることを特徴とする請求項2に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  4. 上記複数個の推定増加バッテリー状態ベクトルを定める段階は、
    メモリーから上記第2所定時間で上記バッテリー状態を示す推定バッテリー状態ベクトルを検索する段階;
    上記推定バッテリー状態ベクトル、予測バッテリー入力ノイズベクトル及び予測センサーノイズベクトルに基づいて上記第2所定時間で上記バッテリー状態、バッテリー入力ノイズ及び上記センサーノイズを示す第1推定増加バッテリー状態ベクトルを定める段階;
    上記メモリーから上記第2所定時間で上記バッテリー状態の不確実性を示す推定バッテリー状態ベクトルスクエア‐ルート共分散マトリックスを検索する段階;
    上記推定バッテリー状態ベクトルスクエア‐ルート共分散マトリックス、上記バッテリー入力ノイズの不確実性を示す共分散マトリックス及び上記センサーノイズの不確実性を示す共分散マトリックスに基づいて上記バッテリー状態の不確実性、上記バッテリー入力ノイズ及び上記センサーノイズを示す第1推定増加バッテリー状態ベクトルスクエア‐ルート共分散マトリックスを定める段階;及び
    上記第1推定増加バッテリー状態ベクトルと上記第1推定増加バッテリー状態ベクトルスクエア‐ルート共分散マトリックスに基づいて上記複数個の推定増加バッテリー状態ベクトルを演算する段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  5. 上記複数個の推定増加バッテリー状態ベクトルを定める段階は、
    上記複数個の推定増加バッテリー状態ベクトルのうち一つのメンバーを上記第1推定増加バッテリー状態ベクトルと同一にセットする段階;
    上記複数個の推定増加バッテリー状態ベクトルのLメンバーの第1追加セットを上記推定増加バッテリー状態ベクトルスクエア‐ルート共分散マトリックスから抽出された各列によってかけられた定数に足された上記第1推定増加バッテリー状態ベクトルと同一にセットする段階;
    上記複数個の推定増加バッテリー状態ベクトルのLメンバーの第2追加セットを上記第1推定増加バッテリー状態ベクトルから引いた上記第1推定増加バッテリー状態ベクトルスクエア‐ルート共分散マトリックスから抽出された各列によってかけられた定数と同一にセットする段階をさらに含み、
    上記Lは、上記推定増加バッテリー状態ベクトルの長さであることを特徴とする請求項4に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  6. 上記複数個の予測バッテリー状態ベクトルを定める段階は、
    第1複数個の推定バッテリー状態ベクトルを獲得するために上記複数個の推定増加バッテリー状態ベクトルから上記第2所定時間で上記バッテリー状態と上記バッテリー状態の不確実性を示す値を抽出する段階;
    複数個の推定バッテリー入力ノイズベクトルを獲得するために上記第1複数個の推定増加バッテリー状態ベクトルから上記第2所定時間で上記バッテリー入力ノイズと上記バッテリー入力ノイズの不確実性を示す値を抽出する段階;
    上記第2所定時間でバッテリー入力変数のうち少なくとも一つの測定値を持つ上記第1バッテリー入力ベクトルを生成する段階;及び
    上記第1複数個の予測バッテリー状態ベクトル、上記第1複数個の予測バッテリー入力ノイズベクトル及び上記第1バッテリー入力ベクトルに基づいて上記複数個の予測バッテリー状態ベクトルを定める段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  7. 上記複数個の予測バッテリーパラメータベクトルを定める段階は、
    メモリーから上記第2所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルを検索する段階;
    上記メモリーから上記第2所定時間で上記バッテリーパラメータの不確実性を示す推定バッテリーパラメータベクトル共分散マトリックスを検索する段階;
    上記第2所定時間でパラメータノイズの不確実性を示す推定バッテリーパラメータノイズ共分散マトリックスを定める段階;及び
    上記推定バッテリーパラメータベクトル、上記推定バッテリーパラメータベクトル共分散マトリックス及び上記推定バッテリーパラメータノイズ共分散マトリックスに基づいて上記複数個の予測バッテリーパラメータベクトルを演算する段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  8. 上記複数個の予測バッテリーパラメータベクトルを演算する段階は、
    上記複数個の予測バッテリーパラメータベクトルのうち一つのメンバーを上記推定バッテリーパラメータベクトルと同一にセットする段階;
    上記推定バッテリーパラメータベクトル共分散マトリックスと上記推定バッテリーパラメータノイズ共分散マトリックスに基づいて上記第1所定時間で上記バッテリーパラメータの不確実性を示す予測バッテリーパラメータベクトル共分散マトリックスを定める段階;
    上記複数個の予測バッテリーパラメータベクトルのLメンバーの第1追加セット(a first additional set)を上記第1予測バッテリーパラメータベクトル共分散マトリックスのマトリックススクエア‐ルートから抽出された各列によってかけられた定数に足された上記第1推定バッテリーパラメータベクトルと同一にセットする段階;及び
    上記複数個の予測バッテリーパラメータベクトルのLメンバーの第2追加セットを上記第1推定増加バッテリーパラメータベクトルから引いた上記第1予測バッテリーパラメータベクトル共分散マトリックスのマトリックススクエア‐ルートから抽出された各列によってかけられた定数と同一にセットする段階をさらに含み、
    上記Lは、上記第1推定バッテリーパラメータベクトルの長さであることを特徴とする請求項7に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  9. 上記複数個の予測バッテリーパラメータベクトルを定める段階は、
    メモリーから上記第1所定時間以前の第2所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルを検索する段階;
    上記メモリーから上記第2所定時間で上記バッテリーパラメータの不確実性を示す推定バッテリーパラメータベクトルスクエア‐ルート共分散マトリックスを検索する段階;
    上記第2所定時間でパラメータノイズの不確実性を示す推定バッテリーパラメータノイズスクエア‐ルート共分散マトリックスを定める段階;及び
    上記推定バッテリーパラメータベクトル、上記推定バッテリーパラメータベクトルスクエア‐ルート共分散マトリックス及び上記推定バッテリーパラメータノイズスクエア‐ルート共分散マトリックスに基づいて上記複数個の予測バッテリーパラメータベクトルを演算する段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  10. 上記複数個の予測バッテリーパラメータベクトルを演算する段階は、
    上記複数個の予測バッテリーパラメータベクトルのうち一つのメンバーを上記推定バッテリーパラメータベクトルと同一にセットする段階;
    上記推定バッテリーパラメータベクトルスクエア‐ルート共分散マトリックスと上記推定バッテリーパラメータノイズ共分散マトリックスに基づいて上記第1所定時間で上記バッテリーパラメータの不確実性を示す予測バッテリーパラメータベクトルスクエア‐ルート共分散マトリックスを定める段階;
    上記複数個の予測バッテリーパラメータベクトルのLメンバーの第1追加セットを上記第1予測バッテリーパラメータベクトルスクエア‐ルート共分散マトリックスから抽出された各列によってかけられた定数に足された上記第1推定バッテリーパラメータベクトルと同一にセットする段階;及び
    上記複数個の予測バッテリーパラメータベクトルのLメンバーの第2追加セットを上記第1推定増加バッテリーパラメータベクトルから引いた上記第1予測バッテリーパラメータベクトルスクエア‐ルート共分散マトリックスから抽出された各列によってかけられた定数と同一にセットする段階をさらに含み、
    上記Lは、上記第1推定バッテリーパラメータベクトルの長さであることを特徴とする請求項9に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  11. 上記第1複数個の予測バッテリー出力ベクトルを定める段階は、
    メモリーから上記第2所定時間で上記バッテリー状態を示す推定バッテリー状態ベクトルを検索する段階;
    上記第1所定時間で平均センサーノイズを示すバッテリーミーンセンサーノイズベクトル(battery mean sensor noise vector)と上記第2所定時間で平均バッテリー入力ノイズを示すバッテリーミーン入力ノイズベクトル(battery mean input noise vector)を生成する段階;
    上記第1所定時間でバッテリー入力変数のうち少なくとも一つの測定値を持つ第1バッテリー入力ベクトルを生成する段階;
    上記第2所定時間でバッテリー入力変数のうち少なくとも一つの測定値を持つ第2バッテリー入力ベクトルを生成する段階;及び
    上記複数個の予測バッテリーパラメータベクトル、上記第1バッテリー入力ベクトル、上記第2バッテリー入力ベクトル、上記バッテリーミーンセンサーノイズベクトル、上記バッテリーミーン入力ノイズベクトル及び上記推定バッテリー状態ベクトルに基づいて上記第1複数個の予測バッテリー出力ベクトルを定める段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  12. 上記第2複数個の予測バッテリー出力ベクトルを定める段階は、
    複数個の推定センサーノイズベクトルを獲得するために上記複数個の推定増加バッテリー状態ベクトルから上記第2所定時間で上記センサーノイズ及び上記センサーノイズの不確実性を示す値を抽出する段階;
    上記第1所定時間でバッテリー入力変数のうち少なくとも一つの測定値を持つ上記第1バッテリー入力ベクトルを生成する段階;及び
    上記複数個の予測バッテリー状態ベクトル、上記第1バッテリー入力ベクトル、上記複数個の推定センサーノイズベクトルに基づいて上記第2複数個の予測バッテリー出力ベクトルを定める段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  13. 上記推定バッテリーパラメータベクトルを定める段階は、
    上記第1複数個の予測バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリーの少なくとも一つの出力変数を示す第1予測バッテリー出力ベクトルを定める段階;
    上記第2所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルに基づいて上記第1所定時間で上記バッテリーパラメータを示す予測バッテリーパラメータベクトルを定める段階;
    上記予測バッテリーパラメータベクトル、上記第1予測バッテリー出力ベクトル、上記複数個の予測バッテリーパラメータベクトル及び上記複数個の予測バッテリー出力ベクトルに基づいてパラメータゲインマトリックスを定める段階;及び
    上記バッテリー出力ベクトル、上記第1予測バッテリー出力ベクトル、上記予測バッテリーパラメータベクトル及び上記パラメータゲインマトリックスに基づいて上記推定バッテリーパラメータベクトルを演算する段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  14. 上記第1予測バッテリー出力ベクトルは、上記複数個の予測バッテリー出力ベクトルの加重値平均を演算する段階を含むことを特徴とする請求項13に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  15. 上記パラメータゲインマトリックスを定める段階は、
    上記第1予測バッテリー出力ベクトルと上記複数個の予測バッテリー出力ベクトルに基づいて予測バッテリー出力ベクトル共分散マトリックスを定める段階;
    上記第1予測バッテリー出力ベクトル、上記複数個の予測バッテリー出力ベクトル、上記予測バッテリーパラメータベクトル及び上記複数個の予測バッテリーパラメータベクトルに基づいてバッテリー出力ベクトルとバッテリーパラメータベクトルの間の予測クロス‐共分散マトリックスを定める段階;及び
    上記予測クロス‐共分散マトリックスと上記予測バッテリー出力ベクトル共分散マトリックスとに基づいて上記パラメータゲインマトリックスを演算する段階を含むことを特徴とする請求項13に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  16. 上記予測バッテリー出力ベクトル共分散マトリックスを定める段階は、
    上記第1予測バッテリー出力ベクトルと上記複数個の予測バッテリー出力ベクトルに基づいて予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスを定める段階;及び
    上記予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスと上記予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスに基づいて上記予測バッテリー出力ベクトル共分散マトリックスを演算する段階を含むことを特徴とする請求項15に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  17. 上記推定バッテリーパラメータベクトルを演算する段階は、
    上記バッテリー出力ベクトルと上記第1予測バッテリー出力ベクトルに基づいてイノベーションベクトル(innovation vector)を定める段階;
    上記パラメータゲインマトリックスと上記イノベーションベクトルに基づいてアップデートベクトル(update vector)を定める段階;及び
    上記第1予測バッテリーパラメータベクトルと上記アップデートベクトルに基づいて上記推定バッテリーパラメータベクトルを演算する段階を含むことを特徴とする請求項13に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  18. 上記第1予測バッテリーパラメータベクトル、上記第1複数個の予測バッテリーパラメータベクトル、上記パラメータゲインマトリックス及び上記第1予測バッテリー出力ベクトル共分散マトリックスに基づいて上記第1所定時間で推定バッテリーパラメータベクトル共分散マトリックスを定める段階をさらに含むことを特徴とする請求項13に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  19. 上記第1予測バッテリーパラメータベクトル、上記第1複数個の予測バッテリーパラメータベクトル、上記パラメータゲインマトリックス及び予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスに基づいて上記第1所定時間で推定バッテリーパラメータベクトルスクエア‐ルート共分散マトリックスを定める段階をさらに含むことを特徴とする請求項13に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  20. 上記推定バッテリー状態ベクトルを定める段階は、
    (1)上記複数個の予測バッテリー状態ベクトルに基づいて上記第1所定時間で上記バッテリー状態を示す予測バッテリー状態ベクトルを定める段階;
    (2)上記複数個の予測出力ベクトルに基づいて上記第1所定時間で上記バッテリーの少なくとも一つの出力変数を示す第2予測バッテリー出力ベクトルを定める段階;
    (3)上記予測バッテリー状態ベクトル、上記第2予測バッテリー出力ベクトル、上記複数個の予測バッテリー状態ベクトル及び上記第2複数個の予測バッテリー出力ベクトルに基づいて状態ゲインマトリックス(state gain matrix)を定める段階;及び
    (4)上記予測バッテリー状態ベクトル、上記第2予測バッテリー出力ベクトル、上記状態ゲインマトリックス及び上記バッテリー出力ベクトルに基づいて上記推定バッテリー状態ベクトルを演算する段階を含むことを特徴とする請求項1に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  21. 上記予測バッテリー状態ベクトルを定める段階は、
    上記複数個の予測バッテリー状態ベクトルの加重値平均を演算する段階を含むことを特徴とする請求項20に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  22. 上記第2予測バッテリー出力ベクトルを定める段階は、
    上記第2複数個の予測バッテリー出力ベクトルの加重値平均を演算する段階を含むことを特徴とする請求項20に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  23. 上記状態ゲインマトリックスを定める段階は、
    上記第2予測バッテリー出力ベクトルと上記第2複数個の予測バッテリー出力ベクトルに基づいて予測バッテリー出力ベクトル共分散マトリックスを定める段階;
    上記第2予測バッテリー出力ベクトル、上記第2複数個の予測バッテリー出力ベクトル、上記予測バッテリー状態ベクトル及び上記複数個の予測バッテリー状態ベクトルに基づいてバッテリー出力ベクトルとバッテリー状態ベクトルの間の予測クロス‐共分散マトリックスを定める段階;及び
    上記予測クロス‐共分散マトリックスと上記予測バッテリー出力ベクトル共分散マトリックスとに基づいて上記状態ゲインマトリックスを演算する段階を含むことを特徴とする請求項20に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  24. 上記予測バッテリー出力ベクトル共分散マトリックスを定める段階は、
    上記第2予測バッテリー出力ベクトルと上記第2複数個の予測バッテリー出力ベクトルに基づいて予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスを定める段階;及び
    上記予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスと上記予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスに基づいて上記予測バッテリー出力ベクトル共分散マトリックスを演算する段階を含むことを特徴とする請求項23に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  25. 上記推定バッテリー状態ベクトルを演算する段階は、
    上記第2バッテリー出力ベクトルと上記第2予測バッテリー出力ベクトルに基づいてイノベーションベクトルを定める段階;
    上記状態ゲインマトリックスと上記イノベーションベクトルに基づいてアップデートベクトルを定める段階;及び
    上記予測バッテリー状態ベクトルと上記アップデートベクトルに基づいて上記推定バッテリー状態ベクトルを演算する段階を含むことを特徴とする請求項20に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  26. 上記予測バッテリー状態ベクトル、上記複数個の予測バッテリー状態ベクトル、上記状態ゲインマトリックス及び上記予測バッテリー出力ベクトル共分散マトリックスに基づいて上記第1所定時間で推定バッテリー状態ベクトル共分散マトリックスを定める段階をさらに含む特徴とする請求項20に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  27. 上記予測バッテリー状態ベクトル、上記複数個の予測バッテリー状態ベクトル、上記状態ゲインマトリックス及び予測バッテリー出力ベクトルスクエア‐ルート共分散マトリックスに基づいて上記第1所定時間で推定バッテリー状態ベクトルスクエア‐ルート共分散マトリックスを定める段階をさらに含む特徴とする請求項20に記載の推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定める方法。
  28. 第1所定時間でバッテリー状態を示す推定バッテリー状態ベクトルと上記第1所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルを定めるシステムであって、
    上記バッテリーの出力変数を示す第1信号を生成するように構成されるセンサー;
    上記センサーと動作可能に連結され、上記第1所定時間以前の第2所定時間でバッテリー状態、バッテリー入力ノイズ、バッテリー出力変数を測定するセンサーに関わるセンサーノイズ、上記バッテリー状態の不確実性、上記バッテリー入力ノイズの不確実性及び上記センサーノイズの不確実性を示す複数個の推定増加バッテリー状態ベクトルを定めるように構成され、上記複数個の推定増加バッテリー状態ベクトルに基づいて上記第1所定時間で上記バッテリー状態と上記バッテリー状態の不確実性を示す複数個の予測バッテリー状態ベクトルを定めるように構成され、上記第1所定時間で上記バッテリーパラメータと上記バッテリーパラメータの不確実性を示す複数個の予測バッテリーパラメータベクトルを定めるように構成され、上記複数個の予測バッテリーパラメータベクトルに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第1複数個の予測バッテリー出力ベクトルを定めるように構成され、上記複数個の推定増加バッテリー状態ベクトルと上記複数個の予測バッテリー状態ベクトルに基づいて上記第1所定時間で少なくとも一つの上記バッテリーの出力変数と上記出力変数の不確実性を示す第2複数個の予測バッテリー出力ベクトルを定めるように構成され、上記第1信号に基づいてバッテリー出力ベクトルを定めるように構成され、上記複数個の予測バッテリーパラメータベクトル、第1複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリーパラメータを示す推定バッテリーパラメータベクトルを定めるように構成され、上記複数個の予測バッテリー状態ベクトル、上記第2複数個の予測バッテリー出力ベクトル及び上記バッテリー出力ベクトルに基づいて上記第1所定時間で上記バッテリー状態を示す推定バッテリー状態ベクトルを定めるように構成されるコンピュータを含むことを特徴とする推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステム。
JP2008557199A 2006-03-02 2007-01-22 推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステム及び方法 Active JP5036733B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US11/366,089 US7521895B2 (en) 2006-03-02 2006-03-02 System and method for determining both an estimated battery state vector and an estimated battery parameter vector
US11/366,089 2006-03-02
PCT/KR2007/000360 WO2007100189A1 (en) 2006-03-02 2007-01-22 System and method for determining both an estimated battery state vector and an estimated battery parameter vector

Publications (2)

Publication Number Publication Date
JP2009528662A JP2009528662A (ja) 2009-08-06
JP5036733B2 true JP5036733B2 (ja) 2012-09-26

Family

ID=38459266

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008557199A Active JP5036733B2 (ja) 2006-03-02 2007-01-22 推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステム及び方法

Country Status (7)

Country Link
US (1) US7521895B2 (ja)
EP (1) EP1989563B1 (ja)
JP (1) JP5036733B2 (ja)
KR (1) KR101355959B1 (ja)
CN (1) CN101395489B (ja)
TW (1) TWI337663B (ja)
WO (1) WO2007100189A1 (ja)

Families Citing this family (57)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7321220B2 (en) * 2003-11-20 2008-01-22 Lg Chem, Ltd. Method for calculating power capability of battery packs using advanced cell model predictive techniques
US7723957B2 (en) * 2005-11-30 2010-05-25 Lg Chem, Ltd. System, method, and article of manufacture for determining an estimated battery parameter vector
US8049460B2 (en) * 2007-07-18 2011-11-01 Tesla Motors, Inc. Voltage dividing vehicle heater system and method
US8628872B2 (en) * 2008-01-18 2014-01-14 Lg Chem, Ltd. Battery cell assembly and method for assembling the battery cell assembly
US7994755B2 (en) 2008-01-30 2011-08-09 Lg Chem, Ltd. System, method, and article of manufacture for determining an estimated battery cell module state
US7883793B2 (en) * 2008-06-30 2011-02-08 Lg Chem, Ltd. Battery module having battery cell assemblies with alignment-coupling features
US9759495B2 (en) 2008-06-30 2017-09-12 Lg Chem, Ltd. Battery cell assembly having heat exchanger with serpentine flow path
US8426050B2 (en) * 2008-06-30 2013-04-23 Lg Chem, Ltd. Battery module having cooling manifold and method for cooling battery module
US8067111B2 (en) * 2008-06-30 2011-11-29 Lg Chem, Ltd. Battery module having battery cell assembly with heat exchanger
US9140501B2 (en) * 2008-06-30 2015-09-22 Lg Chem, Ltd. Battery module having a rubber cooling manifold
KR100927541B1 (ko) * 2008-08-14 2009-11-17 주식회사 엘지화학 배터리 전압 거동을 이용한 배터리 저항 특성 추정 장치 및방법
US8202645B2 (en) 2008-10-06 2012-06-19 Lg Chem, Ltd. Battery cell assembly and method for assembling the battery cell assembly
US9337456B2 (en) 2009-04-20 2016-05-10 Lg Chem, Ltd. Frame member, frame assembly and battery cell assembly made therefrom and methods of making the same
US8663828B2 (en) 2009-04-30 2014-03-04 Lg Chem, Ltd. Battery systems, battery module, and method for cooling the battery module
US8663829B2 (en) 2009-04-30 2014-03-04 Lg Chem, Ltd. Battery systems, battery modules, and method for cooling a battery module
US8403030B2 (en) 2009-04-30 2013-03-26 Lg Chem, Ltd. Cooling manifold
US8852778B2 (en) 2009-04-30 2014-10-07 Lg Chem, Ltd. Battery systems, battery modules, and method for cooling a battery module
US8399118B2 (en) 2009-07-29 2013-03-19 Lg Chem, Ltd. Battery module and method for cooling the battery module
US8703318B2 (en) * 2009-07-29 2014-04-22 Lg Chem, Ltd. Battery module and method for cooling the battery module
US8399119B2 (en) 2009-08-28 2013-03-19 Lg Chem, Ltd. Battery module and method for cooling the battery module
US8041522B2 (en) * 2009-12-02 2011-10-18 American Electric Vehicles, Ind. System and method for recursively estimating battery cell total capacity
US8427105B2 (en) * 2009-12-02 2013-04-23 Gregory L. Plett System and method for equalizing a battery pack during a battery pack charging process
US8918299B2 (en) * 2009-12-02 2014-12-23 American Electric Vehicles, Inc. System and method for maximizing a battery pack total energy metric
RU2491566C1 (ru) * 2010-02-18 2013-08-27 Ниссан Мотор Ко., Лтд. Устройство оценки состояния батареи и способ оценки состояния батареи
JP5691592B2 (ja) 2010-02-18 2015-04-01 日産自動車株式会社 電池状態推定装置
US8341449B2 (en) 2010-04-16 2012-12-25 Lg Chem, Ltd. Battery management system and method for transferring data within the battery management system
US9147916B2 (en) 2010-04-17 2015-09-29 Lg Chem, Ltd. Battery cell assemblies
US8353315B2 (en) 2010-08-23 2013-01-15 Lg Chem, Ltd. End cap
US8920956B2 (en) 2010-08-23 2014-12-30 Lg Chem, Ltd. Battery system and manifold assembly having a manifold member and a connecting fitting
US8758922B2 (en) 2010-08-23 2014-06-24 Lg Chem, Ltd. Battery system and manifold assembly with two manifold members removably coupled together
US8469404B2 (en) 2010-08-23 2013-06-25 Lg Chem, Ltd. Connecting assembly
US9005799B2 (en) 2010-08-25 2015-04-14 Lg Chem, Ltd. Battery module and methods for bonding cell terminals of battery cells together
US8662153B2 (en) 2010-10-04 2014-03-04 Lg Chem, Ltd. Battery cell assembly, heat exchanger, and method for manufacturing the heat exchanger
US8288031B1 (en) 2011-03-28 2012-10-16 Lg Chem, Ltd. Battery disconnect unit and method of assembling the battery disconnect unit
CN102253343B (zh) * 2011-04-21 2013-04-03 北京世纪瑞尔技术股份有限公司 一种蓄电池的健康度与荷电状态的估算方法
US8449998B2 (en) 2011-04-25 2013-05-28 Lg Chem, Ltd. Battery system and method for increasing an operational life of a battery cell
US9178192B2 (en) 2011-05-13 2015-11-03 Lg Chem, Ltd. Battery module and method for manufacturing the battery module
FR2975501B1 (fr) * 2011-05-20 2013-05-31 Renault Sas Procede d'estimation de l'etat de charge d'une batterie electrique
US10234512B2 (en) * 2011-06-11 2019-03-19 Sendyne Corporation Current-based cell modeling
US8880253B2 (en) 2011-06-28 2014-11-04 Ford Global Technologies, Llc Nonlinear adaptive observation approach to battery state of charge estimation
US8974929B2 (en) 2011-06-30 2015-03-10 Lg Chem, Ltd. Heating system for a battery module and method of heating the battery module
US8859119B2 (en) 2011-06-30 2014-10-14 Lg Chem, Ltd. Heating system for a battery module and method of heating the battery module
US8974928B2 (en) 2011-06-30 2015-03-10 Lg Chem, Ltd. Heating system for a battery module and method of heating the battery module
US8993136B2 (en) 2011-06-30 2015-03-31 Lg Chem, Ltd. Heating system for a battery module and method of heating the battery module
US9496544B2 (en) 2011-07-28 2016-11-15 Lg Chem. Ltd. Battery modules having interconnect members with vibration dampening portions
US9182451B2 (en) * 2012-07-30 2015-11-10 Robert Bosch Gmbh System and method for posteriori adaptation of a state of charge model in a battery
US8751086B2 (en) 2012-08-21 2014-06-10 Ford Global Technologies, Llc Online battery capacity estimation
KR101547005B1 (ko) * 2012-10-26 2015-08-24 주식회사 엘지화학 배터리 잔존 용량 추정 장치 및 방법
KR101547006B1 (ko) * 2012-10-26 2015-08-24 주식회사 엘지화학 배터리 잔존 용량 추정 장치 및 방법
KR101767635B1 (ko) 2014-10-24 2017-08-14 주식회사 엘지화학 이차 전지의 충전 상태 추정 장치 및 그 방법
KR101655583B1 (ko) * 2014-11-26 2016-09-07 현대자동차주식회사 배터리 상태 추정 장치 및 방법
FR3030769B1 (fr) 2014-12-22 2018-02-02 Renault S.A.S. Procede d'estimation de grandeurs physiques caracteristiques d'une batterie electrique
US10324135B2 (en) * 2016-06-06 2019-06-18 Mitsubishi Electric Research Laboratories, Inc. Methods and systems for data-driven battery state of charge (SoC) estimation
US10099679B2 (en) 2016-10-20 2018-10-16 Fca Us Llc Battery state and parameter estimation using a mixed sigma-point kalman filtering and recursive least squares technique
JP6815299B2 (ja) * 2017-09-22 2021-01-20 株式会社神戸製鋼所 油圧システムのパラメータ推定方法
KR20220009258A (ko) * 2020-07-15 2022-01-24 주식회사 엘지에너지솔루션 배터리 관리 장치, 배터리 팩, 에너지 저장 시스템 및 배터리 관리 방법
IT202000025441A1 (it) 2020-10-28 2021-01-28 Nesyt Spin Off S R L Alimentatore intelligente per batterie

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4390841A (en) * 1980-10-14 1983-06-28 Purdue Research Foundation Monitoring apparatus and method for battery power supply
AU5555790A (en) * 1989-05-12 1990-11-29 Fraunhofer-Gesellschaft Zur Forderung Der Angewandten Forschung E.V. Process and device for determining physical variables of rechargeable electric accumulators
US5687077A (en) * 1991-07-31 1997-11-11 Universal Dynamics Limited Method and apparatus for adaptive control
US6353815B1 (en) * 1998-11-04 2002-03-05 The United States Of America As Represented By The United States Department Of Energy Statistically qualified neuro-analytic failure detection method and system
DE10106508A1 (de) * 2001-02-13 2002-08-29 Bosch Gmbh Robert Verfahren und Anordnung zur Bestimmung der Leistungsfähigkeit einer Batterie
US6441586B1 (en) * 2001-03-23 2002-08-27 General Motors Corporation State of charge prediction method and apparatus for a battery
EP1417503B1 (de) * 2001-06-29 2010-02-17 Robert Bosch Gmbh Verfahren zur ermittlung des ladezustands und/oder der leistungsfähigkeit eines ladungsspeichers
US6534954B1 (en) * 2002-01-10 2003-03-18 Compact Power Inc. Method and apparatus for a battery state of charge estimator
US20030184307A1 (en) * 2002-02-19 2003-10-02 Kozlowski James D. Model-based predictive diagnostic tool for primary and secondary batteries
US6927554B2 (en) * 2003-08-28 2005-08-09 General Motors Corporation Simple optimal estimator for PbA state of charge

Also Published As

Publication number Publication date
EP1989563A4 (en) 2012-12-19
TW200809242A (en) 2008-02-16
KR20080104296A (ko) 2008-12-02
EP1989563B1 (en) 2014-04-23
WO2007100189A1 (en) 2007-09-07
US7521895B2 (en) 2009-04-21
CN101395489A (zh) 2009-03-25
CN101395489B (zh) 2011-11-02
US20070236182A1 (en) 2007-10-11
WO2007100189A8 (en) 2008-09-12
JP2009528662A (ja) 2009-08-06
EP1989563A1 (en) 2008-11-12
KR101355959B1 (ko) 2014-02-03
TWI337663B (en) 2011-02-21

Similar Documents

Publication Publication Date Title
JP5036733B2 (ja) 推定バッテリー状態ベクトルと推定バッテリーパラメータベクトルとを定めるシステム及び方法
JP4920696B2 (ja) 推定バッテリーパラメータベクトルを定めるシステムと方法及びその製造物
JP4920703B2 (ja) 推定組合せバッテリー状態‐パラメータベクトルを定めるシステム、方法、及びその製造物
JP4920693B2 (ja) 推定バッテリー状態ベクトルを決めるシステムと方法及びその製造物
JP5058814B2 (ja) バッテリーの状態及びパラメーターの推定システム及び方法
KR101021745B1 (ko) 배터리 셀 모듈의 추정 상태를 결정하는 시스템, 방법 및 그 제조물
JP5259190B2 (ja) ジョイントバッテリー状態とパラメーター推定システム及び方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20091225

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120619

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120703

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150713

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5036733

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: R3D02

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250