M = 0; S=0;
for i = 1:Nb_Points
S = 0;
for j=1:K
S = S + 1/((sum((X(i,:)-Centres(j,:)).^2).^0.5))^q;
end
M = M + K/S;
end
T1 = zeros(Nb_Points, K);
for i = 1:Nb_Points
for j = 1:K
if (sum((X(i,:)-Centres(j,:)).^2) == 0)
T1(i,j) = Epsilon;
else
T1(i,j) = (((sum((X(i,:)-Centres(j,:)).^2).^0.5))^(-q-2));
end
end
end
for i = 1:Nb_Points
T2(i) = 0;
for j = 1:K
if (sum((X(i,:)-Centres(j,:)).^2) == 0)
T2(i) = Epsilon;
else
T2(i) = T2(i) + (((sum((X(i,:)-Centres(j,:)).^2).^0.5))^(-q-2));
end
end
end
for j=1:K
T(:,j) = T1(:,j)./T2';
end
L = zeros(Nb_Points,1);
for i = 1:Nb_Points
Denomi = 0; Nomi = 0;
for j = 1:K
if (sum((X(i,:)-Centres(j,:)).^2) == 0)
Denomi = Epsilon;
Nomi = Epsilon;
else
Denomi = Denomi + (((sum((X(i,:)-Centres(j,:)).^2).^0.5))^(-q))^2;
Nomi = Nomi + ((sum((X(i,:)-Centres(j,:)).^2).^0.5))^(-q-2);
end
end
L(i) = Nomi./Denomi;
end