CN102983371B - 计算方法、计算系统、计算装置 - Google Patents

Abstract

本发明涉及计算电池的劣化状态的计算方法、计算系统及计算装置。本实施方式的计算方法的目的在于更正确地计算电池的劣化状态。本实施方式的计算方法实现参照测量值数据库和函数信息数据库，以保存于所述函数信息数据库的函数的活性材料的量为变量，对存储于所述测量值数据库的电池电压进行回归计算的运算功能，所述测量值数据库将电池的充电量和此时产生的电池电压对应着存储，所述函数信息数据库存储表示包含于电池内的多种活性材料的充电量与电压的关系的函数。

Description

计算方法、计算系统、计算装置

本申请基于且主张2011年5月31日申请的在先日本专利申请第2011-122965号的优先权，此申请案的全部内容以引用的方式并入本文。

技术领域

本发明的实施方式涉及一种计算电池的劣化状态的计算方法、计算系统、计算装置。

背景技术

计算电池的劣化状态的计算方法、计算系统、计算装置的研究正在进展中。例如，如日本公开公报P2011-75461所记载，利用电池单元的容量值进行电池单元的评价。

但是，在仅用容量的评价中，无法把握是正极或负极中的哪个发生劣化等更正确的状况。

发明内容

根据上述情况，提供一种实施方式，参照测量值数据库和函数信息数据库，实现以保存于所述函数信息数据库的函数的活性材料的量为变量（変数）对存储于所述测量值数据库的电池电压进行回归计算的运算功能，该测量值数据库使电池的充电量和此时产生的电池电压对应着存储，该函数信息数据库存储表示包含于电池内的多种活性材料的充电量与电压的关系的函数。

根据该实施方式，更正确地计算电池的劣化状态。

附图说明

图1是第一实施方式所涉及的计算系统的方框图。

图2是第一实施方式所涉及的电池装置的方框图。

图3是表示与利用了活性材料A及活性材料B的正极的充电量相对的电动势的一例的图。

图4是表示与正极的充电量相对的电位的一例的图。

图5是表示由电流值的增大引起的与充电量相对的负极电位的变化的一例的图。

图6是表示与充电时间相对的充电电压的测量值的一例的图。

图7是表示第一实施方式所涉及的计算装置的动作的流程图。

图8是表示第一实施方式所涉及的回归计算的一例的流程图。

图9是表示第一实施方式所涉及的充电量分配计算的一例的流程图。

图10是表示第一实施方式所涉及的计算装置的变形例的图。

图11是表示第一实施方式所涉及的计算装置的动作的方框图。

图12是表示第一实施方式所涉及的计算装置的动作的方框图。

图13是表示第一实施方式所涉及的计算装置的动作的方框图。

图14是表示第一实施方式所涉及的计算装置的动作的方框图。

具体实施方式

以下，参照附图对本实施方式进行说明。

（第一实施方式）

图1是本实施方式的计算系统的构成图。

图1所示的计算系统是进行电池装置20的剩余容量的计算的计算机系统。在本实施方式中，作为计算系统的构成要素之一的计算装置10可以由根据处理功能通过LAN、因特网等通信网络将装置组合起来的计算装置群构成。

计算装置10是包括CPU100、RAM（RWM）110、通信IF120、输入IF131、显示IF141、ROM150、存储部160及计时器170的结构。此外，也可以具备装载USB存储器等外部存储装置的IF（接口）。计算装置10是执行程序并进行运算的计算机。

计算装置10通过通信IF120从二次电池装置20收集电流值、电压值等数据，并利用收集到的数据进行各种运算处理。

CPU100是将预先写入到ROM150的各程序读出到RAM110，并进行计算处理的运算处理部（微处理器）。CPU100可以使功能合并而由多个CPU群（微型计算机、微型控制器）构成。另外，CPU内可以具备带有RAM功能的内置存储器。

RAM（RWM）110是CPU100执行程序时使用的存储区域，是用作工作区域的存储器。使处理所需的数据暂时存储是优选的。

通信IF120是与二次电池装置进行数据授受的通信装置、通信单元。例如，有路由器。在本实施方式中，通信IF20与二次电池装置20的连接如有线通信般记载，但可以替换为各种无线通信网。另外，通信IF20与二次电池装置20的连接也可以是通过能够单方向或双方向通信的台数的网络来进行的方式。

