在自动控制原理、信号与系统以及电路分析等工程学科中，幅频特性曲线的相角怎么算是一个核心且常考的问题。许多初学者容易混淆“幅频特性”和“相频特性”的概念，误以为幅频曲线本身直接显示了相角。事实上，幅频特性描述了系统对不同频率信号的增益放大或衰减程度，而相角则包含在相频特性中。然而，两者同源，均取决于系统的传递函数。本文将深入解析如何从传递函数出发，精确计算特定频率下的相角，并探讨其在系统稳定性分析中的关键作用。
当线性定常系统受到正弦信号输入时，其稳态输出也是同频率的正弦信号，但幅值和相位可能发生变化。这种幅值和相位随频率变化的特性，称为系统的频率特性。它可以用复数形式表示为：
G(jω) = A(ω) e^(jφ(ω))
其中：
因此，回答“幅频特性曲线的相角怎么算”这个问题，本质上是在问：如何计算传递函数 G(jω) 的相位角 φ(ω)。
相角反映了系统对输入信号的时间延迟或超前效应。在负反馈控制系统中，相角的累积可能导致系统从负反馈变为正反馈，从而引发振荡甚至不稳定。因此，准确计算相角对于评估系统的相位裕度至关重要。
无论系统多么复杂，计算相角的基本逻辑是一致的。以下是标准的四步计算法：
首先必须明确系统的数学模型，即传递函数 G(s)。它通常是 s 的有理分式，由比例、积分、微分、惯性、振荡等环节串联或并联组成。
将传递函数中的复变量 s 替换为纯虚数 jω。这一步将时域模型转换为频域模型。注意，j 是虚数单位，满足 j² = -1。
将 G(jω) 进行复数运算，化简为标准形式：
G(jω) = Re(ω) + j·Im(ω)
其中 Re(ω) 是实部，Im(ω) 是虚部。这一步通常需要通过分子分母同乘共轭复数来实现。
利用反正切函数计算：
φ(ω) = arctan [ Im(ω) / Re(ω) ]
⚠️ 注意：必须根据 Re(ω) 和 Im(ω) 的正负判断象限，不能直接依赖计算器结果，因为 arctan 的值域通常仅在 (-90°, 90°) 之间。
假设有一个惯性环节：G(s) = 1 / (Ts + 1)
结论：惯性环节的相角范围从 0° (ω→0) 变化到 -90° (ω→∞)。
在实际工程中，大多数系统可以分解为几个典型环节的串联。掌握这些基本环节的相角特性，可以快速估算整个系统的相角。以下是常见环节的相角特性表：
| 环节名称 | 传递函数 G(s) | 频率特性 G(jω) | 相角 φ(ω) | 相角变化范围 |
|---|---|---|---|---|
| 比例环节 | K | K | 0° | 恒为 0° |
| 积分环节 | 1/s | 1/(jω) = -j/ω | -90° | 恒为 -90° |
| 微分环节 | s | jω | +90° | 恒为 +90° |
| 惯性环节 | 1/(Ts+1) | 1/(jTω+1) | -arctan(Tω) | 0° → -90° |
| 一阶微分 | Ts+1 | jTω+1 | +arctan(Tω) | 0° → +90° |
| 振荡环节 | 1/(T²s²+2ξTs+1) | 1/[(1-T²ω²)+j(2ξTω)] | -arctan(2ξTω/(1-T²ω²)) | 0° → -180° |
| 延迟环节 | e^(-τs) | e^(-jτω) | -τω (弧度) | -∞ → 0° (线性下降) |
对于高阶系统或包含延迟环节的系统，相角可能会超出 ±90° 的范围。此时，必须使用 atan2(y, x) 函数（在 MATLAB 中为 angle(G(jω))）来自动处理象限问题，或者手动根据实部和虚部的符号进行调整：
理解了单个环节的相角计算后，绘制整个系统的 Bode 图就变得简单了。Bode 图由两张图组成：上方是幅频特性图（dB），下方是相频特性图（度），横坐标均为对数频率 ω。
虽然渐近线法简便，但在转折频率附近误差较大。若需精确绘制相频曲线，应计算以下关键点：
将这些点连接成光滑曲线，即可得到精确的相频特性曲线。
对于包含非最小相位环节（如右半平面零点）或延迟环节的系统：
计算相角的最终目的是为了分析系统的稳定性。在频域分析中，我们主要关注两个指标：相位裕度和幅值裕度。
定义：在截止频率 ω_c（幅频特性为 0dB 的频率）处，相角 φ(ω_c) 与 -180° 之间的差值。
公式：
γ = 180° + φ(ω_c)
判据：
定义：在相位穿越频率 ω_g（相频特性为 -180° 的频率）处，幅频特性分贝值的负数。
公式：
Kg(dB) = -20lg|G(jω_g)|
判据：
💡 专家提示：在实际工程设计中，相位裕度比幅值裕度更常用，因为它直接反映了系统的相对稳定性及阻尼程度。如果“幅频特性曲线的相角怎么算”是为了评估稳定性，那么重点应放在 ω_c 处的相角计算上。
对于高阶系统，手工计算相角繁琐且易错。MATLAB 提供了强大的工具箱来辅助分析。
% 定义传递函数 G(s) = 10 / (s^2 + 2s + 10)
num = [10];
den = [1 2 10];
sys = tf(num, den);
% 计算频率 w = 2 rad/s 处的频率响应
w = 2;
[G_mag, G_phase] = bode(sys, w);
% G_phase 返回的是弧度，转换为角度
G_phase_deg = rad2deg(G_phase);
fprintf('在 w=%f 时，相角为: %f 度
', w, G_phase_deg);
bode(sys);
grid on;
title('Bode Diagram of G(s)');
% 自动标注相位裕度和幅值裕度
margin(sys);
通过 margin 命令，MATLAB 会自动在图上标出相位裕度和幅值裕度，极大简化了分析过程。
“幅频特性曲线的相角怎么算”并非一个孤立的问题，而是频域分析中的核心技能。关键在于：
通过理论计算与 MATLAB 仿真相结合，工程师可以高效地设计和优化控制系统，确保其在各种工况下的稳定性和性能。