CN107703438A - 一种rc充放电回路下充放电特性曲线拟合算法 - Google Patents
一种rc充放电回路下充放电特性曲线拟合算法 Download PDFInfo
- Publication number
- CN107703438A CN107703438A CN201710849183.8A CN201710849183A CN107703438A CN 107703438 A CN107703438 A CN 107703438A CN 201710849183 A CN201710849183 A CN 201710849183A CN 107703438 A CN107703438 A CN 107703438A
- Authority
- CN
- China
- Prior art keywords
- characteristic curve
- value
- charging
- minizone
- actual
- 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
-
- 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/28—Testing of electronic circuits, e.g. by signal tracer
- G01R31/2832—Specific tests of electronic circuits not provided for elsewhere
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
本发明提供了一种RC充放电回路下充放电特性曲线拟合算法,包括:根据实际充、放电特性曲线,分别假定一条充、放电特性曲线,先通过采集样本数据分别确定充、放电特性曲线系数Vu2_a和E2_a可能分布的取值范围,将各取值范围分别定义为假定取值区间,将各假定取值区间分成N(N为≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为Vu2_a和E2_a当前值,通过算法自动调整系数B2_a或B`2_a,得到与小区间数目相同的若干个SUM值,然后将最小的SUM值对应的小区间作为Vu2_a或E2_a值所在的实际取值区间,在这个小区间分别对充放电特性系数Vu2_a或E2_a进行拟合,从而得到实际充、放电特性曲线特征值。其扩大了拟合技术的应用范围,同时达到高效拟合速度和准确度。
Description
技术领域
本发明涉及一种充放电回路研究,尤其涉及一种RC充放电回路下充放电特性曲线拟合算法。
背景技术
解析RC充放电电路的根本就是找到其对应的物理特性,这其中最重要的一步就是找出其对应的充放电方程。通过曲线拟合,可以找到满足RC充、放电特性曲线对应的物理系数,从而帮助分析。现有充、放电特性曲线高效高速拟合技术是在PC端上实现的,特别是基于MATLAB的拟合技术,都是基于PC强悍的硬件处理能力。但是嵌入式平台硬件处理能力是没法达到那种高度的,这就决定了在嵌入式平台上需要依靠算法优化来达到高速高效的拟合效果。
其次,目前电动汽车采用预估算法求解绝缘电阻值时,由于RC电路存在,导致芯片采集AD样本失准。由于RC电路存在导致电路元器件两端电压失准,电压延迟稳定现象,因而其计算结果经常有误,而且造成极大的工作量。
因而,急需一种拟合算法可应用于嵌入式平台,同时还能解决RC电路存在时采用预估算法对绝缘电阻值的精确求解。
发明内容
基于上述背景,本发明提供了一种RC充放电回路下充放电特性曲线拟合算法,其扩大了RC充放电回路曲线拟合技术的应用范围,同时达到高效充、放电特性曲线拟合速度和准确度。
为实现上述目的,本发明提供了一种RC充放电回路下充放电特性曲线拟合算法,包括:根据实际充电特性曲线1Vt=V0+(Vu1-V0)*[1-exp(-t*B1)]和实际放电特性曲线1Vt=E1*exp(-t*B`1)+Vc1,分别假定一条充电特性曲线2Vt=V0+(Vu2-V0)*[1-exp(-t*B2)]和一条放电特性曲线2Vt=E2*exp(-t*B`2)+Vc2,通过采集样本数据分别确定充电特性曲线系数Vu2_a和放电特性曲线系数E2_a可能分布的取值范围,将各取值范围分别定义为假定取值区间,将各假定取值区间分成N(N取≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为Vu2_a和E2_a当前值,通过算法自动调整系数B2_a或B`2_a,得到与小区间数目相同的若干个SUM值,然后将最小的SUM值对应的小区间作为Vu2_a或E2_a值所在的实际取值区间,在这个小区间分别对充电特性曲线系数Vu2_a或放电特性曲线系数E2_a进行拟合,从而得到实际充电特性曲线特征值(Vu2_b、B2_b)和实际放电特性曲线特征值(E2_b、B`2_b、Vc2_b),从而确定实际充电特性曲线1和实际放电特性曲线1(V0为充电开始时监测对象电压值;Vu1、Vu2为充电完成监测对象终止电压值;Vt为监测对象在任意时刻t电压值;B1、B2为RC充电特性系数;E1、E2为监测对象从开始放电到放电结束总的放电电压差值;B`1、B`2为RC放电特性系数;Vc1、Vc2为监测对象放电结束后本身存在的稳定电压值;exp()为指数函数)。
需要说明的是将各假定取值区间分成N(N取≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为Vu2_a和E2_a当前值,存在下述两种情形:若假定取值区间很小时,N取1,即仅存在一个长度相同的区间,此时针对这个取值范围取中间值即可;若假定取值区间需要再进一步细分,则N取≥2的整数,此时至少存在两个长度相同的小区间,针对至少2个小区间取中间值作为Vu2_a和E2_a当前值。
进一步的,所述通过采集样本数据确定充电特性曲线系数Vu2_a可能分布的取值范围,包括:假设实际充电特性曲线1Vt=V0+(Vu1-V0)*[1-exp(-t*B1)]有n个样本点C0~Cn-1,在充电特性曲线2Vt=V0+(Vu2-V0)*[1-exp(-t*B2)]上选取n个与样本点C0~Cn-1具有相同横坐标的点,通过后台运行采集样本点的数据确定充电特性曲线系数Vu2_a可能分布的取值范围。
进一步的,所述通过采集样本数据确定放电特性曲线系数E2_a可能分布的取值范围,包括:假设实际放电特性曲线1Vt=E1*exp(-t*B`1)+Vc1有n个样本点C0~Cn-1,在放电特性曲线2Vt=E2*exp(-t*B`2)+Vc2上选取n个与样本点C0~Cn-1具有相同横坐标的点,通过后台运行采集样本点的数据确定放电特性曲线系数E2_a可能分布的取值范围。
进一步的,所述将取值范围定义为假定取值区间,将假定取值区间分成N(N为≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为Vu2_a当前值,通过算法自动调整系数B2_a,得到与小区间数目相同的若干个SUM值,为通过最小二乘法,从而得到与区间数目相同的若干个SUM值。
SUM的计算方式为:
(i-表示第i个样本点,T[i]-表示实际充电特性曲线1第i个样本点的横坐标,buf[i]-表示第i个样本点在实际充电特性曲线1上的纵坐标)
进一步的,所述通过最小二乘法从而得到与区间数目相同的若干个SUM值,包括:将实际充电特性曲线1的横坐标皆设定为间隔相同的时刻值T[i](i取0,1,2,...n-1),纵坐标为相应时刻采集出的Vt值,将Vu2_a的可能取值范围定义为假定取值区间[Vu_Cn-1,x*Vu_Cn-1](x取2,...n),Vu_Cn-1为实际充电特征曲线上点Cn-1的Vt值,将假定取值区间再分成N(N取≥1的整数)个长度相同的小区间,分别取各长度相同的小区间的中间值Vum(m取1,2....n),通过最小二乘法的算法自动调整B2_a值计算不同组合(Vum,B2_a)下的最小SUM值,其中最小SUM值对应的小区间就是实际充电特性曲线1的Vu2_a的实际取值区间。
通过最小二乘法的算法自动调整B2值是指以最小二乘法求得的SUM值为判断标准力求SUM值最小,具体的指,当增大B2值时,SUM变大,则减小B2,直到SUM最小为止;当增大B2值时,SUM变小,则继续增大B2,直到SUM最小为止;同理,当减小B2值时,SUM值也变大,则增大B2值,直到SUM最小为止;当减小B2值时,SUM值也变小,则继续减小B2值,直到SUM最小为止。
进一步的,将取值范围定义为假定取值区间,将假定取值区间分成N(N取≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为E2_a当前值,通过算法自动调整系数B`2_a,得到与小区间数目相同的若干个SUM值,为通过最小二乘法,从而得到与区间数目相同的若干个SUM值。
SUM的计算方式为:
(i-表示第i个样本点,T[i]-表示实际放电特性曲线1第i个样本点的横坐标,buf[i]-表示实际放电特性曲线1第i个样本点的纵坐标)
进一步的,所述通过最小二乘法从而得到与区间数目相同的若干个SUM值,包括:将实际放电特性曲线1的横坐标皆设定为间隔相同的时刻值T[i](i取0,1,2,...n-1),纵坐标为相应时刻采集出的Vt值,刚开始放电时的Vt为Vf,因为E2_a=Vf-Vc2_a,将E2_a的可能取值范围定义为假定取值区间[0,Vf],将假定取值区间再分成N(N取≥1的整数)个长度相同的区间,分别取各长度相同的区间的中间值Em(m取1,2....n),通过最小二乘法的算法自动调整B`2_a值计算不同组合(Em,B`2_a,Vf-Em)下的最小SUM值,其中最小SUM值对应的小区间就是实际放电特性曲线1的E2_a的实际取值区间。
通过最小二乘法的算法自动调整B`2_a值是指以最小二乘法求得的SUM值为判断标准力求SUM值最小,具体的指,当增大B`2_a值时,SUM变大,则减小B`2_a,直到SUM最小为止;当增大B`2_a值时,SUM变小,则继续增大B`2_a,直到SUM最小为止;同理,当减小B`2_a值时,SUM值也变大,则增大B`2_a值,直到SUM最小为止;当减小B`2_a值时,SUM值也变小,则继续减小B`2_a值,直到SUM最小为止。
进一步的所述将最小的SUM值对应的小区间作为Vu2_a所在的实际取值区间,在这个小区间对充电特性曲线系数Vu2_a进行拟合,从而得到实际充电特性曲线特征值(Vu2_b、B2_b)为通过最小二乘法,计算实际取值区间内不同组合下的(Vu2_a、B2_a)值的SUM值,找到SUM值最小值,就是最终实际充电特性曲线特征值(Vu2_b、B2_b)。
进一步的,所述将最小的SUM值对应的小区间作为E2_a值所在的实际取值区间,在这个小区间对放电特性曲线系数E2_a进行拟合,从而得到实际放电特性曲线特征值(E2_b、B`2_b、Vc2_b),为通过最小二乘法计算实际取值区间内不同组合下的(E2_a、B`2_a、Vc2_a)值的SUM值,找到SUM值最小值,就是最终实际放电特性曲线特征值(E2_b、B`2_b、Vc2_b)。
本发明就是通过将拟合结果分布的可能区间划分为若干小区间,确定实际Vu或者E分布的小区间后再在这个小区间进行多次的拟合,求得最优解,精确确定最终拟合结果,其准确率高。这样RC充、放电特性曲线的快速拟合,就不再仅仅局限于PC端了,从而促进了RC曲线拟合技术在嵌入式平台上的运用,优化了嵌入式平台的数据处理方法,减小了数据处理量,让嵌入式平台拥有了对满足充放电特性数据的高速高效拟合能力。同样的,本发明新拟合算法充分考虑到由于RC电路存在造成的电路元器件两端电压失准问题采用预估思想,通过分析已有样本点,预估出最终RC电路稳定后电路中元器件两端的准确电压,对具有RC电路特性的现象的绝缘电阻值的预估算法精确求解对分析提供了一种新的方法。
附图说明
图1为充电特性曲线拟合图。
图2为充电特性曲线系数特征值拟合过程图。
图3为放电特性曲线拟合图。
图4为放电特性曲线系数特征值拟合过程图。
具体实施方式
下面通过具体实施方式来进一步说明本发明的技术方案,但不构成对本发明的任何限制。
首先说明充电特性曲线。
对于充电特性曲线方程,我们需要拟合的有两个系数,即充电结束后终止电压Vu和充电物理特性系数B。具体拟合方法如下(如图1所示):假设n个样本点C0~Cn-1所在的曲线Vt=V0+(Vu1-V0)*[1-exp(-t*B1)]为实际充电特性曲线1,即我们需要求解的充电方程系数代表的充电特性曲线。为了求解实际充电曲线系数,我们先给定一条充电特性曲线Vt=V0+(Vu2-V0)*[1-exp(-t*B2)],即假设充电特性曲线2,在假设充电特性曲线2上选取n个与样本点C0~Cn-1具有相同横坐标的点D0~Dn-1。首先通过采集样本数据确定充电特性参数Vu2_a可能分布的取值范围,将该范围定义为假定取值区间[Vu_Cn-1,x*Vu_Cn-1](x为正整数2,3,4,5....,依实际情况而定),将该假定取值区间均分成若干段长度相同的小区间,在该假定取值小区间取中间值作为Vu2_a当前值,然后通过最小二乘法,通过算法自动调整系数B2_a,得到与小区间数相同的若干个SUM值,然后将最小的SUM值对应的小区间[Vu_CX,Vu_CY](注:[Vu_CX,Vu_CY]属于[Vu_Cn-1,x*Vu_Cn-1]的一部分。)作为Vu2_a值所在的实际区间。
找到了充电曲线特征值Vu2_a所在的实际区间[Vu_CX,Vu_CY]后,利用最小二乘法,求解不同组合下(Vu2_a、B2_a)(其中B2_a根据算法自动求解)下所有最小SUM值组,最后选择出组中最小SUM对应的(Vu2_a、B2_a),即最终充电特性曲线参数(Vu2_b、B2_b)。Vu2_a拟合如下:
很明显,Cn-1点处于充电完成或者未完成阶段,如图1所示,对应的Vt值肯定小于或等于充电结束的Vu值。已知Cn-1点对应的Vt值,假设为Vu_Cn-1,这个Vt明显小于或等于最终曲线充电完成后的Vt值,即最终的Vu值。那么首先给定一个Vu的搜索区间[Vu_Cn-1,x*Vu_Cn-1](x为正整数2,3,4,5....,依实际情况而定),将这个区间分成若干相等长度区间,如将[10,20]分成5个长度相等的小区间[10,12)、[12,14)、[14,16)、[16,18)、[18,20]。然后取这若干个小区间的中间值,如[10,20]的5个区间的11、13、15、17、19,让充电曲线Vt=V0+(Vu2-V0)*[1-exp(-t*B2)]中Vu2分别取各个区间的中间值,通过算法自动调整B2值,每个Vu2值和最终的B2值,都对应一种充电曲线的情况,通过不断改变每种充电曲线2的系数B2,使得当前Vu2、B2对应的曲线2与曲线1最接近,同时得到了最小二乘法下的一个最小SUM值。有多少个长度相同的区间就有多少个最终的SUM值,取出这几个中最小的SUM值和对应的Vu2和B2值,就可以得到实际充电特性曲线特征值Vu所在的区间[Vu_CX,Vu_CY](很明显,[Vu_CX,Vu_CY]属于[Vu_Cn-1,x*Vu_Cn-1]的一部分)。
SUM的计算方式为:
(i-表示第i个样本点,T[i]-表示实际充电特性曲线1第i个样本点的横坐标,buf[i]-表示第i个样本点在实际充电特性曲线1上的纵坐标)
然后直接在上面确定的Vu所属区间[Vu_CX,Vu_CY]中,利用最小二乘法,通过不同组合(Vu2、B2)下的SUM值,找到SUM值最小时的(Vu2、B2)值,就是最终充电特性曲线系数(Vu2_b、B2_b)值。
这个过程中样本点横坐标设定为间隔相同的时刻值T[i](i取0,1,2,...n-1),比如10、30、50、70...,纵坐标为相应时刻采集出的Vt值。规定将横坐标相同的两个样本点纵坐标差值的平方和SUM为拟合判断标准,很明显,当SUM值小于一定值时候,说明假设充电特性曲线2和实际充电特性曲线1重合或者很接近,此时的(Vu2、B2)就是实际充电曲线系数(Vu2_b、B2_b)。
我们采用如下方法,实现以最快最高效的方式找到实际充电曲线系数。很明显,充电结束后,T[n-1]时刻采集到的Vt值,最接近充电结束后Vu值,而实际Vu值肯定大于或等于这个Vt值,那么我们就可以在大于这个Vt值的一定范围内搜索实际的充电结束后的Vu值。假设充电结束后Vu取值范围为(Vu0,Vu5),可以取Vu0等于T[n-1]时刻采集到的Vt值,Vu5可以取为T[n-1]时刻采集到的Vt值的整数倍。将这个连续区间均分为N段,这里为了叙述方便,均分为5段,即[Vu0,Vu1)、[Vu1,Vu2)、[Vu2,Vu3)、[Vu3,Vu4)、[Vu4,Vu5),规定所有涉及到SUM的值为了区分都以SUMn(n取0,1,2,...)表示依次存储进缓冲数组C[n]中。首先通过在上面5段区间中分别取区间的中间值作为Vu2_a,然后结合最小二乘法的思路,通过相应算法自动调整B2_a,分别得到5个Vu2_a对应的最小SUM值,然后将5个SUM值中最小SUM值对应的[Vux,Vux+1)(x取值整数0...4)作为实际充电特性参数Vu所在的区间。最后在Vu2_a以前面确定的[Vux,Vux+1)(x取值整数0...4)范围内,通过不同组合(Vu2_a、B2_a)下计算而得到SUM值小于一定值时,此时的Vu2_a、B2_a就是最终的充电特性曲线特征值(Vu2_b、B2_b),充电特性拟合思路如图2所示。
其次,以放电特性曲线为例说明。
对于放电特性曲线方程,需要注意的是多一个拟合参数Vc(Vc为监测对象放电结束后本身存在的稳定电压值),但是Vc与放电电压变化量E之间存在一定的关系,所以从某种程度上讲,还是只有两个拟合系数B`和E。同样假设n个样本点C0~Cn-1所在的曲线Vt=E1*exp(-t*B`1)+Vc1为实际放电特性曲线1,如图3,即我们需要求解的放电方程系数代表的放电特性曲线。为了求解实际放电特性曲线系数,我们先给定一条放电特性曲线Vt=E2*exp(-t*B`2)+Vc2,即假设放电特性曲线2,在假设放电特性曲线2上选取n个与样本点C0~Cn-1具有相同横坐标的点D0~Dn-1,首先通过采集样本数据确定放电特性参数E2_a可能分布的取值范围,将该范围定义为假定取值区间[0,Vf],将假定取值区间均分成若干段长度相同的小区间,在假定取值小区间取中间值作为E2_a当前值,然后通过最小二乘法,通过算法自动调整系数B`2_a,得到与区间数相同的若干个SUM值,然后将最小的SUM值对应的小区间[E_FX,E_FY](很明显,[E_FX,E_FY]属于[0,Vf]的一部分)作为E2_a值所在的实际区间。
找到了放电曲线特征值E2_a所在的实际区间[E_CX,E_CY]后,利用最小二乘法,求解不同组合下(E2_a、B`2_a、Vf-E2_a)(其中B`2_a根据算法自动求解)下所有最小SUM值组,最后选择出组中最小SUM对应的(E2_a、B`2_a、Vf-E2_a),即最终充电特性曲线参数(E2_b、B`2_b、Vc2_b)。E2_a拟合如下:
在假设放电曲线2上选取n个与样本点C0~Cn-1具有相同横坐标的点D0~Dn-1。首先假设刚开始放电0时刻时Vt值为Vf。很明显,Cn-1点处于放电完成或者未完成阶段,如图1所示,对应的Vf-Vt值肯定小于或等于放电结束时E的值。放电分区间拟合算法主要是对E的取值范围进行界定。那么首先给定一个E的搜索区间[0,Vf],将这个区间分成若干相等长度区间,如将[0,20]分成5个长度相等的连续空间[0,4)、[4,8)、[8,12)、[12,16)、[16,20]。然后取这若干个区间的中间值,如[0,20]的5个区间的2、6、10、14、18,让放电曲线Vt=E2*exp(-t*B2)+Vc2中E2分别取各个区间的中间值,则Vc2取值为18、14、10、6、2,通过算法自动调整B`2值,每个E2值和最终的B`2值,都对应一种放电曲线的情况,通过不断改变每种放电曲线2的系数B`2,使得当前E2、B`2对应的曲线2与实际曲线1最接近,同时得到了最小二乘法下的一个最小SUM值。有多少个长度相同的区间就有多少个最终的SUM值,取出这几个中最小的SUM值和对应的E2和B`2值,就可以得到实际放电特性曲线特征值E所在的区间[E_FX,E_FY](很明显,[E_FX,E_FY]属于[0,Vf]的一部分)。
SUM的计算方式为:
(i-表示第i个样本点,T[i]-表示实际放电曲线1第i个样本点的横坐标,buf[i]-表示实际放电曲线1第i个样本点的纵坐标)
然后直接在上面确定的区间[E_FX,E_FY]中,利用最小二乘法,通过不同组合下的(E2、B`2,Vc2)值的SUM值,找到SUM值最小时的(E2、B`2、Vc2)值,就是最终充电特性曲线系数(E2_b、B`2_b、Vc2_b)值。
这个过程中样本点横坐标设定为间隔相同的时刻值T[i](i取0,1,2,...n-1),比如10、30、50、70...,纵坐标为相应时刻采集出的Vt值。我们规定将横坐标相同的两个样本点纵坐标差值的平方和SUM为拟合判断标准,计算不同组合(E,B`,Vc)下的SUM值,很明显,当SUM值小于一定值时候,说明假设放电曲线和实际放电曲线重合,此时的(E,B`,Vc)就是实际放电曲线系数。
我们采用如下方法,实现以最快最高效的方式找到实际放电曲线系数。假设刚开始放电时刻的Vt值为Vf5,那么最终放电结束后E值一定处于[0,Vf5]之间,我们可以取T[0]时刻的Vt值为Vf5。将这个连续区间均分为N段,这里为了叙述方便,均分为5段,即[0,Vf1)、[Vf1,Vf2)、[Vf2,Vf3)、[Vf3,Vf4)、[Vf4,Vf5)。很明显,放电完全结束后E+Vc=Vf5;假设放电结束后Vt值处于[0,Vf1)区间内,那么此时E值范围[Vf5-Vf1,Vf5],此时Vc值范围为[0,Vf1];当Vt值处于[Vf1,Vf2)区间内,那么此时E值范围[Vf5-Vf2,Vf5-Vf1],此时Vc值范围为[Vf1,Vf2],以此类推。规定所有涉及到SUM的值为了区分都以SUMn(n取0,1,2,...)表示;实际放电曲线上样本点纵坐标依次存储进缓冲数组F[n](n为样本点个数)中。首先通过在上面5段区间中分别取区间的中间值作为E2_a,然后结合最小二乘法的思路,通过相应算法自动调整B`2_a,分别得到5个E2_a对应的最小SUM值然后将5个SUM值中最小SUM值对应的[Vfx,Vfx+1)(x取值整数0...4,Vf0=0)作为实际放电特性参数E所在的区间。最后在E2_a以前面确定的[Vfx,Vfx+1)(x取值整数0,...4,Vf0=0)范围内,通过不同组合(E2_a、B`2_a、Vc2_a)下计算而得到SUM值小于一定值时,E2_、B`2、Vc2就是最终放电特性曲线特征值(E2_b、B`2_b、Vc2_b)。放电特性拟合思路如图4所示。
上面以将样本区间分成5段来说明如何对满足RC充、放电特性曲线进行快速拟合的算法思路,以相同的思路,可以将样本区间分得更细,这样使拟合过程变得更加的高效。
Claims (9)
1.一种RC充放电回路下充放电特性曲线拟合算法,包括:根据实际充电特性曲线1Vt=V0+(Vu1-V0)*[1-exp(-t*B1)]和实际放电特性曲线1Vt=E1*exp(-t*B`1)+Vc1,分别假定一条充电特性曲线2Vt=V0+(Vu2-V0)*[1-exp(-t*B2)]和一条放电特性曲线2Vt=E2*exp(-t*B`2)+Vc2,通过采集样本数据分别确定充电特性曲线系数Vu2_a和放电特性曲线系数E2_a可能分布的取值范围,将各取值范围分别定义为假定取值区间,将各假定取值区间分成N(N取≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为Vu2_a和E2_a当前值,通过算法自动调整系数B2_a或B`2_a,得到与小区间数目相同的若干个SUM值,然后将最小的SUM值对应的小区间作为Vu2_a或E2_a值所在的实际取值区间,在这个小区间分别对充电特性曲线系数Vu2_a或放电特性曲线系数E2_a进行拟合,从而得到实际充电特性曲线特征值(Vu2_b、B2_b)和实际放电特性曲线特征值(E2_b、B`2_b、Vc2_b),从而确定实际充电特性曲线1和实际放电特性曲线1(V0为充电开始时监测对象电压值;Vu1、Vu2为充电完成监测对象终止电压值;Vt为监测对象在任意时刻t电压值;B1、B2为RC充电特性系数;E1、E2为监测对象从开始放电到放电结束总的放电电压差值;B`1、B`2为RC放电特性系数;Vc1、Vc2为监测对象放电结束后本身存在的稳定电压值;exp()为指数函数)。
2.根据权利要求1所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,所述通过采集样本数据确定充电特性曲线系数Vu2_a可能分布的取值范围,包括:假设实际充电特性曲线1Vt=V0+(Vu1-V0)*[1-exp(-t*B1)]有n个样本点C0~Cn-1,在充电特性曲线2Vt=V0+(Vu2-V0)*[1-exp(-t*B2)]上选取n个与样本点C0~Cn-1具有相同横坐标的点,通过后台运行采集样本点的数据确定充电特性曲线系数Vu2_a可能分布的取值范围。
3.根据权利要求1所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,所述通过采集样本数据确定放电特性曲线系数E2_a可能分布的取值范围,包括:假设实际放电特性曲线1Vt=E1*exp(-t*B`1)+Vc1有n个样本点C0~Cn-1,在放电特性曲线2Vt=E2*exp(-t*B`2)+Vc2上选取n个与样本点C0~Cn-1具有相同横坐标的点,通过后台运行采集样本点的数据确定放电特性曲线系数E2_a可能分布的取值范围。
4.根据权利要求1所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,所述将取值范围定义为假定取值区间,将假定取值区间分成N(N为≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为Vu2_a当前值,通过算法自动调整系数B2_a,得到与小区间数目相同的若干个SUM值,为通过最小二乘法,从而得到与区间数目相同的若干个SUM值。
5.根据权利要求4所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,所述通过最小二乘法从而得到与区间数目相同的若干个SUM值,包括:将实际充电特性曲线1的横坐标皆设定为间隔相同的时刻值T[i](i取0,1,2,...n-1),纵坐标为相应时刻采集出的Vt值,将Vu2_a的可能取值范围定义为假定取值区间[Vu_Cn-1,x*Vu_Cn-1](x取2,...n),Vu_Cn-1为点Cn-1的Vt值,将假定取值区间再分成N(N取≥1的整数)个长度相同的小区间,分别取各长度相同的小区间的中间值Vum(m取1,2....n),通过最小二乘法的算法自动调整B2_a值计算不同组合(Vum,B2_a)下的最小SUM值,其中最小SUM值对应的小区间就是实际充电特性曲线1的Vu2_a的实际取值区间。
6.根据权利要求1所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,将取值范围定义为假定取值区间,将假定取值区间分成N(N取≥1的整数)个长度相同的小区间,在各假定取值小区间取中间值作为E2_a当前值,通过算法自动调整系数B`2_a,得到与小区间数目相同的若干个SUM值,为通过最小二乘法,从而得到与区间数目相同的若干个SUM值。
7.根据权利要求6所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,所述通过最小二乘法从而得到与区间数目相同的若干个SUM值,包括:将实际放电特性曲线1的横坐标皆设定为间隔相同的时刻值T[i](i取0,1,2,...n-1),纵坐标为相应时刻采集出的Vt值,刚开始放电时的Vt为Vf,因为E2_a=Vf-Vc2_a,将E2_a的可能取值范围定义为假定取值区间[0,Vf],将假定取值区间再分成N(N取≥1的整数)个长度相同的区间,分别取各长度相同的区间的中间值Em(m取1,2....n),通过最小二乘法的算法自动调整B`2_a值计算不同组合(Em,B`2_a,Vf-Em)下的最小SUM值,其中最小SUM值对应的小区间就是实际放电特性曲线1的E2_a的实际取值区间。
8.根据权利要求1所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,所述将最小的SUM值对应的小区间作为Vu2_a所在的实际取值区间,在这个小区间对充电特性曲线系数Vu2_a进行拟合,从而得到实际充电特性曲线特征值(Vu2_b、B2_b)为通过最小二乘法,计算实际取值区间内不同组合下的(Vu2_a、B2_a)值的SUM值,找到SUM值最小值,就是最终实际充电特性曲线特征值(Vu2_b、B2_b)。
9.根据权利要求1所述的RC充放电回路下充放电特性曲线拟合算法,其特征在于,所述将最小的SUM值对应的小区间作为E2_a值所在的实际取值区间,在这个小区间对放电特性曲线系数E2_a进行拟合,从而得到实际放电特性曲线特征值(E2_b、B`2_b、Vc2_b),为通过最小二乘法计算实际取值区间内不同组合下的(E2_a、B`2_a、Vc2_a))值的SUM值,找到SUM值最小值,就是最终实际放电特性曲线特征值(E2_b、B`2_b、Vc2_b)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710849183.8A CN107703438A (zh) | 2017-09-18 | 2017-09-18 | 一种rc充放电回路下充放电特性曲线拟合算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710849183.8A CN107703438A (zh) | 2017-09-18 | 2017-09-18 | 一种rc充放电回路下充放电特性曲线拟合算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107703438A true CN107703438A (zh) | 2018-02-16 |
Family
ID=61172982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710849183.8A Pending CN107703438A (zh) | 2017-09-18 | 2017-09-18 | 一种rc充放电回路下充放电特性曲线拟合算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107703438A (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0980130A (ja) * | 1995-09-18 | 1997-03-28 | Nippon Telegr & Teleph Corp <Ntt> | 蓄電池容量測定方法 |
CN1431521A (zh) * | 2003-01-16 | 2003-07-23 | 华南理工大学 | 锂离子电池电量的测量方法及其装置 |
CN102624017A (zh) * | 2012-03-22 | 2012-08-01 | 清华大学 | 一种基于负荷预测的电池储能系统削峰填谷实时控制方法 |
CN103592579A (zh) * | 2013-07-30 | 2014-02-19 | 中国电力科学研究院 | 一种冲击放电电压计算模型构建方法 |
DE102013113951A1 (de) * | 2013-01-31 | 2014-07-31 | GM Global Technology Operations LLC | Verfahren zum Detektieren von Leerlaufsspannungsverschiebungen mittels Optimierung durch Anpassen der Anodenelektrodenhalbzellspannungskurve |
WO2014185163A1 (ja) * | 2013-05-16 | 2014-11-20 | 日本電気株式会社 | 電池状態推定装置、電池状態管理システム、電池、電池状態推定方法、及び、プログラム |
CN104657520A (zh) * | 2013-11-18 | 2015-05-27 | 国家电网公司 | 一种基于大容量储能锂离子电池的电池建模方法 |
CN104749529A (zh) * | 2015-04-01 | 2015-07-01 | 上海理工大学 | 锂电池充放电特性的标定方法及充放电特性标定仪 |
CN105677979A (zh) * | 2016-01-07 | 2016-06-15 | 北京北交新能科技有限公司 | 基于Butler-Volmer方程的钛酸锂电池改进模型 |
US20170212170A1 (en) * | 2016-01-21 | 2017-07-27 | Yokogawa Electric Corporation | Secondary battery capacity measuring system and secondary battery capacity measuring method |
CN106997028A (zh) * | 2017-05-19 | 2017-08-01 | 江苏大学 | 一种电池差分电压及容量增量曲线高效求解方法 |
-
2017
- 2017-09-18 CN CN201710849183.8A patent/CN107703438A/zh active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0980130A (ja) * | 1995-09-18 | 1997-03-28 | Nippon Telegr & Teleph Corp <Ntt> | 蓄電池容量測定方法 |
CN1431521A (zh) * | 2003-01-16 | 2003-07-23 | 华南理工大学 | 锂离子电池电量的测量方法及其装置 |
CN102624017A (zh) * | 2012-03-22 | 2012-08-01 | 清华大学 | 一种基于负荷预测的电池储能系统削峰填谷实时控制方法 |
DE102013113951A1 (de) * | 2013-01-31 | 2014-07-31 | GM Global Technology Operations LLC | Verfahren zum Detektieren von Leerlaufsspannungsverschiebungen mittels Optimierung durch Anpassen der Anodenelektrodenhalbzellspannungskurve |
WO2014185163A1 (ja) * | 2013-05-16 | 2014-11-20 | 日本電気株式会社 | 電池状態推定装置、電池状態管理システム、電池、電池状態推定方法、及び、プログラム |
CN103592579A (zh) * | 2013-07-30 | 2014-02-19 | 中国电力科学研究院 | 一种冲击放电电压计算模型构建方法 |
CN104657520A (zh) * | 2013-11-18 | 2015-05-27 | 国家电网公司 | 一种基于大容量储能锂离子电池的电池建模方法 |
CN104749529A (zh) * | 2015-04-01 | 2015-07-01 | 上海理工大学 | 锂电池充放电特性的标定方法及充放电特性标定仪 |
CN105677979A (zh) * | 2016-01-07 | 2016-06-15 | 北京北交新能科技有限公司 | 基于Butler-Volmer方程的钛酸锂电池改进模型 |
US20170212170A1 (en) * | 2016-01-21 | 2017-07-27 | Yokogawa Electric Corporation | Secondary battery capacity measuring system and secondary battery capacity measuring method |
CN106997028A (zh) * | 2017-05-19 | 2017-08-01 | 江苏大学 | 一种电池差分电压及容量增量曲线高效求解方法 |
Non-Patent Citations (1)
Title |
---|
张彩萍 等: "用基于遗传优化的扩展卡尔曼滤波", 《吉林大学学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103326598A (zh) | 三相整流器快速模型预测控制方法 | |
CN111428355B (zh) | 一种电力负荷数字统计智能综合的建模方法 | |
CN103259279B (zh) | 基于风电功率min级分量波动的风电场储能容量优化方法 | |
CN102545665A (zh) | 一种三电平pwm整流器直接功率控制方法 | |
CN111222097B (zh) | 光伏直流升压汇集系统可靠性评估及灵敏度分析方法 | |
CN104657520A (zh) | 一种基于大容量储能锂离子电池的电池建模方法 | |
CN110912115B (zh) | 用于混合储能电路平抑功率波动的预测控制方法及系统 | |
CN103970022B (zh) | 一种基于gm(1,n)灰色预测模型的均衡控制方法 | |
Liao et al. | Research on two-stage equalization strategy based on fuzzy logic control for lithium-ion battery packs | |
CN103997241B (zh) | 一种基于h桥模块的多电平变换器电容电压均压方法 | |
Caspar et al. | Optimal cell balancing with model-based cascade control by duty cycle adaption | |
CN113206508A (zh) | 基于预测的微电网储能电池主动均衡方法 | |
CN109889052A (zh) | 一种模块化多电平矩阵换流器电容电压的控制方法及装置 | |
CN110783913A (zh) | 基于群的计及预想事故集最佳电网拓扑在线优化方法 | |
CN107703438A (zh) | 一种rc充放电回路下充放电特性曲线拟合算法 | |
CN110297452A (zh) | 一种蓄电池组相邻型均衡系统及其预测控制方法 | |
CN102184296B (zh) | 基于实测数据的电气化铁路冲击负荷建模方法 | |
CN105958850A (zh) | 基于优化归并排序的模块化多电平换流器电容均压方法 | |
CN110323779B (zh) | 一种分布式发电与储能装置的功率动态聚合的方法和系统 | |
CN110766240A (zh) | 一种不同场景下的快速充电站分层储能配置方法 | |
CN112217226A (zh) | 一种适用于双向dc-dc变换器的改进无模型预测控制方法 | |
Wang et al. | A distributed state-of-charge balancing control scheme for three-port output-series converters in DC hybrid energy storage systems | |
CN113258563B (zh) | 柔直孤岛接入的交流电网继电保护灵敏度确定方法及装置 | |
CN111092425A (zh) | 一种基于拓扑特性的网架聚类分析方法及系统 | |
Wang et al. | An SMPS-based lithium-ion battery test system for internal resistance measurement |
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 | ||
AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20200501 |
|
AD01 | Patent right deemed abandoned |