声明GydF4y2Ba
符号x f=1/(5+4*cos(x));T=泰勒(f,'顺序',8)GydF4y2Ba
返回GydF4y2Ba
T=(49*x^6)/131220+(5*x^4)/1458+(2*x^2)/81+1/9GydF4y2Ba
这是泰勒级数中8阶的所有项,但不包括GydF4y2BaFGydF4y2Ba(GydF4y2Ba十GydF4y2Ba):GydF4y2Ba
从技术上讲,GydF4y2BaTGydF4y2Ba
是麦克劳林系列GydF4y2Ba,因为它的膨胀点是GydF4y2Baa=0GydF4y2Ba
.GydF4y2Ba
这些命令GydF4y2Ba
syms x g=exp(x*sin(x));t=taylor(g,'展开点',2,'顺序',12);GydF4y2Ba
生成泰勒级数的前12个非零项GydF4y2BaGGydF4y2Ba
关于GydF4y2Bax=2GydF4y2Ba
.GydF4y2Ba
TGydF4y2Ba
是一个大的表达;进入GydF4y2Ba
大小(字符(t))GydF4y2Ba
ans=1 99791GydF4y2Ba
找到那个GydF4y2BaTGydF4y2Ba
它的印刷体大约有100000个字符。为了继续使用GydF4y2BaTGydF4y2Ba
,首先简化其表示:GydF4y2Ba
t=简化(t);大小(字符(t))GydF4y2Ba
ans=1 6988GydF4y2Ba
接下来,将这些函数绘制在一起,以查看泰勒近似与实际函数的比较情况GydF4y2BaGGydF4y2Ba
:GydF4y2Ba
xd=1:0.05:3;yd=subs(g,x,xd);fplot(t,[1,3])保持GydF4y2Ba在GydF4y2Ba绘图(xd,yd,GydF4y2Ba'r-.'GydF4y2Ba)职务(GydF4y2Ba'泰勒近似与实际函数'GydF4y2Ba)图例(GydF4y2Ba“泰勒”GydF4y2Ba,GydF4y2Ba'函数'GydF4y2Ba)GydF4y2Ba
特别感谢瑞典UMEA的Gunnar Bäckstrøm教授提供了这个例子。GydF4y2Ba