现代优化方法综述.docx

1、现代优化方法综述1. 引言优化设计英文名是optimization design，从多种方案中选择最佳方案的设计方法。它以数学中的最优化理论为基础，以计算机为手段，根据设计所追求的性能目标，建立目标函数，在满足给定的各种约束条件下，寻求最优的设计方案。第二次世界大战期间，在军事上首先应用了优化技术。1967年，美国的R.L.福克斯等发表了第一篇机构最优化论文。1970年，C.S.贝特勒等用几何规划解决了液体动压轴承的优化设计问题后，优化设计在机械设计中得到应用和发展。随着数学理论和电子计算机技术的进一步发展，优化设计已逐步形成为一门新兴的独立的工程学科，并在生产实践中得到了广泛的应用。通常设计

2、方案可以用一组参数来表示，这些参数有些已经给定，有些没有给定，需要在设计中优选，称为设计变量。如何找到一组最合适的设计变量，在允许的围，能使所设计的产品结构最合理、性能最好、质量最高、成本最低（即技术经济指标最佳），有市场竞争能力，同时设计的时间又不要太长，这就是优化设计所要解决的问题。一般来说，优化设计有以下几个步骤：建立数学模型。选择最优化算法。程序设计。制定目标要求。计算机自动筛选最优设计方案等。2. 数学模型优化设计的数学模型是对优化设计工程问题的数学描述，它包含设计变量、目标函数和设计约束三个基本要素。2.1设计变量2.1.1基本参数a、定义：在设计过程中进行选择变化并最终确定的各项

3、独立参数称为设计变量。b、说明：在设计选择过程中，这些设计变量是变量，但它们一旦被确定后，设计对象也就完全确定了。最优化设计是研究怎样合理地优选这些设计变量的一种现代设计方法。在设计过程中，凡根据设计要求事先给定的，不是设计变量而是设计常量。2.1.2设计方案的表现形式a、设计空间：由n个设计变量为坐标所组成的时空间称作设计空间。b、设计变量的表示法（1）坐标表示法：一维问题一个设计变量数轴上的一个点 二维问题两个设计变量平面直角坐标系上的向量 三维问题三个设计变量空间直角坐标系的向量n维问题n个设计变量n维超越空间的向量一个“设计”方案，可用设计空间中的一点表示，此点可看成是设计变量向量的端

4、点（始点取在坐标原点），称作设计点。也即：在设计空间中的一个点，对应于一组设计变量的值，代表一个设计方案。设计空间包含了该项设计所有可能的设计方案。（2）向量表示法：二维问题二维向量 三维问题三维向量 n维问题n维向量2.1.3设计变量的选取a、维数：设计变量的数目称为最优化问题的维数。如有n个设计变量则称为n维问题。b、常选用的设计变量（1）结构的总体布置尺寸，如中心距。（2）元件的几何尺寸：长度，截面尺寸，某些点的坐标值。（3）材料的力学和物理特性：重量、惯性矩、力或力矩等。通常选择的设计变量都是构件的几个尺寸，因为这不仅可使问题相对简单些，而且由于很多实际结构的几个关系和材料特性已决定的

5、缘故。决定结构布置情况的设计变量的选取要复杂些。较困难的是选取表示材料特性的变量，因为通常所用材料的特性是离散值，选择这些变量时出现了设计变量不连续变化的这一特殊问题。c、设计变量的选择原则（1）对设计影响较大的参数选为设计变量（2）尽量减少设计变量的个数2.2设计约束2.2.1设计约束的种类a、定义：设计空间是所有设计方案的集合，但这些设计方案有些是工程上所不能接受的（例如面积取负等）。如果一个设计满足所有对它提出的要求，就称为可行（或可接受）设计。反之则称为不可行（或不可接受）设计。在设计过程中，为了得到可行的设计方案，必须根据实际要求，对设计变量的取值加以种种限制，这种限制条件称为约束条

6、件。即：一个可行设计必须满足的限制条件称为约束条件。b、分类法一 性能约束：针对性能要求而提出的限制条件称为性能约束。例如：强度条件、刚度或稳定性条件等等。边界约束：对设计变量的取值围加以限制的约束。例如允许选择的尺寸围。法二 等式约束：h（x）=0要求设计点在n维设计空间的约束曲面上不等式约束：g（x）0要求设计点在约束曲面一侧2.2.2可行域与非可行域设计可行域：凡满足所有约束条件的设计点，它在设计空间中的活动围称作可行域。2.3目标函数2.3.1目标函数的定义：a、定义在设计中，设计者总是希望所设计的产品或工程设施具有最好的使用性能，最小的质量或最紧凑的体积和最小的制造成本及最大的经济效

7、益。在最优设计中，可将所追求的设计目标（最优指标）用设计变量的函数形式表达出来。 目标函数是设计中预期要达到的目标，表达为各设计变量的函数表达式： 在优化设计中，用目标函数的大小来衡量设计方案的优劣，故目标函数又叫评价函数。 优化设计中，通常对目标函数求极小值。b、常用的目标函数（1）以成本最低构造目标函数。（2）按最小重量构造目标函数。（3）按几何要求：如最小体积，最小尺寸构造目标函数。（4）按机构的工作精度要求构造（5）按机构的运动轨迹最准确（6）满足应力要求（材料利用最好）（7）振动或噪声最小（齿轮振动，由侧隙产生，寻找一周期啮合点加速度平方根值最小）。（8）平均寿命最长（轴承的寿命计算

8、）。（9）冷却效果最好（轴承的热平衡计算）。（10）可靠性最高。2.4 优化设计数学模型的几何意义2.4.1 优化设计数学模型的一般形式a、模型形式选取设计变量，列出目标函数，给定约束条件后，便可构造最优化设计的数学模型。任何一个最优化问题均可归结为如下描述：在给定的约束条件下，选取适当的设计变量X，使其目标函数f（X）达到最优值，其数学表达式（数学模型）为： b、模型分类：（1）法一 有约束 无约束（2）法二 线性：目标函数和约束函数都是线性的。 非线性：目标函数和约束函数至少有一个为非线性2.4.2最优化问题的几何描述a、约束条件与可行域b、目标函数等值线（1）定义：目标函数是n为变量的函

9、数，它的图象只能在n+1维空间中描述出来。给定一组设计变量的值就相应有一个函数值（并相应在设计空间对应于一个设计点），具有相同函数值的点集在设计空间形成一个曲线或曲面，就是目标函数的等值面或等值线。（c、无约束最优解和约束最优解（1）无约束优化问题：在没有限制条件下，对设计变量求目标函数的极小点，即求等值面中心。（2）约束优化问题：在设计可行域寻求目标函数的极小点。2.4.3局部最优解和全局最优解一、单谷函数和多股函数只有一个极值点的函数称为单谷函数；具有两个以上局部极值点的函数称为多谷函数。二、局部最优解和全局最优解2.5优化设计数学模型大小的分类: n50大型 10n50中型n10小型 3

10、. 经典优化算法小结：3.1无约束优化方法 工程优化问题通常都是多维有约束优化问题，但需从一维无约束问题到多维无约束优化问题再到多维约束优化问题的由简单到复杂的循序渐进的研究过程。 无约束优化问题数学模型： 分类，从是否利用目标函数的导数信息，分直接法和间接法3.1.1 坐标轮换法3.1.1.1 坐标轮换法基本原理将多维无约束优化问题分解、转化为一系列一维优化问题，轮换沿各个坐标轴一维搜索，直到求得最优点。在每次迭代部，要依次沿各坐标轴进行N次（N为优化问题的维数）一维搜索。这种一维搜索是固定其它N-1维变量，视为常量，然后进行一维搜索，对于第k轮迭代，须重复N次该式的一维搜索，搜索的参数为a

11、jk（即要优化的参数是ajk），为相对第j维变量的搜索步长，搜索方向为第j维空间坐标的方向。当k轮迭代结束后，本轮搜索的重点作为下一轮的起点，即，然后投入下一轮迭代。3.1.1.2 该方法特点 不考虑目标函数本身的变化情况（函数特点），简单、效率低、收敛速度慢。3.1.2 共轭方向法3.1.2.1 共轭方向对于N维正定二次函数 (当N=2，为同心椭圆族)，H为函数f的黑塞矩阵(正定对称阵)。若存在两个方向向量，满足，则称与为共轭方向。如何构造共轭方向（二维）？对于某两点，沿方向（不平行）一维搜索得到两个最优点，构成方向，则可以证明与为共轭方向，即当然，这个结论可以从2维推广到N维。同样，说明对

