工程优化设计-无约束直接法.pptVIP

工程优化设计 黄正东 二0一0年九月 内容提要 工程优化问题建模 优化数学理论 一维搜索方法 无约束问题直接搜索方法 无约束问题间接接搜索方法 约束问题直接搜索方法 线性规划与二次规划问题求解 约束问题间接搜索方法 启发式算法 优化软件系统 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 无约束直接搜索方法 Powell法另一形式(保证线性无关) 引理: 设?k =det(S1/|S1|, S2/|S2|,…, Sn/|Sn|), ?k+1=det(S1/|S1|, S2/|S2|,…, Sm-1/|Sm-1|, Sn+1/|Sn+1|, Sm+1/|Sm+1|,…,Sn/|Sn|), 则 ?k+1= ?k *am/|Xn(k)-X0(k)|, 这里 Xm(k)=Xm-1(k)+amSm 证明: ek,i=Si(k)/| Si(k) |, ek+1,i= ek,i,(i!=m), ek+1,m= ek,n+1 ek,n+1=(Xn(k)-X0(k))/rk=(?aiek,i)/rk, rk=|Xn(k)-X0(k)| ?k+1=det(ek,1, ek,2,…, ek,m-1, ek,n+1, ek,m+1,…, ek,n)= ?k *am/rk * 一.坐标轮换法 直接法: 利用迭代过程已有信息和再生信息进行试探和求优, 不需要用到函数导数和分析性质。 (1) 算法思想 轮换地以每一个坐标轴作为一维 搜索方向的优化搜索方法。 一.坐标轮换法 (2) 算法 1. 初始化， ?0，M，n, k=1, x=x(0). 2. 对于i=1,2,…,n， 进行 2.1 Si=ei; 2.2 ?i=min f(x+?Si) 2.3 x= x+?iSi, f=f(x). 3. 如果|?iSi|?, 或者 kM, 转步4; 否则，k=k+1,转步2。 4. 输出x, f. 结束。 一.坐标轮换法 (3) 算法分析 对于维数较高的优化问题，搜索时间过长，一般当n10时， 则不应采用此方法。 2. 算法效率与f(x)形态有关。 二.单纯形方法 (1) 算法思想 单纯形（Simplex）概念: n维空间中的n+1面体 f(xh)=max{f(x0), f(x1), …,f(xn)} ---高点 f(xl)=min{f(x0), f(x1), …,f(xn)} ---低点 x= ---除最高点外的所有点的形心 xh x2 x1 x xr xr=x+a(x-xh) ---反射点，a反射系数，一般取 a=1. 二.单纯形方法 分三种情况选择新点: xe=x+y(xr-x) 若f(xr)f(xe),则xe - xh 否则， xr - xh 1。如果f(xl)f(xr), 进一步扩展 扩展系数y1, 一般y=2. xh x2 xl x xr xe f(xr)比所有单纯形点上值小 不进一步扩展,避免狭窄单纯形 二.单纯形方法 分三种情况选择新点: xe=x+y(xr-x) 若f(xr)f(xe),则xe - xh 否则， xr - xh 1。如果f(xl)f(xr), 进一步扩展 扩展系数y1, 一般y=2. xh x2 xl x xr xe f(xr)比所有单纯形点上值小 二.单纯形方法 分三种情况选择新点: 2。如果max{f(xi),i?h}?f(xr)?f(xl),则 xr - xh xh x2 x1 x xr xe f(xr)在单纯形点上值之间, 但最少比第二最大值小. 分三种情况选择新点: 3。如果f(xr)max{f(xi),i?h},则 f(xh’)=min{f(xh),f(xr)} 压缩xc=x+b(xh’-x), 0b1, 压缩系数，一般b=0.5. 如果f(xc)?f(xh’),xc -xh 否则，以xl为中心压缩整个单纯形： xi=xi+0.5(xl-xi), i=0,1,2,…,n xh x2 xl x2 xh xh x2 x1 x xr xc xh x2 x xr xc x1 f(xr)