线性系统状态空间分析与运动解

控制策略设计
基于状态空间模型，可以设计各种控制策略，如 线性二次调节器（LQR）、最优控制等，以实现 系统的最优控制。
控制系统仿真
通过状态空间模型，可以对控制系统进行仿真， 模拟系统的动态行为，评估控制策略的有效性和 性能。
在信号处理中的应用
信号滤波
状态空间方法可以用于信号滤波， 通过构建状态空间模型来描述信 号的动态变化，实现信号的平滑 和噪声抑制。
状态空间的分解
contents
目录
• 状态空间的基本概念 • 线性系统的状态空间表示 • 状态空间的分解方法 • 状态空间的应用 • 状态空间分解的实例分析
01
状态空间的基本概念
状态变量的定义
状态变量
01
描述系统状态的变量，通常用矢量表示，包含系统的各个独立
变量。
状态变量的选择
02
选择的状态变量应能全面反映系统的动态特性，且便于分析。
线性系统的状态方程
状态方程描述了系统内部状态变 量随时间的变化规律，通常表示 为状态变量的一阶或二阶微分方
程。
对于线性系统，状态方程具有形 式：dx/dt = Ax + Bu，其中x 是状态变量，u是输入，A和B是
系统矩阵。
解状态方程可以得到系统状态变 量的时间响应。
线性系统的输出方程
输出方程描述了系统输出与状态变量和输入之间的关系，通常表示为输出变量与状 态变量的线性组合。
总结词
高阶线性系统的状态空间表示能够精细 地描述系统的动态行为。
VS
详细描述
高阶线性系统是指系统的动态行为需要用 高阶微分方程来描述的系统。其状态空间 表示与一阶和二阶系统类似，但需要更多 的状态变量和方程来描述系统的动态行为 。通过高阶线性系统的状态空间表示，可 以更精确地分析系统的动态性能和稳定性 ，以及设计更有效的控制系统。

称 G ( s ) 为系统的传递函数矩阵。G ( s ) 为的一个有理分式 矩阵。当 g ij ( s ) 除严格真还包含真有理分式时，即 G ( s ) 的一个或一些元传递函数中分母和分子多项式具有相等 的最高幂次时，称为真有理分式矩阵。
7
§1.1-2 传递函数矩阵 当且仅当 G ( s )为真的或严格真的时，它才是物理上可实 现的。当且仅当 lim G ( s ) = 零阵 s →∞ G ( s ) 为严格真的， lim G ( s ) =非零常阵 s →∞ 传递函数矩阵为真的。
8
§1.2 线性定常系统的状态空间描述
§1.2-1 状态和状态空间 系统的状态空间描述是建立在状态和状态空间概念的基 础上的。 定义1.1 动力学系统的状态定义为完全的表征系统时间 域行为的一个最小内部变量组。组成这个变量组的变 xn (t ) 称为系统的状态变量，其中t ≥ t0， "， 量 x1 (t ), x2 (t )， t0 为初始时刻。由状态变量 ⎡ x1 (t ) ⎤ ⎥, t ≥ t 构成的列向量 x(t ) = ⎢ # 0 ⎢ ⎥ 称为系统的状态向量，简称为状态。状态空间则定义为 状态向量取值的一个向量空间。
15
§1.2-2 动态系统的状态空间描述 离散动态过程的状态空间的描述。离散动态过程的一个 重要特点是，系统的各个变量都被处理成为只在离散时 刻取值，其状态空间描述只反映离散时刻的变量组间的 因果关系和转换关系。用k=0,1,2来表示离散的时刻，则 离散时间系统（简称离散系统）的状态方程和输出方程 的最一般形式为：
2
§1.1-1 单变量情形回顾 已知由下列常系数微分方程描述的定常系统
y n + a n −1 y ( n −1) + " + a1 y (1) + a 0 y

