✨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\; =\; \backslash int\_m^n\; f(x)\backslash ,dx$

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

: $S\; =\; f(m\; +\; 1)\; +\; \backslash 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\; =\; \backslash sum\_\{k=1\}^p\; \{\backslash 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ó

: $\backslash sum\_\{i=m\}^n\; f(i)\; =\; \backslash 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

: $\backslash sum\_\{i=m+1\}^n\; f(i)\; =\; \backslash 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,\; \backslash 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\; \backslash ge\; 1$,

: $\backslash 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\; -\; \backslash lfloor\; x\backslash rfloor),$

Ở đâu $\backslash lfloor\; x\; \backslash rfloor$ biểu thị số nguyên lớn nhất nhỏ hơn hoặc bằng $x$ (vậy nên $x\; -\; \backslash lfloor\; x\; \backslash 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\}\backslash int\_m^n\; f^\{(p)\}(x)\; \backslash frac\{P\_p(x)\}\{p!\}\backslash ,dx.$

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

: $|B\_k(x)|\; \backslash le\; \backslash frac\{2\; \backslash cdot\; k!\}\{(2\backslash pi)^k\}\backslash zeta(k),$

Ở đâu $\backslash 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ữ $\backslash 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|\; \backslash leq\; \backslash frac\{2\; \backslash zeta(p)\}\{(2\backslash pi)^p\}\backslash int\_m^n\; |f^\{(p)\}(x)|\backslash ,dx.$