求Matlab大神,急.

求Matlab大神,急.
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);
求解释上述程序.

第2行：获取xo矩阵的行数m和列数n;
第三行：xo的元素个数len
第4行：x是一个列向量,其元素是xo的所有元素
第5行：按照bins指定的中心点和bins的元素个数画统计频率直方图,并返回每个频段上的元素个数xpdf是行向量；
第6行:总的意思是想得出一个累积频率乘以（L-1）,结果四舍五入得到tr.他这个tr是递减的,因为他用xpdf乘上了一个单位阵的上三角阵,所以第一个值包含所有xo的所有元素数,第二个值包含从0到L-1的元素数,第三个值包含从2到L-1的元素数,然后将这些值都乘以L-1.具体什么意思你应该知道.