【正文】
lag=1。else R=B。%循環(huán)計(jì)算R傳遞閉包endend程序四:求截矩陣:function [D k] =jjz4(B)L=unique(B)39。a=size(B)。D=zeros(a)。for m=length(L):1:1 k=L(m)。 for i=1:a for j=1:a if B(i,j)=k D(i,j)=1。 else D(i,j)=0。%求?截距陣,當(dāng)bij≥? 時(shí),bij(?) =1;當(dāng)bij<? 時(shí),bij(?) =0 end end endfprintf(39。當(dāng)分類系數(shù)k=:\n39。)。disp(L(m))。fprintf(39。所得截距陣為:\n39。)。disp(D)。e