JP2012159414A - Battery status estimation device - Google Patents

Battery status estimation device Download PDF

Info

Publication number
JP2012159414A
JP2012159414A JP2011019822A JP2011019822A JP2012159414A JP 2012159414 A JP2012159414 A JP 2012159414A JP 2011019822 A JP2011019822 A JP 2011019822A JP 2011019822 A JP2011019822 A JP 2011019822A JP 2012159414 A JP2012159414 A JP 2012159414A
Authority
JP
Japan
Prior art keywords
battery
voltage
secondary battery
state
cutoff frequency
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
Application number
JP2011019822A
Other languages
Japanese (ja)
Other versions
JP5707982B2 (en
Inventor
Shiho Umeki
志保 梅木
Hideo Nakamura
英夫 中村
Kazuhiko Tazoe
和彦 田添
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.)
Nissan Motor Co Ltd
Original Assignee
Nissan Motor Co 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 Nissan Motor Co Ltd filed Critical Nissan Motor Co Ltd
Priority to JP2011019822A priority Critical patent/JP5707982B2/en
Publication of JP2012159414A publication Critical patent/JP2012159414A/en
Application granted granted Critical
Publication of JP5707982B2 publication Critical patent/JP5707982B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10Energy storage using batteries

Landscapes

  • Tests Of Electric Status Of Batteries (AREA)
  • Secondary Cells (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a battery status estimation device capable of enhancing calculation accuracy of the internal state of a secondary battery including a nonlinear region thereof.SOLUTION: A battery status estimation device includes: current detection means for detecting a current of a secondary battery as a current measuring value; voltage detection means for detecting a terminal voltage of the secondary battery as a voltage measurement value; terminal voltage estimation means that defines a battery model of the secondary battery, and calculates a conversion state quantity by a state quantity conversion of the current measuring value and the voltage measurement value by using a state variable filter based on the battery model in order to estimate a terminal voltage of the secondary battery based on the battery model as a voltage estimation value, from the conversion state quantity; identification means for identifying a parameter of the secondary battery so that a difference between the voltage measurement value and the voltage estimation value converges to zero; and nonlinear region detection means for detecting a presence ratio of a nonlinear region in characteristics between current and voltage in the secondary battery. The terminal voltage estimation means sets a cutoff frequency of the state variable filter according to the presence ratio.

Description

本発明は、二次電池内部の状態を推定する電池状態推定装置に関するものである。   The present invention relates to a battery state estimation device for estimating a state inside a secondary battery.

二次電池の制御装置として、所定の電池モデルを定義し、二次電池の電流および端子電圧の計測値を、電池モデルに基づく状態変数フィルタを用いて状態量変換し、これを用いて、電池モデルに基づく二次電池の端子電圧を推定し、電圧計測値と、電池モデルに基づいて推定された端子電圧との差分がゼロに収束するように、二次電池のパラメータを同定する制御装置が知れている(特許文献1参照)。   As a control device for a secondary battery, a predetermined battery model is defined, and the measured values of the secondary battery current and terminal voltage are converted into state quantities using a state variable filter based on the battery model. A control device that estimates the terminal voltage of the secondary battery based on the model and identifies the parameters of the secondary battery so that the difference between the voltage measurement value and the terminal voltage estimated based on the battery model converges to zero. It is known (see Patent Document 1).

特開2003−185719号公報JP 2003-185719 A

しかしながら、当該従来の制御装置は、二次電池の電流−電圧特性の線形領域に基づく入力検出値を用いて、二次電池の内部状態を推定する装置であった。そのため、従来の制御装置により非線形領域を含む二次電池の内部状態を推定する場合には、当該従来の制御装置は、非線形領域の影響を受け、演算精度が低下するという問題があった。   However, the conventional control device is a device that estimates the internal state of the secondary battery using the input detection value based on the linear region of the current-voltage characteristic of the secondary battery. Therefore, when the internal state of the secondary battery including the nonlinear region is estimated by the conventional control device, the conventional control device is affected by the nonlinear region and has a problem that the calculation accuracy is lowered.

本発明が解決しようとする課題は、非線形領域を含む二次電池の内部状態の演算精度を高めることができる電池状態推定装置を提供することである。   The problem to be solved by the present invention is to provide a battery state estimation device that can improve the calculation accuracy of the internal state of a secondary battery including a nonlinear region.

本発明は、二次電池の電流および端子電圧を検出し、検出した電流と電圧の計測値及び状態変数フィルタを用いて、所定の電池モデルに基づく二次電池の端子電圧を推定し、端子電圧の計測値との端子電圧の推定値との差分がゼロに収束するように、二次電池のパラメータを同定する電池状態推定装置において、当該二次電池の非線形領域の存在割合に応じて、当該状態変数フィルタのカットオフ周波数を設定することにより、上記課題を解決する。   The present invention detects a current and a terminal voltage of a secondary battery, estimates a terminal voltage of a secondary battery based on a predetermined battery model using a measured value of the detected current and voltage, and a state variable filter. In the battery state estimation device that identifies the parameters of the secondary battery so that the difference between the measured value of the terminal value and the estimated value of the terminal voltage converges to zero, according to the existence ratio of the non-linear region of the secondary battery, The above problem is solved by setting the cutoff frequency of the state variable filter.

本発明によれば、状態変数フィルタにおいて、非線形領域の信号を減衰させて線形領域になる信号を用いて、変換状態量を演算しパラメータを同定するため、二次電池の内部状態の演算精度を高めることができる。   According to the present invention, in the state variable filter, the conversion state quantity is calculated and the parameter is identified using the signal that becomes the linear region by attenuating the signal in the nonlinear region. Can be increased.

発明の実施形態に係る電池状態推定装置の構成を示す図である。It is a figure which shows the structure of the battery state estimation apparatus which concerns on embodiment of invention. 図1の電子制御ユニットのブロック図である。It is a block diagram of the electronic control unit of FIG. 図1の二次電池の電池モデルを示す等価回路モデルを示す図である。It is a figure which shows the equivalent circuit model which shows the battery model of the secondary battery of FIG. 図1の二次電池の電流−電圧特性を示すグラフである。2 is a graph showing current-voltage characteristics of the secondary battery of FIG. 1. 図1の電池状態推定装置の制御手順を示すフローチャートである。It is a flowchart which shows the control procedure of the battery state estimation apparatus of FIG. 二次電池の電流−電圧特性(電池温度−25度)を示すグラフである。It is a graph which shows the current-voltage characteristic (battery temperature-25 degree) of a secondary battery. 二次電池の電流−電圧特性(電池温度−0度)を示すグラフである。It is a graph which shows the current-voltage characteristic (battery temperature-0 degree | times) of a secondary battery. 二次電池の電流−電圧特性(電池温度−30度)を示すグラフである。It is a graph which shows the current-voltage characteristic (battery temperature-30 degree | times) of a secondary battery. 二次電池の電流−電圧特性(劣化前)を示すグラフである。It is a graph which shows the current-voltage characteristic (before deterioration) of a secondary battery. 二次電池の電流−電圧特性(劣化後)を示すグラフである。It is a graph which shows the current-voltage characteristic (after deterioration) of a secondary battery. 発明の他の実施形態に係る電池状態推定装置の電子制御ユニット30のブロック図である。It is a block diagram of the electronic control unit 30 of the battery state estimation apparatus which concerns on other embodiment of invention. 実施例1における充電率の推定処理のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the estimation process of the charging rate in Example 1. FIG. 実施例1における充電率の推定処理のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the estimation process of the charging rate in Example 1. FIG. 実施例2における充電率の推定処理のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the estimation process of the charging rate in Example 2. FIG. 実施例2における充電率の推定処理のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the estimation process of the charging rate in Example 2. FIG. 実施例2における充電率の推定処理のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the estimation process of the charging rate in Example 2. FIG. 実施例2における充電率の推定処理のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the estimation process of the charging rate in Example 2. FIG. 実施例2における充電率の推定処理のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the estimation process of the charging rate in Example 2. FIG.

以下、本発明の実施形態を図面に基づいて説明する。   Hereinafter, embodiments of the present invention will be described with reference to the drawings.

《第1実施形態》
図1は、本実施形態に係る二次電池の電池状態推定装置の構成を示す図である。図1に示す電池状態推定装置は、二次電池でモータ等の負荷を駆動したり、モータの回生による電力やエンジンを動力源としてオルタネータで発電した電力で二次電池を充電するシステムに、本発明に係る二次電池の制御装置を適用した例である。
<< First Embodiment >>
FIG. 1 is a diagram illustrating a configuration of a battery state estimation device for a secondary battery according to the present embodiment. The battery state estimation apparatus shown in FIG. 1 is used in a system for driving a load such as a motor with a secondary battery, or for charging a secondary battery with electric power generated by motor regeneration or electric power generated by an alternator using an engine as a power source. It is the example which applied the control apparatus of the secondary battery which concerns on invention.

二次電池10は、複数の単位電池を直列に接続してなるものである。二次電池10を構成する単位電池としては、たとえば、リチウムイオン二次電池などのリチウム系二次電池などが挙げられる。負荷20としては、たとえば、モータなどが挙げられる。   The secondary battery 10 is formed by connecting a plurality of unit batteries in series. Examples of the unit battery constituting the secondary battery 10 include a lithium secondary battery such as a lithium ion secondary battery. An example of the load 20 is a motor.

電流センサ40は、二次電池10に流れる充放電電流を検出するセンサであり、電流センサ40により検出された信号は、電子制御ユニット30へ送出される。また、電圧センサ50は、二次電池10の端子電圧を検出するセンサであり、電圧センサ50により検出された信号は、電子制御ユニット30へ送出される。   The current sensor 40 is a sensor that detects a charging / discharging current flowing through the secondary battery 10, and a signal detected by the current sensor 40 is sent to the electronic control unit 30. The voltage sensor 50 is a sensor that detects the terminal voltage of the secondary battery 10, and a signal detected by the voltage sensor 50 is sent to the electronic control unit 30.

温度センサ60は、二次電池10の電池温度を検出するセンサであり、温度センサ60により検出された信号は、電子制御ユニット30へ送出される。また、劣化センサ70は、二次電池10の劣化度を検出するセンサであり、劣化センサ70により検出された信号は、電子制御ユニット30へ送出される。なお、劣化センサではなく、劣化度を推定する推定部でもよい。   The temperature sensor 60 is a sensor that detects the battery temperature of the secondary battery 10, and a signal detected by the temperature sensor 60 is sent to the electronic control unit 30. The deterioration sensor 70 is a sensor that detects the degree of deterioration of the secondary battery 10, and a signal detected by the deterioration sensor 70 is sent to the electronic control unit 30. Note that an estimation unit that estimates the degree of deterioration may be used instead of the deterioration sensor.

電子制御ユニット30は、二次電池10を制御するための制御ユニットであり、プログラムを演算するCPU、プログラムや演算結果を記憶するROMおよびRAMから構成されるマイクロコンピュータと電子回路等で構成される。   The electronic control unit 30 is a control unit for controlling the secondary battery 10 and includes a CPU that calculates a program, a microcomputer that includes a ROM and RAM that stores programs and calculation results, an electronic circuit, and the like. .

図2に、電子制御ユニット30の機能ブロック図を示す。図2に示すように、電子制御ユニット30は、電流検出部301、電圧検出部302、状態変数フィルタ演算部303、適応同定器304、状態観測機305、SOC変換部306、非線形領域検出部308を備える。   FIG. 2 shows a functional block diagram of the electronic control unit 30. As shown in FIG. 2, the electronic control unit 30 includes a current detector 301, a voltage detector 302, a state variable filter calculator 303, an adaptive identifier 304, a state observer 305, an SOC converter 306, and a nonlinear region detector 308. Is provided.

電流検出部301は、電流計40からの信号を所定周期で取得し、電流計40からの信号に基づき、二次電池10に流れる充放電電流を検出することにより、電流計測値I(k)を取得する。電流検出部301は、取得した電流計測値I(k)を状態変数フィルタ演算部303に送出する。   The current detection unit 301 acquires a signal from the ammeter 40 at a predetermined cycle, and detects a charge / discharge current flowing through the secondary battery 10 based on the signal from the ammeter 40, thereby measuring a current measurement value I (k). To get. The current detection unit 301 sends the acquired current measurement value I (k) to the state variable filter calculation unit 303.

電圧検出部302は、電圧計50からの信号を所定周期で取得し、電圧計50からの信号に基づき、二次電池10の端子電圧を検出することにより、電圧計測値V(k)を取得する。電圧検出部302は、取得した電流計測値V(k)を状態変数フィルタ演算部303に送出する。   The voltage detection unit 302 acquires a signal from the voltmeter 50 at a predetermined cycle, and acquires a voltage measurement value V (k) by detecting a terminal voltage of the secondary battery 10 based on the signal from the voltmeter 50. To do. The voltage detection unit 302 sends the acquired current measurement value V (k) to the state variable filter calculation unit 303.

状態変数フィルタ演算部303は、二次電池10の電池モデルを定義し、電流検出部301により検出された電流計測値I(k)および電圧検出部302により検出された電圧計測値V(k)から、適応デジタルフィルタ演算により、状態量変換して変換状態量ω(k)を演算する。   The state variable filter calculation unit 303 defines a battery model of the secondary battery 10, and the current measurement value I (k) detected by the current detection unit 301 and the voltage measurement value V (k) detected by the voltage detection unit 302. Then, the state quantity is converted by adaptive digital filter calculation to calculate the converted state quantity ω (k).

また適応同定器304は、状態変数フィルタ演算部303の演算結果から、二次電池10の電池モデルの電池パラメータφ^(k)を一括推定し、必要に応じて、当該電池パラメータφ^(k)を変換する。   The adaptive identifier 304 collectively estimates the battery parameter φ ^ (k) of the battery model of the secondary battery 10 from the calculation result of the state variable filter calculation unit 303, and if necessary, the battery parameter φ ^ (k ).

ここで、φ^(k)における右肩に付した「^」は、その値が推定値であることを示す。また、図2中では、推定値である「^」を、それぞれ、φ(k)の「φ」の真上、V(k)の「V」の真上、SOC(k)の「S」の真上としているが、下記式(1)に示すように、これはφ^(k)、V^(k)、SOC^(k)と同義である。以下、V^(k)においても同様である。 Here, “^” attached to the right shoulder in φ ^ (k) indicates that the value is an estimated value. In FIG. 2, “^”, which is an estimated value, is set directly above “φ” of φ (k), directly above “V” of V 0 (k), and “S” of SOC (k). However, as shown in the following equation (1), this is synonymous with φ ^ (k), V 0 ^ (k), and SOC ^ (k). The same applies to V ^ (k) below.

Figure 2012159414
以下、状態変数フィルタ演算部303による変換状態量ω(k)の演算方法及び、適応同定器304による電池パラメータφ^(k)の推定方法について説明する。
Figure 2012159414
Hereinafter, the calculation method of the conversion state quantity ω (k) by the state variable filter calculation unit 303 and the estimation method of the battery parameter φ ^ (k) by the adaptive identifier 304 will be described.

まず、本実施形態で用いる「電池モデル」について、説明する。図3は、二次電池10の電池モデルを示す等価回路モデルである。そして、二次電池10の電池モデルを式(1)のように定義する。   First, the “battery model” used in the present embodiment will be described. FIG. 3 is an equivalent circuit model showing a battery model of the secondary battery 10. And the battery model of the secondary battery 10 is defined like Formula (1).

Figure 2012159414
ここで、A(s)、B(s)の次数を1次とした例で状態変数フィルタの導入を簡単に説明する。図3に示す等価回路モデルは、式(2)の形で表すと、下記の式(3)で表される。
Figure 2012159414
Here, the introduction of the state variable filter will be briefly described with an example in which the orders of A (s) and B (s) are the first order. The equivalent circuit model shown in FIG. 3 is expressed by the following expression (3) when expressed in the form of expression (2).

Figure 2012159414
ここで、モデル入力は電流I[A](正値は充電、負値は放電)、モデル出力は端子電圧V[V]であり、R〔Ω]は電荷移動抵抗、R[Ω]は純抵抗、C[F]は電気二重層容量、V[V]は開路電圧である。また、上記式(2)中、sは微分オペレータである。なお、本実施形態に係る電池モデルは、正極、負極を特に分離していないリダクションモデル(1次)であるが、実際の電池の充放電特性を比較的正確に示すことが可能である。このように本実施形態においては、電池モデルの次数を1次にした構成を例として説明する。
Figure 2012159414
Here, the model input is current I [A] (positive value is charging, negative value is discharging), model output is terminal voltage V [V], R 1 [Ω] is charge transfer resistance, R 2 [Ω] Is a pure resistance, C 1 [F] is an electric double layer capacitance, and V 0 [V] is an open circuit voltage. In the above formula (2), s is a differential operator. Note that the battery model according to the present embodiment is a reduction model (primary) in which the positive electrode and the negative electrode are not particularly separated, but it is possible to show the actual charge / discharge characteristics of the battery relatively accurately. As described above, in the present embodiment, a configuration in which the order of the battery model is first will be described as an example.

