找回密码
 注册
关于网站域名变更的通知
查看: 709|回复: 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
" B% _5 P7 H. H+ c% s.t.4 Q) M" s# ?$ b9 ^$ M
% 0<=xi<=99(i=1,2,...,5)7 a/ t$ A0 t2 }4 h/ v( P" E
% x1+x2+x3+x4+x5<=400
7 Z+ |( E- D# I/ K" U9 @" O% x1+2*x2+2*x3+x4+6*x5<=800
: x1 h% K6 u5 D( i% 2*x1+x2+6*x3<=800
5 p* I7 B. x3 i0 h0 W4 ~3 ?1 z% x3+x4+5*x5<=200' u# R" R, @6 N* x, I1 E9 k
x=intvar(1,5);, z+ C! N5 p: Q2 S* D5 [
f=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';7 P  w  p" \& ~% x0 M5 l
F=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...
$ t& n. ]% _7 j5 }    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
% w% o$ R/ S. ]$ ksolvESDp(F,-f)9 q: A+ U5 ]1 Z2 x# h" \, T6 ^
double(f)
$ v1 t- _' m3 \* _double(x)
! w+ A' G1 _) d! n. s9 G3 p" q* @: M) |5 v+ I9 E, O

+ V2 p* `! l, F5 [$ |% O7 h
  i' Y. o6 D. T- ]( ?) Vans =
" O8 Z& z+ T) M+ S2 ?# @, @6 C. h" K& b7 Y! A
    -5
/ d6 {# C1 {/ s+ J) x2 A0 Uans =% g, Z% Z" I- f2 {) m: i( r1 g

  D) I0 o, h  B6 s7 J5 Y. h     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。: h. j2 Y; l/ U+ Q
                                                             正确解释:  53    99    99    99     0                        80199
% D5 e0 h! G! U7 N2 p3 o

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means
0 @2 }! \  V) N) G6 o; @" [that the branching process is not guaranteed to find a! P! T8 Q% i  N, A8 V6 {7 C2 C
globally optimal solution, since the lower bound can be
+ P6 [( T. g- g% @; ^8 x# r- H# Ginvalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);+ B  ]0 u5 i0 N/ I4 X: C
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);
; c: N1 `. S( @' B8 hF=[];
$ K) R6 c' ]; N3 N, E: DF=[F x>=0;x<=99];
+ Y8 T( c9 k/ S/ zF=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];. O* J  X& o. c' Y, z* U3 l9 w
F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];8 u5 I7 i! \9 n2 `1 |9 J. Q, D
F=[F 2*x(1)+x(2)+6*x(3)<=800];
& ?  r! S4 K2 g; m5 }! k, RF=[F x(3)+x(4)+5*x(5)<=200];
- b1 E2 M9 |/ q' L" ?& E( lops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题
/ T8 H7 S& ~- p: zoptimize(F, -f, ops);. F& B$ \# n1 A! {, G; [9 b
%solvesdp(F,-f)
' ^- `6 L, t9 D4 s1 J  x3 Fdouble(f)1 n, K! u/ {( @' h- O% [8 |
double(x)
1 e! k& A0 |# [# v+ @; D楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-8-19 11:51 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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