第九章 线性系统的状态空间分析与综合在第一章至第七章中，我们曾详细讲解了经典线性系统理论以及用其设计控制系统的方法。

可以看到，经典线性理论的数学基础是拉普拉斯变换和z 变换，系统的基本数学模型是线性定常高阶微分方程、线性常系数差分方程、传递函数和脉冲传递函数，主要的分析和综合方法是时域法、根轨迹法和频域法，分析的主要内容是系统运动的稳定性。

经典线性系统理论对于单输入-单输出线性定常系统的分析和综合是比较有效的，但其显著的缺点是只能揭示输入-输出间的外部特性，难以揭示系统内部的结构特性，也难以有效处理多输入-多输出系统。

在50年代蓬勃兴起的航天技术的推动下，在1960年前后开始了从经典控制理论到现代控制理论的过渡，其中一个重要标志就是卡尔曼系统地将状态空间概念引入到控制理论中来。

现代控制理论正是在引入状态和状态空间概念的基础上发展起来的。

在现代控制理论的发展中，线性系统理论首先得到研究和发展，已形成较为完整成熟的理论。

现代控制理论中的许多分支，如最优控制、最优估计与滤波、系统辨识、随机控制、自适应控制等，均以线性系统理论为基础；非线性系统理论、大系统理论等，也都不同程度地受到了线性系统理论的概念、方法和结果的影响和推动。

现代控制理论中的线性系统理论运用状态空间法描述输入-状态-输出诸变量间的因果关系，不但反映了系统的输入—输出外部特性，而且揭示了系统内部的结构特性，是一种既适用于单输入--单输出系统又适用于多输入—多输出系统，既可用于线性定常系统又可用于线性时变系统的有效分析和综合方法。

在线性系统理论中，根据所采用的数学工具及系统描述方法的不同，又出现了一些平行的分支，目前主要有线性系统的状态空间法、线性系统的几何理论、线性系统的代数理论、线性系统的多变量频域方法等。

由于状态空间法是线性系统理论中最重要和影响最广的分支，加之受篇幅限制，所以本章只介绍线性系统的状态空间法。

9-1 线性系统的状态空间描述1. 系统数学描述的两种基本类型这里所谓的系统是指由一些相互制约的部分构成的整体，它可能是一个由反馈闭合的整体，也可能是某一控制装置或受控对象。

线性系统状态空间分析和运动解状态空间分析方法是一种用来描述线性系统的分析方法。

它将系统的动态特性用一组状态变量来表示，并通过矩阵形式的状态方程进行分析和求解。

状态空间方法是目前广泛应用于自动控制系统设计与分析的一种方法，它可以对系统的稳定性、可控性、可观性以及性能等进行定量分析。

在状态空间分析方法中，首先需要将系统的微分方程表示为矩阵形式的状态方程。

状态方程描述了各个状态变量和它们的变化率之间的关系。

假设系统有n个状态变量x1, x2, ..., xn和m个输入变量u1, u2, ..., um，状态方程可以表示为:dx/dt = Ax + Bu其中，dx/dt是状态变量的变化率，A是状态矩阵，描述状态变量之间的耦合关系，B是输入矩阵，描述输入变量对状态变量的影响。

状态空间分析方法的基本思想是将系统转化为状态空间表达式，然后通过对状态方程进行分析和求解来得到系统的特性和响应。

常见的分析方法包括对系统的稳定性、可控性和可观性进行评估。

稳定性是系统的基本性质之一，用来描述系统在受到扰动时是否能够恢复到平衡状态。

在状态空间方法中，通过研究系统的特征根（或特征值）可以判断系统的稳定性。

特征根是状态方程的解的根，系统的稳定性与特征根的实部有关。

如果特征根的实部都小于零，则系统是稳定的；如果特征根存在实部大于零的情况，则系统是不稳定的。

可控性是指系统是否可以通过输入变量来控制系统的状态变量。

在状态空间方法中，通过可控性矩阵来判断系统的可控性。

