基于谱序列变换的高精度谐波参数估计算法-有色金属在线

分析是电力系统谐波污染治理的前提，由于频谱泄漏的影响，经典的谐波分析算法存在参数估计精度不高、计算复杂度较大等问题，为此，在分析信号谱线衰减特征的基础上，提出一种基于谱序列变换的高精度谐波参数分析算法。该算法通过对信号进行FFT运算得到谱序列并对其实施特定的加权变换，加速非真实频率处谱线的衰减速度，从而达到有效抑制频谱泄漏的目的，在此基础上推导谐波幅值和相位的计算公式。研究结果表明，与经典的加窗插值算法和FFT改进算法相比，所提出的算法所得谐波参数估计精度显著提高，且其在基波频率变化及间谐波条件下均表现出优良的估计性能；此外，与经典的加窗插值算法相比，新算法只需对谱序列进行简单处理，具有计算复杂度低的优点。

关键词：

谐波分析；快速傅里叶变换；谱序列变换；参数估计；频谱泄漏；

中图分类号：TN911.72；TM935.1 文献标志码：A 开放科学(资源服务)标识码(OSID)

文章编号：1672-7207（2020）09-2504-10

High precision harmonic parameter estimation algorithm based on spectral sequence transformation

YANG　Xi^{1, 2}, WANG　Xuming³, CHEN　Bingquan¹, ZHANG　Renmin¹, WANG　Xiangming⁴, LEI　Kejun¹

(1. College of Information Science and Engineering, Jishou University, Jishou 416000, China;

2. School of Information Science and Engineering, Southeast University, Nanjing 210096, China;

3. College of Physics and Electromechanical Engineering, Jishou University, Jishou 416000, China;

4. Shanghai Stock Exchange Technology Co. Ltd., Shanghai 200120, China)

Abstract: High-precision harmonic analysis is the premise of harmonic pollution control in the power system. Due to the influence of spectrum leakage, the classical harmonic parameter estimation algorithms have some disadvantages, such as low accuracy of parameter estimation and high computational complexity.Based on the analysis of spectral line attenuation characteristics, a high-precision harmonic analysis algorithm based on spectral sequence transformation was proposed. The algorithm performed a specific weighted transformation on the spectral sequence obtained via fast Fourier transform(FFT), which accelerated the attenuation rate of spectral lines at the unreal frequencies, and thus effectively suppressed the spectrum leakage. Based on this, an effective estimation method of harmonic amplitude and phase was proposed. The results show that, compared with the classical windowed interpolation algorithm and the improved FFT algorithm, the proposed algorithm in the paper improves the accuracy of parameter estimation significantly, and shows excellent performance when the fundamental frequency changes or inter-harmonics exists. In addition,compared with the classical windowed interpolation algorithms, the proposed algorithm only needs simple processing of spectral sequence, and has the advantage of low computational complexity.

Key words: harmonic analysis; fast Fourier transform(FFT); spectral sequence transformation; parameter estimation; spectrum leakage

