一種綜合的供應(yīng)商評(píng)價(jià)方法

范 琛，王效俐

（同濟(jì)大學(xué) 經(jīng)濟(jì)與管理學(xué)院，上海 200092）

供應(yīng)商選擇問題（VSP）又叫供應(yīng)商評(píng)價(jià)問題，是供應(yīng)鏈管理研究中一個(gè)重要問題.從20世紀(jì)60年代Dickson等學(xué)者對(duì)供應(yīng)商的評(píng)價(jià)準(zhǔn)則進(jìn)行研究開始，到90年代，在全球化分工的背景下該問題得到了廣泛的研究和討論.

劉曉等［1］對(duì) VSP的模型方法進(jìn)行了綜述.Rao［2］詳細(xì)介紹了從20世紀(jì)90年代末到近幾年供應(yīng)商選擇方法的發(fā)展，并詳細(xì)介紹了TOPSIS等算法和算例.Degraeve等［3］從采購(gòu)成本的角度對(duì)幾種方法的案例應(yīng)用結(jié)果進(jìn)行了比較，指出了傳統(tǒng)評(píng)價(jià)模型中忽略的評(píng)價(jià)商務(wù)成本、訂貨成本等間接成本有可能導(dǎo)致在某些情況下評(píng)價(jià)結(jié)果失真.Wan［4］應(yīng)用了數(shù)學(xué)規(guī)劃法來解決VSP問題，利用了類似數(shù)據(jù)包絡(luò)分析（DEA）的變權(quán)思想簡(jiǎn)化并重新構(gòu)造了一個(gè)考慮偏好的線性規(guī)劃模型.Petroni等［5］將主成分分析（PCA）應(yīng)用到供應(yīng)商評(píng)價(jià)中，但未對(duì)評(píng)價(jià)權(quán)重設(shè)置作進(jìn)一步的討論.比較常見的方法還包括1種或多種方法的綜合應(yīng)用，如數(shù)據(jù)包絡(luò)分析同層次分析法（AHP）［6］和受險(xiǎn)價(jià)值法［7］、BP神經(jīng)網(wǎng)絡(luò)［8］等方法結(jié)合的綜合方法，以及其他一些方法如遺傳算法［9］、模糊規(guī)劃［10］、模糊評(píng)價(jià)［11］等.

然而以上各種方法都存著各種各樣的局限性.層次分析法在一定程度上不能精確地給出準(zhǔn)則權(quán)值，并且準(zhǔn)則是人為定義的，無法進(jìn)一步挖掘數(shù)據(jù)之間存在的內(nèi)部聯(lián)系.數(shù)據(jù)包絡(luò)分析無法有效地反應(yīng)決策者對(duì)判斷準(zhǔn)則的偏好.同時(shí)，由于未分析準(zhǔn)則之間的相關(guān)性，導(dǎo)致重復(fù)評(píng)價(jià).主成分分析法僅僅根據(jù)每個(gè)主成分的貢獻(xiàn)率來判斷權(quán)值，不僅無法反應(yīng)決策者的決策偏好，更無法科學(xué)地得到合理的權(quán)值.

本文試圖通過綜合AHP，PCA和DEA 3種方法來克服各方法本身的缺點(diǎn)，從而更好地獲得在VSP問題中的應(yīng)用效果.

1 結(jié)合AHP，PCA和DEA的綜合評(píng)價(jià)法

1.1 步驟1

（1）利用層次分析法［12］畫出層次圖，確定評(píng)估準(zhǔn)則.確定投入性準(zhǔn)則（m個(gè)）和產(chǎn)出性準(zhǔn)則（s個(gè)）.

（2）構(gòu)造判斷矩陣，通過式（1）和式（2）對(duì)判斷矩陣進(jìn)行一致性檢驗(yàn)，算出反映決策者偏好的權(quán)值.其中設(shè)投入準(zhǔn)則的權(quán)值為wxi（i＝1，2，…，m），產(chǎn)出準(zhǔn)則的權(quán)值為wyj（j＝1，2，…，s）.

式中：CI為AHP檢驗(yàn)中的中間檢驗(yàn)參數(shù)；λ為特征值.

式中：CR為AHP中一致性矩陣的檢驗(yàn)參數(shù)；RI的取值見參考文獻(xiàn)［12］.

1.2 步驟2

（1）進(jìn)行主成分分析，將樣本數(shù)據(jù)標(biāo)準(zhǔn)化.

（2）根據(jù)樣本數(shù)據(jù)分別計(jì)算每個(gè)樣本投入準(zhǔn)則和產(chǎn)出準(zhǔn)則矩陣的樣本協(xié)方差的無偏估計(jì)∑x和∑y（見文獻(xiàn)［13］），并分別求出各自的特征值λxi（i＝1，2，…，m）和λyj（j＝1，2，…，s）和對(duì)應(yīng)的特征向量pxi（i＝1，2，…，m）和pyj（j＝1，2，…，s）.對(duì)于λxi和λyj，滿足以下排列：λx1≥λx2≥…≥λxm，λy1≥λy2≥…≥λys.

（3）分別對(duì)投入矩陣和產(chǎn)出矩陣的協(xié)方差的特征值進(jìn)行排序，根據(jù)式（3）計(jì)算累計(jì)貢獻(xiàn)率Ψm，取Ψm達(dá)到一個(gè)較大值時(shí)（通常為85%）的前m′個(gè)和前s′個(gè)特征值＝（λ1，λ2，…，λm′）T（m′≤m）和＝（λ1，λ2，…，λs′）T（s′≤s）所對(duì)應(yīng)的特征向量＝（px1，px2，…，pxm′）和＝（py1，py2，…，pys′）構(gòu)成主成分的投影矩陣，對(duì)各主成分的意義進(jìn)行分析.

（4）分別對(duì)投入準(zhǔn)則的權(quán)重和產(chǎn)出準(zhǔn)則的權(quán)重進(jìn)行投影.

式中：wx，wy分別為投入準(zhǔn)則和產(chǎn)出準(zhǔn)則的權(quán)值向量.

在投入矩陣中希望總投入低，相應(yīng)地就希望權(quán)重高的投入準(zhǔn)則的評(píng)價(jià)值小.投影后，vi≥0表示某樣本第i個(gè)產(chǎn)出主成分評(píng)價(jià)值越大，反映在該主成分中的權(quán)重高的投入準(zhǔn)則下的評(píng)價(jià)值就越高，此時(shí)，希望該投入矩陣的主成分評(píng)價(jià)值越小越好；如果vi≤0，為了滿足投入矩陣評(píng)價(jià)值仍然為越小越好，則令pxi＝-pxi，可得

同理在產(chǎn)出矩陣中希望總產(chǎn)出高，相應(yīng)地就希望權(quán)重高的產(chǎn)出準(zhǔn)則評(píng)價(jià)值高.投影后，如果uj≥0，希望該主成分的評(píng)價(jià)值越高越好；如果uj≤0，為了滿足產(chǎn)出矩陣評(píng)價(jià)值仍然為越高越好，則令pyj＝-pyj可求得

同時(shí)根據(jù)vi和uj的值大小的排序分別獲得投入和產(chǎn)出準(zhǔn)則主成分的偏好約束，vi≥vj（i，j∈1，2，…，m′），ui≥uj（i，j∈1，2，…，s′）.對(duì)投入矩陣和產(chǎn)出矩陣作投影，得投影矩陣cx和cy.

由于DEA方法中隱含著投入x和產(chǎn)出y的線性關(guān)系的假設(shè)，根據(jù)以下性質(zhì)有cx和cy仍滿足線性關(guān)系.

性質(zhì)1 設(shè)x＝（x1，x2，…，xm）T和y＝（y1，y2，…，ym）T滿足線性關(guān)系，y＝f（x）＝Ax＋b，A為s×m階常數(shù)矩陣，b為s×1階常數(shù)向量.分別對(duì)x和y進(jìn)行線性變換，cx＝pxx和cy＝pyy，則對(duì)cx和cy仍滿足線性關(guān)系cy＝Bcx＋d，B為s′×m′階常數(shù)矩陣，d為s′×1階常數(shù)向量.

