找回密码
 注册
关于网站域名变更的通知
查看: 698|回复: 4
打印 上一主题 下一主题

yalmip的问题

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 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
9 O: S3 D5 Q' F6 u1 _1 n5 M% s.t.9 h9 z) F3 a7 Y; ?' z  N' }% M
% 0<=xi<=99(i=1,2,...,5)! W' y& Y: i. x1 Z1 R
% x1+x2+x3+x4+x5<=4002 W0 A1 J6 v7 e& g8 o( o
% x1+2*x2+2*x3+x4+6*x5<=800
( z6 S6 {' f1 y% 2*x1+x2+6*x3<=800/ O9 Z6 X7 K* Z- S& q
% x3+x4+5*x5<=2002 H  [6 r; \+ O# H! a% f
x=intvar(1,5);4 I+ g5 y1 u$ t5 S' Q+ J0 ~
f=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
3 n4 ]+ C- b' `2 ~  AF=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...
. \' C, V! P, l+ B/ Q! |# i1 y2 k    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
" N! {; w% D1 y3 B4 x) ]9 f8 G. m+ FsolvESDp(F,-f)+ U6 T5 H; c- j4 O, i. p
double(f)
! d0 q0 m8 B( M  u. y% ?double(x)/ o* l8 k( o- H2 H  I0 I3 j& z
0 c7 x; [* q% F  f( f) Z& w
: q7 M- K1 r0 t' t' e
" I6 v3 g/ n7 E; _* w- f6 E/ T; ~  Q
ans =
1 P- ~! v: {9 u, A, g# i$ v
8 X$ C, g: \& p0 D/ g. d% ~& I4 t    -5" k/ `" o# @5 a. x- t
ans =
" o2 n6 q& Q7 r7 W+ X& Z& k- P+ n3 `  s+ D% f. i8 k9 g+ I
     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。
% k7 l5 |  h7 {3 ?) @1 y                                                             正确解释:  53    99    99    99     0                        80199
! _: g: B: [. j  [2 S1 b

该用户从未签到

2#
发表于 2020-7-27 10:28 | 只看该作者
帮你顶一下

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means: \( P( }5 G: C: L
that the branching process is not guaranteed to find a- w  b) L/ x" U" d
globally optimal solution, since the lower bound can be
1 R( }4 i% e. ]" ^+ q$ Finvalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
2 `7 }- W; e$ W: K# a8 Rf=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);+ K' P+ y/ h  T9 i9 P) K
F=[];( M0 O+ \' e' S- h$ w
F=[F x>=0;x<=99];% I" Q7 P& s/ k- r( t
F=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];
% I; X6 ]& M' S: q! W6 dF=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];& b# p4 M. i+ \
F=[F 2*x(1)+x(2)+6*x(3)<=800];8 }' u& k9 Q5 C/ {# F6 A* _2 J6 i
F=[F x(3)+x(4)+5*x(5)<=200];/ e  d5 I! r4 W( r! w
ops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题
! K2 s0 o. J1 @2 i' G, v$ t& W1 r1 I0 boptimize(F, -f, ops);
. a! f5 ^) a& U: ?2 S9 E%solvesdp(F,-f)
3 N8 Q' q' p5 R" ~double(f)
) ]4 C, h0 r4 r4 ]! O) ?double(x)
/ a5 e/ `" I2 [- a, ^3 J$ H楼主试试这个程序,可以得出正确的结果

该用户从未签到

5#
发表于 2020-7-30 13:39 | 只看该作者
来学习一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-7-26 20:01 , Processed in 0.125000 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表