prctile

数据集的百分位数

全页折叠

语法

P = prctile(A,p)

P = prctile(A,p,"all")

P = prctile(A,p,dim)

P = prctile(A,p,vecdim)

P = prctile(___,Method=method)

说明

P = prctile(A,p) 返回输入数据 A 在区间 [0,100] 中百分比 p 所对应的百分位数。

如果 A 是向量，则 P 是标量或是与 p 长度相同的向量。P(i) 包含 p(i) 百分位数。
如果 A 是矩阵，则 P 是行向量或矩阵，其中 P 的行数等于 length(p)。P 的第 i 行包含 A 的每列的 p(i) 百分位数。
如果 A 是多维数组，则 P 包含沿其大小不等于 1 的第一个数组维度计算的百分位数。

示例

P = prctile(A,p,"all") 返回 x 中所有元素的百分位数。

示例

P = prctile(A,p,dim) 沿维度 dim 进行运算。例如，如果 A 是矩阵，则 prctile(A,p,2) 对每行中的元素进行运算。

示例

P = prctile(A,p,vecdim) 沿在向量 vecdim 中指定的维度进行运算。例如，如果 A 是矩阵，则 prctile(A,p,[1 2]) 对 A 的所有元素进行运算，因为矩阵的每个元素都包含在由维度 1 和 2 定义的数组切片中。

示例

P = prctile(___,Method=method) 使用指定的方法计算百分位数。除了上述语法中的任何输入参量组合之外，还需要指定方法。

示例

全部折叠

数据向量的百分位数

打开实时脚本

计算给定百分比的数据集的百分位数。

生成一个大小为 7 的数据集。

rng default % for reproducibility
A = randn(1,7)

A = 1×7

    0.5377    1.8339   -2.2588    0.8622    0.3188   -1.3077   -0.4336

计算 A 的元素的第 42 个百分位数。

P = prctile(A,42)

P = 
-0.1026

所有值的百分位数

打开实时脚本

计算一个数组中所有值的百分位数。

创建一个 3×5×2 数组。

rng default % for reproducibility
A = randn(3,5,2)

A = 
A(:,:,1) =

    0.5377    0.8622   -0.4336    2.7694    0.7254
    1.8339    0.3188    0.3426   -1.3499   -0.0631
   -2.2588   -1.3077    3.5784    3.0349    0.7147


A(:,:,2) =

   -0.2050    1.4090   -1.2075    0.4889   -0.3034
   -0.1241    1.4172    0.7172    1.0347    0.2939
    1.4897    0.6715    1.6302    0.7269   -0.7873

计算 A 的所有元素的第 40 个和第 60 个百分位数。

P = prctile(A,[40 60],"all")

P(1) 是 A 的第 40 个百分位数，P(2) 是 A 的第 60 个百分位数。

数据矩阵的百分位数

打开实时脚本

根据指定的百分比计算数据矩阵沿列和行的百分位数。

生成一个 5×5 数据矩阵。

A = (1:5)'*(2:6)

A = 5×5

     2     3     4     5     6
     4     6     8    10    12
     6     9    12    15    18
     8    12    16    20    24
    10    15    20    25    30

计算 A 的各列的第 25 个、第 50 个和第 75 个百分位数。

P = prctile(A,[25 50 75],1)

P = 3×5

    3.5000    5.2500    7.0000    8.7500   10.5000
    6.0000    9.0000   12.0000   15.0000   18.0000
    8.5000   12.7500   17.0000   21.2500   25.5000

矩阵 P 的每列包含矩阵 A 中对应列的三个百分位数。7、12 和 17 是 A 的第三列（元素为 4、8、12、16 与 20）的第 25 个、第 50 个及第 75 个百分位数。P = prctile(A,[25 50 75]) 返回相同的结果。

计算沿 A 的行的第 25 个、第 50 个和第 75 个百分位数。

P = prctile(A,[25 50 75],2)

P = 5×3

    2.7500    4.0000    5.2500
    5.5000    8.0000   10.5000
    8.2500   12.0000   15.7500
   11.0000   16.0000   21.0000
   13.7500   20.0000   26.2500

