CN115062668A - 基于RAdam优化宽度学习的谐波参数检测方法和系统 - Google Patents

基于RAdam优化宽度学习的谐波参数检测方法和系统 Download PDF

Info

Publication number
CN115062668A
CN115062668A CN202210778084.6A CN202210778084A CN115062668A CN 115062668 A CN115062668 A CN 115062668A CN 202210778084 A CN202210778084 A CN 202210778084A CN 115062668 A CN115062668 A CN 115062668A
Authority
CN
China
Prior art keywords
harmonic
parameter detection
parameters
detection model
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210778084.6A
Other languages
English (en)
Inventor
张莹莹
李聪聪
李思齐
孙晨宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202210778084.6A priority Critical patent/CN115062668A/zh
Publication of CN115062668A publication Critical patent/CN115062668A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00002Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by monitoring

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于RAdam优化宽度学习的谐波参数检测方法和系统,该方法包括:获取谐波信号及其谐波参数并构建原始谐波信号数据集,再按比例划分为训练集和测试集;利用训练集训练出谐波参数检测模型;其中,包括初始化设置各项参数,后续再利用RAdam优化器自适应更新宽度学习系统的各项参数直至满足迭代终止条件;最后再输入测试集数据,并利用训练好的谐波参数检测模型估计谐波参数并计算估计误差。通过上述构建的谐波参数检测模型可以估算出电力系统中待检测谐波信号的谐波参数。综上,本发明所述方法应用在由于大量非线性设备的接入而使电流信号产生畸变的电力系统中,可以高效以及准确的估算出谐波信号的谐波参数。

Description