输入IF130是连接输入部131和计算装置10的接口。可以具有将从输入部131发送来的输入信号转换并转换为CPU100能够识别的信号的输入控制功能。本IF作为端子等不是必须的构成要素，也可以直接与计算装置内的布线连接。

输入部131是计算机装置通常具备的各种键盘、按钮等进行输入控制的输入装置、输入单元。此外，也可以具有通过识别人发出的声音而作为输入信号识别或检测的功能。在本实施方式中，设置在计算装置10的外部，但也可以是组装到计算装置中的方式。

显示IF140是连接显示部141和计算装置10的接口。可以从CPU100经由显示IF141进行显示部130的显示控制，也可以通过绘图板（graphic board）等进行描绘处理的LSI（GPU）进行显示控制。作为显示控制功能，例如有将图像数据解密的译码功能。也可以是不使用IF而直接连接至计算装置10内部的方式。

显示部141是液晶显示器、有机EL显示器、等离子显示器等输出装置、输出单元。此外，也可以具有发出声音的功能。在本实施方式中，设置在计算装置10的外部，但也可以组装到计算装置10的内部。

ROM150是保存回归计算程序151和劣化度计算程序152的程序存储器。利用不能进行数据写入的非一次性存储介质是优选的，但也可以是能够随时进行数据的读出、写入的半导体存储器等存储介质。此外，也可以保存将图像数据在显示部141上显示为人能够识别的文字或图案的显示程序、通过通信IF120向终端30配送电池的劣化信息等内容的程序、每隔预先确定的时间将所取得的数据存储到存储部160的信息登记程序等。

回归计算程序151是使CPU100实现计算构成电池装置20的各电池单元中的每个单元或每个电池组的正极的容量值、负极的容量值、内部电阻值的功能的单元。例如，计算（解析）以下的七个值，（1）构成正极的活性材料A的容量，（2）构成正极的活性材料B的容量，（3）负极的容量，（4）构成正极的活性材料A的充电量，（5）构成正极的活性材料B的充电量，（6）负极的充电量，（7）内部电阻值。

（算式1）

利用这些值，计算充电电压相对于时间的变化特性、正极的电位相对于充电量及或负极的电位相对于充电量特性。具体的动作在后阐述。

回归计算程序151由与以下的各算式相对应的程序群构成。另外，关于各程序的顺序，能够进行各种变更。

充电电压V_Ｃ用电池的起电压（起電圧）V_e、内部电阻的电压V_R根据下面的算式2求出。

V_c＝V_e+V_R                （算式2）

电池的起电压V_e用正极的电位E_c及负极的电位E_a根据下面的算式3求出。

V_e＝E_c-E_a                （算式3）

正极和负极的电位用充电量（q）、初始状态下的正极的容量Q_ic及初始状态下的负极的容量Q_ia根据算式4、算式5求出。

E_c＝f_c(q/Q_ic)                        （算式4）

E_a＝f_a(q/Q_ia)                        （算式5）

这里，对以多种活性材料构成正极或负极的情况进行说明。在此情况下，如图3所示，各种活性材料的电动势示出了不同的特性。计算混合了活性材料A（例如，锰酸锂）和活性材料B（例如，钴酸锂）的复合正极的起电压相对于充电量的特性。计算出的特性图示后如图4。

活性材料A的正极的电位E_cA、活性材料B的正极的电位E_cB用初始状态的活性材料A的容量Q_icA、初始状态的活性材料B的容量_QicB、活性材料A的充电量q_A及活性材料B的充电量q_B并存在算式6、算式7、8、9的关系（图3）。

E_cA＝f_cA(q_A/Q_icA)                       （算式6）

E_cB＝f_cB(q_B /Q_icB)                      （算式7）

f_cA(q_A/Q_cA)＝f_cB(q_B/Q_cB)                （算式8）

q＝q_A+q_B                                （算式9）

因此，混合正极的电位E_c用活性材料A的正极的充电开始时的容量q_A、活性材料A的正极的充电量Q_cA或活性材料B的正极的充电开始时的容量q_B及活性材料B的正极的充电量Q_cB根据算式10求出。

