找回密码
 注册
关于网站域名变更的通知
查看: 802|回复: 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
7 q! K) o3 M" M3 o% s.t." q# v# w1 g1 H
% 0<=xi<=99(i=1,2,...,5)0 s4 x- p! Y( X" r% c" G
% x1+x2+x3+x4+x5<=4007 m) E. K) V" d: l) b
% x1+2*x2+2*x3+x4+6*x5<=8007 L  Z' }+ I" P/ t; U1 P' r
% 2*x1+x2+6*x3<=800
, i5 F+ X0 p# Q9 Q/ F6 n3 Y+ ]7 A% x3+x4+5*x5<=200
- c# G9 ^, `$ f$ Dx=intvar(1,5);" v7 a- u# p& g& }! q+ |
f=[1 1 3 4 2]*(x'.^2)-[8 2 3 1 2]*x';
. n4 z3 i8 a7 i# x$ yF=[x>=0;x<=99;[1 1 1 1 1]*x'<=400;[1 2 2 1 6]*x'<=800;...
6 h; |7 O5 o/ Z- z    [2 1 6 0 0]*x'<=800;[0 0 1 1 5]*x'<=200]3 ], Z+ s* r0 h# |
solvESDp(F,-f)
  ~3 S7 ?# }- o- ^double(f)! Y, q. `' W! |$ _% E* P1 A
double(x)2 s! o% n& }, x( Q

) i2 \6 [! x4 b) y" o
" o* [6 N: A  {  \$ W" o7 n& ~  T+ D0 ?# I
ans =: `& T& z+ O! X$ T
! Z/ B5 Z3 q, ?5 E1 a
    -5
+ I2 c2 N; J# s5 Z! s0 e$ ^ans =
1 y8 L+ W* x: x/ H/ t! Y0 G5 @: W% b1 j$ o. d& ^- y* n* l4 q
     1     1     1     1     1                             求助, 这个结果是错的!!但是我真是找不出我错哪里了。
9 l) S0 Z! i1 j/ l& q# `6 E                                                             正确解释:  53    99    99    99     0                        801993 ?' A- ?* x$ |8 _. T& o

该用户从未签到

5#
发表于 2020-7-30 13:39 | 只看该作者
来学习一下

该用户从未签到

4#
发表于 2020-7-30 09:53 | 只看该作者
x=intvar(1,5);
$ d' C3 m+ P! u( T9 }" df=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);- L, D- \$ E- e! B
F=[];
$ f  K$ _9 r$ W3 rF=[F x>=0;x<=99];
5 m9 |1 u4 y. PF=[F x(1)+x(2)+x(3)+x(4)+x(5)<=400];% `8 I1 B7 n2 v
F=[F x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)<=800];
+ h9 Y' g4 L0 r5 i7 ZF=[F 2*x(1)+x(2)+6*x(3)<=800];
2 a% d* c4 p( _& z0 F, SF=[F x(3)+x(4)+5*x(5)<=200];
* F5 @: }. `: rops = sdpsettings('solver', '+Gurobi'); %利用cplex求解问题. q; ~/ c8 [; H, h6 Z: O7 A& {
optimize(F, -f, ops);5 a, {# T2 q( c$ |1 W9 k
%solvesdp(F,-f)+ [# D5 t  j8 A4 e0 F  e% }4 t
double(f)# {0 N1 P4 J7 t$ {
double(x)
. o; ^, X# G( I- ^( ~楼主试试这个程序,可以得出正确的结果

该用户从未签到

3#
发表于 2020-7-29 13:58 | 只看该作者
Warning : The continuous relaxation may be nonconvex. This means$ d! G7 ?# @% j$ o+ t2 n
that the branching process is not guaranteed to find a
' I* m, }" u; Fglobally optimal solution, since the lower bound can be
5 X, G4 |2 x+ `; x+ g4 ]invalid. Hence, do not trust the bound or the gap...

该用户从未签到

2#
发表于 2020-7-27 10:28 | 只看该作者
帮你顶一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

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

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

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

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