矩阵 P 的每行包含矩阵 A 中对应行的三个百分位数。2.75、4 和 5.25 是 A 的第一行（元素为 2、3、4、5 与 6）的第 25 个、第 50 个及第 75 个百分位数。

多维数组的百分位数

打开实时脚本

计算一个多维数组沿多个维度的百分位数。

创建一个 3×5×2 数组。

A = reshape(1:30,[3 5 2])

A = 
A(:,:,1) =

     1     4     7    10    13
     2     5     8    11    14
     3     6     9    12    15


A(:,:,2) =

    16    19    22    25    28
    17    20    23    26    29
    18    21    24    27    30

通过将维度 1 和 2 指定为运算维度，计算 A 的每页的第 40 个和第 60 个百分位数。

Ppage = prctile(A,[40 60],[1 2])

Ppage = 
Ppage(:,:,1) =

    6.5000
    9.5000


Ppage(:,:,2) =

   21.5000
   24.5000

Ppage(1,1,1) 是 A 的第一页的第 40 个百分位数，Ppage(2,1,1) 是 A 的第一页的第 60 个百分位数。

通过将维度 1 和 3 指定为运算维度，计算每个 A(:,i,:) 切片中元素的第 40 个和第 60 个百分位数。

Pcol = prctile(A,[40 60],[1 3])

Pcol = 2×5

    2.9000    5.9000    8.9000   11.9000   14.9000
   16.1000   19.1000   22.1000   25.1000   28.1000

Pcol(1,4) 是 A(:,4,:) 中元素的第 40 个百分位数，Pcol(2,4) 是 A(:,4,:) 中元素的第 60 个百分位数。

给定百分比的 tall 向量的百分位数

打开实时脚本

根据给定百分比计算 tall 列向量的精确和近似百分位数。

对 tall 数组执行计算时，MATLAB® 或者使用并行池（如果有 Parallel Computing Toolbox™，则默认使用并行池），或者使用本地 MATLAB 会话。要在有 Parallel Computing Toolbox 时使用本地 MATLAB 会话运行示例，请使用 mapreducer 函数更改全局执行环境。

mapreducer(0)

为 airlinesmall 数据集创建数据存储。将 NA" 值视为缺失数据，以便 datastore 用 NaN 值替换它们。指定使用 ArrTime 变量。

ds = datastore("airlinesmall.csv","TreatAsMissing","NA", ...
    "SelectedVariableNames","ArrTime");

基于数据存储创建一个 tall 表 tt，并将数据从 tall 表提取到 tall 向量 A 中。

tt = tall(ds)

tt =

  M×1 tall table

    ArrTime
    _______

      735  
     1124  
     2218  
     1431  
      746  
     1547  
     1052  
     1134  
       :
       :

A = tt{:,:}

A =

  M×1 tall double column vector

         735
        1124
        2218
        1431
         746
        1547
        1052
        1134
         :
         :

计算 A 的第 50 个精确百分位数。由于 A 是 tall 列向量且 p 是标量，因此 prctile 默认返回精确的百分位数值。

p = 50;
Pexact = prctile(A,p)

Pexact =

  tall double

    ?

Preview deferred. Learn more.

计算 A 的第 50 个近似百分位数。指定 "approximate" 方法以使用基于 T-Digest 的逼近算法计算百分位数。

Papprox = prctile(A,p,Method="approximate")

Papprox =

  M×N×... tall array

    ?    ?    ?    ...
    ?    ?    ?    ...
    ?    ?    ?    ...
    :    :    :
    :    :    :

Preview deferred. Learn more.

使用 gather 计算 tall 数组并将结果传入内存中。

[Pexact,Papprox] = gather(Pexact,Papprox)

Evaluating tall expression using the Local MATLAB Session:
- Pass 1 of 4: Completed in 0.64 sec
- Pass 2 of 4: Completed in 0.22 sec
- Pass 3 of 4: Completed in 0.36 sec
- Pass 4 of 4: Completed in 0.27 sec
Evaluation completed in 2 sec

