24小时热门版块排行榜    

查看: 4679  |  回复: 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的回帖

lhjyml

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
凝雪小小: 金币+7 2012-05-09 17:29:00
无解原因:再调节相关参数(比如b=0.1;   %ε=b)
3楼2012-02-27 12:16:50
已阅   回复此楼   关注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的回帖

凝雪小小

新虫 (小有名气)

引用回帖:
: 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的回帖

凝雪小小

新虫 (小有名气)

引用回帖:
: Originally posted by lhjyml at 2012-02-27 12:16:50:
无解原因:再调节相关参数(比如b=0.1;   %ε=b)

这些参数是每个都一起调呢 还是每次只调一个呢
我是帝都小愣神
5楼2012-02-29 20:36:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见