Công thức Euler–Maclaurin

Trong toán học, công thức Euler-Maclaurin là một công thức cho sự khác biệt giữa một tích phân và tổng có liên quan chặt chẽ. Nó có thể được sử dụng để tính gần đúng các tích phân bằng các tổng hữu hạn hoặc ngược lại để đánh giá các tổng hữu hạn và chuỗi vô hạn bằng cách sử dụng các tích phân và máy móc tích phân. Ví dụ, nhiều mở rộng tiệm cận có nguồn gốc từ công thức này và công thức Faulhaber cho tổng số lũy thừa là một hệ quả ngay lập tức.

Công thức được Leonhard Euler và Colin Maclaurin phát hiện độc lập vào khoảng năm 1735. Euler cần nó để tính toán chuỗi hội tụ vô hạn một cách chậm chạp trong khi Maclaurin sử dụng nó để tính tích phân. Sau đó, nó đã được khái quát hóa thành công thức Darboux.

Công thức
Nếu m và n là các số tự nhiên và f(x) là một hàm liên tục có giá trị thực hoặc phức cho các số thực x trong khoảng [m, n], thì tích phân

: I = \int_m^n f(x)\,dx

có thể được xấp xỉ bằng tổng (hoặc ngược lại)

: S = f(m + 1) + \cdots + f(n - 1) + f(n)

(xem phương pháp hình chữ nhật). Công thức Eluler-Maclaurin cung cấp các biểu thức cho sự khác biệt giữa tổng và tích phân theo các đạo hàm cao hơn f^{(k)}(x) đượcđánh giá tại các điểm cuối của khoảng, có nghĩa là khi x=m và x=n.

Một cách cụ thể hơn, cho p một số nguyên dương và một hàm f(x) đó là p lần khả vi liên tục trên khoảng [m, n], chúng ta có

: S - I = \sum_{k=1}^p {\frac{B_k}{k!} (f^{(k - 1)}(n) - f^{(k - 1)}(m))} + R_p,

Ở đâu B_k là k số Bernoulli (với B_1=1/2) và R_p là một thuật ngữ lỗi phụ thuộc vào n, m, p và f và thường nhỏ cho các giá trị phù hợp của p.

Công thức thường được viết với chỉ mục con chỉ lấy các giá trị chẵn, vì các số Bernoulli lẻ bằng 0 ngoại trừ B_1. Trong trường hợp này, chúng ta có

: \sum_{i=m}^n f(i) =
\int^n_m f(x)\,dx + \frac{f(n) + f(m)}{2} +
\sum_{k=1}^{\lfloor p/2\rfloor} \frac{B_{2k}}{(2k)!} (f^{(2k - 1)}(n) - f^{(2k - 1)}(m)) + R_p,

Hay cách khác

: \sum_{i=m+1}^n f(i) =
\int^n_m f(x)\,dx + \frac{f(n) - f(m)}{2} +
\sum_{k=1}^{\lfloor p/2\rfloor} \frac{B_{2k}}{(2k)!} (f^{(2k - 1)}(n) - f^{(2k - 1)}(m)) + R_p.

Số hạng còn lại
Số hạng còn lại phát sinh vì tích phân thường không chính xác bằng tổng. Công thức có thể được bắt nguồn bằng cách áp dụng tích hợp lặp đi lặp lại bởi các phần cho các khoảng tiếp theo [r, r+1] cho r = m, m + 1, \ldots, n - 1. Các cận biên trong các tích hợp này dẫn đến các hệ số chính của công thức và các tích phân còn lại tạo thành "phần còn lại".

Phần còn lại có một biểu thức chính xác với các hàm Bernoulli được định kỳ P_k(x). Đa thức Bernoulli có thể được định nghĩa đệ quy bởi B_0(x) = 1 va cho k \ge 1,

: \begin{align}
B_k'(x) &= kB_{k - 1}(x), \\
\int_0^1 B_k(x)\,dx &= 0.
\end{align}

Các hàm Bernoulli định kỳ được định nghĩa là

: P_k(x) = B_k(x - \lfloor x\rfloor),

Ở đâu \lfloor x \rfloor biểu thị số nguyên lớn nhất nhỏ hơn hoặc bằng x (vậy nên x - \lfloor x \rfloor luôn luôn nằm trong khoảng [0,1)).

Với ký hiệu này, phần còn lại R_p bằng

: R_p = (-1)^{p+1}\int_m^n f^{(p)}(x) \frac{P_p(x)}{p!}\,dx.

Khi nào k>0, có thể chỉ ra rằng

: |B_k(x)| \le \frac{2 \cdot k!}{(2\pi)^k}\zeta(k),

Ở đâu \zeta biểu thị hàm zeta Riemann; một cách tiếp cận để chứng minh sự bất bình đẳng này là thu được chuỗi Fourier cho đa thức B_k(x). Các ràng buộc đạt được cho thậm chí k khi nào x bằng không. Thuật ngữ \zeta(k) có thể được bỏ qua cho số lẻ k nhưng chứng minh trong trường hợp này phức tạp hơn (xem Lehmer). Sử dụng bất đẳng thức này, kích thước của số hạng còn lại có thể được ước tính là

: |R_p| \leq \frac{2 \zeta(p)}{(2\pi)^p}\int_m^n |f^{(p)}(x)|\,dx.

Tham khảo