(完整word版)等效电路模型参数在线辨识.docVIP

第四章 等效电路模型参数在线辨识 通过第三章函数拟合的方法可以确定钒电池等效电路模型中的参数， 但是在实际运行过程中模型参数随着工作环境温度、 充放电循环次数、 SOC等因素发生 变化，根据离线试验数据计算得到的参数值估算电池 SOC 可能会造成较大的估计误差。因此，在实际运行时，应对钒电池等效电路模型参数进行在线辨识，做出实时修正，提高基于模型估算 SOC的精度。 4.1 基于遗忘因子的最小二乘算法 参数辨识是根据被测系统的输入输出来， 通过一定的算法， 获得让模型输出值尽量接近系统实际输出值的模型参数估计值。 根据能否实时辨识系统的模型参数，可以将常用的参数辨识方法分为离线和在线两类， 离线辨识只能在数据采集完成后进行，不能对系统模型实时地在线调整参数， 对于具有非线性特性的电池系统往往不能得到满意的辨识结果； 在线辨识方法一般能够根据实时采集到的数据对系统模型进行辨识， 在线调整系统模型参数。 常用的辨识方法有最小二乘法、极大似然估计法和 Kalman 滤波法等。因最小二乘法原理简明、收敛较快、容易理解和掌握、方便编程实现等特点， 在进行电池模型参数辨识时采用了效果较好的含遗忘因子的递推最小二乘法。 4.1.1 批处理最小二乘法简介 假设被辨识的系统模型： y( z)b1z 1 b2 z 2 L bn z n n （ 4-1） G( z) 1 a2 z 2 L an z u( z) 1 a1 z 其相应的差分方程为： n n y(k) ai y( k i ) bi u(k i ) （ 4-2） i 1 i 1 若考虑被辨识系统或观测信息中含有噪声，则被辨识模型式（ 4-2）可改写为： n n z(k) ai y(k i ) bi u( k i ) v(k) （ 4-3） i 1 i 1 式中，z(k) 为系统输出量的第 k 次观测值；y(k ) 为系统输出量的第 k 次真值， y(k i ) 为系统输出量的第 k i 次真值； u(k ) 为系统的第 k 个输入值， u(k i ) 为 系统的第 k i 个输入值； v( k) 为均值为 0 的随机噪声。 令： h( k) [ y( k 1), y(k 2),L , y(k n), u( k 1),u(k 2),L , u(k 3)] （ 4-4） [ a1, a2 ,L an , b1, b2 ,L bn ]T 则式（ 4-3）可变换为： z(k) h( k)θ v(k ) （ 4-5） 式中， θ为待估参数。 k 1,2,L m ，则有 z(1) h(1) y(0) L y(1 n) u(0) L u(1 n) Z m z(2) ,H m h(2) y(1) L y(2 n) u(1) L u(2 n) M M M M M M M M z(m) h(m) y(m 1) L y( m n) u(m 1) L u(m n) θ [a1, a2 ,L an ,b1,b2 ,L bn ]T , Vm [ v(1),v(2), L , v(m)] T 于是，式（ 4-5）的矩阵形式为 Z m H mθ Vm （ 4-6） 最 小 二 乘 法的 思 想 就 是寻 找一 个 θ的估 计 值 θ?， 使 得 各 次 测 量 的 Zi (i 1,2,L ? ? ? , m) 与由估计 θ确定的测量估计 Zi H i θ之差的平方和最小，即 ? (Z m ? T ? min （ 4-7） J ( ) H mθ) (Z m H mθ) 使 J ( ) min 的 估计值记作 ?，称作参数 的最小二乘估计，其值为 ? T 1 T θ (H m H m ) H m Z m （ 4-8） 最小二乘估计虽然不能满足式 （ 4-6）中的每个方程， 使每个方程都有偏差，但是它使所有方程偏差的平方和达到最小， 兼顾了所有方程的近似程度， 使整体误差达到最小，有利于抑制测量误差。 4.1.2 递推最小二乘法简介 前面给出的批处理最小二乘法是拟合工具箱进行数据拟合遵循的基本原理， 可以进行简单的离线辨识。 若每次处理的数据量较大， 应用批处理最小二乘法时， 不仅占用内存大，而且不能用于参数在线实时估计。 而电池系统是强非线性系统， 其参数受工作状态影响较大， 需要利用输入输出数据在线估计模型参数。 为了减 少计算量，减少数据在计算机中所占的内存， 更为了能够实时地辨识出电池系统 的特性，在用最小二乘法进行参数估计时，把它转化成参数递推的估计。 所谓参数递推估计算法就是当辨识系统在运行时， 每获得一次新的观测数据后，就是在上次估计结果的基础上，利用新引入的观测数据对上次估计的结果， 根据递推参数进行修正， 从而递推地得到新的参数估计值。 因此，递推