bins=[0:L-1];[m,n]=size(x0);len=m*n;x=reshape(x0,len,1);xpdf=hist(x,bins);tr=round(xpdf*triu(ones(L))*(L-1)/len); 求解释上述程序。。。