CN110161421A  A kind of method of battery impedance within the scope of online reorganization setpoint frequency  Google Patents
A kind of method of battery impedance within the scope of online reorganization setpoint frequency Download PDFInfo
 Publication number
 CN110161421A CN110161421A CN201910428872.0A CN201910428872A CN110161421A CN 110161421 A CN110161421 A CN 110161421A CN 201910428872 A CN201910428872 A CN 201910428872A CN 110161421 A CN110161421 A CN 110161421A
 Authority
 CN
 China
 Prior art keywords
 voltage
 impedance
 frequency
 scope
 method
 Prior art date
Links
 HBBGRARXTFLTSGUHFFFAOYSAN Lithium Ion Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnID4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTQwLjM2NicgeT0nMTU4LjI1JyBzdHlsZT0nZm9udC1zaXplOjE1cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojMDAwMDAwJyA+PHRzcGFuPkxpPC90c3Bhbj48dHNwYW4gc3R5bGU9J2Jhc2VsaW5lLXNoaWZ0OnN1cGVyO2ZvbnQtc2l6ZToxMS4yNXB4Oyc+KzwvdHNwYW4+PHRzcGFuPjwvdHNwYW4+PC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyA+CjwhLS0gRU5EIE9GIEhFQURFUiAtLT4KPHJlY3Qgc3R5bGU9J29wYWNpdHk6MS4wO2ZpbGw6I0ZGRkZGRjtzdHJva2U6bm9uZScgd2lkdGg9Jzg1JyBoZWlnaHQ9Jzg1JyB4PScwJyB5PScwJz4gPC9yZWN0Pgo8dGV4dCB4PSczMi4zNjYxJyB5PSc1MC4yNScgc3R5bGU9J2ZvbnQtc2l6ZToxNHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzAwMDAwMCcgPjx0c3Bhbj5MaTwvdHNwYW4+PHRzcGFuIHN0eWxlPSdiYXNlbGluZS1zaGlmdDpzdXBlcjtmb250LXNpemU6MTAuNXB4Oyc+KzwvdHNwYW4+PHRzcGFuPjwvdHNwYW4+PC90ZXh0Pgo8L3N2Zz4K [Li+] HBBGRARXTFLTSGUHFFFAOYSAN 0 claims description 4
 238000000034 methods Methods 0 description 12
 238000004458 analytical methods Methods 0 description 2
 238000004364 calculation methods Methods 0 description 2
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
 G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
 G01R31/367—Software therefor, e.g. for battery testing using modelling or lookup tables

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
 G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
 G01R31/385—Arrangements for measuring battery or accumulator variables
 G01R31/387—Determining amperehour charge capacity or SoC

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
 G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
 G01R31/389—Measuring internal impedance, internal conductance or related variables
Abstract
Description
Technical field
The present invention relates to battery impedance testing fields, hinder more particularly, to battery within the scope of a kind of online reorganization setpoint frequency Anti method.
Background technique
Currently, lithium ion battery is big with its capacity density, have extended cycle life, the advantages that charging rate is fast obtained it is extensive Using.But although lithium ion battery possesses plurality of advantages, due to itself belonging to very complicated electrochemical systems, Reaction is made of many complicated electrode processes, this brings great difficulty to its characteristic of research and application.
Electrochemical impedance spectroscopy (Electrochemical Impedance Spectroscopy, EIS) has broadband model The abundant impedance information enclosed can reflect each electrode process of battery, can also characterize temperature, the stateofcharge (State of battery Of Charge, SOC) and the characteristics such as health status (State of Health, SOH), therefore Applied Electrochemistry impedance spectrum is to electricity Pond characteristic, which carries out research, becomes hot spot.
The measuring principle of electrochemical impedance spectroscopy is to work as battery to be under equilibrium state, applies small size simple alternating current excitation letter Number, the impedance under current frequency is calculated by calculating acexcited response under same frequency and simple alternating current excitation.Due to EIS Measuring instrument is costly, and in practical applications without acexcited source, thus traditional EIS measurement usually laboratory into Row.In order to solve this problem, many online researchs for obtaining electrochemical impedance are dedicated to the electrochemistry that design has alternating current source Impedance measurement carmounted device, but this method not only will increase cost, but also will increase battery management system design complexities.This There are some scholars outside using timefrequency convert method, using the voltage and current data computing impedance in battery time domain, but this side Method can be only applied to specific criteria operating condition, however the operating condition of battery is in dynamic change on real vehicle, and institute is in this way There is also certain limitations.
Just because of current impedance Online Measuring Method, there are many limitations, and there is an urgent need to propose that a kind of improvement project is come It is online to obtain battery impedance.
Summary of the invention
It is an object of the present invention to overcome the abovementioned drawbacks of the prior art and provide a kind of online reorganizations to set The method of battery impedance in frequency range.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of method of battery impedance within the scope of online reorganization setpoint frequency, to measure the resistance of vehiclemounted lithium ion battery It is anti, comprising the following steps:
1) voltage data and current data and setpoint frequency range in time domain are obtained online；
2) adding window correction and Fast Fourier Transform (FFT) are carried out to voltage and current data, obtains the voltage and current in frequency domain Value, and choose target frequency section；
3) available voltage and current data are selected using filter criteria；
4) available voltage and current data are fitted to obtain the expression formula of voltage, electric current using least square method, And impedance is obtained accordingly.
The frequency range is set as 0.016Hz.
It is power function variation relation between voltage and current and frequency within the scope of setpoint frequency in the step 4).
The objective function that the voltage and current uses least square method to be fitted are as follows:
Log  U =a_{1}·logf+b_{1}
Log  I =a_{2}·logf+b_{2}
Wherein,  U  it is voltage magnitude,  I  it is current amplitude, f is frequency, a_{1}, a_{2}, b_{1}, b_{2}For fitting constant.
Expression formula in the step 4), between impedance and frequency are as follows:
In the step 2), by adding Hanning window to be corrected.
In the step 2), the target frequency segment limit of selection is 050Hz.
In the step 3), available voltage and electricity are chosen by the fit equation index of standard deviation and impedance magnitude Flow data is then determined as available voltage when standard deviation is less than the fit equation index of given threshold and impedance magnitude less than 0 And current data.
Compared with prior art, the invention has the following advantages that
1) Fast Fourier Transform (FFT) method is utilized, realization obtains battery impedance in the dynamic case, overcomes traditional resistance Anti spectrum can only measure the shortcomings that battery impedance under the offline state of laboratory.
2) present invention carries out windowing process to the voltage and current data of acquisition to reduce spectral leakage, improves impedance width It is worth computational accuracy.
3) the invention proposes screening index, ensure that measurement data availability in computing impedance.
4) least square method approximating method is used, calculating process is simple, and easy to accomplish, calculating speed is fast.
Detailed description of the invention
Fig. 1 is the schematic diagram of 0.01Hz6Hz impedance magnitude under different SOC.
Fig. 2A is the voltage data measured under time domain；Fig. 2 B is the current data measured under time domain.
Fig. 3 A is the voltage amplitudefrequency characteristic schematic diagram in 050Hz frequency range after Fast Fourier Transform (FFT)；Fig. 3 B is electricity Press phasefrequency characteristic schematic diagram；Fig. 3 C is the electric current amplitudefrequency characteristic signal in 050Hz frequency range after Fast Fourier Transform (FFT) Figure；Fig. 3 D is voltage phasefrequency characteristic schematic diagram.
Fig. 4 A is voltage magnitude characteristic and the linear fit result schematic diagram based on least square method；Fig. 4 B is electric current width It is worth characteristic and the linear fit result schematic diagram based on least square method.
Fig. 5 is the method for the present invention flow chart.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
Embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described.Obviously, described embodiment is A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
Firstly, it is necessary to first analyze the impedance magnitude changing rule in target frequency section.Fig. 1 is that impedance magnitude exists Variation tendency in 0.01Hz6Hz frequency range.It can be seen from the figure that impedance magnitude is about within the scope of 0.01Hz6Hz As the raising of frequency quickly reduces in 0.5Hz, after 1Hz6Hz with frequency increase and reduced trend gradually tends to be flat It is slow.From the direct observation of figure it is found that in 0.01Hz6Hz frequency range, power function variation rule are presented in impedance magnitude and frequency Rule.
For sampling thheorem and the indoor equipment sampling precision problem of experiment, the present invention is to standard work each in battery time domain The sampling period of voltage and current data is set as 0.01s under condition, i.e. sample frequency is 100Hz, through fast Fourier (Fast Fourier Transformation, FFT) transformation after the available valuable spectrum information in 50Hz.
Selection for battery floor data in time domain, the present invention choose an example and are illustrated.Fig. 2A and Fig. 2 B institute Be shown as the unified Lightduty Vehicle test procedure in the whole world (Worldwideharmonized Light Vehicles Test Cycle, WLTC) the voltage signal and current signal under operating condition.After selected characteristic frequency band, need to the voltage and current data in time domain Carry out Fast Fourier Transform (FFT).Fourier transformation is a kind of important method of digital processing field, can be by signal from original Beginning time domain is converted to the representation of frequency domain.Any signal continuously measured can be decomposed into the sine wave letter of different frequency Number, to obtain amplitude and the phase angle of sine wave signal, i.e. spectrum information in a frequency domain.Actually should during, be used for Expression and the discrete Fourier transform (Discrete Fourier Transform, DFT) of analysis discrete time signal are more often With.DFT's is defined as:
Wherein, X (k) indicates the data after DFT transform, and x (n) is the analog signal of sampling.
Fast Fourier Transform (FFT) can more rapidly process digital signal, be more widely applied.FFT method is to DFT The improvement of method, this method can greatly reduce multiplication number, improve calculating speed.Based on FFT method, voltage and electricity are obtained The amplitudefrequency characteristic and phasefrequency characteristic of stream are as shown in Fig. 3 A Fig. 3 D.Wherein, Fig. 3 A is voltage amplitudefrequency characteristic；Fig. 3 B is voltage phase frequency Characteristic；Fig. 3 C is electric current amplitudefrequency characteristic；Fig. 3 D is electric current phasefrequency characteristic.
As can be seen that amplitude with the variation of frequency meets Bending influence from Fig. 3 A, Fig. 3 C, this conclusion with The impedance magnitude that experiment measures in 0.01Hz6Hz frequency range is consistent with the variation relation of frequency.The present invention uses Function Fitting Method impedance magnitude online calculation method is designed.
The relationship of voltage, electric current and frequency is respectively as shown in formula (2) and formula (3):
Wherein,  U  it is voltage magnitude,  I  it is current amplitude, f indicates frequency, a_{1}, a_{2},For constant.
In order to simplify fit equation (2) and formula (3), consider to carry out under loglog coordinate using linear least square method Linear fit.When using linear least square fitting method, the objective function of voltage and current fitting is respectively as follows:
Log  U =a_{1}·logf+b_{1}#(4)
Log  I =a_{2}·logf+b_{2}#(5)
Abovementioned equation difference is as shown in Figure 4 A and 4 B shown in FIG..Wherein, Fig. 4 A is the fitting of voltage magnitude logarithm and frequency logarithm As a result, Fig. 4 B is the fitting result of current amplitude logarithm and frequency logarithm.Oblique line in figure is fitting a straight line, and scatterplot is FFT change Data point after changing.
Fig. 5 is finally established battery impedance amplitude online calculation method flow chart.In step 501, the electricity in time domain Pressure and current data can directly obtain in practical applications.In step 503, in order to reduce the influence of spectrum leakage, to adopting Collect datasignal adding window.The premise of FFT assumes that time signal is that the period is unlimited, but one section can only be intercepted in actual operation Signal, which results in spectrum leakages.Therefore to solve this problem it is necessary to carry out windowing process to data.Following Step 505 in, FFT transform is carried out to the signal after adding window, the voltage and current value in available frequency domain.
It, can since the direct current component being superimposed in voltage in time domain, current data can be separated by FFT transform It is rejected with the influence by removal DC component by bias direct current to impedance magnitude, as illustrated at step 507.
According to Shannon's sampling theorem, sample frequency is at least twice of signal highest frequency, otherwise will lead to the height of script Frequency signal is sampled low frequency signal, the aliasing of frequency spectrum occurs.Generally in practical applications, the sample frequency of BMS can be with Reach 100Hz, so range of target frequencies is 050Hz, as indicated at 509.
Found in actual moving process, not all floor data after FFT transform and least square fitting, The impedance magnitude being calculated all meets required precision, i.e. error is maintained within the scope of certain error.Two are primarily present to ask Topic, first is that the impedance magnitude numerical problem calculated, second is that impedance magnitude changing rule Trend Problem.Therefore, the invention proposes Screening conditions, to remove not available data in advance.
For first problem, selection criteria difference is as Rule of judgment.By analysis, not available voltage and current data Have the characteristics that distribution diverging and the order of magnitude is too small.In statistics, standard deviation is a kind of for quantifying the change of one group of data value The measure of change amount.The smaller expression data point of standard deviation tends to the mean value close to set, and standard deviation is larger, indicates data Point is distributed in bigger value range, be easy to cause the inaccuracy of fitting result.It can be with by the comparison of different data index It was found that the standard deviation of voltage and current should limit in a certain range, so as to reject unavailable data.Step 511 embodies This process.
After having selected available data, logarithm can be taken simultaneously to voltage magnitude and current amplitude.In step 515, Function Fitting is carried out using least square method, obtains the functional relation expression between voltage and frequency and electric current and frequency Formula.In step 517, result is obtained to fitting and carries out antilogarithm operation.May finally obtain voltage magnitude  U  with frequency f it Between relationship and current amplitude  I  the relationship between frequency.Based on both sides relation, in step 519, can be calculated Relationship between impedance and frequency:
For Second Problem, according to the available index of Fig. 1 less than zero, therefore second screening conditions can be proposed, That is parameter (a_{1}a_{2}) less than zero.If (a_{1}a_{2}) it is greater than zero, illustrate that the trend of matched curve is incorrect, therefore such Data point needs are removed.Detailed process is as shown in step 521.In conclusion the invention proposes two main filter criterias: 1. the standard deviation of voltage data and current data；2. the index of the fit equation of impedance magnitude.
In addition, systematic error must be taken into consideration in measurement process.In actual measurement, impedance magnitude due to systematic error all Greater than offline measured value, but curve can preferably concentrate in certain section, therefore according to this feature by calculated width Value curve translates downwards a distance, corrects the error.
Claims (8)
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201910428872.0A CN110161421A (en)  20190522  20190522  A kind of method of battery impedance within the scope of online reorganization setpoint frequency 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201910428872.0A CN110161421A (en)  20190522  20190522  A kind of method of battery impedance within the scope of online reorganization setpoint frequency 
Publications (1)
Publication Number  Publication Date 

CN110161421A true CN110161421A (en)  20190823 
Family
ID=67631935
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201910428872.0A CN110161421A (en)  20190522  20190522  A kind of method of battery impedance within the scope of online reorganization setpoint frequency 
Country Status (1)
Country  Link 

CN (1)  CN110161421A (en) 
Citations (8)
Publication number  Priority date  Publication date  Assignee  Title 

JP2006250905A (en) *  20050314  20060921  Fuji Heavy Ind Ltd  Battery management system 
JP4030217B2 (en) *  19990312  20080109  トヨタ自動車株式会社  Abnormality determination device and abnormality determination method for battery pack 
CN104849675A (en) *  20150617  20150819  哈尔滨工业大学  Method for obtaining electrochemical and thermal coupling models of lithium ion battery 
CN105223508A (en) *  20150714  20160106  上海空间电源研究所  Inside lithium ion cell performance state lossless detection method 
CN105891715A (en) *  20141212  20160824  广西大学  Lithium ion battery health state estimation method 
CN106371029A (en) *  20160826  20170201  中国电力科学研究院  Lithium battery AC impedance spectrum online synchronous test method and device 
CN108761341A (en) *  20180601  20181106  哈尔滨工业大学  A kind of lithium ion battery battery chemical modeling parameter acquisition methods 
CN109143074A (en) *  20180628  20190104  中国科学院光电研究院  A kind of power battery model parameter identification method and system 

