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

本发明专利技术公开了一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法。该方法首先依据CSTR连续系统模型，获得了状态分量包含待辨识参数的状态空间表达式；接着，借助欧拉算法对获取的非线性连续状态空间表达式进行了离散化处理，得到了相应的离散迭代模型。最后，运用无迹卡尔曼滤波算法进行多次迭代辨识，获得了准确的辨识结果。该算法收敛性好，且易于与已有软件相结合，具有很好的工程应用前景。

【技术实现步骤摘要】

本专利技术涉及一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，属于系统建模与参数辨识

技术介绍
连续搅拌釜式反应器(CSTR)是化工生产过程中典型的、高度非线性的化学反应系统。在化工生产的核心设备中占有相当重要的地位，在染料、医药试剂、食品及合成材料工业中，CSTR系统得到了广泛的应用。由于CSTR在化工生产中的重要作用，所以有必要对该过程进行详细的研究。在对CSTR建立模型进行分析时，模型的参数有时是未知的，所以对CSTR的模型参数辨识方法进行研究，具有重要的工程意义。然而，已有的方法如扩展卡尔曼滤波等，这些方法辨识结果有时会出现发散，得不到正确的结果。为了提高辨识效率和精度，研究运用无迹卡尔曼滤波算法进行CSTR模型参数辨识，具有重要的意义。
技术实现思路
为了有效的了解化工生产中CSTR化学反应系统，本专利技术提出了一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，有效的实现了CSTR模型的参数辨识。本专利技术的技术解决方案是：一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，其步骤如下：(1)、获取扩维状态向量中包含CSTR模型待辨识参数的状态空间表达式；(2)、运用欧拉算法对连续的状态空间表达式进行离散化，获得离散化的状态空间表达式；(3)、初始化，包括：设定参数辨识的初值初始参数辨识误差协方差以及过程噪声和量测噪声所满足的协方差矩阵Q和R，算法迭代次数最大值L；(4)、选取k-1时刻的sigma点，计算公式为： δ k - 1 = x ^ k - 1 x ^ k - 1 + γ · p ^ x , k - 1 x ^ k - 1 - γ · p ^ x , k - 1 ]]> γ = n + λ , λ = α 2 · ( n + k f ) - n ]]>式中，表示k-1时刻的状态估计值，表示k-1时刻的状态估计误差协方差，γ表示尺度参数，n表示的维数；常数α决定sigma点围绕均值的波动范围，通常α∈[10-4,1]；常数kf是另一尺度参数，用于状态估计和参数辨识时通常取0。(5)、在上一步基础上，计算k-1时刻的sigma点的增值点，计算公式为： ( f k - 1 * ) i = f ( ( δ k - 1 ) i , u k - 1 ) ]]> ( h k - 1 * ) i = h本文档来自技高网...

