一种改进频谱叠加算法的多间谐波识别方法与流程

1.本发明涉及电力电子技术领域，更具体地，涉及一种改进频谱叠加算法的多间谐波识别方法。


背景技术：

2.目前，随着我国大力推进能源改革，风电、光伏等分布式新能源得到了迅速发展，引入大量诸如整流器、逆变器、svg(state voltage generator)、upfc(united power flow controller)等电力电子元件，使得传统电网非线性负载数量增多。系统运行时，非线性负载的电感与回路中的电容引起铁磁谐振，极易在基波及谐波附近产生间谐波。系统中间谐波的含量直接影响了系统及受端设备的稳定运行，增加了滤波装置的复杂程度，在造成大量间谐波损耗的同时也为继电保护的可靠动作带来了困难。此外，间谐波如果造成次同步振荡将会对严重损坏设备，对电力系统造成难以估量的损失。因此，正确识别电力系统信号中间谐波的成分在当下环境中尤为重要。
3.现有技术中，采用单峰谱线插值算法，通过对间谐波附近频谱泄露的幅值进行加权，解出待求间谐波的频率，从而进一步求得间谐波幅值、相位；单峰谱线插值算法能在存在单个间谐波时，准确快速地计算出待求间谐波的频率值，从而得出间谐波的各项参数。实际工程运用中，单间谐波情况并不常见，往往同时存在多个距离相近的间谐波。但当存在多个间谐波时，尤其是当两个间谐波相距较近时，由于存在频谱泄露和主瓣干涉现象，导致离散谱线频点上的频率幅值为各间谐波泄露频谱的叠加，如图1所示，此时，即便采样同步时，传统的单峰谱线插值算法也无法将相互干涉的间谐波分离，也难以得到准确的间谐波频率。
4.现有技术中，还采用传统的加窗插值算法对单个间谐波或多个相距较远的间谐波时进行快速识别，但当间谐波相距较近时，受间谐波频谱泄露、栅栏效应和主瓣干涉的影响识别会出现较大误差。以及采用双峰谱线插值算法对间谐波进行更加精确的识别，但同样不能解决相距较近的多个间谐波的识别问题。此外，还采用adaline方法来简化分析结构、增强跟踪能力，但此方法需要提前知道间谐波的个数频率，具有较大的局限性。还有采用传统小波emd算法来分解信号，以实现对间谐波的精确识别，但无法从根本上解决频谱混叠问题，并且无法准确识别高次频率的间谐波。
5.综上，需要对传统的单峰谱线插值算法进行改进，以满足多间谐波识别的需要。


技术实现要素：

6.为解决现有技术中存在的不足，本发明的目的在于，提供一种改进频谱叠加算法的多间谐波识别方法，结合单峰谱线插值算法与频谱叠加算法，基于韦达定理，提出了多个间谐波存在时的间谐波参数统一求解模型，高效、快速、准确地识别间谐波的频率、幅值、相位参数。
7.本发明采用如下的技术方案。
8.一种改进频谱叠加算法的多间谐波识别方法包括：
9.步骤1，采集电力系统中电压及电流的瞬时信号；
10.步骤2，选择频率k0δf附近多个间谐波的2n个频点；
11.步骤3，基于单峰谱线插值算法和改进频谱叠加算法，利用所采集的瞬时信号，构建所述多个间谐波在2n个频点处的表征参数统一求解模型；
12.步骤4，基于韦达定理，对所述表征参数统一求解模型进行求解，获得各间谐波的频率以及各频率间谐波的幅值和相位。
13.优选地，步骤1中，在10个工频周期内，以5hz为频率分辨率，在基波及谐波处采集电压及电流的幅值、相位的瞬时信号。
14.优选地，步骤2中，频率包括：基波频率和谐波频率；
15.k0表示基波频率或谐波频率转化为频率分辨率的倍数，当k0＝f0/δf时，频率k0δf是基波频率f0；当k0>f0/δf时，频率k0δf是谐波频率k0δf；其中，δf为频率分辨率，f0为基波频率；
16.频率k0δf附近的多个间谐波的频率分别为(k0 i)δf和(k0‑
i)δf，即2n个间谐波的频点；其中i＝1,2,
…
,n，n表示待识别的间谐波的总数。
17.优选地，δf为频率分辨率，满足如下关系式：
18.δf＝f
s
/w
19.式中，f
s
表示采样频率，w表示采样点总数。
20.优选地，步骤3包括：
21.步骤3.1，基于单峰谱线插值算法，构建单个待识别的间谐波在任一个频点处的频谱泄漏模型，满足如下关系式：
[0022][0023]
式中，
[0024]
y
′
(k0‑
i)为第p个待识别的间谐波在间谐波频率(k0‑
i)δf处的频谱泄漏值，
[0025]
y
′
(k0 i)为第p个待识别的间谐波在间谐波频率(k0 i)δf处的频谱泄漏值，
[0026]
其中，i＝1,2,
…
,n，n表示待识别的间谐波的总数，
[0027]
c
p
为第p个待识别的间谐波的幅值相位表征参数，
[0028]
d
p
为第p个待识别的间谐波的频率表征参数，
[0029]
其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；
[0030]
步骤3.2，基于改进频谱叠加算法，构建多个待识别的间谐波在任一个频点处的频谱泄漏模型，满足如下关系式：
[0031][0032]
式中，
[0033]
y(k0‑
i)为n个待识别的间谐波在间谐波频率(k0‑
i)δf处的频谱泄漏值，
[0034]
y(k0 i)为n个待识别的间谐波在间谐波频率(k0 i)δf处的频谱泄漏值，
[0035]
其中，i＝1,2,
…
,n，n表示待识别的间谐波的总数，
[0036]
c1,
…
,c
p
,
…
c
n
分别为n个待识别的间谐波的幅值相位表征参数，
[0037]
d1,
…
,d
p
,
…
d
n
分别为n个待识别的间谐波的频率表征参数，
[0038]
其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；
[0039]
步骤3.3，构建多个待识别的间谐波在2n个频点处的表征参数统一求解模型，满足如下关系式：
[0040][0041]
式中，
[0042]
y(k0‑
n),
…
,y(k0‑
1),y(k0 1),
…
,y(k0 n)，分别为n个待识别的间谐波在2n个频点处的频谱泄漏值，
[0043]
c1,
…
,c
p
,
…
c
n
分别为n个待识别的间谐波的幅值相位表征参数，
[0044]
d1,
…
,d
p
,
…
d
n
分别为n个待识别的间谐波的频率表征参数，
[0045]
其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；
[0046]
任一个离散谱线频点处的谱线幅值，为多个间谐波频谱泄漏幅值的叠加。
[0047]
优选地，步骤3.1包括：
[0048]
步骤3.1.1，基于单峰谱线插值算法，利用瞬时信号经过离散傅里叶变换后得到的信号y(k)，构建单个待识别间谐波的泄露频谱y
′
(k)，满足如下关系式：
[0049][0050]
式中，
[0051]
a
x
为单个待识别间谐波的幅值，
[0052]
为单个待识别间谐波的相位，
[0053]
f
x
为单个待识别间谐波的频率，
[0054]
δf为频率分辨率，
[0055]
k为频谱点，
[0056]
w为采样点总数，
[0057]
其中，定义单个待识别间谐波的幅值相位表征参数，满足如下关系式：
[0058]
[0059]
定义单个待识别间谐波的频率表征参数，满足如下关系式：
[0060]
d
x
＝k
‑
f
x
/δf
[0061]
步骤3.1.2，构建频率为f
x
的间谐波在频率k0δf处的频谱泄露值y
′
(k0)，满足如下关系式：
[0062][0063]
式中，为待识别间谐波的与频率k0f相关的幅值相位表征参数，d
x
为待识别间谐波的频率表征参数；
[0064]
步骤3.1.3，构建频率为f
x
的间谐波分别在频率(k0 1)δf和(k0‑
1)δf处的频谱泄露值y
′
(k0 1)和y
′
(k0‑
1)，满足如下关系式：
[0065][0066]
式中，为待识别间谐波的与频率(k0 1)δf相关的幅值相位表征参数；为待识别间谐波的与频率(k0‑
1)δf相关的幅值相位表征参数；
[0067]
步骤3.1.4，当采样点数w足够大时，频谱点k0的变化不改变幅值相位表征参数，得到单个待识别的间谐波在任一个所述频点处的频谱泄漏模型。
[0068]
进一步，步骤3.1.1中，瞬时信号经过离散傅里叶变换后得到的信号y(k)，包括：
[0069]
以如下关系式表示电力系统中离散采集的瞬时信号y(n)：
[0070][0071]
式中，m为离散采集的信号中含有不同频率的个数，f
s
为采样频率，f
i
为离散采集信号的第i个频率，a
i
为离散采集信号的i个频率对应的幅值，为离散采集信号的i个频率对应的相位，n为采样点，n＝1,2,
…
,w
‑
1，w为采样点总数；
[0072]
以如下关系式对离散采集的瞬时信号y(n)进行离散傅里叶变换：
[0073][0074]
式中，y(k)为离散采集的瞬时信号经过离散傅里叶变换后得到的信号，k为频谱点，k＝1,2,
…
,w
‑
1，w为采样点总数。
[0075]
优选地，步骤4包括：
[0076]
步骤4.1，基于韦达定理，以特定一元高次方程的n个解，作为n个待识别的间谐波的频率表征参数；
[0077]
步骤4.2，利用n个待识别的间谐波的频率表征参数，基于表征参数统一求解模型，求解得到n个待识别的间谐波的幅值相位表征参数；
[0078]
步骤4.3，利用n个待识别的间谐波的频率表征参数和幅值相位表征参数，基于单峰谱线插值算法，得到各间谐波的频率、幅值和相位。
[0079]
进一步，步骤4.3包括：
[0080]
步骤4.3.1，利用n个待识别的间谐波的频率表征参数和幅值相位表征参数，构建多个待识别间谐波的泄露频谱y(k)，满足如下关系式：
[0081][0082]
其中，c表示n个待识别的间谐波的幅值相位表征参数集合，d表示n个待识别的间谐波的频率表征参数集合；
[0083]
步骤4.3.2，将n个待识别的间谐波泄露频谱y(k)分解为单个间谐波频谱y
′
(k)的叠加，满足如下关系式：
[0084][0085]
式中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；
[0086]
步骤4.3.3，对单个待识别的间谐波在2n个频点处的频谱泄漏频谱值进行排序，得到单个间谐波的最大泄露频率幅值|y
′
(k1)|和次最大的泄露频率幅值|y
′
(k2)|；
[0087]
步骤4.3.4，根据单峰谱线插值算法，单个待识别间谐波的频率f
p
满足如下关系式：
[0088][0089]
式中，f
p
为单个待识别间谐波的频率，其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；δf为频率分辨率；k1为最大泄露频率幅值|y
′
(k1)|对应的频谱点；
[0090]
步骤4.3.5，基于多个待识别间谐波的泄露频谱y(k)以及各间谐波频率f
p
，选取n个有关y(k)的方程构成方程组，求解所述方程组得到各间谐波所对应的幅值相位表征参数中的值；
[0091]
步骤4.3.6，对取模长，得到各间谐波对应的幅值，并且中角度为各间谐波对应的相位。
[0092]
进一步，步骤5.3.2中，n个待识别的间谐波泄露频谱y(k)满足如下关系式：
[0093][0094]
式中，
[0095]
a
p
为单个待识别间谐波的幅值，
[0096]
为单个待识别间谐波的相位，
[0097]
f
p
为单个待识别间谐波的频率，
[0098]
δf为频率分辨率，
[0099]
k为频谱点，
[0100]
w为采样点总数。
[0101]
本发明的有益效果在于，与现有技术相比，精确地分离出相距较近的基波及谐波
附近的多个非对称间谐波，解决传统频谱叠加算法对检测间谐波的存在个数和分布限制等问题。同时，本发明模型统一、计算快速、结果精确，能够较好满足实际工程中对间谐波检测的要求。
附图说明
[0102]
图1为现有技术在识别多个间谐波时，各间谐波泄露频谱的叠加示意图；
[0103]
图2为本发明的一种改进频谱叠加算法的多间谐波识别方法的流程框图；
[0104]
图3为本发明一实施例中，改进频谱叠加算法的多间谐波识别方法与传统fft算法，针对频率识别的误差对比图。
具体实施方式
[0105]
下面结合附图对本技术作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本技术的保护范围。
[0106]
如图2，一种改进频谱叠加算法的多间谐波识别方法包括：
[0107]
步骤1，采集电力系统中电压及电流的瞬时信号。
[0108]
具体地，步骤1中，在10个工频周期内，以5hz为频率分辨率，在基波及谐波处采集电压及电流的幅值、相位的瞬时信号。
[0109]
值得注意的是，本发明优选实施例中采用的是iec推荐的10周期采样长度的测量标准，采样频率f
s
为3200hz，采样点数为640点，频率分辨率δf为5hz，是一种非限制性的较优选择，所属领域技术人员可以根据实际应用需要选择不同的测量标准。
[0110]
步骤2，选择频率k0δf附近多个间谐波的2n个频点。
[0111]
具体地，步骤2中，频率包括：基波频率和谐波频率。
[0112]
k0表示基波频率或谐波频率转化为频率分辨率的倍数，当k0＝f0/δf时，频率k0δf是基波频率f0；当k0>f0/δf时，频率k0δf是谐波频率k0δf；其中，δf为频率分辨率，f0为基波频率。
[0113]
频率k0δf附近的多个间谐波的频率分别为(k0 i)δf和(k0‑
i)δf，即2n个间谐波的频点；其中i＝1,2,
…
,n，n表示待识别的间谐波的总数。
[0114]
具体地，δf为频率分辨率，满足如下关系式：
[0115]
δf＝f
s
/w
[0116]
式中，f
s
表示采样频率，w表示采样点总数。
[0117]
步骤3，基于单峰谱线插值算法和改进频谱叠加算法，利用所采集的瞬时信号，构建所述多个间谐波在2n个频点处的表征参数统一求解模型。
[0118]
具体地，步骤3包括：
[0119]
步骤3.1，基于单峰谱线插值算法，构建单个待识别的间谐波在任一个频点处的频谱泄漏模型。
[0120]
基频同步采样前提下，基波频率或者谐波频率仅在自身频率处存在频谱值，不会产生频谱泄露现象。此时，基波频率或谐波频率的频谱值为y
′
(k0)＝y(k0) a，其中a为基波/谐波自身频率的频谱，k0δf为基波频率或者谐波频率。
[0121]
由于基波频率或者谐波频率附近的频率y(k0 i)、y(k0‑
i)不受基波频率或者谐波
频率的频谱泄露的影响，间谐波f
x
在各附近频率处频谱值，满足如下关系式：
[0122][0123]
式中，
[0124]
y
′
(k0‑
i)为第p个待识别的间谐波在间谐波频率(k0‑
i)δf处的频谱泄漏值，
[0125]
y
′
(k0 i)为第p个待识别的间谐波在间谐波频率(k0 i)δf处的频谱泄漏值，
[0126]
其中，i＝1,2,
…
,n，n表示待识别的间谐波的总数，
[0127]
c
p
为第p个待识别的间谐波的幅值相位表征参数，
[0128]
d
p
为第p个待识别的间谐波的频率表征参数，
[0129]
其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数。
[0130]
具体地，步骤3.1包括：
[0131]
步骤3.1.1，基于单峰谱线插值算法，利用瞬时信号经过离散傅里叶变换后得到的信号y(k)，构建单个待识别间谐波的泄露频谱y
′
(k)，满足如下关系式：
[0132][0133]
式中，
[0134]
a
x
为单个待识别间谐波的幅值，
[0135]
为单个待识别间谐波的相位，
[0136]
f
x
为单个待识别间谐波的频率，
[0137]
δf为频率分辨率，
[0138]
k为频谱点，
[0139]
w为采样点总数，
[0140]
其中，定义单个待识别间谐波的幅值相位表征参数，满足如下关系式：
[0141][0142]
定义单个待识别间谐波的频率表征参数，满足如下关系式：
[0143]
d
x
＝k
‑
f
x
/δf
[0144]
步骤3.1.2，构建频率为f
x
的间谐波在频率k0δf处的频谱泄露值y
′
(k0)，满足如下关系式：
[0145][0146]
式中，为待识别间谐波的与频率k0f相关的幅值相位表征参数，d
x
为待识别间谐波的频率表征参数；
[0147]
步骤3.1.3，构建频率为f
x
的间谐波分别在频率(k0 1)δf和(k0‑
1)δf处的频谱泄露值y
′
(k0 1)和y
′
(k0‑
1)，满足如下关系式：
[0148][0149]
式中，为待识别间谐波的与频率(k0 1)δf相关的幅值相位表征参数；为待识别间谐波的与频率(k0‑
1)δf相关的幅值相位表征参数；
[0150]
步骤3.1.4，当采样点数w足够大时，频谱点k0的变化不改变幅值相位表征参数，即因此得到单个待识别的间谐波在任一个所述频点处的频谱泄漏模型。
[0151]
进一步，步骤3.1.1中，瞬时信号经过离散傅里叶变换后得到的信号y(k)，包括：
[0152]
以如下关系式表示电力系统中离散采集的瞬时信号y(n)：
[0153][0154]
式中，m为离散采集的信号中含有不同频率的个数，f
s
为采样频率，f
i
为离散采集信号的第i个频率，a
i
为离散采集信号的i个频率对应的幅值，为离散采集信号的i个频率对应的相位，n为采样点，n＝1,2,
…
,w
‑
1，w为采样点总数；
[0155]
以如下关系式对离散采集的瞬时信号y(n)进行离散傅里叶变换：
[0156][0157]
式中，y(k)为离散采集的瞬时信号经过离散傅里叶变换后得到的信号，k为频谱点，k＝1,2,
…
,w
‑
1，w为采样点总数。
[0158]
考虑到实际采样过程中，对谐波进行同步采样的难度很大，因此所求频率f
x
通常不落在离散谱线频点上。
[0159]
步骤3.2，基于改进频谱叠加算法，构建多个待识别的间谐波在任一个频点处的频谱泄漏模型，满足如下关系式：
[0160][0161]
式中，
[0162]
y(k0‑
i)为n个待识别的间谐波在间谐波频率(k0‑
i)δf处的频谱泄漏值，
[0163]
y(k0 i)为n个待识别的间谐波在间谐波频率(k0 i)δf处的频谱泄漏值，
[0164]
其中，i＝1,2,
…
,n，n表示待识别的间谐波的总数，
[0165]
c1,
…
,c
p
,
…
c
n
分别为n个待识别的间谐波的幅值相位表征参数，
[0166]
d1,
…
,d
p
,
…
d
n
分别为n个待识别的间谐波的频率表征参数，
[0167]
其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数。
[0168]
实际工程运用中，单个间谐波情况并不常见，往往同时存在多个距离相近的间谐波。传统频谱叠加算法在计算时存在很多限制，通常叠加的上限为两个的间谐波对。本发明优选实施例中提出了适用于两个以上简谐波对的频谱叠加算法。
[0169]
步骤3.3，构建多个待识别的间谐波在2n个频点处的表征参数统一求解模型，满足如下关系式：
[0170][0171]
式中，
[0172]
y(k0‑
n),
…
,y(k0‑
1),y(k0 1),
…
,y(k0 n)，分别为n个待识别的间谐波在2n个频点处的频谱泄漏值，
[0173]
c1,
…
,c
p
,
…
c
n
分别为n个待识别的间谐波的幅值相位表征参数，
[0174]
d1,
…
,d
p
,
…
d
n
分别为n个待识别的间谐波的频率表征参数，
[0175]
其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；
[0176]
任一个离散谱线频点处的谱线幅值，为多个间谐波频谱泄漏幅值的叠加。
[0177]
构建求解多个间谐波的表征参数统一求解模型，在工程中具有较高的应用价值，通过求解该模型，可解决对间谐波种类及数量的限制问题。
[0178]
步骤4，基于韦达定理，对所述表征参数统一求解模型进行求解，获得各间谐波的频率以及各频率间谐波的幅值和相位。
[0179]
具体地，步骤4包括：
[0180]
步骤4.1，基于韦达定理，以特定一元高次方程的n个解，作为n个待识别的间谐波的频率表征参数。
[0181]
将表征参数统一求解模型简化为方程组，满足如下关系式：
[0182]
a
·
b＝c
[0183]
其中，矩阵a的表达式如下：
[0184][0185]
其中，矩阵b的表达式如下：
[0186][0187]
其中，矩阵c的表达如下：
[0188]
c＝[(
‑
k)
n
y(k)]
[0189]
在矩阵a、矩阵b、矩阵c的表达式中，y(k)为频谱点k处的频谱泄漏值。
[0190]
通过一元高次方程的韦达定理可知，对于任意一个n次方程∑a
i
x
i
，它的全部n个根x
i
，满足如下关系式：
[0191][0192]
上式所示方程组的左侧与矩阵b的表达式吻合，为简化方程，在此假设a(0)＝1，可以得出该特定一元n次方程的n个系数。
[0193]
该特定方程得出的n个解为上述矩阵b中n个参数d的值，再代入矩阵b的剩余式中，即可求解关于c的线性方程组。
[0194]
步骤4.2，利用n个待识别的间谐波的频率表征参数，基于表征参数统一求解模型，求解得到n个待识别的间谐波的幅值相位表征参数。
[0195]
步骤4.3，利用n个待识别的间谐波的频率表征参数和幅值相位表征参数，基于单峰谱线插值算法，得到各间谐波的频率、幅值和相位。
[0196]
进一步，步骤4.3包括：
[0197]
步骤4.3.1，利用n个待识别的间谐波的频率表征参数和幅值相位表征参数，构建多个待识别间谐波的泄露频谱y(k)，满足如下关系式：
[0198][0199]
其中，c表示n个待识别的间谐波的幅值相位表征参数集合，d表示n个待识别的间谐波的频率表征参数集合；
[0200]
步骤4.3.2，将n个待识别的间谐波泄露频谱y(k)分解为单个间谐波频谱y
′
(k)的叠加，满足如下关系式：
[0201]
[0202]
式中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；
[0203]
步骤4.3.3，对单个待识别的间谐波在2n个频点处的频谱泄漏频谱值进行排序，得到单个间谐波的最大泄露频率幅值|y
′
(k1)|和次最大的泄露频率幅值|y
′
(k2)|；
[0204]
步骤4.3.4，根据单峰谱线插值算法，单个待识别间谐波的频率f
p
满足如下关系式：
[0205][0206]
式中，f
p
为单个待识别间谐波的频率，其中，p＝1,2,
…
,n，n表示待识别的间谐波的总数；δf为频率分辨率；k1为最大泄露频率幅值|y
′
(k1)|对应的频谱点；
[0207]
步骤4.3.5，基于多个待识别间谐波的泄露频谱y(k)以及各间谐波频率f
p
，选取n个有关y(k)的方程构成方程组，求解所述方程组得到各间谐波所对应的幅值相位表征参数中的值；
[0208]
步骤4.3.6，对取模长，得到各间谐波对应的幅值，并且中角度为各间谐波对应的相位。
[0209]
进一步，步骤4.3.2中，n个待识别的间谐波泄露频谱y(k)，满足如下关系式：
[0210][0211]
式中，
[0212]
a
p
为单个待识别间谐波的幅值，
[0213]
为单个待识别间谐波的相位，
[0214]
f
p
为单个待识别间谐波的频率，
[0215]
δf为频率分辨率，
[0216]
k为频谱点，
[0217]
w为采样点总数。
[0218]
实施例1。
[0219]
对于现实情况中多组非对称间谐波共存的电信号，采用iec推荐的10周期采样长度的测量标准，采样频率f
s
为3200hz，采样点数为640点，频率分辨率δf为5hz，信号中存在一组间谐波对，各成分含量如表1所示。
[0220]
表1存在一组间谐波对的实验信号成分含量表
[0221]
频率/hz幅值/％相位/
°
501006045.734554.3345
[0222]
表1中的实验信号为基波和分布在其附近的一对间谐波，此类分布为工程中常见情况，具有实际工程意义。
[0223]
针对包含一组间谐波对的电信号采用传统的fft算法进行识别，其识别内容从频率识别、幅值识别和相位识别三个维度展开，结果如表2所示。
[0224]
表2存在一对间谐波对的实验信号传统fft分析结果
[0225][0226][0227]
根据实验结果可以得知，传统fft对基波的识别能力明显比对间谐波识别能力强，其对基波频率识别误差控制在0.4％，对基波幅值识别误差控制在0.428％，而对基波相位识别的误差较大，相差接近一倍。此外，传统fft对间谐波的识别结果收效甚微，其频率识别的误差分别飙升为3.98％和2.89％，幅值识别的误差之大较难应用于工程实践，分别为30.6％和30.8％，相位识别误差更大。
[0228]
采用本发明所提出的间谐波算法对原模拟信号进行识别，得到识别信号的间谐波的频率、幅值、相位参数，同时结合原信号进行误差分析，得出的实验结果如表3所示。
[0229]
表3存在一对间谐波对的实验信号本文算法分析结果及其误差
[0230]
频率/hz频率误差幅值/％幅值误差相位/
°
相位误差1000100.011.00
×
10
‑459.823.00
×
10
‑345.7003.0258.33
×
10
‑344.882.67
×
10
‑354.3003.000044.992.22
×
10
‑4[0231]
由表3可知，本发明提出的算法能对基波及谐波附近的一对间谐波对进行准确识别，对于一对间谐波对的频率估计几乎零误差，幅值和相位误差都可降低至0.1％数量级。实验结果表明对于工程中常见的存在一对间谐波对的情况，本发明提出的算法能进行高精度识别。
[0232]
为进一步验证本算法的可推广性，本发明优选实施例还设置了基波/谐波附近存在多个间谐波的场景，该信号各成分含量如表4所示。
[0233]
表4存在四个间谐波的实验信号成分含量表
[0234]
频率/hz幅值/％相位/
°
1001006085.533097.92.560104.53120115.75120
[0235]
实验中的信号为二次谐波和分布在其附近的多个临近间谐波，频率、幅值、相位设置参照实际电弧炉信号进行设置，具有现实意义。
[0236]
针对包含四个间谐波的电信号采用传统的fft算法进行识别，其识别内容从频率识别、幅值识别和相位识别三个维度展开，结果如表5所示。
[0237]
表5存在四个间谐波对的实验信号传统fft分析结果
[0238]
频率/hz幅值/％相位/
°
96.6898.81959.6979.12.508100.0290.822.001172.8999.612.193
‑
146.64111.334.723122.15
[0239]
根据实验结果可以得知，相比前一组实验，传统fft对基波频率的识别误差已经从0.4％增加至3.32％，其相对误差量飙升8倍。从幅值的维度观察，其识别误差从0.428％增加至1.181％，相对误差为原来的2.76倍。因此，可以推断当间谐波信号的个数增加后，传统的fft算法对基波的识别能力也有了显著的下降。
[0240]
此外，当间谐波信号个数增加后，传统fft算法对间谐波的识别能力数据支撑如下。传统fft对间谐波的识别结果误差依然很大，其频率识别的误差分别为7.48％、7.23％、4.68％和3.78％，相比前一组实验的3.98％和2.89％，依然没有得到改善。其幅值识别误差分布在5％
‑
27％不等，但都在5％以上，与此同时，其相位识别误差分布在1.79％到222.2％不等，波动较大，难以应用于工程实践。可以推断当间谐波信号的个数增加后，传统的fft算法对间谐波的识别能力也并没有得到改善。
[0241]
综上，无论是在仅有一组间谐波对的情况下，还是存在多个间谐波时，传统的fft对基波和间谐波的识别能力都十分有限。
[0242]
实验采用本发明所提出的宽频谐波识别算法对包含多间谐波的原模拟电信号进行识别，得到识别信号的间谐波的频率、幅值、相位参数，同时结合原信号进行误差分析，得出的实验结果如表6所示。
[0243]
表6存在四个间谐波的实验信号成分含量表
[0244][0245][0246]
其中本发明提出的算法与fft算法针对频率识别的误差对比图如图3所示。
[0247]
根据图3可知，针对当下电力系统中广泛存在的含有间谐波的电信号，本发明提出的算法的频率识别误差远远小于传统fft算法的频率识别误差。
[0248]
此外，由表6可知，本发明提出的算法能对基波及谐波附近的间谐波进行准确识别，其准确度表现不仅仅局限于频率识别，还包括幅值识别，相位识别等不同维度。对于四个间谐波频率误差可降低至0.1％数量级，幅值和相位误差都可降低至1％数量级，相比传统fft得到显著改善。
[0249]
与此同时，相比前一组实验只有一组间谐波对，本实验对于包含多个间谐波的电信号进行三个维度的准确识别也证明了本专利提出的统一模型和算法具有较好的普适性，
能在工程实际中推广。
[0250]
本发明的有益效果在于，与现有技术相比，精确地分离出相距较近的基波及谐波附近的多个非对称间谐波，解决传统频谱叠加算法对检测间谐波的存在个数和分布限制等问题。同时，本发明模型统一、计算快速、结果精确，能够较好满足实际工程中对间谐波检测的要求。
[0251]
本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述，但是本领域技术人员应该理解，以上实施示例仅为本发明的优选实施方案，详尽的说明只是为了帮助读者更好地理解本发明精神，而并非对本发明保护范围的限制，相反，任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。