现代控制理论(第二章)讲解
浙大控制考研-现代控制理论(浙大)第二章

1 A2t 2 2!
1 k
Aktk
)
b0
t 0 x(0) b0
x(t) (I At 1 A2t 2 1 Akt k )x(0)
2!
k
eAt I At 1 A2t 2 1 Akt k
2!
k
矩阵指数函数
Φ(t) 状态转移矩阵
x(t) eAtx(0) 描述了状态向量由初始状态x(0)向任意时 刻状态 x(t)转移的内在特性。
eAt I At 1 A2t 2 1 Akt k
2!
k
1)根据状态转移矩阵的定义求解:
eAt I At 1 A2t 2 1 Akt k
2!
k!
对所有有限的t值来说,这个无穷级数都是收敛的 。
求出的解不是解析形式,适合于计算机求解。
例:求解系统状态方程 解:
x1
x2
0 0
-11
6
-6 -11 5
试计算状态转移矩阵 eAt .
解: 1) 特征值
1 1
I A 6 -11 6 1 2 3 0
6 11 5
1 1,2 2,3 3
2) 计算特征向量:
1 1 1 p1 0, p2 2, p3 6
1 4 9
3) 构造变换阵P:
1 1 1 P 0 2 6
(A B)3 A3 B3 3A2B 3AB 2
(9) x Px Φ(t) P-1Φ(t)P P-1eAtP
证明:非奇异线性变换
x Px
n n非奇异矩阵 另一组状态变量
x Px
x P1AP x x(t) eP1AP x(0)
x Ax APx 新的系统矩阵 新的状态转移矩阵
Ax
eAt x(0) Φ(t)x(0)
现代控制理论(II)-讲稿课件ppt

03
通过具体例子说明最小值原理在最优控制问题中的应
用方法。
06 现代控制理论应用案例
倒立摆系统稳定控制
倒立摆系统模型建立
分析倒立摆系统的物理特性,建立数学模型,包括运动方程和状态 空间表达式。
控制器设计
基于现代控制理论,设计状态反馈控制器,使倒立摆系统实现稳定 控制。
系统仿真与实验
利用MATLAB/Simulink等工具进行系统仿真,验证控制器的有效性; 搭建实际实验平台,进行实时控制实验。
最优控制方法分类
根据性能指标的类型和求解方法, 最优控制可分为线性二次型最优控 制、最小时间控制、最小能量控制 等。
最优控制应用举例
介绍最优控制在航空航天、机器人、 经济管理等领域的应用实例。
05 最优控制理论与方法
最优控制问题描述
控制系统的性能指标
定义控制系统的性能评价标准,如时间最短、能量最小等。
随着网络技术的发展,分布式控制系统逐渐 成为现代控制理论的研究热点,如多智能体 系统、协同控制等。
下一步学习建议
01
02
03
04
深入学习现代控制理论相关知 识,掌握更多先进的控制方法
和技术。
关注现代控制理论在实际系统 中的应用,了解不同领域控制
系统的设计和实现方法。
加强实践环节,通过仿真或实 验验证所学理论知识的正确性
机器人运动学建模
分析机器人的运动学特性, 建立机器人运动学模型, 描述机器人末端执行器的 位置和姿态。
运动规划算法设计
基于现代控制理论,设计 运动规划算法,生成机器 人从起始点到目标点的平 滑运动轨迹。
控制器设计与实现
设计机器人运动控制器, 实现机器人对规划轨迹的 精确跟踪;在实际机器人 平台上进行实验验证。
现代控制理论 2-0

