请高手帮忙用matlab解决线性规划问题

来源:百度知道 编辑:UC知道 时间:2024/09/24 04:19:20
初学线性规划问题,请大家帮帮忙。
用Matlab 求解下列线性规划问题:
min z=|x1|+2|x2|+3|x3|+4|x4|;
s.t. x1-x2-x3+x4=0;
x1-x2+x3-3*x4=1;
x1-x2-2*x3+3*x4=-0.5.
一本书上提供一种转换方法:
对任意的Xi,存在Ui,Vi>0满足
Xi=Ui-Vi,|Xi|=Ui+Vi
Ui=(Xi+|Xi|)/2,Vi=(|Xi|-Xi)/2;
记U=[U1,U2,…Un]^T,V=[V1,V2,……,Vn]^T
原问题可转化为:
min(求和i=1……n,)(Ui+Vi)
s.t.
A(u-v)<=b;
u,v>=0.

不小心的告诉你Matlab的线性规划函数linprog无法求解,非线性的目标函数,因为你的带有绝对值符号,故linprog是无能为力的

但是幸运的是Matlab提供了fmincon

%by dynamic
%see also http://www.matlabsky.com
%2009.2.16
%
%目标函数
objfun=@(x)abs(x(1))+2*abs(x(2))+3*abs(x(3))+4*abs(x(4))
%等式约束
Aeq=[1 -1 -1 1
1 -1 1 -3
1 -1 -2 3];
beq=[0 1 -0.5]';
x0=[0 0 0 0];%给一个初值 很关键 和很重要哦
x=fmincon(objfun,x0,[],[],Aeq,beq)

恩你的想法的确高,佩服,下面是根据的需要改进的程序

min z=[1 2 3 4 1 2 3 4]*[u1 u2 u3 u4 v1 v2 v3 v4]'

s.t.

A=[1 -1 -1 1 1 -1 -1 1
1 -1 1 -3 1 -1 1 -3
1 -1 -2 3 1 -1 -2 3]

x=[u1 u2 u3 u4 v1 v2 v3 v4]'

b=[0 1 -0.5]'

Ax=b
x>=0

%目标函数
f=[1 2 3 4 1 2 3 4];
%等式约束
Aeq=[1 -1 -1 1 1 -1 -1 1
1 -1 1 -3 1 -1 1 -3
1 -1 -2 3 1 -1 -2 3];
beq=[0 1 -0.5]';
%边界条件<