2019
 20190522 CN CN201910428872.0A patent/CN110161421A/en active Search and Examination
Patent Citations (8)
Publication number  Priority date  Publication date  Assignee  Title 

JP4030217B2 (en) *  19990312  20080109  トヨタ自動車株式会社  Abnormality determination device and abnormality determination method for battery pack 
JP2006250905A (en) *  20050314  20060921  Fuji Heavy Ind Ltd  Battery management system 
CN105891715A (en) *  20141212  20160824  广西大学  Lithium ion battery health state estimation method 
CN104849675A (en) *  20150617  20150819  哈尔滨工业大学  Method for obtaining electrochemical and thermal coupling models of lithium ion battery 
CN105223508A (en) *  20150714  20160106  上海空间电源研究所  Inside lithium ion cell performance state lossless detection method 
CN106371029A (en) *  20160826  20170201  中国电力科学研究院  Lithium battery AC impedance spectrum online synchronous test method and device 
CN108761341A (en) *  20180601  20181106  哈尔滨工业大学  A kind of lithium ion battery battery chemical modeling parameter acquisition methods 
CN109143074A (en) *  20180628  20190104  中国科学院光电研究院  A kind of power battery model parameter identification method and system 
NonPatent Citations (2)
Title 

张珺涵 等: "基于时频分析的锂离子电池阻抗计算方法", 《电池》 * 
徐睿: "锂离子电池劣化程度与阻抗模型相关性研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * 
Similar Documents
Publication  Publication Date  Title 

