基于Hanning窗插值FFT的电力谐波分析仪

资 源 简 介

本项目旨在解决电力系统信号在非同步采样情况下产生的频谱泄露和栅栏效应问题，从而实现对电网谐波的高精度分析。系统核心采用了基于Hanning窗的双谱线插值快速傅里叶变换（FFT）算法。由于电力系统频率通常在50Hz附近小幅波动，固定频率采样会导致计算误差。本程序通过对原始采样信号加Hanning窗，有效抑制频率泄露。在获取离散频谱后，利用频谱峰值及其相邻的最大谱线进行比值计算，并代入特定的插值修正公式，重新推算出真实的幅值、频率和初相位。该系统不仅能精确检测基波参数，还能对高次谐波进行定量的比例评估，适用于

详 情 说 明

电力系统谐波分析仪项目说明文档

1. 项目介绍

本项目设计并实现了一个专门用于电力系统信号处理的谐波分析仪。在实际电力系统中，频率往往在50Hz基准值附近波动，导致传统的固定频率采样产生非同步采样现象。非同步采样会引起频谱泄露和栅栏效应，严重影响谐波分析的精度。

为了解决这一问题，本项目采用了基于Hanning窗的加窗插值FFT算法。通过给采样信号加Hanning窗来抑制频谱泄露，并根据离散频谱中峰值谱线及其相邻谱线的幅值关系，利用插值修正公式精确推算出信号的真实频率、幅值和初相位。该方案能够显著提升电网在频率波动情况下的分析可靠性。

2. 功能特性

• 自适应频率校正：能够处理基波频率偏离50Hz的情况，自动修正因离散频率点与实际频率不重合产生的误差。
• 高次谐波定量分析：系统预设支持对基波（1次）及3、5、7、9、11次高次谐波的自动检测。
• 双谱线插值优化：利用Hanning窗的主瓣特性，通过双谱线重心法计算偏差值，实现参数的高精度补偿。
• 电能质量评估：自动计算总谐波畸变率（THD）以及各次谐波的含有率（HR）。
• 可视化交互界面：提供时域采样波形图和频域校正对比图，直观展示分析过程与结果。

3. 系统要求

• 运行环境：MATLAB R2016a 或更高版本。
• 计算资源：标准PC环境，程序采用轻量化矢量计算，响应速度快。

4. 核心功能设计逻辑

程序按照信号处理的标准流程，分为数据构造、预处理、变换、插值修正、评估输出五个环节。

#### 信号生成与环境模拟 系统设定采样频率为6400Hz，计算点数为4096点。模拟生成一个包含基波及多个奇次谐波（3、5、7、9、11次）的复合信号。其中基波频率设定为50.2Hz以模拟非同步采样，并加入适量的高斯白噪声以模拟真实的传感器环境。

#### 窗函数处理与频谱获取 在进行快速傅里叶变换（FFT）之前，对原始时域信号截断并施加Hanning窗。Hanning窗的表达式为 $w(n) = 0.5 times [1 - cos(2pi n / (N-1))]$。加窗能够有效降低截断引起的旁瓣电平，为后续的插值算法提供平滑的主瓣数据支持。

#### 双谱线插值算法实现 这是本程序的技术核心，包含以下步骤：

1. 峰值定位：在每个理论谐波频率附近搜索最大幅值线。
2. 邻域选取：通过比较峰值两侧谱线的大小，确定插值方向，选取包含最大能量的两根相邻谱线。
3. 比值计算：计算两根谱线的幅值比 $alpha$。
4. 修正量推导：利用Hanning窗特定的插值公式 $delta = (2alpha - 1) / (alpha + 1)$ 计算频率偏移系数。
5. 参数恢复：

* 频率修正：将离散频率点向 $delta$ 方向偏移。 * 幅值修正：考虑窗函数的能量衰减和插值损失系数进行补偿。 * 相位修正：修正由FFT复数角偏移和窗中心偏移引起的相位误差，并将结果归一化到正弦信号初相。

5. 关键算法细节分析

#### 频率与幅值修正公式 程序中通过算式 $delta = (2alpha - 1) / (alpha + 1)$ 获取频率偏差量，范围在 -0.5 到 0.5 之间。针对幅值修正，程序使用了一个包含 $sin(pi delta)$ 和 $pi delta (1-delta^2)$ 的修正因子。这种处理方式能够补偿由于离散采样点未落在频谱主线峰值上而导致的幅值衰减，极大地提升了幅值的恢复精度。

#### 相位补偿逻辑 FFT直接得到的是余弦相位。程序通过对结果进行 rad2deg 转换，并叠加 90 度的修正值，将结果统一为符合电力系统习惯的正弦表达式初相位。同时，程序包含了相位周期化处理（-180至180度），确保结果的美观和实用性。

#### 指标计算与报表打印 程序不仅识别物理参数，还通过内部逻辑自动计算：

• THD：基于所有被检测到的谐波有效值平方和的平方根与基波有效值之比。
• 含有率 (HR)：各次谐波幅值相对于基波幅值的百分比。

最终数据以格式化的控制台表格形式输出，包括频率、幅值、相位及含有率的完整清单。

#### 可视化呈现

• 时域图：动态显示采样信号的前10个周期波形，用以观察信号的畸变情况。
• 频域图：在同一坐标系下绘制原始FFT的灰色离散谱线和红色校正后的插值特征点。这种对比能够清晰展示算法如何从“发散”的谱线中精确定位到真实的谐波分布位置。