E_c＝f_c(q/Q_ic)＝f_cA(q_A/Q_cA)＝f_cB(q_B/Q_cB)

（算式10）

另外，活性材料A的正极的电位E_cA和活性材料B的充电量q_B是各活性材料表面的电位。因此，活性材料内的锂离子的分布根据活性材料内的锂离子的扩散电阻而变化，所以认为充电量与起电压的关系根据充电电流而变化。但是，在本实施方式中，在使用于正极的活性材料及使用于负极的碳类活性材料中，活性材料内的扩散电阻小，因此作为即使充电电流变化，充电电流与起电压的关系也不会大幅变化的活性材料对待。

另一方面，当在负极中将利用了如钛酸锂的扩散电阻较大的材料作为活性材料利用于负极时，如图5所示，充电量与起电压的关系根据电流值大幅变化，因此不进行与正极同样的近似。

因此，负极电位E_a以算式11表示。

E_a＝f_a(q/Q_ia，I/Q_ia)                        （算式11）

另外，内部电阻的电压V_R用充电电流I和内部电阻R（q）根据算式12、13求出。

V_R＝R(q)×I                        （算式12）

q＝∫Idt                           （算式13）

也就是说，算式2表示为，

Vc＝fc(q/Qic)-fa(q/Qia，I/Qia)+R(q)×I

（算式2A）

这样，充电电压与活性材料的起电压特性及内部电阻之间存在非线性的相关关系。接下来，以活性材料的容量及内部电阻为变量，关于充电电压相对于充电量的特性曲线进行回归计算，计算并决定活性材料的容量及内部电阻。

劣化度计算程序152是使CPU100实现根据通过执行回归计算程序151求出的活性材料的容量和内部电阻的值来计算电池装置20的劣化度的功能的单元。

以各活性材料的容量和内部电阻为变量，将求出劣化度的函数事先存储于函数信息DB163，从函数信息DB163读出函数并对通过本程序得出的各活性材料的容量和内部电阻进行比较，来计算劣化度。

例如，在事先存储容量和内部电阻的基准值，并且通过本程序得到的容量或内部电阻的值比基准值小时，判断为劣化到无法使用的程度，或者也可以根据使用用途，事先存储根据容量和内部电阻值来运算的劣化指标的运算方法（函数），将利用通过本程序得到的容量和内部电阻运算出的劣化指标作为劣化度提示。

存储部160是硬盘驱动器（HDD）等非易失性的存储装置、存储单元。此外，不限于非易失性的存储单元，可以利用闪存器等半导体存储器，也可以是组合了这些半导体存储器和HDD的方式的存储介质。在本实施方式中，以ROM150和存储部160作为不同的存储介质的方式记载，也可以将ROM和存储部合并而作为一个存储部、存储单元。

存储部160保存计测值DB161、计算结果DB162及函数信息DB163。此外，还存储CPU100的计算处理所需要的数据。

测量结果DB161对计测出的电池装置20的电流值、电压值等测量数据每隔预先确定的测量间隔存储数据。例如，将所保存的数据以表形式保存。若进行图示，则绘出如图6所示的图。与所保存的测量结果有关的数据在执行回归计算程序151时利用。

计算结果DB162保存通过执行回归计算程序151的CPU100计算出的值。保存后的值可以由CPU100读出并通过显示IF140显示在显示部141上。也可以是代替本存储单元，存储在云计算系统内的计算装置10外部的存储介质中的方式。

另外，函数信息DB163存储如图4所示那样表示正极电位和充电量的关系的函数、表示图5所示的负极电位与充电量的关系的函数及值。这些表示图5所示的负极电位与充电量的关系是，随着电流值的增大，负极电位相对于充电量变得恒定的区域变少。本函数信息在评价正极、负极的劣化度时利用。

另外，也保存有求出正极、负极相对于充电量的电位时的数据。在正极或负极由多种材料构成的情况下，保存有与各材料中的每种材料相对应的函数。

计时器170是用于计测时间的时钟。CPU100利用计测出的时刻进行电流值、电压值的计测及保存。计算装置10可以是代替计时器170通过网络取得时间的方式，也可以是根据计时器170的计测和从网络获得的时刻来计算时刻的方式。

