《实验二MATLAB编程单纯形法求解》由会员分享,可在线阅读,更多相关《实验二MATLAB编程单纯形法求解(11页珍藏版)》请在人人文库网上搜索。
1、北京联合大学 实验报告 项目名称: 运筹学专题实验报告 学 院: 自动化 专 业: 物流工程 班 级: 1201B 学 号:2012100358081 姓 名: 管水城 成 绩: 2015 年 5 月 6 日 实验二:MATLAB编程单纯形法求解 一、实验目的: (1)使学生在程序设计方面得到进一步的训练;,掌握Matlab (C或VB)语言进行程序设计中一些常用方法。 (2)使学生对线性规划的单纯形法有更深的理解. 二、实验用仪器设备、器材或软件环境 计算机, Matlab R2006 三、算法步骤、计算框图、计算程序等 本实验主要编写如下线性规划问题的计算程序: mincxAx?b ?s.。
2、t.?x?0,b?0?其中初始可行基为松弛变量对应的列组成. 对于一般标准线性规划问题: mincxAx?b ?s.t.?x?0,b?0?求解上述一般标准线性规划的单纯形算法(修正)步骤如下: 对于一般的标准形式线性规划问题(求极小问题),首先给定一个初始基本可行解。设初始基为B,然后执行如下步骤: ?1b?Bx令x?0,计算目标函数值f?cxbx?B 解(1).,,求得BBNBB ?1b的第i个分量2,.,m)记B(以bi?1, i?1CwB?w?CB,(2).计算单纯形乘子w, 得到,对于非基变量,计算判别BB?1cB1? 令,可直接计算数cA?cp?zc?cB?iiiiiBB?max?,。
3、R 为非基变量集合 kRi?0 ,若判别数则得到一个最优基本可行解,运算结束;否则,转到下一k步 ?1pBy?y0?ypy?B的每个分量均非正数, ,(3).解得到,;若即kkkkkk则停止计算,问题不存在有限最优解,否则,进行步骤(4).确定下标r,使 ? bb0y?min,且trx为离基变量, tkyyBtkrk0t?:y rtkx为进基变量,用p替换p,得到新的基矩阵B,还回步骤(1) ;Bkkr 、计算框图为: 令xB 计算单纯性乘子?令k 开始x B?)(非基变量c,j?R?wp?jjj初始可行基B 1?b,x?Bf?c0,B?N1?B?c,计算判别数wB?Rjmax?,j是?0? 。
4、k 否 得到最优 是1?,?Bp解方程By?p,得到y kkkk?0y? k 否 使得确定下标r,?bbir0y?min|?ikyy?ikrkx为退基变量,Brp代替k? ?不存在有限进基变量,以xk B得到新的基矩阵,pBr图1 :3计算程序(Matlab) );A=input(A= b=);b=input( );c=input(c= %ratformat 可以让结果用分数输出 m,n=size(A); E=1:m;E=E; F=n-m+1:n;F=F; D=E,F; %创建一个一一映射,为了结果能够标准输出 X=zeros(1,n); %初始化X if(n0) & (b1(i)/(A(i,。
5、k)+eps)temp ) temp=b1(i)/A(i,k); %找退基变量 r=i; end end fprintf(x(%d)进基,x(%d)退基n,k,D(r,2); %显示进基变量和退基变量 B(:,r)=A(:,k); cB(r)=c(k); %确定进基退基变量后,相应的基矩阵及新基对应的目标值的c也相应改变 D(r,2)=k; %改变D中的映射关系 end end end end 程序保存为 danchunxin.m 文件 四数值实验及其结果: 打开matlab软件,点击运行danchunxin.m,出现命令符要求输入相应矩阵命令。 1. 求解: 输入数据矩阵如下: A=9 4 。
6、1 0 0;4 5 0 1 0;3 10 0 0 1 b=360 200 300 c=-7 -12 0 0 0 点击运行得如下图: 2图最大由实验结果可知,该问题的最优解为:x1=100,x2=0,x3=540,x4=200,x5=0, 。700值为 2.求解: 输入数据矩阵如下: A=-3 -2 -1 -6 -18 1 0;-1 -0.5 -0.2 -2 -0.5 0 1;0.5 1 0.2 2 0.8 0 0 b=-700 -30 200 c=2 7 4 9 5 0 0 点击运行得如下图: 图3 由实验结果可知,该问题的最优解为:x1=15200/33,x2=0,x3=0,x4=0,x5=。
7、1250/33,x6=0,最小值为:1235/3。 3. 求解: 输入数据矩阵如下: A=1 -2 1 1 0;4 -1 -2 0 1;-2 0 1 0 0 b=11 -3 1 c=3 -1 -1 0 0 点击运行得如下图: 图4 由实验结果可知,该线性问题不存在最优解。 五实验总结 通过该实验,加深了对线性规划问题的理解,并加深了将实际问题进行数学建模转化为可利用MATLAB软件处理的线性规划问题,通过使用该软件,能充分的将复杂的线性规划问题通过编写程序实现相应的逻辑求解线性规划问题,在实际操作过程中,大大的减少了线性规划问题的运算的复杂性。编写函数程序代码是遇到了很多问题,但在老师的帮助下,逐一的解决了,为进行下面的实验奠定了一定的基础。。