4、工艺制度制订中的非线性规划.ppt

4、工艺制度制订中的非线性规划方法 ——若干实用规划方法及计算机算法 4.1 无约束最优化方法 4.1.1 一维收索方法 一维搜索方法实际上是用叠代法求一元函数的极值。叠代法的基本原理是：在给出广的极小点的一个初始点X(o)以后，计算一系列的X(k)(k=1,2,…)，希望点列{X(k)}的极限点X*是f(X)的一个极小点。在点列{X(k)}中相邻两个点之间有如下关系，即： 式中S(k)是一个向量，表示搜索方向；λk是一个实数，称为步长。一个n元函数f(x)在过点x(k)的方向为S(k)的直线的值可以表示为： 一维搜索法有：0.618法（黄金分割法）、分散法、切线法（又称牛顿法）和抛物线法等。 0.618法又称黄金分割法，它通过逐步缩小搜索区间求得一元函数极小点近似值的一种方法。 其基本思路是：在搜索区间[a,b]内适当插入两点x1,x2，它们将[a, b]分为三段，通过比较这两点的函数值而确定删去最左端或最右端，这算是叠代一次；然后在保留下来的区间上作同样的处置。如此叠代下去，可将搜索区间无限缩小。 如何确立两个插入点的位置呢？设区间[a, b]的长为1，在与a点相距分别为α和β的点处插入x1和x2，如图： (a) x1和x2两点在区间[a, b]中的位置是对称的，即 α+β=1 (b) 无论删掉哪一段，例如[x2,b]，在留下的区间[a,x1]里再插入一点x3，使得x3, x2在[a, x2]中的位置与x1, x2在[a, b]中的位置具有相同比例。这样保证了每次量代都以同一β的比率缩短区间 黄金分割法及计算框图 解题步骤： (a) 确定初始搜索区间[a, b]，给定精度要求。（δ0） (b) 计算x2=a+β(b-a)，及f2=f(x2) (c) 计算x1=b-β(b-a)及f1=f(x1) (d) 如果|x2-x1|δ，则输出极小点的近似值 x*=( x2+x1)/2以及f(x*)，叠代终止，否则转入(e) (e) 如果f1f2，则b=x2，x2=x1，f2=f1，转入(c)，否则转(f) (f) 令a=x1，x1=x2，f1=f2，x2= a+β(b-a) ，f2=f(x2)，转(d) 叠代原理是：在给出f(X)极小值的一个初始点 以后，计算一系列的 希望点列 的极限点 是函数f(X)的一个极小值。在点列 中相邻两点之间有如下关系，即 试中 是一个向量，表示搜索方向， 是一个实数，称为步长。一个n元函数f(X)在点 的方向 的直线可以表示为： * * 式中λ为独立变量，每一个λ值对应直线上的一个点，如图，沿直线求出函数的极小值，实际上就是求单变量λ的函数式（4-9）的极小点。 X（K）+λS（K） X（K） λ S（K） β 1 α α a b X1 X3 X2 代入α+β=1 另一负根含去 α=0.382 0.618法的计算步骤如下： β 1 α α a b X1 X3 X2 给立δ0，置 确立初始搜索区间[a, b] X2=a+β(b-a) f2=f(X1) X1=b-β(b-a) f1=f(X1) |X2-X1|δ f1f2 是 是 b=X1 X2=X1 f2=f1 不 a=X1 X1=X2 f1=f2 X2=a+β(b-a) f2=f(X2) 输出X1, f(X*)