实验报告优化设计.docVIP

实验一 一维搜索方法的编程和调试 1、实验目的：（1）明确一维搜索方法的种类和数值寻优的方法 （2）加深对进退法和黄金分割法的基本理论和算法步骤的理解。 （3）培养学生独立编制、调试计算机程序的能力。 （4）掌握黄金分割优化方法程序的使用方法，培养学生解决工程实际问题的能力。 2、实验内容：用进退法确定实例函数的搜素区间，用黄金分割法确定函数的极优点。 3、算法原理及框图: 4、计算程序： 5、应用实例：(任选一个) （1）求函数f（x）=（x-3）^2 的最优解，h =0.1 ε=0.001 （2）求函数f（x）=8x^3-2x^2-7x+3 的最优解，h=0.001 ε=0.001 （3）求函数f（x）=x^2-7x+10 的最优解，h=0.1 ε=0.1 要求通过编制程序生成每步迭代的各个数据a,b,a1,f1,a2,f2,b-a,以及最终求优的结果。 6、计算结果。 实验二 无约束优化方法的编程和调试 1、实验目的：（1）明确无约束优化方法的方法和特点。 （2）加深对鲍威尔法的基本理论和算法步骤的理解。 （3）培养学生独立编制、调试计算机程序的能力。 （4）掌握用鲍威尔法优化方法程序的使用方法，培养学生解决工程实际问题的能力。 2、实验内容：编制鲍威尔法程序来确定实例函数的极优点。 3、算法原理及框图: 4、计算程序： 5、应用实例：(任选一个) （1）求函数f（x1 x2）=x1^2+2*x2^2-4*x1-2*x1*x2 的最优解 （2）求函数f（x1 x2）=(x1-2)^2+(x1-2*x2)^2的最优解 （3）求函数f（x1 x2）=1.5*x1^2+0.5*x2^2-2*x1-x1*x2的最优解 要求通过编制程序生成每步迭代的各个数据以及最终求优的结果。 6、计算结果。  #includeiostream.h double Fun(double x) //待求函数 { double a=(x-3)*(x-3); return a; } void main() { double Fun(double); double a0=0; //初始搜索点 double h=0.1; //搜索步长 double a1,a2,a3,f1,f2,f3; a1=a0; f1=Fun(a1); a2=a0+h; f2=Fun(a2); if(f2=f1) { h=-h; a3=a1;f3=f1; loop1: a1=a2;f1=f2; a2=a3;f2=f3; couta1=a1 a2=a2endl; coutf1=f1 f2=f2endl; loop2: a3=a2+h;f3=Fun(a3); if(f3=f2) { couta=min(a1 ,a2 )endl; coutb=max(a1 ,a2 )endl; coutf1=f1 f2=f2endl; } else goto loop1; } else goto loop2; //黄金分割法求值 double x,a=a1,b=a2,e=0.001; //精度 a1=b-0.618*(b-a);f1=Fun(a1); a2=a+0.618*(b-a);f2=Fun(a2); while((b-a)=e) { if(f1=f2) { b=a2; a2=a1; f2=f1; a1=b-0.618*(b-a); f1=Fun(a1); couta1=a1 a2=a2endl; coutf1=f1 f2=f2endl; } else { a=a1; a1=a2; f1=f2; a2=a+0.618*(b-a); f2=Fun(a2); couta1=a1 a2=a2endl; coutf1=f1 f2=f2endl; } } x=0.5*(b+a); coutx=xendl; }