∫
t
0
e − Aτ f (τ )dτ =
e [ x(0) + ∫ e
At 0 At
t
− Aτ
f (τ )dτ ] + ∫ e A( t −τ ) Bu (τ )dτ
t1 − Aτ
当t = t1时,有 x(t1 ) = e [ x(0) + ∫ e
0
f (τ )dτ ] + ∫ e A( t −τ ) Bu (τ )dτ
λ − 1 0 det[λI − A] = det = (λ − 1)(λ + 3) = 0 λ + 3 2 λ1 = 1, λ2 = −3 0 0 rank [λ1 I − AMb] = rank 2 4 − 4 rank [λ2 I − AMb] = rank 0 系统能控。 1 =2 1 0 1 =2 0 1
0
t1
∫
t1
0
e − Aτ f (τ )dτ为一个确定的值,仅仅相当于把系统
原来的初态改变了一确定的常值。所以在讨论系统 的能控性时,不考虑系统存在的确定性干扰。
第二章 系统的可观性和可控性
(三)能控性判据
判据一: 判据一:若系统能控,则能控性矩阵
Qc = [B AB A 2 B ... A n −1 B ] 满秩,即
第二章 系统的可观性和可控性
现代控制理论基础
主讲人: 主讲人:荣军 mail:rj1219 163. 1219@ E-mail:rj1219@
第二章 系统的可观性和可控性
2-1 能能控性及其判据
-、线性定常系统的能观测性及其判据 -、线性定常系统的能观测性及其判据
线性定常系统状态方程为 x = Ax + Bu 其中x、u分别为n、 r维向量,A、B为满足矩阵运算的常值矩阵。若给定系统的 一个初始状态x0和任一状态x1,如果在的有限时刻tf>0,定义在 时间区间[0,tf]的输入u(t)使状态x(0)=x0转移到x(tf)= x1 ,则称系统状态完全是能控的; 如果系统对任意一个初始状态都能控,则称系统是状态完全 能控的,简称系统是状态能控的或系统是能控的。
第二章 现在控制理论 线性系统的数学描述1

第二章现在控制理论线性系统的数学描述1第二章现在控制理论线性系统的数学描述1第二章线性系统的数学描述数学模型可以存有许多相同的形式,较常用的存有三种:第一种是:把系统的输入量和输出量之间的关系用数学方式表达出来,称之为输入输出描述,或外部描述;比如:微分方程式、传递函数和差分方程。
第二种是:不仅可以描述系统输入、输出之间的关系,而且还可以描述系统的内部特性,称之为状态空间描述或内部描述;它特别适用于于多输出、多输入系统,也适用于于时变系统、非线性系统和随机控制系统。
第三种是:用比较直观的方块图(结构图)和信号流图模型进行描述。
同一系统的数学模型可以表示为不同的形式,需要根据不同的情况对这些模型进行取舍,以利于对控制系统进行有效的分析。
许多表面上全然相同的系统(例如机械系统、电气系统、液压系统和经济学系统)有时却可能将具备完全相同的数学模型;从这个意义上讲,数学模型表达了这些系统的共性,所以只要研究透了一种数学模型,也就能完全了解具有这种数学模型形式的各式各样系统的本质特征。
92.1线性系统的时域数学模型对于单输出、单输入线性定常系统,使用以下微分方程去叙述:c(n)(t)?a1c(n?1)(t)?a2c?b0r(m)(n?2)?(t)?anc(t)(t)an?1c(m?1)(t)?b1r(t)?b2r(m?2)?(t)?bmr(t)(t)bm?1r(2.1)式中,r(t)和c(t)分别是系统的输入信号和输出信号,c(n)(t)为c(t)对时间t的n阶导数;ai(i?1,2,?n)和bj(j?0,1,?m)就是由系统的结构参数同意的系数。
通常情况下,列写控制系统运动方程的步骤就是(建模过程):首先,分析系统的工作原理及其各变量之间的关系,找出系统的输入量和输出量;其次,根据系统运动特性的基本定律,通常从系统的输出端的已经开始依次写下各元件的运动方程,在列写元件运动方程时,须要考量相连元件间的相互作用;最后,由组成系统各元件的运动方程中,消去中间变量,求取只含有系统输入和输出变量及其各阶导数的方程,并将其化为标准形式。
第2章 现代控制理论1PPT课件

