具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

如图1所示，本发明实施例提供一种提高谐波检测精度的自适应VMD检测方法，包括以下步骤：

1)对谐波信号P(t)进行归一化处理，具体如下：

式中，k为采样点，N为采样点总数。

2)将信号P(t)进行k层变分模态分解，k＝2,3…n+1，求取最优分解层数K_best；具体包括：

将信号P(t)进行k层变分模态分解，k＝2,3…n+1，求取不同分解层数各个本征模态函数IMF的形状因数求取信号P(t)的形状因数SFR_P(t)，定义形状因数MEVAR＝mean(SFR)/var(SFR)，var(SFR)越大，mean(SFR)越小，则MEVAR越小，表征着IMF形状规则，且各个IMF之间相关性较小，IMF得到了有效分离。MEVAR最小值MEVAR_min时的分解层数为谐波信号的最优分解层数K_best。其分解过程以及求解公式如下：

a)初始化分解层数k＝2，定义最大分解层数K_max，初始化参数

b)根据式更新u_k；

c)根据式更新

d)依据更新λ；

e)若满足，则停止迭代，若不满则，则返回步骤b)；

f)计算k层本征模态函数IMF的形状因数以及MEVAR_min，其中

g)判断k是否等于K_max，若等于则程序终止，若不等于则k＝k+1，返回步骤b)，最终得到MEVAR最小值MEVAR_min时的分解层数为最优分解层数K_best。

3)利用最优分解层数K_best对信号P(t)进行分解，得到K_best个本征模态分量c_i(t)，i＝1,2…K_best；具体分解如下：

4)对本征模态分量c_i(t)进行Hilbert变换，得到瞬时幅值a_i(t)和瞬时频率f_i(t)；具体包括：

对本征模态分量c_i(t)进行Hilbert变换，构造解析信号：

其中，表示第i个c_i(t)分量的Hilbert变换，τ表示为时间，表示为c_i(t)分量的瞬时幅值函数，表示为c_i(t)的瞬时频率。

5)从瞬时幅值a_i(t)和瞬时频率f_i(t)中找出最符合信号趋势的波形，对其两端进行处理，端点处理后得到的信号为a_i(t)'和f_i(t)'；本实施例以对瞬时幅值a_i(t)左端进行处理为例对自适应波形匹配算法进行详细阐述，其具体过程如下：

a)寻找出除特征波形外所有三角波形的起始点值a_i(i)，其对应的时间点为：

其中m′_i和n′_i分别为信号a_i(t)的极大值和极小值，分别对应的时间为和a_i(t)的左端点为a_i(1)，a_i(1)-m′₁-n′₁为特征波形，a_i(i)-m′_i-n′_i为最匹配波形，当t_a(i)不在取样点时，采用三次样条插值对其进行处理；

b)求取所有特征波形与三角波形的匹配误差，其误差公式为：

e(i)＝|m′_i-m′₁|+|n′_i-n′₁|+|a_i(i)-a_i(1)|

c)得到最小匹配误差mine(i)，将mine(i)所对应的三角波形作为匹配波形，延拓至a_i(t)的左侧，同理得到a_i(t)信号右侧的匹配波形，端点处理后得到的信号为a_i(t)'和f_i(t)'。

6)求取a_i(t)'和f_i(t)'的均值，得到各谐波分量的幅值及频率等特征量。

下面以一个具体实施例对本发明的方法进行详细阐述：

1)首先构建电网谐波信号的信号模型

P(t)＝sin(ωt)+a_k1sin(3ωt)+a_k2sin(5ωt)++a_k3sin(7ωt)

式中，ω＝2πf₀，f₀＝50Hz，t为采样时长为0.2s，采样频率f_s＝3200Hz，共采样640点，a_k1＝0.3，a_k2＝0.2，a_k2＝0.15，为了更好的模拟实际电网信号以及检验本算法的抗干扰性，在上述信号模型中加入了35dB的白噪声。

2)将信号P(t)进行归一化处理得到信号P_g(t)，对P_g(t)进行VMD分解，其中α＝1500，τ＝0，ε＝1×10^-6，K_max＝15,得到MEVAR_min。如图2所示，MEVAR_min最小时的分解层数为4，其值为2476.9911，图2c)为在最佳分解层数Kbest＝4时VMD分解得到的第一个分量(IMF1)，IMF1为基频分量,第二个分量IMF2为5倍频分量，第三个分量IMF3为3倍频分量，第四个分量IMF4为7倍频分量，由图可知，本发明所提算法自适应选取最优分解层数K，能够有效抑制噪声的干扰，准确提取出各谐波分量。

3)利用Hilbert变换对各个模态分量求取瞬时幅值和瞬时频率，图3a)和图3c)为未经处理的频率和幅值，由图可知，端点数据出现严重失真，对检测效果影响较大。图3b)和图3d)为利用自适应波形匹配算法对端点经过处理后得到的瞬时幅值和频率，由图可知，处理后的信号波形与原信号的自然趋势相符且连接光滑。

4)求取瞬时幅值以及频率的均值，得到谐波分量的幅值及频率信息。由表1可知，自适应波形匹配算法较好的抑制了谐波频率和幅值端点失真的问题，提高谐波分量检测的精确度。

表1

基于同一发明构思，本发明实施例还提供了一种提高谐波检测精度的自适应VMD检测装置，由于该装置所解决问题的原理与前述实施例的方法相似，因此该装置的实施可以参见前述方法的实施，重复之处不再赘述。

如图4所示，为本发明实施例提供的一种提高谐波检测精度的自适应VMD检测装置，可以用于执行上述方法实施例，该装置包括：

归一化处理模块，用于对谐波信号P(t)进行归一化处理；

最优分解层数求取模块，用于将信号P(t)进行k层变分模态分解，k＝2,3…n+1，求取最优分解层数K_best；

变分模态分解模块，用于利用最优分解层数K_best对信号P(t)进行分解，得到K_best个本征模态分量c_i(t)，i＝1,2…K_best；

Hilbert变换模块，用于对本征模态分量c_i(t)进行Hilbert变换，得到瞬时幅值a_i(t)和瞬时频率f_i(t)；

端点处理模块，用于从瞬时幅值a_i(t)和瞬时频率f_i(t)中找出最符合信号趋势的波形，对其两端进行处理，端点处理后得到的信号为a_i(t)'和f_i(t)'；

幅值及频率获取模块，用于求取a_i(t)'和f_i(t)'的均值，得到各谐波分量的幅值及频率。

本发明实施例还提供一种提高谐波检测精度的自适应VMD检测装置，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上任一所述方法的步骤。

本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上任一所述方法的步骤。

综上所述，本发明提供的这种提高谐波检测精度的自适应VMD检测方法、装置及存储介质，基于形状参数均值与方差的比值对VMD进行自适应参数优化选择，选取最优分解层数K，利用自适应波形匹配对谐波信号的幅值、频率信息进行端点处理。通过本发明所做的仿真对比实验表明，其谐波分量幅值误差最低能达到0％，频率误差最低能达到0.49％，比未经处理的所得到的谐波分量幅值精度最大提高了6.67％，频率精度最大提高了1.02％。因此本发明具有良好的噪声鲁棒性，能够自适应选取最优分解层数，还能有效抑制瞬时幅值和频率的端点结果失真的问题，能够准确检测电网谐波信号。

本领域普通技术人员可以理解实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器(ROM，Read Only Memory)、随机存取存储器(RAM，Random AccessMemory)、磁盘或光盘等。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。