[求助]LMI求解中，若k×k'＝y，已知y，怎样求k？

grace_105 · 发表于 2006-6-23 16:36

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

在做输出反馈控制器设计时，对于Xdot＝AX+BU Y=CX 设计一个静态输出反馈控制器U=kY时，得到的k是个常数，而需要的k是个1×2的才满足状态模型，怎样得到1×2维的k值呢？ 程序如下： a=[0 0 0 0;0 -0.088317 0.625 0;0 0 0 1;0 -0.002725 0.3283625 0]; b=[0;-0.883167;0;2.356553]; c=[1 0 0 0;0 0 1 0]; 
 setlmis([]); x=lmivar(1,[4 1]); y=lmivar(1,[1 1]); ew=lmivar(2,[1 1]); lmiterm([1 1 1 x],1,a','s'); % LMI #1: x*a'+a*x lmiterm([1 2 1 0],b'); % LMI #1: b' lmiterm([1 2 2 y],1,-1); % LMI #1: -y lmiterm([1 3 1 x],c,1); % LMI #1: c*x lmiterm([1 3 3 0],-inv(ew)); % LMI #1: -inv(ew) lmiterm([-2 1 1 x],1,1); % LMI #2: x 即令正定x>0，默认为<0，所以用-2表示第二个LMI大于0 lmiterm([-3 1 1 ew],1,1); % LMI #3: ew 即令α>0,是正实数，此处将α用ew表示，-3表示第三个LMI大于0 s1=getlmis; [lp kp]=feasp(s1) %此处用来求s1的可行解，返回的lp<0，则该LMI可行。 x=dec2mat(s1,kp,x) y=dec2mat(s1,kp,y); ew=dec2mat(s1,kp,ew) p=inv(x) % 到此已求出了x,y,ew,p的值，以下要用正交补求k y=y/ew; %由ew×inv（y）除以.ew得到inv（y） y=inv(y) %再求逆得到所求的y 此处y＝0.1505 z=orth(y); t=chol(y); k=(z*t)' ％此处得到k＝0.3879，我想得到一个1×2的k，怎样修改程序呢？

[此贴子已经被作者于2006-6-23 16:42:51编辑过]

realyyy · 发表于 2006-6-23 21:40

共同关注ing （楼主的意思是y已知（y为一常数），且k*k'=y，如何求一个维数为1x2的k）

[此贴子已经被作者于2006-6-23 22:42:11编辑过]

账号		自动登录	找回密码
密码			我要加入

[编程技巧] [求助]LMI求解中，若k×k'＝y，已知y，怎样求k？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。