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

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

Trình giải Phương trình Vi phân Thường cấp hai nhận một phương trình vi phân thường tuyến tính có dạng a·y″ + b·y′ + c·y = g(x) với các hệ số thực hằng số, tự động lập phương trình đặc trưng của nó, phân loại chế độ giảm chấn (quá giảm chấn, giảm chấn tới hạn, thiếu giảm chấn, không giảm chấn hoặc không ổn định), và đưa ra cả nghiệm dạng đóng ký hiệu và nghiệm số học độ chính xác cao. Đầu ra tương tác kết hợp biểu đồ thời gian đường cong kép của y(x) và y′(x) với quỹ đạo mặt phẳng pha của (y, y′) — một góc nhìn tiết lộ chế độ hoạt động ngay lập tức: xoắn ốc vào cho thiếu giảm chấn, nút vào cho quá giảm chấn, vòng lặp kín cho không giảm chấn, xoắn ốc ra cho không ổn định.

Phương trình vi phân thường (ODE) tuyến tính cấp hai với hệ số hằng số là gì?

Một phương trình vi phân thường tuyến tính cấp hai với hệ số hằng số thực là một phương trình có dạng

trong đó a ≠ 0, b, c là các hằng số thực và g(x) là số hạng cưỡng bức. Hai điều kiện đầu y(x₀) = y₀ và y′(x₀) = y′₀ biến nó thành một bài toán giá trị đầu với nghiệm duy nhất trên một lân cận của x₀ — điều này tuân theo định lý Picard-Lindelöf áp dụng cho hệ cấp một tương đương.

Nếu g(x) = 0 phương trình là thuần nhất. Nếu không, nó là không thuần nhất, và nghiệm đầy đủ được phân tích thành

trong đó y_h là nghiệm tổng quát của phương trình thuần nhất liên quan (chứa hai hằng số tự do) và y_p là bất kỳ nghiệm riêng nào của phương trình đầy đủ. Việc áp dụng hai điều kiện đầu sẽ xác định hai hằng số tự do này.

Phương trình đặc trưng

Giả định y = e^(r·x) trong phương trình thuần nhất sẽ cho phương trình đặc trưng (hoặc phương trình trợ giúp)

một phương trình bậc hai có biệt thức Δ = b² − 4ac kiểm soát toàn bộ hành vi định tính:

Ba trường hợp nghiệm & Chế độ giảm chấn

Phương pháp hệ số bất định (Trường hợp không thuần nhất)

Khi g(x) có một trong các dạng đơn giản sau, phương pháp hệ số bất định cung cấp một nghiệm riêng bằng cách giả định một hàm thử có cùng dạng với các hệ số chưa biết và giải tìm chúng:

• Hằng số g(x) = k. Thử: y_p = K. Nếu c = 0 hãy nhân với x; nếu b = 0 nữa, hãy nhân với x một lần nữa.
• Đa thức bậc n. Thử: đa thức tổng quát bậc n. Nhân với x hoặc x² nếu số hạng hằng số hoặc tuyến tính cộng hưởng.
• Hàm mũ g(x) = A·e^(k·x). Thử: y_p = K·e^(k·x). Nếu k trùng với một nghiệm đặc trưng, hãy nhân với x (nghiệm đơn) hoặc x² (nghiệm kép) — đây là cộng hưởng.
• Hình sin g(x) = A·cos(ω·x) + B·sin(ω·x). Thử: y_p = K₁·cos(ω·x) + K₂·sin(ω·x). Nhân với x nếu iω là một nghiệm (cộng hưởng tần số thuần túy).
• Tích và tổng tuân theo tính tuyến tính và quy tắc tích.

Đọc mặt phẳng pha

Hệ cấp một tương đương là u = y, v = y′ với u′ = v và v′ = (g(x) − b·v − c·u)/a. Vẽ v theo u dưới dạng tham số theo x cho ta quỹ đạo mặt phẳng pha. Đối với các hệ tự trị thuần nhất (không có x trong g), các quỹ đạo được xác định duy nhất bởi điểm bắt đầu (y₀, y′₀) và tiết lộ chế độ hoạt động ngay lập tức:

• Thiếu giảm chấn: quỹ đạo xoắn ốc hướng vào trong về phía gốc tọa độ.
• Quá giảm chấn: quỹ đạo tiếp cận gốc tọa độ dọc theo một đường bất biến (vectơ riêng chậm).
• Giảm chấn tới hạn: nút suy biến, quỹ đạo tiếp tuyến với vectơ riêng duy nhất.
• Không giảm chấn: hình elip đóng bao quanh gốc tọa độ — dao động vĩnh cửu.
• Không ổn định: quỹ đạo xoắn ốc hoặc chạy ra ngoài vô tận.

Ví dụ thực tế: Máy dao động hài giảm chấn có cưỡng bức

Xét phương trình y″ + 2·y′ + 5·y = 10 với y(0) = 0, y′(0) = 0 — một hệ thống thiếu giảm chấn, bị cưỡng bức.

1. Phương trình đặc trưng: r² + 2r + 5 = 0 → Δ = 4 − 20 = −16 → r = −1 ± 2i.
2. Nghiệm thuần nhất: y_h = e^(−x)·(C₁·cos 2x + C₂·sin 2x).
3. Nghiệm riêng cho lực cưỡng bức hằng số g = 10: thử y_p = K, nên 5K = 10, cho y_p = 2.
4. Áp dụng ICs: y(0) = 0 → C₁ + 2 = 0 → C₁ = −2. y′(0) = 0 → −C₁ + 2C₂ = 0 → C₂ = −1.
5. Đáp án cuối cùng: y(x) = 2 − e^(−x)·(2·cos 2x + sin 2x) — dao động với biên độ suy giảm và giới hạn y → 2.

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

1. Nhập các hệ số a, b, c ở hàng trên cùng. a phải khác không (nếu không phương trình là cấp một).
2. Nhập số hạng cưỡng bức g(x), hoặc để là 0 cho bài toán thuần nhất. Các nghiệm riêng dạng đóng được lập cho các hằng số, đa thức bậc tối đa 2 và hàm mũ đơn A·e^(k·x) bao gồm cả trường hợp cộng hưởng.
3. Cung cấp các điều kiện đầu (x₀, y₀, y′₀). Cả y và y′ tại x₀ đều phải được chỉ định vì phương trình là cấp hai.
4. Chọn phạm vi x cho các biểu đồ. Trình giải tích phân hướng ra ngoài từ x₀ theo cả hai hướng x bằng RK4.
5. Nhấp vào Giải & Trực quan hóa. Bạn nhận được phương trình đặc trưng với các nghiệm của nó trên mặt phẳng phức, phân loại chế độ giảm chấn, các nghiệm dạng đóng thuần nhất và riêng, biểu đồ thời gian đường cong kép của y và y′, và quỹ đạo mặt phẳng pha.

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

• Hệ thống lò xo-khối lượng-giảm chấn cơ học: m·x″ + c·x′ + k·x = F(t). Quá giảm chấn, giảm chấn tới hạn và thiếu giảm chấn tương ứng với các tỷ số giảm chấn khác nhau ζ = c/(2·√(m·k)).
• Mạch điện RLC: mạch RLC nối tiếp tuân theo L·Q″ + R·Q′ + Q/C = V(t) — cấu trúc giống hệt, các ký hiệu khác nhau.
• Con lắc (góc nhỏ): θ″ + (g/L)·θ = 0 tạo ra dao động điều hòa đơn giản; thêm lực cản không khí tạo ra dao động giảm chấn.
• Phản ứng của tòa nhà đối với động đất: cấu trúc một bậc tự do với gia tốc nền là số hạng cưỡng bức.
• Hệ thống servo điều khiển PID: động lực học sai số vòng kín quy về một ODE cấp hai có tỷ số giảm chấn chi phối độ vọt lố.
• Các mô hình quần thể có quán tính: tăng trưởng kinh tế với độ trễ tích lũy vốn, hoặc các mô hình sinh thái có phản hồi trễ.

Phương pháp số học — Runge-Kutta cổ điển (RK4) trên hệ 2D

Công cụ giảm a·y″ + b·y′ + c·y = g(x) về hệ cấp một

với u(x₀) = y₀, v(x₀) = y′₀. Runge-Kutta bốn giai đoạn sau đó được áp dụng cho trạng thái vectơ (u, v). RK4 có sai số cắt địa phương O(h⁵) và sai số toàn cầu O(h⁴); 400 bước phụ mặc định trong mỗi hướng mang lại độ chính xác khoảng sáu chữ số cho các bài toán không cứng.

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