Pexact = 
1522

Papprox = 
1.5220e+03

精确百分位数和近似百分位数的值在所示的四位数字上相同。

tall 矩阵沿不同维度的百分位数

打开实时脚本

根据指定百分比计算 tall 矩阵沿不同维度的精确和近似百分位数。

mapreducer(0)

创建一个 tall 矩阵 A，其中包含来自 airlinesmall 数据集且存储在 varnames 中的变量子集。有关从 tall 数组提取数据的步骤的详细信息，请参阅给定百分比的 tall 向量的百分位数。

varnames = ["ArrDelay","ArrTime","DepTime","ActualElapsedTime"];
ds = datastore("airlinesmall.csv","TreatAsMissing","NA", ...
    "SelectedVariableNames",varnames);
tt = tall(ds);
A = tt{:,varnames}

A =

  M×4 tall double matrix

           8         735         642          53
           8        1124        1021          63
          21        2218        2055          83
          13        1431        1332          59
           4         746         629          77
          59        1547        1446          61
           3        1052         928          84
          11        1134         859         155
          :          :            :           :
          :          :            :           :

当沿非 1 的维度运算时，prctile 函数仅计算精确百分位数，以便它能够使用基于排序的算法（请参阅算法）高效计算，而不是使用基于 T-Digest 的逼近算法。

计算 A 沿第二个维度的第 25 个、第 50 个和第 75 个精确百分位数。

p = [25 50 75];
Pexact = prctile(A,p,2)

Pexact =

  M×N×... tall array

    ?    ?    ?    ...
    ?    ?    ?    ...
    ?    ?    ?    ...
    :    :    :
    :    :    :

Preview deferred. Learn more.

当函数沿第一个维度运算且 p 是百分比向量时，您必须使用基于 T-digest 的逼近算法计算百分位数。如果使用基于排序的算法计算 tall 数组第一个维度的百分位数，计算量会很大。

计算 A 沿第一个维度的第 25 个、第 50 个和第 75 个近似百分位数。由于默认维度是 1，因此无需为 dim 指定值。

Papprox = prctile(A,p,Method="approximate")

Papprox =

  M×N×... tall array

    ?    ?    ?    ...
    ?    ?    ?    ...
    ?    ?    ?    ...
    :    :    :
    :    :    :

Preview deferred. Learn more.

使用 gather 计算 tall 数组并将结果传入内存中。

[Pexact,Papprox] = gather(Pexact,Papprox);

Evaluating tall expression using the Local MATLAB Session:
- Pass 1 of 1: Completed in 1.4 sec
Evaluation completed in 1.9 sec

显示 A 沿第二个维度的第 25 个、第 50 个和第 75 个精确百分位数的前五行。

Pexact(1:5,:)

ans = 5×3
10³ ×

    0.0305    0.3475    0.6885
    0.0355    0.5420    1.0725
    0.0520    1.0690    2.1365
    0.0360    0.6955    1.3815
    0.0405    0.3530    0.6875

矩阵 Pexact 的每行包含 A 中对应行的三个百分位数。30.5、347.5 和 688.5 分别是 A 中第一行的第 25 个、第 50 个和第 75 个百分位数。

显示 A 沿第一个维度的第 25 个、第 50 个和第 75 个近似百分位数。

Papprox

Papprox = 3×4
10³ ×

   -0.0070    1.1149    0.9322    0.0700
         0    1.5220    1.3350    0.1020
    0.0110    1.9180    1.7400    0.1510

矩阵 Papprox 的每列包含 A 中对应列的三个百分位数。Papprox 的第一列包含 A 的第一列的百分位数。

输入参数

全部折叠

`A` — 输入数组
向量 | 矩阵 | 多维数组

输入数组，指定为向量、矩阵或多维数组。

数据类型: double | single | duration

`p` — 要计算其百分位数的百分比
标量 | 向量

要计算其百分位数的百分比，指定为 0 到 100 之间的标量或标量向量。

