816 控制理论攻坚

8.5.2 线性连续定常系统的输出能控性

---

8.5.2 线性连续定常系统的输出能控性

一、核心考点与理论精要

1. 输出能控性的定义与物理意义

💡 【定义与概念】 输出能控性(Output Controllability)

考虑具有直接传输矩阵的线性连续定常系统:

{x˙(t)=Ax(t)+Bu(t)y(t)=Cx(t)+Du(t)\begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = Cx(t) + Du(t) \end{cases}

其中 x(t)Rnx(t) \in \mathbb{R}^n 为状态向量,u(t)Rpu(t) \in \mathbb{R}^p 为控制输入,y(t)Rqy(t) \in \mathbb{R}^q 为输出向量。

定义: 若在有限时间区间 [t0,tf][t_0, t_f] 内,存在容许的(无约束分段连续的)控制输入向量 u(t),t[t0,tf]u(t), t \in [t_0, t_f],能使系统由任意初始输出 y(t0)y(t_0) 转移到指定的任意最终输出 y(tf)y(t_f),则称此系统是输出完全能控的(简称输出能控)。

💡 【技巧与方法】 状态能控与输出能控的核心区别

  1. 关注维度不同:状态能控关注系统内部所有 nn 个状态变量 x1,x2,,xnx_1, x_2, \cdots, x_n 是否均可被控制;输出能控仅关注系统外部可测量的 qq 个输出变量 y1,y2,,yqy_1, y_2, \cdots, y_q 是否可由输入完全操纵。通常工业控制中 qnq \le n
  2. 直接通道的作用:当存在直接传输矩阵 D0D \neq 0 时,控制输入 u(t)u(t) 可直接跨过状态积分器瞬间改变输出 y(t)y(t),因此输出能控矩阵中包含了 DD 矩阵的分块。

2. 输出能控性判据(Kalman 准则)

💡 【THEOREM】 输出能控性判据

线性连续定常系统 (A,B,C,D)(A, B, C, D) 输出完全能控的充要条件是:输出能控性判别矩阵 SoS_o 具有满行秩(行秩等于输出向量维数 qq):

rank(So)=rank[CBCABCA2BCAn1BD]=q\text{rank}(S_o) = \text{rank} \begin{bmatrix} CB & CAB & CA^2B & \cdots & CA^{n-1}B & D \end{bmatrix} = q

其中:

  • SoRq×(np+p)S_o \in \mathbb{R}^{q \times (np + p)}
  • 若系统无直接传输项(即 D=0D = 0),则输出能控性判别矩阵退化为: So=[CBCABCA2BCAn1B]=CMS_o = \begin{bmatrix} CB & CAB & CA^2B & \cdots & CA^{n-1}B \end{bmatrix} = C \cdot M 其中 M=[BABAn1B]M = \begin{bmatrix} B & AB & \cdots & A^{n-1}B \end{bmatrix} 为系统的状态能控性矩阵。

3. 状态能控性与输出能控性的内在关系

💡 【THEOREM】 状态能控性与输出能控性的逻辑关系定理

  1. 状态能控推导输出能控(充分条件): 若系统是状态完全能控的(即 rank(M)=n\text{rank}(M) = n),且输出矩阵 CC 具有满行秩(即 rank(C)=q\text{rank}(C) = q),则系统必然是输出完全能控的
    • 推导依据:由矩阵秩不等式 rank(So)=rank(CM)=rank(C)=q\text{rank}(S_o) = \text{rank}(CM) = \text{rank}(C) = q
  2. 状态不能控 \centernot    \centernot\implies 输出不能控(非必要条件)系统状态不完全能控,系统输出依然可能是完全能控的!
    • 物理机理:当系统中那些不受控的状态分量没有耦合进输出通道 yy(即其在 CC 中对应的系数为 00),或者受控状态子空间的维数已经足以完全张成 qq 维输出空间时,输出变量仍然可以被输入完全操纵。

二、标准解题模型与计算步骤

💡 【METHOD】 输出能控性分析标准解题步骤

  1. 确定系统维度:明确状态维数 nn、输入维数 pp、输出维数 qq
  2. 计算分块矩阵乘积
    • 依次计算 CB,CAB,CA2B,,CAn1BCB, CAB, CA^2B, \cdots, CA^{n-1}B
  3. 构造输出能控性矩阵
    • 若有 DD,构造 So=[CBCABCAn1BD]S_o = \begin{bmatrix} CB & CAB & \cdots & CA^{n-1}B & D \end{bmatrix}
    • D=0D=0,构造 So=[CBCABCAn1B]S_o = \begin{bmatrix} CB & CAB & \cdots & CA^{n-1}B \end{bmatrix}
  4. 求秩与结论判定
    • SoS_o 进行初等行变换,求 rank(So)\text{rank}(S_o)
    • rank(So)=q\text{rank}(S_o) = q,判定系统输出完全能控
    • rank(So)<q\text{rank}(S_o) < q,判定系统输出不完全能控

三、经典例题与真题精析

例题 1:状态不可控而输出可控的典型系统

💡 【PROBLEM】

已知系统的状态方程和输出方程为:

x˙=[0112]x+[11]u,y=[10]x\dot{x} = \begin{bmatrix} 0 & 1 \\ -1 & -2 \end{bmatrix} x + \begin{bmatrix} 1 \\ -1 \end{bmatrix} u, \quad y = \begin{bmatrix} 1 & 0 \end{bmatrix} x

试判别系统的状态能控性与输出能控性,并说明两者的关系。

💡 【SOLUTION】 规范解答

解:

第一步:系统参数与维度确认 系统矩阵为:

A=[0112],B=b=[11],C=c=[10],D=d=0A = \begin{bmatrix} 0 & 1 \\ -1 & -2 \end{bmatrix}, \quad B = b = \begin{bmatrix} 1 \\ -1 \end{bmatrix}, \quad C = c = \begin{bmatrix} 1 & 0 \end{bmatrix}, \quad D = d = 0

状态维数 n=2n = 2,输出维数 q=1q = 1

第二步:判别状态能控性 计算 AbAb

Ab=[0112][11]=[0×1+1×(1)(1)×1+(2)×(1)]=[11]Ab = \begin{bmatrix} 0 & 1 \\ -1 & -2 \end{bmatrix} \begin{bmatrix} 1 \\ -1 \end{bmatrix} = \begin{bmatrix} 0 \times 1 + 1 \times (-1) \\ (-1) \times 1 + (-2) \times (-1) \end{bmatrix} = \begin{bmatrix} -1 \\ 1 \end{bmatrix}

构造状态能控性判别矩阵 SS(即 MM):

S=[bAb]=[1111]S = \begin{bmatrix} b & Ab \end{bmatrix} = \begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix}

计算行列式:

det(S)=1×1(1)×(1)=11=0\det(S) = 1 \times 1 - (-1) \times (-1) = 1 - 1 = 0

故:

rank(S)=1<n=2\text{rank}(S) = 1 < n = 2

结论 1: 系统状态不完全能控

第三步:判别输出能控性 计算矩阵乘积:

cb=[10][11]=1cb = \begin{bmatrix} 1 & 0 \end{bmatrix} \begin{bmatrix} 1 \\ -1 \end{bmatrix} = 1 cAb=[10][11]=1cAb = \begin{bmatrix} 1 & 0 \end{bmatrix} \begin{bmatrix} -1 \\ 1 \end{bmatrix} = -1

构造输出能控性矩阵 SoS_o

So=[cbcAbd]=[110]S_o = \begin{bmatrix} cb & cAb & d \end{bmatrix} = \begin{bmatrix} 1 & -1 & 0 \end{bmatrix}

其行秩为:

rank(So)=1=q\text{rank}(S_o) = 1 = q

结论 2: 系统输出完全能控

第四步:物理机理解析(为什么状态不可控但输出可控?) 展开状态微分方程:

{x˙1=x2+ux˙2=x12x2uy=x1\begin{cases} \dot{x}_1 = x_2 + u \\ \dot{x}_2 = -x_1 - 2x_2 - u \\ y = x_1 \end{cases}

观察状态组合:令 z1=x1+x2z_1 = x_1 + x_2,则:

z˙1=x˙1+x˙2=(x2+u)+(x12x2u)=x1x2=z1\dot{z}_1 = \dot{x}_1 + \dot{x}_2 = (x_2 + u) + (-x_1 - 2x_2 - u) = -x_1 - x_2 = -z_1

发现状态分量 z1(t)=x1(t)+x2(t)z_1(t) = x_1(t) + x_2(t) 的导数中输入项 uu 完全抵消,z˙1=z1\dot{z}_1 = -z_1 是完全自主衰减的不可控模态! 但是输出 y(t)=x1(t)y(t) = x_1(t),输入 u(t)u(t) 可以直接通过 x˙1=x2+u\dot{x}_1 = x_2 + u 驱动 x1x_1。因此,尽管系统内部存在不可控模态 x1+x2x_1 + x_2,输出 yy 本身仍然能够由控制量 u(t)u(t) 自由操纵至任意目标值。


四、易错陷阱与避坑指南

💡 【PITFALL】 考场三大易错点

  1. 满秩目标维度混淆
    • 状态能控性判别要求满秩目标为状态维数 nnrank(M)=n\text{rank}(M) = n);
    • 输出能控性判别要求满秩目标为输出维数 qqrank(So)=q\text{rank}(S_o) = q),切勿将输出能控矩阵与状态维数 nn 进行对比。
  2. 直接传输矩阵 DD 的遗漏
    • 若题目中给出了直接传输矩阵 DD(或标量 dd),构造 SoS_o 矩阵时务必在末尾并入 DD 矩阵(So=[CBCABCAn1BD]S_o = \begin{bmatrix} CB & CAB & \cdots & CA^{n-1}B & D \end{bmatrix})。虽然很多时候不影响秩的大小,但在考研大题中遗漏 DD 将会扣除步骤分。
  3. 因果逻辑颠倒
    • 牢记:“状态能控 ++ CC 满行秩     \implies 输出能控”是单向推导
    • “输出能控 \centernot    \centernot\implies 状态能控”,“状态不能控 \centernot    \centernot\implies 输出不能控”。

Discussion

Comments

Post

Thoughts, corrections, and follow-up notes are welcome here.