三维数据的处理

[phy1987]

请各位高手帮忙！
Z1=T11(:);
zmean=mean(Z1);
zstd=std(Z1);
a=cumsum(T11)*0.1;
b=max(a);
m=max(cumsum(b)*0.1);
c=diag(diag(X1));
mx=max(cumsum(b*c)*0.1);
m1=mx/m;
c2=diag(diag(Y1));
my=max(cumsum(b*c2)*0.1);
m2=my/m;
d=(T11*diag(X1)*0.1);
mxy=d*diag(Y1)*0.1;
covxy=mxy/m-m1*m2
e=b*c^2;
f=max(cumsum(e)*0.1);
bzlcx=sqrt(f/m-m1^2);
cofxy=covxy/bzlcx*bzlcy;
这是小弟要用的的程序，想请高手帮忙！帮我看一下是哪里出现错误，谢谢！我要处理的是231*46矩阵。