matlab 现有一公式,想用matlab实现它,首先有一个矩阵Z,然后Zk表示矩阵的每一行,u表示每一行的均值,n表示

Cr = 0;
n = size(Z, 1);
for i = 1 : n
x = Z( i, : ) - mean( Z(i, :));
Cr = Cr + x * x';
end
Cr = Cr / ( n - 1);