matlab pdist2函数官方详解

pdist2

两组观测之间的成对距离


句法

D = pdist2(X,Y,Distance)

D = pdist2(X,Y,Distance,DistParameter)

D = pdist2(___,Name,Value)

[D,I] = pdist2(___,Name,Value)


描述

D = pdist2(X,Y,Distance)中的每一对观测之间的距离。XY使用Distance.

D = pdist2(X,Y,Distance,DistParameter)指定的度量返回距离。DistanceDistParameter。您可以指定DistParameter只有当Distance'seuclidean''minkowski',或'mahalanobis'.

D = pdist2(___,Name,Value)使用名称-值对参数之一指定附加选项。'Smallest''Largest'除了以前语法中的任何参数之外。

例如,

  • D = pdist2(X,Y,Distance,'Smallest',K)指定的度量计算距离。Distance并返回K与观测的最小成对距离X中的每一次观察Y按升序排列。

  • D = pdist2(X,Y,Distance,DistParameter,'Largest',K)指定的度量计算距离。DistanceDistParameter并返回K按降序排列的最大成对距离。

[D,I] = pdist2(___,Name,Value)还返回矩阵。I。矩阵I中的观测结果的索引。X对应于D.


实例

计算欧氏距离

用三个观察和两个变量创建两个矩阵。

rng('default') % For reproducibility
X = rand(3,2);
Y = rand(3,2);

计算欧氏距离。输入参数的默认值。Distance'euclidean'。当不使用名称-值对参数计算欧几里德距离时,不需要指定Distance.

D = pdist2(X,Y)
D = 3×3

    0.5387    0.8018    0.1538
    0.7100    0.5951    0.3422
    0.8805    0.4242    1.2050

D(i,j)对应于观察之间的成对距离i在……里面X和观察j在……里面Y.

计算Minkowski距离

用三个观察和两个变量创建两个矩阵。

rng('default') % For reproducibility
X = rand(3,2);
Y = rand(3,2);

用默认指数2计算Minkowski距离。

D1 = pdist2(X,Y,'minkowski')
D1 = 3×3

    0.5387    0.8018    0.1538
    0.7100    0.5951    0.3422
    0.8805    0.4242    1.2050

计算Minkowski距离,其指数为1,等于城市块距离。

D2 = pdist2(X,Y,'minkowski',1)
D2 = 3×3

    0.5877    1.0236    0.2000
    0.9598    0.8337    0.3899
    1.0189    0.4800    1.7036

D3 = pdist2(X,Y,'cityblock')
D3 = 3×3

    0.5877    1.0236    0.2000
    0.9598    0.8337    0.3899
    1.0189    0.4800    1.7036

找到两个最小的成对距离

用三个观察和两个变量创建两个矩阵。

rng('default') % For reproducibility
X = rand(3,2);
Y = rand(3,2);

找到两个最小的两两欧几里得距离X中的每一次观察Y.

[D,I] = pdist2(X,Y,'euclidean','Smallest',2)
D = 2×3

    0.5387    0.4242    0.1538
    0.7100    0.5951    0.3422

I = 2×3

     1     3     1
     2     2     2

中的每一次观察Ypdist2通过计算并将距离值与X。的每个列中的距离排序。D按升序排列。I中的观测结果的索引。X对应于D.

使用自定义距离函数计算缺少元素的成对距离

定义一个自定义距离函数,该函数忽略NaN值,并使用自定义距离函数计算成对的距离。

用三个观察和三个变量创建两个矩阵。

rng('default') % For reproducibility
X = rand(3,3)
Y = [X(:,1:2) rand(3,1)]
X =

    0.8147    0.9134    0.2785
    0.9058    0.6324    0.5469
    0.1270    0.0975    0.9575


Y =

    0.8147    0.9134    0.9649
    0.9058    0.6324    0.1576
    0.1270    0.0975    0.9706

X和Y的前两列是相同的。假设X(1,1)丢失了。

X(1,1) = NaN
X =

       NaN    0.9134    0.2785
    0.9058    0.6324    0.5469
    0.1270    0.0975    0.9575

计算汉明距离。

D1 = pdist2(X,Y,'hamming')
D1 =

       NaN       NaN       NaN
    1.0000    0.3333    1.0000
    1.0000    1.0000    0.3333

中频观测i在……里面X或观察j在……里面YNaN值,函数pdist2回报NaN对于两两之间的距离ij。因此,D1(1,1),D1(1,2)和D1(1,3)是NaN价值。

定义自定义距离函数nanhamdist忽略坐标的NaN值并计算Hamming距离。当处理大量的观测时,您可以通过遍历数据的坐标来更快地计算距离。

function D2 = nanhamdist(XI,XJ)  
%NANHAMDIST Hamming distance ignoring coordinates with NaNs
[m,p] = size(XJ);
nesum = zeros(m,1);
pstar = zeros(m,1);
for q = 1:p
    notnan = ~(isnan(XI(q)) | isnan(XJ(:,q)));
    nesum = nesum + ((XI(q) ~= XJ(:,q)) & notnan);
    pstar = pstar + notnan;
end
D2 = nesum./pstar; 

nanhamdist的输入参数传递函数句柄。pdist2.

D2 = pdist2(X,Y,@nanhamdist)
D2 =

    0.5000    1.0000    1.0000
    1.0000    0.3333    1.0000
    1.0000    1.0000    0.3333

向现有集群分配新数据并生成C/C++代码

此示例使用:

执行k-均值聚类

使用三个发行版生成培训数据集。

rng('default') % For reproducibility
X = [randn(100,2)*0.75+ones(100,2);
    randn(100,2)*0.5-ones(100,2);
    randn(100,2)*0.75];

使用kmeans.

[idx,C] = kmeans(X,3);

绘制星系团和星系团质心图。

