现代控制理论第一章02
现代控制理论课件2
38
二、从系统的机理出发建立状态空间表达式
例1、求图示机械系统的状态空间表达式
外力 u(t)
K ---弹性系数 m
牛顿力学定律 my u by ky
阻 尼 系 数
y(t) b
位移 令
b u(t ) ky m y y
x1 y
x2 y
39
动态方程如下
x1 x2
x1 y 1 0 x2
41
例:设有如图所示的机 械系统。它由两个彼 此耦合的平台构成。 并借助于弹簧和阻尼 到达地基。试选择合 适的状态变量,写出 该系统的状态空间模 型。
42
解答:依题意,进行受力分析,可得如下的微分方程:
M1y1 = u -k1 (y1 - y 2 )-f1 (y1 - y 2 ) M2y 2 = k1 (y1 - y 2 ) + f1 (y1 - y 2 )-k 2y 2 -f 2y 2
其中: a11 a12 a1n a a22 a2 n 21 A — 系统内部状态的联系, an1 an 2 ann
18
称为系统矩阵 , 为n n方阵;
多输入——多输出定常系统: 用向量矩阵表示时的状态空间表达式为:
Ax Bu x y Cx Du
其状态变量为: x1 , x2 ,, xn , 则状态方程的一般形式 为:
1 a11x1 a12 x2 a1n xn b11u1 b12u2 b1r ur x 2 a21x1 a22 x2 a2 n xn b21u1 b22u2 b2 r ur x n an1 x1 an 2 x2 ann xn bn1u1 bn 2u2 bnr ur x
现代控制理论基础第一章
Elements of Modern Control Theory主讲:董霞现代控制理论基础西安交通大学机械工程学院Email:xdong@办公地点:西二楼东207参考教材《现代控制工程》王军平董霞主编西安交通大学出版社教材《现代控制理论基础》(机械类)何钺编机械工业出版社《现代控制工程》(第三版)Katsuhiko Ogata著卢伯英、于海勋译电子工业出版社第一章绪论现代控制理论是在20世纪50年代末、60年代初形成的控制理论。
之所以称其为现代控制理论是与经典控制理论相比较而言的。
1.1 控制理论发展简史目前国内外学术界普遍认为控制理论经历了三个发展阶段:经典控制理论现代控制理论智能控制理论这种阶段性发展是由简单到复杂、由量变到质变的辩证发展过程。
并且,这三个阶段不是相互排斥,而是相互补充、相辅相成的,它们各有其应用领域,并还在不同程度地继续发展着。
控制理论中反馈的概念代表性人物:瓦特(J.Watt),于1788年发明了蒸汽机飞球调速器。
这是一个典型的自动调节系统,由此拉开了经典控制理论发展的序幕。
控制理论诞生前,人们对于反馈就有了认识。
经典控制理论的诞生1868年,英国物理学家J.C.Maxwell 发表《论调速器》论文,解决了蒸汽机调速系统中出现的剧烈振荡问题;1877年,英国科学家E.J. Routh 建立了劳斯稳定性判据;1895年,德国数学家A. Hurwitz 提出了胡尔维茨稳定性判据;1892年,俄国数学家A. M.Lyapunov 发表了专著《论运动稳定性的一般问题》;1922年,美国的N. Minorsky 研究出用于船舶驾驶的伺服机构并提出PID 控制方法;1932年,美籍瑞典人H. Nyquist 提出了频域内研究系统稳定性的频率判据;经典控制理论的诞生1940年,H. W.Bode引入了对数坐标,使频域稳定性判据更适合工程应用;1942年,H. Harris引入了传递函数概念;1948年,W.R. Evans提出了根轨迹方法;1948年,N. Wiener发表了著名的《控制论》,标志着经典控制理论的诞生。
现代控制理论习题解答(第一章)
Ra
La
i f = 常数
ua
f ia D J
ω
ML
【解】: 设状态变量为:
题 1-2 图
⎡ x1
⎢ ⎣
x
2
⎤ ⎥ ⎦
=
⎡ia ⎢⎣ω
⎤ ⎥ ⎦
其中 ia 为流过电感上的电流, ω 电动机轴上的角速度。 电动机电枢回路的电压方程为:
eb 为电动机反电势。 电动机力矩平衡方程为
•
ua = La ia + Ra ⋅ ia + eb
(4) y (4) + 3y + 2y = −3u + u
【解】:
5
在零初始条件下,方程两边拉氏变换,得到传递函数,再根据传递函数求状态空间 表达式。 此题多解,一般写成能控标准型、能观标准型或对角标准型,以下解法供参考。 (1)传递函数为:
状态空间表达式为:
G(s) =
2
s3 + 2s2 + 4s + 6
1⎤
R 2 C1 −1
R2C2
⎥ ⎥ ⎥
⎡ ⎢ ⎣
x1 x2
⎥⎦
⎤ ⎥ ⎦
+
⎡ ⎢ ⎢ ⎣
1
R1C1 0
⎤ ⎥⎥u i ⎦
y = u0 = [0
1]⎢⎡
⎣
x1 x2
⎤ ⎥ ⎦
(2)
设状态变量: x1 = iL 、 x2 = uc 而
1
根据基尔霍夫定律得: 整理得
•
iL = C uc
•
ui = R ⋅ iL + LiL + uc
•
M D = J ω + fω + M L
《现代控制理论》课后习题全部答案(最完整打印版)
第一章习题答案1-1试求图1-27系统的模拟结构图,并建立其状态空间表达式。
解:系统的模拟结构图如下:系统的状态方程如下:阿令,则所以,系统的状态空间表达式及输出方程表达式为状态变量的状态方程,和以电阻上的电压作为输出量的输出方程。
解:由图,令,输出量有电路原理可知:既得写成矢量矩阵形式为:1-3参考例子1-3(P19).1-4两输入,,两输出,的系统,其模拟结构图如图1-30所示,试求其状态空间表达式和传递函数阵。
解:系统的状态空间表达式如下所示:1-5系统的动态特性由下列微分方程描述列写其相应的状态空间表达式,并画出相应的模拟结构图。
解:令,则有相应的模拟结构图如下:1-6(2)已知系统传递函数,试求出系统的约旦标准型的实现,并画出相应的模拟结构图解:1-7给定下列状态空间表达式(1)画出其模拟结构图(2)求系统的传递函数解:(2)1-8求下列矩阵的特征矢量(3)解:A的特征方程解之得:当时,解得:令得(或令,得)当时,解得:令得(或令,得)当时,解得:令得1-9将下列状态空间表达式化成约旦标准型(并联分解)(2)解:A的特征方程当时,解之得令得当时,解之得令得当时,解之得令得约旦标准型1-10已知两系统的传递函数分别为W1(s)和W2(s)试求两子系统串联联结和并联连接时,系统的传递函数阵,并讨论所得结果解:(1)串联联结(2)并联联结1-11(第3版教材)已知如图1-22所示的系统,其中子系统1、2的传递函数阵分别为求系统的闭环传递函数解:1-11(第2版教材)已知如图1-22所示的系统,其中子系统1、2的传递函数阵分别为求系统的闭环传递函数解:1-12已知差分方程为试将其用离散状态空间表达式表示,并使驱动函数u的系数b(即控制列阵)为(1)解法1:解法2:求T,使得得所以所以,状态空间表达式为第二章习题答案2-4用三种方法计算以下矩阵指数函数。
(2)A=解:第一种方法:令则,即。
求解得到,当时,特征矢量由,得即,可令当时,特征矢量由,得即,可令则,第二种方法,即拉氏反变换法:第三种方法,即凯莱—哈密顿定理由第一种方法可知,2-5下列矩阵是否满足状态转移矩阵的条件,如果满足,试求与之对应的A阵。
《现代控制理论》课后习题全部答案(最完整打印版)
《现代控制理论》课后习题全部答案(最完整打印版)第⼀章习题答案1-1试求图1-27系统的模拟结构图,并建⽴其状态空间表达式。
解:系统的模拟结构图如下:系统的状态⽅程如下:阿令,则所以,系统的状态空间表达式及输出⽅程表达式为状态变量的状态⽅程,和以电阻上的电压作为输出量的输出⽅程。
解:由图,令,输出量有电路原理可知:既得写成⽮量矩阵形式为:1-4两输⼊,,两输出,的系统,其模拟结构图如图1-30所⽰,试求其状态空间表达式和传递函数阵。
解:系统的状态空间表达式如下所⽰:1-5系统的动态特性由下列微分⽅程描述列写其相应的状态空间表达式,并画出相应的模拟结构图。
解:令,则有相应的模拟结构图如下:1-6(2)已知系统传递函数,试求出系统的约旦标准型的实现,并画出相应的模拟结构图解:1-7给定下列状态空间表达式(1)画出其模拟结构图(2)求系统的传递函数解:(2)1-8求下列矩阵的特征⽮量(3)解:A的特征⽅程解之得:当时,解得:令得(或令,得)当时,解得:令得(或令,得)当时,解得:令得1-9将下列状态空间表达式化成约旦标准型(并联分解)(2)解:A的特征⽅程当时,解之得令得当时,解之得令得当时,解之得令得约旦标准型1-10已知两系统的传递函数分别为W1(s)和W2(s)试求两⼦系统串联联结和并联连接时,系统的传递函数阵,并讨论所得结果解:(1)串联联结(2)并联联结1-11(第3版教材)已知如图1-22所⽰的系统,其中⼦系统1、2的传递函数阵分别为求系统的闭环传递函数解:1-11(第2版教材)已知如图1-22所⽰的系统,其中⼦系统1、2的传递函数阵分别为求系统的闭环传递函数解:1-12已知差分⽅程为试将其⽤离散状态空间表达式表⽰,并使驱动函数u的系数b(即控制列阵)为(1)解法1:解法2:求T,使得得所以所以,状态空间表达式为第⼆章习题答案2-4⽤三种⽅法计算以下矩阵指数函数。
(2)A=解:第⼀种⽅法:令则,即。
现代控制理论(刘豹)第一章
状态变量
状态向量
状态空间
状态方程
状态:表征 系统运动的信 息和行为 状态变量: 能完全表示系 统运动状态的 最小个数的一 组变量
由状态变量 构成的向量 x1(t) x2(t) : xn(t)
以各状态变量 x1(t),x2(t),…… xn(t)为坐标轴 组成的几维空 间。
S nY ( s ) + an −1S n −1Y ( s ) + ... + a0Y ( s ) = bm S mu ( s ) + ... + b0Y ( s )
(bm S m + bm −1S m −1 + ... + b0 ) Y ( s ) Z ( s ) G ( s) = Y ( s) / U ( s) = = ⋅ n n −1 ( S + an −1S + ... + a0 ) Z ( s) U ( s)
& x3 x3
x2 x1
机电工程系
∫
∫
∫
习题2 习题
已知离散系统的差分方程为
y (k + 2) + 3 y (k + 1) + 2 y (k ) = 2u (k + 1) + 3u (k )
试求系统的状态空间表达式,并画出其模拟结构图。
解:假设初始条件为零,系统微分方程的 Z 变换为:
z 2Y ( z ) + 3 zY ( z ) + 2Y ( z ) = 2sU ( z ) + 3U ( z )
S n Z ( s ) + an −1S n −1Z ( s ) + ... + a0 Z ( s ) = U ( s ) Y ( s ) = bn −1S
《现代控制理论》复习提纲()
现代控制理论复习提纲第一章:绪论(1)现代控制理论的根本内容包括:系统辨识、线性系统理论、最优控制、自适应控制、最优滤波(2)现代控制理论与经典控制理论的区别第二章:控制系统的状态空间描述1.状态空间的根本概念;系统、系统变量的组成、外部描述和内部描述、状态变量、状态向量、状态空间、状态方程、状态空间表达式、输出方程2.状态变量图概念、绘制步骤;3.由系统微分方程建立状态空间表达式的建立;第三章:线性控制系统的动态分析1.状态转移矩阵的性质及其计算方法〔1〕状态转移矩阵的根本定义;〔2〕几个特殊的矩阵指数;〔3〕状态转移矩阵的根本性质〔以课本上的5个为主〕;〔4〕状态转移矩阵的计算方法掌握:方法一:定义法方法二:拉普拉斯变换法例题2-2第四章:线性系统的能控性和能观测性(1)状态能控性的概念状态能控、系统能控、系统不完全能控、状态能达(2)线性定常连续系统的状态能控性判别包括;格拉姆矩阵判据、秩判据、约当标准型判据、PBH判据掌握秩判据、PBH判据的计算(3)状态能观测性的概念状态能观测、系统能观测、系统不能观测(4)线性定常连续系统的状态能观测性判别包括;格拉姆矩阵判据、秩判据、约当标准型判据、PBH判据掌握秩判据、PBH判据的计算(5)能控标准型和能观测标准型只有状态完全能控的系统才能变换成能控标准型,掌握能控标准I型和II型的只有状态完全能观测的系统才能变换成能控标准型,掌握能观测标准I型和II 型的计算方法第五章:控制系统的稳定性分析〔1〕平衡状态〔2〕李雅普诺夫稳定性定义:李雅普诺夫意义下的稳定概念、渐进稳定概念、大范围稳定概念、不稳定性概念(3)线性定常连续系统的稳定性分析例4-6第六章线性系统的综合(1)状态反应与输出反应(2)反应控制对能控性与观测性的影响复习题1. 、和统称为系统变量。
2. 系统的状态空间描述由和组成,又称为系统的动态方程。
3. 状态变量图是由、和构成的图形。
4. 计算1001A-⎡⎤=⎢⎥⎣⎦的矩阵指数Ate__________。
第1章现代控制理论
1.1 状态变量及状态空间表达式
1.1.1 状态变量 状态变量是既足以完全确定系统运动状态而个数又是最小的一组变量,
当其在t=t0时刻的值已知时,则在给定t≥t0时刻的输入作用下,便能完全确 定系统在任何t≥t0时刻的行为。
1.1.2 状态矢量
如果 个状态变量用
表示,并把这些状态变量看作
是矢量 的分量,则 就称为状态矢量,记作:
(34)
这种形式下,状态的选取
1.4.3 多输入一多输出系统微分方程的实现 一双输入一双输出的三阶系统为例,设系统的微积分方程为:
(35) 同单输入单输出系统一样,式(35)系统的实现也是非唯一的。现采用 模拟结构图的方法,按高阶导数项求解:
对每一个方程积分:
故得模拟结构图,如下图所示:
取每个积分器的输出为一个状态变量,如上图所示。则式(35)的一种 实现为:
1.1.3 状态方程 以状态变量
为坐标轴所构成的 维空间,称为
状态空间。
1.1.4 状态方程
由系统的状态变量构成的一阶微分方程组称为系统的状态方程。 用图下所示的 网络,说明如何用状态变量描述这一系统。
图一
根据电学原理,容易写出两个含有状态变量的一阶微分方程组:
亦即
(1)
式(1)就是图1.1系统的状态方程,式中若将状态变量用一般符号 ,
第一章 控制系统的状态空间表 达式
回忆:经典控制中学过的控制系统的数学表达式 微分方程、传递函数、方框图、信号流图、频率特性、伯德图等
a 0 ( t ) d d n c t ( n t ) a 1 ( t ) d d n t 1 n c ( 1 t ) a n 1 ( t ) d c d ( t t ) a n ( t ) c ( t ) b 0 ( t ) d d m t r m ( t ) b 1 ( t ) d d m t 1 m r ( 1 t ) b m 1 ( t ) d r d ( t t ) b m ( t ) r ( t )
现代控制理论知识点汇总
1.状态空间表达式n 阶DuCx y Bu Ax x+=+= 1:⨯r u 1:⨯m y n n A ⨯:r n B ⨯:n m C ⨯:rm D ⨯:A 称为系统矩阵,描述系统内部状态之间的联系;B为输入(或控制)矩阵,表示输入对每个状态变量的作用情况;C 输出矩阵,表示输出与每个状态变量间的组成关系,D直接传递矩阵,表示输入对输出的直接传递关系。
2.状态空间描述的特点①考虑了“输入-状态-输出”这一过程,它揭示了问题的本质,即输入引起了状态的变化,而状态决定了输出。
②状态方程和输出方程都是运动方程。
③状态变量个数等于系统包含的独立贮能元件的个数,n 阶系统有n 个状态变量可以选择。
④状态变量的选择不唯一。
⑤从便于控制系统的构成来说,把状态变量选为可测量或可观察的量更为合适。
⑥建立状态空间描述的步骤:a 选择状态变量;b 列写微分方程并化为状态变量的一阶微分方程组;c 将一阶微分方程组化为向量矩阵形式,即为状态空间描述。
⑦状态空间分析法是时域内的一种矩阵运算方法,特别适合于用计算机计算。
3.模拟结构图(积分器加法器比例器)已知状态空间描述,绘制模拟结构图的步骤:积分器的数目应等于状态变量数,将他们画在适当的位置,每个积分器的输出表示相应的某个状态变量,然后根据状态空间表达式画出相应的加法器和比例器,最后用箭头将这些元件连接起来。
4.状态空间表达式的建立1由系统框图建立状态空间表达式:a 将各个环节(放大、积分、惯性等)变成相应的模拟结构图;b 每个积分器的输出选作i x ,输入则为i x;c 由模拟图写出状态方程和输出方程。
2由系统的机理出发建立状态空间表达式:如电路系统。
通常选电容上的电压和电感上的电流作为状态变量。
利用KVL 和KCL 列微分方程,整理。
③由描述系统的输入输出动态方程式(微分方程)或传递函数,建立系统的状态空间表达式,即实现问题。
实现是非唯一的。
方法:微分方程→系统函数→模拟结构图→状态空间表达式。
《现代控制理论》第3版课后习题答案
《现代控制理论参考答案》第一章答案1-1 试求图1-27系统的模拟结构图,并建立其状态空间表达式。
图1-27系统方块结构图解:系统的模拟结构图如下:图1-30双输入--双输出系统模拟结构图系统的状态方程如下:u K K x K K x K K x X K x K x x x x J K x J x J K x J K x x J K x x x pp p p n pb1611166131534615141313322211+--=+-==++--===••••••令y s =)(θ,则1x y =所以,系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡••••••654321165432111111112654321000001000000000000010010000000000010x x x x x x y uK K x x x x x x K K K K K K J K J J K J K J K x x x x x x p p pp n p b1-2有电路如图1-28所示。
以电压)(t u 为输入量,求以电感中的电流和电容上的电压作为状态变量的状态方程,和以电阻2R 上的电压作为输出量的输出方程。
L1L2U图1-28 电路图解:由图,令32211,,x u x i x i c ===,输出量22x R y =有电路原理可知:•••+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=•••写成矢量矩阵形式为:[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡32121321222111321000*********x x x R y u L x x x CCL L R L L R x x x 。
《现代控制理论》刘豹著(第3版)课后习题答案(最完整版)
第一章习题答案1-1 试求图1-27系统的模拟结构图,并建立其状态空间表达式。
11K s K K p +sK s K p 1+s J 11sK n 22s J K b -++-+-)(s θ)(s U 图1-27系统方块结构图解:系统的模拟结构图如下:)(s U )(s θ---+++图1-30双输入--双输出系统模拟结构图1K pK K 1pK K 1+++pK n K ⎰⎰⎰11J ⎰2J K b ⎰⎰-1x 2x 3x 4x 5x 6x系统的状态方程如下:uKK x KK x KK x X K x K x x x x J Kx J x J K x J Kx x J K x x x ppppn pb 1611166131534615141313322211+--=+-==++--===∙∙∙∙∙∙令y s =)(θ,则1x y =所以,系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∙∙∙∙∙∙65432116543211111111265432100000100000000000000010010000000000010x x x x x x y uK K x x x xx x K K K K K K J K J J K J KJ K x x x x x xp p pp n pb1-2有电路如图1-28所示。
以电压)(t u 为输入量,求以电感中的电流和电容上的电压作为状态变量的状态方程,和以电阻2R 上的电压作为输出量的输出方程。
R1L1R2L2CU---------Uc---------i1i2图1-28 电路图解:由图,令32211,,x u x i x i c ===,输出量22x R y =有电路原理可知:∙∙∙+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x Cx Cx x L x L R x uL x L x L R x =+-=+-=+--=∙∙∙写成矢量矩阵形式为:[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡32121321222111321000010111010x x x R y u L x x x CCL L R L L R x x x 。
现代控制理论(1-8讲第1-2章知识点)精品PPT课件
dia dt
Ke
I fD Coபைடு நூலகம்st
n f Const
nDJ , f
其中:Kf 为发电机增益常数;Ke 为电动机反电势常数。
(3).电动机力矩平衡方程:J
d
dt
f
Kmia
(Km
-电动机转矩常数)
以上三式可改写为:
d
dt
f J
Km J
ia
dia dt
Ke Ra
La
La
ia
Kf La
if
试写出其状态空间表达式。
解:选择相变量为系统的状态变量,有
•
•
•• •
x1 y x2 y x1 x3 y x2
故
即
•
x1 x2
•
x2 x3
•
x3
a0 a3
x1
a1 a3
x2
a2 a3
x3
1 a3
u
•
0
x 0
a0
a3
1 0 a1 a3
0
0
1 x 0 u
a2
1
a3 a3
a1 y a0 y
bnu (n)
b u (n1) n 1
b0u
(1)
分为两种情况讨论。
一、输入信号不含有导数项:
此时系统的运动方程为:
•
y(n)
a y(n1) n1
a1 y a0 y b u
故选
x1 y
•
x2 y
..
xn1
y(n2)
xn y(n1)
对左边各式求导一次,即有
18
24
2-3 化系统的频域描述为状态空间描述
现代控制理论--2控制系统的状态空间模型
4. 非线性定常系统:
X (t ) f X (t ) u(t ) Y (t ) g X (t ) u(t )
5.非线性时变系统:
x(t) f x(t), u(t), t y(t) g x(t), u(t), t
6.线性系统状态空间表达式的简便写法:
x1
y
x2 y
令
xn
1
yn2
xn
y n 1
x1 x2
x2
x3
xn1
xn
xn
yn
a1 yn1
=-an x1 an1x2
an1 y a1xn bu
an y
bu
0
K F(t)
y(t)
f
弹簧-质量-阻尼器系统
解:列基本方程:
d2y
dy
m dt 2 f dt ky u t
选择状态变量:取:
x1 (t) y(t)
故得:
x2(t) y(t)
x1(t) x2 (t)
x2 (t)
k m
x1
f m
x2
1 m
u
y(t) x1
将以上方程组写矩阵形式
u(t) 1
La
Ra
+++
x1 L
dt
x1
Ca
L
x2
dt
x2 1 Y(t)
1
Cm
J
+ x3 +
现代控制理论章节习题含答案(大学期末复习资料)
《现代控制理论》第一章习题解答1.1线性定常系统和线性时变系统的区别何在?答:线性系统的状态空间模型为:x = AxBu+y CxDu= +线性定常系统和线性时变系统的区别在于:对于线性定常系统,上述状态空间模型中的系数矩阵A,B,C和中的各分量均为常数,而对线性时变系统,其系数矩阵D A,B,C和D中有时变的元素。
线性定常系统在物理上代表结构和参数都不随时间变化的一类系统,而线性时变系统的参数则随时间的变化而变化。
1.2 现代控制理论中的状态空间模型与经典控制理论中的传递函数有什么区别?答: 传递函数模型与状态空间模型的主要区别如下:1.3 线性系统的状态空间模型有哪几种标准形式?它们分别具有什么特点?答: 线性系统的状态空间模型标准形式有能控标准型、能观标准型和对角线标准型。
对于n 阶传递函数G s( )= b s n−s1nn+−1a s+n−b s1n−n2−1n+−2 + +as a+1 bs b+1 +0 0+d ,分别有⎧⎡0 1 0 0 ⎤⎡⎤0⎪⎢0 0 1 0 ⎥⎥⎢⎥⎢⎥0⎪⎢⎪⎪x =⎢ ⎥x+⎢⎥ u ⑴能控标准型:⎨⎢0 0 0 1 ⎥⎥⎢⎥⎢⎥0⎪⎢⎪⎣⎢−a0 −a1 −a2 −a n−1⎥⎦⎢⎥⎣⎦1⎪⎪⎩y=[b0 b1 b n−2 b n−1]x du+⎧⎡0 0 0 −a0 ⎤⎡b0 ⎤⎪⎪⎢⎢1 0 0 −a1 ⎥⎥⎢⎢b1 ⎥⎥⎪⎪x =⎢0 1 0 −a2 ⎥⎥x+⎢⎢ ⎥⎥u⑵能观标准型:⎨⎢b n−2⎥⎪⎢ ⎥⎢⎪⎣⎢0 0 1 −a n−1⎦⎥⎢⎣b n−1⎥⎦⎪⎪⎩y=[0 0 0 1]x du+⎧⎡p1⎪⎢0⎪x =⎢⎢ 0 p20 0 ⎤⎡1⎤0 ⎥⎢1⎥⎥x+⎢⎥u ⎥⎢ ⎥⎪⑶对角线标准型:⎨⎪⎢⎣0⎪p n⎥⎦⎢⎣1⎥⎦⎪⎩y=[c1 c2 c x du n] + 式中的pp1, 2,, p n和c c1, 2,, c n可由下式给出,G s( )= b s n−s1nn−1a s+n−b s1n−n2−1n+−2 + +as a+1 bs b+1 +0 0 + =d s p−c1 1 + s p−c2 2 + + s p−c n n +d+能控标准型的特点:状态矩阵的最后一行由传递函数的分母多项式系数确定,其余部分具有特定结构,输出矩阵依赖于分子多项式系数,输入矩阵中的元素除了最后一个元素是1 外,其余全为0。
现代控制理论第一章-控制系统数学模型
y b0
b1
bn1
xn
注:如果输入项的导数阶次和输出项导数阶次相同,则有d。
Y (s) R(s)
bn s n an s n
b1s b0 a1s a0
d
bn1sn1 b1s b0 ansn a1s a0
例1-4 已知描述系统的微分方程为 y18y 192y 640y 160u 640u
y bn1z(n1) b1z b0 z b0 x1 b1x2 bn1xn
写成矩阵形式
x1
x2
xn
0
0
0
a0
1 0 0 a1
0 1 0 a2
0 0 0 a3
0
0
0 1 an1
x1 x2
xn
0 u 0
1
x1
第1章 控制系统数学模型
本课程的任务是系统分析和系统设计。而不论是系统分析还是系 统设计,本课程所研究的内容是基于系统的数学模型来进行的。因 此,本章首先介绍控制系统的数学模型。
本章内容为: 1、状态空间表达式 2、由微分方程求出系统状态空间表达式 3、传递函数矩阵 4、离散系统的数学模型 5、线性变换(状态变量选取非唯一)
写成矩阵形式
x1 0 1 0 x1 0
x2
0
0
1
x2
0
u
x3 a0 a1 a2 x3 b0
x1
y 1
0
0
x2
x3
状态图如下:
一般情况下,n 阶微分方程为: y(n) an1 y(n1) a1 y a0 y b0u
选择状态变量如下:
x1 y x1 x2 y x2 x3 y
0
x2
1 M
现代控制理论_第1章
现代控制理论 机械工程硕士研究生学位课
状态空间表达式
ɺ = Ax + bu x y = cx
x1 0 x 0 2 x = ⋮ , A = ⋮ xn −1 0 xn −a0 1 0 ⋮ 0 −a1 0 1 ⋮ ⋯ ⋯ ⋱ 0 0 0 0 ⋮ , b = ⋮ , c = [1 0 ⋯ 0] 1 0 −an −1 b0
现代控制理论 机械工程硕士研究生学位课
实现问题
实现问题:由描述系统输入-输出动态关系 的运动方程式或传递函数,建立系统的状态 空间表达式。 揭示系统的内部关系 讨论单输入单输出线性定常系统
ɺ = Ax + bu x y = cx + du
现代控制理论 机械工程硕士研究生学位课
两类实现问题
di 1 Ri + L + ∫ idt = u dt C
本例子中 1. 输入和输出都已 明确; 2. 选择两个独立的 储能元件作为状 态变量; 3. 根据电路的基本 定律列出方程
现代控制理论 机械工程硕士研究生学位课
1 y = uc = ∫ idt C
系统状态方程的建立
设状态变量为电感器电流和电容器电压,即
现代控制理论 Modern Control Theory
现代控制理论 机械工程硕士研究生学位课
本课程主要内容
系统描述:状态空间表示法 系统分析:状态方程的解、线性系统的能 控和能观测性、稳定性分析 系统设计:状态反馈和状态观测器 最优控制:最优控制系统及其解法
现代控制理论第一章(吴忠强版)
吴忠强
目
录
第一章 控制系统的状态空间表达式 第二章 控制系统状态空间表达式的解 第三章 线性控制系统的能控性与能观性 第四章 控制系统的李亚普诺夫稳定性 第五章 线性定常系统的综合 第六章 最优控制系统设计 参考文献
内容简介
•
本书系统的介绍了现代控制理论的 基本内容,包括控制系统的状态空间描 述、运动分析与离散化、李亚普诺夫稳 定性分析、能控性与能观性、状态反馈 与状态观测器、最优控制系统设计。每 章配有一定的例题和习题.
b11 b 21 B bn1
b12 b 22 bn 2
b1 r b2 r b nr
y1 y2 y ym
——m维输出矢量;
—— n r 输入(或控制)矩阵;
c 11 c 12 c 21 c 22 C c m1 c m 2
1
式(1-3)就是图1-1系统的输出方程,它的矩阵表示为
y 1
T
0
x1 x2
或
y C x
T
y c x
T
(1-4)
式中
c
1
0
六、状态空间表达式
l 状态方程和输出方程总合起来,构成对一个系统完整的动态 描述,称为系统的状态空间表达式, 在经典控制理论中,用指定某个输出量的高阶微分方程来描 述系统的动态过程。如图1-1所示的系统,在以 uc 作输出时, 从式(1-1)消去中间变量i ,得到二阶微分方程为
回到式(1-5)或式(1-6)的二阶系统,若改选 u C 和 u c 作为 两个状态变量,即令 x 1 u C ,
x2 uc
王金城现代控制理论第一章_习题答案
王金城化工出版社第1章习题参考答案:1-1(a )选123123,,,,,y y y v v v 为状态变量,根据牛顿定律,对1M ,有()11112121dv M g K y K y y M dt---= 对2M ,有()()222123232dv M g K y y K y y M dt+---= 对3M ,有()33323433dv M g K y y K y M dt+--= 令312112233415263,,,,,dy dy dyx y x y x y x v x v x v dt dt dt=========,整理得 ()()()122214253641112334233251262322233,,,,,K K K x x x x x x x x xg M M K K K K K x K K xx x g x x x g M M M M M +====-++++=-++=-+()()()1221123222223433300010000001000000010000001100010000K K K M M x x g K K K K M M M K K K M M ⎡⎤⎢⎥⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥+⎢⎥-⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥+⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥+-⎢⎥⎢⎥⎣⎦100000010000001000y x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦(b )选12,12,,y y v v 为状态变量,根据牛顿定律,对1M ,有()11121111dv M g B v v K y M dt+--= 对2M ,有()22221212dv f M g B v B v v M dt+---= 令1211223142,,,dy dyx y x y x v x v dt dt ======,整理得 11113243134111,,K B Bxx x x x x x x g M M M ===--++ ,112434222B B B f x x x g M M M +=-++所以状态空间描述为1111111122220010000001000011100K B B xx g f M M M B B B M M M ⎡⎤⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥=++⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥+⎣⎦⎢⎥-⎢⎥⎣⎦⎣⎦10000100y x ⎡⎤=⎢⎥⎣⎦1-2(a )取电感电流i 和电容电压u 为状态变量,列回路方程122c rc c c u u R (i )u u R di L u u dt u du C dt R ⎧=+++⎪⎪⎪=+⎨⎪⎪=⎪⎩令12c x i,x u,y u ===()1212121212112121211r R R R R L(R R )L(R R )L(R R )xx u R C R RC(R R )C(R R )-⎡⎤⎡⎤⎢⎥⎢⎥+++⎢⎥⎢⎥=+⎢⎥⎢⎥--⎢⎥⎢⎥+++⎣⎦⎣⎦1222121212r R RR R y x u R R R R R R ⎡⎤=--+⎢⎥+++⎣⎦ (b )选择回路电流a i 和电枢角速度ω为状态变量,有aa a a ae di u R i L K dt ω=++ 力矩平衡方程:a a d J B K i ,dtωω+= 其中a K 为转矩常数 1a a e a a a a adi R K i u dt L L L ω=--+a a K d B i dt J J ωω=-- 令12a x i ,x ,ω==有10a e a aa a R K L L L xx u K B JJ -⎡⎤-⎡⎤⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥-⎣⎦⎢⎥⎣⎦ , []01y x ω==1-3 (1)传递函数为3221375Y(s )U(s )s s s =+++将传递函数中的公因子提出,于是有3123211375Y(s )s U(s )s s s----=+++ 按梅逊公式构建系统的状态变量图能控标准形:0100001057131x x ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦u []200y x =能观标准形:0052107001130x x u -⎡⎤⎡⎤⎢⎥⎢⎥=-+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦[]001y =x(2)传递函数为:2332132223123Y(s )s s s U(s )s s s s----++==++++ 按梅逊公式构建系统的状态变量图能控标准形:010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦[]210y x =能观标准形:003210010120x x u -⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦[]001y x =(3)传递函数为:3212332123324515471547Y(s )s s s s s s U(s )s s s s s s------+++---==+++++++ 按梅逊公式构建系统的状态变量图状态空间描述为:010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦[]514y x u =---+(4)①12121221212121b s b b s b s Y(s )U(s )s a s a a s a s ----++==++++ 状态空间描述为:1322140101xx u x a a x ⎡⎤⎡⎤⎡⎤⎡⎤=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦⎣⎦ ,[]21y b b x =②22121201200111s c Z(s )c c Y(s )c s c s c s s c c ---==++++ 状态空间描述为:332144000101x x y c c x x c c ⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦ ,301z x c = 两系统串联,得112122332121440001000001000100x x a a x x u x x c c b b x x c c ⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦(5)由G(z)有,y(k+3)+4y(k+2)+5y(k+1)+2y(k)=u(k)令12312x (k )y(k )x (k )y(k )x (k )y(k )=⎧⎪=+⎨⎪=+⎩ 1230100100102541x (k )x(k )x (k )u(k )x (k )⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥+=+⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦⎣⎦[]123100x (k )y(k )x (k )x (k )⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦(6)由G(z)有,y(k+3)+6y(k+2)+11y(k+1)+6y(k)=2u(k+2)+u(k+1)+2u(k)01001001061161x(k )x(k )u(k )⎡⎤⎡⎤⎢⎥⎢⎥+=+⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦[]123212x (k )y(k )x (k )x (k )⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦1-4 (a )化简系统结构图得系统状态空间描述:1234010000010024220025025x x x u x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦ []0100y x =(b) 化简系统结构图得系统状态空间描述:1112221323255223735353xx u ///x x u ////--⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤=+⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦ []110y x = []201y x =1-5 (1) 传递函数为21233212332322461246s s s s s G(s )s s s s s s ------++++==++++++ 能控标准形:010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦[]231y x =能观标准形:006210430121x x u -⎡⎤⎡⎤⎢⎥⎢⎥=-+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦[]001y x =(2)传递函数为24422431332132s s s G(s )s s s s -----+-+==++++ 能控标准形:01000001000001020301xx u ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦[]1300y x =-能观标准形:00021100030103000100x x u -⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥=+⎢⎥⎢⎥-⎢⎥⎢⎥⎣⎦⎣⎦[]0010y x = 1-6(1) 24512122123123(s )(s )G(s )(s )(s )(s )s s s ++-==++++++++状态空间描述:100102010031x x u -⎡⎤⎡⎤⎢⎥⎢⎥=-+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]12122y x =- (2)223533313313(s )G(s )(s )(s )s s (s )+--==+++++++ 状态空间描述:310003010011x x u -⎡⎤⎡⎤⎢⎥⎢⎥=-+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]333y x =-- 1-7(1)∵31I A ()()λλλ-=++ 1213,λλ=-=-∴1003A -⎡⎤=⎢⎥-⎣⎦11111111,p Ap ,p λλ⎡⎤=-==⎢⎥⎣⎦,22222131,p Ap ,p λλ⎡⎤=-==⎢⎥-⎣⎦∴1111P ⎡⎤=⎢⎥-⎣⎦1111112P ---⎡⎤=-⎢⎥-⎣⎦ 11112B P B -⎡⎤==⎢⎥-⎣⎦∴11020312x x u ⎡⎤⎢⎥-⎡⎤=+⎢⎥⎢⎥-⎣⎦⎢⎥-⎢⎥⎣⎦(2)1230123I A ,,,λλλλ-==-=-=-∴100020003A -⎡⎤⎢⎥=-⎢⎥⎢⎥-⎣⎦111111111,p Ap ,p λλ⎡⎤⎢⎥=-==-⎢⎥⎢⎥-⎣⎦,2222212212,p Ap ,p λλ⎡⎤⎢⎥⎢⎥=-==-⎢⎥⎢⎥⎢⎥⎣⎦333331333,p Ap ,p λλ⎡⎤⎢⎥=-==-⎢⎥⎢⎥⎣⎦∴1111231132P ⎡⎤⎢⎥⎢⎥=---⎢⎥⎢⎥-⎢⎥⎣⎦137272304027162B P B -⎡⎤⎢⎥⎢⎥==--⎢⎥⎢⎥⎢⎥⎣⎦∴37271002020304000327162x x u ⎡⎤⎢⎥-⎡⎤⎢⎥⎢⎥=-+--⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎢⎥⎣⎦(3)519400433030114003433114j x j x j j ⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥-+⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦+⎢⎥⎢⎥⎣⎦1-8 (1)∵A 为友矩阵123012I A ,,λλλλ-====∴ 110010002A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ 101112124P ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ 1111B P B --⎡⎤⎢⎥==-⎢⎥⎢⎥⎣⎦ ∴100101010021x x u -⎡⎤⎡⎤⎢⎥⎢⎥=+-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(2)212331031I A ()(),,λλλλλλ-=--====310030001A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ 120112111P ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ 11335234B P B --⎡⎤⎢⎥==-⎢⎥⎢⎥-⎣⎦ ∴3101330305200134x x u -⎡⎤⎡⎤⎢⎥⎢⎥=+-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦1-9(1)110061031002P -⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦111000062300100111000152020233302100313000222AP AP -⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦1106203502BP B -⎡⎤⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦203640C CP ⎡⎤==⎢⎥⎣⎦ 11000621102203333502022xx u ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦203640y x ⎡⎤=⎢⎥⎣⎦ (2)①1111I A I P AP P P P AP P (I A)P λλλλ-----=-=-=-11I A P I A P P P I A I A λλλλ---=-=-=-∴特征值不变②1111G(s )C(sI A)B CP(sI P AP )P B ----=-=- 111C P(sI P AP )P B ---⎡⎤=-⎣⎦11111C P(sI )P PP APP B C(sI A )B -----⎡⎤=-=-⎣⎦∴传递函数不变1-10证明:11G (s )c(sI A)b -=- 12G (s )c(sI A)b -=- ∵T T TA A ,b c ,c b ===∴12T T T T T T TG (s )b (sI A )c b (sI )A c -⎡⎤=-=-⎣⎦ [11T T T T T Tb (sI A )c b (sI A)c --⎤⎡⎤=-=-⎦⎣⎦11TTc(sI A )b G (s )-⎡⎤=-=⎣⎦ ∵系统为单输入单输出,11T G (s )G (S )= ∴两者传递函数相同。
(完整版)现代控制理论
(完整版)现代控制理论第⼀章线性离散系统第⼀节概述随着微电⼦技术,计算机技术和⽹络技术的发展,采样系统和数字控制系统得到⼴泛的应⽤。
通常把采样系统,数字控制系统统称为离散系统。
⼀、举例⾃动测温,控温系统图;加热⽓体图解:1. 当炉温h变化时,测温电阻R变化→R,电桥失去平衡状态,检流计指针发⽣偏转,其偏转⾓度为)e;(t2. 检流计是个⾼灵敏度的元件,为防磨损不允许有摩擦⼒。
当凸轮转动使指针),接触时间为τ秒;与电位器相接触(凸轮每转的时间为T样偏差)(t e 是连续信号,电位器的输出的e *τ(t)是脉冲信号。
连续信号转变为脉冲信号的过程,成为采样或采样过程。
实现采样的装置成为采样器。
To —采样周期,f s =--To1采样频率,W s =2πf s —采样⾓频率 2.信号复现因接触时间很⼩,τo T ??τ,故可把采样器的输出信号)(t e *近似看成是⼀串强度等于矩形脉冲⾯积的理想脉冲,为了去除采样本⾝带来的⾼额分量,需要把离散信号)(t e *恢复到原信号)(t e 。
实现⽅法:是在采样器之后串联⼀个保持器,及信号复现滤波器。
作⽤:是把)(t e *脉冲信号变成阶梯信号e h (t)3.采样系统结构图r(t),e(t),c(t),y(t)为连续信号,)(t e *为离散信号)(s G h ,)(s G p ,)(s H 分别为保持器,被控对象和反馈环节的传递函数。
(t)r4.采样系统⼯作过程由保持器5. 采样控制⽅式采样周期To ??=≠=?相位不同步采样常数常数6. 采样系统的研究⽅法(或称使⽤的数字⼯具)因运算过程中出现s 的超越函数,故不⽤拉式变换法,⼆采⽤z 变换⽅法,状态空间法。
第⼆节信号的采样和复现第⼀节是定性认识与分析,本节是定量研究。
⼀、采样过程从第3个图形可知,采样器输出信号)(t e *是⼀串理想的脉冲信号,k 瞬时)(t e *的脉冲强度等于此时)(T e 的幅值)(0kT e ,即)0(0T e ,)(0T e ,)2(0T e …. )(0nT e ….采样过程可以看成为⼀个幅值调制过程,采样器如同⼀个幅值调制器。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
其对应的状态空间表达式分别为
1 A x 1 x1 B 1u1 y1 C1 x1 D1u1
2 A2 x 2 B2 u2 x y2 C2 x2 D2 u2
从图可知 u1=u2=u y1+y2=y 故可导出并联联结组合系统的 状态空间模型为
1 A1 x x 2 0
[I+G0(s)F(s)]Y(s)=G0(s)U(s) Y(s)=[I+G0(s)F(s)]-1G0(s)U(s) 反馈联结组合系统的传递函数为 G(s)=[I+G0(s)F(s)]-1G0(s) 或 G(s)=G0(s)[I+F(s)G0(s)]-1
1.6 状态向量的线性变换和状态空间表达式的 特征标准型
因此令x(0)=0,于是由状态方程的拉氏变换式有
X(s)=(sI-A)-1BU(s)
将上述X(s)代入输出方程,有
Y(s)=[C(sI-A)-1B+D]U(s)
线性定常连续系统的传递函数阵为
G(s)=C(sI-A)-1B+D
若对于输入与输出间无直接关联项(即D=0)的系统,则有
G(s)=C(sI-A)-1B
1.4 从状态空间表达式求系统传递函数(阵)
前面已经介绍了SISO系统从传递函数求系统的状态空间表达 式,下面将介绍其逆问题, • 即怎样从状态空间表达式求系统的传递函数阵。
已知MIMO线性定常系统的状态空间表达式为
x Ax Bu y Cx Du
其中x为n维状态向量; u为r维输入向量; y为m维输出向量。
1 1 1 C ( sI A ) B D C ( s I A ) B1 D1 2 2 2 2 1 1
G2 ( s )G1 ( s )
串联联结组合系统的传递函数阵为串联系统各子 系统的传递函数阵的顺序乘积。
应当注意,由于矩阵不满足乘法交换律,故在上式中 G1(s)和G2(s)的位置不能颠倒,它们的顺序与它们在 系统中的串联联结顺序一致。
x2 A2 x2 B2 u2 y2 C2 x2
• 从图可知
u1=u-y2 u2=y1=y 因此可导出反馈组合系统的状 态空间模型为
x1 A1 x1 B1u1 A1 x1 B1 (u y2 ) A1 x1 B1C2 x2 B1u x2 A2 x2 B2 u2 A2 x2 B2 y1 A2 x2 B2C1 x1 y y1 C1 x1
x Ax Bu 对上式取拉氏变换,有 y Cx Du sX ( s ) x(0) AX ( s ) BU ( s ) Y ( s ) CX ( s ) DU ( s )
其中X(s)、U(s)和Y(s)分别为x(t)、u(t)和y(t)的拉氏变换; x(0)为x(t)的在初始时刻t=0的值。 由于传递函数阵描述的是系统输入输出间动态传递关系,不考 虑系统初始条件的影响。
x1 y D2C1 C2 ( D2 D1 )u x2
串联连接组合系统的状态变量的维数为子系统的 状态变量的维数之和。 由串联组合系统的状态空间模型可求得组合系统 的传递函数阵为
A G ( s ) D2C1 C2 sI 1 B2C1 D2 C1 0 B1 D2 D1 A2 B2 D1 B1 D2 D1 1 BD sI A2 2 1 0
即有
x1 A1 B1C2 x1 B1 u A2 x2 0 x2 B2 C1 y C 0 x1 1 x 2
反馈联结组合系统的状态变量的维数为子系统的状态变量 的维数之和。 • 由反馈联结组合系统的联结图可知 Y(s)=G0(s)U1(s)=G0(s)[U(s)-Y2(s)]=G0(s)[U(s)-F(s)Y(s)]
1 1
1 sI A 1 C2 1 1 sI A B C sI A 2 2 1 1 1 1
D2 C1 sI A1 B1 C2 sI A2 B2 C1 sI A1 B1 C2 sI A2 B2 D1 D2 D1
1
C1 ( sI A1 ) 1 B1 D1 C2 ( sI A2 ) 1 B2 D2 G1 ( s ) G2 ( s )
并联组合系统的传递函数阵为各并联子系统的传递函数 阵之和。
2) 串联连接
串联联接组合系统方块结构图
设图所示的串联联结的组合系统的两个子系统的传 递函数阵分别和并联连结的结构相同,其对应的状态 空间表达式也分别相同。
对r维输入、m维输出的MIMO系统,若其输入输出的拉氏变换 分别为U(s)和Y(s),则系统的输入输出间的动态关系可表示为 Y(s)=G(s)U(s) 其中G(s)称为传递函数阵,其每个元素为标量传递函数。 • G(s)的形式为 G11 ( s ) G12 ( s ) ... G1r ( s ) G ( s ) G ( s ) ... G ( s ) 22 2r G ( s ) 21 ... ... ... ... Gm1 ( s ) Gm 2 ( s ) ... Gmr ( s ) 其中Gij(s)描述了第i个输出与第j个输入之间的动态传递关系。
• 相应的输出方程为
y y 2 C2 x 2 D2 u2 C2 x 2 D2 (C1 x1 D1u1 ) D2C1 x1 C2 x 2 D2 D1u
1 A1 x x 2 B2C1
0 x1 B1 u A2 x2 B2 D1
0 0 X k k 1 4 T1 k3 T3 1 T2 0
y x1
写成矩阵形式:
y 1 0 0X
x 0 0 1 k2 x2 0 u T2 1 x k1 T1 3 T1
由组合系统的状态空间表达式可求得组合系统的传递函 数阵为
G ( s ) C1 C1 A1 C2 sI 0 sI A1 1 C2 0 0 B1 ( D1 D2 ) A2 B2 B1 0 ( D1 D2 ) 1 B sI A2 2
[例 ] 0 求由 A 0
[解 ]:
0 1 0 1 1 0 所表述系统的W(s) 0 1 ,B 2 1 ,C 2 1 1 6 11 6 0 2 1
由传递函数矩阵公式得:
W( s) C ( sI A) 1 B D s 1 0 1 1 0 0 s 1 2 1 1 6 11 s 6
3) 反馈连接
反馈连接组合系统结构图
• 设对应于图所示的反馈联结 组合系统的两个子系统的传 递函数阵为
G0 (s) C1 (sI A1 )1 B1
F (s) C2 (sI A2 )1 B2
其对应的状态空间模型分别为
x1 A1 x1 B1u1 y1 C1 x1
从图可知 u1=u u2=y1 y2=y 因此可导出串联组合系统的状态空间方程为
1 A1 x1 B1u1 A1 x1 B1u x
2 A2 x 2 B2 u2 A2 x 2 B2 y1 x A2 x 2 B2 (C1 x1 D1u1 ) B2C1 x1 A2 x 2 B2 D1u
例1(见书)
u
k1 T s 1 1
k2 T2 s 1
k3 T3 s
y
k4
3 x
u
k1 T1
1 T1
x3 k 2
T2
2 xΒιβλιοθήκη 1 T2x2
k3 T3
1 x
x1
y
k4
图中有三个积分环节,三阶系统,取三个状态变量如上图(选 择积分环节后的变量为状态变量): 则有:
k3 1 x2 x T1 1 k2 x2 x2 x3 T2 T2 k1 1 k1 3 k4 x1 x3 u x T1 T1 T1
C1
0 x1 B1 u A2 x 2 B2
y C1 x1 D1u1 C 2 x 2 D2 u2 x1 C 2 ( D1 D2 )u x2
因此,由上述状态空间表达式可知,并联组合系统的状态 变量的维数为子系统的状态变量的维数之和。
Ax Bu x y Cx Du
G( s) C ( sI A)1 B D
SISO系统,用传递函数G(s)描述,G(s)是一个元素;
MIMO系统,多个输入对多个输出,故引入传递函数矩阵 G(s) ,G(s)是一个矩阵,可以表征多个输入对系统输出的影 响; 同一系统,不同的状态空间表达式对应的传递函数阵应是相 同的。即描述系统输入与输出间动态传递关系的传递函数阵对 状态变换具有不变性。
1
求得传递函数阵为:
2 2 s 4s 29 s 3s 4 1 W( s) 3 2 2 2 s 6s 11s 6 4s 56s 52 3s 17s 14
1.5 组合系统的状态空间描述和传递函数矩阵
已知两独立子系统的状态空间描述和传递函数如下
状态空间模型不具有唯一性.
第一章
控制系统的状态空间描述
4、根据系统传递函数的方块图建立状态空间表达式
典型环节转换成结构图