CN103675007B - Rbf网络的工业熔融指数软测量仪表及方法 - Google Patents

Rbf网络的工业熔融指数软测量仪表及方法 Download PDF

Info

Publication number
CN103675007B
CN103675007B CN201310433060.8A CN201310433060A CN103675007B CN 103675007 B CN103675007 B CN 103675007B CN 201310433060 A CN201310433060 A CN 201310433060A CN 103675007 B CN103675007 B CN 103675007B
Authority
CN
China
Prior art keywords
sigma
training sample
model
melt index
fuzzifying equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201310433060.8A
Other languages
English (en)
Other versions
CN103675007A (zh
Inventor
刘兴高
张明明
李见会
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201310433060.8A priority Critical patent/CN103675007B/zh
Publication of CN103675007A publication Critical patent/CN103675007A/zh
Application granted granted Critical
Publication of CN103675007B publication Critical patent/CN103675007B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种RBF网络的工业熔融指数软测量仪表及方法。该方法使用多个RBF神经网络对训练样本进行训练,然后对每个RBF神经网络的输出结果进行模糊化处理,提高了系统的预报精度。在本发明中,用于测量易测变量的现场智能仪表、控制站与DCS数据库连接,DCS数据库与所述软测量模型的输入端连接,所述RBF网络的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接;本发明具有在线测量、计算速度快、输入数据集自动更新、抗噪声能力强的特点。

Description

RBF网络的工业熔融指数软测量仪表及方法
技术领域
本发明设计软测量仪表及方法,尤其涉及一种RBF网络的工业熔融指数软测量仪表和方法。
背景技术
聚丙烯是一种由丙烯聚合而成的半结晶的热塑性塑料,具有较高的耐冲击性,机械性质强韧,抗多种有机溶剂和酸碱腐蚀,在工业界有广泛的应用,是平常最常见的高分子材料之一。熔融指数(MI)是聚丙烯生产中确定最终产品牌号的重要质量指标之一,它决定了产品的不同用途。熔融指数的精确、及时的测量,对生产和科研,都有非常重要的作用和指导意义。然而,熔融指数的在线分析测量目前仍然很难做到,缺乏熔融指数的在线分析仪是制约聚丙烯产品质量的一个主要问题。MI只能通过人工取样、离线化验分析获得,而且一般每2-4小时分析一次,时间滞后大,难以满足生产实时控制的要求。
近年来人工神经网络,尤其是RBF神经网络,在系统优化方面取得了不错的效果。神经网络具有很强的自适应、自组织、自学习的能力以及大规模并行运算的能力。但在实际应用中,神经网络也暴露了一些自身固有的缺陷:权值的初始化是随机的,易陷入局部极小;学习过程中隐含层的节点数目和其他参数的选择只能根据经验和实验来选择;收敛时间过长、鲁棒性差等。其次,工业现场采集到的DCS数据也因为噪音、人工操作误差等带有一定的不确定误差,所以使用确定性强的人工神经网络的预报模型一般推广能力不强。
1965年美国数学家L.Zadeh首先提出了模糊集合的概念。随后模糊逻辑以其更接近于日常人们的问题和语意陈述的方式,开始代替坚持所有事物都可以用二元项表示的经典逻辑。模糊逻辑迄今已经成功应用在了工业的多个领域之中,例如家电、工业控制等领域。2003年,Demirci提出了模糊方程的概念,通过使用模糊隶属度矩阵和和其变形构建一个新的输入矩阵,接着在局部方程中以反模糊方法中的重心法得出解析值作为最后的输出。对于丙烯聚合生产过程中熔融指数的软测量,考虑到工业生产过程中的噪音影响以及操作误差,可以使用模糊逻辑的模糊性能降低误差对整个预报精度的影响。
发明内容
为了克服已有的丙烯聚合生产过程的测量精度不高、对噪声敏感度低、参数选取难度大的不足,本发明提供一种在线测量、计算速度快、模型自动更新、抗噪声能力强的RBF网络的工业熔融指数软测量仪表和方法。
一种RBF网络的工业熔融指数软测量仪表,包括用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库以及熔融指数软测量值显示仪,所述现场智能仪表、控制站与DCS数据库连接,所述软测量仪表还包括RBF网络的工业熔融指数软测量模型,所述DCS数据库与所述RBF网络的工业熔融指数软测量模型的输入端连接,所述RBF网络的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接,所述RBF网络的工业熔融指数软测量模型包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值: TX - = 1 N Σ i = 1 N TX i - - - ( 1 )
计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX - ) - - - ( 2 )
标准化: X = TX - TX - σ x - - - ( 3 )
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则第i个标准化后的训练样本Xi对于模糊群k的隶属度μik为:
μ ik = ( Σ j = 1 c * ( | | X i - v k | | | | X i - v j | | ) 2 m - 1 ) - 1 - - - ( 4 )
式中,m为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xiik)=[1 func(μik)Xi]               (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xiik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
RBF神经网络作为模糊方程系统的局部方程,对每个模糊群进行优化拟合。设第k个RBF神经网络模糊方程的的输出是:
y ^ ik = Σ l w lk Φ lk ( | | X i - C lk | | ) - - - ( 6 )
其中Clk和wlk是第k个RBF神经网络模糊方程第l个结点的中心和输出权值,Φlk(||Xi-Clk||)是第k个RBF神经网络模糊方程第l个结点的径向基函数,由下式确定:
Φ lk ( | | X i - C lk | | ) = exp ( - ( | | X i - C lk | | ) 2 2 × σ lk ) - - - ( 7 )
其中σlk是对应的径向基函数的高斯宽度,||·||为范数表达式。最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
y ^ i = Σ k = 1 c * μ ik y ^ ik Σ k = 1 c * μ ik - - - ( 8 )
式中,为第k个RBF神经网络模糊方程的的输出。
作为优选的一种方案,所述RBF网络的工业熔融指数软测量模型还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊方程模型。
一种RBF网络的工业熔融指数软测量方法,所述软测量方法具体实现步骤如下:
1)、对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1。该处理采用以下算式过程来完成:
2.1)计算均值: TX - = 1 N Σ i = 1 N TX i - - - ( 1 )
2.2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX - ) - - - ( 2 )
2.3)标准化: X = TX - TX - σ x - - - ( 3 )
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
3)、对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则第i个标准化后的训练样本Xi对于模糊群k的隶属度μik为:
μ ik = ( Σ j = 1 c * ( | | X i - v k | | | | X i - v j | | ) 2 m - 1 ) - 1 - - - ( 4 )
式中,m为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xiik)=[1 func(μik)Xi]                 (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xiik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
RBF神经网络作为模糊方程系统的局部方程,对每个模糊群进行优化拟合。设第k个RBF神经网络模糊方程的的输出是:
y ^ ik = Σ l w lk Φ lk ( | | X i - C lk | | ) - - - ( 6 )
其中Clk和wlk是第k个RBF神经网络模糊方程第l个结点的中心和输出权值,Φlk(||Xi-Clk||)是第k个RBF神经网络模糊方程第l个结点的径向基函数,由下式确定:
Φ lk ( | | X i - C lk | | ) = exp ( - ( | | X i - C lk | | ) 2 2 × σ lk ) - - - ( 7 )
其中σlk是对应的径向基函数的高斯宽度,||·||为范数表达式。最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
y ^ i = Σ k = 1 c * μ ik y ^ ik Σ k = 1 c * μ ik - - - ( 8 )
式中,为第k个RBF神经网络模糊方程的的输出。
作为优选的一种方案:所述软测量方法还包括以下步骤:4)、定期将离线化验数据输入到训练集中,更新模糊方程模型。
本发明的技术构思为:对丙烯聚合生产过程的重要质量指标熔融指数进行在线软测量,克服已有的聚丙稀熔融指数测量仪表测量精度不高、对噪声敏感度低、推广性能差的不足,以得到较高精度的软测量结果。此模型相对于已有的熔融指数软测量模型有以下优点:此软测量仪表,相对于已有的熔融指数软测量仪表,具有以下优点:(1)减小了噪声和人工操作误差对预报精度的影响;(2)增强了预报性能;(3)提高了稳定性。
本发明的有益效果主要表现在:1、在线测量;2、模型自动更新;3、抗噪声干扰能力强、4、精度高。
附图说明
图1是RBF网络的工业熔融指数软测量仪表及方法的基本结构示意图;
图2是RBF网络的工业熔融指数软测量模型结构示意图。
具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例1
参照图1、图2,一种RBF网络的工业熔融指数软测量仪表,包括丙烯聚合生产过程1、用于测量易测变量的现场智能仪表2、用于测量操作变量的控制站3、存放数据的DCS数据库4以及熔融指数软测量值显示仪6,所述现场智能仪表2、控制站3与丙烯聚合生产过程1连接,所述现场智能仪表2、控制站3与DCS数据库4连接,所述软测量仪表还包括RBF神经网络模糊方程的软测量模型5,所述DCS数据库4与所述RBF网络的工业熔融指数软测量模型5的输入端连接,所述RBF网络的工业熔融指数软测量模型5的输出端与熔融指数软测量值显示仪6连接,所述RBF网络的工业熔融指数软测量模型包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值: TX - = 1 N Σ i = 1 N TX i - - - ( 1 )
计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX - ) - - - ( 2 )
标准化: X = TX - TX - σ x - - - ( 3 )
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则第i个标准化后的训练样本Xi对于模糊群k的隶属度μik为:
μ ik = ( Σ j = 1 c * ( | | X i - v k | | | | X i - v j | | ) 2 m - 1 ) - 1 - - - ( 4 )
式中,m为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xiik)=[1 func(μik)Xi]                   (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xiik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
RBF神经网络作为模糊方程系统的局部方程,对每个模糊群进行优化拟合。设第k个RBF神经网络模糊方程的的输出是:
y ^ ik = Σ l w lk Φ lk ( | | X i - C lk | | ) - - - ( 6 )
其中Clk和wlk是第k个RBF神经网络模糊方程第l个结点的中心和输出权值,Φlk(||Xi-Clk||)是第k个RBF神经网络模糊方程第l个结点的径向基函数,由下式确定:
Φ lk ( | | X i - C lk | | ) = exp ( - ( | | X i - C lk | | ) 2 2 × σ lk ) - - - ( 7 )
其中σlk是对应的径向基函数的高斯宽度,||·||为范数表达式。最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
y ^ i = Σ k = 1 c * μ ik y ^ ik Σ k = 1 c * μ ik - - - ( 8 )
式中,为第k个RBF神经网络模糊方程的的输出。
作为优选的一种方案,所述RBF网络的工业熔融指数软测量模型还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊方程模型。
根据反应机理以及流程工艺分析,考虑到聚丙烯生产过程中对熔融指数产生影响的各种因素,取实际生产过程中常用的九个操作变量和易测变量作为建模变量,有:三股丙稀进料流率,主催化剂流率,辅催化剂流率,釜内温度、压强、液位,釜内氢气体积浓度。表1列出了作为软测量模型5输入的9个建模变量,分别为釜内温度(T)、釜内压力(p)、釜内液位(L)、釜内氢气体积浓度(Xv)、3股丙烯进料流率(第一股丙稀进料流率f1,第二股丙稀进料流率f2,第三股丙稀进料流率f3)、2股催化剂进料流率(主催化剂流率f4,辅催化剂流率f5)。反应釜中的聚合反应是反应物料反复混合后参与反应的,因此模型输入变量涉及物料的过程变量采用前若干时刻的平均值。此例中数据采用前一小时的平均值。熔融指数离线化验值作为软测量模型5的输出变量。通过人工取样、离线化验分析获得,每4小时分析采集一次。
现场智能仪表2及控制站3与丙烯聚合生产过程1相连,与DCS数据库4相连;软测量模型5与DCS数据库及软测量值显示仪6相连。现场智能仪表2测量丙烯聚合生产对象的易测变量,将易测变量传输到DCS数据库4;控制站3控制丙烯聚合生产对象的操作变量,将操作变量传输到DCS数据库4。DCS数据库4中记录的变量数据作为RBF网络的工业熔融指数软测量模型5的输入,软测量值显示仪6用于显示RBF网络的工业熔融指数软测量模型5的输出,即软测量值。
表1RBF网络的工业熔融指数软测量模型所需建模变量
变量符号 变量含义 变量符号 变量含义
T 釜内温度 f1 第一股丙稀进料流率
p 釜内压强 f2 第二股丙稀进料流率
L 釜内液位 f3 第三股丙稀进料流率
Xv 釜内氢气体积浓度 f4 主催化剂流率
f5 辅催化剂流率
RBF网络的工业熔融指数软测量模型5,包括以下3个部分:
数据预处理模块7,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值: TX - = 1 N Σ i = 1 N TX i - - - ( 1 )
计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX - ) - - - ( 2 )
标准化: X = TX - TX - σ x - - - ( 3 )
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
模糊方程模块8,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则第i个标准化后的训练样本Xi对于模糊群k的隶属度μik为:
μ ik = ( Σ j = 1 c * ( | | X i - v k | | | | X i - v j | | ) 2 m - 1 ) - 1 - - - ( 4 )
式中,m为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xiik)=[1 func(μik)Xi]            (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xiik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
RBF神经网络作为模糊方程系统的局部方程,对每个模糊群进行优化拟合。设第k个RBF神经网络模糊方程的的输出是:
y ^ ik = Σ l w lk Φ lk ( | | X i - C lk | | ) - - - ( 6 )
其中Clk和wlk是第k个RBF神经网络模糊方程第l个结点的中心和输出权值,Φlk(||Xi-Clk||)是第k个RBF神经网络模糊方程第l个结点的径向基函数,由下式确定:
Φ lk ( | | X i - C lk | | ) = exp ( - ( | | X i - C lk | | ) 2 2 × σ lk ) - - - ( 7 )
其中σlk是对应的径向基函数的高斯宽度,||·||为范数表达式。最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
y ^ i = Σ k = 1 c * μ ik y ^ ik Σ k = 1 c * μ ik - - - ( 8 )
式中,为第k个RBF神经网络模糊方程的的输出。
模型更新模块9,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊方程模型。
实施例2
参照图1、图2,一种基于RBF神经网络模糊方程模型的工业聚丙烯生产熔融指数软测量方法,所述软测量方法具体实现步骤如下:
1)、对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1。该处理采用以下算式过程来完成:
2.1)计算均值: TX - = 1 N Σ i = 1 N TX i - - - ( 1 )
2.2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX - ) - - - ( 2 )
2.3)标准化: X = TX - TX - σ x - - - ( 3 )
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
3)、对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则第i个标准化后的训练样本Xi对于模糊群k的隶属度μik为:
μ ik = ( Σ j = 1 c * ( | | X i - v k | | | | X i - v j | | ) 2 m - 1 ) - 1 - - - ( 4 )
式中,m为模糊分类过程中需要的分块矩阵指数,通常取作2,||·||为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xiik)=[1 func(μik)Xi]            (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xiik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
RBF神经网络作为模糊方程系统的局部方程,对每个模糊群进行优化拟合。设第k个RBF神经网络模糊方程的的输出是:
y ^ ik = Σ l w lk Φ lk ( | | X i - C lk | | ) - - - ( 6 )
其中Clk和wlk是第k个RBF神经网络模糊方程第l个结点的中心和输出权值,Φlk(||Xi-Clk||)是第k个RBF神经网络模糊方程第l个结点的径向基函数,由下式确定:
Φ lk ( | | X i - C lk | | ) = exp ( - ( | | X i - C lk | | ) 2 2 × σ lk ) - - - ( 7 )
其中σlk是对应的径向基函数的高斯宽度,||·||为范数表达式。最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
y ^ i = Σ k = 1 c * μ ik y ^ ik Σ k = 1 c * μ ik - - - ( 8 )
式中,为第k个RBF神经网络模糊方程的的输出。
作为优选的一种方案:所述软测量方法还包括以下步骤:4)、定期将离线化验数据输入到训练集中,更新模糊方程模型。
本实施例的方法具体实现步骤如下:
步骤1:对丙烯聚合生产过程对象1,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入。操作变量和易测变量由DCS数据库4获得。
步骤2:对样本数据进行预处理,由数据预处理模块7完成。
步骤3:基于模型训练样本数据建立模糊方程模型8。输入数据如步骤2所述获得,输出数据由离线化验获得。
步骤4:模型更新模块9定期将离线化验数据输入到训练集中,更新模糊方程模型,基于RBF神经网络模糊方程模型的软测量模型5建立完成。
步骤5:熔融指数软测量值显示仪6显示RBF网络的工业熔融指数软测量模型5的输出,完成对工业聚丙烯生产熔融指数软测量的显示。

Claims (2)

1.一种RBF网络的工业熔融指数软测量仪表,包括用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库以及熔融指数软测量值显示仪,所述现场智能仪表、控制站与DCS数据库连接,其特征在于:所述软测量仪表还包括RBF网络的软测量模型,所述DCS数据库与所述RBF网络的工业熔融指数软测量模型的输入端连接,所述RBF网络的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接,所述RBF网络的工业熔融指数软测量模型包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值: TX ‾ = 1 N Σ i = 1 N TX i - - - ( 1 )
计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) - - - ( 2 )
标准化: X = TX - TX ‾ σ x - - - ( 3 )
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本;σx表示训练样本的标准差,σ2 x表示训练样本的方差;
模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化;设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则第i个标准化后的训练样本Xi对于模糊群k的隶属度μik为:
μ ik = ( Σ j = 1 c * ( | | X i - v k | | | | X i - v j | | ) 2 m - 1 ) - 1 - - - ( 4 )
式中,m为模糊分类过程中需要的分块矩阵指数,取作2,||·||为范数表达式;
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xiik)=[1 func(μik) Xi]    (5)
其中func(μik)为隶属度值μik的变形函数,取或exp(μik),Φik(Xiik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵;
RBF神经网络作为模糊方程系统的局部方程,对每个模糊群进行优化拟合;设第k个RBF神经网络模糊方程的的输出是:
y ^ ik = Σ l w lk Φ lk ( | | X i - C lk | | ) - - - ( 6 )
其中Clk和wlk是第k个RBF神经网络模糊方程第l个结点的中心和输出权值,Φlk(||Xi-Clk||)是第k个RBF神经网络模糊方程第l个结点的径向基函数,由下式确定:
Φ lk ( | | X i - C lk | | ) = exp ( - ( | | X i - C lk | | ) 2 2 × σ lk ) - - - ( 7 )
其中σlk是对应的径向基函数的高斯宽度,||·||为范数表达式;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
y ^ i = Σ k = 1 c * μ ik y ^ ik Σ k = 1 c * μ ik - - - ( 8 )
式中,为第k个RBF神经网络模糊方程的的输出;
所述的RBF网络的工业熔融指数软测量模型还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊方程模型。
2.一种用如权利要求1所述的RBF网络的工业熔融指数软测量仪表实现的软测量方法,其特征在于:所述软测量方法主要包括以下步骤:
1)、对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1;该处理采用以下算式过程来完成:
2.1)计算均值: TX ‾ = 1 N Σ i = 1 N TX i - - - ( 1 )
2.2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) - - - ( 2 )
2.3)标准化: X = TX - TX ‾ σ x - - - ( 3 )
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本;σx表示训练样本的标准差,σ2 x表示训练样本的方差;
3)、对从数据预处理模块传过来的训练样本,进行模糊化;设模糊方程系统中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则第i个标准化后的训练样本Xi对于模糊群k的隶属度μik为:
μ ik = ( Σ j = 1 c * ( | | X i - v k | | | | X i - v j | | ) 2 m - 1 ) - 1 - - - ( 4 )
式中,m为模糊分类过程中需要的分块矩阵指数,取作2,||·||为范数表达式;
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xiik)=[1 func(μik) Xi]   (5)
其中func(μik)为隶属度值μik的变形函数,取或exp(μik),Φik(Xiik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵;
RBF神经网络作为模糊方程系统的局部方程,对每个模糊群进行优化拟合;设第k个RBF神经网络模糊方程的的输出是:
y ^ ik = Σ l w lk Φ lk ( | | X i - C lk | | ) - - - ( 6 )
其中Clk和wlk是第k个RBF神经网络模糊方程第l个结点的中心和输出权值,Φlk(||Xi-Clk||)是第k个RBF神经网络模糊方程第l个结点的径向基函数,由下式确定:
Φ lk ( | | X i - C lk | | ) = exp ( - ( | | X i - C lk | | ) 2 2 × σ lk ) - - - ( 7 )
其中σlk是对应的径向基函数的高斯宽度,||·||为范数表达式;最后,由反模糊方法中的重心法得到最后的模糊方程系统的输出:
y ^ i = Σ k = 1 c * μ ik y ^ ik Σ k = 1 c * μ ik - - - ( 8 )
式中,为第k个RBF神经网络模糊方程的的输出;
所述软测量方法还包括以下步骤:4)、定期将离线化验数据输入到训练集中,更新模糊方程模型。
CN201310433060.8A 2013-09-22 2013-09-22 Rbf网络的工业熔融指数软测量仪表及方法 Expired - Fee Related CN103675007B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310433060.8A CN103675007B (zh) 2013-09-22 2013-09-22 Rbf网络的工业熔融指数软测量仪表及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310433060.8A CN103675007B (zh) 2013-09-22 2013-09-22 Rbf网络的工业熔融指数软测量仪表及方法

