一种基于完备总体经验模态分解和极限学习机的风速预测方法技术

本发明专利技术公开一种基于完备总体经验模态分解和极限学习机的风速预测方法，该方法首先利用完备总体经验模态分解对不稳定、随机的原始风速序列进行分解得到一系列的稳定的固有模态分量和一个残差序列，然后利用极限学习机对每个固有分量和残差序列进行训练预测得到各自的子预测结果，最后把所有的子预测结果重构得到最后的风速预测结果。通过与其它3个风速预测模型相比较可知，本发明专利技术提出的方法不仅提高了风速预测的精度，还增强了模型的鲁棒性和训练预测速率。

【技术实现步骤摘要】
一种基于完备总体经验模态分解和极限学习机的风速预测方法
本专利技术涉及机器学习和风力发电的领域，特别是涉及一种基于完备总体经验模态分解和极限学习机的风速预测方法。
技术介绍
随着可再生能源的发展和化石能源的不断减少，风电作为一种清洁的可再生能源越来越受到人们的重视。因此，充分地利用风能来发电成为了一项重要的任务。然而，由于风速具有不断变化的行为，所以风速序列具有随机性、非稳定性和非线性的特点。由于风速的不稳定性，将会导致风电电能不稳定，进而导致整个电力系统不稳定。所以，对风速进行预测并提前采取相应措施来保证电力系统的稳定性是很有必要的。近年来，为了提高风速预测的精度，各国专家学者提出了大量的风速预测模型。目前，有物理和统计两种主流的方法；另外，一些混合预测方法也越来越受重视。物理模型通过物理因素去估计未来的风速，数值天气预报是最常见的物理模型；统计方法主要通过使用大量历史数据去构建模型，常见的统计模型有自回归模型、移动平均模型、自回归移动平均模型和自回归整合移动平均模型。由于物理和统计方法的局限性，许多混合方法被提出。一般来说，混合预测模型比单个预测模型有更好的预测性能。最常见的混合策略有两种，第一种混合策略是先通过信号分解算法对原始风速序列进行分解，把不稳定的原始序列分解成一些稳定的子序列，然后再把每个子序列应用到预测方法上；第二种混合策略是通过智能优化算法优化风速预测模型的参数。然而，以上两种混合策略仍存在一些问题，信号分解常常存在模态混叠问题、智能优化算法常常陷入局部最优解和预测模型的预测性能这些问题都会对预测精度产生影响。于是如今的混合趋势是减少信号分解算法的模态混叠、提高智能优化算法的鲁棒性和提高预测模型的预测性能。总体完备经验模态分解是经验模态算法的改进算法。经验模态算法适用于分析非线性和非稳定的序列，它能把原始序列分解成一系列不同频率的稳定子序列，这些子序列称作固有模态分量。固有模态分量在Hilbert转化后能得到瞬时属性的物理意义。然而，经验模态算法也存在模态混叠的缺陷，模态混叠严重影响了分解的精度和破坏了包含在固有模态分量中的物理意义。为了解决这一问题，Torres等人在2011年提出了总体完备经验模态分解算法，该算法在分解的每一步都会给信号添加一个特殊的噪声去平均模态混叠，该改进算法大大提高了信号分解的精度。由于一般的前馈神经网络的学习速率比要求的要慢。为了提高神经网络的学习速率，Huang等人在2006提出了一种新的单隐藏层前馈神经网络叫极限学习机。极限学习机基于经验风险最小化理论并且它的学习过程只需要单迭代，相比于其它神经网络，它具有更好的泛化能力，鲁棒性和快速学习速率。因此，极限学习机会提高预测的精度。
技术实现思路
本专利技术提出一种基于完备总体经验模态分解和极限学习机的风速预测方法，来解决现有的低精度和低鲁棒性的风速预测方法。为实现上述目的，本专利技术采取以下技术方案：一种基于完备总体经验模态分解和极限学习机的风速预测方法包括如下步骤：S1、采集原始风速数据，构建完备总体经验模态分解和极限学习机的风速预测模型，利用完备总体经验模态分解对原始的风速序列进行分解得到多个稳定的固有模态分量；完备总体经验模态分解对原始风速序列的分解步骤为：S11、给原始信号s[τ]添加高斯白噪声得到信号s[τ]+ε0ωi[τ]，其中ωi表示高斯白噪声，ε表示信噪比；S12、对s[τ]+ε0ωi[τ]重复进行I次经验模态分解可得到第1个固有模态分量为：则第1个残差为：其中，E1(·)是由EMD算法得到的第1个固有模态分量；S13、继续分解ρ1[τ]+ε1E1(ω1[τ])，直到满足(2)中得到的固有模态分量的条件时结束，把全部分解信号的均值作为第2个固有模态分量为：于是对于m＝2,…,M，则第m个残差为：S14、分解满足第1个固有模态分量的条件时结束，把全部分解信号的均值作为第m+1个固有模态分量为：S15、重复步骤(4)，残差不能分解时重复，得到最终的残差为：其中，M表示固有模态分量的数量。S2、针对每个固有模态分量确定各自的训练集和测试集；S3、把每个固有模态分量确定的训练集分别作为极限学习机的输入进行训练，得到它们各自的训练特征；对于N个训练样本{(ξi,τi)|ξi∈Rn,τi∈R,i＝1,…,N}，包含m个节点并且激励函数为f(ξ)的极限学习机模型为：其中，ωi表示输入层中的n个节点分别到隐藏层的m个节点的权重向量，βj表示隐藏层的m个节点到输出层的m个节点的权重向量，ξj表示隐藏层的偏置向量，则极限学习机的训练预测的步骤为：S31、随机地设定权值向量ωi,i＝1,2…,m和偏置向量ξj,j＝1,2,…,m；S32、利用激励函数f(x)计算n×m的转移矩阵T；S33、计算输出节点的权值向量为：β＝T+Y其中，T+表示矩阵T的Moore-Pennrose转置。S4、对每个固有模态分量对应的测试集进行预测，得到每个固有模态分量对应的子预测结果；S5、将S4中得到的各个子预测结果进行重构，得到最终的预测值结果输出；本专利技术的技术效果和优点是：本专利技术的基于完备总体经验模态分解和极限学习机的风速预测方法，首先把非平稳的原始风速序列通过完备总体经验模态分解算法分解成一系列平稳的固有模态分量，然后利用极限学习机对每个固有模态函数进行训练预测，最后把所有预测的分量重构得到最终的预测结果。极限学习机提高了学习能力和学习速率，完备总体经验模态分解则是提高了预测模型的预测精度。附图说明图1是本专利技术模型总流程图；图2是极限学习机的结构图；图3(a)是完备总体经验模态算法分解原始序列得到的第一个固有模态分量图；图3(b)是完备总体经验模态算法分解原始序列得到的第二个固有模态分量图；图3(c)是完备总体经验模态算法分解原始序列得到的第三个固有模态分量图；图3(d)是完备总体经验模态算法分解原始序列得到的第四个固有模态分量图；图3(e)是完备总体经验模态算法分解原始序列得到的第五个固有模态分量图；图3(f)是完备总体经验模态算法分解原始序列得到的第六个固有模态分量图；图3(g)是完备总体经验模态算法分解原始序列得到的第七个固有模态分量图；图3(h)是完备总体经验模态算法分解原始序列得到的第八个固有模态分量图；图3(i)是完备总体经验模态算法分解原始序列得到的第九个固有模态分量图；图3(j)是完备总体经验模态算法分解原始序列得到的第十个固有模态分量图；图3(k)是完备总体经验模态算法分解原始序列得到的残差图；图4是完备总体经验模态和极限学习机模型预测的风速与原始风速比较图；图5是单一极限学习机模型预测的风速与原始风速比较图；图6是经验模态算法和极限学习机模型预测的风速与原始风速比较图；图7是集合经验模态和极限学习机模型预测的风速与原始风速比较图。具体实施方式为了使本专利技术的优点更加清除，以下将结合附图以及具体的实施方式，对上述技术方案进行详细的说明。如图1所示，实例一种基于完备总体经验模态分解和极限学习机的风速预测方法包括以下步骤：S1、采集原始风速数据，利用完备总体经验模态分解对原始的风速序列进行分解得到多个稳定的固有模态分量和残差序列；完备总体经验模态分解对原始风速序列的分解步骤为：S11、给原始信号s[τ]添加高斯白噪本文档来自技高网...