證明 因?yàn)閏x＝xpx，則有x＝cx，其中為廣義逆矩陣，關(guān)于廣義逆矩陣定義見文獻(xiàn)［13］，同理有y＝cy，因?yàn)閤和y滿足線性關(guān)系，y＝Ax＋b，則有cy＝Acx＋b，cy＝pyAcx＋pyb.令B＝pyA，d＝pyb，則有cy＝Bcx＋d.以上得證.

1.3 步驟3

將cx和cy分別作為投入矩陣和產(chǎn)出矩陣作為DEA的投入和產(chǎn)出項(xiàng)并考慮偏好約束，將基本的C2R模型規(guī)劃改寫為如下的分式規(guī)劃：

其中，前2項(xiàng)約束表示投入產(chǎn)出比都小于1，同時(shí)總收入和總投入都為正.第3項(xiàng)和第4項(xiàng)約束為偏好約束.進(jìn)行Charnes-Copper變換，令，μ＝tu，則可將原約束轉(zhuǎn)化為

由于cx和cy中各成分的評(píng)價(jià)值并不一定為正，對(duì)比原C2R模型規(guī)劃，規(guī)劃式（10）并沒有滿足cx＞0和cy＞0的先決條件.

1.4 投影后投入矩陣和產(chǎn)出矩陣經(jīng)濟(jì)解釋

由于cx是投入矩陣的投影，所以它仍然反映針對(duì)各供應(yīng)商的投入的評(píng)價(jià)值.但是新的相互獨(dú)立的各個(gè)主成分反映的并不一定是某一個(gè)或幾個(gè)準(zhǔn)則的評(píng)價(jià)值的投影之和，很可能反映的是幾個(gè)準(zhǔn)則的評(píng)價(jià)值同另外幾個(gè)準(zhǔn)則的評(píng)價(jià)值的投影差值.那么從經(jīng)濟(jì)角度考慮，相對(duì)希望其中投入權(quán)值大的準(zhǔn)則的評(píng)價(jià)值小，這樣綜合的總投入才有可能小.

2 實(shí)證分析

某大型工業(yè)分銷企業(yè)在華設(shè)立了采購(gòu)中心，根據(jù)其與其工作人員訪談和問卷調(diào)查得到的該中心有關(guān)供應(yīng)商判斷的準(zhǔn)則，畫出以下層次，如圖1.得到判斷矩陣見表1和表2.

圖1 判斷準(zhǔn)則的層次Fig.1 Hierarchy process chart

表1 準(zhǔn)則層第1層判斷矩陣Tab.1 Criteria judgment matrix of the first layer

表2 準(zhǔn)則層第2層判斷矩陣Tab.2 Criteria judgment matrix of the first layer

進(jìn)行一致性檢驗(yàn)，λ＝4.0261，CI＝0.0087，RI＝0.9，CR＝0.0097＜0.1，可見通過一致性檢驗(yàn).由于第2層的每個(gè)判斷矩陣都是2個(gè)準(zhǔn)則比較，所以無需檢驗(yàn).用幾何平均法求出第1層的權(quán)重，推導(dǎo)出準(zhǔn)則層第2層的權(quán)重，如圖2.

圖2 各準(zhǔn)則層權(quán)重Fig.2 Weight of each criteria

通過調(diào)查獲得一部分金屬加工類產(chǎn)品供應(yīng)商的評(píng)價(jià)數(shù)據(jù)，見表3.

表3 數(shù)據(jù)標(biāo)準(zhǔn)化后的12個(gè)供應(yīng)商評(píng)價(jià)值Tab.3 Normalized evaluation data of 12target vendors

分別對(duì)投入矩陣和產(chǎn)出矩陣進(jìn)行主成分分析的計(jì)算，貢獻(xiàn)率見表4、表5，同時(shí)計(jì)算權(quán)值投影，如表6、表7.獲得偏好約束v2≥v1.獲得產(chǎn)出主成分的偏好約束u1≥u2≥u3.

