Stevens《飞行器控制与仿真》读书笔记（二）：飞行器模态分析——纵向与横航向

本文是 Stevens《Aircraft Control and Simulation》读书笔记系列的第二篇，覆盖教材第2章的模态分析部分。阅读本文前需掌握6-DOF线性化和状态空间方程，以及特征根/阻尼比概念。

---

第二篇：飞行器模态分析——纵向与横航向

一、概念引入——从”稳不稳”到”怎么稳”

1.1 问题

你在PID调参文章中见过”自然频率”和”阻尼比”，但你有没有想过——这些值从哪来？

答案在第一篇的线性化状态空间方程中。系统矩阵 $\mathbf{A}$ 的特征值就是飞行器的固有力学模态。

控制设计的第一步不是”调参数”，而是理解被控对象本身的自然行为——即解耦后的纵向和横航向模态分析。

1.2 历史的五分钟

时间	人物	贡献
1911	布赖恩（Bryan）	首次将飞行器运动分解为纵向和横航向
1927	格劳尔特（Glauert）	给出了短周期和长周期的近似公式
1940s	美国NACA	系统化飞行器模态命名：短周期、荷兰滚等
1950s	航空工业标准	MIL-STD-8785 飞行品质规范——量化了各模态的要求

---

二、纵向模态

2.1 纵向状态空间

纵向近似解耦后的状态向量：

$$\mathbf{x}_{\text{long}} = [u, w, q, \theta]^T \quad \text{或} \quad [V, \alpha, q, \theta]^T$$

系统矩阵 $\mathbf{A}_{\text{long}}$ 的标准形式：

$$\begin{bmatrix} \dot{V} \\ \dot{\alpha} \\ \dot{q} \\ \dot{\theta} \end{bmatrix} = \begin{bmatrix} X_V & X_\alpha & 0 & -g\cos\theta_0 \\ Z_V/V_0 & Z_\alpha/V_0 & 1 & -g\sin\theta_0/V_0 \\ M_V & M_\alpha & M_q & M_\theta \\ 0 & 0 & 1 & 0 \end{bmatrix} \begin{bmatrix} V \\ \alpha \\ q \\ \theta \end{bmatrix}$$

实际上 $\mathbf{A}_{\text{long}}$ 的四个特征值分为两对共轭复根，对应两个模态：

2.2 短周期模态（Short Period）

特征：频率高、阻尼适中、衰减快（几秒内）

特征值形式：$\lambda_{sp} = -\zeta_{sp}\omega_{sp} \pm i\omega_{sp}\sqrt{1-\zeta_{sp}^2}$

近似公式（来自 Stevens 第2章）：

$$\omega_{sp} \approx \sqrt{-\frac{M_\alpha Z_\alpha}{V_0} + M_q \frac{Z_\alpha}{V_0} + M_\alpha}$$ $$\zeta_{sp} \approx -\frac{M_q + Z_\alpha/V_0 + M_{\dot{\alpha}}}{2\omega_{sp}}$$

物理本质：迎角和俯仰速率的快速振荡，主要受升降舵影响。

参数	典型值（通用航空）	典型值（战斗机）	典型值（四旋翼近似）
$\omega_{sp}$	2-5 rad/s	5-10 rad/s	20-40 rad/s
$\zeta_{sp}$	0.3-0.7	0.5-0.8	0.5-0.9

2.3 长周期模态（Phugoid）

特征：频率低、阻尼很弱、衰减慢（几十秒甚至不衰减）

特征值形式：$\lambda_{ph} = -\zeta_{ph}\omega_{ph} \pm i\omega_{ph}\sqrt{1-\zeta_{ph}^2}$

近似公式：

$$\omega_{ph} \approx \sqrt{2}\frac{g}{V_0}$$ $$\zeta_{ph} \approx \frac{1}{\sqrt{2}} \frac{C_D}{C_L}$$

物理本质：速度和俯仰角之间的慢速能量交换——飞机”上下起伏”的缓慢振荡。在高速飞行中，这个模态的阻尼非常弱。

2.4 纵向模态的物理直觉

想象你拉动升降舵让机头抬升：

1. 短周期：机头迅速上仰（$\theta$ 增大），迎角（$\alpha$）迅速增加，升力增大→机头上扬更快……然后俯仰阻尼 $M_q$ 把它拉回来。整个过程2-3秒结束。
2. 长周期：抬头后飞机爬升，速度降低→升力不足→机头下降→速度恢复→再次抬头……这个能量循环持续几十秒。

---

三、横航向模态

3.1 横航向状态空间

横航向近似解耦后的状态向量：

$$\mathbf{x}_{\text{lat}} = [v, p, r, \phi]^T$$

3.2 三种模态

横航向系统有三个特征值（不像纵向是对称的）：一个实根（滚转收敛）、一对共轭复根（荷兰滚）和一个实根（螺旋模态）。

3.2.1 荷兰滚模态（Dutch Roll）

特征：偏航和滚转的耦合振荡

近似公式：

$$\omega_{dr} \approx \sqrt{-\frac{N_v Y_\beta}{V_0} + N_\beta}$$ $$\zeta_{dr} \approx -\frac{Y_\beta/V_0 + N_r}{2\omega_{dr}}$$

物理本质：飞机像”鸭子走路”一样左右摇摆——偏航和滚转交替出现。这是乘客最不舒服的一种模态。

参数	典型值
$\omega_{dr}$	1-4 rad/s
$\zeta_{dr}$	0.05-0.3（经常太弱！）