Kang et al.  A new neural network model for the stateofcharge estimation in the battery degradation process  
CN103257286B (en)  A kind of charging electric vehicle facility automatic test approach and system  
Li et al.  A quick online state of health estimation method for Liion battery with incremental capacity curves processed by Gaussian filter  
Xiong et al.  Towards a smarter battery management system: A critical review on battery state of health monitoring methods  
Howey et al.  Online measurement of battery impedance using motor controller excitation  
KR101331325B1 (en)  Device for calculating impedances of battery cell and battery impedance measuring system  
Karden et al.  A method for measurement and interpretation of impedance spectra for industrial batteries  
CN104502859B (en)  The detection of battery charge amount and cell health state, diagnostic method  
Ran et al.  Prediction of state of charge of lithiumion rechargeable battery with electrochemical impedance spectroscopy theory  
Blanke et al.  Impedance measurements on lead–acid batteries for stateofcharge, stateofhealth and cranking capability prognosis in electric and hybrid electric vehicles  
CN100394208C (en)  System and method for measuring fuel cell voltage and high frequency resistance  
CN102305891B (en)  Online monitoring method of lowfrequency oscillation of power system  
JP2014074682A (en)  Estimation apparatus and estimation method of parameter or the like of battery  
US20110148424A1 (en)  Apparatus for estimating battery state of health  
DE69927058T2 (en)  Apparatus and method for detecting characterizing parameters of load stores  
US8421486B2 (en)  Oildegradation detecting apparatus  
CN101701986B (en)  System for pretesting and diagnosing electro magnetic interference of electronic equipment and method thereof  
Cottis et al.  Measures for the detection of localized corrosion with electrochemical noise  
CN101622547B (en)  Device and method for detecting abnormality of electric storage device  
JP5354480B2 (en)  Fuel cell evaluation system  
CN101629925B (en)  Method and device for measuring the conductivity of a pure or ultrarapture liquid  
JP3162030B2 (en)  Battery capacity measuring method and battery capacity measuring device using voltage response signal of pulse current  
CN105388422B (en)  Secondary battery capacity measuring system and secondary battery capacity measurement method  
WO2013185405A1 (en)  Power transformer iron core looseness fault diagnosis and preliminary fault location system and method  
US6208147B1 (en)  Method of and apparatus for measuring battery capacity by impedance spectrum analysis 
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 