ufc_prold.m

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

function [U,V,J,ll,r]=ufc_pr(X,U,c,r,m,g,pnorm,cluster_error)

[leng,dim] = size(X);
E = zeros(dim,dim,c);
S = zeros(dim,r,c);
Dr = zeros(c,leng);
Dp = zeros(c,leng);
D = zeros(c,leng);
while 1
    a = U * (X.^m);
    b = sum(U')';
    b = repmat(b,1,dim);
    V = a ./ b;
        
    for i = 1:c
        a = X - repmat(V(i,:),leng,1);
        b = a .* repmat(U(i,:)',1,dim);
        E(:,:,i) = b' * a;
        
        pc = pcacov(E(:,:,i));
        S(:,:,i) = pc(:,1:r);
        a = X - repmat(V(i,:), leng, 1);
        Dr(i,:) =  sum(S(:,:,i)' * a')';
        
        Dp(i,:) = sum(a'.^pnorm);
    end
    
    D = Dp + g * Dr;
    U_new = D ./ repmat(sum(D), c, 1);
    U_error = norm(U_new - U)
    if U_error < cluster_error
        break;
    end
    U = U_new;
end