CN102626557B - 基于ga-bp算法的分子蒸馏工艺参数优化方法 - Google Patents

基于ga-bp算法的分子蒸馏工艺参数优化方法 Download PDF

Info

Publication number
CN102626557B
CN102626557B CN 201210107643 CN201210107643A CN102626557B CN 102626557 B CN102626557 B CN 102626557B CN 201210107643 CN201210107643 CN 201210107643 CN 201210107643 A CN201210107643 A CN 201210107643A CN 102626557 B CN102626557 B CN 102626557B
Authority
CN
China
Prior art keywords
value
neural network
yield
output
population
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
CN 201210107643
Other languages
English (en)
Other versions
CN102626557A (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.)
Changchun University of Technology
Original Assignee
Changchun University of Technology
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 Changchun University of Technology filed Critical Changchun University of Technology
Priority to CN 201210107643 priority Critical patent/CN102626557B/zh
Publication of CN102626557A publication Critical patent/CN102626557A/zh
Application granted granted Critical
Publication of CN102626557B publication Critical patent/CN102626557B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明基于GA-BP算法的分子蒸馏工艺参数优化方法涉及分子蒸馏优化领域,该方法是,首先利用GA算法优化BP神经网络的权值和阈值,再利用BP神经网络实现分子蒸馏工艺中由蒸发面温度、蒸发器内真空度、进料速度三个工艺控制参数预测出馏出物的纯度和得率,最后,利用GA算法实现分子蒸馏工艺中由馏出物的期望纯度和得率值得到优化的工艺控制参数输入值的过程。本发明的方法能够及时地掌握馏出物的纯度和得率,并作出精确地预测,可以及时采取调节措施,使影响馏出物产品质量的关键因素得到合理的控制,从而保证产品的纯度和得率在规定范围内,有效地提高企业的生产效率,节约能源,降低企业的生产成本。

Description

