class?point（object）:
????x=0.0
????y=0.0
pointF=[]
pointType=[]#記錄點(diǎn)屬于的類
AverageD=[]#?記錄每個(gè)聚類的均值
ZArray=[]


StdDiff=[]???#?記錄聚類樣本中心標(biāo)準(zhǔn)差值
Std=[]???????#標(biāo)準(zhǔn)聚類中心
Sum=[]???????#求和臨時(shí)
N=[]??????????#記錄每個(gè)聚類書面

StdDistance=[]?#聚類中心之間距離
StdDisMax=[]
StdDisMaxCor=[]

MaxDiff=1????????#標(biāo)準(zhǔn)差判定區(qū)間?
MinDistance=2????#不同聚類中心最小距離
MaxNumStd=4??????#最大的聚類中心數(shù)目
TotalNum=10???????#點(diǎn)數(shù)

SAArray=[[]]
ZDistance=[]
ZDistanceR=[]
ZDistanceC=[]
StdTime=10
Nc=4
step=2?????????????#記錄步驟及當(dāng)前狀態(tài)
CountTime=0
#---------------------------------初始化
for?i?in?range（TotalNum）:
????pointF+=[point（）]
????pointType+=[0]
????StdDiff+=[point（）]
????ZDistance+=[0]
????ZDistanceR+=[0]
????ZDistanceC+=[0]
for?i?in?range（MaxNumStd）:
????AverageD+=[0]
????Std+=[point（）]
????Sum+=[point（）]
????ZArray+=[point（）]
????N+=[0]
for?i?in?range（MaxNumStd）:
????StdDistance+=[point（）]
????StdDisMax+=[0]
????StdDisMaxCor+=[0]
for?i?in?range（TotalNum）:
????SAArray+=[[]]
????for?j?in?range（TotalNum）:
????????SAArray[i]+=[point（）]
?
[pointF[0].xpointF[0].y]=[0.00.0]????
[pointF[1].xpointF[1].y]=[3.08.0]
[pointF[2].xpointF[2].y]=[2.02.0]
[pointF[3].xpointF[3].y]=[1.01.0]
[pointF[4].xpointF[4].y]=[5.03.0]
[pointF[5].xpointF[5].y]=[4.08.0]
[pointF[6].xpointF[6].y]=[6.03.0]
[pointF[7].xpointF[7].y]=[5.04.0]
[pointF[8].xpointF[8].y]=[6.04.0]????
[pointF[9].xpointF[9].y]=[7.05.0]
[ZArray[0].xZArray[0].y]=[00]
[ZArray[1].xZArray[1].y]=[38]
[ZArray[2].xZArray[2].y]=[22]
[ZArray[3].xZArray[3].y]=[11]
#----------------------------------函數(shù)定義區(qū)
def?DistancePoint（x1y1x2y2）:
????return?（（x1-x2）**2+（y1-y2）**2）**0.5
def?DistancePointF（ab）:
????return?（（a.x-b.x）**2+（a.y-b.y）**2）**0.5
def?ComputeNumStd（TypeList）:
????temp=pointType[7]
????for?i?in?range（6）:
????????if?temp????????????temp=pointType[i]
????return?temp+1
while（CountTime<=StdTime）:
????if?step==2:
????????for?i?in?range（Nc）:
????????????N[i]=0
????????print?‘這是第%d次歸類‘%CountTime
????????CountTime=CountTime+1
????????stdtemp=0
????????for?i?in?range（TotalNum）:
????????????dis=65535
????????????for?j?in?range（Nc）:
????????????????ftemp=DistancePointF（pointF[i]ZArray[j]）
????????????????if?ftemp????????????????????stdtemp=j
????????????????????dis=ftemp
????????????SAArray[stdtemp][N[stdtemp]].x=pointF[i].x
????????????SAArray[stdtemp][N[stdtemp]].y=pointF[i].y
????????????N[stdtemp]=N[stdtemp]+1
????????for?i?in?range（Nc）:
????????????print?“第%d個(gè)聚類中心是:（%d%d）擁有%d個(gè)元素???“%（iZArray[i].xZArray[i].yN[i]）
????????????print?“包含的元素有：“
????????????for?j?in?range（N[i]）:
????????????????print?“（%d%d）“%（SAArray[i][j].xSAArray[i][j].y）
????????step=3???#跳轉(zhuǎn)到第三步
????????#break
????if?step==3:
????????print“第%d步，判斷是否可以去掉一些“%step
????????for?i?in?range（Nc）:
????????????if?N[i]<1:?#1也可以為其他形參
????????????????#取消這個(gè)樣本子集
????????????????for?j?in?range（TotalNum）:
????????????????????if?pointType[j]==i:
????????????????????????poin