【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 3370  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

幸运天使

铜虫 (小有名气)

楼主,你好!

我想请问你的问题解决了没有?如果解决了,可否把代码发给我看看,我也遇到了相同的问题,不甚感激!
8楼2014-05-14 14:53:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dc652037717

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
凝雪小小: 金币+5 2012-04-16 10:05:47
求最优的r值应该是用mincx语句,如果只用feasp语句得加上循环语句。
6楼2012-04-15 13:30:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

凝雪小小

新虫 (小有名气)

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

ep22stone

铜虫 (小有名气)

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

大母虫子

银虫 (初入文坛)

LZ问题解决了吗?过去好久了
术业有专攻
9楼2015-01-17 14:58:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guaixiaofu

新虫 (初入文坛)

我有一个问题,怎么才能传上去
10楼2015-08-27 16:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 凝雪小小 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 请问教育部人文社科难度大吗 +9 锦衣卫寒战 2024-04-15 14/700 2024-04-16 23:10 by 锦衣卫寒战
[考研] 资源与环境求调剂 +15 shallowshiyi 2024-04-11 15/750 2024-04-16 15:12 by nanaliuyun
[考研] 调剂 +8 人间十三月@ 2024-04-10 10/500 2024-04-15 22:01 by zzzhonghl
[考研] 322求调剂 +6 本己上岸 2024-04-15 7/350 2024-04-15 14:19 by mthwyj
[材料综合] 请问哪里可以测试低压吸氢的PCT曲线 +3 yunshengcd 2024-04-14 5/250 2024-04-15 13:59 by zqdsb
[考研] 一志愿南京理工大学材料工程282求调剂 +13 abc16785 2024-04-11 13/650 2024-04-15 11:26 by coolfishwll
[考研] 化工学硕调剂276 +13 星星陨落 2024-04-10 14/700 2024-04-15 10:45 by 奇悠悠
[考研] 汕头大学 -资源与环境-需要较大量调剂-报考07海洋、生物、生态三个专业也可以申请 +6 stusam 2024-04-11 10/500 2024-04-13 23:49 by stusam
[考研] 290,一志愿报考重庆大学生物与医药专硕,求调剂 +4 GPX4 2024-04-13 5/250 2024-04-13 17:19 by lincunhui
[基金申请] 求问大神,教育部人文社科还有其他基金项目组成员是不是一定要有教授成员 +3 nanco123 2024-04-11 3/150 2024-04-13 11:56 by ljh9201
[教师之家] 表弟博士毕业,选30万安家费还是副教授??? +10 yingyu427 2024-04-11 12/600 2024-04-13 11:16 by karl_messi
[考研] 生物与医药调剂 +8 lilinxian 2024-04-12 8/400 2024-04-12 23:35 by heyan奋斗
[考研] 281求调剂 +9 矢志不渝上岸 2024-04-11 10/500 2024-04-12 14:02 by 上岸xd
[考研] 0858能源动力求调剂 +5 毛驴子秋裤 2024-04-11 11/550 2024-04-12 11:41 by 想上岸啊a
[考研] 双9 +5 随鱼而安. 2024-04-11 5/250 2024-04-11 16:40 by 公瑾逍遥
[考研] 希望大家都上岸 +5 1537507 2024-04-10 10/500 2024-04-11 09:34 by zhao200155
[考研] 求调剂,双9, +8 随鱼而安. 2024-04-10 8/400 2024-04-10 22:30 by 1464969820
[考研] 322求调剂 +7 梦南洲 2024-04-10 7/350 2024-04-10 22:27 by 1464969820
[考研] 材化专硕,一志愿中科大284求调剂 +16 玉上花 2024-04-10 17/850 2024-04-10 15:08 by qlm5820
[考研] 求调剂 +4 吃了一口德芙 2024-04-10 4/200 2024-04-10 13:48 by 缘子你好
信息提示
请填处理意见