MATLAB基础性实验的报告3.docVIP

学 生 实 验 报 告  一、实验目的 学习最优化技术和基本原理，了解最优化问题的分类；掌握线性规划的建模技巧和求解方法；熟悉MATLAB软件求解线性规划模型的基本命令；通过范例学习，熟悉建立线性规划模型的基本要素和求解方法。 通过该实验的学习，使学生掌握最优化技术，认识面对现实生活中的最优化问题提出假设和建立优化模型，并且学会使用MATLAB软件进行线性规划模型求解的基本命令。 二、三、实验内容1．最优化问题的提出，提出不同的假设可以建立不同的最优化模型； 2．建立线性规划模型的基本要素和步骤； 3．使用MATLAB命令对线性规划模型进行计算。 四、实验步骤 1．开启MATLAB软件平台，开启MATLAB编辑窗口； 2．根据问题，建立无约束或非线性规划模型，并编写求解规划模型的M文件； 3．保存文件并运行； 4．观察运行结果(数值或图形)，并不断地改变参数设置观察运行结果； 5．根据观察到的结果和体会，写出实验报告。 五、实验要求与任务求解无约束优化 1) 画出该曲面图形直观地判断该函数的最优解2) 使用fminunc命令求解，能否求到全局最优解[x,y]=meshgrid(-5:0.01:5,-5:0.01:5); z=(-20)*exp((-0.2)*sqrt(0.5*(x.^2+y.^2)))-exp(0.5*(cos(2*pi*x)+cos(2*pi*y)))+22.713; mesh(x,y,z) 等高线图： [x,y]=meshgrid(-5:0.01:5,-5:0.01:5); z=(-20)*exp((-0.2)*sqrt(0.5*(x.^2+y.^2)))-exp(0.5*(cos(2*pi*x)+cos(2*pi*y)))+22.713; contour(x,y,z,20) drawnow hold on plot(1,1,o); text(1,1,atart point) plot(-1,-1,o) text(-1,-1,solution) 直观地可以看出最优解为:（0，0，0） 2）主函数： function f=fun1(x) f=(-20)*exp((-0.2)*sqrt(0.5*(x(1).^2+x(2).^2)))-exp(0.5*(cos(2*pi*x(1))+cos(2*pi*x(2))))+22.713; 主程序：oldoptions=optimset(fminunc) options=optimset(oldoptions,largescale,off) options11=optimset(options,HessUpdate,dfp) [x11,fval11,exitflag11,output11]=fminunc(fun1,[1 1],options11) pause x11 = 0 0 fval11 = -0.0053 exitflag11 = 5 output11 = iterations: 1 funcCount: 72 stepsize: 0.6107 firstorderopt: 2.8284 algorithm: medium-scale: Quasi-Newton line search message: [1x362 char] 求解非线性规划 试判定你所求到的解是否是最优？function f=fun2(x) f=-(0.201*x(1).^4*x(2)*x(3).^2)/(10.^7) 建立M文件mycon.m定义非线性约束： function [g,ceq] =mycon(x) ceq=[ ]; g=[x(1).^2*x(2)-675;(x(1).^2*x(3).^2)/(10.^7)-0.419]; 主程序为 x0=[18;4;100]; VLB=[ 0 0 0 ];VUB=[ 36 5 125 ]; [x,fval,exitflag,output]=fmincon(fun2,x0,[ ],[ ],[ ],[ ],VLB,VUB,mycon) fval=-fval 结果： x = 20.3378 1.6319 100.6475 fval = -56.8478 exitflag = 4 output = iterations: 4 funcCount: 16 lssteplength: 1 stepsize: 9.0469e-007 algorithm: medium-sca