function?idx?=?AP（S）
N?=?size（S1）;
%初始化為0
A=zeros（NN）;%Availabilities
R=zeros（NN）;?%Responsibility
?
lam=0.9;?%?Set?damping?factor
%前后兩次聚類結(jié)果相同的累計(jì)次數(shù)，初始為-1
same_time?=?-1;
%最多進(jìn)行10000次循環(huán)退出
for?iter=1:10000
????%?Compute?responsibilities
????Rold=R;
????AS=A+S;
????%找到每行中最大的值max（a（ik）+s（ik））Y是一列最大值，I是序號(hào)
????[YI]=max（AS[]2）;
????%每行的最大值設(shè)為-1000，不會(huì)影響第二大值的尋找
????for?i=1:N
????????AS（iI（i））=-1000;
????end
????%找到每行中第二大的值Y2，計(jì)算最大值位置的r（ik）用
????[Y2~]=max（AS[]2）;
????R=S-repmat（Y[1N]）;
????for?i=1:N
????????R（iI（i））=S（iI（i））-Y2（i）;
????end
????R=（1-lam）*R+lam*Rold;?%?Dampen?responsibilities
????%?Compute?availabilities
????Aold=A;
????Rp=max（R0）;
????for?k=1:N
????????Rp（kk）=R（kk）;
????end
????%豎著復(fù)制
????A=repmat（sum（Rp1）[N1]）-Rp;
????dA=diag（A）;
????A=min（A0）;
????for?k=1:N
????????A（kk）=dA（k）;
????end;
????A=（1-lam）*A+lam*Aold;?%?Dampen?availabilities
?
????if（same_time?==?-1）%第一次
????????E=R+A;
????????[~?idx_old]?=?max（E[]2）;
????????same_time?=?0;
????else
????????E=R+A;
????????[~?idx]?=?max（E[]2）;%計(jì)算得到新的idx
?
????????if（sum（abs（idx_old-idx））?==?0）%前后兩次迭代得到的分類結(jié)果完全相同
????????????same_time?=?same_time?+?1;%次數(shù)加一
????????????if（same_time?==?10）%分類結(jié)果穩(wěn)定在當(dāng)前結(jié)果10次
????????????????%輸出循環(huán)次數(shù)
????????????????disp（num2str（iter））;
????????????????break;
????????????end
????????end
????????%尚未完全相同，繼續(xù)迭代
????????idx_old?=?idx;
????end
end

E=R+A;
[~?idx]?=?max（E[]2）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-06-03?23:18??AP\
?????文件????????3487??2016-05-05?16:02??AP\Aggregation.mat
?????文件????????1621??2016-06-02?19:59??AP\AP.m
?????文件????????1344??2016-06-03?01:48??AP\AP_3D.m
?????文件?????????952??2016-06-02?21:58??AP\AP_top.m
?????文件????????1854??2016-05-05?16:02??AP\Compound.mat
?????文件???????32149??2016-05-05?16:02??AP\D31.mat
?????文件???????17625??2016-04-18?21:06??AP\depth.png
?????文件????????1226??2016-05-05?16:02??AP\flame.mat
?????文件????????1810??2016-05-05?16:02??AP\jain.mat
?????文件????????1499??2016-05-05?16:02??AP\pathbased.mat
?????文件????????4373??2016-05-05?16:02??AP\R15.mat
?????文件????????1556??2016-05-05?16:02??AP\spiral.mat