给出lof的matlab代码

以下是一个简单的 MATLAB 实现 Local Outlier Factor (LOF) 算法的代码：

function lof = LOF(X, k)
% X: 数据集，每行代表一个样本
% k: 选取的邻居个数

[n, d] = size(X);
dist = pdist2(X, X); % 计算样本间距离
lof = zeros(n, 1);

for i = 1:n
    neighbors = knnsearch(X, X(i, :), 'K', k+1); % 找到第 i 个样本的 k+1 个邻居
    neighbors = neighbors(2:end); % 去除自身
    reach_dist = max(dist(neighbors, i), [], 1); % 与 k 邻居的可达距离
    lrd = k / sum(dist(neighbors, neighbors), 1); % 局部可达密度
    lof(i) = sum(lrd) / (k * lrd(i)) * reach_dist; % LOF
end

需要注意的是，这段代码中使用了 MATLAB 自带的 `pdist2` 函数计算样本间的距离，并且使用了 MATLAB 自带的 `knnsearch` 函数找到每个样本的 k 个邻居。因此，这段代码可能不是最优化的实现方式，但是可以用于简单的 LOF 算法实验。