Page 117 - 2022年第53卷第9期
P. 117

S - 1
                                                                 T
                                                      =
                                                  Z S - 1 ∑  β i i ω i  i
                                                            g( W+ b)
                                                        i =1                                            (8)
                                                  S = 2,3,4,…,Z = 0
                                                                    0
                                           T
                                                                                              i
              式中:β i  = [ β i,1 ,β i,2 ,…,β i,d ]为第 i个隐含节点的输出权重,d为输出维度;ω i                和 b分别为第 i个
              隐含节点的输入权重和偏置;g(·)为 SCN模型的激活函数。同时,根据式(9)可计算当前网络的误
              差 e :
                  S - 1
                                         e = [Z - Z ] = [e    ,e    ,…,e       ]                        (9)
                                          S - 1    S - 1  S - 1,1  S - 1,2  S - 1,d
                  ( 3)SCN引入监督机制为隐含节点 S分配参数,具体监督机制形式如下:
                                            T             T                  T       T
                                   g= [g( ω S w + b),g( ω S w + b),…,g( ω S  w + b)]                   (10)
                                    S    S     1  S    S     2  S        S     n  S
                                               2
                                                   2
                                    〈e    ,g〉≥b(1 - r - μ S ‖e      2                                  (11)
                                                           )
                                                               S - 1 ,j ‖ ,j = 1 ,2,…,d
                                       S - 1 ,j  S  g
              式中:〈 e      ,g〉为向量 e        与 g的内积;g为第 S个隐含节点的输出;n为样本数;对于任意 g ∈Γ
                                             S
                       S - 1,j
                             S
                                                        S
                                       S - 1,j
                                                           +
                                                  g
              ( Γ表示张成的空间)都有,0 ≤‖g ‖<b,b∈R ,lim μ S                = 0 且 0< μ S ≤1 - r,r为正则化参数,范围在 0
                                                      g
                                                             S →!
                                                             和 b。
              到 1之间。根据监督机制确定隐含节点最佳参数 ω s
                                                                 s
                  (4)利用最小二乘法计算出隐含层输出权重:
                                                                        S
                                                                              2
                                                          ] = argmin Z - ∑  g                          (12)
                                          [ β 1 ,β 2 ,…,β S               β j j
                                                                       j =1
                  不断增加隐藏层节点,重复式( 8)到式(12),直至模型误差‖e‖达到期望的误差容限 χ 或者隐
                                                                            S
              函节点数达到最大隐含节点数 S ,最终输出最优模型。同时,参考文献 [29] 中 SCN超参数设置原
                                           max
              则,文中将输入权重和偏置的比例因子 λ设置为{0.5,1,5,10,30,50,100,150,200,250},正
              则化参数 r设置为{0.9,0.99,0.9999,0.99999,0.999999},误差容限 χ 设置为 0.01。
              3 仿真实验
              3.1 RCMMSDE抗噪性能分析 为验证 RCMMSDE算法的抗噪性能,通过式(13)模拟不同噪声环境下
              水电机组 “外 8”形的转子轴心轨迹。如图 2所示,加入噪声后,机组转子轴心轨迹进一步复杂化。
                                         {       1         ) + Asin(2 ω t + α 2 ) + z(t)               (13)
                                                               2
                                          x(t) =Asin( ω t + α 1
                                          y(t) =Bsin( ω t + β 1 ) + Bsin(2 ω t + β 2 ) + z(t)
                                                  1            2
              式中:x(t)、y(t)为 X和 Y方向的振动位移;ω为角频率;A、A、B 以及 B 为不同方向位移的幅
                                                                           2
                                                                      1
                                                                               1
                                                                                      2
                                       为 不 同 方 向 位 移 的 初 始 相 位;z(t)为 噪 声 信 号,文 中 取 50dB、40dB、
              度值;α 1  、α 2 、β 1  以 及 β 2
              30dB、20dB以及 10dB的高斯白噪声。













                                               图 2 不同状态不对中故障轴心轨迹图



                     3
                —  1 1 0 —
   112   113   114   115   116   117   118   119   120   121   122