12、于N维函数，有N个共轭方向。对于二次函数，只要经过N个一维搜索即可到达最优点（即N维空间完成一轮迭代）。对于大于二次的函数，则可能需要将上一轮迭代的终点作为新一轮迭代的起点。在构造迭代方程式时，可以用二次泰勒展开式来近似目标函数的等值面。3.1.2.2 共轭方向法基本原理第一轮迭代与坐标轮换法相同。将起点和N次一维搜索的末点组成一个新的方向，沿这个方向一维搜索，得到本轮迭代的终点。从第二轮起，舍去前一轮的第一个一维搜索方向，将前一轮的后N个一维搜索方向作为本轮迭代的前N个方向，这N个方向的一维搜索终点与本轮搜索的起点构成第N+1个一维搜索方向，沿这个方向做一维搜索，得到本轮搜索的终点。若不满足

13、精度要求，则重复迭代。3.1.2.3 共轭方向法的特点收敛速度比坐标轮换法有明显的提高，但前提是每次迭代所产生的新的方向与原来的N-1个方向之间要保持线性无关，若这些方向之间线性相关，则降低了搜索空间的维数，导致不能完全穷尽对设计空间每个方向的搜索，从而不能收敛于真正的最优解。3.1.3 梯度法（最速下降法）3.1.3.1 梯度法基本原理无约束优化的直接法（坐标轮换法和共轭方向法、鲍威尔法）没有考虑无约束优化最优解存在的必要条件（梯度为零），使用这一条件，可以设计出更为高效的算法，所谓间接法（梯度法、牛顿法、变尺度法）。梯度方向是函数值变化最快的方向，那么负梯度方向便是函数值下降最快的方向。从

14、这一点受启发，可以使迭代方向沿梯度方向进行一维搜索来再多维空间寻优。即搜索方向为梯度方向：，或，则迭代公式为。3.1.3.2 梯度法的特点前提是梯度存在。优点是算法简单。相邻两次迭代的搜索方向垂直。即证明：，即k轮迭代经过一次一维搜索由k点到达k+1点，那么，对于一维优化有，所以可见，相邻两轮迭代的搜索方向并不一致，为相互垂直的锯齿形过程。剃度法对于迭代出发点目标函数等值面偏心率为零时很有效，但对于有偏心的其效率就低了，随偏心率的增加，迭代终止的难度也在增加。可见这种搜索在接近目标时的收敛是比较慢（缺点）的，效率也就不会高了。剃度法一般并不作为工程中实际应用的方法，常用于其他方法的初始迭代（类

15、似于坐标轮换法）。3.2约束优化方法实际工程优化问题大多数为设计空间多维且带有约束条件的非线性优化问题。其数学模型为根据对约束条件处理方法的不同：直接法（约束坐标轮换法、随机方向法、复合形法、可行方向法）间接法（简约梯度法、惩罚函数法等）。直接法可以直接从可行域中找到最优解；将问题分解为一系列比较简单的子问题，用子问题的解逼近原问题的解。直接法简单直观、对目标函数要求不高；计算量大、收敛慢，因此效率低。3.2.1 约束随机方向搜索法（随机方向法）3.2.1.1 基本原理从可行域某一点出发，沿某一给定步长，并随机产生搜索方向，直到该方向同时满足可行性和下降性要求，沿着这个方向以该步长继续搜索，直到不满足可行性及下降性条件为止。把上述满足要求的终点作为新的起点，重新产生随机方向，如果能够找到一个合适的方向，同时满足条件，则沿该方向以原步长继续搜索；如若找不到适合的方向，则将步长减半，仍以该点为起点随机搜索，如果能找到新的方向，则沿该方向继续，如果不能，步长再减半。直到找不到新的搜索方向，且步长满足精度要求，则以该起点为最优点。一个需要说明的问题：从某一点出发，如何判断