内控复核型关联控制的混合类快速成型架构方法

摘要

本发明公开了一种内控复核型关联控制的混合类快速成型架构方法，利用深度学习策略构造初始控制集合体以保证控制热源体均匀分布在控制空间中；再利用复核型关联控制对最优控制热源体进行局部控制，寻找最优解。本发明在进化过程中只对精英控制热源体进行单纯形确定性局部控制，能较大减少算法适应度函数计算次数；深度学习策略能使控制热源体均匀分布在控制空间中，为全局控制奠定基础；本发明还通过柯西变异算子以避免改方法陷入局部最优。最终通过最短时间，得到全局最优解。

权利要求

1.内控复核型关联控制的混合类快速成型架构方法，其特征在于，利用深度学习策略构造初始控制集合体以保证控制热源体均匀分布在控制空间中；再利用复核型关联控制对最优控制热源体进行局部控制，寻找全局最优解；包括如下步骤：(1)在D维控制空间中随机产生N个候选控制热源体，构造一初始控制集合体RP；

(2)利用深度学习策略，产生所述初始控制集合体RP的对应反向控制集合体OP；

(3)合并所述初始控制集合体RP和所述反向控制集合体OP，将合并后的2N个控制热源体按照其对应的优化控制度函数大小进行升序排序，选取前N个控制热源体作为类快速成型架构方法的初始控制集合体；

(4)计算所述初始控制集合体中每个带电的控制热源体所受作用力，判断是否满足结束条件，若满足，则结束，否则进行下一步；

(5)将所述初始控制集合体中的控制热源体分为精英控制热源体和一般控制热源体，对精英控制热源体进行单纯形控制；

(6)移动所述控制热源体，更新群体中每个所述控制热源体的位置；

(7)当所有所述控制热源体稳定控制到一定程度时，对所述精英控制热源体和所述一般控制热源体进行柯西变异的局部控制，得到全局最优解；

其中，所述一定程度是指控制热源体群中控制热源体的组合，使控制热源体对应的优化控制度函数值在允许误差范围内。

2.如权利要求1所述的混合类快速成型架构方法，其特征在于，步骤(4)所述结束条件包括步骤结束或者优化控制度函数值满足设定值。

3.如权利要求1所述的混合类快速成型架构方法，其特征在于，所述精英控制热源体是指优化控制度函数值小的前N1个控制热源体；N-N1个控制热源体则为一般控制热源体。

说明书

内控复核型关联控制的混合类快速成型架构方法

技术领域

[0001] 本发明涉及计算机领域，尤其涉及内控复核型关联控制的混合类快速成 型架构方法。

背景技术

[0002] 类快速成型架构(CRPA)法是Birbil和Fang于2003年提出的一种新的全 局优化启发式方法。该方法模拟快速成型场中的吸引和排斥架构，将问题的 解看作一个带电控制热源体，通过带电控制热源体间吸引-排斥架构准则产生 群体智能指导优化控制。由于该方法具有良好的控制性能，因此在函数优化、 项目调度、神经网络训练和旅行商问题等领域中得到广泛的应用。但其存在 控制集合体多样性差、对初始控制集合体的分布要求较高、易出现早熟稳定 控制和局部控制能力差等缺点。且该方法中控制热源体所受合力会忽略解空 间中的某些可行控制区域，在一些复杂优化问题中容易出现早熟稳定控制。 为了避免此问题，Birbil和Fang在标准CRPA方法的基础上，对合力计算步 骤进行了修改，提出一种修改的类快速成型架构方法，在计算合力时，对当 前群体中离最优控制热源体xbest最远的个体xp施加扰动，使其可能移动到被忽 略区域，增强算法的全局控制能力。但其控制速度慢，全局的控制能力有限。

发明内容

[0003] 针对上述现有技术存在的缺陷，本发明提供一种内控复核型关联控制的 混合类快速成型架构方法，增强了算法在最优点周围的局部控制能力，以加 快稳定控制速度，从而提高了全局的控制能力。

[0004] 本发明提供的内控复核型关联控制的混合类快速成型架构方法，其改进 之处在于，利用深度学习策略构造初始控制集合体以保证控制热源体均匀分 布在控制空间中；再利用复核型关联控制对最优控制热源体进行局部控制， 寻找全局最优解；包括如下步骤：

