求多元数据的经验分布函数的MATLAB代码

下面是求解多元数据的经验分布函数的MATLAB代码，其中X为n行p列的数据矩阵，表示n个p维随机向量：

function [F, x] = ecdf_multivariate(X)
% 输入：X为n行p列的数据矩阵，表示n个p维随机向量
% 输出：F为多元经验分布函数，x为对应的数据点

[n, p] = size(X); % n个p维随机向量
x = unique(X, 'rows'); % 去重得到数据点
m = size(x, 1); % 数据点个数

F = zeros(m, 1); % 初始化多元经验分布函数
for i = 1:m
    F(i) = sum(all(X<=repmat(x(i,:), n, 1), 2))/n;
end
end

使用方法：将上述函数保存为.m文件，然后在MATLAB命令行中输入以下代码：

X = randn(100,3); % 生成100个3维标准正态分布的随机向量
[F, x] = ecdf_multivariate(X); % 计算多元经验分布函数
plot(x, F, 'o'); % 绘制多元经验分布函数图像

上述代码会生成100个3维标准正态分布的随机向量，并计算它们的多元经验分布函数，最后绘制出图像。