你做的是鲁棒控制?网络控制?线性控制?
把程序发给我,我给你看看
Solving gevp Problem in Matlab30/10/10 13:51 Filed in: Matlab Code
http://arashfotouhi.com/files/92158ff4aad04c216118e30d5a6fbc71-0.html
In order to solve the following optimization problem:
minimize a
subject to
[A’P+PA+C-aP,(BP+CQ)’
BP+CQ,-I] <0
&
P=P’>0
in Matlab use the code below:
%Solving a generalized eigenvalue problem in Matlab
%Date: October 30, 2010
%Author: Arash Fotouhi
setlmis([])
p=lmivar(1,[n 1])
q=lmivar(2,[m m])
y=lmivar(1,[n 1])
lmiterm([1 1 1 0],0)
lmiterm([-1 1 1 p],1,1)
lmiterm([2 1 1 p],1,A,'s')
lmiterm([2 1 1 0],C)
lmiterm([2 2 1 p],B,1)
lmiterm([2 2 1 q],1,C)
lmiterm([2 2 2 0],-eye(m))
lmiterm([-2 1 1 y],1,1)
lmiterm([-2 2 1 0],0)
lmiterm([-2 2 2 0],0)
lmiterm([3 1 1 y],1,1)
lmiterm([-3 1 1 p],1,1)
lmis=getlmis
[lambda,decvars]=gevp(lmis,1)
Popt = dec2mat(lmis,decvars,p)
Qopt = dec2mat(lmis,decvars,q)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
解:(1)加(2)得3X1=-5,所以X1=-5/3.把X1=-5/3代入(2)和(3),得-10/3-X2加2X3=-4,(4)-20/3加X2加4X3=-2.(5).(4)加(5)得6X3=4,所以x3=2/3,代入(4)得X2=2.所以原方程组解为{X1=-5/3,X2=2,x3=2/3.欢迎分享,转载请注明来源:内存溢出
评论列表(0条)