《线性系统理论讲义PPT课件》由会员分享,可在线阅读,更多相关《线性系统理论讲义PPT课件(300页珍藏版)》请在金锄头文库上搜索。
1、第一部分第一部分: 线性系统时间域理论线性系统时间域理论 第二章第二章 线性系统的状态空间描述线性系统的状态空间描述 2.1 状态和状态空间 线性系统时间域理论是以时间域数学模型为系统描述线性系统时间域理论是以时间域数学模型为系统描述,直直接在时间域内分析和综合线性系统的运动和特性的一种理论接在时间域内分析和综合线性系统的运动和特性的一种理论和方法。和方法。 系统动态过程的数学描述 1/4,1/501/4,1/501(1).系统的外部描述 外部描述常被称作为输出输入描述例如.对SISO线性定常系统:时间域的外部描述:复频率域描述即传递函数描述 (2)系统的内部描述 状态空间描述是系统内部描述的
2、基本形式,需要由两个数学方程表征, 状态方程和输出方程 (3)外部描述和内部描述的比较 一般的说外部描述只是对系统的一种不完全描述,不能反映黑箱内部结构的不能控或不能观测的部分. 内部描述则是系统的一种完全的描述,能够完全反映系统的所有动力学特性. 2/4,2/502/4,2/5022.1 基本概念2.1.1 2.1.1 定义定义(1)(1)状态:状态:系统过去、现在和将来的状况系统过去、现在和将来的状况(2)(2)状态变量:状态变量: 能够完全表征系统运动状态的能够完全表征系统运动状态的最小一组变量:最小一组变量:表示系统表示系统 时刻的状态时刻的状态 当当时的输入时的输入给定,且上述给定,
3、且上述时的行为时的行为 状态确定时,状态变量能完全确定系统状态确定时,状态变量能完全确定系统初始初始在在3(5)(5)状态方程状态方程:描述系统状态与输入之间关系:描述系统状态与输入之间关系的、一阶微分方程(组):的、一阶微分方程(组):作为分量的向量,即作为分量的向量,即(3) (3) 状态向量:以系统的状态向量:以系统的 个独立状态变量个独立状态变量为为(4) (4) 状态空间状态空间: : 以状态变量以状态变量 维空间。维空间。标轴构成的标轴构成的坐坐4输出方程:描述系统输出与状态、输入之间关输出方程:描述系统输出与状态、输入之间关系的数学表达式:系的数学表达式:(6)(6)(7)(7)
4、状态空间表达式:状态空间表达式: (5)+ (6). (5)+ (6).状态变量的特点:状态变量的特点:(1)(1)独立性:状态变量之间线性独立独立性:状态变量之间线性独立. . (2) (2)多样性:状态变量的选取并不唯一,实多样性:状态变量的选取并不唯一,实际上存在无穷多种方案际上存在无穷多种方案. .(3)等价性:两个状态向量之间只差一个非奇异等价性:两个状态向量之间只差一个非奇异变换变换. .5(4)(4)现实性:状态变量通常取为涵义明确的物理量现实性:状态变量通常取为涵义明确的物理量. .(5)(5)抽象性:状态变量可以没有直观的物理意义抽象性:状态变量可以没有直观的物理意义. .2
5、.1.2 2.1.2 状态空间表达式的一般形式:状态空间表达式的一般形式:(1)(1)线性系统线性系统62.2 线性系统的状态空间描述 电路系统状态空间描述的列写示例 以上方程可表为形如 描述系统输入、输出和状态变量之间关系的方程组称为系统的状态空间表达式(动态方程或运动方程),包括状态方程(描述输入和状态变量之间的关系)和输出方程(描述输出和输入、状态变量之间的关系)。 1/7,5/501/7,5/507机电系统状态空间描述的列写示例 上式可表为形如 2/7,6/502/7,6/508连续时间线性系统的状态空间描述 动态系统的结构连续时间线性系统的状态空间描述 线性时不变系统 线性时变系统
6、3/7,7/503/7,7/509连续时间线性系统的方块图 4/7,8/504/7,8/5010人口分布问题状态空间描述的列写示例假设某个国家,城市人口为107,乡村人口为9x107,每年4%的城市人口迁移去乡村, 2%的乡村人口迁移去城市,整个国家的人口的自然增长率为1% 设k为离散时间变量, x1(k)、x2(k)为第k年的城市人口和乡村人口, u(k)为第k年所采取的激励性政策控制手段,设一个单位正控制措施可激励5x104城市人口迁移乡村,而一个单位负控制措施会导致5x104乡村人口去城市, y(k)为第k年全国人口数 写成矩阵形式5/7,9/505/7,9/5011离散时间线性系统的状
7、态空间描述状态空间描述形式离散时间线性时不变系统 离散时间线性时变系统6/7,10/506/7,10/5012离散系统状态空间描述的特点:离散系统状态空间描述的特点:一是:状态方程形式上的差分型属性(即:状态方程为差分方程。)二是:描述方程的线性属性。(状态方程和输出方程的右端,对状态x和输入u都 呈现为线性关系。)三是:变量取值时间的离散属性(所有变量只能在离散时刻k取值)。离散时间线性系统的方块图7/7,11/507/7,11/50132.3.连续变量动态系统按状态空间描述的分类 线性系统和非线性系统 设系统的状态空间描述为 向量函数 若f(x,u,t),g(x,u,t)的全部或至少一个组
8、成元素为x、u的非线性函数,该系统称为非线性系统 。若f(x,u,t),g(x,u,t)的全部组成元为x、u的线性函数,该系统称为线性系统 对于线性系统 非线性系统可以用泰勒展开方法化为线性系统。 1/2,12/501/2,12/5014时变系统和时不变系统 若向量f,g不显含时间变量t,即 该系统称为时不变系统 若向量f,g显含时间变量t,即 该系统称为时变系统 连续时间系统和离散时间系统 当且仅当系统的输入变量,状态变量和输出变量取值于连续时间点,反映变量间因果关系的动态过程为时间的连续过程,该系统称为连续时间系统 当且仅当系统的输入变量,状态变量和输出变量只取值于离散时间点,反映变量间因
9、果关系的动态过程为时间的不连续过程,该系统称为离散时间系统.确定性系统和不确定性系统 称一个系统为确定性系统,当且仅当不论是系统的特性和参数还是系统的输入和扰动,都是随时间按确定的规律而变化的. 称一个动态系统为不确定性系统,或者系统的特性和参数中包含某种不确定性,或者作用于系统的输入和扰动是随机变量 2/2,13/502/2,13/50152.4 由系统输入输出描述导出状态空间描述 由输入输出描述导出状态空间描述 对于单输入,单输出线性时不变系统,其微分方程描述 其传递函数描述 可以导出其状态空间描述为 1/18,14/501/18,14/5016结论1 给定单输入,单输出线性时不变系统的输
10、入输出描述,其对应的状态空间描述可按如下两类情况导出 (1)m=n,即系统为真情形证明:设 2/18,15/502/18,15/5017可见 3/18,16/503/18,16/5018令 有 4/18,17/504/18,17/5019(2)mn,即系统为严真情形 写成矩阵形式: 5/18,18/505/18,18/5020结论2 给定单输入,单输出线性时不变系统的输入输出描述,其对应的状态空间描述可按如下两类情况导出 (1)m=0情形此时输入输出描述为: 选取n个状态变量 6/18,19/506/18,19/5021其对应的状态空间描述为: 7/18,20/507/18,20/5022例例
11、2.2 2.2 考虑系统考虑系统试写出其能控标准型状态空间表达式试写出其能控标准型状态空间表达式。则状态空间表达式为:则状态空间表达式为:选择状态变量:选择状态变量:23(2)m0情形此时输入输出描述为: a:8/18,21/508/18,21/5024其对应的状态空间描述为: 其中9/18,22/509/18,22/5025b:改写为 令 10/18,23/5010/18,23/5026结论3 给定单输入单输出线性时不变系统的传递函数描述为: 其极点即分母方程的根 为两两互异实数,则对应的状态空间描述可按如下两类情形导出:(1) m t0的有限时间区间t0,t1内,存在一个无约束的控制矢量u
12、(t),使x(t1)=0,则称系统是状态完全能控的,简称系统是能控的。可见系统的能控性反映了控制矢量u(t)对系统状态的控制性质,与系统的内部结构和参数有关。 定义1/3,1/451/3,1/4593能控性,能达性定义 对连续时间线性时变系统 如果存在一个时刻 以及一个无约束的容许控制u(t) 使系统状态由x(t0)=x0转移到x(t1)=0,则称非零状态X0在t0时刻为能控。 如果存在一个时刻t1J,t1t0,以及一个无约束的容许控制u(t),tt0,t1,使系统状态由x(t0)=0转移到x(t1)=xf0,则称非零状态xf在t0时刻为能达。 *对连续时间线性时不变系统,能控性和能达性等价;
13、对离散时间线性时不变系统和线性时变系统,若系统矩阵为非奇异,则能控性和能达性等价;对连续时间线性系统,能控性和能达性一般为不等价。 定义:对连续时间线性时变系统 和指定初始时刻t0J,如果状态空间中所有非零状态在时刻t0J都为能控/能达,称系统在时刻t0为完全能控/能达。 2/3,2/452/3,2/4594定义:对连续时间线性时变系统 和指定初始时刻t0J,如果状态空间中存在一个非零状态或一个非空状态集合在时刻t0J为不能控/能达,称系统在时刻t0为不完全能控/能达。 定义:若系统的能控/能达性与初始时刻t0的选取无关,或系统在任意初始时刻t0J均为完全能控/能达,则称系统为一致完全能控/能
14、达。能观测性定义对连续时间线性时变系统和指定初始时刻t0J,如果存在一个时刻t1J,t1t0,使系统以x(t0)=x0为初始状态的输出y(t)恒为零,即y(t)0,tt0,t1,则称非零状态x0在时刻t0为不能观测;如果状态空间中所有非零状态在时刻t0都不为不能观测,则称系统在时刻t0为完全能观测;如果状态空间中存在一个非零状态或一个非零状态集合在时刻t0为不能观测,则称系统在时刻t0为不完全能观测;如果系统对任意时刻均为完全能观测,即能观测性与初始时刻t0的选取无关,则称系统为一致完全能观测。 该系统是不完全能观测的由于 可见系统的状态x(t)的能观测性与x(t0)的能观测性是等价的。3/3
15、,3/453/3,3/459542 连续时间线性系统的能控性判据 结论1:线性时变系统在t0时刻是状态完全能控的充分必要条件是下列格兰姆矩阵为非奇异矩阵 证明 充分性 为非奇异时,系统能控 说明系统是能控的 1/8,4/451/8,4/4596反证法 必要性 是奇异的,且系统能控,看能否导出矛盾的结果。由于是奇异的,故的行向量在t0,t1上线性相关,必存在非零的行向量,使在t0,t1区间成立,若选择非零的初始状态x(t0)= T,则 说明=0,矛盾 2/8,5/452/8,5/4597结论2:连续时间线性时不变系统: 完全能控的充分必要条件是,存在时刻t10,使格拉姆矩阵 为非奇异。 结论3:
16、n 维连续时间线性时变系统 设A(t),B(t)对t为n-1阶连续可微,定义 则系统在时刻t0J完全能控的一个充分条件为,存在一个有限时刻t1J,t1t0,,使 3/8,6/453/8,6/4598结论4 对n 维连续时间线性时不变系统,系统完全能控的充分必要条件为能控性判别矩阵 满秩,即rankQ c=n 结论5n 维连续时间线性时不变系统完全能控的充分必要条件为:rankSI-AB=n,或为系统特征值结论6:n 维连续时间线性时不变系统完全能控的充分必要条件为:矩阵A不存在与B所有列正交的非零左特征向量,即对矩阵A所有特征值i,使同时满足TA= i T , TB=0 的左特征向量T=0。
17、4/8,7/454/8,7/4599结论7:对n维线性时不变系统,若A为对角阵,且其特征值两两相异,系统完全能控的充分必要条件是B中不包含零行向量。结论8:对n维线性时不变系统,若A为约当阵,系统完全能控的充分必要条件是:特征值互异的约当块最后一行对应的B阵中,该行元素不全为零。 特征值相同的各约当块最后一行对应的B阵各行向量线性无关。5/8,8/455/8,8/45100例例 图示电路,判断系统能控性条件 解解 选取状态变量x1=iL,x2=uC,得系统的状态方程为: 6/8,9/456/8,9/45101(R1R4=R2R3)时,系统不能控。否则系统能控。 例例 系统能控的充分必要条件是向
18、量组bl11、bl12、bl13线性无关以及bl21线性无关(即不为零)。7/8,10/457/8,10/45102定义:令 对完全能控连续时间线性时不变系统,定义能控性指数为:使“rankQk=n”成立的最小正整数k。 结论9:对完全能控单输入连续时间线性时不变系统,状态维数为n,则系统能控性指数n。结论10:对完全能控多输入连续时间线性时不变系统,状态维数为n,输入维数为p,设rankB=r,则能控性指数满足 设 为矩阵A的最小多项式次数,则 结论11:多输入连续时间线性时不变系统,状态维数为n,输入维数为p,且rankB=r,则系统完全能控的充分必要条件为: 8/8,11/458/8,1
19、1/4510343 连续时间线性系统的能观测性判据 结论1:线性时变系统在t0时刻是状态完全能观测的充分必要条件是下列格兰姆矩阵为非奇异矩阵 结论2:连续时间线性时不变系统完全能观测的充分必要条件是,存在时刻t10,使格拉姆矩阵 为非奇异。 1/5,12/451/5,12/45104结论3:n 维连续时间线性时变系统设A(t),C(t)对t为n-1阶连续可微,定义 则系统在时刻t0J完全能观测的一个充分条件为,存在一个有限时刻t1J,t1t0,,使 2/5,13/452/5,13/45105结论4 对n 维连续时间线性时不变系统,系统完全能观测的充分必要条件为能观测性判别矩阵 满秩,即rank
20、Q o=n 结论5n 维连续时间线性时不变系统完全能观测的充分必要条件为:或为系统特征值结论6:n维连续时间线性时不变系统完全能观测的充分必要条件为:矩阵A不存在与C所有行正交的非零右特征向量,即对矩阵A所有特征值,使同时满足的右特征向量 3/5,14/453/5,14/45106结论7:对n维连续时间线性时不变系统,若A为对角阵,且其特征值两两相异,系统完全能观测的充分必要条件是C阵中不包含零列向量。 结论8:对n维连续时间线性时不变系统,若A为约当阵,系统完全能观测的充分必要条件是:特征值互异的约当块第一列对应的C阵中,该列元素不全为零。特征值相同的约当块第一列对应的C阵中,各列向量线性无
21、关。4/5,15/454/5,15/45107定义:令 完全能观测n维连续时间线性时不变系统的能观测性指数定义为 使“rankQk=n”成立的最小正整数。 结论9:对完全能观测单输出连续时间线性时不变系统,状态维数为n,则能观测性指数为n。 结论10:对完全能观测多输出连续时间线性时不变系统,状态维数为n,输入维数为q,设rankC=m,则 设 为矩阵A的最小多项式次数,则 结论11:对多输出连续时间线性时不变系统,设rankC=m,则系统完全能观测的充分必要条件是: 5/5,16/455/5,16/451084.4 离散时间线性系统的能控性和能观性判据 时变系统的能控性和能观性判据定义 离散
22、时间线性时变系统 如果对初始时刻hJk 和任意非零初始状态X(h)=X0都存在时刻lJk,lh和对应输入u(k),使输入作用下系统状态在时刻lJk达到原点,即有X(l)=0,则称系统在时刻h完全能控; 如果对初始时刻h和任意非零状态Xl,都存在时刻lJk,lh和对应输入u(k),使输入作用下由初始状态X(h)=0出发的系统运动在时刻lJk达到Xl,则称系统在时刻h完全能达。结论1 离散时间线性时变系统在时刻h完全能达的充分必要条件为,存在时刻lJk,lh,使格兰姆矩阵 为非奇异 1/8,17/451/8,17/45109结论2 若系统矩阵G(k)对所有 kh,l-1 非奇异,则离散时间线性时变
23、系统在时刻hJk完全能控的充分必要条件为,存在时刻lJk,lh,使格兰姆矩阵 为非奇异 若系统矩阵G(k)对一个或一些kh,l-1奇异。格兰姆矩非奇异为系统在时刻h完全能控的一个充分条件。 若系统矩阵G(k) 对所有kh,l-1非奇异,则系统能控性和能达性等价。 若离散时间线性时变系统为连续时间线性时变系统的时间离散化,则系统的能控性和能达性等价。2/8,18/452/8,18/45110时不变系统的能控性和能达性判据 结论3 离散时间线性时不变系统 系统完全能达的充分必要条件为,存在时刻l 0,使格兰姆矩阵 为非奇异。若系统矩阵G非奇异,则系统完全能控的充分必要条件为存在时刻l 0,使格兰姆
24、矩阵为非奇异。若系统矩阵G奇异,则上述格兰姆矩阵非奇异为系统完全能控的充分条件。 3/8,19/453/8,19/45111结论4 n维离散时间线性时不变系统 系统完全能达的充分必要条件为矩阵 满秩 若系统矩阵G非奇异,则系统完全能控的充分必要条件为 rankQkc=n。若系统矩阵G奇异,rankQkc=n 为系统完全能控的一个充分条件。结论5 对于单输入离散时间线性时不变系统,当系统完全能控时,可构造如下一组输入控制 则系统必可在n步内由任意非零初态X(0),转移到状态空间原点,通常称这组控制为最小拍控制。 若系统矩阵G非奇异,则离散时间线性时不变系统能控性和能达性等价。 若离散时间线性时不
25、变系统为连续时间线性时不变系统的时间离散化,则系统的能控性和能达性等价。4/8,20/454/8,20/45112例例 设单输入线性离散系统的状态方程为 试判断系统的能控性,若初始状态x(0)=2,1,0T,确定使x(3)=0的控制序列u(0),u(1),u(2);研究x(2)=0的可能性。 解解 系统是能控的 5/8,21/455/8,21/45113令若令无解。即不存在控制序列u(0),u(1)能够使系统从初始状态x(0)=2,1,0T转移到x(2)=0。6/8,22/456/8,22/45114时变系统的能观测性判据结论6 离散时间线性时变系统在时刻hJk完全能观测的充分必要条件为,存在
26、一个离散时刻lJk,l h,使格兰姆矩阵 为非奇异 时不变系统的能观测性判据 结论7 离散时间线性时不变系统完全能观测的充分必要条件为,存在一个离散时刻l0,使格兰姆矩阵 为非奇异 7/8,23/457/8,23/45115结论8 n 维离散时间线性时不变系统完全能观测的充分必要条件为 满秩 结论9 若单输出离散时间线性时不变系统完全能观测,则利用n步输出值就可构造出相应的初始状态 8/8,24/458/8,24/451164.5 对偶性定义:对连续时间线性时变系统 其对偶系统定义为如下形式的一个连续时间线性时变系统其中,状态X为n维行向量,协状态为n维行向量 输入u为p维列向量,输入为q 维
27、行向量 输出Y为q维列向量,输出为p 维行向量 结论10 :原构系统的状态转移矩阵 与对偶系统的状态转移矩阵 之间满足如下关系 1/2,25/451/2,25/45117结论11 设为原构线性系统, d为对偶线性系统,则有 完全能控 d 完全能观测 完全能观测 d 完全能控 2/2,26/452/2,26/451184.6离散化线性系统保持能控性和能观测性的条件 设连续时间线性时不变系统 对应的时间离散化系统 其中G=eAT H=A的特征值 结论12 如果连续系统(A、B、C)不能控(不能观测),则对任意采样周期T离散化后的系统(G、H、C)也是不能控(不能观测)的。 证明 用反证法 设连续系
28、统不能控,而对于某采样T离散化后的系统却是能控的。则rankH、GH、G2H、Gn-1H=n 1/3,27/451/3,27/45119容易验证 为可交换阵,故 由于eAiT可用I、A、A2、An-1线性表示,故 连续系统是能控的,矛盾。本定理也可叙述为:如果离散化后的系统是能控(能观测)的,则离散化前的连续系统一定是能控(能观测)的。2/3,28/452/3,28/45120结论13 :设连续系统(A、B、C)能控(能观测),则离散化后的系统也能控(能观测)的必要条件是: 不是A的特征值。其中k为非零整数 结论14 对时间离散化,使采样周期T的值 则时间离散化系统能控的充分必要条件是 eAT
29、B为行线性无关 结论15 连续时间线性时不变系统,其时间离散化系统保持完全能控/完全能观测的一个充分条件为,采样周期T满足如下条件:对A的任意两个特征值1、2,不存在非零整数k,使成立对于单输入单输出系统,本定理是充分必要的。3/3,29/453/3,29/451214.7能控性、能观测性与传递函数的关系 结论16 如果A的特征值互不相同,则系统(A、B、C)为能控且能观测的充分必要条件是:传递矩阵G(s)的分母|sI-A|与分子之间不发生因子相消 结论17 单输入、单输出系统(A、b、c)是能控且能观测的充分必要条件是:传递函数G(s)的分母|sI-A|与分子之间不发生因子相消。 结论18
30、单输入、单输出系统(A、b、c),如果A的特征值互不相同,若传递函数存在零、极点对消,则系统或是状态不能控或是状态不能观测的;若传递函数不存在零、极点对消,则系统是状态完全能控且完全能观测的。 证明:单输入、单输出系统动态方程为 如果A的特征值互不相同,则一定可利用非奇异线性变换,使A成为对角阵。即: 1/4,30/451/4,30/45122状态方程可写为:在初始条件为零的情况下,拉氏变换得 对输出方程拉氏变换 此式即为传递函数的部分分式 2/4,31/452/4,31/45123若传递函数存在零、极点对消,传递函数的部分分式中应缺少相应项。如传递函数中相消的零、极点为s-k,则说明fkk=
31、0,k=0,fk 0系统是不能控的;fk=0,k0,系统是不能观测的;k=0,fk=0,系统是既不能控也不能观测的。若传递函数不存在零、极点对消,传递函数的部分分式中,应有fkk0(k=1、2、n)系统是既能控又能观测的。3/4,32/453/4,32/45124例 设单输入、单输出系统的传递函数 由于存在零、极点对消,系统不可能是既能控又能观测的。结论19 如果多输入、多输出系统的状态向量与输入向量之间的传递矩阵 的各行在复数域上线性无关,则系统是能控的。(充分必要条件) 结论20 如果多输入、多输出系统的输出向量与初始状态向量X(0)之间的传递矩阵 的各列在复数域上线性无关,则系统是能观测
32、的。(充分必要条件) 4/4,33/454/4,33/4512548能控规范形和能观测规范形:SISO情形 结论21:连续时间线性时不变系统的能控性和能观测性在线性非奇异变换下保持不变。能控性指数,能观测性指数也保持不变。 定义 一个单输入系统,如果其A、b阵具有如下形式: 则系统一定能控。这种形式的A、b阵称为能控标准形 1/5,34/451/5,34/45126结论22:对完全能控n维单输入单输出连续时间线性时不变系统 则通过变换矩阵 2/5,35/452/5,35/45127可将系统变换成能控规范形,即导出 3/5,36/453/5,36/45128定义 一个单输出系统,如果其A、c阵具
33、有如下形式: 则系统一定能观测,此时的A、c阵称为能观测标准形 结论23:对完全能观测的n 维单输入单输出连续时间线性时不变系统,其能观测规范形可基于线性非奇异变换 导出 4/5,37/454/5,37/45129其中5/5,38/455/5,38/4513049 能控规范形和能观测规范形MIMO情形 旺纳姆能控规范形,旺纳姆能观测规范形龙伯格能控规范形,龙伯格能观测规范形 1/1,39/451/1,39/45131410连续时间线性时不变系统的结构分解 系统按能控性分解 设不能控系统的动态方程为 其能控性矩阵的秩为 rn,选出其中r个线性无关列,再加任意n-r个列,构成非奇异变换T-1 其中
34、 1/6,40/451/6,40/45132经非奇异变换后,系统的动态方程写为 于是可得能控子系统动态方程为: 不能控子系统动态方程为 2/6,41/452/6,41/45133例 已知 试按能控性进行规范分解 解 系统不完全能控,取 能控子系统动态方程为 不能控子系统动态方程为 3/6,42/453/6,42/45134系统按能观测性分解 设不能观测系统的动态方程为 其能观测性矩阵的秩为ln,选出其中l个线性无关行,再加任意n-l个行,构成非奇异变换T 能观测子系统动态方程为 不能观测子系统动态方程为 4/6,43/454/6,43/45135系统按能控性和能观测性的标准分解 设系统(A、B
35、、C)不能控、不能观测,可先对系统按能控性分解,即令 再分别对能控子系统、不能控子系统按能观测性分解 最后得到 5/6,44/455/6,44/45136经T-1变换后,系统的动态方程为 能控、能观测子系统动态方程为:能控、不能观测子系统动态方程为 不能控、能观测子系统动态方程为 不能控、不能观测子系统动态方程为 6/6,45/456/6,45/45137第第5 5章章 系统运动的稳定性系统运动的稳定性 51 外部稳定性和内部稳定性定义:称一个系统的外部稳定(BIBO)是指对任何一个有界输入u(t),即:u(t)1 的任意输入u(t),对应的输出y(t)均为有界,即 结论1:对零初始条件p维输
36、入和q维输出连续时间线性时变系统,tt0,+)则t0时刻系统BIBO稳定的充分必要条件为,存在一个有限正常数,使对一切tt0,+)脉冲响应矩阵H(t,)所有元均满足关系式 证明考虑SISO情形充分性1/4,1/181/4,1/18138必要性采用反证法,即系统BIBO稳定,却存在某个t1使可以取有矛盾结论2:对零初始条件p维输入和q维输出连续时间线性时不变系统,令t0=0,则系统BIBO稳定的充分必要条件为:存在一个有限正常数,使脉冲响应矩阵H(t)所有元均满足关系式 2/4,2/182/4,2/18139结论3:对零初始条件p维输入和q维输出连续时间线性时不变系统,令初始时刻t0=0,则系统
37、BIBO稳定的充分必要条件为:真或严真传递函数矩阵G(s)的所有极点均具有负实部。定义:称连续时间线性时不变系统在t0为内部稳定,是指由时刻t0任意非零初始状态引起的零输入响应Xou(t)对tt0,+)有界,并满足渐近属性,即: 结论4:设n维连续时间线性时变自治系统 系统在t0时刻内部稳定的充分必要条件为:状态转移矩阵(t,t0)对所有tt0,+为有界,并满足: 结论5:对n维连续时间线性时不变自治系统 内部稳定的充分必要条件为 或矩阵A所有特征值均具有负实部,即:Rei(A)0,都对应存在另一位赖于和t0的实数(,t0)0,使得满足不等式X0-Xe(,t0)的任一初始状态x0出发的受扰运动
38、(t;x0,t0)都满足不等式: (t;x0,t0)-Xe 稳定的几何解释李亚普诺夫意义下一致稳定时不变系统的稳定属性李亚普诺夫意义下稳定的实质 1/2,5/181/2,5/18142渐近稳定称自治系统 的孤立平衡状态Xe=0在时刻t0为渐近稳定,如果) Xe=0在时刻t0为李亚普诺夫意义下稳定,)对实数(,t0)0和任给实数0使得满足不等式X0-Xe(,t0)的任一初始状态x0出发的受扰运动(t;x0,t0)满足不等式(t;x0,t0)-Xe, 不稳定 称自治系统 的孤立平衡状态Xe=0在时刻t0为不稳定,如果不管取实数0为多么大,都不存在对应一个实数(,t0)0,使得满足不等式X0-Xe(
39、,t0)的任一初始状态x0出发的受扰运动(t;x0,t0)满足不等式(t;x0,t)-Xe, 不管初始偏差有多大,系统总是稳定的,则称系统是大范围稳定的。不管初始偏差有多大,系统总是渐近稳定的,则称系统是大范围渐近稳定的。大范围渐近稳定的系统只能有一个平衡状态。为了满足稳定条件,初始偏差有一定限制,则称系统是小范围稳定的。对于线性系统,若在小范围稳定,则必大范围稳定;若在小范围渐近稳定,则必大范围渐近稳定 2/2,6/182/2,6/1814353李亚普诺夫第二方法的主要定理 结论7:对连续时间非线性时变自治系统 X=0为系统平衡状态,若可构造对x和t具有连续一阶偏倒数的标量函数V(x,t),
40、V(0,t)=0,且对状态空间中所有非零状态X满足如下条件: )V(x,t)正定且有界,即存在两个连续的非减标量函数(x)和(x),(0)0,(0)0,使对所有tt0,)有:(x)V(x,t)(x)0) V(x,t)对时间t的导数负定且有界。)当x,有V(x,t) 则系统的原点平衡状态x=0为大范围一致渐近稳定。结论8:对连续时间非线性时不变自治系统 X=0为系统平衡状态,若可构造对x具有连续一阶偏倒数的标量函数V(x),V(0)=0,且对状态空间中所有非零状态X满足如下条件: )V(x)为正定) 为负定)当x,有V(x) 则系统原点的平衡状态x=0为大范围一致渐近稳定。 1/4,7/181/
41、4,7/18144例设系统状态方程为坐标原点是系统的一个平衡状态,试确定该系统的稳定性 解 取一正定的标量函数 为一负定的标量函数,且系统是大范围渐近稳定的。 2/4,8/182/4,8/18145结论9 小范围渐近稳定性定理 对连续时间非线性时变自治系统,若可构造对x和t具有连续一阶偏导数的一个标量函数V(x,t),V(0,t)=0,以及围绕状态空间原点的一个吸引区,使对所有非零状态x和所有tt0,)满足如下条件: V(x,t)为正定且有界; 为负定且有界; 则系统原点平衡状态x=0在域内为一致渐近稳定。结论10小范围渐近稳定性定理 对连续时间非线性时不变自治系统,若可构造对x具有连续一阶偏
42、导数的一个标量函数V(x),V(0)=0,以及围绕状态空间原点的一个吸引区,使对所有非零状态x满足如下条件:V(x)为正定; 为负定 则系统原点平衡状态x=0在域内为渐近稳定 3/4,9/183/4,9/18146结论11 小范围渐近稳定性定理 对连续时间非线性时不变自治系统,若可构造对x具有连续一阶偏导数的一个标量函数V(x),V(0)=0,以及围绕状态空间原点的一个吸引区,使对所有非零状态x满足如下条件: V(x)为正定; 为负半定 对任意非零x0 则原点平衡状态x=0在域内为渐近稳定 结论 12不稳定性定理 对连续时间非线性时变自治系统,若可构造对x和t具有连续一阶偏导数的一个标量函数V
43、(x,t),V(0,t)=0,以及围绕状态空间原点的一个吸引区域,使对所有非零状态x和所有tt0,)满足如下条件:()V(x,t)为正定且有界;() 为正定且有界;则系统原点平衡状态x=0为不稳定。 4/4,10/184/4,10/1814754 构造李亚普诺夫函数的规则化方法 变量梯度法设连续时间非线性时不变系统 Xe=0为系统孤立平衡状态, (1)设V(x)的梯度为 (2)设梯度V(x)对应于有势场,则旋度rotV(x)=0,即 (3)由 (4)由(2),(3)定出V(x) (5)1/3,11/181/3,11/18148(6)判断V(x)计算结果的正定性克拉索夫斯基方法设连续时间非线性时
44、不变系统 Xe=0为系统孤立平衡状态, 系统雅可比矩阵 克拉索夫斯基指出:如果存在一个对称正定矩阵B,使对称阵S(x)=BJ(x)+ BJ(x)T是负定的,那么平衡状态x=0是渐近稳定的,系统的李雅普诺夫函数为: V(x)=f(x)TBf(x)如果 则平衡状态x=0是大范围渐近稳定的。结论13:对连续时间线性时不变系统,矩阵A为非奇异,若A+AT为负定,则原点平衡状态x=0为大范围渐近稳定。 2/3,12/182/3,12/18149例 确定平衡状态x=0的稳定性 解 取B=I 为对称负定阵,所以平衡状态x=0是渐近稳定的。平衡状态x=0是大范围渐近稳定的 3/3,13/183/3,13/18
45、15055 连续时间线性系统的状态运动稳定性判据 线性时不变系统的稳定性判据 结论14特征值判据 对连续时间线性时不变系统,原点平衡状态即x=0是李亚普诺夫意义下稳定的充分必要条件为,矩阵A的特征值均具有非正实部即实部为零或负,且零实部特征值只能为A的最小多项式的单根。结论15特征值判据 对连续时间线性时不变系统,原点平衡状态x=0是渐近稳定的充分必要条件为,矩阵A的特征值均具有负实部 结论16李亚普诺夫判据 对n维连续时间线性时不变系统,原点平衡状态xe=0是渐近稳定的充分必要条件为,对任给一个nn正定对称矩阵Q ,李亚普诺夫方程ATP+PA=-Q有唯一nn正定对称解阵P。结论17李亚普诺夫
46、判据推广形式 对n维连续时间线性时不变系统和任给实数0,令矩阵A特征值为i (A), i=1,2,n,则系统所有特征值均位于s平面的直线j左半开平面上,即成立Rei (A)0,使成立: (t,t0)(t0)0使上式成立,系统原点平衡状态xe=0为李亚普诺夫意义下一致稳定。 结论19基于状态转移矩阵的判据 对连续时间线性时变系统,表(t,t0)为系统状态转移矩阵,则系统唯一平衡状态xe=0在时刻t0是渐近稳定的充分必要条件为,存在依赖于t0的一个实数(t0)0,使同时成立:(t,t0)(t0)0和20使成立:(t,t0)1e-2(t-t0)系统原点平衡状态xe=0为一致渐近稳定。2/2,15/1
47、82/2,15/1815256 连续时间线性时不变系统稳定自由运动的衰减性能估计 对渐近稳定的连续时间线性时不变自治系统 衰减系数定义为 最小衰减系数 设 则 1/1,16/181/1,16/1815357 离散时间系统状态运动的稳定性及其判据 结论20大范围渐近稳定判据 对离散时间非线性时不变自治系统,若存在一个相对于离散状态x(k)的标量函数V(x(k),使对任意x(k)n满足:()V(x(k)为正定;()表V(x(k)V(x(k1)V(x(k),V(x(k)为负定;()当x(k),有V(x(k) 。则原点平衡状态即x=0为大范围渐近稳定。结论21大范围渐近稳定判据 对离散时间非线性时不变
48、系统,若存在一个相对于离散状态x(k)的标量函数V(x(k),使对任意x(k)n满足:()V(x(k)为正定;()表V(x(k)V(x(k1)V(x(k),V(x(k)为负半定;()对由任意非零初始状态x(0) n确定的所有自由运动x(k)的轨线,V(x(k)不恒为零;()当x(k),有V(x(k) 。则原点平衡状态即x=0为大范围渐近稳定。1/2,17/181/2,17/18154结论22特征值判据 对离散时间线性时不变自治系统,原点平衡状态即xe=0是李亚普诺夫意义下稳定的充分必要条件为,G的全部特征值i (G)(i=1,2,n)的幅值均等于或小于1,且幅值等于1的特征值只能为G的最小多项
49、式的单根。结论23特征值判据 对离散时间线性时不变自治系统,原点平衡状态是渐近稳定的充分必要条件为,G的全部特征值的幅值均小于1。结论24李亚普诺夫判据 对n维离散时间线性时不变自治系统,原点平衡状态渐近稳定的充分必要条件为,对任一给定nn正定对称矩阵Q,离散型李亚普诺夫方程GTPG-P=-Q有唯一nn正定对称解阵P。 结论25扩展李亚普诺夫判据 对n维离散时间线性时不变自治系统,原点平衡状态以实数0为幂指数稳定,即G的特征值满足:i (G),01, i=1,2,n的充分必要条件为:对任一给定nn正定对称矩阵Q,扩展离散型李亚普诺夫方程 (1/)2GTPG-P=-Q 有唯一nn正定对称解阵P。
50、2/2,18/182/2,18/18155第六章线性反馈系统的时间域综合 61 引言综合问题的提法系统的综合问题由受控系统,性能指标和控制输入三个要素组成 所谓系统综合,就是对给定受控系统,确定反馈形式的控制,使所导出闭环系统的运动行为达到或优于指定的期望性能指标 性能指标的类型 1/1,1/401/1,1/4015662 状态反馈和输出反馈 状态反馈设连续时间线性时不变系统 状态反馈下受控系统的输入为:u=-Kx+,KRpn,反馈系统xf的状态空间描述为: BCAK结论1:对连续时间线性时不变系统,状态反馈保持能控性,不保持能观测性。 1/3,2/401/3,2/40157输出反馈 设连续时
51、间线性时不变系统 输出反馈下受控系统输入u=-Fy+ ,FRpq BCAF输出反馈系统的状态空间描述为: 结论2:对连续时间线性时不变系统,输出反馈保持能控性和能观测性。2/3,3/402/3,3/40158状态反馈和输出反馈的比较 反馈原理反馈原理:状态反馈为系统结构信息的完全反馈,输出反馈则是系统结构信息的不完全反馈。反馈功能反馈功能:状态反馈在功能上远优于输出反馈。改善输出反馈的途径改善输出反馈的途径:扩展输出反馈(动态输出反馈)反馈实现上反馈实现上,输出反馈要优越于状态反馈。解决状态反馈物理实现的途径解决状态反馈物理实现的途径:引入状态观测器扩展状态反馈和扩展输出反馈的等价性等价性。
52、3/3,4/403/3,4/4015963 状态反馈极点配置:单输入情形 极点配置定理 对单输入n 维连续时间线性时不变受控系统 系统全部n个极点可任意配置的充分必要条件为(A,b)完全能控。 极点配置算法step1. 判别(A,b)能控性 step2.计算矩阵A特征多项式det(SI-A)=(s)=Sn+n-1Sn-1+1S+0step3.计算由期望闭环特征值 决定的特征多项式 step4 1/4,5/401/4,5/40160step5 step6. Q=P-1step7 Step8 .停止计算 2/4,6/402/4,6/40161例例1连续时间线性时不变状态方程为期望闭环极点为计算状态
53、反馈阵K解解容易判断 系统能控计算由期望闭环极点组决定的特征多项式 3/4,7/403/4,7/40162计算 4/4,8/404/4,8/4016364 状态反馈极点配置:多输入情况 极点配置定理:对多输入n维连续时间线性时不变系统 系统可通过状态反馈任意配置全部n个极点的充分必要条件为A,B完全能控。 极点配置算法:给定n维多输入连续时间线性时不变受控系统A,B和一组任意期望闭环特征值 要求确定一个pn状态反馈矩阵K,使 step1.判断A的循环性,若非循环,选取一个pn实常阵K1,使 为循环;若循环,表 Step2:选取一个p1实常量,有b=B使 为完全能控 Step3.对等价单输入系统
54、 利用单输入情形极点配置算法,计算状态反馈向量k。 Step4.对A为循环,Kk 对A为非循环,KkK1 1/2,9/401/2,9/40164状态反馈对系统传递函数矩阵零点的影响 结论:对完全能控n维单输入单输出连续时间线性时不变系统,引入状态反馈任意配置传递函数全部n个极点的同时,一般不影响其零点。 结论:对完全能控n维多输入多输出连续时间线性时不变系统,状态反馈在配置传递函数矩阵全部n个极点同时,一般不影响其零点。定义:设完全能控多输入多输出连续时间线性时不变系统 其传递函数矩阵G(s)=C(SI-A) -1 B, G(s)的极点为其特征方程式的根。 零点定义使得 的所有s值 2/2,1
55、0/402/2,10/401656.5 输出反馈极点配置输出反馈极点配置结论:对完全能控连续时间线性时不变受控系统 采用输出反馈 一般不能任意配置系统全部极点。 结论:对完全能控n维单输入单输出连续时间线性时不变受控系统 采用输出反馈 只能使用闭环系统极点配置到根轨迹上,而不能任意配置到根轨迹以外位置上。 1/1,11/401/1,11/401666.6状态反馈镇定状态反馈镇定所谓状态镇定问题就是:对给定时间线性时不变受控系统,找到一个状态反馈型控制律 使所导出的状态反馈型闭环系统 为渐近稳定,即系统闭环特征值均具有负实部。 结论:连续时间线性时不变系统可由状态反馈镇定,当且仅当系统不能控部分
56、为为渐近稳定 结论:连续时间线性时不变系统可由状态反馈镇定的一个充分条件是系统完全能控 状态反馈镇定算法:Step1 判断(A.B)能控性,若完全能控,去Step4。Step2 对(A.B)按能控性分解Step3 对能控部分进行极点配置Step4 计算镇定状态反馈矩阵Step5 计算停止。1/1,12/401/1,12/401676.7状态反馈动态解耦状态反馈动态解耦问题的提法:设多输入多输出连续时间线性时不变系统采用包含输入变换的状态反馈系统BCAKL1/9,13/401/9,13/40168则系统状态空间描述为: 所谓动态解耦控制,就是寻找输入变换 和状态反馈矩阵 使得所导出的闭环传递函数
57、矩阵为非奇异对角有理分式矩阵 系统的结构特征量 输出矩阵 传递函数矩阵 2/9,14/402/9,14/40169对连续时间线性时不变受控系统,结构特性指数定义为: 对连续时间线性时不变受控系统,结构特性向量定义为: 可解耦条件结论:对方连续时间线性时不变受控系统,使包含输入变换状态反馈系统可实现动态解耦的充分必要条件是:基于结构特征向量组成的pp矩阵E非奇异 3/9,15/403/9,15/40170 则可导出包含输入变换状态反馈系统 这种解耦称为积分型解耦系统。4/9,16/404/9,16/40171解耦控制综合算法 给定n维方连续时间线性时不变受控系统 要求综合一个输入变换和状态反馈矩
58、阵对L,K,使系统实现动态解耦,并使解耦后每个单输入单输出系统实现期望极点配置。Step1:计算受控系统(A,B,C)的结构特征量 Step 2:组成并判断矩阵E的非奇异性 若E为非奇异,即能解耦 若E为奇异,则不能解耦。Step3: 5/9,17/405/9,17/40172Step 4:取 导出积分型解耦系统 Step5:判断 的能观测性,若不完全能观测,计算 Step6:引入线性非奇异变换 化积分型解耦系统为解耦规范型。 对完全能观测 6/9,18/406/9,18/401737/9,29/407/9,29/40174Step7 求 Step8:对解耦规范型 选取 状态反馈矩阵 的结构
59、对完全能观测 对不完全能观测 8/9,20/408/9,20/40175Step9:根据指定期望极点组按单输入情形极点配置法,定出状态反馈矩阵 Step10:最后得 例6.4 p2989/9,21/409/9,21/401766.8状态反馈静态解耦状态反馈静态解耦问题的提法: 设多输入多输出连续时间线性时不变受控系统 所谓静态解耦,就是综合一个输入变换和状态反馈矩阵对 使导出的包含输入变换状态反馈系统及其传递函数矩阵 满足:i)闭环控制系统渐近稳定, 即ii)闭环传递函数矩阵当S=0时为非奇异对角常阵,即有 1/2,22/401/2,22/40177可解耦条件 存在输入变换和状态反馈矩阵对L,
60、K,其中 可使方n维连续时间线性时不变受控系统实现静态解耦,当且仅当 2/2,23/402/2,23/40178612 全维状态观测器 一:开环状态观测器为了实现状态反馈,有时需要对状态进行估计,开环估计方法如下:二:全维观测器全维观测器是指重构状态向量的维数与原系统相同 事实上,已知的信息为u(t)和y(t),只有当系统完全能观测时,才能从u(t)和y(t)及其导数的线性组合中获得状态向量x(t)的估计值此时存在状态观测器。状态观测器利用观测器实现状态反馈的系统为: 1/5,24/401/5,24/40179在观测器的设计中,为使 尽快地接近x(t),可利用y(t)和 之间的差作为误差反馈信
61、息,观测器结构如下:BACH写出观测器动态方程为 原系统的状态方程: 定义状态向量的真实值与估计值之间的偏差为误差状态向量,即:2/5,25/402/5,25/40180定理:若系统(A、B、C)是能观测的,其状态可用n维状态观测器进行估计 矩阵H可以按给定极点的位置来选择,所定极点的位置,将决定误差向量趋于零的速率。3/5,26/403/5,26/40181例 设系统动态方程为 试设计一个状态观测器,其中矩阵A-hc的特征值(观测器极点)为-10,-10。解 希望的特征多项式 观测器方程 4/5,27/404/5,27/40182原系统及其状态观测器结构图如下 5/5,28/405/5,28
62、/401836.13降维观测器一个n维的能观测系统由于y可以直接提供一部分状态,故只需要估计其余的状态即可。1:建立n-m维子系统动态方程 设ARnn,BRnr,CRmn,系统(A、B、C)能观测,令:为一个nn矩阵,D的选择应使Q可逆,考虑到 系统的动态方程为 1/6,29/401/6,29/40184可直接有y提供,只须估计 2:降维观测器设计 方程改写为 故降维观测器方程为 2/6,30/402/6,30/40185令 这是一个n-m维观测器,整个状态向量的估计值为: 而系统原状态向量x的估计值为 3:H阵的选择 通过H阵的选择,使 的极点任意配置,极点的位置决定误差向量 衰减到零的速率
63、,而 直接有y提供,不存在估值误差。 3/6,31/403/6,31/40186定理:有m个输出的任一m维能观测系统(A、B、C),可通过状态变换而写成如下形式:其状态可用n-m维龙伯格观测器进行估计 (n-m)m矩阵H可以选得使 极点的位置决定误差向量衰减到零的速率。观测器结构图如下: 的极点任意配置4/6,32/404/6,32/40187例 已知系统: 试构造一降维观测器解 系统完全能观测 令 设降维观测器的特征值为-10,H=h希望的特征多项式为+10,故H=10,降维观测器为: 5/6,33/405/6,33/40188原系统状态向量估计值为 原系统及其降维观测器如下 原系统降维观测
64、器6/6,34/406/6,34/401896.14带状态观测器的状态反馈系统 带状态观测器的状态反馈系统: 现提出两个问题:1,用状态估计进行状态反馈和用x进行状态反馈对系统特性的影响是否一致,或者说系统的闭环传递矩阵是否一致?2,进行状态反馈设计时的K阵和观测器设计时的H阵能否分开设计?显然利用x进行状态反馈时,控制系统的传递矩阵为 状态估计进行状态反馈(A、B、C)状态观测器K1/6,35/401/6,35/40190为了计算传递矩阵,作坐标变换 变换前 变换后 考虑到当R、T可逆时,有2/6,36/402/6,36/40191传递矩阵为 系统的传递矩阵与用x进行状态反馈时的传递矩阵相同
65、。 另外,系统的特征多项式 它由(A-BK)、(A-HC)的特征多项式的乘积组成,可见只要系统(A、B、C)能控、能观测,则可按极点配置的需要选择K,按观测器动态特性的需要H,两者可分开进行设计,这个原理称为分离定理 3/6,37/403/6,37/40192例 设系统的传递函数为 希望利用状态反馈使闭环极点为-4j6,并求实现这个反馈的二维及一维观测器。 解 1:建立能观测标准形实现 系统也是能控的 2:求状态反馈阵K,设K=k1,k2,系统特征方程式:希望的特征方程式K=2,40 3:求二维观测器,设其极点为s1=s2=-10,H=h1,h2T 希望的特征方程式 H=100,14T4/6,
66、38/404/6,38/40193观测器方程 系统结构图 5/6,39/405/6,39/401944:求一维观测器,设其极点为s1=-10,H=h 希望的特征方程式s+10=0H=10 观测器方程系统结构图 yuv-6/6,40/406/6,40/40195第二部分第二部分 线性系统的复频率域理论线性系统的复频率域理论 线性系统的复频率域理论,是以传递矩阵作为系统描述,并在复频率域内分析和综合线性时不变系统的一种理论和方法。 19681 矩阵分式描述矩阵分式描述第第8 8章章 传递函数矩阵的矩阵分式描述传递函数矩阵的矩阵分式描述 矩阵分析描述实质上就是把有理分式矩阵形式的传递函数矩阵G(s)
67、表为两个多项式矩阵之“比”。右MFD和左MFD设p维和n维输出连续时间线性时不变系统其传递函数矩阵 为G(s)的一个右矩阵分式描述 为G(s)的一个左矩阵分式描述 其中 为多项式矩阵 1/4,1/161/4,1/16197例如 上式即为G(s)的一个右MFD 把G(s)按各行通分,可以写出G(s)的左MFD MDF的特性 结论:对传递函数矩阵G(s)的一个右MFD,规定对传递函数矩阵G(s)的一个左MFD,规定 对给定一个G(s),其右MFD和左MFD在次数上一般不相等。 2/4,2/162/4,2/16198结论:对传递函数矩阵G(s),其右MFD和左MFD为不唯一,且不同的MFD可能具有不
68、同的次数。结论:对qp传递函数矩阵G(s),设 W(s)为pp非奇异多项式矩阵,令 为其一个右MFD则 也是G(s)的一个右MFD,且 若W(s)为单模矩阵,则 结论:对qp传递函数矩阵G(s),设 WL(s)为任一qq非奇异多项式矩阵。 为其一个左MFD,则 也是G(s)的一个左MFD,且 若WL(s)为单模矩阵,则 3/4,3/163/4,3/16199结论:对qp传递函数矩阵G(s),设 为其 一个右MFD和一个左MFD 则有 *最小阶MFD也不是唯一的 *称最小阶MFD为不可简约MFD 4/4,4/164/4,4/162008.2矩阵分式描述的真性和严真性矩阵分式描述的真性和严真性 设
69、多输入多输出连续时间线性时不变系统,传递函数矩阵G(s)为 1/3,5/161/3,5/16201结论:定义 真性和严真性的判别准则 结论:对右MFD D(s)为pp阵且列既约 则 例 容易判断D(s)为列既约,且 可知 为真 2/3,6/162/3,6/16202例 给定12右MFD D(s)为非列既约,尽管 但 非真 结论:对左MFD 为qq阵且行既约,则 * 若D(s)或DL(s)为非列既约或行既约,则引入一个单模矩阵,化D(s)或DL(s)为列既约或行既约,进行判断。 3/3,7/163/3,7/162038.3从非真矩阵分式描述导出严真矩阵分式描述从非真矩阵分式描述导出严真矩阵分式描
70、述 结论:对非真右MFD,N(s)D-1(s),D(s) 为pp多项式矩阵,N(s)为qp多项式矩阵Q(s)和R(s),使且R(s)D-1(s)为非真N(s)D-1(s)导出的严真右MFD。 确定严真MFD的算法 Step1:计算给定N(s)D-1(s)的有理分式矩阵G(s) Step2:通过多项式除法,得 Step 3 Step4 其中R(s)D-1(s)为非真右MFD N(s)D-1(s)的严真部分,Q(s)为多项式矩阵部分。 1/3,8/161/3,8/16204结论:对非真左MFD,DL-1(s)NL(s),唯一存在两个多项式矩阵 使 一类特殊情形的多项式矩阵除法问题 在连续时间线性时
71、不变系统中,除式矩阵通常为sI-A 结论:对pp矩阵sI-A和多项式矩阵N(s),唯一存在一个常阵Nr(A)和多项式矩阵Qr(s)满足 其中 显然Nr(A)(sI-A) -1为N(s)(sI-A) -1所导出的严真右MFD 2/3,9/162/3,9/16205结论:对qq矩阵sI-A和多项式矩阵NL(s),唯一存在一个常阵NL(A)和多项式矩阵QL(s)满足 其中 显然(sI-A) -1 NL(A) 为 (sI-A) -1 NL(s)所导出的严真左MFD 3/3,10/163/3,10/1620684 不可简约矩阵分式描述不可简约MFD实质上是系统传递函数矩阵的一类最简约MFD,通常也称为最
72、小阶MFD。 定义:右不可简约D(s)和N(s)为右互质 左不可简约DL(s)和NL(s)为左互质 不可简约MFD的基本特性结论:对qp传递函数矩阵,其右不可简约MFD和左不可简约MFD均为不惟一 结论:设 为qp传递函数矩阵G(s)的任意两个右不可简约MFD,则必存在单模阵U(s)满足: 1/3,11/161/3,11/16207结论:设 为传递函数矩阵G(s)的任意两个左不可简约MFD 则必存在单模阵V(s),满足 结论:传递函数矩阵G(s)的右不可简约MFD满足广义惟一性。 传递函数矩阵G(s)的左不可简约MFD满足广义惟一性。结论:对qp传递函数矩阵G(s)的任一右不可简约MFD N(
73、s)D -1 (s)和任一右可简约MFD ,必存在非奇异多项式矩阵T(s),满足:结论:对qp传递函数矩阵G(s)的任一左不可简约MFD DL -1 (s) NL(s) 和任一左可简约MFD ,必存在非奇异多项式矩阵TL(s),满足:2/3,12/162/3,12/16208结论:对qp传递函数矩阵G(s)的所有右不可简约MFD 必有:1,Ni(s)具有相同史密斯形 2,Di(s)具有相同不变多项式结论:对qp传递函数矩阵G(s)的所有左不可简约MFD 必有:1,NLi(s)具有相同史密斯形 2,DLi(s)具有相同不变多项式结论:对qp传递函数矩阵的任一左不可简约MFD,和任一右不可简约MF
74、D 必有 3/3,13/163/3,13/1620985确定不可简约矩阵分式描述的算法结论:对qp传递函数矩阵G(s)设 为任一右可简约MFD pp多项式矩阵R(s)为 的一个最大右公因子且为非奇异,取 为G(s)的一个右不可简约MFD。 结论:对qp传递函数矩阵G(s)设 为任一左可简约MFD RL(s)为 的一个最大左公因子且为非奇异,取 为G(s)的一个左不可简约MFD。 1/1,14/161/1,14/1621086规范矩阵分式描述规范矩阵分式描述传递函数矩阵的可简约MFD和不可简约MFD具有不惟一性。其惟一化的途径是对MFD分母矩阵限定为规范形而得到规范MFD。埃尔米特形MFD称qp
75、的NH(s)DH-1(s)为传递函数矩阵G(s)的列埃尔米特形MFD,是指分母矩阵具有列埃尔米特形 其中:1) 为首1多项式 2)若 则 3)若 为含S多项式,则 1/2,15/161/2,15/16211称qp的 为传递函数矩阵G(s)的行埃尔米特形MFD,是指分母矩阵 具有行埃尔米特形 其中:1) 为首1多项式 2)若 则 3)若 为含S多项式,则 结论:对qp传递函数矩阵G(s),其所有不可简约右MFD具有相同列埃尔米特形MFD 其所有不可简约左MFD具有相同行埃尔米特形MFD 2/2,16/162/2,16/16212第9章 传递函数矩阵的结构特性91史密斯-麦克米伦形称秩为r的有理分
76、式矩阵为史密斯-麦克米伦形,当且仅当具有形式 其中,1) 为互质,i=1,2,r 2)满足整除性 1/4,1/121/4,1/12213结论:对qp有理分式矩阵G(s),设则必存在qq和pp单模矩阵U(s)和V(s)使变换后传递函数矩阵U(s)G(s)V(s)为史密斯-麦克米伦形史密斯-麦克米伦形基本特性结论:有理分式矩阵G(s)的史密斯-麦克米伦形M(s)为惟一结论:化有理分式矩阵G(s)为史密斯-麦克米伦形M(s)的单模变换阵对U(s),V(s)不惟一。结论:严格有理分式矩阵G(s)的史密斯-麦克米伦形M(s)不具有保持严真属性,M(s)甚至可能为非真。结论:对qq非奇异有理分式矩阵G(s
77、)其中a为非零常数2/4,2/122/4,2/12214结论:对秩为r的qp传递函数矩阵G(s),其史密斯-麦克米伦形M(s)为令 则M(s)表为右MFD 令 则M(s)表为左MFD 3/4,3/123/4,3/12215结论:对qp传递函数矩阵G(s),其史密斯-麦克米伦形为M(s)。单模变换阵对为U(s),V(s)若取 则 为G(s)的不可简约右MFD 若取 则 为G(s)的不可简约左MFD 4/4,4/124/4,4/1221692传递函数矩阵的有限零点和有限极点定义:对秩为r的qp传递函数矩阵G(s),其史密斯-麦克米伦形M(s),则G(s)有限极点=M(s)中 的根,i=1,2,rG
78、(s)有限零点=M(s)中 的根, i=1,2,r 定义:对qp传递函数矩阵G(s),设 和 为G(s)的任一不可简约右MFD和任一不可简约左MFD,则 G(s)有限极点=det(D(s)=0根或 的根 G(s)有限零点= 的s值或 的s值 定义:对qp传递函数矩阵G(s),设其状态空间描述为(A,B,C),且(A,B)全能控,(A,C)完全能观测,则有: G(s)有限极点= 的根 G(s)有限零点=使 降秩的s值 1/2,5/121/2,5/12217结论:对qp严真传递函数矩阵G(s),其能控和能观测状态空间描述为(A,B,C),z0为任一零点,则对满足关系式 的所有非零初始状态x0和输入
79、 系统输出具有阻塞作用,即其能引起的系统输出y(t)强制恒为零。 2/2,6/122/2,6/1221893传递函数矩阵的结构指数 对qp传递函数矩阵G(s), 有限零点和有限极点的集合。 那么,若对任一 导出对应的rr对角矩阵 则称 为G(s)在 的一组结构指数 1/3,7/121/3,7/12219例定出 的结构指数 史密斯-麦克米伦形为 G(s)极点零点集合 =正整数G(s)在S= 有 个零点 =负整数G(s)在S= 有 个极点 =0 G(s)在S= 0 处无零点和极点2/3,8/122/3,8/12220结论:G(s)在 极点重数= 中负指数之和绝对值 结论:G(s)在 零点重数= 中
80、正指数之和 结论:传递函数矩阵在非极点零点处的结构指数必恒为零。结论:对qp传递函数矩阵G(s)设 则G(s)的史密斯-麦克米伦形M(s)为 上式表明,一旦定出G(s)各个极点零点及其结构指数组,便可构造出G(s)的史密斯-麦克米伦形M(s)。3/3,9/123/3,9/1222194传递函数矩阵在无穷远处的极点和零点 确定s=处极点零点的思路对qp传递函数矩阵G(s), 则直接基于G(s)的史密斯-麦克米伦形M(s)不能定义G(s)在无穷远处的极点和零点,若引入变换 则有G(s)在s=处的极点/零点=H()在=0处的极点/零点。结论:对qp传递函数矩阵G(s),设 再基于变换 由G(-1)导
81、出H() 引入单模变换阵。导出其史密斯-麦克米伦形 G(s)在s=处的极点重数= 中 的 根重数i=1,2,r 则有G(s)在s=处的零点重数= 中 的 根重数i=1,2,r 1/3,10/121/3,10/12222例:设 史密斯-麦克米伦形 基于此,可以定出G(s)在s=处极点重数=2G(s)在s=处零点重数=1 2/3,11/122/3,11/12223无穷远处的结构指数对qp传递函数矩阵G(s) 则G(s)在s=处结构指数 在=0处结构指数 3/3,12/123/3,12/12224第10章 传递函数矩阵的状态空间实现101实现的基本概念和基本属性定义101实现对真或严真连续时间线性时
82、不变系统,称一个状态空间描述 或简写为(A,B,C,E)是其传递函数矩阵G(s)的一个实现,如果两者为外部等价即成立关系式:C (sIA) -1 B+E=G(s) 结论101实现维数传递函数矩阵G(s)的实现(A,B,C,E)的结构复杂程度可由其维数表征。一个实现的维数规定为其系统矩阵A的维数,即有 实现维数=dimA 结论102不惟一性传递函数矩阵G(s)的实现(A,B,C,E)满足强不惟一性。即对传递函数矩阵G(s),不仅其实现结果为不惟一,而且其实现维数也为不惟一。 结论103最小实现最小实现定义为传递函数矩阵G(s)的所有实现(A,B,C,E)中维数最小的一类实现。实质上,最小实现就是
83、外部等价于G(s)的一个结构最简状态空间模型。 1/5,1/391/5,1/39225结论104实现间关系对传递函数矩阵G(s),其不同实现间一般不存在代数等价关系,但其所有最小实现间必有代数等价关系。 结论105实现物理本质物理直观上,传递函数矩阵G(s)的实现就是对具有“黑箱”形式的真实系统在状态空间领域寻找一个外部等价的内部假想结构,内部假想结构对真实系统的可否完全表征性依赖于系统的是否能控和能观测。结论106实现形式传递函数矩阵G(s)的实现形式取决于其真性或严真性属性。当G(s)为严真,其实现对应地具有形式(A,B,C)即E=0;当G(s)为真,其实现对应地具有形式(A,B,C,E)
84、即E0,且有设状态空间描述(A,B,C,E)为传递函数矩阵G(s)的一个实现,dimA=n,则对任一nn非奇异阵T,状态空间描述(TAT-1,TB,CT-1,E)必也为G(s)的一个同维实现。其他实现构造2/5,2/392/5,2/39226能控类实现和能观测类实现是两类基本的典型实现定义102能控类实现 称状态空间描述(A,B,C,E)为传递函数矩阵G(s)的一个能控类实现,当且仅当 C(sIA) -1 B+E=G(s) (A,B)能控且有指定形式 定义103能观测类实现称状态空间描述(A,B,C,E)为传递函数矩阵G(s)的一个能观测类实现,当且仅当 C(sIA) -1 B+E=G(s)
85、(A,C)能观测且有指定形式 最小实现是传递函数矩阵G(s)的一类最为重要的实现。最小实现是G(s)的所有实现中结构为最简的实现,即从外部等价的角度实现中不包含任何多余的部分,因此通常也称最小实现为不可简约实现。 结论108设(A,B,C)为严真传递函数矩阵G(s)的一个实现,则其为最小实现的充分必要条件是(A,B)完全能控,(A,C)完全能观测 最小实现判据3/5,3/393/5,3/39227实现的最小维数结论1010实现最小维数 对严真传递函数矩阵G(s),表其为幂级数表达式:为马尔柯夫(Markov)参数矩阵,并基此组成汉克尔(Hankel)矩阵 则G(s)的状态空间实现的最小维数为
86、nmin=rankH 结论109最小实现广义惟一性严真传递函数矩阵G(s)的最小实现为不惟一但满足广义惟一性。即若(A,B,C)和 为G(s)的任意两个n维最小实现,则必可基此构造出一个nn非奇异常阵T使成立: 4/5,4/394/5,4/39228结论1011实现最小维数 对qp传递函数矩阵G(s),rankG(s) =r,其史密斯麦克米伦形为 其中,U(s)和V(s)为qq和pp单模阵。那么,G(s)的状态空间实现的最小维数为 5/5,5/395/5,5/3922910.2标量传递函数的典型实现不失一般性,考虑真标量传递函数g(s),并通过严真化先将其表为常数e和严真有理分式n(s)/d(
87、s)之和,即有 那么,对g(s)的各类典型实现就归结为对严真传递函数n(s)/d(s)导出相应的实现,而常数e为各类实现中的输入输出直接传递函数。1/5,6/391/5,6/39230几点讨论 真标量传递函数g(s)的能控规范形实现 实现形式惟一性 维数非最小性 (Ac,bc,cc)为最小实现条件 结论1012能控规范形实现 标量传递函数g(s)的严真部分n(s)/d(s)的能控规范形实现具有形式:2/5,7/392/5,7/39231结论1018能观测规范形实现 标量传递函数g(s)的严真部分n(s)/d(s)的能观测规范形实现具有形式:几点讨论 真标量传递函数g(s)的能观测规范形实现 实
88、现形式惟一性 维数非最小性 (Ac,bc,cc)为最小实现条件 结论10.24对偶性 严真标量传递函数n(s)/ d(s)的能控规范形实现(Ac,bc,cc)和能观测规范形实现(A0,b0,c0)满足对偶关系,即有A0= AcT,b0= ccT,c0= bcT 3/5,8/393/5,8/39232结论10.25并联形实现 设传递函数g(s)及其严真部分n(s)/ d(s),极点为1(1重),2(2重),m(m重), 表 则严真传递函数n(s)/ d(s)的并联形实现为 4/5,9/394/5,9/39233几点解释 并联形实现为约当型规范形实现 并联形实现在构成上的难点 对极点中包含共轭复数
89、情形的处理 表n(s)/ d(s)为 则严真传递函数n(s)/ d(s) 的串联形实现为 几点解释 (1)串联形实现的优点(2)串联形实现在构成上的难点(3)对极零点中包含共轭复数情形的处理 5/5,10/395/5,10/3923410.3基于有理分式矩阵描述的典型实现:能控形实现和能观测形实现考虑以有理分式矩阵描述给出的真qp传递函数矩阵G(s)G(s)=(gij(s),i=1,q j=1,q 进而,表G(s)为“严真qp传递函数矩阵”和“qp常阵E”之和,即G(s)=(gij(s)=(eij)+(gijsp(s)=E+Gsp(s) 且有E=G()。再表Gsp(s)诸元即G(s)诸元的最小
90、公分母d(s)为d(s)=sl+l-1sl-1+1s+0基此,严真qp传递函数矩阵Gsp(s)可进而表为其中,Pk(k=0,1,l-1)为qp常阵 1/3,11/391/3,11/39235结论10.35能控形实现对以有理分式矩阵描述给出的严真传递函数矩阵Gsp(s),其能控形实现具有形式 而真传递函数矩阵G(s)的能控形实现为 2/3,12/392/3,12/39236结论10.36能观测形实现 对以有理分式矩阵描述给出的严真传递函数矩阵Gsp(s),其能观测形实现具有形式 而真传递函数矩阵G(s)的能观测形实现为 3/3,13/393/3,13/3923710.4基于矩阵分式描述的典型实现
91、:控制器实现和观测器形实现 右MFD的控制器形实现不失一般性,考虑qp右MFD 和D(s)为qp和pp的多项式矩阵,设为D(s)列既约 首先,对真 导出其严真右MFD。为此,通过矩阵除法,可以得到其中,qp常阵E为“商阵”,qp多项式矩阵N(s)为“余式阵”。基此,将上式右乘D-1(s),可以导出 下面的问题就是,对qp严真右MFD N(s)D-1(s),D(s)列既约构造其控制器形实现。1/22,14/391/22,14/39238(1)控制器形实现的定义定义10.4控制器形实现对qp严真右MFD N(s)D-1(s),D(s)列既约,表列次数ciD(s)=kci,i=1,2,p,称一个状态
92、空间描述 为其控制器形实现,其中如果满足Cc(sIAc) -1 Bc= N(s)D-1(s) (AC,BC)为完全能控且具有特定形式2/22,15/392/22,15/39239Dhc为D(s)的列次系数,且detDhc0 DLc为D(s)的低次系数阵 NLc为N(s)的低次系数阵 结论10.37构造(AC,BC,CC)的结构图对qp严真右MFD N(s)D-1(s),D(s)列既约,表列次数ciD(s)=kci,i=1,2,p, 再引入列次表达式:D(s)=DhcSc(s)+DLcC(s) N(s)=NLcC(s)其中3/22,16/393/22,16/39240那么,基此可导出构造(AC,
93、BC,CC)的结构图 称c(s)Sc-1(s)为核心右MFD -uyu0y0图10.5结论10.38构造(AC,BC,CC)的思路给定qp严真右MFD N(s)D-1(s),D(s)列既约,则在图10.5所示构造(AC,BC,CC)的结构图基础上,对(AC,BC,CC)的构造可分为两步进行:首先,对核心右MFDc(s)Sc-1(s)构造实现(Ac0,Bc0,Cc0),称其为N(s)D-1(s)的核实现。进而,用核实现置换图10.5所示结构图中的核心右MFD,再通过结构图化简导出N(s)D-1(s)的控制器形实现。4/22,17/394/22,17/39241(3)核实现(Ac0,Bc0,Cc0
94、)的构造 先来引入积分链组模型。相对于qp右MFD N(s)D-1(s),D(s)列既约,列次数ciD(s)=kci,i=1,2,p,其积分链组的组成如图所示。图中,为使组成表达整齐起见,已经非实质性地假定列次数满足非降性,即成立kc1kc2kcp。 积分链组的输入uch取为 积分链组的输出ych取为各个积分链的输出构成的向量 5/22,18/395/22,18/392426/22,19/396/22,19/39243结论10.40积分链组的状态空间描述 相对于qp右MFD N(s)D-1(s)的积分链组模型,取状态Xch输出Ych和输入uch为 7/22,20/397/22,20/39244
95、8/22,21/398/22,21/39245(4)控制器形实现的构造结论10.42控制器形实现 对真qp右MFD,其严真右MFD为N(s)D-1(s),D(s)列既约,列次数ciD(s)=kci,i=1,2,p,再引入列次表达式:D(s)=DhcSC(s)+DLcc(s) N(s)=NLcc(s)且知核MFDc(s)Sc-1(s)的实现为(Ac0,Bc0,Cc0),则严真N(s)D-1(s)的控制器形实现(AC,BC,CC)的系数矩阵为Ac= Ac0Bc0Dhc-1DLc,Bc= Bc0Dhc-1,Cc=NLc 而真右MFD的控制器形实现为(AC,BC,CC,E) Bc0Cc0NLcD-1h
96、cD-1hcDLcAc09/22,22/399/22,22/39246例101 定出给定22右MFD N(s)D-1(s)的控制器形实现(Ac,Bc,Cc),其中 容易判断,D(s)为列既约,且N(s)D-1(s)为严真,进而,定出列次数 kc1=c1D(s)=2,kc2=c2D(s)=3 基此,又可定出10/22,23/3910/22,23/39247核实现 可导出控制器形实现 11/22,24/3911/22,24/39248控制器形实现的性质结论10.43控制器形实现 对严真右MFD N(s)D-1(s),D(s)列既约,由核实现(Ac0,Bc0,Cc0)的结构所决定,其控制器形实现(A
97、C,BC,CC)具有形式:12/22,25/3912/22,25/39249(2)控制器形实现和列次表达式在系数阵间的对应关系结论10.44对应关系 对严真右MFD N(s)D-1(s),D(s)列既约,控制器形实现系数矩阵(AC,BC,CC)和D(s)列次表达式系数阵之间具有直观关系 Ac的第i个*行Dhc-1DLc的第i行Bc的第i个*行Dhc-1的第i行其中,i=1,2,p。 例102 定出给定22右MFD N(s)D-1(s)的控制器形实现(Ac,Bc,Cc),其中 容易判断,D(s)为列既约,且N(s)D-1(s)为严真,进而,定出列次数 kc1=c1D(s)=2,kc2=c2D(s
98、)=3 13/22,26/3913/22,26/39250基此,又可定出结论10.45不完全能观测属性 对严真右MFD N(s)D-1(s),D(s)列既约的控制器形实现(AC,BC,CC),(AC,BC)为完全能控,但(AC,CC)一般为不完全能观测。14/22,27/3914/22,27/39251结论10.46系数矩阵间关系 对严真右MFD N(s)D-1(s),D(s)列既约,其控制器形实现(AC,BC,CC)和N(s)D-1(s)在系数矩阵之间具有关系:结论10.47系数矩阵行列式间关系 对严真右MFD N(s)D-1(s),D(s)列既约,其控制器形实现(AC,BC,CC)和N(s
99、)D-1(s)在系数矩阵行列式之间具有关系:det(sI-Ac)=(detDhc) -1 detD(s)dim(Ac)=deg(det D(s) 结论10.48实现和N(s)关系 对严真右MFD N(s)D-1(s),D(s)列既约,其控制器形实现(AC,BC,CC)和MFD分子矩阵N(s)之间具有关系 结论10.49联合能控能观测条件 对严真右MFD N(s)D-1(s),D(s)列既约,其控制器形实现(AC,BC,CC) 联合能控和能观测的一个充分条件为,对所有s,qp矩阵N(s)为列满秩即rankN(s)=p 15/22,28/3915/22,28/39252左MFD的观测形实现 考虑真
100、qp左MFD 为多项式矩阵,为既约行列式。为对真 导出严真左MFD,引入矩阵左除法可以得到 其中,DL-1 (s)NL(s)为严真左MFD。下面的问题就是,对qp严真左MFDDL-1 (s)NL(s),DL(s)行既约,构造观测器形实现 定义10.5观测器形实现对qp严真左MFD DL-1 (s)NL(s),DL(s)行既约,表行次数rjDL(s)=krj,j=1,2,q,则称一个状态空间描述 16/22,29/3916/22,29/39253(2)核实现(A00 B00 C00)对严真DL-1 (s)NL(s),行次数rjDL(s)=krj,j=1,2,q,引入行次数表达式 DL(s)=Sr
101、(s)Dhr+r(s)DLrNL(s)= r(s) NLr其中Dhr为DL(s)的行次系数矩阵,且detDhr0DLr为DL(s)的低次系数阵NLr为N(s)的低次系数阵 17/22,30/3917/22,30/39254结论10.51核实现对qp左MFD DL-1 (s)NL(s),其核心MFD Sr-1 (s)r(s)的实现即DL-1 (s)NL(s)的核实现为严真DL 1 (s)NL(s)的观测器形实现(A0,B0,C0)的系数矩阵关系式为 A0=A00DLr(s)Dhr-1C00 B0=NL , C0=Dhr-1C00 18/22,31/3918/22,31/39255观测器形实现的性
102、质结论10.53观测器形实现 对严真左MFD DL1 (s)NL(s),DL(s)行既约,其观测器形实现(A0,B0,C0)具有形式:19/22,32/3919/22,32/39256结论10.54对应关系 对严真左MFD DL1 (s)NL(s),DL(s)行既约,观测器形实现(A0,B0,C0)系数矩阵和DL(s)列次表达式系数矩阵之间具有直观关系:A0的第j个*列=DLrDhr-1的第j列C0的第j个*列=Dhr-1的第j列j=1,2,q。结论10.55不完全能控属性 对严真左MFD DL1 (s)NL(s),DL(s)行既约,则其观测器形实现(A0,B0,C0)中,(A0,C0)为完全
103、能观测,但(A0,B0)一般为不完全能控。20/22,33/3920/22,33/39257结论10.56系数矩阵间关系 对严真左MFD DL1 (s)NL(s),DL(s)行既约,其观测器形实现(A0,B0,C0) 和DL1 (s)NL(s)在系数矩阵之间具有直观关系: 结论10.57系数矩阵行列式间关系 对严真左MFD DL1 (s)NL(s),DL(s)行既约,其观测器形实现(A0,B0,C0) 和DL1 (s)NL(s)在系数矩阵的行列式之间具有直观关系:det(sI-A0)=(detDhr) -1 detDL(s) dim(A0)=deg detDL(s) 结论10.58实现和NL(
104、s)关系 对严真左MFD DL1 (s)NL(s),DL(s)行既约,其观测器形实现(A0,B0,C0) 和MFD 的分子矩阵NL(s)之间具有关系: 21/22,34/3921/22,34/39258结论10.59联合能控能观测条件 对严真左MFD DL1 (s)NL(s),DL(s)行既约,其观测器形实现(A0,B0,C0) 联合能控和能观测的一个充分条件为,对所有s,qp矩阵NL(s)为行满秩即rank NL(s)=q。 结论10.60对偶性设(A0,B0,C0)为“严真左MFD DL1 (s)NL(s),DL(s)行既约”的观测器形实现,(Ac,Bc,Cc)为“严真右MFD N(s)
105、D1 (s) ,D(s)列既约”的控制器形实现,则(A0,B0,C0)和(Ac,Bc,Cc) 形式为对偶,即A0(=)AcT ,C0(=)BcT 22/22,35/3922/22,35/3925910.5基于矩阵分式描述的典型实现:能控性形实现和能观测性形实现基于矩阵分式描述的实现按“右或左MFD”和“分母矩阵列既约或行既约”共有四种可能的组合。上节已就“右MFD N (s) D1 (s),D (s)列既约”和“左MFD DL1 (s)NL(s),DL(s)行既约”构造“控制器形实现”和“观测器形实现”。本节讨论“右MFD N (s) D1 (s),D (s)列既约”和“左MFD DL1 (s
106、)NL(s),DL(s)列既约”构造对应的“能控性形实现”和“能观测性形实现”。 1/1,36/391/1,36/3926010.6不可简约矩阵分式描述的最小实现最小实现也称为不可简约实现。最小实现是传递函数矩阵的维数最小即结构最简约的一类实现。结论10.78不可简约右MFD最小实现 对qp严真右MFD N (s)D1 (s),设n=deg detD(s),表(Ac,Bc,Cc)为“N (s) D1 (s),D (s)列既约”的 n维控制器形实现,则有(Ac,Bc,Cc)为最小实现 N (s) D1 (s)不可简约表(Aco,Bco,Cco)为“N (s) D1 (s),D (s)行既约”的
107、n维能控性形实现,则有 (Ac0,Bc0,Cc0) 为最小实现 N (s) D1 (s)不可简约 需要指出,尽管上述结论为由右MFD确定最小实现提供了一条易于计算的途径,但这并不意味着由右MFD的最小实现只可能有控制器形或能控形的形式。下面,给出右MFD的最小实现的更具普遍性的结论。 结论10.79不可简约右MFD最小实现 对qp严真右MFD N (s)D1 (s),D (s)列或行既约,表(A,B,C)为其任意形式的n维实现,n=deg detD(s),则有(A,B,C)为最小实现 N (s)D1 (s)不可简约 1/3,37/391/3,37/39261结论10.80不可简约左MFD最小实
108、现 对qp严真左MFD DL1 (s)NL(s),设n=deg detDL(s),表(A0,B0,C0)为“DL1 (s)NL(s),DL(s)行既约”的 n维观测器形实现,表(A0b,B0b,C0b)为“DL1 (s)NL(s),DL(s)列既约”的 n维能观测性形实现,则有(A0,B0,C0)为最小实现 DL1 (s)NL(s)不可简约 (A0b,B0b,C0b) 为最小实现 DL1 (s)NL(s)不可简约 结论10.81不可简约左MFD最小实现 对qp严真左MFD DL1 (s)NL(s),DL(s) 行或列既约,表 其任意形式的n维实现,n=deg detDL(s),则有 为最小实现
109、 DL1 (s)NL(s)不可简约 结论10.82狭义惟一性尽管严真不可简约右MFD 或严真不可简约左MFD 的最小实现为不惟一,但其特定形式最小实现则为惟一,如控制器形最小实现、观测器形最小实现、能控性形最小实现和能观测性形最小实现等。 结论10.83不惟一性对严真传递函数矩阵G(s),由不可简约MFD的不惟一性所决定,上述基于MFD的特定形式最小实现也为不惟一。结论10.84维数惟一性对严真传递函数矩阵G(s),不管表为哪种类型的不可简约MFD,也不管导出的为哪种类型的最小实现,最小实现的维数均为相同,且有最小实现维数MFD分母矩阵行列式的次数 2/3,38/392/3,38/39262结
110、论10.85代数等价性 对严真传递函数矩阵G(s)或矩阵分式描述MFD,其各种形式的最小实现之间为代数等价。结论10.86确定最小实现途径 对严真可简约 MFD,确定最小实现的途径可有频率方法和时间域方法两类。频率途径为:严真可简约MFD,分母矩阵为列既约或行既约=导出不可简约MFD,分母矩阵列既约或行既约=导出“控制器形实现/能控性形实现”或“观测器形实现/能观测器性形实现”=所得实现为最小实现,且维数等于分母矩阵行列式的次数时间域途径为:严真可简约MFD,分母矩阵为列既约或行既约 =导出能控能观测部分(Ac0,Bc0,Cc0) =导出能观测能控部分(Ac0,Bc0,Cc0)=最小实现即为(
111、Ac0,Bc0,Cc0)3/3,39/393/3,39/39263第第11章线性时不变系统的多项式矩阵描述章线性时不变系统的多项式矩阵描述 11.1 多项式矩阵描述多项式矩阵描述(polynomial matrix descriptions)简称为PMD,是对线性时不变系统引入的具有更广普遍性的一类内部描述 多项式矩阵描述的形式现在,推广讨论一般形式的多输入多输出线性时不变系统,定义 那么,可以导出系统的多项式矩阵描述为 PMD和其他描述的关系结论11.1PMD的传递函数矩阵 对线性时不变系统,由给出的PMD的传递函数矩阵G(s)为 G(s)=R(s)P-1(s)Q(s)+W(s) 1/6,1
112、/221/6,1/22264结论11.2状态空间描述的PMD 给定线性时不变系统的状态空间描述: 其中,E(p)为多项式矩阵,p=d/dt为微分算子,x(0)=0,且E(p)的存在反映系统的非真性。那么,状态空间描述的等价的PMD为 其中, 为n1广义状态,PMD的各个系数矩阵为 P(s)=(sI-A),Q(s)=B,R(s)=C,W(s)=E(s) 结论11.3MFD的PMD 给定qp线性时不变系统的右MFD N(s)D-1(s)+E(s)和左MFD DL-1(s)NL(s)+E(s),其中N(s)D-1(s)和DL-1(s)NL(s)为严真MFD,E(s)为多项式矩阵。那么,等价于N(s)
113、D-1(s)+E(s)的PMD为 其中, 为p1广义状态,PMD的各个系数矩阵为 P(s)=D(s), Q(s)=I, R(s)=N(s), W(s)=E(s) 2/6,2/222/6,2/22265等价于DL-1(s)NL(s)+E(s)的PMD为 其中, 为q1广义状态,PMD的各个系数矩阵为 P(s)=DL(s), Q(s)=NL(s), R(s)=I, W(s)=E(s) 不可简约PMD 不可简约PMD是线性时不变系统的最为基本和应用最广的一类PMD。 定义11.1不可简约PMD 称(P(s),Q(s),R(s),W(s)为不可简约PMD,当且仅当 P(s),Q(s)左互质, P(s)
114、,R(s)右互质 把可简约PMD化为不可简约PMD是复频率域方法中经常面临的一个问题。 3/6,3/223/6,3/22266情形 P(s),R(s)右互质, P(s),Q(s)非左互质 结论11.5构造不可简约PMD 对“ P(s),R(s)右互质, P(s),Q(s)非左互质”型可简约PMD,表mm多项式矩阵H(s)为非左互质 P(s),Q(s)的任一最大左公因子,再取 则可简约PMD的一个不可简约PMD为 情形 P(s),R(s)非右互质, P(s),Q(s)左互质 结论11.6构造不可简约PMD 对“ P(s),R(s)非右互质, P(s),Q(s)左互质”型可简约PMD,表mm多项式
115、矩阵F(s)为右互质 P(s),R(s)的任一最大右公因子,再取 即有 则可简约PMD的一个不可简约PMD为 4/6,4/224/6,4/22267情形 P(s),R(s)非右互质, P(s),Q(s)非左互质 结论11.7构造不可简约PMD 对“ P(s),R(s)非右互质, P(s),Q(s)非左互质”型可简约PMD,表mm多项式矩阵H(s)为非左互质 P(s),Q(s)的任一最大左公因子,取 mm多项式矩阵 为 的任一最大右公因子,取 则可简约PMD的一个不可简约PMD为 5/6,5/225/6,5/22268结论11.8不可简约PMD不唯一性 设(P(s),Q(s),R(s),W(s)
116、为线性时不变系统的一个不可简约PMD,P(s)为mm多项式矩阵,Q(s)、R(s)和W(s)为mp、qm和qp多项式矩阵。表U(s)和V(s)为任意两个mm单模阵,取则 也为系统的一个不可简约PMD 6/6,6/226/6,6/2226911.2 多项式矩阵描述的状态空间实现 PMD的实现 考虑线性时不变系统,其多项式矩阵描述即PMD为 其中,P(s)为mm多项式矩阵;Q(s)、R(s)和W(s)为mp、qm和qp多项式矩阵。定义11.2PMD的实现 称状态空间描述为PMD (P(s),Q(s),R(s),W(s))的一个实现,如果两者的传递函数矩阵为相等,即成立:R(s)P-1(s)Q(s)
117、+W(s)=C(sI-A)-1B+E(s) 其中,E(s)=E(p)p=s 注 PMD的实现具有强不唯一性,即不仅实现的结果不唯一,且实现的维数也不唯一。 1/1,7/221/1,7/2227011.3 多项式矩阵描述的互质性和状态空间描述的能控性与能观测性 左互质性与能控性 考虑线性时不变系统,其多项式矩阵描述为 其中,P(s)为mm多项式矩阵,Q(s)、R(s)和W(s)为mp、qm和qp多项式矩阵。系统的状态空间描述即PMD的一个实现为其中,A为nn常阵,B和C为np和qn常阵,E(p)为qp多项式矩阵。下面,给出能控性和左互质性间关系的结论 结论11.16左互质性和能控性 对线性时不变
118、系统的PMD及其状态空间实现,有(P(s),Q(s))左互质=(A,B)完全能控 结论11.17右互质性和能观测性 对线性时不变系统的PMD及其状态空间实现有(P(s),R(s))右互质=(A,C)完全能观测 1/4,8/221/4,8/22271结论11.18不可简约PMD的最小描述性 对线性时不变系统,如同称(A,B)完全能控和(A,C)完全能观测的状态空间描述(A,B,C,E(p))为最小描述一样,也称(P(s),Q(s)左互质和(P(s),R(s)右互质的不可简约PMD(P(s),Q(s),R(s),W(s)为最小描述。 结论11.19 MFD右互质性和能观测性 考虑线性时不变系统的右
119、MFD 为严真,其能控类实现为 其中,dim(Ac)=deg detD(s)。则有 D(s),N(s)右互质=(Ac,Cc)完全能观测 2/4,9/222/4,9/22272结论11.20 MFD左互质性和能控性 考虑线性时不变系统的左MFD 为严真,其能观测类实现为 其中,dim(A0)=deg detDL(s)。则有 DL-1(s),NL(s)左互质=(A0,B0)完全能控 结论11.21状态空间描述的互质性 考虑线性时不变系统,其状态空间描述为(A,B,C,E(p)),传递函数矩阵G(s)的关系式为 G(s)=C(sI-A) -1 B+E(s)=R(s)P-1(s)Q(s)+W(s) 则
120、由PMD左右互质性和状态空间描述能控性能观测性的等价关系,可知 sI-A,B左互质=(A,B)完全能控 sI-A,C右互质=(A,C)完全能观测 3/4,10/223/4,10/22273结论11.22 SISO系统互质性 考虑单输入单输出即SISO线性时不变系统,表其传递函数g(s)为 其中,P(s)为mm多项式矩阵,r(s)和q(s)为1m和m1多项式项量,W(s)为多项式,(s)为P(s)的最小多项式。则有 P(s),r(s)右互质=(s)和r(s)H(s)不含相消因子 P(s),q(s)左互质=(s)和H(s)q(s)不含相消因子 P(s),r(s)和P(s),q(s)均互质=g(s)
121、严真部分不含零点极点对消 4/4,11/224/4,11/2227411.4 传输零点和解耦零点 PMD的极点 考虑线性时不变系统的PMD(P(s),Q(s),R(s),W(s),其传递函数矩阵为 G(s)=R(s)P-1(s)Q(s)+W(s) 定义11.3 PMD的极点 对PMD(P(s),Q(s),R(s),W(s),定义: PMD的极点“R(s)P-1(s)Q(s)+W(s)”的极点 结论11.23 PMD的极点 表(A,B,C,E(p))为PMD(P(s),Q(s),R(s),W(s)的一个最小实现,则有PMD的极点“det(sI-A)=0”的根 结论11.24 PMD的极点 若PM
122、D(P(s),Q(s),R(s),W(s)为不可简约,则有PMD的极点“detP(s)=0的根” 定义11.4 PMD的传输零点 对PMD(P(s),Q(s),R(s),W(s),其传递函数矩阵为G(s)=R(s)P-1(s)Q(s)+W(s),则定义: PMD的传输零点“R(s)P-1(s)Q(s)+W(s)”的零点 1/2,12/221/2,12/22275结论11.25 PMD的传输零点 表(A,B,C,E(p))为PMD(P(s),Q(s), R(s),W(s)的任一最小实现,则有 PMD的传输零点使 降秩的s值 结论11.26 PMD的传输零点 若PMD(P(s),Q(s),R(s)
123、,W(s)为不可简约,则有PMD的传输零点使 降秩的s值 2/2,13/222/2,13/2227611.5 系统矩阵 考虑线性时不变系统,其多项式矩阵描述为其中,P(s)为mm非奇异多项式矩阵;Q(s)、R(s)和W(s)为mp、qm和qp多项式矩阵。进而,表上式为增广变量方程形式,有定义11.8 PMD系统矩阵 线性时不变系统PMD的系统矩阵定义为其增广变量方程(11.143)的系数矩阵,即 结论11.35 状态空间描述系统矩阵 线性时不变系统状态空间描述的系统矩阵为 1/4,14/221/4,14/22277结论11.36 MFD系统矩阵 对qp线性时不变系统的MFD,右N(s)D-1(
124、s)的系统矩阵为 左DL -1 (s)NL(s)的系统矩阵为 结论11.37判断不可简约性 线性时不变系统PMD的系统矩阵为S(s),有 PMD不可简约=S(s)的前m行和前m列分别满秩, 结论11.38 PMD的极点零点 线性时不变系统PMD的系统矩阵为S(s),若PMD为不可简约,则有PMD的极点使S(s)左上方mm块矩阵降秩s值 PMD的传输零点使S(s)降秩s值 2/4,15/222/4,15/22278增广系统矩阵通常,一个线性时不变系统的不同类型描述的系统矩阵在维数上为不同。进而,同一类型不同描述的系统矩阵在维数上也常为不同。增广系统矩阵正是为克服由此而引起的不便而在系统矩阵基础上
125、导出的一类广义系统矩阵。定义11.9 PMD增广系统矩阵 线性时不变系统PMD的增广系统矩阵定义为 其中,为正整数且可按需要任取 结论11.42 不可简约性相同 对线性时不变系统的系统矩阵S(s)及其增广系统矩阵Se(s),有Se(s)不可简约=S(s)不可简约 结论11.43互质性相同 对线性时不变系统的系统矩阵S(s)及其增广系统矩阵Se(s),有Pe(s),Qe(s)左互质=P(s),Q(s)左互质 Pe(s),Re(s)右互质=P(s),R(s)右互质 3/4,16/223/4,16/22279结论11.44 极点和传输零点相同 对线性时不变系统的不可简约系统矩阵S(s)及其增广系统矩
126、阵Se(s),有Se(s)的极点S(s)的极点 Se(s)的传输零点S(s)的传输零点 结论11.45解偶零点相同 对线性时不变系统的可简约系统矩阵S(s)及其增广系统矩阵Se(s),有Se(s)的输入解偶零点S(s)的输入解偶零点 Se(s)的输出解偶零点S(s)的输出解偶零点 结论11.46传递函数矩阵相同 线性时不变系统的系统矩阵S(s)及其增广系统矩阵Se(s)具有相同的传递函数矩阵,即有Re(s)Pe-1(s)Qe(s)+W(s)=R(s)P-1(s)Q(s)+W(s) 结论11.47分母矩阵行列式相同 线性时不变系统的系统矩阵S(s)及其增广系统矩阵Se(s)具有相同的分母矩阵行列
127、式,即有detPe(s)=detP(s) 结论11.48特性关系属性相同 对线性时不变系统,引入增广系统矩阵Se(s)代替系统矩阵S(s)以讨论不同描述间关系,不会损失不同描述在特性上的关系属性,如互质性、能控性能观测性、稳定性等。4/4,17/224/4,17/2228011.6 严格系统等价 对线性时不变系统,考虑相同输入和相同输出的两个PMD的系统矩阵S1(s)和S2(s),它们既可属于同一系统也可属于不同系统,并表S1(s) 和S2(s)分别为 其中,Pi(s)为mimi非奇异多项式矩阵,Ri(s)、Qi(s)和Wi(s)为mip、qmi和qp多项式矩阵,i=1,2。进而,不妨设m1m
128、2m。 定义11.10严格系统等价 称两个PMD型系统矩阵S1(s) 和S2(s)为严格系统等价,当且仅当存在mm单模阵U(s)和V(s),以及qm和mp多项式矩阵X(s)和Y(s),使成立: 并且,记为S1(s)S2(s) 1/5,18/221/5,18/22281三点说明:1:严格系统等价式一种变换关系2:严格系统等价变换式一类特定的左右单模变换3:严格系统等价变换满足对称性、自反性和传递性对称性:若S1(s)S2(s),则S2(s)S1(s)。自反性:S1(s)S1(s)。传递性:若S1(s)S2(s),S2(s)S3(s),则S1(s)S3(s)。严格系统等价变换的性质 线性时不变系统
129、的两个PMD型系统矩阵S1(s)和S2(s),若S1(s)S2(s)即严格系统等价,则两者分母矩阵P2(s)和P1(s)具有等同的不变多项式,即有detP2(s)=0detP1(s)其中,0为非零常数 严格系统等价变换下传递函数矩阵保持不变 2/5,19/222/5,19/22282对线性时不变系统,表两个多项式矩阵描述 其系统矩阵为S1(s)和S2(s),再令(A1,B1,C1,E1(p))PMD1的任一能控类或能观测类实现 (A2,B2,C2,E2(p))PMD2的任一能控类或能观测类实现若S1(s)S2(s)即严格系统等价,则两个同类实现具有相同维数和相同特征多项式,即有 dim(A1)
130、= dim(A2) det(sI-A1)=det(sI-A2) 严格系统等价变换下系统同类实现在维数和特征多项式上的等同性 3/5,20/223/5,20/22283左互质性和右互质性在严格系统等价变换下的不变性 对线性时不变系统,PMD的互质性在严格等价变换下保持不变 若S1(s)S2(s)即严格系统等价,则有 P2(s),Q2(s)左互质=P1(s),Q1(s)左互质 P2(s),R2(s)右互质=P1(s),R1(s)右互质 能控性和能观测性在严格系统等价变换下的不变性“状态空间描述代数等价”和“系统矩阵严格系统等价”的等价性 对线性时不变系统,表两个状态空间描述为(A1,B1,C1,E
131、1(p))和(A2,B2,C2,E2(p)) 系统矩阵为 则有(A2,B2,C2,E2(p))代数等价(A1,B1,C1,E1(p))=S2(s)S1(s) 4/5,21/224/5,21/22284传递函数矩阵的所有不可简约MFD的严格系统等价 对线性时不变系统的qp传递函数矩阵G(s),且不要求为严真,则G(s)的所有不可简约MFD必都为严格系统等价 结论11.57不可简约PMD的严格系统等价 对线性时不变系统的qp传递函数矩阵G(s),G(s)的所有不可简约MFD为严格系统等价 严格系统等价描述在结构性质和运动行为上的等同性由严格系统等价性保证,在不可简约的前提下,线性时不变系统的三类描
132、述即状态空间描述、右或左MFD以及PMD在用于系统的分析和综合时的结果为完全等价,不会出现丢失系统结构信息的情况。5/5,22/225/5,22/22285第12章线性时不变系统的复频率域分析12.1并联系统的能控性和能观测性 并联系统并联系统是以“输入相同”和“输出相加”为特征的一类组合系统 首先,给出对子系统的两个基本假定。一是,S1和S2可由其传递函数矩阵G1(s)和G2(s)完全表征,即其相应的状态空间描述为完全能控和完全能观测。二是,子系统传递函数矩阵Gi(s),i=1,2为qipi有理分式矩阵。且表为不可简约右和左MFD:Gi(s)=Ni(s)Di-1(s)=DLi-1(s)NLi
133、(s),i=1,2 S1S2u=u1=u2,y=y1+y2 p1=p2=p,q1=q2=q 结论12.1能控性条件 由线性时不变子系统S1和S2组成的并联系统Sp,若取G1(s)=不可简约右MFD N1(s)D1-1(s)G2(s)=不可简约右MFD N2(s)D2-1(s)则有Sp完全能控D1(s),D2(s)左互质 1/2,1/151/2,1/15286结论12.2能观测性条件 线性时不变子系统S1和S2组成的并联系统Sp,若取G1(s)=不可简约左MFD DL1-1(s) NL1(s)G2(s)=不可简约左MFD DL2-1(s) NL2(s)则有Sp完全能观测DL1(s),DL2(s)
134、右互质 结论12.3不可简约性条件 由线性时不变子系统S1和S2组成的并联系统Sp,若取G1(s)和G2(s)为“不可简约右MFD N1(s)D1-1(s)与N2(s)D2-1(s)”和“不可简约左MFD DL1-1(s) NL1(s)与DL2-1(s)NL2(s)”,则有Sp不可简约,即可用G1(s)+G2(s)完全表征D1(s),D2(s)左互质,DL1(s),DL2(s)右互质 结论12.4能控性和能观测性条件 由线性时不变子系统S1和S2组成的多输入多输出并联系统Sp,则Sp保持完全能控和完全能观测的一个充分条件是,qp传递函数矩阵G1(s)和G2(s)不包含公共极点。结论12.6能控
135、性和能观测性条件 由线性时不变子系统S1和S2组成的单输入单输出并联系统Sp,则Sp保持为完全能控和完全能观测的充分必要条件是,标量传递函数g1(s)和g2(s)不包含公共极点。 2/2,2/152/2,2/1528712.2串联系统的能控性和能观测性串联系统是由子系统按串联方式顺序联接的组合系统 首先,对子系统引入两个基本假定。一是,S1和S2可由其传递函数矩阵G1(s)和G2(s)所完全表征,即其状态空间描述为完全能控和完全能观测。二是,Gi(s),i=1,2,为qipi有理分式矩阵,且表为不可简约右和左MFD:Gi(s)=Ni(s)Di-1(s)=DLi-1NLi(s),i=1,2 S1
136、 S2uu1y1u2y2y进而,由子系统的S1S2串联特征,可以给出系统组成上的相应约束条件为 u=u1,y1=u2,y=y2 p1=p,q1=p2,q2=q 结论12.7能控性条件 由线性时不变子系统S1和S2组成的串联系统ST,若取G1(s)=不可简约右MFD N1(s)D1-1(s)G2(s)=不可简约右MFD N2(s)D2-1(s)则有ST完全能控D2(s),N1(s)左互质 1/4,3/151/4,3/15288结论12.8能控性条件 由线性时不变子系统按S1S2顺序组成的串联系统ST,若取G1(s)=不可简约右MFD N1(s)D1-1(s)G2(s)=不可简约左MFD DL2-
137、1(s)NL2(s)则有ST完全能控 DL2(s),NL2(s)N1(s)左互质 结论12.9能控性条件 由线性时不变子系统按S1S2顺序组成的串联系统ST,若取G1(s)=不可简约左MFD DL1-1(s)NL1(s)G2(s)=不可简约右MFD N2(s)D2-1(s) 则有ST完全能控 DL1(s) D2(s),NL1(s)左互质 结论12.10能观测性条件 由线性时不变子系统按S1S2顺序组成的串联系统ST,若取G1(s)=不可简约左MFD DL1-1(s) NL1(s)G2(s)=不可简约左MFD DL2-1(s) NL2(s)则有ST完全能观测DL1(s),NL2(s)右互质 2/
138、4,4/152/4,4/15289结论12.11能观测性保持条件 由线性时不变子系统按S1S2顺序组成的串联系统ST,若取G1(s)=不可简约左MFD DL1-1(s) NL1(s)G2(s)=不可简约右MFD N2(s)D2-1(s) 则有ST完全能观测DL1(s) D2(s),N2(s)右互质 结论12.12能观测性保持条件 由线性时不变子系统按S1S2顺序组成的串联系统ST,若取G1(s)=不可简约右MFD N1(s) D1-1(s)G2(s)=不可简约左MFD DL2-1(s)NL2(s) 则有ST完全能观测D1(s),NL2(s)N1(s)右互质 结论12.13能控性条件 由线性时不
139、变子系统按S1S2顺序组成的多输入多输出串联系统ST,设p=p1q1=p2,传递函数矩阵G1(s)为满秩,则ST保持完全能控的一个充分条件是,没有G2(s)极点等同于G1(s)传输零点 3/4,5/153/4,5/15290结论 12.15能观测性条件 由线性时不变子系统按S1S2顺序组成的多输入多输出串联系统ST,设p2=q1q2=q,传递函数矩阵G2(s)为满秩,则ST保持完全能观测的一个充分条件是,没有G1(s)极点等同于G2(s)传输零点。结论12.17能控性条件 由线性时不变子系统按S1S2顺序组成的单输入单输出串联系统ST,则ST保持完全能控的充分必要条件是,没有g2(s)极点为g
140、1(s)零点所对消。结论12.18能观测性条件 由线性时不变子系统按S1S2顺序组成的单输入单输出串联系统ST,则ST保持完全能观测的充分必要条件是,没有g1(s)极点为g2(s) 零点所对消。结论12.19完全表征条件 由线性时不变子系统按S1S2顺序组成的单输入单输出串联系统ST,则ST可用g2(s)g1(s)完全表征的充分必要条件是,g1(s)和g2(s)没有极点零点对消现象。 4/4,6/154/4,6/1529112.3状态反馈系统的能控性和能观测性 结论12.21状态反馈系统复频率域形式 对线性时不变受控系统,状态反馈系统复频率域结构基本形式如图所示,并被采用作为复频率域方法中分析
141、和综合状态反馈的基本模型。 Dhc-1(s)S-1(s)NLcDhc-1DLcKD-1(s)N(s)K(s)结论12.22状态反馈系统的MFD 复频率域结构图表征的线性时不变状态反馈系统,闭环传递函数矩阵的右MFD为GK(S)=N(s)DK-1(s)闭环分母矩阵DK(s)为DK(s)=DhcS(s)+(DLc+K)(s) 1/2,7/151/2,7/15292结论12.23能控性 复频率域结构图表征的线性时不变状态反馈系统K和开环受控系统0有K完全能控0完全能控 结论12.24能观测性 复频率域结构图表征的线性时不变状态反馈系统K和开环受控系统0,0完全观测不能保证K必为完全能观测。2/2,8
142、/152/2,8/1529312.4输出反馈系统的能控性和能观测性考虑图示结构组成的线性时不变输出反馈系统F。首先对输出反馈系统F引入三个基本约定:(i)子系统S1和S2为真或严真,且可由传递函数矩阵G1(s)和G2(s)分别完全表征。(ii)为保证输出反馈系统F传递函数矩阵的真性或严真性,令detI+G1(s)G2(s)s= =detI+G2(s)G1(s)s= 0 (iii) 对输出反馈系统F中包含的串联系统,表 进而,由输出反馈的联接特征,可以导出系统组成上的约束关系式为 u1=uy2,y=y1=u2 S1S2结论12.27F的传递函数矩阵 对线性时不变输出反馈系统F,有GF(s)=G1
143、(s)I+G2(s)G1(s) -1 GF(s)=I+G1(s)G2(s) -1 G1(s) 1/2,9/151/2,9/15294结论12.28完全能控条件 线性时不变输出反馈系统F,有 F完全能控S12完全能控 结论12.29完全能观测条件 线性时不变输出反馈系统F,有F完全能观测S21完全能观测 结论12.32联合能控能观测条件 单输入单输出线性时不变输出反馈系统F,有 F完全能控和完全能观测g2(s)极点和g1(s)零点间不存在对消 结论12.33能控和能观测条件 多输入多输出线性时不变输出反馈系统F,若对子系统S2有G2(s)=F(常阵),则有F完全能控S1完全能控 F完全能观测S1
144、完全能观测 2/2,10/152/2,10/1529512.5直接输出反馈系统的稳定性分析系统运动的稳定性可以区分为两种基本类型,“内部稳定性”和“外部稳定性”,前者就是渐进稳定性,后者即为有界输入有界输出稳定性或BIBO稳定性 考虑图所示的线性时不变直接输出反馈系统DF,其特点是反馈通道中子系统的传递函数矩阵G2(s)=I。通常,也称直接输出反馈系统为单位输出反馈系统。进而。约定子系统S1可由传递函数矩阵G1(s)完全表征,G1(s)为方的真有理分式矩阵,detI+G1(s)s= 0。G1(s)结论12.34 两类稳定等价性 线性时不变直接输出反馈系统DF,子系统S1可由G1(s)完全表征,
145、则有DFBIBO稳定DF渐近稳定 1/2,11/151/2,11/15296结论12.35 特征多项式 线性时不变直接输出反馈系统DF,子系统S1由G1(s)完全表征,1(s)为G1(s)的特征多项式,则有DF的特征多项式1(s) detI+G1(s) 其中,为非零常数 结论12.36 稳定条件 线性时不变直接输出反馈系统DF,若G1(s)以有理分式矩阵表征,则有DF为渐近稳定和BIBO稳定“1(s) detI+G1(s)0根”均具有负实部 其中,1(s)为G1(s)的特征多项式。结论12.37 稳定条件 线性时不变直接输出反馈系统DF,若G1(s)以不可简约右MFD N1(s)D1-1(s)
146、表征,则有DF为渐近稳定和BIBO稳定“detD1(s)+N1(s)0根”均具有负实部 结论12.38 稳定条件 线性时不变直接输出反馈系统DF,若G1(s)以不可简约左MFD DL1-1(s)NL1(s)表征,则有DF为渐近稳定和BIBO稳定“detDL1(s)+NL1(s)0根”均具有负实部 2/2,12/152/2,12/1529712.6具有补偿器的输出反馈系统的稳定性分析 具有补偿器的输出反馈系统简表为CF。CF是更为一般的形式的输出反馈系统,其构成如图所示。G1(s)G2(s)同样,对CF引入如下三个基本约定子系统S1和S2为完全能控和完全能观测,即可由其传递函数矩阵G1(s)和G
147、2(s)完全表征,且G1(s)和G2(s)为qp和pq真有理分式矩阵。子系统S1和S2采用有理分式矩阵和不可简约MFD描述:Gi(s)=不可简约右MFD Ni(s)Di-1(s)=不可简约左MFD DLi-1(s)NLi(s) 为保证CF的传递函数矩阵GCF(s)为真,设detI+G1(s)G2(s)s= 0 1/3,13/151/3,13/15298结论12.39 等价条件 具有补偿器的线性时不变输出反馈系统CF,表 若满足“S12完全能控,S21完全能观测”条件,则有CF渐近稳定CF BIBO稳定 若不满足“S12完全能控,S21完全能观测”条件,则有CF渐近稳定=CF BIBO稳定 结论
148、12.40特征多项式 具有补偿器的线性时不变输出反馈系统CF,表CF状态空间描述ACF,BCF,CCF,ECFS1状态空间描述A1,B1,C1,E1S2状态空间描述A2,B2,C2,E2则有CF的多项式det(sIACF)1(s)2(s)detI+G1(s)G2(s) 结论12.41渐近稳定条件 具有补偿器的线性时不变输出反馈系统CF,若G1(s)和G2(s)以有理分式矩阵表征,则有CF渐近稳定“1(s)2(s)detI+G1(s)G2(s)0根”均具有负实部 2/3,14/152/3,14/15299结论12.42渐近稳定条件 具有补偿器的线性时不变输出反馈系统CF,若G1(s)和G2(s)
149、以不可简约左MFD DL1-1(s) NL1(s)和不可简约右MFD N2(s)D2-1(s)表征,则有CF渐近稳定“detDL1(s)D2(s)+NL1(s)N2(s)0根”均具有负实部 结论12.43渐进稳定条件 具有补偿器的线性时不变输出反馈系统CF,若G1(s)和G2(s)以不可简约右MFD N1(s)D1-1(s)和不可简约左MFD DL2-1(s) NL2(s)表征,则有CF渐近稳定“detDL2(s)D1(s)+NL2(s)N1(s)0根”均具有负实部 结论12.45特殊CF稳定条件 具有补偿器的线性时不变输出反馈系统CF,若子系统S2的传递函数矩阵G2(s)F(常阵),则由结论12.41结论12.43给出的CF为渐近稳定的条件也是CF为BIBO稳定的充分必要条件 结论12.41渐近稳定条件 具有补偿器的线性时不变输出反馈系统CF,若G1(s)和G2(s)以有理分式矩阵表征,则有CF渐近稳定“1(s)2(s)detI+G1(s)G2(s)0根”均具有负实部 3/3,15/153/3,15/15300