级数与泰勒展式.ppt

级数与泰勒展式 1、求数项级数的和 2、求函数项级数的和 分析：按定义，如果数项级数的部分和序列收敛，则称该级数收敛，称该极限值为级数的和。 函数项级数在固定其通项定义域内每一点后，就是一个数项级数，级数的收敛点集构成和函数的定义域。 在MATLAB中，用于级数求和的命令为：symsum(expr,v,a,b) expr是级数的通项表达式；v是通项中的求和变量，缺省时对默认变量求和； a和b分别为求和变量的起点和终点，缺省时…？ syms n I1=symsum(1/n^2,1,inf) I1 = pi^2/6 syms n x f=(-1)^(n-1)*x^n/n; I2=symsum(f,n,1,inf) %此处的变量“n”不可省略，表明对变量n求和 I2 = piecewise([x == -1, -Inf], [abs(x) = 1 and x ~= -1, log(x + 1)]) %这是一个分段函数，当x=-1时，级数和为-?；当时|x|?1且x?-1，级数收敛于和函数 ln(1+x)。 % 在计算函数项级数时，计算机一般是给出级数收敛时的结果，除此之外级数发散。 syms x n f=sin(x*n)/n^2; sf=symsum(f,n,1,inf); df=[]; for k=0:0.1:4*pi y=subs(sf,k); df=[df y]; end t=0:0.1:4*pi; plot(t,df,’*’) 和函数是以2?为周期的连续且逐项可微函数。 syms x n f=sin(x*n)/n; sf=symsum(f,n,1,inf); df=[]; for k=0:0.1:4*pi y=subs(sf,k); df=[df y]; end t=0:0.1:4*pi; plot(t,df,’*’) 和函数是以2?为周期的间断函数。 1、求函数sin(x）在x=0点的泰勒展式 2、求函数sin3(x） 在x=1点的泰勒展式 在MATLAB中，用于求泰勒展式的指令为 taylor(function,n,x,a) function是待展开的函数表达式，n为展开项数，缺省时展开至5次幂，即6项，x是function中的变量，a为函数的展开点，缺省时为0，即为麦克劳林展开。 syms x taylor(sin(x)) ans = x^5/120 - x^3/6 + x taylor((sin(x))^3,3,x,1) ans = sin(1)^3 + 3*cos(1)*sin(1)^2*(x - 1) - ((3*sin(1)^4 - 6*cos(1)^2*sin(1)^2)*(x - 1)^2)/(2*sin(1)) syms x; f=sin(x); a=taylor(f,x,Expansionpoint,0,Order,6) a = x^5/120 - x^3/6 + x b=taylor(f,x,Expansionpoint,0,Order,20) b = - x^19/121645100408832000 + x^17/355687428096000 - x^15/1307674368000 + x^13/6227020800 - x^11+ x^9/362880 - x^7/5040 + x^5/120 - x^3/6 + x syms x f=sin(x); a=taylor(f,6); b=taylor(f,10); plotf=ezplot(f,[0,2*pi]); hold on plota=ezplot(a,[0,2*pi]); set(plota,color,red) plotb=ezplot(b,[0,2*pi]); set(plotb,color,cyan); hold off legend(‘sin(x)’,‘a’,‘b’,‘Location’,‘North’)%图标 title(‘Taylor Series Expansion’)%图形名称 结论：在[0]随着展开阶数的增加，泰勒多项式逼近的效果更好。在非常接近展开点（零点）附近，低阶泰勒展式的逼近效果已经很好。 syms x %符号运算 f=sin(x); a=taylor(f,6); b=taylor(f,10); t=-6.4:0.2:6.4; %赋值，转入数值运算 f1=subs(f,t); a1=subs(a,t); b1=subs(b,t); plot(t,f1,’*’,t,a1,‘-r’,t,b1,‘c’), axis([-6.5,6.5,-