Page 45 - 2024年第55卷第9期
P. 45

图 2 BiLSTM网络结构

                                                i = σ (W ·[h ,x] + b)
                                                 t     i    t - 1  t  i
                                                c=        i    t - 1  t  c                             (11)
                                                珓 tanh (W ·[h ,x] + b)
                                                 t
                                                              c
                                                c= f·c + i· 珓
                                                 t  t   t - 1  t  t
              式中:W 、i、b分别为输入门权重、遗忘率与偏置;c 与 c分别为前一状态与当前状态下的隐藏层
                      i   t  i                                   t - 1  t
                           c为中间细胞状态值;W 和 b 为细胞状态层的权重和偏置。
              细胞状态值; 珓                            c   c
                            t
                  ( 3)通过运算 “输出门” 的 Sigmoid函数与 tanh函数得到细胞状态 c与 o并向下传递,生成最终
                                                                                 t   t
              模型的输出结果,计算公式如式(12)(13)所示:
                                                o= σ (W ·[h ,x] + b)                                   (12)
                                                  t     o     t - 1  t  o
                                                     h= o·tanh(c)                                      (13)
                                                          t
                                                      t
                                                                  t
              式中:W 为输出门权重矩阵;b 为输出门偏置。
                      o                     o
                  ( 4)由前一时刻前向层输出 h、反向层输出 r组合叠加为当前时刻 BiLSTM网络每一级的计算结
                                             t              t
              果,计算公式如式( 14)所示:
                                                    Y= g(W h+ W r)                                     (14)
                                                      t     q t   h t
              式中:Y 为 BiLSTM当前时刻的输出;g为输出层的激活函数;W 与 W 分别为前向层与反向层的输
                      t
                                                                                h
                                                                          q
              出映射至输出层的权重矩阵。
              2.4 大坝变形预测模型构建
              2.4.1 模型构建 基于上述理论,构建基于自适应时序分解筛选的大坝变形预测模型,建模实现流程
              如图 3所示,主要步骤如下:
                  ( 1)利用 GWO对 CEEMDAN分解参数中的分解层数 max_imf、添加次数 trials及噪声幅值 epsilon
              进行寻优,利用皮尔逊相关系数法与相关系数原理确定噪声分量并进行相应阈值降噪重构,以含噪信
              号与降噪误差比为每头狼的个体适应度。
                  ( 2)以寻优所得最佳参数组合对原始数据序列进行 CEEMDAN阈值降噪。
                  (3)利用交叉验证的递归特征选择法在原始数据集中寻找最优特征集合,选取广泛使用且效果较
              好的回归方法作为元学习器,对所有可能的特征数量进行遍历;同时以 k折交叉试验所得评价指标的
              均值作为判断依据来选取最优学习器与最优特征集合。
                  ( 4)分析滤波降噪后变形序列特性,设置 BiLSTM模型与 Adam优化算法参数构建预测模型,以经
              递归特征选择法筛选后的最优特征集合为输入因子,通过迭代运行完成构建 BiLSTM预测模型。
              2.4.2 模型评价指标 以拟合的良好性与检验的有效性等为评价原则                             [26] ,选取以下指标进行模型评价:

                                          1  m                      1  m
                                                        2
                                RMSE=       ∑  (Y -Y′) ,     MAE=     ∑   Y -Y′
                                       槡                            m i =1
                                                 i
                                                                                i
                                                                           i
                                                      i
                                          m i =1
                                         m
                                        ∑  (Y -Y′)  2                   m
                                                  i
                                              i
                                                                             i
                                                                                 i
                                 2
                                R =1-   i =1         ,       MAPE=   1 ∑   Y -Y′   ×100%
                                         m
                                                                     m i =1   Y i
                                                    2
                                        ∑  ( 珔    i
                                             Y-Y′)
                                         i =1
              式中:Y 为实测值;Y′为预测值; 珔
                                              Y为实测数据序列均值。
                      i
                                  i
                —  1 0 0 —
                     5
   40   41   42   43   44   45   46   47   48   49   50