そして、R、R、Cを下記式(4)のように表すと、上記式(3)は、下記式(5)で表されることとなる。 When the representative of R 1, R 2, C 1 as shown in Equation (4), the equation (3), and thus represented by the following formula (5).

Figure 2012159414
Figure 2012159414

Figure 2012159414
そして、本実施形態においては、上記式(5)に示される電池モデルから、状態変数フィルタ演算部303によって、変換状態量ω(k)を演算し、適応同定器304は電池パラメータφ^(k)の推定を行なう。
Figure 2012159414
In this embodiment, the state variable filter calculation unit 303 calculates the conversion state quantity ω (k) from the battery model shown in the above formula (5), and the adaptive identifier 304 uses the battery parameter φ ^ (k ) Is estimated.

まず、開路電圧V(t)は、電流I(t)に可変なパラメータhを乗じたものをある初期状態から積分したものと考えれば、開路電圧V(t)は、下記式(6)で表すことができる。 First, assuming that the open circuit voltage V 0 (t) is obtained by integrating a current I (t) multiplied by a variable parameter h from an initial state, the open circuit voltage V 0 (t) is expressed by the following formula (6 ).

Figure 2012159414
そして、上記式(5)に、上記式(6)を代入すると、下記式(7)となり、これを整理すると下記式(8)となる。
Figure 2012159414
When the above formula (6) is substituted into the above formula (5), the following formula (7) is obtained, and when this is arranged, the following formula (8) is obtained.

Figure 2012159414
Figure 2012159414

Figure 2012159414
そして、式(8)より、状態空間方程式を適応同定器304の設計に便利な非最小実現へ変換する。
Figure 2012159414
Then, from Equation (8), the state space equation is converted into a non-minimum realization that is convenient for the design of the adaptive identifier 304.