时不变系统状态转移矩阵Φ tt0或 Φ t是满足如下矩阵微分
方程和初始条件的解,这也是检验一个矩阵是不是状态转移
的条件。
Φ (tt0)AΦ (tt0)或 Φ (t)AΦ (t)
Φቤተ መጻሕፍቲ ባይዱ(0)I
Φ (0)I
(2.5)
1Φ t在 t0的值 lim ΦtI
t0
(2)Φt对t的导 Φ 数 tA Φ tΦ tA
故可求出其解为:
t
X ( t) ( t) X ( 0 ) o ( t ) B () U d ( 2 .2 b )
式中 (t) eAt 为系统的状态转移矩阵。
对于线性时变系统非齐次状态方程,
X ( t) A ( t) X ( t) B ( t) U ( t) ( 2 3 )
类似可求出其解为
x (0 )e a t tb(u )e a (t )d 0
同样,将方程(2.1)写为 X (t)A(X t)B(U t)
在上式两边左乘eAt ,可得:
e A [X t(t) A(t) X ]d[e AX t(t) ]e A B t (tU )
dt
3
将上式由 0 积分到 t ,得
X ( t) e A X t ( 0 ) te A (t )B () U d (2 .2 a ) o
的解,X(t)=Ф (t, t0)X(0) 。 下面不加证明地给出线性时变系统状态转移矩阵的几个
重要性质: 1、 (t,t)I
2 、 ( t 2 ,t 1 ) ( t 1 ,t 0 ) ( t 2 ,t 0 )
3 、 1 (t,t0) (t0 ,t) 4、当A给定后,(t,t0) 唯一
5、计算时变系统状态转移矩阵的公式
令 x (t) b 0 b 1 t b 2 t2 b iti b iti,t 0
第二章现代控制理论状态空间表达式

即
(2-11)
(3) 列出状态空间描述iL 1 − ( R + R )C 1 2 R1 L( R1 + R2 ) − R1 1 ( R1 + R2 )C uC ( R1 + R2 )C (2-12) + e(t ) R1 R2 iL R2 − L( R + R ) L( R1 + R2 ) 1 2
§2.1 状态空间描述的概念 2.1.2 控制系统的状态空间描述举例
例2-1 R-L-C系统,求其状态空间描述
R
u
L i
C
uC
解 (1) 确定状态变量 选择电容两端电压 uC (t )、电感通过的电流 i (t ) (2) 列写微分方程并化为一阶微分方程组 基尔霍夫(Kirchhoff)电压定律,
(2-13)
令
1 − ( R + R )C 1 2 A= R1 L( R + R ) 1 2
1 ( R + R )C 2 b= 1 R2 L( R + R ) 1 2
−
R1 ( R1 + R2 )C R1 R2 − L( R1 + R2 )
n 维列向量,状态向量
a12 a1n a22 a2 n an 2 ann
n×n方阵,系统矩阵(或状态矩阵), 反映系统状态的内在联系
§2.1 状态空间描述的概念
现代控制理论-第二章 控制系统的状态空间描述

DgXu
2.2.1.由物理机理直接建立状态空间表达式: 例2.2.1 系统如图所示
L
R2
u
iL
R1
uc
选择状态变量:
x1 iL , x2 uC ,
13 中南大C diL 1 iL (u L ) C dt R1 dt duC diL L uC C R2 u dt dt
y(s) [C(sI A) B D]U (s)
1
1
得
9
G(s) C (sI A) B D
命题得证
中南大学信息学院自动化系
1
DgXu
例2.1.3
已知系统的状态空间描述为
x1 0 1 0 x1 0 x 0 1 1 x 1 u 2 2 x3 0 0 3 x3 1
28 中南大学信息学院自动化系
DgXu
故有(n-1) 个状态方程:
对xl求导数且考虑式 (2.3.12),经整理有:
则式 (2.3.12) bn=0 时的动态方程为:
(2.3.16)
式中:
29 中南大学信息学院自动化系
DgXu
30 中南大学信息学院自动化系
DgXu
3)
化输入-输出描述为状态空间描述
11 中南大学信息学院自动化系
DgXu
2.3. 线性定常连续系统状态空间表达式的建立
建立状态空间表达式的方法主要有两种: 一是直接根据系统的机理建立相应的微分方程或差分方 程,继而选择有关的物理量作为状态变量,从而导出其状态 空间表达式; 二是由已知的系统其它数学模型经过转化而得到状态达 式。由于微分方程和传递函数是描述线性定常连续系统常用 的数学模型,故我们将介绍已知 n 阶系统微分方程或传递函 数时导出状态空间表达式的一般方法,以便建立统一的研究 理论,揭示系统内部固有的重要结构特性。
现代控制理论第二章