电池装置20是具备电池和计测电池的电流值、电压值等的控制电路的装置。电池中能够利用各种二次电池。例如，通过VTM21进行电压检测，通过电压检测电路602检测电压。

终端30是装载了上述电池装置20的车辆或定置形电池装置。是包括CPU300、RAM（RWM）310、通信IF320、输入IF331、显示IF341、ROM350、存储部360、计时器370及蓄电池380的结构。另外，也可以具备对蓄电池进行充电的充电器381、在从蓄电池对马达383送电时进行直流-交流转换的变流器382。另外，也可以代替对终端30内进行控制的CPU100而使用可编程控制器385。

各构成要素的功能与计算装置的各构成要素是同样的，所以省略。另外，也可以通过CPU100和CPU300分割着执行各程序的执行。蓄电池380是电池装置20。

接着，主要参照图7、图8对具体的动作的一例进行说明。

为了简单，充电电流值结合定电流充电的情况进行记载。充电电流值无需是一定值，但出于测量误差及计算误差较少的观点，定电流充电下的测量是优选的。

另外，内部电阻值作为相对于充电量q不变化的电阻值对待。

另外，在本实施方式中，记载对各单元中的每个单元测量了电压的情况。在对如电池组般具有多个单元的电池模块使用本计算方法的情况下，因为各单元中的每个单元的劣化的进展状况都不同，所以对每个单元测量充电电压的时间变化是优选的。（步骤S101）。

这样，将测量到的充电电压保存在测量结果DB161中。在从充电开始到到达充电终止电压为止的充电时间t_c期间得到N个测量值。刚刚充电之后是过渡期，充电电压尚未稳定，因此在该过渡期ts期间得到的ns个数据未用在用于计算活性材料的容量及内部电阻的回归计算中。

计算装置10取得各单元的充电电压的测量值。将所取得的值暂时保存在RAM110或存储在存储部160中。

CPU100根据ROM150执行回归计算程序151，并通过回归计算对非线性微分方程式的解进行解析（步骤S102）。

因为是进行定电流充电，所以从充电开始起的充电量q_c根据算式13求出：

q_c＝I·t

这样，得到与充电量q对应的测量电压、测量值（V₁，q_c1），（V₂，q_c2），…，（V_N，q_cN）。

CPU100将得到的值暂时保存在RAM110中或存储在存储部160中。

用上述测量值进行回归计算。在进行回归计算时使用的残差平方和以下式表示。

$S=\underset{n=\mathrm{ns}+1}{\overset{N}{\mathrm{\Σ}}}{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/\mathrm{Qc})-f_a((q_{\mathrm{cn}}+q_0^a)/\mathrm{Qa},I/\mathrm{Qa})))}^2$ （算式14）

充电开始时的正极的充电量

充电开始时的负极的充电量

充电开始时的充电量在回归计算时是未知数，所以上述充电开始时的正极的充电量及负极的充电量也是未知数。

在本实施方式中，在正极为活性材料A和活性材料B的复合正极的情况下，回归计算的未知数为，

（算式15）。

因此，求出下面的联立方程式的解（步骤S1022）。

