现代控制理论第一章 控制系统数学模型

uD
系统运动方程式为
KmiD
f
JD
d
dt
（式中， Ke 为电动势常数； Km 为转矩常数； J D 为折合到电动
机轴上的转动惯量； f 为折合到电动机轴上的粘性摩擦系数。）
可选择电枢电流 iD 和角速度 为状态变量，电动机的电 枢电压 uD为输入量，角速度 为输出量。
状态空间表达式 状态图如下：
返回到微分方程形式： z(n) an1z(n1) a1z a0 z u
以及
b z (n1) n1
b1z b0 z
y
选择状态变量如下：
x1 z
x1 x2 z
x2 x3 z
┆ xn1 xn z(n1)
xn z(n) a0 x1 a1x2 an1xn u
y bn1z(n1) b1z b0 z b0 x1 b1x2 bn1xn
x1
x
x2
xn
u1
u
u2
urym
a11
A
an1
a1n ann nn
c11
C
cm1
c1n cmn mn
b11
B
bn1
b1r anr nr
d11
D
dm1
d1r dmr mr
如果矩阵A, B, C, D中的所有元素都是实常数时，则称这样 的系统为线性定常（LTI，即：Linear Time-Invariant）系统。
其中
0 b3 0 1 b2 a20 0 2 b1 a10 a01 160 192 0 640 0 160 3 b0 a00 a11 a22 640 18160 2240
于是系统的状态空间表达式为
x1 0
x2
0
x3 640
1 0 192
0 x1 0
化简后，得
(M m)y ml u
my ml mg
求解得： y mg 1 u
MM
(M m)g 1 u
Ml
Ml
选择状态变量 x1 y ，x2 x1 y，x3 ，x4 x3
u 为系统输入， y 为系统输出
x1 0 1 0 0 x1 1
x2
0
0
mg M
0
x2
1
x2
160
u
18 x3 2240
x1
y 1
0
0
x2
x3
（2）辅助变量法 引入辅助变量z
z18z192z 640z u y 160z 640z
选择状态变量 x1 z x2 z x1 x3 z x2
于是系统的状态空间表达式为
x1 0
x2
0
x3 640
1 0 192
0 x1 0
1
x2
0u
18 x3 1
x1
y 640
160
0
x2
x3
1.3 传递函数矩阵
传递函数——系统初始松弛（即：初始条件为零）时，输出量的 拉氏变换式与输入量的拉氏变换式之比。
1.3.1 传递函数
单入-单出线性定常系统的状态空间表达式为
x Ax Bu
y Cx du
在初始松弛时，求Laplace变换，并且化简
diD dt
d
KRLmDD
dt J D
Ke LD f
JD
iD
1
LD 0
uD
y 0
1iD
例1-3 建立单极倒立摆系统的状态空间表达式。 单级倒立摆系统是许多重要的宇宙空间应用的一个简单模型。
在水平方向，应用牛顿第二定律：
M
d2 y dt2
m
d2 dt2
( y l sin ) u
x2
2
u
Ax
bu
x3 a0 a1 a2 x3 3
x1
y x1 0u 1
0
0
x2
0u
Cx
du
x3
系统的状态图
一般情况下，n 阶微分方程为： y(n) an1 y(n1) a1 y a0 y bnu(n) bn1u(n1) b1u b0u
选择 n 个状态变量为 系统方程为
y 1 0
x1
0
xn
1.2.2 微分方程中含有输入信号导数项
（一）待定系数法 首先考察三阶系统，其微分方程为
y a2 y a1 y a0 y b3u b2u b1u b0u
选择状态变量： x1 y 0u x2 y 0u 1u x1 1u x3 y 0u 1u 2u x2 2u
bn s n an s n
b1s b0 a1s a0
d
bn1sn1 b1s b0 ansn a1s a0
例1-4 已知描述系统的微分方程为 y18y192y 640y 160u 640u
试求系统的状态空间表达式。
解 （1）待定系数法
选择状态变量如下
x1 y 0u x2 x1 1u x3 x2 2u
┆ xn1 xn y(n1) xn y(n) a0 x1 a1x2 an1xn b0u
写成矩阵形式：
x1
x2
0 0
xn
0
a0
10 01 00 a1 a2
0 0 0 a3
0
0
0
x1 x2
0
u
1 an1
xn
0
b0
系统的状态图如下：
k m
y
f m
dy dt
1 m
F
k m
x1
f m
x2
1 m
F
机械系统的系统方程为
x1 x2
0
k m
1 f
m
x1 x2
0 1
m
F
y 1
0
x1 x2
该系统的状态图如下
例1-2 建立电枢控制直流他励电动机的状态空间表达式
电枢回路的电压方程为
LD
diD dt
RDiD
Ke
状态变量对输入量的传递函数
Gxu (s) sI
A 1b
adjsI detsI
A A
b
输出量对输入量的传递函数（即：传递函数）
g yu (s)
CsI
A 1b
d
C
adjsI detsI
i(t) 和 uC (t) 可以表征该电路系统的行为，就是该系统的一组状态
变量
1.1.2 状态空间表达式
前面电路的微分方程组可以改写如下，并且写成矩阵形式：
di(t) R i(t) uC (t) u(t)
dt L
LL
di(t)
dt duC (t)
1RL
dt C
1 L 0
状态空间——以所选择的一组状态变量为坐标轴而构成的正交线 性空间，称为状态空间。
例：如下图所示电路， u(t) 为输入量， uC (t) 为输出量。
建立方程：
L
di(t dt
)
Ri(t
)
uC
(t
)
u(t
)
i C duC (t) dt
初始条件：
i(t) t t0
i(t0 )
uC (t) tt0 uC (t0 )
x1 y 0u x2 x1 1u x3 x2 2u
xn xn1 n1u
x1
x2
0 0
xn
0
a0
10 01 00 a1 a2
0 0 0 a3
0 0
x1
x2
1
2
u
1 an1
xn
n1
n
y 1 0
系统状态图如下
x1
6、组合系统的数学描述 7、利用MATLAB进行模型之间的变换
1.1 状态空间表达式
1.1.1 状态、状态变量和状态空间 状态——动态系统的状态是一个可以确定该系统行为的信息集合。
状态变量——确定系统状态的最小一组变量，如果知道这些变量
在任意初始时刻 t0 的值以及 t ≥ t0 的系统输入，便能够完整地 确定系统在任意时刻 t 的状态。（状态变量的选择可以不同）
写成矩阵形式
x1
x2
0 0
xn
0
a0
10 01 00 a1 a2
0 0 0 a3
0
0
0
x1 x2
0 u
1 an1
xn
0 1
y b0
b1
x1
bn1
xn
注：如果输入项的导数阶次和输出项导数阶次相同，则有d。
Y (s) R(s)
系统的状态方程和输出方程一起，称为系统状态空间表达式，或称 为系统动态方程，或称系统方程。
设： x1 i(t) x2 uC (t)
C 0 1
x
x1
x2
A
R
L 1
-
1 L
0
C
x Ax Bu
则可以写成状态空间表达式：
y Cx
1
B
L 0
推广到一般形式：
x Ax Bu y Cx Du
2、微分方程中含有输入信号导数项，（即1.2.2 中的内容）
1.2.1 微分方程中不含有输入信号导数项
首先考察三阶系统，其微分方程为
y a2 y a1 y a0 y b0u
选取状态变量 x1 y
x2 y x3 y
则有 x1 x2 x2 x3 x3 a0 x1 a1x2 a2 x3 b0u
第1章 控制系统数学模型
本课程的任务是系统分析和系统设计。而不论是系统分析还是系 统设计，本课程所研究的内容是基于系统的数学模型来进行的。因 此，本章首先介绍控制系统的数学模型。