09实验九 无穷级数.pptVIP

图9.8 MATLAB 高等数学实验 实验九 无穷级数 实验目的 掌握用MATLAB求无穷级数的和、求幂级数的收敛域、展开函数为幂级数以及展开周期函数为傅里叶级数的方法。 9.1 学习MATLAB命令 9.1.1 符号表达式求和函数 symsum(S(k)) 返回符号表达式S(k)中的符号变量k从0到k-1的和值 symsum(S(k),v) 返回符号表达式S(k)中指定的符号变量k由v代替, 再对v从0到v-1的和值 symsum(S(k),a,b) 返回符号表达式S(k)中的符号变量k从a到b的和值 symsum(S(k),v,a,b) 返回符号表达式S(k)中的符号变量k由v代替, 再从a到b的和值 9.1.2 符号函数的泰勒级数展开函数 taylor(f(x))或taylor(y,x) 求函数y=f(x)的5阶麦克劳林展开式 taylor(f(x),n) 求函数f(x)的n-1阶麦克劳林展开式 taylor(f(x),a) 求函数f(x)在x=a处的5阶泰勒展开式 taylor(f(x),n,a) 求函数f(x)在x=a处的n-1阶泰勒展开式 taylor(f,n,v,a) 返回符号表达式f在a点泰勒展开到n-1次式, 自变量为v taylortool(f) 对指定的函数f, 用图形用户界面显示出泰勒展开式 9.1.3 泰勒级数计算器函数 subs(S,old,new) 将符号表达式S中的符号变量old用new代替 9.1.4 在符号表达式或矩阵中 进行符号替换的函数 simplify(expr)与simple(expr) 用于化简符号表达式expr 9.1.5 符号表达式的化简函数 9.2 实验内容 9.2.1 级数求和 当符号变量的和存在时, 可以用symsum命令来求无穷级数和。 【例1】求 的值。 输入: syms n s1=symsum(1/(4*n^2+8*n+3),n,1,inf) 得到该级数的和为s1=1/6。 【例2】求级数 的和。 输入: syms k x s2=symsum(x^(3*k),k,1,inf) 得到和函数 s2 =piecewise([x^3 = 1, Inf], [x^3 1 and abs(x) in Dom::Interval(0, 1), -x^3/(x^3 - 1)], [x^3 1 and (1 abs(x) and not 1 x or abs(x) = 1 and ((abs(x) 1 or 0 x) and x = 1 or not 0 x) and x 1), -(limit(x^(3*k)/(x^3 - 1), k = Inf) - x^3*limit(x^(3*k)/(x^3 - 1), k = Inf) + x^3)/(x^3 - 1)]) 【例3】设 , 求 。 首先输入: for n=1:25 a=1; for m=1:n a=a*10/m; end plot(n,a,*) hold on end 得到输出(见图9.1)。从散点图可见 的变化趋势。 syms n; symsum(10^n/gamma(n+1),n,1,inf) 输出： ans = exp(10) - 1 图9.1 9.2.2 求幂级数的收敛域. 【例4】求 的收敛域与和函数。 输入: clear; syms n x a1=4^(2*n)*(x-3)^n/(n+1); a2=subs(a1,n,n+1); p=limit(a2/a1,n,inf) 输出为： p= 16*x-48 注意, 这里对a2和a1都没有加绝对值。因此上式的绝对值小于1时, 幂级数收敛, 大于1时发散。 为了求出收敛区间的端点, 输入: x1=solve(16*x-48=1) x2=solve(16*x-48=-1) 输出为: x1=49/16 x2=47/16 由此可知 时收敛, 或 时发散。 为了判断端点的敛散性, 输入: simplify(subs(a1,x,49/16)) 得到x为右端点时幂级数的一般项为: ans= 1/(n+1) 因此当x=49/16时发散。 再输入: simplify(subs(a1,x,47/16)) 输出结果为: ans= (-1)^n/(n+1) 因此当x=47/16时,级数收敛