数组的质心，可以看为将数组值作为像素点灰度值时的重心。

在文章Matlab质心算法中提到了质心的算法。

第一种算法源自文章Matlab质心算法：

% B为待求质心的数组
x = 0;
y = 0;
sum1 = 0;
for i = 1 : size(B, 1)
	for j = 1 : size(B, 2)
		y = y + B(i, j) * i;
		x = x + B(i, j) * j;
		sum1 = sum1 + B(i, j);
	end
end
x = x / sum1;
y = y / sum1;

文章Matlab质心算法中也给出了另一种算法：求每一行/列的质心，然后再求数组的质心。

但在本文中有一个更简单的算法，代码更简洁：

% B为待求质心的数组
x1 = sum(sum(B .* [1 : size(B, 2)])) / sum(B(:));
y1 = sum(sum(B' .* [1 : size(B, 1)])) / sum(B(:));

欢迎大家批评指正。