如何将相位角（phase degrees）转换为弧度？

Question

在信号处理或控制系统设计中，常需将相位角从度数（degrees）转换为弧度（radians），但许多初学者在使用MATLAB或Python等工具进行频域分析时，误将度数直接用于三角函数计算，导致结果错误。问题在于：这些函数默认输入为弧度而非度数。那么，如何正确地将相位角（如45°、90°）转换为弧度？转换公式是什么？是否可以借助编程语言内置函数实现自动转换？掌握这一基本转换对Bode图、相量分析和傅里叶变换等应用至关重要。

Answer 1

相位角从度数到弧度的转换：在信号处理与控制系统中的关键实践

1. 基础概念：角度单位的本质区别

在数学和工程计算中，角度有两种主要表示方式：度数（degrees）和弧度（radians）。度数将一个圆周划分为360等份，而弧度则是基于圆的半径长度定义的——1弧度是弧长等于半径时所对应的圆心角。

许多初学者在使用MATLAB或Python进行信号处理时，误将度数直接代入sin、cos等三角函数，导致频域分析结果出现偏差。这是因为这些函数默认接受的是弧度输入。

例如，sin(90) 在Python中计算的是 sin(90 radians) ≈ 0.894，而非预期的 sin(90°) = 1。

2. 转换公式：从度数到弧度的数学基础

度数与弧度之间的转换关系由以下公式给出：

• 弧度 = 度数 × (π / 180)
• 度数 = 弧度 × (180 / π)

以常见角度为例：

度数 (°)	弧度 (rad)
0	0
30	π/6 ≈ 0.524
45	π/4 ≈ 0.785
60	π/3 ≈ 1.047
90	π/2 ≈ 1.571
120	2π/3 ≈ 2.094
135	3π/4 ≈ 2.356
150	5π/6 ≈ 2.618
180	π ≈ 3.142
270	3π/2 ≈ 4.712
360	2π ≈ 6.283

3. 编程实现：利用内置函数自动转换

现代科学计算语言提供了便捷的内置函数来处理单位转换，避免手动计算错误。

在Python中（使用NumPy）：

import numpy as np

# 示例：将45度转换为弧度
deg = 45
rad = np.radians(deg)
print(f"{deg}° = {rad:.4f} rad")

# 反向转换
rad_back = np.degrees(rad)
print(f"{rad:.4f} rad = {rad_back}°")

# 正确使用三角函数
result = np.sin(np.radians(90))  # 输出 1.0

在MATLAB中：

% 示例：将90度转换为弧度并求正弦值
deg = 90;
rad = deg2rad(deg);
result = sin(rad);  % 正确结果为 1.0

% 或使用专门的函数
result_sind = sind(90);  % 直接处理度数的sin函数

4. 实际应用场景：Bode图与相量分析中的关键作用

在控制系统设计中，Bode图用于分析系统的频率响应，其中相位裕度常以度数表示。但在构建传递函数或进行数值仿真时，所有角度必须统一为弧度。

例如，在绘制一个二阶系统的Bode图时，若相位峰值为-90°，需将其转换为 -π/2 弧度参与计算。

在傅里叶变换中，复数表示法如 X(ω) = |X|∠θ 中的相位θ通常以弧度存储，尤其是在IFFT或相位校正算法中。

相量分析中，若多个信号的相位以度数给出（如 V₁=10∠30°, V₂=5∠-45°），在进行向量加法前，必须先转换为弧度以便于使用exp(1j*θ)形式表达。

5. 错误案例分析与调试建议

常见错误模式包括：

1. 直接调用 np.sin(90) 而未转换单位
2. 在滤波器设计中误设相位延迟为整数度数而未转为弧度
3. Bode图绘制时相位轴单位混乱，导致稳定性判断失误

调试建议：

• 始终检查输入参数单位是否符合API文档要求
• 在关键节点添加断言验证角度单位：assert abs(np.sin(np.radians(90)) - 1.0) < 1e-10
• 封装通用转换函数以提高代码可读性

6. 高级实践：构建鲁棒的相位处理模块

为提升工程健壮性，可设计如下Python工具类：

class PhaseConverter:
    @staticmethod
    def deg_to_rad(degrees):
        return degrees * np.pi / 180.0
    
    @staticmethod
    def rad_to_deg(radians):
        return radians * 180.0 / np.pi
    
    @staticmethod
    def safe_sin_deg(degrees):
        return np.sin(PhaseConverter.deg_to_rad(degrees))
    
    @staticmethod
    def phase_difference_deg(ph1_deg, ph2_deg):
        diff_rad = np.angle(np.exp(1j * PhaseConverter.deg_to_rad(ph1_deg - ph2_deg)))
        return PhaseConverter.rad_to_deg(diff_rad)

7. 流程图：相位角处理的标准工作流

graph TD
    A[原始相位数据] --> B{单位是度数?}
    B -- 是 --> C[调用 deg2rad 或 radians()]
    B -- 否 --> D[直接使用]
    C --> E[执行三角运算或频域分析]
    D --> E
    E --> F[Bode图/相量合成/FFT逆变换]
    F --> G[输出结果并验证单位一致性]