2. 飞行器建模基础
飞行控制律设计的本质是“基于模型的控制”，因此一个准确、简洁且适合控制综合的飞行器模型，是整个设计流程的基石。本章将系统梳理从牛顿-欧拉方程出发建立刚体六自由度（6-DOF）运动方程的过程，深入剖析气动力与力矩模型的工程化表达方式，并重点讲解线性化与小扰动假设、以及配平（Trim）计算的核心方法。这些内容将直接为后续的控制律结构选择、增益调度设计和仿真验证提供数学与物理基础。
2.1 刚体六自由度运动方程：从物理定律到数学模型
核心原理：飞行器在空间中的运动可分解为质心的平移运动和绕质心的旋转运动，共六个自由度。建立运动方程的基本依据是牛顿第二定律和欧拉方程。在机体坐标系（Body Frame）下，力方程和力矩方程呈现出最简洁且便于控制设计的形式。需要特别注意的是，机体坐标系是一个动坐标系，因此必须引入由于坐标系旋转而产生的“表观力”项（如科氏力、离心力）。
力方程的矢量形式为：
m · (dV/dt + Ω × V) = F_aero + F_gravity + F_thrust其中，V = [u, v, w]ᵀ 为机体坐标系下的速度分量，Ω = [p, q, r]ᵀ 为角速度分量，F_aero 为气动力，F_gravity 为重力在机体轴上的投影，F_thrust 为发动机推力。力矩方程则表示为：
I · Ω̇ + Ω × (I · Ω) = M_aero + M_thrust这里 I 为惯性张量矩阵，对于对称布局的飞行器，通常可简化为对角矩阵（忽略惯性积）。运动学方程则描述了姿态角（滚转φ、俯仰θ、偏航ψ）与角速度之间的非线性转换关系，以及位置（北东地坐标系）与速度之间的积分关系。
2.2 气动力与力矩模型：工程化的系数表达
核心原理：气动力与力矩是飞行器所受外力的主要来源，其精确建模极为困难。工程上普遍采用“系数建模法”，即通过无量纲气动系数（CL, CD, CY, Cl, Cm, Cn）与动压 q̄ = ½ρV²、参考面积 S、参考弦长 c̄ 和翼展 b 的组合来描述气动力与力矩。这些系数通常是迎角 α、侧滑角 β、舵面偏度（升降舵 δe、副翼 δa、方向舵 δr）以及角速率（p, q, r）的非线性函数。
以纵向通道为例，升力系数和俯仰力矩系数可展开为：
CL = CL0 + CLα·α + CLδe·δe + (c̄/2V)·CLq·q
Cm = Cm0 + Cmα·α + Cmδe·δe + (c̄/2V)·Cmq·q其中，CLα 和 Cmα 是静稳定性导数，Cmδe 是操纵导数，CLq 和 Cmq 是动导数。这些导数的数值决定了飞行器的静稳定性和动稳定性特性。对于大迎角或高机动状态，还需要引入非线性修正项，如 α² 项或失速模型。
| 气动导数 | 物理含义 | 对稳定性的影响 | 典型获取途径 |
|---|---|---|---|
| Cmα | 俯仰力矩对迎角的导数 | 负值提供静稳定 | CFD / 风洞试验 |
| CLα | 升力线斜率 | 影响短周期频率 | 理论估算 / 试验 |
| Cmδe | 俯仰力矩对升降舵的导数 | 决定操纵效能 | 风洞试验 |
| Cnβ | 偏航力矩对侧滑角的导数 | 航向静稳定性 | CFD / 试验 |
| Clp | 滚转阻尼导数 | 抑制滚转速率 | 经验公式 / 试验 |
2.3 线性化与小扰动假设：从非线性到状态空间
操作流程：线性化的核心思想是在一个稳定的配平状态（x₀, u₀）附近，将非线性方程 f(x, u) 展开为一阶泰勒级数，忽略高阶小量。具体步骤如下：
- 确定配平点：求解 f(x₀, u₀) = 0，得到定直平飞或定常机动的状态和输入。
- 定义小扰动变量：令 Δx = x - x₀，Δu = u - u₀，并假设 |Δx| 和 |Δu| 足够小。
- 计算雅可比矩阵：求取状态方程对状态量和输入量的偏导数，得到系统矩阵 A 和控制矩阵 B。
- 输出方程线性化：对输出方程 y = g(x, u) 同样求偏导，得到 C 和 D 矩阵。
- 整理为标准状态空间形式：Δẋ = A·Δx + B·Δu，Δy = C·Δx + D·Δu。
对于常规固定翼飞行器，线性化后纵向与横航向通常可解耦，分别得到4阶和5阶的状态方程。纵向状态量通常选取 [ΔV, Δα, Δq, Δθ]，横航向选取 [Δβ, Δp, Δr, Δφ, Δψ]。
2.4 配平计算：寻找力的平衡点
核心原理：配平是线性化的前置条件，其本质是求解一组状态量和输入量，使得飞行器在特定飞行条件下满足力和力矩的平衡。对于定直平飞，需满足升力等于重力、推力等于阻力、俯仰力矩为零。配平计算通常采用数值迭代方法，如牛顿-拉夫森法。
以一个典型的定直平飞配平为例，已知飞行高度 H = 5000 m，马赫数 Ma = 0.6，需要求解的未知量为：迎角 α₀、升降舵偏角 δe₀、油门开度 δT₀。配平方程组为：
ΣFz = 0: L(α₀, δe₀, δT₀) - W = 0
ΣFx = 0: T(δT₀) - D(α₀, δe₀) = 0
ΣMy = 0: M(α₀, δe₀) = 0通过迭代求解上述方程组，即可获得配平状态。配平结果不仅是线性化的基准点，也是飞行性能分析、操纵面偏度校核的重要依据。
背景：某小型电动无人机，质量 m = 12 kg，翼展 b = 3.2 m，平均气动弦长 c̄ = 0.35 m，飞行速度 V = 30 m/s，飞行高度 H = 1000 m（ρ = 1.112 kg/m³）。
初始现象：在开环仿真中，飞行器无法保持高度，出现长周期振荡发散。
分析过程：首先进行配平计算，得到 α₀ = 4.2°，δe₀ = -2.5°，δT₀ = 0.35（归一化油门）。随后在配平点进行线性化，得到纵向状态矩阵 A 的特征值为：λ₁,₂ = -0.45 ± 2.1i（短周期），λ₃,₄ = 0.012 ± 0.18i（长周期）。发现长周期模态具有正实部，说明飞行器纵向静不稳定。
操作步骤：① 检查重心位置，发现重心过于靠后（X_cg = 0.32 c̄），导致静稳定裕度不足；② 将重心前移至 0.28 c̄；③ 重新配平与线性化，得到长周期特征值变为 -0.015 ± 0.16i，系统稳定。
结果对比：重心调整后，配平攻角变为 5.1°，升降舵偏角变为 -3.2°，但系统动态特性满足要求。
复盘经验：配平与线性化不仅是数学计算，更是发现物理问题的有效手段。特征值分析能快速暴露稳定性缺陷，指导总体参数调整。
2.5 常见误区与总结复盘
常见误区一：忽视惯性积的影响。对于非对称布局或载荷不对称的飞行器，惯性张量中的非对角项不可忽略，否则会导致横航向模型失真。
常见误区二：线性化点选取不当。若配平点计算错误或选取的飞行状态本身不满足平衡条件，则线性化模型将包含非零的“常数项”，导致状态空间模型不准确。
常见误区三：气动导数随马赫数变化处理不当。在跨声速区域，气动导数变化剧烈，若仍使用低速定值，将导致控制律在高马赫数下失效。
常见误区四：忽略执行机构动态。线性化模型通常只包含刚体动力学，但实际舵机、发动机存在带宽限制和延迟，若不在模型中体现，控制律设计可能过于乐观。