ここで、本例の状態変数フィルタ演算部303は、二次電池10の電流−電圧特性の非線形領域の存在割合に応じて、上記の状態変数フィルタのカットオフ周波数を設定する。状態変数フィルタ演算部303には、予めカットオフ周波数(λa)とカットオフ周波数(λb)が設定されている。λaはλbより高い周波数である。なお、非線形領域の存在割合は、後述する非線形領域存在割合検出部308に基づき特定される。そして、本例の状態変数フィルタ演算部303は、二次電池10の非線形領域の存在割合に応じて、カットオフ周波数(λa)の状態変数フィルタとカットオフ周波数(λb)の状態変数フィルタとを使い分けて、下記の演算を行う。   Here, the state variable filter calculation unit 303 of this example sets the cutoff frequency of the state variable filter according to the existence ratio of the non-linear region of the current-voltage characteristic of the secondary battery 10. In the state variable filter calculation unit 303, a cutoff frequency (λa) and a cutoff frequency (λb) are set in advance. λa is a higher frequency than λb. The existence ratio of the non-linear region is specified based on a non-linear region existence ratio detection unit 308 described later. Then, the state variable filter calculation unit 303 of this example performs a state variable filter with a cutoff frequency (λa) and a state variable filter with a cutoff frequency (λb) according to the existence ratio of the nonlinear region of the secondary battery 10. The following calculation is performed properly.

状態変数フィルタ演算部303は、下記の式(9)と式(10)、及び、式(11)と式(12)により表される、それぞれの状態変数フィルタを有している。式(9)〜式(12)は、式(8)に、既知定数k(i=1,2、・・・,n)を導入することにより、得られる。 The state variable filter calculation unit 303 includes state variable filters represented by the following expressions (9) and (10), and expressions (11) and (12). Expressions (9) to (12) are obtained by introducing known constants k i (i = 1, 2,..., N) into Expression (8).

Figure 2012159414
Figure 2012159414

Figure 2012159414
Figure 2012159414

Figure 2012159414
Figure 2012159414

Figure 2012159414
式(9)及び式(10)は、カットオフ周波数がλaである場合の状態変数フィルタを示し、式(11)及び式(12)は、カットオフ周波数がλbである場合の状態変数フィルタを示す。
Figure 2012159414
Equations (9) and (10) show the state variable filter when the cutoff frequency is λa, and Equations (11) and (12) show the state variable filter when the cutoff frequency is λb. Show.

なお、上記式(10)及び式(12)において、Ii_1、b0i_1、Ii_2、b0i_2は、未知パラメータ(T,T,K,h)を含むパラメータであり、fVi_1、fIi_1、fVi_2、fIi_2は、電流計40および電圧計50により計測可能な値であるI(k)、V(k)を状態変数フィルタによりフィルタ処理を施した変換状態量である。そして、上記式(10)及び式(12)は、これらの積和式になっているため、適応デジタルフィルタの標準形である下記式(13)と一致する。 In the above formulas (10) and (12), I i_1 , b 0i_1 , I i_2 , b 0i_2 are parameters including unknown parameters (T 1 , T 2 , K, h), and f Vi — 1 , f Ii_1 , f Vi_2 , and f Ii_2 are conversion state quantities obtained by filtering I (k) and V (k), which are values measurable by the ammeter 40 and the voltmeter 50, with a state variable filter. And since the said Formula (10) and Formula (12) are these product-sum formulas, it corresponds with following formula (13) which is a standard form of an adaptive digital filter.

Figure 2012159414
上記式(13)中、φ =[Ii_1,b0i_1]、ω=[fVi_1,fIi_1] であり、φ =[Ii_2,b0i_2]、ω=[fVi_2,fIi_2] である。
Figure 2012159414
In the above formula (13), φ 1 T = [I i_1, b 0i_1], ω 1 = [f Vi_1, f Ii_1] a, φ 2 T = [I i_2 , b 0i_2], ω 2 = [f Vi_2 , F Ii — 2 ].

そして、変換状態量であるω(k)又はω(k)から、上述した電池モデルから推定される二次電池10の端子電圧の推定値である電圧推定値V^(k)と、電圧計50で検出され電圧検出部302により取得された実際の計測値である電圧計測値V(k)と、の差分がゼロに収束するように、適応調整則により、下記式(14)に示すアルゴリズムに基づいて、適応同定器304は、電池モデルの電池パラメータφ^(k)の同定を行なう。なお、この際において、本実施形態では、単純な「最小二乗法による適応デジタルフィルタ」の論理的な欠点(一度推定値が収束すると、その後パラメータが変化しても再度正確な推定ができないこと)を改善した「両限トレースゲイン方式」を用いることができる。 Then, from the conversion state quantity ω 1 (k) or ω 2 (k), a voltage estimated value V ^ (k) that is an estimated value of the terminal voltage of the secondary battery 10 estimated from the battery model described above, The following equation (14) is obtained by the adaptive adjustment law so that the difference between the actual measured value detected by the voltmeter 50 and the voltage measured value V (k) obtained by the voltage detector 302 converges to zero. Based on the algorithm shown, the adaptive identifier 304 identifies the battery parameter φ ^ (k) of the battery model. At this time, in this embodiment, the logical disadvantage of a simple “adaptive digital filter based on the least square method” (once the estimated value converges, accurate estimation cannot be performed again even if the parameter changes thereafter). It is possible to use a “both limit trace gain method” that improves the above.

Figure 2012159414
上記式(14)は、電池パラメータφ^(k)を適応的に求める逐次式であり、γ(k)、Γ(k−1)は、共に適応ゲインであり、これらのうち、γ(k)はスカラゲイン(誤差ゲイン)であり、Γ(k−1)は行列ゲイン(信号ゲイン)である。そして、上記式(13)により、k時点における状態量ζ(k)が得られた際には、電池モデルから推定される二次電池10の端子電圧の推定値である電圧推定値V^(k)と、電圧計50で検出され電圧検出部302により取得された電圧計測値V(k)との差分であるe(k)を求めることができ、このe(k)をゼロに収束させることにより、電池パラメータφ^(k)を逐次的に算出することができる。ここで、λaの状態変数フィルタにより算出された電池パラメータを、φ^(k)、λbの状態変数フィルタにより算出された電池パラメータを、φ^(k)とする。
Figure 2012159414
The above equation (14) is a sequential equation for adaptively obtaining the battery parameter φ ^ (k), and γ (k) and Γ (k−1) are both adaptive gains, and among these, γ (k ) Is a scalar gain (error gain), and Γ (k−1) is a matrix gain (signal gain). Then, when the state quantity ζ (k) at time point k is obtained by the above equation (13), the estimated voltage value V ^ () that is the estimated value of the terminal voltage of the secondary battery 10 estimated from the battery model. k) and a difference e (k) between the voltage measurement value V (k) detected by the voltmeter 50 and acquired by the voltage detector 302 can be obtained, and this e (k) is converged to zero. Accordingly, the battery parameter φ ^ (k) can be calculated sequentially. Here, the battery parameter calculated by the state variable filter of λa is φ 1 ^ (k), and the battery parameter calculated by the state variable filter of λb is φ 2 ^ (k).

電池パラメータφ^(k)、φ^(k)は、上述したように、未知パラメータ(T,T,K,h)を含むパラメータIi_1,b0i_1、Ii_2,b0i_2にそれぞれ相当する。そのため、状態観測器305は、状態変数フィルタ演算部303により算出された電池パラメータφ^(k)および変換状態量ω(k)を用い、これらを、上記式(4)に代入することで、開路電圧推定値V^(k)を求めることができる。このようにして得られた開路電圧推定値V^(k)は、状態観測器305により、SOC変換部306に送出される。 Cell parameters φ 1 ^ (k), φ 2 ^ (k) , as described above, the parameter I i_1 containing unknown parameters (T 1, T 2, K , h), b 0i_1, I i_2, the b 0I_2 Each corresponds. Therefore, the state observer 305 uses the battery parameter φ ^ (k) and the converted state quantity ω (k) calculated by the state variable filter calculation unit 303, and substitutes them into the above equation (4). An open circuit voltage estimated value V 0 ^ (k) can be obtained. The open circuit voltage estimated value V 0 ^ (k) obtained in this way is sent to the SOC converter 306 by the state observer 305.

SOC変換部306は、開路電圧推定値V^(k)に基づき、二次電池10の充電状態であるSOCに変換し演算する。開路電圧とSOCには所定の関係があるため、例えば、SOC変換部306は、予め開路電圧とSOCとの関係を示すテーブルを格納し、当該テーブルを参照して、SOCを演算する。これにより、二次電池10の内部状態が推定される。 The SOC conversion unit 306 converts and calculates the SOC that is the state of charge of the secondary battery 10 based on the open circuit voltage estimated value V 0 ^ (k). Since there is a predetermined relationship between the open circuit voltage and the SOC, for example, the SOC conversion unit 306 stores a table indicating the relationship between the open circuit voltage and the SOC in advance, and calculates the SOC with reference to the table. Thereby, the internal state of the secondary battery 10 is estimated.

次に、非線形領域存在割合検出部308について、図4を用いて説明する。図4は、二次電池10の電流−電圧特性を示すグラフである。   Next, the non-linear region existence ratio detection unit 308 will be described with reference to FIG. FIG. 4 is a graph showing the current-voltage characteristics of the secondary battery 10.

