[发明专利]一种基于randls模型的电池参数在线辨识方法有效
申请号: | 201710504706.5 | 申请日: | 2017-06-28 |
公开(公告)号: | CN107462836B | 公开(公告)日: | 2020-04-17 |
发明(设计)人: | 王立业;王丽芳;廖承林;刘伟龙;陶成轩;张志刚 | 申请(专利权)人: | 中国科学院电工研究所 |
主分类号: | G01R31/367 | 分类号: | G01R31/367;G01R31/387;G01R31/392 |
代理公司: | 北京科迪生专利代理有限责任公司 11251 | 代理人: | 关玲 |
地址: | 100190 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 randls 模型 电池 参数 在线 辨识 方法 | ||
1.一种基于randls模型的电池参数在线辨识方法,其特征在于:所述的辨识方法根据randls电池模型拓扑结构,建立输入输出传递函数,然后对randls电池模型离散化;最后采用递推增广最小二乘法(RELS)对模型参数进行迭代辨识,所述根据randls电池模型拓扑结构建立输入输出传递函数的方法如下:
根据randls电池模型,极化内阻Rn与容量电容Cn串联连接后跨接在m、n两端,极化电容Cb跨接在m、n两端,欧姆内阻Ro一端接m端,一端接输出uout;randls模型的输入为电流i,randls模型输出为randls电池模型两端的电压uout,randls电池模型s域形式的输入输出传递函数H(s)为:
公式(1)中,U(s)为电池端电压s域形式,I(s)为输入电流s域形式,a,b,c,d的含义如公式(2)所示:
对randls电池模型传递函数H(s)进行离散化的方法如下:
采用脉冲响应不变变换对式(1)所述的randls电池模型传递函数H(s)进行离散化,得到如下z域形式的传递函数:
根据z域和时域变换原理,公式(3)进而可以写成如下离散迭代形式:
其中,T为采样周期,将采样周期设为10秒;uk+1为k+1时刻randls电池模型两端的电压,uk为k时刻电池模型两端的电压,uk-1为k-1时刻randls电池模型两端的电压,ik为k 时刻randls电池模型输入电流,ik-1为k-1时刻randls电池模型输入电流,ik-2为k-2时刻randls电池模型输入电流;
randls电池模型的迭代参数矩阵表示为θ=[A1,A2,B0,B1,B2,C0,C1,C2]T,其中A1、A2、B0、B1、B2为公式(4)中的参数,C0、C1、C2为系统噪声,进而实现对randls电池模型物理参数的估算;randls电池模型物理参数与迭代参数之间的关系如下:
采用递推增广最小二乘法对randls电池模型物理参数进行迭代辨识的过程如下:
步骤1:对迭代过程中所需要的未知迭代参数和Pk进行初始化;
其中,Pk为迭代过程的中间矩阵,为randls电池模型迭代参数矩阵θ的估计值,分别为randls电池模型的迭代参数A1、A2、B0、B1、B2、C0、C1、C2的估计值,
将的初始值设为零矩阵,将Pk的初始值设为106I,其中I为单位矩阵;
步骤2:根据电池组k时刻的电压信号uk、k-1时刻的电压uk-1、k时刻的电流信号ik、k-1时刻的电流信号ik-1、k-2时刻的电流信号ik-2,以及k时刻的环境有色噪声估计值k-1时刻的环境有色噪声估计值k-2时刻的环境有色噪声估计值建立k+1时刻的系统观测值矩阵
步骤3:计算k+1时刻的迭代增益Kk+1:
其中为k+1时刻的系统观测值矩阵的转置矩阵;
步骤4:计算k+1时刻randls电池模型两端的电压预测值
其中为k时刻的系统观测值矩阵的转置矩阵;
步骤5:计算k+1时刻的迭代过程的中间矩阵Pk+1:
步骤6:计算k+1时刻的迭代参数矩阵
其中为k+1时刻randls电池模型两端电压的测量值;
步骤7:根据式公式(5)和辨识得到的k+1时刻迭代参数矩阵计算randls电池模型物理参数Ro、Rn、Cb、Cn;
步骤8:重复步骤:2至步骤7,直至程序止工作。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院电工研究所,未经中国科学院电工研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201710504706.5/1.html,转载请声明来源钻瓜专利网。