数学建模算法与应用习题 1-4 解析 MATLAB 换一下思路做题

数学建模算法与应用 习题1-4

我这里解答方案和答案上有一点点不一样，我是设每个货物的总量为1，这样每个货物在三个仓加起来就是小于等于1。

MatLab代码：

clc;clear;
zero = [0,0,0,0];
%% 这里是利润的判断 一共是 12个数 3个仓*4个货物
c = [3100*18,3800*15,3500*23,2850*12];
c = [c,c,c];
%% 这里是判断每个仓所有的重量加起来不会超过限定量
A1 = [18,15,23,12];
A1 = [A1,zero,zero;zero,A1,zero;zero,zero,A1];
%% 这里是判断每个仓所有的体积加起来不会超过限定量
A2 = [480*18,650*15,580*23,390*12];
A2 = [A2,zero,zero;zero,A2,zero;zero,zero,A2];
%% 这里是判断每个仓的所有的值加起来不会超过限定
% 因为我设定的是分数 所以加起来不会超过1
A3 = zeros(4,12);
for i=1:4
   A3(i,i:4:12) = 1; 
end
%% 这里合并一下
A = [A1;A2;A3];
b = [10,16,8,6800,8700,5300,1,1,1,1]';
%% 这里是判断货物重量和最大限定成比例
A1 = [18,15,23,12];
aeq = [A1,A1,zero;zero,A1,A1];
aeq(1,1:4) = aeq(1,1:4)/10;
aeq(1,5:8) = aeq(1,5:8)/(-16);
aeq(2,5:8) = aeq(2,5:8)/16;
aeq(2,9:12) = aeq(2,9:12)/(-8);
beq = zeros(2,1);
%%
options = optimoptions('linprog','Algorithm','dual-simplex');
[x,y] = linprog(-c,A,b,aeq,beq,zeros(12,1),ones(12,1),options);
%disp(x);
for i = 1:4
    fprintf('货物%d = %.4ft\n',i,(x(i)+x(i+4)+x(i+8))*A1(i));
end
fprintf('最优利润是%.4f元\n',-y);
%disp(-y);

输出：

Optimal solution found.

货物1 = 0.0000t
货物2 = 15.0000t
货物3 = 15.9474t
货物4 = 3.0526t
最优利润是121515.7895元