实验五2015在简约的世界里使收益最大数学规划.docVIP

开课学院、实验室：计算机学院 ds1422 实验时间 ：2016 年 4 月 10 日 课程 名称 数学实验 实验项目 名 称 实验五 在简约的世界里使收益最大——数学规划 实验项目类型 验证 演示 综合 设计 其他 指导 教师 成 绩 一、实验目的及意义 [1] 学习最优化技术和基本原理，了解最优化问题的分类； [2] 掌握线性规划的建模技巧和求解方法； [3] 学习灵敏度分析问题的思维方法； [4] 熟悉MATLAB软件求解线性规划模型的基本命令； [5] 通过范例学习，熟悉建立线性规划模型的基本要素和求解方法。 MATLAB软件进行线性规划模型求解的基本命令，并进行灵敏度分析。解决现实生活中的最优化问题是本科生学习阶段中一门重要的课程，因此，本实验对学生的学习尤为重要。 1．最优化问题的提出，提出不同的假设可以建立不同的最优化模型； 2．建立线性规划模型的基本要素和步骤； 3．使用MATLAB命令对线性规划模型进行计算与灵敏度分析； 4．利用优化数值解与图形解对最优化特征作定性与定量分析； 1．开启MATLAB软件平台，开启MATLAB编辑窗口； 2．根据问题，建立的线性规划模型，并编写求解规划模型的M文件； 3．保存文件并运行； 4．观察运行结果(数值或图形)，并不断地改变参数设置观察运行结果； 5．根据观察到的结果和体会，写出实验报告。 根据实验内容和步骤，完成以下实验，要求写出实验报告（实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论） 应用实验 1. 投资策略 某部门现有资金10万元，五年内有以下投资项目可供选择： A：从第一年到第四年每年初投资，次年末收回本金且获利15%； B：第三年初投资，第五年末收回本金且获利25%，最大投资额为4万元； C：第二年初投资，第五年末收回本金且获利40%，最大投资额为3万元； D：6%； 2求解无约束优化 1) 画出该曲面图形, 直观地判断该函数的最优解; 2) 使用fminunc命令求解, 能否求到全局最优解?. 求解非线性规划 问题的解决及源代码： 1. 1.1问题分析 用表示第i年对第j个项目的投资金额 要使第五年年末本息总额最大，应当在每年将所有可用资金都用于投资，以确保资金的充分利用，由于项目投资均发生在年初，故以下只讨论年初的投资情况： 第一年： 第二年：手上资金（即第一年年末收回资金）为，全部用来对可投资项目投资，则有= 第三年：同理，有= 第四年：= 第五年：= 第五年年末本息和为（即第五年所能收回的所有资金） 1.2建立模型： = = = = ， 1.3求解模型： 可编写lingo程序如下： Model: max=1.06*x54+1.15*x41+1.25*x32+1.4*x23; x11+x14=10; 1.06*x14=x21+x23+x24; 1.15*x11+1.06*x24=x31+x32+x34; 1.15*x21+1.06*x34=x41+x44; 1.15*x31+1.06*x44=x54; x23=3; x32=4; end 运行结果如下： 所得最优值为14.375万元，对应的最优解为: x11=7.169811, x14=2.830189, x23=3, x32=4, x34=4.245283, x41=4.5, 其余值为0 即第一年对A项目投资7.169811万元， 对D项目投资2.830189万元； 第二年对C项目投资3万元； 第三年对B项目投资4万元， 对D项目投资4.245283万元； 第四年对A项目投资4.5万元。 2. (1). X=-5:0.1:5; Y=X; [x1,x2]=meshgrid(X,Y); R=-0.2*sqrt(0.5*(x1.^2+x2.^2)); P=0.5*(cos(2*pi.*x1)+cos(2*pi.*x2)); r=-20*exp(R);p=-1*exp(P); yf=r+p+22.713; mesh(x1,x2,yf);title(函数曲面图形) (2). 子函数： function f=hw2(x) R=-0.2*sqrt(0.5*(x(1)^2+x(2)^2)); P=0.5*(cos(2*pi*x(1))+cos(2*pi*x(2))); r=-20*exp(R); p=-1*exp(P); f=r+p+22.713 主函数： x1=[-1 -1]; options=optimset(display,iter,tolfun,1e-10); [x,fval]=fminunc(hw2,x1,options) 结果: x = 0 0 fval =