3、频域分析基础:传递函数与频率响应、伯德图与奈奎斯特图、幅值裕度与相位裕度
各位同学,咱们今天聊聊频域分析。说实话,我刚入行那会儿,觉得时域分析挺直观的——给个阶跃信号,看响应曲线就完了。但后来做飞控调试,发现光看时域根本不够。为什么?因为飞机在天上飞,遇到的是各种频率的扰动,比如阵风、结构振动、舵面抖振。这时候,频域分析就成了我们的“透视镜”。
3.1 传递函数与频率响应
先说说传递函数。说白了,它就是系统在复频域(s域)的数学模型。你想想看,一个线性时不变系统,输入输出之间的关系,用微分方程描述很麻烦,但拉普拉斯变换一搞,就变成了一个简单的代数分式。
比如一个典型的二阶系统:
G(s) = ωn² / (s² + 2ζωn·s + ωn²)
这里ωn是自然频率,ζ是阻尼比。我在项目中遇到过,很多新手只盯着阻尼比看,觉得阻尼比调到0.7就万事大吉。其实不然,自然频率决定了系统的响应速度,两者要一起看。
那频率响应是什么呢?很简单,你把传递函数中的s换成jω(j是虚数单位,ω是角频率),就得到了频率响应函数:
G(jω) = |G(jω)| · e^(j·∠G(jω))
这里|G(jω)|是幅值,∠G(jω)是相位。嗯,这里要注意:频率响应描述的是系统对不同频率正弦输入的稳态响应。你给系统输入一个正弦信号,输出也是同频率的正弦,只是幅值和相位变了。
核心理解:传递函数是系统的“基因”,频率响应是系统的“体检报告”。基因决定了你能长多高,体检报告告诉你现在状态如何。
3.2 伯德图与奈奎斯特图
有了频率响应,怎么直观地看?两个工具:伯德图和奈奎斯特图。
伯德图由两张图组成:幅频特性图和相频特性图。横轴都是频率(对数坐标),纵轴分别是幅值(dB)和相位(度)。我个人习惯,调试飞控时先看伯德图,因为它能直观告诉你:
- 低频增益够不够(决定了稳态精度)
- 穿越频率在哪里(决定了响应速度)
- 高频衰减特性(决定了抗噪声能力)
举个例子,一个典型的PID控制器伯德图:
Gc(s) = Kp + Ki/s + Kd·s
低频段积分项起主导,增益很高,保证稳态精度;中频段比例项主导,决定穿越频率;高频段微分项起主导,提升相位。我曾经调一个四旋翼的角速度环,光看时域曲线怎么调都有超调,后来画出伯德图才发现,穿越频率附近相位裕度只有15度,太低了。
奈奎斯特图则是把频率响应画在复平面上。横轴是实部,纵轴是虚部,频率从0到无穷大,画出一条曲线。为什么要用这个?因为它能直接判断闭环系统的稳定性——看曲线是否包围(-1, j0)点。
我的小技巧:伯德图适合“设计”,奈奎斯特图适合“验证”。先用伯德图调参数,再用奈奎斯特图确认稳定性。两个工具配合使用,事半功倍。
3.3 幅值裕度与相位裕度
这两个指标,是频域分析中最实用的两个数。说白了,它们告诉你:你的系统离“不稳定”还有多远。
幅值裕度(Gain Margin, GM):在相位穿越频率(相位为-180°的频率点)处,系统还能增加多少增益才会变得不稳定。一般要求GM > 6dB。
相位裕度(Phase Margin, PM):在幅值穿越频率(幅值为0dB的频率点)处,系统的相位离-180°还有多远。一般要求PM在30°~60°之间。
我建议你记住这个经验值:
| 指标 | 推荐值 | 说明 |
|---|---|---|
| 幅值裕度 GM | > 6 dB | 太小容易振荡,太大响应太慢 |
| 相位裕度 PM | 30° ~ 60° | 30°以下超调大,60°以上响应迟钝 |
| 穿越频率 ωc | 系统带宽的0.5~0.8倍 | 决定了响应速度 |
避坑指南:我曾经调一个固定翼的俯仰通道,相位裕度调到45度,觉得稳了。结果试飞时遇到强阵风,飞机开始剧烈振荡。后来分析发现,我忽略了舵机执行器的延迟——它给系统额外引入了相位滞后。所以,建模时一定要把执行器、传感器延迟都考虑进去,否则你的裕度全是“虚的”。
最后总结一下频域分析的三个步骤:
- 建模:写出系统的传递函数,包括被控对象、执行器、传感器
- 绘图:画出伯德图,找到穿越频率和相位裕度
- 校验:用奈奎斯特图确认稳定性,调整参数直到裕度满足要求
记住,频域分析不是纸上谈兵。你调出来的每一个参数,最终都要在实物上验证。但有了频域这个工具,你至少知道“往哪个方向调”,而不是瞎蒙。嗯,今天就到这里,下节课我们讲如何用MATLAB工具进行频域分析实战。