function[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
  maxM=maxM+1;
  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;end
        fm=fm+x;
      elseif(js==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=lnf;%％倒數距離法
    for(i=1:n)for(j=i+1:n)d=0;
    if(cs==11)for(k=1:m)d=d+(X(j,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