数学实验报告-非线性规划与多目标规划实验.doc

PAGE PAGE 6 重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1421 学 院 年级 专业班 学 生 姓 名 学 号 开 课 时 间 至 学年第 学期 总 成 绩 教师签名 数学与统计学院制 开课学院、实验室：数学与统计学院DS1421实验时间 ： 2013 年 5 月 15 日 课程 名称 数学实验 实验项目 名 称 非线性规划与多目标规划实验 实验项目类型 验证 演示 综合 设计 其他 指导 教师 成 绩 实验目的 [1] 学习非线性规划模型的标准形式和建模方法； [2] 掌握建立非线性规划模型的基本要素和求解方法； [3] 熟悉MATLAB软件求解非线性规划模型的基本命令； [4] 通过范例学习，了解建立非线性规划模型的全过程，与线性规划比较其难点何在。 本实验包括基础实验、应用实验和创新实验，基础实验和应用实验要求独立完成，创新实验要求合作完成。通过该实验的学习，使学生掌握最优化技术，认识面对什么样的实际问题，提出假设和建立优化模型，并且使学生学会使用MATLAB软件和Lingo软件求解非线性规划模型，注意初始解的选择不同会导致软件求出的解的变化（是局部最优解还是整体最优解）。解决现实生活中的最优化问题是本科生学习阶段中一门重要的课程，因此，本实验对学生的学习尤为重要。 基础实验一 一、实验内容 1求解无约束优化 1) 画出该曲面图形, 直观地判断该函数的最优解; 2) 使用fminunc命令求解, 能否求到全局最优解? 二、实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等， 实验结果及分析） 解：1.（1）编写绘图命令 x=-5:0.01:5; y=-5:0.01:5; [X,Y]=meshgrid(x,y); Z=-20.*exp(-0.2.*(0.5.*(X.^2+Y.^2).^0.5))-exp(0.5.*(cos(2.*pi.*X)+cos(2.*pi.*Y)))+22.713; mesh(X, Y, Z); 结果如下： 最优解近似在x=0，y=0处。 （2）用fminunc命令求解 编写函数文件： 1）建立函数M文件： function f=fun(x) f=-20*exp(-0.2*(0.5*(x(1)^2+x(2)^2)^0.5))-exp(0.5*(cos(2*pi*x(1))+cos(2*pi*x(2))))+22.713; 2） x0=[0,0]; options=optimset(display, iter, tolfun,1e-10); [x,fval]= fminunc(fun,x0, options) 结果： ，fval =-0.0053 与图形结果相对吻合 基础实验二 一、实验内容 2. 求解非线性规划,试判定你所求到的解是否是最优? 二、实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等， 实验结果及分析） 解： 编写函数文件： ①function f=fun(x) f=-1*0.201*x(1)^4*x(2)*x(3)^2/10^7; ②function [c,ceq]=nlcon(x) c(1)= x(1)^2*x(2)-675;c(2)= 1e-007*x(1)^2*x(2)^2-0.419; ceq=[]; 给定初值求解 ③x0=[10 10 2]; L=[0 0 0];U=[36,5,125]; [x,fval]=fmincon(fun2,x0,[],[],[],[],L,U,nlcon) fmax=-fval 结果： fval =-274.7419,fmax =274.7419。 应用实验（或综合实验） 实验内容 4. 组合投资问题 设有8种投资选择:5支股票,2种债券,黄金. 投资者收集到这些投资项目的年收益率的历史数据 (见表6.1), 投资者应如何分配他的投资资金,即需要确定这8种投资的最佳投资分配比例. 表6.1 8种投资项目的年收益率历史数据 项目 年份 债券1 债券2 股票1 股票2 股票3 股票4 股票5 黄金 1973 0.075 -0.058 -0.148 -0.185 -0.302 0.023 -0.149 0.677 1974 0.084 0.02 -0.265 -0.284 -0.338 0.002 -0.232 0.722 1975 0.061 0.