基于GA-BP算法的分子蒸馏工艺参数优化方法
技术领域
本发明涉及分子蒸馏工艺条件优化领域,具体涉及一种基于GA(Genetic Algorithm,遗传算法)和BP(Back Propagation,误差反向传播)网络的分子蒸馏工艺参数优化方法。
背景技术
目前,药厂、保健品生产等厂家都非常关注精油的提取技术,以提取五味子精油为例,在药厂中主要采用浸泡法和水蒸气蒸馏法,这两种方法虽然成本较低,但成分为混合物,精油成分不高,纯度很低,且提取时间较长,浸泡的时间多达几小时到十几小时甚至更长;近些年发展了超临界萃取、分子蒸馏等适合于生产中应用的提取方法,来实现精油提取,其中超临界提取仍为粗提的方法,提取的精油中混有水份、溶剂、杂质等,产品的纯度和得率较低,企业生产效率和经济效益仍然得不到显著提高,所以分子蒸馏成为精制精油提取的最佳方法,分子蒸馏法相对其它提取方法虽然效率高,但是,该方法馏出物的纯度和得率受蒸发面温度、蒸发器内真空度、进料速率三个工艺参数的影响最大,能否准确确定这些工艺参数值,将直接影响产品的质量和生产效率。
分子蒸馏生产过程中工艺参数大多采用经验参数,缺少科学依据,且生产过程复杂,无法建立准确的数学模型,目前较多地用显式的统计分析方法来解决,但是统计分析模型缺乏自适应性,模型更新困难,参数确定不准确;也有采用正交试验进行工艺参数优化的方法,如当真空度不变,通过多次试验得到最佳温度,在此温度下再通过多次蒸馏试验摸索最佳真空度,而分子蒸馏蒸发器的温度与真空度具有耦合关系,当真空度改变,物质的沸点发生变化,即真空度变化,最佳蒸馏温度随之变化,所以需要应用更科学有效的实验参数优化方法来确定最佳工艺参数。本发明采用GA-BP算法实现分子蒸馏工艺参数优化,根据产品指标要求(得率、纯度),运用该算法获得最佳工艺参数值。在实际蒸馏过程中,由于蒸发面温度、蒸发器内真空度、进料速度是影响蒸馏产品指标最重要的工艺参数,也是蒸馏过程需要精确控制的三个量,所以本发明以这三个量作为优化对象,设计了GA-BP蒸馏工艺参数优化方法,对于指导生产具有很好的实际意义。
发明内容
为了解决传统方法确定的分子蒸馏工艺参数不准确,馏出物的纯度和得率低的问题,本发明提供一种基于GA-BP算法的分子蒸馏工艺参数优化方法。
该方法首先建立以蒸发面温度、蒸发器内真空度、进料速度作为输入量的BP神经网络预测模型,通过遗传算法优化该预测模型的权值和阈值,然后对该预测模型进行训练,得到稳定的预测模型;对训练好的BP神经网络预测模型在一定范围内随机给定多组输入向量(蒸发面温度、蒸发器内真空度、进料速度),利用BP神经网络预测模型预测出相应的输出量(纯度和得率),通过个体适应度函数,计算个体适应度值,再利用遗传算子调整输入向量,产生具有更好适应性的新的种群(输入向量),重复以上优化过程,直到适应度函数值收敛,与该值对应的输入向量即为一组最优工艺参数。
本发明基于GA-BP算法的分子蒸馏工艺参数优化方法包括以下步骤:
步骤一:将分子蒸馏系统在同一时刻下的一组蒸发面温度、蒸发器内真空度、进料速度三个对蒸馏过程影响最大的工艺控制参数的值作为输入量Pk=(a1,a2,a3),其中,a1、a2、a3分别是蒸发面温度、蒸发器内真空度、进料速度,k是n组实验样本中的任意一组,k=1,2,…,n,将此时分子蒸馏系统的馏出物的纯度值和得率值作为输出量Ck=(c1,c2),其中,c1、c2分别是纯度值和得率值,确立一个在同一时刻下输入量Pk与输出量Ck具有映射对应关系的实验数据组;在分子蒸馏过程稳态进行的情况下,在不同时刻提取并记录n个具有上述映射关系的实验数据组,将这些实验数据组作为实验样本,存入实验样本数据库中;
步骤二:根据步骤一所述实验样本中输入量与输出量的映射对应关系,建立三层拓扑结构的BP神经网络,设定BP神经网络的输入层神经元为3个,输出层神经元为2个,隐含层神经元为7个;其中输入层和隐含层的激活函数都选取Log-sigmoid型函数,输出层的激活函数选取Pureline型函数;Log-sigmoid型函数定义为
Figure GDA0000368290600000021
其中sj是指第j个隐含层节点的输入,即
Figure GDA0000368290600000022
其中j=1,2…7;wij为第i个输入节点与第j个输出节点的连接权值,θj为隐含层第j节点的阈值,ai为蒸发面温度、蒸发器内真空度、进料速度;Pureline型函数是线性函数,自变量为隐含层节点的输出,因变量是输出层节点的输出;
步骤三:以步骤二所述的BP神经网络作为黑箱模型,建立以步骤一所述实验样本中的输入量作为BP神经网络的输入变量、以对应馏出物纯度和得率的预测值作为输出变量的BP神经网络预测模型;
步骤四:对步骤三所述的BP神经网络预测模型中的输入层与隐含层之间的连接权值wij、阈值θij,隐含层与输出层之间的连接权值vjq、阈值γq进行随机初始化,使BP神经网络预测模型具备最基本的预测条件,其中q为输出层节点数,q=1,2;
步骤五:将从步骤一所述的实验样本数据库中调取的各组工艺控制参数的值,作为步骤四所述的具备最基本预测条件的BP神经网络预测模型的输入变量Pk,从而得到与每组输入变量Pk=(a1,a2,a3)映射对应的BP神经网络预测模型输出的纯度值和得率值的预测输出量Yk的数据组;
步骤六:在由步骤一所述实验样本确定的阈值和权值范围内,根据GA遗传算法的原理优化BP神经网络预测模型的初始阈值和权值,完成BP神经网络预测模型的训练学习过程;
步骤6-1:随机初始化多组BP神经网络预测模型的阈值θm和权值wz作为GA遗传算法的初始种群,种群大小根据经验预先设定,并对该初始种群进行实数编码,设定最大进化代数为100代;
步骤6-2:构造用于优化BP神经网络预测模型初始阈值和权值的GA遗传算法的初始种群中第k个样本的的个体适应度函数
Figure GDA0000368290600000031
根据个体适应度函数公式计算出初始种群中每个个体的适应度f的数值;
式中,
Figure GDA0000368290600000032
为根据实际生产需要预先设定的分子蒸馏后的纯度值和得率值的期望输出量;
Figure GDA0000368290600000033
为步骤五所述的纯度值和得率值的预测输出量,q为输出节点数;
步骤6-3:依据轮盘赌的选择法,从步骤6-2所述初始种群的多个个体中选择对应个体适应度值f数值相对较大的多个待优化个体,进行交叉、变异操作,产生的新个体作为子一代阈值θ'm和权值w'z种群的个体;
步骤6-4:用步骤6-3所述子一代阈值θ'm和权值w'z的种群的新个体替换步骤6-1所述的初始种群的个体,重复步骤6-1至步骤6-3的GA遗传算法对子一代种群个体的优化过程,直到步骤6-2所述的个体适应度函数基本稳定不变或者达到设定的进化代数时,结束优化过程,并得到初步优化完毕的最后一代种群的个体,即得到了初步优化完毕的BP神经网络初始阈值θ"′m和权值w"′z
步骤6-5:从步骤一所述的实验样本数据库中调取由分子蒸馏实验获得的实验样本,此实验样本中的输入量Pk与输出量Yk具有真实的映射对应关系;以相同的实验样本中的输入量Pk代入到BP神经网络预测模型中,并将步骤6-4所述初步优化完毕的阈值θ"′m和权值w"′z作为BP神经网络预测模型的新初始阈值θ和权值w,此时由BP神经网络预测模型得到预测输出量Y'k,此预测输出量Y'k与分子蒸馏期望输出量Tk之间存在误差;
步骤6-6:建立BP神经网络的误差反向传播模型,构造一个由分子蒸馏期望输出量Tk与BP神经网络预测模型的预测输出量Y'k的差的平方和的目标函数,该函数第k个样本的平均误差为
Figure GDA0000368290600000041
Figure GDA0000368290600000042
为根据实际生产需要预先设定的分子蒸馏后的纯度值和得率值的期望输出量;为步骤6-5所述的纯度值和得率值的预测输出量,q为输出节点数,使这个目标函数最小化的过程就是BP神经网络误差反向传播的过程;
步骤6-7:用梯度下降法对步骤6-6所述的BP神经网络误差反向传播的目标函数进行运算,使得步骤6-5中所述BP神经网络预测模型的初始阈值θ和权值w进行进一步优化;按照梯度下降法,伴随迭代次数的增加,目标函数的误差将逐步减小,直至误差满足预先设定的精度要求,结束初始阈值θ和权值w的优化过程;
步骤6-8:经过步骤6-7的反复计算和缩小误差,当误差最终满足预先设定的精度要求时,将会得到一组最优的权值w0和阈值θ0,将此时的权值w0和阈值θ0作为BP神经网络预测模型最终的权值和阈值,不再需要调整和改变,BP神经网络预测模型的训练学习过程至此结束;
步骤七:利用遗传算法实现通过给定多个分子蒸馏系统期望的馏出物的纯度值和对应此时的得率值,得到对应的分子蒸馏系统的工艺控制参数的输入值即蒸发面温度值、蒸发器内真空度值、进料速度值,该遗传算法的实现过程如下:
步骤7-1:给定多个分子蒸馏系统期望的馏出物的纯度值和得率值Tj,从步骤一所述的实验样本数据库中选取与期望的纯度值和得率值Tj相对接近的分子蒸馏实验真实输出量Cj所对应的多个工艺控制参数输入量Pj作为遗传算法的初始种群;
步骤7-2:对步骤7-1所述的多个工艺控制参数输入量Pj的初始种群进行实数编码,设定最大进化代数为100代;
步骤7-3:构造遗传算法的目标函数:设经过BP神经网络预测模型得到的纯度和得率的值分别为y1、y2,期望的纯度和得率的值分别为T1和T2,纯度和得率的权重分别为W1和W2,纯度的最大、最小值分别为T1max和T1min,得率最大、最小值分别为T2max和T2min,则纯度和得率两个期望目标组合加权得到的目标函数f(X)为:
f ( X ) = W 1 * | y 1 - T 1 | | T 1 max - T 1 min | + W 2 * | y 2 - T 2 | | T 2 max - T 2 min | T 1 max > T 1 min , T 2 max > T 2 min , W 1 + W 2 = 1 , 式中目标函数值越小,则表明BP神经网络预测模型预测的纯度值y1、得率值y2越接近期望的纯度值T1、得率值T2
步骤7-4:选取步骤7-3所述的目标函数f(X)的倒数作为遗传算法的初始种群的个体适应度函数f(Y),即f(Y)=1/f(X),则个体适应度函数f(Y)的值越大,表明个体适应能力越强;
步骤7-5:针对步骤7-1所述的遗传算法的初始种群,根据步骤7-4中的个体适应度函数f(Y),计算该初始种群中每个工艺控制参数输入量Pj个体的适应度值;
步骤7-6:依据轮盘赌的选择方法,选出步骤7-5中个体适应度值相对较大的多个工艺控制参数输入量Pj个体,进行交叉、变异操作,产生的新工艺控制参数输入量Pj'作为子一代的初始种群个体;
步骤7-7:将经步骤7-6产生的新工艺控制参数输入量Pj'初始种群替换步骤7-1所述的工艺控制参数输入量Pj初始种群;
步骤7-8:将步骤7-6产生的每个子一代新工艺控制参数输入量Pj'种群的个体都代入步骤6-8所述的训练好的BP神经网络预测模型中,得出新一组的预测结果y1'和y2';用该新一组的预测结果替代步骤7-3中纯度和得率的预测值y1、y2
步骤7-9:将步骤7-8中得出的新一组的预测结果y1'和y2'代入步骤7-4的个体适应度函数f(Y),计算新种群个体对应下的新一组的工艺控制参数输入量Pj种群个体的适应度函数值;
步骤7-10:重复步骤7-4至步骤7-9所述更新下一代种群个体、计算对应新种群下的新的预测值和计算新的个体适应度函数f(Y)值的过程;在该重复过程中个体适应度函数f(Y)值将逐步增大;直至个体适应度函数f(Y)值基本稳定不再明显变化或达到预设的进化代数时,终止遗传算法的优化过程;此时稳定的个体适应度函数值对应的工艺控制参数输入量Pj就是对应期望纯度值和得率值所要输入的蒸发面温度值、蒸发器内真空度值、进料速度值。
本发明的有益效果是:该方法能够预测馏出物的纯度和得率,并根据预测值使影响馏出物产品质量的关键因素得到合理的控制,从而保证产品的纯度和得率在规定范围内,有效地提高企业的生产效率,节约能源,降低企业的生产成本。
具体实施方式
下面结合实例对本发明的技术方案做进一步详细说明。
本发明以分子蒸馏法提取五味子精油的生产实验为例,实验以干燥果实五味子为原料,首先,通过超临界CO2萃取装置萃取得到五味子粗油,然后以该五味子粗油为原料,通过长春工业大学的DCH-300自动控制三级刮膜式分子蒸馏装置进行分子蒸馏,提取五味子精油的实验。
人工神经网络的学习是对数据进行归纳学习的过程,通过对数据样本进行反复学习,来不断修正各个神经元之间相互连接的权值和阈值,最终实现将人工神经网络的权值和阈值收敛于一个稳定的范围。以三层BP神经网络为例,三层BP神经网络算法的核心是为了实现或逼近输入信号与输出信号之间的映射关系,BP网络一方面进行信号的正向传播,一方面进行误差的反向传播,并根据误差修正网络的权值和阈值,这一过程是反复进行的,通过不断的修正权值和阈值,使网络的整体误差达到允许范围或者训练次数限制为止。
遗传算法的基本步骤为:初始群体、编码、计算适应度函数、交叉变异操作和解码。
本发明针对BP网络在实际应用中表现出来的不足,采用三层BP神经网络,利用遗传算法优化BP神经网络各层神经元的连接权值和阈值以及BP网络的构造,从大量数据中选取训练样本并设计,通过反复实验确定隐含层节点数。遗传算法通过这些方法使模型得到不断改进和优化,基于以上遗传算法的BP神经网络预测模型被建立。结合生产实验,充分考虑分子蒸馏生产过程的实际情况,实时采集数据、周期变化等特点,对BP网络的样本进行设计和参数选取,再反复通过网络训练和检验得到较为合理的预测结果。
本发明的方法是以实验采集的真实数据为样本数据,对样本数据进行统计、分析、整合,并在此基础上建立基于GA-BP算法的多级分子蒸馏系统中五味子精油纯度和得率的预测模型;然后,在预测模型的基础上建立基于遗传算法的参数优化模型,提出了分子蒸馏系统中优化的对象,根据优化目标值确定出优化量值(蒸发面温度、蒸发器内真空度、进料速度),为企业工业化生产提供工艺优化策略。
本发明基于GA-BP算法的分子蒸馏工艺参数优化方法包括以下步骤:
步骤一:以25组真实实验数据建立实验样本,即,将分子蒸馏系统的工艺控制参数蒸发面温度值、蒸发器内真空度值、进料速度值作为输入量,将分子蒸馏系统得到的馏出物的纯度值和得率值作为输出量,建立具有映射对应关系的25组实验数据组作为实验样本;
步骤二:建立三层拓扑结构的BP神经网络,设定BP神经网络的输入层神经元为3个,分别对应蒸发面温度值、蒸发器内真空度值和进料速度值,输出层神经元为2个,分别对应馏出物的纯度值和得率值,隐含层神经元为7个;
步骤三:以步骤二所述的BP神经网络作为黑箱模型,建立以步骤一所述实验样本中的输入量作为输入变量、以对应馏出物纯度和得率的预测值作为输出变量的BP神经网络预测模型;
步骤四:对BP神经网络预测模型中的输入层与隐含层之间的连接权值、阈值,隐含层与输出层之间的连接权值、阈值进行随机初始化,使BP神经网络预测模型具备最基本的预测条件;
步骤五:将实验样本中的25组工艺控制参数的值作为BP神经网络预测模型的输入变量,得到与每组输入变量映射对应的纯度和得率的预测输出量的数据组;
步骤六:在由25组实验样本确定的阈值和权值范围内,首先利用GA遗传算法得到初步优化完毕的BP神经网络预测模型的初始阈值和权值,然后再通过BP神经网络的误差反向传播模型的目标函数及梯度下降法,优化上述初步优化完毕的BP神经网络预测模型的初始阈值和权值,进而完成BP神经网络预测模型的训练学习过程。其中,将实验样本中的16组数据作为BP神经网络预测模型的训练样本,9组数据作为BP神经网络预测模型的检验样本以验证BP神经网络预测模型已训练好;
按照梯度下降法,伴随迭代次数的增加,BP神经网络的误差反向传播模型的目标函数的误差将逐步减小,直至误差满足预先设定的纯度和得率的精度要求,BP网络的初始阈值和权值的优化过程结束。
以BP神经网络的隐含层权值wij的修正为例,权值修正的过程就是使平均误差E逐渐减小的过程。具体如下:
wij(N+1)=wij(N)-η×g(N),于是
Figure GDA0000368290600000071
将平均误差
Figure GDA0000368290600000072
的表达式代入
Figure GDA0000368290600000073
其中,c表示输出单元个数,令ej=tj-yj,它是关于wij的隐函数,其中yj=f(netj),yj表示第j个神经元的输出,
Figure GDA0000368290600000074
表示第j个神经元的输入,f表示传递函数;需要进行逐层微分得到
Figure GDA0000368290600000081
于是有
Figure GDA0000368290600000082
其中δj=ejf'(netj),η为BP网络的学习效率,为预先设定的常数。
步骤七:以训练好的BP神经网络预测模型为基础,利用遗传算法实现通过给定多个分子蒸馏系统期望的馏出物的纯度值和得率值,得到对应的分子蒸馏系统的工艺控制参数的输入值即蒸发面温度值、蒸发器内真空度值、进料速度值。
在分子蒸馏实验过程中,本发明充分利用了非线性拟合能力较强的神经网络和具有全局寻优特性的遗传算法,将两种方法相结合,使各自的优点充分发挥,能在很大的程度上克服传统方法的缺陷,完成网络的训练、检验和最优评价,为生产过程的决策和控制提供可靠的依据。本发明利用BP神经网络算法建立了工艺控制参数即蒸发面温度值、蒸发器内真空度值、进料速度值对产品的纯度与得率影响的预测模型,基于遗传算法以产品纯度和得率为综合目标函数,将遗传算法引入到分子蒸馏提纯工艺控制参数的优化中,获得了优化的工艺控制参数。本发明的方法与传统的方法相比,算法简单,计算精度高,有很强的自学习能力。

