应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 运筹学与控制论 » matlab线性矩阵不等式求解做优化问题
51/1返回列表
查看: 3902  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

凝雪小小

新虫 (小有名气)

[求助] 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
麻烦各位指出其中有错误的地方 不胜感激！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
我是帝都小愣神
1楼2012-02-26 22:12:31
已阅   回复此楼   关注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的回帖
查看全部 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的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