gusucode.com > 数学建模竞赛：露天矿的开采问题，运用lingo与matlab编程，内附代码 > 数学建模竞赛：露天矿的开采问题，运用lingo与matlab编程，内附代码/work4/mmax.m

    f=[-1 -1 -1 -1 -1 -1 -1 -1]        %目标函数系数矩阵
a=[0.0826    0.0478    0.0528    0.1199   -0.0257    0.1216    0.1220 -0.0015
   -0.0547    0.1275   -0.0001    0.0332    0.0867   -0.1127   -0.0186 0.0985
   -0.0694    0.0620   -0.1565   -0.0099    0.1247    0.0024   -0.0028 -0.2012
   -0.0346   -0.1028    0.2050   -0.0209   -0.0120    0.0057    0.1452 0.0763
   0.0003    0.2428   -0.0647   -0.0412   -0.0655    0.0700   -0.0039 -0.0092
   0.2376   -0.0607   -0.0781    0.0929    0.0466   -0.0003    0.1664 0.0004]
%a为拟合线性方程关系的系数矩阵
aa=[a;-a]   %限制条件的系数矩阵
ri=[0.13 0.18 0.09 0.11 0.15 0.14]  %安全裕度
y1=[165 150 160 155 132 162]-[110.4775  131.3521 -108.9928  77.6116   133.1334   120.8481]
y2=[165 150 160 155 132 162]+[110.4775  131.3521 -108.9928  77.6116   133.1334   120.8481]
%y1,y2为限制条件的常数项
b=[y1,y2]
xi=[120	73	180	80	125	125	81.1 90]
%xi为0方案的各机组出力
vi=[2.2 1 3.2 1.3 1.8 2 1.4 1.8]
%vi爬坡率
vlb=xi-15*vi%限制条件下限
vub=xi+15*vi%限制条件上限
z=linprog(f,aa,b,[],[],vlb,vub)
%求最优线性规划
sum(z)%求出力总和