Claims (1)

1.基于GA-BP算法的分子蒸馏工艺参数优化方法,其特征在于,该方法包括以下步骤:
步骤一:将分子蒸馏系统在同一时刻下的一组蒸发面温度、蒸发器内真空度、进料速度三个对蒸馏过程影响最大的工艺控制参数的值作为输入量Pk=(a1,a2,a3),其中,a1、a2、a3分别是蒸发面温度、蒸发器内真空度、进料速度,k是n组实验样本中的任意一组,k=1,2,…,n,将此时分子蒸馏系统的馏出物的纯度值和得率值作为输出量Ck=(c1,c2),其中,c1、c2分别是纯度值和得率值,确立一个在同一时刻下输入量Pk与输出量Ck具有映射对应关系的实验数据组;在分子蒸馏过程稳态进行的情况下,在不同时刻提取并记录n个具有上述映射关系的实验数据组,将这些实验数据组作为实验样本,存入实验样本数据库中;
步骤二:根据步骤一所述实验样本中输入量与输出量的映射对应关系,建立三层拓扑结构的BP神经网络,设定BP神经网络的输入层神经元为3个,输出层神经元为2个,隐含层神经元为7个;其中输入层和隐含层的激活函数都选取Log-sigmoid型函数,输出层的激活函数选取Pureline型函数;Log-sigmoid型函数定义为
Figure FDA0000368290590000011
其中sj是指第j个隐含层节点的输入,即
Figure FDA0000368290590000012
其中j=1,2…7;wij为第i个输入节点与第j个输出节点的连接权值,θj为隐含层第j节点的阈值,ai为蒸发面温度、蒸发器内真空度、进料速度;Pureline型函数是线性函数,自变量为隐含层节点的输出,因变量是输出层节点的输出;
步骤三:以步骤二所述的BP神经网络作为黑箱模型,建立以步骤一所述实验样本中的输入量作为BP神经网络的输入变量、以对应馏出物纯度和得率的预测值作为输出变量的BP神经网络预测模型;
步骤四:对步骤三所述的BP神经网络预测模型中的输入层与隐含层之间的连接权值wij、阈值θij,隐含层与输出层之间的连接权值vjq、阈值γq进行随机初始化,使BP神经网络预测模型具备最基本的预测条件,其中q为输出层节点数,q=1,2;
步骤五:将从步骤一所述的实验样本数据库中调取的各组工艺控制参数的值,作为步骤四所述的具备最基本预测条件的BP神经网络预测模型的输入变量Pk,从而得到与每组输入变量Pk=(a1,a2,a3)映射对应的BP神经网络预测模型输出的纯度值和得率值的预测输出量Yk的数据组;
步骤六:在由步骤一所述实验样本确定的阈值和权值范围内,根据GA遗传算法的原理优化BP神经网络预测模型的初始阈值和权值,完成BP神经网络预测模型的训练学习过程;
步骤6-1:随机初始化多组BP神经网络预测模型的阈值θm和权值wz作为GA遗传算法的初始种群,种群大小根据经验预先设定,并对该初始种群进行实数编码,设定最大进化代数为100代;
步骤6-2:构造用于优化BP神经网络预测模型初始阈值和权值的GA遗传算法的初始种群中第k个样本的个体适应度函数
Figure FDA0000368290590000021
根据个体适应度函数公式计算出初始种群中每个个体的适应度f的数值;
式中,
Figure FDA0000368290590000022
为根据实际生产需要预先设定的分子蒸馏后的纯度值和得率值的期望输出量;为步骤五所述的纯度值和得率值的预测输出量,q为输出节点数;
步骤6-3:依据轮盘赌的选择法,从步骤6-2所述初始种群的多个个体中选择对应个体适应度值f数值相对较大的多个待优化个体,进行交叉、变异操作,产生的新个体作为子一代阈值θ'm和权值w'z种群的个体;
步骤6-4:用步骤6-3所述子一代阈值θ'm和权值w'z的种群新个体替换步骤6-1所述的初始种群的个体,重复步骤6-1至步骤6-3的GA遗传算法对子一代种群个体的优化过程,直到步骤6-2所述的个体适应度函数基本稳定不变或者达到设定的进化代数时,结束优化过程,并得到初步优化完毕的最后一代种群的个体,即得到了初步优化完毕的BP神经网络初始阈值θ"′m和权值w"′z
步骤6-5:从步骤一所述的实验样本数据库中调取由分子蒸馏实验获得的实验样本,此实验样本中的输入量Pk与输出量Yk具有真实的映射对应关系;以相同的实验样本中的输入量Pk代入到BP神经网络预测模型中,并将步骤6-4所述初步优化完毕的阈值θ"′m和权值w"′z作为BP神经网络预测模型的新初始阈值θ和权值w,此时由BP神经网络预测模型得到预测输出量Y'k,此预测输出量Y'k与分子蒸馏期望输出量Tk之间存在误差;
步骤6-6:建立BP神经网络的误差反向传播模型,构造一个由分子蒸馏期望输出量Tk与BP神经网络预测模型的预测输出量Y'k的差的平方和的目标函数,该函数第k个样本的平均误差为
Figure FDA0000368290590000024
Figure FDA0000368290590000025
为根据实际生产需要预先设定的分子蒸馏后的纯度值和得率值的期望输出量;
Figure FDA0000368290590000026
为步骤6-5所述的纯度值和得率值的预测输出量,q为输出节点数,使这个目标函数最小化的过程就是BP神经网络误差反向传播的过程;
步骤6-7:用梯度下降法对步骤6-6所述的BP神经网络误差反向传播的目标函数进行运算,使得步骤6-5中所述BP神经网络预测模型的初始阈值θ和权值w进行进一步优化;按照梯度下降法,伴随迭代次数的增加,目标函数的误差将逐步减小,直至误差满足预先设定的精度要求,结束初始阈值θ和权值w的优化过程;
步骤6-8:经过步骤6-7的反复计算和缩小误差,当误差最终满足预先设定的精度要求时,将会得到一组最优的权值w0和阈值θ0,将此时的权值w0和阈值θ0作为BP神经网络预测模型最终的权值和阈值,不再需要调整和改变,BP神经网络预测模型的训练学习过程至此结束;
步骤七:利用遗传算法实现通过给定多个分子蒸馏系统期望的馏出物的纯度值和对应此时的得率值,得到对应的分子蒸馏系统的工艺控制参数的输入值即蒸发面温度值、蒸发器内真空度值、进料速度值,该遗传算法的实现过程如下:
步骤7-1:给定多个分子蒸馏系统期望的馏出物的纯度值和得率值Tj,从步骤一所述的实验样本数据库中选取与期望的纯度值和得率值Tj相对接近的分子蒸馏实验真实输出量Cj所对应的多个工艺控制参数输入量Pj作为遗传算法的初始种群;
步骤7-2:对步骤7-1所述的多个工艺控制参数输入量Pj的初始种群进行实数编码,设定最大进化代数为100代;
步骤7-3:构造遗传算法的目标函数:设经过BP神经网络预测模型得到的纯度和得率的值分别为y1、y2,期望的纯度和得率的值分别为T1和T2,纯度和得率的权重分别为W1和W2,纯度的最大、最小值分别为T1max和T1min,得率最大、最小值分别为T2max和T2min,则纯度和得率两个期望目标组合加权得到的目标函数f(X)为:
f ( X ) = W 1 * | y 1 - T 1 | | T 1 max - T 1 min | + W 2 * | y 2 - T 2 | | T 2 max - T 2 min | T 1 max > T 1 min , T 2 max > T 2 min , W 1 + W 2 = 1 , 式中目标函数值越小,则表明BP神经网络预测模型预测的纯度值y1、得率值y2越接近期望的纯度值T1、得率值T2
步骤7-4:选取步骤7-3所述的目标函数f(X)的倒数作为遗传算法的初始种群的个体适应度函数f(Y),即f(Y)=1/f(X),则个体适应度函数f(Y)的值越大,表明个体适应能力越强;
步骤7-5:针对步骤7-1所述的遗传算法的初始种群,根据步骤7-4中的个体适应度函数f(Y),计算该初始种群中每个工艺控制参数输入量Pj个体的适应度值;
步骤7-6:依据轮盘赌的选择方法,选出步骤7-5中个体适应度值相对较大的多个工艺控制参数输入量Pj个体,进行交叉、变异操作,产生的新工艺控制参数输入量Pj'个体,作为子一代的初始种群个体;
步骤7-7:将经步骤7-6产生的新工艺控制参数输入量Pj′初始种群替换步骤7-1所述的工艺控制参数输入量Pj初始种群;
步骤7-8:将步骤7-6产生的每个子一代新工艺控制参数输入量Pj′种群的个体都代入步骤6-8所述的训练好的BP神经网络预测模型中,得出新一组的预测结果y1'和y2';用该新一组的预测结果替代步骤7-3中纯度和得率的预测值y1、y2
步骤7-9:将步骤7-8中得出的新一组的预测结果y1'和y2'代入步骤7-4的个体适应度函数f(Y),计算新种群个体对应下的新一组的工艺控制参数输入量Pj种群个体的适应度函数值;
步骤7-10:重复步骤7-4至步骤7-9所述更新下一代种群个体、计算对应新种群下的新的预测值和计算新的个体适应度函数f(Y)值的过程;在该重复过程中个体适应度函数f(Y)值将逐步增大;直至个体适应度函数f(Y)值基本稳定不再明显变化或达到预设的进化代数时,终止遗传算法的优化过程;此时稳定的个体适应度函数值对应的工艺控制参数输入量Pj就是对应期望纯度值和得率值所要输入的蒸发面温度值、蒸发器内真空度值、进料速度值。
CN 201210107643 2012-04-13 2012-04-13 基于ga-bp算法的分子蒸馏工艺参数优化方法 Expired - Fee Related CN102626557B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210107643 CN102626557B (zh) 2012-04-13 2012-04-13 基于ga-bp算法的分子蒸馏工艺参数优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210107643 CN102626557B (zh) 2012-04-13 2012-04-13 基于ga-bp算法的分子蒸馏工艺参数优化方法

