fuzzydbindex.m

来自「ufc无监督优化模糊聚类用于彩色图像分割」· M 代码 · 共 24 行

function [FR,DB] = FuzzyDBIndex(X,U,V,c,m)

[leng,dim] = size(X);
dp = zeros(1,c);
dv = zeros(c,c);
FR = zeros(c,c);

P = (U == repmat(max(U),c,1));
for i = 1:c
    Pi = find(P(i,:));
    n(i) = length(Pi);
    Vi = repmat(V(i,:),leng,1);
    dp(i) = sqrt(U(i,:).^m * sum((X-Vi)'.^2)' / n(i));
    %dp(i) = sqrt(sum(sum((X(Pi,:)-Vi)'.^2)') / n(i));
    Vi = Vi(1:c,:);
    dv(i,:) = sqrt(sum((Vi-V)'.^2)) + realmin;
end

for i = 1:c
    dpi = repmat(dp(i),1,c);   
    FR(i,:) = (dpi + dp) ./ dv(i,:);
end
FR = FR - diag(diag(FR));
DB = mean(max(FR));