CN105608267A - 一种多变量全局优化算法 - Google Patents
一种多变量全局优化算法 Download PDFInfo
- Publication number
- CN105608267A CN105608267A CN201510958096.7A CN201510958096A CN105608267A CN 105608267 A CN105608267 A CN 105608267A CN 201510958096 A CN201510958096 A CN 201510958096A CN 105608267 A CN105608267 A CN 105608267A
- Authority
- CN
- China
- Prior art keywords
- meta
- sampling
- function
- points
- model
- 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 47
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 40
- 238000005070 sampling Methods 0.000 claims abstract description 72
- 238000000034 method Methods 0.000 claims abstract description 42
- 239000011159 matrix material Substances 0.000 claims abstract description 31
- 238000013461 design Methods 0.000 claims abstract description 26
- 239000013598 vector Substances 0.000 claims abstract description 19
- 230000004044 response Effects 0.000 claims description 22
- 238000004458 analytical method Methods 0.000 claims description 19
- 230000003044 adaptive effect Effects 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 238000010998 test method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 12
- 230000006870 function Effects 0.000 description 46
- 238000009826 distribution Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013433 optimization analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种多变量全局优化算法,包括以下步骤:初始采样、构建初始元模型集、在元模型集上搜索Pareto点集、更新适存度集合、判断新点集是否满足收敛条件、以增量RBF方法更新元模型集。本发明针对现有基于元模型的多变量优化方法采用的大都为单值元模型,一个变量函数须对应一个元模型,存在复杂度高、计算量大等问题,提出元模型集的概念,根据RBF函数的特点,将其系数向量转化为系数矩阵,使得多个变量可在同一元模型中进行仿真,大大优化了多变量优化算法。大大减少了计算工作量,降低了优化算法难度,节省宝贵计算资源,进而为工业设计中复杂机电产品结构优化设计提供一种新的思路与方法。
Description
【技术领域】
本发明涉及机电产品优化设计的技术领域,特别是多输入多输出问题优化设计的技术领域
【背景技术】
由于现代机电产品结构越来越复杂,对其进行优化仿真时需耗费很多计算机资源,这在实际中是无法接受的。因此,在过去近二十年中,响应面法应运而生并在工业中得到广泛应用,它是利用最小二乘回归法来实现目标模型的拟合,主要通过计算机实验产生的采样点来构造与目标模型相近似的数学模型进行仿真优化分析,因而能节约计算资源,降低计算量。
多输入多输出优化问题常常涉及至源模型的近似化问题,即仿真模型,计算模型、元模型等。这些模型促进了优化和概念搜索的发展,降低函数优化迭代次数,从而减少了计算机资源的消耗。使用这些仿真模型的关键之处在于近似模型的精度问题,如果精度不够准确,则求得的Pareto解集不能有效逼近Pareto边界。因此,选择合适的多输入多输出优化方法对机电产品***的优化仿真是非常重要和迫切的研究之一。
【发明内容】
本发明的目的在于引入元模型集的概念,利用RBF函数线性化的特点,将其系数向量转化为系数矩阵,从而使得多输入输出优化过程在同一元模型中完成,并提出一种与元模型集方法相匹配的增量Pareto适存度算法,从而降低优化算法难度,节省宝贵计算资源,进而为工业设计中复杂机电产品结构优化设计提供一种新的思路与方法。
为实现上述目的,本发明提出了多变量全局优化算法,包括以下步骤;
步骤一:初始采样:通过拉丁超立方采样方法采样构建初始样本点集X(i),主要依据设计变量个数p和采样点个数p2;
步骤二:构造初始元模型集:调用目标函数仿真,获得样本点集X(i)的响应值Y(i),利用X(i)、Y(i)来构建初始元模型集,适存度集合F(i)通过Pareto适存度函数公式来计算;
步骤三:在元模型集上搜索Pareto点集:在元模型集上搜索非支配点集Pa (k),以改进的增量拉丁超立方采样方法增加采样点,更新样本点集适存度值,在元模型集上搜索近似Pareto最优解,以全部边界点适存度值趋于零作为收敛准则;
步骤四:更新适存度集合:令X(i+1)=X(i)∪Pe (i),在目标函数进行仿真时,获得样本点集X(i+1)对应值Y(i+1),然后根据增量Pareto适存度算法来计算样本点X(i+1),更新适存度集合Fe (i+1);
步骤五:新点集是否满足收敛条件:判断经过精确分析的样本点的非支配点集Fe (i+1)是否满足收敛条件,如满足则退出循环,否则转向步骤三;
步骤六:以增量RBF方法更新元模型集,令i=i+1,再次循环。
优选的,在步骤一中,拉丁超立方采样的通用公式为:Xij=[πj(i)-Uij]/n,1≤i≤n,1≤j≤d,其中,πj(1),...,πj(n)是一个随机排列,区间为[1,n],Uij,…,U[0,1]是区间[0,1]的一个均匀随机分布,d为采样点的维度,n为采样点个数,具体采样步骤为:
1)设采样水平为m个点,维度空间为n,将设计空间每个维度进行m等分,每个维度分成m个区间:(0,1/m),(1/m,2/m),…,(1-1/m,1);
2)在每个维度m个区间内随机取值,得值且 以次类推;
3)根据每个维度值随机选取配对,形成m个采样点的n维数据。
优选的,在步骤二中,初始元模型集的具体构造步骤为:
1)选择合适的径向基函数φ,一般常用Cubic基函数,其表达式为φ(r)=(r+c)3,0<c<1;
2)采样点的生成:主要通过计算机试验方法在设计空间里生成,然后对其进行函数估算,从而得到对应响应值;
3)将m个采样点xi及对应响应值f代入公式Aλ=f*得出线性方程组,m×m阶矩阵A中元素Aij=φ(r),式中,r—采样点xi和xj的距离范数;
4)对线性方程组求解,得出扩展RBF函数中基函数对应权重系数λ=[λ1,λ2,...,λn]T,再次将采样点xi代入对应方程组即可得出RBF响应面确切表达式;
5)同时进行多点估值,则 相应的,RBF元模型的响应值f*(x)相应变为矩阵,这样,RBF元模型就被转变为同一元模型下的多值模型,即元模型集合。
优选的,在步骤3)中,公式Aλ=f*中,f*=[f*(x1),f*(x2),...,f*(xn)]T,
Ai,j=φ(||xi-xj||)(i,j=1,2...,n),λ=[λ1,λ2,...,λn]T,由于未知参数个数多于由采样点构造的方程个数,则Cubic基函数为欠定的,所以,系数λ还必须满足式因此可得到n+p个方程,其矩阵方程式为 式中,Gi,j=g(xi)(i=1,2...,n;j=1,2,...,p),c=[c1,c2,...,cp]T,因此,即可求解参数λ和c。
优选的,在步骤二中,适存度集合为 式中,A={xi,i=1,2,…,m}为设计点集合,Fi—第i个给定点xi在A中的适存度表达式,—第i个给定点映射的第k个子目标函数比例缩放值;k=i,2,…,m;l—边界指数,取值为1;max—计算给定点xi对应的目标函数向量和其它给定点xj(j=1,2,…,m,j≠i)对应目标函数向量之差集合的最大值。
优选的,在步骤二中,为了计算Pareto适存度集合,设函数fs *为第i个给定点函数值与其它点j函数值之差的最小值,那么: 式中,j=1,…n,且i≠j;k—子目标函数,k=1,…m,则上式的行向量矩阵为:
对上式取最大值即为Pareto适存度函数公式值。因 设 则上式变换为:
优选的,在步骤三中,改进的增量拉丁超立方采样方法的具体算法为:
1)建造(M+N)×K分析矩阵,将矩阵中所有元素清0;
2)计算原来N个采样点在分析矩阵中的位置,将相对应元素置1;
3)查找矩阵中元素为0的点,将其收集至空余坐标集中;
4)用标准拉丁超立方采样方法生成(M+N)个采样点,按空余坐标集中的位置生成新增采样点坐标。
优选的,在步骤三中,收敛标准为:其中,ε1—大小1的正数,
优选的,在步骤四中,增量Pareto适存度算法具体步骤为:
1)设向量init_def为上一次Pareto适存度值,fit为当前给定点对应的函数值矩阵,并且将fit按比例缩放至[1,0];
2)构造向量def为当前Pareto适存度值域,如果init_def非空,将其加入def前面部分;
3)设i为迭代过程中当前给定点对应函数值cur_fit,将对所有给定点进行循环,设j迭代过程中与i不同的给定点对应函数值oth_fit,对前i-1个给定点进行循环;
4)设f*=cur_fit-oth_fit,那么 。如def第i个元素小于则将第i个元素替换为如def第i个元素小于则将第i个元素替换为
5)令i++,j++,再次循环。
优选的,在步骤五中,经过精确分析的样本点的非支配点集是否满足的收敛条件为::其中,ε2—大小1的正数。
有益效果:本发明针对现有基于元模型的多变量优化方法采用的大都为单值元模型,一个变量函数须对应一个元模型,存在复杂度高、计算量大等问题,提出元模型集的概念,根据RBF函数的特点,将其系数向量转化为系数矩阵,使得多个变量可在同一元模型中进行仿真,大大优化了多变量优化算法。然后,针对通过采样点计算Pareto适存度矩阵困难等问题,提出了增量迭代式Pareto适存度计算方法,利用上一次迭代产生的适存度值直接更新,大大减少了计算工作量,降低了优化算法难度,节省宝贵计算资源,进而为工业设计中复杂机电产品结构优化设计提供一种新的思路与方法。
【附图说明】
图1是本发明的元模型集上搜索近似Pareto流程图;
图2a是本发明中的Pareto边界分布图;
图2b是本发明中的Pareto精确分析点分布图;
图2c是本发明中的Pareto解集分布图;
图2d是本发明中的Pareto收敛曲线图;
图3是本发明中的五杆平面桁架结构图;
图4a是本发明中的Pareto边界分布图二;
图4b是本发明中的Pareto精确分析点分布图二;
图4c是本发明中的Pareto解集分布图二;
图4d是本发明中的Pareto收敛曲线图二;
【具体实施方式】
以下各实施例的说明是参考附加的图式,用以例示本发明可用以实施的特定实施例。本发明所提到的方向用语,例如「上」、「下」、「前」、「后」、「左」、「右」、「内」、「外」、「侧面」等,仅是参考附加图式的方向。因此,使用的方向用语是用以说明及理解本发明,而非用以限制本发明。在图中,结构相似的单元是以相同标号表示。
下面结合附图和具体实施例,对本发明进行详细说明。
如图1至图4d所示,本发明的多变量全局优化算法,包括以下步骤:
步骤一:初始采样:通过拉丁超立方采样方法采样构建初始样本点集X(i),主要依据设计变量个数p和采样点个数p2;
步骤二:构造初始元模型集:调用目标函数仿真,获得样本点集X(i)的响应值Y(i),利用X(i)、Y(i)来构建初始元模型集,适存度集合F(i)通过Pareto适存度函数公式来计算;
步骤三:在元模型集上搜索Pareto点集:在元模型集上搜索非支配点集Pa (k),以改进的增量拉丁超立方采样方法增加采样点,更新样本点集适存度值,在元模型集上搜索近似Pareto最优解,以全部边界点适存度值趋于零作为收敛准则;
步骤四:更新适存度集合:令X(i+1)=X(i)∪Pe (i),在目标函数进行仿真时,获得样本点集X(i+1)对应值Y(i+1),然后根据增量Pareto适存度算法来计算样本点X(i+1),更新适存度集合Fe (i+1);
步骤五:新点集是否满足收敛条件:判断经过精确分析的样本点的非支配点集Fe (i+1)是否满足收敛条件,如满足则退出循环,否则转向步骤三;
步骤六:以增量RBF方法更新元模型集,令i=i+1,再次循环。
进一步的,在步骤一中,拉丁超立方采样的通用公式为:Xij=[πj(i)-Uij]/n,1≤i≤n,1≤j≤d,其中,πj(1),...,πj(n)是一个随机排列,区间为[1,n],Uij,…,U[0,1]是区间[0,1]的一个均匀随机分布,d为采样点的维度,n为采样点个数,具体采样步骤为:
1)设采样水平为m个点,维度空间为n,将设计空间每个维度进行m等分,每个维度分成m个区间:(0,1/m),(1/m,2/m),…,(1-1/m,1);
2)在每个维度m个区间内随机取值,得值且 以次类推;
3)根据每个维度值随机选取配对,形成m个采样点的n维数据。
进一步的,在步骤二中,初始元模型集的具体构造步骤为:
1)选择合适的径向基函数φ,一般常用Cubic基函数,其表达式为φ(r)=(r+c)3,0<c<1;
2)采样点的生成:主要通过计算机试验方法在设计空间里生成,然后对其进行函数估算,从而得到对应响应值;
3)将m个采样点xi及对应响应值f代入公式Aλ=f*得出线性方程组,m×m阶矩阵A中元素Aij=φ(r),式中,r—采样点xi和xj的距离范数;
4)对线性方程组求解,得出扩展RBF函数中基函数对应权重系数λ=[λ1,λ2,...,λn]T,再次将采样点xi代入对应方程组即可得出RBF响应面确切表达式;
5)同时进行多点估值,则 相应的,RBF元模型的响应值f*(x)相应变为矩阵。
进一步的,在步骤二的3)中,公式Aλ=f*中,f*=[f*(x1),f*(x2),...,f*(xn)]T,Ai,j=φ(||xi-xj||)(i,j=1,2...,n),λ=[λ1,λ2,...,λn]T,由于未知参数个数多于由采样点构造的方程个数,则Cubic基函数为欠定的,所以,系数λ还必须满足式因此可得到n+p个方程,其矩阵方程式为 式中,Gi,j=g(xi)(i=1,2...,n;j=1,2,...,p),c=[c1,c2,...,cp]T,因此,即可求解参数λ和c。
进一步的,在步骤二中,适存度集合为
进一步的,在步骤二中,为了计算Pareto适存度集合,设函数fs *为第i个给定点函数值与其它点j函数值之差的最小值,那么: 式中,j=1,…n,且i≠j;k—子目标函数,k=1,…m,则上式的行向量矩阵为:
对上式取最大值即为Pareto适存度函数公式值。因 设 则上式变换为:
进一步的,在步骤三中,改进的增量拉丁超立方采样方法的具体算法为:
1)建造(M+N)×K分析矩阵,将矩阵中所有元素清0;
2)计算原来N个采样点在分析矩阵中的位置,将相对应元素置1;
3)查找矩阵中元素为0的点,将其收集至空余坐标集中;
4)用标准拉丁超立方采样方法生成(M+N)个采样点,按空余坐标集中的位置生成新增采样点坐标。
进一步的,在步骤三中,收敛标准为:其中,ε1—大小1的正数,具体搜索流程如图1所示。
进一步的,在步骤四中,增量Pareto适存度算法具体步骤为:
1)设向量init_def为上一次Pareto适存度值,fit为当前给定点对应的函数值矩阵,并且将fit按比例缩放至[1,0];
2)构造向量def为当前Pareto适存度值域,如果init_def非空,将其加入def前面部分;
3)设i为迭代过程中当前给定点对应函数值cur_fit,将对所有给定点进行循环,设j迭代过程中与i不同的给定点对应函数值oth_fit,对前i-1个给定点进行循环;
4)设f*=cur_fit-oth_fit,那么 如def第i个元素小于则将第i个元素替换为如def第i个元素小于则将第i个元素替换为
5)令i++,j++,再次循环。
进一步的,在步骤五中,经过精确分析的样本点的非支配点集Fe (i+1)是否满足的收敛条件为::其中,ε2—大小1的正数。
为验证本发明提出的算法,采用凸边界优化问题和五杆平面桁架结构设计优化两个例子进行验证。
(1)凸边界优化问题:
测试函数:
收敛标准值ε1=1.02,ε2=1.008,采用本文所述算法,经12次迭代后满足收敛条件,结果如图2a、b、c、d所示,统计结果如表1所示。
表1例1结果分析表
从图2(a)可知,采用上述方法得到的Pareto边界是均匀分布的,获得的近似Pareto边界总体上趋向于实际的Pareto边界,从而使得由近似模型获得的非支配点能直接收敛于实际Pareto边界。由图2(b)可知,逐次增量分析的大部分设计点分布于Pareto边界上,且比例高达61%,而且分析次数也大大降低。
五杆平面桁架结构设计优:
主要特性参数:杨氏模量E=7030.8×107N/m2,密度ρ=27680.37N/m3,许用应力σa=±1.7577×108N,5个设计变量为每根杆的截面区域,其定义范围为:
{64.516mm2<Ai<6451.6mm2}i=1,2,...,5
在满足五杆平面桁架结构屈服强度和稳定性条件的前提下,优化目标为最小化桁架的承载值。经分析可知,其函数为:
{-1.7577×108N<σi<1.7577×108N}i=1,2,...,5
设收敛标准值ε1=1.02,ε2=1.008,经12次迭代后得到94个Pareto边界点,占精确分析点总量的48%,且精确分析次数为249次。例二Pareto边界分布、精确分析点分布、解集分布及收敛曲线如图4a、b、c、d所示,计算结果如表2所示。
表2例2结果分析表
结合两个例子的测试结果可以得出,基于元模型集的多变量全局优化算法相对于多目标遗传算法来讲,获取Pareto边界迭代次数少,进行精确分析的次数也大大降低,这对于节约计算机资源,提高仿真优化效率具有重要意义。
本发明公开了一种多变量全局优化算法,包括以下步骤:初始采样、构建初始元模型集、在元模型集上搜索Pareto点集、更新适存度集合、判断新点集是否满足收敛条件、以增量RBF方法更新元模型集。本发明针对现有基于元模型的多变量优化方法采用的大都为单值元模型,一个变量函数须对应一个元模型,存在复杂度高、计算量大等问题,提出元模型集的概念,根据RBF函数的特点,将其系数向量转化为系数矩阵,使得多个变量可在同一元模型中进行仿真,大大优化了多变量优化算法。然后,针对通过采样点计算Pareto适存度矩阵困难等问题,提出了增量迭代式Pareto适存度计算方法,利用上一次迭代产生的适存度值直接更新,大大减少了计算工作量,降低了优化算法难度,节省宝贵计算资源,进而为工业设计中复杂机电产品结构优化设计提供一种新的思路与方法。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内,本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (10)
1.一种多变量全局优化算法,其特征在于:包括以下步骤:
步骤一:初始采样:通过拉丁超立方采样方法采样构建初始样本点集X(i),主要依据设计变量个数p和采样点个数p2;
步骤二:构造初始元模型集:调用目标函数仿真,获得样本点集X(i)的响应值Y(i),利用X(i)、Y(i)来构建初始元模型集,适存度集合F(i)通过Pareto适存度函数公式来计算;
步骤三:在元模型集上搜索Pareto点集:在元模型集上搜索非支配点集以改进的增量拉丁超立方采样方法增加采样点,更新样本点集适存度值,在元模型集上搜索近似Pareto最优解,以全部边界点适存度值趋于零作为收敛准则;
步骤四:更新适存度集合:令在目标函数进行仿真时,获得样本点集X(i+1)对应值Y(i+1),然后根据增量Pareto适存度算法来计算样本点X(i+1),更新适存度集合
步骤五:新点集是否满足收敛条件:判断经过精确分析的样本点的非支配点集是否满足收敛条件,如满足则退出循环,否则转向步骤三;
步骤六:以增量RBF方法更新元模型集,令i=i+1,再次循环。
2.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤一中,拉丁超立方采样的通用公式为:Xij=[πj(i)-Uij]/n,1≤i≤n,1≤j≤d,其中,πj(1),...,πj(n)是一个随机排列,区间为[1,n],Uij,…,U[0,1]是区间[0,1]的一个均匀随机分布,d为采样点的维度,n为采样点个数,具体采样步骤为:
1)设采样水平为m个点,维度空间为n,将设计空间每个维度进行m等分,每个维度分成m个区间:(0,1/m),(1/m,2/m),…,(1-1/m,1);
2)在每个维度m个区间内随机取值,得值且 以次类推;
3)根据每个维度值随机选取配对,形成m个采样点的n维数据。
3.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤二中,初始元模型集的具体构造步骤为:
1)选择合适的径向基函数φ,一般常用Cubic基函数,其表达式为φ(r)=(r+c)3,0<c<1;
2)采样点的生成:主要通过计算机试验方法在设计空间里生成,然后对其进行函数估算,从而得到对应响应值;
3)将m个采样点xi及对应响应值f代入公式Aλ=f*得出线性方程组,m×m阶矩阵A中元素Aij=φ(r),式中,r—采样点xi和xj的距离范数;
4)对线性方程组求解,得出扩展RBF函数中基函数对应权重系数λ=[λ1,λ2,...,λn]T,再次将采样点xi代入对应方程组即可得出RBF响应面确切表达式;
5)同时进行多点估值,则相应的,RBF元模型的响应值f*(x)相应变为矩阵,这样,RBF元模型就被转变为同一元模型下的多值模型,即元模型集合。
4.根据权利要求3所述的多变量全局优化算法,其特征在于,在步骤3)中,公式Aλ=f*中,f*=[f*(x1),f*(x2),...,f*(xn)]T,
Ai,j=φ(||xi-xj||)(i,j=1,2...,n),λ=[λ1,λ2,...,λn]T,由于未知参数个数多于由采样点构造的方程个数,则Cubic基函数为欠定的,所以,系数λ还必须满足式因此可得到n+p个方程,其矩阵方程式为式中,Gi,j=g(xi)(i=1,2...,n;j=1,2,...,p),c=[c1,c2,...,cp]T,因此,即可求解参数λ和c。
5.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤二中,适存度集合为式中,A={xi,i=1,2,...,m}为设计点集合,Fi—第i个给定点xi在A中的适存度表达式,—第i个给定点映射的第k个子目标函数比例缩放值;k=i,2,...,m;l—边界指数,取值为1;max—计算给定点xi对应的目标函数向量和其它给定点xj(j=1,2,...,m,j≠i)对应目标函数向量之差集合的最大值。
6.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤二中,为了计算Pareto适存度集合,设函数为第i个给定点函数值与其它点j函数值之差的最小值,那么:式中,j=1,...n,且i≠j;k—子目标函数,k=1,...m,则上式的行向量矩阵为:
对上式取最大值即为Pareto适存度函数公式值。因设则上式变换为:
7.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤三中,改进的增量拉丁超立方采样方法的具体算法为:
1)建造(M+N)×K分析矩阵,将矩阵中所有元素清0;
2)计算原来N个采样点在分析矩阵中的位置,将相对应元素置1;
3)查找矩阵中元素为0的点,将其收集至空余坐标集中;
4)用标准拉丁超立方采样方法生成(M+N)个采样点,按空余坐标集中的位置生成新增采样点坐标。
8.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤三中,收敛标准为:其中,ε1—大小1的正数。
9.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤四中,增量Pareto适存度算法具体步骤为:
1)设向量init_def为上一次Pareto适存度值,fit为当前给定点对应的函数值矩阵,并且将fit按比例缩放至[1,0];
2)构造向量def为当前Pareto适存度值域,如果init_def非空,将其加入def前面部分;
3)设i为迭代过程中当前给定点对应函数值cur_fit,将对所有给定点进行循环,设j迭代过程中与i不同的给定点对应函数值oth_fit,对前i-1个给定点进行循环;
4)设f*=cur_fit-oth_fit,那么 如def第i个元素小于则将第i个元素替换为如def第i个元素小于则将第i个元素替换为
5)令i++,j++,再次循环。
10.根据权利要求1所述的多变量全局优化算法,其特征在于,在步骤五中,经过精确分析的样本点的非支配点集是否满足的收敛条件为::其中,ε2—大小1的正数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510958096.7A CN105608267A (zh) | 2015-12-21 | 2015-12-21 | 一种多变量全局优化算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510958096.7A CN105608267A (zh) | 2015-12-21 | 2015-12-21 | 一种多变量全局优化算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105608267A true CN105608267A (zh) | 2016-05-25 |
Family
ID=55988199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510958096.7A Pending CN105608267A (zh) | 2015-12-21 | 2015-12-21 | 一种多变量全局优化算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105608267A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106406087A (zh) * | 2016-07-08 | 2017-02-15 | 上海大学 | 基于变换函数和填充函数的智能控制***优化算法 |
CN109002724A (zh) * | 2018-06-07 | 2018-12-14 | 南京师范大学 | 一种基于紧支撑径向基函数的dem局部脱密与恢复方法 |
CN111985092A (zh) * | 2020-07-30 | 2020-11-24 | 哈尔滨工业大学 | 一种智能汽车仿真测试矩阵生成方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102024082A (zh) * | 2010-12-15 | 2011-04-20 | 同济大学 | 实现汽车仪表板结构***多学科多目标优化的方法 |
US20150019173A1 (en) * | 2013-07-09 | 2015-01-15 | International Business Machines Corporation | Multiobjective optimization through user interactive navigation in a design space |
-
2015
- 2015-12-21 CN CN201510958096.7A patent/CN105608267A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102024082A (zh) * | 2010-12-15 | 2011-04-20 | 同济大学 | 实现汽车仪表板结构***多学科多目标优化的方法 |
US20150019173A1 (en) * | 2013-07-09 | 2015-01-15 | International Business Machines Corporation | Multiobjective optimization through user interactive navigation in a design space |
Non-Patent Citations (1)
Title |
---|
魏昕: "基于元模型的全局优化算法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106406087A (zh) * | 2016-07-08 | 2017-02-15 | 上海大学 | 基于变换函数和填充函数的智能控制***优化算法 |
CN109002724A (zh) * | 2018-06-07 | 2018-12-14 | 南京师范大学 | 一种基于紧支撑径向基函数的dem局部脱密与恢复方法 |
CN109002724B (zh) * | 2018-06-07 | 2021-03-23 | 南京师范大学 | 一种基于紧支撑径向基函数的dem局部脱密与恢复方法 |
CN111985092A (zh) * | 2020-07-30 | 2020-11-24 | 哈尔滨工业大学 | 一种智能汽车仿真测试矩阵生成方法 |
CN111985092B (zh) * | 2020-07-30 | 2024-05-31 | 哈尔滨工业大学 | 一种智能汽车仿真测试矩阵生成方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dey et al. | Metamodel based high-fidelity stochastic analysis of composite laminates: A concise review with critical comparative assessment | |
CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
CN104866692B (zh) | 一种基于自适应代理模型的飞行器多目标优化方法 | |
Iuliano | Global optimization of benchmark aerodynamic cases using physics-based surrogate models | |
CN106372278A (zh) | 一种联合考虑输入参数不确定性和代理模型不确定性的灵敏度分析方法 | |
Kubalík et al. | Symbolic regression driven by training data and prior knowledge | |
CN110262245A (zh) | 基于事件触发机制的多智能体***的控制器设计方法 | |
Zhang et al. | A new extension of newton algorithm for nonlinear system modelling using RBF neural networks | |
Li et al. | Reanalysis assisted metaheuristic optimization for free vibration problems of composite laminates | |
CN105608267A (zh) | 一种多变量全局优化算法 | |
Fan et al. | A novel hybrid algorithm for solving multiobjective optimization problems with engineering applications | |
CN114996995A (zh) | 一种超材料隔振单元性能预报方法和*** | |
CN107846022B (zh) | 基于ilutp预处理并行迭代法的大规模配电网潮流分析方法 | |
Petrone et al. | A probabilistic non-dominated sorting GA for optimization under uncertainty | |
CN103065015A (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
Hofwing et al. | D-optimality of non-regular design spaces by using a Bayesian modification and a hybrid method | |
CN112749517A (zh) | 一种公差优化分配方法 | |
Jo et al. | Variable-fidelity aerodynamic design using gradient-enhanced kriging surrogate model with regression | |
CN106055742B (zh) | 基于设计熵的产品自适应设计更改方法 | |
Cheng et al. | A hybrid artificial neural network method with uniform design for structural optimization | |
CN110727569B (zh) | 一种多因素综合环境下的***性能评估方法 | |
CN112488248A (zh) | 一种基于卷积神经网络的代理模型的构建方法 | |
CN107797966A (zh) | 交互式动态影响图的最优k模型求解方法 | |
Sun et al. | Optimization via simulation using Gaussian process-based search | |
Li et al. | Comparisons of three meta-models for structural reliability analysis: RSM, ANN and SVR |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160525 |