散点矩阵之间和内部

Between and within scatter matrix

有人可以帮我为这两个矩阵创建一个 matlab 代码吗?

this are the formulas

w = 0;

for i = 1:k
    for j = 1:N(i)
        x_diff = x(i,j) - mean(x(i,:));
        w = w + (x_diff * x_diff^T);
    end
end

b = 0;

for i = 1:k
    x_diff = mean(x(i,:)) - mean(x);
    b = b + (N(i) * x_diff * x_diff^T);
end