JP2003075517A - 二次電池の充電率推定装置 - Google Patents
二次電池の充電率推定装置Info
- Publication number
- JP2003075517A JP2003075517A JP2001268314A JP2001268314A JP2003075517A JP 2003075517 A JP2003075517 A JP 2003075517A JP 2001268314 A JP2001268314 A JP 2001268314A JP 2001268314 A JP2001268314 A JP 2001268314A JP 2003075517 A JP2003075517 A JP 2003075517A
- Authority
- JP
- Japan
- Prior art keywords
- value
- terminal voltage
- equation
- charging rate
- open circuit
- 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.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Tests Of Electric Status Of Batteries (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Secondary Cells (AREA)
Abstract
(57)【要約】
【課題】充電率ＳＯＣを正確に推定できる二次電池の充
電率推定装置を提供する。 【解決手段】電池モデルを用いた適応デジタルフィルタ
演算で開路電圧を推定し、その結果を用いて充電率を演
算する充電率推定装置において、端子電圧検出値と端子
電圧初期値との差分値を求め、その端子電圧差分値と電
流検出値とを適応デジタルフィルタの入力としている。
端子電圧差分値は端子電圧検出値と端子電圧初期値との
差分値なので、演算開始時には、端子電圧検出値＝端子
電圧初期値のため、差分値＝０となる。したがって、適
応デジタルフィルタ内部の推定パラメータの初期状態が
約０である推定演算開始時点においては入力値が全て０
となるので、演算開始直後に推定パラメータ値が乱れて
真値に収束しないという問題点を解決することができ、
充電率ＳＯＣを正確に推定することが出来る。
電率推定装置を提供する。 【解決手段】電池モデルを用いた適応デジタルフィルタ
演算で開路電圧を推定し、その結果を用いて充電率を演
算する充電率推定装置において、端子電圧検出値と端子
電圧初期値との差分値を求め、その端子電圧差分値と電
流検出値とを適応デジタルフィルタの入力としている。
端子電圧差分値は端子電圧検出値と端子電圧初期値との
差分値なので、演算開始時には、端子電圧検出値＝端子
電圧初期値のため、差分値＝０となる。したがって、適
応デジタルフィルタ内部の推定パラメータの初期状態が
約０である推定演算開始時点においては入力値が全て０
となるので、演算開始直後に推定パラメータ値が乱れて
真値に収束しないという問題点を解決することができ、
充電率ＳＯＣを正確に推定することが出来る。
Description
【０００１】
【発明の属する技術分野】本発明は、二次電池の充電率
（ＳＯＣ）を推定する装置に関する。
（ＳＯＣ）を推定する装置に関する。
【０００２】
【従来の技術】二次電池の充電率ＳＯＣ（充電状態とも
言う）は開路電圧Ｖ０（通電遮断時の電池端子電圧であ
り、起電力、開放電圧とも言う）と相関があるので、開
路電圧Ｖ０を求めれば充電率を推定することが出来る。
しかし、二次電池の端子電圧は、通電を遮断（充放電を
終了）した後も安定するまでに時間を要するので、正確
な開路電圧Ｖ０を求めるには、充放電を終了してから所
定の時間が必要である。したがって充放電中や充放電直
後では、正確な開路電圧Ｖ０を求めることが出来ないの
で、上記の方法で充電率ＳＯＣを求めることが出来な
い。そのため、従来は、下記のような方法を用いて開路
電圧Ｖ０を推定している。
言う）は開路電圧Ｖ０（通電遮断時の電池端子電圧であ
り、起電力、開放電圧とも言う）と相関があるので、開
路電圧Ｖ０を求めれば充電率を推定することが出来る。
しかし、二次電池の端子電圧は、通電を遮断（充放電を
終了）した後も安定するまでに時間を要するので、正確
な開路電圧Ｖ０を求めるには、充放電を終了してから所
定の時間が必要である。したがって充放電中や充放電直
後では、正確な開路電圧Ｖ０を求めることが出来ないの
で、上記の方法で充電率ＳＯＣを求めることが出来な
い。そのため、従来は、下記のような方法を用いて開路
電圧Ｖ０を推定している。
【０００３】二次電池の充電率（ＳＯＣ）を推定する技
術に関する従来例としては、「論文“適応デジタルフィ
ルタを用いた鉛電池の開路電圧と残存容量の推定”四国
総研、四国電力、湯浅電池 T.IEEE Japan Vol.112-C,N
o.4 1992」に記載されたものがある。この電池状態検出
手法は、通電中の二次電池（鉛電池やリチウムイオン電
池等の充放電可能な電池）の端子電圧と電流の計測デー
タに、「適応デジタルフィルタ」を用いて開路電圧Ｖ０
を推定（パラメータ同定）して、この値から電池の充電
率ＳＯＣを推定するものである。上記の従来例において
は、「非回帰型の電池モデル（出力値が入力値の現在値
および過去値だけで決るモデル）」に相当する下記（数
４）式に、「適応デジタルフィルタ（逐次型のモデルパ
ラメータ同定アルゴリズム）」を用いて、下記（数４）
式中のパラメータの一つである開路電圧Ｃｊを算出し
て、この値から充電率ＳＯＣを算出している。
術に関する従来例としては、「論文“適応デジタルフィ
ルタを用いた鉛電池の開路電圧と残存容量の推定”四国
総研、四国電力、湯浅電池 T.IEEE Japan Vol.112-C,N
o.4 1992」に記載されたものがある。この電池状態検出
手法は、通電中の二次電池（鉛電池やリチウムイオン電
池等の充放電可能な電池）の端子電圧と電流の計測デー
タに、「適応デジタルフィルタ」を用いて開路電圧Ｖ０
を推定（パラメータ同定）して、この値から電池の充電
率ＳＯＣを推定するものである。上記の従来例において
は、「非回帰型の電池モデル（出力値が入力値の現在値
および過去値だけで決るモデル）」に相当する下記（数
４）式に、「適応デジタルフィルタ（逐次型のモデルパ
ラメータ同定アルゴリズム）」を用いて、下記（数４）
式中のパラメータの一つである開路電圧Ｃｊを算出し
て、この値から充電率ＳＯＣを算出している。
【０００４】
【数４】
ただし、Ｖｊ：端子電圧 Ｉｊ−ｋ：ｋサンプル周期
前の電流 Ｃｊ：開路電圧 ｂｋ，ｊ：適応デジタルフィルタの
係数 Ｎ：次数 Ｃｊ：開路電圧 そして、電池はＩｊ＝０の時（通電遮断時）、Ｖｊ＝Ｃ
ｊである。
前の電流 Ｃｊ：開路電圧 ｂｋ，ｊ：適応デジタルフィルタの
係数 Ｎ：次数 Ｃｊ：開路電圧 そして、電池はＩｊ＝０の時（通電遮断時）、Ｖｊ＝Ｃ
ｊである。
【０００５】
【発明が解決しようとする課題】上記のように、従来例
においては、推定パラメータの初期値が全て０である適
応デジタルフィルタを用いて演算しているが、推定演算
開始時の実際の開路電圧等は０ではない。つまり、推定
演算開始時に、適応デジタルフィルタ内部の推定パラメ
ータの初期状態が約０であるにもかかわらず、Ｖ０＝Ｃ
０≠０という値が入力されるので、推定開始直後に推定
パラメータ値が乱れて真値に収束しない、という問題が
あった。
においては、推定パラメータの初期値が全て０である適
応デジタルフィルタを用いて演算しているが、推定演算
開始時の実際の開路電圧等は０ではない。つまり、推定
演算開始時に、適応デジタルフィルタ内部の推定パラメ
ータの初期状態が約０であるにもかかわらず、Ｖ０＝Ｃ
０≠０という値が入力されるので、推定開始直後に推定
パラメータ値が乱れて真値に収束しない、という問題が
あった。
【０００６】本発明は上記のごとき従来技術の問題を解
決するためになされたものであり、充電率ＳＯＣおよび
その他のパラメータを正確に推定することの出来る二次
電池の充電率推定装置を提供することを目的とする。
決するためになされたものであり、充電率ＳＯＣおよび
その他のパラメータを正確に推定することの出来る二次
電池の充電率推定装置を提供することを目的とする。
【０００７】
【課題を解決するための手段】上記の目的を達成するた
め本発明においては、特許請求の範囲に記載するように
構成している。すなわち、請求項１においては、電流遮
断時である推定演算開始時点における端子電圧初期値を
演算し、端子電圧検出値と端子電圧初期値の差分値を演
算し、電流計測値と前記差分値を、電池モデルを用いた
適応デジタルフィルタに入力し、推定演算された開路電
圧の変化分に端子電圧初期値を加算して開路電圧Ｖ０の
推定値を演算し、予め求めた開路電圧Ｖ０と充電率ＳＯ
Ｃとの関係に基づいて上記の推定した開路電圧Ｖ０を用
いて充電率を推定するように構成している。なお、請求
項１の構成において、適応デジタルフィルタとしては、
前記従来例に記載のごときものを用いることも出来る
が、請求項２または請求項３に記載のごとき構成に適用
すれば、さらに正確な値を求めることが出来る。
め本発明においては、特許請求の範囲に記載するように
構成している。すなわち、請求項１においては、電流遮
断時である推定演算開始時点における端子電圧初期値を
演算し、端子電圧検出値と端子電圧初期値の差分値を演
算し、電流計測値と前記差分値を、電池モデルを用いた
適応デジタルフィルタに入力し、推定演算された開路電
圧の変化分に端子電圧初期値を加算して開路電圧Ｖ０の
推定値を演算し、予め求めた開路電圧Ｖ０と充電率ＳＯ
Ｃとの関係に基づいて上記の推定した開路電圧Ｖ０を用
いて充電率を推定するように構成している。なお、請求
項１の構成において、適応デジタルフィルタとしては、
前記従来例に記載のごときものを用いることも出来る
が、請求項２または請求項３に記載のごとき構成に適用
すれば、さらに正確な値を求めることが出来る。
【０００８】また、請求項２においては、請求項１の演
算に用いる適応デジタルフィルタとして、（数１）式に
示す連続時間系の電池モデルを用いて適応デジタルフィ
ルタ演算を行い、（数１）式中のオフセット項である開
路電圧Ｖ０および過渡項であるＡ(ｓ）、Ｂ(ｓ）に対応
するパラメータを一括推定するように構成している。
算に用いる適応デジタルフィルタとして、（数１）式に
示す連続時間系の電池モデルを用いて適応デジタルフィ
ルタ演算を行い、（数１）式中のオフセット項である開
路電圧Ｖ０および過渡項であるＡ(ｓ）、Ｂ(ｓ）に対応
するパラメータを一括推定するように構成している。
【０００９】また、請求項３においては、前記（数１）
式に示した連続時間系の電池モデルを、離散時間系表現
に変換した（数２）式で示される「自己回帰型モデル」
を用いて、適応デジタルフィルタ演算を行い、各パラメ
ータ（ａｋ、ｂｋ、ｃｋ、ｒ ｋ）を一括推定し、さら
に、（数３）式を用いて開路電圧Ｖ０を推定演算するよ
うに構成している。上記請求項１および請求項２の構成
は、例えば後記実施例１に相当する。
式に示した連続時間系の電池モデルを、離散時間系表現
に変換した（数２）式で示される「自己回帰型モデル」
を用いて、適応デジタルフィルタ演算を行い、各パラメ
ータ（ａｋ、ｂｋ、ｃｋ、ｒ ｋ）を一括推定し、さら
に、（数３）式を用いて開路電圧Ｖ０を推定演算するよ
うに構成している。上記請求項１および請求項２の構成
は、例えば後記実施例１に相当する。
【００１０】
【発明の効果】請求項１に記載の発明においては、検出
した電流と端子電圧差分値を適応デジタルフィルタに入
力している。端子電圧差分値は端子電圧検出値と端子電
圧初期値との差分値なので、演算開始時には、端子電圧
検出値＝端子電圧初期値のため、差分値＝０となる。し
たがって、適応デジタルフィルタ内部の推定パラメータ
の初期状態が約０である推定演算開始時点においては入
力値が全て０となるので、演算開始直後に推定パラメー
タ値が乱れて真値に収束しないという問題点を解決する
ことができる。
した電流と端子電圧差分値を適応デジタルフィルタに入
力している。端子電圧差分値は端子電圧検出値と端子電
圧初期値との差分値なので、演算開始時には、端子電圧
検出値＝端子電圧初期値のため、差分値＝０となる。し
たがって、適応デジタルフィルタ内部の推定パラメータ
の初期状態が約０である推定演算開始時点においては入
力値が全て０となるので、演算開始直後に推定パラメー
タ値が乱れて真値に収束しないという問題点を解決する
ことができる。
【００１１】請求項２に記載の発明においては、二次電
池の電流Ｉと端子電圧Ｖと開路電圧Ｖ０の関係を、（数
１）式に示す伝達関数で近似して、開路電圧Ｖ０の項を
定常項（オフセット項）とみなす（定式化する）こと
で、「最小二乗法」等の「適応デジタルフィルタ」（公
知の推定アルゴリズム）を適用することが可能となる。
その結果、(数１）式中のパラメータ（オフセット項で
あるＶ０、および過渡項であるＡ(ｓ）やＢ(ｓ）に対応
するパラメータ）を一括推定することが可能となる。こ
れらパラメータは、充電率ＳＯＣや温度や劣化度などに
影響され、時々刻々変化するものであるが、適応デジタ
ルフィルタにより精度良く逐次推定することが可能であ
る。なお、例えば後記図５に示す開路電圧Ｖ０と充電率
ＳＯＣの関係は、温度や電池の劣化度に影響されにくく
一定の相関関係にあるので、この特性を予め記憶してお
けば、開路電圧Ｖ０から充電率ＳＯＣが直接算出でき
る。したがって、充電率ＳＯＣについても開路電圧Ｖ０
と同様に、条件によらず正確な推定が可能である、とい
う効果がある。
池の電流Ｉと端子電圧Ｖと開路電圧Ｖ０の関係を、（数
１）式に示す伝達関数で近似して、開路電圧Ｖ０の項を
定常項（オフセット項）とみなす（定式化する）こと
で、「最小二乗法」等の「適応デジタルフィルタ」（公
知の推定アルゴリズム）を適用することが可能となる。
その結果、(数１）式中のパラメータ（オフセット項で
あるＶ０、および過渡項であるＡ(ｓ）やＢ(ｓ）に対応
するパラメータ）を一括推定することが可能となる。こ
れらパラメータは、充電率ＳＯＣや温度や劣化度などに
影響され、時々刻々変化するものであるが、適応デジタ
ルフィルタにより精度良く逐次推定することが可能であ
る。なお、例えば後記図５に示す開路電圧Ｖ０と充電率
ＳＯＣの関係は、温度や電池の劣化度に影響されにくく
一定の相関関係にあるので、この特性を予め記憶してお
けば、開路電圧Ｖ０から充電率ＳＯＣが直接算出でき
る。したがって、充電率ＳＯＣについても開路電圧Ｖ０
と同様に、条件によらず正確な推定が可能である、とい
う効果がある。
【００１２】請求項３に記載の発明においては、具体的
な方法の一つとして、連続時間系で記述された（数１）
式の等価回路モデルを、離散時間系に変換する際に、オ
フセット項である開路電圧Ｖ０を（数２）式の右辺第３
項のように展開することで、拡大最小二乗法などの適応
デジタルフィルタを適用することが可能となる。その結
果、マイコン等による演算で推定処理が実行可能にな
る、という効果がある。
な方法の一つとして、連続時間系で記述された（数１）
式の等価回路モデルを、離散時間系に変換する際に、オ
フセット項である開路電圧Ｖ０を（数２）式の右辺第３
項のように展開することで、拡大最小二乗法などの適応
デジタルフィルタを適用することが可能となる。その結
果、マイコン等による演算で推定処理が実行可能にな
る、という効果がある。
【００１３】
【発明の実施の形態】（実施例１）図１は、本発明の実
施例を機能ブロックで表した図である。図１において、
１は二次電池の電流を検出する電流検出手段、２は二次
電池の端子電圧を検出する端子電圧検出手段、３は演算
開始時における端子電圧の初期値を演算する端子電圧初
期値演算手段、４は端子電圧検出手段２で検出した端子
電圧検出値と端子電圧初期値演算手段３で求めた端子電
圧初期値との差分値を演算する端子電圧差分値演算手
段、５は電流検出値と端子電圧差分値とを入力し、適応
デジタルフィルタによって開路電圧を推定する開路電圧
推定手段、６は開路電圧推定手段５で求めた開路電圧値
を用いて、例えば後記図５の特性から充電率ＳＯＣを求
める充電率推定手段である。
施例を機能ブロックで表した図である。図１において、
１は二次電池の電流を検出する電流検出手段、２は二次
電池の端子電圧を検出する端子電圧検出手段、３は演算
開始時における端子電圧の初期値を演算する端子電圧初
期値演算手段、４は端子電圧検出手段２で検出した端子
電圧検出値と端子電圧初期値演算手段３で求めた端子電
圧初期値との差分値を演算する端子電圧差分値演算手
段、５は電流検出値と端子電圧差分値とを入力し、適応
デジタルフィルタによって開路電圧を推定する開路電圧
推定手段、６は開路電圧推定手段５で求めた開路電圧値
を用いて、例えば後記図５の特性から充電率ＳＯＣを求
める充電率推定手段である。
【００１４】図２は、実施例の具体的な構成を示すブロ
ック図である。この実施例は、二次電池でモータ等の負
荷を駆動したり、モータの回生電力で二次電池を充電す
るシステムに、二次電池の充電率推定装置を設けた例を
示す。図２において、１０は二次電池（単に電池とも言
う）、２０はモータ等の負荷、３０は電池の充電状態を
推定する電子制御ユニットで、プログラムを演算するＣ
ＰＵやプログラムを記憶したＲＯＭや演算結果を記憶す
るＲＡＭから成るマイクロコンピュータと電子回路等で
構成される。４０は電池から充放電される電流を検出す
る電流計、５０は電池の端子電圧を検出する電圧計であ
り、それぞれ電子制御ユニット３０に接続される。上記
の電子制御ユニット３０は前記図１の端子電圧初期値演
算手段３、端子電圧差分値演算手段４、開路電圧演算手
段５および充電率推定手段６の部分に相当する。また、
電流計４０は電流検出手段１に、電圧計５０は端子電圧
検出手段に、それぞれ相当する。
ック図である。この実施例は、二次電池でモータ等の負
荷を駆動したり、モータの回生電力で二次電池を充電す
るシステムに、二次電池の充電率推定装置を設けた例を
示す。図２において、１０は二次電池（単に電池とも言
う）、２０はモータ等の負荷、３０は電池の充電状態を
推定する電子制御ユニットで、プログラムを演算するＣ
ＰＵやプログラムを記憶したＲＯＭや演算結果を記憶す
るＲＡＭから成るマイクロコンピュータと電子回路等で
構成される。４０は電池から充放電される電流を検出す
る電流計、５０は電池の端子電圧を検出する電圧計であ
り、それぞれ電子制御ユニット３０に接続される。上記
の電子制御ユニット３０は前記図１の端子電圧初期値演
算手段３、端子電圧差分値演算手段４、開路電圧演算手
段５および充電率推定手段６の部分に相当する。また、
電流計４０は電流検出手段１に、電圧計５０は端子電圧
検出手段に、それぞれ相当する。
【００１５】まず、本実施例で用いる「電池モデル」を
説明する。図３は、二次電池の等価回路モデルを示す図
であり、下記（数５）式で示される。（数５）式におい
て、モデル入力は電流Ｉ［Ａ］（正値は充電、負値は放
電）、モデル出力は端子電圧Ｖ［Ｖ］、Ｒ１〔Ω］は電
荷移動抵抗、Ｒ２［Ω］は純抵抗、Ｃ１［Ｆ］は電気二
重層容量、Ｖ０［Ｖ］は開路電圧である。なお、ｓはラ
プラス演算子である。本モデルは、正極、負極を特に分
離していないリダクションモデル（一次）であるが、実
際の電池の充放電特性を比較的正確に示すことが可能で
ある。
説明する。図３は、二次電池の等価回路モデルを示す図
であり、下記（数５）式で示される。（数５）式におい
て、モデル入力は電流Ｉ［Ａ］（正値は充電、負値は放
電）、モデル出力は端子電圧Ｖ［Ｖ］、Ｒ１〔Ω］は電
荷移動抵抗、Ｒ２［Ω］は純抵抗、Ｃ１［Ｆ］は電気二
重層容量、Ｖ０［Ｖ］は開路電圧である。なお、ｓはラ
プラス演算子である。本モデルは、正極、負極を特に分
離していないリダクションモデル（一次）であるが、実
際の電池の充放電特性を比較的正確に示すことが可能で
ある。
【００１６】
【数５】
上記（数５）式に零次ホールドを付加してＺ変換するこ
とで、線形離散時間システム（数６）式を得る。
とで、線形離散時間システム（数６）式を得る。
【００１７】
【数６】
ただし、ａ１、ｂ０、ｂ１は定数、ｚ−１は遅延演算子
Ｖ(ｋ）は現時点の端子電圧、Ｖ(ｋ−ｎ）はｎサンプル
周期前の端子電圧を示すものであり、以後、他の変数も
これに準ずる。
周期前の端子電圧を示すものであり、以後、他の変数も
これに準ずる。
【００１８】次に、「拡大最小二乗法」と呼ばれる公知
の「適応デジタルフィルタ（逐次型同定アルゴリズ
ム）」を一般形でまず説明する。線形離散時間システム
で記述されるプラントモデルを（数７）式とする。
の「適応デジタルフィルタ（逐次型同定アルゴリズ
ム）」を一般形でまず説明する。線形離散時間システム
で記述されるプラントモデルを（数７）式とする。
【００１９】
【数７】
ただし、Ａ(ｚ−１）、Ｂ(ｚ−１）はｚ−１の多項式
ｎはＡ(ｚ−１）、Ｂ(ｚ−１）の次数
ｙ(ｋ）は出力、ｕ(ｋ）は入力、ｒ(ｋ）は式誤差（雑
音） 上記のｒ(ｋ）は一般に白色雑音ではなく、通常の「最
小二乗法」を用いると推定値に偏り（バイアス）が生じ
る（電池の場合、開路電圧が数７式の右辺第２項の定常
値に相当）。この問題に対応する幾つかの改良手法の一
つとして、「拡大最小二乗法」がある。この手法では、
（数８）式のように式誤差を改めて定義する。
音） 上記のｒ(ｋ）は一般に白色雑音ではなく、通常の「最
小二乗法」を用いると推定値に偏り（バイアス）が生じ
る（電池の場合、開路電圧が数７式の右辺第２項の定常
値に相当）。この問題に対応する幾つかの改良手法の一
つとして、「拡大最小二乗法」がある。この手法では、
（数８）式のように式誤差を改めて定義する。
【００２０】
【数８】
ただし、ｅ(ｋ）は平均値零で分散σ２の白色雑音
Ｃ(ｚ−１）はｚ−１の多項式
ｐはＣ(ｚ−１）の次数
上記の（数７）式、（数８）式を用いて下記（数９）式
が求まる。
が求まる。
【００２１】
【数９】
更に（数１１）式の定義を用いて（数９）式を（数１
０）式に変形する。
０）式に変形する。
【００２２】
【数１０】
なお、（数１０）式中のＴは行列の配置を示す。
【００２３】
【数１１】
ｋ時点のパラメータ推定値をθ＊(ｋ）とすると、逐次
推定アルゴリズムは（数１２）式に示すようになる。
推定アルゴリズムは（数１２）式に示すようになる。
【００２４】ただし、（数１１）式のｎ＝１、ｐ＝１。
【００２５】
【数１２】
次に、前記（数６）式の電池モデルに以上の適応フィル
タアルゴリズムを適用する。（数１０）式において下記
（数１３）式のように定義すると、（数６）式と（数
９）式のモデルが一致するので前述の同定アルゴリズム
（数１２）式を適用することができ、パラメータθ
(ｋ）を推定することができる。
タアルゴリズムを適用する。（数１０）式において下記
（数１３）式のように定義すると、（数６）式と（数
９）式のモデルが一致するので前述の同定アルゴリズム
（数１２）式を適用することができ、パラメータθ
(ｋ）を推定することができる。
【００２６】
【数１３】
なお、ｒ＊(ｋ）の導出には下記（数１４）式を用い
る。
る。
【００２７】
【数１４】
上記のように「適応フィルタ」で同定されたパラメータ
θ＊(ｋ）から下記（数１５）式に示すように、電池の
重要パラメータである開路電圧Ｖ０、内部抵抗Ｒ１＋Ｒ
２、時定数Ｒ１・Ｃ１が求まる。
θ＊(ｋ）から下記（数１５）式に示すように、電池の
重要パラメータである開路電圧Ｖ０、内部抵抗Ｒ１＋Ｒ
２、時定数Ｒ１・Ｃ１が求まる。
【００２８】
【数１５】
図４は、電子制御ユニット３０のマイクロコンピュータ
が行う処理のフローチャートである。図４のルーチンは
一定周期Ｔ０毎に実施される。なお、（ｋ）は今回の
値、（ｋ−ｎ）はｎサンプル周期前（ｎ周期前の演算）
の値を意味する。
が行う処理のフローチャートである。図４のルーチンは
一定周期Ｔ０毎に実施される。なお、（ｋ）は今回の
値、（ｋ−ｎ）はｎサンプル周期前（ｎ周期前の演算）
の値を意味する。
【００２９】ステップＳ１０では、電流Ｉ(ｋ）、端子
電圧Ｖ(ｋ）を計測する。ステップＳ２０では、推定演
算の開始を判断する。電子制御ユニット３０は二次電池
の動作を全て管理しているので、電源投入直後は電流Ｉ
＝０であるから開始可能（ＹＥＳ）として、唯１回のみ
ステップＳ３０へ進む。初期値演算以降はステップＳ４
０へ進む。
電圧Ｖ(ｋ）を計測する。ステップＳ２０では、推定演
算の開始を判断する。電子制御ユニット３０は二次電池
の動作を全て管理しているので、電源投入直後は電流Ｉ
＝０であるから開始可能（ＹＥＳ）として、唯１回のみ
ステップＳ３０へ進む。初期値演算以降はステップＳ４
０へ進む。
【００３０】ステップＳ３０では、適応デジタルフィル
タ（同定アルゴリズムとも言う）の演算開始時点におけ
る端子電圧初期値Ｖ ｉｎｉを算出し、ステップＳ４０
へ進む。演算開始時は電流Ｉ＝０であるから、開路電圧
初期値Ｖ０ ｉｎｉは端子電圧初期値Ｖ ｉｎｉに等し
い。 Ｖ ｉｎｉ＝Ｖ(０)， Ｖ０ ｉｎｉ＝Ｖ(０) なお、Ｖ(０)は今回の端子電圧を示す。
タ（同定アルゴリズムとも言う）の演算開始時点におけ
る端子電圧初期値Ｖ ｉｎｉを算出し、ステップＳ４０
へ進む。演算開始時は電流Ｉ＝０であるから、開路電圧
初期値Ｖ０ ｉｎｉは端子電圧初期値Ｖ ｉｎｉに等し
い。 Ｖ ｉｎｉ＝Ｖ(０)， Ｖ０ ｉｎｉ＝Ｖ(０) なお、Ｖ(０)は今回の端子電圧を示す。
【００３１】ステップＳ４０では、検出した端子電圧Ｖ
(ｋ）と上記の端子電圧初期値Ｖ ｉｎｉとの差分値であ
る端子電圧差分値△Ｖ(ｋ)を算出する。 △Ｖ(ｋ)＝Ｖ(ｋ)−Ｖ ｉｎｉ なお、演算開始時は、Ｖ(ｋ)＝Ｖ ｉｎｉであるから、
△Ｖ(ｋ)＝０である。
(ｋ）と上記の端子電圧初期値Ｖ ｉｎｉとの差分値であ
る端子電圧差分値△Ｖ(ｋ)を算出する。 △Ｖ(ｋ)＝Ｖ(ｋ)−Ｖ ｉｎｉ なお、演算開始時は、Ｖ(ｋ)＝Ｖ ｉｎｉであるから、
△Ｖ(ｋ)＝０である。
【００３２】ステップＳ５０では、ノイズ除去の為に、
電流Ｉ(ｋ）および端子電圧差分値△Ｖ(ｋ）にローパス
フィルタ処理をし、処理後に改めて電流Ｉ(ｋ）および
端子電圧差分値△Ｖ(ｋ）とする。
電流Ｉ(ｋ）および端子電圧差分値△Ｖ(ｋ）にローパス
フィルタ処理をし、処理後に改めて電流Ｉ(ｋ）および
端子電圧差分値△Ｖ(ｋ）とする。
【００３３】ステップＳ６０では、適応デジタルフィル
タに、検出した電流値Ｉ(ｋ)と上記の求めた端子電圧差
分値△Ｖ(ｋ)とを入力して（数１３）式を行う。ここで
入力として端子電圧差分値を用いるのは、推定演算開始
時に適応デジタルフィルタ内の推定パラメータの初期値
を約０としているので、推定パラメータが発散しないよ
うに、入力を全て０とするためである。従来例では推定
演算開始時に電流Ｉ(０)＝０と端子電圧Ｖ(０)＝Ｖ ｉ
ｎｉ≠０を入力するので、推定パラメータが発散してし
まう。なお、電源投入直後から適応フィルタ演算を行う
ので、電流値Ｉ(ｋ）はゼロを初期値として入力され
る。したがって、電流値Ｉ(ｋ）は、必ずしも偏差を用
いる必要性はない。 ｕ(ｋ)＝△Ｉ(ｋ)、 ｙ(ｋ)＝△Ｖ(ｋ) θＴ(ｋ)＝{−ａ１(ｋ)，ｂ０(ｋ)，ｂ１(ｋ)，−ｃ
１(ｋ)} ωＴ(ｋ)＝{ｙ(ｋ−１)，ｕ(ｋ)，ｕ(ｋ−１)，ｒ(ｋ−
１)} なお、ｒ(ｋ−１)は（数１４）式で算出する。
タに、検出した電流値Ｉ(ｋ)と上記の求めた端子電圧差
分値△Ｖ(ｋ)とを入力して（数１３）式を行う。ここで
入力として端子電圧差分値を用いるのは、推定演算開始
時に適応デジタルフィルタ内の推定パラメータの初期値
を約０としているので、推定パラメータが発散しないよ
うに、入力を全て０とするためである。従来例では推定
演算開始時に電流Ｉ(０)＝０と端子電圧Ｖ(０)＝Ｖ ｉ
ｎｉ≠０を入力するので、推定パラメータが発散してし
まう。なお、電源投入直後から適応フィルタ演算を行う
ので、電流値Ｉ(ｋ）はゼロを初期値として入力され
る。したがって、電流値Ｉ(ｋ）は、必ずしも偏差を用
いる必要性はない。 ｕ(ｋ)＝△Ｉ(ｋ)、 ｙ(ｋ)＝△Ｖ(ｋ) θＴ(ｋ)＝{−ａ１(ｋ)，ｂ０(ｋ)，ｂ１(ｋ)，−ｃ
１(ｋ)} ωＴ(ｋ)＝{ｙ(ｋ−１)，ｕ(ｋ)，ｕ(ｋ−１)，ｒ(ｋ−
１)} なお、ｒ(ｋ−１)は（数１４）式で算出する。
【００３４】ステップＳ７０では、（数１２−１）式、
（数１２−２）式、（数１４）式を演算する。ステップ
Ｓ８０では、（数１５）式を演算し、下記のように、△
Ｖ０(ｋ）、Ｒ１＋Ｒ２、Ｒ１・Ｃ１、△Ｖ(ｋ)を求め
る。ただし、△Ｖ０(ｋ）は推定演算開始時からの開路
電圧推定値の変化分、△Ｖ(ｋ）は推定演算開始時から
の端子電圧推定値の変化分である。 △Ｖ０(ｋ)＝−ｃ１(ｋ)・ｒ(ｋ−１)／{１＋ａ１(ｋ)} Ｒ１＋Ｒ２＝{ｂ０(ｋ)＋ｂ１(ｋ)}／{１＋ａ１(ｋ)} Ｒ１・Ｃ１＝−Ｔ０／Ｉｎ{−ａ１(ｋ)} △Ｖ(ｋ)＝−ａ１(ｋ−１)・ｙ(ｋ−１)＋ｂ０(ｋ−１)
・ｕ(ｋ)＋ｂ１(ｋ−１)・ｕ(ｋ−１)−ｃ１(ｋ−１)・
ｒ(ｋ−１) ステップＳ９０では、下記のように、開路電圧推定値Ｖ
０(ｋ）と端子電圧推定値Ｖ(ｋ）を算出する。つまり、
ステップＳ８０で算出した△Ｖ０(ｋ）は同定アルゴリ
ズム開始時からの開路電圧の変化分であるから、それに
開路電圧初期値Ｖ０ ｉｎｉを加算して開路電圧推定値
Ｖ０(ｋ）を算出する。また、△Ｖ(ｋ）は同定アルゴリ
ズム開始時からの端子電圧の変化分であるから、それに
端子電圧初期値Ｖ ｉｎｉを加算して端子電圧推定値Ｖ
(ｋ）を算出する。 Ｖ０(ｋ)＝△Ｖ０(ｋ)＋Ｖ０ ｉｎｉ Ｖ(ｋ)＝△Ｖ(ｋ)＋Ｖ ｉｎｉ ステップＳ１００では、図５に示すような開路電圧Ｖ０
と充電率ＳＯＣの相関マップを用いて、上記で算出した
Ｖ０(ｋ）から充電率ＳＯＣ(ｋ）を算出する。なお、図
５中のＶＬはＳＯＣ＝０％に、ＶＨはＳＯＣ＝１００％
に相当する開路電圧である。ステップＳ１１０では、次
回の演算用に各数値を保存して一連の処理を終了する。
（数１２−２）式、（数１４）式を演算する。ステップ
Ｓ８０では、（数１５）式を演算し、下記のように、△
Ｖ０(ｋ）、Ｒ１＋Ｒ２、Ｒ１・Ｃ１、△Ｖ(ｋ)を求め
る。ただし、△Ｖ０(ｋ）は推定演算開始時からの開路
電圧推定値の変化分、△Ｖ(ｋ）は推定演算開始時から
の端子電圧推定値の変化分である。 △Ｖ０(ｋ)＝−ｃ１(ｋ)・ｒ(ｋ−１)／{１＋ａ１(ｋ)} Ｒ１＋Ｒ２＝{ｂ０(ｋ)＋ｂ１(ｋ)}／{１＋ａ１(ｋ)} Ｒ１・Ｃ１＝−Ｔ０／Ｉｎ{−ａ１(ｋ)} △Ｖ(ｋ)＝−ａ１(ｋ−１)・ｙ(ｋ−１)＋ｂ０(ｋ−１)
・ｕ(ｋ)＋ｂ１(ｋ−１)・ｕ(ｋ−１)−ｃ１(ｋ−１)・
ｒ(ｋ−１) ステップＳ９０では、下記のように、開路電圧推定値Ｖ
０(ｋ）と端子電圧推定値Ｖ(ｋ）を算出する。つまり、
ステップＳ８０で算出した△Ｖ０(ｋ）は同定アルゴリ
ズム開始時からの開路電圧の変化分であるから、それに
開路電圧初期値Ｖ０ ｉｎｉを加算して開路電圧推定値
Ｖ０(ｋ）を算出する。また、△Ｖ(ｋ）は同定アルゴリ
ズム開始時からの端子電圧の変化分であるから、それに
端子電圧初期値Ｖ ｉｎｉを加算して端子電圧推定値Ｖ
(ｋ）を算出する。 Ｖ０(ｋ)＝△Ｖ０(ｋ)＋Ｖ０ ｉｎｉ Ｖ(ｋ)＝△Ｖ(ｋ)＋Ｖ ｉｎｉ ステップＳ１００では、図５に示すような開路電圧Ｖ０
と充電率ＳＯＣの相関マップを用いて、上記で算出した
Ｖ０(ｋ）から充電率ＳＯＣ(ｋ）を算出する。なお、図
５中のＶＬはＳＯＣ＝０％に、ＶＨはＳＯＣ＝１００％
に相当する開路電圧である。ステップＳ１１０では、次
回の演算用に各数値を保存して一連の処理を終了する。
【００３５】上記のように、本発明においては、ステッ
プＳ４０で、端子電圧検出値と端子電圧初期値との差分
値を求め、その端子電圧差分値と電流検出値とを適応デ
ジタルフィルタの入力としている。端子電圧差分値は端
子電圧検出値と端子電圧初期値との差分値なので、演算
開始時には、端子電圧検出値＝端子電圧初期値のため、
差分値＝０となる。したがって、適応デジタルフィルタ
内部の推定パラメータの初期状態が約０である推定演算
開始時点においては入力値が全て０となるので、演算開
始直後に推定パラメータ値が乱れて真値に収束しないと
いう問題点を解決することができる。
プＳ４０で、端子電圧検出値と端子電圧初期値との差分
値を求め、その端子電圧差分値と電流検出値とを適応デ
ジタルフィルタの入力としている。端子電圧差分値は端
子電圧検出値と端子電圧初期値との差分値なので、演算
開始時には、端子電圧検出値＝端子電圧初期値のため、
差分値＝０となる。したがって、適応デジタルフィルタ
内部の推定パラメータの初期状態が約０である推定演算
開始時点においては入力値が全て０となるので、演算開
始直後に推定パラメータ値が乱れて真値に収束しないと
いう問題点を解決することができる。
【００３６】以下、シミュレーションで求めた本発明の
効果を説明する。図６および図７は、電池モデルに真値
のパラメータを与えて、パラメータ同定したシミュレー
ション結果を示す図であり、図６は従来例、図７は本発
明の結果を示す。なお、図６、図７において、点線また
は破線は真値、実線は推定値を示す。
効果を説明する。図６および図７は、電池モデルに真値
のパラメータを与えて、パラメータ同定したシミュレー
ション結果を示す図であり、図６は従来例、図７は本発
明の結果を示す。なお、図６、図７において、点線また
は破線は真値、実線は推定値を示す。
【００３７】従来例では、図６のに示すように、推定
演算開始直後に端子電圧推定値が乱れる。また、図６の
に示すように、その後の内部抵抗の推定値も真値でな
く負の値に収束している。そのため、図６のに示すよ
うに、開路電圧推定値は真値とはかけ離れた値になって
いる。それに対して実施例では、図７のに示すよう
に、推定演算期始直後に端子電圧推定値が乱れないの
で、図７のに示すように、その後の内部抵抗の推定値
も真値に収束する。そのため、図７のに示すように、
開路電圧推定値は真値と滑らかに一致している。
演算開始直後に端子電圧推定値が乱れる。また、図６の
に示すように、その後の内部抵抗の推定値も真値でな
く負の値に収束している。そのため、図６のに示すよ
うに、開路電圧推定値は真値とはかけ離れた値になって
いる。それに対して実施例では、図７のに示すよう
に、推定演算期始直後に端子電圧推定値が乱れないの
で、図７のに示すように、その後の内部抵抗の推定値
も真値に収束する。そのため、図７のに示すように、
開路電圧推定値は真値と滑らかに一致している。
【００３８】また、本発明においては、二次電池の電流
Ｉと端子電圧Ｖと開路電圧Ｖ０の関係を、（数１）式に
示す伝達関数で近似して、開路電圧Ｖ０の項を定常項
（オフセット項）とみなす（定式化する）ことで、「最
小二乗法」等の「適応デジタルフィルタ」（公知の推定
アルゴリズム）を適用することが可能となる。その結
果、(数１）式中のパラメータ（オフセット項であるＶ
０、および過渡項であるＡ(ｓ）やＢ(ｓ）に対応するパ
ラメータ）を一括推定することが可能となる。これらパ
ラメータは、充電率ＳＯＣや温度や劣化度などに影響さ
れ、時々刻々変化するものであるが、適応デジタルフィ
ルタにより精度良く逐次推定することが可能である。な
お、例えば後記図５に示す開路電圧Ｖ０と充電率ＳＯＣ
の関係は、温度や電池の劣化度に影響されにくく一定の
相関関係にあるので、この特性を予め記憶しておけば、
開路電圧Ｖ０から充電率ＳＯＣが直接算出できる。した
がって、充電率ＳＯＣについても開路電圧Ｖ０と同様
に、条件によらず正確な推定が可能である。
Ｉと端子電圧Ｖと開路電圧Ｖ０の関係を、（数１）式に
示す伝達関数で近似して、開路電圧Ｖ０の項を定常項
（オフセット項）とみなす（定式化する）ことで、「最
小二乗法」等の「適応デジタルフィルタ」（公知の推定
アルゴリズム）を適用することが可能となる。その結
果、(数１）式中のパラメータ（オフセット項であるＶ
０、および過渡項であるＡ(ｓ）やＢ(ｓ）に対応するパ
ラメータ）を一括推定することが可能となる。これらパ
ラメータは、充電率ＳＯＣや温度や劣化度などに影響さ
れ、時々刻々変化するものであるが、適応デジタルフィ
ルタにより精度良く逐次推定することが可能である。な
お、例えば後記図５に示す開路電圧Ｖ０と充電率ＳＯＣ
の関係は、温度や電池の劣化度に影響されにくく一定の
相関関係にあるので、この特性を予め記憶しておけば、
開路電圧Ｖ０から充電率ＳＯＣが直接算出できる。した
がって、充電率ＳＯＣについても開路電圧Ｖ０と同様
に、条件によらず正確な推定が可能である。
【００３９】また、具体的な方法の一つとして、連続時
間系で記述された（数１）式の等価回路モデルを、離散
時間系に変換する際に、オフセット項である開路電圧Ｖ
０を（数２）式の右辺第３項のように展開することで、
拡大最小二乗法などの適応デジタルフィルタを適用する
ことが可能となる。その結果、マイコン等による演算で
推定処理が実行可能になる。
間系で記述された（数１）式の等価回路モデルを、離散
時間系に変換する際に、オフセット項である開路電圧Ｖ
０を（数２）式の右辺第３項のように展開することで、
拡大最小二乗法などの適応デジタルフィルタを適用する
ことが可能となる。その結果、マイコン等による演算で
推定処理が実行可能になる。
【図１】本発明の実施例の基本構成を機能ブロックで表
した図。
した図。
【図２】本発明の具体的構成を示すブロック図。
【図３】二次電池の等価回路モデルを示す図。
【図４】電子制御ユニット３０のマイクロコンピュータ
が行う処理のフローチャート。
が行う処理のフローチャート。
【図５】開路電圧Ｖ０と充電率ＳＯＣの相関マップ。
【図６】従来例における各数値の推定シミュレーション
結果の一例を示す図。
結果の一例を示す図。
【図７】実施例における各数値の推定シミュレーション
結果の一例を示す図。
結果の一例を示す図。
【符号の説明】
１…電流検出手段 ２…端子電圧検出手
段 ３…端子電圧初期値演算手段 ４…端子電圧差分値
演算手段 ５…開路電圧推定手段 ６…充電率推定手段 １０…二次電池 ２０…モータ等の負
荷 ３０…電子制御ユニット ４０…電流計 ５０…電圧計
段 ３…端子電圧初期値演算手段 ４…端子電圧差分値
演算手段 ５…開路電圧推定手段 ６…充電率推定手段 １０…二次電池 ２０…モータ等の負
荷 ３０…電子制御ユニット ４０…電流計 ５０…電圧計
─────────────────────────────────────────────────────
フロントページの続き
Ｆターム(参考） 2G016 CB06 CB11 CB13 CC04 CC24
CC27 CC28 CD18
5G003 BA01 DA02 EA05 GC05
5H030 AA03 AA06 AS08 FF42 FF44
Claims (3)
- 【請求項１】二次電池の電流を検出する手段と、 二次電池の端子電圧を検出する手段と、 電流遮断時である推定演算開始時点における端子電圧初
期値を演算する手段と、 前記端子電圧検出値と前記端子電圧初期値との差分値を
演算する手段と、 前記電流計測値と前記差分値を、電池モデルを用いた適
応デジタルフィルタに入力し、その結果として推定演算
された開路電圧の変化分に端子電圧初期値を加算して開
路電圧の推定値を演算する開路電圧推定手段と、 予め求めた開路電圧と充電率との関係に基づいて前記の
推定した開路電圧を用いて充電率を推定する充電率推定
手段と、 を備えたことを特徴とする二次電池の充電率推定装置。 - 【請求項２】下記（数１）式に示す連続時間系の電池モ
デルを用いて、適応デジタルフィルタ演算を行い、（数
１）式中のオフセット項であるＶ０および過渡項である
Ａ(ｓ）、Ｂ(ｓ）に対応するパラメータを一括推定する
ことを特徴とする請求項１に記載の二次電池の充電率推
定装置。 【数１】 ただし、ｓはラプラス演算子、Ａ(ｓ）、Ｂ(ｓ）は、ｓ
の多項式関数 - 【請求項３】前記（数１）式に示す連続時間系の電池モ
デルを、離散時間系表現に変換した下記（数２）式で示
される自己回帰型モデルを用いて、適応デジタルフィル
タ演算を行い、各パラメータ（ａｋ、ｂｋ、ｃｋ、
ｒｋ）を一括推定し、さらに、下記（数３）式を用いて
開路電圧Ｖ０を推定演算することを特徴とする請求項２
に記載の二次電池の充電率推定装置。 【数２】 【数３】 ただし、（数２）式、（数３）式において、 Ｖｊ−ｋはｋサンプリング周期前の端子電圧 Ｉｊ−ｋはｋサンプリング周期前の電流 ａｋ、ｂｋ、ｃｋは定数 ｒｊ−ｋはオフセット
項の仮変数 Ｎ，Ｐは次数 ｅｊは白色雑音
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001268314A JP3714214B2 (ja) | 2001-09-05 | 2001-09-05 | 二次電池の充電率推定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001268314A JP3714214B2 (ja) | 2001-09-05 | 2001-09-05 | 二次電池の充電率推定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003075517A true JP2003075517A (ja) | 2003-03-12 |
| JP3714214B2 JP3714214B2 (ja) | 2005-11-09 |
Family
ID=19094316
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001268314A Expired - Fee Related JP3714214B2 (ja) | 2001-09-05 | 2001-09-05 | 二次電池の充電率推定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3714214B2 (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7358704B2 (en) | 2003-12-25 | 2008-04-15 | Nissan Motor Co., Ltd. | Critical state estimation system and method for secondary cells |
| US7486079B2 (en) | 2004-06-11 | 2009-02-03 | Nissan Motor Co., Ltd. | Available input-output power estimating device for secondary battery |
| JP2010127729A (ja) * | 2008-11-27 | 2010-06-10 | Calsonic Kansei Corp | バッテリの劣化推定方法及び装置 |
| JP2010239705A (ja) * | 2009-03-30 | 2010-10-21 | Primearth Ev Energy Co Ltd | 二次電池の充電制御方法及び充電器 |
| JP2011530697A (ja) * | 2008-08-08 | 2011-12-22 | エルジー・ケム・リミテッド | バッテリーセルの電圧変化挙動を用いたセルバランス装置及び方法 |
| WO2013125118A1 (ja) * | 2012-02-22 | 2013-08-29 | カルソニックカンセイ株式会社 | パラメータ推定装置 |
| US20150219726A1 (en) * | 2014-02-03 | 2015-08-06 | GM Global Technology Operations LLC | Systems and methods for battery state estimation |
| CN105527581A (zh) * | 2016-01-30 | 2016-04-27 | 清华大学 | 混合型正极材料锂离子电池关键参数及容量衰减机理的辨识方法 |
| CN105659428A (zh) * | 2013-10-21 | 2016-06-08 | 丰田自动车株式会社 | 电池系统 |
| CN112514138A (zh) * | 2018-07-25 | 2021-03-16 | 松下知识产权经营株式会社 | 管理装置以及电源系统 |
-
2001
- 2001-09-05 JP JP2001268314A patent/JP3714214B2/ja not_active Expired - Fee Related
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7358704B2 (en) | 2003-12-25 | 2008-04-15 | Nissan Motor Co., Ltd. | Critical state estimation system and method for secondary cells |
| US7486079B2 (en) | 2004-06-11 | 2009-02-03 | Nissan Motor Co., Ltd. | Available input-output power estimating device for secondary battery |
| JP2011530697A (ja) * | 2008-08-08 | 2011-12-22 | エルジー・ケム・リミテッド | バッテリーセルの電圧変化挙動を用いたセルバランス装置及び方法 |
| JP2010127729A (ja) * | 2008-11-27 | 2010-06-10 | Calsonic Kansei Corp | バッテリの劣化推定方法及び装置 |
| JP2010239705A (ja) * | 2009-03-30 | 2010-10-21 | Primearth Ev Energy Co Ltd | 二次電池の充電制御方法及び充電器 |
| JP5319854B1 (ja) * | 2012-02-22 | 2013-10-16 | カルソニックカンセイ株式会社 | パラメータ推定装置 |
| WO2013125118A1 (ja) * | 2012-02-22 | 2013-08-29 | カルソニックカンセイ株式会社 | パラメータ推定装置 |
| CN105659428A (zh) * | 2013-10-21 | 2016-06-08 | 丰田自动车株式会社 | 电池系统 |
| KR101776546B1 (ko) * | 2013-10-21 | 2017-09-07 | 도요타지도샤가부시키가이샤 | 전지 시스템 |
| US10181622B2 (en) | 2013-10-21 | 2019-01-15 | Toyota Jidosha Kabushiki Kaisha | Cell system |
| US20150219726A1 (en) * | 2014-02-03 | 2015-08-06 | GM Global Technology Operations LLC | Systems and methods for battery state estimation |
| CN105527581A (zh) * | 2016-01-30 | 2016-04-27 | 清华大学 | 混合型正极材料锂离子电池关键参数及容量衰减机理的辨识方法 |
| CN112514138A (zh) * | 2018-07-25 | 2021-03-16 | 松下知识产权经营株式会社 | 管理装置以及电源系统 |
| CN112514138B (zh) * | 2018-07-25 | 2024-04-23 | 松下知识产权经营株式会社 | 管理装置以及电源系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3714214B2 (ja) | 2005-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3714321B2 (ja) | 二次電池の充電率推定装置 | |
| JP3714333B2 (ja) | 二次電池の入出力可能電力推定装置 | |
| JP3747826B2 (ja) | 二次電池の充電率推定装置 | |
| JP5842421B2 (ja) | 電池状態推定装置 | |
| JP2023518778A (ja) | 電池充電状態を決定する方法及び装置、電池管理システム | |
| EP2089731A2 (en) | Apparatus and method for determination of the state-of-charge of a battery when the battery is not in equilibrium | |
| JP2007024703A (ja) | 二次電池の充電率推定装置 | |
| JP2007147487A (ja) | 二次電池の入出力可能電力推定装置 | |
| JP2010127729A (ja) | バッテリの劣化推定方法及び装置 | |
| JP2006098135A (ja) | バッテリの劣化度推定装置 | |
| JP2006292492A (ja) | 二次電池の満充電容量推定装置 | |
| JP2007139536A (ja) | 二次電池の入出力可能電力推定装置 | |
| JP7183576B2 (ja) | 二次電池パラメータ推定装置、二次電池パラメータ推定方法及びプログラム | |
| JP2015224927A (ja) | 蓄電残量推定装置、蓄電池の蓄電残量を推定する方法、及びコンピュータプログラム | |
| JP2003075517A (ja) | 二次電池の充電率推定装置 | |
| JP3714314B2 (ja) | 二次電池の充電率推定装置 | |
| JP3714246B2 (ja) | 二次電池の充電率推定装置 | |
| JP2007033112A (ja) | 二次電池の充電率推定装置 | |
| JP3714330B2 (ja) | 二次電池の充電率推定装置 | |
| CN116699415A (zh) | 一种动态可重构电池系统电量估计方法、系统及电子设备 | |
| JP3852372B2 (ja) | 二次電池の充電率推定装置 | |
| JP3852371B2 (ja) | 二次電池の充電率推定装置 | |
| JP4103569B2 (ja) | 二次電池の充電率推定装置 | |
| JP4720364B2 (ja) | 二次電池の内部抵抗推定装置 | |
| JP4666149B2 (ja) | 二次電池の入出力可能電力推定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050714 |
|
| 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: 20050802 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050815 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |