运筹MATLAB实验报告~~.docVIP

《运筹学》Matlab上机 实验报告 2009-2010学年第二学期 学院（部） 管理学院 学 号 031508107 姓 名 张思思 老 师 张伯生 成 绩 实验（一）（二） 1.实验目的： （1）加强学生分析问题的能力，锻炼数学建模的能力。 （2）掌握MATLAB软件中线性规划问题的编程、灵敏度问题的编程、求解和分析。 （3）利用所学MATLAB语言，完成线性规划问题的程序设计。 2。实验内容： （1）线性规划的单纯形法。 （2）线性规划的灵敏度分析。 3 实验步骤与方法（即案例分析）： 题目：例1.11 上海地铁公司正准备增加其人民广场的往来班次，因此需要雇佣更多的工作人员，但是不知道到底雇佣多少数量的工作人员。管理层意识到在向中介公司的客户提供令人满意的服务的同时必须进行成本管理，因此，必须寻求成本与收益之间合意的平衡。于是，要求管理团队研究如何规划人员才能以最小的成本提供令人满意的服务。 时段 x1 x2 x3 x4 x5 x6 x7 x8 6:00-8:00 x11 0 0 0 0 0 0 0 51 8:00-10:00 x21 x22 0 0 0 0 0 0 89 10:00-12:00 x31 x32 x33 0 0 0 0 0 75 12:00-14:00 x41 x42 x43 x44 0 0 0 0 97 14:00-16:00 x51 x51 x53 x54 x55 0 0 0 84 16:00-18:00 0 x62 x63 x64 x65 x66 x67 0 92 18:00-20:00 0 0 x73 x74 x75 x76 x77 x78 88 20:00-22:00 0 0 0 0 0 x86 x87 x88 75 22:00-24:00 0 0 0 0 0 0 x97 x98 43 24:00-0:00 0 0 0 0 0 0 0 x10.841 代理商成本 200 190 180 175 170 170 180 190 模型和截图如下： %m1.m文件：输入并化标准型 m=input(请输入约束条件个数：); n=input(请输入变量个数：); f=input(请输入目标函数系数矩阵：); a=input(请输入约束条件系数矩阵：); b=input(请输入约束条件右端常数矩阵:); f1=[f,zeros(1,m)]; a1=[a,eye(m,m)]; cb=a1(:,n+1:n+m); c=f1(:,n+1:n+m); mm=zeros(1,m+n); jie=[zeros(1,n),b]; disp(初始基可行解为:); jie dd=1:n+m; ee=dd(:,n+1:n+m); m2 %m2.m文件：计算检验数 q=zeros(1,n+m); for i=1:m+n ff=f1(:,i); ss=sum(c*a1(:,i)) t=ff-ss; q(:,i)=t; end disp (各检验数为); disp (q) m3 %m3.m文件判断函数是否有界 for i=1:m+n if q(:,i)0 t=i; if max(a1(:,t))=0 disp(此问题无界);break; else disp(此问题有界); end end end m4 %m4.m判断最优解 g=q(:,1); for i=1:n+m qq=q(:,i) if g=qq g=qq;z=i; end end if g=0 A=目标函数已得到最优值; disp(A); mm=mm; u=sum(f1*mm); disp(u); disp(最优解为：); mm=mm; break; else g=g; end disp(检验数最大值为：); disp(g); disp(换出变量x为：); disp(z); m5 %m5.m文件计算检验数 rr=a1(:,z); ct=b./rr; disp(检验数为); ct m6 %m6.m确定换出变量 h=100000; for j=1:m ctt=ct(j,:) if ctt0 if hctt h=ct(j,:);y=j+n; end end end disp(换出变量为); y m7 %m7.m迭代产生新的单纯形表 if yn cb(:,y-n)=