黄金分割法程序.doc

黄金分割法黄金分割法是通过不断缩短搜索区间的长度来寻求一维函数的极小点，这种方法的基本原理是：在搜索区间a,b内按如下规则对称地取两点a1和a2a1=a+0.382(b-a); a2=a+0.618(b-a);黄金分割法的搜索过程是：1） 给出初始搜索区间 a,b 及收敛精度 e ，将 赋以0.6182） 计算a1 和 a2，并计算起对应的函数值 f(a1),f(a2); ， 3） 根据期间消去法原理缩短搜索区间，为了能用原来的坐标点计算公式，需进行区间名城的代换，并在保留区间中计算一个新的试验点及其函数值。4） 检查区间是否缩短到足够小和函数值收敛到足够近，如果条件不满足则返回到步骤2。5） 如果条件满足，则取最后两试验点的平均值作为极小点的数值近似解。黄金分割法的流程图及程序清单需要说明的是搜索区间a,b不需要给定，只需输入搜索精度e；程序由四个子程序构成；（1）：输入输出子程序io（）；（2）：float fc (float x)求输入函数在某一点的值；（3）void findqujian(float a3,float f3)确定搜索区间；（4）：float xunyou(float *value)寻找最小值 黄金分割法程序运行截图 #include iostream.h#include math.h#include stdio.h#include conio.h#define steplength 0.01#define n 5float e;float a,b,c,d,g;float q5;void io()cout假设多项式的最高次幂为4endlendl;cout设多项式的一般形式为f=a*x4+b*x3+c*x2+d*x+eendlendl;cout请输入要求解的目标多项式的系数endlendl;printf(a=);scanf(%f,&a);q4=a; printf(b=); scanf(%f,&b);q3=b;printf(c=); scanf(%f,&c);q2=c; printf(d=); scanf(%f,&d);q1=d; printf(g=);scanf(%f,&g);q0=g; coutendl; cout请输入搜索精度eendlendl; scanf(%f,&e); cout=0;i-) u=u*x+qi;return u; void findqujian(float a3,float f3)float t=float(steplength), a1,f1,ia;a0=0;f0=fc(a0);for(int i=0;i+)a1=a0+t;f1=fc(a1);if(f1=e)t=-t;a0=a1;f0=f1;else if(ia=1)return;t=t/2;ia=1;for(i=0;i+)a2=a1+t;f2=fc(a2);if(f2f1) break;t=2*t;a0=a1; f0=f1;a1=a2; f1=f2;if(a0a2)a1=a0;f1=f0;a0=a2;f0=f2;a2=a1; f2=f1;return;float xunyou(float *value)float a13,f13,a4,f4;float aa;findqujian(a1,f1);a0=a10;f0=f10;a3=a12;f3=f12;a1=a0+float(0.382)*(a3-a0);a2=a0+float(0.618)*(a3-a0);f1=fc(a1);f2=fc(a2);for(int i=0;i+)if(f1=f2)a0=a1;f0=f1;a1=a2;f1=f2;a2=a0+float(0.618)*(a3-a0);f2=fc(a2);elsea3=a2;f3=f2;a2=a1;f2=f1;a1=a0+float(0.382)*(a3-a0);f1=fc(a1);if(fabs(a3-a0)e)aa=(a1+a2)/2; *value=fc(aa);break;return(aa);void main()float xx,value;io(); xx=xunyou(&value);printf(f(x)=%2.1f*x4+%2.1f*x3+%2.1f*x2+%2.1f*x+%2.1f,a,b,c,d,g); cout取得最小值的坐标为endlendl;printf(nzuobiao*=%fnn,xx);cout函数的最小值是endl;printf(nminmum*=%fn,value);getch();坐标轮换发#include #include #define EP 0.001main()float M1,M2,min,H,x1,x2;int n=0;x1=100;x2=25;do n=n+1; M1=x1;x1=2+x2; /* 直接用数学方法求，最小值 ：X=-B/2/A；Y=X1*X1-2*(2+X2)*X1+2*X2*X2; */ M2=x2;x2=x1/2; /* 直接用数学方法求，最小值 ：X=-B/2/A；Y=2(X2*X2-X2*X1)-X1*X1-4*X1; */ H=x2+x1-M1-M2; printf(n X1 is %f,X2 is %f, H is %f.,x1,x2,H); while(fabs(H)EP);min=pow(x1,2)+2*pow(x2,