馮曉蒲,張鐵峰
(華北電力大學 電氣與電子工程學院,河北 保定 071003)
聚類分析是一種重要的人類行為,早在孩提時代,一個人就通過不斷改進下意識中的聚類模式來學會如何區分貓狗、動物植物。目前在許多領域都得到了廣泛的研究和成功的應用,如用于模式識別、數據分析、圖像處理、市場研究、客戶分割、Web文檔分類等[1]。
聚類就是按照某個特定標準(如距離準則)把一個數據集分割成不同的類或簇,使得同一個簇內的數據對象的相似性盡可能大,同時不在同一個簇中的數據對象的差異性也盡可能地大。即聚類后同一類的數據盡可能聚集到一起,不同數據盡量分離。
聚類技術[2]正在蓬勃發展,對此有貢獻的研究領域包括數據挖掘、統計學、機器學習、空間數據庫技術、生物學以及市場營銷等。各種聚類方法也被不斷提出和改進,而不同的方法適合于不同類型的數據,因此對各種聚類方法、聚類效果的比較成為值得研究的課題。
目前,有大量的聚類算法[3]。而對于具體應用,聚類算法的選擇取決于數據的類型、聚類的目的。如果聚類分析被用作描述或探查的工具,可以對同樣的數據嘗試多種算法,以發現數據可能揭示的結果。
主要的聚類算法可以劃分為如下幾類:劃分方法、層次方法、基于密度的方法、基于網格的方法以及基于模型的方法[4-6]。
每一類中都存在著得到廣泛應用的算法,例如:劃分方法中的k-means[7]聚類算法、層次方法中的凝聚型層次聚類算法[8]、基于模型方法中的神經網絡[9]聚類算法等。
目前,聚類問題的研究不僅僅局限于上述的硬聚類,即每一個數據只能被歸為一類,模糊聚類[10]也是聚類分析中研究較為廣泛的一個分支。模糊聚類通過隸屬函數來確定每個數據隸屬于各個簇的程度,而不是將一個數據對象硬性地歸類到某一簇中。目前已有很多關于模糊聚類的算法被提出,如著名的FCM算法等。
本文主要對k-means聚類算法、凝聚型層次聚類算法、神經網絡聚類算法之SOM,以及模糊聚類的FCM算法通過通用測試數據集進行聚類效果的比較和分析。
k-means是劃分方法中較經典的聚類算法之一。由于該算法的效率高,所以在對大規模數據進行聚類時被廣泛應用。目前,許多算法均圍繞著該算法進行擴展和改進。
k-means算法以k為參數,把n個對象分成k個簇,使簇內具有較高的相似度,而簇間的相似度較低。kmeans算法的處理過程如下:首先,隨機地選擇k個對象,每個對象初始地代表了一個簇的平均值或中心;對剩余的每個對象,根據其與各簇中心的距離,將它賦給最近的簇;然后重新計算每個簇的平均值。這個過程不斷重復,直到準則函數收斂。通常,采用平方誤差準則,其定義如下:

這里E是數據庫中所有對象的平方誤差的總和,p是空間中的點,mi是簇Ci的平均值[9]。該目標函數使生成的簇盡可能緊湊獨立,使用的距離度量是歐幾里得距離,當然也可以用其他距離度量。k-means聚類算法的算法流程如下:
輸入:包含n個對象的數據庫和簇的數目k;
輸出:k個簇,使平方誤差準則最小。
步驟:
(1)任意選擇k個對象作為初始的簇中心;
(2)repeat;
(3)根據簇中對象的平均值,將每個對象(重新)賦予最類似的簇;
(4)更新簇的平均值,即計算每個簇中對象的平均值;
(5)until不再發生變化。
根據層次分解的順序是自底向上的還是自上向下的,層次聚類算法分為凝聚的層次聚類算法和分裂的層次聚類算法。
凝聚型層次聚類的策略是先將每個對象作為一個簇,然后合并這些原子簇為越來越大的簇,直到所有對象都在一個簇中,或者某個終結條件被滿足。絕大多數層次聚類屬于凝聚型層次聚類,它們只是在簇間相似度的定義上有所不同。四種廣泛采用的簇間距離度量方法如下:
最小距離:

最大距離:

平均值的距離:

平均距離:

這里給出采用最小距離的凝聚層次聚類算法流程:
(1)將每個對象看作一類,計算兩兩之間的最小距離;
(2)將距離最小的兩個類合并成一個新類;
(3)重新計算新類與所有類之間的距離;
(4)重復(2)、(3),直到所有類最后合并成一類。
SOM神經網絡[11]是由芬蘭神經網絡專家Kohonen教授提出的,該算法假設在輸入對象中存在一些拓撲結構或順序,可以實現從輸入空間(n維)到輸出平面(2維)的降維映射,其映射具有拓撲特征保持性質,與實際的大腦處理有很強的理論聯系。
SOM網絡包含輸入層和輸出層。輸入層對應一個高維的輸入向量,輸出層由一系列組織在2維網格上的有序節點構成,輸入節點與輸出節點通過權重向量連接。學習過程中,找到與之距離最短的輸出層單元,即獲勝單元,對其更新。同時,將鄰近區域的權值更新,使輸出節點保持輸入向量的拓撲特征。
算法流程:
(1)網絡初始化,對輸出層每個節點權重賦初值;
(2)將輸入樣本中隨機選取輸入向量,找到與輸入向量距離最小的權重向量;
(3)定義獲勝單元,在獲勝單元的鄰近區域調整權重使其向輸入向量靠攏;
(4)提供新樣本、進行訓練;
(5)收縮鄰域半徑、減小學習率、重復,直到小于允許值,輸出聚類結果。
1965年美國加州大學柏克萊分校的扎德教授第一次提出了‘集合’的概念。經過十多年的發展,模糊集合理論漸漸被應用到各個實際應用方面。為克服非此即彼的分類缺點,出現了以模糊集合論為數學基礎的聚類分析。用模糊數學的方法進行聚類分析,就是模糊聚類分析[12]。
FCM算法是一種以隸屬度來確定每個數據點屬于某個聚類程度的算法。該聚類算法是傳統硬聚類算法的一種改進。
設數據集X={x1,x2,…xn},它的模糊c劃分可用模糊矩陣U=[uij]表示,矩陣U的元素uij表示第 j(j=1,2,…,n)個數據點屬于第 i(i=1,2,…,c)類的隸屬度,uij滿足如下條件:

目前被廣泛使用的聚類準則為取類內加權誤差平方和的極小值,即:

其中V為聚類中心,m為加權指數,

算法流程:
(1)標準化數據矩陣;
(2)建立模糊相似矩陣,初始化隸屬矩陣;
(3)算法開始迭代,直到目標函數收斂到極小值;
(4)根據迭代結果,由最后的隸屬矩陣確定數據所屬的類,顯示最后的聚類結果。
實驗中,選取專門用于測試分類、聚類算法的國際通用的UCI數據庫中的IRIS[13]數據集,IRIS數據集包含150個樣本數據,分別取自三種不同的鶯尾屬植物setosa、versicolor和virginica的花朵樣本,每個數據含有4個屬性,即萼片長度、萼片寬度、花瓣長度,單位為cm。在數據集上執行不同的聚類算法,可以得到不同精度的聚類結果。
文中基于前面所述各算法原理及算法流程,用matlab進行編程運算,得到表1所示聚類結果。

表1 三種聚類方法的實驗對比結果
如表1所示,對于四種聚類算法,按三方面進行比較:(1)聚錯樣本數:總的聚錯的樣本數,即各類中聚錯的樣本數的和;(2)運行時間:即聚類整個過程所耗費的時間,單位為 s;(3)平均準確度:設原數據集有 k個類,用 ci表示第 i類,ni為 ci中樣本的個數,mi為聚類正確的個數,則mi/ni為第i類中的精度,則平均精度為:

四種聚類算法中,在運行時間及準確度方面綜合考慮,k-means和FCM相對優于其他。但是,各個算法還是存在固定缺點:k-means聚類算法的初始點選擇不穩定,是隨機選取的,這就引起聚類結果的不穩定,本實驗中雖是經過多次實驗取的平均值,但是具體初始點的選擇方法還需進一步研究;層次聚類雖然不需要確定分類數,但是一旦一個分裂或者合并被執行,就不能修正,聚類質量受限制;FCM對初始聚類中心敏感,需要人為確定聚類數,容易陷入局部最優解;SOM與實際大腦處理有很強的理論聯系。但是處理時間較長,需要進一步研究使其適應大型數據庫。
聚類分析因其在許多領域的成功應用而展現出誘人的應用前景,除經典聚類算法外,各種新的聚類方法正被不斷被提出。
[1]HAN Jia Wei,KAMBER M.數據挖掘概念與技術[M].范明,孟曉峰,譯.北京:機械工業出版社,2001.
[2]楊小兵.聚類分析中若干關鍵技術的研究[D].杭州:浙江大學,2005.
[3]XU Rui,Donald Wunsch 1 1.survey of clustering algorithm[J].IEEE.Transactions on Neural Networks,2005,16(3):645-67 8.
[4]YI Hong,SAM K.Learning assignment order of instances for the constrained k-means clustering algorithm[J].IEEE Transactions on Systems,Man,and Cybernetics,Part B:Cybernetics,2009,39(2):568-574.
[5]賀玲,吳玲達,蔡益朝.數據挖掘中的聚類算法綜述[J].計算機應用研究,2007,24(1):10-13.
[6]孫吉貴,劉杰,趙連宇.聚類算法研究[J].軟件學報,2008,19(1):48-61.
[7]孔英會,苑津莎,張鐵峰,等.基于數據流管理技術的配變負荷分類方法研究.中國國際供電會議,CICED2006.
[8]馬曉艷,唐雁.層次聚類算法研究[J].計算機科學,2008,34(7):34-36.
[9]汪海波,張海臣,段雪麗.基于MATLAB的自組織競爭神經網絡聚類研究[J].邢臺職業技術學院學報,2005,22(1):45-47.
[10]呂曉燕,羅立民,李祥生.FCM算法的改進及仿真實驗研究[J].計算機工程與應用,2009,45(20):144-147.
[11]李戈,邵峰晶,朱本浩.基于神經網絡聚類的研究[J].青島大學學報,2001,16(4):21-24.
[12]戈國華,肖海波,張敏.基于 FCM的數據聚類分析及matlab 實現[J].福建電腦,2007,4:89-90.
[13]FISHER R A.Iris Plants Database//http://www.ics.uci.edu/~mlearn/MLRepository.Html.Authorized license.