安徽师范大学 数学建模实验课 (实验报告) 7-1无约束优化.docVIP

安徽师范大学 数计 学院实验报告 专业名称 数学与应用数学 实 验 室 6#203 实验课程 建模实验 实验名称 无约束优化 姓 名 夏 康 学 号 同组人员 无 实验日期 2014.5.13 注：实验报告应包含（实验目的，实验原理，主要仪器设备和材料，实验过程和步骤，实验原始数据记录和处理，实验结果和分析，成绩评定）等七项内容。具体内容可根据专业特点和实验性质略作调整，页面不够可附页。 实验目的 1、了解无约束优化的基本算法； 2、掌握Matlab优化工具箱的基本用法； 3、掌握用Matlab求解无约束优化实际问题。 实验原理 1、无约束优化的基本操作步骤 2、Matlab的基本操作步骤 主要仪器设备和材料 PC机 实验过程和步骤 2求解，不同初值，不同算法（搜索方向）、步长搜索、数值梯度与分析梯度等）的结果分析、比较。 根据求解梯度算法的函数，我们建立不同的matlab文件进行试验。首先建立函数文件 fun1.m %用于数值梯度 function f=fun1(x,a,b) f=x(1)^2/a+x(2)^2/b; fun2.m %用于分析梯度 function [f,g]=fun2(x,a,b) f=x(1)^2/a+x(2)^2/b; if nargout1 g(1)=2*x(1)/a;g(2)=2*x(2)/b; end 建立主程序 %jizhi3.m a=2;b=1;x0=[1,1]; format short e ---largescale [x,f,ef,out]=fminunc(@fun1,x0,[],a,b) ---case1:bfg---- opt1=optimset(LargeScale,off) [x1,v1,exit1,out1]=fminunc(fun1,x0,opt1,a,b); ---case2:dfp fopt=optimset(opt1,HessUpdate,dfp) [x2,v2,exit2,out2]=fminunc(fun1,x0,fopt,a,b) ---case3:steep fopt=optimset(opt1,HessUpdate,steepdesc); [x3,v3,exit3,out3]=fminunc(fun1,x0,fopt,a,b) ---case4:bfg1 opt2=optimset(LargeScale,off,LineSearchType,cubicpoly) [x4,v4,exit4,out4]=fminunc(fun1,x0,opt2,a,b) ---case5:dfp fopt2=optimset(opt2,HessUpdate,dfp,LineSearchType,cubicpoly) [x5,v5,exit5,out5]=fminunc(fun1,x0,fopt2,a,b) ---case6:steep fopt2=optimset(opt2,HessUpdate,steepdesc,LineSearchType,cubicpoly); [x6,v6,exit6,out6]=fminunc(fun1,x0,fopt2,a,b) +++++results of solutions++++ solutions=[x1;x2;x3;x4;x5;x6]; funvalues=[v1;v2;v3;v4;v5;v6]; iterations=[out1.funcCount;out2.funcCount;out3.funcCount;out4.funcCount;out5.funcCount;out6.funcCount]; [solutions,funvalues,iterations] ---case7:bfg---- opt1=optimset(LargeScale,off,GradObj,on) [x1,v1,exit1,out1]=fminunc(fun2,x0,opt1,a,b) ---case8:dfp fopt=optimset(opt1,HessUpdate,dfp) [x2,v2,exit2,out2]=fminunc(fun2,x0,fopt,a,b) ---case9:steep fopt=optimset(opt1,HessUpdate,steepdesc); [x3,v3,exit3,out3]=fminunc(fun2,x0,fopt,a,b) ---case10:bfg1 opt2=optimset(LargeScale,off,LineSearchType,cubicpoly,G