如果可控性矩阵的秩等于系统的状态变量个数，则系统是可控的；如果可控性矩阵的秩小于系统的状态变量个数，则系统是不可控的。

可观性是指系统的状态变量是否可以通过观测变量来测量得到。

在状态空间方法中，通过可观性矩阵来判断系统的可观性。

如果可观性矩阵的秩等于系统的状态变量个数，则系统是可观的；如果可观性矩阵的秩小于系统的状态变量个数，则系统是不可观的。

除了稳定性、可控性和可观性外，状态空间分析方法还可以用来分析系统的性能指标，如系统的响应时间、稳态误差和系统的最大误差等。

： x A(t ) x B(t )u y C (t ) x
（1）
线性时变系统的对偶系统的状态空间描述为：
d ： T AT (t ) T C T (t ) T T BT (t ) T
（2）
式中： —协状态, n维行向量； —输出, p维行向量；
如果其状态空间描述具有如下形式
ˆ ˆ ˆ ˆ x Ao x bou
其中：
0 0 0 1 1 ˆ Ao 1 n-1
ˆ ˆ y co x
ˆ co 0 0 1
则称此状态空间描述为能观测规范形。
25
总复习：现代控制理论
2．PBH秩判据
i I A rank n; C
i 1, 2, , n
3．对角线规范型判据
4．约当规范型判据
13
总复习：现代控制理论
3. 对角线规范型判据（※）
当矩阵A的特征值 1 , 2 ,, n 为两两相异时， 线性定常连续系统 x Ax x(0) x0 t0 y Cx
x (t ) L1 X ( s ) L1 (s A) 1[ x0 +B U ( s )]
9
总复习：现代控制理论
第4章 线性系统的可控性与可观测性
一、线性定常连续系统的可控性判据（※） 1．秩判据
rankQc rank B AB An 1 B n
2．PBH秩判据
rank i I A B n
i 1, 2, , n
3．对角线规范型判据 4．约当规范型判据
10
总复习：现代控制理论
3．对角线规范型判据（※）
当矩阵A的特征值 1 , 2 ,, n 为两两相异时， 线性定常连续系统 x(t ) Ax(t ) Bu (t ) x(0) x0 t 0 完全能控的充分必要条件是：其对角线规范型

1

2 n
n 1 n
t e n
1

0 1
21
n 1 2
(n 1)1 (n 1)(n 2) n 3 1 2! n2 (n 1)1 n 1 1 1
矩阵指数函数的算法 1：定义法
e At I At
1 2 2 A t 2!
只能得到eAt的数值结果，难以获得eAt解析表达式，但用计算机计算，具 有编程简单和算法迭代的优点。 2：特征值法
A P 1 AP
A PA P 1
e At Pe A t P 1
P为变换A为约当规范型的变换矩阵 1）若A的特征值为两两互异
如果系统矩阵A(t),B(t)的所有元在时间定义区间[t0,tα]上为时间t的连续实函数，输 入u(t)的所有元为时间t的连续实函数，那么状态方程的解x(t)存在且唯一。 从数学观点，上述条件可减弱为： ①系统矩阵A(t)的各个元aij(t)在时间区间[t0,tα]上为绝对可积，即：

t
t0
| aij (t ) | dt ,
－1
te1t 1t e e3t
0 2tet e 2t 1 3tet 2et 2e 2t 2 tet et e 2t
e At 0 I 1 A 2 A2 (2tet e 2t ) I (3tet 2et 2e 2t ) A (tet et e 2t ) A2 2et e 2t 0 e t e 2t 0 et 0 2et 2e 2t 0 et 2e 2t
s3 ( s 1)( s 2) 2 ( s 1)( s 2)

