CN115130749A - 融合nsga-ⅲ和topsis的数据驱动多目标优化方法 - Google Patents
融合nsga-ⅲ和topsis的数据驱动多目标优化方法 Download PDFInfo
- Publication number
- CN115130749A CN115130749A CN202210751216.6A CN202210751216A CN115130749A CN 115130749 A CN115130749 A CN 115130749A CN 202210751216 A CN202210751216 A CN 202210751216A CN 115130749 A CN115130749 A CN 115130749A
- Authority
- CN
- China
- Prior art keywords
- evaluation
- optimization
- data
- distance
- 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.)
- Pending
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims abstract description 106
- 230000006870 function Effects 0.000 claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 25
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000013135 deep learning Methods 0.000 claims abstract description 10
- 230000002068 genetic effect Effects 0.000 claims abstract description 9
- 238000007781 pre-processing Methods 0.000 claims abstract description 8
- 230000002159 abnormal effect Effects 0.000 claims abstract description 4
- 238000011156 evaluation Methods 0.000 claims description 81
- 239000011159 matrix material Substances 0.000 claims description 49
- 239000013598 vector Substances 0.000 claims description 14
- 238000013527 convolutional neural network Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 8
- 239000003795 chemical substances by application Substances 0.000 claims description 7
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 abstract description 13
- 210000004027 cell Anatomy 0.000 description 14
- 238000005520 cutting process Methods 0.000 description 13
- 238000002474 experimental method Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- 238000011176 pooling Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000005265 energy consumption Methods 0.000 description 6
- 230000003746 surface roughness Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000013400 design of experiment Methods 0.000 description 3
- 238000003801 milling Methods 0.000 description 3
- 230000004913 activation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003050 experimental design method Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000012803 optimization experiment Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Genetics & Genomics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种融合NSGA‑Ⅲ和TOPSIS的数据驱动多目标优化方法,包括如下步骤:步骤1):历史数据采集:采集与优化目标和工艺参数相关的数据;步骤2):数据预处理:对采集得到的历史数据进行预处理,剔除原始数据中的异常值;步骤3):生成代理模型:利用深度学习方法,根据历史数据生成优化目标的预测模型,并将预测模型转换成代理模型;步骤4):多目标优化:将代理模型作为目标函数,利用非支配排序遗传算法(NSGA‑Ⅲ)对目标函数进行优化,生成Pareto前沿;步骤5):选取最优解:利用逼近理想解排序法(TOPSIS)对Pareto解集进行排序,确定最优解,能准确的得到不同的生产任务和生产要求对应的最优方案,解决生产过程中的多目标优化问题。
Description
技术领域
本发明属于多目标优化技术领域,具体的为一种融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法。
背景技术
多目标优化是指某个情景中在需要达到多个目标时,由于目标间存在内在冲突,一个目标的优化是以其他目标劣化为代价,因此很难出现唯一最优解,而只能在它们中间进行协调和折中处理,使各个子目标都尽可能地达到最优。多目标优化问题在工程应用等现实生活中非常普遍并且处于非常重要的地位,相较于单目标优化问题,多目标优化问题的求解难度更大,解决多目标优化问题具有非常重要的科研价值和实际意义。
根据国内外学者关于多目标优化问题的研究,目前解决多目标优化问题可分为两步。第一步是通过实验技术和理论方法建立参数和不同优化目标之间的数学关系。第二步是开发一种多目标优化算法或组合算法,用于解决参数的优化问题。
第一步中,实验技术主要采用试验设计(DoE)的方法来构造统计曲线拟合函数,并将其作为预测模型来实现优化过程。DoE的有效使用已被证明在筛选最有候选参数和确定优化的因果效应方面至关重要,从而提高产品的质量和可靠性。在实际制造操作中,使用精确的实验方法构建曲线拟合函数通常需要昂贵的实验加工数据采集,导致人工成本和材料成本增加。采用理论方法进行优化主要是基于物理-数学模型或有限元模拟。理论方法在实际应用前就可以得到合理的参数,能提高加工工艺性能和降低生产成本,但是该方法依赖于专家知识去构建多目标优化函数。
第二步中,“Multi-objective parameter optimization of turbine impellerbased on RBF neural network and NSGA-II genetic algorithm”(Energy Reports,2021,7(S7),Ji Yunguang,Yang Zhikuo,Ran Jingyu,Li Hongtao)内记载了一种基于径向基函数神经网络和NSGA-II的多目标优化方法,用于优化涡轮的扬程和效率。但是该方法针对扬程和效率这两个相互冲突的目标,需要额外的时间和人力资源用于筛选出最优的参数组合。“Hybrid NSGA-II based decision-making in fuzzy multi-objectivereliability optimization problem”(SN Applied Sciences,2019,1(11),HemantKumar,Shiv Prasad Yadav)内记载了一种利用NSGA-II优化燃气轮机超速保护系统的系统可靠性、系统成本和系统重量的技术方案,利用模糊排序方法从Pareto解集中得到最优折衷解。但是该方法依赖于专家知识去确定每一个目标的隶属度函数,不能快速应用于具有新目标的优化任务。公开号为CN114358409A的中国专利中公开了一种多目标优化结果的排序方法及相关装置,通过多次多目标优化将多目标优化结果的排序问题转换为单目标排序问题,降低了多目标优化排序问题的复杂度。但是该方法需要人为的去干预,优化步骤复杂且优化时间长。
综上所述,现有的多目标优化方法存在以下缺陷:
(i)现有方法的优化时间长、成本高。试验设计方法针对新的优化目标需要进行新的实验,并依赖于昂贵的实验数据采集,然后手动获得每个目标的曲线拟合函数,导致人工成本和材料成本增加。采用理论方法时物理预测模型和仿真模型的计算效率往往较低,导致迭代过程的时间成本较高,优化耗时几个小时到几天不等。
(ii)现有方法的依赖于专家知识,方法的适应性差。现有的方法依赖于专家知识去构建目标函数,但是针对多目标优化任务,由于专家知识仅在特定领域发挥作用,对于新的优化目标仅依赖于专家知识去开发不同目标的目标函数是很困难的。在实际的操作中存在各种类型的任务,不同任务的目标和相关的参数也不相同,由于专家知识的有限性,传统的方法不适用于解决不同领域的多优化问题。
(iii)现有方法不能准确获得最优方案。多目标优化求解是从基于物理预测模型的目标函数获得优化后的参数,然而物理建模过程中经常进行许多假设和简化,导致目标函数的精度较低,优化后的参数可能不满足约束条件。此外,在生成的Pareto解集中只有少数几个解代表了优化目标之间的最优权衡,导致决策者很难根据目标函数从Pareto解集中做出最优选择。
发明内容
有鉴于此,本发明的目的在于提供一种融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,能准确的得到不同的生产任务和生产要求对应的最优方案,解决生产过程中的多目标优化问题,以快速、准确地获得最优方案。
为达到上述目的,本发明提供如下技术方案:
一种融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,包括如下步骤:
步骤1):历史数据采集:采集与优化目标和工艺参数相关的数据;
步骤2):数据预处理:对采集得到的历史数据进行预处理,剔除原始数据中的异常值;
步骤3):生成代理模型:利用深度学习方法,根据历史数据生成优化目标的预测模型,并将预测模型转换成代理模型;
步骤4):多目标优化:将代理模型作为目标函数,利用非支配排序遗传算法(NSGA-Ⅲ)对目标函数进行优化,生成Pareto前沿;
步骤5):选取最优解:利用逼近理想解排序法(TOPSIS)对Pareto解集进行排序,确定最优解。
进一步,所述步骤2)中,利用局部离群因子算法(LOF)对历史数据进行预处理,方法如下:
在样本集D中,定义某样本点o与离其第k远的近邻点间的距离为点o的第k距离dk(o),并定义点o的第k距离邻域Nk(o)为离点o的距离不超过dk(o)的所有样本点,即:
Nk(o)={o'|o'∈D,d(o,o')≤dk(o)}
其中,d(o,o')为样本点o与o'间的距离;若d(o,o')<dk(o),则定义样本点o的可达距离dreach(o,o')为dk(o),否则为d(o,o'),即:
dreach(o,o')=max{d(o,o'),dk(o)}
其中,||Nk(o)||为点o的第k距离邻域内包含的样本点个数;
定义点o的局部可达密度与其第k距离邻域Nk(o)内所有点的局部可达密度之比的平均值为局部离群因子LOFk(o),即:
样本点o的局部可达密度相比其k近邻的局部可达密度越低,则LOFk(o)值越大,说明点o是离群点的可能性也越大。
进一步,所述步骤3)中,若历史数据为离散数据,则采用卷积神经网络(CNN)生成优化目标的预测模型;若历史数据为时序数据,则采用长短期记忆人工神经网络(LSTM)生成优化目标的预测模型。
进一步,所述步骤4)中,利用非支配排序遗传算法(NSGA-Ⅲ)对目标函数进行优化的方法为:
41)通过边界交叉构造权重的方法,基于优化目标数M和各维度目标划分数p,均匀生成H(M,p)个参考点:
42)将第t代种群Pt通过交叉变异操作得到子代种群Qt,每代种群规模为S,将父代种群和子代种群合并形成新种群Rt;
43)将合并种群Rt中的个体进行快速非支配排序,得到若干非支配层F1,F2,F3,...,Fl;
44)将优先级较高的非支配层存入下一代种群Pt+1,临界层为Fτ,满足:
45)得到当前种群的所有个体在每一维目标上的最优值,构成当前种群在每一维目标上的原点。将步骤41)中生成的参考点与原点相连,连接线组成参考向量;计算当前Pt+1中所有个体到参考向量的垂直距离,并与最近的参考向量相关联;
判断迭代次数Gen是否达到设定阈值maxGen;若是,则输出当下种群,得到工艺参数优化解集;若否,则执行步骤42)。
进一步,所述步骤5)中,利用逼近理想解排序法(TOPSIS)确定最优解的方法为:
51)建立初始判断矩阵
设评价方案数量为m,所有评价方案构成方案集A={A1,A2,…,Am};每个评价方案有n个评价指标,单一评价方案评判集Ai={Xi1,Xi2,…,Xin},评判指标Xij为第i个评价方案的第j个评价指标值,i∈[1,m],j∈[1,n];构建相应的初始判断矩阵为:
52)建立标准化决策矩阵
将负向指标以倒数方式完成正向化,使所有指标具有同趋势性;将原始数据进行归一化处理以消除量纲:
最终得到标准化决策矩阵Z=(zij)m×n;
53)构建加权评价矩阵
对评价方案的n维评价指标分配权重,构建的评价指标权重矩阵对标准化决策矩阵进行加权得到加权评价矩阵Z’:
其中,Z'表示加权评价矩阵;W表示评价指标的权重矩阵;wj表示第j个评价指标的权重;
54)确定理想解
正理想解Z+是由加权评价矩阵Z’中每列最优值组成:
负理想解Z-是由加权评价矩阵Z’中每列最劣值组成:
55)计算距离尺度与优劣排序
(1)计算距离尺度,确定评价方案与正、负理想解的欧式距离:
其中,Di +表示第i个评价方案到正理想解的欧式距离;Di -表示第i个评价方案到正理想解的欧式距离;z’ij表示加权评价矩阵中的第i个方案得到的第j个指标的值;
(2)计算评价方案与正理想解的接近度:
其中,Ci表示第i个方案与正理想解的接近度;
根据Ci高低进行优劣排序,Ci值越接近1表明评价方案越优。
本发明的有益效果在于:
本发明融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法具有以下优点:
1)基于深度学习的代理模型构建方法可以自适应的构建优化参数和优化目标之间的复杂函数关系。与传统的方法相比,本方法不需要构建新的模型或实验也不依赖于专家知识,减少了优化时间和成本。
2)使用历史数据搭建目标函数,与传统的方法相比,深度学习建立代理模型,采用自动化、集成的方式,人为干扰最小,不需要理解生产过程中的物理规律或化学规律,避免了非最优隐藏节点问题,保证了模型精度的稳定。
3)可提高从Pareto非劣解集中寻找最优方案的搜索性能,与传统的方法相比,可以根据预先设置的目标数量和目标间的权重,自适应的获得符合生产过程要求的最优方案。
4)可以针对不同的生产任务和生产要求,对目标的个数和不同目标的权重进行动态调整,解决传统方法不能快速适应不同任务的优化问题。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法的流程图;
图2为CNN模型的结构图;
图3为LSTM模型的结构图;
图4为利用非支配排序遗传算法(NSGA-Ⅲ)对目标函数进行优化的流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
如图1所示,本实施例融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,包括如下步骤。
步骤1):历史数据采集:采集与优化目标和工艺参数相关的数据。
根据优化的目标和优化的参数,利用传感器和实验设备采集相关的历史数据。对于在生产过程中不能从设备中直接获取的状态信息,可以在需要采集信息的部位安装相应的传感器,通过数据采集系统获取历史数据。
步骤2):数据预处理:对采集得到的历史数据进行预处理,剔除原始数据中的异常值。
本实施例利用局部离群因子算法(LOF)对历史数据进行预处理,提高原始数据的准确性,方法如下:
在样本集D中,定义某样本点o与离其第k远的近邻点间的距离为点o的第k距离dk(o),并定义点o的第k距离邻域Nk(o)为离点o的距离不超过dk(o)的所有样本点,即:
Nk(o)={o'|o'∈D,d(o,o')≤dk(o)}
其中,d(o,o')为样本点o与o'间的距离;若d(o,o')<dk(o),则定义样本点o的可达距离dreach(o,o')为dk(o),否则为d(o,o'),即:
dreach(o,o')=max{d(o,o'),dk(o)}
其中,||Nk(o)||为点o的第k距离邻域内包含的样本点个数;
为了避免样本分布形状对离群点选取的影响,LOF算法采用局部可达密度来表征样本点o的密度。定义点o的局部可达密度lrdk(o)为:
其中,||Nk(o)||为点o的第k距离邻域内包含的样本点个数。
定义点o的局部可达密度与其第k距离邻域Nk(o)内所有点的局部可达密度之比的平均值为局部离群因子LOFk(o),即:
样本点o的局部可达密度相比其k近邻的局部可达密度越低,则LOFk(o)值越大,说明点o是离群点的可能性也越大。
步骤3):生成代理模型:利用深度学习方法,根据历史数据生成优化目标的预测模型,并将预测模型转换成代理模型。
本实施例中,代理模型是指利用各种预测模型所开发的目标函数。通常,代理模型是一个快速运行的逼近函数,可用于全局或局部逼近昂贵的函数值。由于复杂的工作条件阻碍了物理模型的构建来描述加工过程,使其优化极其困难。深度学习是机器学习中最重要的分支领域,由于具有从原始数据中提取特征的能力,能够自动建立输入参数和输出目标之间的关系,利用深度学习开发的预测模型可以在较少人为干预的情况下运行,所有模型参数,包括特征和回归模块,都可以联合训练。具体的,若历史数据为离散数据,则采用卷积神经网络(CNN)生成优化目标的预测模型;若历史数据为时序数据,则采用长短期记忆人工神经网络(LSTM)生成优化目标的预测模型。
CNN具有一般的一维结构,如图2所示。在一维CNN建模过程中,考虑了五个不同的层:输入层、卷积层、池化层、全连接层和输出层。采集到的数据,包括工艺参数和输出目标,送到输入层进行处理,再送到卷积层进行卷积运算。
在卷积层之后,池化层的作用是减少特征映射级别,同时保留之前输入中隐藏的信息。通常,多个池操作是必要的。在本实施例中,采用最大池化,可以表示为:
f(x)=max(0)
在实施例的池化操作为:
其中,x为池化层的输入,f(x)为池化层的输出。在卷积和池化操作后,后续为全连接层,以池化层输出的1D向量作为输入。全连接层由许多神经元组成,这些神经元与前一层中的所有节点聚集在一起。在输出层,激活函数通常用于完成回归或分类。由于预测问题是一项回归任务,因此采用线性激活函数。
LSTM作为特殊RNN,在循环网络的基础上添加记忆细胞(cell)单元,用于保存长期状态,并增加遗忘门、输入门和输出门共同控制信息的遗忘和记忆,使之更擅长处理时间序列数据。
如图3,LSTM的单元结构xt和ht分别为t时刻的输入和输出,ht-1为上一时刻LSTM单元的输出,ct-1为上一时刻记忆细胞单元的输出,ct为t时刻记忆细胞单元的输出。LSTM各单元计算公式如下。
遗忘门决定上一时刻有多少信息保留或遗忘在记忆细胞ct,即:
ft=σ(Wfhht-1+Wfxxt+bf)
输入门决定当前时刻有多少信息保留在记忆细胞ct,即:
it=σ(Witht-1+Wixxt+bi)
当前时刻的候选记忆细胞单元ct’,即:
ct'=tanh(Wchht-1+Wcxxt+bc)
记忆细胞由ft点乘上一时刻记忆细胞单元ct-1的值与it点乘当前时刻临时储存单元ct’的值相加更新,即:
ct=ft⊙ct-1+it⊙ct'
输出门决定当前记忆细胞ct有多少输出到当前隐藏层输出值ht,即:
ot=σ(Wohht-1+Woxxt+bo)
将记忆细胞ct经过tanh函数处理后与输出门ot相乘,得到最终输出:
ht=ot⊙tanh(ct)
其中,Wfh、Wih、Woh分别为遗忘门、输入门、输出门与上一时刻隐藏层ht-1的权重矩阵;Wfx、Wix、Wox分别为遗忘门、输入门、输出门与输入向量xt的权重矩阵;bf、bi、bc、bo分别为遗忘门、输入门、记忆细胞、输出门的偏执项;σ代表sigmoid非线性函数;tanh为双曲正切函数;⊙表示向量间点乘。
步骤4):多目标优化:将代理模型作为目标函数,利用非支配排序遗传算法(NSGA-Ⅲ)对目标函数进行优化,生成Pareto前沿。
本实施例使用NSGA-III来对工艺参数进行多目标优化,开发的基于深度学习的数据驱动代理模型作为目标函数。本实施例可以避免计算复杂的数学建模和昂贵的实验。当要构建一个新的种群St,利用基于Pareto支配的非支配排序将原始种群分为若干不同的非支配层(F1,F2…)。从F1开始逐次将各非支配层的解加入到St中。这种分布干扰了最优方案的选择,直到达到形成下一代所需的非支配解的数量。如果要插入的最后一个非支配层FL包含的个体多于需要的,那么必须在FL内根据参考线进行选择。基准线从FL的中心开始并通过基准点。与基准线垂直距离最小的解决方案被优先选择。利用下列公式计算得到工艺参数的帕累托(Pareto)前沿。
具体的,如图4所示,本实施例利用非支配排序遗传算法(NSGA-Ⅲ)对目标函数进行优化的方法如下。
41)通过边界交叉构造权重的方法,基于优化目标数M和各维度目标划分数p,均匀生成H(M,p)个参考点:
42)将第t代种群Pt通过交叉变异操作得到子代种群Qt,每代种群规模为S,将父代种群和子代种群合并形成新种群Rt;
43)将合并种群Rt中的个体进行快速非支配排序,得到若干非支配层F1,F2,F3,...,Fl;
44)将优先级较高的非支配层存入下一代种群Pt+1,临界层为Fτ,满足:
45)得到当前种群的所有个体在每一维目标上的最优值,构成当前种群在每一维目标上的原点。将步骤41)中生成的参考点与原点相连,连接线组成参考向量;计算当前Pt+1中所有个体到参考向量的垂直距离,并与最近的参考向量相关联;
判断迭代次数Gen是否达到设定阈值maxGen;若是,则输出当下种群,得到工艺参数优化解集;若否,则执行步骤42)。
步骤5):选取最优解:利用逼近理想解排序法(TOPSIS)对Pareto解集进行排序,确定最优解。
NSGA-III产生的优化结果信息量较大,同时可供决策者选择的解决方案数目多。从NSGA-III算法得到的Pareto解中选取的代表解可以看作是一个多目标优化问题。为了实现这一目标,利用理想解相似度排序技术对最优有效解进行排序。提出的理想相似度排序方法的基本机制是计算在n维空间中定义的正理想解(PIS)和负理想解(NIS)的每个备选方案之间的距离,其中n为决策问题中准则的数量。最佳的方案与PIS之间的矢量距离最小,与NIS之间的矢量距离最大。
具体的,本实施例利用逼近理想解排序法(TOPSIS)确定最优解的方法如下。
51)建立初始判断矩阵
设评价方案数量为m,所有评价方案构成方案集A={A1,A2,…,Am};每个评价方案有n个评价指标,单一评价方案评判集Ai={Xi1,Xi2,…,Xin},评判指标Xij为第i个评价方案的第j个评价指标值,i∈[1,m],j∈[1,n];构建相应的初始判断矩阵为:
52)建立标准化决策矩阵
评价指标可分为正向指标与负向指标,以最小化为目的的目标通常被认为是负向指标,而以最大化为目的的目标通常被认为是正向指标。将负向指标以倒数方式完成正向化,使所有指标具有同趋势性;将原始数据进行归一化处理以消除量纲:
最终得到标准化决策矩阵Z=(zij)m×n;
53)构建加权评价矩阵
对评价方案的n维评价指标分配权重,构建的评价指标权重矩阵对标准化决策矩阵进行加权得到加权评价矩阵Z`:
其中,Z'表示加权评价矩阵;W表示评价指标的权重矩阵;wj表示第j个评价指标的权重;
54)确定理想解
正理想解Z+是由加权评价矩阵Z`中每列最优值组成:
负理想解Z-是由加权评价矩阵Z`中每列最劣值组成:
55)计算距离尺度与优劣排序
(1)计算距离尺度,确定评价方案与正、负理想解的欧式距离:
其中,Di +表示第i个评价方案到正理想解的欧式距离;Di -表示第i个评价方案到正理想解的欧式距离;z’ij表示加权评价矩阵中的第i个方案得到的第j个指标的值;
(2)计算评价方案与正理想解的接近度:
其中,Ci表示第i个方案与正理想解的接近度;
根据Ci高低进行优劣排序,Ci值越接近1表明评价方案越优。
下面结合具体实例对本实施例融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法进行进一步说明。
1、铣削加工实验环境
为了证明本实施例融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法的可行性,本实施例以铣削加工实验为例,实验所用主要设备如表1所示。
表1.实验主要设备
实验在VCML850数控铣床上进行,数控系统为Siemens 828D,数据处理多目标优化求解采用一台带有4个NVIDIA GeForce RTX 2080Ti显卡的Windows系统计算机。本实施例采用Kistler 9139AA切削力计测量切削过程中的切削力,利用HIOKI PW3337功率仪采集加工过程中的能耗数据,利用Rtec MFT-5000测量工件已加工表面的表面粗糙度。
本实施例中采用的刀具直径为8mm,齿数为4。为了进一步验证多目标优化方法对切削参数优化的有效性,在固定刀具轨迹上进行了优化实验。为了说明在不同优化目标上的有效性和适应性,本实施例进行了两个案例研究。在第一种情况下,有四个相互冲突的目标:生产率(MRR)、表面粗糙度、最大切削力和能耗。在第二种情况下,有三个相互冲突的目标:生产率、最大切削力和能耗,以表面粗糙度(Ra)为约束条件。相应的公式如下:
0<ap≤4
0<ae≤10
0.16≤f*N≤0.5
30≤vc≤100
Ra≤0.5μm
其中,ap为切削深度,ae为切削宽度,f为进给速度,N为刀具的齿数,vc为切削速度。
在本实施例中,由于切削参数和目标具有离散性,所以采用了卷积神经网络(CNN)开发最大切削力、表面粗糙度和能耗的预测模型。数据采集系统采集的整个数据集分为三个部分,用于训练、验证和测试的比例分别为70%、15%和15%。利用平均绝对误差(MAE)、均方根误差(RMSE)和预测精度(PA)对预测模型的性能进行评价。计算公式如下:
在本实施例中,NSGA-Ⅲ算法的突变概率和交叉概率分别设为1/55和1。此外,代数和种群大小都设置为100。在第一种情况下,能耗、最大切削力和表面粗糙度被认为是负向指标,因为目的是使其最小化。生产率(即MRR)被认为是一个正向指标,因为优化目的是使其最大化。本案例在第一种情况中假设所有目标的权重相等,即为(0.25,0.25,0.25,0.25)。在第二种情况下,将能耗和最大切削力作为负指标,将生产率(即MRR)作为正指标。为了说明本实施例提出的方法在不同权重上的有效性,在第二种情况下,所有目标的权重被分为四组。第一组、第二组、第三组、第四组的权重分别设为(0.2,0.2,0.6)、(0.3,0.3,0.4)、(0.5,0.4,0.1)、(0.6,0.2,0.2)。
2、实验结果分析
由表2可以看出,通过基于深度学习的数据驱动代理模型得到的目标函数具有较高的精度和共性,能满足多目标优化对模型精度的要求,可以大大提高优化结果的可靠性。
表2预测性能的统计结果
在第一种情况下利用NSGA-Ⅲ和TOPSIS对工艺参数进行多目标优化和求解,结果如表3。
表3第一种情况多目标优化及排序结果
由表3可知,75号方案效果最好,得分为0.611,工艺参数分别为97.85、11.15×10-2、9.98、3.99,生成的目标值分别为13.94、966.13、1151×10-3、0.97,这些值满足加工过程的约束条件。结果最差的是26号方案,工艺参数分别为95.30、5.19×10-2、1.63、3.89,对应的目标值分别为40.27、432.30、83.18×10-3、0.52,得到的TOPSIS得分为0.373。
表4第二种情况多目标优化及排序结果
由表4可知,不同分组的工艺参数和相应的目标变化较大,这主要是由于赋予不同的权值的造成的。这表明,最佳工艺参数确实随着需求的不同而改变。而且往往很难找到一个合适的解决方案,满足当前的工作要求。但是利用本实施例可以很容易地找出最满足现有工作要求的工艺参数的组合。
第一种情况下优化过程的计算时间汇总如表5所示。在第一种情况下,需要优化四个相互冲突的目标,用于在优化上的时间比第二种情况要长。因此,本实施例用第一种情况的时间来说明本方法的效率。优化总时间为160.4s,模型训练总时间为914.2s,优化总时间为1074.6s,这比传统方法用的时间更短,进一步说明了本实施例方法能解决传统方法优化时间长的问题。
表5第一种情况下优化过程的计算时间
优化过程已在两个案例中成功实施,这证明了本实施例所提出的方法在优化过程中对不同目标的有效性。为了进一步说明这种差异,我们对两种情况下的最佳优化结果进行了比较。对于第二种情况,由于分配给不同目标的权重不同,得到了4个最佳优化结果。通过对比可以发现,第一种情况下的切削力比第一组和第二组在第二种情况下要小,而其他目标则相反。可以发现,虽然得到了最佳的参数组合,但不同的优化目标在不同场合下的响应值差异很大,这进一步说明最佳优化结果往往依赖于不同的工作要求。综上所述,本实施例提出的方法在不同权重和优化目标下的有效性和可行性,证明了提出的方法可以适用于不同的工况,同时操作人员可以很容易地根据他们的实际加工要求调整不同目标的权重,这有助于从复杂的加工参数空间中设计和找到最佳的近似加工参数组合。
总的来说,本实施例所提出的方法的主要优势在于通过使用我们提出的方法,可以在不需要额外的成本或理解物理规律的情况下,轻松地构建高精度和通用的目标函数,更具有成本效益和工业实践。本实施例所提出的优化方法是端到端结构,目标函数的构造、优化计算和Pareto集求解均采用自动化、集成化的方式,人为干扰最小。与传统的优化方法相比,我们提出的方法具有快速的优化速度,只需几分钟就可以完成优化过程。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
Claims (5)
1.一种融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,其特征在于:包括如下步骤:
步骤1):历史数据采集:采集与优化目标和工艺参数相关的数据;
步骤2):数据预处理:对采集得到的历史数据进行预处理,剔除原始数据中的异常值;
步骤3):生成代理模型:利用深度学习方法,根据历史数据生成优化目标的预测模型,并将预测模型转换成代理模型;
步骤4):多目标优化:将代理模型作为目标函数,利用非支配排序遗传算法(NSGA-Ⅲ)对目标函数进行优化,生成Pareto前沿;
步骤5):选取最优解:利用逼近理想解排序法(TOPSIS)对Pareto解集进行排序,确定最优解。
2.根据权利要求1所述融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,其特征在于:所述步骤2)中,利用局部离群因子算法(LOF)对历史数据进行预处理,方法如下:
在样本集D中,定义某样本点o与离其第k远的近邻点间的距离为点o的第k距离dk(o),并定义点o的第k距离邻域Nk(o)为离点o的距离不超过dk(o)的所有样本点,即:
Nk(o)={o'|o'∈D,d(o,o')≤dk(o)}
其中,d(o,o')为样本点o与o'间的距离;若d(o,o')<dk(o),则定义样本点o的可达距离dreach(o,o')为dk(o),否则为d(o,o'),即:
dreach(o,o')=max{d(o,o'),dk(o)}
其中,||Nk(o)||为点o的第k距离邻域内包含的样本点个数;
定义点o的局部可达密度与其第k距离邻域Nk(o)内所有点的局部可达密度之比的平均值为局部离群因子LOFk(o),即:
样本点o的局部可达密度相比其k近邻的局部可达密度越低,则LOFk(o)值越大,说明点o是离群点的可能性也越大。
3.根据权利要求1所述融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,其特征在于:所述步骤3)中,若历史数据为离散数据,则采用卷积神经网络(CNN)生成优化目标的预测模型;若历史数据为时序数据,则采用长短期记忆人工神经网络(LSTM)生成优化目标的预测模型。
4.根据权利要求1所述融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,其特征在于:所述步骤4)中,利用非支配排序遗传算法(NSGA-Ⅲ)对目标函数进行优化的方法为:
41)通过边界交叉构造权重的方法,基于优化目标数M和各维度目标划分数p,均匀生成H(M,p)个参考点:
42)将第t代种群Pt通过交叉变异操作得到子代种群Qt,每代种群规模为S,将父代种群和子代种群合并形成新种群Rt;
43)将合并种群Rt中的个体进行快速非支配排序,得到若干非支配层F1,F2,F3,...,Fl;
44)将优先级较高的非支配层存入下一代种群Pt+1,临界层为Fτ,满足:
45)得到当前种群的所有个体在每一维目标上的最优值,构成当前种群在每一维目标上的原点。将步骤41)中生成的参考点与原点相连,连接线组成参考向量;计算当前Pt+1中所有个体到参考向量的垂直距离,并与最近的参考向量相关联;
判断迭代次数Gen是否达到设定阈值maxGen;若是,则输出当下种群,得到工艺参数优化解集;若否,则执行步骤42)。
5.根据权利要求1所述融合NSGA-Ⅲ和TOPSIS的数据驱动多目标优化方法,其特征在于:所述步骤5)中,利用逼近理想解排序法(TOPSIS)确定最优解的方法为:
51)建立初始判断矩阵
设评价方案数量为m,所有评价方案构成方案集A={A1,A2,…,Am};每个评价方案有n个评价指标,单一评价方案评判集Ai={Xi1,Xi2,…,Xin},评判指标Xij为第i个评价方案的第j个评价指标值,i∈[1,m],j∈[1,n];构建相应的初始判断矩阵为:
52)建立标准化决策矩阵
将负向指标以倒数方式完成正向化,使所有指标具有同趋势性;将原始数据进行归一化处理以消除量纲:
最终得到标准化决策矩阵Z=(zij)m×n;
53)构建加权评价矩阵
对评价方案的n维评价指标分配权重,构建的评价指标权重矩阵对标准化决策矩阵进行加权得到加权评价矩阵Z’:
其中,Z'表示加权评价矩阵;W表示评价指标的权重矩阵;wj表示第j个评价指标的权重;
54)确定理想解
正理想解Z+是由加权评价矩阵Z’中每列最优值组成:
其中,z'j +表示加权评价矩阵Z’中第j列的最优值;
负理想解Z-是由加权评价矩阵Z’中每列最劣值组成:
其中,z'j -表示加权评价矩阵Z’中第j列的最劣值;
55)计算距离尺度与优劣排序
(1)计算距离尺度,确定评价方案与正、负理想解的欧式距离:
其中,Di +表示第i个评价方案到正理想解的欧式距离;Di -表示第i个评价方案到正理想解的欧式距离;z′ij表示加权评价矩阵中的第i个方案得到的第j个指标的值;
(2)计算评价方案与正理想解的接近度:
其中,Ci表示第i个方案与正理想解的接近度;
根据Ci高低进行优劣排序,Ci值越接近1表明评价方案越优。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210751216.6A CN115130749A (zh) | 2022-06-29 | 2022-06-29 | 融合nsga-ⅲ和topsis的数据驱动多目标优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210751216.6A CN115130749A (zh) | 2022-06-29 | 2022-06-29 | 融合nsga-ⅲ和topsis的数据驱动多目标优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115130749A true CN115130749A (zh) | 2022-09-30 |
Family
ID=83379801
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210751216.6A Pending CN115130749A (zh) | 2022-06-29 | 2022-06-29 | 融合nsga-ⅲ和topsis的数据驱动多目标优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115130749A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116362518A (zh) * | 2023-04-04 | 2023-06-30 | 无锡物联网创新中心有限公司 | 一种纺纱配棉方法及相关装置 |
CN117113589A (zh) * | 2023-10-23 | 2023-11-24 | 泉州装备制造研究所 | 一种弱刚性薄壁结构预连接工艺多目标智能优化方法 |
CN117709168A (zh) * | 2024-02-02 | 2024-03-15 | 北京理工大学深圳汽车研究院(电动车辆国家工程实验室深圳研究院) | 一种强制风冷散热系统的优化设计方法及其装置 |
CN117993893A (zh) * | 2024-04-03 | 2024-05-07 | 华东交通大学 | 一种铁路接触网维修方案生成方法 |
CN118230838A (zh) * | 2024-03-18 | 2024-06-21 | 重庆大学 | 一种融合甲醇生产系统机理的多目标自适应智能优化算法 |
-
2022
- 2022-06-29 CN CN202210751216.6A patent/CN115130749A/zh active Pending
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116362518A (zh) * | 2023-04-04 | 2023-06-30 | 无锡物联网创新中心有限公司 | 一种纺纱配棉方法及相关装置 |
CN116362518B (zh) * | 2023-04-04 | 2023-12-19 | 无锡物联网创新中心有限公司 | 一种纺纱配棉方法及相关装置 |
CN117113589A (zh) * | 2023-10-23 | 2023-11-24 | 泉州装备制造研究所 | 一种弱刚性薄壁结构预连接工艺多目标智能优化方法 |
CN117113589B (zh) * | 2023-10-23 | 2024-02-13 | 泉州装备制造研究所 | 一种弱刚性薄壁结构预连接工艺多目标智能优化方法 |
CN117709168A (zh) * | 2024-02-02 | 2024-03-15 | 北京理工大学深圳汽车研究院(电动车辆国家工程实验室深圳研究院) | 一种强制风冷散热系统的优化设计方法及其装置 |
CN117709168B (zh) * | 2024-02-02 | 2024-04-30 | 北京理工大学深圳汽车研究院(电动车辆国家工程实验室深圳研究院) | 一种强制风冷散热系统的优化设计方法及其装置 |
CN118230838A (zh) * | 2024-03-18 | 2024-06-21 | 重庆大学 | 一种融合甲醇生产系统机理的多目标自适应智能优化算法 |
CN117993893A (zh) * | 2024-04-03 | 2024-05-07 | 华东交通大学 | 一种铁路接触网维修方案生成方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115130749A (zh) | 融合nsga-ⅲ和topsis的数据驱动多目标优化方法 | |
US20220011740A9 (en) | Energy consumption prediction system and method based on the decision tree for CNC lathe turning | |
CN112966954A (zh) | 一种基于时间卷积网络的防洪调度方案优选方法 | |
CN106446478B (zh) | 一种切削加工工艺优选方法 | |
CN113780153B (zh) | 一种刀具磨损监测和预测方法 | |
CN114254915B (zh) | 轴类零件全流程加工质量合格状态决策及优化方法 | |
CN105975701A (zh) | 一种基于混合模糊模型的并行调度拆卸路径生成方法 | |
CN114611379A (zh) | 一种基于数据驱动的加工工艺节能规划方法 | |
CN115907195A (zh) | 一种光伏发电功率预测方法、系统、电子设备及介质 | |
CN113947332B (zh) | 一种地下工程综合保障能力的评估方法及系统 | |
Wang et al. | An intelligent process parameters determination method based on multi-algorithm fusion: a case study in five-axis milling | |
CN116187836A (zh) | 一种电力物资质量评价及抽检方法 | |
CN112990601B (zh) | 基于数据挖掘的蜗轮加工精度自愈系统及方法 | |
CN113450013A (zh) | 基于改进nsga-ⅲ算法求解车间节能调度问题的方法 | |
CN117458480A (zh) | 基于改进lof的光伏发电功率短期预测方法及系统 | |
Guan et al. | Machining scheme selection of digital manufacturing based on genetic algorithm and AHP | |
CN115081551A (zh) | 基于K-Means聚类和优化RVM线损模型建立方法及系统 | |
CN115034070A (zh) | 一种基于多目标优化和vikor方法的复杂机械产品的选择装配优化和决策方法 | |
CN112365168A (zh) | 一种基于主分量分析评价环境空气质量的方法 | |
CN111626465A (zh) | 一种新能源功率短期区间预测方法及系统 | |
CN114297802B (zh) | 一种薄板激光切割工艺参数的多目标优化方法 | |
CN116736796B (zh) | 一种钢材尺寸的高精度工控方法 | |
CN113988391B (zh) | 一种基于相似日理论和神经网络的光伏发电功率预测方法 | |
CN114943304B (zh) | 一种基于贝叶斯的3c电子产品装配误差传递方法 | |
Li et al. | Machining process energy consumption modelling using response surface methodology and neural network |
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 |