CN103336869B - 一种基于高斯过程联立mimo模型的多目标优化方法 - Google Patents
一种基于高斯过程联立mimo模型的多目标优化方法 Download PDFInfo
- Publication number
- CN103336869B CN103336869B CN201310281958.8A CN201310281958A CN103336869B CN 103336869 B CN103336869 B CN 103336869B CN 201310281958 A CN201310281958 A CN 201310281958A CN 103336869 B CN103336869 B CN 103336869B
- Authority
- CN
- China
- Prior art keywords
- gaussian process
- model
- matrix
- formula
- phi
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 150
- 230000008569 process Effects 0.000 title claims abstract description 88
- 238000013461 design Methods 0.000 claims abstract description 29
- 238000005457 optimization Methods 0.000 claims abstract description 29
- 230000004044 response Effects 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000013401 experimental design Methods 0.000 claims abstract description 6
- 230000002068 genetic effect Effects 0.000 claims abstract description 6
- 230000035772 mutation Effects 0.000 claims abstract description 6
- 238000002474 experimental method Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 76
- 238000005070 sampling Methods 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 5
- 238000003050 experimental design method Methods 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000002372 labelling Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000013507 mapping Methods 0.000 abstract description 2
- 230000009467 reduction Effects 0.000 abstract description 2
- 239000003795 chemical substances by application Substances 0.000 description 11
- 238000000354 decomposition reaction Methods 0.000 description 7
- 239000000463 material Substances 0.000 description 7
- 238000013178 mathematical model Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 235000015429 Mirabilis expansa Nutrition 0.000 description 1
- 244000294411 Mirabilis expansa Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229960000074 biopharmaceutical Drugs 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004941 influx Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 235000013536 miso Nutrition 0.000 description 1
- 230000004952 protein activity Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于高斯过程联立MIMO模型的多目标优化方法,该方法通过实验设计方法得到样本,利用高斯过程联立MIMO模型近似建立设计变量与待考察响应之间的映射关系,采用高斯变异混合遗传算法、蚁群算法等对高斯过程联立MIMO模型进行多目标优化,得到关于设计变量组合的Pareto前沿,进而判断Pareto前沿中解集的分布质量,最后根据具体要求选择Pareto前沿中某一设计变量组合进行高精度分析求解,当得到的结果满意,则进行物理实验。本发明将实验设计、高精度分析求解、代理模型技术以及优化算法集成应用于优化设计中,能够大大缩减优化设计的耗时和计算成本,显著提高工作效率。
Description
技术领域
本发明涉及一种代理模型技术结合优化算法的多目标优化方法,尤其涉及一种高斯过程联立建模结合优化算法的多目标优化设计方法。
背景技术
代理模型技术的发展为优化技术的引入提供了有效的桥梁,通过先进实验设计方法得到样本数据,利用代理模型技术近似建立设计变量与待考察响应之间的映射关系,最后运用优化算法求解该关系模型得到最优设计变量组合。这种将实验设计、高精度分析求解、代理模型技术以及优化算法集成应用于优化设计中的方法,不仅能够极大减少变量修改中的人为因素,而且能够在复杂的设计方案中得到尽可能完善的设计变量组合。
现有基于代理模型技术的多目标优化方法中,有以下几种:
(1)基于Kriging模型的多目标优化方法,如KaymazI在:“Applicationofkrigingmethodtostructuralreliabilityproblems(Kriging模型在结构可靠性问题中的应用)”(StructuralSafety,2005,27(2):133-151)(结构安全性,2005年27期第2卷,133-151页)进行优化时,核函数参数的寻优过程会造成大量的计算成本,这在高维的各向异性问题中尤为明显,这一缺陷导致建立Kriging模型所耗费的时间相对其他建模方法而言要多很多。
(2)基于传统独立高斯过程回归模型的多目标优化方法,如JinYuan在:“Reliablemulti-objectiveoptimizationofhigh-speedWEDMprocessbasedonGaussianprocessregression(基于高斯过程回归的高速线切割加工可靠性多目标优化)”(InternationalJournalofMachineTools&Manufacture,2008(48):47–60)(机床与制造国际期刊,2008年48期,47-60页)进行多目标优化时,实际上是采用传统高斯过程回归模型分别独立地对每一个输出端建立高斯过程模型,在对同一组未知输入向量的各个响应进行拟合或预测时,不同输出端对应的各模型所利用的估计参数也是不同的,因此,传统高斯过程模型本质上可以看作一组相对独立的多输入—单输出(multipleinputsandsingleoutput,MISO)模型,无法充分利用数据信息。
(3)基于限定性相关矩阵的代理模型的多目标优化方法,如McMillianNJ在:“AnalysisofProteinActivityDatabyGaussianStochasticProcessModels(基于高斯随机过程模型的蛋白质活动分析)”(JournalofBiopharmaceuticalStatistics,1999(9):145-160)(生物制药统计学期刊,1999年9期,145-160页)预测结果中,限定性相关矩阵赋予类别输入变量某些特定的假设参数化关系,虽然能够大大简化相关计算的复杂程度,但是限定性相关矩阵通常只能用于特定的过程或者系统,因此缺乏捕捉各输出端之间不同类型相关性的柔性,如乘积相关矩阵只能捕捉类别输入变量水平间的正相关性,而缺乏捕捉其负相关性的能力。
综上所述,现有代理模型多目标优化方法中,不足主要体现在以下四方面:1、存在计算耗时长的问题;2、在假设所有的输入信息仅有量化了的设计变量(数值输入变量)的前提下,通过计算未知点与已知点的相对位置关系,确定未知点在数值输入变量空间中的具体位置(数值输入变量的空间相关性),继而进行模型的拟合或预测,然而事实上,需要考察的输出端个数也是已知的,且每个输出端在性质上是分类的(类别输入变量)并能决定代理模型的一些重要性质,如果把系统中的已知条件都看作建模时的输入信息,那么传统多目标建模方法没有综合考虑这两种类型的输入变量,因此常常未能充分利用数据信息;3、许多工程问题的多目标模型中的输出端通常来说都承载了彼此之间的一些相关性,传统多目标建模方法忽略了在同一组数值输入向量下,不同输出端的响应之间存在的相关性(类别输入变量的交叉相关性),由于没有对类别输入变量的交叉相关性进行建模,因此无法利用该相关性来提高模型对未知点的拟合或预测精度;4、缺乏捕捉各输出端之间不同类型相关性的柔性。
发明内容
本发明的目的在于提出一种基于高斯过程联立MIMO(multipleinputsandmultipleoutputs)模型的多目标优化方法,凭借高斯过程联立MIMO模型的模型相关参数较少、易实现、样本规模较小的优势,解决优化计算耗时长、系统信息利用不充分、多目标代理模型无法利用类别输入变量的交叉相关性来提高模型对未知点拟合或预测的精度以及无法捕捉输出端之间不同类型相关性的问题。
本发明实现上述目的所采取的技术方案是:一种基于高斯过程联立MIMO模型的多目标优化方法,包括实验设计方法抽样,标准高斯过程建模及多目标优化算法寻优与验证,具体按照如下步骤进行:
步骤1,根据用户给定的初始条件,通过增强型平移传播拉丁超立方ETPLHD抽样方法进行实验设计,得到设计变量组合X=[x1,x2,…,xN]T,并对设计变量组合分别进行高精度分析求解,得到相应的响应值y=[y1,y2…,yN]T,与设计变量组合一并组成模型的训练数据集(X,y);
步骤2,以步骤1得到的训练数据集,建立基于非限定性相关矩阵的高斯过程联立MIMO模型,确定基于非限定性相关矩阵的高斯过程联立MIMO模型中的协方差函数,该协方差函数的计算公式如下公式(1):
式中h=(xt,zc)t,其中xt是数值输入变量,z为类别输入变量,并以c=1,…,m作为标记,c表示模型输出端的个数,是输出端zr和zs之间的交叉关联度大小,定义非限定性相关矩阵T={τr,s},设θ=[w1,…,wD,v0,v1]T是指数平方平稳协方差函数的超参数,v0是白噪音估计方差,v1是纵坐标上的估计方差,D是设计变量的维数,δpq是Kronecker算子;
并对所述的协方差函数中的模型参数θ、v0和v1进行估计和优化,估计和优化的方法是通过最大化边界似然函数法,转化为求解最小化边界log-似然函数;
步骤3,采用高斯变异混合遗传算法GGA-h、多目标粒子群优化算法和蚁群算法对步骤2所建立的基于非限定性相关矩阵的高斯过程联立MIMO模型进行多目标优化;
步骤4,判断步骤3中多目标优化得到的Pareto前沿中解的分布是否均匀,解的个数是否足够以及边界处的解是否接近单目标极值而具有良好拓展性,当前沿质量达标时,保存该Pareto前沿,否则返回重新设定多目标优化的算法参数;
步骤5,在步骤4得到符合要求的Pareto前沿的前提下,根据用户具体要求选择Pareto前沿中某一设计变量组合进行高精度分析求解,当得到的结果满意时,则为优化结果,否则返回步骤1中重新进行ETPLHD抽样实验设计。
步骤2所述的确定基于非限定性相关矩阵的高斯过程联立MIMO模型中的协方差函数,包括有效定义的数值输入变量空间相关性和类别输入变量交叉相关性,首先根据用户要求的具体设计变量维数D赋予相应超参数的初值个数为D+2,即增加v0和v1,然后根据具体输出端个数c确定非限定性相关矩阵T的规模,即c×c,并赋予初值,具体方法为该矩阵主对角线元素的自相关系数为1,其余矩阵元素为0.01,高斯过程联立MIMO模型相关参数个数为(D+2)+c(c-1)/2。
所述的非限定性相关矩阵T,其超球面分解参数化步骤如下:
步骤A,对一个非限定性相关矩阵T进行乔里斯基分解(Cholesky-typedecomposition),得到公式(2):
T=LLt(2)
式中L={lr,s}是一个主对角线元素严格正定的下三角矩阵;
步骤B,每个行向量(lr,1,…,lr,r)根据一个r-维单位超球面体的曲面坐标点来定义,根据球坐标系,当r=1时,l1,1=1,r=2,…,m,得到公式(3):
式中φr,s∈(0,π),所有的φr,s记作Φ,且φr,s的取值约束在(0,π)之间,即将优化问题中的矩阵正定性不等式约束参数化转变到易于计算的边界约束。
步骤2所述的将模型参数θ、v0和v1进行估计和优化的具体步骤如下:
步骤a,将高斯过程联立MIMO模型中的协方差矩阵Km表示成Kronecker形式K为标准高斯过程协方差矩阵,根据克罗内克积的特性,Km转变为一个Nm×Nm的分块矩阵;
步骤b,按照公式(4)、(5)和(6)对协方差矩阵Km进行一系列矩阵计算和变换;
从公式(6)可知,非限定性相关矩阵T和标准高斯过程协方差矩阵K是相互独立的,而T经过超球面分解与Φ一一对应,通过公式(7)和(8)可分别得到参数θ和Φ。
与现有技术相比,本发明具备的有益效果是:提出一种简单灵活的基于高斯过程联立MIMO模型的多目标优化方法,以标准高斯过程回归模型的重要概念和模型的公式实现为基础,分别以距离基函数定义数值输入变量的空间相关性,以相关矩阵定义类别输入变量的交叉相关性,充分利用各类过程或系统信息,构建基于非限定性相关矩阵的高斯过程联立MIMO模型中协方差函数的一般方法;进而采用超球面分解参数化了非限定性相关矩阵,保持应用柔性的同时,将优化问题中复杂的正定性约束转化为较为简单的边界约束,大大降低计算复杂度,显著减少计算耗时。
附图说明
图1为本发明具体实施例1中IG建模方法得到的Pareto前沿。
图2为本发明具体实施例1中SG建模方法得到的Pareto前沿。
图3为本发明具体实施例1中解析式传回法得到的Pareto前沿。
图4为本发明具体实施例中二维超球面分解示意图。
图5为本发明具体实施例中三维超球面分解示意图。
图6为本发明具体实施例二中汽车桥壳后盖剖视图。
图7为本发明具体实施例二中汽车桥壳后盖正二测视图。
图8为本发明具体实施例二中优化流程图。
图9为本发明具体实施例二中本基于高斯过程联立MIMO模型的多目标优化方法得到的Pareto前沿。
图10为本发明具体实施例二中优化后的成形极限图。
图11为本发明具体实施例二中材料流入情况图。
图12为本发明具体实施例二中料厚的样件值与模拟值。
具体实施方式
下面结合实施例及其附图详细说明基于高斯过程联立MIMO模型的多目标优化方法的技术方案,但本发明的具体实施方式不局限于下述的实施例。
在实施例1中将通过一个数值算例来证明本发明所提出的基于高斯过程联立MIMO模型的多目标优化方法的有效性和可行性。为了达到比较的目的,将考虑下述两种建模方法:
A)独立高斯过程模型(individualGaussianprocess),记作IG。该方法实际上就是通过标准高斯过程模型分别对每一个需要考察的目标建模。
B)联立高斯过程模型(simultaneousGaussianprocess),记作SG。即通过本发明所述的高斯过程联立MIMO模型对目标进行多目标联立建模。
实施例1
解决多目标问题的数值算例
一种基于高斯过程联立MIMO模型的多目标优化方法,包括实验设计方法抽样,标准高斯过程联立MIMO模型及多目标优化算法寻优,包括如下步骤:
步骤1,根据用户给定的数值多目标问题,其数学描述如公式(9),通过拉丁超立方实验设计方法进行实验设计,在变量x∈[0,1]中抽样得到20个样本点,并代入公式(9)的解析式f1,g(x1,x2)中,得到相应的响应值,与变量x的样本集一并组成模型的训练数据集;
f1=x1
f2=g(x1,x2)h(x1,x2)
g(x1,x2)=11+x2 2-10cos(2πx2)(9)
(x∈(0,1))
步骤2,以步骤1得到的训练数据集,分别通过IG和SG建模方法对本实施例中的解析式f1,g(x1,x2)构建代理模型,其中包括分别确定模型中的协方差函数,IG和SG建模方法的协方差函数分别如公式(10)和公式(11)所示,并对两个协方差函数中的模型相关参数估计和优化,估计和优化的方法是通过最大化边界似然函数法,转化为求解最小化边界log-似然函数;
式中h=(xt,z2)t,其中xt是数值输入变量,z为类别输入变量,并以c=2作为标记,c表示模型输出端的个数,是输出端z1和z2之间的交叉关联度大小,定义T={τ1,2},设θ=[w1,v0,v1]T是指数平方平稳协方差函数的超参数,v0是白噪音估计方差,v1是纵坐标上的估计方差,变量的维数为1,δpq是Kronecker算子;
步骤3,采用第二代遗传算法NSGAII分别对步骤2中IG和SG建模方法所建立的代理模型进行多目标优化,算法参数设定如下:种群大小为50,交叉概率为0.8,变异概率为0.3,迭代次数为200,如图1为IG建模方法得到的Pareto前沿,图2为SG建模方法得到的Pareto前沿;
步骤4,在使用第二代遗传算法解决本实施例中的多目标问题时,算法寻优过程中,是将每一次迭代产生的新种群传回到公式(9)的解析式f1,g(x1,x2)中,根据响应值进行相应的适应度初值计算并进行非劣排序(本算例中称为解析式传回法),使用解析式传回法得到的标准Pareto前沿,如图3。对比图1、图2、图3,可知步骤3中应用IG建模方法时得到的图1中Pareto前沿在边界处出现了很大的误差,而且解的分布性也不均匀;而图2中应用SG建模方法时得到的Pareto前沿边界附近的解比较准确,虽然边界附近的分布性稍稍差于图3,但是前沿整体的吻合程度要比IG高得多,且解的分布更加均匀,在Pareto前沿中最常使用的前沿中段部分,其解的分布性甚至优于图3;
步骤5,由步骤4得到的Pareto前沿结果分析中,可以证明本发明提出的基于高斯过程联立MIMO模型的多目标优化方法比基于标准高斯过程建模方法具有更高的预测精度;且能够充分利用数据中的信息;并在解决多目标问题中有着更加良好的表现。
在步骤2中,SG建模方法的高斯过程联立MIMO模型中的协方差函数,此时的具体变量维数为1赋予相应超参数的初值个数为1+2,即增加v0和v1,且具体输出端个数2,可确定非限定性相关矩阵T2的规模,即2×2,并赋予初值,如公式(12),具体方法为该矩阵主对角线元素的自相关系数为1,其余矩阵元素为0.01,高斯过程联立MIMO模型相关参数个数为4。
非限定性相关矩阵T2,其超球面分解参数化步骤如下:
步骤A,对一个非限定性相关矩阵T进行乔里斯基分解(Cholesky-typedecomposition),得到公式(13)
式中L={lr,s}是一个主对角线元素严格正定的下三角矩阵;
步骤B,每个行向量(l2,1,l2,2)根据一个2-维单位超球面体的曲面坐标点来定义,此时输出端为2维,则其分解结果是一个半球面,如图4所示;若输出端为3维,则其分解结果是一个半球体,如图5所示;若输出端超过3维,则其分解结果为超球面,结果无法以图形形式表示。本实施例中根据如下球坐标系,当r=1,l1,1=1,r=2,得到公式(14)
式中φr,s∈(0,π),所有的φr,s记作Φ,且φr,s的取值约束在(0,π)之间,即将优化问题中的矩阵正定性不等式约束参数化转变到易于计算的边界约束。
在步骤2所述的模型相关参数估计和优化步骤如下:
步骤a,将高斯过程联立MIMO模型中的协方差矩阵K2表示成Kronecker形式K为标准高斯过程协方差矩阵,根据克罗内克积的特性,K2转变为一个40×40的分块矩阵;
步骤b,对K2进行一系列矩阵计算和变换,如公式(15)、(16)和(17)所示:
从公式(17)可知,非限定性相关矩阵T和标准高斯过程协方差矩阵K是相互独立的,而T经过超球面分解与Φ一一对应,则有公式(18)和(19):
因此模型相关参数参数θ和Φ可以通过(18)和(19)两个公式在估计和优化过程中分别得到。
实施例2
基于高斯过程联立MIMO模型的多目标优化方法在板料拉深中的应用
采用基于高斯过程联立MIMO模型的多目标优化方法优化某汽车桥壳后盖,如图6和图7所示,拉深工艺参数多目标优化,其优化过程流程图如图8所示,包括如下步骤:
步骤1,根据客户制订的工艺文件,选定待优化的成形工艺参数,并确定其取值范围如下:80≤压边力BHF(kN)≤500,275≤坯料直径D(mm)≤325,0.10≤板料与凹模的静摩擦系数μ1≤0.15,0.10≤板料与压边圈的静摩擦系数μ2≤0.15,通过增强型平移传播拉丁超立方(ETPLHD)抽样方法进行实验设计,得到30组工艺参数组合,并利用板料成形CAE软件根据每一组样本所对应的工艺参数组合分别建立有限元模型进行数值模拟,根据成形极限图FLD,定义得到拉裂与起皱目标函数,如公式(20)(21)所示,再根据每一组工艺参数组合下的数值模拟结果数据分别计算目标函数值DC和DW,各组响应值与工艺参数组合一并组成高斯过程联立MIMO模型的训练数据集,见表1;
表1ETPLHD抽样样本及其响应
步骤2,以步骤1得到的训练数据集,建立基于非限定性相关矩阵的高斯过程联立MIMO模型,其中包括确定基于非限定性相关矩阵的高斯过程联立MIMO模型中的协方差函数,如公式(22)所示,并对该协方差函数中的模型相关参数估计和优化,估计和优化的方法是通过最大化边界似然函数法,转化为求解最小化边界log-似然函数;
式中h=(xt,z2)t,其中xt是数值输入变量,z为类别输入变量,并以c=2作为标记,c表示模型输出端的个数,是输出端z1和z2之间的交叉关联度大小,定义T={τ1,2},设θ=[w1,…,w4,v0,v1]T是指数平方平稳协方差函数的超参数,v0是白噪音估计方差,v1是纵坐标上的估计方差,设计变量的维数为4,δpq是Kronecker算子;
步骤3,采用高斯变异混合遗传算法(GGA-h)算法参数设定如下:种群大小为50,交叉概率为0.9,变异概率为0.1,迭代次数为300,如图9所示,为基于高斯过程联立MIMO模型的多目标优化方法得到的Pareto前沿;
步骤4,由步骤3中的Pareto前沿,解的分布均匀情况良好,边界处的解比较接近单目标极值,说明解集有良好的拓展性。与此同时,能够知道本例中拉裂和起皱存在非线性关系,而交叉关联度则说明两者存在负相关,但不是完全负相关,即并非拉裂目标函数值减小越快,起皱目标函数值就增大越快,当受压达到一定程度时,板料成形失稳,由面内应变转变为弯曲变形时,材料流动受阻,也会导致拉裂缺陷增加,这也与实际情况相符。为了验证本发明中的优化方法得到的结果的准确性,表2为图9所示Pareto前沿中的某3组工艺参数组合的模型算法预测值和计算机辅助工程(CAE)验证值,由对比数据可知模型算法预测值与CAE验证值十分接近,虽有少许偏差,但其变化趋势一致,因此优化结果可以指导该桥壳后盖的生产,保存该Pareto前沿;
表2对比数据验证
步骤5,在步骤4得到了符合要求的Pareto前沿的前提下,表3是从该Pareto前沿中选择的与初始样本点中选出的拉裂目标函数值相近的3组数据作对比,Pareto前沿中非劣解的起皱目标函数值均比初始样本中的该值有不同程度的降低。当从该Pareto前沿和初始样本中选择起皱目标函数值相近的数据时,也可得到相同的结论。这表明本发明的优化方法能够同时优化拉裂和起皱两个质量指标。根据用户制订的工艺文件,选取如图9所示的P点作为满意解,其DC和DW值分别为0.01307和4.518,设计变量BHF、D、μ1和μ2所对应的值分别为114.7kN,286mm,0.127和0.119,利用该工艺参数组合进行数值模拟,根据得到的成形极限图图10所示,能够看出零件的拉裂和起皱趋势都很小。材料的流入情况如图11所示,可以看出拉深成形过程中,料流状况良好,材料流入均匀,最大最小流入量差值为3.22mm。根据优化的如图9所示的P点工艺参数组合,做相应的微调,使其满足实际生产设置要求,样件成形质量良好。对样件在对称轴处进行线切割,量取样件对称轴上圆角、直边等特征起始位置以及中段均分点处24个点的料厚值,并在数值模拟结果中提取相应位置的料厚值,样件值与模拟值如图12,两者的分布和变化趋势基本是一致的,同时说明了样件的料厚分布情况良好。根据用户制订的检测计划,对试制的5个样件进行尺寸测量,结果见表4。
表3拉裂目标函数值DC相当时,起皱目标函数DW优化值和样本值对比数据
表4拉深成形样件检验报告
在步骤2中,本发明的优化方法中的高斯过程联立MIMO模型中的协方差函数,此时的具体设计变量维数为4,赋予相应超参数的初值个数为4+2,即增加v0和v1,且具体输出端个数2,可确定非限定性相关矩阵T2的规模,即2×2,并赋予初值,如公式(12),具体方法为该矩阵主对角线元素的自相关系数为1,其余矩阵元素为0.01,高斯过程联立MIMO模型相关参数个数为7。
非限定性相关矩阵T2,其超球面分解参数化步骤如下:
步骤A,对一个非限定性相关矩阵T进行乔里斯基分解(Cholesky-typedecomposition),得到公式(13),式中L={lr,s}是一个主对角线元素严格正定的下三角矩阵;
步骤B,每个行向量(l2,1,l2,2)根据一个2-维单位超球面体的曲面坐标点来定义,此时输出端为2维,则其分解结果是一个半球面,如图4所示;若输出端为3维,则其分解结果是一个半球体,如图5所示;若输出端超过3维,则其分解结果为超球面,结果无法以图形形式表示。本实施例中根据如下球坐标系,当r=1,l1,1=1,r=2,得到公式(14),式中φr,s∈(0,π),所有的φr,s记作Φ,且φr,s的取值约束在(0,π)之间,即将优化问题中的矩阵正定性不等式约束参数化转变到易于计算的边界约束。
在步骤2所述的模型相关参数估计和优化步骤如下:
步骤a,将高斯过程联立MIMO模型中的协方差矩阵K2表示成Kronecker形式K为标准高斯过程协方差矩阵,根据克罗内克积的特性,K2转变为一个40×40的分块矩阵;
步骤b,对K2进行一系列矩阵计算和变换,如公式(15)(16)(17)所示,从公式(17)可知,非限定性相关矩阵T和标准高斯过程协方差矩阵K是相互独立的,而T经过超球面分解与Φ一一对应,则有公式(18)(19),因此模型相关参数参数θ和Φ可以通过(18)(19)两个公式在估计和优化过程中分别得到。
由此可见,高斯过程模型有严格的统计学理论基础,对小样本、非线性以及高维问题有着良好的适应性,且实现过程相对简单。因此本发明在标准高斯过程回归模型的基础上,提出了一种简单灵活的基于高斯过程联立MIMO模型的多目标优化方法。分别以距离基函数定义了数值输入变量的空间相关性,以相关矩阵定义了类别输入变量的交叉相关性,构建了基于非限定性相关矩阵的高斯过程联立MIMO模型中协方差函数的一般方法;进而采用超球面参数化分解方法构建非限定性相关矩阵,保持了应用柔性的同时,将优化问题中复杂的正定性约束转化为较为简单的边界约束,大大降低计算复杂度。通过两个实施例,证明基于高斯过程联立MIMO模型的多目标优化方法在模型精度、输出端相关性捕捉以及在解决多目标优化问题中具有非常明显的优势。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于高斯过程联立MIMO模型的多目标优化方法,包括实验设计方法抽样,标准高斯过程建模及多目标优化算法寻优与验证,其特征在于,具体按照如下步骤进行:
步骤1,根据用户给定的初始条件,通过增强型平移传播拉丁超立方ETPLHD抽样方法进行实验设计,得到设计变量组合X=[x1,x2,…,xN]T,并对设计变量组合分别进行高精度分析求解,得到相应的响应值y=[y1,y2…,yN]T,与设计变量组合一并组成模型的训练数据集(X,y);
步骤2,以步骤1得到的训练数据集,建立基于非限定性相关矩阵的高斯过程联立MIMO模型,确定基于非限定性相关矩阵的高斯过程联立MIMO模型中的协方差函数,该协方差函数的计算公式如下公式(1):
式中h=(xt,zc)t,其中xt是数值输入变量,z为类别输入变量,并以c=1,…,m作为标记,c表示模型输出端的个数,是输出端zr和zs之间的交叉关联度大小,定义非限定性相关矩阵T={τr,s},设θ=[w1,…,wD,v0,v1]T是指数平方平稳协方差函数的超参数,v0是白噪音估计方差,v1是纵坐标上的估计方差,D是设计变量的维数,δpq是Kronecker算子;
并对所述的协方差函数中的模型参数θ、v0和v1进行估计和优化,估计和优化的方法是通过最大化边界似然函数法,转化为求解最小化边界log-似然函数;
步骤3,采用高斯变异混合遗传算法GGA-h、多目标粒子群优化算法和蚁群算法对步骤2所建立的基于非限定性相关矩阵的高斯过程联立MIMO模型进行多目标优化;
步骤4,判断步骤3中多目标优化得到的Pareto前沿中解的分布是否均匀,解的个数是否足够以及边界处的解是否接近单目标极值而具有良好拓展性,当前沿质量达标时,保存该Pareto前沿,否则返回重新设定多目标优化的算法参数;
步骤5,在步骤4得到符合要求的Pareto前沿的前提下,根据用户具体要求选择Pareto前沿中某一设计变量组合进行高精度分析求解,当得到的结果满意时,则为优化结果,否则返回步骤1中重新进行ETPLHD抽样实验设计。
2.根据权利要求1所述的基于高斯过程联立MIMO模型的多目标优化方法,其特征在于,步骤2所述的确定基于非限定性相关矩阵的高斯过程联立MIMO模型中的协方差函数,包括有效定义的数值输入变量空间相关性和类别输入变量交叉相关性,首先根据用户要求的具体设计变量维数D赋予相应超参数的初值个数为D+2,即增加v0和v1,然后根据具体输出端个数c确定非限定性相关矩阵T的规模,即c×c,并赋予初值,具体方法为该矩阵主对角线元素的自相关系数为1,其余矩阵元素为0.01,高斯过程联立MIMO模型相关参数个数为(D+2)+c(c-1)/2。
3.根据权利要求2所述的基于高斯过程联立MIMO模型的多目标优化方法,其特征在于,所述的非限定性相关矩阵T,其超球面分解参数化步骤如下:
步骤A,对一个非限定性相关矩阵T进行乔里斯基分解,得到公式(2):
T=LLt(2)
式中L={lr,s}是一个主对角线元素严格正定的下三角矩阵;
步骤B,每个行向量(lr,1,…,lr,r)根据一个r-维单位超球面体的曲面坐标点来定义,根据球坐标系,当r=1时,l1,1=1,r=2,…,m,得到公式(3):
式中φr,s∈(0,π),所有的φr,s记作Φ,且φr,s的取值约束在(0,π)之间,即将优化问题中的矩阵正定性不等式约束参数化转变到易于计算的边界约束。
4.根据权利要求1或2所述的基于高斯过程联立MIMO模型的多目标优化方法,其特征在于,步骤2所述的将模型参数θ、v0和v1进行估计和优化的具体步骤如下:
步骤a,将高斯过程联立MIMO模型中的协方差矩阵Km表示成Kronecker形式K为标准高斯过程协方差矩阵,根据克罗内克积的特性,Km转变为一个Nm×Nm的分块矩阵;
步骤b,按照公式(4)、(5)和(6)对协方差矩阵Km进行一系列矩阵计算和变换;
从公式(6)可知,非限定性相关矩阵T和标准高斯过程协方差矩阵K是相互独立的,而T经过超球面分解与Φ一一对应,通过公式(7)和(8)可分别得到参数θ和Φ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310281958.8A CN103336869B (zh) | 2013-07-05 | 2013-07-05 | 一种基于高斯过程联立mimo模型的多目标优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310281958.8A CN103336869B (zh) | 2013-07-05 | 2013-07-05 | 一种基于高斯过程联立mimo模型的多目标优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103336869A CN103336869A (zh) | 2013-10-02 |
CN103336869B true CN103336869B (zh) | 2016-07-06 |
Family
ID=49245033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310281958.8A Expired - Fee Related CN103336869B (zh) | 2013-07-05 | 2013-07-05 | 一种基于高斯过程联立mimo模型的多目标优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103336869B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103488851B (zh) * | 2013-10-15 | 2016-07-06 | 浙江大学 | 基于几何结构信息的多目标优化方法 |
CN104008248B (zh) * | 2014-06-04 | 2017-10-10 | 广西大学 | 基于高斯过程的注塑成形工艺稳健设计与容差设计方法 |
JP6371870B2 (ja) * | 2014-06-30 | 2018-08-08 | アマゾン・テクノロジーズ・インコーポレーテッド | 機械学習サービス |
CN107783998A (zh) * | 2016-08-26 | 2018-03-09 | 华为技术有限公司 | 一种数据处理的方法以及装置 |
CN106682772A (zh) * | 2016-12-14 | 2017-05-17 | 贵州电网有限责任公司电力科学研究院 | 一种基于多维云模型的输电线路故障率综合预测方法 |
CN107038303B (zh) * | 2017-04-13 | 2020-08-18 | 西北工业大学 | 用于机械可靠性分析与设计的基于代理模型的双层实验设计方法 |
CN107844855A (zh) * | 2017-10-23 | 2018-03-27 | 深圳大学 | 一种基于昂贵多目标优化问题的优化方法及系统 |
CN107844659B (zh) * | 2017-11-16 | 2021-05-14 | 华东理工大学 | 水煤浆气化过程的代理模型建模方法 |
CN108197738B (zh) * | 2017-12-29 | 2020-09-15 | 东华大学 | 一种涤纶长丝熔体输送过程的工艺参数优化方法 |
CN111382210B (zh) * | 2018-12-27 | 2023-11-10 | 中国移动通信集团山西有限公司 | 一种分类方法、装置及设备 |
CN110046378B (zh) * | 2019-02-28 | 2022-09-13 | 昆明理工大学 | 一种基于进化多目标优化的选择性分层集成高斯过程回归软测量建模方法 |
CN110046434B (zh) * | 2019-04-18 | 2021-03-16 | 电子科技大学 | 基于高斯过程回归和遗传算法的机械零件可靠性设计方法 |
CN114861458B (zh) * | 2022-05-26 | 2023-07-14 | 西南交通大学 | 一种基于映射关系代理模型的桥上行车安全快速评价方法 |
CN116117818A (zh) * | 2023-02-28 | 2023-05-16 | 哈尔滨工业大学 | 一种机器人传感器手眼标定方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101944141A (zh) * | 2010-08-18 | 2011-01-12 | 北京理工大学 | 一种高效的基于模糊聚类自适应径向基全局优化方法 |
-
2013
- 2013-07-05 CN CN201310281958.8A patent/CN103336869B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103336869A (zh) | 2013-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103336869B (zh) | 一种基于高斯过程联立mimo模型的多目标优化方法 | |
Lukaczyk et al. | Active subspaces for shape optimization | |
Cheung et al. | Bayesian uncertainty analysis with applications to turbulence modeling | |
Huang et al. | Research on multi-fidelity aerodynamic optimization methods | |
Park et al. | Verification of unstructured grid adaptation components | |
Le Quilliec et al. | A manifold learning-based reduced order model for springback shape characterization and optimization in sheet metal forming | |
Seshadri et al. | Dimension reduction via gaussian ridge functions | |
Liu | Development of gradient-enhanced kriging approximations for multidisciplinary design optimization | |
Wang et al. | Development of a prediction model based on RBF neural network for sheet metal fixture locating layout design and optimization | |
Vicario et al. | Meta‐models in computer experiments: Kriging versus artificial neural networks | |
Kwon et al. | A trended Kriging model with R2 indicator and application to design optimization | |
Kahhal et al. | Multi-objective optimization of sheet metal forming die using genetic algorithm coupled with RSM and FEA | |
Boullé et al. | A mathematical guide to operator learning | |
Morse et al. | A multi-fidelity boundary element method for structural reliability analysis with higher-order sensitivities | |
Zheng et al. | Interval model updating using universal grey mathematics and Gaussian process regression model | |
Koziel et al. | Variable-fidelity CFD models and co-Kriging for expedited multi-objective aerodynamic design optimization | |
Yi et al. | Efficient global optimization using a multi-point and multi-objective infill sampling criteria | |
Toal et al. | An adjoint for likelihood maximization | |
Likeng et al. | Wing-body optimization based on multi-fidelity surrogate model | |
Lukaczyk et al. | Managing gradient inaccuracies while enhancing optimal shape design methods | |
Matsui et al. | Calibration of Spalart-Allmaras model for simulation of corner flow separation in linear compressor cascade | |
Xie | Gradient-based optimum aerodynamic design using adjoint methods | |
El-Beltagy et al. | Gaussian processes for model fusion | |
Liang et al. | Three‐stage non‐Gaussian homogeneous random field representation over manifolds | |
Kowalczyk et al. | Numpress–integrated computer system for analysis and optimization of industrial sheet metal forming processes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 |