Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
SYSTEM AND METHOD FOR STATE OF POWER ESTIMATION OF A BATTERY USING IMPEDANCE MEASUREMENTS
Document Type and Number:
WIPO Patent Application WO/2023/212699
Kind Code:
A9
Abstract:
A method is provided for pretraining a hyper model configured for use in predicting a state of power (SoP) of a vehicle battery. The method includes performing electrochemical impedance spectroscopy (EIS) scans on a plurality of batteries having a set of similar operating characteristics to the vehicle battery. The EIS scans are performed across various states of the vehicle battery. The method further includes fitting parameters of the hyper model by applying an optimization technique to results of the EIS scans. The hyper model includes a family of models that each define a voltage response of a respective cell from among a plurality of cells of the vehicle battery to a current profile over the various states of the vehicle battery.

Inventors:
TRAA JOHANNES (US)
TANOVIC OMER (US)
GULLAPALLI HEMTEJ (US)
Application Number:
PCT/US2023/066370
Publication Date:
March 28, 2024
Filing Date:
April 28, 2023
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
ANALOG DEVICES INC (US)
International Classes:
G01R31/367; B60L3/00; B60L3/12; B60L58/18; G01R31/3842; G01R31/389; G01R31/396
Attorney, Agent or Firm:
ARORA, Suneel et al. (US)
Download PDF:
Claims:
CLAIMS: 1. A method for pretraining a hyper model configured for use in predicting a state of power (SoP) of a vehicle battery, the method comprising: performing electrochemical impedance spectroscopy (EIS) scans on a plurality of batteries having a set of similar operating characteristics to the vehicle battery, the EIS scans performed across various states of the vehicle battery; and fitting parameters of the hyper model by applying an optimization technique to results of the EIS scans, the hyper model comprising a family of models that each define a voltage response of a respective cell from among a plurality of cells of the vehicle battery to a current profile over the various states of the vehicle battery. 2. The method according to claim 1, further comprising estimating the state of power (SoP) of the vehicle battery using the hyper model, the results of the EIS scans, operational constraints of the vehicle battery, and the various states of the vehicle battery used in the EIS scans. 3. The method according to claim 1, wherein the various states of the vehicle battery include at least some of: different temperature ranges of the vehicle battery; different states of charge (SoCs) of the vehicle battery; age of the vehicle battery; and a nature of a current load the vehicle battery is subjected to. 4. The method according to claim 1, wherein the hyper model is an equivalent circuit model (ECM) that maps various circuit elements such as resistors, capacitors, inductors, and Warburg impedance of the ECM to the various states of the vehicle battery used in the EIS scans. 5. The method according to claim 1, wherein the hyper model is a model reduced from a physics based model that maps physics representations of battery elements to the various states of the vehicle battery used in the EIS scans. 6. The method according to claim 1, wherein the hyper model is an adaptive filter inferred using a frequency response of an equivalent impedance of the vehicle battery learned using the EIS scans under the various states of the vehicle battery. 7. The method according to any of claims 1-6, wherein the fitting is performed offline. 8. The method according to any of claims 1-7, wherein a multidimensional look up table is used to map the hyper model, an expected battery state, the results of the EIS scans, and battery system current and voltage limitations into a corresponding SoP value. 9. The method according to claim 1, wherein a mapping from the EIS scans, the various states of the vehicle battery, and operational constraints of the vehicle battery to the state of power (SoP) of the battery is performed using a regression function learned under the various states of the vehicle battery. 10. The method according to claim 1, wherein the results of the EIS scans are used to infer a hyper model using an optimization technique selected from the group consisting of a gradient-based linear optimization method, a non-negative least squares (NNLS) method, and a convex optimization method. 11. The method according to claim 1, wherein the hyper model is a family of equivalent circuit models (ECMs), and the method further comprises performing a smart initialization of a parameter fitting method by setting a series resistor R in an R—RC ECM model having the series resistor R in series with one or more RC parallel sub-circuits to a smallest observed impedance value, setting model parameters to determined values and holding the determined values fixed while scanning over a range based on a value of the series resistor R to identify a determined value that minimizes an objective function. 12. A method for predicting a state of power (SoP) of a battery, the method comprising: performing a plurality of electrochemical impedance spectroscopy (EIS) scans on the battery prior to an initial use of the battery in a vehicle; calibrating a pretrained hyper model using results of the plurality of EIS scans, the hyper model comprising a family of models that each define a voltage response of a respective cell from among a plurality of cells of the battery to a current profile over various states of the batterys; performing a plurality of additional EIS scans on the battery subsequent to the initial use of the battery in the vehicle; recalibrating the pretrained hyper model using results of the plurality of additional EIS scans; predicting the SoP of each of multiple cells of the vehicle battery responsive to a current battery state of each of the multiple cells; combining the SoP of each of multiple cells into a battery SoP; and controlling an amount of current extracted from or put into the battery responsive to a SoP value. 13. The method according to claim 12, wherein the plurality of additional EIS scans are periodically performed subsequent to the initial use of the battery in the vehicle. 14. The method according to claim 12, wherein the various states of the battery include at least some of different temperature ranges of the battery, different states of charge (SoCs) of the battery, age of the battery, and a nature of a current load the battery is subjected to. 15. The method according to claim 12, wherein the hyper model is an equivalent circuit model (ECM) that maps various circuit elements such as resistors, capacitors, inductors, and Warburg impedance of the ECM to the various states of the battery used in the EIS scans. 16. The method according to claim 12, wherein the hyper model is a model reduced from a physics based model that maps physics representations of battery elements to the various states of the battery used in the EIS scans 17. The method according to claim 12, wherein the hyper model is an adaptive filter inferred using a frequency response of an equivalent impedance of the battery learned using the EIS scans under the various states of the battery used in the EIS scans. 18. The method according to claim 12, wherein the results of the EIS scans are used to infer a hyper model using an optimization technique selected from the group consisting of a gradient-based linear optimization method, a non-negative least squares (NNLS) method, and a convex optimization method. 19. The method according to claim 12, wherein the hyper model is a family of equivalent circuit models (ECMs), and the method further comprises performing a smart initialization of a parameter fitting method by setting a series resistor R in an R—RC ECM model having the series resistor R in series with one or more RC parallel sub-circuits to a smallest observed impedance value, setting model parameters to determined values and holding the determined values fixed while scanning over a range based on a value of the series resistor R to identify a determined value that minimizes an objective function. 20. The method according to claim 12, wherein the plurality of additional EIS scans is performed while the battery is in a key-on condition. 21. The method according to claims 12 or 20, further comprising, during a key- on condition of the battery: measuring a current output from each of multiple cells of the battery; and determining the current battery state of each of the multiple cells of the battery responsive to the current output from each of the multiple cells of the battery. 22. The method according to claim 12, further comprising predicting the SoP of at least one of the cells of the battery, wherein at least one constraint comprises at least one of a terminal voltage of the battery, a current of the battery, a temperature of the battery, and a state of charge (SoC) of the battery. 23. The method according to claim 12, wherein the plurality of additional EIS scans is performed while the battery is in a key-off condition to update the hyper model. 24. The method according to claim 12, wherein in a key-off condition, the battery SoP is estimated using a regression algorithm given EIS measurements and a current state of the battery. 25. The method according to claim 12, wherein the SoP comprises a maximum allowable static current that can be sustained for a given time period. 26. The method according to claim 25, wherein the maximum allowable static current comprises a level of current that does not cause a constraint violation. 27. The method according to claim 25, wherein the maximum allowable static current comprises a level of current that does not result in a temporary loss of capacity of more than a specified fraction over a given time period. 28. The method according to claim 12, wherein the SoP comprises a maximum allowable constant power that can be sustained for a given time period. 29. The method according to claim 20, further comprising tracking a battery state starting from a known initial condition responsive to measurable inputs comprising the results of the EIS scans. 30. The method according to claim 20, further comprising tracking a battery state including a terminal voltage of the battery using an Extended Kalman filter (EKF) technique that comprises predicting a battery state vector using a state space model and predicting the terminal voltage of the battery response over a desired period of time. 31. The method according to claim 30, wherein the EKF technique comprises forming battery system matrices expressed as a function of a step index that is dependent on parameters of the hyper model that, in turn, are dependent on the various states of the battery. 32. The method according to claim 12, in a key-on condition, the method further comprises estimating the battery state including a terminal voltage of the battery using a closed form solution to discrete approximations of a state space model and predicting the terminal voltage of the battery response over a desired period of time. 33. The method according to claim 20, further comprising generating a prediction for a current SoP of the battery responsive to a present prediction of the battery state, measurable inputs, and system limitations relating to current and voltage maximum and minimum values. 34. The method according to claim 33, wherein the battery state comprises a present battery current, a present terminal voltage, a present open circuit voltage (OCV), and a present impedance measurement obtained using EIS. 35. The method according to claim 12, wherein, for a recalibration, respective complexities of models in the family of models comprised in the hyper model increase with increasing battery cell age. 36. The method according to claim 12, wherein the SoP of the battery is equal to one of the SoPs of each of multiple cells. 37. The method according to claim 12, wherein the SoP value of a battery is equal to a lowest value from among the multiple cell SoPs for the multiple cells that constitute the battery. 38. The method according to claim 12, wherein the SoP value is equal to the battery SoP. 39: A system for predicting a state of power (SoP) of a battery, the system comprising: an electrochemical impedance spectroscopy (EIS) system for performing a plurality of EIS scans on the battery prior to an initial use of the battery in a vehicle, and a plurality of additional EIS scans on the battery in the vehicle; a memory device for storing program code; and a processing device operatively coupled to the EIS system and the memory device for running the program code to: calibrate a pretrained hyper model using results of the plurality of EIS scans, the pretrained hyper model comprising a family of models that each define a voltage response of a respective cell from among a plurality of cells of the battery to a current profile over various states of the battery; recalibrate the pretrained hyper model using results of the plurality of additional EIS scans; predict the SoP of each of multiple cells of the vehicle battery responsive to a current battery state of each of the multiple cells; combine the SoPs of each of multiple cells into a battery SoP; and control an amount of current extracted from or put into the battery responsive to at least one of the SoPs of each of multiple cells or the battery SoP. 40. The system according to claim 39, wherein the system is comprised in a battery management system. 41. The system according to claim 39, wherein the plurality of additional EIS scans are periodically performed subsequent to the initial use of the battery in the vehicle.
Description:
SYSTEM AND METHOD FOR STATE OF POWER ESTIMATION OF A BATTERY USING IMPEDANCE MEASUREMENTS CROSS REFERENCE TO RELATED APPLICATION [0001] This application is a non-provisional utility patent application based on U.S. Provisional Application Serial No. 63/336,724, entitled “SYSTEM AND METHOD FOR STATE-OF-POWER ESTIMATION OF A BATTERY USING IMPEDANCE MEASUREMENTS”, filed on April 29, 2022, the disclosure of which is incorporated herein by reference in its entirety. FIELD OF THE DISCLOSURE [0002] This disclosure relates generally to battery monitoring and, more particularly, to a system and method for estimating battery state of power (SoP) using battery impedance measurements. [0003] As the development of electric vehicles grows at rapid pace, battery management systems (BMS) are expected to track the state of battery packs to enable effective operation of the corresponding control algorithms. One such state to monitor is the maximum amount of power that can be drawn from (during rapid acceleration) or put back (during regenerative breaking) into the pack, which translates directly to vehicle performance and user experience. Currently, BMS’ use simplistic and conservative estimates computed as the product of the minimal terminal voltage with the maximum current that can be drawn from it for a given period of time. However, this significantly underestimates the maximal power metric. [0004] State of Power (SoP) quantifies the maximum amount of power that a battery can deliver over a short period of time. Traditional battery management systems infer SoP-like information by training a model based on time-domain voltage and current data, often in the form of a current pulse schedule (i.e., Direct Current Internal Resistance (DCIR) measurements), where SoP estimation is accomplished through lengthy DCIR measurement tests involving current pulses with long relaxation times. This can be very time-consuming and is impractical for several reasons, one being that the resulting equivalent circuit model (ECM) must be recalibrated after the cell has aged somewhat. Thus, there is a need for improved SoP measurement and improved ECM recalibration. SUMMARY [0005] The following presents a simplified summary of one or more aspects in order to provide a basic understanding of such aspects. This summary is not an extensive overview of all contemplated aspects and is intended to neither identify key or critical elements of all aspects nor delineate the scope of any or all aspects. Its sole purpose is to present some concepts of one or more aspects in a simplified form as a prelude to the more detailed description that is presented later. [0006] In this exposition we define a hyper model as a family of models that each characterize the voltage response of a respective cell from among a plurality of cells of the vehicle battery to a current profile under varying battery states. Thus, for any given battery state, we have an equivalent state space model that represents the electrical characteristics of the battery accurately. According to an aspect, a method is provided for pretraining a hyper model that may be used to predict the state of power (SoP) of a vehicle battery. The method includes performing electrochemical impedance spectroscopy (EIS) scans on a plurality of batteries similar to the vehicle battery under various states of the battery. In one aspect, the various battery states can include at least some of different temperature ranges, different states of charge (SoCs), age of the battery, and the nature of the current load the battery is subjected to. The method further includes fitting parameters of the hyper model by applying an optimization technique to results of the EIS scans. [0007] According to another aspect, a method is provided for predicting a state of power (SoP) of a battery. The method includes performing a plurality of electrochemical impedance spectroscopy (EIS) scans on the battery prior to an initial use of the battery in a vehicle. The method further includes calibrating a pretrained hyper model using results of the EIS scans. The pretrained hyper model includes a family of models that each define a voltage response of a respective cell from among a plurality of cells of the battery to a current profile over various states of the battery. The method further includes performing a plurality of additional EIS scans on the battery subsequent to the initial use of the battery in the vehicle. The method also includes recalibrating the pretrained hyper model using results of the additional EIS scans. [0008] According to a further aspect, a system is provided for predicting a state of power (SoP) of a battery. The system includes an electrochemical impedance spectroscopy (EIS) system for performing a plurality of EIS scans on the battery prior to an initial use of the battery in a vehicle, and a plurality of additional EIS scans on the battery subsequent to the initial use of the battery in the vehicle. The system further includes a memory device for storing program code. The system also includes a processing device operatively coupled to the EIS system and the memory device for running the program code to calibrate a pretrained hyper model using results of the plurality of EIS scans. The pretrained hyper model includes a family of models that each define a voltage response of a respective cell from among a plurality of cells of the battery to a current profile over various states of the battery. The processor device additionally runs the program code to recalibrate the pretrained hyper model using results of the plurality of additional EIS scans. [0009] To the accomplishment of the foregoing and related ends, the one or more aspects comprise the features hereinafter fully described and particularly pointed out in the claims. The following description and the annexed drawings set forth in detail certain illustrative features of the one or more aspects. These features are indicative, however, of but a few of the various ways in which the principles of various aspects may be employed, and this description is intended to include all such aspects and their equivalents. BRIEF DESCRIPTION OF THE DRAWINGS [0010] The disclosed aspects will hereinafter be described in conjunction with the appended drawings, provided to illustrate and not to limit the disclosed aspects, wherein like designations denote like elements, wherein dashed lines may indicate optional elements, and in which: [0011] FIG. 1 is a diagram of a battery system within a vehicle, in accordance with an exemplary aspect; [0012] FIG.2 is a schematic diagram showing a battery equivalent circuit model (ECM), in accordance with an exemplary aspect; [0013] FIG. 3 is a diagram showing an example Nyquist plot of the negative of the imaginary part of the EIS versus real part of the impedance, showing the measured impedance and that estimated from the ECM of FIG. 2 for a battery, in accordance with an exemplary aspect; [0014] FIG.4 is a diagram showing a plot of terminal voltage versus time for the battery ECM of FIG.2, in accordance with an exemplary aspect; [0015] FIG.5 is a schematic diagram showing another battery ECM, in accordance with an exemplary aspect; [0016] FIGS.6-8 and 10-11 are plots showing the estimated parameters of the ECM in FIG.5 computed using non-linear least squares (NLLS) as a function of states of charge (SoC) and temperature (T), in accordance with an exemplary aspect; [0017] FIG.9 is a diagram showing the Nyquist plots of the negative of the imaginary part of the cell impedance versus the real part of the cell impedance; [0018] FIG. 12 is a plot of open circuit voltage (OCV) versus SoC, in accordance with an exemplary aspect; [0019] FIG. 13 is a diagram showing an extended Kalman Filter (EKF) functioning, in accordance with an exemplary aspect; [0020] FIG. 14 is a block diagram showing a key-off use case, in accordance with an exemplary aspect; [0021] FIG. 15 is a block diagram showing a key-on use case, in accordance with an exemplary aspect; [0022] FIG.16 is a diagram showing a plot of current versus time for the battery ECM of FIG.2; [0023] FIG.17 is a diagram showing a plot of voltage versus time for the battery ECM of FIG.2; [0024] FIG.18 is a flow diagram illustrating a method for offline pretraining of hyper model parameters, in accordance with an exemplary aspect; [0025] FIG. 19 is a flow diagram showing a method for calibration and recalibration of a new cell, in accordance with an exemplary aspect; and

