MathLabs
Định lýĐã chứng minh

Chặn sai số toàn cục của phương pháp Euler

Phát biểu

Giả sử f(t,y)f(t,y) liên tục Lipschitz theo yy với hằng số LL, và nghiệm chính xác thỏa ∣y′′(t)∣≤M|y''(t) | \le M trên [t0,T][t_0, T]. Khi đó sai số toàn cục en=y(tn)−yne_n = y(t_n) - y_n sau nn bước với bước nhảy hh thỏa ∣en∣≤hM2L(eL(tn−t0)−1)|e_n| \le \frac{hM}{2L}\left(e^{L(t_n-t_0)}-1\right).

Vì sao đúng?

Mỗi bước Euler mắc một sai số cục bộ nhỏ cỡ O(h2)O(h^2) do cắt bỏ chuỗi Taylor sau số hạng bậc nhất, nhưng các sai số cục bộ này có thể cộng dồn qua khoảng (T−t0)/h(T-t_0)/h bước cần thiết để tới thời điểm TT cố định. Định lý cho thấy sự cộng dồn này không quá tệ: điều kiện Lipschitz ngăn các sai số cục bộ nhỏ bị khuếch đại nhanh hơn hàm mũ, nên sai số toàn cục chỉ còn bậc nhất, O(h)O(h), theo bước nhảy — kém một bậc so với sai số cục bộ, đúng theo quy luật thường gặp ở các phương pháp một bước.

Phác thảo chứng minh

Viết sai số cục bộ mắc phải trong một bước là hiệu giữa nghiệm chính xác khai triển bằng chuỗi Taylor và bước cập nhật Euler. Định lý Taylor cho y(tn+1)=y(tn)+hf(tn,y(tn))+h22y′′(ξn)y(t_{n+1}) = y(t_n) + h f(t_n, y(t_n)) + \frac{h^2}{2} y''(\xi_n) với ξn\xi_n nào đó nằm giữa tnt_n và tn+1t_{n+1}, nên sai số cắt cụt cục bộ là τn=h22y′′(ξn)\tau_n = \frac{h^2}{2} y''(\xi_n), bị chặn bởi h2M2\frac{h^2 M}{2} nhờ giả thiết ∣y′′(t)∣≤M|y''(t) | \le M.

Trừ bước cập nhật Euler yn+1=yn+hf(tn,yn)y_{n+1} = y_n + h f(t_n, y_n) khỏi khai triển Taylor chính xác này. Viết en=y(tn)−yne_n = y(t_n) - y_n, hiệu của hai vế phải cho en+1=en+h[f(tn,y(tn))−f(tn,yn)]+τne_{n+1} = e_n + h\left[f(t_n, y(t_n)) - f(t_n, y_n)\right] + \tau_n. Điều kiện Lipschitz chặn số hạng trong ngoặc vuông bởi L∣en∣L|e_n|, nên ∣en+1∣≤(1+hL)∣en∣+h2M2|e_{n+1}| \le (1+hL)|e_n| + \frac{h^2 M}{2}.

Vì e0=0e_0 = 0 (giá trị ban đầu là chính xác), khai triển đệ quy này bằng quy nạp cho ∣en∣≤h2M2∑k=0n−1(1+hL)k=hM2L[(1+hL)n−1]|e_n| \le \frac{h^2M}{2}\sum_{k=0}^{n-1}(1+hL)^k = \frac{hM}{2L}\left[(1+hL)^n - 1\right], dùng đẳng thức chuỗi hình học. Cuối cùng, (1+hL)n≤eLhn=eL(tn−t0)(1+hL)^n \le e^{Lhn} = e^{L(t_n-t_0)} vì 1+x≤ex1+x\le e^x với mọi số thực xx, biến chặn trên thành ∣en∣≤hM2L(eL(tn−t0)−1)|e_n| \le \frac{hM}{2L}\left(e^{L(t_n-t_0)}-1\right) — đúng là bất đẳng thức cần chứng minh, và rõ ràng là O(h)O(h) khi tnt_n cố định.

Chủ đề chứa định lý này

Chứng minh từng bước

Chưa có chứng minh từng bước cho định lý này.

Tài liệu tham khảo

  1. J. C. Butcher (2016). Numerical Methods for Ordinary Differential Equations
  2. E. Hairer, S. P. Norsett, G. Wanner (1993). Solving Ordinary Differential Equations I: Nonstiff Problems