3、谐波检测与测量技术：傅里叶变换（FFT）原理、加窗插值FFT、小波变换在谐波检测中的应用、谐波测量仪器选型

各位工程师朋友，大家好。今天我们来聊聊谐波检测与测量技术。说实话，这块内容是我在PCS项目里踩坑最多的地方。你想想看，谐波都测不准，后面的抑制措施就是瞎忙活。所以这一章，我把自己这些年摸爬滚打的经验都掏出来，咱们一起把这块硬骨头啃下来。

3.1 傅里叶变换（FFT）原理——谐波检测的基石

FFT，说白了就是把时域信号掰开揉碎，看看里面到底藏了哪些频率成分。我习惯把它比作「光谱分析」——白光经过三棱镜，就分出七色光；电流信号经过FFT，就分出基波和各次谐波。

数学上，离散傅里叶变换（DFT）的公式是：

X(k) = Σ[n=0 to N-1] x(n) * e^(-j*2π*k*n/N)

但实际工程中没人手算这个。我们直接用FFT，它只是DFT的一种快速算法。N点FFT的计算量从N²降到了N*log₂(N)。举个例子，1024点FFT，直接算要100多万次乘法，用FFT只要1万次左右。嗯，差距就是这么夸张。

我在项目中遇到过一个问题：采样率设得太高，点数又不够，结果50Hz和55Hz的谐波根本分不开。后来我学乖了，先算好Δf再设参数。

3.2 加窗插值FFT——对付频谱泄漏的利器

直接做FFT有个大坑——频谱泄漏。为什么会这样？因为FFT假设信号是周期无限的，但我们只能截取一段有限长的数据。截断就相当于加了个矩形窗，矩形窗的旁瓣会「污染」相邻频率。

我建议用加窗FFT。常用的窗函数有这些：

窗函数	主瓣宽度	旁瓣衰减	适用场景
矩形窗	最窄	-13dB	瞬态信号，或频率成分间隔很大
汉宁窗	较宽	-32dB	一般谐波检测，我最常用
海明窗	较宽	-43dB	旁瓣要求更高的场合
布莱克曼窗	最宽	-58dB	高精度测量，但频率分辨率牺牲大

加窗之后还有个问题：幅值会变。这时候就需要插值修正。我常用的方法是双谱线插值——找到峰值附近的两根谱线，用它们的比值来估算真实频率和幅值。

加窗插值FFT的流程大致如下：

1. 采集N点数据
2. 乘以窗函数（如汉宁窗）
3. 做N点FFT
4. 找峰值附近的几根谱线
5. 用插值公式修正频率、幅值、相位

3.3 小波变换在谐波检测中的应用

FFT有个硬伤：它假设信号是平稳的。但实际电网里的谐波，经常是突变的——比如大功率设备突然启动，谐波含量瞬间飙升。这时候FFT就抓瞎了，因为它只能告诉你「这段时间平均有多少谐波」，却说不清「谐波具体什么时候变的」。

小波变换就是来解决这个问题的。它用可伸缩的「小波基」去匹配信号，既能看全局，又能看局部。我习惯这么理解：FFT像一把固定刻度的尺子，小波像一把可以放大缩小的显微镜。

在谐波检测中，我主要用离散小波变换（DWT）。它的分解过程是这样的：

原始信号 → 低频近似分量(cA1) + 高频细节分量(cD1)
cA1 → cA2 + cD2
cA2 → cA3 + cD3
...

每一层分解，cD分量就对应一个频段的信号。比如采样率10kHz，第一层cD1对应2500-5000Hz，第二层cD2对应1250-2500Hz，以此类推。这样就能把各次谐波分到不同的频段里。

我曾经用db4小波做过一个PCS的暂态谐波检测。设备启动瞬间的谐波突变，小波变换抓得清清楚楚，FFT根本看不到这个细节。不过小波变换的计算量比FFT大不少，选型时要注意处理器的性能。

3.4 谐波测量仪器选型

讲完算法，咱们聊聊硬件。仪器选型这块，我踩过的坑最多，给大家分享几点实在的。

先看这张选型决策图，是我自己总结的：

选型时，我主要看这几个指标：

1. 精度等级：IEC 61000-4-30把仪器分为A、B、C三级。A级最高（0.5%），用于实验室和标准检测；B级（1%）用于一般现场；C级（5%）用于初步筛查。我个人建议，PCS项目至少用B级，关键场合用A级。
2. 测量带宽：至少要能测到50次谐波（2500Hz）。有些场合需要测到100次，比如有高频开关设备的场景。
3. 同步方式：锁相环同步是必须的。没有同步，相位测量就是扯淡。我见过有人用普通万用表测谐波，那结果根本没法看。
4. 数据处理能力：实时FFT、数据记录、通信接口（Modbus、IEC 61850）这些功能，按需选配。别为了省钱买了个「裸机」，后面加功能加不上去。

最后说一句，仪器只是工具，算法才是灵魂。再好的仪器，如果FFT参数设不对、窗函数选错了，测出来的数据也是垃圾。所以，先把前面讲的FFT和小波变换搞明白，再谈仪器选型。