Publications (2)

Publication Number Publication Date
CN102626557A CN102626557A (zh) 2012-08-08
CN102626557B true CN102626557B (zh) 2013-10-23

Family

ID=46585060

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210107643 Expired - Fee Related CN102626557B (zh) 2012-04-13 2012-04-13 基于ga-bp算法的分子蒸馏工艺参数优化方法

Country Status (1)

Country Link
CN (1) CN102626557B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103500280A (zh) * 2013-09-30 2014-01-08 沈阳化工大学 谷氨酸发酵过程菌体浓度在线软测量方法
CN104298112B (zh) * 2014-09-28 2017-02-15 长春工业大学 分子蒸馏过程的多变量系统辨识方法
CN104778338A (zh) * 2015-05-03 2015-07-15 长春工业大学 一种低能耗分子蒸馏过程给定值的优化方法
CN106842948B (zh) * 2017-03-01 2019-08-06 长春工业大学 基于bp网络的hdp分子蒸馏系统的最优控制方法
CN107977707B (zh) * 2017-11-23 2020-11-06 厦门美图之家科技有限公司 一种对抗蒸馏神经网络模型的方法及计算设备
CN108073075A (zh) * 2017-12-21 2018-05-25 苏州大学 基于ga优化bp神经网络的硅微加速度计温度补偿方法、系统
CN109100975A (zh) * 2018-09-03 2018-12-28 深圳市智物联网络有限公司 一种参数优化方法及系统
CN110348037B (zh) * 2019-04-22 2020-08-25 武汉理工大学 汽车尾气热电转化装置电气拓扑结构的优化方法
CN113433906B (zh) * 2021-06-24 2022-06-10 之江实验室 一种蒸馏装置产物预测及蒸馏操作参数优化的方法
CN114493057B (zh) * 2022-04-18 2022-07-22 希望知舟技术(深圳)有限公司 基于异常工况的生产工艺参数推荐方法及相关设备
CN115231525B (zh) * 2022-07-29 2022-12-23 福建德尔科技股份有限公司 电子级三氟化氯的智能分离纯化系统
CN116661402B (zh) * 2023-07-27 2023-10-03 北京盛世利科技有限公司 一种化工材料的生产管控方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101333961B (zh) * 2008-08-07 2011-12-28 清华大学 氢气天然气混合燃料发动机的优化方法
US8586365B2 (en) * 2009-12-15 2013-11-19 Exxonmobil Research And Engineering Company Methods for analyzing and optimizing biofuel compositions
CN101833314B (zh) * 2010-03-30 2012-07-25 深圳达实智能股份有限公司 污水处理控制系统及污水处理控制方法