[0026] FIGS. 20-21 are flow diagrams further showing a step the method of FIG.19, in accordance with an exemplary aspect. DETAILED DESCRIPTION OF EXAMPLE EMBODIMENTS [0027] In various aspects, systems and methods are provided for estimating the state of power (SoP) of a battery using electrochemical impedance spectroscopy (EIS) measurement technology. In various aspects, the battery state of power (SoP) is estimated using impedance measurements taken at multiple frequencies. Battery SoP is defined as the maximal power that battery can generate or absorb at any point in time without exceeding manufacturers limits, such as those on maximal current, maximal and minimal terminal voltage and temperature, etc. In various aspects, the impedance measurements may be used to train a family of equivalent circuit models (ECM), labeled an “ECM hyper model”. The hyper model defines the voltage response of each cell to a respective current profile throughout varying conditions (for example, varying combinations of temperatures and states of charge, or in another example, varying combinations of ages of the cell). When the cell is in operation, the model is used to track the internal state of the cell, giving access to a SoP estimate of that cell at any point in time through a prediction mechanism. Alternatively, a static model can be defined for relevant inputs that are measured in situ. [0028] Thus, the hyper model may be an equivalent circuit model (ECM) that may map various resistive and capacitive components of the ECM to the varying combinations of temperatures and states of charge used in the EIS scans. While one or more aspects may be directed to the use of a battery equivalent circuit model (ECM) to model a battery, other representations can be used such as, for example, and not limited to, a model reduced from a physics based model that maps physics representations of battery elements to various combinations of temperatures and states of charge used in the EIS scans, and so forth. It is to be appreciated that aspects of the present disclosure may use any type of battery model to obtain a SoP estimate in accordance with the present disclosure. [0029] Aspects described herein are designed to report how much power a battery can generate or absorb at any point in time using information from electrochemical impedance spectroscopy (EIS) because exceeding certain limits on the power applied to/from an EV battery during charge/discharge can lead to premature aging, electrochemical faults, or serious problems like thermal runaway. Additionally, traditional BMS systems infer SoP-like information by training a model based on time- domain voltage and current data, often in the form of a current pulse schedule (i.e., Direct Current Internal Resistance (DCIR) measurement). In another example, traditional BMS systems infer SoP-like information by training a model based on time- domain voltage and current data from hybrid pulse power characterization (HPPC) current schedules. Moreover, EIS scan measurements enable more complete and precise modeling of a battery’s behavior by training models on impedance spectra. [0030] EIS is an emerging technology in Battery Management Systems (BMS’). In EIS, a cell is stimulated with a sinusoidal current and the resulting voltage is measured (or vice versa). This is repeated at many different frequencies to produce an impedance spectrum that can be used to fit models characterizing the electrochemical properties of the cell. [0031] Aspects described herein include collecting impedance spectra while a car is at rest (“key-off”) at various state-of-charge (SoC) levels and temperatures and using this data to train a family of equivalent circuit models (ECM), referred to herein as an “ECM hyper model”. Impedance spectra measurements can also be collected on a corpus of batteries before they are placed in the car, while the in situ measurements (“key-off”) are then used to calibrate the model. The hyper model defines the voltage response of the cell to a current profile throughout varying conditions. When the cell is in use, the model is used to track the internal state of the cell, enabling SoP estimation at any point in time through a prediction mechanism. Alternatively, a static model can be defined for relevant inputs that are measured in situ. [0032] SoP estimation can be accomplished through lengthy DCIR tests involving current pulses with long relaxation times. This can be very time-consuming and is impractical for several reasons, one being that the resulting ECM must be calibrated after the cell has aged somewhat. The use of EIS spectra is much more practical because it is relatively fast and does not involve stimulating the cell with a large current pulse. Recalibration of the ECM throughout the life of the cell is much easier because EIS measurements may be taken while the car is at rest (a.k.a., “key- off”). The concept of EIS-based ECMs for SoP estimation is represented in FIGs.1-4. [0033] Referring to FIG.1, a battery system 102 within a vehicle 100 is shown. The vehicle 100 may be an electric vehicle (EV) or a hybrid vehicle such as a plug-in hybrid electric vehicle. The vehicle 100 may be considered a host of the battery system 102. The battery system 102 may include a BMS controller 110, a wireless battery management system (wBMS) 120, and a plurality of battery modules 130. [0034] The BMS controller 110 may be a component of the vehicle 100 configured to interface with the wBMS 120. For example, the BMS controller 110 may be an electronic control unit (ECU) of the vehicle 100. The BMS controller 110 may execute a BMS controller application, which may be referred to as a safety application. For instance, the BMS controller 110 may communicate with the wBMS 120 to receive information about the battery system such as state of charge, voltage, temperatures, and any faults that have occurred. The BMS controller 110 may also communicate with other vehicle components such as an inverter or charger. [0035] The wBMS 120 may be configured to interface between the battery modules 130 and the BMS controller 110. For example, the wBMS 120 may receive packets including measurement messages and fault messages from the battery modules 130. The wBMS 120 may aggregate the messages from the individual battery modules to provide system level information to the BMS controller 110. The wBMS 120 may include a wireless manager 122 and a wireless radio 124. The wireless manager 122 may be configured to generate messages for transmission to the battery modules 130 and receive messages from the battery modules 130. The wireless manager 122 may include a radio protocol stack. The wireless radio 124 may include one or more radios configured to transmit radio-frequency (RF) signals to the battery modules 130. The wireless radio 124 may be referred to as a head radio and may control (e.g., schedule) the communications with the battery modules 130. [0036] The battery module 130 may include a wireless radio 140, a safety processor 150, a battery monitoring system 160, and a plurality of cells 170. The cells 170 may be battery cells that store power. In some implementations, each battery module 130 may include between 3 and 24 individual cells 170. [0037] The battery monitoring system 160 may be configured to monitor one or more parameters of the plurality of battery cells. For example, the battery monitoring system 160 may monitor voltage and temperature of each cell. The battery monitoring system 160 may provide the measurements to the safety processor 150. The battery monitoring system 160 may be referred to as a battery monitoring integrated circuit (BMIC). [0038] In various aspects, battery monitoring system 160 includes an electrochemical impedance spectroscopy (EIS) system 125 for taking EIS measurements at various frequencies. The measurements may include, among other variables, voltage and current from which an impedance spectrum of each cell can be calculated. Such calculation can be performed by the safety processor 150 in order to calculate SoP using the measurements to control how the battery is used. In this way, premature battery damage among other issues can be avoided. [0039] The safety processor 150 may be a computer processor configured to execute computer code such as a script. In some implementations, the safety processor 150 is a separate processor connected to the wireless radio 140 and the battery monitoring system 160 via interfaces such as a serial peripheral interface (SPI). In other implementations, the safety processor may be a processor of the wireless radio 140. The safety processor 150 may be configured to perform various tasks with respect to managing the battery module 130. For example, the safety processor 150 may be configured to execute a battery monitoring script to generate a battery information payload defined by a battery information payload format. The battery information payload may be based on the SoP calculated using the EIS measurements from the EIS system 125. The safety processor may also receive commands (e.g., a cell balancing command) from the wBMS 120, and write commands to the battery monitoring system 160. The safety processor 150 may include a script engine 152, a data processing engine 154, a scheduler 156, and a parser. The script engine 152 may execute a script 191 received from wBMS 120. The data processing engine 154 may perform data processing commands defined by the script to write commands to the battery monitoring system 160 and process data received from the battery monitoring system 160. In some implementations, the data processing engine may be a component of the script engine. The scheduler 156 may be configured to generate a plurality of packets according to a battery information payload format. For example, the scheduler 156 may operate according to an initialization schedule at initialization to potentially generate fault packs and operate according to an active mode schedule to generate measurement messages within a timed subloop. The parser 158 may be configured to parse various components of the script. [0040] The wireless radio 140 may be configured to communicate with the wBMS 120 and/or a service device. The wireless radio 140 may be a radio similar to the wireless radio 124 configured to transmit and receive RF signals. The wireless radio 140 may be referred to a tail radio because the wireless radio 140 may be managed by the wireless radio 124. The wireless radio 140 is configured to receive a container file 190 from the wBMS 120. The container file 190 includes metadata 192 defining the battery information payload format and the battery monitoring script. The wireless radio 140 is configured to transmit at least the metadata 192 of the container file 190 to the service device. The wireless radio 140 is also configured to transmit packets 181 including a battery information payload to either the wBMS 120 or the service device. [0041] To calculate SoP, a model of the battery is used. In some aspects, the model is a battery equivalent circuit model (ECM). In other aspects, the model is a model reduced from a physics-based model. Other model representations of a battery can also be used, while maintaining the spirit of the present disclosure. [0042] Referring to FIG.2, a battery equivalent circuit model (ECM) 200 is shown, in accordance with an exemplary aspect. The battery ECM 200 includes an R section 210 and an RC section 220. While one RC section 220 is shown, other circuit models can have more than one RC section 220, as shown in FIG.5. [0043] The R section 210 of ECM 200 includes resistor R 0 211. The RC section 220 includes resistor R1221 and capacitor C1222 in parallel and having voltage VC1 across them. A positive terminal of a voltage VOCV 230 is connected to one side the resistor R 0 211. A positive terminal of a battery terminal voltage V term is connected to one side of the RC section 220. The other side of RC section 220 is connected to another side of the resister R0211. A negative terminal of the battery terminal voltage Vterm is connected to a negative terminal of the voltage V OCV 230. [0044] Referring to FIG. 3, a Nyquist plot 300 of the negative of the imaginary part of the electrochemical impedance spectrum of a battery versus the real part of the impedance for the battery ECM 200 of FIG.2 is shown, in accordance with an exemplary aspect. Curve 310 denotes impedance as computed using the estimated ECM model of the battery and curve 320 denotes the measured impedance of the battery. [0045] Referring to FIG.4, a plot 400 of terminal voltage versus time for the battery ECM 200 of FIG.2 is shown, in accordance with an exemplary aspect. In particular, time is shown in the x-axis and terminal voltage (V term ) is shown in the y axis. Two curves are shown, a predicted curve 420 and a true curve 410. [0046] Referring to FIG. 5, another battery equivalent circuit model (ECM) 500 is shown, in accordance with an exemplary aspect. The battery ECM 500 can be used to accommodate edge devices, e.g., in a wireless BMS such as that shown and described with respect to FIG.1. The battery ECM 500 includes an R section 510, a first RC section 520, and a second RC section 530. The R section 510 includes resistor R0 511. The first RC section 520 includes resistor R1 521 and capacitor C1 522 in parallel and having voltage V C1 across them. The second RC section 530 includes resistor R2531 and capacitor C2532 in parallel and having voltage VC2 across them. A positive terminal of a voltage V OCV 540 is connected to one side of the resistor R o 511. A positive terminal of a battery terminal voltage V term Is connected to one side of the second RC section 530. The other side of the second RC section 530 is connected to another side of first RC section 520. The other side of first RC section 520 is connected to another side of the resister R 0 511. A negative terminal of the battery terminal voltage Vterm is connected to a negative terminal of the voltage VOCV 540. [0047] For the ECM 500 of FIG. 5, battery impedance is given by the following: [0048] Using EIS, this impedance is measured at various radial frequencies (ωk) and the ECM can be inferred from the constrained non-linear least squares estimate of the parameters. [0049] Referring to FIGS. 6-8 and FIGS. 10-11, estimates of ECM parameters as a function of temperature and state of charge (SoC) for an R-RC-RC approximation of a 3Ah Li-ion battery are shown. The estimates were derived using EIS sweeps made at different temperatures and SoCs. [0050] The plot in FIG. 9, shows the estimated battery impedance (the Nyquist plot) in comparison to the measured EIS at various temperatures at 50% SoC. This clearly shows that the estimated ECM approximates the battery impedance well. In particular, FIG. 9 shows the Nyquist plots of the negative of the imaginary part of the cell impedance versus the real part of the cell impedance which, in turn, shows the comparison of the measured impedance through EIS and that estimated from the ECM of FIG.4 where estimated parameters of the ECM in FIG.5 were computed using non- linear least squares (NLLS) as a function of temperature (T), in accordance with an exemplary aspect; [0051] EIS analysis uses a computer to find the model parameters that give the best agreement between a model's impedance spectrum and a measured spectrum. A non-linear least squares fitting (NLLS) algorithm may be used. NLLS starts with initial estimates for all the model’s parameters. Starting from this initial point, the algorithm makes changes in several or all of the parameter values and evaluates the resulting fit. If the change improves the fit, the new parameter value is accepted. If the change worsens the fit, the old parameter value is retained. Next, a different parameter value is changed and the test is repeated. Each trial with new values is called an iteration. Iterations continue until the goodness of fit exceeds an acceptance criterion, or until the number of iterations reaches a limit. [0052] FIGS. 6-11 show an ECM fit using NLLS as a function of SoC and temperature, and also show that the circuit parameters are a function of the battery state (T, SoC, Age). For instance, the cell resistances decrease with increasing temperature, whereas the opposite is broadly true of the capacitances. To account for this variation, it is possible to infer an equivalent hyper model of a battery which maps the battery state to the ECM such as the following: [0053] Such a hyper model also helps us interpolate effectively across different battery states without requiring explicit experimentation. [0054] Alternatively, projected gradient descent may be used to track the ECM parameters by using EIS measurements. The gradients for the algorithm may be derived based on the partial derivatives of the estimate ^^ as follows: [0055] A time-domain voltage prediction technique will now be described, in accordance with an exemplary aspect, according to an exemplary aspect. To track the cell terminal voltage, the mathematical model is first derived as a function of the ECM parameters and the current applied to the battery as follows: [0056] Here, the component voltage dynamics are described by the following: . [0057] Thus, tracking equivalently translates to tracking the OCV and The OCV is typically a non-linear function of the SoC, such as shown in FIG. 12 which shows the OCV-SoC for a battery, and the evolves as follows: where ^^ ^^ ^^ is the battery capacity. The ECM parameters also vary as non-linear functions of the battery state. Accounting for these, an Extended Kalman Filter (EKF) is described to track [0058] A description will now be given regarding EKF-based terminal voltage tracking, in accordance with an exemplary aspect. [0059] The EKF technique includes forming battery system matrices expressed as a function of a step index that is dependent on parameters of the hyper model that, in turn, are dependent on the varying conditions. [0060] Discretizing the differential equations, the state space model for the ECM 500 in FIG.5 is given by the following: where The measurement is given by the following: [ ] [ ] [ ]([ ] [ ]) [ ] [ ] [0061] Thus, state estimates are corrected when measurements are available, and when measurements are unavailable, the state vector can be predicted using the state space model. Consider the state and measurement models as follows: ( ) where the noise processes are normal random vectors given by Then an EKF, functioning 1300 as shown in Fig. 13, can track the terminal voltage of the battery. [0062] The EKF functioning 1300 includes a predict portion 1310, a linearize portion 1320, and a correct portion 1330. [0063] It may be noted here that alternate approaches other than an EKF can be used to track and predict the battery state including adaptive filtering and autoregressive models that predict the terminal voltage of the battery given the current state of the system. The methodology described here to compute the state of power using the hyper model of a battery is compatible with any such method. [0064] A description will now be given regarding estimating the SoP of a battery, in accordance with an exemplary aspect. [0065] The SoP of a battery measures the maximum amount of power that can be delivered by the battery. Here, two possible definitions are described and methods to calculate the SoP. One definition is referred to as constant current SoP, and the other definition is referred to as constant power SoP. [0066] Regarding constant current SoP, a conservative estimate of the maximal power that can be drawn from the battery may be defined as follows: where is the maximum current pulse that can be drawn from the battery. The idea is that the maximal cell current constrained by the minimum terminal voltage, results in the maximal power that can be drawn from the battery. [0067] be the terminal voltage estimate as a function of current ^^, over time, where is the time period over which SoP is to be computed. Then, the maximum value of i is the solution to the following: [0068] In an aspect, the Newton-Raphson method can be used to find the root of this equation, and hence compute the SoP. This however is a conservative estimate as it does not account for instantaneous power fluctuations on account of the variation in [0069] Regarding constant power SoP, the maximal constant power that can be drawn from a battery amounts to more net energy output than the constant current estimate as it adapts to the variation in Thus, a more accurate estimate of the maximal energy that can be drawn from a battery is derived by solving for the instantaneous currents that results in a constant power over the time period of interest. This power is computed by using Powell’s conjugate direction method while constraining to the safety limits of the battery [0070] Unlike aspects described herein, conventional techniques for fitting ECMs and corresponding methods are current pulse-based measurement techniques. [0071] In accordance with features of aspects described herein, EIS measurements may be used for training of parameters of a general-purpose algorithmic model to be used for in-situ (while an electric vehicle (EV) is in-use) SoP prediction. Additionally, EIS measurements are used at key-off (for example, when an EV is parked in a garage) to recalibrate parameters of the SoP algorithmic model. Moreover, EIS measurements are used during EV in-use scenario to dynamically recalibrate parameters of the SoP model. [0072] A primary purpose of aspects described herein is to estimate how much power a cell can output at any given time. SoP can be quantified indirectly as the maximum allowable static current that can be sustained for Δt seconds (or some prespecified time period). "Maximum allowable" is defined as the largest current that does not cause a constraint violation. Constraints can be defined for terminal voltage, current, state-of-charge, cell temperature, etc. Alternatively, SoP can be quantified indirectly as the maximum allowable static current that can be sustained for Δt seconds (or some prespecified time period). "Maximum allowable" is defined as desired, e.g., an amount that does not result in a temporary loss of capacity of more than a specified fraction (e.g., 5%). SoP (or max allowable current) is a dynamic variable and depends on the initial state of the battery at the moment of prediction. In other words, any system that maps measurable parameters to SoP would need to include memory in order to achieve acceptable accuracy, unless battery state is supplied as an additional input. [0073] SoP estimation use cases include key-off and key-on. For key- off, a model is developed that generates a prediction of SoP for an a priori given “battery state.” Battery state could be defined, for example, as the most recent battery current, terminal voltage, open circuit voltage (OCV), state of charge (SoC), temperature, and impedance measurement (this might be different for individual batteries and the whole pack). The SoP model is a multidimensional look-up table that maps each expected battery state into a corresponding SoP value (could also be implemented as a function that performs this mapping/calculation when called – that is, not necessarily a hard encoded table that is just read from memory when called). [0074] For the key-on use case, the model may be equipped with a tracker which starts from a known initial condition and then tracks “battery state.” Whenever SoP prediction is needed, the model can use the most recent battery state supplied by the tracker. The SoP output of the model is instantaneous (i.e., evaluated in real-time) and provides a prediction for the current state of the battery only. [0075] Referring to FIG. 14 a key-off use case 1400 is shown, in accordance with an exemplary aspect. The key-off use case 1400 involves a State of Power (SoP) model 1410 which inputs measurable inputs (EIS), expected battery state, and system limitations (e.g., maximum/minimum current/voltage (I/V) values), and outputs a predicted SoP. [0076] In the key-off use-case, a SoP mapping is generated offline and so there is no need for measurables at the time of generating the mapping (except maybe EIS to correct for ageing factors, as shown above); however, even in the key-off use case, it is implicitly assumed that the SoP information will be eventually used for key-on purposes. [0077] Referring to FIG. 15, a key-on use case 1500 is shown, in accordance with an exemplary aspect. The key-on use case 1500 involves a SoP model 1510 and a state tracker 1520. The SoP model 1510 and the state tracker 1520 both input measurable inputs (EIS, V, I, time (T), SoC, state of health (SoH), and so forth). The SoP model 1510 further inputs predicted battery state output from the state tracker 1520 and system limitations (e.g., maximum/minimum current/voltage (I/V) values), and outputs a predicted SoP. [0078] In the key-on use case, battery state is tracked with a state tracker. In this case, more measurables, as well as system limitations, are provided to the model. Battery state, which affects SoP at a given time, is tracked over time. It will be recognized that SoP is relevant to both discharging and fast charging the battery. [0079] It will be recognized that the proposed model provides battery module- level info, while SoP is a battery pack-level metric; therefore, the individual module- level SoP measurements must be combined appropriately to produce a pack-level SoP. [0080] The following are SoP estimation approaches for key-off and key-on use cases. For the key-off use case, a regression model approach includes training a baseline regression model when at rest (offline) and generating SoP prediction by evaluating the regression function given predicted observable quantities and battery state. An ECM approach includes fitting an ECM to EIS measurements when at rest and generating SoP prediction by evaluating the ECM given predicted observable quantities and battery state. [0081] For the key-on use case, a regression model approach includes training a baseline regression model when at rest (offline) and evaluating real-time SoP by evaluating the regression function given observable quantities. An ECM with Kalman Filter (KF) tracking approach includes fitting an ECM to EIS measurements when at rest (offline), using KF to track battery state, and evaluating a prediction based on the ECM (with constraints) to infer real-time SoP. [0082] Referring back to FIG. 2 and further to FIGS.16-17, the use of ECMs to perform SoP estimation is shown, in accordance with an exemplary aspect. [0083] FIG. 2 as described above relates to battery ECM 200 used by the present disclosure in one exemplary aspect. Other battery ECMs can also be used while maintaining the spirit of the present disclosure such as that shown in FIG.5 and others as described herein and as readily envisioned by one of ordinary skill in the art given the teachings of the present disclosure provided herein. [0084] FIG. 16 shows a plot 1600 of current versus time for the battery ECM 200 of FIG. 2. In particular, time is represented in the x-axis, and current (including Icharge and Imax) is represented in the y-axis. [0085] FIG.17 shows a plot 1700 of voltage versus time for time for the battery ECM 200 of FIG. 2. In particular, time is represented in the x-axis, and voltage (including V charge and V max ) is represented in the y-axis. [0086] First, the battery is modeled with R section 210 and RC section 220 and a voltage source V OCV 230 if the ECM 200 of FIG. 2 is used. Alternatively, another model or representation can be used such as, but not limited to the ECM 500 of FIG.5. The model is fit to the data taking into account parameter variations with temperature, SoC, and age. The parameterization of the model can be in the form of a look-up-table or regression function. Voltage response to a current step of desired duration (e.g., 10 secs, 60 secs, etc.) is calculated from an impedance spectrum that is determined at many frequencies and that can be used to fit models characterizing the electrochemical properties of the cell. The parameterization can be performed in closed-form or via sequential inference with a KF algorithm. For increased accuracy, ECM parameters can be adapted (closed-form approximations or on-line system identification in a KF framework). Constraints are set on current, SoC, voltage, etc., to find max current; this gives P = I * V. [0087] Referring to FIGs.18-19, a method 1800 for offline pretraining of hyper model parameters is shown, in accordance with an exemplary aspect. [0088] At step 1810, perform EIS scans on a corpus of batteries at selected temperatures and SoCs. This may result in the generation of an impedance spectrum including current and voltage measurements taken at different frequencies for each of the cells of a vehicle battery. The batteries in the corpus have a set of similar operating characteristics to the vehicle battery. The corpus of batteries having a set of similar operating characteristics to the vehicle battery refers to batteries having comparable battery chemistry, form factor, battery capacity, and operating conditions. [0089] In an aspect, step 1810 may include one or more of step 1810A through 1810D. [0090] Step 1810A includes, when the hyper model is an equivalent circuit model (ECM) hyper model, performing a smart initialization of the parameter fitting method by setting a series resistor R in an R—RC ECM model have the series resistor R in series with one or more RC parallel sub-circuits to a smallest observed impedance value, setting model parameters to determined values and holding the determined values fixed while scanning over a range based on a value of the series resistor R to identify a determined value that minimizes an objective function. [0091] At step 1810B use, as the hyper model, an ECM that maps various circuit elements such as resistors, capacitors, inductors, and Warburg impedance of the ECM to the various states of the vehicle battery used in the EIS scans. [0092] At step 1810C use, as the hyper model, a model reduced from a physics based model that maps physics representations of battery elements to the various states of the vehicle battery used in the EIS scans. [0093] At step 1810D, use, as the hyper model, an adaptive filter inferred using a frequency response of an equivalent impedance of the vehicle battery learned using the EIS scans under the various states of the vehicle battery. In an aspect, the various states of the battery can include at least some of different temperature ranges of the battery, different states of charge (SoCs) of the battery, age of the battery, and a nature of a current load the battery is subjected to. [0094] At step 1820, fit parameters of a hyper model by applying an optimization technique to results of the EIS scans. The hyper model includes a family of models, where each of the models define a voltage response of a respective cell of the vehicle battery to a current profile over varying conditions. The family of models may include equivalent circuit models (ECMs), models from a physics model or some other type of battery model. The fitting may be performing using a multidimensional look up table and/or a regression function. The results of the EIS scans may be applied to any one of more of the following optimization techniques: a gradient-based linear optimization method; a non-negative least squares (NNLS) method; and a convex optimization method. [0095] At 1830, estimate the state of power (SoP) of the vehicle battery using the hyper model, the results of the EIS scans, operational constraints of the vehicle battery, and the various states of the vehicle battery used in the EIS scans. [0096] Referring to FIG.20, a method 2000 for calibration and recalibration of a new cell is shown, in accordance with an exemplary aspect. [0097] At step 2010, perform a plurality of electrochemical impedance spectroscopy (EIS) scans on a battery prior to an initial use of the battery in a vehicle (e.g., at assembly). This may result in the generation of an impedance spectrum including current and voltage measurements taken at different frequencies for each of the cells of the battery. [0098] At step 2020, calibrate a pretrained hyper model of the battery using results of the EIS scans. The pretrained hyper model includes a family of models, where each of the models define a voltage response of a respective cell of the vehicle battery to a current profile over varying conditions. The family of models may include equivalent circuit models (ECMs), models from a physics model or some other type of battery model. [0099] At step 2030, cycle the battery (normal use). [00100] At step 2040, collect additional EIS scans in in-situ at a key-off condition (that is, subsequent to the initial use of the battery in the vehicle). [00101] At step 2050, recalibrate the pretrained hyper model using results of the additional EIS scans. In an aspect, for a recalibration, respective complexities of models in the family of models included in the hyper model may increase with increasing battery age. In this way, the individual differences between each of the cells forming a battery can be accounted for and their respective aging and corresponding affects can be considered. [00102] It will be noted that EIS scans can also be collected during in- use (i.e., key-on) rather than key-off in accordance with certain aspects. [00103] Referring to FIGS.21-22, step 2030 of method 2000 of FIG.20 is further shown, in accordance with an exemplary aspect. [00104] At step 2131, measure a current which is output from each cell of the battery. [00105] At step 2132, step forward in the battery state tracker (if key-on use case) to predict a current state of each cell of the battery. [00106] At step 2133, predict the SoP of each cell of the battery responsive to the current state of each cell of the battery. [00107] In an aspect, step 2133 can include step 2133A. [00108] At step 2133A, consider at least one constraint in the SoP prediction of each cell of the battery. In an aspect, at least one constraint may include at least one of a terminal volage of the battery, a current of the battery, a temperature of the battery, and a state of charge (SoC) of the battery. [00109] At step 2134, combine the SoPs of each of the cells of the battery to determine a battery SoP. It is to be appreciated that various methods can be used to combine the individual cell SoPs into a battery SoP. A simple sum or more advanced combining techniques can be used including weighted sums and so forth while maintaining the spirit of the present disclosure. [00110] At step 2135, perform an action with respect to the battery responsive to the individual cell SoPs and/or the battery SoP. [00111] In an aspect, step 2135 can include one or more of steps 2035A and 2035B. [00112] At step 2135A, control an amount of current extracted from or put into the battery responsive to the battery SoP and/or a SoP having a lowest value from among multiple cell SoPs for the multiple cells that constitute the battery. [00113] At step 2135B, perform a service level action on the battery including replacing the battery responsive to the SoP being below a threshold value and/or replacing individual cells having an individual SoP below the same or another threshold value to optimize the performance of a given battery by replacing it’s weak link(s) (cell(s)). [00114] A further description of model fitting will now be given, in accordance with an exemplary aspect. [00115] To fit this model to impedance measurements ^^ ^^ at multiple frequencies, an objective function is defined as follows: where is the complex-valued error of the fit and re and im subscripts denote real and imaginary parts. Its derivatives with respect to the (real-valued) model parameters are as follows: [00116] To simplify the derivations, the following shorthand notation is used: [00117] The individual partial derivatives of the EIS model with respect to each parameter is as follows: [00118] The impedance model formula can also be abbreviated as follows: [00119] Given these formulas for the error and its derivatives, a non- linear optimization technique can be applied to tune the parameters to a dataset. It has to be ensured that the parameters remain non-negative. [00120] A description will now be given regarding log-scaled parameters, in accordance with an exemplary aspect. [00121] It is more natural to adjust the model parameters on a logarithmic scale, which is as simple as defining for logarithmic-scale parameter counterparts This also has the added benefit of ensuring non-negativity. Conveniently, so applying the chain rules gives:

[00122] A description will now be given regarding R-T parameterization, in accordance with an exemplary aspect.

[00123] The model can be re-parameterized via r ( = the impedance model is as follows: and its derivatives are as follows:

[00124] A description will now be given regarding additional considerations, in accordance with an exemplary aspect.

[00125] A line search at each iteration ensures that the error is never increased. Also, well-known measures to improve convergence behavior can be implemented like adaptive step size and momentum.

[00126] An initialization can be used. The initialization is particularly applicable to R-RC-type models in order to get in the vicinity of an optimum. R o is set to the smallest observed impedance value. Then the τ parameters are set to a reasonable intermediate value (empirically / visually determined). Holding those, fixed, a scan is performed over a range of R i i > 0 values with all of them held equal to each other to find the one that minimizes the objective. Finally, the Ri and r ( (i > 0) values are deterministically spread over a small range around the initial values to avoid ambiguities. This seems to work quite well in practice. This results in a smart initialization of the R-RC type models.

[00127] A further description will now be given of the hyper model, in accordance with an exemplary aspect.

[00128] In practice, a model of impedance is wanted that generalizes to many temperatures and SoCs because it is known that the model parameters vary quite a lot with respect to these states. If an RC model is fit with a standardized, deterministic initialization procedure, a permutation problem can be avoided between the RC pairs that could severely complicate the generalization procedure. Once all the individual models are fit, a 3 rd -order polynomial can be fit via least-squares regression to each parameter. [00129] A description will now be given regarding time-domain voltage prediction, in accordance with an exemplary aspect. [00130] Given an ECM, the time-domain voltage response to a current input can be expressed. For an R-RC-type model, the following equations from voltage and current laws can be used: [00131] The latter can be re-written as follows: [00132] Thus, only the voltages over the RC pairs need to be solved to be able to predict and it is also desired that those voltage functions be discretized. [00133] A description will now be given regarding Euler approximations, in accordance with an exemplary aspect. [00134] Applying the forward Euler approximation, the following is obtained: and applying backward Euler, the following is obtained: [00135] Both of these solutions describe a convex combination of the voltage and the voltage but differ in the mixing coefficients and whether to use the present or previous step’s current. Difference equations of this form, i.e.: (34) have the following closed-form solution in the time-domain for constant u[n]: (35)

[00136] This equation may be used to predict what the voltage will be at any time in the future given a constant current and initial voltage condition.

[00137] A description will now be given regarding closed-form voltage prediction for the discretized model, in accordance with an exemplary aspect.

[00138] Substituting the solution for the RC pair voltage from the backward Euler approximation into the terminal voltage formula, the following is obtained: for a constant current i* using and the shorthand:

[00139] There are two issues with this expression:

[00140] SoC dependence: It does not take into account the dependence of the

ECM parameters on SoC, which could change non-negligibly during the prediction horizon depending on the current magnitude. This is addressed by an SoC-dependent ECM, and as Cj is a constant, this equation does not work for varying SoCs.

[00141] Initial conditions: Initial values have to be provided for the voltages across each of the RC pairs. This is handled by an appropriate tracking scheme.

[00142] A description will now be given regarding explicit solution to the continuous-time model, in accordance with an exemplary aspect.

[00143] Since the full system model is completely decoupled, and each ODE is linear and of first order, this system can be solved explicitly to get expressions for voltages Namely, from (27), the following is directly obtained:

[00144] Assuming that i(t) = i 0 is constant in the interval of interest, the following is obtained:

[00145] Similarly, if the system parameters R j and τ j are time-varying (e.g., if they depend on SOC or cell temperature), the explicit solution for V c . can be written as follows:

Where is the antiderivative of If model parameters can be explicitly described as functions of time then the above integrals can be either explicitly solved, or numerically integrated. For example, given starting SOC and cell temperature, it should be possible to approximate model parameters with low order polynomials (e.g., linear or quadratic) over a small time-window under consideration.

[00146] For example may be modeled as: where e is some small positive constant.

[00147] A description will now be given regarding prediction based on the explicit solution to a continuous-time state-space model, in accordance with an exemplary aspect.

[00148] While the battery is in use, the SoP algorithm has to be capable of providing, at any moment in time, an estimate of the maximal power that can be put into or drawn from the battery in some pre-defined small time window. Therefore, in order to accurately estimate P(t) it is necessary to have an accurate prediction of the cell terminal voltage waveform over the time-horizon The battery performance can be evaluated by the discharge or charge energy E(t) while maintaining an SoP. The energy is expressed as follows:

[00149] The energy can be approximated using the composite Simpson’s rule

Riemann sum shown below where , , , , even subintervals and the ^^-th time step is [00150] A description will now be given regarding a conservative SoP Estimation, in accordance with an exemplary aspect. [00151] Applying a constant discharge or charge at does not result in constant power because of varying battery voltage. A constant power discharge or charge using the average power as the SoP might not be possible since voltages near will require higher currents than which can be constrained to the manufacturers’ limits. The peak power is described in (45) which is an underestimate of the maximum power. (45) [00152] A constant current pulse would imply an infinite current change rate at the beginning of the pulse. This can be considered the worst-case scenario and is used because of the following reasons: [00153] - Current can be easily controlled by power electronics. [00154] - Voltage change during constant current pulses are negligible. [00155] A description will now be given regarding a Newton-Raphson Method to find ^^ ^^ ^^ ^^ for the Conservative SoP Estimate, in accordance with an exemplary aspect. [00156] The max cell current results in the maximal power that can be drawn from the battery and is constrained by the minimum terminal voltage, The problem statement becomes: [00157] Problem: Solve for such that over the time-horizon [00158] The Newton-Raphson method is a powerful technique to solve root- finding problems numerically using linear approximations. With 2 initial guesses for ^^ ^^ ^^ ^^ being ^^ 0 and ^^ 1 the secant method is used to approximate the derivative and follows (46)until the tolerance for termination is satisfied for

[00159] The function returns the final terminal voltage at t + T p . This functions uses the state space model from (56) without observations with constant current i such that it is predicting the state for time T p .

[00160] The main issues with using this method include:

[00161] Bad initial guesses where the function doesn’t converge;

[00162] Sub-quadratic convergence rate; and

[00163] Solutions with noisy voltages that can reach below in previous time steps t + T p .

[00164] Possible improvements:

[00165] Use Brent’ s Method that combines the bisection method, the secant method and inverse quadratic interpolation to bound current within an interval;

[00166] - Log updates to restrict the current to be non-negative; and

[00167] Use a termination condition for instead of i n for steep slopes

[00168] A description will now be given of a constant power method to find in accordance with an exemplary aspect.

[00169] The max constant power results in more energy output than the conservative estimate for SoP. The maximum constant power discharge that results in the battery specified minimum terminal voltage at t + T p is found numerically utilizing Powell’s method. Batteries have maximum current discharge and charge pulses that need to be considered to restrict the current profile. This is accomplished by solving for the current at t + T p that is equal to the maximum discharge pulse current up to T p if the current reaches higher than the specification to reach the cut-off voltage. However, searching for this current profile that results in a constant power discharge is required with changing voltage, temperature, and SoC.

[00170] A description will now be given regarding a constant power simulation, in accordance with an exemplary aspect.

[00171] To find the max constant power discharge or charge a current profile needs to be solved to compensate for the changing voltage of the model. The current at each time step is computed numerically using Powell’s methods. With shorthand form

[00172] The current profile can be derived using 26 and 33 at the next time step n. [00173] V ocv is dependent on SoC and temperature. To solve for the current i[n] that results in a constant power discharge, Vocv [«] can be expressed in terms of i[n] and V ocv [n — 1] using the OCV-SoC inverse function relationship. However, since the OCV-SOC function is not defined yet an approximation of V ocv [n] « V ocv [n — 1] can be made with a small enough time step since the SoC change is insignificant. The current i[n] can be solved directly using (51) with the quadratic equation with a negative sign for positive overall cell voltage.

[00174] Results show that with a small enough time step (~ 100 ms) the desired power matches very closely to the calculated power. This function could be useful for more efficient computing of the SoP on an embedded system.

[00175] A description will now be given regarding tracking, in accordance with an exemplary aspect. The description will address a state space model and EKF approximation.

[00176] A description of a state space model will now be given, in accordance with an exemplary aspect.

[00177] The full state-space model corresponding to the backward Euler discretization for 2 RC pairs has state transition equation: which can be abbreviate to the following: for J x1 state vector x, J × J state transition matrix A , and J x1 control input vectob b. The following shorthand is used: [00178] The measurement equation is as follows: (56) which is abbreviated to the following: (57) for non-linear function and control input scalar d. The dependence of S on the OCV at time step n has been included. The OCV-SoC function can be estimated empirically and approximated with a 9 th -order polynomial, for example, to allow for deriving an EKF approximation. [00179] The presence of the control input ^^[ ^^] translates to very simple adjustments to the resulting Kalman filter’s predict and correct steps. [00180] A description will now be given of an EKF approximation, in accordance with an exemplary aspect. [00181] The Extended Kalman Filter (EKF) is a useful and simple approach for approximating non-linearities in a dynamical model. Given the state space model defined above, the following generative model is assumed: with noise processes: [00182] The corresponding EKF equations are as follows and as shown in FIG.13: [00183] Predict [00184] Linearize [00185] Correct where use is made of the Jacobian of the state-dependent measurement function.

[00186] The system matrices are expressed as a function of step index because they depend on the ECM parameters, which themselves are functions of SoC and temperature. SoC is tracked in the model and temperature is provided externally, so it is assumed that the non-linearity introduced by that interdependence is negligible. The system matrices also depend on the step duration, which can vary over time.

[00187] When no measurements are available, the EKF will only perform the predict step. This allows it to continue “tracking” when measurements are not available.

[00188] Additional Aspects

[00189] The present disclosure may additionally include one or more of the following aspects.

[00190] Aspect 1. A method for pretraining a hyper model configured for use in predicting a state of power (SoP) of a vehicle battery, the method comprising: performing electrochemical impedance spectroscopy (EIS) scans on a plurality of batteries having a set of similar operating characteristics to the vehicle battery, the EIS scans performed across various states of the vehicle battery; and fitting parameters of the hyper model by applying an optimization technique to results of the EIS scans, the hyper model comprising a family of models that each define a voltage response of a respective cell from among a plurality of cells of the vehicle battery to a current profile over the various states of the vehicle battery. [00191] Aspect 2. The method according to aspect 1, further comprising estimating the state of power (SoP) of the vehicle battery using the hyper model, the results of the EIS scans, operational constraints of the vehicle battery, and the various states of the vehicle battery used in the EIS scans. [00192] Aspect 3. The method according to aspect 1, wherein the various states of the vehicle battery include at least some of: different temperature ranges of the vehicle battery; different states of charge (SoCs) of the vehicle battery; age of the vehicle battery; and a nature of a current load the vehicle battery is subjected to. [00193] Aspect 4. The method according to aspect 1, wherein the hyper model is an equivalent circuit model (ECM) that maps various circuit elements such as resistors, capacitors, inductors, and Warburg impedance of the ECM to the various states of the vehicle battery used in the EIS scans. [00194] Aspect 5. The method according to aspect 1, wherein the hyper model is a model reduced from a physics based model that maps physics representations of battery elements to the various states of the vehicle battery used in the EIS scans. [00195] Aspect 6. The method according to aspect 1, wherein the hyper model is an adaptive filter inferred using a frequency response of an equivalent impedance of the vehicle battery learned using the EIS scans under the various states of the vehicle battery. [00196] Aspect 7. The method according to any of aspects 1-6, wherein the fitting is performed offline. [00197] Aspect 8. The method according to any of aspects 1-7, wherein a multidimensional look up table is used to map the hyper model, an expected battery state, the results of the EIS scans, and battery system current and voltage limitations into a corresponding SoP value. [00198] Aspect 9. The method according to aspect 1, wherein a mapping from the EIS scans, the various states of the vehicle battery, and operational constraints of the vehicle battery to the state of power (SoP) of the battery is performed using a regression function learned under the various states of the vehicle battery. [00199] Aspect 10. The method according to aspect 1, wherein the results of the EIS scans are used to infer a hyper model using an optimization technique selected from the group consisting of a gradient-based linear optimization method, a non-negative least squares (NNLS) method, and a convex optimization method. [00200] Aspect 11. The method according to aspect 1, wherein the hyper model is a family of equivalent circuit models (ECMs), and the method further comprises performing a smart initialization of a parameter fitting method by setting a series resistor R in an R—RC ECM model having the series resistor R in series with one or more RC parallel sub-circuits to a smallest observed impedance value, setting model parameters to determined values and holding the determined values fixed while scanning over a range based on a value of the series resistor R to identify a determined value that minimizes an objective function. [00201] Aspect 12. A method for predicting a state of power (SoP) of a battery, the method comprising: performing a plurality of electrochemical impedance spectroscopy (EIS) scans on the battery prior to an initial use of the battery in a vehicle; calibrating a pretrained hyper model using results of the plurality of EIS scans, the hyper model comprising a family of models that each define a voltage response of a respective cell from among a plurality of cells of the battery to a current profile over various states of the batterys; performing a plurality of additional EIS scans on the battery subsequent to the initial use of the battery in the vehicle; recalibrating the pretrained hyper model using results of the plurality of additional EIS scans; predicting the SoP of each of multiple cells of the vehicle battery responsive to a current battery state of each of the multiple cells; combining the SoP of each of multiple cells into a battery SoP; and controlling an amount of current extracted from or put into the battery responsive to a SoP value. [00202] Aspect 13. The method according to aspect 12, wherein the plurality of additional EIS scans are periodically performed subsequent to the initial use of the battery in the vehicle. [00203] Aspect 14. The method according to aspect 12, wherein the various states of the battery include at least some of different temperature ranges of the battery, different states of charge (SoCs) of the battery, age of the battery, and a nature of a current load the battery is subjected to. [00204] Aspect 15. The method according to aspect 12, wherein the hyper model is an equivalent circuit model (ECM) that maps various circuit elements such as resistors, capacitors, inductors, and Warburg impedance of the ECM to the various states of the battery used in the EIS scans. [00205] Aspect 16. The method according to aspect 12, wherein the hyper model is a model reduced from a physics based model that maps physics representations of battery elements to the various states of the battery used in the EIS scans. [00206] Aspect 17. The method according to aspect 12, wherein the hyper model is an adaptive filter inferred using a frequency response of an equivalent impedance of the battery learned using the EIS scans under the various states of the battery used in the EIS scans. [00207] Aspect 18. The method according to aspect 12, wherein the results of the EIS scans are used to infer a hyper model using an optimization technique selected from the group consisting of a gradient-based linear optimization method, a non-negative least squares (NNLS) method, and a convex optimization method. [00208] Aspect 19. The method according to aspect 12, wherein the hyper model is a family of equivalent circuit models (ECMs), and the method further comprises performing a smart initialization of a parameter fitting method by setting a series resistor R in an R—RC ECM model having the series resistor R in series with one or more RC parallel sub-circuits to a smallest observed impedance value, setting model parameters to determined values and holding the determined values fixed while scanning over a range based on a value of the series resistor R to identify a determined value that minimizes an objective function. [00209] Aspect 20. The method according to aspect 12, wherein the plurality of additional EIS scans is performed while the battery is in a key-on condition. [00210] Aspect 21. The method according to aspects 12 or 20, further comprising, during a key-on condition of the battery: measuring a current output from each of multiple cells of the battery; and determining the current battery state of each of the multiple cells of the battery responsive to the current output from each of the multiple cells of the battery. [00211] Aspect 22. The method according to aspect 12, further comprising predicting the SoP of at least one of the cells of the battery, wherein at least one constraint comprises at least one of a terminal voltage of the battery, a current of the battery, a temperature of the battery, and a state of charge (SoC) of the battery. [00212] Aspect 23. The method according to aspect 12, wherein the plurality of additional EIS scans is performed while the battery is in a key-off condition to update the hyper model. [00213] Aspect 24. The method according to aspect 12, wherein in a key-off condition, the battery SoP is estimated using a regression algorithm given EIS measurements and a current state of the battery. [00214] Aspect 25. The method according to aspect 12, wherein the SoP comprises a maximum allowable static current that can be sustained for a given time period. [00215] Aspect 26. The method according to aspect 25, wherein the maximum allowable static current comprises a level of current that does not cause a constraint violation. [00216] Aspect 27. The method according to aspect 25, wherein the maximum allowable static current comprises a level of current that does not result in a temporary loss of capacity of more than a specified fraction over a given time period. [00217] Aspect 28. The method according to aspect 12, wherein the SoP comprises a maximum allowable constant power that can be sustained for a given time period. [00218] Aspect 29. The method according to aspect 20, further comprising tracking a battery state starting from a known initial condition responsive to measurable inputs comprising the results of the EIS scans. [00219] Aspect 30. The method according to aspect 20, further comprising tracking a battery state including a terminal voltage of the battery using an Extended Kalman filter (EKF) technique that comprises predicting a battery state vector using a state space model and predicting the terminal voltage of the battery response over a desired period of time. [00220] Aspect 31. The method according to aspect 30, wherein the EKF technique comprises forming battery system matrices expressed as a function of a step index that is dependent on parameters of the hyper model that, in turn, are dependent on the various states of the battery. [00221] Aspect 32. The method according to aspect 20, further comprising estimating the battery state including a terminal voltage of the battery using a closed form solution to discrete approximations of a state space model and predicting the terminal voltage of the battery response over a desired period of time. [00222] Aspect 33. The method according to aspect 20, further comprising generating a prediction for a current SoP of the battery responsive to a present prediction of the battery state, measurable inputs, and system limitations relating to current and voltage maximum and minimum values. [00223] Aspect 34. The method according to aspect 33, wherein the battery state comprises a present battery current, a present terminal voltage, a present open circuit voltage (OCV), and a present impedance measurement obtained using EIS. [00224] Aspect 35. The method according to aspect 12, wherein, for a recalibration, respective complexities of models in the family of models comprised in the hyper model increase with increasing battery cell age. [00225] Aspect 36. The method according to aspect 12, wherein the SoP of the battery is equal to one of the SoPs of each of multiple cells. [00226] Aspect 37. The method according to aspect 12, wherein the SoP value of a battery is equal to a lowest value from among the multiple cell SoPs for the multiple cells that constitute the battery.

[00227] Aspect 38. The method according to aspect 12, wherein the SoP value is equal to the battery SoP.

[00228] Aspect 39: A system for predicting a state of power (SoP) of a battery, the system comprising: an electrochemical impedance spectroscopy (EIS) system for performing a plurality of EIS scans on the battery prior to an initial use of the battery in a vehicle, and a plurality of additional EIS scans on the battery in the vehicle; a memory device for storing program code; and a processing device operatively coupled to the EIS system and the memory device for running the program code to: calibrate a pretrained hyper model using results of the plurality of EIS scans, the pretrained hyper model comprising a family of models that each define a voltage response of a respective cell from among a plurality of cells of the battery to a current profile over various states of the battery; recalibrate the pretrained hyper model using results of the plurality of additional EIS scans; predict the SoP of each of multiple cells of the vehicle battery responsive to a current battery state of each of the multiple cells; combine the SoPs of each of multiple cells into a battery SoP; and control an amount of current extracted from or put into the battery responsive to at least one of the SoPs of each of multiple cells or the battery SoP.

[00229] Aspect 40. A system having one or more components configured to perform the function of any of aspects 13 to 38.

[00230] It should be noted that all of the specifications, dimensions, and relationships outlined herein (e.g., the number of elements, operations, steps, etc.) have only been offered for purposes of example and teaching only. Such information may be varied considerably without departing from the spirit of the present disclosure, or the scope of the appended claims. The specifications apply only to one non-limiting example and, accordingly, they should be construed as such. In the foregoing description, exemplary aspects have been described with reference to particular component arrangements. Various modifications and changes may be made to such aspects without departing from the scope of the appended claims. The description and drawings are, accordingly, to be regarded in an illustrative rather than in a restrictive sense. [00231] Note that with the numerous examples provided herein, interaction may be described in terms of two, three, four, or more electrical components. However, this has been done for purposes of clarity and example only. It should be appreciated that the system may be consolidated in any suitable manner. Along similar design alternatives, any of the illustrated components, modules, and elements of the FIGURES may be combined in various possible configurations, all of which are clearly within the broad scope of this Specification. In certain cases, it may be easier to describe one or more of the functionalities of a given set of flows by only referencing a limited number of electrical elements. It should be appreciated that the electrical circuits of the FIGURES and its teachings are readily scalable and may accommodate a large number of components, as well as more complicated/sophisticated arrangements and configurations. Accordingly, the examples provided should not limit the scope or inhibit the broad teachings of the electrical circuits as potentially applied to myriad other architectures. [00232] It should also be noted that in this Specification, references to various features (e.g., elements, structures, modules, components, steps, operations, characteristics, etc.) included in "one aspect", "exemplary aspect", "an aspect", "another aspect", "some aspects", "various aspects", "other aspects", "alternative aspect", and the like are intended to mean that any such features are included in one or more aspects of the present disclosure, but may or may not necessarily be combined in the same aspects. [00233] It is to be appreciated that the use of any of the following “/”, “and/or”, and “at least one of”, for example, in the cases of “A/B”, “A and/or B” and “at least one of A and B”, is intended to encompass the selection of the first listed option (A) only, or the selection of the second listed option (B) only, or the selection of both options (A and B). As a further example, in the cases of “A, B, and/or C” and “at least one of A, B, and C”, such phrasing is intended to encompass the selection of the first listed option (A) only, or the selection of the second listed option (B) only or the selection of the third listed option (C) only, or the selection of the first and the second listed options (A and B) only, or the selection of the first and third listed options (A and C) only, or the selection of the second and third listed options (B and C) only, or the selection of all three options (A and B and C). This may be extended, as readily apparent by one of ordinary skill in this and related arts, for as many items listed. [00234] It should also be noted that the functions related to circuit architectures illustrate only some of the possible circuit architecture functions that may be executed by, or within, systems illustrated in the FIGURES. Some of these operations may be deleted or removed where appropriate, or these operations may be modified or changed considerably without departing from the scope of the present disclosure. In addition, the timing of these operations may be altered considerably. The preceding operational flows have been offered for purposes of example and discussion. Substantial flexibility is provided by aspects described herein in that any suitable arrangements, chronologies, configurations, and timing mechanisms may be provided without departing from the teachings of the present disclosure. [00235] Numerous other changes, substitutions, variations, alterations, and modifications may be ascertained to one skilled in the art and it is intended that the present disclosure encompass all such changes, substitutions, variations, alterations, and modifications as falling within the scope of the appended claims. [00236] Note that all optional features of the device and system described above may also be implemented with respect to the method or process described herein and specifics in the examples may be used anywhere in one or more aspects. [00237] The “means for” in these instances (above) may include (but is not limited to) using any suitable component discussed herein, along with any suitable software, circuitry, hub, computer code, logic, algorithms, hardware, controller, interface, link, bus, communication pathway, etc. [00238] Note that with the example provided above, as well as numerous other examples provided herein, interaction may be described in terms of two, three, or four network elements. However, this has been done for purposes of clarity and example only. In certain cases, it may be easier to describe one or more of the functionalities of a given set of flows by only referencing a limited number of network elements. It should be appreciated that topologies illustrated in and described with reference to the accompanying FIGURES (and their teachings) are readily scalable and may accommodate a large number of components, as well as more complicated/sophisticated arrangements and configurations. Accordingly, the examples provided should not limit the scope or inhibit the broad teachings of the illustrated topologies as potentially applied to myriad other architectures. [00239] It is also important to note that the steps in the preceding flow diagrams illustrate only some of the possible signaling scenarios and patterns that may be executed by, or within, communication systems shown in the FIGURES. Some of these steps may be deleted or removed where appropriate, or these steps may be modified or changed considerably without departing from the scope of the present disclosure. In addition, a number of these operations have been described as being executed concurrently with, or in parallel to, one or more additional operations. However, the timing of these operations may be altered considerably. The preceding operational flows have been offered for purposes of example and discussion. Substantial flexibility is provided by communication systems shown in the FIGURES in that any suitable arrangements, chronologies, configurations, and timing mechanisms may be provided without departing from the teachings of the present disclosure. [00240] Although the present disclosure has been described in detail with reference to particular arrangements and configurations, these example configurations and arrangements may be changed significantly without departing from the scope of the present disclosure. For example, although the present disclosure has been described with reference to particular communication exchanges, aspects described herein may be applicable to other architectures. [00241] Numerous other changes, substitutions, variations, alterations, and modifications may be ascertained to one skilled in the art and it is intended that the present disclosure encompass all such changes, substitutions, variations, alterations, and modifications as falling within the scope of the appended claims. In order to assist the United States Patent and Trademark Office (USPTO) and, additionally, any readers of any patent issued on this application in interpreting the claims appended hereto, Applicant wishes to note that the Applicant: (a) does not intend any of the appended claims to invoke paragraph six (6) of 35 U.S.C. section 142 as it exists on the date of the filing hereof unless the words “means for” or “step for” are specifically used in the particular claims; and (b) does not intend, by any statement in the specification, to limit this disclosure in any way that is not otherwise reflected in the appended claims.