原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
黄金分割法的优化设计
实验报告
学院：机电工程
机制自动化11—03班
学号：541102010326
姓名：刘点点
1，黄金分割法的程序流程图
2，对应流程图的C语言程序
下面应用C语言程序利用黄金分割法求一元函数F=x^2+2*x的最优解,已知初始区间为[—3，
5］,取收敛精度e=10-4。
C语言程序如下:
#include〈stdio.h
＃include〈math.h
＃definef(x）pow（x，2)+2*x
＃defineM0。618
voidmain(）
{
doubley1，y2,x1，x2,x，a，b,e;
intn；
n=1;
printf(”请输入收敛精度”请输入收敛精度）;
scanf(”％％，e);
请输入区间左值请输入区间左值）;
scanf（”%lf”,＆a)；
printf（（请输入区间右值b=”）;
scanf（（＆b)；
printf（”nabx1x2y1y2y2n”）;
x1=b—M＊(b—a）；
x2=a+M*（b-a）;
y1=f（x1)；
y2=f（x2）;
printf(”%d％.4lf％。4lf％.4lf%。4lf%。4lfn”，n,a,b,x1，x2,y1,y2)；
n=n++；
do
｛
if（y1y2)
{
b=x2;
x2=x1；
y2=y1;
x1=b—M＊（b—a）；
y1=f（x1)；
printf（”%d％。4lf%。4lf%.4lf%。4lf%.4lfn”，n,a,b，x1,x2,y1，
y2）；
n=n++;
｝
else
｛
a=x1;
x1=x2;
y1=y2；
x2=a+M＊(b-a）；
y2=f（x2);
％d%。4lf%。4lf%.4lf%.4lf％.4lf%。。n”，n，a，b，x1，x2,y1，
y2）;
n=n++;
｝
}while(fabs((b—a)/b）=efabs（（y2—y1)/y2)〉=e）；
x=（a+b）*0。5;
n”,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,与解析法求得的精确值相同
文档评论（0）