基于ＶＢ與ＭａｔｒｉｘＶＢ的模糊Ｃ均值方法實(shí)現(xiàn)

劉　超　李茂寬　王偉華

摘 要：針對Visual Basic數(shù)值計算能力和圖像處理能力的不足，不利于系統(tǒng)開發(fā)，介紹了模糊C均值法的原理及其基于VB和MatrixVB的實(shí)現(xiàn);該方法將Matlab的強(qiáng)大計算功能與VB的Windows用戶界面開發(fā)方面優(yōu)勢結(jié)合起來，充分發(fā)揮各自優(yōu)勢，縮短了軟件的開發(fā)周期。軟件測試結(jié)果表明，該計算方法正確，軟件界面友好，計算速度快，系統(tǒng)資源消耗少，操作簡便易行，能滿足數(shù)據(jù)分類的要求。

關(guān)鍵詞：Visual Basic;MatrixVB;模糊C均值;模式識別

中圖分類號：TQ15.9;TP311 文獻(xiàn)標(biāo)識碼：B 文章編號：1004 373X(2009)02 102 02

Realization of Fuzzy C－Means Based on Visual Basic and MatrixVB

LIU Chao1,LI Maokuan2,WANG Weihua3

(1.PLA General Staff Meteorological and Hydrological Department,Beijing,100081,China;

2.Naval Aeronautical and Astronautical University,Yantai,264001,China;3.Unit 61828 of PLA,Kashi,844200,China)

Abstract：According to the situation that Visual Basic is insufficient in numerical computation,image processing ability and is weak in algorithm development.The theory of fuzzy C－means and the realization of fuzzy C－means by incorporating VBand MatrixVB are introduced.The incorporated advances of Matlab′s powerful calculation function and Visual Basic′s friendly graphic user interface shorten the period of software development by taking the advantages of both.It has been demonstrated by the instance that the calculating procedure is accurate.The software interface is friendly.Beside this,resource expenditure of computer is low and calculating speed is quick,the software is convenient to use.It can fulfill the data clustering problem well.

Keywords：Visual Basic;MatrixVB;fuzzy C－means;pattern recognition

0 引 言

聚類是將一組給定的未知類標(biāo)號的樣本分成內(nèi)在的多個類別，使得同一類中的樣本具有較高的相似度，而不同類中的樣本差別大。聚類分析的目的是揭示和刻畫數(shù)據(jù)的內(nèi)在結(jié)構(gòu)，其內(nèi)容涉及統(tǒng)計學(xué)、生物學(xué)、以及機(jī)器學(xué)習(xí)等研究領(lǐng)域，并在模式識別、數(shù)據(jù)分析和挖掘、圖像處理等領(lǐng)域獲得了廣泛的應(yīng)用［1］。模糊C均值（Fuzzy C－Means，F(xiàn)CM）算法是一種基于劃分的聚類算法，它的思想就是使得被劃分到同一類的對象之間相似度最大，而不同類之間的相似度最小。模糊C均值算法是普通C均值算法的改進(jìn)，普通C均值算法對于數(shù)據(jù)的劃分是硬性的，而FCM則是一種柔性的模糊劃分。

目前大多數(shù)針對模式識別的軟件都是基于Matlab，模糊C均值就已經(jīng)集成在Matlab語言中。Matlab是集數(shù)值分析、矩陣運(yùn)算、信號處理和圖形顯示于一體的高性能數(shù)學(xué)軟件，但Matlab也存在其自身的局限性：一般基于Matlab的應(yīng)用程序不能脫離Matlab集成環(huán)境工作，而且編寫界面的功能相對較弱，Matlab本身是一種解釋性的語言，運(yùn)行速度非常慢，大大降低了程序的效率。因此，將Matlab與VB有效地集成在一起，實(shí)現(xiàn)應(yīng)用系統(tǒng)的無縫集成，對于有效縮短開發(fā)周期，優(yōu)化系統(tǒng)性能是十分有意義的。介紹用VB與MatrixVB的混編來實(shí)現(xiàn)模糊C均值聚類方法。

1 關(guān)于MatrixVB 的簡介

MatrixVB［2］ 是Mathworks 公司針對VB 提供的一個Matlab 組件庫(Component Object Model Library) ，它提供了600 多個函數(shù)，包括基本的數(shù)學(xué)運(yùn)算和功能強(qiáng)大的信號處理、線性代數(shù)、串運(yùn)算及圖形圖像處理功能等，用來彌補(bǔ)VB 內(nèi)建函數(shù)的不足，主要是針對數(shù)學(xué)函數(shù)而言，它具有Matlab繪圖的強(qiáng)大功能，讓VB 能很輕易地畫出一些數(shù)學(xué)函數(shù)圖形，為VB 提供了強(qiáng)大的功能擴(kuò)展，讓程序員更容易地去開發(fā)計算應(yīng)用方面的程序代碼。在VB中使用該數(shù)學(xué)工具包可以避免重復(fù)性勞動，可以減少開發(fā)人員實(shí)現(xiàn)算法和界面設(shè)計方面的困難。

2 模糊C均值原理及其實(shí)現(xiàn)

模糊C均值（FCM），是用隸屬度確定每個數(shù)據(jù)點(diǎn)屬于某個聚類的程度的一種聚類算法。由Bezdek［3］提出了該算法，作為早期硬C均值聚類（HCM）方法的一種改進(jìn)。

FCM把n個向量x璱（i=1，2，…，n）分為c個模糊組，并求每組的聚類中心，使得非相似性指標(biāo)的價值函數(shù)達(dá)到最小。FCM與HCM的主要區(qū)別在于FCM用模糊劃分，使得每個給定數(shù)據(jù)點(diǎn)用值在［0，1］間的隸屬度來確定其屬于各個組的程度。與引入模糊劃分相適應(yīng)，隸屬矩陣U允許有取值在［0，1］間的元素。不過，加上歸一化規(guī)定，一個數(shù)據(jù)集的隸屬度和總等于1：

∑ci=1u璱j=1,衘=1,2,…,n(1)

則FCM的價值函數(shù)（或目標(biāo)函數(shù)）就是一般化形式：

J(U,c1,c2,…,c璫)=∑ci=1J璱=∑ci=1∑njum璱jd2璱j(2)

式中，c璱為模糊組i的聚類中心;d璱j=||c璱-x璲||為第i個聚類中心與第j個數(shù)據(jù)點(diǎn)間的歐幾里德距離；且m是一個加權(quán)指數(shù)。構(gòu)造如下新的目標(biāo)函數(shù)，可求得使式（2）達(dá)到最小值的必要條件：

(U,c1 ,c2,…,c璫 ,λ1,λ2,…,λ璶 )=

J(U,c1,c2,…,c璫 ) + ∑nj = 1λ璲 (∑ci = 1u璱j -1)=

∑ci = 1∑njum璱jd2璱j + ∑nj = 1(∑ci = 1u璱j -1)(3)

式中，λ璲，j=1,2,…,n是式（1）的n個約束式拉格朗日乘子。對所有輸入?yún)⒘壳髮?dǎo)，使式（2）達(dá)到最小的必要條件為：

c璱=∑nj=1um璱jx璲/∑nj=1um璱j(4)

和:

u璱j=1/∑ck=1d璱jd璳kj2/(m-1)(5)

根據(jù)上述兩個必要條件，模糊C均值聚類算法是一個簡單的迭代過程。在批處理方式運(yùn)行時，F(xiàn)CM用下列步驟確定聚類中心c璱和隸屬矩陣U ：

步驟1 用值在0，1間的隨機(jī)數(shù)初始化隸屬矩陣U，使其滿足式（1）中的約束條件。

步驟 2 用式（4）計算c個聚類中心c璱，i=1,2,…,c。

步驟3 用式（2）計算價值函數(shù)。如果它小于某個確定的閾值，或它相對上次價值函數(shù)值的改變量小于某個閾值，則算法停止。

步驟4 用式（5）計算新的U矩陣。返回步驟2。

通過上面的討論，不難看出FCM算法需要兩個參數(shù)，一個是聚類數(shù)目C；另一個是參數(shù)m。以下利用VB和MatrixVB實(shí)現(xiàn)上述算法。首先讀入待分類的數(shù)據(jù)，數(shù)據(jù)分布如圖1所示，是200組二維數(shù)據(jù)。設(shè)置參數(shù)：將數(shù)據(jù)分為3類，初始化隸屬度矩陣U，m=2。程序如下所示:

Cluster_N = 3

data_N = 200

Data = vbload("c:\data.mat")‘讀入數(shù)據(jù)并顯示

ax = vbaxes(Form1.hWnd)

aa = Data(linspace(1,200,200),1)

bb = Data(linspace(1,200,200),2)

Call scatter(aa,bb,".")

U = randM(Cluster_N,data_N) ‘初始化隸屬矩陣

tmp = Sum(U)

tmp = mtimes(ones(Cluster_N,1),tmp)

U = rdivide(U,tmp)

dist = zeros(3,200) ‘初始距離矩陣

算法的輸出是C個聚類中心點(diǎn)向量和大小為C×N的一個模糊劃分矩陣，這個矩陣表示的是每個樣本點(diǎn)屬于每個類的隸屬度。根據(jù)這個劃分矩陣，按照模糊集合中的最大隸屬原則就能夠確定每個樣本點(diǎn)歸為哪個類。聚類中心表示的是每個類的平均特征，可以認(rèn)為是這個類的代表點(diǎn)。以下是迭代求取聚類中心和隸屬度矩陣的程序：

mf = power(U,2)

center = rdivide(mtimes(mf,Data),

transpose(mtimes(ones(2,1),Sum(transpose(mf)))))

For j = 1 To 200

For i = 1 To 3

a1 = center(i,linspace(1,2,2))

a2 = Data(j,linspace(1,2,2))

dist(i,j) = Sum(power(minus(a1,a2),2))

Next i

Next j

obj_fuc = Sum(Sum(times(power(dist,2),mf)))

tmp = power(dist,-2)

U = rdivide(tmp,mtimes(ones(3,1),Sum(tmp)))

圖2是對未知類別的待分類數(shù)據(jù)進(jìn)行模糊C均值聚類的最終結(jié)果，不同的圖形表示不同的類別。

圖1 待分類數(shù)據(jù)

圖2 分類后數(shù)據(jù)

3 結(jié) 語

在此介紹VB與MatrixVB混合編程的模糊C均值方法實(shí)現(xiàn)，程序在 Windows XPProfessional 版上用MatrixVB4.5和Visual Basic6.0中文版實(shí)現(xiàn)，具有效率高，占用系統(tǒng)資源少等優(yōu)點(diǎn)。文中給出了部分基本代碼，可以根據(jù)自己的需要進(jìn)行擴(kuò)充，并應(yīng)用于分類處理。這種編程方式既節(jié)省了時間，又提高了軟件的性能，對于優(yōu)化系統(tǒng)，縮短軟件開發(fā)周期很有意義。

參考文獻(xiàn)

［1］邊肇祺，張學(xué)工.模式識別［M］.2版.北京：清華大學(xué)出版社,2000.

［2］MathWorks.Inc.MatrixVB User′s Guide［Z］.2000.

［3］Bezdek J C.Pattern Recognition with Objective Function Algorithms ［M］.New York:Plenum Press,1981.

［4］王躍強(qiáng),王紀(jì)龍,王云才.VB 程序中實(shí)現(xiàn)調(diào)用Matlab 的方法［J］.計算機(jī)應(yīng)用,2002,2(21):95－96.

［5］何強(qiáng),何英.Matlab 擴(kuò)展編程［M］.北京:清華大學(xué)出版社,2002.

［6］劉炳文.精通Visual Basic6.0中文版［M］.北京:電子工業(yè)出版社,1999.

［7］黃錫泉.VB與Matlab 無縫接口編程［J］.微計算機(jī)應(yīng)用,2005,26(2):238－239.

［8］董樑,李光軍,賀仲雄.復(fù)雜系統(tǒng)的集對模式識別方法.現(xiàn)代電子技術(shù),2007(2):68－70.

作者簡介 劉 超 男，1981年出生，山東青島人，碩士研究生。研究方向?yàn)橛嬎銠C(jī)仿真。

李茂寬 男，1979年出生，山東即墨人，博士研究生。研究方向?yàn)閳D像處理與模式識別。

王偉華 男，1974年出生，河南人。研究方向?yàn)橛嬎銠C(jī)仿真。