一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN105975747A(43)申请公布日2016.09.28(21)申请号201610272832.8(22)申请日2016.04.27(71)申请人渤海大学地址121000辽宁省锦州市松山新区科技路19号(72)发明人周琪于占东王焕清王巍(74)专利代理机构锦州辽西专利事务所21225代理人李辉(51)Int.Cl.G06F19/00(2011.01)权利要求书2页说明书5页附图2页(54)发明名称一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法(57)摘要本发明公开了一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法。该方法首先依据CSTR连续系统模型，获得了状态分量包含待辨识参数的状态空间表达式；接着，借助欧拉算法对获取的非线性连续状态空间表达式进行了离散化处理，得到了相应的离散迭代模型。最后，运用无迹卡尔曼滤波算法进行多次迭代辨识，获得了准确的辨识结果。该算法收敛性好，且易于与已有软件相结合，具有很好的工程应用前景。CN105975747ACN105975747A权利要求书1/2页1.一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，其特征在于，包含如下步骤：(1)、获取状态分量中包含CSTR模型待辨识参数的状态空间表达式。(2)、运用欧拉算法对连续的状态空间表达式进行离散化，获得离散的状态空间表达式。(3)、初始化，包括：设定参数辨识的初值初始参数辨识误差协方差以及过程噪声和量测噪声所满足的协方差矩阵Q和R，算法迭代次数最大值L。(4)、选取k-1时刻的sigma点，计算公式为：式中，表示k-1时刻的状态估计值，表示k-1时刻的状态估计误差协方差，γ表示尺度参数，n表示的维数；常数α决定sigma点围绕均值的波动范围，通常α∈[10-4,1]；常数kf是另一尺度参数，用于状态估计和参数辨识时通常取0。(5)、在上一步基础上，计算k-1时刻的sigma点的增值点，计算公式为：式中，f(·)是对应具体问题系统方程的非线性函数，h(·)是对应具体问题输出方程中的非线性函数，uk-1是k-1时刻输入控制矩阵，下标i表示对应于第i个sigma点的相关取值，i＝0…2n。(6)、计算k-1时刻的状态向量均值和协方差，计算公式为：式中，表示k-1时刻的状态向量均值，表示k-1时刻的状态协方差，Qk-1表示k-1时刻系统噪声所满足的协方差矩阵，权重系数和取值的计算公式如下：式中β通常是包含x分布的先验知识，对高斯分布来说，其最优值一般取2。(7)、计算k-1时刻量测向量均值和协方差，计算公式为：式中表示k-1时刻量测向量均值，表示k-1量测向量的协方差，Rk-1表示k-1时刻2CN105975747A权利要求书2/2页的量测噪声所满足的协方差矩阵。(8)、计算交互协方差，计算公式如下：式中，表示k-1时刻的交互协方差。(9)、在上一步的基础上，计算k-1时刻的卡尔曼滤波增益，其遵循的计算公式为：式中，Kk-1表示k-1时刻的卡尔曼滤波增益。(10)、运用无迹卡尔曼滤波更新步，获得k时刻的状态估计值和协方差，计算公式为：式中，表示k时刻的状态估计值，yk-1表示k-1时刻量测输出真实值，表示k时刻的估计协方差。(11)、依据上述步骤进行多次迭代辨识，直至k≥L时，结束迭代过程，输出辨识结果。3CN105975747A说明书1/5页一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法技术领域[0001]本发明涉及一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，属于系统建模与参数辨识技术领域。背景技术[0002]连续搅拌釜式反应器(CSTR)是化工生产过程中典型的、高度非线性的化学反应系统。在化工生产的核心设备中占有相当重要的地位，在染料、医药试剂、食品及合成材料工业中，CSTR系统得到了广泛的应用。[0003]由于CSTR在化工生产中的重要作用，所以有必要对该过程进行详细的研究。在对CSTR建立模型进行分析时，模型的参数有时是未知的，所以对CSTR的模型参数辨识方法进行研究，具有重要的工程意义。然而，已有的方法如扩展卡尔曼滤波等，这些方法辨识结果有时会出现发散，得不到正确的结果。为了提高辨识效率和精度，研究运用无迹卡尔曼滤波算法进行CSTR模型参数辨识，具有重要的意义。发明内容[0004]为了有效的了解化工生产中CSTR化学反应系统，本发明提出了一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，有效的实现了CSTR模型的参数辨识。[0005]本发明的技术解决方案是：一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，其步骤如下：[0006](1)、获取扩维状态向量中包含CSTR模型待辨识参数的状态空间表达式；[0007](2)、运用欧拉算法