找回密码
 注册
关于网站域名变更的通知
查看: 692|回复: 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; o0 R( l6 ], k2 a
% s.t.3 Q. ~) R1 W2 U2 U# }& ~6 T! `
% 0<=xi<=99(i=1,2,...,5)
; q$ \5 Y- u9 ]+ e% B5 V# c8 [% x1+x2+x3+x4+x5<=4004 b6 Z+ h8 w9 x; Z6 G% H7 D: C
% x1+2*x2+2*x3+x4+6*x5<=800' @. ]3 j$ k* E8 [6 R& P! z' z' E& i) ]/ C
% 2*x1+x2+6*x3<=800$ i- x' Q+ `; d
% x3+x4+5*x5<=200
% d0 V- `! n% S7 G& a, g, Cx=intvar(1,5);
+ C" u8 X' T+ bf=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
8 w; I0 q# R3 x: QF=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...& ~3 I" f1 \1 ?; Y0 ~
    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
- B! R: F& n  E1 d& I0 L) AsolvESDp(F,-f)
7 I9 ]3 P6 H7 D; ~double(f): ^" U5 E  Z8 h$ D! _; D
double(x)
4 a4 [6 P. ~: p: J# K6 m. P
3 Y7 z1 N) z7 U5 A5 k( J; r3 J' F3 F/ [$ y0 _

5 M1 e% t) R5 d2 m7 }  K( ?! B' Zans =. s; B2 N2 P! l8 l
9 k3 l9 U7 _' z/ B* r) H1 o5 {
    -5
# A3 I1 m. P. `- y" i& aans =% ?4 U; W! C2 [* f2 M
5 ]% E7 L- c0 F
     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。. K+ P% i* C* v7 [( ]
                                                             正确解释:  53    99    99    99     0                        80199
2 R5 x$ w4 j/ H0 Q

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means
0 Y: ?4 y# N/ U: t, G5 Ithat the branching process is not guaranteed to find a
7 H: @# `7 }' m( Y) Y$ U, Mglobally optimal solution, since the lower bound can be; [# R, U' R9 ~8 M% V
invalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
# j1 O5 [0 F0 L0 D5 Nf=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);
% q; M0 ?" V5 r1 m1 R. mF=[];
  K- E2 I: r# U' q8 r- _8 y( J+ EF=[F x>=0;x<=99];) P4 D& o/ z# d% i8 c) V. ?. X
F=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];
. [4 A5 o! t! _0 w, xF=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];
/ |3 T. T0 H: d7 Y8 d) MF=[F 2*x(1)+x(2)+6*x(3)<=800];
7 H3 w% z# d: K+ D3 C, QF=[F x(3)+x(4)+5*x(5)<=200];7 H; u! S* ?: q" u
ops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题# g- h( }" M9 Z3 u5 c
optimize(F, -f, ops);" b" t) v- l, H5 n% o; q0 g
%solvesdp(F,-f)" R5 n$ T+ t2 s4 G4 L0 `( [
double(f): {- x* o4 z* F  G
double(x)1 B1 g, ~% ]5 o
楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-6-23 12:10 , Processed in 0.093750 second(s), 23 queries , Gzip On.

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

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

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