Máy Tính Phương Pháp Euler

Giới thiệu về Máy Tính Phương Pháp Euler

Máy tính Phương pháp Euler giải bằng số bất kỳ bài toán giá trị ban đầu cấp một nào có dạng \( y' = f(x, y), \; y(x_0) = y_0 \) bằng phương pháp Euler cổ điển (tiến). Nó cung cấp một bảng lặp đầy đủ, vẽ đa giác Euler trên trường độ dốc trực tiếp, so sánh lời giải tại ba kích thước bước khác nhau để bạn có thể quan sát trực quan sự hội tụ của phương pháp và — nếu bạn cung cấp lời giải chính xác dạng đóng — nó sẽ thực hiện phân tích sai số từng bước.

Phương pháp Euler là gì?

Phương pháp Euler là thuật toán đơn giản nhất để xấp xỉ lời giải của một bài toán giá trị ban đầu. Bắt đầu từ một điểm đã biết \( (x_0, y_0) \) trên đường cong lời giải, nó lặp lại việc tiến thêm một bước nhỏ có kích thước h dọc theo độ dốc cục bộ \( f(x, y) \):

Về mặt hình học, mỗi bước là một đoạn thẳng ngắn có độ dốc bằng giá trị của phương trình vi phân tại điểm hiện tại. Đường gấp khúc thu được — đa giác Euler — là một xấp xỉ cho lời giải thực (thường là đường cong).

Nó chính xác đến mức nào?

Phương pháp Euler là một phương pháp cấp một. Sai số cắt cụt cục bộ tại mỗi bước là \( O(h^2) \) và sai số toàn cục sau khi tích phân trên một khoảng cố định là \( O(h) \). Trong thực tế:

• Giảm một nửa kích thước bước sẽ làm giảm khoảng một nửa sai số toàn cục.
• Sai số tăng tuyến tính với chiều dài của khoảng tích phân.
• Sai số tồi tệ nhất ở những nơi lời giải có độ cong lớn.

Việc so sánh kích thước bước tích hợp (h, h/2, h/4) cho phép bạn thấy trực tiếp sự hội tụ tuyến tính này: hãy bật tùy chọn này và kiểm tra xem ba giá trị cuối cùng có tiến tới một giới hạn chung với mỗi giá trị cách giới hạn khoảng một nửa so với giá trị trước đó hay không.

Cách đọc biểu đồ

Hình ảnh trực quan hóa lớp bốn loại thông tin trên một mặt phẳng tọa độ duy nhất:

• Trường độ dốc màu xám — các đoạn thẳng ngắn có độ nghiêng bằng \( f(x, y) \) tại điểm đó. Hãy coi nó là "dòng chảy mà ODE quy định". Bất kỳ đường cong lời giải nào cũng phải tiếp tuyến với trường tại mọi điểm.
• Đa giác Euler màu chàm — lời giải số theo từng bước. Mỗi đoạn bắt đầu tại điểm lưới trước đó và hướng theo \( f(x_n, y_n) \) trong một khoảng cách h.
• Đường cong chính xác màu xanh lá đứt nét — chỉ xuất hiện khi bạn cung cấp lời giải dạng đóng. Các đoạn gạch dọc màu cam đứt nét là các sai số cục bộ có dấu \( y_n - y_{\text{chính xác}}(x_n) \).
• Đường cong so sánh màu cam và xanh lá — cùng một bài toán được chạy lại tại h/2 và h/4, được hiển thị khi bật so sánh kích thước bước.

Cách sử dụng máy tính này

1. Nhập vế phải của ODE vào trường được đánh dấu y' =. Sử dụng x và y làm biến. Các toán tử được hỗ trợ là + − × ÷ ^, và các hàm được hỗ trợ bao gồm sin, cos, tan, asin, acos, atan, exp, ln, log, log10, log2, sqrt, abs.
2. Thiết lập các điều kiện ban đầu: giá trị bắt đầu x₀, giá trị y₀ ban đầu tại điểm đó, kích thước bước h (dương để tích phân tiến, âm để tích phân lùi) và số bước n.
3. (Tùy chọn) Cung cấp lời giải chính xác y(x) nếu bạn biết. Máy tính sẽ tính toán \( |y_n - y(x_n)| \) tại mỗi bước và báo cáo sai số tối đa và sai số cuối cùng.
4. Chuyển đổi các tùy chọn trực quan hóa: trường độ dốc được bật theo mặc định; so sánh kích thước bước chồng thêm hai đường cong bổ sung tại h/2 và h/4.
5. Nhấp vào Chạy. Phần kết quả hiển thị số liệu thống kê tóm tắt, biểu đồ, bảng so sánh hội tụ và bảng lặp đầy đủ. Di chuột qua một hàng sẽ làm nổi bật điểm tương ứng trên biểu đồ (và ngược lại).

Ví dụ minh họa

Xét \( y' = x + y, \; y(0) = 1 \) với h = 0.1 và 10 bước. Lời giải chính xác là \( y(x) = -x - 1 + 2e^x \). Áp dụng công thức Euler ta có:

Sai số cuối cùng là khoảng 0.249. Giảm một nửa h xuống 0.05 làm sai số cuối cùng giảm xuống còn khoảng 0.13, và giảm một nửa lần nữa xuống 0.025 làm nó giảm xuống còn khoảng 0.067 — sự hội tụ tuyến tính rõ ràng, đúng như lý thuyết dự đoán.

Phương pháp Euler so với các phương pháp số khác

Khi Euler sai

Phương pháp Euler tiến có thể hoạt động không tốt trong ba tình huống:

• Kích thước bước quá lớn — đa giác dao động hoặc phân kỳ. Cách khắc phục là giảm h; việc so sánh h, h/2, h/4 giúp điều này có thể quan sát được ngay lập tức.
• ODE cứng (Stiff ODEs) — các phương trình có đồng thời các chế độ suy giảm nhanh và suy giảm chậm buộc h phải cực nhỏ để ổn định. Hãy chuyển sang phương pháp ẩn (Euler ngược) hoặc phương pháp BDF.
• Các điểm kỳ dị trong f(x, y) — chia cho không, sqrt của số âm hoặc ln của số không dương sẽ dừng quá trình tích phân. Máy tính báo cáo rõ ràng bước vi phạm.

Các ứng dụng phổ biến

• Vật lý — Định luật thứ hai của Newton dưới dạng hệ cấp một, phân rã phóng xạ \( \dot{N} = -\lambda N \), định luật làm lạnh của Newton.
• Sinh học & dịch tễ học — tăng trưởng logistic \( \dot{y} = r\,y(1 - y/K) \), các mô hình SIR theo ngăn.
• Kinh tế học — lãi kép liên tục, các mô hình tăng trưởng Solow đơn giản.
• Hóa học — động học phản ứng bậc nhất \( \dot{c} = -k c \).
• Giảng dạy — giới thiệu khái niệm tích phân số trước khi chuyển sang RK4 hoặc các trình giải thích ứng.

Câu hỏi thường gặp

Phương pháp Euler là gì?

Phương pháp Euler là quy trình số đơn giản nhất để giải bài toán giá trị ban đầu y' = f(x, y), y(x0) = y0. Tại mỗi bước, nó tiến triển lời giải bằng y_{n+1} = y_n + h · f(x_n, y_n), thực tế là đi theo độ dốc tại điểm hiện tại trong một khoảng cách ngắn h. Nó có độ chính xác cấp một, nghĩa là sai số toàn cục là O(h).

Phương pháp Euler chính xác đến mức nào?

Phương pháp Euler có sai số cắt cụt cục bộ O(h²) và sai số toàn cục O(h). Việc giảm một nửa kích thước bước sẽ làm giảm khoảng một nửa sai số toàn cục. Đây là lý do tại sao việc so sánh hội tụ tại h, h/2 và h/4 trong máy tính này rất hữu ích: bạn có thể thấy sai số thu hẹp xấp xỉ tuyến tính với h.

Khi nào phương pháp Euler thất bại?

Phương pháp Euler có thể trở nên không ổn định đối với các bài toán cứng hoặc khi kích thước bước quá lớn so với độ cong cục bộ của lời giải. Bạn có thể thấy lời giải số dao động, bùng nổ lên vô hạn hoặc trôi ra xa khỏi lời giải thực. Giảm h thường giúp ích; đối với các phương trình cứng, các phương pháp ẩn như Backward Euler được ưu tiên hơn.

Làm thế nào để chọn kích thước bước?

Bắt đầu với h sao cho có khoảng 10 đến 50 bước trong khoảng quan tâm. Nếu đa giác Euler lệch rõ rệt so với trường độ dốc hoặc từ lời giải chính xác của bạn, hãy giảm một nửa h và chạy lại. Sử dụng so sánh h, h/2, h/4 tích hợp để kiểm tra xem ba đường cong có đang hội tụ về phía nhau không.

Sự khác biệt giữa phương pháp Euler và Runge-Kutta (RK4) là gì?

Runge-Kutta bậc bốn đánh giá độ dốc tại bốn điểm mỗi bước và kết hợp chúng với các trọng số (1, 2, 2, 1)/6, mang lại sai số toàn cục O(h⁴) — tốt hơn nhiều bậc so với O(h) của Euler cho cùng một số bước. Euler vẫn có giá trị trong việc giảng dạy khái niệm tích phân số và cho các ứng dụng rất đơn giản hoặc độ chính xác thấp.

Tôi có thể sử dụng công cụ này cho các hệ ODE không?

Máy tính này xử lý một ODE cấp một vô hướng đơn lẻ y' = f(x, y). Đối với các hệ hoặc ODE bậc cao hơn, bạn có thể viết lại phương trình thành một hệ cấp một và sử dụng trình giải hệ chuyên dụng, hoặc chuyển đổi một phương trình bậc hai thành hai phương trình cấp một và giải chúng theo từng thành phần.

Tôi có thể tích phân ngược thời gian không?

Có — hãy nhập kích thước bước h âm. Máy tính sẽ tiến từ x₀ theo hướng âm trong n bước. Điều này hữu ích để tái cấu trúc quá khứ từ một trạng thái hiện tại đã biết.

Đọc thêm

• Phương pháp Euler — Wikipedia
• Phương pháp Runge–Kutta — Wikipedia
• Trường độ dốc — Wikipedia
• Phương trình cứng — Wikipedia