Công thức xấp xỉ Stirling cơ bản nhất có dạng biểu thức như thế nào?

Công thức Stirling cơ bản xác định giai thừa n! xấp xỉ tiệm cận với căn bậc hai của 2 pi n nhân với (n chia e) tất cả lũy thừa n, với tỷ số giữa hai vế tiến về 1 khi n tiến ra vô cùng.

Bất đẳng thức kép Robbins cải tiến công thức Stirling như thế nào?

Herbert Robbins thiết lập bất đẳng thức kép kẹp chặt giá trị giai thừa n! với mọi số nguyên n >= 1 bằng cách bổ sung số hạng e^(1/(12n+1)) ở cận dưới và e^(1/(12n)) ở cận trên.

Tại sao xấp xỉ Stirling lại giữ vai trò cốt lõi trong cơ học thống kê?

Trong cơ học thống kê hạt Boltzmann, entropy vi mô được tính qua logarit tự nhiên của số trạng thái vi mô chứa tổ hợp giai thừa phân bố hạt; với hệ vĩ mô cỡ số Avogadro, dạng rút gọn ln n! ~ n ln n - n là công cụ giải tích then chốt.

Chuỗi tiệm cận Stirling có phải là chuỗi hội tụ vô hạn thông thường không?

Chuỗi Stirling là chuỗi tiệm cận theo nghĩa Poincaré; với một giá trị n cố định, chuỗi sẽ phân kỳ nếu cộng vô hạn số hạng, do đó trong thực tế người ta phải cắt cụt chuỗi ở số hạng tối ưu.

Xấp xỉ Stirling là gì? Công thức, chứng minh và ứng dụng

