最优化实验报告.docVIP

广东金融学院实验报告 课程名称：最优化方法 实验编号 及实验名 称 共轨梯度法及二分法Mat lab编程求解实验 系另U 应用数学系 姓名 陈格 学 号 111613225 111613226 班级 1116132 实验地点 实验楼504 实验日期 2013年12月25日 实验时数 2 指导教师 刘伟 同组其他成 员 陈宇强 成绩 一、 实验目的及耍求 运用Mat lab软件编程实现共轨梯度法求解书本173页第10小题及二分法求解137页笫5小题。 二、 实验环境及相关情况（包含使用软件、实验设备、主要仪器及材料等） 使用软件:Mat lab 三、实验内容及步骤(包含简耍的实验步骤流程) 第一题用一般的二分法求解，第二题由于所求的函数是非二次函数，需要非二次 函数的共辘梯度法，程序中也需要optimset函数来求解。 1、 用二分法解 min k + 2：x 呂” — 3 ￡ Ji ￡ 6 最后长度*0.2。 实验步骤：1、先对x: + 2x求导，W-f(x) = 2x+ 2 2、 初始区间取[-3,6],最后区间长度为6 = 0.20即％ = —3 ‘加=6 d = 0.2 3、 每次区间长度记录为哲=务迭代终止条件兔￡ 0.2o 4、 取c* = + 代入f(x)屮，通过判断c： 大于0赋值给b, + 小于零赋值给a, 1 = 1 + 1； 等于0结束迭代，此吋c；为最优解。 2、 用共轨梯度法极小化函数 f(x) = 100(^2 一 兀』)2 + (1- = (-1切片 实验步骤：1、先对f(x)中的H，2求偏导分别为g⑴、g(2)o 得： g(l) = -400(^ X x2 - x^) - 2 + 3 X xt; S(2) = 200(xj -琉》 2、再调用optimset函数求最小值。 四、 实验结果（包括程序或图表、结论陈述、数据记录及分析等，可附页） 1、按步骤编写二分法mat lab程序及其结果: syms x; Z二X 2+2*x; % 口标原函数 y 二 diff(z); f二inline( 2*x+2‘，‘ x); d二0. 2; %要求误差 a=-3; %区间下界 b=6;%区间上界 e=b-a; k二0 ; while ed c 二(a+b)/2; if f(c)0 a=c; elseif f (c)0 b二c; else a二 c; b二 c end e=e/2; k=k+l; end x二(a+b)/2; x k 运行结果： x = -0. 9609 k = 6 迭代6次后得到最优值x二-0. 9609o 2、按步骤川共轨梯度法极小化函数： 首先建立y函数 function [f,g]=y(x) f= 100*(x(2)-x(l)A2)A2+(l -x(l))A2; g(l)=-400*x(2)*x(l )+400*x( 1 )A3-2+2*x( 1); g(2)=200*x(2)-200*x( 1)A2; 调川optimset函数： options=optimset(GradObj,on);[x fval]=fminunc(y1.2 1],options) x = 1.0000 1 .()()(X) fval = 4.0728e-013 与题Fl极小化的函数值0较为接近，兀=口4厂。 五、 实验总结（包括心得体会、问题回答及实验改进意见，可附页） 学会了用二分法与共辘梯度法通过niatlab程序，极小化函数，可以看到得岀的解不一定是 精确的，有可能只是1?分接近精确解（如第二题的用用共辘梯度法极小化函数），毕竞这些方法 只是迭代过程，是需要控制在一定的精度里，当小于精度吋会停止迭代，得出最终解。 六、教师评语 评语 评语等级 优 中 及格 不格 1 ?实验态度认真，实验目的明确 2.实验方案、程序设计合理 3.实验过程（实验步骤详细，记录完整，数 据合理） 4.实验结论正确，分析透彻 5?实验报告独立完成，无抄袭现彖，并按时 提交，格式规范，文字叙述流畅，逻辑性强 附录: