CN115219902B - A method and system for quickly testing the life of a power battery - Google Patents

A method and system for quickly testing the life of a power battery Download PDF

Info

Publication number
CN115219902B
CN115219902B CN202210789672.XA CN202210789672A CN115219902B CN 115219902 B CN115219902 B CN 115219902B CN 202210789672 A CN202210789672 A CN 202210789672A CN 115219902 B CN115219902 B CN 115219902B
Authority
CN
China
Prior art keywords
battery
variable
sequence
cycle life
life
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210789672.XA
Other languages
Chinese (zh)
Other versions
CN115219902A (en
Inventor
张承慧
李世鹏
商云龙
朱昱豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
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 Shandong University filed Critical Shandong University
Priority to CN202210789672.XA priority Critical patent/CN115219902B/en
Publication of CN115219902A publication Critical patent/CN115219902A/en
Application granted granted Critical
Publication of CN115219902B publication Critical patent/CN115219902B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Tests Of Electric Status Of Batteries (AREA)
  • Secondary Cells (AREA)

Abstract

本发明属于动力电池领域,公开了一种动力电池寿命快速测试方法及系统,包括:获取待测试动力电池的变量曲线,基于获取的变量曲线得到横向和纵向变量序列;将横向和纵向变量序列作为自变量,基于核典型相关分析方法将自变量拟合到电池循环寿命,引入递归核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数;将最大相关系数高于阈值的电池循环寿命对应的序列结合机器学习模型进行测试得到动力电池寿命测试结果。本发明避免了对先验知识的要求,与传统测试方法相比,在电池循环寿命测试中表现出良好的有效性并提高了鲁棒性,且极大缩短了测试时间。

The present invention belongs to the field of power batteries, and discloses a method and system for rapid testing of power battery life, including: obtaining a variable curve of a power battery to be tested, and obtaining a horizontal and vertical variable sequence based on the obtained variable curve; using the horizontal and vertical variable sequence as independent variables, fitting the independent variables to the battery cycle life based on the kernel canonical correlation analysis method, and introducing a recursive kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variable and the battery cycle life; combining the sequence corresponding to the battery cycle life with a maximum correlation coefficient higher than a threshold with a machine learning model to obtain a power battery life test result. The present invention avoids the requirement for prior knowledge, and compared with traditional testing methods, it shows good effectiveness and improved robustness in battery cycle life testing, and greatly shortens the testing time.

Description

一种动力电池寿命快速测试方法及系统A method and system for quickly testing the life of a power battery

技术领域Technical Field

本发明属于动力电池领域,尤其涉及一种动力电池寿命快速测试方法及系统。The present invention belongs to the field of power batteries, and in particular relates to a method and system for quickly testing the life of a power battery.

背景技术Background Art

电池寿命是指可用容量下降至初始值80%时的充放电次数。对电池进行循环寿命测试,能够进一步了解电池特性,论证电池是否达到设计目标,在使用电池的过程中实现更好的管控,同时评估电池满足不同应用场景需求的能力。Battery life refers to the number of charge and discharge cycles until the available capacity drops to 80% of the initial value. Cycle life testing of batteries can further understand battery characteristics, prove whether the battery meets the design goals, achieve better control during battery use, and evaluate the battery's ability to meet the needs of different application scenarios.

目前,锂离子电池循环寿命测试主要采用国标GB/T 31484-2015《电动汽车用动力蓄电池循环寿命要求及试验方法》。该标准中规定,进行电池标准循环寿命测试时,需将电池连续循环充放电至其寿命终止条件,即“一测到底”。锂离子电池充放电循环次数可达几千甚至一万次以上,测试时间可达1年乃至更长时间,这项测试通常十分耗时。现有技术公开了一种电池测试方法:控制电池进行N轮次充放电操作,每轮次充放电操作包括多次充电操作和多次放电操作,多次充电操作至少基于两种电流规格进行,多次放电操作至少基于两种电流规格进行;在控制电池进行N轮次充放电操作后,将测试次数记录加1;重新执行确定电池的电池容量值的步骤,直至电池容量值小于预设容量值。现有技术公开了一种电池测试方法:获取电流充放电曲线,基于所述电池使用电压上下限,基于该电流充放电曲线进行电压区间划分;对相同的电芯样品平行在不同电压区间内进行循环测试,对应获得不同电压区间的循环衰减;对所述不同电压区间的循环衰减进行融合,获得循环总衰减。然而以上方法获得电池循环寿命的需要通过全周期、长时间、高成本的测试。At present, the cycle life test of lithium-ion batteries mainly adopts the national standard GB/T 31484-2015 "Requirements and Test Methods for Cycle Life of Power Storage Batteries for Electric Vehicles". The standard stipulates that when conducting a standard battery cycle life test, the battery needs to be continuously charged and discharged until its life end condition, that is, "test to the end". The number of charge and discharge cycles of lithium-ion batteries can reach thousands or even more than 10,000 times, and the test time can reach 1 year or even longer. This test is usually very time-consuming. The prior art discloses a battery testing method: controlling the battery to perform N rounds of charge and discharge operations, each round of charge and discharge operations includes multiple charging operations and multiple discharging operations, multiple charging operations are performed based on at least two current specifications, and multiple discharging operations are performed based on at least two current specifications; after controlling the battery to perform N rounds of charge and discharge operations, the number of tests is recorded plus 1; re-execute the step of determining the battery capacity value of the battery until the battery capacity value is less than the preset capacity value. The prior art discloses a battery testing method: obtaining a current charge and discharge curve, dividing the voltage interval based on the upper and lower limits of the battery voltage, and performing cycle tests on the same battery sample in parallel in different voltage intervals to obtain the cycle attenuation of the different voltage intervals; and fusing the cycle attenuation of the different voltage intervals to obtain the total cycle attenuation. However, the above method requires full-cycle, long-term, and high-cost testing to obtain the battery cycle life.

基于预测的电池寿命快速测试方法是解决上述问题优秀方案,即“以估代测”。通过以预测代替全周期测试,可省去大量时间,能够实现电池循环寿命的高效准确评估。现有技术公开了一种预测方法,通过将电池进行不同循环次数的短期循环性能测试,记录循环次数和容量保持率,然后将不同循环次数后的电池进行拆解,利用X射线衍射法测试石墨负极材料的石墨化度,根据循环次数、容量保持率和石墨化度三种数据进行测试。此种方法需要拆解破坏电池,存在一定局限性。现有技术公开了一种锂离子电池循环寿命的测试方法,在电池表面安装压力传感器,记录电池一定循环次数内的容量信息,根据循环次数、放电容量和压力传感器中的电压数据进行拟合,实现对电池寿命的测试。此方法需要额外的压力传感器,实际中会增加成本,存在一定局限性。现有技术还公开了一种锂离子循环寿命的快速测试方法,通过将电池置于不同倍率下进行500次循环测试,拟合方程得到测试方程实现电池的寿命测试。此方法无需精密的测试设备和复杂的理论计算,但依然需要较长时间的循环测试,实用性不高。电池寿命测试方法有助于电池的生产、优化和开发,但是动力电池是一个强非线性系统,对环境高度敏感,循环寿命测试极其困难。新能源汽车产业要求动力电池能够快速进行技术迭代和产品升级,但现有电池的循环寿命快速测试方法缺失,严重制约新能源汽车行业的快速高质量发展。因此,如何缩短寿命测试时间,快速准确高效地测试并评估电池循环寿命,已成为突破动力电池及其相关产业快速发展关键技术瓶颈的重要手段。The prediction-based rapid battery life test method is an excellent solution to the above problems, that is, "estimate instead of test". By replacing the full cycle test with prediction, a lot of time can be saved, and the efficient and accurate evaluation of the battery cycle life can be achieved. The prior art discloses a prediction method, which performs a short-term cycle performance test on the battery with different cycle times, records the cycle times and capacity retention rate, and then disassembles the battery after different cycle times, uses X-ray diffraction to test the graphitization degree of the graphite negative electrode material, and tests according to the three data of cycle times, capacity retention rate and graphitization degree. This method requires disassembly and destruction of the battery, and has certain limitations. The prior art discloses a test method for the cycle life of a lithium-ion battery, installs a pressure sensor on the surface of the battery, records the capacity information of the battery within a certain number of cycles, and fits the cycle times, discharge capacity and voltage data in the pressure sensor to achieve the test of the battery life. This method requires an additional pressure sensor, which will increase the cost in practice and has certain limitations. The prior art also discloses a rapid test method for the cycle life of a lithium-ion battery, which performs 500 cycle tests on the battery at different rates, and obtains a test equation to achieve the battery life test by fitting the equation. This method does not require sophisticated testing equipment and complex theoretical calculations, but it still requires a long cycle test, so it is not very practical. Battery life test methods are helpful for battery production, optimization and development, but power batteries are a strongly nonlinear system that is highly sensitive to the environment, and cycle life testing is extremely difficult. The new energy vehicle industry requires that power batteries be able to quickly iterate technology and upgrade products, but the lack of rapid cycle life test methods for existing batteries seriously restricts the rapid and high-quality development of the new energy vehicle industry. Therefore, how to shorten the life test time and quickly, accurately and efficiently test and evaluate the battery cycle life has become an important means to break through the key technical bottlenecks in the rapid development of power batteries and related industries.

目前已存在许多用于锂电池循环寿命测试的模型测试方法,通常可以分为机理模型、经验模型和机器学习模型。目前用于动力电池寿命测试的许多机理模型为基于伪二维的电化学模型,其考虑了由副反应引起的容量退化,但它们无法考虑电池的不一致性。与机理模型不同,经验模型直接拟合容量下降轨迹而忽略了电池内部机制;为了根据最新的电池状态和工作条件更新模型参数,卡尔曼滤波器(KF)和局部滤波器(PF)被广泛使用;然而,由于这些经验模型与收集的老化曲线拟合,它们通常只能在相似条件下应用,这意味着泛化能力较差。机器学习模型通过将一组提取的电池特征映射到基于大量数据的电池循环寿命来进行测试。然而,大部分机器学习模型通常输入的特征较少,在实际应用中可能会导致鲁棒性较弱;此外,这些特征通常集中在整个数据集的几个点上,数据集的其余部分没有被充分利用,导致准确性降低;此外,手动找到合适的数据特征集通常需要花费大量精力。There are many model testing methods for lithium battery cycle life testing, which can generally be divided into mechanism models, empirical models and machine learning models. Many mechanism models currently used for power battery life testing are based on pseudo-two-dimensional electrochemical models, which take into account capacity degradation caused by side reactions, but they cannot take into account battery inconsistencies. Unlike mechanism models, empirical models directly fit the capacity decline trajectory and ignore the internal mechanism of the battery; in order to update model parameters according to the latest battery status and working conditions, Kalman filters (KF) and local filters (PF) are widely used; however, since these empirical models fit the collected aging curves, they can usually only be applied under similar conditions, which means that the generalization ability is poor. Machine learning models are tested by mapping a set of extracted battery features to the battery cycle life based on a large amount of data. However, most machine learning models usually have fewer input features, which may lead to weak robustness in practical applications; in addition, these features are usually concentrated on a few points in the entire data set, and the rest of the data set is not fully utilized, resulting in reduced accuracy; in addition, manually finding a suitable data feature set usually takes a lot of effort.

发明内容Summary of the invention

为了解决传统测试方法的周期长、成本高的技术问题,本发明提出了一种动力电池寿命快速测试方法及系统,基于获取的变量曲线得到横向和纵向变量序列;将横向和纵向变量序列作为自变量,基于核典型相关分析方法将自变量拟合到电池循环寿命,引入递归核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数以及对应的核相关系数分析解决方案,在电池循环寿命测试中表现出良好的有效性,极大缩短了测试时间。In order to solve the technical problems of long cycle and high cost of traditional testing methods, the present invention proposes a rapid testing method and system for power battery life, which obtains horizontal and vertical variable sequences based on the acquired variable curves; takes the horizontal and vertical variable sequences as independent variables, fits the independent variables to the battery cycle life based on the kernel canonical correlation analysis method, introduces the recursive kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variable and the battery cycle life and the corresponding kernel correlation coefficient analysis solution, which shows good effectiveness in battery cycle life testing and greatly shortens the test time.

为了实现上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solution:

本发明的第一个方面提供一种动力电池寿命快速测试方法,采用估计代替待测试的方法,包括如下步骤:The first aspect of the present invention provides a method for quickly testing the life of a power battery, which adopts a method of estimating instead of testing, and comprises the following steps:

获取待测试动力电池的变量曲线;Obtaining a variable curve of the power battery to be tested;

基于获取的变量曲线得到横向和纵向变量序列;Based on the obtained variable curves, horizontal and vertical variable sequences are obtained;

将横向和纵向变量序列作为自变量,基于核典型相关分析方法将自变量拟合到电池循环寿命,引入递归核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数;The horizontal and vertical variable sequences are used as independent variables, and the independent variables are fitted to the battery cycle life based on the kernel canonical correlation analysis method. The recursive kernel canonical correlation analysis method is introduced to obtain the maximum correlation coefficient between the independent variable and the battery cycle life.

将最大相关系数高于阈值的电池循环寿命对应的序列结合机器学习模型进行测试得到动力电池寿命测试结果。The sequence corresponding to the battery cycle life with the maximum correlation coefficient higher than the threshold is tested with the machine learning model to obtain the power battery life test result.

本发明的第二个方面提供一种动力电池寿命快速测试系统,包括:A second aspect of the present invention provides a power battery life rapid testing system, comprising:

数据获取模块,用于获取待测试动力电池的变量曲线;A data acquisition module, used to acquire a variable curve of the power battery to be tested;

序列构建模块,用于基于获取的变量曲线得到横向和纵向变量序列;A sequence building module is used to obtain horizontal and vertical variable sequences based on the obtained variable curves;

数据拟合模块,用于将横向和纵向变量序列作为自变量,基于核典型相关分析方法将自变量拟合到电池循环寿命,引入递归核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数;A data fitting module is used to take the horizontal and vertical variable sequences as independent variables, fit the independent variables to the battery cycle life based on the kernel canonical correlation analysis method, and introduce the recursive kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variable and the battery cycle life;

电池寿命测试模块,将最大相关系数高于阈值的电池循环寿命对应的序列结合机器学习模型进行测试得到动力电池寿命测试结果。The battery life test module combines the sequence corresponding to the battery cycle life with the maximum correlation coefficient higher than the threshold with the machine learning model to obtain the power battery life test result.

本发明的第三个方面提供一种计算机可读存储介质。A third aspect of the present invention provides a computer-readable storage medium.

一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如上述所述的一种动力电池寿命快速测试方法中的步骤。A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps in the method for quickly testing the life of a power battery as described above.

本发明的第四个方面提供一种计算机设备。A fourth aspect of the present invention provides a computer device.

一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述所述的一种动力电池寿命快速测试方法中的步骤。A computer device comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps in the method for rapid testing the life of a power battery as described above are implemented.

与现有技术相比,本发明的有益效果是:Compared with the prior art, the present invention has the following beneficial effects:

(1)与传统测试方法相比,本发明采用了以估计代替测试的方法,在电池循环寿命测试中表现出良好的有效性,极大缩短了测试时间,依靠精准预测即可完成寿命测试,极大减少寿命测试时间,加速电池升级换代。(1) Compared with traditional testing methods, the present invention adopts a method of replacing testing with estimation, which shows good effectiveness in battery cycle life testing and greatly shortens the testing time. The life test can be completed by relying on accurate prediction, which greatly reduces the life test time and accelerates battery upgrading.

(2)本发明通过比较不同电池在同一循环中的状态来测试循环寿命和不同电池在同一点的变化趋势来预测循环寿命,引入核典型相关分析方法,通过将变量序列拟合到几乎完美的循环寿命,解决了过度拟合的问题,且可与传统特征相结合,与传统测试方法相比,在电池循环寿命测试中表现出良好的有效性并提高了鲁棒性。(2) The present invention tests the cycle life by comparing the states of different batteries in the same cycle and predicts the cycle life by the changing trends of different batteries at the same point. The kernel canonical correlation analysis method is introduced to solve the problem of overfitting by fitting the variable sequence to an almost perfect cycle life. The method can be combined with traditional features and shows good effectiveness and improved robustness in battery cycle life testing compared with traditional testing methods.

(3)本发明基于核典型相关分析方法将自变量拟合到电池循环寿命,引入核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数,并提出递归典型相关分析求解算法以提高求解运算速度;将最大相关系数高于阈值的电池循环寿命对应的序列结合机器学习模型进行测试得到动力电池寿命测试结果,避免了对先验知识的要求。(3) The present invention fits the independent variable to the battery cycle life based on the kernel canonical correlation analysis method, introduces the kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variable and the battery cycle life, and proposes a recursive canonical correlation analysis solution algorithm to improve the solution operation speed; the sequence corresponding to the battery cycle life with a maximum correlation coefficient higher than the threshold is combined with the machine learning model to obtain the power battery life test result, avoiding the requirement for prior knowledge.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

构成本发明的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。The accompanying drawings in the specification, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

图1是本发明实施例动力电池寿命快速测试方法的整体流程示意图;FIG1 is a schematic diagram of the overall process of a method for rapidly testing the life of a power battery according to an embodiment of the present invention;

图2是本发明实施例电池横向与纵向变量序列示意图。FIG. 2 is a schematic diagram of the transverse and longitudinal variable sequences of a battery according to an embodiment of the present invention.

图3是本发明实施例利用前50循环数据获得的测试结果。FIG. 3 is a test result obtained using the first 50 cycle data according to an embodiment of the present invention.

图4是本发明实施例利用前100循环数据获得的测试结果。FIG. 4 is a test result obtained using the first 100 cycle data according to an embodiment of the present invention.

具体实施方式DETAILED DESCRIPTION

下面结合附图与实施例对本发明作进一步说明。The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

应该指出,以下详细说明都是例示性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。It should be noted that the following detailed descriptions are all illustrative and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.

需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprising" and/or "including" are used in this specification, it indicates the presence of features, steps, operations, devices, components and/or combinations thereof.

本发明的整体思路为:首先,从每个选定的电池变量中导出一个横向序列集和一个纵向序列集,通过核典型相关分析将这些序列拟合到电池循环寿命。利用递归核典型相关分析,减少计算负载;另外,提出递归典型相关分析求解算法以提高求解运算速度。通过上述分析方法方式,生成适当的拟合结果作为机器学习模型的可选输入。最后指定相关系数的阀值,选择高质量的机器学习算法输入。The overall idea of the present invention is as follows: First, a horizontal sequence set and a vertical sequence set are derived from each selected battery variable, and these sequences are fitted to the battery cycle life through kernel canonical correlation analysis. Recursive kernel canonical correlation analysis is used to reduce the computational load; in addition, a recursive canonical correlation analysis solution algorithm is proposed to improve the solution operation speed. Through the above analysis method, appropriate fitting results are generated as optional inputs to the machine learning model. Finally, a threshold of the correlation coefficient is specified to select high-quality machine learning algorithm inputs.

实施例一Embodiment 1

如图1所示,本实施例提供一种动力电池寿命快速测试方法,采用估计代替待测试的方法,包括如下步骤:As shown in FIG1 , this embodiment provides a method for quickly testing the life of a power battery, which uses estimation instead of a method to be tested, and includes the following steps:

S1:获取待测试动力电池的变量曲线;S1: Obtaining a variable curve of the power battery to be tested;

作为一种或多种实施例,S1中,本实施例选取麻省理工学院(MIT)收集的数据集。As one or more embodiments, in S1, this embodiment selects a data set collected by the Massachusetts Institute of Technology (MIT).

需要说明的是,本实施例中选取的数据集中电池的数量和电池的变量曲线可以根据本领域技术人员进行设置。It should be noted that the number of batteries and the variable curves of the batteries in the data set selected in this embodiment can be set according to those skilled in the art.

例如,本实施例选取共123节电池,7个电池变量曲线;For example, in this embodiment, a total of 123 batteries and 7 battery variable curves are selected;

待测试动力电池的变量曲线包括:放电容量-电压曲线(Q(V))、放电温度-电压曲线(T(V))、电压增量-容量(IC)、放电电压-时间曲线(V(t))、放电温度-时间曲线(T(t)))、放电电压/电流-时间曲线(V/I(t))和整个CV(constant voltage的缩写,即充电的恒压阶段)步骤期间的充电电流插值曲线(ICC)。The variable curves of the power battery to be tested include: discharge capacity-voltage curve (Q(V)), discharge temperature-voltage curve (T(V)), voltage increment-capacity (IC), discharge voltage-time curve (V(t)), discharge temperature-time curve (T(t)), discharge voltage/current-time curve (V/I(t)) and the charging current interpolation curve (ICC) during the entire CV (abbreviation for constant voltage, i.e. the constant voltage stage of charging) step.

