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

2744557306

这段 MATLAB 代码实现了使用二阶和四阶的泰勒展开法（Taylor's method）来数值求解常微分方程初值问题，并将结果与精确解进行比较。以下是代码的主要解释：
function y = Taylor(a, b, N, af)
    h = (b - a) / N;
    x(1) = a;
    y2(1) = af;
% w9 Y7 f  O3 O0 b. a    y4(1) = af;
! @" [% c/ ]# k' Z. j    jqj(1) = af;
7 C& {9 m% J3 d& e  a' m, A
    for i = 2:N
        % 二阶Taylor法
' ]3 H1 \, O' O5 {5 m) a& f        y2(i) = y2(i-1) + h * ((1 - h/2) * (x(i-1) - y2(i-1)) + 1);
( k& G; z8 L1 n" }' |0 V6 o
        % 四阶Taylor法
2 A3 Y% v# a/ ~; e$ L        y4(i) = y4(i-1) + h * ((1 - h/2 + h^2/6 - h^3/24) * (x(i-1) - y4(i-1)) + 1);
4 M7 S( a% U* V8 v' c
        x(i) = a + (i-1) * h;
        jqj(i) = x(i) + exp(-x(i));
& E7 U; x: h4 C$ r) ]    end
% n' Q& O% T6 v& O0 L: x1 y
7 J8 D4 L/ l2 i3 n$ |    % 将结果输出为矩阵形式
    result = [x', y2', y4', jqj'];

    % 绘制曲线
9 Y1 D+ P# h- Z- F. s9 o" j1 w# ~    plot(x, y2, 'r', x, y4, 'b', x, jqj, 'g');
    legend('Taylor2法', 'Taylor4法', '精确解');
! a% K) _9 T5 i+ Y! V& Mend
* |& k. x; K) b) N
该函数的输入参数包括求解区间 [a, b]、步数 N 以及初始条件 af。函数使用二阶和四阶的泰勒展开法来逐步求解微分方程。最后，结果以矩阵形式输出，并绘制了二阶Taylor法、四阶Taylor法和精确解的曲线。
3 c" @( b4 j- ?2 Y8 e9 U* f你可以调用这个函数，例如 Taylor(0, 1, 100, 1)，其中参数表示解在区间 [0, 1] 上，步数为 100，初始条件为 1。

$ S. B9 o* o2 m
3 y9 b3 F% P" }0 k! ^" S