利用for循環(huán)計算三維矩陣的平方和，其計算需要數(shù)小時如何解決？

我應(yīng)該

f = X^2+Y^2+Z^2;
grid_x = ones(281,281,151);
grid_y  =   ones(281,281,151);
grid_z  =   ones(281,281,151);
for i = 1:151   
    X = grid_x(:,:,i);
    Y = grid_y(:,:,i);
    Z = grid_z(:,:,i);
    C(:,:,i) = vpa(subs(f),20);
    clear X;
    clear Y;
    clear Z;
end