黄金分割法求极值(c编程)

合集下载

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

流程图：

C程序:

#include <>

double f(double x) {double y;

y=*x**x*x+;

return(y);

}

void main()

{double a,b,x1,x2,w,s,min,x0;

int i=1;

double h;

w=(sqrt(5)-1)/2;

s=1e-5;

printf("input a,h:\n");

scanf("%lf,%lf",&a,&h);

do /*确定初始区域[a,b]*/

{b=a+i*h;

i++;

printf("i=%d,b=%lf\n",i,b);

}while (f(a)>f(b));

x1=a+(1-w)*(b-a);

x2=a+w*(b-a);

printf("[%lf,%lf]\n",a,b);

loop:if(f(x1)

if(fabs(b-a)>s)

{x2=x1;

x1=a+(1-w)*(b-a);

printf("x1=%lf,x2=%lf\n",x1,x2);

goto loop;

}

else f((a+b)/2);

}

else

{a=x1;

if(fabs(b-a)>s)

{x1=x2;

x2=a+w*(b-a);

printf("x1=%lf,x2=%lf\n",x1,x2);

goto loop;

}

else f((a+b)/2);

}

x0=(a+b)/2;

min=f(x0);

printf("x0=%lf,min=%lf\n",(a+b)/2,f((a+b)/2));

}

运算结果：