figure
gscatter(X(:,1),X(:,2),idx,'bgm')
hold on
plot(C(:,1),C(:,2),'kx')
legend('Cluster 1','Cluster 2','Cluster 3','Cluster Centroid')

matlab pdist2函数官方详解

将新数据分配给现有集群

生成测试数据集。

Xtest = [randn(10,2)*0.75+ones(10,2);
    randn(10,2)*0.5-ones(10,2);
    randn(10,2)*0.75];

使用现有集群对测试数据集进行分类。从每个测试数据点查找最近的质心,使用pdist2.

[~,idx_test] = pdist2(C,Xtest,'euclidean','Smallest',1);

绘制测试数据并使用idx_testgscatter.

gscatter(Xtest(:,1),Xtest(:,2),idx_test,'bgm','ooo')
legend('Cluster 1','Cluster 2','Cluster 3','Cluster Centroid', ...
    'Data classified to Cluster 1','Data classified to Cluster 2', ...
    'Data classified to Cluster 3')

matlab pdist2函数官方详解

 

输入参数

X,Y — 输入数据
数值矩阵

输入数据,指定为数字矩阵。XMX-按-n矩阵和Y我的-按-n矩阵。行对应于单独的观察,列对应于单个变量。

数据类型:single | double

Distance — 距离度量
字符向量 | 串标量 | 功能手柄

距离度量,指定为字符向量、字符串标量或函数句柄,如下表所述。

价值 描述
'euclidean'

欧几里德距离(默认)。

'squaredeuclidean'

平方欧氏距离(本选项仅为提高效率而提供。它不满足三角形不等式。

'seuclidean'

标准化欧氏距离观测值之间的每个坐标差除以标准差的对应元素,S = nanstd(X)。使用DistParameter指定另一个值S.

'mahalanobis'

使用样本协方差的Mahalanobis距离XC = nancov(X)。使用DistParameter指定另一个值C,其中的矩阵C是对称的正定的。

'cityblock'

城市街区距离。

'minkowski'

明考斯基距离。默认指数为2。DistParameter指定不同的指数P,在哪里P是指数的正标量值。

'chebychev'

切比切夫距离(最大坐标差)

'cosine'

1减去点间夹角的余弦(作为向量处理)。

'correlation'

1减去点之间的样本相关性(作为值的序列处理)。

'hamming'

Hamming距离,即不同坐标的百分比。

'jaccard'

1减去Jaccard系数,这是不同的非零坐标的百分比。

'spearman'

其中一项减去样本Spearman之间的秩相关性(作为值的序列处理)。

@distfun

自定义距离函数句柄。距离函数具有以下形式

function D2 = distfun(ZI,ZJ)
% calculation of distance
...
哪里

 

  • ZI1-按-n包含单一观测的向量。

  • ZJm2-按-n包含多个观察的矩阵。distfun必须接受矩阵ZJ任意数量的观察。

  • D2m2-按-1距离向量,以及D2(k)是观察之间的距离ZIZJ(k,:).

 

如果数据不稀疏,通常可以使用内置距离而不是函数句柄来更快地计算距离。

有关定义,请参见距离度量.

当你使用'seuclidean''minkowski',或'mahalanobis',您可以指定一个附加的输入参数。DistParameter来控制这些指标。您还可以使用这些指标的方式与默认值为DistParameter.

例子: 'minkowski'

DistParameter — 距离度量参数值
正标量 | 数值向量 | 数值矩阵

距离度量参数值,指定为正标量、数字向量或数字矩阵。此参数仅在您指定Distance'seuclidean''minkowski',或'mahalanobis'.

 

  • 如果Distance'seuclidean'DistParameter是每个维度的缩放因子的向量,指定为正向量。默认值是nanstd(X).

  • 如果Distance'minkowski'DistParameter是Minkowski距离的指数,表示为正标量。默认值为2。

  • 如果Distance'mahalanobis'DistParameter是一个协方差矩阵,指定为数字矩阵。默认值是nancov(X)DistParameter必须是对称的和正定的。

 

例子: 'minkowski',3

数据类型:single | double

名称值对参数

指定可选的逗号分隔对Name,Value争论。Name是参数名和Value对应的值。Name必须出现在引号中。可以将多个名称和值对参数按任何顺序指定为Name1,Value1,...,NameN,ValueN.

例子:任一'Smallest',K'Largest',K。你不能同时使用'Smallest''Largest'.

'Smallest' — 要查找的最小距离数
正整数

要查找的最小距离数,指定为由'Smallest'以及一个正整数。如果您指定'Smallest',然后pdist2的每一列中的距离排序。D按升序排列。

例子:'Smallest',3

数据类型:single | double

'Largest' — 要查找的最大距离数
正整数

要查找的最大距离数,指定为由'Largest'以及一个正整数。如果您指定'Largest',然后pdist2的每一列中的距离排序。D按降序排列。

例子:'Largest',3

数据类型:single | double

D-成对距离
数值矩阵

成对距离,返回为数字矩阵。

如果您没有指定'Smallest''Largest',然后DMX-按-我的矩阵,其中MX我的中的观察数XY分别。D(i,j)是观察之间的距离i在……里面X和观察j在……里面Y。中频观测i在……里面X或观察j在……里面YNaN,然后D(i,j)NaN内建距离函数。

如果您指定'Smallest''Largest'K,然后DK-按-我的矩阵。D包含K最小的或最小的K与观测的最大成对距离X中的每一次观察Y。中的每一次观察Ypdist2发现K的最小或最大距离,通过计算和比较的距离值与所有的观测X。如果K大于MXpdist2返回MX-按-我的矩阵。

I-排序索引
正整数矩阵

排序索引,作为正整数矩阵返回。I是和DI中的观测结果的索引。X对应于D.