无穷级数的计算.PDFVIP

无穷级数的计算 MATLAB 求解数列部分和符号运算指令为 symsum(s,t,a,b)—表示s 中的符号变量t 从a 到b 的级数和（t 缺省时设定为x 或 最接近 x 的字母） k 例1：计算∑n 。 n 1 syms k symsum(k) 结果显示 ans = 1/2*k^2-1/2*k 3 1 ∞ 1 ∞ 1 ∞ 1 例2 ： 、 、 和 ∑ ∑ ∑ ∑ k 1 k k 1 k 2 k 1 k 3 k 0 x k syms x k s=symsum(1/x,1,3) 结果显示 ans s=11/6 s1=symsum(1/k^2,k,1,inf) 结果显示 s1 =1/6*pi^2 s3=symsum(1/k^3,k,1,inf) 结果显示 s3 = zeta(3) vpa(zeta(3)) ans = 1.2020569031595942366408280577161 批注： help zeta ZETA Symbolic Riemann zeta function. ZETA(z) = sum(1/k^z,k,1,inf). ZETA(n,z) = n-th derivative of ZETA(z) Overloaded methods help sym/zeta.m help vpa VPA Variable precision arithmetic. R = VPA(S) numerically evaluates each element of the double matrix S using variable precision floating point arithmetic with D decimal digit accuracy, where D is the current setting of DIGITS. The resulting R is a SYM. VPA(S,D) uses D digits, instead of the current setting of DIGITS. D is an integer or the SYM representation of a number. Examples: vpa(pi,780) shows six consecutive 9s near digit 770 in the decimal expansion of pi. vpa(hilb(2),5) returns [ 1., .50000] [.50000, .33333] See also DOUBLE, DIGITS. Overloaded methods help sym/vpa.m S4=symsum(x^k,k,0,inf) 结果显示 S4 =-1/(x-1) ∞ 1 例 3：通过编写MATLAB 程序计算p-级数∑ p 的部分和数列，观测部分和数列的变化趋 n 1 n 势。 （1）当p 1 时，级数收敛。设p 2 ，p 3 ，观察部分和数列的变化趋势。 for n=1:50 s1=0; s2=0 for i=1:n s1=s1+1/i^2; s2=s2+1/i^3; ss(n)=s1; sss(n)=s2; end end ss 结果显示 ss= 1 .0000 1.2500 1.3611 1.4236 1.4636 1.4914 1.5118 1.5274 1.5398 1.5498 1.5580