黄金分割法-C语言

黄金分割法C语言代码实现：

#include<stdio.h>

#include<stdlib.h>
#include<math.h>
#define cs1 0.382
#define cs2 0.618
int main()                                    
{
double fc(double);
double lmd[10],miu[10],flmd[10],fmiu[10],a[10],b[10];
int k=0;
double e;
printf("请输入区间a,b和精度e，并用逗号隔开。\n");
scanf("%lf,%lf,%lf",&a[0],&b[0],&e);
lmd[0]=a[0]+cs1*(b[0]-a[0]);
miu[0]=a[0]+cs2*(b[0]-a[0]);
flmd[0]=fc(lmd[0]),fmiu[0]=fc(miu[0]);
pd1:
printf("a[k]=%lf\t b[k]=%lf\t lmd[k]=%lf\t miu[k]=%lf\t flmd[k]=%lf\t fmiu[k]=%lf\t k=%d\n",a[k],b[k],lmd[k],miu[k],flmd[k],fmiu[k],k);
if(flmd[k]>fmiu[k])
{
if(b[k]-lmd[k]<=e)
{
printf("函数近似最优解x=%lf\n",miu[k]);
exit(0);
}
else
{
k++;
a[k]=lmd[k-1];
b[k]=b[k-1];
lmd[k]=miu[k-1];
fmiu[k-1]=fc(miu[k-1]);
flmd[k]=fmiu[k-1];
miu[k]=a[k]+cs2*(b[k]-a[k]);
fmiu[k]=fc(miu[k]);
goto pd1;
}
}
else
{
if(miu[k]-a[k]<=e)
{
printf("函数近似最优解为x=%lf\n",miu[k]);
exit(0);
}
else
{
k++;
a[k]=a[k-1];
b[k]=miu[k-1];
miu[k]=lmd[k-1];
flmd[k-1]=fc(lmd[k-1]);
fmiu[k]=flmd[k-1];
lmd[k]=a[k]+cs1*(b[k]-a[k]);
flmd[k]=fc(lmd[k]);
goto pd1;
}
}


return 0;
}
double fc(double x)
{
double y;
y=pow(x,3)-2*x+1;
return y;
}

a,b分别为搜索区间左右段值，e为区间精度。本次算法的函数为f(x)=x^3-2*x+1。

design by zam.