飞行品质要求（MIL-STD-8785）：$\zeta_{dr} \ge 0.08$（1级），$\zeta_{dr} \ge 0.02$（2级）。低于这个值的飞机需要加装偏航阻尼器。

3.2.2 滚转收敛模态（Roll Convergence）

特征：纯粹的实根，指数衰减，没有振荡

特征值 $\lambda_{roll} = L_p$（滚转阻尼导数），总是负的。

时间常数：$\tau_{roll} = -1/L_p$

飞机类型	$\tau_{roll}$
大型客机	1-3 秒
战斗机	0.3-1 秒
小型无人机	0.1-0.3 秒

3.2.3 螺旋模态（Spiral Mode）

特征：也是实根，但可能为负（稳定）、接近零（中立）或为正（不稳定！）

物理本质：如图机稍微倾斜→侧滑→偏航→进一步倾斜…最终要么恢复、要么卷入螺旋俯冲。

大多数飞机在洗流配置下的螺旋模态是轻微不稳定的（特征值实部为正，但非常接近0）。飞行员可以通过微量副翼输入轻易控制。

---

四、模态分析的飞控意义

4.1 模态与控制器的关系

模态	主要受影响的状态	主要控制手段	设计者需要保证
短周期	$\alpha, q, \theta$	升降舵	$\zeta_{sp} \ge 0.3$
长周期	$V, \theta$	油门+升降舵	$\zeta_{ph} \ge 0.04$
荷兰滚	$v, p, r$	方向舵（偏航阻尼器）	$\zeta_{dr} \ge 0.08$
滚转收敛	$p$	副翼	$\tau_{roll} \le 1.0s$（1级）
螺旋	$\phi, \psi$	副翼（自动修正）	时间双倍 $\ge 20s$

4.2 与 PID 调参的联系

你在PX4 PID调参文章中调的是：

• 速率环（角速度）：直接影响短周期和滚转收敛模态的阻尼

  • $P_{\text{rate}}$ 增大 → $\omega_{sp}$ 增大（响应更快）
  • $D_{\text{rate}}$ 增大 → $\zeta_{sp}$ 增大（阻尼更强，但过大则响应迟钝）

• 姿态环（角度）：影响长周期和荷兰滚的阻尼

  • $P_{\text{att}}$ 增大 → 长周期模态频率增大

• 位置环：把飞行器的12状态系统作为一个跟踪控制器设计

模态分析告诉你：为什么飞控的分层结构（速率→姿态→速度→位置）是自然的——因为每一层对应的物理模态的时间尺度本来就不同。

---

五、完整概念地图

12状态系统
     ↓ 小扰动线性化
    A矩阵 (12×12)
     ↓ 纵向/横航向解耦
    
纵向 (4-5阶)                         横航向 (4-5阶)
    │                                      │
    ├── 短周期 (高频, 2-5rad/s)            ├── 荷兰滚 (1-4rad/s, 弱阻尼)
    │   α, q 振荡 → 升降舵可调             │   v, p, r 耦合 → 偏航阻尼器
    │                                      │
    ├── 长周期 (低频, 0.3-1rad/s)          ├── 滚转收敛 (实根, 快)
    │   V, θ 能量交换 → 弱阻尼             │   p 指数衰减 → 滚转响应
    │                                      │
    └── (高度: 积分不改变特征值)            └── 螺旋 (实根, 通常略不稳)
                                              φ, ψ → 飞行员可修正

控制设计含义:
    速率环 → 短周期 + 滚转收敛
    姿态环 → 长周期 + 荷兰滚
    位置环 → 外回路, 时间尺度更大

---

六、核心公式速查卡

模态	特征值形式	估算公式	设计要求
短周期	$-\zeta_{sp}\omega_{sp} \pm i\omega_{sp}\sqrt{1-\zeta_{sp}^2}$	$\omega_{sp} \approx \sqrt{-M_\alpha Z_\alpha/V_0}$	$\zeta \ge 0.3$
长周期	$-\zeta_{ph}\omega_{ph} \pm i\omega_{ph}\sqrt{1-\zeta_{ph}^2}$	$\omega_{ph} \approx \sqrt{2}g/V_0$	$\zeta \ge 0.04$
荷兰滚	$-\zeta_{dr}\omega_{dr} \pm i\omega_{dr}\sqrt{1-\zeta_{dr}^2}$	$\omega_{dr} \approx \sqrt{-N_v Y_\beta/V_0 + N_\beta}$	$\zeta \ge 0.08$
滚转收敛	$\lambda = L_p$（实根）	$\tau_{roll} = -1/L_p$	$\tau \le 1.0s$
螺旋	$\lambda$（实根，可正可负）	—	时间双倍 $\ge 20s$

---

参考文献

1. Stevens, B. L., Lewis, F. L., & Johnson, E. N. (2015). Aircraft Control and Simulation (3rd ed., Chapter 2). Wiley.
2. MIL-STD-8785C (1980). Military Standard: Flying Qualities of Piloted Aircraft.
3. Cook, M. V. (2012). Flight Dynamics Principles (3rd ed., Chapter 6-8). Butterworth-Heinemann.
4. Etkin, B., & Reid, L. D. (1995). Dynamics of Flight (3rd ed., Chapter 7-8). Wiley.

---

下一节：PID控制与飞控实现

从模态分析到控制器设计：单回路PID、串级PID、内环/外环结构——将你已有的PID调参文章放在Stevens的框架下重新审视，并与真实的PX4飞控代码对应。