非線形領域存在割合検出部308は、例えば、本例の推定装置の起動時であって、二次電池10が無負荷な状態の時に、二次電池10の開放電圧VOCVを演算する。または、非線形領域存在割合検出部308は、電流検出部301の計測値を積算することで算出されるSOCからVOCVを演算する。なお、計測値を積算する際には、積算誤差が蓄積されないように定期的に初期化を行ってもよい。そして、非線形領域存在割合検出部308は、電流検出部301及び電圧検出部302の所定周期のサンプリングデータを用いて、電流−電圧の特性上に、サンプリングデータ及び開放電圧VOCVをプロットする。そして、プロット密度の高いエリアを抽出する。抽出されたエリア、当該電流−電圧の特性上において、予め決められた所定の領域内に収まるか否かによって、非線形領域の存在割合が判定される。 For example, the non-linear region existence ratio detection unit 308 calculates the open-circuit voltage V OCV of the secondary battery 10 when the estimation device of this example is activated and the secondary battery 10 is in a no-load state. Alternatively, the non-linear region presence ratio detection unit 308 calculates V OCV from the SOC calculated by integrating the measurement values of the current detection unit 301. When integrating the measurement values, initialization may be performed periodically so that integration errors are not accumulated. Then, the non-linear region existence ratio detection unit 308 plots the sampling data and the open-circuit voltage V OCV on the current-voltage characteristics using the sampling data of the predetermined period of the current detection unit 301 and the voltage detection unit 302. Then, an area having a high plot density is extracted. The existence ratio of the non-linear region is determined depending on whether or not the extracted area and the current-voltage characteristics fall within a predetermined region.

ここで、予め決められた所定の領域とは、図4の斜線部で示すように、演算された開放電圧VOCVを切片とするIV基準直線(P)を中心とした、所定の幅で特定される領域である。なお、IV基準直線(P)の傾きは、例えば、二次電池10の初期の状態における内部抵抗が用いられる。 Here, the predetermined area is specified by a predetermined width centered on the IV reference straight line (P) with the calculated open-circuit voltage V OCV as an intercept, as indicated by the hatched portion in FIG. It is an area to be done. For example, the internal resistance in the initial state of the secondary battery 10 is used for the inclination of the IV reference straight line (P).

例えば、二次電池10が低中領域の電流域で使用されている場合には、電流−電圧の特性は線形性が保たれるため、サンプリングデータは、当該斜線部の中に収まる(図4のエリアa)。一方、二次電池10が低中領域の電流域で使用されている場合(図4のエリアb)には、電流−電圧の特性は非線形性をもつようになり、サンプリングデータは、当該斜線部から外れる。   For example, when the secondary battery 10 is used in a low-middle current region, the current-voltage characteristics are kept linear, so that the sampling data is within the shaded portion (FIG. 4). Area a). On the other hand, when the secondary battery 10 is used in the current region in the low to middle region (area b in FIG. 4), the current-voltage characteristic has nonlinearity, and the sampling data is the shaded portion. Deviate from.

非線形領域存在割合検出部308は、サンプリングデータで抽出されるエリアが斜線部内にある時(図4のエリアa)には、二次電池10は線形性が保たれており、非線形領域の存在割合は所定の割合より低いと判断する。一方、サンプリングデータで抽出されるエリアが斜線部外にある時(図4のエリアb)には、二次電池10は線形性が保たれておらず、非線形領域の存在割合が所定の割合より高いと判断する。そして、非線形領域存在割合検出部308は、当該判断結果に基づく信号を状態変数フィルタ演算部303に送信し、状態変数フィルタ演算部303は、当該存在割合に基づいて、存在割合が小さい場合には、カットオフ周波数をλaに、存在割合が大きい場合には、カットオフ周波数をλbとする。   When the area extracted by the sampling data is within the shaded area (area a in FIG. 4), the non-linear area existence ratio detection unit 308 maintains the linearity, and the non-linear area existence ratio. Is determined to be lower than a predetermined ratio. On the other hand, when the area extracted by the sampling data is outside the shaded area (area b in FIG. 4), the secondary battery 10 does not maintain linearity, and the existence ratio of the non-linear region is higher than the predetermined ratio. Judged to be high. Then, the nonlinear region existence ratio detection unit 308 transmits a signal based on the determination result to the state variable filter calculation unit 303, and the state variable filter calculation unit 303 determines that the existence ratio is small based on the existence ratio. When the cutoff frequency is λa and the existence ratio is large, the cutoff frequency is λb.

次に、本例の電池状態推定装置の制御手順を、図5を用いて説明する。図5は、本例の電池状態推定装置の制御手順を示すフローチャートである。   Next, the control procedure of the battery state estimation apparatus of this example will be described with reference to FIG. FIG. 5 is a flowchart showing a control procedure of the battery state estimation device of this example.

まず、ステップS1にて、電流検出部301、および電圧検出部302により、電流計測値I(k)、および電圧計測値V(k)の取得が行われ、状態変数フィルタ演算部303に送出される。   First, in step S <b> 1, the current measurement value I (k) and the voltage measurement value V (k) are acquired by the current detection unit 301 and the voltage detection unit 302 and sent to the state variable filter calculation unit 303. The

ステップS2にて、非線形領域存在割合検出部308は、開路電圧VOCVと、ステップS1の計測値に基づいて、バッテリー12の非線形領域の存在割合を検出する。次に、ステップS3にて、状態変数フィルタ演算部303は、ステップS2の存在割合に応じて、状態変数フィルタのカットオフ周波数をλa又はλbに設定する。 In step S2, the non-linear region existence ratio detection unit 308 detects the non-linear region existence ratio of the battery 12 based on the open circuit voltage V OCV and the measured value in step S1. Next, in step S3, the state variable filter calculation unit 303 sets the cutoff frequency of the state variable filter to λa or λb according to the existence ratio in step S2.

ステップS4にて、電流計測値I(k)および電圧計測値V(k)について、状態変数フィルタ演算部303により、上記式(9)及び(10)、又は、(11)及び(12)に従い、設定されたカットオフ周波数に基づく状態変数フィルタを用いたフィルタ処理が行なわれ、変換状態量ω(k)が算出される。   In step S4, the current variable value I (k) and the voltage measurement value V (k) are calculated by the state variable filter calculation unit 303 according to the above formulas (9) and (10) or (11) and (12). Then, a filter process using a state variable filter based on the set cutoff frequency is performed, and a conversion state quantity ω (k) is calculated.

ステップS5にて、ステップS4において算出された変換状態量ω(k)及び電圧計測値V(k)を用いて、適応同定器304により、上記式(14)に従い、電池モデルの電池パラメータφ^(k)の同定が行なわれる。   In step S5, using the conversion state quantity ω (k) and the voltage measurement value V (k) calculated in step S4, the adaptive identifier 304 performs battery parameter φ ^ of the battery model according to the above equation (14). Identification of (k) is performed.

ステップS6にて、状態観測器305により、状態変数フィルタ演算部303及び適応同定器304により算出された電池パラメータφ^(k)および変換状態量ω(k)に基づいて、開路電圧推定値V^(k)の算出が行なわれる。 In step S6, based on the battery parameter φ ^ (k) and the converted state quantity ω (k) calculated by the state observer 305 by the state variable filter calculation unit 303 and the adaptive identifier 304, the open circuit voltage estimated value V Calculation of 0 ^ (k) is performed.

そしてステップS7にて、SOC変換部306は、ステップS5の開路電圧推定値V^(k)を用いて、予め定められた二次電池10の開路電圧−SOC特性に基づいて、充電率推定値SOC^(k)の算出が行なわれる。 In step S7, SOC conversion unit 306 estimates the charging rate based on the predetermined open circuit voltage-SOC characteristic of secondary battery 10 using open circuit voltage estimated value V 0 ^ (k) in step S5. The value SOC ^ (k) is calculated.

本実施形態では、以上のようにして二次電池10の電池モデルの電池パラメータφ^(k)および充電率推定値SOC^(k)の推定が行われる。   In the present embodiment, the battery parameter φ ^ (k) and the estimated charge rate SOC ^ (k) of the battery model of the secondary battery 10 are estimated as described above.

上記のように本例は、電流計測値I(k)および電圧計測値V(k)を用いて、二次電池10の電池モデルに基づく状態変数フィルタを用いて変換状態量を演算し、電圧推定値V^(k)を算出し、電圧計測値V(k)および電圧推定値V^(k)の差分e(k)がゼロに収束するように、電池モデルの電池パラメータφ^(k)を推定する際に、二次電池10の非線形領域の存在割合に応じて、状態変種フィルタのカットオフ周波数を設定する。これにより、非線形領域にある電流計測値I(k)および電圧計測値V(k)を減衰させて、線形領域にある計測値を効率的に用いて、変換状態量を演算し、パラメータφ^(k)を推定することができる。   As described above, this example calculates the conversion state quantity using the state variable filter based on the battery model of the secondary battery 10 using the current measurement value I (k) and the voltage measurement value V (k), and the voltage The estimated value V ^ (k) is calculated, and the battery parameter φ ^ (k) of the battery model is set such that the difference e (k) between the measured voltage value V (k) and the estimated voltage value V ^ (k) converges to zero. ) Is set, the cutoff frequency of the state variant filter is set according to the existence ratio of the non-linear region of the secondary battery 10. As a result, the current measurement value I (k) and the voltage measurement value V (k) in the non-linear region are attenuated, and the conversion state quantity is calculated using the measurement value in the linear region efficiently, and the parameter φ ^ (K) can be estimated.

