4. Ziegler-Nichols频域整定法：基于临界增益K_u和临界周期T_u的整定公式与步骤
说到PID整定，Ziegler-Nichols频域法是我最早学会的“实战派”方法。那时候我刚入行，师傅丢给我一句话：“先把系统搞振荡，再算参数。”我当时心里直犯嘀咕——把系统搞振荡？这不是找死吗？
后来我才明白，这方法的核心就是找到系统的“临界点”。说白了，就是让系统在纯比例控制下刚好开始等幅振荡。这个状态下的增益叫临界增益K_u，振荡周期叫临界周期T_u。有了这两个数，剩下的就是套公式了。
4.1 方法的核心思想
Z-N频域法的本质，是通过实验获取系统的频域特征。你想想看，一个闭环系统在临界稳定时，它的开环幅频特性刚好穿过0dB线，相频特性刚好穿过-180°线。这个点，就是系统的幅值裕度=1、相位裕度=0°的状态。
我个人习惯把这种方法叫做“振荡法”。因为它确实需要你故意让系统振荡起来。不过别怕，只要操作得当，这个振荡是可控的。
核心公式一览：
P控制器：K_p = 0.5 × K_u
PI控制器：K_p = 0.45 × K_u，T_i = 0.85 × T_u
PID控制器：K_p = 0.6 × K_u，T_i = 0.5 × T_u，T_d = 0.125 × T_u
4.2 实验步骤详解
嗯，这里要注意，整个实验过程必须保证系统安全。我在项目中遇到过一位同事，上来就把比例增益调到很大，结果系统直接发散，差点把执行器烧了。所以，请一定按步骤来。
- 准备工作：将PID控制器设置为纯比例模式（即T_i = ∞，T_d = 0）。初始比例增益K_p设为一个很小的值，比如0.1或0.2。
- 施加阶跃信号：给系统一个小的设定值阶跃，观察响应曲线。
- 逐步增大K_p：每次增加20%-30%，重复步骤2。直到系统出现等幅振荡。
- 记录临界参数：此时的K_p就是K_u。用秒表或示波器测量振荡周期，记为T_u。
- 套用公式：根据你需要的控制器类型，代入上面的公式计算参数。
⚠️ 重要警告：
我曾经见过有人把系统振荡幅度搞到满量程的80%，结果机械限位直接撞坏。记住：振荡幅度控制在满量程的5%-10%就足够了。一旦发现振荡幅度快速增大，立即减小K_p或切断控制信号。
4.3 整定公式的由来
你可能会问：为什么是0.6倍？为什么积分时间是0.5倍周期？
其实Z-N当年是通过大量实验总结出来的经验值。他们发现，对于大多数工业过程，这个比例能让系统获得约1/4衰减比的响应——也就是系统振荡三次后，振幅衰减到原来的1/4。这个衰减比在工程上被认为是一个不错的折中方案。
我个人的理解是：0.6倍的K_u给了你足够的“劲”，0.5倍的T_i帮你消除稳态误差，0.125倍的T_d则提前“刹车”防止超调过大。这三个数配合起来，就像老司机开车——油门、刹车、方向盘配合得刚刚好。
4.4 适用场景与局限性
这个方法不是万能的。我在实际项目中总结了几条经验：
- 适用：自衡对象（如温度、液位、压力），尤其是那些滞后时间较小的系统。
- 不适用：积分对象（如电机速度）、大滞后系统（滞后时间超过时间常数的1/3）。
- 慎用：不允许超调的场合（如化学反应器温度控制）。
说白了，这个方法适合“皮实”的系统。如果你的系统很娇贵，或者安全要求极高，我建议你改用后面章节会讲到的Cohen-Coon法或ISTE最优整定法。
4.5 一个实际案例
记得有一次，我给一个恒温箱做温控。加热器功率5kW，温度范围0-200°C。我按步骤操作：
- 初始K_p = 0.5，系统响应很慢，没有振荡
- 逐步增加到K_p = 3.2时，温度开始出现±1°C的等幅振荡
- 测量周期T_u ≈ 45秒
- 于是K_u = 3.2，T_u = 45s
- 套用PID公式：K_p = 0.6×3.2 = 1.92，T_i = 0.5×45 = 22.5s，T_d = 0.125×45 = 5.625s
整定完成后，系统在设定值50°C时，超调量约8%，调节时间约90秒，稳态误差±0.3°C。说实话，这个结果对于恒温箱来说已经够用了。但如果你需要更快的响应或更小的超调，可以在这个基础上微调。
4.6 知识体系总览
下面这张图是我自己整理的Z-N频域法知识框架，你可以把它当作一个快速索引：
💡 我的小技巧：
如果你手头有示波器或数据采集卡，建议把响应曲线录下来。这样你可以回放分析，精确测量T_u。我习惯在曲线上标记相邻两个波峰的时间差，取3-5个周期的平均值，这样更准。
另外，如果系统噪声比较大，可以在反馈通道加一个低通滤波器，截止频率设为1/T_u左右。这样能滤掉高频噪声，避免误判振荡。
好了，Z-N频域法的核心内容就这些。这个方法虽然老，但胜在简单实用。你只要记住：找到临界点，套公式，微调。三步走，搞定大部分常规系统。
公众号:蓝海资料掘金营，微信deep3321