1、程序 1：求傳遞閉包的程序:unction t=bibao(X,c) I=ones(size(X,1);Y=zeros(size(X,1); for i=1:size(X,1)for j=1:size(X,1)for k=1:size(X,2)Y(i,j)=abs(X(i,k)-X(j,k)+Y(i,j); endendenddisp(模糊相似矩陣為：)R=I-c*Yflag=0; C=R;B=R;while flag=0for i=1:size(X,1)for j=1:size(X,1)for k=1:size(X,1)B(i,j)=max(min(C(i,k),C(k,j),B(i,j);

2、endendendif B=Cflag=1;elseC=B;endend t=B;程序 2： R 的K 截矩陣 function Z=fenlei(t,k) Z=t;for i=1:size(t)for j=1:size(t)if t(i,j)=k Z(i,j)=1; else Z(i,j)=0;endendend end程序 3：樣本的距離矩陣 function Y=dmatrix(X,a) Y=zeros(size(X,1);for i=1:size(X,1)for j=1:size(X,1)for k=1:size(X,2)Y(i,j)= Y(i,j)+a(1,k)*(X(i,k)-X(

3、j,k)2;endY(i,j)=sqrt( Y(i,j);endend end程序 4：尋找能合并的樣品的標號D=z;a=min(z*eye(size(D);for i=1:size(D)for j=i:size(D)if D(i,j)=adisp(i j)endendend程序 5：主成分分析functioA(X)p,n=size(X);for j=1:nmju(j)=mean(X(:,j);sigma(j)=sqrt(cov(X(:,j);endfor i=1:pfor j=1:nY(i,j)=(X(i,j)-mju(j)/sigma(j);endendsigmaY=cov(Y);T,la

4、mbda=eig(sigmaY); disp(特征根(由小到大):)disp(lambda)disp(特征向量:) disp(T)Xsum=sum(sum(lambda,2),1); for i=1:nfai(i)=lambda(i,i)/Xsum;endfor i=1:npsai(i)= sum(sum(lambda(1:i,1:i),2),1)/Xsum;enddisp(方差貢獻率:) disp(fai)disp(累計方差貢獻率:)disp(psai)程序 6：多元線性回歸 x=d;Y=m;X=ones(length(Y),1) x;b br rb,b,ssss=regress(Y,X);

5、rcoplot(r,r)程序 7：層次分析法 C+程序 #include#includeusing namespatd;const double m=1e-5;constMax=11; /定義最大階const double RI10=0,0,0.58,0.90,1.12,1.24,1.32,1.41,1.45,1.49;class Aprivate:unsigned x,y,n;doubMaxMax;double vMax,wMax,uMax,fMax,hMax;public:void creat(); void strum(); void sum();void open();void A:c

6、reat()coutn;x=n; y=n;for(unsigned i=0;ix;i+)cout輸入第i+1行 for(unsigned j=0;jBij;void A:strum()double s,e=0,t=0,b,r; unsigned i,j,c; for(c=0;cn-1;c+) uc=0;un-1=1;do b=0; r=0;for(i=0;i=n-1;i+)vi=0;wi=ui; for(j=0;jb)b=fabs(vi);for(i=0;ir)r=fabs(ui-wi);while(r=m);cout主特征根 T=bn;for(i=0;i=x-1;i+)cout特征向量i+1

7、為: ;coutuicoutn; s=(b-n)/(n-1);coutnCI=sn;coutCR=s/RIn-1n; for(i=0;i=n-1;i+)for(j=0;j=n-1;j+) t+=Bij ;for(i=0;i=n-1;i+) hi=0;fi=0;for(j=0;j=n-1;j+) hi+=Bij;fi=hi/t;e+=fi;for(i=0;i=n-1;i+)coutn權重向量i+1為: fi/e; coutn;void A:sum() unsigned i,j;coutn 請輸入所要乘的一維數組: ; for(i=0;ivi;for(i=0;ix;i+)ui=0;for(j=0;jy;j+)ui+=Bij*vj;coutn 總排序向量i+1為: ui;coutn;void A:open()coutx;couty;for(unsigned i=0;ix;i+)cout輸入第i+1行:; for(unsigned j=0;jBij;main()k;A a1;cout求主