| 查看: 4281 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
凝雪小小新虫 (小有名气)
|
[求助]
matlab线性矩阵不等式求解做优化问题
|
|
|
在求解如下凸优化问题中遇到了问题： 1.不知道如何编程得到最优的γ² 2.求解提示为没有可行解 是否是我退出的不等式有问题 或者其他方面有问题 A=[0,1,0,-1;-882,-28.4,0,28.4;0,0,0,1;1696.15,54.62,-657.12,-2680.39;]; B=[0,0,0,-0.039]'; Bw=[0,0,-1,2625.77]'; C1=[-882,-28.4,0,28.4]; C2=[0,0,10,0]; Ea=[0,0,-65.71,262.58] %Eb=0 L=[0,0,0,1]' I=eye(1) gam=4.2; %γ=gam a=0.01; %ρ=a b=0.1; %ε=b c=inv(b); umax=3000 setlmis([]); X=lmivar(1,[4 1]); %定义决策变量 Z=lmivar(2,[1 4]); lmiterm([1 1 1 X],A,1,'s'); lmiterm([1 1 1 Z],B,1,'s'); lmiterm([1 1 2 0],L); lmiterm([1 1 3 X],1,Ea'); %Ea→Ea' lmiterm([1 1 4 0],Bw); lmiterm([1 1 5 X],1,C1'); lmiterm([1 2 2 0],-c); lmiterm([1 3 3 0],-b); lmiterm([1 4 4 0],-gam^2); lmiterm([1 5 5 0],-1); lmiterm([-2 1 1 X],1,1); lmiterm([3 1 1 0],-1); lmiterm([3 1 2 Z],0.95,1); lmiterm([3 2 2 X],-3000,1); lmiterm([4 1 1 0],-1); lmiterm([4 1 2 X],C2,0.95); lmiterm([4 2 2 X],-1,1); lmisys=getlmis; %完成LTI框架的设设置 [tmin,xfeas]=feasp(lmisys); %求解可行解问题 X=dec2mat(lmisys,xfeas,X); %提取解矩阵 把决策变量转化为矩阵形式 Z=dec2mat(lmisys,xfeas,Z); P=inv(X); K=Z*P 运行提示为： Solver for LMI feasibility problems L(x) < R(x) This solver minimizes t subject to L(x) < R(x) + t*I The best value of t should be negative for feasibility Iteration : Best value of t so far 1 2487.161836 2 1661.789005 3 1200.565677 4 542.424422 5 311.999933 6 311.999933 7 279.917289 8 279.917289 *** new lower bound: 36.447556 Result: best value of t: 279.917289 f-radius saturation: 0.000% of R = 1.00e+009 These LMI constraints were found infeasible K = 2.0861 0.1139 -16.8718 3.5340 麻烦各位指出其中有错误的地方 不胜感激！  |
回复此楼» 猜你喜欢
国自然面上和省基金B类撒花
已经有9人回复
-
华师大读博
已经有5人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有5人回复
-
急需审稿人！！！
已经有3人回复
-
申博/考博
已经有8人回复
-
295分求调剂
已经有6人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有13人回复
-
085600材料与化工调剂
已经有6人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
-
337求调剂
已经有3人回复
11楼2016-05-18 10:56:30
2楼2012-02-27 12:15:32
3楼2012-02-27 12:16:50
凝雪小小
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 421.8
- 散金: 10
- 帖子: 105
- 在线: 61.7小时
- 虫号: 1394131
- 注册: 2011-09-07
- 性别: GG
- 专业: 控制理论与方法
4楼2012-02-27 14:54:28
5