+

CN111948563A - A method for predicting the remaining life of electric forklift lithium battery based on multi-neural network coupling - Google Patents

A method for predicting the remaining life of electric forklift lithium battery based on multi-neural network coupling Download PDF

Info

Publication number
CN111948563A
CN111948563A CN202010568146.1A CN202010568146A CN111948563A CN 111948563 A CN111948563 A CN 111948563A CN 202010568146 A CN202010568146 A CN 202010568146A CN 111948563 A CN111948563 A CN 111948563A
Authority
CN
China
Prior art keywords
lithium battery
prediction
input
open
time
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
CN202010568146.1A
Other languages
Chinese (zh)
Other versions
CN111948563B (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.)
Zhejiang University ZJU
Hangcha Group Co Ltd
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202010568146.1A priority Critical patent/CN111948563B/en
Publication of CN111948563A publication Critical patent/CN111948563A/en
Application granted granted Critical
Publication of CN111948563B publication Critical patent/CN111948563B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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]
    • G01R31/396Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
    • 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]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • 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]
    • G01R31/392Determining battery ageing or deterioration, e.g. state of health
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/049Temporal neural networks, e.g. delay elements, oscillating neurons or pulsed inputs
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Secondary Cells (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

本发明是为了克服现有技术预测锂电池剩余寿命时,计算复杂耗时长,预测精度低的问题,提供一种基于多神经网络耦合的电动叉车锂电池剩余寿命预测方法,提高了预测计算精度,减少了预测模型训练时间,包括以下步骤:建立基于长短时记忆神经网络的开路电压预测模型,采用RMSprop算法和dropout正则化方法对网络进行优化,从而预测锂电池在放电循环中的开路电压值VOC;把预测结果按顺序划分成多个放电循环,统计每个放电循环中从初始电压至最小电压间的开路电压样本个数NS,利用采样时间TS相同,得到每个放电循环中放电至最小电压的时间Tmin;建立基于人工神经网络的容量预测模型,以预测锂电池容量C,从而得到锂电池剩余寿命预测值RUL。

Figure 202010568146

In order to overcome the problems of complicated calculation and long time-consuming and low prediction accuracy when predicting the remaining life of a lithium battery in the prior art, the invention provides a method for predicting the remaining life of an electric forklift lithium battery based on multi-neural network coupling, which improves the prediction calculation accuracy. The training time of the prediction model is reduced, including the following steps: establishing an open-circuit voltage prediction model based on a long-short-term memory neural network, and using the RMSprop algorithm and the dropout regularization method to optimize the network to predict the open-circuit voltage value V of the lithium battery during the discharge cycle. OC ; Divide the prediction results into multiple discharge cycles in sequence, count the number of open-circuit voltage samples N S from the initial voltage to the minimum voltage in each discharge cycle, and use the same sampling time T S to obtain the discharge in each discharge cycle. Time T min to the minimum voltage; establish a capacity prediction model based on an artificial neural network to predict the lithium battery capacity C, thereby obtaining the predicted value RUL of the remaining life of the lithium battery.

Figure 202010568146

Description

一种基于多神经网络耦合的电动叉车锂电池剩余寿命预测 方法A prediction of remaining life of electric forklift lithium battery based on multi-neural network coupling method

技术领域technical field

本发明涉及电池技术领域,尤其是涉及一种基于多神经网络耦合的电池剩余寿命预测方法。The invention relates to the technical field of batteries, in particular to a method for predicting the remaining life of a battery based on multi-neural network coupling.

背景技术Background technique

锂电池的能量密度高、循环寿命长、放电平台高,可为多种工程机械提供可靠的动力源。目前电动叉车多采用锂电池作为其动力源,锂电池在不断充放电的情况下,其电化学性能会逐渐下降,容量慢慢出现衰减,在达到最大寿命后锂电池性能衰减速度加快,发热量增加,安全性降低,容易引起设备运行不稳定造成经济的损失,严重的甚至会引发安全事故。因此,锂电池需要及时进行维护和更换,而剩余寿命预测方法利用能够提前给出电池容量的衰减情况,从而可以提前制定出合理的维护方案,从而保证电动叉车的稳定运行。剩余寿命预测是指在已知部分容量衰减数据的情况下,去预测电池在达到寿命终止前仍能运行的循环数。Lithium batteries have high energy density, long cycle life and high discharge platform, which can provide a reliable power source for a variety of construction machinery. At present, most electric forklift trucks use lithium batteries as their power source. When the lithium battery is continuously charged and discharged, its electrochemical performance will gradually decrease, and the capacity will gradually decay. Increase, reduce safety, easily lead to unstable equipment operation and cause economic losses, and even lead to serious safety accidents. Therefore, the lithium battery needs to be maintained and replaced in time, and the remaining life prediction method can give the attenuation of the battery capacity in advance, so that a reasonable maintenance plan can be formulated in advance to ensure the stable operation of the electric forklift. Remaining life prediction refers to predicting the number of cycles that a battery can still run before reaching the end of its life, given the partial capacity decay data.

传统的剩余寿命预测方法是基于模型的方法,该方法基于电池失效原理和电化学反应建立模型,并对电池寿命衰减情况进行预测,该方法的缺点是计算复杂,精度低。The traditional remaining life prediction method is a model-based method, which establishes a model based on the battery failure principle and electrochemical reaction, and predicts the battery life decay. The disadvantage of this method is that the calculation is complicated and the accuracy is low.

发明内容SUMMARY OF THE INVENTION

本发明是为了克服现有技术预测锂电池剩余寿命时,计算复杂耗时长,预测精度低的问题,提供一种基于多神经网络耦合的电动叉车锂电池剩余寿命预测方法,提高了预测计算精度,减少了预测模型训练时间。In order to overcome the problems of complicated calculation and long time-consuming and low prediction accuracy when predicting the remaining life of a lithium battery in the prior art, the invention provides a method for predicting the remaining life of an electric forklift lithium battery based on multi-neural network coupling, which improves the prediction calculation accuracy. Reduced predictive model training time.

为了实现上述目的,本发明采用以下技术方案,包括以下步骤:In order to achieve the above object, the present invention adopts the following technical solutions, comprising the following steps:

S1:建立基于长短时记忆神经网络的开路电压预测模型,采用RMSprop算法和dropout正则化方法对网络进行优化,从而预测锂电池在放电循环中的开路电压值VOCS1: establish an open-circuit voltage prediction model based on a long-short-term memory neural network, and use the RMSprop algorithm and the dropout regularization method to optimize the network to predict the open-circuit voltage value V OC of the lithium battery during the discharge cycle;

S2:把预测结果按顺序划分成多个放电循环,统计每个放电循环中从初始电压至最小电压间的开路电压样本个数NS,利用采样时间TS相同,得到每个放电循环中放电至最小电压的时间TminS2: Divide the prediction results into multiple discharge cycles in sequence, count the number of open-circuit voltage samples N S from the initial voltage to the minimum voltage in each discharge cycle, and use the same sampling time T S to obtain the discharge cycle in each discharge cycle. the time T min to the minimum voltage;

S3:建立基于人工神经网络的容量预测模型,以预测锂电池容量C,从而得到锂电池剩余寿命预测值RUL。S3: Establish a capacity prediction model based on an artificial neural network to predict the capacity C of the lithium battery, thereby obtaining the predicted value RUL of the remaining life of the lithium battery.

通过结合长短时记忆神经网络和人工神经网络建立电动叉车锂电池剩余寿命预测模型,充分利用了长短时记忆神经网络对非线性数据的拟合能力和预测能力,利用了人工神经网络结构简洁的优点,又利用RMSprop算法和dropout正则化方法对预测模型进行优化,提高模型计算精度,减少模型训练时间。By combining the long-short-term memory neural network and artificial neural network, a prediction model for the remaining life of electric forklift lithium batteries is established, which makes full use of the long-short-term memory neural network's ability to fit and predict nonlinear data, and takes advantage of the simple structure of the artificial neural network. , and use the RMSprop algorithm and the dropout regularization method to optimize the prediction model, improve the calculation accuracy of the model, and reduce the model training time.

作为优选,所述S1的具体步骤如下:Preferably, the specific steps of the S1 are as follows:

S11:设定输入神经元个数为NI1,设定神经隐藏层层数为LH1和隐藏层神经元个数NH1,利用滑移窗口建立输入参数与输出参数之间的关系映射,设滑移窗口长度为L,设训练数据长度为M,则训练时建立输入Vinput到输出分别为Voutput的映射关系,如式(1)所示:S11: Set the number of input neurons to N I1 , set the number of neural hidden layers to L H1 and the number of hidden layer neurons to N H1 , use the sliding window to establish the relationship mapping between input parameters and output parameters, set The length of the sliding window is L, and the length of the training data is M, then the mapping relationship between the input V input and the output V output is established during training, as shown in formula (1):

Figure BDA0002548270080000021
Figure BDA0002548270080000021

其中,Vinput中的每一行对应一个输入样本,长度为L;Among them, each row in V input corresponds to an input sample, and the length is L;

S12:利用S11训练得到的模型对开路电压VOC采用如下方法进行预测,输入与输出如式(2)所示:S12: Use the model trained in S11 to predict the open-circuit voltage V OC using the following method, and the input and output are shown in formula (2):

Figure BDA0002548270080000031
Figure BDA0002548270080000031

其中,VOC中的每一行对应一次预测的输入样本,每次预测得到一个开路电压预测值

Figure BDA0002548270080000032
利用得到的预测值
Figure BDA0002548270080000033
对下一次预测的输入样本序列进行更新,即将预测值
Figure BDA0002548270080000034
加入到输入样本序列的末尾,并删掉输入样本序列的第一个数据,从而完成窗口的滑移操作,得到全部的开路电压预测值
Figure BDA0002548270080000035
作为优选,所述S2的具体步骤如下:Among them, each row in V OC corresponds to a predicted input sample, and each prediction obtains a predicted value of open circuit voltage
Figure BDA0002548270080000032
use the predicted value
Figure BDA0002548270080000033
Update the input sample sequence for the next prediction, that is, the predicted value
Figure BDA0002548270080000034
Add to the end of the input sample sequence, and delete the first data of the input sample sequence, so as to complete the sliding operation of the window and obtain all the predicted values of open circuit voltage
Figure BDA0002548270080000035
As preferably, the concrete steps of described S2 are as follows:

S21:把预测结果按顺序划分成多个放电循环的方法是:S21: The method of dividing the prediction result into multiple discharge cycles in sequence is:

若时刻Tj使得式(15)成立:If time T j makes equation (15) true:

Figure BDA0002548270080000036
Figure BDA0002548270080000036

则Tj(j=0,1,2,…,n;T0=1)为第j个循环的终止时刻,第j个循环的开路电压序列为

Figure BDA0002548270080000037
Then T j (j=0, 1, 2, ..., n; T 0 =1) is the termination time of the jth cycle, and the open-circuit voltage sequence of the jth cycle is
Figure BDA0002548270080000037

S22:统计每个放电循环中从初始电压至最小电压间的开路电压样本个数NS,利用采样时间TS相同,得到每个放电循环中放电至最小电压的时间Tmin,如式(16)所示:S22: Count the number of open-circuit voltage samples N S from the initial voltage to the minimum voltage in each discharge cycle, and use the same sampling time T S to obtain the time T min from discharging to the minimum voltage in each discharge cycle, as shown in formula (16 ) as shown:

Tmin=NS×TS (16)T min =N S ×T S (16)

作为优选,所述S3的具体步骤如下:As preferably, the concrete steps of described S3 are as follows:

S31:设定输入神经元个数为NI2,设定神经隐藏层层数为LH2和隐藏层神经元个数NH2,建立从放电循环中放电至最小电压的时间Tmin至锂电池容量C的映射关系;S31: Set the number of input neurons as N I2 , set the number of neural hidden layers as L H2 and the number of hidden layer neurons as N H2 , establish the time T min from the discharge cycle to the minimum voltage to the capacity of the lithium battery The mapping relationship of C;

S32:隐藏层各神经元计算过程如式(17)所示:S32: The calculation process of each neuron in the hidden layer is shown in formula (17):

HT=σ(WHTTmin+bH); (17)H T =σ(W HT T min +b H ); (17)

HT表示隐藏层计算函数,其中WHT和bH为该函数输入值的权值矩阵和偏置参数矩阵;迭代20次后结束训练,保存模型。将预测得到的每个放电循环中放电至最小电压的时间Tmin输入模型,得到容量预测结果;H T represents the hidden layer calculation function, where W HT and b H are the weight matrix and bias parameter matrix of the input value of the function; after 20 iterations, the training ends and the model is saved. Input the predicted time T min from discharging to the minimum voltage in each discharge cycle into the model to obtain the capacity prediction result;

S33:定义锂电池寿命终止点EOL为容量Cj等于额定容量C0的70%时所对应的循环数j;定义并计算锂电池剩余寿命RUL如式(18)所示:S33: Define the end point EOL of the lithium battery life as the cycle number j corresponding to when the capacity C j is equal to 70% of the rated capacity C 0 ; define and calculate the remaining life RUL of the lithium battery as shown in formula (18):

RUL=UL-EOL; (18)RUL=UL-EOL; (18)

其中,UL为当前锂电池已运行的循环数,EOL为锂电池的寿命终止点。Among them, UL is the number of cycles that the current lithium battery has run, and EOL is the end-of-life point of the lithium battery.

本发明的有益效果The beneficial effects of the present invention

本发明建立了长短时记忆-人工神经网络的多神经网络耦合模型,有效提高了预测的精度。利用RMSprop算法进行训练,提高了模型预测的准确度;本模型利用dropout方法对模型进行了优化,避免了过拟合问题,减少了模型训练的时间,提高了模型预测的准确度。利用了开路电压作为预测输入数据,与采用放电时间作为特征参数的长短时记忆-人工神经网络的耦合模型与简单长短时记忆神经网络模型相比预测结果的精度明显提高,能够准确预测电动叉车锂电池的剩余寿命。The invention establishes a long-short-term memory-artificial neural network multi-neural network coupling model, which effectively improves the prediction accuracy. The RMSprop algorithm is used for training, which improves the accuracy of model prediction; this model uses the dropout method to optimize the model, which avoids the problem of overfitting, reduces the time of model training, and improves the accuracy of model prediction. Using the open circuit voltage as the prediction input data, the accuracy of the prediction results is significantly improved compared with the long-short-term memory-artificial neural network coupling model using the discharge time as the characteristic parameter and the simple long-short-term memory neural network model, which can accurately predict the lithium ion of the electric forklift. remaining battery life.

附图说明Description of drawings

图1是本发明的总流程框图;Fig. 1 is the general flow chart of the present invention;

图2是本发明的第一组预测结果的分析图;Fig. 2 is the analysis diagram of the first group of prediction results of the present invention;

图3是本发明的第二组预测结果的分析图;Fig. 3 is the analysis diagram of the second group of prediction results of the present invention;

图4为剩余寿命预测的绝对误差;Figure 4 shows the absolute error of remaining life prediction;

图5为容量衰减预测的均方根误差。Figure 5 shows the root mean square error of capacity decay predictions.

具体实施方式Detailed ways

下面结合附图和具体实施方式对本发明作详细说明。The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.

本发明的基于多神经网络耦合的电动叉车锂电池剩余寿命预测方法,总流程框图如图1所示,具体操作过程包括:The overall flow diagram of the method for predicting the remaining life of an electric forklift lithium battery based on multi-neural network coupling of the present invention is shown in Figure 1, and the specific operation process includes:

S1:建立基于长短时记忆神经网络的开路电压预测模型,采用RMSprop算法和dropout正则化方法对网络进行优化,从而预测锂电池在放电循环中的开路电压值VOC,具体步骤如下:S1: Establish an open-circuit voltage prediction model based on a long-short-term memory neural network, and use the RMSprop algorithm and dropout regularization method to optimize the network to predict the open-circuit voltage value V OC of the lithium battery during the discharge cycle. The specific steps are as follows:

S11:设定输入神经元个数为NI1,设定神经隐藏层层数为LH1和隐藏层神经元个数NH1,利用滑移窗口建立输入参数与输出参数之间的关系映射,设滑移窗口长度为L,设训练数据长度为M,则训练时建立输入Vinput到输出分别为Voutput的映射关系,如式(1)所示:S11: Set the number of input neurons to N I1 , set the number of neural hidden layers to L H1 and the number of hidden layer neurons to N H1 , use the sliding window to establish the relationship mapping between input parameters and output parameters, set The length of the sliding window is L, and the length of the training data is M, then the mapping relationship between the input V input and the output V output is established during training, as shown in formula (1):

Figure BDA0002548270080000051
Figure BDA0002548270080000051

其中,Vinput中的每一行对应一个输入样本,长度为L;Among them, each row in V input corresponds to an input sample, and the length is L;

S12:利用S11训练得到的模型对开路电压VOC采用如下方法进行预测,输入与输出如式(2)所示:S12: Use the model trained in S11 to predict the open-circuit voltage V OC using the following method, and the input and output are shown in formula (2):

Figure BDA0002548270080000052
Figure BDA0002548270080000052

其中,VOC中的每一行对应一次预测的输入样本,每次预测得到一个开路电压预测值

Figure BDA0002548270080000053
利用得到的预测值
Figure BDA0002548270080000054
对下一次预测的输入样本序列进行更新,即将预测值
Figure BDA0002548270080000055
加入到输入样本序列的末尾,并删掉输入样本序列的第一个数据,从而完成窗口的滑移操作,得到全部的开路电压预测值
Figure BDA0002548270080000061
Among them, each row in V OC corresponds to a predicted input sample, and each prediction obtains a predicted value of open circuit voltage
Figure BDA0002548270080000053
use the predicted value
Figure BDA0002548270080000054
Update the input sample sequence for the next prediction, that is, the predicted value
Figure BDA0002548270080000055
Add to the end of the input sample sequence, and delete the first data of the input sample sequence, so as to complete the sliding operation of the window and obtain all the predicted values of open circuit voltage
Figure BDA0002548270080000061

长短时记忆神经网络对于时序数据具有很好的预测及拟合能力,为了防止过拟合,对该网络进行dropout正则化处理,即随机产生0和1向量,从而随机丢弃网络中的部分神经元,如式(3)和(4)所示:The long-short-term memory neural network has good prediction and fitting ability for time series data. In order to prevent over-fitting, dropout regularization is performed on the network, that is, 0 and 1 vectors are randomly generated, thereby randomly discarding some neurons in the network. , as shown in equations (3) and (4):

Figure BDA0002548270080000062
Figure BDA0002548270080000062

Figure BDA0002548270080000063
Figure BDA0002548270080000063

其中RT服从伯努利(Bernoulli)分布函数,p为dropout正则化比率;每个长短时记忆神经网络的神经元中包括输入门、遗忘门和输出门三种信息的选取及转换机制,实现过程如下:Among them, R T obeys the Bernoulli distribution function, and p is the dropout regularization ratio; the neurons of each long and short-term memory neural network include three kinds of information selection and conversion mechanisms: input gate, forget gate and output gate. The process is as follows:

遗忘门决定来自于前一个神经元的信息是否被保留,取值为1或0,前者代表保留,后者代表丢弃,如式(5)所示:The forget gate determines whether the information from the previous neuron is retained, and takes a value of 1 or 0. The former represents retention and the latter represents discard, as shown in formula (5):

Figure BDA0002548270080000064
Figure BDA0002548270080000064

其中Wfv和Wfh分别为遗忘门函数的当前神经元的输入值和前一神经元的输出值的权值矩阵,bf为偏置参数矩阵;where W fv and W f h are the weight matrix of the input value of the current neuron of the forget gate function and the output value of the previous neuron, respectively, and b f is the bias parameter matrix;

输入门计算当前神经元输入的保留信息rT和新的状态信息sT,如式(6)和(7)所示:The input gate calculates the retained information r T and the new state information s T of the current neuron input, as shown in equations (6) and (7):

Figure BDA0002548270080000065
Figure BDA0002548270080000065

Figure BDA0002548270080000066
Figure BDA0002548270080000066

其中,Wrv和Wsv分别为输入门中函数rT和sT当前神经元输入值的权值矩阵,Wrh和Wsh分别为输入门中函数rT和sT前一神经元输出值的权值矩阵,br和bs分别为输入门中函数rT和sT的偏置参数矩阵;Among them, W rv and W sv are the weight matrices of the current neuron input values of functions r T and s T in the input gate, respectively, and W rh and W sh are the output values of the previous neuron of functions r T and s T in the input gate, respectively The weight matrix of , b r and b s are the bias parameter matrices of the functions r T and s T in the input gate, respectively;

之后对长短时记忆神经网络神经元的状态cT做如式(8)更新:Afterwards, the state c T of the long-short-term memory neural network neurons is updated as in Equation (8):

cT=rT⊙sT+cT-1⊙fT (8)c T =r T ⊙s T +c T-1 ⊙f T (8)

输出门计算当前神经元的输出信息hT,如式(9)和(10)所示:The output gate calculates the output information h T of the current neuron, as shown in equations (9) and (10):

Figure BDA0002548270080000071
Figure BDA0002548270080000071

hT=oT⊙tanh(cT) (10)h T =o T ⊙tanh(c T ) (10)

其中Wov和Woh为输出门中函数oT当前神经元输入值和前一神经元输出值的权值矩阵,bo为函数oT的偏置参数矩阵;where W ov and W oh are the weight matrix of the input value of the current neuron and the output value of the previous neuron in the function o T in the output gate, and b o is the bias parameter matrix of the function o T ;

利用RMSprop算法对长短时记忆神经网络的权值矩阵ω和偏置参数矩阵b做如式(11-14)更新:Use the RMSprop algorithm to update the weight matrix ω and bias parameter matrix b of the long-short-term memory neural network as shown in Equation (11-14):

Figure BDA0002548270080000072
Figure BDA0002548270080000072

Figure BDA0002548270080000073
Figure BDA0002548270080000073

Figure BDA0002548270080000074
Figure BDA0002548270080000074

Figure BDA0002548270080000075
Figure BDA0002548270080000075

其中J1为计算权值矩阵的成本函数,J2为计算偏置参数矩阵的成本函数,k代表当前迭代的次数,β为系数,通常取0.999,ε用于防止除数为零,通常取10-8,α为学习速率,该算法属于加快梯度下降法,因此经20次迭代后长短时记忆神经网络的训练结束,将模型保存;输入已知数据,得到开路电压预测序列。where J 1 is the cost function for calculating the weight matrix, J 2 is the cost function for calculating the bias parameter matrix, k represents the number of current iterations, β is the coefficient, usually 0.999, ε is used to prevent division by zero, usually 10 -8 , α is the learning rate, the algorithm belongs to the accelerated gradient descent method, so after 20 iterations, the training of the long-short-term memory neural network is over, and the model is saved; input the known data to get the open-circuit voltage prediction sequence.

S2:把预测结果按顺序划分成多个放电循环,统计每个放电循环中从初始电压至最小电压间的开路电压样本个数NS,利用采样时间TS相同,得到每个放电循环中放电至最小电压的时间Tmin,具体步骤如下:S2: Divide the prediction results into multiple discharge cycles in sequence, count the number of open-circuit voltage samples N S from the initial voltage to the minimum voltage in each discharge cycle, and use the same sampling time T S to obtain the discharge cycle in each discharge cycle. The time T min to the minimum voltage, the specific steps are as follows:

S21:把预测结果按顺序划分成多个放电循环的方法是:S21: The method of dividing the prediction result into multiple discharge cycles in sequence is:

若时刻Tj使得式(15)成立:If time T j makes equation (15) true:

Figure BDA0002548270080000076
Figure BDA0002548270080000076

则Tj(j=0,1,2,…,n;T0=1)为第j个循环的终止时刻,第j个循环的开路电压序列为

Figure BDA0002548270080000081
Then T j (j=0, 1, 2, ..., n; T 0 =1) is the termination time of the jth cycle, and the open-circuit voltage sequence of the jth cycle is
Figure BDA0002548270080000081

S22:统计每个放电循环中从初始电压至最小电压间的开路电压样本个数NS,利用采样时间TS相同,得到每个放电循环中放电至最小电压的时间Tmin,如式(16)所示:S22: Count the number of open-circuit voltage samples N S from the initial voltage to the minimum voltage in each discharge cycle, and use the same sampling time T S to obtain the time T min from discharging to the minimum voltage in each discharge cycle, as shown in formula (16 ) as shown:

Tmin=NS×TS (16)T min =N S ×T S (16)

S3:建立基于人工神经网络的容量预测模型,以预测锂电池容量C,从而得到锂电池剩余寿命预测值RUL,具体步骤如下:S3: Establish a capacity prediction model based on an artificial neural network to predict the capacity C of the lithium battery, thereby obtaining the predicted value RUL of the remaining life of the lithium battery. The specific steps are as follows:

S31:设定输入神经元个数为NI2,设定神经隐藏层层数为LH2和隐藏层神经元个数NH2,建立从放电循环中放电至最小电压的时间Tmin至锂电池容量C的映射关系;S31: Set the number of input neurons as N I2 , set the number of neural hidden layers as L H2 and the number of hidden layer neurons as N H2 , establish the time T min from the discharge cycle to the minimum voltage to the capacity of the lithium battery The mapping relationship of C;

S32:隐藏层各神经元计算过程如式(17)所示:S32: The calculation process of each neuron in the hidden layer is shown in formula (17):

HT=σ(WHTTmin+bH); (17)H T =σ(W HT T min +b H ); (17)

HT表示隐藏层计算函数,其中WHT和bH为该函数输入值的权值矩阵和偏置参数矩阵;迭代20次后结束训练,保存模型。将预测得到的每个放电循环中放电至最小电压的时间Tmin输入模型,得到容量预测结果;H T represents the hidden layer calculation function, where W HT and b H are the weight matrix and bias parameter matrix of the input value of the function; after 20 iterations, the training ends and the model is saved. Input the predicted time T min from discharging to the minimum voltage in each discharge cycle into the model to obtain the capacity prediction result;

S33:定义锂电池寿命终止点EOL为容量Cj等于额定容量C0的70%时所对应的循环数j;定义并计算锂电池剩余寿命RUL,如式(18)所示:S33: Define the end point EOL of the lithium battery life as the cycle number j corresponding to when the capacity C j is equal to 70% of the rated capacity C 0 ; define and calculate the remaining life RUL of the lithium battery, as shown in formula (18):

RUL=UL-EOL; (18)RUL=UL-EOL; (18)

其中,UL为当前锂电池已运行的循环数,EOL为锂电池的寿命终止。Among them, UL is the number of cycles that the current lithium battery has run, and EOL is the end of life of the lithium battery.

本发明利用Python语言进行建模,分别验证了上述方法对于NASA AmesPrognostics Center of Excellence(PCoE)和Center for Advanced Life CycleEngineering(CALCE)提供的锂电池寿命衰减数据的预测结果,并与采用放电时间作为特征参数的长短时记忆-人工神经网络的耦合模型和简单长短时记忆神经网络模型的预测结果进行了对比,各组实验均采用30%数据用于训练,预测结果如图2和图3所示(图2与图3所用锂电池样本不同),其中M1为以开路电压作为特征参数的长短时记忆-人工神经网络耦合模型的预测结果,M2为以放电时间作为特征参数的长短时记忆-人工神经网络耦合模型的预测结果,M3为简单长短时记忆神经网络模型的预测结果,可以看出本发明建立的以开路电压作为特征参数的长短时记忆-人工神经网络的耦合模型和其他模型相比具有更好的预测效果,通过对不同类型的电池进行预测,可以看出本发明建立的以开路电压作为特征参数的长短时记忆-人工神经网络的耦合模型具有很好的泛化性能。为进一步说明,利用绝对误差(AE,absolute error)和均方根误差(RMSE,root mean squared error)对预测结果进行评价,计算过程如式(19)和(20)所示:The present invention uses Python language for modeling, respectively verifies the prediction results of the above methods for the lithium battery life decay data provided by NASA AmesPrognostics Center of Excellence (PCoE) and Center for Advanced Life CycleEngineering (CALCE), and uses the discharge time as a feature. The prediction results of the long-short-term memory-artificial neural network coupling model of parameters and the simple long-short-term memory neural network model were compared. Each group of experiments used 30% data for training, and the prediction results are shown in Figure 2 and Figure 3 ( Figure 2 is different from the lithium battery sample used in Figure 3), where M1 is the prediction result of the long-short-term memory-artificial neural network coupling model with open circuit voltage as the characteristic parameter, M2 is the long-short-term memory-artificial neural network with discharge time as the characteristic parameter The prediction result of the network coupling model, M3 is the prediction result of a simple long-short-term memory neural network model, it can be seen that the long-short-term memory-artificial neural network coupling model established by the present invention with open circuit voltage as a characteristic parameter has a better performance than other models. Better prediction effect. By predicting different types of batteries, it can be seen that the long-short-term memory-artificial neural network coupling model with open circuit voltage as a characteristic parameter established by the present invention has good generalization performance. For further explanation, the prediction results are evaluated by absolute error (AE, absolute error) and root mean squared error (RMSE, root mean squared error). The calculation process is shown in equations (19) and (20):

Figure BDA0002548270080000091
Figure BDA0002548270080000091

Figure BDA0002548270080000092
Figure BDA0002548270080000092

其中,RUL为电池剩余循环寿命的真实值,

Figure BDA0002548270080000093
为电池剩余循环寿命的预测值,Cj为第j个循环电池容量的真实值,
Figure BDA0002548270080000094
为第j个循环电池容量的预测值。计算结果如图4和图5所示(图4为剩余寿命预测的绝对误差,图5为容量衰减预测的均方根误差),Among them, RUL is the true value of the remaining cycle life of the battery,
Figure BDA0002548270080000093
is the predicted value of the remaining cycle life of the battery, C j is the actual value of the battery capacity in the jth cycle,
Figure BDA0002548270080000094
is the predicted value of the battery capacity for the jth cycle. The calculation results are shown in Figure 4 and Figure 5 (Figure 4 is the absolute error of the remaining life prediction, Figure 5 is the root mean square error of the capacity decay prediction),

本发明建立的以开路电压作为特征参数的长短时记忆-人工神经网络的耦合模型,对锂电池剩余寿命预测的误差较小,且对实际容量衰减曲线拟合精度更高,具有非常高的剩余寿命预测精度。The long-short-term memory-artificial neural network coupling model with open circuit voltage as the characteristic parameter established by the present invention has smaller error in predicting the remaining life of lithium battery, and has higher fitting accuracy for the actual capacity decay curve, and has a very high residual life. Lifetime prediction accuracy.

上述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。The foregoing embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that the foregoing embodiments can still be used for The recorded technical solutions are modified, or some or all of the technical features thereof are equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.

Claims (4)

1. A method for predicting the remaining life of a lithium battery of an electric forklift based on multi-neural network coupling is characterized by comprising the following steps:
s1: establishing an open-circuit voltage prediction model based on a long-time and short-time memory neural network, and adopting an RMSprop algorithm and a dOptimizing the network by using a ropout regularization method, thereby predicting the open-circuit voltage value V of the lithium battery in the discharge cycleOC
S2: dividing the prediction result into a plurality of discharge cycles in sequence, and counting the number N of open-circuit voltage samples from the initial voltage to the minimum voltage in each discharge cycleSUsing the sampling time TSThe same, the time T to discharge to the minimum voltage in each discharge cycle is obtainedmin
S3: and establishing a capacity prediction model based on an artificial neural network to predict the capacity C of the lithium battery so as to obtain a predicted value RUL of the residual life of the lithium battery.
2. The method for predicting the residual life of the lithium battery of the electric forklift based on the multi-neural-network coupling as claimed in claim 1, wherein the specific steps of S1 are as follows:
s11: setting the number of input neurons to NI1Setting the number of the nerve hiding layers to be LH1And number of hidden layer neurons NH1Establishing a relation mapping between input parameters and output parameters by using a sliding window, and establishing an input V during training if the length of the sliding window is L and the length of training data is MinputTo outputs of V respectivelyoutputThe mapping relation of (2) is shown in formula (1):
Figure FDA0002548270070000011
wherein, VinputEach row in the set corresponds to an input sample with a length of L;
s12: model obtained by training S11 is utilized to split open-circuit voltage VOCThe prediction is carried out by adopting the following method, wherein the input and the output are shown as the formula (2):
Figure FDA0002548270070000021
wherein, VOCOne for each row inEach time of prediction to obtain a predicted value of open-circuit voltage
Figure FDA0002548270070000022
(t ═ M +1, M +2, …, N), using the predicted values obtained
Figure FDA0002548270070000023
Updating the input sample sequence for the next prediction, i.e. the predicted value
Figure FDA0002548270070000024
Adding the data to the end of the input sample sequence, and deleting the first data of the input sample sequence, thereby completing the window sliding operation and obtaining all predicted values of the open-circuit voltage
Figure 1
3. The method for predicting the residual life of the lithium battery of the electric forklift based on the multi-neural-network coupling as claimed in claim 1, wherein the specific steps of S2 are as follows:
s21: the method for sequentially dividing the prediction result into a plurality of discharge cycles is as follows:
if the time TjSo that equation (15) holds:
Figure FDA0002548270070000026
then Tj(j=0,1,2,…,n;T01) is the end time of the jth cycle, and the open circuit voltage sequence of the jth cycle is
Figure FDA0002548270070000027
S22: counting the number N of open-circuit voltage samples from the initial voltage to the minimum voltage in each discharge cycleSUsing the sampling time TSSame, get each discharge cycleTime T for discharging to minimum voltageminAs shown in formula (16):
Tmin=NS×TS (16)
4. the method for predicting the residual life of the lithium battery of the electric forklift based on the multi-neural-network coupling as claimed in claim 1, wherein the specific steps of S3 are as follows:
s31: setting the number of input neurons to NI2Setting the number of the nerve hiding layers to be LH2And number of hidden layer neurons NH2Establishing the time T from discharge to minimum voltage in the discharge cycleminMapping relation to lithium battery capacity C;
s32: the calculation process of each neuron of the hidden layer is shown as the formula (17):
HT=σ(WHTTmin+bH); (17)
HTrepresenting a hidden layer computation function, where WHTAnd bHInputting a weight matrix and a bias parameter matrix of the function input value; after 20 times of iteration, the training is finished, and the model is saved. The predicted time T for discharging to the minimum voltage in each discharge cycleminInputting a model to obtain a capacity prediction result;
s33: defining EOL as capacity CjEqual to rated capacity C0The number j of cycles corresponding to 70%; defining and calculating the residual life RUL of the lithium battery, as shown in formula (18):
RUL=UL-EOL; (18)
wherein, UL is the number of the operating cycles of the current lithium battery, and EOL is the end point of the service life of the lithium battery.
CN202010568146.1A 2020-06-19 2020-06-19 A method for predicting the remaining life of electric forklift lithium battery based on multi-neural network coupling Active CN111948563B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010568146.1A CN111948563B (en) 2020-06-19 2020-06-19 A method for predicting the remaining life of electric forklift lithium battery based on multi-neural network coupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010568146.1A CN111948563B (en) 2020-06-19 2020-06-19 A method for predicting the remaining life of electric forklift lithium battery based on multi-neural network coupling

Publications (2)

Publication Number Publication Date
CN111948563A true CN111948563A (en) 2020-11-17
CN111948563B CN111948563B (en) 2021-12-17

Family

ID=73337098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010568146.1A Active CN111948563B (en) 2020-06-19 2020-06-19 A method for predicting the remaining life of electric forklift lithium battery based on multi-neural network coupling

Country Status (1)

Country Link
CN (1) CN111948563B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731184A (en) * 2020-12-28 2021-04-30 深圳供电局有限公司 Battery service life detection method and system
CN112986827A (en) * 2021-04-12 2021-06-18 山东凯格瑞森能源科技有限公司 Fuel cell residual life prediction method based on deep learning
CN113359048A (en) * 2021-04-28 2021-09-07 中国矿业大学 Indirect prediction method for remaining service life of lithium ion battery
CN113777496A (en) * 2021-09-06 2021-12-10 北京化工大学 Remaining life prediction method of lithium-ion battery based on temporal convolutional neural network
CN114565008A (en) * 2022-01-12 2022-05-31 西安理工大学 Anomaly detection method of lithium battery based on long-short-term memory autoencoder
CN114859249A (en) * 2022-07-06 2022-08-05 苏州清研精准汽车科技有限公司 Method and device for detecting battery pack capacity
CN115186370A (en) * 2022-05-18 2022-10-14 广东海洋大学 Engineering forklift transfer learning system based on deep learning training model
CN118861604A (en) * 2024-09-23 2024-10-29 西安君扬红芯电子股份有限公司 A life prediction method for pulse power supply system

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09243716A (en) * 1996-03-08 1997-09-19 Sanken Electric Co Ltd Method of judging state of secondary battery
JP2014206499A (en) * 2013-04-15 2014-10-30 新電元工業株式会社 Secondary battery life prediction system, secondary battery characteristic evaluation device, secondary battery mounting device, and secondary battery life prediction method
CN104849671A (en) * 2015-05-22 2015-08-19 大连理工大学 Battery pack capacity detection system based on combined neural network
CN107167741A (en) * 2017-06-06 2017-09-15 浙江大学 A kind of lithium battery SOC observation procedures based on neutral net
CN109343505A (en) * 2018-09-19 2019-02-15 太原科技大学 Prediction method of gear remaining life based on long short-term memory network
CN109993270A (en) * 2019-03-27 2019-07-09 东北大学 Remaining life prediction method of lithium-ion battery based on gray wolf pack optimization LSTM network
CN110542866A (en) * 2019-10-12 2019-12-06 上海新微技术研发中心有限公司 method for estimating residual electric quantity parameter of battery
CN110687460A (en) * 2019-10-12 2020-01-14 上海新微技术研发中心有限公司 Soc estimation method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09243716A (en) * 1996-03-08 1997-09-19 Sanken Electric Co Ltd Method of judging state of secondary battery
JP2014206499A (en) * 2013-04-15 2014-10-30 新電元工業株式会社 Secondary battery life prediction system, secondary battery characteristic evaluation device, secondary battery mounting device, and secondary battery life prediction method
CN104849671A (en) * 2015-05-22 2015-08-19 大连理工大学 Battery pack capacity detection system based on combined neural network
CN107167741A (en) * 2017-06-06 2017-09-15 浙江大学 A kind of lithium battery SOC observation procedures based on neutral net
CN109343505A (en) * 2018-09-19 2019-02-15 太原科技大学 Prediction method of gear remaining life based on long short-term memory network
CN109993270A (en) * 2019-03-27 2019-07-09 东北大学 Remaining life prediction method of lithium-ion battery based on gray wolf pack optimization LSTM network
CN110542866A (en) * 2019-10-12 2019-12-06 上海新微技术研发中心有限公司 method for estimating residual electric quantity parameter of battery
CN110687460A (en) * 2019-10-12 2020-01-14 上海新微技术研发中心有限公司 Soc estimation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YONGZHI ZHANG 等: "Long Short-Term Memory Recurrent Neural Network for Remaining Useful Life Prediction of Lithium-Ion Batteries", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》 *
李练兵 等: "锂离子电池组可用剩余容量的估算方法", 《电源技术》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731184A (en) * 2020-12-28 2021-04-30 深圳供电局有限公司 Battery service life detection method and system
CN112731184B (en) * 2020-12-28 2023-03-03 深圳供电局有限公司 Battery service life detection method and system
CN112986827A (en) * 2021-04-12 2021-06-18 山东凯格瑞森能源科技有限公司 Fuel cell residual life prediction method based on deep learning
CN112986827B (en) * 2021-04-12 2022-06-03 山东凯格瑞森能源科技有限公司 Fuel cell residual life prediction method based on deep learning
CN113359048A (en) * 2021-04-28 2021-09-07 中国矿业大学 Indirect prediction method for remaining service life of lithium ion battery
CN113777496A (en) * 2021-09-06 2021-12-10 北京化工大学 Remaining life prediction method of lithium-ion battery based on temporal convolutional neural network
CN113777496B (en) * 2021-09-06 2023-10-24 北京化工大学 Lithium ion battery residual life prediction method based on time convolution neural network
CN114565008A (en) * 2022-01-12 2022-05-31 西安理工大学 Anomaly detection method of lithium battery based on long-short-term memory autoencoder
CN115186370A (en) * 2022-05-18 2022-10-14 广东海洋大学 Engineering forklift transfer learning system based on deep learning training model
CN114859249A (en) * 2022-07-06 2022-08-05 苏州清研精准汽车科技有限公司 Method and device for detecting battery pack capacity
CN114859249B (en) * 2022-07-06 2022-10-18 苏州清研精准汽车科技有限公司 Method and device for detecting battery pack capacity
CN118861604A (en) * 2024-09-23 2024-10-29 西安君扬红芯电子股份有限公司 A life prediction method for pulse power supply system

Also Published As

Publication number Publication date
CN111948563B (en) 2021-12-17

Similar Documents

Publication Publication Date Title
CN111948563A (en) A method for predicting the remaining life of electric forklift lithium battery based on multi-neural network coupling
CN112241608B (en) Lithium battery life prediction method based on LSTM network and transfer learning
CN110824364B (en) Lithium battery SOH estimation and RUL prediction method based on AST-LSTM neural network
CN109993270B (en) Lithium-ion battery remaining life prediction method based on gray wolf pack optimized LSTM network
CN109991542B (en) Lithium ion battery residual life prediction method based on WDE optimization LSTM network
CN114861527A (en) Lithium battery life prediction method based on time series characteristics
CN111308375B (en) LSTM-FFNN-based electric forklift lithium ion battery health state prediction method
CN112036084B (en) A kind of similar product lifetime migration screening method and system
CN111736084B (en) Valve-regulated lead-acid storage battery health state prediction method based on improved LSTM neural network
CN110568359B (en) Lithium battery residual life prediction method
CN114325450A (en) Lithium ion battery health state prediction method based on CNN-BilSTM-AT hybrid model
CN108519556A (en) A Lithium-ion Battery SOC Prediction Method Based on Recurrent Neural Network
CN113064093A (en) Energy storage battery state of charge and state of health joint estimation method and system
Dai Nguyen et al. Equivalent circuit simulated deep network architecture and transfer learning for remaining useful life prediction of lithium-ion batteries
CN115453399B (en) Battery pack SOH estimation method considering inconsistency
CN111999648A (en) A method for predicting the remaining life of lithium batteries based on long short-term memory network
CN108549036A (en) Ferric phosphate lithium cell life-span prediction method based on MIV and SVM models
CN113361692B (en) Lithium battery remaining life combined prediction method
CN111812519B (en) A battery parameter identification method and system
CN115308606B (en) A method for estimating the health state of lithium-ion batteries based on neighboring features
CN115047350B (en) Digital-analog linkage based lithium ion battery remaining service life prediction method
CN115032544B (en) Method for predicting remaining life of power lithium battery under incomplete data scene
CN116643196A (en) Battery health state estimation method integrating mechanism and data driving model
CN112734002A (en) Service life prediction method based on data layer and model layer joint transfer learning
CN118169582B (en) Lithium ion battery health state and residual life prediction method

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Tong Shuiguang

Inventor after: Tong Zheming

Inventor after: Li Yuansong

Inventor after: Miao Jiazhi

Inventor before: Tong Zheming

Inventor before: Miao Jiazhi

Inventor before: Tong Shuiguang

Inventor before: Li Yuansong

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20211122

Address after: 310000 Yuhang Tang Road, Xihu District, Hangzhou, Zhejiang 866

Applicant after: ZHEJIANG University

Applicant after: Hangcha Group Co., Ltd

Address before: 310058 Yuhang Tang Road, Xihu District, Hangzhou, Zhejiang 866

Applicant before: ZHEJIANG University

GR01 Patent grant
GR01 Patent grant
点击 这是indexloc提供的php浏览器服务,不要输入任何密码和下载