中科大《优化设计》课程大作业之一维搜索优化设计实验报告.doc





《中科大《优化设计》课程大作业之一维搜索优化设计实验报告.doc》由会员分享,可在线阅读,更多相关《中科大《优化设计》课程大作业之一维搜索优化设计实验报告.doc(15页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、中科大优化设计课程大作业之一维搜索优化设计实验报告优化设计实验报告力学系型号:联想y470CPU:i52450M内存:2GB系统:win7-64位如下是三个目标函数(包括自定义函数)以及初值和精度选取:1、min f(x)=xx*xx4*xxx-6*x*x16x+4 初值=3; 精度=0。0000012、min f(x)=x*x+exp(-x) 初值=100; 精度=0.0000013、min f(x)=(x-1)*(x-1)+1初值= 3; 精度=0。000001如下是各个方法运算结果:1.黄金分割法函数时间次数极值点极值10.007099334.0-156。020。006005380。35
2、1810.8271830。000718311.01.02。平分法函数时间次数极值点极值10.13946104。0156。021.520886260。351730.8271831。319126211。01。03.成功失败法函数时间次数极值点极值10。015909124。0-156.020.002930390.351730.8271830。001100421。01.04.牛顿法函数时间次数极值点极值10.52764344。0156。020.63695730。351730。8271830。15020401。01.05。三点二次插值法函数时间次数极值点极值10.01416434。0-156。020.0
3、10109200。351730。8271830。00133901。01。0简要总结:通过上面的5个表格可以综合看出,牛顿法是这几个方法中最有效的方法。同时,根据编程来看,牛顿法的算法简单,所以程序相对其他几种简洁,故牛顿法可作为常用方法使用。在程序的运行方面,分别设置了可变的函数选择、起始点、方法选择、精度这四个输入,故可以在命令窗口运行主程序main,再根据提示分别输入这四个参数的所需值,就可以得到运行结果.程序如下:1、 主程序clear;global k;k=0;disp(1。f(x)=xx*xx-4*xxx6x*x-16*x+4);disp(2.f(x)=x*x+exp(x);disp
4、(3.f(x)=(x-1)(x-1)+1);while 1 no=input(请输入上面所想选择函数的编号(1、2、3):); if no=1|no=2|no=3 break; end disp(此次输入无效.);end disp( );x0=input(请输入初始值:);disp( );disp(1。黄金分割法);disp(2。平分法);disp(3。成功失败法);disp(4.牛顿法);disp(5。三点二次插值法);disp(6.三次插值法);while 1 m1=input(请输入上面所想选择方法的编号(1、2、3、4、5、6):); if m1=1m1=2|m1=3m1=4|m1=5
5、|m1=6 break; end disp(此次输入无效.);end disp( );while 1e=input(请输入精度(建议0。001或0.000001):); if e0 break; end disp(此次输入无效.);end disp( );disp();x,y=fmin(no,x0,m1,e);fprintf(迭代次数为: %8.0fn, k);fprintf(所求极值的横坐标x为: %16.5fn, x);fprintf(所求极值的纵坐标y为: %16。5fn, y);2、 调用函数function x,y=fmin(no,x0,m1,e)amin,amax,c=range1
6、(no,x0);if m1=1 tic;x,y=gold(no,amin,amax,e);toc;elseif m1=2 tic;x,y=pingfen(no,amin,amax,e);toc;elseif m1=3 tic;x,y=chenggong(no,amin,e);toc;elseif m1=4 tic;x,y=newton(no,amax,e);toc;elseif m1=5 tic;x,y=chazhi(no,amin,amax,c,e);toc;elseif m1=6 tic;x,y=sanci(no,amin,amax,e);toc;endend3、三点二次插值functio
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 优化设计 中科大 优化 设计 课程 作业 之一 搜索 实验 报告

限制150内