また本発明は、電流計測値I(k)および電圧計測値V(k)に基づき、非線形領域の存在割合を検出し、存在割合が所定の存在割合より大きい場合には、状態変種フィルタのカットオフ周波数を低いカットオフ周波数(λb)にする。これにより、二次電池10の電流−電圧特性は、環境や使用条件により経時的に変化する。そのため、二次電池10の現在の状態を、電流計測値I(k)および電圧計測値V(k)から把握し、二次電池10の非線形領域の存在割合を検出することにより、二次電池10の非線形領域が多くなった場合には、カットオフ周波数を下げて、同定性能を向上させる。これにより演算精度を高めることができる。   The present invention also detects the existence ratio of the nonlinear region based on the measured current value I (k) and the measured voltage value V (k), and if the existence ratio is larger than the predetermined existence ratio, the state variant filter is cut. The off frequency is set to a low cutoff frequency (λb). Thereby, the current-voltage characteristics of the secondary battery 10 change over time depending on the environment and use conditions. Therefore, the secondary battery 10 is grasped from the current measurement value I (k) and the voltage measurement value V (k), and the secondary battery 10 is detected by detecting the existence ratio of the non-linear region of the secondary battery 10. When 10 nonlinear regions increase, the cutoff frequency is lowered to improve the identification performance. Thereby, calculation accuracy can be improved.

なお、本例は、二次電池10の電池温度に基づいて、二次電池10の非線形領域の存在割合を検出し、カットオフ周波数を設定してもよい。図6a〜cは、各温度における電流に対する端子間電圧の特性(IV特性)の一例を示し、図6aは電池温度25度の場合のIV特性を、図6bは電池温度0度の場合のIV特性を、図6cは電池温度−30度の場合のIV特性を示す。図6a〜cに示すように、IV特性は、電池温度を相関性があり、電池温度が低くなるにつれて、非線形領域の存在割合が大きくなる。そして、図6a〜cに示す特性は、二次電池10に使用される電池の特性により予め決まる。そのため、本例において、非線形領域存在割合検出部308は、電池温度と非線形領域の存在割合との相関表を予め記録し、温度センサ60の検出温度に基づいて、存在割合を特定する。そして、状態変数フィルタ演算部303は、非線形領域存在割合検出部308により特定された、非線形領域の存在割合に応じて、二次電池10の非線形領域が多くなった場合には、カットオフ周波数を下げて、同定性能を向上させる。これにより演算精度を高めることができる。   In the present example, the cut-off frequency may be set by detecting the existence ratio of the non-linear region of the secondary battery 10 based on the battery temperature of the secondary battery 10. 6A to 6C show an example of the characteristics (IV characteristics) of the voltage between terminals with respect to the current at each temperature, FIG. 6A shows the IV characteristics when the battery temperature is 25 degrees, and FIG. 6B shows the IV characteristics when the battery temperature is 0 degrees. FIG. 6 c shows the IV characteristics when the battery temperature is −30 degrees. As shown in FIGS. 6a to 6c, the IV characteristics correlate with the battery temperature, and the non-linear region existence ratio increases as the battery temperature decreases. 6a to 6c are determined in advance by the characteristics of the battery used for the secondary battery 10. Therefore, in this example, the non-linear region existence ratio detection unit 308 previously records a correlation table between the battery temperature and the non-linear region existence ratio, and specifies the existence ratio based on the temperature detected by the temperature sensor 60. Then, the state variable filter calculation unit 303 determines the cutoff frequency when the non-linear area of the secondary battery 10 increases according to the non-linear area existence ratio specified by the non-linear area existence ratio detection unit 308. To improve identification performance. Thereby, calculation accuracy can be improved.

なお、本例は、二次電池10の劣化度に基づいて、二次電池10の非線形領域の存在割合を検出し、カットオフ周波数を設定してもよい。図7a及びbは、劣化における電流に対する端子間電圧の特性(IV特性)の一例を示し、図7aは初期の電池のIV特性を、図7bは劣化後の電池のIV特性を示す。図7a及び図7bに示すように、IV特性は、劣化度と相関性があり、劣化するにつれて、非線形領域の存在割合が大きくなる。そして、図7a及び図7bに示す特性は、二次電池10に使用される電池の特性により予め決まる。そのため、本例において、非線形領域存在割合検出部308は、劣化度と非線形領域の存在割合との相関表を予め記録し、劣化度センサ70の検出温度に基づいて、存在割合を特定する。そして、状態変数フィルタ演算部303は、非線形領域存在割合検出部308により特定された、非線形領域の存在割合に応じて、二次電池10の非線形領域が多くなった場合には、カットオフ周波数を下げて、同定性能を向上させる。これにより演算精度を高めることができる。   In the present example, the cut-off frequency may be set by detecting the existence ratio of the non-linear region of the secondary battery 10 based on the degree of deterioration of the secondary battery 10. 7A and 7B show an example of the characteristics (IV characteristics) of the voltage between the terminals with respect to the current during deterioration, FIG. 7A shows the IV characteristics of the initial battery, and FIG. 7B shows the IV characteristics of the battery after deterioration. As shown in FIGS. 7a and 7b, the IV characteristic has a correlation with the degree of deterioration, and the existence ratio of the non-linear region increases as it deteriorates. The characteristics shown in FIGS. 7 a and 7 b are determined in advance by the characteristics of the battery used for the secondary battery 10. Therefore, in this example, the non-linear region existence ratio detection unit 308 records a correlation table between the deterioration degree and the non-linear area existence ratio in advance, and specifies the existence ratio based on the detected temperature of the deterioration degree sensor 70. Then, the state variable filter calculation unit 303 determines the cutoff frequency when the non-linear area of the secondary battery 10 increases according to the non-linear area existence ratio specified by the non-linear area existence ratio detection unit 308. To improve identification performance. Thereby, calculation accuracy can be improved.

なお本例は、電流及び電圧の計測値、電池温度又は電池の劣化度に応じて、非線形領域の存在割合が大きくなった場合に、現在の電流計測値が所定の電流値より小さい場合には、カットオフ周波数を高いカットオフ周波数(λa)に設定してもよい。二次電池10の非線形領域の存在割合が大きくなった場合でも、電流計測値の絶対値の低い領域(小中電流領域)では、二次電池10の線形性が保たれている。一方、絶対値の高い領域(高電流領域)では、二次電池10の線形性が保たれていない。そのため、状態変数フィルタ演算部303は、非線形領域の存在割合が所定の割合よりも大きいと判断された場合であり、電流検出部301の測定値に基づき、小中電流領域が使用されている場合には、カットオフ周波数を高いカットオフ周波数(λa)に設定する。これにより、本例は、演算精度を高めることができる。   In this example, the current current measurement value is smaller than the predetermined current value when the non-linear region exists in proportion to the current and voltage measurement values, battery temperature, or battery deterioration level. The cutoff frequency may be set to a high cutoff frequency (λa). Even when the existence ratio of the non-linear region of the secondary battery 10 increases, the linearity of the secondary battery 10 is maintained in the region where the absolute value of the current measurement value is low (small and medium current region). On the other hand, in the region having a high absolute value (high current region), the linearity of the secondary battery 10 is not maintained. Therefore, the state variable filter calculation unit 303 is a case where the existence ratio of the non-linear region is determined to be larger than a predetermined rate, and the small / medium current region is used based on the measurement value of the current detection unit 301. In this case, the cutoff frequency is set to a high cutoff frequency (λa). Thereby, this example can improve calculation accuracy.

なお、電流の大小が頻繁に切り替わり、カットオフ周波数を毎回切り換えると、パラメータの収束性が悪くなる場合がある。そのため、例えば、ある時間内の平均電流に基づいて、カットオフ周波数を変更してもよい。   Note that if the magnitude of the current is frequently switched and the cut-off frequency is switched each time, the parameter convergence may be deteriorated. Therefore, for example, the cutoff frequency may be changed based on an average current within a certain time.

なお、本例の電流センサ40は本発明の「電流検出手段」に相当し、電圧センサ50が「電圧検出手段」に、状態変数フィルタ演算部303が「端子電圧推定手段」に、適応同定器304が「同定手段」に、状態観測器305が「開路電圧演算手段」に、非線形領域存在割合検出部308が「非線形領域存在割合検出手段」に相当する。   The current sensor 40 of this example corresponds to the “current detection means” of the present invention, the voltage sensor 50 as the “voltage detection means”, the state variable filter calculation unit 303 as the “terminal voltage estimation means”, and the adaptive identifier. Reference numeral 304 corresponds to “identification means”, the state observer 305 corresponds to “open circuit voltage calculation means”, and the non-linear area existence ratio detection unit 308 corresponds to “non-linear area existence ratio detection means”.

《第2実施形態》
図8は、発明の他の実施形態に係る電池状態推定装置の電子制御ユニット30のブロック図である。本例では上述した第1実施形態に対して、電子制御ユニット30の一部が異なる。これ以外の構成は上述した第1実施形態と同じであるため、その記載を援用する。
<< Second Embodiment >>
FIG. 8 is a block diagram of the electronic control unit 30 of the battery state estimation device according to another embodiment of the invention. In this example, a part of the electronic control unit 30 is different from the first embodiment described above. Since the other configuration is the same as that of the first embodiment described above, the description thereof is incorporated.

図8に示すように、電子制御ユニット30はパラメータ変換部307を有し、状態変数フィルタ演算部303は、第1演算部3031と第2演算部3032を含む。状態変数フィルタ演算部303は、複数の状態変数フィルタを用いて、変換状態量ω(k)を演算する。第1演算部3031は、式(9)及び式(10)により表される、カットオフ周波数(λa)の状態変数フィルタ(第1の状態変数フィルタ)を用いて演算し、第2演算部3032は、式(11)及び式(12)により表される、カットオフ周波数(λb)の状態変数フィルタ(第2の状態変数フィルタ)を用いて演算する。   As shown in FIG. 8, the electronic control unit 30 includes a parameter conversion unit 307, and the state variable filter calculation unit 303 includes a first calculation unit 3031 and a second calculation unit 3032. The state variable filter calculation unit 303 calculates the converted state quantity ω (k) using a plurality of state variable filters. The first calculation unit 3031 performs calculation using a state variable filter (first state variable filter) having a cutoff frequency (λa) represented by Expression (9) and Expression (10), and the second calculation unit 3032. Is calculated using a state variable filter (second state variable filter) having a cut-off frequency (λb) represented by Expression (11) and Expression (12).

第1演算部3031は、電池モデルに基づく第1の状態変数フィルタを用いて、変換状態量ω(k)を演算し、第2演算部3032は、電池モデルに基づく第2の状態変数フィルタを用いて、変換状態量ω(k)を演算する。そして、適応同定器304は、変換状態量であるω(k)から、電圧推定値V^(k)と、電圧計50で検出され電圧検出部302により取得された実際の計測値である電圧計測値V(k)と、の差分がゼロに収束するように、適応調整則により、下記式(14)に示すアルゴリズムに基づいて、適応同定器304は、電池モデルの電池パラメータφ^(k)の同定を行なう。 The first calculation unit 3031 calculates the converted state quantity ω 1 (k) using the first state variable filter based on the battery model, and the second calculation unit 3032 uses the second state variable filter based on the battery model. Is used to calculate the conversion state quantity ω 2 (k). The adaptive identifier 304 is a voltage estimation value V ^ (k) from the conversion state quantity ω 1 (k) and an actual measurement value detected by the voltmeter 50 and acquired by the voltage detection unit 302. Based on the algorithm shown in the following equation (14), the adaptive identifier 304 uses the battery parameter φ 1 ^ of the battery model by the adaptive adjustment rule so that the difference from the voltage measurement value V (k) converges to zero. (K) is identified.

パラメータ変換部307は、予め設定されているテーブルまたは数式を参照し、λa及びλbを用いて、適応同定器304により同定されたパラメータφ^(k)を、パラメータφ^(k)に変換する。そして、状態観測器305は、パラメータ変換部307により変換されたパラメータφ^(k)と、第2演算部3032により演算された変換状態量ω(k)を用いて、開路電圧推定値V^(k)を算出する。 The parameter conversion unit 307 refers to a preset table or mathematical expression, and uses λa and λb to change the parameter φ 1 ^ (k) identified by the adaptive identifier 304 to the parameter φ 2 ^ (k). Convert. Then, the state observer 305 uses the parameter φ 2 ^ (k) converted by the parameter conversion unit 307 and the converted state quantity ω 2 (k) calculated by the second calculation unit 3032, so that the open circuit voltage estimated value is obtained. V 0 ^ (k) is calculated.

上記のように本例は、λaの状態変数フィルタを用いて演算する第1演算部3031とλbの状態変数フィルタを用いて演算する第2演算部3032を設ける。そして、本例は、第1演算部3031の演算結果に基づき適応同定器304により演算し、第2演算部3032の演算結果に基づき第2演算部3032により演算する構成をもつ。共通のカットオフ周波数をもつ状態変数フィルタを用いて演算した結果を、適応同定器304及び状態観測器305にそれぞれ出力する場合には、状態変数フィルタの周波数特性がパラメータ同定性能へ与える影響と、状態量推定性能へ与える影響がトレードオフになり、演算精度が充分に発揮できない可能性があった。本例は、適応同定器304において同定に必要な入力周波数帯を無駄に落すことなく、状態観測器305において非線形性域の入力を減衰できるように、各状態変数フィルタに異なるカットオフ周波数を設定できるため、二次電池10の非線形性の影響を受ける事なく、開放電圧を精度よく推定することができる。   As described above, this example includes the first calculation unit 3031 that calculates using the state variable filter of λa and the second calculation unit 3032 that calculates using the state variable filter of λb. In this example, the adaptive identifier 304 calculates based on the calculation result of the first calculation unit 3031, and the second calculation unit 3032 calculates based on the calculation result of the second calculation unit 3032. When the results calculated using the state variable filter having a common cutoff frequency are output to the adaptive identifier 304 and the state observer 305, respectively, the influence of the frequency characteristics of the state variable filter on the parameter identification performance; The effect on the state quantity estimation performance becomes a trade-off, and there is a possibility that the calculation accuracy cannot be fully exhibited. In this example, different cut-off frequencies are set for each state variable filter so that the input in the nonlinearity region can be attenuated in the state observer 305 without wasting the input frequency band necessary for identification in the adaptive identifier 304. Therefore, the open circuit voltage can be accurately estimated without being affected by the non-linearity of the secondary battery 10.

なお、本例は、非線形領域存在割合検出部308により検出された存在割合に応じて、第1状態変数フィルタのカットオフ周波数又は第2状態変数フィルタのカットオフ周波数の少なくとも一方を設定してもよい。また本例は、存在割合に応じて、第1状態変数フィルタのカットオフ周波数又は第2状態変数フィルタのカットオフ周波数の少なくとも一方を、低い周波数に設定してもよい。二次電池10の非線形領域の割合が小さい場合には、λbのフィルタを用いて、演算する必要は必ずしも必要ないため、本例は、二次電池10の状態に応じて、カットオフ周波数の大きさを調整することにより、演算精度を高めることができる。   In this example, at least one of the cutoff frequency of the first state variable filter or the cutoff frequency of the second state variable filter may be set according to the presence ratio detected by the nonlinear region presence ratio detection unit 308. Good. In this example, at least one of the cutoff frequency of the first state variable filter or the cutoff frequency of the second state variable filter may be set to a low frequency according to the existence ratio. When the ratio of the non-linear region of the secondary battery 10 is small, it is not always necessary to perform the calculation using the filter of λb. Therefore, in this example, the cutoff frequency is increased depending on the state of the secondary battery 10. By adjusting the height, the calculation accuracy can be increased.

なお本例の第1演算部3031は本発明の「第1演算手段」に相当し、第2演算部3032は本発明の「第2演算手段」に相当する。   The first calculation unit 3031 of this example corresponds to the “first calculation unit” of the present invention, and the second calculation unit 3032 corresponds to the “second calculation unit” of the present invention.

《実施例1》
電池モデルを用いたシミュレーションにより、本発明の実施形態1に係る電池状態推定装置の効果を検証したので、図9〜図11を用いて、説明する。図9〜11において、上から電流計測値I(k)の変化を示すプロファイル、電圧計測値V(k)の変化を示すプロファイル、同定パラメータ(T1)の変化を示すプロファイル、SOC誤差(%)を示すプロファイル、SOCの推定値の変化を示すプロファイルである。
Example 1
Since the effect of the battery state estimation device according to the first embodiment of the present invention was verified by simulation using a battery model, description will be made with reference to FIGS. 9 to 11, from the top, a profile showing a change in the current measurement value I (k), a profile showing a change in the voltage measurement value V (k), a profile showing a change in the identification parameter (T1), and an SOC error (%). And a profile showing a change in the estimated value of the SOC.

図9は、電流−電圧特性が線形領域でのシミュレーション結果を示し、図10は電流−電圧特性が非線形領域でのシミュレーション結果である。また図11は、図10の点線部分の拡大図である。また図9〜図11において、一点鎖線は、シミュレーション時に求めた真値を示し、破線は、カットオフ周波数をλ1とした場合を示し、実線はカットオフ周波数をλ2とした場合を示す。ここで示したカットオフ周波数は、λ1>λ2という関係とした。   FIG. 9 shows a simulation result when the current-voltage characteristic is in a linear region, and FIG. 10 shows a simulation result when the current-voltage characteristic is in a non-linear region. FIG. 11 is an enlarged view of a dotted line portion of FIG. 9 to 11, the alternate long and short dash line indicates the true value obtained during the simulation, the broken line indicates the case where the cutoff frequency is λ1, and the solid line indicates the case where the cutoff frequency is λ2. The cutoff frequency shown here has a relationship of λ1> λ2.

図9に示すように、電流−電圧特性が線形領域では、カットオフ周波数λ2に比べ、λ1の方が、同定パラメータ(T1)の収束性能もよく、SOC推定誤差も小さくなった。よって、線形領域においては、カットオフ周波数をλ1にすること演算精度が向上した。   As shown in FIG. 9, in the linear region of the current-voltage characteristics, the convergence performance of the identification parameter (T1) is better and the SOC estimation error is smaller at λ1 than at the cutoff frequency λ2. Therefore, in the linear region, the calculation accuracy is improved by setting the cutoff frequency to λ1.

図10及び図11に示すように、IV特性が非線形領域に至るような場合に、線形領域で設定したカットオフ周波数λ1では、SOC推定差も大きく、非線形性の影響を受けて同定パラメータ(T1)が振動的になる場合があった(図11のX参照)。カットオフ周波数λ2では、同定パラメータの振動を抑制した。よって、非線形領域に置いては、カットオフ周波数をλ2にすること演算精度が向上した。   As shown in FIGS. 10 and 11, when the IV characteristic reaches the non-linear region, the SOC estimation difference is large at the cut-off frequency λ1 set in the linear region, and the identification parameter (T1) is affected by the non-linearity. ) Sometimes vibrates (see X in FIG. 11). At the cutoff frequency λ2, the vibration of the identification parameter was suppressed. Therefore, in the nonlinear region, the calculation accuracy is improved by setting the cutoff frequency to λ2.

《実施例2》
実施例1と同様に、本発明の実施形態2に係る電池状態推定装置の効果を検証したので、図12〜図15を用いて、説明する。図12〜15においては、上から電流計測値I(k)の変化を示すプロファイル、電圧計測値V(k)の変化を示すプロファイル、同定パラメータ(T1)の変化を示すプロファイル、SOC誤差(%)を示すプロファイル、SOCの推定値の変化を示すプロファイルである。
Example 2
Similar to Example 1, the effect of the battery state estimation device according to Embodiment 2 of the present invention was verified, and will be described with reference to FIGS. 12 to 15, from the top, a profile showing a change in the current measurement value I (k), a profile showing a change in the voltage measurement value V (k), a profile showing a change in the identification parameter (T1), and an SOC error (%). ) Indicating the change in the estimated value of the SOC.

図12及び図13は、電流−電圧特性が線形領域でのシミュレーション結果であり、図14及び図15は、電流−電圧特性が非線形領域でのシミュレーション結果である。一点鎖線は、シミュレーション時に求めた真値である。図12と図14の破線は、第1演算部3031と第2演算部3032の状態変数フィルタのカットオフ周波数λ1とした場合を示す。実線は第1演算部3031の状態変数フィルタのカットオフ周波数λ1とし、第2演算部3032の状態変数フィルタのカットオフ周波数λ3とした場合を示す。図13と図15の破線は第1演算部3031と第2演算部3032の状態変数フィルタのカットオフ周波数λ2とした場合を示す。実線は第1演算部3031の状態変数フィルタのカットオフ周波数λ2と第2演算部3032の状態変数フィルタのカットオフ周波数λ3とした場合を示す。ここで示したカットオフ周波数は、λ1>λ2>λ3という関係である。   12 and 13 show simulation results when the current-voltage characteristic is in the linear region, and FIGS. 14 and 15 show simulation results when the current-voltage characteristic is in the non-linear region. A one-dot chain line is a true value obtained at the time of simulation. The broken lines in FIGS. 12 and 14 indicate the case where the cutoff frequency λ1 of the state variable filter of the first calculation unit 3031 and the second calculation unit 3032 is set. A solid line indicates a case where the state variable filter cutoff frequency λ1 of the first calculation unit 3031 is set and the state variable filter cutoff frequency λ3 of the second calculation unit 3032 is set. The broken lines in FIGS. 13 and 15 indicate the case where the cutoff frequency λ2 of the state variable filter of the first calculation unit 3031 and the second calculation unit 3032 is set. A solid line indicates a case where the cutoff frequency λ2 of the state variable filter of the first calculation unit 3031 and the cutoff frequency λ3 of the state variable filter of the second calculation unit 3032 are used. The cut-off frequency shown here has a relationship of λ1> λ2> λ3.

図12及び図13より、電流−電圧特性が線形領域では、カットオフ周波数λ2(図13)に比べ、λ1(図12)の方が、同定パラメータ(T1)の収束性能がよく、SOC推定誤差も小さくなった。よって、線形領域においては、カットオフ周波数をλ1にすることで演算精度が向上した。   12 and 13, when the current-voltage characteristic is in the linear region, the convergence performance of the identification parameter (T1) is better at λ1 (FIG. 12) than the cutoff frequency λ2 (FIG. 13), and the SOC estimation error is increased. Became smaller. Therefore, in the linear region, the calculation accuracy is improved by setting the cutoff frequency to λ1.

図14、図15より、電流−電圧特性が非線形領域に至るような場合に、線形領域で設定したカットオフ周波数λ1(図14)では、SOC推定差が大きく、非線形性の影響を受けて同定パラメータ(T1)が振動的になる場合があった(図14のY部分を参照)。カットオフ周波数λ2(図15)では、同定パラメータの振動を抑制した。よって、非線形領域に置いては、カットオフ周波数をλ2にすることで演算精度が向上した。   14 and 15, when the current-voltage characteristic reaches the nonlinear region, the cutoff frequency λ1 (FIG. 14) set in the linear region has a large SOC estimation difference and is identified by the influence of nonlinearity. There was a case where the parameter (T1) became oscillatory (see the Y portion in FIG. 14). At the cutoff frequency λ2 (FIG. 15), the vibration of the identification parameter was suppressed. Therefore, in the non-linear region, the calculation accuracy is improved by setting the cutoff frequency to λ2.

10…二次電池
20…負荷
30…電子制御ユニット
301…電流検出部
302…電圧検出部
303…状態変数フィルタ演算部
3031…第1演算部
3032…第2フィルタ演算部
304…適応同定器
305…状態観測器
306…SOC変換部
307…パラメータ変換部
308…非線形領域存在割合検出部
40…電流センサ
50…電圧センサ
60…温度センサ
70…劣化度センサ
DESCRIPTION OF SYMBOLS 10 ... Secondary battery 20 ... Load 30 ... Electronic control unit 301 ... Current detection part 302 ... Voltage detection part 303 ... State variable filter calculating part 3031 ... 1st calculating part 3032 ... 2nd filter calculating part 304 ... Adaptive identifier 305 ... State observer 306 ... SOC conversion unit 307 ... parameter conversion unit 308 ... nonlinear region existence ratio detection unit 40 ... current sensor 50 ... voltage sensor 60 ... temperature sensor 70 ... degradation degree sensor

Claims (11)

二次電池の電流を、電流計測値として検出する電流検出手段と、
前記二次電池の端子電圧を、電圧計測値として検出する電圧検出手段と、
前記二次電池の電池モデルを定義し、前記電流計測値および前記電圧計測値を、前記電池モデルに基づく状態変数フィルタを用いて、状態量変換して変換状態量を算出し、前記変換状態量から、前記電池モデルに基づく前記二次電池の端子電圧を電圧推定値として推定する端子電圧推定手段と、
前記電圧計測値と前記電圧推定値との差分がゼロに収束するように、前記二次電池のパラメータを同定する同定手段と、
前記二次電池の電流―電圧と特性のうちの非線形領域の存在割合を検出する非線形領域検出手段とを備え、
前記端子電圧推定手段は、前記存在割合に応じて、前記状態変数フィルタのカットオフ周波数を設定することを特徴とする電池状態推定装置。
Current detection means for detecting the current of the secondary battery as a current measurement value;
Voltage detection means for detecting the terminal voltage of the secondary battery as a voltage measurement value;
A battery model of the secondary battery is defined, and the current measurement value and the voltage measurement value are converted into a state quantity using a state variable filter based on the battery model to calculate a conversion state quantity, and the conversion state quantity From the terminal voltage estimation means for estimating the terminal voltage of the secondary battery based on the battery model as a voltage estimated value,
Identification means for identifying parameters of the secondary battery so that a difference between the voltage measurement value and the voltage estimation value converges to zero;
A non-linear region detecting means for detecting the existence ratio of the non-linear region among the current-voltage and characteristics of the secondary battery,
The terminal voltage estimation means sets a cutoff frequency of the state variable filter according to the existence ratio, and the battery state estimation device.
前記端子電圧推定手段は、前記存在割合が所定の存在割合より大きい場合には、前記カットオフ周波数を所定のカットオフ周波数より低く設定することを特徴とする請求項1記載の電池状態推定装置。 2. The battery state estimation device according to claim 1, wherein the terminal voltage estimation unit sets the cutoff frequency lower than a predetermined cutoff frequency when the presence ratio is larger than a predetermined existence ratio. 前記端子電圧推定手段は、
前記存在割合が所定の存在割合より小さい場合には、前記カットオフ周波数を第1のカットオフ周波数に設定し、
前記存在割合が所定の存在割合より大きい場合には、前記カットオフ周波数を第1のカットオフ周波数より低い第2のカットオフ周波数に設定することを特徴とする請求項1又は2記載の電池状態推定装置。
The terminal voltage estimation means includes
If the abundance ratio is less than a predetermined abundance ratio, the cutoff frequency is set to a first cutoff frequency,
3. The battery state according to claim 1, wherein when the existence ratio is larger than a predetermined existence ratio, the cutoff frequency is set to a second cutoff frequency lower than the first cutoff frequency. Estimating device.
前記端子電圧推定手段は、第1の状態変数フィルタを用いて変換状態量を演算する第1の演算部と、前記第1の状態変数フィルタとは異なる第2の状態変数フィルタを用いて変換状態量を演算する第2の演算部とを含むことを特徴とする請求項1〜3のいずれか一項に記載の電池状態推定装置。 The terminal voltage estimation means uses a first state variable filter to calculate a conversion state quantity and a second state variable filter different from the first state variable filter to convert the state The battery state estimation apparatus according to claim 1, further comprising a second calculation unit that calculates the amount. 前記状態変換量と前記パラメータを用いて、前記二次電池の開路電圧を演算する開路電圧演算手段をさらに有し、
前記同定手段は、前記第1の演算部により演算された変換状態量に基づき、前記パラメータを同定し、
前記開路電圧演算手段は、前記第2の演算部により演算された変換状態量に基づき、前記開路電圧を演算することを特徴とする請求項4記載の電池状態推定装置。
Using the state conversion amount and the parameter, further comprising an open circuit voltage calculating means for calculating an open circuit voltage of the secondary battery,
The identification means identifies the parameter based on the conversion state quantity calculated by the first calculation unit,
5. The battery state estimation device according to claim 4, wherein the open circuit voltage calculation means calculates the open circuit voltage based on a conversion state quantity calculated by the second calculation unit.
前記端子電圧推定手段は、前記存在割合に応じて、前記第1の状態変数フィルタ又は前記第2の状態変数フィルタの少なくとも一方のフィルタの前記カットオフ周波数を設定することを特徴とする請求項4又は5記載の電池状態推定装置。 5. The terminal voltage estimating means sets the cut-off frequency of at least one of the first state variable filter and the second state variable filter according to the existence ratio. Or the battery state estimation apparatus of 5. 前記端子電圧推定手段は、前記存在割合が所定の存在割合より大きい場合には、前記第1の状態変数フィルタのカットオフ周波数又は前記第2の状態変数フィルタのカットオフ周波数の少なくとも一方のカットオフ周波数を所定のカットオフ周波数より低くすることを特徴とする請求項4〜6のいずれか一項に記載の電池状態推定装置。 The terminal voltage estimation means is configured to cut off at least one of a cutoff frequency of the first state variable filter and a cutoff frequency of the second state variable filter when the existence ratio is larger than a predetermined existence ratio. The battery state estimation device according to any one of claims 4 to 6, wherein the frequency is set lower than a predetermined cutoff frequency. 前記非線形領域検出手段は、前記電流計測値及び前記電圧計測値に基づき前記非線形領域の存在割合を検出し、
前記端子電圧推定手段は、前記存在割合が所定の存在割合より大きい場合には、前記カットオフ周波数を所定のカットオフ周波数より低くすることを特徴とする請求項1〜7のいずれか一項に記載の電池状態推定装置。
The nonlinear area detecting means detects the existence ratio of the nonlinear area based on the current measurement value and the voltage measurement value,
The said terminal voltage estimation means makes the said cut-off frequency lower than a predetermined cut-off frequency, when the said presence ratio is larger than a predetermined presence ratio, The Claim 1 characterized by the above-mentioned. The battery state estimation device described.
前記非線形領域検出手段は前記二次電池の温度に基づき、前記前記非線形領域の存在割合を検出し、
前記端子電圧推定手段は、前記存在割合が所定の存在割合より大きい場合には、前記カットオフ周波数を所定のカットオフ周波数より低くすることを特徴とする請求項1〜7のいずれか一項に記載の電池状態推定装置。
The non-linear area detecting means detects the existence ratio of the non-linear area based on the temperature of the secondary battery,
The said terminal voltage estimation means makes the said cut-off frequency lower than a predetermined cut-off frequency, when the said presence ratio is larger than a predetermined presence ratio, The Claim 1 characterized by the above-mentioned. The battery state estimation device described.
前記非線形領域検出手段は前記二次電池の劣化度に基づき、前記前記非線形領域の存在割合を検出し、
前記端子電圧推定手段は、前記存在割合が所定の存在割合より大きい場合には、前記カットオフ周波数を所定のカットオフ周波数より低くすることを特徴とする請求項1〜7のいずれか一項に記載の電池状態推定装置。
The non-linear area detecting means detects the existence ratio of the non-linear area based on the degree of deterioration of the secondary battery,
The said terminal voltage estimation means makes the said cut-off frequency lower than a predetermined cut-off frequency, when the said presence ratio is larger than a predetermined presence ratio, The Claim 1 characterized by the above-mentioned. The battery state estimation device described.
前記端子電圧推定手段は、前記電流計測値が所定の電流値より小さい場合には、前記カットオフ周波数を前記所定のカットオフ周波数より高くすることを特徴とする請求項8〜10のいずれか一項に記載の電池状態推定装置。 11. The terminal voltage estimation unit, wherein the cutoff frequency is higher than the predetermined cutoff frequency when the measured current value is smaller than a predetermined current value. The battery state estimation device according to item.
JP2011019822A 2011-02-01 2011-02-01 Battery state estimation device Active JP5707982B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011019822A JP5707982B2 (en) 2011-02-01 2011-02-01 Battery state estimation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011019822A JP5707982B2 (en) 2011-02-01 2011-02-01 Battery state estimation device

