This application claims priority from Japanese Patent Applications No. 2010-002803, filed on Jan. 8, 2010, the entire contents of which are herein incorporated by reference.
1. Technical Field
The present invention relates to a battery characteristic evaluator.
2. Related Art
The ammeter 3 measures a rising or falling value of output current of the battery 1 varying depending on the turning-on/off of the load 2, and the voltmeter 4 measures a rising or falling value of output voltage of the battery 1 varying depending on the turning-on/off of the load 2. Such a specific measurement procedure is described in JP-A-2003-4780.
A circuit constant optimizing unit 6 includes a voltage calculator 6a and a determination unit 6b, optimizes a circuit constant of an equivalent circuit model of the battery 1 as an identification value FV based on the current value data IM measured by the ammeter 3, the voltage value data VM measured by the voltmeter 4, and the equivalent circuit model data EM of the battery 1 which are input from the input unit 5, and outputs the optimized circuit constant of the equivalent circuit model to an output unit 7.
In the circuit constant optimizing unit 6, the current value data IM measured by the ammeter 3, the equivalent circuit model data EM of the battery 1, and a candidate of the circuit constant CC from the determination unit 6b are input to the voltage calculator 6a, and a calculated voltage value VC is calculated and provided to the determination unit 6b.
The voltage value data VM measured by the voltmeter 4 and the calculated voltage value VC calculated by the voltage calculator 6a are input to the determination unit 6b. The measured voltage value data VM and the calculated voltage value VC are compared with each other and it is determined whether the circuit constant is the optimal value. When it is determined that the circuit constant is not the optimal value, a new circuit constant CC is generated from the comparison result and is input to the voltage calculator 6a, and the voltage is calculated again. These processes are repeatedly performed until it is determined that the circuit constant is the optimal value. The identification value FV optimized as the circuit constant of the equivalent circuit model in this way is provided to the output unit 7.
The output unit 7 generates a characteristic curve of the battery 1 based on the identification value FV of the circuit constant of the equivalent circuit model optimized by the circuit constant optimizing unit 6 and displays the generated characteristic curve on a display unit (not shown).
When circuit data shown in
JP-A-2003-4780 discloses the configuration of method and apparatus for measuring internal impedance of a battery.
JP-A-2005-100969 discloses removing an influence of a response voltage due to polarization at the time of measuring internal impedance of a battery.
In a low-frequency region of impedance of the battery 1, Warburg impedance is exhibited due to the influence of diffusion. The Warburg impedance may be calculated as impedance in a frequency domain as shown in
However, a voltage drop curve based on the Warburg impedance could not be reproduced by the combination of a resistor and a capacitor. Nevertheless, when the identification thereof is performed using a DC source, a resistor, and a capacitor, the resistance value of the resistor or the capacitance value of the capacitor becomes a large value that does not correspond to reality as shown in
Exemplary embodiments of the present invention address the above disadvantages and other disadvantages not described above. However, the present invention is not required to overcome the disadvantages described above, and thus, an exemplary embodiment of the present invention may not overcome any disadvantages.
Accordingly, it is an illustrative aspect of the present invention to provide a battery characteristic evaluator which can improve precision of a circuit constant identification value in an equivalent circuit model of a battery in view of Warburg impedance.
According to one or more illustrative aspects of the present invention, there is provided a battery characteristic evaluator configured to identify a circuit constant of an equivalent circuit model based on a current-voltage characteristic of a battery. The evaluator includes: a current waveform divider configured to divide a certain current waveform into a plurality of step functions with a plurality of infinitesimal time intervals and output the step functions; and a circuit constant optimizing unit configured to calculate the optimized circuit constant of the equivalent circuit model, based on the step functions, a measured voltage value, and equivalent circuit model data.
Other aspects and advantages of the present invention will be apparent from the following description, the drawings and the claims.
Hereinafter, exemplary embodiments of the present invention will be described with reference to the accompanying drawings.
In
In the circuit constant optimizing unit 6, a step response calculator 6c and a voltage adder 6d adding the response calculation results V1 to Vn+m of the step response calculator 6c are provided instead of the voltage calculator 6a of
Equivalent circuit model data EM, a candidate of a circuit constant CC from a determination unit 6b, and the step functions I1 to In+m corresponding to the current from the current waveform divider 8 are input to the step response calculator 6c. Accordingly, the step response calculator 6c calculates step response voltages V1 to Vn+m for the current given as the step functions I1 to In+m and inputs the step response voltages V1 to Vn+m as the calculation results to an input terminal of the voltage adder 6d.
The voltage adder 6d adds the step response voltages V1 to Vn+m as the calculation results of the step response calculator 6c to obtain a calculated voltage value VC. Then, the calculated voltage value VC is provided to the determination unit 6b.
Voltage value data VM measured by a voltmeter 4 and the calculated voltage value VC calculated by the voltage adder 6d are input to the determination unit 6b. The measured voltage value VM and the calculated voltage value VC are compared to determine whether the circuit constant is the optimal value as the comparison result. When it is determined that the circuit constant is not the optimal value, a new circuit constant CC is generated from the comparison result and is input to the step response calculator 6c so as to calculate a voltage again. These processes are repeatedly performed until it is determined that the circuit constant is the optimal value. An identification value FV optimized as the circuit constant of the equivalent circuit model in this way is provided to an output unit 7.
The output unit 7 generates a characteristic curve of the battery 1 based on the identification value FV of the circuit constant of the equivalent circuit model optimized by the circuit constant optimizing unit 6 and displays the generated characteristic curve on a display unit (not shown).
The details shown in
I(t)=Ii·u(t−b1)+I2·u(t−b2)+I3·u(t−b3)+ . . . +In·u(t−bn)−In+1·u(t−bn+i)−In+2·u(t−bn+2)− . . . −In+m·u(t−bn+m)=I1·u(t1)+I2·u(t2)+I3·u(t3)+ . . . +In·u(tn)−In+1·u(tn+1)+In+2·u(tn+2)+ . . . +In+m·u(tn+m), (1)
where u(t) is set so that u(ti)=0 (if ti<0) and 1 (if ti≧0) at time ti (where i=1 to n+m).
In Expression 1, Ii (ti) (where I=1 to n) can be expressed as follows by the Laplace transform.
I
i(s)=L(Ii·u(t−bi))=Ii·(1/s) (2)
Similarly, Ii(ti) (where i=n+1 to n+m) can be also expressed as follows by the Laplace transform.
I
i(s)=−L(Ii·u(t−bi))=−Ii·(1/s) (3)
Since these current signals flow in impedance Z(s) and are thus converted into voltages, the voltages Vi(s) (where i=1 to n+m) based on the currents are expressed as follows.
V
i(s)=Z(s)·Ii·1/s (if i=1 to n)
V
i(s)=−Z(s)·Ii·1/s (if i=n+1 to m) (4)
When step current flows in impedance Z, transient voltage response signals Vi(ti) are obtained as follows by the Laplace-transforming Expression (4).
V
i(ti)=L[Vi(s)]=Ii·L[Z(s)·1/s] (if i=1 to n)
V
i(ti)=L[Vi(s)]=−I·L[Z(s)·1/s] (if i=n+1 to m) (5)
Therefore, by recombining the step responses divided into (n+m) steps, a transient voltage response waveform V(t) when a certain current waveform flows in the impedance Z can be expressed by Expression (6).
V(t)=V1(t1)+V2(t2)+V3(t3)+ . . . +Vn(tn)−Vn+1(tn+1)−Vn+2(tn+2) . . . −Vn+m(tn+m) (6)
Accordingly, even when a certain current waveform is input, it is possible to calculate a voltage response of the battery.
According to this configuration, the Warburg impedance can be included in the equivalent circuit and the identification precision of the battery increases, thereby making the current-voltage characteristic closer to reality. Realistic values can be obtained for the circuit constants other than the Warburg impedance.
Although the equivalent circuit model in which the Warburg impedance is connected in parallel has been described in the above-mentioned embodiment, an equivalent circuit in which the Warburg impedance W1 is singly connected in series as shown in
In
In this case, the voltage Vw in the time domain of the Warburg impedance block W1 can be calculated as follows and thus the calculation is simplified.
Vw=(δ√2t)×Ip/Γ(3/2), (7)
where δ represents a constant of diffusion and Γ represents a gamma function.
The total voltage of the equivalent circuit shown in
The method according to the invention can be applied when the input current has a rectangular waveform.
In the above-mentioned embodiment, the current is changed and identified with the measured response voltage, but the voltage may be changed and identified with the measured current value.
According to the above-mentioned invention, it is possible to provide a battery characteristic evaluator which can identify a circuit constant with high precision in an equivalent circuit model of a battery in consideration of the Warburg impedance so as to evaluate a battery characteristic with high precision, and can be suitably used to efficiently analyze various parameters of a battery.
Number | Date | Country | Kind |
---|---|---|---|
2010-002803 | Jan 2010 | JP | national |