e At Te AtT1
A 有二种标准形式： 对角线矩阵、约当矩阵
（1）当A的特征值 1, 2,, n 为两两相异时：对角线标准型
e1t
0
e At Te AtT1 T
T 1
0
ent
其中： T为使A化为对角线标准型的非奇异变换矩阵。
求状态转移矩阵的步骤：
1） 先求得A阵的特征值 i 。
2） 求对应于 i 的特征向量 pi ，并得到T阵及T的逆阵。
j0
其中： a0(t), a1(t),, an1(t) 为t的标量函数，可按A的特 征值确定。
1） A的特征值 1, 2 ,, n 两两相异时，
注意求逆
a0(t)
1 1
12
n1 1
1
e 1t
a1 ( t )
1
1
22
n1 2
e
2t
an1 ( t )
1
1
2n
n1 n
e
nt
三、几个特殊的矩阵指数函数
（1）设 A diag 1 n ，即A为对角阵且具有互异元素时，
有
e1t
Φt
e2t
0
0
ent
（2）若A能通过非奇异变换为对角阵时，即
T -1AT Λ
e1t
ΦΦ(t)t
T
e2t
0
0
T
-1
(2-9)
ent
1
0
（3）设A为 (n n) 约当阵，即 A
a0
a1
(t (t
) )
1 1
1
1
e
1t
2
e
2t
在第3种方法中已经求得特征根，所以得：

m1m2 y (4) + [(k1 + k 2 )m1 + k1m2 ]̇ẏ + k1k2 y = k1 F (t)
3)
(2) 对式 3）进行拉氏变换可得
Y(s) =
k1
4)
F (s) m1m2 s 4 + [(k1 + k2 )m1 + k1m2 ]s2 + k1k2
·258·
(3) 对式（1）进行拉氏变换可得
(1) 令
Z (s) Y (s) G(s) = •
2)
U (s) Z(s)
式中
Z(s) =
1
3)
U (s) s3 + 7s2 + 14s + 8
Y ( s ) = s 2 + 8 s + 15
4)
Z (s)
由式 3)
̇ż̇ + 7̇ż + 14ż + 8z = u
令
x1 = z ẋ1 = x2 = ż ẋ2 = x3 = ̇ż
⎡−1 0 0 ⎤ ⎡1⎤
故有
ẋ
=
⎢ ⎢
0
−2
0
⎥ ⎥
x
+
⎢⎢1⎥⎥u
⎢⎣ 0 0 − 4⎥⎦ ⎢⎣1⎥⎦
y
=
⎡8 ⎢⎣ 3
−3 2
−
1 6
⎤ ⎥⎦
x
·261·
两种形式的状态空间表达式所对应的状态图分别如图 9-4（a），（b）所示。
8
u
1 x3
1x
1
x1
1
y
s
s
2
s
－7
－

2．状态空间的基本概念 （1）状态：系统在时间域中的行为或运动信息的集合。 （2）状态变量：能够完全表征系统运动状态的一组独立的变量，常用符号 x1（t），x2 （t），…，xn（t）表示。 （3）状态向量：由 n 个用来描述系统状态的状态变量 x1（t），x2（t），…，xn（t）组 成的向量 x（t）称为 n 维状态向量，表示为 x（t）＝[x1（t），x2（t），…，xn（t）]T。 （4）状态空间：以 n 个状态变量为基底所组成的 n 维空间。 （5）状态轨迹：系统状态在状态空间中随时间变化而形成的轨迹，又称状态轨迹。 （6）线性系统的状态空间表达式：又称为动态方程。
具有非正（负或零）实部，且具有零实部的特征值为 A 的最小多项式单根。
（2）系统的唯一平衡状态 xe＝0 是渐近稳定的充分必要条件：A 的所有特征根均具有
3．线性定常连续系统状态方程的解 （1）齐次方程求解方法：幂级数法；拉普拉斯变换法。 （2）非齐次方程求解方法：积分法；拉普拉斯变换法。
4．传递函数矩阵 表达式：G（s）＝C（sI－A）－1B＋D
二、线性系统的可控性与可观测性 1．可控性 如果系统的每一个状态变量的运动都可由输入来影响和控制，而由任意的始点达到原点， 则该系统是完全可控系统，简称为系统可控。 （1）可控标准形
5 / 75
圣才电子书 十万种考研考证电子书、题库视频学习平台

的任意初始态 x0 出发的运动轨迹 x（t；x0，t0），在 t→∞都满足：||x（t；x0，t0）－xe||≤ε，
t≥t0，则称 xe 是李雅普诺夫意义下稳定的。
（3）渐近稳定
系统不仅满足李氏意义下的稳定，且
（2）可观测性判据
3 / 75
圣才电子书 十万种考研考证电子书、题库视频学习平台

•
x2=y(1)
X n 1 X n
•
·········
X n Y (n)
xn=y(n-1)
6
其状态空间体现式为
X AX BU
Y CX
X1 X2 X= …… Xn-1 Xn
0 0 B= … 0 β0
0 1 0 ··· 0
0 0 1 ··· 0
A=
…………
0 0 0 ··· 1
-an -an-1 -an-2 ··· -a1

I )]}
S
2.Sylve假seteA如tr重P(根A)d定是d理方[e:阵stA旳Ad任j(意A多 项s式I )],且方s 阵A
具有s个相同特征值,则有:
14
求下列矩阵旳矩阵指数。
61 A=
-4 2
-3 2 A=
1 -2
01 A=
-2 -3
11 A=
-4 5
15
eAt旳某些性质
(1). d dt
17
1.无输入的情况（即U（t）=0）
X AX
1. ①只有一种状态变量A→a,且有X(0)=X0
2.
X(t)= eat ·X0
3. ②推广到多种变量,n个
4.
X(t)= eAt ·X0
5.
18
2.有输入的情况
X AX BU
(1). X aX bU
X (t) ea(tt0 ) X (t0 )
x1,x2,…,xn做轴构成旳n维空间称。 4 状态方程：状态变量旳一阶导数与状态变量、
输入量旳关系。
2
2 线性定常系统旳状态空间描述
一状态变量选用
u1
x1
输入 被控过程 状态变量 输出装置

从系统的机理出发
例 建立如右图所示机械系统的状
态空间表达式，并画出系统的状态图。
k F
根据牛顿第二定理有
m
或表示成
F ky f dy m d 2 y
dt
dt 2
d 2 y dy m dt2 f dt ky F
y f
机械位移系统
选择位移 y 和速度dy / dt为状态变量，令 x1 y, x2 dy / dt
同一个系统，状态变量的选取不是惟一的。 对于一般的物理系统，状态变量的个数应等于储能元 件的个数。
基本概念
用状态变量来表征系统时，还有如下基本概念：
状态向量 把描述系统的 n个状态变量 x1(t), x2 (t), , xn (t) 看作向
量x(t) 的分量，则x(t) 称为 n 维状态向量，记作
状态空间
x1(t)
x1
x(t
)
x2
(t
),
简记为
x
x2
xn
(t
)
xn
以状态变量 x1, x2, , xn为坐标轴所张成的n维空间，
称为状态空间。系统在任意时刻的状态，在状态空间中是一个
点，随时间推移，状态在变化，在状态空间中绘出一条轨迹，
称为状态轨线。
基本概念
状态方程 由系统的状态变量构成的一阶微分方程组，
则
x&1
x2
,
x&2
k m
x1
f m
x2
1 m
F,
y
x1
从系统的机理出发
用向量—矩阵表示的状态空间表达式为
x&1 x&2
0
k
m
1 f

4、状态轨线 的向量称为 维状态向量。
10、线性系统的状态空间表达式
状态：系统的状态是一个在时域中可以确定该系统行为或运动信息的集合。
线性系系统的统状态状空间描态述 向量在状态空间中随时间变化的轨迹。
10、线性系统的状态空间表达式
系统内部描述：状态空间描述，对系统的一种完全描述，表征系统所有动力学特征。
二、系统描述中常用的基本概念
系统状态向输量在入状态空和间中输随时出间变化：的轨由迹。外部施加到系统上的全部激励称
10、线性系统的状态空间表达式
系统状态空为间表输达式中入， 和，均是能线性从函数外。 部测量到的来自系统的信息称为
的向量称为 维状态向量。
输出。 Bellman提出了最优控制的动态规划方法。
x(t)f[x(t),u(t),t]
y(t)g[x(t),u(t),t]
x(tk1)f[x(tk),u(tk),tk] y(tk)g[x(tk),u(tk),tk]
8、自治系统
系统状态空间表达式中，函数 f 和 g不显含 t 或t k 。 x(t) f[x(t),u(t)]
y(t) g[x(t),u(t)]
谢谢观看
由输入
唯一确定，则称系统在 时刻是松弛的。
x (tk 1 ) f[x (tk)u ,(tk)tk ,]
6、输出方程 描述系统输出变量与状态变量和输入变量之间函数 关系的代数方程组。
y(t)g [x(t)u ,(t)t],
y (tk)g [x (tk)u ,(tk)tk ,]
7、状态空间表达式 状态方程与输出方程的组合称为状态空间表达式。
研究对象：多输入、多输出系统，线性、非线性、定常或时变、连续或离散系统。
线性系统的状态空间描述

第二章 线性系统的状态空间描述§2-1 状态空间的基本概念1、状态：系统的状态，是指系统的过去、现在和将来的状况。

（如：一个质点作直线运动，它的状态就是它每个时刻的位置和速度）2、状态变量：能完全表征系统运行状态的最小数目的一组变量。

（如果用最少的n 个变量x 1(t)， x 2(t)，……， x n (t)就能完全描述系统的状态，那么这n 个变量就是一组状态变量。

）3、状态向量：设一个系统有n 个状态变量，即x 1(t)，x 2(t)，……，x n (t)，用这n 个状态变量作为分量构成的向量x(t)称为该系统的状态向量。