【技术保护点】
一种基于完备总体经验模态分解和极限学习机的风速预测方法，其特征在于，包括以下步骤：S1、采集原始风速数据，构建完备总体经验模态分解和极限学习机的风速预测模型，利用完备总体经验模态分解对原始的风速序列进行分解得到多个稳定的固有模态分量；S2、针对每个固有模态分量确定各自的训练集和测试集；S3、把每个固有模态分量确定的训练集分别作为极限学习机的输入进行训练，得到它们各自的训练特征；S4、对每个固有模态分量对应的测试集进行预测，得到每个固有模态分量对应的子预测结果；S5、将S4中得到的各个子预测结果进行重构，得到最终的预测值结果输出；S6、上述模型对比单个极限学习机和经验模态分解与极限学习机这两个模型的性能指标，得到性能最优的风速预测模型。

【技术特征摘要】
1.一种基于完备总体经验模态分解和极限学习机的风速预测方法，其特征在于，包括以下步骤：S1、采集原始风速数据，构建完备总体经验模态分解和极限学习机的风速预测模型，利用完备总体经验模态分解对原始的风速序列进行分解得到多个稳定的固有模态分量；S2、针对每个固有模态分量确定各自的训练集和测试集；S3、把每个固有模态分量确定的训练集分别作为极限学习机的输入进行训练，得到它们各自的训练特征；S4、对每个固有模态分量对应的测试集进行预测，得到每个固有模态分量对应的子预测结果；S5、将S4中得到的各个子预测结果进行重构，得到最终的预测值结果输出；S6、上述模型对比单个极限学习机和经验模态分解与极限学习机这两个模型的性能指标，得到性能最优的风速预测模型。2.根据权利要求1所述的基于完备总体经验模态分解和极限学习机的风速预测方法，其特征在于，步骤S1中的完备总体经验模态分解对原始风速序列的分解步骤为：S1、给原始信号s[τ]添加高斯白噪声得到信号s[τ]+ε0ωi[τ]，其中ωi表示高斯白噪声，ε表示信噪比；S2、对s[τ]+ε0ωi[τ]重复进行I次经验模态分解可得到第1个固有模态分量为：则第1个残差为：其中，E1(·)是由EMD算法得到的第1个固有模态分量；S3、继续分解ρ1[τ]+ε1E1(ω1[τ])，直到满足(2)中得到的固有模态分量的条件时结束，把全部分解信号的均值作为第2个固有模态分量为：

【专利技术属性】
技术研发人员：颜宏文，卢格宇，
申请(专利权)人：长沙理工大学，
类型：发明
国别省市：湖南,43