黄金分割法.cpp

c写的黄金分割法（0.618）最优化程序

#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<<"假设多项式的最高次幂为4"<<endl<<endl;
	cout<<"设多项式的一般形式为f=a*x^4+b*x^3+c*x^2+d*x+e"<<endl<<endl;
	cout<<"请输入要求解的目标多项式的系数"<<endl<<endl;
	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<<endl;
  cout<<"请输入搜索精度e"<<endl<<endl;
   scanf("%f",&e);
     cout<<endl;
}
 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]<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]=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];
a[1]=a[2];f[1]=f[2];
a[2]=a[0]+float(0.618)*(a[3]-a[0]);f[2]=fc(a[2]);
}
else
{a[3]=a[2];f[3]=f[2];
a[2]=a[1];f[2]=f[1];
a[1]=a[0]+float(0.382)*(a[3]-a[0]);f[1]=fc(a[1]);
}
if(fabs(a[3]-a[0])<e/2)
{aa=(a[1]+a[2])/2; *value=fc(aa);
break;}
}return(aa);
}


void main()
{float xx,value;
	io();
 xx=xunyou(&value);
printf("f(x)=%2.1f*x^4+%2.1f*x^3+%2.1f*x^2+%2.1f*x+%2.1f",a,b,c,d,g);
 cout<<"取得最小值的坐标为"<<endl<<endl;
printf("\nzuobiao*=%f\n\n",xx);
cout<<"函数的最小值是"<<endl;
printf("\nminmum*=%f\n",value);
getch();
}