CN107145066A - 一种多参数优化方法 - Google Patents

一种多参数优化方法 Download PDF

Info

Publication number
CN107145066A
CN107145066A CN201710223294.8A CN201710223294A CN107145066A CN 107145066 A CN107145066 A CN 107145066A CN 201710223294 A CN201710223294 A CN 201710223294A CN 107145066 A CN107145066 A CN 107145066A
Authority
CN
China
Prior art keywords
particle
value
population
fitness value
vector
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
Application number
CN201710223294.8A
Other languages
English (en)
Other versions
CN107145066B (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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201710223294.8A priority Critical patent/CN107145066B/zh
Publication of CN107145066A publication Critical patent/CN107145066A/zh
Application granted granted Critical
Publication of CN107145066B publication Critical patent/CN107145066B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0265—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion

Landscapes

  • Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种多参数优化方法，能够加快学习效率，降低时间复杂度、空间复杂度。所述方法包括：初始化种群中每个粒子的位置和迭代次数，粒子每个维度对应一个待求解的参数；根据计算得到的每个粒子对目标函数的适应度值，更新每个粒子的历史最优位置和种群的最优位置，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量，并根据每个粒子的位置及学习向量确定下一代粒子的位置；若进行反向学习，则将每个粒子反向得到反向粒子，若反向粒子对目标函数的适应度值大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子。本发明适用于工业参数优化技术领域。

Description

一种多参数优化方法
技术领域
本发明涉及工业参数优化技术领域，特别是指一种多参数优化方法。
背景技术
近年来，越来越多的工业、企业开始重视大数据在优化产品性能上的重要之处，例如，使用大数据技术对工业生产中的工业参数进行设置，并进行相应的优化改进。然而，由于工业参数都存在着数据量大，目标参数维度高等问题，决定了工业参数优化是一个运算时间和空间都极高的机器学习问题。
发明内容
本发明要解决的技术问题是提供一种多参数优化方法，以解决现有技术所存在的时间复杂度和空间复杂度高的问题。
为解决上述技术问题，本发明实施例提供一种多参数优化方法，包括：
S1，初始化种群中每个粒子的位置和迭代次数，所述粒子处于多维空间中，粒子每个维度对应一个待求解的参数；
S2，计算每个粒子对目标函数的适应度值，根据计算得到的每个粒子对目标函数的适应度值，更新每个粒子的历史最优位置和种群的最优位置，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量，并根据每个粒子的位置及学习向量确定下一代粒子的位置，迭代次数加1；
S3，判断是否进行反向学习，若进行反向学习，则将每个粒子反向得到反向粒子，判断反向粒子对目标函数的适应度值是否大于原粒子对目标函数的适应度值，若大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子；
S4，检查是否满足预设的迭代结束条件，若满足预设的迭代结束条件，则输出结果；否则，则返回S2继续执行。
进一步地，所述S1包括：
S11，对待求解的多参数进行预处理，根据预处理结果，创建数学模型，并根据创建的数学模型，确定目标函数；其中，所述预处理包括：归一化处理，所述目标函数用于计算适应度值；
S12，在搜索空间中随机产生均匀分布的种群，记录种群中每个粒子的坐标向量，得到种群中每个粒子的位置；
S13，设置种群粒子的惯性系数、学习向量的学习系数，迭代次数初始值；
S14，设置迭代结束条件。
进一步地，所述S2包括：
S21，根据目标函数计算每个粒子的适应度值，将计算得到的适应度值与所述粒子的历史最优适应度值进行比较，若计算得到的适应度值优于历史最优适应度值，则用所述粒子的位置更新所述粒子的历史最优位置，并用所述粒子的适应度值更新所述粒子的历史最优适应度值；
S22，根据更新后的每个粒子的历史最优适应度值，获取适应度值最优的粒子的位置作为种群的最优位置；
S23，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量；
S24，根据种群粒子的惯性系数、学习向量的学习系数、每个粒子的当前位置及学习向量确定下一代粒子的位置；
S25，将迭代次数加1。
进一步地，所述S23包括：
S231，根据第一公式，计算当前迭代得到的每个粒子的每一种取值的平均适应度值，所示第一公式表示为：
其中，n表示第n个参数，q表示第q个取值，M表示第q个取值出现了M次，fm表示当前迭代得到的粒子的对应于所述取值的适应度值，Zmnq表示第q个取值的数值；
S232，根据得到的每个粒子的每一种取值的平均适应度值，从每个粒子的历史最优位置与种群的最优位置的每个维度取值中选取值较优的作为相应粒子的学习向量在相应维度的值。
进一步地，所述S24包括：
通过第二公式，确定下一代粒子的位置，所述第二公式表示为：
Xi+1＝α*Xi+β*Xlearni
其中，α表示种群粒子的惯性系数，β表示学习向量的学习系数，Xi表示第i次迭代中粒子的位置，Xlearni表示第i个粒子的学习向量，Xi+1表示第i+1次迭代中粒子的位置。
进一步地，所述S3包括：
S31，生成一个随机数，若随机数小于预设的第一阈值，则跳过S32，否则，则执行S32；
S32，将每个粒子反向得到反向粒子，计算反向粒子对目标函数的适应度值，判断反向粒子对目标函数的适应度值是否大于原粒子对目标函数的适应度值，若大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子；若不大于原粒子对目标函数的适应度值，则丢弃反向粒子。
进一步地，所述位置为向量；
所述将每个粒子反向得到反向粒子包括：
对粒子位置向量的每一个维度，取搜索空间在相应维度的中点为对称中心，则反向值为粒子原位置向量在相应维度关于对称中心的对称值；
将粒子位置向量的每个维度都取相应的对称值，得到所述粒子位置向量的反向向量；
根据所述粒子位置向量的反向向量，得到反向粒子。
进一步地，所述方法还包括：
若不进行反向学习，则执行S4。
进一步地，所述S4包括：
判断当前次迭代和上一次迭代得到的种群最优位置之间的差值是否小于预设的第二阈值；
若小于预设的第二阈值，则输出当前种群最优位置作为待求解的多参数的最优解，算法结束；
若不小于预设的第二阈值，则返回S2继续执行。
进一步地，所述S4还包括：
判断迭代次数是否达到预设的最大迭代次数；
若达到预设的最大迭代次数，则算法结束，求解失败；
若没有达到预设的最大迭代次数，则返回S2继续执行。
本发明的上述技术方案的有益效果如下：
上述方案中，通过正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量，并根据每个粒子的位置及学习向量确定下一代粒子的位置；并通过反向学习将每个粒子反向得到反向粒子，若反向粒子对目标函数的适应度值大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子。这样，通过正交算法可以加快进化速度及提高收敛速度，减少迭代的次数，而反向学习可以使迭代不容易陷入局部最优解，同时使用正交算法和反向学习能够加快学习效率，降低时间复杂度、空间复杂度。
附图说明
图1为本发明实施例提供的多参数优化方法的流程示意图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
本发明针对现有的时间复杂度和空间复杂度高的问题，提供一种多参数优化方法。
如图1所示，本发明实施例提供的多参数优化方法，包括：
S1，初始化种群中每个粒子的位置和迭代次数，所述粒子处于多维空间中，粒子每个维度对应一个待求解的参数；
S2，计算每个粒子对目标函数的适应度值，根据计算得到的每个粒子对目标函数的适应度值，更新每个粒子的历史最优位置和种群的最优位置，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量，并根据每个粒子的位置及学习向量确定下一代粒子的位置，迭代次数加1；
S3，判断是否进行反向学习，若进行反向学习，则将每个粒子反向得到反向粒子，判断反向粒子对目标函数的适应度值是否大于原粒子对目标函数的适应度值，若大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子；
S4，检查是否满足预设的迭代结束条件，若满足预设的迭代结束条件，则输出结果；否则，则返回S2继续执行。
本发明实施例所述的多参数优化方法，通过正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量，并根据每个粒子的位置及学习向量确定下一代粒子的位置；并通过反向学习将每个粒子反向得到反向粒子，若反向粒子对目标函数的适应度值大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子。这样，通过正交算法可以加快进化速度及提高收敛速度，减少迭代的次数，而反向学习可以使迭代不容易陷入局部最优解，同时使用正交算法和反向学习能够加快学习效率，降低时间复杂度、空间复杂度。
本实施例中，种群中的每个粒子处于多维空间中，每个粒子有多个维度/多维坐标，每个维度/每维坐标代表一个参数，因此，每个粒子的位置可以用于表示待求解的多个参数。粒子在搜索空间中运动，在每个位置都计算一个目标函数值，当目标函数值达到最小时，就是粒子的位置为最优解，粒子的多维坐标就是待求解的多个参数的最优解。
本实施例中，搜索空间是指求解过程中粒子的运动空间，也就是带求解参数的可能空间。搜索空间是一个多维空间，也是最开始粒子初始化时随机分布的空间。
在前述多参数优化方法的具体实施方式中，进一步地，所述初始化种群中每个粒子的位置和迭代次数(S1)包括：
S11，对待求解的多参数进行预处理，根据预处理结果，创建数学模型，并根据创建的数学模型，确定目标函数；其中，所述预处理包括：归一化处理，所述目标函数用于计算适应度值；
S12，在搜索空间中随机产生均匀分布的种群，记录种群中每个粒子的坐标向量，得到种群中每个粒子的位置；
S13，设置种群粒子的惯性系数、学习向量的学习系数，迭代次数初始值；
S14，设置迭代结束条件。
本实施例中，第一次迭代时，需对待求解的多个参数进行归一化等预处理操作，并根据预处理后的结果创建数学模型，根据创建的数学模型，找到目标函数f，所述目标函数f用于计算适应度值；进而，在搜索空间中随机产生均匀分布的种群，所述种群也可称为粒子群，记录粒子群中每个粒子的坐标向量，得到粒子群中每个粒子的位置Xi，其中，粒子的位置是向量，向量的维度等于待求解的参数的数目，粒子群规模根据所求参数的复杂程度设置，例如，粒子群规模可以在20～1000之间；在第一次迭代时，还需设置种群粒子的惯性系数α、学习向量的学习系数β，迭代次数iter初始值(iter＝0)及迭代结束条件。
在前述多参数优化方法的具体实施方式中，进一步地，所述计算每个粒子对目标函数的适应度值，根据计算得到的每个粒子对目标函数的适应度值，更新每个粒子的历史最优位置和种群的最优位置，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量，并根据每个粒子的位置及学习向量确定下一代粒子的位置，迭代次数加1(S2)包括：
S21，根据目标函数计算每个粒子的适应度值，将计算得到的适应度值与所述粒子的历史最优适应度值进行比较，若计算得到的适应度值优于历史最优适应度值，则用所述粒子的位置更新所述粒子的历史最优位置，并用所述粒子的适应度值更新所述粒子的历史最优适应度值；
S22，根据更新后的每个粒子的历史最优适应度值，获取适应度值最优的粒子的位置作为种群的最优位置；
S23，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量；
S24，根据种群粒子的惯性系数、学习向量的学习系数、每个粒子的当前位置及学习向量确定下一代粒子的位置；
S25，将迭代次数加1。
本实施例中，在第g次迭代时，首选根据目标函数f计算每个粒子的适应度值Pfit，将计算得到的当前适应度值Pfit与所述粒子的历史最优适应度值进行比较，若所述粒子的当前适应度值Pfit优于历史最优适应度值，则用所述粒子的当前位置更新所述粒子的历史最优位置Pbest，并用所述粒子的当前适应度值Pfit更新所述粒子的历史最优适应度值；按照更新后的每个粒子的历史最优适应度值对所有粒子的历史最优位置进行排序，获取适应度值最优的粒子的位置作为种群的最优位置Gbest，其中，适应度值越小越优；利用正交算法对每个粒子的历史最优位置Pbest与种群的最优位置Gbest进行正交计算，得到每个粒子的学习向量Xlearni；根据种群粒子的惯性系数、学习向量的学习系数、每个粒子的当前位置及学习向量确定下一代粒子的位置，每次更新子代都会替换原来的粒子；并将迭代次数iter加1。
本实施例中，历史最优位置Pbest与种群的最优位置Gbest都是向量，每个向量都有多个维度，每个维度表示一个参数，且每个维度上有多个取值。
在前述多参数优化方法的具体实施方式中，进一步地，所述利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量(S23)包括：
S231，因为每个维度上有多个取值，可以根据第一公式，计算当前迭代得到的每个粒子的每一种取值的平均适应度值，所示第一公式表示为：
其中，n表示第n个参数，q表示第q个取值，M表示第q个取值出现了M次，fm表示当前迭代得到的粒子的对应于所述取值的适应度值，Zmnq表示第q个取值的数值；
S232，根据得到的每个粒子的每一种取值的平均适应度值，从每个粒子的历史最优位置与种群的最优位置的每个维度取值中选取值较优的作为相应粒子的学习向量在相应维度的值。
本实施例中，带求解的参数以钢铁轧制的控制参数为例，对正交算法进行举例说明：
假设，钢铁轧制的控制参数包括：温度和压力，其中，温度和压力的数值组成一个二维向量，就是粒子的位置坐标；例如：温度可以有100℃，200℃，300℃；压力有100N，200N。
这样两两有6种组合，并根据目标函数计算每个粒子的适应度值，得到6个适应度值，如表1所示：
表1适应度值
温度 压力 适应度值(越小越优)
100℃ 100N 10
100℃ 200N 14
200℃ 100N 16
200℃ 200N 26
300℃ 100N 22
300℃ 200N 24
计算正交时，对温度这个参数，需要分别计算100℃，200℃，300℃的适应度平均值。
100℃平均适应度值为(10+14)/2＝12；
200℃平均适应度值为(16+26)/2＝21；
300℃平均适应度值为(22+24)/2＝23；
100N平均适应度值为(10+16+22)/3＝16；
200N平均适应度值为(14+26+24)/3＝18；
接着，对得到的每个粒子的每一种取值的平均适应度值进行排序，平均适应度值越小越优，则可知对温度参数而言，温度值越小越优；对压力而言，压力值越小越优；根据得到的每个粒子的每一种取值的平均适应度值，从每个粒子的历史最优位置与种群的最优位置的每个维度取值中选取值较优的作为相应粒子的学习向量在相应维度的值，从而得到当前迭代每个粒子的学习向量，进而可以根据得到的每个粒子的学习向量指导相应粒子的下一次飞行；
本实施例中，假设粒子1的历史最优位置是(300℃，100N)；
粒子2的历史最优位置是(100℃，300N)；
种群的最优位置是(200℃，200N)；
则选择学习向量时，温度上粒子1选择全局的，压力选自身历史的，组成的粒子1的学习向量就是(200℃，100N)；同理，粒子2的学习向量就是(100℃，200N)。
在前述多参数优化方法的具体实施方式中，进一步地，所述根据种群粒子的惯性系数、学习向量的学习系数、每个粒子的当前位置及学习向量确定下一代粒子的位置(S24)包括：
通过第二公式，确定下一代粒子的位置，所述第二公式表示为：
Xi+1＝α*Xi+β*Xlearni
其中，α表示种群粒子的惯性系数，β表示学习向量的学习系数，Xi表示第i次迭代中粒子的位置，Xlearni表示第i个粒子的学习向量，Xi+1表示第i+1次迭代中粒子的位置。
在前述多参数优化方法的具体实施方式中，进一步地，所述判断是否进行反向学习，若进行反向学习，则将每个粒子反向得到反向粒子，判断反向粒子对目标函数的适应度值是否大于原粒子对目标函数的适应度值，若大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子(S3)包括：
S31，生成一个随机数，若随机数小于预设的第一阈值，则跳过S32，否则，则执行S32；
S32，将每个粒子反向得到反向粒子，计算反向粒子对目标函数的适应度值，判断反向粒子对目标函数的适应度值是否大于原粒子对目标函数的适应度值，若大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子；若不大于原粒子对目标函数的适应度值，则丢弃反向粒子。
本实施例中，例如，可以生成一个随机取值为(0,1)的随机数randiter，并假设所述预设的第一阈值为0.3，即：若随机数randiter<0.3，则跳过步骤32，否则执行步骤32。
本实施例中，在执行步骤32时，每个粒子反向得到反向粒子Preverse，计算反向粒子Preverse对目标函数f的适应度值P′fit，判断反向粒子Preverse对目标函数f的适应度值P′fit是否大于原粒子(反向之前的该粒子)对目标函数f的适应度值Pfit，若大于原粒子对目标函数f的适应度值Pfit，则利用反向粒子Preverse替换原粒子；若不大于原粒子对目标函数f的适应度值Pfit，则丢弃反向粒子Preverse。
本实施例中，粒子的属性有空间位置和速度矢量，所以替换粒子就是更新这两个属性的值。由于反向学习不是一次独立的迭代(因为：不是每个粒子都会进行反向学习，只有当生成的随机数小于预设的第一阈值时，相应的粒子才会进行反向学习)，所以粒子的下标不变，就相当于是把原粒子替换掉了，原粒子被替换后：位置就是反向后的反向粒子Preverse的位置，速度变成0。
在前述多参数优化方法的具体实施方式中，进一步地，所述位置为向量；
所述将每个粒子反向得到反向粒子包括：
对粒子位置向量的每一个维度，取搜索空间在相应维度的中点为对称中心，则反向值为粒子原位置向量在相应维度关于对称中心的对称值；
将粒子位置向量的每个维度都取相应的对称值，得到所述粒子位置向量的反向向量；
根据所述粒子位置向量的反向向量，得到反向粒子。
本实施例中，用Xi表示第i个粒子的位置，粒子位置Xi是向量，对第i个粒子位置向量Xi的每一个维度，取搜索空间在相应维度的中点为对称中心，则反向值为粒子原位置向量Xi在相应维度关于对称中心的对称值，将粒子的位置向量Xi每个维度都取相应的对称值，则得到所述粒子的位置向量Xi的反向向量Xireverse；根据反向向量Xireverse，得到反向粒子Preverse。
在前述多参数优化方法的具体实施方式中，进一步地，所述方法还包括：
若不进行反向学习，则执行S4。
在前述多参数优化方法的具体实施方式中，进一步地，所述检查是否满足预设的迭代结束条件，若满足预设的迭代结束条件，则输出结果；否则，则返回S2继续执行(S4)包括：
判断当前次迭代和上一次迭代得到的种群最优位置之间的差值是否小于预设的第二阈值；
若小于预设的第二阈值，则输出当前种群最优位置作为待求解的多参数的最优解，算法结束；
若不小于预设的第二阈值，则返回S2继续执行。
本实施例中，检查最后两次(当前次迭代和上一次迭代)迭代得到的种群最优位置Gbest之间的差值，若差值小于预设的第二阈值，例如，所述预设的第二阈值为10-6，则停止迭代，输出当前种群最优位置Gbest作为待求解的多参数的最优解，算法结束；若差值不小于预设的第二阈值，则返回S2继续执行。
在前述多参数优化方法的具体实施方式中，进一步地，所述检查是否满足预设的迭代结束条件，若满足预设的迭代结束条件，则输出结果；否则，则返回S2继续执行(S4)还包括：
判断迭代次数是否达到预设的最大迭代次数；
若达到预设的最大迭代次数，则算法结束，求解失败；
若没有达到预设的最大迭代次数，则返回S2继续执行。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种多参数优化方法，其特征在于，包括：
S1，初始化种群中每个粒子的位置和迭代次数，所述粒子处于多维空间中，粒子每个维度对应一个待求解的参数；
S2，计算每个粒子对目标函数的适应度值，根据计算得到的每个粒子对目标函数的适应度值，更新每个粒子的历史最优位置和种群的最优位置，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量，并根据每个粒子的位置及学习向量确定下一代粒子的位置，迭代次数加1；
S3，判断是否进行反向学习，若进行反向学习，则将每个粒子反向得到反向粒子，判断反向粒子对目标函数的适应度值是否大于原粒子对目标函数的适应度值，若大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子；
S4，检查是否满足预设的迭代结束条件，若满足预设的迭代结束条件，则输出结果；否则，则返回S2继续执行。
2.根据权利要求1所述的多参数优化方法，其特征在于，所述S1包括：
S11，对待求解的多参数进行预处理，根据预处理结果，创建数学模型，并根据创建的数学模型，确定目标函数；其中，所述预处理包括：归一化处理，所述目标函数用于计算适应度值；
S12，在搜索空间中随机产生均匀分布的种群，记录种群中每个粒子的坐标向量，得到种群中每个粒子的位置；
S13，设置种群粒子的惯性系数、学习向量的学习系数，迭代次数初始值；
S14，设置迭代结束条件。
3.根据权利要求1所述的多参数优化方法，其特征在于，所述S2包括：
S21，根据目标函数计算每个粒子的适应度值，将计算得到的适应度值与所述粒子的历史最优适应度值进行比较，若计算得到的适应度值优于历史最优适应度值，则用所述粒子的位置更新所述粒子的历史最优位置，并用所述粒子的适应度值更新所述粒子的历史最优适应度值；
S22，根据更新后的每个粒子的历史最优适应度值，获取适应度值最优的粒子的位置作为种群的最优位置；
S23，利用正交算法对每个粒子的历史最优位置与种群的最优位置进行正交计算，得到每个粒子的学习向量；
S24，根据种群粒子的惯性系数、学习向量的学习系数、每个粒子的当前位置及学习向量确定下一代粒子的位置；
S25，将迭代次数加1。
4.根据权利要求3所述的多参数优化方法，其特征在于，所述S23包括：
S231，根据第一公式，计算当前迭代得到的每个粒子的每一种取值的平均适应度值，所示第一公式表示为：
<mrow> <msub> <mi>S</mi> <mrow> <mi>n</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </msubsup> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>*</mo> <msub> <mi>Z</mi> <mrow> <mi>m</mi> <mi>n</mi> <mi>q</mi> </mrow> </msub> </mrow> <mrow> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </msubsup> <msub> <mi>Z</mi> <mrow> <mi>m</mi> <mi>n</mi> <mi>q</mi> </mrow> </msub> </mrow> </mfrac> </mrow>
其中，n表示第n个参数，q表示第q个取值，M表示第q个取值出现了M次，fm表示当前迭代得到的粒子的对应于所述取值的适应度值，Zmnq表示第q个取值的数值；
S232，根据得到的每个粒子的每一种取值的平均适应度值，从每个粒子的历史最优位置与种群的最优位置的每个维度取值中选取值较优的作为相应粒子的学习向量在相应维度的值。
5.根据权利要求3所述的多参数优化方法，其特征在于，所述S24包括：
通过第二公式，确定下一代粒子的位置，所述第二公式表示为：
Xi+1＝α*Xi+β*Xlearni
其中，α表示种群粒子的惯性系数，β表示学习向量的学习系数，Xi表示第i次迭代中粒子的位置，Xlearni表示第i个粒子的学习向量，Xi+1表示第i+1次迭代中粒子的位置。
6.根据权利要求1所述的多参数优化方法，其特征在于，所述S3包括：
S31，生成一个随机数，若随机数小于预设的第一阈值，则跳过S32，否则，则执行S32；
S32，将每个粒子反向得到反向粒子，计算反向粒子对目标函数的适应度值，判断反向粒子对目标函数的适应度值是否大于原粒子对目标函数的适应度值，若大于原粒子对目标函数的适应度值，则利用反向粒子替换原粒子；若不大于原粒子对目标函数的适应度值，则丢弃反向粒子。
7.根据权利要求1或6所述的多参数优化方法，其特征在于，所述位置为向量；
所述将每个粒子反向得到反向粒子包括：
对粒子位置向量的每一个维度，取搜索空间在相应维度的中点为对称中心，则反向值为粒子原位置向量在相应维度关于对称中心的对称值；
将粒子位置向量的每个维度都取相应的对称值，得到所述粒子位置向量的反向向量；
根据所述粒子位置向量的反向向量，得到反向粒子。
8.根据权利要求1所述的多参数优化方法，其特征在于，所述方法还包括：
若不进行反向学习，则执行S4。
9.根据权利要求1所述的多参数优化方法，其特征在于，所述S4包括：
判断当前次迭代和上一次迭代得到的种群最优位置之间的差值是否小于预设的第二阈值；
若小于预设的第二阈值，则输出当前种群最优位置作为待求解的多参数的最优解，算法结束；
若不小于预设的第二阈值，则返回S2继续执行。
10.根据权利要求1或9所述的多参数优化方法，其特征在于，所述S4还包括：
判断迭代次数是否达到预设的最大迭代次数；
若达到预设的最大迭代次数，则算法结束，求解失败；
若没有达到预设的最大迭代次数，则返回S2继续执行。
CN201710223294.8A 2017-04-07 2017-04-07 一种多参数优化方法 Active CN107145066B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710223294.8A CN107145066B (zh) 2017-04-07 2017-04-07 一种多参数优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710223294.8A CN107145066B (zh) 2017-04-07 2017-04-07 一种多参数优化方法

Publications (2)

Publication Number Publication Date
CN107145066A true CN107145066A (zh) 2017-09-08
CN107145066B CN107145066B (zh) 2020-09-11

Family

ID=59774694

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710223294.8A Active CN107145066B (zh) 2017-04-07 2017-04-07 一种多参数优化方法

Country Status (1)

Country Link
CN (1) CN107145066B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111950802A (zh) * 2020-08-21 2020-11-17 上海明略人工智能(集团)有限公司 一种生产调度控制方法和装置
CN113378468A (zh) * 2021-06-18 2021-09-10 中国科学院地理科学与资源研究所 多维地学参量的权重优化方法及系统
CN114599004A (zh) * 2022-01-28 2022-06-07 北京邮电大学 一种基站布局方法和装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB202110582D0 (en) * 2021-07-22 2021-09-08 Owlstone Med Ltd Breath sampling device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05173625A (ja) * 1991-12-20 1993-07-13 Komatsu Ltd Ｆａシステムの下位モジュール識別方法
CN103544526A (zh) * 2013-11-05 2014-01-29 辽宁大学 一种改进粒子群算法及其应用
CN106371316A (zh) * 2016-10-25 2017-02-01 大唐（北京）水务工程技术有限公司 基于pso‑lssvm的水岛加药在线控制方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05173625A (ja) * 1991-12-20 1993-07-13 Komatsu Ltd Ｆａシステムの下位モジュール識別方法
CN103544526A (zh) * 2013-11-05 2014-01-29 辽宁大学 一种改进粒子群算法及其应用
CN106371316A (zh) * 2016-10-25 2017-02-01 大唐（北京）水务工程技术有限公司 基于pso‑lssvm的水岛加药在线控制方法和装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
WEI HONG LIM 等: "A Self-Adaptive Topologically Connected-Based Particle Swarm Optimization", 《IEEE》 *
ZHI-HUI ZHAN 等: "Orthogonal Learning Particle Swarm Optimization", 《IEEE TRANSACTIONS ON EVOLUTIONARY COMPUTATION》 *
王利霞: "改进的粒子群优化算法及其在数据分类中的应用", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
肖文显 等: "一种融合反向学习和量子优化的粒子群算法", 《微电子学与计算机》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111950802A (zh) * 2020-08-21 2020-11-17 上海明略人工智能(集团)有限公司 一种生产调度控制方法和装置
CN113378468A (zh) * 2021-06-18 2021-09-10 中国科学院地理科学与资源研究所 多维地学参量的权重优化方法及系统
CN113378468B (zh) * 2021-06-18 2024-03-29 中国科学院地理科学与资源研究所 多维地学参量的权重优化方法及系统
CN114599004A (zh) * 2022-01-28 2022-06-07 北京邮电大学 一种基站布局方法和装置
CN114599004B (zh) * 2022-01-28 2024-01-05 北京邮电大学 一种基站布局方法和装置

Also Published As

Publication number Publication date
CN107145066B (zh) 2020-09-11

Similar Documents

Publication Publication Date Title
CN107145066A (zh) 一种多参数优化方法
CN104834215B (zh) 一种变异粒子群优化的bp神经网络pid控制算法
Li et al. A hybrid assembly sequence planning approach based on discrete particle swarm optimization and evolutionary direction operation
CN105138717A (zh) 一种动态变异粒子群优化神经网络的变压器状态评估方法
CN104317195B (zh) 一种基于改进极限学习机的非线性逆模型控制方法
CN111079891A (zh) 一种基于双隐含层bp神经网络的离心泵性能预测方法
CN104615856B (zh) 基于热风炉组的煤气消耗量预测模型建立方法及装置
WO2023207139A1 (zh) 利用交流型樽海鞘算法求解拉伸/压缩弹簧参数的方法
CN107578093A (zh) 滑坡变形的Elman神经网络动态预测方法
CN111522233A (zh) 参数自整定的mimo异因子全格式无模型控制方法
CN105469142A (zh) 一种基于样本增量驱动的神经网络增量型前馈算法
CN115048804A (zh) 一种城市管网多传感器优化部署方法
CN109193807A (zh) 一种电力系统经济调度方法及系统
JP2018180799A5 (zh)
Masrom et al. Hybridization of particle swarm optimization with adaptive genetic algorithm operators
CN109074348A (zh) 用于对输入数据集进行迭代聚类的设备和迭代方法
CN109782586B (zh) 参数自整定的miso异因子紧格式无模型控制方法
CN114186709A (zh) 基于帝王蝶算法优化灰色模型关键参数的能源预测方法
CN106292289B (zh) 流程工业控制回路对象的混合精英随机搜索优化方法
Hamza et al. Differential evolution with multi-constraint consensus methods for constrained optimization
CN111279276B (zh) 用于控制复杂系统的随机化加强学习
CN103105775A (zh) 基于序优化与在线核极限学习机的分层迭代优化调度方法
Jun et al. Models and algorithms of production scheduling in tandem cold rolling
CN106371321A (zh) 一种焦化炉炉膛压力系统模糊网络优化pid控制方法
CN116979516A (zh) 基于rngo算法的多机电力系统稳定器参数优化方法

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