莫紅枝

（玉林師范學院 教育技術(shù)中心，廣西 玉林 537000）

0 引言

神經(jīng)網(wǎng)絡算法是一種非線性計算模型，近年來成為模式識別中常用的工具之一。在多層神經(jīng)網(wǎng)絡中，系統(tǒng)的性能不僅受到隱含層數(shù)、隱含層神經(jīng)元數(shù)量的影響，而且還與激勵函數(shù)的選取和訓練算法直接相關。

在目前的研究中，采用最為廣泛的為S型激勵函數(shù)[1-2]，S型函數(shù)容易減慢網(wǎng)絡的收斂速度，甚至可能導致陷入局部最小值[3]。針對這一問題，近幾年采用自適應激勵函數(shù)對神經(jīng)元輸入的加權(quán)和進行計算已經(jīng)成為一種趨勢，并應用于股票預測[4]、文字識別[5]等方面。本文針對常見的S型函數(shù),改進了自適應性激勵函數(shù)神經(jīng)網(wǎng)絡系統(tǒng)框架，提出基于馬爾可夫鏈的學習算法，并將其應用到故障診斷領域，取得了比較好的實驗結(jié)果。

1 自適應性神經(jīng)網(wǎng)絡

1.1 多層神經(jīng)網(wǎng)絡

多層神經(jīng)網(wǎng)絡一般包括一個輸入層、一個輸出層、一個或者多個隱藏層。隱藏層每一層網(wǎng)絡中都包含多個神經(jīng)元，對每個神經(jīng)元的輸入都是由上一層輸出的加權(quán)和，例如對于第n組樣本數(shù)據(jù)，第k層上的第j個神經(jīng)元的輸入可以通過計算上一層的加權(quán)和來計算，計算公式為：

其中，wkji表示上一層網(wǎng)絡中的第i個神經(jīng)元對于第k層上的第j個神經(jīng)元之間的權(quán)重值。yki(n)表示第 k層網(wǎng)絡上第i個神經(jīng)元的輸出值。當k=1時，表示該層為輸入層，其輸出值就是對于整個神經(jīng)網(wǎng)絡的輸入值，即y1i(n)=xi(n)，神經(jīng)元的輸出由輸入經(jīng)過激勵函數(shù)計算得到，即：

其中，φkj表示第k層第j個神經(jīng)元的激勵函數(shù)。為了增加神經(jīng)網(wǎng)絡算法的非線性映射能力，隱含層的激勵函數(shù)可以選取非線性函數(shù)，常見的包括雙正切函數(shù)和S型函數(shù)。輸出層上神經(jīng)元的輸出就是整個神經(jīng)網(wǎng)絡系統(tǒng)的輸出。為了獲得最優(yōu)的權(quán)重值，需要對神經(jīng)網(wǎng)絡進行訓練。在給定訓練樣本后，對神經(jīng)網(wǎng)絡進行訓練過程的實質(zhì)是不斷調(diào)整權(quán)重值，使神經(jīng)網(wǎng)絡計算的輸出值與理論上的輸出值之間的誤差最小，即：

其中，din和 yin分別為第i個神經(jīng)元上真實情況下的輸出值和實際計算結(jié)果的輸出值，N為訓練樣本的個數(shù)。

1.2 激勵函數(shù)

本文以傳統(tǒng)的三層神經(jīng)網(wǎng)絡結(jié)構(gòu)，設計出基于改進的S函數(shù)的自適應性神經(jīng)網(wǎng)絡。其特點在于隱含層的激勵函數(shù)不再是固定的函數(shù)，而是包含了可變參數(shù)的激勵函數(shù)，這種神經(jīng)網(wǎng)絡系統(tǒng)框架如圖1所示。其中輸入層的神經(jīng)元數(shù)量由選取的特征的個數(shù)決定(1，2，…，N)，輸出層神經(jīng)元個數(shù)為 4 個（S1，S2，S3，S4），隱含層神經(jīng)元的輸入是輸入層各個神經(jīng)元輸出的加權(quán)和，并采用了自適應性的激勵函數(shù)對隱含層的輸入進行計算(I1，I2…Ik)。計算結(jié)果通過加權(quán)求和作為輸出層的輸入。在輸入層并未采用任何激勵函數(shù)，輸出層采用經(jīng)典的S型激勵函數(shù)，如：

