《第1章系统的状态空间法》由会员分享,可在线阅读,更多相关《第1章系统的状态空间法(101页珍藏版)》请在金锄头文库上搜索。
1、绪绪 论论一、现代控制理论的性质及发展一、现代控制理论的性质及发展 控制理论研究的问题:控制理论研究的问题:控制理论研究的问题:控制理论研究的问题:如何改进系统的动如何改进系统的动态性能,达到所需的性能指标。态性能,达到所需的性能指标。 控制系统中两个重要的概念:控制系统中两个重要的概念:控制系统中两个重要的概念:控制系统中两个重要的概念:1、反馈的概念、反馈的概念 2、最优控制的概念、最优控制的概念控制理论发展的三个时期:控制理论发展的三个时期:控制理论发展的三个时期:控制理论发展的三个时期: 1、经典控制理论时期(二十世纪、经典控制理论时期(二十世纪3050年代)年代)研究对象主要是线性系
2、统,以拉氏变换为数学工具。研究对象主要是线性系统,以拉氏变换为数学工具。较好的解决了单输入单输出反馈控制问题。较好的解决了单输入单输出反馈控制问题。 2、现代控制理论时期(二十世纪、现代控制理论时期(二十世纪5070年代)年代)研究对象为多变量、非线性、时变、离散系统。研究对象为多变量、非线性、时变、离散系统。 以线性代数和微分方程为主要的数学工具,以以线性代数和微分方程为主要的数学工具,以状态空间法为基础,分析和设计控制系统。状态空间法为基础,分析和设计控制系统。 3、大系统理论和智能控制理论时期(二十世、大系统理论和智能控制理论时期(二十世纪纪70年代至今)年代至今)二、现代控制理论基础主
3、要内容二、现代控制理论基础主要内容1、线性系统理论、线性系统理论2、系统辨识、系统辨识3、最优控制、最优控制4、最优估计、最优估计5、自适应控制、自适应控制1、线性系统理论、线性系统理论 建立系统的状态方程,系统的响应特性,系建立系统的状态方程,系统的响应特性,系统的稳定性,能控性,能观测性,状态反馈,状统的稳定性,能控性,能观测性,状态反馈,状态观测器态观测器2、系统辨识、系统辨识包括结构辨识和参数辨识包括结构辨识和参数辨识3、最优控制、最优控制 通过观测一个系统的输入输出关系来确定其数学通过观测一个系统的输入输出关系来确定其数学模型的方法。模型的方法。 在已知系统状态方程、初始条件及某些约
4、束条在已知系统状态方程、初始条件及某些约束条件下,寻找一个最优控制量,使系统的状态或输出件下,寻找一个最优控制量,使系统的状态或输出在控制向量作用下,使某一指标达到最优值。在控制向量作用下,使某一指标达到最优值。4、最优估计、最优估计在在通讯工程中,接受到的信号为:通讯工程中,接受到的信号为:Y(t)=S(t)+ (t)有用信号有用信号干扰躁声干扰躁声由由Y(t)求)求S(t)的估计的估计S(t)5、自适应控制、自适应控制 自适应控制一般分为两类:模型参考自适应自适应控制一般分为两类:模型参考自适应控制,自校正自适应控制。控制,自校正自适应控制。 当控制对象的结构或参数随环境条件的变化而有大当
5、控制对象的结构或参数随环境条件的变化而有大的变化时,为了保证控制系统在整个控制过程中都满足的变化时,为了保证控制系统在整个控制过程中都满足某一最优准则,则最优控制器的参数就要随之加以调节,某一最优准则,则最优控制器的参数就要随之加以调节,这类控制为自适应控制。这类控制为自适应控制。四、本课程主要内容四、本课程主要内容四、本课程主要内容四、本课程主要内容1、状态空间法、状态空间法2、动态分析、动态分析3、能控、能控性与性与能观测性能观测性三、控制理论的应用三、控制理论的应用航天与航空、电机械、化工、冶金、交通、医疗航天与航空、电机械、化工、冶金、交通、医疗4、结构分解与实现、结构分解与实现5、稳
6、定性分析、稳定性分析6、状态反馈、状态反馈7、最优控制、最优控制8、最小值原理、最小值原理五、参考书五、参考书1、现代控制理论基础现代控制理论基础机械工业出版社机械工业出版社 常春馨编常春馨编2、现代控制理论基础现代控制理论基础北京工业大学出版社北京工业大学出版社 谢克明编谢克明编3、现代控制理论现代控制理论机械工业出版社机械工业出版社 刘豹编刘豹编4、现代控制理论基础现代控制理论基础电子工业出版社电子工业出版社 尤昌德编尤昌德编第第1 1章章 控制系统的状态空间表达式控制系统的状态空间表达式1.1 概述概述1.2 控制系统的状态空间表达式控制系统的状态空间表达式1.3 状态空间表达式的建立状
7、态空间表达式的建立1.4 状态方程的线性变换状态方程的线性变换1.5 系统的传递函数阵系统的传递函数阵1.6 离散系统的状态空间表达式离散系统的状态空间表达式1.7 时变系统和非线性系统的状时变系统和非线性系统的状 态空间表达式态空间表达式1.1 概述概述古典控制理论是基于传递函数来分析与设计系统。古典控制理论是基于传递函数来分析与设计系统。现代控制理论是建立在状态空间法基础上。现代控制理论是建立在状态空间法基础上。1.2 控制系统的状态空间表达式控制系统的状态空间表达式1 1. .2 2 . . 1 1 基本概念基本概念基本概念基本概念1、系统的状态:系统运动信息的集合,表示系统、系统的状态
8、:系统运动信息的集合,表示系统 过去、现在、将来的运动状况。过去、现在、将来的运动状况。2、系统的状态变量:唯一确定系统状态的一组独、系统的状态变量:唯一确定系统状态的一组独立变量。能够完全描述系统时域行为的最小变量立变量。能够完全描述系统时域行为的最小变量组。状态变量的选取不唯一。组。状态变量的选取不唯一。1 1. .2 2 . . 1 1 基本概念基本概念基本概念基本概念3、状态矢量:以、状态矢量:以n个状态变量为分量,构成一个个状态变量为分量,构成一个n维矢量。维矢量。 X(t)=x1(t)x2(t) : : xn(t)4、状态空间、状态空间 :以以n个状态变量为坐标轴所构成的空个状态变
9、量为坐标轴所构成的空间,称为间,称为n维状态空间。维状态空间。 5、状态方程、状态方程 :状态变量的一阶导数与输入变量及:状态变量的一阶导数与输入变量及状态变量的关系式。状态变量的关系式。dx1dt=f1(x1 x2 u1 u2)一阶微分方程一阶微分方程6、输出方程、输出方程 :输出变量与输入变量及状态变量的:输出变量与输入变量及状态变量的关系式。关系式。1 1. .2 2 . . 1 1 基本概念基本概念基本概念基本概念 代数方程代数方程7、状态空间表达式、状态空间表达式 :状态方程和输出方程。:状态方程和输出方程。1 1. .2 2 . . 2 2 控制系统状态空间表达式控制系统状态空间表
10、达式控制系统状态空间表达式控制系统状态空间表达式例:某机械运动系统的物理模型,它是一个弹簧例:某机械运动系统的物理模型,它是一个弹簧质质量量阻尼系统,试建立输入的外力阻尼系统,试建立输入的外力u (t),输出为位移输出为位移 y (t)的的状态空间表达式。状态空间表达式。fmkuy y1= f1(x1 x2 u1 u2)K:弹性系数弹性系数f:阻尼系数阻尼系数1 1. .2 2 . . 2 2 控制系统状态空间表达式控制系统状态空间表达式控制系统状态空间表达式控制系统状态空间表达式fmk解:系统的运动方程:解:系统的运动方程:d2ydt 2m=ufdydtkyd2ydt 2m+ fdydt+k
11、y = u系统的状态变量:系统的状态变量:x1=yuy x2=y=x1x2=y系统的状态方程:系统的状态方程:x1 = x2系统的输出方程:系统的输出方程: y = x1x2=kmx1fmx21mu+1 1. .2 2 . . 2 2 控制系统状态空间表达式控制系统状态空间表达式控制系统状态空间表达式控制系统状态空间表达式fmkuy矩阵形式:矩阵形式:x1 = x2y = x1y =1 0x1x2x1x2=u1m0 1kmfmx1x2+x1x20简写为:简写为:X=AX+buY=CXx2=kmx1fmx21mu+1 1. .2 2 . . 2 2 控制系统状态空间表达式控制系统状态空间表达式控
12、制系统状态空间表达式控制系统状态空间表达式多多输入多输出线性定常系统:输入多输出线性定常系统:x1x2 : :xn=Xu1u2 : :ur=Uy1y2 : :ym=YX=AX+BUY=CX+DUA=a11a12. a1n a21a22. a2n . . . . . an1an2. ann B=b11b12. b1r b21b22. b2r . . . . . bn1bn2. bnr1 1. .2 2 . . 2 2 控制系统状态空间表达式控制系统状态空间表达式控制系统状态空间表达式控制系统状态空间表达式X=AX+BUY=CX+DUC=c11c12. c1n c21c22. c2n . . .
13、. . cm1cm2. cmn B=b11b12. b1r b21b22. b2r . . . . . bn1bn2. bnr A=a11a12. a1n a21a22. a2n . . . . . an1an2. ann D=d11d12. d1rd21d22. d2r . . . . . dm1dm2. dmr 1 1. .2 2 . . 3 3 控制系统状态空间的一般表达式控制系统状态空间的一般表达式控制系统状态空间的一般表达式控制系统状态空间的一般表达式X=F(X U t)Y=G(X U t)1 1. .2 2 . . 4 4 线性系统状态空间表达式的模拟结构图和线性系统状态空间表达式
14、的模拟结构图和线性系统状态空间表达式的模拟结构图和线性系统状态空间表达式的模拟结构图和 信号流图信号流图信号流图信号流图X=AX+BUY=CX+DU BDCAU(t)+ Y(t)DUAXCXX比例器比例器加法器加法器积分器积分器1、结构图、结构图BUX1 1. .2 2 . . 4 4 线性系统状态空间表达式的模拟结构图和线性系统状态空间表达式的模拟结构图和线性系统状态空间表达式的模拟结构图和线性系统状态空间表达式的模拟结构图和 信号流图信号流图信号流图信号流图例:线性系统的状态空间表达式为例:线性系统的状态空间表达式为x1=x2x2=x3x3=8x114x27x3+uy =x1+2x2 试画
15、出它的系统结构图。试画出它的系统结构图。解:这是一个三阶系统,需解:这是一个三阶系统,需3个积分器个积分器例:线性系统的状态方程为例:线性系统的状态方程为x1=x2x2=x3x3=8x114x27x3+uy =x1+2x2 试画出它的系统结构图。试画出它的系统结构图。解:这是一个三阶系统,需解:这是一个三阶系统,需3个积分器个积分器x1x3x3=x2x2=x18147+2+ uyx1x3x3=x2x2=x18147+2+ uy2、信号流图、信号流图X=AX+BUY=CX+DUUBAC DYXX将将上例中的结构图用信号流图表示上例中的结构图用信号流图表示x2=x1x3=x2x3x18147+2+
16、 uy2、信号流图、信号流图将将上例中的结构图用信号流图表示上例中的结构图用信号流图表示u 2 y1x3x2x1x1114871.3 状态空间表达式的建立状态空间表达式的建立1 1. .3 3 .1 .1由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式方框图方框图结构图结构图例:试建立系统的状态空间表达式例:试建立系统的状态空间表达式解:将惯性环节变为积分环节解:将惯性环节变为积分环节 k1T1S+1 k2T2S+1k4 k1T1S+1uy+ k3 T3S k1 T1 1 S+1 T11 1. .3 3 .1 .1由系统方
17、框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式 k1 T1解:将惯性环节变为积分环节解:将惯性环节变为积分环节 k1T1S+1 1 S+1 T1 k1 T1 1 S 1 T1+ k1 T1 1 T1+x3x31 1. .3 3 .1 .1 由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式 k2T2S+1k4 k1T1S+1uy+ k3 T3S 1 T1 +x3x3 k1 T1 k4 1 T2 +x2x2 k2 T2 k3 T3 x1x1=y+u1 1. .3
18、 3 .1 .1 由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式 1 T1 +x3x3 k1 T1 k4 1 T2 +x2x2 k2 T2 k3 T3 x1x1=y+ux1= x2k3 T3x3= x3+ (u k4x1)k1 T11 T1x2= x2+ x31 T2k2 T2y =x11 1. .3 3 .1 .1 由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式例:含有零点例:含有零点 kSSZS+Puy+ 1S+aZ +PS+PSZS+P=1Z
19、 +PS+P kSuy+ 1S+a+1 1. .3 3 .2 .2 由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式例:由例:由RLC组成的系统如图,组成的系统如图,u为输入变量,为输入变量,y为输出变为输出变 量,试建立它的状态空间表达式。量,试建立它的状态空间表达式。CRLuRuLuCiu+y解:解:u= uR + uL + uC uL =Ldt di u= RCduCdt+LCd2uCdt2+uC设状态变量为:设状态变量为:x1=uC、x2= x1= uCy =1 0x1x2 x1x2=u1LC0 1
20、 1LCRLx1x2+x1x20i =Cdt duC1 1. .3 3 .2 .2 由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式例:试建立图中所示的机械旋转运动的状态空间表例:试建立图中所示的机械旋转运动的状态空间表 达式。设转动惯量为达式。设转动惯量为J。设状态变量为:设状态变量为:x1= 、x2= = BKTB:粘性阻尼系数。粘性阻尼系数。 K: 扭转轴的刚性系数。扭转轴的刚性系数。T:施加于扭转轴上的力矩。施加于扭转轴上的力矩。 :转动的角度。:转动的角度。解:设扭转轴的转动角度解:设扭转轴的转动
21、角度 及其角速度及其角速度 为状态变量。为状态变量。u=T根据牛顿定律:根据牛顿定律:1 1. .3 3 .2 .2 由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式例:试建立图中所示的机械旋转运动的状态空间表例:试建立图中所示的机械旋转运动的状态空间表 达式。设转动惯量为达式。设转动惯量为J。y =1 0x1x2 x1x2=u1 J0 1 kJBJx1x2+x1x20 BKTB:粘性阻尼系数。粘性阻尼系数。 K: 扭转轴的刚性系数。扭转轴的刚性系数。T:施加于扭转轴上的力矩。施加于扭转轴上的力矩。 :转动
22、的角度。:转动的角度。解:设扭转轴的转动角度解:设扭转轴的转动角度 及其角速度及其角速度 为状态变量。为状态变量。1 1. .3 3 .3 .3 由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式1、输入函数不包含导数项时、输入函数不包含导数项时设设系统的微分方程:系统的微分方程:y(n)+a1y(n1)+an1 y+any=bu变换为:变换为:X=AX+BUY=CX设设y、 yy(n1)为为系统的状态变量。系统的状态变量。令:令:x1=yx2=yxn1=y(n2)xn=y(n1)系统状态方程:系统状态方程:
23、x1=x2x2=x3xn1= xnxn= y(n)=a1xn a2xn1 anx1+bu1 1. .3 3 .3 .3 由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式1、输入函数不包含导数项时、输入函数不包含导数项时系统状态方程:系统状态方程:x1=x2x2=x3xn1= xnxn= y(n)=a1xn a2xn1 anx1+buy =x1 =0 1 0 0x1x2xnu0x1x2+x1x20xn0b0 0 1 0an an1 an2 a1 y =x11 1. .3 3 .3 .3 由系统的微分方程建立状
24、态空间表达式由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式y =10 0x1x2x1x2xna3a2a1x1x2x3=u0 1 0+x1x2x30060 0 16 11 5y = 1 0 0x1x2x3例:将例:将y+5y+11y+6y=6u变换为状态空间表达式变换为状态空间表达式 2、输入函数包含导数项时、输入函数包含导数项时设设系统的微分方程:系统的微分方程:y(n)+a1y(n1)+an1 y+any=b0u(n)+b1u(n1) +bn1u +bnu设:设:y+a1y+a2y+a3y=b0 u+ b1 u+b2u+ b3u 状态空
25、间表达式状态空间表达式x1x2x3=u0 1 0+x1x2x3c1c2c30 0 1a3 a2a1选择待定系数选择待定系数c1 、c2 、c3使使状态方程中不含导数项状态方程中不含导数项2、输入函数包含导数项时、输入函数包含导数项时x1x2x3=u0 1 0+x1x2x3c1c2c30 0 1a3a2a1x1= x2+ c1 ux2= x3+ c2 ux3=a3 x1 a2 x2 a1 x3 + c3 u将将上式展开:上式展开:求求c1 、c2 、c32、输入函数包含导数项时、输入函数包含导数项时令:令:y= x1+c0u (1)x1= x2+ c1 ux2= x3+ c2 ux3=a3 x1
26、 a2 x2 a1 x3 + c3 uy= x1+c0u = x2+c1u + c0u (2)y= x2+c1u + c0u= x3+c2u + c1u+c0u (3)y= x3+c2u+c1u + c0u = a3 x1 a2 x2 a1 x3+ c3 u +c2u + c1u+c0u (4)a1 (3)+ a2 (2)+ a3 (1)+(4)即:即: y+a1y+a2y+a3y=b0 u+ b1 u+b2u+ b3u 左左式式=c0u +(c1+a1c0)u+(c2+a1c1+ a2c0 )u+(c3+a1c2+ a2c1+a3c0)u 2、输入函数包含导数项时、输入函数包含导数项时y+a
27、1y+a2y+a3y=b0 u+ b1 u+b2u+ b3u 左左式式=c0u +(c1+a1c0)u+(c2+a1c1+ a2c0 )u+(c3+a1c2+ a2c1+a3c0)u 比较系数得:比较系数得:c0= b0c1=b1a1c0c2=b2a1c1 a2c0c3=b3a1c2 a2c1 a3c0对于对于n阶系统:阶系统:cn=bna1cn1a2c n2 aic ni anc02、输入函数包含导数项时、输入函数包含导数项时求系统的状态变量求系统的状态变量y= x1+c0u (1)y= x1+c0u = x2+c1u + c0u (2)y= x2+c1u + c0u= x3+c2u + c
28、1u+c0u (3)x1 = y c0u (1)x2 = y c1u c0u (2)x3 = y c2u c1u c0u(3)因为:因为:所以:所以:状态变量是由状态变量是由y、u及它的各价导数组成。及它的各价导数组成。解:解:c0=0a3a2a1例:将例:将y+4y+2y+y=u+u+3u变换为状态空间表达式变换为状态空间表达式 b1b2b3b0=0c1=b1a1c0=140=1c2=b2a1c1 a2c0=1 41= 3c3=b3a1c2 a2c1 a3c0=3 4( 3)21=13x1x2x3=u0 1 0+x1x2x313130 0 11 2 4y= x1作业:作业:1-1试求系统的模
29、拟结构图,并建立状态空间表达式。试求系统的模拟结构图,并建立状态空间表达式。 k1T1S+1 k2 Suy+ 1 T2S+1+ 1 S1-2 将将y+2y+4y+6y=2u变换为状态空间表达式。变换为状态空间表达式。 1-3 将将变换为状态空间表达式。变换为状态空间表达式。1-3 试建立图中所示的机械旋转运动的状态空间表试建立图中所示的机械旋转运动的状态空间表 达式。设转动惯量为达式。设转动惯量为J。 BKTB:粘性阻尼系数。粘性阻尼系数。 K: 扭转轴的刚性系数。扭转轴的刚性系数。T:施加于扭转轴上的力矩。施加于扭转轴上的力矩。 :转动的角度。:转动的角度。1 1. .3 3 .4 .4 由
30、系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式已知系统的传递函数已知系统的传递函数U(S)Y(S)=G (S)=b0Sn + b1Sn1+ bn1 S+bn a0Sn + a1Sn1+ an1 S+an G (S)=b1Sn1+ b2Sn2+ + bn1 S+bnSn + a1Sn1+ an1 S+an+d=G(S)+d化为真分式:化为真分式: 输出与输入之间的直接传递关系输出与输入之间的直接传递关系首先讨论首先讨论G(S)1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表
31、达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式b1Sn1+ b2Sn2+ + bn1 S+bnSn + a1Sn1+ an1 S+anG(S) =1、G(S)特征方程的特征方程的n个极点互异个极点互异用用部分分式法部分分式法G(S) =k1SS1+ +k2SS2+knSSnS1、 S2、 Sn:特征方程的极点特征方程的极点k1、 k2、 kn:待定系数待定系数ki=Lim (SSi) G(S)SSiLim (SSi) SSik1SS1+ +k2SS2+knSSn因为因为 ki=1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式
32、由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式设第设第i个状态变量的拉氏变换为个状态变量的拉氏变换为xi(S)= 1SSiU(S)(SSi) xi(S)= U(S) Sxi (S)=Sixi(S)+U(S) 由拉氏反由拉氏反 变换得状态方程:变换得状态方程:xi(t)= Sixi(t)+ux1(t)= S1x1(t)+ux2(t)= S2x2(t)+uxi(t)= Sixi(t)+uxn(t)= Snxn(t)+u1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式求输
33、出方程:求输出方程:G(S) =k1SS1+ +k2SS2+knSSn=k1x1(S)+ k2x2(S)+ + knxn(S)y(t)=k1x1(t)+ k2x2(t)+ + knxn(t)y(t)=k1x1(t)+ k2x2(t)+ + knxn(t)+du计入计入d的影响的影响Y(S) =k1SS1+ +U(S) +k2SS2U(S) knSSnU(S) 1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式矩阵形式:矩阵形式:y(t)=k1x1(t)+ k2x2(t)+ + knxn
34、(t)+dux1(t)= S1x1(t)+ux2(t)= S2x2(t)+uxi(t)= Sixi(t)+u=S1 0 0 0x1x2xnu1x1x2+x1x21xn110 S2 0 00 0 0 Sny =k1k2 knx1x2x1x2xn对角线标准形对角线标准形+du1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式信号流图:信号流图: x1x1k1S1 xnxnknSny(t)=k1x1(t)+ k2x2(t)+ + knxn(t)+duxi(t)= Sixi(t)+u x2x2
35、k2S2u111111yd例:已知系统传递函数为:例:已知系统传递函数为:G(S) = 2S+1S3+7S2+14S+8试用部分分式法写出状态空间表达式。试用部分分式法写出状态空间表达式。解:由解:由 S3+7S2+14S+8=0求得:求得:S1= 1、 S2= 2、 S3= 4则则 G(S) =k1S+1+k2S+2+k3S+4k1 = Lim (S+1) G(S)= Lim (S+1) S 1S 1 2S+1(S+1)(S+2)(S+4)=1 3k2= Lim (S+2) G(S)= Lim (S+2) S 2S 2 2S+1(S+1)(S+2)(S+4)=3 2k3= Lim (S+4)
36、 G(S)= Lim (S+4) S 4S 4 2S+1(S+1)(S+2)(S+4)=7 6x1x2x3=u1 0 0+x1x2x31110 2 00 04y=x1x2x31 3 3 27 6 x1x11 x3x34 x2x22u111111y1 3 3 26 71 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式2、G(S)特征方程有相重极点特征方程有相重极点设设系统有系统有5个特征根:个特征根:S1、 S1、 S1、 S4、 S5。G(S) =k11(SS1)3+k12(SS1)2
37、+k5SS5k13(SS1)+k4SS4+Kij= Lim 1(j1)!SS1dj1 G(S) (SS1)m dSj1 重极点系数:重极点系数:单极点系数:单极点系数:ki=Lim (SSi) G(S)SSim:重极点的个数重极点的个数1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式2、G(S)特征方程有相重极点特征方程有相重极点x1(S)= 1(SS1)3U(S)设设状态变量的拉氏变换为状态变量的拉氏变换为x2(S)= 1(SS1)2U(S)x3(S)= 1(SS1)U(S)x4(
38、S)= 1(SS4)U(S)x5(S)= 1(SS5)U(S)则:则: x1(S)= 1(SS1)x2(S)x2(S)= 1(SS1)x3(S)x3(S)= 1(SS1)U(S)1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式2、G(S)特征方程有相重极点特征方程有相重极点x4(S)= 1(SS4)U(S)x5(S)= 1(SS5)U(S)x1(S)= 1(SS1)x2(S)x2(S)= 1(SS1)x3(S)x3(S)= 1(SS1)U(S)整理后得:整理后得:Sx1(S)=S1x
39、1(S)+x2(S) Sx2(S)=S1x2(S)+x3(S) Sx3(S)=S1x3(S)+U(S) Sx4(S)=S4x4(S)+U(S) Sx5(S)=S5x5(S)+U(S) 1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式2、G(S)特征方程有相重极点特征方程有相重极点取拉氏反取拉氏反 变换,得系统状态方程:变换,得系统状态方程:Sx1(S)=S1x1(S)+x2(S) Sx2(S)=S1x2(S)+x3(S) Sx3(S)=S1x3(S)+U(S) Sx4(S)=S4x4
40、(S)+U(S) Sx5(S)=S5x5(S)+U(S) x1(t)= S1x1(t)+ x2(t)x4(t)= S4x4(t)+ux2(t)= S1x2(t)+ x3(t)x3(t)= S1x3(t)+ ux5(t)= S5x5(t)+u1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式2、G(S)特征方程有相重极点特征方程有相重极点x1(t)= S1x1(t)+ x2(t)x4(t)= S4x4(t)+ux2(t)= S1x2(t)+ x3(t)x3(t)= S1x3(t)+ ux
41、5(t)= S5x5(t)+u=S1 1 0 0 0x1x2x3x4x5u+001110 S1 1 0 00 0 0 0 S50 0 S1 0 00 0 0 S4 0x1x2x3x4x5系统状态方程:系统状态方程:约当约当标准形标准形2、G(S)特征方程有相重极点特征方程有相重极点Y(S) =+k11(SS1)3U(S) +k5SS5U(S) k13(SS1)U(S) +k12(SS1)2U(S) +k4SS4U(S) =k11x1(S)+ k12x2(S)+ k13x3(S) + k4x4(S)+ k5x5(S)x2(S)= 1(SS1)2U(S)x3(S)= 1(SS1)U(S)x4(S)
42、= 1(SS4)U(S)x5(S)= 1(SS5)U(S)x1(S)= 1(SS1)3U(S)求求输出方程:输出方程:1 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式2、G(S)特征方程有相重极点特征方程有相重极点Y(S)=k11x1(S)+ k12x2(S)+ k13x3(S) + k4x4(S)+ k5x5(S)y(t) =k11x1(t)+ k12x2(t)+ k13x3(t) + k4x4(t)+ k5x5(t)系统输出方程:系统输出方程:y =k11 k12 k13 k4
43、k5x1x2x3x4x51 1. .3 3 .4 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式y(t) =k11x1(t)+ k12x2(t)+ k13x2(t) + k4x4(t)+ k5x5(t)x1(t)= S1x1(t)+ x2(t)x4(t)= S4x4(t)+ux2(t)= S1x2(t)+ x3(t)x3(t)= S1x3(t)+ ux5(t)= S5x5(t)+u信号流图:信号流图: x5x5S5 x4x4S4uk4111 x3x3S11k5y x2x2S11 x1x1S1k11k12k1
44、3例:已知系统传递函数为:例:已知系统传递函数为:G(S) = 4S2+17S+16S3+7S2+16S+12试用部分分式法写出状态空间表达式。试用部分分式法写出状态空间表达式。解:由解:由 S3+7S2+16S+12=0求得:求得:S1= 2、 S2= 2、 S3= 3则则 G(S) =k11(S+2)2+k12S+2+k3S+3(S+2)2(S+3)=0Kij= Lim 1(j1)!SS1dj1 G(S) (SS1)m dSj1 G(S) =k11(S+2)2+k12S+2+k3S+3S2K11= Lim 1(11)!G(S) (S+2)2 Kij= Lim 1(j1)!SS1dj1 G(
45、S) (SS1)m dSj1 K12= Lim 1(21)!S 2d21 G(S) (S+2)2 dS21 1(21)!S 2d 4S2+17S+16dS (S+3) K12= Lim =3= LimS2(S+2)2(S+3)4S2+17S+16(S+2)2=2k3=Lim (S+3) G(S)S3=Lim 4S2+17S+16 (S+2)2 S3=1x1x2x3=u2 1 0+x1x2x30110 2 00 03(S+2)2(S+3)4S2+17S+16G(S)=y=x1x2x32 3 1信号流图或结构图?信号流图或结构图?作业作业1-4 已知系统传递函数为:已知系统传递函数为: S2+S+
46、1S3+6S2+11S+6试用部分分式法写出状态空间表达式,画出系统的试用部分分式法写出状态空间表达式,画出系统的模拟结构图或信号流图。模拟结构图或信号流图。1.3 .1 1.3 .1 由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式由系统方框图建立状态空间表达式1.3 .2 1.3 .2 由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式由系统的工作机理建立状态空间表达式1.3 .3 1.3 .3 由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式由系统的微分方程建立状态空间表达式
47、由系统的微分方程建立状态空间表达式1.3 .4 1.3 .4 由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式由系统传递函数建立状态空间表达式1.4 状态方程的线性变换状态方程的线性变换 特征矢量线性变换法,把状态方程化为对角线标准形特征矢量线性变换法,把状态方程化为对角线标准形或约当标准形。或约当标准形。1 1. .4 4 .1 .1 系统状态的线性变换系统状态的线性变换系统状态的线性变换系统状态的线性变换设:设:X=x1 x2 x3 xnTX=x1 x2 x3 xnTX=x1 x2 x3 xnT它们之间的线性变换:它们之间的线性变换:X=P1
48、XX=PXP:nn非奇异变换阵非奇异变换阵已知线性系统为:已知线性系统为:X=AX+BUY=CX+DUX=PX令:令:X=PX代入状态方程代入状态方程Y=CPX+DUPX=APX+BUX=P1APX+ P1BUY=CPX+DU1 1. .4 4 .1 .1 系统状态的线性变换系统状态的线性变换系统状态的线性变换系统状态的线性变换X=P1APX+ P1BUY=CPX+DUPX=APX+BUY=CPX+DUY=CX+DUX=AX+BU 式中式中A= P1AP B= P1B C=CP D=D1 1. .4 4 .2 .2 系统的特征值系统的特征值系统的特征值系统的特征值1、特征值及特征矢量、特征值及
49、特征矢量AP= PA的特征值的特征值 P AP=0( I A)P =0有非零解的必要条件:有非零解的必要条件:| I A | =0A的特征方程的特征方程1 1. .4 4 .2 .2 系统的特征值系统的特征值系统的特征值系统的特征值| I A | = n +a1 n1+a n1 + a n求得求得A的特征值:的特征值: 1、 2、 n Api= ipi对应于对应于 i的一个特征矢量的一个特征矢量由由全部全部 所所对应的特征矢量:对应的特征矢量:P=p1 p2 pi pn =p11p12. p1n p21p22. P2n . . . . . Pn1Pn2. pnn 线性变换阵线性变换阵2、特征值
50、不变性、特征值不变性特征多项式特征多项式经线性变换后,其特征值不变。经线性变换后,其特征值不变。| I A |2、特征值不变性、特征值不变性X=AX+BU对于对于X=AX+BU对于对于| I A |A= P1AP B= P1B| I A |=| I A |证明:证明:| I A |=| I P1AP| =| P1P P1AP| |P1( I A) P|= |P1| I A| P |=|P1| P| I A|所以:所以: = 1 1. .4 4 .3 .3 化状态方程为对角线标准形化状态方程为对角线标准形化状态方程为对角线标准形化状态方程为对角线标准形设设X=AX+BU,若,若A的特征值的特征值
51、 1、 2、 n互异,互异,则必则必存在非奇异变换阵存在非奇异变换阵P,使其进行使其进行X=PX的变换后,其的变换后,其状态方程状态方程X=AX+BU将为将为对角线标准形,即对角线标准形,即 1 0 0 00 2 0 00 0 0 nA=且且P=p1 p2 pi pn =p11p12. p1n p21p22. P2n . . . . . Pn1Pn2. pnn A= P1APB= P1BC=CPD=DX=P1X设:设:A的特征值:的特征值: 1、 2、 n 特征矢量:特征矢量:p1 p2 pi pn证明证明A= P1AP 1 0 0 00 2 0 00 0 0 nA=证明:证明:X=P1APX
52、+ P1BUX=AX+BU若若Pi是对应于是对应于 i的一个特征矢量的一个特征矢量则必则必满足满足( iI A)pi =0Ap2= 2p2Api= ipiAp1= 1p1证明:证明:Ap2= 2p2Api= ipiAp1= 1p1Api= ipiApn= npn写成矩阵形式:写成矩阵形式:Ap1 Ap2 Api Apn= 1p1 2p2 ipi npn Ap1 p2 pi pn= 1p1 2p2 ipi npn AP= p1 p2 pi pn 1 0 0 00 2 0 00 0 0 n1 1. .4 4 .3 .3 化状态方程为对角线标准形化状态方程为对角线标准形化状态方程为对角线标准形化状态
53、方程为对角线标准形 1 0 0 00 2 0 00 0 0 n=P证明:证明: AP= p1 p2 pi pn 1 0 0 00 2 0 00 0 0 nA= P1AP 1 0 0 00 2 0 00 0 0 n=证毕证毕例:试将状态方程例:试将状态方程变换为对角线标准形。变换为对角线标准形。x1x2x3=u0 1 1+x1x2x3001 6 11 6 6 11 5 解:解:(1)求系统特征值求系统特征值| I A | =0 1 1 6 11 6 6 11 5 0 0 0 0 0 0 = 1 1 6 +11 6 6 11 5 | I A |= 1 1 6 +11 6 6 11 5 = 3+6
54、2+11 +6=0( +1)( +2) ( +3)=0(2)求特征矢量求特征矢量对应于对应于 1= 1 的特征矢量为的特征矢量为p1,Ap1= 1p10 1 1 6 11 6 6 11 5p11p21p31= p11 p21 p31p21p31= p11 6p11 11 p21+6p31= p21 6p11 11 p21+5p31= p31p21p31= p11 6p11 10 p21+6p31= 0 6p11 11 p21+6p31= 0可以看出:可以看出: p21=0 p11=p31令令 p11=1 p31 =1求得:求得: 1= 1、 2= 2、 3= 3同理,将同理,将 2= 2、 3
55、= 3分别代入分别代入Ap2= 2p2Ap3= 3p3求得:求得: 1 0 1=p11p21p31p1= 1 2 4=p12p22p32p2= 1 6 9=p13p23p33p3=变换阵变换阵P=p1 p2 p3=3 5/2 1 3 4 31 3/2 1P1=1 1 10 2 61 4 9A= P1AP则则= 1 0 0 0 2 0 0 0 3B= P1B = 2 3 1A= P1AP则则= 1 0 0 0 2 0 0 0 3B= P1B = 2 3 1=+x1x2x3x1x2x3 1 0 0 0 2 0 0 0 3 2 3 1u1 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程
56、为约当标准形化状态方程为约当标准形化状态方程为约当标准形当当A有相重特征值时;有相重特征值时;1、A的线性独立特征矢量数的线性独立特征矢量数等于等于它的阶数它的阶数n,这时这时A仍可以化为仍可以化为对角线标准形对角线标准形;2、A的线性独立特征矢量数的线性独立特征矢量数小于小于它的阶数它的阶数n,这时这时A不能化为对角线标准形,只能化为不能化为对角线标准形,只能化为约当标准形约当标准形;1 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形例:例: A=对应特征值:对应特征值: 1= 1、 2= 1、 3= 2对应于对应于 1=
57、 1 的特征矢量为的特征矢量为p1,0 0 10 0 00 0 1 1I A p1 =p11p21p31= 01 0 10 1 00 0 2显然显然 1I A 的秩是的秩是1, p1有有两个独立的解,对应两个两个独立的解,对应两个独立的特征矢量,即独立的特征矢量,即A= P1AP =1 0 0 0 1 0 0 0 21 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形 1 0 0p1= 1 0 1p3=对应于对应于 3= 3 的特征矢量为的特征矢量为p3,则则 P=1 0 1 0 1 0 0 0 1 0 1 0p2=0 0 1
58、0 0 00 0 1 1I A p1 =p11p21p31= 01 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形若:若: A=对应特征值:对应特征值: 1= 1、 2= 1、 3= 21 1 20 1 30 0 1但但rank 1I A =2,独立的特征矢量只有一个。独立的特征矢量只有一个。约当标准形:约当标准形:Y=CX+DUX=JX+BU 式中式中J= Q1AQ B= Q1B C=CQ D=D1 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形设设A是
59、是55的方阵,其特征值为的方阵,其特征值为 1、 1、 1、 4和和 5,存在一个变换阵存在一个变换阵Q,使得使得J= Q1AQA1 0 0 0 A2 0 0 0 A3=A的约当标准形的约当标准形J由三个约当块组成。若由三个约当块组成。若 1 只有一个只有一个独立的特征矢量,则独立的特征矢量,则 1 1 0 0 00 1 1 0 00 0 0 0 50 0 1 0 00 0 0 4 0J= Q1AQ=1 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形J= Q1AQ= 1 1 0 0 00 1 0 0 00 0 0 0 50
60、0 1 0 00 0 0 4 0若若 1 有两个独立的特征矢量,则有两个独立的特征矢量,则若若 1 有三个独立的特征矢量,则有三个独立的特征矢量,则J= Q1AQ= 1 0 0 0 00 1 0 0 00 0 0 0 50 0 1 0 00 0 0 4 01 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形设设 1 有一个独立的特征矢量,求有一个独立的特征矢量,求Q。J= Q1AQQ J=AQ令令Q=q11q12. q15q21q22. q25. . . . . q51q52. q55=Q1 Q2 Q3 Q4 Q5 1 1 0
61、 0 00 1 1 0 00 0 0 0 50 0 1 0 00 0 0 4 0Q1 Q2 Q3 Q4 Q5 =A Q1 Q2 Q3 Q4 Q51 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形 1 1 0 0 00 1 1 0 00 0 0 0 50 0 1 0 00 0 0 4 0Q1 Q2 Q3 Q4 Q5 =A Q1 Q2 Q3 Q4 Q5将将上式展开得:上式展开得: 1 Q1=A Q1Q1 + 1 Q2=A Q2Q2 + 1 Q3=A Q3 5Q5=A Q5 4 Q4=A Q4( 1 IA)Q1= 0( 1 IA)
62、Q2= Q1( 1 IA)Q3= Q2( 4 IA)Q4= 0( 5 IA)Q5= 0Q1、Q4、Q5为独立特征矢量,为独立特征矢量,Q2、Q3为非独立特征矢量为非独立特征矢量1 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形( 1 IA)Q1= 0( 1 IA)Q2= Q1( 1 IA)Q3= Q2( 4 IA)Q4= 0( 5 IA)Q5= 0解此解此解此解此方程组得变换阵方程组得变换阵方程组得变换阵方程组得变换阵QQ例:例: A=0 6 51 0 23 2 4化化A为标准形。为标准形。 解:解:(1)求系统特征值求系统
63、特征值| I A | =0 6 5 1 0 2 3 2 4 0 0 0 0 0 0 = 6 5 1 2 3 2 4 = ( 1)2( 2) 1 1. .4 4 .4 .4 化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形化状态方程为约当标准形 | I A | =( 1)2( 2) 求得:求得: 1= 1、 2= 1、 3= 2将将 1= 1 代入代入( 1 IA)Q1= 0中中 6 5 1 2 3 2 4 q11q21q311 6 5 1 1 2 3 2 3 =0rank 1I A =2,独立的特征矢量只有一个。独立的特征矢量只有一个。任取任取q11=1解得解得q21= 3/
64、7, q31= 5/7再将再将Q1代入代入( 1 IA)Q2= Q1中中 1 3/7 5/7Q1=rank 2I A =2,独立的特征矢量只有一个。独立的特征矢量只有一个。任取任取q12=1解得解得q22= 22/49, q32= 46/49再将再将Q1代入代入( 1 IA)Q2= Q1中中要要保证保证Q阵非奇异阵非奇异将将 3=2代入代入( 3 IA)Q3= 0中中2 6 5 1 2 2 3 2 2 q13q23q33=0rank 2I A =2,独立的特征矢量只有一个。独立的特征矢量只有一个。令令q13=2,则,则q23= 1,q33= 2 1 22/49 46/49Q2=q12q22q3
65、2=13/75/71 6 5 1 1 2 3 2 3 2 1 2Q3= 1 22/49 46/49Q2= 1 3/7 5/7Q1=Q=1 1 2 3/7 22/49 1 5/7 46/49 21 1 00 1 00 0 2J= Q1AQ=作业作业1-5 已知已知Y= 1 0x1x2x1x2x1x2=+u11 0 1 2 3试试化为标准形并求其传递函数。化为标准形并求其传递函数。1.5 系统的传递函数阵系统的传递函数阵1 1. .5 5 .1 .1 传递函数阵的概念传递函数阵的概念传递函数阵的概念传递函数阵的概念G11(S)G21(S)G12 (S)G22 (S)+ +Y2(S)U1(S)U2(
66、S)Y1(S)当初始条件为零时当初始条件为零时Y1(S)=G11(S) U1(S)+ G12(S) U2(S)Y2(S)=G21(S) U1(S)+ G22(S) U2(S)双输入双输出双输入双输出1 1. .5 5 .1 .1 传递函数阵的概念传递函数阵的概念传递函数阵的概念传递函数阵的概念Y1(S)=G11(S) U1(S)+ G12(S) U2(S)Y2(S)=G21(S) U1(S)+ G22(S) U2(S)Y1(S)Y2(S)G11(S) G12(S)G21(S) G22(S) U1(S)U2(S)=简记为:简记为:Y (S) =G(S) U (S)G11G12G1nG21G22
67、G2n Gn1Gn2 Gnnn个输入个输入n个输出个输出G11 0 0 00 G22 0 0 0 0 0Gnn分离式控制分离式控制分离式控制分离式控制1 1. .5 5 .2 .2 闭环系统的传递函数阵闭环系统的传递函数阵闭环系统的传递函数阵闭环系统的传递函数阵Y (S)=G0(S) E(S)= G0(S) U(S)F(S)G0(S)H(S)E(S)F(S)U(S)Y(S)+=G0(S) U(S)H(S) Y(S) I+G0(S) H(S) Y(S)= G0(S) U(S)Y (S)= I+G0(S) H(S) 1 G0(S) U(S)闭环传递函数阵:闭环传递函数阵:G(S)= I+G0(S)
68、 H(S) 1 G0(S) 1 1. .5 5 .3 .3 由状态空间表达式求传递函数阵由状态空间表达式求传递函数阵由状态空间表达式求传递函数阵由状态空间表达式求传递函数阵X=AX+BUY=CX+DUX:n维,维,Y:m维,维,U:r维,维,对对上式取拉氏变换上式取拉氏变换SX(S)X(0)=AX(S)+BU(S)X(0)=0Y(S)=CX(S)+DU(S)SX(S) =AX(S)+BU(S)X(S) =SIA1BU(S)Y(S)= CSIA1B +D U(S)传递函数阵:传递函数阵:G(S)= CSIA1B +D例:已知系统的状态空间表达式为例:已知系统的状态空间表达式为x1=x2+u1x2
69、=x3+2u1u2x3=6x111x26x3+2u2y1=x1 x2y2 =2x1+x2 x3 试求其传递函数阵。试求其传递函数阵。解:解:A=0 1 00 0 16 11 61 02 10 2B=C=11 02 1 1 G(S)= CSIA1B +D11 02 1 1 S 1 0 0 S 1 6 11 S+ 6 11 02 10 2G(S)=S3+6S2+11S+61S2 4S+29S2 + 3S 44S2 + 56S+523S2 17S 14=1 14 3S2+4 356 1729 452 14S+S3+6S2+11S+61 1. .5 5 .4 .4 传递函数阵的不变性传递函数阵的不变性
70、传递函数阵的不变性传递函数阵的不变性 对对状态方程进行线性变换后,其对应的传递函数阵不变状态方程进行线性变换后,其对应的传递函数阵不变.证明:证明:G(S):原原系统的传递函数阵系统的传递函数阵 G(S):经过线性变换的状态空间表达式所对应的经过线性变换的状态空间表达式所对应的 传递函数阵。传递函数阵。G(S)= CSIA1B +DG(S)= CSIA1B +D设设P是非奇异变换阵,由于是非奇异变换阵,由于A= P1APB= P1BC=CPD=D= CPS P1PP1AP1 P1B +D= CPP1( S IA) P1 P1B +D= CP P1 S IA 1 PP1B +D= C S IA
71、1 B +D故故G(S)= CPSIP1AP1 P1B +D1 1. .5 5 .5 .5 组合系统的状态空间表达式与组合系统的状态空间表达式与组合系统的状态空间表达式与组合系统的状态空间表达式与 传递函数阵传递函数阵传递函数阵传递函数阵1 1、子系统的并联联结、子系统的并联联结、子系统的并联联结、子系统的并联联结设设子系统子系统S1、S2分别为分别为n1、n2维,其组合系统的示意图维,其组合系统的示意图 B1D1C1A1U(t)+Y1(t)X1X1 B2D2C2A2+Y2(t)X2X2+Y (t)X1=A1X1+B1UY1=C1X1+D1UX2=A2X2+B2UY2=C2X2+D2UX1=A
72、1X1+B1UY1=C1X1+D1UX2=A2X2+B2UY2=C2X2+D2U组合系统组合系统Y=Y1+ Y2= C1X1+D1U+ C2X2+D2UY= C1 C2D1+ D2X1X2+U传递函数阵:传递函数阵:G1(S)= C1SIA11B1 +D1G2(S)= C2SIA21B2 +D2G(S)= CSIA1B +DX1X2X1X2=+UB1B2A1 0 0 A21 1、子系统的并联联结、子系统的并联联结、子系统的并联联结、子系统的并联联结传递函数阵:传递函数阵:G1(S)= C1SIA11B1 +D1G2(S)= C2SIA21B2 +D2G(S)= CSIA1B +DB1B2SIA
73、1 0 0 SI A2C1 C21+ D1+ D2=B1B2(SIA1) 0 0 (SI A2)C1 C21+ D1+ D2=1= C1SIA11B1 + C2SIA21B2 + D1+ D2= C1SIA11B1 +D1+ C2SIA21B2 +D2G(S) = G1(S)+ G2(S) 1 1、子系统的并联联结、子系统的并联联结、子系统的并联联结、子系统的并联联结 1 1. .5 5 .5 .5 组合系统的状态空间表达式与组合系统的状态空间表达式与组合系统的状态空间表达式与组合系统的状态空间表达式与 传递函数阵传递函数阵传递函数阵传递函数阵2 2、子系统的串联联结、子系统的串联联结、子系统
74、的串联联结、子系统的串联联结X1=A1X1+B1UY1=C1X1+D1U=U2X2=A2X2+B2U2 =A2X2+B2(C1X1+D1U) =A2X2+B2C1X1+ B2 D1UY2=C2X2+D2(C1X1+D1U)= C2X2+D2C1X1+D2D1U B1D1C1A1U(t)+Y1(t)X1X1 B2D2C2A2+Y2(t)X2X2Y (t)U2(t)1 1. .5 5 .5 .5 组合系统的状态空间表达式与传递函数阵组合系统的状态空间表达式与传递函数阵组合系统的状态空间表达式与传递函数阵组合系统的状态空间表达式与传递函数阵2 2、子系统的串联联结、子系统的串联联结、子系统的串联联结
75、、子系统的串联联结X1=A1X1+B1UY1=C1X1+D1U=U2X2=A2X2+B2U2 =A2X2+B2(C1X1+D1U) =A2X2+B2C1X1+ B2 D1UY2=C2X2+D2(C1X1+D1U)= C2X2+D2C1X1+D2D1UY= D2C1 C2D2 D1 UX1X2+X1X2X1X2=+UB1B2D1A1 0B2C1 A22 2、子系统的串联联结、子系统的串联联结、子系统的串联联结、子系统的串联联结X1=A1X1+B1UY1=C1X1+D1U=U2X2=A2X2+B2C1X1+ B2 D1UY2=C2X2+D2C1X1+D2D1UY= D2C1 C2D2 D1 UX1
76、X2+X1X2X1X2=+UB1B2D1A1 0B2C1 A2传递函数阵:传递函数阵: G(S)= CSIA1B +DB1B2D1D2C1 C2SIA1 0 B2C1 SI A21+ D2 D1=2 2、子系统的串联联结、子系统的串联联结、子系统的串联联结、子系统的串联联结SIA1 0 B2C1 SI A21=(SIA1)1(SIA2)1 B2C1(SIA1)1(SIA2)10G(S)= C2SIA21B2 + D2 C1SIA11B1 + D1G(S)= G 2(S) G1(S) G(S)B1B2D1D2C1 C2SIA1 0 B2C1 SI A21+ D2 D1=3 3、子系统的反馈联结、
77、子系统的反馈联结、子系统的反馈联结、子系统的反馈联结Y (t) B1C1A1U(t)+Y1(t)X1X1 B2C2A2+Y2(t)X2X2U2(t)U1(t)X1=A1X1+B1U1=A1X1+B1(UC2X2)= A1X1B1C2X2 +B1UY1=C1X1=YX2=B2U2+A2X2= B2 C1X1+ A2X2Y2=C2X23 3、子系统的反馈联结、子系统的反馈联结、子系统的反馈联结、子系统的反馈联结X1=A1X1+B1U1=A1X1+B1(UC2X2)= A1X1B1C2X2 +B1UY1=C1X1=YX2=B2U2+A2X2= B2 C1X1+ A2X2Y= C1 0X1X2X1X2
78、X1X2=+UB1 0A1 B1C2B2C1 A2G1(S)= C1SIA11B1 =G0(S)G2(S)= C2SIA21B2=H(S)G(S)= CSIA1B传递函数阵:传递函数阵:G(S)=I+ G0(S) H(S)1G0(S)y2= 2 1x3x4x1x2x1x2=+u01 0 1 2 3作业作业1-6 设子系统设子系统S1、S2的状态空间表达式分别为的状态空间表达式分别为x3x4x3x4=+u01 0 1 12 7y1= 1 0x1x2求求并联系统的状态空间表达式及它的传递函数。并联系统的状态空间表达式及它的传递函数。1.6 离散系统的状态空间表达式离散系统的状态空间表达式X(k+1
79、)T=G(kT) X(kT)+ H(kT)u(kT) (1)Y (kT)=C(kT) X(kT)+ D(kT)u(kT) (2)T:采样周期。采样周期。 (1)式是一阶差分方程组式是一阶差分方程组(2)式是代数方程组式是代数方程组离散系统的状态空间表达式由系统的高阶差分方程或离散系统的状态空间表达式由系统的高阶差分方程或脉冲传递函数获得。脉冲传递函数获得。1.7 时变系统和非线性系统的状时变系统和非线性系统的状 态空间表达式态空间表达式一、线性时变系统一、线性时变系统X=A(t)X+B (t) UY=C (t) X+D (t) UA(t)=a11 (t) a12 (t). a1n (t) a2
80、1 (t) a22 (t). a2n (t) . . . . . an1 (t) an2 (t). ann (t) B(t)=b11 (t) b12 (t). b1r (t) b21 (t) b22 (t). b2r (t) . . . . . bn1 (t) bn2 (t). bnr (t) C(t)=c11 (t) c12 (t). c1n (t) c21 (t) c22 (t). c2n (t) . . . . . cm1 (t) cm2 (t). cmn (t) D(t)=d11 (t) d12 (t). d1r (t) d21 (t) d22 (t). d2r (t) . . . . . dm1 (t) dm2 (t). dmr (t) 二、非线性系统二、非线性系统X=F(X U t)Y=G(X U t)