找回密码
 注册
关于网站域名变更的通知
查看: 699|回复: 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
2 \8 Q% l7 ?* O) Z; c% s.t.
" O. ^$ j! P; S8 `5 q" D/ N& m% 0<=xi<=99(i=1,2,...,5)" ?# Z7 r! l( q/ V: _
% x1+x2+x3+x4+x5<=400
3 I- q# g/ U* ?8 F6 M) Y% x1+2*x2+2*x3+x4+6*x5<=800
: K- y8 r1 i7 ?! M" v" V( Y  a% 2*x1+x2+6*x3<=8003 d3 ^0 g2 @0 f
% x3+x4+5*x5<=2001 c( o8 v3 a; D* a# R
x=intvar(1,5);4 O$ \  N" z* r/ m
f=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
" V/ ]; ?3 U! J9 t1 e( y- XF=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...
3 _, u7 s; r" ]    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
+ A% {  D. r' QsolvESDp(F,-f)- T  }! s/ Y/ C$ ?
double(f)2 j: P1 I4 r+ g% D& {# ~$ s+ P
double(x)  ~6 M9 C- U- k( D

3 I1 E" O8 X( H& i" _; X4 Q5 @6 }  Y
4 t/ U) |8 u/ S3 ^% ?6 `
3 f* a) L# }4 _8 {+ y# I- X+ P  a; Ians =
1 |0 E5 C8 Q; y, G, H
1 l) m8 u5 w- S3 q9 M& B    -5
, w9 G8 S6 g$ E" E6 Cans =0 p8 m! S8 {0 L7 o! A# U

/ P0 t/ p2 _# h$ z, n     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。
+ U- d0 f; i/ A                                                             正确解释:  53    99    99    99     0                        80199
, L4 S2 P& J* V! i+ E$ k

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means
9 C( \$ B6 w- {. ^. n8 dthat the branching process is not guaranteed to find a
7 H4 @( F9 t; v. O1 I# pglobally optimal solution, since the lower bound can be
5 X4 V  Q; d, ?4 I' V8 P/ l( l2 D! linvalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
0 s% o6 F. d3 c9 ?: S+ mf=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);
# u( ~1 y7 |( q) C) bF=[];
2 h: B) }8 E1 e2 H! ZF=[F x>=0;x<=99];+ i& R6 I& _3 S
F=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];9 G1 R+ Q9 e3 x( r) T- g$ f
F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];
2 D/ Y' W4 @2 c' s' gF=[F 2*x(1)+x(2)+6*x(3)<=800];& @; {6 v6 s- W% E
F=[F x(3)+x(4)+5*x(5)<=200];
7 y8 @. o1 k3 @ops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题
/ f& X: K+ p! n  m" c; M2 poptimize(F, -f, ops);
- H3 T  B4 f4 D  E9 C# \$ }%solvesdp(F,-f)
' _% }( |; v+ z& b1 F3 S9 e; ?double(f)+ }- J0 u8 g* u- Q) @" ^  R8 G
double(x)
3 e( x5 ^7 [! ^2 j% Y" P# s楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-7-27 08:02 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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