k-means python代码

k-means python代码

 

 

def fun_dis (x,y,n):
    return sum (map (lambda v1,v2:pow(abs(v1-v2),n), x,y))

def distance (x,y):
    #return fun_dis (x,y,2)
    return fun_dis (x,y,1)
    

def min_dis_center(center, node):
    min_index = 0
    min_dista = distance (center[0][1],node)
    for i in range (1,len(center)):
        tmp = distance (center[i][1],node)
        if (tmp < min_dista):
            min_dista = tmp
            min_index = i
    return min_index


# input [[x1,y1,z1..],[x2,y2,z2..]..]
# return 
def k_means (info,k=5):
    center = [[1,info[i]] for i in range(k)]
    result = [[i] for i in range(k)]
    width  = len (info[0])
    for i in range(k,len(info)):
        min_center = min_dis_center (center,info[i])
        for j in range(width):
            center[min_center][1][j] = (center[min_center][1][j] * center[min_center][0] + info[i][j])/ (1.0+center[min_center][0])
        center[min_center][0] += 1
        result[min_center].append (i)
    return result,center


test_data = [[1,100,300,9],
             [9,200,400,20],
             [2,102,305,10],
             [3,103,299,9],
             [3,105,299,11],
             [3,105,299,11],
             [9,200,400,20],
             [2,102,305,10],
             [11,220,409,20]]


data = k_means (test_data,2)
print data[0]
print data[1]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值