yalmip的问题

shuddkk · 发表于 2020-7-27 09:49

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

% Max z=x1^2+x2^2+3*x3^2+4*x4^2+2*x5^2-8*x1-2*x2-3*x3-x4-2*x5
% s.t.
% 0<=xi<=99(i=1,2,...,5)
% x1+x2+x3+x4+x5<=400
% x1+2*x2+2*x3+x4+6*x5<=800
% 2*x1+x2+6*x3<=800
% x3+x4+5*x5<=200
x=intvar(1,5);
f=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
F=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...
[2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
solvESDp(F,-f)
double(f)
double(x)

ans =

-5
ans =

1 1 1 1 1 求助，这个结果是错的！！但是我真是找不出我错哪里了。
正确解释： 53 99 99 99 0 80199

大小的小 · 发表于 2020-7-27 10:28

帮你顶一下

小小鲁班 · 发表于 2020-7-29 13:58

Warning : The continuous relaxation may be nonconvex. This means
that the branching process is not guaranteed to find a
globally optimal solution, since the lower bound can be
invalid. Hence, do not trust the bound or the gap...

lupkpu9 · 发表于 2020-7-30 09:53

x=intvar(1,5);
f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)^2-8*x(1)-2*x(2)-3*x(3)-x(4)-2*x(5);
F=[];
F=[F x>=0;x<=99];
F=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];
F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];
F=[F 2*x(1)+x(2)+6*x(3)<=800];
F=[F x(3)+x(4)+5*x(5)<=200];
ops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题
optimize(F, -f, ops);
%solvesdp(F,-f)
double(f)
double(x)
楼主试试这个程序，可以得出正确的结果

nkkopd · 发表于 2020-7-30 13:39

来学习一下

帐号		自动登录	找回密码
密码			注册

yalmip的问题

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1