黄金分割法程序
黄金分割法
黄金分割法是通过不断缩短搜索区间的长度来寻求一维函数的极小点,这种方法的基本原理是:在搜索区间[a,b]内按如下规则对称地取两点a1和a2
a1=a+0.382(b-a); a2=a+0.618(b-a);
黄金分割法的搜索过程是:
给出初始搜索区间 [a,b] 及收敛精度 e ,将 赋以0.618
计算a1 和 a2,并计算起对应的函数值 f(a1),f(a2); ,
根据期间消去法原理缩短搜索区间,为了能用原来的坐标点计算公式,需进行区间名城的代换,并在保留区间中计算一个新的试验点及其函数值。
检查区间是否缩短到足够小和函数值收敛到足够近,如果条件不满足则返回到步骤2。
如果条件满足,则取最后两试验点的平均值作为极小点的数值近似解。
黄金分割法的流程图及程序清单
需要说明的是搜索区间[a,b]不需要给定,只需输入搜索精度e;
程序由四个子程序构成;
(1):输入输出子程序io();
(2):float fc (float x)求输入函数在某一点的值;
(3)void findqujian(float a[3],float f[3])确定搜索区间;
(4):float xunyou(float *value)寻找最小值
黄金分割法程序运行截图
#include "iostream.h"
#include "math.h"
#include "stdio.h"
#include "conio.h"
#define steplength 0.01
#define n 5
float e;float a,b,c,d,g;
float q[5];
void io()
{
cout<
cout<
cout<
printf("a=");
scanf("%f",&a);q[4]=a;
printf("b=");
scanf("%f",&b);q[3]=b;
printf("c=");
scanf("%f",&c);q[2]=c;
printf("d=");
scanf("%f",&d);q[1]=d;
printf("g=");scanf("%f",&g);q[0]=g;
cout<
cout<
scanf("%f",&e);
cout<
}
float fc(float x)
{int i;
float u=q[n-1];
for(i=n-2;i>=0;i--)
u=u*x+q[i];
return u;
}
void findqujian(float a[3],float f[3])
{
float t=float(steplength), a1,f1,ia;
a[0]=0;
f[0]=fc(a[0]);
for(int i=0;;i++)
{
a[1]=a[0]+t;f[1]=fc(a[1]);
if(f[1]
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]=fc(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 xunyou(float *value)
{float a1[3],f1[3],a[4],f[4];
float aa;
findqujian(a1,f1);
a[0]=a1[0];f[0]=f1[0];
a[3]=a1[2];f[3]=f1[2];
a[1]=a[0]+float(0.382)*(a[3]-a[0]);a[2]=a[0]+float(0.618)*(a[3]-a[0]);
f[1]=fc(a[1]);f[2]=fc(a[2]);
for(int i=0;;i++)
{if(f[1]>=f[2])
{a[0]=a[1];f[0]=f[1];