式（4）由S型函數(shù)演化而來，是一種常用的自適應性激勵函數(shù)，式中的α和β為可變參數(shù)。雖然該函數(shù)已經(jīng)應用到神經(jīng)網(wǎng)絡算法中，卻很少有文獻將其應用在解決機械設備故障分類問題中。

圖1 多層神經(jīng)網(wǎng)絡算法系統(tǒng)架構(gòu)

2 基于馬爾可夫鏈的訓練算法

2.1 算法描述

樣本訓練即是在給定一定數(shù)量的樣本時，利用式(3)對所有的權(quán)重進行最優(yōu)化估計的過程[6-7]。當樣本數(shù)據(jù)中含有噪聲時，會造成程序魯棒性很差，給傳統(tǒng)的訓練方法帶來困難，本文假設式(3)中的誤差服從于高斯分布，然后根據(jù)后驗概率構(gòu)造出馬爾可夫鏈，完成對權(quán)重的訓練，可以有效避免噪聲對訓練結(jié)果造成的影響，具有收斂速度快的優(yōu)勢。假設樣本中含有噪聲，因此實際輸出與理想輸出之間的關系為：

式(5)的含義是對權(quán)重和自適應性參數(shù)進行估計，首先建立出的最大似然估計為：

其中 θ={w，α，β}為要估計的參數(shù)向量，xi為第 i個含噪樣本。根據(jù)Hammersley-Clifford理論，在給定樣本X時，利用條件分布 pi(θi|θ{j≠i}，X，E)可以從聯(lián)合分布 p(θ|X，E)中產(chǎn)生足夠的點，趨近使得誤差E最小的最小二乘估計值。因此本文通過條件分布不斷調(diào)整權(quán)重和可變參數(shù)的值，使其得到訓練：

假設樣本數(shù)據(jù)中的噪聲符合正態(tài)分布，則：

因此：

在得到θ的條件分布后，則各個權(quán)重及可變參數(shù)可以通過以下的算法進行更新。

算法一：

輸入：樣本 X，迭代次數(shù) I，初始值 θ(0)={w0，α0，β0}

從上面算法可以看到，通過不斷對各個參數(shù)進行更新，形成了馬爾可夫鏈，最終可以得到最小二乘估計。

2.2 參數(shù)分析

下面以式（4）為例給出條件分布的計算公式：

（1）對于權(quán)重 wk

求取其分布時只需要將其他變量看作固定值，則可以得到其分布：

算法二：

第二步：u～U(0，1)

第三步：wk～(wk)

圖2 正常軸承振動信號

（2）對于參數(shù) α

通過簡單的推導可以得出參數(shù)α的條件分布仍然服從于正態(tài)分布：

（3）對于參數(shù) β

式（12）中的概率也是很難處理的，為簡化程序，同樣采用拒絕性采樣算法對β進行更新。

3 實驗結(jié)果仿真

3.1 數(shù)據(jù)準備

為對神經(jīng)網(wǎng)絡性能進行驗證，利用本文設計的自適應性神經(jīng)網(wǎng)絡設計出了分類器，應用于軸承故障診斷當中。選取的樣本數(shù)據(jù)來自于美國凱斯西儲大學股東軸承數(shù)據(jù)中心。軸承型號為SKF公司的6205-2RS型的深溝球軸承?？紤]了4種軸承故障，分別為內(nèi)圈單點故障、外圈點蝕及滾動體點蝕和正常工作信號。4種信號的波形分別如圖2～圖5所示。

訓練樣本空間總共選取了1 136個個體，每個個體包含512個采樣點。通過小波分解提取出了20個小波系數(shù)作為分類器的輸入。

3.2 訓練結(jié)果

圖3 滾動體點蝕振動信號

圖4 內(nèi)圈單點故障振動信號

圖5 外圈點蝕故障振動信號

由于每段數(shù)據(jù)提取的特征個數(shù)為20個，因此將分類器的輸入層神經(jīng)元個數(shù)設置為了20個。通過實驗得到了最佳的隱含層神經(jīng)元數(shù)量。輸出層神經(jīng)元個數(shù)對應于4種故障，最終的神經(jīng)網(wǎng)絡架構(gòu)和參數(shù)設置如表1所示。

表1 測試的神經(jīng)網(wǎng)絡架構(gòu)及其參數(shù)設置

