欢迎来到天天文库
浏览记录
ID:51342819
大小:346.50 KB
页数:10页
时间:2020-03-22
《机械优化设计课程论文.doc》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、一、一维搜索F(x)=0.8x-5x-7x-l2x+170方法:黄金分割法编程语言:C++实现平台:VS2010给定h°y
2、黃金分割法框图收敛精度:0.0000001Ymin=-208.04609313进退法框图运算结果:进退法初始点:1()迭代次数:38源程序如下:步长:0.01极小点Xmin=5.59013095#include^include^include^include#include3、-5*x*x*x-*x*x-12*x+l;voidfind(double&*,double&b)//进]层法{doublexl,x2,x3zy1zy2fy3zStep,xOzh;inti=l;cout<<"•、J•找小lb筆区4、X5、”<>xO>>Step;cout<6、T<7、=y3){h=2*h;xl=x2;x2=x3;yi=y2;y2=y3;x3=x2+h;y3=f(x3);coutvvi++vvsetw()<8、(xlzx3);coutv区间为:"«""«a<>EPS;cout<<'!N,'<9、10、recision(4)<EPS){if(yl>=y2){a=xl;xl=x2;yl=y2;x2=a+k*(b-a);y2=f(x2);}else{b=x2;x2=xl;y2=y1;xl=b-k*(b-a);yl=f(xl);}err=fabs(b-a);couvj++vvsetw(0)<
3、-5*x*x*x-*x*x-12*x+l;voidfind(double&*,double&b)//进]层法{doublexl,x2,x3zy1zy2fy3zStep,xOzh;inti=l;cout<<"•、J•找小lb筆区
4、X
5、”<>xO>>Step;cout<6、T<7、=y3){h=2*h;xl=x2;x2=x3;yi=y2;y2=y3;x3=x2+h;y3=f(x3);coutvvi++vvsetw()<8、(xlzx3);coutv区间为:"«""«a<>EPS;cout<<'!N,'<9、10、recision(4)<EPS){if(yl>=y2){a=xl;xl=x2;yl=y2;x2=a+k*(b-a);y2=f(x2);}else{b=x2;x2=xl;y2=y1;xl=b-k*(b-a);yl=f(xl);}err=fabs(b-a);couvj++vvsetw(0)<
6、T<7、=y3){h=2*h;xl=x2;x2=x3;yi=y2;y2=y3;x3=x2+h;y3=f(x3);coutvvi++vvsetw()<8、(xlzx3);coutv区间为:"«""«a<>EPS;cout<<'!N,'<9、10、recision(4)<EPS){if(yl>=y2){a=xl;xl=x2;yl=y2;x2=a+k*(b-a);y2=f(x2);}else{b=x2;x2=xl;y2=y1;xl=b-k*(b-a);yl=f(xl);}err=fabs(b-a);couvj++vvsetw(0)<
7、=y3){h=2*h;xl=x2;x2=x3;yi=y2;y2=y3;x3=x2+h;y3=f(x3);coutvvi++vvsetw()<8、(xlzx3);coutv区间为:"«""«a<>EPS;cout<<'!N,'<9、10、recision(4)<EPS){if(yl>=y2){a=xl;xl=x2;yl=y2;x2=a+k*(b-a);y2=f(x2);}else{b=x2;x2=xl;y2=y1;xl=b-k*(b-a);yl=f(xl);}err=fabs(b-a);couvj++vvsetw(0)<
8、(xlzx3);coutv区间为:"«""«a<>EPS;cout<<'!N,'<9、10、recision(4)<EPS){if(yl>=y2){a=xl;xl=x2;yl=y2;x2=a+k*(b-a);y2=f(x2);}else{b=x2;x2=xl;y2=y1;xl=b-k*(b-a);yl=f(xl);}err=fabs(b-a);couvj++vvsetw(0)<
9、10、recision(4)<EPS){if(yl>=y2){a=xl;xl=x2;yl=y2;x2=a+k*(b-a);y2=f(x2);}else{b=x2;x2=xl;y2=y1;xl=b-k*(b-a);yl=f(xl);}err=fabs(b-a);couvj++vvsetw(0)<
10、recision(4)<EPS){if(yl>=y2){a=xl;xl=x2;yl=y2;x2=a+k*(b-a);y2=f(x2);}else{b=x2;x2=xl;y2=y1;xl=b-k*(b-a);yl=f(xl);}err=fabs(b-a);couvj++vvsetw(0)<
此文档下载收益归作者所有