随着现代电力电子技术的快速发展，越来越多的非线性元器件被应用于电力系统^[1]，在电网中产生大量谐波污染。谐波污染使得电能质量下降，从而严重影响电力系统的安全和经济运行。为了更好地进行谐波污染治理，必须对谐波参数进行高精度估计^[2]。快速傅里叶变换(fast Fourier transform, FFT)是一种经典的谐波参数估计方法。然而，由于在实际运行过程中电网基波处于波动状态，因此，实际的采样过程通常处于非同步状态。直接采用FFT算法对谐波信号进行参数估计，其过程相当于对原始信号序列施加矩形窗。矩形窗实现方式虽然简单，但这种窗所引起的频谱泄露非常严重，其旁瓣峰值衰减仅为13 dB，由此导致FFT算法对应的估计结果精度较低。为了减小频谱泄漏在谐波参数估计过程中带来的负面影响，采用在时域加窗处理和频域谱线插值处理相结合的方法来提高估计精度^[3^-13]。在前期的时域加窗处理过程中，常用的窗函数包括Hanning窗、Blackman窗、Nuttall窗、Rife-Vincent窗和MSD窗。这些窗函数相对简单，但其旁瓣衰减速度有限而不利于抑制频谱泄露。为进一步加速窗函数旁瓣衰减，近年来，一些更加复杂的优化窗相继被提出，如乘法窗^[11]、卷积窗^[12^-13]等。然而，当采用一些复杂的窗函数时，频谱泄露的抑制效果尽管有所提高，但另一方面也将不可避免地导致更高的计算复杂度而不利于信号参数的实时估计。通过加窗处理可以在一定程度上抑制长距离频谱泄漏的影响，但其对短距离范围内频谱泄露的抑制效果一般，因此，在随后的频域谱线插值处理过程中通过选择峰值谱线附近不同数量的谱线进行插值，以达到减小短距离频谱泄漏的影响，从而进一步提高参数估计精度。例如，文献[3，4，8]提出了加窗双谱线插值算法，文献[5，9]提出了加窗三谱线插值算法，文献[6，10]提出了加窗四谱线插值算法，文献[7]则提出了一种加窗六谱线插值算法。这些插值方法在一定程度上都可以提高谐波参数估计的精度。值得注意的是，上述基于不同数量谱线的插值方法有着各自的特点：对于利用较少谱线进行插值的方法，其插值公式推导和相关参数数值计算的复杂度较低，但其谐波参数估计的精度也相应较低；而对于利用较多谱线进行插值的方法，虽然参数估计的精度有所提高，但其计算复杂度也相应提高。相比较而言，文献[9]提出的基于Nuttall窗的三峰插值谐波算法较好地兼顾了参数估计精度和计算复杂度。与此同时，考虑到FFT方法在工程上简单易行，一些研究者提出了基于FFT的改进谐波参数分析方法，如文献[14-15]提出了基于FFT序列进行变换处理的谐波参数估计算法，在此基础上，文献[16]通过在变换处理之前引入时域加窗处理环节以进一步减小频谱泄露，从而提高后继处理过程中谐波参数估计精度。整体而言，基于FFT的改进算法虽然比前述基于加窗与插值处理相结合的算法其实现复杂度有所降低，但在实际应用过程中，谐波参数估计的精度仍有待进一步提高。本文在分析信号谱序列衰减特征的基础上，提出一种新的基于谱序列变换的高精度谐波参数估计算法。该算法通过对谱序列进行特定的加权变换处理，减少了频谱泄漏造成的谐波间干扰，极大地提高了谐波参数估计的精度。与此同时，由于该方法只需对信号FFT变换后的谱序列直接进行简单的变换处理，因此，具有计算复杂度低的特点。

1 信号频谱分析

设单频率无限长余弦信号为

(1)

其中：，和分别为信号的幅值、频率、相位。信号对应的抽样序列为

(2)

其中：；为抽样频率。将截断成1个长为N点的序列，这相当于将乘上1个长为N点的窗函数：

(3)

注意到对进行FFT变换时相当于预先对加了1个矩形窗，此时，可以表示为

(4)

原信号的复指数序列形式可以表示为

(5)

其中：j为参数单位，j=；为数字域角频率。相应地，可以表示为

(6)

利用频域卷积定理，的离散时间傅里叶变换可表示为

(7)

式中：表示原序列的离散时间傅里叶变换，

(8)

其中：为冲激函数。表示的离散时间傅里叶变换，

(9)

因此，原序列的频谱为

(10)

利用冲激函数的性质可以得到

(11)

式(11)为连续复谱函数，由于该复谱函数的正、负频率点具有对称性，且负频点对谐波参数估计的影响较小，所以，在分析过程中一般忽略负频率成分的影响，而只对信号的正频率部分加以分析和利用。在频域内取N点，以采样间隔对连续谱进行离散化。相应地，在式(11)中令，即得到离散后的谱序列为

(12)

又由于，故谱序列可进一步近似为

(13)

2 基于谱序列变换的谐波分析

2.1　谱序列变换原理

由式(13)可知，当处于同步采样时仅当时，有。此时，谱线的峰值恰好等于准确值，而在其他谱线上有。电网中基波频率处于波动状态使得实际采样总是处于非同步状态，即为一个非整数，故设，其中为的整数部分，为其小数部分，故有