其中,ICC采用线性插值进行标准化。Among them, ICC is standardized using linear interpolation.

S2:基于获取的变量曲线得到横向和纵向变量序列;S2: Obtain horizontal and vertical variable sequences based on the obtained variable curves;

如图2所示,其中,将变量曲线中一个周期内变化的电池变量构成的序列作为横向变量序列,不同周期内同一采样点变化的电池变量构成的序列作为纵向变量序列。As shown in FIG2 , the sequence of battery variables that change within one cycle in the variable curve is taken as the horizontal variable sequence, and the sequence of battery variables that change at the same sampling point in different cycles is taken as the vertical variable sequence.

不同横向变量序列曲线在后期的循环中更加明显,所以通常表现更好。The curves of different horizontal variable series are more obvious in the later cycles, so they usually perform better.

相反,在早期的循环中,电池之间的差异很小且不规则,纵向序列的性能与采样点的位置相关,通常在变量变化更大、更规律的情况下表现更好。In contrast, in early cycles, where differences between cells are small and irregular, the performance of longitudinal series is related to the location of the sampling points, generally performing better where the variables vary more widely and more regularly.

例如,基于S1获取的数据的基础上,每个变量派生出一个横向和纵向变量序列集,即从每个选择的变量中导出横向和纵向变量序列集。For example, based on the data obtained by S1, a set of horizontal and vertical variable sequences is derived for each variable, that is, a set of horizontal and vertical variable sequences is derived from each selected variable.

按照此方式,可生成14种变量序列作为后续采用递归的核典型相关分析方法输入的自变量。In this way, 14 variable sequences can be generated as independent variables for subsequent input of the recursive kernel canonical correlation analysis method.

S3:将横向和纵向变量序列作为自变量,基于核典型相关分析方法将自变量拟合到电池循环寿命,引入递归核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数;S3: Taking the horizontal and vertical variable sequences as independent variables, fitting the independent variables to the battery cycle life based on the kernel canonical correlation analysis method, and introducing the recursive kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variables and the battery cycle life;

核典型相关分析(Kernel Canonical Correlation Analysis,KCCA)可以通过横向变量序列比较不同电池在同一循环中的状态来测试电池的循环寿命,通过纵向变量序列比较不同电池在同一点的变化趋势来测试电池的循环寿命。Kernel Canonical Correlation Analysis (KCCA) can test the cycle life of batteries by comparing the states of different batteries in the same cycle through a horizontal variable sequence, and can test the cycle life of batteries by comparing the change trends of different batteries at the same point through a vertical variable sequence.

由于KCCA可以在训练步骤中将变量序列拟合到几乎完美的循环寿命,因此必须解决过度拟合的问题,KCCA在训练组上表现良好的超参数可能不适合测试组。Since KCCA can fit the variable sequence to almost perfect cycle life in the training step, the problem of overfitting must be addressed. The hyperparameters for which KCCA performs well on the training set may not be suitable for the test set.

