在做分类时常常需要估算不同样本之间的相似性度量(Similarity Measurement),
这时通常采用的方法就是计算样本间的“距离”(Distance)。
采用什么样的方法计算距离是很讲究,甚至关系到分类的正确与否。
本文的目的就是对常用的相似性度量作一个总结
本文目录:
1. 欧氏距离
2. 曼哈顿距离
3. 切比雪夫距离
4. 闵可夫斯基距离
5. 标准化欧氏距离
6. 马氏距离
7. 夹角余弦
8. 汉明距离
9. 杰卡德距离 & 杰卡德相似系数
10. 相关系数 & 相关距离
11. 信息熵
1. 欧氏距离(Euclidean Distance)
欧氏距离是最易于理解的一种距离计算方法,源自欧氏空间中两点间的距离公式。
(1)二维平面上两点a(x1,y1)与b(x2,y2)间的欧氏距离:
(2)三维空间两点a(x1,y1,z1)与b(x2,y2,z2)间的欧氏距离:
(3)两个n维向量a(x11,x12,…,x1n)与 b(x21,x22,…,x2n)间的欧氏距离:
也可以用表示成向量运算的形式:
(4)计算欧氏距离
java
* 欧几里得度量算法实现 public class Euclidean { /** * 两个向量可以为任意维度,但必须保持维度相同,表示n维度中的两点 */ public static void main(String[] args) { double distance = 0; double[] vector1= {1.1,1.2}; double[] vector2= {2.1,2.2}; if (vector1.length == vector2.length) { for (int i = 0; i < vector1.length; i++) { double temp = Math.pow((vector1[i] - vector2[i]), 2); distance += temp; } distance = Math.sqrt(distance); } System.out.println(distance); } }
python
import numpy as np
vector1 = np.mat([1,2,3])
vector2= np.mat([4,5,6])
print(np.sqrt((vector1-vector2)*((vector1-vector2).T))) #向量
print(np.power(np.sum(np.power(np.abs(vector1-vector2),2)),1/2)) #原始公式
2. 曼哈顿距离(Manhattan Distance)
从名字就可以猜出这种距离的计算方法了。想象你在曼哈顿要从一个十字路口开车到另外一个十字路口,驾驶距离是两点间的直线距离吗?显然不是,除非你能穿越大楼。实际驾驶距离就是这个“曼哈顿距离”。而这也是曼哈顿距离名称的来源, 曼哈顿距离也称为城市街区距离(City Block distance)。
(1)二维平面两点a(x1,y1)与b(x2,y2)间的曼哈顿距离
(2)两个n维向量a(x11,x12,…,x1n)与 b(x21,x22,…,x2n)间的曼哈顿距离
(3)计算曼哈顿距离
java
public class Manhattan { public static void main(String[] args) { double distance = 0; double[] vector1= {1,2,3}; double[] vector2= {4,5,6}; if (vector1.length == vector2.length) { for (int i = 0; i < vector1.length; i++) { double temp = Math.abs(vector1[i] - vector2[i]); distance += temp; } } System.out.println(distance); } }
python
import numpy as np vector1=np. mat([1,2,3]) vector2 = np.mat([4,5,6]) print(np.sum(np.abs(vector1-vector2)))
3. 切比雪夫距离 ( Chebyshev Distance )
国际象棋玩过么?国王走一步能够移动到相邻的8个方格中的任意一个。那么国王从格子(x1,y1)走到格子(x2,y2)最少需要多少步?自己走走试试。你会发现最少步数总是max( | x2-x1 | , | y2-y1 | ) 步 。有一种类似的一种距离度量方法叫切比雪夫距离。
(1)二维平面两点a(x1,y1)与b(x2,y2)间的切比雪夫距离
(2)两个n维向量a(x11,x12,…,x1n)与 b(x21,x22,…,x2n)间的切比雪夫距离
这个公式的另一种等价形式是
看不出两个公式是等价的?提示一下:试试用放缩法和夹逼法则来证明。
(3c
public class Chebyshev{ public static void main(String[] args) { double distance = 0; double[] vector1= {1,2,3}; double[] vector2= {4,5,7}; if (vector1.length == vector2.length) { for (int i = 0; i < vector1.length; i++) { double temp = Math.abs(vector1[i] - vector2[i]); if (temp>distance) { distance = temp; } } } System.out.println(distance); } }
import numpy as np vector1=np. mat([1,2,3]) vector2 = np.mat([4,5,6]) print(np.max(np.abs(vector1-vector2)))
4. 闵可夫斯基距离(Minkowski Distance)
闵氏距离不是一种距离,而是一组距离的定义。
(1) 闵氏距离的定义
两个n维变量a(x11,x12,…,x1n)与 b(x21,x22,…,x2n)间的闵可夫斯基距离定义为:
其中p是一个变参数。
当p=1时,就是曼哈顿距离
当p=2时,就是欧氏距离
当p→∞时,就是切比雪夫距离
根据变参数的不同,闵氏距离可以表示一类的距离。
(2)闵氏距离的缺点
闵氏距离,包括曼哈顿距离、欧氏距离和切比雪夫距离都存在明显的缺点。
举个例子:二维样本(身高,体重),其中身高范围是150~190,体重范围是50~60,有三个样本:a(180,50),b(190,50),c(180,60)。那么a与b之间的闵氏距离(无论是曼哈顿距离、欧氏距离或切比雪夫距离)等于a与c之间的闵氏距离,但是身高的10cm真的等价于体重的10kg么?因此用闵氏距离来衡量这些样本间的相似度很有问题。
简单说来,闵氏距离的缺点主要有两个:(1)将各个分量的量纲(scale),也就是“单位”当作相同的看待了。(2)没有考虑各个分量的分布(期望,方差等)可能是不同的。
(3)计算闵氏距离
public class Minkowski{ public static void main(String[] args) { double distance = 0; double[] vector1= {1,2,3}; double[] vector2= {4,5,6}; if (vector1.length == vector2.length) { for(int p=1;p<=10;p++) { if(p==1){ for (int i = 0; i < vector1.length; i++) { double temp = Math.abs(vector1[i] - vector2[i]); distance += temp; } } if(p>1) { distance = 0; for (int i = 0; i < vector1.length; i++) { double temp = Math.pow(Math.abs(vector1[i] - vector2[i]),p); distance += temp; } distance = Math.pow(distance,1.0/p); } System.out.println(distance); } } } }
import numpy as np vector1=np. mat([1,2,3]) vector2 = np.mat([4,5,6]) for p in range(1,10) : if p == 1 : print(np.sum(np.abs(vector1-vector2))) # elif p == 2 : # print(np.sqrt((vector1-vector2)*((vector1-vector2).T))) else: print(np.power(np.sum(np.power(np.abs(vector1-vector2),p)),1/p))
5. 标准化欧氏距离 (Standardized Euclidean distance )
(1)标准欧氏距离的定义
标准化欧氏距离是针对简单欧氏距离的缺点而作的一种改进方案。标准欧氏距离的思路:既然数据各维分量的分布不一样,好吧!那我先将各个分量都“标准化”到均值、方差相等吧。均值和方差标准化到多少呢?这里先复习点统计学知识吧,假设样本集X的均值(mean)为m,标准差(standard deviation)为s,那么X的“标准化变量”表示为:
而且标准化变量的数学期望为0,方差为1。因此样本集的标准化过程(standardization)用公式描述就是:
标准化后的值 = ( 标准化前的值 - 分量的均值 ) /分量的标准差
经过简单的推导就可以得到两个n维向量a(x11,x12,…,x1n)与 b(x21,x22,…,x2n)间的标准化欧氏距离的公式:
如果将方差的倒数看成是一个权重,这个公式可以看成是一种加权欧氏距离(Weighted Euclidean distance)。
(2)计算标准化欧氏距离
public class StandardizedEuclideanDistance { public static void main(String[] args) { double distance = 0; double[] vector1= {1,2,3}; double[] vector2= {4,5,6}; double []s=new double[vector2.length+1]; double []avg=new double[vector2.length]; // 均值 if (vector1.length == vector2.length) { for(int i=0;i<vector2.length;i++){ avg[i]=( vector1[i]+vector2[i])/2; //均值 s[i]= Math.sqrt((Math.pow( vector1[i]-avg[i],2)+Math.pow( vector2[i]-avg[i],2))/2); //标准差 // s[i]= (Math.pow( vector1[i]-avg[i],2)+Math.pow( vector2[i]-avg[i],2))/2; 也可以使用方差 //标准化欧氏距离 double temp= Math.pow((vector1[i] - vector2[i])/s[i],2); distance = distance+temp; //System.out.println( s[i]+","+temp+","+distance); } } double y = distance; // System.out.println(Math.pow(y,0.5)); System.out.println( Math.sqrt(y)); } }
import numpy as np x=np. mat([1,2,3]) y=np.mat([4,5,7]) #合并 """ stack()函数 函数原型:stack(arrays, axis=0) axis=0时"增加一维,新维度的下标为0 即按行合并; axis=0时"增加一维,新维度的下标为1" 即按行合并; 等同于hstack()函数 函数原型:hstack(tup) ,参数tup可以是元组,列表, 或者numpy数组,返回结果为numpy的数组 作用:在水平方向把元素堆起来 vstack()函数 函数原型:vstack(tup) ,参数tup可以是元组,列表,或者numpy数组,返回结果为numpy的数组 作用:在垂直方向把元素堆叠起来 """ X=np.vstack([x,y]) print(X) #方法一:根据公式求解 sk=np.var(X,axis=0,ddof=0)#方差 #numpy.var(a, axis=None, dtype=None, out=None, ddof=0, keepdims=<class numpy._globals._NoValue>) print(sk) d1=np.sqrt(((np.power((x-y),2)/sk).sum())) print(d1) #方法二:根据scipy库求解 """ pdist """ from scipy.spatial.distance import pdist d2=pdist(X,'seuclidean',V=None) print(d3) """ scipy.spatial.distance.seuclidean scipy.spatial.distance.seuclidean(u, v, V) 三个参数 Returns the standardized Euclidean distance between two 1-D arrays. The standardized Euclidean distance between u and v. Parameters: u : (N,) array_like Input array. v : (N,) array_like Input array. V : (N,) array_like 方差 V is an 1-D array of component variances. It is usually computed among a larger collection vectors. Returns: seuclidean : double The standardized Euclidean distance between vectors u and v. """ from scipy.spatial.distance import seuclidean d3= seuclidean(x, y, sk) print(d3)
6. 马氏距离(Mahalanobis Distance)
(1)马氏距离定义
有M个样本向量X1~Xm,协方差矩阵记为S,均值记为向量μ,则其中样本向量X到u的马氏距离表示为:
而其中向量Xi与Xj之间的马氏距离定义为:
若协方差矩阵是单位矩阵(各个样本向量之间独立同分布),则公式就成了:
也就是欧氏距离了。
若协方差矩阵是对角矩阵,公式变成了标准化欧氏距离。
(2)马氏距离的优缺点:量纲无关,排除变量之间的相关性的干扰。
(3)
package aa; import org.apache.commons.math3.linear.Array2DRowRealMatrix; import org.apache.commons.math3.linear.LUDecomposition; import org.apache.commons.math3.linear.RealMatrix; public class MahalanobisDistance { public static void main(String[] args){ double[] x0={10,2,3}; double[] x1={1,30,1}; double[] x2={200,200,5000}; double[] x3={3000,1000,1000}; double[][] x = {{10,2,3},{1,30,1},{200,200,5000},{3000,1000,1000},{400000,400000,8}}; double[][] cov = new double[5][5] ; double [] mean = new double[5] ; for (int i = 0; i < 5; i++) { for(int j = 0;j<3;j++) { mean[i] += x[i][j]/3;//均值 } for(int j = 0;j<5;j++) { for(int k = 0;k<3;k++) { cov[i][j] += (x[i][k]-mean[i])*(x[j][k]-mean[j])/2;//协方差 } System.out.println(i+",,"+j+":"+cov[i][j]); } } double[] x01 =new double [3]; double[] x02 =new double [3]; double[] x03 =new double [3]; for (int i = 0; i < 3; i++) { x01[i]=x0[i]-x1[i]; x02[i]=x0[i]-x2[i]; x03[i]=x0[i]-x3[i]; } RealMatrix matrix = new Array2DRowRealMatrix(cov); //将数组转化为矩阵 System.out.println(matrix); RealMatrix matrix1 = new Array2DRowRealMatrix(x01); //将数组转化为矩阵 RealMatrix matrix2 = matrix1.transpose(); //矩阵转置 RealMatrix inversetestMatrix = inverseMatrix(matrix);//求逆 //求MahalanobisDistance 并将矩阵转化为数组 RealMatrix detNum = matrix2.multiply(matrix1).multiply(inversetestMatrix); System.out.println("MahalanobisDistance:\t"+detNum); } public static RealMatrix inverseMatrix(RealMatrix A) { RealMatrix result = new LUDecomposition(A).getSolver().getInverse(); return result; } }注:Java 运行过程中报错,目前问题 没有解决!!!后期更新维护
X1=np.mat([[10,2,3],[1,30,1],[200,200,5000],[3000,1000,1000],[400000,400000,8]]) Y = pdist(X1,'mahal',V=None) y2=np.cov(X1) print(y2) print(Y)
7. 夹角余弦(Cosine)
有没有搞错,又不是学几何,怎么扯到夹角余弦了?各位看官稍安勿躁。几何中夹角余弦可用来衡量两个向量方向的差异,机器学习中借用这一概念来衡量样本向量之间的差异。
(1)在二维空间中向量A(x1,y1)与向量B(x2,y2)的夹角余弦公式:
(2) 两个n维样本点a(x11,x12,…,x1n)和b(x21,x22,…,x2n)的夹角余弦
类似的,对于两个n维样本点a(x11,x12,…,x1n)和b(x21,x22,…,x2n),可以使用类似于夹角余弦的概念来衡量它们间的相似程度。
即:
夹角余弦取值范围为[-1,1]。夹角余弦越大表示两个向量的夹角越小,夹角余弦越小表示两向量的夹角越大。当两个向量的方向重合时夹角余弦取最大值1,当两个向量的方向完全相反夹角余弦取最小值-1。
夹角余弦的具体应用可以参阅参考文献[1]。
(3)计算夹角余弦
8. 汉明距离(Hamming distance)
(1)汉明距离的定义
两个等长字符串s1与s2之间的汉明距离定义为将其中一个变为另外一个所需要作的最小替换次数。例如字符串“1111”与“1001”之间的汉明距离为2。
应用:信息编码(为了增强容错性,应使得编码间的最小汉明距离尽可能大)。
(2)计算汉明距离
9. 杰卡德相似系数(Jaccard similarity coefficient)
(1) 杰卡德相似系数
两个集合A和B的交集元素在A,B的并集中所占的比例,称为两个集合的杰卡德相似系数,用符号J(A,B)表示。
杰卡德相似系数是衡量两个集合的相似度一种指标。
(2) 杰卡德距离
与杰卡德相似系数相反的概念是杰卡德距离(Jaccard distance)。杰卡德距离可用如下公式表示:
杰卡德距离用两个集合中不同元素占所有元素的比例来衡量两个集合的区分度。
(3) 杰卡德相似系数与杰卡德距离的应用
可将杰卡德相似系数用在衡量样本的相似度上。
样本A与样本B是两个n维向量,而且所有维度的取值都是0或1。例如:A(0111)和B(1011)。我们将样本看成是一个集合,1表示集合包含该元素,0表示集合不包含该元素。
p :样本A与B都是1的维度的个数
q :样本A是1,样本B是0的维度的个数
r :样本A是0,样本B是1的维度的个数
s :样本A与B都是0的维度的个数
那么样本A与B的杰卡德相似系数可以表示为:
这里p+q+r可理解为A与B的并集的元素个数,而p是A与B的交集的元素个数。
而样本A与B的杰卡德距离表示为:
(4)计算杰卡德距离
10. 相关系数 ( Correlation coefficient )与相关距离(Correlation distance)
(1) 相关系数的定义
相关系数是衡量随机变量X与Y相关程度的一种方法,相关系数的取值范围是[-1,1]。相关系数的绝对值越大,则表明X与Y相关度越高。当X与Y线性相关时,相关系数取值为1(正线性相关)或-1(负线性相关)。
(2)相关距离的定义
(3)计算相关系数与相关距离
11. 信息熵(Information Entropy)
信息熵是衡量分布的混乱程度或分散程度的一种度量。分布越分散(或者说分布越平均),信息熵就越大。分布越有序(或者说分布越集中),信息熵就越小。
计算给定的样本集X的信息熵的公式:
参数的含义:
n:样本集X的分类数
pi:X中第i类元素出现的概率
信息熵越大表明样本集S分类越分散,信息熵越小则表明样本集X分类越集中。。当S中n个分类出现的概率一样大时(都是1/n),信息熵取最大值log2(n)。当X只有一个分类时,信息熵取最小值0