实例化矩阵变量并评估涉及此矩阵变量的所有LMI术语
mnewsys = setmvar(lmisys,x,xval)
setmvar.
设置矩阵变量X使用标识符X
价值xval.
。所有涉及的条款X评估,相应地更新常数术语,并且X从矩阵变量列表中删除。返回所得LMI系统的描述新闻
。
整数X
是返回的标识符LMIVAR.
什么时候X被宣布为。实例化X和setmvar.
不会更改剩余矩阵变量的标识符。
功能setmvar.
冻结某些矩阵变量并针对剩余的矩阵优化是有用的。它通过避免局部或完全重新定义该组LMI约束来节省时间。
考虑系统
X=斧头+BU.
以及寻找稳定国家反馈法的问题你=kx.在哪里K.是一个未知的矩阵。
由Lyapunov定理,这相当于找到P.> 0和K.这样
(A + BK.)P.+P.(A + BK.T.)+一世<0。
随着变量的变化y:=kp.,这种情况降低到LMI
AP+paT.+经过+yT.B.T.+一世<0。
此LMI由命令输入
n =尺寸(a,1)%状态ncon = size(b,2)%的输入setlmis([])p = lmivar(1,[n 1])%p全对称y = lmivar(2,[ncon n])%y矩形Lmiterm([111 1 p],a,1,'s')%ap + pa'lmiterm([1 1 1y],b,1,'s')%by +Y'B'LMITERM([1 1 1 0],1)%I LMIS = GETLMIS
要了解此问题是否有解决方案K.对于特定的Lyapunov矩阵P.=一世, 放P.至一世通过打字
News = setmvar(LMI,P,1)
由此产生的LMI系统消息
只有一个变量y=K.。通过呼叫来评估其可行性屁股
:
[Tmin,XFEAS] = FEAP(新闻)Y = DEC2MAT(新闻,XFEAS,Y)
计算的y是可行的Tmin.
<0。