CN113011669B - 一种生猪月度存栏量预测方法 - Google Patents
一种生猪月度存栏量预测方法 Download PDFInfo
- Publication number
- CN113011669B CN113011669B CN202110339114.9A CN202110339114A CN113011669B CN 113011669 B CN113011669 B CN 113011669B CN 202110339114 A CN202110339114 A CN 202110339114A CN 113011669 B CN113011669 B CN 113011669B
- Authority
- CN
- China
- Prior art keywords
- mode
- algorithm
- decomposition
- predicted value
- modes
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 241000282887 Suidae Species 0.000 title claims abstract description 26
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 50
- 238000012549 training Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000010606 normalization Methods 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 abstract 1
- 239000010410 layer Substances 0.000 description 8
- 238000013528 artificial neural network Methods 0.000 description 7
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 1
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Marine Sciences & Fisheries (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Animal Husbandry (AREA)
- Educational Administration (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Primary Health Care (AREA)
- Artificial Intelligence (AREA)
- Agronomy & Crop Science (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种生猪月度存栏量预测方法,包括:通过变分模态分解算法,将生猪存栏量的原始时间序列数据x(t)分解为K个模态;对于每一个由变分模态分解得到的模态,先进行归一化处理,再按照预设比例划分为训练集和测试集;采用训练集数据对极限学习机算法进行训练,确定算法的最优参数;以步长为V的滑动窗口选取极限学习机算法的输入集数据;对于每一个模态,分别将测试集的输入集数据输入到训练好的极限学习机算法中,输出对下一时点的预测值,将其反归一化处理后得到预测值序列uk(t);对所有模态的预测值进行相加重构,得到最终的预测值结果。本发明方法能够提高预测效率和预测精度,并在实验中验证了预测的有效性。
Description
技术领域
本发明涉及生猪月度存栏量预测技术领域,特别涉及一种基于自适应变分模态分解和极限学习机的生猪月度存栏量预测方法。
背景技术
通过挖掘时间序列数据的内在统计特征对未来时点的数据进行预测,目前已存在于诸多科学领域中。时间序列数据由一系列取自不同时点的观测值构成。挖掘时间序列数据内在统计特征的目的在于找出历史样本点数据之间的变化规律,并构建时间序列模型对未来时点的数据进行样本外预测。一般地,时间序列中不同时点的数据具有口径一致的特点。根据数据的截取方式不同,时序数据既可以是时点数,也可以是时期数。
对我国生猪月度存栏量进行精准预测,将为制定阶段性畜牧业政策提供重要参考。但受经济、环境等多方面因素共同作用的影响,生猪月度存栏量序列一般都是非平稳、非线性的,这两个特征为其精准预测提出了更高的要求。目前被广泛应用到这类时间序列数据预测的方法主要可以分为三种:计量统计模型、基于人工智能方法的预测模型和混合预测模型。计量统计模型主要包括自差分整合移动平均自回归模型(Autoregressiveintegrated moving average,ARIMA)、广义自回归条件异方差模型(Generalizedautoregressive conditional heteroskedasticity model,GARCH)等模型,这类模型在处理时间序列数据的线性特征时取得了很好的效果,但在处理非平稳、非线性的复杂序列数据时的表现往往不佳。近年来,基于人工智能的预测模型逐渐应用到复杂时间序列数据的处理中,具有代表性的有人工神经网络(Artificial neural network,ANN)、支持向量机(Support vector regression,SVM)和极限学习机(Extreme learning machine,ELM)等,其原理在于将输入集数据与输出集数据之间的内在联系映射到高维度进行训练,并通过预设的有限次迭代达到最优参数。
然而,传统ANN采用梯度下降的迭代算法调整权重参数,通常会伴随着训练速度慢,易陷入局部最优和出现过拟合等现象,导致预测结果鲁棒性较差。相较于ANN,SVR在训练速度方面和预测结果鲁棒性得到了有效提升,但其泛化能力对核函数及其参数的选择具有较高的敏感性,需要不断调整参数才能实现较高精度的预测结果。与ANN和SVR相比,ELM是单层前向反馈神经网络学习算法,兼具运算速度快和泛化能力强的优点,大幅提高了运算效率和收敛速度,使预测结果更具鲁棒性。目前,ELM在复杂时间序列数据预测中表现出了良好的性能。
混合预测模型制作有多种实现形式,该类模型同时将至少两个算法结合在一起,以实现更优的预测效果。一方面,混合预测模型可将多个预测模型结合在一起,从不同的角度对复杂时间序列数据进行特征提取并分别作出预测,结合各个地方预测结果得到最终预测结果;另一方面,混合预测模型可将数据预处理算法与预测模型相结合,该类模型基于“分而治之”的思想,先对复杂的原始数据进行降维、分解等预处理,对处理后的数据进行建模预测,再对预测结果进行重构得到最终的预测结果。目前,常见的分解算法包括小波变换、经验模态分解(Empirical mode decomposition,EMD)、集成经验模态分解(Ensembleempirical mode decomposition,EEMD)、自适应噪声完备集成经验模态分解(CEEMDAN),改进的CEEMDAN(Improved complete ensemble empirical model decomposition withadaptive noise,ICEEMDAN)和变分模态分解(Variational mode decomposition,VMD)等。
EMD是一种信号时频分析方法,能实现对非平稳非线性数据的平稳化,对金融时间序列数据的分析起到了很大作用。然而,该方法也存在一些缺陷,主要是模边界效应、模式重叠、对噪声敏感等,这可能会对分解精度产生负面影响,从而导致最终预测结果失真。为了降低以上问题对分解结果的不良影响,EEMD将高斯白噪声加入原始信号中再进行分解,取得了不错的效果。但由于人为选择白噪声的不同,分解得到的本征模函数也有差异,这使得EEMD方法不稳定,而且EEMD方法难以完全消除由加入的高斯白噪声引起的重构误差,这也将导致建立精确预测模型难度的增大。在此基础上,通过将具有自适应特征等性质的白噪声加入EMD,学者陆续提出了互补集成经验模式分解(CEEMDAN)和改进的CEEMDAN(ICEEMDAN),逐步降低了模态混叠问题对分解结果精度的不利影响。然而EEMD等改良的算法仍然存在无法充分分解频率相近分量等问题,导致其在分解高频数据时效果有限。
VMD是一种起源于信号处理的多分辨率技术。与EMD不同,VMD是一种完全非递归算法,可以将原始信号或序列数据分解为频谱中具有指定带宽的多个成分(子序列)。研究表明,VMD在抑制噪声和提高分解精度等方面远远优于同类模型。但是,VMD分解得到的模态个数需要预先设定,但对这一数值的预设缺乏统一的指导,因此将对最终分解和预测结果造成明显的影响。分解模态数量过多,将导致计算资源的浪费,而分解模态数量过少,将导致模态中存在较多的噪声数据,导致精准预测的难度大,最终导致预测精度低的后果。
发明内容
本发明的目的在于提供一种生猪月度存栏量预测方法,通过引入模态分解标准改良变分模态分解(VMD),并基于“分解-集成”框架将其与极限学习机(ELM)相结合,实现对生猪月度存栏量的精准预测。
为解决上述技术问题,本发明的实施例提供如下方案:
一种生猪月度存栏量预测方法,包括以下步骤:
通过变分模态分解算法,将生猪存栏量的原始时间序列数据x(t)分解为K个模态;
对于每一个由变分模态分解得到的模态,先进行归一化处理,再按照预设比例划分为训练集和测试集;
采用训练集数据对极限学习机算法进行训练,确定算法的最优参数;
以步长为V的滑动窗口选取极限学习机算法的输入集数据;
对于每一个模态,分别将测试集的输入集数据输入到训练好的极限学习机算法中,输出对下一时点的预测值,将其反归一化处理后得到预测值序列uk(t);
对所有模态的预测值进行相加重构,得到最终的预测值结果。
优选地,K根据指标rres自适应确定,公式如下:
其中,x(t)表示原始时间序列,N表示总时点数;当rres开始达到小于0.01时,即确定模态个数K。
优选地,对于每一个由变分模态分解得到的模态,进行归一化处理的公式如下:
优选地,所述预设比例为8:2。
优选地,对所有模态的预测值进行相加重构的公式如下:
优选地,所述变分模态分解算法的过程如下:
初始化模态分量的带宽uk,中心频率ωk以及循环次数n;
当ω大于等于0时,循环更新每个模态分量的带宽uk和中心频率ωk,自适应分解成以中心频率ωk为中心扩散的模态分量:
其中,uk为分解后对模态分量,ωk为模态分量对应的中心频率,λ为拉格朗日乘子,当满足下述条件时,终止该循环:
其中,∈为常数;
在计算每个分量时,将信号分布到变分模型中进行分解,使用希尔伯特变换和高斯平滑,通过寻找约束变分模型的最优解实现复杂信号的分解;求解每一个模态的希尔伯特变换,然后将每一个模态的频谱移到基带上,之后使用解调信号的H高斯平滑,最小化模态带宽的总和:
引入二次惩罚因子α和拉格朗日乘子λ,将约束问题转化为非约束问题,采用交替乘子方向法进行求解:
优选地,所述极限学习机算法的实现过程如下:
给定一个样本(xi,yi),其中xi=[xi1,x2,...,xiN]T,表示N维输入集,yi=[yi1,yi2,...,yiM]T,表示M维输出集,i=1,2,...,N表示样本标签;极限学习机算法的计算过程如下:
其中,L表示隐含层节点数量,βl表示第l个节点的输出权重矩阵,F为执行运算的激活函数,wl和bl分别表示隐含层第l个节点的输入权重向量和偏置向量;因此,上述运算又可以表述为Hβ=Y,其中β=[β1,β2,...,βL]T,H表示隐含层输出权重矩阵,表示如下:
因此,将β表示为:βT=H+T=HT(HHT)-1T,其中H+=HT(HHT)-1表示隐含层输出矩阵H的Moore-Penrose广义逆矩阵。
优选地,所述方法还包括:
通过比较预测值序列和实际值序列,对预测效果进行验证。
本发明实施例提供的技术方案带来的有益效果至少包括:
本发明通过引入模态分解标准改良变分模态分解(VMD)算法,并基于“分解-集成”框架将其与极限学习机(ELM)算法相结合,实现对生猪月度存栏量的精准预测。
当前在变分模态分解中,模态个数K需要预先设定,但对这一数值的预设缺乏统一的指导,因此将对最终分解和预测结果造成明显的影响。分解模态数量过多,将导致计算资源的浪费,而分解模态数量过少,将导致模态中存在较多的噪声数据,导致精准预测的难度大,最终导致预测精度低的后果。本发明提出了一个适用于复杂时间序列数据的分解标准,使变分模态分解能够在该标准下自适应地确定最佳分解模态个数K,在节省计算资源的同时实现对原始序列的显著降噪,为后一阶段的预测工作奠定基础。
考虑到传统ANN在采用梯度下降迭代算法调整权重参数时存在的训练速度慢,易陷入局部最优和出现过拟合等问题,以及SVR在泛化能力方面对核函数及其参数的选择具有较高的敏感性,本发明采用了兼具较快运算速度和较强泛化能力的ELM作为主要预测模型,大幅提高了运算效率和收敛速度,使预测结果更具鲁棒性。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种生猪月度存栏量预测方法的流程图;
图2是对未来1阶的数据进行单步预测的示意图;
图3是ELM、EMD-ELM、ICEEMDAN-ELM和VMD-ELM之间的预测结果对比示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明的实施例提供了一种生猪月度存栏量预测方法,如图1所示,所述方法包括以下步骤:
通过变分模态分解(VMD)算法,将生猪存栏量的原始时间序列数据x(t)分解为K个模态,也称本征模函数(IMF)。
其中,K根据指标rres自适应确定,公式如下:
式中,x(t)表示原始时间序列,N表示总时点数;当rres开始达到小于0.01的水平且不存在明显的下降趋势时,即可确定模态个数K。
对于每一个由变分模态分解(VMD)得到的模态,先进行归一化处理,再按照预设比例划分为训练集和测试集。
其中,进行归一化处理的公式如下:
作为本发明的一种具体实施方式,归一化处理后按照8:2的比例将其划分为训练集和测试集。
采用训练集数据对极限学习机算法进行训练,确定算法的最优参数,训练过程中,网络中的参数不断迭代更新。
采用步长为V的滑动窗口,即以V阶的历史数据作为极限学习机(ELM)算法的输入集数据,对未来一阶的数据进行单步预测,如图2所示。
对于每一个模态,分别将测试集的输入集数据输入到训练好的极限学习机(ELM)算法中,输出对下一时点的预测值,将其反归一化处理后得到预测值序列uk(t)。
对所有模态的预测值进行相加重构,得到最终的预测值结果,计算公式如下:
进一步地,所述方法还包括:
通过比较预测值序列和实际值序列,对该模型的预测效果进行验证。
进一步地,在本发明的实施例中,所述变分模态分解(VMD)算法的过程如下:
初始化模态分量的带宽uk,中心频率ωk以及循环次数n;
当ω大于等于0时,循环更新每个模态分量的带宽uk和中心频率ωk,自适应分解成以中心频率ωk为中心扩散的模态分量:
其中,uk为分解后对模态分量,ωk为模态分量对应的中心频率,λ为拉格朗日乘子,当满足下述条件时,终止该循环:
其中,∈为常数;
VMD是一种新型多分量信号分解算法,在计算每个分量时,将信号分布到变分模型中进行分解,使用希尔伯特变换和高斯平滑,通过寻找约束变分模型的最优解实现复杂信号的分解;求解每一个模态的希尔伯特变换,然后将每一个模态的频谱移到基带上,之后使用解调信号的H高斯平滑,最小化模态带宽的总和:
引入二次惩罚因子α和拉格朗日乘子λ,将约束问题转化为非约束问题,采用交替乘子方向法进行求解:
进一步地,在本发明的实施例中,所述极限学习机(ELM)算法的实现过程如下:
给定一个样本(xi,yi),其中xi=[xi1,x2,...,xiN]T,表示N维输入集,yi=[yi1,yi2,...,yiM]T,表示M维输出集,i=1,2,...,N表示样本标签;极限学习机算法的计算过程如下:
其中,L表示隐含层节点数量,βl表示第l个节点的输出权重矩阵,F为执行运算的激活函数,wl和bl分别表示隐含层第l个节点的输入权重向量和偏置向量;因此,上述运算又可以表述为Hβ=Y,其中β=[β1,β2,...,βL]T,H表示隐含层输出权重矩阵,表示如下:
因此,可以将β表示为:βT=H+T=HT(HHT)-1T,其中H+=HT(HHT)-1表示隐含层输出矩阵H的Moore-Penrose广义逆矩阵。
为验证本发明所述混合模型VMD-ELM算法的优越性,将其与其他一些方法做对比,包括单一预测模型SVR、BPNN和ELM,混合预测模型EMD-SVR、EMD-BPNN和EMD-ELM,ICEEMDAN-SVR、ICEEMDAN-BPNN和ICEEMDAN-ELM,VMD-SVR和VMD-BPNN。所有模型均采用同样的数据集,原始的生猪存栏量数据将直接作为SVR、BPNN和ELM三个单一模型的输入数据集进行模型训练,其余模型均分别使用EMD,ICEEMDAN和VMD方法分解得到的数据来训练。表1给出了本发明提出的混合模型(VMD-ELM)与其他模型在中国生猪月度存栏量数据中的预测误差。
表1不同预测方法的预测误差对比
由表1可知,本发明提出的混合模型在生猪月度存栏量数据上的四项预测误差指标MAE、RMSE、MAPE和TIC均要显著小于其他模型。四种评价标准的计算公式如下:
图3进一步给出了ELM、EMD-ELM、ICEEMDAN-ELM和VMD-ELM之间的预测结果对比示意图。结果表明,单一预测模型ELM的预测结果存在显著的“滞后预测”效应,即对未来一期的预测值容易受到最近一期观测值的影响,从而表现出预测值序列曲线近似于是将当前观测值曲线向后平移一期得到,最终导致了单一预测模型的误差水平相对较大的结果。出现该现象的原因是生猪月度存栏量数据具有典型的非线性、非平稳特征,序列中存在大量噪声,导致模型训练困难,泛化能力相对较差,而在三类混合预测模型中,尤其是本发明所提出的VMD-ELM模型中,该问题得到了有效解决,预测精度也得到了有效的提升。该数据的测试集样本为2019年7月至2021年11月的数据,包含29个月的月度数据,数据来源于中国政府网(http://www.gov.cn)和前瞻数据库(d.qianzhan.com)。其中2021年的数据为官方预测数据。实验结果表明,本发明所提出的VMD-ELM混合模型预测方法在中国生猪月度存栏量序列数据预测中能够取得更高的预测精度。
综上所述,本发明基于变分模态分解(VMD)和极限学习机(ELM)提供了一种自适应预测方法用于对生猪月度存栏量进行精准预测,与现有技术相比,本发明方法能够提高预测效率和预测精度,并在实验中验证了预测的有效性。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种生猪月度存栏量预测方法,其特征在于,包括以下步骤:
通过变分模态分解算法,将生猪存栏量的原始时间序列数据x(t)分解为K个模态;
对于每一个由变分模态分解得到的模态,先进行归一化处理,再按照预设比例划分为训练集和测试集;
采用训练集数据对极限学习机算法进行训练,确定算法的最优参数;
以步长为V的滑动窗口选取极限学习机算法的输入集数据;
对于每一个模态,分别将测试集的输入集数据输入到训练好的极限学习机算法中,输出对下一时点的预测值,将其反归一化处理后得到预测值序列uk(t);
对所有模态的预测值进行相加重构,得到最终的预测值结果;
其中,所述极限学习机算法的实现过程如下:
给定一个样本(xi,yi),其中xi=[xi1,x2,...,xiN]T,表示N维输入集,yi=[yi1,yi2,...,yiM]T,表示M维输出集,i=1,2,...,N表示样本标签;极限学习机算法的计算过程如下:
其中,L表示隐含层节点数量,βl表示第l个节点的输出权重矩阵,F为执行运算的激活函数,wl和bl分别表示隐含层第l个节点的输入权重向量和偏置向量;因此,上述计算过程又可以表述为Hβ=Y,其中β=[β1,β2,...,βL]T,H表示隐含层输出权重矩阵,表示如下:
因此,将β表示为:βT=H+T=HT(HHT)-1T,其中H+=HT(HHT)-1表示隐含层输出矩阵H的Moore-Penrose广义逆矩阵。
4.根据权利要求1所述的生猪月度存栏量预测方法,其特征在于,所述预设比例为8:2。
6.根据权利要求1所述的生猪月度存栏量预测方法,其特征在于,所述变分模态分解算法的过程如下:
初始化模态分量的带宽uk,中心频率ωk以及循环次数n;
当ω大于等于0时,循环更新每个模态分量的带宽uk和中心频率ωk,自适应分解成以中心频率ωk为中心扩散的模态分量:
其中,uk为分解后对模态分量,ωk为模态分量对应的中心频率,λ为拉格朗日乘子,当满足下述条件时,终止该循环:
其中,∈为常数;
在计算每个分量时,将信号分布到变分模型中进行分解,使用希尔伯特变换和高斯平滑,通过寻找约束变分模型的最优解实现复杂信号的分解;求解每一个模态的希尔伯特变换,然后将每一个模态的频谱移到基带上,之后使用解调信号的H高斯平滑,最小化模态带宽的总和:
引入二次惩罚因子α和拉格朗日乘子λ,将约束问题转化为非约束问题,采用交替乘子方向法进行求解:
7.根据权利要求1所述的生猪月度存栏量预测方法,其特征在于,所述方法还包括:
通过比较预测值序列和实际值序列,对预测效果进行验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110339114.9A CN113011669B (zh) | 2021-03-30 | 2021-03-30 | 一种生猪月度存栏量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110339114.9A CN113011669B (zh) | 2021-03-30 | 2021-03-30 | 一种生猪月度存栏量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113011669A CN113011669A (zh) | 2021-06-22 |
CN113011669B true CN113011669B (zh) | 2022-08-02 |
Family
ID=76409205
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110339114.9A Expired - Fee Related CN113011669B (zh) | 2021-03-30 | 2021-03-30 | 一种生猪月度存栏量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113011669B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113764108A (zh) * | 2021-08-16 | 2021-12-07 | 淮阴工学院 | 一种面向covid-19的经验模态分解模糊预测方法 |
CN114037130A (zh) * | 2021-10-28 | 2022-02-11 | 重庆川南环保科技有限公司 | 一种家电维修订单预测方法、装置及终端设备 |
CN115438301B (zh) * | 2022-10-26 | 2023-06-20 | 北京建筑大学 | 一种基于ICEEMDAN二次分解耦合informer模型的设备运行趋势预测方法 |
CN118485186A (zh) * | 2024-07-16 | 2024-08-13 | 安徽省(水利部淮河水利委员会)水利科学研究院(安徽省水利工程质量检测中心站) | 城市用水量数据分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106897794A (zh) * | 2017-01-12 | 2017-06-27 | 长沙理工大学 | 一种基于完备总体经验模态分解和极限学习机的风速预测方法 |
CN109299430A (zh) * | 2018-09-30 | 2019-02-01 | 淮阴工学院 | 基于两阶段分解与极限学习机的短期风速预测方法 |
CN109816236A (zh) * | 2019-01-18 | 2019-05-28 | 中国农业科学院农业信息研究所 | 一种生猪养殖的规划方法及规划系统 |
CN111193256A (zh) * | 2020-01-07 | 2020-05-22 | 上海电力大学 | 基于变分模态分解和门控循环单元的电力负荷预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190188611A1 (en) * | 2017-12-14 | 2019-06-20 | Business Objects Software Limited | Multi-step time series forecasting with residual learning |
-
2021
- 2021-03-30 CN CN202110339114.9A patent/CN113011669B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106897794A (zh) * | 2017-01-12 | 2017-06-27 | 长沙理工大学 | 一种基于完备总体经验模态分解和极限学习机的风速预测方法 |
CN109299430A (zh) * | 2018-09-30 | 2019-02-01 | 淮阴工学院 | 基于两阶段分解与极限学习机的短期风速预测方法 |
CN109816236A (zh) * | 2019-01-18 | 2019-05-28 | 中国农业科学院农业信息研究所 | 一种生猪养殖的规划方法及规划系统 |
CN111193256A (zh) * | 2020-01-07 | 2020-05-22 | 上海电力大学 | 基于变分模态分解和门控循环单元的电力负荷预测方法 |
Non-Patent Citations (1)
Title |
---|
基于HP分离ARIMA—M...v模型的我国生猪存栏量预测;吴雪;《统计与决策》;20130630(第6期);第102-104页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113011669A (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113011669B (zh) | 一种生猪月度存栏量预测方法 | |
Li et al. | Smart wind speed forecasting approach using various boosting algorithms, big multi-step forecasting strategy | |
CN109871977A (zh) | 基于离散小波变换和优化最小支持向量机的负荷预测方法 | |
Huang et al. | A hybrid model for short-term dissolved oxygen content prediction | |
CN112884236B (zh) | 一种基于vdm分解与lstm改进的短期负荷预测方法及系统 | |
CN113836823A (zh) | 一种基于负荷分解和优化双向长短期记忆网络的负荷组合预测方法 | |
CN111353534B (zh) | 一种基于自适应分数阶梯度的图数据类别预测方法 | |
Sun et al. | Spatiotemporal wind power forecasting approach based on multi-factor extraction method and an indirect strategy | |
CN112668611A (zh) | 一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法 | |
CN114970345A (zh) | 短期负荷预测模型构建方法、装置、设备及可读存储介质 | |
Liu et al. | A wind speed forcasting model based on rime optimization based VMD and multi-headed self-attention-LSTM | |
Li et al. | Short-term wind power prediction based on modal reconstruction and CNN-BiLSTM | |
CN113128666A (zh) | 基于Mo-S-LSTMs模型的时间序列多步预测方法 | |
CN115344927A (zh) | 基于工程历史数据的高拱坝施工仿真参数动态更新方法 | |
CN114676645A (zh) | 一种非平稳时间序列预测方法及系统 | |
Sun et al. | Multi-step wind speed forecasting based on a hybrid decomposition technique and an improved back-propagation neural network | |
CN113222234B (zh) | 基于集成模态分解的用气需求预测方法及其系统 | |
Xinxin et al. | Short-term wind speed forecasting based on a hybrid model of ICEEMDAN, MFE, LSTM and informer | |
Wang et al. | A deterministic and probabilistic hybrid model for wind power forecasting based improved feature screening and optimal Gaussian mixed kernel function | |
Ao et al. | Entropy estimation via normalizing flow | |
CN112232565A (zh) | 基于两阶段的时间序列预测方法、预测系统、终端及介质 | |
CN116865255A (zh) | 基于改进熵权法和seceemd的短期风电功率预测方法 | |
Zhang et al. | A Novel Combined Model Based on Hybrid Data Decomposition, MSWOA and ENN for Short-Term Wind Speed Forecasting | |
CN115619563A (zh) | 一种基于神经网络的股票价格分析方法 | |
Ma et al. | Predicting Stock Prices Using Hybrid LSTM and ARIMA Model. |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220802 |