LOF异常值剔除算法

%LOF算法
function lof =  LOF（dist）
K=10;
m=size（dist1）;                 %m为对象数，dist为两两之间的距离
distance = zeros（mm）;
num = zeros（mm）;               %distance 和num用来记录排序后的顺序，和对象编号顺序
kdistance =  zeros（m1）;         %计算每个对象的kdistance
count  = zeros（m1）;            %k邻域的对象数
reachdist =  zeros（mm）;         %计算两两之间的reachable-distance
lrd = zeros（m1）;
lof = zeros（m1）;
%计算k-距离
for i=1:m
     [distance（i:）num（i:）]=sort（dist（i:）‘ascend‘）;
    kdistance（i）=distance（iK+1）;
    count（i） = -1;%自己的距离为0，要去掉自己
    for j = 1:m
        if dist（i