例如,如果在训练步骤中简单地选择小的η以使相关系数变大,则拟合结果在测试组可能不满意。For example, if a small η is simply chosen in the training step to make the correlation coefficient larger, the fitting results may not be satisfactory in the test set.

为了解决这个问题,本实施例创建了一个KCCA验证组来测试和选择KCCA超参数。考虑到批次之间日历老化和休息时间之间的差异,将每批次中大致相同数量的电池分配给训练、验证和测试组。To address this issue, this example creates a KCCA validation group to test and select KCCA hyperparameters. Considering the differences in calendar aging and rest time between batches, approximately the same number of cells in each batch are allocated to the training, validation, and test groups.

如图2所示,将数据集分为三组,KCCA训练组、KCCA验证组和KCCA测试组,每组41个电池,KCCA训练组和KCCA验证组都构成了机器学习训练组。As shown in Figure 2, the dataset is divided into three groups, KCCA training group, KCCA validation group and KCCA test group, each with 41 batteries. The KCCA training group and the KCCA validation group both constitute the machine learning training group.

KCCA训练组使用不同的超参数集对变量序列进行训练,不同的超参数集通过网格搜索方法有效识别。The KCCA training group trains the variable sequence using different hyperparameter sets, which are effectively identified by the grid search method.

本实施例基于核典型相关分析(KCCA),提出了递归KCCA(RKCCA)求解算法,以方便计算。This embodiment proposes a recursive KCCA (RKCCA) solution algorithm based on kernel canonical correlation analysis (KCCA) to facilitate calculation.

传统的典型相关分析(CCA)旨在寻找两组多维变量之间的线性关系。Traditional canonical correlation analysis (CCA) aims to find the linear relationship between two sets of multidimensional variables.

两组随机变量,形式为(x,y),均值为零,假设有一个动力电池数据实例S=((x1,y1),…,(xn,yn))的样本(x,y),其中,x表示由动力电池横向和纵向变量序列组成的序列,即自变量,y表示对应的动力电池的寿命,即因变量。Two sets of random variables, of the form (x, y), The mean is zero. Suppose there is a sample (x, y) of a power battery data instance S = ((x 1 , y 1 ), …, (x n , yn )), where x represents a sequence composed of horizontal and vertical variable sequences of the power battery, that is, the independent variable, and y represents the life of the corresponding power battery, that is, the dependent variable.

其中,用X表示为(x1,…,xn)T,Y表示为(y1,…,yn)THere, X is represented by (x 1 , …, x n ) T , and Y is represented by (y 1 , …, yn ) T .

假设保证n≥rank(X),rank(Y),且X和Y都具有零均值。定义一个新的向量w,命名为x的加载向量并将x投影到该方向,通过选择加载向量v对y来做同样的事情,获得y的投影。Assume that n ≥ rank(X), rank(Y), and that both X and Y have zero mean. Define a new vector w, named the loading vector of x and project x into that direction, and do the same for y by choosing a loading vector v to obtain the projection of y.

将投影t=Xw和u=Yv分别定义为X和Y的得分向量。CCA确定w和v以最大化相关系数在t=Xw和u=Yv之间。Define the projections t=Xw and u=Yv as the score vectors of X and Y, respectively. CCA determines w and v to maximize the correlation coefficient between t=Xw and u=Yv.

在本实施例中,所有||·||表示欧几里得范数。In this embodiment, all ||·|| represent Euclidean norms.

所述将横向和纵向变量序列作为自变量,基于核典型相关分析方法将自变量拟合到电池循环寿命,具体包括:The method of using the horizontal and vertical variable sequences as independent variables and fitting the independent variables to the battery cycle life based on the kernel canonical correlation analysis method specifically includes:

S301:采用核函数化解将自变量映射至高维空间包括:S301: Using kernel function to map independent variables to high-dimensional space includes:

将每个xi映射到一个具有无限维的高维向量空间F:Map each xi into a high-dimensional vector space F with infinite dimensions:

其中,xi表示自变量序列中的每个元素,Φ(xi)表示将自变量序列中的每个元素映射至高维空间的结果;Wherein, xi represents each element in the independent variable sequence, and Φ( xi ) represents the result of mapping each element in the independent variable sequence to a high-dimensional space;

Φ(X)=(Φ(x1),Φ(x2),…,Φ(xn))T (2)Φ(X)=(Φ(x 1 ), Φ(x 2 ),…, Φ(x n )) T (2)

其中,Φ(X)的均值也应该为零:Among them, the mean of Φ(X) should also be zero:

通过应用核技巧来避免求解具体Φ(xi):Avoid solving the specific Φ( xi ) by applying the kernel trick:

K(X)i,j=Φ(xi)Φ(xj)T=k(xi,xj) (4)K(X) i,j =Φ(x i )Φ(x j ) T =k(x i ,x j ) (4)

K(X)=Φ(X)Φ(X)T (5)K(X)=Φ(X)Φ(X) T (5)

其中,核函数k(xi,xj)是一个对称函数。因此,K(X)是一个对称矩阵,表示自变量核矩阵。Among them, the kernel function k( xi , xj ) is a symmetric function. Therefore, K(X) is a symmetric matrix, which represents the independent variable kernel matrix.

S302:中心化S302: Centralization

通常k(xi,xj)不能确保Φ(X)的均值为零,因此需要中心化:Usually k( xi , xj ) cannot ensure that the mean of Φ(X) is zero, so centering is required:

S303:计算映射至高维空间自变量对应的序列和电池循环寿命的投影;以上述投影之间的相关系数最大化为目标构建优化函数。S303: Calculate the projection of the sequence and battery cycle life corresponding to the independent variable mapped to the high-dimensional space; and construct an optimization function with the goal of maximizing the correlation coefficient between the above projections.

需要说明的是,为了方便,本实施例的其余部分,分别使用Φ和K代替Φ(X)和K(X)。It should be noted that, for convenience, in the rest of this embodiment, Φ and K are used to replace Φ(X) and K(X) respectively.

优化问题变成最大化相关系数ρ1The optimization problem becomes maximizing the correlation coefficient ρ 1 :

其中,Φ和Φ(X)含义一致,表示映射到高维空间的自变量矩阵,Y表示因变量矩阵,w1表示映射到高维空间的自变量矩阵的投影向量,v1表示因变量矩阵的投影向量。Among them, Φ and Φ(X) have the same meaning, representing the independent variable matrix mapped to the high-dimensional space, Y represents the dependent variable matrix, w 1 represents the projection vector of the independent variable matrix mapped to the high-dimensional space, and v 1 represents the projection vector of the dependent variable matrix.

Φ通常满秩,所以wi始终是Φ(xi)的线性组合:Φ is usually full rank, so wi is always a linear combination of Φ( xi ):

wi=ΦTαi (8)w iT α i (8)

与CCA类似,可以得到:Similar to CCA, we can get:

(KTK)-1KTY(YTY)-1YT1=ρ1 2α1 (9)(K T K) -1 K T Y(Y T Y) -1 Y T11 2 α 1 (9)

(YTY)-1YTK(KTK)-1KTYv1=ρ1 2v1 (10)(Y T Y) -1 Y T K(K T K) -1 K T Yv 11 2 v 1 (10)

但是,K的对称的,这导致:However, K is symmetric, which leads to:

(KTK)-1KTY(YTY)-1YTK=I (11)(K T K) -1 K T Y(Y T Y) -1 Y T K=I (11)

其中,I是单位矩阵,所有的特征值都等于1。在任意选择αi时|ρi||=1总成立,总可以形成完美的相关性。Where I is the identity matrix, and all eigenvalues are equal to 1. When α i is arbitrarily selected, |ρ i ||=1 always holds, and a perfect correlation can always be formed.

核方法中需要正则化,ηα1 T1项被添加到优化方程:Regularization is required in kernel methods, and the term ηα 1 T1 is added to the optimization equation:

其中,η是手动选择的超参数,通常η很小。αTKα是核偏最小二乘(KPLS)方法的优化问题中的项。where η is a manually chosen hyperparameter, usually small, and α T Kα is a term in the optimization problem of the kernel partial least squares (KPLS) method.

正则化问题可以看作是非正则化KCCA和KPLS的过渡形式。The regularization problem can be viewed as a transitional form of non-regularized KCCA and KPLS.

在这一步之后,可以确保方程是可解的且矩阵特征值不总为1:After this step, it is guaranteed that the equation is solvable and that the matrix eigenvalues are not always 1:

(KTK+η1I)-1KTY(YTY)-1YT1=λ1 2α1: (13)(K T K+η 1 I) -1 K T Y(Y T Y) -1 Y T11 2 α 1 : (13)

(YTY)-1YTK(KTK+η1I)-1KTYv1=λ1 2v1 (14)(Y T Y) -1 Y T K(K T K+η 1 I) -1 K T Yv 11 2 v 1 (14)

需要强调的是,由于限制已更改,λ1 2与ρ1 2并不完全相同,但它们的差异应该很小。It is important to emphasize that λ 1 2 is not exactly the same as ρ 1 2 because the restrictions have changed, but their differences should be small.

其中,所述引入递归核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数包括:The maximum correlation coefficient between the independent variable and the battery cycle life obtained by introducing the recursive kernel canonical correlation analysis method includes:

基于乘幂法,基于给定的可对角化矩阵和随机非零向量,采用迭代求解方法找到可对角化矩阵最大特征值,根据特征值和相关系数之间的关系得到最大特征值对应的最大相关系数。Based on the power method, based on a given diagonalizable matrix and a random non-zero vector, an iterative solution method is used to find the maximum eigenvalue of the diagonalizable matrix, and the maximum correlation coefficient corresponding to the maximum eigenvalue is obtained according to the relationship between the eigenvalue and the correlation coefficient.

与CCA类似,α和v是(KTK+ηI)-1KTY(YTY)-1YT1和(YTY)-1YTK(KTK+η1I)-1KTYv1的对应特征向量的解。Similar to CCA, α and v are the corresponding eigenvectors of the solution (K T K+ηI) -1 K T Y(Y T Y) -1 Y T1 and (Y T Y) -1 Y T K(K T K+η 1 I) -1 K T Yv 1 .

因为只需要找到最大特征值,使用乘幂法来完成这项任务是一种更好的方法。Since we only need to find the largest eigenvalue, using the power method is a better approach to accomplish this task.

给定一个可对角化矩阵D和一个随机非零向量b0,当k→∞时,D的最大特征值对应的特征向量bk可以通过以下递推关系得到:Given a diagonalizable matrix D and a random non-zero vector b 0 , when k→∞, the eigenvector b k corresponding to the maximum eigenvalue of D can be obtained by the following recursive relation:

可得α和v之间的关系:The relationship between α and v can be obtained:

v∝(YYT)-1YKTα (16)v∝(YY T ) -1 YK T α (16)

α∝(KKT+ηI)-1KYTv (17)α∝(KK T +ηI) -1 KY T v (17)

本实施例中,根据典型相关分析CCA的性质,所述特征值和相关系数之间的关系为:特征值的平方为对应的相关系数。In this embodiment, according to the properties of canonical correlation analysis CCA, the relationship between the eigenvalue and the correlation coefficient is: the square of the eigenvalue is the corresponding correlation coefficient.

针对这一特点,本实施例基于给定的可对角化矩阵和随机非零向量,采用迭代求解方法找到可对角化矩阵最大特征值,包括:In view of this feature, this embodiment uses an iterative solution method to find the maximum eigenvalue of the diagonalizable matrix based on a given diagonalizable matrix and a random non-zero vector, including:

(1)选择一个随机的非零向量作为α的初始值,通常是K的第一列(对核矩阵的权重向量alpha标准化);(1) Select a random non-zero vector as the initial value of α, usually the first column of K (alpha normalize the weight vector of the kernel matrix);

然后执行以下迭代:Then perform the following iterations:

(2)Φ的得分向量:t=Kα,(收敛时映射到高维空间的自变量矩阵的得分向量t为自变量核矩阵K与权重向量alpha的积,对t标准化)(2) The score vector of Φ: t = Kα, (The score vector t of the independent variable matrix mapped to the high-dimensional space at convergence is the product of the independent variable kernel matrix K and the weight vector alpha, and t is standardized)

(3)Y载荷向量:v=(YTY)-1YTt,(收敛时因变量的载荷向量v与因变量矩阵Y、映射到高维空间的自变量矩阵的得分向量t的关系,对v标准化)(3) Y load vector: v = (Y T Y) -1 Y T t, (The relationship between the load vector v of the dependent variable and the dependent variable matrix Y and the score vector t of the independent variable matrix mapped to the high-dimensional space at convergence, and v is standardized)

(4)Y的得分向量:u=Yv;(收敛时因变量的得分向量u为因变量矩阵Y与因变量的载荷向量v的积)(4) Score vector of Y: u = Yv; (When converged, the score vector u of the dependent variable is the product of the dependent variable matrix Y and the dependent variable loading vector v)

(5)权重向量:α=(KKT+ηI)-1KYTv,(收敛时核矩阵的权重向量alpha与核矩阵、因变量矩阵Y、因变量的载荷向量v的关系,对核矩阵的权重向量alpha标准化)(5) Weight vector: α = (KK T + ηI) -1 KY T v, (The relationship between the weight vector alpha of the kernel matrix and the kernel matrix, the dependent variable matrix Y, and the load vector v of the dependent variable at convergence, and the weight vector alpha of the kernel matrix is standardized)

通过上述方式,本实施例可以以微小的精度损失为代价,以更方便的方式找到最大特征值和对应的特征向量,而不是对两个矩阵进行完整的特征向量分解。In the above manner, this embodiment can find the maximum eigenvalue and the corresponding eigenvector in a more convenient way at the cost of a slight loss of accuracy, rather than performing a complete eigenvector decomposition on the two matrices.

在KCCA的测试算法中,x被认为是独立的变量,y作为因变量。目标是根据先前捕获的关系用新的给定自变量Xnew测试Ynew。Xnew的核矩阵Knew预处理如下:In the KCCA testing algorithm, x is considered as an independent variable and y as a dependent variable. The goal is to test Y new with a new given independent variable X new based on the previously captured relationship. The kernel matrix K new of X new is preprocessed as follows:

Knew=Φ(Xnew)Φ(X)T (18)K new =Φ(X new )Φ(X) T (18)

然后中心化KnewThen centralize K new :

通过假设Xnew和Ynew的得分矩阵相同来进行测试Ynew。Xnew的得分矩阵为:The test Y new is performed by assuming that the score matrices of X new and Y new are the same. The score matrix of X new is:

Tnew=KnewA (20)T new = K new A (20)

将Tnew视为Ynew的得分矩阵:Think of T new as the score matrix of Y new :

可得估计的YnewThe estimated Y new can be obtained:

S4:将最大相关系数高于阈值的电池循环寿命对应的序列结合机器学习模型进行测试得到动力电池寿命测试结果。S4: The sequence corresponding to the battery cycle life with a maximum correlation coefficient higher than the threshold is tested in combination with the machine learning model to obtain the power battery life test result.

本实施例考虑到测试模型输入的数据的质量比数量更重要,同时深度学习方法的性能优于传统机器学习方法(GPR和Elastic net)。This embodiment takes into account that the quality of the data input to the test model is more important than the quantity, and the performance of the deep learning method is better than the traditional machine learning method (GPR and Elastic net).

本实施例,使用4个具有代表性的机器学习模型测试了这些相关系数高于阈值的拟合结果对应的序列,所述机器学习模型包括人工神经网络(ANN)、随机森林(RF)、高斯过程回归(GPR)和弹性网络(Elastic net)。In this embodiment, four representative machine learning models are used to test the sequences corresponding to the fitting results with correlation coefficients higher than the threshold, and the machine learning models include artificial neural network (ANN), random forest (RF), Gaussian process regression (GPR) and elastic net.

利用前100个周期的数据,ANN、RF、GPR和Elastic net的最佳平均绝对百分比误差(MAPE)分别为8.4%、8.6%、11.1%和14.0%,使用前50个周期的数据,最佳结果范围为9.5%到14.9%,其中RF实现了最佳性能。通过上述实验验证了ANN和RF的性能优于传统机器学习算法(GPR和Elastic网络)。Using the data from the first 100 cycles, the best mean absolute percentage errors (MAPE) for ANN, RF, GPR, and Elastic net were 8.4%, 8.6%, 11.1%, and 14.0%, respectively. Using the data from the first 50 cycles, the best results ranged from 9.5% to 14.9%, with RF achieving the best performance. The above experiments verified that the performance of ANN and RF is better than that of traditional machine learning algorithms (GPR and Elastic network).

对于阈值的设定,并非是越高的相关系数阀值总是更好,尤其是对于深度学习方法,较低的阈值意味着更多的拟合结果可以输入到机器学习模型中,这也有助于测试。Regarding threshold setting, a higher correlation coefficient threshold is not always better, especially for deep learning methods. A lower threshold means that more fitting results can be input into the machine learning model, which is also helpful for testing.

因此阈值的设定根据实际选取的机器学习模型进行设置。Therefore, the threshold is set according to the actual selected machine learning model.

例如,通常对于ANN,随着阈值的降低,MAPE先降低后升高,当阈值设置为0.8时出现最低点,在这些输入的数量和质量之间存在平衡。对于GPR和Elastic net,它们的最佳性能出现在阈值高时。For example, generally for ANN, as the threshold decreases, MAPE decreases first and then increases, with the lowest point occurring when the threshold is set to 0.8, where there is a balance between the quantity and quality of these inputs. For GPR and Elastic net, their best performance occurs when the threshold is high.