[0005] (1)在D维控制空间中随机产生N个候选控制热源体，构造一始控制集 合体RP；

[0006] (2)利用深度学习策略，产生所述初始控制集合体RP的对应反向控制集 合体OP；

[0007] (3)合并所述一始控制集合体RP和所述反向控制集合体OP，将合并后 的2N个控制热源体按照其对应的优化控制度函数大小进行升序排序，选取前N个控制热源体作为类快速成型架构方法的初始控制集合体；

[0008] (4)计算所述初始控制集合体中每个带电的控制热源体所受作用力，判 断是否满足结束条件，若满足，则结束，否则进行下一步；

[0009] (5)将所述初始控制集合体中的控制热源体分为精英控制热源体和一般 控制热源体，对精英控制热源体进行单纯形控制；

[0010] (6)移动所述控制热源体，更新群体中每个所述控制热源体的位置；

[0011] (7)当所有所述控制热源体稳定控制到一定程度时，对所述精英控制热 源体和所述一般控制热源体进行柯西变异的局部控制，得到全局最优解。

[0012] 优选的，在D维控制空间中：

[0013] 设任一控制热源体为x∈[l,u]，则：

[0014] 其反向控制热源体x'为：x'＝l+u-x；

[0015] 将所述反向控制热源体x'扩展到D维控制空间，设P为D维控制空间的一 个候选解，则：

[0016] P＝(x1,x2,…,xD)；

[0017] 其中：xi∈[ai,bi],i＝1,2,…,D；

[0018] 所述候选解P对应的反向解为：

[0019] P'＝(x1',x'2,…,x'D)；

[0020] 其中：xi'＝li+ui-xi；

[0021] 在所述一始控制集合体RP中，第i个控制热源体的位置为：

[0022] Xi＝(xi1,xi2,…,xiD),xij∈[li,uj],i＝1,2,…,N,j＝1,2,…,D；

[0023] 在所述反向控制集合体OP中，第i个控制热源体的位置为：

[0024] Xi'＝(xi'1,xi'2,…,xi'D),i＝1,2,…,N；

[0025] 在所述反向控制集合体OP中，第i个控制热源体的第j维分量为：

[0026] xi'j＝lj+uj-xij。

[0027] 较优选的，采用最小化方法进行所述2N个控制热源体的升序排序。

[0028] 较优选的，步骤(4)计算所述初始控制集合体中每个带电控制热源体所 受作用力，包括：

[0029] 计算控制热源体xi的电荷量qi的公式为：

[0030]

[0031] 根据库仑定律，对控制集合体中任意两个控制热源体xi和xj，计算它们之 间的作用力Fij为：

[0032]

[0033] 其中i,j＝1,2,…,N；每个控制热源体xi所受到的作用力Fi为：

[0034]

[0035] 较优选的，步骤(4)所述结束条件包括步骤结束或者优化控制度函数值 满足设定值。

[0036] 较优选的，所述精英控制热源体是指优化控制度函数值小的前N1个控制 热源体；N-N1个控制热源体则为一般控制热源体。

[0037] 较优选的，所述对精英控制热源体进行单纯形控制的步骤包括：

[0038] 1)给定初始单纯形参数：

[0039] 顶点的控制热源体xi∈Rn，i＝1,2,…,n+1；

[0040] 反射系数α，α>0；

[0041] 扩展系数γ，γ>1；

[0042] 压缩系数β，β∈(0,1)；

[0043] 允许误差ε，ε>0；

[0044] 计算顶点的优化控制度函数 i＝1,2,…,n+1，置k＝0；

[0045] 2)确定位置最高点的控制热源体k次方为 位置次高点的控制热源体k 次方为位置最低点的控制热源体k次方为 确定优化控制度函数为：

[0046]

[0047]

[0048]

[0049] 计算n个点的形心 并计算形心处的优化控制度函数 其中：

[0050]

[0051] 3)进行反射，令优化控制率r型的控制热源体k次方为：

[0052]

[0053] 计算优化控制率r型的优化控制度函数

[0054] 4)判断优化控制率r型的优化控制度函数值若小于位置最低点的优化控 制度函数值，即 则进行扩展，令优化控制率e型的控制热源体 为：

