MATLAB频域分析：控制系统稳定性评估实战指南

1. 项目背景与核心需求

在控制系统工程领域，频域分析是评估系统动态特性的重要手段。SISO（单输入单输出）模型的频域特性数值化表征，能直观反映系统的稳定性、响应速度等关键指标。这个项目源于我在帮同事调试一个工业温度控制系统时，发现MATLAB官方文档中关于频域特性数值提取的部分描述较为分散，而中文社区缺乏系统性的实操指南。

频域特性数值主要包括：幅值裕度(Gain Margin)、相位裕度(Phase Margin)、截止频率(Crossover Frequency)等。这些参数对判断控制系统鲁棒性至关重要。例如在无人机飞控系统设计中，相位裕度低于30°就可能出现振荡现象。通过MATLAB编程提取这些数值，可以快速评估数百个设计方案的稳定性，大幅提升开发效率。

2. 关键MATLAB函数解析

2.1 margin函数深度使用

margin函数是获取频域特性的核心工具，其完整语法为：

matlab复制[Gm,Pm,Wcg,Wcp] = margin(sys)

输出参数含义：

• Gm：幅值裕度（单位：dB）
• Pm：相位裕度（单位：度）
• Wcg：相位穿越频率（单位：rad/sec）
• Wcp：增益穿越频率（单位：rad/sec)

实际工程中建议使用带输出参数的调用方式，而非直接绘制Bode图。我曾遇到一个案例：某伺服系统在20Hz附近出现异常峰值，通过自动记录Wcg数值变化，最终定位到是编码器信号线受到变频器干扰。

2.2 多系统批量处理方法

对于需要比较多个控制器设计的场景，推荐使用cell数组存储结果：

matlab复制systems = {sys1, sys2, sys3}; 
results = cellfun(@margin, systems, 'UniformOutput', false);

这种写法比for循环效率高40%左右（实测i7-1185G7处理器处理100个系统耗时从0.8s降至0.47s）。注意要设置'UniformOutput'为false以处理不同阶次的系统。

3. 工程实践中的进阶技巧

3.1 不确定系统的鲁棒性评估

面对存在参数波动的系统（如电池模型中的内阻变化），可以结合Robust Control Toolbox进行蒙特卡洛分析：
`