Giải Hệ Phương Trình Vi Phân Thường

Giới thiệu về Giải Hệ Phương Trình Vi Phân Thường

Trình giải Hệ Phương trình Vi phân Thường là một hộp công cụ vi phân tất cả trong một cho các hệ tuyến tính và phi tuyến kết hợp. Nhập ma trận hệ số 2×2 hoặc 3×3 và công cụ sẽ thực hiện phân tích đầy đủ giá trị riêng / vectơ riêng, viết nghiệm tổng quát và nghiệm riêng dạng đóng bằng LaTeX, phân loại điểm cân bằng tại gốc tọa độ là điểm yên ngựa, nút, xoắn ốc hoặc tâm, và vẽ biểu đồ pha tương tác với quỹ đạo hoạt hình. Đối với các hệ phi tuyến phẳng, bạn có thể nhập các vế phải tùy ý \(f(x,y)\) và \(g(x,y)\) và công cụ sẽ tạo ra biểu đồ pha RK4 với độ chính xác cao.

Hệ Phương trình Vi phân Thường (ODEs) là gì?

Một hệ phương trình vi phân thường liên kết nhiều hàm chưa biết của một biến duy nhất — thường là thời gian \(t\) — thông qua các đạo hàm của chúng. Ở dạng gọn nhất,

Khi \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\) với ma trận hằng số \(A\), hệ thống là tuyến tính và tự trị — và đây là nơi lý thuyết đạt đến vẻ đẹp rực rỡ nhất: toàn bộ hành vi dài hạn được xác định bởi các giá trị riêng của \(A\).

Công thức Giải Hệ Tuyến tính bằng Giá trị riêng

Đối với \(\mathbf{x}' = A\mathbf{x}\), phương pháp tiêu chuẩn là:

1. Tính đa thức đặc trưng \(\det(\lambda I - A) = 0\).
2. Giải tìm các giá trị riêng \(\lambda_1, \lambda_2, \dots\).
3. Với mỗi giá trị riêng, tìm một vectơ riêng \(v\) bằng cách giải \((A - \lambda I) v = 0\).
4. Lập nghiệm tổng quát dưới dạng tổ hợp tuyến tính: \(\mathbf{x}(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} + \cdots\).
5. Xác định các hằng số \(c_i\) bằng cách thay điều kiện ban đầu \(\mathbf{x}(0)\) vào nghiệm tổng quát.

Ba trường hợp cho Hệ 2×2

Mặt phẳng Vết-Định thức

Đối với ma trận 2×2 có vết \(T = a_{11} + a_{22}\) và định thức \(D = a_{11} a_{22} - a_{12} a_{21}\), toàn bộ phân loại khớp trong một sơ đồ duy nhất:

Đây là lý do tại sao bảng kết quả hiển thị nổi bật \(T\), \(D\), và \(\Delta = T^2 - 4D\) — ba con số là đủ để đặt tên cho điểm cân bằng.

Hệ Phi tuyến và Biểu đồ Pha

Hầu hết các ODE thực tế là phi tuyến và không có nghiệm dạng đóng. Công cụ xử lý chúng bằng cách tích phân các phương trình theo phương pháp số Runge–Kutta bậc 4 (RK4), có sai số cắt cục bộ \(O(h^5)\) và là phương pháp mặc định cho các trường vectơ trơn.

Lớp phủ biểu đồ pha bao gồm:

• Một trường vectơ được lấy mẫu trên lưới 13×13, hiển thị hướng dòng chảy tại mọi điểm.
• Quỹ đạo từ điều kiện ban đầu của bạn, được vẽ màu đỏ với một chấm cam hoạt hình hiển thị hướng thời gian.
• Một vài đường dòng mầm từ một vòng tròn các điểm bắt đầu, đưa ra một bức tranh tổng thể về động lực học.
• Đối với các hệ tuyến tính 2×2, các trục vectơ riêng (màu xanh lơ đứt đoạn) — đây là các hướng bất biến mà các nghiệm trượt theo hàm mũ.

Cách sử dụng Trình giải này

1. Chọn một chế độ — Tuyến tính 2×2, Tuyến tính 3×3, hoặc Phi tuyến 2D — thông qua các tab ở đầu biểu mẫu.
2. Điền vào các hệ số hoặc phương trình. Nhấp vào bất kỳ Ví dụ nhanh nào để điền trước một hệ tiêu chuẩn (nút ổn định, tâm, yên ngựa, con lắc, Van der Pol, v.v.).
3. Nhập điều kiện ban đầu \((x_0, y_0)\) và một khoảng thời gian \(T\). Các giá trị \(T\) điển hình là 6–20 cho các bộ dao động và 3–6 cho các hệ ổn định suy giảm nhanh.
4. Nhấp Giải. Trang kết quả đầy đủ sẽ xuất hiện với phân loại, giá trị riêng, vectơ riêng, nghiệm dạng đóng (chế độ tuyến tính), biểu đồ pha hoạt hình và biểu đồ chuỗi thời gian.
5. Xem lại quỹ đạo bằng cách sử dụng nút dưới biểu đồ pha nếu bạn muốn xem lại quá trình di chuyển trên đường cong IC.

Ví dụ Thực tế — Bộ dao động điều hòa có giảm chấn

Bộ dao động có giảm chấn \(\ddot{x} + 2\zeta \omega \dot{x} + \omega^2 x = 0\) có thể được viết lại dưới dạng hệ 2D bằng cách đặt \(y = \dot{x}\):

Với \(\omega = 1\) và \(\zeta = 0.2\) (giảm chấn dưới mức), ma trận là \(A = \begin{pmatrix} 0 & 1 \\ -1 & -0.4 \end{pmatrix}\). Vết \(T = -0.4\), định thức \(D = 1\), biệt số \(\Delta = 0.16 - 4 = -3.84 < 0\), vì vậy chúng ta có một xoắn ốc ổn định với các giá trị riêng \(-0.2 \pm 0.9798\,i\). Quỹ đạo xoắn ốc vào gốc tọa độ, và chuỗi thời gian hiển thị các hàm sin suy giảm theo hàm mũ.

Ứng dụng

• Hệ thống cơ khí — hệ thống lò xo-khối lượng kết hợp, con lắc, con quay hồi chuyển.
• Mạch điện — mạng RLC, bộ lọc op-amp, điều khiển không gian trạng thái.
• Động lực học quần thể — mô hình thú săn mồi-con mồi Lotka–Volterra, các loài cạnh tranh, dịch tễ học (SIR, SIS).
• Động học hóa học — mạng lưới phản ứng, bộ dao động Belousov–Zhabotinsky.
• Khoa học thần kinh — mô hình nơ-ron FitzHugh–Nagumo, rút gọn Hodgkin–Huxley.
• Lý thuyết điều khiển — mô hình nhà máy tuyến tính hóa, thiết kế bộ quan sát, biên độ ổn định.

Mẹo & Lưu ý

• Nếu quỹ đạo của bạn tăng vọt quá nhanh, hãy giảm khoảng thời gian T — một hệ thống không ổn định có thể tràn khỏi khung hình chỉ trong vài đơn vị thời gian.
• Đối với giá trị riêng lặp lại, trình giải sẽ tự động tìm vectơ riêng tổng quát \(w\) bằng cách giải \((A - \lambda I)w = v\), do đó bạn nhận được số hạng \(tv\) mà không cần tính toán thủ công.
• Đối với hệ phi tuyến, các mũi tên trường vectơ cũng tiết lộ các điểm cân bằng khác gốc tọa độ dưới dạng các chấm màu xanh lơ — hãy quan sát biểu đồ để tìm các vùng có độ lớn bằng không.
• Đối với hệ 3×3 không có biểu đồ pha (3D khó hiển thị trên trang 2D), nhưng nghiệm dạng đóng và kết luận về độ ổn định vẫn được áp dụng.
• Điều kiện ban đầu và khoảng thời gian tách biệt với việc phân loại: thay đổi chúng chỉ làm di chuyển quỹ đạo màu đỏ, không làm thay đổi kết luận về giá trị riêng.

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

Hệ phương trình vi phân thường là gì?

Hệ phương trình vi phân thường (ODEs) là một tập hợp các phương trình liên kết với nhau liên quan đến đạo hàm của nhiều hàm chưa biết của một biến độc lập duy nhất, thường là thời gian. Dạng cổ điển là \( \mathbf{x}'(t) = F(t, \mathbf{x}(t)) \), trong đó \( \mathbf{x} \) là một vectơ trạng thái và \(F\) là trường vectơ. Các hệ tuyến tính có thể được viết gọn dưới dạng \( \mathbf{x}' = A\mathbf{x} + \mathbf{b} \), và hành vi của chúng được xác định gần như hoàn toàn bởi các giá trị riêng của ma trận hệ số \(A\).

Làm thế nào để các giá trị riêng phân loại điểm cân bằng của hệ tuyến tính 2×2?

Đối với hệ 2×2 \( \mathbf{x}' = A\mathbf{x} \), gốc tọa độ được phân loại theo vết \(T\) và định thức \(D\) của \(A\): \(D < 0\) cho điểm yên ngựa (không ổn định); \(D > 0\) với \(T^2 > 4D\) cho một nút (ổn định nếu \(T < 0\), không ổn định nếu \(T > 0\)); \(D > 0\) với \(T^2 < 4D\) cho một đường xoắn ốc (ổn định nếu \(T < 0\), không ổn định nếu \(T > 0\), một tâm thuần túy nếu \(T = 0\)). Đường biên \(T^2 = 4D\) tạo ra một nút suy biến.

Nghiệm dạng đóng trông như thế nào khi các giá trị riêng là số phức?

Nếu \(A\) có các giá trị riêng phức liên hợp \( \alpha \pm i\beta \) với vectơ riêng phức \( v = p + iq \), nghiệm tổng quát thực là \( \mathbf{x}(t) = e^{\alpha t} \left[ c_1 (p \cos\beta t - q \sin\beta t) + c_2 (p \sin\beta t + q \cos\beta t) \right] \). Hàm mũ \(e^{\alpha t}\) kiểm soát biên độ (tăng, giảm hoặc hằng số), trong khi sin và cos xử lý chuyển động quay.

Điều gì xảy ra khi ma trận có giá trị riêng lặp lại?

Nếu ma trận có giá trị riêng lặp lại \(\lambda\) nhưng chỉ có một vectơ riêng độc lập tuyến tính \(v\), bạn cũng cần một vectơ riêng tổng quát \(w\) thỏa mãn \( (A - \lambda I) w = v \). Nghiệm tổng quát khi đó có dạng \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \). Nếu không gian riêng có hai chiều, ma trận là một bội số vô hướng của ma trận đơn vị trên không gian con bất biến đó và nghiệm rút gọn thành \( \mathbf{x}(t) = (c_1 v_1 + c_2 v_2) e^{\lambda t} \).

Công cụ này có thể giải các hệ phi tuyến dưới dạng biểu tượng không?

Chế độ phi tuyến giải hệ thống bằng phương pháp số bằng cách sử dụng bộ tích phân Runge–Kutta bậc 4 (RK4) và vẽ biểu đồ pha. Hầu hết các hệ phi tuyến không có nghiệm dạng đóng, vì vậy đây là cách tiếp cận tiêu chuẩn. Bạn vẫn có thể đọc ra hành vi cục bộ gần các điểm cân bằng bằng cách tuyến tính hóa, điều mà chế độ tuyến tính 2×2 xử lý được — hãy tính ma trận Jacobian tại điểm cố định và nhập nó vào dưới dạng \(A\).

Biểu đồ pha là gì?

Biểu đồ pha là một hình ảnh hình học về các nghiệm của một hệ 2D trong mặt phẳng \(x\)–\(y\). Mỗi nghiệm vạch ra một đường cong gọi là quỹ đạo, và tập hợp các quỹ đạo cùng với các mũi tên trường vectơ tiết lộ hành vi định tính: liệu các nghiệm xoắn ốc vào trong, tách ra tại điểm yên ngựa, dao động hay ổn định tại các điểm cân bằng. Biểu đồ pha giúp cấu trúc toàn cục của một hệ thống có thể nhìn thấy được trong nháy mắt.

Đọc thêm

• Hệ phương trình vi phân — Wikipedia
• Biểu đồ pha — Wikipedia
• Giá trị riêng và vectơ riêng — Wikipedia
• Các phương pháp Runge–Kutta — Wikipedia
• Bộ dao động Van der Pol — Wikipedia