实验二黄金分割法实验指导书.docVIP

实验二 黄金分割法实验指导书 一、实验目的： 1、加深对黄金分割法的基本理论和算法步骤的理解。 2、培养学生独立编制、调试机械优化算法程序的能力。 3、培养学生灵活运用优化设计方法解决工程实际问题的能力。 二、主要设备及软件配置 硬件：计算机（1台/人） 软件：VC6.0（c++） 三、算法程序框图及算法步骤 图1-2 黄金分割法程序框图 算法程序框图：如图1-2所示。 算法步骤： 1）给出初始搜索区间[a，b]及收敛精度ε，将λ赋以0.618。 2）按坐标点计算公式计算a1,a2；并计算其对应的函数值y1,y2。 3）根据区间消去法原理缩短搜索区间。为了能用原来的坐标点计算公式，需进行区间名称的代换，并在保留区间中计算一个新的试验点及其函数值。 如果y1≥y2，则新区间＝[a1 b]，a=a1,a1=a2,y1=y2, 插入新点a2=a+λ*(b-a) ,y2=f(a2); 如果y1＜y2，则新区间＝[a a2]，b=a2,a2=a1,y2=y1 插入新点a1=b-λ*(b-a) ,y1=f(a1); 4）检查区间是否缩短到足够小和函数值收敛到足够精度，如果收敛条件满足，则取最后两试验点的平均值作为极小点的数值近似解。如果条件不满足则转向步骤3）进行新的区间缩小。 四、实验内容与结果分析 1、根据算法程序框图和算法步骤编写计算机程序; 2、求解函数f(x)=x+20/x的极小点，初始搜索区间为[0.2，1], 收敛精度ε=0.01; 3、如果初始搜索区间为[0.1，10]，收敛精度ε=0.001,结果又如何？ 4、试分析初始搜索区间和收敛精度对搜索计算的影响。 附录：黄金分割法源程序 #include stdio.h #include conio.h #include math.h #include iostream.h #define e 0.001 #define tt 0.01 float function(float x) { // float y=xa*pow(x,3)-b*pow(x,2)-c*x+d; //求解的一维函数 float y=1*pow(x,4)-5*pow(x,3)-4*pow(x,2)-6*x+60; //求解的一维函数 return(y); } void finding(float a[3], float f[3]) { float t=tt,a1,f1,ia; //初始区间的下界值 a[0]=0; f[0]=function(a[0]); for (int i=0;;i++) { a[1]=a[0]+t; f[1]=function(a[1]); if ((f[1]f[0])) break; if ((fabs(f[1]-f[0])=e)) { t=-t; a[0]=a[1]; f[0]=f[1]; } else { if (ia==1) return; t=t/2; ia=1; } } for (i=0;;i++) { a[2]=a[1]+t; f[2]=function(a[2]); if(f[2]f[1]) break; t=2*t; a[0]=a[1]; f[0]=f[1]; a[1]=a[2]; f[1]=f[2]; } if (a[0]a[2]) { a1=a[0]; f1=f[0]; a[0]=a[2]; f[0]=f[2]; a[2]=a1; f[2]=f1; } return; } float gold(float *ff) { float a1[3],f1[3],a[4],f[4]; float aa; finding (a1,f1); a[0]=a1[0]; f[0]=f1[0]; a[3]=a1[2]; f[3]=f1[2]; a[1]=a[0]+0.382*(a[3]-a[0]); a[2]=a[0]+0.618*(a[3]-a[0]); f[1]=function(a[1]); f[2]=function(a[2]); for (int i=0;;i++) { if (f[1]=f[2]) { a[0]=a[1]; f[0]=