基本信息
源码名称:状态反馈(state feedback.m)
源码大小:4.14KB
文件格式:.m
开发语言:MATLAB
更新时间:2021-01-06
   友情提示:(无需注册或充值,赞助后即可获取资源下载链接)

     嘿,亲!知识可是无价之宝呢,但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下,绝对物超所值哦!如有下载和支付问题,请联系我们QQ(微信同号):813200300

本次赞助数额为: 2 元 
   源码介绍

%方法2.2:
%将ebxl1,ebxl2,ebxl3看成变量
setlmis([])

X=lmivar(1,[4 1]);
Z=lmivar(2,[2 4]);
rs=lmivar(1,[1 0]);
ebxl1=lmivar(1,[1 0]);
ebxl2=lmivar(1,[1 0]);
ebxl3=lmivar(1,[1 0]);

%first LMI
lmiterm([1 1 1 X],A,1,'s');
lmiterm([1 1 1 Z],B,1,'s');
lmiterm([1 1 1 rs],1,1);
lmiterm([1 1 1 ebxl2],rA^2,1);
lmiterm([1 1 1 ebxl3],rB^2,1);

lmiterm([1 2 1 X],Fa,1);
lmiterm([1 2 1 Z],Fb,1);

lmiterm([1 3 1 X],1,1);
lmiterm([1 4 1 Z],1,1);
lmiterm([1 2 2 ebxl1],-1,1);
lmiterm([1 3 3 ebxl2],-1,1);
lmiterm([1 4 4 ebxl3],-1,1);

%second LMI X>0
lmiterm([-2 1 1 X],1,1);

%third LMI  rs>0
lmiterm([-3 1 1 rs],1,1);

%forth LMI ebxl1>0
lmiterm([-4 1 1 ebxl1],1,1);

%fifth LMI ebxl2>0
lmiterm([-5 1 1 ebxl2],1,1);

%sixth LMI ebxl3>0
lmiterm([-6 1 1 ebxl3],1,1);

lmis=getlmis;
[tmin,xfeas]=feasp(lmis)
xx=dec2mat(lmis,xfeas,X);
zz=dec2mat(lmis,xfeas,Z);
rr=dec2mat(lmis,xfeas,rs);
EBXL1=dec2mat(lmis,xfeas,ebxl1);
EBXL2=dec2mat(lmis,xfeas,ebxl2);
EBXL3=dec2mat(lmis,xfeas,ebxl3);

K=zz*inv(xx);
region=(rr/EBXL1)^0.5;