陶立權(quán)，楊坤

（中國民航大學(xué)天津市民用航空器適航與維修重點實驗室，天津　300300）

?

基于AHP-BP神經(jīng)網(wǎng)絡(luò)的航空材料供應(yīng)商選擇

陶立權(quán)，楊坤

（中國民航大學(xué)天津市民用航空器適航與維修重點實驗室，天津300300）

摘要：根據(jù)適航規(guī)章要求，提出了航空材料供應(yīng)商選擇的特色指標，構(gòu)建了航空材料供應(yīng)商選擇的AHP-BP神經(jīng)網(wǎng)絡(luò)模型，通過AHP確定各指標權(quán)重，再結(jié)合BP神經(jīng)網(wǎng)絡(luò)，從訓(xùn)練數(shù)據(jù)中提取隱含的知識和規(guī)律，能夠方便地用于新供應(yīng)商的選擇。該模型求解算法為動量梯度下降反向傳播算法，具有良好的可擴展性，從而增加了評價的動態(tài)性。算例驗證結(jié)果表明，將AHP-BP神經(jīng)網(wǎng)絡(luò)模型用于航空材料供應(yīng)商選擇具有較強的實用性。

關(guān)鍵詞：供應(yīng)商；適航；AHP；BP神經(jīng)網(wǎng)絡(luò)

適航性是航空器在預(yù)期的環(huán)境和運行限制下安全飛行的固有品質(zhì)。一架飛機大約由600萬個零件組成，每個零件都是由金屬或非金屬材料經(jīng)過一系列工藝程序制成的，所以材料是影響飛機安全、質(zhì)量和壽命的重要因素。

然而，目前國內(nèi)材料大多是由航空部門以外的材料供應(yīng)商提供的，這些供應(yīng)商大多缺乏甚至沒有適航經(jīng)驗，因此對材料供應(yīng)商的適航管理顯得尤為重要。根據(jù)適航規(guī)章要求，飛機制造商對供應(yīng)商提供零部件、設(shè)備及材料的適航性負責。供應(yīng)商的選擇對飛機制造商來說至關(guān)重要。

目前，關(guān)于航空供應(yīng)商管理的國外資料較少，國內(nèi)航空供應(yīng)商管理研究多集中于航材、備件供應(yīng)商的研究［1-6］。文獻［7］針對航空發(fā)動機企業(yè)采購成本高、供應(yīng)商快速反應(yīng)能力低等問題，構(gòu)造了航空發(fā)動機企業(yè)供應(yīng)商評價指標體系，依據(jù)供應(yīng)商評定系統(tǒng)優(yōu)選供應(yīng)商，最終得以有效控制成本。文獻［8］利用層次分析法研究了客機項目中的供應(yīng)商選擇問題，根據(jù)項目評價的一般原則，結(jié)合大飛機項目的特點，建立了大飛機項目的評價指標體系。文獻［9］在層次分析法（AHP）的基礎(chǔ)上，提出了民機供應(yīng)商質(zhì)量保證評價指標體系，并構(gòu)建了民機供應(yīng)商質(zhì)量保證體系的評價模型。上述研究從傳統(tǒng)的供應(yīng)商管理角度，對供應(yīng)商的選擇方法進行改進，尚未結(jié)合飛機型號取證的適航要求對供應(yīng)商的選擇進行完善。

本文根據(jù)適航規(guī)章要求，提出航空材料供應(yīng)商選擇的特色指標，構(gòu)建AHP-BP神經(jīng)網(wǎng)絡(luò)模型，實現(xiàn)航空材料供應(yīng)商的合理選擇。

1　基于AHP-BP神經(jīng)網(wǎng)絡(luò)的航材供應(yīng)商模型

AHP-BP神經(jīng)網(wǎng)絡(luò)方法將AHP與BP神經(jīng)網(wǎng)絡(luò)相結(jié)合，即利用AHP的計算結(jié)果來訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)，既避免了繁雜計算又簡化了網(wǎng)絡(luò)訓(xùn)練過程。根據(jù)以往專家或業(yè)內(nèi)知名人士對供應(yīng)商的評價，應(yīng)用神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)并積累其評價的經(jīng)驗，通過權(quán)值和閾值的學(xué)習(xí)，以及對網(wǎng)絡(luò)的訓(xùn)練，得到訓(xùn)練后的網(wǎng)絡(luò)模型，用來對供應(yīng)商進行綜合性評價，并預(yù)測供應(yīng)商的綜合評價得分，通過綜合得分的排序，選擇出最優(yōu)供應(yīng)商。

1.1航空材料供應(yīng)商選擇指標體系

