数学实验(MATLAB版旱墨明版)2.5-3.4部分答案.doc

练习2.5 判断下列级数的敛散性，若收敛，求出其收敛值. ，，，，，，. 解：编程： （1） syms n s=symsum(1/n^(2^n),n,1,inf) s1 =sum(1/((n^2)^n),n = 1 .. Inf) （2） syms n s=symsum(sin(1/n),n,1,inf) s =sum(sin(1/n),n = 1 .. Inf) （3） syms n s=symsum(log(n)/n^3,n,1,inf) s = -zeta(1,3) （4） syms n s=symsum(1/(log(n))^n,n,3,inf) s = sum(1/(log(n)^n),n = 3 .. Inf) （5） syms n s=symsum(1/(n*log(n)),n,2,inf) s = sum(1/n/log(n),n = 2 .. Inf) （6） syms n s=symsum((-1)^n*n/(1+n^2),n,1,inf) s = -1/2*hypergeom([2, 1+i, 1-i],[2-i, 2+i],-1) （7） 显然，上面级数（1）-（6）都收敛，分别等于：sum(1/(n^(2^n)),n = 1 .. Inf)、sum(sin(1/n),n = 1 .. Inf)、-zeta(1,3)、sum(1/(log(n)^n),n = 3 .. Inf)、sum(1/n/log(n),n = 2 .. Inf)、symsum((-1)^n*n/(1+n^2),n,1,inf)。 求当k=4,5,6,7,8时公式中t的值。 解：当k=4时，编程： syms n s=symsum(1/n^8,n,1,inf) s = 1/9450*pi^8 此时，，，再次编程： syms t solve(1/9450*pi^8=pi^8/t,t) ans =9450 解得：t=9450; 当k=5时，编程： syms n s=symsum(1/n^10,n,1,inf) s =1/93555*pi^10 此时，，，再次编程： syms t solve(1/93555*pi^10=pi^10/t,t) ans = 93555 解得：t=93555; 当k=6时，利用上面步骤直接编程有： syms n s=symsum(1/n^12,n,1,inf) s = 691/638512875*pi^12 syms t solve(691/638512875*pi^12=pi^12/t,t) ans =638512875/691 解得：t=638512875/691; 当k=7时，利用上面步骤直接编程有： syms n s=symsum(1/n^14,n,1,inf) s = 2pi^14 syms t solve(2pi^14=pi^14/t,t) ans 2 解得：t2; 当t=8时，利用上面步骤直接编程有： syms n s=symsum(1/n^16,n,1,inf) s = 3617/325641566250*pi^16 syms t solve(3617/325641566250*pi^16=pi^16/t,t) ans =325641566250/3617 解得：t=325641566250/3617. 用Taylor命令观测函数的Machlaurin展开式的前几项，然后在同一坐标系里作出函数和它的Taylor展开式的前几项构成的多项式函数的图形，观测这些多项式函数的图形的图形逼近的情况. ;(2);(3); ;(5);(6). 解：（1）编程： syms x f=asin(x); t1=taylor(f,1);t2=taylor(f,2);t3=taylor(f,3);t4=taylor(f,4); t5=taylor(f,5);t6=taylor(f,6); t1,t2,t3,t4,t5,t6 t1 = 0 t2 = x t3 =x t4 =x+1/6*x^3 t5 =x+1/6*x^3 t6 =x+1/6*x^3+3/40*x^5 然后作图观察图形逼近情况： x=-1:0.1:1; f=asin(x);t2=x;t4=x+1/6*x.^3;t6=x+1/6*x.^3+3/40*x.^5; plot(x,f,x,t2,*,x,t4,+,x,t6,o) grid on 图像： 通过图像可知当：它们之间逼近程度最高。 （2）编程： syms x f=atan(x);t1=taylor(f,1);t2=taylor(