CN108647807A - 河流流量的预测方法 - Google Patents
河流流量的预测方法 Download PDFInfo
- Publication number
- CN108647807A CN108647807A CN201810317991.4A CN201810317991A CN108647807A CN 108647807 A CN108647807 A CN 108647807A CN 201810317991 A CN201810317991 A CN 201810317991A CN 108647807 A CN108647807 A CN 108647807A
- Authority
- CN
- China
- Prior art keywords
- data
- component
- vmd
- input
- indicate
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000013528 artificial neural network Methods 0.000 claims abstract description 40
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 35
- 238000012549 training Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000003416 augmentation Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 8
- 206010068052 Mosaicism Diseases 0.000 abstract 1
- 210000003765 sex chromosome Anatomy 0.000 abstract 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 16
- 238000002474 experimental method Methods 0.000 description 4
- 238000007726 management method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007728 cost analysis Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 210000004218 nerve net Anatomy 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Primary Health Care (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种河流流量的预测方法,包括如下步骤,1:数据预处理;2:VMD模型分解;3:VMD分解结果的分量重构;4:VMD‑BP模型的集成及河流流量预测。变分模态分解是处理非平稳信号的有效方法,结合BP神经网络在处理非线性函数拟合的优势,提出和建立了基于VMD‑BP模型的河流流量预测方法。该方法将原始数据分解为多个反映数据特征的本征模态函数,进行数据的平稳化处理,进而解决水位流量数据的非线性和波动性问题,提高了预测精度。
Description
技术领域
本发明涉及水资源管理技术领域,具体涉及一种河流流量的预测方法。
背景技术
河流流量是水文监测和水资源管理的重要指标,对于水利工程建设、流量计算、航运规划等方面具有重要的指导意义和应用价值。在水利工程效益成本分析中,水文、工程和经济学的结合已成为普遍现象。与水资源系统各组成部分的规划和运作有关的许多活动都需要对未来事件进行预测。合理的流量预测有助于洪水预报及管理、水库蓄水管理、农田灌溉等,将带来积极的经济效益和社会效益。同时,流量的测量技术相对复杂,测量设备比较昂贵,难以实现长期监测,而水位的测量相对简单,因此,通常利用水位流量的拟合曲线,结合水位流量历史资料进行河流流量预测。
河流流量预测主要有三类方法:基于物理过程的模型、基于数理统计的回归模型和基于仿真模拟的智能方法。物理模型法是通过分析河流流动的过程,构建实际模型进行流量模拟预测,但河流流动过程复杂,难以构建合适的物理模型。回归模型法基于河流水文数据构建回归方程进行预测,但回归模型难以处理复杂的水位流量关系。
近年来,由于人工智能方法不必考虑复杂的流动过程,广泛应用于水文预报,包括模型树、支持向量机、基因表达式编程、人工神经网络等。模型树通过将数据分类为多组数据进行线性拟合进行预测,但水位流量并非简单的线性关系,预测精度不高。支持向量机通过内积函数所定义的非线性变换,将输入空间变换到一个高维空间,在高维空间中寻找输出变量和输入变量之间的非线性关系,适合多支支流流入的流域研究,但当训练样本过多会导致运算效率降低。基因表达式编程的公式是透明的,可以表达为一个简单的数学函数,但后期收敛慢,且容易陷入局部最优。
神经网络有很强的非线性拟合能力,并且能够在不解析内部结构的情况下提取出输入和输出之间的关系,适用于不同的水文环境,因此,逐渐应用到河流流量预测。神经网络在河流流量预测的重点主要是输入变量的确定以及模型的优化,如Sivapragasam在密西西比河的研究将降雨量、蒸发量、温度、流量等因素作为输入变量;Yaseen等人将前期流量作为输入变量在马来西亚Johor河流域进行前馈反向传播神经网络和径向基函数神经网络的对比研究;Shiri混合小波变换和模糊系统对神经网络进行优化,在土耳其Filyos河进行了模型的验证;Kalteh结合小波分解优化神经网络,但小波分解存在需要选取小波基、确定分解尺度的问题。
发明内容
本发明的目的在于提供一种河流流量的预测方法,该方法运用变分模态分解将水位、流量数据分解为反映数据特征的本征模态函数,降低数据的波动性,之后再重构分量进行BP(back propagation)神经网络训练预测,叠加得到最终的预测值,实现河流流量的预测。将两种模型结合起来,可以充分发挥各自的特点和优势,为河流流量预测提供更加精确的方法。
为解决上述技术问题,本发明公开的一种河流流量的预测方法,其特征在于,它包括如下步骤:
步骤1:数据预处理,将待预测河流区域的历史水位数据和历史流量数据进行多种随机组合,作为输入条件进行BP神经网络进行BP神经网络训练及预测,确定预测相对误差最小的一组输入数据{u1(t),u2(t)...un(t)}和输出数据o(t),其中un表示第n个输入变量,t表示第t个数据;
步骤2:VMD(变分模态分解)模型分解,利用VMD模型对输入数据{u1(t),u2(t)...un(t)}和输出数据o(t)分别进行平稳化处理,分解为表示输入数据{u1(t),u2(t)...un(t)}输出数据o(t)的多个分量,VMD模型分解的用于寻找输入数据{u1(t),u2(t)...un(t)}的K个分量unk(t)(k=1,2,…,K),使得K个分量unk(t)(k=1,2,…,K)的宽带之和最小;得到输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok},其中unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量;
步骤3:利用VMD模型将输入数据{u1(t),u2(t)...un(t)}和输出数据o(t)分解后,按其对应的分解量进行重构,即将输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok}分别按第一分量、第二分量一直到第k分量的形式进行组合,形成如以下公式所示的k组数据组合,unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量;
步骤4:VMD-BP模型的集成及河流流量预测,将步骤3中VMD分解的分量重构得到的k组数据组合
作为输入条件,分别进行BP神经网络的训练及预测,得到各个分量组合的预测值Dn,之后将每个分量组合的预测值叠加得到最终的待预测河流流量预测值D;D=D1+D2+…+Dn,即实现了基于VMD-BP模型的集成及河流流量预测。
本发明中,BP神经网络是一种按照误差逆向传播算法训练的多层前馈神经网络,具有很强的非线性映射能力,在处理水位流量的非线性拟合上有很大的优势。
本发明中,VMD模型是一种完全非递归的方法,可以将信号分解为一定数量的有限宽带的本征模态函数(BLIMF),每个BLIMF包含原信号的特征信息。它克服了离散小波变换、经验模态分解、局域均值分解等模式中存在的模态混合、不能正确地消除附加噪声和确定本征模态函数(IMF)个数的问题,适合处理具有波动性的水位流量数据。该方法融合VMD模型和BP神经网络,建立VMD-BP预测模型,以解决基于水位流量多值关系的河流流量预测问题。
本发明的有益效果为:
1、将前期水位、前期流量作为模型的输入条件,有效的解决了水位流量之间复杂的绳套曲线关系。
2、VMD模型是一种处理非线性、非平稳信号的谐波分析方法,在处理具有波动性的水位流量数据上有很大的优势,BP神经网络具有很强的非线性拟合能力,结合两种模型的优势,提高了河流流量预测的精度。
附图说明
图1为本发明的流程图;
图2为输出数据VMD模型分解结果。
图2中,横轴表示样本序列。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的详细说明:
本发明的一种河流流量的预测方法,如图1所示,它包括如下步骤:
步骤1:数据预处理,将待预测河流区域的历史水位数据和历史流量数据进行多种随机组合,作为输入条件进行BP神经网络进行BP神经网络训练及预测,确定预测相对误差最小的一组输入数据{u1(t),u2(t)...un(t)}和输出数据o(t),其中un表示第n个输入变量,t表示第t个数据;
步骤2:VMD模型分解,利用VMD模型对输入数据{u1(t),u2(t)...un(t)}和输出数据o(t)分别进行平稳化处理,分解为表示输入数据{u1(t),u2(t)...un(t)}输出数据o(t)的多个分量,VMD模型分解的用于寻找输入数据{u1(t),u2(t)...un(t)}的K个分量unk(t)(k=1,2,…,K),使得K个分量unk(t)(k=1,2,…,K)的宽带之和最小;得到输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok},其中unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量;
步骤3:利用VMD模型将输入数据{u1(t),u2(t)...un(t)}和输出数据o(t)分解后,按其对应的分解量进行重构,即将输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok}分别按第一分量、第二分量一直到第k分量的形式进行组合,形成如以下公式所示的k组数据组合,unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量;
步骤4:VMD-BP模型的集成及河流流量预测,将步骤3中VMD分解的分量重构得到的k组数据组合
作为输入条件,分别进行BP神经网络的训练及预测,得到各个分量组合的预测值Dn,之后将每个分量组合的预测值叠加得到最终的待预测河流流量预测值D;D=D1+D2+…+Dn,即实现了基于VMD-BP模型的集成及河流流量预测。
上述技术方案的步骤2中,得到输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok}的具体方法为:
步骤201:首先对每个输入数据分量unk(t)通过希尔伯特变换
计算相关联的分析信号从而获得单边频谱:
步骤202:通过将得到的所述分析信号与希尔伯特变换得到的分析信号的指数调谐项相乘得到解调信号,将每个有限宽带的本征模态函数(BLIMF)的频谱移位到基带,再根据解调信号的梯度的二范数平方估计每个有限宽带的本征模态函数的宽带,最后得到变量unk(t)的约束变分公式:
式中,{unk}:={un1,...,unk}是第n个输入数据分解的k个分量的集合,{ωnk}:={ωn1,...,ωnk}是第n个输入数据分解的分量对应的中心频率的集合,表示所有分量的和,f(t)是原始数据,即输入数据un,表示对t求偏导数,t表示分析信号的个数,δ(t)表示脉冲函数,表示相位,表示希尔伯特变换得到的分析信号的指数谐调项;
步骤203:引入二次惩罚项α和拉格朗日乘子λn(t)计算得到增广拉格朗日公式:
上述的目的是将约束变分公式变为无约束问题,方便求解,式中,二次惩罚项α为足够大的正数,用于保证信号的重构精度,拉格朗日乘子λn(t)使约束条件保持严格性,表示对t求偏导数;
步骤204:利用交替方向乘子法,通过迭代更新未知的值,直到,迭代更新之后得到的值的精度E满足E<ε停止迭代,ε表示给定的判别精度,ε>0,比如10-7;
式中,m表示迭代次数,ω为输入信号的中心频率, 分别为对应的傅立叶变换,表示第n个输入变量的第k个分量迭代m次之后的值,λn(t)表示拉格朗日乘子,α表示二次惩罚项,ωnk表示第n个输入数据分解的第k个分量对应的中心频率;
步骤205:最终得到输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok},其中unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量。
上述技术方案的步骤4中得到最终的待预测河流流量预测值D的方法为:
步骤401:对重构的分量组合进行归一化处理,消除单位差异和异常数据的影响,同时归一化处理可加快网络的收敛,减少运行时间,因为进行VMD分解后得到的有限宽带的本征模态函数中既有正数也有负数,因此将分解后得到的有限宽带的本征模态函数归一化映射到[-1,1]
y=2*(x-min)/(max-min)-1
式中,x是重构的分量组合,min、max分别为重构的分量组合的最小值和最大值,y为归一化输出;
步骤402:第n个输入变量的第k个分量unk进入BP神经网络的输入层,累加unk与神经网络输入层和隐藏层之间的权值wnh、神经网络输入层和隐藏层之间的偏置值θn的乘积传入BP神经网络的隐藏层,经过BP神经网络的传递函数fS(x)的计算后得到隐藏层输出,之后再累计隐藏层输出与神经网络网络隐藏层与输出层之间权值wh、神经网络网络隐藏层与输出层之间偏置值θh的乘积后传入BP神经网络的输出层,得到输出信号yk;所述传递函数fS(x)为双极S型函数;
式中,x为输入隐藏层的数据,e为自然对数;
步骤403:根据神经网络在训练的过程中产生的预测输出信号yk和对应的期望输出Ok计算得到误差E,判断误差是否小于给定值ε,如果大于则将该误差反向传播给隐藏层求其误差梯度,更新权值和偏置值;并返回正向传播中进行计算,直到误差小于期望误差或者达到最大迭代次数,得到最终训练的网络。
步骤404:重复上述步骤401~403分别得到每个分量组合训练后的网络;
步骤405:根据步骤1中确定的输入数据方案确定河流流量预测数据的输入数据,输入对应的训练后的网络,得到每个组合的预测值Dn;
步骤406:将上述得到的预测值叠加得到最终的预测值;
D=D1+D2+…+Dn。
本实施例以长江螺山观测站1998、1999两年的日水位和日流量共604个数据作为研究数据,随机选取60个(约10%)数据作为测试数据,剩余的数据作为训练数据。
采用相对误差和确定性系数DC来衡量预测值和实际值的吻合程度,其中Qi表示实测值,表示预测值,表示实测平均值。较好的预测结果应该是DC>0.9,且相对误差小于5%。
1数据预处理
因为螺山水文站受到上游洪水涨落及下游回水顶托等的影响,水位流量存在多值关系,呈现复杂的绳套曲线,单独使用水位来描述流量的效果不佳,因此采用多输入的方式用前后几天的水位流量数据来反映水势,避免多值关系对于拟合精度的影响。为选取最佳的输入方案,进行表1中各个方案的测试。其中输入数据包含Ht(当天水位)、Ht-2(前两天水位)、Qt-2(前两天流量)、Ht-1(前一天水位)、Qt-1(前一天流量)、Ht+1(后一天水位)等变量;输出数据为Qt(当天流量)。
表1输入方案及其相对误差
表1输入方案及其相对误差
根据表1各方案相对误差的比较,确定模型最终的输入条件为Ht、Ht-1、Qt-1、Ht+1,输出条件为Qt。
2VMD模型分解
将输入数据及输出数据进行VMD分解,其中输出数据的分解结果如附图2。为避免过多分解带来的误差以及减少后期BP神经网络的构建量,VMD分解中选择将数据分解为三个分量,分别为反映原始数据基本趋势的趋势量、频率波动较低的低频量以及频率波动较大的高频量。惩罚参数C采用VMD默认值2 000;τ选取0.3以保证数据分解的保真度。
3VMD分解结果的分量重构
将四个输入数据及一个输出数据的VMD分解结果按照趋势量、低频量、高频量进行分量组合,形成三组组合。
4VMD-BP模型的集成及河流流量预测
水位、流量数据的计量单位以及数量级不相同,有可能由于其数据范围的不同导致其对于最终结果的影响程度不同,从而导致结果预测误差增大。因此,对重构的分量组合进行归一化处理,消除单位差异和异常数据的影响。同时归一化处理可加快网络的收敛,减少运行时间。因为进行VMD分解后得到的BLIMF中既有正数也有负数,因此将数据归一化映射到[-1,1]。
将归一化后的三组组合分别进行BP神经网络的训练及预测,分别得到趋势量、低频量、高频量的预测值,然后将三者的结果叠加得到最终的预测值,预测结果如表2。
表2 VMD各分量误差
为更加具体的描述VMD-BP模型的预测精度,进行三个实验,单输入单输出BP神经网络、多输入单输出BP神经网络以及VMD-BP模型。三个实验的预测结果如表3。
表3三个模型误差对比表
第一个和第二个实验的结果对比体现多输入的BP神经网络可以通过输入多个变量描述水势来有效的解决水位流量的多值关系;比较第二个和第三个实验的结果反映了VMD-BP模型的有效性和精确性。VMD-BP模型在三个实验中相对误差最小,且确定性系数最接近于1,是更好的河流流量预测方法。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (5)
1.一种河流流量的预测方法,其特征在于,它包括如下步骤:
步骤1:将待预测河流区域的历史水位数据和历史流量数据进行多种随机组合,作为输入条件进行BP神经网络进行BP神经网络训练及预测,确定预测相对误差最小的一组输入数据{u1(t),u2(t)...un(t)}和输出数据o(t),其中un表示第n个输入变量,t表示第t个数据;
步骤2:利用VMD模型对输入数据{u1(t),u2(t)...un(t)}和输出数据o(t)分别进行平稳化处理,分解为表示输入数据{u1(t),u2(t)...un(t)}输出数据o(t)的多个分量,VMD模型分解的用于寻找输入数据{u1(t),u2(t)...un(t)}的K个分量unk(t)(k=1,2,…,K),使得K个分量unk(t)(k=1,2,…,K)的宽带之和最小;得到输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok},其中unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量;
步骤3:利用VMD模型将输入数据{u1(t),u2(t)...un(t)}和输出数据o(t)分解后,按其对应的分解量进行重构,即将输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok}分别按第一分量、第二分量一直到第k分量的形式进行组合,形成如以下公式所示的k组数据组合,unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量;
步骤4:将步骤3中VMD分解的分量重构得到的k组数据组合
作为输入条件,分别进行BP神经网络的训练及预测,得到各个分量组合的预测值Dn,之后将每个分量组合的预测值叠加得到最终的待预测河流流量预测值D;D=D1+D2+…+Dn。
2.根据权利要求1所述的河流流量的预测方法,其特征在于:
步骤2中,得到输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok}的具体方法为:
步骤201:首先对每个输入数据分量unk(t)通过希尔伯特变换
计算相关联的分析信号从而获得单边频谱:
步骤202:通过将得到的所述分析信号与希尔伯特变换得到的分析信号的指数调谐项相乘得到解调信号,将每个有限宽带的本征模态函数的频谱移位到基带,再根据解调信号的梯度的二范数平方估计每个有限宽带的本征模态函数的宽带,最后得到变量unk(t)的约束变分公式:
式中,{unk}:={un1,...,unk}是第n个输入数据分解的k个分量的集合,{ωnk}:={ωn1,...,ωnk}是第n个输入数据分解的分量对应的中心频率的集合,表示所有分量的和,f(t)是原始数据,即输入数据un,表示对t求偏导数,t表示分析信号的个数,δ(t)表示脉冲函数,表示相位,表示希尔伯特变换得到的分析信号的指数谐调项;
步骤203:引入二次惩罚项α和拉格朗日乘子λn(t)计算得到增广拉格朗日公式:
上述的目的是将约束变分公式变为无约束问题,方便求解,式中,二次惩罚项α为足够大的正数,用于保证信号的重构精度,拉格朗日乘子λn(t)使约束条件保持严格性,表示对t求偏导数;
步骤204:利用交替方向乘子法,通过迭代更新未知的值,直到,迭代更新之后得到的值的精度E满足E<ε停止迭代,ε表示给定的判别精度;
式中,m表示迭代次数,ω为输入信号的中心频率, 分别为λn(t)、f(t)、uni(t)、对应的傅立叶变换,表示第n个输入变量的第k个分量迭代m次之后的值,λn(t)表示拉格朗日乘子,α表示二次惩罚项,ωnk表示第n个输入数据分解的第k个分量对应的中心频率;
步骤205:最终得到输入数据VMD分解的分量组合{un1,un2,...,unk}及输出数据VMD分解的分量组合{o1,o2,...,ok},其中unk表示第n个输入变量的第k个分量,ok表示输出变量的第k个分量。
3.根据权利要求1所述的河流流量的预测方法,其特征在于:步骤4中得到最终的待预测河流流量预测值D的方法为:
步骤401:对重构的分量组合进行归一化处理,消除单位差异和异常数据的影响,同时归一化处理可加快网络的收敛,减少运行时间,因为进行VMD分解后得到的有限宽带的本征模态函数中既有正数也有负数,因此将分解后得到的有限宽带的本征模态函数归一化映射到[-1,1]
y=2*(x-min)/(max-min)-1
式中,x是重构的分量组合,min、max分别为重构的分量组合的最小值和最大值,y为归一化输出;
步骤402:第n个输入变量的第k个分量unk进入BP神经网络的输入层,累加unk与神经网络输入层和隐藏层之间的权值wnh、神经网络输入层和隐藏层之间的偏置值θn的乘积传入BP神经网络的隐藏层,经过BP神经网络的传递函数fS(x)的计算后得到隐藏层输出,之后再累计隐藏层输出与神经网络网络隐藏层与输出层之间权值wh、神经网络网络隐藏层与输出层之间偏置值θh的乘积后传入BP神经网络的输出层,得到输出信号yk;
步骤403:根据神经网络在训练的过程中产生的预测输出信号yk和对应的期望输出Ok计算得到误差E,判断误差是否小于给定值ε,如果大于则将该误差反向传播给隐藏层求其误差梯度,更新权值和偏置值;并返回正向传播中进行计算,直到误差小于期望误差或者达到最大迭代次数,得到最终训练的网络。
4.根据权利要求3所述的权利要求1所述的河流流量的预测方法,其特征在于:所述步骤403后还步骤404:重复上述步骤401~403分别得到每个分量组合训练后的网络;
步骤405:根据步骤1中确定的输入数据方案确定河流流量预测数据的输入数据,输入对应的训练后的网络,得到每个组合的预测值Dn;
步骤406:将上述得到的预测值叠加得到最终的预测值;
D=D1+D2+…+Dn。
5.根据权利要求3所述的权利要求1所述的河流流量的预测方法,其特征在于:所述BP神经网络的传递函数fS(x)为双极S型函数;
式中,x为输入隐藏层的数据,e为自然对数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810317991.4A CN108647807B (zh) | 2018-04-10 | 2018-04-10 | 河流流量的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810317991.4A CN108647807B (zh) | 2018-04-10 | 2018-04-10 | 河流流量的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108647807A true CN108647807A (zh) | 2018-10-12 |
CN108647807B CN108647807B (zh) | 2021-08-13 |
Family
ID=63745790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810317991.4A Expired - Fee Related CN108647807B (zh) | 2018-04-10 | 2018-04-10 | 河流流量的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108647807B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110921A (zh) * | 2019-04-30 | 2019-08-09 | 武汉理工大学 | 一种考虑时滞效应的河流水位预测方法 |
CN110197306A (zh) * | 2019-05-31 | 2019-09-03 | 东北大学 | 一种基于vmd-pacf-bp模型的动态血糖在线预测方法 |
CN111753965A (zh) * | 2020-06-30 | 2020-10-09 | 长江水利委员会水文局 | 一种基于深度学习的河流流量自动整编方法及系统 |
CN111860970A (zh) * | 2020-06-29 | 2020-10-30 | 浪潮电子信息产业股份有限公司 | 一种河流流量预测方法、装置及电子设备和存储介质 |
CN115687308A (zh) * | 2022-12-08 | 2023-02-03 | 长江水利委员会水文局 | 一种水位流量关系智能定线方法 |
CN117744884A (zh) * | 2023-12-29 | 2024-03-22 | 南方电网调峰调频发电有限公司鲁布革水力发电厂 | 水库水流量预测模型的构建方法和水库水流量预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846753A (zh) * | 2010-04-29 | 2010-09-29 | 南京信息工程大学 | 基于经验模态分解和支持向量机的气候时序预测方法 |
CN102147273A (zh) * | 2010-01-29 | 2011-08-10 | 大连理工大学 | 一种基于数据的冶金企业高炉煤气动态预测方法 |
CN105978037A (zh) * | 2016-08-03 | 2016-09-28 | 河海大学 | 基于风速预测的多时段电气互联系统最优潮流计算方法 |
CN106777775A (zh) * | 2017-01-10 | 2017-05-31 | 清华大学 | 一种基于多断面水位预测河流流量的神经网络方法 |
CN106845641A (zh) * | 2017-03-03 | 2017-06-13 | 东南大学 | 基于经验模态分解和bp神经网络的地铁沉降预测方法 |
CN107315892A (zh) * | 2017-08-10 | 2017-11-03 | 北京交通大学 | 一种基于极限学习机的轴承故障诊断方法 |
-
2018
- 2018-04-10 CN CN201810317991.4A patent/CN108647807B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102147273A (zh) * | 2010-01-29 | 2011-08-10 | 大连理工大学 | 一种基于数据的冶金企业高炉煤气动态预测方法 |
CN101846753A (zh) * | 2010-04-29 | 2010-09-29 | 南京信息工程大学 | 基于经验模态分解和支持向量机的气候时序预测方法 |
CN105978037A (zh) * | 2016-08-03 | 2016-09-28 | 河海大学 | 基于风速预测的多时段电气互联系统最优潮流计算方法 |
CN106777775A (zh) * | 2017-01-10 | 2017-05-31 | 清华大学 | 一种基于多断面水位预测河流流量的神经网络方法 |
CN106845641A (zh) * | 2017-03-03 | 2017-06-13 | 东南大学 | 基于经验模态分解和bp神经网络的地铁沉降预测方法 |
CN107315892A (zh) * | 2017-08-10 | 2017-11-03 | 北京交通大学 | 一种基于极限学习机的轴承故障诊断方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110921A (zh) * | 2019-04-30 | 2019-08-09 | 武汉理工大学 | 一种考虑时滞效应的河流水位预测方法 |
CN110197306A (zh) * | 2019-05-31 | 2019-09-03 | 东北大学 | 一种基于vmd-pacf-bp模型的动态血糖在线预测方法 |
CN111860970A (zh) * | 2020-06-29 | 2020-10-30 | 浪潮电子信息产业股份有限公司 | 一种河流流量预测方法、装置及电子设备和存储介质 |
CN111860970B (zh) * | 2020-06-29 | 2022-03-22 | 浪潮电子信息产业股份有限公司 | 一种河流流量预测方法、装置及电子设备和存储介质 |
CN111753965A (zh) * | 2020-06-30 | 2020-10-09 | 长江水利委员会水文局 | 一种基于深度学习的河流流量自动整编方法及系统 |
CN115687308A (zh) * | 2022-12-08 | 2023-02-03 | 长江水利委员会水文局 | 一种水位流量关系智能定线方法 |
CN117744884A (zh) * | 2023-12-29 | 2024-03-22 | 南方电网调峰调频发电有限公司鲁布革水力发电厂 | 水库水流量预测模型的构建方法和水库水流量预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108647807B (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108647807A (zh) | 河流流量的预测方法 | |
Olyaie et al. | A comparison of various artificial intelligence approaches performance for estimating suspended sediment load of river systems: a case study in United States | |
Kisi et al. | Precipitation forecasting using wavelet-genetic programming and wavelet-neuro-fuzzy conjunction models | |
Talebizadeh et al. | Uncertainty analysis for the forecast of lake level fluctuations using ensembles of ANN and ANFIS models | |
Yu et al. | Forecasting crude oil price with an EMD-based neural network ensemble learning paradigm | |
CN109886498A (zh) | 一种基于特征选择的emd-gru短期电力负荷预测方法 | |
CN108038580A (zh) | 基于同步挤压小波变换的光伏功率多模型综合预测方法 | |
CN107886161A (zh) | 一种提高复杂信息系统效能的全局敏感性分析方法 | |
CN110334875A (zh) | 计及评估指标冲突的风电功率组合概率预测方法 | |
CN109886464A (zh) | 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法 | |
Han et al. | Network traffic prediction using variational mode decomposition and multi-reservoirs echo state network | |
CN103268525A (zh) | 一种基于wd-rbf的水文时间序列模拟预测方法 | |
CN110490366A (zh) | 基于变分模态分解和迭代决策树的径流量预测方法 | |
Li et al. | Decomposition-ANN methods for long-term discharge prediction based on Fisher’s ordered clustering with MESA | |
Wu et al. | Local and global Bayesian network based model for flood prediction | |
Rajaee et al. | Utilization of WGEP and WDT models by wavelet denoising to predict water quality parameters in rivers | |
Mazandaranizadeh et al. | Development of a PSO-ANN model for rainfall-runoff response in basins, Case Study: Karaj Basin | |
CN116822743A (zh) | 一种基于两阶段分解重构与误差校正的风电功率预测方法 | |
CN110110921A (zh) | 一种考虑时滞效应的河流水位预测方法 | |
Guo et al. | Groundwater depth forecasting using configurational entropy spectral analyses with the optimal input | |
Dalavi et al. | Modeling runoff in Bhima River catchment, India: A comparison of artificial neural networks and empirical models | |
CN116628442B (zh) | 一种基于人工神经网络的地下水储量变化空间降尺度方法 | |
CN107704944A (zh) | 一种基于信息论学习的股市波动区间预测方法 | |
Dayal et al. | Application of hybrid artificial neural network algorithm for the prediction of standardized precipitation index | |
Tao et al. | A combined model of wavelet and neural network for short term load forecasting |
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 |
Granted publication date: 20210813 |
|
CF01 | Termination of patent right due to non-payment of annual fee |