表4 投入矩陣協(xié)方差陣的特征值和累計(jì)貢獻(xiàn)率Tab.4 The accumulated contribution rate of the covariance matrix of i_n_p_u_t___m_a_t_r_ix

對(duì)原評(píng)價(jià)矩陣進(jìn)行投影，如表8.

根據(jù)式（11）進(jìn)行計(jì)算，得到表9.

表5 產(chǎn)出矩陣協(xié)方差的特征值和累計(jì)貢獻(xiàn)率Tab.5 The accumulated contribution rate of the covariance matrix of output matrix

綜合方法同帶偏好約束的DEA相比，結(jié)果見表10.

分析原評(píng)價(jià)矩陣中的投入準(zhǔn)則和產(chǎn)出準(zhǔn)則之間的相關(guān)性，求得相關(guān)性矩陣（方法見文獻(xiàn)［14］）見表11和表12.

表6 投入矩陣主成分分析和權(quán)重投影Tab.6 Principle components of input matrix and we ght projection

表7 產(chǎn)出矩陣的主成分分析和權(quán)重投影Tab.7 Principle components of output matrix and weight projection

表8 12家供應(yīng)商投影后各主成分的投影評(píng)價(jià)值Tab.8 Evaluation value of each principle component

表9 12家供應(yīng)商DEA得分和權(quán)重Tab.9 DEA score and weight of 12vendors

從表11可以看出，報(bào)價(jià)同響應(yīng)時(shí)間存在著很小的正相關(guān)性，報(bào)價(jià)和距離存在著很小的負(fù)相關(guān)性，而距離和響應(yīng)時(shí)間之間存在著一定的負(fù)相關(guān)性.

表10 綜合方法同帶偏好約束的DEA的計(jì)算結(jié)果的比較Tab.10 Result comparison between the synthesis method of AHP，PCA ＆ DEA and the method of DEA with preference constrain

表11 投入矩陣相關(guān)系數(shù)矩陣Tab.11 Correlation coefficient input matrix

表12 產(chǎn)出矩陣相關(guān)系數(shù)矩陣Tab.12 Correlation coefficient output matrix

從表12可以看出產(chǎn)能庫存、管理水平、技術(shù)能力之間都存在著較大的正相關(guān)性.對(duì)比表11，可以看出絕大部分的工廠得分在綜合方法下相比帶偏好的DEA得分要低，這部分是因?yàn)楫a(chǎn)出之間的相關(guān)性出現(xiàn)了重復(fù)評(píng)價(jià)，一定程度上夸大了工廠的實(shí)際產(chǎn)出評(píng)價(jià).

從投入角度來說，距離和響應(yīng)時(shí)間之間的較高負(fù)相關(guān)性導(dǎo)致原投入評(píng)價(jià)準(zhǔn)則下獲得的總評(píng)價(jià)值比獨(dú)立情況下高，也就是原準(zhǔn)則的總投入被少評(píng)估了一部分.

對(duì)比表4和表11，工廠1雖然在原產(chǎn)出準(zhǔn)則中各項(xiàng)得分都很高，由于原產(chǎn)出準(zhǔn)則的相互之間的高相關(guān)性產(chǎn)出評(píng)價(jià)被夸大了，同時(shí)相比工廠3，工廠1投入準(zhǔn)則的報(bào)價(jià)高出很多造成在相互獨(dú)立的準(zhǔn)則評(píng)價(jià)下得分比帶偏好約束的DEA方法下小得較多，得分排名也相差較大.而工廠1、工廠3、工廠4、工廠5在帶偏好約束的DEA方法下得分非常相近，而在本文綜合方法下得分被更加清晰地區(qū)分了出來.

3 結(jié)語

