频域控制实战：基于MATLAB的伯德图与奈奎斯特图设计指南

1. 频域控制与MATLAB工具入门

第一次接触频域控制时，我也曾被那些复杂的图表搞得头晕眼花。直到真正用MATLAB画出了第一个伯德图，才突然明白这些曲线背后隐藏的系统秘密。频域分析就像给控制系统做"心电图"，通过观察系统在不同频率下的反应，我们能精准诊断其稳定性与性能。

MATLAB在这个领域就像医生的听诊器。它的控制系统工具箱提供了bode()、nyquist()等函数，让我们能快速绘制专业图表。比如下面这个简单的二阶系统分析：

s = tf('s');
G = 1/(s^2 + 0.5*s + 1);
bode(G)
grid on

短短四行代码就能生成完整的幅频和相频特性曲线。新手常犯的错误是直接看曲线却不关注坐标尺度，建议第一次使用时先右键图表选择"Properties"，将频率轴改为线性刻度更容易理解。

2. 伯德图实战解析

2.1 从理论到实践

记得我第一次调试无人机飞控时，发现系统总是轻微振荡。导师让我画出伯德图，这才发现相位裕度只有30度。伯德图由两幅子图组成：上方是幅频特性（dB vs Hz），下方是相频特性（度 vs Hz）。关键要关注三个特征点：

• 增益穿越频率：幅值曲线穿越0dB线的位置
• 相位穿越频率：相位曲线穿越-180度的位置
• 带宽频率：幅值下降3dB对应的频率

用margin()函数可以直接获取这些参数：

[Gm,Pm,Wcg,Wcp] = margin(G);
disp(['增益裕度:',num2str(20*log10(Gm)),'dB']);
disp(['相位裕度:',num2str(Pm),'度']);

2.2 设计实例：水箱液位控制

假设要设计一个水箱控制系统，要求超调量小于20%