Matlab新手必看：5分钟搞定THD计算（附电力系统谐波分析代码）

从“唱歌跑调”到精准量化：用Matlab驾驭电力谐波分析的实战指南

想象一下，你正在欣赏一场音乐会，歌手的声音本该纯净悠扬，但偶尔会冒出几个刺耳的音符，破坏了整体的和谐。在电力电子和电气工程的世界里，这种“跑调”现象同样存在，我们称之为谐波失真。电网或电力电子设备输出的理想正弦波，如果混入了频率是基波整数倍的“杂音”分量，就会导致波形畸变，影响电能质量，甚至损坏敏感设备。总谐波失真（Total Harmonic Distiation, THD）就是衡量这种“跑调”程度的黄金指标。

对于电力电子工程师、电气专业的学生，或是任何需要评估逆变器、变频器、UPS电源输出质量的从业者来说，快速、准确地计算THD是一项核心技能。过去，这可能需要复杂的公式推导和繁琐的频谱分析。但现在，借助Matlab强大的信号处理工具箱，我们可以在几分钟内完成从数据导入到结果可视化的全过程。本文将带你绕过理论深水区，直击实战核心，通过生活化的类比和可直接套用的代码模板，让你迅速掌握用Matlab进行电力系统谐波分析的要领。无论你是刚接触Matlab的新手，还是希望优化现有工作流的工程师，这里都有你需要的“干货”。

1. 理解THD：不仅仅是数学公式

在深入代码之前，我们有必要厘清THD究竟在衡量什么。简单来说，THD量化了信号中所有谐波分量总有效值与基波分量有效值的比值。一个完美的50Hz正弦波（基波）其THD为0%。一旦电路中存在非线性负载（如整流器、开关电源），就会产生100Hz、150Hz、200Hz等谐波，THD值随之升高。

THD的两种常见表达形式：

• 百分比形式：THD(%) = (sqrt(谐波功率总和) / 基波功率) * 100%。这是最直观的表示，例如“THD<3%”意味着谐波总含量低于基波的3%。
• 分贝形式：THD(dB) = 20 * log10(THD)。在音频和某些精密测量领域更为常见，能更好地表征大动态范围。

为什么THD如此重要？高THD意味着：

• 设备过热：谐波电流会导致电机、变压器等设备产生额外的铁损和铜损。
• 误动作：可能干扰基于过零检测的精密电子设备。
• 共振风险：可能与系统电容、电感形成谐振，放大特定次谐波，造成灾难性后果。

注意：在电力系统分析中，我们通常关注电流THD（THDi）和电压THD（THDv）。两者的成因和影响有所不同，但计算方法相通。本文示例主要围绕电压波形展开。

为了更清晰地对比不同场景下的THD典型值及其影响，可以参考下表：

2. 五分钟核心实战：调用Matlab内置thd()函数

Matlab的信号处理工具箱（Signal Processing Toolbox）提供了一个非常便捷的函数——thd()。对于初学者，这是最快速的上手途径。假设我们已经从仿真软件（如Simulink）或实际数据采集卡中，获得了一个逆变器输出电压的采样序列。

2.1 数据准备与基础分析

首先，我们需要将数据加载到Matlab工作区。数据可能来自.mat文件、.txt文件或.csv文件。

% 示例1：从 .mat 文件加载数据
load('inverter_output_voltage.mat'); % 假设文件包含变量 'voltage' 和 'Fs'
% 此时变量 voltage 和 Fs (采样频率) 应已存在于工作区

% 示例2：从文本文件加载数据（常见于从示波器导出）
data = load('oscillos