找回密码
 注册
关于网站域名变更的通知
查看: 768|回复: 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
% _5 a4 a+ \  P3 d1 ?9 |/ |, Z% s.t.4 I5 l( c5 s. m3 i
% 0<=xi<=99(i=1,2,...,5)  o% @6 Q; ~- v0 f; z# f4 f2 n
% x1+x2+x3+x4+x5<=400
! h5 K; k2 o  f% {$ ^$ A! j% x1+2*x2+2*x3+x4+6*x5<=8002 U. w& [1 \8 _3 U, r( i
% 2*x1+x2+6*x3<=800% P  N% v6 q( b1 Y
% x3+x4+5*x5<=200( q! ~( i) f+ w/ M* A; M) D* A! v7 Q7 w4 V
x=intvar(1,5);6 i& X6 s+ X* t6 J& {. ~
f=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
3 R" Z# y# S4 ~% O. hF=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;.... L4 o) _2 Z4 |1 C; k
    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]. P% w) B7 k  l( Z, B
solvESDp(F,-f); P. V% U! y8 f# N7 |9 O
double(f)5 [  W, h; ?. o3 h
double(x)
0 M3 U( b3 @! u
0 K* {8 }1 }. R  Y4 L
# S. r1 m+ v  o; e# L& A6 [/ w2 y4 p0 \
( {9 Y9 e6 R, i9 O( o& N2 kans =; {4 D- [6 Q! @! g8 v) t
& c- F: C+ ?2 V/ v, r! @; X, S
    -5& c( h4 w8 [9 l0 ?$ p
ans =
/ R. W2 N6 M5 n/ W* p9 P* o5 M7 t% ]4 _4 @
     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。+ X6 |( w& b% G. v4 k) e
                                                             正确解释:  53    99    99    99     0                        80199, \; y/ Z. [: W. d8 V

该用户从未签到

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

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means* \% l* \$ h* V
that the branching process is not guaranteed to find a
' o! }2 b0 X4 o" v; X8 B% w) Iglobally optimal solution, since the lower bound can be
; I$ S1 l* P. S" c. w' Oinvalid. Hence, do not trust the bound or the gap...

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
3 P/ J+ k5 ~0 \3 ?% h6 \4 v! kf=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);
$ L7 _( y# {" C" YF=[];9 x8 C! l/ g6 _2 ^# c) J
F=[F x>=0;x<=99];
; @6 _: K: ]' ]6 z% u8 ~+ X4 ZF=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];+ o6 V0 l0 M( J; G$ T8 h% _
F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];( K# a; y8 G; }0 ]4 B+ Y4 S4 ?
F=[F 2*x(1)+x(2)+6*x(3)<=800];
4 ]) `/ K* o2 UF=[F x(3)+x(4)+5*x(5)<=200];
/ A. u% W3 j" A% A( h. D7 Cops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题
& G& m1 U9 S$ o3 o* b5 |* h5 qoptimize(F, -f, ops);9 u& `1 }! n4 a) L: N; @# R
%solvesdp(F,-f)6 y/ n, g0 z4 C  g" [# U! q# W. ?
double(f)  J. R, |( P. H8 }; ]
double(x)
/ n9 J' ?  i9 X/ ?" a8 i楼主试试这个程序,可以得出正确的结果

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-5 02:45 , Processed in 0.140625 second(s), 24 queries , Gzip On.

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

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

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