欢迎来到天天文库
浏览记录
ID:27829004
大小:195.70 KB
页数:6页
时间:2018-12-06
《黄金分割法调试程序1》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、黄金分割法求解/=8x3-2x:-7x+3的]=J小值,区间[-10,10],精度:0.001C程序:#include"stdio.h"#include"math.h"doubleeq(doublex){doubley;y二&0*x*x*x・2.0*x*x・7.0*x+3.0;returny;}voidmain(){doubleal,a2,a3,a4,f2,f3,x;al=-10.0;//左边界a4=10.0;//右边界Labi:a3=al+0.618*(a4・al);f3=eq(a3);Lab2:a2
2、=al+0.382*(a4-al);f2=eq(a2);Lab3:if(fabs(a4-al)<0.001){x=(al+a4)/2.0;printf("x=%gfmin=%e",x,eq(x));}else{if(f2⑶{al=a2;a2=a3;f2=f3;a3=al+0.618*(a4-al);f3=eq(a3);gotoLab3;};};}运行结果
3、:编程安装CYuYanbinwte>p.exex=-9.99959fmin=-8.126OO5e+OO3Pressanykeytocontinue」J——C078/宙11新建晶保存剪切复制粘贴查找注释睡撤销重复运行调试校正对憑VC6.0筆试机试教7c,C/C++程序设计学习与实验系统文件0)编辑0)运行®工具©帮助0)信X-我的程序0.618c.c+较件应用问题解答+简单的入门程序实例+C语言学习方法指导+C/C卄语法参考+VC卄错误代码及程序实+函数查询+典型源程序+常用自定义函数算法+课程设
4、计源程序TurboC2.0常见错误令TurboC卄3.0常见错乍VisualC++6・0常见惜C语言常见专业词汇分Z运算符的优先级别次序常见控需悖符表(免费)ASCII码宇符对照表(竟nonameO.c01020304050607080910/*Note:VourchoiceisCIDE*/^includei,stdio.h"ttinclude,inath.hiidoubleeq(doublex){doubley;y=8.0*x»x*x-2.0«xreturny;uoidnain()1213141516
5、17181920212223242526272829303132ccdoublea1=a4=Labi:a3-f3=Lab2:a2-F2・Lab3:a1,a2»a3,a4,F2,f3,x;-10.0;〃左边界10.0;//右边界a1•0.618»(a4-a1);eq(a3);a1♦0.382•(aU-a1);eq(a2);if(Fabs(a4-a1)<0.001)6、F2(F2(F2a3二F2;gotoLab2;};F3)7、10;%左边界a4=10;%右边界e=0.001;%精度r=(sqrt(5)-1)/2;%艮卩0.618a2=al+(l・r)*(a4・al);a3=al+r*(a4-al);f2=func(a2);f3=func(a3);k=1;while(abs(a4-al)>=e)讦f28、=(al+a4)/2f=func(x)a3=al+0.618*(a4-al)f3=func(a3)Na4・al9、>=eY黄金分割算法框图运行结果:f=-96.8658f=80.5746f=-1.1925e+003f=-2.8893e+003f=_4.4795e+003f二_5.7048e+003f=-6.5633e+003f=-7.1345e+003f=-7.5035e+003f=-7.7377e+003f=-7.8849e+003f=-7
6、F2(F2(F2a3二F2;gotoLab2;};F3)7、10;%左边界a4=10;%右边界e=0.001;%精度r=(sqrt(5)-1)/2;%艮卩0.618a2=al+(l・r)*(a4・al);a3=al+r*(a4-al);f2=func(a2);f3=func(a3);k=1;while(abs(a4-al)>=e)讦f28、=(al+a4)/2f=func(x)a3=al+0.618*(a4-al)f3=func(a3)Na4・al9、>=eY黄金分割算法框图运行结果:f=-96.8658f=80.5746f=-1.1925e+003f=-2.8893e+003f=_4.4795e+003f二_5.7048e+003f=-6.5633e+003f=-7.1345e+003f=-7.5035e+003f=-7.7377e+003f=-7.8849e+003f=-7
7、10;%左边界a4=10;%右边界e=0.001;%精度r=(sqrt(5)-1)/2;%艮卩0.618a2=al+(l・r)*(a4・al);a3=al+r*(a4-al);f2=func(a2);f3=func(a3);k=1;while(abs(a4-al)>=e)讦f28、=(al+a4)/2f=func(x)a3=al+0.618*(a4-al)f3=func(a3)Na4・al9、>=eY黄金分割算法框图运行结果:f=-96.8658f=80.5746f=-1.1925e+003f=-2.8893e+003f=_4.4795e+003f二_5.7048e+003f=-6.5633e+003f=-7.1345e+003f=-7.5035e+003f=-7.7377e+003f=-7.8849e+003f=-7
8、=(al+a4)/2f=func(x)a3=al+0.618*(a4-al)f3=func(a3)Na4・al
9、>=eY黄金分割算法框图运行结果:f=-96.8658f=80.5746f=-1.1925e+003f=-2.8893e+003f=_4.4795e+003f二_5.7048e+003f=-6.5633e+003f=-7.1345e+003f=-7.5035e+003f=-7.7377e+003f=-7.8849e+003f=-7
此文档下载收益归作者所有