結(jié)合AHP，PCA和DEA 3種方法提出了一種綜合的供應(yīng)商評(píng)價(jià)法，利用AHP獲得評(píng)價(jià)準(zhǔn)則和決策者對(duì)準(zhǔn)則的偏好.利用PCA方法分別對(duì)評(píng)價(jià)樣本的投入準(zhǔn)則得分和產(chǎn)出準(zhǔn)則得分進(jìn)行正交投影，從而獲得相互獨(dú)立評(píng)價(jià)主成分得分，對(duì)準(zhǔn)則的主觀權(quán)值進(jìn)行投影，獲得關(guān)于主成分的偏好約束.采用DEA方法在傳統(tǒng)的C2R中加入了總投入和總產(chǎn)出為正的約束，同時(shí)增加了偏好約束以針對(duì)投影后的評(píng)價(jià)矩陣構(gòu)造合適的規(guī)劃.

通過實(shí)際案例將以上方法同帶偏好約束的DEA方法進(jìn)行對(duì)比，可以看出綜合方法能有效解決在供應(yīng)商評(píng)價(jià)過程中容易出現(xiàn)的數(shù)據(jù)冗余、重復(fù)評(píng)價(jià)夸大供應(yīng)商評(píng)價(jià)的問題，同時(shí)又集合了各個(gè)方法的優(yōu)點(diǎn)，從而能更有效地幫助決策者在眾多的供應(yīng)商中作出識(shí)別和判斷.

［1］劉曉，李海越，王成恩，等.供應(yīng)商選擇模型與方法綜述［J］.中國(guó)管理科學(xué)，2004（1）：139.LIU Xiao，LI Haiyue，WANG Cheng’en，et al.Vendor selection module and method review［J］.China Management Science，2004（1）：139.

［2］Rao R V.Decision making in the manufacturing environment［M］.London：Springer-Verlag London Ltd，2007.

［3］Degraeve Ze-Ger，Labro Eva，Roodhooft Filip.An evaluation of vendor selection models from a total cost of ownershipperspective［J］.European Journal of Operational Research，2000，125（1）：34.

［4］Wan LUNG Ng.An efficient and simple model for multiple criteria supplier selection problem［J］.European Journal of Operational Research，2008，186（3）：1059.

［5］Petroni Braglia Mar-Cello.Vendor Selection using principal component analysis ［J］. Journal of Supply Chain Management，2000，36（2）：63.

［6］Zhang H，Liu X，Liu W H.An AHP／DEA methodology for 3PL vendor selection in 4PL ［J］.Computer Supported Cooperative Work in Design，2006，LNCS 3865：646.

［7］WU D S，Olson D.Enterprise risk management：a DEA VaR approach in vendor selection［J］.International Journal of Production Research，2010，48（16）：4919.

［8］史成東，陳菊紅，張雅琪.基于BP神經(jīng)網(wǎng)絡(luò)和DEA的物流供應(yīng)商選擇［J］.工業(yè)工程，2010，13（4）：112.SHI Chengdong，CHEN Juhong，ZHANG Yaqi.Logistic supplier selection based on BP neural networks and DEA［J］.Industry Engineering，2010，13（4）：112.

［9］Bas-Net C，Weintraub A.A genetic algorithm for a bi-criteria supplier selection problem［J］.International Transactions in Operational Research，2009，16（2）：173.

［10］WANG J Y，ZHAO R Q，TANG W S.Fuzzy programming models for vendor selection problem in a supply chain［J］.Tsinghua Science ＆ Technology，2008，13（1）：106.

［11］趙小惠，趙小苗.基于模糊決策的供應(yīng)商選擇方法［J］.工業(yè)工程與管理，2002（4）：27.ZHAO Xiaohui，ZHAO Xiaomiao.Vendor selection based on fuzzy decision making ［J］.Industrial Engineering and Management，2002（4）：27.

［12］劉新憲，朱道立.選擇與判斷［M］.上海：上海科學(xué)普及出版社，1990.LIU Xinxian，ZHU Daoli.Selection and judgement［M］.Shanghai：Shanghai Popular Science Press，1990.

［13］王靜龍.多元統(tǒng)計(jì)分析［M］.北京：科學(xué)出版社，2008.WANG Jinglong. Multivariate statistics analysis ［M ］.Beijing：Science Press，2008.