《《状态空间》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《状态空间》PPT课件.ppt（46页珍藏版）》请在三一办公上搜索。
1、第八章 系统状态空间分析法,内容,系统特征方程及解关于系统可观性、可控性判别的状态反馈极点配置状态观测器的设计,.1 系统状态方程的解,状态转移矩阵,若状态方程是齐次的,即有:,EX1,a=0 1 0;0 0 1;-6-11-6;x0=1;1;1;t=0:0.1:10;for i=1:length(t)x(:,i)=expm(a*t(i)*x0;endplot3(x(1,:),x(2,:),x(3,:);grid on,系统的特征方程、特征值及特征向量,特征方程：|sI-A|=0特征值及特征向量：,V，D=eig（A）,特征向量矩阵,特征值矩阵,A*V=V*D,EX2 已知控制系统求控制系统的
2、特征方程,A=2 1-1;1 2-1;-1-1 2;I=1 0 0;0 1 0;0 0 1;syms s%符号计算det(s*I-A)s=solve(det(s*I-A)%求解,ans=s3-6*s2+9*s-4,s=4 1 1,EX2 求控制系统的特征值及特征向量,符号计算Symbolic Toolbox,EIGENSYS Obsolete Symbolic Toolbox function.V,D=EIGENSYS(A)is the same as V,D=eig(sym(A),.2 传递矩阵G,CsI-A-1B+D,A=0 1;0-2;B=1 0;0 1;C=1 0;0 1;D=0;sy
3、ms sI=1 0;0 1;G=C*inv(s*I-A)*B,G=1/s,1/s/(s+2)0,1/(s+2),.3 线性变换,状态方程的线性变换 ss2ss(sys,T),EX3,A=0-2;1-3;B=2 0;C=0 3;P=6 2;2 0;%变换矩阵x=PzP1=inv(P);A1=P1*A*P%z坐标系的模型B1=P1*BC1=C*P,A1=0 1-2-3B1=0 1C1=6 0,The eigenvalues of system are unchanged by the linear transformation:(线性变换不改变系统的特征值）,约当标准形,canon(sys,mod
4、el)canon(sys,companion),EX4利用特征值及范德蒙特矩阵求约当阵,A=0 1 0;0 0 1;2-5 4;V,D=eig(A)P=1 0 1;1 1 2;1 2 4P1=inv(P)；J=P1*A*P,V=-0.5774 0.5774-0.2182-0.5774 0.5774-0.4364-0.5774 0.5774-0.8729D=1.0000 0 0 0 1.0000 0 0 0 2.0000P=1 0 1 1 1 2 1 2 4J=1 1 0 0 1 0 0 0 2,符号计算,Jo=jordan(A),Jo=2 0 0 0 1 1 0 0 1,.4 系统的可控性和可
5、观性,MATLAB提供函数分别计算能控性矩阵和能观测性矩阵可控性矩阵CO=ctrb(A,B)可观测性矩阵OB=obsv(A,C),可控性判定,A=1 1 0;0 1 0;0 1 1;B=0 1;1 0;0 1;n=length(A)CO=ctrb(A,B);rCO=rank(CO);if rCO=n disp(System is controllable)elseif rCOn disp(System is uncontrollable)end,n=3 CO=0 1 1 1 2 1 1 0 1 0 1 0 0 1 1 1 2 1 rCO=2 System is uncontrollable,可
6、观测性判定,A=-3 1;1-3;B=1 1;1 1;C=1 1;1 1;D=0;n=length(A);OB=obsv(A,C);rOB=rank(OB)if rOB=n disp(System is observable)elseif rOBn disp(System is unobservable)end,OB=1 1 1 1-2-2-2-2rOB=1System is unobservable,可控标准形,若S为非奇异，逆矩阵存在，设为,则，变换矩阵为P,A=-2 2-1;0-2 0;1-4 0;B=0 1 1;n=length(A);CAM=ctrb(A,B);if det(CAM)
7、=0 CAM1=inv(CAM);endP=CAM1(3,:);CAM1(3,:)*A;CAM1(3,:)*A*A;P1=inv(P);A1=P*A*P1B1=P*B运行后,可知系统的可控标准形矩阵对为:A1=0 1 0 0 0 1-2-5-4,B1=0 0 1,可观测标准形,则，变换矩阵为M=PT,若V为非奇异，逆矩阵存在，设为,.5 系统状态反馈与状态观测器,利用反馈结构，研究在什么条件下能实现闭环系统极点的任意配置，以达到预期要求。状态反馈与状态观测器原理参见线性控制系统工程Module24，25,24.1 The Structure of State Space Feedback Co
8、ntrol（状态反馈控制的结构）,1.State Variable Feedback Control System,n the number of state variable,If the desired location of the closed-loop poles are,the desired characteristic equation will be,We can obtained,to make the closed-loop poles to be located in desired position.,The principle of designing a stat
9、e space controller,2.The sufficient and necessary condition of state feedback for closed-loop placement:(状态反馈实现极点配置的充要条件）,极点配置:,MATLAB直接用于系统极点配置计算的函数有acker和place,A,B为系统矩阵,K=acker(A,B,P),P为期望极点向量,K反馈增益向量,K=place(A,B,P)K,prec,message=place(A,B,P)Pre为实际极点偏离期望极点位置的误差Message是当系统某一非零极点偏离期望值位置大于10%时给出的警告信息
10、.,极点配置步骤:,(1)获得系统闭环的状态空间方程.(2)根据系统性能要求,确定期望极点分布P.(3)利用极点配置设计函数求取系统反馈增益K.(4)检验系统性能.,例:已知控制系统的系数矩阵为:A=-2-2.5-0.5;1 0 0;0 1 0,B=1 0 0,闭环系统的极点为S=-1,-2,-3,对其进行极点配置.,%极点配置A=-2,-2.5,-0.5;1,0,0;0,1,0;B=1,0,0;P=-1,-2,-3;K=acker(A,B,P)%闭环系统矩阵Ac=A-B*K%验证闭环系统的特征值eig(Ac)K=4.0000 8.5000 5.5000验证后发现配置结果正确,所以反馈控制器为
11、K=4 8.5 5.5,例:已知控制系统的传递函数为:判断系统的可控性并设计反馈控制器,使得闭环系统的极点为-2,-1 i,.,%判断系统的可控性n1=10;d1=conv(conv(1,0,1 1),1,2);a,b,c,d=tf2ss(n1,d1);n=3;CAM=ctrb(a,b);if det(CAM)=0;rcam=rank(CAM);if rcam=n,disp(系统可控)elseif rcamn,disp(系统不可控)endelseif det(CAM)=0,disp(系统不可控)End程序运行结果为:系统可控.,由于系统是可控的,所以可以任意配置系统的极点.%极点配置P=-2,
12、-1+i,-1-i;K=place(a,b,P);%显示极点配置信息K,Prec,Mes=place(a,b,P)%绘制系统阶跃响应曲线sys=ss(a-b*K,b,c,d);poles=pole(sys);%求极点,验证step(sys/dcgain(sys),10);%t为10秒的阶跃响应Axis(0,10,0,1.5),运行结果:,K=1.0000 4.0000 4.0000Prec=15Mes=,由运算结果知道,配置过程中没有出错和警告信息.由反馈系统的阶跃响应曲线可知,闭环系统的动态性能良好.所以状态反馈控制器为:K=1 4 4,25.1 Observer A model of th
13、e system under study,The approach taken to solve the problem is as following:To construct a model of the system under study;Assume(subject to certain restrictions)that the computed state variables are good approximations to the true state variables;From these computed state variables,a suitable cont
14、roller for the actual system may be constructed using the techniques described in Module 24.,Where,x are assumed to be unmeasured directly.,状态观测器设计,To decrease the error,(that is error),we take to correct to make well approach:,Select the matrix K to make the solution of this equation on error be co
15、nvergent(收敛的),then,The gain matrix K is written as:,The closed-loop poles of this model(observer)can be selected by selecting the gain matrix K,so that the state variables will be same as in the end.Hence,we can use as the state variables in the state variable feedback system.,The closed-loop system
16、 with observer,B,C,A,+,+,u,x,y,B,C,A,+,+,K,-,G,-,r,+,+,状态观测器,状态反馈,25.2 The sufficient and necessary condition of constructing a state variable observer,Observability criterion:A system A,C is state observable if and only if,状态观测器设计一般原理归结为使用极点配置法求观测器的增益矩阵G,A,C为系统矩阵,G=acker(A,C,P),P为观测器的期望极点向量,G为观测器增益
17、向量,K=place(A,C,P),例:设系统的状态空间表达式为:,试设计一个状态观测器,使极点为-3,-4,-5.解:(1)判断系统的能观性.首先检验系统是否完全能观,程序如下:A=0,0,2;1,0,9;0,1,0;B=3,2,1;C=0,0,1;ob=obsv(A,C);roam=rank(ob);if roam=n disp(系统可观)elseif roam=n disp(系统不可观)end,程序结果为:系统可观表示该系统是可观的,存在状态观测器.,(2)设计状态观测器.使用MATLAB的acker()函数来设计状态观测器,程序如下:P=-3,-4,-5;A1=A;B1=C;K=ack
18、er(A1,B1,P);H=(K)ahc=A-H*C程序结果如下:,H=62 56 12ahc=0 0-60 1 0-47 0 1-12即系统的状态观测器为:,已知离散控制系统如下:,采样周期为0.01秒,试设计状态观测器,使极点为0.1 0.1i.解:对于离散控制系统,函数acker和 place的调用方法基本相同.(1)判断系统的可观测性程序如下:,%判断系统的能观性A=0,1;-2,-1;B=0,1;C=1,1;D=0;T=0.01;sys=ss(A,B,C,D,T);ob=obsv(sys);unob=length(A)-rank(ob);if unob=0 disp(系统可观)elseif roam=0 disp(系统不可观)end运行结果为:系统可观,(2)设计状态观测器.%极点配置P=0.1+0.1*i;0.1-0.1*i;A1=A;B1=C;K=acker(A1,B1,P);H=(K)ahc=A-H*C由运行结果可知:系统的状态观测器为:,总结,
