CN103675012B - Bp粒子群优化的工业熔融指数软测量仪表和方法 - Google Patents
Bp粒子群优化的工业熔融指数软测量仪表和方法 Download PDFInfo
- Publication number
- CN103675012B CN103675012B CN201310436856.9A CN201310436856A CN103675012B CN 103675012 B CN103675012 B CN 103675012B CN 201310436856 A CN201310436856 A CN 201310436856A CN 103675012 B CN103675012 B CN 103675012B
- Authority
- CN
- China
- Prior art keywords
- particle
- equation
- iter
- fuzzy
- training sample
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种BP粒子群优化的工业熔融指数软测量仪表及方法。该软测量方法使用BP神经网络作为模糊方程***的局部方程,随后引入粒子群算法对整个软测量模型进行进一步优化。在本发明中,用于测量易测变量的现场智能仪表、控制站与DCS数据库连接,软测量值显示仪包括BP粒子群优化的工业熔融指数软测量模型,DCS数据库与软测量模型的输入端连接,所述BP粒子群优化的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接;本发明具有在线优化、计算速度快、模型自动更新、抗噪声能力强的特点。
Description
技术领域
本发明涉及软测量仪表及方法,尤其涉及一种BP粒子群优化的工业熔融指数软测量仪表及方法。
背景技术
聚丙烯是一种由丙烯聚合而成的半结晶的热塑性塑料,具有较高的耐冲击性,机械性质强韧,抗多种有机溶剂和酸碱腐蚀,在工业界有广泛的应用,是平常最常见的高分子材料之一。熔融指数(MI)是聚丙烯生产中确定最终产品牌号的重要质量指标之一,它决定了产品的不同用途。熔融指数的精确、及时的测量,对生产和科研,都有非常重要的作用和指导意义。然而,熔融指数的在线分析测量目前仍然很难做到,缺乏熔融指数的在线分析仪是制约聚丙烯产品质量的一个主要问题。MI只能通过人工取样、离线化验分析获得,而且一般每2-4小时分析一次,时间滞后大,难以满足生产实时控制的要求。
近年来人工神经网络,尤其是BP神经网络,在***优化方面取得了不错的效果。神经网络具有很强的自适应、自组织、自学习的能力以及大规模并行运算的能力。但在实际应用中,神经网络也暴露了一些自身固有的缺陷:权值的初始化是随机的,易陷入局部极小;学习过程中隐含层的节点数目和其他参数的选择只能根据经验和实验来选择;收敛时间过长、鲁棒性差等。其次,工业现场采集到的DCS数据也因为噪音、人工操作误差等带有一定的不确定误差,所以使用确定性强的人工神经网络的预报模型一般推广能力不强。
1965年美国数学家L.Zadeh首先提出了模糊集合的概念。随后模糊逻辑以其更接近于日常人们的问题和语意陈述的方式,开始代替坚持所有事物都可以用二元项表示的经典逻辑。模糊逻辑迄今已经成功应用在了工业的多个领域之中,例如家电、工业控制等领域。2003年,Demirci提出了模糊方程的概念,通过使用模糊隶属度矩阵和和其变形构建一个新的输入矩阵,接着在局部方程中以反模糊方法中的重心法得出解析值作为最后的输出。对于丙烯聚合生产过程中熔融指数的软测量,考虑到工业生产过程中的噪音影响以及操作误差,可以使用模糊逻辑的模糊性能降低误差对整个预报精度的影响。
粒子群算法,即Particle Swarm Optimization,是由Kennedy和Eberhart教授提出来的一种通过模仿鸟类飞行行为来寻求全局最优的一种生物智能寻优算法,简称PSO。该算法通过群体中粒子间的相互影响,减少了搜索算法陷入局部最优解的风险,具有很好的全局搜索性能。粒子群算法被用来搜索BP神经网络模糊方程的最优参数组合,以达到优化模型的目的。
发明内容
为了克服已有的丙烯聚合生产过程的测量精度不高、对噪声敏感度低、参数选取难度大的不足,本发明提供一种在线测量、计算速度快、模型自动更新、抗噪声能力强的BP粒子群优化的工业熔融指数软测量仪表及方法。
一种BP粒子群优化的工业熔融指数软测量仪表,包括用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库以及熔融指数软测量值显示仪,所述现场智能仪表、控制站与DCS数据库连接,所述软测量仪表还包括BP粒子群优化的工业熔融指数软测量模型,所述DCS数据库与所述BP粒子群优化的工业熔融指数软测量模型的输入端连接,所述BP粒子群优化的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接,所述BP粒子群优化的工业熔融指数软测量模型包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值:
计算方差:
标准化:
其中,TXi为第i个训练样本,,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程***中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,‖·‖为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xi,μik)=[1func(μik)Xi] (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
BP神经网络作为模糊方程***的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为sl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程***的输出:
其中,为第k个BP神经网络模糊方程输出层的预测输出。
粒子群算法优化模块,用于采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体实现步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程***的误差函数,表示为:
式中,是模糊方程***的预测输出,Oi为模糊方程***的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为整个粒子群的全局最优值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
作为优选的一种方案,所述BP粒子群优化的工业熔融指数软测量模型还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊方程模型。
一种BP粒子群优化的工业熔融指数丙烯聚合生产过程软测量方法,所述软测量方法具体实现步骤如下:
1)、对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1。该处理采用以下算式过程来完成:
2.1)计算均值:
2.2)计算方差:
2.3)标准化: 其中,TXi为第i个训练样本,,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
3)、对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程*** 中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,‖·‖为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xi,μik)=[1func(μiz)Xi] (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
BP神经网络作为模糊方程***的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为sl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程***的输出:
其中,为第k个BP神经网络模糊方程输出层的预测输出。
4)、采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体步骤如下:①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误 差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程***的误差函数,表示为:
式中,是模糊方程***的预测输出,Oi为模糊方程***的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为整个粒子群的全局最优值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
作为优选的一种方案:所述软测量方法还包括以下步骤:5)、定期将离线化验数据输入到训练集中,更新模糊方程模型。
本发明的技术构思为:对丙烯聚合生产过程的重要质量指标熔融指数进行在线软测量,克服已有的聚丙稀熔融指数测量仪表测量精度不高、对噪声敏感度低、推广性能差的不足,引入粒子群算法对模糊方程模型进行自动优化,不需要人为经验来多次调整模糊方程中BP神经网络局部方程的参数。此模型相对于已有的熔融指数软测量模型有以下优点:(1)减小了噪声和人工操作误差对模型预报精度的影响;(2)增强了模型的预报性能;(3)对模型的参数进行自动寻优提高了模型的稳定性,降低了模型陷入局部最优的可能性。
本发明的有益效果主要表现在:1、在线测量;2、在线参数自动优化;3、模型自动更 新;4、抗噪声干扰能力强、5、精度高。
附图说明
图1是BP粒子群优化的工业熔融指数软测量仪表及方法的基本结构示意图;
图2是BP粒子群优化的工业熔融指数软测量模型结构示意图。
具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例1
参照图1、图2,一种BP粒子群优化的工业熔融指数软测量仪表,包括丙烯聚合生产过程1、用于测量易测变量的现场智能仪表2、用于测量操作变量的控制站3、存放数据的DCS数据库4以及熔融指数软测量值显示仪6,所述现场智能仪表2、控制站3与丙烯聚合生产过程1连接,所述现场智能仪表2、控制站3与DCS数据库4连接,所述软测量仪表还包括粒子群算法优化BP神经网络模糊方程的软测量模型5,所述DCS数据库4与所述BP粒子群优化的工业熔融指数软测量模型5的输入端连接,所述BP粒子群优化的工业熔融指数软测量模型5的输出端与熔融指数软测量值显示仪6连接,所述BP粒子群优化的工业熔融指数软测量模型包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值:
计算方差:
标准化:
其中,TXi为第i个训练样本,,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程***中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,‖·‖为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xi,μik)=[1func(μik)Xi] (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
BP神经网络作为模糊方程***的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为sl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程***的输出:
其中,为第k个BP神经网络模糊方程输出层的预测输出。
粒子群算法优化模块,用于采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体实现步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程***的误差函数,表示为:
式中,是模糊方程***的预测输出,Oi为模糊方程***的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为整个粒子群的全局最优值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
作为优选的一种方案,所述BP粒子群优化的工业熔融指数软测量模型还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新神经网络模型。
根据反应机理以及流程工艺分析,考虑到聚丙烯生产过程中对熔融指数产生影响的各种因素,取实际生产过程中常用的九个操作变量和易测变量作为建模变量,有:三股丙稀进料流率,主催化剂流率,辅催化剂流率,釜内温度、压强、液位,釜内氢气体积浓度。表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中记录的变量数据作为BP粒子群优化的工业熔融指数软测量模型5的输入,软测量值显示仪6用于显示BP粒子群优化的工业熔融指数软测量模型5的输出,即软测量值。
表1:BP粒子群优化的工业熔融指数软测量模型所需建模变量
变量符号 | 变量含义 | 变量符号 | 变量含义 |
T | 釜内温度 | f1 | 第一股丙稀进料流率 |
p | 釜内压强 | f2 | 第二股丙稀进料流率 |
L | 釜内液位 | f3 | 第三股丙稀进料流率 |
Xv | 釜内氢气体积浓度 | f4 | 主催化剂流率 |
f5 | 辅催化剂流率 |
BP粒子群优化的工业熔融指数软测量模型5,包括以下4个部分:
数据预处理模块7,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值:
计算方差:
标准化:
其中,TXi为第i个训练样本,,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
模糊方程模块8,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程***中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,‖·‖为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xi,μik)=[1func(μik)Xi] (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
BP神经网络作为模糊方程***的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为yl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程***的输出:
其中,为第k个BP神经网络模糊方程输出层的预测输出。
粒子群算法优化模块9,用于采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体实现步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程***的误差函数,表示为:
式中,是模糊方程***的预测输出,Oi为模糊方程***的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为整个粒子群的全局最优值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
模型更新模块10,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊方程模型。
实施例2
参照图1、图2,一种基于粒子群算法优化BP神经网络模糊方程模型的丙烯聚合生产过程软测量方法,所述软测量方法具体实现步骤如下:
1)、对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1。该处理采用以下算式过程来完成:
2.1)计算均值:
2.2)计算方差:
2.3)标准化: 其中,TXi为第i个训练样本,,N为训练样本数,为训练样本的均值,X为标准化后的训练样本。σx表示训练样本的标准差,σ2 x表示训练样本的方差。
3)、对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化。设模糊方程***中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,通常取作2,‖·‖为范数表达式。
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xi,μik)=[1func(μik)Xi] (5)
其中func(μik)为隶属度值μik的变形函数,一般取exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵。
BP神经网络作为模糊方程***的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为sl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程***的输出:
其中,为第k个BP神经网络模糊方程输出层的预测输出。
4)、采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体实现步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest。
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程***的误差函数,表示为:
式中,是模糊方程***的预测输出,Oi为模糊方程***的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为整个粒子群的全局最优值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax。
作为优选的一种方案:所述软测量方法还包括以下步骤:5)、定期将离线化验数据输入到训练集中,更新模糊方程模型。
本实施例的方法具体实施步骤如下:
步骤1:对丙烯聚合生产过程对象1,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入。操作变量和易测变量由DCS数据库4获得。
步骤2:对样本数据进行预处理,由数据预处理模块7完成。
步骤3:基于模型训练样本数据建立初始模糊方程模型8。输入数据如步骤2所述获得,输出数据由离线化验获得。
步骤4:由粒子群算法9优化初始模糊方程模型8的局部BP神经网络方程参数。
步骤5:模型更新模块10定期将离线化验数据输入到训练集中,更新模糊方程模型,基于粒子群算法优化BP神经网络模糊方程模型的软测量模型5建立完成。
步骤6:熔融指数软测量值显示仪6显示BP粒子群优化的工业熔融指数软测量模型5的输出,完成对工业聚丙烯生产熔融指数软测量的显示。
Claims (2)
1.一种BP粒子群优化的工业熔融指数软测量仪表,包括用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库以及熔融指数软测量值显示仪,所述现场智能仪表、控制站与DCS数据库连接,其特征在于:所述软测量仪表还包括BP粒子群优化的工业熔融指数软测量模型,所述DCS数据库与所述BP粒子群优化的工业熔融指数软测量模型的输入端连接,所述BP粒子群优化的工业熔融指数软测量模型的输出端与熔融指数软测量值显示仪连接,所述BP粒子群优化的工业熔融指数软测量模型包括:
数据预处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:
计算均值:
计算方差:
标准化:
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本;σx表示训练样本的标准差,σ2 x表示训练样本的方差;
模糊方程模块,对从数据预处理模块传过来的标准化后的训练样本X,进行模糊化;设模糊方程***中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,取作2,||·||为范数表达式;
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xi,μik)=[1 func(μik) Xi] (5)
其中func(μik)为隶属度值μik的变形函数,取或exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵;
BP神经网络作为模糊方程***的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为sl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程***的输出:
其中,为第k个BP神经网络模糊方程输出层的预测输出;
粒子群算法优化模块,用于采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体实现步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest;
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程***的误差函数,表示为:
式中,是模糊方程***的预测输出,Oi为模糊方程***的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为整个粒子群的全局最优值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax;
所述BP粒子群优化的工业熔融指数软测量模型还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新模糊方程模型。
2.一种用如权利要求1所述的BP粒子群优化的工业熔融指数软测量仪表实现的软测量方法,其特征在于:所述软测量方法具体实现步骤如下:
1)、对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
2)、将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化,使得其均值为0,方差为1;该处理采用以下算式过程来完成:
2.1)计算均值:
2.2)计算方差:
2.3)标准化:
其中,TXi为第i个训练样本,N为训练样本数,为训练样本的均值,X为标准化后的训练样本;σx表示训练样本的标准差,σ2 x表示训练样本的方差;
3)、对从数据预处理模块传过来的训练样本,进行模糊化;设模糊方程***中有c*个模糊群,模糊群k、j的中心分别为vk、vj,则标准化后的第i个训练样本Xi对于模糊群k的隶属度μik为:
式中,n为模糊分类过程中需要的分块矩阵指数,取作2,||·||为范数表达式;
使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:
Φik(Xi,μik)=[1 func(μik) Xi] (5)
其中func(μik)为隶属度值μik的变形函数,取或exp(μik)等,Φik(Xi,μik)表示第i个输入变量Xi及其模糊群k的隶属度μik所对应的新的输入矩阵;
BP神经网络作为模糊方程***的局部方程,设第k个BP神经网络模糊方程输出层的预测输出为输入为net,与此层相邻的隐含层中任一神经元l的输出为sl,则有:
式中,wlk是隐层神经元l与输出层神经元之间的连接权,f(·)为神经元的输出函数,通常取为Sigmoid函数,表示为:
式中,h为输出层神经元的阈值,θ0为陡度参数,用以调节Sigmoid函数的陡度;最后,由反模糊方法中的重心法得到最后的模糊方程***的输出:
其中,为第k个BP神经网络模糊方程输出层的预测输出;
4)、采用粒子群算法对模糊方程中BP神经网络局部方程的wlk进行优化,具体实现步骤如下:
①确定粒子群的优化参数为BP神经网络局部方程的wlk、粒子群个体数目popsize、最大循环寻优次数itermax、第p个粒子的初始位置rp、初始速度vp、局部最优值Lbestp以及整个粒子群的全局最优值Gbest;
②设定优化目标函数,将其转换为适应度,对每个局部模糊方程进行评价;通过相应的误差函数计算适应度函数,并认为误差大的粒子适应度小,粒子p的适应度函数表示为:
fp=1/(Ep+1) (10)
式中,Ep是模糊方程***的误差函数,表示为:
式中,是模糊方程***的预测输出,Oi为模糊方程***的目标输出;
③按照如下公式,循环更新每个粒子的速度和位置,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(12)
rp(iter+1)=rp(iter)+vp(iter+1) (13)
式中,vp表示更新粒子p的速度,rp表示更新粒子p的位置,Lbestp表示更新粒子p的个体最优值,Gbest即为整个粒子群的全局最优值,iter表示循环次数,ω是粒子群算法中的惯性权重,m1、m2是对应的加速系数,a1、a2是[0,1]之间的随机数;
④对于粒子p,如果新的适应度大于原来的个体最优值,更新粒子的个体最优值:
Lbestp=fp (14)
⑤如果粒子p的个体最优值Lbestp大于原来的粒子群全局最优值Gbest,更新原来的粒子群全局最优值Gbest:
Gbest=Lbestp (15)
⑥判断是否满足性能要求,若是,结束寻优,得到一组优化的模糊方程的局部方程参数;否则返回步骤③,继续迭代寻优,直至达到最大迭代次数itermax;
所述软测量方法还包括以下步骤:5)、定期将离线化验数据输入到训练集中,更新模糊方程模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310436856.9A CN103675012B (zh) | 2013-09-22 | 2013-09-22 | Bp粒子群优化的工业熔融指数软测量仪表和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310436856.9A CN103675012B (zh) | 2013-09-22 | 2013-09-22 | Bp粒子群优化的工业熔融指数软测量仪表和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675012A CN103675012A (zh) | 2014-03-26 |
CN103675012B true CN103675012B (zh) | 2015-08-19 |
Family
ID=50313217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310436856.9A Expired - Fee Related CN103675012B (zh) | 2013-09-22 | 2013-09-22 | Bp粒子群优化的工业熔融指数软测量仪表和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675012B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104408518B (zh) * | 2014-11-12 | 2015-08-26 | 山东地纬数码科技有限公司 | 基于粒子群优化算法的神经网络学习优化方法 |
CN111382933B (zh) * | 2020-03-04 | 2023-05-09 | 海南金盘智能科技股份有限公司 | 一种变压器排程方案的生成方法及*** |
CN113425273A (zh) * | 2021-06-21 | 2021-09-24 | 合肥工业大学 | 一种建模个体脉搏波特征的血压测量*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101315557A (zh) * | 2008-06-25 | 2008-12-03 | 浙江大学 | 基于遗传算法优化bp神经网络的丙烯聚合生产过程最优软测量仪表及方法 |
CN101382801A (zh) * | 2008-06-25 | 2009-03-11 | 浙江大学 | 基于ega优化的丙烯聚合生产过程最优软测量仪表及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102008036968A1 (de) * | 2008-08-08 | 2010-02-11 | Endress + Hauser Gmbh + Co. Kg | Diagnoseverfahren eines Prozessautomatisierungssystem |
-
2013
- 2013-09-22 CN CN201310436856.9A patent/CN103675012B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101315557A (zh) * | 2008-06-25 | 2008-12-03 | 浙江大学 | 基于遗传算法优化bp神经网络的丙烯聚合生产过程最优软测量仪表及方法 |
CN101382801A (zh) * | 2008-06-25 | 2009-03-11 | 浙江大学 | 基于ega优化的丙烯聚合生产过程最优软测量仪表及方法 |
Non-Patent Citations (3)
Title |
---|
Comparison of particle swarm optimization and backpropagation as training algorithms for neural networks;Gudise V G et al.;《Proceedings of the 2003 IEEE》;20031231;全文 * |
基于改进混沌粒子群优化算法的聚丙烯生产调优;夏陆岳等;《计算机与应用化学》;20120728;第29卷(第7期);全文 * |
应用混沌粒子群优化训练的BP神经网络预报高炉铁水含硅量;徐生林等;《冶金分析》;20101231;第30卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103675012A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103674778B (zh) | Rbf粒子群优化的工业熔融指数软测量仪表和方法 | |
CN101315557B (zh) | 基于遗传算法优化bp神经网络的丙烯聚合生产过程最优软测量仪表及方法 | |
CN103675006B (zh) | 最小二乘的工业熔融指数软测量仪表和方法 | |
CN103675011B (zh) | 最优支持向量机的工业熔融指数软测量仪表和方法 | |
CN103472865B (zh) | 智能最小二乘的农药废液焚烧炉炉温最佳化***及方法 | |
CN103472866B (zh) | 智能模糊***的农药废液焚烧炉炉温最佳化***及方法 | |
CN103675012B (zh) | Bp粒子群优化的工业熔融指数软测量仪表和方法 | |
CN103675005B (zh) | 最优模糊网络的工业熔融指数软测量仪表及方法 | |
CN103675010B (zh) | 支持向量机的工业熔融指数软测量仪表及方法 | |
CN103675009B (zh) | 模糊方程的工业熔融指数软测量仪表及方法 | |
CN103839103A (zh) | 丙烯聚合生产过程bp最优预报***和方法 | |
CN103630568B (zh) | Bp网络的工业熔融指数软测量仪表及方法 | |
CN103675007B (zh) | Rbf网络的工业熔融指数软测量仪表及方法 | |
CN103838205B (zh) | Bp全局最优丙烯聚合生产过程最优软测量仪表和方法 | |
CN103838209B (zh) | 丙烯聚合生产过程自适应最优预报***和方法 | |
CN109856971A (zh) | 基于群智寻优的丙烯聚合生产过程最优在线预报*** | |
CN103675008B (zh) | 加权模糊的工业熔融指数软测量仪表及方法 | |
CN103838958A (zh) | 模糊智能最优丙烯聚合生产过程最优软测量仪表和方法 | |
CN103838957A (zh) | 丙烯聚合生产过程径向基最优软测量仪表和方法 | |
CN103472729B (zh) | 群智的农药废液焚烧炉有害物排放达标控制***及方法 | |
CN103488206B (zh) | 智能径向基的农药生产废液焚烧炉炉温最佳化***及方法 | |
CN103823963A (zh) | 一种鲁棒聚丙烯熔融指数软测量***和方法 | |
CN103823369A (zh) | 基于bp多模网络的丙烯聚合生产过程预报***和方法 | |
CN103488209B (zh) | 智能支持向量机的农药废液焚烧炉炉温最佳化***及方法 | |
CN103488207B (zh) | 模糊***的农药生产废液焚烧炉炉温最佳化***及方法 |
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: 20150819 Termination date: 20180922 |
|
CF01 | Termination of patent right due to non-payment of annual fee |