线性系统状态空间分析报告与运动解
[计算机软件及应用]2 线性系统运动分析
![[计算机软件及应用]2 线性系统运动分析](https://img.taocdn.com/s3/m/3f8e017ebe23482fb4da4c7e.png)
1 0 2t e 0
1 1 2 1 0 2
1 (1 e 2t ) 2 e 2t
方法三
由于
s sI A 0 0 0 s 0 1 s 2 0
1 1 1
0 1 2
0 1 2
方法四：化eAt为A的有限项法(Caley-Hamilton定理法) 利用凯莱-哈密尔顿定理，化为A的有限项，然后 通过求待定时间函数获得的方法。 必须指出，这种方法相当系统，而且计算过程简单。 设A的最小多项式阶数为m。可以证明，采用赛尔维斯特 内插公式，通过求解行列式 1 e
2.1 状态方程的齐次解 所谓系统的自由解,是指系统输入为零时,由初始 状态X(0)引起的自由运动.即零输入响应.
x(t ) Ax(t )
若
x(t 0 ) x 0
则
x(t ) (t; t 0 , x 0 ,0) e A(t t 0 ) , t t 0
x(0) x 0 若 则 x(t ) (t;0, x 0 ,0) e At , t 0 证明: 当 x(0) x 0 时,对矩阵方程进行拉氏变换,得
方法一：直接计算法(矩阵指数函数)
e
At
A2t 2 A3 t 3 I At 2! 3!
1 k k A t k 0 k!
(2.9)
可以证明，对所有常数矩阵A和有限的t值来说， 这个无穷级数都是收敛的。
方法二：对角线标准形与Jordan标准形法 1)若可将矩阵A变换为对角线标准形，
s sI A 0
( sI A)1
(t ) e At L1[(sI A) 1 ]
线性系统理论3线性系统的运动分析
THANKS
伯德图判据
通过观察系统开环伯德图（对数幅频特性和相频特性曲线）来判断系统的稳定性。若开环伯 德图在穿越频率处的相位裕度大于0，则系统是稳定的。
不稳定系统的分析与处理
不稳定原因分析
不稳定系统可能由于系统内部参数摄动、外部扰动或控 制器设计不当等原因导致。需要对系统进行详细分析， 找出不稳定的原因。
不稳定系统处理
线性微分方程
01
描述线性系统动态行为的数学工具，通过求解微分方程可以得
到系统的输出响应。
传递函数
02
在频域中描述线性系统输入输出关系的数学表达式，常用于控
制系统的分析和设计。
状态空间方程
03
描述线性系统状态变量和输入输出关系的数学方程组，适用于
多输入多输出系统和时变系统。
线性系统的建模方法
1 2
机理建模
运动方程的物理意义
描述系统运动状态
运动方程描述了线性系统的运动状态，包括位置、速度和 加速度等物理量。通过求解运动方程，可以得到这些物理 量的时域解和频域解。
预测系统响应
根据已知输入和初始条件，通过求解运动方程可以预测线 性系统的响应。这对于控制系统的设计和分析具有重要意 义。
分析系统稳定性
通过分析运动方程的解的性质，可以判断线性系统的稳定 性。例如，如果解是收敛的，则系统是稳定的；如果解是 发散的，则系统是不稳定的。
对求解结果进行可视化展示和数据分 析，研究电路系统的动态响应特性， 如谐振频率、阻尼振荡等。
建立模型
运动方程
求解方法
结果分析
根据电路元件的连接方式和电气特性， 建立电路系统的数学模型，如RLC串 联或并联电路。
采用解析法或数值法求解运动方程， 得到电路中各元件的电压、电流等电 气参数。
第三章线性系统的运动分析
Chapter 3 Analysis of Linear System3.1 INTRODUCTION运动分析的数学实质：从数学的角度，运动分析的实质就是求解系统的状态方程。
以解析形式或数值分析形式，建立系统状态随输入和初始状态的演化规律。
(Solving the time-invariant state equation)3.2 连续时间线性时不变系统的运动分析SOLVING THE TIME-INVARIANT STATE EQUATION系统响应=系统的零输入响应+系统的零状态响应System response=a term consisting of the transition of the initial state +a term arising from the input vector零输入响应：自由运动，由系统矩阵决定，不受外输入影响。
零状态响应：强迫运动，响应稳态时具有和输入相同的函数形态。
01!k k ∞−+=∑0k k b t ∞=+=∑2012Ab Ab t Ab t +=+++b k 0)b +Equating the coefficients of the equal powers of t, we obtain By substituting this assumed solution in to Equation (1)解的说明：1.零输入响应是状态空间中由初始状态经线性变换矩阵所导出的一个变换点。
2.自由运动3.自由运动的轨迹由唯一决定。
4.当自由运动轨迹趋于平衡状态时，则系统是渐近稳定的。
At e0x Ate 0=x若初始时间取为t 0≠0则0)(,)(0t t x e t x t t A ou ≥=−00)(x t x =01!k k ∞−+=∑+232322332323332)()2!3!F F I Ft t t F t A t A Ft AF t F t ++++++0+=0,1,2,))AtAt Ae A e A ++=+=利用性质+λ)neλ)n t0000i i λλ⎤⎥⎥⎥⎥⎥⎥⎦12)l J t J tJ t e e 0i i t t e e e λλ⎤⎥⎥⎥⎥⎥⎥⎥⎦系统状态运动规律的基本表达式设系统的状态空间描述为有表达式⎰⎰≥−+=+=−t A Att t A At t d t Bu e x e d Bu e x e t x 000)(00,)(,)()(ττττττ⎰≥+=−−t t t A t t A t t d Bu e x e t x 000)(0)(,)()(τττ对初始时刻t 0=0 情形有表达式注意：物理意义解的讨论:(1)卷积特征；(2)零初始响应的几何特征；(3)可达性；(4)任意时刻的表达式00≥,=)(),(+=t t x t x t Bu Ax x3.3连续时间线性时不变系统的状态转移矩阵State-Transition Matrix设连续时间线性时不变系统，状态方程为：as To verify this, note thatWe thus confirm that Equation (2) is the solution of Equation (1))2()0()()(x t t x Φ=where )(Φt is n n ⨯Matrix and is the unique solution of)0()0()0()0(x x x =Φ=Ate t =)(Φ)(=)0()(Φ=)0()(Φ=)(t Ax x t A x t t xI t A t =)0(Φ)(Φ=)(Φ )1(=Ax x and状态转移矩阵的形式为()()()0000,0000t t e t t t t e t t t t A At ≥=−Φ≠≥=Φ=−时，时，基于状态转移矩阵的系统响应表达式()()()()()()()()()⎰⎰−Φ+−Φ=≥−Φ=−Φ=tt t t ox ou d Bu t x t t t x t t d Bu t t x x t t t x 0000000ττττττ。
第1章线性系统的状态空间描述
x&(t) Ax(t) Bu(t) y(t) Cx(t) Du(t)
• 情况1：输入u不含导数
y(n) an1y(n1) L a1y& a0 y bu
自主技术与智能控制研究中心
二、状态空间模型的建立
输入u不含导数 y(n) an1y(n1) L a1y& a0 y bu
选取状态变量 x1 y x2 y x3 y
I ml2 ml
自主技术与智能控制研究中心
ml M m
二、状态空间模型的建立
用一阶微分方程组表示系统模型！
&x& 1m2l2 g 1(I ml2 )u && 1(M m)mgl 1mlu
引入新的变量
x1 x x2 x&
x3 x4 &
x&1 x2
x&2 x&3
{1m2l x4
x&% Ax% Bu%
y%
Cx%
Du%
f1
A
f x
|x0
,u0
x1
M
fn x1
L O L
f1 xn
M
fn xn
B
f u
|x0 ,u0
,C
g x
|x0
,u0
,
D
g u
|x0 ,u0
自主技术与智能控制研究中心
二、状态空间模型的建立 例3：质量-弹簧-阻力器系统
自主技术与智能控制研究中心
u
线性化 0 V mg
m
d2 dt 2
(x
l
sin )
H
I&& Vl Hl
m d 2 (l cos) V mg
第三章线性系统状态方程的解
第三章 线性系统的运动分析§3-1线性连续定常齐次方程求解一、齐次方程和状态转移矩阵的定义１、齐次方程状态方程的齐次方程部分反映系统自由运动的状况（即没有输入作用的状况），设系统的状态方程的齐次部分为：)()(t Ax t x= 线性定常连续系统：Ax x =２、状态转移矩阵的定义齐次状态方程Ax x = 有两种常见解法：（1）幂级数法；（2）拉氏变换法。
其解为)0()(x e t x At ⋅=。
其中Ate 称为状态转移矩阵（或矩阵指数函数、矩阵指数），记为：At e t =)(φ。
若初始条件为)(0t x ，则状态转移矩阵记为：)(00)(t t A e t t -=-Φ对于线性时变系统，状态转移矩阵写为),(0t t φ，它是时刻t ，t 0的函数。
但它一般不能写成指数形式。
（１）幂级数法设Ax x= 的解是t 的向量幂级数 +++++=kk t b t b t b b t x 2210)(式中 ，，，，，k b b b b 210都是n 维向量，则+++++=-1232132)(k k t kb t b t b b t x)(2210 +++++=kk t b t b t b b A故而有：⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧======00323021201!1!31312121b A k b b A Ab b b A Ab b Ab b K K2且有0)0(b x =。
故+++++=kk t b t b t b b t x 2210)(+++++=k k t b A k t b A t Ab b 020200!1!21 )0()!1!21(22x t A k t A At I kk +++++=定义：∑∞==+++++=022!1!1!21K k k k k Att A k t A k t A At I e则)0()(x e t x At ⋅=。
（２）拉氏变换解法将Ax x= 两端取拉氏变换，有 )()0()(s Ax x s sx =- )0()()(x s x A sI =- )0()()(1x A sI s x ⋅-=-拉氏反变换，有)0(])[()(11x A sI L t x ⋅-=-- 则])[()(11---==A sI L e t At φ【例3.1.1】 已知系统的状态方程为x x⎥⎦⎤⎢⎣⎡=0010 ，初始条件为)0(x ，试求状态转移矩阵和状态方程的解。
线性系统理论 第2章 线性系统的状态空间描述
u(k )
H (k )
x(k 1)
x(k )
单位延迟
C (k )
y(k )
G (k )
7/7,11/50
2.3.连续变量动态系统按状态空间描述的分类
线性系统和非线性系统
设系统的状态空间描述为 x f ( x, u, t ) y g ( x, u, t )
向量函数
g1 ( x, u, t ) f1 ( x, u, t ) g ( x, u , t ) f ( x, u , t ) ，g ( x, u, t ) 2 f ( x, u , t ) 2 g q ( x, u , t ) f n ( x, u , t )
和t≥t0 各时刻的任意输入变量组 u1 (t ),u2 t ,, u p (t ) 那么系统的任何一个内部变量在t≥t0各时刻的运动行为也就随之而完全确定
3/4,3/50
(2).状态变量组最小性的物理特征： 少一个不行，多一个没用 (3). 状态变量组最小性的数学特征：极大线性无关变量组 (4). 状态变量组的不唯一性 ：任意
1/18,14/50
结论1
给定单输入,单输出线性时不变系统的输入输出描述,
y ( n) an1 y ( n1) a1 y (1) a0 y bmu ( m) bm1u ( m1) b1u (1) b0u
Y (s) bm s m bm1 s m1 b1 s1 b0 g ( s) U ( s) s n an1 s n1 a1 s a0
时变系统和时不变系统
f f ( x, u ) 若向量f,g不显含时间变量t,即 g g ( x, u )
第一章 线性定常系统的状态空间描述及运动分析
称 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
线性系统状态空间分析和运动解
线性系统状态空间分析和运动解状态空间分析方法是一种用来描述线性系统的分析方法。
它将系统的动态特性用一组状态变量来表示，并通过矩阵形式的状态方程进行分析和求解。
状态空间方法是目前广泛应用于自动控制系统设计与分析的一种方法，它可以对系统的稳定性、可控性、可观性以及性能等进行定量分析。
在状态空间分析方法中，首先需要将系统的微分方程表示为矩阵形式的状态方程。
状态方程描述了各个状态变量和它们的变化率之间的关系。
假设系统有n个状态变量x1, x2, ..., xn和m个输入变量u1, u2, ..., um，状态方程可以表示为:dx/dt = Ax + Bu其中，dx/dt是状态变量的变化率，A是状态矩阵，描述状态变量之间的耦合关系，B是输入矩阵，描述输入变量对状态变量的影响。
状态空间分析方法的基本思想是将系统转化为状态空间表达式，然后通过对状态方程进行分析和求解来得到系统的特性和响应。
常见的分析方法包括对系统的稳定性、可控性和可观性进行评估。
稳定性是系统的基本性质之一，用来描述系统在受到扰动时是否能够恢复到平衡状态。
在状态空间方法中，通过研究系统的特征根（或特征值）可以判断系统的稳定性。
特征根是状态方程的解的根，系统的稳定性与特征根的实部有关。
如果特征根的实部都小于零，则系统是稳定的；如果特征根存在实部大于零的情况，则系统是不稳定的。
可控性是指系统是否可以通过输入变量来控制系统的状态变量。
在状态空间方法中，通过可控性矩阵来判断系统的可控性。
如果可控性矩阵的秩等于系统的状态变量个数，则系统是可控的；如果可控性矩阵的秩小于系统的状态变量个数，则系统是不可控的。
可观性是指系统的状态变量是否可以通过观测变量来测量得到。
在状态空间方法中，通过可观性矩阵来判断系统的可观性。
如果可观性矩阵的秩等于系统的状态变量个数，则系统是可观的；如果可观性矩阵的秩小于系统的状态变量个数，则系统是不可观的。
除了稳定性、可控性和可观性外，状态空间分析方法还可以用来分析系统的性能指标，如系统的响应时间、稳态误差和系统的最大误差等。
线性系统的运动分析
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种方法中已经求得特征根，所以得：
现代控制工程-第二章线性系统的状态空间描述
1 x3 s
1 s
1 x1 s
y(t )
2
3
8 64
解：第一步：化简方框图，使得整个系统只有标准积分器（1/s）、 比例器（k）及加法器组成。 第二步：将上述调整过的结构图中的每个标准积分器（1/s） 的输出作为一个独立的状态变量xi，积分器的输入端就是状态变 量的一阶导数dxi/dt。 第三步：写出每个状态变量的一阶微分方程，从而写出系统 的状态方程。
y Cx Du
图2－2 系统动态方程的方块图结构
状态空间分析法具有下列优越之处：
便于在数字计算机上求解；
容易考虑初始条件； 能了解并利用处于系统内部的状态信息； 数学描述简化；
适于描述多输入－多输出、时变、非线性、随机、离散等各类 系统，是最优控制、最优估计、辨识、自适应控制等现代控制系 统的基本描述方法。
例2.2.3求如图所示系统的动态方程。
（a）系统方块图
u(t )
s 1 s2
1 s3
1 s 2 8s 64
y(t )
（b）第一次等效变换
1 s3
u(t )
1 s2
1 s( s 8)
y(t )
64
（c）由标准积分器组成的等效方块图
u(t )
1 x4 s
(2-5)
y t cx t du(t )
，cn ，d为直接联系输入量、输出量 其中 c c1，c2， 的前向传递（前馈）系数，又称前馈系数。
多输入－多输出（含q个输出变量）线性定 常连续系统的输出方程一般表达形式为：
y1 c11 x1 c1n xn d11u1 d1 pu p yq cq1 x1 cqn xn d q1u1 d qp u p
现代控制理论第2章 线性系统的运动
定义矩阵向量eAt为状态转移矩阵
于是齐次状态方程的解为:
x(t ) e x(0) e x0
At At
(2 7)
另用拉氏变换法求解齐次微分方程:
x(t ) Ax(t ) sx(s) x(0) Ax(s)
x(s) (sI A) x(0)
1
拉氏反变换后得到齐次状态方程的解:
2. 状态转移矩阵的计算。
a. 直接求取；
b. 拉普拉斯变换；
c. 化矩阵A为对角型或约当型；
d.化矩阵指数 e At
为A的有限项。
（1）线性系统状态转移矩阵的基本性质
1 2 2 1 i i Φ(t ) e I At A t A t (2 9 1) 2! i!
表明 (t ) 具有分段组合的性质。
④ Φ1 (t ) Φ(t ) , Φ1 (t ) Φ(t ) 证：根据性质①和③及逆矩阵定义，有
Φt t Φt Φ t I Φ t t Φ t Φt I
Φ 1 (t ) Φ(t )
At
t0 0
x(t ) (t t0 )x(t0 )
x(t ) (t )x(0)
(2 10)
(2 9)
状态转移矩阵 (t )包含了系统自由运动的全部信息， 完全表征了系统的动态特性,A的状态矩阵唯一。
几点解释：
① 如果t为某给定常数T，那么零输入响应 x(t ) 就是状 态空间中由初始状态 x 0 经线性变换常数阵 (t ) 所致。
（一）齐次状态方程解的一般表达式
x(t ) A(t )x(t ),
x(t0 ) x0 ,
2
t t0
i
b 2t bi t
线性系统的状态空间分析与综合PPT资料(正式版)
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、状态空间表达式 状态方程与输出方程的组合称为状态空间表达式。
研究对象：多输入、多输出系统，线性、非线性、定常或时变、连续或离散系统。
线性系统的状态空间描述
现代控制理论实验报告
现代控制理论实验指导书实验一：线性系统状态空间分析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)。
现代控制理论-线性系统的状态空间描述
c11(t) c12 (t) c1n (t)
C
(t)
c21
(
t
)
c22 (t)
c2n (
t
)
,
m n维输出矩阵 表 征 输 出 和 每 个 状 态 量 变 的 关 系
cm1(t) cm2 (t) cmn (t)
d11(t)
D(t)
d 21 ( t )
d12 (t)
d22 (t)
最小个数：意味着这组变量是互相独立的。一个用n阶微分方
程描述的含有n个独立变量的系统，当求得n个独立变量随时
间变化的规律时，系统状态可完全确定。若变量数目多于n，
必有变量不独立；若少于n，又不足以描述系统状态。
2021/8/24
电气信息学院《现代控制理论课程》
12
状态变量的选取具有非唯一性，即可 用某一组、也可用另一组数目最少的变量 （状态变量不唯一）。状态变量不一定要 象系统输出量那样，在物理上是可测量或 可观察的量，但在实用上毕竟还是选择容 易测量的一些量，以便满足实现状态反馈、 改善系统性能的需要。
常用符号： 积分器
比例器
ki
注：有几个状态变量，就建几个积分器
加法器
注:负反馈时为－
系统框图：
U
B
D
•
X
A
X C Y
X•
AX
BU
Y CX D U
2021/8/24
电气信息学院《现代控制理论课程》
22
线性时变系统状态空间描述：x A(t)x B(t)u y C(t)x D(t)u
t)
b11(t) b12 (t) b1r (t)
B(t)
b21 ( t )
线性系统的状态空间分析与综合
线性系统的状态空间分析与综合第九章线性系统的状态空间分析与综合⼀、教学⽬的与要求：通过本章内容的学习，使学⽣建⽴起状态变量和状态空间的概念，掌握线性定常系统状态空间模型的建⽴⽅法，状态空间表达式的线性变换，状态完全能控或状态完全能观测的定义，及其多种判据⽅法，状态转移矩阵的求法，传递函数矩阵与状态空间表达式的关系。
⼆、授课主要内容：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．试简述数学模型各种表达式之间的对应关系。
线性系统的状态空间描述
状态向量：是由状态变量所构成的向量，即向量
x(t)x1(t),x2(t),L,xn(t)T称为n维状态向量。
状态空间：以n个线性无关的状态变量作为基底所组 成的 n 维空间称为状态空间Rn。
状态轨线：随着时间推移，系统状态x(t)在状态空间 所留下的轨迹称为状态轨线或状态轨迹。
连续系统：
x&(t)f [x(t),u(t),t] y(t)g[x(t),u(t),t]
离散系统：
xy(t(ktk1))gf[[xx(t(ktk),),uu(t(ktk),)t,ktk]] 或 x(yk(k)1)g[fx[(xk()k,)u,(uk()k,)k,]k]
4．线性系统状态空间表达式：状态方程与输出方 程都是线性方程的系统是线性系统。线性系统的状态方 程是一阶向量线性微分方程或一阶向量线性差分方程。
关于状态的几点说明
系统的状态空间描述
状态变量组选取上的不唯一性： 由于系统中变量的个数必大于n，而其中仅有n个
是线性无关的，因此决定了状态变量组在选取上的不 唯一性。
➢状态变量不是所有变量的总和。 ➢输出量可以选作状态变量。 ➢输入量不允许选作状态变量。 ➢状态变量有时是不可测量的。 ➢状态变量是时间域的。
对于控制工程而言，它可能是被控对象、控 制装置，也可能是某些部件的串联、并联和反馈 组合。
图1-1 系统的方块图表示
✓ 图中方块以外的部分为系统环境; ✓ 环境对系统施加的作用或激励称为系统输入，
用向量 u[u1,u2,Lup]T表示; ✓ 系统对环境的作用（即从外部量测到的系统信
息）称为系统输出，用向量 y[y1,y2,Lyq]T表示； ✓ 系统输入和输出统称为系统的外部变量。 ✓ 描述系统内部状况的变量称为系统的状态变量，
线性系统的状态空间描述
状态向量：是由状态变量所构成的向量，即向量 称为n维状态向量。 状态空间：以n个线性无关的状态变量作为基底所组成的 n 维空间称为状态空间Rn。 状态轨线：随着时间推移，系统状态x(t)在状态空间所留下的轨迹称为状态轨线或状态轨迹。
状态方程（※）：描述系统状态变量与输入变量之间关系的一阶微分方程组(连续时间系统)或一阶差分方程组(离散时间系统)称为系统的状态方程。 状态方程表征了系统由输入所引起的内部状态变化，其一般形式为： 或 线性系统的状态空间描述
对角型实现和约当标准型实现，需要计算系统的极点(特征值)和特征向量，很不方便。
总结：
由系统微分方程建立状态空间表达式（自学P405-409）
01
由系统微分方程建立状态空间表达式的整个思路与由系统传递函数建立状态空间表达式的思路是类似的，所以这里不再详细介绍，请参看教材P405-407。 另外，当给定系统微分方程时，可先求出其传递函数，然后按照前面推导的公式直接写出其可控标准型和可观测标准型实现，例如我们在例1-2种所做的那样。
状态变量组选取上的不唯一性： 由于系统中变量的个数必大于n，而其中仅有n个是线性无关的，因此决定了状态变量组在选取上的不唯一性。
系统的状态空间描述
系统的任意选取的两个状态变量组之间为线性非奇异变换的关系。
状态变量是时间域的。
状态变量有时是不可测量的。
状态变量不是所有变量的总和。
1.5 组合系统的状态空间描述
1.4 线性系统等价的状态空间描述
1.1 线性系统状态空间描述
2021
2023
1.1 线性系统状态空间描述
一．系统数学描述的基本类型
1．几个基本定义
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
【实验地点】课外(宿舍) 【实验目的】1、学会利用MATLAB 实现离散系统传递函数模型的生成2、学会利用MATLAB 将连续系统离散化 【实验设备与软件】1、MATLAB/Simulink 数值分析软件2、计算机一台 【实验原理】1、求矩阵特征值和特征向量命令格式[V J]=eig （A ） Cv=eig(A)说明：V 特征向量，J 是Jordan 型，cv 是特征值列向量 2、求运动的方法（1）利用Laplace 逆变换----适合于连续/离散线性系统采用ilaplace/iztrans 对传递函数求逆，这种方法一般是零输入情况下求响应。
（2）用连续（离散）状态转移矩阵表示系统解析解----适合于线性定常系统 对连续定常系统有：假设初始时刻为零，LTI 系统的解析解为dt Bu ee x e t x tAtAtAt⎰⎰+=0)()0()(τ。
若u （t ）是单位阶跃输入，则上述解可写成dtBu ee x e t x tAtAt At ⎰⎰+=0)()0()(τ。
进一步简化为：Bu A Bu A x e t x At 11))0(()(---+=对离散线性定常系统有：∑---+=11)()0()(k i kki Hu G x G k x（3）状态方程的数值分析方法----适合于连续线性系统和非线性系统采用直接数值积分很容易的处理各种定常/时变和线性/非线性系统。
有很多数值积分方法，其中有一类预测-修正数值积分方法+自适应步长调整的算法比较有效。
在MATLAB/Simulink 中包含的多种有效的、适用于不同类型的ODE 求解算法，典型的是Runge-Ktuta 算法，其通常使用如下的函数格式：[t,x]=ode45(odefun,[ti,tf],x0,options)----采用四阶、五阶Runge-Ktuta 算法 [t,x]=ode23(odefun,[ti,tf],x0,options)----采用二阶、三阶Runge-Ktuta 算法 说明：a.这两个函数是求解非刚性常微分方程的函数。
b.参数options 为积分的误差设置，取值为相对误差‘reltol ’和绝对误差‘abstol ’;[ti,tf]求解的时间围；x0是初值是初值向量；[t,x]是解。
（4）利用CotrolToolBox 的离散化求解函数----适合于TLI 系统 用step （）/impulse()函数求取阶跃输入/冲激输入时系统的状态响应： 当系统G 是连续的情况下：调用[y,t,x]=step/impulse(G )会自动对连续系统G 选取采样时间围和周期；调用[y,t,x]=step/impulse(G ，ti:Ts:tf)由用户自己定义对连续系统G 的样时间围和周期； 当系统G 是离散的情况下:调用[y,t,x]=step/impulse(G )会按离散系统G 给出的采样周期计算；调用[y,t,x]=step/impulse(G ，ti:Ts:tf)是Ts 必须与离散系统G 的采样时间围和周期一致。
另外lsim()函数调用格式：[y,x,t]=lsim(G,u,ti,TS,tf,x0) 零输入响应调用函数initial （）,格式：[y,x,t]=(G,x0) (5)利用simulink 环境求取响应----适用于所有系统求取响应 使用simulink 求取线性或非线性系统的响应，调用格式如下：[t,x,y]=sim(‘XX.mdl ’,ti:Ts:tf,options,u)【实验容】已知线性系统：])(201)()(210)(404040202119201921)(t x t y t u t x t x +-----•已知线性系统1、利用Matlab 求零状态下的阶跃响应（包括状态和输出），生成两幅图：第一幅绘制各状态响应曲线并标注；第二幅绘制输出响应曲线。
状态响应曲线：A=[-21 19 -20;19 -21 20;40 -40 -40]; B=[0;1;2]; C=[1 0 2];D=[0]; %输入状态空间模型各矩阵,若没有相应值,可赋空矩阵 X0=[0;0;0]; % 输入初始状态 sys=ss(A,B,C,D); %构造传递函数[y,x,t]=step(sys); % 绘以时间为横坐标的状态响应曲线图 plot(t,x); grid;title('状态响应曲线') 输出响应程序：A=[-21 19 -20;19 -21 20;40 -40 -40]; B=[0;1;2]; C=[1 0 2]; D=0; X0=[0;0;0][num,den]=ss2tf(A,B,C,D,1); sys=tf(num,den); step(sys) gridtitle('输出响应曲线')图一（状态响应曲线）图二（输出响应曲线）2、利用Matlab 求零状态下的冲激响应（包括状态和输出），生成两幅图：第一幅绘制各状态响应曲线并标注；第二幅绘制输出响应曲线。
状态响应曲线程序:A=[-21 19 -20;19 -21 20;40 -40 -40]; B=[0;1;2]; C=[1 0 2];D=[]; %输入状态空间模型各矩阵,若没有相应值,可赋空矩阵 x0=[0;0;0]; % 输入初始状态 sys=ss(A,B,C,D); %构造传递函数 [y,x,t]= impulse(sys); plot(t,x); grid;title('状态响应曲线') 输出响应曲线程序:A=[-21 19 -20;19 -21 20;40 -40 -40]; B=[0;1;2]; C=[1 0 2]; D=0; X0=[0;0;0][num,den]=ss2tf(A,B,C,D,1); sys=tf(num,den); impulse(sys); grid; title('')图三（状态响应曲线'）图四（输出响应曲线）3、若控制输入为，且初始状态为，求系统的响应，要求a.在simulink只能够画出模型求响应，生成两幅图：第一幅绘制各状态响应曲线并标注；第二幅绘制输出响应曲线。
程序如下：t=[0:0.01:5];u=(1+exp(-t).*cos(5*t)).*(t<3)+1*(t>=3);t=t';u=u';ut=[t,u];[t1,x,y]=sim('shiyan5.mdl',t,[],ut);plot(t1,x)figure(2);plot(t1,y)创建的模型图如下：图五（模型图）b.编写.m文件求响应，生成两幅图：第一幅绘制各状态响应曲线并标注；第二幅绘制输出响应曲线。
状态响应曲线：t=[0:0.02:5];u=(1+exp(-t).*cos(5*t)).*(t<3)+1*(t>=3);t=t';u=u';A=[-21 19 -20;19 -21 20;40 -40 -40];B=[0;1;2];C=[1 0 2];D=[0]; %输入状态空间模型各矩阵,若没有相应值,可赋空矩阵 X0=[0.2;0.2;0.2]; % 输入初始状态u=(t==0); %就是个条件判断，只有t=0的时候，u才为“1”sys=ss(A,B,C,D); %构造传递函数plot(t,x);grid;title('状态响应曲线')输出响应曲线：plot(t,y);grid;title('输出响应曲线')图六（状态响应曲线）图七（输出响应曲线）4、以阶跃输入情况下的，分析各模块对响应有什么影响。
图八（阶跃输入时）阶跃输入的图像到答稳定时间快，曲线平滑5、求系统的传递函数在MATLAB软件Command Window窗口中输入以下程序A=[-21 19 -20;19 -21 20;40 -40 -40];B=[0;1;2];C=[1 0 2];D=0;[num,den]=ss2tf(A,B,C,D,1);printsys(num,den)程序运行结果为图七6、若采用K增益负反馈，绘制闭环根轨迹图，并对根轨迹加以描述说明。
A=[-21 19 -20;19 -21 20;40 -40 -40];B=[0;1;2];C=[1 0 2];D=0;[num,den]=ss2tf(A,B,C,D,1);rlocus(num,den);gridtitle('K增益负反馈闭环根轨迹图')图九（K增益负反馈闭环根轨迹图）采用K增益负反馈，画出如图所示的根轨迹图。
由图可知，共有3条根轨迹，第一条最终趋于原点；第二条收敛在20~60之间；第三条最终趋于无穷远处。
7、在Matlab中绘制Bode图和Nyquist图，并对图给予说明。
绘制Bode图：A=[-21 19 -20;19 -21 20;40 -40 -40];B=[0;1;2];C=[1 0 2];D=0;sys=tf(num,den) bode(num,den) gridtitle('Bode 图')汇出的波特图如图所示，由图可知，对复制响应分析可得，交越频率在转折频率之后，故复制的变化主要发生在低频段。
对相频特性进行分析，可知此系统的相频特性角度均为负值，并且最后的相角是趋于-90度的。
绘制Nyquist 图： nyquist(sys) title('Nyquist 图 ')图十（波特图）图十一（Nyquist 图）画出的奈奎斯特图如上所示，根据此图可知，此系统是稳定的系统，由奈奎斯特曲线可以分析出此系统的稳定性。
【实验结论与总结】经过本次试验我们了解了线性系统状态空间，知道如何在simulink 中用状态方程搭建系统仿真模型，以及调用模型。
不同输入对对系统有着一定的影响。
如何把方程转换成simulink 中的方块模型是非常重要的，在本次实验中，最困难的地方在于运用Simulink 生成图形，在这过程中由于要利用到sim 函数，不熟悉导致出错后来经过检查才改正错误。