[发明专利]最优模糊网络的工业熔融指数软测量仪表及方法

权利要求书

1.一种最优模糊网络的工业熔融指数软测量仪表，包括用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库以及熔融指数软测量值显示仪，所述现场智能仪表、控制站与DCS数据库连接，其特征在于：所述软测量仪表还包括最优模糊网络的工业熔融指数软测量模型，所述DCS数据库与所述最优模糊网络的工业熔融指数软测量模型的输入端连接，所述最优模糊网络的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接，所述最优模糊网络的工业熔融指数软测量模型包括：

数据预处理模块，用于将从DCS数据库输入的模型训练样本进行预处理，使得训练样本的均值为0，方差为1，该处理采用以下算式过程来完成：

计算均值：TX‾=1NΣi=1NTXi---(1)]]>

计算方差：σx2=1N-1Σi=1N(TXi-TX‾)---(2)]]>

标准化：X=TX-TX‾σx---(3)]]>

其中，TX_i为第i个训练样本，N为训练样本数，为训练样本的均值，X为标准化后的训练样本。σ_x表示训练样本的标准差，σ²_x表示训练样本的方差。

模糊神经网络模块，对从数据预处理模块传过来的输入变量，进行模糊推理和建立模糊规则。对从数据预处理模块传过来的经过预处理过的训练样本X进行模糊分类，得到模糊规则库中每个模糊聚类的中心和宽度。设第p个标准化后的训练样本X_p=[X_p1,…,X_pn]，其中n是输入变量的个数。

设模糊神经网络有R个模糊规则，为了求得每个模糊规则对于训练样本X_p的每个输入变量X_pj,j=1,…,n，下面的模糊化方程将求出其对第i个模糊规则的隶属度：

Mij=exp{-(Xpj-mij)2σij2}---(4)]]>

其中m_ij和σ_ij分别表示第i个模糊规则的第j个高斯成员函数的中心和宽度，由模糊聚类求得。

设标准化后的训练样本X_p对模糊规则i的适应度为μ⁽ⁱ⁾(X_p)，则μ⁽ⁱ⁾(X_p)的大小可由下式决定：

μ(i)(Xp)=Πj=1nMij(Xp)=exp{-Σj=1n(Xpj-mij)2σij2}---(5)]]>

求得输入训练样本对于每个规则的适应度之后，模糊神经网络对模糊规则输出进行推导以得到最后的解析解。在常用的模糊神经网络结构中，每个模糊规则推导的过程都可以表示为如下：首先求得训练样本中所有输入变量的线性乘积和，然后用此线性乘积和与规则的适用度μⁱ(X_p)相乘，得到最终的每条模糊规则的输出。模糊规则i的推导输出可以表示如下：

f(i)=μ(i)(Xp)×(Σj=1naij×Xpj+ai0)---(6)]]>

y^p=Σi=1Rf(i)+b=Σi=1R[μi(Xp)×(Σj=1naij×Xpj+ai0)]+b---(7)]]>

式中，f⁽ⁱ⁾为第i条模糊规则的输出，是模糊神经网络模型对第p个训练样本的预测输出，a_ij,j=1,…,n是第i条模糊规则中第j个变量的线性系数，a_i0是第i条模糊规则中输入变量线性乘积和的常数项，b是输出偏置量。

支持向量机优化模块，在式（7）中，输入变量线性乘积和中的参数的确定是模糊神经网络使用中用到的一个主要问题，这里我们采用把原有的模糊规则推导输出形式转换为支持向量机优化问题，再使用支持向量机进行线性优化，转换过程如下：

y^p=Σi=1Rf(i)+b]]>

=Σi=1R[μ(i)(Xp)×(Σj=1naij×Xpj+ai0)]+b---(8)]]>

=Σi=1RΣj=0naij×μ(i)(Xp)×Xpj+b]]>

其中X_p0为常数项且恒等于1。令

φ→(Xp)=[μ(1)×Xp0,...,μ(1)×Xpn,......,μ(R)×Xp0,...,μ(R)×Xpn]---(9)]]>

其中，表示原训练样本的转化形式，即把原来的训练样本转换为如上式形式，作为支持向量机的训练样本：

S={(φ→(X1),y1),(φ→(X2),y2),...,(φ→(XN),yN),}---(10)]]>

其中y₁,…,y_N是训练样本的目标输出，取S作为新的输入训练样本集合，那么原有问题可以转化为如下的支持向量机对偶优化问题：

R(ω,b)=γ1NΣp=1NLϵ(yp,f(Xp))+12ωTω---(11)]]>

其中y_p是输入训练样本X_p的目标输出，ω是支持向量机超平面的法向量，f(X_p)是对应于X_p的模型输出，γ是支持向量机的惩罚因子，R(ω,b)是优化问题的目标函数，N是训练样本数，L_ε(y_p,f(X_p))表达式如下：

其中ε是优化问题的误差容限，接下来使用支持向量机求得模糊神经网络的模糊规则最优推导线性参数和对偶优化问题的预报输出：

aij=Σk=1N(αk*-αk)μ(i)Xkj=Σk∈SVN(αk*-αk)μ(i)Xkji=1,...,R;j=0,...,n---(13)]]>

y^p=Σk=1N(αk*-αk)<φ→(X),φ→(Xk)>+b---(14)]]>

其中α_k，(k=1,…,N)分别是y_p-f(X_p)大于0和小于0时对应的拉格朗日乘子,即为对应于第p个标准化后的训练样本X_p的MI预报值。

所述最优模糊网络的工业熔融指数软测量模型还包括：

