s154.m

来自「该文件主要包括了matlab的一些主要应用例程!对学习matlab很有帮助」· M 代码 · 共 27 行

function [x,I]=s154(m,n,a,b,c)
% m为生产地点数目，正整数
% n为销售地点数目，正整数
% a为m维产量行向量
% b为n维销售量行向量
% c为m*n阶单位运价矩阵
% x为m*n阶最优调运方案矩阵
% I为最小总运价
B=[a,b];
d=reshape(c',1,m*n);
F=[zeros(m,n);eye(n,n)];
A=F;
for i=1:m-1
   for j=(i-1)*n+1:i*n
      A(i,j)=A(i,j)+1;
   end
   A=[A,F];
end
for j=(m-1)*n+1:m*n
   A(m,j)=A(m,j)+1;
end
v1=zeros(m*n,1);
options(13)=m+n;
y=linprog(d,[],[],A,B,v1);
I=d*y;
x=reshape(y,n,m);
x=x';