Phương trình vi phân thường (ODE) tuyến tính cấp hai với hệ số hằng số là gì?

Một ODE tuyến tính cấp hai với hệ số hằng số có dạng a·y″ + b·y′ + c·y = g(x), trong đó a, b, c là các hằng số thực và g(x) là số hạng cưỡng bức (không thuần nhất). Với hai điều kiện đầu y(x₀) = y₀ và y′(x₀) = y′₀ nghiệm là duy nhất. Trường hợp thuần nhất g(x) = 0 luôn có nghiệm dạng đóng thông qua phương trình đặc trưng a·r² + b·r + c = 0; trường hợp không thuần nhất được giải dưới dạng y(x) = y_h(x) + y_p(x).

Phương trình đặc trưng là gì?

Đối với a·y″ + b·y′ + c·y = 0, thay thế giả định y = e^(r·x) dẫn đến a·r² + b·r + c = 0 — phương trình đặc trưng hoặc phương trình trợ giúp. Các nghiệm của nó quyết định dạng của nghiệm thuần nhất: hai nghiệm thực phân biệt cho y_h = C₁·e^(r₁·x) + C₂·e^(r₂·x); một nghiệm kép r cho y_h = (C₁ + C₂·x)·e^(r·x); các nghiệm phức liên hợp α ± β·i cho y_h = e^(α·x)·(C₁·cos(β·x) + C₂·sin(β·x)).

Thiếu giảm chấn, giảm chấn tới hạn và quá giảm chấn có nghĩa là gì?

Các thuật ngữ này đến từ mô hình lò xo-khối lượng-giảm chấn m·x″ + c·x′ + k·x = 0. Quá giảm chấn (biệt thức > 0, hai nghiệm thực) có nghĩa là hệ thống trở về trạng thái cân bằng chậm mà không có dao động. Giảm chấn tới hạn (biệt thức = 0, nghiệm kép) là sự trở về nhanh nhất mà không có vọt lố. Thiếu giảm chấn (biệt thức < 0, nghiệm phức) tạo ra dao động suy giảm. Không giảm chấn (b = 0, c/a > 0) tạo ra dao động hình sin thuần túy mãi mãi.

Phương pháp hệ số bất định là gì?

Đối với lực cưỡng bức g(x) đơn giản — hằng số, đa thức, hàm mũ, sin, cos và tích của chúng — nghiệm riêng y_p được giả định có cùng dạng với g với các hệ số chưa biết, được xác định bằng cách thay vào ODE và so khớp các số hạng. Hàm thử phải được nhân với x (hoặc x² cho nghiệm kép) khi g(x) cộng hưởng với một nghiệm đặc trưng.

Mặt phẳng pha là gì?

Đối với một phương trình cấp hai được quy về hệ 2D (y, y′) mặt phẳng pha vẽ đồ thị y′ theo y khi x tiến tới. Các đường cong nghiệm trong mặt phẳng pha tiết lộ chế độ hoạt động ngay lập tức: các xoắn ốc suy giảm cho thiếu giảm chấn, các nút hướng vào trong cho quá giảm chấn, các hình elip đóng cho chuyển động hài không giảm chấn, và các xoắn ốc hướng ra ngoài cho dao động không ổn định. Nó là bản đối chiếu hình học của sơ đồ nghiệm phương trình đặc trưng.

Công cụ này sử dụng phương pháp số học nào?

Phương pháp Runge-Kutta bậc bốn (RK4) cổ điển được áp dụng cho hệ cấp một tương đương u = y, v = y′, với u′ = v và v′ = (g(x) − b·v − c·u)/a. RK4 có sai số cắt địa phương O(h⁵) và 400 bước phụ mặc định cho mỗi hướng mang lại độ chính xác khoảng sáu chữ số cho các phương trình không cứng trên cửa sổ đã chọn.

Đọc thêm

• Phương trình vi phân tuyến tính — Wikipedia
• Phương trình đặc trưng — Wikipedia
• Phương pháp hệ số bất định — Wikipedia
• Máy dao động hài — Wikipedia
• Mặt phẳng pha — Wikipedia
• Phương pháp Runge-Kutta — Wikipedia