找回密码
 注册
关于网站域名变更的通知
查看: 799|回复: 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
3 V& ?( r% G( ^" }  x) K% s.t.
8 |; E! X; N6 J- m% 0<=xi<=99(i=1,2,...,5); T: C5 w/ D9 h
% x1+x2+x3+x4+x5<=400, d; l) `* @( ]  G7 o8 Q  h
% x1+2*x2+2*x3+x4+6*x5<=8002 I( x/ K, ]$ {1 H! c
% 2*x1+x2+6*x3<=800" \& h) X; \& `' G$ z0 {* R9 r# b
% x3+x4+5*x5<=2006 j7 k+ s/ i/ ^9 s9 h0 t
x=intvar(1,5);
4 ]$ d+ C& K8 X0 P, Z; Tf=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
: q4 c! z: e; Q" C1 ~F=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...; ~! w; u' y6 k1 V, m
    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
8 E! p/ P% }) x8 w4 y! x# ]8 M- ]solvESDp(F,-f). w, R7 j% e' Q  g! K5 M
double(f)
: x7 u$ @, _1 c8 qdouble(x)
- a$ y/ m! T4 Y& c( X. D! ~9 {
5 }5 h  A* d7 m% @' g, X& ?( v- J& u9 b: P4 f7 z  g% b1 S' J
$ k6 [5 w. ~" t& \) d
ans =
, a* |8 q; V' I" U4 m+ H; K1 a8 F6 ^- {, g) D
    -5% a: J- X; A' |( h8 H
ans =
) t7 B, e$ [) a9 E/ ~: s2 j+ _/ i
2 M& T# n9 i; Z$ B$ \     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。
7 |5 v6 J3 L( w* Y                                                             正确解释:  53    99    99    99     0                        80199
1 ]3 ?  X* m0 L# T# x

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means- s! a8 G9 o0 r8 r! L; l
that the branching process is not guaranteed to find a
: r' l# S1 y- {& W6 J. m. zglobally optimal solution, since the lower bound can be
& {" r' k' [# q4 o; O9 O/ y# J. ^invalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
5 P0 }' H) h* g0 S  \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);
  [' }! ~! v  Q* r! N2 SF=[];! j+ E5 \, }2 \# q
F=[F x>=0;x<=99];2 ?/ G! M3 g' L* H
F=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];1 u. x0 s( P0 D% {3 U
F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];
) }" e& X3 {% O3 G" Y/ iF=[F 2*x(1)+x(2)+6*x(3)<=800];  |6 y. J9 l/ M+ E) e2 s# ?5 `3 w
F=[F x(3)+x(4)+5*x(5)<=200];
9 N- O1 ^2 P) o- A9 Tops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题) W4 I$ ^# h6 _5 }+ `+ t
optimize(F, -f, ops);
- Q5 e( F8 U0 u7 v  H%solvesdp(F,-f). u6 ^$ [! ?- ]) o
double(f)4 X, ~9 O' W! P5 z) `
double(x)* E2 h- V/ l- @, t) G/ U
楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-24 17:59 , Processed in 0.156250 second(s), 24 queries , Gzip On.

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

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

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