黄金分割法机械的优化设计C语言程序

1、黄金分割法的优化设计实验报告学院：机电工程 机制自动化11-03班 学号：541102010326 姓名：刘点点1，黄金分割法的程序流程图开始rt?c2| 口2儿T 2i1口"2 口十久(bcf)2 /(2)是b 勺c?2否1-(a + b)结束）是否給定u、b、E1A 0.618ib A(b d)儿一于（4）+A(/> a)夕"2 /（勺）2,对应流程图的 C语言程序F面应用 C语言程序利用黄金分割法求一元函数F=xA2+2*x的最优解，已知初始区间为-3,5,取收敛精度e=10-4。C语言程序如下: #in clude<stdio.h>#in clud

2、e<math.h>#defi ne f(x) pow(x,2)+2*x#defi ne M 0.618void mai n()double y1,y2,x1,x2,x,a,b,e;int n;n=1;prin tf("请输入收敛精度e=");scan f("%lf", &e);prin tf("请输入区间左值 a=");scan f("%lf",&a);prin tf("请输入区间右值 b=");scan f("%lf",&b);prin

3、tf("na b x1 x2 y1 y2n");x仁 b-M*(b-a);x2=a+M*(b-a);y1= f(x1);y2=f(x2);prin tf("%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lfn", n,a,b,x1,x2,y1,y2); n=n+;doif(y1<y2)b=x2;x2=x1;y2=y1;x仁 b-M*(b-a);y1= f(x1);printf("%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lfn",n,a,b,x1,x2,y1,y2); n=

4、n+;elsea=x1;x1=x2;y仁 y2;x2=a+M*(b-a);y2=f(x2);prin tf("%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lfn", n,a,b,x1,x2,y1,y2); n=n+;while(fabs(b-a)/b)>=e&& fabs(y2-y1)/y2)>=e);x=(a+b)*0.5;prin tf("x=%.5lfn",x);getchar();3.运行结果:区区AAA有青主冃-日.日'b青左右e =0.0301 a=-3b=5xl03000.05

5、60x2口1-94400.11517,66?1-3.00001.94400.1151-3.00000-0560-1_8326-0_3068-0.9876-1 -83260-05&0-0-987&-0-8881-1.8326-0.6654-1.3868-0.«504-0.9876-1-3868-M.9410-0-9876-1.11140 9410-0.8358-0,?730-1.1114-0.8358-1.0061-0.9410-1-0000-0.9965-1 .1114-0.9410-1.0463-0-997911llll-1-0463-1-0214-0.9410-0.9812-0-9812-1.0061-0.981209996-1.0214-1.0061-0.9995-1.0000-1.0061-0.9966-1.0000-1.0000reas any keu to continue假定经十二次迭代后已满足收敛精度要求，则得x*=1/2(a