[TTT_IOU]

对于任意的矩阵或向量A：
如：
A=magic(5);
[i,j]=find(A<=5);
L=length(x);
B=diag(A(x,y));
希望有帮助