1、黄金分割法的优化设计实验报告学院:机电工程机制自动化 11-03 班学号:541102010326姓名:刘点点1,黄金分割法的程序流程图2,对应流程图的 C 语言程序下面应用 C 语言程序利用黄金分割法求一元函数 F=x2+2*x 的最优解,已知初始区间为-3,5 ,取收敛精度 e=10-4。C 语言程序如下:#include #include#define f(x) pow(x,2)+2*x #define M 0.618 void main() double y1,y2,x1,x2,x,a,b,e; int n; n=1; printf(“请输入收敛精度 e=“); scanf(“%lf“
2、, printf(“请输入区间左值 a=“); scanf(“%lf“, printf(“请输入区间右值 b=“); scanf(“%lf“, printf(“n a b x1 x2 y1 y2n“);x1=b-M*(b-a);x2=a+M*(b-a);y1=f(x1);y2=f(x2); printf(“%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lfn“,n,a,b,x1,x2,y1,y2); n=n+;do if(y1=e x=(a+b)*0.5; printf(“x=%.5lfn“,x); getchar(); 3.运行结果:假定经十二次迭代后已满足收敛精度要求,则得 x*=1/2(a+b)=1/2(-1.0214-0.9812)=-1.0013,相应的函数极值 f(x*)=-0.9999;近似精确值 x*=-1,f(x*)=-1,与解析法求得的精确值相同