$\left\{\begin{array}{c}\∂S/{\∂Q}_{\mathrm{cA}}=0\\ \∂S/{\∂Q}_{\mathrm{cB}}=0\\ \∂S/{\∂Q}_a=0\\ \∂S/{\∂q}_0^{\mathrm{cA}}=0\\ \∂S/{\∂q}_0^{\mathrm{cB}}=0\\ \∂S/{\∂q}_0^a=0\\ \∂S/\∂R=0\end{array}\right.$

（算式16）

作为初始值，利用适当的值，例如上次测量时的值（步骤S1021）。

$\left\{\begin{array}{c}\∂S/{\∂Q}_{\mathrm{cA}}=0\\ \∂S/{\∂Q}_{\mathrm{cB}}=0\\ \∂S/{\∂Q}_a=0\\ \∂S/{\∂q}_0^{\mathrm{cA}}=0\\ \∂S/{\∂q}_0^{\mathrm{cB}}=0\\ \∂S/{\∂q}_0^a=0\\ \∂S/\∂R=0\end{array}\right.$ （算式17）

下面的步骤的各值以下式求出（步骤S1023）。

$\left\{\begin{array}{c}{Q}_{\mathrm{cA}}\←Q_{\mathrm{cA}}+{\mathrm{\δQ}}_{\mathrm{cA}}\\ Q_{\mathrm{cB}}\←Q_{\mathrm{cB}}+{\mathrm{\δQ}}_{\mathrm{cB}}\\ Q_a\←Q_a+{\mathrm{\δQ}}_a\\ q_0^{\mathrm{cA}}\←q_0^{\mathrm{cA}}+{\mathrm{\δq}}_0^{\mathrm{cA}}\\ q_0^{\mathrm{cB}}\←q_0^{\mathrm{cB}}+{\mathrm{\δq}}_0^{\mathrm{cB}}\\ q_0^a\←q_0^a+{\mathrm{\δq}}_0^a\\ R\←R+\mathrm{\δR}\end{array}\right.$ （算式18）

此时，

（算式19）

求解下式而获得。在本实施方式中，用牛顿法进行说明，但可以用马奎特（Levenberg-Marquardt）法等其他数值解析法代替牛顿法。

（算式20）

反复运算，直到求得的值满足以下的收敛条件（收敛半径）（步骤S1024，步骤S1025）。

$\left\{\begin{array}{c}{\mathrm{\&delta;Q}}_{\mathrm{cA}}<{\mathrm{\&delta;}}_{\mathrm{QcA}}\\ {\mathrm{\&delta;Q}}_{\mathrm{cB}}<{\mathrm{\&delta;}}_{\mathrm{QcB}}\\ {\mathrm{\&delta;Q}}_a<{\mathrm{\&delta;}}_{\mathrm{Qa}}\\ {\mathrm{\&delta;q}}_0^{\mathrm{cA}}<{\mathrm{\&delta;}}_{q_0^{\mathrm{cA}}}\\ {\mathrm{\&delta;q}}_0^{\mathrm{cB}}<{\mathrm{\&delta;}}_{q_0^{\mathrm{cB}}}\\ {\mathrm{\&delta;q}}_0^a<{\mathrm{\&epsiv;}}_{q_0^a}\\ \mathrm{\&delta;R}<{\mathrm{\&delta;}}_R\end{array}\right.$ （算式21）

这样，反复运算后的结果保存在存储部160中（步骤S103）。

接着，CPU100接下来执行劣化度计算程序152。

通过对由上述回归计算程序151计算出的容量及内部电阻和预先存储的劣化基准进行比较或通过预先存储的运算方法（函数），用容量和内部电阻来计算劣化度（步骤S104）。

CPU100执行显示程序等显示控制单元并将计算出的结果显示在显示部141上。

另外，也可以是将在步骤102中运算处理过的充电电压和内部电阻直接原样地显示在显示部141上的方式。

接着，例示出步骤S1022的运算处理。

第一，在正极或负极为单极的情况下，例如对算式27的第3列第3行求解后，成为

$\frac{{\&PartialD;}^{2}S}{{{\&PartialD;Q}_a}^2}=\frac{{\&PartialD;}^2}{{{\&PartialD;Q}_a}^2}\underset{n=\mathrm{ns}+1}{\overset{N}{\mathrm{\&Sigma;}}}{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/Q_a,I/Q_a))-R\&times;I)}^2$

$=\underset{n=\mathrm{ns}+1}{\overset{N}{\mathrm{\&Sigma;}}}\frac{{\&PartialD;}^2}{{{\&PartialD;Q}_a}^2}{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/Q_a,I/Q_a))-R\&times;I)}^2$ （算式22），利用小变动(小変異)ΔQc计算微分，而成为，

$\underset{n=\mathrm{ns}+1}{\overset{N}{\mathrm{\&Sigma;}}}\frac{{\&PartialD;}^2}{{\&PartialD;Q}_a}{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/Q_a,I/Q_a))-R\&times;I)}^2$

$=\underset{n=\mathrm{ns}+1}{\overset{N}{\mathrm{\&Sigma;}}}(\frac{{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/(Q_a+{\mathrm{\&delta;Q}}_a),I/(Q_a+{\mathrm{\&delta;Q}}_a)))-R\&times;I)}^2}{{{\mathrm{\&delta;Q}}_a}^2}$

$-2(\frac{{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/\left(Q_a\right),I/\left(Q_a\right)))-R\&times;I)}^2}{{{\mathrm{\&delta;Q}}_a}^2}$

$\frac{{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/(Q_a-{\mathrm{\&delta;Q}}_a),I/(Q_a-{\mathrm{\&delta;Q}}_a)))-R\&times;I)}^2}{{{\mathrm{\&delta;Q}}_a}^2})$ （算式23）。

同样地，对第1行第7列的算式计算后，表示为

$\underset{n=\mathrm{ns}+1}{\overset{N}{\mathrm{\&Sigma;}}}\frac{{\&PartialD;}^2}{{\&PartialD;Q}_a\&PartialD;R}{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/Q_a,I/Q_a))-R\&times;I)}^2$

$=\underset{n=\mathrm{ns}+1}{\overset{N}{\mathrm{\&Sigma;}}}(\frac{{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/(Q_a+{\mathrm{\&delta;Q}}_a),I/(Q_a+{\mathrm{\&delta;Q}}_a)))-(R+\mathrm{\&delta;R})\&times;I)}^2}{{4\mathrm{\&delta;Q}}_a\mathrm{\&delta;R}}$

$-(\frac{{(V_n-f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/(Q_a+{\mathrm{\&delta;Q}}_a),I/(Q_a+{\mathrm{\&delta;Q}}_a)))-(R-\mathrm{\&delta;R})\&times;I)}^2}{{4\mathrm{\&delta;Q}}_a\mathrm{\&delta;R}}$

$-(\frac{{(V_n-f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/(Q_a+{\mathrm{\&delta;Q}}_a),I/(Q_a-{\mathrm{\&delta;Q}}_a)))-(R+\mathrm{\&delta;R})\&times;I)}^2}{{4\mathrm{\&delta;Q}}_a\mathrm{\&delta;R}}$

$\left(\frac{{(V_n-(f_c((q_{\mathrm{cn}}+q_0^c)/Q_c)-f_a((q_{\mathrm{cn}}+q_0^a)/(Q_a-{\mathrm{\&delta;Q}}_a),I/(Q_a-{\mathrm{\&delta;Q}}_a)))-(R-\mathrm{\&delta;R})\&times;I)}^2}{{4\mathrm{\&delta;Q}}_a\mathrm{\&delta;R}}\right)$ （算式24）。

对其他项也同样进行计算。

接着，利用图10对正极为复合件的情况的步骤S1022的运算的前处理进行说明。

首先，设定初始值（步骤S10211）。

q_cn←q_cns+1、q_A，q_B

接着，对q_A、q_B生成联立方程式，并求出解（步骤S10212，S10213）。

$\frac{{\&PartialD;f}_{\mathrm{cA}}(q_{\mathrm{An}}/Q_{\mathrm{cA}})}{\&PartialD;q}{\mathrm{\&delta;q}}_A-\frac{{\&PartialD;f}_{\mathrm{cB}}(q_{\mathrm{Bn}}/Q_{\mathrm{cB}})}{\&PartialD;q}{\mathrm{\&delta;q}}_B=-f_{\mathrm{cA}}(q_{\mathrm{An}}/Q_{\mathrm{cA}})+f_{\mathrm{cB}}(q_{\mathrm{Bn}}/Q_{\mathrm{cB}})$

$\left[A_c\right]=\left(\begin{array}{cc}\frac{{\&PartialD;f}_{\mathrm{cA}}((q_{\mathrm{An}}+q_0^{\mathrm{cA}})/Q_{\mathrm{cA}})}{\&PartialD;q}& \frac{{\&PartialD;f}_{\mathrm{cB}}((q_{\mathrm{Bn}}+q_0^{\mathrm{cB}})/Q_{\mathrm{cB}})}{\&PartialD;q}\\ 1& 1\end{array}\right)$

