【Python经典题目】闵可夫斯基距问题

题目

定义一个高维空间样本点集类HDPoints，须包含以下数据属性与方法属性：

(a)数据属性self.points：类型为列表，由多个子列表构成，每个子列表表示高维空间中的一个数据点，且数据维度可以任意，并通过初始化构造函数获得。

(b)方法属性centerpoint(self)：计算点集的中心点。

(c)方法属性minkowski (self, x, y, p)：计算两点x和y之间的闵可夫斯基距离，p为非负整数，用p=0情形表示切比雪夫距离。由此定义的距离称为p-闵氏距离，其数学定义如下：

(d)方法属性farthestpoint(self, p)：找出离中心点p-闵氏距离最远的点，返回在self.points中的下标以及最大距离。

(e)方法属性farthest2points(self, p)：找出点集self.points中p-闵氏距离最远的两点，返回两点在self.points中的下标及其最大距离。

要求

接下来，实例化类HDPoints，利用random模块，随机产生至少50个高维空间数据点，样本点的维度至少在5以上，且每个分量取值服从区间[0,1]上的均匀分布。同时，随机产生一个0~5之间的一个非负整数，赋值传递给p-闵氏距离函数中的参数p，对HDPoints实例对象的全部自定义方法属性（即centerpoint()、minkowski()、farthestpoint()和farthest2points()）进行功能测试。

代码

import numpy as np

class HDPoints(object):
    def __init__(self,list_HDPoints):
        self.points=list_HDPoints

    def centerpoint(self):
        # 转化为高维数组,数组要求长度一致
        m=np.array(self.points)
        return sum(m) / len(m)

    def minkowski (self, x, y, p):
        if (p==0):p=1
        return sum([abs(x - y) ** p for x, y in zip(x, y)])**(1/p)

    def farthestpoint(self, p):
        centerPoint=self.centerpoint()
        list_instance=[self.minkowski(centerPoint,i,p)for i in self.points]
        max1=max(list_instance)
        return list_instance.index(max1),max1

    def farthest2points(self, p):
        max=0
        l1,l2=[],[]
        for i in self.points:
            for j in self.points[self.points.index(i)+1::]:
                if self.minkowski(i,j,p)>max:
                    max=self.minkowski(i,j,p)
                    l1,l2=i,j
        return max,self.points.index(l1),self.points.index(l2)

#产生50个高维点
a=np.random.uniform(0,5,(50,5))
m=[]
for i in a:
    m1=[]
    for j in i:
        m1.append(int(j))
    m.append(m1)
print(m)
HD_Point=HDPoints(m)

print("中心点是:",HD_Point.centerpoint())
print("两点之间的闵可夫斯基距离是:",HD_Point.minkowski([3,4,5,6],[2,3,4,5],4))
print("当p=5时,离中心点距离最远的点和下标是:",HD_Point.farthestpoint(5))
print("当p=5时，距离最远的两个点下标和距离是:",HD_Point.farthest2points(5))