智能支持向量机的农药废液焚烧炉炉温最佳化系统及方法

摘要

本发明公开了一种智能支持向量机的农药废液焚烧炉炉温最佳化系统及方法。训练样本经过加权最小二乘支持向量机的处理，得到的多个支持向量机输出结果与模糊系统结合得到炉温预报值和使炉温最佳的操作变量值。在本发明中，标准化处理模块与模糊系统模块相连，用于对训练样本中心化；模糊系统模块输出的炉温预报值和使炉温最佳的操作变量值作为结果显示模块的输入；智能优化模块，用于求解模型输出误差最小问题；模型更新模块，用于按设定的采样时间间隔，采集现场智能仪表信号。本发明实现了炉温的准确控制、抑制了系统中的误差，避免出现炉温过低或过高。

主权项

一种智能支持向量机的农药废液焚烧炉炉温最佳化系统，包括与工业过程对象连接的现场智能仪表、DCS系统以及上位机，所述的DCS系统包括控制站和数据库；所述现场智能仪表与DCS系统连接，所述DCS系统与上位机连接，其特征在于：所述的上位机包括：标准化处理模块，用于将从DCS数据库输入的模型训练样本进行预处理，对训练样本中心化，即减去样本的平均值，然后对其进行标准化：计算均值：$\overline{TX}=\frac{1}{N}\underset{i=1}{\overset{N}{\Sigma }}{\mathrm{TX}}_i---\left(1\right)$计算方差：${\sigma }_x^2=\frac{1}{N-1}\underset{i=1}{\overset{N}{\Sigma }}({\mathrm{TX}}_i-\overline{TX})---\left(2\right)$标准化：$X=\frac{TX-\overline{TX}}{{\sigma }_x}---\left(3\right)$其中，TX_i为第i个训练样本，是从DCS数据库中采集的生产正常时的关键变量、炉温和使炉温最佳化的操作变量的数据，N为训练样本数，为训练样本的均值，X为标准化后的训练样本；σ_x表示训练样本的标准差，σ²_x表示训练样本的方差；模糊系统模块，对从标准化处理模块传过来的标准化后的训练样本X，进行模糊化；设模糊系统中有c^*个模糊群，模糊群k、j的中心分别为v_k、v_j，则第i个标准化后的训练样本X_i对于模糊群k的隶属度μ_ik为：${\mu }_{ik}={\left(\underset{j=1}{\overset{c^{*}}{\Sigma }}{\left(\frac{\left|\right|X_i-v_k\left|\right|}{\left|\right|X_i-v_j\left|\right|}\right)}^{\frac{2}{n-1}}\right)}^{-1}---\left(4\right)$式中，n为模糊分类过程中需要的分块矩阵指数，取作2，||·||为范数表达式；使用以上隶属度值或者它的变形以获得新的输入矩阵，对于模糊群k，其输入矩阵变形为：Φ_ik(X_i,μ_ik)＝[1 func(μ_ik) X_i]             (5)其中func(μ_ik)为隶属度值μ_ik的变形函数，取或exp(μ_ik)，Φ_ik(X_i,μ_ik)表示第i个输入变量X_i及其模糊群k的隶属度μ_ik所对应的新的输入矩阵；支持向量机作为模糊系统的局部方程；模糊系统的局部方程对每个模糊群进行优化拟合；设模型训练样本的第i个目标输出为O_i，加权重的支持向量机通过变换把拟合问题等价于如下二次规划问题：$\min R(w,\xi )=\frac{1}{2}{w}^{T}w+\frac{1}{2}{\mathrm{\gamma \Sigma }}_{i=1}^N{\omega }_i{\xi }_i^2---\left(6\right)$同时定义拉格朗日函数：其中，R(w,ξ)是优化问题的目标函数，minR(w,ξ)是优化问题的目标函数的最小值，N是训练样本数，ξ＝{ξ₁,…,ξ_N}是松弛变量，ξ_i是松弛变量的第i个分量，w是支持向量机超平面的法向量，b是相应的偏移量，而ω_i,i＝1,…,N和γ分别是加权支持向量机的权重和惩罚因子，α_i，i＝1,…,N是对应的拉格朗日乘子的第i个分量，是非线性映射函数，μ_ik表示第i个标准化后的训练样本X_i对于模糊群k的隶属度，Φ_ik(X_i,μ_ik)表示第i个输入变量X_i及其模糊群k的隶属度μ_ik所对应的新的输入矩阵；由(6)(7)(8)式可推导出模糊群k在训练样本i的输出为：${\hat{y}}_{ik}=\underset{m=1}{\overset{N}{\Sigma }}{\alpha }_m\times K<{\Phi }_{im}(X_m,{\mu }_{mk}),{\Phi }_{ik}(X_i,{\mu }_{ik})>+b---\left(9\right)$其中，为模糊群k在训练样本i的输出，K<·>是加权支持向量机的核函数，这里K<·>取线性核函数；α_m，m＝1,…,N是对应的拉格朗日乘子；μ_mk表示第m个训练样本X_m对于模糊群k的隶属度，Φ_mk(X_m,μ_mk)表示第m个输入变量X_m及其模糊群k的隶属度μ_mk所对应的新的输入矩阵；由反模糊方法中的重心法得到最后的模糊系统的输出：${\hat{y}}_i=\frac{{\Sigma }_{k=1}^{c^{*}}{\mu }_{ik}{\hat{y}}_{ik}}{{\Sigma }_{k=1}^{c^{*}}{\mu }_{ik}}---\left(10\right)$式中，为模糊系统的输出，为模糊群k在训练样本i的输出智能优化模块，用于采用粒子群算法对模糊系统的局部方程的惩罚因子和误差容限值进行优化，具体实现步骤如下：①确定粒子群的优化参数为模糊系统的局部方程的惩罚因子和误差容限值、粒子群个体数目popsize、最大循环寻优次数iter_max、第p个粒子的初始位置r_p、初始速度v_p、局部最优值Lbest_p以及整个粒子群的全局最优值Gbest；②设定优化目标函数，将其转换为适应度，对每个模糊系统的局部方程进行评价；通过相应的误差函数计算适应度函数，并认为误差大的粒子适应度小，粒子p的适应度函数表示为：f_p＝1/(E_p+1)                        (11)式中，E_p是模糊系统的误差函数，表示为：$E_p=\sqrt{\frac{1}{N}\underset{i=1}{\overset{N}{\Sigma }}{({\hat{y}}_i-O_i)}^2}---\left(12\right)$式中，是模糊系统的预测输出，O_i为模糊系统的目标输出；③按照如下公式，循环更新每个粒子的速度和位置，v_p(iter+1)＝ω×v_p(iter)+m₁a₁(Lbest_p‑r_p(iter))+m₂a₂(Gbest‑r_p(iter))                                                        (13)r_p(iter+1)＝r_p(iter)+v_p(iter+1)                       (14)式中，v_p表示更新粒子p的速度，r_p表示更新粒子p的位置，Lbest_p表示第p粒子的局部最优值，Gbest为整个粒子群的全局最优值，iter表示循环次数，ω是粒子群算法中的惯性权重，m₁、m₂是对应的加速系数，a₁、a₂是[0,1]之间的随机数；④对于粒子p，如果新的适应度大于原来的局部最优值，更新粒子的局部最优值：Lbest_p＝f_p                               (15)⑤如果粒子p的局部最优值Lbest_p大于原来的粒子群全局最优值Gbest，更新原来的粒子群全局最优值Gbest:Gbest＝Lbest_p                           (16)⑥判断是否满足性能要求，若是，结束寻优，得到一组优化的模糊系统的局部方程的参数；否则返回步骤③，继续迭代寻优，直至达到最大循环寻优次数iter_max；Gbest即为对应于第i个标准化后的训练样本X_i的炉温预报值和使炉温最佳的操作变量值；所述的上位机还包括：模型更新模块，用于按设定的采样时间间隔，采集现场智能仪表信号，将得到的实测炉温与系统预报值比较，如果相对误差大于10％或炉温超出生产正常上下限范围，则将DCS数据库中生产正常时的使炉温最佳的新数据加入训练样本数据，更新软测量模型；结果显示模块，用于将得到的炉温预报值和使炉温最佳的操作变量值传给DCS系统，在DCS的控制站显示，并通过DCS系统和现场总线传递到现场操作站进行显示；同时，DCS系统将所得到的使炉温最佳的操作变量值作为新的操作变量设定值，自动执行炉温最佳化操作；信号采集模块，用于依照设定的每次采样的时间间隔，从数据库中采集数据；所述的关键变量包括进入焚烧炉的废液流量、进入焚烧炉的空气流量和进入焚烧炉的燃料流量；所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。