CN110393954A - 一种基于强化学习的浓密机在线控制方法 - Google Patents
一种基于强化学习的浓密机在线控制方法 Download PDFInfo
- Publication number
- CN110393954A CN110393954A CN201910636652.7A CN201910636652A CN110393954A CN 110393954 A CN110393954 A CN 110393954A CN 201910636652 A CN201910636652 A CN 201910636652A CN 110393954 A CN110393954 A CN 110393954A
- Authority
- CN
- China
- Prior art keywords
- network
- value
- concentrator
- underflow
- pump speed
- 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.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D21/00—Separation of suspended solid particles from liquids by sedimentation
- B01D21/01—Separation of suspended solid particles from liquids by sedimentation using flocculating agents
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D21/00—Separation of suspended solid particles from liquids by sedimentation
- B01D21/30—Control equipment
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D21/00—Separation of suspended solid particles from liquids by sedimentation
- B01D21/30—Control equipment
- B01D21/32—Density control of clear liquid or sediment, e.g. optical control ; Control of physical properties
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D21/00—Separation of suspended solid particles from liquids by sedimentation
- B01D21/30—Control equipment
- B01D21/34—Controlling the feed distribution; Controlling the liquid level ; Control of process parameters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Feedback Control In General (AREA)
- Separation Of Suspended Particles By Flocculating Agents (AREA)
Abstract
本发明提供一种基于强化学习的浓密机在线控制方法,能够降低时间消耗,并提高控制精度。所述方法包括:获取生产过程中所监测到的历史记录数据;建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练;通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速。本发明涉及采矿领域。
Description
技术领域
本发明涉及采矿领域,特别是指一种基于强化学习的浓密机在线控制方法。
背景技术
在冶金领域等复杂过程工业场景下,浓密机是一个被广泛应用的大型沉降工具,它通过重力沉降作用可以将低浓度的固液混合物进行浓缩形成高浓度的混合物,起到减水、浓缩的作用。
在实际生产过程中,由于浓密机运行机理复杂,难以建立数学模型,大部分的控制算法都是基于人工设计的专家系统或手工制定模糊控制器中的规则库,并辅助以传统的比例积分控制手段实现对底流泵速、絮凝剂泵速的控制。此类方法过度依赖人工经验、缺乏自适应性。
发明内容
本发明要解决的技术问题是提供一种基于强化学习的浓密机在线控制方法,以解决现有技术所存在的难以建立数学模型,浓密机控制方法过度依赖人工经验、缺乏自适应性的问题。
为解决上述技术问题,本发明实施例提供一种基于强化学习的浓密机在线控制方法,包括:
获取生产过程中所监测到的历史记录数据,所述历史记录数据包括:底流浓度、泥层高度、进料流量、进料浓度、底流泵速以及絮凝剂泵速;
建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练;
通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速。
进一步地,在建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练之前,所述方法还包括:
对获取到的底流浓度、泥层高度、进料流量、进料浓度、底流泵速以及絮凝剂泵速进行归一化放缩,其中,归一化放缩的公式为:
其中,z表示获取到的原始参量值,参量为底流浓度、泥层高度、进料流量、进料浓度、底流泵速或絮凝剂泵速;表示归一化放缩后的参量值;zmin表示参量z的最小值;zmax表示参量z的最大值。
进一步地,模型网络表示为:
其中,k为采样时间;分别表示模型网络预测的k+1时刻的底流浓度、泥层高度;Wm1和Wm2都表示模型网络的权重;tanh为激活函数;分别表示归一化放缩后的底流浓度、泥层高度、进料流量、进料浓度、底流泵速、絮凝剂泵速;表示归一化放缩后的控制动作;上标T表示矩阵转置。
进一步地,在训练模型网络的过程中调整权重Wm1和Wm2,当模型网络的损失函数的值收敛到第一预设值,则表明模型网络训练完成,其中,模型网络的损失函数表示为:
其中,Em(k)表示模型网络的损失函数;em(k)为简写形式, 分别表示实际的k+1时刻的底流浓度、泥层高度。
进一步地,在训练模型网络的过程中,采用梯度下降法来调节权值,表示为:
其中,lm表示学习率,i表示第i个权重。
进一步地,评价网络表示为:
其中,表示累计代价函数值;Wc1和Wc2都表示评价网络的权重。
进一步地,在训练评价网络的过程中,使用基于经验回放技术的损失函数训练评价网络,当评价网络的损失函数的值收敛到第二预设值,则表明评价网络训练完成,其中,评价网络的损失函数表示为:
其中,ec(k)表示评价网络的损失函数;L表示回放点数;U()表示效用函数;γ表示折扣因子。
进一步地,在训练评价网络的过程中,采用梯度下降法来调节权值,表示为:
其中,lc表示学习率,i表示第i个权重。
进一步地,所述通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速包括:
A1,随机选取控制动作 初始化j=0;
A2,通过训练好的模型网络,预测下一时刻的底流浓度、泥层高度:
其中,k为采样时间,即:当前时刻;k+1为下一时刻;
A3,已训练好的评价网络,根据预测到的下一时刻的底流浓度、泥层高度,估计下一时刻的累计代价值:
其中,
A4,根据估计得到的下一时刻的累计代价值,计算当前k时刻的累计代价值:
其中,表示第j次迭代得到的控制动作
A5,根据得到的当前k时刻的累计代价值,利用梯度下降算法对进行更新:
其中,lu为学习率;
A6,j=j+1;
A7,返回执行步骤A1-A6,直到的变化值小于预设的第三阈值,或者当前的迭代次数大于预设的最大迭代次数,此时的为最优控制动作:底流泵速、絮凝剂泵速。
A8,对进行反归一化,得到最终的最优控制动作u:
其中,⊙表示按位乘,umax、umin、umid分别表示控制动作的上限、下限、上限与下限的均值。
进一步地,在建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练之后,所述方法还包括:
获取实际充填过程中所监测到的实时监测数据;
利用所述实时监测数据对评价网络进行学习训练;或,
在建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练之后,所述方法还包括:
建立浓密机仿真模型;
使用建立好的浓密机仿真模型验证所述控制模型的有效性。
本发明的上述技术方案的有益效果如下:
上述方案中,获取生产过程中所监测到的历史记录数据,所述历史记录数据包括:底流浓度、泥层高度、进料流量、进料浓度、底流泵速以及絮凝剂泵速;建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练;通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速。这样的控制模型不仅具有自适应、不依赖于精确数学模型的特点,同时相比较已有的自适应动态规划算法,该控制模型在双网结构中去掉了动作网络,直接采用梯度下降迭代算法求解控制动作,有着更少的时间消耗和更高的控制精度,且结构简单。
附图说明
图1为本发明实施例提供的浓密机系统的结构示意图;
图2为本发明实施例提供的基于强化学习的浓密机在线控制方法的流程示意图;
图3为本发明实施例提供的基于强化学习的浓密机在线控制方法的详细流程示意图;
图4为本发明实施例提供的模型网络的结构示意图;
图5为本发明实施例提供的噪音量变化曲线示意图一;
图6(a)、(b)分别为本发明实施例提供的HCNVI与其他ADP算法的底流浓度、效用值对比示意图;
图7(a)、(b)分别为本发明实施例提供的在HDP算法、HCNVI算法中引入经验回放对效用值的影响示意图;
图8为本发明实施例提供的HDP与HCNVI在时间消耗上的对比示意图一;
图9为本发明实施例提供的噪音量变化曲线示意图二;
图10(a)、(b)分别为本发明实施例提供的HCNVI与其他ADP算法在波动噪声输入下的底流浓度、效用值对比示意图;
图11(a)、(b)为本发明实施例提供的噪音持续变化下经验回放对HCNVI的底流浓度、效用值影响示意图;
图12为本发明实施例提供的HDP与HCNVI在时间消耗上的对比示意图二。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明针对现有的难以建立数学模型,浓密机控制方法过度依赖人工经验、缺乏自适应性的问题,提供一种基于强化学习的浓密机在线控制方法。
为了更好地理解本发明实施例所述的基于强化学习的浓密机在线控制方法,参考图1,先对浓密机系统运行过程进行简要说明:
浓密机在采矿领域是重要的沉降分离设备,上游工段会产生浓度、流量不断波动的低浓度料浆。利用泥沙颗粒的密度大于水的特性以及絮凝剂的絮凝作用,砂粒会不断沉降并在浓密机底部形成高浓度的底流,在底流泵的压力作用下被吸入输送管道。对于浓密沉降过程的性能进行评价,其核心控制指标为底流浓度y(k)。该因素受控制输入、系统状态参量、及其他外部噪音影响,其中,控制输入u(k)为底流泵速u1(k)以及絮凝剂泵速u2(k),其中,k为采样时间,系统状态参量为底流浓度y(k)和泥层高度h(k),进料流量c1(k)、进料浓度c2(k)为外部噪音输入c(k),这是由于在部分工业场景中,上游工序产生的物料浓度、物料流量是不可控的,为了使本发明提出的浓密机控制模型具有通用性,因此将进料状态(进料流量c1(k)、进料浓度c2(k))作为噪音输入量。根据上述定义,其中为可控制输入量,为实数, 为不可控但是可观测的噪音量,为系统状态量,该参量是表征当前浓密机状态的重要参量,它可被间接控制但不作为控制目标。在工业领域中,浓密机进料颗粒大小,进料成分都会对浓密机底流浓度产生影响。不过由于此类变量无法观测且波动较小,为了简化问题,假定其保持恒定的。浓密机系统可表述为一个未知非线性系统:
[y(k+1),h(k+1)]T=f(y(k),u(k),c(k),h(k)) (1)
其中,y(k)表示底流浓度,f(·)为非线性未知函数,k为采样时间。
对于式(1),智能控制的首要目标是使底流浓度y(k),追踪其设定的理想值y*。另外,为了保证浓密机系统运行安全与仪器寿命,控制输入必须满足一定的限制条件。综合上述指标因素,可以将浓密机控制问题转化为有约束的最优化问题:
其中,J(k)表示累计代价值函数,是一种评价值函数;U(k)表示效用函数;J(k)代表在当前状态y(k)下,执行控制输入u(k)需要承受的代价;γ∈(0,1]是折扣因子,代表系统短期控制过程中产生的惩罚值在累计惩罚项所占比重;Q>0,是一个实验参数,用于平衡(y(k)-y*)2和的比重;R是对称正定矩阵,ui min,ui max分别表示ui(k)的下限和上限,u1(k)以及絮凝剂泵速u2(k)的下限统称为:umin,上限统称为:
式(2)可以表示为式(4)贝尔曼方程的形式:
根据贝尔曼最优原则,第k时刻的最优评价值函数J*(k)满足离散哈密顿-雅可比-贝尔曼方程:
第k时刻,最优的控制输入u*(k)可以表示为:
基于式(5)、公式(6),提出一种启发式评价网络值迭代算法(HCNVI),该算法能根据浓密机系统生产过程中产生的实时监测数据x(k)进行在线学习,并产生满足Ωu={u:umin≤u≤umax}约束的控制输入量u(k),且最小化J(k)。
实施例一
如图2所示,本发明实施例提供的基于强化学习的浓密机在线控制方法
S101,获取生产过程中所监测到的历史记录数据,所述历史记录数据包括:底流浓度、泥层高度、进料流量、进料浓度、底流泵速以及絮凝剂泵速;
S102,建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练;
S103,通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速。
本发明实施例所述的基于强化学习的浓密机在线控制方法,获取生产过程中所监测到的历史记录数据,所述历史记录数据包括:底流浓度、泥层高度、进料流量、进料浓度、底流泵速以及絮凝剂泵速;建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练;通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速。这样的控制模型不仅具有自适应、不依赖于精确数学模型的特点,同时相比较已有的自适应动态规划算法,该控制模型在双网结构中去掉了动作网络,直接采用梯度下降迭代算法求解控制动作,有着更少的时间消耗和更高的控制精度,且结构简单。
本实施例中所述的基于强化学习的浓密机在线控制方法是一种新的自适应动态规划(ADP)算法,所述方法是在传统启发式动态规划(HDP)算法的基础上,建立由模型网络和评价网络组成的双网结构的控制模型,可以称为启发式评价网络值迭代算法(HCNVI),所述方法利用评价网络与模型网络分别近似累计代价函数和对浓密机生产过程进行建模与预测,所述方法不依赖于精确的数学模型,并采用梯度下降迭代算法求解式(6),实现对浓密机底流泵速和絮凝剂泵速的稳定控制,因此可以有效实现对浓密机等大型过程工业设备的控制。
为了更好地理解本发明实施例所述的一种基于强化学习的浓密机在线控制方法,对其进行详细说明,如图3所示,具体可以包括以下步骤:
步骤(1)数据获取及预处理:获取实际生产过程中各种传感器所监测到的历史记录数据,并对其进行预处理,具体可以包括以下步骤:
步骤(1-1)数据获取
获取实际充填过程中各种传感器所监测到的历史记录数据,其中,历史记录数据用于模型网络和评价网络的最初训练,模型网络的训练全部离线进行,在控制任务开始后,将不再对模型网络进行调整。
另需要说明的是:
在浓密机系统稳定运行时,还需获取实际充填过程中各种传感器所监测到的实时监测数据,并采用在线学习的方式,用实时监测数据对评价网络进行学习训练。
本实施例中,获取数据主要可以包括以下步骤:
步骤(1-1-1)获取实际充填过程中各种传感器所监测到的历史记录数据,假设,本实施例中的历史记录数据,由矿山的自动化系统记录并存储于企业数据库中,通过使用OPC服务器将企业数据中的数据进行导出,保存在本地。数据记录跨度为1个月,总数据记录约1万余条,记录间隔为1分钟。历史记录数据的属性有底流泵速u1(k)、絮凝剂泵速u2(k)、泥层高度h(k)、进料流量c1(k),进料浓度c2(k)以及底流浓度y(k)。定义x(k)=[y(k),h(k),c(k)]T,上标T表示矩阵转置。
步骤(1-1-2):浓密机系统稳定运行后,将浓密机系统产生的实时监测数据通过OPC服务器实时传送到控制模型中进行评价网络的在线学习与实时控制。
由于本实施例是通过浓密机仿真模型来验证控制模型的有效性,所以在本实施例中,用于在线学习的实时监测数据为仿真模型在运行过程中所产生的数据,也就是说,是将浓密机仿真模型产生的数据实时传送到控制模型中进行评价网络的在线学习与实时控制。
步骤(1-2)数据预处理
由于不同物理量的取值差异很大,这会导致网络无法有效学习并且造成超参数设定困难。因此需要统计由浓密机产生的数据中各参量的极值,对所有数据利用下式进行归一化放缩。
式(7)中,表示归一化放缩后的参量值,z表示获取到的原始参量值,zmin表示参量z的最小值,zmax表示参量z的最大值,参量为底流浓度、泥层高度、进料流量、进料浓度、底流泵速或絮凝剂泵速。
例如,对于底流浓度y(k),进行归一化放缩后的底流浓度值为:其它参量同理。则得到预处理后的数据属性为:底流泵速絮凝剂泵速泥层高度进料流量进料浓度以及底流浓度预处理后的x(k)表示为预处理后的u(k)表示为
步骤(2)建立控制模型
建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练。
需要说明的是:图3中的经验数据就是历史记录数据,本实施例将历史记录数据暂存,在训练模型时批量取出用于控制模型更新,其中,模型网络训练过程是离线进行的,不需要使用经验回放,模型网络在图中只起到预测和的作用。
本实施例中,所述建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练,具体可以包括以下步骤:
步骤(2-1)建立模型网络
建立模型网络用来对浓密机生产过程进行建模,即根据当前系统状态、外部噪音量、控制输入、预测下一时刻底流浓度和泥层高度变化。模型网络的训练全部离线进行,在控制任务开始后,将不再对模型网络进行调整。主要包括以下步骤:
步骤(2-1-1)模型网络的自变量为因变量为下一时刻的底流浓度和泥层高度初始化模型网络的网络参数,设置网络的学习率为lm=0.01,批训练的数据个数predict_batch_size=32;
步骤(2-1-2)构造一个3层神经网络结构的模型网络,其基本结构如图4所示,模型网络可以表示为:
其中,分别表示模型网络预测的k+1时刻的底流浓度、泥层高度;Wm1和Wm2都表示模型网络的权重并且均初始化为0~1之间的随机数;tanh为激活函数;且模型网络的输入层包含6个节点,隐含层包含20个节点,输出层包含2个节点。
步骤(2-1-3)确定模型网络的损失函数,训练模型网络;其中,模型网络的损失函数表示为:
其中,Em(k)表示模型网络的损失函数;em(k)为简写形式, 分别表示实际的k+1时刻的底流浓度、泥层高度。
步骤(2-1-4)设置训练次数predict_epoch=50,使用Adam(adaptive momentestimation,自适应矩估计)优化器作为神经网络的优化器,并采用梯度下降法来调节权重Wm1和Wm2,如式(11)所示,当模型网络的损失函数的值收敛到第一预设值,例如,0.08,则表明模型网络已训练完成。
其中,lm表示学习率;i表示第i个权重,i的取值为1、2。
步骤(2-2)建立评价网络
评价网络用于对累积代价函数进行近似。在训练评价网络时,将短期内的监测数据共同用于评价网络的训练,主要包含以下步骤:
步骤(2-2-1)评价网络的自变量为因变量为累计代价值。初始化模型网络的网络参数,设置网络的学习率为lc=0.01,批训练的数据个数batch_size=1;
步骤(2-2-2)构造一个3层神经网络结构的评价网络,评价网络表示为:
其中,表示累计代价函数值;Wc1和Wc2都表示评价网络的权重且均初始化为0~1之间的随机数;tanh为激活函数,评价网络的输入层包含4个节点,隐含层包含14个节点,输出层为1个节点。
步骤(2-2-3)确定损失函数,训练评价网络。
本实施例中,使用时序差分误差作为评价网络的损失函数,损失函数表示为:
其中,ec(k)表示评价网络的损失函数;U()表示效用函数;γ表示折扣因子。
本实施例中,由于浓密机具有运行缓慢、状态迁移滞后的特点,为了维持评价网络的局部梯度差异以增强控制模型的收敛速度,在训练评价网络的过程中,本实施例利用经验回放技术以对短期内的监测数据进行回放训练,以增加评价网络训练的准确性和收敛速度,且增强了收敛稳定性。
本实施例中,经验回放方法将式(13)所述的损失函数修改为式(14),得到基于经验回放技术的损失函数:
其中,L表示回放点数,一般设置为2。
步骤(2-2-4)设置训练次数Nc=500,使用Adam优化器作为神经网络的优化器,并采用梯度下降法来调节权重Wc1和Wc2,如式(15)所示,当评价网络的损失函数的值收敛到第二预设值∈c,例如,∈c=0.001,则表明评价网络已训练完成。
其中,lc表示学习率;i表示第i个权重,i的取值为1、2;由于浓密机所处环境的外界噪音是不断波动的,lc需设定为固定值以保持学习能力。
步骤(3)利用梯度下降算法计算控制动作
通过输入第k时刻(当前时刻)的底流浓度、泥层高度、进料流量、进料浓度、底流泵速或絮凝剂泵速,得到第k时刻的最优控制动作输出:底流泵速、絮凝剂泵速;主要包含以下步骤:
步骤(3-1)随机选取控制动作 初始化j=0;;
步骤(3-2)预测下一(k+1)时刻的系统状态:
其中,k为采样时间,即:当前时刻;k+1为下一时刻;
步骤(3-3)已训练好的评价网络,根据预测到的下一时刻的底流浓度、泥层高度,估计下一时刻的累计代价值:
其中,
需要说明的是:由于真实工业环境下进料噪音都是连续变化的,很少出现突变,因此本实施例采用当前时刻噪音c(k)来充当下一时刻噪音c(k+1);
步骤(3-4)根据估计得到的k+1时刻的累计代价值,计算当前k时刻的累计代价值:
其中,表示第j次迭代得到的控制动作
步骤(3-5)根据得到的当前k时刻的累计代价值,利用梯度下降算法对进行更新:
其中,lu为学习率,设置为0.4,令j=j+1;
步骤(3-6)一直循环步骤(3-2)、(3-3)、(3-4)、(3-5),直到u的变化值小于预设的第三阈值∈a,例如,∈a=0.0001,或者当前迭代的次数j大于最大迭代次数Na(例如,4000),此时的为最优控制动作;
步骤(3-7)对进行反归一化,得到最终的最优控制动作u:
其中,⊙表示按位乘,umax、umin、umid分别表示控制动作的上限、下限、上限与下限的均值。
步骤(4)建立浓密机仿真模型,并使用建立好的浓密机仿真模型验证本发明提出的控制模型的有效性。由于在真实工业场景下进行浓密机控制实验成本较高,本发明采用浓密机仿真模型验证所提出控制模型的有效性,主要包括以下步骤:
步骤(4-1)根据已有的理论知识建立浓密机仿真模型,该浓密机仿真模型建立在四个假设基础上:(a)进料都是球形颗粒;(b)絮凝剂在浓密机的静态混合器中作用完全;(c)流体的扩散以固液混合物形式进行;(d)忽略颗粒间相互作用、浓密机中把机中轴的影响,主要包含以下步骤:
步骤(4-1-1)确定仿真模型推导中出现的变量,如表1,表2,表3所示:
表1仿真模型常量
变量 | 含义 | 量纲 | 参考值 |
ρ<sub>s</sub> | 干砂密度 | kg/m<sup>3</sup> | 4150 |
ρ<sub>e</sub> | 介质表观密度 | kg/m<sup>3</sup> | 1803 |
μ<sub>e</sub> | 悬浮体系的表观粘度 | Pa·s | 1 |
d<sub>0</sub> | 进料颗粒直径 | m | 0.00008 |
p | 平均浓度系数 | 无 | 0.5 |
A | 浓密机横截面积 | m<sup>2</sup> | 300.5 |
k<sub>s</sub> | 絮凝剂作用系数 | s/m<sup>2</sup> | 0.157 |
k<sub>i</sub> | 压缩层浓度系数 | m<sup>3</sup>/s | 0.0005*3600 |
K<sub>i</sub> | 进料流量与进料泵频的系数 | m<sup>3</sup>/r | 50/3600 |
K<sub>u</sub> | 底流流量与底流泵频的系数 | m<sup>3</sup>/r | 2/3600 |
K<sub>f</sub> | 絮凝剂流量与絮凝剂泵频的系数 | m<sup>3</sup>/r | 0.75/3600 |
θ | 压缩时间 | s | 2300 |
表2参量定义
变量 | 含义 | 量纲 | 初始值 | 补充说明 |
f<sub>i</sub>(t) | 进料泵频 | HZ | 40 | 扰动量 |
f<sub>u</sub>(t) | 底流泵频 | HZ | 85 | 控制量 |
f<sub>f</sub>(t) | 絮凝剂泵频 | HZ | 40 | 控制量 |
c<sub>i</sub>(t) | 进料浓度 | kg/m<sup>3</sup> | 73 | 扰动量 |
h(t) | 泥层高度 | m | 1.48 | 状态量 |
c<sub>u</sub>(t) | 底流浓度 | kg/m<sup>3</sup> | 680 | 目标量 |
表3部分变量计算方法
步骤(4-1-2)根据已有的理论知识,进行仿真模型推导。
首先可得泥层高度与泥层内平均单位体积含固量之间的关系:
根据固体守恒定律,泥层内固体质量变化量等于由进料导致泥层内固体量增加量与底流导致泥层内固体减少量的差,因此可以建立泥层内平均单位体积含固量与粒子沉降速度的关系:
对上式做变形得到下式:
联立式(19)、式(17),可得泥层高度h(t)与底流浓度cu(t)的一阶变化率:
在这个仿真模型中,絮凝剂泵速ff和底流泵速fu是控制输入u=[fu,ff]T,进料泵速fi和进料浓度ci是外部干扰量c=[fi,ci]T,底流浓度cu为系统追踪变量y=cu。理想的控制系统能够在外界干扰量c不断波动下,通过在合理范围内调节u,驱使y追踪其设定值y*。根据真实生产数据对变量取值范围进行限制,umin=[40,30],umax=[120,50],hmin=0.75,hmax=2.8,ymin=280,ymax=1200,cmin=[40,30],cmax=[120,50],y*=680。
步骤(4-2)使用建立好的浓密机仿真模型验证本发明提出的控制模型的有效性。主要在两种类型噪音量c(k)输入下,展示本发明提出的控制模型的控制效果,并与其他ADP算法,主要包括HDP、双启发式动态规划(Dual Heuristic Programming,DHP)、穿插学习策略迭代(Interleaved Learning Policy Iteration,ILPL)进行比较。主要包含以下步骤:
步骤(4-2-1):恒定-阶跃型噪音输入下浓密机控制仿真实验:第一个实验首先设置干扰量输入c为恒定值,并在某一时刻为其增加阶跃突变,噪音输入量如图5所示。该实验用来验证控制模型能否在浓密机处于非最优态下,快速寻找到u*,使被控模型达到理想收敛稳态。主要包含以下步骤:
步骤(4-2-1-1):在这个仿真实验中设置控制器参数如下:迭代轮次T=270,仿真步长Td=12min,Q=0.004,γ=0.6,Na=4000,Nc=500,∈c=0.001,∈a=0.0001,lm=0.01,lc=0.01,lu=0.4。控制效果如图6所示,其中HDP、DHP算法也使用经验回放,回放点数为2。实验中HDP、ILPL、HCNVI的评价网络结构相同,且网络参数初始化为相同数值。
根据实验结果可以发现,对于不同控制算法,由于网络参数初始值均为随机设定值,训练初期底流浓度有较大幅度的波动,且在设定值两侧持续震荡。随着各个控制模型的学习,系统状态与网络参数不断趋于平稳,直到某一时刻,底流浓度开始稳定并与设定值重合,且不再产生波动,此时控制模型参数也不再发生变化,被控系统和控制模型同时收敛到最优态。从效用值变化曲线也可以看出,早期由于底流浓度与其设定值偏差较大,效用值较高,但是随着模型与系统趋于稳态,U(k)不断缩减直到接近于0的位置。到达270分钟时,系统进料浓度、进料流量发生突变,底流浓度无法维持稳态,开始远离设定值。控制模型根据噪音量改变后的系统所产生的轨迹数据重新训练,将底流浓度拉回设定值位置。由于在第一阶段控制模型已经到达过一次稳态,在第二阶段仅需要较少的迭代次数就可以使系统重归理想收敛稳态。通过观察不同控制算法产生的系统轨迹,可以发现不同控制算法到达最优态所需的时间有较大差别,且在收敛到最优态的过程中,底流浓度的波动也有较大差异。在实验第一阶段,为使系统达到稳态,HCNVI算法所需要的迭代次数更少,训练过程中产生的底流浓度振幅也更小。并且在噪音量改变后,HCNVI算法可以迅速地使模型重归最优态,且底流浓度几乎未发生大幅度波动。
步骤(4-2-1-2)为了验证经验回放技术对控制算法性能的影响,本发明分别对比了无经验回放、经验回放数量为2的情况下,HDP、HCNVI的控制性能。对比结果如图7所示。在本实验中,仅比较了两种算法的效用值变化,效用值越快地收敛到0说明算法控制效果越佳。通过观察图7(a)和图7(b)中无经验回放情况下的效用值变化曲线,可以发现曲线波动较大,相比于使用短期经验回放,无经验回放情况下控制模型需要更多的迭代轮次才能够使系统达到收敛。特别是在图7(b)的HCNVI的实验中,270分钟时系统噪音输入量改变,效用值开始剧增,底流浓度开始偏离设定值,被控系统几乎无法收敛,但在增加了经验数据回放后效用函数值可以快速收敛至最低点。该实验结果表明经验回放技术无论对HDP算法还是HCNVI算法都有很好的促进作用。对于在线学习的控制算法,增加对短期历史数据点的经验回放可以通过维护局部梯度准确性,有效地指导评价网络学习正确的评价值,增强模型的收敛能力。
步骤(4-2-1-3)本文进行了十组实验来对比HCNVI算法在时间上的优势。选取HDP算法作为参考对象,T=270,结果如图8所示。由于每次实验中网络初始值不同,系统运行轨迹以及模型训练过程也不同,因此每组实验中模型学习以及控制所需的累积时间略有差异。但是从多次实验结果可以看出,由于HCNVI算法中去掉了动作网络,仅需要训练评价网络,所以模型整体训练时间大大缩减,尽管HCNVI算法中计算控制输入所需时间相比于HDP算法直接利用动作网络前向传播求解控制动作所需时间长,但是HCNVI算法总消耗时间明显少于HDP算法。
步骤(4-2-2)高斯噪音波动输入下浓密机控制仿真实验。真实工业场景下,浓密机的进料浓度和流量是实时波动的。步骤(4-2-1)中仿真模型的进料状态是基本恒定的,只在某一时刻产生突变,其目的是为了更好地观察不同控制算法的收敛速度。在步骤(4-2-2)实验中,进料流量和进料浓度两个噪音量存在持续波动,用来模仿真实工业场景。噪音输入的单步变化量服从高斯分布:
c(k+1)=c(k)+Δc Δc~N(μ=0,∑=diag(0.6,0.6)) (26)
本实施例中,进料波动变化如图9所示。
步骤(4-2-2-1)HCNVI控制器参数与步骤(4-2-1-1)相同,迭代轮次T=270,仿真步长Td=120s。利用该仿真模型再次对比HCNVI与其他算法控制性能的差异,结果如图10(a)、(b)所示。
通过观察实验结果发现在环境噪音连续变化条件下,浓密机底流浓度会发生持续震荡。随着对模型参数的不断训练,各个算法的控制性能趋于平稳,由于进料噪音导致的底流浓度波动稍有减弱。对比不同控制算法的控制性能,可以发现HCNVI相比于其他ADP算法,能够更快地将底流浓度锁定在设定值临域范围内,且浓度振幅小于其他ADP算法。从效用值变化曲线也可以看出,相比于其他算法,HCNVI算法的效用值整体较小,且在训练后期几乎0。
步骤(4-2-2-2)进一步分析不同控制算法之间的性能。表4给出了不同算法在恒定-阶跃型噪音输入下步骤(4-2-1-1)实验(称为:实验一)和高斯噪音波动输入下步骤(4-2-2-1)实验(称为:实验二)中底流浓度控制性能指标对比结果。相比其他算法,HCNVI算法可以更好地控制底流浓度稳定在其设定值附近,其控制总体稳定性(由MSE、IAE体现)、控制鲁棒性(由MAE体现)更佳。在过程工业控制场景中,控制系统的MAE指标尤为重要,某一工序的物料性质发生剧烈波动会使下游物料加工工序出现连带波动,严重影响生产的稳定性和最终产品的质量。HCNVI算法在MAE指标上的优势证实了其在过程工业控制问题中的适用性。
表4不同控制算法之间性能分析
步骤(4-2-2-3)在噪音持续变化下,算法中引入经验回放技术对于控制性能的改善也是很明显的。图11(a)、(b)对比了在HCNVI算法中,不使用经验回放和设置经验回放数量为2两种情况下的模型控制模型性能。实验结果表明,经验回放技术在环境噪音持续变化下仍对模型收敛速度有重要促进作用。
步骤(4-2-2-4)再次对比了在噪音持续变化下,两个算法在时间消耗上的差异性。本发明中仍然重复十次实验,每次实验的迭代轮次为270轮。实验结果如图12所示,在噪音持续变化环境下,HCNVI算法在时间消耗上的优势稍有削减,这是由于当环境存在持续外部扰动时,评价网络的参数也伴随着持续波动,每轮迭代中评价网络训练的时间增加,而这一部分的时间损耗在HCNVI和HDP中是共有的,因此使得HCNVI算法在整体时间消耗的优势减少,但总时间消耗仍少于HDP算法。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种基于强化学习的浓密机在线控制方法,其特征在于,包括:
获取生产过程中所监测到的历史记录数据,所述历史记录数据包括:底流浓度、泥层高度、进料流量、进料浓度、底流泵速以及絮凝剂泵速;
建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练;
通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速。
2.根据权利要求1所述的基于强化学习的浓密机在线控制方法,其特征在于,在建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练之前,所述方法还包括:
对获取到的底流浓度、泥层高度、进料流量、进料浓度、底流泵速以及絮凝剂泵速进行归一化放缩,其中,归一化放缩的公式为:
其中,z表示获取到的原始参量值,参量为底流浓度、泥层高度、进料流量、进料浓度、底流泵速或絮凝剂泵速;表示归一化放缩后的参量值;zmin表示参量z的最小值;zmax表示参量z的最大值。
3.根据权利要求2所述的基于强化学习的浓密机在线控制方法,其特征在于,模型网络表示为:
其中,k为采样时间;分别表示模型网络预测的k+1时刻的底流浓度、泥层高度;Wm1和Wm2都表示模型网络的权重;tanh为激活函数;分别表示归一化放缩后的底流浓度、泥层高度、进料流量、进料浓度、底流泵速、絮凝剂泵速;表示归一化放缩后的控制动作;上标T表示矩阵转置。
4.根据权利要求3所述的基于强化学习的浓密机在线控制方法,其特征在于,在训练模型网络的过程中调整权重Wm1和Wm2,当模型网络的损失函数的值收敛到第一预设值,则表明模型网络训练完成,其中,模型网络的损失函数表示为:
其中,Em(k)表示模型网络的损失函数;em(k)为简写形式, 分别表示实际的k+1时刻的底流浓度、泥层高度。
5.根据权利要求4所述的基于强化学习的浓密机在线控制方法,其特征在于,在训练模型网络的过程中,采用梯度下降法来调节权值,表示为:
其中,lm表示学习率,i表示第i个权重。
6.根据权利要求5所述的基于强化学习的浓密机在线控制方法,其特征在于,评价网络表示为:
其中,表示累计代价函数值;Wc1和Wc2都表示评价网络的权重。
7.根据权利要求6所述的基于强化学习的浓密机在线控制方法,其特征在于,在训练评价网络的过程中,使用基于经验回放技术的损失函数训练评价网络,当评价网络的损失函数的值收敛到第二预设值,则表明评价网络训练完成,其中,评价网络的损失函数表示为:
其中,ec(k)表示评价网络的损失函数;L表示回放点数;U()表示效用函数;γ表示折扣因子。
8.根据权利要求7所述的基于强化学习的浓密机在线控制方法,其特征在于,在训练评价网络的过程中,采用梯度下降法来调节权值,表示为:
其中,lc表示学习率,i表示第i个权重。
9.根据权利要求8所述的基于强化学习的浓密机在线控制方法,其特征在于,所述通过训练好的模型网络预测下一时刻的底流浓度和泥层高度,且训练好的评价网络根据预测到的所述下一时刻的底流浓度、泥层高度,估计所述下一时刻的累计代价值,根据估计得到的下一时刻的累计代价值,计算当前时刻的累计代价值,根据得到的当前时刻的累计代价值,利用梯度下降迭代算法确定当前时刻最优控制动作:底流泵速、絮凝剂泵速包括:
A1,随机选取控制动作 初始化j=0;
A2,通过训练好的模型网络,预测下一时刻的底流浓度、泥层高度:
其中,k为采样时间,即:当前时刻;k+1为下一时刻;
A3,已训练好的评价网络,根据预测到的下一时刻的底流浓度、泥层高度,估计下一时刻的累计代价值:
其中,
A4,根据估计得到的下一时刻的累计代价值,计算当前k时刻的累计代价值:
其中,表示第j次迭代得到的控制动作
A5,根据得到的当前k时刻的累计代价值,利用梯度下降算法对进行更新:
其中,lu为学习率;
A6,j=j+1;
A7,返回执行步骤A1-A6,直到的变化值小于预设的第三阈值,或者当前的迭代次数大于预设的最大迭代次数,此时的为最优控制动作:底流泵速、絮凝剂泵速。
A8,对进行反归一化,得到最终的最优控制动作u:
其中,⊙表示按位乘,umax、umin、umid分别表示控制动作的上限、下限、上限与下限的均值。
10.根据权利要求1所述的基于强化学习的浓密机在线控制方法,其特征在于,在建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练之后,所述方法还包括:
获取实际充填过程中所监测到的实时监测数据;
利用所述实时监测数据对评价网络进行学习训练;或,
在建立由模型网络和评价网络组成的双网结构的控制模型,并利用获取到的历史记录数据对所述模型网络和评价网络进行训练之后,所述方法还包括:
建立浓密机仿真模型;
使用建立好的浓密机仿真模型验证所述控制模型的有效性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910636652.7A CN110393954B (zh) | 2019-07-15 | 2019-07-15 | 一种基于强化学习的浓密机在线控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910636652.7A CN110393954B (zh) | 2019-07-15 | 2019-07-15 | 一种基于强化学习的浓密机在线控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110393954A true CN110393954A (zh) | 2019-11-01 |
CN110393954B CN110393954B (zh) | 2020-05-12 |
Family
ID=68325568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910636652.7A Active CN110393954B (zh) | 2019-07-15 | 2019-07-15 | 一种基于强化学习的浓密机在线控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110393954B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112084645A (zh) * | 2020-09-02 | 2020-12-15 | 沈阳工程学院 | 一种基于混合迭代adp方法的锂离子电池储能系统能量管理方法 |
CN112445136A (zh) * | 2020-12-16 | 2021-03-05 | 北京科技大学 | 一种基于连续时间神经网络的浓密机预测控制方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106933099A (zh) * | 2017-03-06 | 2017-07-07 | 太原理工大学 | 一种选煤厂浓缩机和压滤机药剂添加协同控制系统 |
AU2018100166A4 (en) * | 2017-02-10 | 2018-03-01 | Outotec (Finland) Oy | Method and arrangement for controlling a dewatering process |
CN109460890A (zh) * | 2018-09-21 | 2019-03-12 | 浙江大学 | 一种基于强化学习与控制性能监测的智能自愈方法 |
-
2019
- 2019-07-15 CN CN201910636652.7A patent/CN110393954B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2018100166A4 (en) * | 2017-02-10 | 2018-03-01 | Outotec (Finland) Oy | Method and arrangement for controlling a dewatering process |
CN106933099A (zh) * | 2017-03-06 | 2017-07-07 | 太原理工大学 | 一种选煤厂浓缩机和压滤机药剂添加协同控制系统 |
CN109460890A (zh) * | 2018-09-21 | 2019-03-12 | 浙江大学 | 一种基于强化学习与控制性能监测的智能自愈方法 |
Non-Patent Citations (1)
Title |
---|
王旭等: "浓缩生产过程优化控制系统的研究与应用", 《有色金属(选矿部分)》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112084645A (zh) * | 2020-09-02 | 2020-12-15 | 沈阳工程学院 | 一种基于混合迭代adp方法的锂离子电池储能系统能量管理方法 |
CN112084645B (zh) * | 2020-09-02 | 2023-06-09 | 沈阳工程学院 | 一种基于混合迭代adp方法的锂离子电池储能系统能量管理方法 |
CN112445136A (zh) * | 2020-12-16 | 2021-03-05 | 北京科技大学 | 一种基于连续时间神经网络的浓密机预测控制方法及系统 |
CN112445136B (zh) * | 2020-12-16 | 2022-02-22 | 北京科技大学 | 一种基于连续时间神经网络的浓密机预测控制方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110393954B (zh) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108764540B (zh) | 基于并行lstm串联dnn的供水管网压力预测方法 | |
US11543787B2 (en) | Networked control system time-delay compensation method based on predictive control | |
CN102122132A (zh) | 一种基于模糊神经网络的用于环境模拟系统的智能控制系统 | |
CN103226741B (zh) | 城市供水管网爆管预测方法 | |
CN109343505A (zh) | 基于长短期记忆网络的齿轮剩余寿命预测方法 | |
CN102831269A (zh) | 一种流程工业过程工艺参数的确定方法 | |
CN110090478B (zh) | 一种充填场景下的深锥浓密机智能控制方法 | |
CN112819240B (zh) | 基于物理约束的lstm模型预测页岩油产量的方法 | |
CN110393954A (zh) | 一种基于强化学习的浓密机在线控制方法 | |
CN107967803A (zh) | 基于多源数据及变权重组合预测模型的交通拥堵预测方法 | |
CN109472397B (zh) | 基于粘度变化的聚合工艺参数调节方法 | |
CN106873381A (zh) | 喷氨控制系统 | |
CN106406083A (zh) | 一种模糊控制选煤方法 | |
CN109724657A (zh) | 基于改进的Elman神经网络的水表流量计量方法和系统 | |
CN109254530A (zh) | 基于磨矿过程基础回路的无模型自适应控制方法 | |
CN116430722A (zh) | 一种供水网络系统及其控制方法、装置及存储介质 | |
Dai et al. | Dual-rate adaptive optimal tracking control for dense medium separation process using neural networks | |
CN113326968A (zh) | 基于调整pso惯性权重的母线短期负荷预测方法及装置 | |
CN108319146B (zh) | 一种径向基神经网络基于离散粒子群训练的方法 | |
CN113536685B (zh) | 一种基于神经网络外推的风速概率模型的建模方法 | |
CN103713619B (zh) | 一种磨矿过程运行控制方法的测试方法 | |
CN110609525A (zh) | 一种基于在线adp的非线性时滞系统最优控制方法 | |
Mukesh et al. | Influence of search algorithms on aerodynamic design optimisation of aircraft wings | |
CN109142167B (zh) | 一种磨矿粒度的在线鲁棒正则软测量模型的建立方法 | |
Ding et al. | The application of the fuzzy self-adaptive PID controller to the automatic operation control of water transfer canal system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |