如果你的是matlab5,有lp函数,那么修改如下:
c=[1;1;1;1;1];%用分号
a=[-1,-1,-1,-1,-1];%都要化为小于
b=[-19250];%-13475>-19250,留小的
vlb=[0;0;0;0;0];
vub=[10000;5000;2500;1250;500];
[x,lam]=lp(c,a,b,[],[],vlb,vub)%无等式约束用[]表示出来。
 
如果不是matlab5,没有lp函数,就改用linprog
c=[1;1;1;1;1];%用分号
a=[-1,-1,-1,-1,-1];%都要化为小于
b=[-19250];%-13475>-19250,留小的
vlb=[0;0;0;0;0];
vub=[10000;5000;2500;1250;500];
[x,lam]= linprog(c,a,b,[],[],vlb,vub)
 
计算结果:
 
x =
  1.0e+004 *
    1.0000
    0.5000
    0.2500
    0.1250
    0.0500
lam =
       19250