%LOF算法
%dist為m*m的距離矩陣每一行代表一個(gè)數(shù)據(jù)與其它數(shù)據(jù)距離的列向量，所以該矩陣為
%對(duì)角線為0的，關(guān)于對(duì)角線對(duì)稱的矩陣K為k-近鄰
%?function?lof?=?LOF（distK）

clear;
clc;

A=importdata（‘data.mat‘）;%需要進(jìn)行離群因子算法處理的數(shù)據(jù)集
numData=size（A1）;
KD=[];
for?i=1:1:numData
[~dist]=knnsearch（A（i:）A（::））;
KD=[KD;dist‘];
end

m=size（dist1）;?????????????????%m為對(duì)象數(shù)，dist為兩兩之間的距離
distance?=?zeros（mm）;
num?=?zeros（mm）;???????????????%distance?和num用來記錄排序后的順序，和對(duì)象編號(hào)順序
kdistance?=?zeros（m1）;?????????%計(jì)算每個(gè)對(duì)象的kdistance
count??=?zeros（m1）;????????????%k鄰域的對(duì)象數(shù)
reachdist?=?zeros（mm）;?????????%計(jì)算兩兩之間的reachable-distance
lrd?=?zeros（m1）;
lof?=?zeros（m1）;
%計(jì)算k-distance
for?i=1:m?
????[distance（i:）num（i:）]=sort（dist（i:）‘a(chǎn)scend‘）;%distance按照升序?qū)ist進(jìn)行排序，num記錄排序前各個(gè)dist所在位置信息
????kdistance（i）=distance（iK+1）;%獲得k近鄰距離，因?yàn)榕判蚝蟮谝粋€(gè)值為自身到自身的距離為0，所以k+1才是第k個(gè)最近距離
????count（i）?=?-1;???????????????%自己的距離為0，要去掉自己
????for?j?=?1:m????????????