[0055]

[0056] 计算优化控制率e型的优化控制度函数 且进入步骤5)；

[0057] 判断优化控制率r型的优化控制度函数值若不小于位置最低点的优化控 制度函数值且不大于位置次高点的优化控制度函数值，即 则置位置最高点的控制热源体k次方等于优化控制率r型的控制热源体k次 方，对应的优化控制度函数值也相等，即 且进入步骤 7)；

[0058] 判断优化控制率r型的优化控制度函数值若大于位置最高点的优化控制 度函数值大小，即 则进行压缩，令：

[0059]

[0060] 其中h′∈{h,r}， 为位置最高点中不易控制点的优化控制度函数；且令优化 控制率c型的控制热源体k次方为：

[0061]

[0062] 计算优化控制率c型的优化控制度函数 且进入步骤6)；

[0063] 5)判断优化控制率e型的优化控制度函数值若小于优化控制率r型的优化 控制度函数值，即 则令位置最高点的控制热源体k次方等于优 化控制率e型的控制热源体k次方，对应的优化控制度函数值也相等，即  且进入步骤6)；否则，令位置最高点的控制热源体k 次方等于优化控制率r型的控制热源体k次方，对应的优化控制度函数值也相 等，即 且进入步骤7)；

[0064] 6)判断优化控制率c型的优化控制度函数值若不大于位置最高点中不易控 制点的优化控制度函数值，即 则令位置最高点的控制热源体k 次方等于优化控制率c型的控 制热 源体k次方 ，对应的优化 控制度函数值也相  等 ，即且进入步骤7)；否则进行收缩，令：

[0065]

[0066] 其中i＝1,2,…,n+1，计算顶点的优化控制度函数 再进入步骤7)；

[0067] 7)判断顶点的优化控制度函数值 和形心处控制热源体的优化控制度 函数值 是否满足稳定控制准则，若：

[0068]

[0069] 则停止计算，且xi视为最优解；否则，令k的次数加1，即k＝k+1，并返回步 骤2)。

[0070] 较优选的，步骤(6)移动控制热源体xi后，其移动位置公式为：

[0071]

[0072] 移动的控制热源体xi″将沿着作用力的方向移动。

[0073] 较优选的，步骤(7)中一定程度是指控制热源体群中控制热源体的组合， 使其对应的优化控制度函数值在允许误差范围内。

[0074] 本发明的技术方案中，在进化过程中只对精英控制热源体进行单纯形确 定性局部控制，能较大减少算法适应度函数计算次数；深度学习策略能使控 制热源体均匀分布在控制空间中，为全局控制奠定基础；本发明还通过柯西 变异算子以避免改方法陷入局部最优。最终通过最短时间，得到全局最优解。

附图说明

[0075] 图1为本发明实施例的流程图；

[0076] 图2为本发明实施例的三种方法对四个函数的优化操控结果比较；

[0077] 图3为本发明实施例的Sphere函数的优化操控稳定控制曲线；

[0078] 图4为本发明实施例的Rosenbrock函数的优化操控稳定控制曲线；

[0079] 图5为本发明实施例的Griewank函数的优化操控稳定控制曲线；

[0080] 图6为本发明实施例的Rastrigin函数的优化操控稳定控制曲线；

[0081] 图7为本发明实施例的四种方法对4个函数的优化操控结果比较。

具体实施方式

[0082] 为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举 出优选实施例，对本发明进一步详细说明。然而，需要说明的是，说明书中 列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的 理解，即便没有这些特定的细节也可以实现本发明的这些方面。

[0083] 本实施例提出的一种内控复核型关联控制的混合类快速成型架构方法， 其流程图如图1所示，利用深度学习策略构造初始控制集合体以保证控制热源 体均匀分布在控制空间中；再利用复核型关联控制对最优控制热源体进行局 部控制，寻找最优解。具体包括如下步骤：

[0084] (1)在D维控制空间中随机产生N个候选控制热源体，构造一始控制集 合体RP；

[0085] 在D维控制空间中：

[0086] 设任一控制热源体为x∈[l,u](l为拉普拉斯空间，u为热源空间)，将x 映射到上述两个复合型空间中进行计算。则：

[0087] 其反向控制热源体x'为：x'＝l+u-x；

[0088] 将所述反向控制热源体x'扩展到D维控制空间，设P为D维控制空间的一 个候选解，则：

[0089] P＝(x1,x2,…,xD)；

[0090] 其中：xi∈[ai,bi],i＝1,2,…,D；ai，bi为受控对象的数字化控制向量；D为 向量中的整数；

[0091] 所述候选解P对应的反向解为：

[0092] P'＝(x1',x'2,…,x'D)；

[0093] 其中：xi'＝li+ui-xi；

[0094] 在一始控制集合体RP中，第i个控制热源体的位置为：

[0095] Xi＝(xi1,xi2,…,xiD),xij∈[li,uj],i＝1,2,…,N,j＝1,2,…,D；

[0096] (2)利用深度学习策略，产生所述初始控制集合体RP的对应反向控制集 合体OP；其中，第i个控制热源体的位置为：

[0097] Xi'＝(xi'1,xi'2,…,xi'D),i＝1,2,…,N；

[0098] 式中，N为正整数；

[0099] 第i个控制热源体的第j维分量为：

[0100] xi'j＝lj+uj-xij；

[0101] 式中，lj为l为拉普拉斯空间第j维分量；uj为l为热源空间第j维分量； xij为第i个控制热源体第j个控制部分；

[0102] (3)合并所述一始控制集合体RP和所述反向控制集合体OP，将合并后 的2N个控制热源体按照其对应的优化控制度函数 的值大小采用最小化 方法进行升序排序，选取前N个控制热源体作为类快速成型架构方法的初始 控制集合体；

[0103] (4)计算所述初始控制集合体中每个带电控制热源体所受作用力，判断 是否满足结束条件，即步骤结束或者优化控制度函数值满足用户设定的设定 值，若满足，则结束，否则进行下一步；其中：

[0104] 计算所述初始控制集合体中每个带电控制热源体所受作用力，包括：

[0105] 计算控制热源体xi的电荷量qi的公式为：

[0106]

[0107] 式中，f(xi)为控制热源体xi的优化控制度函数，f(xbest)为最优的控制热源 体的优化控制度函数，n为正整数；

[0108] 根据库仑定律，对控制集合体中任意两个控制热源体xi和xj，计算它们之 间的作用力Fij为：

[0109]

[0110] 其中i,j＝1,2,…,N，N为正整数；每个控制热源体xi所受到的作用力Fi为：

[0111]

[0112] (5)将所述初始控制集合体中控制热源体分为精英控制热源体和一般控 制热源体。精英控制热源体是指优化控制度函数值小的前N1个控制热源体； N-N1个控制热源体则为一般控制热源体。对精英控制热源体进行单纯形控 制。单纯形控制法是一种处理无约束优化问题的直接方法，具有计算量小， 优化速度快，局部控制能力强，且不需要优化函数的梯度信息的特点。本实 施例具体步骤包括：

[0113] 1)给定初始单纯形参数：顶点的控制热源体xi∈Rn，i＝1,2,…,n+1(n为正 整数，其值可等于N或不等于N)；R为实数；反射系数α>0；扩展系数γ>1； 压缩系数β∈(0,1)；允许误差ε>0；

[0114] 计算顶点的优化控制度函数 i＝1,2,…,n+1，置k＝0；(k为次数序号)[0115] 2)确定位置最高点的控制热源体k次方为 位置次高点的控制热源体k 次方为位置最低点的控制热源体k次方为 (h,g,l∈{1,2,…,n+1})，确定优 化控制度函数为：

[0116]

[0117]

[0118]

[0119] 式中， 是位置最高点的优化控制度函数； 是位置次最高点的优 化控制度函数； 是位置最低点的优化控制度函数；

[0120] 计算n个热源体的形心 并计算形心处的优化控制度函数 其中计 算形心的公式为：

[0121]

[0122] 3)进行反射，令优化控制率r型的控制热源体k次方为：

[0123]

[0124] 式中， 是优化控制率r型的控制热源体的k次方； 形心处的控制热源体的 k次方； 是位置最高点的控制热源体的k次方；计算优化控制率r型的优化 控制度函数[0125] 4)判断优化控制率r型的优化控制度函数值若小于位置最低点的优化控 制度函数值，即 则进行扩展，令优化控制率e型的控制热源体 为：