1 s 1 1 s
Y(s)
2.5 由系统方块图导出状态空间描述
方块图导出状态空间模型的步骤 (1)将系统方块图中的每一环节都分解为积分环节和 惯性环节的组合。 (2)以所有惯性环节和积分环节的输出作为状态变量 的拉氏变换。 (3)列出所有惯性环节和积分环节输入输出的拉氏变 换关系式。 (4)对所有(3)中的拉氏变换关系式求拉氏反变换 得到一阶微分方程组。 (5)把(4)中的一阶微分方程组化成向量矩阵表示 的状态方程与输出方程。
例2.7
2.5 由系统方块图导出状态空间描述
(一)方块图方法的思路 当系统的描述以方块图形式给出时,常常无须求 出系统的总传递函数和状态变量图,可以直接由方块 图导出其相应的状态空间模型. 方法:把系统中二阶以上的环节化为由惯性环节 和积分环节组成。
2.5 由系统方块图导出状态空间描述
(二)典型二阶系统状态空间描述
2.6
习题课
2.7 系统状态空间描述与传递函数
设线性连续定常系统的状态空间模型为
x Ax Bu y C x Du
对以上两式分别做拉氏变换,得
sX s AX s BU s
Y s CX s DU s
从以上两式中消去 X s , 则
s a1 s
n 1
k1 s s1
k2 s s2
kn s sn
k i lim W s s s i
例2.5
s si
2.3 系统的频域描述化为状态空间描述
(二)控制系统传递函数的极点为重根 (1) 传递函数的极点为一个重根
W s Y s U s k 11 k 12 k 1n s s1
现代控制理论基础_周军_第二章状态空间分析法

2.1 状态空间描述的基本概念系统一般可用常微分方程在时域内描述,对复杂系统要求解高阶微分方程,这是相当困难的。
经典控制理论中采用拉氏变换法在复频域内描述系统,得到联系输入-输出关系的传递函数,基于传递函数设计单输入-单输出系统极为有效,可从传递函数的零点、极点分布得出系统定性特性,并已建立起一整套图解分析设计法,至今仍得到广泛成功地应用。
但传递函数对系统是一种外部描述,它不能描述处于系统内部的运动变量;且忽略了初始条件。
因此传递函数不能包含系统的所有信息。
由于六十年代以来,控制工程向复杂化、高性能方向发展,所需利用的信息不局限于输入量、输出量、误差等,还需要利用系统内部的状态变化规律,加之利用数字计算机技术进行分析设计及实时控制,因而可能处理复杂的时变、非线性、多输入-多输出系统的问题,但传递函数法在这新领域的应用受到很大限制。
于是需要用新的对系统内部进行描述的新方法-状态空间分析法。
第一节基本概念状态变量指描述系统运动的一组独立(数目最少的)变量。
一个用阶微分方程描述含有个独立变量的系统,当求得个独立变量随时间变化的规律时,系统状态可完全确定。
若变量数目多于,必有变量不独立;若少于,又不足以描述系统状态。
因此,当系统能用最少的个变量完全确定系统状态时,则称这个变量为系统的状态变量。
选取状态变量应满足以下条件:给定时刻的初始值,以及的输入值,可唯一确定系统将来的状态。
而时刻的状态表示时刻以前的系统运动的历史总结,故状态变量是对系统过去、现在和将来行为的描述。
状态变量的选取具有非唯一性,即可用某一组、也可用另一组数目最少的变量。
状态变量不一定要象系统输出量那样,在物理上是可测量或可观察的量,但在实用上毕竟还是选择容易测量的一些量,以便满足实现状态反馈、改善系统性能的需要。
状态向量把描述系统状态的个状态变量看作向量的分量,则称为状态向量,记以,上标为矩阵转置记号。
若状态向量由个分量组成,则称维状态向量。
现代控制理论(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
3
3
100
(sI
A) 1
1 s3
0 1 0 s2
001
010
0 0 1s 000
001 000 000
11 1 s s2 s3 01 1
s s2 00 1
s
11 1
次(u(t) 0 )状态方程的解
x(t) e At x0
At k 0 kk ! k x0
定义矩阵指数: e At
Akt k k 0 k!
I At 1 A2t 2 2
阵。
1
A t ,它仍是一个矩
kk
k!
若初始时间为t 0,则状态方程的解为
x(t ) e x A(t t0 ) 0
Ak (t t0 )k x0
幂零矩阵:存在某一正整数 k ,使得 A 0 称为 k 次“幂零矩阵”。 A 为幂零矩 阵的“充要条件”是 A 的所有特征值为k零: AX A , i 0 i 1,2, , n
特例: A 为数字矩阵,即 A
P31 例 2-1: A
010 0 0 1 , A2 000
001 0 0 0 , A3 000
(2) e0 I ; (3) e A]t 称e A为频e域A(t求)法;或(叫4)Lapl(aecAet )变1 换法e A;t
(5) 若矩阵 A、B 满足交换律 AB BA,则有e At e Bt ( A B)t (;A 、 B 可交换的 e
充要条件是 AB 为反称矩阵, A A 称为对称矩阵, A A 称为反称矩阵)
(9) 传递性:对任意满足t t t ,有e A(t2 t ) e A(t1 t )
1
0
2
1
0
e A(t2 t0 ) 。这表明状态
现代控制理论(第二章)线性系统的状态空间描述

