传递函数与频域分析：从复域模型到系统稳定性判据
飞行控制系统的稳定性分析，本质上是在回答两个问题：系统在受到扰动后能否恢复平衡？以及系统在何种频率的输入下会表现出最危险的响应？传递函数与频域分析正是解决这两个问题的核心工具。传递函数将时域微分方程转化为复域代数方程，使系统特性一目了然；而频域分析则通过频率响应、伯德图和奈奎斯特图，将抽象的极点位置转化为直观的增益与相位关系。本章将系统讲解传递函数的定义与结构、极点零点对稳定性的影响、频率响应的物理意义，以及伯德图与奈奎斯特图的绘制与判读方法，为后续的稳定性判据（如奈奎斯特判据、相位裕度与增益裕度）奠定基础。
4.1 传递函数的定义与物理意义
传递函数是线性时不变系统在零初始条件下，输出量的拉普拉斯变换与输入量的拉普拉斯变换之比。对于飞行控制系统，无论是纵向俯仰通道还是横航向滚转通道，其动力学方程都可以表示为常系数线性微分方程。对微分方程两边取拉氏变换，并整理成输出与输入之比的形式，即可得到系统的传递函数。以典型的俯仰角控制为例，升降舵偏角δe到俯仰角速率q的传递函数通常具有如下形式：
G(s) = q(s)/δ_e(s) = (K(s + 1/T_θ)) / (s² + 2ζω_n s + ω_n²)其中，K为增益，1/Tθ为传递函数零点，分母的二次项对应短周期模态，ωn为自然频率，ζ为阻尼比。传递函数的分母多项式等于零即为特征方程，其根就是系统的极点。极点的实部决定了模态的衰减速率，虚部决定了振荡频率。飞行控制系统的稳定性完全由极点在s平面上的位置决定：所有极点均位于左半平面（实部小于零）时，系统渐进稳定；只要有一个极点位于右半平面，系统就不稳定；极点落在虚轴上则处于临界稳定状态。
4.2 极点零点与系统动态响应
极点与零点在s平面上的分布直接决定了系统的时域响应。对于飞行控制系统的短周期模态，一对共轭复数极点p1,2 = -ζωn ± jωn√(1-ζ²)对应着俯仰角速率的振荡响应。阻尼比ζ决定了振荡衰减的快慢，ζ越小，振荡越剧烈，衰减越慢；ζ大于1时，响应变为过阻尼，无振荡但响应迟缓。飞行品质规范（如MIL-F-8785C）对短周期阻尼比有明确要求，通常要求ζ在0.3到0.8之间，以保证飞行员有良好的操纵感受。
零点的影响则更为微妙。零点并不改变系统的稳定性，但会改变各极点模态在输出中的参与程度。例如，若传递函数中存在一个靠近原点的零点，它会削弱低频模态在响应中的比重，使系统对低频输入不敏感。在飞行控制中，升降舵到迎角的传递函数通常包含一个零点，该零点与短周期极点共同作用，决定了俯仰响应的初始斜率。若零点位于右半平面（非最小相位零点），系统在阶跃输入下会出现反向响应——这是飞控设计中需要特别避免的，因为它会导致飞行员产生错误的操纵直觉。
| 极点位置 | 时域响应特征 | 飞行控制含义 |
|---|---|---|
| 负实轴上（σ<0） | 指数衰减，无振荡 | 滚转模态，快速收敛 |
| 共轭复根（-ζωn±jωd） | 衰减振荡，ζ决定衰减率 | 短周期/荷兰滚模态 |
| 虚轴上（σ=0） | 等幅振荡，不衰减 | 临界稳定，需避免 |
| 右半平面（σ>0） | 发散，指数增长 | 不稳定，必须通过增稳系统修正 |
| 原点处 | 积分特性，阶跃响应斜坡 | 航向积分项，用于消除稳态误差 |
4.3 频率响应与伯德图
频率响应是传递函数在s = jω时的特例，它描述了系统对正弦输入的稳态响应特性。当输入为幅值A、频率ω的正弦信号时，线性系统的稳态输出也是同频率的正弦信号，但幅值变为A·|G(jω)|，相位滞后∠G(jω)。幅频特性|G(jω)|和相频特性∠G(jω)共同构成了系统的频率响应。伯德图将频率响应绘制在半对数坐标上：横轴为角频率ω（对数刻度），纵轴分别为幅值（dB）和相位（度）。这种表示方法的优势在于：
- 乘法运算转化为加法，串联环节的幅频特性直接相加，相频特性也相加。
- 采用渐近线近似，可以快速手绘幅频特性草图，便于初步分析。
- 对数刻度覆盖宽频率范围，从低频到高频都能清晰展示。
绘制伯德图的标准流程如下：首先将传递函数分解为典型环节（比例、积分、惯性、一阶微分、振荡环节）的乘积；然后分别绘制各环节的渐近幅频和相频曲线；最后将所有环节的曲线叠加，得到总特性。对于飞行控制系统，升降舵指令到俯仰角速率的典型伯德图在低频段呈现0dB水平线（比例特性），在中频段因短周期模态出现谐振峰（若ζ较小），在高频段以-40dB/dec斜率下降。穿越频率ωc（幅值曲线穿越0dB的频率）是衡量系统响应速度的关键指标，ωc越大，系统带宽越宽，响应越快，但也会放大高频噪声。
4.4 奈奎斯特图与稳定性判据
奈奎斯特图将频率响应绘制在复平面上，横轴为实部Re[G(jω)]，纵轴为虚部Im[G(jω)]，ω从0变化到∞时，轨迹形成一条曲线。与伯德图相比，奈奎斯特图能够直接应用奈奎斯特稳定性判据，判断闭环系统的稳定性。该判据的核心思想是：闭环系统稳定的充要条件是，开环传递函数G(s)H(s)的奈奎斯特曲线逆时针包围(-1, j0)点的次数N，等于开环传递函数在右半s平面的极点数P。即Z = P - N，当Z=0时闭环稳定。
奈奎斯特判据的工程价值在于：它允许我们仅通过开环频率响应实验数据（无需建立精确模型）来判断闭环稳定性。对于飞行控制系统，这一特性尤为重要——因为飞机在飞行中气动参数会随马赫数、迎角变化，精确建模困难。通过扫频飞行试验获得开环频率响应，绘制奈奎斯特图，即可评估闭环系统的稳定裕度。此外，奈奎斯特图还能直观显示系统对延迟的敏感度：延迟环节e-τs在奈奎斯特图上表现为轨迹的旋转，延迟越大，轨迹越靠近(-1, j0)点，稳定性越差。
4.5 案例拆解：某型无人机俯仰通道稳定性分析
某小型固定翼无人机在巡航状态下，升降舵到俯仰角速率的开环传递函数经辨识为：
G(s) = 12.5(s + 0.8) / (s² + 1.2s + 4)该传递函数对应短周期模态，自然频率ωn=2 rad/s，阻尼比ζ=0.3，零点位于s=-0.8。初始分析发现，该无人机在手动遥控模式下俯仰振荡明显，飞行员反馈"飞机点头不止"。通过频域分析定位问题：
- 绘制伯德图，发现谐振峰幅值达8.5dB（对应ζ=0.3），穿越频率ωc≈3.2 rad/s，相位裕度仅PM=22°，远低于45°的工程要求。
- 绘制奈奎斯特图，轨迹虽未包围(-1, j0)点（开环无右半平面极点，P=0，N=0，Z=0），但轨迹距离(-1, j0)点最近距离仅0.35，增益裕度GM≈9dB，但相位裕度严重不足。
- 进一步分析发现，零点s=-0.8靠近原点，在低频段产生相位超前，但同时也抬高了中频段幅值，加剧了谐振峰。
解决方案：在俯仰角速率反馈回路中增加超前-滞后校正网络，设计校正环节D(s) = (s+2)/(s+8)，将穿越频率处的相位提升约30°。校正后系统开环传递函数变为G(s)D(s)，重新绘制伯德图，相位裕度提升至PM=48°，增益裕度GM=11dB，谐振峰降至3dB以内。飞行试验验证：俯仰振荡明显收敛，飞行员评分从"不满意"提升至"满意"。
| 指标 | 校正前 | 校正后 | 工程要求 |
|---|---|---|---|
| 相位裕度 PM | 22° | 48° | ≥45° |
| 增益裕度 GM | 9 dB | 11 dB | ≥6 dB |
| 穿越频率 ωc | 3.2 rad/s | 2.8 rad/s | 1~4 rad/s |
| 谐振峰 Mr | 8.5 dB | 3.0 dB | ≤4 dB |
复盘经验：本案例中，仅凭时域响应（振荡）难以准确定位问题根源，而伯德图清晰揭示了相位裕度不足这一本质原因。校正网络的设计充分利用了频域方法的直观性——通过调整穿越频率附近的相位曲线来提升稳定裕度。值得注意的是，校正后穿越频率略有降低，响应速度稍有牺牲，但换取了更充足的稳定裕度，这在飞控设计中是合理的权衡。
4.6 常见误区与总结
在传递函数与频域分析的学习和工程应用中，以下误区需要特别警惕：
- 误区一：认为传递函数可以描述所有系统。实际上，传递函数仅适用于线性时不变系统。飞行器在大迎角、失速等非线性状态下，传递函数模型失效，需采用描述函数或非线性分析方法。
- 误区二：将零点与极点混为一谈。零点不影响稳定性，但显著影响动态响应形状。忽略零点可能导致对系统响应的误判，尤其是在非最小相位系统中。
- 误区三：只关注幅频特性而忽视相频特性。伯德图中幅值满足要求但相位裕度不足的情况非常常见，必须同时校验PM和GM。
- 误区四：奈奎斯特图判读时忘记考虑开环右半平面极点。若开环不稳定（P>0），即使奈奎斯特曲线不包围(-1, j0)点，闭环也可能不稳定。
- 误区五：将穿越频率设计得过高以追求快速性。过高的穿越频率会放大传感器噪声和执行机构的高频动态，反而降低系统鲁棒性。
总结而言，传递函数与频域分析为飞行控制系统稳定性分析提供了从复域到频域的完整工具链。传递函数建立了系统的数学模型，极点零点揭示了系统的固有特性，频率响应将抽象模型与可测量的物理量联系起来，伯德图和奈奎斯特图则分别从对数频率特性和极坐标轨迹两个视角给出了稳定性判据。掌握这些工具，工程师就能在设计中预判系统行为，在试验中定位问题根源，在调试中量化稳定裕度，从而确保飞行控制系统在各种飞行条件下都能安全、可靠地运行。