function?LOF
clc
clear
fileOrder=1;
%讀取原始數(shù)據(jù)
x=readFile（fileOrder）;
K=10;
%對(duì)數(shù)據(jù)操作
count=1;%outlier點(diǎn)個(gè)數(shù)
x=opeData（xcountfileOrder）;
[mn]=size（x）;
%%?計(jì)算各個(gè)點(diǎn)之間距離
for?i=1:m
????for?j=i:m
????????dist（ij）=abs（x（i）-x（j））;
????????dist（ji）=dist（ij）;
????end
end

%%?核心部分，對(duì)每個(gè)點(diǎn)分別計(jì)算LOF值
for?i=1:m
????%計(jì)算k-distance
????[distancenum]=sort（dist（i:）‘a(chǎn)scend‘）;
????kdistance=distance（K+1）;
????%集合
????pointArray（i:）=num（2:K+1）;
????%???對(duì)集合中每個(gè)點(diǎn)計(jì)算可達(dá)距離
????n=size（pointArray2）;
????for?j=1:n
????????d（ij）=distance（j+1）;
????????[distempnumtemp]=sort（dist（i:）‘a(chǎn)scend‘）;
????????kdistemp=distemp（K+1）;
????????kdis（j）=kdistemp;
????????%可達(dá)距離
????????reachdis（ij）=max（d（ij）kdis（j））;
????end
????sum_reachdis=0;
????for?j=1:n
????????sum_reachdis=sum_reachdis+reachdis（ij）;
????end
????%計(jì)算每個(gè)點(diǎn)的lrd
????lrd（i）=n/sum_reachdis;
end
%%?得到lof值
for?i=1:m
????sumlrd=0;
????for?j=1:n
????????sumlrd=sumlrd+lrd（j）/lrd（i）;
????end
????lof（i）=sumlrd/n;
end
%%畫(huà)圖
subplot（211）
plot（1:size（x1）x‘-b.‘‘linewidth‘2‘markersize‘14.5）;
legend（‘原始數(shù)據(jù)‘）;
title（‘LOF方法效果圖‘）;
set（legend‘location‘‘best‘）;
subplot（212）
plot（1:size（lof2