justcc2001 发表于 2008-12-18 10:27

关于lmi中mincx求解器的问题

clear all
clc
A=[-0.08 -0.03 -0.157 0;-0.73 -0.377 0 1;0 0 0 1;0 -8.65 0 -0.5]
B=
B1=
C1=
D1=
M0=
M=
J=
q=5
s=4.5
a=0.03
b=0.02
setlmis([])
X=lmivar(1,)
Y=lmivar(2,)
r=lmivar(1,)
e11=lmivar(1,)
e12=lmivar(1,)
e13=lmivar(1,)
e14=lmivar(1,)
e21=lmivar(1,)
e22=lmivar(1,)
e23=lmivar(1,)
e24=lmivar(1,)
e25=lmivar(1,)

lmiterm(,-s,1)
lmiterm(,1,a)
lmiterm(,1,B*J*B')
lmiterm(,1,b)
lmiterm(,1,b)
lmiterm(,A,1)
lmiterm(,q,1)
lmiterm(,B,1)
lmiterm(,-s,1)
lmiterm(,1,1)
lmiterm(,1,1)
lmiterm(,1,1)
lmiterm(,1,1)
lmiterm(,-1,1)
lmiterm(,-1,inv(J))
lmiterm(,-1,1)
lmiterm(,-1,inv(J))
lmiterm(,A,1,'s')
lmiterm(,B,1,'s')
lmiterm(,1,B*J*B')
lmiterm(,1,a)
lmiterm(,1,b)
lmiterm(,1,b)
lmiterm(,B1)
lmiterm(,C1,1)
lmiterm(,D1,1)
lmiterm(,-1,1)
lmiterm(,-1,1)
lmiterm(,1,D1*J*D1')
lmiterm(,1,1)
lmiterm(,1,1)
lmiterm(,1,1)
lmiterm(,1,1)
lmiterm(,1,1)
lmiterm(,-1,inv(J))
lmiterm(,-1,1)
lmiterm(,-1,1)
lmiterm(,-1,inv(J))
lmiterm(,-1,inv(J))

lmis=getlmis
c=mat2dec(lmis,zeros(4,4),zeros(2,4),eye(1),zeros(1,1),zeros(1,1),zeros(1,1),zeros(1,1),zeros(1,1),zeros(1,1),zeros(1,1),zeros(1,1),zeros(1,1))
options=
=mincx(lmis,c,options)
Xopt=dec2mat(lmis,xopt,X)
Yopt=dec2mat(lmis,xopt,Y)
ropt=dec2mat(lmis,xopt,r)
K=inv(M0)*Yopt*inv(Xopt)

运行后出现 Error using ==> dec2mat
DECVARS must be a vector of length 28
Error in ==> hshiyan at 78
Xopt=dec2mat(lmis,xopt,X)

justcc2001 发表于 2008-12-18 10:28

第一个lmi是区域稳定,第二个是H无穷范数,单独的话都能运行,和在一起就报错了。

dadaya 发表于 2008-12-26 17:11

我的也是出现了这个问题,哪位高手可以解决一下:'(

fengsm 发表于 2009-1-7 16:27

nihao

能不能加我,想请教你一下问题

fengsm 发表于 2009-1-7 16:28

justcc2001:qq 357486273
页: [1]
查看完整版本: 关于lmi中mincx求解器的问题