$\left(B_c\right)=\left(\begin{array}{c}-f_{\mathrm{cA}}((q_{\mathrm{An}}+q_0^{\mathrm{cA}})/Q_{\mathrm{cA}})+f_{\mathrm{cB}}((q_{\mathrm{Bn}}+q_0^{\mathrm{cB}})/Q_{\mathrm{cB}})\\ q_{\mathrm{cn}}-q_{\mathrm{An}}-q_{\mathrm{Bn}}\end{array}\right)$ （算式25）

接着，判定解是否落入预先确定的收敛半径内（步骤S10214）。

$\begin{array}{c}{\mathrm{\&delta;q}}_A<{\mathrm{\&epsiv;}}_{q_A}\\ \mathrm{and}\\ {\mathrm{\&delta;q}}_B<{\mathrm{\&epsiv;}}_{q_B}\end{array}$ （算式26）

未落入的情况下，设定新的初始值，并求出q_A、q_B（步骤S10215）。

落入的情况下（步骤S10214为是），进行对下一测量点的计算（步骤S10216，S10217）。

得到的解以下式表示。

$\left(\begin{array}{ccc}{q}_{\mathrm{cns}+1}& q_{\mathrm{Ans}+1}& q_{\mathrm{Bns}+1}\\ q_{\mathrm{cns}+2}& q_{\mathrm{Ans}+2}& q_{\mathrm{Bns}+2}\\ .& .& .\\ .& .& .\\ .& .& .\\ q_{\mathrm{cN}-1}& q_{\mathrm{AN}-1}& q_{\mathrm{BN}+1}\\ q_{\mathrm{cN}}& q_{\mathrm{AN}}& q_{\mathrm{BN}}\end{array}\right)$

（算式27）

这里，正极电压能够取得为：

(q_Ans+1q_Ans+2…q_AN)        （算式28）。

该数列可以记载为

$\left(\begin{array}{cccc}{q}_{\mathrm{Ans}+1}(Q_{\mathrm{cA}},Q_{\mathrm{cB}},q_0^{\mathrm{cA}},q_0^{\mathrm{cB}})& q_{\mathrm{Ans}+2}(Q_{\mathrm{cA}},Q_{\mathrm{cB}},q_0^{\mathrm{cA}},q_0^{\mathrm{cB}})& ...& q_{\mathrm{AN}}(Q_{\mathrm{cA}},Q_{\mathrm{cB}},q_0^{\mathrm{cA}},q_0^{\mathrm{cB}})\end{array}\right)$ （算式29）。

于是，若对进行例示，则成为

（算式30）。对于关系到其他的正极的项也可以同样地求出。

这样，通过本实施方式，可以计算正极的电位、负极的电位及充电电压。

这样，通过更精密地求出正极的电池、负极的电位，能够检测仅用充电电压的劣化判断无法判断的仅正极、负极的劣化在进展的情况，所以能够进行更正确的劣化判定。

进而，通过根据电池的特性使一部分值一定，也能够实现处理的高速化。

<变形例>

参照图10至图14对第一实施方式的变形例进行说明。

在图10、11中，在终端30内具备计算装置10。在此情况下，无需在上位侧设置服务器，就能够执行通信网络1的信息通信量的降低。

在图12中，是计算设备通过因特网及通信机汇集终端30的信息并在计算装置10评价电池状态的实施方式。能够进行数据的一元管理，能够更正确地进行电池状态的评价。另外，对可靠性的保持也是有益的。

在图13中，通过无线通信进行来自终端30的电池信息的收集。在电池特性的计算中能够每隔一定时间进行是优选的。

在图14中，是在快速充电器等充电器中具备计算装置10的实施方式。

虽已对本发明的若干实施方式进行了描述，但这些实施方式仅是作为示例而提出的，并不意欲限定发明的范畴。当然，本文中所描述的新颖的方法和系统可通过其他各种方式具体化，而且，可在不偏离本发明的精神的前提下对本文中所描述的方法和系统进行各种省略、替换和更改。随附的权利要求书及其等效物意在涵盖属于本发明的范畴及精神的这些方式或修改。

Claims (17)

1.一种计算方法，使计算电池电压的计算装置实现运算功能，包括：

