姚頔，王瑛，馬躍飛

1.空軍工程大學(xué)裝備管理與安全工程學(xué)院，陜西西安 710051；2.國家飛行流量監(jiān)控中心，北京 100094 3.國防科學(xué)技術(shù)大學(xué)信息系統(tǒng)與管理學(xué)院，湖南長沙 410073

一種基于多子網(wǎng)與線性目標規(guī)劃的模糊網(wǎng)絡(luò)分析法模型

姚頔1，2，王瑛1，馬躍飛3

1.空軍工程大學(xué)裝備管理與安全工程學(xué)院，陜西西安 710051；2.國家飛行流量監(jiān)控中心，北京 100094 3.國防科學(xué)技術(shù)大學(xué)信息系統(tǒng)與管理學(xué)院，湖南長沙 410073

針對經(jīng)典ANP模型在解決現(xiàn)實復(fù)雜評價決策問題中的困難，通過引入“子網(wǎng)”結(jié)構(gòu)，提出了一個多層遞階分解的網(wǎng)絡(luò)模型，有效降低了建模工作量與復(fù)雜度。同時，考慮比較判斷的不確定性與多樣性，提供了集成多種表達的基于線性目標規(guī)劃的模糊判斷矩陣權(quán)向量求解算法。整合多子網(wǎng)與模糊判斷，構(gòu)建了評估流程。以空管體系為例進行能力評估，驗證了所提模型算法的有效性。

網(wǎng)絡(luò)分析法；模糊判斷；線性目標規(guī)劃；子網(wǎng)

Satty于1996年提出了網(wǎng)絡(luò)分析法（analytic network process，ANP）作為層次分析法（analytic hierarchy process，AHP）的延伸［1］，并獲得廣泛應(yīng)用［2-5］。從網(wǎng)絡(luò)結(jié)構(gòu)看，ANP將元素劃分為控制層和網(wǎng)絡(luò)層，控制層包括目標及控制準則，一般控制準則呈遞階層次、同層相互獨立；控制準則下所有支配元素分簇（或稱元素組）并按簇內(nèi)簇間相互關(guān)系抽象成網(wǎng)絡(luò)結(jié)構(gòu)。Satty和王蓮芬將ANP的主要結(jié)構(gòu)歸納為內(nèi)部獨立的遞階層次結(jié)構(gòu)、內(nèi)部獨立外部循環(huán)系統(tǒng)、內(nèi)部依存的遞階層次結(jié)構(gòu)和內(nèi)部依存外部循環(huán)系統(tǒng)［6-7］。2012年，Navid等提出了多影響樹（multiple-effect tree）和多準則樹（multiple-criteria tree）2種新的網(wǎng)絡(luò)結(jié)構(gòu)［8］。從測度原理看，延續(xù)AHP基于比較對判斷矩陣的影響排序思想，ANP以簇間影響為子塊，以控制層元素為準則，以網(wǎng)絡(luò)層某一元素為次準則，通過間接優(yōu)勢度比較，導(dǎo)出網(wǎng)絡(luò)層元素對次準則的影響排序，逐個考量所有元素相互影響并進行加權(quán)歸一化處理，構(gòu)建反映元素間一步優(yōu)勢度的超矩陣，利用冪法等求取極限相對排序向量。同AHP一樣，作為算法基礎(chǔ)的判斷矩陣假定決策者可以對任意2個元素進行比較判斷，并給出精確的比例標度。

1　多子網(wǎng)模型描述

1.1模型結(jié)構(gòu)

定義1 子網(wǎng) ANP網(wǎng)絡(luò)中由相同類別、屬性或按某種規(guī)則聚合的元素及其相互關(guān)系構(gòu)成的網(wǎng)絡(luò)，是多子網(wǎng)模型的基本組成單元，包括源節(jié)點、控制準則、底層節(jié)點3個層次元素及層內(nèi)層間關(guān)系，由如下五元組表示：