示例: 25

示例: [25, 50, 75]

数据类型: double | single

`dim` — 沿其运算的维度
正整数标量

沿其运算的维度，指定为正整数标量。如果不指定维度，则默认值是大小不等于 1 的第一个数组维度。

以输入矩阵 A 和由百分比 p 组成的向量为例：

P = prctile(A,p,1) 计算 A 中各列在百分比 p 的百分位数。
P = prctile(A,p,2) 计算 A 中各行在百分比 p 的百分位数。

维度 dim 指示 P 中与 p 长度相同的维度。

`vecdim` — 要沿其运算的维度组成的向量
正整数向量

要沿其运算的维度组成的向量，指定为正整数向量。每个元素代表输入数据的一个维度。

输出 P 在最小指定运算维度中的大小等于 p 的长度。P 在 vecdim 中指定的其他运算维度中的大小为 1。P 在 vecdim 中未指定的所有维度中的大小保持与输入数据相同。

假设有一个 2×3×3 输入数组 A 和百分比 p。prctile(A,p,[1 2]) 将返回一个 length(p)×1×3 数组，因为 1 和 2 是运算维度且 min([1 2]) = 1。返回的数组的每页包含 A 的对应页中元素的百分位数。

`method` — 计算百分位数的方法
`"midpoint"` (默认) | `"inclusive"` | `"exclusive"` | `"approximate"`

计算百分位数的方法，指定为以下值之一：

"midpoint" - 通过使用排序的中点算法计算百分位数。
在 R2025a 之前的版本中: 对此方法使用 "exact"。
"inclusive" - 通过使用排序的算法计算百分位数，该算法将第 0 个和第 100 个百分位数包含在数据边界内。 (自 R2025a 起)
"exclusive" - 通过使用排序的算法计算百分位数，该算法将第 0 个和第 100 个百分位数排除在数据边界外。 (自 R2025a 起)
"approximate" - 通过使用 T-Digest 的算法计算近似百分位数（适用于 double 或 single 输入数组）。

有关百分位数计算的详细信息，请参阅算法。

详细信息

全部折叠

线性插值

线性插值使用线性多项式求 y_i = f(x_i)，即底层函数 Y = f(X) 在向量或数组 x 中的点处的值。如有数据点 (x₁, y₁) 和 (x₂, y₂)，其中 y₁ = f(x₁) 且 y₂ = f(x₂)，则线性插值根据下式计算 y₁ 和 x₂ 之间给定 x 处的 x= f(x)：

$y = f (x) = y_{1} + \frac{(x - x_{1})}{(x_{2} - x_{1})} (y_{2} - y_{1}) .$

同样，如果第 100(1.5/n) 个百分位数是 y_1.5/n 且第 100(2.5/n) 个百分位数是 y_2.5/n，则线性插值根据下式计算第 100(2.3/n) 个百分位数 y_2.3/n：

$y_{\frac{2.3}{n}} = y_{\frac{1.5}{n}} + \frac{(\frac{2.3}{n} - \frac{1.5}{n})}{(\frac{2.5}{n} - \frac{1.5}{n})} (y_{\frac{2.5}{n}} - y_{\frac{1.5}{n}}) .$

T-Digest

T-digest [2] 是一种概率数据结构体，它是数据集的经验累积分布函数 (CDF) 的稀疏表示。T-digest 可用于根据在线或分布式数据计算基于秩的统计量（例如百分位数和分位数）的逼近值，所采用的方式支持可控精度，特别是在数据分布的尾部附近。

对于分布在不同分区中的数据，T-digest 会单独计算每个数据分区的分位数估计值（和百分位数估计值），然后合并这些估计值，同时保持恒定的内存界限和恒定的计算相对精度（对于第 q 个分位数，为 $q (1 - q)$ ）。由于这些原因，T-digest 在处理 tall 数组时非常实用。

要估计分布在不同分区中的数组的分位数，首先在数据的每个分区中构建一个 T-digest。T-digest 对分区中的数据进行聚类，并通过质心值和表示对聚类有贡献的采样数量的累积权重来汇总每个聚类。T-digest 使用大簇（质心间距大）表示 q = 0.5 附近的 CDF 区域，使用小簇（质心间距小）表示 q = 0 和 q = 1 附近的 CDF 区域。

T-digest 通过使用缩放函数控制簇大小，该函数使用压缩参数 δ 将分位数 q 映射到索引 k。即

$k (q, δ) = δ \cdot (\frac{\sin^{- 1} (2 q - 1)}{π} + \frac{1}{2}),$

，其中映射 k 呈单调形态，最小值为 k(0,δ) = 0，最大值为 k(1,δ) = δ。下图显示了 δ = 10 的缩放函数。

Scaling function plotted with q on the x-axis and k on the y-axis

缩放函数将分位数 q 转换为缩放因子 k，以便在 q 中给出可变大小的步。因此，簇大小不相等（中心分位数附近的簇较大，q = 0 和 q = 1 附近的簇较小）。较小的簇使得数据边缘附近的准确度更高。

要用具有权重和位置的新观测值更新 T-digest，请找到距离新观测值最近的簇。然后，添加权重，并基于加权平均值更新簇质心，前提是更新后的簇权重不超过大小限制。

您可以合并来自每个数据分区的独立 T-digest，方法是对这些 T-digest 求并集，并合并其质心。要合并 T-digest，请首先对来自所有独立 T-digest 的簇按簇权重降序排序。然后，合并满足大小限制的相邻簇，以形成一个新 T-digest。

一旦形成表示完整数据集的 T-digest，就可以估计 T-digest 中每个簇的端点（或边界），然后使用每个簇的端点之间的插值来计算准确的分位数估计值。

算法

对于 n 元素向量 A，当您选择除 "approximate" 之外的任何方法时，prctile 函数使用基于排序的算法计算百分位数。

A 中的排序元素根据选定方法映射到百分位数，如下表中所述。

百分位数	`Method`
百分位数	`"midpoint"` 在 R2025a 之前的版本中: `"exact"`	`"inclusive"` (自 R2025a 起)	`"exclusive"` (自 R2025a 起)
第一个排序元素的百分位数	50/n	0	100/(n+1)
第二个排序元素的百分位数	150/n	100/(n−1)	200/(n+1)
第三个排序元素的百分位数	250/n	200/(n−1)	300/(n+1)
...	...	...	...
第 k 个排序元素的百分位数	50(2k−1)/n	100(k−1)/(n−1)	100k/(n+1)
...	...	...	...
第 (n−1) 个排序元素的百分位数	50(2n−3)/n	100(n−2)/(n−1)	100(n−1)/(n+1)
第 n 个排序元素的百分位数	50(2n−1)/n	100	100n/(n+1)

例如，如果 A 是 [6 3 2 10 1]，则百分位数如下表所示。

百分位数	`Method`
百分位数	`"midpoint"` 在 R2025a 之前的版本中: `"exact"`	`"inclusive"` (自 R2025a 起)	`"exclusive"` (自 R2025a 起)
`1` 的百分位数	10	0	50/3
`2` 的百分位数	30	25	100/3
`3` 的百分位数	50	50	50
`6` 的百分位数	70	75	200/3
`10` 的百分位数	90	100	250/3

prctile 函数使用线性插值计算 A 的第一个和最后一个排序元素之间的百分比的百分位数。有关详细信息，请参阅线性插值。
例如，如果 A 是 [6 3 2 10 1]，则：
- 对于中点方法，第 40 个百分位数是 2.5。
  在 R2025a 之前的版本中: 对于精确方法，第 40 个百分位数是 2.5。
- 对于包含方法，第 40 个百分位数是 2.6。 (自 R2025a 起)
- 对于排除方法，第 40 个百分位数是 2.4。 (自 R2025a 起)
prctile 函数将 A 中元素的最小值或最大值赋给该范围之外的百分比对应的百分位数。
例如，如果 A 是 [6 3 2 10 1]，则对于中点方法和排除方法，第 5 个百分位数是 1。 (自 R2025a 起)
在 R2025a 之前的版本中: 例如，如果 A 是 [6 3 2 10 1]，则对于精确方法，第 5 个百分位数是 1。

prctile 函数将 NaN 值视为缺失值并将其删除。

参考

[1] Langford, E. “Quartiles in Elementary Statistics”, Journal of Statistics Education. Vol. 14, No. 3, 2006.

[2] Dunning, T., and O. Ertl. “Computing Extremely Accurate Quantiles Using T-Digests.” August 2017.

扩展功能

全部展开

tall 数组
对行数太多而无法放入内存的数组进行计算。

prctile 函数支持 tall 数组，但存在以下使用说明和限制：

仅当 A 是 tall 数值列向量时，P = prctile(A,p) 才返回精确百分位数（使用基于排序的算法）。
仅当满足以下条件之一时，P = prctile(A,p,dim) 才返回精确百分位数：
- A 是 tall 数值列向量。
- A 是 tall 数值数组且 dim 不是 1。例如，prctile(A,p,2) 返回 tall 数组 A 沿行的精确百分位数。
如果 A 是 tall 数值数组且 dim 是 1，则您必须将 method 指定为 "approximate"，以使用基于 T-Digest 的逼近算法计算百分位数。例如，prctile(A,p,1,"Method","approximate") 返回 tall 数组 A 沿列的近似百分位数。
仅当满足以下条件之一时，P = prctile(A,p,vecdim) 才返回精确百分位数：
- A 是 tall 数值列向量。
- A 是 tall 数值数组且 vecdim 不包含 1。例如，如果 A 是 3×5×2 数组，则 prctile(A,p,[2,3]) 返回每个 A(i,:,:) 切片中元素的精确百分位数。
- A 是 tall 数值数组且 vecdim 包含 1 以及 A 中其大小不等于 1 的所有维度。例如，如果 A 是 10×1×4 数组，则 prctile(A,p,[1 3]) 返回 A(:,1,:) 中元素的精确百分位数。
如果 A 是 tall 数值数组且 vecdim 包含 1 但不包含 A 中其大小不等于 1 的所有维度，则必须将 method 指定为 "approximate" 以使用逼近算法。例如，如果 A 是 10×1×4 数组，您可以使用 prctile(A,p,[1 2],"Method","approximate") 计算 A 的每页的近似百分位数。

有关详细信息，请参阅tall 数组。

C/C++ 代码生成
使用 MATLAB® Coder™ 生成 C 代码和 C++ 代码。

用法说明和限制：

不支持 "all" 和 vecdim 输入。
不支持 Method 名称-值参量。
dim 输入参量必须为编译时常量。
如果不指定 dim 输入参量，生成的代码中的工作（或运算）维度可能不同。因此，可能会发生运行时错误。有关详细信息，请参阅在默认维度选择方面与 MATLAB 不兼容 (MATLAB Coder)。
如果输出 P 是向量，当以下所有条件都满足时，P 的方向与 MATLAB^® 不同：
- 未提供 dim。
- A 在编译时是可变大小数组，而不是可变大小向量，但在运行时 A 是向量。
- 向量 A 的方向与向量 p 的方向不匹配。
在这种情况下，输出 P 与 A 的方向匹配，而不是与 p 的方向匹配。

GPU 数组
通过使用 Parallel Computing Toolbox™ 在图形处理单元 (GPU) 上运行来加快代码执行。

prctile 函数支持 GPU 数组输入，但有以下用法说明和限制：

不支持 "all" 和 vecdim 输入。
不支持 Method 名称-值参量。

有关详细信息，请参阅在 GPU 上运行 MATLAB 函数 (Parallel Computing Toolbox)。

分布式数组
使用 Parallel Computing Toolbox™ 在集群的组合内存中对大型数组进行分区。

用法说明和限制：

不支持持续时间输入。

有关详细信息，请参阅使用分布式数组运行 MATLAB 函数 (Parallel Computing Toolbox)。

版本历史记录

