Trình Giải Phương Trình Vi Phân Thường Cấp Một

Giới thiệu về Trình Giải Phương Trình Vi Phân Thường Cấp Một

Trình giải phương trình vi phân thường cấp một tiếp nhận một phương trình vi phân thường dưới dạng dy/dx = f(x, y), tự động phân loại cấu trúc của nó (tách biến, tuyến tính, tự trị, chính xác hoặc tổng quát), và tạo ra cả nghiệm giải tích dạng đóng nếu có thể và nghiệm số độ chính xác cao ở mọi nơi. Trực quan hóa trường hướng trực tiếp với đường cong nghiệm đè lên giúp ý nghĩa hình học của phương trình trở nên rõ ràng ngay lập tức — các nghiệm chính là các đường cong tiếp xúc với mọi mũi tên.

Phương trình vi phân thường cấp một là gì?

Một phương trình vi phân thường cấp một (ODE) chỉ bao gồm một hàm số chưa biết y(x) và đạo hàm cấp một của nó y'(x). Dạng hiện tiêu chuẩn là:

Kết hợp với một điều kiện ban đầu y(x₀) = y₀, điều này định nghĩa một bài toán giá trị ban đầu (IVP). Định lý Picard-Lindelöf đảm bảo một nghiệm duy nhất trong một lân cận nào đó của x₀ miễn là f là liên tục Lipschitz đối với y gần (x₀, y₀). Về mặt hình học, IVP yêu cầu tìm đường cong duy nhất đi qua (x₀, y₀) có độ dốc tại mọi điểm khớp với f tại điểm đó — chính xác là đường cong tiếp xúc với trường hướng.

Sáu lớp phương trình mà bộ giải nhận diện

Phương pháp dạng đóng: Tuyến tính với hệ số hằng số

Khi vế phải rút gọn thành dy/dx = a·y + b với các hằng số a và b, thừa số tích phân μ(x) = e^(-a·x) cho một nghiệm chính xác. Nghiệm tổng quát là:

Áp dụng điều kiện ban đầu y(x₀) = y₀ xác định hằng số C và mang lại nghiệm riêng duy nhất. Lớp phương trình này bao gồm rất nhiều bài toán trong sách giáo khoa:

• Tăng trưởng mũ — dy/dx = k·y, nghiệm riêng y(t) = y₀·e^(k·t).
• Suy giảm mũ — dy/dx = -k·y, chu kỳ bán rã ln 2 / k.
• Định luật làm mát của Newton — dy/dx = -k·(y - T_amb), nhiệt độ cơ thể giảm theo hàm mũ tiến về nhiệt độ môi trường.
• Sạc mạch RC — dV/dt = (1/RC)·(V_in - V), điện áp tụ điện tiếp cận nguồn.
• Đào thải thuốc — dược động học cấp một với tốc độ đào thải k.

Cách đọc một trường hướng

Tại mỗi điểm lưới (x, y), công cụ vẽ một đoạn thẳng ngắn có độ dốc bằng f(x, y). Ba quan sát hữu ích:

• Điểm cân bằng là những điểm mà f(x, y) = 0 — trường hướng nằm ngang. Đối với các phương trình tự trị, đây là các điểm cố định y* thỏa mãn f(y*) = 0; các quỹ đạo lân cận hoặc tiếp cận (ổn định) hoặc chạy xa khỏi (không ổn định) y*.
• Đường đẳng hướng (Isoclines) là những đường cong mà f(x, y) bằng một hằng số c, do đó tất cả các mũi tên dọc theo đường cong đều có cùng độ dốc c.
• Các đường cong nghiệm không bao giờ cắt nhau (khi f là Lipschitz) — rõ ràng về mặt trực quan vì hai đường cong cắt nhau sẽ cần các độ dốc khác nhau tại điểm giao nhau.

Phương pháp số: Runge-Kutta cổ điển (RK4)

Cho trước (x_n, y_n), giá trị tiếp theo được tính bằng cách lấy trung bình của bốn ước tính độ dốc:

RK4 có sai số cắt cụ bộ O(h⁵) và sai số toàn cục O(h⁴), mang lại độ chính xác khoảng sáu chữ số ở số bước mặc định cho các phương trình không cứng. Bộ giải tích phân từ điểm ban đầu theo cả hai hướng x và dừng lại nếu độ lớn của y vượt quá 10¹⁵ — thường thấy ở các nghiệm bùng nổ trong thời gian hữu hạn, như dy/dx = y².

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

1. Nhập vế phải vào trường dy/dx = ... Sử dụng x và y làm biến, * để nhân, ^ hoặc ** cho lũy thừa, và các hàm tiêu chuẩn như sin, cos, exp, log, sqrt. Các hằng số pi và e được công nhận.
2. Chỉ định điều kiện ban đầu (x₀, y₀) — đường cong nghiệm duy nhất sẽ đi qua điểm này.
3. Chọn phạm vi x để vẽ trường hướng và đường cong nghiệm. Phạm vi y được tự động điều chỉnh từ nghiệm tích phân.
4. Nhấp Giải & Trực quan hóa. Bộ phân loại chạy trước; nếu phương trình của bạn khớp với mẫu dạng đóng (tuyến tính với hệ số hằng), bạn sẽ nhận được câu trả lời ký hiệu. Trường hướng và đường cong nghiệm luôn được hiển thị.
5. Bật/tắt trường hướng để tập trung vào đường cong nghiệm, hoặc xem lại hoạt ảnh vẽ đường cong để xem quá trình tích phân tiến triển như thế nào từ điểm ban đầu.

Ví dụ thực tế: Định luật làm mát của Newton

Một tách cà phê ở 80 °C nguội đi trong phòng 20 °C. Tốc độ truyền nhiệt tỷ lệ thuận với sự chênh lệch nhiệt độ:

Đây là phương trình tuyến tính với hệ số hằng số (a = -0.1, b = 2). Dạng đóng là:

Sau 30 phút: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. Chế độ xem trường hướng giúp hành vi giới hạn trở nên rõ ràng — mọi đường cong nghiệm, bất kể nhiệt độ bắt đầu, đều tiệm cận với đường ngang T = 20.

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

• Động lực học quần thể — các mô hình tăng trưởng mũ, logistic, hiệu ứng Allee.
• Dược động học — hấp thụ và đào thải thuốc, tính toán chu kỳ bán rã.
• Truyền nhiệt — định luật làm mát của Newton, mô hình điện dung gộp.
• Mạch RC và RL — các quá trình quá độ điện tuyến tính cấp một.
• Phân rã phóng xạ — chuỗi phân rã đơn đồng vị.
• Bình trộn — nồng độ của chất tan dưới tác động của dòng vào/ra.
• Vật rơi có lực cản — phân tích vận tốc cuối dv/dt = g - kv.

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

Phương trình vi phân thường cấp một là gì?

Phương trình vi phân thường (ODE) cấp một là một phương trình có dạng dy/dx = f(x, y) liên quan đến hàm số chưa biết y(x) và đạo hàm cấp một của nó. Giải ODE nghĩa là tìm hàm số y(x) có đạo hàm khớp với vế phải. Với một điều kiện ban đầu y(x₀) = y₀, nghiệm là duy nhất dưới các giả định quy tắc nhẹ (định lý Picard-Lindelöf).

Trường hướng là gì?

Trường hướng (hoặc trường hướng vectơ) vẽ một đoạn thẳng nhỏ tại mỗi điểm lưới (x, y) có độ dốc bằng f(x, y). Các đường cong nghiệm của ODE chính là các đường tiếp xúc với các đoạn thẳng này tại mọi điểm. Trường hướng cung cấp trực giác hình học tức thì về hành vi tổng thể của các nghiệm mà không cần giải phương trình bằng ký hiệu.

Công cụ này giải được những loại ODE cấp một nào?

Công cụ tự động phân loại phương trình vào một trong các loại: có thể tích phân (chỉ phụ thuộc vào x, giải bằng tích phân trực tiếp), tuyến tính với hệ số hằng số y' = a·y + b (cung cấp nghiệm dạng đóng đầy đủ), tự trị (chỉ phụ thuộc vào y), tách biến (phân tích thành g(x)·h(y)), tuyến tính với hệ số biến đổi (P(x)·y + Q(x)), hoặc tổng quát. Đối với mọi loại, một nghiệm số Runge-Kutta độ chính xác cao và trực quan hóa trường hướng sẽ được tạo ra.

Phương pháp số nào được sử dụng?

Phương pháp Runge-Kutta bậc bốn cổ điển (RK4) được áp dụng với 300 bước phụ trong mỗi hướng từ điểm ban đầu. RK4 có sai số cắt cụ bộ O(h⁵) và là phương pháp tiêu chuẩn cho các ODE không cứng ở quy mô này. Bộ giải phát hiện sự phân kỳ (tràn số hoặc NaN) và dừng tích phân một cách sạch sẽ để biểu đồ vẫn hợp lệ.

Phương pháp thừa số tích phân cho ODE tuyến tính là gì?

Đối với một ODE tuyến tính cấp một y' + P(x)·y = Q(x), nhân cả hai vế với thừa số tích phân μ(x) = e^∫P(x) dx. Vế trái trở thành đạo hàm chính xác d/dx[μ·y], vì vậy y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C). Khi P và Q là hằng số, nó rút gọn thành dạng đóng y = -b/a + C·e^(a·x), mà công cụ sẽ tự động trả về.

Công cụ này có thể xử lý các phương trình cứng hoặc hệ phương trình ODE không?

Bộ giải này dành cho các ODE vô hướng cấp một không cứng. Các vấn đề rất cứng (nơi nghiệm có nhiều thang thời gian khác nhau nhiều bậc độ lớn) có thể cần một phương pháp ẩn như Euler ngược hoặc Rosenbrock; các hệ liên kết yêu cầu một bộ giải giá trị vectơ. Trong những trường hợp đó, hãy sử dụng một gói chuyên dụng như solve_ivp của SciPy hoặc các bộ giải ODE cứng chuyên dụng.

Đọc thêm

• Phương trình vi phân thường — Wikipedia
• Trường hướng (Slope field) — Wikipedia
• Phương pháp Runge-Kutta — Wikipedia
• Thừa số tích phân (Integrating factor) — Wikipedia
• Định lý Picard-Lindelöf — Wikipedia