24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4276  |  回复: 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的回帖

guaixiaofu

新虫 (初入文坛)

我有一个问题,怎么才能传上去
10楼2015-08-27 16:36:26
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 苦命人。。。 2026-04-18 4/200 2026-04-19 02:01 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 291求调剂 +10 关忆北. 2026-04-14 10/500 2026-04-18 23:32 by 路病情
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 化工学硕294分,求导师收留 +33 yzyzx 2026-04-12 37/1850 2026-04-17 23:00 by wunaiy88
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 294求调剂 +7 淡然654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 105500药学求调剂 +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
信息提示
请填处理意见