3. 纵向运动建模：小扰动方程、传递函数与模态分析

各位同学，今天我们来聊聊纵向运动建模。说实话，这部分内容是飞控系统稳定性分析的基石。我个人习惯把纵向运动比作飞机的「呼吸」——它不像横航向那么剧烈，但处理不好，飞机就会像打摆子一样，忽上忽下。

嗯，咱们直接进入正题。

3.1 纵向小扰动方程

为什么要用小扰动？说白了，真实飞行中，飞机状态一直在变。我们不可能用一个方程描述所有情况。但如果我们假设飞机只在某个平衡点附近做微小运动，那问题就简化了。

我在项目中遇到过这样一个坑：有次做某型无人机建模，直接用了全量方程，结果仿真步长设得特别小，跑一次要半小时。后来换成小扰动方程，速度提升了十几倍，精度也完全够用。

Δẋ = A · Δx + B · Δu

其中：
Δx = [ΔV, Δα, Δq, Δθ]ᵀ   // 速度、迎角、俯仰角速率、俯仰角
Δu = [Δδe]                 // 升降舵偏角

A矩阵（纵向稳定性导数矩阵）：
A = [Xv  Xα  0   -g·cosθ₀]
    [Zv  Zα  1   -g·sinθ₀/V₀]
    [Mv  Mα  Mq  0        ]
    [0   0   1   0         ]

你想想看，这个A矩阵里的每个元素都有物理意义。比如Mα，就是迎角变化引起的俯仰力矩变化，我们叫它「纵向静稳定性导数」。如果Mα是负的，说明飞机是静稳定的——迎角增大，低头力矩产生，飞机自动恢复。

3.2 纵向传递函数

从状态方程到传递函数，其实就是拉普拉斯变换。我们关心的是：升降舵偏角Δδe到俯仰角速率Δq的传递函数。

为什么特别关注这个？因为俯仰角速率是飞控系统最常用的反馈信号。你去看任何一架电传飞机的控制律，俯仰速率反馈几乎是标配。

G(s) = Δq(s) / Δδe(s) = (s + 1/T₁)(s + 1/T₂) / [(s² + 2ζ₁ω₁s + ω₁²)(s² + 2ζ₂ω₂s + ω₂²)]

这个传递函数的分母是两个二阶环节的乘积。嗯，这里要注意：这两个二阶环节对应着两种完全不同的运动模态。

模态	典型周期	阻尼比ζ	物理表现
短周期	1~3秒	0.3~0.7	迎角和俯仰速率快速振荡
长周期	20~100秒	0.01~0.1	速度和俯仰角缓慢波动

我曾经在调试某型教练机时，发现短周期阻尼比只有0.15。飞行员反馈说「飞机像坐在弹簧上」，稍微动一下杆就晃个不停。后来我们调整了控制律的增益，把阻尼比提到0.5左右，问题就解决了。

3.3 长周期与短周期模态

这两种模态，说白了就是飞机纵向运动的两种「性格」。

短周期模态

短周期是飞行员最直接感受到的。你拉杆，飞机抬头，迎角增大，然后产生低头力矩——这个过程在几秒内完成。如果短周期阻尼太差，飞机会像「点头」一样，非常难受。

为什么会这样？从物理上看，短周期主要涉及迎角α和俯仰角速率q的耦合。迎角变化产生气动力矩，力矩改变俯仰速率，俯仰速率又反过来影响迎角。这是一个典型的二阶振荡系统。

长周期模态

长周期又叫「沉浮模态」。你想想看，飞机受到扰动后速度变化，升力跟着变，然后飞机开始缓慢地上下起伏。这个过程可能持续一分钟以上。

长周期的阻尼通常很弱。我记得有次做某大型运输机的仿真，长周期阻尼比只有0.02，几乎是无阻尼振荡。飞行员说「飞机像在海上漂」，虽然不会发散，但操纵起来很累。

从数学上看，长周期主要涉及速度V和俯仰角θ的耦合。速度变化导致升力变化，升力变化导致轨迹弯曲，轨迹弯曲又改变俯仰角，俯仰角再影响速度分量。这个循环很慢，所以周期长。

3.4 知识体系总览

为了让大家更直观地理解本章的知识结构，我画了一张图：

从这张图可以看得很清楚：小扰动方程是源头，传递函数是桥梁，最终落脚点是两种模态的分析。而短周期和长周期，又分别服务于控制律设计和飞行品质评估这两个工程目标。

好了，纵向运动建模的核心内容就这些。记住：小扰动方程是工具，传递函数是手段，模态分析才是目的。下次你调试飞控时，如果遇到纵向振荡问题，先判断是短周期还是长周期，然后对症下药。