为了验证本实施例方法的有效性,进行了实验:In order to verify the effectiveness of the method in this embodiment, the following experiments were conducted:

表1中为基于MIT数据集的实验结果,图3是本发明实施例利用前50循环数据获得的测试结果;图4是本发明实施例利用前100循环数据获得的测试结果。Table 1 shows the experimental results based on the MIT data set, FIG. 3 shows the test results obtained by the embodiment of the present invention using the first 50 cycle data, and FIG. 4 shows the test results obtained by the embodiment of the present invention using the first 100 cycle data.

通过引入MAPE和RMSE进行评价实验结果,其中MAPE为平均绝对百分比误差,RMSE为均方根误差。The experimental results are evaluated by introducing MAPE and RMSE, where MAPE is the mean absolute percentage error and RMSE is the root mean square error.

最终实验结果表明,对比传统全寿命周期测试需2000次,本发明以少量精度为代价,依靠精准测试即可完成寿命测试,仅需50至300次循环即可完成估计,无需全面测试,缩短1950个测试循环,速度提高40倍,可极大减少寿命测试时间,加速电池升级换代。The final experimental results show that compared with the traditional full life cycle test which requires 2,000 times, the present invention can complete the life test at the expense of a small amount of accuracy by relying on precise testing. It only takes 50 to 300 cycles to complete the estimation, without the need for comprehensive testing, shortening the test cycles by 1,950 and increasing the speed by 40 times, which can greatly reduce the life test time and accelerate battery upgrades.

表1实验结果Table 1 Experimental results

循环数Number of cycles MAPE(%)MAPE(%) RMSE(循环)RMSE(loop) 缩短循环数Shorten the number of cycles 50(2.5%)50(2.5%) 9.19.1 122122 19501950 100(5%)100(5%) 8.28.2 9898 19001900 200(10%)200(10%) 8.08.0 9393 18001800 300(15%)300(15%) 7.77.7 8989 17001700

实施例二Embodiment 2

本实施例提供一种动力电池寿命快速测试系统,包括:This embodiment provides a power battery life rapid testing system, including:

数据获取模块,用于获取待测试动力电池的变量曲线;A data acquisition module, used to acquire a variable curve of the power battery to be tested;

序列构建模块,用于基于获取的变量曲线得到横向和纵向变量序列;A sequence building module is used to obtain horizontal and vertical variable sequences based on the obtained variable curves;

数据拟合模块,用于将横向和纵向变量序列作为自变量,基于核典型相关分析方法将自变量拟合到电池循环寿命,引入递归核典型相关分析方法得到自变量和电池循环寿命之间的最大相关系数;A data fitting module is used to take the horizontal and vertical variable sequences as independent variables, fit the independent variables to the battery cycle life based on the kernel canonical correlation analysis method, and introduce the recursive kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variable and the battery cycle life;

电池寿命测试模块,将最大相关系数高于阈值的电池循环寿命对应的序列结合机器学习模型进行测试得到动力电池寿命测试结果。The battery life test module combines the sequence corresponding to the battery cycle life with the maximum correlation coefficient higher than the threshold with the machine learning model to obtain the power battery life test result.

实施例三Embodiment 3

本实施例提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如上述所述的一种动力电池寿命快速测试方法中的步骤。This embodiment provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the steps in the method for quickly testing the life of a power battery as described above are implemented.

实施例四Embodiment 4

本实施例提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述所述的一种动力电池寿命快速测试方法中的步骤。This embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the method for rapid testing the life of a power battery as described above are implemented.

本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of hardware embodiments, software embodiments, or embodiments combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage and optical storage, etc.) containing computer-usable program codes.

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention is described with reference to the flowchart and/or block diagram of the method, device (system), and computer program product according to the embodiment of the present invention. It should be understood that each process and/or box in the flowchart and/or block diagram, as well as the combination of the process and/or box in the flowchart and/or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(RandomAccessMemory,RAM)等。A person skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program, and the program can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes of the embodiments of the above-mentioned methods. The storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM), etc.

以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims (9)

1. The rapid power battery life test method is characterized by adopting an estimated replacement test method and comprises the following steps:
Acquiring a variable curve of a power battery to be tested;
the variable curve of the power battery to be tested comprises a discharge capacity-voltage curve, a discharge temperature-voltage curve, a voltage increment-capacity curve, a discharge voltage-time curve, a discharge temperature-time curve, a discharge voltage/current-time curve and a charge current interpolation curve during the constant voltage stage step of the whole charge;
obtaining a transverse variable sequence and a longitudinal variable sequence based on the obtained variable curve;
The transverse and longitudinal variable sequences obtained based on the obtained variable curves are as follows:
taking a sequence formed by battery variables changing in one period in a variable curve as a transverse variable sequence, and taking a sequence formed by battery variables changing in the same sampling point in different periods as a longitudinal variable sequence;
Using the transverse and longitudinal variable sequences as independent variables, fitting the independent variables to the cycle life of the battery based on a kernel canonical correlation analysis method, and introducing a recursive kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variables and the cycle life of the battery;
And combining the sequence corresponding to the battery cycle life with the maximum correlation coefficient higher than the threshold value with a machine learning model to test so as to obtain a power battery life test result.
2. The rapid power cell life test method of claim 1, wherein said fitting the independent variables to the battery cycle life based on a kernel-based canonical correlation analysis method using the sequence of lateral and longitudinal variables as the independent variables comprises:
Mapping the sequence corresponding to the independent variable to a high-dimensional space by adopting a kernel function solution;
Calculating projections of the sequence and the battery cycle life corresponding to the high-dimensional space independent variable;
And constructing an optimization function by maximizing the correlation coefficient between the projections as a target.
3. A rapid power battery life test method according to claim 2, wherein the term in the optimization problem of the hyper-parametric and the kernel-biased least squares method is added to the objective function by regularization when the optimization function is constructed.
4. The rapid power battery life test method of claim 1, wherein said introducing a recursive kernel canonical correlation analysis method to obtain a maximum correlation coefficient between the independent variable and the battery cycle life comprises: based on the exponentiation method, based on a given diagonable matrix and a random non-zero vector, the maximum eigenvalue of the diagonable matrix is found by adopting an iterative solution method, and the maximum correlation coefficient corresponding to the maximum eigenvalue is obtained according to the relation between the eigenvalue and the correlation coefficient.
5. A rapid power battery life test system, characterized in that it is based on the rapid power battery life test method according to any one of claims 1 to 4, further comprising:
The data acquisition module is used for acquiring a variable curve of the power battery to be tested;
The sequence construction module is used for obtaining transverse and longitudinal variable sequences based on the obtained variable curves;
The transverse and longitudinal variable sequences obtained based on the obtained variable curves are as follows:
taking a sequence formed by battery variables changing in one period in a variable curve as a transverse variable sequence, and taking a sequence formed by battery variables changing in the same sampling point in different periods as a longitudinal variable sequence;
the data fitting module is used for taking the transverse and longitudinal variable sequences as independent variables, fitting the independent variables to the cycle life of the battery based on a kernel canonical correlation analysis method, and introducing a recursive kernel canonical correlation analysis method to obtain the maximum correlation coefficient between the independent variables and the cycle life of the battery;
And the battery life test module is used for testing the sequence corresponding to the battery cycle life with the maximum correlation coefficient higher than the threshold value by combining the machine learning model to obtain a power battery life test result.
6. The rapid power cell life test system of claim 5 wherein said fitting the independent variables to the battery cycle life based on a kernel-based canonical correlation analysis method using the transverse and longitudinal variable sequences as independent variables comprises:
Mapping the sequence corresponding to the independent variable to a high-dimensional space by adopting a kernel function solution;
Calculating projections of the sequence and the battery cycle life corresponding to the high-dimensional space independent variable;
And constructing an optimization function by maximizing the correlation coefficient between the projections as a target.
7. The rapid power battery life test system of claim 5 wherein said introducing a recursive kernel canonical correlation analysis method results in a maximum correlation coefficient between the independent variable and the battery cycle life comprising: based on the exponentiation method, based on a given diagonable matrix and a random non-zero vector, the maximum eigenvalue of the diagonable matrix is found by adopting an iterative solution method, and the maximum correlation coefficient corresponding to the maximum eigenvalue is obtained according to the relation between the eigenvalue and the correlation coefficient.
8. A computer readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the steps of a power battery life fast test method according to any one of claims 1-4.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor performs the steps of a method for rapid testing of the life of a power battery as claimed in any one of claims 1 to 4.
CN202210789672.XA 2022-07-06 2022-07-06 A method and system for quickly testing the life of a power battery Active CN115219902B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210789672.XA CN115219902B (en) 2022-07-06 2022-07-06 A method and system for quickly testing the life of a power battery

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210789672.XA CN115219902B (en) 2022-07-06 2022-07-06 A method and system for quickly testing the life of a power battery