中國民用航空局（CAAC）對供應(yīng)商的管理要求主要集中在：CCAR 21部21.13和21.143條款，適航程序AP-21-AA-2011-03、AP-21-04，以及咨詢通告AC-21-04、AC-ARJ21-01之中。針對大型運輸類飛機的適航要求，在CCAR 25部中予以規(guī)定。CCAR25.603 和CCAR25.605兩個條款對航空材料提出了適航要求。

本文在客觀分析影響材料供應(yīng)商選擇因素的基礎(chǔ)上，根據(jù)適航規(guī)章的管理要求和技術(shù)要求，通過專家問卷調(diào)查，確定從研發(fā)能力、產(chǎn)品質(zhì)量、交付能力、生產(chǎn)能力和客戶支持5個因素來評價材料供應(yīng)商，并進一步提出產(chǎn)品適航性、資質(zhì)認證等指標，如圖1所示。

圖1　飛機材料供應(yīng)商選擇指標體系Fig.1　Index system for aeronautical material supplier selection

1.2指標權(quán)重的確定

1）構(gòu)造兩兩比較判斷矩陣，如表1所示。

表1　判斷矩陣Tab.1　Judgment matrix

判斷矩陣中的數(shù)值要以Saaty專門設(shè)計的1～9的比例標度為比較標準，如表2所示。其中bij表示對于A而言，Bi對Bj的相對重要性。

表2　判斷尺度定義表Tab.2　Description of judgment scales

2）計算相對權(quán)重有很多方法［10］，如方根法、和積法、特征根法和最小二乘法。其中方根法最為簡便，在研究中應(yīng)用也最為廣泛，表達式為

對于每個判斷矩陣都需對其進行一致性檢驗。一致性是指判斷矩陣中各要素的重要性判斷是否一致，不能出現(xiàn)邏輯矛盾。首先要計算一致性指標

式中：λmax為判斷矩陣的最大特征值，可根據(jù)公式AW = λW計算出所有的特征值，再找出最大值。下一步是計算隨機一致性比率C.R.，其中R.I.的值由表3給出，即

表3　判斷矩陣的平均隨機性指標R.I.值Tab.3　Average random index R.I. of judgment matrix

當C.R.＜0.1時，則認為判斷矩陣具有滿意的一致性，否則就需調(diào)整判斷矩陣，使其具有滿意的一致性。

3）權(quán)重總排序及總一致性檢驗。權(quán)重總排序如表4所示。其中

向量（W1，W2，W3）為方案層C對目標層A的權(quán)重向量。權(quán)重大的方案優(yōu)先選擇。一致性指標計算公式為

表4　權(quán)重總排序Tab.4　Total rank of weight

1.3供應(yīng)商選擇流程

1）建立判斷矩陣

根據(jù)表2對影響因素B1、B2、B3、B4、B5進行兩兩比較，建立判斷矩陣，如表5所示。影響因素B1、B2、B3、B4、B5的評價指標間的判斷矩陣分別如表6～表10所示。

表5　準則層與目標層之間的判斷矩陣Tab.5 Judgment matrix between criterion layer and target layer

表6　指標B1的判斷矩陣Tab.6 Judgment matrix of index B1

表7　指標B2的判斷矩陣Tab.7 Judgment matrix of index B2

表8　指標B3的判斷矩陣Tab.8 Judgment matrix of index B3

表9　指標B4的判斷矩陣Tab.9 Judgment matrix of index B4

2）權(quán)重計算及一致性檢驗

首先進行權(quán)重單排序及一致性檢驗。指標權(quán)重根據(jù)式（1）和式（2）計算，一致性檢驗指標依據(jù)式（3）和式（4）以及表3進行計算。由于計算復(fù)雜，本文采用Matlab軟件來完成權(quán)重計算及一致性檢驗。對于判斷矩陣A，計算得出權(quán)重向量為W =（0.51 0.263 8 0.129 6 0.063 6 0.033 0），一致性指標C.R. = 0.040 9＜0.1。

表10　指標B5的判斷矩陣Tab.10 Judgment matrix of index B5

同理：對于判斷矩陣B1，權(quán)重向量為W1=（0.148 8 0.785 4 0.065 8），一致性指標C.R.1= 0.027 7＜0.1；對于判斷矩陣B2，權(quán)重向量為W2=（0.75 0.25），二階矩陣自然符合一致性，不需要一致性檢驗；對于判斷矩陣B3，權(quán)重向量為W3=（0.730 6 0.081 0 0.188 4），一致性指標C.R.3= 0.022 4＜0.1；對于判斷矩陣B4，權(quán)重向量為W4=（0.833 3 0.166 7），二階矩陣自然符合一致性，不需要一致性檢驗；對于判斷矩陣B5，權(quán)重向量為W5=（0.649 1 0.279 0 0.071 9），一致性指標C.R.5= 0.022 4＜0.1。

下面進行權(quán)重總排序及一致性檢驗指標層對目標層的權(quán)重（根據(jù)表4）計算，計算結(jié)果如表11所示。

表11　指標層對目標層的綜合權(quán)重Tab.11　Global weight of index layer relative to target layer

因此，指標層各指標對目標層的權(quán)重向量為

一致性檢驗指標根據(jù)式（5）來計算，結(jié)果如下

3）數(shù)據(jù)處理及BP網(wǎng)絡(luò)設(shè)計和訓(xùn)練

首先，需要由行業(yè)專家根據(jù)經(jīng)驗以及候選供應(yīng)商的實際情況，對候選供應(yīng)商的各個指標在1～10之間進行打分，各指標分數(shù)與其權(quán)重乘積之和為總分。

其次，要將計算得到的分數(shù)由（0，10）轉(zhuǎn)化為（-1，1），以便于BP網(wǎng)絡(luò)的訓(xùn)練。計算公式為

最后，采用BP神經(jīng)網(wǎng)絡(luò)通過對樣本的學(xué)習(xí)和訓(xùn)練，獲取相關(guān)領(lǐng)域?qū)＜业闹R和評價經(jīng)驗，當對新的供應(yīng)商做出綜合評價時，利用訓(xùn)練成熟網(wǎng)絡(luò)中存儲的專家經(jīng)驗對其做出綜合評價。

AHP-BP神經(jīng)網(wǎng)絡(luò)模型的建立過程如圖2所示。

圖2　AHP-BP神經(jīng)網(wǎng)絡(luò)模型建立過程Fig.2　Building process of AHP-BP neural network model

2　算例分析

案例中假設(shè)有25家候選材料供應(yīng)商，基于上述評價指標采用十分制對樣本的各項指標進行打分，如表12所示。

由于各項指標的量綱不同，需要做出相應(yīng)的轉(zhuǎn)換。MATLAB R2012a的神經(jīng)網(wǎng)絡(luò)工具箱中所定義的premnmx函數(shù)可以將數(shù)據(jù)歸一化到（-1，1）之間。將數(shù)據(jù)轉(zhuǎn)換后結(jié)合AHP得出的指標權(quán)重進行綜合評價，結(jié)果如表13所示。

本文利用MATLAB R2012a神經(jīng)網(wǎng)絡(luò)工具箱，以13個指標評分作為網(wǎng)絡(luò)的輸入，綜合評價值作為輸出，以規(guī)范化后的數(shù)據(jù)作為訓(xùn)練樣本，網(wǎng)絡(luò)擬合供應(yīng)商的綜合評價結(jié)果及網(wǎng)絡(luò)擬合誤差如圖3～圖5所示。

表12　供應(yīng)商評價指標評分Tab.12　Grade of supplier evaluation index

表13　綜合評分結(jié)果Tab.13　Global grade

圖3　訓(xùn)練樣本、驗證樣本、測試樣本及樣本總體R值Fig.3　R value of training，validation，test and all samples

圖4　網(wǎng)絡(luò)訓(xùn)練中MSE的變化過程Fig.4　MSE change in training

經(jīng)訓(xùn)練后，訓(xùn)練樣本、驗證樣本、測試樣本及樣本總體的R值都在0.95以上，圖3顯示了網(wǎng)絡(luò)性能指標MSE值達到1e-15以下，圖4顯示了網(wǎng)絡(luò)訓(xùn)練過程中3種樣本MSE值的變化，圖5顯示了目標輸出與實際輸出間的誤差絕大多數(shù)落在（-0.129 5，0.184 2），只有兩例不在該范圍，一個為-0.317 7，另一個為0.278 4。這幾個指標說明網(wǎng)絡(luò)已達到網(wǎng)絡(luò)訓(xùn)練的精度要求。

在此基礎(chǔ)上，利用訓(xùn)練好的網(wǎng)絡(luò)對5個新的供應(yīng)商進行綜合評價，其指標評分結(jié)果如表14所示。AHP模型指標權(quán)重綜合評價結(jié)果及BP神經(jīng)網(wǎng)絡(luò)評價結(jié)果對比如表15所示。

圖5　誤差直方圖Fig.5　Error histogram

表14　供應(yīng)商指標評分表Tab.14　Supplier index grades

表15　網(wǎng)絡(luò)擬合結(jié)果Tab.15　Neural network fitting results

表15中目標值是根據(jù)上文得到的各項指標權(quán)重與5組供應(yīng)商的評分無量綱后加權(quán)所得。從中可以看出，第4個供應(yīng)商是最佳供應(yīng)商，其綜合評價值為8.542。5個供應(yīng)商的網(wǎng)絡(luò)擬合值與目標值非常接近，最大絕對誤差只有0.005，因此可以認為，學(xué)習(xí)訓(xùn)練后的BP神經(jīng)網(wǎng)絡(luò)模型對供應(yīng)商的評價具有較好的學(xué)習(xí)效果和非線性映射能力。

上述分析表明，該模型求解算法為動量梯度下降反向傳播算法，具有很好的可擴展性，從而增加了評價的動態(tài)性。

3　結(jié)語

航空材料供應(yīng)商的選擇與一般制造企業(yè)供應(yīng)商的選擇有所不同。本文結(jié)合適航要求，針對航空材料供應(yīng)商建立了評價體系，并建立了AHP-BP神經(jīng)網(wǎng)絡(luò)模型來進行供應(yīng)商的選擇。AHP-BP神經(jīng)網(wǎng)絡(luò)方法將層次分析法與BP神經(jīng)網(wǎng)絡(luò)相結(jié)合，即利用層次分析法的計算結(jié)果來訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)，既避免了繁雜的計算又簡化了網(wǎng)絡(luò)的訓(xùn)練過程。通過算例驗證，AHP-BP神經(jīng)網(wǎng)絡(luò)模型能夠記錄相關(guān)領(lǐng)域?qū)＜业闹R和評價經(jīng)驗，有效減少人工誤差，提高了評價的準確性。對航空制造業(yè)科學(xué)地選擇材料供應(yīng)商具有實際的指導(dǎo)意義和推廣價值。

參考文獻:

［1］王琨.航材供應(yīng)商的評估與選擇［J］.江蘇航空，2002（4）：17-19.

［2］倪現(xiàn)存，左洪福.基于模糊多屬性群決策及OWA算子的民機備件供應(yīng)商評估［J］.機械科學(xué)與技術(shù)，2006，25（12）：1404-1407.

［3］倪現(xiàn)存，左洪福.航材承修商可拓綜合評價研究［J］.哈爾濱工業(yè)大學(xué)學(xué)報，2006，38（7）：1168-1172.

［4］陳玖圣，張曉瑜，陳靜杰.航材供應(yīng)商評估決策過程模型及實現(xiàn)［J］.系統(tǒng)工程理論與實踐，2008.11（28）：149-154.

［5］周媛，左洪福.民機供應(yīng)商綜合評估實證研究［J］.科技管理研究，2010（18）：70-72.

［6］黃勇，張磊.基于集對層次分析法的備件供應(yīng)商評估［J］.江蘇航空，2011，121（1）：21-23.

［7］傅艷華，唐加福.航空發(fā)動機企業(yè)供應(yīng)商評定系統(tǒng)構(gòu)建［J］.儀器儀表學(xué)報，2006，27（6）：938-939.

［8］胡慶江，郭靜雯.層次分析法在民用客機項目供應(yīng)商選擇中的應(yīng)用［J］.國防技術(shù)基礎(chǔ)，2009（10）：54-58.

［9］郭巧榮，徐建新.民機供應(yīng)商質(zhì)量保證體系的評價模型研究［J］.航空維修與工程，2012（3）：93-96.

［10］秦吉，張翼鵬.現(xiàn)代統(tǒng)計信息分析技術(shù)在安全工程方面的應(yīng)用-層次分析法原理［J］.工業(yè)安全與防塵，1999（5）：44-48.

（責任編輯：楊媛媛）

Selection of aeronautical material suppliers based on AHP-BP neural networks

TAO Liquan，YANG Kun

（Civil Aircraft Airworthiness and Maintenance Key Lab of Tianjin，CAUC，Tianjin 300300，China）

Abstract:According to the requirements of airworthiness regulations，a special index system for the selection of aeronautical material suppliers is proposed with the establishment of AHP-BP neural network model. The weight of each index is determined through AHP，and then combined with BP neural networks，the implicit knowledge and discipline from the training data is extracted，so that it can be used for new suppliers conveniently. The momentum gradient descent back-propagation algorithm is used in the model with well extensibility and increasing dynamic nature of the evaluation. Verification shows that AHP-BP neural network model is with strong practicability for the selection of aeronautical material suppliers.

Key words:supplier；airworthiness；AHP；BP neural networks

中圖分類號：V25；F253.2

文獻標志碼：A

文章編號：1674－5590（2016）02－0056－06

收稿日期：2014-11-07；修回日期：2015-03-25基金項目：民用飛機專項科研項目（MJ-J-2012-07）

作者簡介：陶立權(quán)（1978—），男，黑龍江牡丹江人，助理研究員，碩士，研究方向為民用飛機及發(fā)動機適航審定.