技术摘要：
本发明公开了一种基于二次分解与回声状态网络（ESN）的中长期径流预测方法，属于径流预测技术领域。该方法包括：步骤一：获取径流序列x(t)，根据径流序列的数据情况，将其分为训练样本和测试样本；步骤二：利用自适应噪声完备经验模态分解（CEEMDAN）将径流序列分解为  全部
背景技术：
目前径流预测方法使用较为广泛的是主要以人工神经网络为代表的人工智能方 法，但传统的递归神经网络具有记忆渐消的缺点。因此，具有良好短期记忆功能、简单训练 算法以及容易确定网络结构的ESN获得越来越多研究人员的关注。然而，径流序列的复杂性 影响了ESN模型的参数率定及学习效率，单一模型不能充分挖掘并学习径流序列中所包含 的全部信息，预测精度提升空间有限，从而出现了“分解-预测-重构”模式的组合预测方法， 通过将径流序列分解成若干子序列的方式，简化原始径流序列。但这种模式仍有一定缺陷， 分解出的最高频分量包含了原始径流序列的主要特征，具有波动强烈且振动幅值大的特 点，序列仍较为复杂，导致最高频分量的预测精度不足，进而对整个径流序列预测效果产生 较大的影响。 因此，寻求新的组合模式来减小最高频分量预测误差，以期提高整体预测精度，已 成为水文预报领域的热点问题。
技术实现要素：
本发明旨在提供一种基于二次分解与回声状态网络(ESN)的中长期径流预测方 法，解决了一次分解中最高频分量仍波动强烈的问题，进一步提高了性能较好的模型预测 精度，对中长期径流预测具有很好的参考价值。 本发明提供了一种基于二次分解与回声状态网络(ESN)的中长期径流预测方法， 将CEEMDAN-VMD与ESN模型组合使用，首先通过CEEMDAN方法，将径流序列分解出一系列包含 不同尺度局部信息的子序列，以解决径流的复杂性影响模型参数确定及预测效率；再采用 VMD对子序列中频率最高的分量进行二次频域特性分解，降低最高频分量的波动程度，在此 基础上，采用预测性能较好的ESN模型，从而综合提高预测精度，对中长期径流预测具有借 鉴意义。 上述内容使用CEEMDAN可有效地将复杂的径流序列分解为包含不同尺度时域信息 的子序列，可避免EMD分解出现的模态混叠问题，与集合经验模态分解(EEMD)相比，提高了 分解效率，同时解决互补集合经验模态分解(CEEMD)添加不同噪声信号产生不同数量模态 分解的问题。该方法实现了分解信号的近似完美重构，解决了不同加噪信号产生不同数量 模态的缺陷；然而其频率最高的IMF分量包含主要的原始序列频率信息，波动较为强烈，序 列仍较为复杂，导致预测精度较低，借鉴分解方法对原始径流序列的处理思路，采用VMD方 法对最高频分量进行频域特性分解，进一步提取最高频分量的不同频域信息，降低序列复 杂程度，以期进一步提高ESN模型预测精度；而使用的ESN模型，其作为传统循环神经网络 (RNN)的改进方法，具有独特的动态储备池结构，通过调整储备池中神经元状态的数值特 5 CN 111553513 A 说　明　书 2/7 页 征，使得系统具有良好的短期记忆功能及非线性逼近能力，且简化了网络训练过程，提高计 算效率。 上述基于二次分解与回声状态网络的中长期径流预测方法，具体包括以下步骤： 步骤一：获取径流序列x(t)，根据径流序列的数据情况，将其分为训练样本和测试 样本； 步骤二：利用CEEMDAN将径流序列分解为若干IMF分量和一个Res； 步骤三：使用VMD方法对频率最高的IMF分量进行二次分解，得到若干Mode； 步骤四：将两次分解出的子序列分别输入ESN进行预测； 步骤五：将各子序列的预测结果重构，得到最终径流序列预测值。 下面对各步骤具体说明如下： 步骤二所述利用CEEMDAN将径流序列分解为若干IMF分量和一个Res，具体包括以 下步骤： (1)对原始径流序列x(t)添加服从标准正态分布的白噪声序列： xm(t)＝x(t) ε0σm(t)m＝1,2,…,M 式中，xm(t)为加入白噪声后的序列；ε0为适用于xm(t)的信噪比；M为集总次数，即 实现M次经验模态分解EMD；σm(t)为服从标准正态分布的白噪声； (2)对时间序列xm(t)进行M次EMD分解，并计算总体平均值，将其定义为原始信号 的第一个IMF分量，即： 式中，cm为EMD分解得到的第m个IMF分量；c1为第一个IMF分量； (3)计算一阶残差r1(t)： r1(t)＝x(t)-c1(t) (4)对一阶残差添加白噪声后再次进行EMD分解，求取平均值得到第二个IMF分量 c2(t)： 式中，ε1为适用于r1(t)的信噪比；E1(·)为EMD算法产生第1个IMF分量的算子； (5)对后面的每个阶段，k＝2,3,…,K得到的第k个余量信号rk(t)与第(3)步相同， 再重复第(4)步计算第k 1个IMF分量，即： rk(t)＝rk-1(t)-ck(t) 式中，εk为适用于rk(t)的信噪比；Ek(·)为EMD算法产生第k 1个IMF分量的算子； (6)直至残差不能再分解，从而求得最终残差R(t)： (7)则原始序列x(t)的表达式为： 6 CN 111553513 A 说　明　书 3/7 页 式中，K代表在给定精度下原始序列分解出的IMF分量总个数。 步骤三所述使用VMD方法对频率最高的IMF分量进行二次分解，得到若干Mode，具 体包括以下步骤： (1)采用中心频率观察法确定VMD分解出的变分模态Mode的个数r，令r＝1,2,…, R，采用VMD对IMF中的最高频分量进行二次分解，若当r＝R 1时，VMD分解出的Mode中心频率 接近，则r＝R为最终选取的Mode个数； (2)对于IMF分量中的最高频分量采用Hilbert变换，得到原始Mode的解析信号，并 得到单边频谱： 式中，δ(t)为单位脉冲函数；j表示虚数单位；ur(t)为第r个Mode；*表示卷积； (3)将各个Mode的频谱调解到基带，即对 乘以 来调整整个ur (t)的估算中心频率，将各个Mode调解到基频带； (4)利用约束变分模型求解每个Mode分量的带宽： 同时满足： 式中，ωk为各个Mode的角频率；c1(t)为被二次分解的最高频分量； 表示对t求偏 导； (5)通过引入二次惩罚因子α和拉格朗日乘法算子λ(t)求解约束变分模型的最优 解，即可得到最优Mode；其中二次惩罚因子α能保证信号的重构精度，拉格朗日乘法算子λ (t)能加强约束。 步骤四所述将两次分解出的子序列分别输入ESN进行预测，具体包括以下步骤：将 步骤二、步骤三分解出的各子序列，分为与原始径流序列长度相同的训练样本及测试样本， 再将各子序列训练样本输入ESN进行模型训练，寻找最优参数，最后采用建立好的模型预测 各子序列的测试样本。 进一步地，所述的ESN建立过程为： (1)训练阶段 ①初始化ESN网络： 对输入连接权值Win、反馈连接权值Wback和内部状态向量W进行初始化，设定储备池 的内部各个参数，构建动态储备池，并更新内部状态向量： x(n)＝f(Winu(n) Wx(n-1) Wbacky(n))n＝2,3,…,N 7 CN 111553513 A 说　明　书 4/7 页 式中，u(n)＝u1(n) ,u2(n) ,…,uP(n)为ESN的输入向量；x(n)＝x1(n) ,x2(n) ,…,xN (n)为ESN的状态向量；y(n)＝y1(n) ,y2(n) ,…,yL(n)为ESN的输出向量；P为输入向量个数；N 为储备池规模；L为输出向量个数；f为储备池处理单元的内部神经元激活函数，常用tanh函 数； ②收集内部状态向量B与相应输出信号T，以此求得输出权重Wout： Wout＝B-1T ③根据求出的输出权重Wout建立高效的ESN预测模型； (2)测试阶段根据训练阶段得到网络输出权重Wout求得最终预测结果： y(n)＝fout(Woutu(n) ,x(n) ,y(n-1))n＝2,3,…,N 式中，y(n)＝y1(n) ,y2(n) ,…,yL(n)为输出向量；L为输入向量个数；fout为输出处 理单元的内部神经元的激活函数，常用tanh函数。 本发明将CEEMDAN-VMD与ESN模型组合使用，带来的有益效果： (1)使用CEEMDAN可有效地将复杂的径流序列分解为包含不同尺度时域信息的子 序列，可避免EMD分解出现的模态混叠问题，与集合经验模态分解(EEMD)相比，提高了分解 效率，同时解决互补集合经验模态分解(CEEMD)添加不同噪声信号产生不同数量模态分解 的问题。 (2)采用VMD方法对最高频分量进行频域特性分解，进一步提取最高频分量的不同 频域信息，降低序列复杂程度，以期进一步提高ESN模型预测精度； (3)使用的ESN模型，其作为传统循环神经网络(RNN)的改进方法，具有独特的动态 储备池结构，通过调整储备池中神经元状态的数值特征，使得系统具有良好的短期记忆功 能及非线性逼近能力，且简化了网络训练过程，提高计算效率。 附图说明 图1为本发明的流程图； 图2为ESN的结构图。