第1个回答 2019-06-27
符号计算的话用symsum函数。比如求和∑1/n^2(n从1到k)
syms
n
k
>>
s=symsum(1/n^2,n,1,k)
s
=
-Psi(1,k+1)+1/6*pi^2
再比如求和∑n^2(n从1到k)
syms
n
k
>>
s=symsum(n^2,n,1,k)
s
=
1/3*(k+1)^3-1/2*(k+1)^2+1/6*k+1/6
s=simplify(s)
s
=
1/3*k^3+1/2*k^2+1/6*k
而要用数值计算的话,应当用一个矩阵先把要求和的元素存起来,然后用sum函数。比如上面的两个例子
第一个为:
k=1000;
a=1./(1:k).^2%存储1,1/4,1/9,1/16...
s=sum(a)
s
=
1.6439
第二个为:
k=1000;
a=(1:k).^2;%存储1,4,9,16...
s=sum(a)