| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1122 人关注过本帖, 1 人收藏
标题:使用yalmip工具箱求解半定规划问题时显示凸检测失败
取消只看楼主 加入收藏
叁生石0212
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2013-9-12
收藏(1)
 问题点数:0 回复次数:0 
使用yalmip工具箱求解半定规划问题时显示凸检测失败
程序代码如下:
%%定义变量
r=sdpvar(1,1);
t=sdpvar(1,1);
x=sdpvar(6,1);
Q=sdpvar(6,6);
q=sdpvar(6,1);
eit=sdpvar(1,1);
y=sdpvar(1,1);
w=sdpvar(6,1);
sigma0=[0.5 0.6 0.7 0.8 0.9 1;0.6 0.5 0.6 0.7 0.8 0.9;0.7 0.6 0.5 0.6 0.7 0.8;
       0.8 0.7 0.6 0.5 0.6 0.7;0.9 0.8 0.7 0.6 0.5 0.6;1 0.9 0.8 0.7 0.6 0.5];
mju0=[2.8;3;3.2;4;5;5.5];
a=[0 0 0 0 0 0]';
b=[2 1.75 1 3.25 3 3]';
c=[0.00375 0.0175 0.0625 0.00834 0.025 0.025]';
p=[1 1 1 1 1 1]';
%%约束条件
F=[Q>=0];
F=F+[trace((2.3.*sigma0+mju0*mju0')*Q)+mju0'*q+0.2^0.5.*norm(sqrtm(sigma0)*(q+2.*Q*mju0))];
F=F+[[Q q/2+x/2;q'/2+x'/2 r-eit]>=0];
F=F+[eit>=p'*a+b'*x+p'*x+y];
F=F+[y>=w'*w];
F=F+[w>=sqrt(c).*x];
F=F+[w<=sqrt(c).*x];
F=F+[p'*x<=150];
F=F+[[50 20 15 10 10 12]'<=x];
F=F+[x<=[200 80 50 35 30 40]'];
%%目标函数
f=r+t;
%%调用解算器
solvesdp(F,f)

运行后显示为:
ans =
    solvertime: 0
       problem: 14
          info: 'Convexity check failed (, exact model not available.)'
    yalmiptime: 0.0460

如果在上述程序下加double(f),显示为:
ans =
   NaN
 
搜索更多相关主题的帖子: 工具箱 检测 
2013-09-12 20:38
快速回复:使用yalmip工具箱求解半定规划问题时显示凸检测失败
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.032674 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved