Main Content

本页翻译不是最新的。点击此处可查看最新英文版本。

pdf

高斯混合分布的概率密度函数

说明

示例

y = pdf(gm,X) 返回高斯混合分布 gm 的概率密度函数 (pdf),在 X 中的值处计算函数值。

示例

全部折叠

创建一个 gmdistribution 对象并计算其 pdf 值。

定义一个二成分二元高斯混合分布的分布参数(均值和协方差)。

mu = [1 2;-3 -5];
sigma = [1 1]; % shared diagonal covariance matrix

使用 gmdistribution 函数创建一个 gmdistribution 对象。默认情况下,该函数会创建一个等比例混合。

gm = gmdistribution(mu,sigma)
gm = 

Gaussian mixture distribution with 2 components in 2 dimensions
Component 1:
Mixing proportion: 0.500000
Mean:     1     2

Component 2:
Mixing proportion: 0.500000
Mean:    -3    -5

计算 gm 的 pdf 值。

X = [0 0;1 2;3 3;5 3];
pdf(gm,X)
ans = 4×1

    0.0065
    0.0796
    0.0065
    0.0000

创建一个 gmdistribution 对象并对绘制其 pdf。

定义两个二元高斯混合成分的分布参数(均值、协方差和混合比例)。

p = [0.4 0.6];               % Mixing proportions     
mu = [1 2;-3 -5];            % Means
sigma = cat(3,[2 .5],[1 1])  % Covariances 1-by-2-by-2 array
sigma = 
sigma(:,:,1) =

    2.0000    0.5000


sigma(:,:,2) =

     1     1

cat 函数沿第三个数组维度串联协方差。定义的协方差矩阵是对角矩阵。sigma(1,:,i) 包含成分 i 的协方差矩阵的对角线元素。

使用 gmdistribution 函数创建一个 gmdistribution 对象。

gm = gmdistribution(mu,sigma)
gm = 

Gaussian mixture distribution with 2 components in 2 dimensions
Component 1:
Mixing proportion: 0.500000
Mean:     1     2

Component 2:
Mixing proportion: 0.500000
Mean:    -3    -5

使用 fsurf 绘制高斯混合分布的 pdf。

gmPDF = @(x,y) arrayfun(@(x0,y0) pdf(gm,[x0 y0]),x,y);
fsurf(gmPDF,[-10 10])

Figure contains an axes object. The axes object contains an object of type functionsurface.

输入参数

全部折叠

高斯混合分布,也称为高斯混合模型 (GMM),指定为 gmdistribution 对象。

您可以使用 gmdistributionfitgmdist 创建一个 gmdistribution 对象。通过指定分布参数,使用 gmdistribution 函数创建一个 gmdistribution 对象。在给定固定成分数量的情况下,使用 fitgmdist 函数对数据进行 gmdistribution 模型拟合。

用于计算 pdf 的值,指定为 n×m 数值矩阵,其中 n 是观测值数目,m 是每个观测值中的变量数。

数据类型: single | double

输出参数

全部折叠

高斯混合分布 gm 的 pdf 值,在 X 处计算,以 n×1 数值向量形式返回,其中 n 是 X 中的观测值数目。

在给定每个观测值和成分概率的情况下,pdf 函数通过使用每个成分的似然来计算 pdf 值。

y(i)=j=1kL(Cj|Oi)P(Cj),

其中 L(Cj|Oj) 是给定观测值 i 时成分 j 的似然,P(Cj) 是成分 j 的概率。pdf 函数通过使用在观测值 i 处计算的第 j 个高斯混合成分的多元正态 pdf 来计算似然项。成分概率是混合成分的混合比例,即 gmComponentProportion 属性。

版本历史记录

在 R2007b 中推出