记为Tn t x t x t x t x )](,),(),([)(21 =4、状态空间：由n 个状态变量作为坐标轴所构成的n 维空间，称为状态空间。

引入了状态和状态空间的概念之后，就可以建立动力学系统的状态空间描述了。

从结构的角度讲，一个动力学系统可用图2-1所示的方块图来表示。

其中x(t)表征系统的状态变量，u(t)为系统控制量（即输入量），y(t)为系统的输出变量。

与输入—输出描述不同，状态空间描述把系统动态过程的描述考虑为一个更为细致的过程：输入引起系统状态的变化，而状态和输入则决定了输出的变化。

5、状态方程：状态变量的一阶导数与状态变量、输入量的关系，称为系统的状态方程。

例：设单输入线性定常系统(LTI-Linear Time Invariant )的状态变量为x 1(t)，x 2(t)，……，x n (t)，输入为u(t),则一般形式的状态方程为：)()()()()()()()()()()()()()()()()()()()()(2211222221212112121111t u b t x t a t x t a t x a t x t u b t x t a t x t a t x a t x t u b t x t a t x t a t x a t x n n nn n n nn n n n ++++='++++='++++='图2-1 动力学系统结构示意图上式可写成向量—矩阵形式：其中：6、输出方程：在指定系统输出的情况下，该输出与状态变量、输入量之间的函数关系式，称为系统的输出方程。

基此，选取电容端电压uc和流经电感的电流iL作 为电路状态变量组。 显然，uc和iL必满足状态变量定义中所指出的线 性无关极大组属性。
2013/11/22
线性系统理论
23
第2章 线性系统的状态空间描述
广东工业大学 自动化学院 自动控制系 陈玮
1、电路系统状态空间描述的列写示例
R1 C + uC iL R2 iC + -
由此，得 和，
X QX
X PX PQX
X QX QPX
显然， PQ QP I
即矩阵P和Q互逆。
结论：系统的任意选取的两个状态X和 X 之间 为线性非奇异变换的关系。
2013/11/22
线性系统理论
18
第2章 线性系统的状态空间描述
广东工业大学 自动化学院 自动控制系 陈玮
2013/11/22
线性系统理论
2
第2章 线性系统的状态空间描述
广东工业大学 自动化学院 自动控制系 陈玮
一、状态和状态空间
1、系统动态过程的两类数学描述 2、状态和状态空间的定义
2013/11/22
线性系统理论
3
第2章 线性系统的状态空间描述
广东工业大学 自动化学院 自动控制系 陈玮
1、系统动态过程的两类数学描述
线性系统理论
21
第2章 线性系统的状态空间描述
广东工业大学 自动化学院 自动控制系 陈玮
1、电路系统状态空间描述的列写示例
电路系统如图所示，设各组元件的参数值为已 知，取电压源e(t)为输入变量，电阻R2端电压uR2为输 出变量。 C
R1 iC + e (t ) -
L
+ uC iL R2

现代控制理论实验指导书实验一：线性系统状态空间分析1、模型转换图1、模型转换示意图及所用命令传递函数一般形式：)()(1111110nmasasasabsbsbsbsGnnnnmmmm≤++++++++=----MATLAB表示为：G=tf(num,den)，其中num,den分别是上式中分子，分母系数矩阵。

零极点形式：∏∏==--=nijmiipszsKsG11)()()(MATLAB表示为：G=zpk(Z,P,K)，其中Z，P，K分别表示上式中的零点矩阵，极点矩阵和增益。

传递函数向状态空间转换：[A,B,C,D] = TF2SS(NUM,DEN)；状态空间转换向传递函数：[NUM,DEN] = SS2TF(A,B,C,D,iu)---iu表示对系统的第iu个输入量求传递函数；对单输入iu为1;验证教材P438页的例9-6。

求P512的9-6题的状态空间描述。

>> A=[0 1;0 -2];>> B=[1 0;0 1];>> C=[1 0;0 1];>> D=[0 0;0 0];>> [NUM,DEN] = ss2tf(A,B,C,D,1)NUM =0 1 20 0 0DEN =1 2 0>> [NUM,DEN] = ss2tf(A,B,C,D,2)NUM =0 0 10 1 0DEN =1 2 0给出的结果是正确的，是没有约分过的形式P512 9-6>> [A,B,C,D]=tf2ss([1 6 8],[1 4 3])A =-4 -31 0B =1C =2 5D =12、状态方程求解单位阶跃输入作用下的状态响应：G=ss(A,B,C,D);[y,t,x]=step(G);plot(t,x).零输入响应[y,t,x]=initial(G,x0)其中，x0为状态初值。

验证P435的例9-4，P437的例9-5。

9-4A=[0 1;-2 -3];B=[0;0];C=[0 0];D=[0];G=ss(A,B,C,D);[y,t,x]=initial(G,[1;2]);plot(t,x)（设初始状态为[1 ;2]）零输入响应00.20.40.60.81 1.2 1.4 1.6 1.82-1-0.50.511.529-5零输入响应A=[0 1;-2 -3];B=[0;1];C=[0 0];D=[0];G=ss(A,B,C,D);[y,t,x]=initial(G,[1;2]);plot(t,x)00.20.40.60.81 1.2 1.4 1.6 1.82-1-0.50.511.52零状态响应，阶跃信号激励下>> A=[0 1;-2 -3];B=[0;1];C=[0 0];D=[0];>> G=ss(A,B,C,D);[y,t,x]=step(G);plot(t,x)00.20.40.60.81 1.2 1.4 1.6 1.8200.050.10.150.20.250.30.350.4总响应>> A=[0 1;-2 -3];B=[0;1];C=[0 0];D=[0];G=ss(A,B,C,D);[y1,t1,x1]=step(G);[y2,t2,x2]=initial(G,[1;2]);>> x=x1+x2;>> plot(t1,x)00.20.40.60.81 1.2 1.4 1.6 1.82-0.500.511.523、系统可控性和可观测性可控性判断：首先求可控性矩阵：co=ctrb(A ，B)。

线性系统的状态空间分析与综合第九章线性系统的状态空间分析与综合⼀、教学⽬的与要求：通过本章内容的学习，使学⽣建⽴起状态变量和状态空间的概念，掌握线性定常系统状态空间模型的建⽴⽅法，状态空间表达式的线性变换，状态完全能控或状态完全能观测的定义，及其多种判据⽅法，状态转移矩阵的求法，传递函数矩阵与状态空间表达式的关系。

⼆、授课主要内容：1．线性系统的状态空间描述2．线性系统的可控性与可观测性3．线性定常系统的状态反馈与状态观测器（详细内容见讲稿）三、重点、难点及对学⽣的要求（掌握、熟悉、了解、⾃学）1.重点掌握线性定常系统状态空间模型的建⽴⽅法与其他数学描述（微分⽅程、传递函数矩阵）之间的关系。

2.掌握采⽤状态空间表述的系统运动分析⽅法，状态转移矩阵的概念和求解。

3.掌握系统基本性质——能控性和能观测性的定义、有关判据及两种性质之间的对偶性。

4.理解状态空间表达式在线性变换下的性质，对于完全能控或能观测系统，构造能控、能观测标准形的线性变换⽅法，对于不完全能控或不完全能观测系统，基于能控性或能观测性的结构分解⽅法。

5.掌握单变量系统的状态反馈极点配置和全维状态观测器设计⽅法，理解分离定理，带状态观测器的状态反馈控制系统的设计。

重点掌握线性系统的状态空间描述和求解，线性系统的可控性与可观测性及状态反馈与状态观测器。

四、主要外语词汇线性系统 linear system状态空间 state space状态⽅程 state equation状态向量 state vector传递函数矩阵 translation function matrix状态转换矩阵 state-transition matrix可观测标准形 observational standard model可控标准形 manipulative standard model李亚普诺夫⽅程Lyaponov equation状态观测器 state observation machine对偶原理 principle of duality五、辅助教学情况（见课件）六、复习思考题1．什么是系统的状态空间模型？状态空间模型中的状态变量、输⼊变量、输出变量各指什么？2．通过机理分析法建⽴系统状态空间模型的主要步骤有哪些？3．何为多变量系统？如何⽤传递矩阵来描述多变量系统的动态特性？在多变量系统中，环节串联、并联、反馈连接时，如何求取总的传递矩阵？4．试简述数学模型各种表达式之间的对应关系。

则有 x1 x2
x2 x3
x3 a0 x1 a1x2 a2 x3 b0u
写成矩阵形式
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
选择 n 个状态变量为 系统方程为
x1 y 0u x2 x1 1u x3 x2 2u
xn xn1 n1u
x1
x2
xn
0
0
0
a0
1 0 0 a1
0 1 0 a2
0 0 0 a3
0 0 1 an1
x1 x2
xn
1
x3 x4
0 0 0
0 0 0
mg M
0
(M m)g Ml
0
x2
1 M
u
;
1 0
x3 x4
0
1 Ml
x1
y 1
0
0
0
x2
x3 x4
状态图为
1.2 由微分方程求状态空间表达式
一个系统，用线性定常微分方程描述其输入和输出的关系。通过选 择合适的状态变量，就可以得到状态空间表达式。
设小球的重心坐标为： ( yG , zG )
则 yG y l sin
zG l cos
在水平方向，应用牛顿第二定律：
M
d2 y dt2
m
d2 dt2
(y
l
sin )
u
转动方向的力矩平衡方程式：