找回密码
 注册
关于网站域名变更的通知
查看: 770|回复: 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
1 K/ E# P0 }& K6 k. i% s.t.
2 c! g! o# p6 k5 g: T9 p0 Y% 0<=xi<=99(i=1,2,...,5)4 R! @% d: O" i# j: _( `
% x1+x2+x3+x4+x5<=400
+ S7 [( _$ k$ c8 v* A+ b& b. l1 q% x1+2*x2+2*x3+x4+6*x5<=800# O) s9 N" \% E& ~8 p) T
% 2*x1+x2+6*x3<=800" ~0 E! W8 i8 K$ n/ S- x
% x3+x4+5*x5<=200) p$ {5 _/ t9 `4 e$ \) h
x=intvar(1,5);
; _( _) T  ~; A7 Y! U5 x* ^. tf=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
3 V2 Z7 j* W* O8 b  JF=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...9 ]! l. h; _9 E) z3 c8 S
    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
2 y" j1 {! ~" IsolvESDp(F,-f)+ N: b" `" B8 Q
double(f)
- c1 o2 i" x9 f: R4 kdouble(x)" e& ^* a, Z" A. S

1 O' n6 r7 Q- z/ w6 _+ Z8 _! G( I3 R* R1 O5 o  Y- j2 ?" w
! h& a7 ?9 e: c6 d
ans =$ T$ e8 i- M: z$ K, Y, O

2 B. C& B! t+ Q# |8 r9 o7 L: t    -5
8 a; F# u' {# i$ m( yans =
' n+ V( H8 ]- _; g/ t( C! z- N0 k: A: d' E% R& b) v  ]
     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。) O$ t/ e9 w5 C: [
                                                             正确解释:  53    99    99    99     0                        801999 i* k* Z0 c4 _

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means
! ~. ^1 t# q$ F4 H7 @, Wthat the branching process is not guaranteed to find a
; B( O0 n; B" c, }globally optimal solution, since the lower bound can be
" k; R, \8 @$ Pinvalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
1 B2 V) |# g  i0 }" |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);
  h' L' M; p) o, BF=[];
: h9 M0 G2 ^% x5 A. N! IF=[F x>=0;x<=99];
5 h, {9 n+ K, O. sF=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];
4 o! ]: j& _7 Y3 kF=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];- w- O2 Q- I7 L: s( e2 A7 o
F=[F 2*x(1)+x(2)+6*x(3)<=800];
9 x% y  ?4 T1 s# e4 J! `F=[F x(3)+x(4)+5*x(5)<=200];
! @6 x" a/ x; U; lops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题; v. q( b& u; P) n5 c; @" e
optimize(F, -f, ops);
- K4 b. t' F1 h- B%solvesdp(F,-f)' U: c1 W5 ?1 d8 H
double(f)
# h6 E' W  P6 ]) j* c1 r% Gdouble(x)% Z* q4 }; G9 J2 _7 v+ C2 T/ r( `8 j0 I( e
楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-5 20:02 , Processed in 0.156250 second(s), 23 queries , Gzip On.

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

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

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