%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-distance
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（ij）<=kdistance（i）
????????????count（i）?=?count（i）+1;
????????end
????end
end
????
for?i?=?1:m
????for?j=1:i-1
????????reachdist（ij）?=?max（dist（ij）kdistance（j））;
????????reachdist（ji）?=?reachdist（ij）;
????end
end

for?i?=?1:m
????sum_reachdist=0;
????for?j=1:count（i）
????????sum_reachdist=sum_reachdist+reachdist（inum（j+1））;
????end
????%計算每個點的lrd
????lrd（i）=count（i）/sum_reachdist;
end
%?得到lof值
for?i=1:m
????sumlrd=0;
????for?j=1:count（i）
????????sumlrd=sumlrd+lrd（num（j+1））/lrd（i）;
????end
????lof（i）=sumlrd/count（i）;
end

%畫圖
%?subplot（211）
%?plot（1:size（x1）x‘-b.‘‘linewidth‘2‘markersize‘14.5）;
%?legend（‘原始數據‘）;
%?title（‘LOF方法效果圖‘）;
%?set（legend‘location‘‘best‘）;
%?subplot（212）
%?plot（1:size（lof2）lof‘-r.‘‘linewidth‘2‘markersize‘14.5）;
%?legend（‘LOF值‘）;
%?set（legend‘location‘‘best‘）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1519??2017-12-28?15:40??LOF.m
?????文件?????????455??2017-12-28?15:40??testLOF.m