查看: 111|回复: 1|关注: 0

[未答复] 怎么通过fminsearch.m找到LMI中α的局部最优值啊以及怎么用gevp求解器求解广义特征值最小化问题

[复制链接]

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
发表于 2020-10-18 11:28:50 | 显示全部楼层 |阅读模式
本帖最后由 晨落蹊下 于 2020-10-18 21:28 编辑

这是程序
A=[-2 0;0 -0.9];
Ad=[-1 0;-1 -1];
B=[-0.5; 1];
h=0.7;
w=1;
u=0.1;
setlmis([]);
P=lmivar(1,[2 1]);
S=lmivar(1,[2 1]);
R=lmivar(1,[2 1]);
W=lmivar(1,[2 1]);
X=lmivar(1,[2 1]);
Y=lmivar(1,[2 1]);
Z=lmivar(1,[2 1]);
a=lmivar(1,[1 0]);
q=lmivar(1,[1 0]);
r=(1/h)*exp(-a*h);
t=-(1-u)*exp(-a*h);
g=1/q;
lmiterm([1 1 1 P],1,A,'s');
lmiterm([1 1 1 X],1,A,'s');
lmiterm([1 1 1 Y],1,1,'s');
lmiterm([1 1 1 P],a,1);
lmiterm([1 1 1 X],a,1);
lmiterm([1 1 1 S],1,1);
lmiterm([1 1 1 R],h,1);
lmiterm([1 1 2 P],1,Ad);
lmiterm([1 1 2 X],1,Ad);
lmiterm([1 1 2 Y],-1,1);
lmiterm([1 1 3 Y],A',1);
lmiterm([1 1 3 Z],1,1);
lmiterm([1 1 3 Y],a,1);
lmiterm([1 1 4 P],1,B);
lmiterm([1 1 4 X],1,B);
lmiterm([1 1 5 Y],1,1);
lmiterm([1 2 2 S],t,1);
lmiterm([1 2 2 W],u^2,1);
lmiterm([1 2 3 Y],Ad',1);
lmiterm([1 2 3 Z],-1,1);
lmiterm([1 3 3 R],-r,1);
lmiterm([1 3 3 Z],a,1);
lmiterm([1 3 4 -Y],1,B);
lmiterm([1 3 5 Z],1,1);
lmiterm([1 4 4 0],-a/w^2);
lmiterm([1 5 5 W],-1,1);
lmiterm([-2 1 1 X],1,1);
lmiterm([-2 1 2 Y],1,1);
lmiterm([-2 2 2 Z],1,1);
lmiterm([-2 2 2 S],r,1);
lmiterm([-3 1 1 P],1,1);
lmiterm([-4 1 1 S],1,1);
lmiterm([-5 1 1 W],1,1);
lmiterm([-6 1 1 R],1,1);
lmiterm([-7 1 1 0],1,g);
lmiterm([-7 1 2 0],1,1);
lmisterm([-7 2 2 P],1,1);
lmisys=getlmis
[g,xopt]=gevp(lmisys)
pp=dec2mat(lmisys,xfeas,P)
rr=dec2mat(lmisys,xfeas,R)
ss=dec2mat(lmisys,xfeas,S)
ww=dec2mat(lmisys,xfeas,W)

出现的错误如下
错误使用 lmiterm (line 85)
There should be no B matrix for constant terms

我用的是MATLAB 64位R2018a

这是我要求解的LMI

这是我要求解的LMI
回复主题 已获打赏: 0 积分

举报

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
 楼主| 发表于 2020-10-18 21:02:30 | 显示全部楼层
是不是要先通过fminsearch.m找到上面这个LMI中α的局部最优值啊?怎么通过fminsearch.m找到上面这个LMI中α的局部最优值啊
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表