I. Khái niệm cơ bản về phương pháp tuyến tính đa bước
Phương pháp tuyến tính đa bước là một lớp phương pháp số quan trọng trong giải các phương trình vi phân thường. Đây là phương pháp sử dụng thông tin từ nhiều bước tính toán trước đó để xấp xỉ nghiệm tại bước hiện tại. Khác với phương pháp một bước như Euler, phương pháp đa bước tận dụng dữ liệu lịch sử để nâng cao độ chính xác và hiệu quả tính toán. Các phương pháp này được xây dựng dựa trên nội suy đa thức Newton và công thức Simpson, cho phép tính toán với cấp chính xác cao. Việc hiểu rõ cấu trúc và nguyên lý của phương pháp tuyến tính đa bước là nền tảng để áp dụng hiệu quả trong thực tế.
1.1. Định nghĩa và công thức tổng quát
Phương pháp tuyến tính đa bước được định nghĩa bằng công thức: yn+k = a1yn+k-1 + a2yn+k-2 + ... + akyn + h(b0f(xn+k,yn+k) + b1f(xn+k-1,yn+k-1) + ... + bkf(xn,yn)). Các hệ số a và b được xác định dựa trên điều kiện nội suy. Công thức này cho phép tính yn+k từ k giá trị trước đó, giúp cải thiện hiệu suất tính toán.
1.2. Phân loại phương pháp tuyến tính đa bước
Các phương pháp tuyến tính đa bước chia thành hai loại chính: phương pháp hiển (explicit) khi b0 = 0, và phương pháp ẩn (implicit) khi b0 ≠ 0. Phương pháp hiển tính toán dễ dàng nhưng ổn định kém, trong khi phương pháp ẩn có tính ổn định tốt hơn nhưng cần giải phương trình phi tuyến.
II. Phương pháp Adams trong giải số phương trình vi phân
Phương pháp Adams là một trong những phương pháp tuyến tính đa bước phổ biến nhất, được chia thành Adams-Bashforth (hiển) và Adams-Moulton (ẩn). Các phương pháp này xây dựng dựa trên nội suy đa thức các giá trị hàm f tại các điểm trước. Phương pháp Adams có ưu điểm là dễ cài đặt, có cấp chính xác cao, và tính toán hiệu quả. Phương pháp Adams-Bashforth thường được sử dụng để dự đoán giá trị, trong khi Adams-Moulton dùng để hiệu chỉnh. Sự kết hợp hai phương pháp này tạo thành công thức dự đoán-hiệu chỉnh rất hiệu quả trong thực hành.
2.1. Phương pháp Adams Bashforth
Phương pháp Adams-Bashforth là phương pháp hiển sử dụng các giá trị f tại k bước trước. Công thức k-bước: yn+1 = yn + h∑βjfn-j với các hệ số βj xác định bằng nội suy lùi. Ưu điểm là tính toán đơn giản, nhưng ổn định yếu với các bước h lớn.
2.2. Phương pháp Adams Moulton
Phương pháp Adams-Moulton là phương pháp ẩn bao gồm giá trị f tại bước hiện tại. Công thức: yn+1 = yn + h(β₀fn+1 + ∑βjfn-j+1) với β₀ ≠ 0. Phương pháp này có tính ổn định tốt hơn và cấp chính xác cao, thích hợp cho các bài toán khó.
III. Phương pháp Nyström Milne Simpson và BDF
Ngoài các phương pháp Adams, còn có phương pháp Nyström dùng cho phương trình vi phân bậc hai, và phương pháp Milne-Simpson kết hợp hai phương pháp khác nhau. Phương pháp BDF (Backward Differentiation Formula) là lớp phương pháp tuyến tính đa bước ẩn, xây dựng từ nội suy ngược hàm y thay vì f. Phương pháp BDF đặc biệt hiệu quả cho các phương trình vi phân stiff, nhờ tính ổn định A-ổn định của nó. Các phương pháp BDF bậc cao (2, 3 bước) được sử dụng rộng rãi trong các ứng dụng kỹ thuật phức tạp.
3.1. Phương pháp Nyström và ứng dụng
Phương pháp Nyström được thiết kế cho phương trình bậc hai y'' = f(x,y,y') không chứa y'. Công thức: yn+1 = yn-1 + 2h·f(xn,yn) + O(h³). Phương pháp này tiết kiệm tính toán vì không cần đưa về hệ phương trình bậc một.
3.2. Phương pháp BDF và tính ổn định
Phương pháp BDF k-bước định nghĩa bởi: ∑αjyn+j = h·βkfn+k với αj từ nội suy ngược. Công thức BDF 2 bước: (3yn+2 - 4yn+1 + yn)/2h = fn+2 có cấp chính xác 2 và A-ổn định, lý tưởng cho phương trình stiff.
IV. Tính chất toán học của phương pháp tuyến tính đa bước
Để đánh giá hiệu quả của phương pháp tuyến tính đa bước, cần xét các tính chất toán học quan trọng bao gồm cấp chính xác, tính phù hợp, zero-ổn định và sự hội tụ. Cấp chính xác xác định bậc của sai số cắt xén địa phương, được xác định thông qua khai triển Taylor. Tính phù hợp đòi hỏi phương pháp cho kết quả chính xác cho các đa thức bậc thấp. Zero-ổn định đảm bảo sai số không phát triển khi h→0. Cuối cùng, định lý hội tụ Dahlquist khẳng định rằng phương pháp hội tụ khi và chỉ khi tính phù hợp và zero-ổn định được thỏa mãn. Những tính chất này là cơ sở lý thuyết cho sự ứng dụng thành công của phương pháp tuyến tính đa bước.
4.1. Cấp chính xác và sai số cắt xén
Cấp chính xác p của phương pháp tuyến tính đa bước được định nghĩa sao cho sai số cắt xén địa phương τn = O(hp). Sai số toàn cục en = y(xn) - yn đáp ứng en = O(hp). Để đạt cấp chính xác cao, cần xác định chính xác các hệ số a và b từ điều kiện hội tụ.
4.2. Zero ổn định và định lý Dahlquist
Zero-ổn định kiểm tra xem phương pháp có ổn định khi f = 0 (phương trình y' = 0). Điều kiện: các nghiệm của đa thức đặc trưng ρ(r) = rk - a₁rk-1 - ... - ak phải có |r| ≤ 1 với nghiệm đơn |r| = 1. Định lý Dahlquist phát biểu: nếu zero-ổn định và tính phù hợp được thỏa mãn thì phương pháp hội tụ.
