正弦函数的泰勒展开式

2744557306

1. x0=-2*pi:0.01:2*pi; y0=sin(x0); syms x; y=sin(x);
   $ i, @0 Z+ M0 Q* n$ q6 o
2. 
     Y( [+ Y+ K8 Y+ q+ f% _$ l6 l
3.       plot(x0,y0), axis([-2*pi,2*pi,-1.5,1.5]); hold on& t1 Q  d$ K5 e5 N\" n- O
   
4. 
   + ?' J1 u, x5 h6 A% t\" g
5.       for n=[8:2:16]  }5 z; q; t4 P  y& U, w
   
6. 
   / ^' W! A5 k8 a1 {' `
7.          p=taylor(y,x,n), y1=subs(p,x,x0); line(x0,y1)
   & ]3 J7 R% `8 ?8 K8 K; f, [
8. 
   9 |/ z8 Y. w1 c0 O9 U8 i
9.       end

复制代码

这段代码主要是用来绘制正弦函数在给定区间 [-2π, 2π] 上的泰勒展开多项式，并将展开多项式在不同阶数下的近似结果与原始正弦函数的图像进行比较。
# C+ y* v) _5 ]! D4 a7 m
首先，代码中定义了一个变量 `x0`，表示从 -2π 到 2π，间隔为 0.01 的一组 x 值。然后计算了这些 x 值对应的正弦函数值，存储在变量 `y0` 中。
6 k( f% Y# r; O, f9 P
" ^! F% k3 X1 O! O" M2 r/ |接着，定义了符号变量 x 和正弦函数 y=sin(x)。然后绘制了原始正弦函数在区间 [-2π, 2π] 上的图像，并设置了坐标轴的范围为 x 轴从 -2π 到 2π，y 轴从 -1.5 到 1.5。
7 h( a* D$ V) H8 o$ p$ l
接下来，通过循环遍历了 n 从 8 到 16，步长为 2 的一系列值。在每次循环中，使用 `taylor` 函数对正弦函数 y=sin(x) 进行 n 阶的泰勒展开，并将展开多项式存储在变量 p 中。然后通过 `subs` 函数将这个展开多项式在 x0 上求值，得到对应的近似值，存储在变量 y1 中。最后，使用 `line` 函数将 x0 和 y1 连线，表示泰勒展开多项式在不同阶数下的近似结果。

3 y) e! b, t; U0 N通过这段代码，可以直观地观察到随着泰勒展开多项式阶数的增加，对正弦函数的逼近效果会越来越好。
) B) f% ~6 z6 o4 K% t. K


4 ?9 r6 n+ a  S' v) d" P! ^