Advertisement

通过矩阵求两个向量集中元素两两之间的欧氏距离(python实现)

阅读量:

通过矩阵求两个向量集中元素两两之间的欧氏距离(python实现)

在许多算法中都涉及计算向量间的欧氏距离。例如,在机器学习领域中常用的KNN算法就需要考虑从训练集A和测试集B的所有向量对(Ai, Bi)进行分析,并计算这些向量对之间的欧氏距离。这样的处理会导致需要执行两个层次的循环结构,在处理大量数据时可能会降低效率。

这里介绍如何将这一过程用矩阵运算实现。

假设有两个三维向量集,用矩阵表示:

要求A,B两个集合中的元素两两间欧氏距离。

先求出AB^T

然后对AB^T分别求其中每个向量的模平方,并扩展为2*3矩阵:

然后:

将上面这个矩阵一开平方,就得到了A,B向量集两两间的欧式距离了。

下面是python实现:

复制代码
    import numpy
    def EuclideanDistances(A, B):
    BT = B.transpose()
    vecProd = A * BT
    SqA =  A.getA()**2
    sumSqA = numpy.matrix(numpy.sum(SqA, axis=1))
    sumSqAEx = numpy.tile(sumSqA.transpose(), (1, vecProd.shape[1]))    
    SqB = B.getA()**2
    sumSqB = numpy.sum(SqB, axis=1)
    sumSqBEx = numpy.tile(sumSqB, (vecProd.shape[0], 1))    
    SqED = sumSqBEx + sumSqAEx - 2*vecProd   
    ED = (SqED.getA())**0.5
    return numpy.matrix(ED)

全部评论 (0)

还没有任何评论哟~