3、数学预备知识(二):动态系统理论——状态空间表示、传递函数、稳定性与可控可观性

各位同学,欢迎来到动态系统理论这一讲。

说实话,搞飞行器控制这么多年,我最大的体会就是:你如果不懂动态系统,你连飞机为什么抖都说不清楚。今天咱们就把这块硬骨头啃下来。

3.1 状态空间表示——把系统装进一个“盒子”里

先问大家一个问题:一个飞行器在空中飞,你怎么描述它此刻的状态?

位置?速度?姿态角?角速度?没错,这些就是状态变量

状态空间表示法,说白了就是用一个一阶微分方程组,把系统的全部动态信息装进去。数学形式长这样:

ẋ(t) = A·x(t) + B·u(t)
y(t) = C·x(t) + D·u(t)

这里:

  • x(t) —— 状态向量,比如 [位置, 速度, 俯仰角, 俯仰角速率]
  • u(t) —— 输入向量,比如舵面偏角、油门
  • y(t) —— 输出向量,比如传感器测到的数据
  • A, B, C, D —— 系统矩阵,决定了系统的“脾气”

我个人习惯:拿到一个新模型,第一件事就是写出它的状态空间方程。哪怕后面要转传递函数,也得先有这个“骨架”。

举个例子,一个简单的俯仰通道短周期运动,状态变量取攻角α和俯仰角速率q,输入是升降舵偏角δe:

[α̇]   [Zα  1+Zq] [α]   [Zδe]
[q̇] = [Mα  Mq ] [q] + [Mδe]·δe

你看,就这么简洁。我在项目中遇到过好几次,有人拿着几十页的微分方程推导报告,我让他先写成状态空间形式,瞬间就清晰了。

3.2 传递函数——从“时域”跳到“频域”

状态空间是时域的描述,但很多时候我们想知道:输入一个正弦波,输出会怎样?这时候就需要传递函数了。

传递函数怎么来的?对状态空间方程做拉普拉斯变换,假设初始条件为零:

G(s) = C·(sI - A)⁻¹·B + D

嗯,这个公式看着有点吓人。但你想想看,它本质上就是输出/输入的比值,只不过是在复频域里。

传递函数的好处是直观:

  • 极点 —— 决定系统稳定性(后面细说)
  • 零点 —— 影响系统响应形状
  • 增益 —— 决定稳态输出大小

避坑指南:我曾经在分析一个飞翼布局的横航向稳定性时,直接用传递函数算,结果发现算出来的零点位置和实际试飞数据对不上。后来才发现,是忽略了非最小相位特性。所以记住:传递函数只适用于线性时不变系统,而且不能处理初始条件

3.3 稳定性——系统会不会“炸”?

搞飞控的,最怕什么?最怕飞机发散。稳定性就是回答这个问题:给系统一个扰动,它能自己回来吗?

判断稳定性,有几种常用方法:

方法 适用场景 核心思想
特征值法 状态空间模型 矩阵A的特征值实部全部为负 → 稳定
极点法 传递函数模型 所有极点位于左半平面 → 稳定
李雅普诺夫法 非线性系统 找一个能量函数,看它是否衰减
奈奎斯特判据 频域分析 开环频率响应绕(-1, j0)点的圈数

我个人最常用的是特征值法。为什么?因为状态空间模型里,A矩阵的特征值直接告诉你每个模态的阻尼比和自然频率。比如:

  • 特征值 λ = -0.5 ± 2i → 阻尼比0.24,频率2.06 rad/s → 短周期振荡,阻尼偏弱
  • 特征值 λ = -0.02 → 长周期模态,收敛很慢
  • 特征值 λ = 0.1 → 不稳定!赶紧改控制律

注意:特征值实部为负只是“渐近稳定”的条件。如果实部为零,那是临界稳定,实际工程中绝对不能接受——稍微有点非线性因素就发散。

3.4 可控性与可观性——你能管住它吗?你能看到它吗?

这两个概念,是卡尔曼老爷子提出来的。我当年学的时候觉得抽象,后来做项目才明白:这是决定你能不能设计控制器的前提

可控性

问:给定任意初始状态,能不能通过输入u(t)在有限时间内把系统拉到任意目标状态?

判断方法:构造可控性矩阵

Qc = [B, AB, A²B, ..., Aⁿ⁻¹B]

如果Qc满秩(秩等于状态维数n),则系统可控。

举个例子,我之前做倾转旋翼机时,发现俯仰通道和滚转通道在某些飞行模式下耦合严重,可控性矩阵的秩下降了一阶。这意味着有一个状态方向是控制不了的。后来我们调整了舵面布局才解决。

可观性

问:给定一段时间内的输出y(t)和输入u(t),能不能唯一确定系统的初始状态?

判断方法:构造可观性矩阵

Qo = [C; CA; CA²; ...; CAⁿ⁻¹]

满秩即可观。

我建议:每次设计状态观测器之前,先算一下可观性矩阵的秩。我曾经见过一个同事,花了两周设计卡尔曼滤波器,结果发现系统不可观——白干了。

3.5 知识体系总览

下面这张图,是我自己总结的动态系统理论核心逻辑,画出来方便大家理解:

动态系统理论核心知识体系 动态系统模型 时域:状态空间表示 频域:传递函数 A矩阵特征值 可控性矩阵 可观性矩阵 极点位置 零点位置 奈奎斯特判据 稳定性 → 可控性 → 可观性 三者共同决定控制器设计可行性

这张图把今天讲的内容串起来了。你看,状态空间和传递函数是两条主线,最终都指向稳定性、可控性、可观性这三个核心问题

搞飞行器控制,说白了就是:先建模(状态空间或传递函数),再分析(稳定性),最后设计控制器(依赖可控可观性)。每一步都离不开今天讲的这些基础。

最后分享一个经验:我每次拿到一个新机型的数据,都会先做三件事——算特征值看稳定性,算可控性矩阵看能不能控,算可观性矩阵看能不能观。这三步走完,心里就有底了。你也试试这个习惯,保证少走弯路。


公众号:蓝海资料掘金营,微信deep3321