Publications (2)

Publication Number Publication Date
CN103675007A CN103675007A (zh) 2014-03-26
CN103675007B true CN103675007B (zh) 2015-08-26

Family

ID=50313212

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310433060.8A Expired - Fee Related CN103675007B (zh) 2013-09-22 2013-09-22 Rbf网络的工业熔融指数软测量仪表及方法

Country Status (1)

Country Link
CN (1) CN103675007B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113030156B (zh) * 2021-03-13 2023-02-24 宁波大学科学技术学院 一种基于非线性慢特征模型的聚丙烯熔融指数软测量方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101315556A (zh) * 2008-06-25 2008-12-03 浙江大学 基于Chaos混沌寻优的丙烯聚合生产过程最优软测量仪表及方法
WO2010015465A1 (de) * 2008-08-08 2010-02-11 Endress+Hauser Gmbh+Co.Kg Diagnoseverfahren eines prozessautomatisierungssystem

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101315556A (zh) * 2008-06-25 2008-12-03 浙江大学 基于Chaos混沌寻优的丙烯聚合生产过程最优软测量仪表及方法
WO2010015465A1 (de) * 2008-08-08 2010-02-11 Endress+Hauser Gmbh+Co.Kg Diagnoseverfahren eines prozessautomatisierungssystem

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MUSTAFA DEMIRCI.FOUNDATIONS OF FUZZY FUNCTIONS AND VAGUE ALGEBRA BASED ON MANY-VALUED EQUIVALENCE RELATIONS, PART I: FUZZY FUNCTIONS AND THEIR APPLICATIONS.《International Journal of General Systems》.2003,第32卷(第2期),123-155. *
张志猛等.遗传算法的聚丙烯熔融指数最优预报.《计算机与应用化学》.2011,第28卷(第7期),843-846. *
李九宝等.基于PSO-SA算法的聚丙烯熔融指数预报.《化工学报》.2010,第61卷(第8期),1955-1959. *
李炜等.基于模糊RBF 神经网络的管道泄漏检测方法.《兰州理工大学学报》.2009,第35卷(第6期),79-83. *
车荣杰.聚丙烯熔融指数软测量方法研究.《CNKI中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》.2010,(第03期),41-60. *

Also Published As

Publication number Publication date
CN103675007A (zh) 2014-03-26

Similar Documents

Publication Publication Date Title
CN101315557B (zh) 基于遗传算法优化bp神经网络的丙烯聚合生产过程最优软测量仪表及方法
CN103674778B (zh) Rbf粒子群优化的工业熔融指数软测量仪表和方法
CN101863088B (zh) 一种橡胶混炼过程中门尼粘度的预报方法
CN101382801B (zh) 基于ega优化的丙烯聚合生产过程最优软测量仪表及方法
CN103675011A (zh) 最优支持向量机的工业熔融指数软测量仪表和方法
CN103675006A (zh) 最小二乘的工业熔融指数软测量仪表和方法
CN103823430B (zh) 智能加权丙烯聚合生产过程最优软测量系统和方法
CN103839103B (zh) 丙烯聚合生产过程bp最优预报系统和方法
CN102621953B (zh) 一种橡胶硬度的在线自动质量监控和预测模型更新的方法
CN103675005B (zh) 最优模糊网络的工业熔融指数软测量仪表及方法
CN103675009A (zh) 模糊方程的工业熔融指数软测量仪表及方法
CN103675007B (zh) Rbf网络的工业熔融指数软测量仪表及方法
CN103675012A (zh) Bp粒子群优化的工业熔融指数软测量仪表和方法
CN103675010A (zh) 支持向量机的工业熔融指数软测量仪表及方法
CN103824121A (zh) 基于多模群智的丙烯聚合生产过程最优预报系统和方法
CN103838209B (zh) 丙烯聚合生产过程自适应最优预报系统和方法
CN109856971A (zh) 基于群智寻优的丙烯聚合生产过程最优在线预报系统
CN103838958A (zh) 模糊智能最优丙烯聚合生产过程最优软测量仪表和方法
CN103630568A (zh) Bp网络的工业熔融指数软测量仪表及方法
CN103838957A (zh) 丙烯聚合生产过程径向基最优软测量仪表和方法
CN103823369B (zh) 基于bp多模网络的丙烯聚合生产过程预报系统和方法
CN103678953A (zh) 基于贝叶斯组合神经网络的生物发酵产量在线预报方法
CN103675008A (zh) 加权模糊的工业熔融指数软测量仪表及方法
CN103472867A (zh) 支持向量机的农药生产废液焚烧炉炉温最佳化系统及方法
CN103838205A (zh) Bp全局最优丙烯聚合生产过程最优软测量仪表和方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150826

Termination date: 20180922

CF01 Termination of patent right due to non-payment of annual fee