Xấp xỉ Stirling (Stirling's approximation), hay công thức Stirling, là một hệ thức tiệm cận toán học dùng để xấp xỉ giá trị của giai thừa các số nguyên dương lớn hoặc mở rộng sang hàm Gamma trên trường số phức. Trong giải tích toán học, tổ hợp, lý thuyết xác suất và cơ học thống kê, phép tính giai thừa trực tiếp đối với các đại lượng lớn đòi hỏi chi phí tính toán bùng nổ, do đó xấp xỉ Stirling đóng vai trò là một công cụ phân tích giải tích không thể thiếu. Bài viết này trình bày chi tiết lịch sử khám phá, các phương pháp chứng minh giải tích, chuỗi tiệm cận hiệu chỉnh sai số, các bất đẳng thức chặn kẹp chính xác và các miền ứng dụng tiêu biểu của công thức Stirling.

Bối cảnh lịch sử và nguồn gốc khám phá

Lịch sử ra đời của công thức Stirling gắn liền với sự phát triển của lý thuyết xác suất cổ điển vào đầu thế kỷ mười tám, đặc biệt là bài toán ước tính hệ số nhị thức trong phân phối nhị thức với số lần thử rất lớn. Nhà toán học người Pháp Abraham de Moivre trong ấn bản năm 1730 của công trình kinh điển Miscellanea Analytica đã lần đầu tiên phát hiện ra rằng khi số nguyên $n$ tăng lên vô hạn, giai thừa $n!$ tăng trưởng tiệm cận tỷ lệ thuận với biểu thức $n^{n + 1/2} e^{-n}$ . Tuy nhiên, de Moivre chưa thể xác định được giá trị giải tích chính xác của hằng số tỷ lệ mà chỉ ước tính được giá trị xấp xỉ số học của nó.

Cũng vào năm 1730, nhà toán học người Scotland James Stirling đã công bố chuyên khảo Methodus Differentialis. Bằng cách kết hợp phép tính vi phân hữu hạn với công thức tích phân tích số vô hạn nổi tiếng của John Wallis, Stirling đã chứng minh xuất sắc rằng hằng số tỷ lệ còn thiếu trong công thức của de Moivre chính là căn bậc hai của tích 2 nhân số pi (xấp xỉ 2,5066). Như được hệ thống hóa và phân tích lịch sử tường minh trong công trình của Tweddle (2003), sự hợp tác và hoàn thiện này đã chính thức khai sinh ra công thức mang tên xấp xỉ Stirling:

n! \sim \sqrt{2\pi n} \left(\frac{n}{e}\right)^n

Trong toán học giải tích, ký hiệu ngã $\sim$ biểu thị quan hệ tương đương tiệm cận, nghĩa là giới hạn của tỷ số giữa hai vế tiến về 1 khi biến số tiến dần ra vô cùng:

\lim_{n \to \infty} \frac{n!}{\sqrt{2\pi n} (n/e)^n} = 1

Phương pháp chứng minh giải tích thông qua hàm Gamma

Để mở rộng phạm vi của xấp xỉ Stirling từ tập số tự nhiên sang tập số thực dương và số phức, các nhà toán học sử dụng biểu diễn giải tích của hàm Gamma của Leonhard Euler:

\Gamma(n+1) = \int_0^\infty t^n e^{-t} \, dt = n!

Phương pháp chứng minh hiện đại và tổng quát nhất cho công thức Stirling dựa trên phương pháp Laplace (hay phương pháp điểm yên ngựa). Biểu thức dưới dấu tích phân có thể được viết lại dưới dạng hàm mũ:

t^n e^{-t} = \exp(n \ln t - t)

Hàm số ở số mũ đạt giá trị cực đại duy nhất tại điểm $t = n$ . Thực hiện phép đổi biến số $t = n(1 + x)$ , tích phân Euler chuyển thành:

n! = n^{n+1} e^{-n} \int_{-1}^\infty \exp\left( n [\ln(1+x) - x] \right) \, dx

Khai triển chuỗi Taylor cho hàm số logarit quanh lân cận điểm gốc $x = 0$ cho ta:

\ln(1+x) - x = -\frac{x^2}{2} + \frac{x^3}{3} - \dots

Khi $n$ rất lớn, toàn bộ giá trị của tích phân tập trung cục bộ xung quanh điểm cực đại $x = 0$ . Bỏ qua các số hạng bậc cao hơn và lấy xấp xỉ tích phân Gaussian chuẩn trên toàn bộ trục thực:

\int_{-\infty}^\infty \exp\left(-\frac{n x^2}{2}\right) \, dx = \sqrt{\frac{2\pi}{n}}

Nhân kết quả này với hệ số bên ngoài tích phân trực tiếp dẫn tới công thức xấp xỉ bậc nhất Stirling.

Chuỗi tiệm cận Stirling và số Bernoulli

Để đạt độ chính xác cao hơn khi tính toán cho các giá trị $n$ hữu hạn, công thức Stirling được mở rộng thành một chuỗi tiệm cận đầy đủ thông qua công thức tổng Euler-Maclaurin. Lấy logarit tự nhiên của hàm giai thừa, chuỗi khai triển Stirling có dạng:

\ln n! = n \ln n - n + \frac{1}{2} \ln(2\pi n) + \sum_{k=1}^m \frac{B_{2k}}{2k(2k-1) n^{2k-1}} + \mathcal{O}\left(\frac{1}{n^{2m+1}}\right)

Trong công thức trên, các hệ số $B_{2k}$ đại diện cho các số Bernoulli chẵn. Thay thế các giá trị cụ thể của các số Bernoulli, nghiên cứu của Nemes (2010) trình bày các số hạng điều chỉnh đầu tiên của chuỗi tiệm cận Stirling:

\ln n! = n \ln n - n + \frac{1}{2} \ln(2\pi n) + \frac{1}{12n} - \frac{1}{360n^3} + \frac{1}{1260n^5} - \dots

Dạng lũy thừa tương ứng của chuỗi này đối với hàm giai thừa được biểu diễn qua khai triển:

n! \sim \sqrt{2\pi n} \left(\frac{n}{e}\right)^n \left( 1 + \frac{1}{12n} + \mathcal{O}\left(\frac{1}{n^2}\right) \right)

Bất đẳng thức chặn kẹp Robbins

Một điểm hạn chế của chuỗi tiệm cận là chuỗi này phân kỳ đối với mọi giá trị $n$ cố định khi số lượng số hạng tiến ra vô cùng. Do đó, trong giải tích số và xác suất ứng dụng, các nhà toán học cần những bất đẳng thức chặn trên và chặn dưới có hiệu lực phổ quát cho mọi số nguyên dương. Năm 1955, nhà toán học Herbert Robbins đã công bố một bất đẳng thức kép kinh điển trên tạp chí The American Mathematical Monthly:

\sqrt{2\pi n} \left(\frac{n}{e}\right)^n \exp\left(\frac{1}{12n+1}\right) < n! < \sqrt{2\pi n} \left(\frac{n}{e}\right)^n \exp\left(\frac{1}{12n}\right)

Bất đẳng thức của Robbins (1955) có giá trị lý thuyết và ứng dụng to lớn vì nó đúng tuyệt đối với mọi số nguyên $n \ge 1$ . Công thức này cho thấy sai số tương đối của xấp xỉ Stirling giảm nhanh theo quy luật bậc nghịch đảo của $n$ , cung cấp một công cụ kiểm soát sai số giải tích tin cậy mà không đòi hỏi phải tính toán chuỗi phức tạp.

Các công thức cải tiến hiện đại và đóng góp của Ramanujan

Trong cuốn sổ tay toán học của mình, thiên tài toán học Ấn Độ Srinivasa Ramanujan đã đề xuất một công thức xấp xỉ giai thừa có tốc độ hội tụ cực nhanh. Trong công trình khảo cứu chi tiết của Hirschhorn và Villarino (2013), công thức xấp xỉ của Ramanujan được chứng minh với dạng biểu thức giải tích tường minh:

n! \approx \sqrt{\pi} \left(\frac{n}{e}\right)^n \sqrt[6]{8n^3 + 4n^2 + n + \frac{1}{30}}

Nghiên cứu của Hirschhorn và Villarino (2013) đã khẳng định rằng công thức Ramanujan cho độ chính xác phi thường ngay cả với trường hợp biên nhỏ nhất khi $n = 1$ , với sai số tương đối nhỏ hơn rất nhiều so với công thức Stirling sơ cấp.

Tiếp nối hướng nghiên cứu này, các công trình của Mortici (2010) đã phát triển các lớp xấp xỉ tiệm cận tối ưu mới cho hàm Gamma, cho phép kiểm soát sai số với tốc độ tiệm cận đạt bậc cao hơn bằng cách tham số hóa các phân thức hữu tỷ bổ trợ trong số mũ của công thức Stirling.

Công thức xấp xỉ	Tác giả và thời gian	Dạng biểu thức xấp xỉ	Bậc hội tụ sai số
Stirling sơ cấp	Stirling & de Moivre (1730)	Căn bậc hai của 2 pi n nhân (n/e)^n	Bậc một O(1/n)
Khai triển Euler-Maclaurin	Nemes (2010) hệ thống	Thêm số hạng 1/(12n) - 1/(360n^3) + 1/(1260n^5)	Tùy ý theo số hạng cắt cụt
Chặn kẹp Robbins	Robbins (1955)	Kẹp giữa exp(1/(12n+1)) và exp(1/(12n))	Chặn chặt phổ quát với mọi n >= 1
Công thức Ramanujan	Hirschhorn & Villarino (2013)	Căn bậc 6 của đa thức (8n^3 + 4n^2 + n + 1/30)	Bậc siêu cao O(1/n^4)
Xấp xỉ Mortici	Mortici (2010)	Tối ưu hóa đa thức tham số trong hàm mũ	Tăng tốc hội tụ tiệm cận

Miền ứng dụng trong khoa học và kỹ thuật

Công thức xấp xỉ Stirling là một trong những công cụ toán học được vận dụng nhiều nhất trong khoa học tự nhiên và tính toán hiện đại:

Cơ học thống kê và nhiệt động lực học: Trong cơ học thống kê Boltzmann, entropy vi mô của một hệ nhiệt động chứa $N$ hạt được xác định qua công thức liên hệ với trọng lượng thống kê: $S = k_B \ln \Omega$ . Đại lượng $\Omega$ thường bao gồm các tổ hợp giai thừa phân bố hạt. Khi $N$ đạt tới số Avogadro (ở quy mô số Avogadro với bậc độ lớn hạt khổng lồ), xấp xỉ dạng logarit $\ln N! \approx N \ln N - N$ là công cụ duy nhất cho phép suy dẫn các phương trình trạng thái khí lý tưởng và phân bố vi chính tắc.
Lý thuyết xác suất và định lý giới hạn: Định lý giới hạn trung tâm de Moivre-Laplace là nền tảng xấp xỉ phân phối nhị thức bằng phân phối chuẩn Gaussian khi số phép thử độc lập rất lớn. Việc chứng minh sự hội tụ của hệ số nhị thức sang hàm mật độ chuẩn dựa trực tiếp vào xấp xỉ Stirling.
Khoa học máy tính và lý thuyết thuật toán: Trong phân tích độ phức tạp thuật toán, cận dưới lý thuyết của thuật toán sắp xếp so sánh là $\log_2(n!)$ . Ứng dụng xấp xỉ Stirling cho phép thiết lập ngay lập tức định lý cận dưới tiệm cận $\Omega(n \log n)$ , chứng minh tính tối ưu của các thuật toán sắp xếp như Merge Sort hay Heapsort.

Hạn chế toán học và lưu ý khi sử dụng

Mặc dù là công cụ xấp xỉ mạnh mẽ, việc áp dụng công thức Stirling cần lưu ý các giới hạn giải tích quan trọng:

Bản chất phân kỳ của chuỗi tiệm cận: Chuỗi Stirling là một chuỗi tiệm cận theo nghĩa Poincaré, không phải là một chuỗi hội tụ theo nghĩa Cauchy thông thường. Với một giá trị $n$ cố định, nếu cộng thêm vô hạn số hạng trong chuỗi, tổng chuỗi sẽ phân kỳ ra vô cùng. Do đó, trong tính toán số học thực tế, người ta phải cắt cụt chuỗi ở số hạng tối ưu có độ lớn nhỏ nhất.
Sai số tương đối ở giá trị n nhỏ: Công thức Stirling sơ cấp cho sai số tương đối khoảng vài phần trăm khi $n$ nhỏ. Nếu không áp dụng các hệ số hiệu chỉnh của Robbins (1955) hay công thức Ramanujan, sai số này có thể dẫn đến chệch hướng nghiêm trọng trong các mô hình tổ hợp quy mô nhỏ.
Điểm kỳ dị trên mặt phẳng phức: Khi mở rộng sang hàm Gamma $\Gamma(z)$ , công thức Stirling chỉ hợp lệ trong miền góc phức không chứa trục thực âm, bởi hàm Gamma sở hữu các cực điểm đơn tại các số nguyên không dương $z = 0, -1, -2, \dots$ .