提升技术已成为采样系统设计的主要工具,但提升变换不能保证频率响应不变,故不适用于加权H∞设计。为此,针对采样控制系统的结构特点,将频率响应分为两个通道进行计算,采用经典的多速率感知环境下的状态观测器,既可以得到频率响应,又可以求得系统的L2诱导范数。方法简单直观,物理概念清楚。讨论了一类不确定时变时滞切换系统基于状态观测器的非脆弱H∞控制问题。给出基于观测器的非脆弱H∞控制器形式,在同时限定时滞上界和时滞导数上界的条件下,利用单李雅普诺夫函数方法和凸组合技术,构造李雅普诺夫泛函和其对应的切换规则;利用不等式引理对求导后的李雅普诺夫泛函进行不等式放缩进而消除不等式中的时变时滞项;再次,引入J函数得到满足H∞性能指标的非线性矩阵不等式,并利用Schur补引理将非线性矩阵不等式转化为等价的2个线性矩阵不等式,得到一类不确定变时滞切换系统存在基于状态观测器的H∞控制器且满足H∞性能指标的充分条件;给出由2个子系统构成的一个切换系统,利用LMI工具箱进行数值例子仿真说明定理的有效性和实用性。
????2运行结果
部分代码:
set;set;set;set;
al1=0.01;al2=0.01;
A=[0.950.50.2;-0.9-0.2;00.85];%StablePlant%A=[20.50.2;-0.11-0.2;00.11];%UnstablePlant
B=[1;5;2];C=[110;011];D=[10;01];
S0=[10;01];S1=[00;00];S2=[10;00];S3=[00;01];S4=[10;00];S5=[00;00];
C0=S0*C;C1=S1*C;C2=S2*C;C3=S3*C;C4=S4*C;C5=S5*C;D0=S0*D;D1=S1*D;D2=S2*D;D3=S3*D;D4=S4*D;D5=S5*D;
setlmis
[gamma2,n,sgamma2]=lmivar;
[P1,n,sP1]=lmivar;[P2,n,sP2]=lmivar;[P3,n,sP3]=lmivar;[P4,n,sP4]=lmivar;[P5,n,sP5]=lmivar;[P0,n,sP0]=lmivar;
[Y1,n,sY1]=lmivar;[Y2,n,sY2]=lmivar;[Y3,n,sY3]=lmivar;[Y4,n,sY4]=lmivar;[Y5,n,sY5]=lmivar;[Y0,n,sY0]=lmivar;
S1=newlmi;lmitermlmitermlmitermlmitermlmitermlmitermlmiterm)lmitermlmiterm)S2=newlmi;lmitermlmitermlmitermlmitermlmitermlmitermlmiterm)lmitermlmiterm)S3=newlmi;lmitermlmitermlmitermlmitermlmitermlmitermlmiterm)lmitermlmiterm)S4=newlmi;lmitermlmitermlmitermlmitermlmitermlmitermlmiterm)lmitermlmiterm)S5=newlmi;lmitermlmitermlmitermlmitermlmitermlmitermlmiterm)lmitermlmiterm)S6=newlmi;lmitermlmitermlmitermlmitermlmitermlmitermlmiterm)lmitermlmiterm)
LMIs=getlmis;c=[1zeros];
[copt,xopt]=mincx;P0opt=dec2mat;Y0opt=dec2mat;P1opt=dec2mat;Y1opt=dec2mat;P2opt=dec2mat;Y2opt=dec2mat;P3opt=dec2mat;Y3opt=dec2mat;P4opt=dec2mat;Y4opt=dec2mat;P5opt=dec2mat;Y5opt=dec2mat;
L0=-P0opt^*Y0opt;L1=-P1opt^*Y1opt;L2=-P2opt^*Y2opt;L3=-P3opt^*Y3opt;L4=-P4opt^*Y4opt;L5=-P5opt^*Y5opt;
x=0;x=0;x=0;xp=0;xp=0;xp=0;
Ts=0.01;t=0:1:600;a=0:0.01:2*pi;u=sin;
fork=1:600;ifk<=600u_proposed=al1*randn;y_proposed=al2*[randn;randn];elseu_proposed=0;y_proposed=0;endxp=A*xp+B*u+B*u_proposed;yp=C*xp+D*y_proposed;ifrem==0x=*x+B*u+L0*yp;y=C*x;elseifrem==1x=*x+B*u+L1*yp;y=C*x;elseifrem==2x=*x+B*u+L2*yp;y=C*x;elseifrem==3x=*x+B*u+L3*yp;y=C*x;elseifrem==4x=*x+B*u+L4*yp;y=C*x;elseifrem==5x=*x+B*u+L5*yp;y=C*x;enduy_proposed_sum=uy_proposed_sum+u_proposed^2+norm^2;e_proposed_sum=e_proposed_sum+norm-xp)^2;end
figure,stairs,"k","LineStyle","--"),holdon,stairs,"k","LineStyle","--"),holdon,stairs,"k","LineStyle","--")figure,stairs,"k"),holdon,stairs,"k"),holdon,stairs,"k")
ylabelxlabel
%figure,stairs-x,"k"),holdon,stairs-x,"k"),holdon,stairs-x,"k"),axis;%%ylabel%xlabel
figure%z,zQ1subplot;holdon;%gridonstairs-x,"k","LineWidth",0.;axis%set;xlabelylabelboxon%%figure%v1subplot;holdon;%gridonstairs-x,"k","LineWidth",0.;axis%set;xlabelylabelboxon
subplot;holdon;%gridonstairs-x,"k","LineWidth",0.;axis%set;xlabelylabelboxon
刘彦文,王广雄,何朕.采样系统的频率响应和L_2诱导范数[J].控制与决策,2005:1133-113DOI:13195/j.cd.2005liuyw.01
????4Matlab代码实现
文章为作者独立观点,不代表 股票程序化软件自动交易接口观点