function[R]=F_JlR(cs,X)%定義函數
%模糊聚類分析建立模糊相似矩陣: [R]=F_JlR(cs,X)
%X,數據矩陣
%cs=1,數量積法
%cs=2,夾角余弦法
%cs=3,相關系數法
%cs=4,指數相似系數法
%cs=5,最大最小法
%cs=6,算術平均最小法
%cs=7,幾何平均最小法
%cs=8,一般歐式距離法
%cs=9,一般海明距離法
%cs=10,一般切比雪夫距離法
%cs=11,倒數歐式距離法
%cs=12,倒數海明距離法
%cs=13,倒數切比雪夫距離法
%cs=14,指數歐式距離法
%cs=15,指數海明距離法
%cs=16,指數切比雪夫距離法
[n,m]=size(X);%獲得矩陣的行列數
R=[];
if(cs==1)maxM=0;pd=0;%數量積法
	for(i=1:n)for(j=1:n)if(j~=i)x=0;
		for(k=1:m)x=x+X(i,k)*X(j,k);end
		if(maxM<x)maxM=x;end
    end;end;end
	if(maxM<0.000001)return;end
	for(i=1:n)for(j=1:n)
		if(i==j)R(i,j)=1;
		else R(i,j)=0;
			for(k=1:m)R(i,j)=R(i,j)+X(i,k)*X(j,k);end
			R(i,j)=R(i,j)/maxM;
			if(R(i,j)<0)pd=1;end
		end
    end;end
	if(pd)for(i=1:n)for(j=1:n)R(i,j)=(R(i,j)+1)/2;end;end;end
elseif(cs==2)%夾角余弦法
	for(i=1:n)for(j=1:n)xi=0;xj=0;
		for(k=1:m)xi=xi+X(i,k)^2;xj=xj+X(j,k)^2;end
		s=sqrt(xi*xj);R(i,j)=0;
		for(k=1:m)R(i,j)=R(i,j)+X(i,k)*X(j,k);end
		R(i,j)=R(i,j)/s;
    end;end
elseif(cs==3)%相關系數法
	for(i=1:n)for(j=1:n)xi=0;xj=0;
		for(k=1:m)xi=xi+X(i,k);xj=xj+X(j,k);end
		xi=xi/m;xj=xj/m;xis=0;xjs=0;
		for(k=1:m)xis=xis+(X(i,k)-xi)^2;xjs=xjs+(X(j,k)-xj)^2;end
		s=sqrt(xis*xjs);R(i,j)=0;
		for(k=1:m)R(i,j)=R(i,j)+abs((X(i,k)-xi)*(X(j,k)-xj));end
		R(i,j)=R(i,j)/s;
    end;end
elseif(cs==4)%指數相似系數法
	for(i=1:n)for(j=1:n)R(i,j)=0;
		for(k=1:m)xk=0;
			for(z=1:n)xk=xk+X(z,k);end
			xk=xk/n;sk=0;
			for(z=1:n)sk=sk+(X(z,k)-xk)^2;end
			sk=sk/n;R(i,j)=R(i,j)+exp(-0.75*((X(i,k)-X(j,k))/sk)^2);
        end
		R(i,j)=R(i,j)/m;
    end;end
elseif(cs<=7)%最大最小法 算術平均最小法 幾何平均最小法
	for(i=1:n)for(j=1:n)fz=0;fm=0;
		for(k=1:m)
			if(X(j,k)<0)R=[];return;end
			if(X(j,k)>X(i,k))x=X(i,k);
			else x=X(j,k);end
			fz=fz+x;
        end
		if(cs==5)%最大最小法
			for(k=1:m)if(X(i,k)>X(j,k))x=X(i,k);else x=X(j,k);end
			fm=fm+x;end
		elseif(cs==6)for(k=1:m)fm=fm+(X(i,k)+X(j,k))/2;end%算術平均最小法
		else for(k=1:m)fm=fm+sqrt(X(i,k)*X(j,k));end;end%幾何平均最小法
		R(i,j)=fz/fm;
    end;end
elseif(cs<=10)C=0;%一般距離法
	for(i=1:n)for(j=i+1:n)d=0;
		if(cs==8)for(k=1:m)d=d+(X(i,k)-X(j,k))^2;end
			d=sqrt(d);%歐式距離
		elseif(cs==9)for(k=1:m)d=d+abs(X(i,k)-X(j,k));end%海明距離
		else for(k=1:m)if(d<abs(X(i,k)-X(j,k)))d=abs(X(i,k)-X(j,k));end;end;end%切比雪夫距離
		if(C<d)C=d;end
	end;end
	C=1/(1+C);
	for(i=1:n)for(j=1:n)d=0;
		if(cs==8)for(k=1:m)d=d+(X(i,k)-X(j,k))^2;end
			d=sqrt(d);%歐式距離
		elseif(cs==9)for(k=1:m)d=d+abs(X(i,k)-X(j,k));end%海明距離
		else for(k=1:m)if(d<abs(X(i,k)-X(j,k)))d=abs(X(i,k)-X(j,k));end;end;end%切比雪夫距離
		R(i,j)=1-C*d;
	end;end
elseif(cs<=13)minM=Inf;%倒數距離法
	for(i=1:n)for(j=i+1:n)d=0;
		if(cs==11)for(k=1:m)d=d+(X(i,k)-X(j,k))^2;end
			d=sqrt(d);%歐式距離
		elseif(cs==12)for(k=1:m)d=d+abs(X(i,k)-X(j,k));end%海明距離
		else for(k=1:m)if(d<abs(X(i,k)-X(j,k)))d=abs(X(i,k)-X(j,k));end;end;end%切比雪夫距離
		if(minM>d)minM=d;end
	end;end
	minM=0.9999*minM;
	if(minM<0.000001)return;end
	for(i=1:n)for(j=1:n)d=0;
		if(j==i)R(i,j)=1;continue;end
		if(cs==11)for(k=1:m)d=d+(X(i,k)-X(j,k))^2;end
			d=sqrt(d);%歐式距離
		elseif(cs==12)for(k=1:m)d=d+abs(X(i,k)-X(j,k));end%海明距離
		else for(k=1:m)if(d<abs(X(i,k)-X(j,k)))d=abs(X(i,k)-X(j,k));end;end;end%切比雪夫距離
		R(i,j)=minM/d;
	end;end
else for(i=1:n)for(j=1:n)d=0;%指數距離法
	if(cs==14)for(k=1:m)d=d+(X(i,k)-X(j,k))^2;end
		d=sqrt(d);%歐式距離
	elseif(cs==15)for(k=1:m)d=d+abs(X(i,k)-X(j,k));end%海明距離
	else for(k=1:m)if(d<abs(X(i,k)-X(j,k)))d=abs(X(i,k)-X(j,k));end;end;end%切比雪夫距離
	R(i,j)=exp(-d);
end;end;end