Phương pháp đa bước tuyến tính
Các phương pháp đa bước tuyến tính tính toán mỗi giá trị nghiệm mới từ một tổ hợp tuyến tính của một số giá trị nghiệm và đạo hàm trước đó, tái sử dụng các kết quả đã có để đạt được bậc cao với chi phí thấp cho mỗi bước.
Definition
Phương pháp đa bước tuyến tính là một phương pháp cho các phương trình vi phân thường, xác định giá trị nghiệm tiếp theo thông qua một quan hệ tuyến tính cố định giữa một số giá trị nghiệm trước đó và các đánh giá vế phải.
Scope
Chủ đề này bao gồm các họ Adams-Bashforth (hiện) và Adams-Moulton (ẩn), các công thức vi phân ngược cho các bài toán cứng, triển khai dự đoán-hiệu chỉnh, các đa thức đặc trưng và điều kiện nghiệm xác định tính ổn định zero, và các rào cản bậc của Dahlquist giới hạn những gì các phương pháp như vậy có thể đạt được.
Core questions
- Các phương pháp đa bước tái sử dụng các giá trị trước đó như thế nào để đạt được bậc cao với một lần đánh giá hàm mới cho mỗi bước?
- Ổn định zero là gì, và điều kiện nghiệm trên đa thức đặc trưng thể hiện nó như thế nào?
- Các cặp dự đoán-hiệu chỉnh kết hợp các công thức hiện và ẩn trong thực tế như thế nào?
- Các rào cản bậc của Dahlquist nói gì về giới hạn của độ chính xác và ổn định của phương pháp đa bước?
Key theories
- Ổn định zero và điều kiện nghiệm
- Một phương pháp đa bước là ổn định zero, và do đó hội tụ khi nhất quán, chính xác khi các nghiệm của đa thức đặc trưng đầu tiên của nó nằm trong đĩa đơn vị đóng với chỉ các nghiệm đơn trên biên; điều kiện nghiệm này là tương tự đa bước của tính ổn định.
- Các rào cản Dahlquist
- Rào cản đầu tiên của Dahlquist giới hạn bậc của một phương pháp k-bước ổn định zero, và rào cản thứ hai của ông cho thấy không có phương pháp đa bước tuyến tính ổn định A nào có thể có bậc lớn hơn hai, đó là lý do tại sao các bộ giải cứng bậc cao dựa vào sự thỏa hiệp BDF về ổn định tương đối hơn là ổn định tuyệt đối.
Mechanisms
Các phương pháp Adams tích phân một đa thức nội suy thông qua các giá trị đạo hàm trước đó: Adams-Bashforth chỉ sử dụng các giá trị đã biết (hiện), Adams-Moulton bao gồm giá trị mới chưa biết (ẩn) để đạt được độ chính xác và ổn định cao hơn. Trong thực tế, hai phương pháp này được ghép nối thành một cặp dự đoán-hiệu chỉnh: công thức hiện dự đoán, công thức ẩn hiệu chỉnh, thường trong một hoặc hai lần lặp. Các công thức vi phân ngược thay vào đó lấy hiệu các giá trị nghiệm trước đó để xấp xỉ đạo hàm tại điểm mới, tạo ra các phương pháp ổn định cứng là cốt lõi của các mã ODE cứng. Bởi vì các phương pháp đa bước cần một số giá trị khởi đầu, chúng được khởi động bằng một phương pháp một bước.
Clinical relevance
Các phương pháp đa bước tuyến tính, đặc biệt là các công thức vi phân ngược, là nền tảng cho các bộ giải ODE cứng được sử dụng trong động học hóa học, mô phỏng mạch điện tử và các hệ thống vi phân-đại số lớn, nơi việc đánh giá vế phải tốn kém và việc tái sử dụng các đánh giá trước đó thông qua các công thức đa bước mang lại hiệu quả đáng kể.
History
Adams và Bashforth đã giới thiệu các công thức đa bước vào thế kỷ XIX, với Moulton bổ sung các biến thể ẩn; phân tích của Dahlquist vào những năm 1950-1960 đã thiết lập lý thuyết ổn định và các rào cản bậc chi phối lĩnh vực này, và công trình của C. William Gear vào những năm 1970 đã biến các mã công thức vi phân ngược thành tiêu chuẩn cho các bài toán cứng.
Key figures
- John Couch Adams
- Francis Bashforth
- Forest Ray Moulton
- Germund Dahlquist
- C. William Gear
Related topics
Seminal works
- hairer1993
- iserles2008
Frequently asked questions
- Các phương pháp đa bước khác với các phương pháp Runge-Kutta như thế nào?
- Các phương pháp Runge-Kutta thực hiện một số đánh giá đạo hàm mới trong mỗi bước nhưng loại bỏ chúng sau đó, trong khi các phương pháp đa bước tái sử dụng các giá trị đạo hàm từ các bước trước. Do đó, các phương pháp đa bước rẻ hơn cho mỗi bước nhưng cần các giá trị khởi đầu bổ sung và xử lý đặc biệt khi thay đổi kích thước bước.
- Điều kiện nghiệm là gì?
- Đó là yêu cầu các nghiệm của đa thức đặc trưng đầu tiên của phương pháp phải nằm bên trong hoặc trên đường tròn đơn vị, với các nghiệm biên là đơn. Nó đảm bảo rằng các lỗi nhỏ không bị khuếch đại khi các bước tích lũy, đảm bảo phương pháp ổn định zero và do đó hội tụ.