现代控制理论Chapter4

Controllability and Observability

能控性与能观测性

在有限时间内,控制作用能否使系统从初始状态转移到要求的状态?
• 指控制作用对状态变量的影响或控制能力,称之为状态能控性问题决定能否实现最优控制

 在有限时间内,能否通过对系统输出的测定来估计系统的各个状态?
• 指系统的输出量对系统状态的识别能力,称之为状态能观性问题决定能否实现状态反馈控制

Controllability 能控性

如果对于任意初始状态x(0) = x0和任意目标状态x1,存在一个输入,能在有限时间内将x0驱动到x1,就说系统(A, B)称为能控的

在Ch3中,我们已经知道系统的状态响应方程:

为了满足定义,我们希望在输入项的作用下,能够抵消掉初始状态的值,并且产生一个目标状态的值,为此,我们构造如下

我们希望是一个方阵,能够方便的求逆,所以对其乘以
因此可以令输入项为:
抵消掉 矩阵,产生一个单位阵。
代入原式之后可以得到:

因此得到以下定理:

能控性的判据1

(A, B) is controllable if and only if the following matrix

is for every .

也即,只有当 矩阵非奇异时,可以对它求逆矩阵,这个时候系统才稳定。说明所有状态方向都可以被控制输入影响到。
如果奇异,说明有某些状态方向无论如何施加控制输入,都无法到达,也就是说系统不可控。

如果矩阵非奇异,即可逆时,根据我们的构造过程,显然系统能控,充分性得证。

只要有某一个时刻 非奇异,那么对于任意 的时刻该矩阵都非奇异。

必要性:
即证明系统能控时, 矩阵必然非奇异(可逆)

采用 反证法 的思想:
假设 不可逆,这意味着W矩阵有零特征值和特征向量,所以

which implies

加入系统此时仍然可控,我们可以把一个非零的状态 转化为 .
因此

.
然而作为特征向量
出现了矛盾,这说明必须是非奇异的,必要性得证。

若系统能控,在控制律
的作用下,可在任意短时间内将非零状态转移到零状态,我们把下面的矩阵叫做格拉姆矩阵

能控性的判据2

是可控的,当且仅当如下的 矩阵
的秩为 (即满秩)。

该判据基于对 的泰勒展开得到,当 矩阵满秩时,可以推出 矩阵奇异,故而系统能控。
证明参考:过程证明 > 能观性判据2 (rank条件)

能控性的判据3

当且仅当对于 的任意特征值 ,有
时, 是可控的。

对于证明过程中提到的控制输入:

能够将 转换到
这是一个开环控制器,非常的敏感,

若系统可控,必定可以得出系统稳定。

Observability 能观性

能观性表示的是根据系统输出 和输入 反映状态矢量 的能力。

PixPin_2025-04-21_10-21-08.png
如果对于任何未知的初始状态 ,存在一个有限的 ,从而可以根据 上精确地确定 ,则称该系统为能观系统。

由于输入 是可以任意给定的,可以令 ,从而能观性只依赖于矩阵 和矩阵

前的系数构造为方阵

进行积分

这与可控性的形式是很像的,所以得到
其中
能观的条件:

能观性的判据1

非奇异

同样,我们期待一种更简单的方法,将 进行泰勒展开,表示为
的线性组合。

可以得到

能观性的判据2

是可观测的,当且仅当以下矩阵 的秩为 n(满列秩)。

等价于

能观性的判据3

对于的任意特征值 ,有

基于能观性,可以得到一个 的估计器
这是一个开环估计器,在理想情况下估计是准确的。

能控能观具有鲜明的对偶性
对于能控

  • x:系统状态
  • u:输入(控制作用)
  • y:输出(观测量)
  • A, B, C:系统矩阵
    对于能观
  • z:对偶系统的状态
  • v:对偶系统的输入(对应原系统的“观测”)
  • w:对偶系统的输出(对应原系统的“控制”)
判据类别 能控性 (Controllability) 等价条件 能观性 (Observability) 等价条件
基本定义 系统对 是状态可控的。 系统对 是状态能观的。
秩条件 (Kalman) 可控性矩阵 满列秩,即 能观性矩阵 满列秩,即
PBH 判据 (Hautus) 的任意特征值 ,矩阵 满行秩,即 的任意特征值 ,矩阵 满列秩,即
格拉姆矩阵判据 (连续时间) 可控性格拉姆矩阵 对任意 是满秩(正定)的。 能观性格拉姆矩阵 对任意 是满秩(正定)的。
格拉姆矩阵判据 (离散时间) (离散) 可控性格拉姆矩阵 是满秩(非奇异)的。 (离散) 能观性格拉姆矩阵 是满秩(非奇异)的。

Observability and Controllability of Equivalent Systems 等价系统的能观能控性

定理: 考虑线性时不变系统 。若存在一个非奇异矩阵 (即 满秩,可逆),进行状态变换 ,得到等价系统 ,其中 。那么,原系统 是状态可控的,当且仅当变换后的系统 是状态可控的。

只要 非奇异,就有

能观性的证明同理。

等价变换并不改变系统的能控性/能观性。(不改变传递函数,不改变极点,不改变稳定性)

对于任意能控/能观系统,都可以通过线性变换 ,转化为
系统被分解为Jordan Block Systems

对于对角形式
PixPin_2025-04-21_15-12-53.png
B矩阵的每一行都不能全零,确保输入通道能够影响该状态变量。
C矩阵的每一列都不能全零,使得每个状态分量都能在输出中留下痕迹。

对于Jordan Form
,有一个若当块,该若当块的输入尾行和输出首列需不全0
PixPin_2025-04-21_15-41-29.png

,对于如图所示的情况,有两个若当块,每一个的输入尾行和输出首行都不全0,并且要求这些行之间和这些列之间线性无关
PixPin_2025-04-21_15-30-25 1.png

一般定理 (General Theorem)

  • 不同 特征值
  • 对每个 ,在 Jordan 标准型里有 个 Jordan 块:
  • 定义
    • 中对应 最后一行的行向量;
    • 中对应 第一列的列向量。

能控 ⇔ 对每个 ,集合

在输入空间中线性无关

能观 ⇔ 对每个 ,集合

在输出空间中线性无关

直观理解
输入只能立刻作用在“链尾”,输出只能最先“看见”链首,所以对链首和链尾有对应要求。

对于一个SISO系统:

  • 可控 当且仅当对于每个不同的特征值,只有一个Jordan块,并且B中对应Jordan块最后一行的每个元素都不为零。
  • 可观 当且仅当对于每个不同的特征值,只有一个Jordan块,并且C中对应Jordan块第一行的每个元素都不为零。
    若单个特征值对应多个Jordan块,则不可能满足线性无关这一约束,这时系统无法能观/能控。

Decomposition 分解

Controllability Decomposition 能控性分解

对于任意系统 ,总存在一个等价变换 ,使得变换后的系统 具有如下结构:

这样一个变换将系统分解为能控部分和不能控的部分。
PixPin_2025-04-29_00-52-28.png

其中:

  • 子系统 是完全能控的,其维数为
  • , ,
  • 对应系统的不能控部分。

矩阵可以看出,输入 无法影响 对应的状态,故而不能控。

我们已经知道,等价变换不改变系统的传递函数。当我们对系统做能控性分解后,通过简单的矩阵运算不难发现系统的传递函数仅由能控部分决定

系统状态空间被划分为能控子空间和不能控子空间。

如何寻找的变换矩阵

  1. 写出原系统的能控性矩阵

这里 就是可控子空间的维数。

  1. 中选出一组最大线性无关列

它们张成可控子空间

  1. 补齐基底
    再任取 个向量 ,作为不可控子空间的方向

使得

构成 的一组基底。

  1. 构造变换矩阵

可逆。定义

代入可得

即完成了能控性分解。

注意这里的顺序不能搞反了

Observability Decomposition 能观性分解

能观性分解可以看作能控性分解的对偶(Duality)

类似地,对于任意系统 ,总存在一个等价变换 ,使得变换后的系统 具有:

这样一个变换将系统分解为能观部分和不能观的部分。
PixPin_2025-04-29_01-23-40.png

其中:

  • 子系统 是完全能观的,其维数为
  • , ,
  • 对应系统的不能观部分。

从矩阵不难看出,状态 (对应 的状态) 并不会影响输出 ,故而是不能观的。

系统的传递函数仅由能观部分决定,证明类似过程证明 > 系统的传递函数仅由能控部分决定,不再赘述。

系统状态空间被划分为能观子空间和不能观子空间。

注意,能控性分解里补全后的矩阵 ,而能控性补全后的矩阵直接就是P

Kalman Decomposition 卡尔曼分解

先对一个系统进行能控性分解,再进行能观性分解,就可以得到卡尔曼分解。

任一线性时不变系统 都可以通过等价变换 分解为如下形式:

其中:

  • : 能控且能观 () 的状态子空间
  • : 能控但不能观 () 的状态子空间
  • : 不能控但能观 () 的状态子空间
  • : 不能控且不能观 () 的状态子空间

PixPin_2025-04-29_20-39-54.png

并且可以得出,系统的传递函数仅由能控且能观的部分决定

只有能控且能观的部分才对系统的传递函数(即输入输出特性)有贡献。

Minimal Realization 最小实现

状态方程表达形式的多样性与传递函数的唯一性反映了系统描述中的冗余可能性。通过能控性与能观性分解的卡尔曼分解过程,我们实际上是在识别并消除这种冗余,因为分解后得到的矩阵维度通常会减小,而系统的输入输出特性(即传递函数)保持不变。这一过程本质上是寻找系统的最小实现,即使用最小维度的状态空间来完整描述系统动态特性,从而在保持系统传递函数不变的前提下,获得最为精简且无冗余的状态空间表达。

对于一个给定的传递函数(或传递矩阵) ,可以找到多个状态空间表示 使得:

这些不同的状态空间表示称为 实现 (Realization)

定义:
如果一个实现 状态空间的维数 (即矩阵 的阶数)在所有可能的实现中是最小的,则称该实现为 最小实现 (Minimal Realization)

如果一个系统(A,B,C,D) 是最小实现,当且仅当 (A,B) 能控且 (A,C) 能观

给定一个传递函数

其中,分母多项式是 n 阶(系统阶数为 n),分子多项式是 m 阶(m ≤ n)。

写成能控标准型(CCF)

  • A 矩阵(n × n):

    A 的主对角线上是 0,次对角线上是 1,最后一行是分母多项式的负系数

  • B 矩阵(n × 1):

只有一个元素为 1,在最后一行。

  • C 矩阵(1 × n):基于分子多项式。如果分子阶数 m < n,则用 0 填充:

其中,填充到 n 个元素(如果 m < n-1,则末尾补 0)。

  • D 矩阵(标量):通常 D = 0,除非系统有直接馈通项。

这种形式确保了系统的能控矩阵是满秩的,从而系统是完全能控的。

写成能观标准型(OCF)

  • A 矩阵(n × n):是 CCF 的 A 矩阵的“转置式”变体,结构如下:

第一列是 ,最后一列是

  • B 矩阵(n × 1):基于分子多项式:

如果 m < n,则末尾补 0,使其成为 n × 1 向量。

  • C 矩阵(1 × n):
    即,只有一个元素为 1,在第一位。
  • D 矩阵:通常 D = 0。

这种形式确保了系统的能观矩阵是满秩的,从而系统是完全能观的。

对于单输入单输出 (SISO) 系统,其传递函数 的最小实现维数等于最简分式形式下分母 的阶数。当且仅当 prime(互质 ,即没有零极点相消)时,其标准能控型或能观型实现是最小实现