运筹学C-实验报告.docVIP

实 验 报 告 《运筹学C》 2012～2013学年第 二 学期 学院（部） 管 理 学 院 指导教师 陈克东 班级代号 031B112 姓名/学号 同组人 无 提交时间 2013.06.07 成绩评定 实验目的： 加强学生分析问题的能力，锻炼数学建模的能力。 掌握MATLAB软件中线性规划、灵敏度问题的编程、求解和分析。 利用所学MATLAB语言，完成运输问题程序设计。 用 WORD 书写实验报告：包括详细规划模型、试验步骤和结果分析。 实验内容： 题1： 某厂的一个车间有，两个工段可以生产三种产品，各工段开工一天生产三种产品的数量和成本，以及合同对三种产品的每周最低需求量由下表给出。问每周各工段对该生产任务应开工几天，可使生产合同的要求得到满足，并使成本最低。建立模型。 建立模型： 程序代码：A=[1,1;3,1;1,3]; B=[5,9,9]; C=[1000,2000]; [x,liuqiangning]=linprog(C,-A,-B,[ ],[ ],zeros(2,1)) 运行结果： 分析结果： 题2： 某工厂生产两种产品，分别需在A、B、C、D设备上加工，有关数据如下左表所示。它的标准型的最优表（、、、为松弛变量）如下右表所示。假定B设备增加10个台时所需费用为12，问增加B设备10个台时是否合算？若能增加利润，能增加多少？ 建立模型： 程序代码：A=[1,0,0,0;-1,0,-2,1/2;-1/4,1/4,1/2,-1/8;0,0,1,0] B=[22,18,16,28] C=A*B 运行结果： 分析结果： 题3 已知运输问题的产销平衡表与单位运价表如表所示，问怎样调配使总运费最小？试建立模型并进行求解。 销地 产地 B1 B2 B3 B4 产量 A1 A2 A3 3 11 3 10 1 9 2 8 7 4 10 5 7 4 9 销量 3 6 5 6 建立模型： 程序代码： %e运算主程序段 m=input(请输入产地数：); n=input(请输入销地数：); e=input(请输入单位运价表：); b1=input(请输入产量矩阵：); b2=input(请输入销量矩阵：); num1=sum(b1); num2=sum(b2); if num1==num2 [xx,fva]=trap1(m,n,e,b1,b2); elseif num1num2 [xx,fva]=trap2(m,n,e,b1,b2); else [xx,fva]=trap3(m,n,e,b1,b2); end p=round(xx); beq=[b1;b2]; A=调运运输表 disp(p); B=总运费： disp(fva); %trap1.产销平衡 function[xx,fv]=trap1(m,n,e,b1,b2) p=e;f=p(:); [aeq1,aeq2]=trb(m,n); aeq=[aeq1;aeq2]; beq=[b1;b2]; lb=zeros(m*n,1); [x,fval,exitflag,output,lambda]=linprog(f,[],[],aeq,beq,lb); fv=fval; xx=reshape(x,n,m); xx=xx; %trap2.产大于销 function[xx,fv]=trap(m,n,e,b1,b2) p=e;f=p(:); [aq1,aq2]=trb(m,n); a=aq1;b=b1; aeq=aq2;beq=b2; lb=zeros(m*n,1); [x,fval,exitflag]=linprog(f,a,b,aeq,beq,lb); fv