Also Published As

Publication number Publication date
CN102626557A (zh) 2012-08-08

Similar Documents

Publication Publication Date Title
CN102626557B (zh) 基于ga-bp算法的分子蒸馏工艺参数优化方法
CN108053061B (zh) 一种基于改进卷积神经网络的太阳能辐照度预测方法
CN104765346B (zh) 一种炼油过程全流程建模方法
KR102139358B1 (ko) 머신러닝 기반 플랫폼을 이용한 공정제어방법, 그를 수행하기 위한 컴퓨터 프로그램 매체 및 공정제어장치
CN106354017B (zh) 一种稀土萃取分离过程组分含量区间控制方法
CN110210495B (zh) 基于并行LSTM自编码器动态特征提取的XGBoost软测量建模方法
CN103927412A (zh) 基于高斯混合模型的即时学习脱丁烷塔软测量建模方法
CN103177088A (zh) 一种生物医学空缺数据弥补方法
CN108667069B (zh) 一种基于偏最小二乘法回归的短期风电功率预测方法
CN108021773B (zh) 基于dss数据库的分布式水文模型多场次洪水参数率定方法
CN108647373A (zh) 一种基于xgboost模型的工业过程软测量方法
CN102880905B (zh) 一种常顶油干点在线软测量方法
CN105512745A (zh) 基于粒子群—bp神经网络的风电功率区间预测方法
CN109657411B (zh) 一种基于数据驱动的溶剂脱沥青装置建模及优化方法
CN109146186A (zh) 一种基于双重分解的短期风电功率预测方法
CN107868979A (zh) 一种基于恒拉速控制结构的硅单晶直径控制方法
CN105260532A (zh) 基于序列近似优化的薄板拉伸变压边力不确定性设计方法
CN109816167A (zh) 径流预报方法及径流预报装置
CN108876094A (zh) 一种清洁能源接入的电网运行综合评价方法
CN108537366A (zh) 基于最优卷积二维化的水库调度方法
CN105608295A (zh) 焦化炉压力的多目标遗传算法与rbf神经网络优化建模方法
CN107729988B (zh) 基于动态深度置信网络的蓝藻水华预测方法
CN105184400A (zh) 一种烟田土壤水分预测方法
Wang et al. A new approach of obtaining reservoir operation rules: Artificial immune recognition system
CN104732067A (zh) 一种面向流程对象的工业过程建模预测方法

Legal Events

Date Code Title Description
C06 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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131023

Termination date: 20140413