【技术保护点】
一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，其特征在于，包含如下步骤：(1)、获取状态分量中包含CSTR模型待辨识参数的状态空间表达式。(2)、运用欧拉算法对连续的状态空间表达式进行离散化，获得离散的状态空间表达式。(3)、初始化，包括：设定参数辨识的初值初始参数辨识误差协方差以及过程噪声和量测噪声所满足的协方差矩阵Q和R，算法迭代次数最大值L。(4)、选取k‑1时刻的sigma点，计算公式为：δk-1=x^k-1x^k-1+γ·p^x,k-1x^k-1-γ·p^x,k-1]]>γ=n+λ,λ=α2·(n+kf)-n]]>式中，表示k‑1时刻的状态估计值，表示k‑1时刻的状态估计误差协方差，γ表示尺度参数，n表示的维数；常数α决定sigma点围绕均值的波动范围，通常α∈[10‑4,1]；常数kf是另一尺度参数，用于状态估计和参数辨识时通常取0。(5)、在上一步基础上，计算k‑1时刻的sigma点的增值点，计算公式为：(fk-1*)i=f((δk-1)i,uk-1)]]>(hk-1*)i=h((δk-1)i,uk-1)]]>式中，f(·)是对应具体问题系统方程的非线性函数，h(·)是对应具体问题输出方程中的非线性函数，uk‑1是k‑1时刻输入控制矩阵，下标i表示对应于第i个sigma点的相关取值，i＝0…2n。(6)、计算k‑1时刻的状态向量均值和协方差，计算公式为：x^k-1-=Σi=02nwim(fk-1*)i]]>p^x,k-1-=Σi=02nwic((fk-1*)i-x^k-1-)((fk-1*)i-x^k-1-)T+Qk-1]]>式中，表示k‑1时刻的状态向量均值，表示k‑1时刻的状态协方差，Qk‑1表示k‑1时刻系统噪声所满足的协方差矩阵，权重系数和取值的计算公式如下：w0m=λn+λ,w0c=λn+λ+(1-α2+β)]]>wim=wic=12(n+λ),i=1...2n]]>式中β通常是包含x分布的先验知识，对高斯分布来说，其最优值一般取2。(7)、计算k‑1时刻量测向量均值和协方差，计算公式为：yk-1-=Σi=02nwim(hk-1*)i]]>p^y,k-1-=Σi=02nwic((hk-1*)i-yk-1-)((hk-1*)i-yk-1-)T+Rk-1]]>式中表示k‑1时刻量测向量均值，表示k‑1量测向量的协方差，Rk‑1表示k‑1时刻的量测噪声所满足的协方差矩阵。(8)、计算交互协方差，计算公式如下：p^xy,k-1-=Σi=02nwic((fk-1*)i-x^k-1-)((hk-1*)i-yk-1-)T]]>式中，表示k‑1时刻的交互协方差。(9)、在上一步的基础上，计算k‑1时刻的卡尔曼滤波增益，其遵循的计算公式为：Kk-1=p^xy,k-1-·(p^y,k-1-)-1]]>式中，Kk‑1表示k‑1时刻的卡尔曼滤波增益。(10)、运用无迹卡尔曼滤波更新步，获得k时刻的状态估计值和协方差，计算公式为：x^k=x^k-1-+Kk-1·(yk-1-yk-1-)]]>p^x,k=p^x,k-1--Kk-1·p^y,k-1-·Kk-1T]]>式中，表示k时刻的状态估计值，yk‑1表示k‑1时刻量测输出真实值，表示k时刻的估计协方差。(11)、依据上述步骤进行多次迭代辨识，直至k≥L时，结束迭代过程，输出辨识结果。...

【技术特征摘要】
1.一种基于无迹卡尔曼滤波算法的CSTR模型参数辨识方法，其特征在于，包含如下步骤：(1)、获取状态分量中包含CSTR模型待辨识参数的状态空间表达式。(2)、运用欧拉算法对连续的状态空间表达式进行离散化，获得离散的状态空间表达式。(3)、初始化，包括：设定参数辨识的初值初始参数辨识误差协方差以及过程噪声和量测噪声所满足的协方差矩阵Q和R，算法迭代次数最大值L。(4)、选取k-1时刻的sigma点，计算公式为： δ k - 1 = x ^ k - 1 x ^ k - 1 + γ · p ^ x , k - 1 x ^ k - 1 - γ · p ^ x , k - 1 ]]> γ = n + λ , λ = α 2 · ( n + k f ) - n ]]>式中，表示k-1时刻的状态估计值，表示k-1时刻的状态估计误差协方差，γ表示尺度参数，n表示的维数；常数α决定sigma点围绕均值的波动范围，通常α∈[10-4,1]；常数kf是另一尺度参数，用于状态估计和参数辨识时通常取0。(5)、在上一步基础上，计算k-1时刻的sigma点的增值点，计算公式为： ( f k - 1 * ) i = f ( ( δ k - 1 ) i , u k - 1 ) ]]> ( h k - 1 * ) i = h ( ( δ k - 1 ) i , u k - 1 ) ]]>式中，f(·)是对应具体问题系统方程的非线性函数，h(·)是对应具体问题输出方程中的非线性函数，uk-1是k-1时刻输入控制矩阵，下标i表示对应于第i个sigma点的相关取值，i＝0…2n。(6)、计算k-1时刻的状态向量均值和协方差，计算公式为： x ^ k - 1 - = Σ i = 0 2 n w i m ( f k - 1 * ) i ]]> p ^ x , k - 1 - = Σ i = 0 2 n w i c ( ( f k - 1 * ) i - x ^ k - 1 - ) ( ( f k - 1 * ) i - x ^ k - 1 - ) T + Q k - 1 ]]>式中，表示k-1时刻的状态向量均值，表示k-1时...

【专利技术属性】
技术研发人员：周琪，于占东，王焕清，王巍，
申请(专利权)人：渤海大学，
类型：发明
国别省市：辽宁;21