我有以下循环
syms k m
summ = symsum(k^2,k,1,m)
for j=1:100
summ = subs(symsum((m-j)^2*summ,m,j+1,k),k,m);
end
在R2014a、W7、64b上计算大约需要600s。Mathematica 中的等效计算需要 60 秒。结果是相同的,但有区别。所有和的属性之一(在这种情况下全部为一百个)是它可以以一种半好的方式分解。Mathematica 默认不分解它,因为它没有被告知这样做。但是 symsum 通过检查总和来做到这一点,我认为这是放缓的根源。
你会同意吗?如果是,有没有办法告诉 symsum 不要简化输出?如果没有,有没有办法加快速度?我真的很想在 Matlab 中做而不是在 Mathematica 中。