[发明专利]一种化学放热反应的分布式主元分析神经网络建模方法

摘要

本发明公开了一种化学放热反应的分布式主元分析神经网络建模方法。本发明通过采集对象的输入和输出数据，利用主元分析法通过降低维数将分布式参数系统分为自回归线性模型和非线性模型，通过最小二乘法对自回归线性模型进行辨识。非线性模型利用最小二乘法建立对象的RBF神经网络模型，然后通过遗传算法优化RBF‑神经网络模型的参数。本发明建立的模型具有较高的精确性，能够很好的描述过程对象的动态性能。

主权项

1.一种化学放热反应的分布式主元分析神经网络建模方法，其特征在于该方法的具体步骤是：步骤1、采集催化棒过程的实时运行数据，建立催化棒的分布式参数模型1.1以为催化棒输入的时空数据，为催化棒采集的输出数据，和相应的催化棒状态变量其中t是时间序列，L为时间序列的长度，z_i为采集的第i组催化棒输出数据所处的空间位置，N为采集的输出数据的总数；1.2催化棒时空变量通过傅里叶变换得：转换为有限空间得：其中是n次的逼近，是傅里叶变换得到的正交基函数，是时间系数；1.3由于基函数是正交向量，得到：其中(.，.)是向量的内积运算，则步骤1.2中的时空变量进一步变换得：1.4通过求解函数的最小值得形式如下；其中，定义||f(z)||＝(f(z)，f(z))^1/2，总体平均值其中f(z)为函数；求解的满足如下条件；其中是空间两点的相关函数，正交基函数线性表示为：1.5由步骤1.4可得：时间两点的相关函数定义为经过化简表示为：Cγ_it＝λ_iγ_it其中γit＝[γi1，...，γiL]是第i个特征向量，C是时间两点相关函数的矩阵；把特征值按照大小的顺序进行排列，即λ1＞λ2＞…＞λL，通过求取前面n个特征值使其满足前面n项的Ei之和大于98％，由此来确定n的大小；1.6把希尔伯特空间划分为H_s和H_f两个子空间其中系统表示为：其中As和Bs是需要辨识的矩阵，fs(Xs，0)是u≡0的非线性部分，Y为系统的输出变量，Cs为单位矩阵；步骤2.催化棒的线性时不变自回归模型的解耦辨识的具体步骤；2.1通过步骤1.6，空间的低维时间序列的输出yi(t)表示为：2.2忽略fs对催化棒系统造成的影响，将催化棒系统近似为线性时不变自回归模型；其中A(q‑1)，B(q‑1)分别为n×n，n×p项对角矩阵多项式，ny、nu分别为输入输出的阶次，q‑1是后移算子，进而变换得：其中y(t‑1)，…y(t‑n_y)分别表示t‑1时刻…t‑n_y时刻的输出；2.3通过使用最小二乘法进行辨识系统中θ；其中0＜μ＜1是遗忘因子，K(t)矩阵的权重系数，P(t)一个正定协方差矩阵，计算出Ai，Bj；i＝1，...，ny，j＝1，...，nu步骤3.基于RBF神经网络辨识催化棒非线性部分具体步骤如下；3.1通过一系列催化棒输出y(t)，并且在u(t)≡0的情况下，经过ΔTs采样间隔时间后，进行采样获得催化棒的输出y(t+ΔTs)，得到非线性部分为：fs(t)＝y(t+ΔTs)‑Ay(t)，fs＝[f1，f2，...，fn]3.2 RBF神经网络由输入层、隐含层、输出层3部分组成；其中包含ny个输入，nr个隐含层的节点，一个输出节点f1；隐含层的基函数选取为薄板样条函数，则第i层隐含节点可以表示为：hi(X)＝||X‑Ci||2lg(||X‑Ci||)，i＝1，2，...，nr其中Ci为隐含层的中心3.3 RBF神经网络的第j层的输出，表示为：其中Wij是隐含层到输出层的权重系数，通过最小二乘法获得：3.4由步骤3.2和步骤3.3进一步得到RBF神经网络的形式如下：步骤4.利用遗传算法优化RBF神经网络；4.1设定种群染色体数量为M，进化的最大迭代次数为N1，随机初始化种群；通过适应度函数优化径向基的个数和径向基中心数目可得目标函数：f＝es+ω·nr其中ω为权重系数，es是时间[t1 tL]内实际的输出与RBF神经网络的输出误差总和：4.2采用十进制编码方式对染色体进行编码，第i个的染色体可以表示为：其中i＝1，2，…，M，M是种群染色体数量，nr是隐含节点的个数，ny是输入节点的个数，1≤nr≤D，D是隐含层节点的最大个数，因此Ci′为m×D的矩阵；染色体Ci′中的元素为：cij＝fmin+r(fmax‑fmin)，1≤i≤nr，1≤j≤ny其中r是位于[0.1，1]之间的随机变量，fmax，fmin分别为RBF神经网络的输出的最大值和最小值；4.3选取染色体的操作算子的具体步骤为：4.3.1染色体交叉运算；选取交叉算子Pc，使染色体Ci′和下一个染色体Ci+1′以概率Pc进行交叉运算，产生下一代染色体；4.3.2染色体校正运算；为了产生新的径向基神经网络结构，以校正算子Pr的概率增加或者减少隐含层节点的个数，从而产生新的隐含层点个数nr＝nr+r1，r1∈[‑D/4，D/4]随机的整数，进而改变染色体Ci′中元素cij；4.3.3染色体变异运算；为了得到搜索空间中更好的解，以Pm的概率改变隐含层节点的个数nr，进而改变染色体中的元素Ci′，从而完成染色个体变异的操作；4.4依照步骤4.2到4.3中的步骤进行循环重复优化搜索，如果达到最大的进化次数N1结束优化搜索计算，得到经过遗传算法优化后的染色体Ci′，进一步通过解码得到RBF神经网络的参数。