第一章
状态方程转对角标准型
- 求特征方程 ∣λI−A∣=0 得出 λ
- 得出的解分别代入 λI−A 求解出 对应的特征向量p
- 线性非奇异变换矩阵 T=[p1,p2,...,pn]
- A=T−1AT,B=T−1B
无重根系统的标准型状态空间表达式
1. 能控标准型
状态空间表达式:
A=0 0 ⋮ 0 −a010⋮0−a101⋮0−a2⋯⋯⋱⋯⋯00⋮1−an−1
B=0 0 ⋮ 0 1
C=[b0b1⋯bn−1]
2. 能观标准型
状态空间表达式:
A=0 1 0 ⋮ 0001⋮0⋯⋯⋯⋱⋯000⋮1−a0−a1−a2⋮−an−1
B=b0 b1 b2 ⋮ bn−1
C=[00⋯01]
第二章
状态空间表达式 >> 传递函数
w(s)=C∗(sI−A)−1∗B+D
判断矩阵是否满足状态转移矩阵
Φ(0)=I
状态转移矩阵 >> 系统矩阵
Φ˙(t)∣t=0=A
已知 Φ,B,U 求输出响应 x(t)
x(t)=Φ(t−t0)x(t0)+∫t0tΦ(t−τ)BU(τ)dτ
初始时刻t0=0时
x(t)=L−1[(sI−A)−1]x(0)+L−1[(sI−A)−1Bu(s)]
第三章
线性定常连续系统能控性判据
秩判据
状态方程:x˙=Ax+Bu
有能控型判别矩阵Qc=[B,AB,A2B,...,An−1B] 满秩
PBH判据
所有特征值均满足
rank[λiI−A,B]=n
约当标准型判据
看书
线性定常连续系统输出能控性
能控性判别矩阵的秩等于输出向量的维数
rankQm=rank[CB,CAB,...,CAn−1B,D]=m
线性时变连续系统能控性判据
rank[M0(tf),M1(tf),...,Mn−1(tf)]=n
其中,
M0(tf) M1(tf) M2(tf) Mn−1(tf)=B(tf)=[−A(t)M0(t)+dtdM0(t)]t=tf=[−A(t)M1(t)+dtdM1(t)]t=tf=[−A(t)Mn−2(t)+dtdMn−2(t)]t=tf
例如,将t = 3s代入,发现成立。那么0 ~ 3s内都能控
线性定常连续系统能观测性判据
秩判据
系统方程:x˙=Ax+Bu, y=Cx+Du
有能观测性判别矩阵 Qo=C CA CA2 ⋮ CAn−1 满秩
PBH判据
所有特征值均满足
rank[λiI−A C]=n
约当标准型判据
看书
线性时变连续系统能观测性判据
系统在时刻 t0 能观测的充要条件是存在一个有限时刻 tf>t0,使得如下能观测性判别矩阵 满秩,即:
rankN0(tf) N1(tf) ⋮ Nn−1(tf)=n
其中,
N0(t) N1(t) N2(t) Nn−1(t)=C(t)=N0(t)A(t)+dtdN0(t)=N1(t)A(t)+dtdN1(t) ⋮=Nn−2(t)A(t)+dtdNn−2(t)
说明:判断时需将 t=t0 代入以上各矩阵。例如,在 t0=0 时判别矩阵满秩,则系统在 [0,t1] 区间内是能观测的。
第五章
闭环系统极点配置
被控系统∑o(A,B,C)的状态空间表达式为{x˙=[1 03−1]x+[0 1]u,求反馈增益矩阵F,使得闭环系统极点配置为−1+j,−1−j,并画出状态变量图
- 判断系统能控
- 确定闭环系统期望特征多项式
因为题给信息要求−1+j,−1−j,
有(s+i+j)(s+i−j)=s2+2s+2
所以a2∗=2,a1∗=2
- 求状态反馈增益阵F=[f1f2]
系统的特征多项式为(s−1)(s+1)=s2−1
有a2=1,a1=0
则能控标准型下的状态反馈增益阵为 F=[a2∗−a2a1∗−a1]=[32]
对于原式转成能控标准型的转换矩阵,有Tcc−1=([BAB⋯An−1B]an−1 ⋮ a1 1⋯1a11)
对于题目,可得Tcc−1=([3 −101][0 110])−1=31 3101
则状态反馈增益阵为 F=[f1f2]=FTcc−1=[32][31 3101]=[352]
- 画图
状态观测器
全维观测器
被控系统∑0(A,B,C)的状态空间表达式为⎩⎨⎧x˙=[1 03−1]x+[0 1]u y=[11]x,试设计全维观测器,使其极点为 −3,−3
- 判断系统能观性
- 确认闭环状态观测器系统矩阵的期望特征多项式
对于要求的极点为 −3,−3,有闭环特征多项式:
p∗(s)=∏i=1n(s−λi∗)=sn+a1∗sn−1+⋯+an−1∗s+an∗
对于题给信息,有 (s+3)(s+3)=s2+6s+9
有 a2∗=9,a1∗=6
- 求观测器偏差反馈增益矩阵G=[g1g2]T
先求当前系统的特征多项式:(s−1)(s+1),得a2=−1,a1=0
则能观标准型下的状态观测器增益矩阵为:Go=[a2∗−a2 a1∗−a1]=[10 6]
需要一个将原式转换成状态标准型的矩阵Toc=(an−1 ⋮ a1 1⋯1a11C CA ⋮ CAn−1)−1
对于原式,有Toc=([a1 110][C CA])−1=[10 6]
那么G=TocGo=[2 4]
- 列写观测器的状态方程
x^˙=Ax^+Bu−G(y^−y)
- 画图
降维观测器
设系统∑o(A,B,C)的状态空间表达式为⎩⎨⎧x˙=4 −7 10014−81x+1 0 −1u y=[101]x,试设计一个极点为 -4,-4的降维状态观测器
- 判断系统能观性
- 作线性变换
构造n∗n阶非奇异矩阵T=[Tn−m C]=1 0 1010001
T−1=1 0 −1010001
得 ⎩⎨⎧x˙=TAT−1+TBu=Ax+Bu=0 1 00014−85x+1 0 0u y=CT−1x=Cx=[001]x
Comments NOTHING