一种原料的调度方法及装置的制造方法

一种原料的调度方法及装置的制造方法
【专利摘要】本发明提供了一种原料的调度方法及装置，该方法包括建立原料不确定调度模型；对原料不确定调度模型进行确定化处理，得到确定的带约束的多目标调度模型；对确定的带约束的多目标调度模型进行求解，得到原料调度信息，根据原料调度信息执行对混匀过程中原料的调度。通过建立原料不确定调度模型可以准确地描述实际生产过程，通过确定化和带变异的CMOSTA算法可以求取准确度很高的原料调度信息，原料调度信息可执行性很强，根据该原料调度信息执行原料的调度，可以大大提高生产效率，以及提高产品的质量。
【专利说明】
一种原料的调度方法及装置
技术领域
[0001] 本发明涉及物料调度技术领域，尤其涉及一种原料的调度方法及装置。
【背景技术】
[0002] 目前，炼铁厂在炼铁过程中，需要根据炼铁混匀过程的需求制定出原料调度方案， 通过调度方案对炼铁原料进行调度，以此提高炼铁效率及保证钢铁的质量。
[0003] 当前，在制定原料调度方案时，技术人员根据上级部门下达的原料的计划指标、各 矿物原料的总下料量以及生产工艺约束来设定堆配比方案，并根据以往的配料经验计算出 各混匀料堆所用物料种类及相应的配比，从而确定出原料的配料方案，以满足各车间对各 混匀料堆的基本要求，并通过生成的炼铁混匀原料成分情况调整原料的基本配料方案，当 在生产过程中执行该原料的配料方案时，如果某个混匀料堆的某个成分的指标不符合要求 时，技术人员则根据以往的配比经验，调整该料堆中该成分的含量，并调整其它各料堆相应 的配比以满足生产中原料的配比要求，最后根据调整的数据对原料的配比方案进行修改， 并根据修改后的原料的配比方案执行生产过程中原料的堆配比。
[0004] 由于炼铁混匀过程所需的原料种类很多，确定原料的调度方案需要很多的人工干 预，通过技术人员的经验确定原料的配料方案时，无法考虑来料时间的不确定性，如此，导 致确定出的调度方案中的原料的配比数据与实际炼铁的配比需求之间存在很大的误差，生 产效率低且炼铁质量差。

【发明内容】