H[t0 ,)
yc
1
yc
u
t t0 0
容易得到其解
yc
(t )
e
1t
yc
(0)
t
e1
(t
)u(
)d
显然,若其初始条件
yc
0
(0)
不能确定,则不能
唯一地确定其输出。
1.非零初始条件与脉冲输入
零初始条件:系统的初始条件为零是指系统在初 始时刻没有能量储备。
注意:在建立线性系统的输入—输出描述时, 必须假设系统的初始条件为零。
单变量线性时变系统输入-输出关系: y L(u)
用符号 g(t,τ) 表示该系统的单位脉冲响应,即
g(t,τ)L( (t ))
注意: g(t,τ) 是双变量函数; τ— 代表δ函数作用于系统的时刻; t — 代表观测其输出响应的时刻。
结论1:对单变量线性时变系统,u(t)为其输 入变量,g(t,τ)为其单位脉冲响应,在初始
y
kp
u
s3 1s 2 2s 3
若对其参数一无所知,它的控制律设计就会复 杂得多,而稳定性的分析事实上是无法进行的。
系统的输入—输出描述仅在松弛的条件下才能采用。
若系统在t0时刻是非松弛的,输出 y[t0 ,) 并不能单
单由 u[t0 ,) 所决定,即关系式 不成立。考察简单的一阶系统:
y[t0 ,)
初始条件不为零时,可以将非零的初始条件等 效成在初始时刻的一个脉冲输入。
单位脉冲函数(δ函数 )
令
0
(t
t1
)
1
0
t t1 t1 t t1 t t1
当Δ→0时, (t t1) 的极限函数,即
现代控制理论基础 第2章 控制系统的状态空间描述

