黄金分割法-C语言

本文介绍了一种使用C语言实现的黄金分割法求解一元函数最小值的方法。通过输入搜索区间的左右端值及精度,算法逐步逼近最优解。示例中采用的函数为f(x) = x^3 - 2*x + 1。

黄金分割法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.


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值