[0005] 鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上 述问题的一种原料的调度方法及装置。
[0006] 依据本发明的一个方面，提供了 一种原料的调度方法，所述方法包括：
[0007] 建立原料不确定调度模型；
[0008] 对所述原料不确定调度模型进行确定化处理，得到确定的带约束的多目标调度模 型；
[0009] 对所述确定的带约束的多目标调度模型进行求解，得到原料调度信息；
[0010] 根据所述原料调度信息执行对混匀过程中原料的调度。
[0011] 可选地，所述建立原料不确定调度模型，包括：
[0012] 获取混匀过程的配料基本信息；
[0013] 采用三角模糊数，建立异常调度对应的惩罚函数；
[0014] 根据所述配料基本信息和所述异常调度对应的惩罚函数，确定混匀过程中的目标 函数；
[0015] 确定混匀过程中各种原料的约束条件；
[0016] 根据所述目标函数和所述约束条件建立原料不确定调度模型，所述原料不确定调 度模型如下式(1)所示，
[0020] 式（1)中，/U_)为异常调度对应的惩罚函数;xu为决策变量，代表第j堆中所需第 i种原料的相应下料量;C1;1为第i中原料的采购价格;&(x)为生产成本，代表完成堆计划所 需原料的总花费，fi(X)为目标函数，代表最小的生产成本;m表示原料数;η表示堆数；
[0021] Ci>2为第i种原料每吨每天所需的保存费用，_为额外的存储成本，C2(x)为总存储 成本，代表完成堆计划期间所有原料存储所消耗的费用，f 2(X)为目标函数，代表最小的总 存储成本；
[0022] s.t.为各种原料的约束条件，Qy为所配第j堆混合料堆中第τ种化学成分的含量 指标，Ay为计划第j堆混合料堆中第τ种化学成分的含量指标，qy为计划第j堆混合料堆中 第τ种化学成分含量指标的最大波动范围， Sl(t)为第i种原料当前的库存量，h为配比第j堆 所堆堆重，D」为计划中第j堆的标准堆重，为计划中相对应仏指标的最大波动范围， X1>J为 特殊物料在j堆中掺配量;为第j堆所需第i种原料相应的下料量。
[0023] 可选地，所述对所述原料不确定调度模型进行确定化处理，得到确定的带约束的 多目标调度模型，包括：
[0024] 通过模糊模拟技术对所述原料不确定调度模型中每个含有模糊函数方程的目标 函数进行处理，得到所述目标函数确定的输入输出样本集；
[0025] 对所述确定的输入输出数据样本集进行训练，得到确定的带约束的多目标调度模 型。
[0026] 可选地，所述通过模糊模拟技术对所述原料不确定调度模型中每个含有模糊函数 方程的目标函数进行处理，得到所述目标函数对应的确定的输入输出样本集，包括：
[0027] 将每个目标函数中的惩罚函数表示为预设的模糊函数方程表示式，所述预设的模 糊函数方程表示式如下式(2)所示，
[0028] ,n(X，4) (2)
[0029] 式⑵中，多为目标函数中惩罚函数对应的模糊函数，紅,"__,04表示，与父和#的 函数关系，X为一组决策变量，表示原料的下料量，#为三角模糊数表示的函数参数；
[0030] 获取三角模糊数的隶属度函数；
[0031] 确定目标函数的一组决策变量X1，对所述三角模糊数在所述三角模糊数的隶属度 函数内抽样，得到I的一系列抽样值取& = 1，2,...，/〇，1(为抽样数；
[0032] 根据所述决策变量X1和每个抽样值0，通过公式（2)得到K个确定的输出 乂 〇 = 1,2,...，夂），并求取1(个){〇 = 1,2，...，[)的平均值，得到目标函数对应于决策变量父1 的确定输出值/;
[0033] 根据所述目标函数输出值尹的计算方法，计算N次，得到N组输入数据XHlil， 2，. . .，N)和对应的N组输出数据芦（/ = 1,2,...,TV )，将所述N组输入数据和对应的所述N组输 出数据确定为目标函数确定的输入输出样本集。
[0034] 可选地，所述对所述确定的输入输出数据样本集进行训练，得到确定的带约束的 多目标模型，包括：
[0035]将所述确定的输入输出数据样本集排列为训练数据集样本；
[0036]对所述训练数据集样本，通过状态转移算法STA训练神经元网络学习过程的权重 向量，逼近所述不确定调度模型中对应的目标函数，得到确定的生产成本函数和/存储成本 函数；
[0037] 根据训练后得到的所述确定的生产成本函数和存储成本函数建立确定的带约束 的多目标调度模型。
[0038] 可选地，所述对所述确定的带约束的多目标调度模型进行求解，得到原料调度信 息，包括：
[0039]初始化设置CM0STA算法的相关参数、Pareto最优搜索解集种群Ps、Pareto非劣解 集种群Pb、最大迭代次数Gmax，最大个体违约程度值ε以及水平截集pc;
[0040] 从当前种群出发，开始当前迭代操作，对所述确定的带约束的多目标调度模型采 用CM0STA搜索策略，得到新的种群；
[0041] 判断所述新的种群是否陷入局部最优，若是，对所述种群参数采用变异操作进行 种群更新，所述局部最优为当前的Pareto非劣解集种群Ps与上一轮迭代得到的Pareto非劣 解集种群P s相等时的状态；
[0042] 判断当前迭代次数是否达到最大迭代次数Gmax，若没达到，将所有中间种群置空， 对CM0STA搜索中得到的Pareto最优搜索解集种群Ps和Pareto非劣解集种群P B执行新一轮迭 代操作，即重新搜索新种群并判断是否需要种群更新;若达到，则输出Pareto非劣解集种群 Pb作为搜索结果的非劣最优解集；
[0043] 通过T0PSIS(Technique for Order Preference by Similarity to an Ideal Solution，逼近理想解排序法)从所述搜索结果输出的非劣最优解集中选取一组非劣解，将 选取的所述一组非劣解确定为原料调度信息。
[0044] 依据本发明的另一个方面，提供了一种原料的调度装置，所述装置包括：
[0045] 建立模块，用于建立原料不确定调度模型；
[0046]处理模块，用于对所述原料不确定调度模型进行确定化处理，得到确定的带约束 的多目标调度模型；
[0047]计算模块，用于对所述确定的带约束的多目标调度模型进行求解，得到原料调度 信息；
[0048] 执行模块，用于根据所述原料调度信息执行对混匀过程中原料的调度。
[0049] 可选地，所述建立模块包括：
[0050] 获取单元，用于获取混匀过程的配料基本信息；
[0051 ]第一建立单元，用于采用三角模糊数，建立异常调度对应的惩罚函数；
[0052]第一确定单元，根据所述配料基本信息和所述异常调度对应的惩罚函数，确定混 匀过程中的目标函数；
[0053]第二确定单元，用于确定混匀过程中各种原料的约束条件；
[0054]第二建立单元，用于根据所述目标函数和所述约束条件建立原料不确定调度模 型，所述原料不确定调度模型如下式(1)所示，
[0058] 式（1)中，/(尤)为异常调度对应的惩罚函数;Xi,j为决策变量，代表第j堆中所需第 i种原料的相应下料量;C1;1为第i中原料的采购价格;&(x)为生产成本，代表完成堆计划所 需原料的总花费，fi(X)为目标函数，代表最小的生产成本;m表示原料数;η表示堆数；
[0059] Ci,2为第i种原料每吨每天所需的保存费用，爲为额外的存储成本，C2(x)为总存储 成本，代表完成堆计划期间所有原料存储所消耗的费用，f 2(X)为目标函数，代表最小的总 存储成本；
[0060] S.t.为各种原料的约束条件，Qy为所配第j堆混合料堆中第τ种化学成分的含量 指标，Ay为计划第j堆混合料堆中第τ种化学成分的含量指标，qy为计划第j堆混合料堆中 第τ种化学成分含量指标的最大波动范围， Sl(t)为第i种原料当前的库存量，h为配比第j堆 所堆堆重，D」为计划中第j堆的标准堆重，为计划中相对应仏指标的最大波动范围， X1>J为 特殊物料在j堆中掺配量;为第j堆所需第i种原料相应的下料量。
[0061 ] 可选地，所述处理模块包括：
[0062]处理单元，用于通过模糊模拟技术对所述原料不确定调度模型中每个含有模糊函 数方程的目标函数进行处理，得到所述目标函数确定的输入输出样本集；
[0063]训练单元，用于对所述确定的输入输出数据样本集进行训练，得到确定的带约束 的多目标调度模型。
[0064]可选地，所述计算模块包括：
[0065]设置单元，用于初始化设置CM0STA算法的相关参数、Pareto最优搜索解集种群Ps、 Pareto非劣解集种群Pb、最大迭代次数Gmax，最大个体违约程度值ε以及水平截集pc;
[0066] 搜索单元，从当前种群出发，开始当前迭代操作，对所述确定的带约束的多目标模 型采用CM0STA搜索策略，得到新的种群；
[0067] 更新单元，用于判断所述新的种群是否陷入局部最优，若是，对所述种群采用变异 操作进行种群更新，所述局部最优为当前的Pareto非劣解集种群P s与上一轮迭代得到的 Pareto非劣解集种群Ps相等时的状态；
[0068] 迭代判断单元，用于判断搜索单元和更新单元是否需要迭代进行，若当前迭代次 数没达到最大迭代次数Gmax，则重新启动搜索单元和更新单元;若达到，则迭代过程终止， 输出Pareto非劣解集种群Pb作为搜索结果的非劣最优解集；
[0069] 选取单元，用于通过T0PSIS(Technique for Order Preference by Similarity to an Ideal Solution，逼近理想解排序法)法从所述搜索结果输出的非劣最优解集中选 取一组非劣解，将选取的所述一组非劣解确定为原料调度信息。
[0070] 本发明提供的原料的调度方法及装置，实现建立的原料不确定调度模型，可以准 确地描述实际生产过程，通过确定化和带变异的CM0STA算法可以求取准确度很高的原料调 度信息，原料调度信息可执行性很强，根据该原料调度信息执行原料的调度，可以大大提高 生产效率，以及提高产品的质量。
【附图说明】
[0071] 通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通 技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明 的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
[0072] 图1是本发明一种实施方式的一种原料的调度方法的流程图；
[0073] 图2是本发明一种实施方式的建立原料不确定调度模型的流程图；
[0074] 图3是本发明一种实施方式的三角模糊数的隶属度函数图形；
[0075] 图4是本发明一种实施方式的对确定的带约束的多目标调度模型进行求解的流程 图；
[0076] 图5是本发明一种实施方式的获取新的种群的流程图；
[0077] 图6是本发明一种实施方式的一种原料的调度装置的结构示意图。
【具体实施方式】
[0078] 下面结合附图和实施例对本发明的实施方式作进一步详细描述。以下实施例用于 说明本发明，但不能用来限制本发明的范围。
[0079] 本发明实施例提供的原料的调度方法，针对原料场混匀过程中某些原料种类、到 达时间等不确定性问题，对不确定条件下约束优化问题进行建模，对模型确定化处理后，并 利用带变异的CM0STA(Constrained Multi-Objective State Transition Algorithm，约 束多目标状态转移算法)算法对建立的模型求解，得到准确度很高的原料的调度信息，根据 该原料的调度信息执行生产过程中的原料调度，具体过程如下。
[0080]图1是本发明一种实施方式的一种原料的调度方法的流程图，参照图1，该方法包 括以下S101-S104步骤。
[0081 ]步骤S101:建立原料不确定调度模型。
[0082]在本发明实施例中，上述原料可以是炼铁所需的原料，建立上述原料不确定模型， 可以真实反映出炼铁所需的原料种类、品位、库存量、来料时间及配料时间的范围，其中，在 建立上述原料不确定调度模型时，具体可通过以下S1011-S1015步骤来完成，如图2所示。 [0083]步骤S1011:获取混匀过程的配料基本信息。
[0084] 上述获取混匀过程的配料基本信息时，首先针对原料场混匀过程的实际问题，分 析现场工况，包括分析现场混匀所需的原料种类，计划内混匀料堆数，每种原料各个化学成 分含量，每种原料的采购价格，每种原料计划内进料量，计划内混匀料堆中各个化学成分含 量指标，根据分析结果建立原料分析表，该原料分析表示意性地如表1所示。
[0085] 表 1
[0087] 建立上述原料分析表后，可以通过该原料分析表及实际生产需求，确定出计划内 购买原料的种类、每种原料的成分指标、每种原料的购买量、每种原料的来料时间，每种原 料的不确定配料时间、混匀料堆的每种物料含量范围、各个混匀料堆应堆重范围、每堆物料 特殊物料使用信息及堆料场库存信息，最终生成混匀过程的配料基本信息。
[0088] 上述计划可以是周计划、月计划、季度计划或年计划，具体可根据实际生产的需求 确定，为了便于清楚地说明本发明实施例提供的数据，无特殊说明的情况下，上述计划均指 月计划。
[0089] 在原料的调度过程中，异常调度原料都会对生产过程产生影响，以至带来额外的 损失，因此，需要在调度原料时考虑异常调度的情况。其中，异常调度主要包括提前调度或 拖延调度。
[0090] 步骤S1012:采用三角模糊数，建立异常调度对应的惩罚函数。
[0091]在建立异常调度对应的惩罚函数前，由于混匀过程中存在来料时间的不确定性， 因此，引入三角模糊数，描述原料不确定调度模型中所有的模糊函数方程的不确定参数。 [0092]假设，用三角模糊数1=( &1，&[?，&2)"来表示生产过程中原料的不确定参数，则八对 应的隶属度函数为
[0094]上述三角模糊数对应的隶属度函数图形如图3所示。
[0095]在建立异常调度对应的惩罚函数时，首先需要确定提前来料的惩罚费用权重，延 迟来料的惩罚费用权重，以及确定提前来料的惩罚费用，延迟来料的惩罚费用;然后根据提 前来料的惩罚费用权重、延迟来料的惩罚费用权重、提前来料的惩罚费用、延迟来料的惩罚 费用建立异常调度对应的惩罚函数/〇〇,惩罚函数如下式(4)所示，
[0097] 式(4)中，《:为提前来料的惩罚费用权重，w2为延迟来料的惩罚费用权重，式为需额 外提供的存储成本为需额外提供的生产成本。其中，定义区间[e^U]为原料i(i = l， 2，...，m)的交货期窗口，其中ei为原料i的最早交货期为原料i的最晚交货期；用三角模 糊变量(? 表示原料i(i = l，2,…，m)的配料时间，Cf，，0/依次为原料 i(i = l，2, . . .，m)的最早配料时间、最可能配料时间以及最晚配料时间；若$ 称原料i 提前，保证了原料在其交货期窗口内完工，但需提供额外的存储成本尾，焉]ηκυφ,Α-Γ,}; 若G</,，称原料i滞后，影响了正常的原料堆料过程，则需提供不可估量损失成本?；, % = maxio, ^ - Ο. }0
[0098] 在保证炼铁效率以及钢铁质量的前提下，控制炼铁成本也是制定原料调度方案所 必需考虑的关键指标。根据炼铁厂对成本控制的要求，可以确定出炼铁过程的目标函数。 [0099]步骤S1013:根据所述配料基本信息和所述异常调度对应的惩罚函数，确定混匀过 程中的目标函数。
[0100] 上述目标函数主要指调度物料所花费的成本控制函数，在本发明中，优选地，上述 目标函数主要指调度物料所花费的生产成本函数和/或存储成本函数。
[0101] 在实际混匀过程中，由于各种原料参数受到生产条件的影响，因此，在建立原料的 调度方案时需要考虑各种原料的约束条件。
[0102] 步骤S1014:确定混匀过程中各种原料的约束条件。
[0103] 在确定混匀过程中各种原料的约束条件时，需要确定出混匀过程中的原料的成分 约束条件、堆重约束条件、库存约束条件及特殊原料约束条件，每种约束条件的具体的确定 过程如下。
[0104]根据混匀过程中每种原料的下料量和每种原料中各化学成分的含量，计算混匀料 堆中每种物料的成分含量，要求配完以后的每一个混匀料堆任何一个时刻的成分都必须满 足指标要求，确定出原料的成分约束条件，原料的成分约束条件如下式(5)所示；
[0106] 式(5)中，QCPm为第j堆混匀料中某个化学成分所占的百分比，ai_cpn(i = l，2,…， m)为第i种原料中某种化学成分的含量;Xu为第j堆所需第i种原料相应的下料量;为计 划中某种原料中某种化学成分的含量指标 ;qi为计划中某种原料的某种化学成分含量指标 的最大波动范围
表示第j堆混匀料中某种化学成分的总含量；
表示第j 堆混匀料的总堆重。
[0107] 上述混匀料堆中通常包括铁、硅、钙、硫及磷等化学成分，当混匀料堆中铁比、硅 比、钙比、硫比及磷比等达到炼铁所需的比例要求时，生产出的钢铁才可以满足质量要求。
[0108] 在实际炼铁过程中，受原料场空间的限制，各混匀料堆堆重也受到一定限制。根据 原料场物料容量，计算原料场各堆物料的堆重范围，可以确定出堆重约束条件，堆重物料约 束条件如式(6)所示；
[0109] bj彡Dj土pj，即 |bj-Dj | 彡pj，j = l，2, · · ·，n (6)
[0110]
为实际第j堆所堆堆重，Dj为计划中第j堆的标准堆重，pj为计 划中相对应Dj指标的最大波动范围。
[0111] 根据原料场各原料所需最大用量和各堆中所需每种原料的下料量，计算每种原料 当前的库存量，确定出库存约束条件，库存约束条件如下式(7)所示；
[0113] 式(7)中，xu为决策变量，代表第j堆中所需第i种原料的下料量，Sl(t)为第i种原 料当前的库存量。
[0114] 统计可重复利用的烧结物料，确定出特殊原料约束条件，假设只有第一种物料为 特殊原料，则特殊原料约束条件如式(8)所示；
[0115] 返矿的约束条件是xi,j>0,其余Xi,j彡0，i = 2,3, · · ·，m (8)
[0116] 对于一个实际的混匀堆料过程，考虑到一些过程参数同样受到生产工艺条件的约 束的情况，因此，根据上述确定出的原料的各种约束条件可以生成该原料不确定调度模型 的约束条件，如式(9)所示；
[0118] 式(9)中，Qy为所配第j堆混合料堆中第τ种化学成分的含量指标，Ay为计划第j 堆混合料堆中第τ种化学成分的含量指标，qy为计划第j堆混合料堆中第τ种化学成分含量 指标的最大波动范围，Sl(t)为第i种原料当前的库存量，h为配比第j堆所堆堆重，D」为计划 中第j堆的标准堆重，Pj为计划中相对应Dj指标的最大波动范围，XI, j为特殊物料在j堆中掺 配量;xu为第j堆所需第i种原料相应的下料量;m表示原料数;η表示堆数。
[0119] 上述特殊原料约束条件用以维持返矿平衡，提高资源的有效利用率，从而保证炼 铁资源的合理利用。通过上述原料的成分约束条件、上述堆重约束条件、上述库存约束条件 及上述特殊原料约束条件，可以真实反映出原料场各堆物料的堆重范围，每堆料堆中各个 化学成分的比例范围以及每种原料的库存可用量，如此，在生产过程中，可以根据各个原料 的约束条件及时填充原料，避免料堆无料、出现生产中断的情况，从而可以提高产品的生产 效率和产品质量。
[0120] 步骤S1015:根据目标函数和约束条件建立原料不确定调度模型，原料不确定调度 模型如式(1)所示，
[0124] 式（1)中，/(I)为异常调度对应的惩罚函数;Xu为决策变量，代表第j堆中所需第 i种原料的相应下料量;C1;1为第i中原料的采购价格;&(x)为生产成本，代表完成堆计划所 需原料的总花费，fi(X)为目标函数，代表最小的生产成本;m表示原料数;η表示堆数；
[0125] Ci,2为第i种原料每吨每天所需的保存费用，我为额外的存储成本，C2(x)为总存储 成本，代表完成堆计划期间所有原料存储所消耗的费用，f 2(X)为目标函数，代表最小的总 存储成本；
[0126] s.t.为各种原料的约束条件。
[0127] 在实际炼铁过程中，由于炼铁所需的含铁原料来料时间存在不确定性，导致根据 上述原料不确定调度模型难以进行原料的调度，因此，需要通过以下S102及S103步骤确定 出原料的调度信息，以使根据该原料的调度信息可以准确地执行炼铁过程中原料的调度。
[0128] 步骤S102:对原料不确定调度模型进行确定化处理，得到确定的带约束的多目标 调度模型。
[0129] 上述对原料不确定调度模型进行确定化处理时，具体可通过以下S1021-S1022步 骤进行确定化。
[0130] 步骤S1021:通过模糊模拟技术对原料不确定调度模型中每个含有模糊函数方程 的目标函数进行处理，得到目标函数确定的输入输出样本集。
[0131] 当上述建立的原料不确定调度模型包括的目标函数为fKX)及f2(X)时，由于异常 调度对应的惩罚函数Λ幻中存在不确定参数，即/(幻是不确定配料时间的函数，因此 /(X)是一个不确定函数，当采用三角模糊数表示/(X)中存在的不确定参数时，得到的 /(Υ)是一个关于三角模糊数的模糊函数方程，从而目标函数心0)及f 2(X)是含有模糊函数 方程的不确定函数。
[0132] 在本发明实施中，通过模糊模拟技术对原料不确定调度模型中每个含有模糊函数 方程的目标函数进行处理，得到目标函数确定的输入输出样本集时，具体可通过以下过程 进行处理。
[0133] 第一步:将每个目标函数中的惩罚函数表示为预设的模糊函数方程表示式，预设 的模糊函数方程表示式如下式(2)所示，
[0134] (2)
[0135] 式(2)中，罗为目标函数中惩罚函数对应的模糊函数，心__,.1>，0表示夕与父和卢的 函数关系，X为一组决策变量，表示原料的下料量，#为三角模糊数表示的函数参数；
[0136] 第二步:获取三角模糊数的隶属度函数；
[0137] 第三步：确定目标函数的一组决策变量X1，对三角模糊数在三角模糊数的隶属度 函数内抽样，得到#的一系列抽样值劣0=1,2,... 5[)，K为抽样数；
[0138] 第四步：根据决策变量X1和每个抽样值贫，通过公式（2)得到Κ个确定的输出 （々=1,2,...,/〇，并求取1(个.以々=1,2,...,《）的平均值，得到目标函数对应于决策变量父 1 的确定输出值;^ ;
[0139] 第五步:通过上述第三步和第四步的提供的方法，根据目标函数输出值歹的计算 方法，计算N次，得至IjN组输入数据X1 (1 = 1，2，. . .，N)和对应的N组输出数据文（/ = 1,2,...,iV )， 将N组输入数据和对应的N组输出数据确定为目标函数确定的输入输出样本集。
[0140] 当通过以上方法确定出目标函数对应的确定的输入输出样本集后，根据该确定的 输入输出样本集，通过以下S1022步骤得到确定的带约束的多目标调度模型，以使对该确定 的带约束的多目标调度模型求解后，得到确定的、符合实际的调度要求的原料的调度信息。
[0141] 步骤S1022:对确定的输入输出数据样本集进行训练，得到确定的带约束的多目标 调度模型。
[0142] 在本发明实施例提供的原料的调度方法中，可以通过多层前向神经元网络对上述 确定的输入输出数据样本集进行训练，该多层前向神经元网络具有可以逼近连续非线性函 数的优点，从而使得训练得到的数据与实际混匀过程中需要的数据之间的误差大大减小， 提高训练后得到的数据的准确性。
[0143] 在通过多层前向神经元网络对上述确定的输入输出数据样本集进行训练时，首 先，将确定的输入输出数据样本集排列为训练数据集样本{(&丨)1 = 1，2，...，《，设置该 原料的不确定调度模型对应的多层前向神经元网络，包括设置输入神经元个数、隐层神经 元个数、输出神经元个数、预设的最小误差值;然后，对训练数据集样本，通过STA(State Transition Algorithm，状态转移算法)算法根据式（10)训练神经元网络学习过程的权重 向量，直至实际输出矩阵H(X，W)在允许误差内接近其训练数据集样本Y，从而得到极小化误 差函数Err(W)，从而逼近原料不确定调度模型中对应的目标函数，得到确定的生产成本函 数和存储成本函数;最后，根据训练后得到的确定的生产成本函数和存储成本函数得到确 定的带约束的多目标调度模型。
[0145]式（10)中，W为训练权重向量，Yi (i = 1，2，…，N)训练数据集输出样本向量，Xi (i = 1，2，…，N)训练数据集输入样本向量。
[0146] 当训练后得到的权重向量符合一定要求时，通过上述训练后得到的生产成本函数 计算出的生产成本与实际生产成本之间的误差会大大减小，以及通过上述训练后得到的存 储成本函数计算出的存储成本与实际存储成本之间的误差会大大减小，如此，对训练后得 到的确定的带约束的多目标调度模型进行求解后，得到的原料调度信息准确度很高，该原 料调度?目息具有很强的可执彳丁性，从而可以大大提尚生广效率，以及提尚广品质量。
[0147] 步骤S103:对确定的带约束的多目标调度模型进行求解，得到原料调度信息；
[0148] 在本发明实施例提供的原料的调度方法中，上述原料调度信息包括原料的配比信 息以及相应的处理措施等。其中，配比信息包括每个车间每堆的配料时间、配料种类、配料 量、每种原料中的掺配量等，处理措施规定了调度任务，以及异常调度的处理等，通过上述 原料调度信息可以对生产中原料的配比及成本进行控制。上述对确定的带约束的多目标调 度模型进行求解，得到原料调度信息，参见图4,具体可通过以下S1031-S1037步骤进行求 解。
[0149] 上述对带约束的多目标模型优化求解前，需要对多目标和多约束进行处理。首先 将上述确定的带约束的多目标调度模型中的约束条件转化为一个违约程度函数疗，具 体可通过以下步骤转化。
[0150] 第一步:将确定的带约束的多目标调度模型中的等式约束条件，通过式（11)将其 转化为不等式约束条件。
[0151] hi(X) |-μ^0,? = ρ+1,ρ+2, . . . ,q (11)
[0152] 式（11)中，|hi(X) I为确定性函数的等式约束条件，X为决策变量，μ通常取一个较 小的正数，代表等式约束能够接受的偏离范围。
[0153]第二步:将某一个体对第i个约束条件的违约程度定义为Gi(X)，该Gi(X)如式（12) 所示，
[0155] 式（12)中，Gi(X)为第i个约束条件的违约程度，gi(X)表示模型中的不等式约束，i 为第i个约束条件，第1~p个约束条件为不等式约束条件，第p+1~q个约束条件为等式约 束，X为决策变量。
[0156] 第三步:将某一个体对所有约束条件的违约程度定义为该个体违反约束条件的程 度，同时由于各约束条件的类型不同，因而具有不同的量纲和特性，需对每一约束条件的违 约程度进行标准化处理，得到标准化违约程度泛(XL该&U )如式(13)所示，
[0158] 式（13)中，q为约束条件个数，X为决策变量，Gi(X)为第i个约束条件的违约程度， Gi_max为第i个约束条件的最大违约程度。
[0159] 当得到标准化违约程度函数后，将违约程度函数)和确定的带约束的多 目标调度模型中的确定性函数进行融合，得到融合函数F(X)。该融合函数F(X)如式（14)所 示，
[0161] 式（14)中，f/ (X)为确定的生产成本函数，f21X)为确定的存储成本函数，可义)为 标准化违约程度。
[0162] 得到融合函数F(X)后，通过基于带变异的CM0STA(约束多目标状态转移算法)算法 对该融合函数求解，得到原料调度信息。
[0163] 步骤S1031:初始化设置CM0STA算法的相关参数、Pareto最优搜索解集种群Ps、 Pareto非劣解集种群Pb、最大迭代次数Gmax，最大个体违约程度值ε以及水平截集pc。
[0164] 上述设置CM0STA算法的相关参数，包括设置STA的搜索种群参数、扩展因子α、转移 因子β、扩展因子γ和平移因子S的值;上述设置Pareto最优搜索解集种群Ps以及Pareto非 劣解集种群Pb，包括设置种群个数SE、随机初始化并均匀产生数目为num存放Pareto最优搜 索解集种群Ps，同时初始化搜索中间种群为空，且设置Pareto非劣解集种群P B = Ps;设置水 平截集pc>〇。
[0165] 步骤S1032:从当前种群出发，开始当前迭代操作，对带约束的多目标模型采用 CM0STA搜索策略，得到新的种群。
[0166] 当设置好上述参数后，具体可通过执行以下S10321-S10323步骤得到新的种群，如 图5所示。
[0167] 步骤S10321:对种群Ps中的每个个体执行状态转移操作，得到中间种群PM1。
[0168] 在本发明实施例中，在对种群Ps中的每个个体执行状态转移操作，时，即依次进行 扩展操作、旋转操作和转移操作，将三次操作中非劣等级排序为1的个体放入中间种群Pmi 中，得到中间种群Pmi。
[0169] 步骤S10322:对种群PB和得到的中间种群?(《进行融合，以及进行Pareto非劣等级 排序，得到新的种群Pm。
[0170] 当将种群Ps中的所有个体都进行状态转移操作后，得到中间种群PM1，对PM1和Pb进 行融合，形成初始的新种群Pm;利用训练好的融合函数F(X)检验初始的新种群Pm的可行性， 将得到的可行解进行Pareto非劣等级排序后放入中间种群，将不可行解进行违约程度 升序排序后放入中间种群Pj中，将排序后的Pi和Pj融合得到新种群Pm，其中，Pm= [Pi，Pj]。
[0171] 在本发明实施例中，更新种群参数时，利用训练好的神经元网络检验搜索种群的 可行性。其中，对Pm中&(.〇 = 0的个体，即为可行解，用Pareto非劣等级排序可行解PI;对PM中 ?(.Λ^>?的个体，即为不可行解，用违约程度排序不可行解巧。此时，在F(X)所得的可行解种 群PdP不可行解种群Pj个体中，最优解的选取采取以下准则进行选取：
[0172] 当两个个体1和1'都可行时，将适应值F(X)比较小的个体作为最优个体；当两个 个体f和都不可行时，将违约程度值比较小的个体作为最优个体；当f可行且 不可行时，如果浮@，且，将f'个体作为最优个体，如果或 ，将i个体作为最优个体。
[0173] 步骤S10323:根据新的种群Pm对种群参数进行更新。
[0174]更新种群参数时，主要更新非劣解集种群PB、最优解搜索集种群Ps以及最大个体违 约程度值ε。在对种群Ps和Pb进行赋值时，将Pm中排在前num个个体赋予Ps，即Ps = Pm( 1 :num); 种群Pm中非劣等级为1的个体集Pmi放入到Pb中，即Pb = Pmi。同时对最大个体违约程度值ε做 自适应调整。
[0175]在本发明实施例中，为了将不可行解的比例保持在一个固定的水平pc，根据得到 中间种群Pm以及自适应调整公式对最大个体违约程度值ε自适应调整，其中，自适应调整公 式如下式(15)所示，
[0177] 式（15)中，pk为每进化Κ代后不可行解在群体中所占比例的平均值，pc为水平截集， ε为最大个体违约程度值，Other为除pk<p。以及除pk>p。外的情况，ε '为自适应调整后的最 大个体违约程度值。
[0178] 上述pk可通过以下式（16)计算出，
[0180] 式(16)中，Ni表示第i代群体中不可行解数量;num是迭代搜索中个体数量。
[0181] 步骤S1033:判断新的种群是否陷入局部最优，若陷入局部最优，则执行以下S1034 步骤，若没有陷入局部最优，则执行以下S1035步骤，其中，局部最优为当前的Pareto非劣解 集种群P s与上一轮迭代得到的Pareto非劣解集种群Ps相等时的状态。
[0182]在本发明实施例中，首先判断迭代次数k是否为50的整数倍且Ps中的第k代Ps和第 k-Ι代Ps是否相等，若是，则说明迭代陷入局部最优。
[0183] 步骤S1034:执行全局更新策略，对种群采用变异操作进行种群更新，即执行最优 解更新公式，得到变异后的父代中种群Ps。上述最优解更新公式如下（17)式所示，
[0184] xk = xk+rand(0,1) X (xbest-xk-1) (17)
[0185] 式（17)中，xk为第k代种群，xh为第k-1代种群，rand(0，l)为在(0，1)范围内的随 机数， Xbe3St为第k-Ι代种群中的最优种群。
[0186]步骤S1035:判断当前迭代次数是否达到最大迭代次数Gmax，若没达到，将所有中 间种群置空，则执行以上步骤S1032，对CM0STA搜索中得到的Pareto最优搜索解集种群Ps和 Pareto非劣解集种群P B执行新一轮迭代操作，即重新搜索新种群并判断是否需要种群更 新;若达到，则执行以下S1036步骤。
[0187]步骤S1036:输出Pareto非劣解集种群PB作为搜索结果的非劣最优解集；
[0188]步骤S1037:通过T0PSIS法从搜索结果输出的非劣最优解集中选取一组非劣解，将 选取的一组非劣解确定为原料调度信息。
[0189] 步骤S104:根据原料调度信息执行对混匀过程中原料的调度。
[0190] 上述确定出的原料调度信息与实际中原料的调度之间的误差很小，原料调度信息 的准确度很高，根据该原料调度信息执行原料的调度，可以大大提高生产效率，以及提高产 品的质量。
[0191] 本发明实施例提供的原料的调度方法，通过建立的原料不确定调度模型可以准确 地描述实际生产过程，通过确定化和带变异的CM0STA算法可以求取准确度很高的原料调度 信息，原料调度信息可执行性很强，根据该原料调度信息执行原料的调度，可以大大提高生 广效率，以及提尚广品的质量。
[0192] 图6是本发明一种实施方式的一种原料的调度装置的示意图，参照图6,该装置包 括：
[0193] 建立模块S1，用于建立原料不确定调度模型；
[0194] 处理模块S2,用于对原料不确定调度模型进行确定化处理，得到确定的带约束的 多目标调度模型；
[0195] 计算模块S3,用于对确定的带约束的多目标调度模型进行求解，得到原料调度信 息；
[0196] 执行模块S4,用于根据原料调度信息执行对混匀过程中原料的调度。
[0197] 在本发明实施例中，由于混匀过程中存在来料时间的不确定性，因此，引入三角模 糊数来描述原料场混匀过程堆计划原料不确定调度模型中的不确定参数，从而建立以各混 匀堆的最小生产成本和/或最小存储成本为目标的原料的不确定调度模型。
[0198] 在本发明实施例中，上述原料可以是炼铁所需的原料，建立上述原料不确定模型， 可以真实反映出出炼铁所需的原料种类、品位、库存量、来料时间及配料时间的范围，其中， 上述建立模块S1在建立上述原料不确定调度模型时，具体可通过以下获取单元、第一建立 单元、第一确定单元和第二确定单元来建立。
[0199] 上述建立模块S1包括获取单元、第一建立单元、第一确定单元、第二确定单元和第 二建立单元。
[0200] 获取单元，用于获取混匀过程的配料基本信息；
[0201] 上述获取单元获取混匀过程的配料基本信息时，首先针对原料场混匀过程的实际 问题，分析现场工况，其中包括分析现场混匀所需的原料种类，计划内混匀料堆数，每种原 料各个化学成分含量，每种原料的采购价格，每种原料计划内进料量，计划内混匀料堆中各 个化学成分含量指标，根据分析结果建立原料分析表，该原料分析表与上述实施例中提供 的一种原料调度方法中的原料分析表相同，在此不再赘述。获取单元建立上述原料分析表 后，可以通过该原料分析表及实际生产的需求，确定出计划内购买原料的种类、每种原料的 成分指标、每种原料的购买量、每种原料的来料时间，每种原料的不确定配料时间、混勾料 堆的每种物料含量范围、各个混匀料堆应堆重范围、每堆物料特殊物料使用信息及堆料场 库存信息，生成混匀过程的配料基本信息。
[0202] 上述计划可以是周计划、月计划、季度计划或年计划，具体可根据实际生产的需求 确定，为了便于清楚地说明本发明实施例提供的数据，无特殊说明的情况下，上述计划均指 月计划。
[0203] 在原料的调度过程中，异常调度原料都会对生产过程产生影响，以至带来额外的 损失，因此，需要在调度原料时考虑异常调度的情况。其中，异常调度主要包括提前调度或 拖延调度。
[0204] 第一建立单元，用于采用三角模糊数建立异常调度对应的惩罚函数；
[0205] 上述第一建立单元可通过上述实施例一种炼铁原料的调度方法中提供的惩罚函 数的建立方法，建立异常调度对应的惩罚函数，在此不再赘述。
[0206] 在保证炼铁效率以及钢铁质量的前提下，控制炼铁成本也是制定原料调度方案所 必需考虑的关键指标。根据炼铁厂对成本控制的要求，通过以下第一确定单元可以确定出 炼铁过程的目标函数。
[0207] 第一确定单元，根据所述配料基本信息和所述异常调度对应的惩罚函数，确定混 匀过程中的目标函数。
[0208] 上述目标函数主要指调度物料所花费的成本控制函数，在本发明中，优选地，上述 目标函数主要指调度物料所花费的生成成本函数和/或存储成本函数。
[0209] 在实际混匀过程中，由于各种原料参数受到生产条件的影响，因此，在建立原料的 调度方案时需要考虑各种原料的约束条件。
[0210] 第二确定单元，用于确定混匀过程中各种原料的约束条件；
[0211 ]上述第二确定单元，在确定混勾过程中各种原料的约束条件时，需要确定出混勾 过程中的原料的成分约束条件、堆重约束条件、库存约束条件及特殊原料约束条件，每种约 束条件的具体确定过程，可通过上述一种原料的调度方法中提供的约束条件的确定方法来 确定，在此不再赘述。
[0212] 在本发明实施例提供的原料的调度装置中，通过上述原料的成分约束条件、上述 堆重约束条件、上述库存约束条件及上述特殊原料约束条件，可以真实反映出原料场各堆 物料的堆重范围，每堆料堆中各个化学成分的比例范围以及每种原料的库存可用量，如此， 在生产过程中，可以根据各个原料的约束条件及时填充原料，避免料堆无料、出现生产中断 的情况，从而可以提尚广品的生广效率和广品质量。
[0213] 第二建立单元，用于根据所述目标函数和所述约束条件建立原料不确定调度模 型，原料不确定调度模型如下式(1)所示，
[0217] 式（1)中，/pQ为异常调度对应的惩罚函数;xu为决策变量，代表第j堆中所需第 i种原料的相应下料量;C1;1为第i中原料的采购价格;&(x)为生产成本，代表完成堆计划所 需原料的总花费，fi(X)为目标函数，代表最小的生产成本;m表示原料数;η表示堆数；
[0218] C1>2为第i种原料每吨每天所需的保存费用，爲为额外的存储成本，C2(x)为总存储 成本，代表完成堆计划期间所有原料存储所消耗的费用，f 2(X)为目标函数，代表最小的总 存储成本；
[0219] s.t.为各种原料的约束条件，Qy为所配第j堆混合料堆中第τ种化学成分的含量 指标，Ay为计划第j堆混合料堆中第τ种化学成分的含量指标，qy为计划第j堆混合料堆中 第τ种化学成分含量指标的最大波动范围， Sl(t)为第i种原料当前的库存量，h为配比第j堆 所堆堆重，D」为计划中第j堆的标准堆重，为计划中相对应仏指标的最大波动范围， X1>J为 特殊物料在j堆中掺配量;为第j堆所需第i种原料相应的下料量。
[0220] 在实际炼铁过程中，由于炼铁所需的含铁原料来料时间存在不确定性，导致根据 上述原料不确定调度模型难以进行原料的调度，因此，需要通过上述处理模块S2确定出原 料的调度信息，以使根据该原料的调度信息可以准确执行炼铁过程中原料的调度。
[0221] 上述处理模块S2包括处理单元和训练单元。
[0222] 处理单元，用于通过模糊模拟技术对原料不确定调度模型中每个含有模糊函数方 程的目标函数进行处理，得到目标函数确定的输入输出样本集。
[0223] 在本发明实施例中，由于混匀过程中存在来料时间的不确定性，因此，上述处理单 元在确定输入值后，通过三角模糊数的隶属度函数对不确定参数进行抽样，并根据基于抽 样得到的抽样值，得到目标函数确定的输出样本值。具体的确定过程与上述实施例一种炼 铁原料的调度方法中提供的确定过程类同，在此不再赘述。
[0224]当通过上述处理单元处理得到目标函数对应的确定的输入输出样本集后，根据该 确定的输入输出样本集，通过以下训练单元得到确定的带约束的多目标调度模型，以使对 该确定的带约束的多目标调度模型求解后，得到确定的、符合实际的调度要求的原料的调 度信息。
[0225] 训练单元，用于对确定的输入输出数据样本集进行训练，得到确定的带约束的多 目标调度模型。
[0226] 上述训练单元，可通过上述实施例一种原料的调度方法中提供的训练方法，对确 定的输入输出数据样本集进行训练，具体的训练过程在此不再赘述。
[0227] 当上述处理模块确定出确定的带约束的多目标调度模型后，可通过以下计算模块 得到原料调度信息。
[0228] 上述计算模块S3包括设置单元、搜索单元、更新单元、迭代判断单元和选取单元。 [0229]设置单元，用于初始化设置CM0STA算法的相关参数、Pareto最优搜索解集种群P s、 Pareto非劣解集种群Pb、最大迭代次数Gmax，最大个体违约程度值ε以及水平截集pc;
[0230] 搜索单元，用于从当前种群出发，开始当前迭代操作，对带约束的多目标优化模型 采用CM0STA搜索策略，得到新的种群；
[0231] 更新单元，用于判断新的种群是否陷入局部最优，若是，对种群采用变异操作进行 种群更新，局部最优为当前的Pareto非劣解集种群P s与上一轮迭代得到的Pareto非劣解集 种群Ps相等时的状态；
[0232] 迭代判断单元，用于判断搜索单元和更新单元是否需要迭代进行，若当前迭代次 数没达到最大迭代次数Gmax，则重新启动搜索单元和更新单元;若达到，则迭代过程终止， 输出Pareto非劣解集种群Pb作为搜索结果的非劣最优解集；
[0233] 选取单元，用于通过T0PSIS(Technique for Order Preference by Similarity to an Ideal Solution，逼近理想解排序法)法从搜索结果输出的非劣最优解集中选取一 组非劣解，将选取的一组非劣解确定为原料调度信息。
[0234] 在本发明实施例中，上述设置单元可通过上述实施例一种原料的调度方法中提供 的设置方法设置上述各种参数;上述搜索单元，可通过上述实施例一种原料的调度方法中 提供的搜索方法，得到新的种群;上述更新单元，可通过上述实施例一种原料的调度方法中 提供的判断准则，判断是否需要用变异操作对种群更新;上述迭代判断单元，可通过上述实 施例一种原料的调度方法中提供的迭代终止判断准则，判断是否迭代终止;上述选取单元， 可通过上述实施例一种原料的调度方法中提供的选取方法，从搜索结果输出的非劣最优解 集中选取一组非劣解，将选取的一组非劣解确定为原料调度信息，在此均不再赘述。
[0235] 当上述选取单元将选取的一组非劣解确定为原料调度信息，该原料调度信息与实 际中原料的调度之间的误差很小，原料调度信息的准确度很高，上述执行模块S4根据该原 料调度?目息执彳丁原料的调度，可以大大提尚生广效率，以及提尚广品的质量。
[0236] 综上所述，本实施例提供了一种原料的调度方法及装置，通过建立的原料不确定 调度模型可以准确地描述实际生产过程，通过确定化和带变异的CM0STA算法可以求取准确 度很高的原料调度信息，原料调度信息可执行性很强，根据该原料调度信息执行原料的调 度，可以大大提尚生广效率，以及提尚广品的质量。
[0237] 本发明的实施例是为了示例和描述起见而给出的，而并不是无遗漏的或者将本发 明限于所公开的形式。很多修改和变化对于本领域的普通技术人员而言是显而易见的。选 择和描述实施例是为了更好说明本发明的原理和实际应用，并且使本领域的普通技术人员 能够理解本发明从而设计适于特定用途的带有各种修改的各种实施例。
【主权项】
1. 一种原料的调度方法，其特征在于，所述方法包括： 建立原料不确定调度模型； 对所述原料不确定调度模型进行确定化处理，得到确定的带约束的多目标调度模型； 对所述确定的带约束的多目标调度模型进行求解，得到原料调度信息； 根据所述原料调度信息执行对混匀过程中原料的调度。2. 根据权利要求1所述方法，其特征在于，所述建立原料不确定调度模型，包括： 获取混匀过程的配料基本信息； 采用Ξ角模糊数，建立异常调度对应的惩罚函数； 根据所述配料基本信息和所述异常调度对应的惩罚函数，确定混匀过程中的目标函 数； 确定混匀过程中各种原料的约束条件； 根据所述目标函数和所述约束条件建立原料不确定调度模型，所述原料不确定调度模 型如下式(1)所示，式(1)中，./'(.V)为异常调度对应的惩罚函数;XI,功决策变量，代表第j堆中所需第i种原 料的相应下料量;Cl,1为第i中原料的采购价格;Ci(x)为生产成本，代表完成堆计划所需原 料的总花费;fi(X)为目标函数，代表最小的生产成本;m表示原料数;η表示堆数； Cl,2为第i种原料每吨每天所需的保存费用，4为额外的存储成本，C2(X)为总存储成本， 代表完成堆计划期间所有原料存储所消耗的费用，f2(X)为目标函数，代表最小的总存储成 本； s.t.为各种原料的约束条件，为所配第j堆混合料堆中第τ种化学成分的含量指标， At,j为计划第j堆混合料堆中第τ种化学成分的含量指标，为计划第j堆混合料堆中第τ种 化学成分含量指标的最大波动范围，Si(t)为第i种原料当前的库存量，bj为配比第j堆所堆 堆重，Dj为计划中第j堆的标准堆重，为计划中相对应町指标的最大波动范围，XU为特殊 物料在j堆中渗配量;XI,J为第j堆所需第i种原料相应的下料量。3. 根据权利要求1所述方法，其特征在于，所述对所述原料不确定调度模型进行确定化 处理，得到确定的带约束的多目标调度模型，包括： 通过模糊模拟技术对所述原料不确定调度模型中每个含有模糊函数方程的目标函数 进行处理，得到所述目标函数确定的输入输出样本集； 对所述确定的输入输出数据样本集进行训练，得到确定的带约束的多目标调度模型。4. 根据权利要求3所述方法，其特征在于，所述通过模糊模拟技术对所述原料不确定调 度模型中每个含有模糊函数方程的目标函数进行处理，得到所述目标函数对应的确定的输 入输出样本集，包括： 将每个目标函数中的惩罚函数表示为预设的模糊函数方程表示式，所述预设的模糊函 数方程表示式如下式(2)所示， (X，浸） （2) 式(2)中，夫为目标函数中惩罚函数对应的模糊函数，表示，与乂和易的函数 关系，X为一组决策变量，表示原料的下料量，#为Ξ角模糊数表示的函数参数； 获取Ξ角模糊数的隶属度函数； 确定目标函数的一组决策变量χ?，对所述Ξ角模糊数在所述Ξ角模糊数的隶属度函数 内抽样，得到I的一系列抽样值线二1，2，...乂），Κ为抽样数； 根据所述决策变量Χ?和每个抽样值边，通过公式（2)得到Κ个确定的输出 乂（& = ]，2,...，/〇，并求取1(个乂（&二1，2，".，。的平均值，得到目标函数对应于决策变量乂1 的确定输出值户； 根据所述目标函数输出值V的计算方法，计算Ν次，得到腺且输入数据χ?(1 = 1，2,.. .，Ν) 和对应的Ν组输出数据户(/ = 1，2，,..，Λ〇,将所述Ν组输入数据和对应的所述Ν组输出数据确 定为目标函数确定的输入输出样本集。5. 根据权利要求3所述的方法，其特征在于，所述对所述确定的输入输出数据样本集进 行训练，得到确定的带约束的多目标模型，包括： 将所述确定的输入输出数据样本集排列为训练数据集样本； 对所述训练数据集样本，通过状态转移算法STA训练神经元网络学习过程的权重向量， 逼近所述原料不确定调度模型中对应的目标函数，得到确定的生产成本函数和/存储成本 函数； 根据训练后得到的所述确定的生产成本函数和存储成本函数建立确定的带约束的多 目标调度模型。6. 根据权利要求1所述的方法，其特征在于，所述对所述确定的带约束的多目标调度模 型进行求解，得到原料调度信息，包括： 初始化设置CMOSTA算法的相关参数、Pareto最优捜索解集种群Ps、Pareto非劣解集种群 时、最大迭代次数Gmax，最大个体违约程度值εΚ及水平截集Pc; 从当前种群出发，开始当前迭代操作，对所述确定的带约束的多目标调度模型采用 CMOSTA捜索策略，得到新的种群； 判断所述新的种群是否陷入局部最优，若是，对所述种群采用变异操作进行种群更新， 所述局部最优为当前的化reto非劣解集种群Ps与上一轮迭代得到的化reto非劣解集种群Ps 相等时的状态； 判断当前迭代次数是否达到最大迭代次数Gmax，若没达到，将所有中间种群置空，对 CMOSTA捜索中得到的化reto最优捜索解集种群Ps和化reto非劣解集种群Pb执行新一轮迭代 操作，即重新捜索新种群并判断是否需要种群更新;若达到，则输出化reto非劣解集种群Pb 作为捜索结果的非劣最优解集； 通过TOPSIS(Technique for Order Preference by Similarity to an Ideal Solution,逼近理想解排序法)从所述捜索结果输出的非劣最优解集中选取一组非劣解，将 选取的所述一组非劣解确定为原料调度信息。7. -种原料的调度装置，其特征在于，所述装置包括： 建立模块，用于建立原料不确定调度模型； 处理模块，用于对所述原料不确定调度模型进行确定化处理，得到确定的带约束的多 目标调度模型； 计算模块，用于对所述确定的带约束的多目标调度模型进行求解，得到原料调度信息； 执行模块，用于根据所述原料调度信息执行对混匀过程中原料的调度。8. 根据权利要求7所述装置，其特征在于，所述建立模块包括： 获取单元，用于获取混匀过程的配料基本信息； 第一建立单元，用于采用Ξ角模糊数，建立异常调度对应的惩罚函数； 第一确定单元，根据所述配料基本信息和所述异常调度对应的惩罚函数，确定混匀过 程中的目标函数； 第二确定单元，用于确定混匀过程中各种原料的约束条件； 第二建立单元，用于根据所述目标函数和所述约束条件建立原料不确定调度模型，所 述原料不确定调度模型如下式(1)所示，式（1)中，ΛΧ)为异常调度对应的惩罚函数;xiu为决策变量，代表第j堆中所需第i种 原料的相应下料量;Cl,1为第i中原料的采购价格;Ci(x)为生产成本，代表完成堆计划所需 原料的总花费，fi(X)为目标函数，代表最小的生产成本;m表示原料数;η表示堆数； Cl,2为第i种原料每吨每天所需的保存费用，巧为额外的存储成本，C2(x)为总存储成本， 代表完成堆计划期间所有原料存储所消耗的费用，f2(X)为目标函数，代表最小的总存储成 本； s.t.为各种原料的约束条件，为所配第j堆混合料堆中第τ种化学成分的含量指标， At,j为计划第j堆混合料堆中第τ种化学成分的含量指标，为计划第j堆混合料堆中第τ种 化学成分含量指标的最大波动范围，Si(t)为第i种原料当前的库存量，bj为配比第j堆所堆 堆重，Dj为计划中第j堆的标准堆重，为计划中相对应町指标的最大波动范围，XU为特殊 物料在j堆中渗配量;XI,J为第j堆所需第i种原料相应的下料量。9. 根据权利要求7所述装置，其特征在于，所述处理模块包括： 处理单元，用于通过模糊模拟技术对所述原料不确定调度模型中每个含有模糊函数方 程的目标函数进行处理，得到所述目标函数确定的输入输出样本集； 训练单元，用于对所述确定的输入输出数据样本集进行训练，得到确定的带约束的多 目标调度模型。10. 根据权利要求7所述的装置，其特征在于，所述计算模块包括： 设置单元，用于初始化设置CMOSTA算法的相关参数、Pareto最优捜索解集种群Ps、 化reto非劣解集种群Pb、最大迭代次数Gmax，最大个体违约程度值εΚ及水平截集Pc; 捜索单元，用于从当前种群出发，开始当前迭代操作，对所述确定的带约束的多目标模 型采用CMOSTA捜索策略，得到新的种群； 更新单元，用于判断所述新的种群是否陷入局部最优，若是，对所述种群采用变异操作 进行种群更新，所述局部最优为当前的化reto非劣解集种群Ps与上一轮迭代得到的化reto 非劣解集种群Ps相等时的状态； 迭代判断单元，用于判断捜索单元和更新单元是否需要迭代进行，若当前迭代次数没 达到最大迭代次数Gmax，则重新启动捜索单元和更新单元;若达到，则迭代过程终止，输出 化reto非劣解集种群Pb作为捜索结果的非劣最优解集； 选取单元，用于通过逼近理想解排序法T0PSIS法从所述捜索结果输出的非劣最优解集 中选取一组非劣解，将选取的所述一组非劣解确定为原料调度信息。
【文档编号】G05B19/418GK106094751SQ201610453502
【公开日】2016年11月9日
【申请日】2016年6月21日
【发明人】王雅琳, 孙克楠, 何海明, 陈晓方, 曹跃, 孙备, 桂卫华
【申请人】中南大学