找回密码
 注册
关于网站域名变更的通知
查看: 758|回复: 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. I+ b  V) _" _* T. P" I
% s.t.
% o8 ]+ U1 v$ \% 0<=xi<=99(i=1,2,...,5)( S4 v; H/ _% V5 I
% x1+x2+x3+x4+x5<=400" _  k8 `' j4 R$ d2 N3 i7 P4 U, O/ x
% x1+2*x2+2*x3+x4+6*x5<=8001 }  C# k8 J# k' d6 H9 ^  c3 @0 n
% 2*x1+x2+6*x3<=800; G3 j: c. ^" j0 ~( t
% x3+x4+5*x5<=200+ A: X9 A* }% e& T3 N" G! l& A& G
x=intvar(1,5);
) t6 z& K5 {" I/ _5 _. g7 }, U% Cf=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';/ Z. O& ]& I0 e0 q3 k5 ~  P8 M0 S4 L" u
F=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...  Q1 ?: u/ X& T$ }4 S1 ^
    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
- j& h$ x4 J! a4 {+ S# r( jsolvESDp(F,-f)3 y) U& X2 o2 ~& \" M
double(f)
. C  t+ u' s$ g8 h4 s: Udouble(x)
- h. i% H) e+ Z' e3 B
0 q& Z5 {! d: H% X9 w$ b! D+ V* U$ |
% ?/ R+ S2 h. c! j2 ~6 h
ans =% Z" q$ A! W4 }

# W) o. _9 @" c. `9 Z    -5
: z4 D+ m3 |0 G, s3 \% nans =- J# r% G/ e9 H1 c* W* o& |; u
; E) z! p- K( v! m( B9 S
     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。7 A9 r# m$ V* ^! {
                                                             正确解释:  53    99    99    99     0                        80199
$ k  ~( F' x& B) o9 q- m0 F* Z; i2 S

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means4 F* |' r2 K" j3 A+ Z5 U$ J
that the branching process is not guaranteed to find a  y) c/ C6 P: O  E6 ?' j& ]/ s5 V
globally optimal solution, since the lower bound can be' \& I5 W  c9 r- i% R3 z
invalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
6 k3 q- i6 ]" L/ _" tf=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);
2 W2 L3 b- Q) F. X( k8 gF=[];7 d5 }" q% `  r0 m
F=[F x>=0;x<=99];0 S3 X$ F; s/ D( j2 Y$ H" [
F=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];7 ]7 F' |8 m! x  F7 q0 n
F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];- v0 I) M3 q6 o$ V7 L! ]* }
F=[F 2*x(1)+x(2)+6*x(3)<=800];3 l. e" u5 V' U9 U! k% v$ i9 p$ h
F=[F x(3)+x(4)+5*x(5)<=200];
- A9 w6 ]/ L. `/ {ops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题( |9 n7 N% t( r
optimize(F, -f, ops);
: _) U* m1 R( V, y; K5 K/ K2 C%solvesdp(F,-f)' f! c8 d( A+ V+ |' D8 n5 {; }  L
double(f)
! C! L# ?) D0 @3 ddouble(x)# m* M; V1 ?3 r) \; t7 l
楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-10-8 01:44 , Processed in 0.109375 second(s), 23 queries , Gzip On.

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

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

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