线性规划（matlab篇）

线性规划

1、所有线性规划问题化成matlab中规定的标准形式：

f, x, b, beq, lb, ub为列向量，其中f为价值向量，b为资源向量，A，Aeq为矩阵。

2、linprog函数

上述对应Matlab中的求解命令为：

[x,fval] = linprog(f,A,b)
[x,fval] = linprog(f,A,b,Aeq,beq)
[x,fval] = linprog(f,A,b,Aeq,beq,lb,ub)

x返回决策向量的取值；fval返回目标函数的最优值；f为价值向量；

A、b对应不等式约束

Aeq和beq对应线性等式约束

lb和ub分别对应决策向量的下界和上界

3、Easy example:

$$\begin{gathered} maxz=2x_1+3x_2-5x_3, \\ \{\begin{array}{l}{x}_1+x_2+x_3=7,\\ 2x_1-5x_2+x_3\ge 10,\\ x_1+3x_2+x_3\le 12,\\ x_1,x_2,x_3\ge 0\end{array} \end{gathered}$$

①化为matlab标准型：
$$\begin{gathered} minw=-2x_1-3x_2+5x_3, \\ \left(\begin{array}{ccc}-2& 5& -1\\ 1& 3& 1\end{array}\right)\left(\begin{array}{c}{x}_1\\ x_2\\ x_3\end{array}\right)\le \left(\begin{array}{c}-10\\ 12\end{array}\right), \\ \left[1,1,1\right]\cdot {\left[x_1,x_2,x_3\right]}^T=7, \\ {\left[x_1,x_2,x_3\right]}^T\ge [0,0,0{]}^T \end{gathered}$$

>> f=[-2;-3;5];			% 价值向量
>> a=[-2,5,-1;1,3,1];	% a、b对应不等式的左边和右边
>> b=[-10;12];
>> aeq=[1,1,1];			% aeq和beq对应等式的左边和右边
>> beq=7;
>> [x,y]=linprog(f,a,b,aeq,beq,zeros(3,1));

Optimal solution found.

>> x					% x=[x1;x2;x3]最优解

x =

    6.4286
    0.5714
         0

>> y=-y					% 最优值（该题所求为max，但linprog函数求得的是min，故只需加个-）

y =

   14.5714

4、example2:

某架货机有三个货舱：前仓、中仓、后仓。三个货舱所能装载的货物的最大重量和体积都有限制，如表2所示。并且，为了保持飞机的平衡，三个货舱中实际装载货物的重量必须与其最大容许重量成比例。

现有四类货物用该货机进行装运：

​ 货物规格及利润表

假设：

（1）每种货物可以无限细分

（2）每种货物可以分布在一个或者多个货舱内

（3）不同的货物可以放在同一个货舱内，并且可以保证不留空隙

问：应如何装运，使货机飞行利润最大？

解：

用 i = 1,2,3,4表示货物1、2、3、4；用 j =1,2,3表示前舱、中仓、后仓。
$$\begin{gathered} 设x_{ij}表示第i种货物在第j个货舱内的重量,w_j和v_j表示可以运输的第j个舱的重量限制和提及限制,a_i,b_i,c_i表示 \\ 第i种货物的重量，单位重量所占的空间体积,单位货物的利润。 \end{gathered}$$
（1）目标函数：
$$z=c_1\sum _{j=1}^3{x}_{1j}+c_2\sum _{j=1}^3{x}_{2j}+c_3\sum _{j=1}^3{x}_{3j}+c_4\sum _{j=1}^3{x}_{4j}=\sum _{i=1}^4{c}_i\sum _{j=1}^3{x}_{ij}$$
（2）约束条件：

① 四种货物的重量约束：
$$\sum _{j=1}^3{x}_{ij}\le a_i,i=1,2,3,4$$
② 三个货舱的重量限制：
$$\sum _{i=1}^4{x}_{ij}\le w_j,j=1,2,3$$
③ 三个货舱的体积限制：
$$\sum _{i=1}^4{b}_i{x}_{ij}\le v_j,j=1,2,3$$
④ 三个货舱的平衡限制:
$$\frac{{\sum }_{i=1}^4{x}_{i1}}{10}=\frac{{\sum }_{i=1}^4{x}_{i2}}{16}=\frac{{\sum }_{i=1}^4{x}_{i3}}{8}$$
综上，建立线性规划模型：
$$\begin{gathered} maxz=\sum _{i=1}^4{c}_i\sum _{j=1}^3{x}_{ij} \\ \{\begin{array}{l}{\sum }_{j=1}^3{x}_{ij}\le a_i,i=1,2,3,4\\ {\sum }_{i=1}^4{x}_{ij}\le w_j,j=1,2,3\\ {\sum }_{i=1}^4{b}_i{x}_{ij}\le v_j,j=1,2,3\\ \frac{{\sum }_{i=1}^4{x}_{i1}}{10}=\frac{{\sum }_{i=1}^4{x}_{i2}}{16}=\frac{{\sum }_{i=1}^4{x}_{i3}}{8}\end{array} \end{gathered}$$

% 目标函数
c=[3100;3800;3500;2850];     % 由目标函数可知有12个变量
c=c*ones(1,3);               % 使c变成4行3列，对应目标函数未知量前面的系数
c=c(:);                    % c每列合并成一个长的列向量
% 约束条件1,需要构造如下数组
a1=zeros(4,12);
for i=1:4                             
    a1(i,i:4:12)=1;                    
end
b1=[18 15 23 12]';
% 约束条件2
a2=zeros(3,12);
for i=1:3
    a2(i,4*i-3:4*i)=1;          
end
b2=[10 16 8]';
% 约束条件3
bb=[480;650;580;390];
a3=zeros(3,12);
for i=1:3
    a3(i,4*i-3:4*i)=bb;           
end
b3=[6800 8700 5300]';
a=[a1;a2;a3];
b=[b1;b2;b3];
%----------------------------
% 等式条件
aeq=zeros(2,12);             
aeq(1,1:4)=1/10;
aeq(1,5:8)=-1/16;
aeq(2,5:8)=1/16;
aeq(2,9:12)=-1/8;
beq=[0;0];
lb=zeros(12,1);
[x,y]=linprog(-c,a,b,aeq,beq,lb,[])
x=reshape(x,[4,3])         % 化为3*4的数组
x=sum(x')                  % 舱1,2,3分别放入货物的总重量
y=-y;                  

y =

  -1.2152e+05


x =

         0         0         0
    7.0000         0    8.0000
    3.0000   12.9474         0
         0    3.0526         0


x =

         0   15.0000   15.9474    3.0526

答： 求得的四种货物的吨数分别为 0， 15 ，15.9474 ，3.0526

​ 总利润为 1.2152×10^5元。