【例3】建立图2-1所示RLC电路的状态方程。
取电容上的电压uC (t)和电感中的电流i(t)作为状态变量, 根据电路原理有
C duc (t) i(t) dt
di(t) L dt Ri(t) uc (t) u(t)
将上式中状态变量的一阶导数放在方程左边,其余项 移至方程右边,整理得一阶微分方程组为
状态空间法具备如下优点: (1)在数字计算机上求解一阶微分方程组或者差分方程
组,比求解与它相当的高阶微分方程或差分方程要容易。
(2)状态空间法引入了向量矩阵,大大简化了一阶微分方 程组的数学表示法。
(3)在控制系统的分析中,系统的初始条件对经典法感 到困难的问题,采用状态空间法就迎刃而解了。
(4)状态空间法能同时给出系统的全部独立变量的响应, 不但反映了系统的输入输出外部特性,而且揭示了系统 内部的结构特性,既适用单输入单输出系统又适用多输 入多输出系统。
x = A(t)x B(t)u
y
C (t ) x
D(t)u
式中,各个系数矩阵分别为
(2-8)
a11 (t)
A(t)
an1 (t)
c11 (t)
C
(t)
cm1 (t)
a1n (t)
b11 (t)
,
B(t)
ann (t)
bn1 (t)
c1n (t)
d11 (t)
,
D(t)
cmn (t)
述把系统的输出取为系统外部输入的直接响应, 显然这种描述回避了表征系统内部的动态过程 即把系统当成一个“黑匣”,认为系统的内部 结构和内部信息全然不知,系统描述直接反映 了输出变量与输入变量间的动态因果关系。
考察图2-1所示的n级RC网络。图中虚线框内 为具有放大器隔离的n级RC电路,设放大器的输入阻
现代控制理论_2-1_线性系统的状态空间描述

第二章 线性系统的状态空间分析法c§1 线性系统的状态空间描述§2 线性定常连续系统的分析e a第二章 线性系统的状态空间分析法§3 线性定常离散系统的分析 §4 系统的传递函数矩阵ty cc§1 线性系统的状态空间描述§2 线性定常连续系统的分析e a§3 线性定常离散系统的分析 §4 系统的传递函数矩阵ty c叠加 原理一、系统描述的基本概念 一、系统描述的基本概念1,输入、输出2,松弛性:若系统的输出y[t0,∞) 由输入u[t0,∞)唯一确定,则称系统在t0时刻是松弛的。
系统在t0 时刻不存储能量,初始条件为零!ce a4,线性: H (u1 + u 2 ) = Hu1 + Hu2 可加性5,定常性: Qa 为位移算子y = Hu算子,如传递函数3,因果性:系统在t 时刻的输出仅取决于t时刻和t 时刻之前的输入,与t 时刻之后的输入无关。
ty cy1 y2 yq输入延迟输出相应延迟 y = Hu = HQa u = Qa Hu = Qa y = y (t − a )u (t ) y (t ) u (t ) y (t )cH (αu1 ) = αH (u1 )e at齐次性a为实数u (t ) = Qa u (t ) = u (t − a )ty ct系统 数学描述cu1 u2 up外部描述(输入—输出描述)不完全描述微分方程、传递函数e aM二、状态空间的基本概念 二、状态空间的基本概念Mx 例:机械位移系统 依据牛顿定律: ∑ F = m&&(t )系 统x1 , x2 , L , xn内部描述(状态空间描述)状态方程+输出方程ty c完全描述ckmF (t )e ax (t )& ∴ F (t ) − kx(t ) − fx(t ) = m&&(t ) x微分方程:& m&&(t ) + fx(t ) + kx(t ) = F (t ) xf传递函数:X (s ) 1 = (s ) ms 2 + fs + k F经典控制理论中的数学模型(外部描述),反映 了输入输出的关系,不能反映内部变量的关系。
现代控制理论课件chapter2