(14)

利用三角函数和差化积公式，式(14)可以表示为

(15)

由式(15)可知，当信号处于非同步采样时，频谱的峰值不等于准确值，而是在其他谱线上产生了分量，即产生了频谱泄漏，并且各次谐波之间将产生相互干扰，从而使得对谐波参数的准确估计十分困难，因此，如何降低频谱泄漏成为高精度谐波参数估计的关键。进一步令

(16)

则式(15)可以简洁地表示为

(17)

式(17)表明信号的幅度谱与成反比。这意味着随着增大，远离真实频率位置的谱线幅值只是以的速度衰减。由于这些谱线的幅值衰减速度不够快，频谱泄露将对谐波参数的估计造成很大干扰，极大地影响谐波参数估计的精度。为了减小频谱泄漏对谐波参数估计精度带来的负面影响，需要尽量加快偏离真实频率位置谱线幅值的衰减速度，为此，本文新设计出一个衰减速度为的新序列：

(18)

需指出的是，尽管可以选择更多的谱线来构造衰减速度更快的新序列，但考虑到算法复杂度和估计精度的平衡，在式(18)所构造的新序列中只是选择了真实频率位置左右各5根谱线参与谐波参数的估计。注意到可以看作1个有理多项式，因此，根据代数理论，新序列可以进一步展开为如下表达式：

(19)

联立式(18)和式(19)可求得

(20)

由式(15)~(17)可得到

(21)

式中：。所以，可以等价地表示为

(22)

可见新的谱序列实际上是将真实频率附近的11个点进行加权处理，由此频谱序列的衰减速度由原来的变为，衰减速度得到极大提高，使得离真实谱线较远处频点的幅值急剧减小，从而有效地抑制频谱泄漏。为验证本文算法对频谱泄漏的抑制效果，现对1个由基波和三次谐波组成的复合信号进行谱分析，信号为

(23)

其中：；。对信号直接进行FFT处理和基于谱序列变换后得到的频谱特征分别如图1和图2所示(其中，频点为频域中离散谱线序号)。对比图1和图2可见：利用FFT算法对信号直接处理得到的谱序列频谱泄漏严重，而采用本文所提方法对谱序列进行变换后，新的谱序列中具有较大幅值的谱线主要集中在峰值谱线附近，减少了谐波之间的干扰，从而达到了进一步抑制频谱泄漏的目的，进而提高了谐波参数估计的精度。

FX_GRP_ID80006236

图1　直接FFT频谱图

Fig. 1　Spectrum graph of direct FFT

FX_GRP_ID8000414C

图2　基于谱序列变换后的频谱图

Fig. 2　Spectrum graph of spectral sequence transformation

2.2　谐波参数估计方法

非同步采样使得谐波真实频率的位置与峰值谱线的位置并不对应，而是位于相邻的2根谱线之间。设真实频率偏离峰值谱线的频率为，若能求得，则可准确计算真实谱线的位置，进而估计相关参数。为此，先在变换后的序列中搜寻幅值最大的2根谱线和。由式(16)和式(18)得：

(24)

(25)

设和的比值为，则

(26)

可以求出偏移量为

(27)

结合式(16)和式(18)不难得到谐波幅值和相位的计算表达式：

(28)

(29)

其中：为的相位。经归纳，谐波参数估计的具体步骤如下。

1) 对序列进行FFT运算，得到谱序列；

2) 由原谱序列根据式(22)计算得到新的谱序列；

3) 在变换后的序列中搜索幅值最大的2根谱线和；

4) 求出和的比值；

5) 通过与之间的关系式(27)计算偏移量；

6) 根据式(28)和(29)计算出幅值和相位。

基于谱序列变换的谐波参数估计流程图如图3所示。

FX_GRP_ID8000478D

图3　谐波参数估计的流程图

Fig. 3　Flow chart of harmonic parameter estimation

上面以单频信号为例介绍了谐波参数估计的方法，该方法可以方便地推广到用于由基波和谐波组成的复合信号的参数估计。具体而言，在复合信号参数估计中，先计算出基波频率，再将第次谐波的频域搜索范围设置为，并在该范围内搜索出幅值最大的2根谱线，根据式(24)~(27)即可依次计算出各次谐波的幅值和相位。

3 仿真与分析

为了验证本文所提算法的有效性，利用本文算法、文献[9]中算法和文献[16]中算法分别对单频信号和复合信号及相关的重要场景进行仿真实验对比。

3.1　单频信号场景

对单频信号作参数估计：

(30)

其中：；；；；采样点数N为64。

利用文献[9]中算法计算的幅值A和相位φ为：

A=24.999 999 879 984 617 V；

φ=1.499 999 859 847 008 rad。

利用文献[16]中算法计算的幅值A和相位φ为：

A=24.999 999 789 256 243 V；

=1.499 999 934 913 475 rad。

利用本文算法计算的幅值A和相位φ为：

A=24.999 999 999 113 207 V；

=1.499 999 999 514 312 rad。

由计算结果可知，利用文献[9]中算法估计幅值和相位的相对误差分别为4.80×10^-9和9.34×10^-8，利用文献[16]中算法估计幅值和相位的相对误差分别为8.43×10^-9和4.33×10^-8，而利用本文算法估计幅值和相位的相对误差分别为3.55×10^-11和3.24×10^-10。显然，对于单频信号而言，本文算法比文献[9]和文献[16]中所提算法具有更高的参数估计精度。

3.2　多频信号场景

设1个由基波和2~21次谐波组成的复合信号，其基波频率为，表达式为

(31)

式中：为基波的幅值；为各次谐波的幅值；为基波的相位；为各次谐波的相位。在仿真过程中，具体参数设置如表1所示。

表1　基波和各次谐波的幅值和相位

Table 1　Amplitude and phase of the fundamental and harmonics

设置采样频率，采样点数N=1 024，对信号分别采用本文算法、文献[9]中算法和文献[16]中算法进行参数估计，这3种算法的幅值和相位估计相对误差对比结果如图4所示。从图4(a)可见：对于基波的幅值估计相对误差，本文算法所得结果比文献[9]中算法所得结果高近7个数量级，比文献[16]中算法结果高3个以上数量级；对于基波以外的谐波，本文算法得到的估计相对误差范围为7.39×10^-11~3.03×10^-14，文献[9]中算法结果相对误差范围为8.74×10^-8~4.82×10^-10，文献[16]中算法结果相对误差范围为7.18×10^-9~6.10×10^-11。从图4(b)可见：对于基波的相位估计相对误差，本文算法所得结果比文献[9]中算法结果高将近4个数量级，比文献[16]中算法结果高近3个数量级；进一步，对于基波以外的各次谐波，本文所提算法的相对误差范围为1.28×10^-9~8.14×10^-13，而文献[9]中算法结果的相对误差范围为2.29×10^-6~2.08×10^-9，文献[16]中算法结果的相对误差范围则为6.21×10^-7~6.34×10^-10。总体来看，对于基波和谐波组成的复合信号，本文算法在幅值和相位参数的估计结果精度与文献[9]和文献[16]中算法结果精度相比显著提高。

FX_GRP_ID80000539

图4　不同算法下谐波幅值和相位估计相对误差对比

Fig. 4　Comparison of harmonic amplitude and phase estimation relative error under different algorithms

在采样点数N为512，1 024和2 048时，幅值、相位估计相对误差仿真结果见图5。从图5可以看出：当采样点数由512增加到1 024时，幅值估计相对误差至少提高了2个数量级，对4次和11次谐波则提高了近5个数量级；相位估计相对误差至少提高了3个数量级，对7次和14次谐波则提高了近4个数量级；当采样点数由512增加到2 048时，幅值估计相对误差至少提高了3个数量级，对2次和15次谐波则提高了近6个数量级；同时，相位估计相对误差至少提高了3个数量级，对5次和9次谐波则提高了近7个数量级。显然，采样点数的增加可以有效地提高算法的估计精度。

FX_GRP_ID8000344F

图5　采样点数N不同时谐波幅值和相位估计相对误差对比

Fig. 5　Comparison of harmonic amplitude and phase estimation relative error when sampling points are different

表2　不同基波频率下幅值估计相对误差

Table 2　Amplitude estimation relative error at different fundamental frequencies

表3　不同基波频率相位估计相对误差

Table 3　Phase estimation relative error at different fundamental frequencies

3.3　基波频率变化场景

电网实际运行过程中的基波频率并不是固定不变的，通常在50 Hz左右波动，因此，考察算法在基波频率变化情况下参数估计性能的稳定性非常必要。GB/T 15945—2008规定电力系统基波频率偏差最大范围为±0.5 Hz，故在仿真过程中基波频率设定在49.5~50.5 Hz^[17^-18]。信号时域表达式如式(31)所示，仿真过程中设置采样频率，采样点数N=1 024。表2和表3所示仿真结果表明：在基波频率变化时，本文所提算法仍然表现出优良的估计性能和良好的稳定性。

3.4　间谐波场景分析

电网在实际运行过程中可能会产生间谐波，其存在会造成电压的闪变并对电网产生危害，因此，在处理谐波过程中，对间谐波参数的高精度估计非常必要^[19^-20]。为了验证本文算法对间谐波参数估计的有效性，在仿真过程中采用1个由7次间谐波组成的复合信号(其基波频率)：

(32)

具体参数设置如表4所示。另外，在仿真过程中设置采样频率，采样点数1 024。采用文献[9]中算法、文献[16]中算法和本文算法分别对上述间谐波信号的幅值和相位参数进行估计，得到各次间谐波的幅值、相位估计的相对误差如表5和表6所示。从仿真结果可以看出：文献[9]和文献[16]中算法在一定程度上提高了间谐波参数估计的精度，但本文算法与文献[9]和[16]中算法相比，估计精度显著提高。

表4　各次间谐波的幅值和相位

Table 4　Amplitude and phase of each inter-harmonic

表5　采用不同算法的时间谐波幅值估计相对误差

Table 5　Amplitude estimation relative error of inter-harmonics under different algorithms

表6　采用不同算法的时间谐波相位估计相对误差

Table 6　Phase estimation relative error of inter-harmonics under different algorithms

4 结论

1）在对谐波信号FFT谱序列的衰减特征进行分析的基础上，通过对原谱序列实施特定加权变换，加快非真实频率幅值的衰减速度，从而达到有效抑制频谱泄漏的目的，并在此基础上给出了一种谐波幅值和相位的高精度估计方法。

2）该算法能有效抑制频谱泄漏。与经典的加窗插值算法和FFT改进算法所得结果相比，该算法估计精度显著提高，且该算法在基波频率变化及间谐波场景下均表现出优良的估计性能。

3）与经典的加窗插值算法相比，本文所提出的新算法无需复杂的加窗处理环节，也无需复杂的插值处理操作，实现过程更简单，计算复杂度低。

参考文献：

[1] MEDINA A, SEGUNDO J, RIBEIRO P, et al. Harmonic analysis in frequency and time domain[J]. IEEE Transactions on Power Delivery, 2013, 28(3): 1813-1821.

[2] SHUAI Zhikang,HUANG Wen,SHEN Chao,et al.Characteristics and restraining method of fast transient inrush fault currents in synchronverters[J]. IEEE Transactions on Industrial Electronics, 2017, 64(9): 7487-7497.

[3] 庞浩, 李东霞, 俎云霄, 等. 应用FFT进行电力系统谐波分析的改进算法[J]. 中国电机工程学报, 2003, 23(6): 50-54.

PANG Hao, LI Dongxia, ZU Yunxiao, et al. An improved algorithm for harmonic analysis of power system using FFT technique[J]. Proceedings of the CSEE, 2003, 23(6): 50-54.

[4] 周西峰, 赵蓉, 郭前岗. Blackman-Harris窗的插值FFT谐波分析与应用[J]. 电测与仪表, 2014, 51(11): 81-85.

ZHOU Xifeng, ZHAO Rong, Guo Qiangang. Blackman-Harris window based interpolation FFT harmonic analysis and its application[J]. Electrical Measurement & Instrumentation, 2014, 51(11): 81-85.

[5] 牛胜锁, 梁志瑞, 张建华, 等. 基于三谱线插值FFT的电力谐波分析算法[J]. 中国电机工程学报, 2012, 32(16): 130-136.

NIU Shengsuo, LIANG Zhirui, ZHANG Jianhua, et al. An algorithm for electrical harmonic analysis based on triple-spectrum-line interpolation FFT[J]. Proceedings of the CSEE, 2012, 32(16): 130-136.

[6] 张俊敏, 刘开培, 汪立, 等. 基于四谱线插值FFT的谐波分析快速算法[J]. 电力系统保护与控制, 2017, 45(1): 139-145.

ZHANG Junmin, LIU Kaipei, WANG Li, et al. A rapid algorithm for harmonic analysis based on four-spectrum-line interpolation FFT[J]. Power System Protection and Control, 2017, 45(1): 139-145.

[7] 华敏, 陈剑云. 一种高精度六谱线插值FFT谐波与间谐波分析方法[J]. 电力系统保护与控制, 2019, 47(11): 9-15.

HUA Min, CHEN Jianyun. A high precision approach for harmonic and interharmonic analysis based on six-spectrum-line interpolation FFT[J]. Power System Protection and Control, 2019, 47(11): 9-15.

[8] 卿柏元, 滕召胜, 高云鹏, 等. 基于Nuttall窗双谱线插值FFT的电力谐波分析方法[J]. 中国电机工程学报, 2008, 28(25): 153-158.

QING Baiyuan, TENG Zhaosheng, GAO Yunpeng, et al. An approach for electrical harmonic analysis based on Nuttall window double-spectrum-line interpolation FFT[J]. Proceedings of the CSEE, 2008, 28(25): 153-158.

[9] 翟瑞淼, 英超, 任国臣, 等. 基于Nuttall窗的三峰插值谐波算法分析[J]. 电力系统保护与控制, 2015, 43(10): 38-43.

ZHAI Ruimiao, YING Chao, REN Guochen, et al. An approach for harmonic analysis based on Nuttall window and triple-spectral-line interpolation[J]. Power System Protection and Control, 2015, 43(10): 38-43.

[10] 李得民, 何怡刚. 基于Nuttall窗四谱线插值FFT的电力谐波分析[J]. 电力系统保护与控制, 2016, 44(3): 64-71.

LI Demin, HE Yigang. Harmonic analysis of power system based on Nuttall window four-spectrum-line interpolation FFT[J]. Power System Protection and Control, 2016, 44(3): 64-71.

[11] 张俊敏, 刘开培, 汪立, 等. 基于乘法窗函数的插值FFT的谐波分析方法[J]. 电力系统保护与控制, 2016, 44(13): 1-5.

ZHANG Junmin, LIU Kaipei, WANG Li, et al. An algorithm for harmonic analysis based on multiplication window function[J]. Power System Protection and Control, 2016, 44(13): 1-5.

[12] WEN He, ZHANG Junhao, MENG Zhuo, et al. Harmonic estimation using symmetrical interpolation FFT based on triangular self-convolution window[J]. IEEE Transactions on Industrial Informatics, 2015, 11(1): 16-26.

[13] 刘冬梅, 杨重良, 何怡刚, 等. 基于Rife-Vincent自卷积窗三谱线插值FFT电力谐波分析[J]. 电子测量与仪器学报, 2016, 30(9): 1351-1360.

LIU Dongmei, YANG Zhongliang, HE Yigang, et al. Harmonic analysis of power system based on Rife-Vincentself-convolution window triple-spectral-line interpolation FFT[J]. Journal of Electronic Measurement and Instrument, 2016, 30(9): 1351-1360.

[14] 薛蕙, 杨仁刚. 基于FFT的高精度谐波检测算法[J]. 中国电机工程学报, 2002, 22(12): 106-110.

XUE Hui, YANG Rengang. Precise algorithms for harmonic analysis based on FFT algorithm[J]. Proceedings of the CSEE, 2002, 22(12): 106-110.

[15] 张红瑛, 滕召胜, 温和, 等. 九点变换改进FFT高精度谐波分析方法[J]. 电力系统及其自动化学报, 2010, 22(6): 38-43, 61.

ZHANG Hongying, TENG Zhaosheng, WEN He, et al. Nine points' polynomial transform improved FFT for high precise harmonic analysis[J]. Proceedings of the CSU-EPSA, 2010, 22(6): 38-43, 61.

[16] 王刘旺, 黄建才, 孙建新, 等. 基于加汉宁窗的FFT高精度谐波检测改进算法[J]. 电力系统保护与控制, 2012, 40(24): 28-33.

WANG Liuwang, HUANG Jiancai, SUN Jianxin, et al. An improved precise algorithm for harmonic analysis based on Hanning-windowed FFT[J]. Power System Protection and Control, 2012, 40(24): 28-33.

[17] 曾博, 滕召胜, 周毅波. 改进FFT非稳态电力谐波分析及应用[J]. 电工技术学报, 2012, 27(11): 256-262.

ZENG Bo, TENG Zhaosheng, ZHOU Yibo. Power system harmonic analysis under nonstationary situations and its application based on improved FFT method[J]. Transactions of China Electrotechnical Society, 2012, 27(11): 256-262.

[18] TAN Tianyuan, CHEN Wenjuan, LIU Kaipei, et al. Harmonic analysis based on time domain mutual-multiplication window[J]. Journal of Modern Power Systems and Clean Energy, 2016, 4(1): 47-53.

[19] 段晨, 杨洪耕. 邻近基波/谐波的间谐波检测算法[J]. 电网技术, 2019, 43(05): 1818-1825.

DUAN Chen, YANG Honggeng. An algorithm for detecting interharmonics adjacent to fundamental/harmonic components[J]. Power System Technology, 2019, 43(05): 1818-1825.

[20] 李宁, 左培丽, 王新刚, 等. 基于改进DFT和时域准同步的间谐波检测算法[J]. 电力自动化设备, 2017, 37(4): 170-178.

LI Ning, ZUO Peili, WANG Xingang, et al. Inter-harmonic detection based on improved DFT and time-domain quasi-synchronization[J]. Electric Power Automation Equipment, 2017, 37(4): 170-178.

(编辑陈灿华)

收稿日期： 2019 -12 -10; 修回日期： 2020 -03 -05

基金项目(Foundation item)：国家自然科学基金资助项目(61861019, 61362018)；湖南省自然科学基金资助项目(2019JJ50483)；湖南省教育厅优秀青年基金资助项目(18B316)；吉首大学研究生科研创新项目(JGY201852)；吉首大学研究生校级科研项目(Jdy19042) (Projects(61861019, 61362018) supported by the National Natural Science Foundation of China; Project(2019JJ50483) supported by the Natural Science Foundation of Hunan Province; Project(18B316) supported by the Department of Education of Hunan Province; Project(JGY201852) supported by the Research Innovation Funds for the Postgraduates of Jishou University; Project(Jdy19042) supported by the School-level Scientific Research Funds for the Postgraduates of Jishou University)

通信作者：雷可君，博士，副教授，从事智能电网、数字信号处理技术与应用等研究；E-mail: leikejun-123@163.com

摘要：基于高精度谐波分析是电力系统谐波污染治理的前提，由于频谱泄漏的影响，经典的谐波分析算法存在参数估计精度不高、计算复杂度较大等问题，为此，在分析信号谱线衰减特征的基础上，提出一种基于谱序列变换的高精度谐波参数分析算法。该算法通过对信号进行FFT运算得到谱序列并对其实施特定的加权变换，加速非真实频率处谱线的衰减速度，从而达到有效抑制频谱泄漏的目的，在此基础上推导谐波幅值和相位的计算公式。研究结果表明，与经典的加窗插值算法和FFT改进算法相比，所提出的算法所得谐波参数估计精度显著提高，且其在基波频率变化及间谐波条件下均表现出优良的估计性能；此外，与经典的加窗插值算法相比，新算法只需对谱序列进行简单处理，具有计算复杂度低的优点。