[0126]

[0127] 计算优化控制率e型的优化控制度函数 且进入步骤5)；

[0128] 判断优化控制率r型的优化控制度函数值若不小于位置最低点的优化控 制度函数值且不大于位置次高点的优化控制度函数值，即 则置位置最高点的控制热源体k次方等于优化控制率r型的控制热源体k次 方，对应的优化控制度函数值也相等，即 且进入步骤 7)；

[0129] 判断优化控制率r型的优化控制度函数值若大于位置最高点的优化控制 度函数值大小，即 则进行压缩，令：

[0130]

[0131] 其中h′∈{h,r}， 为位置最高点中不易控制点的优化控制度函数；且令 优化控制率c型的控制热源体k次方为：

[0132]

[0133] 计算优化控制率c型的优化控制度函数 且进入步骤6)；需要说明 的是，在位置次高点和位置最高点中，有些不易控制的热源体点(由用户确 定)，其控制时可能采用h+r型控制率进行控制、h+e型控制率进行控制或者 h+c型控制率进行控制。

[0134] 5)判断优化控制率e型的优化控制度函数值若小于优化控制率r型的优化 控制度函数值，即 则令位置最高点的控制热源体k次方等于优 化控制率e型的控制热源体k次方，对应的优化控制度函数值也相等，即  且进入步骤6)；否则，令位置最高点的控制热源体k 次方等于优化控制率r型的控制热源体k次方，对应的优化控制度函数值也相 等，即 且进入步骤7)；

[0135] 6)判断优化控制率c型的优化控制度函数值若不大于位置最高点中不易 控制点的优化控制度函数值，即 则令位置最高点的控制热源体 k次方等于优化控制率c型的控 制热 源体k次方 ，对应的优化 控制度函数值也  相等 ，即且进入步骤7)；否则进行收缩，令：

[0136]

[0137] 其中i＝1,2,…,n+1，计算顶点的优化控制度函数 再进入步骤7)；

[0138] 7)判断顶点的优化控制度函数值 和形心处控制热源体的优化控制度 函数值 是否满足稳定控制准则，若：

[0139]

[0140] 则停止计算，且xi视为最优解；否则，令k的次数加1，即k＝k+1，并返 回步骤2)。

[0141] (6)移动控制热源体，更新群体中各控制热源体的位置；其公式为：

[0142]

[0143] 移动的控制热源体xi″将沿着作用力的方向移动。

[0144] (7)当控制热源体群稳定控制到一定程度，即控制热源体群中控制热源 体的组合，使其对应的优化控制度函数值在允许误差范围内时，对所述精英 控制热源体和所述一般控制热源体进行柯西变异的局部控制，得到全局最优 解。

[0145] 本实施例选取4个标准测试函数Sphere函数、Rosenbrock函数、Griewank 函数和Rastrigin函数进行测试。

[0146] 在4个测试函数中，Sphere函数是一个简单的单峰值函数，其他3个函数均 为复杂的多峰值函数，4个函数的全局最优值均为0。

[0147] 将本实施例提出的方法与修改的类快速成型架构方法、单纯形(SM)方 法得到的结果进行比较。在比较中，4个函数最大迭代次数设置为1000，稳定 控制精度设置为1e-05。在修改的类快速成型架构方法中，v＝0.5；在单纯形方 法中，α＝1.2,γ＝2,β＝0.3。在本实施例提供的方法中，v＝0.5，α＝1.2,γ＝2,β＝0.3。 变异概率pm＝0.05。每个测试问题在相同条件下独立运行20次实验，记录其最 优值和平均最优值，如图2中的数据所示，给出了在上述参数设置下，3种方 法对4个函数的优化操控结果比较。从数据中可知，对4个高维测试函数，本 实施例的方法对Sphere函数和Grienwank函数20次实验中一致地找到全局最优 解。

[0148] 图3至图6给出了本实施例的方法与单纯形方法和比较类快速成型架构方 法对4个测试问题的优化操控稳定控制曲线。从图3至图6可以清晰地看出，本 实施例提出的方法能快速地稳定控制到问题的全局最优解。

[0149] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普 通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润 饰，这些改进和润饰也应视为本发明的保护范围。