式中，Si是第i個子網(wǎng)；Oi是該子網(wǎng)的源節(jié)點；Ci是控制準則集；Bi是底層節(jié)點集；CRi是底層節(jié)點和控制準則間隸屬關(guān)系的集合，即CRi=｛〈u，v〉｜u∈Ci，v∈Bi｝；IRi是底層節(jié)點之間影響關(guān)系的集合，即IRi=｛〈x，y〉｜x，y∈Bi｝。設(shè)連接函數(shù)f滿足：若底層節(jié)點v在控制準則u下，則f〈u，v〉=1，反之為0；若底層節(jié)點y影響x，則f〈x，y〉=1（即x指向y），若無相互影響，則f〈x，y〉=0。

子網(wǎng)與經(jīng)典模型結(jié)構(gòu)相似，將元素劃分為控制層和受控層，控制層包括目標及控制準則，受控層為底層節(jié)點。子網(wǎng)的源節(jié)點即目標，源節(jié)點之下按控制準則將底層節(jié)點分為相互獨立的部分，每個控制準則支配若干底層節(jié)點，底層節(jié)點按相互影響關(guān)系分為依賴與獨立2類，依賴底層節(jié)點構(gòu)成網(wǎng)絡(luò)結(jié)構(gòu)，獨立底層節(jié)點構(gòu)成遞階層次結(jié)構(gòu)。可以沒有準則，但必須有目標和底層節(jié)點。

與經(jīng)典模型不同在于：①為控制子網(wǎng)規(guī)模，建立的子網(wǎng)控制準則僅為1層，其主要功能是將底層指標分類。②子網(wǎng)中沒有經(jīng)典模型中簇的概念，建立簇旨在降低判斷矩陣維度，通過將元素分配到各個子網(wǎng)中，得到的元素規(guī)模一般不會很大，因此簇的構(gòu)建是不必要的。下圖是一個典型的子網(wǎng)結(jié)構(gòu)：

圖1　子網(wǎng)的典型結(jié)構(gòu)

圖中，root是該子網(wǎng)的源節(jié)點；a、b是控制準則；c、d和e是具有相互影響關(guān)系的底層節(jié)點；f、g是相互獨立的底層節(jié)點。

定義2 子網(wǎng)間關(guān)系分為遞階層次關(guān)系和影響關(guān)系2類。遞階層次關(guān)系描述子網(wǎng)間的層次結(jié)構(gòu)，即?SP，Sq（p，q∈I，p≠q），若?z滿足：z∈Bp且z=Oq，則稱Sq是由SP派生的子網(wǎng)，SP和Sq呈遞階層次關(guān)系，子網(wǎng)間所有遞階層次關(guān)系的集合構(gòu)成H=｛〈Sp，Sq〉｜Sqis derived from Sp｝；影響關(guān)系描述源節(jié)點在同一個子網(wǎng)中的不同子網(wǎng)間的影響關(guān)系，即若Sm和Sn是由Sl派生的子網(wǎng)（m≠n），則顯然有：〈Sm，Sn〉∈IRl。

定義3 多子網(wǎng)模型（multi-subnets analytic network process，MSANP） 由子網(wǎng)及網(wǎng)間關(guān)系構(gòu)成的遞階層次結(jié)構(gòu)模型，由如下二元組表示：

式中，S=｛Si｝是子網(wǎng)集合，H是子網(wǎng)遞階層次關(guān)系集合。模型中除處于頂層的源節(jié)點外，每個節(jié)點僅能有一個源節(jié)點。

圖2給出了一個多子網(wǎng)模型示例，圖中共有3個子網(wǎng)Si、Sj和Sk。其中，Sj和Sk是由Si派生的子網(wǎng)，其源節(jié)點b、c分別為Si中源節(jié)點a下的底層節(jié)點，b 和c的依賴關(guān)系反映了Sj和Sk的相互影響；Sj中準則d下底層節(jié)點h、i、j具有依賴關(guān)系，并與e相互獨立，顯然b節(jié)點有自影響關(guān)系；Sk中底層節(jié)點f、g相互獨立。

圖2　一個MSANP模型示例

1.2模型復(fù)雜度分析

證明：

證畢。

相較ANP模型，MSANP模型有較少的元素間影響關(guān)系，相應(yīng)超矩陣的規(guī)模也變小了，復(fù)雜度被顯著降低。

2　基于線性目標規(guī)劃的模糊判斷

AHP/ANP模型的基礎(chǔ)信息源于構(gòu)建的判斷矩陣，一般矩陣單元表現(xiàn)形式一致，對于不同表達的判斷，傳統(tǒng)方法無法求得權(quán)向量，為此采用如下模型求解。

設(shè)存在歸一化模糊權(quán)向量

若該向量是歸一的，當且僅當如下條件成立［9］：

由于

如果求得的權(quán)向量完全一致，即滿足

式中

在實際問題中，由于判斷的不一致性，上式存在一定偏差，即

構(gòu)建如下線性目標規(guī)劃模型求解模糊權(quán)向量［10］：

式中，eT=［1，1，…，1］。

K（aij）＞1，即矩陣存在猶豫數(shù)時，將分解為個矩陣，f=1，2，…，ψ。對于任意，按（6）式求得歸一的模糊權(quán)向量及其偏移量。由此可以近似得到對于專家判斷的隸屬度為

t的值根據(jù)實際情況選取，一般可設(shè)t=1。為得到最符合專家判斷的模糊權(quán)向量，拓展式（6）得：

通過此模型求解一致性最大的模糊權(quán)向量。當判斷矩陣不一致性超出設(shè)定閾值η，即最小偏差量Jf?≥η時，需要專家對判斷矩陣重新調(diào)整。

3　MSFANP算法

以網(wǎng)絡(luò)類比連通圖，節(jié)點的影響分為子網(wǎng)影響、單分支全網(wǎng)影響和全網(wǎng)影響。子網(wǎng)影響（local value）是節(jié)點在所處子網(wǎng)中的局部影響值，包括遞階層次結(jié)構(gòu)下的聚合權(quán)重和網(wǎng)絡(luò)結(jié)構(gòu)下的極限相對權(quán)重；單分支全網(wǎng)影響（absolute value）是節(jié)點在從頂層源節(jié)點出發(fā)的一條鏈上，相對頂層源節(jié)點的絕對影響值；節(jié)點可能被包含在多個子網(wǎng)內(nèi)（例如待評估的方案），合成全網(wǎng)下節(jié)點在各鏈上的絕對影響，得到全網(wǎng)影響（global/synthetic value）。算法流程如圖3所示。

步驟1 根據(jù)多子網(wǎng)模型剖分思想，劃分網(wǎng)絡(luò)結(jié)構(gòu)，將各子網(wǎng)底層節(jié)點按控制準則分類，確定相互獨立和依賴關(guān)系（如無控制準則，以子網(wǎng)源節(jié)點為控制準則確定關(guān)系）；

步驟2 如果子網(wǎng)中有控制準則或獨立底層節(jié)點，對所有控制準則或?qū)ⅹ毩⒌讓庸?jié)點按所屬控制準則構(gòu)建判斷矩陣，由（2）～（11）式得相應(yīng)權(quán)向量；

步驟3 如果子網(wǎng)中有依賴底層節(jié)點，依次以某一依賴節(jié)點為次準則，進行依賴底層節(jié)點間間接優(yōu)勢度比較，由（2）～（11）式導(dǎo)出相對次準則的影響排序，最終形成加權(quán)超矩陣Mi；

步驟4 設(shè)子網(wǎng)源節(jié)點的影響值為1，控制準則影響值則為判斷矩陣導(dǎo)出的自身權(quán)重，則獨立底層節(jié)點的子網(wǎng)影響為：

式中，Vib為獨立底層節(jié)點的子網(wǎng)排序向量，Wib為獨立底層節(jié)點相對其控制準則的歸一化權(quán)向量，Vc為控制準則權(quán)重；

步驟5 利用冪法求取Mi的極限相對排序：

式中，t為一較大正數(shù)，Mi每一列均為依賴底層節(jié)點相對控制準則的極限相對排序，取第1列設(shè)為Wdb，取代Wib代入（10）式得Vdb，若無控制準則，Wdb= Vdb；

步驟6 由于每個節(jié)點僅有一個源節(jié)點，則按遞階層次得節(jié)點的單分支全網(wǎng)影響值：

式中，R（j）為從頂層源節(jié)點出發(fā)的第j條鏈，為 R（j）上第k個節(jié)點，為的單分支全網(wǎng)影響，為的子網(wǎng)影響；

步驟7 合成不同子網(wǎng)中相同節(jié)點的單分支影響得節(jié)點全網(wǎng)影響：

式中，el為網(wǎng)絡(luò)節(jié)點，l=｛1，2，…，L｝，為el的全網(wǎng)影響，M為頂層源節(jié)點出發(fā)的鏈數(shù)，指示函數(shù)εj滿足：若為el，則εj=1，否則εj=0。

圖3　MSFANP算法流程

4　計算實例

為驗證多子網(wǎng)模型的有效性，以國家空管體系為例說明所述算法的評估流程。根據(jù)SHEL模型剖分，空管體系能力自頂向下評估層次見圖4，分為3層指標體系，底層為方案層，分別為待評估比較國家。按MSFANP模型構(gòu)建子網(wǎng)見圖5，人機環(huán)管緊密耦合，形成相對體系能力的子網(wǎng)；其下支配要素分別張成子網(wǎng)，并相互影響；方案自成子網(wǎng)，與每個底層指標相連，方案間沒有明顯影響關(guān)系，方案對目標的反饋也不明確。

圖4　空管體系能力評估層次模型

構(gòu)建各子網(wǎng)下指標間的判斷矩陣，統(tǒng)計輸入群專家的模糊判斷信息，以圖5a）子網(wǎng)為例見表1：

圖5　空管體系能力評估子網(wǎng)構(gòu)成

表1　以H為次準則的判斷矩陣

表2　以S為次準則的判斷矩陣

表3　以E為次準則的判斷矩陣

表4　以M為次準則的判斷矩陣

求得表1判斷矩陣模糊權(quán)向量為

去模糊化得：（0.056，0.224，0.168，0.552）T

表2判斷矩陣的去模糊化權(quán)向量為

表3判斷矩陣的去模糊化權(quán)向量為

表4判斷矩陣的去模糊化權(quán)向量為

以求得的各判斷矩陣權(quán)向量為列向量，形成超矩陣，通過冪法求得該子網(wǎng)元素極限排序向量

作為各元素的子網(wǎng)影響值。

同理，求得圖5b）～圖5e）各子網(wǎng)中底層指標的子網(wǎng)影響值和方案在各底層能力下得分（即方案子網(wǎng)在各底層能力指標下的子網(wǎng)影響值）。依次計算方案節(jié)點的單分支全網(wǎng)影響和全網(wǎng)影響，最終以最優(yōu)方案國的全網(wǎng)影響值為1，美、歐、我國空管體系能力相對評估值為（1，0.925，0.765）T。計算結(jié)果與文獻［13］以主要宏觀績效領(lǐng)域?qū)嵶C數(shù)據(jù)，采用熵權(quán)-ANP得出的體系能力評估結(jié)論較為相近，從定性判斷定量化的解析角度予以了佐證。

5　結(jié) 論

本文提出的MSFANP模型，拓展了網(wǎng)絡(luò)分析法的描述能力。通過引入子網(wǎng)，基于遞階分解思想給出了一個多層網(wǎng)絡(luò)框架，由于僅需考慮子網(wǎng)內(nèi)的影響關(guān)系，模型被大為簡化。同時，為處理各種不確定性判斷形式提供了統(tǒng)一范式和基于線性目標規(guī)劃的模糊權(quán)向量求解算法。通過逐步計算網(wǎng)絡(luò)節(jié)點的子網(wǎng)影響、單分支全網(wǎng)影響，得出節(jié)點全網(wǎng)影響，從而構(gòu)建了MSFANP的算法流程。一個空管體系能力評估的實例證明了所提模型更符合人們的分析習(xí)慣，也驗證了模型方法的可行性與有效性。下一步將根據(jù)具體問題深化判斷標度研究，進一步提升評估精細化水平。

［1］ Saaty T L.Decision Making with Dependence and Feedback：The Analytic Network Process［M］.RWS Publications，Pittsburgh，1996

［2］ Leung L C，Lam K C，Cao D.Implementing the Balanced Score Card Using the Analytic Hierarchy Process&the Analytic Network Process［J］.Journal of the Operational Research Society，2005，57（6）：682-691

［3］ Chung S H，Lee A H，Pearn W L.Analytic Network Process（ANP）Approach for Product Mix Planning in Semiconductor Fabricator［J］.International Journal of Production Economics，2005，96（1）：15-36

［4］ Guneri A F，Gengiz M，Seker S.A Fuzzy ANP Approach to Shipyard Location Selection［J］.Expert Systems with Applications，2009，36（4）：7992-7999

［5］ Aznar J，F(xiàn)errís-O?ate J，Guijarro F.An ANP Framework for Property Pricing Combining Quantitative and Qualitative Attributes ［J］.Journal of the Operational Research Society，2010，61（5）：740-755

［6］ Saaty T L.Fundamentals of the Analytic Network Process-Dependence and Feedback in Decision-Making with a Single Network ［J］.Journal of Systems Science and Systems Engineering，2004，13（2）：129-157

［7］ 王蓮芬.網(wǎng)絡(luò)分析法（ANP）的理論與算法［J］.系統(tǒng)工程理論與實踐，2001，21（3）：44-50

Wang Lianfen.The Theory and Algorithm of Analytic Network Process［J］.Systems Engineering-Theroy&Practice，2001，21 （3）：44-50（in Chinese）

［8］ Navid K，Afshin S M，Jalil S，Seghir Z.An Algorithm for the Analytic Network Process（ANP）Structure Design［J］.Journal of Multi-Criteria Decision Analysis，2012，19（1/2）：33-55

［9］ Wang Y M，Elhag T M S.On the Normalization of Interval and Fuzzy Weights［J］.Fuzzy Sets and Systems，2006，157：2456-2471

［10］Wang Y M，Kwai-Sang C.A Linear Goal Programming Priority Method for Fuzzy Analytic Hierarchy Process and Its Applications in New Product Screening［J］.International Journal of Approximate Reasoning，2008，49：451-465

［11］Vahdani B，Mousavi S M，Tavakkoli-Moghaddam R.Group Decision Making Based on Novel Fuzzy Modified TOPSIS Method ［J］.Applied Mathematical Modelling，2011，37（9）：72-93

［12］Lin H F，Lee G G.A Study of Service Quality Evalution Model for Virtual Knowledge Communities［J］.E-Commerce Study，2006，4（2）：211-234

［13］國家飛行流量監(jiān)控中心.我國與航空發(fā)達國家空管體系能力比較研究［R］.北京：國家空管委辦公室，2015

China State Air Traffic Flow Management Center.The Comparative Study on the Capability of Air Traffic Management System-of-Systems between China and Developed Countries［R］.Beijing，CATC，2015（in Chinese）

A Model of Fuzzy Analytic Network Process Based on Multi-Subnets and Linear Goal Programming

Yao Di1，2，Wang Ying1，Ma Yuefei3

1.College of Equipment Management&Safety Engineering，Air Force Engineering University，Xi′an 710051，China 2.State Air Traffic Flow Management Center，Beijing 100094，China 3.College of Information System&Management，National University of Defense Technology，Changsha 410073，China

Aiming at the difficulty of using the classical ANP model to solve complex decision problems，we propose a network model in multiple levels based on hierarchical decomposition，through introducing‘subnet’structure，which reduces the workload and complexity effectively.At the same time，considering the uncertainty and diversity of the judgment，we provide a linear goal programming model for fuzzy pairwise comparison matrices with multiple expressions.By adopting the extended ANP model with subnets and fuzzy judgments，we construct the evaluation process.Finally，an example for the capability evaluation of air traffic management system-of-systems indicates that the method is effective.

analytic hierarchy process，decision making，electric network analysis，flow charting，linear programming，mathematical models，matrix algebra，membership functions，vectors；analytic network process （ANP），fuzzy judgment，linear goal programming，subnet

N945

A

1000-2758（2015）06-1041-08

2015-04-02

國家自然科學(xué)基金（71171199）與國家空管“十二五”科研專項課題（GKG201401003）資助

姚頔（1984—），空軍工程大學(xué)博士研究生，主要從事信息系統(tǒng)工程與智能決策、空域與飛行流量管理的研究。