Publications (2)

Publication Number Publication Date
JP2012159414A true JP2012159414A (en) 2012-08-23
JP5707982B2 JP5707982B2 (en) 2015-04-30

Family

ID=46840068

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011019822A Active JP5707982B2 (en) 2011-02-01 2011-02-01 Battery state estimation device

Country Status (1)

Country Link
JP (1) JP5707982B2 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015035839A (en) * 2013-08-07 2015-02-19 Tdk株式会社 Power storage system, and electric vehicle including the same
JP2015053826A (en) * 2013-09-09 2015-03-19 パナソニック株式会社 Charge control device, charging apparatus, vehicle and charge control method
US20150219726A1 (en) * 2014-02-03 2015-08-06 GM Global Technology Operations LLC Systems and methods for battery state estimation
DE102015117171A1 (en) 2014-10-09 2016-04-14 Denso Corporation BATTERIEZUSTANDSABSCHÄTZVORRICHTUNG
JP2016167357A (en) * 2015-03-09 2016-09-15 株式会社デンソー Secondary battery control device
CN106896329A (en) * 2017-03-23 2017-06-27 东软集团股份有限公司 The Forecasting Methodology and device of a kind of battery terminal voltage
CN110573892A (en) * 2017-11-02 2019-12-13 株式会社Lg化学 method, apparatus and recording medium for estimating parameters of battery equivalent circuit model
US10663524B2 (en) 2015-03-06 2020-05-26 Denso Corporation Battery state estimating apparatus

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010286445A (en) * 2009-06-15 2010-12-24 Honda Motor Co Ltd Battery status estimation device

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010286445A (en) * 2009-06-15 2010-12-24 Honda Motor Co Ltd Battery status estimation device

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015035839A (en) * 2013-08-07 2015-02-19 Tdk株式会社 Power storage system, and electric vehicle including the same
JP2015053826A (en) * 2013-09-09 2015-03-19 パナソニック株式会社 Charge control device, charging apparatus, vehicle and charge control method
US20150219726A1 (en) * 2014-02-03 2015-08-06 GM Global Technology Operations LLC Systems and methods for battery state estimation
US10527677B2 (en) 2014-10-09 2020-01-07 Denso Corporation Battery state estimation apparatus
DE102015117171B4 (en) 2014-10-09 2019-03-21 Denso Corporation BATTERIEZUSTANDSABSCHÄTZVORRICHTUNG
DE102015117171A1 (en) 2014-10-09 2016-04-14 Denso Corporation BATTERIEZUSTANDSABSCHÄTZVORRICHTUNG
US10663524B2 (en) 2015-03-06 2020-05-26 Denso Corporation Battery state estimating apparatus
JP2016167357A (en) * 2015-03-09 2016-09-15 株式会社デンソー Secondary battery control device
US9827868B2 (en) 2015-03-09 2017-11-28 Denso Corporation Secondary battery control unit
CN106896329A (en) * 2017-03-23 2017-06-27 东软集团股份有限公司 The Forecasting Methodology and device of a kind of battery terminal voltage
CN110573892A (en) * 2017-11-02 2019-12-13 株式会社Lg化学 method, apparatus and recording medium for estimating parameters of battery equivalent circuit model
EP3627167A4 (en) * 2017-11-02 2020-07-29 LG Chem, Ltd. Method, apparatus, and recording medium for estimating parameter of battery equivalent circuit model
US11119157B2 (en) 2017-11-02 2021-09-14 Lg Chem, Ltd. Method, apparatus and recording medium for estimating parameters of battery equivalent circuit model
CN110573892B (en) * 2017-11-02 2021-10-22 株式会社Lg化学 Method, apparatus and recording medium for estimating parameters of battery equivalent circuit model

Also Published As

Publication number Publication date
JP5707982B2 (en) 2015-04-30

Similar Documents

Publication Publication Date Title
JP5707982B2 (en) Battery state estimation device
JP5691592B2 (en) Battery state estimation device
JP5842421B2 (en) Battery state estimation device
JP6507375B2 (en) Battery state estimation device and battery state estimation method
JP5944291B2 (en) Battery parameter estimation apparatus and method
JP6657967B2 (en) State estimation device and state estimation method
US8918300B2 (en) Apparatus and method for battery state of charge estimation
JP4638195B2 (en) Battery degradation degree estimation device
JP5163542B2 (en) Secondary battery input / output possible power estimation device
KR20170138488A (en) METHOD AND APPARATUS FOR DETERMINING CHARGE AND HEALTH CONDITION OF LITHIUM SULFIDE
JP2018096953A (en) Battery state estimation device
JP2015052482A (en) Soundness estimation device and soundness estimation method for battery
JP2016529478A (en) Monitoring the charge stored in the battery
JP2012057998A (en) Charge rate calculation apparatus for secondary battery and charge rate calculation method
JP2006292492A (en) Full charge capacity estimation system for secondary battery
JP2010203854A (en) Device of estimating internal state of secondary battery
JP2012225713A (en) Charge rate estimation device
KR20160002336A (en) Battery status estimation method, system and recording medium for performing the method
JP5029416B2 (en) Secondary battery charging rate estimation device and charging rate estimation method
JP5504657B2 (en) Secondary battery total capacity estimation device
JP5549449B2 (en) Battery state estimation device
JP6541412B2 (en) Charging rate calculation method and charging rate calculation device
JP2010203935A (en) Device of estimating inputtable/outputtable power of secondary battery
JP2010210457A (en) Control device of secondary battery
RU2491566C1 (en) Device for estimation of battery state and method for estimation of battery state

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20131224

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140618

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140624

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140822

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: 20150203

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20150216

R151 Written notification of patent or utility model registration

Ref document number: 5707982

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151