Modern Control Theory
L06
Chapter 2 State Space Analysis of Control Systems
2.1 线性定常系统齐次状态方程的解
矩阵指数法
(t ) Ax(t ), 对应于 t 的同次幂系数相等 x
1 1 2 1 2 b1 Ab0 ,2b2 Ab1 b2 Ab1 A b0 A b0 2 2 2! 1 1 1 3 3 3b3 Ab2 b3 Ab2 A b0 A b0, , b0 x(0) 3 3 2! 3! 1 k bk A b0 k!
Modern Control Theory
L06
Chapter 2 State Space Analysis of Control Systems
2.1 线性定常系统齐次状态方程的解
矩阵指数法
1 2 2 1 k k 所以 x(t ) [ I At A t A t ]x(0) 2! k! 1 2 2 1 k k 因为 e 1 at a t a t 2! k!
2.2状态转移矩阵
转移矩阵的计算
Modern Control Theory
2.1 线性定常系统齐次状态 方程的解 (自由解)
L06
Chapter 2 State Space Analysis of Control Systems
2.1 线性定常系统齐次状态方程的解
定义:自由解 矩阵指数法 拉氏变换法
为一阶齐次微分方程组。
自由解:系统在没有输入的情况下,由初始状 态引起的自由运动。
Modern Control Theory
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
sI
A 1
s 2
s3
1 1 s 3
(s
1)(s 2
2)
(s 1)(s 2)
1
(s
1)(s s
2)
(s 1)(s 2)
s3
e At
L1
(s
1)( s 2
2)
(s 1)(s 2)
EAST CHINA INSTITUTE OF TECHNOLOgy
第二章 控制系统状态空间表达式的解
2.1 线性定常齐次状态方程的解(自由解) 2.2 矩阵指数函数——状态转移矩阵 2.3 线性定常系统非齐次方程的解 2.4 * 线性时变系统的解 2.5 * 离散时间系统状态方程的解 2.6* 连续时间状态空间表达式的离散化
(s
1)( s 2
2)
(s 1)(s 2)
1
(s
1)( s s
2)
(s 1)(s 2)
eAt L1
sI A 1
2et e2t 2et 2e2t
et e2t
et
2e2t
et
2e2t
例2-6,利用凯莱-哈密顿定理— -----------------自学! 例2-3与例2-7也请注意自学!
EAST CHINA INSTITUTE OF TECHNOLOgy
2.3 线性定常系统非齐次方程的解
现在讨论线性定常系统在控制作用 方程为非齐次矩阵微分方程:
作用下的强制运动。此时状态
1 1 1 1
P 1
2
1
2
P1 12
1
1
2 1
1 1
e At
Pe At P1
e1t
P
0
0 e2t
P
1
1 11 et Nhomakorabea2
0
et
e2t 2 1 2et e2t
证明: 级数形式
和标量微分方程求解类似,先假设式(1)的解
(3) 为 的矢量幂
EAST CHINA INSTITUTE OF TECHNOLOgy
(4) 代入式(1)得:
(5)
既然式(4)是式(1)的解,则式(5)对任意时刻 都成立,故 的同次 幂项的系数应相等,有:
EAST CHINA INSTITUTE OF TECHNOLOgy
EAST CHINA INSTITUTE OF TECHNOLOgy
在特定控制作用下,如脉冲函数、阶跃函数和斜坡函数的激励下,则 系统的解式(2)可以简化为以下公式:
1.脉冲响应
即当
时
2.阶跃响应
即当
时
3.斜坡响应
即当
时
(6) (7)
(8)
例2-8 要求掌握!
EAST CHINA INSTITUTE OF TECHNOLOgy
sI A 1
2et e2t 2et e2t
et e2t
et
2e2t
x(t) L1 (sI A)1 x(0) L1 (sI A)1BU(s)
EAST CHINA INSTITUTE OF TECHNOLOgy
在式(4)中,令
,可得:
将以上结果代入式(4),故得: (6)
EAST CHINA INSTITUTE OF TECHNOLOgy
等式右边括号内的展开式是 即
于是式(6)可表示为:
矩阵,它是一个矩阵指数函数,记为 , (7)
再用 的正确性。
代替
即在代替 的情况下,同样可以证明式2)
EAST CHINA INSTITUTE OF TECHNOLOgy
(0)
et et
1 e2t 2 e2t
1
2
EAST CHINA INSTITUTE OF TECHNOLOgy
例2-8:已知系统状态方程中
0 1 0
试求解该系统的单位阶跃响应。A 2 3,b 1
解法二:拉氏变换法
或
(1)
这就是组合性质,它意味着从 转移到0,再从0转移到 的组合。
2.性质二
或
(2)
注:本性质可用于判断矩阵是否符合状态转移矩阵的条件 3.性质三
或
(3)
EAST CHINA INSTITUTE OF TECHNOLOgy
4.性质四
或
(4)
这个性质说明,
矩阵与A矩阵是可以交换的。
注:本性质还表明,由状态转移矩阵
EAST CHINA INSTITUTE OF TECHNOLOgy
2.1 线性定常齐次状态方程的解(自由解)
所谓系统的自由解,是指系统输入为零时,由初始状态引起的自由 运动。此时,状态方程为齐次微分方程:
(1)
若初始时刻 时的状态给定为
则式(1)有唯一确定解:
若初始时刻从
开始,即
(2) 则其解为:
x(t) Φ(t)x(0) x(t) Φ(t t0 )x(t0 )
线性定常系统的状态转移矩阵
注:状态矩阵一般不是常数,而是时间的函数 起始矢量可以任意取,系统求解区间可任意选定—状态空间法的优点
EAST CHINA INSTITUTE OF TECHNOLOgy
2.2.2 状态转移矩阵(矩阵指数函数)的基本性质 1.性质一
可反推A!
5.性质五
对于
方阵A和B,当且仅当AB=BA时,有
而当AB≠BA是,则
这个性质说明,除非距阵A与B是可交换的,它们各目的矩阵指数函 数之积与其和的矩阵指数函数不等价。这与标量指数函数的性质是不同的。
EAST CHINA INSTITUTE OF TECHNOLOgy
2.2.3 几个特殊的矩阵指数函数 1.若 A 为对角线矩阵,即
对式(4)在
上间积分,有:
整理后可得式(2):
同理,若对式(4)在
上积分,即可证明式(3)。
式(2)也可从拉氏变换法求得,对式(1)进行拉氏变换,有:
即
EAST CHINA INSTITUTE OF TECHNOLOgy
上式左乘
,得:
注意式(5)等式右边第二项,其中:
(5)
两个拉氏变换函数的积是一个卷积的拉氏变换,即 以此代入式(5),并取拉氏反变换,即得 :
et
2e2t
1
1 2et 2e2t
0 2 1
e
2t
1
1
et e2t
et
2e2t
3)用拉氏变换法求解 e At L1 (sI A)1
EAST CHINA INSTITUTE OF TECHNOLOgy
(s (s
1 1)( s
s 1)( s
2) 2)
L1
2
s 1 2
s 1
1
s2 2
s2
1
s 1 1
1
s2 2
s 1 s 2
2et e2t
et e2t
2et 2e2t
2et 2et
e2t 2e2t
et e2t
et
2e2t
u(t) 1
t Φ( )Bu(t )d 0
t 2e e2 0 2e 2e2
e e
e2 2e2
10d
2.2 矩阵指数函数——状态转移矩阵
2.2.1 状态转移矩阵 齐次微分方程(1)的自由解为:
或
令Φ(t) eAt, 反应了由初始状态到时 间t的运动规律
该式反应了状态矢量由初始状态到任意时刻的矢量变换关系,反应了 状态矢量在空间随时间转移的规律,因此称为状态转移矩阵。
EAST CHINA INSTITUTE OF TECHNOLOgy
求解,记得公式(13)。
EAST CHINA INSTITUTE OF TECHNOLOgy
例2-1,2-2,2-4:求以下矩阵A的状态转移矩阵
[解]: 1)直接算法(略)
2)用标准型法求解
A
0 2
1 3
特征值: 1 1, 2 2
特征值互异 ,转化成对角标准型,且A为友矩阵
2.2.4 1.根据
的计算 的定义直接计算
编程,用计算机算,最终能得到收敛解。但很难得到解析解。例2-1 2.变换 A 为约旦标准型 (1)A 特征根互异
其中 T 是使 A 变换为对角线矩阵的变换阵。由式(7),有:
EAST CHINA INSTITUTE OF TECHNOLOgy
3.利用拉氏反变换法求 (10)
证明 齐次微分方程
两边取拉氏变换
即 故
EAST CHINA INSTITUTE OF TECHNOLOgy
对上式两边取拉氏反变换,从而得到齐次微分方程的解:
4.应用凯莱—哈密顿定理求 (1)由凯莱—哈密顿定理,方阵A满足其自身的特征方程,即
所以有
它是 同理
的线性组合。
EAST CHINA INSTITUTE OF TECHNOLOgy