Publications (2)

Publication Number Publication Date
CN115219902A CN115219902A (en) 2022-10-21
CN115219902B true CN115219902B (en) 2024-08-16

Family

ID=83610385

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210789672.XA Active CN115219902B (en) 2022-07-06 2022-07-06 A method and system for quickly testing the life of a power battery

Country Status (1)

Country Link
CN (1) CN115219902B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115825755B (en) * 2022-12-30 2023-09-19 哈尔滨昆宇新能源有限公司 Method for evaluating consistency of voltages of battery cells of energy storage battery

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102778653A (en) * 2012-06-20 2012-11-14 哈尔滨工业大学 Data-driven lithium ion battery cycle life prediction method based on AR (Autoregressive) model and RPF (Regularized Particle Filtering) algorithm
CN103954913A (en) * 2014-05-05 2014-07-30 哈尔滨工业大学深圳研究生院 Predication method of electric vehicle power battery service life

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7591818B2 (en) * 2001-12-04 2009-09-22 Endoscopic Technologies, Inc. Cardiac ablation devices and methods
US7411400B2 (en) * 2003-12-19 2008-08-12 Battery Control Corp. Method for testing battery condition
US10393813B2 (en) * 2013-08-27 2019-08-27 The Regents Of The University Of Michigan On-board state of health monitoring of batteries using incremental capacity analysis
AT519335B1 (en) * 2016-11-04 2022-08-15 Avl List Gmbh Method for diagnosing a technical system
CN107024663B (en) * 2017-04-01 2019-04-05 湖南银杏数据科技有限公司 Lithium battery screening technique based on charging curve feature KPCA cluster
CN107728072A (en) * 2017-10-10 2018-02-23 合肥国轩高科动力能源有限公司 Method for rapidly predicting cycle life of lithium ion battery
CN107942261B (en) * 2017-12-29 2021-04-02 上海电气集团股份有限公司 Method and system for estimating state of charge of battery
KR20220011487A (en) * 2020-07-21 2022-01-28 삼성전자주식회사 Method for identifying status of battery, the electronic device and storage medium therefor
GB2600757A (en) * 2020-11-09 2022-05-11 Horiba Mira Ltd Battery performance optimisation
CN113533988B (en) * 2021-06-04 2022-09-27 上海空间电源研究所 Long-term circulation capacity attenuation analysis method for lithium ion battery

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102778653A (en) * 2012-06-20 2012-11-14 哈尔滨工业大学 Data-driven lithium ion battery cycle life prediction method based on AR (Autoregressive) model and RPF (Regularized Particle Filtering) algorithm
CN103954913A (en) * 2014-05-05 2014-07-30 哈尔滨工业大学深圳研究生院 Predication method of electric vehicle power battery service life

Also Published As

Publication number Publication date
CN115219902A (en) 2022-10-21

Similar Documents

Publication Publication Date Title
CN113030764B (en) A method and system for estimating the state of health of a battery pack
Lipu et al. Real-time state of charge estimation of lithium-ion batteries using optimized random forest regression algorithm
CN110888059B (en) Charge state estimation algorithm based on improved random forest combined volume Kalman
CN113219343A (en) Lithium battery health state prediction method, system, equipment and medium based on elastic network
CN111007417A (en) Battery pack SOH and RUL prediction method and system based on inconsistency assessment
CN111448467B (en) Method and system for modeling and estimating battery capacity
CN110443377B (en) Support vector machine storage battery life prediction method based on immune algorithm optimization
JP2023501995A (en) Battery performance prediction
CN107132490B (en) A method for realizing state of charge estimation of lithium battery pack
CN109001640B (en) Data processing method and device for power battery
WO2022213789A1 (en) Lithium battery soc estimation method and apparatus, and computer-readable storage medium
Chen et al. State of health estimation of lithium-ion batteries based on equivalent circuit model and data-driven method
CN110095731B (en) A direct prediction method for remaining life of long-life space lithium-ion batteries
JP2021150290A (en) Battery diagnostic system for estimating remaining useful life (rul) of battery
CN113189490B (en) Lithium battery health state estimation method based on feature screening and Gaussian process regression
WO2024016496A1 (en) Method and apparatus for estimating soh of lithium battery
Mazzi et al. Lithium-ion battery state of health estimation using a hybrid model based on a convolutional neural network and bidirectional gated recurrent unit
CN113721149A (en) Lithium battery capacity prediction method based on semi-supervised transfer learning
CN115219902B (en) A method and system for quickly testing the life of a power battery
CN110806540A (en) Cell test data processing method, device, system and storage medium
Wu et al. Remaining useful life prediction of Lithium-ion batteries based on PSO-RF algorithm
TWI465939B (en) A Multi - stage Fast Charge Method for Optimizing Lithium Batteries
CN115291126B (en) Lithium battery SOC estimation method, device and medium based on SGPR and adaptive infinite filtering
Cai et al. Prediction of lithium-ion battery remaining useful life based on hybrid data-driven method with optimized parameter
CN115586452A (en) Lithium ion battery health state estimation method based on novel health characteristics

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant