楊家豪，歐陽(yáng)森，石怡理，黃瑞藝，劉子文

(華南理工大學(xué)電力學(xué)院，廣東 廣州510640)

1 引言

模糊綜合評(píng)價(jià)方法是國(guó)內(nèi)外的電能質(zhì)量評(píng)價(jià)的主要手段之一［1-8］。其運(yùn)用于電能質(zhì)量的評(píng)價(jià)，關(guān)鍵在于隸屬度函數(shù)的選擇和設(shè)計(jì)。在實(shí)際評(píng)價(jià)過(guò)程中，選擇不同的隸屬度函數(shù)會(huì)導(dǎo)致評(píng)價(jià)結(jié)果出現(xiàn)差異，因此設(shè)計(jì)一種具有普遍適用性的隸屬度函數(shù)來(lái)提高評(píng)價(jià)結(jié)果的可信度具有重要意義。

文獻(xiàn)［2-5］中隸屬度函數(shù)的參數(shù)的確定需要依靠實(shí)際情況和個(gè)人的經(jīng)驗(yàn)進(jìn)行反復(fù)設(shè)計(jì)，使評(píng)價(jià)結(jié)果具有人為主觀性而導(dǎo)致可信度降低;文獻(xiàn)［6，10］選擇了單一的嶺形或三角形隸屬度函數(shù)，根據(jù)最大隸屬度原則和極限原則確定左右零點(diǎn)，沒有考慮隸屬度函數(shù)的選擇對(duì)評(píng)價(jià)結(jié)果的影響;文獻(xiàn)［7］中采用數(shù)學(xué)擬合并修正后得到隸屬度函數(shù)中參數(shù)，這種方法的可信度受到曲線擬合精度的限制。文獻(xiàn)［9］中采用Bezier曲線逼近方式，當(dāng)曲線的冪次較低時(shí)效果較差，此時(shí)靈活性受到限制，且控制點(diǎn)數(shù)目會(huì)對(duì)擬合出的隸屬函數(shù)精度造成影響。

優(yōu)選組合方法在預(yù)測(cè)領(lǐng)域得到了廣泛的應(yīng)用，在許多文獻(xiàn)中已經(jīng)驗(yàn)證“組合預(yù)測(cè)模型的預(yù)測(cè)誤差平方和不大于參與組合的各個(gè)單一模型的預(yù)測(cè)誤差平方和的最小者”［11，12］。因此組合的結(jié)果可以有效減小波動(dòng)，提高穩(wěn)定性，同時(shí)其適用性必然是最好的。

針對(duì)不同隸屬度函數(shù)會(huì)導(dǎo)致評(píng)價(jià)結(jié)果出現(xiàn)差異的問(wèn)題，本文擬選取五種典型隸屬度函數(shù)進(jìn)行模糊評(píng)價(jià)，再引入方差—協(xié)方差優(yōu)選組合法［12］，對(duì)五種典型的隸屬度函數(shù)進(jìn)行賦權(quán)組合，構(gòu)造出一種組合隸屬度函數(shù)。方差—協(xié)方差優(yōu)選組合法對(duì)所選取的隸屬度函數(shù)的評(píng)價(jià)信息實(shí)現(xiàn)了充分利用，對(duì)多種單一模糊評(píng)價(jià)模型所包含的信息進(jìn)行最優(yōu)組合，從而使組合隸屬度函數(shù)具有普遍適用性，最終使評(píng)價(jià)結(jié)果充分可信。文章最后給出一個(gè)算例，通過(guò)實(shí)例驗(yàn)證了本算法的有效性。

2 典型隸屬度函數(shù)分析

設(shè)在論域U上給定一個(gè)映射:

A:U→［0，1］

u a A(u)

則稱A為U上的模糊集;A(u)稱為A的隸屬度函數(shù)(或稱為u對(duì)A的隸屬度)。

隸屬度函數(shù)又可以分為:偏小型、中間型、偏大型。偏小型和偏大型只是單側(cè)區(qū)間有值，而中間型的特點(diǎn)是區(qū)間對(duì)稱有值。由于電能質(zhì)量評(píng)價(jià)涉及到的指標(biāo)較多，為了使選取的隸屬度函數(shù)均能適用于各個(gè)指標(biāo)的評(píng)價(jià)，提高普遍適用性，故本文均采用中間型隸屬度函數(shù)，保證其對(duì)稱性和平衡性［13］。

本文選取五種典型的隸屬度函數(shù)作為構(gòu)造組合隸屬度函數(shù)的基礎(chǔ)，如表1所示。由表1可知，每一種隸屬度函數(shù)都具有自身的特點(diǎn)，當(dāng)選用不同的隸屬度函數(shù)進(jìn)行模糊評(píng)價(jià)時(shí)，得到的隸屬度必然不同，從而導(dǎo)致評(píng)價(jià)結(jié)果出現(xiàn)差異。同時(shí)，每個(gè)隸屬度函數(shù)中又含有一個(gè)或兩個(gè)待定的參數(shù)(表達(dá)式中a、b均為待定參數(shù))，同樣的隸屬度函數(shù)，只要參數(shù)的設(shè)計(jì)不同，亦會(huì)導(dǎo)致隸屬度出現(xiàn)差異。因此，在單一隸屬度函數(shù)下得到的評(píng)價(jià)結(jié)果可信度不高。

表1 五種典型的隸屬度函數(shù)Tab．1 Fivetypicalmembershipfunctions

3 基于優(yōu)選組合法的隸屬度函數(shù)構(gòu)造

3.1 同一評(píng)價(jià)指標(biāo)的隸屬度平均值

先求取同一指標(biāo)在不同隸屬度函數(shù)下的隸屬度，再求取這些隸屬度的平均值:

式中，fi表示該指標(biāo)在第i個(gè)隸屬度函數(shù)下的隸屬度值;m表示總共有m種隸屬度函數(shù)，本文中m=5。

3.2 方差—協(xié)方差優(yōu)選組合法

設(shè)f1，f2是在兩個(gè)不同隸屬度函數(shù)下關(guān)于f－的無(wú)偏隸屬度值，fc是加權(quán)平均的組合隸屬度值，偏離f－的誤差分別為e1，e2和ec。w1和w2是相應(yīng)的權(quán)系數(shù)，且w1+w2=1，則有:fc=w1f1+w2f2，若fc也為無(wú)偏的，則可以證明:ec=e1+e2，從而有fc的方差為:

關(guān)于w1對(duì)Var(ec)求極小值，可得:

又w2=1－w1，記方差Var(e2)=σ22，Var(e1)=σ11，協(xié)方差cov(e1，e2)=σ12。由于e1、e2相互獨(dú)立，則σ12=0，有:

把兩種隸屬度函數(shù)的組合結(jié)果一般化。設(shè)m個(gè)無(wú)偏隸屬度值分別為f1，f2，．．．，fm，各自相對(duì)于f－的誤差的方差為σ11，σ22，．．．，σkk。在許多實(shí)際情況下，不同隸屬度函數(shù)之間的誤差是不相關(guān)的。關(guān)于wi(i=1，2，．．．，m)對(duì)Var(ec)求極小值是個(gè)求條件極值的問(wèn)題，在約束條件∑wi=1的條件下，對(duì)Var(ec)引入Lagrange乘子來(lái)求極小值，可以得到各個(gè)隸屬度函數(shù)的權(quán)重:

式中，i=1，2，…，m。

由此可見，隸屬度值越靠近平均值f－的隸屬度函數(shù)，即方差越小其權(quán)重就越大。

3.3 組合隸屬度函數(shù)

根據(jù)上面求得的各個(gè)隸屬度函數(shù)的權(quán)重，構(gòu)造組合隸屬度函數(shù)，如下:

從式(6)可知，新構(gòu)造的組合隸屬度函數(shù)f對(duì)五種隸屬度函數(shù)進(jìn)行優(yōu)化整合，突出評(píng)價(jià)結(jié)論一致性較高的部分隸屬度函數(shù)的作用，同時(shí)兼顧一致性程度較低的另一部分函數(shù)所承載的評(píng)價(jià)信息，即僅弱化其作用，而不完全剔除該信息，從而既對(duì)評(píng)價(jià)結(jié)論相對(duì)一致的多評(píng)價(jià)結(jié)果進(jìn)行了優(yōu)化組合，又充分利用了所有的評(píng)價(jià)方法信息，最終使應(yīng)用組合隸屬度函數(shù)進(jìn)行評(píng)價(jià)所得的評(píng)價(jià)結(jié)果的可信度更高。

4 電能質(zhì)量模糊評(píng)價(jià)算法設(shè)計(jì)

4.1 電能質(zhì)量指標(biāo)權(quán)重的確定

本文側(cè)重于使用組合隸屬度函數(shù)實(shí)現(xiàn)對(duì)評(píng)價(jià)結(jié)果的優(yōu)化，故在權(quán)重的確定上并沒有強(qiáng)制性的要求，為不失一般性，擬采用G1法(序關(guān)系分析法)［14］來(lái)確定電能質(zhì)量指標(biāo)的權(quán)重，該方法過(guò)程清晰簡(jiǎn)便，能充分體現(xiàn)專家的意愿。主要指標(biāo)包括電壓偏差、頻率偏差、諧波電壓總畸變率、三相電壓不平衡度以及電壓閃變5項(xiàng)。步驟如下:

(1)首先確定序關(guān)系。設(shè)對(duì)于某評(píng)價(jià)準(zhǔn)則指標(biāo)的重要程度具有的序關(guān)系x1?x2?···?xm。

(2)然后確定相鄰指標(biāo)之間的相對(duì)重要程度。設(shè)專家關(guān)于評(píng)價(jià)指標(biāo)xk－1與xk的重要程度之比wk－1/wk的理性判斷分別為:

其中，對(duì)xk－1與xk的重要程度之比的判斷分為五種情況，分別對(duì)應(yīng)不同的rk，其賦值可參考表2。

(3)最后計(jì)算權(quán)重系數(shù)wk:

由wk得出其他指標(biāo)的權(quán)重:

表2 rk賦值參考表Tab．2 Assignment reference of rk

4.2 電能質(zhì)量評(píng)價(jià)

在選定的隸屬度函數(shù)下，可得到對(duì)應(yīng)的模糊隸屬矩陣:

式中，n為電能質(zhì)量指標(biāo)數(shù);m為質(zhì)量等級(jí)數(shù);Ri為第i項(xiàng)指標(biāo)的單因素評(píng)價(jià)。

矩陣?yán)锏摩蘨，j須進(jìn)行歸一化處理，即為:

最后的評(píng)價(jià)結(jié)果為:其中，W為電能質(zhì)量指標(biāo)權(quán)重矢量，W=(w1，w2，…，wm);B為相對(duì)于各質(zhì)量等級(jí)隸屬度向量。

對(duì)5個(gè)質(zhì)量等級(jí)依次賦以分值c1=5、c2=4、c3=3、c4=2、c5=1，分別代表“優(yōu)質(zhì)”、“良好”、“中”、“合格”、“不合格”，最后采用加權(quán)平均得到評(píng)價(jià)點(diǎn)的電能質(zhì)量評(píng)價(jià)值:

4.3 算法步驟

利用方差—協(xié)方差法對(duì)選定的五種隸屬度函數(shù)進(jìn)行優(yōu)化組合，從而構(gòu)造出組合隸屬度函數(shù)并進(jìn)行模糊評(píng)價(jià)。算法包括五個(gè)步驟:

(1)進(jìn)行電能質(zhì)量等級(jí)劃分。

(2)利用G1法來(lái)確定電能質(zhì)量指標(biāo)的權(quán)重。

(3)分別利用五種典型的隸屬度函數(shù)計(jì)算電能質(zhì)量各指標(biāo)數(shù)據(jù)的隸屬度值，利用式(1)計(jì)算各指標(biāo)的隸屬度值的平均值。

(4)利用式(5)計(jì)算各隸屬度函數(shù)對(duì)于各電能質(zhì)量指標(biāo)的權(quán)重。

(5)利用式(6)構(gòu)造出相對(duì)于各個(gè)指標(biāo)的組合隸屬度函數(shù);同時(shí)利用組合隸屬度函數(shù)計(jì)算各電能質(zhì)量指標(biāo)數(shù)據(jù)的隸屬度值。

(6)利用式(10)、式(11)分別計(jì)算在五種隸屬度函數(shù)和組合隸屬度函數(shù)下的模糊評(píng)價(jià)矩陣R。

(7)利用式(7)～式(9)計(jì)算電能質(zhì)量指標(biāo)權(quán)重W。再利用式(12)得到隸屬度向量;再通過(guò)式(13)進(jìn)行加權(quán)平均得到評(píng)價(jià)值。

5 實(shí)例分析

5.1 計(jì)算過(guò)程

應(yīng)用五種典型隸屬度函數(shù)和組合隸屬度函數(shù)對(duì)某變電站電能質(zhì)量數(shù)據(jù)進(jìn)行模糊評(píng)價(jià)，并進(jìn)行對(duì)比分析。本文劃分的電能質(zhì)量等級(jí)如表3所示。

表3 電能質(zhì)量指標(biāo)及其質(zhì)量等級(jí)劃分準(zhǔn)則Tab．3 Power quality indexes and their division norms of quality level

本文數(shù)據(jù)來(lái)源于廣州供電局某110kV變電站1～12月的電能質(zhì)量數(shù)據(jù)，如表4所示。根據(jù)表3中的電能質(zhì)量數(shù)據(jù)，利用5種典型的隸屬度函數(shù)計(jì)算每一個(gè)數(shù)據(jù)關(guān)于各電能質(zhì)量等級(jí)的隸屬度值，再利用式(1)和式(5)計(jì)算各個(gè)隸屬度函數(shù)在每一個(gè)電能質(zhì)量指標(biāo)下的權(quán)重，如表5所示。

表4 變電站電能質(zhì)量數(shù)據(jù)Tab．4 Power quality data of substation

表5 各電能指標(biāo)下5種隸屬度函數(shù)的權(quán)重Tab．5 Weights of fivemembership functions of each power indicator

從而根據(jù)式(6)計(jì)算各電能質(zhì)量指標(biāo)的組合隸屬度函數(shù)，再利用組合隸屬度函數(shù)計(jì)算每一個(gè)數(shù)據(jù)關(guān)于各質(zhì)量等級(jí)的隸屬度值。

電能質(zhì)量指標(biāo)按照“諧波電壓總畸變率?電壓偏差?電壓閃變?三相電壓不平衡度?頻率偏差”的序關(guān)系，利用G1法得出各電能質(zhì)量指標(biāo)權(quán)重為W=(0.2728，0.248，0.206667，0.15897，0.1135);分別把5種隸屬度函數(shù)的隸屬度值和組合函數(shù)算得的隸屬度值代入式(10)中得到每月關(guān)于不同隸屬度函數(shù)的模糊評(píng)判矩陣，再利用式(11)～式(13)得到每個(gè)月的關(guān)于不同隸屬度函數(shù)的模糊評(píng)價(jià)結(jié)果，如表6和圖1所示。

其中，Y1在嶺型分布隸屬度函數(shù)下12個(gè)月的評(píng)價(jià)值，Y2在尖Γ分布隸屬度函數(shù)下12個(gè)月的評(píng)價(jià)值，Y3在柯西分布隸屬度函數(shù)下12個(gè)月的評(píng)價(jià)值，Y4在梯形分布隸屬度函數(shù)下12個(gè)月的評(píng)價(jià)值，Y5在正態(tài)分布隸屬度函數(shù)下12個(gè)月的評(píng)價(jià)值，Y6在組合隸屬度函數(shù)下12個(gè)月的評(píng)價(jià)值。

5.2 結(jié)果分析

從表6和圖1可知，五種典型的隸屬度函數(shù)算出來(lái)的12個(gè)月模糊評(píng)價(jià)值分布在2.8～3.5之間;而評(píng)價(jià)值等于3時(shí)，是電能質(zhì)量等級(jí)“中”和“合格”的分界線，這時(shí)就會(huì)出現(xiàn)在某個(gè)隸屬度函數(shù)下該月的電能質(zhì)量等級(jí)為“中”，而在另一隸屬度函數(shù)下為合格的現(xiàn)象，因而導(dǎo)致評(píng)價(jià)結(jié)果具有不確定性。此時(shí)，基于優(yōu)選組合的隸屬度函數(shù)，對(duì)各隸屬度函數(shù)進(jìn)行了最優(yōu)組合，必然比單一隸屬度函數(shù)的評(píng)價(jià)結(jié)果更具有可信性。

表6 對(duì)應(yīng)各個(gè)隸屬度函數(shù)的12個(gè)月電能質(zhì)量綜合評(píng)價(jià)值YTab．6 Power quality comprehensive evaluation value Y of12 months corresponding to each membership function

(1)從表6可以看出，4月份的評(píng)價(jià)值中，隸屬度函數(shù)為柯西分布和正態(tài)分別時(shí)分別比3小0.0205和0.0959，而隸屬度函數(shù)為嶺型分布、尖Γ分布和梯形分布分別比3大0.1078、0.163和0.124，綜合來(lái)看，4月份的評(píng)價(jià)值更趨向于大于3，即電能質(zhì)量等級(jí)為中;利用組合隸屬度函數(shù)算出的評(píng)價(jià)值3.041399，驗(yàn)證了上述分析，因此有理由相信4月份的電能質(zhì)量的質(zhì)量等級(jí)為中更可信。同樣，對(duì)于7月份，有三個(gè)評(píng)價(jià)值小于3，兩個(gè)評(píng)價(jià)值大于3，且小于3的評(píng)價(jià)值偏離“合格”的程度高于大于3的評(píng)價(jià)值，因此，類似4月份的分析過(guò)程，7月份評(píng)價(jià)值應(yīng)小于3，即電能質(zhì)量等級(jí)為“合格”;利用組合隸屬度函數(shù)算出的評(píng)價(jià)值2.974526，也驗(yàn)證了上述分析。對(duì)于其他月份，進(jìn)行類似的分析，可以得到類似的分析結(jié)果。

圖1 對(duì)應(yīng)各個(gè)隸屬度函數(shù)的12個(gè)月電能質(zhì)量綜合評(píng)價(jià)值YFig．1 Power quality comprehensive evaluation value Y of 12 months corresponding to eachmembership function

(2)由圖1可以看出，優(yōu)選組合分布的評(píng)價(jià)結(jié)果分布在五種隸屬度函數(shù)評(píng)價(jià)結(jié)果中間，類似于評(píng)價(jià)結(jié)果取了平均值，使評(píng)價(jià)結(jié)果的可信度提高。必須指出，這與單純計(jì)算評(píng)價(jià)值的平均值存在本質(zhì)區(qū)別，因?yàn)榉讲睢獏f(xié)方差優(yōu)選組合法是從隸屬度層面進(jìn)行賦權(quán)的，當(dāng)某一隸屬度函數(shù)算出的隸屬度值偏離平均隸屬度值很大時(shí)，認(rèn)為該隸屬度函數(shù)不夠精確，因而賦予的權(quán)重相應(yīng)就小，反之亦然。通過(guò)對(duì)隸屬度函數(shù)加權(quán)得到了組合隸屬度函數(shù)更準(zhǔn)確，從而提高了評(píng)價(jià)結(jié)果的可信度。

6 結(jié)論

在電能質(zhì)量模糊評(píng)價(jià)過(guò)程中，隸屬度函數(shù)的選擇問(wèn)題由于人為的主觀性及經(jīng)驗(yàn)性太強(qiáng)，要得到一個(gè)讓所有人都認(rèn)可的結(jié)果近乎不可能。將基于多種隸屬度函數(shù)的評(píng)價(jià)結(jié)果通過(guò)賦予不同的權(quán)重來(lái)綜合分析是非常有益的嘗試。

本文引入方差—協(xié)方差最優(yōu)組合法，對(duì)五種典型隸屬度函數(shù)進(jìn)行最優(yōu)組合，構(gòu)造了具有普遍適用性的組合隸屬度函數(shù)，其中突出了評(píng)價(jià)結(jié)論一致性較高的那部分隸屬度函數(shù)的作用，同時(shí)兼顧一致性程度較低的另一部分函數(shù)，僅起將其弱化的作用，而不完全剔除該方法，充分利用了五種典型函數(shù)的評(píng)價(jià)信息，使模糊評(píng)價(jià)結(jié)果的可信度得到提高。

本文所提出的算法很好地處理了單一隸屬度函數(shù)導(dǎo)致評(píng)價(jià)結(jié)果出現(xiàn)的差異性。通過(guò)實(shí)例分析驗(yàn)證了本文算法的有效性。

［1］Farghal SA，Kandil M S，Elmitwally A．Quantifying electric power quality via fuzzymodeling and analytic hierarchy processing［J］．IEEE Proceedings Generation Transmission and Distribution，2002，149(1):44-49．

［2］賈清泉，宋家驊，蘭華，等(Jia Qingquan，Song Jiahua，Lan Hua，etal．)．電能質(zhì)量及其模糊方法評(píng)價(jià)(Power quality and fuzzymethod evaluation)［J］．電網(wǎng)技術(shù)(Power System Technology)，2000，24(6):46-49．

［3］趙 霞，趙 成 勇，賈 秀 芳，等(Zhao Xia，Zhao Chengyong，Jia Xiufang，et al．)．基于可變權(quán)重的電能質(zhì)量模糊綜合評(píng)價(jià)(Fuzzy synthetic evaluation of power quality based on changeable weight)［J］．電網(wǎng)技術(shù)(Power System Technology)，2005，29(6):11-16．

［4］李連結(jié)，姚建剛，龍立波，等(Li Lianjie，Yao Jiangang，Long Libo，et al．)．組合賦權(quán)法在電能質(zhì)量模糊綜合評(píng)價(jià)中的應(yīng)用(Application of combination weighting method in fuzzy synthetic evaluation of power quality)［J］．電力系統(tǒng)自動(dòng)化(Automation of Electric Power Systems)，2007，31(4):56-60．

［5］蔣金良，袁金晶，歐陽(yáng)森(Jiang Jinliang，Yuan Jinjing，Ouyang Sen)．基于改進(jìn)隸屬度函數(shù)的電能質(zhì)量模糊綜合評(píng)價(jià)方法研究(Research of fuzzy synthetic evaluation method of power quality based on improved membership function)［J］．華南理工大學(xué)學(xué)報(bào)(Journal of South China University of Technology)，2012，40(11):107-112．

［6］譚家茂，黃少先(Tan Jiamao，Huang Shaoxian)．基于模糊理論的電能質(zhì)量綜合評(píng)價(jià)方法研究(Study on power quality comprehensive evaluation based on fuzzy theory)［J］．繼電器(Relay)，2006，34(3):55-59．

［7］唐會(huì)智，彭建春(Tang Huizhi，Peng Jianchun)．基于模糊理論的電能質(zhì)量綜合量化指標(biāo)研究(Research on synthetic and quantificated appraisal index of power quality based on fuzzy theory)［J］．電網(wǎng)技術(shù)(Power System Technology)，2003，27(12):85-88．

［8］陶順，肖湘寧(Tao Shun，Xiao Xiangning)．基于短板效應(yīng)的電能質(zhì)量綜合等級(jí)評(píng)價(jià)(Synthetic grading evaluation of power quality based on incentivemechanism)［J］．電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2008，27(2):16-20．

［9］Medaglia A L，F(xiàn)ang SC，Nuttle H LW，et al．An efficient and flexiblemechanism for constructing membership function［J］．European J of Operational Research，2002，139(1):84-95．

［10］汪金剛，張作鵬，龍震澤，等(Wang Jingang，Zhang Zuopeng，Long Zhenze，etal．)．采用模糊推理對(duì)絕緣子放電紫外檢測(cè)結(jié)果進(jìn)行外絕緣放電狀況評(píng)估(Discharge status evaluation on external insulation of insulator based on results of ultra-violet radiation measurement and fuzzy reasoning)［J］．電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2011，30(2):79-83．

［11］汪同三，張濤(Wang Tongsan，Zhang Tao)．組合預(yù)測(cè):理論、方法及其應(yīng)用(Forecast combination:theory，method and application)［M］．北京:社會(huì)科學(xué)文獻(xiàn)出版社(Beijing:Social Sciences Academic Press)，2008．

［12］康重慶，夏清，劉梅(Kang Chongqing，Xia Qing，Liu Mei)．電力系統(tǒng)負(fù)荷預(yù)測(cè)(Power system load forecasting)［M］．北京:中國(guó)電力出版社(Beijing:China Electric Power Press)，2007．

［13］王季方，盧正鼎(Wang Jifang，Lu Zhengding)．模糊控制中隸屬度函數(shù)的確定方法(The determine method ofmembership function in fuzzy control)［J］．河南科學(xué)(Henan Science)，2000，18(4):348-351．

［14］郭亞軍(Guo Yajun)．綜合評(píng)價(jià)理論、方法及應(yīng)用(Comprehensive evaluation theory，methods and application)［M］．北京:科學(xué)出版社(Beijing:Science Press)，2007．