模型更新模块，用于模型的在线更新，定期将离线化验数据输入到训练集中，更新模糊神经网络模型。

2.一种用如权利要求1所述的基于支持向量机优化模糊神经网络的丙烯聚合生产过程软测量仪表实现的软测量方法，其特征在于：所述软测量方法具体实现步骤如下：

1）、对丙烯聚合生产过程对象，根据工艺分析和操作分析，选择操作变量和易测变量作为模型的输入，操作变量和易测变量由DCS数据库获得；

2）、用于将从DCS数据库输入的模型训练样本进行预处理，使得训练样本的均值为0，方差为1，该处理采用以下算式过程来完成：

计算均值：TX‾=1NΣi=1NTXi---(1)]]>

计算方差：σx2=1N-1Σi=1N(TXi-TX‾)---(2)]]>

标准化：X=TX-TX‾σx---(3)]]>

其中，TX_i为第i个训练样本，N为训练样本数，为训练样本的均值，X为标准化后的训练样本。σ_x表示训练样本的标准差，σ²_x表示训练样本的方差。

3）、对从数据预处理模块传过来的输入变量，进行模糊推理和建立模糊规则。对从数据预处理模块传过来的经过预处理过的训练样本X进行模糊分类，得到模糊规则库中每个模糊聚类的中心和宽度。设第p个标准化后的训练样本X_p=[X_p1,…,X_pn]，其中n是输入变量的个数。

设模糊神经网络有R个模糊规则，为了求得每个模糊规则对于训练样本X_p的每个输入变量X_pj,j=1,…,n，下面的模糊化方程将求出其对第i个模糊规则的隶属度：

Mij=exp{-(Xpj-mij)2σij2}---(4)]]>

其中m_ij和σ_ij分别表示第i个模糊规则的第j个高斯成员函数的中心和宽度，由模糊聚类求得。

设标准化后的训练样本X_p对模糊规则i的适应度为μ⁽ⁱ⁾(X_p)，则μ⁽ⁱ⁾(X_p)的大小可由下式决定：

μ(i)(Xp)=Πj=1nMij(Xp)=exp{-Σj=1n(Xpj-Mij)2σij2}---(5)]]>

求得输入训练样本对于每个规则的适应度之后，模糊神经网络对模糊规则输出进行推导以得到最后的解析解。在常用的模糊神经网络结构中，每个模糊规则推导的过程都可以表示为如下：首先求得训练样本中所有输入变量的线性乘积和，然后用此线性乘积和与规则的适用度μ⁽ⁱ⁾(X_p)相乘，得到最终的每条模糊规则的输出。模糊规则i的推导输出可以表示如下：

f(i)=μ(i)(Xp)×(Σj=1naij×Xpj+ai0)---(6)]]>

y^p=Σi=1Rf(i)+b=Σi=1R[μi(Xp)×(Σj=1naij×Xpj+ai0)]+b---(7)]]>

式中，f⁽ⁱ⁾为第i条模糊规则的输出，是模糊神经网络模型对第p个训练样本的预测输出，a_ij,j=1,…,n是第i条模糊规则中第j个变量的线性系数，a_i0是第i条模糊规则中输入变量线性乘积和的常数项，b是输出偏置量。

4）、在式（7）中，输入变量线性乘积和中的参数的确定是模糊神经网络使用中用到的一个主要问题，这里我们采用把原有的模糊规则推导输出形式转换为支持向量机优化问题，再使用支持向量机进行线性优化，转换过程如下：

y^p=Σi=1Rf(i)+b]]>

=Σi=1R[μ(i)(Xp)×(Σj=1naij×Xpj+ai0)]+b---(8)]]>

=Σi=1RΣj=0naij×μ(i)(Xp)×Xpj+b]]>

其中X_p0为常数项且恒等于1。令

φ→(Xp)=[μ(1)×Xp0,...,μ(1)×Xpn,......,μ(R)×Xp0,...,μ(R)×Xpn]---(9)]]>

其中，表示原训练样本的转化形式，即把原来的训练样本转换为如上式形式，作为支持向量机的训练样本：

S={(φ→(X1),y1),(φ→(X2),y2),...,(φ→(XN),yN),}---(10)]]>

其中y₁,…,y_N是训练样本的目标输出，取S作为新的输入训练样本集合，那么原有问题可以转化为如下的支持向量机对偶优化问题：

R(ω,b)=γ1NΣp=1NLϵ(yp,f(Xp))+12ωTω---(11)]]>

其中y_p是输入训练样本X_p的目标输出，ω是支持向量机超平面的法向量，f(X_p)是对应于X_p的模型输出，γ是支持向量机的惩罚因子，R(ω,b)是优化问题的目标函数，N是训练样本数，L_ε(y_p,f(X_p))表达式如下：

其中ε是优化问题的误差容限，接下来使用支持向量机求得模糊神经网络的模糊规则最优推导线性参数和对偶优化问题的预报输出：

aij=Σk=1N(αk*-αk)μ(i)Xkj=Σk∈SVN(αk*-αk)μ(i)Xkji=1,...,R;j=0,...,n---(13)]]>

y^p=Σk=1N(αk*-αk)<φ→(X),φ→(Xk)>+b---(14)]]>

其中α_k，(k=1,…,N)分别是y_p-f(X_p)大于0和小于0时对应的拉格朗日乘子,即为对应于第p个标准化后的训练样本X_p的MI预报值。

所述软测量方法还包括以下步骤：5）、定期将离线化验数据输入到训练样本集中，更新模糊神经网络模型。