参照测量值数据库和函数信息数据库，以保存于所述函数信息数据库的函数的活性材料的量为变量，对存储于所述测量值数据库的电池电压进行回归计算，所述测量值数据库将电池的充电量和此时产生的电池电压对应着存储，所述函数信息数据库存储表示包含于电池内的多种活性材料的充电量与电压的关系的函数。

2.根据权利要求1所述的计算方法，还以所述电池的内部电阻值和根据充电电流值求出的从所述电池的充电开始起的到测量时为止的充电量为变量，计算所述电池的充电电压的时间特性。

3.根据权利要求1所述的计算方法，还以活性材料的充电量的变化量为变量，计算所述电池的充电电压的时间特性。

4.根据权利要求2所述的计算方法，其特征在于，

所述充电电流值是定电流充电时的电流值。

5.根据权利要求1所述的计算方法，其特征在于，

所述电池包括含有多种活性材料的正极及/或含有多种活性材料的负极，

参照保存所述活性材料的电位-充电量特性的函数信息数据库，计算含有所述多种活性材料的正极的正极电压及/或含有所述多种活性材料的负极的负极电压。

6.根据权利要求1所述的计算方法，其特征在于，

还计算构成电池的各电池单元中每个电池单元或每个电池组的正极的容量值、负极的容量值、内部电阻值。

7.根据权利要求1所述的计算方法，其特征在于，

计算出构成正极的第一活性材料的容量、构成正极的第二活性材料的容量、负极的容量、构成正极的第一活性材料的充电量、构成正极的第二活性材料的充电量及负极的充电量的值，并利用这些值计算充电电压相对于时间的变化特性、正极相对于充电量的电位特性及/或负极相对于充电量的电位特性。

8.根据权利要求1所述的计算方法，其特征在于，

在负极中用钛酸锂作为活性材料。

9.一种计算装置，包括：

测量值数据库，将电池的充电量和此时产生的电池电压对应着存储；

函数信息数据库，存储表示包含于电池内的多种活性材料的充电量与电压的关系的函数；以及

运算单元，参照所述测量值数据库和所述函数信息数据库，以保存于所述函数信息数据库的函数的活性材料的量为变量，对存储于所述测量值数据库的电池电压进行回归计算。

10.根据权利要求9所述的计算装置，其特征在于，

所述电池包括含有多种活性材料的正极及/或含有多种活性材料的负极，

所述运算单元还参照保存所述活性材料的电位-充电量特性的函数信息数据库，计算含有所述多种活性材料的正极的正极电压及/或含有所述多种活性材料的负极的负极电压。

11.根据权利要求9所述的计算装置，其特征在于，

在负极中用钛酸锂作为活性材料。

12.根据权利要求9所述的计算装置，其特征在于，

所述运算单元计算构成正极的第一活性材料的容量、构成正极的第二活性材料的容量、负极的容量、构成正极的第一活性材料的充电量、构成正极的第二活性材料的充电量及负极的充电量的值，并用这些值计算充电电压相对于时间的变化特性、正极相对于充电量的电位特性及/或负极相对于充电量的电位特性。

13.一种电池状态评价用的计算装置，其特征在于，具备：

测量值存储单元，将电池的充电量和此时产生的电池电压对应着存储；

函数存储单元，存储表示包含于电池内的多种活性材料的充电量与电压的关系的函数；以及

运算单元，以存储于所述函数存储单元的函数的活性材料的量为变量，对存储于所述测量值存储单元的电池电压进行回归计算。

14.根据权利要求13所述的电池状态评价用的计算装置，其特征在于，

所述测量值存储单元包括根据电池的通电电流求出所述电池的充电量的单元，所述运算单元将电池的内部电阻作为未知数处理。

15.根据权利要求13所述的电池状态评价用的计算装置，其特征在于，

存储于所述函数存储单元的函数除了活性材料的充电量以外，还以活性材料的充电量的变化量为变量。

16.根据权利要求14所述的电池状态评价用的计算装置，其特征在于，

所述测量值存储单元存储对电池进行定电流充电时的电池电压。

17.根据权利要求13所述的电池状态评价用的计算装置，其特征在于，

所述函数存储单元关于在正极或负极任一个中含有多种活性材料的电池，按其每种活性材料存储函数。