二阶和四阶的泰勒展开法（Taylor's method）来数值求解常微分方程初值问题

2744557306

这段 MATLAB 代码实现了使用二阶和四阶的泰勒展开法（Taylor's method）来数值求解常微分方程初值问题，并将结果与精确解进行比较。以下是代码的主要解释：
/ @/ j& ?1 p7 n# R, x) c7 z- {function y = Taylor(a, b, N, af)
    h = (b - a) / N;
" T# r/ P( k; Z- d$ r2 @    x(1) = a;
9 w) A9 E0 a: S, F6 n* V9 g# v8 F    y2(1) = af;
    y4(1) = af;
    jqj(1) = af;
1 V: G5 D9 Y: A. ~$ X$ U9 J
2 J5 u) S* `" A( B6 ~5 U& G    for i = 2:N
5 I1 Y& F0 [" ~  w        % 二阶Taylor法
        y2(i) = y2(i-1) + h * ((1 - h/2) * (x(i-1) - y2(i-1)) + 1);

        % 四阶Taylor法
! F; ]' x/ Q2 h' |. Z        y4(i) = y4(i-1) + h * ((1 - h/2 + h^2/6 - h^3/24) * (x(i-1) - y4(i-1)) + 1);
2 @: Q- }+ q. d7 H7 V8 G7 r' @
' I+ r. x& L* h1 x( N        x(i) = a + (i-1) * h;
/ I6 G  y, w2 Z1 H        jqj(i) = x(i) + exp(-x(i));
& y5 t; V& C. O* {0 B    end
2 \8 ~; R' p  @! l
    % 将结果输出为矩阵形式
- R7 g# r3 d( }6 `7 }3 W% n6 x  q# H    result = [x', y2', y4', jqj'];
- h7 ]/ r9 V0 {+ c1 G0 s0 C  [* o
  p7 ]: [+ L: a6 K, A    % 绘制曲线
    plot(x, y2, 'r', x, y4, 'b', x, jqj, 'g');
    legend('Taylor2法', 'Taylor4法', '精确解');
$ C( k4 M' c( r' ^. n$ B: N8 H9 Nend
9 f5 p) J+ Y( D9 H# L- {
+ I+ u: }$ X; w* F1 y该函数的输入参数包括求解区间 [a, b]、步数 N 以及初始条件 af。函数使用二阶和四阶的泰勒展开法来逐步求解微分方程。最后，结果以矩阵形式输出，并绘制了二阶Taylor法、四阶Taylor法和精确解的曲线。
' \, j7 z3 s4 g" ?# M6 w# W# R你可以调用这个函数，例如 Taylor(0, 1, 100, 1)，其中参数表示解在区间 [0, 1] 上，步数为 100，初始条件为 1。

: B- K: c% C$ Y* G, i  w