优化设计程序.pdfVIP

《机械优化设计》 实验报告 班级 : 机械设计（ 5 ）班 姓名： 林斯宇 学号： 0901102028 1 目标函数方程式： F=ax^4+bx^3+cx^2+dx+e 其中 a=1.2,b=-1,c=-1,d=-10,c=181, 初始点 x0=0，步长 h=0.1 ，收 敛精度 e=0.0001 2 实验目的 ：通过上机编程，理解进退法、黄金分割法的原理，了 解计算机在优化设计中的应用。 3 黄金分割法的基本原理 一维收索是解函数极小值的方法之一， 其解法思想为沿某已知方向求目标函数的极小值点 。这种方法的适应面非常广。黄金分割法也是建立在区间削去法则原理基础上的试探方法。 黄金分割法是用于一元函数 f(x) 在给定初始区间 [a,b] 内搜索极小点 α* 的一种方法。它是优 化计算中的经典算法，以算法简单、 收敛速度均匀、效果较好而著称， 是许多优化算法的基 础，但它只适用于一维区间上的凸函数 [6] ，即只在单峰区间内才能进行一维寻优，其收敛效 率较低。其基本原理是：依照 “去劣存优 ”原则、对称原则、以及等比收缩原则来逐步缩小搜 索区间 [7] 。具体步骤是：在区间 [a,b] 内取点： a1 ，a2 把 [a,b] 分为三段。如果 f(a1)f(a2) ， 令 a=a1,a1=a2,a2=a+r*(b-a) ；如果 f(a1)f(a2) ，令 b=a2 ，a2=a1,a1=b-r*(b-a), 如果｜ (b-a)/b ｜ 和｜ (y1-y2)/y2 ｜都大于收敛精度 ε重新开始。因为 [a,b] 为单峰区间，这样每次可将搜索区 间缩小 0.618 倍或 0.382 倍，处理后的区间都将包含极小点的区间缩小，然后在保留下来的 区间上作同样的处理，如此迭代下去，将使搜索区 [a,b] 逐步缩小，直到满足预先给定的精度 时 ， 即 获 得 一 维 优 化 问 题 的 近 似 最优解。黄金分割法原理如图所示， 4 实验所编程序框图 （1 ）进退法计算框图 （2） 黄金分割法计算框图 5 程序源代码 6 程序 #includemath.h #includestdio.h #define f(x) 1.2*pow(x,4)+(-1)*pow(x,3)+(-1)*pow(x,2)+(-10)*pow(x,1 )+181 main() { double F1,F2,F3,W,x1,x2,x3,q,a,b,c; int n; n=1; printf(c=); scanf(%lf,c); printf(a=); scanf(%lf,a); printf(b=); scanf(%lf,b); printf(n a b x1 x2 F1 F2\n); x1=a+0.382*(b-a); x2=a+0.618*(b-a); F1=f(x1); F2=f(x2); printf(%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4 lf\n,n,a,b,x1,x2,F1,F2); n=n+1; do { if (F1=F2) { a=x1; x1=x2;