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

2744557306

这段 MATLAB 代码实现了使用二阶和四阶的泰勒展开法（Taylor's method）来数值求解常微分方程初值问题，并将结果与精确解进行比较。以下是代码的主要解释：
: {9 f1 ~2 Q8 B0 O$ N$ b+ cfunction y = Taylor(a, b, N, af)
    h = (b - a) / N;
    x(1) = a;
/ [) }, N: \. J    y2(1) = af;
    y4(1) = af;
    jqj(1) = af;
7 I+ Q, A1 Z0 V& l6 y% H$ u" t
    for i = 2:N
: i" w4 Q- [8 p$ a( X        % 二阶Taylor法
        y2(i) = y2(i-1) + h * ((1 - h/2) * (x(i-1) - y2(i-1)) + 1);
9 r+ f; P  H' \8 _, x/ ^" V) ~
        % 四阶Taylor法
& S8 R+ }  |: w# [- K2 R        y4(i) = y4(i-1) + h * ((1 - h/2 + h^2/6 - h^3/24) * (x(i-1) - y4(i-1)) + 1);
# w) r# p4 J  z  v! q# L8 f0 C$ ?
        x(i) = a + (i-1) * h;
        jqj(i) = x(i) + exp(-x(i));
    end
: ?9 }; ?4 \5 Z# q% K
    % 将结果输出为矩阵形式
    result = [x', y2', y4', jqj'];
1 ~8 q7 f3 C5 ~' J9 F$ d5 r
    % 绘制曲线
    plot(x, y2, 'r', x, y4, 'b', x, jqj, 'g');
" l9 X5 \& _4 d. ~2 E    legend('Taylor2法', 'Taylor4法', '精确解');
* l' X8 q; j' Y. \end

6 P; p* m4 d6 [! ]  ~& @该函数的输入参数包括求解区间 [a, b]、步数 N 以及初始条件 af。函数使用二阶和四阶的泰勒展开法来逐步求解微分方程。最后，结果以矩阵形式输出，并绘制了二阶Taylor法、四阶Taylor法和精确解的曲线。
3 |2 D1 t7 t6 z5 n你可以调用这个函数，例如 Taylor(0, 1, 100, 1)，其中参数表示解在区间 [0, 1] 上，步数为 100，初始条件为 1。
3 }3 Y, i; D# F2 s9 I$ c' V
  \: O- f$ u4 k