实验4 运输问题求解的编程实现.doc

实验4 运输问题求解的编程实现成绩专业班级 信息121 学号 201212030120 姓名 刘帅 报告日期 .实验类型：验证性实验 综合性实验 设计性实验实验目的：熟练运输问题求解。实验内容：运输问题求解表上作业法2个（产销平衡、不平衡各一个）实验原理 运用运输问题表上作业法求解，首先确定初始调运方案，然后计算空格的检验数，判定是否最优方案，如果不是，进行调整，然后再重复上述过程，直到求出最优解。实验步骤1 要求上机实验前先编写出程序代码 2 编辑录入程序3 调试程序并记录调试过程中出现的问题及修改程序的过程4 经反复调试后，运行程序并验证程序运行是否正确。5 记录运行时的输入和输出。 预习编写程序代码：实验报告：根据实验情况和结果撰写并递交实验报告。实验总结：参考程序 function metr X = LP_transport(Cost, Sc, De ) if nargin = 3 error(There must be 3 arguments );end nS nD = size(Cost) ;nDi = length(De) ;nSi = length(Sc);if nD = nDi error(Unmatched column number );elseif nS = nSi error(Unmatched row number );end clear nDi nSi ; if abs(sum(De) - sum(Sc) 1e-10 disp(unbalanced transportation problem is not supported.); X = 0; metr = 0; return;end Aeq = conv_constr_matrix( nS, nD ); X Bflag=northwest(Sc,De);metr = sum(sum(X.* Cost ); while (1) optflag entB= IsOptimal( Bflag, Cost, Aeq ) ; if optflag break; end Y Bout=cycle(X, entB(1), entB(2), Bflag); Bflag = Bout; X = Yend metr = sum(sum(X.* Cost ); function optflag entB= IsOptimal( Bflag, Cost, Aeq ) nS nD = size(Cost);cnum = nS * nD; Bflaglst = reshape(Bflag, 1, cnum );Clst = reshape(Cost, 1, cnum );Aeq = Aeq(1:(end-1), :); B = Aeq(:, logical(Bflaglst); %BaseCb = Clst( logical(Bflaglst) ); nidx = find( Bflaglst ) ;N = Aeq(:, nidx);Cn = Clst( nidx ); U = Cb/B; R = Cn - U*N; NegFlag = R Cost=4 7 3 8;2 9 2 8;9 4 10 5; Sc=9 4 9; De=4 6 5 6; LP_transport(Cost, Sc, De )X = 2 0 5 0 1 0 0 3 0 6 0 3ans =85最优解矩阵为：X。此时最小运费为85。二 所选例题：产销不平衡课本185(5-27)例题 需求地区化肥厂产量（万吨）A2613221760B1413191550C192023-60最低需求（万吨）最高需求（万吨）3050707003010不限 De=40 60 0 10;40 80 30 1000; Sc=60 50 60; Cost=26 13 12 17;14 13 19 15;19 20 23 1000; LP_transport(Cost,Sc,De)unbalanced transportation problem is not supported.不支持不平衡运输问题。ans =0需要改进为：假定生产厂D，将产销不平衡问题转化为产销平衡问题在进行求解。 销地 产地产量（万吨）A16161322171750B14141319151560C19192023MM50DM0M0M050销量302070301050 Cost=16 16 13 22 17 17;14 14 13 19 15 15;19 19 20 23 1000 1000;1000 0 1000 0 1000 0； Sc=50 60 50 50； De=30 20 70 30 10 50； LP_transport(Cost,Sc,De)X = 0 0 50 0 0 0 0 0 20 0 10 30 30 20 0 0 0 0 0 0 0 30 0 20ans = 2460最优解矩阵为：X。此时最