Máy Giải Phương Trình Vi Phân Bernoulli

Giới thiệu về Máy Giải Phương Trình Vi Phân Bernoulli

Máy giải phương trình vi phân Bernoulli giải quyết một trong những phương trình vi phân bậc nhất phi tuyến nổi tiếng nhất — phương trình Bernoulli y' + P(x)y = Q(x)yⁿ — và biến phép đạo hàm kinh điển trong sách giáo khoa thành một hướng dẫn từng bước mang tính tương tác. Nó tuyến tính hóa phương trình thông qua phép thế v = y¹⁻ⁿ, xây dựng thừa số tích phân μ(x) và ghi đè đường cong dạng đóng kết quả lên nghiệm số RK4 và trường độ dốc để bạn có thể xem mọi chi tiết cùng một lúc.

Phương trình vi phân Bernoulli là gì?

Được giới thiệu bởi Jacob Bernoulli vào năm 1695, phương trình Bernoulli là một phương trình vi phân bậc nhất có dạng

Khi n = 0, phương trình đã là tuyến tính; khi n = 1, nó là phương trình tách biến. Đối với mọi số thực n khác, phương trình là phi tuyến, nhưng phép thế kinh điển v = y¹⁻ⁿ chuyển đổi nó thành phương trình vi phân tuyến tính theo v, có thể giải được bằng mẹo thừa số tích phân tiêu chuẩn.

Phương pháp Bernoulli sáu bước

Bắt đầu từ y' + P(x)y = Q(x)yⁿ:

