Invention content
It is a primary object of the present invention in view of the deficiencies of the prior art, provide a kind of part charge and discharge based on secondary cell
The method and system of the electric complete charge and discharge data of data reconstruction, complete charge and discharge can be rebuild using part charge and discharge electrical measurement data
Data and curves shorten the time of the techniques such as partial volume, stepping, improve production efficiency.
To achieve the above object, the present invention uses following technical scheme:
A method of based on the complete charge and discharge data of battery portion charge and discharge data reconstruction, the battery especially lithium electricity
Pond the described method comprises the following steps:
A, the part charge and discharge electrical measurement data curve for including multiple discrete sampling point datas of the battery is obtained;
B, the relatively preset Global center reference data curve of the part charge and discharge electrical measurement data curve is determined
Starting point;
C, using the part charge and discharge electrical measurement data curve and preset a plurality of stepping center reference data and curves,
By linearly rebuild or least square rebuild, determine curve of approximation, and according to the curve of approximation and the part charge and discharge
The starting point of the relatively preset Global center reference data curve of measurement data curve, polishing remove the part charge and discharge
The charge and discharge electrical measurement data curve of remainder except electrical measurement data curve, to rebuild complete charge and discharge electrical measurement data
Curve.
Further, step A includes being pre-processed to the part charge and discharge electrical measurement data curve, the pretreatment packet
Include following steps:
A1, fast Fourier (FFT) transformation is carried out to the multiple discrete sampling point data, time domain data is transformed into frequency
Domain;
A2, windowing process progress smothing filtering is carried out to frequency domain data;
A3, zero padding interpolation is carried out to filtered frequency domain data, keeps its sample rate bent with the Global center reference data
Line and the data sampling rate of the stepping center reference data and curves are consistent;
A4, by the transform frequency domain data after the smothing filtering and the interpolation to time domain.
Further, in step B, by finding the part charge and discharge electrical measurement data in preset sliding window
The least absolute value of curve and the difference of the Global center reference data curve, to determine the part charge and discharge electrical measurement data
The starting point of the relatively described Global center reference data curve of curve.
Further, step B includes:
B1, using the set point on the Global center reference data curve as starting point, correspond to sampling sequence number w, take and set
Corresponding measurement data subtraction calculations difference on the one piece of data of measured length, with the part charge and discharge electrical measurement data curve, will
Difference after subtracting each other adds up and calculates the absolute value of accumulated value, and wherein w is initialized as 1;
B2, by the absolute value of the accumulated value calculated in step B1 and the difference minimum value e_min for being initially set maximal number
Compare, if it is less than e_min, then using e_min as the absolute value of the accumulated value calculated in step B1;
W plus 1 if w is more than the size of the sliding window, and is terminated, is otherwise returned by B3, the update starting point
Step B1.
Further, step C includes the following steps:
C11. by the part charge and discharge electrical measurement data curve respectively in a plurality of stepping center reference data and curves
Each matched, determine and the stepping center reference data of part charge and discharge electrical measurement data curve difference value minimum song
Line;
C12. use the stepping center reference data that step C11 is determined will be described according to the starting point as curve of approximation
The charge and discharge electrical measurement data curve of remainder described in polishing, it is bent to rebuild complete charge and discharge electrical measurement data after curve of approximation translation
Line.
Further, step C includes the following steps:
C21. it is gone to approach the part charge and discharge electrical measurement data curve with a plurality of stepping center reference data and curves, really
Surely the optimal coefficient for minimizing the quadratic sum of approximate error;
C22. according to the optimal coefficient, curve of approximation is established using a plurality of stepping center reference data and curves, it will
The charge and discharge electrical measurement data curve of remainder described in polishing, rebuilds complete charge and discharge electrical measurement number after the curve of approximation translation
According to curve.
Further, in step C21,
If a plurality of stepping center reference data and curves are L curve, if every stepping center reference data and curves have N
A discrete sampling point data,
Including the part charge and discharge electrical measurement data of N number of discrete sampling point data is with vector
Completely charge and discharge data are M of l stepping center reference data and curves with vector
rl=[rl(1) … rl(N-1)rl(N) … rl(M)]T
The matrix of L × M of L stepping reference data composition is
L stepping reference data is respectively with pstdFor starting point, the matrix of the L × N for the data composition that length is N is
If undetermined coefficient vector is h=[h(1) … h(L)]T, least square approximation is carried out with following formula
Optimal coefficient h is solved to
After determining optimal coefficient, the approximating curve rebuild using L stepping center reference data and curves is
Further, the Global center reference data curve and a plurality of stepping center reference data and curves are according to electricity
The material and production technology in pond are preset.
A kind of computer readable storage medium, is stored thereon with computer program, and the computer program is held by processor
The method based on the complete charge and discharge data of battery portion charge and discharge data reconstruction is realized when row.
A kind of system based on the complete charge and discharge data of battery portion charge and discharge data reconstruction, including:
Processor;And
Memory, the executable instruction for storing the processor;
Wherein, the processor is configured to execute via the executable instruction is executed described is filled based on battery portion
The method that discharge data rebuilds complete charge and discharge data.
Beneficial effects of the present invention:
In order to shorten the time of the techniques such as partial volume, inner walkway, production efficiency is improved, the present invention provides a kind of based on part
The method of the complete charge and discharge data of charge and discharge data reconstruction.The method of the present invention is based on preset reference data curve, only
Charge and discharge electrical measurement data in part need to be utilized high-precision can to rebuild complete charge and discharge data and curves.The method of the present invention is applied to
When partial volume technique, the lithium cell of band part electricity need to only be charged to fullcharging electricity condition, subsequent section electric discharge is required up to manufacture, i.e.,
Using the complete charge and discharge data and curves of data reconstruction of part charging and partial discharge, to complete subsequent charge/discharge capacity
Estimation and battery stepping eliminate the cycle of intermediate primary complete charge and discharge, shorten for 50% or more time, substantially increase
Production efficiency.
The aforementioned feature and technical advantage for quite widely having elaborated the present invention, so as to more fully understand this hair
Bright detailed description.The other feature and advantage of the present invention will be described below.
Specific implementation mode
It elaborates below to embodiments of the present invention.It is emphasized that following the description is only exemplary,
The range being not intended to be limiting of the invention and its application.
In conjunction with attached drawing, the present invention is further described in detail by the following examples.It is emphasized that following
Explanation is only exemplary, the range being not intended to be limiting of the invention and its application.It should be appreciated by those skilled in the art that
The concept and specific embodiment of disclosure are used for changing or design other structures to complete based on being easily used as
The identical purpose of the present invention.Those skilled in the art are not it should also be appreciated that this equivalent construction offsets the present invention's
Spirit and scope.Be considered as the novel features of feature of the present invention, structure and operation method, and further purpose and
Advantage will be better understood when from description below and in conjunction with attached drawing.However, it should profoundly recognize, each spy provided
Sign is all merely to description and explanation, the definition without being intended to the limitation present invention.
Refering to fig. 1, in one embodiment, a kind of based on the complete charge and discharge data of battery portion charge and discharge data reconstruction
Method, the battery especially lithium battery, the described method comprises the following steps:
Step A, the part charge and discharge electrical measurement data curve for including multiple discrete sampling point datas of the battery is obtained;
Step B, determine that the relatively preset Global center reference data of the part charge and discharge electrical measurement data curve is bent
The starting point of line;
Step C, the part charge and discharge electrical measurement data curve and preset a plurality of stepping center reference data are utilized
Curve, by linearly rebuild or least square rebuild, determine curve of approximation, and fill according to the curve of approximation and the part
The starting point of the relatively preset Global center reference data curve of discharge measuring data and curves, polishing remove the part
The charge and discharge electrical measurement data curve of remainder except charge and discharge electrical measurement data curve, to rebuild complete charge and discharge electrical measurement
Data and curves.
In various embodiments, charge and discharge electrical measurement data curve can be but not limited to constant-current charge voltage-time curve,
Constant-current discharge voltage-time curve, constant-voltage charge current versus time curve, constant-current charge current versus time curve etc..
Refering to fig. 1 and Fig. 2, in a preferred embodiment, step A includes being carried out to the part charge and discharge electrical measurement data curve
Pretreatment, the pretreatment include the following steps:
Step A1, it carries out fast Fourier (FFT) to the multiple discrete sampling point data to convert, time domain data is converted
To frequency domain;
Step A2, windowing process is carried out to frequency domain data and carries out smothing filtering;
Step A3, zero padding interpolation is carried out to filtered frequency domain data, makes its sample rate and the Global center reference number
Data sampling rate according to curve and the stepping center reference data and curves is consistent;
Step A4, by the transform frequency domain data after the smothing filtering and the interpolation to time domain.
In a preferred embodiment, in step B, by finding the part charge and discharge electrical measurement in preset sliding window
The least absolute value for measuring data and curves and the difference of the Global center reference data curve, to determine the part charge and discharge electrical measurement
Measure the starting point of the relatively described Global center reference data curve of data and curves.
In a further preferred embodiment, step B includes:
Step B1, using the set point on the Global center reference data curve as starting point, correspond to sampling sequence number w,
The one piece of data of setting length is taken, it is poor with the corresponding measurement data subtraction calculations on the part charge and discharge electrical measurement data curve
Value, the difference after subtracting each other is added up and calculates the absolute value of accumulated value, wherein w is initialized as 1;
Step B2, by the absolute value of the accumulated value calculated in step B1 and the difference minimum value e_ for being initially set maximal number
Min compares, if it is less than e_min, then using e_min as the absolute value of the accumulated value calculated in step B1;
Step B3, the starting point is updated, w plus 1 if w is more than the size of the sliding window, and terminates, otherwise
Return to step B1.
In an advantageous embodiment, step C includes the following steps:
Step C11. is bent with a plurality of stepping center reference data respectively by the part charge and discharge electrical measurement data curve
Each in line is matched, and determines the stepping center reference number with the part charge and discharge electrical measurement data curve difference value minimum
According to curve;
Step C12. uses the stepping center reference data that step C11 is determined will according to the starting point as curve of approximation
The charge and discharge electrical measurement data curve of remainder described in polishing, rebuilds complete charge and discharge electrical measurement number after the curve of approximation translation
According to curve.
In another preferred embodiment, step C includes the following steps:
Step C21. is gone to approach the part charge and discharge electrical measurement data song with a plurality of stepping center reference data and curves
Line, determination make the optimal coefficient that the quadratic sum of approximate error minimizes;
Step C22. is established approximate bent according to the optimal coefficient using a plurality of stepping center reference data and curves
Line, the charge and discharge electrical measurement data curve of remainder described in polishing, rebuilds complete charge and discharge after the curve of approximation is translated
Measurement data curve.
In a further preferred embodiment, in step C21,
If a plurality of stepping center reference data and curves are L curve, if every stepping center reference data and curves have N
A discrete sampling point data,
Including the part charge and discharge electrical measurement data of N number of discrete sampling point data is with vector
Completely charge and discharge data are M of l stepping center reference data and curves with vector
rl=[rl(1) … rl(N-1)rl(N) … rl(M)]T
The matrix of L × M of L stepping reference data composition is
L stepping reference data is respectively with pstdFor starting point, the matrix of the L × N for the data composition that length is N is
If undetermined coefficient vector is h=[h(1) … h(L)]T, least square approximation is carried out with following formula
Optimal coefficient h is solved to
After determining optimal coefficient, the approximating curve rebuild using L stepping center reference data and curves is
In a preferred embodiment, the Global center reference data curve and a plurality of stepping center reference data and curves
It is preset according to the material of battery and production technology.
Refering to Fig. 8, in another embodiment, a kind of computer readable storage medium, such as memory are stored thereon with meter
Calculation machine program, when the computer program is executed by processor realize according to aforementioned any embodiment based on battery portion
The method of the complete charge and discharge data of charge and discharge data reconstruction.
Refering to Fig. 8, in another embodiment, one kind being based on the complete charge and discharge data of battery portion charge and discharge data reconstruction
System, including:
Processor;And
Memory, the executable instruction for storing the processor;
Wherein, the processor is configured to execute according to aforementioned any embodiment institute via the executable instruction is executed
The method based on the complete charge and discharge data of battery portion charge and discharge data reconstruction stated.
In order to shorten the time of the techniques such as partial volume, stepping, production efficiency is improved, the embodiment of the present invention proposes above-mentioned base
In the method and system of the complete charge and discharge data of part charge and discharge data reconstruction.The present invention is based on preset reference data songs
Line need to only utilize charge and discharge electrical measurement data in part high-precision can rebuild complete charge and discharge data and curves.The system is applied to
When partial volume technique, the lithium cell of band part electricity need to only be charged to 100% state-of-charge, subsequent section electric discharge is wanted up to manufacture
It asks, you can using the complete charge and discharge data and curves of data reconstruction of part charging and partial discharge, to complete subsequent charge and discharge
Capacitance estimates and battery stepping, eliminates the cycle of intermediate primary complete charge and discharge, shortens for 50% or more time, significantly
Improve production efficiency.
The feature and its advantage of the specific embodiment of the invention elaborated further below.
1, part charge and discharge electrical measurement data is pre-processed
If charge and discharge electrical measurement data in part includes K discrete sampling point { sk, k=0,2 ... K-1
1) it carries out fast Fourier (FFT) to K measurement data to convert, time-domain signal is transformed into frequency domain
Wherein k is discrete frequency point serial number.
2) windowing process is carried out to frequency-region signal and carries out smothing filtering, be such as superimposed raised cosine (raised cosine) window
The frequency spectrum of raised cosine window, β are cutoff frequency
Rate, T are amplitude.
3) sample rate is improved to filtered frequency-region signal zero padding interpolation, it is consistent with reference data sample rate
4) time domain is transformed into the frequency domain data after smothing filtering and interpolation
2, to the reconstruction of part charge and discharge cycles measurement data
The reconstruction of part charge and discharge electrical measurement data needs to be based on several preset reference data curves, including the model
The stepping center reference data and curves of each stepping of Global center reference data curve and lithium battery of lithium battery.These reference numbers
It is closely related according to the material and leading portion production technology of curve and lithium battery, it is set in advance according to material and technique in process of production
It is fixed.It can it will be understood by those skilled in the art that presetting battery charging and discharging reference data curve according to the material of battery and technique
To be realized using ordinary skill in the art means.
Starting point of the 2.1 estimation part charge and discharge measurement data curves with respect to reference data curve
The nuance of state-of-charge when starting due to charging and discharging lithium battery first has to determine measurement data relative to the overall situation
The starting point of center reference data and curves.The present invention proposes that the method for finding minimal difference in one section of sliding window calculates and measures
Curve starting point.By the minimum for finding measurement data curve and the difference of Global center reference data curve in sliding window
Absolute value determines starting point of the experiment curv relative to reference curve.The specific method is as follows:
Step S1., using set point as starting point, corresponds to sampling sequence number w, takes setting in Global center reference data curve
Difference after subtracting each other is added up with measurement data subtraction calculations difference and calculates the absolute of accumulated value by the one piece of data of length N
Value.
Step S2. is by the absolute value of the accumulated value calculated in step S1 and the difference minimum value for being initially set maximal number
eminCompare, if it is less than emin, then by eminAs the absolute value of the accumulated value calculated in step S1, and record w.
W is added 1 by step S3., if w is more than the size W of sliding window, and is terminated, otherwise return to step S1.
Parameter declaration:
w:Sampling corresponding to Global center reference data curve for the starting point of the data segment compared with measurement data
Serial number is initialized as 1;
N:Measurement data sampled point number;
N-th of sampled point of measurement data;
r0(n):N-th of sampled point of Global center reference data curve;
W:Sliding window size;
emin:The least absolute value for recording measurement data curve and the difference of Global center reference data curve, is initialized as
Maximal number;
pstd:Minimal difference starting point is recorded in the sampled point serial number of Global center reference data curve, is initialized as 0.
The pseudocode of algorithm is as follows:
2.2. complete measurement data and curves are rebuild
1) linear to rebuild
In one embodiment, the reconstruction of complete measurement data is completed by the method linearly rebuild.Linear differential weight
Building is matched respectively with each stepping center reference data and curves using part charge and discharge electrical measurement data curve, finds difference most
Approximation of the small stepping center reference data and curves as measurement data, to rebuild complete measurement data.
It is as follows:
Step 1. matches part charge and discharge electrical measurement data curve with each stepping center reference data and curves respectively, seeks
Look for the stepping center reference data and curves with part charge and discharge electrical measurement data difference minimum;
Step 2. uses the stepping center reference data chosen as curve of approximation, and polishing after approximate curvilinear translation is measured number
According to curve remainder, to rebuild complete measurement data curve.
Parameter declaration:
L:Battery stepping number;
p:It is directed toward with part charge and discharge electrical measurement data curve difference value and minimum stepping center reference data, p initialization
For maximum;
Dmin:Record the minimal difference of part charge and discharge measurement data curve and stepping center reference data and curves and,
DminIt is initialized as maximum.
Pseudo-code of the algorithm is as follows:
Complete measurement data after reconstruction are represented by
WhereinTo use m-th of sampled point of the partial data curve after the reconstruction of linear reconstruction method,For portion
Divide m-th of sampled point of charge and discharge electrical measurement data, rp(m) it is and part charge and discharge electrical measurement data curve difference value and minimum point
M-th of sampled point of shelves center reference data.
2) least squares are rebuild
In another embodiment, the reconstruction of complete measurement data is completed by least square method for reconstructing.Minimum two
It is that the song of measurement data is rebuild with the criterion of least square using existing stepping center reference data as substrate to multiply reconstruction
Line.
Least square reconstruction is gone to approach the part charge and discharge electrical measurement number got with L stepping center reference data and curves
According to determining optimal coefficient with the quadratic sum of the approximate error of part charge and discharge electrical measurement data by minimizing.Determine most major clique
After number, approximating curve, polishing measurement data curve after approximating curve is translated are rebuild using L stepping center reference data and curves
Remainder is to rebuild complete measurement data and curves.
It is as follows:
If including the part charge and discharge electrical measurement data of N number of sampled point is expressed as with vector
If M complete charge and discharge data of l stepping center reference data and curves are expressed as with vector
rl=[rl(1) … rl(N-1)rl(N) … rl(M)]T
The matrix of L × M of L stepping reference data composition is represented by
L stepping reference data is respectively with pstdMatrix for starting point, the L × N for the data composition that length is N is represented by
If undetermined coefficient vector is h=[h(1) … h(L)]T, the mathematic(al) representation of least square approximation is as follows
Optimal coefficient h can be solved to
After determining optimal coefficient, the approximating curve rebuild using L stepping center reference data and curves is represented by
Polishing measurement data curve remainder is to rebuild complete measurement data and curves after approximating curve is translated, after reconstruction
Complete measurement data be represented by
WhereinTo use m-th of sampled point of the partial data curve after the reconstruction of least square reconstruction method,For m-th of sampled point of part charge and discharge electrical measurement data,The approximating curve established for least square method
M-th of sampled point.
Test result
1200 battery samples use same specification lithium cell, 5 grades of battery point, each stepping center reference curve such as Fig. 4 institutes
Show.
Component capacitive equipment uses negative pressure forming probe bed and 5V-300A power cabinets, is carried out at the appointed time with 60A electric currents
Constant-current discharge measures complete constant-current discharge voltage curve.Using the first half discharge curve of 1200 batteries, this hair is used
Bright method carries out Curve Reconstruction, by reconstructed results and experiment curv results contrast, calculates the evaluated error of reconstruction, and utilize weight
Curve after building carries out stepping.
In terms of result, the method for reconstructing in the present invention has accurately rebuild complete discharge data using partial discharge data,
When least square reconstruction method higher using complexity, mean square error reaches 5.6E-7, and no battery is by wrong stepping.It is relatively common
Stepping technique, entire stepping cycle time 50% or more, substantially increases production efficiency, has high practical value.
1 mean square error of table and wrong sorting rate
|
Mean square error |
Mistake sorting rate |
Linear reconstruction method |
1.4E-4 |
0.2% |
Minimum two reconstruction methods |
5.6E-7 |
0% |
It will be understood by those skilled in the art that various illustrative components, blocks, module, circuit and algorithm in conjunction with described in present disclosure
Step may be implemented as electronic hardware, computer software or both combination.In order to clearly demonstrate this of hardware and software
Interchangeability, it is above that various Illustrative components, module, module, circuit and step are generally described in terms of its function.
This function is to depend on applying specific application over the whole system and design constraint by hardware or software implementation.Technology
Personnel can be implemented in various ways the function for each specific application, but this implementation decision is not necessarily to be construed as partially
From the scope of the present invention.Those skilled in the art also will readily appreciate that, component described herein, method or interaction
Sequence or combination are only example, and component, method or the interaction of the various aspects of present disclosure can be to be different from the present invention
Those of described mode is combined or implements.
Functional block and module may include processor, electronic equipment, hardware device, electronic building brick, logic circuit, storage
Device, software code, firmware code etc. or its arbitrary combination.With it is aforementioned consistent, can use general processor, digital signal
Processor (DSP), application-specific integrated circuit (ASIC), field programmable gate array (FPGA) or other programmable logic devices, solely
The various explanations in conjunction with described in present disclosure are realized or are executed in vertical door or transistor logic, independent hardware component or its arbitrary combination
Property logical block, module and circuit for executing function of the present invention.General processor can be microprocessor, but replace
Dai Di, processor can be any conventional processors, controller, microcontroller or state machine.Processor can also be implemented as
The combination of computing device, for example, the combination of DSP and microprocessor, multi-microprocessor, in conjunction with DSP core one or more it is micro-
Processor or any other such configuration.
In conjunction with the step of present disclosure the method or algorithm can directly with hardware, by processor execute software module,
Or both combination mode embody.Software module may reside within RAM memory, flash memory, ROM memory, eprom memory,
The storage medium of eeprom memory, register, hard disk, removable disk, any other formats of CD-ROM or known in the art
In.Storage medium is connected to processor so that letter can be written from read information and to storage medium in processor
Breath.In alternative solution, storage medium can be integrated into processor.Pocessor and storage media may reside in ASIC.
ASIC may reside in user terminal, base station, sensor or any other communication equipment.In alternative solution, processor and
Storage medium can be used as stand-alone assembly resident in the user terminal.
In one or more exemplary designs, the function can be with hardware, software, firmware, or any combination thereof
Implement.If implemented with software mode, it is computer-readable that the function can be used as one or more instructions or code to be stored in
It transmits on medium or on a computer-readable medium.Computer-readable medium includes computer storage media and communication media,
Including any medium for helping for computer program to be transmitted to another place from a place.Computer readable storage medium
It can be any usable medium that can be accessed by general or specialized computer.As an example, not a limit, this computer can
It may include that RAM, ROM, EEPROM, CD-ROM or other disk storage, magnetic disk storage or other magnetic storages are set to read medium
Any other medium that is standby or can be used for it is expected with instruction or data structure format carrying or storage program code devices, can
With by general or specialized computer or the access of general or specialized processor.In addition, connection be properly termed as computer can
Read medium.For example, if using coaxial cable, fiber optic cables, twisted-pair feeder or Digital Subscriber Line (DSL) from website, server or
Other remote source softwares, then coaxial cable, fiber optic cables, twisted-pair feeder or DSL are comprised in the definition of medium.The present invention
The disk and CD used includes CD (CD), laser disk, CD, digital versatile disc (DVD), floppy disk and Blu-ray disc, wherein
Disk usually magnetically reproduce data, and CD using laser optics reproduce data.Combination of the above also should by comprising
In the range of computer-readable medium.
Although embodiments of the present invention and their advantages have been described in detail, it is to be understood that, is not departing from by appended power
In the case that profit requires the embodiment spirit and scope limited, various changes, replacement and change can be carried out herein.This
Outside, the scope of the present invention is not limited to process, machine, manufacture, material composition, means, method and step described in the specification
Rapid specific embodiment.It is executed and corresponding embodiment described herein one ordinarily skilled in the art will readily appreciate that can utilize
Essentially identical function or obtain with embodiment described herein essentially identical result presently, there are or to be developed later above-mentioned
Disclosure, process, machine, manufacture, material composition, means, method or step.Therefore, appended claims be intended to by these processes,
Machine, manufacture, material composition, means, method or step are embraced within their scope.