在 R2006a 之前推出

全部展开

R2025a: 使用包含法和排除法计算百分位数

您可以使用包含法或排除法计算百分位数。将方法名称分别指定为 "inclusive" 或 "exclusive"。"inclusive" 方法在数据边界内包含第 0 个和第 100 个百分位数；"exclusive" 方法则排除它们。

此外，默认方法的名称已从 "exact" 更改为 "midpoint"。prctile 函数继续支持 Method="exact" 以实现向后兼容。

以下示例说明默认方法与两种新方法之间的差异。

x = 1:5;
P = prctile(x,25,Method="midpoint")

P = 1.7500

x = 1:5;
P = prctile(x,25,Method="inclusive")

P = 2

x = 1:5;
P = prctile(x,25,Method="exclusive")

P = 1.5000

R2022b: 改进了使用较少的输入数据时的性能

prctile 函数由于输入解析更快而提高了性能。当输入解析在计算时间中占比很大时，性能改进最为显著。这种情况出现在：

输入数据的大小很小时。
要为其计算百分位数的百分比的数量很少。
计算沿默认运算维度进行时。

例如，以下代码为一个包含 3000 个元素的矩阵计算四个百分位数。该代码的执行速度比上一版本大约快 5 倍。

function timingPrctile
A = rand(300,10);
for k = 1:3e3
  P = prctile(A,[20 40 60 80]);
end
end

大致的执行时间是：

R2022a：1.0 秒

R2022b：0.2 秒

代码是在运行 Windows^® 10 的 Intel^® Xeon^® CPU E5-1650 v4 @ 3.60 GHz 测试系统上使用 timeit 函数进行计时的：

timeit(@timingPrctile)

R2022a: 从 Statistics and Machine Learning Toolbox 移至 MATLAB

以前，prctile 需要 Statistics and Machine Learning Toolbox™。

另请参阅

quantile | median | iqr

prctile

语法

说明

示例

数据向量的百分位数

所有值的百分位数

数据矩阵的百分位数

多维数组的百分位数

给定百分比的 tall 向量的百分位数

tall 矩阵沿不同维度的百分位数

输入参数

A — 输入数组 向量 | 矩阵 | 多维数组

p — 要计算其百分位数的百分比 标量 | 向量

dim — 沿其运算的维度 正整数标量

vecdim — 要沿其运算的维度组成的向量 正整数向量

method — 计算百分位数的方法 "midpoint" (默认) | "inclusive" | "exclusive" | "approximate"

详细信息

线性插值

T-Digest

算法

参考

扩展功能

tall 数组 对行数太多而无法放入内存的数组进行计算。

C/C++ 代码生成 使用 MATLAB® Coder™ 生成 C 代码和 C++ 代码。

GPU 数组 通过使用 Parallel Computing Toolbox™ 在图形处理单元 (GPU) 上运行来加快代码执行。

分布式数组 使用 Parallel Computing Toolbox™ 在集群的组合内存中对大型数组进行分区。

版本历史记录

R2025a: 使用包含法和排除法计算百分位数

R2022b: 改进了使用较少的输入数据时的性能

R2022a: 从 Statistics and Machine Learning Toolbox 移至 MATLAB

另请参阅

`A` — 输入数组
向量 | 矩阵 | 多维数组

`p` — 要计算其百分位数的百分比
标量 | 向量

`dim` — 沿其运算的维度
正整数标量

`vecdim` — 要沿其运算的维度组成的向量
正整数向量

`method` — 计算百分位数的方法
`"midpoint"` (默认) | `"inclusive"` | `"exclusive"` | `"approximate"`

tall 数组
对行数太多而无法放入内存的数组进行计算。

C/C++ 代码生成
使用 MATLAB® Coder™ 生成 C 代码和 C++ 代码。

GPU 数组
通过使用 Parallel Computing Toolbox™ 在图形处理单元 (GPU) 上运行来加快代码执行。

分布式数组
使用 Parallel Computing Toolbox™ 在集群的组合内存中对大型数组进行分区。