Method of detecting system function by measuring frequency response

• US 7,395,163 B1
• Filed: 07/05/2007
• Issued: 07/01/2008
• Est. Priority Date: 12/21/2004
• Status: Active Grant

• Current Assignee: Battelle Energy Alliance LLC (Battelle Memorial Institute)
• Original Assignee: Montana Tech of the University of Montana (University of Montana System)
• Inventors: Morrison, John L., Morrison, William H.
• Primary Examiner(s): WACHSMAN, HAL D
• Application Number:
  US11/825,629
• Time in Patent Office: 362 Days
• Field of Search: 702/117, 702/75, 702/79
• US Class Current: 702/75
• CPC Class Codes: G01R 31/389 Measuring internal impedance, internal conductance or related variables

Abstract

Real time battery impedance spectrum is acquired using one time record, Compensated Synchronous Detection (CSD). This parallel method enables battery diagnostics. The excitation current to a test battery is a sum of equal amplitude sin waves of a few frequencies spread over range of interest. The time profile of this signal has duration that is a few periods of the lowest frequency. The voltage response of the battery, average deleted, is the impedance of the battery in the time domain. Since the excitation frequencies are known, synchronous detection processes the time record and each component, both magnitude and phase, is obtained. For compensation, the components, except the one of interest, are reassembled in the time domain. The resulting signal is subtracted from the original signal and the component of interest is synchronously detected. This process is repeated for each component.

Claims

1. A method for detecting system function of a unit under test by measuring frequency response, the method comprising the steps of:

(a) selecting a number of frequencies over which the function of the unit under test will be tested;

(b) assembling an excitation time record that is the sum of the sinusoids of the frequencies and a duration of greater than or equal to a period of the lowest selected frequency;

(c) conditioning the excitation time record to be compatible with the unit under test;

(d) exciting the unit under test with the excitation time record and simultaneously capturing a response time record with a data acquisition system;

(e) processing the response time record using the following synchronous detection equations to obtain estimated frequency components of magnitude and phase for one of the selected frequencies;

2. The method of claim 1, wherein said error is minimized using the following equations which create a corrected time record and a compensated time record:

$\begin{array}{cc}{f}_{\mathrm{RK}}<br><br>\left[j\right]=\underset{p=1,p\ne <br><br>K}{\overset{M}{\sum <br><br>}}<br><br><br><br>\left(F_p<br><br>\sin <br><br>\left({\mathrm{\omega <br><br>}}_p<br><br>\left(j-1\right)<br><br>\mathrm{\Delta <br><br>}<br><br><br><br>t\right)+F_{\mathrm{Qp}}<br><br>\cos <br><br>\left({\mathrm{\omega <br><br>}}_p<br><br>\left(j-1\right)<br><br>\mathrm{\Delta <br><br>}<br><br><br><br>t\right)\right);<br><br>j=R+1<br><br>\text{;<br><br>}<br><br>N& \left(11\right)\end{array}$
CfK_Out[j]=f_Out[j]−

f_RK[j]

Where;

f_Out is the original time record

f_RK is the correction time record

Cf_KOut is the compensated time record

F_p estimated in phase amplitude response at the p^th frequency

F_Qp estimated quadrature amplitude response at the p^th frequency

ω

_p radian frequency of the p^th sinusoid

Δ

t time step of data system

N number of points of the response time record

M number of different sinusoids of the excitation time record

R number of points of the response time record that may be discarded

K is the frequency index for the system function being detected.

3. The method of claim 1, wherein said error is minimized using a neural network.

4. The method of claim 1, wherein said unit under test is a battery.