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

2744557306

这段 MATLAB 代码实现了使用二阶和四阶的泰勒展开法（Taylor's method）来数值求解常微分方程初值问题，并将结果与精确解进行比较。以下是代码的主要解释：
function y = Taylor(a, b, N, af)
    h = (b - a) / N;
; V& m% z# O7 K' }    x(1) = a;
" ^1 y' b: Y% }8 @% n    y2(1) = af;
5 ~! _5 u1 h8 Y1 w    y4(1) = af;
( l3 i' ?9 }; k    jqj(1) = af;
" T& b; ?& d1 u8 c2 R( b
% N( }+ J' Q! j    for i = 2:N
1 W# k/ x( c/ L0 L2 m: N  _        % 二阶Taylor法
        y2(i) = y2(i-1) + h * ((1 - h/2) * (x(i-1) - y2(i-1)) + 1);
- c0 D7 n% E$ Z$ f
        % 四阶Taylor法
" L6 t3 x; `# J: h; L! T+ C" F        y4(i) = y4(i-1) + h * ((1 - h/2 + h^2/6 - h^3/24) * (x(i-1) - y4(i-1)) + 1);
  K. ~" ]8 Y9 f) `
4 }4 Z6 X3 k, r6 e" W6 p! \& Q        x(i) = a + (i-1) * h;
        jqj(i) = x(i) + exp(-x(i));
( K* E. p% D$ h* c: W    end
4 H3 f# {& v, ^  h& v! H. R
2 }5 y3 S% f  C, }; i    % 将结果输出为矩阵形式
    result = [x', y2', y4', jqj'];
/ m0 R( \# ?. z$ E! F+ l* |. h
    % 绘制曲线
    plot(x, y2, 'r', x, y4, 'b', x, jqj, 'g');
% m4 q' X2 E: n  Y1 R    legend('Taylor2法', 'Taylor4法', '精确解');
end
3 z# J+ _" R( A
该函数的输入参数包括求解区间 [a, b]、步数 N 以及初始条件 af。函数使用二阶和四阶的泰勒展开法来逐步求解微分方程。最后，结果以矩阵形式输出，并绘制了二阶Taylor法、四阶Taylor法和精确解的曲线。
+ r6 T9 F: }& x9 |3 J你可以调用这个函数，例如 Taylor(0, 1, 100, 1)，其中参数表示解在区间 [0, 1] 上，步数为 100，初始条件为 1。

. {/ ]8 g: z. E6 F' R# _" R- q