loading
课 程 实 验 报 告
专 业 年 级 13级数学与应用数学 课 程 名 称 运筹学 指 导 教 师 文永松 学 生 姓 名 雷中
学 号 20132211021057 实 验 日 期 2015年12月17日 实 验 地 点 笃行楼B栋407 实 验 成 绩
教务处制 2015 年 12 月7日
实验项目名称 实验 目的及要求 函数作图 （1）熟练掌握matlab程序的基本操作 （2）为数学建模打好基础,提高自己的编程能力。 （3）学会解决实际问题中的最优化问题，为企业管理决策提供支持。 （4）简单的编程求解现实中的线性规划问题问题。 实验 内容 编写一个MATLAB的函数文件，用于求解标准形的线性规划问题： 1、模型： min z=cX .t.AX?bs命令：x=linprog（c，A，b） 2、模型： minz?cXs.t. AX?b AeqX?beq 命令：x=linprog（c，A，b，Aeq,beq） 注意：若没有不等式：AX?b存在，则令A=[ ]，b=[ ]. 若没有等式约束, 则令Aeq=[ ], beq=[ ]. 3、模型： minz?cXs.t. AX?b AeqX?beq VLB?X?VUB 命令：[1] x=linprog（c，A，b，Aeq,beq, VLB，VUB） [2] x=linprog（c，A，b，Aeq,beq, VLB，VUB, X0） 注意：[1] 若没有等式约束, 则令Aeq=[ ], beq=[ ]. [2]其中X0表示初始点 4、命令：[x,fval]=linprog(…) 返回最优解x及x处的目标函数值fv 实际具体题目： 一、用MATLAB做函数图。 （1）.y=sin(x)和y=cos(x),[-4π，4π] （2）.(3). 旋转抛物面 ， 实验步骤 算法：1.运行MATLAB.exe文件。 2. 输入MATLAB作图的程序代码。 对于一般的标准形式线性规划问题(求极小问题),首先给定一个初始基本可行解。设初始基为B,然后执行如下步骤： (1).解BxB?b?1x?Bb，令xN?0,计算目标函数值f?cBxB B,求得以bi(i?1,2,...,m)记B?1b的第i个分量 (2).计算单纯形乘子w, wB?CB?1w?CBB,得到,对于非基变量，zi?ci},R为非基变量集计算判别数?i?zi?ci?cBB?1pi?ci,令 ?k?max{i?R合 若判别数?k?0 ,则得到一个最优基本可行解，运算结束；否则，转到下一步。 实验环境 MATLAB 实验结果与 分析 计算及程序代码 一.用MATLAB做函数图像： 1.编写的M码为： 输出的函数图像： 2．编写的M码为： 输出的函数图像为：