其中S-MPL代表了S型函數(shù)作為隱含層激勵函數(shù)的神經(jīng)網(wǎng)絡系統(tǒng)。F1-MPL代表以式（4）中的函數(shù)作為激勵函數(shù)的神經(jīng)網(wǎng)絡系統(tǒng)。對F1-MPL的訓練過程如圖6和圖7所示。圖6顯示的是利用本文算法的訓練過程，其中σ=0.5，初始值在 0～1之間隨機生成。圖 7展示了利用共軛梯度法作為訓練算法的收斂過程。共軛梯度法是介于最速下降法與牛頓法之間的一個方法，它僅需利用一階導數(shù)信息，但克服了最速下降法收斂慢的缺點，又避免了牛頓法需要存儲和計算Hesse矩陣并求逆的缺點，學習率選為1.2。

圖6 本文算法訓練過程

圖7 共軛梯度算法訓練過程

從圖6和圖7中的對比可以看出，利用本文提出的算法在第5次更新時就基本可以達到穩(wěn)定，具有穩(wěn)定性高、收斂速度快的特點。

3.3 分類結(jié)果

通過訓練后兩種神經(jīng)網(wǎng)絡對4類信號最終的分類結(jié)果如表2所示。

表2 神經(jīng)網(wǎng)絡算法的比較結(jié)果

其中樣本個數(shù)一欄分別表示了4種類型的信號的樣本個數(shù)，S-MPL、F1-MPL分別指的是通過 S-MPL網(wǎng)絡和F1-MPL網(wǎng)絡分類正確的4種信號的數(shù)目。因此可以計算出兩種神經(jīng)網(wǎng)絡算法的分類精度如表3所示。

表3 分類器的分類精度比較結(jié)果

從表1中看出，本文的訓練算法比傳統(tǒng)固定型的S型函數(shù)更加耗時，這是由于將可變參數(shù)引入到激勵函數(shù)中后，系統(tǒng)在訓練時往往需要更多的運算。然而訓練樣本的收斂精度也有所提高，這表明了本文算法訓練精度也更高，因此可以推測，自適應性的神經(jīng)網(wǎng)絡應用到其他問題當中時，比傳統(tǒng)的神經(jīng)網(wǎng)絡更加容易搜索到全局最優(yōu)值。

表2和表3證明了將本文的訓練算法應用在解決滾動軸承故障診斷問題方面的優(yōu)越性，取得了更高的分類精度。對于正常信號、滾動體點蝕振動信號、內(nèi)圈單點故障信號、外圈點蝕故障信號的分類精度分別可以達到99.69%、99.15%、98.90%、99.67%，平均分類精度可以達到99.38%。

4 結(jié)論

本文對傳統(tǒng)的S型激勵函數(shù)進行了改進，提出一種自適應性的神經(jīng)網(wǎng)絡分類器；基于馬爾可夫鏈對神經(jīng)網(wǎng)絡進行訓練，提高了網(wǎng)絡訓練速度；最后，將該分類器應用到滾動軸承故障診斷問題中。結(jié)果證明，使用該分類器可以比傳統(tǒng)的S型神經(jīng)網(wǎng)絡分類器獲得更高的分類精度。

[1]唐貴基，范德功，胡愛軍，等.基于小波包能量特征向量神經(jīng)網(wǎng)絡的旋轉(zhuǎn)機械故障診斷[J].汽輪機技術(shù)，2006(3)：215-217.

[2]張來斌，崔厚璽，王朝暉，等.基于信息熵神經(jīng)網(wǎng)絡的風力發(fā)電機故障診斷方法研究[J].機械強度，2009(1)：132-135.

[3]BURSE K，YADAV R N，SHRIVASTAVA S C.Channel equalization using neural networks：a review[J].IEEE Transactions on Systems，Man，and Cybernetics Part C-Applications and Reviews，2010，40(3)：352-357.

[4]BILDIRICI M，ALP E A，ERSIN O O.TAR-cointegration neural network model：An empirical analysis of exchange rates and stock returns[J].Expert Systems with Applications，2010，37(1)：2-11.

[5]KANG M，PALMER-BROWN D.A modal learning adaptive function neural network applied to handwritten digit recognition[J].Information Sciences，2008，178(20)：3802-3812.

[6]滕輝.一種改進的神經(jīng)網(wǎng)絡學習算法研究[J].科技通報，2012(4)：97-98.

[7]鄭緒枝，夏薇，雷靖.一種改進的Jacobi正交多項式的BP神經(jīng)網(wǎng)絡算法[J].云南大學學報(自然科學版)，2011(S2)：188-191.