当前位置：首页 > news >正文

线性系统理论学懵了？手把手带你推导能控性格拉姆矩阵判据（附详细证明步骤）

news 2026/5/25 7:54:35

线性系统理论学懵了？手把手带你推导能控性格拉姆矩阵判据（附详细证明步骤）

第一次翻开《线性系统理论》教材时，那些抽象的数学符号和跳跃的证明步骤总让人望而生畏。特别是关于能控性的格拉姆矩阵判据，看似简单的定理背后隐藏着精妙的数学构造。本文将从一个学习者的视角，还原证明过程中的每一个关键思考节点，带你真正理解为什么格拉姆矩阵的非奇异性与系统能控性等价。

1. 理解基本概念：从零开始搭建知识框架

1.1 能控性的直观理解

能控性描述的是系统状态能否被外部输入所影响。想象驾驶一辆汽车：

完全能控：方向盘、油门和刹车可以让你到达任何位置（状态）
不完全能控：某些方向上的移动无法通过控制实现

数学上，线性时不变系统的状态方程为：

ẋ(t) = Ax(t) + Bu(t)

其中：

x(t) ∈ ℝⁿ是状态向量
u(t) ∈ ℝᵖ是控制输入
A ∈ ℝⁿˣⁿ和B ∈ ℝⁿˣᵖ是常值矩阵

1.2 格拉姆矩阵的构造原理

格拉姆矩阵在控制理论中扮演着重要角色，其定义为：

W_c(0, t₁) = ∫₀ᵗ¹ e^{-Aτ}BBᵀe^{-Aᵀτ} dτ

这个看似复杂的表达式实际上蕴含深刻的物理意义：

e^{-Aτ}表示系统动态的逆向演化
BBᵀ反映了控制输入的耦合强度
积分操作考虑了时间累积效应

提示：格拉姆矩阵实质上是描述控制能量在状态空间中的分布情况

2. 充分性证明：从格拉姆矩阵到能控性

2.1 证明思路构建

已知格拉姆矩阵非奇异，需要证明系统完全能控。核心思想是构造性证明——直接找到一个控制输入u(t)能将任意初始状态转移到原点。

关键步骤：

假设存在t₁>0使得W_c(0,t₁)可逆
构造控制输入：u(t) = -Bᵀe^{-Aᵀt}W_c⁻¹(0,t₁)x₀
验证该控制能将x₀在t₁时刻转移到0

2.2 详细推导过程

将构造的控制输入代入系统解：

x(t₁) = e^{At₁}x₀ + ∫₀ᵗ¹ e^{A(t₁-τ)}Bu(τ)dτ = e^{At₁}x₀ - ∫₀ᵗ¹ e^{A(t₁-τ)}BBᵀe^{-Aᵀτ}W_c⁻¹x₀ dτ

利用e^{A(t₁-τ)} = e^{At₁}e^{-Aτ}，可得：

x(t₁) = e^{At₁}[I - ∫₀ᵗ¹ e^{-Aτ}BBᵀe^{-Aᵀτ}dτ W_c⁻¹]x₀ = e^{At₁}[I - W_cW_c⁻¹]x₀ = 0

这个漂亮的构造直接验证了充分性条件。

3. 必要性证明：从能控性到格拉姆矩阵非奇异

3.1 反证法的运用

假设系统完全能控但格拉姆矩阵奇异，则存在非零向量η≠0使得：

ηᵀW_cη = ∫₀ᵗ¹ ||Bᵀe^{-Aᵀτ}η||² dτ = 0

这意味着对几乎所有τ∈[0,t₁]，有：

Bᵀe^{-Aᵀτ}η ≡ 0

3.2 能控性矛盾的导出

考虑将x₀ = e^{At₁}η转移到原点。由能控性定义，存在控制u(t)使得：

0 = e^{At₁}e^{At₁}η + ∫₀ᵗ¹ e^{A(t₁-τ)}Bu(τ)dτ

左乘ηᵀe^{-At₁}得：

0 = ηᵀη + ∫₀ᵗ¹ ηᵀe^{-Aτ}Bu(τ)dτ = ||η||²

这与η≠0矛盾，故原假设不成立。

4. 实践视角：MATLAB验证与数值案例

4.1 数值验证方法

通过MATLAB可以直观验证判据的有效性：

A = [1 2; 0 -1]; B = [1; 0]; t1 = 5; % 计算格拉姆矩阵 syms tau; Wc = int(expm(-A*tau)*B*B'*expm(-A'*tau), tau, 0, t1); disp('格拉姆矩阵：'); disp(vpa(Wc,4)); % 检查能控性 rank(ctrb(A,B)) == size(A,1) % 标准判据 det(Wc) ~= 0 % 格拉姆判据