1. Chia cho yⁿ: \( y^{-n}y' + P(x)\,y^{1-n} = Q(x) \).
2. Thay v = y¹⁻ⁿ: lưu ý rằng \( v' = (1-n)y^{-n}y' \), vì vậy \( y^{-n}y' = v'/(1-n) \).
3. Tuyến tính hóa: \( v' + (1-n)P(x)\,v = (1-n)Q(x) \) — một phương trình vi phân tuyến tính bậc nhất theo v.
4. Thừa số tích phân: \( \mu(x) = \exp\!\left(\int (1-n)P(x)\,dx\right) \), vì vậy \( (\mu v)' = \mu(1-n)Q(x) \).
5. Giải tìm v(x): \( v(x) = \frac{1}{\mu(x)}\left[\mu(x_0)v_0 + \int_{x_0}^{x}\mu(t)(1-n)Q(t)\,dt\right] \).
6. Thế ngược lại: \( y(x) = v(x)^{1/(1-n)} \).

Khi các tích phân liên quan là sơ cấp, bạn sẽ thu được một dạng đóng rõ ràng; khi không phải, máy tính sẽ đánh giá chúng bằng số bằng quy tắc Simpson để vẽ đường cong nghiệm.

Các trường hợp đặc biệt được xử lý tự động

Ví dụ minh họa — n = 2, Dạng Logistic

Xét y' + y/x = x·y² với điều kiện ban đầu y(1) = 1. Ở đây P(x) = 1/x, Q(x) = x, và n = 2, nên 1 − n = −1.

1. Thay v = y⁻¹ = 1/y. Khi đó v' = −y⁻²y' và phương trình trở thành v' − (1/x)v = −x.
2. Thừa số tích phân: μ(x) = exp(∫−1/x dx) = 1/x.
3. (μ·v)' = μ·(−x) = −1. Tích phân: (1/x)·v = −x + C, tức là v = −x² + Cx.
4. Áp dụng ĐKBD: tại x = 1, v = 1/1 = 1, nên 1 = −1 + C ⇒ C = 2. Do đó v(x) = −x² + 2x.
5. Thế ngược lại: y(x) = 1/v(x) = 1/(2x − x²) = 1/(x(2 − x)).

Nghiệm dạng đóng y = 1/(x(2−x)) có các tiệm cận đứng tại x = 0 và x = 2 — chính xác là những gì một trường độ dốc làm cho rõ ràng ngay từ cái nhìn đầu tiên.

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

1. Điền vào trình xây dựng phương trình. Nhập P(x) và Q(x) vào các ô màu xanh, và số mũ n vào ô nhỏ phía trên. Bố cục phản ánh dạng chuẩn y' + P(x)y = Q(x)yⁿ.
2. Thiết lập điều kiện ban đầu (x₀, y₀) và phạm vi vẽ đồ thị [x tối thiểu, x tối đa]. Phạm vi nên chứa x₀.
3. Nhấp vào Giải. Máy tính phát hiện xem bạn có đang ở trong trường hợp đặc biệt (n = 0 hoặc n = 1) hay không và hiển thị đạo hàm tương ứng. Nếu không, nó sẽ thực hiện phép thế Bernoulli đầy đủ sáu bước với các phương trình được hiển thị bằng MathJax.
4. Đọc đồ thị. Đường cong màu cam là nghiệm số RK4. Đường cong đứt nét màu xanh là dạng đóng được đánh giá thông qua thừa số tích phân. Trường mũi tên hiển thị y' ở mọi nơi, vì vậy bạn cũng có thể quan sát các nghiệm khác.
5. Sao chép CSV các điểm mẫu nếu bạn muốn nhập quỹ đạo vào một chương trình khác.

Mẹo, cạm bẫy và các trường hợp biên

• n không phải số nguyên yêu cầu y > 0. Trình giải gắn cờ các tổ hợp như n = 1/2 với y₀ ≤ 0, trong đó yⁿ sẽ là số phức.
• y₀ = 0 thường là điểm kỳ dị. Bất kỳ phương trình Bernoulli nào với Q ≠ 0 và n > 0 đều có nghiệm tầm thường y ≡ 0, đây thường không phải là nhánh bạn muốn.
• Tránh các điểm P(x) bùng nổ gần x₀. Các biểu thức như 1/x yêu cầu x₀ ≠ 0; trình giải sẽ xác thực điều này trước khi chạy.
• Số mũ lớn (|n| > 20) bị từ chối để ngăn hiện tượng tràn số. Trong thực tế, các phương trình Bernoulli với n lớn như vậy hầu như không bao giờ xuất hiện trong các bài toán thực tế.
• Tiệm cận đứng. Nếu RK4 phân kỳ, hãy thử thu hẹp phạm vi x về phía x₀ nơi nghiệm vẫn hữu hạn.

Nơi phương trình Bernoulli xuất hiện

• Động lực học dân số — phương trình logistic y' = ry(1 − y/K) là một phương trình Bernoulli ngụy trang (n = 2 sau khi sắp xếp lại).
• Động học hóa học — các phản ứng tự xúc tác thường tuân theo y' ∝ y − y².
• Mạch điện — một số mạch RL điện trở phi tuyến tạo ra dạng Bernoulli.
• Cơ học chất lưu — phương trình lớp biên sau khi giảm độ tương tự.
• Mô hình dịch bệnh — tỷ lệ cá thể nhạy cảm của mô hình SIR có thể được đưa về dạng Bernoulli.
• Tăng trưởng kinh tế — mô hình Solow–Swan với tỷ lệ tiết kiệm không đổi là Bernoulli với n = α.

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

Phương trình vi phân Bernoulli là gì?

Phương trình Bernoulli là một phương trình vi phân bậc nhất có dạng y' + P(x)y = Q(x)yⁿ, trong đó P và Q là các hàm liên tục và n là một số thực bất kỳ. Đây là một ví dụ kinh điển về phương trình vi phân phi tuyến có thể chuyển đổi thành phương trình tuyến tính thông qua phép thế v = y¹⁻ⁿ.

Phép thế v = y¹⁻ⁿ hoạt động như thế nào?

Nhân phương trình ban đầu với y⁻ⁿ để mọi số hạng y trở thành y¹⁻ⁿ hoặc y⁻ⁿy'. Đặt v = y¹⁻ⁿ cho ta v' = (1−n)y⁻ⁿy'. Việc thay thế sẽ biến đổi phương trình Bernoulli thành v' + (1−n)P(x)v = (1−n)Q(x), đây là phương trình tuyến tính đối với v và có thể giải được bằng thừa số tích phân.

Điều gì xảy ra khi n = 0 hoặc n = 1?

Khi n = 0, phương trình đã là tuyến tính bậc nhất, vì vậy không cần thay thế. Khi n = 1, công thức Bernoulli sẽ chia cho 1 − n = 0, vì vậy chúng tôi xử lý riêng: phương trình thu gọn thành y' = (Q(x) − P(x))·y, đây là phương trình tách biến với nghiệm dạng đóng y = y₀·exp(∫(Q−P) dx).

Các phương trình Bernoulli có luôn giải được dưới dạng đóng không?

Về nguyên tắc là có, nhưng các tích phân thu được liên quan đến thừa số tích phân có thể không có nguyên hàm sơ cấp. Khi điều đó xảy ra, máy tính sẽ đánh giá chúng bằng số bằng quy tắc Simpson và vẽ đường cong nghiệm. Bản thân phương pháp này luôn đưa một phương trình Bernoulli về các phép cầu phương.

Tại sao y âm và n không phải số nguyên lại gây rắc rối?

Nếu n không phải là số nguyên, yⁿ được định nghĩa là exp(n·ln y) và chỉ là số thực khi y > 0. Việc nhập y âm sẽ tạo ra một số phức. Trình giải sẽ gắn cờ tình huống này và yêu cầu y₀ > 0 hoặc một số mũ nguyên để nghiệm duy trì giá trị thực.

Trường độ dốc hiển thị điều gì?

Trường độ dốc là một lưới gồm các đoạn tiếp tuyến nhỏ có góc bằng y' tại điểm (x, y) đó. Bất kỳ đường cong nghiệm nào cũng bị buộc phải tuân theo các tiếp tuyến này, vì vậy trường độ dốc cho phép bạn xem hình dạng định tính của tất cả các nghiệm cùng một lúc, với điều kiện ban đầu xác định đường cong cụ thể.

Đọc thêm

• Phương trình vi phân Bernoulli — Wikipedia
• Thừa số tích phân — Wikipedia
• Hàm Logistic — Wikipedia
• Trường độ dốc — Wikipedia