24小时热门版块排行榜    

查看: 4188  |  回复: 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
麻烦各位指出其中有错误的地方 不胜感激!
回复此楼
我是帝都小愣神
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ep22stone

铜虫 (小有名气)

mincx语句如何用,楼主能不能介绍一下
7楼2013-03-23 22:58:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

lhjyml

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
凝雪小小(金币+18): ★★★很有帮助 2012-02-27 14:54:55
1)最优的γ²,需要加入一个循环语句。
2)无解原因:(1)首先仔细检查你的理论推导,若没有问题,检查程序,若还没有问题,在条件参数(比如b=0.1;   %ε=b)
2楼2012-02-27 12:15:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lhjyml

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
凝雪小小: 金币+7 2012-05-09 17:29:00
无解原因:再调节相关参数(比如b=0.1;   %ε=b)
3楼2012-02-27 12:16:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

凝雪小小

新虫 (小有名气)

引用回帖:
: Originally posted by lhjyml at 2012-02-27 12:15:32:
1)最优的γ&sup2,需要加入一个循环语句。
2)无解原因:(1)首先仔细检查你的理论推导,若没有问题,检查程序,若还没有问题,在条件参数(比如b=0.1;   %ε=b)

1)要加一个怎样的循环语句呢  之前没有编过  求高手指点呀
2)不确定的参数有什么试的经验么
我是帝都小愣神
4楼2012-02-27 14:54:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试 317 +6 半拉月丙 2026-03-20 6/300 2026-03-21 18:27 by 学员8dgXkO
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见