基于RAdam优化宽度学习的谐波参数检测方法和系统
技术领域
本发明属于电力系统中的谐波信号智能检测技术,具体涉及一种基于RAdam优化宽度学习的谐波参数检测方法和系统。
背景技术
在现代电力系统迅速发展和能源体系迭代升级的背景下,非线性电力电子设备应用广泛。大量非线性设备的接入可能导致电力系统中电流信号畸变而形成谐波,谐波的存在会降低公共电网中发电和输配电效率,影响用电设备的正常运转且增加额外的损耗,同时还隐藏着巨大的安全隐患。精确、快速的谐波检测和分析方法可以有效地为谐波的抑制和补偿奠定基础,是解决电力系统中非现性设备接入带来的电流信号畸变问题的前提。
谐波检测的经典方法如快速傅里叶变换(FFT)、小波包变换(WPT)等,不过需要完整周期的信号进行检测和分析是这类方法的特点,当采用较少周期数时,FFT和WPT一般只对基波检测比较有效,当存在基频偏差、间谐波和噪声等复杂多干扰因素时,其对高次谐波的检测误差将进一步增加,一般需要10/12个周期才能达到期望的精度。而电力系统中由于非线性设备的接入将导致电力系统内存在基波偏差、间谐波和噪声等复杂多干扰因素,对存在上述多干扰因素的谐波信号进行估计时,谐波次数越高,信号的幅值、频率、相位等参数越容易受干扰因素的影响,因此,经典的FFT以及WPT应用于电力系统的谐波估计将导致较大的检测误差,从而需要更多个周期信号才能达到期望精度。
近年来,神经网络得益于其强大的非线性关系拟合能力,在图像处理、目标检测、故障诊断等领域取得了蓬勃发展和应用。目标应用神经网络进行谐波参数估计技术主要有4种,依次为基于多层感知器神经网络(MLPNN)、径向基网络(RBFNN)、小波神经网络(AWNN)、回声状态网络(ESN)的谐波参数估计方法,然而,MLPNN的多层结构和基于BP算法参数学习机制本身的贪婪特性容易导致算法收敛速度慢或陷入局部最优;基于RBFNN和基于AWNN的方法均采用了梯度下降法参数学习机制,该机制的学习率对方法的估计精度影响较大且需要人工调节,太小容易收敛速度慢,太大容易发散或溢出,基于ESN方法采用元启发式群智能优化算法更新参数,因为需要不断地进行参数寻优过程,在提高精度的同时也具有较高的时间复杂度,且参数最优解的形成依赖于优化算法的初始值,包括初始种群数等。
因此,本发明致力于探索另外一种神经网络,将其应用于电力系统的谐波参数估计。其中,宽度学习系统(BLS)是在随机向量函数链接神经网络(RVFLNN)的基础上被提出,不同其他神经网络通过增加层数提高网络响应精度,宽度学习系统(BLS)通过增加网络节点提升拟合准确率,其具有拓扑结构简单、训练速度快等特点,且网络结构改变后无需重新训练整个网络,采用增量学习算法训练新增节点即可。为此,本发明将探索如何将BLS网络应用于电力系统的谐波参数估计,提供一种全新的谐波参数检测技术,充分发挥以及利用BLS网络优势。
发明内容
本发明为了解决电力系统因非线性设备接入带来的谐波检测的问题,提供了一种全新的技术思路来实现谐波参数检测。具体是提供了一种基于RAdam优化宽度学习的谐波参数检测方法和系统,该谐波参数检测方法将BLS网络引入到电力系统中的谐波信号智能检测技术内,更是进一步将RAdam与BLS网络结合,其中,利用RAdam优化BLS网络的参数,使得基于RAdam与BLS网络构建的谐波参数检测模型的精度更高,可以用于检测谐波信号的幅值、频率、相位等谐波参数。
一方面,本发明提供了一种基于RAdam优化宽度学习的谐波参数检测方法,包括以下步骤:
步骤S1:获取电力系统中的谐波信号以及所述谐波信号对应的谐波参数,进而构建出原始谐波信号数据集;
步骤S2:利用所述原始谐波信号数据集训练出谐波参数检测模型,其中,所述谐波检测模型是基于宽度学习系统构建的,并采用RAdam优化模型参数;
其中,所述谐波参数检测模型的模型输入为谐波信号,模型输出为谐波信号对应的谐波参数;
步骤S3:获取待检测的电力系统中的谐波信号并输入至所述谐波参数检测模型得到谐波信号的谐波参数。
在电力系统中,由于大量非线性负载的接入,导致电力系统的谐波信号相比于其他领域的谐波畸变程度严重得多,因此估计谐波参数的难度也更高;同时电力系统的谐波畸变具有随机性和波动性,运行时谐波状态随时可能发生变化,因此,对电力系统的谐波信号的谐波参数检测,实时性以及准确性是主要考虑点。其中,BLS网络结构简单且训练速度快,同时其本身的增量学习网络更新机制使得在谐波状态发生变化需要保证模型精度且不用重新训练整个网络,因此面对谐波状态突发变化等情况,BLS可快速适应谐波实时状态并准确预测谐波参数,更加符合电力系统的谐波检测需求。为此,本发明技术方案创造性地将BLS网络引入到电力系统的谐波信号检测中,进而提供一种基于RAdam优化宽度学习的谐波参数检测方法,实验结果表明,基于RAdam优化的宽度学习谐波参数检测方法的估计精度高,同时运算速度能够满足谐波实时检测和分析需求。
其次,本发明在BLS网络基础上加入了RAdam优化器。其中,RAdam是Adam的改进算法,其克服了Adam在优化过程中因方差太大而导致的局部最优,其收敛性和稳定性也更佳。再者,RAdam在Adam的基础上加入了方差动态减少机制—“整流项”,克服了Adam在训练初期因数据有限而造成方差过大,也避免了预热阶段手动调整学习率,这种机制保证了RAdam优化宽度学习的谐波参数估计方法在运行过程中不会像Adam等优化器那样可能陷入待优化参数的局部最优解,且能够自适应地调整学习率,因此能够与快速、稳定地收敛至谐波参数估计所需的BLS各层全局最优参数,最终能够准确估计出电力系统的各项谐波参数,更加符合电力系统的谐波参数检测需求。
此外,和经典的FFT、WPT等方法相比,针对
Figure BDA0003719224430000031
周期或
Figure BDA0003719224430000032
周期的谐波信号,应用本发明技术方案提供的谐波参数检测方法也可以实现谐波参数估计;整体而言,本发明技术方案提供的所述谐波参数检测方法具有拓扑结构简单、估计精度高和计算时间短的特点,适用于在大量非线性设备接入的电力系统中,在存在基波偏差、间谐波和噪声等复杂多干扰因素的情况下对谐波参数进行实时高效估计。
进一步可选地,步骤S2中所述谐波参数检测模型的训练过程如下:
步骤S2-1:初始化设置:其中,包括设置所述谐波参数检测模型中宽度学习系统内特征映射层-增强层、增强层-输出层的权值和偏置;
步骤S2-2:将所述原始谐波信号数据集中的谐波信号输入所述宽度学习系统得到映射特征与增强特征之间的联合特征;
步骤S2-3:判断是否需要初始化所述宽度学习系统中输出层的权值;
其中,若是第一次进入训练所述谐波参数检测模型,则利用W=A+Y得到输出层的权值,进入下一次迭代训练,返回步骤S2-2;否则,执行步骤S2-4;W表示输出层的权值,Y表示期望的网络输出,并选用所述原始谐波信号数据集中谐波信号的谐波参数表示;A+表示对联合特征A进行伪逆运算;
步骤S2-4:采用RAdam优化器更新所述宽度学习系统的参数,其中,更新的参数包括:特征映射层-增强层、增强层-输出层的权值和偏置以及输出层的权值;
步骤S2-5:判断是否达到迭代终止条件,若达到,训练结束,得到所述谐波参数检测模型的最优模型参数;否则,返回步骤S2-2进入下一次的迭代。
进一步可选地,所述RAdam优化器的更新策略为:
Figure BDA0003719224430000033
其中,t为时间步数,即参数迭代更新的次数;θtt-1表示更新后、更新前的谐波参数检测模型的模型参数,αt为步长,kt为方差修正项,
Figure BDA0003719224430000034
为偏差修正移动平均线,lt为自适应学习率。
进一步可选地,所述方差修正项,偏差修正移动平均线以及自适应学习率的计算如下:
γmax=2/(1-μ2)-1
Figure BDA0003719224430000041
当γt>γ时,kt
Figure BDA0003719224430000042
和lt满足:
Figure BDA0003719224430000043
否则,取kt=lt=1以及
Figure BDA0003719224430000044
其中,γmax、γt、γ分别表示第t步的近似SMA的最大值、近似SMA值、近似SMA的阈值;μ1为一阶矩估计指数衰减率,μ2为二阶矩估计指数衰减率,
Figure BDA0003719224430000045
为μ2的t次方,mt为有偏一阶矩估计,nt为有偏二阶矩估计,其计算公式如下:
mt=μ1mt-1+(1-μ1)gt
Figure BDA0003719224430000046
其中,gt为待更新参数θ的梯度,其满足:
Figure BDA0003719224430000047
E为目标函数,其与谐波参数的模型输出值以及真实值有关。
本发明的技术方案,其基于参数的梯度信息,建立基于RAdam优化器的参数学习策略,使得待优化各参数的学习率根据预测目标进行自适应调整;此外,RAdam优化器的收敛速度快、收敛精度高,因此可进一步提高模型结果的准确率和稳定性。
进一步可选地,所述宽度学习系统内的增强函数为双曲正切函数,特征映射函数采用线性函数。
进一步可选地,所述原始谐波信号数据集中谐波信号对应的独立样本是将整周期的谐波信号分解为1/2周期或1/4周期得到的。
第二方面,本发明提供一种基于所述谐波参数检测方法的模拟系统,其包括相互连接的谐波信号采集子系统以及谐波信号处理子系统;
其中,所述谐波信号采集子系统用于采集电力系统的谐波信号;所述谐波信号处理子系统用于根据步骤S2方式进行处理谐波信号进而得到谐波参数检测模型,其中,将待检测的电力系统中的谐波信号输入至所述谐波参数检测模型得到谐波信号的谐波参数;
所述谐波信号采集子系统至少包括:输入交流电源、非线性电力电子设备、可编程大功率直流电子负载、以及信号采集单元以及信号存储单元;
其中,输入交流电源与非线性电力电子设备连接,所述输入交流电源向所述谐波信号采集子系统提供可变的输出电压、输出频率和负载电流于非线性电力电子设备,进而产生谐波;
所述可编程大功率直流电子负载用于模拟实际负载且大小可调;
所述信号采集单元用于采集畸变的谐波信号,并与所述信号存储单元相连,将采样后的信号数据传输到信号存储单元进行存储,以便谐波信号处理子系统进行后续分析。
第三方面,本发明提供一种基于所述谐波参数检测方法的系统,其包括:
原始谐波信号数据集构建模块,用于获取电力系统中的谐波信号以及所述谐波信号对应的谐波参数,进而构建出原始谐波信号数据集;
谐波参数检测模型构建模块,用于利用所述原始谐波信号数据集训练出谐波参数检测模型,其中,所述谐波检测模型是基于宽度学习系统构建的,并采用RAdam优化模型参数;
其中,所述谐波参数检测模型的模型输入为谐波信号,模型输出为谐波信号对应的谐波参数;
检测模块,用于将待检测的电力系统中的谐波信号输入至所述谐波参数检测模型得到谐波信号的谐波参数。
第四方面,本发明提供一种电子终端其包括:
一个或多个处理器;
存储了一个或多个计算机程序的存储器;
所述处理器调用所述计算机程序以实现:
一种基于RAdam优化宽度学习的谐波参数检测方法的步骤。
第五方面,本发明提供一种可读存储介质,其存储了计算机程序,所述计算机程序被处理器调用以实现:
一种基于RAdam优化宽度学习的谐波参数检测方法的步骤。
有益效果
本发明技术方案提供的谐波参数检测方法,将拓扑结构简单、运算速度快的宽度学习系统(BLS)应用于电力系统中谐波信号的参数估计,充分利用BLS网络结构简单、训练速度快以及增量学习网络更新机制,能够快速适应谐波实时状态并准确预测谐波参数。此外,本发明所述谐波参数检测方法还在BLS网络基础上加入了RAdam优化器,克服了Adam在训练初期因数据有限而造成方差过大,也避免了预热阶段手动调整学习率,保证了RAdam优化的宽度学习谐波参数估计方法运行过程中不会像Adam等优化器那样可能陷入待优化参数的局部最优解,且能够自适应地调整学习率,最终能够与快速、稳定地收敛至谐波参数估计所需的BLS各层全局最优参数,进而准确估计出电力系统的各项谐波参数。综上,本发明技术方案提供的谐波参数检测方法能够准确快速估算出各项谐波参数,同时,针对
Figure BDA0003719224430000061
周期甚至
Figure BDA0003719224430000062
周期的谐波信号,也能准确估算,更符合电力系统的谐波参数检测需求。
附图说明
图1是本发明提供的基于RAdam优化宽度学习的谐波参数检测方法的流程图;
图2是本发明提供的谐波参数检测模型的训练流程图;
图3是本发明提供的谐波参数检测方法的模拟系统的结构框图;
图4是本发明实例中采用的谐波信号波形;
图5是本发明实例中12个周期的谐波信号FFT分析结果;
图6是本发明实例中1—15次奇次谐波幅值参数;
图7是本发明实例中基于RAdam优化的宽度学习系统总均方误差(平均值)变化过程;
图8是本发明实例中测试集相对误差最大值;
图9是本发明实例中测试集相对误差平均值。
具体实施方式
本发明技术方案提供的一种基于RAdam优化宽度学习的谐波参数检测方法适用于电力系统中谐波信号的参数估计。其核心是将BLS网络引入到电力系统中的谐波信号智能检测技术内,提供了一种全新的技术思路来实现谐波参数检测;其次是,在BLS网络的基础上,引入RAdam,利用RAdam优化BLS网络的参数,使得基于RAdam与BLS网络构建的谐波参数检测模型的精度更高。具体的,本发明从谐波参数估计的实时性和精确性角度考虑,从而选择BLS网络以及RAdam。一方面,在实际的谐波估计场景中,基频偏差、间谐波和噪声等多种干扰因素的存在给电力系统的谐波参数估计模型带来了挑战,这也是经典的方法如快速傅里叶变换(FFT)、小波包变换(WPT)只对基波估计比较有效的原因之一,而宽度学习系统(BLS)具有强大的非线性拟合功能,对所有次数的谐波估计都有效,进一步采用迭代速度快、寻优能力强的RAadm优化器,进而建立基于参数梯度信息的参数自适应学习策略,可以明显提升模型的估计精度。二方面,宽度学习系统(BLS)具有拓扑结构简单、训练周期短和运算速度快的特点,RAadm优化器的强大寻优能力使得模型能够快速收敛,满足谐波估计的实时性需求。因此,本发明的谐波参数检测方法可以准确、高效地估算出电力系统内谐波信号的谐波参数。
应当理解,本发明技术方案提供的谐波参数检测方法可以实现幅值、频率、相位等谐波参数的估计,其中,幅值参数受容易系统中非线性负载的影响且畸变较大,因此下述实施例将以包含幅值为模型响应为例进行详细陈述,下面将结合实施例对本发明做进一步的说明。
实施例1:
本实施例提供一种基于RAdam优化宽度学习的谐波参数检测方法,包括以下步骤:
步骤1:获取电力系统的谐波信号,进而构建出原始谐波信号数据集,并按比例将所述原始谐波信号数据集划分为训练集和测试集。
其中,电力系统的谐波信号是因基波偏差、间谐波和噪声等复杂多干扰因素影响而产生畸变的谐波失真信号,即在电力系统中因大量非线性设备接入后产生失真现象的电流信号。譬如,响应输出包括幅值时,则采集的数据至少包括电力系统的谐波信号以及谐波幅值,以此为基础构建样本。其中谐波幅值是可以采用现有技术手段来获取的,在满足精度要求的基础上,本发明并不约束其技术手段。若包括其他响应,则同样采用现有技术手段来获取满足精度要求的谐波参数。
步骤2:利用所述原始谐波信号数据集训练出谐波参数检测模型。具体包括以下步骤:
步骤2-1:初始化设置。
本实施例中,宽度学习系统的参数:特征映射层节点组数a=8,每组节点数b=8,增强层节点数m=8,输出层节点数N=8,特征映射层—增强层的权值Wf和偏置βf以及增强层—输出层的权值We和偏置βe随机初始化,输出层的权值则在后续步骤中采用伪逆算法初始化。RAdam优化器的参数包括初始步长α0=6×10-7、近似简单移动平均线(SMA)阈值γ=4、一阶矩估计指数衰减率μ1=0.9,二阶矩估计指数衰减率μ2=0.995;其他可行的实施例中,μ1和μ2的范围均为[0,1),初始步长α0可设为其他定值。后续每一步的步长αt也可以设置为定值或者设置为随迭代次数减小。应当注意的是,在不同的谐波参数的估计应用中,以上参数的设定可根据具体的任务需求/精度进行调整。
步骤2-2:将训练集中的谐波信号输入所述宽度学习系统得到映射特征与增强特征之间的联合特征。其中,利用训练集以及初始化设置执行训练过程,具体如下:
S201:将训练集中的原始谐波信号X映射为a组映射特征,第i组映射特征表示为Fi,其计算公式如下:
Figure BDA0003719224430000071
其中,X表示由训练集中的原始谐波信号构成的L个独立样本,每个独立样本为D维,该D维对应一个独立样本的采样点数,i=1,2,...,a,Wfi
Figure BDA0003719224430000072
为第i组特征映射层—增强层的权值和偏置项,ζi为第i个特征映射函数,令Wf=[Wf1,...,Wfa]和
Figure BDA0003719224430000081
本实施例中独立样本数目为800个(由400个整周期信号分解而来),若采用
Figure BDA0003719224430000082
周期,每个独立样本包含62个采样点,若采用
Figure BDA0003719224430000083
周期,每个独立样本包含31个采样点。
本实施例中,优选在输入宽度学习系统前采用公式:
Figure BDA0003719224430000084
对谐波信号进行了归一化处理,式中,x为归一化后的采样点幅值数据,x0为原始谐波信号,xmax和xmin分别为原始谐波信号的最大值和最小值。本实施例中优选特征映射函数ζi采用线性函数。其他可行的实施例中,根据任务需求和具体的预测精度确定特征映射函数的类型,只是基于目前的验证情况来看,选择线性函数的效果最佳。应当注意,本实例中训练时每次输入的信号为单个独立样本。
S202:将所有映射特征串联组成映射特征矩阵F=[F1,F2,...,Fa];
S203:将所述映射特征矩阵F映射为m个增强特征,第j个增强特征表示为Hj,其计算公式如下:
Figure BDA0003719224430000085
其中,i=1,2,...,m,Wej
Figure BDA0003719224430000086
为第j个增强层—输出层的权值和偏置项,ξj为第j个增强函数,记We=[We1,...,Web]以及
Figure BDA0003719224430000087
本实施例中增强函数采用双曲正切函数,其表达式和导数表达式如下:
Figure BDA0003719224430000088
Figure BDA0003719224430000089
其中,v表示输入到增强函数的采样点。应当理解,其他可行的实施例中,根据具体的需求以及精度可以选择其他类型的激活函数。
S204:所有增强特征串联组成增强特征矩阵H=[H1,H2,...,Hm];
S205:映射特征矩阵F和增强特征矩阵H的联合特征构成宽度学习系统的输入特征A=[F|H]。
步骤2-3:判断是否需要初始化所述宽度学习系统中输出层的权值。其中,若是第一次进入谐波参数估计模型,则需要初始化,初始化后进入下一次的迭代训练;否则,不需要初始化输出层权值,而是采用RAdam优化器更新宽度学习系统参数,执行步骤S2-4。
本实施例中,如果需要初始化输出层的权值,期望的网络输出Y=AW,进而输出层矩阵W=A+Y,其中,A+表示对联合特征A进行伪逆运算,期望的网络输出选择是训练集中谐波信号对应的幅值。
步骤2-4:采用RAdam优化器更新所述宽度学习系统的参数,其中,更新的参数包括:特征映射层-增强层、增强层-输出层的权值和偏置以及输出层的权值。具体是基于待更新参数的梯度信息进行优化,梯度信息为基于目标函数求待优化/更新参数的梯度,因此,RAdam的本质还是基于梯度下降思路。
其中,基于RAdam优化器的参数学习更新策略为:
Figure BDA0003719224430000091
其中,t为时间步数,即参数迭代更新的次数;θ为需要学习的参数,具体为Wf,Weef,WF,WH;gt-1为更新前参数,gt为更新后参数,αt为步长,kt为方差修正项,
Figure BDA0003719224430000092
为偏差修正移动平均线,lt为自适应学习率。
本实施例中,第t步的近似SMA值γt和其最大值γmax的计算公式如下:
γmax=2/(1-μ2)-1
Figure BDA0003719224430000093
当γt>γ时,kt
Figure BDA0003719224430000094
和lt的计算方式如下:
Figure BDA0003719224430000095
否则,取kt=lt=1;
其中,mt为有偏一阶矩估计,nt为有偏二阶矩估计,其计算公式如下:
mt=μ1mt-1+(1-μ1)gt
Figure BDA0003719224430000096
其中,gt为待更新参数的梯度信息。其满足:
Figure BDA0003719224430000097
E为目标函数,其与谐波参数的模型输出值以及真实值有关。应当理解,目标函数不同,待更新参数的梯度的计算方式也不同。本实施例中选择真实值和响应值之间的总均方误差E为目标函数,其他可行的实施例中,可以选择其他类型的目标函数。总均方误差E表示为:
Figure BDA0003719224430000098
其中,n=1,...,N,N表示输出节点个数,
Figure BDA0003719224430000099
是谐波参数的期望值,Yn为系统响应值。本实施例中,系统响应值为模型输出,具体如下:
取X的第l个独立样本,记
Figure BDA0003719224430000101
第l个独立样本所对应的参数记为:F=(fi)1×(a×b),H=(hi)1×m,
Figure BDA0003719224430000102
输出层权值W分解为两部分
Figure BDA0003719224430000103
其中
Figure BDA0003719224430000104
连接F到输出层,
Figure BDA0003719224430000105
连接H到输出层,则宽度学习系统的第n个响应输出(譬如,本实施例中包括为幅值的响应输出n)可表示为:
Figure BDA0003719224430000106
其中,
Figure BDA0003719224430000107
Figure BDA0003719224430000108
为第l个独立样本中的第d个采样点。
按照上述目标函数,各类待更新的谐波参数检测模型的模型参数的梯度如下所示:
Figure BDA0003719224430000109
其中,fi,hi,
Figure BDA00037192244300001010
均表示第l个独立样本所对应的映射特征矩阵F、增强特征矩阵H、特征映射层—增强层的权值Wf和偏置βf、增强层—输出层的权值We和偏置βe、输出层权值W中WF、WH中的元素。h'j是hj的导数。
步骤S2-5:判断是否达到迭代终止条件,若达到,训练结束,得到所述谐波参数检测模型的最优模型参数;否则,返回步骤S2-2进入下一次的迭代。
其中,训练迭代终止条件是达到最高训练次数或者最低误差阈值,本实例中最高训练次数是100轮,最低误差阈值为1×10-7
步骤S2-6:输入测试集数据,利用训练好的宽度学习系统估计谐波参数并计算估计误差。
综上,本发明的训练和测试样本是由输入数据和输出数据组成的,输入数据就是采集的谐波信号数据,输出数据是期望的谐波信号参数。因此训练过程的简述为:以训练样本中每个独立样本(谐波信号数据)对应映射特征和增强特征的联合特征为输入,得到输出结果(BLS预测的谐波信号参数)与训练样本的输出数据相对应并衡量误差,从而以此为标准进行训练,确定特征映射层参数、增强层参数和输出层的参数。
且本实施例中优选训练阶段是将独立样本一个一个或一组一组输入宽度学习系统中进行训练,即下一次迭代使用下一个或下一组独立样本,若达到了最后一个样本,则返回从第一个样本重新开始,直至达到迭代终止条件。
步骤S3:获取待检测的电力系统中的谐波信号并输入至所述谐波参数检测模型得到谐波信号的谐波参数。
应当理解,本实施例的谐波参数检测方法采用特征映射层和增强层的联合特征组成网络输入,以及基于宽度学习系统(BLS)学习参数的梯度信息,建立基于梯度下降策略的参数学习机制,提高了宽度学习系统的网络性能;网络参数初始化时,输出层的权值初始化采用伪逆算法,减少因参数随机生成而训练过程耗时,同时也可以提高模型的估计精度。
实施例2:
基于上述谐波参数检测方法,本发明为了验证该方法的性能,搭建了如图3所示的谐波参数检测方法的模拟系统,是为了验证本发明提出的谐波参数检测方法性能而搭建的,其包括:谐波信号采集子系统和谐波信号处理子系统。该谐波信号采集子系统获取谐波信号,包括:输入交流电源、多脉冲整流电源(非线性电力电子设备)、电流探头、可编程大功率直流电子负载、示波器、数字万用表、功率分析仪以及信号采集单元以及信号存储单元。
其中,输入交流电源与多脉冲整流电源连接,所述输入交流电源向整个信号采集子系统提供可变的输出电压、输出频率和负载电流,输出正弦电压加压于多脉冲整流电源,使得基波电流畸变而产生谐波;
其中,所述电流探头是将待测信号完整可靠地传输到示波器的信号链路;
其中,所述示波器用于可视化电路中输入和输出信号波形;
其中,所述数字万用表用于测量电压、电流和电阻值;
其中,所述功率分析仪用于监测谐波数量、功率因数等参数;
其中,所述可编程大功率直流电子负载用于模拟实际负载且大小可调;
其中,所述信号采集单元用于采集畸变谐波信号,并与所述信号存储单元相连,将采样后的信号数据传输到信号存储单元进行存储,以便谐波信号处理子系统进行后续分析;
其中,所述谐波信号处理子系统设有信号处理分析单元,用于根据本发明技术方案提供的所述谐波参数估计方法估计各次谐波参数。
基于上述系统,本实施例中可变频交流电源输出值设定为电压115V、频率400Hz,可编程大功率直流电子负载的设定电压为270V、电流为5A,电力电子设备为400Hz的18脉冲整流电源,并采用上述系统以49.6kHz的采样率获取了400组整周期的谐波信号,训练集和测试集比例设定为7:3。应当注意的是,在本发明的谐波参数估计中,400个整周期的谐波信号将根据估计需要进一步拆分为
Figure BDA0003719224430000111
周期和
Figure BDA0003719224430000112
周期信号,原始信号数据集的输出值,则采用12个周期的已获取真实谐波信号FFT分析值。
为了验证本发明所述方法的可靠性,本实例中采集了400组整周期谐波信号,其波形如图4所示,上述信号进一步分解为800组半周期信号,本发明提出的基于RAdam优化宽度学习的谐波参数检测方法采用
Figure BDA0003719224430000121
周期和
Figure BDA0003719224430000122
周期两种输入信号,
Figure BDA0003719224430000123
周期的信号选用
Figure BDA0003719224430000124
周期的任意一半数据组成,按照70%的训练集比例,训练数据为560组,测试数据为240组。谐波参数估计模型的真实值采用12个整周期的快速傅里叶变换(FFT)分析值,FFT的分析结果如图5所示,本实例中选择1、3、5、7、9、11、13、15次谐波进行估计,其幅值参数进一步如图6所示。首先,基于RAdam优化的宽度学习谐波估计方法进行参数初始化,具体地,本实例中需要优化的参数包括特征映射层—增强层的权值Wf和偏置βf、增强层—输出层的权值We和偏置βe、输出层的权值WF、WH,以总均方误差作为优化目标函数,图7展示了100次迭代过程中每轮迭代结束谐波参数估计模型响应值和测试集真实值之间的总均方误差的平均值,由图可知,采用伪逆算法初始化输出层权值使得模型初始误差值不高,且随着迭代次数的增加,误差值越来越小,并逐渐趋于稳定,
Figure BDA0003719224430000125
周期和
Figure BDA0003719224430000126
周期的输入数据相比,前者的误差更小且参数寻优效果更好。训练过程结束后,即得到了基于RAdam优化宽度学习的谐波参数检测模型的最优权值和阈值,测试过程中误差以测试集240组输入信号和系统响应值之间的相对误差平均值和最大值衡量,并与经典的FFT、WPT两种方法的估计结果进行比较,其中FFT和WPT采用1个和2个整周期的输入信号,其误差结果比较依次如图8和图9所示,可以看出,FFT和WPT对基波的参数估计仍然是相当有效的,当信号周期数由1变成2时,FFT和WPT两种方法的误差显著减小,但因为基频偏差、间谐波和噪声等多复杂因素的干扰,对于高次谐波的估计,FFT、WPT和本发明提出的方法相比误差较大,当输入信号由
Figure BDA0003719224430000127
周期减少为
Figure BDA0003719224430000128
周期,本发明提出方法的误差虽有细微的增加,但并无大的波动,因此在谐波参数估计应用中,从实时性和精确度两个角度比较,实时性需求占优者可采用
Figure BDA0003719224430000129
周期输入信号,精确度需求占优者可采用
Figure BDA00037192244300001210
周期输入信号。通过本实例中的初步测试,本发明的最初的想法也得到验证,即在大量非线性设备接入的电力系统中信号畸变严重的情况下,且存在基波偏差、间谐波和噪声等多复杂因素的干扰时,采用基于RAdam优化宽度学习的谐波参数检测方法能够实现谐波信号参数的实时、准确估计。
实施例3:
本实施例提供一种基于所示谐波参数检测方法的系统,其包括:相互连接或依次连接的原始谐波信号数据集构建模块、谐波参数检测模型构建模块以及检测模块。
其中,原始谐波信号数据集构建模块,用于获取电力系统中的谐波信号以及所述谐波信号对应的谐波参数,进而构建出原始谐波信号数据集;
谐波参数检测模型构建模块,用于利用所述原始谐波信号数据集训练出谐波参数检测模型,其中,所述谐波检测模型是基于宽度学习系统构建的,并采用RAdam优化模型参数;
检测模块,用于将待检测的电力系统中的谐波信号输入至所述谐波参数检测模型得到谐波信号的谐波参数。
具体各个模块的实现过程请参照上述方法的内容,在此不再赘述。应该理解到,上述功能模块的划分仅仅是一种逻辑功能的划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。同时,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
实施例4:
本实施例提供一种电子终端,其包括一个或多个处理器以及存储了一个或多个计算机程序的存储器,该处理器调用计算机程序以实现:一种基于RAdam优化宽度学习的谐波参数检测方法的步骤。具体执行:
步骤S1:获取电力系统中的谐波信号以及所述谐波信号对应的谐波参数,进而构建出原始谐波信号数据集;
步骤S2:利用所述原始谐波信号数据集训练出谐波参数检测模型,其中,所述谐波检测模型是基于宽度学习系统构建的,并采用RAdam优化模型参数;
步骤S3:获取待检测的电力系统中的谐波信号并输入至所述谐波参数检测模型得到谐波信号的谐波参数。
其中,存储器可能包含高速RAM存储器,也可能还包括非易失性除颤器,例如至少一个磁盘存储器。
如果存储器、处理器独立实现,则存储器、处理器和通信接口可以通过总线相互连接并完成相互间的通信。所述总线可以是工业标准体系结构总线,外部设备互联总线或扩展工业标准体系结构总线等。所述总线可以分为地址总线、数据总线、控制总线等。
可选的,在具体实现上,如果存储器、处理器集成在一块芯片上,则存储器、处理器可以通过内部接口完成相互之间的通信。
各个步骤的具体实现过程请参照前述方法的阐述。
应当理解,在本发明实施例中,所称处理器可以是中央处理单元(CentralProcessing Unit,CPU),该处理器还可以是其他通用处理器、数字信号处理器(DigitalSignal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。存储器可以包括只读存储器和随机存取存储器,并向处理器提供指令和数据。存储器的一部分还可以包括非易失性随机存取存储器。例如,存储器还可以存储设备类型的信息。
实施例5:
本实施例提供一种可读存储介质,其存储了计算机程序,存储了计算机程序,所述计算机程序被处理器调用以实现:一种基于RAdam优化宽度学习的谐波参数检测方法的步骤。具体执行:
步骤S1:获取电力系统中的谐波信号以及所述谐波信号对应的谐波参数,进而构建出原始谐波信号数据集;
步骤S2:利用所述原始谐波信号数据集训练出谐波参数检测模型,其中,所述谐波检测模型是基于宽度学习系统构建的,并采用RAdam优化模型参数;
步骤S3:获取待检测的电力系统中的谐波信号并输入至所述谐波参数检测模型得到谐波信号的谐波参数。
各个步骤的具体实现过程请参照前述方法的阐述。
所述可读存储介质为计算机可读存储介质,其可以是前述任一实施例所述的控制器的内部存储单元,例如控制器的硬盘或内存。所述可读存储介质也可以是所述控制器的外部存储设备,例如所述控制器上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述可读存储介质还可以既包括所述控制器的内部存储单元也包括外部存储设备。所述可读存储介质用于存储所述计算机程序以及所述控制器所需的其他程序和数据。所述可读存储介质还可以用于暂时地存储已经输出或者将要输出的数据。
基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分,或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的可读存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,RandomAccess Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
需要强调的是,本发明所述的实例是说明性的,而不是限定性的,因此本发明不限于具体实施方式中所述的实例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,不脱离本发明宗旨和范围的,不论是修改还是替换,同样属于本发明的保护范围。

Claims (10)

1.一种基于RAdam优化宽度学习的谐波参数检测方法,其特征在于:包括以下步骤:
步骤S1:获取电力系统中的谐波信号以及所述谐波信号对应的谐波参数,进而构建出原始谐波信号数据集;
步骤S2:利用所述原始谐波信号数据集训练出谐波参数检测模型,其中,所述谐波检测模型是基于宽度学习系统构建的,并采用RAdam优化模型参数;
其中,所述谐波参数检测模型的模型输入为谐波信号,模型输出为谐波信号对应的谐波参数;
步骤S3:获取待检测的电力系统中的谐波信号并输入至所述谐波参数检测模型得到谐波信号的谐波参数。
2.根据权利要求1所述的谐波参数检测方法,其特征在于:步骤S2中所述谐波参数检测模型的训练过程如下:
步骤S2-1:初始化设置;
设置所述谐波参数检测模型中宽度学习系统内特征映射层-增强层、增强层-输出层的权值和偏置;
步骤S2-2:将所述原始谐波信号数据集中的谐波信号输入所述宽度学习系统得到映射特征与增强特征之间的联合特征;
步骤S2-3:判断是否需要初始化所述宽度学习系统中输出层的权值;
其中,若是第一次进入训练所述谐波参数检测模型,则利用W=A+Y得到输出层的权值,进入下一次迭代训练,返回步骤S2-2;否则,执行步骤S2-4;W表示输出层的权值,Y表示期望的网络输出,并选用所述原始谐波信号数据集中谐波信号的谐波参数表示;A+表示对联合特征A进行伪逆运算;
步骤S2-4:采用RAdam优化器更新所述宽度学习系统的参数,其中,更新的参数包括:特征映射层-增强层、增强层-输出层的权值和偏置以及输出层的权值;
步骤S2-5:判断是否达到迭代终止条件,若达到,训练结束,得到所述谐波参数检测模型的最优模型参数;否则,返回步骤S2-2进入下一次的迭代。
3.根据权利要求2所述的谐波参数检测方法,其特征在于:所述RAdam优化器的更新策略为:
Figure FDA0003719224420000011
其中,t为时间步数,即参数迭代更新的次数;θtt-1表示更新后、更新前的谐波参数检测模型的模型参数,αt为步长,kt为方差修正项,
Figure FDA0003719224420000012
为偏差修正移动平均线,lt为自适应学习率。
4.根据权利要求3所述的谐波参数检测方法,其特征在于:所述方差修正项,偏差修正移动平均线以及自适应学习率的计算如下:
γmax=2/(1-μ2)-1
Figure FDA0003719224420000021
当γt>γ时,kt
Figure FDA0003719224420000022
和lt满足:
Figure FDA0003719224420000023
否则,取kt=lt=1以及
Figure FDA0003719224420000027
其中,γmax、γt、γ分别表示第t步的近似SMA的最大值、近似SMA值、近似SMA的阈值;μ1为一阶矩估计指数衰减率,μ2为二阶矩估计指数衰减率,
Figure FDA0003719224420000024
为μ2的t次方,mt为有偏一阶矩估计,nt为有偏二阶矩估计,其计算公式如下:
mt=μ1mt-1+(1-μ1)gt
Figure FDA0003719224420000025
其中,gt为待更新参数θ的梯度,其满足:
Figure FDA0003719224420000026
E为目标函数,其与谐波参数的模型输出值以及真实值有关。
5.根据权利要求2所述的谐波参数检测方法,其特征在于:所述宽度学习系统内的增强函数为双曲正切函数,特征映射函数采用线性函数。
6.根据权利要求1所述的谐波参数检测方法,其特征在于:所述原始谐波信号数据集中谐波信号对应的独立样本是将整周期的谐波信号分解为1/2周期或1/4周期得到的。
7.一种基于权利要求1-6任一项所述谐波参数检测方法的模拟系统,其特征在于:包括相互连接的谐波信号采集子系统以及谐波信号处理子系统;
其中,所述谐波信号采集子系统用于采集电力系统的谐波信号;所述谐波信号处理子系统用于根据步骤S2方式进行处理谐波信号进而得到谐波参数检测模型,其中,将待检测的电力系统中的谐波信号输入至所述谐波参数检测模型得到谐波信号的谐波参数;
所述谐波信号采集子系统至少包括:输入交流电源、非线性电力电子设备、可编程大功率直流电子负载、以及信号采集单元以及信号存储单元;
其中,输入交流电源与非线性电力电子设备连接,所述输入交流电源向所述谐波信号采集子系统提供可变的输出电压、输出频率和负载电流,并作用于非线性电力电子设备进而产生谐波;
所述可编程大功率直流电子负载用于模拟实际负载且大小可调;
所述信号采集单元用于采集畸变的谐波信号,并与所述信号存储单元相连,将采样后的信号数据传输到信号存储单元进行存储,以便谐波信号处理子系统进行后续分析。
8.一种基于权利要求1-6任一项所述谐波参数检测方法的系统,其特征在于:包括:
原始谐波信号数据集构建模块,用于获取电力系统中的谐波信号以及所述谐波信号对应的谐波参数,进而构建出原始谐波信号数据集;
谐波参数检测模型构建模块,用于利用所述原始谐波信号数据集训练出谐波参数检测模型,其中,所述谐波检测模型是基于宽度学习系统构建的,并采用RAdam优化模型参数;
其中,所述谐波参数检测模型的模型输入为谐波信号,模型输出为谐波信号对应的谐波参数;
检测模块,用于将待检测的电力系统中的谐波信号输入至所述谐波参数检测模型得到谐波信号的谐波参数。
9.一种电子终端,其特征在于:包括:
一个或多个处理器;
存储了一个或多个计算机程序的存储器;
所述处理器调用所述计算机程序以实现:
权利要求1-6任一项所述谐波检测方法的步骤。
10.一种可读存储介质,其特征在于:存储了计算机程序,所述计算机程序被处理器调用以实现:
权利要求1-6任一项所述谐波检测方法的步骤。
CN202210778084.6A 2022-06-28 2022-06-28 基于RAdam优化宽度学习的谐波参数检测方法和系统 Pending CN115062668A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210778084.6A CN115062668A (zh) 2022-06-28 2022-06-28 基于RAdam优化宽度学习的谐波参数检测方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210778084.6A CN115062668A (zh) 2022-06-28 2022-06-28 基于RAdam优化宽度学习的谐波参数检测方法和系统

Publications (1)

Publication Number Publication Date
CN115062668A true CN115062668A (zh) 2022-09-16

Family

ID=83204523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210778084.6A Pending CN115062668A (zh) 2022-06-28 2022-06-28 基于RAdam优化宽度学习的谐波参数检测方法和系统

Country Status (1)

Country Link
CN (1) CN115062668A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116596061A (zh) * 2023-05-16 2023-08-15 北京航空航天大学 一种基于迁移学习的功率放大器谐波预测方法
CN116595368A (zh) * 2023-05-16 2023-08-15 北京航空航天大学 一种基于非线性建模的功率放大器谐波预测方法
CN117096956A (zh) * 2023-10-20 2023-11-21 江苏力普电子科技有限公司 一种高压变频器的谐波控制方法及系统

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116596061A (zh) * 2023-05-16 2023-08-15 北京航空航天大学 一种基于迁移学习的功率放大器谐波预测方法
CN116595368A (zh) * 2023-05-16 2023-08-15 北京航空航天大学 一种基于非线性建模的功率放大器谐波预测方法
CN116596061B (zh) * 2023-05-16 2023-12-12 北京航空航天大学 一种基于迁移学习的功率放大器谐波预测方法
CN116595368B (zh) * 2023-05-16 2024-01-26 北京航空航天大学 一种基于非线性建模的功率放大器谐波预测方法
CN117096956A (zh) * 2023-10-20 2023-11-21 江苏力普电子科技有限公司 一种高压变频器的谐波控制方法及系统
CN117096956B (zh) * 2023-10-20 2023-12-26 江苏力普电子科技有限公司 一种高压变频器的谐波控制方法及系统

Similar Documents

Publication Publication Date Title
US11409347B2 (en) Method, system and storage medium for predicting power load probability density based on deep learning
CN115062668A (zh) 基于RAdam优化宽度学习的谐波参数检测方法和系统
CN111860982A (zh) 一种基于vmd-fcm-gru的风电场短期风电功率预测方法
CN110443724B (zh) 一种基于深度学习的电力系统快速状态估计方法
CN110413657B (zh) 面向季节型非平稳并发量的平均响应时间评估方法
CN116316591A (zh) 基于混合双向门控循环的短期光伏功率预测方法及系统
CN112884213B (zh) 小波分解与动态混合深度学习的燃煤锅炉NOx预测方法
CN116707331B (zh) 基于模型预测的逆变器输出电压高精度调节方法及系统
CN113095598A (zh) 一种多能负荷预测方法、系统、设备和介质
CN110879927A (zh) 一种用于海目标检测的海杂波幅度统计分布现场建模方法
CN114490065A (zh) 一种负载预测方法、装置及设备
CN115201563A (zh) 一种基于联合熵的多谐波源定位方法及系统
US20220243347A1 (en) Determination method and determination apparatus for conversion efficiency of hydrogen production by wind-solar hybrid electrolysis of water
CN114611676A (zh) 基于神经网络的新能源发电系统阻抗模型辨识方法及系统
CN107563511B (zh) 一种实时系统可用时间快速估算与优化方法
CN112287605B (zh) 一种基于图卷积网络加速的潮流校核方法
CN113553771B (zh) 基于rnn网络的动态x参数核计算方法
CN113642784B (zh) 一种计及风机状态的风电功率超短期预测的方法
CN112685900B (zh) 一种表征冲击负荷功率特性的电力负荷模拟方法
CN115169740A (zh) 基于压缩感知的池化回声状态网络的序列预测方法及系统
CN115632970A (zh) 非高斯噪声下通信干扰信号带宽估计方法、设备及存储介质
CN115189990A (zh) 一种基于循环储蓄学习网络的ofdm时域信道预测方法及系统
CN114298300A (zh) 不确定度预测方法、装置、电子设备及存储介质
Wang et al. Power Quality Disturbance Recognition Method in Park Distribution Network Based on One-Dimensional VGGNet and Multi-label Classification
CN111310907A (zh) 一种微波组件故障诊断方法、装置及设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination