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

: .
黄金分割法的优化设计
实验报告
学院：机电工程
机制自动化11-03班
学号：541102: .
黄金分割法的优化设计
实验报告
学院：机电工程
机制自动化11-03班
学号：541102010326
姓名：刘点点
1，黄金分割法的程序流程图
开始给定Q、h、总
A—
b——d)
儿一/2i）
n+乂(h—£/)
J‘2J\先)
1、
a
1
1
b—口2
y2—V]
1t
O+/(/?—tf)y2-na2)
o[—力_/（丹_&）儿一心】）
I
是
b言⑺+/0
2,对应流程图的C语言程序F面应用C语言程序利用黄金分割法求一元函数
F=xA2+2*x的最优解，已知初始区间为
[-3,5],取收敛精度e=10-4。
C语言程序如下:
#include<>
#include<>
#definef(x)pow(x,2)+2*x
#
voidmain()
{doubley1,y2,x1,x2,x,a,b,e;
intn;n=1;
printf("请输入收敛精度e=");scanf("%lf",&e);printf("请输入区间左值a=");scanf("%lf",&a);
printf("请输入区间右值b=");scanf("%lf",&b);y2\n");
y2\n");
%.4lf\n",n,a,b,x1,x2,y1,y2);
printf("nabx1x2y1x1=b-M*(b-a);x2=a+M*(b-a);
y1=f(x1);y2=f(x2);
printf("%d%.4lf%.4lf%.4lf%.4lf%.4lfn=n++;do
{if(y1<y2)
{b=x2;
x2=x1;y2=y1;x1=b-M*(b-a);
y1=f(x1);printf("%d%.4lf%.4lf%.4lf%.4lf%.4lf%.4lf\n",n,a,b,x1,x2,y1,y2);
n=n++;}
else{
a=x1;x1=x2;
y1=y2;x2=a+M*(b-a);
y2=f(x2);printf("%d%.4lf%.4lf%.4lf%.4lf%.4lf%.4lf\n",n,a,b,x1,x2,y1,y2);
n=n++;}
}while(fabs((b-a)/b)>=e&&fabs((y2-y1)/y2)>=e);x=(a+b)*;
printf("x=%.5lf\n",x);getchar();}
AAAb二--Hu.---UJJ冃丄戸主fl
:
凤昭yr吕绅左区囘右x2yi
i・94400・1151
x2yi
i・9440