MATLAB频域分析：SISO模型特征提取与工程应用

1. 项目背景与核心价值

在控制系统工程领域，频域分析是评估系统动态特性的重要手段。作为一名长期从事自动控制系统设计的工程师，我经常需要从MATLAB的SISO（单输入单输出）模型频域响应数据中提取关键数值特征。这些特征包括幅值裕度、相位裕度、带宽频率等，它们直接决定了控制系统的稳定性和响应速度。

MATLAB的官方帮助文档虽然详尽，但对于非英语母语的工程师来说，直接阅读英文文档存在理解门槛。DeepSeek的翻译功能为技术文档的本地化提供了新思路，但如何准确提取和解读频域特征数值，仍是实际工程中的痛点。本文将结合我十年来的MATLAB使用经验，详解SISO模型频域特征值的提取方法与实际应用技巧。

2. 频域特征参数体系解析

2.1 基础特征参数定义

频域分析的核心在于以下几个关键参数：

1. 增益裕度(Gain Margin)：系统达到临界稳定时允许增加的增益量，通常以分贝(dB)表示。在伯德图中对应相位为-180°时的幅值。

2. 相位裕度(Phase Margin)：系统达到临界稳定时允许增加的相位滞后量，单位为度(°)。对应幅值为0dB时的相位与-180°的差值。

3. 带宽频率(Bandwidth Frequency)：系统幅频响应下降至-3dB时对应的频率，反映系统对输入信号的跟踪能力。

4. 谐振峰值(Resonant Peak)：幅频响应的最大值，表征系统相对稳定性。

2.2 MATLAB中的计算原理

MATLAB通过margin()函数实现这些参数的计算，其底层算法基于以下步骤：

1. 通过模型传递函数计算频率响应数据
2. 寻找相位穿越-180°的点（增益裕度点）
3. 寻找幅值穿越0dB的点（相位裕度点）
4. 通过插值法精确确定临界频率
5. 计算谐振峰值需要扫描整个频率响应曲线

注意：对于高阶系统，可能存在多个相位穿越点，MATLAB默认返回最保守（最小）的裕度值。

3. 实操：从SISO模型提取频域参数

3.1 基础数据获取方法

以典型二阶系统为例，演示完整操作流程：

matlab复制% 创建示例系统
sys = tf([1],[1 1.414 1]); % 自然频率=1rad/s，阻尼比=0.707

% 获取频域特征