找回密码
 注册
关于网站域名变更的通知
查看: 805|回复: 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
/ o5 z! ~9 q: O; k& h1 F" ^3 ?% s.t.& ?7 ]0 R( \& f
% 0<=xi<=99(i=1,2,...,5)$ z4 t. @  r, y/ X! P
% x1+x2+x3+x4+x5<=4000 C# {& v+ o3 M  G: _4 C7 y7 I
% x1+2*x2+2*x3+x4+6*x5<=800# r  w, Q% c5 ^& e7 i; V  A  N
% 2*x1+x2+6*x3<=800
* q: C5 U+ a( b% x3+x4+5*x5<=200
( l+ r! q' Z% k2 yx=intvar(1,5);: z9 r) l1 A% |7 n; d$ t
f=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
% O) u% E& c5 R- xF=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...  M) X# R% q' N* O% \& S8 w
    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]
& E3 r6 g! c+ w7 _3 k1 O$ I* s/ y+ MsolvESDp(F,-f)/ u5 M  _/ v$ _+ O9 z! o1 p
double(f)
8 I# s8 s/ ]  F( {double(x)
$ |4 X. r! X; c2 r/ ?5 ]! R6 |- t4 I# U+ i# }! h" K0 g2 k  v
$ K7 j0 V) j$ @
$ v7 h  @' j2 V9 M3 _1 w; ?- m
ans =
& J' j: P* B9 L, j3 Z3 c# O( n+ i  S( m+ d$ `+ R
    -5
7 z$ Q  k. P% q9 Bans =- A% A& P& U( x) j1 m& F6 w

$ ?( \# Y" d3 w: d+ W  h# `+ k     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。
0 |0 T+ `. U. E+ L; Q( l                                                             正确解释:  53    99    99    99     0                        80199
/ J8 @1 Z: W7 t) T

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means# H# [0 g, f2 ?$ @! m
that the branching process is not guaranteed to find a5 `( {* q( k/ F" z
globally optimal solution, since the lower bound can be
4 O( |9 |- r1 D3 H+ s/ C( E! q, _0 j  Yinvalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);1 T1 F& w) b0 G( e7 j
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);, D5 ^1 J' _# Y' V4 l7 v* K) b* J( ^
F=[];
4 e  X# y" Q% d9 jF=[F x>=0;x<=99];: p) k) d! H! s
F=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];
9 \# v5 Z1 ~- \F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];- @) ~1 a) G& V+ f7 S
F=[F 2*x(1)+x(2)+6*x(3)<=800];
; r1 K0 e5 T7 l& b/ |F=[F x(3)+x(4)+5*x(5)<=200];
6 p  M6 E0 R2 Z( `6 cops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题
4 {5 Z1 U+ O2 Zoptimize(F, -f, ops);
* q) e) V! u4 i( H%solvesdp(F,-f)
/ O. G- Z8 k3 ^# m7 fdouble(f)
$ M4 E: [$ i- X: j+ f( gdouble(x)6 H1 p9 {1 s" M: n* Z9 T
楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-24 22:42 , Processed in 0.171875 second(s), 23 queries , Gzip On.

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

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

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