CN107219466A - 一种混合扩展卡尔曼滤波的锂电池soc估算方法 - Google Patents
一种混合扩展卡尔曼滤波的锂电池soc估算方法 Download PDFInfo
- Publication number
- CN107219466A CN107219466A CN201710436071.XA CN201710436071A CN107219466A CN 107219466 A CN107219466 A CN 107219466A CN 201710436071 A CN201710436071 A CN 201710436071A CN 107219466 A CN107219466 A CN 107219466A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mover
- mtd
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L58/00—Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles
- B60L58/10—Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries
- B60L58/12—Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries responding to state of charge [SoC]
- B60L58/13—Maintaining the SoC within a determined range
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/392—Determining battery ageing or deterioration, e.g. state of health
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/70—Energy storage systems for electromobility, e.g. batteries
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Sustainable Energy (AREA)
- Sustainable Development (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Secondary Cells (AREA)
Abstract
一种混合扩展卡尔曼滤波的锂电池SOC估算方法,包括如下步骤:建立锂电池等效模型:锂电池等效模型,电池的开路电压OCV表示为电化学模型,R指电池内阻,该值在充电和放电时的值是不同的,充电时将其设为R+,放电时设为R‑,p1、p2、p3、p4、p5、p6、p7为符合锂电池模型的拟合参数,根据电池充放电实验数据,可以估算出电池等效电路模型参数;建立电池***离散的状态空间模型,采用改进的EKF算法对电池SOC进行估算。本发明克服了EKF算法中由于电池模型的不确定性、噪声统计特性未知及泰勒展开线性化可能造成滤波精度降低甚至产生滤波发散的问题。
Description
【技术领域】
本发明属于锂电池技术领域,具体涉及一种混合扩展卡尔曼滤波的锂电池SOC估算方法。
【背景技术】
电池是电动汽车的主要能量载体和动力来源,也是电动汽车整车车体的主要组成部分,准确的估算电池SOC,不仅可以提高电池的容量利用效率,还可以延长电池使用寿命。由于电池内部复杂的电化学性质及物理反应,SOC不能直接测量,一般是通过电池外部参数,如电压及工作电流等进行估算。
SOC是电池管理***最重要的参数之一,准确估计SOC不仅为电动汽车驾驶员提供准确的剩余电量,也为电池管理***的管理和控制提供依据。目前,最为常用的SOC估算算法通常为卡尔曼滤波算法。传统的卡尔曼滤波是在标准条件下获得的,是一种无偏的线性最小方差估计算法,在动态***的数学模型和噪声统计特性已知的情况下,卡尔曼滤波通过测量值对预测估计进行修正,可以得到状态的精确估计。但是在实际应用中很难得到电池精确的数学模型和噪声的统计特性,使滤波精度降低甚至产生滤波发散的现象。
中国发明专利CN 103744027A公开了一种基于卡尔曼滤波的自校正电池SOC。算法首先初始化状态真实值,状态估计值,状态真实值可根据开路电压法确定;然后利用安时积分法的SOC计算式,作为卡尔曼滤波算法中状态方程,计算Xk+1;然后观察方程计算真实值Yk+1;最后通过卡尔曼递推原理,不断对当前时刻SOC值进行修正。该发明具有下列优点:(1)有效解决开路电压法、安时积分法等传统单一方法估算SOC存在的弊端,(2)不断对电池总容量和可用容量进行校正,可有效弥补SOC估算中初始误差和累积误差的优点。由于在实际应用中很难得到***精确的数学模型和噪声的统计特性,这就可能使滤波精度降低甚至产生滤波发散的现象。
中国发明专利CN105510829 A公开了一种新型锂离子动力电池SOC估计方法。通过采用强跟踪滤波器估计SOC,克服了扩展卡尔曼滤波器由于模型不确定性造成SOC估计不准的缺点,强跟踪滤波器由扩展卡尔曼滤波器改造而来,主要针对***模型不确定性导致滤波器估计不准及发散问题,具有以下的优点:(1)对模型不确定性具有较强的鲁棒性;(2)对突变状态的跟踪能力极强,甚至在***达到平衡状态时,仍保持对缓变状态与突变状态的跟踪能力;(3)适中的计算复杂度。由于强跟踪算法要求测量噪声统计是已知的,但大多数应用问题中噪声统计特性是未知的,用错误的噪声统计设计滤波器将使滤波器的误差增大。
中国发明专利CN105093122 A公开了一种基于强跟踪自适应SQKF的应急灯电池SOC估计方法。通过建立应急灯电池的二阶RC等效模型;其次在电池等效模型的基础上,建立相应的离散状态空间模型方程;最后采用强跟踪自适应SQKF对电池的SOC值进行滤波估计。该发明具有下列优点:(1)比现有卡尔曼滤波框架下的SOC估计方法更高的估计精度。(2)通过引入时变渐消因子和在线估计***噪声方差,有效抑制了电池***建模时由于噪声统计特性未知时变而导致的滤波发散问题。由于每次滤波都要估噪声R,增加了滤波的复杂度。
目前SOC估算方法主要存在以下问题:
1、模型过于简化或***参数发生变动。实际中的***模型一般比较复杂,若要精确的描述***,模型中的状态变量需要达到较高的维数,这不利于对***状态进行重构,因此大多数情况下都要使用模型简化的办法,即忽略***中某些不重要的因素,从而使用相对较少的状态变量来描述***的主要特征,这就可能造成应用模型与实际***之前的不匹配。***在实际运行过程中可能会出现元件损坏、老化问题,这些问题均会对***模型造成影响,进而使模型参数发生变化,与原模型的匹配程度降低。
2、噪声统计特性不准确。在模型中都需要考虑***噪声和测量噪声,在大多数模型中应用的噪声统计特性都比较理想,从而导致实际过程噪声的统计特性与理论特性不能达到一致性。在实际情况中,***的统计特性可能因为***干扰而发生变化,这将会造成噪声统计特性的不准确性。
3、泰勒展开线性化误差。当线性化展开点为***状态的真值时,线性化误差最小。但实际运用中,***状态真值是未知的,因此必定存在线性化误差,且线性化展开点与真值相差越远,误差越大,从而会导致滤波结果不稳定。
【发明内容】
本发明提出了一种混合扩展卡尔曼滤波的锂电池SOC估算方法,克服了EKF算法中由于电池模型的不确定性、噪声统计特性未知及泰勒展开线性化可能造成滤波精度降低甚至产生滤波发散的问题。
本发明是这样实现的:
一种混合扩展卡尔曼滤波的锂电池SOC估算方法,包括如下步骤:
步骤一:建立锂电池等效模型:
锂电池等效模型,电池的开路电压OCV表示为电化学模型,R指电池内阻,该值在充电和放电时的值是不同的,充电时将其设为R+,放电时设为R-,p1、p2、p3、p4、p5、p6、p7为符合锂电池模型的拟合参数,根据电池充放电实验数据,可以估算出电池等效电路模型参数;
步骤二:建立电池***离散的状态空间模型:
Zk=OCV(k)+i(t)R=f(SOC(k))=g(xk,uK)+vk (2)
式中Cn是额定容量;i为电池电流;η为库伦效率,其中E[Wk]=0, E[vk]=0,为***的随机干扰,Qk为过程激励噪声协方差矩阵,Rk为观测噪声协方差矩阵。f(wk.,uK)是一个非线性状态转换函数,g(xk,uk)是一个非线性测量函数;
定义
步骤三:采用改进的EKF算法对电池SOC进行估算;
1)k=0,选择初始值P0/0
2)将式(1)在处一阶Taylor展开
3)计算残差及时变渐消因子λk+1:
残差:
式中γk为k时刻残差,γ(1)表示k=0时刻的残差,Sk表示残差的协方差矩阵,0<ρ<1为遗忘因子,通常在0.95-0.99之间选取,β≥1为弱化因子,目的是使状态估计值更加平滑;
4)计算预测值
5)预测误差协方差矩阵
6)将式(2)在处在处一阶Taylor展开,有
7)计算卡尔曼滤波增益矩阵
8)计算滤波估计值
9)计算后向平滑值
10)计算SOC更新值及量测噪声,即
量测噪声计算:
首先判断滤波是否异常:
γ(k)γT(k)>κtr[γ(k)γT(k)] (13)
利用式(13)对滤波的状态进判断,若式(13)成立,则说明出现滤波异常,原来的动力学模型已不适应当前的滤波,应采用更加重视当前观测数据的方法来使其适应模型变化对滤波产生的影响,需要对噪声R(k)进行估计,使其适应当前的滤波;反之,式(18)成立,说明滤波无异常,不需要估计R(k),其中R(k)的计算公式如下:
dk=(1-b)/(1-bk+1)
R(k)=(1-dk)Rk-1+dk[γk+1γT k+1HkPk-1HT k]
式中,κ为储备系数,κ>1;0<b<1为遗忘因子,通常在0.95-0.99之间选取;tr表示矩阵的迹;γ(k)为新息系列;
11)判断滤波是否执行,如果是,返回步骤2);否则,结束算法。
本发明的优点在于:考虑到电池模型、算法中泰勒展开线性化及噪声统计特性对状态估算的影响,为了使电池的状态估计更加准确,提出了一种集成STF与Sage-Husa自适应滤波的混合EKF算法,在EKF算法的基础上,首先对观测方程线性化方式进行改进,考虑线性化展开点和雅可比矩阵取值点对线性化逼近的影响,采用基于中值定理的线性化处理方式,有效降低了非线性函数的线性化误差;并通过强跟踪滤波算法原理,引入渐消因子λk+1,使得残差序列在每一步互相正交,减少了EKF算法对电池模型不确定性的影响;根据简化的Sage-Husa自适应滤波算法估算噪声原理,采用协方差匹配技术,确定量测噪声协方差的值,减少了对噪声估算的运算量。与传统的EKF算法相比,该算法降低了EKF算法中泰勒展开线性化可能造成的误差,克服了电池初始SOC值及EKF算法中由于电池模型的不确定性及噪声统计特性未知,而造成的滤波精度降低的缺点,算法在复杂度上会略有增加,但滤波精度、稳定性及实时性有所提高。
【附图说明】
下面参照附图结合实施例对本发明作进一步的描述。
图1是本发明的电池模型图。
图2是本发明的方法流程示意图。
【具体实施方式】
一种混合扩展卡尔曼滤波的锂电池SOC估算方法,包括如下步骤:
步骤一:建立锂电池等效模型:
锂电池等效模型,电池的开路电压OCV表示为电化学模型,R指电池内阻,该值在充电和放电时的值是不同的,充电时将其设为R+,放电时设为R-,p1、p2、p3、p4、p5、p6、p7为符合锂电池模型的拟合参数,根据电池充放电实验数据,可以估算出电池等效电路模型参数,如图1所示;
步骤二:建立电池***离散的状态空间模型:
Zk=OCV(k)+i(t)R=f(SOC(k))=g(xk,uK)+vk (2)
式中Cn是额定容量;i为电池电流;η为库伦效率,其中E[Wk]=0, E[vk]=0,为***的随机干扰,Qk为过程激励噪声协方差矩阵,Rk为观测噪声协方差矩阵。f(wk.,uK)是一个非线性状态转换函数,g(xk,uk)是一个非线性测量函数;
定义
步骤三:采用改进的EKF算法对电池SOC进行估算;
1)k=0,选择初始值P0/0
2)将式(1)在处一阶Taylor展开
3)计算残差及时变渐消因子λk+1:
残差:
式中γk为k时刻残差,γ(1)表示k=0时刻的残差,Sk表示残差的协方差矩阵,0<ρ<1为遗忘因子,通常在0.95-0.99之间选取,β≥1为弱化因子,目的是使状态估计值更加平滑;
4)计算预测值
5)预测误差协方差矩阵
6)将式(2)在处在处一阶Taylor展开,有
7)计算卡尔曼滤波增益矩阵
8)计算滤波估计值
9)计算后向平滑值
10)计算SOC更新值及量测噪声,即
量测噪声计算:
首先判断滤波是否异常:
γ(k)γT(k)>κtr[γ(k)γT(k)] (13)
利用式(13)对滤波的状态进判断,若式(13)成立,则说明出现滤波异常,原来的动力学模型已不适应当前的滤波,应采用更加重视当前观测数据的方法来使其适应模型变化对滤波产生的影响,需要对噪声R(k)进行估计,使其适应当前的滤波;反之,式(18)成立,说明滤波无异常,不需要估计R(k),其中R(k)的计算公式如下:
dk=(1-b)/(1-bk+1)
R(k)=(1-dk)Rk-1+dk[γk+1γT k+1HkPk-1HT k] (14)
式中,κ为储备系数,κ>1;0<b<1为遗忘因子,通常在0.95-0.99之间选取;tr表示矩阵的迹;γ(k)为新息系列;
11)判断滤波是否执行,如果是,返回步骤2);否则,结束算法。
考虑到电池模型、算法中泰勒展开线性化及噪声统计特性对状态估算的影响,为了使电池的状态估计更加准确,本发明提出了一种集成强跟踪滤波与Sage-Husa自适应滤波的混合EKF算法,在EKF算法的基础上,首先对观测方程线性化方式进行改进,考虑线性化展开点和雅可比矩阵取值点对线性化逼近的影响,采用基于中值定理的线性化处理方式,有效降低了非线性函数的线性化误差;并通过强跟踪滤波算法原理,引入渐消因子λk+1,使得残差序列在每一步互相正交,减少了EKF算法对电池模型不确定性的影响;根据简化的Sage-Husa自适应滤波算法估算噪声原理,采用协方差匹配技术,确定量测噪声协方差的值,减少了对噪声估算的运算量。与传统的EKF算法相比,该算法降低了EKF算法中泰勒展开线性化可能造成的误差,克服了电池初始SOC值及EKF算法中由于电池模型的不确定性及噪声统计特性未知,而造成的滤波精度降低的缺点,算法在复杂度上会略有增加,但滤波精度、稳定性及实时性有所提高。
以上所述仅为本发明的较佳实施用例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换以及改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种混合扩展卡尔曼滤波的锂电池SOC估算方法,其特征在于:包括如下步骤:
步骤一:建立锂电池等效模型:
锂电池等效模型,电池的开路电压OCV表示为电化学模型,R指电池内阻,该值在充电和放电时的值是不同的,充电时将其设为R+,放电时设为R-,p1、p2、p3、p4、p5、p6、p7为符合锂电池模型的拟合参数,根据电池充放电实验数据,可以估算出电池等效电路模型参数;
步骤二:建立电池***离散的状态空间模型:
(1)
(2) Zk=OCV(k)+i(t)R=f(SOC(k))=g(xk,uK)+vk
式中Cn是额定容量;i为电池电流;η为库伦效率,其中E[Wk]=0,E[vk]=0,E[vk vT j]=Rkδk,j为***的随机干扰,Qk为过程激励噪声协方差矩阵,Rk为观测噪声协方差矩阵。f(wk.,uK)是一个非线性状态转换函数,g(xk,uk)是一个非线性测量函数;
定义
步骤三:采用改进的EKF算法对电池SOC进行估算;
1)k=0,选择初始值P0/0
2)将式(1)在处一阶Taylor展开
<mrow>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&Gamma;w</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
3)计算残差及时变渐消因子λk+1:
残差:
<mrow>
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
<mo>&GreaterEqual;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
<mo><</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mrow>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>&rho;S</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>&gamma;</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&rho;</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>&gamma;</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>&GreaterEqual;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>S</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<msub>
<mi>Q</mi>
<mi>k</mi>
</msub>
<msub>
<msup>
<mi>H</mi>
<mi>T</mi>
</msup>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;R</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<msub>
<mi>F</mi>
<mi>k</mi>
</msub>
<msub>
<mi>P</mi>
<mi>k</mi>
</msub>
<msub>
<msup>
<mi>F</mi>
<mi>T</mi>
</msup>
<mi>k</mi>
</msub>
<msub>
<msup>
<mi>H</mi>
<mi>T</mi>
</msup>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中γk为k时刻残差,γ(1)表示k=0时刻的残差,Sk表示残差的协方差矩阵,0<ρ<1为遗忘因子,通常在0.95-0.99之间选取,β≥1为弱化因子,目的是使状态估计值更加平滑;
4)计算预测值
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
5)预测误差协方差矩阵
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>|</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&lambda;</mi>
<mi>k</mi>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>F</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>Q&Gamma;</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
6)将式(2)在处在处一阶Taylor展开,有
<mrow>
<msub>
<mi>Z</mi>
<mi>k</mi>
</msub>
<mo>&ap;</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
7)计算卡尔曼滤波增益矩阵
<mrow>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
8)计算滤波估计值
<mrow>
<msub>
<mover>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>k</mi>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
9)计算后向平滑值
<mrow>
<msub>
<mover>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>F</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<msup>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
10)计算SOC更新值及量测噪声,即
<mrow>
<mfenced open = "" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<mo>=</mo>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>|</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<msub>
<mover>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>K</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<msub>
<msup>
<mi>P</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mover>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>^</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<msup>
<mi>K</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Z</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<msub>
<msup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>I</mi>
<mo>-</mo>
<msub>
<msup>
<mi>K</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<msub>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>I</mi>
<mo>-</mo>
<msub>
<msup>
<mi>K</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<msub>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>+</mo>
<msub>
<msup>
<mi>K</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
</msub>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<msubsup>
<msup>
<mi>K</mi>
<mo>&prime;</mo>
</msup>
<mi>k</mi>
<mi>T</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
量测噪声计算:
首先判断滤波是否异常:
γ(k)γT(k)>κtr[γ(k)γT(k)] (13)
利用式(13)对滤波的状态进判断,若式(13)成立,则说明出现滤波异常,原来的动力学模型已不适应当前的滤波,应采用更加重视当前观测数据的方法来使其适应模型变化对滤波产生的影响,需要对噪声R(k)进行估计,使其适应当前的滤波;反之,式(18)成立,说明滤波无异常,不需要估计R(k),其中R(k)的计算公式如下:
dk=(1-b)/(1-bk+1)
R(k)=(1-dk)Rk-1+dk[γk+1γT k+1HkPk-1HT k] (14)
式中,κ为储备系数,κ>1;0<b<1为遗忘因子,通常在0.95-0.99之间选取;tr表示矩阵的迹;γ(k)为新息系列;
11)判断滤波是否执行,如果是,返回步骤2);否则,结束算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710436071.XA CN107219466A (zh) | 2017-06-12 | 2017-06-12 | 一种混合扩展卡尔曼滤波的锂电池soc估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710436071.XA CN107219466A (zh) | 2017-06-12 | 2017-06-12 | 一种混合扩展卡尔曼滤波的锂电池soc估算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107219466A true CN107219466A (zh) | 2017-09-29 |
Family
ID=59947773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710436071.XA Pending CN107219466A (zh) | 2017-06-12 | 2017-06-12 | 一种混合扩展卡尔曼滤波的锂电池soc估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107219466A (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
CN108594135A (zh) * | 2018-06-28 | 2018-09-28 | 南京理工大学 | 一种用于锂电池均衡充放电控制的soc估算方法 |
CN108732509A (zh) * | 2018-06-06 | 2018-11-02 | 哈尔滨工业大学 | 一种面向空间应用的锂离子电池荷电状态的在线估计方法 |
CN108872865A (zh) * | 2018-05-29 | 2018-11-23 | 太原理工大学 | 一种防滤波发散的锂电池soc估算方法 |
CN109000663A (zh) * | 2018-08-07 | 2018-12-14 | 中国人民解放军海军航空大学 | 基于时变自调整的滚动时域估计到达代价计算方法 |
CN109782181A (zh) * | 2018-12-20 | 2019-05-21 | 宁波飞拓电器有限公司 | 一种基于组合滤波的应急灯电池soc估计方法 |
CN109870651A (zh) * | 2019-01-22 | 2019-06-11 | 重庆邮电大学 | 一种电动汽车动力电池***soc和soh联合在线估算方法 |
CN109885849A (zh) * | 2018-05-07 | 2019-06-14 | 长春工业大学 | 基于强跟踪滤波的轨道客车微动开关剩余寿命预测方法 |
CN111007400A (zh) * | 2019-11-22 | 2020-04-14 | 西安工程大学 | 基于自适应双扩展卡尔曼滤波法的锂电池soc估算方法 |
CN111376793A (zh) * | 2018-12-29 | 2020-07-07 | 观致汽车有限公司 | 用于管理电池的方法、设备和计算机可读介质 |
CN111650517A (zh) * | 2020-04-20 | 2020-09-11 | 北京理工大学 | 一种电池荷电状态估算方法 |
CN111781503A (zh) * | 2020-06-15 | 2020-10-16 | 国网江苏省电力有限公司无锡供电分公司 | 一种锂离子储能电池soc在线估算方法 |
CN112003527A (zh) * | 2020-07-22 | 2020-11-27 | 西安理工大学 | 用于异步电机转速辨识的迭代扩展卡尔曼滤波的改进方法 |
CN113219830A (zh) * | 2021-05-06 | 2021-08-06 | 无锡灵鸽机械科技股份有限公司 | 一种失重式喂料机控制方法 |
CN113884914A (zh) * | 2021-10-26 | 2022-01-04 | 河南理工大学 | 一种渐消并行卡尔曼滤波动力电池荷电状态估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6534954B1 (en) * | 2002-01-10 | 2003-03-18 | Compact Power Inc. | Method and apparatus for a battery state of charge estimator |
CN103675683A (zh) * | 2012-09-02 | 2014-03-26 | 东莞市振华新能源科技有限公司 | 一种锂电池荷电状态(soc)估算方法 |
CN104217107A (zh) * | 2014-08-27 | 2014-12-17 | 华南理工大学 | 一种基于多传感器信息仿人机器人跌倒状态检测方法 |
CN105510829A (zh) * | 2014-09-29 | 2016-04-20 | 山东大学 | 一种新型锂离子动力电池soc估计方法 |
-
2017
- 2017-06-12 CN CN201710436071.XA patent/CN107219466A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6534954B1 (en) * | 2002-01-10 | 2003-03-18 | Compact Power Inc. | Method and apparatus for a battery state of charge estimator |
CN103675683A (zh) * | 2012-09-02 | 2014-03-26 | 东莞市振华新能源科技有限公司 | 一种锂电池荷电状态(soc)估算方法 |
CN104217107A (zh) * | 2014-08-27 | 2014-12-17 | 华南理工大学 | 一种基于多传感器信息仿人机器人跌倒状态检测方法 |
CN105510829A (zh) * | 2014-09-29 | 2016-04-20 | 山东大学 | 一种新型锂离子动力电池soc估计方法 |
Non-Patent Citations (2)
Title |
---|
刘浩: "基于EKF的电动汽车用锂离子电池SOC估算方法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
王晴晴等: "改进EKF算法在固定单站被动目标跟踪中的应用", 《空军预警学院学报》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
CN109885849A (zh) * | 2018-05-07 | 2019-06-14 | 长春工业大学 | 基于强跟踪滤波的轨道客车微动开关剩余寿命预测方法 |
CN108872865A (zh) * | 2018-05-29 | 2018-11-23 | 太原理工大学 | 一种防滤波发散的锂电池soc估算方法 |
CN108732509A (zh) * | 2018-06-06 | 2018-11-02 | 哈尔滨工业大学 | 一种面向空间应用的锂离子电池荷电状态的在线估计方法 |
CN108594135A (zh) * | 2018-06-28 | 2018-09-28 | 南京理工大学 | 一种用于锂电池均衡充放电控制的soc估算方法 |
CN109000663A (zh) * | 2018-08-07 | 2018-12-14 | 中国人民解放军海军航空大学 | 基于时变自调整的滚动时域估计到达代价计算方法 |
CN109782181A (zh) * | 2018-12-20 | 2019-05-21 | 宁波飞拓电器有限公司 | 一种基于组合滤波的应急灯电池soc估计方法 |
CN111376793A (zh) * | 2018-12-29 | 2020-07-07 | 观致汽车有限公司 | 用于管理电池的方法、设备和计算机可读介质 |
CN109870651A (zh) * | 2019-01-22 | 2019-06-11 | 重庆邮电大学 | 一种电动汽车动力电池***soc和soh联合在线估算方法 |
CN111007400A (zh) * | 2019-11-22 | 2020-04-14 | 西安工程大学 | 基于自适应双扩展卡尔曼滤波法的锂电池soc估算方法 |
CN111650517A (zh) * | 2020-04-20 | 2020-09-11 | 北京理工大学 | 一种电池荷电状态估算方法 |
CN111650517B (zh) * | 2020-04-20 | 2021-08-17 | 北京理工大学 | 一种电池荷电状态估算方法 |
CN111781503A (zh) * | 2020-06-15 | 2020-10-16 | 国网江苏省电力有限公司无锡供电分公司 | 一种锂离子储能电池soc在线估算方法 |
CN111781503B (zh) * | 2020-06-15 | 2023-09-22 | 国网江苏省电力有限公司无锡供电分公司 | 一种锂离子储能电池soc在线估算方法 |
CN112003527A (zh) * | 2020-07-22 | 2020-11-27 | 西安理工大学 | 用于异步电机转速辨识的迭代扩展卡尔曼滤波的改进方法 |
CN112003527B (zh) * | 2020-07-22 | 2022-08-16 | 西安理工大学 | 用于异步电机转速辨识的迭代扩展卡尔曼滤波的改进方法 |
CN113219830A (zh) * | 2021-05-06 | 2021-08-06 | 无锡灵鸽机械科技股份有限公司 | 一种失重式喂料机控制方法 |
CN113884914A (zh) * | 2021-10-26 | 2022-01-04 | 河南理工大学 | 一种渐消并行卡尔曼滤波动力电池荷电状态估计方法 |
CN113884914B (zh) * | 2021-10-26 | 2024-03-12 | 河南理工大学 | 一种渐消并行卡尔曼滤波动力电池荷电状态估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107219466A (zh) | 一种混合扩展卡尔曼滤波的锂电池soc估算方法 | |
CN104617623B (zh) | 一种电动汽车动力电池组均衡控制方法 | |
CN107368619B (zh) | 一种扩展卡尔曼滤波soc估算方法 | |
CN103852727B (zh) | 动力电池荷电状态在线估算方法和装置 | |
CN102998629B (zh) | 一种锂电池荷电状态的估计方法 | |
CN103020445B (zh) | 一种电动车车载磷酸铁锂电池的soc与soh预测方法 | |
EP2321663B1 (en) | Apparatus and method for estimating state of health of battery based on battery voltage variation pattern | |
CN109946623A (zh) | 一种锂电池的soc在线估测方法 | |
CN103616647B (zh) | 一种用于电动汽车电池管理***的电池剩余电量估计方法 | |
CN106909716A (zh) | 计及容量损耗的磷酸铁锂电池建模及soc估计方法 | |
CN103744030B (zh) | 电池组健康状态和荷电状态在线估算装置及估算方法 | |
CN101625397B (zh) | 一种电池剩余电量的混合快速估计方法 | |
CN109633479B (zh) | 基于嵌入式容积卡尔曼滤波的锂电池soc在线估算方法 | |
CN103797374A (zh) | 用于电池监控的***和方法 | |
CN104122504A (zh) | 一种电池的soc估算方法 | |
CN104678316A (zh) | 锂离子电池荷电状态估算方法和装置 | |
CN107153163A (zh) | 一种基于自适应ukf的锂电池soc估算方法 | |
CN106443478A (zh) | 基于闭环混合算法的磷酸铁锂电池剩余电量的估算方法 | |
CN102788957A (zh) | 一种动力电池荷电状态估算方法 | |
CN102645637A (zh) | 一种均衡电池的soc估算方法 | |
KR100901252B1 (ko) | 슬라이딩 모드 관측기를 이용한 2차 전지 soc 예측방법 및 장치 | |
CN105699910A (zh) | 一种锂电池剩余电量在线估计方法 | |
CN109342950A (zh) | 一种用于锂电池荷电状态的评估方法、装置及其设备 | |
CN108646199A (zh) | 一种动力电池soc估算方法和*** | |
CN113625174B (zh) | 一种锂离子电池soc与容量联合估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170929 |
|
RJ01 | Rejection of invention patent application after publication |